隔振设计与传递率

分类: 構造解析 | 综合版 2026-04-06
CAE visualization for vibration isolation theory - technical simulation diagram
防振設計と伝達率

理论与物理

防振(振动隔离)是什么

🧑‍🎓

老师,“防振”是停止振动吗?


🎓

不是“停止”振动,而是减少传递。在振源和被振体之间插入弹簧(防振橡胶、防振支架)来抑制振动的传递。


传递率

🎓

传递率(transmissibility)$T$ 是输出与输入的比值:


$$ T(\omega) = \frac{|X_{out}|}{|X_{in}|} = \sqrt{\frac{1 + (2\zeta r)^2}{(1-r^2)^2 + (2\zeta r)^2}} $$

其中 $r = \omega / \omega_n$(频率比),$\zeta$ 是阻尼比。


🧑‍🎓

$r = 1$(共振)时传递率达到峰值,$r > \sqrt{2}$ 时 $T < 1$(防振效果),对吧。


🎓

完全正确。$r > \sqrt{2}$(即 $f > \sqrt{2} f_n$)的区域是防振区域。在此区域内输出小于输入。


🎓

设计要点:

  • 降低 $f_n$ — 防振区域变宽。使支架更柔软
  • 但过于柔软则静挠度会变大 — 存在实际约束
  • 若运行过程会通过共振点,则需要阻尼 — 用 $\zeta$ 抑制峰值

防振支架的选型

🎓
支架类型弹簧常数阻尼用途
橡胶支架中等中等($\zeta$ 5〜15%)发动机支架、设备支架
螺旋弹簧低($\zeta$ < 1%)精密设备防振
空气弹簧非常低半导体制造设备
钢丝支架中等中等(摩擦阻尼)军用设备
🧑‍🎓

精密设备用空气弹簧…能将 $f_n$ 降低到约 0.5 Hz 左右呢。


🎓

空气弹簧的 $f_n = 0.5 \sim 2$ Hz。几乎可以隔绝所有外部振动。半导体曝光设备和激光设备中空气弹簧是标准配置。


FEM中的防振设计

🎓

FEM中的防振设计:


1. 构建设备+支架+基础的FEM模型

2. 将支架建模为弹簧单元(+阻尼器)

3. 对基础施加输入振动(频率响应或时程)

4. 计算设备的响应(位移、加速度)

5. 绘制传递率 $T = |X_{out}| / |X_{in}|$

6. 确认 $T < T_{target}$


总结

🧑‍🎓

我来整理一下防振设计和传递率。


🎓

要点:


  • 传递率 $T$ 是设计的核心指标 — $T < 1$ 表示有防振效果
  • $f > \sqrt{2} f_n$ 为防振区域 — $f_n$ 越低效果越好
  • 共振峰值需用阻尼抑制 — 合理设置 $\zeta$
  • 支架的选型 — 橡胶、螺旋、空气、钢丝
  • 用FEM计算传递率 — 弹簧单元+频率响应分析

Coffee Break 闲谈

防振的黄金比例:固有频率为激励频率的1/3以下

防振设计的基本法则是“支架的固有频率 fn ≤ 激励频率 f0 / √2 ≈ f0 × 0.7以下”,若将 fn 降至 f0/3(3倍法则),则传递率可降至1/8以下。该法则源自J.P. Den Hartog于1934年在其著作《Mechanical Vibrations》中展示的传递率曲线。电子显微镜(SEM/TEM)的安装中,fn ≤ 1Hz的超低刚度空气支架是标准配置,用以防止百万倍放大倍率下的图像模糊。

各项的物理意义
  • 惯性项(质量项):$\rho \ddot{u}$,即“质量×加速度”。您是否有过急刹车时身体向前冲的经验?那种“被带走的感觉”正是惯性力。物体越重越难启动,一旦启动也越难停止。地震时建筑物摇晃,也是因为地面突然移动,而建筑物的质量“被落下”。静力分析中此项设为零,这是基于“缓慢施力,加速度可忽略”的假设。对于冲击载荷或振动问题,此项绝不能省略。
  • 刚度项(弹性恢复力):$Ku$ 或 $\nabla \cdot \sigma$。拉伸弹簧时能感觉到“试图恢复原状的力”吧?那就是胡克定律 $F=kx$,也是刚度项的本质。那么提问——用相同的力拉伸铁棒和橡皮筋,哪个伸长更多?当然是橡皮筋。这种“难以伸长的程度”就是杨氏模量 $E$,它决定了刚度。常见的误解是:“刚度高=强度高”。刚度是“抵抗变形的能力”,强度是“抵抗破坏的能力”,两者是不同的概念。
  • 外力项(载荷项):体积力 $f_b$(如重力)和表面力 $f_s$(如压力、接触力)。可以这样理解——桥上卡车的重量是“作用于整个内部体积的力”(体积力),轮胎压路面的力是“仅作用于表面的力”(表面力)。风压、水压、螺栓预紧力…这些都是外力。这里容易犯的错误是:弄错载荷方向。本想施加“拉力”却成了“压力”——听起来像笑话,但在3D空间中坐标系发生旋转时确实会发生。
  • 阻尼项:瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他弦。声音会一直持续吗?不,会逐渐变小。这是因为振动能量通过空气阻力和弦的内部摩擦转化为热能。汽车的减震器也是同样原理——故意吸收振动能量以改善乘坐舒适性。如果阻尼为零会怎样?建筑物在地震后会一直摇晃不停。现实中不会如此,因此合理设置阻尼很重要。
假设条件与适用范围
  • 连续体假设:将材料视为连续介质,忽略微观不均匀性
  • 小变形假设(线性分析时):变形相对于初始尺寸足够小,应力-应变关系呈线性
  • 各向同性材料(除非特别指定):材料特性不依赖于方向(各向异性材料需另行定义张量)
  • 准静态假设(静力分析时):忽略惯性力·阻尼力,仅考虑外力与内力的平衡
  • 不适用的情形:大变形·大旋转问题需要考虑几何非线性。塑性·蠕变等非线性材料行为需要扩展本构关系
量纲分析与单位制
变量SI单位注意事项·换算备忘
位移 $u$m(米)以mm输入时,载荷·弹性模量也需统一为MPa/N系
应力 $\sigma$Pa(帕斯卡)= N/m²MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致
应变 $\varepsilon$无量纲(m/m)注意工程应变与对数应变的区别(大变形时)
弹性模量 $E$Pa钢:约210 GPa,铝:约70 GPa。注意温度依赖性
密度 $\rho$kg/m³mm系中为tonne/mm³(钢约为 10⁻⁹ tonne/mm³)
力 $F$N(牛顿)mm系和m系均统一用N

数值解法与实现

FEM中的防振支架建模

🧑‍🎓

防振支架在FEM中如何建模?


🎓

弹簧单元+阻尼器(粘性单元)并联表示。设置三个方向的弹簧常数和阻尼系数。


Nastran

```

CBUSH, 100, 200, 1000, 2000 $ 衬套单元

PBUSH, 200, K, 1000., 1000., 5000. $ kx, ky, kz

, B, 10., 10., 50. $ cx, cy, cz

```

Abaqus

```

*CONNECTOR SECTION, BEHAVIOR=mount

BUSHING,

*CONNECTOR BEHAVIOR, NAME=mount

*CONNECTOR ELASTICITY

1000., 1000., 5000.

*CONNECTOR DAMPING

10., 10., 50.

```

橡胶支架的非线性特性

🎓

橡胶支架具有频率依赖的刚度和阻尼(粘弹性特性)。


  • 静态刚度 — 低频下的弹簧常数
  • 动态刚度 — 高频下的弹簧常数(比静态高20〜50%)
  • 损耗系数 $\eta$ — 频率依赖的阻尼

🧑‍🎓

动态刚度比静态刚度高?


🎓

橡胶是粘弹性材料,振动频率越高会变得越硬。若将静态试验得到的弹簧常数直接用于动态分析,会高估防振效果。应通过动态试验(DMA:动态力学分析)测量频率依赖特性。


传递率的计算

🎓

```

$ 传递率 = 输出点的加速度 / 输入点的加速度

T(f) = |a_output(f)| / |a_input(f)|

```


通过FEM的频率响应分析输出输入点和输出点的加速度,然后取比值。


总结

🧑‍🎓

我来整理一下防振设计的数值方法。


🎓

要点:


  • 用 CBUSH(Nastran)/ CONNECTOR(Abaqus)表示支架
  • 设置三个方向的弹簧常数+阻尼 — 也可设置各向异性
  • 橡胶的动态刚度高于静态刚度 — 使用DMA数据
  • 传递率 = 输出/输入的比值 — 从频率响应分析中计算得出

Coffee Break 闲谈

空气弹簧的固有频率随配管长度变化

空气弹簧的固有频率与封闭空气体积V的-1/2次方成正比,因此通过增加辅助储气罐体积可以降低固有频率(最低可达0.5〜1Hz左右)。仅用螺旋弹簧难以实现3Hz以下的固有频率,因此半导体制造设备(如ASML TWINSCAN等)全部采用空气弹簧。可将制造车间的地板振动(主要为2〜10Hz)衰减至1/100以下。

线性单元(一阶单元)

节点间线性插值。计算成本低,但应力精度低。注意剪切自锁(可通过减缩积分或B-bar法缓解)。

二阶单元(带中间节点)

可表现曲线变形。应力精度大幅提高,但自由度约增加2〜3倍。推荐:应力评估重要的场合。

完全积分 vs 减缩积分

完全积分:有过约束(自锁)风险。减缩积分:有沙漏模式(零能量模式)风险。根据情况选择合适方法。

自适应网格

基于误差指标(如ZZ估计量)的自动细化。高效提高应力集中区域的精度。有h法(单元细分)和p法(增加阶次)两种。

牛顿-拉弗森法

非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性,但计算成本高。

修正牛顿-拉弗森法

切线刚度矩阵使用初始值或每隔数次迭代更新。每次迭代成本低,但收敛速度为线性。

收敛判定准则

力残差范数: $||R|| / ||F_{ext}|| < \epsilon$(通常 $\epsilon = 10^{-3}$〜$10^{-6}$)。位移增量范数: $||\Delta u|| / ||u|| < \epsilon$。能量范数: $\Delta u \cdot R < \epsilon$

载荷增量法

不一次性施加全部载荷,而是分小步增加。弧长法(Riks法)可以超越载荷-位移关系的极值点进行追踪。

直接法 vs 迭代法的比喻

直接法是“用笔算精确解联立方程”的方法——可靠但大规模问题耗时过长。迭代法是“反复猜测逼近正确答案”的方法——最初答案粗略,但每次迭代精度提高。就像查字典时,从第一页开始顺序查找(直接法)不如先估计大致位置翻开,再前后调整(迭代法)来得高效,原理相同。

网格阶次与精度的关系

一阶单元如同“用直尺近似曲线”——用直线折线表现,精度有限。二阶单元如同“柔性曲线”——可以表现曲线变化,即使网格密度相同,精度也显著提高。但每个单元的计算成本增加,需根据总体的成本效益来判断。

实践指南

防振设计的实务

🧑‍🎓

防振设计在实际工作中如何进行?


STEP 1 振动环境的确定

関連シミュレーター

この分野のインタラクティブシミュレーターで理論を体感しよう

シミュレーター一覧

関連する分野

熱解析製造プロセス解析V&V・品質保証
この記事の評価
ご回答ありがとうございます!
参考に
なった
もっと
詳しく
誤りを
報告
参考になった
0
もっと詳しく
0
誤りを報告
0
Written by NovaSolver Contributors
Anonymous Engineers & AI — サイトマップ
关于作者