<4D F736F F D FD0B1C0ADC7C5CAA9B9A4B9FDB3CCCBF7C1A6BEABC8B7C4A3C4E2BCBCCAF52E646F63>

Size: px

Start display at page:

Download "<4D F736F F D FD0B1C0ADC7C5CAA9B9A4B9FDB3CCCBF7C1A6BEABC8B7C4A3C4E2BCBCCAF52E646F63>"

侍翟
5 years ago
Views:

1 浅谈斜拉桥施工过程中索力精确模拟技术 [ 桑登峰戴宇文胡若邻李治学 ] [ 中交四航工程研究院有限公司广州 ] [ 摘要 ] 斜拉桥在国内外高速公路及高速铁路桥梁上的应用越来越广泛, 而在斜拉桥施工过程中斜拉索是结构的关键部位, 也是建造过程中控制的要点, 如何精确模拟计算施工过程中斜拉桥索力是斜拉桥的建设其中一项关键技术由于非线性等因素的影响, 采用常规计算方法在计算过程中施加索力与目标索力往往有很大偏差, 难以满足施工精度要求本文将引入迭代法和影响矩阵法这两种方法来模拟计算施工过程的索力, 并结合工程实例验证这两种方法在斜拉桥施工过程索力模拟计算中达到较高的精确度, 将施加索力与目标索力的偏差控制在 1% 以内, 提高了斜拉桥施工过程中索力模拟计算的效率与精度 [ 关键词 ] 施工过程 ; 影响矩阵 ; 迭代 ; 索力 Precise Simulation of Cable Force for Construction State of Cable-Stayed Bridge [Sang Dengfeng Dai Yuwen Hu Ruolin Li Zhixue] [CCCC fourth harbor engineering institute Co.,ltd.,Guangzhou,510230] [ Abstract ] Cable_Stayed Bridge is used more and more on highway and high-speed railway,enter abstract information here. whereas the cables are the vital structures of cable-stayed bridge,so precise simulation of cable force for construction state of cable-stayed bridge is a key technology,because of non-linearity, the simple methods can t simulate the cable force precisely.two methods based on iteration and influence matrix are introduced, The feasibility is verified by cable force simulation of one cable-stayed bridge on construction.the error is less than 1%.so the precise simulation of cable force is improved obviously. [ Keyword ] construction state; influence Matrix; iteration; cable force 1 前言斜拉桥在国内外高速公路及高速铁路桥梁上的应用越来越广泛, 斜拉索是斜拉桥的生命线, 斜拉桥结构内力对斜拉索索力变化比较敏感, 如何精确模拟计算施工过程中斜拉桥索力是斜拉桥的建设的关键技术目前国内外众多学者提出了不少算法, 包括初始应变法和温度荷载法模拟索力 [1-2], 其公式简单实用, 但与目标索力具有一定的偏差, 而且这种偏差并会逐步累积, 可能会使计算最终结果的产生较大的偏差, 导致理论计算难以指导施工

2 [3-4] 为了减少这种施加索力与目标索力的偏差, 本文将引入迭代法和影响矩阵法这两种方法来模拟计算施工过程的索力, 并结合工程实例验证这两种方法在斜拉桥施工过程索力模拟计算中达到较高的精确度, 有效减少了施加索力与目标索力的偏差, 大大提高斜拉桥施工过程中索力模拟计算的效率与精度 2 斜拉索索力计算非线性影响因素斜拉桥是索塔梁组成的组合结构, 属于高次超静定的柔性体系, 复杂桥梁结构在荷载作用下受有轴力弯矩剪力扭矩等作用, 其截面应力处于复杂应力状态, 其受力状态很难用数学公式来表达, 结构体系呈现明显的几何非线性 [5], 斜拉桥结构内力和变形对斜拉索索力变化比较敏感, 斜拉索索力的变化对整个结构体系影响很大, 索力模拟受斜拉索垂度效应和斜拉索水平分力引起的主梁 P- 效应及柔性结构大变形等非线性因素的影响 1) 斜拉索垂度效应, 用杆单元近似模拟索类构件, 由于斜拉索垂度的影响导致单元刚度变化, 特别是在施工阶段斜拉索初张索力较小的情况下, 非线性影响越发明显, 索结构自重对结构平衡的影响不可忽略, 计算中可根据斜拉索索力变化使拉索单元在工程结构中表现的弹性模量用 Ernst 公式不断修正, 不断更新拉索单元刚度矩阵 Ernst 公式表达式即 : 其中, 为 Ernst 修正后的有效弹性模量, 为拉索材料弹性模量, 为拉索单位体积质量,L 为拉索的水平投影长度, 为拉索初应力 Ernst 公式修正弹性模量, 在小位移高应力水平下, 精度较高, 但是索在大位移状态或者应力水平不高的情况下会产生较大的误差 [6] 2) 斜拉索水平分力引起的主梁 P- 效应, 斜拉索索力使主梁及桥塔等构件处于弯矩和轴力组合作用下, 在这种组合作用下, 即使这些构件的材料满足胡克定律处于线性范围内, 但由于其单元几何刚度矩阵发生变化, 平衡方程是建立在变形后的位置上的, 因此整个结构体系也会呈现非线性特性 3) 结构大变形, 由于斜拉索等柔性结构的应用, 斜拉桥结构具有大变形小应变特性, 即当荷载作用在斜拉桥结构的某节点上, 该节点将发生位移, 荷载也随之移动, 荷载的作用方向也随之改变, 荷载对其他节点的弯矩也将改变, 如果位移过大, 结构刚度矩阵是几何非线性的因此其平衡方程将不再是线性关系, 在几何非线性分析时, 应考虑大变形的影响, 平衡方程是建立在结构变形后的位置上的 (1) 3 常用索力模拟法及其局限性结论索力模拟可以对索定义初应变来模拟索力, 对于索力的初应变计算方法为 : (2)

3 其中为斜拉索等效弹性模量, 可用 Ernst 公式不断修正,A 为斜拉索单元截面积, F 为索力, 为初应变还可以对斜拉索施加温度荷载法来模拟索力, 对于温度荷载施加的计算方法为 : T 为需施加的温度荷载, 为斜拉索线膨胀系数以上两种方法虽然通过 Ernst 公式修正其垂度弹性模量, 而且公式简单明了, 操作方便, 但是其目标索力及其平衡公式是建立在结构变形前的位置, 而索力的大小与其锚固两端的位移相关, 由此带来的非线性影响因素较大, 特别是在施工阶段中, 模拟索力初张拉时, 由于此时索力较小, 索力应力水平较低,Ernst 公式修正误差会进一步加大导致模拟索力与目标索力有更大偏差, 并且这种偏差会在后面施工阶段中逐步积累, 导致最终计算结果有较大偏差, 直接影响工程质量 (3) 4 索力精确模拟计算方法由于施工过程仿真分析时常用的模拟索力方法无法精确模拟索力, 很难有效模拟实际施工过程, 笔者参考前人关于索力模拟优化的方法, 引入一般迭代法和基于影响矩阵的迭代法, 用于施工过程初张索力的模拟, 对比两种方法优劣 4.1 循环迭代法由于非线性因素的影响, 采用施加单元初应变或者温度荷载来模拟张拉力一次模拟误差较大, 可以基于上述两种方法通过迭代的方法提高模拟精度, 即基于温度荷载的迭代法和基于初应变的迭代法, 基于初应变的迭代法的表达式为 : 式中, 为目标索力 ; 为第 n 次迭代完成时的索力, 为迭代的精度要求当不满足式 (5) 收敛准则时, 通过公式 (4) 继续迭代计算, 直到满足公式 (5) 的收敛条件基于温度荷载的迭代法的表达式为 : (4) (5) 同样根据公式 (5) 定义的收敛准则进行迭代, 直到满足公式 (5) 的收敛条件 (6) 4.2 基于影响矩阵的迭代法影响矩阵法优化索力也是建立在结构满足线性叠加的基础上, 在非线性结构体系中也是通过线性计算原理并通过线性迭代不断逼近目标值, 达到模拟精度需求因此, 非线性结构也采用线性计算原理, 通过不断的迭代来获取最终数值, 根据线性叠加原理 [7] :

4 (7) 则 : (8) 式中为影响矩阵 ; 为施调向量 ; 为受调向量, 为的逆矩阵在斜拉桥施工模拟过程中, 当模拟到第 n 束斜拉索张拉阶段时, 该斜拉索索力初次模拟之后, 实际索力与目标索力的偏差, 也即受调向量为, 为调整向量, 就是该工况下的索力影响矩阵, 如果调整向量为应变, 影响矩阵中的元素为第 j 束斜拉索单位应变增量对第 i 束斜拉索索力的改变量 ; 如果调整向量中的元素因此 : 为第 j 束斜拉索单位温度增量对第 i 束斜拉索力的改变量为温度, 影响矩阵由于影响矩阵以满足线性叠加原理为前提条件, 影响矩阵是索力调整之前的影响矩阵, 由于几何非线性的影响, 当索力变化时, 结构体系也发生变化, 影响矩阵也会随之变化, 采用影响矩阵法一次模拟也难以达到高精度, 因此, 也需要通过迭代法不断逼近目标索力同样可以通过设定公式 (5) 的收敛准则, 直到满足收敛条件 ; 基于影响矩阵迭代公式表示为 : (9) (10) (11) 5 工程实例某大跨度钢桁梁斜拉桥跨径布置为 ( )m, 主塔为花瓶形钢筋混凝土索塔, 在锚索区下端设置上横梁, 在主梁下端设置下横梁, 塔高米, 斜拉索共 9 4 2=72 根, 主梁为钢桁梁, 主桁采用整体节点, 上弦杆下弦杆腹杆都采用箱型截面采用 ANSYS 软件, 建立全桥模型其中, 主梁和主塔均为 Beam188 单元, 斜拉索为只受拉的 Link10 单元, 桥面板采用 Shell63 单元模拟拉索与主梁节点之间通过刚性杆连接 [8], 建立的有限元模型如下图所示拉索编号以主塔为界, 边跨侧从索塔处向梁端依次为 S0~S8, 中跨侧从索塔处向跨中依次为 M0~M8

5 图 1 斜拉桥有限元分析模型该桥跨度大, 非线性影响大, 拉索之间相互影响大施工方法从主塔位置向两边对称悬拼, 施工过程中随着钢梁的拼装和拉索的张拉, 结构体系在不断的变化, 而且初张索力较小, 索力应力水平较小, 施工过程的初张索力采用不同方法模拟的结果见表 1 和图 2. 其中, 方法一是直接施加初应变 ; 方法二是循环迭代法, 迭代 10 次 ; 方法三是基于影响矩阵的迭代法, 迭代 10 次表 1 索力模拟结果对比拉索编号目标索力 (kn) 方法一 (kn) 误差 % 方法二 (kn) 误差 % 方法三 (kn) 误差 % S S S S S S S S S M M M M M M M M M

6 目标索力 (kn) 5000 目标索力 (kn) N k µ Ö Á Ë 方法一 (kn) 方法二 (kn) 方法三 (kn) N k µ Ö Á Ë 方法一 (kn) 方法二 (kn) 方法三 (kn) S0 S1 S2 S3 S4 S5 S6 S7 S M0 M1 M2 M3 M4 M5 M6 M7 M8 斜拉索编号斜拉索编号 A)S0~S8 拉索模拟结果对比 B)M0~M8 拉索模拟结果对比图 2 索力模拟结果对比分析表 1 和图 2 可见, 直接施加初应变的方法模拟拉索索力, 误差很大, 随着施工过程的积累, 误差不断增大, 超过了 30%, 可见施工过程中此种方法很难模拟到精确索力, 不能满足工程精度需要 ; 而循环迭代法和基于影响矩阵的迭代法模拟索力收敛较快, 误差都小于 1%, 都能够满足工程精度需要 ; 对比两种方法, 一般的迭代法收敛速度相对慢一些, 误差在 0.9% 以内, 基于影响矩阵的迭代法收敛精度更高误差基本控制在 0.4% 以内, 基于影响矩阵的迭代法收敛更快 6 结论循环迭代法和基于影响矩阵迭代法是两种较精确模拟索力计算的方法, 从应用情况来看其精度基本能控制在 1% 以内, 满足工程需要两种方法各有特点, 前者公式简单方法明了, 但是迭代收敛速度不如后者 ; 后者收敛速度快, 但其影响矩阵的计算相对繁琐因此, 在斜拉桥斜拉索模拟计算过程中要根据实际工程情况采用相应的计算方法进行计算现今, 由于计算机的计算速度日益加快, 采用循环迭代法模拟更有利于工程人员的理解和掌握, 适用于一般跨径结构的斜拉桥而对于超大跨径复杂的斜拉桥, 一般迭代法收敛速度过慢, 造成过多的耗用计算资源从而加大了计算的难度, 在这种情况下建议采用影响矩阵迭代法进行模拟计算, 可加快收敛速度, 优化计算资源, 有利提高计算精度与速度满足工程需要 [ 参考文献 ] [1] [1] 周强, 杨文兵, 杨新华. 斜拉桥索力调整在 ANSYS 中的实现 [J], 华中科技大学学报 ( 城市科学版 ), 2005,22(5):81-83 [2] [2] 叶梅新, 韩衍群, 张敏,ANSYS 二次开发技术在确定斜拉桥初始恒载索力中的应用 [J], 铁道科学与工程学报,2005,27(5):56-59 [3] [3] 李义强, 张彦兵, 杨丽. ANSYS 中准确施加斜拉桥索力方法的研究 [J]. 国防交通工程与技术. 2006(1): [4] [4] 张杨永, 吴万忠, 周云岗. 斜拉桥索力精确模拟的矩阵分析法 [J]. 重庆交通大学学报 ( 自然科学版 ). 2009, 28(6):

7 [5] [5] 卫星, 强士中, 利用 ANSYS 实现斜拉桥非线性分析, 四川建筑科学研究,2003,12:47 [6] [6] 肖汝诚, 项海帆. 斜拉桥索力优化的影响矩阵法 [J]. 同济大学学报. 1998, 26(3): [7] [7] 项海帆, 姚玲森, 高等桥梁结构理论, 人民交通出版社.2002:247 [8] [8] 包龙生, 杨炳成, 于玲. 吊拉组合索桥的模型建立与初始索力确定 [J]. 沈阳建筑大学学报 ( 自然科学版 ).2004,20(4):

Fig1 Theforceappliedtothetrainwhenrunning :w = w j +w q (3) :w = w = w 0 +w j (4) w i 121 基本阻力 w r = 600 R ( N/kN) (8) :R : [2] w s [3] w s =0

Fig1 Theforceappliedtothetrainwhenrunning :w = w j +w q (3) :w = w = w 0 +w j (4) w i 121 基本阻力 w r = 600 R ( N/kN) (8) :R : [2] w s [3] w s =0 31 4 2012 8 JournalofLanzhouJiaotongUniversity Vol31No4 Aug2012 :1001-4373(2012)04-0097-07 * 张友兵张波 ( 100073) : 分析了列车运行过程中的受力情况给出了制动过程中减速度的计算方法并采用正向反向两种迭代方式计算列车制动曲线两种方式计算出的制动曲线一致证明了计算制动曲线的方法是正确的