钢-混凝土组合梁非线性分析ANSYS命令流

这是一个带托梁的钢――混凝土组合梁模型，实验时采用的是弯筋连接件，我用COMBIN39单元模拟。
/clear
/prep7
!朱聘儒的组合梁试验模拟
ET ,1 ,SOLID65
R,1
ET,2,SHELL43
R,2, 11.4!钢梁上下翼缘厚度
R,3, 7!钢梁腹板厚度
ET,3,pipe20
R,5,8, 3.99!受拉纵筋
R,6,8, 3.99!箍筋
et,4,combin39
keyopt,4,3,0
keyopt,4,3,2
R,7,-19.47,-35325,-1.288,-32499,-0.18,-21902
RMORE,,,0.18,21902,1.288,32499
RMORE,19.47,35325,,,,,
mp,ex,1,3e4!混凝土
mp,prxy,1,0.2
TB,KINH,1,1,10,
TBTEMP,0
TBPT,,0.0001356,4.068
TBPT,,0.000678,11.12
TBPT,,0.000949,13.211
TBPT,,0.001085,13.83
TBPT,,0.001356,14.3
TBPT,,0.001627,14.13
TBPT,,0.00203,13.48
TBPT,,0.00271,12.1
TBPT,,0.004068,9.63
TBPT,,0.0053,7.15
TB,CONCR,1,1,9,
TBTEMP,0
TBDATA,,0.5,0.9,21.5,2.15,,
TBDATA,,,,0.6,,,
mp,ex,2,2.1e5!钢板
mp,prxy,2,0.3
TB,BKIN,2,1,2
TBTEMP,0
TBDATA,,215E6,2.1E3
mp,ex,3,2.1e5!纵筋及箍筋
mp,prxy,3,0.3

k,1,,305,-285
k,2,,305,285
k,3,,335,285
k,4,,335,-285
k,5,2500,305,-285
k,6,2500,305,285
k,7,2500,335,285
k,8,2500,335,-285
k,9,,305,-100
k,10,,305,100
k,11,,335,100
k,12,,335,-100
k,13,2500,305,-100
k,14,2500,305,100
k,15,2500,335,100
k,16,2500,335,-100
l,1,2
l,2,3
l,3,4
l,4,1
lgen,9,1,4,1,300,
l,1,5
l,2,6
l,3,7
l,4,8
l,9,13
l,10,14
l,11,15
l,12,16
latt,3,6,3
lsel,s,length,,30
lesize,all,,,1
lsel,a,length,,570
lesize,all,,,5
lmesh,all
lsel,all
latt,3,5,3
lsel,s,length,,2500
lesize,all,,,25
lmesh,all
allsel,all
k,201,0,290,-300
k,202,0,290,-140
k,203,0,200,-50
k,204,0,200,50
k,205,0,290,140
k,206,0,290,300
k,207,0,350,300
k,208,0,350,140
k,209,0,350,-140
k,300,0,350,-300
a,201,202,209,300
a,202,205,208,209
a,205,206,207,208
a,202,203,204,205
asel,all
esize,80
mshape,0,2d
mshkey,1
amesh,all
asel,s,,,1,4,1
mat,1
REAL,1
TYPE,1
extopt,esize,25
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0
VEXT,1,4 ,1 ,2500
k,3001,0,200,-50
k,3002,0,200,0
k,3003,0,200,50
k,3004,2500,200,50
k,3005,2500,200,0
k,3006,2500,200,-50
k,3007,0,0,-50
k,3008,0,0,0
k,3009,0,0,50
k,4000,2500,0,50
k,4001,2500,0,0
k,4002,2500,0,-50
a,3001,3003,3004,3006
a,3007,3009,4000,4002
a,3002,3008,4001,3005
mshape,0,2d
mshkey,1
esize,100
asel,s,,,22,23
lsel,s,length,,100
lesize,all,,,4
aatt, 2,2,2
amesh,all
asel,all
asel,s,,,24
aatt,2,3,2
amesh,all
allsel,all
numcmp,all
type,4
real,7
*do,i,1,8
e,1097+3i,871+3i
e,1145+3i,919+3i
*enddo

/solu
nsel,s,loc,x,0
nsel,r,loc,y,0
d,all,ux
d,all,uy
allsel,all
nsel,s,loc,x,2500
dsym,symm!设置跨中面为对称面约束
allsel,all
nsel,s,loc,y,350
sf,all,pres,0.08!弹性极限状态时的荷载

antype,0
NCNV,2,
NSUBST,5
autots,on
lnsrch,1
outres,all,all
time,1
CNVTOL,F,,0.05,2
neqit,50
solve
finish
运行后，系统出现如下提示：
ERROR：Real const

ant set 7 undefined but refenced by element 963。
我定义的COMBIN39的实常数难道不对么？要是不对的话，不这样定义又该如何定义啊？
大侠们看一看，帮帮在下，我都快急死了！
TBDATA,,0.5,0.9,21.5,2.15,,

solid65的Material Data里的顺序是
1 Shear transfer coefficients for an open crack.
2 Shear transfer coefficients for a closed crack.
3 Uniaxial tensile cracking stress.
4 Uniaxial crushing stress (positive).
5 Biaxial crushing stress (positive).
6 Ambient hydrostatic stress state for use with constants 7 and 8.
7 Biaxial crushing stress (positive) under the ambient hydrostatic stress state (constant 6).
8 Uniaxial crushing stress (positive) under the ambient hydrostatic stress state (constant 6).
9 Stiffness multiplier for cracked tensile condition, used if KEYOPT(7) = 1 (defaults to 0.6).
那是不是应该把21.5和2.15调换位置啊？
楼上的amoisonic兄：
谢谢阁下，2.15和21.5的顺序确实搞错了。而且我也知道我的命令流哪里有问题了，是实常数没有按顺序编号，系统自动排序后就没有了7号实常数。我改过之后就没有大问题了，但是计算结果还是有些小问题。我正在努力寻找！
唉，发点牢骚，ANSYS学起来真是不顺手！一大堆的命令，也记不住，遇到一个就翻书找。真是难！