OpenSees的高强钢筋混凝土柱抗震性能数值模拟.docVIP

基于OpenSees的高强钢筋混凝土柱抗震性能数值模拟 摘要：基于OpenSees平台，采用梁柱纤维单元建立了钢筋混凝土单柱的纤维模型，对循环荷载下的配置HRB500纵筋和箍筋的混凝土单柱非线性滞回反应进行了数值模拟。通过对模拟结果和试验结果的对比分析，验证了OpenSees平台对于配置高强钢筋混凝土柱的抗震性能模拟的可行性和精确性。 关键词：OpenSees;高强钢筋混凝土柱；纤维模型；数值模拟；抗震性能 0 引言 我国目前的用钢水平远低于国际先进水平，根据混凝土结构的发展趋势，推广应用高强钢筋将促进我国结构用钢水平的提高，在建筑结构中推广应用 HRB500钢筋，具有巨大的经济效益和社会效益。2010年《混凝土结构设计规范》正式将HRB500级钢筋纳入主力钢筋并用于混凝土梁、柱和剪力墙的设计。 目前，高强钢筋用于的试验研究，pa级纵筋抗压强度、高强约束箍筋抗拉强度和高强钢筋对柱的位移延性的影响等问题上存在争议，有必要进行进一步的试验研究和理论分析，为高强钢筋的推广和应用提供科研依据。该单元允许刚度沿杆长变化，首先通过结点位移得到相应的单元杆端位移，然后根据位移插值型函数求得截面的变形，再根据截面的恢复力关系得到相应截面抗力与截面切线刚度矩阵，最后按照Gauss—Legendre积分方法沿杆长积分计算出整个单元的抗力与切线刚度矩阵。 式中，为单元位移插值型函数，为截面切线刚度矩阵。 2 材料本构关系 2.1 混凝土材料本构模型 混凝土模型的应力-应变滞回关系曲线，如图1所示。混凝土模型没有考虑混凝土受拉力学性能，即受拉区混凝土的应力和刚度为零。 图 时， 时， 时， 其中： 各个参数代表的含义如下：为混凝土峰值应力对应的应变，为考虑约束引起的强度和峰值应变增大系数，是应变下降段的斜率，为混凝土圆柱体抗压强度，是核芯区箍筋的屈服强度，是横向约束箍筋的体积配箍率，是混凝土核芯区的宽度，是箍筋间距。 对于混凝土受箍筋约束的情况，Scott等建议混凝土极限压应变可偏保守地按下式确定： 2.2 钢筋材料本构模型 钢筋纤维采用Giuffre一Menegotto—Pinto[4]钢筋本构关系，其应力一应变关系如图2所示。 图 钢筋骨架曲线为双折线型，该模型由于采用了应变显函数表达形式，因而计算方便，并且保持了与钢筋反复加载试验结果较好的一致性，可以反映筋的Bauschinger[5]效应。模型在数学上可以表示为 其中：与分别为钢筋屈服点处的应力与应变；与分别为钢筋反向点处的应力与应变；为钢筋的硬化刚度与原点切线模量的比值；R为反映Bauschinger效应的参数；R0，与均为材料常数；为上一循环塑性应变的绝对值。 3 单柱模拟算例 模拟了同济大学2009年所做的4根配置500Mpa高强钢筋的混凝土柱的低周反复荷载试验表试件参数 试件 编号 (Mpa) 试件配筋 轴压力 (kN) (Mpa) (Mpa) (%) 纵向钢筋 箍筋 Z5F-1 42.2 4DF16 DF8@100 115.9 0.11 527.3 568.4 1.51 0.19 3.16 Z5F-2 42.2 4DF16 DF 8@60 231.9 0.22 527.3 568.4 2.52 0.31 3.16 Z5F-3 55.2 4DF16 DF 8@60 188.7 0.14 527.3 568.4 2.52 0.24 3.16 Z5F-4 55.2 4DF16 DF8@100 377.4 0.28 527.3 568.4 1.51 0.14 3.16 注：DF表示HRB500级钢筋，B表示HRB335级钢筋；为混凝土立方体抗压强度实测值； 为箍筋屈服强度实测值；为加密区配箍特征值；为加密区体积配箍率；为纵筋实测屈服强度；为试验轴压比；为剪跨比。试件 (a)Z5F-1 (b)Z5F-2 (c)Z5F-3 (d)Z5F-4 图5骨架曲线对比图 (a)Z5F-1 (b)Z5F-2 (c)Z5F-3 (d)Z5F-4 由上图可以看出试验滞回曲线与模拟滞回曲线基本吻合，相较之下模拟滞回曲线更为饱满，由于试验过程中竖向力施加没有维持恒定，导致模拟结果有所偏差。骨架曲线的特征值基本相同，误差不超过6%。 5 结束语 用OpenSees程序采取合理的单元类型和本构关系，合理的划分单