墩梁固接区设计空间局部有限元方法分析.docVIP

墩梁固接区设计空间局部有限元方法分析摘要：刚构桥中墩梁固接处的结构构造和应力分布情况都比较复杂，属于刚构桥中的关键部位。本文以柳林大桥主桥为例，运用有限元分析方法对墩梁固接处进行空间分析，以便指导工程的设计和施工。 关键词：墩梁固接处；有限元分析；模型 1工程概况 柳林大桥主桥中墩处，钢主梁、主翼和主墩刚性连接(如图 1所示)，该固接区构造极其复杂，地震荷载作用下的传力途径和局部应力分布情况十分不明朗。为此，本文基于偏于安全地认为，模态阻尼比0.02和50年超越概率2％的全桥反应谱作用是该墩梁固接区最不利的地震荷载作用工况，该工况下通过反应谱分析获得的墩梁固接区内力，可作为局部应力分析的作用荷载。 2有限元计算模型的建立 2.1材料 墩梁固接区除承台和桩基础外全采用钢材，且钢板所用材料统一为Q345qD。材料物理参数为：材料密度7.8 t/m3； 材料弹模2.1×1011 N/m；材料泊松比0.3，根据不同板厚对应不同的屈服强度和极限抗拉强度。 2.2单元介绍 该固接区局部应力分析采用大型计算程序ANSYS。因该区域设计均采用Q345qD钢板，故本分析采用 ANSYS单元库中的Shell63单元。Shell63单元为四节点实体单元，单元既具有弯曲能力又具有膜力，可以承受单元平面内荷载和垂直单元平面的法向荷载。该单元每个节点具有六个自由度――三个平动自由度和三个转动自由度。单元需要定义4个节点和4个厚度。 在单元面内，其节点厚度为输入的四个厚度，单元的厚度假定为均匀变化。如果单元厚度不变，只需输入参数TK(I)即可；如果厚度是变化的，则四个节点的厚度均需输入。 2.3有限元模型描述 固接区计算模型网格采用自由网格划分，根据结构设计所采用的不同钢板厚度，分批分区域形成计算网格，并最终融合形成整体有限元模型。模型由7578个关键点，13658条线和5506块面元构成，全部考虑了设计图纸中的主梁断面横坡、加劲肋、横隔板、过人孔，以及U肋、横隔板连接欠焊口等构造细节。最终计算模型包括496653个节点和499360个单元，得到的有限元分析模型见图2和图3。 因模型复杂，本报告建立有限元模型时，根据板材厚度确定不同的实常数，从而能将不同厚度的结构部位单元分开编组。相同厚度的单元归属于同一个单元组，方便计算结果处理、显示和分析。 由于结构比较复杂，计算网格的划分需要适用应力梯度的变化，本报告在固接区域内可能出现应力集中的部位进行单元局部细分。 网格检查表明，整个有限元模型网格划分质量满足ANSYS计算要求。 3有限元模型分析 3.1有限元模型验证 为了验证所建立的墩梁固接区有限元实体模型的正确性，通过比较关键节点位移和关键截面应力响应，分别验证主梁、主翼模拟和中墩模拟的正确性。 (1)主梁、主翼模型验证 本报告基于结构自重作用下的模型内力进行验证。在结构自重作用下，通过全桥有限元静力分析，获得的墩梁固接区边界内力如表1所示。这些内力包括主翼边界内力和主梁边界内力。 节间2边界主体结构的四个相交面弯矩扭矩正负号不同，以对称方式布载；主翼根部与主结构的两个相交面弯矩扭矩正负号不同，以对称方式布载。 将自重作用下的内力施加到墩梁固接区实体模型上，可得到边界面上的等效应力，计算结果见下表，下表也列出了墩梁固接区边界在全桥杆系模型中的等效应力值。 可见，节间2边界箱梁顶面应力二者吻合非常好。在节间2边界箱梁底面和主翼边界顶、底面，由于实体模型箱型截面的剪力滞效应和局部刚度变化，应力分布在某一范围变化，这一应力范围是以杆系模型相应位置的应力值为基准波动，其合理性是可以理解和接受的。 与全桥模型应力对比可知，本报告建立的实体模型之主梁、主翼部分较好地反映了全桥有限元模型中墩梁固接区主梁、主翼特征。 (2)中墩验证 为验证所建立的实体模型之中墩模拟合理性，本报告分别在实体模型和杆系模型取出的子结构对应点处加一相同的力，以比较两模型相应点处的位移。选取的集中力加载点为主翼边界面中心，力方向为顺桥向，大小为1×107N，具体加载方式见图4和图5。 由于在实体模型中不存在与杆系模型主翼单元形心位置相同的节点，这里采用主翼边界顶底板中心位移的平均值作为比较对象。选取的比较参数为该点顺桥向转动角度，主要考虑角度值受局部翘曲的影响较小。两种模型对应的顺桥向转动角度如表2所示。 可见，两种模型位移大体相等，由于整体建模时并未考虑过人孔对截面削弱的影响，故实体模型的整体刚度应小于整体模型的刚度，加之实体模型的空间和局部效应，因此局部实体模型的位移值大于整体模型的位移值是可以理解的，二