• <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中文av极速乱| 看非洲黑人一级黄片| 国产淫语在线视频| 国产黄色视频一区二区在线观看| 国产亚洲最大av| 香蕉丝袜av| a级片在线免费高清观看视频| 免费看av在线观看网站| 在现免费观看毛片| 国产综合精华液| 日韩大片免费观看网站| 免费高清在线观看视频在线观看| 亚洲精华国产精华液的使用体验| 一区二区三区激情视频| 两个人看的免费小视频| 日本色播在线视频| 日韩,欧美,国产一区二区三区| 97在线人人人人妻| 人妻人人澡人人爽人人| 午夜久久久在线观看| 午夜日本视频在线| 99re6热这里在线精品视频| 欧美 日韩 精品 国产| 男女国产视频网站| 最新中文字幕久久久久| 国产成人精品久久久久久| 亚洲三级黄色毛片| 麻豆精品久久久久久蜜桃| 午夜福利视频精品| 亚洲av福利一区| 欧美xxⅹ黑人| 亚洲欧美成人精品一区二区| 女人高潮潮喷娇喘18禁视频| 99久久精品国产国产毛片| 激情五月婷婷亚洲| 波多野结衣一区麻豆| 在线天堂最新版资源| 国产xxxxx性猛交| 在现免费观看毛片| 少妇的逼水好多| 另类亚洲欧美激情| 亚洲精华国产精华液的使用体验| 成人国产麻豆网| 国产一区二区三区综合在线观看| 欧美在线黄色| 999精品在线视频| 一本久久精品| 久久99热这里只频精品6学生| 美女国产视频在线观看| 色94色欧美一区二区| 国产欧美日韩一区二区三区在线| 日本av手机在线免费观看| 亚洲欧洲国产日韩| 精品国产一区二区三区久久久樱花| 亚洲人成网站在线观看播放| 国产成人欧美| 亚洲欧洲精品一区二区精品久久久 | 丝瓜视频免费看黄片| 91精品国产国语对白视频| 欧美另类一区| 国产欧美日韩综合在线一区二区| 久久久国产精品麻豆| 高清在线视频一区二区三区| av不卡在线播放| av网站在线播放免费| 一级,二级,三级黄色视频| 日本欧美视频一区| 亚洲av综合色区一区| av不卡在线播放| 丰满乱子伦码专区| 18禁裸乳无遮挡动漫免费视频| 亚洲激情五月婷婷啪啪| 午夜福利在线观看免费完整高清在| 秋霞伦理黄片| 天天躁夜夜躁狠狠躁躁| 欧美成人午夜精品| 不卡av一区二区三区| 啦啦啦视频在线资源免费观看| 欧美bdsm另类| freevideosex欧美| 中文字幕人妻丝袜制服| 国产精品国产av在线观看| 免费在线观看完整版高清| 免费黄频网站在线观看国产| 国产爽快片一区二区三区| 一区二区三区激情视频| av女优亚洲男人天堂| 午夜免费鲁丝| 欧美中文综合在线视频| 精品国产乱码久久久久久小说| 亚洲熟女精品中文字幕| 欧美精品一区二区免费开放| 欧美精品人与动牲交sv欧美| 日韩,欧美,国产一区二区三区| 国产一区有黄有色的免费视频| 亚洲av综合色区一区| 国产成人aa在线观看| 国产精品久久久久久精品电影小说| 国语对白做爰xxxⅹ性视频网站| 在线观看美女被高潮喷水网站| 黑人欧美特级aaaaaa片| 久久综合国产亚洲精品| 国产一区二区 视频在线| 国产精品欧美亚洲77777| 可以免费在线观看a视频的电影网站 | 老司机影院成人| 在现免费观看毛片| 亚洲欧美一区二区三区久久| 丰满少妇做爰视频| 午夜91福利影院| 春色校园在线视频观看| 亚洲国产欧美网| 人妻人人澡人人爽人人| 一本久久精品| 婷婷色综合www| 欧美日韩国产mv在线观看视频| 日日爽夜夜爽网站| 久久久久久久精品精品| 国产人伦9x9x在线观看 | 纯流量卡能插随身wifi吗| 黄网站色视频无遮挡免费观看| 黄网站色视频无遮挡免费观看| 建设人人有责人人尽责人人享有的| 国产精品国产三级国产专区5o| a 毛片基地| 亚洲国产色片| 91精品伊人久久大香线蕉| 热re99久久国产66热| 好男人视频免费观看在线| 亚洲精品国产av成人精品| 99久国产av精品国产电影| 精品卡一卡二卡四卡免费| 日日爽夜夜爽网站| videossex国产| 我要看黄色一级片免费的| 激情视频va一区二区三区| 人妻少妇偷人精品九色| 亚洲综合精品二区| 校园人妻丝袜中文字幕| 国产成人av激情在线播放| 国产综合精华液| 久久人人爽av亚洲精品天堂| 男女边吃奶边做爰视频| 国产精品欧美亚洲77777| 日韩精品有码人妻一区| 又大又黄又爽视频免费| 777米奇影视久久| 国产成人精品婷婷| 精品午夜福利在线看| av免费在线看不卡| 十分钟在线观看高清视频www| 又大又黄又爽视频免费| 精品久久久久久电影网| 久久久久精品久久久久真实原创| 欧美bdsm另类| 又大又黄又爽视频免费| 又大又黄又爽视频免费| 26uuu在线亚洲综合色| 最黄视频免费看| 亚洲美女黄色视频免费看| 天天躁夜夜躁狠狠久久av| 国产熟女欧美一区二区| 欧美xxⅹ黑人| 久久精品国产亚洲av高清一级| 成人毛片a级毛片在线播放| 精品人妻偷拍中文字幕| 纵有疾风起免费观看全集完整版| 精品国产乱码久久久久久男人| 制服丝袜香蕉在线| 五月开心婷婷网| 午夜福利在线免费观看网站| 国产免费现黄频在线看| 精品国产乱码久久久久久男人| 久久久久久久久久人人人人人人| 寂寞人妻少妇视频99o| a级毛片在线看网站| 中文天堂在线官网| 一级,二级,三级黄色视频| 亚洲一区二区三区欧美精品| 麻豆av在线久日| 亚洲第一区二区三区不卡| 母亲3免费完整高清在线观看 | 国产精品免费大片| 久久ye,这里只有精品| 91在线精品国自产拍蜜月| 99久久综合免费| 热99国产精品久久久久久7| 亚洲欧美中文字幕日韩二区| 国产精品久久久久久精品古装| 女的被弄到高潮叫床怎么办| 国产一级毛片在线| 亚洲少妇的诱惑av| 波多野结衣av一区二区av| 80岁老熟妇乱子伦牲交| 捣出白浆h1v1| 欧美激情高清一区二区三区 | 国产熟女欧美一区二区| 男女啪啪激烈高潮av片| freevideosex欧美| 最黄视频免费看| 精品久久蜜臀av无| 国产片内射在线| 免费在线观看视频国产中文字幕亚洲 | 超色免费av| 亚洲av成人精品一二三区| 少妇人妻久久综合中文| 夫妻午夜视频| 亚洲婷婷狠狠爱综合网| 少妇 在线观看| 丝瓜视频免费看黄片| 午夜免费观看性视频| 国产精品.久久久| 满18在线观看网站| 观看av在线不卡| 国产男女内射视频| 久久国产亚洲av麻豆专区| 欧美国产精品va在线观看不卡| 欧美成人午夜精品| 一个人免费看片子| 深夜精品福利| 免费女性裸体啪啪无遮挡网站| 成人手机av| 久久久久国产一级毛片高清牌| 男女高潮啪啪啪动态图| 在线观看免费高清a一片| 下体分泌物呈黄色| 久久久久久久国产电影| 亚洲三级黄色毛片| 九草在线视频观看| 寂寞人妻少妇视频99o| 日本色播在线视频| 国产精品免费大片| 久久精品久久精品一区二区三区| 女性被躁到高潮视频| 久久精品久久久久久久性| 国产又色又爽无遮挡免| 久久女婷五月综合色啪小说| 亚洲精品视频女| 国产极品粉嫩免费观看在线| √禁漫天堂资源中文www| 国产免费视频播放在线视频| a 毛片基地| 黑人巨大精品欧美一区二区蜜桃| 美女中出高潮动态图| 成人亚洲精品一区在线观看| 天天操日日干夜夜撸| 成人国产麻豆网| 国产亚洲最大av| 欧美日韩精品网址| 久久久久国产一级毛片高清牌| 十分钟在线观看高清视频www| 校园人妻丝袜中文字幕| 99久久精品国产国产毛片| 久久久久国产精品人妻一区二区| 国产爽快片一区二区三区| 欧美变态另类bdsm刘玥| 熟女电影av网| 国产 一区精品| 有码 亚洲区| 久久精品国产综合久久久| 最近2019中文字幕mv第一页| 乱人伦中国视频| 午夜免费男女啪啪视频观看| 99国产综合亚洲精品| 如日韩欧美国产精品一区二区三区| 爱豆传媒免费全集在线观看| 青青草视频在线视频观看| 久久久久精品性色| 欧美成人午夜免费资源| 午夜免费男女啪啪视频观看| 美女脱内裤让男人舔精品视频| 香蕉国产在线看| 午夜久久久在线观看| 一区二区三区乱码不卡18| 女人被躁到高潮嗷嗷叫费观| 伦精品一区二区三区| 亚洲av国产av综合av卡| 777米奇影视久久| 王馨瑶露胸无遮挡在线观看| 国产精品一二三区在线看| videos熟女内射| 性高湖久久久久久久久免费观看| 久久精品亚洲av国产电影网| 欧美国产精品一级二级三级| 国产精品久久久久成人av| av线在线观看网站| 午夜福利,免费看| 日韩av不卡免费在线播放| 国产精品人妻久久久影院| 亚洲国产最新在线播放| 亚洲三区欧美一区| 极品少妇高潮喷水抽搐| 国产精品.久久久| 两性夫妻黄色片| 国产精品欧美亚洲77777| 国产精品熟女久久久久浪| 亚洲国产精品国产精品| 久久精品国产自在天天线| 久久久久人妻精品一区果冻| 中文乱码字字幕精品一区二区三区| 久久久久久人妻| 国产一区二区激情短视频 | 国产成人精品久久二区二区91 | 蜜桃在线观看..| 男的添女的下面高潮视频| 2018国产大陆天天弄谢| 人妻人人澡人人爽人人| 亚洲av在线观看美女高潮| 九草在线视频观看| 国产精品成人在线| 搡老乐熟女国产| 免费日韩欧美在线观看| 精品亚洲成国产av| www.av在线官网国产| 欧美激情高清一区二区三区 | www.自偷自拍.com| 免费人妻精品一区二区三区视频| 亚洲精品成人av观看孕妇| 男人爽女人下面视频在线观看| 韩国精品一区二区三区| 国产在线免费精品| 日韩不卡一区二区三区视频在线| 亚洲色图综合在线观看| 亚洲成人一二三区av| 老汉色av国产亚洲站长工具| 一区二区三区精品91| 国产一区二区 视频在线| 久久精品久久久久久久性| 日韩在线高清观看一区二区三区| 麻豆av在线久日| 在线观看国产h片| 三级国产精品片| 伦精品一区二区三区| 大香蕉久久成人网| 国产精品一区二区在线不卡| 亚洲精品成人av观看孕妇| 亚洲精品一二三| 亚洲av国产av综合av卡| www.熟女人妻精品国产| 国产乱来视频区| 深夜精品福利| 国产精品成人在线| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩另类电影网站| 午夜免费男女啪啪视频观看| 亚洲精品视频女| 免费av中文字幕在线| 免费黄色在线免费观看| 亚洲视频免费观看视频| 这个男人来自地球电影免费观看 | 亚洲av福利一区| av女优亚洲男人天堂| 久久婷婷青草| 亚洲av中文av极速乱| 欧美国产精品一级二级三级| 亚洲欧美中文字幕日韩二区| 成年女人在线观看亚洲视频| 一级片免费观看大全| 多毛熟女@视频| 久久久久国产精品人妻一区二区| 国产成人午夜福利电影在线观看| 国产成人91sexporn| 欧美日韩av久久| 91精品三级在线观看| 国产国语露脸激情在线看| 最近手机中文字幕大全| 老女人水多毛片| 纵有疾风起免费观看全集完整版| 夫妻午夜视频| 亚洲精品成人av观看孕妇| 免费在线观看完整版高清| 国产熟女午夜一区二区三区| 久久99热这里只频精品6学生| 男人舔女人的私密视频| 久久精品人人爽人人爽视色| 日本vs欧美在线观看视频| 亚洲国产最新在线播放| 一级爰片在线观看| 男女啪啪激烈高潮av片| 久久久国产一区二区| 最近的中文字幕免费完整| 国产精品国产av在线观看| av在线观看视频网站免费| 久久影院123| 一级黄片播放器| 人人妻人人爽人人添夜夜欢视频| 少妇被粗大猛烈的视频| 日韩欧美精品免费久久| 久久狼人影院| 只有这里有精品99| 黄片播放在线免费| 久久久久久人妻| 黄频高清免费视频| 夫妻午夜视频| 老汉色∧v一级毛片| 国产在线视频一区二区| 国产一区二区在线观看av| 老司机影院成人| 免费看不卡的av| 成人黄色视频免费在线看| 91精品伊人久久大香线蕉| 天天影视国产精品| 日韩av在线免费看完整版不卡| 亚洲精品久久午夜乱码| 18禁观看日本| 亚洲第一区二区三区不卡| √禁漫天堂资源中文www| 久久久久久免费高清国产稀缺| 女性生殖器流出的白浆| 久久久久久伊人网av| 日韩电影二区| 日韩人妻精品一区2区三区| 在线观看美女被高潮喷水网站| 免费黄网站久久成人精品| av.在线天堂| 男女高潮啪啪啪动态图| 另类精品久久| 国产成人精品在线电影| 国产成人精品婷婷| 精品久久久精品久久久| 欧美精品高潮呻吟av久久| av网站免费在线观看视频| 久久人人爽av亚洲精品天堂| 久久久久久久大尺度免费视频| 熟女电影av网| 国产片内射在线| 精品99又大又爽又粗少妇毛片| 亚洲国产日韩一区二区| 又大又黄又爽视频免费| 高清av免费在线| 亚洲av电影在线观看一区二区三区| 人妻系列 视频| kizo精华| 亚洲伊人久久精品综合| 91在线精品国自产拍蜜月| 美女国产高潮福利片在线看| 在线观看人妻少妇| 国产在线免费精品| 中国三级夫妇交换| 国产片特级美女逼逼视频| 日韩在线高清观看一区二区三区| 青草久久国产| 99九九在线精品视频| 午夜福利乱码中文字幕| 欧美日韩视频精品一区| 亚洲,一卡二卡三卡| videos熟女内射| 2022亚洲国产成人精品| 精品少妇久久久久久888优播| 亚洲少妇的诱惑av| 一级片免费观看大全| 精品少妇久久久久久888优播| 国产亚洲午夜精品一区二区久久| 午夜激情av网站| 久久久久国产网址| 免费久久久久久久精品成人欧美视频| 黄网站色视频无遮挡免费观看| 人体艺术视频欧美日本| 色哟哟·www| 韩国精品一区二区三区| 男的添女的下面高潮视频| 日日摸夜夜添夜夜爱| 男女无遮挡免费网站观看| 国产精品无大码| 不卡av一区二区三区| 国产精品二区激情视频| 精品久久蜜臀av无| 九色亚洲精品在线播放| 深夜精品福利| 国产精品人妻久久久影院| 欧美亚洲日本最大视频资源| 人体艺术视频欧美日本| 9色porny在线观看| 在线观看一区二区三区激情| 91aial.com中文字幕在线观看| 国语对白做爰xxxⅹ性视频网站| a级片在线免费高清观看视频| 少妇人妻 视频| 亚洲欧美日韩另类电影网站| 国产av码专区亚洲av| 国产日韩一区二区三区精品不卡| 精品亚洲成国产av| 国产黄频视频在线观看| 好男人视频免费观看在线| 亚洲人成网站在线观看播放| 精品一区二区免费观看| 亚洲三区欧美一区| 亚洲一区二区三区欧美精品| 欧美人与性动交α欧美软件| 夫妻性生交免费视频一级片| 校园人妻丝袜中文字幕| 日韩一区二区视频免费看| 少妇人妻 视频| 亚洲精品成人av观看孕妇| 欧美最新免费一区二区三区| 久久热在线av| 精品国产一区二区三区四区第35| 免费观看a级毛片全部| 亚洲欧美一区二区三区国产| 亚洲经典国产精华液单| 黑人欧美特级aaaaaa片| 成年女人在线观看亚洲视频| 最黄视频免费看| 少妇人妻精品综合一区二区| 最近中文字幕2019免费版| 青青草视频在线视频观看| 久久精品熟女亚洲av麻豆精品| 免费不卡的大黄色大毛片视频在线观看| 纯流量卡能插随身wifi吗| 男女啪啪激烈高潮av片| 男男h啪啪无遮挡| 深夜精品福利| 日韩av免费高清视频| 波多野结衣av一区二区av| 国产极品粉嫩免费观看在线| 97精品久久久久久久久久精品| 亚洲成人av在线免费| 91在线精品国自产拍蜜月| 一区二区日韩欧美中文字幕| 最近最新中文字幕免费大全7| 女人被躁到高潮嗷嗷叫费观| 久久久亚洲精品成人影院| 亚洲欧美一区二区三区久久| 在线天堂中文资源库| 制服人妻中文乱码| 少妇的逼水好多| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲美女搞黄在线观看| 青草久久国产| 国产成人aa在线观看| 男女免费视频国产| 热re99久久国产66热| 捣出白浆h1v1| 黑人猛操日本美女一级片| 国产精品麻豆人妻色哟哟久久| 人人妻人人澡人人看| 黄色 视频免费看| 18在线观看网站| 人妻少妇偷人精品九色| 天堂中文最新版在线下载| 天天躁日日躁夜夜躁夜夜| 精品少妇黑人巨大在线播放| 国产片内射在线| 80岁老熟妇乱子伦牲交| 国产男人的电影天堂91| 精品亚洲成a人片在线观看| 26uuu在线亚洲综合色| 国产精品久久久久久精品古装| 成人毛片a级毛片在线播放| 不卡av一区二区三区| 亚洲成色77777| 一二三四中文在线观看免费高清| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产一区二区在线观看av| 香蕉精品网在线| 免费人妻精品一区二区三区视频| 亚洲精品第二区| 久久免费观看电影| 老司机亚洲免费影院| 熟女电影av网| 精品卡一卡二卡四卡免费| 国产精品一区二区在线观看99| 999精品在线视频| 亚洲国产精品一区二区三区在线| 性色avwww在线观看| 一本久久精品| 午夜免费男女啪啪视频观看| 自拍欧美九色日韩亚洲蝌蚪91| tube8黄色片| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产精品国产精品| 看非洲黑人一级黄片| 久久久久久久国产电影| 久久久久精品久久久久真实原创| 久久久久久人人人人人| 国产精品无大码| 在线观看www视频免费| 欧美日韩精品成人综合77777| 亚洲欧美日韩另类电影网站| 男女边摸边吃奶| 一级a爱视频在线免费观看| 欧美日韩一区二区视频在线观看视频在线| 国产成人a∨麻豆精品| 人人妻人人澡人人爽人人夜夜| 高清不卡的av网站| 亚洲精品美女久久av网站| 久久97久久精品| 欧美激情极品国产一区二区三区| 午夜免费鲁丝| 国产在线视频一区二区| 最近中文字幕高清免费大全6| 99re6热这里在线精品视频| 成人漫画全彩无遮挡| 久久亚洲国产成人精品v| 国产日韩一区二区三区精品不卡| 亚洲第一青青草原| 精品少妇久久久久久888优播| 亚洲精品美女久久av网站| 亚洲精品av麻豆狂野| 色94色欧美一区二区| 亚洲婷婷狠狠爱综合网| 黑丝袜美女国产一区| 2021少妇久久久久久久久久久| 国产成人精品久久久久久| 丝袜脚勾引网站| 天天操日日干夜夜撸| 中文乱码字字幕精品一区二区三区| 日本-黄色视频高清免费观看| 免费日韩欧美在线观看| 欧美精品一区二区大全| 97人妻天天添夜夜摸| 最新的欧美精品一区二区| 在线观看www视频免费| 美女主播在线视频| 久久精品国产亚洲av天美|