索桥结构的非线性有限元法.pdf

索桥结构的非线性有限元法 z 屈本宁 空 —1趸 力学系) 摘 要 结合新 型倒 张拱 索桥 结构特 点，考 虑索桥非线性 大变形性质，提 出三 节点= 次 由线元，直梁单元和杆 单元共存 的混合 单元有 限元模 式，从 结构整体上模拟索 桥 ．由Reissner广 义变分 原理 导 出由巍元，梁单元的刚度 矩阵，井得 到用节点坐 标表示 的大住牡空间杆 的单元剐度矩 阵，实现 7索、梁 杆单元的位移协调，得 到整体 规模下 的有限元运 动方程．用 Newton-Raphson方法求解静 力 问题 ．实例 的计算结果与实测位对 比吻合较好 ，表 明方法是正确和可行的． 关键词 倒张境索桥：非缉些·太变形；位移协调 — ’ - 0引言 暇元t，京强 倒张拱索桥是在索道桥基础上发展起来的新型柔性桥梁．它与原索道桥在结构上的区 别是在桥面下方增加了预拉倒张索．如图 l所示．从而大大增加了桥的稳定性，克服了原 索道桥变形过大，稳定性差的致命弱点 索桥主要由钢索、横梁 拉杆及桥板组成，是复杂的大变形组合结构，受力及稳定性 分析均较复杂．本文应用非线性有限元理论，结合索桥的结构特点，提出索 梁、杆单元 相结合能从整体上模拟桥体结构的混合有限元模式_由Reissner变分原理导出索．粱单元剐 度矩阵，并得到用坐标表示 的大位移空间杆单元刚度据 阵．在 3个假设的前提下实 现了三种单元 的位移协调， 从而建立了整体规模下的有 限元运动方程 ．应用求解静 横粱 力问题的Newton-Raphson 田1倒张拱襄桥 方法编制了计算程序．以跨度 105m，载重 30t的索桥为算例，计算结果与实测结果对 比具 有较好的一致性．表明混合的有限元模式是正确而实用的． 本文的理论和方法对索网类结构计算具有重要意义， 1各单元刚度矩阵，载荷列阵疑质量矩阵 索桥受力后的变形特征是大位移小应变，这类问题的求解可归结为求 Reissner广义 变分原理提出的泛函的变分，此泛函为 n=J。 {盯((，+u+(：『(，：)一；]一}dI，一』G．d一』Gd (1) 式中， 为广义应力张量 ， 为位移， “： ／。 ，c为体膨胀系数， 为温度 变化，E 为余应变 比能，V。为物体体积，S。为表面积． 1．1柔索单元剐度矩阵，载荷列阵殛质量矩阵 [1=孚J [{】Ⅳ【Ⅳ }一1—2筇 N[c]d (2) 索单元几何刚度矩阵为 ， [，]= 』 ([ⅣⅣ[]f}{}Ⅳ【lⅣ【】)d (3) 载荷列阵为 {R}={r】+J [Ⅳ l，}山+，J ![_v]f^}d (4) 质量矩阵为 ’ M【I=J ．，m[N] 【Id (5) 以 各 式 中 -Ⅳ『 0 0N；0 0Ⅳ 0 0] [Ⅳ]=Ⅳ【]：d；loⅣC。ooⅣ；ooⅣol l0 0_v 0 0Ⅳ；0 0NJ 而Ⅳ—一+等，Ⅳ=一等ⅣC=；+等 S为索单元拉 氏坐标， 为材料弹模；A为索单元截面积，，，为索单元变形前、后长 度；m为单位索长质量；g为熏力加速度；{ }‘为索单元节点坐标列阵，{r} 为集 中载荷列阵 (，}为变形后单位索长上的分布力 (含分布惯性力)列阵；f^} 为重力方 向与坐标的关系列阵． 1．2 梁单元刚度矩阵，载荷列阵及质量