基于全尺寸自主开发材料的桥梁钢筋混凝土墩柱施工的ABAQUS动力分析.docxVIP

基于全尺寸自主开发材料的桥梁钢筋混凝土

墩柱施工的ABAQUS动力分析

万鹏

江西省交通运输厅

摘要：本文提出两个纤维梁单元混凝土模型。基于全尺寸自主开发材料，模拟了桥梁钢筋

混凝土墩柱施工的ABAQUS动力分析，开发了两个适用于纤维梁单元的混凝土模型用户子程

序，通过对比全尺寸振动台试验各项结果，验证了本文开发的子程序的可靠性，结果表明：

用户子程序能很好地反映受荷载作用。

关键词混凝土结构，地震荷载，材料模型，数值分析，ABAQUS

1.引言

钢筋混凝土的非线性性能越来越受到行业研究者的重视，试验研究的成果也越来越广

泛。除钢筋混凝土结构的动力试验越来越被研究者重视外，钢筋混凝土结构的动力非线性数

值分析也在学术界渐渐受到关注。各种钢筋混凝土结构动力分析的理论研究与材料模型研究

也越来越受到学者们的青睐，ABAQUS软件具有成熟高效的非线性算法和功能强大的前后处

理能力，结构非线性分析的研究者更愿意在这类通用平台上二次开发自己的材料模型或单元

来进行结构的非线性分析。本文也是在ABAQUS通用平台上开发出用于纤维单元的混凝土

材料模型以及钢筋模型，并用这些开发的模型对钢筋混凝土桥墩柱进行非线性动力分析。

2.混凝土模型理论模型

本文混凝土受压受拉骨架曲线如图1所示，其中受压曲线分成4段，受拉曲线分成2段。

受压曲线第一段从坐标原点到压应力峰值的40%为结束点，为直线段，反应混凝土的弹性特

性。第二段是应变从最大弹性应变经压应力峰值对应的峰值应变再到极限压应变为二次曲线

段，反应混凝土的非线性特性。第二段的应力应变关系，采用EN1992欧盟混凝土规范中

建议的应力应变关系公式。

图1混凝土模型受压受拉骨架曲线

本文在ABAQUS中开发的F1-CON和F2-CON模型，在没有实测弹性模量和峰值应变

时，输入参数为3或者4个，前3个参数分别对应图1：其中F1-CON模型不计混凝土拉应

力的影响。两种模型的卸载和重加载路径分别与OpenSees中CONCRETE01和CONCRETE02

模型一致，此处不再赘述。文中计算ABAQUS自带混凝土两种模型所需要的参数时，混凝

土拉压应力应变关系均使用图1所示的骨架曲线，这样无论用ABAQUS自带混凝土模型还

是本文开发模型，进行数据对比时，就不会因混凝土骨架曲线不同而导致没有可比性。本文

的钢筋模型采用作者开发的钢筋混凝土结构的钢筋滞回模型[7]，在ABAQUS上开发材料模

型的接口程序同文[8]。

3.桥墩柱试验简介

该试验是PEER于加州大学举办的一个钢筋混凝土桥墩柱动力响应试验测试竞赛[1][2][3]

(如图2照片所示)，桥墩柱的截面配筋如图3。

图2试验照片图3柱截面尺寸、配筋和单元分割

图3中，纵筋为直径36mm，箍筋2根直径16mm，间距153mm，箍筋外的混凝土保

护层厚度51mm，具体材料力学性能参数参见文[1][2][3]。试验时柱顶配有一个质心位于柱顶

的混凝土质量块，重2322kN，即质量236.7吨。

4.模态分析

首先进行模态分析，确定系统的固有频率。试验时由于与振动台激励方向垂直的方向上，

试件位移受到约束，因此，可以认为系统在作平面振动。在ABAQUS平台上系统的模态分

析，整个柱子长7.315米，平分成15个节点，每2个节点间采用线性平面铁摩辛柯梁单元

B21。考虑到柱子箍筋对混凝土的约束作用，同时ABAQUS对圆形截面平面梁单元的截面积

分最大积分点数只能是9，也即，从轴线到梁的受拉或受压边缘的辛普森积分点只有5点，

这样的积分点数很难准确描述如此巨大截面的混凝土非线性特性。因此，本文分析时，每2

个节点之间采用3个B21单元来描述柱子截面(如图3所示)。其中最外围的圆环截面表示

柱子混凝土保护层截面，环的厚度51mm，ABAQUS梁单元截面特性选择SECTION=PIPE。

图3中间的环，取以钢筋为中心位置的圆环，其厚度为67.5mm，ABAQUS截面选择同保

护层一样。这两个圆环的积分点都选择11点，这样相当于将环分成20等分。由于外围的

环是非约束混凝土，同时，其混凝土的弹性模量和峰值应变均为已知，因此，ABAQUS的输

入命令如下：

*MATERIAL,NAME=F1-CONe

*Depvar

10,

*UserMaterial,constants=6

40.9,.2,0.004,0.01,22900.,0.0026

上述命令中第5行的参数分别为：混凝土峰值应力，残存应力比值，最大压应