參考文獻 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