基于纤维模型及柔度法钢管混凝土压弯剪构件力学性能模拟及灵敏度分析.docVIP

基于纤维模型及柔度法钢管混凝土压弯剪构件力学性能模拟及灵敏度分析摘 要：本文基于有限元柔度法和纤维模型法的基本思想，研究了往复荷载作用下的CFT框架柱滞回性能的建模方法，借助非线性梁柱单元，编制了钢管混凝土压弯剪构件的有限元分析程序。分析了影响模型计算精度的因素。建立了不同的钢管混凝土框架柱计算模型，讨论了材料本构关系、截面纤维划分、积分点数量等因素对模拟结果的影响。以试验为依据，运用OPENSEES软件，通过计算结果和试验结果的对比，分析了影响数值模拟准确性的因素。结果表明，受约束混凝土本构关系的下降段主要影响到滞回曲线的退化，钢材本构关系的强化段主要影响到水平荷载的大小，二者是影响模型准确性的重要因素，而网格划分则主要体现在计算效率上。本研究为钢管混凝土框架柱滞回性能的数值模拟分析提供了理论依据。 1. 引 言 目前研究结构的非弹性静力和动力反应的主要方法有两种。其一是结构静力和动力试验的方法，一般采用小比例尺的结构模型试验，以模拟原型结构的性能。小比例尺的模型试验往往带有较大的尺寸效应[1]，而大比例尺以及足尺的模型试验往往受到试验条件以及研究费用的限制，难以得到推广。 其二是计算机仿真试验方法或数值模型试验方法。该方法不受时空条件的限制，节省研究费用，易于追踪结构受力全过程中的性能变化，适用于对结构进行大量的参数分析等等。与传统的有限元法相比，纤维模型法能够在保证计算精度的同时耗费较小的机时，因此纤维模型法更适合进行复杂结构的非线性分析。目前国内外基于纤维模型的数值模拟越来越多的应用于试验验证，因此有必要深入分析影响模拟精度的因素。本文基于纤维模型的基本原理借助Opensees[3]开发平台研究了钢管混凝土框架柱的数值模拟方法，探讨了材料本构关系、截面纤维划分、积分点数量等因素对模型计算精度的影响。 2.纤维模型和有限单元柔度法的计算原理 2.1 纤维截面 纤维截面以其高精度模拟材料的非线性而得到学者的广泛认可。基本理论是通过一系列纤维截面将梁沿x轴方向划分为数段，每一段的特性以中间截面来代表。将中间截面离散为纤维网格，每个网格可定义不同的材料本构关系，截面特性由材料的本构关系通过截面积分确定。然后由截面的力与变形关系得到杆端力与位移关系。作为一种有效的截面离散化方法，纤维截面法广泛适用于各种宏观有限元程序。 2.2 有限单元柔度法 基于柔度法建立非线性梁柱单元时，需要对单元节点力（本文将单个杆件设定为一个单元则节点力为杆端力）引入力的插值函数，由静力平衡关系得到截面力。对于假定的内力分布，在不考虑单元荷载分布任意变化的前提下，无论杆件处于何种状态，作为单元控制方程之一的平衡条件都能严格满足。对比刚度法，柔度法应用的难点在于单元抗力难以直接从截面抗力推导得到，只有通过迭代逼近的方式获得。虽然这需要加大计算量，但其优势体现在通过引入力的插值函数可在任意计算时刻由单元节点力得到精确的单元内力（截面力），迭代收敛时单元的内力和抗力一致，单元的变形和单元节点位移协调，单元变形可以呈现复杂的分布形式，模拟精度比刚度法有很大的提高。 3. 钢管混凝土压弯剪构件的建模方法 3.1 材料模型 材料的本构模型是影响数值计算精度的最主要因素之一。钢管混凝土中的混凝土由于被钢管所约束，强度可以得到大幅度的提高，进而提升构件的综合性能。为了研究材料本构关系对计算结果的影响，本文基于不同的材料本构模型，进行了CFT柱在循环荷载作用下的数值模拟。 3.1.1 混凝土材料模型 约束混凝土的本构关系中最为典型的为mander[6]模型和韩林海[7]模型，二者的曲线特征点的计算存在着一定差异。本构关系的数学表达式分别见文献[6]和文献[7]。 3.1.2 钢材材料模型 钢材硬化段对计算结果的影响较为明显，此为本研究考察的对象之一。本文选择Giuffre-Mengegotto-Pinto模型，分考虑强化和不考虑强化两种情况。当考虑钢材强化时，强化段后的硬化模量取为0.01，为钢材的弹性模量，下文简称（模型S1）；当不考虑强化时，硬化模量取为0，下文简称（模型S2）。 以上两种钢材模型在反复荷载的作用下考虑材料刚度退化，加卸载过程均考虑包辛格效应，采用随动强化理论模型。 采用Mander JB[6]等人和韩林海[7]等提出的两种混凝土本构模型组合两种钢材本构模型，应用纤维模型法建模计算得到的滞回曲线名称如表1所示。 表1滞回曲线简称 3.2 纤维截面的划分和积分点数设定 纤维截面网格的划分方式、疏密程度和积分点个数会影响到程序求解的精度和效率，为了分析截面的划分疏密程度对模拟的影响，采用适合于圆形截面特性的划分方式，将截面按环向