• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    一種汽輪機(jī)基礎(chǔ)的并行多目標(biāo)優(yōu)化方法

    2015-02-15 06:56:14軍,征,,
    關(guān)鍵詞:黑箱算例汽輪機(jī)

    李 召 軍, 李 征, , 李 克 秋

    ( 1.大連理工大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院, 遼寧 大連 116024;2.大連理工大學(xué) 工業(yè)裝備結(jié)構(gòu)分析國家重點(diǎn)實(shí)驗(yàn)室, 遼寧 大連 116024 )

    工程力學(xué)

    一種汽輪機(jī)基礎(chǔ)的并行多目標(biāo)優(yōu)化方法

    ( 1.大連理工大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院, 遼寧 大連 116024;2.大連理工大學(xué) 工業(yè)裝備結(jié)構(gòu)分析國家重點(diǎn)實(shí)驗(yàn)室, 遼寧 大連 116024 )

    提出一種汽輪機(jī)基礎(chǔ)的并行多目標(biāo)優(yōu)化方法,旨在搜尋能夠使汽輪機(jī)基礎(chǔ)的質(zhì)量和最大動位移取得最小值的設(shè)計(jì)變量的解.基礎(chǔ)的質(zhì)量可用其梁柱截面積的線性表達(dá)式表達(dá),而其在外加擾力作用下產(chǎn)生的動位移可由黑箱優(yōu)化方法計(jì)算獲?。贙riging替代函數(shù)模型建立的基礎(chǔ)梁柱截面積與動位移的近似表達(dá)式基礎(chǔ)上,發(fā)展了基于熵的期望提高準(zhǔn)則用于提高優(yōu)化的精度,并使用并行計(jì)算策略提高優(yōu)化的效率.對兩個實(shí)際汽輪機(jī)基礎(chǔ)算例優(yōu)化的結(jié)果表明,該方法可以有效解決該類工程優(yōu)化問題,得到比較理想的優(yōu)化解.

    汽輪機(jī)基礎(chǔ);多目標(biāo)優(yōu)化;Kriging;并行計(jì)算;期望提高

    0 引 言

    汽輪機(jī)基礎(chǔ)通常是包含梁柱結(jié)構(gòu)的混凝土基礎(chǔ)框架,其荷載一般包括基礎(chǔ)的自重以及基礎(chǔ)之上電機(jī)轉(zhuǎn)動所產(chǎn)生的動擾力.汽輪機(jī)基礎(chǔ)的設(shè)計(jì)原則是在降低基礎(chǔ)質(zhì)量的同時,盡量減少其產(chǎn)生的動位移.在傳統(tǒng)的嚴(yán)格設(shè)計(jì)中,柱的位置是固定的并且基礎(chǔ)頻率比擾力頻率要高很多;但現(xiàn)代的靈活設(shè)計(jì)中,柱的位置是可變的,并且基礎(chǔ)頻率也可以小于擾力頻率.現(xiàn)階段汽輪機(jī)基礎(chǔ)的體積越來越大,傳統(tǒng)意義上的嚴(yán)格設(shè)計(jì)已經(jīng)不再適用,靈活的基礎(chǔ)設(shè)計(jì)模式在當(dāng)今的汽輪機(jī)基礎(chǔ)設(shè)計(jì)領(lǐng)域占據(jù)了主導(dǎo)地位[1].

    近年來,工程優(yōu)化設(shè)計(jì)發(fā)展迅猛.將優(yōu)化技術(shù)應(yīng)用在汽輪機(jī)基礎(chǔ)設(shè)計(jì)時,梁柱的截面積及柱的幾何位置都可以作為設(shè)計(jì)變量,最終實(shí)現(xiàn)靈活設(shè)計(jì)的目的.通常這類工程優(yōu)化問題都是帶約束的多設(shè)計(jì)變量多目標(biāo)問題.而黑箱優(yōu)化方法的出現(xiàn),無疑為解決這類問題提供了全新的思路[2].

    20世紀(jì)70年代以來快速發(fā)展的高性能計(jì)算技術(shù)為資源密集型的工程優(yōu)化問題提供了便利.并行計(jì)算技術(shù)作為代表,在幾十年的發(fā)展過程中得到了長足的進(jìn)步,并出現(xiàn)了許多通用的編程框架和計(jì)算環(huán)境.如何將這些技術(shù)應(yīng)用在工程優(yōu)化問題特別是汽輪機(jī)基礎(chǔ)的優(yōu)化設(shè)計(jì)當(dāng)中以獲取最大的效率提高,是一個亟須解決的實(shí)際應(yīng)用問題[3].

    本文提出的并行多目標(biāo)優(yōu)化方法建立在黑箱優(yōu)化的基礎(chǔ)上,通過ANSYS軟件作為黑箱方法計(jì)算基礎(chǔ)在擾力下的動位移,使用Kriging代理模型模擬動位移和設(shè)計(jì)截面變量的函數(shù)關(guān)系;同時,設(shè)計(jì)一種基于熵的期望提高準(zhǔn)則用來添加優(yōu)化解,并使用并行計(jì)算技術(shù)完成整個優(yōu)化過程;最后對由兩個實(shí)際基礎(chǔ)組成的優(yōu)化算例進(jìn)行測試以驗(yàn)證其有效性.

    1 研究背景

    1.1 黑箱優(yōu)化方法

    通常情況下工程優(yōu)化應(yīng)用包含幾個相似的過程:(1)抽樣,采用抽樣算法根據(jù)設(shè)計(jì)變量上下限獲取一定規(guī)模的樣本組群;(2)分析,使用通用或特定的工程應(yīng)用軟件得到所有樣本的響應(yīng)值;(3)建模,在樣本和響應(yīng)值之間建立一個近似的數(shù)學(xué)模型(如果不使用替代函數(shù)法,該過程可省略);(4)優(yōu)化,采用某種優(yōu)化算法在給定的約束和設(shè)計(jì)變量上下限范圍內(nèi)得到最優(yōu)值;(5)計(jì)算最優(yōu)值樣本的目標(biāo)函數(shù)值,使用同樣的軟件計(jì)算該優(yōu)化值的響應(yīng)值;(6)收斂檢測,若優(yōu)化值滿足所有約束條件且其目標(biāo)函數(shù)值符合優(yōu)化要求,則該優(yōu)化值即為最終的優(yōu)化解,否則將該優(yōu)化值與其響應(yīng)值分別加入樣本組群和目標(biāo)函數(shù)組群,轉(zhuǎn)過程(3)[2,4].

    而黑箱優(yōu)化方法,顧名思義,就是將上述各個階段均采用相對獨(dú)立的模塊完成對應(yīng)的功能,上下模塊間使用標(biāo)準(zhǔn)的輸入輸出進(jìn)行參數(shù)的傳遞,并使用通用型的控制程序完成整個優(yōu)化過程.整個優(yōu)化應(yīng)用執(zhí)行期間,子模塊各司其職,就如同一個個不透明的黑箱,因此這種序列化和標(biāo)準(zhǔn)化的工程優(yōu)化過程也被廣泛稱為黑箱優(yōu)化.因抽樣、分析、建模等黑箱過程均在實(shí)際問題中有許多可以利用的商業(yè)或開源軟件及程序,黑箱優(yōu)化方法也被設(shè)計(jì)人員大量采用和推廣[4-5].

    1.2 通用工程分析軟件

    工程領(lǐng)域在應(yīng)用分析軟件的使用上逐漸形成了一系列通用的選擇:結(jié)構(gòu)設(shè)計(jì)中的ANSYS和ABAQUS;模具設(shè)計(jì)中的Moldflow和Unigraphics UX;建筑設(shè)計(jì)中的Revit、3DMax和AutoCAD等.近年來,越來越多的數(shù)據(jù)中心開始以集中服務(wù)形式(如云平臺)提供給用戶使用這些軟件的各項(xiàng)功能,在方便用戶使用的同時,也極大地節(jié)約了能源和購買成本.這些軟件豐富的應(yīng)用編程接口(API)進(jìn)一步使得用戶可針對具體的應(yīng)用開發(fā)出特定的優(yōu)化邏輯和優(yōu)化模型,極大地促進(jìn)了工程優(yōu)化領(lǐng)域的發(fā)展[6].

    1.3 優(yōu)化算法

    蟻群算法(ant colony optimization)、粒子群算法(particle swarm optimization)、遺傳算法(genetic algorithm)和序列二次規(guī)劃(sequential quadratic programming)、人工神經(jīng)網(wǎng)絡(luò)(artificial neural networks)算法等是近些年比較流行的優(yōu)化算法.具體問題不同,算法選擇也不一.這些算法對工程優(yōu)化領(lǐng)域的理論研究及實(shí)際應(yīng)用產(chǎn)生了重要的影響.對這些優(yōu)化算法的具體分析已經(jīng)超過了本文的范圍[7].

    1.4 并行計(jì)算技術(shù)

    20世紀(jì)70年代之后,諸如多核及多處理機(jī)計(jì)算機(jī)廣泛發(fā)展,并行計(jì)算技術(shù)也因此產(chǎn)生,用來充分利用這些硬件提升所帶來的計(jì)算便利.這期間出現(xiàn)了諸多并行算法與編程語言,消息處理接口(MPI)便是比較流行的一種[3].

    隨后,集群計(jì)算開始出現(xiàn):將硬件和軟件環(huán)境類似的物理機(jī)通過高性能網(wǎng)絡(luò)設(shè)備加以互聯(lián),輔以并行計(jì)算環(huán)境的配置,可被用來處理更多更復(fù)雜和更耗時的計(jì)算任務(wù).通常,集群通過批處理作業(yè)管理系統(tǒng)進(jìn)行計(jì)算任務(wù)的管理,如OpenPBS.組成一個集群的計(jì)算機(jī)硬件一般包括磁盤處理單元、控制節(jié)點(diǎn)、計(jì)算節(jié)點(diǎn)及網(wǎng)絡(luò)設(shè)備等[8].

    1.5 Kriging替代函數(shù)模型

    諸如梁柱截面等汽輪機(jī)基礎(chǔ)優(yōu)化中的設(shè)計(jì)變量值(樣本)和在動擾力作用下這些樣本產(chǎn)生的動位移(響應(yīng)值)之間存在某種隱函數(shù)關(guān)系.基于這些樣本和它們的響應(yīng)值,Kriging可以計(jì)算出一種替代的多項(xiàng)式函數(shù)模型.假設(shè)樣本組為x=(x1x2…xn),對應(yīng)的響應(yīng)值為y,Kriging模型可構(gòu)建為[9]

    (1)

    該模型包含線性回歸部分和隨機(jī)部分.β為回歸系數(shù);f(x)為確定性漂移,通常使用x的多項(xiàng)式形式對整個設(shè)計(jì)域進(jìn)行全局近似;隨機(jī)函數(shù)z(x)為波動,對整個設(shè)計(jì)域進(jìn)行局部近似,其均值為0,方差可表示為σ2.

    工程領(lǐng)域常使用稱為DACE的開源Matlab優(yōu)化工具箱來建立樣本與響應(yīng)的Kriging替代函數(shù)模型.使用者僅需提供樣本及響應(yīng)值數(shù)據(jù)和一個用來定義模型建立精確度的θ函數(shù)值,創(chuàng)建最基本的Kriging代理模型[10].

    1.6 期望提高準(zhǔn)則

    積分后可得

    (2)

    其中Φ和φ分別代表正態(tài)分布函數(shù)和正態(tài)密度函數(shù).而優(yōu)化解的搜索方向也即為式(2)的最大值方向[11].

    EI準(zhǔn)則也存在一種加權(quán)形式的表述方式,記為

    E[I(x)]=σ(x)[λvΦ(v)+(1-λ)φ(v)];λ∈[0,1]

    (3)

    權(quán)系數(shù)λ代表著搜索時對局部和全局的重視程度,λ值越小,EI準(zhǔn)則的優(yōu)化解越集中在全局的解空間內(nèi),反之則將搜索過程集中在當(dāng)前的局部解空間內(nèi).

    2 汽輪機(jī)基礎(chǔ)的并行多目標(biāo)優(yōu)化方法

    2.1 汽輪機(jī)基礎(chǔ)的多目標(biāo)優(yōu)化模型

    考慮決定基礎(chǔ)造價成本的基礎(chǔ)質(zhì)量以及決定基礎(chǔ)力學(xué)性能指標(biāo)之一的動位移兩個目標(biāo),汽輪機(jī)基礎(chǔ)的多目標(biāo)優(yōu)化模型可創(chuàng)建為

    w0和u0為使該問題量綱一化所引入的基礎(chǔ)參考質(zhì)量和參考動位移,計(jì)算時采用初始設(shè)計(jì)對應(yīng)的值進(jìn)行賦值.進(jìn)一步采用K-S函數(shù)對優(yōu)化模型進(jìn)行處理后可得

    (5)

    2.2 基于熵的期望提高(EEI)準(zhǔn)則

    為獲取一個最優(yōu)的權(quán)系數(shù)求解式(3)的問題,特引入信息熵模型:

    (6)

    使用拉格朗日乘子法對上述模型轉(zhuǎn)化,可得

    (7)

    η和μ分別為多目標(biāo)優(yōu)化的加權(quán)系數(shù)及拉格朗日乘子,此處η=0.5.在最優(yōu)解處對式(7)求導(dǎo),有

    (8)

    可計(jì)算出

    (9)

    其中r=(η-1)/η,稱為準(zhǔn)權(quán)系數(shù),而式(9)可計(jì)算出最優(yōu)的加權(quán)系數(shù),并進(jìn)而根據(jù)這組加權(quán)系數(shù)預(yù)測出一個最優(yōu)的Pareto解.

    2.3 收斂準(zhǔn)則

    同時考慮EI準(zhǔn)則的精確程度對優(yōu)化的影響以及多目標(biāo)問題目標(biāo)函數(shù)的精確程度,特引入兩個收斂準(zhǔn)則:

    (10)

    (11)

    式(10)表示預(yù)測的最優(yōu)值的提高程度已經(jīng)非常小,因此沒有進(jìn)一步迭代的必要;式(11)表示當(dāng)前模型下Kriging的近似程度已經(jīng)接近真實(shí)情況,故可以停止優(yōu)化.Fn為第n次迭代時的實(shí)際目標(biāo)函數(shù)值.ε1和ε2是收斂系數(shù),均定義為0.05.

    2.4 并行優(yōu)化方法

    本文使用一個由22個計(jì)算節(jié)點(diǎn)和一個控制節(jié)點(diǎn)組成的天梭10000型計(jì)算集群作為并行程序的測試環(huán)境.每個計(jì)算節(jié)點(diǎn)有2個CPU,每個CPU含4個物理核心.集群采用英特爾MPI軟件配置了并行計(jì)算環(huán)境.假定計(jì)算節(jié)點(diǎn)上可利用的物理核心共m個(本文將每個核心對應(yīng)一個并行MPI進(jìn)程),則并行優(yōu)化的方法可描述為

    (1)在控制節(jié)點(diǎn)上使用拉丁超立方采樣方法生成一個N(本文中等于截面向量維度的5倍)個的樣本群.

    (2)使用ANSYS軟件作為黑箱分析模塊對樣本進(jìn)行分析,前m-1個MPI進(jìn)程獲得個樣本,最后一個MPI進(jìn)程獲得剩余的樣本.

    (3)在控制節(jié)點(diǎn)上對所有的樣本分析數(shù)據(jù)進(jìn)行匯總,按式(5)求出所有樣本的目標(biāo)函數(shù)值.

    (4)在控制節(jié)點(diǎn)上對樣本和目標(biāo)函數(shù)建立Kriging代理模型,并根據(jù)式(9)計(jì)算出一個最優(yōu)加權(quán)系數(shù).

    (5)在計(jì)算節(jié)點(diǎn)上,前m-1個MPI進(jìn)程使用加權(quán)系數(shù)λ=j/m計(jì)算式(3)(j為進(jìn)程的編號),得出m-1個優(yōu)化解;最后一個MPI進(jìn)程使用最優(yōu)加權(quán)系數(shù),得到一個優(yōu)化解.

    (6)使用2.3節(jié)的收斂準(zhǔn)則對優(yōu)化過程進(jìn)行控制,同時滿足時優(yōu)化結(jié)束,否則置N=m,轉(zhuǎn)第(2)個優(yōu)化過程.

    3 測試算例

    本文采用兩個實(shí)際工程算例對上述并行多目標(biāo)優(yōu)化方法進(jìn)行測試.分別使用EI準(zhǔn)則及EEI準(zhǔn)則對算例進(jìn)行計(jì)算,并定義計(jì)算節(jié)點(diǎn)可利用的物理核心數(shù)為1、2、4、8、16.

    3.1 測試算例1

    第1個測試算例為一個600MW的汽輪機(jī)基礎(chǔ)模型,包含15個截面設(shè)計(jì)變量,其FEM[12]和優(yōu)化結(jié)果分別如圖1和表1所示.

    圖1 600 MW汽輪機(jī)基礎(chǔ)FEMFig.1 FEM of 600 MW turbine foundation

    3.2 測試算例2

    第2個測試算例為一個300 MW的汽輪機(jī)基礎(chǔ)模型,包含17個截面設(shè)計(jì)變量,其FEM和優(yōu)化結(jié)果分別如圖2和表2所示.

    3.3 算例的加速比與并行效率

    上述算例的并行計(jì)算性能使用表3的數(shù)據(jù)加以描述.

    3.4 結(jié)果分析與討論

    和初始設(shè)計(jì)相比,優(yōu)化后(EI和EEI準(zhǔn)則)的基礎(chǔ)設(shè)計(jì)均有更好的目標(biāo)函數(shù)值(F),對應(yīng)地也會使基礎(chǔ)的自重有所減輕,進(jìn)而減少建造基礎(chǔ)時的經(jīng)濟(jì)投入;在符合式(4)的質(zhì)量和動位移約束條件下,基礎(chǔ)的梁柱截面積部分會有相應(yīng)的減少或者增加,但基礎(chǔ)整體會在迭代過程中朝質(zhì)量的下限和動位移的上限靠近.

    與傳統(tǒng)的EI準(zhǔn)則相比,EEI準(zhǔn)則在對加權(quán)形式的EI準(zhǔn)則進(jìn)一步改進(jìn)后,可以獲得一個更加優(yōu)良的加權(quán)系數(shù),用以對最優(yōu)值進(jìn)行預(yù)測.從表1和表2所示的數(shù)據(jù)來看,使用EEI準(zhǔn)則可以提高最優(yōu)值的精確度,獲取更優(yōu)的目標(biāo)函數(shù)值,進(jìn)而對工程造價等產(chǎn)生積極作用.

    從表1和表2中EEI準(zhǔn)則的優(yōu)化結(jié)果可看出,隨著投入計(jì)算的MPI進(jìn)程的數(shù)量增多,用于預(yù)測最優(yōu)值的EEI加權(quán)系數(shù)也會相應(yīng)增多,在同一個迭代過程中,同時獲得的最優(yōu)值也會增加.樣本群及響應(yīng)值的豐富,使得Kriging替代模型的預(yù)測準(zhǔn)確度增加,進(jìn)而也會增加在最優(yōu)值計(jì)算時的準(zhǔn)確度,提高優(yōu)化結(jié)果的精度.

    以黑箱方法對優(yōu)化進(jìn)行模塊化分割,能夠使得各模塊之間相互獨(dú)立,進(jìn)而使用特定的程序和軟件作最優(yōu)化的分析和計(jì)算.將這種算例使用并行環(huán)境運(yùn)算后,可最大程度地發(fā)揮并行計(jì)算的優(yōu)勢,各進(jìn)程間無頻繁的信息交互,每個進(jìn)程有自己相對匹配的數(shù)據(jù)和結(jié)果.如表3數(shù)據(jù)所示,黑箱形式的優(yōu)化模式可獲得比較高的并行加速比和效率.

    表1 600 MW汽輪機(jī)基礎(chǔ)優(yōu)化結(jié)果Tab.1 Optimized result for 600 MW turbine foundation

    圖2 300 MW汽輪機(jī)基礎(chǔ)FEMFig.2 FEM of 300 MW turbine foundation

    表2 300 MW汽輪機(jī)基礎(chǔ)優(yōu)化結(jié)果Tab.1 Optimized result for 300 MW turbine foundation

    表3 加速比和效率Tab.3 Speed up ratio and efficiency

    4 結(jié) 語

    本文提出的汽輪機(jī)基礎(chǔ)的并行多目標(biāo)優(yōu)化方法,在對傳統(tǒng)的EI準(zhǔn)則進(jìn)行改進(jìn)的基礎(chǔ)上,將熵原理加入加權(quán)形式的EI準(zhǔn)則中(EEI準(zhǔn)則):通過最優(yōu)權(quán)系數(shù)的作用,增加最優(yōu)值的預(yù)測準(zhǔn)度;對使用EEI準(zhǔn)則的優(yōu)化方法進(jìn)行了并行化,在獲得較好的優(yōu)化解的同時,得到了比較滿意的并行加速效果.鑒于黑箱方法在工程優(yōu)化領(lǐng)域的大量采用,本文所發(fā)展的并行多目標(biāo)優(yōu)化方法可進(jìn)而應(yīng)用到其他實(shí)際的工程問題中.

    [1] 馬曉光. 基于Kriging的汽輪機(jī)基礎(chǔ)穩(wěn)健優(yōu)化[D]. 大連:大連理工大學(xué), 2011.

    MA Xiao-guang. Robustness design for turbine engine foundation based on Kriging [D]. Dalian:Dalian University of Technology, 2011. (in Chinese)

    [2]Rao S S. Engineering Optimization:Theory and Practice [M]. 4th ed. New Jersey:John Wiley & Sons, 2009.

    [3]Trobec R, Vajtersic M, Zinterhof P. Parallel Computing:Numerics, Applications, and Trends [M]. London:Springer, 2009.

    [4]Jones D R, Schonlau M, Welch W J. Efficient global optimization of expensive black-box functions [J]. Journal of Global optimization, 1998, 13(4):455-492.

    [5]Deb K. Multi-objective optimization [M] // Burke E K, Kendall G, eds. Search Methodologies. New York:Springer US, 2005:273-316.

    [6]XU Wei. Application of commercial FEA software [M] // GENG Jian-ping, YAN Wei-qi, XU Wei, eds. Application of the Finite Element Method in Implant Dentistry. Berlin:Springer Berlin Heidelberg, 2008:92-134.

    [7]Deb K. Multi-objective optimisation using evolutionary algorithms:an introduction [M] // WANG Li-hui, Ng A H C, Deb K, eds. Multi-objective Evolutionary Optimisation for Product Design and Manufacturing. London:Springer, 2011:3-34.

    [8]Bauke H, Mertens S. Cluster Computing [M]. Berlin:Springer Berlin Heidelberg, 2006.

    [9]Gao Y, Wang X. An effective warpage optimization method in injection molding based on the Kriging model [J]. The International Journal of Advanced Manufacturing Technology, 2008, 37(9-10):953-960.

    [10]Lophaven S N, Nielsen H B, S?ndergaard J. DACE-A Matlab Kriging toolbox, version 2. 0 [R]. Copenhagen:Technical University of Denmark, 2002.

    [11]Kleijnen J P, van Beers W, van Nieuwenhuyse I. Expected improvement in efficient global optimization through bootstrapped Kriging [J]. Journal of Global Optimization, 2012, 54(1):59-73.

    [12]Madenci E, Guven I. The Finite Element Method and Applications in Engineering Using ANSYS [M]. New York:Springer US, 2006.

    A parallel multi-objective optimization method for turbine foundation

    ( 1.School of Computer Science and Technology, Dalian University of Technology, Dalian 116024, China;2.State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China )

    A parallel multi-objective optimization method for turbine foundation is proposed. The problem is described as to find the design variables which get the minimum values of a turbine foundation′s mass and maximum vibration amplitude. The mass of the turbine foundation is a linear expression of the beams and columns′ cross-sectional area, and the vibration amplitude of the foundation acted by additional disturbing force can be calculated by black-box functions. Based on the simulation model which Kriging surrogate function established between the cross-sectional area and the vibration amplitude of the foundation, an entropy-based expected improvement criterion and a parallel strategy are developed to improve the computing accuracy and the efficiency, separately. Two practical turbine foundations are investigated as test cases, and the experimental results demonstrate that the proposed method could be effectively used in this kind of engineering optimization problems with satisfactory optimization solutions achieved.

    turbine foundation; multi-objective optimization; Kriging; parallel computing; expected improvement

    1000-8608(2015)03-0229-07

    2014-10-28;

    2015-03-26.

    “九七三”國家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃資助項(xiàng)目(2012CB025905);高等學(xué)校學(xué)科創(chuàng)新引智計(jì)劃資助項(xiàng)目(B14013);國家自然科學(xué)基金資助項(xiàng)目(11202049);中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(DUT14RC(3)060).

    李召軍(1986-),男,博士生,E-mail:lizhaojun@mail.dlut.edu.cn;李 征*(1982-),男,講師,E-mail:lizheng@dlut.edu.cn.

    TP393;TU323.3

    A

    10.7511/dllgxb201503001

    猜你喜歡
    黑箱算例汽輪機(jī)
    拒絕黑箱
    東汽百萬等級汽輪機(jī)低壓軸承偏載治理研究
    能源工程(2020年5期)2021-01-04 01:29:00
    失蹤的黑箱
    淺析給水泵汽輪機(jī)跳閘回路改造
    廣西電力(2016年4期)2016-07-10 10:23:38
    汽輪機(jī)排汽缸噴水量計(jì)算
    基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
    互補(bǔ)問題算例分析
    基于CYMDIST的配電網(wǎng)運(yùn)行優(yōu)化技術(shù)及算例分析
    燃煤PM10湍流聚并GDE方程算法及算例分析
    135MW汽輪機(jī)組滑壓運(yùn)行經(jīng)濟(jì)性分析
    河南科技(2014年11期)2014-02-27 14:09:55
    大陆偷拍与自拍| 人妻制服诱惑在线中文字幕| 国产精品精品国产色婷婷| 国产av在哪里看| 日韩欧美一区视频在线观看 | 精品一区二区免费观看| 欧美xxⅹ黑人| 免费av观看视频| 99热这里只有精品一区| 啦啦啦啦在线视频资源| 好男人在线观看高清免费视频| 欧美激情在线99| 人妻制服诱惑在线中文字幕| 成人美女网站在线观看视频| 午夜爱爱视频在线播放| 看非洲黑人一级黄片| 免费人成在线观看视频色| 成人午夜高清在线视频| 国内精品美女久久久久久| 又爽又黄无遮挡网站| 99久久中文字幕三级久久日本| 黄色一级大片看看| 乱系列少妇在线播放| 中国国产av一级| 一级毛片 在线播放| 日韩av不卡免费在线播放| 欧美激情国产日韩精品一区| 联通29元200g的流量卡| 高清在线视频一区二区三区| 国产精品熟女久久久久浪| 我的女老师完整版在线观看| 欧美日韩综合久久久久久| 69人妻影院| 高清视频免费观看一区二区 | a级一级毛片免费在线观看| 最近的中文字幕免费完整| 国产精品日韩av在线免费观看| 寂寞人妻少妇视频99o| 熟妇人妻不卡中文字幕| 韩国av在线不卡| 亚洲欧洲日产国产| 美女主播在线视频| 日韩在线高清观看一区二区三区| 2018国产大陆天天弄谢| 精品午夜福利在线看| 1000部很黄的大片| 午夜精品一区二区三区免费看| 97精品久久久久久久久久精品| 欧美性猛交╳xxx乱大交人| 国产精品国产三级国产av玫瑰| 免费看a级黄色片| 少妇高潮的动态图| 成年女人在线观看亚洲视频 | 嫩草影院入口| 欧美一区二区亚洲| 国产国拍精品亚洲av在线观看| 日韩不卡一区二区三区视频在线| 热99在线观看视频| 亚洲欧美精品专区久久| 边亲边吃奶的免费视频| 日本免费在线观看一区| 久久鲁丝午夜福利片| 一级毛片我不卡| 精品人妻视频免费看| 精品国内亚洲2022精品成人| 久久国内精品自在自线图片| 菩萨蛮人人尽说江南好唐韦庄| 在线播放无遮挡| 国产精品久久久久久精品电影小说 | 久久精品国产亚洲av天美| 美女高潮的动态| 免费黄色在线免费观看| 午夜免费激情av| 色播亚洲综合网| 少妇裸体淫交视频免费看高清| 国产精品爽爽va在线观看网站| 日韩av免费高清视频| 精品熟女少妇av免费看| 亚洲av在线观看美女高潮| 天天躁日日操中文字幕| 免费播放大片免费观看视频在线观看| 99久久九九国产精品国产免费| 亚洲国产精品国产精品| 麻豆乱淫一区二区| 精品欧美国产一区二区三| 精品久久久久久久久av| 欧美精品一区二区大全| 特级一级黄色大片| 国产成人精品一,二区| 男人狂女人下面高潮的视频| 亚洲av免费高清在线观看| 天美传媒精品一区二区| 国国产精品蜜臀av免费| 国产精品久久久久久精品电影| 天堂网av新在线| 亚洲av电影不卡..在线观看| 丰满少妇做爰视频| 欧美成人a在线观看| 一区二区三区四区激情视频| 91在线精品国自产拍蜜月| 美女xxoo啪啪120秒动态图| 国产精品女同一区二区软件| 国产黄a三级三级三级人| 22中文网久久字幕| 亚洲国产日韩欧美精品在线观看| 观看免费一级毛片| 免费黄网站久久成人精品| 久久精品国产鲁丝片午夜精品| 青春草国产在线视频| 丰满少妇做爰视频| 女人十人毛片免费观看3o分钟| 少妇裸体淫交视频免费看高清| 欧美性猛交╳xxx乱大交人| 国产亚洲最大av| 亚洲欧洲日产国产| 久久久精品94久久精品| 亚洲国产欧美人成| 亚洲最大成人av| 青青草视频在线视频观看| 乱人视频在线观看| 久久久久久久久久黄片| 女人久久www免费人成看片| 一个人观看的视频www高清免费观看| av免费在线看不卡| 免费播放大片免费观看视频在线观看| 久久久久久久久久成人| 欧美3d第一页| 久久精品久久精品一区二区三区| 色综合亚洲欧美另类图片| 亚洲成人中文字幕在线播放| 国产精品久久久久久久电影| 精品国产一区二区三区久久久樱花 | 国产成人aa在线观看| 国产人妻一区二区三区在| 三级毛片av免费| 精品一区二区三区视频在线| 国产精品嫩草影院av在线观看| 91精品伊人久久大香线蕉| 免费不卡的大黄色大毛片视频在线观看 | 国产精品综合久久久久久久免费| 中国美白少妇内射xxxbb| 老女人水多毛片| 干丝袜人妻中文字幕| 亚洲四区av| 国产精品久久视频播放| 日本免费a在线| 夫妻午夜视频| 国产国拍精品亚洲av在线观看| 免费不卡的大黄色大毛片视频在线观看 | 777米奇影视久久| 欧美精品国产亚洲| 日本wwww免费看| 最新中文字幕久久久久| 97人妻精品一区二区三区麻豆| 亚洲国产精品sss在线观看| 人人妻人人澡欧美一区二区| 亚洲激情五月婷婷啪啪| 亚洲乱码一区二区免费版| 国产三级在线视频| 亚洲成人中文字幕在线播放| 免费av毛片视频| 一级片'在线观看视频| 男女国产视频网站| 免费观看性生交大片5| 一级毛片我不卡| 成人漫画全彩无遮挡| 国产男人的电影天堂91| 国产精品国产三级国产av玫瑰| xxx大片免费视频| 免费大片18禁| 国产精品女同一区二区软件| 禁无遮挡网站| 久久人人爽人人片av| 亚洲av不卡在线观看| 最近的中文字幕免费完整| 欧美丝袜亚洲另类| 久久99热6这里只有精品| 亚洲精品国产av蜜桃| 老师上课跳d突然被开到最大视频| 美女被艹到高潮喷水动态| 亚洲无线观看免费| 欧美不卡视频在线免费观看| 最近中文字幕2019免费版| 国产乱人视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲自偷自拍三级| 青青草视频在线视频观看| 精品国产三级普通话版| 亚洲欧美精品专区久久| 波野结衣二区三区在线| 欧美人与善性xxx| 99视频精品全部免费 在线| 美女大奶头视频| 可以在线观看毛片的网站| 91狼人影院| 国产一级毛片七仙女欲春2| 久久精品国产自在天天线| 综合色av麻豆| 国产精品嫩草影院av在线观看| 蜜桃亚洲精品一区二区三区| 内地一区二区视频在线| 一级毛片久久久久久久久女| 在线观看一区二区三区| 亚洲激情五月婷婷啪啪| 日韩欧美 国产精品| 少妇熟女aⅴ在线视频| 日韩,欧美,国产一区二区三区| 日韩在线高清观看一区二区三区| 国产精品99久久久久久久久| 老司机影院毛片| 黄片wwwwww| 亚洲内射少妇av| 一区二区三区四区激情视频| 亚洲四区av| 男人狂女人下面高潮的视频| 欧美一级a爱片免费观看看| 亚洲欧美日韩东京热| 丝瓜视频免费看黄片| 亚洲av不卡在线观看| 91精品国产九色| 99热这里只有精品一区| 精品国内亚洲2022精品成人| 成年女人在线观看亚洲视频 | 最后的刺客免费高清国语| 丰满少妇做爰视频| 大香蕉97超碰在线| 成人av在线播放网站| 免费在线观看成人毛片| 中文字幕制服av| 国产成人a区在线观看| 小蜜桃在线观看免费完整版高清| 国产伦一二天堂av在线观看| 不卡视频在线观看欧美| 伊人久久精品亚洲午夜| 一级毛片黄色毛片免费观看视频| 亚洲精品亚洲一区二区| 婷婷色综合大香蕉| 久久人人爽人人爽人人片va| 久久久久久久午夜电影| 国产精品女同一区二区软件| 天天躁日日操中文字幕| 亚洲精品国产成人久久av| av在线天堂中文字幕| 婷婷色综合www| 亚洲精品456在线播放app| 国产精品麻豆人妻色哟哟久久 | 黄片wwwwww| 国产亚洲最大av| 国产精品久久视频播放| 少妇熟女欧美另类| 色尼玛亚洲综合影院| 国产精品99久久久久久久久| kizo精华| 天堂av国产一区二区熟女人妻| 国产高清不卡午夜福利| 一级毛片黄色毛片免费观看视频| 亚洲在久久综合| 成人亚洲欧美一区二区av| 欧美人与善性xxx| 欧美人与善性xxx| 一个人看视频在线观看www免费| 赤兔流量卡办理| 国产探花极品一区二区| 亚洲国产日韩欧美精品在线观看| 精品酒店卫生间| 欧美三级亚洲精品| 午夜福利在线观看免费完整高清在| 亚洲真实伦在线观看| 久久久久久久久久黄片| 日韩成人av中文字幕在线观看| 最后的刺客免费高清国语| 欧美日韩亚洲高清精品| 日韩一区二区视频免费看| 亚洲色图av天堂| 亚洲国产精品sss在线观看| 国产成人福利小说| 日韩av在线大香蕉| 欧美zozozo另类| 久久精品夜色国产| or卡值多少钱| 最近2019中文字幕mv第一页| 91狼人影院| 精品国产三级普通话版| 天美传媒精品一区二区| 99久久人妻综合| 亚洲电影在线观看av| 精品久久久久久久久av| 久久人人爽人人片av| 免费av观看视频| av福利片在线观看| 日韩欧美 国产精品| 亚洲欧洲国产日韩| 久久久久久久午夜电影| 日本av手机在线免费观看| 婷婷色综合www| 超碰97精品在线观看| 国产综合懂色| 国产午夜精品论理片| .国产精品久久| 亚洲精品乱码久久久久久按摩| 国产一区二区在线观看日韩| 亚洲av免费在线观看| 国产黄色小视频在线观看| 国产人妻一区二区三区在| 久久精品国产鲁丝片午夜精品| 中文资源天堂在线| 两个人的视频大全免费| 亚洲国产av新网站| 久久99热这里只频精品6学生| 成人性生交大片免费视频hd| 久久6这里有精品| 一本一本综合久久| 伊人久久国产一区二区| 少妇的逼水好多| 国产午夜精品一二区理论片| 日韩,欧美,国产一区二区三区| 欧美人与善性xxx| 自拍偷自拍亚洲精品老妇| 亚洲av成人av| 国语对白做爰xxxⅹ性视频网站| 久久久久久久大尺度免费视频| 欧美+日韩+精品| 只有这里有精品99| 国产欧美另类精品又又久久亚洲欧美| 亚洲高清免费不卡视频| 成人亚洲精品av一区二区| 久久鲁丝午夜福利片| 中文字幕av成人在线电影| 亚洲怡红院男人天堂| 丰满乱子伦码专区| 成年人午夜在线观看视频 | 国产v大片淫在线免费观看| 久久亚洲国产成人精品v| 日本免费在线观看一区| 欧美bdsm另类| 搡老妇女老女人老熟妇| 成年免费大片在线观看| 日韩欧美一区视频在线观看 | 久久亚洲国产成人精品v| 成年女人看的毛片在线观看| 99久久精品热视频| 国产精品久久久久久av不卡| 欧美成人午夜免费资源| 久久久久久久久久久免费av| 亚洲av成人精品一二三区| 欧美另类一区| 校园人妻丝袜中文字幕| 99久国产av精品| 色尼玛亚洲综合影院| 黄色日韩在线| 国产精品综合久久久久久久免费| 国产 亚洲一区二区三区 | 日韩欧美三级三区| 亚洲精品国产av成人精品| 狂野欧美白嫩少妇大欣赏| av在线老鸭窝| 久久久欧美国产精品| 插逼视频在线观看| 亚洲精品色激情综合| av在线播放精品| 91精品国产九色| 麻豆成人av视频| 国产亚洲5aaaaa淫片| 激情 狠狠 欧美| 久久人人爽人人片av| 国产亚洲精品av在线| 国内少妇人妻偷人精品xxx网站| 国产毛片a区久久久久| 91久久精品电影网| 2018国产大陆天天弄谢| 亚洲国产色片| 成人欧美大片| 亚洲精品视频女| 亚洲欧洲日产国产| 国产精品熟女久久久久浪| 亚洲精品456在线播放app| 夜夜爽夜夜爽视频| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 少妇猛男粗大的猛烈进出视频 | 欧美xxxx性猛交bbbb| 大又大粗又爽又黄少妇毛片口| 国产有黄有色有爽视频| 我的女老师完整版在线观看| 午夜福利在线观看吧| 亚洲av在线观看美女高潮| 亚洲成色77777| 久久这里只有精品中国| 91精品伊人久久大香线蕉| 九九爱精品视频在线观看| 三级男女做爰猛烈吃奶摸视频| 国产一区二区在线观看日韩| 久久6这里有精品| 蜜桃久久精品国产亚洲av| 99久久九九国产精品国产免费| 日韩不卡一区二区三区视频在线| 男女视频在线观看网站免费| 亚洲国产精品国产精品| 亚洲欧美一区二区三区黑人 | 亚洲精品国产av成人精品| 中文字幕制服av| 欧美最新免费一区二区三区| 国产不卡一卡二| 亚洲国产欧美在线一区| 久久久久久国产a免费观看| 搡老妇女老女人老熟妇| 男人和女人高潮做爰伦理| 女人被狂操c到高潮| 国产色爽女视频免费观看| av在线亚洲专区| 大话2 男鬼变身卡| 国产av在哪里看| 我的女老师完整版在线观看| 狂野欧美激情性xxxx在线观看| 欧美xxxx黑人xx丫x性爽| 国产亚洲精品久久久com| 老女人水多毛片| 在线a可以看的网站| 亚洲综合色惰| 亚洲综合精品二区| 天堂影院成人在线观看| 大陆偷拍与自拍| 国产欧美另类精品又又久久亚洲欧美| 久久草成人影院| 欧美精品一区二区大全| 色播亚洲综合网| 国产在线男女| 国产精品国产三级国产av玫瑰| 99久国产av精品| 男女啪啪激烈高潮av片| videossex国产| 99re6热这里在线精品视频| 又粗又硬又长又爽又黄的视频| av线在线观看网站| 亚洲欧美日韩无卡精品| 国产黄色视频一区二区在线观看| 国产成人精品福利久久| 噜噜噜噜噜久久久久久91| 免费少妇av软件| 亚洲欧美成人综合另类久久久| 如何舔出高潮| 免费av毛片视频| 久久久a久久爽久久v久久| 久久人人爽人人爽人人片va| 亚洲精华国产精华液的使用体验| 亚洲av电影在线观看一区二区三区 | 国产精品美女特级片免费视频播放器| 人人妻人人澡欧美一区二区| 国产精品久久久久久久久免| 久久久成人免费电影| 久久久a久久爽久久v久久| 午夜激情久久久久久久| 亚洲欧美成人精品一区二区| 亚洲欧美日韩卡通动漫| 在线天堂最新版资源| 精品一区二区三区视频在线| 只有这里有精品99| 亚洲人成网站在线播| 美女国产视频在线观看| 久久久久久久大尺度免费视频| 26uuu在线亚洲综合色| 精品人妻一区二区三区麻豆| 2018国产大陆天天弄谢| 国产午夜精品论理片| 免费观看无遮挡的男女| 欧美日韩一区二区视频在线观看视频在线 | 日韩av免费高清视频| 成人综合一区亚洲| 一级二级三级毛片免费看| 黄色配什么色好看| 最近中文字幕高清免费大全6| 国产精品一区二区在线观看99 | 白带黄色成豆腐渣| 日本欧美国产在线视频| 国产精品国产三级专区第一集| 欧美成人精品欧美一级黄| 亚洲美女视频黄频| 免费少妇av软件| 亚洲精品国产成人久久av| 十八禁网站网址无遮挡 | 精品国产一区二区三区久久久樱花 | 国产av码专区亚洲av| 好男人视频免费观看在线| 国内精品美女久久久久久| 99热全是精品| 成人国产麻豆网| 99热这里只有是精品50| 舔av片在线| 国产成人精品福利久久| 亚洲最大成人手机在线| 男人舔奶头视频| 女的被弄到高潮叫床怎么办| 在线免费观看的www视频| 波多野结衣巨乳人妻| 日本爱情动作片www.在线观看| 久久久久久久久久人人人人人人| 国产淫片久久久久久久久| 亚洲精品,欧美精品| 亚洲精品国产av蜜桃| 色吧在线观看| 免费av观看视频| 色播亚洲综合网| 日本黄大片高清| 青春草国产在线视频| 久久精品夜夜夜夜夜久久蜜豆| 国产成人a区在线观看| 搡老乐熟女国产| 人妻制服诱惑在线中文字幕| 毛片女人毛片| 日韩av不卡免费在线播放| 97精品久久久久久久久久精品| av网站免费在线观看视频 | 国产成人精品一,二区| 欧美日韩精品成人综合77777| 18禁动态无遮挡网站| 亚洲乱码一区二区免费版| 色视频www国产| ponron亚洲| 人妻系列 视频| 国产精品av视频在线免费观看| 免费观看精品视频网站| 日本爱情动作片www.在线观看| 在线天堂最新版资源| 亚洲精品第二区| 日韩不卡一区二区三区视频在线| 国产淫语在线视频| 69人妻影院| 国产精品.久久久| 女人十人毛片免费观看3o分钟| 亚洲欧美日韩东京热| 成年版毛片免费区| 青春草亚洲视频在线观看| 午夜免费男女啪啪视频观看| 激情 狠狠 欧美| 在线播放无遮挡| 免费看美女性在线毛片视频| 中国美白少妇内射xxxbb| or卡值多少钱| 91精品一卡2卡3卡4卡| 18禁在线播放成人免费| 搞女人的毛片| 亚洲欧美日韩卡通动漫| 欧美高清成人免费视频www| 亚洲精品国产成人久久av| 国产探花在线观看一区二区| av又黄又爽大尺度在线免费看| 狂野欧美白嫩少妇大欣赏| 中文字幕av成人在线电影| 三级毛片av免费| 国产高潮美女av| 日韩欧美精品免费久久| 人妻制服诱惑在线中文字幕| 一级毛片电影观看| 丝袜喷水一区| 爱豆传媒免费全集在线观看| 人妻制服诱惑在线中文字幕| 日本av手机在线免费观看| 国产精品一区二区三区四区免费观看| 日本免费a在线| 国产探花极品一区二区| 欧美激情在线99| 国产高清不卡午夜福利| 春色校园在线视频观看| 全区人妻精品视频| 婷婷色综合www| 十八禁网站网址无遮挡 | 亚洲国产精品专区欧美| 国产在视频线在精品| 在线免费观看的www视频| 欧美高清性xxxxhd video| 久久这里只有精品中国| 国产在视频线在精品| 亚洲av免费在线观看| 成年女人在线观看亚洲视频 | 偷拍熟女少妇极品色| 天天一区二区日本电影三级| 卡戴珊不雅视频在线播放| 黄色日韩在线| 韩国av在线不卡| 日韩大片免费观看网站| 熟妇人妻久久中文字幕3abv| 亚洲精品日本国产第一区| 国产不卡一卡二| a级一级毛片免费在线观看| 少妇猛男粗大的猛烈进出视频 | 久久精品综合一区二区三区| 老司机影院毛片| 天堂俺去俺来也www色官网 | 亚洲人成网站在线播| 久久精品综合一区二区三区| 大陆偷拍与自拍| 女人被狂操c到高潮| 天堂√8在线中文| 国产白丝娇喘喷水9色精品| 久久人人爽人人片av| 丰满人妻一区二区三区视频av| 精品少妇黑人巨大在线播放| 久久久色成人| 久久精品夜夜夜夜夜久久蜜豆| 午夜免费男女啪啪视频观看| 午夜免费激情av| 日本午夜av视频| 可以在线观看毛片的网站| 精品久久久久久久末码| 蜜臀久久99精品久久宅男| 亚洲在线自拍视频| av在线播放精品| 男女边摸边吃奶| 免费黄频网站在线观看国产| 日本av手机在线免费观看| 午夜福利高清视频| 亚洲天堂国产精品一区在线| 超碰av人人做人人爽久久| 看十八女毛片水多多多| 国产精品一二三区在线看|