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

    隨機風(fēng)沙沖蝕葉片涂層的數(shù)值模擬研究

    2023-05-30 14:45:42劉士毅岑海堂劉光友
    機械科學(xué)與技術(shù) 2023年5期
    關(guān)鍵詞:沙粒沖蝕風(fēng)沙

    劉士毅,岑海堂,劉光友

    (內(nèi)蒙古工業(yè)大學(xué) 機械工程學(xué)院,呼和浩特 010051)

    中國幅員遼闊,但其地貌復(fù)雜、各地區(qū)氣候存在很大差異。西北華北地區(qū)沙化嚴(yán)重,并形成中國八大沙漠[1-3]、東北冰雪多、東南沿海地區(qū)暴雨頻發(fā),然而中國的風(fēng)能資源又主要聚積于此,所以風(fēng)力機的工作環(huán)境差,常年遭受風(fēng)沙沖蝕、陽光曝曬以及極端溫差作用[4],導(dǎo)致葉片涂層損傷嚴(yán)重、修復(fù)困難、維修成本高,嚴(yán)重影響葉片的氣動性能以及風(fēng)力機組的輸出功率,極大地降低了風(fēng)力機的經(jīng)濟性[5-7]。張永等[8]基于自制試驗臺,采用氣流挾沙噴射法研究了風(fēng)沙作用下風(fēng)力機葉片涂層沖蝕特性,進而得到涂層沖蝕磨損量與沖蝕速度的關(guān)系,并建立了涂層沖蝕磨損程度的評價公式。徐偉勝等[9]對航空發(fā)動機葉片表面硬質(zhì)涂層受沙粒高角度沖擊而出現(xiàn)裂紋的問題進行了研究,通過實驗?zāi)M硬質(zhì)顆粒以高角度重復(fù)沖擊TiN/Ti 硬質(zhì)涂層,并研究了調(diào)制比、層數(shù)結(jié)構(gòu)對硬質(zhì)涂層沖擊損傷的影響,結(jié)果表明硬質(zhì)顆粒重復(fù)沖擊作用引起涂層疲勞剝落和圓周疲勞裂紋的產(chǎn)生。戴麗萍等[10]針對含沙氣流對風(fēng)力機葉片沖蝕作用問題進行了數(shù)值模擬仿真,并基于實驗數(shù)據(jù)建立了玻璃鋼材料的沙蝕沖蝕率模型,進而基于風(fēng)力發(fā)電機葉輪三維流場及沙粒運動特點,計算得到了葉片表面在各工況下的沖蝕率,結(jié)果表明葉片中、外葉展前緣區(qū)域沖蝕損傷最為嚴(yán)重。本文基于已有試驗結(jié)果,利用EDEM 與Fluent 耦合進行風(fēng)沙沖蝕進行數(shù)值模擬,得到隨機載荷數(shù)據(jù),并對葉片涂層進行有限元分析,為涂層疲勞強度計算和后期維護提供理論基礎(chǔ)。

    1 風(fēng)機葉片涂層三維模型

    本文參照3 MW 風(fēng)力機結(jié)合Glauert 法確定風(fēng)機主要參數(shù),通過氣動分析軟件Profili 和三維建模軟件UG 建立涂層三維模型。

    該風(fēng)機葉片模型中,風(fēng)輪直徑D= 90 m,取輪轂半徑rhub= 1.5 m,則葉片長度L= 43.5 m,風(fēng)輪轉(zhuǎn)速n= 19.3 r/min,沿葉片展向選取6 個翼型截面,其占位分別為0.17、0.42、0.57、0.75、0.95 和1。本文結(jié)合風(fēng)力機葉片專用翼型的相關(guān)特性選用FX77-343、S818、S830 及S831 這4 種翼型,并由氣動分析軟件Profili 計算得到各翼型在最大升阻比時的攻角和升力系數(shù),并通過葛勞渥設(shè)計法的相關(guān)公式計算得到葉片各翼型的攻角αi、入流角φi、安裝角βi及弦長ci等建模參數(shù),如表1 所示。

    表1 葉片建模參數(shù)

    通過Profili 生成上述各類型翼型輪廓數(shù)據(jù),將數(shù)據(jù)導(dǎo)入UG 建立各翼型曲線串,建立的葉片三維模型,如圖1 所示;葉片表面建立的涂層三維模型,厚度為2 mm,如圖2 所示。

    圖1 葉片三維模型

    圖2 涂層三維模型

    2 隨機風(fēng)速模型建立

    本文以內(nèi)蒙古地區(qū)為例,對大型風(fēng)電場風(fēng)力進行風(fēng)速數(shù)據(jù)監(jiān)測如圖3 所示。得到實測風(fēng)速樣本總均值和總方差分別為12.34 和5.24。

    圖3 內(nèi)蒙古某大型風(fēng)電場實測風(fēng)速

    自然風(fēng)受大氣壓強的影響而瞬息萬變,導(dǎo)致風(fēng)速的變化具有明顯的隨機性,工程實踐表明威布爾分布是目前最能描述自然風(fēng)速的數(shù)學(xué)模型,其結(jié)果與風(fēng)場實際監(jiān)測數(shù)據(jù)較為吻合,被廣泛應(yīng)用于風(fēng)能相關(guān)領(lǐng)域的計算中。

    求解方程為:

    式中:V為風(fēng)速;U為[0,1]之間的隨機數(shù)。

    對于給定樣本的均值和方差,求解威布爾分布的形狀參數(shù)k和尺度參數(shù)c[11],即

    最后可得到基于威布爾分布模型的風(fēng)速數(shù)學(xué)為

    通過MATLAB 自帶隨機函數(shù)可得威布爾風(fēng)速模型概率密度如圖4 所示,威布爾風(fēng)速模型的分布如圖5所示,威布爾風(fēng)速模型的隨機風(fēng)速時間序列如圖6 所示。

    圖4 威布爾風(fēng)速模型的概率密度圖

    圖5 威布爾風(fēng)速模型的分布圖

    圖6 威布爾風(fēng)速模型的時間序列圖

    3 隨機風(fēng)沙沖蝕仿真

    由于自然風(fēng)速不可避免的呈威布爾分布隨機變化,因此對于內(nèi)蒙古地區(qū)服務(wù)于沙塵環(huán)境下的風(fēng)力機葉片涂層耐久性的研究,應(yīng)充分考慮變量(風(fēng)速大小、沙粒直徑等)的隨機性,進而生成葉片涂層及流體域的三維網(wǎng)格模型、離散元分析軟件EDEM 計算沙粒數(shù)據(jù)以及流體分析軟件Fluent 計算流場數(shù)據(jù),得到隨機風(fēng)沙沖蝕葉片涂層的隨機載荷數(shù)據(jù)。

    3.1 EDEM 和Fluent 耦合仿真的基本流程

    Fluent 自帶的離散項也能進行數(shù)值模擬仿真,但是其限制離散項體積分?jǐn)?shù)不超過10%,只能進行稀疏的低濃度兩項流的仿真,這為工程實例的仿真帶來極大的不便,因此本文通過Fluent 與EDEM 耦合仿真,以此實現(xiàn)對風(fēng)沙兩項流沖蝕葉片涂層的模擬仿真。基于拉格朗日耦合仿真的基本思路可得到EDEM 與Fluent 耦合仿真的基本流程,如圖7 所示。

    圖7 EDEM 與Fluent 耦合仿真流程

    3.2 網(wǎng)格劃分

    創(chuàng)建計算模型流體域,設(shè)置流體域的入口邊界inlet、出口邊界outlet、壁面wall 以及流固耦合面FSI,設(shè)置旋轉(zhuǎn)的流體進行旋轉(zhuǎn)速度與方向風(fēng)機葉片的旋轉(zhuǎn)方向為(0,0,1),旋轉(zhuǎn)速度為19.3 r/min。建立流體域四面體網(wǎng)格,其節(jié)點為270 136,單元數(shù)為1 154 386。網(wǎng)格劃分后的流體域模型如圖8 所示。

    圖8 流體網(wǎng)格劃分

    3.3 隨機風(fēng)沙沖蝕仿真

    本文以內(nèi)蒙古某地為例,其沙粒直徑分布及其含量分布如表2 所示。為了最大限度貼合實際工況,本次模擬仿真分析中采用的沙粒直徑數(shù)據(jù)直接參照表2 選取,結(jié)果如表3 所示。

    表2 庫布齊沙漠沙粒粒徑分布及其含量分布

    表3 沙粒粒徑分布及含量分布

    由于沙粒中90%以上為二氧化硅,因此直接設(shè)置沙粒材料為二氧化硅物性參數(shù)如表4 所示,風(fēng)力機葉片涂層材料為聚氨酯物性參數(shù)如表5 所示[12]。

    表4 二氧化硅物性參數(shù)

    表5 聚氨酯涂層物性參數(shù)

    將涂層三維模型加載到EDEM,并按表5 設(shè)置涂層物性參數(shù);由于不同形狀沙粒對葉片涂層最大磨損率的變化規(guī)律基本一致,將沙塵顆粒簡化為球形顆粒,忽略非球形顆粒的影響[13]。建立球形沙粒模型,通過顆粒粒徑自定義函數(shù)按表3 設(shè)置沙粒粒徑分布參數(shù),設(shè)置接觸模型為Hertz-Mindlin,沙粒靜摩擦因數(shù)為0.5、動摩擦因數(shù)為0.15、恢復(fù)系數(shù)為0.5;沙粒靜摩擦因數(shù)為0.4、動摩擦因數(shù)為0.1、恢復(fù)系數(shù)為0.3。在流場入口建立四邊形顆粒工廠邊界,進而在邊界平面上建立動態(tài)顆粒工廠,并設(shè)置顆粒生成速率為1.2 kg/s,總質(zhì)量為5 kg。網(wǎng)格大小為最小顆粒半徑的20 倍,激活EDEM 耦合器,等待與Fluent 耦合仿真分析。

    3.4 Fluent 隨機風(fēng)蝕仿真分析

    假定接觸線周圍流場流體不可壓縮,而且湍流流動發(fā)展足夠充分,選擇標(biāo)準(zhǔn)大渦模擬模型,并設(shè)置模型參數(shù)為缺省值。采用基于壓力的求解器,流體域求解算法為SIMPLE 算法,流體域離散方法為二階迎風(fēng)離散格式[14]。

    3.5 控制方程

    3.5.1 連續(xù)性方程

    所有流動問題都必須遵循質(zhì)量守恒定律,即:單位時間內(nèi)流體微元質(zhì)量的增量等于該時間間隔內(nèi)流入該微元體的凈質(zhì)量,由此可推出連續(xù)性方程,其微分形式為

    式中: ρ為密度;t為時間;u、v、w分別為速度矢量u在x、y、z方向上的分量。

    對于不可壓縮流體,密度ρ 為常數(shù),則有

    3.5.2 動量守恒方程

    動量守恒定律是所有流動體系都必須遵循的基本定律,即:給定流體系統(tǒng)的動量對時間的變化率等于其所受外力的總和,根據(jù)該定律推導(dǎo)可得動量守恒方程,也稱為運動方程,或N-S 方程[15]:

    式中:p為靜壓; τxx、 τyx、 τzx、 τxy、 τyy、 τzy、 τxz、 τyz及τzz為 應(yīng)力張量;Fx、Fy及Fz為體力,若體力為重力,且Z軸豎直向上,則Fx=0,Fy=0,Fz=-ρg。

    3.6 流體域的離散

    以上論述表明CFD 數(shù)學(xué)模型由偏微分方程組組成,因此很難得到其解析解。目前主要通過對流體域進行離散化,進而在離散域節(jié)點上建立代數(shù)方程,然后用有限差分法、有限元法及有限體積法等數(shù)值方法對所得的代數(shù)方程進行求解。

    3.6.1 計算區(qū)域的離散

    對于有限體積法,其計算域的離散實質(zhì)是將計算域分解為多個相互獨立的子區(qū)域,并確定各子區(qū)域的節(jié)點位置和該節(jié)點所代表的控制體積,計算域離散后可得到節(jié)點、控制體積、界面及網(wǎng)格線4 種幾何要素。

    3.6.2 控制方程的離散

    對于三維瞬態(tài)流場,其離散方程為

    其中:

    4 仿真結(jié)果分析

    4.1 EDEM 和Fluent 耦合仿真結(jié)果分析

    流場仿真開始后,F(xiàn)luent 會通過耦合程序關(guān)聯(lián)啟動EDEM 進行顆粒系統(tǒng)的仿真計算,由于隨機風(fēng)沙沖蝕葉片涂層的模擬仿真為瞬態(tài)流場問題,所以要求流場在每個時間步都需迭代至收斂狀態(tài)(Fluent默認(rèn)殘差下降到10-3為收斂)才會進行下一個時間步的求解。

    整個仿真過程結(jié)束后可以得到入口速度、流固耦合面FSI 動態(tài)壓力的時間歷程數(shù)據(jù)以及沙粒沖擊葉片涂層的沖擊載荷,分別如圖9、圖10 以及圖11 所示。

    圖9 入口速度

    圖10 流固耦合面動壓

    圖11 隨機沖擊載荷

    由圖9 可知入口速度介于5~25 m/s 之間與風(fēng)場實測風(fēng)速數(shù)據(jù)較為吻合。

    由圖10 可知流固耦合面FSI 的動壓介于20~140 Pa 之間。

    由圖11 可知沙粒的沖擊載荷介于0.4~7.9 N 之間,三者均具有顯著的隨機性。

    4.2 涂層仿真結(jié)果分析

    通過上文對葉片表面的風(fēng)沙沖蝕涂層的隨機載荷分析,得到外載荷主要為沙粒的沖擊載荷和風(fēng)載荷。通過有限元軟件對隨機載荷作用到葉片涂層進行應(yīng)力響應(yīng)分析,得到葉片涂層在隨機載荷作用下應(yīng)力響應(yīng)的時間歷程數(shù)據(jù)以及葉片涂層應(yīng)力響應(yīng)最大時的應(yīng)力云圖,分別如圖12 和圖13 所示。

    圖12 隨機載荷時間歷程曲線

    圖13 0.86 s 時葉片涂層的應(yīng)力云圖

    由圖12 可知:風(fēng)沙環(huán)境下,葉片涂層的應(yīng)力響應(yīng)具有顯著的隨機性,最大應(yīng)力響應(yīng)介于0.4~2.79 MPa之間且主要集中在0.8~1.6 MPa 之間。0.86 s 時葉片涂層的應(yīng)力響應(yīng)在葉片尖端附近達到最大為2.79 MPa,如圖13 所示。

    5 結(jié)論

    基于EDEM 與Fluent 的耦合仿真分析得到了隨機風(fēng)沙沖蝕作用下葉片涂層的載荷數(shù)據(jù)自然風(fēng)速的隨機性以及沙粒直徑的不確定性使得風(fēng)沙沖蝕葉片涂層的過程具有顯著的隨機性。

    沖蝕過程中:風(fēng)速V=13.9×[-ln(1-U)]1/2.54,葉片涂層所受的外載荷介于0.4~7.9 N 之間,葉片涂層的應(yīng)力響應(yīng)介于0.4~2.79 MPa之間且主要集中在0.8~1.6 MPa 之間。

    葉片涂層的最大應(yīng)力響應(yīng)為2.79 MPa,位于葉片尖端前緣及其側(cè)面附近,由于最大應(yīng)力小于涂層材料的屈服強度,所以葉片涂層受載后未發(fā)生明顯的塑性變形,因此風(fēng)沙環(huán)境下葉片涂層的失效形式表現(xiàn)為隨機脈動應(yīng)力下的高周疲勞破壞,而且葉片涂層的疲勞損傷主要集中在葉片尖端前緣及其側(cè)面附近。

    根據(jù)計算結(jié)果,確定葉片涂層疲勞壽命,建立涂層精細化設(shè)計方法, 并提出預(yù)防性維修方案,提高葉片涂層防沙粒侵蝕性能。

    猜你喜歡
    沙粒沖蝕風(fēng)沙
    140MPa井口壓裂四通管道沖蝕分析
    沙粒和水珠
    風(fēng)沙地區(qū)公路路基防護措施初探
    想看山的小沙粒
    想看山的小沙粒
    東明縣風(fēng)沙化土地監(jiān)測與治理
    輸氣管道砂沖蝕的模擬實驗
    高考試題中“風(fēng)沙、河沙”知識點的解讀
    環(huán)氧樹脂及其復(fù)合材料的固體顆粒沖蝕磨損
    三種不銹鋼材料抗固體顆粒沖蝕性能研究
    亚洲成av人片免费观看| 欧美一区二区国产精品久久精品| 亚洲五月天丁香| 99久久精品国产亚洲精品| 欧美最黄视频在线播放免费| 欧美三级亚洲精品| 国产高清激情床上av| 色综合站精品国产| 非洲黑人性xxxx精品又粗又长| 一卡2卡三卡四卡精品乱码亚洲| 啦啦啦韩国在线观看视频| 国产大屁股一区二区在线视频| 国产精品亚洲美女久久久| 国产老妇女一区| 美女免费视频网站| av在线天堂中文字幕| 久9热在线精品视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 五月伊人婷婷丁香| 免费大片18禁| 一级av片app| 亚洲中文日韩欧美视频| 欧美又色又爽又黄视频| 我要看日韩黄色一级片| 老司机午夜十八禁免费视频| 国产精品野战在线观看| 女同久久另类99精品国产91| 国产精品一及| 男插女下体视频免费在线播放| 欧美极品一区二区三区四区| 一个人看视频在线观看www免费| 中文字幕av成人在线电影| 久久人人爽人人爽人人片va | 9191精品国产免费久久| 97碰自拍视频| 国产一区二区三区在线臀色熟女| 少妇人妻精品综合一区二区 | 国产高清激情床上av| 69人妻影院| 欧美日韩福利视频一区二区| av专区在线播放| 中文亚洲av片在线观看爽| 久久精品国产自在天天线| 国内毛片毛片毛片毛片毛片| 首页视频小说图片口味搜索| 波多野结衣巨乳人妻| 欧美三级亚洲精品| 免费高清视频大片| 亚洲国产色片| 成人精品一区二区免费| 日本免费一区二区三区高清不卡| 亚洲av美国av| 久久午夜亚洲精品久久| 性欧美人与动物交配| 高清在线国产一区| 国产伦在线观看视频一区| 国内精品一区二区在线观看| 婷婷色综合大香蕉| 欧美黑人巨大hd| 蜜桃久久精品国产亚洲av| 熟妇人妻久久中文字幕3abv| www日本黄色视频网| 舔av片在线| 国产成人av教育| 高潮久久久久久久久久久不卡| 全区人妻精品视频| 久久中文看片网| 天天躁日日操中文字幕| 久久久久精品国产欧美久久久| 人妻制服诱惑在线中文字幕| 婷婷六月久久综合丁香| 69av精品久久久久久| 亚洲av不卡在线观看| 欧美绝顶高潮抽搐喷水| 自拍偷自拍亚洲精品老妇| 欧美一区二区国产精品久久精品| 99热这里只有是精品在线观看 | 女人十人毛片免费观看3o分钟| 又爽又黄无遮挡网站| 91狼人影院| 亚洲精品亚洲一区二区| 偷拍熟女少妇极品色| 亚洲人成伊人成综合网2020| 国产伦在线观看视频一区| 午夜福利在线在线| 免费看美女性在线毛片视频| 99久国产av精品| 91狼人影院| 精品久久久久久,| 欧美一区二区亚洲| 人妻久久中文字幕网| 亚洲va日本ⅴa欧美va伊人久久| 亚洲欧美日韩无卡精品| 精品久久久久久久末码| 国产 一区 欧美 日韩| 变态另类成人亚洲欧美熟女| 99热只有精品国产| 国产淫片久久久久久久久 | 精品国产三级普通话版| 一本久久中文字幕| av黄色大香蕉| 欧美日韩黄片免| 香蕉av资源在线| 亚洲精品影视一区二区三区av| 色av中文字幕| 欧美区成人在线视频| 色综合亚洲欧美另类图片| 18禁裸乳无遮挡免费网站照片| 乱码一卡2卡4卡精品| 久久香蕉精品热| 88av欧美| 美女高潮的动态| 亚洲一区二区三区不卡视频| 少妇人妻一区二区三区视频| 18美女黄网站色大片免费观看| 黄色丝袜av网址大全| 黄色日韩在线| 天堂√8在线中文| 好男人在线观看高清免费视频| 中文字幕久久专区| 12—13女人毛片做爰片一| 亚洲18禁久久av| 亚洲av免费高清在线观看| 国产色婷婷99| 12—13女人毛片做爰片一| 丰满乱子伦码专区| 亚洲精品456在线播放app | 久久精品国产自在天天线| 99国产精品一区二区蜜桃av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 男女视频在线观看网站免费| 天堂网av新在线| 国产高清视频在线播放一区| 国产伦人伦偷精品视频| 老司机福利观看| 国产乱人伦免费视频| 俺也久久电影网| 怎么达到女性高潮| 久久国产乱子伦精品免费另类| 九九热线精品视视频播放| 午夜激情福利司机影院| 国产精品乱码一区二三区的特点| 黄片小视频在线播放| 禁无遮挡网站| 麻豆一二三区av精品| 此物有八面人人有两片| 美女xxoo啪啪120秒动态图 | 久久久久久国产a免费观看| 桃色一区二区三区在线观看| 亚洲av二区三区四区| 国产精品一及| 美女 人体艺术 gogo| 久久久久免费精品人妻一区二区| 国产野战对白在线观看| 啪啪无遮挡十八禁网站| 日日夜夜操网爽| 又粗又爽又猛毛片免费看| 淫秽高清视频在线观看| 老司机福利观看| 一级黄色大片毛片| 成年女人永久免费观看视频| 蜜桃亚洲精品一区二区三区| 久久香蕉精品热| 亚洲国产欧洲综合997久久,| 美女高潮的动态| 激情在线观看视频在线高清| 天堂√8在线中文| 99热6这里只有精品| 好看av亚洲va欧美ⅴa在| 亚洲精品成人久久久久久| 精品久久久久久久末码| 看片在线看免费视频| 2021天堂中文幕一二区在线观| 日日干狠狠操夜夜爽| 国产aⅴ精品一区二区三区波| 免费在线观看影片大全网站| 欧美性猛交黑人性爽| 日本黄色片子视频| 在线天堂最新版资源| АⅤ资源中文在线天堂| 国产成人欧美在线观看| 欧美日韩黄片免| 国产精品爽爽va在线观看网站| 麻豆av噜噜一区二区三区| 伦理电影大哥的女人| 香蕉av资源在线| 国产真实乱freesex| 国产精品亚洲美女久久久| 99久久精品热视频| 国产亚洲欧美98| 国产乱人伦免费视频| 国产欧美日韩精品亚洲av| 美女高潮的动态| 亚洲精品一卡2卡三卡4卡5卡| 免费高清视频大片| 成熟少妇高潮喷水视频| 国内少妇人妻偷人精品xxx网站| 久久精品综合一区二区三区| 免费黄网站久久成人精品 | 99热这里只有是精品在线观看 | 国产人妻一区二区三区在| 国产亚洲精品综合一区在线观看| 午夜影院日韩av| 淫秽高清视频在线观看| 久久久久国内视频| 国产淫片久久久久久久久 | 精品人妻偷拍中文字幕| 亚洲熟妇中文字幕五十中出| 丝袜美腿在线中文| 国产成人欧美在线观看| 亚洲午夜理论影院| 亚洲第一区二区三区不卡| 欧美+亚洲+日韩+国产| 丰满的人妻完整版| av视频在线观看入口| x7x7x7水蜜桃| 国产在线男女| 两个人的视频大全免费| 宅男免费午夜| 日本 欧美在线| 国产亚洲欧美98| 人人妻人人看人人澡| 一区二区三区高清视频在线| 一a级毛片在线观看| 一级毛片久久久久久久久女| 熟女人妻精品中文字幕| 99精品在免费线老司机午夜| 久久亚洲精品不卡| 热99re8久久精品国产| 久久精品国产亚洲av涩爱 | 成人欧美大片| 久久精品国产亚洲av香蕉五月| 午夜激情福利司机影院| 亚洲中文字幕一区二区三区有码在线看| 免费看美女性在线毛片视频| 国产欧美日韩精品一区二区| 亚洲av免费在线观看| 99热这里只有是精品在线观看 | 51国产日韩欧美| 亚洲欧美日韩东京热| 欧美黑人欧美精品刺激| 国产单亲对白刺激| 不卡一级毛片| 在线观看美女被高潮喷水网站 | 少妇人妻一区二区三区视频| 99精品久久久久人妻精品| 国产成+人综合+亚洲专区| 精品一区二区免费观看| 性欧美人与动物交配| 夜夜看夜夜爽夜夜摸| 51午夜福利影视在线观看| 怎么达到女性高潮| 一进一出好大好爽视频| 亚洲第一欧美日韩一区二区三区| 亚洲精品在线美女| 久久久久久久精品吃奶| 成人精品一区二区免费| 88av欧美| 亚洲国产精品久久男人天堂| 久久人人精品亚洲av| 久久久久国内视频| 人妻制服诱惑在线中文字幕| 国产一区二区三区在线臀色熟女| 亚洲在线观看片| 久久国产乱子伦精品免费另类| 久久性视频一级片| 757午夜福利合集在线观看| 亚洲自偷自拍三级| 国产精品永久免费网站| 日韩欧美精品v在线| 国产精华一区二区三区| 久久久久国产精品人妻aⅴ院| 十八禁网站免费在线| 人妻久久中文字幕网| 欧美一区二区精品小视频在线| 啪啪无遮挡十八禁网站| 亚洲,欧美,日韩| 亚洲18禁久久av| 91麻豆av在线| 亚洲真实伦在线观看| 亚洲久久久久久中文字幕| 国产高清视频在线播放一区| 99国产精品一区二区蜜桃av| netflix在线观看网站| 深夜精品福利| 国产成人啪精品午夜网站| 亚洲av成人不卡在线观看播放网| 男插女下体视频免费在线播放| 内射极品少妇av片p| 深夜精品福利| 又爽又黄a免费视频| 午夜影院日韩av| 日韩亚洲欧美综合| 国产视频一区二区在线看| 精品免费久久久久久久清纯| 在线播放无遮挡| 国产精品,欧美在线| xxxwww97欧美| 在线观看av片永久免费下载| 一本一本综合久久| 麻豆成人av在线观看| 极品教师在线免费播放| 欧美高清成人免费视频www| 男女床上黄色一级片免费看| 尤物成人国产欧美一区二区三区| 午夜福利欧美成人| 麻豆av噜噜一区二区三区| 老熟妇仑乱视频hdxx| 啪啪无遮挡十八禁网站| 日本熟妇午夜| 超碰av人人做人人爽久久| 国产aⅴ精品一区二区三区波| 久久国产乱子免费精品| 国产中年淑女户外野战色| 99视频精品全部免费 在线| 观看美女的网站| 免费一级毛片在线播放高清视频| 在线免费观看的www视频| 色av中文字幕| 婷婷色综合大香蕉| 中出人妻视频一区二区| 午夜福利在线在线| 18禁黄网站禁片午夜丰满| 九九在线视频观看精品| 久久久久国内视频| 日本成人三级电影网站| 天天躁日日操中文字幕| 网址你懂的国产日韩在线| 深爱激情五月婷婷| 国产精品不卡视频一区二区 | 国产精品99久久久久久久久| 久久久精品欧美日韩精品| 又黄又爽又刺激的免费视频.| 五月伊人婷婷丁香| 亚洲一区高清亚洲精品| 久久欧美精品欧美久久欧美| 国产人妻一区二区三区在| 国产精品久久久久久久电影| 欧美在线一区亚洲| 少妇熟女aⅴ在线视频| 国产精品爽爽va在线观看网站| 直男gayav资源| 久久人人精品亚洲av| 麻豆国产av国片精品| h日本视频在线播放| 亚洲熟妇熟女久久| 成年女人毛片免费观看观看9| 婷婷亚洲欧美| 欧美日本亚洲视频在线播放| 香蕉av资源在线| 欧美xxxx黑人xx丫x性爽| 欧美国产日韩亚洲一区| 少妇的逼水好多| 久久国产乱子免费精品| 国模一区二区三区四区视频| 高潮久久久久久久久久久不卡| 日韩人妻高清精品专区| 国产黄色小视频在线观看| 高清在线国产一区| 熟女电影av网| 黄片小视频在线播放| 欧美日韩福利视频一区二区| 国产精品不卡视频一区二区 | 久久精品国产亚洲av天美| 亚洲精品在线观看二区| 51午夜福利影视在线观看| 别揉我奶头 嗯啊视频| 欧美不卡视频在线免费观看| 亚洲狠狠婷婷综合久久图片| 午夜福利成人在线免费观看| 国产精品综合久久久久久久免费| 亚洲七黄色美女视频| 中文在线观看免费www的网站| 此物有八面人人有两片| 日韩欧美精品免费久久 | 日韩欧美免费精品| 最新在线观看一区二区三区| 国产午夜福利久久久久久| 嫩草影院新地址| 日韩大尺度精品在线看网址| 国产不卡一卡二| 男人和女人高潮做爰伦理| 亚洲国产精品合色在线| 国产亚洲欧美在线一区二区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美一区二区亚洲| 国产色婷婷99| 最新中文字幕久久久久| 国产精华一区二区三区| 搡老熟女国产l中国老女人| 久久精品影院6| 亚洲精品乱码久久久v下载方式| 又粗又爽又猛毛片免费看| 亚洲av第一区精品v没综合| 国产在线精品亚洲第一网站| 色噜噜av男人的天堂激情| 中国美女看黄片| 午夜视频国产福利| 欧美+亚洲+日韩+国产| 两个人视频免费观看高清| 国内精品久久久久精免费| 久久欧美精品欧美久久欧美| 色播亚洲综合网| 久久久久久久亚洲中文字幕 | 桃色一区二区三区在线观看| 又紧又爽又黄一区二区| 午夜福利免费观看在线| 亚洲av电影不卡..在线观看| 99久久久亚洲精品蜜臀av| 日韩欧美免费精品| 亚洲在线观看片| 97人妻精品一区二区三区麻豆| 午夜免费激情av| 91在线观看av| 国产91精品成人一区二区三区| 成年免费大片在线观看| 激情在线观看视频在线高清| 久久久久久国产a免费观看| 国产久久久一区二区三区| 久久欧美精品欧美久久欧美| 国产亚洲欧美98| 亚洲美女黄片视频| 亚洲人成网站在线播| 国产伦精品一区二区三区视频9| 又爽又黄a免费视频| 黄片小视频在线播放| 国产 一区 欧美 日韩| 99热只有精品国产| 免费观看精品视频网站| 欧美一区二区国产精品久久精品| 91在线观看av| 男女之事视频高清在线观看| 黄色日韩在线| 男人的好看免费观看在线视频| 日韩中文字幕欧美一区二区| 国产精品嫩草影院av在线观看 | 2021天堂中文幕一二区在线观| 日本黄色片子视频| 美女免费视频网站| 免费黄网站久久成人精品 | a级一级毛片免费在线观看| 国产国拍精品亚洲av在线观看| a级一级毛片免费在线观看| 99精品久久久久人妻精品| 1000部很黄的大片| 美女cb高潮喷水在线观看| 亚洲avbb在线观看| 亚洲 国产 在线| 老司机午夜十八禁免费视频| 亚洲国产精品999在线| 国产中年淑女户外野战色| 午夜影院日韩av| 国产精品久久久久久久久免 | 欧美黄色淫秽网站| 亚洲片人在线观看| 中文在线观看免费www的网站| 99热只有精品国产| 美女被艹到高潮喷水动态| 嫩草影院精品99| 91九色精品人成在线观看| 一个人免费在线观看的高清视频| 91九色精品人成在线观看| 一夜夜www| 免费无遮挡裸体视频| 国模一区二区三区四区视频| 亚洲成人久久性| 久久久久性生活片| av国产免费在线观看| 成人亚洲精品av一区二区| 久久精品夜夜夜夜夜久久蜜豆| 在线播放国产精品三级| 十八禁网站免费在线| 久久久久久久久中文| 我要搜黄色片| 少妇高潮的动态图| 一本久久中文字幕| 日韩免费av在线播放| 国内精品久久久久久久电影| 精品人妻视频免费看| 精品久久久久久成人av| 高清日韩中文字幕在线| 欧美潮喷喷水| 一个人免费在线观看电影| 看十八女毛片水多多多| 一本久久中文字幕| 麻豆久久精品国产亚洲av| 国产人妻一区二区三区在| 亚洲乱码一区二区免费版| 免费黄网站久久成人精品 | 免费在线观看日本一区| 综合色av麻豆| 757午夜福利合集在线观看| 18美女黄网站色大片免费观看| 国产极品精品免费视频能看的| 99精品久久久久人妻精品| 精品人妻1区二区| 极品教师在线视频| 日本 av在线| 久久草成人影院| 免费看a级黄色片| 成熟少妇高潮喷水视频| 久久这里只有精品中国| 99热这里只有精品一区| 日本 av在线| 中文字幕久久专区| 观看美女的网站| 免费av不卡在线播放| 桃红色精品国产亚洲av| 成人国产一区最新在线观看| 婷婷丁香在线五月| 五月玫瑰六月丁香| 国产成+人综合+亚洲专区| 中国美女看黄片| 日本免费一区二区三区高清不卡| 男人狂女人下面高潮的视频| 国产精品av视频在线免费观看| 亚洲在线观看片| 国产一区二区三区在线臀色熟女| 亚洲最大成人中文| 一边摸一边抽搐一进一小说| 成人av一区二区三区在线看| 婷婷丁香在线五月| 久久久久久久久久黄片| 午夜日韩欧美国产| а√天堂www在线а√下载| 国产一级毛片七仙女欲春2| 午夜激情欧美在线| 日韩欧美一区二区三区在线观看| 精品久久久久久久人妻蜜臀av| 日日干狠狠操夜夜爽| 我的老师免费观看完整版| 欧美丝袜亚洲另类 | 日韩av在线大香蕉| av在线蜜桃| 亚洲欧美精品综合久久99| 亚洲av成人av| 婷婷亚洲欧美| 51国产日韩欧美| 午夜a级毛片| 午夜视频国产福利| 欧美xxxx黑人xx丫x性爽| 日韩欧美国产一区二区入口| 级片在线观看| 偷拍熟女少妇极品色| 夜夜看夜夜爽夜夜摸| 成人鲁丝片一二三区免费| 亚州av有码| 日韩国内少妇激情av| 1024手机看黄色片| 婷婷色综合大香蕉| 美女被艹到高潮喷水动态| 国产精品久久久久久人妻精品电影| 免费在线观看影片大全网站| 精品久久久久久成人av| a级毛片免费高清观看在线播放| 免费av不卡在线播放| 久久久久性生活片| 99热只有精品国产| 伊人久久精品亚洲午夜| 国产精品嫩草影院av在线观看 | 天堂动漫精品| 一个人看视频在线观看www免费| 午夜福利视频1000在线观看| 亚洲自拍偷在线| 日本一本二区三区精品| 岛国在线免费视频观看| 国产高潮美女av| 亚洲不卡免费看| 乱人视频在线观看| 精品一区二区三区视频在线观看免费| 亚洲第一欧美日韩一区二区三区| 国产国拍精品亚洲av在线观看| 白带黄色成豆腐渣| 亚洲欧美激情综合另类| 深夜精品福利| 国产探花极品一区二区| 欧美日韩福利视频一区二区| 热99在线观看视频| 在线国产一区二区在线| 禁无遮挡网站| 欧美bdsm另类| 赤兔流量卡办理| 国模一区二区三区四区视频| 国产午夜精品久久久久久一区二区三区 | 91在线观看av| 国产精品,欧美在线| 亚洲无线观看免费| 露出奶头的视频| 欧美激情久久久久久爽电影| 男人和女人高潮做爰伦理| 国产精品综合久久久久久久免费| 特大巨黑吊av在线直播| 少妇的逼好多水| 亚洲av二区三区四区| 桃色一区二区三区在线观看| 亚洲av成人av| 亚洲专区国产一区二区| 麻豆av噜噜一区二区三区| 免费高清视频大片| .国产精品久久| 99热这里只有是精品在线观看 | 美女 人体艺术 gogo| 亚洲,欧美精品.| 人妻制服诱惑在线中文字幕| 少妇人妻精品综合一区二区 | 色综合欧美亚洲国产小说| 国产欧美日韩一区二区三| 欧美乱色亚洲激情| 欧美色视频一区免费| 亚洲精品在线观看二区| 高潮久久久久久久久久久不卡| 日本 欧美在线| 亚洲成人久久性| 亚洲天堂国产精品一区在线| 宅男免费午夜|