请教高手一个关于非线性屈曲分析的问题

中升肖慧丽

请教大家一个非线性屈曲分析问题:'( :'(
模型边界条件：模型边界条件：两端固定，长度为17800毫米
有限元模型:长度为8900毫米，边界为一端固定、一端自由
单元类型beam188，
最后画变形-应力曲线时，变形取自由端端点，应力应该取哪一点的作用反力？以下是命令流，能否帮我看看什么毛病？
!单位：N,M
FINISH
/CLEAR
/UNITS,SI
/FILNAME,cellguide
！进入前处理器
/PREP7
!指定单元类型
ET,1,BEAM188
!定义材料属性
MP,EX,1,2.1E11
MP,PRXY,1,0.3
TB,BKIN,1
TBDATA,1,235e6,0
！定制梁截面
SECTYPE,1,beam,hrec,fangxing,3
SECOFFSET,CENT
SECDATA,0.189,0.2,0.02,0.02,0.012,0.012
!创建有限元模型
K,1,0,0,0              !创建关键点
K,2,0,8.9,0
K,3,0.25,4.45,0
LSTR,1,2                !生成线
LATT,1,,1,,3,,1                    !指定关键点4作为方向点
LESIZE,ALL,,,89,,1,,,1             !指定单元尺寸,划分178个网格
LMESH,1    !划分网格
FINISH


!求解
/SOLU
!施加固定约束
D,1,ALL,O
!施加单位集中载荷
F,2,FY,-1

ANTYPE,0
PSTRES,ON                  ！打开预应力
SOLVE
FINISH

!进行屈曲分析
/SOLU
ANTYPE,BUCKLE              ！指定分析类型为特征值屈曲分析
BUCOPT,LANB,1,             !激活模态扩展计算
MXPAND,1                   ！设置模态扩展选项
SOLVE
FINISH

!进入通用后处理器
/POST1
SET,FIRST                     ！显示一阶模态
NSORT,U,SUM                   ！储存节点位移
*GET,DMAX,SORT, ,MAX          !保存最大位移
*GET,PINLV,ACTIVE,,SET,FREQ   ！把FREQ的值赋给PINLV并保存
FINISH


!添加初始缺陷
/PREP7
UPGEOM, 0.2/DMAX,1,1,cellguide,RST
FINISH

! 非线性求解
/SOLU
D,1,ALL,0           !施加固定约束
F,2,FY,-PINLV*1.5   !施加载荷
ANTYPE,0          !静力分析
NLGEOM,1          !打开大变形
NSUBST,20         !设置步长为20步
ARCLEN,ON,5        !打开弧长法求解
ARCTRM,U,2,2,UX   
OUTRES,ALL,ALL
OUTPR,ALL,ALL
SOLVE

! 绘制荷载位移曲线
FINISH
/POST26
FILE
FILE,'cellguide','rst','.'
NUMVAR,200
NSOL,2,2,U,X,UX2             !定义变量UX2储存节点2的的位移UX
RFORCE,3,1,F,Y,FY3           !定义变量FY3储存节点3的的位移反作用力
PRVAR,2,3                    !列出变量
XVAR,2
PLVAR,3
/AXLAB,X,UX2(m)
/AXLAB,Y,FY3(N)
/REPLOT

qlp529

也刚开始学习这方面的
看了楼主的帖子启发很大