1

*CONTROL_ENERGY 2 2 2 2 *CONTROL_SHELL 20.0 1 -1 1 2 2 1 *CONTROL_TIMESTEP

0.0000 0.9000 0 0.00 0.00 *CONTROL_TERMINATION

$1000.0000 0 0.00000 0.00000 0.00000 0.800E+05 0 0.00000 0.00000 0.00000 *CONTROL_SPH 2 0 $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ TIME HISTORY $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*DATABASE_BINARY_D3PLOT 200.0E+00 $0.500E+00

*DATABASE_BINARY_D3THDT 0.8000E+02

*DATABASE_EXTENT_BINARY 0 0 3 1 0 0 0 0 0 0 4 0 0 0 $ $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ SECTION DEFINITIONS $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*SECTION_SOLID 2 1

*SECTION_SPH 1 $

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ MATERIAL DEFINITIONS $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*MAT_RIGID

2 7.80 2.10 0.300000 0.0 0.0 0.0 1.00 7.00 7.00

*MAT_RIGID

3 7.80 2.10 0.300000 0.0 0.0 0.0 1.00 6.00 7.00

*MAT_RIGID

4 7.80 2.10 0.300000 0.0 0.0 0.0 1.00 6.00 7.00

*MAT_NULL 1 1.00 *EOS_GRUNEISEN

1 .1484000 1.9790000 .0000000 .0000000 .1100000 3.0000000 .0000000 .0000000 $ $ $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ PARTS DEFINITIONS $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $ *PART

Part 1 for Mat 4 and Elem Type 1 1 1 1 1 0 0 0 $ *PART

Part 2 for Mat 2 and Elem Type 1

2 2 2 0 0 0 0 $ *PART

Part 3 for Mat 3 and Elem Type 1 3 2 3 0 0 0 0 *PART

Part 3 for Mat 3 and Elem Type 1 4 2 4 0 0 0 0 $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ RIGID BOUNDRIES $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*LOAD_BODY_Y 2,1.0

*DEFINE_CURVE 2

0.0,9.8E-10 1.0,9.8E-10 $

*DEFINE_CURVE

1 0 1.000 1.000 0.000 0.000

0.000000000000E+00 1.000000000000E-04 1.000000000000E+05 1.000000000000E-04 *BOUNDARY_PRESCRIBED_MOTION_RIGID 3 2 0 1 -1.00 0 0.000 0.000 $ $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ NODE DEFINITIONS $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*CONTACT_AUTOMATIC_NODES_TO_SURFACE 1 1 3 2 0 0 0 0

0.000 0.000 0.000 0.000 0.000 0 0.000 0.0000E+08 0.000 0.000 0.100 0.000 0.000 0.000 0.000 0.000

1 0.1000000 3 *SET_PART_LIST 1 2,3,4 *INCLUDE mesh.k *END

3、 ALE（接触算法）

采用接触算法分析流固耦合问题也是一种选择，在液面波动幅度较小时可以采用此种方法进行分析，流体用ALE算法描述，结构采用Lagrange算法；需要注意的一点：对ALE网格要进行滑顺处理，以控制网格形态，保证求解精度。 下面是某算例的部分命令流：

*KEYWORD *TITLE ALE $

*DATABASE_FORMAT 0 $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ CONTROL OPTIONS $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*CONTROL_PARALLEL 1

*CONTROL_ENERGY 2 2 2 2 *CONTROL_ALE

3 1 2 1.0000000 1.0000000 0.000000 1.0000000 1.0000e+9 0.000000 0.000000 2 *CONTROL_TIMESTEP

0.0000 0.9000 0 0.00 0.00 *CONTROL_CONTACT

0.0000000 0.0000000 1 0 2 0 0 0 0 0 0

*CONTROL_TERMINATION