1. 课程导入与LS-DYNA概览

碰撞与冲击是工程结构最常见的极端载荷之一。汽车正面碰撞、手机跌落、电池包挤压、飞机鸟撞、防护头盔吸能,这些场景的共同特点是:作用时间极短(毫秒级)、变形极大(材料进入塑性甚至失效)、接触状态不断变化(部件之间反复碰撞与分离)。用静力学方法分析这类问题会得到完全错误的结论,因为惯性效应和应力波传播才是主导因素。本课程的目标,就是让学习者掌握用LS-DYNA求解这类瞬态强非线性问题的完整能力。

学习目标:从"会用软件"到"会判断结果"

完成本课程后,你应当能够独立完成以下任务:第一,根据问题特征判断是否适合显式动力学求解,并合理选择单位制、时间步长与求解控制参数;第二,建立包含壳单元、实体单元、焊点与连接关系的碰撞有限元模型,正确设置接触、材料本构与边界条件;第三,读懂d3plot结果文件中的能量曲线、接触力曲线与变形云图,判断仿真是否物理可信;第四,识别常见的数值问题——如负体积、沙漏能过大、初始穿透、质量缩放过度——并给出针对性的修正方案。

需要强调的是,本课程不追求"点完按钮出结果",而是强调结果判读能力。一个能跑通的模型不等于一个正确的模型。能量守恒是否满足、沙漏能是否控制在总能量的5%以内、接触力曲线是否出现非物理振荡,这些才是判断仿真可信度的关键依据。

LS-DYNA在显式动力学中的定位

LS-DYNA由Livermore Software Technology Corporation开发,现属Ansys旗下,是显式动力学领域事实上的工业标准。它的核心求解器采用中心差分法进行时间积分,无需组装整体刚度矩阵,也不进行迭代求解,因此单步计算成本极低,适合处理大规模接触与高度非线性问题。

与隐式求解器(如Abaqus/Standard、Ansys Mechanical)相比,显式方法的时间步长受Courant条件限制:稳定时间步长近似为最小单元特征长度除以材料声速。这意味着网格越小,时间步越短,计算步数越多。以钢材为例,声速约5000 m/s,若最小单元尺寸为5 mm,则稳定时间步约为1微秒。模拟100毫秒的碰撞过程需要约10万步,这正是显式求解适合短时瞬态问题的原因。

LS-DYNA的另一个优势是丰富的接触算法库。它提供单面接触、面面接触、点面接触等多种类型,支持自动搜索接触对、处理自接触与侵蚀接触,这对碰撞中部件大变形、材料失效后新表面暴露的场景至关重要。

LS-DYNA 显式动力学求解流程与关键控制点 前处理 几何清理 / 网格划分 材料 / 单元 / 连接 接触 / 边界 / 初始条件 求解 中心差分时间积分 时间步 / 质量缩放 接触搜索 / 沙漏控制 后处理 能量曲线 / 接触力 变形 / 应力云图 结果可信度判断 求解流程中的关键控制点 单位制一致性 mm-ton-s-N-MPa 或 m-kg-s-N-Pa 时间步控制 *CONTROL_TIMESTEP TSSFAC / DT2MS 接触定义 *CONTACT_AUTOMATIC SOFT / SURFACE_TO_SURFACE 沙漏与能量 *CONTROL_HOURGLASS 沙漏能 < 5% 总能量 结果可信度判读依据 能量守恒检查 总能量 / 动能 / 内能 / 沙漏能 接触力曲线 峰值 / 振荡 / 穿透量 变形与失效模式 塑性铰 / 撕裂 / 侵入量

软件界面与求解流程总览

LS-DYNA本身是一个求解器,通常配合前处理器(如LS-PrePost、HyperMesh、Ansys Workbench)和后处理器(LS-PrePost)使用。LS-PrePost是LSTC官方免费工具,既能做简单建模,也能完成全部后处理工作,建议初学者优先掌握。

一个典型的碰撞分析流程包含以下环节:

  1. 几何准备与网格划分:清理CAD模型中的小圆角、小孔等对碰撞响应影响不大的特征,划分壳单元或实体单元。碰撞分析中壳单元使用尤为广泛,因为薄壁结构在碰撞中主要通过弯曲和膜力吸能。
  2. 材料与单元属性定义:为每个部件指定材料本构(如*MAT_PIECEWISE_LINEAR_PLASTICITY用于金属塑性,*MAT_PLASTIC_KINEMATIC用于考虑应变率效应)和单元公式(如壳单元ELFORM=2为Belytschko-Tsay,计算效率高)。
  3. 连接与接触设置:焊点用*CONSTRAINED_SPOTWELD或*MAT_SPOTWELD模拟,螺栓连接可用*CONSTRAINED_NODAL_RIGID_BODY。接触通常用*CONTACT_AUTOMATIC_SINGLE_SURFACE处理自接触,用*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE处理部件间接触。
  4. 边界条件与初始速度:固定刚性墙或约束节点,通过*INITIAL_VELOCITY给运动部件赋初速度。例如整车正面碰撞中,车辆以50 km/h(约13.89 m/s)撞击刚性墙。
  5. 求解控制与输出设置:设置终止时间(*CONTROL_TERMINATION)、时间步缩放因子(TSSFAC,默认0.9)、输出间隔(*DATABASE_BINARY_D3PLOT),提交求解器计算。
  6. 结果查看与判断:在LS-PrePost中打开d3plot文件,查看能量曲线、接触力、变形动画,判断仿真是否可信。
初学者常见误区:直接使用CAD原始模型划分网格,导致大量小单元使时间步急剧减小;忽略单位制统一,导致材料参数数量级错误;接触定义时未检查初始穿透,使计算一开始就出现异常接触力。

本课程的组织方式

后续章节将按照"建模—求解—判读—修正"的主线展开。我们会用具体的教学示例贯穿:一个薄壁方管以10 m/s轴向撞击刚性墙,通过这个简单但典型的模型,逐步引入单元选择、材料定义、接触设置、能量分析与结果判读。所有示例参数均为教学假设值,用于说明方法,不代表特定产品的实测数据。

学习过程中建议同步打开LS-PrePost,跟随章节内容动手操作。碰撞仿真的很多细节——比如接触厚度的设置、沙漏控制的类型选择——只有亲手做过、看过结果差异,才能真正理解其影响。