馬睿佳，魏正英，陳雪麗，馬勝利(西安交通大學(xué)機(jī)械制造系統(tǒng)工程國家重點實驗室，西安 710049)

迷宮型灌水器具有結(jié)構(gòu)相對簡單、制造方便、成本低，又具有一定的壓力補(bǔ)償性等諸多優(yōu)點而得到廣泛應(yīng)[1]。迷宮型灌水器主要由3部分組成：進(jìn)口柵格、迷宮流道和出口。其中迷宮流道利用復(fù)雜的壁面邊界使其中的水流紊亂，從而引起能量損失，水流在流道內(nèi)的流動近似表現(xiàn)為紊流，是此類灌水器消能的主要部位，它的結(jié)構(gòu)特性在很大程度上決定著灌水器的2大性能：水力性能和抗堵性能。

傳統(tǒng)的求解方法需要對多個目標(biāo)進(jìn)行合并,需要選取適當(dāng)?shù)恼壑邢禂?shù),有很強(qiáng)的先驗性,增加了決策者決策的難度，如灰靶理論[2]等。為了更好地解決復(fù)雜的多目標(biāo)優(yōu)化問題，研究人員引進(jìn)了多目標(biāo)進(jìn)化算法，其中，粒子群算法，也稱粒子群優(yōu)化算法(Particle Swarm Optimization 簡稱PSO)， 是近年來發(fā)展起來的一種新的進(jìn)化算法。PSO 算法是從隨機(jī)解出發(fā)，通過多次迭代尋找最優(yōu)解也稱為非劣解，它也是通過適應(yīng)度來評價解的品質(zhì)，它通過追隨當(dāng)前搜索到的最優(yōu)值來尋找全局最優(yōu)[3]。本文通過對兩大性能的定量分析，在正交實驗的基礎(chǔ)上，采用粒子群算法對彎道環(huán)流式流道進(jìn)行多目標(biāo)優(yōu)化。

1 彎道環(huán)流式流道結(jié)構(gòu)設(shè)計

1.1 設(shè)計依據(jù)

彎道環(huán)流現(xiàn)象是指在2個連續(xù)且相反的90°圓弧連接處，由于前后兩端壓強(qiáng)的作用剛好相反，會形成2個轉(zhuǎn)向相反的螺旋狀前進(jìn)流，在兩個圓弧連接處垂向截面有雙環(huán)流現(xiàn)象[4]。流體的能量損失包括局部損失和沿程損失，對于灌水器的流道而言，由于內(nèi)部流道的長度相對較短，能量損失主要為局部損失。利用彎道環(huán)流現(xiàn)象可以增加流動內(nèi)部的碰撞和漩渦，增加流態(tài)內(nèi)部的紊亂程度，從而有效耗能。借鑒上述現(xiàn)象設(shè)計彎道環(huán)流式流道如圖1所示。

圖1 彎道環(huán)流式流道Fig.1 Bend circulating current emitter 注：1-出口；2-進(jìn)口；3-流道；4-直線過渡段；5-寬彎段；6-圓弧漸縮段；7-細(xì)彎段。

1.2 流道結(jié)構(gòu)參數(shù)化

流道單元包括直線過渡段、寬彎段、圓弧漸縮段及細(xì)彎段4部分，如圖1所示。寬彎段和細(xì)彎段分別為2個同心圓弧構(gòu)成，且2個部分的圓心位于同一水平上，直線過渡段與2段圓弧外圍進(jìn)行相同大小的圓弧連接，圓弧漸縮段為過渡輔助。流道單元設(shè)計為寬進(jìn)窄出，在圓弧漸縮段產(chǎn)生漸縮、碰撞及本身的螺旋渦旋損失，進(jìn)入直線過渡段后又產(chǎn)生漸擴(kuò)損失。如圖2所示，提取流道中結(jié)構(gòu)參數(shù)，寬彎段的內(nèi)徑為R1，直線過渡段兩端的過渡圓弧半徑都為R2，細(xì)彎段的內(nèi)徑為R3，直線過渡段寬度為H，2段圓弧的外徑保持不變，寬彎段外徑為1.3 mm，細(xì)彎段的外徑為1.4 mm，漸縮段為與寬彎段相切的且必須通過A點的圓弧。

圖2 流道單元結(jié)構(gòu)Fig.2 Channel unit of emitter

2 流道結(jié)構(gòu)正交實驗設(shè)計

2.1 正交實驗設(shè)計

選擇流道參數(shù)R1、R2、R3、H共4個因素，不考慮個因素之間的交互作用，加入誤差項，選用5因素4水平L16(45)正交表，試驗因素水平表見表1，正交實驗設(shè)計見表2。共16組，采用數(shù)值模擬的方法，通過建模-網(wǎng)格劃分-FLUENT計算的方式對流道性能參數(shù)進(jìn)行提取分析，建立流道結(jié)構(gòu)參數(shù)與灌水器2大性能之間的目標(biāo)函數(shù)。

表1 正交實驗因素水平 mmTab.1 Factors and levels of orthogonal experiment

表2 L16(45)設(shè)計及試驗結(jié)果Tab.2 L16(45) design table and experiment results

2.2 數(shù)值模型

灌水器流道深度為1 mm,流道進(jìn)出口采用六面體結(jié)構(gòu)網(wǎng)格劃分，流道部分由于多圓弧部分，采用四面體網(wǎng)格劃分，網(wǎng)格密度為0.1 mm,總網(wǎng)格數(shù)為58萬個，細(xì)化網(wǎng)格密度分別為0.08和0.06 mm，統(tǒng)計質(zhì)量流誤差不超過1%，但計算耗時分別增長至2～3倍。流道進(jìn)口為壓力進(jìn)口，出口為大氣壓。模型采用標(biāo)準(zhǔn)k-ε模型，近壁面選用標(biāo)準(zhǔn)壁面函數(shù)，采用SIMPLE算法，二階離散格式進(jìn)行離散，收斂精度為10-4。采用離散相模型進(jìn)行流道內(nèi)部固液兩相流的研究[4]。

2.3 正交實驗結(jié)果分析

取流態(tài)指數(shù)作為水力性能的衡量指標(biāo)，粒子通過率作為抗堵性能的指標(biāo)，進(jìn)行16組數(shù)值模擬，分別采用直觀比較分析、方差分析、線性回歸分析等方法對表2進(jìn)行分析，并得到結(jié)構(gòu)參數(shù)與2個指標(biāo)的回歸模型。

(1)流態(tài)指數(shù)分析。通過直觀分析可知各結(jié)構(gòu)對流態(tài)參數(shù)影響為R1>R3>R2>H，根據(jù)流態(tài)指數(shù)因素-指標(biāo)圖3，選擇流態(tài)指數(shù)的極小值，可知各因素的最佳組合為R1=0.6 mm，R2=0.65 mm，R3=0.71 mm，H=0.85 mm，此組合沒有在正交組合中出現(xiàn)。對其進(jìn)行方差分析，對α=0.1，查表可知F0.1(3,3)=5.39，結(jié)果如表3所示，可知R1對結(jié)果的影響一般顯著，其他參數(shù)對結(jié)果影響不顯著。用SPSS軟件對數(shù)據(jù)進(jìn)行回歸分析，得到流態(tài)指數(shù)與各結(jié)構(gòu)參數(shù)的回歸模型為：

X=3.022-1.511R1-0.989R2-2.325R3-2.239H+

1.243R21+0.781R22+1.667R23+1.319H2

(1)

表3 流態(tài)指數(shù)方差分析結(jié)果Tab.3 Variance analysis results of flow index

圖3 流態(tài)指數(shù)因素-指標(biāo)圖Fig.3 The factor-index chart of flow index

模型的相關(guān)系數(shù)R=0.965，R2=0.931，線性關(guān)系顯著。將通過回歸模型計算的流態(tài)指數(shù)與正交試驗的結(jié)果進(jìn)行比較，二者的相對誤差在2.6%以內(nèi)，回歸模型精度較高。

(2)粒子通過率分析。同理，通過直觀分析可知流道結(jié)構(gòu)參數(shù)對對粒子通過率的影響順序為H>R1>R2>R3。具有最大的粒子通過率的結(jié)構(gòu)組合為：R1=0.6 mm，R2=0.60 mm，R3=0.74 mm，H=0.85 mm。在相關(guān)系數(shù)R方為0.878的顯著水平下，回歸分析得到的粒子通過率與流道結(jié)構(gòu)參數(shù)的回歸模型為:

K=3.909-3.215R1-0.067R2+3.482R3-8.602H+

3.09R21+0.133R22-2.243R23+5.168H2

(2)

3 基于粒子群算法的流道結(jié)構(gòu)優(yōu)化

3.1 多目標(biāo)粒子群算法

粒子群算法是一種隨機(jī)的基于種群的優(yōu)化策略，通過模擬鳥群覓食的過程建立的一種優(yōu)化算法。粒子群算法中，每只鳥作為一個粒子，每個粒子有位置和速度兩個屬性，在飛行過程中根據(jù)自身找到的距離食物最近的位置和整個群體找到最近位置去不斷改變自身的前進(jìn)方向，最后整個群體都會趨向同一個地方，即食物的位置。在原有的粒子群算法基礎(chǔ)上，Shi[6]等提出了慣性權(quán)重粒子群算法，也稱標(biāo)準(zhǔn)的粒子群算法，包括了3個部分，粒子自身的慣性部分、自身認(rèn)知部分和社會部分。該算法的速度更新公式如下：

vid(t+1)=ωvid(t)+c1ri1,t[Xpbid(t)-xid(t)]+

c2ri2,t[Xgbid(t)-xid(t)]

(3)

式中：ω為慣性權(quán)重，描述粒子慣性對當(dāng)前速度的影響；c1、c2為加速系數(shù)；ri1,t、ri2,t為0到1之間的隨機(jī)數(shù)。位置更新公式如下：

xid(t+1)=xid(t)+vid(t+1)

(4)

之后，Clerc[7]等提出了一種帶有收縮因子的粒子群算法，在保證算法收斂的同時提高算法的收斂速度。本文采用這種算法對速度進(jìn)行更新，該算法的描述如下：

vid(t+1)=χ{vid(t)+c1ri1,t[Xpbid(t)-xid(t)]+

c2ri2,t[Xgbid(t)-xid(t)]}

計算取c1=c2=2.05，則φ=4.1，χ=0.729。由式(3)和式(5)可知，2種算法從表達(dá)形式來說基本一致，式(5)是式(3)的變形。

多目標(biāo)粒子群算法(multiobjective particle swarm optimization,MOPSO)中以2004年Coello[8]提出基于精英級策略的MOPSO算法最具有里程碑意義，本文采用該算法進(jìn)行灌水器結(jié)構(gòu)的優(yōu)化計算。該算法的具體步驟如下。

(1)隨機(jī)初始化種群，計算各個粒子的適應(yīng)度值，將其中的非劣解及其位置信息存入外部檔案庫REP。

(2)確定粒子初始的個體最優(yōu)Xpbi(0)和全局最優(yōu)Xgbi(0)。

(3)迭代次數(shù)增加，每迭代一次，更新粒子的速度和位置，其中個體最優(yōu)位置通過比較單個個體經(jīng)歷過的歷史最優(yōu)位置，同時為每一個粒子選取全局最優(yōu)，全局最優(yōu)采用自適應(yīng)網(wǎng)格法和輪盤選擇法[9]從外部檔案庫REP中選取一個。

(4)迭代次數(shù)增加，直至滿足終止條件。最終外部檔案庫REP保存的就是所需要的最優(yōu)解。

外部檔案庫REP的作用機(jī)制包括以下幾個方面：首先，每一次迭代得到的全部粒子中的非劣解與上一次迭代得到的REP內(nèi)的值進(jìn)行比較，若新的非劣解被REP內(nèi)原有的任一解支配，則無法進(jìn)入REP內(nèi)；若REP內(nèi)的所有值都無法支配新的非劣解，則新的非劣解進(jìn)入；若新的非劣解可以支配其中的某個非劣解，則原來的非劣解丟棄，新的非劣解進(jìn)入REP內(nèi)。當(dāng)外部檔案庫REP達(dá)到數(shù)量上限時，采取二次準(zhǔn)則對里面的元素進(jìn)行保留，即網(wǎng)格內(nèi)稀疏位置的元素具有優(yōu)先保留權(quán)。

3.2 基于MOPSO的彎道環(huán)流式流道結(jié)構(gòu)優(yōu)化

灌水器流道多目標(biāo)優(yōu)化的目標(biāo)函數(shù)由上文得到如下：

MinX=3.022-1.511R1-0.989R2-2.325R3-

2.239H+1.243R21+0.781R22+1.667R23+1.319H2

MaxK=3.909-3.215R1-0.067R2+3.482R3-

8.602H+3.09R21+0.133R22-2.243R23+5.168H2

其中，結(jié)構(gòu)參數(shù)的約束范圍為：

(6)

由于優(yōu)化目標(biāo)是找到具有優(yōu)良的水力性能及抗堵塞性能的小流量灌水器，故將流量作為額外的約束標(biāo)準(zhǔn)，流量回歸方程由上述正交實驗獲得：

Q=2.319-0.803R1-0.253R2-0.343R3+0.354H

1.5

(7)

其中設(shè)置種群規(guī)模為100個，迭代次數(shù)為100次，分別比較經(jīng)過初始代、30次、70次、100次迭代后的非劣解集如圖4所示，可知粒子通過率K隨著流態(tài)指數(shù)x的增大而增大，基本上保持為線性關(guān)系。而理想的結(jié)果為粒子通過率越大且流態(tài)指數(shù)越小越好，故可知兩者是相互抵制的，在得到非劣解集后，仍需要人為選擇兩者關(guān)系的權(quán)重來確定相應(yīng)的組合。分別在迭代第100次得到非劣解集中選取5個點，如圖5所示，并列出對應(yīng)的結(jié)構(gòu)參數(shù)的組合，如表4所示。

由表4可知，并對其中非劣解集中的點的結(jié)構(gòu)參數(shù)及目標(biāo)函數(shù)值分析可知，R1的大小基本為0.6，其他3個參數(shù)的可浮動量的大小分別為H最大，R2次之，R3最小。說明R1的取值對2大性能的影響很大，只有在取值為0.6的時候才可以得到最優(yōu)值，而非劣解集的前沿分布主要由H和R2的大小決定。同時，對比前面的正交實驗結(jié)果可知，前面的最佳組合基本圍繞在參數(shù)變量的附近，可知正交分析與多目標(biāo)粒子群算法的正確性。得到非劣解集后，應(yīng)當(dāng)根據(jù)設(shè)計需求 期、抽穗期、灌漿期的適宜灌溉控制下限標(biāo)準(zhǔn)(占田持%)分別為45%、65%、70%、65%，且在生產(chǎn)中盡量避免出現(xiàn)長期連續(xù)干旱現(xiàn)象。

圖4 迭代過程中非劣解集的進(jìn)化過程Fig.4 Evolutionary process of non dominated solution set in the iterative process

圖5 取點位置圖Fig.5 Location of point

表4 非劣解參數(shù)信息Tab.4 Parameters of non dominated solution

□

[1] 路振廣，孟春紅，王曉麗.基于遺傳算法的冬小麥優(yōu)化灌溉制度研究[J].人民黃河，2009,31(11)：76-78.

[2] 劉 鈺，汪 林，倪廣恒，等.中國主要作物灌溉需水量空間分布特征[J].農(nóng)業(yè)工程學(xué)報，2009,25(12)：6-12.

[3] 高 陽，黃 玲，李新強(qiáng)，等.開花后水分脅迫對冬小麥旗葉光合作用和保護(hù)酶活性的影響[J].水土保持學(xué)報，2013，27(4)：201-206.

[4] 申孝軍，孫景生，劉祖貴，等.灌水控制下限對冬小麥產(chǎn)量和品質(zhì)的影響[J].農(nóng)業(yè)工程學(xué)報，2011，26(12)：58-65.

[5] 徐建文，梅旭榮，居 輝，等.黃淮海地區(qū)冬小麥關(guān)鍵生育期不同灌溉水平對產(chǎn)量影響的模擬[J].作物學(xué)報，2014，40(8)：1 485-1 492.

[6] 姚 寧，宋利兵，劉 健，等.不同生長階段水分脅迫對旱區(qū)冬小麥生長發(fā)育和產(chǎn)量的影響[J].中國農(nóng)業(yè)科學(xué)，2015，48(12)：2 379-2 389.

[7] 王俊儒，李生秀.不同生育時期水分有限虧缺對冬小麥產(chǎn)量及其構(gòu)成因素的影響[J].西北植物學(xué)報，2000，20(2)：193-200.

[8] 閆永鑾，郝衛(wèi)平，梅旭榮，等.拔節(jié)期水分脅迫-復(fù)水對冬小麥干物質(zhì)積累和水分利用效率的影響[J].中國農(nóng)業(yè)氣象，2011，32(2)：190-195.

[9] 霍治國，李世奎，白月明，等.冬小麥中輕度水分脅迫的增產(chǎn)節(jié)水效應(yīng)研究[J].自然資源學(xué)報，2003，18(1)：58-66.

[10] 銀敏華，李援農(nóng)，周昌明，等.調(diào)虧灌水和分蘗干擾對冬小麥生長的補(bǔ)償效應(yīng)[J].應(yīng)用生態(tài)學(xué)報，2015，26(10)：3 011-3 019.