4. 结构动力学基础：单自由度体系振动方程、固有频率与振型、阻尼比的概念、共振原理

各位工程师朋友，咱们今天聊点硬核的。结构动力学，说白了就是研究结构在动荷载下怎么“动”的学问。你想想看，塔筒在风里摇摆，它到底是怎么振起来的？频率是多少？会不会越晃越大最后垮掉？这些问题的答案，都藏在这一章的基础里。

我个人习惯，讲动力学一定先从最简单的模型入手。别一上来就搞几十个自由度的有限元模型，那会把人绕晕。咱们先看一个弹簧上挂个铁球，这就是经典的单自由度体系。

4.1 单自由度体系振动方程

一个质量块，一根弹簧，一个阻尼器，外加一个外荷载。这就是单自由度体系的全部家当。它的运动方程长这样：

m * ẍ(t) + c * ẋ(t) + k * x(t) = F(t)

解释一下：

• m：质量，单位kg。塔筒的等效质量。
• c：阻尼系数，单位N·s/m。说白了就是阻碍运动的力。
• k：刚度，单位N/m。结构抵抗变形的能力。
• ẍ(t)：加速度，位移对时间的二阶导。
• ẋ(t)：速度，位移对时间的一阶导。
• x(t)：位移，我们最关心的量。
• F(t)：外荷载，比如风荷载。

这个方程，是结构动力学的“牛顿第二定律”。我在项目中遇到过不少年轻工程师，拿着有限元软件算了一堆结果，但问他这个方程怎么写，答不上来。嗯，这里要注意：不懂这个方程，你永远只是软件的操作工，不是工程师。

4.2 固有频率与振型

如果外荷载F(t)=0，结构会怎么动？这就是自由振动。这时候方程变成：

m * ẍ(t) + c * ẋ(t) + k * x(t) = 0

先忽略阻尼（c=0），方程简化为：

m * ẍ(t) + k * x(t) = 0

这个方程的解，是一个简谐振动：

x(t) = A * sin(ωₙ * t + φ)

其中：

• ωₙ：无阻尼固有圆频率，单位rad/s。ωₙ = √(k/m)。
• fₙ：固有频率，单位Hz。fₙ = ωₙ / (2π)。
• Tₙ：固有周期，单位s。Tₙ = 1 / fₙ。

为什么会这样？因为结构一旦被扰动，它就会按照自己的“脾气”振动。这个“脾气”就是固有频率。刚度越大、质量越小，频率越高；反之越低。

对于多自由度体系（比如塔筒），每个自由度都有自己的固有频率和对应的振型。振型，说白了就是结构在某个频率下振动的“形状”。一阶振型像摆锤，二阶振型像S形，三阶振型像W形。

4.3 阻尼比的概念

现实中没有永动机。结构振动总会慢慢停下来，为什么？因为有阻尼。阻尼消耗能量，让振动衰减。

阻尼比ζ（zeta）是衡量阻尼大小的无量纲参数：

ζ = c / (2 * √(m * k))

或者写成：

ζ = c / c_cr

其中c_cr = 2√(mk)是临界阻尼系数。

阻尼比的意义：

• ζ = 0：无阻尼，振动永不停止。
• 0 < ζ < 1：欠阻尼，实际结构都属此类。振动逐渐衰减。
• ζ = 1：临界阻尼，结构刚好不振动，直接回到平衡位置。
• ζ > 1：过阻尼，缓慢回到平衡位置，没有振动。

钢结构塔筒的阻尼比通常很小，大约在0.5%~2%之间。混凝土塔筒稍大，约2%~5%。

4.4 共振原理

共振，是结构工程师的噩梦，也是动力学里最核心的概念。

当外荷载的频率接近结构的固有频率时，即使荷载很小，结构也会产生巨大的振动响应。这就是共振。

用数学表达：对于简谐荷载F(t) = F₀ * sin(ωt)，结构的稳态响应幅值为：

X = F₀ / k * 1 / √[(1 - r²)² + (2ζr)²]

其中：

• r：频率比，r = ω / ωₙ。
• F₀/k：静位移，即把荷载当静力算的位移。
• 放大系数：D = 1 / √[(1 - r²)² + (2ζr)²]。

当r=1（即ω=ωₙ）时，放大系数达到最大值：

D_max = 1 / (2ζ)

如果阻尼比ζ=0.01，那么D_max = 50。也就是说，共振时的位移是静力位移的50倍！

你想想看，塔筒顶部如果装了风机，叶轮旋转频率（1P、3P）也会激励塔筒。如果1P频率恰好等于塔筒一阶频率，那后果不堪设想。这就是为什么大型风机塔筒的一阶频率通常设计在0.2~0.4Hz之间，避开叶轮旋转频率。

好了，这一章的内容就到这里。结构动力学是塔筒分析的根基，单自由度体系虽然简单，但它的概念贯穿整个课程。后面的多自由度体系、风荷载响应分析，都是在这个基础上扩展的。

---