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

    基于PB-NSGA-Ⅲ算法的高速離心泵葉輪性能優(yōu)化研究*

    2023-12-20 12:39:10賴喜德陳小明劉升波
    機(jī)電工程 2023年12期
    關(guān)鍵詞:性能參數(shù)揚(yáng)程葉輪

    張 翔,賴喜德*,陳小明,劉升波

    (1.西華大學(xué) 能源與動力工程學(xué)院,四川 成都 610039;2.成都成發(fā)科能動力工程有限公司,四川 成都 610052)

    0 引 言

    隨著高速離心泵在航空航天、石油化工等行業(yè)應(yīng)用范圍的逐步擴(kuò)大[1],人們對高速泵的性能也提出了更高的要求;但國內(nèi)高速泵的設(shè)計研發(fā)起步較晚[2],現(xiàn)有產(chǎn)品的水力性能達(dá)不到使用要求。

    為提高高速泵的水力性能,國內(nèi)外學(xué)者針對高速泵優(yōu)化設(shè)計開展了大量的研究,并已取得了一定的成果。

    袁軍婭等人[3]建立了高速泵的性能預(yù)測模型,推導(dǎo)出了高速泵效率的數(shù)學(xué)表達(dá)式,以效率為優(yōu)化目標(biāo),對高速泵進(jìn)行了優(yōu)化設(shè)計,優(yōu)化后高速泵的水力效率提高了9%。

    相較于理論計算,數(shù)值模擬計算可以更加精確地預(yù)測性能特性。由此,數(shù)值模擬結(jié)合正交試驗在高速泵的優(yōu)化設(shè)計中逐步增多。

    賀青等人[4]以空化余量為目標(biāo),采用正交優(yōu)化方法得出優(yōu)化方案,對航空燃油離心泵進(jìn)行了抗汽蝕優(yōu)化設(shè)計,優(yōu)化后其蒸汽質(zhì)量分?jǐn)?shù)降低了19.98%;但其空化余量下降不明顯。陳建華等人[5]以揚(yáng)程和水力效率為目標(biāo),采用正交優(yōu)化方法,對高速井泵進(jìn)行了優(yōu)化設(shè)計,優(yōu)化后高速井泵的揚(yáng)程上升了36 m;但井泵的水力效率并未得到提高。

    正交試驗得出的優(yōu)化方案并不是在變量范圍內(nèi)對數(shù)值進(jìn)行全面搜索的結(jié)果,且優(yōu)化后,部分性能目標(biāo)的提升也不明顯。

    嚴(yán)俊峰等人[6]將葉輪幾何參數(shù)作為優(yōu)化變量,僅以效率為目標(biāo),采用遺傳算法在變量的約束范圍內(nèi)全面搜索最優(yōu)參數(shù)組合,為降低優(yōu)化難度而構(gòu)建出優(yōu)化變量間的映射關(guān)系,對低比轉(zhuǎn)速高速泵進(jìn)行了優(yōu)化設(shè)計;然而其在優(yōu)化中需要反復(fù)修正映射的系數(shù),這加大了優(yōu)化設(shè)計的工作量。為克服繁瑣的系數(shù)修正過程,XU Z K等人[7]結(jié)合近似模型與多島遺傳算法,以水力效率和揚(yáng)程為優(yōu)化目標(biāo),對高速電磁泵進(jìn)行了性能優(yōu)化,經(jīng)過104次迭代得出了最優(yōu)參數(shù)組合,優(yōu)化后電磁泵的水力效率上升了6.23%;但該方法迭代次數(shù)多,計算時間長。為減少計算時長,白永明等人[8]基于NSGA-Ⅱ算法對高速泵進(jìn)行了多目標(biāo)優(yōu)化設(shè)計,消除了揚(yáng)程駝峰;但優(yōu)化結(jié)果中可供設(shè)計人員選擇的方案過于單一。

    因此,結(jié)合高速泵特點,有效縮減優(yōu)化變量的數(shù)量以提高效率,減少迭代次數(shù),探索優(yōu)化變量和優(yōu)化目標(biāo)間的聯(lián)系,使優(yōu)化后的性能參數(shù)組合在較大范圍內(nèi)均勻分布,是掌握高速泵多目標(biāo)優(yōu)化設(shè)計的關(guān)鍵。

    近年來,研究人員采用傅里葉幅度靈敏度檢驗擴(kuò)展法(extended Fourier amplitude sensitivity test, EFAST)確定影響高速泵性能的主要幾何參數(shù)[9],并將其作為優(yōu)化變量,以有效縮減變量個數(shù);改進(jìn)的智能算法在葉輪機(jī)械多目標(biāo)優(yōu)化方面[10-12]應(yīng)用廣泛并取得了不錯的效果,這均為高速泵優(yōu)化設(shè)計提供了參考。

    高速泵的性能參數(shù)對幾何參數(shù)比較敏感[13],使得高速泵優(yōu)化難度加大;且針對多目標(biāo)優(yōu)化問題,提升算法預(yù)測能力也是多目標(biāo)優(yōu)化的一個難點。已有研究并未利用改進(jìn)遺傳算法,并結(jié)合高速泵特點對泵進(jìn)行有效地優(yōu)化設(shè)計;且現(xiàn)有研究中應(yīng)用的算法預(yù)測能力較差,算法結(jié)果可供設(shè)計人員選擇的方案較少,高速泵的綜合性能提升不明顯。

    針對高速泵的特點,首先,筆者應(yīng)用EFAST法有效篩選優(yōu)化變量;然后,采用支持向量機(jī)高效構(gòu)建近似模型,在此基礎(chǔ)上采用PB-NSGA-Ⅲ算法[14]進(jìn)行尋優(yōu);最終,形成一套針對高速泵的高效優(yōu)化方法,以有效提升優(yōu)化設(shè)計的能力。

    1 高速泵葉輪多目標(biāo)優(yōu)化設(shè)計流程

    某工程擬采用IS-60-50-158型高速泵。該泵主要由吸入管、誘導(dǎo)輪、葉輪和蝸殼等過流部件組成。但該泵的空化性能、水力效率以及揚(yáng)程不能滿足使用需求。

    為綜合提高該高速泵水力性能,在其他過流部件基本不變的條件下,可以通過優(yōu)化葉輪來達(dá)到優(yōu)化目標(biāo)。筆者以IS-60-50-158型高速泵葉輪的幾何模型為基礎(chǔ),將其作為初始模型進(jìn)行多目標(biāo)優(yōu)化設(shè)計。

    其中,葉輪的優(yōu)化流程如圖1所示。

    圖1 高速泵葉輪的優(yōu)化流程Fig.1 Optimization process of high-speed pump impeller

    在圖1的優(yōu)化流程中,空化余量是空化性能的重要指標(biāo)之一,空化余量需按現(xiàn)有誘導(dǎo)輪揚(yáng)程考慮。

    筆者以水力效率η、空化余量NPSHr和揚(yáng)程H為優(yōu)化目標(biāo),將泵性能參數(shù)組合為{η,NPSHr,H}。

    優(yōu)化流程分為兩步:

    1)建立高速泵葉輪幾何參數(shù)與性能參數(shù)間的近似模型。其中,結(jié)合高速泵的結(jié)構(gòu)特點以及設(shè)計要求,確定幾何參數(shù)的取值范圍;為提高優(yōu)化效率,減少優(yōu)化變量個數(shù),通過靈敏度分析篩選出影響高速泵性能的主要幾何參數(shù),將其作為優(yōu)化變量;利用拉丁超立方試驗在樣本空間中生成優(yōu)化變量的60組幾何參數(shù)組合;在CFX軟件中完成計算流體動力學(xué)(computational fluid dynamics,CFD)數(shù)值計算,獲取60組幾何參數(shù)組合對應(yīng)的性能參數(shù)組合,采用支持向量機(jī)(support vector machine,SVM)建立高速泵幾何參數(shù)組合與性能參數(shù)組合間的近似函數(shù),完成近似模型的構(gòu)建。

    2)利用所建立的近似函數(shù),以PB-NSGA-Ⅲ算法為尋優(yōu)基礎(chǔ),找出性能最優(yōu)的參數(shù)組合,完成多目標(biāo)優(yōu)化設(shè)計。

    2 基于數(shù)值模擬的高速泵性能預(yù)測計算方法

    筆者簡述了初始的高速泵模型,通過數(shù)值計算獲取初始高速泵的性能參數(shù)組合,同時確定數(shù)值計算方案,為后續(xù)的數(shù)值計算提供參考。

    2.1 初始的高速泵模型

    筆者的研究對象是IS-60-50-158型高速泵。該泵為單級懸臂式結(jié)構(gòu),轉(zhuǎn)速n=8 000 r/min,設(shè)計工況Qd=50 m3/h。

    初始高速泵幾何參數(shù)如表1所示。

    表1 初始高速泵的幾何參數(shù)

    2.2 數(shù)值模擬方法

    根據(jù)IS-60-50-158型高速泵的幾何參數(shù),筆者運用CFturbo軟件建立了高速泵幾何模型并導(dǎo)出流體域。

    高速泵的流體域如圖2所示。

    圖2 高速泵流體域Fig.2 High-speed pump fluid domain

    圖2中,高速泵流體域包括葉輪進(jìn)出口段、葉輪和蝸殼。

    為減少數(shù)值計算誤差,筆者在UG/NX中修改了流體域的幾何模型,進(jìn)口段的長度為葉輪進(jìn)口直徑的4倍。筆者采用MESH軟件將高速泵的流體域劃分為非結(jié)構(gòu)網(wǎng)格,交界面和葉片頭部網(wǎng)格需要做加密處理,網(wǎng)格質(zhì)量達(dá)到0.2以上;進(jìn)口采用靜壓,出口采用質(zhì)量流量;采用CFX進(jìn)行數(shù)值模擬并分析其內(nèi)部流場。

    為了解初始模型水力性能,筆者通過數(shù)值模擬獲取了初始模型的性能參數(shù),如表2所示。

    表2 初始模型的性能參數(shù)

    2.3 高速泵的優(yōu)化設(shè)計目標(biāo)參數(shù)

    結(jié)合項目要求可知,高速泵的水力性能不能滿足使用需求。為此,筆者提出的設(shè)計參數(shù)如表3所示。

    表3 高速泵的設(shè)計參數(shù)

    對比表2和表3可以看到:初始的高速泵模型不能滿足設(shè)計要求,需探討其優(yōu)化方法。

    3 幾何參數(shù)與性能參數(shù)間的近似模型

    筆者以提高高速泵水力性能為目標(biāo),應(yīng)用PB-NSGA-Ⅲ算法尋優(yōu)時,需要對目標(biāo)中的性能特性值進(jìn)行預(yù)測;但很難準(zhǔn)確預(yù)測性能特性值。

    因此,筆者利用近似模型建立高速泵葉輪性能參數(shù)組合與幾何參數(shù)組合的一一對應(yīng)關(guān)系,預(yù)測幾何參數(shù)對應(yīng)的性能參數(shù)。近似模型的本質(zhì)是根據(jù)離散幾何參數(shù)樣本和性能參數(shù)構(gòu)建近似函數(shù)。

    筆者建立近似模型包括優(yōu)化變量的選取、拉丁超立方試驗,以及利用支持向量機(jī)建立的幾何參數(shù)與性能參數(shù)間的近似函數(shù)。

    3.1 優(yōu)化變量的選取

    為提高優(yōu)化效率,減少優(yōu)化變量個數(shù),筆者結(jié)合幾何參數(shù)的取值范圍,通過靈敏度分析,篩選出影響高速泵性能的主要幾何參數(shù),并將其作為優(yōu)化變量。

    由于幾何參數(shù)的約束范圍較寬,而全局敏感性分析中的EFAST法在約束范圍較寬的情形中具有優(yōu)越性,其能夠準(zhǔn)確分析幾何參數(shù)的敏感性。因此,筆者采用EFAST法對幾何參數(shù)進(jìn)行全局靈敏度分析。

    首先,需要明確幾何參數(shù)的約束條件。根據(jù)設(shè)計參數(shù)中的比轉(zhuǎn)速,結(jié)合速度系數(shù)法[16]的統(tǒng)計,可選取葉輪進(jìn)口直徑Dj、葉輪直徑D2、葉輪出口寬度b2的取值范圍。袁壽其[17]在大量試驗和理論推導(dǎo)的基礎(chǔ)上,給出了葉片進(jìn)口角β1、葉片出口角β2的推薦值。

    筆者適度放寬取值范圍,得到約束條件,如表4所示。

    表4 幾何參數(shù)的約束條件

    根據(jù)水力損失模型[18],計算水力損失,可得揚(yáng)程和水力效率,可計算泵空化余量[19-21],公式如下:

    H=HT-Δh

    (1)

    (2)

    (3)

    式中:H為揚(yáng)程;h為水力效率;HT為理論揚(yáng)程;Δh為水力損失之和(水力損失包括吸入式水力損失、葉輪進(jìn)口水力損失、葉輪流道摩擦損失、進(jìn)口液流變向水力損失、葉輪出水口水力損失、蝸殼流道摩擦損失、蝸殼內(nèi)擴(kuò)散損失等);l為葉片進(jìn)口稍前的絕對速度;v0為葉片進(jìn)口稍前的相對速度;w0為進(jìn)口相對速度。

    筆者的主要目標(biāo)是提高高速泵的水力效率和空化性能,其次是讓揚(yáng)程達(dá)到設(shè)計要求,因為初始高速泵的揚(yáng)程未達(dá)到設(shè)計要求但相差不大。

    為找到影響高速泵性能的主因,全局靈敏度分析主要針對水力效率與空化余量。筆者結(jié)合幾何參數(shù)約束條件,根據(jù)水力損失模型及泵空化余量計算模型,采用EFAST法進(jìn)行了靈敏度分析,得到一階和全局敏感性指數(shù),如圖3所示。

    圖3 幾何參數(shù)的敏感性指數(shù)Fig.3 Sensitivity index of geometric parameters

    圖3中,從幾何參數(shù)的敏感性指數(shù)中可以篩選出主要幾何參數(shù),將其作為優(yōu)化變量。

    一階敏感性指數(shù)代表單個變量對性能參數(shù)的敏感性;全局敏感性指數(shù)表示變量間的耦合作用。

    在葉輪的幾何參數(shù)中,對空化余量影響最大的是葉輪進(jìn)口直徑Dj,對水力效率影響最大的是葉輪出口寬度b2,其余3個幾何參數(shù)中,對水力效率以及空化余量影響最大的是葉片出口角β2。

    因此,筆者選擇進(jìn)口直徑Dj、葉輪出口寬度b2以及葉片出口角β2作為優(yōu)化變量,記作幾何參數(shù)組合為{Dj,b2,β2}。

    3.2 優(yōu)化變量均勻樣本的生成

    為提高數(shù)值模擬效率,優(yōu)化變量的幾何參數(shù)組合需要均勻分布,從而在優(yōu)化變量的約束范圍內(nèi),使建立的近似模型能夠預(yù)測對應(yīng)的性能參數(shù)值。

    在拉丁超立方抽樣框架下,筆者采用連續(xù)局部枚舉方法和平移傳播算法在抽樣空間內(nèi)生成均勻樣本,生成的60組幾何參數(shù)組合如圖4所示。

    圖4 樣本點分布圖Fig.4 Sample point distribution map

    3.3 幾何參數(shù)與性能參數(shù)間近似函數(shù)的建立

    筆者應(yīng)用PB-NSGA-Ⅲ算法對高速泵葉輪進(jìn)行多目標(biāo)優(yōu)化設(shè)計,迭代中需要確定各性能參數(shù)與幾何參數(shù)組合{Dj,b2,β2}間的函數(shù)表達(dá)式。但對于高速泵的各性能參數(shù),實際上很難使用數(shù)學(xué)表達(dá)式來描述。

    在CFturbo中,筆者利用60組幾何參數(shù)組合生成幾何模型,采用數(shù)值計算獲取60組幾何模型對應(yīng)的性能參數(shù)值組合,構(gòu)建出幾何參數(shù)組合與性能參數(shù)值組合間的一一對應(yīng)關(guān)系,即建立幾何參數(shù)組合與性能參數(shù)值組合間的近似函數(shù),便于尋優(yōu)時利用近似函數(shù)關(guān)系。

    為高效構(gòu)建近似函數(shù),筆者采用林智仁基于支持向量機(jī)開發(fā)的LIBSVM工具箱[22-24],其在構(gòu)建小樣本問題的近似函數(shù)中具有良好的擬合預(yù)測效果。

    支持向量機(jī)的應(yīng)用分為訓(xùn)練集和測試集兩個部分。在60組樣本中,隨機(jī)抽取50個樣本作為訓(xùn)練集,剩余10個作為測試集。支持向量機(jī)構(gòu)造近似函數(shù)的過程較為簡單,在此不加贅述。

    為使近似模型的精度更高,減小模型參數(shù)對模型的影響,同時抑制過擬合和欠擬,筆者應(yīng)用LIBSVM工具箱中的RBF核函數(shù),采用交叉驗證策略選取懲罰因子(c)和方差(g)的最佳參數(shù)組合,利用最佳參數(shù)組合構(gòu)建出訓(xùn)練模型。

    因為近似函數(shù)的精度對尋優(yōu)有較大的影響,所以筆者采用回歸分析對比CFD計算值與支持向量機(jī)預(yù)測值,以評估近似函精度,如圖5所示。

    圖5 支持向量機(jī)測試集的回歸分析Fig.5 Regression analysis of support vector machine test set

    在圖5中,測試集中的性能參數(shù)組合包含支持向量機(jī)所建立的近似函數(shù)的預(yù)測結(jié)果與CFD計算結(jié)果,回歸分析中的支持向量機(jī)預(yù)測值為支持向量機(jī)預(yù)測幾何參數(shù)得出的性能參數(shù)值;回歸分析中的CFD計算值為采用CFX軟件獲取的幾何參數(shù)對應(yīng)的性能參數(shù)值,兩者形成對比,得到?jīng)Q定系數(shù)。

    圖中的直線為45°對角線,其含義為支持向量機(jī)的預(yù)測值與CFD計算的值相等。決定系數(shù)R2越大,說明兩者吻合性越好。決定系數(shù)大于0.9,說明支持向量機(jī)建立的性能參數(shù)與幾何參數(shù)間的近似函數(shù)精度較高,預(yù)測結(jié)果有很高的精度。

    4 高速泵葉輪的多目標(biāo)優(yōu)化

    筆者將構(gòu)造出的近似模型與改進(jìn)的智能算法相結(jié)合。其能夠找出性能最優(yōu)的幾何參數(shù)組合。但改進(jìn)的智能算法在3個目標(biāo)以上的優(yōu)化問題中較難收斂。因此,在高速泵優(yōu)化中,需要采用近似模型預(yù)測幾何參數(shù)對應(yīng)的性能參數(shù),應(yīng)用PB-NSGA-Ⅲ算法找到最優(yōu)性能的參數(shù)組合。

    由于高速泵轉(zhuǎn)速高,其性能對幾何參數(shù)很敏感,子群的適應(yīng)度范圍較寬,從圖5中支持向量機(jī)的預(yù)測值可以看到揚(yáng)程范圍較寬。為此,筆者采取自適應(yīng)策略,避免“早熟”現(xiàn)象。此外,為減少迭代次數(shù),筆者采用PB-NSGA-Ⅲ算法,并和基于參考點的非支配遺傳算法(NSGA-Ⅲ算法)作對比。

    PB-NSGA-Ⅲ與NSGA-Ⅲ算法的流程如圖6所示。

    圖6 PB-NSGA-Ⅲ算法與NSGA-Ⅲ算法流程Fig.6 PB-NSGA-III algorithm and NSGA-III algorithm process

    從圖6中可以看到:PB-NSGA-Ⅲ算法與NSGA-Ⅲ算法不同的是其采用了雙指標(biāo)準(zhǔn)則,在程序中第一個指標(biāo)為個體與帕累托前沿的距離;第二個指標(biāo)為個體到理想點的距離,其中理想點為每代子群中每一個目標(biāo)最小值的組合,筆者取兩個指標(biāo)中的較小值作為選擇操作的依據(jù)[14]。

    5 優(yōu)化結(jié)果及分析

    結(jié)合高速泵的特點,筆者應(yīng)用多學(xué)科優(yōu)化方法以提高其水力性能。為驗證該方法的正確性,筆者利用基于數(shù)學(xué)和計算機(jī)的NSGA-Ⅲ算法和PB-NSGA-Ⅲ算法,找出性能最優(yōu)的參數(shù)組合。參數(shù)組合包含幾何參數(shù)組合與性能參數(shù)組合,對比分析兩種算法結(jié)果,從更優(yōu)秀的算法結(jié)果中選出優(yōu)化后的參數(shù)組合,將其作為優(yōu)化方案。

    筆者結(jié)合數(shù)值模型并運用葉輪機(jī)械相關(guān)理論,分析優(yōu)化方案的參數(shù)組合,確認(rèn)葉輪幾何參數(shù)組合是否合理;并進(jìn)一步對比優(yōu)化前后葉輪的流場,分析優(yōu)化前后流場內(nèi)流線、靜壓及湍動能的變化,驗證優(yōu)化效果。

    5.1 優(yōu)化結(jié)果分析

    筆者應(yīng)用NSGA-Ⅲ算法和PB-NSGA-Ⅲ算法在MATLAB中完成相應(yīng)程序。兩種算法均設(shè)定參考點數(shù)量為200個,種群數(shù)量為200個,迭代次數(shù)為400。

    經(jīng)過程序迭代得到200組參數(shù)組合,200組性能參數(shù)組合{η,NPSHr,η}在空間的分布能夠反映多目標(biāo)優(yōu)化的效果,如圖7所示。

    圖7 帕累托最優(yōu)前沿Fig.7 Pareto optimal front

    從圖7中可以看到:PB-NSGA-Ⅲ算法得到的性能參數(shù)組合集中在狹長窄帶上,分布更為均勻,優(yōu)化結(jié)果中最高水力效率為87.24%,最低效率為82.62%,其中180組性能參數(shù)組合的水力效率大于83.4%;NSGA-Ⅲ算法得到的結(jié)果較為集中,優(yōu)化結(jié)果中最高水力效率為85.24%,最低效率為82.04%,其中114組性能參數(shù)組合的水力效率大于83.4%。

    針對高速泵的優(yōu)化,PB-NSGA-Ⅲ算法的多目標(biāo)優(yōu)化能力更強(qiáng)。水力效率相對較高的解,空化余量也相對較高;但其揚(yáng)程相對較低。

    在PB-NSGA-Ⅲ算法的200組性能參數(shù)組合中,兼顧水力效率與空化余量,同時在揚(yáng)程滿足設(shè)計要求的情況下,筆者篩選出水力效率≥85.5%、空化余量≤9.6 m、揚(yáng)程≥220 m的性能參數(shù)組合,共有17組,然后選擇水力效率最高的一組作為優(yōu)化方案。

    優(yōu)化前后幾何參數(shù)如表5所示。

    表5 初始方案與優(yōu)化方案幾何參數(shù)對比

    由表5可知:不同數(shù)值可反映優(yōu)化前后幾何模型的變化。

    優(yōu)化前后高速泵的水力性能如表6所示。

    表6 初始方案與優(yōu)化方案性能對比

    在表6中,筆者結(jié)合數(shù)值模擬驗證了優(yōu)化方案中的參數(shù)組合。其中,水力效率提高了2.22%,空化余量降低了3.20%,揚(yáng)程提高了6.60%;在設(shè)計工況下,對比支持向量機(jī)的預(yù)測值和數(shù)值計算獲取的結(jié)果,水力效率偏差為0.96%,空化余量偏差為1.06%,揚(yáng)程偏差為1.77%,支持向量機(jī)的預(yù)測值和數(shù)值計算獲取的結(jié)果偏差在2%以內(nèi)。

    以上結(jié)果再次證明近似模型有很高的精度。

    5.2 優(yōu)化前后外特性分析

    優(yōu)化方案中的性能參數(shù)組合是基于數(shù)學(xué)和計算機(jī)知識得出的設(shè)計工況點的水力性能,需運用葉輪機(jī)械理論進(jìn)一步驗證優(yōu)化方案中的幾何參數(shù)組合是否合理。數(shù)值計算可獲取優(yōu)化前后高速泵在各運行工況下的水力效率以及揚(yáng)程。

    優(yōu)化前后水力效率和揚(yáng)程如圖8所示。

    圖8 初始方案與優(yōu)化方案外特性曲線對比Fig.8 Comparison of the external characteristic curves of the initial scheme and the optimization scheme

    從圖8中可以看到:最高水力效率出現(xiàn)在所設(shè)計的工況附近。優(yōu)化前后在0.7Qd~0.9Qd工況的水力效率有幾處相等,但優(yōu)化后在1.0Qd、1.2Qd、1.4Qd工況下,水力效率分別提高了2.22%、5.73%、30.04%;在0.6Qd、0.8Qd、1.0Qd、1.2Qd、1.4Qd工況下,揚(yáng)程分別提高了0.96%、1.97%、6.60%、16.09%、42.83%。

    所以,優(yōu)化后高速泵在大流量工況下水力效率更高,高效率區(qū)更寬;優(yōu)化后高速泵在各個工況點的揚(yáng)程均有提高。

    5.3 高速泵的流場分析

    優(yōu)化前后0.1倍葉高處的流線圖如圖9所示。

    圖9 優(yōu)化前后葉輪中間截面的流線圖Fig.9 The streamline diagram of the middle section of the impeller before and after optimization

    從圖9中可以看到:優(yōu)化前后大多數(shù)區(qū)域的流線方向是穩(wěn)定的,不同的是,優(yōu)化前靠近葉片尾部出現(xiàn)了漩渦,其會導(dǎo)致能量損失,對高速泵的穩(wěn)定性以及水力效率均有一定的影響;而優(yōu)化后漩渦消失,說明優(yōu)化后的葉輪內(nèi)部流態(tài)更為合理,流動損失減小。

    優(yōu)化前后設(shè)計工況點葉片靜壓分布如圖10所示。

    圖10 優(yōu)化前后葉片靜壓圖Fig.10 Blade static pressure diagram before and after optimization

    從圖10中可以看到:優(yōu)化后,葉片背面的靜壓分布更加均勻。在葉片頭部,消除了較小的(-0.514 MPa~-0.256 MPa)的靜壓,葉片頭部低壓區(qū)獲得了一定的控制。這是因為優(yōu)化后葉輪進(jìn)口直徑相對較小,葉輪進(jìn)口、輪轂以及葉片進(jìn)口邊所形成的空間更小,葉輪進(jìn)口到葉片進(jìn)口區(qū)域的壓差增大,葉片頭部的靜壓提高,從而使葉片頭部的低壓區(qū)得到一定控制。

    優(yōu)化前后設(shè)計工況葉輪中間截面湍動能分布如圖11所示。

    圖11 優(yōu)化前后中間截面湍動能分布圖Fig.11 The turbulent kinetic energy distribution map of the middle section before and after optimization

    從圖11中可以看到:湍動能是描述流場中湍流發(fā)展的一個重要物理量,湍動能值越小,則湍流發(fā)散的程度越小,葉片對流體的做功能力越強(qiáng)。從葉片頭部一直到葉片尾部的流道內(nèi),湍動能值較小;但在葉輪與蝸殼交界面附近的區(qū)域湍動能較大,該區(qū)域的湍流程度較大,湍流的脈動長度以及時間長度較大。

    優(yōu)化后,葉輪與蝸殼交界面附近湍流值較大的區(qū)域明顯減小,該區(qū)域的湍流耗散程度減小,說明優(yōu)化后葉輪對流體做功更強(qiáng)。

    6 結(jié)束語

    筆者針對高速泵葉輪進(jìn)行了多目標(biāo)優(yōu)化設(shè)計,優(yōu)化后的高速泵在水力性能上有大幅提高,其可為高速泵的優(yōu)化設(shè)計提供參考。

    研究結(jié)論如下:

    1)筆者基于近似模型和PB-NSGA-Ⅲ算法,對高速泵葉輪進(jìn)行了多目標(biāo)優(yōu)化設(shè)計。優(yōu)化方案中,近似模型的預(yù)測值和數(shù)值計算獲取的結(jié)果偏差在2%以內(nèi),決定系數(shù)R2大于0.9,表明近似模型的預(yù)測精度很高。優(yōu)化后,設(shè)計工況點的水力效率提高了2.22%,揚(yáng)程滿足設(shè)計要求,在大流量工況下水力效率提升明顯,也表明該優(yōu)化方法是可靠的;

    2)優(yōu)化后在0.1倍葉高處,靠近葉片尾部的漩渦消失,使泵在運行過程中的穩(wěn)定性得到了提高。葉片頭部的靜壓有所提高,泵的空化余量下降,空化性能在一定程度上有所提高;

    3)相對于其他優(yōu)化設(shè)計方法,運用PB-NSGA-Ⅲ算法結(jié)合近似模型能夠在變量取值范圍內(nèi)更全面地搜索最優(yōu)解,優(yōu)化結(jié)果中的性能參數(shù)組合在狹長窄帶上分布均勻,且沒有出現(xiàn)大量偏離狹長窄帶的解。因此,采用雙指標(biāo)的PB-NSGA-Ⅲ算法對于高速泵的優(yōu)化設(shè)計是十分有利的。

    筆者研究了高速泵葉輪的優(yōu)化設(shè)計,但是該研究中不包括含誘導(dǎo)輪的高速泵。因此,在后續(xù)研究中,筆者將對含誘導(dǎo)輪的高速泵進(jìn)行優(yōu)化設(shè)計,以便更好地提高高速泵的水力性能。

    猜你喜歡
    性能參數(shù)揚(yáng)程葉輪
    紡織品吸濕速干性能參數(shù)研究
    1.4317 QT2鋼在高能泵葉輪上的應(yīng)用
    管路受力誘發(fā)高揚(yáng)程離心泵振動加劇原因分析
    應(yīng)用石膏型快速精密鑄造技術(shù)制造葉輪
    離心泵葉輪切割方法
    時間數(shù)字轉(zhuǎn)換器性能參數(shù)及誤差分析
    天線電性能參數(shù)的測量
    基于CFD/CSD耦合的葉輪機(jī)葉片失速顫振計算
    新型中高揚(yáng)程大流量水錘泵結(jié)構(gòu)技術(shù)改進(jìn)研究
    中國水利(2015年7期)2015-02-28 15:12:58
    高揚(yáng)程長距離輸水管道系統(tǒng)水錘防護(hù)的模擬分析
    偷拍熟女少妇极品色| 男的添女的下面高潮视频| 亚洲丝袜综合中文字幕| 男人狂女人下面高潮的视频| 午夜免费激情av| 久久久久久久午夜电影| 热99在线观看视频| 亚洲欧美精品自产自拍| 国语自产精品视频在线第100页| 青春草国产在线视频| 欧美一区二区亚洲| 日产精品乱码卡一卡2卡三| 欧美日本视频| 国产一级毛片在线| 联通29元200g的流量卡| 欧美区成人在线视频| 亚洲真实伦在线观看| 亚洲怡红院男人天堂| 美女黄网站色视频| 狠狠狠狠99中文字幕| 男人和女人高潮做爰伦理| 久久草成人影院| 天堂网av新在线| 国产成人一区二区在线| 插逼视频在线观看| 国内揄拍国产精品人妻在线| 18禁在线播放成人免费| 亚洲欧美一区二区三区国产| 成人美女网站在线观看视频| 国产精品1区2区在线观看.| 国产爱豆传媒在线观看| 五月玫瑰六月丁香| 黄色一级大片看看| 欧美激情久久久久久爽电影| 亚洲精品乱码久久久久久按摩| 亚洲成人中文字幕在线播放| 听说在线观看完整版免费高清| 嫩草影院新地址| 麻豆成人午夜福利视频| 亚洲av电影不卡..在线观看| www.av在线官网国产| 大话2 男鬼变身卡| 亚洲国产欧美在线一区| 亚洲高清免费不卡视频| 晚上一个人看的免费电影| 亚洲丝袜综合中文字幕| 老司机影院毛片| 国产一区二区在线观看日韩| 国产亚洲最大av| 一本久久精品| 国产亚洲av片在线观看秒播厂 | 99九九线精品视频在线观看视频| av女优亚洲男人天堂| 日本三级黄在线观看| 日韩在线高清观看一区二区三区| 又粗又爽又猛毛片免费看| 免费不卡的大黄色大毛片视频在线观看 | 一卡2卡三卡四卡精品乱码亚洲| 成人一区二区视频在线观看| 伦理电影大哥的女人| av福利片在线观看| 亚洲成人精品中文字幕电影| 亚洲怡红院男人天堂| 国产黄色小视频在线观看| 能在线免费观看的黄片| 久久久精品大字幕| 日韩欧美 国产精品| 国产精品久久久久久精品电影| 日本爱情动作片www.在线观看| 亚洲欧美一区二区三区国产| 久久这里有精品视频免费| 免费观看在线日韩| 91aial.com中文字幕在线观看| АⅤ资源中文在线天堂| 国产精品嫩草影院av在线观看| 内地一区二区视频在线| 久久鲁丝午夜福利片| 中文天堂在线官网| 亚洲欧洲日产国产| 九草在线视频观看| 99久久成人亚洲精品观看| 国产免费又黄又爽又色| 国产v大片淫在线免费观看| 中文字幕久久专区| 一区二区三区高清视频在线| av国产免费在线观看| 视频中文字幕在线观看| 国产人妻一区二区三区在| 99热这里只有是精品50| 久久久久国产网址| 亚洲欧美精品专区久久| 免费大片18禁| 免费在线观看成人毛片| 亚洲精华国产精华液的使用体验| 中文天堂在线官网| 久久久精品94久久精品| 久久午夜福利片| 国产一级毛片在线| 久久久久久久久久成人| 男人舔女人下体高潮全视频| 搡女人真爽免费视频火全软件| 我要看日韩黄色一级片| 少妇熟女欧美另类| 免费看日本二区| av国产久精品久网站免费入址| 国产亚洲精品久久久com| av国产久精品久网站免费入址| 亚洲一区高清亚洲精品| 国产伦理片在线播放av一区| 亚洲真实伦在线观看| 精品国产一区二区三区久久久樱花 | 国产三级在线视频| 色哟哟·www| 毛片女人毛片| 老女人水多毛片| 日日摸夜夜添夜夜添av毛片| 少妇高潮的动态图| 国产一区二区在线观看日韩| 久久午夜福利片| 中文字幕免费在线视频6| 日韩欧美在线乱码| 1000部很黄的大片| 97超视频在线观看视频| 一区二区三区乱码不卡18| 久久久久九九精品影院| 亚洲av中文字字幕乱码综合| 日本色播在线视频| 深爱激情五月婷婷| 一本一本综合久久| 美女xxoo啪啪120秒动态图| 亚洲婷婷狠狠爱综合网| 精品久久久久久久久久久久久| 色综合色国产| 久久亚洲精品不卡| av在线播放精品| 国产人妻一区二区三区在| 又爽又黄a免费视频| av专区在线播放| 亚洲av日韩在线播放| 国产高清视频在线观看网站| 99在线视频只有这里精品首页| 高清在线视频一区二区三区 | 夜夜看夜夜爽夜夜摸| 日韩av在线免费看完整版不卡| 日韩av不卡免费在线播放| 在线观看一区二区三区| 欧美激情在线99| 超碰97精品在线观看| 一级毛片我不卡| 亚洲美女搞黄在线观看| 在线观看66精品国产| 亚洲中文字幕日韩| 久久精品综合一区二区三区| 久99久视频精品免费| 国产精品女同一区二区软件| 欧美xxxx性猛交bbbb| 亚洲av福利一区| 国产高清国产精品国产三级 | 天堂√8在线中文| 欧美色视频一区免费| 国产综合懂色| 国产免费福利视频在线观看| 国产黄a三级三级三级人| 欧美精品国产亚洲| 非洲黑人性xxxx精品又粗又长| 欧美一级a爱片免费观看看| 亚洲欧美精品自产自拍| 少妇的逼好多水| 亚洲国产精品成人久久小说| 亚洲无线观看免费| 久久国产乱子免费精品| 白带黄色成豆腐渣| 黄片wwwwww| 日日干狠狠操夜夜爽| 晚上一个人看的免费电影| 黄色日韩在线| 高清日韩中文字幕在线| 18禁在线无遮挡免费观看视频| 蜜桃亚洲精品一区二区三区| 九草在线视频观看| 婷婷六月久久综合丁香| 国产亚洲精品久久久com| 老师上课跳d突然被开到最大视频| www.av在线官网国产| 国产高潮美女av| 变态另类丝袜制服| 久久韩国三级中文字幕| 亚洲精品一区蜜桃| 白带黄色成豆腐渣| 国产精品1区2区在线观看.| 国产极品天堂在线| 有码 亚洲区| 久久久久网色| 亚洲在线观看片| 日本-黄色视频高清免费观看| 久久久亚洲精品成人影院| 亚洲av福利一区| 淫秽高清视频在线观看| 男人和女人高潮做爰伦理| 男女那种视频在线观看| 边亲边吃奶的免费视频| 色噜噜av男人的天堂激情| 国产免费视频播放在线视频 | 日本色播在线视频| 69av精品久久久久久| 晚上一个人看的免费电影| 国产一区亚洲一区在线观看| 免费搜索国产男女视频| 国产精品永久免费网站| av在线蜜桃| 精品久久久久久电影网 | 美女xxoo啪啪120秒动态图| 美女高潮的动态| 一个人看视频在线观看www免费| 精品久久久久久久久亚洲| 欧美一区二区亚洲| 久久精品国产亚洲网站| 国产探花极品一区二区| 大香蕉久久网| 久久亚洲国产成人精品v| 精品国产一区二区三区久久久樱花 | 97超视频在线观看视频| 亚洲婷婷狠狠爱综合网| 日本三级黄在线观看| 免费搜索国产男女视频| 成人美女网站在线观看视频| 自拍偷自拍亚洲精品老妇| 日韩精品有码人妻一区| 国产精品嫩草影院av在线观看| 亚洲在线观看片| 午夜福利高清视频| 亚洲av福利一区| 国产片特级美女逼逼视频| 韩国高清视频一区二区三区| 欧美潮喷喷水| 男人舔奶头视频| 色综合站精品国产| 99久国产av精品| 三级经典国产精品| 日本黄大片高清| 久久久久精品久久久久真实原创| 精品久久久噜噜| 国产成人freesex在线| 精品国产三级普通话版| av视频在线观看入口| 晚上一个人看的免费电影| 久久久久国产网址| 中文欧美无线码| 亚洲欧美日韩卡通动漫| 18禁动态无遮挡网站| 在线a可以看的网站| 我的女老师完整版在线观看| 欧美zozozo另类| 久久久久性生活片| 久久久久国产网址| 日韩三级伦理在线观看| 高清av免费在线| 亚洲欧美精品综合久久99| 精品国产一区二区三区久久久樱花 | 乱人视频在线观看| 免费黄色在线免费观看| 久久久久国产网址| 久久久成人免费电影| 麻豆成人午夜福利视频| 亚洲,欧美,日韩| 亚洲怡红院男人天堂| 欧美激情在线99| 老师上课跳d突然被开到最大视频| 成人毛片a级毛片在线播放| 久久人人爽人人爽人人片va| 亚洲精品国产av成人精品| 国产色爽女视频免费观看| 欧美区成人在线视频| 99热网站在线观看| 免费搜索国产男女视频| 最近最新中文字幕大全电影3| 午夜福利网站1000一区二区三区| 免费大片18禁| 久热久热在线精品观看| 国产在视频线精品| 成人午夜高清在线视频| 久热久热在线精品观看| 边亲边吃奶的免费视频| 黄片wwwwww| 2021少妇久久久久久久久久久| 99热精品在线国产| 日本色播在线视频| 一区二区三区免费毛片| 国产欧美日韩精品一区二区| 高清视频免费观看一区二区 | 最近中文字幕高清免费大全6| 晚上一个人看的免费电影| 国产高清不卡午夜福利| av福利片在线观看| 在线观看美女被高潮喷水网站| 在线播放国产精品三级| 亚洲精品亚洲一区二区| 天天躁日日操中文字幕| 看片在线看免费视频| 97超碰精品成人国产| 亚洲色图av天堂| 亚洲成人av在线免费| 美女cb高潮喷水在线观看| 亚洲国产精品久久男人天堂| 永久免费av网站大全| 日本av手机在线免费观看| 少妇熟女欧美另类| 少妇的逼水好多| 日韩中字成人| 亚洲av免费高清在线观看| 91狼人影院| 色5月婷婷丁香| 中文天堂在线官网| 丰满人妻一区二区三区视频av| 少妇熟女欧美另类| 国产成人精品久久久久久| 少妇被粗大猛烈的视频| 成人午夜高清在线视频| 可以在线观看毛片的网站| 少妇人妻一区二区三区视频| 久久久久网色| 看免费成人av毛片| 亚洲怡红院男人天堂| 春色校园在线视频观看| 99久久九九国产精品国产免费| 一边亲一边摸免费视频| 成年av动漫网址| 精品熟女少妇av免费看| 青青草视频在线视频观看| 男女那种视频在线观看| 少妇丰满av| 亚洲最大成人手机在线| 国产精品嫩草影院av在线观看| 黄色配什么色好看| 精品国产三级普通话版| 国产中年淑女户外野战色| 看十八女毛片水多多多| 亚洲av一区综合| 久久久色成人| 国产精品一区二区在线观看99 | 精品久久久久久久久久久久久| 国产亚洲5aaaaa淫片| 一本久久精品| 中文字幕久久专区| 看非洲黑人一级黄片| 亚洲av男天堂| 久久久久久久久久黄片| 91精品国产九色| 国产精品不卡视频一区二区| av免费在线看不卡| 国产一级毛片在线| 中文字幕免费在线视频6| 国产高潮美女av| 日韩av不卡免费在线播放| 日韩人妻高清精品专区| 国产亚洲91精品色在线| 国产精品久久视频播放| 国产精品乱码一区二三区的特点| АⅤ资源中文在线天堂| 日韩欧美国产在线观看| 久久草成人影院| 亚洲中文字幕一区二区三区有码在线看| 免费看美女性在线毛片视频| 青春草国产在线视频| 亚洲美女搞黄在线观看| 少妇被粗大猛烈的视频| 韩国高清视频一区二区三区| 亚洲国产高清在线一区二区三| 亚洲欧美清纯卡通| 内射极品少妇av片p| 亚洲av免费在线观看| 国产一级毛片在线| 2022亚洲国产成人精品| 国产精品熟女久久久久浪| 亚洲国产精品合色在线| 男女边吃奶边做爰视频| 亚洲人成网站在线观看播放| 七月丁香在线播放| 久热久热在线精品观看| 两性午夜刺激爽爽歪歪视频在线观看| 丰满人妻一区二区三区视频av| 晚上一个人看的免费电影| 秋霞伦理黄片| 美女黄网站色视频| 成人性生交大片免费视频hd| 老司机影院毛片| 国产三级中文精品| 精品免费久久久久久久清纯| 中文字幕制服av| 日本黄色视频三级网站网址| 免费大片18禁| 欧美色视频一区免费| 伦理电影大哥的女人| 久久久亚洲精品成人影院| 精品久久国产蜜桃| 美女高潮的动态| 天天躁夜夜躁狠狠久久av| 亚洲精品乱码久久久久久按摩| 波多野结衣高清无吗| 日韩精品青青久久久久久| 久久久精品欧美日韩精品| 18禁在线无遮挡免费观看视频| 乱码一卡2卡4卡精品| 精品人妻偷拍中文字幕| 午夜免费激情av| 亚洲经典国产精华液单| 亚洲国产欧洲综合997久久,| 一卡2卡三卡四卡精品乱码亚洲| 九九在线视频观看精品| av又黄又爽大尺度在线免费看 | 精品午夜福利在线看| 日韩一区二区三区影片| 国产一级毛片七仙女欲春2| 国产精品久久久久久av不卡| 男人的好看免费观看在线视频| 国产老妇女一区| 国产黄色小视频在线观看| 少妇的逼水好多| 久久99蜜桃精品久久| 最近的中文字幕免费完整| 国产69精品久久久久777片| 免费看av在线观看网站| 麻豆成人av视频| 少妇的逼好多水| av黄色大香蕉| 少妇丰满av| 久久草成人影院| 一个人看的www免费观看视频| 我要搜黄色片| videos熟女内射| 日韩中字成人| 久久久久网色| 国产精品.久久久| 久久久精品94久久精品| 成年女人永久免费观看视频| 99热这里只有是精品在线观看| 日韩av在线免费看完整版不卡| 在线免费十八禁| 国语自产精品视频在线第100页| 女人十人毛片免费观看3o分钟| 欧美激情国产日韩精品一区| 九九久久精品国产亚洲av麻豆| 久久草成人影院| 国产精品永久免费网站| 亚洲经典国产精华液单| 国产成人a∨麻豆精品| 国产精品不卡视频一区二区| 国产精品爽爽va在线观看网站| 久久精品熟女亚洲av麻豆精品 | 国产精品女同一区二区软件| 久久精品国产亚洲网站| 成人毛片60女人毛片免费| 日本色播在线视频| 午夜亚洲福利在线播放| 国产大屁股一区二区在线视频| 国产伦精品一区二区三区四那| 日本色播在线视频| 免费看日本二区| 国产精品无大码| 最近2019中文字幕mv第一页| 亚洲国产高清在线一区二区三| 99久久精品国产国产毛片| 精品国产露脸久久av麻豆 | 国产在视频线在精品| 九九久久精品国产亚洲av麻豆| 高清毛片免费看| 久久国产乱子免费精品| 色视频www国产| 国产日韩欧美在线精品| 亚洲欧美清纯卡通| 一级毛片我不卡| 亚洲经典国产精华液单| 91久久精品国产一区二区成人| 成人性生交大片免费视频hd| 日韩一本色道免费dvd| 国产国拍精品亚洲av在线观看| 国产日韩欧美在线精品| 国产高清视频在线观看网站| 夜夜爽夜夜爽视频| 一边亲一边摸免费视频| 免费观看精品视频网站| 亚洲欧洲日产国产| 亚洲精华国产精华液的使用体验| 乱人视频在线观看| 内射极品少妇av片p| 校园人妻丝袜中文字幕| 国产91av在线免费观看| 日韩制服骚丝袜av| 国产免费视频播放在线视频 | 久久久久久大精品| 国产精品不卡视频一区二区| 美女xxoo啪啪120秒动态图| 中文字幕av在线有码专区| 国产午夜福利久久久久久| 综合色av麻豆| 亚洲美女搞黄在线观看| 久久精品国产自在天天线| 桃色一区二区三区在线观看| 日本五十路高清| av国产久精品久网站免费入址| 精品国产一区二区三区久久久樱花 | 亚洲欧美清纯卡通| 国产午夜福利久久久久久| 国产人妻一区二区三区在| 精品人妻偷拍中文字幕| 色视频www国产| 少妇人妻一区二区三区视频| 亚洲av不卡在线观看| 色吧在线观看| 亚洲精华国产精华液的使用体验| 观看美女的网站| 看免费成人av毛片| 一个人观看的视频www高清免费观看| 白带黄色成豆腐渣| 激情 狠狠 欧美| 国产午夜精品论理片| 观看美女的网站| 特大巨黑吊av在线直播| 91午夜精品亚洲一区二区三区| 免费黄色在线免费观看| 亚洲色图av天堂| 一边摸一边抽搐一进一小说| 亚洲第一区二区三区不卡| 天天躁日日操中文字幕| 色哟哟·www| 亚洲成av人片在线播放无| 国产精品人妻久久久影院| 特级一级黄色大片| 一二三四中文在线观看免费高清| 1000部很黄的大片| 少妇人妻精品综合一区二区| 精品一区二区免费观看| 亚洲在线自拍视频| 一夜夜www| 日本黄大片高清| 午夜亚洲福利在线播放| 一个人看视频在线观看www免费| 一区二区三区四区激情视频| 亚洲中文字幕一区二区三区有码在线看| 麻豆成人午夜福利视频| 在线观看一区二区三区| 午夜久久久久精精品| 亚洲av电影在线观看一区二区三区 | 九九久久精品国产亚洲av麻豆| 日韩欧美国产在线观看| 深爱激情五月婷婷| 国产精品一区二区三区四区免费观看| 亚洲国产色片| 汤姆久久久久久久影院中文字幕 | 国国产精品蜜臀av免费| 国产精品麻豆人妻色哟哟久久 | 麻豆国产97在线/欧美| 亚洲av中文字字幕乱码综合| 国内揄拍国产精品人妻在线| 看黄色毛片网站| 观看美女的网站| 1000部很黄的大片| 99热全是精品| 97超视频在线观看视频| 国产探花在线观看一区二区| 少妇猛男粗大的猛烈进出视频 | 老女人水多毛片| 97超视频在线观看视频| 不卡视频在线观看欧美| 国产精品久久视频播放| 2021天堂中文幕一二区在线观| 亚洲国产成人一精品久久久| 久久久久久久午夜电影| 亚洲,欧美,日韩| 国产老妇伦熟女老妇高清| 亚洲精品久久久久久婷婷小说 | 建设人人有责人人尽责人人享有的 | 国产精品乱码一区二三区的特点| 少妇猛男粗大的猛烈进出视频 | 欧美成人午夜免费资源| 久久99热这里只有精品18| 免费黄网站久久成人精品| 日韩三级伦理在线观看| 欧美变态另类bdsm刘玥| 日日撸夜夜添| 熟妇人妻久久中文字幕3abv| 一个人观看的视频www高清免费观看| 日本免费在线观看一区| 女的被弄到高潮叫床怎么办| 国产精品,欧美在线| 国产伦精品一区二区三区四那| 欧美精品国产亚洲| 伊人久久精品亚洲午夜| 精品久久久久久久末码| 国产一级毛片七仙女欲春2| 国产中年淑女户外野战色| 欧美日韩精品成人综合77777| av在线老鸭窝| 国产高清视频在线观看网站| av在线亚洲专区| 亚洲婷婷狠狠爱综合网| 亚洲成色77777| 精品人妻熟女av久视频| 久久99蜜桃精品久久| 欧美另类亚洲清纯唯美| 中国国产av一级| 三级经典国产精品| 春色校园在线视频观看| 久久精品久久久久久噜噜老黄 | 国产精品综合久久久久久久免费| 1024手机看黄色片| 国产精品.久久久| 插阴视频在线观看视频| 精品人妻偷拍中文字幕| 成人性生交大片免费视频hd| 精品人妻视频免费看| 久久久久久久久久久丰满| 日韩欧美三级三区| 日韩,欧美,国产一区二区三区 | 色网站视频免费|