一种分析复杂桥梁结构的新单元.pdfVIP

一种分析复杂桥梁结构的新单元 项贻强 徐兴 (浙江大学交通工程研究所) 摘要：本文提出了一种分析复杂桥梁结构的新单元，它突破了传统的空问有限元数值分析 法对结构的不同材料越须采用不同单元的限制，用较少的单元就能够方便模拟和描 述各种复杂的结构。可广泛适用于具有任意支承条件的复杂桥梁的空间整体分析、 静动力分析及非线性问题的分析。 关键词：屡舍单元理论等参元空阃分析 1问题的提出 对于具有复杂截面形状(如变高度的箱形面、T形截面和Ⅱ形截面粱)、支承条件和线形 变化的斜弯桥梁、斜拉桥，若要进行精确的分析，理论上都应归结为空间问题，应用解决空 间问题的方法进行计算。传统的考虑荷载横向分布的平面设计理论和弹性理论解析解法在这 方面几乎很难奏效，传统的三维板壳单元和块体单元等有限元数值分析法和通用程序，理论 上可通过对结构进行细致的剖分，来解决工程各种复杂空间计算问题，但由于其建模时输入 数据量大、单元划分多(有时多达十万，甚上百万个单元)，不同单元有时拼装较为困难，而 且有些夹在结构物中的材料特性难以描述和模拟，致使某些问题的分析受到了限制。 本文在传统的空间有限元法基础上，假定结构的单元为由若干的层合或段组成的单元， 将单元的刚度矩阵、单元质量矩阵、单元结点力矩阵在单元区域内的积分用单元内若干层 或段的小区域的积分之和来代替，而每一层和段可以分别描述不同的几何形状和材料特性， 这样就可以方便地将结构物的整体形状和特性描述出来。但是，仅将单元分层和分段描述 进行积分还不能避免解题单元细分、计算量庞大的困难。因此，在此基础上，进一步引进 “虚拟层”和“虚拟段”的概念(必要时为描述几何形状、变形特性还应引入虚拟并非真 实的结点)。所谓“虚拟层”或“虚拟段”，是指该“层”或该“段”实际上可以没有材料， 或假想该“层”或该“段”的某些材料特性参数(如弹性模量E、G等)为零，如箱梁结 构的中空的箱室部分等，其材料的弹性模量E、G等可以为零。即该“层”或该“段”为 空。这样通过层间夹“空层”、“空段”或层外补“空层”、“空段”来构造单元，一方面， 一个单元可有多种分层的不同材料，一层中还可以有多种材料，甚至可以没有材料，使原 来需要多个单元才能进行分析和多层不同材料归为一个单元内部的层，从而减少了单元数， 提高了计算效率；另一方面，使一些具有复杂几何形状的结构可以由单元内部的分层或分 段来处理，从而避免了单元局部细节的结点连接带来的一系列问题。这就是“虚拟层合单 元”的独到之处，根据这一思路，它不仅可以粗放地划分单元，并且可以方便地考虑各种 材料的差异及分布情况。 2层合单元的基本理论 2．1坐标变换式和位移插值模式 根据一般的有限元理论，对8～20结点的空间等参单元，可由边长为2的立方体母单元 的坐标变换而得到，其坐标变换式和位移模式可统一写成如下形式： H。 ixi ∽掣鲁m。搓、●●●●●●●●L，●●●●●●』 肿严肌f惶 (2) 式中：n为单元的结点数；五、Y；、zi和“；、vi、Wi分别为结点的整体笛卡尔坐标和 位移，Ⅳ；为形函数。当n=8时，其形函数为： (f=1,2，…8) (3) N{=去(1+毒{专)(1+17i叼)(1+f{f) 其中：f目f为边长2的立方体母单元的坐标系，其原点在母单元的形心处。量，叩f，f。 的取值范围为母单元的边界，即磊，仇，色=+l或-1。 当n=20时，其边中结点的形函数为： ⅣI=丢(1髻)(1+W：D(IKf)(f=9，lo，1112) 4’l516) (4) (Ⅲ3’1 Ⅳl=三(1可)(1幢跏Kf) 其角结点的形函数由如下的线性组合得出： 12 N。