采用梁单元，分析一个弹性地基梁的问题，ANSYS提示约束不满足要示，那位大侠能帮求...

xieyangbo

/com,structural
/title,beam analysis
/prep7
et,1,beam188
keyopt,1,1,1
et,2,COMBIN14
!keyopt,2,1,1
keyopt,2,2,2
mp,ex,1,3e10
mp,prxy,1,0.15
mp,dens,1,2500
sectype,1,beam,rect,,0
SECOFFSET, CENT
secdata,1.5,1,2,2
!secnum,1
!mp,ex,2,3e10
!mp,prxy,2,0.15
k,1,,,,
k,2,,7,,,

l,1,2

latt,1,,1,,,,1
lesize,all,1
lmesh,all
esel,all,

cm,beam,elem

*do,i,1,5
n,100+i,,i-1,0
n,200+i*2,-1,i-1,
*enddo

*do,i,1,5
type,2
!mat,2
real,10
esys,0           !!!!!
en,500+i,100+i,200+i*2
*enddo
numcmp,all

nummrg,all,   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!有问题

esel,all
!cmsel,u,beam
esel,s,type,,2
cm,spring,elem
cmsel,s,spring
nsle,s
cmsel,s,beam
nsle,u
d,all,all,,,,,

h=1         !!!!定义系数
b=1         !!!!定义系数
m1=50e6     !!!!定义系数

cmsel,s,spring
*DO,i,0,4
NSEL,S,LOC,y,1*i                   !从当前集选择一组节点子集
ESLN,S                            !从已选集中选择
ESEL,R,TYPE,,2                   !选择单元类型号为2的单元
R,i+1,m1*(4-i)*b*h,,,               !单位面积内受的力随深度增加而增加   !!!!!!!!!!!!!!!!!!!!!!!!!!!!采用m法计算         
EMODIF, ALL, REAL, i+1,             !对已存在单元进行修改
*enddo

secplot,1   !!!!!显示截面参数

!acel,,9.8
*set,qk,384.4  
allsel,all
cmsel,s,beam
nsle,all
f,all,fx,-qk
!sfcum,pres,,
!sfbeam,all,,pres,-qk
allsel,all

cmsel,s,beam
nsle,s
nsel,r,loc,y,0

d,all,uy,0,,,,
d,all,ux,0,,,,
d,all,uy,0,,,,
allsel,all

!cmsel,s,beam

!r,1,1*1.5,1*1.5**3/12
!emodif,all,real,1

allsel,all
/solu
solve