Устойчивость в нелинейной постановке

Добрый день!

Путаюсь определить устойчивость в нелинейной постановке (Nonlineal buckling), однако результат расчета удручает. Подскажите, пожалуйста, где в коде проблемы. Спасибо.

finish
/clear,start
/prep7
!
! Beam-web, loaded with transverse load
! (размеры в см)
! geometrically nonlineal
!
FNL=0 !1 ! Учет физической нелинейности (0 - упругая работа, 1 - пластика)
PIN=0 ! simply suported
QMOD=1 !2 !1 !Loaded with one force
QBM=40 ! Value of load , N
BU=1 ! 1 ! Buckling
!--------------------
HBM=75 !64 ! Depth of beam
!---------------------
TWB=1 ! Thickness of web, mm
ELSZ=10 ! Element size, mm
LBM=800 ! Length of beam, mm
HW=HBM! Depth of beam

! Material
MP,EX,1,2.1e5
MP,NUXY,1,0.3
ET,1,SHELL181
!----------------
R,1,TWB !  Thickness
!
k,1,0,0,0
k,2,LBM,0,0
k,3,LBM,0,HW
k,4,(LBM/2),0,HW
k,5,0,0,HW
!---------------------
a,1,2,3,4,5
asel,all
!----------------------------------------------------------------
asel,s,loc,z,0,HW
aatt,1,1
!-------------
!*go,:end
!-------------
esize,ELSZ
amesh,all
!=================
*if,PIN,eq,0,then
nsel,s,loc,x,0
nsel,r,loc,z,0
d,all,ux
d,all,uy
d,all,uz
!-------------
nsel,s,loc,x,LBM
nsel,r,loc,z,0
d,all,uy
d,all,uz
!-------------
nsel,all
*endif
!-----условие прямолинейности вертикальных кромок----------
nsel,s,loc,y,0
nsel,r,loc,x,0
d,all,uy
!---------------
nsel,s,loc,y,0
nsel,r,loc,x,LBM
d,all,uy
!---------------
nsel,all
!-----условие прямолинейности горизонтальных кромок----------
nsel,s,loc,y,0
nsel,r,loc,z,HW
d,all,uy
nsel,s,loc,y,0
nsel,r,loc,z,0
d,all,uy
nsel,all
!--------конец условия прямолинейности-------
*if,QMOD,eq,1,then
fk,4,fz,-QBM     
!------------------------
f,NODE(LBM/4,0,HW/2),fy,(QBM*0.1)
f,NODE(LBM/2,0,HW/2),fy,(QBM*0.1)
!------------------------
*endif
!-------------------------
!*go,:end
!-----------------
! Решение задачи
allsel
!===============
! Решение задачи
!============
/SOLU
!NLGEOM,ON
!ANTYPE,0
!-----------------
PSTRES,ON   
solve
finish
*if,BU,eq,1,then
/SOLU
antype,1
BUCOPT,lanb,5,0,0  
MXPAND,5,0,0,0  
!-----------------
SOLVE   
FINISH
*endif
/POST1  
SET,FIRST   
PLDISP,0
:end