多层中空圆筒的屈曲分析(压杆稳定)可不可以用轴对称方法分析(vm128)
以下是我的命令 (运行到最后时,结果是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) 为什么没人回帖子?我的问题太没水平? 多层中空圆筒的屈曲分析(压杆稳定)可不可以用轴对称方法分析.
肯定是不太合适的,VM128的例子和你的有本质区别,例子上bar用plane单元来模拟是没问题的。你的圆筒用轴对称是不合适的,去找本屈曲方面的书看看基础知识。圆筒分长圆筒和短圆筒,屈曲模态是不相同的。对圆筒,不说你轴对称模型,就是1/4,1/2模型在没有确定把握的情况下也尽量不要用,因为屈曲模态不一定是对称的。所以屈曲分析在一般情况下最好用整体模型。 这是我做的体模型,网格很大,只是初步的计算。
/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.最后四句命令用的对不对?
页:
[1]