LSDYNA小球入水分析
参考文献 Modelling strategies for numerical simulation of aircraft ditching
International Journal of Crashworthiness 4区 影响因子2.055
FEM→ELMENT TOOLS→Duplicate Nodes→选取上下的重合节点，将其合并
FEM→MEASURE→选取球，测量其体积
质量3.97，除以球的体积，得到四分之一球的密度
采用Element Editing→Area→Delete删去多余的四分之三球
材料量纲为mm-kg-ms-kN-GPa 转化为m-kg-s 1000kg/m3=1000*kg/109mm3=10-6kg/mm3
定义真空Vacuum，材料编号1，密度1.0e-12
定义Water，材料编号2，密度1.0e-6，动力黏度系数1e-8
定义刚体Rigid sphere，编号3，密度7.02e-7，弹性模量E=200GPa，泊松比0.3，约束X和Y的平动自由度和所有转动自由度
Define SECTION SOLID
MM ALE 编号1 ELFORM 11
LAGRANGE 编号3 ELFORM 1
定义状态方程EOS_LINEAR_POLYNOMIAL
EOSID=2，对应Water的编号，C1=2.723，C2=7.727，C3=14.66
至此材料定义完毕
定义沙漏HOURGLASS
HGID=1，JHQ=1，QM=1e-6
定义PART
此处注意Vacuum和Water的SECID=1，里面设置了ELFORM=11的算法，且Hourglass=1，里面设置了QM=1e-6，而Rigid Sphere的SECID=3，里面设置了ELFORM=1的算法，且Hourglass=0
Part1 Vacuum，SECID=1，Hourglass=1
Part2 Water，SECID=1，EOSID=2，Hourglass=1
Part3 Rigid Sphere SECID=3 Hourglass=0
定义part集SET_PART_LIST
集合编号1，里面包含part1和part2
定义CONSTRAINED_LAGRANGE_IN_SOLID
SLAVE=3 MASTER=1 SSTYP=1 MSTYP=0 意味着SLAVE是part ID 3 MASTER是part set ID 1，NQUAD=3，CTYPE=4，DIREC=2 意味着如果球落水则水会被压缩到上面，而如果球提升之后，水不会被球所拉升，MCOUP=1，只与有最高密度的材料相耦合
定义ALE_MULTI_MATERIAL_GROUP
Vacuum SID=1 IDTYPE=1
Water SID=2 IDTYPE=1
定义边界条件
Create Entity→Cre→Boundary→Spc→By element 勾选prop并选择水的底面，约束所有方向的平动和转动→Apply
Create Entity→Cre→Boundary→Spc→By element 勾选prop并选择对称面1，限制Y向平动和X和Z向的转动
Create Entity→Cre→Boundary→Spc→By element 勾选prop并选择对称面2，限制X向平动和Y和Z向的转动
定义CONTROL_ALE
DCT=-1，NADV=1，METH=2，AFAC=-1
其他采用默认值
定义CONTROL_TERMINATION
定义ENDTIM=10 结束时间10ms
定义CONTROL_TIMESTEP
TSSFAC=0.67 计算时间步的缩放因子，默认0.9，如果使用了高爆物，默认值降为0.67
定义DATABASE_ASCII_option
Default DT=0.001 勾选RBDOUT DT=0.001
定义DATABASE_BINARY_FSIFOR
DT=0.001 PSETID=1 选取输出FSI力，选择的Part set ID为1
定义DATABASE_BINARY_D3PLOT
NPLTC=100 控制输出步数为100
定义DATABASE_FSI
Model Checking→Keyword Check
定义INITIAL_VELOCITY_GENERATION
Assembly and Select Part
由四分之一模型转换成全部模型，Settings→Post Settings→Reflections
最后生成k文件计算时，如果提示CONTROL_ALE读取错误，将最后一行删除之后可以正常计算
*CONTROL_ALE
-1 1 2 -1.0 0.0 0.0 0.0 0.0
0.01.00000E20 1.01.00000E-6 0 0 0.0 0
1 50 0 0.0 0.0 0 0.0 0.0
0 0 01.00000E-5
