FINISH
/CLEAR
/TITLE,YLRQ
/FILN,E653
!ϲ,λΪʵλ
*AFUN,DEG                                                ! Ƕȵλ
!ģͲ  
*SET,DI,3660E-3                                          ! Ͳھ
*SET,L,43000E-3                                          ! 
*SET,TN,26E-3                                            ! Ͳ幫ƺ 
*SET,A,6860E-3                                           ! ߵľ
*SET,B,455E-3                                            ! 
*SET,THETA,135                                           ! 
*SET,C,762E-3                                            ! 
*SET,TD,40E-3                                            ! 
*SET,H1,500E-3                                           ! H1 
*SET,H2,1500E-3                                          ! H2
*SET,H3,1750E-3                                             ! 
!ϲ
*SET,DENSITY,7850                                        ! ܶ  
*SET,EXX,2.01E11                                         ! ϵģ
*SET,PRXY,0.3                                            ! ϲɱ
*SET,P1,1E6                                              ! ѹ
!************************
!*******ǰ***********
!************************
/PREP7
!嵥Ԫ
ET,1,SOLID95                                             ! 95ŵԪ
MP,EX,1,EXX                                              ! 嵯ģ
MP,NUXY,1,PRXY                                           ! 岴ɱ
MP,DENS,1,DENSITY                                        ! ܶ
!ģ
CYLIND,DI/2,DI/2+TN,0,L/2,-90,90                         ! Ͳ
SPHERE,DI/2+TN,DI/2,-90,90,                              ! ͷ
VSBW,ALL                                                 ! ùƽз
VDELE,4, , ,1                                            ! ɾõͷ
CYLIND,DI/2+TN,DI/2+TN+TD,A-C/2,A+C/2,-90,-90+THETA/2    ! 
WPOFF,0,0,A-B/2                                          ! ƽƽǵ
BLC4,0,-DI/2-TN-TD-H1,H3,H2,B                            ! 
VPTN,ALL                              ! 廥
ALLS                                  ! ȫѡ
LGLUE,ALL                             ! ճϸ
VDELE,5,,,1                           ! ɾ
CYLIND,0,DI/4,0,-DI*4,-90,90          ! Ͳ
VPTN,ALL                              ! 廥
VDELE,2,3,1,1                         ! ɾ弰丽
!ùƽʷ
WPROT,0,90,0                          ! תƽ
VSEL,S,VOLU,,4,5,1,                   ! ѡзͷ
VSBW,ALL                              ! ùƽзַͷ
ALLSEL                                ! ȫѡ
WPROT,0,-90,0                         ! תƽ
KWPAVE,20                             ! ƽ涨λڹؼ20
VSBW,ALL                              ! ùƽз
KWPAVE,22                             ! ƽ涨λڹؼ22
VSBW,ALL                              ! ùƽз
KWPAVE,27                             ! ƽ涨λڹؼ27
VSBW,ALL                              ! ùƽз
KWPAVE,31                             ! ƽ涨λڹؼ31
VSBW,ALL                              ! ùƽз
WPROT,0,0,90                          ! תƽ
KWPAVE,43                             ! ƽ涨λڹؼ43
VSBW,9                                ! ùƽз֧Ų
VGLUE,ALL                             ! ճ

!,ѡ
VSEL,S,VOLU,,14,17,3
VSEL,A,VOLU,,6,16,10
CM,INTANK,VOLU                        ! ӴͰڱ
VSEL,S,VOLU,,13,15,2
CM,BACKSTOP,VOLU                      ! 
VSEL,S,VOLU,,5,7,2
VSEL,A,VOLU,,10,12,2
CM,GASKET,VOLU                        ! 
VSEL,S,VOLU,,1,2,1
VSEL,A,VOLU,,3,8,5
CM,CLOSE,VOLU                         ! ͷ
VSEL,S,VOLU,,4,11,7
CM,TANK,VOLU                          ! ѹͰ

!
ALLSEL                                ! ȫѡ
CMSEL,S,INTANK                        ! ѡӴͰڱ
ESIZE,0.02                            ! 趨ԪС
VSWEEP,ALL                            ! ӳ񻮷
CMSEL,S,GASKET                        ! ѡе
ESIZE,0.02                            ! 趨ԪС
VSWEEP,ALL                            ! ӳ񻮷
CMSEL,S,BACKSTOP                      ! ѡа
ESIZE,0.06                            ! 趨ԪС
VSWEEP,ALL                            ! ӳ񻮷
ESIZE,0.2                             ! 趨ԪС
CMSEL,S,CLOSE                         ! ѡзͷ
VSWEEP,ALL                            ! ӳ񻮷
CMSEL,S,TANK                          ! ѡѹͰ
ESIZE,0.4                             ! 趨ԪС
VSWEEP,ALL                            ! ӳ񻮷
ALLSEL                                ! ȫѡ
WPCSYS,-1,0                           ! תԭʼϵ
LOCAL,11,0,0,0,L/2,                   ! ֲϵ
CSYS,11                               ! ֲϵ
VSYMM,Z,ALL, , , ,0,0                 ! 
NUMMRG,ALL                            ! ϲ
NUMCMP,ALL                            ! ѹ
FINISH

!************************
!**************
!************************
/SOLU
ALLSEL
ACEL,0,9.8,0,                         ! ʩٶ 
!*********ע⣬ANSYSʩӵΪٶȷʵ෴*************
NSEL,S,LOC,Z,0,L
NSEL,R,LOC,Y,-DI/2-TN-TD-H1           ! ѡǰ˰
D,ALL,ALL                             ! ʩӹ̶Լ
ALLSEL
NSEL,S,LOC,Z,-L,0
NSEL,R,LOC,Y,-DI/2-TN-TD-H1           ! ѡк˰
D,ALL,UY,0                            ! ̶YԼ
ALLSEL
NSEL,S,LOC,X,0                        ! ѡжԳнڵ
DSYM,SYMM,X,                          ! ʩӶԳλԼ
CSYS,0                                ! ѿϵ           
LSEL,S,LOC,Z,0,L                      ! ѡгͷ֮Ͱ弰߶
CSYS,1                                ! ϵ
LSEL,R,LOC,X,0,DI/2                   ! ѡͰڱ߶
ASLL,S,1                              ! ѡͲڱ
ASEL,A,AREA,,10,11,1
ASEL,A,AREA,,28,29,1
ASEL,A,AREA,,92,118,26
ASEL,A,AREA,,83,86,3                  ! ѡͰڱڼͷڱ
NSLA,S,0                              ! ѡڱнڵ
SF,ALL,PRES,P1                        ! ʩӾѹ
SF,ALL,PRES,%P2%                      ! ʩӾˮѹ
ALLSEL                                ! ȫѡ
SOLVE                                 ! 
FINISH

!************************
!******************
!************************
/POST1
!*******************************************
PLNSOL,S,EQV,0,1                      ! ʾЧӦֲ
PLNSOL,U,Y,2,1                        ! ʾYλ
!*************ֲ******************************
RSYS,1                                ! ϵΪϵ
ALLSEL                                ! ȫѡ
PLNSOL,S,SY,0,1                       ! ʾӦֲ
CMSEL,S,GASKET                        ! ѡе
CMSEL,A,BACKSTOP                      ! ѡа
ESLV,S                                ! ѡڵ+ϵĵԪ
PLNSOL,S,SY,0,1                       ! ʾӦֲ

!**************Ӧ**********************************
CSYS,0                                ! ѿϵ
RSYS,0                                ! ѿϵ
ALLSEL                                ! ȫѡ
!****************·1****************
NODE1=NODE(0,DI/2,0)                  ! ڵNODE1
NODE2=NODE(0,DI/2,L)                  ! ڵNODE2
PATH,PATH_1,2,30,40,                  ! ·PATH_1
PPATH,1,NODE1                         ! ͨڵλö·ؼNODE1
PPATH,2,NODE2                         ! ͨڵλö·ؼNODE2
PDEF, ,S,EQV,AVG                      ! ЧӦӳ䵽·
PDEF,SY,S,Y,AVG                       ! Ӧӳ䵽·
PDEF,SX,S,X,AVG                       ! Ӧӳ䵽·
PDEF,SZ,S,Z,AVG                       ! Ӧӳ䵽· 
/AXLAB,X,DISTANCE                     ! Xʾ
/AXLAB,Y,STRESS                       ! YʾӦ
PLPATH,SEQV,SY,SX,SZ                  ! ʾ·ϵӦֲ

!****************·2****************
NODE3=NODE(0,-DI/2,0)                 ! ڵNODE3
NODE4=NODE(0,-DI/2,L)                 ! ڵNODE4
PATH,PATH_2,2,30,40,                  ! ·PATH_2
PPATH,1,NODE3                         ! ͨڵλö·ؼNODE3
PPATH,2,NODE4                         ! ͨڵλö·ؼNODE4
PDEF, ,S,EQV,AVG                      ! ЧӦӳ䵽·
PDEF,SY,S,Y,AVG                       ! Ӧӳ䵽·
PDEF,SX,S,X,AVG                       ! Ӧӳ䵽·
PDEF,SZ,S,Z,AVG                       ! Ӧӳ䵽· 
/AXLAB,X,DISTANCE                     ! Xʾ
/AXLAB,Y,STRESS                       ! YʾӦ
PLPATH,SEQV,SY,SX,SZ                  ! ʾ·ϵӦֲ

!****************·3****************
TOP_NODE=NODE (0,DI/2,A)              ! ڵTOP_NODE
BOT_NODE=NODE (0,-DI/2,A)             ! ڵBOT_NODE
CSYS,1                                ! ϵ
RSYS,1                                ! ϵΪϵ
PATH,PATH_3,2,30,20,                  ! ·PATH_3  
PPATH,1,TOP_NODE                      ! ͨڵλö·ؼ1
PPATH,2,BOT_NODE                      ! ͨڵλö·ؼ2
PDEF,SY,S,Y,AVG                       ! Ӧӳ䵽·
PDEF,SX,S,X,AVG                       ! Ӧӳ䵽· 
PDEF,SXY,S,XY,AVG                     ! Ӧӳ䵽· 
/AXLAB,X,DISTANCE                     ! Xʾ
/AXLAB,Y,STRESS                       ! YʾӦ
PLPATH,SY,SX,SXY                      ! ʾ·ϵӦֲ

!****************·4****************
CSYS,0                                ! ѿϵ
TOP_NODE2=NODE (0,DI/2,L-A)           ! ڵTOP_NODE
BOT_NODE2=NODE (0,-DI/2,L-A)          ! ڵBOT_NODE
CSYS,1                                ! ϵ
RSYS,1                                ! ϵΪϵ
PATH,PATH_4,2,30,20,                  ! ·PATH_4  
PPATH,1,TOP_NODE2                     ! ͨڵλö·ؼ1
PPATH,2,BOT_NODE2                     ! ͨڵλö·ؼ2
PDEF,SY,S,Y,AVG                       ! Ӧӳ䵽·
PDEF,SX,S,X,AVG                       ! Ӧӳ䵽· 
PDEF,SXY,S,XY,AVG                     ! Ӧӳ䵽· 
/AXLAB,X,DISTANCE                     ! Xʾ
/AXLAB,Y,STRESS                       ! YʾӦ
PLPATH,SY,SX,SXY                      ! ʾ·ϵӦֲ

!**************ӦԻ·1****
CSYS,0                                ! ѿϵ
NODE5=NODE(0,-DI/2,A)                 ! ڵNODE5
NODE6=NODE(0,-DI/2-TN,A)              ! ڵNODE6
PATH,PATH_5,2,30,20,                  ! ·PATH_5  
PPATH,1,NODE5                         ! ͨڵλö·ؼNODE5
PPATH,2,NODE6                         ! ͨڵλö·ؼNODE6
PLSECT,S,X,-1,0                       ! ӦԻ

!**************ӦԻ·2****
CSYS,0                                ! ѿϵ
NODE7=NODE(0,-DI/2,A/2)               ! ڵNODE7
NODE8=NODE(0,-DI/2-TN,A/2)            ! ڵNODE8
PATH,PATH_6,2,30,20,                  ! ·PATH_6  
PPATH,1,NODE7                         ! ͨڵλö·ؼNODE7
PPATH,2,NODE8                         ! ͨڵλö·ؼNODE8
PLSECT,S,X,-1,0                       ! ӦԻ
