空間索面懸索橋空纜線形分析

張永水，吳章旭，馮 偉

(重慶交通大學 土木工程學院，重慶 400074)

空間索面懸索橋空纜線形分析

張永水，吳章旭，馮 偉

(重慶交通大學 土木工程學院，重慶 400074)

空間索面懸索橋的鞍座曲線對空纜線形分析有較大的影響，若還是采用虛焦點法分析空纜線形將會得到不準確的結(jié)果。采用兩種方法分析空纜線形的鞍座影響：解析迭代法在鞍座曲線上找切點分析空纜線形；基于ANSYS有限元軟件準確模擬主纜與鞍座的接觸關(guān)系分析空纜線形。兩種方法的計算結(jié)果完全一致，均可作為空間索面懸索橋空纜線形的計算方法。

橋梁工程；空間索面；鞍座曲線；空纜線形；解析法；有限元

0 引 言

空間索面懸索橋主纜在塔頂之間的橫向間距較小，因此主塔不必像平行索面懸索橋那樣設(shè)計成H型框架結(jié)構(gòu)，主塔的設(shè)計有了更大的自由度。設(shè)計時不僅能考慮結(jié)構(gòu)的安全可靠而且能夠與該地區(qū)的人文風俗、環(huán)境相融合，目前這一新的結(jié)構(gòu)體系，主要用于人行橋和小跨徑橋梁中。

主索鞍的鞍槽曲線設(shè)計成空間弧線主要是為了成橋狀態(tài)主纜的走向在鞍槽出口處與曲線相切，由于空間索面懸索橋具有立面矢跨比和平面矢跨比，因此主索鞍曲線在設(shè)計上一般采用平面圓弧和立面圓弧組合而成的曲線要素?；鶞仕鞴傻募茉O(shè)若還是采用MIDAS一次成橋倒拆模型的空纜線形再修正到基準索股存在以下問題：①索股的橫向受到主索鞍曲線的約束，無法考慮兩根基準索股的橫向間距和實際不符合；②虛焦點的計算方法[1]無法考慮主纜與鞍座的接觸關(guān)系；③主索鞍預偏量的計算只考慮虛焦點縱向位移量，再以中跨與邊跨縱向水平力是否相等作為判據(jù)，而沒有考慮虛焦點在高程上的變化。只有真實的模擬主索鞍曲線進行空纜線形分析才能夠解決MIDAS模型中存在的各種問題。采用精確模擬主索鞍曲線的計算方法分為兩個方向：一是采用解析法模擬主索鞍找切點[2-3]，二是采用有限元法模擬主索鞍與主纜的接觸關(guān)系[4-5]。

1 解析法模擬主索鞍

1.1 計算假定

1)主纜為柔性結(jié)構(gòu)只能受拉，且沒有抗彎剛度[6]；

2)主索鞍的空間圓弧曲線要求素立面圓弧和平面圓弧已知；

3)以立面圓弧最高點為界限，邊跨到散索點與中跨的無應力長度已知；

4)主索鞍切點間的曲線弧長近似等于索股的無應力長度。

1.2 計算方法

空纜線形分析采用經(jīng)典彈性懸鏈線理論，主纜的自重集度按主纜變形前的無應力長度分布。計算見圖1，其解析表達式見式(1)和式(2)[7]。中跨與邊跨的不平衡力決定了主索鞍圓心偏移量的大小，主索鞍的位移量與縱向不平衡力的調(diào)節(jié)由式(1)的微分增量關(guān)系進行迭代，能夠得到較快的收斂，見式(3)。

(1)

(2)

(3)

式中：H，V，dH分別為空纜水平力、豎向力、邊跨與中跨縱向不平衡力；q，s0，E，A0分別為空纜自重集度、計算無應力索長、彈性模量、截面面積；l，h，dl分別為計算跨徑，計算高差，索鞍圓心偏移量，空纜計算見圖1。

圖1 空纜計算Fig.1 Hollow cable calculation chart

主索鞍曲線立面圓心與水平圓心不一定在同一里程上，因此主要以立面圓弧圓心確定切點坐標，再由散索點的坐標來分析計算跨徑和計算高差，平面圓弧作為橫向偏距的計算。解析法計算程序在MATLABR2013b平臺編制腳本文件實現(xiàn)，計算步驟如下，解析迭代流程圖見圖2。

圖2 解析迭代計算流程Fig.2 Analytical iterative calculation flow chart

1)根據(jù)設(shè)計圖紙計算兩個塔的空間主索鞍曲線立面圓心坐標和水平面圓心坐標，假設(shè)中跨切點處縱向水平力Hc和Vc之間的夾角為α1，邊跨縱向水平力Hl和豎向力Vl的夾角為β1。

2)由主索鞍曲線要素計算中跨與邊跨的計算跨徑l，相對高差h以及切點間無應力索長s0。

4)判斷α2-α1與β2-β1的二范數(shù)是否小于收斂值eps1若小于eps1計算收斂，進入步驟5)若大于eps1則賦值α1=α2,β1=β2回到步驟2)。

6)根據(jù)dH=Hl-Hc的差值由式(3)計算主索鞍圓心坐標的偏移量得到兩主索鞍新的圓心坐標并進入步驟1)。

7)得到主索鞍圓心的偏移量，計算中跨和邊跨在不同溫度下里程對應的高程或無應力長度對應的里程和高程。

2 有限元法模擬主索鞍

2.1 計算假定

1)索股在主索鞍的索槽內(nèi)橫向受到剛性約束位移為0；

2)采用只受拉的桿單元能較準確地模擬主纜；

3)主索鞍的空間圓弧曲線要素平面圓弧和立面圓弧已知；

4)立面圓弧最高點為界限，邊跨到散索點與中跨的無應力長度已知。

2.2 計算方法

主索鞍建模以通用有限元分析軟件ANSYS作為計算平臺進行分析，主索鞍的空間曲線具有精確的解析表達式，因此在模型中能夠準確模擬主索鞍的位置和空間曲線長度。模擬主索鞍的懸索橋空纜線形分析主要采用了三類單元[8-9]：①只受拉的Link10桿單元模擬所有纜長；②只受壓的Link10桿單元模擬主索鞍與索股的接觸關(guān)系[10]；③采用MPC184剛性梁單元模擬主索鞍的底座。

建模過程需要注意以下幾點:

1)中跨與邊跨主纜的建模長度均以已知的無應力索長控制；

2)在主纜與鞍座可能接觸的地方鞍座桿單元需進行加密對應曲線長0.5～1cm可以較精確模擬主纜與鞍座的接觸關(guān)系；

3)施加位移邊界條件應包含散索點的回退量。

求解過程中同樣以邊跨與中跨的縱向水平力相等得到主索鞍圓心的偏移量，求解結(jié)束后更新模型的節(jié)點坐標即為空纜線形，見圖3。

圖3 有限元分析空纜線形Fig.3 Finite element analysis of hollow cable shape

3 算例分析

天河大橋北叉橋跨徑為100 m+266 m+100 m雙塔三跨的空間索面自錨式懸索橋。橋跨完全對稱布置，主索鞍空間曲線為平面圓弧和立面圓弧的組合形式，其構(gòu)造如圖4，立面圓弧頂與錨箱之間的偏距為13.9 m空間性非常明顯。利用文中解析法編制的腳本文件和ANSYS編寫的APDL宏文件進行模擬主索鞍的空纜線形分析，考慮主塔與主梁的壓縮量后計算分析所需的初始數(shù)據(jù)列于表1。

表1 初始參數(shù)

由初始參數(shù)分別用筆者介紹的解析法和有限元法得到主索鞍圓心的偏移量分別為48.4 m和48.2 m。為了對兩種方法的線形進行對比將有限元分析中的主索鞍圓心強制移動48.4 m，得到解析法與有限元立面投影線形對比見圖3，平面投影線形對比見圖4，相對誤差見圖5。從分析結(jié)果可以看出解析法和有限元得到的線形計算結(jié)果十分接近，相對誤差在1 mm內(nèi)，兩種方法的結(jié)果相互驗證。

圖4 主索鞍構(gòu)造Fig.4 Main cable saddle configuration

圖5 立面空纜線形Fig.5 Vertical hollow cable shape

圖6 水平面空纜線形Fig.6 Horizontal hollow cable shape

圖7 解析法與有限元誤差對比Fig.7 Error comparison between analytic method and finite element

4 結(jié) 語

基于彈性懸鏈線理論的解析法和基于ANSYS的有限單元法，同時模擬具有空間鞍座曲線的空纜架設(shè)線形分析，編制了相應的腳本文件和APDL宏文件。利用文中的程序?qū)μ旌哟髽螂p塔空間自錨式懸索橋基準索股架設(shè)進行分析，解析法和有限元法進行相互校核均能夠得到較高的計算精度，可作為類似空間索面懸索橋空纜線形計算分析的兩種參考方法。

[1] 黃強. 混凝土自錨式懸索橋施工控制研究[D]. 哈爾濱：哈爾濱工業(yè)大學, 2013.

HUANG Qiang.StudyonConstructionControlofConcreteSelf-anchorSuspensionBridge[D]. Harbin: Harbin Institute of Technology, 2013.

[2] 齊東春,沈銳利,郭永成. 懸索橋空間纜索線形的解析計算方法[J]. 武漢理工大學學報, 2013, 35(12): 109-113. QI Dongchun, SHEN Ruili, GUO Yongcheng. Analytical calculation method of main cable shape-finding of suspension bridge with spatial cables[J].JournalofWuhanUniversityofTechnology, 2013, 35(12):109-113.

[3] 羅喜恒. 懸索橋主纜線形的鞍座影響[J]. 公路交通科技, 2005, 22(8): 36-39. LUO Xiheng. Effect of saddle on cable shape of suspension bridges[J].JournalofHighwayandTransportonResearchandDevelopment, 2005, 22(8): 36-39.

[4] 潘永仁. 懸索橋結(jié)構(gòu)非線性分析理論與方法[M]. 北京：人民交通出版社, 2004. PAN Yongren.No-linerAnalysisTheoryMethodforSuspensionBridgeStructure[M]. Beijing: China Communications Press，2004.

[5] 羅喜恒,肖汝誠,項海帆. 用于懸索橋非線性分析的鞍座-索單元[J]. 土木工程學報, 2005, 38(6): 47-53. LUO Xiheng, XIAO Rucheng, XIANG Haifan. Saddle-cable elements for nonlinear analysis of suspension bridges[J].ChinaCivilEngineeringJournal, 2005, 38(6): 47-53.

[6] 占維. 懸索橋主纜施工計算的解析迭代方法[D].武漢：武漢理工大學2007. ZHAN Wei.AnalyticalIterationMethodforMainCableConstructionofSuspensionBridges[D]. Wuhan: Wuhan University of Technology, 2007.

[7] 唐茂林. 大跨度懸索橋空間幾何非線性分析與軟件開發(fā)[D]. 成都：西南交通大學, 2003. TANG Maolin. 3DGeometricNo-linearAnalysisofLong-spanSuspensionBridgeandItsSoftwareDevelopment[D].Chengdu: Southwest Jiaotong University, 2003.

[8] 王新敏，李義強，許宏偉. ANSYS 結(jié)構(gòu)分析單元與應用[M]. 北京：人民交通出版社, 2011. WANG Xinmin, LI Yiqiang, XU Hongwei.ANSYSStructuralAnalysisUnitandApplication[M]. Beijing: China Communications Press,2011.

[9] 王新敏. ANSYS 工程結(jié)構(gòu)數(shù)值分析[M]. 北京：人民交通出版社, 2007. WANG Xinmin.ANSYSNumericalAnalysisofEngineeringStructures[M]. Beijing: China Communications Press，2007.

[10] 許琳. 懸索橋施工過程中纜索系統(tǒng)的精細模擬[D]. 大連：大連理工大學, 2013. XU Lin.AccurateSimulationofCableSysteminConstructionProcessoftheSuspensionBridge[D]. Dalian: Dalian University of Technology,2013.

(責任編輯 朱漢容)

Hollow Cable Shape Analysis of Suspension Bridge with Spatial Cable Plane

ZHANG Yongshui, WU Zhangxu, FENG Wei

(School of Civil Engineering, Chongqing Jiaotong University, Chongqing 400074, P.R.China)

Spatial cables suspension with saddle curve has a great effect on the hollow cable shape analysis, while still analysis hollow cable shape is obtained by using the virtual focus method will be inaccurate results. Two kinds of methods were adopted to analyze the hollow cable shape with saddle effect:analytic iteration method on the saddle curve wasused to find point of tangency in analysis empty cable shape. Based on the ANSYS finite element software accurate simulation of the contact relationship of main cable and saddle was made to analysis of hollow cable shape. The calculation results of the two methods were completely consistent, both can be used as suspension hollow cable shape with spatial cables calculation method.

bridge engineering; spatial cable plane; saddle curve; hollow cable shape; analytical method; finite element

10.3969/j.issn.1674- 0696.2017.02.02

2016- 01- 03；

2016- 03-12

重慶市研究生科研創(chuàng)新項目(CYS16184)

張永水(1965—)，男，山西臨汾人，教授，主要從事橋梁工程控制理論與智能監(jiān)控技術(shù)工作。Email:zhang_yongshui@163.com。

吳章旭(1992—)，男，重慶人，碩士研究生，主要從事橋梁工程控制理論與智能監(jiān)控技術(shù)工作。Email:13883200742@163.com。

U448.25

A

1674- 0696(2017)02- 007- 04