!===============================================================================

!                        强夯冲击动应力分析

!                         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

刚体不发生运动，是怎么回事