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

    翼傘充氣過程的流固耦合方法數(shù)值仿真

    2020-07-31 09:48:48張思宇余莉劉鑫
    關(guān)鍵詞:傘衣氣室前緣

    張思宇,余莉,,*,劉鑫

    (1.南京航空航天大學(xué) 飛行器環(huán)境控制與生命保障工業(yè)和信息化部重點(diǎn)實(shí)驗(yàn)室,南京210016;2.南京航空航天大學(xué) 航空學(xué)院,南京210016)

    沖壓式翼傘相比于傳統(tǒng)降落傘,升阻比高、滑翔穩(wěn)定性能和操作性良好,廣泛應(yīng)用于航空航天等領(lǐng)域[1-3],近年來成為各國(guó)研究重點(diǎn)。翼傘充氣成功與否決定著翼傘系統(tǒng)工作的成敗。

    翼傘和常規(guī)降落傘一樣,在充氣展開過程中,由于本身的柔性特性,在氣流作用下柔性傘衣極易變形,工作特性發(fā)生改變,涉及到的流固耦合問題具有高度非線性、非定常特點(diǎn),理論研究難度很大。隨著計(jì)算機(jī)的發(fā)展和數(shù)值精度的提高,數(shù)值模擬逐漸成為研究降落傘流固耦合性能的重要手段。Purvis[4]、潘星和曹義華[5]建立了降落傘結(jié)構(gòu)的彈簧阻尼多節(jié)點(diǎn)模型,與內(nèi)部簡(jiǎn)化流場(chǎng)進(jìn)行耦合計(jì)算?;跁r(shí)空離散的MSD-CFD(Mass Spring Damper-Computational Fluid Dynamics)耦合方法[6-8]側(cè)重于傘衣外形的計(jì)算,無法得到結(jié)構(gòu)應(yīng)力信息。Takizawa[9]、Stein[10]和Tezduyar[11]等將DSD/SST(Deform ing-Spatial-Domain/Stabilized Space-Time)有限元計(jì)算方法用于三維降落傘流固耦合仿真計(jì)算研究中,該方法側(cè)重于外部流場(chǎng)的數(shù)值計(jì)算,對(duì)初始充氣段缺少研究。Kim和Peskin[12]采用浸入邊界法進(jìn)行降落傘充氣過程的流固耦合模擬。Tutt[13]和Yu[14]等采用任意拉格朗日-歐拉(Arbitrary Lagrange-Euler,ALE)方法模擬了三維降落傘開傘過程。

    但是,目前降落傘的流固耦合研究大多集中于平面圓形傘,針對(duì)翼傘充氣過程的三維流固耦合數(shù)值研究一直沒有突破性進(jìn)展,常常做出一定簡(jiǎn)化,主要面向剛性翼傘[15-16]或者基于松散耦合方法進(jìn)行穩(wěn)態(tài)定常計(jì)算。例如,F(xiàn)ogell等[17]采用弱耦合方法分析了翼傘在穩(wěn)定滑翔階段的氣動(dòng)外形。Kalro等[18]采用并行有限元計(jì)算對(duì)翼傘初始充氣階段進(jìn)行了數(shù)值模擬。汪龍芳和賀衛(wèi)亮[19]對(duì)定常狀況下翼傘的流固耦合變形問題進(jìn)行了三維數(shù)值模擬。張春和曹義華[20]基于弱耦合方法分析了翼傘的氣動(dòng)變形。

    翼傘充氣過程大變形、大位移、多尺度流動(dòng)的非定常流固耦合研究基本處于空白階段。這主要是由于:翼傘為非展平曲面構(gòu)成的多氣室-雙翼面復(fù)雜結(jié)構(gòu),之間由開孔翼肋隔成若干連通氣室,前緣充氣,后緣封閉,結(jié)構(gòu)的復(fù)雜性導(dǎo)致其無法采用常規(guī)方法建立折疊模型。此外,多氣室翼傘在充氣過程中多方向同時(shí)展開,是一個(gè)典型的非穩(wěn)態(tài)流場(chǎng)-柔性傘衣結(jié)構(gòu)的劇烈耦合問題,規(guī)律復(fù)雜。上述問題作為翼傘研究領(lǐng)域的難點(diǎn),一直備受學(xué)者關(guān)注。

    為了解決上述問題,本文提出了基于自由曲面變形(Free Form Deformation,F(xiàn)FD)理論的多氣室柔性沖壓翼傘展向折疊建模方法,流體域通過時(shí)步更新實(shí)現(xiàn)隨傘載系統(tǒng)運(yùn)動(dòng),基于ALE方法開展折疊翼傘充氣過程的非定常流固耦合數(shù)值計(jì)算,分析翼傘的非線性動(dòng)力學(xué)行為,系統(tǒng)地探究非穩(wěn)態(tài)流場(chǎng)-柔性傘衣結(jié)構(gòu)之間的耦合機(jī)理,為翼傘安全設(shè)計(jì)提供理論依據(jù)。

    1 計(jì)算模型

    1.1 自由曲面變形理論

    翼傘為非展平曲面所構(gòu)成的多氣室-雙翼面結(jié)構(gòu),為解決其展向折疊問題,本文根據(jù)FFD理論[21],提出了多氣室柔性沖壓翼傘的展向折疊建模方法。將待變形的翼傘模型嵌入FFD控制晶格中,通過移動(dòng)晶格節(jié)點(diǎn),將晶格變形傳遞給嵌入其中的翼傘模型,達(dá)到對(duì)其折疊建模的目的。

    本文采用六面體控制晶格,如圖1所示,建立局部坐標(biāo)系WSH,W、S、H是局部坐標(biāo)系軸矢量,沿坐標(biāo)軸分別均勻布置l、m、n個(gè)控制節(jié)點(diǎn),局部坐標(biāo)系原點(diǎn)在全局坐標(biāo)系下的坐標(biāo)記為Q0。

    FFD控制晶格各節(jié)點(diǎn)的全局坐標(biāo)為

    圖1 FFD六面體控制晶格Fig.1 FFD hexahedral control lattice

    式中:(i,j,k)為控制晶格各節(jié)點(diǎn)的局部標(biāo)識(shí),i=1,2,…,l;j=1,2,…,m;k=1,2,…,n。

    在FFD控制晶格中,翼傘任意一點(diǎn)的局部坐標(biāo)是(w,s,h),0≤w,s,h≤1,有

    該點(diǎn)的全局坐標(biāo)為Q=Q0+w W+s S+h H。本文采用Bernstein基函數(shù)建立控制晶格各節(jié)點(diǎn)位移ΔPi,j,k與ΔQ之間的映射函數(shù)關(guān)系:

    1.2 流固耦合控制方程

    流場(chǎng)域控制方程為

    流場(chǎng)域采用無反射邊界條件,以消除流場(chǎng)邊界處的單元速度反射波干擾,防止邊界對(duì)計(jì)算產(chǎn)生影響。

    翼傘傘衣、傘繩和流場(chǎng)發(fā)生耦合,采用拉格朗日方法計(jì)算,控制方程為

    ALE網(wǎng)格運(yùn)動(dòng)方程為

    式中:Li為拉格朗日坐標(biāo)。

    對(duì)控制方程式(4)~式(6)進(jìn)行全耦合計(jì)算,采用中心差分時(shí)間顯式方法求解。對(duì)于每一個(gè)節(jié)點(diǎn),流場(chǎng)和結(jié)構(gòu)的位移x、速度u按式(7)進(jìn)行更新:

    式中:Fint、Fext分別為內(nèi)、外力矢量;M 為質(zhì)量對(duì)角矩陣。

    1.3 流場(chǎng)域時(shí)步更新

    翼傘充氣過程是有限質(zhì)量充氣,為減小計(jì)算消耗,本文在每一時(shí)間步都對(duì)流場(chǎng)域進(jìn)行更新,使其跟隨載荷對(duì)象一起進(jìn)行空間運(yùn)動(dòng)。首先在結(jié)構(gòu)網(wǎng)格上任意選擇不共線的3個(gè)節(jié)點(diǎn)A、B、C,坐標(biāo)分別為KA、KB和KC,建立一個(gè)以A為原點(diǎn)的局部坐標(biāo)系(見圖2),其3個(gè)軸矢量分別為

    圖2 流場(chǎng)域時(shí)間步更新Fig.2 Time-step update of flow field

    每個(gè)時(shí)間步后,局部坐標(biāo)將跟隨結(jié)構(gòu)網(wǎng)格發(fā)生位移,根據(jù)局部坐標(biāo)在每個(gè)時(shí)間步前后的位移獲得變換矩陣T,則流場(chǎng)網(wǎng)格新的坐標(biāo)系為

    式中:[e*1e*2e*31]為更新后的齊次坐標(biāo);[e1e2e31]為更新前的齊次坐標(biāo)。流體域通過時(shí)間步更新技術(shù)實(shí)現(xiàn)運(yùn)動(dòng)與重構(gòu)。

    2 研究對(duì)象與模型驗(yàn)證

    2.1 研究對(duì)象與計(jì)算工況

    本文采用文獻(xiàn)[22]的沖壓式翼傘模型,如圖3所示,翼梢兩側(cè)有穩(wěn)定幅,32根傘繩,平鋪展長(zhǎng)5.48m,弦長(zhǎng)2.74m,繩長(zhǎng)3.36m,有7個(gè)氣室,每個(gè)氣室由一個(gè)非承載肋片分割成2個(gè)半氣室。翼傘的折疊模型如圖4所示。

    本文采用ALE流固耦合數(shù)值方法,將結(jié)構(gòu)網(wǎng)格單元穿插于流場(chǎng)網(wǎng)格中,其流固耦合數(shù)值計(jì)算模型如圖5所示(規(guī)模為6WS×8WS×10WS,其中WS為傘衣弦長(zhǎng)),其中結(jié)構(gòu)附近流場(chǎng)進(jìn)行加密處理,參數(shù)設(shè)置如表1所示。翼傘以0°迎角充氣,初始充氣速度為17m/s。本文基于ANSYS平臺(tái)計(jì)算,采用8核處理器,仿真0.6 s耗時(shí)約172 h。

    圖3 沖壓式翼傘結(jié)構(gòu)示意圖Fig.3 Schematic diagram of rammed parafoil structure

    圖5 流固耦合數(shù)值計(jì)算模型Fig.5 Numerical calculation model of fluid-structure interaction

    表1 流固耦合過程的數(shù)值模型信息Table 1 Numerical model information of fluid-structure interaction process

    2.2 模型驗(yàn)證

    為驗(yàn)證本文模型的合理性和準(zhǔn)確性,采用文獻(xiàn)[22]的空投試驗(yàn)數(shù)據(jù)進(jìn)行驗(yàn)證計(jì)算,同時(shí)進(jìn)行流場(chǎng)網(wǎng)格的無關(guān)性驗(yàn)證。本文采用3種不同網(wǎng)格密度的流場(chǎng)進(jìn)行數(shù)值計(jì)算,對(duì)比結(jié)果如表2所示。

    隨網(wǎng)格數(shù)量增加,計(jì)算誤差減小,但消耗時(shí)長(zhǎng)增加,因此綜合考慮,本文選擇網(wǎng)格數(shù)量為849 000的模型進(jìn)行數(shù)值計(jì)算,得到的翼傘充滿外形和開傘動(dòng)載F(充氣過程中作用在傘衣上的縱向載荷)與空投試驗(yàn)的對(duì)比情況如圖6、圖7所示??梢钥闯?,翼傘充滿外形與試驗(yàn)外形基本一致。本文開傘動(dòng)載計(jì)算結(jié)果的變化趨勢(shì)與空投試驗(yàn)變化規(guī)律相近,數(shù)值計(jì)算得到的開傘動(dòng)載峰值出現(xiàn)時(shí)刻略早于文獻(xiàn)[22],數(shù)值為4 760.1 N,和文獻(xiàn)[22]試驗(yàn)數(shù)據(jù)(4 963.2 N)比較,誤差為4.09%,該誤差滿足柔性翼傘開傘計(jì)算的誤差要求。綜上認(rèn)為,本文基于多氣室雙翼面沖壓翼傘建立的數(shù)值計(jì)算模型可以有效模擬翼傘充氣階段工作過程。

    表2 網(wǎng)格無關(guān)性對(duì)比結(jié)果Table 2 Comparison results of grid independence

    圖6 數(shù)值計(jì)算結(jié)果與空投試驗(yàn)充滿外形對(duì)比Fig.6 Comparison of canopy shape between numerical calculation result and airdrop test

    圖7 開傘動(dòng)載隨時(shí)間的變化曲線Fig.7 Curves of opening load over time

    3 研究結(jié)果

    3.1 傘衣外形及流場(chǎng)變化

    翼傘充氣展開過程為有限質(zhì)量充氣情況,許多隨機(jī)參數(shù)都會(huì)對(duì)充氣性能產(chǎn)生影響,本文未考慮風(fēng)場(chǎng)等隨機(jī)環(huán)境因素對(duì)充氣過程的影響。圖8為翼傘充氣過程的外形變化。在充氣初期,翼傘前緣切口進(jìn)氣,尾部最先膨脹展開。之后下翼面高壓氣流作用于穩(wěn)定幅,拉動(dòng)傘衣展向伸展,穩(wěn)定幅產(chǎn)生應(yīng)力集中現(xiàn)象。充氣過程中,前緣切口約束較少,且在氣室充滿后由于切口角度的存在,始終承受氣流吹襲,震顫頻率大,產(chǎn)生持續(xù)的應(yīng)力集中現(xiàn)象。因此,在翼傘設(shè)計(jì)中應(yīng)加強(qiáng)穩(wěn)定幅與前緣切口處的材料強(qiáng)度,避免傘衣破損。

    為進(jìn)一步探究翼傘充氣過程中的結(jié)構(gòu)及流場(chǎng)的非定常變化規(guī)律,選取a、b、c 3個(gè)截面(見圖9)進(jìn)行研究,各個(gè)截面的結(jié)構(gòu)及流場(chǎng)變化如圖10所示(各個(gè)時(shí)刻從上到下依次為a、b、c截面),為刻畫細(xì)節(jié)情況,各個(gè)時(shí)刻的圖形縮比尺寸不同。

    圖8 充氣過程翼傘外形變化(范式等效應(yīng)力云圖)Fig.8 Parafoil shape changes during inflation process(von M ises equivalent stress contours)

    圖9 選取截面示意圖Fig.9 Schematic diagram of selected sections

    首先,通過觀察圖10所示翼傘充氣過程的弦向變化(a、b截面):在翼傘初始折疊狀態(tài)時(shí),a、b截面的速度場(chǎng)分布相似,前緣產(chǎn)生渦。隨翼傘逐漸充氣,上翼面的渦逐漸后移至脫體,前緣處不斷產(chǎn)生新的渦。在充氣過程中,a截面隨氣室充氣逐漸形成飽滿翼型,內(nèi)腔壓力增加,傘衣剛度增加,類似剛性翼,外部流暢穩(wěn)定。而b截面處于最外側(cè),傘衣由于約束較少變形明顯,在翼傘充滿后翼型飽滿度不足,外部流場(chǎng)不穩(wěn)定,傘衣變形狀態(tài)具有一定的隨機(jī)性。

    圖10中翼傘充氣過程的展向變化(c截面)表明:在翼傘初始折疊狀態(tài)時(shí),c截面周圍流場(chǎng)產(chǎn)生對(duì)稱繞流現(xiàn)象,此時(shí)傘衣折疊壓縮率大,下翼面形成高壓區(qū)(壓力云圖見圖11)。之后翼傘底部對(duì)稱渦開始破裂,氣流繞翼尖流出,翼傘展開充氣。該過程上翼面流速快,下翼面流速慢,上下翼面形成的壓力造成“翼尖上翹,中部凹陷”的翼傘尾流再附現(xiàn)象。因此在翼傘開傘過程中應(yīng)加入引導(dǎo)傘改善充氣狀況,引導(dǎo)傘拉出傘衣中部減輕氣室塌陷。此外,設(shè)計(jì)合理的穩(wěn)定幅,阻擋下翼面到上翼面的繞流,消弱翼尖渦強(qiáng)度,減輕干擾影響。此后,傘衣隨氣室充滿保持穩(wěn)定充滿翼型。

    圖11 充氣過程壓力云圖Fig.11 Pressure contours of inflation process

    3.2 氣室充氣規(guī)律

    為進(jìn)一步研究翼傘的各氣室充氣規(guī)律,將翼傘沿弦向分為前(前緣)、中(中部)、后(尾緣)三部分,沿展向?qū)?個(gè)氣室逐次編號(hào),兩側(cè)為氣室1和氣室7,中央氣室為4號(hào)。圖12為各氣室充滿時(shí)間和充滿寬度的變化規(guī)律??梢钥闯觯?/p>

    1)各氣室的充氣規(guī)律關(guān)于中央氣室具有對(duì)稱性,中央氣室充氣快且充滿外形飽滿,越靠近兩側(cè)充滿速率越慢,充氣效果越差。

    2)在高壓氣流推動(dòng)作用下,進(jìn)入氣室的氣體涌入尾部,翼傘尾緣最先展向打開,充氣速率明顯快于其他部分,前緣展開最慢。充滿后中部充氣效果最好,氣室最為飽滿。

    圖12 各氣室充滿時(shí)間和充滿寬度Fig.12 Inflation time and width of each air chamber

    3)由于尾緣上下翼面縫合,翼尖繞流導(dǎo)致傘衣變形,內(nèi)外氣壓不穩(wěn)定,兩側(cè)氣室的后部未完全充滿。

    充氣過程的翼展變化見圖13。翼傘的前緣、中部、尾緣前、中、后三部分于0.39、0.37、0.25 s分別達(dá)到翼展峰值5.05、5.40、5.20m。之后由于傘繩和傘衣的彈性作用,各部位充滿后產(chǎn)生一定的回彈和氣室“鼓包”現(xiàn)象,中部充滿翼展穩(wěn)定在5.0m左右,略小于設(shè)計(jì)值。

    圖13 充氣過程的翼展變化Fig.13 Changes in wingspan during inflation process

    3.3 氣動(dòng)特性分析

    翼傘的氣動(dòng)特性參數(shù)變化情況如圖14所示。0.1 s之前,計(jì)算域尚未形成穩(wěn)定繞流流場(chǎng),氣動(dòng)參數(shù)不具有參考意義。隨繞流流場(chǎng)形成,翼傘的升阻力特性參數(shù)變化趨勢(shì)一致,在0.28 s時(shí)達(dá)到峰值,此時(shí)傘衣投影面積最大。之后由于呼吸回彈現(xiàn)象,傘衣面積減小,氣動(dòng)系數(shù)減小,翼傘充滿穩(wěn)定后氣動(dòng)系數(shù)略有回升。充滿后傘載系統(tǒng)存在一定迎角,滑翔比CL/CD穩(wěn)定在2.24,升力系數(shù)CL為1.56,阻力系數(shù)CD為0.70。

    圖14 氣動(dòng)特性參數(shù)的變化Fig.14 Changes in aerodynamic characteristic parameter

    4 結(jié) 論

    本文針對(duì)非定常充氣展開過程,數(shù)值模擬了翼傘的三維流固耦合動(dòng)力學(xué)特性,得出:

    1)翼傘充氣展開過程中穩(wěn)定幅和前緣切口處產(chǎn)生應(yīng)力集中現(xiàn)象,因此,在翼傘設(shè)計(jì)中應(yīng)加強(qiáng)穩(wěn)定幅與前緣切口處的材料強(qiáng)度。

    2)翼傘充氣展開過程存在“翼尖上翹,中部凹陷”的翼傘尾流再附現(xiàn)象,造成傘衣內(nèi)塌,導(dǎo)致系統(tǒng)穩(wěn)定性受影響,因此在翼傘開傘過程中應(yīng)加入引導(dǎo)傘以改善充氣狀況,設(shè)計(jì)合理的穩(wěn)定幅可以減小翼尖渦的影響。

    3)各氣室的充氣規(guī)律關(guān)于中央氣室具有對(duì)稱性,中央氣室充氣快且充滿外形飽滿,兩側(cè)充滿速率越慢,氣室飽滿度下降;翼傘尾部最先膨脹展開,前緣展開最慢,充滿后存在一定的回彈和氣室“鼓包”現(xiàn)象,充滿后展長(zhǎng)小于設(shè)計(jì)值。

    4)升阻力特性參數(shù)變化趨勢(shì)一致,傘衣充滿后翼傘滑翔比穩(wěn)定在2.24,升力系數(shù)為1.56,阻力系數(shù)為0.70。

    猜你喜歡
    傘衣氣室前緣
    小型降落傘充氣和傘衣塌陷過程
    基于Hyperworks的重卡氣室支架結(jié)構(gòu)優(yōu)化設(shè)計(jì)研究
    設(shè)計(jì)參數(shù)及大氣參數(shù)對(duì)降落傘充氣性能的影響
    一種飛機(jī)尾翼前緣除冰套安裝方式
    某重型特種車制動(dòng)氣室支架優(yōu)化設(shè)計(jì)
    汽車零部件(2019年5期)2019-06-13 06:25:50
    傘衣透氣性對(duì)翼傘氣動(dòng)特性的影響
    一起220kV GIS設(shè)備斷路器氣室罐體放電原因分析
    廣西電力(2016年5期)2016-07-10 09:16:44
    氣室裝置對(duì)變威力發(fā)射影響分析
    深水沉積研究進(jìn)展及前緣問題
    成人无遮挡网站| kizo精华| 国产精品久久视频播放| 最近视频中文字幕2019在线8| 精华霜和精华液先用哪个| 男女边摸边吃奶| 永久网站在线| 欧美日韩亚洲高清精品| 女人十人毛片免费观看3o分钟| 日本wwww免费看| 亚洲欧美一区二区三区黑人 | 十八禁国产超污无遮挡网站| 免费大片18禁| 天天躁夜夜躁狠狠久久av| 精品酒店卫生间| 天堂av国产一区二区熟女人妻| 国产一级毛片在线| 一本一本综合久久| 亚洲婷婷狠狠爱综合网| 乱人视频在线观看| 国产亚洲5aaaaa淫片| 国产午夜精品论理片| 亚洲精品456在线播放app| 久久国内精品自在自线图片| 亚洲最大成人av| 久久精品人妻少妇| 日韩成人伦理影院| 乱系列少妇在线播放| 乱人视频在线观看| 亚洲乱码一区二区免费版| 三级经典国产精品| 久久久久久久亚洲中文字幕| 精品久久久精品久久久| 国产亚洲午夜精品一区二区久久 | 国产精品99久久久久久久久| 99久国产av精品| 毛片一级片免费看久久久久| 五月伊人婷婷丁香| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av在线观看美女高潮| 成年人午夜在线观看视频 | 午夜老司机福利剧场| 黑人高潮一二区| 日韩欧美 国产精品| av在线亚洲专区| 观看免费一级毛片| 国产在视频线精品| 99九九线精品视频在线观看视频| 91aial.com中文字幕在线观看| 网址你懂的国产日韩在线| 国产午夜精品一二区理论片| 成年女人在线观看亚洲视频 | 青青草视频在线视频观看| 你懂的网址亚洲精品在线观看| 亚洲丝袜综合中文字幕| 九九爱精品视频在线观看| 韩国av在线不卡| 韩国高清视频一区二区三区| 成人鲁丝片一二三区免费| 午夜福利在线观看吧| 九色成人免费人妻av| av免费观看日本| 欧美激情久久久久久爽电影| www.av在线官网国产| 国产伦理片在线播放av一区| 成人美女网站在线观看视频| 亚洲av不卡在线观看| 国产片特级美女逼逼视频| 国产淫片久久久久久久久| 少妇裸体淫交视频免费看高清| 菩萨蛮人人尽说江南好唐韦庄| 女人久久www免费人成看片| 中文在线观看免费www的网站| 色视频www国产| 婷婷色综合www| 亚洲精品国产av成人精品| 在线观看美女被高潮喷水网站| 国产不卡一卡二| 国产探花极品一区二区| 麻豆精品久久久久久蜜桃| 精品人妻熟女av久视频| 又爽又黄无遮挡网站| 日日啪夜夜爽| 亚洲国产精品国产精品| 久久久久久久国产电影| 国产视频内射| 久久草成人影院| 免费看av在线观看网站| 免费观看无遮挡的男女| 国产精品国产三级专区第一集| 搞女人的毛片| 91精品一卡2卡3卡4卡| 亚洲精品亚洲一区二区| 日韩不卡一区二区三区视频在线| 亚洲精品一区蜜桃| 国产精品蜜桃在线观看| 精品久久久久久成人av| 亚洲不卡免费看| 中文字幕人妻熟人妻熟丝袜美| 国产成人精品一,二区| 青春草国产在线视频| 亚洲最大成人手机在线| a级毛片免费高清观看在线播放| eeuss影院久久| 欧美日韩一区二区视频在线观看视频在线 | 内地一区二区视频在线| 卡戴珊不雅视频在线播放| 国产亚洲最大av| av免费观看日本| 久久精品人妻少妇| 免费av观看视频| 欧美高清性xxxxhd video| 精品欧美国产一区二区三| 深爱激情五月婷婷| 亚州av有码| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最新中文字幕久久久久| 九草在线视频观看| av在线亚洲专区| 日韩欧美一区视频在线观看 | 久久午夜福利片| 午夜激情福利司机影院| 深爱激情五月婷婷| 97在线视频观看| 亚洲av.av天堂| 成人av在线播放网站| 日本一二三区视频观看| 最近中文字幕高清免费大全6| 一级毛片久久久久久久久女| av一本久久久久| 夫妻午夜视频| 亚洲av福利一区| 亚洲自偷自拍三级| 搡女人真爽免费视频火全软件| 婷婷色综合大香蕉| 亚洲精品自拍成人| 精品人妻偷拍中文字幕| 老司机影院毛片| 九草在线视频观看| 色播亚洲综合网| a级毛片免费高清观看在线播放| 爱豆传媒免费全集在线观看| 国产精品久久久久久精品电影| 神马国产精品三级电影在线观看| 亚洲av日韩在线播放| 亚洲,欧美,日韩| 国产成人精品一,二区| 麻豆精品久久久久久蜜桃| 免费观看的影片在线观看| 69人妻影院| 色哟哟·www| 成人国产麻豆网| 午夜激情欧美在线| 97在线视频观看| 中文乱码字字幕精品一区二区三区 | 肉色欧美久久久久久久蜜桃 | 亚洲国产精品成人久久小说| 一夜夜www| 激情五月婷婷亚洲| 成人国产麻豆网| 男人舔奶头视频| 成年女人看的毛片在线观看| 久久久久免费精品人妻一区二区| 日韩电影二区| 天堂影院成人在线观看| 好男人在线观看高清免费视频| 最近视频中文字幕2019在线8| 国产男人的电影天堂91| 搡女人真爽免费视频火全软件| 日韩中字成人| 亚洲成人中文字幕在线播放| 男女边摸边吃奶| 成年女人看的毛片在线观看| 纵有疾风起免费观看全集完整版 | 国产探花在线观看一区二区| 晚上一个人看的免费电影| 能在线免费看毛片的网站| 综合色丁香网| 免费不卡的大黄色大毛片视频在线观看 | 天天一区二区日本电影三级| 日韩欧美一区视频在线观看 | 亚洲成人一二三区av| 国产亚洲一区二区精品| 麻豆成人午夜福利视频| 亚洲精品成人av观看孕妇| 久久精品夜色国产| 中文字幕免费在线视频6| 乱码一卡2卡4卡精品| 亚洲av中文字字幕乱码综合| 岛国毛片在线播放| av.在线天堂| 你懂的网址亚洲精品在线观看| 日韩欧美 国产精品| 亚洲第一区二区三区不卡| 精品国产一区二区三区久久久樱花 | 日韩一区二区三区影片| 亚洲欧美中文字幕日韩二区| 日本欧美国产在线视频| 麻豆国产97在线/欧美| 婷婷色麻豆天堂久久| 亚洲av中文av极速乱| 两个人的视频大全免费| 亚洲国产高清在线一区二区三| 中文天堂在线官网| 国产精品福利在线免费观看| 搡女人真爽免费视频火全软件| 又大又黄又爽视频免费| 国产av在哪里看| 亚洲国产精品sss在线观看| 精品一区二区免费观看| 亚洲欧美成人综合另类久久久| 高清日韩中文字幕在线| 草草在线视频免费看| 日本wwww免费看| 亚洲国产精品专区欧美| 婷婷六月久久综合丁香| 国产一区亚洲一区在线观看| av免费在线看不卡| 99热网站在线观看| 国产成人精品婷婷| 91久久精品国产一区二区三区| 亚洲av福利一区| 国产亚洲91精品色在线| 一级毛片电影观看| 国产精品日韩av在线免费观看| 国产伦精品一区二区三区四那| 欧美变态另类bdsm刘玥| 国产激情偷乱视频一区二区| 中文资源天堂在线| 免费看不卡的av| 欧美+日韩+精品| 久久久久久久久久成人| 九草在线视频观看| 亚洲av成人精品一区久久| 亚洲av电影不卡..在线观看| 亚洲欧美中文字幕日韩二区| 婷婷色av中文字幕| av黄色大香蕉| 赤兔流量卡办理| 国产毛片a区久久久久| 午夜视频国产福利| 伦精品一区二区三区| 看十八女毛片水多多多| 你懂的网址亚洲精品在线观看| 永久免费av网站大全| 国产黄色免费在线视频| 国产不卡一卡二| 精品久久久久久电影网| 高清日韩中文字幕在线| 久久久久久伊人网av| 寂寞人妻少妇视频99o| 国产女主播在线喷水免费视频网站 | 直男gayav资源| 亚洲人成网站在线观看播放| kizo精华| 国产男女超爽视频在线观看| 亚洲欧美精品专区久久| 欧美日韩精品成人综合77777| 国产三级在线视频| 2021天堂中文幕一二区在线观| 国产成人福利小说| 国产精品久久久久久av不卡| 精品久久久精品久久久| 丰满少妇做爰视频| 韩国av在线不卡| 亚洲精品中文字幕在线视频 | 精品一区在线观看国产| 亚洲av福利一区| 91精品一卡2卡3卡4卡| 久久久久久久久久人人人人人人| 男女啪啪激烈高潮av片| 久久久成人免费电影| 乱人视频在线观看| 少妇熟女欧美另类| 老师上课跳d突然被开到最大视频| 伦精品一区二区三区| 欧美日韩综合久久久久久| 哪个播放器可以免费观看大片| 狂野欧美激情性xxxx在线观看| 久久久久网色| av黄色大香蕉| 免费看光身美女| 久久久久久久久久久免费av| 国产高潮美女av| 午夜激情久久久久久久| 亚洲av一区综合| 卡戴珊不雅视频在线播放| 麻豆成人午夜福利视频| 高清毛片免费看| 亚洲最大成人中文| 精品熟女少妇av免费看| 久久久久精品久久久久真实原创| 中国国产av一级| 激情 狠狠 欧美| av卡一久久| 婷婷色av中文字幕| 国产成人午夜福利电影在线观看| 夫妻性生交免费视频一级片| 欧美xxⅹ黑人| 全区人妻精品视频| 一个人看视频在线观看www免费| 国产爱豆传媒在线观看| 中文精品一卡2卡3卡4更新| 69人妻影院| 亚洲欧美一区二区三区国产| 亚洲国产精品国产精品| 精品一区二区三卡| 男的添女的下面高潮视频| 又粗又硬又长又爽又黄的视频| 美女内射精品一级片tv| 男人狂女人下面高潮的视频| 搞女人的毛片| 久久99精品国语久久久| 最近视频中文字幕2019在线8| 久久人人爽人人片av| 久久久久免费精品人妻一区二区| 中文欧美无线码| 亚洲精品国产成人久久av| 日日干狠狠操夜夜爽| 精品一区二区三区视频在线| 亚洲精品一区蜜桃| 国产午夜精品一二区理论片| 欧美日韩综合久久久久久| 国产真实伦视频高清在线观看| 国产成人福利小说| 观看美女的网站| 亚洲国产精品成人久久小说| 久久久久久久久久久丰满| 97超碰精品成人国产| 国产在线男女| www.av在线官网国产| av网站免费在线观看视频 | 九九在线视频观看精品| 亚洲真实伦在线观看| 丰满乱子伦码专区| 国内精品宾馆在线| 免费黄网站久久成人精品| 全区人妻精品视频| 国产伦精品一区二区三区视频9| 两个人视频免费观看高清| 麻豆av噜噜一区二区三区| 高清毛片免费看| 波多野结衣巨乳人妻| 大话2 男鬼变身卡| 99久久中文字幕三级久久日本| 国产精品精品国产色婷婷| 欧美日韩国产mv在线观看视频 | 人妻制服诱惑在线中文字幕| 纵有疾风起免费观看全集完整版 | 欧美日韩国产mv在线观看视频 | 日韩欧美 国产精品| 亚洲欧洲国产日韩| 精品久久久久久久久av| 免费看美女性在线毛片视频| 非洲黑人性xxxx精品又粗又长| 久热久热在线精品观看| 国产日韩欧美在线精品| 国产男女超爽视频在线观看| 少妇人妻一区二区三区视频| 天美传媒精品一区二区| 免费在线观看成人毛片| 日韩在线高清观看一区二区三区| 日韩精品青青久久久久久| 亚洲av电影在线观看一区二区三区 | 国产精品.久久久| 国产高潮美女av| 汤姆久久久久久久影院中文字幕 | 国产一区二区亚洲精品在线观看| 一级毛片我不卡| 熟妇人妻不卡中文字幕| 黄色配什么色好看| 少妇丰满av| 卡戴珊不雅视频在线播放| 超碰97精品在线观看| 午夜久久久久精精品| 看十八女毛片水多多多| eeuss影院久久| 搡老乐熟女国产| 水蜜桃什么品种好| 国产又色又爽无遮挡免| 久久精品人妻少妇| 亚洲精品乱久久久久久| 久久97久久精品| 免费看av在线观看网站| 亚洲在线自拍视频| 免费高清在线观看视频在线观看| 国产综合精华液| 午夜福利成人在线免费观看| 国产伦在线观看视频一区| 亚洲精品亚洲一区二区| 激情 狠狠 欧美| 亚洲成人一二三区av| 在线a可以看的网站| 高清日韩中文字幕在线| 亚洲婷婷狠狠爱综合网| 亚洲久久久久久中文字幕| 亚洲欧美成人精品一区二区| 欧美三级亚洲精品| 日韩中字成人| 亚洲av电影在线观看一区二区三区 | 国产老妇女一区| 久久久久精品久久久久真实原创| 亚洲国产精品成人综合色| 尤物成人国产欧美一区二区三区| ponron亚洲| 日本猛色少妇xxxxx猛交久久| 99久久精品国产国产毛片| 亚洲av成人精品一区久久| 精品熟女少妇av免费看| 亚洲精品色激情综合| 毛片女人毛片| 大又大粗又爽又黄少妇毛片口| 国产免费一级a男人的天堂| 亚洲欧洲日产国产| 亚洲国产精品成人久久小说| 国产淫语在线视频| 日韩av在线免费看完整版不卡| 狂野欧美激情性xxxx在线观看| 女人被狂操c到高潮| 久久精品国产亚洲av涩爱| 男女边摸边吃奶| 国产国拍精品亚洲av在线观看| 国产麻豆成人av免费视频| 熟妇人妻久久中文字幕3abv| 国产精品国产三级国产av玫瑰| 天天躁夜夜躁狠狠久久av| 久久久精品免费免费高清| 亚洲欧美一区二区三区黑人 | 男人爽女人下面视频在线观看| 国产探花在线观看一区二区| 国产av在哪里看| 女人十人毛片免费观看3o分钟| 精品人妻偷拍中文字幕| 午夜激情久久久久久久| 最近的中文字幕免费完整| 亚洲精华国产精华液的使用体验| 亚洲乱码一区二区免费版| 一个人免费在线观看电影| 女人被狂操c到高潮| 国内精品一区二区在线观看| 久久综合国产亚洲精品| 婷婷六月久久综合丁香| 日本一二三区视频观看| 日韩视频在线欧美| 国产黄频视频在线观看| 午夜老司机福利剧场| 中国美白少妇内射xxxbb| 五月伊人婷婷丁香| 国产在视频线精品| 日日啪夜夜爽| 成人毛片60女人毛片免费| 日韩av免费高清视频| 最近的中文字幕免费完整| 岛国毛片在线播放| 精品午夜福利在线看| 国产乱来视频区| 人人妻人人澡人人爽人人夜夜 | 日韩制服骚丝袜av| 亚洲欧美清纯卡通| 国产激情偷乱视频一区二区| 日韩欧美精品v在线| 一区二区三区乱码不卡18| 亚洲av成人精品一区久久| 日产精品乱码卡一卡2卡三| av在线播放精品| 国产免费又黄又爽又色| 免费观看无遮挡的男女| 超碰av人人做人人爽久久| 成年av动漫网址| 国产黄色免费在线视频| 最近视频中文字幕2019在线8| 亚洲精品456在线播放app| 不卡视频在线观看欧美| 美女国产视频在线观看| 美女高潮的动态| 能在线免费观看的黄片| 在线观看人妻少妇| 少妇人妻一区二区三区视频| 91久久精品国产一区二区成人| 可以在线观看毛片的网站| 一级毛片黄色毛片免费观看视频| 丰满人妻一区二区三区视频av| 亚洲va在线va天堂va国产| 日韩成人伦理影院| 亚洲成人中文字幕在线播放| 国产亚洲av嫩草精品影院| 欧美人与善性xxx| 爱豆传媒免费全集在线观看| 久久精品久久久久久久性| 肉色欧美久久久久久久蜜桃 | 老司机影院毛片| 婷婷色综合大香蕉| 国产黄片视频在线免费观看| 最近最新中文字幕大全电影3| 久久精品人妻少妇| 六月丁香七月| 国产激情偷乱视频一区二区| 最新中文字幕久久久久| 国产精品久久久久久精品电影小说 | 内射极品少妇av片p| 免费av毛片视频| 少妇熟女欧美另类| 国产一级毛片七仙女欲春2| 一级毛片久久久久久久久女| 精品久久久噜噜| 少妇被粗大猛烈的视频| 久久久色成人| 麻豆成人av视频| 亚洲精品色激情综合| 女的被弄到高潮叫床怎么办| 日本免费在线观看一区| 久久久久久久午夜电影| 狠狠精品人妻久久久久久综合| 成人毛片60女人毛片免费| 综合色丁香网| 亚洲人成网站高清观看| 3wmmmm亚洲av在线观看| 午夜爱爱视频在线播放| 伊人久久精品亚洲午夜| 嫩草影院精品99| 免费少妇av软件| 国产一区二区在线观看日韩| 麻豆成人午夜福利视频| 欧美日韩综合久久久久久| 69人妻影院| 日韩在线高清观看一区二区三区| xxx大片免费视频| 最近最新中文字幕大全电影3| 汤姆久久久久久久影院中文字幕 | 丰满乱子伦码专区| 99热这里只有精品一区| 国产av在哪里看| 亚洲国产色片| 亚洲电影在线观看av| 欧美成人一区二区免费高清观看| 免费大片黄手机在线观看| 免费大片18禁| 中文字幕av在线有码专区| 亚洲成人精品中文字幕电影| 美女内射精品一级片tv| 少妇熟女aⅴ在线视频| 一级a做视频免费观看| 日日啪夜夜爽| 亚洲乱码一区二区免费版| 乱人视频在线观看| 尤物成人国产欧美一区二区三区| 亚洲激情五月婷婷啪啪| 成人午夜高清在线视频| 国产激情偷乱视频一区二区| 一级毛片 在线播放| 亚洲久久久久久中文字幕| 久久国产乱子免费精品| 少妇人妻精品综合一区二区| 人妻夜夜爽99麻豆av| 成人无遮挡网站| 国产白丝娇喘喷水9色精品| 校园人妻丝袜中文字幕| 国产探花极品一区二区| 午夜亚洲福利在线播放| 亚洲人成网站在线观看播放| 免费不卡的大黄色大毛片视频在线观看 | 国产精品伦人一区二区| 国产黄色免费在线视频| 成人性生交大片免费视频hd| 亚洲国产高清在线一区二区三| 91精品伊人久久大香线蕉| 国产精品99久久久久久久久| 男女国产视频网站| 美女被艹到高潮喷水动态| 久久人人爽人人爽人人片va| 国产一区有黄有色的免费视频 | 国产综合精华液| 亚洲性久久影院| 免费看美女性在线毛片视频| 91狼人影院| 超碰av人人做人人爽久久| 中文资源天堂在线| 久久久精品94久久精品| 精品国内亚洲2022精品成人| av免费在线看不卡| 午夜福利成人在线免费观看| 国产极品天堂在线| 亚洲四区av| 插阴视频在线观看视频| 国产精品.久久久| 天堂√8在线中文| 亚洲美女搞黄在线观看| 少妇被粗大猛烈的视频| 亚洲成人一二三区av| 午夜福利网站1000一区二区三区| 最近手机中文字幕大全| 尾随美女入室| 久久精品综合一区二区三区| av线在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 久久精品夜色国产| 蜜臀久久99精品久久宅男| 国产精品国产三级国产专区5o| 精品不卡国产一区二区三区| 美女cb高潮喷水在线观看| 一个人看的www免费观看视频| 久久99热这里只频精品6学生| 欧美日本视频| 国产亚洲91精品色在线| 免费少妇av软件| 免费播放大片免费观看视频在线观看| 九九爱精品视频在线观看| 亚洲av成人精品一二三区| 免费看av在线观看网站| 国产不卡一卡二| 亚洲精品456在线播放app| 中国美白少妇内射xxxbb| 日韩欧美三级三区| 中文字幕免费在线视频6| 亚洲伊人久久精品综合|