求助：厚壁圆筒的自增强

Leeyajiang

我做的是厚壁圆筒自增强模拟。求解完成后，从路径图上看结果是正确的了，但从云图上看好像不符合轴对称厚壁圆筒的应力分布规律，可能我的加载约束有问题。那位高手做过厚壁圆筒的分析，帮帮忙看看怎么改进一下。
　　附上我求解的图和命令流。

Leeyajiang

/BATCH  
/COM,ANSYS RELEASE  9.0    UP20041104       10:12:59    11/18/2006
/input,menust,tmp,'',,,,,,,,,,,,,,,,1   
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1   
/REPLOT,RESIZE  
WPSTYLE,,,,,,,,0
/PREP7  
!*  
ET,1,SOLID45
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.024e5
MPDATA,PRXY,1,,0.3  
TB,BISO,1,2,2,  
TBTEMP,20   
TBDATA,,960,2.1e4,,,,   
TBTEMP,400  
TBDATA,,720,1.99e4,,,,  
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,400
MPDE,EX,1   
MPDE,PRXY,1
MPDATA,EX,1,,2.024E+005
MPDATA,EX,1,,1.775e5
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
TBDE,BISO,1,,,  
TB,BISO,1,2,2,  
TBTEMP,20   
TBDATA,,960,21000,,,,   
TBTEMP,400  
TBDATA,,720,19900,,,,      !define material

*SET,h,0.5  
*SET,R1       , 0.150   
*SET,R2       , 0.255   
CYLIND,r2,r1,0,h,-90,90,   !builded modle


FINISH  
/SOL
FINISH  
/PREP7  
FLST,5,2,4,ORDE,2   
FITEM,5,6   
FITEM,5,8   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,15, , , , ,1
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,11  
FITEM,5,-12
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,40, , , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,7   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,40, , , , ,1
!*  
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*                                 !meshed



FINISH  
/SOL
FLST,2,1,4,ORDE,1   
FITEM,2,10  
!*  
/GO
DL,P51X, ,ALL,  
FLST,2,2,5,ORDE,2   
FITEM,2,5   
FITEM,2,-6  
DA,P51X,SYMM                    !constaint


!load 1
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,313.91
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-166.09
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,1  
LSWRITE,1,  


!load 2
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,329.61  
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-174.39
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,2  
LSWRITE,2,  



!load 3
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,346.09
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-183.12
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,3  
LSWRITE,3,  


!load 4
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,363.39
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-192.27
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,4  
LSWRITE,4,  


!load 5
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,381.56
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-207.9
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,5  
LSWRITE,5,  


!load 6
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,400.64  
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-211.98
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,6  
LSWRITE,6,  


!load 7
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,420.67
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-222.58
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,7  
LSWRITE,7,  


!load 8
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,441.7
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-233.71
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,8  
LSWRITE,8,  


!load 9
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,470.19  
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-245.39
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,9  
LSWRITE,9,  


!load 10
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,329.13
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-174.14
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,10  
LSWRITE,10,  



!load 11
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,230.39
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-121.9
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,11  
LSWRITE,11,  


!load 12
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,161.27
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-85.33
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,12
LSWRITE,12,  



!load 13
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,112.89
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-59.73
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,13
LSWRITE,13,  


!load 14
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,79.02
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,-41.81
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,14
LSWRITE,14,  


!load 15
FLST,2,1,5,ORDE,1   
FITEM,2,4   
/GO
!*  
SFA,P51X,1,PRES,0
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
/GO
!*  
SFA,P51X,1,PRES,0
ANTYPE,0
NLGEOM,1
NSUBST,10,20,5  
AUTOTS,1
TIME,15
LSWRITE,15,

jwooh

做的挺不错的，可惜我不太懂塑性方面的东西，请别人给你解答吧！