LSDYNA中ALE与SALE的新旧关键字

本文链接：https://blog.csdn.net/xiaocai13/article/details/142975130

1.流固耦合：

旧：*CONSTRAINED_LARGRANGE_IN_SOLID

新：*ALE_STRUCTURED_FSI

• 耦合类型： *ALE_STRUCTURED_FSI 只有罚函数公式耦合，类似于 *CONSTRAINED_LAGRANGE_IN_SOLID 中 CTYPE=4 和 5 的组合

• 耦合点数量。对于每个拉格朗日段，都有一定数量的耦合点，均匀分布在段面上。这些耦合点连接着惩罚弹簧。在使用* CONSTRAINED_LAGRANGE_IN_SOLID 时，需要指定数量，而使用此关键字时，则无需指定，因为 LS-DYNA 会在初始化阶段自动确定耦合点数量。

• 泄漏控制。使用 *ALE_STRUCTURED_FSI 时，泄漏控制是自动的。流体泄漏被自动检测并修复，无需用户预先设置。

• 正态类型。正态类型选择是基于本地几何图形自动进行的。用户不需要在节点 / 段正态之间进行选择。

• 边缘耦合。边缘耦合是自动的。壳单元段中暴露的边缘会耦合。不需要使用* CONSTRAINED LAGRANGE IN SOLID_-EDGE 。

• 侵蚀耦合。侵蚀的固体元素将改变耦合段。属于侵蚀元素的段需要删除新暴露的段需要添加。使用* CONSTRAINED LAGRANGEIN SOLID 时，需要指定 CTYPE=5 才能激活此段修改算法。在新的* ALE_STRUCTURED_FSI 中，该过程始终处于打开状态，不需要任何标志。

在SALE中官方更推荐*ALE_STRUCTURED_FSI，新的关键词对泄露控制的更好。两种关键词对计算结果影响较大。

2.体积填充：

旧：*INITIAL_VOLUME_FRACTION_GEOMETRY

新：*ALE_STRUCTURED_MESH_VOLUME_FILLING

经对比，两者不同对计算结果有较小区别

官方手册中认为*INITIAL_VOLUME_FRACTION_GEOMETRY过于老旧，很少使用。

3.多物质：

旧：*ALE_MUTLI-MATERIAL_GROUP

新：*ALE_STRUCTURED_MUTLI-MATERIAL_GROUP

*ALE_STRUCTURED_MUTLI-MATERIAL_GROUP是为SALE新增的关键字，好像是ANSYS 2024R2中的LSPP才有，新的关键字可以添加材料和状态方程，以及基准压力，取代了*CONTROL_ALE中的PREF。

PLNEPS和AXISYM可以模拟2D平面和轴对称的情况

对计算结果的影响未验证

4.多物质组集

旧：*SET_MULTI

新：*SET_MULTI_MATERIAL_GROUP_LIST、*SET_MULTI_MATERIAL_GROUP_LIST_GPNAME

两者功能几乎相同，官方手册中提到*SET_MULTI是旧关键字，两者使用后对计算结果有较大区别。在最新版2024R2中，*SET_MULTI_MATERIAL_GROUP_LIST又消失了，只剩*SET_MULTI_MATERIAL_GROUP_LIST_GPNAME，暂未查到两者的区别。

2024.10.17：经研究，*SET_MULTI是新版本中把以前的*SET_MULTI_MATERIAL_GROUP_LIST名字缩减了。

5.*CONTROL_ALE

参考B站教程及其他相关的教程，爆炸流固耦合仿真常使用以下参数

DTC-调用备用对流逻辑，默认设置为“-1”，更好模拟炸药；

NADV-循环次数，默认设置为“1”

METH-(欧拉网格)移流算法（聚能射流一般选择METH=2 Van Leer with HIS)

METH=-2与METH=2的区别：

METH=2时，单元前后上下左右六个方向上的连续3个单元（共18个单元）的变形程度被限制

METH=-2时，变形程度限制在周围26个单元上（不考虑方向）。能更好地保留*MAT_HIGH_EXPLOSIVE_BURN材料界面

METH=1与METH=3的区别：

METH=1，每个网格移流步骤仅保留单元内能

METH=3，每个网格移流步骤保留全部能量

但在LSDYNA发布的Structured ALE Tutorial教程（https://ftp.lstc.com/anonymous/outgoing/hao/sale/tutorials）中，

在SALE中，官方希望我们尽可能极简*CONTROL_ALE的卡片

对比使用后发现两者的区别对计算结果影响较大。

6.流固耦合中的MCOUP

在*CONSTRAINED_LARGRANGE_IN_SOLID和*ALE_STRUCTURED_FSI中，都有MCOUP。对应*SET_MULTI或者*SET_MULTI_MATERIAL_GROUP。如果填0，是所有物质都考虑，但是官方手册中提到，如果填0耦合只考虑50%。所以应该指定材料，但也有两种情况：