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