有没有做钢管混凝土柱的小伙伴啊？？？

/BATCH
/input,menust,tmp,''
WPSTYLE,,,,,,,,0
/PREP7
ET,1,SOLID185
ET,2,SOLID65
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2e5
MPDATA,PRXY,1,,0.3
TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,235,,,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,2.8e4
MPDATA,PRXY,2,,0.2
TB,CONC,2,1,9,
TBTEMP,0
TBDATA,,0.35,0.9,2.35,-1,,
TBDATA,,,,,,,
TBDE,MISO,2
TB,MISO,2,1,21
TBTEMP
TBPT,,0.0001,2.8
TBPT,,0.0003,7.2049
TBPT,,0.0004,9.4398
TBPT,,0.0005,1.1589E+001
TBPT,,0.0006,1.365E+001
TBPT,,0.0007,1.5621E+001
TBPT,,0.0008,1.7501E+001
TBPT,,0.0009,1.9287E+001
TBPT,,0.001,2.0976E+001
TBPT,,0.0011,2.2566E+001
TBPT,,0.0012,2.4053E+001
TBPT,,0.0013,2.5434E+001
TBPT,,0.0014,2.6704E+001
TBPT,,0.0015,2.7859E+001
TBPT,,0.0016,2.8891E+001
TBPT,,0.0017,2.9795E+001
TBPT,,0.0018,3.0559E+001
TBPT,,0.0019,3.1168E+001
TBPT,,0.002,3.1601E+001
TBPT,,0.0021,3.18E+001
TBPT,,0.0031,3.18E+001
TBDE,MISO,2
TB,MISO,2,2,14,0 !多线性各向同性强化模型
TBTEMP,0
TBPT,,0.0001,2.8 !应力应变数据表
TBPT,,0.0003,5.5782
TBPT,,0.0006,10.5336
TBPT,,0.0009,14.8662
TBPT,,0.0012,18.576
TBPT,,0.0016,22.5536
TBPT,,0.002,25.4239
TBPT,,0.0033,27.739
TBPT,,0.005,22.267
TBPT,,0.01,10.7514
TBPT,,0.015,7.911
TBPT,,0.02,6.658
TBPT,,0.025,5.915
TBPT,,0.03,5.40255
TBDE,MISO,2,,,
TB,MISO,2,1,14,0
TBTEMP,0
TBPT,,0.0001,2.8
TBPT,,0.0003,5.5782
TBPT,,0.0006,10.534
TBPT,,0.0009,14.866
TBPT,,0.0012,18.576
TBPT,,0.0016,22.554
TBPT,,0.002,25.424
TBPT,,0.0033,27.739
TBPT,,0.005,22.267
TBPT,,0.01,10.751
TBPT,,0.015,7.911
TBPT,,0.02,6.658
TBPT,,0.025,5.915
TBPT,,0.03,5.4025
TBDE,MISO,2
TB,MISO,2,1,8,0 !多线性各向同性强化模型
TBTEMP,0
TBPT,,0.0001,2.8
TBPT,,0.0003,2.664
TBPT,,0.0006,4.896
TBPT,,0.0009,6.696
TBPT,,0.0012,8.064
TBPT,,0.0016,9.216
TBPT,,0.002,9.6
TBPT,,0.033,9.6
TBDE,MISO,2
TBDE,MISO,2
TB,MISO,2,1,21
TBTEMP
TBPT,,0.0001,2.8
TBPT,,0.0003,7.2049
TBPT,,0.0004,9.4398
TBPT,,0.0005,1.1589E+001
TBPT,,0.0006,1.365E+001
TBPT,,0.0007,1.5621E+001
TBPT,,0.0008,1.7501E+001
TBPT,,0.0009,1.9287E+001
TBPT,,0.001,2.0976E+001
TBPT,,0.0011,2.2566E+001
TBPT,,0.0012,2.4053E+001
TBPT,,0.0013,2.5434E+001
TBPT,,0.0014,2.6704E+001
TBPT,,0.0015,2.7859E+001
TBPT,,0.0016,2.8891E+001
TBPT,,0.0017,2.9795E+001
TBPT,,0.0018,3.0559E+001
TBPT,,0.0019,3.1168E+001
TBPT,,0.002,3.1601E+001
TBPT,,0.0021,3.18E+001
TBPT,,0.0031,3.18E+001
BLOCK,0,150,0,150,0,1500,
BLOCK,3,147,3,147,0,1500,
vsbv,1,2
BLOCK,3,147,3,147,0,1500,
vglue,all
allsel,all
lsel,s,lengeh,,,150
lsel,a,length,,,147
ALLSEL,ALL
GPLOT
lsel,s,length,,150
lsel,a,length,,147
lsel,a,length,,,147
lesize,all,,,60
ALLSEL,ALL
GPLOT
lsel,s,length,,1500
lesize,all,,,10
allsel,all
mtplot
GPLOT
VSEL,S,,,1
type,2
mat,2
vatt,2,,2
vsweep,all
VCLEAR, 1
GPLOT
lsel,s,length,,1500
lesize,all,,,30
GPLOT
FLST,2,8,4,ORDE,4
FITEM,2,9
FITEM,2,-12
FITEM,2,33
FITEM,2,-36
LCLEAR,P51X
lsel,s,length,,1500
FLST,5,8,4,ORDE,4
FITEM,5,9
FITEM,5,-12
FITEM,5,33
FITEM,5,-36
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,60, , , , ,1
ALLSEL,ALL
GPLOT
VSEL,S,,,1
type,2
mat,2
vatt,2,,2
vsweep,all
ALLSEL,ALL
GPLOT
vsel,,2
vsel,s,,,2
type,1
mat,1
vatt,1,,1
vsweep,all
ALLSEL,ALL
GPLOT
FLST,5,2,5,ORDE,2
FITEM,5,1
FITEM,5,17
ASEL,S, , ,P51X
NSLA,S,1
FINISH
/SOL
FLST,2,4049,1,ORDE,4
FITEM,2,7142
FITEM,2,-10862
FITEM,2,227070
FITEM,2,-227397
/GO
D,P51X, , , , , ,ALL, , , , ,
ALLSEL,ALL
GPLOT
FLST,2,2,5,ORDE,2
FITEM,2,2
FITEM,2,18
/GO
SFA,P51X,1,PRES,15
NSUBST,100,10,1
OUTRES,ERASE
OUTRES,ALL,ALL
TIME,100
NSUBST,10,100,1
SOLVE
FINISH