落下冲击试验仿真
落下冲击试验的理论基础
自由落下的力学
老师，我想用FEM预测智能手机跌落时的冲击，能从基础的落下物理开始讲解吗？
很好的问题。落下试验的出发点是忽略空气阻力的自由落下。从高度 $h$ 静静释放的物体，在撞击地面瞬间的速度可由能量守恒导出：
例如，智能手机从 1.2 m 高度（从口袋掏出的高度）跌落时：
这个速度称为"初速度"，是落下FEM的基本设置参数。质量 $m$ 在公式中消去了，但碰撞后的变形和应力显然受其影响。
4.85 m/s，换成时速大概17 km/h左右吧。自行车碰撞地面的速度，想想就很快…。
没错，而且接触时间极短。智能手机的角与混凝土地面碰撞时，接触时间仅为 0.3〜0.5 ms。在这极短的时间内，动量变化 $\Delta p = mv$ 产生的平均冲击力为：
200克的物体承受约2400 N的力，对应 约1200 G的加速度。这就是摧毁电子器件焊接接头的原因。
冲突的接触力学
冲突瞬间的接触力用什么模型来计算？简单地用弹簧模型可以吗？
弹性碰撞时，赫兹接触理论（Hertz Contact Theory）是基础。两个弹性体接触时，接触力 $F$ 与压入量（压陷量）$\delta$ 的关系为：
其中 $K_H$ 是赫兹接触刚度，由两个物体的弹性模量、泊松比和曲率半径决定。但落下试验中会发生 塑性变形和破损，仅用赫兹理论不够。FEM在单元级直接计算弹塑性本构关系，得到真实的非线性接触力。
FEM的接触具体怎么处理？
主要有 惩罚法（Penalty Method）和 拉格朗日乘数法（Lagrange Multiplier Method）两种。实务中的落下分析绝大多数用惩罚法。原理很简单：两个面试图穿透时，虚拟弹簧推回：
$k_p$ 是惩罚刚度，$g_n$ 是间隙（负值表示穿透量）。$k_p$ 太小会产品"穿过"地面，太大会让时间步长极端减小导致计算爆炸。先用默认值跑，再微调是实务做法。
能量收支和吸收机制
落下的能量都到哪去了？没坏的话有时候不反弹，有时候会弹起来…。
极其重要的问题。碰撞前的全能量是动能 $E_k = \frac{1}{2}mv^2$，碰撞后分配到多种形式：
- 弹性应变能 $E_{elastic}$ — 变形后能恢复的部分。是反弹的能量来源
- 塑性耗散能 $E_{plastic}$ — 永久变形消耗。筐体凹陷、裂纹
- 摩擦能 $E_{friction}$ — 接触面滑动消耗
- 阻尼能 $E_{damping}$ — 材料粘性吸收
- 反弹能 $E_{rebound}$ — 弹起后的动能
用反发系数（COR）$e$ 表示，则 $E_{rebound} = e^2 \cdot E_k$。智能手机玻璃面撞瓷砖时，$e \approx 0.2\text{〜}0.4$，大部分能量被塑性变形和内部阻尼消耗。
落下姿态和应力集中
听说落下角度对结果影响很大，哪种姿态最危险？
落下姿态影响巨大。相同高度落下，平面接触和角接触的局部应力可能相差 10倍以上：
| 落下姿态 | 初始接触面积 | 峰值加速度 | 破坏风险 | 典型损伤 |
|---|---|---|---|---|
| 平面落下（Flat） | 大（全面） | 200〜500 G | 低 | 液晶破裂、基板挠度 |
| 边落下（Edge） | 中（线状） | 500〜1,500 G | 中 | 框架变形、按钮破损 |
| 角落下（Corner） | 极小（点状） | 1,000〜5,000 G | 最大 | 筐体裂纹、IC剥离 |
IEC 60068-2-31要求 6个面、12条边、8个角共26个方向的落下试验。FEM至少要做面、边、角三个代表性方向，余量评估的话理想是全26个方向。
26个方向，一次30分钟的话得13小时…。现实吗？
这时候HPC集群就派上用场了。三星、苹果这样的大厂用几十台工作站并行跑26个方向，一天内全出结果。中小企业的做法是先聚焦角落下3个方向+平面落下1个方向共4个工况，找出最严苛条件，这样更现实。
应变速率效应
教科书上说冲击问题中应变速率影响很大，落下试验也要考虑吗？
必须的。落下冲击的典型应变速率是 $10^1 \sim 10^3$ /s，比准静态试验（$10^{-3}$ /s）高万倍以上。很多金属和塑料在高应变速率下 屈服应力显著提升。代表性模型是Cowper-Symonds本构：
$\sigma_0$ 是准静态屈服应力，$D$ 和 $q$ 是材料常数。比如聚碳酸酯（PC）：$D = 10$ /s、$q = 2$，应变速率100 /s时屈服应力约为准静态的3倍。忽视这点会导致FEM高估变形，误判"会坏"。
iPhone落下试验的秘辛
初代iPhone开发时，苹果每天让数十台试作机跌落，据传说。iPhone 12之后采用Ceramic Shield玻璃，落下耐性提升4倍，这个"4倍"背后是LS-DYNA数千工况参数化分析的结果。用内聚单元模拟玻璃微裂纹扩展，把边缘倒角曲率精细优化到0.1 mm级别。
落下冲击试验的数值计算方法
显式法的时间积分
落下分析几乎都用显式法，隐式法不行吗？
落下冲击是 数毫秒超短时间现象，同时发生接触、大变形、材料塑性化。隐式法每步都要求解联立方程，接触条件一变收敛就困难。显式法基于中心差分：
质量矩阵 $M$ 对角化（lumped mass）后，无需求解方程，每个节点加速度独立算。非线性很强的问题也能稳定推进。缺点是条件稳定，安定时间步为：
钢的1 mm六面体单元约 $\Delta t \approx 0.2$ μs。5 ms冲击现象需约25,000步的计算量。
25,000步看起来很多，但每步轻的话整体也快吧…。
完全对。智能手机规模模型（50万单元）用16核LS-DYNA跑下来 30分钟到1小时就完成了。同样问题用隐式法，接触切换导致不收敛，要花多久谁也说不好。落下分析用显式法不是"快"，是"能解"。
接触算法
接触设置是最难的吧。具体怎么做？
落下分析的接触分两类：
- 产品 vs 地面 — 外部接触。把地面设为刚体面就不需要网格，效率高。混凝土或瓷砖非常硬，用刚体近似完全没问题
- 产品内部自接触 — 筐体变形时内部基板或电池接触。常被忽视但对电子产品很关键
LS-DYNA用 *CONTACT_AUTOMATIC_SINGLE_SURFACE 是标准，全外表面自动检测接触。惩罚刚度比例因子（SFS/SFM）默认1.0先跑，看穿透情况再上调到2.0〜5.0。
LS-DYNA 输入示例
LS-DYNA的具体设置能示范一下吗？1 m高平面落下的情况。
用mm-ms-ton制式。1 m落下冲速度 $v = \sqrt{2 \times 9.81 \times 1.0} = 4.43$ m/s = 4.43 mm/ms。这样设置：
$ --- 初速度（产品组全体向下）---
*INITIAL_VELOCITY_SET
$ nsid vx vy vz
1, 0.0, 0.0, -4.43
$ --- 重力（追踪反弹时需要）---
*LOAD_BODY_Z
$ lcid sf
1, 1, 9.81e-3 $ mm/ms^2
$ --- 刚体壁（无限硬地面）---
*RIGIDWALL_PLANAR
$ xt yt zt xh yh zh
0., 0., 0., 0., 0., 1.
$ --- 接触（自动单面）---
*CONTACT_AUTOMATIC_SINGLE_SURFACE
$ ssid msid sstyp mstyp sboxid mboxid spr mpr
1, 0, 3, 0, 0, 0, 1, 1
$ fs fd dc vc
0.3, 0.2, 0., 0.
$ --- 时间控制 ---
*CONTROL_TIMESTEP
$ dtinit tssfac
0., 0.9关键是 TSSFAC=0.9 用稳定时间步的90%。通常0.9就可以，接触激烈可降到0.67。
Abaqus/Explicit 输入示例
Abaqus/Explicit怎么写？公司是用Abaqus的…。
Abaqus一般用SI单位。相同的1 m落下：
** --- 初速度 ---
*INITIAL CONDITIONS, TYPE=VELOCITY
product_set, 3, -4.43
** --- 重力 ---
*DLOAD
product_set, GRAV, 9.81, 0., 0., -1.
** --- 刚体地面（分析型刚体面）---
*RIGID BODY, REF NODE=floor_rp,
ANALYTICAL SURFACE=floor_surf
*SURFACE, TYPE=PLANAR, NAME=floor_surf
DATA LINE: 0., 0., 0., 0., 0., 1.
** --- 通用接触（推荐）---
*CONTACT, OP=NEW
*CONTACT INCLUSIONS, ALL EXTERIOR
** --- 步（5 ms分析）---
*DYNAMIC, EXPLICIT
, 0.005Abaqus的 General Contact 很方便，全外表自动标记接触候选。分析型刚体面定义地面，省去网格。效果和LS-DYNA的 AUTOMATIC_SINGLE_SURFACE 一样。
单元选择和网格策略
单元类型选固体还是壳体？
按零件厚度和变形性质选：
| 零件 | 推荐单元 | LS-DYNA | Abaqus | 理由 |
|---|---|---|---|---|
| 筐体（1〜3 mm厚） | 壳体 | ELFORM=2（BT） | S4R | 板厚方向5积分点保证弯曲精度 |
| 基板（0.8〜1.6 mm） | 壳体 | ELFORM=16（完全） | S4 | 反翘变形精确捕捉 |
| IC/BGA（焊点） | 固体 | ELFORM=1（低减） | C3D8R | 需要三维应力评估 |
| 橡胶垫片 | 固体 | ELFORM=-1（完全） | C3D8H（混合） | 回避非压缩橡胶体积锁定 |
| 泡沫缓冲 | 固体 | ELFORM=1 | C3D8R | 适应大压缩变形 |
网格尺寸：接触面附近 0.5〜1.0 mm，其他2〜5 mm。壳体沿厚度至少5积分点，固体厚度方向至少3层。
低减单元出沙漏模式了怎么办？
沙漏能量超内部能的 5%警戒，超10%必须重做。对策三选：
- 沙漏控制调参 — LS-DYNA用
*HOURGLASS试IHQ=4（Flanagan-Belytschko粘性）或IHQ=5（刚性） - 改完全积分单元 — 沙漏彻底消除，但计算成本2〜3倍
- 网格优化 — 扭曲单元（纵横比>5）是沙漏温床，改善单元品质能大幅缓解
接触刚度调参是手艺活
落下分析现场，"惩罚刚度调整"是经验差异显现的地方。Abaqus通用接触默认用有效单元刚度的10倍设惩罚，但带橡胶垫片的筐体容易过刚性，冲击力峰值成了实测的2倍。反过来降刚度又穿透。某大家电厂CAE团队规定：新材料用单元压入模拟预先做"接触标定"，每次材料变都验证一遍接触参数。
落下冲击试验的实务应用
行业标准和试验条件
实务中落下试验按什么标准定？产品不同吧。
按产品类型汇总主要标准：
| 产品 | 标准 | 落下高度 | 地面 | 姿态 |
|---|---|---|---|---|
| 智能手机 | 各企业内标 | 1.0〜1.8 m | 混凝土/钢板 | 6面+4角 |
| 笔记本电脑 | MIL-STD-810H Method 516.8 | 0.75 m（运作中） | 胶合板 | 26方向 |
| 平板电脑 | IEC 60068-2-31 | 0.5〜1.0 m | 钢板 | 6面+12边+8角 |
| 包装货物 | ISTA 3A / ASTM D5276 | 0.3〜1.0 m | 混凝土 | 角+边+面 |
| 军用电子器材 | MIL-STD-810H Method 516.8 | 1.22 m | 胶合板/混凝土 | 26方向 |
| 医疗器械 | IEC 60601-1 | 0.5〜1.0 m | 混凝土 | 厂商规定 |
| EV电池包 | UN 38.3 Test T.7 | 1.0 m | 混凝土 | 平面落下 |
关键是 地面材质。混凝土比钢板硬得多，冲击力峰值可差1.5〜2倍。FEM要准确设地面弹性模量，或干脆用刚体面。
建模工作流程
收到CAD，从建模到出结果的整个流程是什么？
标准工作流如下：
- CAD清理 — 螺纹孔、小倒角（R < 0.3 mm）、微小台阶简化。但冲击受力的倒角 一定保留（应力集中关键）
- 网格生成 — 冲击面0.5〜1.0 mm，其他2〜5 mm。板金筐用壳体，焊点用固体
- 材料定义 — 弹塑性+应变速率依存。树脂还要加温度依存
- 接触定义 — 自动单面+产品-地面接触对。摩擦系数用实测值（通常0.2〜0.5）
- 初速度、重力设置 — $v = \sqrt{2gh}$ 竖直向下
- 输出设置 — 加速度（重心、IC位置）、应力、能量高频输出（0.01〜0.1 ms间隔）
- 求解运算 — 冲击+反弹周期（通常5〜20 ms）
- 结果验证 — 能量守恒检查 → 加速度、应力评价 → 变形动画确认
结果评价和合否判定
分析结果怎么判断"设计合格"？看什么？
合否有四个主要评估轴：
- 结构完整性 — 筐体有没有裂。FEM看塑性应变是否超材料破断应变（如PC约0.10〜0.15）
- 冲击加速度 — 内部IC、BGA受的加速度在规范内吗。典型BGA焊点许容约1,500〜3,000 G（脉宽0.5 ms）
- 永久变形 — 筐体间隙、高度变化在许容范围。如"落后间隙变化 < 0.3 mm"
- 功能部件完好 — 按钮、接头、铰链等活动部件能否正常动作
"不坏=合格"不对，正确的是"坏的方式在许容内"。设计前要和产品团队明确界线：允许些微凹陷，但电池短路、液晶破裂NG。
能量平衡验证
前辈说"先看能量"，具体怎么看？
落下分析品质保证的 最关键指标。检查项目：
- 全能量守恒 — $E_{total}(t)$ 对初值（$= E_k = \frac{1}{2}mv^2$）变动在 ±5%以内。大幅波动说时间步太大或接触设置不当
- 沙漏能比 — $E_{hourglass} / E_{internal} < 5\%$。超这个是低减单元"乱跑"的证据
- 冲击前后能量转移 — 冲击前动能100%，冲击中转内部能（弹性+塑性），弹起后回到动能。物理合理吗？
- 负内部能 — 如果出现说单元"凭空生能"，完全破裂。立即停止分析
亚马逊纸箱包装优化
亚马逊2019年启用"APASS（Amazon Packaging Academic & Scientific Services）"包装落下仿真系统，年执行数十万工况虚拟试验。用LS-DYNA的CRUSHABLE_FOAM材料模型再现纸板缓冲特性，用DOE（实验计划法）自动优化产品与包装组合。结果包材成本降15%，物流破损率降3成。
落下冲击试验的软件比较
主要求解器的特点
落下分析有什么求解器可选？各有什么优劣？
主要显式法求解器对比：
| 求解器 | 开发方 | 优势 | 落下分析特色 | 许可形式 |
|---|---|---|---|---|
| LS-DYNA | ANSYS（原LSTC） | 冲撞分析业界标准 | RIGIDWALL、CONSTRAINED_NODE_SET、SPOTWELD单元完备 | 节点锁/浮动 |
| Abaqus/Explicit | Dassault Systemes | 材料库丰富 | General Contact易用、超弹性材料精度高 | 令牌/浮动 |
| Ansys Explicit (AUTODYN) | ANSYS | Workbench集成GUI | Drop Test Analysis System模板配置简便 | 令牌 |
| RADIOSS | Altair | 优化联动 | 与HyperStudy集成实现落下设计优化 | HWU（共通令牌） |
| PAM-CRASH | ESI Group | 流程仿真 | 成型→落下联成分析（残余应力继承） | 浮动 |
用途别选型指南
公司主要做智能手机落下试验，选哪个最好？
手机用 LS-DYNA或Abaqus/Explicit 两强。按用途分：
- 电子产品（手机、PC） → LS-DYNA首选。SPOTWELD单元高效建模螺钉接合，显式计算快。三星、联想、小米都用
- 橡胶树脂多的产品（保护壳、医疗器械） → Abaqus/Explicit。Mullins效应（橡胶应力软化）、Arruda-Boyce超弹性标配
- 初学（重GUI） → Ansys Explicit Dynamics。Drop Test Analysis System模板，输入落高和姿态就自动设好
- 包材优化 → LS-DYNA + LS-OPT。CRUSHABLE_FOAM泡沫模型+优化工具完善
- 连续分析（成型→落下） → PAM-CRASH。注塑残余应力直接继承
开源选项没有吗？许可费很贵…。
有个 OpenRadioss（Altair 2022年开源RADIOSS源码）。LS-DYNA兼容关键字也能读，落下基本功能齐全。但前处理、后处理、技术支持跟商用版差距大，教学和原型验证用还行，产品认证级别的生产解析还是商用求解器更放心。
ThinkPad的MIL-STD-810H认证背后
联想ThinkPad通过MIL-STD-810H Method 516.8落下试验的名声，背后是LS-DYNA对26方向0.75 m落下的全仿真。镁合金框架的壁厚和加强筋通过FEM预优化到90%以上合格率，才进实物试验。从设计到认证的"前期装载"策略使开发周期降为传统的一半。
落下冲击试验的故障排除
穿透和接触不良
模型跑起来，产品就穿过地面飞出去了…。啥情况？
最常见的落下分析故障。逐项检查：
- 缺接触定义 — 产品外面和地面的接触对没定。LS-DYNA检查
*CONTACT卡有没有，Abaqus检查*CONTACT INCLUSIONS - 法线反向 — 壳体法线朝内，接触检测反了。LS-DYNA看
*CONTACT的SST/MST，Abaqus在CAE看法线可视化 - 初始穿透 — 网格阶段产品和地面就重叠了。留1〜2 mm间隙，再赋速度
- 惩罚刚度不足 — 默认值穿透就上调比例因子到2〜5倍。LS-DYNA的
SFS/SFM，Abaqus的*CONTACT PROPERTY的PENALTY
加速度尖峰
加速度时间曲线上有物理上不可能的尖峰（10万G超）。滤波能消除吗？
滤波消除不行。得找根本原因：
- 接触刚度太高 — 接触瞬间反力过大。降低比例因子（0.1〜0.5）
- 节点接触-离脱反复 — 打颤现象。加接触阻尼参数，LS-DYNA的
VDC、Abaqus的DAMPING设20% - 质量缩放过度 — 用质量放大扩大时间步，反而改变冲击力。额外加的质量超全体5%以上的话就是问题
- 输出加速度未滤波 — 对比实测数据要用SAE CFC滤波（CFC-600或CFC-1000）
角落下的单元崩溃
角落下时碰撞后单元奇怪变形，"negative volume"错误停掉了。
"negative volume"是六面体单元翻转（雅可比行列式负值）。角落下接触极小，局部应变爆炸。对策：
- 细化角部网格 — 冲击点周围细到0.3〜0.5 mm。代价是时间步也缩短，计算成本增
- 单元删除（侵蚀） — 达到破断应变的单元自动移除。LS-DYNA的
*MAT_ADD_EROSION、Abaqus的*ELEMENT DELETION。要确保删除单元数 < 全体1%以保质量、能量 - 加倒角 — CAD里角上加0.5〜1.0 mm倒角。实产品也有，物理合理
- ALE法 — 用任意拉格朗日-欧拉法在极变形区。单元跟随变形重新划分，回避翻转
能量发散
全能量曲线碰撞后一直涨，啥毛病？
能量被"生成"了，说明 数值上有能量注入。这是致命bug，结果不可信。原因：
- 接触不稳定 — 节点在接触面来回振荡，惩罚力"泵送"能量。降接触刚度或加阻尼解决
- 沙漏模式失控 — 低减单元零能量模式成长。改沙漏控制参数（IHQ）或换完全积分
- 时间步超稳定极限 — 质量缩放扩大时间步。降
TSSFAC到0.67试试
必须让 能量变动收束到初始动能的±5%，否则应力加速度无信度。
与实验相关性差
能量守恒OK，但加速度波形和实测数据完全对不上。脉宽、幅度都不对。
多因素耦合。逐项检查清单：
- 材料应变速率特性 — 只用准静态参数跑的话结果就偏。要加Cowper-Symonds参数
- 缓冲材（泡沫、橡胶）的压缩曲线 — 测试机实测应力-应变，不能用目录值。发泡料批间波动大
- 接触摩擦系数 — 摩擦大脉宽短峰值高。没实测就用敏感性分析估范围
- 传感器位置对应 — FEM输出节点和试验传感器位置一致吗？一般对比重心位置
- 滤波一致性 — 试验数据和FEM数据都用同一SAE CFC滤波吗？
- 落下姿态偏差 — 试验很难完全平面，通常有几度倾角。FEM也扫±2°的敏感性
调试铁律：从简单开始
结果不合时，资深工程师先做"加法调试"：用均质直方体（铝块）落刚体面，解析解能对上吗？不对就接触设置根本错。再试空筐体。最后加内部器件。"减法调试"比加全部再删更高效，复杂问题尤其管用。