finish
/clear
/prep7
NUMSTR,KP,1
NUMSTR,LINE,1
NUMSTR,AREA,1
NUMSTR,VOLU,1
CSYS,0
k, 1  ,84.375,67.956,.000,
k, 2  ,-202.715,36.029,.000,
l, 1  ,2
k, 3  ,400.507,301.831,.000,
k, 4  ,439.590,238.471,.000,
l, 3  ,4
k, 5  ,704.649,237.944,.000,
l, 4  ,5
k, 6  ,704.649,-104.637,.000,
l, 5  ,6
k, 7  ,-202.715,-104.637,.000,
l, 6  ,7
l, 7  ,2
k, 8  ,272.611,301.831,.000,
l, 3  ,8
k, 9  ,268.768,287.044,.000,
k, 10  ,251.496,229.523,.000,
k, 11  ,194.633,138.969,.000,
k, 12  ,120.527,91.240,.000,
FLST,3,6,3
FITEM,3,8
FITEM,3,9
FITEM,3,10
FITEM,3,11
FITEM,3,12
FITEM,3,1
BSPLIN, ,P51X
k, 13  ,219.517,301.831,.000,
l, 8  ,13
k, 14  ,704.649,67.956,.000,
l, 1  ,14
k, 15  ,150.217,181.901,.000,
l, 15  ,4
k, 16  ,161.517,224.791,.000,
FLST,3,4,3
FITEM,3,13
FITEM,3,16
FITEM,3,15
FITEM,3,1
BSPLIN, ,P51X
alls
lovlap,all
lglue,all
alls
nummrg,all,,,,low
numcmp,all
alls
 


FLST,2,4,3  
FITEM,2,13  
FITEM,2,8   
FITEM,2,17  
FITEM,2,15  
A,P51X  
FLST,2,3,3  
FITEM,2,15  
FITEM,2,17  
FITEM,2,1   
A,P51X  
FLST,2,4,3  
FITEM,2,8   
FITEM,2,3   
FITEM,2,4   
FITEM,2,17  
A,P51X  
FLST,2,5,3  
FITEM,2,17  
FITEM,2,4   
FITEM,2,5   
FITEM,2,14  
FITEM,2,1   
A,P51X  
FLST,2,5,3  
FITEM,2,2   
FITEM,2,1   
FITEM,2,14  
FITEM,2,6   
FITEM,2,7   
A,P51X 
alls
lglue,all
aglue,all
nummrg,all
numcmp,all

mp,ex,1,20e9     !!
mp,prxy,1,0.3
mp,dens,1,2100 
tb,dp,1
tbmodif,1,1,1e6
tbmodif,1,2,50

mp,ex,2,49e7     !!
mp,prxy,2,0.3
mp,dens,2,2100
tb,dp,2
tbmodif,1,1,65e3
tbmodif,1,2,48

mp,ex,3,48e7     !!ͶΪ15%
mp,prxy,3,0.3
mp,dens,3,2140
tb,dp,3
tbmodif,1,1,50e3
tbmodif,1,2,43.2

mp,ex,4,47e7     !!ͶΪ30%
mp,prxy,4,0.3
mp,dens,4,2180
tb,dp,4
tbmodif,1,1,42.5e3
tbmodif,1,2,39.6

mp,ex,5,46e7     !!ͶΪ45%
mp,prxy,5,0.3
mp,dens,5,2210
tb,dp,5
tbmodif,1,1,38.6e3
tbmodif,1,2,38.4

mp,ex,6,45e7     !!ͶΪ60%
mp,prxy,6,0.3
mp,dens,6,2240
tb,dp,6
tbmodif,1,1,37.4e3
tbmodif,1,2,37.2

mp,ex,7,44e7     !!ͶΪ80%
mp,prxy,7,0.3
mp,dens,7,2270
tb,dp,7
tbmodif,1,1,35e3
tbmodif,1,2,36

mp,ex,8,43e7     !!ͶΪ100%
mp,prxy,8,0.3
mp,dens,8,2300
tb,dp,8
tbmodif,1,1,33.8e3
tbmodif,1,2,34.8

et,1,plane42
keyopt,1,3,2

alls
lsel,s,,,2
lsel,a,,,6,7
lsel,a,,,9,11
lsel,a,,,14,16
lesize,all,,,10

lsel,s,,,1,3,2
lsel,a,,,12
lesize,all,,,10,2

lsel,s,,,5,13,8
lesize,all,,,8,0.5

lsel,s,,,8
lesize,all,,,20
lsel,s,,,4
lesize,all,,,30

type,1
mat,2
asel,s,,,2
amesh,all
alls
AMAP,1,13,8,17,15 

alls
type,1
mat,1
 
AMAP,3,8,3,4,17 
AMAP,4,17,5,14,1
AMAP,5,2,14,6,7 

nsel,s,loc,x,704.6490
d,all,ux,0
nsel,s,loc,x,-202.7150
d,all,ux,0
nsel,s,loc,y,-104.6370
d,all,uy,0
allsel
acel,,9.8
finish

/solu
antype,static
alls
fcum,add,, 
deltim,0.2,0.05,0.4 
autots,on !ʹԶʱ䲽 
pred,on !ʱ䲽Ԥ 
!lnsrch,on ! 
nlgeom,on !򿪴λЧ 
nropt,full,,on !趨ţ-ɭѡ 
KBC,0               !Ӻ
TIME,1 
OUTRES,basic,1  !ѡƣÿһӲrescontrol,,all,1,100
cnvtol,f,,0.02,2,1
solve
save
rescontrol,file_summany

/solu
antype,,rest,1,4
asel,s,,,1,2
esla,s
mpchg,3,all
allsel,all
TIME,2
deltim,0.2,0.05,0.4  
autots,on !ʹԶʱ䲽 
pred,on !ʱ䲽Ԥ 
lnsrch,on ! 
nlgeom,on !򿪴λЧ 
nropt,full !趨ţ-ɭѡ 
cutcontrol,plslimit,0.20
KBC,0               !Ӻ
cnvtol,f,,0.02,2,1 
!cnvtol,f,,0.7,2
!NSUBST,30
OUTRES,basic,1   
solve 
save
rescontrol,file_summany

/solu
antype,,rest,2,4
asel,s,,,1,2
esla,s
mpchg,4,all
allsel,all
TIME,3
deltim,0.2,0.05,0.4  
autots,on !ʹԶʱ䲽 
pred,on !ʱ䲽Ԥ 
lnsrch,on ! 
nlgeom,on !򿪴λЧ 
nropt,full !趨ţ-ɭѡ 
cutcontrol,plslimit,0.20
KBC,0               !Ӻ
cnvtol,f,,0.02,2,1 
!cnvtol,f,,0.7,2
!NSUBST,30
OUTRES,basic,1   
solve 
save
rescontrol,file_summany

/solu
antype,,rest,3,4
asel,s,,,1,2
esla,s
mpchg,5,all
allsel,all
TIME,4
deltim,0.2,0.05,0.4  
autots,on !ʹԶʱ䲽 
pred,on !ʱ䲽Ԥ 
lnsrch,on ! 
nlgeom,on !򿪴λЧ 
nropt,full !趨ţ-ɭѡ 
cutcontrol,plslimit,0.20
KBC,0               !Ӻ
cnvtol,f,,0.02,2,1 
!cnvtol,f,,0.7,2
!NSUBST,30
OUTRES,basic,1   
solve 
save
rescontrol,file_summany


/solu
antype,,rest,4,4
asel,s,,,1,2
esla,s
mpchg,6,all
allsel,all
TIME,5
deltim,0.2,0.05,0.4  
autots,on !ʹԶʱ䲽 
pred,on !ʱ䲽Ԥ 
lnsrch,on ! 
nlgeom,on !򿪴λЧ 
nropt,full !趨ţ-ɭѡ 
cutcontrol,plslimit,0.20
KBC,0               !Ӻ
cnvtol,f,,0.02,2,1 
!cnvtol,f,,0.7,2
!NSUBST,30
OUTRES,basic,1   
solve 
save
rescontrol,file_summany


/solu
antype,,rest,5,4
asel,s,,,1,2
esla,s
mpchg,7,all
allsel,all
TIME,6
deltim,0.2,0.05,0.4  
autots,on !ʹԶʱ䲽 
pred,on !ʱ䲽Ԥ 
lnsrch,on ! 
nlgeom,on !򿪴λЧ 
nropt,full !趨ţ-ɭѡ 
cutcontrol,plslimit,0.20
KBC,0               !Ӻ
cnvtol,f,,0.02,2,1 
!cnvtol,f,,0.7,2
!NSUBST,30
OUTRES,basic,1   
solve 
save
rescontrol,file_summany


/solu
antype,,rest,6,4
asel,s,,,1,2
esla,s
mpchg,8,all
allsel,all
TIME,7
deltim,0.2,0.05,0.4  
autots,on !ʹԶʱ䲽 
pred,on !ʱ䲽Ԥ 
lnsrch,on ! 
nlgeom,on !򿪴λЧ 
nropt,full !趨ţ-ɭѡ 
cutcontrol,plslimit,0.20
KBC,0               !Ӻ
cnvtol,f,,0.02,2,1 
!cnvtol,f,,0.7,2
!NSUBST,30
OUTRES,basic,1   
solve 
save
rescontrol,file_summany









9.5  
9.5.1  ʼ״̬±ȶ
ɺ󣬶һڸñµӦλƼӦͼӦͼڽ/POST1еģʶ֮ǰͨúЩͼΪ
/POST1
SET,1
PLNSOL,S,Y,0,1          !ʼ״̬µӦ
PLNSOL,S,XY,0,1         !ʼ״̬µļӦ
PLNSOL,EPTO,XY,0,1     !ʼ״̬µļӦ
    plnsol,eppl,EQV,0,1
    plnsol,s,eqv,0,1  !ЧӦ

/POST1
SET,2
    plnsol,eppl,EQV,0,1

/POST1
SET,3
    plnsol,eppl,EQV,0,1
    
/POST1
SET,4
    plnsol,eppl,EQV,0,1    
    
/POST1
SET,5
    plnsol,eppl,EQV,0,1   
    
/POST1
SET,6
    plnsol,eppl,EQV,0,1      
    
/POST1
SET,7
    plnsol,eppl,EQV,0,1        
PLNSOL,S,XY,0,1 
PLNSOL,EPTO,XY,0,1     !ʼ״̬µļӦ
    plnsol,eppl,EQV,0,1    
    plnsol,s,eqv,0,1  !ЧӦ
    






/POST1
LCDEF,1,1, ,          !1㲽Ϊ1
LCDEF,7,7, ,          !6㲽(ͶΪ100%)Ϊع2
LCASE,7,             !ع2
LCOPER,SUB,1, , ,     !ع2ȥع1
PLNSOL,U,sum,0,1          !ʼ״̬µλƳ
PLNSOL,EPTO,XY,0,1     !ʼ״̬µļӦ
    plnsol,eppl,EQV,0,1




9.5.2  ˮ״̬±ȶ
жԲͬͶµıȶԶ⣬ƪޣȫ״̬ıȶԽз뷽ӦõͶµȶԷ
Ҫ㱥ˮĸĸıֵȻҪǰ6.4.4еĺعӦ޸ģʾ
!*******loadcase.mac*******
/POST1
LCDEF,1,1, ,          !1㲽Ϊ1
LCDEF,2,7, ,          !6㲽(ͶΪ100%)Ϊع2
LCASE,2,             !ع2
LCOPER,SUB,1, , ,     !ع2ȥع1
ȻλơͼӦͼΪ
PLNSOL,U,sum,0,1          !ʼ״̬µλƳ


                
Mohr-Column׼Աˮ״̬µֲм㣬бˮ״̬µԲΪC=50kpa=30?£
!********Mohr.mac********
/post1
asel,s,,,1,2,1                     !ѡԪ
esla,s
etable,S1_1,S,1                  !ѡԪĵ13ӦΪԪ
etable,S3_1,S,3
sadd,S1,S3_1,,-1.00,1
sadd,S3,S1_1,,-1.00,1
sadd,smf_over,S1,S3,1.00,-1.00
cos_fai=cos(30/180*3.1415926)     !Ħǵ
sin_fai=sin(30/180*3.1415926)
constant=2*50000*cos_fai         
sadd,S1_S3,S1,S3,1,1
sexp,S1__S3,S1_S3,,1
sadd,smf_down,S1__S3,,sin_fai,,constant
sexp,s_m_f,smf_over,smf_down,1,-1
esel,all
/contour
pletab,s_m_f,avg                  !Ԫͼ
