輕型堤塊系統(tǒng)概念重建與優(yōu)化

侯武杰，刁景華，同建輝

(1.海軍勤務(wù)學(xué)院，天津 300450；2.陽光大地(天津)環(huán)保工程有限公司，天津 300393)

所謂堤道系統(tǒng)，就是通過快速部署到軍事形勢較為嚴(yán)峻或受災(zāi)地區(qū)，將人員和物資運送上岸的一種橋接系統(tǒng)[1]。

1 堤道系統(tǒng)發(fā)展概況

第二次世界大戰(zhàn)期間，諾曼底登陸時開設(shè)的“桑樹”人工港首次使用了堤道系統(tǒng)。隨后分別研制了能夠適應(yīng)2級海況的海軍駁船系統(tǒng)[2]、模塊化堤道系統(tǒng)[3]和聯(lián)合模塊化駁船系統(tǒng)[4-5]以及能夠適應(yīng)3級海況的高架堤道系統(tǒng)[4-5]。其中海軍駁船系統(tǒng)由長寬分別為27.43 m和6.40 m的鋼模塊組成，但該設(shè)計并不符合國際標(biāo)準(zhǔn)，需要特殊的起重設(shè)備進(jìn)行安放；模塊化堤道系統(tǒng)的設(shè)計符合國際標(biāo)準(zhǔn)，可以采用浮動堤道、滾動或滾降安放設(shè)施、駁運門橋、絞灘拖船等多種不同的方式靈活配置；聯(lián)合模塊化駁船系統(tǒng)由長寬高分別為12.19 m、2.44 m和2.44 m的鋼模塊連接而成，連接和配置方式都很靈活，但存在模塊間出現(xiàn)裂紋、側(cè)向連接系統(tǒng)不穩(wěn)定等風(fēng)險和不足；高架堤道系統(tǒng)由模塊化部件焊接而成，通過轉(zhuǎn)盤實現(xiàn)了物資的雙向傳輸，是目前最有效的海岸間物資運輸?shù)牡痰老到y(tǒng)。

上述堤道系統(tǒng)在特定的歷史環(huán)境下都發(fā)揮了一定的作用，隨著時代的發(fā)展，堤道系統(tǒng)需要戰(zhàn)略性海運資源(大型運輸船等)運輸和部署的矛盾也逐漸凸顯，為此設(shè)計了輕型堤塊系統(tǒng)。該系統(tǒng)由長3.05 m、寬6.10 m的12個模塊單元組成，上部為鋁制結(jié)構(gòu)，下部為氣動浮筒，模塊之間采用端到端進(jìn)行連接[6]，便捷高效。

該系統(tǒng)的設(shè)計主要考慮了四個方面的因素：一是確保吃水較淺的船只能夠在惡劣環(huán)境(如泥灘或濕地)或港口損壞(淺水或軟土條件)的情況下迅速展開堤道部署；二是設(shè)計符合國際標(biāo)準(zhǔn)配置，可綜合運用陸運、海運、空運等多種方式實施戰(zhàn)略投送；三是便于組裝；四是模塊單元重量輕，體積小。其最大優(yōu)勢在于可以使用聯(lián)合高速船、后勤支援艦或軍用駁船運輸，通過浮筒充氣、模塊連接等方式，使堤道系統(tǒng)快速部署到水深較淺的前沿軍事地域，并且能夠支撐多輛主戰(zhàn)坦克同時通過。

2 輕型堤塊系統(tǒng)概念重建的思路

輕型堤塊系統(tǒng)在現(xiàn)場實驗和評估中顯示出了良好的前景和性能，但連接設(shè)計比較復(fù)雜，連接部件重量占到每個模塊自身重量的一半，需要重新定義現(xiàn)有的設(shè)計，選擇現(xiàn)有鉸鏈的旋轉(zhuǎn)柔度，充分利用下部浮筒的浮力實現(xiàn)旋轉(zhuǎn)。

圖1為輕型堤塊系統(tǒng)的三種不同配置形式，曲線表示水位。若系統(tǒng)是剛性的，即模塊之間進(jìn)行完全力矩連接，那么該結(jié)構(gòu)就不能充分利用波浪間的浮力(如圖1-a所示)。若能采用柔性連接(如圖1-b所示)，就可以更好地利用波浪間的浮力(如圖1-c所示)。

圖1 堤道系統(tǒng)配置Fig.1 Configurations for causeway system

圖2-a為輕型堤塊系統(tǒng)內(nèi)部的具體連接方式，這種連接方式既繁重又復(fù)雜。概念重建的目的就是通過將旋轉(zhuǎn)柔度傳遞到甲板自身來消除剛性連接，在給定載荷條件下，可通過將甲板剛度降低到目標(biāo)值來實現(xiàn)特定的旋轉(zhuǎn)。設(shè)計目標(biāo)是在1 140 kN·m荷載條件下，使甲板在3.05 m長的模塊上旋轉(zhuǎn)0.08 rad(堤道支撐一個坦克和斜坡負(fù)載時的設(shè)計彎矩)。圖2-b則是一種更輕型、更簡單的模塊連接方式，且模塊可以分開進(jìn)行運輸。此外，模塊采用重量輕、耐腐蝕的玻璃纖維增強聚合材料設(shè)計，并將其橫截面改為箱梁結(jié)構(gòu)。

2-a 輕型堤塊系統(tǒng)原設(shè)計的連接方式 2-b 用于概念重構(gòu)系統(tǒng)的力矩連接圖2 鉸鏈配置Fig.2 Hinge configurations

3 輕型堤塊系統(tǒng)的設(shè)計優(yōu)化

堤道系統(tǒng)的自身重量和體積對于可部署結(jié)構(gòu)設(shè)計來講非常重要。對于概念重建的輕型堤塊系統(tǒng)而言，在符合美國鋼結(jié)構(gòu)規(guī)范中關(guān)于拉擠玻璃鋼結(jié)構(gòu)條件和幾何要求前提下，運用多目標(biāo)結(jié)構(gòu)優(yōu)化算法，以最小目標(biāo)剛度值為約束條件，求得系統(tǒng)的最小自重和剛度。目前，多目標(biāo)結(jié)構(gòu)優(yōu)化算法被廣泛應(yīng)用于浮動甲板設(shè)計、永久橋面鋪裝系統(tǒng)設(shè)計以及橋梁修復(fù)工程中，其中遺傳算法是應(yīng)用最廣泛的一種優(yōu)化算法，其不僅能避免局部最小值，還可以處理離散和不可微分優(yōu)化函數(shù)。

3.1 問題分析

概念重建就是探索一種具有目標(biāo)剛度的輕量化設(shè)計方法。即在符合鋼結(jié)構(gòu)規(guī)范和封裝標(biāo)準(zhǔn)的基礎(chǔ)上，以慣性矩不低于目標(biāo)值為約束條件，對最小自重和最小慣性矩進(jìn)行多目標(biāo)結(jié)構(gòu)優(yōu)化。

其中，設(shè)計變量與箱梁截面有關(guān)，包括整個截面深度H、上下翼緣厚度tft和tfb、下翼緣寬度wfb、外腹板厚度twe(圖3)。

圖3 一般的橫截面Fig.3 General cross section

這些變量從9.5 mm開始以3.2 mm的增量不斷增加到6 096 mm。圖3中wft為上翼緣寬度、hw為腹板高度、bwe為外翼緣水平寬度，按照設(shè)計要求，上翼緣寬度wft=6 096 mm為固定值，內(nèi)部腹板的數(shù)量n和厚度twi在算法中確定。最后兩個維度與設(shè)計變量的幾何關(guān)系如下

hw=H-tft-tfb

(1)

(2)

這是簡化后的方程。該優(yōu)化問題的正式定義如下

使得

(3)

式中：W為上部結(jié)構(gòu)的自重，即每個元素的面積乘以密度p的總和。轉(zhuǎn)動慣量目標(biāo)函數(shù)Ix由平行軸定理得到，其中d為單個構(gòu)件質(zhì)心到整個上部構(gòu)件質(zhì)心的距離。最后，算法生成的截面必須符合u個約束條件，其中c1代表最小轉(zhuǎn)動慣量值，c2～c6為鋼結(jié)構(gòu)規(guī)范手冊的約束條件，幾何約束條件c7、c8為封裝標(biāo)準(zhǔn)。表1列出了軍事橋接系統(tǒng)的負(fù)載組合作用于系統(tǒng)的設(shè)計彎矩、剪力和扭矩的具體數(shù)值，并將其分別乘以荷載系數(shù)1.33以滿足鋼結(jié)構(gòu)規(guī)范設(shè)計。

表1 設(shè)計彎矩、剪力和扭轉(zhuǎn)Tab.1 Design moment, shear and torsion

第一個約束條件c1通過設(shè)置一個較低的限制，以確保上部結(jié)構(gòu)的慣性矩不超出預(yù)期的靈活性。即

(4)

式中：p為所需的曲率半徑；MS為工作荷載作用下的彎矩(如表1所示)；EL為縱向彈性模量。在運用彎矩與曲率關(guān)系時，假定截面為平面，材料為有限彈性材料。均勻彎矩下，在3.05 m長的模塊上期望旋轉(zhuǎn)0.08 rad，目標(biāo)慣性矩為0.003 503 m4。

鋼結(jié)構(gòu)規(guī)范手冊中定義的結(jié)構(gòu)性限制標(biāo)準(zhǔn)，與材料切變(材料在剪切破壞、腹板剪切失穩(wěn)和腹板橫向穩(wěn)定性等因素作用下而引起的構(gòu)件強度)、彎曲(材料斷裂和局部不穩(wěn)定性導(dǎo)致的構(gòu)件強度)、扭矩(強度控制下的扭轉(zhuǎn)能力)和集中負(fù)荷需求(腹板拉伸斷裂、腹板破壞、腹板壓縮屈曲和翼緣彎曲破壞而引起的構(gòu)件強度)等因素有關(guān)。在所有的計算中，假定玻璃纖維增強聚合材料構(gòu)件具有相同的性能，且不發(fā)生分層或分離。

第二個約束條件是基于材料剪切破壞與腹板剪切屈曲之間的控制行為。剪切能力Vn必須超過剪力需求Vu，關(guān)系如下所示

c2=Vu-λφVn≤0

(5)

式中：λ=0.8為時間效應(yīng)因子；φ為剪切阻力系數(shù)，系統(tǒng)的剪力需求為表1給出的值乘以荷載系數(shù)1.33。對于剪切過程中的材料斷裂，φ=0.65，Vn=FLTAs，其中FLT為平面抗剪強度特性、As為剪切面積(As=Hbwe)。此處的剪切面積指的是單腹板的面積，實際上剪力會通過幾個或全部腹板傳遞到截面上，由于每個腹板所承載的剪力需求只能通過詳細(xì)的有限元分析得出，為了計算方便，假設(shè)任何單個腹板必須能夠承載全部的剪力需求，最終設(shè)計時，通過詳細(xì)分析每塊腹板的實際需求來減少腹板的寬度。對于腹板剪切屈曲，φ=0.80，Vn=FcrAs，其中Fcr為臨界剪切屈曲應(yīng)力，為設(shè)計變量twe和H的函數(shù)。

第三、四個約束條件用來檢查該部分是否滿足彎曲要求。約束條件c3是確保材料在彎曲過程中不會斷裂。彎曲能力Mn必須超過需求Mu，關(guān)系如下所示

c3=Mu-λφMn≤0

(6)

式中：阻力系數(shù)φ=0.65，彎曲能力的計算公式為Mn=FLIx/y，其中FL為縱向抗彎強度，y為截面中性軸到構(gòu)件極端纖維的距離。Mu為因式矩(1.33乘以表1設(shè)計值)。

第二個彎曲約束是用來檢查外部腹板是否足夠厚，防止局部不穩(wěn)定。即

c4=twer-twe≤0

(7)

twer為所需的腹板厚度，通過如下公式求得

(8)

式中：φ=0.65；EL,w和ET,w分別表示腹板的縱向彈性模量和橫向彈性模量；vLT為縱向泊松比；GLT為平面剪切模量。

第五個約束條件c5與構(gòu)件的扭轉(zhuǎn)能力有關(guān)，定義為

c5=Tu-λφTn≤0

(9)

式中：Tu為所需扭轉(zhuǎn)量(1.33乘以表1設(shè)計值)；φ=0.70；Tn為抗扭承載能力。封閉截面的抗扭承載力為Tn=2tFLTAo，其中t為任意外部單元的最小厚度，Ao為外部單元中心線圍成的面積。對于集中負(fù)荷，負(fù)荷需求Ru必須超過容量Rn，如下所示

c6=Ru-λφRn≤0

(10)

系統(tǒng)的負(fù)荷需求主要是坦克的集中載荷以及船與堤道相連的坡道重量，坦克的集中載荷可以看成是其重量除以車輪的數(shù)量(Ru=68 kN，包括負(fù)載系數(shù))。該約束條件主要考慮了腹板拉伸斷裂、腹板破壞和腹板壓縮屈曲等三種極限狀態(tài)的情況，分別對內(nèi)外腹板進(jìn)行了詳細(xì)分析，并使得φRn最小。對于腹板拉伸斷裂的情況，Rn=ltenFT,wtw，其中l(wèi)ten為腹板高度，F(xiàn)T,w為腹板的橫向抗彎強度，同樣假設(shè)需要一個腹板來承載全部負(fù)載，因此tw的值為內(nèi)外腹板中最薄腹板的厚度(twe或twi)，在這種極限狀態(tài)下φ=0.65。對于腹板破壞嚴(yán)重的情況，φ=0.70，Rn=0.7hwtwFsh,int{1+[(2k+6tplate+bplate)/dw]}，式中Fsh,int為層間剪切強度，k為構(gòu)件頂部到圓角底部的距離，由于截面上沒有圓角，所以取上翼緣的厚度tft代替，tplate和bplate分別為承重板的厚度和長度，由于不存在承重板，取上翼緣的厚度和寬度代替，最大值為102 mm，dw為腹板深度，在本案例中就是腹板高度hw。對于腹板壓縮屈曲的情況，φ=0.80，Rn=fcrAeff，其中Aeff為有效面積(Aeff=lefftw，其中l(wèi)eff為腹板深度dw或上下翼緣之間垂直距離的較小者)，則fcr為

(11)

第七個、第八個約束條件是與功能需求、包裝需求和物理約束相關(guān)的幾何約束。上、下翼緣之間的夾角必須大于45%，可正式表示為

c7=wft-wfb-2H≤0

(12)

對于包裝需求來說，橫截面的深度必須足夠淺，使得4個模塊可以安裝在一個標(biāo)準(zhǔn)容器中(最大深度為400 mm)。

c8=H-400 mm≤0

(13)

第九個、第十個幾何約束是確保模塊的橫截面如實存在，具體來說，上、下翼緣的厚度不能大于整個構(gòu)件的高度。

c9=tft+tfb-H≤0

(14)

并且外部腹板的厚度不能大于橫截面的寬度。

c10=2twe-wft≤0

(15)

內(nèi)部腹板的數(shù)量和厚度沒有設(shè)置為變量，而是通過計算以滿足各種極限狀態(tài)。受壓翼緣構(gòu)件的強度是由整個截面彎曲所引起的局部不穩(wěn)定性和翼緣抗壓彎曲破壞造成的，而內(nèi)腹板的數(shù)量正是按照保持受壓翼緣構(gòu)件穩(wěn)定的標(biāo)準(zhǔn)設(shè)計的，由于截面的上翼緣和下翼緣都能受壓，因此有必要對它們分別進(jìn)行分析，并提出基于集中荷載要求nwi的計算方法

(16)

式中：FT，f為翼緣橫向抗彎強度，φ=0.65，nwi為這兩個計算之間的控制值。假設(shè)內(nèi)部腹板在整個截面上均勻間隔，并且集中力作用在最關(guān)鍵的位置，即兩個腹板之間的中間位置，腹板厚度是由材料在剪切破壞、腹板剪切失穩(wěn)和彎曲局部失穩(wěn)之間的控制計算確定的，與外部腹板計算的情況一樣，仍然假設(shè)任何單個腹板必須能夠承載全部的剪力需求，材料剪切破壞與腹板剪切屈曲仍然由c2中描述的相同方程導(dǎo)出，重新排列后得到twi(代替了c2方程中的bwe和twe)，根據(jù)彎曲不穩(wěn)定性檢驗，所要求的腹板厚度與c4中描述的相同。

如前所述，對A709 HPS 70鋼(屈曲強度為483 MPa)和A709 HPS 100鋼(屈曲強度為689 MPa)兩種高強度鋼進(jìn)行了類似的處理。設(shè)計變量和目標(biāo)函數(shù)保持不變，但是結(jié)構(gòu)約束c2～c6由美國公路橋梁設(shè)計規(guī)范中的約束代替，用于鋼箱梁的設(shè)計。目標(biāo)轉(zhuǎn)動慣量c1也進(jìn)行了改變，以獲得相同的剛度和不同的彈性模量0.000 217 m4。

3.2 模擬退火算法

模擬退火算法是一種基于晶體形成類比的迭代改進(jìn)算法，在物理降溫過程中，隨著熔化體的緩慢冷卻，系統(tǒng)的能量逐漸降低，在此過程中存在可能發(fā)生更高能量配置的概率P，這最終會導(dǎo)致更低的能量配置。概率公式P=e(-ΔE)/T，其中ΔE為能量構(gòu)型差，T為物體溫度，當(dāng)溫度降低時，發(fā)生高能量狀態(tài)的可能性就會降低。這一過程可以擴展到結(jié)構(gòu)優(yōu)化，其中概率與接受高值函數(shù)的概率有關(guān)，能量與現(xiàn)有函數(shù)有關(guān)，T為用戶可以控制的變量，從而使算法能夠擺脫局部最小值。

對于單目標(biāo)模擬退火算法(如最小權(quán)值優(yōu)化)，該算法首先從用戶指定的離散值數(shù)據(jù)庫中隨機生成一組可行的設(shè)計變量，選擇初始解。這個集合成為初始的最佳解。然后隨機擾動一個或多個變量，生成一個新的解，如果該解符合所有約束條件，并能得到一個較低權(quán)值的解，它就成為新的通用解，當(dāng)然也存在一定概率沒有最佳解，但它仍然可以作為當(dāng)前的解決方案接受，算法在此基礎(chǔ)上繼續(xù)迭代。

這種對單個設(shè)計變量進(jìn)行優(yōu)化的過程可以推廣到多目標(biāo)優(yōu)化，不是使用加權(quán)平均將目標(biāo)函數(shù)組合在一起，而是將目標(biāo)函數(shù)分離以產(chǎn)生帕累托最優(yōu)解集，如果在任何一個目標(biāo)函數(shù)中，該解都不被其他解所掩蓋，那么該解就是帕累托最優(yōu)解。在收斂過程中，帕累托最優(yōu)解集為設(shè)計者提供了一系列可能的解，這些解跨越了通過單目標(biāo)優(yōu)化可以找到的兩個極端，根據(jù)設(shè)計優(yōu)先級，設(shè)計師可以選擇最終的解決方案。

3.3 優(yōu)化設(shè)計

從本質(zhì)上講，啟發(fā)式算法并不能保證每次使用時都收斂于同一個解。用戶定義了幾個參數(shù)，包括v(一次變量變化的最大數(shù)量)、pm(允許值數(shù)據(jù)庫允許的擾動量)、r(降溫因子)、m(冷卻周期的長度)和n(算法收斂的冷卻循環(huán)次數(shù))。算法的質(zhì)量和魯棒性取決于這些參數(shù)的選擇(魯棒性是指同時具有較低的平均值和標(biāo)準(zhǔn)差)。為了確定這些參數(shù)對多目標(biāo)優(yōu)化的魯棒性選擇，對每個目標(biāo)函數(shù)使用16種不同的參數(shù)組合v、pm、r、m和n進(jìn)行單目標(biāo)優(yōu)化，對每種組合進(jìn)行20次數(shù)值模擬。表2給出了考慮的組合以及玻璃纖維增強聚合材料設(shè)計的平均結(jié)果μ、標(biāo)準(zhǔn)差σ和變異系數(shù)cv。對于兩個目標(biāo)函數(shù)，最穩(wěn)健的組合是SA4(v=1；pm=10；r=0.8；m=10 000；n=2；在表2中以粗體突出顯示)。

表2 玻璃纖維增強聚合材料設(shè)計的模擬退火數(shù)值實驗Tab.2 Simulated annealing numerical tests for GFRP design

通過這些參數(shù)的選擇，對玻璃纖維增強聚合材料進(jìn)行多目標(biāo)模擬退火算法設(shè)計，找到了帕累托最優(yōu)結(jié)果集，設(shè)計者可以從中選擇最終的設(shè)計(如圖4中空心圓所示)。菱形和方形分別代表了重量和慣性矩的單目標(biāo)優(yōu)化的最佳結(jié)果，代表了帕累托最優(yōu)集跨越的兩個極值，虛線表示目標(biāo)轉(zhuǎn)動慣量值，轉(zhuǎn)動慣量和重量的最佳組合由一個填充圓所表示的橫截面來表示。該方案的慣性矩接近目標(biāo)值(0.003 504 m4)，重量也可以接受(包括浮體重量在內(nèi)，為1 208 kg/m)，對應(yīng)的截面如圖5-a所示。對這兩種等級的鋼進(jìn)行了類似的處理，得到的優(yōu)化截面如圖5-b和5-c所示。

圖4 多目標(biāo)模擬退火算法數(shù)值模擬得出的帕累托最優(yōu)解集Fig.4 Pareto-optimal (PO) set of solutions from one multiobjective simulated annealing (MOSA) numerical simulation

5-a 優(yōu)化的玻璃纖維增強聚合材料截面

表3對優(yōu)化結(jié)果進(jìn)行了總結(jié)，并與現(xiàn)有的輕型堤塊系統(tǒng)進(jìn)行了比較，所采用的三種優(yōu)化設(shè)計均達(dá)到或接近所需的系統(tǒng)目標(biāo)剛度，因此可以根據(jù)需要設(shè)計鉸鏈的連接方式，可采用圖2-b所示的簡化連接替換。除了簡化連接之外，最重要的任務(wù)之一就是減少系統(tǒng)的總重量，表3給出了單獨的模塊重量W(不含鉸鏈和浮漂)、有連接的模塊重量(含鉸鏈，不含浮漂)以及包含浮漂的系統(tǒng)總重量。但由于頂部翼緣的局部約束，玻璃纖維增強聚合材料的設(shè)計并不能使系統(tǒng)的重量降低，而兩種鋼的設(shè)計都比玻璃纖維增強聚合材料設(shè)計的方案重量輕，其中A709 HPS 100鋼的設(shè)計方案使得堤塊系統(tǒng)重量較現(xiàn)有系統(tǒng)重量輕。此外，系統(tǒng)的上部結(jié)構(gòu)高度H被降低，使得一個標(biāo)準(zhǔn)容器可以裝運更多的模塊，減少了運輸時間和包裝成本。

表3 輕型堤塊系統(tǒng)原型與設(shè)計方案的比較Tab.3 Comparison between original prototype and proposed designs

4 結(jié)論和展望

本文研究提出了一種基于玻璃纖維增強聚合材料和兩種高強度鋼的多目標(biāo)結(jié)構(gòu)優(yōu)化的輕型堤塊系統(tǒng)，所得到的設(shè)計結(jié)果均達(dá)到或接近系統(tǒng)所需的目標(biāo)剛度，以取代繁重復(fù)雜的連接，并建議簡化連接。這比原系統(tǒng)的兼容連接更簡單、更輕便、更容易部署和維護(hù)，也使得系統(tǒng)上部結(jié)構(gòu)設(shè)計制造和部署更為簡單。每一個優(yōu)化設(shè)計降低了截面深度，從而減少了運輸成本和時間。雖然玻璃纖維增強聚合材料和A709 HPS 70鋼的設(shè)計并不能降低整個系統(tǒng)的重量，但是A709 HPS 100鋼的設(shè)計達(dá)到了降低系統(tǒng)自重的目標(biāo)，從而實現(xiàn)了該項目的設(shè)計優(yōu)先級。

在對重新定義的輕型堤塊系統(tǒng)進(jìn)行初步優(yōu)化設(shè)計的基礎(chǔ)上，下一步需要繼續(xù)深入研究詳細(xì)的連接設(shè)計，分析通過減小腹板寬度來減輕系統(tǒng)自重的可能性，并要考慮材料和制造成本，以及橫截面的可構(gòu)造性和部件之間的連接細(xì)節(jié)等。