怎样取出圆柱壳的变形后的径向最大变形

bluecmy

各位高手，请问一个多层圆柱壳，
怎样取出各层的变形后的径向最大变形？使用什么命令
模型如下：
/units,si
T1=0.056
T2=0.284
layer=3
R=0.678-0.5*T1
L=12*(R+3*T1+2*T2)
d=0.5*(R+3*T1+2*T2)
opint=1500000


/prep7
et,1,mesh200
et,2,solid95

keyopt,1,1,6
keyopt,1,2,0

MP,EX,1,2.1e9
MP,PRXY,1,0.3
MP,EX,2,0.7e9
MP,PRXY,2,0.3

cyl4,0,0,R   
cyl4,0,0,R+T1
cyl4,0,0,R+T1+T2
cyl4,0,0,R+2*T1+T2
cyl4,0,0,R+2*T1+2*T2
cyl4,0,0,R+3*T1+2*T2
APTN,all
adele,1

LSTR,      23,       3  
LSTR,       1,      21  
LSTR,      22,       2  
LSTR,       4,      24  
FLST,2,5,5,ORDE,2   
FITEM,2,7   
FITEM,2,-11
FLST,3,4,4,ORDE,2   
FITEM,3,25  
FITEM,3,-28
ASBL,P51X,P51X  
aglue,all
NUMCMP,all

allsel
*do,i,1,2*layer-1,2
*do,j,0,3
aSEL,s,,,4*i-j
aATT,1, , 1, 0
*enddo
*enddo

allsel
*do,i,2,2*layer-1,2
*do,j,0,3
aSEL,s,,,4*i-j
aATT,2, , 1, 0
*enddo
*enddo

allsel
/PNUM,MAT,1
/REPLOT


FLST,5,24,4,ORDE,2  
FITEM,5,1   
FITEM,5,-24  
LSEL, , , ,P51X  
LESIZE,all, , ,8, , , , ,1       !环向

allsel
lsel,s,,,25
*do,j,0,3
*do,i,25+j,25+j+8*(layer-1),8
lsel,a,,,i
*enddo
*enddo
LESIZE,all, , ,1, , , , ,1                     !硬层
  
allsel
lsel,s,,,29
*do,j,0,3
*do,i,29+j,29+j+8*(layer-2),8
lsel,a,,,i
*enddo
*enddo
LESIZE,all, , ,2, , , , ,1                     !软层  
allsel
AMESH,all

K,100,0,0,0
K,101,0,0,L
LSTR,      100,      101
type,2
EXTOPT,ESIZE,80,0,
EXTOPT,ACLEAR,1
EXTOPT,ATTR,1,1,0   
ESYS,0  
VDRAG,all, , , , , ,45

jxxansys

可以用圆柱坐标系做为结果坐标系