懸索橋計(jì)算數(shù)值方法研究
2018-01-22
前言
對于大跨度懸索橋的施工控制需要確定懸索橋主纜的初始理想狀態(tài)以及成橋狀態(tài)，通常計(jì)算都采用有限元法和解析法，有限元法一般根據(jù)成橋的線形和受力情況，迭代出空纜狀態(tài)的線形和受力；解析法則根據(jù)成橋設(shè)計(jì)線形計(jì)算主纜無應(yīng)力長度，利用任何情況下主纜的無應(yīng)力長度不變的原理計(jì)算結(jié)構(gòu)參數(shù)，一般在解析法在數(shù)學(xué)方法上采用牛頓迭代法或擬牛頓法進(jìn)行計(jì)算。
本文結(jié)合主纜的實(shí)際情況：采用分段懸鏈線法計(jì)算，此法是考慮除主纜外的一期恒載及二期恒載作為多個(gè)集中力作用在各吊點(diǎn)處，主纜在各吊點(diǎn)之間線形為懸鏈線。并在計(jì)算結(jié)構(gòu)參數(shù)時(shí)考慮主纜的自重約束方程，即荷載集度在施工過程中不斷的變化，主纜總的質(zhì)量不變；根據(jù)主纜在主索鞍處的受力情況，以及中邊跨空纜與成橋狀態(tài)下無應(yīng)力長度相等；在數(shù)學(xué)方法上進(jìn)行改進(jìn)，采用精度更高、收斂更快、計(jì)算更穩(wěn)定的新Aitken迭代方法和同倫算法計(jì)算懸索橋的結(jié)構(gòu)參數(shù)。
1.分段懸鏈線的計(jì)算方法
1.1 基本假定：
?。?）主纜材料為線彈性，符合胡克定律；
（2）主纜是是理想柔性的，只能承受拉力，不能受壓，截面抗彎剛度對主纜線形影響忽略不計(jì)；
（3）忽略主纜橫截面在變形前后的變化；
通過以上假定，主纜的自重恒載集度沿主纜索長為常量，但變形前后可以不一樣。
1.2 分段懸鏈線法原理
考慮加勁梁的一期恒載及二期恒載作為多個(gè)集中力作用在各吊點(diǎn)處，此時(shí)主纜線形在各吊點(diǎn)之間為懸鏈線，取任意兩吊點(diǎn)間自由懸索建立坐標(biāo)系，以豎向?yàn)?方向，向下為正，水平向?yàn)?方向，向右為正，主纜上任意一點(diǎn)的拉格朗日坐標(biāo)為 ，對應(yīng)的笛卡兒坐標(biāo)為( , )，如圖1所示。
式中： 為成橋狀態(tài)主纜集度， 、 分別為塔頂主纜水平力和豎直力；
1.3 利用數(shù)學(xué)方法Aitken迭代計(jì)算豎向力
工程中常常會遇到許多非線性方程求根的問題，對于這一類的問題，一般不能用解析方法求得其解，而只能利用數(shù)值方法求得其近似解，目前常用的是牛頓迭代法，本文將Aitken迭代法結(jié)合實(shí)際工程，運(yùn)用其求解，與牛頓迭代法進(jìn)行比較。
本文采用了一種新的Aitken迭代法，新算法將二分法和迭代法結(jié)合起來，先用二分法預(yù)報(bào)初值，當(dāng)區(qū)間縮小到一定程度時(shí)再用改進(jìn)的Aitken算法迭代。
1.3.1 計(jì)算豎向力
假定索鞍水平力初值 、豎向力 和 可以通過式(2)求得 ，式(3)可求得 ，再由 通過式(2)求得 ，式(3)求得 ，按類似的方法進(jìn)行計(jì)算，一直計(jì)算到 ，并對 和 進(jìn)行修正，計(jì)算 。直到 , 為收斂精度，且同時(shí)滿足跨中斜率為0。
1.3.2Aitken迭代法求解
?。?）一般Aitken算法
1）簡單迭代格式：= ， ；
2）加速迭代格式： ；
如果 則將 賦值給 ，并修正 ( ＝ + ）；重復(fù)上面的式（1），（2）直到滿足精度要求，迭代停止，輸出 、 。
觀察圖2在相同初始迭代值、相同精度要求下，通過三種不同的數(shù)值迭代方法計(jì)算，新Aitken算法二階收斂，與牛頓法相比，新算法在迭代過程中，不需要計(jì)算導(dǎo)數(shù)值，減少了工作量，計(jì)算結(jié)果表明新的Aitken算法相比一般的Aitken算法和牛頓算法更快收斂，計(jì)算時(shí)可以節(jié)約計(jì)算時(shí)間。
2.空纜狀態(tài)下結(jié)構(gòu)參數(shù)計(jì)算
2.1 主索鞍預(yù)偏量計(jì)算
主索鞍預(yù)偏量的計(jì)算原則:保證各跨主纜無應(yīng)力索長空掛于主索鞍上，主纜在主索鞍槽內(nèi)不發(fā)生相對滑動。下面針對大跨度單跨不對稱懸索橋, 根據(jù)無應(yīng)力長度不變的原則列出下面的變形協(xié)調(diào)條件，自重約束條件和力學(xué)平衡條件:
?。?）幾何變形協(xié)調(diào)條件：中、邊跨主纜無應(yīng)力長度在任何階段、任何狀態(tài)都保持不變。
?。?）自重約束條件：主纜的重量在任何狀態(tài)下都保持不變。設(shè)主纜無應(yīng)力狀態(tài)時(shí)集度為 ，無應(yīng)力長度為 ，成橋狀態(tài)時(shí)集度為 ，主纜索長為 ，空纜狀態(tài)時(shí)集度為 ，主纜索長為 則有：
?。?） 力學(xué)平衡條件：空纜在索鞍兩側(cè)的水平力相等，設(shè)左索鞍預(yù)偏量為 ，右索鞍預(yù)偏量為 。
?。?）根據(jù)式（1）、（2）、（3）得到以下方程組：
式中： , 表示左、右邊跨空纜時(shí)索長， , 表示左、右邊跨成橋時(shí)索長， 表示左、右邊跨無應(yīng)力自重集度， 分別為中跨、左邊跨、右邊跨主纜在空纜、成橋時(shí)伸長量， 、 表示左、右邊跨成橋時(shí)自重集度。
式（4）中第一式和第二式可以采用新Aitken迭代法求解左、右邊跨主纜成橋集度，聯(lián)立第三式至第八式得到非線性方程組。工程計(jì)算中，經(jīng)常遇到需要求解非線性方程組的問題，非線性方程組的收斂速度及收斂性都比線性方程組要差。在求解非線性方程組時(shí)，牛頓迭代法是比較經(jīng)典的方法，其在局部收斂點(diǎn)附近是平方收斂的，但其解與初始解有關(guān)，且迭代過程中需要求導(dǎo)，計(jì)算量非常大且有時(shí)會出現(xiàn)計(jì)算困難。本文提出采用同倫算法，是基于其在大范圍收斂，并對初始值沒有嚴(yán)格限制，其思想是從容易求解的方程組開始，逐步過渡到原方程組的求解，最終求得方程組的近似解，下面簡單介紹一下同倫算法的相關(guān)內(nèi)容：
設(shè)非線性方程組為： ，其解為 。
?。?）構(gòu)造泛函 ：； 定義為：(其中： 為任意給的初值，假定為 函數(shù)（ )；
?。?）對于 的方程 ，當(dāng) 時(shí)， ； 是方程的解；當(dāng) 時(shí)， ； 是方程的解，即 = ；
?。?）基于這個(gè)思想最后得如下關(guān)系式： ( ，對初始值 )； 為雅可比矩陣，對 在 上積分，就可得到 = ；上面的非線性方程組問題就轉(zhuǎn)化為數(shù)值積分問題。
本文采用同倫算法求解非線性方程組，結(jié)合Matlab語言編寫求解程序，使用時(shí)只需修改相應(yīng)的雅可比矩陣和非線性方程組表達(dá)式，便可求解。
3.算例
本文編寫懸索橋通用計(jì)算程序，為說明本文程序正確性，下面通過算例進(jìn)行驗(yàn)算。
算例1. 江蘇江陰長江大橋?yàn)橹骺?385m的單跨鋼箱梁懸索橋，中跨主纜的設(shè)計(jì)參數(shù)：E=1.9×1011Pa， A=0.9027m2，無應(yīng)力狀態(tài)下的沿索長的均布荷載集度q0=77.70kN/m；邊跨主纜的設(shè)計(jì)參數(shù)：Eb=1.9×1011Pa，Ab =0.9526m2，q0=81.99kN/m；
算例2. 貴州某大橋?yàn)橹骺?36m的單跨簡支鋼桁梁懸索橋，結(jié)構(gòu)布置及幾何尺寸如圖3所示，成橋中跨主纜（單纜）的設(shè)計(jì)參數(shù)：E=2.0×1011Pa，A=0.16916m2，荷載集度q0=15.755kN/m；邊跨主纜荷載集度q0=15.01kN/m；
4.結(jié)語
?。?）在懸索橋的結(jié)構(gòu)參數(shù)計(jì)算中，本文考慮主纜自重約束方程，將主纜集度在不同的施工階段變化考慮進(jìn)去，計(jì)算更為精確。
?。?）新Aitken迭代在相同精度要求下能減少迭代次數(shù)，降低計(jì)算量，并能克服發(fā)散現(xiàn)象，收斂速度快、計(jì)算精度高，初值選擇范圍大。
?。?）同倫算法能解決擬牛頓法是在一定條件下計(jì)算時(shí)超收斂的，穩(wěn)定性差，迭代效果不理想的問題，其收斂范圍也比牛頓算法擴(kuò)大了，而且精確度也比原算法提高了。本文在解析迭代計(jì)算中采用新Aitken迭代法和同倫算法結(jié)合實(shí)際工程，運(yùn)用其求解，通過比較發(fā)現(xiàn)其是非常有效的求解數(shù)值方法。
