ansys钢筋混凝土锈蚀模拟过程中的问题？
ansys初学者，根据导师的论文想要做一个关于混凝土锈蚀的模拟，计算的时候出现如下提示，能力有限，自己检查了很多遍也没发现问题，希望各位大神帮忙解答，后附命令流，谢谢大家
finish
/clear
/filname,rcbeam
/prep7
et,1,solid65 !混凝土
et,2,solid45 !垫块
et,3,link8 !钢筋
et,4,combin39 !弹簧
keyopt,1,1,1 !不考虑形函数的附加项
keyopt,1,3,2 !无筋单元的压碎后行为
keyopt,1,7,1 !考虑开裂后的应力释放
!=====================
!参数设置
x=0.015 !锈蚀率
d1=28
d2=6
d3=12 !三种钢筋的直径
*if,x,gt,1.5e-2,then,p=1.346*2.71828^(-0.198*x)
*else,p=1
fck=36.45 !C40混凝土抗压强度标准值
ftk=2.77 !C40混凝土抗拉强度标准值
a=sqrt(fck)
b=sqrt(40.7)
!====================================
!混凝土本构模型定义
mp,ex,1,fck*0.19/0.0002
mp,prxy,1,0.2
tb,concr,1
tbdata,,0.5,0.95,ftk,-1
tb,miso,1,,11
tbpt,,0.0002,fck*0.19
tbpt,,0.0004,fck*0.36
tbpt,,0.0006,fck*0.51
tbpt,,0.0008,fck*0.64
tbpt,,0.001,fck*0.75
tbpt,,0.0012,fck*0.84
tbpt,,0.0014,fck*0.91
tbpt,,0.0016,fck*0.96
tbpt,,0.0018,fck*0.99
tbpt,,0.002,fck
tbpt,,0.0033,fck
!钢筋本构模型定义
mp,ex,2,2.1e5
mp,prxy,2,0.3
tb,biso,2
tbdata,,360,0
!垫块定义
mp,ex,3,2.1e6
mp,prxy,3,0.3
!=================================
!65单元实常数定义
!素混凝土
r,1
pi=acos(-1)
r,2,28*28*pi/4
r,6,12*12*pi/4
r,7,6*6*pi/4
k1=pi*d1*3200 !表面积*周长1
k2=pi*d2*440*2+190*2 !表面积*周长2
k3=pi*d3*12*3200 !表面积*周长3
!弹簧实常数定义
r,3,0,0,0.5E-2,a/b*k1*p*0.0054,1.0E-2,a/b*k1*p*0.014
rmore,1.5E-2,a/b*k1*p*0.0154,2.0E-2,a/b*k1*p*0.024,2.5E-2,a/b*k1*p*0.0254
rmore,3.0E-2,a/b*k1*p*0.034,3.5E-2,a/b*k1*p*0.0354,4.0E-2,a/b*k1*p*0.044
rmore,4.5E-2,a/b*k1*p*0.0454,5.0E-2,a/b*k1*p*0.054
r,4,0,0,0.5E-2,a/b*k2*p*0.0054,1.0E-2,a/b*k2*p*0.014
rmore,1.5E-2,a/b*k2*p*0.0154,2.0E-2,a/b*k2*p*0.024,2.5E-2,a/b*k2*p*0.0254
rmore,3.0E-2,a/b*k2*p*0.034,3.5E-2,a/b*k2*p*0.0354,4.0E-2,a/b*k2*p*0.044
rmore,4.5E-2,a/b*k2*p*0.0454,5.0E-2,a/b*k2*p*0.054
r,5,0,0,0.5E-2,a/b*k3*p*0.0054,1.0E-2,a/b*k3*p*0.014
rmore,1.5E-2,a/b*k3*p*0.0154,2.0E-2,a/b*k3*p*0.024,2.5E-2,a/b*k3*p*0.0254
rmore,3.0E-2,a/b*k3*p*0.034,3.5E-2,a/b*k3*p*0.0354,4.0E-2,a/b*k3*p*0.044
rmore,4.5E-2,a/b*k3*p*0.0454,5.0E-2,a/b*k3*p*0.054
!=====================
blc4,,,3600,250,500
wpcsys,-1
wpoff,,,-50
blc4,150,0,100,250,50
vgen,2,2,,,3200
wpcsys,-1
wpoff,0,0,500
blc4,1000,0,100,250,50
vgen,2,4,,,1600
wpcsys,-1
allsel
wpoff,,,30
vsbw,all
wpoff,,,440
vsbw,all
wprota,,,90
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpoff,,,200
vsbw,all
wpcsys,-1
wprota,0,-90
wpoff,0,0,30
vsbw,all
wpoff,0,0,190
vsbw,all
wpcsys,-1
wprota,,,90
wpoff,,,150
vsbw,all
wpoff,,,100
vsbw,all
wpoff,,,850
vsbw,all
wpoff,,,1600
vsbw,all
wpoff,,,3350-2700
vsbw,all
wpoff,,,100
vsbw,all
lsel,s,loc,y,29,221
lsel,r,loc,z,29,471
lsel,r,tan1,x
lsel,r,loc,x,1,3599
cm,stirrup,line !圈梁
allsel
wprota,-90
wpoff,-125
wprota,,,90
vsbw,all
allsel
wpcsys,-1
lsel,s,loc,z,30
lsel,r,tan1,y
lsel,r,loc,y,1,249
cm,long,line !受拉钢筋
lsel,s,loc,z,470
lsel,r,tan1,y
lsel,r,loc,y,1,249
lsel,u,loc,y,124,126
cm,hanger,line !构造钢筋
vsel,s,loc,z,499,1000
vsel,a,loc,z,-100,1
allsel
vglue,all
!===========================================
lsel,s,length,,30
lesize,all,,,6
lsel,s,length,,50
lesize,all,,,10
lsel,s,length,,100
lesize,all,,,20
lsel,s,length,,95
lesize,all,,,95/5
lsel,s,length,,440
lesize,all,,,440/5
lsel,s,length,,200
lesize,all,,,200/5
vsel,s,loc,z,0,500
vatt,1,,1
mshkey,1
esize,5
vmesh,all
allsel,all
vsel,u,mat,,1
vatt,3,,2
esize,5
vmesh,all
k,1001,200,30,30
k,1017,3400,30,30
kfill,1001,1017
k,1021,200,220,30
k,1037,3400,220,30
kfill,1021,1037
k,1041,200,125,30
k,1057,3400,125,30
kfill,1041,1057
l,1001,1017
l,1021,1037
l,1041,1057
k,2001,200,30,470
k,2017,3400,30,470
kfill,2001,2017
k,2021,200,220,470
k,2037,3400,220,470
kfill,2021,2037
l,2001,2017
l,2021,2037
*do,i,1001,1017
l,i,i+20
*enddo
*do,i,2001,2017
l,i,i+20
*enddo
*do,i,1001,1017
l,i,i+1000
*enddo
*do,i,1021,1037
l,i,i+1000
*enddo
lsel,s,length,,3200
lsel,a,loc,z,0,50
latt,2,2,3
lesize,all,,50
lmesh,all
!================================
lsel,s,loc,z,30
lsel,r,tan1,y
lsel,r,loc,y,1,249
cpintf,ux,0.01 !创建节点的耦合自由度UX
*get,cpmax,cp,,max !获取藕合组总数
type,4
real,3
*do,i,1,cpmax
*get,n1,cp,i,term,1,node !获取第i组的第一个节点号
*get,n2,cp,i,term,2,node !获取第i组的第二个节点号
e,n1,n2
*enddo
cpintf,uy
!===============================================
lsel,s,loc,y,29,221
lsel,r,loc,z,29,471
lsel,r,tan1,x
lsel,r,loc,x,1,3599
cpintf,ux,0.01 !创建节点的耦合自由度UX
*get,cpmax,cp,,max !获取藕合组总数
type,4
real,4
*do,i,1,cpmax
*get,n1,cp,i,term,1,node !获取第i组的第一个节点号
*get,n2,cp,i,term,2,node !获取第i组的第二个节点号
e,n1,n2
*enddo
cpintf,uy
!===============================================
lsel,s,loc,z,470
lsel,r,tan1,y
lsel,r,loc,y,1,249
lsel,u,loc,y,124,126
cpintf,ux,0.01 !创建节点的耦合自由度UX
*get,cpmax,cp,,max !获取藕合组总数
type,4
real,5
*do,i,1,cpmax
*get,n1,cp,i,term,1,node !获取第i组的第一个节点号
*get,n2,cp,i,term,2,node !获取第i组的第二个节点号
e,n1,n2
*enddo
cpintf,uy
!=======================
allsel,all
!约束
asel,s,loc,z,-50
da,all,all,0
!====================
asel,s,loc,z,550
sfa,all,1,pres,200e3/200/100
allsel,all
save
!=====================
/solu
antype,0
outres,all,all
autots,on
nsubst,500,,50
neqit,50
cnvtol,u,,0.05
cnvtol,f,,0.05
solve
张小 