于 潔 賈 璞 肖超妮 鄭曉暉

（西北大學(xué)生命科學(xué)學(xué)院）

優(yōu)化問(wèn)題是實(shí)際的工程和日常生活中經(jīng)常會(huì)遇到的問(wèn)題，隨著社會(huì)發(fā)展和技術(shù)進(jìn)步，優(yōu)化問(wèn)題所涉及到的因素和需求越來(lái)越多，越來(lái)越復(fù)雜。 許多優(yōu)化問(wèn)題往往需要同時(shí)考慮多個(gè)目標(biāo)或?qū)傩?，屬于多目?biāo)優(yōu)化問(wèn)題。 多目標(biāo)優(yōu)化問(wèn)題，在一定程度上應(yīng)該歸屬于系統(tǒng)科學(xué)類問(wèn)題。 以系統(tǒng)論的觀點(diǎn)，將多目標(biāo)問(wèn)題視為一個(gè)系統(tǒng)，進(jìn)行整體優(yōu)化，才能使系統(tǒng)在整體上達(dá)到最優(yōu)，而各個(gè)目標(biāo)之間在系統(tǒng)的整體優(yōu)化的旗幟下達(dá)到某種協(xié)同或者折中。

目前，常用的多目標(biāo)算法，主要有線性加權(quán)法、帕累托解法和ε-約束法，以及機(jī)械設(shè)計(jì)和選材的層次分析法 （AHP）， 多準(zhǔn)則妥協(xié)解排序法（VIKOR），理想解相似度排序法（TOPSIS），基于比值分析的多目標(biāo)優(yōu)化法（MOORA）等。 但這些方法仍然存在一些基本問(wèn)題，如：線性加權(quán)法中涉及到的權(quán)重因子和歸一化因子的確定依據(jù)不明確； 帕累托解法只能給出解的一個(gè)集合；ε-約束法只挑選出其中某一目標(biāo)進(jìn)行優(yōu)化，而其他目標(biāo)則作為約束條件進(jìn)行處理，也就是將多目標(biāo)優(yōu)化問(wèn)題降級(jí)為單一目標(biāo)的優(yōu)化問(wèn)題等。 此外，從集合論和概率論的角度來(lái)看，線性加權(quán)法就是對(duì)幾個(gè)目標(biāo)的集合進(jìn)行并集運(yùn)算，這顯然不符合多個(gè)目標(biāo)同時(shí)優(yōu)化應(yīng)該是幾個(gè)目標(biāo)的集合之交集運(yùn)算的科學(xué)內(nèi)涵。 因此，要解決多目標(biāo)優(yōu)化問(wèn)題，就需要以系統(tǒng)論的觀點(diǎn)從整體上深刻揭示多目標(biāo)優(yōu)化的本征內(nèi)涵，以各個(gè)目標(biāo)之間的內(nèi)在關(guān)系為導(dǎo)向，從系統(tǒng)工程的角度，探討多目標(biāo)優(yōu)化的異質(zhì)同型的系統(tǒng)性，建立普適的各個(gè)多目標(biāo)之間統(tǒng)籌兼顧的系統(tǒng)模型，并在應(yīng)用中加以發(fā)展。

近年來(lái)，ZHENG M S等以系統(tǒng)論的觀點(diǎn)出發(fā)，從集合論和概率論的角度分析了多目標(biāo)優(yōu)化問(wèn)題，建立了全新的概率基多目標(biāo)優(yōu)化方法［1］。筆者將概率基多目標(biāo)優(yōu)化方法應(yīng)用于化工過(guò)程及機(jī)械設(shè)計(jì)的多目標(biāo)優(yōu)化問(wèn)題，必將有助于恰當(dāng)解決化工機(jī)械設(shè)計(jì)和加工方面的問(wèn)題，得到科學(xué)和理性的結(jié)果。

1 概率基多目標(biāo)同時(shí)優(yōu)化方法的正交試驗(yàn)設(shè)計(jì)步驟

1.1 概率基多目標(biāo)優(yōu)化方法簡(jiǎn)介

概率基多目標(biāo)優(yōu)化方法（Probability-based Multi-Objective Optimization，PMOO） 企望解決以往多目標(biāo)優(yōu)化方法中存在人為（或主觀）因素的固有問(wèn)題［1］。該方法以系統(tǒng)論的觀點(diǎn)，引入了青睞概率這一概念，來(lái)反映候選對(duì)象在優(yōu)化過(guò)程中被青睞的程度，并將候選對(duì)象的所有性能效用指標(biāo)劃分為效益型和成本型兩種基本類型。 效益型指標(biāo)具有越大越好的特征，成本型指標(biāo)具有越小越好的特征。 各候選對(duì)象的每個(gè)性能效用指標(biāo)都定量地貢獻(xiàn)出一個(gè)分青睞概率。 從概率論的觀點(diǎn)來(lái)看，候選對(duì)象的總青睞概率就是其所有分青睞概率的乘積，這是它在優(yōu)選過(guò)程中整體性和唯一的決定性指標(biāo)。

在新的概率基多目標(biāo)優(yōu)化方法中，候選對(duì)象的效益型性能指標(biāo)的效用值以正線性相關(guān)的方式貢獻(xiàn)其分青睞概率，即：

其中，Xij是第i個(gè)候選對(duì)象的第j個(gè) （效益型）性能指標(biāo)的效用值；Pij為該效益型性能指標(biāo)的效用值Xij的分青睞概率；n是相關(guān)候選對(duì)象組中候選對(duì)象的總數(shù)；m是該組中每個(gè)候選對(duì)象的性能指標(biāo)（目標(biāo)）的總數(shù)；αj是候選對(duì)象的第j個(gè)效用性能指標(biāo)的歸一化因子，αj=1/（nμj），μj是該效益型性能指標(biāo)的效用值在所涉及的對(duì)象組中的算術(shù)平均值。

對(duì)于成本型性能指標(biāo)， 可以進(jìn)行對(duì)等地處理，以負(fù)線性相關(guān)的方式貢獻(xiàn)其分青睞概率，即：

其中，Xjmax和Xjmin分別表示對(duì)象組中該性能指標(biāo)效用值的最大值和最小值，βj是第j個(gè)性能指標(biāo)的效用值的歸一化因子，且βj=1/［n（Xjmax+Xjminμj）］。

此外，從概率論的角度看，對(duì)于多目標(biāo)的“同時(shí)優(yōu)化”，第i個(gè)候選對(duì)象的總（綜合）青睞概率Pi是其所有分青睞概率Pij的乘積，即：

候選對(duì)象的總青睞概率在整個(gè)優(yōu)選過(guò)程中作為唯一的決定性指標(biāo)，通過(guò)它就將多目標(biāo)優(yōu)化問(wèn)題轉(zhuǎn)化為一個(gè)單目標(biāo)優(yōu)化問(wèn)題。 新的概率基多目標(biāo)優(yōu)化的主要特點(diǎn)是對(duì)于效益型效用指標(biāo)和成本型效用指標(biāo)進(jìn)行對(duì)等處理，而且在整個(gè)處理過(guò)程中不涉及人為（或主觀）的標(biāo)度因子。

1.2 概率基多目標(biāo)優(yōu)化方法與正交試驗(yàn)設(shè)計(jì)的結(jié)合

由于通過(guò)候選方案的總青睞概率能夠?qū)⒍嗄繕?biāo)優(yōu)化問(wèn)題轉(zhuǎn)化為單目標(biāo)優(yōu)化問(wèn)題，而且按照概率論，總青睞概率就是多個(gè)目標(biāo)同時(shí)優(yōu)化中唯一的、全面的決定性指標(biāo)，因此就可以對(duì)總青睞概率進(jìn)行正交試驗(yàn)設(shè)計(jì)的極差分析，再由極差分析所得到的結(jié)果來(lái)確定相應(yīng)的最佳配置。 這個(gè)操作程序就構(gòu)建了概率基多個(gè)目標(biāo)同時(shí)優(yōu)化的正交試驗(yàn)設(shè)計(jì)方法。

1.3 概率基多目標(biāo)優(yōu)化方法的穩(wěn)健性設(shè)計(jì)

早在1950年，田口玄一就意識(shí)到可以通過(guò)穩(wěn)健性（Robust，魯棒）設(shè)計(jì)來(lái)提高產(chǎn)品的質(zhì)量，并且提出可以通過(guò)降低噪聲影響的方式來(lái)進(jìn)行優(yōu)化，這就是田口方法［2］。在田口方法中，將影響因素分為可控和不可控因素，他指出可以通過(guò)設(shè)計(jì)實(shí)驗(yàn)來(lái)研究可控和不可控因素對(duì)響應(yīng)的影響。 田口玄一所稱的不可控因素就是噪聲因素。 穩(wěn)健性設(shè)計(jì)的思想是關(guān)于一組可控因素的設(shè)計(jì)，在最優(yōu)點(diǎn)處產(chǎn)品的質(zhì)量對(duì)所謂的噪聲因素具有不敏感性，或具有最小的噪聲影響。 在田口方法中，還進(jìn)一步假設(shè)可控因素主要指那些實(shí)驗(yàn)者或生產(chǎn)者容易控制的因素， 比如在注射成型過(guò)程中材料的選擇、注射成型時(shí)間及模具溫度等，而噪聲因素是那些不容易或非常昂貴或不可能控制的因素。 因此， 穩(wěn)健性設(shè)計(jì)就是尋找到一套產(chǎn)品和工藝參數(shù)，使在該參數(shù)處產(chǎn)品的最終質(zhì)量響應(yīng)指標(biāo)對(duì)不可控因素的變化具有最小敏感性，而不需要設(shè)法消除不可控因素。 進(jìn)一步，田口玄一引入了“信噪比（SNR）”一詞來(lái)實(shí)施其穩(wěn)健性設(shè)計(jì)，可控因素的最佳設(shè)計(jì)值就對(duì)應(yīng)于信噪比最大值的狀態(tài)。 田口玄一還給出了3種標(biāo)準(zhǔn)類型的信噪比。 對(duì)于以目標(biāo)值的期望值為優(yōu)化點(diǎn)的情形有：

然而，文獻(xiàn)［3～7］的研究表明，在實(shí)際的實(shí)驗(yàn)過(guò)程中， 通常實(shí)驗(yàn)測(cè)試到的結(jié)果的平均值μ和標(biāo)準(zhǔn)差σ是一組互相獨(dú)立的響應(yīng)。 但是，式（4）中信噪比SNR的表達(dá)式將這兩個(gè)響應(yīng)（μ和σ）固化成為一個(gè)信噪比（SNR）響應(yīng)，這是不合理的。 并且，SNR的最大值的優(yōu)化不等價(jià)于同時(shí)使σ達(dá)到最小、使μ接近于其目標(biāo)的優(yōu)化。 更為嚴(yán)重的是，在具有“越小越好”和“越大越好”特征問(wèn)題的表達(dá)式中，甚至標(biāo)準(zhǔn)差σ并未出現(xiàn)。這一點(diǎn)經(jīng)常受到統(tǒng)計(jì)學(xué)家的批評(píng)［3～7］。

另外，統(tǒng)計(jì)學(xué)家們建議，可以通過(guò)使用兩個(gè)單獨(dú)的模型來(lái)考慮均值μ和標(biāo)準(zhǔn)差σ的響應(yīng)。 因此， 對(duì)于σ的最小值和μ接近于其目標(biāo)的優(yōu)化，應(yīng)該同時(shí)采用單獨(dú)的模型來(lái)處理，以達(dá)到合理的穩(wěn)健性優(yōu)化設(shè)計(jì)之目的。

這里，按照概率論的方法，以候選對(duì)象（方案）性能指標(biāo)的算術(shù)平均值及其標(biāo)準(zhǔn)差作為處理中的兩個(gè)獨(dú)立響應(yīng)，便可實(shí)施穩(wěn)健性設(shè)計(jì)和產(chǎn)品的工藝優(yōu)化［1］。

2 概率基多目標(biāo)優(yōu)化方法在化工機(jī)械設(shè)計(jì)中的應(yīng)用實(shí)例

2.1 應(yīng)用于管殼式換熱器折流板結(jié)構(gòu)參數(shù)優(yōu)化的正交試驗(yàn)設(shè)計(jì)

管殼式換熱器廣泛應(yīng)用于石油、化工及冶煉等領(lǐng)域，具有諸多優(yōu)點(diǎn)，如結(jié)構(gòu)簡(jiǎn)單、安全可靠及傳熱效率高等。 管殼式換熱器中的折流板是不可缺少的工作部件，它會(huì)影響換熱器的傳熱、流動(dòng)特性，以及殼程壓降等行為。

李德濤等對(duì)折流板的缺口高度h在殼程圓筒內(nèi)徑的20%～45%范圍內(nèi)進(jìn)行設(shè)計(jì)，發(fā)現(xiàn)折流板缺口高度的變化會(huì)對(duì)流體的流動(dòng)狀態(tài)產(chǎn)生影響，進(jìn)而對(duì)其傳熱效率造成影響，而且換熱系數(shù)和殼程壓降會(huì)隨著折流板缺口高度的增大而減小，折流板的最小間距應(yīng)不小于殼體內(nèi)徑的1/5 也不應(yīng)小于50 mm。研究表明，折流板間距越大殼程傳熱系數(shù)越小、壓降越小。 實(shí)際上，折流板間距應(yīng)該有其最佳的參數(shù)值，故針對(duì)管殼式換熱器進(jìn)行了優(yōu)化設(shè)計(jì)［8］。 管殼式換熱器殼程內(nèi)徑為260 mm、殼程總長(zhǎng)1 500 mm、進(jìn)出口管內(nèi)徑90 mm、換熱管外徑24 mm、換熱管數(shù)量24根，換熱管間距36 mm。 將折流板缺口高度h、折流板數(shù)量N和折流板開(kāi)孔孔徑d作為輸入自變量， 每個(gè)因素有3個(gè)設(shè)計(jì)水平，進(jìn)行正交試驗(yàn)設(shè)計(jì)，以換熱器的殼程壓降和表面換熱系數(shù)作為目標(biāo)，進(jìn)行優(yōu)化設(shè)計(jì)。 對(duì)管殼式換熱器， 采用三維簡(jiǎn)化模型以及CFD進(jìn)行了模擬計(jì)算。 當(dāng)換熱介質(zhì)的入口流速為4 m/s時(shí)，其正交試驗(yàn)設(shè)計(jì)結(jié)果列于表1。

表1 管殼式換熱器多因素正交試驗(yàn)方案L9（34）及結(jié)果

按照概率基多目標(biāo)優(yōu)化方法對(duì)表1中數(shù)據(jù)進(jìn)行處理，得到表2的結(jié)果。 其中，壓降△P具有越低越好的特征，屬于成本型指標(biāo)；表面換熱系數(shù)η具有越大越好的特征，屬于效益型指標(biāo)。 從表2結(jié)果可以看出， 方案2具有最大的總青睞概率，因此，優(yōu)化結(jié)果應(yīng)該在方案2附近。 進(jìn)一步，對(duì)表2給出的總青睞概率進(jìn)行極差分析， 得到表3的結(jié)果。

表2 管殼式換熱器多目標(biāo)試驗(yàn)結(jié)果的評(píng)價(jià)情況

表3 管殼式換熱器概率基多目標(biāo)評(píng)價(jià)極差分析

極差分析表明，因素的重要性次序?yàn)閔＞N＞d。而最終優(yōu)化的結(jié)果為h1N1d3， 即折流板缺口高度h=95 mm、折流板數(shù)量N=4塊、折流板開(kāi)孔孔徑d=10 mm時(shí)，換熱器換熱效果可達(dá)到最優(yōu)。 然而，李德濤等采用綜合性能評(píng)價(jià)因子法，獲得的結(jié)果是h=95 mm、折流板數(shù)量N=4塊、折流板開(kāi)孔孔徑d=8 mm［8］，異于本文的結(jié)果。

2.2 應(yīng)用于降黏聚結(jié)一體化裝置的優(yōu)化設(shè)計(jì)

卜凡熙等針對(duì)旋流器處理含聚介質(zhì)的“消化不良”現(xiàn)象，設(shè)計(jì)出一種降黏聚結(jié)一體化裝置。 將該裝置串聯(lián)于聚介質(zhì)的旋流處理器之前端，在降低聚合物溶液黏度的同時(shí)， 還增大了油滴粒徑，從而降低了含聚介質(zhì)黏度和油滴乳化程度對(duì)旋流器分離性能的影響。 同時(shí)，為了提高裝置的降黏聚結(jié)效果，采用正交試驗(yàn)方法對(duì)其進(jìn)行了結(jié)構(gòu)尺寸優(yōu)化。試驗(yàn)選取A、B、C、D、E、F、G7個(gè)因素，每個(gè)因素設(shè)定3個(gè)水平，并以油滴粒徑、黏度和壓降為3個(gè)優(yōu)化目標(biāo)，進(jìn)行了優(yōu)化設(shè)計(jì)。 采用了三維簡(jiǎn)化模型以及CFD進(jìn)行模擬計(jì)算， 其正交設(shè)計(jì)結(jié)果列于表4［9］。

表4 降黏聚結(jié)一體化裝置的正交試驗(yàn)L18（2137）方案及結(jié)果

按照概率基多目標(biāo)優(yōu)化方法對(duì)表4中數(shù)據(jù)進(jìn)行處理，得到的結(jié)果列于表5，其中，油滴粒徑為效益型指標(biāo)，黏度和壓降為成本型指標(biāo)。從表5可以看出， 試驗(yàn)方案6具有最大的總青睞概率，因此，優(yōu)化結(jié)果應(yīng)該在方案6附近。 進(jìn)一步，對(duì)表5給出的總青睞概率進(jìn)行極差分析， 得到表6的結(jié)果。

表5 降黏聚結(jié)一體化裝置多目標(biāo)試驗(yàn)結(jié)果的評(píng)價(jià)情況

表6 降黏聚結(jié)一體化裝置概率基多目標(biāo)評(píng)價(jià)極差分析

極差分析結(jié)果表明， 因素的重要性次序?yàn)镃＞F＞A＞D＞G＞E＞B， 最終優(yōu)化結(jié)果為A2B3C3D2E2F2G3時(shí)，降黏聚結(jié)一體化裝置可達(dá)到最優(yōu)。卜凡熙等采用“綜合平衡法”處理了這一問(wèn)題［9］，得到的最終優(yōu)化結(jié)果為A2B2C3D2E2F2G1，和本文結(jié)果有點(diǎn)差異。

2.3 應(yīng)用于等距螺旋錐齒輪MIM工藝參數(shù)多目標(biāo)穩(wěn)健性優(yōu)化

劉贛華等針對(duì)等距螺旋錐齒輪， 利用Moldex3D軟件和正交試驗(yàn)設(shè)計(jì)進(jìn)行了金屬粉末注射成型（MIM）工藝成形的仿真分析。 以體積收縮率和粉末濃度為目標(biāo)函數(shù)，試圖通過(guò)對(duì)目標(biāo)函數(shù)均值分析和極差的分析，并利用BP神經(jīng)網(wǎng)絡(luò)結(jié)合非支配排序遺傳算法（NSGA-Ⅱ）進(jìn)行MIM工藝參數(shù)的優(yōu)化。 MIM工藝參數(shù)優(yōu)化考慮了5個(gè)影響因子，即保壓壓力A、充填時(shí)間B、澆口直徑C、熔體溫度D和模具溫度E，每個(gè)因子設(shè)計(jì)了5個(gè)水平，并采用L25（55）進(jìn)行正交試驗(yàn)設(shè)計(jì)［10］。 其正交設(shè)計(jì)結(jié)果列于表7。其中，μ和σ分別表示體積收縮率與粉末濃度的平均值和標(biāo)準(zhǔn)差；合金粉末的輸入體積分?jǐn)?shù)為60%。

表7 MIM工藝參數(shù)的L25（55）正交試驗(yàn)方案及體積收縮率與粉末濃度結(jié)果

筆者針對(duì)這一問(wèn)題，采用概率基多目標(biāo)穩(wěn)健性評(píng)價(jià)［1］，由于粉末濃度的輸入值為60%，故其效用值應(yīng)為ε=|60%-μv|，具有越小越好的特征；同時(shí)體積收縮率的平均值和標(biāo)準(zhǔn)差、粉末濃度的標(biāo)準(zhǔn)差也具有越小越好的特征。按照概率基多目標(biāo)優(yōu)化方法的穩(wěn)健性評(píng)價(jià)方法，以及粉末濃度的效用值ε的定義，對(duì)表7中數(shù)據(jù)進(jìn)行處理得到表8的結(jié)果。

表8 MIM工藝參數(shù)試驗(yàn)結(jié)果的多目標(biāo)穩(wěn)健性評(píng)價(jià)情況

以Pt=（Pμv·Pσv）×（Pμd·Pσd）為依據(jù)進(jìn)行總青睞概率Pt的評(píng)估［1］。 由表8可知，試驗(yàn)方案25（A5B5C4D3E2）具有最大的總青睞概率，因此，優(yōu)化結(jié)果應(yīng)該在方案25附近。 進(jìn)行極差分析得到表9結(jié)果。

表9 MIM工藝參數(shù)試驗(yàn)結(jié)果的概率基多目標(biāo)評(píng)價(jià)極差分析

極差分析結(jié)果顯示，因素的重要性次序?yàn)镃＞B＞A＞D＞E，最終優(yōu)化結(jié)果為A5B1C4D3E2，即保壓壓力A=106 MPa、 充填時(shí)間B=0.20 s、 澆口直徑C=0.55 mm、熔體溫度D=240 ℃和模具溫度E=40 ℃。劉贛華等采用“BP神經(jīng)網(wǎng)絡(luò)結(jié)合非支配排序遺傳算 法（NSGA-Ⅱ）”處 理 了 這 一 問(wèn) 題［10］，得 到 的Pareto最優(yōu)解為：保壓壓力A=99.183 MPa、充填時(shí)間B=0.241 s、澆口直徑C=0.548 mm、熔體溫度D=240.702 ℃、模具溫度E=63.981 ℃，有異于本文的結(jié)果。

上述3個(gè)案例的計(jì)算結(jié)果均與本文結(jié)果存在差異，其原因在于概率基多目標(biāo)優(yōu)化屬于系統(tǒng)的優(yōu)化方法，期待之后有學(xué)者對(duì)計(jì)算效果進(jìn)行驗(yàn)證評(píng)價(jià)。

3 結(jié)束語(yǔ)

將概率基多目標(biāo)優(yōu)化方法和正交試驗(yàn)設(shè)計(jì)，以及穩(wěn)健性設(shè)計(jì)結(jié)合起來(lái)，可以得到概率基多目標(biāo)正交試驗(yàn)優(yōu)化設(shè)計(jì)和穩(wěn)健性設(shè)計(jì)方法。 將概率基多目標(biāo)正交試驗(yàn)優(yōu)化設(shè)計(jì)方法應(yīng)用于管殼式換熱器折流板結(jié)構(gòu)參數(shù)優(yōu)化的正交試驗(yàn)設(shè)計(jì)和降黏聚結(jié)一體化裝置的優(yōu)化設(shè)計(jì)，概率基多目標(biāo)穩(wěn)健性設(shè)計(jì)方法應(yīng)用于等距螺旋錐齒輪MIM工藝參數(shù)多目標(biāo)的穩(wěn)健性設(shè)計(jì)，得到了多個(gè)目標(biāo)同時(shí)優(yōu)化的結(jié)果。 該方法簡(jiǎn)便易行，可望在廣泛領(lǐng)域得到應(yīng)用。