1. 跌落冲击仿真概述

跌落冲击是消费电子、家电、汽车电子、包装运输等行业最常见的结构失效工况之一。手机从口袋滑落到地面、路由器在运输途中被抛掷、电池包在整车碰撞中受到挤压,本质上都是结构在极短时间内承受巨大加速度或位移冲击。Ansys Mechanical 提供了完整的显式动力学求解能力,可以在设计阶段预测这类瞬态响应,从而减少实物跌落试验次数、缩短迭代周期。

1.1 跌落冲击的物理过程

一次典型的跌落冲击可以拆成三个阶段:自由落体、接触碰撞、振荡衰减。以一台 200 g 的手机从 1.2 m 高度跌落到刚性地面为例,落地瞬间速度约为 4.85 m/s(由 v=√(2gh) 估算)。碰撞持续时间通常在 0.5~3 ms 量级,接触力峰值可达数百到数千牛,结构内部产生应力波并以弹性波形式传播。冲击结束后,结构进入阻尼振荡阶段,能量通过材料阻尼、接触摩擦和塑性变形耗散。

理解这一过程的关键在于:冲击载荷不是静态力,而是随时间剧烈变化的瞬态载荷;结构的响应不仅取决于刚度,还强烈依赖质量分布、接触刚度、材料应变率效应和阻尼。因此,跌落仿真必须使用能够处理大变形、接触和瞬态效应的求解器。

① 自由落体 v = √(2gh) 手机 200 g 刚性地面 h = 1.2 m ② 接触碰撞 变形 接触力峰值 500~3000 N 持续时间 0.5~3 ms ③ 振荡衰减 应力波传播与阻尼耗散 能量通过塑性变形/摩擦耗散 关注:最大应力、塑性应变、加速度 t=0 t≈0.5 ms t≈5 ms

1.2 显式动力学与隐式静力学的区别

Ansys Mechanical 中同时存在隐式求解器(Static Structural、Transient Structural)和显式求解器(Explicit Dynamics)。两者在数学本质、时间步长、适用场景上有根本差异。

隐式方法求解的是平衡方程 K·u = F,每个时间步都需要组装刚度矩阵并迭代求解,时间步长可以取得较大,但每一步计算成本高,且在大变形、复杂接触下容易不收敛。显式方法则基于中心差分法,直接由当前时刻的加速度推算下一时刻的速度和位移,不需要求解联立方程组,因此每一步计算极快,但时间步长受 Courant 条件限制,必须小于应力波穿过最小单元的时间。

对比维度 隐式静力学 / 隐式瞬态 显式动力学
求解方程 K·u = F,需迭代 M·a = F_ext − F_int,显式推进
时间步长 较大,可自适应 极小,受 Courant 条件限制
单步成本 高(矩阵分解) 低(向量运算)
大变形/接触 易不收敛 天然适合
典型工况 静载、模态、热应力 跌落、碰撞、冲击、爆炸
结果输出 平衡态应力/位移 时间历程:应力波、加速度、能量

需要强调的是:显式动力学并非“更高级”的求解器,而是针对瞬态、强非线性问题的专用工具。对于静态或准静态问题,使用显式求解反而会因时间步长过小导致计算时间不可接受。选择求解器的依据是物理过程的时间尺度和非线性程度,而不是软件模块的“新旧”。

1.3 Ansys Mechanical 显式动力学模块的适用场景

在 Ansys Mechanical 中,显式动力学通过 Explicit Dynamics 分析系统实现,底层求解器为 Autodyn 或 LS-DYNA(取决于配置)。其典型适用场景包括:

  • 电子产品跌落:手机、平板、路由器从 0.5~1.5 m 高度跌落到不同地面(刚性、混凝土、木地板),关注屏幕玻璃破裂、焊点失效、外壳塑性变形。
  • 包装运输冲击:产品带缓冲泡沫跌落,关注泡沫压缩量、产品承受的加速度峰值是否超过许用值。
  • 汽车碰撞与电池包挤压:低速碰撞、电池包底部球击,关注结构侵入量、电芯是否受损。
  • 冲击波与爆炸:压力容器、防护结构在爆炸载荷下的响应。

不适用或需谨慎使用的场景:长期蠕变、疲劳裂纹扩展、纯静态强度校核、模态分析。这些应交给隐式求解器或专用模块。

1.4 显式动力学工作流

在 Workbench 中搭建一个跌落仿真,通常遵循以下流程:

  1. 几何准备:简化非关键特征(倒角、小孔、螺纹),保留影响刚度与质量分布的主要结构。对薄壁件可抽取中面或使用壳单元。
  2. 材料定义:至少输入密度、弹性模量、泊松比。若关注塑性变形,需补充双线性或多线性硬化曲线;应变率敏感材料还需 Cowper-Symonds 或 Johnson-Cook 参数。
  3. 接触设置:跌落仿真中接触是核心。常用 Frictional 或 Frictionless 接触,接触刚度、穿透容差需根据网格尺寸调整。地面通常设为刚性体。
  4. 网格划分:显式求解对网格质量敏感。单元尺寸决定时间步长,通常要求最小单元尺寸不小于 0.5 mm,否则计算时间急剧增加。接触区域需局部加密。
  5. 初始条件与载荷:给定初速度(由跌落高度换算)或直接施加重力场+初始高度。注意方向定义:通常 Z 轴负方向为重力方向。
  6. 求解设置:设定终止时间(通常 3~10 ms)、结果输出频率、能量控制选项。开启沙漏控制以避免零能模式。
  7. 后处理:查看等效应力、塑性应变、加速度时间历程、接触力、能量曲线。判断依据包括:最大应力是否超过材料强度、塑性应变是否超过许用值、加速度峰值是否超过元器件耐受阈值。
示例参数参考:某塑料外壳手机从 1.2 m 跌落,初速度约 4.85 m/s,求解终止时间设为 5 ms,最小单元尺寸 0.8 mm,时间步长约 1.6×10⁻⁷ s,总步数约 31000 步。接触刚度因子取 0.1,沙漏控制采用 Flanagan-Belytschko 刚度形式。这些数值属于教学示例,实际项目需根据模型规模和精度要求调整。
常见误区:把显式动力学当成“万能求解器”,对静态问题也使用显式求解;或者忽略应变率效应,用准静态材料参数直接做跌落仿真,导致应力预测偏低。跌落仿真中材料在高应变率下屈服强度可能提升 20%~50%,必须使用动态材料参数或应变率模型。