多层中空圆筒的屈曲分析（压杆稳定）可不可以用轴对称方法分析（vm128）

bluecmy

以下是我的命令 （运行到最后时，结果是ANSYS自动关闭:@@ 郁闷）

/prep7  
et,1,plane42
keyopt,1,3,1

MP,Ex,1,,6.33e12
MP,prxy,1,,0.3
   

MP,EX,2,,0.021e12  
MP,PRXY,2,,0.3  


k,1,0.17,,0   
k,2,0.312,,0
k,3,0.368,,0
k,4,0.652,,0
k,5,0.708,,0
k,6,0.992,,0
k,7,1.048,,0
k,8,1.19,,0
k,9,0.17,11.9,
k,10,0.312,11.900,
k,11,0.368,11.900,
k,12,0.652,11.900,
k,13,0.708,11.900,
k,14,0.992,11.900,
k,15,1.048,11.900,
k,16,1.19,11.900,

k,100,0,0,0
k,101,0,11.900,

a,1,2,10,9   
a,3,11,10,2   
a,3,4,12,11
a,4,5,13,12
a,5,6,14,13
a,6,7,15,14
a,7,8,16,15

GPLOT   
  
aglue,all


asel,s,loc,x,0.312,0.368
asel,a,loc,x,0.652,0.708
asel,a,loc,x,0.992,1.048
AATT,1, ,1,0,                   !!材料属性

allsel,all

asel,s,loc,x,0.17,0.312
asel,a,loc,x,0.368,0.652
asel,a,loc,x,0.708,0.992
asel,a,loc,x,1.048,1.190
AATT,2, ,1,0,                    !!材料属性

allsel,all
/PNUM,MAT,1
/REPLOT

LSEL,S,LOC,y,1,11
LESIZE,all,0.084, , , , , , ,1
esize,0.056
type,1
mshape,0,2d
mshkey,1
amesh,all                       !!!!!生成网格


                        
/SOLU   
PSTRES,ON                 ! CALCULATE PRESTRESS EFFECTS

nsel,s,loc,y,0
D,all,ALL                 ! 约束

allsel
nsel,s,loc,y,11.900
*get,no_nodes,node,,count
F,all,FY,-1/no_nodes          !加单位载荷
OUTPR,,1
SOLVE
FINISH
/SOLU   
ANTYPE,BUCKLE             ! 屈曲分析
BUCOPT,SUBSP,1            !
MXPAND,1                  !
SOLVE
*GET,FCR,MODE,1,FREQ
*status,parm
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'Fcr '
LABEL(1,2) = 'lb '
*VFILL,VALUE(1,1),DATA,38.553
*VFILL,VALUE(1,2),DATA,FCR
*VFILL,VALUE(1,3),DATA,ABS(FCR/38.553)

bluecmy

为什么没人回帖子？我的问题太没水平？

venture

多层中空圆筒的屈曲分析（压杆稳定）可不可以用轴对称方法分析.
肯定是不太合适的，VM128的例子和你的有本质区别，例子上bar用plane单元来模拟是没问题的。你的圆筒用轴对称是不合适的，去找本屈曲方面的书看看基础知识。圆筒分长圆筒和短圆筒，屈曲模态是不相同的。对圆筒，不说你轴对称模型，就是1/4，1/2模型在没有确定把握的情况下也尽量不要用，因为屈曲模态不一定是对称的。所以屈曲分析在一般情况下最好用整体模型。

bluecmy

这是我做的体模型，网格很大，只是初步的计算。
/prep7  
et,1,solid45
MP,Ex,1,6.33e12
MP,prxy,1,0.3
   
MP,EX,2,0.021e12  
MP,PRXY,2,0.3  
  
k,1,0.17,,0   
k,2,0.312,,0
k,3,0.368,,0
k,4,0.652,,0
k,5,0.708,,0
k,6,0.992,,0
k,7,1.048,,0
k,8,1.19,,0
k,9,0.17,,11.9   
k,10,0.312,,11.9
k,11,0.368,,11.9
k,12,0.652,,11.9
k,13,0.708,,11.9
k,14,0.992,,11.9
k,15,1.048,,11.9
k,16,1.19,,11.9
k,100,0,0,0
k,101,0,0,11.9
a,1,2,10,9   
a,3,11,10,2   
a,3,4,12,11
a,4,5,13,12
a,5,6,14,13
a,6,7,15,14
a,7,8,16,15

GPLOT   
vrotat,all,,,,,,100,101,360,2               !面旋转生成体         
   
vglue,all                                 

FLST,5,8,6,ORDE,8   
FITEM,5,1   
FITEM,5,3   
FITEM,5,5   
FITEM,5,7   
FITEM,5,-8  
FITEM,5,10  
FITEM,5,12  
FITEM,5,14  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       2, ,   1,       0    !!指定材料属性
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,6,6,ORDE,6   
FITEM,5,2   
FITEM,5,4   
FITEM,5,6   
FITEM,5,9   
FITEM,5,11  
FITEM,5,13  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       1, ,   1,       0     !!指定材料属性
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
/PNUM,MAT,1
/REPLOT                     
lsel,s,loc,z,0
lsel,a,loc,z,11.9
lsel,u,loc,y,0
LESIZE,all,,15, , , , , ,1
allsel

LSEL,S,LOC,Z,1,11
LESIZE,all,0.595, , , , , , ,1

esize,0.238
VSWEEP,all                          !!!!!生成网格
/solu
antype,static
asel,s,loc,z,0
lsel,r,ext
nsll,r,1
d,all,all                               !截面边缘加约束           
allsel
nsel,s,loc,z,11.9
*get,no_nodes,node,,count
f,all,fz,-1/no_nodes                    
allsel
OUTPR,,1
SOLVE
FINISH
/SOLU                                   
ANTYPE,BUCKLE                           
BUCOPT,LANB
MXPAND,1
SOLVE

1.solution是根据vm127做的，但是vm127的单元类型是beam3,我这里是solid45,不知道会对结果产生什么影响。
2.最后四句命令用的对不对？