ale2d方法模拟楔形体入水
Wb前处理，2d分析
1）模型：2个面1个part下+楔形体
2）单元，平面应变对称
3）材料空气水及及钢【刚体,WB下2D不支持刚体】
*MAT_NULL
$水
1,1000
*EOS_LINEAR_POLYNOMIAL
1,101325,2.25E9
0
*MAT_NULL
$空气
2,1.2929
*EOS_LINEAR_POLYNOMIAL
$环境压力不是吧下边的0改为1e5，而是使用control_ale控制；
2,0,,,,0.4,0.4,0
2.5E5,1
4）网格
5）Part【1-水，2-空气，3-楔形体】
6）6m/s速度；
7）输出设置，单位制mkgs【建议用mkgs标准单位制】，时间0.05s下降0.3m左右【模型0.8m】
8）k文件修改：
$
*SECTION_SHELL
1 13 1.0000 1.0 0.0 0.0 0 1
0.00 0.00 0.00 0.00 0.00
$
*SECTION_ALE2D
$流体单元，可用13号算法（13平面应变14轴对称）
2,11,,13
$
*CONTROL_ALE
$ale控制，NE.-1: Use default advection logic.，1为1个循环，EQ.-2 Modified Van Leer ；
-1,1,-2,0
$1.013e5为环境压力， PREF :=Reference pressure applied to the free surfaces of the ALE mesh boundary.
,,,,,0,1.013e5
$
*ALE_MULTI-MATERIAL_GROUP
$流体物质组
1,1
2,1
$
*SET_PART_LIST
1
1,2
*SET_PART_LIST
2
3
*SET_MULTI-MATERIAL_GROUP_LIST
$设置多物质组集 合，对应mcoup的-101
101
$1为水，只与水耦合；
1
*CONSTRAINED_LAGRANGE_IN_SOLID
$# slave master sstyp mstyp nquad ctype direc mcoup
2 1 0 0 2 4 3 -101
$# start end pfac fric frcmin norm normtyp damp
0.0 0.0 0.00 0.0 0.3 0 0 0.0
$# cq hmin hmax ileak pleak lcidpor nvent blockage
0.0 0.0 0.0 0.3 0 0 0
$# iboxid ipenchk intforc ialesof lagmul pfacmm thkf
0 0 0 0 0.0 0 0.0
$*CONSTRAINED_LAGRANGE_IN_SOLID
$ 1 2 0 0 2 5 2 0
$ 0 0 0.15 0.1
$ 2 0.3
$
*initial_VELOCITY_GENERATION
1,2,,0,-0.01
$
$
*DATABASE_TRACER
$TIME=0为粒子开始记录时间；TRACK=0表示粒子跟随物质运动，=1表示粒子在空间固定不动；xyz表示粒子初始坐标（定义在炸药的外径上）；ammgid为物质组号，=0表示ale单元内物质按照体积分数加权；
0,0,0.5,0.5,0,0
*DATABASE_TRACER
$TRACK=1表示粒子在空间固定不动；xyz表示粒子初始坐标=78，看x=78固定位置的压力；
0,1,0.5,0.5,0,0
*DATABASE_TRHIST
1
*DATABASE_FSI
$# dt binary
0.001 1
$#dbfsi_id sid sidtype输出楔形体part3与水101的留固耦合作用力
101 3 1
$
