钢筋砼柱抗震性能模拟例子.ppt

柱截面弯矩-曲率关系的分析与梁相似（柱存在轴力）压弯构件荷载-挠度单调曲线非线性全过程分析框图卸载Clough恢复力模型武藤恢复力模型武田恢复力模型三折线恢复力骨架曲线模型的验证三折线恢复力滞回模型的验证第一次程序大作业钢筋混凝土圆柱抗震性能试验和有限元模拟1.4试验结果2.OpenSEES有限元分析2.1.2非线性分析2.2OpenSEESTCL语言编程2.2.2建立模型3.定义结点质量mass$nodeTag(ndf$massValus)例：mass120;#定义结点1的x方向质量为2(此命令在此程序中未应用)当进行时程反应分析时需要指定质量。钢筋混凝土纤维模型示意图箍筋约束混凝土采用Mander模型Concrete01单调与滞回本构关系uniaxialMaterialConcrete01$matTag$fpc$epsc0$fpcu$epsu$lambda$ft$Ets(1)定义箍筋约束混凝土的应力-应变关系uniaxialMaterialConcrete01$matTag$fpc$epsc0$fpcu$epsu$lambda$ft$EtssetIDreinf2;#约束混凝土区编号setfc1f[expr–($fc0*(-1.245+2.245*sqrt(1+7.94*$y2/$fc0)-2*$y2/$fc0];#峰值应力seteps1f[expr2.0*$fc1f/$Ec];#峰值应变seteps2f[expr3.0*$eps1f];#极限应变setfc2f[expr($fc1f*($eps2f/$eps1f)*$R)/($R-1+pow($eps2f/$eps1f,$R)];#极限应力setlambda0.15;#极限应变的卸载斜率和初始斜率比setftu[expr-$fc1u/10.0];#抗拉强度setEts[expr-$Ec/2];#抗拉软化刚度uniaxialMaterialConcrete02$IDreinf2$fc1f$eps1f$fp2f$eps2f$lambda$ftu$Ets;注意：模型中的各参数设置需按照命令格式中的顺序，而参数符号可自行定义。实际上是主程序与子程序之间的实参虚参调用(2)定义未约束混凝土的应力-应变关系uniaxialMaterialConcrete01$matTag$fpc$epsc0$fpcu$epsu$lambda$ft$Ets;setIDconcu1;#未约束砼编号setfc1u$fc;#峰值应力seteps1u-0.002;#峰值应变setfc2u[expr0.2*$fc];#极限应力seteps2u-0.005;#极限应变uniaxialMaterialConcrete01$IDconcu$fc1u$eps1u$fc2u$eps2u$lambda$ft$Ets;Steel02单调与滞回本构关系(3)定义钢筋的应力-应变关系uniaxialMaterialSteel02$matTag$Fy$E$b$R0$cR1$cR2$a$b$c$d;setIDsteel3;#钢筋区编号setFy[expr353.4*$MPa];#钢筋屈服应力setEs[expr2.0e5*$MPa];#初始弹性模量setBs0.00085;#应变强化比setR018.5;setcR10.925;setcR20.15;#R0、cR1、cR2是钢筋由弹性阶段向塑性阶段的控制点6.纤维截面单元的划分纤维单元的划分sectionFiber$ColSecTag;#纤维截面patchcirc$matTag$numSubdivCirc$numSubdivRad$yCenter$zCenter$intRad$extRad<$startAng$endAng>;setnfcoreC25;#核心砼沿周长方向的纤维数setnfcoreR21;#核心砼沿半径方向的纤维数setnfcoverC25;#保护层砼沿周长方向的纤维数setnfcoverR2;#保护层砼沿半径方向的纤维数patchcirc$IDconcu$nfcoverC$nfcoverR0.0.$Rcore$Rcol0.360.;#对未约束砼区进行纤维元划分patchcirc$Idreinf$nfcoreC$nfcoreR0.0.0.$Rcore0.360.;#对约束砼区进行纤维元划分layercirc$IDsteel$numBarsCol$Abar0.0.$Rcore0.[expr360-360/$numBarsCol];#箍筋单元划分7.几何坐标转换9.结果输出11.轴向力的分析参量——力控制2.2.3非线性分析3结果输出及后期处理