结构阻尼（滞后阻尼）——复数刚度、损耗系数、FEM实现

分类: 構造解析 > 周波数応答解析 | 更新 2026-04-11
Structural damping complex stiffness model visualization showing hysteresis loop and frequency response function with loss factor
構造減衰モデル — 複素剛性 $K^* = K(1+i\eta)$ による周波数応答と損失係数の可視化

理论与物理

复刚度与损失系数的定义

🧑‍🎓

老师，振动分析中出现了“结构阻尼”，它和普通的阻尼有什么区别呢？是和阻尼器之类的东西不一样吧？

🎓

问得好。阻尼器（粘性阻尼器）的阻尼力与速度成正比——$F_d = c\dot{x}$。这就是粘性阻尼。另一方面，结构阻尼（迟滞阻尼）表示材料内部微观摩擦或滑移引起的能量耗散。在数学上用复刚度表示：

$$ K^* = K(1 + i\eta) $$

这里的 $\eta$ 就是损失系数（loss factor）。实部 $K$ 表示通常的弹性刚度，虚部 $K\eta$ 表示能量耗散的大小。

🧑‍🎓

用复数表示刚度，在物理上是什么意思呢？很难想象虚数的弹簧……

🎓

可以这样理解。假设简谐振动 $x = X e^{i\omega t}$，则恢复力为：

$$ F = K^* x = K(1+i\eta) X e^{i\omega t} $$

实部是与位移同相的弹性恢复力，虚部是比位移超前90°的力——即对应于速度分量的耗散力。与粘性阻尼力 $c\dot{x} = ci\omega x$ 比较，结构阻尼的耗散力与 $\omega$ 无关。这是决定性的区别。

与粘性阻尼的本质区别

🧑‍🎓

与频率无关具体会产生什么差异呢？比如在汽车振动分析中。

🎓

在粘性阻尼中，等效阻尼比与频率成正比：

$$ \zeta_{\text{viscous}} = \frac{c\omega}{2k} $$

所以20 Hz和200 Hz时阻尼比会相差10倍之多。在汽车仪表板振动分析处理50〜500 Hz的宽频带时，如果只用粘性阻尼，会导致低频阻尼不足，高频阻尼过大。

如果是结构阻尼，等效阻尼比为：

$$ \zeta_{\text{eq}} = \frac{\eta}{2} $$

在所有频率下保持恒定。实际的金属或FRP结构的阻尼特性，与结构阻尼模型要接近得多。所以在频率响应分析（NVH分析）中，结构阻尼是标准。

🧑‍🎓

那么 $\eta$ 和 $\zeta$ 的关系就是 $\eta \approx 2\zeta$ 对吧。如果知道模态分析的阻尼比，可以转换吗？

🎓

没错。作为共振频率下的等效关系：

$$ \eta = 2\zeta \quad (\text{在共振点附近严格成立}) $$

例如，如果实验模态分析识别出 $\zeta = 1\%$，则设定 $\eta = 0.02$。但需要注意的是，这个等效关系仅在共振频率下严格成立。在远离共振的频率下，粘性阻尼和结构阻尼的响应会有所不同。

能量耗散与迟滞回线

🧑‍🎓

之所以称为“迟滞阻尼”，是因为与应力-应变曲线的迟滞回线有关吗？

🎓

正是如此。对材料施加重复载荷时，应力-应变曲线会描绘出一个闭合的回线。这个回线的面积对应于能量耗散量 $\Delta W$。损失系数可以定义为每个周期的耗散能量与最大弹性能量之比：

$$ \eta = \frac{\Delta W}{2\pi W_{\max}} $$

其中 $W_{\max} = \frac{1}{2}K X^2$ 是最大弹性应变能。例如，钢材的回线面积非常小，所以 $\eta \approx 0.001$〜$0.01$；橡胶的回线面积大，所以 $\eta \approx 0.2$〜$1.0$。字面意思就是“迟滞的大小”决定了阻尼的强弱。

半功率带宽法识别损失系数

🧑‍🎓

实验中如何测量损失系数呢？可以通过冲击锤试验求得吗？

🎓

最常用的是半功率带宽法。从FRF的共振峰读取损失系数。步骤如下：

  1. 确定FRF的共振频率 $f_n$ 和最大振幅 $|H|_{\max}$
  2. 找出振幅为 $|H|_{\max}/\sqrt{2}$（= -3 dB）的两个频率 $f_1, f_2$
  3. 计算损失系数：
$$ \eta = \frac{f_2 - f_1}{f_n} = \frac{\Delta f}{f_n} \approx 2\zeta $$

例如，冲击锤试验中，共振在200 Hz，-3 dB带宽为2 Hz，则 $\eta = 2/200 = 0.01$。这相当于 $\zeta = 0.5\%$。作为钢结构，这是焊接结构的典型值。

🧑‍🎓

半功率带宽法在模态密集的情况下似乎效果不好吧？

🎓

很敏锐。模态密集时，-3 dB线会与相邻的峰值重叠，无法准确读取。这种情况下可以使用圆拟合法（在奈奎斯特图上进行圆弧拟合）或曲线拟合法（RFP法、LSCE法等）。特别是像CFRP层压板这样阻尼大、模态容易重叠的材料，应避免使用简单的半功率带宽法。

Coffee Break 闲谈

复刚度诞生于飞机颤振分析

用复刚度 $K^* = K(1+i\eta)$ 表示结构阻尼的模型，诞生于1930年代的飞机机翼颤振研究。为了求解机翼弹性变形与空气动力的耦合振动，需要一种简洁的方法来表达材料的耗散特性。1960年，加州理工学院的 T.K. Caughey 严格整理了复刚度与粘性阻尼的等效关系，确立了 $\zeta = \eta/2$ 这个实用公式。至今仍作为Nastran的GE参数被沿用。

数值解法与实现

直接法与模态法中的结构阻尼

🧑‍🎓

结构阻尼是在频率响应分析中使用的吧？是用直接法还是模态法呢？

🎓

两种方法都可以使用。直接法（Nastran SOL 108）中，将复刚度直接代入运动方程：

$$ \left[-\omega^2 M + i\omega C + K(1+i\eta)\right] X = F $$

模态法（Nastran SOL 111）中，先展开到固有模态，再对每个模态应用阻尼：

$$ -\omega^2 q_r + i\omega \cdot 2\zeta_r \omega_r q_r + \omega_r^2(1+i\eta_r) q_r = \phi_r^T F $$

这里 $q_r$ 是模态坐标，$\phi_r$ 是模态向量。在实际工作中，模态法占绝大多数。因为计算成本只有直接法的1/10〜1/100。

🧑‍🎓

模态法可以对每个模态设置不同的 $\eta_r$ 吗？实验中各模态的阻尼比不同是很常见的吧？

🎓

使用Nastran的TABDMP1卡片，可以定义阻尼相对于频率的表格。它会根据各模态的固有频率自动插值 $\eta$。例如，如果实验模态分析识别出一阶模态 $\zeta_1=0.5\%$，二阶模态 $\zeta_2=1.2\%$，则在表格中输入 $\eta_1=0.01$，$\eta_2=0.024$。

频率响应函数（FRF）的公式化

🧑‍🎓

请告诉我包含结构阻尼的FRF公式。我想知道共振点的振幅如何变化。

🎓

单自由度系统包含结构阻尼的FRF（导纳）如下：

$$ H(\omega) = \frac{1}{K(1+i\eta) - \omega^2 M} = \frac{1}{K\left[(1-r^2) + i\eta\right]} $$

其中 $r = \omega/\omega_n$（频率比）。在共振点 $r=1$ 处：

$$ |H(\omega_n)| = \frac{1}{K\eta} $$

也就是说共振振幅与 $1/(K\eta)$ 成正比，与损失系数成反比。如果 $\eta$ 的数值错一位，共振振幅就会差10倍，所以设置值的位数一定要确认。在汽车转向柱振动分析中，把 $\eta=0.02$ 误输入为 $\eta=0.002$，导致FRF峰值大了10倍而引起大麻烦——这是实际发生过的案例。

向瑞利阻尼的转换方法

🧑‍🎓

瞬态响应分析中不能使用结构阻尼对吧？转换为瑞利阻尼时，如何确定 $\alpha$ 和 $\beta$ 呢？

🎓

瑞利阻尼为 $C = \alpha M + \beta K$，阻尼比为：

$$ \zeta(\omega) = \frac{\alpha}{2\omega} + \frac{\beta \omega}{2} $$

结构阻尼的 $\zeta = \eta/2$ 是恒定的，所以联立方程，使在两个频率 $\omega_1, \omega_2$ 处满足 $\zeta(\omega_i) = \eta/2$：

$$ \alpha = \eta \cdot \frac{\omega_1 \omega_2}{\omega_1 + \omega_2}, \quad \beta = \frac{\eta}{\omega_1 + \omega_2} $$

例如，$\eta = 0.02$，关注频带为50〜300 Hz（$\omega_1 = 100\pi, \omega_2 = 600\pi$），则计算 $\alpha, \beta$ 使在这两点处 $\zeta = 1\%$。需要注意的是，这两点之间的阻尼会偏小，之外的阻尼会偏大。频带的选择会显著影响结果，需要谨慎。

因果律约束与在时域中的处理

🧑‍🎓

请再详细解释一下为什么结构阻尼不能在时域中使用。“违反因果律”具体是什么意思？

🎓

将复刚度 $K(1+i\eta)$ 进行傅里叶逆变换回到时域，会发现脉冲响应函数 $h(t)$ 在 $t < 0$ 时不等于零。也就是说输入还没到来，响应就已经出现——因果律（原因→结果的时间顺序）被破坏了。

这是因为损失系数 $\eta$ 在所有频率下恒定这一假设本身，违反了Kramers-Kronig关系（满足因果律所需的实部与虚部约束条件）。物理上严格的阻尼模型必然具有频率依赖性。结构阻尼应作为“频域中的近似模型”来使用。

Coffee Break 闲谈

半功率带宽法已持续使用60多年

半功率带宽法是通过FRF峰值的 $1/\sqrt{2}$ 振幅对应的两个频率，计算损失系数 $\eta = \Delta f / f_n \approx 2\zeta$ 的方法。自D.J. Ewins在1960年代将其体系化以来，至今仍作为振动试验的标准方法使用。钢材的 $\eta$ 约为0.001〜0.01，CFRP复合材料约为0.005〜0.02，减振钢板约为0.05〜0.2，不同材料间存在两个数量级以上的差异。

实践指南

材料/结构对应的典型损失系数

🧑‍🎓

损失系数 $\eta$ 的值一般是多少？分析输入时，应该参考什么？

🎓

我来总结一下典型值。首先从材料本身的损失系数开始：

材料损失系数 $\eta$等效阻尼比 $\zeta$备注
纯铝0.0001〜0.0010.005〜0.05%取决于纯度、晶体结构
铝合金（A6061等）0.001〜0.0050.05〜0.25%航空结构中常见
软钢（SS400等）0.001〜0.0060.05〜0.3%材料本身的值
钢结构（焊接）0.005〜0.0150.25〜0.75%包含连接部位的摩擦
钢结构（螺栓连接）0.01〜0.050.5〜2.5%螺栓面的微小滑移起主导作用
CFRP层压板0.005〜0.020.25〜1%取决于铺层构成、基体
混凝土0.02〜0.11〜5%裂纹状态变化时差异大
橡胶（天然橡胶）0.1〜0.55〜25%强烈依赖于温度、频率
减振橡胶/丁基橡胶0.3〜1.515〜75%用作减振材料
减振钢板（夹层结构）0.05〜0.32.5〜15%约束型减振材料
🧑‍🎓

螺栓连接的结构比材料本身的阻尼要大呢。为什么？

🎓

注意到了关键点。根据Beards（1983）的研究，机械结构总阻尼的60〜90%源于连接部位的摩擦。螺栓连接面在振动时会发生微小滑移，通过库仑摩擦耗散能量。所以即使是相同的钢材，螺栓连接结构（$\eta \approx 0.03$）的阻尼也比焊接结构（$\eta \approx 0.01$）大近3倍。FEM中只输入材料阻尼而与实测不符，大多是因为这个原因。

汽车NVH分析中的应用案例

🧑‍🎓

在汽车NVH分析中使用结构阻尼时，具体如何设置？

🎓

汽车车身（BIW: Body In White）的频率响应分析是结构阻尼的典型应用案例。典型的设置总结如下：

  • 钢板车身面板：$\eta = 0.01$〜$0.02$（包含点焊部位的摩擦）
  • 仪表板减振材料粘贴部位：$\eta = 0.02$〜$0.05$（约束型减振材料的效果）
  • 橡胶悬置（发动机/副车架）：$\eta = 0.1$〜$0.3$
  • 玻璃：$\eta = 0.002$〜$0.005$

在Nastran的SOL 111中进行模态法分析，使用MAT1的GE字段为每种材料设置 $\eta$ 是标准做法。用PARAM, G设置全局统一的值是比较粗糙的方法，如果需要精度，则必须按材料分别设置。

航空航天/工厂管道中的应用

🧑‍🎓

航空航天领域呢？和汽车不同吗？

🎓

在航空航天领域，由于阻尼较小，设置值更为关键。举几个具体例子：

  • 人造卫星面板（CFRP/铝蜂窝）：$\eta = 0.005$〜$0.01$。用于运载火箭的振动环境分析。如果 $\eta$ 设置过大，会低估峰值应力，有在轨道上损坏的危险
  • 喷气发动机叶片：$\eta = 0.001$〜$0.003$。高周疲劳评估必须准确预测共振振幅
  • 工厂管道系统：$\eta = 0.01$〜$0.03$（包含支架、卡箍的摩擦）。ASME/JSME标准推荐管道系统的设计阻尼比为 $\zeta = 1$〜$2\%$（$\eta = 0.02$〜$0.04$）

实践检查清单

🎓

总结一下结构阻尼设置时需要检查的项目：