opensees模拟：有粘结实体桥墩分析算例

1、# Written by wyy# Units: N, m, secwipe;#wipe memery for new calculationmodel basic -ndm 2 -ndf 3;#2-d problem, 3-degree of freedom per node# Define constantsetPI expr2*asin(1.0)#等于 nsetg 9.81;# Data inputsethCol 0.36;# Height of ColumnsetbCol 0.50;# Width of ColumnsetcolAspect 4.86;# Column aspect r

2、atio: L/D 柱子的纵横比 setv expr0.1;# Column axial load rate -> P = v*fc*Ag轴向荷载率 setLCol expr$hCol*$colAspect;# Column length 柱长setfc -4.16E7;# Concrete compressive strengthsetfpc -5.747E+7setepsc0 -0.006328setfpcu -3.302E+7setepsu -0.0215#set ratio 0.3#set ft 5.474E+5setEc 2.43E10setcolArea expr$hCol*

3、$bCol;# Column cross-section area横截面积 setP expr2*$v*$fc*$colArea;# Column axial load 轴向荷载 #set P expr -576000-691484setm expr-$P/$g;# Beam no dal mass节点质量setbeamArea $colArea;# Beam cross-section areasetGJ 1.0e15;# Torsio nal stiffnes扭转刚度setnp 3;# Number of integration points# Define uniaxial materi

4、als# Core concrete （confined 约束混凝土 ）uniaxialMaterialConcrete011 expr$fpc $epsc0 $fpcu $epsu#uniaxialMaterial Concrete02 1 expr $fpc $epsc0 $fpcu $epsu $ratio$ft $Ets# Cover con crete （unconfined 无约束）uniaxialMaterialConcrete012 expr$fc -0.002 0 -0.005# STEEL# Cover steel Tube#set fy 3.67e8;# Yield st

5、ress#set E 2.0E11;# Young's modulus E=200000000kN/m# tag fy E0 b#uniaxialMaterial Steel01 3 $fy $E 0.01# STEELsetfy3.2E8setfu4.59E8setEs2.0E11setEsh6E9;# 0.03*Essetesh0.016;# 8*eyseteult0.08;# 40*eysetlsr 6.0setbeta 1.0setr 0.8setgama 0.8setCf 0.7setalpha 0.506setCd 0.7uniaxialMaterialReinforcin

6、gSteel 3 $fy $fu $Es $Esh $esh $eult-GABuck $lsr $beta $r $gama -CMFatigue $Cf $alpha $Cd# Post-tension threaded rod tag E0 eyt eye initial strain以负为正)setbeta1 1.0setr1 0.5setgama1 0.5setCf1 0.4setalpha1 0.3setCd1 0.3uniaxialMaterialReinforeingSteel 4 1.25E9 1.5E9 1.95E11 $Esh 0.060.3 -GABuek $lsr $

7、beta1 $r1 $gama1 -CMFatigue $Cf1 $alpha1 $Cd1#-# Post-tension threaded rod tag E0 eyt eye initial strain以负为正，通过输入初 始应变的方式施加预应力 )#uniaxialMaterial ElastiePP 4 1.95e11 expr 0.0086+0.0043 expr0.0086-0.0043 -0.0043# Define nodes# tag X Ynode1 0.0 0.0node2 0.0 expr0.2*$LColnode3 0.0 expr0.4*$LColnode4 0.

8、0 expr0.6*$LColnode5 0.0 expr0.8*$LColnode6 0.0 $LCol# Nodal mass# node MX MY RZmass6 $m $m 0.0# Boundary conditions# node DX DY RZfix1 1 1 1#fix 2 0 0 0#fix 3 0 0 0#fix 4 0 0 0#fix 5 0 0 0fix6 0 0 0# Define cross-section for nonlinear columns# set some parameterssetcolWidth 0.36setcolDepth 0.5setco

9、ver 0.024setAs expr$PI*pow(0.01,2)/4;# area of no. 7 bars n *D2/4PI= n setAPs expr7*$PI*pow(0.005,2)/4;# area of no. 7 bars# some variables derived from the parameterssety1 expr$colWidth/2.0setz1 expr$colDepth/2.0section Fiber 2 # Create the concrete core fiberspatchrect12020expr$cover-$y1expr$cover

10、-$z1expr$y1-$coverexpr$z1- $cover# Create the concrete cover fibers (top, bottom, left, right)patchrect 2 20 1 expr-$y1 expr$z1-$cover $y1 $z1patchrect 2 20 1 expr-$y1 expr-$z1 $y1 expr$cover-$z1patchrect 2 2 1 expr-$y1 expr$cover-$z1 expr$cover-$y1 expr$z1-$coverpatchrect221expr$y1-$coverexpr$cover

11、-$z1$y1expr$z1-$cover# Create the reinforcing fibers (left, right,top,bottom)layerstraight3 9 $As expr$cover-$y1 expr$cover-$z1 expr$cover-$y1 expr$z1-$coverlayerstraight3 9 $As expr$y1-$cover expr$cover-$z1 expr$y1-$cover expr$z1-$coverlayerstraight34$Asexpr$y1-$cover-(0.360-0.048)/5expr$z1-$covere

12、xpr$cover+(0.360-0.048)/5-$y1 expr$z1-$coverlayerstraight34$Asexpr$y1-$cover-(0.360-0.048)/5expr$cover-$z1expr$cover+(0.360-0.048)/5-$y1 expr$cover-$z1layerstraight4 3 $APs -0.13 -0.20 -0.13 0.20layerstraight4 3 $APs 0.13 -0.20 0.13 0.20# Define torsional stiffness and attach it to sectionuniaxialMa

13、terialElastic10 $GJ# attach torsion to flexure and create a new section IDtag# tag uniTag uniCode secTagsection Aggregator 1 10 T -section 2# Define column elements# tag vecxzgeomTransfLinear1# Plastic region# tag ndI ndJ nsecs secIDtransfTagelementnonlinearBeamColumn1 1 2 $np 11elementnonlinearBeam

14、Column2 2 3 $np 1 1elementnonlinearBeamColumn3 3 4 $np 1 1elementnonlinearBeamColumn4 4 5 $np 1 1elementnonlinearBeamColumn5 5 6 $np 1 1# Define beam element# tag $iNode $jNode $secTag# Record nodal displacementsrecorderNode -file S3-1.out -time -node 6 -dof 1 disprecorderplotS3-1.out nodeDispbpc 10

15、 10 300 600 -columns 2 1recorderElement -file date/steel-p.out -time -ele 1 section 1 fiberexpr$cover-$y1 expr$cover-$z1 stressStrainrecorderElement -file date/steel-l.out -time -ele 1 section 1 fiberexpr- $cover+$y1 expr$cover-$z1 stressStrainrecorderElement -file date/confined-concrete-p.out -time

16、 -ele 1section 1 fiber expr2*$cover-$y1 expr2*$cover-$z1 stressStrainrecorderElement -file date/confined-concrete-l.out -time -ele 1section 1 fiber expr-2*$cover+$y1 expr2*$cover-$z1 stressStrainrecorderElement -filedate/cover-concrete-p.out- time-ele 1 section1 fiber expr0.5*$cover-$y1 expr0.5*$cov

17、er-$z1 stressStrainrecorderElement -filedate/cover-concrete-l.out-time-ele 1 section1 fiber expr-0.5*$cover+$y1 expr0.5*$cover-$z1 stressStrain# # Constant gravity loadpatternPlain1 Linear # FX FY FZ MX MY MZload6 0.0 expr$P 0.0systemBandGeneralconstraintsTransformationnumbererRCMtestNormDispIncr1.0

18、e-12 1000 3algorithmKrylovNewtonintegratorLoadControl0.1analysisStaticinitializeanalyze10puts"Finish gravity load"loadConst-time 0.0#cyclic loadingsetH 1.0;# Reference lateral load# Set lateral load patternpatternPlain2 Linear load6 $H 0.0 0.0setDincr 0.001setDmaxPush 100integratorDisplace

19、mentControl6 1 $Dincrputs"开始 Pushover分析”setIDctrlNode 6setDxPush expr$DincrsetNDxPush expr-1*$Dincrdof init Jd min $DxPush ;analyze2;puts"2#nodemaxintegratorDisplacementControl$IDctrlNode 1 reached"integratorDisplacementControl$IDctrlNode 1 $DxPush ;analyze5;puts"3 reached"#

20、3Dx#-Dx#3DxintegratorDisplacementControl$IDctrlNode 1$NDxPush ;analyze6;puts"-3 reached"integratorDisplacementControl$IDctrlNode reached"integratorDisplacementControl$IDctrlNode5 reached"integratorDisplacementControl$IDctrlNode reached"integratorDisplacementControl$IDctrlNod

21、eDx#3Dx#-Dxanalyze14;puts"-7 reached"integratorDisplacementControl$IDctrlNode reached"integratorDisplacementControl$IDctrlNode10 reached"integratorDisplacementControl$IDctrlNode reached"integratorDisplacementControl$IDctrlNode15 reached"integratorDisplacementControl$IDc

22、trlNode reached"integratorDisplacementControl$IDctrlNode20 reached"integratorDisplacementControl$IDctrlNode reached"1 $DxPush ;analyze8;puts"51 $NDxPush ;analyze10;puts"-1 $DxPush ;analyze12;puts"71 $NDxPush ;#-Dx#3Dx#-1 $DxPush ;analyze17;puts"101 $NDxPush ;analyz

23、e20;puts"-1 $DxPush ;analyze25;puts"151 $NDxPush ;analyze30;puts"-1 $DxPush ;analyze35;puts"201 $NDxPush ;analyze40;puts"-1 $DxPush ;analyze45;puts"251 $NDxPush ;analyze50;puts"-integratorDisplacementControl$IDctrlNode reached"1 $DxPush ;analyze55;puts"30

24、integratorDisplacementControl$IDctrlNode30 reached"1 $NDxPush ;analyze60;puts"-integratorDisplacementControl$IDctrlNode reached"1 $DxPush ;analyze65;puts"35integratorDisplacementControl$IDctrlNode35 reached"1 $NDxPush ;analyze70;puts"-integratorDisplacementControl$IDctr

25、lNode reached"1 $DxPush ;analyze75;puts"40integratorDisplacementControl$IDctrlNode40 reached"1 $NDxPush ;analyze80;puts"-integratorDisplacementControl$IDctrlNode reached"1 $DxPush ;analyze85;puts"45integratorDisplacementControl$IDctrlNode45 reached"1 $NDxPush ;anal

26、yze90;puts"-integratorDisplacementControl$IDctrlNode1 $DxPush ;#3Dx#-Dx#3Dx#-Dx#3Dx#-Dx#3Dx#-Dx#3Dx#-Dx#3Dx#-Dx#3Dx#-Dx#3Dx#-Dx#3Dxanalyze95;puts"50reached"integratorDisplacementControl$IDctrlNode 1 $NDxPush ;analyze100;puts"-50 reached"integratorDisplacementControl$IDctrlNo

27、de 1 $DxPush ;analyze105;puts"55 reached"integratorDisplacementControl$IDctrlNode 1$NDxPush ;analyze110;puts"-55 reached"1 $DxPush ;#-Dx#3Dx#-11 $DxPush;integratorDisplacementControl$IDctrlNode$NDxPush ;analyze120;puts"-60 reached"integratorDisplacementControl$IDctrlNod

28、e 1$DxPush ;analyze125;puts"65 reached"integratorDisplacementControl$IDctrlNode 1$NDxPush ;analyze130;puts"-65 reached"integratorDisplacementControl$IDctrlNodeanalyze135;puts"70 reached"integratorDisplacementControl$IDctrlNode 1$NDxPush ;analyze140;puts"-70 reached"integratorDisplacementControl$IDctrlNode 1 $DxPush ;analyze145;puts"75 reached"integratorDisplacemen