翻译 新桥梁顶梁元.docVIP

新桥梁顶推梁元 摘要 一种新梁单元已被用于桥梁顶推过程中的静态和动态分析。该单元可以描述在梁上的任意点处的横向约束变形，从而模拟桥梁顶推过程中主梁支座位置的连续变化过程。此梁单元简化了桥梁顶推施工过程的有限元模型，并且大大降低了计算成本。给出的一些数值算例显示了建议单元的优越性。 关键词：梁，动态分析，桥梁，变形 引言 混凝土结构的顶推施工法被广泛应用于世界各地，因高强混凝土的质量能得到保证且其成本很低（Grant 1975）。使用顶推法施工时,需先在桥后台分节段预制梁体，然后将桥台上的整个梁段逐步抬高向前顶推为下一节段让出空间如图.1(Zellner and Svensson 1983)。对于桥梁设计和施工控制来说在顶推过程中结构的详细分析是必要的，该有限元方法（FEM）优选这样的任务，由于它功能的多样性。然而，使用商业软件是不方便的，因为在推出过程中没有特定的可用元素可以直接模拟连续变化的梁体支撑位置。因此，在不同的施工阶段需要用不同的有限元模型（FEM）使支撑位置的连续变化得以被考虑。此外，这种间接的方法在桥梁顶推过程中的动力分析是不合适的。在文献中，传递矩阵法（MTT）已用于分析桥梁的顶推过程(Sasmal et al.2004;Sasmaland Ramanjaneyulu 2006),但前面提到的缺点仍然存在。 这项工作的目的是制定一个适合静态和动态分析的顶推桥梁的梁单元。这种单元包括梁的纵轴线上任意位置的横向约束，因此在桥梁顶推过程中连续变化的支撑位置可以很容易的被合并。介绍一些元素验证已被证明的数值例子。 公式 梁单元横向支撑约束形函数 考虑到如图2.所示平面梁单元在x=d处的横向支撑约束，由于在顶推的过程中支点位置的连续变化，约束点不能在该位置处设置。对于共同的梁单元，通常假设以三项多项式坐标(Bathe and Wilson 1976) 的方式表示挠度v。然而，在支撑位置处一个三次多项式不能描述整个单元的不连续剪力。因此，假设挠度v是在支承点的两侧上的两个三次多项式的形式。即， 式中，为节点位移向量，其分量 和 分别为节点i、j的挠度和转角位移，并且形函数的矩阵为 式中 梁单元横向支承约束刚度矩阵与质量矩阵 根据有限元法的原理，可以得到刚度矩阵 其中EIZ表示梁弯曲刚度，并且 将式（3）~（5）代入式（6）和（7），可以得到刚度矩阵的确切表达试，这是在附录中给出的。以类似的方式，元素的质量矩阵可以一致表示为 其中ρ和A分别为元素的材料密度和横截面面积，并且 将式（3）代入（8）和（9）得出明确的质量矩阵形式，这也是在附录中给出的。 重力的等效节点力 在顶推的过程中，通常应考虑梁的重力，相应的等效节点力由下式给出 式中，g为重力加速度。 在实践中，梁单元将被连接以形成平面（2D）或空间（3D）的结构，其中可能会发生轴向以及扭转变形，。对于轴向延伸的平面结构元素，刚度矩阵 其中As=EA/l表示元素的轴向刚度且Kij(i,j=1,2,3,4)是单元的刚度矩阵由附录中式（13）~（16）给定。至于平面单元的质量矩阵， 式中，Am=ρAl表示单元的质量，Kij(i,j=1,2,3,4)是单元的刚度矩阵由式（8）给定，也在附录中明确给出。 对于一个空间结构单元，在x-z平面的也可能存在弯曲以及绕x轴的扭转。在顶推的过程中，通常提供一个方向上的支撑约束，也就是说y方向。然后，可以得到支承约束并在附录中给出的空间结构单元刚度矩阵和质量矩阵。 数值例子 为了验证本文提出的单元，FE程序使用MATLAB语言。已经研究过的几个数值例子。第一个研究对象为如图3.的两跨连续梁。建立2个有限元模型，采用传统的梁单元和目前梁单元的横向支承约束，分别为如图3（b和c）。采用商业软件程序ANSYS进行传统梁单元的结构分析。表1列出了有限元模型和弯矩区模型1-7位置处的连续梁挠度和转角。如果在位置4,的支承约束被移位到位置3，相应的挠度和转角表列于表2中。表3列出了连续梁的前三个频率，以及那些的解析方法的计算（Weaver等人，1990）。在位置3的连续梁横向约束结果也可以通过这个单元和所述解析方法得到并列于表4。 图.3.两跨连续梁的有限元模型 另一个例子用来证明本元素的有效性是一个在顶推过程中的拱形桥的动态分析。一座正在中国杭州某河上建造的拱桥。图.4给出了的FE模型，顶推过程的速度为2米/秒。顶推的过程的示意图是在显示于图.5。在启动过程，拱桥轴承约束沿着鼻导梁和两个边界梁（图4