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

    基于勢(shì)流理論的尾空泡對(duì)航行體表面壓力影響研究

    2021-10-11 07:23:30尤天慶權(quán)曉波劉元清鮑文春王占瑩
    船舶力學(xué) 2021年9期
    關(guān)鍵詞:物面空泡壁面

    尤天慶,權(quán)曉波,劉元清,鮑文春,王占瑩

    (1.北京宇航系統(tǒng)工程研究所,北京 100076;2.中國(guó)運(yùn)載火箭技術(shù)研究院,北京 100076)

    0 引 言

    航行體水下垂直干發(fā)射情況下,彈射氣體會(huì)附著在尾部低壓區(qū),形成附著尾空泡[1]。考慮不可壓縮流體下游的擾動(dòng)會(huì)對(duì)上游流動(dòng)產(chǎn)生影響,尾空泡狀態(tài)及演化過(guò)程必然會(huì)對(duì)上游航行體表面壓力產(chǎn)生不同程度的干擾[2],進(jìn)而對(duì)航行體整體流體動(dòng)力特性產(chǎn)生影響[3]。

    目前對(duì)尾空泡的研究主要集中于空泡形態(tài)和泡內(nèi)壓力的預(yù)示,其計(jì)算分析方法主要分為簡(jiǎn)化工程計(jì)算[4]以及基于雷諾平均N-S方程的流場(chǎng)數(shù)值模擬兩類[5]。針對(duì)流線型航行體及尾空泡組合體,在其流動(dòng)分離不明顯的情況下,基于勢(shì)流的邊界元計(jì)算可在有限計(jì)算資源下分析流場(chǎng)特性,從計(jì)算效率和流場(chǎng)信息獲取等方面比較,是一種介于工程簡(jiǎn)化計(jì)算和基于雷諾平均N-S方程仿真計(jì)算的方法。

    運(yùn)用勢(shì)流理論對(duì)空泡發(fā)展演化的研究,一方面集中于水中孤立氣泡的生長(zhǎng)和潰滅[6],最早見(jiàn)于Blake和Gibson的工作[7],另一方面集中于航行體肩部附著空泡研究[8-9]。相關(guān)研究表明基于勢(shì)流的邊界元方法,在空泡形態(tài)、流場(chǎng)結(jié)構(gòu)以及物面整體受力等方面的研究分析工作中發(fā)揮了重要作用。

    本文基于勢(shì)流理論,利用邊界元數(shù)值方法,建立了二維軸對(duì)稱航行體及附著尾空泡的流場(chǎng)計(jì)算方法,并形成了計(jì)算程序。在此基礎(chǔ)上計(jì)算分析了不同空泡數(shù)條件下尾空泡對(duì)上游航行體物面壓力分布的影響,并以此為基礎(chǔ)進(jìn)一步分析了瞬態(tài)情況下尾空泡收縮過(guò)程對(duì)上游壓力的擾動(dòng),基于數(shù)值計(jì)算結(jié)果結(jié)合實(shí)驗(yàn)數(shù)據(jù)說(shuō)明了壓力擾動(dòng)的衰減規(guī)律。

    1 數(shù)值計(jì)算方法

    1.1 控制方程

    假設(shè)空泡附著在航行體后部擾流低壓區(qū),在無(wú)限大流體域內(nèi)發(fā)展演化。并假設(shè)流體無(wú)粘、無(wú)旋且不可壓,則N-S方程可簡(jiǎn)化為拉普拉斯方程。

    式中,V為速度向量,U∞為來(lái)流速度(設(shè)來(lái)流為x軸正向),Φ為總速度勢(shì),x為位置坐標(biāo),φ為擾動(dòng)速度勢(shì)。

    在物面,流體需滿足不可穿透條件,即壁面法向速度為零:

    式中,n為物面法向方向向量。

    在空泡壁面,流體壓力在流線方向上為恒定,在計(jì)算中主要通過(guò)伯努利方程計(jì)算空泡切向速度,進(jìn)而通過(guò)延空泡壁面流線方向,進(jìn)行積分獲得空泡壁面上的速度勢(shì)分布。t為物面切向方向向量,lc為空泡壁面流線方向長(zhǎng)度。

    對(duì)于空泡周圍流場(chǎng)的瞬態(tài)變化,速度勢(shì)Φ、速度V及壓力p需滿足伯努利方程:

    1.2 數(shù)值算法

    設(shè)空間場(chǎng)點(diǎn)為p,物面和空泡壁面上分布的點(diǎn)源位置為q,其強(qiáng)度為λ(q),則其在p點(diǎn)產(chǎn)生的速度勢(shì)和速度分別為

    式中,G(p,q)為格林函數(shù)。對(duì)于拉普拉斯方程,三維情況下格林函數(shù)為G(p,q)=-1/4πr,其中r=|p-q|。

    從拉普拉斯方程邊值問(wèn)題角度,本文為第三類邊界條件,即混合邊界條件。航行體表面為流動(dòng)不可穿透條件,即第二類邊界條件(Newmann條件)。數(shù)值計(jì)算中,航行體與空泡表面被劃分為N個(gè)單元,對(duì)航行體表面的第i個(gè)單元,考慮在邊界上源匯自身對(duì)速度的影響,擾動(dòng)勢(shì)所滿足的邊界元積分控制方程如式(8)所示,式中Sj為單元面積,ni為單元法向向量。

    空泡壁面為第一類邊界條件(Dirichlet條件),考慮沿航行體子午線方向,尾空泡分離點(diǎn)位于第d個(gè)單元,其擾動(dòng)勢(shì)所滿足的邊界元積分控制方程如式(9)所示,式中ti為單元切向向量。

    式(8)和(9)即為邊界元計(jì)算的基本方程,通過(guò)這兩個(gè)方程即可求得物面及空泡壁面的λ分布,進(jìn)而確定流場(chǎng)中速度勢(shì)、速度及壓力。

    為節(jié)省計(jì)算資源提高計(jì)算效率,本文采用二維軸對(duì)稱計(jì)算模型,在航行體及空泡表面的邊界單元為圓環(huán)狀。與速度勢(shì)及速度相關(guān)的源匯均勻沿圓環(huán)分布,如圖1所示。對(duì)于子午線長(zhǎng)度為l圓環(huán)邊界單元,其影響速度勢(shì)和速度由圓環(huán)上源匯積分獲得,其表達(dá)式如下:

    圖1 軸對(duì)稱單元積分示意圖Fig.1 Configuration of axisymmetric element integration

    由于第一類完整橢圓積分函數(shù)的奇異性,當(dāng)R趨近于1時(shí),式(9)~(11)所示的邊界單元速度勢(shì)和速度影響函數(shù)值趨近于無(wú)窮大,但其積分是有限的。本文利用偶數(shù)積分點(diǎn)的Gauss積分進(jìn)行邊界單元積分計(jì)算。

    為驗(yàn)證數(shù)值計(jì)算代碼的正確性,計(jì)算一無(wú)粘無(wú)旋理想流體條件下軸對(duì)稱旋成體零攻角物面的壓力分布,并與水洞測(cè)壓數(shù)據(jù)對(duì)比。對(duì)比表明數(shù)值計(jì)算與數(shù)據(jù)吻合較好,如圖2所示。

    圖2 數(shù)值計(jì)算與試驗(yàn)數(shù)據(jù)對(duì)比Fig.2 Comparison between calculation and experiment results

    1.3 尾空泡閉合模式

    針對(duì)物面階梯后部低壓區(qū)附著空泡尾部呈現(xiàn)出很明顯的湍流,如果利用勢(shì)流理論對(duì)該問(wèn)題進(jìn)行計(jì)算,就必須將這一區(qū)域進(jìn)行理想簡(jiǎn)化處理。目前針對(duì)不同類型的空泡,已有多種空泡閉合模式[10]。本文根據(jù)附著空泡的形態(tài)特征,考慮數(shù)值計(jì)算的穩(wěn)定性,采用了圓球壁面條件對(duì)空泡尾部流場(chǎng)進(jìn)行了簡(jiǎn)化處理,如圖3所示。空泡壁面邊界條件利用式(8)計(jì)算,圓球壁面利用式(9)計(jì)算。

    2 數(shù)值計(jì)算結(jié)果分析

    2.1 穩(wěn)態(tài)計(jì)算

    采用上述數(shù)值計(jì)算方法,對(duì)不同來(lái)流條件下軸對(duì)稱半球頭型航行體后部附著尾空泡的二維軸對(duì)稱流場(chǎng)進(jìn)行了迭代計(jì)算,獲得收斂后的空泡形態(tài)如圖4所示。σ=(p∞-pc)/0.5ρU2∞為空化數(shù),p∞為環(huán)境壓力,pc為尾空泡內(nèi)部壓力,ρ為流體密度,U∞為來(lái)流速度大小,X/L和Y/L分別為無(wú)量綱的軸向和徑向坐標(biāo),L為航行體長(zhǎng)度。X/L≤0部分為航行體物面,X/L>0部分為不同空化數(shù)下的附著尾空泡形態(tài)。計(jì)算結(jié)果表明穩(wěn)態(tài)空泡形態(tài)與空化數(shù)緊密對(duì)應(yīng),即針對(duì)不同的環(huán)境壓力、泡內(nèi)壓力或來(lái)流速度,只要空泡數(shù)一致,勢(shì)流條件下穩(wěn)態(tài)空泡形態(tài)即確定。不同空泡數(shù)計(jì)算結(jié)果表明,隨著空泡數(shù)的降低,尾空泡長(zhǎng)度增加,尾空泡輪廓曲線曲率變化趨緩。

    圖4 穩(wěn)態(tài)計(jì)算所形成的附著空泡Fig.4 Tail bubble profile of steady calculation

    不同空化數(shù)條件下航行體和空泡表面的壓力系數(shù)分布如圖5所示。與穩(wěn)態(tài)空泡形態(tài)和空泡數(shù)的對(duì)應(yīng)關(guān)系相同,穩(wěn)態(tài)計(jì)算的壓力系數(shù)分布也與空泡數(shù)緊密相關(guān)。不同空泡數(shù)的數(shù)值計(jì)算表明,隨著空泡數(shù)降低,尾空泡尺寸增大,壓力系數(shù)平直段也相應(yīng)地延長(zhǎng)。不同空化數(shù)條件下,遠(yuǎn)離尾空泡的航行體沾濕區(qū)壓力保持一致。尾空泡主要對(duì)上游臨近尾空泡的物面壓力產(chǎn)生較大影響,如圖5中放大圖所示。數(shù)值計(jì)算表明,隨著空化數(shù)提高,泡內(nèi)壓力系數(shù)逐漸降低,導(dǎo)致臨近尾空泡的物面壓力越低,其影響范圍也相應(yīng)擴(kuò)大。

    圖5 沿航行體及空泡表面壓力分布Fig.5 Pressure distribution along vehicle and bubble meridian

    2.2 瞬態(tài)計(jì)算

    為分析空泡瞬態(tài)變化對(duì)物面壓力的影響,采用與穩(wěn)態(tài)計(jì)算相同的航行體外形,在穩(wěn)態(tài)計(jì)算結(jié)果上,降低空泡內(nèi)部壓力,計(jì)算空泡收縮過(guò)程流場(chǎng)變化??张菔湛s過(guò)程形態(tài)變化如圖6所示。由于空泡內(nèi)部壓力降低,空泡收縮長(zhǎng)度縮短,空泡輪廓曲線變化趨于劇烈??张葑兓鸬牧鲌?chǎng)速度及加速度改變,進(jìn)一步導(dǎo)致了上游物面壓力的變化。下面主要針對(duì)空泡收縮過(guò)程的t0+dt和t0+4dt兩時(shí)刻流場(chǎng)參數(shù)變化進(jìn)行分析。

    圖6 空泡收縮過(guò)程形態(tài)變化Fig.6 Cavity profile change while shrinking

    由于空泡收縮過(guò)程中,整個(gè)航行體物面和空泡組合體形態(tài)發(fā)生變化,其長(zhǎng)度更短,流場(chǎng)下游空泡輪廓線曲率變化更加劇烈,因此在接近尾空泡的航行體表面流體流速增大。由此動(dòng)壓系數(shù)1-q/q∞在X/L=0附近明顯降低,如圖7所示。其中,q為流體當(dāng)?shù)靥巹?dòng)壓,q∞為無(wú)窮遠(yuǎn)處動(dòng)壓。

    圖7 收縮過(guò)程不同時(shí)刻物面及空泡表面動(dòng)壓分布對(duì)比Fig.7 Comparison of dynamic pressure distribution

    針對(duì)尾空泡收縮瞬態(tài)變化過(guò)程,本文流場(chǎng)壓力計(jì)算采用非穩(wěn)態(tài)伯努利方程,如式(5)所示。基于此,壓力系數(shù)為Cp=1-q/q∞-ρΦ˙/q∞。式中,ρ為流體密度,Φ˙=?Φ/?t為速度勢(shì)變化率??梢?jiàn)引起壓力變化的因素為兩部分:1-q/q∞為動(dòng)壓變化所導(dǎo)致的壓力系數(shù)變化,其分布如圖7所示;-ρΦ˙/q∞為速度勢(shì)變化率導(dǎo)致的流場(chǎng)壓力變化。兩者綜合作用表現(xiàn)為物面壓力系數(shù)的分布,如圖8所示。

    對(duì)比圖7和圖8,在靠近尾空泡的航行體物面區(qū)域X/L=0附近,動(dòng)壓導(dǎo)致的壓力明顯降低,而實(shí)際此段壓力在空泡收縮過(guò)程中壓力明顯升高,兩者變化趨勢(shì)截然相反。導(dǎo)致這種差異的原因是速度勢(shì)變化率導(dǎo)致的流場(chǎng)壓力變化,即壓力計(jì)算中的-ρΦ˙/q∞,其分布如圖9所示。圖9表明,在空泡收縮過(guò)程中,其所引起的流場(chǎng)中速度勢(shì)變化率對(duì)空泡上游物面壓力有著較大影響,其影響超過(guò)流場(chǎng)中速度變化的影響,導(dǎo)致靠近空泡部分物面壓力的上升。速度勢(shì)變化率Φ˙=?Φ/?t為歐拉方程中流體加速度項(xiàng)沿流線方向的積分。計(jì)算中速度勢(shì)變化率導(dǎo)致的流場(chǎng)壓力變化,本質(zhì)上是由流體加速度所引起。

    圖8 收縮過(guò)程壓力分布對(duì)比Fig.8 Comparison of static pressure distribution

    圖9 收縮過(guò)程速度勢(shì)變化率引起的壓力分布對(duì)比Fig.9 Time derivative of the velocity potential

    對(duì)于圖10顯示的空泡形態(tài)變化,初始時(shí)刻泡內(nèi)壓力越低,尾空泡附近流體質(zhì)點(diǎn)運(yùn)動(dòng)加速度就越大。在相同時(shí)間內(nèi)速度變化量也越大,尾空泡收縮程度越劇烈。

    在圖10顯示的壓力變化方面,由于尾空泡內(nèi)壓力在空間上是均勻的,在其收縮過(guò)程中空泡壁面壓力整體升高,即X/L>0部分壓力變化ΔCp在空間上也均勻分布。尾空泡收縮過(guò)程引起流場(chǎng)上游航行體壁面壓力變化,這種對(duì)流場(chǎng)的擾動(dòng)在接近尾空泡處最大,在遠(yuǎn)離尾空泡過(guò)程中逐漸衰減,衰減趨勢(shì)近似呈雙指數(shù)函數(shù)。

    圖10 不同初始?jí)毫ξ部张菔湛s引起的物面壓力變化Fig.10 Pressure change during the tail cavity shrinking process with different initialpressures

    為驗(yàn)證以上壓力擾動(dòng)空間分布規(guī)律的正確性,提取了航行體水下附著尾空泡收縮過(guò)程上游物面壓力測(cè)點(diǎn)變化情況,如圖10中的實(shí)驗(yàn)數(shù)據(jù)點(diǎn)。對(duì)比結(jié)果表明,數(shù)值計(jì)算壓力空間衰減趨勢(shì)與試驗(yàn)吻合較好。

    不同尾空泡收縮初始泡內(nèi)壓力會(huì)導(dǎo)致流體運(yùn)動(dòng)加速度產(chǎn)生明顯差異,進(jìn)而影響尾空泡瞬態(tài)收縮過(guò)程中的物面壓力,如圖10所示。尾空泡收縮初始泡內(nèi)壓力越小,流場(chǎng)中流體運(yùn)動(dòng)加速度就越大,對(duì)上游航行體物面壓力的擾動(dòng)也越大。

    3 結(jié) 論

    針對(duì)水下垂直干發(fā)射尾空泡對(duì)上游物面壓力影響問(wèn)題,本文基于勢(shì)流理論,考慮航行體物面及空泡壁面混合邊界條件,建立了適用于尾空泡流場(chǎng)發(fā)展演化的二維軸對(duì)稱邊界元數(shù)值計(jì)算方法,對(duì)航行體及尾空泡流場(chǎng)進(jìn)行了計(jì)算分析,得到的結(jié)論如下:

    (1)不同狀態(tài)尾空泡流場(chǎng)穩(wěn)態(tài)數(shù)值計(jì)算表明,隨著空化數(shù)的提高,尾空泡壓力系數(shù)逐漸降低,導(dǎo)致臨近尾空泡的物面壓力系數(shù)下降,尾空泡對(duì)上游物面壓力影響范圍也相應(yīng)擴(kuò)大;

    (2)尾空泡收縮瞬態(tài)流場(chǎng)演化計(jì)算表明,收縮引起的流場(chǎng)加速度對(duì)空泡上游物面壓力有著較大影響,導(dǎo)致流場(chǎng)中靠近尾空泡物面壓力的上升;

    (3)不同尾空泡收縮過(guò)程初始?jí)毫τ?jì)算表明,尾空泡收縮初始泡內(nèi)壓力越小,流場(chǎng)中流體運(yùn)動(dòng)加速度就越大,對(duì)上游航行體物面壓力的擾動(dòng)也越大,這種壓力擾動(dòng)在遠(yuǎn)離尾空泡過(guò)程中呈雙指數(shù)函數(shù)逐漸衰減。

    猜你喜歡
    物面空泡壁面
    二維有限長(zhǎng)度柔性壁面上T-S波演化的數(shù)值研究
    激波/湍流邊界層干擾壓力脈動(dòng)特性數(shù)值研究1)
    水下航行體雙空泡相互作用數(shù)值模擬研究
    讓吸盤掛鉤更牢固
    壁面溫度對(duì)微型內(nèi)燃機(jī)燃燒特性的影響
    基于LPV的超空泡航行體H∞抗飽和控制
    基于CFD的對(duì)轉(zhuǎn)槳無(wú)空泡噪聲的仿真預(yù)報(bào)
    船海工程(2015年4期)2016-01-05 15:53:28
    新型單面陣自由曲面光學(xué)測(cè)量方法成像特性仿真
    彎曲網(wǎng)格上的間斷有限元湍流數(shù)值解法研究
    SPH在水下高速物體空泡發(fā)展模擬中的應(yīng)用
    中文字幕另类日韩欧美亚洲嫩草| 一级a爱片免费观看的视频| 国产精品久久久久久亚洲av鲁大| 国产男靠女视频免费网站| 亚洲精品国产一区二区精华液| 又大又爽又粗| 国产免费av片在线观看野外av| 一本综合久久免费| 欧美日韩瑟瑟在线播放| 久久精品亚洲精品国产色婷小说| 免费在线观看影片大全网站| 国产av一区二区精品久久| 久久中文看片网| 丁香欧美五月| 久久久久久久午夜电影| 色精品久久人妻99蜜桃| 非洲黑人性xxxx精品又粗又长| 美女 人体艺术 gogo| 欧美另类亚洲清纯唯美| 可以在线观看毛片的网站| 免费在线观看视频国产中文字幕亚洲| 在线观看舔阴道视频| 50天的宝宝边吃奶边哭怎么回事| 日韩精品中文字幕看吧| 麻豆久久精品国产亚洲av| 国产av又大| 欧美丝袜亚洲另类 | 欧美黄色淫秽网站| 日韩一卡2卡3卡4卡2021年| 琪琪午夜伦伦电影理论片6080| 欧美乱码精品一区二区三区| 精品卡一卡二卡四卡免费| 精品第一国产精品| 亚洲欧美精品综合一区二区三区| 欧美+亚洲+日韩+国产| 亚洲中文字幕一区二区三区有码在线看 | 色精品久久人妻99蜜桃| 国产成年人精品一区二区| 免费少妇av软件| 村上凉子中文字幕在线| 超碰成人久久| 国产精品永久免费网站| 欧洲精品卡2卡3卡4卡5卡区| 国产国语露脸激情在线看| 国产精品久久久久久亚洲av鲁大| 亚洲国产欧美日韩在线播放| 久热爱精品视频在线9| 97人妻天天添夜夜摸| 在线十欧美十亚洲十日本专区| 欧美乱妇无乱码| 久久九九热精品免费| 嫁个100分男人电影在线观看| 99久久精品国产亚洲精品| 亚洲国产中文字幕在线视频| 中文字幕精品免费在线观看视频| 好看av亚洲va欧美ⅴa在| 手机成人av网站| 亚洲熟女毛片儿| 制服诱惑二区| av有码第一页| 久久精品人人爽人人爽视色| 亚洲av第一区精品v没综合| 18禁黄网站禁片午夜丰满| 色尼玛亚洲综合影院| 免费女性裸体啪啪无遮挡网站| 伊人久久大香线蕉亚洲五| 亚洲情色 制服丝袜| 97超级碰碰碰精品色视频在线观看| 丝袜人妻中文字幕| 好男人在线观看高清免费视频 | 国产成人精品久久二区二区91| 制服丝袜大香蕉在线| 在线天堂中文资源库| 欧美日韩瑟瑟在线播放| 亚洲欧美日韩另类电影网站| 少妇 在线观看| 亚洲精华国产精华精| 欧美成人性av电影在线观看| 51午夜福利影视在线观看| 国产高清有码在线观看视频 | 美女 人体艺术 gogo| 精品国产一区二区久久| 大香蕉久久成人网| 好男人在线观看高清免费视频 | 一级作爱视频免费观看| 九色亚洲精品在线播放| 极品人妻少妇av视频| 午夜福利在线观看吧| 亚洲七黄色美女视频| 高清黄色对白视频在线免费看| 亚洲精品美女久久久久99蜜臀| 在线播放国产精品三级| 精品国产一区二区三区四区第35| 日韩国内少妇激情av| 精品国产国语对白av| 精品国产国语对白av| 国产1区2区3区精品| 一级毛片女人18水好多| 国产单亲对白刺激| 亚洲五月天丁香| 免费搜索国产男女视频| 亚洲第一欧美日韩一区二区三区| 精品久久久久久久久久免费视频| 国产日韩一区二区三区精品不卡| 日韩精品中文字幕看吧| 日本黄色视频三级网站网址| 欧美老熟妇乱子伦牲交| 亚洲五月婷婷丁香| www.熟女人妻精品国产| 777久久人妻少妇嫩草av网站| 咕卡用的链子| 少妇粗大呻吟视频| 777久久人妻少妇嫩草av网站| 亚洲avbb在线观看| 亚洲色图综合在线观看| 亚洲精品国产区一区二| 乱人伦中国视频| 亚洲精品美女久久av网站| 一级,二级,三级黄色视频| 日韩大码丰满熟妇| 欧美日韩福利视频一区二区| 一级毛片女人18水好多| 欧美日韩黄片免| 亚洲中文字幕一区二区三区有码在线看 | 久久亚洲精品不卡| 亚洲国产欧美日韩在线播放| www.999成人在线观看| 亚洲aⅴ乱码一区二区在线播放 | 人成视频在线观看免费观看| 亚洲欧美精品综合久久99| 黄色毛片三级朝国网站| 在线视频色国产色| 搡老妇女老女人老熟妇| 亚洲人成伊人成综合网2020| 国产精品野战在线观看| 老熟妇仑乱视频hdxx| 久久婷婷人人爽人人干人人爱 | 亚洲中文字幕一区二区三区有码在线看 | 成人av一区二区三区在线看| 国产精品香港三级国产av潘金莲| 亚洲七黄色美女视频| 亚洲人成伊人成综合网2020| 国产欧美日韩一区二区精品| 国产欧美日韩综合在线一区二区| 午夜免费鲁丝| 欧美 亚洲 国产 日韩一| 韩国精品一区二区三区| 成在线人永久免费视频| 88av欧美| 日本 欧美在线| 女警被强在线播放| 亚洲欧美日韩高清在线视频| 黄色a级毛片大全视频| 欧美乱色亚洲激情| 国产男靠女视频免费网站| 国产欧美日韩精品亚洲av| 久久亚洲真实| 精品久久蜜臀av无| 欧美性长视频在线观看| 91老司机精品| 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| 在线观看午夜福利视频| 国产熟女xx| 久久亚洲精品不卡| 这个男人来自地球电影免费观看| 亚洲av成人一区二区三| 村上凉子中文字幕在线| 亚洲av第一区精品v没综合| 国产高清videossex| 女性被躁到高潮视频| 精品人妻1区二区| 亚洲熟妇中文字幕五十中出| 国产高清有码在线观看视频 | 国产精品一区二区精品视频观看| 国产亚洲精品综合一区在线观看 | 亚洲av电影在线进入| 高清毛片免费观看视频网站| 成人18禁高潮啪啪吃奶动态图| 精品国产国语对白av| 国产单亲对白刺激| 热99re8久久精品国产| 免费在线观看日本一区| 精品久久蜜臀av无| 免费久久久久久久精品成人欧美视频| 亚洲中文日韩欧美视频| 伊人久久大香线蕉亚洲五| 亚洲成av片中文字幕在线观看| 欧美日韩亚洲综合一区二区三区_| 国产蜜桃级精品一区二区三区| 亚洲av日韩精品久久久久久密| 久久精品人人爽人人爽视色| 制服诱惑二区| 免费久久久久久久精品成人欧美视频| 国产精品野战在线观看| 涩涩av久久男人的天堂| 欧美激情高清一区二区三区| 亚洲avbb在线观看| 欧美成人免费av一区二区三区| 久久精品亚洲精品国产色婷小说| 午夜福利欧美成人| 午夜亚洲福利在线播放| 制服诱惑二区| 久久 成人 亚洲| 男人舔女人的私密视频| 亚洲欧美一区二区三区黑人| 久久午夜综合久久蜜桃| 久久久久国产一级毛片高清牌| 不卡av一区二区三区| 亚洲国产看品久久| 国产亚洲欧美在线一区二区| 最新在线观看一区二区三区| 欧美国产日韩亚洲一区| 亚洲欧美精品综合一区二区三区| 午夜福利18| 中文字幕高清在线视频| 亚洲国产高清在线一区二区三 | 国产欧美日韩一区二区精品| 中文字幕av电影在线播放| 日日夜夜操网爽| 日韩精品青青久久久久久| 亚洲av电影在线进入| 国产伦人伦偷精品视频| 久久亚洲真实| bbb黄色大片| 脱女人内裤的视频| 亚洲专区国产一区二区| 欧美成人午夜精品| 久久精品国产清高在天天线| 婷婷丁香在线五月| 又黄又爽又免费观看的视频| 亚洲久久久国产精品| 欧美人与性动交α欧美精品济南到| 老司机午夜十八禁免费视频| 丁香六月欧美| 麻豆国产av国片精品| 久久久久国内视频| 国产成人av激情在线播放| 国产激情久久老熟女| 99香蕉大伊视频| 亚洲五月色婷婷综合| 757午夜福利合集在线观看| 亚洲少妇的诱惑av| 老汉色av国产亚洲站长工具| 亚洲五月婷婷丁香| 精品国内亚洲2022精品成人| 天堂影院成人在线观看| 999精品在线视频| 欧美黄色片欧美黄色片| 久99久视频精品免费| 99久久国产精品久久久| 色播在线永久视频| 一a级毛片在线观看| 免费在线观看影片大全网站| or卡值多少钱| 黄色片一级片一级黄色片| 高潮久久久久久久久久久不卡| 亚洲色图 男人天堂 中文字幕| 九色国产91popny在线| 淫妇啪啪啪对白视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品永久免费网站| 欧美一区二区精品小视频在线| 黑人巨大精品欧美一区二区mp4| 午夜福利一区二区在线看| 18美女黄网站色大片免费观看| 天堂影院成人在线观看| 久久久久久人人人人人| 欧美激情高清一区二区三区| bbb黄色大片| 亚洲一区二区三区不卡视频| av视频免费观看在线观看| 黄频高清免费视频| tocl精华| 看片在线看免费视频| 国产片内射在线| 狂野欧美激情性xxxx| 少妇熟女aⅴ在线视频| 欧美乱码精品一区二区三区| 午夜视频精品福利| 久久中文字幕一级| 一进一出抽搐动态| 777久久人妻少妇嫩草av网站| 91字幕亚洲| 国产av精品麻豆| 久久这里只有精品19| 国产精品永久免费网站| 日韩中文字幕欧美一区二区| or卡值多少钱| x7x7x7水蜜桃| 9热在线视频观看99| 国产一卡二卡三卡精品| 日韩有码中文字幕| 欧美乱色亚洲激情| 黄色女人牲交| 19禁男女啪啪无遮挡网站| 最新在线观看一区二区三区| 久久国产精品影院| 亚洲欧美激情在线| 男人舔女人的私密视频| 国产成人影院久久av| 精品免费久久久久久久清纯| 极品人妻少妇av视频| 两性午夜刺激爽爽歪歪视频在线观看 | 涩涩av久久男人的天堂| 夜夜爽天天搞| 国产精品电影一区二区三区| 脱女人内裤的视频| 国产男靠女视频免费网站| 麻豆成人av在线观看| 97人妻天天添夜夜摸| 咕卡用的链子| 9热在线视频观看99| 亚洲片人在线观看| 悠悠久久av| 日本精品一区二区三区蜜桃| 欧美色视频一区免费| 99精品在免费线老司机午夜| 国产成人系列免费观看| 国产精品久久视频播放| 亚洲av片天天在线观看| 久久婷婷成人综合色麻豆| 亚洲 国产 在线| 久久中文看片网| 国产精品秋霞免费鲁丝片| 欧美国产精品va在线观看不卡| 精品人妻在线不人妻| 18禁国产床啪视频网站| 久久久国产欧美日韩av| 日本a在线网址| 欧美另类亚洲清纯唯美| 亚洲中文字幕日韩| 亚洲自偷自拍图片 自拍| 亚洲全国av大片| 一个人观看的视频www高清免费观看 | 亚洲av第一区精品v没综合| 欧美成人免费av一区二区三区| 国产成人欧美| 天堂√8在线中文| 三级毛片av免费| 身体一侧抽搐| 日日干狠狠操夜夜爽| 久久久水蜜桃国产精品网| 大型av网站在线播放| 亚洲自偷自拍图片 自拍| 黄色视频,在线免费观看| 久久天躁狠狠躁夜夜2o2o| 国产1区2区3区精品| 欧美黑人欧美精品刺激| 色综合欧美亚洲国产小说| 老司机在亚洲福利影院| 免费搜索国产男女视频| 国产激情欧美一区二区| 丁香欧美五月| 真人一进一出gif抽搐免费| 国产1区2区3区精品| 一区福利在线观看| 亚洲人成伊人成综合网2020| 国产亚洲精品久久久久久毛片| 成人免费观看视频高清| 亚洲国产精品久久男人天堂| 成人国产一区最新在线观看| 麻豆av在线久日| 亚洲色图av天堂| 亚洲七黄色美女视频| 久久精品亚洲熟妇少妇任你| 99国产精品一区二区蜜桃av| av免费在线观看网站| 一级片免费观看大全| 久久久久国产精品人妻aⅴ院| 亚洲欧美一区二区三区黑人| 免费无遮挡裸体视频| 亚洲黑人精品在线| 黄色成人免费大全| 在线观看免费日韩欧美大片| 每晚都被弄得嗷嗷叫到高潮| 在线观看日韩欧美| 天天躁夜夜躁狠狠躁躁| 免费不卡黄色视频| 国产精品久久久av美女十八| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产毛片av蜜桃av| 国产激情久久老熟女| 一区二区日韩欧美中文字幕| 国产日韩一区二区三区精品不卡| 亚洲av成人av| 国产精品野战在线观看| av在线天堂中文字幕| 精品卡一卡二卡四卡免费| 国产成人精品在线电影| 免费在线观看视频国产中文字幕亚洲| 亚洲国产精品成人综合色| 国产精品99久久99久久久不卡| 亚洲五月天丁香| 黑人巨大精品欧美一区二区mp4| 国产成人影院久久av| 少妇熟女aⅴ在线视频| 女生性感内裤真人,穿戴方法视频| 午夜久久久在线观看| 亚洲国产精品合色在线| 国产欧美日韩精品亚洲av| 亚洲欧美精品综合久久99| 每晚都被弄得嗷嗷叫到高潮| 亚洲黑人精品在线| 热re99久久国产66热| 成人av一区二区三区在线看| 国产一区二区激情短视频| 婷婷丁香在线五月| 精品人妻1区二区| 欧美成狂野欧美在线观看| 女性被躁到高潮视频| 19禁男女啪啪无遮挡网站| 亚洲国产精品久久男人天堂| 日韩欧美国产一区二区入口| 十八禁人妻一区二区| 在线国产一区二区在线| 亚洲成av片中文字幕在线观看| 午夜精品久久久久久毛片777| 少妇的丰满在线观看| 日韩免费av在线播放| 国产精品久久久久久人妻精品电影| 久久青草综合色| 久久婷婷人人爽人人干人人爱 | 免费观看精品视频网站| 久久国产乱子伦精品免费另类| 精品免费久久久久久久清纯| 黄色 视频免费看| 国内毛片毛片毛片毛片毛片| 亚洲情色 制服丝袜| 亚洲精品在线观看二区| 亚洲免费av在线视频| 身体一侧抽搐| 精品国产一区二区久久| 成年版毛片免费区| 男女之事视频高清在线观看| 国产aⅴ精品一区二区三区波| 欧美午夜高清在线| 一进一出抽搐动态| 妹子高潮喷水视频| 香蕉丝袜av| 欧美绝顶高潮抽搐喷水| 91麻豆精品激情在线观看国产| 多毛熟女@视频| 日本五十路高清| 亚洲五月婷婷丁香| 大型黄色视频在线免费观看| 久久香蕉激情| 在线观看www视频免费| 长腿黑丝高跟| 国产成+人综合+亚洲专区| 人成视频在线观看免费观看| 动漫黄色视频在线观看| 久久久久久久久久久久大奶| 亚洲三区欧美一区| 母亲3免费完整高清在线观看| 在线观看舔阴道视频| 人人妻人人爽人人添夜夜欢视频| 国产亚洲欧美精品永久| 亚洲视频免费观看视频| 久久久久国产一级毛片高清牌| 国产精品久久久久久亚洲av鲁大| 如日韩欧美国产精品一区二区三区| 亚洲国产欧美日韩在线播放| 欧美成狂野欧美在线观看| 精品不卡国产一区二区三区| 国产成人精品久久二区二区91| 亚洲精品一卡2卡三卡4卡5卡| 亚洲aⅴ乱码一区二区在线播放 | 黄色视频不卡| 亚洲男人的天堂狠狠| 欧美+亚洲+日韩+国产| 国产熟女午夜一区二区三区| 亚洲国产高清在线一区二区三 | 亚洲中文日韩欧美视频| 亚洲激情在线av| 色av中文字幕| 一级a爱片免费观看的视频| 99精品久久久久人妻精品| 老熟妇乱子伦视频在线观看| 18禁美女被吸乳视频| 香蕉久久夜色| 精品久久久久久久毛片微露脸| 99国产精品99久久久久| 美女高潮到喷水免费观看| 丰满人妻熟妇乱又伦精品不卡| 18禁观看日本| 久久午夜亚洲精品久久| 一边摸一边做爽爽视频免费| 啪啪无遮挡十八禁网站| 色综合婷婷激情| 日韩av在线大香蕉| 亚洲最大成人中文| 叶爱在线成人免费视频播放| 一级a爱视频在线免费观看| 色精品久久人妻99蜜桃| 欧美性长视频在线观看| 亚洲成人国产一区在线观看| 国产成人免费无遮挡视频| 99精品欧美一区二区三区四区| 中文字幕最新亚洲高清| 欧美在线黄色| 国产麻豆69| 自线自在国产av| 免费在线观看黄色视频的| 亚洲国产精品合色在线| 真人一进一出gif抽搐免费| 精品一区二区三区视频在线观看免费| 在线观看免费日韩欧美大片| 曰老女人黄片| tocl精华| 免费在线观看影片大全网站| 亚洲一卡2卡3卡4卡5卡精品中文| 淫妇啪啪啪对白视频| 亚洲第一青青草原| 脱女人内裤的视频| 国产一区二区三区综合在线观看| 亚洲男人的天堂狠狠| 19禁男女啪啪无遮挡网站| 精品国内亚洲2022精品成人| АⅤ资源中文在线天堂| 精品卡一卡二卡四卡免费| 亚洲av电影在线进入| 欧美久久黑人一区二区| 99re在线观看精品视频| 国产亚洲欧美精品永久| 久久亚洲精品不卡| 国产欧美日韩综合在线一区二区| 男人舔女人下体高潮全视频| www国产在线视频色| 最好的美女福利视频网| 丝袜人妻中文字幕| 岛国在线观看网站| 亚洲久久久国产精品| 桃色一区二区三区在线观看| 又黄又爽又免费观看的视频| 久久久久久久久免费视频了| 精品国产一区二区三区四区第35| 制服诱惑二区| 国产精品永久免费网站| 国产精品1区2区在线观看.| 国产午夜精品久久久久久| 日日爽夜夜爽网站| 国产精品乱码一区二三区的特点 | 丰满的人妻完整版| 久久久国产成人免费| 亚洲国产精品合色在线| 在线天堂中文资源库| 多毛熟女@视频| 嫩草影视91久久| 一区二区三区精品91| 午夜日韩欧美国产| 亚洲全国av大片| 一a级毛片在线观看| 久久久国产成人免费| 亚洲欧美精品综合一区二区三区| 涩涩av久久男人的天堂| 欧美日韩福利视频一区二区| 黄色成人免费大全| 亚洲视频免费观看视频| 性欧美人与动物交配| 又黄又粗又硬又大视频| 成人国产一区最新在线观看| 国产一区在线观看成人免费| 丝袜人妻中文字幕| 亚洲国产日韩欧美精品在线观看 | 亚洲国产精品合色在线| 天堂动漫精品| 可以免费在线观看a视频的电影网站| 黑人巨大精品欧美一区二区mp4| 国产成人精品在线电影| 国产高清videossex| 久久国产精品男人的天堂亚洲| 国产欧美日韩综合在线一区二区| 国产av精品麻豆| 脱女人内裤的视频| 国产av在哪里看| 色精品久久人妻99蜜桃| 真人做人爱边吃奶动态| 国内毛片毛片毛片毛片毛片| 天天一区二区日本电影三级 | 精品第一国产精品| 男女床上黄色一级片免费看| 天天添夜夜摸| 久久久久久亚洲精品国产蜜桃av| 国产又爽黄色视频| 欧美日韩亚洲国产一区二区在线观看| 美女免费视频网站| √禁漫天堂资源中文www| 亚洲男人天堂网一区| 脱女人内裤的视频| 91成年电影在线观看| 亚洲va日本ⅴa欧美va伊人久久| 中国美女看黄片| 老司机靠b影院| 国产极品粉嫩免费观看在线| 午夜影院日韩av| 日韩欧美在线二视频| 长腿黑丝高跟| 精品久久久精品久久久| 亚洲七黄色美女视频| 亚洲三区欧美一区| 午夜老司机福利片| 亚洲国产看品久久| 久久热在线av| 视频区欧美日本亚洲| 欧美性长视频在线观看| 久久久久九九精品影院| ponron亚洲| 99久久综合精品五月天人人| 99热只有精品国产| 一级毛片高清免费大全| 国产高清有码在线观看视频 | 视频区欧美日本亚洲| 成人永久免费在线观看视频| 国产成人精品久久二区二区91|