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

    逆向與正向結(jié)合的導(dǎo)葉式混流泵葉片系統(tǒng)反求設(shè)計(jì)

    2020-07-20 05:14:46唐立新賴喜德陳小明
    關(guān)鍵詞:核主泵葉輪曲面

    杜 江,唐立新,程 海,賴喜德,張 翔,陳小明

    (1.東方法馬通核泵有限責(zé)任公司,四川 德陽(yáng) 618000;2.西華大學(xué)能源與動(dòng)力工程學(xué)院,四川 成都 610039)

    混流泵廣泛用于化工、發(fā)電、水利、市政等行業(yè),特別是核電和核動(dòng)力設(shè)備行業(yè)。核主泵是核電站中最關(guān)鍵的設(shè)備之一,為滿足現(xiàn)有的壓水堆核電站反應(yīng)堆堆芯冷卻循環(huán)的參數(shù)要求,一般都采用導(dǎo)葉式混流泵[1-2]?,F(xiàn)有的和在建的核電站的核主泵都是從國(guó)外進(jìn)口或者國(guó)外公司在國(guó)內(nèi)合資廠家制造。其核心技術(shù)嚴(yán)格控制在國(guó)外幾家公司,特別是水力單元的設(shè)計(jì)核心技術(shù)國(guó)內(nèi)仍未完全掌握。這使得在核主泵制造、維修及運(yùn)行過(guò)程中故障分析等方面都非常困難[3]。如何從已有的產(chǎn)品中獲得準(zhǔn)確的水力單元流道幾何數(shù)據(jù)顯得非常重要。逆向工程技術(shù)是解決該類需求的有效技術(shù)手段[4],已被應(yīng)用于水輪機(jī)葉片、汽輪機(jī)葉片、壓縮機(jī)葉片等反求與三維CAD 模型重構(gòu)[5-8],但是在三維測(cè)量的方便性、葉片曲面重構(gòu)的準(zhǔn)確性、后續(xù)工程應(yīng)用的易用性等方面都有待于提高。這些應(yīng)用多是針對(duì)單個(gè)葉片進(jìn)行三維測(cè)量或者掃描可獲得完整的葉片曲面三維數(shù)據(jù)。

    核電站采用的導(dǎo)葉混流式核主泵的葉輪葉片和空間導(dǎo)葉均為三維雕塑曲面。流道空間狹窄,采用非接觸式光學(xué)三維測(cè)量只能獲取一部分?jǐn)?shù)據(jù),采用接觸式三維測(cè)量更為困難。在導(dǎo)葉式混流泵產(chǎn)品不被破壞的前提下,采用只依靠測(cè)量或掃描的逆向反求技術(shù)很難準(zhǔn)確地獲取葉片系統(tǒng)的曲面幾何數(shù)據(jù)。另外,即使通過(guò)測(cè)量獲取了葉片系統(tǒng)曲面三維數(shù)據(jù)點(diǎn),在沒(méi)有原始設(shè)計(jì)數(shù)據(jù)的情況下如何驗(yàn)證這些數(shù)據(jù)的準(zhǔn)確性也是非常困難的。本文結(jié)合工程實(shí)際需求,以逆向設(shè)計(jì)為基礎(chǔ),采用正向設(shè)計(jì)補(bǔ)充逆向設(shè)計(jì)中信息缺失部分,在獲取葉片系統(tǒng)的不完整曲面幾何數(shù)據(jù)的情況下,根據(jù)有限數(shù)據(jù)點(diǎn)采用葉片數(shù)字化設(shè)計(jì)軟件建立初步葉片系統(tǒng)曲面三維曲面幾何模型,再通過(guò)泵全流道數(shù)值模擬計(jì)算預(yù)測(cè)曲線與原泵外特性曲線的對(duì)比,來(lái)驗(yàn)證反求的葉片系統(tǒng)流道曲面的準(zhǔn)確性,以探索一條基于逆向和正向結(jié)合的導(dǎo)葉式混流泵葉片系統(tǒng)反求設(shè)計(jì)方法及技術(shù)途徑。

    1 混流泵葉片系統(tǒng)的近似測(cè)量

    以某型混流式核主泵為例,其流道包括進(jìn)口管段、閉式葉輪、空間導(dǎo)水機(jī)構(gòu)、類球狀殼體。其進(jìn)口管段和類球狀殼體流道幾何形狀簡(jiǎn)單,可以準(zhǔn)確地測(cè)量到數(shù)據(jù),但其閉式葉輪和空間導(dǎo)水機(jī)構(gòu)的葉片都是雕塑曲面,在不破壞產(chǎn)品的前提下很難準(zhǔn)確地測(cè)量和獲得完整數(shù)據(jù)。通過(guò)現(xiàn)有一些簡(jiǎn)易尺寸測(cè)量方法,可以準(zhǔn)確地獲得葉輪和空間導(dǎo)水機(jī)構(gòu)外輪廓的幾何數(shù)據(jù),但葉輪葉片和導(dǎo)葉葉片的幾何形狀數(shù)據(jù)很難準(zhǔn)確地測(cè)量。參考葉片泵水力設(shè)計(jì)中的水平面葉片剪裁圖繪制思路[4],如圖1 所示,本文采用平面樣板測(cè)繪出葉輪葉片和導(dǎo)葉葉片的等高面截線近似數(shù)據(jù),并在如UG?等三維CAD 軟件中進(jìn)行適當(dāng)?shù)墓饣幚恚垣@得葉輪和空間導(dǎo)水機(jī)構(gòu)的葉片的初步曲面數(shù)據(jù)。該數(shù)據(jù)較為粗糙且與實(shí)際葉片曲面存在偏差,一些曲率較小的區(qū)域偏差較小,曲率變化大的和開(kāi)敞性狹窄的區(qū)域偏差較大。

    圖1 葉片系統(tǒng)測(cè)量的等高面截線

    2 葉片系統(tǒng)的曲面重構(gòu)方法

    葉輪和空間導(dǎo)水機(jī)構(gòu)的葉片都是復(fù)雜雕塑曲面組成的曲面零件,按等高面截線取得初步曲面數(shù)據(jù)不便于修改和控制曲面形狀[9]。為了實(shí)現(xiàn)后續(xù)的反求設(shè)計(jì),需要按葉片曲面體的特征劃分為多個(gè)曲面片來(lái)重構(gòu)。曲面重構(gòu)過(guò)程是根據(jù)實(shí)物模型的幾何拓?fù)湫畔⒅貥?gòu),再現(xiàn)曲面產(chǎn)品特征的過(guò)程。重建曲面片的幾何建模采用NURBS 模型[4]。為了便于基于流場(chǎng)數(shù)值模擬計(jì)算預(yù)測(cè)性能與試驗(yàn)的對(duì)比,以驗(yàn)證反求設(shè)計(jì)過(guò)程中控制葉片形狀,并考慮到葉片水力設(shè)計(jì)過(guò)程是在回轉(zhuǎn)流面上設(shè)計(jì)葉型,所以對(duì)初步曲面數(shù)據(jù)沿回轉(zhuǎn)流面上求交線,再沿回轉(zhuǎn)流面上的交線用蒙皮法構(gòu)造葉片曲面。

    葉片曲面重構(gòu)的主要過(guò)程如下。

    1)將按等高面截線取得的初步曲面數(shù)據(jù)沿等高面截線用三維CAD 軟件(如UG?)按分區(qū)構(gòu)造曲面片,并縫合分別建立葉片正、背面曲面。

    2)利用葉片建模軟件CFturbo?計(jì)算確定回轉(zhuǎn)流面的軸面流線。首先按葉輪前后蓋板線在CFturbo?定義軸面流道,然后根據(jù)軸面流道寬度計(jì)算5~7 條軸面流線,本項(xiàng)目取5 條。

    3)將軸面流線數(shù)據(jù)導(dǎo)入U(xiǎn)G?軟件并繞泵軸心線旋轉(zhuǎn)成回轉(zhuǎn)流面,將葉片正、背面曲面與回轉(zhuǎn)流面求交得到各回轉(zhuǎn)流面上葉片正、背面的交線并進(jìn)行適當(dāng)光滑處理。

    4)將各回轉(zhuǎn)流面上葉片正、背面的交線數(shù)據(jù)導(dǎo)入CFturbo?,在CFturbo?中在回轉(zhuǎn)流面上重建葉型。遵循葉片水力設(shè)計(jì)的相關(guān)準(zhǔn)則和要求,補(bǔ)齊葉片的進(jìn)、出口邊等信息。

    5)對(duì)一些特征數(shù)據(jù)如葉片進(jìn)、出口角,特征點(diǎn)的厚度,與前后蓋板的交點(diǎn)坐標(biāo)尺寸等進(jìn)行校核和必要的調(diào)整和光滑。

    6)在CFturbo?中沿回轉(zhuǎn)流面上重建葉片曲面。

    圖2 為葉片系統(tǒng)的重構(gòu)過(guò)程。葉片的重構(gòu)過(guò)程中根據(jù)葉片水力設(shè)計(jì)方法的相關(guān)要求分析葉片的各流面骨線沿周向保角變換、葉片角變化以及從進(jìn)水邊至出水邊的厚度變化,在CFturbo?中適當(dāng)調(diào)整這些變化曲線,可方便調(diào)整和修改葉片的正、背面曲面。對(duì)一些已知的特征點(diǎn)按偏差控制在小于±0.05 mm 校核調(diào)整后,采用NURBS 模型重建葉片曲面,將重建葉片三維模型導(dǎo)入U(xiǎn)G?。在UG?中將前后蓋板三維模型和葉片三維模型組合生成反求重建的葉片系統(tǒng),如圖3 所示。

    3 基于全流道流場(chǎng)數(shù)值模擬的外特性評(píng)估

    3.1 驗(yàn)證反求設(shè)計(jì)的葉片系統(tǒng)準(zhǔn)確性的策略

    圖2 葉片系統(tǒng)的重構(gòu)過(guò)程

    在沒(méi)有原泵的葉片設(shè)計(jì)數(shù)據(jù)情況下,很難判斷反求設(shè)計(jì)結(jié)果是否與原設(shè)計(jì)一致,但從泵的外特性曲線與廠家提供的試驗(yàn)曲線比較,可以間接驗(yàn)證反求設(shè)計(jì)結(jié)果的準(zhǔn)確性。如果泵的外特性曲線與原泵的試驗(yàn)曲線一致,在0.2~1.2Qd運(yùn)行流量范圍其外特性值與原設(shè)計(jì)試驗(yàn)值的最大偏差在±2%內(nèi),就認(rèn)可反求設(shè)計(jì)的結(jié)果符合工程要求的精度。試驗(yàn)驗(yàn)證不僅成本高,而且周期長(zhǎng),也不便于通過(guò)正向數(shù)值模擬分析修改設(shè)計(jì),因此,本項(xiàng)目采用基于全流道流場(chǎng)數(shù)值模擬預(yù)測(cè)外特性的方法來(lái)預(yù)測(cè)在0.2~1.2Qd運(yùn)行流量范圍的外特性曲線。

    圖3 反求的葉片系統(tǒng)

    3.2 全流道流場(chǎng)數(shù)值模擬

    3.2.1 全流道及離散

    如圖4 所示,將反求得到的葉片系統(tǒng)與進(jìn)口段、殼體和出口延伸段組成該泵的全流道。沿流動(dòng)方向,計(jì)算域依次為進(jìn)口段、葉輪、導(dǎo)葉、殼體和出口延伸段。其中出口延伸段長(zhǎng)度為殼體出口直徑的6 倍。采用CFD 軟件進(jìn)行全流道流場(chǎng)數(shù)值模擬。

    圖4 核主泵全流道計(jì)算域

    對(duì)圖4 所示的整個(gè)流體域的網(wǎng)格用ANSYS/ICEM CFD?進(jìn)行離散,其中進(jìn)口段、葉輪和空間導(dǎo)葉流道內(nèi)的網(wǎng)格為六面體結(jié)構(gòu)化網(wǎng)格,殼體和出口延伸段內(nèi)網(wǎng)格為混合網(wǎng)格。由于整個(gè)流體域中有空間三維扭曲的葉輪和導(dǎo)葉的葉片,網(wǎng)格很難滿足y+都小于20??紤]到本文研究的重點(diǎn)是外特性預(yù)測(cè),而不是邊界層內(nèi)的流動(dòng)規(guī)律模擬,因此,在CFX?中一般要求y+不大于90 即可滿足非邊界層內(nèi)流動(dòng)的計(jì)算要求[10]。

    3.2.2 流動(dòng)模型及求解方法

    在實(shí)際運(yùn)行中,核主泵內(nèi)部為非定常的三維流動(dòng),為了便于采用有限體積法求解,其流動(dòng)控制方程可以寫(xiě)成通用形式[10],為

    式中:φ為描述流場(chǎng)的通用變量,可以表示p、ui、T等物理量;Г為廣義擴(kuò)散系數(shù);S為廣義源項(xiàng)。關(guān)于φ不同,對(duì)應(yīng)于連續(xù)方程、動(dòng)量守恒方程、能量方程等流體狀態(tài)方程參見(jiàn)文獻(xiàn)[10]。

    核主泵內(nèi)部流動(dòng)是復(fù)雜的非線性湍流流動(dòng),再加之邊界形狀復(fù)雜,目前直接求解式(1)非常困難。對(duì)式(1)進(jìn)行雷諾時(shí)均化,采用SSTk-ω湍流模型封閉方程組[10]。SSTk-ω模型的雙方程為:

    式中:Pk表示湍流脈動(dòng)動(dòng)能k的生成項(xiàng);Pω表示湍流脈動(dòng)頻率ω的生成項(xiàng);Γk、Γω表示k和ω的有效擴(kuò)散系數(shù);Yk、Yω表示k和ω的耗散項(xiàng);Dω表示正交擴(kuò)散項(xiàng)。

    SSTk-ω模型是一個(gè)自適應(yīng)湍流模型。其巧妙之處在于它在完全湍流區(qū)運(yùn)用k-ε進(jìn)行求解,在近壁區(qū)運(yùn)用k-ω求解。關(guān)于SST 模型的雙方程中各變量和常數(shù)的詳細(xì)定義參見(jiàn)文獻(xiàn)[10]。對(duì)水力單元過(guò)流部件流道幾何離散化,采用ANSYS/CFX?求解。流體在壁面處的流動(dòng)模式采用無(wú)滑移邊界條件,在CFX?中選用基于k-ω方程的自動(dòng)壁面處理模型[10]。該模型會(huì)根據(jù)網(wǎng)格密度的大小自動(dòng)判斷和選擇處理邊界層內(nèi)的流動(dòng),不需要像標(biāo)準(zhǔn)壁面函數(shù)必須將y+嚴(yán)格控制在小于20 一樣,其y+值只是提供近壁分辨率信息。

    為了驗(yàn)證網(wǎng)格精度和數(shù)值模擬計(jì)算方法的正確性,采用基于流場(chǎng)數(shù)值模擬預(yù)測(cè)出的外特性曲線與該核主泵的冷態(tài)試驗(yàn)曲線比較,來(lái)進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證。其結(jié)果表明,當(dāng)全流道單元總數(shù)達(dá)到280 萬(wàn)后,在設(shè)計(jì)工況下其揚(yáng)程的相對(duì)偏差為±0.6%,效率的相對(duì)偏差為±0.85%,計(jì)算結(jié)果已足以滿足工程精度的需要。在后述流場(chǎng)數(shù)值模擬過(guò)程中,單元數(shù)總數(shù)約285 萬(wàn)。

    3.2.3 邊界條件設(shè)置

    全流道流場(chǎng)數(shù)值模擬為所有過(guò)流部件耦合計(jì)算,葉輪與導(dǎo)葉采用STAGE 模型級(jí)間連接。對(duì)流項(xiàng)和湍流項(xiàng)數(shù)值精度均為2 階。計(jì)算時(shí)按工況以質(zhì)量流量作為進(jìn)口邊界條件。

    4 反求設(shè)計(jì)的葉片曲面精度驗(yàn)證

    4.1 核主泵水力單元主要設(shè)計(jì)參數(shù)

    該泵設(shè)計(jì)工況主要參數(shù)為:體積流量Qd=23 790 m3/h,揚(yáng)程Hd=98 m,轉(zhuǎn)速nd=1 485 r/min,進(jìn)口壓力pin=15.16 MPa,工作溫度T=293 ℃。冷卻劑介質(zhì)為:冷態(tài)密度1 000 kg/m3,冷態(tài)運(yùn)動(dòng)黏度8.93×10-7m2/s;熱態(tài)密度742 kg/m3,熱態(tài)運(yùn)動(dòng)黏度9.42×10-8m2/s。其流道主要幾何參數(shù)為:葉輪出口最大直徑D2max=800 mm,葉輪葉片數(shù)Zb=7,導(dǎo)葉數(shù)Zg=12。制造廠提供了冷態(tài)的試驗(yàn)曲線。

    4.2 外特性預(yù)測(cè)計(jì)算方法

    采用驗(yàn)證了的方法在冷態(tài)工況下對(duì)不同流量的核主泵內(nèi)流場(chǎng)進(jìn)行三維非定常流動(dòng)數(shù)值模擬,得到全流道內(nèi)流場(chǎng)的速度矢量和壓力分布,在此基礎(chǔ)上提取相關(guān)數(shù)據(jù),核主泵外特性預(yù)測(cè)計(jì)算為:

    式中:Pin、Pout分別為泵進(jìn)、出口質(zhì)量加權(quán)平均總壓值,根據(jù)三維流場(chǎng)數(shù)值模擬提取計(jì)算Pin、Pout;M為壓力或黏性力對(duì)葉輪回轉(zhuǎn)軸的力矩,葉輪扭矩,N·m,其計(jì)算方法參見(jiàn)文獻(xiàn)[10];ω為葉輪的角速度,rad/s;ρ為流體密度,kg/m3;g為重力加速度。

    4.3 預(yù)測(cè)的外特性與試驗(yàn)曲線對(duì)比

    根據(jù)該泵的設(shè)計(jì)參數(shù),參考制造廠提供的冷態(tài)試驗(yàn)曲線的流量范圍,采用3.2 節(jié)中所述的全流道流場(chǎng)數(shù)值模擬方法,對(duì)冷態(tài)下0.2~1.2Qd流量范圍,每隔0.2Qd進(jìn)行三維穩(wěn)態(tài)流場(chǎng)數(shù)值模擬,共6 個(gè)工況。根據(jù)數(shù)值模擬結(jié)果按4.2 節(jié)中預(yù)測(cè)計(jì)算揚(yáng)程、效率,將預(yù)測(cè)曲線與制造廠提供的冷態(tài)試驗(yàn)曲線對(duì)比,其結(jié)果如圖5 所示??梢钥闯?,預(yù)測(cè)結(jié)果與試驗(yàn)具有較高的吻合程度,預(yù)測(cè)的揚(yáng)程和效率變化趨勢(shì)幾乎與試驗(yàn)結(jié)果一致。設(shè)計(jì)流量下,預(yù)測(cè)計(jì)算揚(yáng)程為96.94 m,與原泵設(shè)計(jì)揚(yáng)程(Hd=98 m)吻合程度很高。在0.2~1.2Qd流量范圍,其揚(yáng)程的最大相對(duì)偏差為1.38%,效率的最大相對(duì)偏差為1.41%,主要是在小流量工況下偏差較大,這是因?yàn)樵谛×髁抗r下采用三維穩(wěn)態(tài)流場(chǎng)數(shù)值計(jì)算不能準(zhǔn)確模擬泵內(nèi)流場(chǎng)。

    圖5 水力單元預(yù)測(cè)與試驗(yàn)外特性對(duì)比

    在0.2~1.2Qd范圍預(yù)測(cè)的外特性曲線與原設(shè)計(jì)的揚(yáng)程、效率試驗(yàn)值的最大偏差都在±2%以內(nèi),在設(shè)計(jì)工況下其揚(yáng)程的相對(duì)偏差為0.6%,效率的相對(duì)偏差為0.85%。從工程的觀點(diǎn)來(lái)看,反求設(shè)計(jì)的結(jié)果已滿足工程要求的精度。

    5 基于正向流場(chǎng)分析修改設(shè)計(jì)

    根據(jù)前述的驗(yàn)證反求設(shè)計(jì)葉片系統(tǒng)的策略,如果在0.2~1.2Qd范圍預(yù)測(cè)的外特性曲線與原設(shè)計(jì)的揚(yáng)程、效率試驗(yàn)值的最大偏差超過(guò)±2%,應(yīng)該按第2 節(jié)中所述方法在CFturbo?中進(jìn)一步修改。在修改設(shè)計(jì)過(guò)程中,分析葉片系統(tǒng)中如圖6 所示各流面上的流速分布、壓力分布,根據(jù)泵葉片設(shè)計(jì)方法在CFturbo?中通過(guò)修改各流面上的葉片角變化來(lái)適當(dāng)修改葉片曲面;再將修改的葉片系統(tǒng)三維幾何模型導(dǎo)入CFD 軟件,基于正向設(shè)計(jì)分析思路,采用第4 節(jié)中的流場(chǎng)數(shù)值模擬方法進(jìn)行三維穩(wěn)態(tài)流場(chǎng)數(shù)值模擬和外特性預(yù)測(cè),按4.3 節(jié)中所述的過(guò)程進(jìn)行反求設(shè)計(jì)的葉片系統(tǒng)準(zhǔn)確性評(píng)估。如果反求設(shè)計(jì)葉片系統(tǒng)的外特性預(yù)測(cè)結(jié)果沒(méi)有滿足驗(yàn)證策略的偏差控制要求,再按前述過(guò)程迭代修改。

    6 結(jié)論

    本文以某型混流式核主泵為研究對(duì)象,采用以逆向設(shè)計(jì)為基礎(chǔ),正向設(shè)計(jì)補(bǔ)充逆向設(shè)計(jì)中信息缺失部分的技術(shù)途徑來(lái)探索產(chǎn)品不被破壞的前提下如何反求導(dǎo)葉式混流泵的葉片系統(tǒng)曲面幾何數(shù)據(jù)方法。

    圖6 數(shù)值模擬中提取流面上的流場(chǎng)分布

    1)基于文中提出的逆向與正向相結(jié)合的葉片反求設(shè)計(jì)方法,在產(chǎn)品不被破壞的前提下,可以較好地解決具有復(fù)雜曲面的混流泵在不能獲取葉片系統(tǒng)的完整曲面幾何數(shù)據(jù)的情況下的流道反求設(shè)計(jì)問(wèn)題。

    2)采用基于全流道流場(chǎng)數(shù)值模擬預(yù)測(cè)泵外特性結(jié)果與原廠家提供的試驗(yàn)曲線比較的方法來(lái)評(píng)估過(guò)流部件反求設(shè)計(jì)的結(jié)果,這種替代試驗(yàn)評(píng)估的技術(shù)路線是可行的,能夠滿足工程精度要求。

    3)文中提出的葉片反求設(shè)計(jì)的方法可推廣用于其他葉片式流體機(jī)械的葉片系統(tǒng)反求設(shè)計(jì)。

    猜你喜歡
    核主泵葉輪曲面
    1.4317 QT2鋼在高能泵葉輪上的應(yīng)用
    相交移動(dòng)超曲面的亞純映射的唯一性
    圓環(huán)上的覆蓋曲面不等式及其應(yīng)用
    應(yīng)用石膏型快速精密鑄造技術(shù)制造葉輪
    核主泵水潤(rùn)滑推力軸承不銹鋼螺紋鎖死分析及改進(jìn)
    誰(shuí)在分割中國(guó)核主泵巨型蛋糕?
    能源(2016年1期)2016-12-01 05:10:05
    基于曲面展開(kāi)的自由曲面網(wǎng)格劃分
    離心泵葉輪切割方法
    核主泵端面密封類型的可靠性分析
    基于CFD/CSD耦合的葉輪機(jī)葉片失速顫振計(jì)算
    纯流量卡能插随身wifi吗| 日韩欧美三级三区| 久久精品亚洲精品国产色婷小说| 后天国语完整版免费观看| or卡值多少钱| 国产一区二区三区在线臀色熟女| 欧美激情极品国产一区二区三区| 一区二区日韩欧美中文字幕| 亚洲全国av大片| 一级a爱视频在线免费观看| 好看av亚洲va欧美ⅴa在| 真人做人爱边吃奶动态| 国内久久婷婷六月综合欲色啪| 男人操女人黄网站| 一级片免费观看大全| 国产精品av久久久久免费| 大陆偷拍与自拍| 欧美日本中文国产一区发布| 老熟妇仑乱视频hdxx| 国产伦一二天堂av在线观看| 女人精品久久久久毛片| 亚洲黑人精品在线| 亚洲熟妇中文字幕五十中出| 久久久久久人人人人人| 精品国产国语对白av| 精品久久久久久成人av| 桃红色精品国产亚洲av| 国产欧美日韩精品亚洲av| 午夜免费成人在线视频| 亚洲专区国产一区二区| 国产精品一区二区免费欧美| svipshipincom国产片| 69精品国产乱码久久久| 亚洲狠狠婷婷综合久久图片| 欧美黄色片欧美黄色片| 午夜福利高清视频| 亚洲最大成人中文| 日韩大尺度精品在线看网址 | 久久人妻av系列| 婷婷丁香在线五月| 午夜老司机福利片| 久9热在线精品视频| 电影成人av| 香蕉丝袜av| 欧美日韩中文字幕国产精品一区二区三区 | 大码成人一级视频| 日韩有码中文字幕| 欧美激情高清一区二区三区| 欧美老熟妇乱子伦牲交| 国产精品亚洲一级av第二区| 欧美激情久久久久久爽电影 | 精品国产一区二区三区四区第35| 亚洲 国产 在线| 妹子高潮喷水视频| 99国产精品99久久久久| 最新美女视频免费是黄的| 日韩一卡2卡3卡4卡2021年| 一区二区三区精品91| 99国产综合亚洲精品| 99热只有精品国产| 国产av一区二区精品久久| 日韩精品青青久久久久久| 动漫黄色视频在线观看| 777久久人妻少妇嫩草av网站| 成人精品一区二区免费| 免费一级毛片在线播放高清视频 | 久久精品91蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 9热在线视频观看99| 欧美精品亚洲一区二区| 欧美黄色淫秽网站| 一边摸一边抽搐一进一小说| 国产精品久久久久久人妻精品电影| 亚洲全国av大片| 久久久久国产精品人妻aⅴ院| 久久久久精品国产欧美久久久| 亚洲avbb在线观看| 成人三级做爰电影| 性少妇av在线| 午夜免费激情av| 亚洲aⅴ乱码一区二区在线播放 | 精品久久蜜臀av无| 中文字幕人妻熟女乱码| 久久精品91无色码中文字幕| 这个男人来自地球电影免费观看| 最近最新中文字幕大全免费视频| 免费人成视频x8x8入口观看| 极品教师在线免费播放| 免费观看人在逋| 亚洲精品美女久久av网站| 日韩欧美一区二区三区在线观看| 国产三级在线视频| 99国产极品粉嫩在线观看| 亚洲熟女毛片儿| 亚洲国产精品合色在线| 咕卡用的链子| 日韩中文字幕欧美一区二区| 久久久久久人人人人人| 一进一出抽搐gif免费好疼| 免费无遮挡裸体视频| 岛国在线观看网站| 久热爱精品视频在线9| 妹子高潮喷水视频| 亚洲美女黄片视频| 91字幕亚洲| 色哟哟哟哟哟哟| 精品不卡国产一区二区三区| 99国产精品一区二区三区| 丝袜美足系列| 可以在线观看毛片的网站| 国产高清videossex| 日韩欧美在线二视频| 十分钟在线观看高清视频www| 精品国产一区二区久久| www.www免费av| 色播亚洲综合网| 高清黄色对白视频在线免费看| 俄罗斯特黄特色一大片| 欧美日韩中文字幕国产精品一区二区三区 | 99国产精品一区二区三区| 好看av亚洲va欧美ⅴa在| 亚洲av片天天在线观看| 人人澡人人妻人| 伦理电影免费视频| 国产激情久久老熟女| 亚洲第一电影网av| 色av中文字幕| 久久久精品国产亚洲av高清涩受| 一区二区三区精品91| 色精品久久人妻99蜜桃| 成人手机av| 桃红色精品国产亚洲av| www.www免费av| 免费看十八禁软件| 久久精品91无色码中文字幕| 在线观看午夜福利视频| 久久亚洲真实| 18美女黄网站色大片免费观看| 久久天躁狠狠躁夜夜2o2o| 啪啪无遮挡十八禁网站| 99国产综合亚洲精品| or卡值多少钱| 成在线人永久免费视频| 黄片小视频在线播放| 久久中文看片网| 欧美黄色片欧美黄色片| 99精品欧美一区二区三区四区| 亚洲成人国产一区在线观看| 午夜福利一区二区在线看| 久久狼人影院| 一级a爱片免费观看的视频| 搡老岳熟女国产| 美女免费视频网站| 国产在线观看jvid| 97人妻精品一区二区三区麻豆 | 熟女少妇亚洲综合色aaa.| 性欧美人与动物交配| 十分钟在线观看高清视频www| 极品人妻少妇av视频| 欧美色欧美亚洲另类二区 | 亚洲人成网站在线播放欧美日韩| 老熟妇乱子伦视频在线观看| 亚洲最大成人中文| 色婷婷久久久亚洲欧美| 欧美日韩精品网址| 亚洲专区国产一区二区| 夜夜看夜夜爽夜夜摸| 亚洲精品在线美女| 不卡一级毛片| 12—13女人毛片做爰片一| 午夜福利高清视频| 亚洲全国av大片| 色在线成人网| 免费av毛片视频| 国产一区二区激情短视频| 黑丝袜美女国产一区| а√天堂www在线а√下载| 国产精品久久视频播放| 可以在线观看的亚洲视频| 国产区一区二久久| 一卡2卡三卡四卡精品乱码亚洲| 欧美最黄视频在线播放免费| 中文字幕av电影在线播放| 真人做人爱边吃奶动态| 99国产精品99久久久久| 精品国产乱子伦一区二区三区| 久久青草综合色| 香蕉久久夜色| 亚洲专区国产一区二区| 亚洲一区二区三区色噜噜| 女人被躁到高潮嗷嗷叫费观| 少妇粗大呻吟视频| 99re在线观看精品视频| 亚洲男人天堂网一区| 91av网站免费观看| 黄频高清免费视频| 亚洲欧美一区二区三区黑人| 人成视频在线观看免费观看| av视频在线观看入口| 日日干狠狠操夜夜爽| 国产熟女午夜一区二区三区| av福利片在线| 国产精品亚洲av一区麻豆| 人妻久久中文字幕网| 久久热在线av| 热re99久久国产66热| 一级a爱视频在线免费观看| 国产一区二区三区视频了| 国产一区二区三区视频了| 波多野结衣高清无吗| 亚洲精品久久成人aⅴ小说| 亚洲精品美女久久久久99蜜臀| 欧美黑人精品巨大| 国产精品自产拍在线观看55亚洲| 黄频高清免费视频| 亚洲精华国产精华精| 欧美色视频一区免费| 亚洲男人的天堂狠狠| 日韩大尺度精品在线看网址 | 一区二区三区高清视频在线| 此物有八面人人有两片| 久久天堂一区二区三区四区| 天天添夜夜摸| 久久精品91无色码中文字幕| 亚洲五月婷婷丁香| 欧美黄色淫秽网站| 50天的宝宝边吃奶边哭怎么回事| 热re99久久国产66热| 老司机在亚洲福利影院| 国产成人欧美| 欧美日本亚洲视频在线播放| 久久精品国产清高在天天线| 国产黄a三级三级三级人| 免费人成视频x8x8入口观看| 亚洲人成电影免费在线| 国内久久婷婷六月综合欲色啪| 淫妇啪啪啪对白视频| 91在线观看av| 国产亚洲av高清不卡| 午夜精品国产一区二区电影| av天堂在线播放| 国产成年人精品一区二区| 免费看美女性在线毛片视频| 一区二区三区精品91| 免费在线观看亚洲国产| 免费久久久久久久精品成人欧美视频| 国产亚洲精品av在线| 亚洲熟女毛片儿| 久久国产乱子伦精品免费另类| 91精品三级在线观看| 老鸭窝网址在线观看| 两人在一起打扑克的视频| 欧美久久黑人一区二区| 无限看片的www在线观看| 亚洲熟女毛片儿| 91成人精品电影| 国产私拍福利视频在线观看| 一个人免费在线观看的高清视频| 老鸭窝网址在线观看| 人妻久久中文字幕网| 久久精品亚洲熟妇少妇任你| 欧美+亚洲+日韩+国产| 18禁国产床啪视频网站| 69av精品久久久久久| 精品一区二区三区四区五区乱码| 黄片播放在线免费| 亚洲色图av天堂| 免费看a级黄色片| 亚洲精品美女久久久久99蜜臀| 久久精品国产亚洲av高清一级| 日韩大码丰满熟妇| 免费在线观看亚洲国产| 人妻丰满熟妇av一区二区三区| 亚洲一区中文字幕在线| 色在线成人网| 免费看美女性在线毛片视频| 真人一进一出gif抽搐免费| 多毛熟女@视频| 亚洲午夜理论影院| 国产成人欧美| 精品国产美女av久久久久小说| 91av网站免费观看| 无人区码免费观看不卡| 一级毛片女人18水好多| 欧美久久黑人一区二区| 久久青草综合色| av免费在线观看网站| 亚洲精品一卡2卡三卡4卡5卡| 中文亚洲av片在线观看爽| 精品福利观看| 国产99白浆流出| 国产在线观看jvid| 黄色视频不卡| 老鸭窝网址在线观看| 美女 人体艺术 gogo| 亚洲成a人片在线一区二区| 精品无人区乱码1区二区| 伦理电影免费视频| 国产伦一二天堂av在线观看| 久久精品91蜜桃| 国产xxxxx性猛交| 国产一区二区在线av高清观看| 老鸭窝网址在线观看| 国产麻豆69| 午夜久久久久精精品| 久久香蕉精品热| 男女床上黄色一级片免费看| 成人免费观看视频高清| 1024视频免费在线观看| 欧美黄色淫秽网站| 脱女人内裤的视频| 国产亚洲av高清不卡| 成人三级黄色视频| 亚洲av第一区精品v没综合| 午夜老司机福利片| 久久人妻av系列| 国产精品1区2区在线观看.| 免费少妇av软件| 久久久精品欧美日韩精品| 69精品国产乱码久久久| 我的亚洲天堂| 久久国产精品影院| 午夜a级毛片| 欧美在线一区亚洲| www.熟女人妻精品国产| 久久精品亚洲熟妇少妇任你| 中文字幕人成人乱码亚洲影| 99国产综合亚洲精品| 无人区码免费观看不卡| 真人做人爱边吃奶动态| 午夜视频精品福利| 19禁男女啪啪无遮挡网站| 十八禁人妻一区二区| 午夜福利成人在线免费观看| 天天一区二区日本电影三级 | 国产伦人伦偷精品视频| 99久久国产精品久久久| 一区二区三区国产精品乱码| 日韩精品青青久久久久久| 女性生殖器流出的白浆| 久久九九热精品免费| 法律面前人人平等表现在哪些方面| 国产免费男女视频| 人人妻,人人澡人人爽秒播| 天堂动漫精品| 久久午夜亚洲精品久久| 91大片在线观看| 国产亚洲精品久久久久久毛片| 日韩欧美三级三区| 欧美成人一区二区免费高清观看 | av片东京热男人的天堂| 精品久久久久久久人妻蜜臀av | 18禁国产床啪视频网站| 欧美乱妇无乱码| 一级毛片女人18水好多| 亚洲国产欧美日韩在线播放| 亚洲av日韩精品久久久久久密| 成年女人毛片免费观看观看9| 午夜激情av网站| 亚洲久久久国产精品| 一进一出抽搐动态| 国产激情久久老熟女| 少妇 在线观看| 在线免费观看的www视频| 亚洲人成电影观看| 人人澡人人妻人| 999精品在线视频| 亚洲av电影在线进入| 精品日产1卡2卡| 999精品在线视频| 午夜亚洲福利在线播放| 国产人伦9x9x在线观看| 国产三级在线视频| 亚洲五月婷婷丁香| 丰满人妻熟妇乱又伦精品不卡| 一级a爱视频在线免费观看| www.999成人在线观看| 国产精品亚洲av一区麻豆| 欧美日韩精品网址| 国产欧美日韩一区二区三| 18禁美女被吸乳视频| 国产熟女午夜一区二区三区| 高清在线国产一区| 欧美精品啪啪一区二区三区| www.自偷自拍.com| АⅤ资源中文在线天堂| 久久午夜亚洲精品久久| 亚洲天堂国产精品一区在线| 老熟妇乱子伦视频在线观看| 美女大奶头视频| 亚洲av日韩精品久久久久久密| 麻豆国产av国片精品| 黑人欧美特级aaaaaa片| 精品少妇一区二区三区视频日本电影| 黄色视频不卡| 亚洲欧洲精品一区二区精品久久久| www.自偷自拍.com| 欧美成狂野欧美在线观看| 99国产精品一区二区三区| 十八禁网站免费在线| 成人三级做爰电影| 免费看a级黄色片| 香蕉国产在线看| 亚洲国产精品合色在线| 香蕉丝袜av| 日韩三级视频一区二区三区| 中文字幕久久专区| 淫秽高清视频在线观看| 亚洲欧洲精品一区二区精品久久久| 国产私拍福利视频在线观看| 人人澡人人妻人| 国产一区二区三区综合在线观看| 亚洲欧美日韩高清在线视频| 亚洲精品粉嫩美女一区| 咕卡用的链子| cao死你这个sao货| 又紧又爽又黄一区二区| 亚洲色图 男人天堂 中文字幕| 12—13女人毛片做爰片一| av福利片在线| 91成人精品电影| 欧美成人免费av一区二区三区| 精品一品国产午夜福利视频| 亚洲国产精品久久男人天堂| 久久中文看片网| 美女国产高潮福利片在线看| 黄频高清免费视频| 国产精品1区2区在线观看.| 一边摸一边抽搐一进一小说| 久久久久久免费高清国产稀缺| 大香蕉久久成人网| 午夜视频精品福利| 精品无人区乱码1区二区| 久久人人精品亚洲av| 91麻豆av在线| 国产一卡二卡三卡精品| 欧美最黄视频在线播放免费| 99精品欧美一区二区三区四区| 精品国内亚洲2022精品成人| 精品久久久久久久毛片微露脸| 国产熟女xx| 最近最新中文字幕大全免费视频| 搞女人的毛片| 国产黄a三级三级三级人| 久久久久精品国产欧美久久久| 色在线成人网| 国产精品美女特级片免费视频播放器 | av欧美777| 一级作爱视频免费观看| 精品一区二区三区视频在线观看免费| 乱人伦中国视频| 欧美午夜高清在线| 亚洲第一欧美日韩一区二区三区| 91字幕亚洲| 精品欧美国产一区二区三| 久久香蕉精品热| 午夜福利在线观看吧| 久久国产精品影院| 人成视频在线观看免费观看| 无限看片的www在线观看| 天堂√8在线中文| 成人18禁高潮啪啪吃奶动态图| 黑人巨大精品欧美一区二区蜜桃| 日韩视频一区二区在线观看| 啦啦啦韩国在线观看视频| 黄片小视频在线播放| 国产片内射在线| 日本 av在线| 久久久久精品国产欧美久久久| 婷婷六月久久综合丁香| 国产精品永久免费网站| 久久精品国产亚洲av香蕉五月| 成人18禁在线播放| 久久狼人影院| 国产欧美日韩综合在线一区二区| 亚洲 欧美 日韩 在线 免费| 夜夜爽天天搞| 757午夜福利合集在线观看| 久久中文字幕一级| 色播在线永久视频| 久久久久久久午夜电影| 欧美最黄视频在线播放免费| 亚洲成av片中文字幕在线观看| 成人18禁高潮啪啪吃奶动态图| 两个人免费观看高清视频| 美女扒开内裤让男人捅视频| 欧美av亚洲av综合av国产av| 日本撒尿小便嘘嘘汇集6| 此物有八面人人有两片| 欧美日韩一级在线毛片| 麻豆久久精品国产亚洲av| 男人的好看免费观看在线视频 | 亚洲人成电影观看| 日韩精品青青久久久久久| 欧美成人免费av一区二区三区| 亚洲第一青青草原| 久久精品国产综合久久久| 人成视频在线观看免费观看| 悠悠久久av| 十分钟在线观看高清视频www| 日韩欧美一区视频在线观看| 亚洲色图综合在线观看| 久久人妻熟女aⅴ| 女同久久另类99精品国产91| 黄色片一级片一级黄色片| 啦啦啦免费观看视频1| 亚洲情色 制服丝袜| 久久久久久大精品| av免费在线观看网站| 欧美日韩一级在线毛片| 国产又色又爽无遮挡免费看| 精品国内亚洲2022精品成人| 国产高清有码在线观看视频 | 精品久久久久久久久久免费视频| 成在线人永久免费视频| 亚洲五月色婷婷综合| 精品国产一区二区久久| 男女做爰动态图高潮gif福利片 | 国产欧美日韩一区二区三| 91精品国产国语对白视频| 69av精品久久久久久| 亚洲成人久久性| 中文字幕色久视频| 人人妻人人爽人人添夜夜欢视频| 成人永久免费在线观看视频| 桃红色精品国产亚洲av| 人人妻人人澡欧美一区二区 | 亚洲中文av在线| 岛国在线观看网站| 午夜精品久久久久久毛片777| 97人妻精品一区二区三区麻豆 | 精品国产国语对白av| 色精品久久人妻99蜜桃| 此物有八面人人有两片| 久久久久精品国产欧美久久久| 久久久久久大精品| 亚洲欧美日韩另类电影网站| 国产极品粉嫩免费观看在线| 男女下面进入的视频免费午夜 | 超碰成人久久| 后天国语完整版免费观看| 午夜福利18| 9色porny在线观看| 香蕉丝袜av| 精品国产超薄肉色丝袜足j| 成人国语在线视频| 亚洲久久久国产精品| 又黄又粗又硬又大视频| 青草久久国产| 丁香欧美五月| 啦啦啦韩国在线观看视频| 亚洲aⅴ乱码一区二区在线播放 | 女同久久另类99精品国产91| 日本vs欧美在线观看视频| avwww免费| 黄色毛片三级朝国网站| 亚洲精华国产精华精| 国产精品影院久久| 最近最新中文字幕大全电影3 | 一个人观看的视频www高清免费观看 | 深夜精品福利| 欧美性长视频在线观看| 欧美色视频一区免费| 午夜福利欧美成人| 99国产精品99久久久久| 久久久久国内视频| 亚洲美女黄片视频| 亚洲欧美日韩高清在线视频| 99国产精品一区二区蜜桃av| ponron亚洲| 丰满人妻熟妇乱又伦精品不卡| 天天添夜夜摸| 国产99白浆流出| 老司机福利观看| 久久久久久国产a免费观看| 国语自产精品视频在线第100页| 久久久精品国产亚洲av高清涩受| 欧美午夜高清在线| 国产免费男女视频| 国产熟女午夜一区二区三区| 日本 av在线| 国产亚洲精品综合一区在线观看 | а√天堂www在线а√下载| 国产高清视频在线播放一区| 久久伊人香网站| 99国产综合亚洲精品| 国产一级毛片七仙女欲春2 | 国产精品一区二区精品视频观看| 免费在线观看亚洲国产| 午夜福利,免费看| 色婷婷久久久亚洲欧美| 人人妻人人澡人人看| 久久天堂一区二区三区四区| 国产免费男女视频| 最近最新中文字幕大全电影3 | www.www免费av| 亚洲成人精品中文字幕电影| 亚洲 欧美一区二区三区| 亚洲五月天丁香| 久久久国产欧美日韩av| 夜夜躁狠狠躁天天躁| 母亲3免费完整高清在线观看| 嫩草影视91久久| 国产日韩一区二区三区精品不卡| а√天堂www在线а√下载| 亚洲 国产 在线| 国产aⅴ精品一区二区三区波| а√天堂www在线а√下载| 久久人人97超碰香蕉20202| 日韩欧美国产在线观看| 亚洲av第一区精品v没综合| 男人舔女人的私密视频| 久久国产精品人妻蜜桃| 两个人视频免费观看高清| 波多野结衣高清无吗|