1. 课程导入与工程失效问题界定
在结构强度与安全评估中,材料失效与断裂往往不是“算不出来”的问题,而是“算得不对、判得不准”的问题。LS-DYNA 作为显式动力学有限元程序,其核心优势在于处理高应变率、大变形、接触碰撞与多物理场耦合,而这些工况恰恰是材料最容易发生失效与断裂的场景。因此,本课程的第一件事不是立刻打开软件点按钮,而是先界定:我们要解决的工程失效问题究竟是什么,LS-DYNA 在其中扮演什么角色,以及后续 29 章的学习路径如何组织。
1.1 LS-DYNA 在失效与断裂分析中的定位
LS-DYNA 的求解器以显式时间积分(中心差分)为主,配合隐式求解能力,适合求解瞬态、非线性、接触主导的问题。在材料失效与断裂分析中,它的定位可以概括为三个层面:
- 承载失效过程的模拟器:从弹性变形、塑性屈服、损伤演化到单元删除或裂纹扩展,LS-DYNA 通过材料本构与失效准则的组合,把“失效”表达为可计算的力学过程。
- 多物理场耦合的平台:热-力耦合、流-固耦合、爆炸冲击等场景下,失效往往由温度、压力、应变率共同驱动,LS-DYNA 提供对应的材料模型与耦合算法。
- 工程判据的验证工具:通过仿真获得应力三轴度、等效塑性应变、损伤变量等场量,与试验或标准判据对比,支撑结构安全评估。
需要明确的是,LS-DYNA 不会自动“知道”材料何时断裂。断裂判据、失效参数、单元删除阈值都需要用户根据材料试验、标准规范或工程经验来定义。因此,本课程的重点不是软件操作本身,而是如何把工程失效问题转化为可计算的模型与判据。
1.2 典型工程失效案例与问题类型
为了建立直观认识,先看几类在 LS-DYNA 中高频出现的失效问题。它们对应不同的材料行为、失效模式与建模策略:
| 工程场景 | 典型失效模式 | 关键影响因素 | LS-DYNA 关注点 |
|---|---|---|---|
| 汽车碰撞吸能盒 | 塑性屈曲、撕裂、焊缝开裂 | 应变率、应力三轴度、连接方式 | 材料失效准则、焊点失效、接触 |
| 电池包壳体挤压 | 局部凹陷、裂纹萌生、密封失效 | 挤压速度、壳体厚度、材料延性 | 壳单元失效、损伤演化、接触穿透 |
| 航空发动机叶片外物损伤 | 缺口裂纹、高周疲劳、局部断裂 | 冲击速度、缺口敏感性、温度 | 高应变率本构、裂纹扩展、热力耦合 |
| 管道爆炸撕裂 | 延性撕裂、脆性断裂、碎片飞散 | 内压、壁厚、材料韧性、温度 | 流固耦合、失效判据、碎片处理 |
| 复合材料层合板冲击 | 分层、纤维断裂、基体开裂 | 铺层角度、冲击能量、界面强度 | 层间失效、渐进损伤、单元删除 |
这些案例的共同点是:失效不是单一时刻的“开关”,而是损伤累积、局部化、裂纹萌生与扩展的连续过程。LS-DYNA 的建模难点也正在于此——如何选择材料模型、如何标定失效参数、如何控制单元删除后的数值稳定性。
1.3 课程目标与学习路径
本课程的目标不是让读者记住所有材料卡片的关键字,而是建立一套可复用的分析流程:从工程失效现象出发,选择合理的材料模型与失效判据,完成建模与求解,并对结果进行工程判读。围绕这一目标,课程按“基础—材料—判据—应用—排障”五个层次展开:
- 基础层(第1~4章):失效与断裂的基本概念、LS-DYNA 显式求解框架、单位制与坐标系约定、材料数据来源与试验基础。
- 材料层(第5~12章):金属塑性、应变率相关本构、损伤力学、复合材料、混凝土、泡沫与聚合物等材料模型的适用条件与参数标定。
- 判据层(第13~18章):等效塑性应变准则、应力三轴度相关准则、损伤演化与单元删除、内聚力模型、XFEM 与裂纹扩展方法。
- 应用层(第19~26章):碰撞吸能、电池包挤压、爆炸冲击、管道撕裂、复合材料冲击、焊点与连接失效等典型工程案例。
- 排障层(第27~30章):负体积、沙漏、能量异常、接触穿透、单元删除导致的数值不稳定等常见问题的诊断与处理。
1.4 一个最小示例:失效判据如何影响结果
为了说明“判据选择”对结果的影响,考虑一个简单的拉伸试件示例。材料采用 *MAT_PLASTIC_KINEMATIC,弹性模量 210 GPa,泊松比 0.3,屈服强度 350 MPa,切线模量 1.2 GPa。分别设置两种失效方式:
*MAT_PLASTIC_KINEMATIC
$ MID RO E PR SIGY ETAN BETA
1 7.85e-9 210000 0.3 350 1200 0.0
$ 失效方式A:等效塑性应变达到0.25时删除单元
*MAT_ADD_EROSION
$ MID EXCL MXPRES MNEPS EFFEPS VOLEPS NUMFIP NCS
1 0 0 0 0.25 0 1 0
$ 失效方式B:应力三轴度相关准则,等效塑性应变阈值随三轴度变化
$ 此处省略具体曲线定义,仅示意判据差异
在相同网格与加载条件下,方式A会在试件标距段均匀出现单元删除,断裂位置偏随机;方式B由于考虑了应力三轴度,裂纹更倾向于从应力三轴度较高的中心区域萌生,断裂位置更接近实际试件的杯锥状断口。这个示例说明:失效判据不是“越复杂越好”,而是要与材料实际失效机制匹配。后续章节会详细展开各类判据的标定方法。
本章完成了课程的问题界定与路径规划。接下来将从显式动力学的基本原理出发,逐步建立失效与断裂分析所需的力学基础与软件操作能力。