有效质量比与参与系数
理论与物理
什么是有效质量
老师，什么是“有效质量”（模态有效质量）？
有效质量是衡量各振动模态对外部惯性力（地震、加速度）响应程度的指标。有效质量大的模态主导结构的动态响应。
并不是所有模态都以相同方式响应吗？
是的。例如悬臂梁的一阶弯曲模态具有较大的有效质量，但二阶模态的有效质量比一阶小。有效质量量化了各模态的“重要性”。
数学定义
第 $i$ 阶模态在 $x$ 方向的有效质量：
其中 $\{1_x\}$ 是 $x$ 方向的单位向量（所有节点的 $x$ 自由度为1，其他为零）。
分子是模态形状与“整体沿 $x$ 方向运动”模式的内积…是模态与惯性力的“相似度”吧。
理解得很完美。整体沿单一方向运动的模态（平动模态）具有较大的有效质量，而局部模态或旋转模态的有效质量较小。
有效质量比
有效质量的总和等于总质量：
通常用各模态的有效质量比（有效质量/总质量）来表示。
所有模态的有效质量比总和是100%对吧。
是的。在实际工程中，通常要求“累积有效质量比达到90%所需的模态数”。这是确定所需模态数的判断标准。
参与系数
“参与系数”（模态参与因子）是什么？
相当于有效质量分子的平方根：
在Nastran中使用 PARAM,EFFMASS 输出有效质量。在Abaqus中使用History Output中的 MODAL EFFECTIVE MASS 输出。
实际应用方法
| 用途 | 有效质量的用法 |
|---|---|
| 地震响应分析的模态数确定 | 各方向覆盖90%所需的模态数 |
| 主导模态的识别 | 有效质量最大的模态是主要响应 |
| 模态叠加法的精度确认 | 通过有效质量的覆盖率判断模态数的合理性 |
| 质量分布的确认 | 有效质量为零的模态是对称模态（不响应非对称输入） |
地震响应分析中，“覆盖90%”是基准对吧。
建筑基准法和欧洲规范8中规定了模态叠加法的有效质量比需达到90%。满足此条件的模态数就是最低限度的必要模态数。
总结
我来整理一下有效质量比和参与系数。
要点：
- 有效质量 = 衡量模态对惯性力响应程度的指标
- 所有模态的总和 = 总质量 —— 各模态的“重要性比率”
- 覆盖90%是实际工程基准 —— 判断所需模态数
- 有效质量大的模态是主导模态 —— 设计中应关注的模态
- 各方向的有效质量不同 —— 需分别检查 $x, y, z$ 各方向
一看有效质量就能明白“哪个模态重要”。不仅要看固有频率，也要养成看有效质量的习惯。
固有频率和模态形状之后，下一个要看的就是有效质量。这三者（频率、模态形状、有效质量）是振动分析的三大支柱。
有效质量分数的“90%规则”
模态有效质量分数用于确认模态分析中是否提取了足够数量的模态。ASME和欧洲抗震规范要求计算到主要三个方向各自的累积有效质量达到总质量90%以上的模态数。这被称为“90%规则”，是从1970年代NRC（美国核管理委员会）规范开始确立的惯例。
各项的物理意义
- 惯性项（质量项）：$\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系用N，m系也用N统一 |
数值解法与实现
有效质量的输出方法
在FEM中如何输出有效质量？
Nastran
```
PARAM, EFFMASS, YES
```
f06文件中会输出各模态的有效质量比和累积值。
Abaqus
```
*OUTPUT, HISTORY
*MODAL DYNAMIC
PARTICIPATION FACTOR, EFFECTIVE MODAL MASS
```
Ansys
```
MODOPT, LANB, 50
MXPAND, 50, , , YES ! 展开模态形状以输出有效质量
```
或在Workbench的“Effective Mass Summary”中确认。
如何读取Nastran的f06文件中的有效质量表？
```
MODE FREQ(HZ) T1 FRACTION T2 FRACTION T3 FRACTION R1 FRACTION ...
1 15.2 0.4521 0.0000 0.0000 0.0000
2 23.8 0.0000 0.3856 0.0000 0.0000
3 45.1 0.1823 0.0000 0.0000 0.0012
...
TOTAL 0.9234 0.9156 0.8912 ...
```
T1, T2, T3 是 $x, y, z$ 方向的平动有效质量比。R1, R2, R3 是旋转。TOTAL行是累积值。
$x$ 方向的累积值为0.92（92%），意味着50个模态覆盖了 $x$ 方向响应的92%，对吧。
是的。如果未达到90%，则需要增加模态数。
有效质量小的模态的解释
有很多有效质量几乎为零的模态。
有效质量为零并不意味着可以忽略，对吧。
它们对地震响应没有贡献，但与机械振动（旋转体不平衡等）的共振有关。根据用途，“重要的模态”是不同的。
总结
我来整理一下有效质量的数值方法。
要点：
- 使用PARAM,EFFMASS（Nastran）可轻松输出 —— 检查f06文件中的表格
- 分别检查各方向的有效质量比 —— $x, y, z, \theta_x, \theta_y, \theta_z$
- 累积90%是模态数的判断基准 —— 设计规范的要求
- 有效质量为零的模态 = 局部模态 —— 对整体响应无贡献但需注意局部共振
有效质量的计算公式与参与系数
第n阶模态的有效质量Mn,eff = (Γn)²/Mn（Γn=模态参与系数，Mn=广义质量）。理论上保证所有模态的有效质量之和等于总质量，利用这一性质可以估算“计算中遗漏的模态”的质量。FEM软件通常在特征值分析后自动输出这些数据。
线性单元（一阶单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切自锁（可通过减缩积分或B-bar法缓解）。
二阶单元（带中间节点）
可以表现曲线变形。应力精度大幅提高，但自由度约增加2~3倍。推荐：应力评估很重要的情况。
完全积分 v
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告
なった
詳しく
報告