ansys鋼筋混凝土銹蝕模擬過程中的問題？
ansys初學(xué)者，根據(jù)導(dǎo)師的論文想要做一個關(guān)于混凝土銹蝕的模擬，計算的時候出現(xiàn)如下提示，能力有限，自己檢查了很多遍也沒發(fā)現(xiàn)問題，希望各位大神幫忙解答，后附命令流，謝謝大家
finish
/clear
/filname,rcbeam
/prep7
et,1,solid65 !混凝土
et,2,solid45 !墊塊
et,3,link8 !鋼筋
et,4,combin39 !彈簧
keyopt,1,1,1 !不考慮形函數(shù)的附加項
keyopt,1,3,2 !無筋單元的壓碎后行為
keyopt,1,7,1 !考慮開裂后的應(yīng)力釋放
!=====================
!參數(shù)設(shè)置
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混凝土抗壓強度標(biāo)準(zhǔn)值
ftk=2.77 !C40混凝土抗拉強度標(biāo)準(zhǔn)值
a=sqrt(fck)
b=sqrt(40.7)
!====================================
!混凝土本構(gòu)模型定義
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
!鋼筋本構(gòu)模型定義
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單元實常數(shù)定義
!素混凝土
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
!彈簧實常數(shù)定義
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 !構(gòu)造鋼筋
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 !創(chuàng)建節(jié)點的耦合自由度UX
*get,cpmax,cp,,max !獲取藕合組總數(shù)
type,4
real,3
*do,i,1,cpmax
*get,n1,cp,i,term,1,node !獲取第i組的第一個節(jié)點號
*get,n2,cp,i,term,2,node !獲取第i組的第二個節(jié)點號
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 !創(chuàng)建節(jié)點的耦合自由度UX
*get,cpmax,cp,,max !獲取藕合組總數(shù)
type,4
real,4
*do,i,1,cpmax
*get,n1,cp,i,term,1,node !獲取第i組的第一個節(jié)點號
*get,n2,cp,i,term,2,node !獲取第i組的第二個節(jié)點號
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 !創(chuàng)建節(jié)點的耦合自由度UX
*get,cpmax,cp,,max !獲取藕合組總數(shù)
type,4
real,5
*do,i,1,cpmax
*get,n1,cp,i,term,1,node !獲取第i組的第一個節(jié)點號
*get,n2,cp,i,term,2,node !獲取第i組的第二個節(jié)點號
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
張小 
