单自由度系统受迫振动与共振现象仿真分析

1. 项目概述：单自由度系统受迫振动与共振现象

在结构动力学领域，单自由度系统的受迫振动是最基础却最具代表性的研究课题之一。作为一名长期从事结构仿真分析的工程师，我发现许多实际工程振动问题（从桥梁风振到机械设备异常）都可以简化为单自由度系统模型来理解其本质。本次仿真主题聚焦于系统在简谐激励下的响应特性，特别是当激励频率接近系统固有频率时出现的共振现象——这个看似简单的物理现象，曾导致过塔科马海峡大桥坍塌等著名工程事故。

受迫振动分析的核心价值在于：通过掌握幅频特性曲线的变化规律，我们能够预测结构在不同频率激励下的动态响应，进而避免共振带来的灾难性后果。在航空航天、机械制造和土木工程等领域，这种分析是结构安全性评估的必修课。现代工程仿真软件（如ANSYS、ADINA等）虽然能自动完成计算，但只有理解单自由度系统的理论基础，才能正确解读仿真结果并指导设计。

2. 核心理论解析

2.1 系统动力学方程建立

典型的单自由度系统由质量块m、弹簧k和阻尼器c组成（如图1所示）。当系统受到简谐激励力P(t)=P₀sin(ωt)作用时，其运动微分方程为：

code复制mẍ + cẋ + kx = P₀sin(ωt)

这个二阶常微分方程的解由齐次解（自由振动）和特解（稳态受迫振动）组成。在实际工程中，我们更关注稳态响应，其解析解为：

code复制x(t) = Asin(ωt - φ)

其中振幅A和相位角φ的表达式揭示了系统响应的本质规律：

code复制A = (P₀/k) / √[(1-(ω/ωₙ)²)² + (2ζω/ωₙ)²]
φ = arctan[2ζ(ω/ωₙ)/(1-(ω/ωₙ)²)]

关键参数说明：

• ωₙ=√(k/m)：系统固有频率
• ζ=c/(2√(km))：阻尼比
• β=A/(P₀/k)：动态放大系数

2.2 幅频特性曲线特征

通过改变激励频率ω并记录振幅A，可以得到图2所示的典型幅频曲线。曲线呈现三个特征区域：

1. 准静态区（ω/ωₙ << 1）：
   响应振幅接近静态变形P₀/k，相位滞后可忽略
   （此时结构刚度起主导作用）

2. 共振区（ω≈ωₙ）：
   振幅出现峰值，且阻尼比ζ显著影响峰值高度
   （实测中常用半功率带宽法测定ζ）