集中質量要素
理论与物理
什么是集中质量单元
老师，“集中质量单元”是什么？静力分析中会需要用到质量吗？
集中质量单元是在单个节点上附加质量的单元。在静力分析中用于计算自重（重力载荷），但其主要用途是动态分析（固有振动、时程响应、冲击）。
质量矩阵
集中质量单元的质量矩阵是对角矩阵：
其中 $m$ 是平移质量，$I_{xx}, I_{yy}, I_{zz}$ 是转动惯量。
原来转动惯量也可以设置啊。
当用质量单元替代电机或大型设备等旋转部件时，如果没有转动惯量，固有振型就会不准确。不仅平移质量，正确设置转动惯量也很重要。
用途
集中质量单元的主要用途：
| 用途 | 说明 |
|---|---|
| 附加非结构质量 | 管道内的流体、设备重量 |
| 质量平衡调整 | 将FEM模型的整体质量调整至设计值 |
| 设备建模 | 用质点替代电机、阀门等 |
| 重心调整 | 使重心位置与实际结构一致 |
| 质量缩放（显式解法） | 增大稳定时间增量 |
管道内流体的质量怎么添加呢？
在管道单元的每个节点上，以集中质量的形式添加流体单位长度的质量。在Abaqus中，NONSTRUCTURAL MASS 很方便，可以向单元添加面密度（kg/m²）或线密度（kg/m）。
各求解器的单元名称
| 求解器 | 单元名称 | 备注 |
|---|---|---|
| Nastran | CONM2 | 6自由度集中质量。可偏移 |
| Abaqus | MASS / ROTARY INERTIA | 平移和转动分开定义 |
| Ansys | MASS21 | 通过KEYOPT(3)选择有无转动惯量 |
Nastran的CONM2“可偏移”？
CONM2可以将质量的质心从节点位置偏移。这样即使质心不在节点位置，也能表现正确的惯性效应。对于大型设备的安装点与质心不重合的情况很有用。
总结
我来整理一下集中质量单元。
要点：
- 在单个节点上附加质量的单元 — 不具有刚度
- 主要用于动态分析 — 固有振动、时程响应
- 平移质量与转动惯量 — 两者都应设置（别忘了转动惯量）
- 附加非结构质量 — 流体、设备、附属物的质量
- CONM2（Nastran）的偏移 — 质心位置的精确建模
在静力分析中只与自重有关，但在动态分析中支配着固有频率，所以超级重要啊。
没错。如果结构的质量分布错了，固有频率也会错。质量单元的设置是动态分析最基础的部分。
集中质量的理论背景
集中质量矩阵（Lumped Mass Matrix）相对于一致质量矩阵（Consistent Mass Matrix），是一种将单元质量集中分配到节点的近似方法。1968年，Hinton等人证明了集中质量矩阵数值衰减更少，更适合显式时间积分。Abaqus Explicit在所有单元中都采用了这种方式，支撑着现代爆炸、冲击分析的计算效率。
各项的物理意义
- 惯性项（质量项）：$\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中有“集中质量”和“分布质量（一致质量）”，有什么区别呢？
通常的FEM单元（梁、壳、实体）会根据材料密度 $\rho$ 自动生成质量矩阵。这就是一致质量矩阵（consistent mass matrix）。而集中质量单元是向特定节点赋予附加质量。
此外，还有将一致质量矩阵对角化后的集中质量矩阵（lumped mass matrix）。
| 质量矩阵 | 特点 | 用途 |
|---|---|---|
| 一致质量 | 有非对角成分。精度高 | 隐式解法的动态分析 |
| 对角集中质量 | 仅对角成分。计算效率高 | 显式解法（Explicit） |
| 附加集中质量（CONM2等） | 用户添加的质量 | 附加非结构质量 |
显式解法中必须使用对角集中质量吗？
是的。显式解法的时间积分（中心差分法）需要使用质量矩阵的逆矩阵，如果不是对角矩阵就无法高效计算。LS-DYNA和Abaqus/Explicit会自动将质量矩阵对角化。
质量确认方法
如何确认模型的整体质量是否正确呢？
确认求解器的质量摘要输出：
确认项目：
- 整体质量是否与设计值（图纸值、实测值）一致
- 重心位置（$X_{CG}, Y_{CG}, Z_{CG}$）是否合理
- 各方向的转动惯量是否合理
重心位置的确认很容易被忽略呢。
即使整体质量一致，如果重心偏移，动态响应（特别是倾覆力矩、偏心效应）也会出错。质量摘要必须在分析开始时确认。
质量缩放
“质量缩放”是什么？
显式解法的稳定时间增量由 $\Delta t \propto L / c$（$L$: 最小单元尺寸，$c$: 声速）决定。如果存在很小的单元，$\Delta t$ 会变得极小。人为增加质量以降低 $c$，从而增大 $\Delta t$ 就是质量缩放。
增加质量不会改变惯性效应吗？
会改变。因此质量缩放只应适用于准静态问题，并且需要确认增加的质量（动能）占总能量的5%以下。
总结
我来整理一下集中质量的数值方法。
なった
詳しく
報告