• <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
    国产精品一区二区在线观看99| 老汉色av国产亚洲站长工具| 午夜免费鲁丝| 国产精品.久久久| 亚洲精品日本国产第一区| 国产亚洲欧美精品永久| 国产精品一国产av| 日韩精品免费视频一区二区三区| 美女国产高潮福利片在线看| av又黄又爽大尺度在线免费看| 成人亚洲精品一区在线观看| 亚洲图色成人| 一级毛片我不卡| 自拍欧美九色日韩亚洲蝌蚪91| 日本黄色日本黄色录像| 亚洲国产欧美一区二区综合| 在线亚洲精品国产二区图片欧美| 欧美日韩一级在线毛片| 午夜免费观看性视频| 亚洲成人手机| 亚洲色图 男人天堂 中文字幕| 欧美在线黄色| 午夜免费鲁丝| 丁香六月天网| 18在线观看网站| 91麻豆精品激情在线观看国产 | 狠狠婷婷综合久久久久久88av| 亚洲国产成人一精品久久久| 久久99热这里只频精品6学生| 男人操女人黄网站| 国产免费现黄频在线看| 曰老女人黄片| 首页视频小说图片口味搜索 | 1024视频免费在线观看| tube8黄色片| 亚洲精品av麻豆狂野| 又大又爽又粗| 在线观看一区二区三区激情| 中文字幕av电影在线播放| 国产福利在线免费观看视频| 亚洲 国产 在线| 50天的宝宝边吃奶边哭怎么回事| 久久精品成人免费网站| 国产日韩欧美在线精品| 精品国产国语对白av| 香蕉国产在线看| 国产一区有黄有色的免费视频| 久久综合国产亚洲精品| 免费观看av网站的网址| av欧美777| 在线观看www视频免费| 精品欧美一区二区三区在线| 丝袜人妻中文字幕| 欧美日本中文国产一区发布| 老熟女久久久| 黄色毛片三级朝国网站| 亚洲国产av影院在线观看| a级片在线免费高清观看视频| 97精品久久久久久久久久精品| 久久ye,这里只有精品| 免费人妻精品一区二区三区视频| 啦啦啦 在线观看视频| 国产在线一区二区三区精| 午夜福利影视在线免费观看| 如日韩欧美国产精品一区二区三区| 深夜精品福利| 秋霞在线观看毛片| 欧美日韩亚洲国产一区二区在线观看 | 国产三级黄色录像| 丝袜在线中文字幕| 亚洲精品久久午夜乱码| 一级毛片 在线播放| 亚洲一区二区三区欧美精品| 国产精品九九99| 欧美亚洲日本最大视频资源| 人人澡人人妻人| 在线观看免费日韩欧美大片| 亚洲,欧美,日韩| 久久国产精品男人的天堂亚洲| 亚洲欧洲精品一区二区精品久久久| 精品国产乱码久久久久久男人| 国产在线观看jvid| 日韩一区二区三区影片| 精品国产国语对白av| 欧美日韩一级在线毛片| 免费女性裸体啪啪无遮挡网站| 两个人免费观看高清视频| 欧美乱码精品一区二区三区| svipshipincom国产片| 亚洲三区欧美一区| 久久精品国产亚洲av高清一级| 成年av动漫网址| 在线观看免费高清a一片| 99精品久久久久人妻精品| 黄色一级大片看看| 91字幕亚洲| 丰满少妇做爰视频| 狂野欧美激情性xxxx| 只有这里有精品99| 一区二区三区精品91| 久热这里只有精品99| 欧美 日韩 精品 国产| 在现免费观看毛片| 亚洲,一卡二卡三卡| 十分钟在线观看高清视频www| 国产成人精品久久二区二区91| 国产成人系列免费观看| 热99国产精品久久久久久7| 爱豆传媒免费全集在线观看| 国产精品成人在线| 黄频高清免费视频| 99国产精品免费福利视频| 男女无遮挡免费网站观看| 好男人电影高清在线观看| 亚洲精品一卡2卡三卡4卡5卡 | kizo精华| 精品国产一区二区三区四区第35| 晚上一个人看的免费电影| 免费日韩欧美在线观看| 啦啦啦中文免费视频观看日本| 只有这里有精品99| 亚洲国产日韩一区二区| 日韩欧美一区视频在线观看| 国产爽快片一区二区三区| 亚洲精品日本国产第一区| 久久久国产一区二区| 国产成人影院久久av| 久久狼人影院| 日韩欧美一区视频在线观看| 欧美日韩亚洲高清精品| 国产高清视频在线播放一区 | 欧美激情高清一区二区三区| 丝瓜视频免费看黄片| 中文字幕精品免费在线观看视频| 国产成人av教育| 精品免费久久久久久久清纯 | 99香蕉大伊视频| 性高湖久久久久久久久免费观看| 中文字幕制服av| 亚洲av片天天在线观看| 各种免费的搞黄视频| 亚洲七黄色美女视频| 在线观看免费视频网站a站| 欧美精品亚洲一区二区| 亚洲欧美日韩高清在线视频 | 国产亚洲午夜精品一区二区久久| 美女大奶头黄色视频| 国产一区二区激情短视频 | 夜夜骑夜夜射夜夜干| 校园人妻丝袜中文字幕| 久久久久国产一级毛片高清牌| 久久人妻熟女aⅴ| 亚洲国产看品久久| 国产黄色视频一区二区在线观看| 久久久久精品国产欧美久久久 | 别揉我奶头~嗯~啊~动态视频 | 亚洲欧美清纯卡通| svipshipincom国产片| 极品人妻少妇av视频| 最新的欧美精品一区二区| 精品少妇黑人巨大在线播放| 精品亚洲成a人片在线观看| 曰老女人黄片| 青青草视频在线视频观看| 国产一区亚洲一区在线观看| 日本午夜av视频| 丝袜美腿诱惑在线| 国产亚洲精品第一综合不卡| 女人爽到高潮嗷嗷叫在线视频| 婷婷丁香在线五月| 波多野结衣一区麻豆| 精品国产乱码久久久久久小说| 最新的欧美精品一区二区| 精品人妻熟女毛片av久久网站| 国产高清视频在线播放一区 | 国产亚洲av片在线观看秒播厂| 日韩一区二区三区影片| 一边摸一边抽搐一进一出视频| 99久久综合免费| 啦啦啦中文免费视频观看日本| 欧美中文综合在线视频| 国产成人av教育| 日本午夜av视频| 中国美女看黄片| 国产91精品成人一区二区三区 | 一级毛片我不卡| 亚洲国产中文字幕在线视频| 日韩大码丰满熟妇| 51午夜福利影视在线观看| 国产精品国产av在线观看| 一区二区三区激情视频| 又黄又粗又硬又大视频| www.熟女人妻精品国产| 精品国产一区二区久久| 一级片免费观看大全| 国产精品.久久久| 国产精品香港三级国产av潘金莲 | 2021少妇久久久久久久久久久| 2018国产大陆天天弄谢| 99精品久久久久人妻精品| 国产免费福利视频在线观看| 亚洲成色77777| 亚洲国产日韩一区二区| 亚洲久久久国产精品| 精品国产一区二区三区久久久樱花| 亚洲中文av在线| 黄色视频不卡| 国产在线观看jvid| av线在线观看网站| 老司机靠b影院| 亚洲中文字幕日韩| 97精品久久久久久久久久精品| 日本欧美国产在线视频| 精品第一国产精品| 国产97色在线日韩免费| 视频在线观看一区二区三区| 黑丝袜美女国产一区| 熟女av电影| 日本欧美视频一区| 国产无遮挡羞羞视频在线观看| 老司机深夜福利视频在线观看 | 丝袜美足系列| 亚洲情色 制服丝袜| 久久99精品国语久久久| 免费观看人在逋| 午夜福利,免费看| 久久国产精品男人的天堂亚洲| 在线观看人妻少妇| 自线自在国产av| 久久久久久久大尺度免费视频| 亚洲熟女毛片儿| 丝袜脚勾引网站| 超碰成人久久| 免费在线观看视频国产中文字幕亚洲 | 老司机影院成人| av有码第一页| 精品人妻在线不人妻| 亚洲国产精品国产精品| 97精品久久久久久久久久精品| 精品高清国产在线一区| 欧美日韩亚洲国产一区二区在线观看 | 国产麻豆69| 极品人妻少妇av视频| 亚洲伊人久久精品综合| 亚洲伊人久久精品综合| 日本vs欧美在线观看视频| 国产精品秋霞免费鲁丝片| av又黄又爽大尺度在线免费看| 日韩伦理黄色片| 久久这里只有精品19| av片东京热男人的天堂| 十八禁人妻一区二区| 交换朋友夫妻互换小说| 大香蕉久久网| 看十八女毛片水多多多| 国产淫语在线视频| 精品人妻一区二区三区麻豆| 黄色 视频免费看| 99热全是精品| 高清欧美精品videossex| 免费在线观看日本一区| 一级片'在线观看视频| 亚洲国产av影院在线观看| 日本vs欧美在线观看视频| 18在线观看网站| 日韩一本色道免费dvd| 一区二区av电影网| 99久久综合免费| 国产亚洲精品第一综合不卡| 80岁老熟妇乱子伦牲交| 中文字幕人妻丝袜一区二区| 久久中文字幕一级| 国产精品国产三级专区第一集| 亚洲欧洲国产日韩| 各种免费的搞黄视频| 日韩中文字幕视频在线看片| av网站在线播放免费| 国产成人啪精品午夜网站| 久久ye,这里只有精品| 我的亚洲天堂| av在线播放精品| 91精品国产国语对白视频| 91九色精品人成在线观看| 高清欧美精品videossex| 成年人午夜在线观看视频| 久久精品亚洲熟妇少妇任你| 国产精品久久久久成人av| 午夜激情av网站| 国产精品久久久久久精品电影小说| 午夜免费鲁丝| 欧美中文综合在线视频| 精品欧美一区二区三区在线| 桃花免费在线播放| 亚洲伊人久久精品综合| 亚洲av成人精品一二三区| 久9热在线精品视频| 一区在线观看完整版| 中文欧美无线码| 日本色播在线视频| 99精品久久久久人妻精品| a级毛片在线看网站| 永久免费av网站大全| 亚洲av欧美aⅴ国产| 久久天躁狠狠躁夜夜2o2o | 美女高潮到喷水免费观看| 夫妻性生交免费视频一级片| 国产女主播在线喷水免费视频网站| 99国产精品免费福利视频| av网站免费在线观看视频| 亚洲免费av在线视频| 99久久综合免费| 欧美精品亚洲一区二区| 国产精品国产三级国产专区5o| 一区二区三区乱码不卡18| 日韩免费高清中文字幕av| 免费在线观看视频国产中文字幕亚洲 | 免费av中文字幕在线| 久久午夜综合久久蜜桃| 久久人人爽av亚洲精品天堂| 校园人妻丝袜中文字幕| 国产av国产精品国产| 丝袜美足系列| 一本综合久久免费| 亚洲午夜精品一区,二区,三区| 97人妻天天添夜夜摸| 丝袜美足系列| 亚洲成人手机| 99国产精品一区二区三区| 男人舔女人的私密视频| 午夜激情av网站| 手机成人av网站| 午夜精品国产一区二区电影| 国产淫语在线视频| 99国产精品免费福利视频| 亚洲欧美色中文字幕在线| 丝袜喷水一区| 欧美激情 高清一区二区三区| www.av在线官网国产| 国产成人精品久久久久久| 中国美女看黄片| 操美女的视频在线观看| 2018国产大陆天天弄谢| 亚洲国产中文字幕在线视频| 国产精品一区二区在线不卡| av国产精品久久久久影院| 大片免费播放器 马上看| 免费在线观看视频国产中文字幕亚洲 | 成人黄色视频免费在线看| 色综合亚洲欧美另类图片| 在线观看免费视频日本深夜| 国产三级在线视频| 露出奶头的视频| 中文字幕最新亚洲高清| 亚洲电影在线观看av| 深夜精品福利| 国产av又大| 国产99久久九九免费精品| 久久久久久九九精品二区国产 | 日韩精品免费视频一区二区三区| 老鸭窝网址在线观看| 级片在线观看| 久久久久久免费高清国产稀缺| 国产欧美日韩精品亚洲av| 国产熟女午夜一区二区三区| 亚洲人成网站高清观看| 丝袜在线中文字幕| 国产精品九九99| 亚洲免费av在线视频| 欧美成人午夜精品| 亚洲第一青青草原| 国产熟女午夜一区二区三区| 亚洲av中文字字幕乱码综合 | 午夜福利在线观看吧| www.精华液| 婷婷精品国产亚洲av| 啦啦啦免费观看视频1| 91av网站免费观看| 欧美黄色淫秽网站| 色综合婷婷激情| 亚洲精品在线观看二区| 露出奶头的视频| 国产极品粉嫩免费观看在线| 一级毛片高清免费大全| 国产爱豆传媒在线观看 | 日韩 欧美 亚洲 中文字幕| а√天堂www在线а√下载| 少妇的丰满在线观看| 国产精品亚洲美女久久久| 一a级毛片在线观看| 国产黄色小视频在线观看| 性欧美人与动物交配| 手机成人av网站| 久久狼人影院| 亚洲国产精品成人综合色| 亚洲国产中文字幕在线视频| 侵犯人妻中文字幕一二三四区| 亚洲国产欧洲综合997久久, | 免费在线观看影片大全网站| 一个人免费在线观看的高清视频| 可以免费在线观看a视频的电影网站| 欧美不卡视频在线免费观看 | 校园春色视频在线观看| 女同久久另类99精品国产91| 变态另类丝袜制服| 婷婷丁香在线五月| 国产精品1区2区在线观看.| 日本精品一区二区三区蜜桃| 99热6这里只有精品| 美女国产高潮福利片在线看| 午夜福利欧美成人| 视频在线观看一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 不卡av一区二区三区| 日本 av在线| 久久亚洲真实| 久久久久精品国产欧美久久久| 成人特级黄色片久久久久久久| 久久久久久久久久黄片| 欧美人与性动交α欧美精品济南到| 精品高清国产在线一区| 神马国产精品三级电影在线观看 | 岛国视频午夜一区免费看| 欧美人与性动交α欧美精品济南到| or卡值多少钱| 美女扒开内裤让男人捅视频| 老司机靠b影院| 制服丝袜大香蕉在线| 国产高清videossex| 好男人电影高清在线观看| 性色av乱码一区二区三区2| 久久热在线av| 亚洲精品一区av在线观看| 18美女黄网站色大片免费观看| 怎么达到女性高潮| 国产又爽黄色视频| 香蕉国产在线看| 亚洲专区字幕在线| 亚洲,欧美精品.| 亚洲精品国产一区二区精华液| 嫩草影视91久久| 欧美zozozo另类| 啦啦啦观看免费观看视频高清| 波多野结衣av一区二区av| 丝袜人妻中文字幕| 国产激情欧美一区二区| www.自偷自拍.com| 极品教师在线免费播放| 成熟少妇高潮喷水视频| 亚洲人成网站在线播放欧美日韩| 在线免费观看的www视频| 日韩有码中文字幕| 成人欧美大片| 此物有八面人人有两片| 国产成人精品无人区| 亚洲国产精品久久男人天堂| av中文乱码字幕在线| 日日干狠狠操夜夜爽| 国产亚洲精品av在线| 在线永久观看黄色视频| 在线观看免费午夜福利视频| 日韩精品中文字幕看吧| 哪里可以看免费的av片| 亚洲欧美精品综合久久99| 午夜视频精品福利| 成人午夜高清在线视频 | 美女大奶头视频| 中文字幕人妻熟女乱码| 免费无遮挡裸体视频| 久久久久久国产a免费观看| 亚洲成国产人片在线观看| 在线观看免费视频日本深夜| 色综合站精品国产| 校园春色视频在线观看| 亚洲精品中文字幕在线视频| 国产精品精品国产色婷婷| 久久久久精品国产欧美久久久| 黄色 视频免费看| 国产伦人伦偷精品视频| 制服人妻中文乱码| 久久欧美精品欧美久久欧美| 啪啪无遮挡十八禁网站| 欧美最黄视频在线播放免费| 国产又爽黄色视频| 欧美成狂野欧美在线观看| 欧美 亚洲 国产 日韩一| 午夜a级毛片| 国产国语露脸激情在线看| 在线十欧美十亚洲十日本专区| 亚洲av成人av| 99久久综合精品五月天人人| 妹子高潮喷水视频| 一夜夜www| 91老司机精品| 国产不卡一卡二| 国产欧美日韩一区二区精品| 天天一区二区日本电影三级| 久久天堂一区二区三区四区| 夜夜夜夜夜久久久久| 亚洲 欧美 日韩 在线 免费| 99国产精品一区二区蜜桃av| 久久青草综合色| 少妇被粗大的猛进出69影院| 免费在线观看成人毛片| 色哟哟哟哟哟哟| a级毛片在线看网站| 中文字幕精品免费在线观看视频| 天堂√8在线中文| 伦理电影免费视频| 97人妻精品一区二区三区麻豆 | 精品一区二区三区av网在线观看| 久久久久久人人人人人| 欧美在线一区亚洲| 中文字幕av电影在线播放| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲成人国产一区在线观看| 亚洲欧美激情综合另类| 麻豆久久精品国产亚洲av| 免费女性裸体啪啪无遮挡网站| 精品国内亚洲2022精品成人| 一a级毛片在线观看| 亚洲欧美激情综合另类| 国产亚洲精品久久久久久毛片| 国产伦人伦偷精品视频| 中文字幕久久专区| 国产av一区二区精品久久| 国内精品久久久久久久电影| 黄色成人免费大全| 免费看美女性在线毛片视频| 中文字幕久久专区| 国产精品久久久人人做人人爽| 欧美丝袜亚洲另类 | 久久国产精品男人的天堂亚洲| 免费在线观看亚洲国产| 悠悠久久av| 一本一本综合久久| 亚洲人成网站在线播放欧美日韩| 人成视频在线观看免费观看| 久久精品国产综合久久久| 黄色视频不卡| 欧美一级a爱片免费观看看 | 国产私拍福利视频在线观看| 美女扒开内裤让男人捅视频| 波多野结衣av一区二区av| xxxwww97欧美| 日韩欧美国产一区二区入口| 午夜精品在线福利| 国产精品一区二区精品视频观看| 亚洲无线在线观看| 国内精品久久久久精免费| 日韩欧美一区二区三区在线观看| 少妇的丰满在线观看| 此物有八面人人有两片| 亚洲成av人片免费观看| 黄色成人免费大全| 日日干狠狠操夜夜爽| 国产熟女xx| 久久久久久国产a免费观看| 91老司机精品| 国产av一区二区精品久久| 特大巨黑吊av在线直播 | 制服诱惑二区| 少妇 在线观看| 一个人观看的视频www高清免费观看 | 欧美黄色淫秽网站| 国产精品影院久久| 午夜精品久久久久久毛片777| 亚洲国产高清在线一区二区三 | 国产区一区二久久| 午夜福利在线在线| 99精品久久久久人妻精品| 亚洲电影在线观看av| 嫁个100分男人电影在线观看| 亚洲人成网站在线播放欧美日韩| 在线观看午夜福利视频| 亚洲专区国产一区二区| 久久精品成人免费网站| 婷婷精品国产亚洲av| 最近最新免费中文字幕在线| 国产精品免费视频内射| 老鸭窝网址在线观看| 亚洲专区国产一区二区| 麻豆av在线久日| 国产亚洲av高清不卡| www.www免费av| 久久香蕉精品热| 成人手机av| 免费在线观看影片大全网站| 男人操女人黄网站| 老熟妇仑乱视频hdxx| 国产国语露脸激情在线看| 日韩欧美国产一区二区入口| aaaaa片日本免费| 欧美+亚洲+日韩+国产| 1024手机看黄色片| 精品一区二区三区视频在线观看免费| 18禁裸乳无遮挡免费网站照片 | 一级黄色大片毛片| 欧美zozozo另类| 亚洲午夜精品一区,二区,三区| 精品免费久久久久久久清纯| 日韩国内少妇激情av| 欧美成人一区二区免费高清观看 | 1024视频免费在线观看| 日本黄色视频三级网站网址| 一本久久中文字幕| 亚洲成人免费电影在线观看| 在线观看免费日韩欧美大片| 免费无遮挡裸体视频| 一级片免费观看大全| 国产黄片美女视频| 日韩高清综合在线| 一区二区三区激情视频| 99在线视频只有这里精品首页| 色哟哟哟哟哟哟| 亚洲人成电影免费在线|