我在用ansys建一个钢筋混凝土梁碳纤维加固，但是加载以后就没变化，因为初学者，不会，但又很急很急。如果有的话就谁好心发发我啊。这个是我建的，好像完全不行的。梁高宽200*200，梁的长为1500
!碳纤维布加固钢筋混凝土梁有限元分析
/com,structural
/prep7
!定义单元类型
et,1,link8
et,2,solid65
et,3,plane42
et,4,shell41

!定义实常数
r,1,28.3!箍筋和受压钢筋
r,2,78.5!受拉钢筋
r,3!混凝土
r,4,0.12!碳纤维布

!定义材料属性
mp,ex,1,2e5!箍筋和受压钢筋
mp,prxy,1,0.25
tb,bkin,1,1,2,1
tbdata,,200

mp,ex,2,2e5!受拉钢筋
mp,prxy,2,0.27
tb,bkin,2,1,2,1
tbdata,,300

mp,ex,3,28000!混凝土
mp,prxy,3,0.2
tb,conc,3,1,9
tbdata,,0.5,1,1.78,-1
mp,ex,4,225000!碳纤维布
mp,prxy,4,0
tb,bkin,4,1,2,1
tbdata,,4300

/pnum,node,1
/pnum,elem,1

!产生所有的节点
n,1
n,9,200
fill,1,9
ngen,9,9,1,9,1,,25
ngen,301,1000,1,81,1,,,-5
/view,1,1,1,1
nplot

!箍筋和受压钢筋单元建立
type,1
real,1
mat,1
*do,ii,20011,20016,1
e,ii,ii+1
*enddo
*do,ii,20017,20062,9
e,ii,ii+9
*enddo
*do,ii,20065,20070,1
e,ii,ii+1
*enddo
*do,ii,20011,20056,9
e,ii,ii+9
*enddo
*do,ii,36011,36016,1
e,ii,ii+1
*enddo
*do,ii,36017,36062,9
e,ii,ii+9
*enddo
*do,ii,36065,36070,1
e,ii,ii+1
*enddo
*do,ii,36011,36056,9
e,ii,ii+9
*enddo
*do,ii,52011,52016,1
e,ii,ii+1
*enddo
*do,ii,52017,52062,9
e,ii,ii+9
*enddo
*do,ii,52065,52070,1
e,ii,ii+1
*enddo
*do,ii,52011,52056,9
e,ii,ii+9
*enddo
*do,ii,68011,68016,1
e,ii,ii+1
*enddo
*do,ii,68017,68062,9
e,ii,ii+9
*enddo
*do,ii,68065,68070,1
e,ii,ii+1
*enddo
*do,ii,68011,68056,9
e,ii,ii+9
*enddo
*do,ii,84011,84016,1
e,ii,ii+1
*enddo
*do,ii,84017,84062,9
e,ii,ii+9
*enddo
*do,ii,84065,84070,1
e,ii,ii+1
*enddo
*do,ii,84011,84056,9
e,ii,ii+9
*enddo
*do,ii,100011,100016,1
e,ii,ii+1
*enddo
*do,ii,100017,100062,9
e,ii,ii+9
*enddo
*do,ii,100065,100070,1
e,ii,ii+1
*enddo
*do,ii,100011,100056,9
e,ii,ii+9
*enddo
*do,ii,125011,125016,1
e,ii,ii+1
*enddo
*do,ii,125017,125062,9
e,ii,ii+9
*enddo
*do,ii,125065,125070,1
e,ii,ii+1
*enddo
*do,ii,125011,125056,9
e,ii,ii+9
*enddo
*do,ii,150011,150016,1
e,ii,ii+1
*enddo
*do,ii,150017,150062,9
e,ii,ii+9
*enddo
*do,ii,150065,150070,1
e,ii,ii+1
*enddo
*do,ii,150011,150056,9
e,ii,ii+9
*enddo
*do,ii,175011,175016,1
e,ii,ii+1
*enddo
*do,ii,175017,175062,9
e,ii,ii+9
*enddo
*do,ii,175065,175070,1
e,ii,ii+1
*enddo
*do,ii,175011,175056,9
e,ii,ii+9
*enddo
*do,ii,200011,200016,1
e,ii,ii+1
*enddo
*do,ii,200017,200062,9
e,ii,ii+9
*enddo
*do,ii,200065,200070,1
e,ii,ii+1
*enddo
*do,ii,200011,200056,9
e,ii,ii+9
*enddo
*do,ii,216011,216016,1
e,ii,ii+1
*enddo
*do,ii,216017,216062,9
e,ii,ii+9
*enddo
*do,ii,216065,216070,1
e,ii,ii+1
*enddo
*do,ii,216011,216056,9
e,ii,ii+9
*enddo
*do,ii,232011,232016,1
e,ii,ii+1
*enddo
*do,ii,232017,232062,9
e,ii,ii+9
*enddo
*do,ii,232065,232070,1
e,ii,ii+1
*enddo
*do,ii,232011,232056,9
e,ii,ii+9
*enddo
*do,ii,248011,248016,1
e,ii,ii+1
*enddo
*do,ii,248017,248062,9
e,ii,ii+9
*enddo
*do,ii,248065,248070,1
e,ii,ii+1
*enddo
*do,ii,248011,248056,9
e,ii,ii+9
*enddo
*do,ii,264011,264016,1
e,ii,ii+1
*enddo
*do,ii,264017,264062,9
e,ii,ii+9
*enddo
*do,ii,264065,264070,1
e,ii,ii+1
*enddo
*do,ii,264011,264056,9
e,ii,ii+9
*enddo
*do,ii,280011,280016,1
e,ii,ii+1
*enddo
*do,ii,280017,280062,9
e,ii,ii+9
*enddo
*do,ii,280065,280070,1
e,ii,ii+1
*enddo
*do,ii,280011,280056,9
e,ii,ii+9
*enddo


!受压钢筋
*do,ii,65,299065,1000
e,ii,ii+1000
*enddo
*do,ii,71,299071,1000
e,ii,ii+1000
*enddo

!纵向受拉钢筋单元建立
type,1
real,2
mat,2
*do,ii,11,299011,1000
e,ii,ii+1000
*enddo
*do,ii,17,299017,1000
e,ii,ii+1000
*enddo
!碳纤维布单元的建立
type,4
real,4
mat,4
*do,ii,1,8,1
*do,jj,ii,299000+ii,1000
e,jj,jj+1,jj+1000
e,jj+1,jj+1000,jj+1000+1
*enddo
*enddo

/view,1,1,1,1
/pnum,elem,0
/pnum,node,0
/eshape,1
eplot

!混凝土单元的建立
k,1
k,2,200
k,3,200,200
k,4,,200
k,9,100,200,-500
k,10,100,200,-1000
a,1,2,3,4

lsel,s,loc,y,0
lsel,a,loc,y,200
lesize,all,,,5
lsel,all
lsel,s,loc,x,0
lsel,a,loc,x,200
lesize,all,,,8
type,3
amesh,all

type,2
real,3
mat,3
extopt,esize,32
extopt,aclear,1
vext,all,,,,,-1500

/pnum,mat,1
/pnum,node,0
/pnum,elem,0
/view,1,1,1,1
eplot
allsel
nummrg,all
numcmp,all
eplot

/solu
nsel,s,loc,y,0
nsel,r,loc,z,-100
d,all,uy,ux

nsel,s,loc,z,-1400
nsel,r,loc,y,0
d,all,uy

allsel
gplot


fk,9,fy,-200000
fk,10,fy,-200000
allsel
nlgeom,on
nsubst,200
outres,all,all
neqit,50
pred,on
!ncnv,2
cnvtol,u,,0.015
solve