基于組合代理模型的箱形梁輕量化設(shè)計*

陳 楊 秦義校 楊俊樂

太原科技大學(xué)機械工程學(xué)院 太原 030024

0 引言

由于復(fù)雜工程問題的有限元建模時間和求解時間過長，且在尋優(yōu)過程中需要大量的仿真數(shù)據(jù)，會使復(fù)雜工程問題的尋優(yōu)耗費大量時間。針對此問題，研究人員迫切需要一種能夠高效求解工程問題的代理模型，使復(fù)雜工程問題的尋優(yōu)得到合理的解決方法。

代理模型又稱為元模型或響應(yīng)面模型，是基于少量數(shù)據(jù)而構(gòu)建的一種真實模型的模型，即模型的模型，具有結(jié)構(gòu)簡單、計算效率高等特點[1]。最常用的代理模型包括多項式響應(yīng)面（Respones Surface methodology，RSM）、Kriging 模型、徑向基（Radial Bssis Function，RBF）函數(shù)等，很多學(xué)者針對代理模型的精確問題進行了研究[2]，并對RSM、Kriging 和RBF這3 種代理模型進行了對比[3-5]，在實際工程中驗證了其優(yōu)越性[6]。

針對箱形梁、薄板結(jié)構(gòu)的優(yōu)化主要集中在輕量化設(shè)計和對加強筋的優(yōu)化（如利用有限元軟件對橋式起重機合理、輕量化、節(jié)能化的箱形梁結(jié)構(gòu)等）[7]。Pinca C B等[8]利用Cosmos/M 有限元軟件對冶金起重機的箱形梁截面進行輕量化設(shè)計；Abid M 等[9]針對橋式起重機箱形梁中加強筋的個數(shù)、形狀、位置進行優(yōu)化，并優(yōu)化了箱形梁的截面形狀，以獲得最大的允許彎曲應(yīng)力；Zhang H 等[10]將函數(shù)逼近法與梯度下降法結(jié)合，通過對加強筋布局位置的優(yōu)化，提出了一種非對稱的加強筋布局形式；Putra G L 等[11]提出了一種混合遺傳算法（GA）對船舶結(jié)構(gòu)的加強筋數(shù)量、類型以及加強筋間距和板厚進行優(yōu)化，得到了更加輕量化的船舶結(jié)構(gòu)；Ehsani A 等[12]提出針對薄板結(jié)構(gòu)的輕量化問題，提出了一種格柵結(jié)構(gòu)的可變肋板模型（VRM）找到網(wǎng)格板的最佳布局方式。

本文針對橋式起重機箱形梁結(jié)構(gòu)進行輕量化設(shè)計，首先建立了有限元模型確定箱形梁各部位的優(yōu)化變量[13]，利用最優(yōu)拉丁超立方對設(shè)計變量抽樣，并對RSM、RBF 和Kriging 這3 種代理模型進行擬合和插值[14]。其次，對得到的3 種代理模型進行誤差分析。最后，提出了一種基于代理模型和自適應(yīng)模擬退火算法（Adaptive Simulated Annealing，ASA）以及序列二次規(guī)劃法（Sequence Quadratic Programing，SQP）相結(jié)合的優(yōu)化策略，對3 種模型在滿足強度、剛度和穩(wěn)定性條件下進行輕量化設(shè)計。

1 目標(biāo)函數(shù)及約束條件

橋式起重機的箱形梁通常在高溫、重載等惡劣環(huán)境中工作。由于箱形梁一般為細(xì)長結(jié)構(gòu)，在彎扭、屈曲變形作用下容易產(chǎn)生局部失穩(wěn)和開裂，從而對橋式起重機的正常運行產(chǎn)生影響。因此，對箱形梁的結(jié)構(gòu)進行優(yōu)化布局，進而得到更好的力學(xué)性能是十分必要的[15，16]。

為了得到更加合理的布局形式，本文以箱形梁跨中最大應(yīng)力為驗證條件，即目標(biāo)函數(shù)為最小化跨中最大應(yīng)力。在實際情況下，布局方式與箱形梁跨中最大應(yīng)力是高度非線性方程，代理模型能夠模擬這種黑盒特性及計算復(fù)雜性。為便于箱形梁的各個部位優(yōu)化，本文對各個優(yōu)化部位進行歸一化操作，優(yōu)化結(jié)束后再反歸一化得到對應(yīng)的優(yōu)化結(jié)果。

如圖1 所示，本文的設(shè)計變量共有9 個，各個優(yōu)化部位用 替代整個優(yōu)化過程，其中X＝[x1，x2，x3，…，x9]T。x1和x2為主腹板和副腹板的厚度，x3為隔板的厚度，x4和x5為上翼緣板和下翼緣板的厚度，x6為上下翼緣板的間距，x7為腹板的間距，x8和x9為隔板孔寬和隔板孔長。

圖1 橋式起重機箱形梁變量設(shè)計圖

本文用G1（X）表示設(shè)計變量與跨中最大應(yīng)力，箱形梁的變形同樣采用代理模型進行模擬，用G2（X）表示設(shè)計變量與跨中最大變形的關(guān)系，優(yōu)化前箱形梁滿載時跨中最大應(yīng)力σmax＝200.69 MPa，跨中最大變形Ymax＝20.02 mm，由起重機設(shè)計準(zhǔn)則可得約束條件。

1）強度條件

式中：G1（X）為變量跨中最大應(yīng)力，σs為材料屈服強度，σb為材料抗拉強度。

2）剛度條件

式中：G2（X）為變量跨中最大變形，[Y]為箱形梁的許用變形。

最終得到的優(yōu)化方程為

式中：G1（X）為目標(biāo)函數(shù)，G2（X）為約束條件，Lb[X]為變量的下界，Ub[X]為變量的上界。

為了在優(yōu)化布局方式的同時不使模型質(zhì)量增加，加入質(zhì)量條件約束，設(shè)優(yōu)化前的質(zhì)量為Mt，優(yōu)化后的質(zhì)量為M，則有

2 組合代理模型全局優(yōu)化算法

本文采用自適應(yīng)模擬退火算法（ASA）與序列二次規(guī)劃法（SQP）相結(jié)合的算法（即ASA-SQP 算法），以此得到箱形梁的最佳布局方式。模擬退火算法是模擬金屬熱態(tài)加工過程而發(fā)展的算法，與遺傳算法有很多相似之處。模擬退火算法比遺傳算法更易于實現(xiàn)，因為其在每次搜索空間只檢查1 個設(shè)計點，而遺傳算法對1 組尋優(yōu)點進行尋優(yōu)。模擬退火算法的基本思路是從初始點開始對目標(biāo)函數(shù)進行評估，函數(shù)值下降，新的尋優(yōu)點即可被接受，如此反復(fù)直至找到最優(yōu)點。為了避免陷入局部最優(yōu)，根據(jù)Metropolis 判斷準(zhǔn)則函數(shù)值上升的點亦可被接受。自適應(yīng)模擬退火算法具有比模擬退火算法更強大的全局尋優(yōu)能力和計算效率，為了加速其在找到全局最優(yōu)解最近的尋優(yōu)速度，引入了梯度算法（序列二次規(guī)劃法）有更強的局部尋優(yōu)能力。

序列二次規(guī)劃法（SQP）是將目標(biāo)函數(shù)進行二次泰勒級數(shù)展開，并將約束條件線性化，通過解二次規(guī)劃得到下一個設(shè)計點，然后根據(jù)2 個可供選擇的優(yōu)化函數(shù)執(zhí)行一次線性搜索，其中Hessian 矩陣由BFGS 公式更新。

代理模型的建模首先需要對樣本空間進行抽樣，目前常見的抽樣方法為拉丁超立方設(shè)計(Latin Hypercube Design，LHD)，其原理是在n維設(shè)計空間中，將每一維坐標(biāo)區(qū)間［xkmin，xkmax］，k∈［1，n］均勻地等分成m個區(qū)間，每個小區(qū)間被記為隨機選取m個點，保證一個因子的每個水平只被研究1 次，即構(gòu)成維度為n，樣本數(shù)為m的拉丁超立方設(shè)計，記為m×nLHD。

拉丁超立方設(shè)計比傳統(tǒng)試驗設(shè)計能夠更好地布滿整個設(shè)計空間，但仍可能存在試驗點分布不夠均勻、丟失設(shè)計區(qū)域的情況。最優(yōu)拉丁超立方設(shè)計(Optimal Latin Hypercube Design，OLHD)改善了拉丁超立方設(shè)計的均勻性，使因子和響應(yīng)的擬合更精確真實，圖2 為最優(yōu)拉丁超立方設(shè)計能生成更均勻樣本點的分布情況。

圖2 最優(yōu)拉丁超立方設(shè)計對比示意圖

本文的優(yōu)化流程分為建模和優(yōu)化2 個模塊。首先確定優(yōu)化變量并對其進行參數(shù)化建模；其次利用最優(yōu)拉丁超立方設(shè)計對設(shè)計變量進行優(yōu)化，為了得到更好擬合效果，對優(yōu)化變量進行歸一化操作，并在設(shè)計空間內(nèi)隨機生成200 個樣本點；然后對3 種代理模型分別進行擬合、插值，在設(shè)計空間中隨機生成50 個誤差樣本點，對誤差進行分析；最后，利用ASA-SQP 算法對3 種代理模型進行優(yōu)化，并將優(yōu)化結(jié)果與數(shù)值仿真進行對比驗證，得到最終優(yōu)化結(jié)果。

為了更好地評估3 種不同代理模型的建模精度，采用相同的樣本點和誤差點對3 種代理模型進行誤差評估，工程上常用的代理模型誤差分析準(zhǔn)則為確定系數(shù)R2，即

式中：m為新隨機生成的驗證點個數(shù)，yi為驗證點上的真實值，為代理模型生成的替代值，為平均值。

確定系數(shù)R2是一種直觀反映模型預(yù)測能力的評價標(biāo)準(zhǔn)，也是使用最多的一種誤差判斷準(zhǔn)則，其值越接近1 表示預(yù)測模型與真實值越接近。一般地，當(dāng)R2＞0.9 時，代理模型的擬合效果有效。

在約束條件中，共有3 個輸出變量，分別是質(zhì)量Mass、最大應(yīng)力Maxs 和跨中最大位移Maxu，此時3個輸出變量在3 個代理模型下的確定系數(shù)如圖3 所示。

圖3 3 個輸出變量的確定系數(shù)

由圖3 可知，Kriging 模型對3 個變量的代理均不符合確定系數(shù)的要求，而代理模型RBF、RSM 在Mass和Maxu 上的代理效果趨近于1，而在最大應(yīng)力上也不符合確定系數(shù)的要求。

由有限元軟件中改變輸入變量可知，最大應(yīng)力隨變量的改變其位置也發(fā)生改變，而代理模型的代理效果往往只針對于某一個區(qū)域的最大應(yīng)力，如箱形梁的最大變形總發(fā)生于跨中，故其代理效果往往尚佳。對有限元模型進一步分析，可知箱形梁的最大應(yīng)力基本發(fā)生在跨中下蓋板及2 個端梁處，如圖4 所示。

圖4 優(yōu)化前應(yīng)力云圖

綜上所述，為了得到更好代理效果，將Maxs 變量替換為跨中下蓋板及2 個端梁處的最大應(yīng)力，分別表示為Maxs_1、Maxs_2、Maxs_3，重復(fù)上述代理模型舍棄掉Kriging 模型，得到的確定系數(shù)R2如圖5 所示。

圖5 確定系數(shù)R2

由圖5 可知，通過添加代理變量，RBF 和RSM 在5 個代理變量下均達(dá)到確定系數(shù)許用要求。然而為了得到更好代理效果，從5 個變量確定系數(shù)的絕對大小分別看，提出一種組合代理模型，即Mass 和Maxs_1 采用RSM 代理，而Maxs_2、Maxs_3 和Maxu 變量采用RBF 代理，得到一種組合代理模型RSM-RBF。

3 優(yōu)化算法

利用ASA-SQP 算法對3 種代理模型進行優(yōu)化，并將優(yōu)化結(jié)果與數(shù)值仿真進行對比驗證，得到最終優(yōu)化結(jié)果。對RSM、RBF 和RSM-RBF 這3 種代理模型進行優(yōu)化得到迭代圖像如圖6 所示，3 種代理模型優(yōu)化結(jié)果與有限元的對比如表1 所示。

表1 各代理模型優(yōu)化對比

圖6 收斂迭代圖

由表1 可知，基于RBF 得到的最小質(zhì)量最大，但其在Maxs_2、Maxs_3 以及Maxu 這3 個變量的誤差要顯著小于RBF代理模型，這與確定系數(shù)的誤差判斷一致。RSM 以及RSM-RBF 這2 種代理模型優(yōu)化得到的質(zhì)量相差無幾，但RSM-RBF 代理5 個變量與有限元結(jié)果進行對比所得到的綜合誤差要優(yōu)于組合變量之前的兩種單獨代理模型的誤差。綜上，本文最終采用RSM-RBF 組合代理模型作為本文最終優(yōu)化所用的模型。

4 優(yōu)化結(jié)果

各部分的優(yōu)化結(jié)果如表2 所示，所得到優(yōu)化后的質(zhì)量為55.33 t，相比于優(yōu)化前的65.54 t，箱形梁的質(zhì)量降低了15.59%。

表2 各代理模型優(yōu)化對比

其優(yōu)化后的應(yīng)力云圖7 所示，結(jié)果符合強度、剛度條件。

圖7 優(yōu)化后的應(yīng)力云

5 結(jié)論

本研究針對工程優(yōu)化中常用的3 種代理模型，即多項式響應(yīng)面（RSM）、RBF 神經(jīng)網(wǎng)絡(luò)和Kriging 模型進行了比較研究，對橋式起重機的進行輕量化設(shè)計。通過對3 種代理模型在確定系數(shù)R2的誤差分析，按子結(jié)構(gòu)進行分塊代理達(dá)到了更好的代理效果，根據(jù)不同代理模型的代理效果，提出了RSM-RBF 組合代理模型。最后，利用一種智能與梯度組合算法對3 種代理模型進行優(yōu)化，并將優(yōu)化結(jié)果與有限元分析結(jié)果進行對比.結(jié)果顯示，RSM-RBF 組合代理模型有更好的誤差代理效果，且其尋優(yōu)結(jié)果使得質(zhì)量降低了15.59%，證明了改優(yōu)化策略具有與一定的工程價值。