摘要 懸索橋主纜計算理論有很多，但是傳統(tǒng)的方法存在諸多的缺點，特別是懸索橋主纜成橋線形分析求解過程收斂性問題，是傳統(tǒng)方法所必須面對的問題。為解決懸索橋主纜成橋線形分析求解過程收斂性問題，文章研究了外荷載與主纜線形的幾何關(guān)系，并構(gòu)造出基于外荷載—線形關(guān)系的解析法基本方程，最后通過算例驗證了新方法的可行性、有效性及精度。所提出的方法對初值要求低，并且能快速收斂得到主纜線形精確解，適用于懸索橋主纜成橋線形分析。

關(guān)鍵詞 懸索橋；主纜線形；解析法；找形

中圖分類號 U412.352.1 文獻(xiàn)標(biāo)識碼 A 文章編號 2096-8949（2024）04-0041-03

0 引言

目前懸索橋主纜計算理論主要包括有限元法、解析法。

有限元法發(fā)展歷史較為久遠(yuǎn)，1966年Brotton首次發(fā)表一種以矩陣位移法求解懸索橋結(jié)構(gòu)的通用分析方法^[1]。后續(xù)很多學(xué)者圍繞著幾何非線形梁單元剛度矩陣做了許多工作，如陳政清等的空間桿系UL行列式法等^[2]。

解析法最早是基于沿跨徑水平方向分布荷載假定的傳統(tǒng)拋物線理論，隨著跨度的增加，恒載沿跨徑方向分布的不均勻性趨于明顯，很多專家深入分析研究，逐步發(fā)展并形成了基于分段懸鏈線的解析法^[3]。分段懸鏈線法也納入了現(xiàn)行的懸索橋規(guī)范^[4]。

基于分段懸鏈線的解析法在研究過程中逐步與有限元法有機結(jié)合，大多以分段懸鏈線基本方程推導(dǎo)得到索單元剛度矩陣，納入當(dāng)前主流懸索橋計算軟件。

基于分段懸鏈線的數(shù)值解析法參數(shù)少，收斂速度快，是目前計算精度最高、應(yīng)用最普遍的方法，但是這類方法原理上需要引入包含索段無應(yīng)力長度及變形相容條件的索段狀態(tài)方程，由于纜索結(jié)構(gòu)強烈的非線性使得初值的選擇對迭代收斂性及收斂速度影響較大。

為了克服分段懸鏈線解析法的初值要求苛刻的不足，多位專家圍繞著力素與主纜線形幾何關(guān)系開展研究：孫遠(yuǎn)等^[5]提出“合理軸線”思路，基于桁架單元和結(jié)構(gòu)力學(xué)的基本原理直接建立單元坐標(biāo)與桿端力的關(guān)系式，劉超等^[6]提出基于整體力學(xué)分析的主纜找形方法，以理想柔索為假定，對主纜控制點取矩構(gòu)成基本方程，并求解支點反力，再通過迭代求解出主纜線形坐標(biāo)。這些思路另辟蹊徑，不再使用分段懸鏈線解析法的基本方程，避開了容易導(dǎo)致迭代發(fā)散的無應(yīng)力長度及變形相容條件，大大降低了對迭代初值的要求。

該文提出一種基于外荷載與主纜幾何線形的關(guān)系，用數(shù)值解析法迭代求解主纜線形的新方法，旨在降低初值標(biāo)準(zhǔn)，提高算法普適性，并能得到主纜線形的精確解。

1 基于外荷載—線形關(guān)系的解析法

1.1 基本假定

懸索橋主纜線形計算基本假定如下：

（1） 主纜為理想柔索，只能受拉，不考慮其抗彎剛度及抗扭剛度。

（2） 主纜的應(yīng)力—應(yīng)變符合胡克定律。

（3） 受力前后主纜的橫截面面積及抗拉強度均保持不變。

（4） 集中荷載均作用在單元節(jié)點上。

如圖1所示，以主纜地錨點為原點O，主纜塔頂錨點為A，中點為C點，B（x，z）點為主纜上任意點，主纜跨徑為L，兩端點Z方向高差為h，跨中垂度為f。

不失一般性，假設(shè)主纜上作用沿主纜索長方向分布的纜索自重等分布荷載q_s（x）和若干集中荷載P（x）如吊桿荷載等，考慮可能存在的斜吊桿，集中荷載P（x）可分解為水平方向的P_x（x）和豎直方向的P_z（x），地錨點O支點反力為F_x和F_z，分布荷載、集中力及支點反力方向均按與坐標(biāo)軸同向為正，反向為負(fù)。

1.2 外荷載節(jié)點力素

為方便表述，引入外荷載節(jié)點力素定義如下：主纜上任意節(jié)點B左側(cè)外荷載在B點產(chǎn)生的節(jié)點力素定義為外荷載節(jié)點力素。外荷載節(jié)點力素包含外荷載節(jié)點力和外荷載節(jié)點力矩，分別用函數(shù)F（x）和M（x）表示。

如圖2所示，取任意主纜單元，單元左右節(jié)點分別為i和j，已知左節(jié)點外荷載節(jié)點力素為F_x（x_i）、F_z（x_i） 和M（x_i），單元上作用分布荷載q（x）和集中荷載Px（xj）、Pz（xj）。

外荷載節(jié)點力：

外荷載節(jié)點力矩迭代表達(dá)式為：

與常見的主纜線形顯式方程不同，基于外荷載—線形關(guān)系的解析法實際是采用外荷載節(jié)點力素方程（1）（2）（3）（6）（7）和線形方程（12）以迭代方式構(gòu)建的隱式表達(dá)式。并通過引入“外荷載節(jié)點力素”概念，消除了對荷載分布形式的限定，應(yīng)用范圍擴展為能夠求得“外荷載節(jié)點力素”的任意荷載分布形式都可以應(yīng)用公式求解線形。

常規(guī)的無集中荷載拋物線法可以認(rèn)為是基于外荷載—線形關(guān)系的解析法在跨徑方向均布荷載且無吊桿張力等集中荷載作用下的特例，均布荷載彎矩代入（12）可得到常見的主纜拋物線方程解析式：

主纜無應(yīng)力長度S₀計算公式如下：

式中，T_i——任意主纜單元的張力；L_i——主纜單元長度；S——主纜的形狀長度。

2 算例

為方便對比，算例采用張志國等^[7]和沈銳利^[8]獨立計算過的某兩支點等高懸索橋，具體參數(shù)如下：

懸索橋跨度L=888 m，吊索間距12 m，主纜恒載集度q=54 kN/m，加勁梁等其余恒載集度W=200 kN/m，主纜截面積A=0.6 m²，索材彈性模量E=2.0×10⁵MPa，跨中矢高f=60 m，70 m，80 m，90 m，100 m。應(yīng)用該文算法迭代計算了主纜水平張力H，主纜線形坐標(biāo)及無應(yīng)力長度，列表對比如表1～3。

由表1～3可知，該文算法的計算結(jié)果與文獻(xiàn)^[7-8]基本一致，說明該文提出來的基于外荷載—線形關(guān)系的解析法正確可行。

3 結(jié)論

該文提出了一種求解懸索橋主纜線形的新方法即基于外荷載與主纜線形關(guān)系的解析法。該方法以隱函數(shù)及迭代方式構(gòu)造線形基本方程，力學(xué)概念清晰，迭代過程簡單，且計算得到的主纜線形坐標(biāo)均為精確解。

該文的方法通過后置無應(yīng)力索長及變形相容條件求解，簡化了求解難度，降低了迭代初值要求，提高了算法普適性。

參考文獻(xiàn)

[1]D. M. Brotton. A General Computer Program for the Solution of Suspension Bridge Problems. Lisbon， 1966， 27-49.

[2]陳政清， 曾慶元， 顏全勝. 空間桿系結(jié)構(gòu)大撓度問題內(nèi)力分析的UL列式法[J]. 土木工程學(xué)報， 1992（5）： 34-44.

[3]李傳習(xí). 混合梁懸索橋非線性精細(xì)計算理論及其應(yīng)用[D]. 長沙：湖南大學(xué)， 2006.

[4]公路懸索橋設(shè)計規(guī)范：JTG/T D65-05—2015[S]. 北京：人民交通出版社， 2015.

[5]孫遠(yuǎn)， 羅文孝， 劉梅. 基于有限元列式的懸索橋主纜找形算法[J]. 土木工程與管理學(xué)報， 2018（2）： 116-123.

[6]劉超， 高展. 基于整體力學(xué)分析的自錨式懸索橋主纜找形法[J]. 同濟大學(xué)學(xué)報（自然科學(xué)版）， 2020（1）： 1-6.

[7]張志國， 鄒振祝， 趙玉成， 等．懸索橋主纜線形解析方程解及應(yīng)用[J]. 工程力學(xué)， 2005（3）： 172-176.

[8]沈銳利. 懸索橋主纜系統(tǒng)設(shè)計及架設(shè)計算方法研究[J]. 土木工程學(xué)報， 1996（2）： 3-9.

收稿日期：2023-12-07

作者簡介：楊勇智（1976—），男，本科，高級工程師，研究方向：橋梁結(jié)構(gòu)設(shè)計。