复合材料圆柱壳弹性屈曲求解

笑天蛋

目前需要计算一个复合圆柱壳的弹性屈曲载荷，但是自己的程序一输入ANSYS就跳出来了，不知道是什么原因。请问各位路过的高手帮我看看，下面是命令流
finish
/clear
/filn, PSP250-3D
/title, Buckling analysis of PSP
! ********* 参数设定*********
rst=1.2/2
tt=15.5
ratio=0.6
Ro=500/2
L=6*Ro     
N=320/2   
Ri=Ro-tt
R=ratio*(Ro-Ri)+Ri
Rj=R-2*rst
Rk=R+2*rst
tm=Ro-Rk      
tc=Rk-Rj        
tn=Rj-Ri            
nhoop=20                 
naxial=80               
p=1                     
!**************************前处理啦**********************************
/PREP7  
!一号材料，复合层材料
ET,1,SOLID46
*SET,_RC_SET,1,
R,1
RMODIF,1,1,1,0,0,0  
RMODIF,1,7,0
RMODIF,1,13,1,0,2.4,
TB,ANEL,1,1,21,1
TBTEMP,0
TBDATA,,0.0003,-0.0004,0,0,0,0  
TBDATA,,0.0009,-0.0001,0,0,0,0.0004
TBDATA,,0,0,0,0.0017,0,0
TBDATA,,0.0011,0,0.0001,,,  
!****************二号材料
!****************HDPE材料***********************
ET,2,SOLID45
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,2,,1002   
MPDATA,PRXY,2,,0.45
!******************建模***************************
K,1,0,0,0,                       !定义筒体轴心线关键点
K,2,0,0,L,                       !定义筒体轴心线另一关键点
L,1,2                            !生成轴心线
LESIZE,all, , ,naxial, , , , ,1  !设定筒体轴向剖分数
CYL4,0,0,Ri,360,Rj,360
CYL4,0,0,Rj,360,Rk,360
CYL4,0,0,Rk,360,Ro,360
VDRAG,all, , , , , ,1           !利用拖拉命令生成圆柱壳
VGLUE, all,, , , , , , ,
!******************glue mission is accomplished*********************
lsel,s,loc,z,0                   !选择圆周线各线段
LESIZE,all, , ,nhoop,, , , ,1    !设定圆周线剖分数

VSWEEP,all,,
!*********************施加约束了*************************
nsel,s,loc,z,0                   !选择z=0端面 !这就加约束了
D,all,uz,                        !约束z=0端面的一个点轴向位移
nsel,s,loc,z,0                   !选择z=0端面各节点
nsel,a,loc,z,L                   !选择z=L端面各节点
csys,1                           !激活总体柱坐标
NROTAT,all                       !旋转所有已选节点的节点坐标系
D,all,UY                         !施加环向位移约束，已保证该截面保持为圆
fini
/SOLU
ANTYPE,STATIC                  !静态分析
PSTRES,ON                      !激活预应力影响
SFA,22,1,PRES,p               !在圆柱壳的外表面施加外压
SFA,23,1,PRES,p
SFA,24,1,PRES,p
SFA,25,1,PRES,p
Alls                           !全选
SOLVE                          !求解
FINISH                           !结束并退出求解器
/SOLU                          !重新进入求解器
ANTYPE,BUCKLE                  !指定分析类型为屈曲分析
BUCOPT,LANB,1                  !利用Block Lanczos法进行结算，提取一节模态
MXPAND,1                       !设定模态扩展数为1
SOLVE                          !求解

是想求弹性屈曲载荷，程序中关于复合层的材料定义是不是有问题，不太懂
输入柔度矩阵时，D11，D22，D33中1,2,3分别和系统坐标系的X,Y,Z对应吗？
定义时是希望把1设为环向，2设为z向，3设为r向，这个怎么弄？......
谢谢大家啦！

笑天蛋

无人来看，自顶一下吧
差不多已经发现问题了，程序里面掉了两行。。。。。。
汗

fruitboy

这个我没弄出来呢？