第二章桥梁结构有限元法及可视化软件的开发1.docxVIP

第二章 桥梁结构分析有限元法及可视化软件开发 2.1 前言 有限元法是伴随着计算机技术进步而发展起来的一种数值分析方法， 广泛应用于分析各类复杂科学与工程问题。 由于它采用矩阵形式表达基本公式， 能灵活地处理和求解各种复杂问题， 便于应用于计算机编程运算等优点， 有限元法显示了其强大的生命力。 目前，各种有限元软件广泛应用于科学研究和工程应用领域。 但许多大型通用软件涉及领域广泛， 在具体问题的操作中难于满足工程设计人员各方面的要求； 一些有限元应用软件程序庞大， 不利于普通技术人员的掌握和使用。为此，结合工程实际编制易学、易用的可视化专用有限元软件能极大的提高 工作效率。本文开发的基于 Visual Basic 语言的大跨桥梁结构施工分析系统(BSCAS)，采用桥梁结构施工控制中普遍采用的正装法 (前进分析法 )和倒装法 (倒拆分析法 )实现了程序的结构化分析和施工模拟。 VB 友好的用户界面， 使设计过程生动直观， 具有有限元结构分析初步知识的工程技术人员在界面的引导下， 可以畅通无阻地进行应用设计。 2.2 桥梁结构有限元法概述 有限元法是将连续体分成有限个单元， 单元间相互有结点连接的理想结点系统。分析时，先进行单元分析，用结点位移表示单元内力，然后将单元再合成结构，进行整体分析，建立整体平衡关系，由此求出结点位移。 有限元求解程序的内部过程如图 2.2.1 所示： 2.2.1 桥梁结构有限元分析一般步骤 采用有限元法进行桥梁施工控制中的结构分析计算， 首先要建立数据文件。 数据文件准备一般分为四步： 桥梁结构的模型化。桥梁结构的模型化就是将实际结构理想化为有限个单元的集合。根据结构的受力特性与工作行为选择恰当的单元形式来模拟实际结构以及选择正确的约束模拟形式， 选择正确的约束模拟形式尤为重要。 就结构分析模型来看，与一般的已成桥梁分析不同的是施工控制中的结构分析模型一般是随着施工的不断推进而不断变化的， 这是由于桥梁在形成过程中的结构体系是在不断变化的。实际工作中， 通常可考虑建立一个统一的模型， 而对某个施工状态的结构模拟则可通过某些单元的是否激活来实现。 计算模型中单元的选择应以能准确描述施工过程中结构受力与变形状态为 准。对于桥梁结构而言， 最主要的是结构纵向的受力分析计算， 而将纵向分析模 型近似地处理成杆件系统， 往往是一种可取的简化模型。 一个模型可以是由一种单元组成，也可是由几种单元组成， 对一些特殊的施工工艺需要采用特殊的单元 结构离散化，输入或生成有限元网格 计算单元刚度矩阵形成总刚度矩阵 形成结点荷载向量 引入约束条件 解线性代数方程组 输出结点位移 计算输出单元的应力 图 2.2.1 有限元程序图 来描述。如连续刚构桥的墩梁固结处， 则假定墩顶零号块是一个具有无限大刚度的刚性杆件，应用程序的刚臂单元来处理。 桥梁结构的离散化。 桥梁结构的离散化就是在模型化处理后， 将结构离散为带有有限个自由度的结构， 这些单元在结点上相互联结， 构成一个与真实结构等价的计算模型， 即将真实的复杂空间结构模拟为平面的杆件系统。 单元大小 与结点位置确定应充分考虑结构受力情况与施工单元的划分。 结构离散化必须遵循以下原则 [5] ：1) 计算模型应尽量符合实际结构的构造特点和受力特点，以保证计算分析结果的真实性； 2) 保证体系的几何不变性，特别是在复杂的体系转换过程中更应注意， 同时要避免与实际结构受力不符的多余联结； 3) 在合理模拟的前提下，应尽量减少结点数目，减少未知量数目，以缩小计算规模，节省时间和计算机空间； 4) 本构关系的确定应尽量反映材料的性质。 杆件单元的划分，应根据结构的构造特点、实际问题的需要以及计算精度 的要求来决定，用来划分单元的结点应包括构件的转折点、 交接点、截面变化点。由于杆件结构有限元程序的计算结果是以结点位移及内力形式输出的， 故在进行结点划分时， 应根据设计需要验算的截面以及求影响线的作用点的要求， 来确定所需增加的中间结点。 按所用软件的输入要求形成数据文件 检查校正数据文件。 计算模型最终体现为数据文件， 数据文件正确方能保证计算模型的正确，乃至才能保证计算结果的正确性。 其次，运行分析软件。对于桥梁施工控制中的结构分析，由于计算模型随 着施工过程的改变，同时要求分析跟踪进行， 应采用具有施工控制跟踪、 仿真分析功能的软件。 最后，对分析结果进行分析和处理。 2.2.2 桥梁结构几何非线性有限元分析的单元刚度矩阵 在结构发生大位移时， 单元的外力和位移之间的线性关系不再适用。 可以将荷载分成几步逐渐施加， 从而把一非线性问题用一系列线性步骤来处理， 求出每一荷载增量引起的位移增量，即增量法。 在每一步中刚度矩阵可以分为两部分， 一部分是在每一步按当时的几何构造计算