混凝土斜拉橋的優(yōu)化設(shè)計(jì)
2018-04-16
混凝土斜拉橋優(yōu)化被歸結(jié)為多目標(biāo)優(yōu)化問(wèn)題，其目的是獲得最低成本，最小偏移和應(yīng)力，并尋求帕累托解決方案。通過(guò)凸標(biāo)量函數(shù)的最小化發(fā)現(xiàn)了該解決方案。數(shù)值方法允許纜索安裝力、纜索調(diào)整力和斜拉索區(qū)域、橋面和橋塔橫截面的計(jì)算，以滿足整個(gè)結(jié)構(gòu)位移和應(yīng)力，無(wú)論是在架設(shè)期間或竣工。因此本文所提出的優(yōu)化算法可以被用來(lái)研究拓?fù)浜蛶缀卧O(shè)計(jì)變量的效應(yīng)。
在MATLAB開(kāi)發(fā)的計(jì)算機(jī)程序被用于結(jié)構(gòu)分析、靈敏性分析和優(yōu)化。結(jié)構(gòu)分析采用一種有限元模型，其包括由于混凝土施工順序、幾何非線性和時(shí)間依存性效應(yīng)造成的負(fù)載歷史和幾何變化。施工過(guò)程和混凝土流變行為顯著影響了混凝土斜拉橋的應(yīng)力和變形。在分析中，同樣應(yīng)考慮處理纜索或大型靈活結(jié)構(gòu)時(shí)產(chǎn)生的幾何非線性情況。
1.時(shí)間依存性效應(yīng)建模
在這項(xiàng)研究中，根據(jù)歐洲規(guī)范2公式進(jìn)行評(píng)估混凝土的蠕變、收縮和老化的時(shí)間依存賴性效應(yīng)。蠕變模型基于線性粘彈性及考慮老化影響。收縮應(yīng)變依存于時(shí)間，但與應(yīng)力無(wú)關(guān)。由作者在最近的研究工作中提出了有關(guān)時(shí)間依存性效應(yīng)建模的詳細(xì)考慮。
1. 1混凝土老化
由于固化的結(jié)果，混凝土的強(qiáng)度和彈性系數(shù)隨時(shí)間增加。在早期，強(qiáng)度和彈性系數(shù)迅速增加，然后增加逐漸停滯，但不會(huì)完全停止。根據(jù)EN1992-1-1（2010），天數(shù)t，混凝土的彈性系數(shù)由下式給出（1）。
其中，Ecm是混凝土彈性系數(shù)在28天內(nèi)的平均值，t是混凝土固化天數(shù)，s是取決于水泥類型的系數(shù)。
1.2混凝土蠕變
在時(shí)刻t0施加單軸應(yīng)力σc的混凝土試件，其在時(shí)刻t的總應(yīng)變可以寫(xiě)成相關(guān)應(yīng)力？ca（t，t0）和與無(wú)關(guān)應(yīng)力？cn（t）的總和，應(yīng)變?yōu)椋?）。
其中J（t，t0）是蠕變功能；如果應(yīng)力小于45%的混凝土抗壓強(qiáng)度標(biāo)準(zhǔn)值，則疊加原理是有效的，并且蠕變應(yīng)變隨著所施加的應(yīng)力而線性變化。
在施工階段和結(jié)構(gòu)使用壽命期間，斜拉橋的應(yīng)力不斷變化。根據(jù)變量的應(yīng)力以及使用疊加原理，方程（2）可改寫(xiě)為（3）。
已經(jīng)提出了幾種方法來(lái)解決這個(gè)方程，包括簡(jiǎn)化方法，逐步數(shù)值積分以及蠕變函數(shù)的粗略估算。在這篇文章中，由Dirichlet級(jí)數(shù)（Bazant，1988）概略估算該蠕變函數(shù)，從而得到（4）。
其中，n是Dirichlet級(jí)數(shù)的項(xiàng)數(shù)，系數(shù)aj通過(guò)使用最小二乘方法從曲線擬合得出。系數(shù)1/aj被稱為延遲時(shí)間，并被選擇來(lái)覆蓋用于計(jì)算蠕變系數(shù)的時(shí)間值范圍以。
1.3混凝土收縮
根據(jù)EN1992-1-1（2010），在時(shí)刻t的總收縮應(yīng)變？cs（t），是自收縮（？ca）和干燥收縮（？cd）的總和。時(shí)刻t的干燥收縮被定義為（5）。
其中，系數(shù)βds（t，ts）和kh取決于構(gòu)件理論尺寸和混凝土干燥收縮開(kāi)始時(shí)刻。參數(shù)εcd（t）取決于環(huán)境相對(duì)濕度、水泥類型和混凝土抗壓強(qiáng)度。
自收縮是混凝土早期硬化期間化學(xué)反應(yīng)的結(jié)果，它可以通過(guò)時(shí)刻t由以下公式表示（6）。
其中，εca（∞）是自收縮應(yīng)力的長(zhǎng)期值，βas（t）是自收縮隨時(shí)間演變的函數(shù)。
1.4時(shí)間依存性效應(yīng)模擬
在結(jié)構(gòu)分析中，通過(guò)產(chǎn)生相同位移場(chǎng)的等效節(jié)點(diǎn)力模擬時(shí)間依存效應(yīng)，并將其作為時(shí)間依存性效應(yīng)及從中計(jì)算出實(shí)際變形狀態(tài)。每個(gè)時(shí)間間隔，使用有限元公式和相應(yīng)的蠕變值以及根據(jù)先前提出公式計(jì)算出的收縮應(yīng)變，將這些計(jì)算的應(yīng)力作為初始變形。于是，僅用應(yīng)力和機(jī)械原點(diǎn)變形之間的彈性本構(gòu)關(guān)系計(jì)算該應(yīng)力。
1.5幾何非線性效應(yīng)
斜拉橋中幾何非線性有三個(gè)主要來(lái)源：由于自身重量造成下垂的非線性軸向力-延伸關(guān)聯(lián)斜拉索；根據(jù)聯(lián)合彎曲和軸向力，非線性軸向力和彎矩-變形關(guān)聯(lián)橋塔和橋面；以及大位移引起的幾何變化。在這篇文章中，通過(guò)二階彈性分析方法考慮幾何非線性效應(yīng)。
考慮斜拉索中幾何非線性的一種廣泛使用方法是用于考慮能描述纜索懸索效應(yīng)并具有等效彈性系數(shù)（Ernst 1965年）的等效直弦構(gòu)件。
纜索等效彈性系數(shù)值由下式給出（7）。
其中，Eeq是等效的纜索彈性系數(shù)，E是等效的纜索材料彈性系數(shù)，y是纜索材料的比重，L是弦的長(zhǎng)度，α是索弦與水平方向之間的角度，σ是纜索的張力應(yīng)力。
二階效應(yīng)被認(rèn)為是采用等效側(cè)力法，也被稱為虛擬側(cè)向荷載法或P-A迭代方法，通常用于建筑結(jié)構(gòu)二階分析。
1.6包括架設(shè)階段的結(jié)構(gòu)分析
繼施工階段后，使用正向分析程序，進(jìn)行架設(shè)階段的建模和分析，然后獲得整個(gè)施工的應(yīng)力和位移，使其能直接考慮時(shí)間依存性效應(yīng)。這適合于結(jié)構(gòu)優(yōu)化目的，因?yàn)樵诟麟A段結(jié)束時(shí)和要求優(yōu)化模塊之前，可同時(shí)使用所有有關(guān)應(yīng)力、位移及其靈敏度的信息。
在斜拉橋施工中，采用最普遍和廣泛使用的方法――平衡懸臂施工法進(jìn)行架設(shè)。施工一開(kāi)始，要建造橋塔，并在橋塔的兩側(cè)開(kāi)始懸臂施工。在隨后的階段，架設(shè)其他橋面節(jié)段和斜拉索，直到橋面合攏。
使用一個(gè)有限元計(jì)算機(jī)程序進(jìn)行結(jié)構(gòu)分析，該橋梁被構(gòu)建成一個(gè)二維構(gòu)架結(jié)構(gòu)。該模型不考慮橋面扭轉(zhuǎn)和結(jié)構(gòu)的三維特性。然而，由于該研究以架設(shè)和使用情況為中心，因而利用二維模型是足夠用的。 2.優(yōu)化設(shè)計(jì)公式
在斜拉橋的優(yōu)化設(shè)計(jì)中，由于設(shè)計(jì)變量數(shù)量高且設(shè)計(jì)目標(biāo)非線性特征也會(huì)產(chǎn)生沖突，因此搜索空間很復(fù)雜。這里作為一個(gè)多目標(biāo)優(yōu)化問(wèn)題被提出，從中獲得帕累托最優(yōu)解向量。這意味著在不增加至少一個(gè)目標(biāo)的情況下，不存在其他可行的能降低一個(gè)目標(biāo)的向量。它涉及設(shè)計(jì)變量、設(shè)計(jì)目標(biāo)和目標(biāo)函數(shù)的定義。
2. 1設(shè)計(jì)變量
考慮的設(shè)計(jì)變量是斜拉索區(qū)域、預(yù)應(yīng)力、橋面和橋塔的橫截面尺寸。設(shè)計(jì)變量用xi來(lái)表示，而全球設(shè)計(jì)變量向量是（8）。
斜拉索由0.6英寸直徑的鋼索制成（15.7mm的公稱直徑和1.5cm2的橫截面積）。對(duì)于橋塔的橫截面，考慮使用矩形空心型鋼，對(duì)于橋面考慮使用三種橫截面類型：梁板、單室箱和三室箱。橋面和橋塔的橫截面設(shè)計(jì)變量對(duì)減重（或降低成本）有直接的影響。斜拉索區(qū)域和斜拉索索力在整個(gè)結(jié)構(gòu)的應(yīng)力分布中發(fā)揮極為重要的作用，因?yàn)樾崩鲄^(qū)域和斜拉索索力對(duì)橋面的梁型特性的范圍作出了定義。此外，斜拉索區(qū)域和斜拉索索力是調(diào)整橋梁幾何形狀和擾度控制的基礎(chǔ)，否則可以只通過(guò)橋面的嚴(yán)格加勁來(lái)實(shí)現(xiàn)，這與預(yù)期的減少材料相反。
2.2設(shè)計(jì)目標(biāo)
斜拉橋的設(shè)計(jì)涉及到實(shí)現(xiàn)一些設(shè)計(jì)目標(biāo)，以檢驗(yàn)使用和強(qiáng)度標(biāo)準(zhǔn)。該目標(biāo)應(yīng)以規(guī)范化的形式來(lái)配置。這些目標(biāo)會(huì)因處于施工階段及恒載和活荷載下的整個(gè)橋梁位移和應(yīng)力極限值而產(chǎn)生。該設(shè)計(jì)還應(yīng)力求將該結(jié)構(gòu)的成本降至最低。考慮到這一點(diǎn)，第一個(gè)目標(biāo)可表示為（9）。
其中，C是該結(jié)構(gòu)的現(xiàn)時(shí)成本，C0為參比成本，這對(duì)應(yīng)于每個(gè)分析和優(yōu)化周期的初始成本。這確保了在每個(gè)周期，成本始終是優(yōu)化算法的主要目標(biāo)之一。該結(jié)構(gòu)的成本被當(dāng)做材料成本（混凝土、加強(qiáng)鋼和預(yù)應(yīng)力鋼）。這些材料的單價(jià)可以通過(guò)咨詢葡萄牙供應(yīng)商公司而獲得。第二組目標(biāo)產(chǎn)生于限制橋面的垂直位移和橋塔的水平位移，以獲得要求的最終橋面縱剖面，并使橋塔彎曲擾度最小化（10）。
其中，δ和δ0分別是位移值和控制下的位移限值。
第三組的目標(biāo)會(huì)因處于施工階段和恒載下的整個(gè)橋梁橋面位移和橋塔應(yīng)力極限值而產(chǎn)生。這些目標(biāo)都與使用情況有關(guān)。
根據(jù)歐洲規(guī)范2（EN1992-1-1 2010）的建議，混凝土的抗壓應(yīng)力限定為45%的混凝土抗壓強(qiáng)度（fck）的特征值，從而使混凝土保持在線性蠕變范圍內(nèi)，并預(yù)防縱向開(kāi)裂?；炷晾瓚?yīng)力限定為5%的混凝土典型軸向抗拉強(qiáng)度的分位數(shù)（fctk，0.05），以避免開(kāi)裂，確保耐久性，得（11）、（12）。
其中，σc是在混凝土構(gòu)件的作用應(yīng)力。使用抗拉和抗壓應(yīng)力及相應(yīng)的容許應(yīng)力時(shí)，要考慮相應(yīng)的信號(hào)。
也應(yīng)對(duì)混凝土構(gòu)件的最大應(yīng)力進(jìn)行檢查。該目標(biāo)可以表示為（13）。
其中，σc是在混凝土構(gòu)件的作用應(yīng)力，σallow是抗拉或抗壓中相應(yīng)的容許應(yīng)力。作用應(yīng)力通過(guò)作用軸向力（NEd）和彎矩（MEd）計(jì)算得出。容許值被確定為一個(gè)應(yīng)力，等同于橫截面組合的軸向力-彎矩設(shè)計(jì)阻力（NRd；MRd）。此值通過(guò)混凝土構(gòu)件相應(yīng)橫截面產(chǎn)生的無(wú)量綱交互圖獲得。對(duì)于這些計(jì)算結(jié)果，加強(qiáng)鋼區(qū)域不是一個(gè)設(shè)計(jì)變量，被認(rèn)為是混凝土截面積的 2%。這被用來(lái)作為混凝土構(gòu)件的加強(qiáng)鋼區(qū)域的平均值，表示普遍的實(shí)用價(jià)值。
斜拉索中應(yīng)力的剩余目標(biāo)為（14）、（15）。
其中，σ和fpk分別是斜拉索中的作用應(yīng)力和預(yù)應(yīng)力鋼抗拉強(qiáng)度的特征值。如果作用應(yīng)力大于0.1fpk，則方程式（14）適用；如果作用應(yīng)力小于或等于0.1fpk，則方程式（15）適用。斜拉索拉伸應(yīng)力的0.1fpk下限值被考慮，以確保其結(jié)構(gòu)效能。方程式（14）中的k值等于架設(shè)期間的0.55，使用情況下的0.50和強(qiáng)度校核的0.74。
2.3目標(biāo)函數(shù)
多目標(biāo)優(yōu)化的目標(biāo)是最大限度地減少該組所有目標(biāo)的設(shè)計(jì)變量。這通過(guò)最小最大值優(yōu)化來(lái)實(shí)現(xiàn)，得（16）、（17）。
該問(wèn)題是不連續(xù)、不可微的問(wèn)題，因此不好解決。然而，如Simoes和Templeman所述（1989），可以通過(guò)一個(gè)不受約束的凸標(biāo)量最小函數(shù)來(lái)間接獲得該求解方法，結(jié)果是Kreisselmeier-Steinhauser函數(shù)（17）。這是連續(xù)、可微的問(wèn)題，因此很容易解決。該函數(shù)只取決于一個(gè)控制參數(shù)p，它不能被減小。
我們的目標(biāo)函數(shù)gj（x）沒(méi)有一個(gè)明確的代數(shù)形式，只能通過(guò)一個(gè)特定的設(shè)計(jì)變量向量的結(jié)構(gòu)分析結(jié)果得出數(shù)值。采取的策略是通過(guò)明確的近似模型迭代序列手段，來(lái)解方程式（17）。一個(gè)明確的近似值可以通過(guò)線性項(xiàng)后截?cái)嗟乃心繕?biāo)函數(shù)gj（x）的泰勒級(jí)數(shù)展開(kāi)式制定。該算法得（18）。
其中，N和M分別是設(shè)計(jì)變量數(shù)量和目標(biāo)數(shù)量。goj（x）和dg0j（x）/dxi是根據(jù)當(dāng)前設(shè)計(jì)變量向量（x0）計(jì)算的目標(biāo)及其靈敏度，其中制定了泰勒級(jí)數(shù)展開(kāi)式。g0j（x）特定數(shù)值求解方程（18）僅構(gòu)成該問(wèn)題完整求解方法的一個(gè)迭代法。使用MATLAB函數(shù)fmincon，進(jìn)行目標(biāo)函數(shù)的最小值計(jì)算，采用一系列二次問(wèn)題，使處于約束范圍內(nèi)的各變量標(biāo)量函數(shù)最小化。
該算法與多個(gè)起點(diǎn)相關(guān)聯(lián)，提供一個(gè)優(yōu)化的工程求解方法，通過(guò)重新排列斜拉索的支承剛度及橋面和橋塔之間的質(zhì)量分布變得合理。目標(biāo)函數(shù)使過(guò)程從一個(gè)可行或不可行的求解方法開(kāi)始。
2.4敏感性分析
考慮到源代碼的可用性和大量的目標(biāo)（應(yīng)力和位移）與設(shè)計(jì)變量數(shù)量相對(duì)，則采用分析直接離散方法，進(jìn)行靈敏度分析。相對(duì)于設(shè)計(jì)變量的位移靈敏度通過(guò)區(qū)分平衡方程獲得（19）～（21）。
可在形式上被重寫(xiě)，其中，Qvi是相對(duì)于ith設(shè)計(jì)變量的系統(tǒng)虛擬假載向量。位移靈敏度可表示為（22）。其中需要存儲(chǔ)剛度矩陣，給剛度矩陣和右側(cè)導(dǎo)數(shù)預(yù)編程序，以使位移導(dǎo)數(shù)可以通過(guò)N假載右側(cè)求解方法來(lái)計(jì)算。應(yīng)力靈敏度通過(guò)有限元應(yīng)力位移關(guān)系的鏈求導(dǎo)來(lái)確定（22）、（23）。 在對(duì)導(dǎo)數(shù)表達(dá)式預(yù)編程序并請(qǐng)求這一階段的條件下，右側(cè)的第一項(xiàng)可在全球系統(tǒng)元素貢獻(xiàn)計(jì)算期間直接計(jì)算出來(lái)。由于位移導(dǎo)數(shù)已知，則很容易計(jì)算出來(lái)右側(cè)第二項(xiàng)。
3.數(shù)值例子和結(jié)果
3. 1數(shù)值模型的說(shuō)明
數(shù)值模型涉及一個(gè)總跨徑為284m，中跨距為148m的對(duì)稱混凝土斜拉橋。橋塔的總高度為52m，橋面置于基礎(chǔ)以上15m的位置。圖1說(shuō)明了作為示例橋的幾何結(jié)構(gòu)。
使用平衡懸臂施工法的按工序劃分施工階段。有關(guān)完整的橋，橋面到橋塔連接僅僅是垂直位置的連接，但要保證結(jié)構(gòu)穩(wěn)定，在假設(shè)階段，該連接固定。
在施工階段考慮采取的行動(dòng)包括橋面自重和1.0kN/m2（這取決于人員及手動(dòng)工具），0.5kN/m2（這取決于臨時(shí)設(shè)備），和懸臂邊緣上400kN的集中載荷（EN1991-1-6 2005）。
在第1階段，在橋塔各側(cè)執(zhí)行第一橋面節(jié)段，并且安裝第一對(duì)斜拉索。在以后的每個(gè)階段，在橋塔各側(cè)上對(duì)稱地增加一個(gè)橋面節(jié)段和一個(gè)斜拉索。在第9階段，在橫向跨距末端處增加一個(gè)簡(jiǎn)單的支座；而第10階段，在中跨處增加控制彎矩和軸向力的連接件，以模擬橋梁合攏。
在兩個(gè)步驟（安裝施力和調(diào)整力）中計(jì)算斜拉索施力。使用第3節(jié)中所描述的優(yōu)化方法，來(lái)確定安裝施力，以確保在第9階段中橋面要求的幾何結(jié)構(gòu)（在橫向跨距中增加支座并在主跨距中實(shí)現(xiàn)橋梁合攏），并且施工階段期間該結(jié)構(gòu)的應(yīng)力要保持在容許限值內(nèi)。要做到這一點(diǎn)，設(shè)計(jì)目標(biāo)gj包括在第9階段的位移及第9階段和施工階段期間的應(yīng)力。在第10階段，對(duì)索力進(jìn)行了調(diào)整，以獲得在自重和和額外恒載為2.5kN/m2（橋面、人行道、安全防護(hù)欄和護(hù)欄）的條件下，橋梁施工結(jié)束時(shí)要求的橋面縱剖面。在第10階段后，對(duì)永久狀態(tài)且沒(méi)有時(shí)間效應(yīng)影響下的橋梁進(jìn)行了三種荷載情況的分析，要考慮最大應(yīng)力的設(shè)計(jì)目標(biāo)。荷載情況包括施加于整個(gè)橋面或僅施加于中跨或邊跨上的恒載加上4 kN/m2（道路交通）的活荷載。
3.2設(shè)計(jì)變量
在所分析例子的數(shù)值計(jì)算法中，橋的整體幾何形狀和橋面寬度（19.0m）為預(yù)先指定的恒定設(shè)計(jì)參數(shù)。為確定定徑設(shè)計(jì)變量，考慮了橋塔的兩個(gè)區(qū)域和橋面的三個(gè)區(qū)域。橋面橫截面的幾何形狀設(shè)計(jì)變量分析有三個(gè)選項(xiàng)：梁板、單室箱和三室箱截面。橋塔的設(shè)計(jì)變量，考慮了矩形空心型材。
3.3優(yōu)化設(shè)計(jì)結(jié)果
實(shí)施例1：梁板橫截面，例2：?jiǎn)问蚁浣孛?，?shí)施例3：三室箱截面。索力從橋塔到中跨增加。對(duì)稱分布的斜拉索力相似。最大力發(fā)生在平衡與邊跨相比最大的中央跨度負(fù)荷支索上。支索力基本控制橋塔頂部的水平位移，從而控制橋塔彎曲變形和應(yīng)力。這是斜拉橋充分結(jié)構(gòu)行為的必要特征。實(shí)施例1、2、3的優(yōu)化結(jié)果是成本分別降低17.8%、15.7%和19.1%。這是由于在新設(shè)計(jì)中，減少了橋面、斜拉索和橋塔的橫截面。
永久荷載作用下的全橋橋塔頂部的水平位移約為1cm。結(jié)果表明，在結(jié)構(gòu)收尾橋面的正應(yīng)力仍在使用條件允許的范圍之內(nèi)，并且橋面幾乎完全受壓?？偨Y(jié)了總成本和橋面成本、橋塔和拉索成本。可以說(shuō)，實(shí)施例2的成本最小。橋面占橋成本最大，為總成本的68.8%-71.1%。拉索約占橋梁成本的12%，而橋塔占其余成本。
4.結(jié)論
綜上所述，混凝土斜拉橋的設(shè)計(jì)作為一個(gè)多目標(biāo)優(yōu)化問(wèn)題提出，涉及到拉索預(yù)應(yīng)力和設(shè)計(jì)變量即拉索、橋塔和橋面的橫截面尺寸。優(yōu)化程序?yàn)檎麄€(gè)結(jié)構(gòu)的足夠位移值和應(yīng)力值提供了結(jié)構(gòu)效率解決方案，證明離散直接法靈敏度分析（每次迭代只需要進(jìn)行一次結(jié)構(gòu)分析），在預(yù)測(cè)橋梁結(jié)構(gòu)性能上是有效的。
鑒于相當(dāng)數(shù)量的目標(biāo)和設(shè)計(jì)變量，解決方案只由少量迭代后獲得，因此本文中使用的優(yōu)化算法是有效的。
