• <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ù)合材料的固體顆粒沖蝕磨損
    三種不銹鋼材料抗固體顆粒沖蝕性能研究
    国产欧美日韩综合在线一区二区| 国产熟女欧美一区二区| 一本—道久久a久久精品蜜桃钙片| 丝袜脚勾引网站| 精品国产乱码久久久久久小说| 久久av网站| 日日啪夜夜爽| 欧美 亚洲 国产 日韩一| 黑人巨大精品欧美一区二区蜜桃 | 亚洲性久久影院| 久久精品国产自在天天线| 精品久久久精品久久久| 久久毛片免费看一区二区三区| 啦啦啦啦在线视频资源| 久久久久久久久久成人| av免费观看日本| 中文字幕av电影在线播放| 国产一级毛片在线| 交换朋友夫妻互换小说| 国产免费一区二区三区四区乱码| 国产免费又黄又爽又色| 最新的欧美精品一区二区| 男女国产视频网站| 91成人精品电影| 边亲边吃奶的免费视频| xxxhd国产人妻xxx| 一级毛片aaaaaa免费看小| 搡女人真爽免费视频火全软件| 最近中文字幕高清免费大全6| 日韩亚洲欧美综合| 日韩欧美一区视频在线观看| 日韩av在线免费看完整版不卡| 国模一区二区三区四区视频| 在线观看免费视频网站a站| 亚洲国产av新网站| 亚洲婷婷狠狠爱综合网| 18+在线观看网站| 亚州av有码| 欧美3d第一页| 亚洲欧美中文字幕日韩二区| 国产成人freesex在线| 亚洲成人手机| 亚洲欧美精品自产自拍| 人妻少妇偷人精品九色| 一区二区三区精品91| 亚洲av不卡在线观看| 日韩大片免费观看网站| av专区在线播放| 综合色丁香网| 欧美最新免费一区二区三区| 成年人午夜在线观看视频| freevideosex欧美| 欧美日韩亚洲高清精品| 一本色道久久久久久精品综合| 天天影视国产精品| 黄片播放在线免费| 亚洲成色77777| 免费日韩欧美在线观看| 最新的欧美精品一区二区| 高清毛片免费看| 国产国拍精品亚洲av在线观看| 另类亚洲欧美激情| 国精品久久久久久国模美| 欧美日韩一区二区视频在线观看视频在线| 在线 av 中文字幕| av播播在线观看一区| 久久久久久久久久久丰满| 多毛熟女@视频| 五月开心婷婷网| 亚洲国产精品成人久久小说| 国产午夜精品久久久久久一区二区三区| 精品亚洲乱码少妇综合久久| 大片电影免费在线观看免费| 亚洲无线观看免费| 能在线免费看毛片的网站| a级毛片在线看网站| av线在线观看网站| 99久久人妻综合| 男人添女人高潮全过程视频| 午夜免费鲁丝| 国产 一区精品| 永久免费av网站大全| 最近中文字幕2019免费版| 人妻系列 视频| 一本色道久久久久久精品综合| 午夜福利,免费看| 亚洲精品国产av成人精品| 亚洲色图 男人天堂 中文字幕 | 各种免费的搞黄视频| 亚洲,欧美,日韩| 国产精品熟女久久久久浪| 免费观看性生交大片5| 亚洲av免费高清在线观看| 国产av国产精品国产| 亚洲欧美成人综合另类久久久| 两个人的视频大全免费| 成人亚洲精品一区在线观看| 十八禁网站网址无遮挡| 狂野欧美白嫩少妇大欣赏| 亚州av有码| 亚洲成人av在线免费| 18禁动态无遮挡网站| 国产精品欧美亚洲77777| 亚洲经典国产精华液单| 麻豆乱淫一区二区| av又黄又爽大尺度在线免费看| 久久久久国产网址| 亚洲人与动物交配视频| 亚洲av.av天堂| 亚洲国产成人一精品久久久| 久久久久久久久大av| 欧美老熟妇乱子伦牲交| videossex国产| 亚洲熟女精品中文字幕| 日韩电影二区| 18禁在线无遮挡免费观看视频| 亚洲精品自拍成人| 又黄又爽又刺激的免费视频.| 天美传媒精品一区二区| 国产av码专区亚洲av| 久久99热6这里只有精品| 国产男人的电影天堂91| videosex国产| 国产精品蜜桃在线观看| 中国国产av一级| 我的老师免费观看完整版| 99久久人妻综合| 久久国产亚洲av麻豆专区| 欧美精品一区二区免费开放| 久久精品熟女亚洲av麻豆精品| 日本av手机在线免费观看| 亚洲精品视频女| 欧美97在线视频| 亚洲综合精品二区| 一区在线观看完整版| av专区在线播放| 中文乱码字字幕精品一区二区三区| 亚洲欧美成人综合另类久久久| 人妻一区二区av| av女优亚洲男人天堂| 97超视频在线观看视频| 国产av精品麻豆| 黑人巨大精品欧美一区二区蜜桃 | 99热国产这里只有精品6| 91精品三级在线观看| 人妻系列 视频| 亚洲欧洲精品一区二区精品久久久 | 日日摸夜夜添夜夜添av毛片| 美女脱内裤让男人舔精品视频| 51国产日韩欧美| 五月玫瑰六月丁香| 亚洲成人av在线免费| av线在线观看网站| 少妇猛男粗大的猛烈进出视频| 国产一区二区三区av在线| 欧美+日韩+精品| 久久99热6这里只有精品| 亚洲av电影在线观看一区二区三区| 91午夜精品亚洲一区二区三区| 男人添女人高潮全过程视频| 亚洲精品日本国产第一区| 国产精品人妻久久久影院| 精品国产露脸久久av麻豆| 在线观看人妻少妇| av国产久精品久网站免费入址| 亚洲国产av影院在线观看| 高清视频免费观看一区二区| 在线免费观看不下载黄p国产| 夫妻午夜视频| 亚洲第一区二区三区不卡| 亚洲国产毛片av蜜桃av| 老司机影院成人| 一边亲一边摸免费视频| 永久免费av网站大全| 高清午夜精品一区二区三区| 人人妻人人澡人人看| 国产欧美日韩综合在线一区二区| 日本午夜av视频| 成人黄色视频免费在线看| 国产精品一国产av| 国产精品国产三级国产av玫瑰| 成年美女黄网站色视频大全免费 | 日韩中字成人| 久久久久久久久久人人人人人人| 少妇精品久久久久久久| 狂野欧美激情性bbbbbb| 观看美女的网站| 26uuu在线亚洲综合色| 欧美最新免费一区二区三区| 亚洲欧美色中文字幕在线| 久久久久网色| 美女内射精品一级片tv| 久久精品夜色国产| 妹子高潮喷水视频| 日本黄色片子视频| 纯流量卡能插随身wifi吗| 自线自在国产av| 中文字幕人妻熟人妻熟丝袜美| www.av在线官网国产| 久久久a久久爽久久v久久| 在线免费观看不下载黄p国产| av在线播放精品| 久久久久视频综合| 国产精品偷伦视频观看了| 国产成人精品福利久久| 国产高清三级在线| 亚洲人与动物交配视频| 高清不卡的av网站| 十分钟在线观看高清视频www| 日韩人妻高清精品专区| 91国产中文字幕| 天美传媒精品一区二区| 国产日韩欧美视频二区| 国产伦理片在线播放av一区| 欧美日韩视频精品一区| 亚洲婷婷狠狠爱综合网| 久久精品国产鲁丝片午夜精品| 亚洲成人一二三区av| 一区二区三区乱码不卡18| 国国产精品蜜臀av免费| 欧美人与性动交α欧美精品济南到 | 51国产日韩欧美| 91精品国产国语对白视频| 国产欧美另类精品又又久久亚洲欧美| 国产极品天堂在线| 国产高清不卡午夜福利| 一区二区日韩欧美中文字幕 | 3wmmmm亚洲av在线观看| av.在线天堂| 91精品三级在线观看| 一级爰片在线观看| 超碰97精品在线观看| 国产精品久久久久成人av| 人人澡人人妻人| 亚洲,欧美,日韩| 少妇精品久久久久久久| 婷婷色av中文字幕| 能在线免费看毛片的网站| 成人影院久久| 97超视频在线观看视频| 如日韩欧美国产精品一区二区三区 | 精品亚洲成国产av| 在线观看www视频免费| 一本久久精品| 视频中文字幕在线观看| 国产成人91sexporn| 日本与韩国留学比较| 老司机影院毛片| 亚洲第一av免费看| 亚洲av综合色区一区| 在线 av 中文字幕| 在线观看免费高清a一片| 成人影院久久| 黑人高潮一二区| 只有这里有精品99| 制服丝袜香蕉在线| 日韩熟女老妇一区二区性免费视频| av视频免费观看在线观看| 乱人伦中国视频| 啦啦啦啦在线视频资源| 国产精品嫩草影院av在线观看| 日日摸夜夜添夜夜爱| 黑人猛操日本美女一级片| 一级毛片 在线播放| 欧美精品一区二区免费开放| 两个人的视频大全免费| 成人亚洲精品一区在线观看| 国产成人一区二区在线| 精品久久久噜噜| 亚洲人成网站在线播| 亚洲五月色婷婷综合| 一本一本综合久久| 亚洲高清免费不卡视频| 中国国产av一级| 色哟哟·www| av免费在线看不卡| 99热国产这里只有精品6| 大香蕉久久成人网| 国产精品 国内视频| 欧美xxⅹ黑人| 黑丝袜美女国产一区| 久久鲁丝午夜福利片| 熟女av电影| 另类亚洲欧美激情| 日韩视频在线欧美| 午夜福利视频在线观看免费| 亚洲国产成人一精品久久久| 青春草国产在线视频| 欧美精品亚洲一区二区| 一本久久精品| 精品国产国语对白av| 最近中文字幕高清免费大全6| av天堂久久9| 亚洲欧洲日产国产| 日韩三级伦理在线观看| 精品国产一区二区久久| 国产一区二区三区av在线| freevideosex欧美| 久久久亚洲精品成人影院| 人妻夜夜爽99麻豆av| 女人精品久久久久毛片| 亚洲欧美中文字幕日韩二区| 亚洲国产精品国产精品| 在线播放无遮挡| 成年美女黄网站色视频大全免费 | 一边摸一边做爽爽视频免费| 18+在线观看网站| 免费看不卡的av| a级片在线免费高清观看视频| 另类亚洲欧美激情| 欧美老熟妇乱子伦牲交| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品自拍成人| 人人妻人人添人人爽欧美一区卜| 亚洲精品一区蜜桃| 永久免费av网站大全| 91久久精品电影网| 亚洲激情五月婷婷啪啪| 天天躁夜夜躁狠狠久久av| 男女边摸边吃奶| 国产极品天堂在线| 国产精品无大码| 色婷婷av一区二区三区视频| 国产老妇伦熟女老妇高清| 精品亚洲成a人片在线观看| 欧美人与善性xxx| freevideosex欧美| 日韩av免费高清视频| 最近手机中文字幕大全| 18禁动态无遮挡网站| www.av在线官网国产| 国产免费一级a男人的天堂| 中文天堂在线官网| 汤姆久久久久久久影院中文字幕| 亚洲精品久久久久久婷婷小说| 国产黄频视频在线观看| 男人操女人黄网站| 一区二区av电影网| 亚洲一级一片aⅴ在线观看| 亚洲精品乱久久久久久| 2018国产大陆天天弄谢| 伦理电影免费视频| 午夜免费鲁丝| 18禁在线无遮挡免费观看视频| 美女cb高潮喷水在线观看| av女优亚洲男人天堂| 一级黄片播放器| 亚洲国产欧美日韩在线播放| 国产在线一区二区三区精| 国产精品麻豆人妻色哟哟久久| 亚洲丝袜综合中文字幕| 久久av网站| 亚洲成色77777| 天堂中文最新版在线下载| 国产永久视频网站| 亚洲美女黄色视频免费看| 3wmmmm亚洲av在线观看| 最近最新中文字幕免费大全7| 免费久久久久久久精品成人欧美视频 | 国产精品嫩草影院av在线观看| 欧美精品亚洲一区二区| 欧美3d第一页| av在线app专区| 人妻 亚洲 视频| 国产探花极品一区二区| 免费高清在线观看日韩| 边亲边吃奶的免费视频| 国产午夜精品久久久久久一区二区三区| 又大又黄又爽视频免费| 一本—道久久a久久精品蜜桃钙片| 色哟哟·www| 91在线精品国自产拍蜜月| 一级毛片我不卡| 亚洲精品第二区| 亚洲图色成人| 在线观看免费视频网站a站| 欧美+日韩+精品| 高清午夜精品一区二区三区| 亚洲av中文av极速乱| 欧美日韩精品成人综合77777| 免费不卡的大黄色大毛片视频在线观看| 国产男女内射视频| 久久97久久精品| 成年人午夜在线观看视频| 国产精品久久久久久精品古装| videos熟女内射| a级片在线免费高清观看视频| 欧美最新免费一区二区三区| 天天躁夜夜躁狠狠久久av| av福利片在线| av又黄又爽大尺度在线免费看| 欧美性感艳星| 能在线免费看毛片的网站| 99热6这里只有精品| 久久久久久久久久久丰满| 国产精品国产av在线观看| 久久久久久久久久人人人人人人| 亚洲激情五月婷婷啪啪| 菩萨蛮人人尽说江南好唐韦庄| 国产精品免费大片| 毛片一级片免费看久久久久| 下体分泌物呈黄色| 九九爱精品视频在线观看| 80岁老熟妇乱子伦牲交| 波野结衣二区三区在线| 男女边摸边吃奶| 久久久久久久久大av| 人人澡人人妻人| 在线 av 中文字幕| 国产爽快片一区二区三区| 亚洲人成网站在线观看播放| 久久女婷五月综合色啪小说| 久久免费观看电影| 女人精品久久久久毛片| 视频中文字幕在线观看| 丝袜脚勾引网站| 精品亚洲乱码少妇综合久久| 久久av网站| 人成视频在线观看免费观看| 99热国产这里只有精品6| 狂野欧美激情性xxxx在线观看| 夜夜骑夜夜射夜夜干| 久热久热在线精品观看| 国产av一区二区精品久久| 国产熟女欧美一区二区| 哪个播放器可以免费观看大片| 大片电影免费在线观看免费| 亚洲婷婷狠狠爱综合网| 妹子高潮喷水视频| 久久久精品免费免费高清| 大陆偷拍与自拍| 亚洲天堂av无毛| 成人影院久久| av女优亚洲男人天堂| 精品亚洲成国产av| 性高湖久久久久久久久免费观看| 国产极品天堂在线| 国产男女内射视频| 亚洲欧美清纯卡通| 丝瓜视频免费看黄片| 99国产精品免费福利视频| 午夜日本视频在线| 九色成人免费人妻av| 99热这里只有精品一区| 日本av免费视频播放| 2021少妇久久久久久久久久久| 高清毛片免费看| 在线精品无人区一区二区三| 永久免费av网站大全| 亚洲av国产av综合av卡| 国产伦精品一区二区三区视频9| 十分钟在线观看高清视频www| 欧美少妇被猛烈插入视频| 精品一区二区三区视频在线| 午夜福利视频在线观看免费| 久久精品人人爽人人爽视色| a级毛片免费高清观看在线播放| 欧美一级a爱片免费观看看| 在线观看免费日韩欧美大片 | 国产精品99久久久久久久久| 色婷婷av一区二区三区视频| 熟女人妻精品中文字幕| 看免费成人av毛片| 精品一区在线观看国产| 日韩 亚洲 欧美在线| 在线观看免费视频网站a站| 纵有疾风起免费观看全集完整版| 女人精品久久久久毛片| 最后的刺客免费高清国语| 菩萨蛮人人尽说江南好唐韦庄| 免费日韩欧美在线观看| 伦精品一区二区三区| 最近最新中文字幕免费大全7| 99热全是精品| 精品国产露脸久久av麻豆| 亚洲精品av麻豆狂野| 久久久久久久国产电影| 成人亚洲欧美一区二区av| 十分钟在线观看高清视频www| 婷婷色综合www| 欧美一级a爱片免费观看看| 人妻夜夜爽99麻豆av| 精品久久国产蜜桃| 97精品久久久久久久久久精品| 另类亚洲欧美激情| 亚洲精品,欧美精品| xxxhd国产人妻xxx| 精品99又大又爽又粗少妇毛片| videosex国产| 一级毛片 在线播放| 免费日韩欧美在线观看| 精品视频人人做人人爽| 午夜日本视频在线| 黄片播放在线免费| 狠狠精品人妻久久久久久综合| 亚洲av二区三区四区| 国产一区二区在线观看av| 亚洲综合精品二区| 丝袜喷水一区| 我的老师免费观看完整版| 女人精品久久久久毛片| 久久影院123| 亚洲精品中文字幕在线视频| 国产一区有黄有色的免费视频| 99国产综合亚洲精品| 国产成人精品婷婷| 99热全是精品| 岛国毛片在线播放| 九九久久精品国产亚洲av麻豆| 夜夜骑夜夜射夜夜干| 亚洲丝袜综合中文字幕| 蜜臀久久99精品久久宅男| xxxhd国产人妻xxx| 久久久久久久久久久丰满| 亚洲欧美一区二区三区黑人 | 亚洲欧美日韩另类电影网站| 99re6热这里在线精品视频| 最新中文字幕久久久久| 亚洲欧美清纯卡通| 另类亚洲欧美激情| 日韩中字成人| 夜夜骑夜夜射夜夜干| 赤兔流量卡办理| 看非洲黑人一级黄片| 婷婷色综合www| 国产精品国产av在线观看| 国产一区二区在线观看日韩| 最后的刺客免费高清国语| 久久这里有精品视频免费| 亚洲色图 男人天堂 中文字幕 | av.在线天堂| 久久午夜综合久久蜜桃| 亚洲美女搞黄在线观看| 国产成人精品无人区| 日韩成人av中文字幕在线观看| 亚洲国产精品国产精品| 男的添女的下面高潮视频| 久久人人爽av亚洲精品天堂| 亚洲国产精品国产精品| 91午夜精品亚洲一区二区三区| 91在线精品国自产拍蜜月| 国产一区有黄有色的免费视频| av在线观看视频网站免费| 好男人视频免费观看在线| 精品人妻在线不人妻| 人人妻人人澡人人看| 国产精品99久久99久久久不卡 | 国产一级毛片在线| 成年美女黄网站色视频大全免费 | 丝袜喷水一区| 一区二区三区精品91| 亚洲综合精品二区| 欧美人与善性xxx| xxx大片免费视频| 晚上一个人看的免费电影| 国产成人精品无人区| 久久97久久精品| 秋霞伦理黄片| 国产极品粉嫩免费观看在线 | 晚上一个人看的免费电影| av卡一久久| 高清午夜精品一区二区三区| 亚洲欧美中文字幕日韩二区| 免费人妻精品一区二区三区视频| 久久久久久久国产电影| 黑人高潮一二区| 亚洲精华国产精华液的使用体验| 日韩成人av中文字幕在线观看| 免费高清在线观看日韩| 男人爽女人下面视频在线观看| 亚洲精品亚洲一区二区| 国产精品蜜桃在线观看| 成人国语在线视频| 精品人妻在线不人妻| 久久久久久久精品精品| 国产成人精品在线电影| av电影中文网址| 亚洲精品aⅴ在线观看| 精品久久久精品久久久| 夜夜爽夜夜爽视频| 日本wwww免费看| 欧美丝袜亚洲另类| 男的添女的下面高潮视频| 精品国产国语对白av| 久久精品国产亚洲av涩爱| 国产精品免费大片| 中文字幕人妻熟人妻熟丝袜美| 日本色播在线视频| 少妇 在线观看| 大话2 男鬼变身卡| av国产久精品久网站免费入址| 欧美激情国产日韩精品一区| 多毛熟女@视频| 午夜影院在线不卡| 亚洲伊人久久精品综合| 久热久热在线精品观看| 亚洲精品视频女| 国产老妇伦熟女老妇高清| 免费av不卡在线播放| 啦啦啦中文免费视频观看日本| 国产又色又爽无遮挡免| 中文精品一卡2卡3卡4更新| 97在线人人人人妻| 久久精品国产自在天天线| 丰满饥渴人妻一区二区三| 人成视频在线观看免费观看| 夜夜骑夜夜射夜夜干| 国产男女超爽视频在线观看| 亚洲色图综合在线观看| 欧美+日韩+精品| 久久久久人妻精品一区果冻| h视频一区二区三区| 啦啦啦在线观看免费高清www| 老女人水多毛片| 亚洲综合色惰|