1. 跌落冲击仿真概述
跌落冲击是消费电子、家电、汽车电子、包装运输等行业最常见的结构失效工况之一。手机从口袋滑落到地面、路由器在运输途中被抛掷、电池包在整车碰撞中受到挤压,本质上都是结构在极短时间内承受巨大加速度或位移冲击。Ansys Mechanical 提供了完整的显式动力学求解能力,可以在设计阶段预测这类瞬态响应,从而减少实物跌落试验次数、缩短迭代周期。
1.1 跌落冲击的物理过程
一次典型的跌落冲击可以拆成三个阶段:自由落体、接触碰撞、振荡衰减。以一台 200 g 的手机从 1.2 m 高度跌落到刚性地面为例,落地瞬间速度约为 4.85 m/s(由 v=√(2gh) 估算)。碰撞持续时间通常在 0.5~3 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 中搭建一个跌落仿真,通常遵循以下流程:
- 几何准备:简化非关键特征(倒角、小孔、螺纹),保留影响刚度与质量分布的主要结构。对薄壁件可抽取中面或使用壳单元。
- 材料定义:至少输入密度、弹性模量、泊松比。若关注塑性变形,需补充双线性或多线性硬化曲线;应变率敏感材料还需 Cowper-Symonds 或 Johnson-Cook 参数。
- 接触设置:跌落仿真中接触是核心。常用 Frictional 或 Frictionless 接触,接触刚度、穿透容差需根据网格尺寸调整。地面通常设为刚性体。
- 网格划分:显式求解对网格质量敏感。单元尺寸决定时间步长,通常要求最小单元尺寸不小于 0.5 mm,否则计算时间急剧增加。接触区域需局部加密。
- 初始条件与载荷:给定初速度(由跌落高度换算)或直接施加重力场+初始高度。注意方向定义:通常 Z 轴负方向为重力方向。
- 求解设置:设定终止时间(通常 3~10 ms)、结果输出频率、能量控制选项。开启沙漏控制以避免零能模式。
- 后处理:查看等效应力、塑性应变、加速度时间历程、接触力、能量曲线。判断依据包括:最大应力是否超过材料强度、塑性应变是否超过许用值、加速度峰值是否超过元器件耐受阈值。