自动控制——过阻尼、欠阻尼、临界阻尼及无阻尼

引言

在自动控制系统和振动分析中，系统的阻尼特性对于系统的动态响应至关重要。阻尼决定了系统在受到扰动或输入信号后，如何恢复到稳定状态。本文将详细介绍过阻尼（overdamped）、欠阻尼（underdamped）、临界阻尼（critically damped）和无阻尼（undamped）的定义、公式，并通过Python代码演示不同阻尼条件下系统的响应。

系统的阻尼特性

考虑一个典型的二阶线性系统，其传递函数可以表示为：

$$H(s)=\frac{{\omega }_n^2}{s^2+2\zeta {\omega }_{n}s+{\omega }_n^2}$$

其中：

-${\omega }_n$是系统的自然频率。
-$\zeta$是阻尼比，定义为$\zeta =\frac{c}{2\sqrt{mk}}$，其中$c$是阻尼系数，$m$是质量，$k$是弹簧常数。

根据阻尼比$\zeta$的不同，系统的动态响应分为以下几种类型：

1. 无阻尼（Undamped）

当$\zeta =0$时，系统没有阻尼力的影响，其响应为纯粹的正弦波动，永不衰减。无阻尼系统在自然频率${\omega }_n$下振荡。

$$x(t)=A\cos ({\omega }_{n}t)+B\sin ({\omega }_{n}t)$$

2. 欠阻尼（Underdamped）

当$0<\zeta <1$时，系统处于欠阻尼状态。系统响应在频率${\omega }_d={\omega }_n\sqrt{1-{\zeta }^2}$下振荡，并随着时间指数衰减。

$$x(t)=e^{-\zeta {\omega }_{n}t}\left(A\cos ({\omega }_{d}t)+B\sin ({\omega }_{d}t)\right)$$

欠阻尼系统的特征是系统会振荡几次，然后逐渐趋于平衡。

3. 临界阻尼（Critically Damped）

当$\zeta =1$时，系统处于临界阻尼状态。此时系统响应速度最快，不会出现振荡，但也不会像过阻尼那样慢慢回到平衡状态。

$$x(t)=(A+Bt)e^{-{\omega }_{n}t}$$

临界阻尼是最理想的阻尼状态，它能确保系统迅速回到平衡而不出现振荡。

4. 过阻尼（Overdamped）

当$\zeta >1$时，系统处于过阻尼状态。系统响应不再振荡，但相比于临界阻尼系统，它恢复到平衡所需的时间更长。

$$x(t)=A{e}^{({\lambda }_{1}t)}+B{e}^{({\lambda }_{2}t)}$$

其中，${\lambda }_1$和${\lambda }_2$是系统的特征根，且${\lambda }_1,{\lambda }_2<0$。

Python代码实现

下面的Python代码演示了一个简单的二阶系统在不同阻尼条件下的时间响应。

import numpy as np
import matplotlib.pyplot as plt
from scipy.signal import lti, step

# 定义系统的自然频率和时间范围
omega_n = 5.0  # 自然频率
t = np.linspace(0, 5, 500)

# 定义不同阻尼比的系统
systems = {
    'Undamped (ζ=0)': lti([omega_n**2], [1, 0, omega_n**2]),
    'Underdamped (ζ=0.5)': lti([omega_n**2], [1, 2*0.5*omega_n, omega_n**2]),
    'Critically Damped (ζ=1)': lti([omega_n**2], [1, 2*1.0*omega_n, omega_n**2]),
    'Overdamped (ζ=2)': lti([omega_n**2], [1, 2*2.0*omega_n, omega_n**2])
}

# 绘制每个系统的阶跃响应
plt.figure(figsize=(10, 6))
for label, system in systems.items():
    t, y = step(system, T=t)
    plt.plot(t, y, label=label)

# 图表设置
plt.title('不同阻尼条件下的系统阶跃响应')
plt.xlabel('时间 [s]')
plt.ylabel('响应')
plt.grid(True)
plt.legend()
plt.show()

代码解释

• lti: scipy.signal.lti 函数用于定义线性时不变（LTI）系统。
• step: step 函数用于计算系统的阶跃响应。
• systems: 定义了四种不同阻尼比$\zeta$的系统：无阻尼、欠阻尼、临界阻尼和过阻尼。
  

结论

通过上述分析和代码演示，我们可以看到阻尼比$\zeta$对系统动态响应的影响：

• 无阻尼系统（$\zeta =0$）会永远振荡，不会衰减。
• 欠阻尼系统（$0<\zeta <1$）会振荡几次，随着时间的推移逐渐平稳。
• 临界阻尼系统（$\zeta =1$）是理想的，它在不振荡的情况下快速回到平衡。
• 过阻尼系统（$\zeta >1$）不会振荡，但比临界阻尼系统更慢地回到平衡状态。

理解和设计适当的阻尼对于工程应用至关重要，能够有效避免系统的过冲、延迟或持续振荡。在实际应用中，通常会设计系统使其接近临界阻尼状态，以确保快速而稳定的响应。