Ansys/ls-dyna 顯示動(dòng)力分析，施加初始速度，剛體未產(chǎn)生運(yùn)動(dòng)。
!===============================================================================
! 強(qiáng)夯沖擊動(dòng)應(yīng)力分析
! m-N-s-Pa-kg-kg/m3
!===============================================================================
!----分析環(huán)境設(shè)置---------------------------------------------------------------
FINISH
/CLEAR
/FILNAME,Pounder,0
/TITLE,Analysis Of Dynamic Compaction
/UNITS,SI
!----定義幾何、材料參數(shù)---------------------------------------------------------
/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
!----網(wǎng)格劃分-------------------------------------------------------------------
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
!----施加無(wú)反射邊界條件-----------------------------------------------------------
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
!----求解選項(xiàng)設(shè)置---------------------------------------------------------------
/SOLU
!設(shè)置能量選項(xiàng)
EDENERGY,1,1,1,1
!設(shè)置時(shí)間步長(zhǎng)因子
EDCTS,0,0.4
!設(shè)置求解時(shí)間
TIME,0.2
!設(shè)置結(jié)果文件輸出步數(shù)
EDRST,100
EDHTIME,1
!設(shè)置結(jié)果文件輸出類型
EDOPT,ADD,blank,LSDYNA
!輸出關(guān)鍵字文件
EDWRITE,LSDYNA,Pounder,K
剛體不發(fā)生運(yùn)動(dòng)，是怎么回事
藍(lán)牙 
