ansys做一个卷筒的屈曲分析,为什么一阶屈曲分析的变形那么大!都1米了都

ansys做一个卷筒的屈曲分析,为什么一阶屈曲分析的变形那么大!都1米了都
我要用ansys做一个卷筒的屈曲分析,可是做出来了以后,为什么一阶屈曲分析的变形那么大!都1米了都!
我定义了国际单位制,静力分析的结果是对的!为什么屈曲分析的结果就这么差呢?代码如下,电脑好一点运行也就几分钟,希望大虾解决下啊!
解答出来有用+50分啊!
FINISH
/CLEAR,START
/FILNAME,DRUM-SG,1
/TITLE,THE STABILITY OF DRUM
/COM,STRUCTURAL
------进入前处理器
/PREP7
/UNITS,SI
ET,1,SOLID45
MP,EX,1,2.06E11
MP,PRXY,1,0.3
------ 建立几何模型
CYLIND,0.300,0.333,2.000,
CYLIND,0.225,0.305,0.138,0.178,
CYLIND,0.225,0.305,1.910,1.950,
CYLIND,0.326,0.340,0.8772,0.8962,0,180
CYLIND,0.326,0.340,0.8772,0.8962,180,360
CYLIND,0.326,0.340,1.1038,1.1228,0,180
CYLIND,0.326,0.340,1.1038,1.1228,180,360
CYLIND,0.326,0.340,0.019,
VGEN,34,8,0.0258,
CYLIND,0.326,0.340,1.1296,1.1486,
VGEN,34,42,0.0258,
CYLIND,0.295,0.305,0.103,
VSEL,S,4
VSEL,A,5,76
CM,subvol,VOLU 所有部分接合
ALLSEL 全部选中
VSBV,1,subvol 布尔减运算
VADD,all 全部布尔加运算
-----------划分网格
WPCSYS,-1,0
WPOFFS,1,
KWPAVE,3
WPROT,0,90
VSBW,1 将圆柱分为两部分
VPLOT 体 显示
ESIZE,0.020,0
VSWEEP,2,下半部分扫略划分
ESIZE,0.020,0
VSWEEP,3,下半部分扫略划分
WPCSYS,1,0 恢复坐标系
ALLSEL,ALL
FINISH
-----------静态求解
/SOLU
-----------载荷与约束
CSYS,1 转化成极坐标
ASEL,S,LOC,X,0.225 选中极坐标半径X
ASEL,R,LOC,Z,0.138,0.178 选中极坐标高度Z
DA,ALL,ALL 全部约束
ASEL,S,LOC,Z,0 选中极坐标高度Z=0
DA,ALL,UX 全部约束
DA,ALL,UY
ALLSEL,ALL 全部恢复
ASEL,S,LOC,X,0.225 选中极坐标半径X
ASEL,R,LOC,Z,1.910,1.950 选中极坐标高度Z
DA,ALL,UX,UY 全部约束
ASEL,S,LOC,Z,2.000 选中极坐标高度Z
DA,ALL,ALL 全部约束
ALLSEL,ALL 全部恢复
ASEL,S,LOC,X,0.326 选中极坐标半径X
ASEL,R,LOC,Y,0,180 选中极坐标角度Y
ASEL,R,LOC,Z,0.8772,0.8962 选中极坐标高度Z
SFA,ALL,1,PRES,9670000
ALLSEL,ALL 全部恢复
ASEL,S,LOC,X,0.326 选中极坐标半径X
ASEL,R,LOC,Y,0,180 选中极坐标角度Y
ASEL,R,LOC,Z,1.1038,1.1228 选中极坐标高度Z
SFA,ALL,1,PRES,9670000
ALLSEL,ALL
-----------求解选项定义
ANTYPE,0
NLGEOM,0 忽略大变形
EQSLV,SPAR 设定稀疏方程求解
PSTRES,ON,预应力效果打开
SOLVE
FINISH
-----------进入线性动态求解
/SOLU 动态解
ANTYPE,BUCKLE 分析类型是特征屈曲分析
BUCOPT,LANB,3 设定block lanczos算法,模态提取方法与数量
MXPAND,3,0,0,YES 模态扩张方法与模态扩张数量
SOLVE

仔细看了看命令流,问题可能出在这句,SFA,ALL,1,PRES,9670000,施加载荷受力面积很小,压力（9670000）却相当大,变形能不大吗?改变一下压力载荷尝试一下.