在汽车碰撞仿真中，为什么吸能盒的设计要让塑性变形能尽量大？

根据能量守恒原理，碰撞的初始动能会转化为弹性变形能、塑性变形能和热量。设计良好的吸能盒旨在最大化塑性变形能，因为这部分能量被结构永久吸收并耗散，从而显著减少传递到乘员舱的剩余动能。这样能更有效地保护乘员安全，是汽车被动安全设计的核心目标之一。

有限元分析中的沙漏能是什么？为什么需要控制它？

沙漏能是使用一阶减缩积分单元（如LS-DYNA中某些壳单元）时可能出现的一种虚假能量。它对应于一种不产生实际应力的变形模式（如沙漏状畸变）。为了抑制这种非物理的数值现象，需引入人为的沙漏控制刚度，其做功即为沙漏能。通常要求沙漏能不超过内能的5%，以确保仿真结果的可靠性和准确性。

做完显式动力学仿真后，为什么首先要检查能量曲线图？

检查能量曲线是显式动力学分析后的关键步骤。通过观察总能量是否守恒（无跳变）、沙漏能占比是否过高（通常<5%）、以及动能与内能的转化是否符合物理过程（如碰撞中动能应转化为内能），可以快速判断模型的整体稳定性与可靠性。能量曲线的异常往往比应力云图更早地揭示出模型设置或接触定义中的根本性错误。

在有限元仿真中，如何判断沙漏能是否在可接受范围内？

判断标准主要看沙漏能与内能（即结构的真实变形能）的比值。行业普遍经验是，沙漏能应控制在总内能的5%以内。如果占比过高，说明虚假的沙漏模态过于显著，可能严重扭曲应力、应变结果，导致仿真不可信。此时需要调整沙漏控制参数、网格或单元类型来抑制沙漏现象。

功与能量 — 弹性应变能与有限元虚功原理

Q: 在有限元仿真中，如何判断沙漏能是否在可接受范围内？

判断标准主要看沙漏能与内能（即结构的真实变形能）的比值。行业普遍经验是，沙漏能应控制在总内能的5%以内。如果占比过高，说明虚假的沙漏模态过于显著，可能严重扭曲应力、应变结果，导致仿真不可信。此时需要调整沙漏控制参数、网格或单元类型来抑制沙漏现象。

分类：物理基础 | 2026-03-25 | サイトマップ

NovaSolver Contributors

能量方法在CAE中的地位
功的定义与计算
动能与势能
弹性应变能
能量守恒在结构分析中的应用
虚功原理：FEM的弱形式基础
从虚功原理推导FEM方程
工程实例：汽车碰撞能量吸收
能量平衡作为FEM质量检验

1. 能量方法在CAE中的地位

🧑‍🎓

碰撞仿真里经常说"吸收能量"——汽车碰撞的动能到底去哪里了？是变成了什么？

🎓

主要变成了材料的塑性变形能（永久变形存储在内部）和热量（金属塑性变形约有5%到10%转化为热）。从能量守恒角度：初始动能 = 弹性变形能 + 塑性变形能（耗散）+ 热量。设计良好的吸能盒，就是让前车头的塑性变形能尽量大，使传到乘员舱的动能尽量少。

在结构力学中，能量方法提供了一种与直接用方程求解力和位移截然不同的视角。虚功原理（Principle of Virtual Work）是大多数FEM软件底层数学的基础，理解它就理解了FEM为什么"工作"。

2. 功的定义与计算

力 $\vec{F}$ 在位移 $d\vec{u}$ 上做的功：

$$dW = \vec{F} \cdot d\vec{u} = F \, du \cos\theta$$

其中 $\theta$ 是力与位移方向的夹角。沿路径的总功：

$$W = \int_{\vec{u}_1}^{\vec{u}_2} \vec{F} \cdot d\vec{u}$$

几种特殊情况：

力垂直于位移（$\theta = 90°$）：做功为零（如向心力对圆周运动）
力与位移同向（$\theta = 0°$）：做正功
力与位移反向（$\theta = 180°$）：做负功（阻力、弹性恢复力）

在FEM中，外力做的功（外功）等于内力（应力）做的功（内功），这就是虚功原理的本质。

3. 动能与势能

动能（Kinetic Energy）：

$$T = \frac{1}{2}mv^2 = \frac{1}{2}\int_V \rho \dot{u}_i \dot{u}_i \, dV = \frac{1}{2}\{\dot{u}\}^T [M] \{\dot{u}\}$$

弹性势能（Elastic Potential Energy / Strain Energy）：

$$\Pi = \frac{1}{2}\int_V \sigma_{ij}\varepsilon_{ij} \, dV = \frac{1}{2}\{u\}^T [K] \{u\}$$

外力势能：

$$V_{\text{ext}} = -\{u\}^T \{F\}$$

系统总势能（Potential Energy）：

$$\Pi_{\text{total}} = \frac{1}{2}\{u\}^T [K] \{u\} - \{u\}^T \{F\}$$

FEM的静力分析其实是在求 $\Pi_{\text{total}}$ 对位移的极小值：

$$\frac{\partial \Pi_{\text{total}}}{\partial \{u\}} = [K]\{u\} - \{F\} = 0 \implies [K]\{u\} = \{F\}$$

这就是最小势能原理，它与牛顿定律等价，但从能量角度给出了FEM的另一种推导路径。

4. 弹性应变能

弹性应变能密度（单位体积的应变能）：

$$u = \frac{1}{2}\sigma_{ij}\varepsilon_{ij} = \frac{1}{2}\{\sigma\}^T \{\varepsilon\}$$

对于单轴线弹性情况（胡克定律 $\sigma = E\varepsilon$）：

$$u = \frac{1}{2}E\varepsilon^2 = \frac{\sigma^2}{2E}$$

弹性应变能在工程中的应用：

应用场景	应变能的意义	工程指标
弹簧设计	储能容量	应变能密度 J/m³
弹道防护	装甲吸收弹丸动能	比吸收能（SEA）J/kg
汽车吸能盒	塑性变形吸能	平均碰撞力 × 压溃行程
复合材料损伤	分层起始判据（GIC）	临界应变能释放率 J/m²
断裂力学	裂纹扩展驱动力	J积分 J/m²

应变能在疲劳分析中也有直接应用——Coffin-Manson关系中塑性应变能密度与疲劳寿命直接相关：

$$\Delta W_p = \sigma_f' \varepsilon_f' (2N_f)^{b+c}$$

5. 能量守恒在结构分析中的应用

封闭系统的能量守恒：

$$E_{\text{kinetic}} + E_{\text{internal}} + E_{\text{contact}} + E_{\text{hourglass}} = E_{\text{external}} + E_{\text{initial}}$$

在显式动力学仿真中，各能量项的物理含义：

$E_{\text{kinetic}}$：系统动能（$\frac{1}{2}\dot{u}^T M \dot{u}$）
$E_{\text{internal}}$：内能（弹性应变能 + 塑性耗散能）
$E_{\text{contact}}$：接触能（接触界面的摩擦耗散）
$E_{\text{hourglass}}$：沙漏控制能（人为引入，应保持极小）
$E_{\text{external}}$：外力做的功

🧑‍🎓

沙漏能是什么？从来没听说过这个词……

🎓

这是FEM里特有的数值现象。一阶减缩积分单元（如LS-DYNA默认的shell单元）计算效率很高，但存在一种不产生任何应力的变形模式——就像沙漏的形状那样，节点位移有了但积分点处应变为零。为了抑制这种假模态，需要加入人为的沙漏控制刚度，做的功就是沙漏能。通常要求沙漏能不超过内能的5%，否则结果不可信。

6. 虚功原理：FEM的弱形式基础

虚功原理（Principle of Virtual Work）：对于处于平衡状态的弹性体，对任意满足位移边界条件的虚位移 $\delta u$，内力虚功等于外力虚功：

$$\delta W_{\text{int}} = \delta W_{\text{ext}}$$

展开：

$$\int_V \delta\varepsilon_{ij} \sigma_{ij} \, dV = \int_V \delta u_i f_i \, dV + \int_S \delta u_i t_i \, dS$$

其中：

左边：虚应变 × 应力，对整个体积积分 = 内力虚功
右边第一项：虚位移 × 体积力 = 重力虚功
右边第二项：虚位移 × 面力（边界上的力） = 外载荷虚功

虚功原理与强形式（直接写平衡微分方程）等价，但它：

对位移的可微性要求更低（弱形式），允许C⁰连续性
自然地处理了力的边界条件（自然边界条件）
天然地适合有限元离散化处理

7. 从虚功原理推导FEM方程

将位移用形函数近似 $u \approx [N]\{u_e\}$，代入虚功原理：

$$\delta u_e^T \int_V [B]^T \{\sigma\} \, dV = \delta u_e^T \left(\int_V [N]^T \{f\} \, dV + \int_S [N]^T \{t\} \, dS\right)$$

由于 $\delta u_e$ 任意，两边可消去，得到单元方程：

$$\underbrace{\int_V [B]^T [D] [B] \, dV}_{[K]^{(e)}} \{u_e\} = \underbrace{\int_V [N]^T \{f\} \, dV + \int_S [N]^T \{t\} \, dS}_{\{F\}^{(e)}}$$

这正是 $[K]^{(e)}\{u_e\} = \{F\}^{(e)}$！将所有单元组装后得到全局方程 $[K]\{u\} = \{F\}$。

虚功原理推导的关键步骤：

选定位移场近似（形函数 $[N]$）
计算应变-位移矩阵 $[B] = \partial[N]$
选定本构关系（材料矩阵 $[D]$）
对单元体积积分（高斯积分）
组装为全局矩阵

8. 工程实例：汽车碰撞能量吸收

以NCAP正面碰撞测试（56 km/h = 15.56 m/s）为例，评估能量吸收设计：

车辆初始动能：

$$E_k = \frac{1}{2}mv^2 = \frac{1}{2} \times 1500 \times 15.56^2 \approx 181 \text{ kJ}$$

设计目标：前车头吸能区（吸能盒+前纵梁）吸收约70%动能，即127 kJ；乘员舱最大侵入量 < 100 mm。

能量吸收部件	吸能量（估算）	吸能机制
吸能盒（前端）	15-25 kJ	轴向屈曲折叠（Progressive buckling）
前纵梁	40-60 kJ	弯曲 + 折叠变形
发动机舱横梁	20-35 kJ	弯曲变形
门槛梁	10-20 kJ	弯曲强化
其他结构	5-10 kJ	残余变形

衡量吸能效率的指标——比吸能（SEA, Specific Energy Absorption）：

$$SEA = \frac{E_{\text{absorbed}}}{m_{\text{structure}}} \quad \text{[J/kg]}$$

高SEA材料和结构对汽车轻量化设计至关重要：

材料/结构	SEA (kJ/kg)	优点
铝合金挤出型材	30-60	轻量化，可设计性好
高强钢（AHSS）	15-35	成本低，成形性好
CFRP（碳纤维）	60-120	极轻，高端应用
铝蜂窝夹芯	40-80	可控变形，空间利用率高
泡沫铝	10-25	各向同性，填充复杂形状

9. 能量平衡作为FEM质量检验

在LS-DYNA中，能量平衡检查是验证仿真模型正确性的重要手段。完整的能量方程：

$$E_k(t) + E_{\text{int}}(t) = E_{\text{ext}}(t) + E_{\text{init}} - E_{\text{contact}}(t) - E_{\text{viscous}}(t)$$

常见的能量异常和对应的诊断：

能量异常现象	可能原因	排查方法
能量突增（跳变）	单元严重扭曲（负体积）	检查最小单元体积，启用侵蚀
沙漏能 > 内能 × 10%	减缩积分单元沙漏失控	加密网格或使用全积分单元
总能量不守恒	接触穿透，质量缩放过大	检查接触设置，减小质量缩放比
能量一直增大（静力分析）	未施加足够约束（刚体运动）	检查边界条件，确认约束充分

🧑‍🎓

所以看能量曲线是判断仿真结果可靠性的一个快速手段？

🎓

对，这是显式动力学工程师的基本功。每次跑完仿真，第一件事就是看energies plot：总能量曲线有没有跳变？沙漏能占比多少？动能和内能的演化是否符合物理直觉？能量曲线的形状本身就会告诉你模型有没有问题，比看应力云图更早发现错误。

总结

能量方法是FEM的灵魂。从功-能定理到虚功原理，从弹性应变能到塑性耗散——这些物理概念直接决定了FEM方程的推导方式，也提供了验证仿真结果正确性的物理工具。汽车碰撞吸能设计的背后，就是对这些能量守恒和转化规律的精确控制。

关联话题： 牛顿运动定律动量与冲量胡克定律与刚度矩阵应力与应变入门结构力学基础