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 的建模难点也正在于此——如何选择材料模型、如何标定失效参数、如何控制单元删除后的数值稳定性。

LS-DYNA 失效与断裂分析:从工程问题到计算判据 工程失效问题 碰撞/冲击 高应变率、大变形 挤压/压溃 屈曲、撕裂、接触 爆炸/冲击波 流固耦合、碎片 疲劳/损伤累积 循环载荷、裂纹萌生 复合材料分层 层间失效、渐进损伤 LS-DYNA 建模与判据 材料本构选择 *MAT_PLASTIC_KINEMATIC / *MAT_JOHNSON_COOK / *MAT_CDM 失效准则定义 等效塑性应变 / 应力三轴度 / 损伤变量 单元删除与裂纹模拟 *MAT_ADD_EROSION / 内聚力单元 / XFEM 接触与连接失效 焊点失效 / 螺栓预紧 / 接触穿透控制 求解控制与稳定性 时间步 / 质量缩放 / 沙漏控制 / 能量平衡 结果判读与验证 应力/应变场 峰值位置、分布范围 损伤演化曲线 损伤变量随时间变化 失效模式对比 与试验/标准对照 能量平衡检查 沙漏能、内能、动能 工程结论 安全裕度、改进方向 参数标定与模型修正反馈

1.3 课程目标与学习路径

本课程的目标不是让读者记住所有材料卡片的关键字,而是建立一套可复用的分析流程:从工程失效现象出发,选择合理的材料模型与失效判据,完成建模与求解,并对结果进行工程判读。围绕这一目标,课程按“基础—材料—判据—应用—排障”五个层次展开:

  1. 基础层(第1~4章):失效与断裂的基本概念、LS-DYNA 显式求解框架、单位制与坐标系约定、材料数据来源与试验基础。
  2. 材料层(第5~12章):金属塑性、应变率相关本构、损伤力学、复合材料、混凝土、泡沫与聚合物等材料模型的适用条件与参数标定。
  3. 判据层(第13~18章):等效塑性应变准则、应力三轴度相关准则、损伤演化与单元删除、内聚力模型、XFEM 与裂纹扩展方法。
  4. 应用层(第19~26章):碰撞吸能、电池包挤压、爆炸冲击、管道撕裂、复合材料冲击、焊点与连接失效等典型工程案例。
  5. 排障层(第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由于考虑了应力三轴度,裂纹更倾向于从应力三轴度较高的中心区域萌生,断裂位置更接近实际试件的杯锥状断口。这个示例说明:失效判据不是“越复杂越好”,而是要与材料实际失效机制匹配。后续章节会详细展开各类判据的标定方法。

常见误区:把单元删除当作“断裂模拟”的全部。单元删除只是失效后的数值处理手段,真正的断裂分析需要关注损伤萌生位置、扩展路径、能量耗散以及删除后对整体刚度与接触的影响。

本章完成了课程的问题界定与路径规划。接下来将从显式动力学的基本原理出发,逐步建立失效与断裂分析所需的力学基础与软件操作能力。