基于代理優(yōu)化算法的岸橋軌道夾優(yōu)化設(shè)計

余 鎮(zhèn),樊志華,李志華

(杭州電子科技大學(xué)機械工程學(xué)院,浙江 杭州 310018)

0 引 言

隨著經(jīng)濟全球化逐步深入，港口貿(mào)易發(fā)展迅猛。岸橋作為裝卸集裝箱的主要設(shè)備，具有結(jié)構(gòu)高、體型大等特點，也易受臺風(fēng)等惡劣天氣的影響[1]。軌道夾是固定在岸橋底部的機械組件，為了防止岸橋滑動，軌道夾鉗口直接與軌道接觸，是軌道夾中受應(yīng)力最大的零件。代理模型方法是解決復(fù)雜產(chǎn)品和系統(tǒng)優(yōu)化設(shè)計問題的最有效途徑之一，可以提高機械產(chǎn)品質(zhì)量，節(jié)約成本。常見代理模型包括多項式回歸模型、徑向基回歸模型、Kriging模型等。其中，Kriging模型能在給出預(yù)測值的同時給出預(yù)測偏差，成為最具代表性的代理模型[2]。張文豐等[3]采用Kriging代理模型方法修正了岸橋整體結(jié)構(gòu)有限元模型，提高了岸橋的性能校核功能。LI等[4]提出一種關(guān)于岸橋結(jié)構(gòu)模型驗證和可靠性評估的綜合框架，為岸橋可靠性評估提供科學(xué)依據(jù)。劉立成等[5]提出一種合理的湍流模型，提高了岸橋風(fēng)載荷數(shù)值模擬的準確性，更好地模擬了優(yōu)化岸橋結(jié)構(gòu)在惡劣變化風(fēng)載荷下的穩(wěn)定性。設(shè)計合理高效的加點策略是代理優(yōu)化算法的關(guān)鍵。JONES等[6]提出一種基于期望改進(Expected Improvement, EI)準則的有效全局優(yōu)化(Efficient Global Optimization, EGO)算法，是最常用的加點準則。DONG等[7]使用多起點方法找到了代理模型的局部最優(yōu)值，并在這些局部最優(yōu)位置附近搜索全局最優(yōu)解。HAO等[8]提出一種自適應(yīng)多保真預(yù)期改進準則，以較少的分析調(diào)用提供最佳的優(yōu)化結(jié)果。目前，現(xiàn)有的代理優(yōu)化算法大多數(shù)關(guān)注無約束優(yōu)化問題，應(yīng)用范圍受到一定的限制，對可行域邊界的開發(fā)能力較弱[9]。SCHONLAU等[10]使用可行性概率(Probability of Feasibility, PoF)作為EI權(quán)重的約束加點準則，用于增強對可行域的探索。文獻[11]結(jié)合PoF準則和EI準則進行雙目標優(yōu)化，開發(fā)了一類針對黑箱約束問題的并行優(yōu)化算法。

岸橋軌道夾結(jié)構(gòu)優(yōu)化研究中，傳統(tǒng)的結(jié)構(gòu)優(yōu)化方法一般采用經(jīng)驗公式進行計算，在定義復(fù)雜形狀設(shè)計變量和防止網(wǎng)狀失真方面效率較低，為此，本文基于Kriging代理模型，提出一種自適應(yīng)約束代理優(yōu)化算法，通過一個動態(tài)更新Kriging方法對軌道夾結(jié)構(gòu)進行優(yōu)化。

1 Kriging模型與代理優(yōu)化算法

1.1 Kriging代理模型

Kriging代理模型為：

y(x)=μ(x)+Z(x)

(1)

式中，μ(x)為設(shè)計空間的全局模擬，反應(yīng)總體變化趨勢；Z(x)為均值為0、方差為σ2、協(xié)方差非零的隨機過程，表示全局模擬基礎(chǔ)上的局部偏差。

協(xié)方差矩陣和高斯函數(shù)如下：

(2)

Kriging代理模型未知點處的預(yù)測響應(yīng)值為預(yù)測方差可以通過最小二乘法得到，有關(guān)Kriging代理模型的具體描述可參考文獻[2,12]。

1.2 加點準則

1.2.1EI準則

Kriging代理模型不僅可以得到未知點處的預(yù)測響應(yīng)值，還可以得到未知點處的預(yù)測標準差，EI加點準則結(jié)合了Kriging代理模型的這2個特點，表達式為：

(3)

1.2.2 并行加點準則

PEI函數(shù)表達式如下：

(4)

P函數(shù)表達式為：

(5)

式中，n表示加點的個數(shù)，本文取n=2；Nd表示設(shè)計變量的數(shù)目；Tc表示當前迭代步驟的相對時間，是先前迭代和當前迭代所添加樣本數(shù)目與總預(yù)期樣本數(shù)目之比；e-2.3026Tc表示優(yōu)化過程中比例因子的衰減。

2 約束代理優(yōu)化算法

本文提出的自適應(yīng)代理優(yōu)化算法結(jié)合改進的約束并行期望改進準則(Constrained Parallel Expected Improvement, CPEI)、重要邊界采樣準則(Importance Boundary Sampling, IBS)和約束函數(shù)Kriging模型更新策略。

2.1 約束處理

2.1.1PoF準則及CPEI準則

PoF準則計算樣本點落在可行區(qū)域的概率，表達式如下：

(6)

由式(6)可知，PoF(·)是一個取值范圍在0～1的聯(lián)合概率密度。

將PoF準則作為EI準則權(quán)重的約束應(yīng)對方法是現(xiàn)有文獻中較常見的做法。同理，結(jié)合PoF準則，將式(4)中的ymin用最優(yōu)可行解ymin_fea代替，根據(jù)PEI加點準則可以得到約束并行期望改進準則CPEI為：

CPEI(x)=PoF(x)PEI(x)

(7)

2.1.2IBS準則

IBS準則旨在解決黑箱約束函數(shù)問題，在距離當前設(shè)計點較小的約束區(qū)域內(nèi)選擇采樣點。IBS值越大，當前采樣點與約束函數(shù)之間的距離越小。在每個循環(huán)中，IBS準則都會選擇該區(qū)域中IBS值最大的采樣點。IBS函數(shù)表達式如下：

(8)

2.2 冗余點的刪除

對于樣本點的重復(fù)問題，可以通過高斯相關(guān)函數(shù)式(2)來判斷。建立Kriging模型后，可獲得2個樣本點的相關(guān)性參數(shù)R，R∈[0,1]。R越接近0，2個點距離越大；R越接近1，2個點距離越?。划擱=1時，2個點重合。

加入的樣本點過度集中，容易造成樣本點冗余，計算效率降低，為此，本文引入Euclidean距離E來刪除相同樣本點。

(9)

式中，xA表示通過加點準則選取的新樣本點，xN表示加點前的樣本點。若E<0.1，說明新樣本點過度集中，刪除加入的新樣本點，否則可以添加新樣本點。

2.3 約束函數(shù)Kriging模型更新判斷方法

約束函數(shù)Kriging模型插值的不確定性可能會得到不可行的最優(yōu)解，為此，本文引入一種客觀的約束函數(shù)Kriging模型更新策略[13]。約束函數(shù)Kriging模型更新判斷條件如下：

(10)

(11)

2.4 自適應(yīng)約束代理優(yōu)化算法流程

本文提出的自適應(yīng)代理優(yōu)化算法流程包括3個階段。(1)第1階段確定設(shè)計空間。使用拉丁超立方抽樣方法得到樣本庫L，如果初始樣本庫L中沒有滿足約束條件的樣本點，通過PoF準則添加1個可行樣本點。(2)第2階段通過IBS準則尋找離約束區(qū)域邊界較近的樣本點。受PoF準則的影響，CPEI準則更加傾向于在約束區(qū)域的中心區(qū)域取點，對約束邊界區(qū)域的搜索能力較弱，故通過引入IBS準則來增強對約束邊界的探索。同時，為了防止CPEI準則和IBS準則加入的樣本點過于集中，通過Euclidean距離來判斷是否刪除新加入的樣本點。兩種加點準則均使用遺傳算法尋優(yōu)。(3)第3階段在新增樣本點附近，采用信任域法[12]進行取點，調(diào)用Kriging代理模型獲得新增樣本點目標函數(shù)、約束函數(shù)的預(yù)測響應(yīng)值和預(yù)測誤差，得到各信任域內(nèi)滿足約束函數(shù)Kriging模型更新條件的目標函數(shù)最小樣本點，并通過仿真得到所需樣本點的目標函數(shù)、約束函數(shù)真實值，加入樣本庫L，更新Kriging模型，直到滿足終止條件，輸出最優(yōu)結(jié)果。

圖1 一維函數(shù)Kriging曲線

3 算例分析

3.1 PEI加點準則的驗證

為了驗證PEI加點準則具有取得EI函數(shù)多個峰值點的特性，用一維函數(shù)進行驗證。一維函數(shù)如下：

f(x)=(6x-2)2sin(12x-4),x∈[0,1]

(12)

初始采樣點x分別取0.1，0.2，0.3，0.4，0.9，得到一維函數(shù)Kriging模型曲線如圖1所示，P=1和P=2時的PEI函數(shù)分布情況如圖2所示。

圖2 PEI函數(shù)值分布曲線

從圖2(a)可以看出，當P=1時，PEI函數(shù)與EI函數(shù)表達式相同，在x=0.14處得到最大峰值；從圖2(b)可以看出，當P=2時，在x=0.63處得到最大峰值，并且和EI函數(shù)的第2個峰值點十分接近，接近圖1中真實函數(shù)全局最優(yōu)點。這是因為當P=2時，所有PEI函數(shù)值都被放大，而第2個峰值區(qū)域經(jīng)放大后超過第1個峰值區(qū)域。綜上分析可知，PEI準則具有尋找EI函數(shù)多個峰值近似區(qū)域的能力。

3.2 基于約束代理優(yōu)化算法的函數(shù)優(yōu)化

約束條件下的數(shù)學(xué)算例函數(shù)表達式為：

(13)

式中，x1,x2∈[0,1]。該函數(shù)的真實全局最優(yōu)解為x*=[0.195 4,0.404 4]，對應(yīng)真實最小值為0.599 8。函數(shù)在較大連通可行域內(nèi)，存在一個全局最優(yōu)點和多個局部最優(yōu)解[3]，容易陷入局部最優(yōu)區(qū)域，尋優(yōu)困難。因此，這類函數(shù)需要更多的函數(shù)估值，算法的不同，對這類函數(shù)的尋優(yōu)影響較大。

圖3 3種優(yōu)化算法的收斂曲線

通過拉丁超立方抽樣得到10個初始樣本點，分別使用本文算法、CEI單點加點準則和EI-PoF兩目標約束并行加點準則優(yōu)化算例(13)，一直加點至總樣本數(shù)為60。每種算法分別重復(fù)運行20次，得到平均收斂曲線如圖3所示，橫坐標為樣本點數(shù)，縱坐標為絕對誤差A(yù)，A=|fmin(x)-freal(x)|，fmin(x)為函數(shù)最小值,freal(x)為真實最小值。

從圖3中可以看出，本文算法得到的曲線收斂最快，到達收斂近似區(qū)域所用點數(shù)較小，迭代次數(shù)最少，且優(yōu)化結(jié)果精度最高。

3種優(yōu)化算法的函數(shù)優(yōu)化結(jié)果如表1所示。

表1 3種優(yōu)化算法的函數(shù)優(yōu)化結(jié)果

從表1可以看出，和CEI及EI-PoF算法相比，本文算法的迭代次數(shù)分別減少了68.40%和36.80%；達到收斂所需樣本數(shù)相差不大；絕對誤差分別減少了48.17%和46.03%；計算時間最少，和CEI算法的差別比較明顯。由此可見，本文算法的精度和效率最高。因為本文算法并行加點，減少了迭代次數(shù)，效率得到提高；加入CPEI準則和IBS準則探索未知點，精度也得到提高，驗證了本文算法的可行性和有效性。

4 工程實例

4.1 港口集裝箱起重機軌道夾模型

本文主要針對嚴重風(fēng)負載條件下的鉗口型軌道夾進行結(jié)構(gòu)優(yōu)化。軌道夾與軌道直接接觸，使軌道夾與軌道夾緊或分離。文獻[14]確定了對軌道夾鉗口有重要影響的關(guān)鍵性結(jié)構(gòu)參數(shù)，本文通過Solidworks軟件，得到如圖4、圖5所示的三維模型，將鉗口三維模型導(dǎo)入Ansys軟件進行仿真分析，并通過靜力學(xué)分析處理后，得到如圖6所示的鉗口有限元分析圖。

圖5中，變量的初始值設(shè)置為：結(jié)構(gòu)厚度t1=30.0 mm，t2=30.0 mm，t3=30.0 mm，孔中心和曲率中心之間長度l1=30.0 mm，鉗口的重量w=43.5 kg。圖6中，初始設(shè)計區(qū)域產(chǎn)生最大應(yīng)力σmax=545 MPa。鉗口材料為42CrMo，屈服極限不低于823 MPa，根據(jù)歐洲起重機設(shè)計規(guī)范安全系數(shù)取為1.5，許用應(yīng)力σa=549 MPa。

圖4 楔形軌道夾三維模型

圖5 軌道夾鉗口三維模型

圖6 鉗口應(yīng)力輪廓圖

4.2 軌道夾的鉗口結(jié)構(gòu)設(shè)計優(yōu)化

4.2.1 優(yōu)化模型

選取t1，t2，t3，l1作為設(shè)計空間，以鉗口受到的最大應(yīng)力σmax為約束函數(shù)，以鉗口質(zhì)量w為目標函數(shù)，得到鉗口結(jié)構(gòu)優(yōu)化數(shù)學(xué)模型：

minw(t1,t2,t3,l1)

(14)

4.2.2 鉗口結(jié)構(gòu)優(yōu)化結(jié)果

通過拉丁超立方抽樣得到21組樣本點，采用Ansys軟件進行仿真分析，部分結(jié)果如表2所示。

表2 21組軌道夾鉗口結(jié)構(gòu)樣本數(shù)據(jù)

利用21組樣本建立軌道夾鉗口質(zhì)量和應(yīng)力的初始Kriging代理模型，分別采用本文算法、CEI算法和EI-PoF算法對Kriging模型不斷優(yōu)化加點，并及時更新Kriging模型，優(yōu)化加點后獲得最優(yōu)解如表3所示。

表3 優(yōu)化前后軌道夾鉗口參數(shù)對比

從表3可以看出，采用本文提出的自適應(yīng)代理優(yōu)化算法對岸橋軌道夾鉗口進行優(yōu)化后，質(zhì)量減少了13.1%，滿足強度約束條件；對比CEI算法和EI-PoF算法，本文算法的優(yōu)化次數(shù)最少，質(zhì)量最輕，故優(yōu)化效率更高，精度更好。

5 結(jié)束語

針對岸橋軌道夾優(yōu)化過程中存在的問題，本文提出一種自適應(yīng)約束代理優(yōu)化算法?；贙riging模型，結(jié)合IBS準則和CPEI準則，增加了對可行域邊界的刻畫，使得產(chǎn)生的樣本點精確收斂于約束優(yōu)化問題的真實最優(yōu)解，有效減少了約束函數(shù)仿真次數(shù)，提高仿真效率；并在滿足強度約束條件下，對岸橋軌道夾結(jié)構(gòu)進行輕量化設(shè)計，解決了優(yōu)化效率低的問題。但是，本文算法對復(fù)雜約束函數(shù)下最優(yōu)解的搜索能力有所欠缺，后續(xù)將針對復(fù)雜約束下的函數(shù)最優(yōu)解問題展開進一步研究。