FINISH$/CLEAR$/CONFIG,NRES,2000     !Ӳֵ
/prep7   
!/nerr,nmerr,nmabt                  !رվϢ
et,1,plane42
et,2,solid45
et,3,combin14
et,4,surf154

!һŲ
mp,ex,1,1.2e9
MP,PRXY,1,0.35
mp,dens,1,2613
!Ų
mp,ex,2,1e9
MP,PRXY,2,0.35
mp,dens,2,2000


!ģ
k,1,-0.14,0
k,2,0.14,0
k,3,0.14,0.6
k,4,-0.14,0.6
l,1,2
l,2,3
l,3,4
l,4,1
al,1,2,3,4
rectng,0.14,3,0,0.6
rectng,-0.14,-3,0,0.6
aplot
nummrg,all               
numcmp,all
aplot

!1
lsel,s,,,1,3,2
lesize,all,,,4
lsel,s,,,2,4,2
lesize,all,,,6
amesh,1
eplot

!2

lsel,s,,,6
lesize,all,,,6
lsel,s,,,5
lesize,5,,,8,4
lsel,s,,,7
lesize,7,,,8,0.25                 
amesh,2
eplot

!3
lsel,s,,,8
lesize,8,,,8,0.25
lsel,s,,,9
lesize,9,,,8,4  
lsel,s,,,10
lesize,all,,,6
amesh,3
eplot
!
extopt,esize,200,0
type,2
mat,1
vext,1,3,1,,,10
/view,1,1,1,1
vplot
eplot
!ı
esel,s,mat,,1
nsle,s
nplot
nsel,r,loc,y,0.401,0.6
nsel,r,loc,x,-3,3
nplot
esln,s
mpchg,2,all
/replot
allsel
eplot

allsel                           !ϲѹ   
nummrg,all
numcmp,all
/pnum,mat,1                     !ֲͬɫʾ  
eplot

!ģ,1
allsel
csys,0
ngen,2,500000,1,5,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,1,5,1
  e,i,500000+i
*enddo
gplot

!2
allsel
csys,0
ngen,2,500006,36,43,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,36,43,1
  e,i,500006+i
*enddo
gplot

!3
allsel
csys,0
ngen,2,500014,92,99,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,92,99,1
  e,i,500014+i
*enddo
gplot

!4
allsel
csys,0
ngen,2,500021,148,152,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,148,152,1
  e,i,500021+i
*enddo
gplot


!5
allsel
csys,0
ngen,2,500026,183,580,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,183,580,1
  e,i,500026+i
*enddo
gplot

!6
allsel
csys,0
ngen,2,500435,979,1575,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,979,1575,1
  e,i,500435+i
*enddo
gplot


!7
allsel
csys,0
ngen,2,501433,7148,7155,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,7148,7155,1
  e,i,501433+i
*enddo
gplot



!8
allsel
csys,0
ngen,2,501442,7204,7402,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,7204,7402,1
  e,i,501442+i
*enddo
gplot

!9
allsel
csys,0
ngen,2,501542,7602,8994,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,7602,8994,1
  e,i,501542+i
*enddo
gplot

!10
allsel
csys,0
ngen,2,503236,18348,18355,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,18348,18355,1
  e,i,503236+i
*enddo
gplot

!11
allsel
csys,0
ngen,2,503245,18404,18602,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,18404,18602,1
  e,i,503245+i
*enddo
gplot

!12
allsel
csys,0
ngen,2,503600,18802,20194,1,,-2
type,3
r,3,3e7,1e5,
real,3
*do,i,18802,20194,1
  e,i,503600+i
*enddo
gplot







/solu
csys,0
dsys,0
nsel,s,loc,x,-3                 !ʩӱ߽ 
nsel,a,loc,x,3
nplot
d,all,ux
allsel
nsel,s,loc,z,10                !ʩӱ߽ 
nsel,a,loc,z,0 
nplot
d,all,UZ
allsel
nsel,s,loc,y,-2               !ʩӱ߽ 
nplot
d,all,all	
allsel
gplot

alphad,2.6907                   !ϵ 
betad,0.0009                    !նϵ  

anty,trans                  !ѡȡ˲̬ѧ
trnopt,full                !ѡȡȫ
eqslv iterative            !ѡԶ

allsel,all

*do,ii,1,20,1               !м  
outres,all,none  
outres,strs,,mycare  
outres,epel,,mycare  
outres,nsol  
outres,a  
outres,v         
time,ii*0.1/25  

auto,on  
sfedele,all,all,all      
nsel,s,loc,x,-0.14-0.005,0.14+0.005  
nsel,r,loc,y,0.6
nsel,r,loc,z,(90+ii-1)*0.05+0.005,(92+ii+1)*0.05-0.005
esln,s                          !ѡЩԪصĽڵ
/psf,pres,norm,2,0,1            !úʾʽΪͷ
sfe,all,4,pres,,700000           !   0.7MPa

ALLSEL
SOLVE
*ENDDO












