finish
/clear
/prep7 
!嵥ԪԵ 
rd0=20.0 !ֱֽ 
et,1,solid65 
et,2,link8 
mp,ex,1,33e3 
mp,prxy,1,0.20 
r,1 
hntra=28 
hntrl=2.6 
tb,concr,1 
tbdata,,0.7,1.0,hntrl,-1 
mp,ex,2,2.1e5 
mp,prxy,2,0.30 
r,2,acos(-1)*0.25*rd0*rd0 
!弴Ԫ 
blc4, , ,100,200,3000 
/view,1,1,1,1 
/ang,1 
gplot 
!ʱ߳ 
lsel,s,loc,z,1,2999 
lsel,r,loc,y,0 
latt,2,2,2 
lesize,all,,,20 !ֽĿ 
lmesh,all 
lsel,s,loc,z,0 
lesize,all,,,4 !ϵĿ4x4 
vsel,all 
vatt,1,1,1 
mshape,0,3d 
mshkey,1 
vmesh,all 
allsel,all 
finish 
/solu 
!ʩԼ 
lsel,s,loc,z,0 
lsel,r,loc,y,0 
dl,all,,uy 
dl,all,,uz 
lsel,all 
lsel,s,loc,z,3000 
lsel,r,loc,y,0 
dl,all,,uy 
lsel,all 
ksel,s,loc,x,0 
ksel,r,loc,y,0 
dk,all,ux 
allsel,all 
!ʩӺ 
qmz=0.3 
asel,s,loc,y,200 
sfa,all,1,pres,qmz 
allsel,all 
nsubst,40 
outres,all,all 
time,qmz*10 
neqit,40 
solve 
finish 


/post1 
pldisp,1
plnsol,s,1,0,1
plnsol,s,3,0,1
plnsol,s,yz,0,1
/DEVICE,VECTOR,1 
set,last 
plcrack 
fini
