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