Ls-dyna动力分析，DP材料模型

可以参见关键字手册，严格执行就行

修改是右对齐，或者参数之间加逗号

先发现原因是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