最近分别建了素混凝土和FRP包裹混凝土的模型，计算后有点问题，还望高手指点下

jslyzz

素混凝土轴压
/prep7
et,1,65,,,,,,,1
mp,ex,1,13585$mp,prxy,1,0.2
fc=14.3$ft=1.43
tb,concr,1$tbdata,,0.5,0.95,ft,-1
tb,miso,1,,11
tbpt,,0.0002,fc*0.19$tbpt,,0.0004,fc*0.36$tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64$tbpt,,0.0010,fc*0.75$tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91$tbpt,,0.0016,fc*0.96$tbpt,,0.0018,fc*0.99
tbpt,,0.002,fc$tbpt,,0.0033,fc*0.85
cyl4,,,150,0,,90,450
elemsiz=50
vatt,1,,1$mshkey,1$esize,elemsiz$vmesh,all
csys,0
asel,s,loc,x,0
da,all,symm
asel,s,loc,y,0
da,all,symm
nsel,s,loc,z,0
d,all,all
nsel,s,loc,z,450
d,all,ux,,,,,uy
cp,1,uz,all
D,38,uz,-1
nsel,all
/solu
antype,0$nsubst,50$outres,all,all$autos,off
neqit,50
cnvtol,f,,0.015,0
nsel,all
solve
/POST26
NSOL,2,38,U,z,Uz
RFORCE,3,38,F,z,Fz
ADD,4,2,,,STRAIN,,,-1/450
ADD,5,3,,,STRESS,,,-1/17662.5
/AXLAB,X,STRAIN
/AXLAB,Y,STRESS
XVAR,4
PLVAR,5


碳纤维约束圆柱混凝土轴压
/prep7
et,1,65,,,,,,,1$et,2,41,2
mp,ex,1,13585$mp,prxy,1,0.2
fc=14.3$ft=1.43
tb,concr,1$tbdata,,0.5,0.95,ft,-1
tb,miso,1,,11
tbpt,,0.0002,fc*0.19$tbpt,,0.0004,fc*0.36$tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64$tbpt,,0.0010,fc*0.75$tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91$tbpt,,0.0016,fc*0.96$tbpt,,0.0018,fc*0.99
tbpt,,0.002,fc$tbpt,,0.0033,fc*0.85
mp,ex,2,400e3$mp,prxy,2,3e-4
R,1,2
cyl4,,,151,0,,90,450
vatt,1,,1
csys,1$asel,s,loc,x,151$aatt,2,1,2
esize,50
vsel,s,,,1$vmesh,all
mshkey,0$asel,s,,,3$amesh,all
allsel,all
nummrg,all$numcmp,all
allsel,all
csys,0
nsel,s,loc,z,0
d,all,all
asel,s,loc,y,0
da,all,symm
asel,s,loc,x,0
da,all,symm
nsel,s,loc,z,450
d,all,ux,,,,,uy
cp,next,uz,all
D,38,uz,-1.2
/solu
ANTYPE,0
autos,off
nsubst,200
neqit,100
CNVTOL,f, ,0.015,2
OUTRES,ALL,1
allsel,all
solve
/POST26
NSOL,2,38,U,z,Uz
RFORCE,3,38,F,z,Fz
ADD,4,2,,,STRAIN,,,-1/450
ADD,5,3,,,STRESS,,,-1/17662.5
/AXLAB,X,STRAIN
/AXLAB,Y,STRESS
XVAR,4
PLVAR,5
FRP包裹混凝土的那个算下来FRP的应力才400Mpa跟以前看的文献中最终FRP拉坏不符，而且极限承载力较素混凝土提高不大  还望高手指点迷津  小弟不胜感激。