王海軍， 王子山， 劉紅釗， 魏 華， 吳德強(qiáng)

(1. 沈陽(yáng)工業(yè)大學(xué) 建筑與土木工程學(xué)院， 沈陽(yáng) 110870； 2. 中鐵十五局集團(tuán) 第五工程有限公司， 天津 300000)

ANSYS軟件提供了二次開(kāi)發(fā)的參數(shù)化語(yǔ)言，可以模擬拉索的單元主要有Link8、Link10與Link180[1].拉索是斜拉橋的關(guān)鍵構(gòu)件，對(duì)于拉索初始恒載索力的準(zhǔn)確施加成為后續(xù)斜拉橋分析的重要基礎(chǔ)[2-4].初始索力的施加可分為不同的方式，目前，ANSYS中初始索力的施加方法主要有三種：初應(yīng)變法、生死單元法與循環(huán)迭代法.

初應(yīng)變法適用于拉索點(diǎn)固定不變的情況，將拉索索力換算成為拉索初應(yīng)變，將初應(yīng)變賦予單元，相當(dāng)于給單元施加了已知的索力，此時(shí)索力是相當(dāng)準(zhǔn)確的.該方法的缺點(diǎn)是在幾何非線性的影響下，拉索錨固點(diǎn)必然會(huì)產(chǎn)生移動(dòng)，此時(shí)并不能得到準(zhǔn)確的拉索索力.

生死單元法通過(guò)對(duì)拉索節(jié)點(diǎn)施加集中力，利用軟件自帶的生死單元功能將其殺死，除去集中力后將其激活，此時(shí)便得到了初始索力，對(duì)多根索力的施加其原理一致[5].該方法的缺點(diǎn)是拉索索力在變形前后始終與拉索端點(diǎn)的連線一致，并不是方向不變，因此，該方法并不能得到準(zhǔn)確的初始索力.

循環(huán)迭代法分為初應(yīng)變迭代法和初始溫度迭代法.基于初應(yīng)變的迭代法將拉索索力換算為初應(yīng)變，在ANSYS中需要將其分解為豎向與橫向的初應(yīng)變，計(jì)算結(jié)果很準(zhǔn)確，該方法的缺點(diǎn)是計(jì)算量龐大且容易產(chǎn)生計(jì)算錯(cuò)誤.基于初始溫度的迭代法將拉索索力換算為初始溫度的形式進(jìn)行迭代計(jì)算，計(jì)算過(guò)程相對(duì)初應(yīng)變法簡(jiǎn)單，且不容易產(chǎn)生計(jì)算誤差，克服了上述方法中的不足，使計(jì)算結(jié)果更精確有效.

1 幾何非線性

(1)

式中：Eeq為Ernst修正后的彈性模量；E0為拉索的彈性模量；γ為單位體積斜拉索的重量；σ0為恒載作用下斜拉索的應(yīng)力；l為斜拉索水平投影的長(zhǎng)度.

本文假定設(shè)計(jì)院提供的設(shè)計(jì)斜拉索成橋索力為成橋恒載作用下的斜拉索索力，考慮斜拉橋的幾何非線性效應(yīng)，反復(fù)計(jì)算并調(diào)整施加斜拉索的初始索力值，直至在恒載作用下計(jì)算所得到的斜拉索索力值與設(shè)計(jì)提供的成橋索力值的相對(duì)誤差在工程允許的范圍內(nèi).實(shí)例證明，以設(shè)計(jì)院提供的斜拉橋恒載作用下的成橋索力通過(guò)溫度循環(huán)迭代進(jìn)行索力的迭代計(jì)算，可以得到較為精確的結(jié)果，加快成橋索力的準(zhǔn)確施加速度，并保證其在工程誤差5%的范圍內(nèi)，索力的計(jì)算公式為

(2)

式中：w為索力的迭代誤差；T迭代為經(jīng)過(guò)迭代計(jì)算后的索力；T設(shè)計(jì)為設(shè)計(jì)院提供的初始索力.

2 基于溫度循環(huán)迭代法索力的施加過(guò)程

利用ANSYS提供的參數(shù)化設(shè)計(jì)語(yǔ)言對(duì)其進(jìn)行二次開(kāi)發(fā)，在ANSYS中實(shí)現(xiàn)斜拉橋初始索力的確定方法[8-10].由于ANSYS不能直接輸入初始索力，因此，初始索力按照考慮垂度效應(yīng)Ernst公式修正彈性模量求出相應(yīng)溫度變化，以改變初始溫度的形式加到斜拉索上，流程圖如圖1所示.圖1中，T0為初始溫度；F0為目標(biāo)索力；F為提取索力；A為拉索斷面；α為溫度膨脹系數(shù).

圖1 溫度迭代流程圖Fig.1 Flow chart for temperature iteration

具體實(shí)現(xiàn)過(guò)程如下：

1) 在ANSYS通用軟件平臺(tái)上建立橋梁的有限元模型，定義施加重力場(chǎng)與二期恒載以及邊界約束條件.

2) 定義并輸入可變更參數(shù)數(shù)組，將初始索力對(duì)應(yīng)的初始溫度差導(dǎo)入軟件平臺(tái)中，作為可修改的參數(shù)存儲(chǔ)在數(shù)組內(nèi).

4) 利用ANSYS平臺(tái)的循環(huán)命令開(kāi)關(guān)進(jìn)行迭代求解計(jì)算，并將每一次求解的索力結(jié)果與提供的成橋恒載索力進(jìn)行誤差分析，將其誤差控制在5%的工程接受誤差范圍內(nèi).如果誤差超限，則把索力差值換算為溫差ΔT，令T=T+ΔT，繼續(xù)計(jì)算，直至計(jì)算索力滿足要求.

以通用軟件ANSYS平臺(tái)為基礎(chǔ)，每次經(jīng)過(guò)靜力計(jì)算，提取對(duì)應(yīng)斜拉索索力值，并且每次提取的索力值都會(huì)與設(shè)計(jì)院提供的初始索力值進(jìn)行對(duì)比，將每次提取得到的索力與初始索力的差值轉(zhuǎn)換為等大的索力值繼續(xù)施加到斜拉索上，通過(guò)循環(huán)計(jì)算的方式達(dá)到準(zhǔn)確施加初始索力的目的.

3 計(jì)算實(shí)例

贛江二橋工程起于S324省道(清宜公路)與工業(yè)四路交叉口，經(jīng)工業(yè)四路，在荷湖館張家村跨袁河合流處贛江北支，經(jīng)譽(yù)家洲后跨贛江南支，在傘頭村方家西側(cè)跨贛東大堤上岸，經(jīng)湖頭村，終點(diǎn)與藥都大道相交，路線全長(zhǎng)4.188 km，其中橋梁長(zhǎng) 2 529 m，主橋?yàn)殡p塔400 m跨斜拉橋.

主橋采用雙塔雙索面結(jié)合梁斜拉橋方案，邊跨設(shè)置一個(gè)輔助墩.主橋鋼結(jié)構(gòu)橋全長(zhǎng)(54+114)×2+400 m，梁底橫向水平設(shè)置，橋面設(shè)置雙向2%橫坡，道路中心線為橫向最高點(diǎn)，梁高3 m，橋面總寬29.2 m，為鋼框架結(jié)構(gòu).鋼橋總重6 920.2 t.斜拉索采用空間扇形雙索面布置，全橋共21×4×2=168根斜拉索，斜拉索采用PESC7-109、PESC7-121、PESC7-139、PESC7-151、PESC7-187、PESC7-199、PESC7-253共七種規(guī)格.斜拉索鋼絲的抗拉標(biāo)準(zhǔn)強(qiáng)度Rby=1 670 MPa，最長(zhǎng)索214.808 m(無(wú)應(yīng)力索長(zhǎng))，最短索55.656 m，單根索最大重量17.62 t(含PE和錨具重量)，順橋向梁上基本索距9 m，塔上索距1.3～2.2 m，采用扭絞型平行鋼絲拉索.塔端采用張拉端冷鑄錨具，梁端采用固定端冷鑄錨具.拉索布置圖如圖2所示(單位：m).

圖2 拉索布置Fig.2 Cable location

該斜拉橋?yàn)殡p塔雙索面結(jié)構(gòu)，復(fù)雜的結(jié)構(gòu)形式使得斜拉索之間相互影響程度加大，通過(guò)往常以手動(dòng)的形式來(lái)調(diào)整施加的索力，不但工作任務(wù)量大，而且還不一定能夠得到滿意的結(jié)果.由于斜拉橋結(jié)構(gòu)復(fù)雜，有限元模型規(guī)模巨大，累計(jì)節(jié)點(diǎn)數(shù)多達(dá)上萬(wàn)個(gè)，ANSYS完成一次響應(yīng)的索力施加計(jì)算往往就需要較長(zhǎng)時(shí)間，計(jì)算人員不能離開(kāi)，比較耗費(fèi)時(shí)間.按照本文提出的準(zhǔn)確施加初始索力的方法，利用APDL參數(shù)化設(shè)計(jì)語(yǔ)言進(jìn)行二次開(kāi)發(fā)，很方便地使ANSYS自動(dòng)完成該斜拉橋初始索力的施加過(guò)程.施加索力用到的1/2斜拉橋有限元模型與全橋有限元模型如圖3、4所示.

索力計(jì)算結(jié)果對(duì)比如表1所示.表1中，B為邊跨，Z為中跨.由于斜拉橋結(jié)構(gòu)復(fù)雜，斜拉索數(shù)量較多，表1中只給出了1/2斜拉索初始索力施加的迭代計(jì)算結(jié)果.ANSYS計(jì)算結(jié)果與設(shè)計(jì)索力的對(duì)比柱狀圖如圖5所示.圖5中拉索從左側(cè)邊跨到跨中的索號(hào)即為表1中索號(hào).

由表1和圖5可知，15次循環(huán)迭代計(jì)算后，計(jì)算所得到的索力值足夠接近設(shè)計(jì)院提供的設(shè)計(jì)索力值，并且相對(duì)誤差很小，所得到的成橋狀態(tài)時(shí)恒載作用下的迭代索力與設(shè)計(jì)院提供的設(shè)計(jì)成橋索力誤差在5%以內(nèi)，證明了本文方法具備迭代計(jì)算速度快、計(jì)算結(jié)果精度高、人工參與程度小的優(yōu)勢(shì).

表1 索力計(jì)算結(jié)果Tab.1 Calculation results of cable force

圖5 索力柱狀圖Fig.5 Histogram of cable force

4 結(jié) 論

本文通過(guò)ANSYS平臺(tái)的APDL參數(shù)化設(shè)計(jì)語(yǔ)言進(jìn)行二次開(kāi)發(fā)，實(shí)現(xiàn)了初始索力的迭代施加方法，通過(guò)利用本文初始索力迭代施加的計(jì)算方法，使斜拉橋繁瑣的調(diào)索過(guò)程在ANSYS軟件中自動(dòng)完成.算例表明，通過(guò)該方法可以方便快捷地迭代施加、計(jì)算斜拉橋拉索索力，從而使索力的施加變得快捷、準(zhǔn)確、有效.在參數(shù)化的分析過(guò)程中，可以修改其中的參數(shù)用來(lái)反復(fù)分析各種不同尺寸、不同恒荷載大小的多種斜拉索索力的施加，為以后其他斜拉橋初始索力的準(zhǔn)確施加提供了可利用的命令文件，為相關(guān)問(wèn)題的二次開(kāi)發(fā)提供了參考.

[1] 陳顏，高明大.基于ANSYS平臺(tái)的斜拉橋成橋索力優(yōu)化方法 [J].山西建筑，2011，36(16)：183-184.

(CHEN Yan，GAO Ming-da.Optimization method of cable tensions for stayedcable bridge based on ANSYS [J].Journal of Shanxi Architecture，2011，36(16)：183-184.)

[2] 王海軍，孔祥宇.格構(gòu)式鋼管混凝土風(fēng)力機(jī)塔架的設(shè)計(jì) [J].沈陽(yáng)工業(yè)大學(xué)學(xué)報(bào)，2015，37(3)：253-259.

(WANG Hai-jun，KONG Xiang-yu.Design of lattece wind turbine tower using concerte filled steel tube [J].Journal of Shenyang University of Technology，2015，37(3)：253-259.)

[3] 周強(qiáng)，楊文兵，楊新華.斜拉橋索力調(diào)整在ANSYS中的實(shí)現(xiàn) [J].華中科技大學(xué)學(xué)報(bào)(城市科學(xué)版)，2005，22(增刊)：81-83.

(ZHOU Qiang，YANG Wen-bing，YANG Xin-hua.The realization of cable force adjustment of cable stayed bridge in ANSYS [J].Journal of Huazhong University of Science and Technology，2005，22(Sup)：81-83.)

[4] 劉杰.曲線寬箱梁橋梁格模型建立與分析 [D].長(zhǎng)沙：中南林業(yè)科技大學(xué)，2014.

(LIU Jie.The grillage modeling and force analysis of wide curved box girder bridge [D].Changsha：Central South University of Forestry and Technology，2014.)

[5] 張楊永，周云崗，姜海西.基于ANSYS的超大跨度斜拉橋的索力模擬 [J].沈陽(yáng)建筑大學(xué)學(xué)報(bào)，2009，25(5)：909-913.

(ZHANG Yang-yong，ZHOU Yun-gang，JIANG Hai-xi.Cable force simulation for super long-span cable-stayed bridge based on ANSYS [J].Journal of Shenyang Jianzhu University，2009，25(5)：909-913.)

[6] 鄭瑩.自定義截面在ANSYS分析復(fù)雜鋼架結(jié)構(gòu)中的應(yīng)用 [J].化學(xué)工程與裝備，2015(9)：134-136.

(ZHENG Ying.Study on the technology of the composite board of the nickel iron chromium alloy and stainless steel [J].Chemical Engineering & Equipment，2015(9)：134-136.)

[7] 汪峰，劉沐宇.基于ANSYS平臺(tái)大跨度斜拉橋非線性成橋初始索力確定 [J].土木工程與管理學(xué)報(bào)，2013，30(2)：37-40.

(WANG Feng，LIU Mu-yu.Analysis of nonlinear bridge-completing cable force of long span cable-stayed bridge based on ANSYS [J].Journal of Civil Engineering and Management，2013，30(2)：37-40.)

[8] 朱云仙，真虹.國(guó)外并行工程研究與應(yīng)用進(jìn)展綜述 [J].機(jī)械設(shè)計(jì)與研究，2005，24(4)：66-69.

(ZHU Yun-xian，ZHEN Hong.Review on research and application of concurrent engineering in foreign countries [J].Mechanical Design and Research，2005，24(4)：66-69.)

[9] 李黎明.ANSYS 有限元分析適用教程 [M].北京：清華大學(xué)出版社，2005.

(LI Li-ming.Application tutorial of ANSYS finite element analysis [M].Beijing：Tsinghua University Press，2005.)

[10]金晶.斜拉橋靜動(dòng)力有限元分析 [D].合肥：合肥工業(yè)大學(xué)，2009.

(JIN Jing.Finite element analysis of the static and dynamic characteristics of cable-stayed bridge [D].Hefei：Hefei University of Technology，2009.)