/filname,C3

/title,C3

/PREP7  

! FRP 厚度, FRP thickness

Fthick=5

! 混凝土弹模, Concrete Elastic Modulus

Ec=33109.7651

! 混凝土泊松比, Concrete Poisson ratio

EMUc=0.2          

! 混凝土抗压抗拉强度

! Concrete Compressive and Tensile Strength

fc=45

ft=2.00

! FRP弹模和泊松比

! Elastic Modulue and Poisson ratio for FRP

Ef=22000

EMUf=0.2

!型钢弹性模量

Es=2.06e5

!型钢泊松比

Esmuf=0.3

!型钢抗压抗拉强度

fy=235         !常温下屈服应力

exx=2.06e5             !常温下杨氏模量

! Height and Radius of Tube

! 圆柱高度和半径

H=1400

Rf=95

! 选择混凝土单元, Concrete element

ET,1,SOLID65

keyopt,1,5,1         !混凝土的线性解输出控制：给出每个积分点的解

keyopt,1,6,3         !混凝土的线非性解输出控制：给出每个积分点的解

keyopt,1,7,1         !考虑应力松弛

! 选择FRP单元

ET,2,solid45   

! 输入混凝土材料

MPTEMP,,,,,,,,  

MPTEMP,1,0

MPDATA,EX,1,,Ec

MPDATA,PRXY,1,,EMUc

! 输入应力应变曲线, Concrete stress-stain curve

TB,kinh,1,1,

TBTEMP,0

tbpt,,164e-6,5.43

tbpt,,493e-6,15.61

tbpt,,821e-6,24.84

tbpt,,1149e-6,33.15

tbpt,,1478e-6,40.54

tbpt,,1806e-6,46.99

tbpt,,2134e-6,52.5

tbpt,,2462e-6,57.09

tbpt,,2790e-6,60.75

tbpt,,3118e-6,63.49

tbpt,,3446e-6,65.35

tbpt,,3774e-6,66.7

tbpt,,4102e-6,68.06

tbpt,,4430e-6,69.42

tbpt,,4758e-6,70.78

! 输入破坏曲面, Failure envelop surface

TB,CONC,1,1,9,

TBTEMP,0

TBDATA,,0.5,1,ft,-1,,

TBDATA,,,,1,,,

MPTEMP,,,,,,,,

MPTEMP,1,0

! 输入FRP材料参数

MPDATA,EX,2,,Ef

MPDATA,PRXY,2,,EMUf

TB,bkin,2,1,2

TBTEMP,0

TBDATA,,430,0

! 输入型钢材料参数

MPDATA,EX,3,,Es

MPDATA,PRXY,3,,Esmuf             

TB,BISO,3,1,2,                 

TBTEMP,0                      

TBDATA,,235,2.06E4,,,,            

                             

! 建立组合柱模型

k,1,-34,-50,0

k,2,34,-50,0

k,3,34,-42.4,

k,4,2.3,-42.4,0

k,5,2.3,42.4,0

k,6,34,42.4,0

k,7,34,50,0

k,8,-34,50,0

k,9,-34,42.4,0

k,10,-2.3,42.4,0

k,11,-2.3,-42.4,0

k,12,-34,-42.4,0

a,1,2,3,4,5,6,7,8,9,10,11,12

pcirc,Rf,Rf+Fthick,0,360

pcirc,0,Rf,0,360

asba,3,1,sepo,,keep,

k,100,100,100,0

k,101,100,100,H

l,100,101

vdrag,all,,,,,,21

LDELE,21,,,1

numcmp,all

vglue,all

numcmp,all

! 赋予材性

VSEL,S,,,2                  !FRP管    2

VATT,2,0,2

VSEL,S,,,3                  !混凝土   1

VATT,1,0,1

VSEL,S,,,1                 !型钢      3

VATT,3,0,2

allsel,all

lesize,35,,,35$lesize,37,,,35$lesize,47,,,35$lesize,49,,,35$lesize,57,,,35

lesize,59,,,35$lesize,29,,,35$lesize,39,,,35$lesize,31,,,35$lesize,33,,,35

lesize,22,,,35$lesize,43,,,35$lesize,27,,,35$lesize,41,,,35$lesize,51,,,35

lesize,55,,,35$lesize,54,,,35$lesize,46,,,35$lesize,23,,,35$lesize,25,,,35

！划分环向

lesize,48,,,10$lesize,58,,,10$lesize,45,,,10$lesize,60,,,10$lesize,52,,,10

lesize,53,,,10$lesize,50,,,10$lesize,56,,,10$lesize,14,,,10$lesize,19,,,10

lesize,20,,,10$lesize,13,,,10$lesize,16,,,10$lesize,17,,,10$lesize,15,,,10

lesize,18,,,10

！划分工字钢截面

lesize,34,,,6$lesize,21,,,6$lesize,7,,,6$lesize,1,,,6

lesize,40,,,6$lesize,28,,,6$lesize,4,,,6$lesize,10,,,6

lesize,9,,,2.5$lesize,5,,,2.5$lesize,11,,,2.5$lesize,3,,,2.5

lesize,38,,,2.5$lesize,30,,,2.5$lesize,26,,,2.5$lesize,42,,,2.5

lesize,32,,,1$lesize,36,,,1$lesize,44,,,1$lesize,24,,,1

lesize,6,,,1$lesize,8,,,1$lesize,2,,,1$lesize,12,,,1

vsweep,1,1,15,1                 !型钢

vsweep,3,25,26,1                !混凝土      

vsweep,2,2,24,1                 !FRP管      

numcmp,all                     !压缩节点编号

nummrg,all

ALLSEL,ALL

nsel,s,loc,Z,H

Sf,all,pres,-0.0001

nsel,s,loc,Z,0

d,all,ux,0,,,,uy,uz

allsel,all

finish

/SOLU

eqslv,spar

pstres,on

Solve

finish

/SOLU

antype,1

bucopt,lanb,6,0,0

mxpand,6,0,0,1,0.001

Outres,all,all

solve

finish

/PREP7  

Upgeom,0.01,1,1,’C3’,’rst’,’’

finish

/solu

antype,0

nlgeom,on

outres,all,all

nsel,s,loc,Z,H

Sf,all,pres,70

allsel,all

nsubst,200

arclen,on

NEQIT,20

solve