钢筋混凝土T梁模型改来改去还是那个错误，请高手指点修改

The value of UY at node 3585 is 11757596.2.  It is greater than the  current limit of 1000000.  This generally indicates rigid body motion   as a result of an
unconstrained model.  Verify that your model is   properly constrained.我把箍筋去掉还是出现类似这种错误，不知道模型哪里出问题了，请高手指点修改，下面是命令流


/clear   
/prep7  
et,1,solid45
et,2,link8  
et,3,plane42
mp,ex,1,3.0E4
mp,prxy,1,0.2   


mp,ex,2,2.0e5   
mp,prxy,2,0.3   
tb,biso,2   
tbdata,,300,0   
*SET,pi,acos(-1)
R,1,0.25*pi*32*32   
R,2,0.25*pi*22*22   
R,3,0.25*pi*8*8  

K,1,90  
K,2,90,1300
K,3,1100,1350   
K,4,1100,1500   
K,5,-1100,1500  
K,6,-1100,1350  
K,7,-90,1300
K,8,-90
A,1,2,3,4,5,6,7,8  

wprota,,,90
wpoffs,,,90
asbw,all
wpoffs,,,-180  
asbw,all
   
wpcsys,-1     
wprota,,-90
wpoffs,,,1300   
asbw,all

wpcsys,-1
wprota,,-90
wpoff,,,48
asbw,all

wpoff,,,72
asbw,all

wpoff,,,72
asbw,all

wpoff,,,72
asbw,all

wpoff,,,72
asbw,all

wpoff,,,1114
asbw,all

wpcsys,-1
wpoff,41
wprota,,,90
asbw,all


wpcsys,-1
wpoff,-41
wprota,,,90
asbw,all

esize,100
mshape,0,2d
mshkey,1
type,3
amesh,all
numcmp,all

*get,amax,area,,num,max
*do,i,1,amax
voffs,i,20000
*enddo
vglue,all
vsweep,all


lsel,s,loc,x,41
lsel,a,loc,x,-41
lsel,r,tan1,y
latt,2,2,2,
lesize,all,100,
lmesh,all

asel,s,loc,z,0
lsla,s
lsel,u,loc,x,-90
lsel,u,loc,x,90
cm,gujin,line
latt,2,3,2
lesize,all,100,
lmesh,all
lgen,2,all,,,,,250
cmsel,u,gujin
cm,djgj,line
lsel,a,,,gujin
lgen,65,djgj,,,,,300
lgen,2,gujin,,,,,19750
lgen,2,gujin,,,,,20000
allsel,all
nummrg,KP
numcmp,all
save

lsel,s,loc,z,249.9,250.1
lsel,a,loc,z,19749.9,19750.1
lsel,r,loc,y,0
esll,s
nsle,s
d,all,uy,,,,,uz
nod1=node(90,0,250)
d,nod1,ux
allsel,all
nod2=node(90,0,19750)
d,nod1,ux
allsel,all

!加载
nod3=node(0,1500,1600)
f,nod3,fy,-140

nod4=node(0,1500,3000)
f,nod4,fy,-140

nod5=node(0,1500,10000)
f,nod5,fy,-120

nod6=node(0,1500,11400)
f,nod6,fy,-120

nod7=node(0,1500,14400)
f,nod7,fy,-300

/solu
antype,0
nsubst,20
outres,all,all
autos,on
NEQIT,10
CNVTOL,U,,0.015
solve

/POST1
SET,LAST
PLDISP,1

根据以往的经验，是两个方面的原因。一个可能是边界约束不够，另一个是加载步设置问题。

你的模型虚线表示钢筋吧，你的约束怎么加到钢筋上了？应该加到板单元的节点上啊。

nod2=node(90,0,19750)
d,nod1,ux
这个地方有问题吧，应该约束nod2，uy