!˳Ϊµеس
/PREP7
!嵥Ԫ
ET,1,shell181    
!                             
MP,EX,1,2.5e5            !嵯ģ
MP,PRXY,1,0.289          !岴ɱ
MP,DENS,1,7.85e-6        !ܶ
TB,BISO,1     
TBDATA,,512,0            !ǿ
!ʵ
R,1,16                   !ܺ
R,2,80                   !
R,3,12                   !м֧ܺ
R,4,10                   !֧ܺ
/PNUM,kp,1               !򿪹ؼɫ
/PNUM,line,1             !߲ɫ  
/PNUM,area,1             !ɫ
!
K,1,0,0,0                !ܹؼ1
K,2,1059,0,0             !ܹؼ2
K,3,-665,0,0             !ܹؼ3 
L,2,3                    !ӹؼ2͹ؼ3һ
LWPLAN,,1,0              !תƽ1Ϊz
PCIRC,169.5, ,0,360      !һԲ뾶Ϊ169.5
ADEL,1                   !ɾ1
ADRAG,2,3,4,5,,,1        !12345
NUMMRG,ALL               !ϲ
NUMCMP,ALL               !ѹ
ALLSEL,ALL               !ѡ
AGLUE,ALL                !ճе
!֧1
K,101,-131.5,0,0         !ؼ101
K,102,-769,127,-803      !ؼ102
L,101,102                !ӹؼ101͹ؼ102һ
LWPLAN,-1,14,1           !תƽ1Ϊz
PCIRC,116.5, ,0,360      !һԲ뾶Ϊ116.5
ADEL,5                   !ɾ5
ADRAG,15,16,17,18,,,14   !1415161718
NUMMRG,ALL               !ϲ
NUMCMP,ALL               !ѹ
ALLSEL,ALL               !ѡ
APTN,ALL                 !ཻ߷ָཻ       
ADELE,12,13              !ɾ1213
ADELE,16,17              !ɾ1617
AGLUE,ALL                !ճе
NUMMRG,ALL               !ϲ
NUMCMP,ALL               !ѹ
ALLSEL,ALL               !ѡ
!֧2
K,101,9,140,-889         !ؼ101
L,1,101                  !ӹؼ1͹ؼ101һ
LWPLAN,-1,39,1           !תƽ39Ϊz
PCIRC,116.5, ,0,360      !һԲ뾶Ϊ116.5
ADEL,11                  !ɾ11 
ADRAG,40,41,42,43,,,39   !3940414243
NUMMRG,ALL               !ϲ
NUMCMP,ALL               !ѹ
ALLSEL,ALL               !ѡ
APTN,ALL                 !ཻ߷ָཻ   
ADELE,18,19              !ɾ1819
ADELE,22,23              !ɾ2223  
AGLUE,ALL                !ճе
NUMMRG,ALL               !ϲ
NUMCMP,ALL               !ѹ
ALLSEL,ALL               !ѡ
!֧3
K,101,131.5,0,0              
K,102,643.5,99,-625
L,101,102
LWPLAN,-1,64,1
PCIRC,116.5, ,0,360
ADEL,17
ADRAG,65,66,67,68,,,64
NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
APTN,ALL
ADELE,24,25
ADELE,28,29
AGLUE,ALL
NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
!֧4
K,1001,104,840,-859  
L,1,1001
LWPLAN,-1,89,1
PCIRC,85, ,0,360
ADEL,23
ADRAG,90,91,92,93,,,89
NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
APTN,ALL
ADELE,28
ADELE,31,32
ADELE,35
ADELE,38,39
AGLUE,ALL
NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
!֧5
K,1001,103,-800.5,-850
L,1,1001
LWPLAN,-1,119,1
PCIRC,85, ,0,360
ADEL,31
ADRAG,120,121,122,123,,,119
NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
APTN,ALL
ADELE,36
ADELE,39,40
ADELE,43
ADELE,47,48
AGLUE,ALL
NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
!̣и֧
K,200,-42,-8,-587      
K,201,-38,-8,-616.5
L,201,200
LWPLAN,-1,149,1
PCIRC,1000, ,0,360
APTN,ALL
ADEL,40,43
ADEL,46,53
ADEL,56,59
ADEL,62,65
ADEL,87
ADEL,68,72
NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
!
K,1001,-615,0,-169.5  
K,1002,-615,0,386.5
K,1003,1008.7,0,953
K,1004,1008.7,0,-169.5
L,1001,1002
L,1002,1003
L,1003,1004
L,1001,1004
LFILLT,194,195,300
LFILLT,195,196,250
K,2001,761.5,0,604.5  
KWPLAN,-1,1,2001,121
KWPAVE,2001
PCIRC,68, ,0,360
AL,194,198,195,199,196,197
ASBA,40,39
LDIV,200, , ,2,0
LDIV,203, , ,2,0
APTN,ALL
AGLUE,ALL
NUMMRG,ALL
NUMCMP,ALL
ALLSEL,ALL
!񻮷
TYPE,1                        !õԪͱΪ1
MAT,1                         !ͱΪ1
REAL,1                        !ʵΪ1
ASEL,S,AREA,,2,6              !ѡ2-6
ASEL,A,AREA,,9,10             !ѡ9-10
ASEL,A,AREA,,12               !ѡ12
ASEL,A,AREA,,14,16            !ѡ14-16
ASEL,A,AREA,,36,37            !ѡ36-37
ASEL,A,AREA,,39               !ѡ39
SMRT,3                        !ܻ񻮷֣óߴȼΪ3
AMESH,ALL                     !ѡ
!м֧123񻮷
TYPE,1                        !õԪͱΪ1
MAT,1                         !ͱΪ1
REAL,3                        !ʵΪ3
ASEL,S,AREA,,7,8              !ѡ7-8
ASEL,A,AREA,,11,13,2          !ѡ1113
ASEL,A,AREA,,17,20            !ѡ17-20
ASEL,A,AREA,,25,29            !ѡ25-29
ASEL,A,AREA,,31,35,2          !ѡ313335
SMRT,3                        !ܻ񻮷֣óߴȼΪ3
AMESH,ALL                     !ѡ
!֧45񻮷
TYPE,1                   
MAT,1
REAL,4
ASEL,S,AREA,,1
ASEL,A,AREA,,21,24
ASEL,A,AREA,,30,34,2
SMRT,3                       !ܻ񻮷֣óߴȼΪ3
AMESH,ALL                    !ѡ
!񻮷
TYPE,1                   
MAT,1
REAL,2
ASEL,S,AREA,,38,40,2
SMRT,3                       !ܻ񻮷֣óߴȼΪ3
AMESH,ALL                    !ѡ
LSEL,S,LINE,,22,25           !(ϸֲ)ѡ22-25
LSEL,A,LINE,,29,30           !ѡ29-30
LSEL,A,LINE,,45,46           !ѡ45-46
LSEL,A,LINE,,59,62           !ѡ59-62
LSEL,A,LINE,,66,67           !ѡ66-67
LSEL,A,LINE,,81,85,2         !ѡ818385
LSEL,A,LINE,,86              !ѡ86
LSEL,A,LINE,,91,93,2         !ѡ9193
LSEL,A,LINE,,94,96,2         !ѡ9496
LSEL,A,LINE,,97              !ѡ97
LSEL,A,LINE,,109             !ѡ109
LSEL,A,LINE,,111,113,2       !ѡ111113
LSEL,A,LINE,,114             !ѡ114
LSEL,A,LINE,,120,123         !ѡ120-123
LSEL,A,LINE,,125             !ѡ125
LSEL,A,LINE,,202             !ѡ202
LSEL,A,LINE,,194,199         !ѡ194-199
LREFINE,ALL                  !ѡߵΧϸ
NUMMRG,ALL                   !ϲ
NUMCMP,ALL                   !ѹ
ALLSEL,ALL                   !ѡ
LSEL,S,LINE,,133,137,2       !ʩԼ
LSEL,A,LINE,,138
LSEL,A,LINE,,144,148,2
LSEL,A,LINE,,149
LSEL,A,LINE,,152,156,2
LSEL,A,LINE,,157
LSEL,A,LINE,,163,167,2
LSEL,A,LINE,,168
LSEL,A,LINE,,174,178,2
LSEL,A,LINE,,179
NSLL                        !ѡѡнڵ
D,ALL,ALL                   !ѡڵʩuxuyuzԼ
ALLSEL,ALL                  !ѡ
!ӳ
*GET,line6,line,6,leng           !6ĳֵֵline6
*GET,line8,line,8,leng           !8ĳֵֵline8
*GET,line10,line,10,leng         !10ĳֵֵline10
*GET,line11,line,11,leng         !11ĳֵֵline11
lleng=line6+line8+line10+line11  !ֵ:llengΪline6line8line10line11֮
!ӳ
*GET,line2,line,2,leng           !2ĳֵֵline2
*GET,line3,line,3,leng           !3ĳֵֵline3
*GET,line4,line,4,leng           !4ĳֵֵline4
*GET,line5,line,5,leng           !5ĳֵֵline5
rleng=line2+line3+line4+line5    !ֵ:rlengΪline2line3line4line5֮  
!  
*GET,line194,line,194,leng      !194ĳֵֵline194
*GET,line199,line,199,leng      !199ĳֵֵline199
erleng=line194+line199          !ֵ:erlengΪline194line199֮
SAVE                            !ģ
FINISH                          !˳ǰ
/SOLU                           !
ANTYPE,static                   !ѡ
NLGEOM,1                        !򿪴λЧ
ACEL,,,9.8                      !ʩٶ
!DOѭ彥ߺ
*DO,i,1,60
    pl=(27000+27000*(i-1))/lleng     !ӳܵѹpl
    LSEL,S,LINE,, 6,8,2
    LSEL,A,LINE,, 10,11
    SFL,ALL,PRES,pl                  !ѡʩѹpl
    ALLSEL,ALL                       !ѡ
    pr=(145000+145000*(i-1))/rleng   !ӳܵѹpr
    LSEL,S,LINE,,2,5            
    SFL,all,pres,pr                  !ѡʩѹpr
    ALLSEL,all
    per=(150000+150000*(i-1))/erleng !ܵѹpr

    LSEL,S,LINE,,194,199,5
    SFL,ALL,PRES,per                 !ѡʩѹper
    ALLSEL,ALL                       !ѡ
    OUTRES,ALL,59                    !Ŀÿ59һΰһӲ
    TIME,i                           !ضӦյʱΪi
    NSUBST,50,,,1                    !غӲΪ50,һֹʱΪһصʼʱ
    KBC,0                            !غɰ
    AUTOTS,ON                        !ʹԶʱ䲽
    NCNV,2,0,0,0,0                   !òʱֹansys,˳
    NEQIT,50                         !ÿһӲƽ
    LSWRITE,i                        !дiغɲļ
*ENDDO
LSSOLVE,1,60,1                       !iغɲļִ
FINISH
SAVE,'jiedian',,,all                 !ļ

/POST1
SET,list
SET,,,,,, ,25                        !SET=25
PLNSOL,U,SUM,0,1                     !۲λƵֵͼ
PLNSOL,S,EQV,0,1                     !۲Von MisesЧӦֵͼ
SET,,,,,, ,47                        !SET=47
PLNSOL,U,SUM,0,1                     !۲λƵֵͼ
PLNSOL,S,EQV,0,1                     !۲Von MisesЧӦֵͼ
ESEL,S,REAL,,1                       !ѡ
ETABLE,zhuguan,S,EQV                 !嵥Ԫ洢Von MisesЧӦ
ESORT,ETAB,zhuguan,0,0, ,            !ԵԪzhuguan洢Von MisesЧӦ
*GET,n1,SORT,0,IMAX                  !ܵVon MisesЧӦֵıŸN1
ESEL,S,REAL,,2                       !ѡ
ETABLE,erban,S,EQV                   !嵥Ԫ洢Von MisesЧӦ
ESORT,ETAB,erban,0,0,,               !ԵԪerban洢Von MisesЧӦ
*GET,n2,SORT,0,IMAX                  !Von MisesЧӦֵıŸN2
ESEL,S,REAL,,3                       !ѡм֧
ETABLE,midzhiguan,S,EQV              !嵥Ԫ洢м֧ܵVon MisesЧӦ
ESORT,ETAB,midzhiguan,0,0, ,         !ԵԪmidzhiguan洢Von MisesЧӦ
*GET,n3,SORT,0,IMAX                  !м֧ܵVon MisesЧӦֵıŸN3
eSEL,S,REAL,,4                       !ѡ֧
ETABLE,dszhiguan,S,EQV               !嵥Ԫ洢֧ܵVon MisesЧӦ
ESORT,ETAB,dszhiguan,0,0, ,          !ԵԪdszhiguan洢Von MisesЧӦ
*GET,n4,SORT,0,IMAX                  !֧ܵVon MisesЧӦֵıŸN4
ALLSEL,ALL
/POST26                              !ʱ̺
File,jidian,rst,.        !򿪽ļ
ANSOL,2,10884,S,EQV,zhuguans         !2洢N1Von MisesЧӦ
ANSOL,3,8991,S,EQV,erbans            !3洢N3Von MisesЧӦ
ANSOL,4,20532,S,EQV,midzhiguans      !4洢N3ĵVon MisesЧӦ
ANSOL,5,13949,S,EQV,dszhiguans       !5洢N4Von MisesЧӦ
XVAR,1                               !Ϊʱ
PLVAR,2,3,4,5                          !ܡ塢м֧ܡ֧Von MisesЧӦʱı仯ͼ

