刚体单元
刚体单元的理论基础
刚体单元是什么
老师，FEM中的"刚体单元"是什么？不变形的单元？
是的。刚体单元不变形。将两个或以上节点「刚体结合」，使一方的位移·旋转传递到另一方。
物理上表现什么样的结构？
数学定式化
刚体单元在FEM中作为多点约束（MPC: Multi-Point Constraint）实现。对于独立节点（主节点）的位移 $\{u_m\}$，从属节点（从节点）的位移为：
其中 $[T]$ 是位移变换矩阵，基于刚体运动的公式。
从属节点是从独立节点的位移"自动确定"的吧。
没错。从属节点的自由度被消去（缩约）。实际上DOF会减少。
刚体单元的种类
| 种类 | Nastran | Abaqus | 特征 |
|---|---|---|---|
| 刚体结合（全DOF） | RBE2 | RIGID BODY / COUPLING | 全自由度刚体结合 |
| 加权平均（荷载分配） | RBE3 | *DISTRIBUTING COUPLING | 分配力。不添加刚性 |
| 刚体面 | RBAR | *RIGID BODY | 2个节点的刚体链接 |
RBE2和RBE3完全不同吗？
本质上不同。RBE2是"刚体结合"，向结构增加刚性。RBE3是"荷载分配"，不增加刚性。如果理解不到位，就会导致严重的模型化错误。详细说明见RBE2、RBE3页面。
刚体单元的注意事项
刚体单元很方便，但用法不当会严重扭曲结果：
1. 刚体单元导致局部刚性过大 — 连接处变得异常坚硬
2. 产生应力奇点 — 刚体单元与变形单元的边界处应力不精确
3. 过约束 — 无意中约束了不必要的自由度
刚体单元边界处的应力不可信吗？
与刚体单元相连接的变形单元的最初1~2层应力仅供参考。应评估远处的应力。
总结
整理一下刚体单元的理论。
要点：
- 不变形单元 — 节点间刚体结合
- MPC方式实现 — 从属节点的DOF被消去
- RBE2（刚体结合）和RBE3（荷载分配）完全不同 — 严禁混淆
- 注意局部刚性过大 — 连接处应力不精确
- 方便用于不同单元间的连接 — 梁-实体、壳-实体
RBE2 vs. RBE3的理解是FEM建模最重要的技能之一呀。
正是如此。能否正确区分这两者，能看出一个FEM工程师的能力。
刚体单元的力学基础
刚体单元是没有内部变形的理想刚体，作为有限单元出现于20世纪60年代NASA的太空结构分析中，由特纳（M.J. Turner）等人提出。从数学上看，通过刚体变换矩阵T，从独立节点的6个自由度可以确定性地计算从属节点的全部6个自由度。现实材料具有弹性变形，因此刚体单元可以理解为模拟真实刚性的弹性率至少是实际材料的10~100倍。
刚体单元的数值计算方法
MPC（多点约束）的实现
刚体单元的MPC怎样实现？
有两种方法：
1. 从属DOF的消去（变换法）
将从属节点的DOF用独立节点的DOF表示，从全局方程中消去。联立方程的规模减小。Nastran的RBE2采用此方法。
2. 惩罚法
将刚体约束近似为"非常硬的弹簧"。不是完全刚体，但精度上足以使用。惩罚值过大会恶化条件数。Abaqus的KINEMATIC选项是变换法，PENALTY选项是惩罚法。
哪种更好？
变换法（KINEMATIC）是基本推荐。惩罚法在与接触组合的情况下使用（为了与接触的惩罚法保持一致）。
求解器别的实现
| 功能 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 刚体结合 | RBE2 | RIGID BODY / COUPLING, KINEMATIC | CERIG |
| 荷载分配 | RBE3 | *COUPLING, DISTRIBUTING | RBE3(MPC) |
| 刚体链接（2节点） | RBAR | *RIGID BODY | MPC184 |
| 刚体面 | RBE2（多节点） | *RIGID BODY | TARGE170 |
刚体单元的连接模式
典型的连接模式：
梁单元和实体单元的连接
梁单元有旋转DOF，实体单元没有。用RBE2/RBE3连接：
- RBE2 — 梁端点作为主节点，实体面的节点作为从节点。截面刚体运动（容易太硬）
- RBE3 — 梁端点作为从节点（参考点），实体面的节点作为主节点（荷载分配点）。不添加刚性（更现实）
RBE3的连接"更温和"吧。
完全同意。RBE2是"硬焊缝"式的连接，RBE3是"平衡荷载传递"式的连接。实际的螺栓连接往往更接近RBE3。
总结
整理刚体单元的数值方法。
要点：
- 变换法（KINEMATIC）是基本推荐 — DOF消去，精确
- 惩罚法用于与接触的结合 — 近似的刚体
- 梁-实体的连接 — RBE2（硬）vs. RBE3（软）
- 连接处应力不精确 — 1~2单元远处评估
刚体单元的实现方式比较
刚体单元的实现主要有①惩罚法、②拉格朗日乘数法、③直接刚性缩约三种。惩罚法实现简单但罚则系数的选择影响解精度。拉格朗日乘数法数学上严格但增加联立方程规模。Nastran采用直接缩约法，消去从属DOF的同时保持系统方程规模。
刚体单元的实务应用
刚体单元的实务应用
请告诉我刚体单元的实务用法。
整理最常见的应用。
螺栓接合的建模
简化表现螺栓接合部分时：
RBE2方式 — 螺栓孔周围的节点用RBE2结合。"螺栓孔部分是刚体"处理。简单但接合处太硬。
RBE3+弹簧方式 — 螺栓位置用RBE3分配荷载，用弹簧单元表现螺栓的轴刚度。更现实但设置复杂。
荷载的分配
"在一点施加荷载，分配到面上"的场景中RBE3很有用。例如吊机荷载分配到法兰、千斤顶荷载分配到结构等。
边界条件的表现
"面的平均位移为零""面保持平面"等边界条件用RBE2/RBE3表现。压力容器的端面条件（面保持平面+轴向自由移动）用RBE2设置。
实务检查清单
请提供刚体单元的检查清单。
"RBE2 vs. RBE3的判断"是首要检查项。这决定了全部。
应该在RBE2的地方用RBE3会导致刚性不足，在RBE3的地方用RBE2会导致刚性过大。这个判断错误会根本改变结果。
刚体单元在发动机悬置分析中的应用
汽车发动机悬置分析中，将发动机缸体用一体刚体单元（RBE2或RBAR集合）模型化，大幅降低计算成本。本田公开技术文献（2015年）报告称，从详细FEM替换为刚体置换可以将模型规模缩减约60%，同时保持悬置部反力精度在5%以内。
刚体单元的软件比较
刚体单元的工具比较
请比较各求解器的刚体单元。
Nastran的RBE2/RBE3最有实绩吗？
可以说RBE2/RBE3是由Nastran发明的。20世纪70年代以来的历史，航空航天的大规模模型有膨大的使用实绩。Abaqus和Ansys的同等功能本质上相同，只是名称和设置方法不同。
选型指南
"仅仅是名称不同"，但设置方法差异很大吧。
Nastran的RBE2/RBE3卡片记法和Abaqus的*COUPLING记法完全不同。很重要的是熟悉各求解器的手册，能用不同的记法设置相同的物理现象。
刚体单元各求解器实现概览
刚体单元的主要实现包括Nastran的RBAR/RBE1/RBE2/RBE3/RROD/RTRPLT、Abaqus的Rigid Body/Kinematic Coupling、Ansys的CE/CP/RBE3指令组。Siemens NX Nastran 2022增加了AUTORBE自动刚体单元生成功能，能从CAD实体直接提取刚体子结构。LS-DYNA多采用基于接触的约束代替刚体单元。
刚体单元的先端研究
刚体单元的先端话题
刚体单元有先端研究吗？
刚体单元本身很简单，但在子结构法和多尺度分析中起着关键作用。
超单元与RBE2/RBE3
超单元（子结构法）中，将结构的一部分缩约成"一个大的单元"。缩约部分的接口用RBE2/RBE3连接。航空航天的大规模模型（数千万DOF）中不可或缺的技术。
数字样机
汽车的数字样机（全车模型）中使用数千个RBE2/RBE3。车身、底盘、发动机、悬置的连接全部用刚体单元表现。连接的建模支配全车振动特性。
与接触的整合
Abaqus的*RIGID BODY可用作接触的刚体面（模具、冲头、工具）。冲压成形中，将模具设置为刚体面，毛坯为变形体建模。
总结
整理刚体单元的先端话题。
刚体单元是看不见的幕后英雄，但是大规模分析的基础设施不可缺少。
刚体单元与柔性体的耦合模型
多体动力学（MBD）与FEM的耦合中，刚体单元作为柔性体（CMS：部件模式综合法）的惯性特性参考节点起作用。20世纪90年代SimPack与ADAMS/Flex的耦合分析普及，柔性机械手臂末端振动分析中基于刚体单元的Craig-Bampton缩约成为业界标准。
刚体单元的故障处理
刚体单元的故障
请告诉我刚体单元常见的故障。
刚体单元的故障主要原因是无意中的刚性增加和DOF过约束。
RBE2局部刚性过大
用RBE2后，接合处的应力异常高了。
RBE2完全刚体化了连接部。原本应该变形的部分变成刚体，导致应力集中在相邻变形单元。
对策：
- 切换到RBE3 — 仅分配荷载，不增加刚性
- 不评估接合处的应力 — 1~2单元远处评估
- 用弹簧单元替换 — 有限刚性连接
RBE3结构不稳定
用RBE3后出现"奇异刚性矩阵"错误。
RBE3不增加刚性。只用RBE3支持结构的话，结构会不稳定（刚体移动没约束）。
对策：
- 在RBE3的参考点（从节点）另外添加约束
- RBE3仅用于"荷载分配"，"结构支持"用RBE2
DOF过约束
刚体单元太多，导致过约束。
用RBE2约束全部6个自由度会完全固定连接部的变形。如果某些自由度（例如旋转）应该自由，就约束必需的DOF。
Nastran的RBE2中用CM（约束成分）字段指定123456中的一部分。全DOF约束就变成"完全固定"。
总结
整理刚体单元的故障处理。
记住这个根本原则，大部分故障就能防止。
完全同意。理解RBE2和RBE3的本质区别是FEM建模最重要的技能。
刚体单元过约束的检查方法
多个刚体单元在同一节点约束不同DOF导致的"过约束"是分析的致命错误。MSC Nastran中含有MPC/RBE单元时常出现Fatal 9050。诊断步骤是①检查AUTOSPC输出的SPC力，②用PARAM,CHECKOUT,YES输出MPC约束矩阵，③搜索重复DOF号，共3步。Abaqus中会表现为局部位移突变。
价值
更详细
错误