1. 单自由度系统受迫振动与共振现象概述
在结构动力学领域，单自由度系统的受迫振动是最基础却至关重要的研究课题。这类系统由一个质量块、一个弹簧和一个阻尼器组成，当受到外部周期性激励时，系统会产生特定的振动响应。共振现象则是当激励频率接近系统固有频率时，振幅急剧增大的特殊状态。
作为一名长期从事结构仿真分析的工程师，我发现许多实际工程问题（如桥梁风振、机械振动等）都可以简化为单自由度系统来研究其核心动力学特性。通过仿真分析，我们能够预测结构在动态载荷下的行为，这对工程安全评估和优化设计具有重要意义。
2. 理论基础与数学模型
2.1 系统动力学方程
单自由度系统的受迫振动可以用二阶常微分方程描述：
code复制mẍ + cẋ + kx = F0sin(ωt)
其中：
- m：系统质量(kg)
- c：阻尼系数(N·s/m)
- k：弹簧刚度(N/m)
- F0：激励力幅值(N)
- ω：激励频率(rad/s)
这个方程揭示了质量惯性力、阻尼力、弹性恢复力与外部激励力之间的动态平衡关系。
2.2 稳态响应解
系统达到稳态时的位移响应为：
code复制x(t) = Asin(ωt - φ)
振幅A和相位角φ的表达式为：
code复制A = (F0/k) / √[(1-(ω/ωn)^2)^2 + (2ζω/ωn)^2]
φ = arctan[2ζ(ω/ωn)/(1-(ω/ωn)^2)]
其中：
- ωn = √(k/m)：系统固有频率
- ζ = c/(2√(mk))：阻尼比
3. 共振现象解析
3.1 共振条件
当激励频率ω接近系统固有频率ωn时，系统进入共振状态。此时：
- 振幅达到最大值：
code复制Amax = F0/(2kζ√(1-ζ^2)) - 相位差φ接近90°
- 能量输入效率最高
3.2 频率响应特性
通过绘制幅频特性曲线（如图1所示），可以观察到：
- 当ω/ωn << 1时：响应由刚度控制
- 当ω/ωn ≈ 1时：出现共振峰
- 当ω/ωn >> 1时：响应由质量控制
4. 仿真实现步骤
4.1 建模方法
-
质量-弹簧-阻尼系统建模：
python复制import numpy as np from scipy.integrate import odeint def system(y, t, m, c, k, F0, omega): x, v = y dxdt = v dvdt = (F0*np.sin(omega*t) - c*v - k*x)/m return [dxdt, dvdt] -
参数设置：
python复制m = 1.0 # 质量(kg) k = 100.0 # 刚度(N/m) c = 0.5 # 阻尼系数(N·s/m) F0 = 10.0 # 激励幅值(N) omega_n = np.sqrt(k/m) # 固有频率
4.2 求解与可视化
-
时域响应求解：
python复制t = np.linspace(0, 10, 1000) y0 = [0, 0] # 初始条件 omega = 0.8*omega_n # 激励频率 sol = odeint(system, y0, t, args=(m,c,k,F0,omega)) -
结果可视化：
python复制import matplotlib.pyplot as plt plt.plot(t, sol[:,0]) plt.xlabel('Time (s)') plt.ylabel('Displacement (m)') plt.title('Forced Vibration Response') plt.grid() plt.show()
5. 工程应用与注意事项
5.1 典型应用场景
-
机械振动分析：
- 旋转机械的临界转速确定
- 减震器设计优化
-
建筑抗震设计：
- 结构基频计算
- 地震响应谱分析
-
航空航天：
- 飞行器颤振分析
- 有效载荷隔振设计
5.2 仿真实践要点
-
阻尼比选择：
- 钢结构：ζ≈0.01-0.05
- 混凝土结构：ζ≈0.05-0.10
- 带阻尼器结构：ζ可达0.20以上
-
频率扫描技巧：
- 先用粗扫确定共振区
- 在共振区附近采用细扫
- 建议频率步长Δω ≤ 0.1ωn
-
结果验证方法：
- 检查能量平衡
- 验证稳态响应幅值
- 对比解析解与数值解
6. 常见问题与解决方案
6.1 数值不稳定问题
现象：高频激励时解算发散
解决方法：
- 减小时间步长
- 使用隐式积分方法
- 检查阻尼系数设置
6.2 共振峰识别困难
现象：幅频曲线峰值不明显
可能原因：
- 阻尼过大（ζ > 0.3）
- 频率分辨率不足
- 激励幅值太小
6.3 工程应用建议
-
安全裕度设置：
- 工作频率应避开0.7ωn ~ 1.3ωn范围
- 对于关键结构，建议保持ω/ωn < 0.5或>1.5
-
减振措施：
- 增加阻尼（阻尼器、粘弹性材料）
- 调频设计（改变质量或刚度）
- 动力吸振器应用
7. 进阶分析技巧
7.1 参数化研究
通过系统化改变参数，研究其对响应的影响：
python复制zeta_values = [0.01, 0.05, 0.1, 0.2]
for zeta in zeta_values:
c = 2*zeta*np.sqrt(m*k)
# 运行仿真并绘制曲线
7.2 频响函数计算
计算并绘制频响函数(FRF)：
python复制omega_range = np.linspace(0, 2*omega_n, 200)
Amp = []
for omega in omega_range:
sol = odeint(system, y0, t, args=(m,c,k,F0,omega))
Amp.append(np.max(sol[-100:,0])) # 取稳态振幅
7.3 实验与仿真对比
建议将仿真结果与实验数据进行对比验证：
- 测量实际结构的频率响应
- 修正仿真模型参数
- 进行相关性分析
通过这种系统的仿真分析方法，工程师可以深入理解单自由度系统的动力学特性，为更复杂的多自由度系统分析奠定基础。在实际应用中，建议结合理论计算、数值仿真和实验验证三种手段，以获得可靠的分析结果。