LS-DYNA — CAE术语解说

分类:术语表 | 2026-01-15
CAE visualization for ls dyna - technical simulation diagram

LS-DYNA

🧑‍🎓

老师,LS-DYNA是什么样的求解器?听说在碰撞分析中很有名。

LS-DYNA的理论基础

LS-DYNA的基本概念和控制方程

🧑‍🎓

LS-DYNA中经常提到的"显式法"和"隐式法"的根本区别是什么?不只是计算快慢的问题吧?

🎓

说得对,本质区别在于时间积分方法。显式法仅用当前时刻

$$ t_n $$
的信息直接计算下一时刻
$$ t_{n+1} $$
的位移和速度。例如中心差分法,从运动方程
$$ M \ddot{u}_n = F_n^{\text{ext}} - F_n^{\text{int}} $$
出发,通过对质量矩阵
$$ M $$
进行对角化,无需求解联立方程组就能更新。这正是它高速的原因。

🧑‍🎓

那隐式法反过来需要求解联立方程组所以慢?但为什么对碰撞这样的短时间现象,显式法更合适呢?

🎓

是的,隐式法在

$$ t_{n+1} $$
包含的方程中求解,需要迭代计算(如Newton-Raphson法),所以成本很高。显式法适合碰撞和爆炸的原因是,时间步长必须非常小(例如1微秒)。用隐式法来求解这么小的步长,迭代次数会爆炸式增长。而显式法每步计算轻,虽然步数多,但总耗时仍然较少。但要注意的是,显式法的稳定性受CFL条件制约。时间步长
$$ \Delta t \propto L_e / c $$
(其中
$$ L_e $$
是单元特征长度,
$$ c $$
是音速),例如1mm单元在钢材中约为0.2微秒的上限。

🧑‍🎓

支配方程应该就是"运动方程"吧?但我听说它也能解流体和热传导,这些怎么处理呢?

🎓

LS-DYNA的核心确实是拉格朗日形式的运动方程,但为了支持多物理场,它引入了"任意拉格朗日-欧拉(ALE)"法和"光滑粒子流体动力学(SPH)"。例如,在流固耦合(FSI)分析中,对流体应用ALE法,对结构体应用拉格朗日法,然后在界面处交换数据。流体的支配方程包括压缩性欧拉方程,如

$$ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 $$
等。

🧑‍🎓

"拉格朗日"和"欧拉"的视点区别能否用LS-DYNA的网格具体解释?

🎓

好的。拉格朗日网格随材料变形而变形。汽车碰撞时,车身被压扁,网格也跟着被压扁。欧拉网格则固定在空间中,材料在其中流动。ALE则是混合体,网格可以自由移动到一定程度,这样可以缓解大变形时网格扭曲的问题。在LS-DYNA中,通过*SECTION_SOLID的ELFORM参数来指定,比如1表示拉格朗日,11表示ALE。

LS-DYNA的数值计算方法

离散化、单元、求解器设置

🧑‍🎓

关键词文件中经常看到的"沙漏控制(Hourglass control)"是什么意思?为什么需要它?

🎓

沙漏是单点积分低阶单元(计算快)特有的数值不稳定模式。单元虽然不变形,却能以零能量随意扭曲,这就叫沙漏模式。沙漏控制就是要抑制这种现象。在LS-DYNA中,通过*HOURGLASS卡片指定类型。例如IHQ=4(Flanagan-Belytschko粘性型)适合冲击问题,IHQ=5(Flanagan-Belytschko刚性型)适合弹性问题。设置不当的话,控制力可能过大,导致结果变硬。

🧑‍🎓

"接触"算法中的"节点到面接触"和"面到面接触"怎样使用才合适?

🎓

"节点到面接触"(NTS)检查一个面(主面)的节点是否与另一个面(从面)的单元段相接触。计算成本低,但非对称。用在主面网格细密的部件上。"面到面接触"(STS)检查两个面的单元段之间是否相接触。对称性好,接触力传递更准确,但成本约为NTS的2倍。在板金成形分析等两表面都需要精确接触的场景中,推荐用STS。卡片中如*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE这样指定。

🧑‍🎓

材料模型很多,比如处理钢铁高应变变形的"MAT_PIECEWISE_LINEAR_PLASTICITY"内部是怎样的?

🎓

MAT_24(通常这么叫)是基于等向硬化的弹塑性模型。输入包括密度、杨氏模量、泊松比、屈服应力,以及应力-塑性应变曲线(LCID)。硬化规律由输入的曲线进行分段线性插值。从公式来说,屈服条件用von Mises准则表示为

$$ \sqrt{\frac{3}{2} s_{ij} s_{ij}} = \sigma_y(\bar{\varepsilon}^p) $$
,其中
$$ s_{ij} $$
是偏应力张量,
$$ \sigma_y $$
是硬化曲线给出的屈服应力。也可以考虑应变率依赖性,通过动态屈服应力的缩放系数来实现。

🧑‍🎓

"质量缩放"具体用什么卡片,改哪个参数?为什么说这是危险的手段?

🎓

通过*CONTROL_TIMESTEP卡片的DT2MS参数,人为增加具有最小时间步的单元的质量,从而增大步长。比如DT2MS=-0.9能让步长扩大10倍。危险之处有两点。第一,质量增加会高估惯性力,动态响应变得不精确。第二,增加的质量会影响结构刚性,固有频率下降。应该仅在静力分析或准静力问题的初期阶段谨慎使用。实际工作中,应该优先改进网格品质。

LS-DYNA的实务应用

工作流程和检查清单

🧑‍🎓

开始碰撞安全分析时,前处理一定要检查的"模型检查"项目有哪些?

🎓

最少要用LS-PrePost或网格工具的检查功能确认以下几点:

1. **自由边和重复节点**:确认没有无意的孔或双重网格。 2. **单元宽纵比和歪斜度**:宽纵比超过10、内角小于10度的极端单元要删除或重新网格划分。 3. **接触面法向量方向**:所有接触面的法向量是否指向正确方向(通常向外)。 4. **材料特性单位一致性**:密度是
$$ ton/mm^3 $$
还是
$$ kg/m^3 $$
,力的单位是否一致。 5. **初始接触**:部件之间初始是否有穿透(负间隙)。
🧑‍🎓

计算中途因"Terminated due to mass increase"停止了。最先要怀疑哪些地方?

🎓

首先怀疑网格大变形产生的极小单元。显式法的时间步长

$$ \Delta t \propto L_e / c $$
$$ L_e $$
是单元特征长度,
$$ c $$
是音速),当
$$ L_e $$
趋近0时,
$$ \Delta t $$
也趋近0。求解器为防止这种情况,会增加该单元的质量,超过一定次数(默认10次)就停止计算。对策是:(1) 粗化该单元所在区域的网格,(2) 强化沙漏控制(改变IHQ),(3) 设置材料破坏准则(*MAT_ADD_EROSION)让单元删除。

🧑‍🎓

后处理时查看能量历史(GLSTAT),什么样的能量平衡才是理想的?

🎓

能量守恒是大前提。即 总能量 = 动能 + 内能 + 滑动界面能 + 沙漏能 应该基本等于初始总能量(主要是初始动能)。特别要注意:

- **沙漏能**:应控制在内能的5%以下。超过10%表示控制过度,结果可能被扭曲。 - **滑动界面能**:接触摩擦和分离消耗的能量。异常过大要检查接触参数(摩擦系数、惩罚刚性)。 - **数值误差能**:越小越好,是数值精度的指标。
🧑‍🎓

实际工作中"网格尺寸"怎么确定?特别是碰撞分析。

🎓

经验法则加权衡。汽车碰撞安全分析中,车身整体平均单元尺寸一般为5-10mm。但重要的变形和折叠区域(前部缓冲盒、支柱)要细化到2-5mm。根据是单元长度应至少为板厚的2倍以上(板厚1mm则单元尺寸2mm以上)。更细的网格会大幅增加计算成本。另外接触的两部件网格尺寸应该接近(比例在3倍以内),有助于接触力的稳定性。

LS-DYNA软件对比

Ansys/LS-DYNA、Abaqus、RADIOSS等的区别

🧑‍🎓

Ansys Mechanical里也有"显式动力学"模块,和LS-DYNA有什么区别?都是Ansys的产品吧?

🎓

历史和架构完全不同。Ansys显式动力学(基于前AUTODYN技术)是"Ansys原生"显式求解器,深度集成到Ansys Workbench环境。而Ansys LS-DYNA是Ansys收购LSTC(Livermore Software Technology Corporation)后得到的求解器,通过Workbench或Mechanical APDL调用。前者配置简单,适合中等非线性问题。后者在极限非线性(爆炸、超高速碰撞、复杂接触)和行业标准材料模型(MAT_***系列)上有压倒性的实绩,是尖端研发的首选。

🧑‍🎓

Abaqus/Explicit和LS-DYNA都是显式法的顶级产品,强项有什么不同?

🎓

好问题。Abaqus/Explicit以Dassault Systèmes的一体化CAE环境(从CAD关联到疲劳分析)为优势,特别在制造业(板金成形、跌落测试、缓冲包装)的应用强。用户界面(工程设计工具)整合得好,容易上手。而LS-DYNA作为"求解器"本身性能更强,特别是"高速碰撞""爆炸""兵器安全(军用规范MIL-STD-810的模拟)""生物伤害分析"等极限复杂、苛刻的物理现象,在建模能力上遥遥领先。另外用户自定义子程序的自由度也远高于Abaqus。

🧑‍🎓

Altair Radioss在碰撞分析上也有名,和LS-DYNA怎么选?

🎓

Radioss也是优秀的显式求解器,在汽车工业(特别是欧洲)广泛使用。选择往往取决于企业历史用途和统合环境的依赖度。若已在主力使用Altair HyperWorks生态(HyperMesh、HyperView、OptiStruct),自然倾向Radioss。技术上,Radioss的"块"并行架构对超大规模模型的并行效率据称更优。而LS-DYNA的材料模型数量(约300种)和对极端物理现象(电磁成形、ICF等)的支持范围更广。

🧑‍🎓

免费开源求解器(如CalculiX、Code_Aster)能代替LS-DYNA吗?

🎓

用途决定一切。CalculiX和Code_Aster主要是隐式法,显式法实现不成熟或缺失。复杂接触算法、高级材料破坏模型、ALE和SPH等多物理手段的完成度,以及验证和实绩的积累,都相差数量级。研究基础冲击问题可用开源工具试试。但若要通过汽车碰撞安全法规认证(如Euro NCAP)或产品可信度担保,业界标准的LS-DYNA、Abaqus、Radioss是必选。

LS-DYNA故障排除

常见错误及对策

🧑‍🎓

计算中出现"negative volume in solid element"错误就停止了。这是什么情况,怎么修?

🎓

这表示单元被极度压缩或剪切,雅可比行列式(体积的度量)变成负数,数值上已破裂。原因主要三个:

1. **过度变形**:材料模型不能表现现实的大变形(如缺少破坏准则)。 2. **不当接触设置**:部件间异常大的力导致单元被压扁。 3. **网格品质差**:初始就有宽纵比极端的单元。 对策是,先用d3plot可视化问题单元的历史,确定发生地点。加入*MAT_ADD_EROSION让材料破坏而删除单元,降低接触的惩罚刚性(*CONTACT的SFS、SFM),或在问题区域强化网格划分。
🧑‍🎓

"termination due to shell element thinning"是什么意思?壳单元厚度变成零了?

🎓

正是。壳单元的厚度相对初始厚度降到允许比率(默认0.1或0.01等)以下就停止计算。这其实是安全装置,物理上对应材料断裂、剥离。如果是正当的现象(如板金被拉断),改*CONTROL_SHELL的IRNXX参数自动删除变薄的单元。如果不切实际地变薄,要检视材料模型(硬化规律)或网格尺寸,缓和局部集中变形。

文章评价
感谢您的回答!
有帮助
想要
更详细
报告
错误
有帮助
0
想要更详细
0
报告错误
0
编写者:NovaSolver贡献者
匿名工程师 & AI — 网站地图
查看个人资料