3. 材料模型选择与定义

跌落冲击仿真的可信度，很大程度上不取决于网格有多细，而取决于材料模型是否选对。跌落过程通常伴随大变形、接触碰撞、局部塑性、应变率强化甚至材料失效，如果仍用一条线弹性曲线描述全部行为，计算结果往往表现为“应力虚高、变形偏小、能量不守恒”。本节围绕五类常用材料模型，说明它们各自适合描述什么、在 Ansys Mechanical 中如何落地，以及参数从哪里来。

一、先判断材料在冲击中处于什么状态

选择材料模型前，先回答三个问题：结构是否允许永久变形？变形量是否超过 5%～10%？应变率是否显著高于准静态？例如手机中框铝合金在 1.5 m 跌落中局部进入塑性，应变率可达 10²～10³ s⁻¹；而橡胶缓冲垫属于超弹性大变形；发泡 EPE 则依靠泡孔坍塌吸能。判断清楚后，模型选择基本就确定了。

二、五类模型的适用性与关键参数

线弹性模型只定义弹性模量 E 和泊松比 ν，密度 ρ 用于惯性。它适合两类场景：一是刚性体或远未屈服的部件，如玻璃盖板在低能量冲击下的初始响应；二是作为复杂模型的“占位”材料，用于快速验证接触和边界。跌落中若某区域应力已超过屈服强度，继续用线弹性会高估回弹、低估吸能，因此它不适合作为主承力金属件的最终模型。

弹塑性模型是金属件跌落的主力。Ansys Mechanical 中通过 Engineering Data 添加 Plasticity→Bilinear Isotropic Hardening 或 Multilinear Isotropic Hardening。双线性只需屈服强度 σy 和切线模量 Et；多线性则输入真实应力—塑性应变曲线。跌落常用多线性，因为高应变率下曲线明显上翘。注意输入的是真实应力与塑性应变，若只有工程应力—应变曲线，需先转换：σ_true=σ_eng(1+ε_eng)，ε_pl=ln(1+ε_eng)−σ_true/E。

超弹性模型用于橡胶、硅胶、TPU 等。Ansys 提供 Mooney-Rivlin、Neo-Hookean、Ogden 等。跌落中橡胶缓冲垫常经历 50%～200% 应变，必须用超弹性而非线弹性。参数通常由单轴拉伸、等双轴拉伸、平面剪切试验拟合得到；若只有单轴数据，可用 Neo-Hookean 两参数近似，但大变形下误差较大。设置时在 Engineering Data 选择 Hyperelastic→Ogden，输入 μi、αi、Di，其中 Di 控制可压缩性，橡胶近似不可压时 Di 取小值。

泡沫材料（EPE、EPS、PU 泡沫）的典型特征是应力—应变曲线分三段：线弹性、平台坍塌、密实化。Ansys 中可用 Crushable Foam 配合 Plasticity，或直接用 Multilinear Isotropic Hardening 输入平台段。关键参数是屈服应力—塑性应变曲线和拉伸/压缩强度比。泡沫吸能主要发生在平台段，若只给弹性模量，仿真会完全丢失缓冲效果。

应变率相关模型用于应变率敏感材料，如铝合金、PC/ABS、部分钢材。Ansys 中常用 Johnson-Cook 或 Cowper-Symonds。Johnson-Cook 屈服应力形式为 σ=(A+Bεⁿ)(1+C ln ε̇*)(1−T*ᵐ)，其中 A、B、n 描述硬化，C 描述率敏感，m 描述温度软化。跌落仿真若忽略率效应，铝合金在 10³ s⁻¹ 下的屈服强度可能被低估 20%～40%，导致塑性区偏大、吸能偏多。设置时在 Engineering Data 选择 Plasticity→Johnson-Cook，逐项填入 A、B、n、C、m 及参考应变率。

三、在 Ansys Mechanical 中的具体设置方法

打开 Engineering Data，新建材料后按以下路径操作。线弹性：Physical Properties→Density 填 ρ，Linear Elastic→Isotropic Elasticity 填 E、ν。弹塑性：在 Plasticity 下选 Bilinear Isotropic Hardening，Yield Strength 填 σy，Tangent Modulus 填 Et；若选 Multilinear Isotropic Hardening，则在表格中逐行输入 Plastic Strain 与 Stress，第一行塑性应变必须为 0，应力等于 σy。超弹性：Hyperelastic→Ogden，输入 μ1、α1、D1，必要时增加 μ2、α2、D2。泡沫：Plasticity→Crushable Foam，或直接用 Multilinear Isotropic Hardening 输入平台段。应变率相关：Plasticity→Johnson-Cook，依次填 A、B、n、C、m，Reference Strain Rate 通常取 1 s⁻¹。

设置完成后，在 Model 中把材料指派给对应几何体。跌落分析还需注意：密度必须正确，否则惯性力错误；接触对中目标面与接触面刚度差异大时，可适当调整法向刚度因子；求解设置中打开大变形（Large Deflection），并控制时间步长使冲击阶段至少有 10～20 个时间点。

四、参数来源与验证

材料参数优先来自材料供应商数据表、文献或材料试验。若只有工程应力—应变曲线，按前述公式转换。验证方法：先做单轴拉伸仿真，与试验曲线对比；再做简单落球冲击，检查接触力峰值和塑性区是否合理。若仿真中塑性应变远超材料延伸率，应检查是否遗漏失效准则或应变率强化。

最后给出一个可参考的参数配置示例，便于直接对照填写：

以上示例值属于教学合理假设，实际项目应以试验或供应商数据为准。材料模型选对并正确输入参数后，跌落仿真的冲击力峰值、塑性区和吸能分配才具备工程参考价值。