粘弾性减衰材的周波数应答

分类:结构分析 | 统合版 2026-04-06
CAE visualization for viscoelastic damping theory - technical simulation diagram
粘弾性减衰材的周波数应答

粘弾性减衰材周波数应答的理论基础

粘弾性材料

🧑‍🎓

老师,粘弾性材料的减衰与结构减衰有什么不同?


🎓

结构减衰假设 $g$ 与周波数无关,但粘弾性材料的减衰强烈依赖周波数和温度。橡胶、聚合物、粘合剂和制振材是典型例子。


复杂弹性率

🎓

粘弾性材料的应力-应变关系:


$$ \sigma(\omega) = E^*(\omega) \varepsilon(\omega) $$
$$ E^*(\omega) = E'(\omega) + iE''(\omega) = E'(\omega)(1 + i\eta(\omega)) $$

  • $E'(\omega)$ — 储能弹性率(刚性。储存能量)
  • $E''(\omega)$ — 损耗弹性率(减衰。散逸能量)
  • $\eta(\omega) = E''/E'$ — 损失因子

🧑‍🎓

刚性和减衰都随周波数变化。


🎓

典型行为:

  • 低周波 — $E'$ 较低(柔软)、$\eta$ 中等
  • 玻璃转变区域 — $E'$ 急增、$\eta$ 达峰值(最大减衰)
  • 高周波 — $E'$ 较高(硬质)、$\eta$ 下降

温度-周波数等价规则(WLF规则)

🎓

粘弾性材料表现为温度上升时向低周波侧移位温度下降时向高周波侧移位Williams-Landel-Ferry(WLF)式描述这种等价性:


$$ \log a_T = -\frac{C_1(T-T_{ref})}{C_2 + (T-T_{ref})} $$

$a_T$ 为移位因子。由此可从单一温度的数据推估全温度的特性。


🧑‍🎓

通过DMA(动态粘弾性试验)在一个温度测量数据,就能用WLF式外推到其他温度吗?


🎓

建立主曲线。将基准温度下的DMA数据按WLF规则移位,获得宽频率范围的 $E'(\omega), \eta(\omega)$。然后输入FEM。


总结

🎓

要点:


  • $E^*(\omega) = E'(1+i\eta)$ — 周波数依存的复杂弹性率
  • $E', \eta$ 均依赖周波数和温度 — 玻璃转变处达峰值
  • WLF式实现温度-周波数等价 — 主曲线的构建
  • 直接法(周波数应答)输入FEM — 模态法处理困难

Coffee Break 闲话

粘弾性的记忆效应由Maxwell在1867年发现

粘弾性材料"应变依赖于应力历史"的记忆效应由詹姆斯·克拉克·麦克斯韦在1867年气体粘性研究中发现。后来Kelvin-Voigt(1890年)提出并联模型,Maxwell串联+Kelvin并联组合形成的"Zener(标准线性固体)模型"成为结构减衰FEM实现的基础。3M阻尼片的特性评估至今仍使用这一模型。

粘弾性减衰材周波数应答的数值计算手法

FEM中的粘弾性材料

🧑‍🎓

如何在FEM中处理粘弾性材料?


🎓

用Prony级数对 $E'(\omega), E''(\omega)$ 建模:


$$ E'(\omega) = E_\infty + \sum_{i=1}^{N} E_i \frac{\omega^2 \tau_i^2}{1 + \omega^2 \tau_i^2} $$
$$ E''(\omega) = \sum_{i=1}^{N} E_i \frac{\omega \tau_i}{1 + \omega^2 \tau_i^2} $$

$E_i, \tau_i$ 为Prony级数参数。$N = 5 \sim 15$ 项覆盖宽频率范围。


Abaqus

```

*VISCOELASTIC, FREQUENCY=PRONY

g_1, k_1, tau_1

g_2, k_2, tau_2

...

```

Nastran

```

MAT1, 1, ...

TABLEM1, 100, ...

$ 周波数依存的E'和η通过表格定义

```

Nastran可直接通过表格(周波数 vs. E', η)输入。

🧑‍🎓

Abaqus的Prony级数能直接从DMA数据拟合吗?


🎓

Abaqus/CAE的*VISCOELASTIC, TEST DATA选项支持输入实测的$E'(\omega), E''(\omega)$,并自动拟合Prony级数。非常便利。


制振材(CLD)的分析

🎓

CLD(约束层阻尼)的分析:

1. 建模基板(钢板等)+粘弾性层+约束层的层叠结构

2. 粘弾性层设置周波数依存的Prony级数

3. 直接法周波数应答分析评估减衰效果

4. 确认FRF峰值降低


总结

🎓
  • 用Prony级数建模 — N项 $E_i, \tau_i$ 覆盖宽频率范围
  • Abaqus的TEST DATA自动拟合 — 直接从DMA数据
  • 必须用直接法 — 周波数依存材料模态法精度不足
  • CLD分析 — 基板+粘弾性+约束层层叠评估制振效果

  • Coffee Break 闲话

    时间-温度换算规则将数据压缩到1/100

    粘弾性材料的损失因子随温度和周波数的组合变化。利用Williams-Landel-Ferry(WLF)式(1955年)的时间-温度换算规则,仅需基准温度Tr下的测量数据就能预测任意温度、周波数条件下的特性。车内内装材(乙丙橡胶系)的周波数特性可用20℃基准的单条主曲线表示-40℃~100℃、1Hz~10kHz的全部条件。

    粘弾性减衰材周波数应答的实务应用

    粘弾性减衰的实务

    🎓

    制振材设计在汽车NVH中是最重要的应用。


    制振材的设计流程

    🎓

    1. DMA试验测量粘弾性材料的$E'(\omega, T), \eta(\omega, T)$

    2. 用WLF规则建立主曲线

    3. 拟合Prony级数

    4. FEM优化制振材配置和厚度

    5. 确认FRF共振峰值降低


    实务检查清单

    🎓
    • [ ] DMA数据是否覆盖目标温度和周波数范围
    • [ ] Prony级数拟合是否良好(与实测 $E', \eta$ 对比)
    • [ ] 是否在全使用温度范围内评估(注意橡胶转变温度附近)
    • [ ] 是否用直接法计算周波数应答(模态法精度不足)
    • [ ] 制振材配置是否在最大应变位置(曲折最大应变点贴附)

    • 🧑‍🎓

      制振材要贴在"最大应变的位置"?


      🎓

      制振材通过剪切变形散逸能量。基板弯曲应变最大的位置(固定端附近)贴附制振材效果最佳。应变较小的位置(自由端附近)贴附效果有限。


      Coffee Break 闲话

      3M制振片将新干线噪声降低5dB

      3M公司的制振材"Dynamat"系列为钢板积层粘弾性聚合物(丙烯酸系)结构,通过熔化100~1000Hz曲折波将其转换为热能。JR东日本E5系新干线(2011年)车底钢板贴附3M Dynamat相当品,将轨道噪声源车内固体音降低约5dB。仅厚2mm,面密度2.4kg/m²。由于成本优势,已广泛用于量产车。

      粘弾性减衰材周波数应答的软件比较

      粘弾性分析工具

      🎓
      工具特点
      Abaqus *VISCOELASTICProny级数+TEST DATA。最灵活
      Nastran TABLEM1周波数表直接输入
      Ansys支持Prony级数
      Actran声学分析中粘弾性支持
      Comet/VA OneSEA中制振材优化

      选定指南

      🎓
      • 制振材设计(FEM → Abaqus直接法 + *VISCOELASTIC
      • NVH制振材优化Nastran + 汽车OEM内部工具
      • SEA级制振材设计 → VA One

      • Coffee Break 闲话

        日东电工与ESI联动的粘弾性材料数据库

        日东电工(胶带系最大企业)与ESI Group(VA One)合作,将主要制振材(Rethertex系列等)的动态特性数据库集成到VA One中。这使设计者能在同一软件中完成材料选择和分析。3M也向全球公开发布自家制振材的Abaqus Material Card,大幅减少材料数据的输入工时。截至2022年,主要制振材30种以上已数据库化。

        粘弾性减衰材周波数应答的先端研究

        分数导数粘弾性模型

        🎓

        作为Prony级数的替代方案,分数导数模型研究在进行中。用少量参数(3~4个)覆盖宽频率范围。FEM实现仍在发展中。


        拓扑优化制振材配置

        🎓

        制振材贴在哪里用拓扑优化自动决定的研究。目的函数:共振峰值最小化。设计变量:制振材配置(0/1)。


        总结

        🎓
        • 分数导数模型 — 少参数宽频带
        • 制振材配置拓扑优化 — 自动确定最优贴附位置

        • Coffee Break 闲话

          约束型制振比自由型损失因子高5倍

          制振材结构分为①自由层(单体贴附)和②约束层(金属板夹心)两类。约束型通过集中剪切应变,损失因子高5~10倍。Nashif(1985年)实验数据表明1mm厚约束型制振处理比同质量自由型振动低减效果高3倍。汽车仪表板采用约束型为主流。设计可将3M、日东电工的CAE数据库直接导入Abaqus/Ansys。

          粘弾性减衰材周波数应答的故障处理

          Prony级数拟合不佳

          🎓
          • 增加项数(5→10→15)
          • 调整松弛时间 $\tau_i$ 范围以覆盖目标周波数
          • 使用Abaqus的TEST DATA自动拟合

          • 忽略温度依存性

            🎓

            粘弾性随温度剧烈变化。玻璃转变温度附近$E'$变化10倍,$\eta$达峰值。应在全使用温度范围内评估。


            在模态法中使用粘弾性

            🎓

            模态法在固有模式(与材料特性无关的基函数)上展开,无法准确处理周波数依存材料。应切换为直接法。


            总结

            🎓
            • 拟合不良 → 增加项数,调整 $\tau_i$ 范围
            • 温度依存性 → 全使用温度范围评估
            • 模态法误用 → 切换为直接法
            • 粘弾性"同时依赖温度和周波数" — 两者都不能忽视

            • Coffee Break 闲话

              夏冬季减衰效果相差3倍

              粘弾性材料的温度依存性强,同一制振片在0℃(冬季)和40℃(夏季)的损失因子相差2~3倍。北海道铁路车辆夏季试验通过的制振对策在冬季运行中失效的事例曾在2008年报告。对策是在全使用温度范围内获取主曲线,在最恶劣温度条件下确认设计余裕。ANSYS的Viscoelastic Material Model中设置温度表是标准流程。

              相关模拟器

              通过本领域的交互式模拟器体验理论

              模拟器列表

              相关领域

              热分析制造工艺分析V&V·质量保证
              本文评价
              感谢您的回答!
              有参考
              价值
              需要
              更多细节
              报告
              错误
              有参考价值
              0
              需要更多细节
              0
              报告错误
              0
              由NovaSolver贡献者撰写
              匿名工程师与AI — 网站地图
              查看档案