9.2  
2ǰĵԪ͡ڸԣģֻPlane 42һֵԪڱԶάƽӦ⣬轫Plane 42KEYOPT(3)Ϊ2Ӧ£
/prep7
et,1,plane42
KEYOPT,1,3,2     !ƽӦѡ
3ģеĲϲģй7ֲϣҺͲͬͶµĻĲΡĲε£
mp,ex,1,20e9                               !嵯ģ
mp,prxy,1,0.22                              !岴ɱ
mp,dens,1,2600                             !ܶ
mp,ex,2,3e9                                !ʼ
mp,prxy,2,0.35
mp,dens,2,2140
mp,ex,3,2.7e9                               !ͶΪ35%
mp,prxy,3,0.35
mp,dens,3,2180
mp,ex,4,2.4e9                               !ͶΪ50%
mp,prxy,4,0.35
mp,dens,4,2210
mp,ex,5,2.1e9                               !ͶΪ65%
mp,prxy,5,0.35
mp,dens,5,2240
mp,ex,6,1.8e9                               !ͶΪ80%
mp,prxy,6,0.35
mp,dens,6,2270
mp,ex,7,1.5e9                               !ͶΪ100%ĵģ
mp,prxy,7,0.35                              !ͶΪ100%Ĳɱ
mp,dens,7,2300                             !ͶΪ100%ܶ
save                                       !洢ݿ

4ģͲֵԪģ͵ļ״ΪӣʲԵϵĽģʽɵ㣬Ȼ棬󻮷ֵԪ
k,1,0,0,0                                  !ɹؼ
k,2,0,505,0
k,3,562,505,0
k,,621,459,0
k,,658,433,0
k,,693,393,0
k,,802,313,0
k,,850,303,0
k,,892,273,0
k,,913,253,0
k,,930,243,0
k,,1034,233,0
k,,1186,228,0
k,,1216,223,0
k,,1232,208,0
*do,i,1,14                                !ѭķʽֱߣ14
  l,i,i+1
*enddo
k,16,660,387,0                            !ٴؼ
k,,770,248,0
k,,827,230,0
k,,888,217,0
k,,930,216,0
k,,1034,213,0
k,,1186,209,0
l,4,16                                    !ٸݹؼ㴴ֱ
l,16,17
l,17,18
l,18,19
l,19,20
l,20,21
l,21,22
l,22,15
nummrg,all                                 !ϲԪ
numcmp,all                                 !ѹԪر
k,23,1243,198,0                             !ٴؼ
k,,612,380,0
k,,745,214,0
k,,884,181,0
k,,1242,174,0
l,3,24                                     !ٸݹؼ㴴ֱ
l,24,25
l,25,26
l,26,27
l,27,23
l,23,15
nummrg,all                                 !ϲԪ
numcmp,all                                 !ѹԪر
k,28,1610,198,0                             !ٴؼ
k,,1610,0,0
k,,1243,0,0
k,,884,0,0
k,,745,0,0
k,,612,0,0
l,6,16                                     !ٸݹؼ㴴ֱ
l,7,17
l,8,18
l,9,19
l,11,20
l,13,22
l,16,24
l,17,25
l,19,26
l,24,33
l,25,32
l,26,31
l,27,30
l,1,33
l,30,31
l,31,32
l,32,33
l,23,28
l,28,29
l,29,30
save
al,4,5,15,29                                 !ߴ
al,29,6,16,30
al,7,17,30,31
al,8,18,31,32
al,9,10,19,32,33
al,11,12,20,21,33,34
al,13,14,22,34
al,3,15,35,23
al,35,16,24,36
al,36,17,18,37,25
al,19,20,21,22,26,27,28,37
al,1,2,23,38,42
al,38,39,24,45
al,39,40,25,44
al,40,41,26,43
al,46,47,48,27,41
allsel,all                                     !ѡԪ
type,1                                      !ѡҪֵԪ
mat,1                                       !ѡҪֵԪĲģ
mshape,1,2d                                 !óƽεԪ
mshkey,0                                    !óɷʽ
ESIZE,10,0,                                  !õԪС
amesh,1,11,1                                 !111ֵԪ
ESIZE,,,                                     !õԪСΪĬֵ
allsel,all                                     !ѡԪ
amesh,12,16,1                                !1216
save
finish
!ֺúԪͼ9-2ʾ

   
9.4  

1ʩӱ߽á
/solu                                       !
eplot                                       !ʾԪͼ
save                                        
nsubst,20                                    !ָӲ
nlgeom,on                                   !ָ
nropt,full                                    !ónewton-lapnace
pred,on                                     !÷Ԥ
lnsrch,on                                    !
outres,all,last                                 !ֻһ
!ʩԼ
nsel,s,loc,x,0                                  !ѡڵ
d,all,ux                                      !ԼڵXλ
nsel,s,loc,x,1610                               !ѡڵ
d,all,ux                                      !ԼڵXλ
nsel,s,loc,y,0                                  !ѡڵ
d,all,all                                      !ԼڵXYλ
allsel,all                                      !ѡԪ
acel,0,10,0                                    !ʩٶ
allsel,all                                      !ѡԪ
save
finish
2Ӧ
/solu                                   !
time,1                                  !һⲽ
asel,s,,,1,7,1                             !ѡ
esla,s                                  !ѡϵԪ
mpchg,2,all                             !˲ֻ嵥ԪĲı
allsel,all                                !ѡԪ
solve                                   !
save                                    !ݿ
finish                                   !һ˵˳
3ⲻͬͶµӦλƳ
!******ͶΪ35ʱ********
/solu                                   !
antype,,rest                              !
time,2                                  !õ2ⲽ
asel,s,,,1,7,1                             !ѡ
esla,s                                   !ѡϵԪ
mpchg,3,all                              !ıѡĵԪĲϲ
allsel,all                                 !ѡԪ
solve                                   !
save                                    !ݿ
finish                                   !һ˵˳
!******ͶΪ50ʱ********
!עͬͶΪ35ʱ
/solu
antype,,rest
time,3
asel,s,,,1,7,1
esla,s
mpchg,4,all
allsel,all
solve 
save
finish
!******ͶΪ65ʱ********
/solu
antype,,rest
time,4
asel,s,,,1,7,1
esla,s
mpchg,5,all
allsel,all
solve 
save
finish
!******ͶΪ80ʱ********
/solu
antype,,rest
time,5
asel,s,,,1,7,1
esla,s
mpchg,6,all
allsel,all
solve
save
finish
!******ͶΪ100ʱ********
/solu
antype,,rest
time,6
asel,s,,,1,7,1
esla,s
mpchg,7,all
allsel,all
solve
save
finish
9.5  
9.5.1  ʼ״̬±ȶ
ɺ󣬶һڸñµӦλƼӦͼӦͼڽ/POST1еģʶ֮ǰͨúЩͼΪ
/POST1
SET,1
PLNSOL,S,Y,0,1          !ʼ״̬µӦ
PLNSOL,U,Y,0,1          !ʼ״̬µλƳ
PLNSOL,S,XY,0,1         !ʼ״̬µļӦ
PLNSOL,EPTO,XY,0,1     !ʼ״̬µļӦ

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

                
Mohr-Column׼Աˮ״̬µֲм㣬бˮ״̬µԲΪC=50kpa=30?£
!********Mohr.mac********
/post1
asel,s,,,1,7,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                  !Ԫͼ
