• <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)展及前緣問題
    精品久久久噜噜| 日本wwww免费看| 亚洲精品国产色婷婷电影| 亚洲婷婷狠狠爱综合网| 精品99又大又爽又粗少妇毛片| 中国美白少妇内射xxxbb| 日本-黄色视频高清免费观看| 人人妻人人澡人人看| 久久人人爽av亚洲精品天堂| 欧美丝袜亚洲另类| 国产不卡av网站在线观看| 亚洲欧美精品自产自拍| 女人精品久久久久毛片| 一区二区三区免费毛片| 在线观看人妻少妇| 免费观看av网站的网址| 国产精品秋霞免费鲁丝片| 久久人人爽av亚洲精品天堂| 美女视频免费永久观看网站| www.av在线官网国产| 中文欧美无线码| 夫妻午夜视频| 亚洲av免费高清在线观看| 高清毛片免费看| 一区二区三区四区激情视频| 国产乱来视频区| 特大巨黑吊av在线直播| 亚洲美女搞黄在线观看| 18禁观看日本| 少妇被粗大的猛进出69影院 | 亚州av有码| 夫妻午夜视频| 视频中文字幕在线观看| 大片免费播放器 马上看| 热99国产精品久久久久久7| 少妇 在线观看| 欧美 日韩 精品 国产| 国产精品麻豆人妻色哟哟久久| 国产又色又爽无遮挡免| 岛国毛片在线播放| 国产免费一级a男人的天堂| 亚洲av国产av综合av卡| 一级二级三级毛片免费看| 少妇熟女欧美另类| 精品人妻熟女毛片av久久网站| 最新的欧美精品一区二区| 一边亲一边摸免费视频| 秋霞在线观看毛片| 91精品三级在线观看| 亚洲精品,欧美精品| 国产片内射在线| 下体分泌物呈黄色| 免费看不卡的av| 91午夜精品亚洲一区二区三区| 国产精品免费大片| 国产黄色视频一区二区在线观看| 久久亚洲国产成人精品v| 久久久亚洲精品成人影院| 国产精品免费大片| 建设人人有责人人尽责人人享有的| 欧美变态另类bdsm刘玥| 国产 一区精品| 日韩成人av中文字幕在线观看| 国产精品国产三级国产专区5o| 免费看不卡的av| 2022亚洲国产成人精品| 久久99热这里只频精品6学生| 建设人人有责人人尽责人人享有的| 97超碰精品成人国产| 18在线观看网站| 激情五月婷婷亚洲| 国产 精品1| 久久婷婷青草| 国产极品天堂在线| 最近中文字幕高清免费大全6| 欧美丝袜亚洲另类| 国产成人精品在线电影| 久久久久久久久大av| 日韩欧美精品免费久久| 日韩电影二区| 制服诱惑二区| 亚洲国产最新在线播放| 一区二区日韩欧美中文字幕 | 91aial.com中文字幕在线观看| 久久ye,这里只有精品| 亚洲五月色婷婷综合| 80岁老熟妇乱子伦牲交| 男人操女人黄网站| 中文字幕av电影在线播放| 国产爽快片一区二区三区| 亚洲中文av在线| 日韩不卡一区二区三区视频在线| 欧美 亚洲 国产 日韩一| 又黄又爽又刺激的免费视频.| 欧美三级亚洲精品| √禁漫天堂资源中文www| 精品人妻熟女av久视频| 蜜桃在线观看..| 只有这里有精品99| 男男h啪啪无遮挡| av在线app专区| 国产精品久久久久久久电影| 久久久精品94久久精品| 国产一级毛片在线| 99视频精品全部免费 在线| 日韩大片免费观看网站| 97精品久久久久久久久久精品| 国产精品人妻久久久久久| 欧美亚洲 丝袜 人妻 在线| 黄片播放在线免费| 久久精品国产鲁丝片午夜精品| 亚洲成人av在线免费| 各种免费的搞黄视频| 热99国产精品久久久久久7| h视频一区二区三区| 视频在线观看一区二区三区| 91久久精品国产一区二区三区| 麻豆精品久久久久久蜜桃| 亚洲,欧美,日韩| 日韩中字成人| 我要看黄色一级片免费的| 欧美精品高潮呻吟av久久| 两个人的视频大全免费| 色视频在线一区二区三区| 久久久久久久国产电影| 久久久国产精品麻豆| 欧美最新免费一区二区三区| 成人二区视频| 久久精品国产a三级三级三级| 国产精品国产三级国产专区5o| 极品少妇高潮喷水抽搐| 两个人免费观看高清视频| 成人免费观看视频高清| 最近最新中文字幕免费大全7| 热99国产精品久久久久久7| 少妇的逼水好多| 街头女战士在线观看网站| 91精品伊人久久大香线蕉| tube8黄色片| 母亲3免费完整高清在线观看 | 99精国产麻豆久久婷婷| 免费av中文字幕在线| 韩国高清视频一区二区三区| 女人久久www免费人成看片| 久久久久久久久大av| 久久久久网色| 日韩成人伦理影院| 日日撸夜夜添| 一边摸一边做爽爽视频免费| 国产av码专区亚洲av| 免费久久久久久久精品成人欧美视频 | av黄色大香蕉| 高清不卡的av网站| 国产亚洲精品久久久com| videos熟女内射| 永久免费av网站大全| 中国国产av一级| 国产精品蜜桃在线观看| 成人亚洲精品一区在线观看| 看非洲黑人一级黄片| 免费播放大片免费观看视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 蜜臀久久99精品久久宅男| 日韩精品有码人妻一区| 日韩强制内射视频| 美女cb高潮喷水在线观看| 久久久久国产精品人妻一区二区| 新久久久久国产一级毛片| 国产在线免费精品| 国国产精品蜜臀av免费| 国产深夜福利视频在线观看| 亚洲一区二区三区欧美精品| 欧美日韩国产mv在线观看视频| 国产淫语在线视频| 亚洲熟女精品中文字幕| 精品一区二区三区视频在线| 亚洲成色77777| 一边亲一边摸免费视频| 亚洲精品456在线播放app| 美女脱内裤让男人舔精品视频| 国产亚洲午夜精品一区二区久久| 丝瓜视频免费看黄片| 国产成人精品无人区| 国产精品女同一区二区软件| 十八禁高潮呻吟视频| 一级,二级,三级黄色视频| 亚洲色图综合在线观看| 久久久精品免费免费高清| 国产 一区精品| 一本久久精品| 国产一区二区三区综合在线观看 | 国产亚洲欧美精品永久| 老熟女久久久| 国产精品一国产av| 国产精品一区www在线观看| 亚洲久久久国产精品| 国产精品国产三级国产av玫瑰| 久久精品夜色国产| 中文字幕人妻丝袜制服| 成人漫画全彩无遮挡| 国产乱人偷精品视频| 一级,二级,三级黄色视频| 国产成人免费无遮挡视频| 中文字幕久久专区| 国产精品秋霞免费鲁丝片| 狂野欧美白嫩少妇大欣赏| 色94色欧美一区二区| √禁漫天堂资源中文www| 国产爽快片一区二区三区| 国产成人精品无人区| 啦啦啦中文免费视频观看日本| 免费黄频网站在线观看国产| 国产精品无大码| av不卡在线播放| 国产毛片在线视频| 少妇的逼好多水| 亚洲av不卡在线观看| 亚洲国产精品一区二区三区在线| 丝袜美足系列| 91精品三级在线观看| 亚洲av福利一区| 欧美亚洲日本最大视频资源| 国产精品 国内视频| 亚洲成人手机| 精品少妇黑人巨大在线播放| 多毛熟女@视频| 久久精品国产a三级三级三级| 这个男人来自地球电影免费观看 | 汤姆久久久久久久影院中文字幕| 99久久精品国产国产毛片| 十八禁网站网址无遮挡| 国产一区亚洲一区在线观看| 亚洲精品乱码久久久v下载方式| 亚洲精品亚洲一区二区| 亚洲精品国产av成人精品| 伦精品一区二区三区| 亚洲综合色网址| 午夜精品国产一区二区电影| 日本与韩国留学比较| 国内精品宾馆在线| 国产高清不卡午夜福利| 久久 成人 亚洲| 久久久欧美国产精品| 少妇人妻 视频| av国产精品久久久久影院| 成人二区视频| 一个人看视频在线观看www免费| 亚洲精品亚洲一区二区| 亚洲欧洲精品一区二区精品久久久 | 汤姆久久久久久久影院中文字幕| 丰满迷人的少妇在线观看| 久久久精品区二区三区| 国产探花极品一区二区| 少妇丰满av| 能在线免费看毛片的网站| 国产亚洲av片在线观看秒播厂| 日本免费在线观看一区| 五月玫瑰六月丁香| 男人爽女人下面视频在线观看| 80岁老熟妇乱子伦牲交| 国产精品一国产av| 好男人视频免费观看在线| 久久久a久久爽久久v久久| 男人操女人黄网站| 丰满饥渴人妻一区二区三| 男女免费视频国产| 国语对白做爰xxxⅹ性视频网站| 精品少妇黑人巨大在线播放| 亚洲国产精品一区三区| 国产成人午夜福利电影在线观看| 亚洲av福利一区| 亚洲图色成人| 亚洲欧洲国产日韩| 一边摸一边做爽爽视频免费| a级毛片在线看网站| 插阴视频在线观看视频| 久久久久久久久久人人人人人人| av国产精品久久久久影院| 国产午夜精品一二区理论片| 又粗又硬又长又爽又黄的视频| 蜜桃久久精品国产亚洲av| av播播在线观看一区| 国产 精品1| 九九在线视频观看精品| 久久午夜福利片| 十八禁网站网址无遮挡| 综合色丁香网| 日韩在线高清观看一区二区三区| 欧美精品国产亚洲| 天天影视国产精品| 18禁裸乳无遮挡动漫免费视频| 欧美日韩亚洲高清精品| 国产精品无大码| 少妇人妻精品综合一区二区| 国产成人精品久久久久久| 国产免费视频播放在线视频| 国产在线免费精品| 中文字幕av电影在线播放| 日韩成人伦理影院| 2021少妇久久久久久久久久久| 国产白丝娇喘喷水9色精品| 一级二级三级毛片免费看| 99久久精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 国产欧美日韩一区二区三区在线 | 亚洲不卡免费看| 久久99蜜桃精品久久| 51国产日韩欧美| 寂寞人妻少妇视频99o| 99久久综合免费| 亚洲精品美女久久av网站| 中文字幕最新亚洲高清| 久久久欧美国产精品| 国产男女内射视频| 免费人成在线观看视频色| 18禁裸乳无遮挡动漫免费视频| 高清毛片免费看| 日本午夜av视频| 在线观看www视频免费| 伊人亚洲综合成人网| 免费日韩欧美在线观看| 免费观看av网站的网址| 2021少妇久久久久久久久久久| 国产精品一区二区在线观看99| 国产av国产精品国产| 色5月婷婷丁香| 亚洲精品国产色婷婷电影| 成年女人在线观看亚洲视频| 亚洲内射少妇av| 国产一区有黄有色的免费视频| 久久久久久伊人网av| 亚洲精品乱码久久久v下载方式| 美女主播在线视频| 如何舔出高潮| 精品久久久久久久久亚洲| 国产亚洲一区二区精品| 街头女战士在线观看网站| 99热6这里只有精品| 亚洲精品aⅴ在线观看| 在现免费观看毛片| 91久久精品电影网| 欧美xxxx性猛交bbbb| 欧美精品一区二区大全| 久久久久网色| a级毛片在线看网站| 人妻少妇偷人精品九色| 欧美精品人与动牲交sv欧美| 国产黄色视频一区二区在线观看| 日韩一本色道免费dvd| 久久久久精品性色| 男女边摸边吃奶| 国产 精品1| 天天操日日干夜夜撸| 亚洲精品日本国产第一区| 免费观看在线日韩| 国模一区二区三区四区视频| 欧美日韩一区二区视频在线观看视频在线| 国内精品宾馆在线| 中文天堂在线官网| 日韩 亚洲 欧美在线| 成人毛片60女人毛片免费| 色婷婷久久久亚洲欧美| 22中文网久久字幕| 亚洲精品第二区| 午夜日本视频在线| 三级国产精品片| 日产精品乱码卡一卡2卡三| 少妇 在线观看| 视频在线观看一区二区三区| 久久精品国产亚洲网站| 汤姆久久久久久久影院中文字幕| 另类精品久久| 日本午夜av视频| 飞空精品影院首页| 草草在线视频免费看| 老熟女久久久| 各种免费的搞黄视频| 午夜影院在线不卡| 免费久久久久久久精品成人欧美视频 | 久久久久久久大尺度免费视频| 观看av在线不卡| 久热久热在线精品观看| 视频区图区小说| 免费黄网站久久成人精品| 国产精品一国产av| 亚洲人成网站在线播| 日韩精品有码人妻一区| 丁香六月天网| 高清av免费在线| 一区二区日韩欧美中文字幕 | av卡一久久| 亚洲av二区三区四区| 18禁在线无遮挡免费观看视频| 老司机影院毛片| 丰满少妇做爰视频| 国产成人精品无人区| 国产精品一区www在线观看| 91午夜精品亚洲一区二区三区| 99热6这里只有精品| 一区二区三区精品91| 久久久精品区二区三区| 欧美bdsm另类| 伊人久久国产一区二区| 亚洲av福利一区| 午夜免费观看性视频| 丰满乱子伦码专区| 亚洲性久久影院| 校园人妻丝袜中文字幕| 女性生殖器流出的白浆| 夜夜骑夜夜射夜夜干| 三级国产精品欧美在线观看| 一级黄片播放器| 大香蕉久久成人网| 国产欧美日韩一区二区三区在线 | 精品人妻在线不人妻| 成人影院久久| 国产在线视频一区二区| 欧美日韩国产mv在线观看视频| 免费观看av网站的网址| 国产精品99久久99久久久不卡 | 国产一区二区在线观看日韩| 精品人妻一区二区三区麻豆| 少妇的逼好多水| 久久久欧美国产精品| 欧美日韩视频高清一区二区三区二| 国产成人freesex在线| 卡戴珊不雅视频在线播放| 亚洲av二区三区四区| 亚洲精品国产av成人精品| 国产免费一级a男人的天堂| a 毛片基地| 亚洲av中文av极速乱| 国产精品一二三区在线看| 久久99蜜桃精品久久| 久久人人爽人人片av| 纵有疾风起免费观看全集完整版| av在线播放精品| 在线观看一区二区三区激情| 最黄视频免费看| 九色亚洲精品在线播放| 精品国产一区二区三区久久久樱花| 观看美女的网站| 热99久久久久精品小说推荐| 国产一区二区在线观看av| 99热6这里只有精品| 在线看a的网站| 国产精品久久久久久av不卡| 全区人妻精品视频| 亚洲欧美一区二区三区国产| 国产熟女午夜一区二区三区 | 五月玫瑰六月丁香| 久久亚洲国产成人精品v| 一级毛片电影观看| 水蜜桃什么品种好| 九色亚洲精品在线播放| 色94色欧美一区二区| 黑人猛操日本美女一级片| 18禁动态无遮挡网站| 女人精品久久久久毛片| videos熟女内射| 九色成人免费人妻av| 亚洲精品久久成人aⅴ小说 | 国产精品蜜桃在线观看| 亚洲精品乱久久久久久| 国产午夜精品一二区理论片| 亚洲精品av麻豆狂野| 九九在线视频观看精品| 51国产日韩欧美| 哪个播放器可以免费观看大片| 精品人妻一区二区三区麻豆| 一级毛片电影观看| 国产精品一区二区三区四区免费观看| 插阴视频在线观看视频| 婷婷色麻豆天堂久久| 特大巨黑吊av在线直播| 亚洲精品色激情综合| 国产精品国产av在线观看| 91在线精品国自产拍蜜月| 国产熟女欧美一区二区| 亚洲欧美日韩另类电影网站| 国产日韩欧美亚洲二区| 男女高潮啪啪啪动态图| 午夜福利,免费看| 卡戴珊不雅视频在线播放| 永久网站在线| 国产视频内射| 美女国产视频在线观看| 国产精品国产av在线观看| 亚洲精品av麻豆狂野| 成年美女黄网站色视频大全免费 | 少妇丰满av| 国产日韩欧美在线精品| 国产综合精华液| 人体艺术视频欧美日本| 免费黄色在线免费观看| 欧美bdsm另类| 精品视频人人做人人爽| 一区二区三区四区激情视频| 亚洲国产精品一区二区三区在线| 狂野欧美激情性bbbbbb| 少妇丰满av| 久久狼人影院| 免费大片18禁| av.在线天堂| 99热这里只有是精品在线观看| 日韩制服骚丝袜av| 亚洲人成网站在线观看播放| 嫩草影院入口| 日韩一本色道免费dvd| 中文字幕亚洲精品专区| 自拍欧美九色日韩亚洲蝌蚪91| 欧美+日韩+精品| 性色avwww在线观看| a级毛片在线看网站| 三上悠亚av全集在线观看| 中文精品一卡2卡3卡4更新| 久久久久久久久久成人| 久久国内精品自在自线图片| 欧美日韩在线观看h| 老女人水多毛片| 一级二级三级毛片免费看| 精品国产一区二区三区久久久樱花| 少妇 在线观看| 亚洲av成人精品一二三区| 一区二区日韩欧美中文字幕 | 久久久久网色| 婷婷色综合大香蕉| 日韩电影二区| 卡戴珊不雅视频在线播放| 美女大奶头黄色视频| 亚洲精品国产av蜜桃| 有码 亚洲区| 女人久久www免费人成看片| 黄色视频在线播放观看不卡| 26uuu在线亚洲综合色| 亚洲丝袜综合中文字幕| 国产免费视频播放在线视频| 日韩av不卡免费在线播放| 亚洲精品乱码久久久v下载方式| 国国产精品蜜臀av免费| .国产精品久久| 午夜激情福利司机影院| 亚洲怡红院男人天堂| 人妻 亚洲 视频| 亚洲婷婷狠狠爱综合网| 日产精品乱码卡一卡2卡三| 国产日韩欧美亚洲二区| 久久人人爽av亚洲精品天堂| 亚洲欧美一区二区三区黑人 | 亚洲精品成人av观看孕妇| 亚洲国产av新网站| 在线播放无遮挡| 日韩欧美精品免费久久| 国产精品成人在线| 午夜福利视频在线观看免费| 黄片播放在线免费| 大话2 男鬼变身卡| 午夜精品国产一区二区电影| 一级毛片电影观看| 欧美日韩综合久久久久久| 一级黄片播放器| 成人国产av品久久久| 搡老乐熟女国产| 大陆偷拍与自拍| 日韩一区二区三区影片| 亚州av有码| 久久久久国产网址| 国产成人aa在线观看| 少妇人妻久久综合中文| 日日摸夜夜添夜夜添av毛片| 亚洲欧美清纯卡通| 欧美人与性动交α欧美精品济南到 | 男女国产视频网站| 国产熟女欧美一区二区| 999精品在线视频| 2018国产大陆天天弄谢| 日韩熟女老妇一区二区性免费视频| 亚洲精品第二区| 97超视频在线观看视频| 十八禁网站网址无遮挡| 国产精品一区www在线观看| 日韩精品有码人妻一区| 蜜桃在线观看..| 大香蕉97超碰在线| 精品国产乱码久久久久久小说| 久久毛片免费看一区二区三区| 最后的刺客免费高清国语| 国产精品成人在线| 大又大粗又爽又黄少妇毛片口| 日日摸夜夜添夜夜爱| 亚洲av中文av极速乱| 国产一区有黄有色的免费视频| 亚洲精品av麻豆狂野| 亚洲精品日韩在线中文字幕| av在线老鸭窝| 亚洲精品av麻豆狂野| 亚洲av中文av极速乱| 在线观看一区二区三区激情| 卡戴珊不雅视频在线播放| 免费高清在线观看视频在线观看| 午夜影院在线不卡| 男男h啪啪无遮挡| 日本与韩国留学比较| 久久精品熟女亚洲av麻豆精品| 大片免费播放器 马上看| 欧美精品一区二区免费开放| 男女无遮挡免费网站观看| 国产免费一区二区三区四区乱码| 91aial.com中文字幕在线观看| 在线观看免费视频网站a站| 丰满少妇做爰视频| 国产熟女午夜一区二区三区 | 热re99久久国产66热| 国产黄频视频在线观看| 美女中出高潮动态图|