3. 零部件冲击跌落仿真建模流程总览

零部件跌落仿真的本质，是把一次真实跌落事件在计算机中"重演"一遍：零件从某个高度以某个姿态落到刚性或柔性地面，经历接触、反弹、塑性变形甚至断裂。要让这次"重演"可信，必须依次完成几何清理、网格划分、材料定义、接触设置、边界条件、求解与后处理七个环节。它们不是彼此独立的操作，而是一条前后耦合的链路——前一步的输入质量直接决定后一步能否成立。本节先建立这条链路的整体视图，明确每一步的输入、输出与关键判断，后续章节再逐项展开。

链路全景与数据传递关系

下图给出七个环节的串联关系，并标注了每一步的核心输入与输出。可以看到，几何清理的输出是"可用于网格划分的干净几何"，网格划分的输出是"带单元与节点的有限元模型"，材料定义把单元与真实材料行为绑定，接触设置定义零件与地面之间的相互作用，边界条件给出初始速度与约束，求解器据此推进时间积分，后处理再把结果翻译成工程判断。任何一环的输入不成立，后续结果都不可信。

各环节的输入、输出与关键判断

几何清理是整条链路的起点。CAD 模型常带有微小圆角、倒角、螺纹、装配间隙和重复面，这些特征在跌落仿真中既不影响刚度，又会把网格尺寸压到极小，导致时间步长骤降。清理的目标是保留承载与吸能结构，删除或简化非承载细节。例如某塑料外壳上的 0.3 mm 装饰圆角，若保留则单元尺寸被迫取 0.15 mm，临界时间步长可能从 1e-6 s 降到 1e-7 s，计算量增加十倍以上。清理后应检查几何是否封闭、有无自由边和重复面，这是后续网格能否生成的前提。

网格划分把几何离散为单元与节点。跌落仿真以显式动力学为主，常用六面体或壳单元，四面体单元在接触和沙漏控制上更敏感。网格尺寸需与最小特征和波传播需求匹配：一般要求沿厚度方向至少 2～3 层单元，接触区域适当加密。划分完成后必须检查单元质量，重点关注雅可比、翘曲度、长宽比和负体积。质量差的单元会直接导致求解中止或结果失真。

材料定义决定零件在冲击下的真实响应。金属常用 *MAT_PLASTIC_KINEMATIC 或 *MAT_JOHNSON_COOK，塑料常用 *MAT_PLASTIC_KINEMATIC 配合失效应变，泡沫或橡胶则用 *MAT_CRUSHABLE_FOAM、*MAT_MOONEY-RIVLIN 等。关键参数包括密度、弹性模量、泊松比、屈服强度、切线模量和失效应变。以某 ABS 外壳为例，可取密度 1.05e-9 t/mm³、弹性模量 2200 MPa、泊松比 0.35、屈服强度 40 MPa、失效应变 0.15。这些数值应来自材料试验或可靠手册，不能凭经验随意填写。

接触设置定义零件与地面、零件与零件之间的相互作用。跌落仿真中地面通常用 *RIGIDWALL_PLANAR 或刚性板模拟，零件自身用 *CONTACT_AUTOMATIC_SINGLE_SURFACE 处理自接触。需要设置静摩擦与动摩擦系数，塑料对钢地面一般取 0.2～0.4。接触刚度、穿透容差和接触厚度是影响结果稳定性的关键参数，设置不当会出现穿透或接触力振荡。

边界条件给出初始速度、姿态和约束。跌落高度 h 与落地速度 v 的关系为 v = √(2gh)，取 g = 9.81 m/s²，h = 1.0 m 时 v ≈ 4.43 m/s。姿态决定冲击位置，是角跌、面跌还是棱跌，直接影响应力分布。求解时间一般取 5～20 ms，需覆盖零件完成首次冲击与主要反弹过程。若只关心首次峰值，可适当缩短。

求解阶段把上述设置组织成完整的 K 文件，提交 LS-DYNA 计算。显式求解的稳定性由临界时间步长控制，Δt ≈ L/c，其中 L 为最小单元特征长度，c 为材料声速。求解过程中应监控能量曲线：总能量应基本守恒，沙漏能应小于内能的 5%，否则需调整网格或沙漏控制。若出现负体积或节点速度异常，通常意味着网格或接触参数需要回溯修改。

后处理把 d3plot、binout 等结果文件翻译成工程判断。常用观察量包括：关键位置的等效应力峰值、塑性应变分布、零件最大变形量、接触力时程曲线以及各能量分量随时间的变化。判断依据应结合设计指标，例如外壳最大应力是否超过材料屈服强度、卡扣是否发生断裂、屏幕区域变形是否超过允许值。若结果异常，应回到对应环节排查，而不是直接修改结果。

链路中的常见耦合与回溯

七个环节并非严格单向。网格尺寸会影响临界时间步长，进而影响求解时间；材料失效应变会影响接触区域的单元删除，进而改变接触行为；接触刚度设置不当会引入额外能量，干扰能量守恒判断。因此实际建模中常需要多次迭代：先粗算一遍确认链路通畅，再逐步细化网格和材料参数。把这条链路理解为一个可回溯的整体，比记住七个孤立步骤更有价值。