惯性释放

定义

惯性释放处理的对象是约束不足、允许发生刚体运动的结构。这类结构的整体刚度矩阵是奇异的，常规静力求解无法直接算出位移^[2]。惯性释放并不去掉刚体运动，而是在外载荷之外补充一组与之相平衡的惯性载荷，使结构在等效静力意义上处于平衡，从而解出弹性变形与应力^[1]。

计算结果中的位移，是各节点相对于参考点（虚支座）的相对位移，参考点选在哪里只改变位移数值，不改变结构内部的变形与应力。它既能用于外载荷本身已经平衡的情形，也能用于载荷合力不为零、结构整体带有加速度的情形。

原理

惯性释放的依据是达朗贝尔原理：给加速运动的物体补上惯性力与惯性力矩之后，就可以把它当作等效的静力系统求解。设 $K$ 为无约束结构的刚度矩阵、$u$ 为位移、$F$ 为外载荷的等效节点力、$M$ 为质量矩阵、$a_r$ 为外载荷引起的刚体加速度，平衡关系可写成 $K u = F - M a_r$^[4]。外载荷与惯性载荷叠加后构成自平衡力系，因此虚支座处的反力理论上为零^[5]。

求解通常分四步：先由有限元模型组装整体质量矩阵、刚度矩阵和外载荷等效节点力；再选定参考点，把质量矩阵与载荷向量缩聚到该点，解出参考点在 3 个平动和 3 个转动方向上的刚体加速度；随后按刚体运动条件把加速度分布到全部节点，得到节点惯性载荷；最后把惯性载荷与外载荷叠加，用虚支座约束参考点的 6 个自由度，按常规静力问题求出位移、应变与应力^[6][7]。

mermaid\nflowchart LR\n A[有限元模型：质量矩阵、刚度矩阵与外载荷] --> B[选取参考点并缩聚]\n B --> C[求解参考点刚体加速度]\n C --> D[分布惯性载荷并与外载荷叠加]\n D --> E[虚支座约束参考点6个自由度]\n E --> F[求解位移、应变与应力]\n

求解完成后可用两项指标检验虚支座是否选得合理：一是约束点处反力之和趋于零，二是与单位刚体位移对应的应变能为零，后者说明模型可以无应变地作刚体运动，实际计算中受舍入误差影响会出现很小的非零值^[7][8]。惯性载荷算得准不准，取决于质量矩阵是否准确，所以模型中必须定义材料的质量密度^[6]。

发展历程

惯性释放的方法在 20 世纪 70 年代已见于工程文献。1977 年在美国密歇根州南菲尔德举行的第二届车辆结构力学国际会议上，有研究者报告了用惯性释放估算结构冲击载荷的做法：先取力的时间历程中的代表值作为静载施加于自由-自由结构，算出刚体加速度，再由加速度与质量得到各点惯性力，最后约束结构并作常规静力分析^[9]。该研究指出，当载荷周期远大于结构被激振模态的周期时，这种做法给出的结果是精确的，两者接近时精度则明显下降^[9]。

此后这一功能进入航空航天用有限元软件。MSC/NASTRAN 提供专门的惯性释放求解序列，用惯性载荷平衡施加在自由-自由结构上的稳态载荷^[10]；其中 Solution 91 即为含惯性释放的超单元静力求解序列，全模型由同一个参考点统一处理，并允许直接输入刚体加速度^[11]。20 世纪 90 年代，美国国家航空航天局刘易斯研究中心与 MacNeal-Schwendler 公司的技术报告还给出了用程序改动把随位移变化的载荷并入惯性释放静力分析的闭式解法^[12]。

随后各主流软件陆续集成该功能。MSC Nastran 提供手动与自动两种方式：手动方式需指定虚支座节点并设置相应参数，自动方式由程序自行选取参考系并把约束分配到有质量的节点，后者更适合存在建模误差、结构缺陷或由不连通部件组成的模型，但不能用于迭代求解器^[7]。ANSYS 通过 IRLF 命令启用惯性释放，只能作线性静力分析，且不支持非线性、轴对称与广义平面应变单元；Abaqus 也提供惯性释放载荷，并可结合几何非线性处理大变形问题^[7]。OptiStruct 自 8.0 版起可直接通过控制卡片激活该功能。

近年的研究转向求解效率与适用范围。德国航空航天中心的研究人员 2024 年在美国航空航天学会科技论坛上提出了一种面向大变形的非线性惯性释放列式，用于把结构响应与飞行力学方程耦合求解^[13]。《计算力学学报》刊出的研究把惯性释放有限元系统归结为增广矩阵法、基于广义逆的直接解法和 Krylov 子空间迭代法三种求解思路，结果显示前两者的计算时间和内存占用呈超线性增长，只适合中小规模问题，而迭代法接近线性增长，在大规模工程问题中更有优势^[2]。在结构优化方面，也出现了把惯性释放与自由漂浮弹性体拓扑优化相结合的方法研究^[4]。

应用

航空航天是惯性释放应用最早、也最普遍的方向。运载器与航天器在准静态飞行段被当作自由-自由系统，其位移常借助带惯性释放的静力分析求得^[10]。典型对象包括承受气动载荷的火箭整流罩、运载火箭的起飞与着陆工况^[1][8]，以及航空发动机中介机匣、风扇轴等主承力构件的强度分析^[3]。

在汽车行业，惯性释放用于车架等部件的结构优化与强度计算，也用于白车身关键接附点的刚度分析：由于约束区域的大小和位置不同会带来较大差异，引入惯性释放可以减小这类影响^[3]。整车行驶工况中，轮胎接地力与惯性力共同实现平衡^[1]。有研究以某高过载干扰弹头为例，比较了惯性释放与限制约束两种做法，二者得到的最大应力与应力集中位置差别明显，且试验结果支持惯性释放的仿真结论^[14]。

船舶结构分析中，船体静止于水中时重力与浮力相抵，整体处于无约束状态，正适合用惯性释放模拟。全船结构强度直接计算既可采用支座法，也可采用惯性释放虚支座法，后者只需约束某一个节点的 6 个方向自由度。已有应用案例包括某观光潜艇与江海直达推船的全船结构强度直接计算^[15]，以及某型舰载升降机升降平台的强度与刚度有限元计算、某型无人直升机的整机刚度仿真等^[3]。

局限

惯性释放只适用于线性静力分析。以 ANSYS 为例，其惯性释放不允许非线性分析，也不支持轴对称单元与广义平面应变单元^[7]。同时，模型必须具有真实的质量分布并定义材料密度，否则算出的惯性载荷不可用^[7]。

惯性载荷的分布建立在刚体加速度为常量的假设之上，平动加速度由外载荷合力决定，转动加速度则与求和参考点的选取有关^[7]。部分软件只计入平动加速度场，忽略由此产生的转动分量。还有软件明确不建议用惯性释放去平衡重力、离心力或热载荷，这样处理可能得到错误结果。

用惯性释放处理动力问题时只是一种近似。只有当载荷周期远大于结构被激振模态的周期时结果才精确，载荷周期与结构周期接近时结果明显变差；在串联的质量系统中，离载荷施加位置越远，误差越大^[9]。此外，虚支座的自由度必须与模型实际的自由度匹配，选得不合适会影响惯性载荷的计算精度；自动惯性释放虽然免去了人工选取约束的麻烦，但不适用于迭代求解器^[7]。由于结构本身没有约束，得到的位移还需换算才能对应真实位移。

参见

• 达朗贝尔原理 —— 惯性释放所依据的力学原理，用惯性力把动力学问题转化为等效静力问题

• 有限元法 —— 惯性释放通常在有限元框架内实现，是处理无约束结构的手段之一

• 刚度矩阵 —— 无约束结构的该矩阵奇异，是必须采用惯性释放的直接原因

• MSC Nastran —— 较早提供惯性释放功能的有限元软件，设有专门的求解序列

• 模态分析 —— 惯性释放处理动力载荷时的适用条件与结构固有频率有关