用实体退化板壳单元对立交桥进行的空间分析.docVIP

用实体退化板壳单元对立交桥进行的空间分析 　　摘要:采用实体退化板壳单元对杭州复兴立交桥的单箱多室箱梁的结构进行空间分析,得出应力和挠度与试验结果相符,证明了实体退化板壳单元的正确性和有效性。 　　关键词:实体退化板壳单元;单箱多室;有限元;空间分析 　　 　　1 概述 　　 　　箱型截面梁具有自重轻、抗扭刚度大、便于现代化施工等优点,因而广泛应用于桥梁建设上,立交桥中应用更为普遍。虽然近20年来,单箱单室静定箱梁扭转、畸变、翘曲、剪力滞后等方面的研究已经取得一些比较好的解决方法,但这几个影响因素一般被分开来研究。然而,大型、多室、多箱的超静定梁的使用越来越多,对多箱多室结构仍运用单箱单室箱梁的分析方法会产生较大的误差。50年代出现的有限单元法在理论上成功的解决了箱梁在多种作用下的受力问题。鉴于箱梁的外形特征,一般用传统的板壳单元来分析,如4节点12(或24)自由度矩形协调单元;用二次分片插值法的4节点12自由度四边形协调单元;用条件极值法建立的协调单元;用分项插值法建立的协调单元等。对于三广义位移的板壳理论,最常用的是8节点24自由度的等参元。但传统的板壳单元是不协调单元,且有转角自由度,因而单元间连接处理较麻烦,对于变厚度的处理也十分粗糙,而最近提出的实体退化板壳单元能很好的解决以上问题。本文便采用这一种单元对单箱多室的杭州复兴立交桥进行空间分析,取得了令人满意的结果。 　　 　　2 实体退化单元 　　 　　由于用位移法对于经典板壳理论构造协调的板壳单元存在一定困难,板壳单元的构造一直是有限元领域中研究得最多的问题,基于板壳理论得到的板壳单元种类繁多,且层出不穷,又由于对简洁明了的等参单元的偏爱,所以希望把板壳结构当作一种特殊的块体结构,采用三维实体单元的等参单元来进行分析。徐兴和凌道盛通过引入相应的板壳假定,得到一系列适用于薄板壳、中厚板壳、膜、梁结构分析的实体退化单元。 　　如图1所示的16节点板壳单元,每个节点有u、v、w三个自由度,共48个自由度 　　单元坐标和位移插值形函数和三维实体单元相同。记16节点三维等参元的节点位移矢量为 　　 (1) 　　其形函数、应力向量、应变向量、应变矩阵均与8～20节点三维等参元的相同。 　　2.1.基本假定 　　经典板壳理论是一般的三维弹性理论根据板壳结构特殊的几何形状引入一定的简化假定后得到的,可以认为板壳理论是一种特定条件下简化的三维弹性理论。从三维弹性理论直接导出的是三维实体等参元,由此可以看出,板壳单元其实是一种特定条件下简化的三维实体等参元,只要在三维实体等参元中引入必要的简化的假定即可。取z坐标方向为板法线方向。 　　2.1.1 中厚板单元基本假定 　　 　　引入以上24个约束,厚壳单元由16节点40自由度单元变成16节点24自由度厚板单元,与通常的8节点24自由度厚板元相同。 　　2.1.2 中厚板壳单元的材料本构关系 　　根据?滓z=0的假定,引入?棕上=?棕下的约束条件和罚系数的办法,采用与三维等参元相同的应力应变向量,各向同性材料的本构关系为 　　 　　式中Ni为三维块体等参元的形函数,?籽为质量密度,|J|为Jacobi行列式值。 　　2.3 相对位移的引入和刚度矩阵的修改 　　如果按(5)式的本构关系,则必须引入 　　 的约束条件,为此要将绝对位移?棕变换为相对位移?驻?棕,这种变换可在单元中进行。 　　考虑如下形式的单元平衡方程 　　 　　3 杭州市复兴立交桥的空间分析和试验研究 　　 　　杭州市复兴立交工程桥梁部分A主线跨铁路(即25#～28#桥墩处)的结构为33+65+46m的三跨变高度预应力砼连续梁(图2),并且A主线在26#桥墩处与H、J两条匝道相连(图3)。桥梁宽度为变宽度,变化范围为19m～35.22m。该桥设计荷载为汽-20、挂-100,人群4.0kN/m2。 　　 　　3.1 计算模型与分析方法 　　根据复兴立交工程A25～A28跨连续梁的施工设计图,采用三维退化板壳单元及浙江大学交通工程研究所的程序USAP,建立了该桥空间三维有限元模型(图4)。通过计算求出各试验控制截面的位置(图3),并且求得各控制截面的影响线。分析中没有考虑普通钢筋以及预应力束的影响;支座采用刚性支座,A25支座的纵向位移为零。 　　3.2 计算与分析比较 　　本次检测利用三维退化板壳单元及USAP,计算各控制截面的应力分布,通过比较得到各控制截面的最大应力点,试验中对这些最大应力点附近的点作为试验控制点,各截面最大控制点位置及控制应力见表1;另I-I截面理论挠度值与实际值比较见表2。 　　从表中可以看出,试验结果与理论值符合较好。理论值较试验值总体偏大的主要原因在于截面实际刚度大于计算刚度