楊恒，李元松，祁超，周小龍

武漢工程大學資源與土木工程學院，湖北武漢430074

大跨徑自錨式懸索橋主纜的線形計算及誤差分析

楊恒，李元松*，祁超，周小龍

武漢工程大學資源與土木工程學院，湖北武漢430074

針對自錨式懸索橋的主纜線形在施工過程中的變化特性，自編分段懸鏈線法和拋物線法程序?qū)崿F(xiàn)了大跨徑自錨式懸索橋主纜線形、空纜線形、各節(jié)段無應(yīng)力長度、吊點坐標和索鞍預偏量等幾何特征參數(shù)的精確計算；結(jié)合工程實例，對自編程序與Midas civil軟件計算結(jié)果進行了比較分析.結(jié)果表明：與有限單元法相比，分段懸鏈線法計算得到的索鞍預偏量相對誤差為11.3%，拋物線法計算得到的索鞍預偏量相對誤差為23%.在施工中，分段懸鏈線法和拋物線法計算結(jié)果需要反復修正迭代以綜合確定主纜線形.

自錨式懸索橋；主纜線形；分段懸鏈線法；拋物線法

1 引言

自錨式懸索橋因其造型新穎、跨越能力強、對地質(zhì)條件要求低等特點，現(xiàn)已成為大、中跨徑橋梁中具有較強競爭力的一種橋型，尤其是在軟土地區(qū)和城市景觀橋梁中，自錨式懸索橋越來越受到青睞［1-3］.隨著國民經(jīng)濟和交通事業(yè)的快速發(fā)展，人們審美需求的不斷提高，橋梁的跨越能力和靚麗外觀逐漸占據(jù)重要地位，相應(yīng)地伴隨跨度的不斷增大，橋梁的受力與變形特性更加復雜，其中幾何非線性特征也愈顯突出［4］.精確計算自錨式懸索橋的成橋線形與施工過程中的變化特性已成為該類橋梁設(shè)計與施工控制的關(guān)鍵技術(shù)問題，而主纜線形又是成橋線形的決定因素之一.

本文針對柔性結(jié)構(gòu)的受力特點，基于靜力平衡理論，詳細討論了主纜線形的計算公式與參數(shù)取值方法，自編Fortran程序?qū)崿F(xiàn)了自錨式懸索橋主纜線形的精確計算，結(jié)合工程實例，給出自編程序與Midas civil軟件計算結(jié)果，并進行了比較分析.

2 計算方法

2.1 基本平衡方程

現(xiàn)代大跨度懸索橋的主纜一般由鋼絲集束而成，抗彎剛度很小，可按完全柔性處理.如圖1所示，對于固定于A、B兩端，受均布荷載w作用的索段，基本假定如下［5］：

1）柔索僅能承受拉力而不能承受彎矩；

2）成橋時吊桿為豎直方向；

3）索的橫截面在荷載作用下不發(fā)生變化.

圖1 均布荷載作用下索段力學模型（a）索段AB和（b）微索段1-2Fig.1Mechanical model of cable segments under uniform load（a）segment of cable AB and（b）micro-segment of cable 1-2

對于長度為dx的曲線段1和曲線段2，由力的平衡條件，得單邊懸索的基本平衡方程［6-7］：

式（1）中：H為主纜拉力的水平分量，ω為均布荷載集度.

按荷載沿水平方向均勻分布或沿主纜長度均勻分布假設(shè)，式（1）的解有拋物線和懸鏈線2種.

2.2 拋物線法

拋物線法［8-9］假定主纜自重qc和橋面荷載qb沿跨長均勻分布，ω=qc+qb，如圖2所示.

2.2.1 成橋主纜線形

基本平衡方程（1）的解為：

根據(jù)邊界條件可求得主纜線形方程.

對兩端等高的主纜、中跨：

圖2 拋物線法主纜受力模式Fig.2Force model of main cable by parabola method

式（3）和式（4）中L1、L2、f1、h的定義見圖2.

2.2.2 成橋狀態(tài)索長

分別將式（3）、式（4）對弧長積分可得各節(jié)段變形后的索長［10］：

中跨的彈性伸長量：

邊跨的彈性伸長量：

2.2.3 空纜線形和索鞍預偏量

設(shè)索鞍預偏量為d，空纜狀態(tài)的水平分力為H#.主纜中跨、邊跨跨度分別為：

中跨、邊跨矢高分別為：

將上述幾何參數(shù)分別代入式（5）～式（7）中，可以得空纜狀態(tài)各節(jié)段的無應(yīng)力長度.

由于預偏量d和水平力H#均為未知量，不能直接求出定解，根據(jù)各跨主纜在成橋恒載狀態(tài)和空纜狀態(tài)下的無應(yīng)力索長均應(yīng)相等，經(jīng)迭代可求得空纜狀態(tài)d和H#的解.

同理，根據(jù)吊桿間各索段的無應(yīng)力長度不變原則，可依次迭代計算出各吊點的縱橋向坐標，代入空纜線形方程，從而求出空纜線形.

2.3 多段懸鏈線法

多段懸鏈線法的主纜受力如圖3所示，假設(shè)主纜自重qc沿主纜均勻分布，橋面恒載qb轉(zhuǎn)化為吊桿集中力Pi.此時平衡微分方程為［11-12］：

圖3 多段懸鏈線法的主纜受力圖示Fig.3Force diagram of main cable by segmented catenary method

2.3.1 成橋恒載狀態(tài)對第i段進行受力分析，微分方程（10）的解為［13-15］：

第i段主纜的高度：

第i段主纜的索長：

忽略主纜彈性模量隨應(yīng)力大小的變化，設(shè)主纜的彈性模量為Ec，截面積為Ac，則第i段索的彈性伸長量為：

2.3.2 空纜狀態(tài)

2.3.2.1 索鞍預偏量

設(shè)索鞍預偏量為d#，空纜狀態(tài)主纜水平力為H#，主纜的中跨、邊跨跨度分別為

根據(jù)各跨主纜在成橋恒載和空纜狀態(tài)下的無應(yīng)力索長相等的條件進行迭代求解.此時，求解d和H#.

2.3.2.2 主纜線形

以主纜狀態(tài)下的中跨任意第i段主纜為例，設(shè)左右兩吊點的坐標分別為(xi'yi)和(xi+1'yi+1).

有應(yīng)力索長為：

彈性伸長量為：

根據(jù)成橋和空纜兩種狀態(tài)下，吊桿之間各段的無應(yīng)力索長相等的原則進行迭代運算.在已知

xi的基礎(chǔ)上，求得各吊點的縱橋向坐標xi+1，將

xi+1代入懸鏈線方程即可得各吊點的空纜狀態(tài)坐標，從而求得空纜線形.

2.3.3 計算流程及程序?qū)崿F(xiàn)歸納前述各工況的計算公式，形成如圖4所示計算流程，并自編Fortran程序.

圖4 多段懸鏈線法計算流程圖Fig.4Flowchart of calculation for main cable by segmented catenary method

3 工程實例

3.1 基本數(shù)據(jù)

某橋為雙索面鋼箱加勁梁自錨式懸索橋，立面布置如圖5所示.橋梁基本參數(shù)：中跨吊桿間距L1=600 m（40×15 m），f1=60 m，邊跨吊桿間距L2=195 m（11×15 m+30 m），高度h=65.2 m；單根主纜截面Ac=0.257 8 m2，Ec＝2.0×105MPa，單位長度重qc=20.436kN/m；單根吊桿截面Ad=5.439×10-3m2，Ed=2.0×105MPa，單位長度重qd=0.424 kN/m；加勁梁截面Ab=3.066 7 m2，Eb=2.06×105MPa，單位長度重qb=165.2 kN/m（加二期恒載）.

圖5 自錨式懸索橋結(jié)構(gòu)體系簡圖（單位：m）Fig.5Structural profile of self-anchored suspension bridge（unit:m）

3.2 計算結(jié)果

3.2.1 主纜坐標3種方法計算的成橋恒載狀態(tài)下主纜坐標、空纜狀態(tài)下主纜坐標，分別列于表1、表2中.

表1 成橋恒載狀態(tài)下的主纜坐標（單位：米）Tab.1Coordinates of main cable under dead load state（unit：m）

表2 空纜狀態(tài)下的主纜坐標（單位：米）Tab.2Coordinates of main cable under empty load state（unit：m）

從表1和表2可以看出，成橋恒載作用下，3種方法計算的吊點坐標較為一致，最大相對誤差為0.2%.空纜狀態(tài)下，最大相對誤差為0.14%.分段懸鏈線與有限單元法計算結(jié)果較為吻合.

3.2.2 索長和索鞍預偏量主纜索長和索鞍預偏量對橋梁的幾何構(gòu)形和初始受力狀態(tài)具有較大影響，若存在較大誤差，將會對施工階段的索塔受力帶來安全隱患，須加以重視并解決.索長列于表3、表4中，索鞍預偏量列于表4中.

表3 成橋恒載狀態(tài)的索長Tab.3Cable under dead load state

表4 空纜狀態(tài)的索鞍預偏量和索長Tab.4Pre-bias of cable saddle and cable length under empty load state

從表3和表4可以看出，成橋狀態(tài)和空纜狀態(tài)下，有應(yīng)力索長和無應(yīng)力索長，多段懸鏈線法與Midas civil計算結(jié)果誤差較??；索鞍預偏量存在誤差，多段懸鏈線法更接近有限單元法，拋物線法誤差較大，誤差為0.15 m，相對誤差為23%.

4 結(jié)語

1）基于分段懸鏈線法自編程序計算結(jié)果與Midas Civil軟件計算結(jié)果基本吻合，用于大跨度自錨式懸索橋成橋、空纜線形分析，能夠達到設(shè)計和施工要求的精度.

2）對于成橋恒載狀態(tài)，分段懸鏈線法和拋物線法均可滿足工程設(shè)計和施工的精度要求，但空纜的吊點坐標和索鞍預偏量結(jié)果表明，拋物線法誤差較大.

3）分段懸鏈線與拋物線法均沒有考慮散索點的約束對線形的影響，且計算精度決定于對各吊桿力的計算與分配，因此仍存在一定的近似性，正確的處理方法應(yīng)與有限元法計算結(jié)果對比分析，經(jīng)反復修正迭代后，綜合防治確定.

［1］WOLLMANN G P，OCHSENDORF J A，ILLINGTON D P.Self-anchored suspension bridges［J］.Journal of bridge engineering，2001，6（2）：156-158.

［2］張哲.混凝土自錨式懸索橋［M］.北京：人民交通出版社，2005.

［3］OCHSENDORF J A，BILLINGTON D P.Self-anchored suspension bridges［J］.Journal of bridge engineering，1999，4（3）：151-156.

［4］王楨.超大跨徑自錨式懸索橋幾何非線性行為理論與實驗性能研究［D］.重慶：重慶交通大學，2013.

［5］邱文亮.自錨式懸索橋非線性分析與實驗研究［D］.大連：大連理工大學，2003.

［6］劉釗，劉厚軍.懸索橋主纜變形及重力剛度新算法［J］.工程力學，2009，26（6）：127-132.

LIU Z，LIU H J.New arithmetic for cable deflection and gravity stiffness of suspension bridges［J］.Engineering mechanics，2009，26（6）：127-132.

［7］肖海波.自錨式懸索橋的力學特性分析［D］.杭州：浙江大學，2004.

［8］張志國，張慶芳，鄒振祝.懸索橋成橋狀態(tài)計算方法［J］.長安大學學報（自然科學版），2006，26（1）：59-62，66.

ZHANG Z G，ZHANG Q F，ZOU Z Z.Analytical methods of suspension bride geometry under dead loads［J］.Journal of Chang'an university（natural science edition），2006，26（1）：59-62，66.

［9］王國興.自錨式懸索橋主纜線形控制與吊索張拉分析［D］.西安：長安大學，2008.

［10］張炎.自錨式懸索橋主纜線形及吊索張拉方案研究［D］.重慶：重慶交通大學，2014.

［11］楊勇.非對稱懸索橋主纜線形程序開發(fā)與參數(shù)分析［D］.重慶：重慶交通大學，2010.

［12］張佳文，趙彬.懸索橋主纜成橋線形的解析迭代計算方法［J］.湘潭大學自然科學學報，2011，33（1）：47-52.

ZHANGJW，ZHAOB.Calculationmethodof iterative analysis on main cable shape of suspension bridges［J］.Natural science journal of Xiangtan university，2011，33（1）：47-52.

［13］KIM H K，LEE M J，CHANG S P.Non-linear shapefinding analysis of a self-anchored suspension bridge［J］.Engineering structures，2002，24（12）：1547-1559.

［14］肖汝誠，項海帆.大跨徑懸索橋結(jié)構(gòu)分析理論及其專用程序系統(tǒng)的研究［J］.中國公路學報，1998，11（4）：42-50.

XIAO R C，XIANG H F.Research on the structural analysis theoryandspecial programof long-span suspension［J］.China journal of highway transport，1998，11（4）：42-50.

［15］邱順冬.橋梁工程軟件midas_civil應(yīng)用工程實例［M］.北京：人民交通出版社，2011.

本文編輯：苗變

Computation and Error Analysis of Main Cable Configuration of Large-Span Self-Anchored Suspension Bridges

YANG Heng，LI Yuansong*，QI Chao，ZHOU Xiaolong
School of Resource and Civil Engineering，Wuhan Institute of Technology，Wuhan 430074，China

Aimed at the change characteristics of the main cable configuration of self-anchored suspension bridge in construction process，we developed the program of segmental catenary method and parabola method to realize the accuracy computation of geometrical parameters of large-span self-anchored suspension bridges，such as main cable configuration，unloaded cable configuration，stress-free length of individual segments，hoist point coordinates，and pre-bias of cable saddle，etc.The results of the self-developed program and Midas civil software were compared with a project.The relative error of pre-bias of cable saddle is 11.3%by the segmental catenary method and is 23%by the parable method compared with that by the finite element method.We suggest that the results of segmental catenary method and parabola method should be modified and iterated repeatedly to determine the main configuration synthetically in construction process.

self-anchored suspension bridge；main cable configuration；segmental catenary method；parabola method

U448.25

A

10.3969/j.issn.1674-2869.2016.06.010

1674-2869（2016）06-0565-06

2016-05-12

武漢工程大學研究生創(chuàng)新基金（CX2015040）

楊恒，碩士研究生.E-mail：358031860@qq.com

*通訊作者：李元松，博士，教授.E-mail：li_yuan_song.@126.com