可以參見關鍵字手冊，嚴格執行就行

修改是右對齊，或者參數之間加逗號

先發現原因是ansys/ls-dyna中沒有D-P材料模型，需要修改k文件自定義土壤D-P材料模型，但是不知道該材料模型的k文件中參數格式。

@周游

@lc @落岸天使 @crystal_w

已附上命令流，求各位給看看什么原因！

!=============================================================================== ! 強夯沖擊動應力分析 ! m-N-s-Pa-kg-kg/m3 !=============================================================================== !----分析環境設置--------------------------------------------------------------- FINISH /CLEAR /FILNAME,Pounder,0 /TITLE,Analysis Of Dynamic Compaction !----定義幾何、材料參數--------------------------------------------------------- /PREP7 J=1200E3 !夯擊能 H0=10 !夯錘落距 R_HC=1 !夯錘半徑 H_HC=0.5 !夯錘高 S_TR=6 !土壤水平方向尺寸 H_TR=10 !土壤豎向尺寸 DENS_HC=7500 E_HC=2.1E11 u_HC=0.23 DENS_TR=1600 E_TR=3E6 C_TR=21E3 FI_TR=13 u_TR=0.29 !----定義單元類型--------------------------------------------------------------- ET,1,SOLID164 !----定義材料模型--------------------------------------------------------------- EDMP,RIGI,1,7,7 MP,DENS,1,DENS_HC MP,EX,1,E_HC MP,NUXY,1,u_HC MP,DENS,2,DENS_TR MP,EX,2,E_TR MP,NUXY,2,u_TR TB,DP,2 TBDATA,2,C_TR,FI_TR,0 !----建立1/4幾何模型------------------------------------------------------------ BLC4,0.0,0.0,2*S_TR,H_TR,-S_TR WPOFFS,,H_TR WPOFFS,S_TR WPROTA,,-90 CYL4,,,R_HC,180,,,H_HC WPCSYS,-1 /WINDOW,1,LTOP !Window 1 左上角 /REPLOT /WINDOW,2,RTOP !Window 2 右上角 /REPLOT /WINDOW,3,LBOT !Window 3 左下角 /REPLOT /WINDOW,4,RBOT !Window 4 右下角 /REPLOT /VIEW,1,-1 !Window 1 左視圖 /REP,FAST /VIEW,2,,1 !Window 2 俯視圖 /REP,FAST /VIEW,3,,,1 !Window 3 前視圖 /REP,FAST /VIEW,4,1,1,1 !Window 4 Iso視圖 /REP,FAST /WINDOW,2,OFF /WINDOW,3,OFF /WINDOW,4,OFF /WINDOW,1,FULL /VIEW,1,1,1,1 /REPLOT !----網格劃分------------------------------------------------------------------- WPOFFS,S_TR WPROTA,,,90 VSBW,ALL WPOFFS,,,1/3*S_TR VSBW,ALL WPOFFS,,,-2/3*S_TR VSBW,ALL NUMCMP,ALL VSEL,S,LOC,Y,-0.1,H_TR+0.1 VPLOT CM,VTR,VOLU ALLSEL,ALL VSEL,S,LOC,Y,H_TR-0.1,H_TR+H_HC VPLOT CM,VHC,VOLU ALLSEL,ALL LSEL,S,LINE,,13,14,1 LSEL,A,LINE,,17,18,1 LSEL,A,LINE,,22,25,1 LSEL,A,LINE,,27,28 LESIZE,ALL,,,6 LSEL,ALL LSEL,S,LINE,,19,21,1 LSEL,A,LINE,,26 LESIZE,ALL,,,3 LSEL,ALL CMSEL,S,VHC MAT,1 VMESH,VHC VSEL,ALL LSEL,ALL LSEL,S,LENGTH,,S_TR LESIZE,ALL,,,30 LSEL,ALL LSEL,S,LENGTH,,H_TR LESIZE,ALL,,,50 LSEL,S,LENGTH,,1/3*S_TR LESIZE,ALL,,,20 LSEL,S,LENGTH,,2/3*S_TR LESIZE,ALL,,,20 CMSEL,S,VTR MAT,2 VMESH,VTR ALLSEL WPCSYS,-1 EPLOT /WINDOW,1,LTOP !Window 1 左上角 /REPLOT /WINDOW,2,RTOP !Window 2 右上角 /REPLOT /WINDOW,3,LBOT !Window 3 左下角 /REPLOT /WINDOW,4,RBOT !Window 4 右下角 /REPLOT /VIEW,1,-1 !Window 1 左視圖 /REP,FAST /VIEW,2,,1 !Window 2 俯視圖 /REP,FAST /VIEW,3,,,1 !Window 3 前視圖 /REP,FAST /VIEW,4,1,1,1 !Window 4 Iso視圖 /REP,FAST EDPART,CREATE PARTSEL,'PLOT' /WINDOW,2,OFF /WINDOW,3,OFF /WINDOW,4,OFF /WINDOW,1,FULL /VIEW,1,1,1,1 /REPLOT !----定義表面接觸信息----------------------------------------------------------- EDCGEN,ESTS,2,1 EDCONTACT,0.5,,,,,,1,,0 !----施加位移邊界條件----------------------------------------------------------- ASEL,S,LOC,X,-0.01,0.01 ASEL,A,LOC,X,2*S_TR-0.01,2*S_TR+0.01 APLOT NSLA NPLOT D,ALL,UX ASEL,ALL ASEL,S,LOC,Z,-0.01,0.01 ASEL,A,LOC,Z,-S_TR-0.01,-S_TR+0.01 APLOT NSLA NPLOT D,ALL,UZ ASEL,ALL ASEL,S,LOC,Y,-0.01,0.01 APLOT NSLA NPLOT D,ALL,UY ALLSEL,ALL EPLOT !----施加無反射邊界條件----------------------------------------------------------- NSEL,S,LOC,X,-0.01,0.01 NSEL,A,LOC,X,2*S_TR-0.01,2*S_TR+0.01 NSEL,A,LOC,Z,-S_TR-0.01,-S_TR+0.01 NPLOT CM,NONR,NODE EDNB,ADD,NONR,1,1 ALLSEL,ALL !----定義夯錘的初始速度--------------------------------------------------------- EDPVEL,VGEN,1,0,-14,0.000 ACEL,0,-9.8,0 FINISH !----求解選項設置--------------------------------------------------------------- /SOLU !設置能量選項 EDENERGY,1,1,1,1 !設置時間步長因子 EDCTS,0,0.002 !設置求解時間 TIME,0.2 !設置結果文件輸出步數 EDRST,50 EDHTIME,1 !設置結果文件輸出類型 EDOPT,ADD,blank,LSDYNA !輸出關鍵字文件 EDWRITE,LSDYNA,Pounder,K