Ansys/ls-dyna 显示动力分析，施加初始速度，刚体未产生运动。
!===============================================================================
! 强夯冲击动应力分析
! m-N-s-Pa-kg-kg/m3
!===============================================================================
!----分析环境设置---------------------------------------------------------------
FINISH
/CLEAR
/FILNAME,Pounder,0
/TITLE,Analysis Of Dynamic Compaction
/UNITS,SI
!----定义几何、材料参数---------------------------------------------------------
/PREP7
J=1200E3 !夯击能
H0=10 !夯锤落距
R_HC=1 !夯锤半径
H_HC=0.5 !夯锤高
S_TR=6 !土壤水平方向尺寸
H_TR=10 !土壤竖向尺寸
DENS_HC=6370
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,6370
MP,EX,1,2.1E11
MP,PRXY,1,0.23
EDMP,HGLS,2,4,0.1,1.5,0.06
MP,DENS,2,1600
MP,EX,2,3E6
MP,PRXY,2,0.29
!----建立1/4几何模型------------------------------------------------------------
BLC4,0.0,0.0,S_TR,H_TR,-S_TR
WPOFFS,,H_TR+0.1
WPROTA,,-90
CYL4,,,R_HC,90,,,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
!----网格划分-------------------------------------------------------------------
WPROTA,,,90
WPOFFS,,,1/3*S_TR
VSBW,ALL
WPOFFS,1/3*S_TR
WPROTA,,,90
VSBW,ALL
WPOFFS,,4/5*H_TR
WPROTA,,90
VSBW,ALL
NUMCMP,ALL
VSEL,S,LOC,Y,-0.1,H_TR
VPLOT
CM,VTR,VOLU
ALLSEL,ALL
VSEL,S,LOC,Y,H_TR+0.001,H_TR+H_HC
VPLOT
CM,VHC,VOLU
ALLSEL,ALL
LSEL,S,LINE,,8,13,1
LESIZE,ALL,,,6
LSEL,ALL
LSEL,S,LINE,,14,16,1
LESIZE,ALL,,,3
LSEL,ALL
CMSEL,S,VHC
TYPE,1
MAT,1
VMESH,VHC
VSEL,ALL
LSEL,ALL
LSEL,S,LENGTH,,1/5*H_TR
LESIZE,ALL,,,20
LSEL,S,LENGTH,,4/5*H_TR
LESIZE,ALL,,,40
LSEL,S,LENGTH,,1/3*S_TR
LESIZE,ALL,,,20
LSEL,S,LENGTH,,2/3*S_TR
LESIZE,ALL,,,20
CMSEL,S,VTR
TYPE,1
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
EDPART,LIST
PARTSEL,'PLOT'
/WINDOW,2,OFF
/WINDOW,3,OFF
/WINDOW,4,OFF
/WINDOW,1,FULL
/VIEW,1,1,1,1
/REPLOT
!----定义表面接触信息-----------------------------------------------------------
EDCGEN,ESTS,2,1,0.15,0.15
!----施加位移边界条件-----------------------------------------------------------
CMSEL,S,VTR
ASLV
ASEL,R,LOC,Z,-0.01,0.01
APLOT
NSLA
NPLOT
D,ALL,UZ
CMSEL,S,VTR
ASLV
ASEL,R,LOC,X,-0.01,0.01
APLOT
NSLA
NPLOT
D,ALL,UX
ASEL,ALL
ASEL,S,LOC,Y,-0.01,0.01
APLOT
NSLA
NPLOT
D,ALL,UY
ALLSEL,ALL
EPLOT
!----施加无反射边界条件-----------------------------------------------------------
NSEL,S,LOC,X,S_TR-0.001,S_TR+0.001
NSEL,A,LOC,Z,-S_TR-0.001,-S_TR+0.001
NPLOT
CM,NONR,NODE
EDNB,ADD,NONR,1,1
ALLSEL,ALL
EPLOT
!----定义夯锤的初始速度---------------------------------------------------------
!CMSEL,S,VHC
!NSLV,S
!CM,VNODE,NODE
!EDVEL,VELO,VNODE,0,-19.796,0
EDPVEL,VGEN,1,0,-19.796,0.000
ALLSEL,ALL
EPLOT
FINISH
!----求解选项设置---------------------------------------------------------------
/SOLU
!设置能量选项
EDENERGY,1,1,1,1
!设置时间步长因子
EDCTS,0,0.4
!设置求解时间
TIME,0.2
!设置结果文件输出步数
EDRST,100
EDHTIME,1
!设置结果文件输出类型
EDOPT,ADD,blank,LSDYNA
!输出关键字文件
EDWRITE,LSDYNA,Pounder,K
刚体不发生运动，是怎么回事
蓝牙 
