别让参数选错毁了SPH模拟！LS-DYNA三维建模要点详解（附赠模型）

详细信息

作者优秀: 优秀教师/意见领袖/博士学历/特邀专家/独家讲师
平台推荐: 内容稀缺

8月前浏览810

平滑粒子流体动力学（Smoothed Particle Hydrodynamics，SPH）方法因其高度的灵活性和对复杂几何与大变形问题的良好适应能力，在工程力学、冲击动力学以及流固耦合等领域得到了广泛应用。

作为一种典型的无网格方法，SPH 从根本上避免了传统基于网格方法在大变形条件下面临的网格畸变问题，从而在一定程度上提升了计算的稳定性与适用范围。然而，SPH方法在实际使用过程中对建模策略和参数设置具有较高的敏感性，若相关参数选取不当，极易引发数值噪声、粒子非物理运动甚至计算发散。因此，在开展 SPH 数值模拟之前，合理的建模方法与参数理解尤为关键。

本文将重点围绕 LS-DYNA 中 SPH 三维模型的建模流程展开介绍，并结合工程实践经验，对建模过程中常见的问题与注意事项进行说明。

LS-DYNA中SPH粒子的生成方式

在 LS-PrePost 中，SPH 粒子的建立主要通过 SPH Generation 功能完成。该功能提供了多种粒子生成方式，主要包括：

基于标准几何体形状生成粒子；
在由壳单元或梁单元所围成的空间内生成粒子；
基于实体单元的节点或单元中心生成粒子。

在 SPH 建模过程中，必须遵循的一个基本原则是粒子分布的协调性。即粒子在空间中的分布应尽可能均匀、连续，避免出现局部过度加密、稀疏不均或杂乱无章的情况。粒子分布不均往往会引入较大的数值噪声，显著降低计算稳定性，甚至直接导致计算失败。基于上述原因，并不推荐采用实体单元节点生成粒子的方式，尤其是基于四面体或四边形实体单元节点生成粒子。该方式极易造成粒子间距不一致，从而破坏 SPH 方法的基本假设。此外，虽然可以通过实体单元中心点生成 SPH 粒子，但该方法在实际应用中灵活性较差，不利于后续统一调整粒子分辨率。相比之下，直接通过指定空间范围生成粒子或基于壳单元空间生成粒子在工程应用中更为高效和可靠。

复杂几何模型的粒子生成策略

对于几何形状较为复杂的模型，尤其是存在明显起伏、非规则边界或微观结构特征的情况，推荐采用基于壳单元所围空间生成 SPH 粒子的方法。

例如，如图所示的微观研磨垫表面模型，其表面形貌复杂，若采用规则几何体直接生成粒子，将难以准确表征真实几何边界。因此，该类模型通常通过壳单元定义计算域，并在其围成的空间内生成 SPH 粒子。该方法不仅能够较好地保持复杂几何的完整性，同时也有助于保证粒子分布的均匀性，从而提高后续计算的稳定性与结果可靠性。

粒子间距的定义与常见误区

在确定 SPH 粒子的生成方式之后，需要进一步设置粒子之间的间距，该参数通过 Num Particles Definition 选项进行定义。需要特别注意的是，该选项在不同粒子生成方式下，其物理含义并不完全一致。

当采用标准几何体形状生成粒子时，该选项表示在三个坐标方向上生成的粒子数量；

当采用基于壳单元空间生成粒子时，该选项则表示三个方向上的粒子间距。

这一差异是 LS-DYNA 使用过程中较为典型、且极易引发误操作的问题之一。相同的界面选项在不同功能模式下对应不同的定义方式，若未加以区分，极易导致粒子尺度设置错误，从而影响计算精度甚至引发数值不稳定。因此，在 SPH 建模过程中，务必结合所选粒子生成策略，明确相关参数的实际物理含义。

SPH 粒子质量的设置与校核

在 SPH 建模过程中，除粒子的几何参数外，还必须为粒子赋予合理的质量属性。粒子质量通常通过材料密度与几何体总体积相乘得到系统总质量，再将该质量均匀分配至所有 SPH 粒子的方式进行确定。

粒子质量参数在数值计算中具有至关重要的作用，其大小将直接影响粒子的动量、惯性响应以及体系的整体动力学行为。质量设置不当往往会导致计算结果出现明显偏差，甚至引发数值不稳定问题。

在实际建模中，这里提供一个较为实用的技巧：在定义 SPH 粒子属性时，可将 Density 参数设置为 -1 。此时程序将自动调用所关联材料模型中的密度值。采用该方式后，即使后续对材料参数进行修改，也无需重新生成 SPH 粒子模型，从而显著提高建模效率并减少人为错误。

在完成 SPH 粒子建模后，可通过 SPH Generation 中的 Modify 功能对模型质量进行检查。该功能本质上用于修改粒子质量，但在实际工程应用中更适合作为一种校核工具，用于确认当前 SPH 模型的总质量及单粒子质量是否合理。

LS-DYNA后续内容安排

本文主要围绕 SPH 三维模型的建模方法进行了系统介绍。对于二维平面 SPH 模型，其基本建模思想与参数设置原则与三维情况基本一致，因此在此不再赘述。在后续文章中，将进一步围绕 SPH 模型中的边界条件处理方法以及关键控制参数的设置技巧展开讨论，重点分析不同参数对计算稳定性与结果精度的影响，并结合工程算例给出相应的设置建议。

温馨提示：我把本文相关的模型资料放到仿真学习资料包，大家可以在资料包搜索文章标题的文件夹下载。

为了帮助读者理解并掌握SPH分析方法理论及应用，笔者在仿真秀官网独家首发了《 LS-DYNA ALE流固耦合专题18讲：鸟撞/船舶/爆炸/激光融化/射流仿真》视频课程，通过课程的学习可快速掌握LS_DYNA中SPH法的使用及相关技巧，并且通过案例的演示帮助理解各个关键字以及相关参数的设置方式及含义。欢迎识别下方二维码试看。

可回放，开发票，奖学金，加餐

提供vip交流群/答疑/学习资料

《 LS-DYNA ALE流固耦合专题18讲：鸟撞/船舶/爆炸/激光融化/射流仿真》

课程配套资料如下，提供VIP群答疑服务。

1、LS-Dyna2023R1基础入门

该视频教程基于ANSYS LS-dyna最新的版本2023 R1进行讲解，手把手全流程界面操作演示。具有完整的案例介绍，从建模、设置边界条件、求解设置、后处理全流程讲解；包括19个最新仿真设计案例讲解与下载，涵盖包括航空、汽车等多个工业制造行业，及多个常见的模拟应用场景演示，工程实践性强；适合学生、老师以及有意愿学习显式求解算法的初学者。