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

    基于Fluent/ISIGHT 貫流式潛水泵流場與結(jié)構(gòu)優(yōu)化研究

    2022-11-03 06:11:00劉國勇徐培智王志明尹忠俊郝鵬鵬朱冬梅
    關(guān)鍵詞:優(yōu)化結(jié)構(gòu)

    劉國勇,徐培智,王志明,尹忠俊,郝鵬鵬,朱冬梅

    (1.北京科技大學(xué)機(jī)械工程學(xué)院,北京 100083;2.湖北航天化學(xué)技術(shù)研究所,湖北襄陽 441001)

    礦用潛水泵[1]作為機(jī)電一體化新型排水設(shè)備,具有不怕水淹、高揚(yáng)程、大流量等優(yōu)勢,在處理透水事故中發(fā)揮積極作用,是應(yīng)對礦井水災(zāi)的重要搶險工具[2].潛水泵的結(jié)構(gòu)優(yōu)化設(shè)計方法發(fā)展至今,從最初的速度系數(shù)法優(yōu)化設(shè)計、損失極值法優(yōu)化設(shè)計,衍化到如今結(jié)合計算機(jī)技術(shù)的計算流體力學(xué)(Computational Fluid Dynamics,CFD)優(yōu)化設(shè)計和多工況優(yōu)化設(shè)計法[3-4],大大縮短了泵的設(shè)計周期.董敏等[5]選擇多目標(biāo)優(yōu)化與數(shù)值模擬相結(jié)合的方法,利用遺傳算法對目標(biāo)函數(shù)進(jìn)行多目標(biāo)優(yōu)化求解.王玉勤等[6]建立以泵汽蝕余量為目標(biāo)函數(shù)的優(yōu)化數(shù)學(xué)模型,采用正交試驗對離心泵進(jìn)行優(yōu)化設(shè)計.Wang等[7]采用直接數(shù)值模擬方法,研究了單向耦合條件下平面湍流壁面射流中慣性粒子的分布和運(yùn)動規(guī)律.邱勇等[8]采用Mixture 多相流模型對雙葉片泵內(nèi)部流動進(jìn)行三維非定常模擬計算,分析了不同工況和顆粒物性下泵徑向力特性的變化規(guī)律.曹衛(wèi)東等[9]基于CFX軟件,采用標(biāo)準(zhǔn)k-ε模型對裝配了不同尺寸后蓋板的葉輪、導(dǎo)葉、進(jìn)出水段的礦用搶險多級排水泵的全流場模型進(jìn)行了數(shù)值模擬.21 世紀(jì)以來,隨著計算機(jī)技術(shù)的高速發(fā)展,CFD 也隨之日趨成熟[10].Wang 等[11]采用混合多相流模型、RNGk-ε湍流模型和協(xié)調(diào)一致壓力耦合方程的半隱式方法(Semi-Implicit Method for Pressure Linked Equations Consistent,SIMPLEC)算法對離心泵的兩相流進(jìn)行了模擬,分析了顆粒尺寸對內(nèi)部壓力分布、固體體積分布和外部特性的影響.以上關(guān)于潛水泵的結(jié)構(gòu)優(yōu)化的研究仍然停留在對單一結(jié)構(gòu)參數(shù)影響分析,對于泵內(nèi)流場的復(fù)雜流域結(jié)構(gòu),必須綜合多種結(jié)構(gòu)參數(shù)整體分析結(jié)構(gòu)影響特點(diǎn),從而得出最優(yōu)化的結(jié)果.

    本文針對貫流式潛水泵的流場特性進(jìn)行了系統(tǒng)分析,提出一種綜合揚(yáng)程、效率和輕量化目標(biāo)的多結(jié)構(gòu)參數(shù)整體優(yōu)化方案,基于Fluent 流體仿真軟件建立貫流式潛水泵全流場仿真模型.分析葉輪結(jié)構(gòu)參數(shù)對潛水泵流場流動特性的影響,通過ISIGHT 優(yōu)化平臺實(shí)現(xiàn)參數(shù)化潛水泵建模仿真過程,采用智能優(yōu)化算法對潛水泵結(jié)構(gòu)參數(shù)進(jìn)行多目標(biāo)多因素優(yōu)化,確定具體的優(yōu)化結(jié)構(gòu).本文研究旨在優(yōu)化出一種能滿足礦井水突事故搶險救災(zāi)需求的大流量輕量化高效礦用潛水泵.

    1 CFD仿真模型建立

    1.1 流場模型建立

    貫流式潛水泵吸水方式是由下至上,在搶險吸水的過程中能夠為電動機(jī)發(fā)熱段提供持續(xù)冷卻散熱,保證電動機(jī)長時間工作的穩(wěn)定性和安全性.在應(yīng)對不同復(fù)雜搶險救災(zāi)環(huán)境時可以根據(jù)需求將多個潛水泵串聯(lián),實(shí)現(xiàn)增大揚(yáng)程的效果,避免多級潛水泵的體型笨重和運(yùn)輸困難問題,可以柔性應(yīng)用于各種復(fù)雜搶險情況.

    貫流式潛水泵結(jié)構(gòu)示意圖如圖1 所示,主要包括導(dǎo)葉體、葉輪、進(jìn)水段、電動機(jī)以及主軸等機(jī)械部件.基于這種貫流式潛水泵結(jié)構(gòu)通過SpaceClaim 中的體積抽取工具獲取水體結(jié)構(gòu)模型,體積抽取工具能夠完好地獲取流體機(jī)械內(nèi)部的封閉流體域,然后對支撐結(jié)構(gòu)進(jìn)行簡化,得到潛水泵的流場水體模型.貫流式潛水泵流場三維模型和截面圖如圖2所示.

    圖1 貫流式潛水泵結(jié)構(gòu)示意圖Fig.1 Schematic diagram of the structure of a tubular submersible pump

    圖2 貫流式潛水泵流場三維模型和截面圖Fig.2 Three-dimensional model of flow field and sectional view of tubular submersible pump

    在潛水泵流體域中,為了引導(dǎo)葉輪水流旋轉(zhuǎn)方向,在貫流式蝸殼結(jié)構(gòu)與葉輪出水口外添加軸向?qū)~,將流體引流成軸向方向,提高流體結(jié)構(gòu)的葉輪旋轉(zhuǎn)能量利用效率,從而起到提高揚(yáng)程和效率的作用.數(shù)值計算過程中將整體貫流式潛水泵水體結(jié)構(gòu)劃分為進(jìn)水段、葉輪和外流道三部分,如圖3所示.

    圖3 貫流式潛水泵水體結(jié)構(gòu)三維模型Fig.3 Three-dimensional model of water structure of tubular submersible pump

    1.2 網(wǎng)格劃分與無關(guān)性驗證

    采用ICEM-CFD 軟件分別對進(jìn)水段、葉輪以及外流道流域進(jìn)行網(wǎng)格劃分.進(jìn)水段采用O-Grid 方法劃分為六面體結(jié)構(gòu)網(wǎng)格,因葉輪和外流道流域為不規(guī)則的空間曲面結(jié)構(gòu),故采用非結(jié)構(gòu)四面體網(wǎng)格,然后通過Make Polyhedra 功能轉(zhuǎn)化為多面體網(wǎng)格.各水體部分網(wǎng)格模型如圖4所示.

    圖4 貫流式潛水泵流場網(wǎng)格模型Fig.4 Grid model of flow field of tubular submersible pump

    流場數(shù)值模擬精度受網(wǎng)格數(shù)量和質(zhì)量的影響,而網(wǎng)格數(shù)量和質(zhì)量又受到計算機(jī)資源配置的限制.因此,需要進(jìn)行網(wǎng)格獨(dú)立性驗證.采用Fluent 流體仿真軟件對貫流式潛水泵進(jìn)行穩(wěn)態(tài)定常計算,根據(jù)網(wǎng)格數(shù)量與揚(yáng)程關(guān)系來驗證網(wǎng)格無關(guān)性.結(jié)果如表1所示.

    由表1 可知,當(dāng)網(wǎng)格數(shù)量大于146 萬時,潛水泵揚(yáng)程的誤差值在1.0%以內(nèi),基本可以認(rèn)為網(wǎng)格數(shù)量達(dá)到146 萬后仿真結(jié)果與網(wǎng)格數(shù)量無關(guān).考慮計算準(zhǔn)確性和計算資源配置合理性,選取方案d 作為后續(xù)仿真模擬的網(wǎng)格數(shù)量.

    表1 網(wǎng)格無關(guān)性驗證Tab.1 Grid independence validation

    1.3 邊界條件及求解器設(shè)置

    采用Fluent 軟件對貫流式潛水泵進(jìn)行全流場仿真,全局設(shè)置重力加速度9.81 m/s2,方向為軸方向,在穩(wěn)態(tài)計算中,葉輪流體區(qū)域設(shè)置為框架運(yùn)動中的旋轉(zhuǎn)運(yùn)動,旋轉(zhuǎn)方向由右手定則確定,轉(zhuǎn)速N為3 600 r/min.

    進(jìn)水口:本文的潛水泵設(shè)計額定流量為240 m3/h,采用速度入口,速度大小v根據(jù)流量和入口截面確定:

    式中:Q為泵的流量;ρ為流體的密度;S為截面面積;R為截面半徑.

    出水口采用壓力出口.靜壁面采用無滑移壁面,選用標(biāo)準(zhǔn)粗糙模型,粗糙度值設(shè)置為0.025 mm.旋轉(zhuǎn)壁面的葉輪部分選用旋轉(zhuǎn)運(yùn)動,采用無滑移壁面.交界面選擇進(jìn)水段與葉輪進(jìn)口的交界面及葉輪出口與貫流式外流道入口的交界面.計算定常過程的流體主要為液態(tài)水,選用壓力基礎(chǔ)求解器,求解算法選用壓力耦合方程的半隱式方法(Semi-Implicit Method for Pressure Linked Equations,SIMPLE)算法,插值計算方法平衡考慮精度和計算量后選用二階迎風(fēng)格式.初始化進(jìn)水口,設(shè)定計算步數(shù),開始計算直到結(jié)果收斂.

    2 計算結(jié)果與分析

    2.1 葉輪結(jié)構(gòu)參數(shù)定義

    本文潛水泵水力性能的主要指標(biāo)為揚(yáng)程與水力效率,葉輪為中心對稱結(jié)構(gòu).葉輪參數(shù)結(jié)構(gòu)定義如圖5所示.

    由圖5 可知,潛水泵的葉輪結(jié)構(gòu)參數(shù)主要有6個:①葉輪直徑d2,葉輪的主要尺寸指標(biāo);②葉片包角φ,葉片兩端與軸心連線形成的夾角;③出口安放角β2,葉片出口切線與葉輪圓周法線夾角;④葉片厚度s,葉片厚度設(shè)為均勻厚度;⑤出口寬度b2,葉輪出口的寬度,也可以認(rèn)為葉輪出口的厚度;⑥葉片數(shù)量n.

    圖5 葉輪結(jié)構(gòu)參數(shù)定義Fig.5 Impeller structure parameters definition

    CFturbo 旋轉(zhuǎn)機(jī)械建模軟件按照實(shí)際潛水泵流量240 m3/h、揚(yáng)程120 m、轉(zhuǎn)速3 600 r/min要求得到葉輪結(jié)構(gòu)經(jīng)驗參數(shù),如表2 所示.后續(xù)的葉輪結(jié)構(gòu)參數(shù)在此基礎(chǔ)上進(jìn)行優(yōu)化和分析.

    表2 CFturbo葉輪結(jié)構(gòu)經(jīng)驗參數(shù)Tab.2 CFturbo experiential impeller structure parameters

    2.2 潛水泵性能參數(shù)

    潛水泵的主要水力性能指標(biāo)為揚(yáng)程與水力效率,其詳細(xì)說明如下.

    1)揚(yáng)程是指水泵能夠揚(yáng)水的高度,通常以符號H來表示,單位為m,計算公式為:

    式中:P1、P2分別為潛水泵進(jìn)、出口處液體總壓力;v1、v2分別為潛水泵進(jìn)、出口處液體的絕對速度;Z1、Z2分別為潛水泵進(jìn)、出口處相對任一基準(zhǔn)垂直方向距離的位置;g為重力加速度.

    2)水力效率為有效功率和軸功率之比,用η表示,其公式表達(dá)如下:

    式中:M為葉輪軸輸入力矩,N·m;ω為葉輪角速度,rad/s.

    2.3 葉輪直徑對潛水泵性能的影響

    建立葉輪直徑d2分別為266 mm、276 mm、286 mm 的仿真模型,通過潛水泵定常仿真模擬,分析不同葉輪直徑的流場影響規(guī)律.

    圖6 為不同流量工況下葉輪直徑對揚(yáng)程和水力效率影響曲線.由圖6 可知,當(dāng)流量為160~440 m3/h時,隨著流量增大,流速增大,導(dǎo)致水泵內(nèi)部水流消耗揚(yáng)程增加,有效揚(yáng)程降低,所以揚(yáng)程隨著葉輪直徑的增大而增大.當(dāng)流量為160~360 m3/h 時,水力效率隨著葉輪直徑增大而降低.當(dāng)流量為360~440 m3/h時,水力效率基本隨著葉輪直徑增大而增大,水力效率極值點(diǎn)隨著葉輪直徑的增大逐漸向小流量區(qū)偏移.隨著葉輪直徑增大,葉片尺寸也增大,葉片掃掠面積增大,對于流體的驅(qū)動力更大,所以揚(yáng)程會提升.當(dāng)潛水泵中的流量逐漸增大時,葉輪中的流體速度增大,單位時間內(nèi)流體和葉片表面接觸的摩擦力增大,導(dǎo)致流體局部摩擦能量損失增加.由圖6(a)可以看出,公式(3)中流量和揚(yáng)程的乘積大小是隨著流量的增大而減小的,所以,葉輪在大流量范圍的水力效率降低.

    圖6 不同流量工況下葉輪直徑對揚(yáng)程和水力效率影響曲線Fig.6 Influence curve of impeller diameter on head and efficiency under different flow conditions

    2.4 葉輪包角對潛水泵性能的影響

    葉片包角是影響葉輪對流體控制能力的重要結(jié)構(gòu)參數(shù),為了分析葉輪結(jié)構(gòu)的影響規(guī)律,建立葉片包角φ分別為116°、120°、124°的3 組不同潛水泵全流場仿真模型.

    圖7 為不同流量工況下葉片包角對揚(yáng)程和水力效率影響曲線.由圖7 可知,當(dāng)流量為160~440 m3/h時,揚(yáng)程隨著葉片包角的增大而減小,隨著流量增大而減小.這是因為葉片包角影響了葉輪出口液體的滑移,包角的增大使液體受到葉片更強(qiáng)的約束,滑移程度降低,隨著流量增大,揚(yáng)程下降速度加快.

    圖7 不同流量工況下葉片包角對揚(yáng)程和水力效率影響曲線Fig.7 Influence curve of impeller blade wrap angles on head and efficiency under different flow conditions

    當(dāng)流量為160~360 m3/h 時,水力效率隨著葉片包角增大而升高;當(dāng)流量為360~440 m3/h 時,水力效率基本隨著葉片包角增大而降低,效率極值點(diǎn)隨著葉片包角的增大逐漸向小流量區(qū)偏移.這主要是因為葉片掃掠接觸面積增大,對進(jìn)出口流體的控制力更強(qiáng),但是摩擦阻力耗能增加更明顯.在小流量范圍內(nèi),由于包角的增大,葉輪流道內(nèi)的流動擴(kuò)散減小,流動更貼近葉片形狀,葉片控制力作用明顯,水利損失較小,所以隨著包角增大,能夠提高潛水泵整體的效率;但是隨著流量增大,單位時間內(nèi)流體與葉片表面接觸產(chǎn)生的摩擦力增大,摩擦能量損失增速加快,導(dǎo)致水力效率開始降低.

    2.5 葉輪出口寬度對潛水泵性能的影響

    葉輪出口寬度是控制葉輪流場空間的重要參數(shù).建立出口寬度b2分別為20 mm、22 mm、24 mm 的3組不同葉輪流場仿真模型,對葉輪的出口寬度影響規(guī)律進(jìn)行分析.

    圖8 為不同流量工況下葉輪出口寬度對揚(yáng)程和效率影響曲線.由圖8 可知,當(dāng)流量為160~440 m3/h時,揚(yáng)程隨著葉輪出口寬度的增大而增大,隨著流量的增加而減小,葉輪出口寬度越大,揚(yáng)程降低速率越小.當(dāng)流量為160~440 m3/h 時,效率隨著葉輪出口寬度增大而減小,效率極值點(diǎn)基本保持在300 m3/h 左右.葉輪出口寬度的大小決定葉輪內(nèi)部的流體空間和葉片軸向方向的寬度,流體空間增大,葉片與流體接觸面積就會增大,導(dǎo)致葉片表面的沿程阻力能量損失提高,加上葉輪出口處的局部阻力能量損失提高,產(chǎn)生了出口寬度越大,潛水泵效率越低的現(xiàn)象.

    圖8 不同流量工況下葉輪出口寬度對揚(yáng)程和效率影響曲線Fig.8 Influence curve of impeller outlet width on head and efficiency under different flow conditions

    2.6 葉片數(shù)量對潛水泵性能的影響

    葉片數(shù)量的變化能夠很明顯地改變?nèi)~輪內(nèi)部的流場結(jié)構(gòu).建立葉片數(shù)量n分別為4、5、6 的3 組不同的葉輪流場仿真模型.

    圖9 為不同葉片數(shù)量下?lián)P程和水力效率隨流量變化曲線.由圖9可知,當(dāng)流量為160~440 m3/h時,揚(yáng)程隨葉片數(shù)量的增加而增大,但隨著流量的增大,揚(yáng)程增加的幅度減小,符合壓力云圖的入口負(fù)壓值增大的流場特性.當(dāng)流量為160~440 m3/h 時,水力效率隨葉片數(shù)量增加而減小,效率極值點(diǎn)隨著葉片數(shù)量的增加逐漸向小流量區(qū)偏移.葉片數(shù)量的增加能夠顯著提高葉輪對流體的做功能力,提高潛水泵的揚(yáng)程特性.但是,由于增加了葉片數(shù)量,導(dǎo)致葉輪中葉片與流體接觸面積增大,產(chǎn)生的沿程阻力耗能增加,從而導(dǎo)致潛水泵整體的效率下降.由圖9 可知,較多的葉片能夠提高揚(yáng)程,但是只適用于流量工況不是很大的情況.

    圖9 不同葉片數(shù)量下?lián)P程和水力效率隨流量變化曲線Fig.9 Influence curve of impeller number on head nd efficiency under different flow conditions

    3 葉輪結(jié)構(gòu)參數(shù)優(yōu)化

    3.1 ISIGHT優(yōu)化流程

    ISIGHT 優(yōu)化平臺是一個模塊化的流程軟件[12],可以通過調(diào)用其中的模塊將多個外部軟件集合到優(yōu)化平臺中.本文選用CFturbo 和PumpLinx 作為ISIGHT 優(yōu)化平臺中潛水泵葉輪結(jié)構(gòu)建模和仿真軟件.主要原因是,CFturbo 和PumpLinx 軟件均能夠把參數(shù)結(jié)構(gòu)設(shè)置記錄下來,生成批處理文件(batch 文件),可以在結(jié)構(gòu)優(yōu)化過程中節(jié)省大量編寫程式的時間.同時,CFturbo 與PumpLinx 之間具有非常優(yōu)秀的接口關(guān)系,CFturbo 建立的潛水泵模型能夠直接從CFturbo建模軟件中啟動PumpLinx流體仿真軟件,并且自動選擇相應(yīng)的泵模板和求解模型,直接進(jìn)入PumpLinx 的計算階段.ISIGHT 可以將其中的關(guān)鍵結(jié)構(gòu)參數(shù)化,便于智能算法進(jìn)行結(jié)構(gòu)優(yōu)化和更改.圖10為ISIGHT優(yōu)化流程.

    圖10 ISIGHT優(yōu)化流程Fig.10 ISIGHT optimization flow

    由圖10 可知,計算器模塊統(tǒng)一整個流程中的參數(shù)變量;CFturbo 讀取寫入模塊負(fù)責(zé)batch 文件的讀取和寫入;CFturbo 軟件模塊主要負(fù)責(zé)調(diào)用軟件;PumpLinx 讀取寫入模塊負(fù)責(zé)將spro 文件讀取和寫入,將其中的迭代步數(shù)、收斂條件和網(wǎng)格尺寸設(shè)置成參數(shù)變量;PumpLinx軟件模塊調(diào)用軟件;仿真結(jié)果數(shù)據(jù)提取模塊負(fù)責(zé)計算揚(yáng)程和效率所需的進(jìn)出口總壓和葉輪力矩;輸出結(jié)果讀取模塊將輸出的總壓和力矩結(jié)果讀取并參數(shù)化寫入優(yōu)化流程中;計算器模塊2將結(jié)果轉(zhuǎn)化為揚(yáng)程和效率;優(yōu)化算法模塊能監(jiān)測所有參數(shù)化的數(shù)據(jù),并且設(shè)置葉輪結(jié)構(gòu)參數(shù)的優(yōu)化范圍、約束條件和結(jié)果優(yōu)化目標(biāo)[13].

    整個優(yōu)化流程通過優(yōu)化算法模塊進(jìn)行調(diào)控,仿真結(jié)束后提取關(guān)鍵數(shù)據(jù)計算揚(yáng)程和效率,然后將性能結(jié)果輸入優(yōu)化模塊中,通過智能算法計算下一步優(yōu)化參數(shù),從而繼續(xù)修改葉輪結(jié)構(gòu)參數(shù),實(shí)現(xiàn)自動化、智能化的高效優(yōu)化過程.

    3.2 優(yōu)化數(shù)值模型

    3.2.1 目標(biāo)函數(shù)

    本文有3個優(yōu)化目標(biāo):目標(biāo)1是提高潛水泵設(shè)計工況下的揚(yáng)程f1(x),目標(biāo)2 是提高潛水泵設(shè)計工況下的效率f2(x),目標(biāo)3 是減小潛水泵葉輪結(jié)構(gòu)的體積從而實(shí)現(xiàn)輕量化f3(x).因此,可以將整個優(yōu)化流程簡化為求極值的問題[14]:

    3.2.2 設(shè)計變量

    設(shè)計變量也就是潛水泵結(jié)構(gòu)優(yōu)化流程中的結(jié)構(gòu)參數(shù),包括葉輪直徑、葉片包角、葉片出口安放角、葉片厚度、葉輪出口寬度和葉片數(shù)量.設(shè)計變量表達(dá)式為:

    3.2.3 約束條件

    約束條件是優(yōu)化過程中的必要限制,可以分為性能約束和側(cè)面約束.在本文中,性能約束為潛水泵的基本性能要求,即揚(yáng)程120 m、效率65%;側(cè)面約束則是設(shè)計變量的優(yōu)化范圍,也就是給定的結(jié)構(gòu)參數(shù)的優(yōu)化范圍.設(shè)計變量取值范圍如表3所示.

    表3 設(shè)計變量取值范圍Tab.3 Value range of design variables

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

    多島遺傳算法能夠明顯提升潛水泵輕量化,并且能夠明顯優(yōu)化揚(yáng)程和效率.基于本文輕量化優(yōu)化目標(biāo),確定選用多島遺傳算法進(jìn)行潛水泵性能優(yōu)化和研究.

    通過多島遺傳算法設(shè)置5 個島,6 次遺傳,10 次迭代優(yōu)化,最終優(yōu)化次數(shù)為300 次.潛水泵優(yōu)化結(jié)果散點(diǎn)分布圖如圖11所示.

    圖11 潛水泵優(yōu)化結(jié)果散點(diǎn)分布圖Fig.11 Scatter distribution diagram of submersible pump optimization results

    優(yōu)化迭代后得到貫流式潛水泵優(yōu)化仿真結(jié)果,揚(yáng)程為124.4 m、水力效率為68.1%、葉輪直徑為270.7 mm.優(yōu)化前后設(shè)計變量對比如表4 所示.額定流量工況下,貫流式潛水泵揚(yáng)程、水力效率、葉輪直徑對比結(jié)果如表5 所示,生成的外特性性能曲線如圖12所示.

    圖12 優(yōu)化前后潛水泵外特性性能曲線Fig.12 The external characteristic curve of submersible pump before and after optimization

    表4 優(yōu)化前后設(shè)計變量對比Tab.4 Comparison of design variables before and after optimization

    表5 優(yōu)化前后貫流式潛水泵目標(biāo)參數(shù)對比Tab.5 Comparison of target parameters of tubular submersible pump before and after optimization

    由表5 可知,葉輪參數(shù)優(yōu)化前后潛水泵額定流量工況的揚(yáng)程提升5.1%,水力效率提升2.1%,葉輪直徑減小1.9%,3 個優(yōu)化目標(biāo)都有非常明顯的優(yōu)化提升.由圖12 可知,優(yōu)化后的揚(yáng)程在所有流量工況下都明顯優(yōu)于優(yōu)化前,但是優(yōu)化后揚(yáng)程曲線在高流量工況下的下降趨勢更快.優(yōu)化后的水力效率曲線的極值點(diǎn)在280 m3/h 流量附近,優(yōu)化前的水力效率曲線的極值點(diǎn)在280~320 m3/h 流量內(nèi);優(yōu)化后的水力效率總體優(yōu)于優(yōu)化前,并且水力效率曲線的極值點(diǎn)更靠近240 m3/h的額定流量工況.

    3.4 試驗驗證

    為了驗證數(shù)值模擬的精度及準(zhǔn)確性,由潛水泵廠家根據(jù)優(yōu)化結(jié)果設(shè)計生產(chǎn)了QW-240-120 排水搶險潛水泵.使用流量計和壓力表對該潛水泵進(jìn)行現(xiàn)場測試(圖13),在水溫25 ℃時,測得流量為240.65 m3/h,進(jìn)出口壓差值為1.213 kPa,可計算得出揚(yáng)程為124.78 m.優(yōu)化后潛水泵結(jié)果與現(xiàn)場測試值對比如表6所示.

    圖13 QW-240-120潛水泵測量數(shù)據(jù)Fig.13 QW-240-120 submersible pump measurement data

    由表6 可知,在流量誤差為0.27%的情況下,揚(yáng)程誤差為0.31%,差值較小,可以說明優(yōu)化結(jié)果良好,潛水泵性能提升明顯.

    表6 試驗值與優(yōu)化值對比Tab.6 Comparison of experimental values and optimized values

    4 結(jié)論

    基于Fluent 流體軟件建立貫流式潛水泵全流場仿真模型,在CFturbo 旋轉(zhuǎn)機(jī)械建模軟件生成的經(jīng)驗葉輪結(jié)構(gòu)基礎(chǔ)上,對葉輪結(jié)構(gòu)引起的流場特性和性能特性進(jìn)行定常仿真模擬分析,得出結(jié)論如下:

    1)葉輪直徑增大能夠顯著提高潛水泵的揚(yáng)程,提升大流量工況下潛水泵的水力效率,在保證潛水泵輕量化要求下,可以選用較大的葉輪直徑.葉片包角增大會降低潛水泵的揚(yáng)程,但是,在小流量工況下能夠明顯提升潛水泵的水力效率,并且最優(yōu)效率點(diǎn)向小流量區(qū)偏移.由于葉輪出口寬度增大會增大潛水泵的揚(yáng)程,降低水力效率,所以需要按照實(shí)際需求來選擇合適的葉輪出口寬度.葉片數(shù)量增加能夠明顯提高潛水泵的揚(yáng)程,但是,水力效率會降低,必須綜合揚(yáng)程和水力效率要求選擇合適的葉片數(shù)量.

    2)實(shí)現(xiàn)ISIGHT 優(yōu)化平臺的集成參數(shù)化建模和智能算法多因素多目標(biāo)參數(shù)優(yōu)化方案,通過多島遺傳算法5 個島,6 次遺傳,10 次迭代優(yōu)化,總次數(shù)300次的參數(shù)優(yōu)化過程,得到最優(yōu)的葉輪參數(shù)結(jié)構(gòu)為:葉輪直徑270.67 mm、出口寬度24.6 mm、葉片數(shù)量6個、葉片厚度4.5 mm、出口安放角27.6°、葉片包角116.4°.優(yōu)化后結(jié)果在額定流量工況下?lián)P程提升5.1%,水力效率提升2.1%,葉輪直徑減小1.9%,不同流量工況下,優(yōu)化后的潛水泵揚(yáng)程和水力效率總體優(yōu)于優(yōu)化前.在額定流量為240 m3/h 工況下,通過現(xiàn)場測定,得到進(jìn)出口壓差值為1.213 kPa、揚(yáng)程為124.78 m,優(yōu)化仿真揚(yáng)程為124.4 m,差值僅為0.31%.說明通過ISIGHT 優(yōu)化平臺得出的優(yōu)化方案可靠,能夠提升潛水泵性能.

    猜你喜歡
    優(yōu)化結(jié)構(gòu)
    超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    民用建筑防煙排煙設(shè)計優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    論結(jié)構(gòu)
    中華詩詞(2019年7期)2019-11-25 01:43:04
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    論《日出》的結(jié)構(gòu)
    創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長
    欧美色视频一区免费| 国产精品一区二区三区四区久久| 国产激情偷乱视频一区二区| 国产人妻一区二区三区在| eeuss影院久久| 精品福利观看| 日韩欧美在线二视频| 国产黄色小视频在线观看| 中文亚洲av片在线观看爽| 国产高清视频在线观看网站| 一级av片app| 免费高清视频大片| 精品一区二区三区人妻视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜影院日韩av| 一进一出抽搐gif免费好疼| 亚洲 国产 在线| 亚洲激情在线av| 欧美性感艳星| 亚洲人与动物交配视频| 99热这里只有是精品50| 亚洲中文字幕日韩| av中文乱码字幕在线| 亚洲熟妇中文字幕五十中出| 高清毛片免费观看视频网站| 国内久久婷婷六月综合欲色啪| 中文资源天堂在线| 人人妻,人人澡人人爽秒播| 在线免费观看的www视频| 亚洲av第一区精品v没综合| 简卡轻食公司| 丁香六月欧美| ponron亚洲| 日韩欧美一区二区三区在线观看| 天美传媒精品一区二区| 两人在一起打扑克的视频| 高潮久久久久久久久久久不卡| 成人欧美大片| 99久久久亚洲精品蜜臀av| 免费在线观看亚洲国产| 久久亚洲精品不卡| 最后的刺客免费高清国语| 成人永久免费在线观看视频| 好男人在线观看高清免费视频| 一区福利在线观看| 日韩中文字幕欧美一区二区| 午夜福利成人在线免费观看| 搡老妇女老女人老熟妇| 亚洲内射少妇av| 神马国产精品三级电影在线观看| 看黄色毛片网站| 国产亚洲欧美98| 最近最新中文字幕大全电影3| 亚洲av中文字字幕乱码综合| 日本免费一区二区三区高清不卡| 麻豆国产97在线/欧美| 一本久久中文字幕| 国产成人a区在线观看| 网址你懂的国产日韩在线| 人人妻,人人澡人人爽秒播| 宅男免费午夜| 午夜精品在线福利| 免费观看精品视频网站| 日本免费一区二区三区高清不卡| 淫秽高清视频在线观看| 老司机午夜十八禁免费视频| 中文资源天堂在线| 日韩成人在线观看一区二区三区| 欧美激情在线99| 夜夜夜夜夜久久久久| 免费一级毛片在线播放高清视频| 亚洲人成网站在线播放欧美日韩| 亚洲激情在线av| 亚洲黑人精品在线| 欧洲精品卡2卡3卡4卡5卡区| 国产精品98久久久久久宅男小说| 欧美午夜高清在线| 午夜免费激情av| 欧美另类亚洲清纯唯美| 亚洲色图av天堂| 精品人妻视频免费看| 久久久久久久午夜电影| 99热6这里只有精品| 国产色婷婷99| 欧美一区二区精品小视频在线| 99国产极品粉嫩在线观看| 亚洲精品粉嫩美女一区| 国产麻豆成人av免费视频| 精品久久久久久久久久免费视频| 国内精品美女久久久久久| 性色avwww在线观看| 一本一本综合久久| 国产精品,欧美在线| 精品国产亚洲在线| 小说图片视频综合网站| 夜夜躁狠狠躁天天躁| 成人av在线播放网站| 精品人妻熟女av久视频| 亚洲国产色片| 国内毛片毛片毛片毛片毛片| 免费人成在线观看视频色| 久久久国产成人免费| 亚洲精品亚洲一区二区| 国产精品不卡视频一区二区 | 老熟妇仑乱视频hdxx| 亚洲 欧美 日韩 在线 免费| 欧美丝袜亚洲另类 | 欧美+亚洲+日韩+国产| 成人毛片a级毛片在线播放| 美女免费视频网站| 狂野欧美白嫩少妇大欣赏| 亚洲av日韩精品久久久久久密| 一个人免费在线观看电影| 男女之事视频高清在线观看| 69人妻影院| 久久久久国内视频| 一级作爱视频免费观看| 国产日本99.免费观看| 亚洲不卡免费看| 亚洲三级黄色毛片| 国产精品野战在线观看| 国产免费男女视频| 亚洲国产高清在线一区二区三| 亚洲 国产 在线| 国产欧美日韩一区二区三| 亚洲熟妇熟女久久| 美女xxoo啪啪120秒动态图 | 午夜视频国产福利| 免费看日本二区| 久久午夜福利片| 国产成人aa在线观看| 极品教师在线视频| 亚洲avbb在线观看| 国产视频内射| a级毛片a级免费在线| 女生性感内裤真人,穿戴方法视频| 亚洲av日韩精品久久久久久密| 国内毛片毛片毛片毛片毛片| 午夜福利在线观看免费完整高清在 | 99久久精品一区二区三区| 午夜免费男女啪啪视频观看 | 最近最新免费中文字幕在线| 中亚洲国语对白在线视频| 欧美黑人巨大hd| 91久久精品国产一区二区成人| 岛国在线免费视频观看| 欧美黑人巨大hd| 中国美女看黄片| 麻豆国产av国片精品| 久久久久九九精品影院| 十八禁人妻一区二区| 最近中文字幕高清免费大全6 | 久久精品国产自在天天线| 国产真实伦视频高清在线观看 | 一二三四社区在线视频社区8| 美女xxoo啪啪120秒动态图 | 网址你懂的国产日韩在线| 久久精品91蜜桃| 97碰自拍视频| 亚洲av.av天堂| 国产亚洲欧美在线一区二区| 舔av片在线| 日本 欧美在线| 最后的刺客免费高清国语| av女优亚洲男人天堂| bbb黄色大片| 国产精品久久久久久久电影| 人人妻,人人澡人人爽秒播| 九九久久精品国产亚洲av麻豆| 免费人成视频x8x8入口观看| 婷婷丁香在线五月| 国产伦一二天堂av在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕熟女人妻在线| 99国产综合亚洲精品| 精品久久久久久久末码| 午夜福利视频1000在线观看| 国内精品久久久久久久电影| 又爽又黄a免费视频| 色5月婷婷丁香| 国产高清视频在线播放一区| 亚洲av.av天堂| 久久久久久久午夜电影| 俺也久久电影网| 国产亚洲欧美98| 99热精品在线国产| 国产精品98久久久久久宅男小说| 国产精品影院久久| 久久这里只有精品中国| 琪琪午夜伦伦电影理论片6080| 少妇裸体淫交视频免费看高清| 一区二区三区免费毛片| 最新在线观看一区二区三区| 又爽又黄a免费视频| 中文资源天堂在线| 国产毛片a区久久久久| 国产三级中文精品| 99久久精品热视频| 制服丝袜大香蕉在线| 性色avwww在线观看| 欧美丝袜亚洲另类 | 日日夜夜操网爽| 成年女人看的毛片在线观看| 美女高潮喷水抽搐中文字幕| 久久精品国产99精品国产亚洲性色| 国产精品1区2区在线观看.| 十八禁网站免费在线| 在线免费观看不下载黄p国产 | 草草在线视频免费看| 久久久国产成人免费| 老司机深夜福利视频在线观看| www.色视频.com| 男插女下体视频免费在线播放| 真实男女啪啪啪动态图| 欧美日本亚洲视频在线播放| 久久性视频一级片| 亚洲最大成人中文| 亚洲熟妇中文字幕五十中出| 亚洲不卡免费看| 三级国产精品欧美在线观看| 亚洲黑人精品在线| 久久久久免费精品人妻一区二区| 一区福利在线观看| 国产中年淑女户外野战色| 又黄又爽又免费观看的视频| 欧美一区二区国产精品久久精品| 尤物成人国产欧美一区二区三区| 无遮挡黄片免费观看| 日韩欧美精品免费久久 | 69人妻影院| 亚洲欧美日韩高清在线视频| 国产三级中文精品| 美女cb高潮喷水在线观看| 在线播放无遮挡| 亚洲av.av天堂| 久久久久性生活片| 亚洲欧美日韩东京热| 久久亚洲精品不卡| 午夜激情欧美在线| 成人特级黄色片久久久久久久| 最近视频中文字幕2019在线8| 一本精品99久久精品77| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕精品亚洲无线码一区| 一个人看的www免费观看视频| 欧美日韩中文字幕国产精品一区二区三区| 观看免费一级毛片| 亚洲片人在线观看| 非洲黑人性xxxx精品又粗又长| 国产黄片美女视频| 国产高清有码在线观看视频| 在线播放无遮挡| 精品午夜福利在线看| 午夜激情福利司机影院| 老司机深夜福利视频在线观看| АⅤ资源中文在线天堂| 蜜桃久久精品国产亚洲av| 欧美又色又爽又黄视频| 精品欧美国产一区二区三| 欧美精品啪啪一区二区三区| 成人特级av手机在线观看| 国产av一区在线观看免费| av天堂在线播放| 久久九九热精品免费| 日本熟妇午夜| 麻豆国产97在线/欧美| 精品久久久久久久久久久久久| 国产真实伦视频高清在线观看 | 精品国内亚洲2022精品成人| 美女高潮喷水抽搐中文字幕| 中文在线观看免费www的网站| 成人三级黄色视频| 国产视频内射| 人人妻,人人澡人人爽秒播| 亚洲成av人片免费观看| 国产精品亚洲一级av第二区| 一个人观看的视频www高清免费观看| 亚洲成av人片在线播放无| 97超视频在线观看视频| 在线播放国产精品三级| avwww免费| 桃色一区二区三区在线观看| 人人妻,人人澡人人爽秒播| 日本a在线网址| www.熟女人妻精品国产| 亚洲专区中文字幕在线| 精品欧美国产一区二区三| 亚洲专区国产一区二区| 亚洲av第一区精品v没综合| 国产在线男女| 日本与韩国留学比较| 最近最新免费中文字幕在线| 一个人看视频在线观看www免费| 国产三级在线视频| 午夜福利在线观看吧| 99热精品在线国产| 国产伦精品一区二区三区四那| 久久国产乱子伦精品免费另类| 国产亚洲av嫩草精品影院| 黄色一级大片看看| 免费观看人在逋| 直男gayav资源| www日本黄色视频网| 精品熟女少妇八av免费久了| 少妇高潮的动态图| 亚洲国产精品成人综合色| 欧美中文日本在线观看视频| 色哟哟哟哟哟哟| 一个人免费在线观看电影| 老女人水多毛片| 国产不卡一卡二| 欧美日韩中文字幕国产精品一区二区三区| 窝窝影院91人妻| 丰满人妻一区二区三区视频av| 三级毛片av免费| 两个人的视频大全免费| 69av精品久久久久久| 搞女人的毛片| 国产av不卡久久| 一本精品99久久精品77| 此物有八面人人有两片| 日韩亚洲欧美综合| 久久久久久久久久成人| a级毛片a级免费在线| 中文字幕免费在线视频6| 欧美日本亚洲视频在线播放| 亚洲最大成人手机在线| 全区人妻精品视频| 亚洲av一区综合| 看片在线看免费视频| eeuss影院久久| 日韩高清综合在线| 精品不卡国产一区二区三区| 丰满的人妻完整版| 日日干狠狠操夜夜爽| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美日韩高清在线视频| 国产不卡一卡二| 亚洲成av人片免费观看| 国产精品一区二区性色av| 日本免费a在线| 1024手机看黄色片| 不卡一级毛片| 成人午夜高清在线视频| 99热精品在线国产| 久久精品综合一区二区三区| 别揉我奶头 嗯啊视频| eeuss影院久久| 麻豆成人午夜福利视频| 搞女人的毛片| 我要搜黄色片| 国内精品久久久久精免费| 亚洲精品在线观看二区| 精品人妻熟女av久视频| 婷婷精品国产亚洲av| 欧美精品国产亚洲| 看免费av毛片| 一区福利在线观看| 午夜老司机福利剧场| 国产精品98久久久久久宅男小说| 亚洲最大成人av| 欧美丝袜亚洲另类 | 久久国产精品影院| 午夜老司机福利剧场| 亚洲第一电影网av| 婷婷色综合大香蕉| 在现免费观看毛片| 日韩欧美国产在线观看| 国产精品亚洲一级av第二区| 一二三四社区在线视频社区8| 久久天躁狠狠躁夜夜2o2o| 免费搜索国产男女视频| 欧美国产日韩亚洲一区| 日韩国内少妇激情av| av在线观看视频网站免费| 亚洲黑人精品在线| 老司机午夜福利在线观看视频| 特级一级黄色大片| 亚洲,欧美精品.| 青草久久国产| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 91久久精品电影网| 日韩欧美三级三区| 久久午夜亚洲精品久久| 日本黄色片子视频| 麻豆成人av在线观看| 在线观看一区二区三区| 久久6这里有精品| 国产亚洲欧美98| 亚洲精品456在线播放app | 亚洲av不卡在线观看| 91久久精品电影网| 亚洲一区二区三区色噜噜| 九九热线精品视视频播放| 99在线视频只有这里精品首页| 国产精品电影一区二区三区| 欧美色视频一区免费| 高清日韩中文字幕在线| 高清毛片免费观看视频网站| 1024手机看黄色片| 桃色一区二区三区在线观看| 精品熟女少妇八av免费久了| av专区在线播放| 国产精品久久久久久精品电影| 国产一区二区在线观看日韩| 深爱激情五月婷婷| bbb黄色大片| 午夜日韩欧美国产| 熟妇人妻久久中文字幕3abv| 黄色女人牲交| 最近中文字幕高清免费大全6 | 老熟妇仑乱视频hdxx| 亚洲av成人av| 国产人妻一区二区三区在| 757午夜福利合集在线观看| 婷婷色综合大香蕉| 一个人观看的视频www高清免费观看| 国产欧美日韩精品一区二区| 91字幕亚洲| 亚洲无线在线观看| 精品一区二区三区av网在线观看| 日韩中字成人| 国产成人啪精品午夜网站| 日韩人妻高清精品专区| 99热这里只有是精品在线观看 | 一级作爱视频免费观看| 一个人看的www免费观看视频| a在线观看视频网站| 好看av亚洲va欧美ⅴa在| 在线免费观看的www视频| 日日夜夜操网爽| 日韩欧美一区二区三区在线观看| 又爽又黄无遮挡网站| 日韩中字成人| 国产亚洲欧美98| 人妻夜夜爽99麻豆av| 日韩 亚洲 欧美在线| 又紧又爽又黄一区二区| 亚洲av第一区精品v没综合| 伦理电影大哥的女人| 露出奶头的视频| 国产成+人综合+亚洲专区| 日本与韩国留学比较| 在线观看午夜福利视频| 欧美激情国产日韩精品一区| 亚洲性夜色夜夜综合| 国产高清三级在线| 欧美午夜高清在线| 久久久精品大字幕| 国产精品av视频在线免费观看| 一级黄片播放器| 欧美又色又爽又黄视频| 欧美性感艳星| 赤兔流量卡办理| 亚洲国产高清在线一区二区三| 最近在线观看免费完整版| 国产毛片a区久久久久| 欧美激情在线99| www.熟女人妻精品国产| 久久久久国产精品人妻aⅴ院| 乱码一卡2卡4卡精品| 国产老妇女一区| or卡值多少钱| 91麻豆av在线| 在线国产一区二区在线| 国产精品影院久久| 精品免费久久久久久久清纯| 九色国产91popny在线| 麻豆成人午夜福利视频| 首页视频小说图片口味搜索| 又爽又黄a免费视频| 特级一级黄色大片| 久久精品人妻少妇| 久久这里只有精品中国| 美女xxoo啪啪120秒动态图 | 成年女人永久免费观看视频| 久久午夜福利片| 国产视频内射| 一个人看的www免费观看视频| 欧美黑人欧美精品刺激| 日韩高清综合在线| 脱女人内裤的视频| 亚洲熟妇中文字幕五十中出| 午夜老司机福利剧场| 亚洲av不卡在线观看| 国产午夜福利久久久久久| 亚洲片人在线观看| 舔av片在线| 精品久久久久久久人妻蜜臀av| 日韩欧美精品v在线| 欧美在线一区亚洲| 成人一区二区视频在线观看| 亚洲精品日韩av片在线观看| 国产视频内射| 中文字幕av成人在线电影| a在线观看视频网站| 亚洲专区国产一区二区| 麻豆国产av国片精品| 成年女人看的毛片在线观看| 99久久无色码亚洲精品果冻| 亚洲中文字幕一区二区三区有码在线看| 亚洲精华国产精华精| 午夜日韩欧美国产| 亚洲国产精品成人综合色| 色综合婷婷激情| 欧美黄色片欧美黄色片| 美女被艹到高潮喷水动态| 亚洲自偷自拍三级| 热99re8久久精品国产| 亚洲专区中文字幕在线| av在线天堂中文字幕| 亚洲在线观看片| 国产精品三级大全| 男插女下体视频免费在线播放| 一区二区三区激情视频| 国产69精品久久久久777片| 亚洲美女视频黄频| 欧美国产日韩亚洲一区| 午夜福利欧美成人| 毛片女人毛片| 久久久久免费精品人妻一区二区| 黄色配什么色好看| 欧美一区二区亚洲| 精品不卡国产一区二区三区| 热99re8久久精品国产| 婷婷精品国产亚洲av在线| 亚洲久久久久久中文字幕| 午夜免费成人在线视频| 别揉我奶头 嗯啊视频| 亚洲成人中文字幕在线播放| 日韩大尺度精品在线看网址| 亚洲综合色惰| 一本精品99久久精品77| 99国产综合亚洲精品| 免费av观看视频| 在线观看午夜福利视频| 日韩欧美在线乱码| 日韩欧美精品免费久久 | 欧美xxxx黑人xx丫x性爽| 午夜福利高清视频| 给我免费播放毛片高清在线观看| 亚洲国产精品sss在线观看| 男插女下体视频免费在线播放| 日本与韩国留学比较| 18禁黄网站禁片午夜丰满| 国产精品av视频在线免费观看| 国产精品伦人一区二区| 日韩中字成人| 国产精品野战在线观看| 亚洲无线在线观看| 欧美区成人在线视频| 嫩草影院精品99| 在线天堂最新版资源| 亚洲国产精品成人综合色| 成人午夜高清在线视频| 好看av亚洲va欧美ⅴa在| 欧美性猛交╳xxx乱大交人| 小蜜桃在线观看免费完整版高清| 久久精品人妻少妇| 国产成人影院久久av| 久久久精品大字幕| 熟女电影av网| 深爱激情五月婷婷| 亚洲av成人精品一区久久| 国产三级在线视频| 成人美女网站在线观看视频| 麻豆国产97在线/欧美| 国产精品1区2区在线观看.| 亚洲精品影视一区二区三区av| 久久精品人妻少妇| 色视频www国产| 男女床上黄色一级片免费看| 精品久久久久久成人av| 精品久久久久久久久亚洲 | 高清日韩中文字幕在线| 成人亚洲精品av一区二区| 丁香欧美五月| av专区在线播放| 有码 亚洲区| av天堂中文字幕网| 一夜夜www| 熟妇人妻久久中文字幕3abv| 亚洲国产色片| 在线观看av片永久免费下载| 精品久久久久久久久久久久久| 99热只有精品国产| 日韩精品中文字幕看吧| 91在线精品国自产拍蜜月| 两个人的视频大全免费| 中文字幕精品亚洲无线码一区| 国产精品免费一区二区三区在线| 一级作爱视频免费观看| 又紧又爽又黄一区二区| 亚洲七黄色美女视频| 在线国产一区二区在线| 欧美一区二区国产精品久久精品| 国产精品一区二区免费欧美| 长腿黑丝高跟| 久久国产精品人妻蜜桃| 亚洲av.av天堂| 女人被狂操c到高潮| 12—13女人毛片做爰片一| 国产亚洲精品久久久com| 国产精品人妻久久久久久| 久久久久久久午夜电影| 精品午夜福利视频在线观看一区| 两个人的视频大全免费| 嫩草影院精品99| 别揉我奶头 嗯啊视频| 18+在线观看网站| 久久人人爽人人爽人人片va | 欧美激情在线99| 国产高清三级在线| 婷婷精品国产亚洲av在线|