/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