LS-DYNA深水爆炸水介质参数初始化方法

在借助数值仿真方法研究爆炸载荷随水深变化规律时，需设定流场的初始热力学参数。在《水介质初始参数设置对水下爆炸载荷的影响》一文中（注：郑永辉, 魏继锋. 爆炸与冲击, 2022, 42(05): 63-72），已经给出了水介质初始参数的设置方法，即根据等温过程或者等内能过程假设计算不同水深下水介质的初始密度和内能或者直接使用LS-DYNA提供的关键字INITIAL_EOS_ALE，此处不再赘述。不过，文中只给出了流体计算域中流体静压力处处相同的工况（简化为一维轴对称模型），当考虑水介质参数在水深方向的连续变化时，还需要在上述基础上进行拓展，对于LS-DYNA程序，本文给出了两种设置方法，各有优劣。

一、关键字INITIAL_HYDROSTATIC_ALE和ALE_AMBIENT_HYDROSTATIC

根据相关资料，采用关键字INITIAL_HYDROSTATIC_ALE和ALE_ AMBIENT_HYDROSTATIC来设置水介质参数水深方向的连续变化是一种较新的方法，也是官方较为推荐的方法。关键字手册中对关键字INITIAL_ HYDROSTATIC_ALE的描述如下：“When an ALE model contains one or more regular (not reservoir-type) ALE parts (ELFORM = 11 and AET = 0), this command may be used to initialize the hydrostatic pressure field in the regular ALE domain due to gravity .The *LOAD_-BODY_(OPTION) keyword must be defined.”对关键字ALE_AMBIENT_HYDROSTATIC的描述为“When an ALE model contains one or more ambient (or reservoir-type) ALE parts (ELFORM = 11 and AET = 4), this command may be used to initialize the hydrostatic pressure field in the ambient ALE domain due to gravity. The *LOAD_-BODY_{OPTION} keyword must be defined. The associated *INITIAL_HYDROSTAT-IC_ALE keyword may be used to define a similar initial hydrostatic pressure field for the regular ALE domain (not reservoir-type region).”从中可以看出，两个关键字分别适用于不同类型的Part，其中，第二个关键字的本质作用在于施加边界条件以维持边界处的压力梯度，使其不随时间变化；此外，在定义这两个关键字的同时，也必须同时定义加载体力的关键字，本质上是在流场施加重力加速度。关于关键字的具体使用，可参考关键字手册中的相关说明及案例。
虽然该组关键字较为简单，但也存在一定的限制，包括：（1）关键字对水介质初始参数的求解算法未知，即水介质参数水深的变化规律未知，只能根据计算结果去反推；（2）对水介质状态方程有限制，R11版本的手册指出，两个关键字暂时只支持EOS_GRUNEISEN和EOS_LINEAR_POLYNOMIAL形式的状态方程，而且对于EOS_LINEAR_POLYNOMIAL，进一步限制如下
   c 3 = c 4 = c 5 = c 6 = 0 E 0 = 0 c 4 = c 5 > 0 , c 1