钢筋混凝土非线性分析第一次大作业

Harbin Institute of Technology

钢筋混凝土结构非线性分

析课程作业

院 系： 姓 名： 学 号: 时 间：

?哈尔滨工业大学

基于Opensees钢混柱结构滞回曲线比较分析

一、试验资料

本实验结构如下图所示，柱净高850mm，截面尺寸为250mm*250mm，采用C40混凝土材料，受力纵筋采用HRB335，箍筋采用HPB235。柱子采用6根直径为12mm的HRB335级钢筋，柱身箍筋采用直径6mm的HPB235级钢筋，箍筋间距为50mm。混凝土保护层厚度为25mm。竖向施加760KN荷载，轴压比为0.268。

混凝土强度：

钢筋力学性能：

试验测得的滞回曲线：

二、Opensees建模过程

1、主程序

wipe;

source Units.tcl; source GeometricParameters.tcl; source Material.tcl; source FiberSection.tcl; source Elements.tcl; source RecorderRC.tcl; source PointGravityLoad.tcl;

source Ex4.Portal2D.analyze.Static.Cycle.tcl; 2、定义量纲

set NT 1.0; set mm 1.0; set sec 1.0;

set kN [expr 1000.0*$NT]; set MPa [expr 1.0*$NT/pow($mm,2)]; set m [expr 1000.0*$mm]; set mm2 [expr $mm*$mm];

set mm4 [expr $mm*$mm*$mm*$mm]; set cm [expr 10.0*$mm]; set PI [expr 2*asin(1.0)]; set Ubig 1.e10; set Usmall [expr 1/$Ubig]; puts \3、定义节点

model BasicBuilder -ndm 2 -ndf 3; node 1 0 0; # node#, X, Y node 2 0 $Lb;

fix 1 1 1 1; fix 2 0 0 0; # node DX DY RZ puts \completely ||\

4、定义几何参数

set B [expr 250.*$mm];

set H [expr 250.*$mm]; set Ab [expr $B*$H]; set Izb [expr 1./12.*$B*pow($H,3)];

set Lb [expr 875.*$mm]; set Cover [expr 25.*$mm]; set Bcore [expr ($B-2*$Cover)]; set Hcore [expr ($H-2*$Cover)]; set Acor [expr $Bcore*$Hcore]; set NmBar_12 4; set Dlbar_12 [expr 12*$mm]; set Albar_12 [expr $PI*pow($Dlbar_12,2)/4.0]; set CAlbar_12 [expr $NmBar_12*$Albar_12]; set CAlbar [expr $CAlbar_12]; set Fy_12 [expr 380.2*$MPa]; set Es_12 [expr 200000.*$MPa];

set Dhbar [expr 6.*$mm]; set Ash [expr $PI*pow($Dhbar,2)/4.0]; set Fyh [expr 588.3*$MPa];

set Esh [expr 210000*$MPa]; set xyt [expr $Fyh/$Esh];

set LDhb [expr 2*($Bcore-$Dhbar)+2*($Hcore-$Dhbar)]; set SV [expr 50*$mm]; set SVc [expr $SV-$Dhbar]; puts \

5、定义材料本构关系

set IDsteel_12 1;

set IDcoverC 3; set IDcoreC 4;

set fc [expr 45.4*$MPa]; set x0 0.002; set ft [expr $fc/10.]; set Et [expr $ft/0.002]; set fc0 -$fc; set xc0 -$x0;

set fcu0 [expr 0.2*$fc0]; set xcu0 -0.006;

set lambda 0.10;

uniaxialMaterial Concrete01 $IDcoverC $fc0 $xc0 $fcu0 $xcu0;

set ps [expr $CAlbar/$Ab]; set pcc [expr $CAlbar/$Acor];

set pst [expr $Ash*$LDhb*1.5/($Acor*$SV)]; set lmda [expr $pst*$Fyh/$fc]; set fcc [expr $fc0*(1+0.5*$lmda)]; set xcc [expr $xc0*(1+2.5*$lmda)]; set fcu [expr 0.35*$fcc]; set xcu [expr 20*$xcc];

uniaxialMaterial Concrete01 $IDcoreC $fcc $xcc $fcu $xcu;

set R0 18.5; set cR1 0.925; set cR2 0.15;

set haRatio 0.0037;

uniaxialMaterial Steel02 $IDsteel_12 $Fy_12 $Es_12 $haRatio $R0 $cR1 $cR2;

puts \

6、定义纤维截面

set RCSecTag 1; set b1 20; set b2 1; set h1 20; set h2 1;

set hs1 [expr $Hcore/2-$Dhbar-$Dlbar_12/2]; set bs1 [expr $Bcore/2-$Dhbar-$Dlbar_12/2]; set coverY [expr $H/2]; set coverZ [expr $B/2]; set coreY [expr $Hcore/2]; set coreZ [expr $Bcore/2]; section Fiber $RCSecTag {;

patch quad $IDcoreC $b1 $h1 -$coreY $coreZ -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ;

patch quad $IDcoverC $b1 $h2 -$coverY $coverZ -$coverY -$coverZ -$coreY -$coverZ -$coreY $coverZ;

patch quad $IDcoverC $b1 $h2 $coreY $coverZ $coreY -$coverZ $coverY -$coverZ $coverY $coverZ;

patch quad $IDcoverC $b2 $h1 -$coreY $coverZ -$coreY $coreZ $coreY $coreZ $coreY $coverZ;

patch quad $IDcoverC $b2 $h1 -$coreY -$coreZ -$coreY -$coverZ $coreY -$coverZ $coreY -$coreZ;

layer straight $IDsteel_12 3 $Albar_12 $hs1 $bs1 $hs1 -$bs1;

本文来源：https://www.bwwdw.com/article/fpjd.html