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

    微重力下板式貯箱內(nèi)液體晃動(dòng)性能研究*

    2017-11-21 04:28:30梁軍強(qiáng)劉錦濤
    航天控制 2017年5期
    關(guān)鍵詞:貯箱表面張力液面

    胡 齊 李 永 梁軍強(qiáng) 劉錦濤 宋 濤

    1.北京控制工程研究所,北京 100190 2.國家體育總局體育科學(xué)研究所, 北京100061

    微重力下板式貯箱內(nèi)液體晃動(dòng)性能研究*

    胡 齊1,2李 永1梁軍強(qiáng)1劉錦濤1宋 濤1

    1.北京控制工程研究所,北京 100190 2.國家體育總局體育科學(xué)研究所, 北京100061

    板式貯箱內(nèi)推進(jìn)劑管理裝置(PMD)以板式結(jié)構(gòu)部件為主,其抑制液體晃動(dòng)性能應(yīng)能夠滿足衛(wèi)星平臺高穩(wěn)定度和快速機(jī)動(dòng)的需求。針對微重力下板式貯箱內(nèi)液體晃動(dòng)性能,建立基于有限體積(VOF)方法的氣液兩相流流動(dòng)模型,通過數(shù)值模擬計(jì)算手段,比較分析研究了微重力環(huán)境下有防晃葉片和無防晃葉片的板式PMD晃動(dòng)性能。計(jì)算分析結(jié)果表明,有防晃葉片的板式PMD抑制液體晃動(dòng)性能明顯優(yōu)于無防晃葉片的板式PMD,使液面快速趨于穩(wěn)定。同時(shí),將計(jì)算分析結(jié)果與液體晃動(dòng)微重力試驗(yàn)結(jié)果進(jìn)行了分析比較,微重力試驗(yàn)結(jié)果與計(jì)算分析結(jié)果非常相近。研究結(jié)果表明,有防晃葉片的板式PMD具有更好的抑制液體晃動(dòng)功能,能夠顯著地抑制貯箱內(nèi)推進(jìn)劑晃動(dòng),以滿足衛(wèi)星平臺的高穩(wěn)定度和快速機(jī)動(dòng)的需求。

    液體晃動(dòng)特性;板式推進(jìn)劑管理裝置;微重力;防晃葉片

    板式貯箱內(nèi)部推進(jìn)劑管理裝置(PMD)以板式結(jié)構(gòu)為主,能夠?qū)崿F(xiàn)推進(jìn)劑在軌全管理,是一種更為先進(jìn)的貯箱。隨著航天技術(shù)的不斷發(fā)展,推進(jìn)劑重量占航天器總重量的比重不斷增大,對航天器的姿態(tài)控制精度與運(yùn)動(dòng)穩(wěn)定性的要求也越來越高,但是航天器貯箱內(nèi)液體晃動(dòng)對航天器姿態(tài)控制與運(yùn)動(dòng)穩(wěn)定性都有很大影響,液體可能出現(xiàn)劇烈晃動(dòng)[1],因此必須考慮貯箱液體晃動(dòng)的影響。但航天器貯箱多在微重力環(huán)境中工作,微重力環(huán)境中的液體晃動(dòng)特性與常重力明顯不同,因此開展微重力條件下板式貯箱內(nèi)液體晃動(dòng)研究具有非常重要的工程價(jià)值。目前工程上常用基于線性理論的空間等效擺模型和彈簧-質(zhì)量塊模型近似描述液體晃動(dòng),在線性假設(shè)條件下,這些模型可較好地模擬小幅晃動(dòng),但無法準(zhǔn)確模擬液體大幅晃動(dòng)[2-4]。為此,本文采用雷諾時(shí)均-流體體積(RANS-VOF)模擬方法,建立基于有限體積(VOF)方法的氣液兩相流流動(dòng)模型,對微重力下某一板式貯箱內(nèi)液體晃動(dòng)進(jìn)行三維數(shù)值模擬,分析比較了有防晃葉片和無防晃葉片的板式PMD液體晃動(dòng)性能的計(jì)算結(jié)果,并將計(jì)算結(jié)果與液體晃動(dòng)微重力試驗(yàn)結(jié)果進(jìn)行了分析比較。

    1 研究方法

    貯箱內(nèi)液體晃動(dòng)過程是氣液兩相流流動(dòng)過程。為求解板式貯箱內(nèi)液體晃動(dòng)過程,本文建立氣液兩相流體流動(dòng)過程的N-S方程,采用雷諾時(shí)均-流體體積(RANS-VOF)計(jì)算方法對貯箱內(nèi)氣液兩相流動(dòng)過程進(jìn)行數(shù)值模擬研究。

    1.1 基于VOF方法的氣液兩相流流動(dòng)數(shù)學(xué)模型

    1.1.1 Navier-Stokes方程

    瞬態(tài)、不可壓、不相溶、等溫、定粘度且存在表面張力的兩相流Navier-Stokes方程:

    (1)

    (2)

    (3)

    式中,u為速度,ρ為流體密度,F(xiàn)為表面張力,g為重力,μ為流體粘度,T代表轉(zhuǎn)置。式(1)為連續(xù)性方程,式(2)為動(dòng)量守恒方程。式(2)考慮了壓力、湍流、表面張力及重力對流動(dòng)的影響。

    1.1.2 VOF方法

    不相溶的氣液兩相,相與相之間存在界面。不同流體分子間作用力的差異將導(dǎo)致相界面上出現(xiàn)表面張力。利用體積分?jǐn)?shù)方法(volume of fluid,VOF)能夠只采用一組N-S方程就能對氣液兩相及相界面進(jìn)行描述,使問題簡化。位于相界面處的單元內(nèi)氣液兩相同時(shí)共存,非相界面單元內(nèi)只存在單相,如圖1所示。因此,引入k相的體積分?jǐn)?shù)εk,某個(gè)單元內(nèi)εk可表示為

    (4)

    整個(gè)流場單元可被劃分為3個(gè)區(qū)域:εk(cell)=0說明單元內(nèi)不存在k相流體;εk(cell)=1說明單元內(nèi)只存在k相流體;0<εk(cell)<1說明單元內(nèi)有相界面存在。

    圖1 氣液界面求解方法

    相界面的跟蹤是不相溶多相流數(shù)值模擬的關(guān)鍵問題之一?;赩OF方法可以對每一相的體積分?jǐn)?shù)建立輸運(yùn)方程來對相界面進(jìn)行跟蹤。

    (5)

    式中,Sεk是描述相間質(zhì)量轉(zhuǎn)移的物理量。本次計(jì)算由于不考慮燃油的蒸發(fā),因此相間不存在質(zhì)量轉(zhuǎn)移,即Sεk=0。

    ρ=εlρl+(1-εl)ρg

    (6)

    μ=εlμl+(1-εl)μg

    (7)

    式中,下標(biāo)l,g分別代表液相和氣相。

    實(shí)際的相界面結(jié)構(gòu)在流場作用下非常復(fù)雜,尤其是要精確求解三維空間上的界面曲率和表面張力是相當(dāng)困難的問題。Brackbill等[5]利用較為簡潔的連續(xù)表面力模型(continuum surface force,CSF)解決了表面曲率和表面張力求解的困難,最重要的是該方法能直接推廣到三維空間,計(jì)算量則增大不多。

    CSF模型將表面張力以體積力的形式附加在流體動(dòng)量方程中,具體的表達(dá)式為:

    (8)

    式中,k是表面曲率,是表面張力系數(shù),n是垂直表面的單位向量,F(xiàn)為液相體積分?jǐn)?shù)。

    1.2 數(shù)值計(jì)算模型

    1.2.1 計(jì)算網(wǎng)格模型

    本文研究的貯箱液體晃動(dòng)過程需要采用整個(gè)貯箱作為計(jì)算域。為簡化問題,去除了原貯箱結(jié)構(gòu)中的導(dǎo)流板、蓄液器等結(jié)構(gòu),保留防晃葉片結(jié)構(gòu),計(jì)算域整體采用六面體為主、五面體為輔的網(wǎng)格劃分方式,如圖2所示。在表面張力影響顯著的計(jì)算區(qū)域內(nèi)應(yīng)盡量使用六面體網(wǎng)格,提高數(shù)值預(yù)測精度。計(jì)算域內(nèi)網(wǎng)格主要由六面體網(wǎng)格構(gòu)成,五面體網(wǎng)格比例<5%。計(jì)算域的網(wǎng)格總數(shù)約80萬,壁面附近的網(wǎng)格尺度在0.5~4mm范圍內(nèi),并對網(wǎng)格無關(guān)性進(jìn)行了驗(yàn)證。

    圖2 貯箱內(nèi)液體晃動(dòng)計(jì)算域及網(wǎng)格劃分

    1.2.2 計(jì)算初始條件

    貯箱液體晃動(dòng)過程的仿真在重定位液面狀態(tài)穩(wěn)定的基礎(chǔ)上開展,因此將微重力環(huán)境下重定位過程的穩(wěn)定狀態(tài)作為貯箱液體晃動(dòng)計(jì)算的初始條件,液體晃動(dòng)計(jì)算的充液比為40%,如圖3所示,數(shù)值模擬貯箱在非沿軌跡工況條件下的液體晃動(dòng)過程,分析有無防晃葉片對板式貯箱液體晃動(dòng)特性的影響。計(jì)算中將貯箱壁面設(shè)為無限光滑,無滑移邊界條件,采用靜態(tài)接觸角來模擬固體壁面對氣液界面的吸附作用,同時(shí)對近壁面區(qū)域流動(dòng)采用標(biāo)準(zhǔn)壁面函數(shù)來模擬,靜態(tài)接觸角設(shè)為1°,液體選用推進(jìn)劑MON,微重力大小為1×10-5g0,方向沿-Y向,同時(shí)貯箱晃動(dòng)的起始點(diǎn)設(shè)置在非沿軌跡工況的起始點(diǎn),持續(xù)時(shí)間為15s,以X向角速度、Y和Z向線速度為主。

    圖3 晃動(dòng)計(jì)算貯箱內(nèi)初始液面分布狀態(tài)

    2 研究結(jié)果與分析

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

    圖4給出了有無防晃葉片情況下貯箱內(nèi)的液體晃動(dòng)氣液界面分布情況。在相同時(shí)刻下,左側(cè)圖為有防晃葉片的貯箱液體晃動(dòng)結(jié)果,右側(cè)圖為無防晃葉片的貯箱液體晃動(dòng)結(jié)果,氣液界面上還顯示了速度值。計(jì)算結(jié)果顯示,在非沿軌跡工況條件下貯箱晃動(dòng)時(shí)(0~15s),葉片的存在會強(qiáng)化其對液體的攪拌作用,并引起氣液界面存在較為復(fù)雜的結(jié)構(gòu)?;蝿?dòng)開始后,氣泡中心區(qū)域速度增加,液體在該區(qū)域存在較顯著的流動(dòng)?;蝿?dòng)開始后6s時(shí),無防晃葉片貯箱內(nèi)的流體流動(dòng)速度顯著高于有葉片貯箱。當(dāng)貯箱停止晃動(dòng)后,在有防晃葉片貯箱內(nèi),氣泡運(yùn)動(dòng)明顯被約束在貯箱頂部區(qū)域,液面基本上沒有出現(xiàn)破裂現(xiàn)象,葉片的存在明顯起到了阻尼和約束液體流動(dòng)的作用,同時(shí)葉片上的通孔有利于葉片之間的液體流通,使得液面不會出現(xiàn)破裂現(xiàn)象,也起到了一定阻尼作用。在無防晃葉片貯箱內(nèi),氣泡在液體慣性力作用下,從后側(cè)逐漸運(yùn)動(dòng)至貯箱底部,然后運(yùn)行經(jīng)過貯箱前部,并最終回到頂部區(qū)域,而且在此過程中液面出現(xiàn)較為明顯的破裂現(xiàn)象。

    圖4 液體晃動(dòng)過程氣液界面分布情況

    從圖5中可以看出,有防晃葉片時(shí),液體質(zhì)心在晃動(dòng)結(jié)束后迅速趨于穩(wěn)定,且在Y和Z向(波動(dòng)幅度遠(yuǎn)大于X向)的液體質(zhì)心分別在18s和30s之后均趨于穩(wěn)定;無防晃葉片時(shí),Y和Z向的液體質(zhì)心在60s時(shí)仍存在較大波動(dòng)。在貯箱運(yùn)動(dòng)結(jié)束后,有防晃葉片貯箱可迅速抑制液體晃動(dòng),而無防晃葉片貯箱液體衰減的速率顯著降低。

    從圖6中可以看出,有防晃葉片時(shí),晃動(dòng)力在晃動(dòng)結(jié)束后迅速趨于穩(wěn)定,且在Y和Z向(晃動(dòng)力遠(yuǎn)大于X向)的晃動(dòng)力分別在20s和18s之后趨于穩(wěn)定;無防晃葉片時(shí),Y和Z向的晃動(dòng)力在20s~45s范圍內(nèi)仍存在較大波動(dòng),尤其是Z向。在貯箱運(yùn)動(dòng)結(jié)束后,有防晃葉片貯箱由于可迅速抑制液體晃動(dòng)使得晃動(dòng)力急劇減小,快速趨于穩(wěn)定,而無防晃葉片貯箱晃動(dòng)力幅值仍比較大,而且晃動(dòng)力衰減的速率顯著降低。

    從圖7中可以看出,有防晃葉片時(shí),晃動(dòng)力矩在晃動(dòng)結(jié)束后迅速趨于穩(wěn)定,且在X向(晃動(dòng)力矩遠(yuǎn)大于Y向和Z向)的晃動(dòng)力矩在20s之后趨于穩(wěn)定;無防晃葉片時(shí),X向的晃動(dòng)力在20~45s范圍內(nèi)仍存在較大波動(dòng)。在貯箱運(yùn)動(dòng)結(jié)束后,有防晃葉片貯箱由于可迅速抑制液體晃動(dòng)使得晃動(dòng)力矩急劇減小而快速區(qū)域穩(wěn)定,無防晃葉片貯箱晃動(dòng)力矩幅值仍比較大,而且晃動(dòng)力衰減的速率顯著降低。

    2.2 計(jì)算結(jié)果與試驗(yàn)結(jié)果對比分析

    圖8顯示了數(shù)值計(jì)算結(jié)果與微重力試驗(yàn)結(jié)果的對照情況??梢钥闯?,微重力試驗(yàn)結(jié)果與計(jì)算分析結(jié)果非常相近,在晃動(dòng)過程中氣液界面分布狀態(tài)基本一樣,而且均未出現(xiàn)液面破裂現(xiàn)象,從而驗(yàn)證了本文中數(shù)值計(jì)算分析方法的可靠性和正確性。還可以看出,無防晃葉片的貯箱模型內(nèi)液面晃動(dòng)比有防晃葉片的貯箱模型內(nèi)液面晃動(dòng)劇烈一些,液面晃動(dòng)幅度也更大一些,這是防晃葉片起抑制液體晃動(dòng)作用的效果。微重力試驗(yàn)結(jié)果表明,有防晃葉片的板式PMD抑制液體晃動(dòng)效果顯著。

    圖5 有無防晃葉片液體質(zhì)心的影響

    圖6 有無防晃葉片晃動(dòng)力的影響

    3 結(jié)論與展望

    通過對微重力下板式貯箱內(nèi)液體晃動(dòng)性能的分析研究工作,得到以下結(jié)論:

    1)微重力試驗(yàn)結(jié)果與計(jì)算分析結(jié)果非常相近,驗(yàn)證了本文數(shù)值計(jì)算分析方法的可靠性和正確性,而且數(shù)值計(jì)算可以模擬試驗(yàn)無法模擬的一些邊界條件及微重力環(huán)境,能直接給出試驗(yàn)中測量較為困難的參數(shù),在工程應(yīng)用上有較大的實(shí)用價(jià)值;

    2)有防晃葉片的板式PMD具有更好的抑制液體晃動(dòng)功能,能夠顯著地抑制貯箱內(nèi)推進(jìn)劑晃動(dòng)以滿足衛(wèi)星平臺的高穩(wěn)定度和快速機(jī)動(dòng)的需求;

    3)防晃葉片對液體晃動(dòng)影響明顯,同時(shí)可在防晃葉片上設(shè)置一些通孔,這些均為將來板式貯箱PMD抑制液體晃動(dòng)的結(jié)構(gòu)設(shè)計(jì)提供有益的參考依據(jù)。

    圖7 有無防晃葉片對液體晃動(dòng)力矩的影響

    圖8 計(jì)算結(jié)果與試驗(yàn)結(jié)果對照

    參 考 文 獻(xiàn)

    [1] 曲廣吉. 衛(wèi)星動(dòng)力學(xué)工程[M]. 北京: 中國科學(xué)技術(shù)出版社, 2000.

    [2] Kana D D.Validated Spherical Pendulum Model or Rotary Liquid Slosh [R]. American Institute of Aeronautics and Astronautics,1988.

    [3] Choi J, Sarigul-Klijn N. Liquid Sloshing in Flexible Tanks under Thruster Firing for Orbit Adjustment [C]. AIAA Joint Propulsion Conference & Exhibit. Washinnton D.C.: AIAA, 2006.

    [4] Mitra S, Sinhamahpatra K P. Slosh Dynamics of Liquid-Filled Container with Submerged Components using Pressure-Based Finite Element Method [J]. Journal of Sound and Vibration, 2007, 304: 361-381.

    [5] Brackbill J U, Kothe D B , Zemach C. A Continuum Method for Modeling Surface Tension [J]. Journal of Computational Physics, 1992(100): 335-354.

    LiquidSloshingPerformanceinVaneTypeTankunderMicrogravity

    Hu Qi1,2,Li Yong1,Liang Junqiang1,Liu Jintao1,Song Tao1

    1.Beijing Institute of Control Engineering, Beijing 100190, China 2.China Institute of Sport Science,Beijing 100061,China

    Propellantmanagementdevice(PMD)invanetypetankmainlycomprisesofvanetypestructureparts,whoseperformanceofrestrainingliquidsloshingshouldsatisfyspacecraftrequirementsofhighstabilizationandfastorbitalmaneuver.Aimingatliquidsloshingperformanceinvanetypetankundermicrogravityenvironment,gas-liquidflowmodelbasedonthevolumeoffluidmethodwasputforward,andvianumericalsimulationliquidsloshingperformancesofvanetypePMDwithanti-sloshingbafflesandwithoutanti-sloshingbafflesinmicrogravitywereanalyzedandcompared.SimulationresultsrevealthatliquidsloshingperformanceofvanetypePMDwithanti-sloshingbafflesismarkedlysuperiorvanetypePMDwithoutanti-sloshingbafflesandthebafflesmakeliquidsurfacebecomestablefast.Thenbycomparingbetween

    *國家自然科學(xué)基金資助項(xiàng)目(51406010)

    resultsofmicrogravityexperimentsandresultsofnumericalsimulations,theyareverysimilar.Accordingtopresentresearch,vanetypePMDwithanti-sloshingbaffleshasbettereffectsonrestrainingliquidsloshingandisabletorestrainobservablypropellantsloshingintanksinordertosatisfyspacecraftrequirementsofhighstabilizationandfastorbitalmaneuver.

    Liquidsloshingperformance;Vanetypepropellantmanagementdevice;Microgravity;Anti-sloshingbaffles

    TP316.2

    A

    1006-3242(2017)05-0019-06

    2016-08-08

    胡齊(1985-),男,江西人,工程師,主要研究方向?yàn)楹教焱七M(jìn)技術(shù)與微重力下流體機(jī)理;李永(1977-),男,山東人,研究員,主要研究方向?yàn)楹教焱七M(jìn)技術(shù)、微重力下流體機(jī)理與先進(jìn)流動(dòng)測量技術(shù);梁軍強(qiáng)(1970-),男,河北人,研究員,主要研究方向?yàn)楹教炱魍七M(jìn)技術(shù);劉錦濤(1986-),男,山東人,工程師,主要研究方向?yàn)楹教炱魍七M(jìn)技術(shù)與微重力流動(dòng)理論;宋濤(1980-),男,河南人,工程師,主要研究方向?yàn)楹教炱魍七M(jìn)技術(shù)。

    猜你喜歡
    貯箱表面張力液面
    運(yùn)載火箭貯箱補(bǔ)償器結(jié)構(gòu)剛度的試驗(yàn)研究
    基于Surface Evolver的推進(jìn)劑貯箱氣液界面分析
    貯箱爆炸碎片初始速度及影響因素
    貯箱輕量化設(shè)計(jì)幾何參數(shù)優(yōu)化方法
    吸管“喝”水的秘密
    基于DCS自動(dòng)控制循環(huán)水液面的改造
    電子測試(2018年6期)2018-05-09 07:31:47
    神奇的表面張力
    小布老虎(2016年4期)2016-12-01 05:46:08
    MgO-B2O3-SiO2三元體系熔渣表面張力計(jì)算
    上海金屬(2016年2期)2016-11-23 05:34:45
    CaF2-CaO-Al2O3-MgO-SiO2渣系表面張力計(jì)算模型
    上海金屬(2014年3期)2014-12-19 13:09:06
    CaO-A12O3-TiO2熔渣表面張力計(jì)算模型
    上海金屬(2014年2期)2014-12-18 06:52:45
    多毛熟女@视频| www日本在线高清视频| 国产精品一区二区精品视频观看| 亚洲专区国产一区二区| 97碰自拍视频| 亚洲成人精品中文字幕电影| 久久中文看片网| 人人妻人人爽人人添夜夜欢视频| 黄色视频不卡| 久久香蕉国产精品| 久久久久久国产a免费观看| 长腿黑丝高跟| 亚洲性夜色夜夜综合| 亚洲欧美一区二区三区黑人| 亚洲七黄色美女视频| 最近最新免费中文字幕在线| 桃红色精品国产亚洲av| 精品欧美国产一区二区三| 精品国产一区二区久久| 中文字幕色久视频| 午夜久久久在线观看| 最近最新中文字幕大全电影3 | av天堂久久9| 国产精华一区二区三区| 亚洲人成电影观看| 高清毛片免费观看视频网站| 一区在线观看完整版| 亚洲成a人片在线一区二区| 成人三级做爰电影| 国产av一区在线观看免费| 国产精品一区二区免费欧美| 免费一级毛片在线播放高清视频 | 好男人电影高清在线观看| 久久久国产欧美日韩av| 啦啦啦免费观看视频1| 非洲黑人性xxxx精品又粗又长| 99精品久久久久人妻精品| 国产精品美女特级片免费视频播放器 | 成人亚洲精品一区在线观看| 久久人妻福利社区极品人妻图片| 国内毛片毛片毛片毛片毛片| 精品久久久久久久毛片微露脸| 国产成年人精品一区二区| 亚洲久久久国产精品| 天堂√8在线中文| 在线天堂中文资源库| 国产一区二区三区在线臀色熟女| 日韩视频一区二区在线观看| 国产精品爽爽va在线观看网站 | 国产av在哪里看| av天堂在线播放| 亚洲伊人色综图| 亚洲精品在线观看二区| 日韩精品青青久久久久久| 午夜福利高清视频| 妹子高潮喷水视频| 国产成人欧美| 亚洲中文字幕一区二区三区有码在线看 | 黄色片一级片一级黄色片| 午夜福利,免费看| 91字幕亚洲| 国产成+人综合+亚洲专区| 国产av又大| www.999成人在线观看| 纯流量卡能插随身wifi吗| 一级毛片精品| 欧美成人一区二区免费高清观看 | 欧美亚洲日本最大视频资源| 欧美绝顶高潮抽搐喷水| 精品欧美一区二区三区在线| 国产私拍福利视频在线观看| 91在线观看av| 十八禁人妻一区二区| 激情视频va一区二区三区| 性欧美人与动物交配| 国产高清视频在线播放一区| aaaaa片日本免费| 精品久久久久久,| 欧美日韩精品网址| 日韩大码丰满熟妇| 天天躁夜夜躁狠狠躁躁| 久久影院123| 18禁美女被吸乳视频| 国产一区在线观看成人免费| 咕卡用的链子| 一本久久中文字幕| 91麻豆精品激情在线观看国产| 亚洲国产高清在线一区二区三 | 精品国产超薄肉色丝袜足j| www.www免费av| 亚洲国产欧美日韩在线播放| 波多野结衣巨乳人妻| 亚洲人成77777在线视频| 嫩草影视91久久| 在线十欧美十亚洲十日本专区| aaaaa片日本免费| 一区二区日韩欧美中文字幕| 国产成人免费无遮挡视频| 满18在线观看网站| 国产xxxxx性猛交| 女人爽到高潮嗷嗷叫在线视频| 一进一出抽搐动态| avwww免费| 国产成人精品久久二区二区免费| 99国产综合亚洲精品| 九色亚洲精品在线播放| 亚洲av成人一区二区三| 亚洲欧美日韩高清在线视频| 一边摸一边抽搐一进一小说| 热re99久久国产66热| 国产精品美女特级片免费视频播放器 | 最新美女视频免费是黄的| 国产亚洲欧美98| 一边摸一边抽搐一进一小说| 欧美成人一区二区免费高清观看 | 一边摸一边做爽爽视频免费| 久久精品国产99精品国产亚洲性色 | 久久午夜综合久久蜜桃| 夜夜爽天天搞| 日本欧美视频一区| 变态另类成人亚洲欧美熟女 | 久久人人97超碰香蕉20202| 咕卡用的链子| 可以在线观看毛片的网站| 国产精品秋霞免费鲁丝片| 后天国语完整版免费观看| 日韩欧美免费精品| 日韩高清综合在线| 操出白浆在线播放| 国产单亲对白刺激| 久久伊人香网站| 精品国内亚洲2022精品成人| 老司机深夜福利视频在线观看| 9热在线视频观看99| 女人被躁到高潮嗷嗷叫费观| 波多野结衣高清无吗| 国产又色又爽无遮挡免费看| 在线免费观看的www视频| 成人亚洲精品一区在线观看| 久久精品国产亚洲av高清一级| 麻豆av在线久日| 精品国产一区二区久久| 狂野欧美激情性xxxx| 国产片内射在线| 波多野结衣巨乳人妻| 成年版毛片免费区| 国产99白浆流出| 成人av一区二区三区在线看| 日本欧美视频一区| 十八禁网站免费在线| 午夜亚洲福利在线播放| 黄片大片在线免费观看| 桃色一区二区三区在线观看| 国产日韩一区二区三区精品不卡| av电影中文网址| 亚洲成人久久性| 日韩精品免费视频一区二区三区| 中文字幕人妻丝袜一区二区| 制服丝袜大香蕉在线| 如日韩欧美国产精品一区二区三区| 午夜老司机福利片| 久久午夜综合久久蜜桃| 亚洲av成人一区二区三| 国产精品一区二区精品视频观看| 欧美一区二区精品小视频在线| 亚洲国产毛片av蜜桃av| 久久中文看片网| 91av网站免费观看| 麻豆av在线久日| 99在线视频只有这里精品首页| 亚洲人成电影免费在线| 免费不卡黄色视频| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美精品综合久久99| 最近最新中文字幕大全电影3 | 别揉我奶头~嗯~啊~动态视频| 亚洲第一电影网av| 国语自产精品视频在线第100页| 人人妻,人人澡人人爽秒播| 久久久国产精品麻豆| 亚洲人成伊人成综合网2020| 亚洲国产精品久久男人天堂| 国产亚洲精品综合一区在线观看 | 久9热在线精品视频| 最近最新中文字幕大全电影3 | 99在线视频只有这里精品首页| 69精品国产乱码久久久| 18禁黄网站禁片午夜丰满| 免费不卡黄色视频| 久久久久久亚洲精品国产蜜桃av| 国产精品久久久人人做人人爽| 无遮挡黄片免费观看| 午夜影院日韩av| 日韩欧美一区视频在线观看| 欧美黄色片欧美黄色片| 亚洲中文av在线| 色精品久久人妻99蜜桃| 成熟少妇高潮喷水视频| 久久精品国产综合久久久| 岛国在线观看网站| 亚洲熟女毛片儿| 悠悠久久av| 日本欧美视频一区| 深夜精品福利| 一个人免费在线观看的高清视频| 高清黄色对白视频在线免费看| 欧美成人一区二区免费高清观看 | 亚洲av熟女| 十八禁人妻一区二区| 久久久国产成人免费| e午夜精品久久久久久久| 国产不卡一卡二| 国产在线观看jvid| av欧美777| 国产精品二区激情视频| 黄网站色视频无遮挡免费观看| 精品一区二区三区视频在线观看免费| 老汉色∧v一级毛片| 免费看a级黄色片| 99re在线观看精品视频| 午夜福利一区二区在线看| www.999成人在线观看| 精品久久久久久久毛片微露脸| 最近最新中文字幕大全免费视频| 亚洲精品久久国产高清桃花| 精品欧美国产一区二区三| bbb黄色大片| 欧美一级毛片孕妇| 精品第一国产精品| 亚洲国产欧美网| 又黄又爽又免费观看的视频| 在线观看www视频免费| 国产精品久久久久久亚洲av鲁大| 人人妻,人人澡人人爽秒播| 啦啦啦 在线观看视频| 日韩有码中文字幕| 91精品国产国语对白视频| 999久久久国产精品视频| 精品久久久久久,| 一区福利在线观看| 变态另类丝袜制服| 日本免费a在线| 欧美成人一区二区免费高清观看 | 久久久久九九精品影院| 久久久精品国产亚洲av高清涩受| 精品久久蜜臀av无| 91老司机精品| 搡老熟女国产l中国老女人| 九色亚洲精品在线播放| 久久人人97超碰香蕉20202| 精品国产一区二区三区四区第35| 深夜精品福利| 欧美日韩乱码在线| 一本久久中文字幕| 十分钟在线观看高清视频www| 国产精品 欧美亚洲| 国产1区2区3区精品| 免费无遮挡裸体视频| 久久婷婷成人综合色麻豆| 午夜免费成人在线视频| 欧美午夜高清在线| 亚洲av日韩精品久久久久久密| 久久亚洲真实| av福利片在线| 多毛熟女@视频| 99久久99久久久精品蜜桃| 亚洲欧美激情在线| 亚洲美女黄片视频| 岛国视频午夜一区免费看| 亚洲性夜色夜夜综合| 麻豆av在线久日| 女人精品久久久久毛片| 淫妇啪啪啪对白视频| 老司机深夜福利视频在线观看| 妹子高潮喷水视频| 美女免费视频网站| 成人免费观看视频高清| 国产精品久久视频播放| 女人被狂操c到高潮| 国产1区2区3区精品| 麻豆成人av在线观看| 亚洲美女黄片视频| 黄色视频不卡| 亚洲国产精品成人综合色| 精品国产一区二区久久| 中文字幕最新亚洲高清| 搞女人的毛片| 国产高清视频在线播放一区| 久久天堂一区二区三区四区| 久久中文字幕一级| 波多野结衣一区麻豆| 亚洲精品国产一区二区精华液| 真人做人爱边吃奶动态| 欧美成人免费av一区二区三区| 伦理电影免费视频| 国产成+人综合+亚洲专区| 身体一侧抽搐| 日韩欧美一区视频在线观看| 国产精品爽爽va在线观看网站 | avwww免费| 满18在线观看网站| 在线播放国产精品三级| 午夜日韩欧美国产| 久久精品91蜜桃| 国产精品香港三级国产av潘金莲| 成人免费观看视频高清| 一区二区三区精品91| 99久久综合精品五月天人人| 亚洲色图综合在线观看| 亚洲中文字幕日韩| 一a级毛片在线观看| 国产av精品麻豆| 国产色视频综合| 在线播放国产精品三级| 亚洲最大成人中文| 嫩草影视91久久| 99在线人妻在线中文字幕| 日韩大码丰满熟妇| 成人三级做爰电影| 一区二区三区国产精品乱码| 国产99白浆流出| 久久天堂一区二区三区四区| 色播亚洲综合网| 老司机午夜福利在线观看视频| 国产亚洲精品综合一区在线观看 | 一级作爱视频免费观看| 国产成人精品久久二区二区免费| 国产欧美日韩一区二区三区在线| 久久精品人人爽人人爽视色| 国产区一区二久久| www日本在线高清视频| 少妇被粗大的猛进出69影院| 久久天躁狠狠躁夜夜2o2o| 欧洲精品卡2卡3卡4卡5卡区| 一二三四在线观看免费中文在| avwww免费| 99riav亚洲国产免费| 我的亚洲天堂| 97碰自拍视频| 国产亚洲精品久久久久5区| 午夜福利成人在线免费观看| 久久精品91蜜桃| 国产午夜精品久久久久久| e午夜精品久久久久久久| 色老头精品视频在线观看| 1024香蕉在线观看| 桃红色精品国产亚洲av| 脱女人内裤的视频| 久久欧美精品欧美久久欧美| 国产精品精品国产色婷婷| 成人国语在线视频| 国产精品免费视频内射| 午夜福利一区二区在线看| 免费看十八禁软件| 两个人看的免费小视频| 日韩精品青青久久久久久| 91成年电影在线观看| 老司机午夜十八禁免费视频| 亚洲av片天天在线观看| 亚洲欧美日韩无卡精品| 中文字幕另类日韩欧美亚洲嫩草| 精品欧美一区二区三区在线| 成人国产综合亚洲| 日韩欧美一区二区三区在线观看| 久久久久九九精品影院| 脱女人内裤的视频| av天堂在线播放| 久久人妻熟女aⅴ| 成人三级做爰电影| 亚洲三区欧美一区| 亚洲熟妇熟女久久| 欧美日韩福利视频一区二区| 一本综合久久免费| 日韩欧美一区二区三区在线观看| 久久九九热精品免费| 中文字幕av电影在线播放| 好男人电影高清在线观看| 国产精品av久久久久免费| 精品无人区乱码1区二区| 精品卡一卡二卡四卡免费| 久久人人精品亚洲av| 国产成人精品无人区| 法律面前人人平等表现在哪些方面| 在线观看66精品国产| 宅男免费午夜| 久9热在线精品视频| 精品卡一卡二卡四卡免费| 丁香欧美五月| 母亲3免费完整高清在线观看| 正在播放国产对白刺激| 纯流量卡能插随身wifi吗| 免费人成视频x8x8入口观看| 国产高清videossex| 亚洲中文字幕日韩| 99精品欧美一区二区三区四区| 很黄的视频免费| 在线观看舔阴道视频| 男人舔女人下体高潮全视频| 99在线视频只有这里精品首页| 国产亚洲精品综合一区在线观看 | bbb黄色大片| 久久精品亚洲精品国产色婷小说| 国产精品香港三级国产av潘金莲| 制服诱惑二区| 看黄色毛片网站| 日韩 欧美 亚洲 中文字幕| 免费看十八禁软件| 亚洲中文日韩欧美视频| 久久国产乱子伦精品免费另类| 男女下面插进去视频免费观看| 淫妇啪啪啪对白视频| 亚洲av成人不卡在线观看播放网| 日日干狠狠操夜夜爽| 人妻丰满熟妇av一区二区三区| 最好的美女福利视频网| 亚洲国产精品成人综合色| 99国产精品99久久久久| 一a级毛片在线观看| 亚洲成av人片免费观看| 国产熟女xx| 欧美激情 高清一区二区三区| 免费不卡黄色视频| 国产亚洲精品久久久久久毛片| 成在线人永久免费视频| 亚洲欧洲精品一区二区精品久久久| 亚洲精品中文字幕在线视频| 中文字幕精品免费在线观看视频| 午夜a级毛片| 大香蕉久久成人网| 亚洲国产欧美网| 国产欧美日韩一区二区三| 国产欧美日韩综合在线一区二区| 国产亚洲精品av在线| 午夜福利欧美成人| 久久久国产成人免费| 精品福利观看| 国内毛片毛片毛片毛片毛片| 亚洲情色 制服丝袜| 美女大奶头视频| 久久精品成人免费网站| 一区二区三区国产精品乱码| 国产xxxxx性猛交| 69精品国产乱码久久久| 免费无遮挡裸体视频| 日韩三级视频一区二区三区| 亚洲欧美精品综合久久99| 午夜福利高清视频| cao死你这个sao货| 看黄色毛片网站| 国产蜜桃级精品一区二区三区| 一区在线观看完整版| 欧美乱色亚洲激情| 男女床上黄色一级片免费看| 国产精品秋霞免费鲁丝片| 亚洲欧美日韩另类电影网站| 午夜福利视频1000在线观看 | 纯流量卡能插随身wifi吗| 亚洲精品中文字幕一二三四区| 免费在线观看日本一区| 中文字幕人成人乱码亚洲影| 在线观看免费午夜福利视频| 午夜激情av网站| 精品久久久久久久毛片微露脸| 国产高清videossex| 黑人巨大精品欧美一区二区蜜桃| 久久天躁狠狠躁夜夜2o2o| netflix在线观看网站| 丰满的人妻完整版| 国产区一区二久久| 久久欧美精品欧美久久欧美| 99精品久久久久人妻精品| 妹子高潮喷水视频| 18美女黄网站色大片免费观看| 黑人操中国人逼视频| 变态另类丝袜制服| 亚洲国产中文字幕在线视频| 97碰自拍视频| 久久香蕉激情| 日本黄色视频三级网站网址| 电影成人av| 一本大道久久a久久精品| 国产精品二区激情视频| 婷婷丁香在线五月| 日本免费一区二区三区高清不卡 | 看黄色毛片网站| 日韩欧美一区二区三区在线观看| 热99re8久久精品国产| 老司机在亚洲福利影院| 久久久久久大精品| 亚洲第一电影网av| 人人妻人人澡欧美一区二区 | 久久久久国产精品人妻aⅴ院| 国产欧美日韩精品亚洲av| 老司机午夜福利在线观看视频| 少妇粗大呻吟视频| 最近最新免费中文字幕在线| 丝袜美足系列| 亚洲成人国产一区在线观看| 最新在线观看一区二区三区| 亚洲电影在线观看av| 亚洲精品中文字幕一二三四区| 亚洲精品中文字幕一二三四区| 老司机福利观看| 婷婷精品国产亚洲av在线| 青草久久国产| 午夜激情av网站| 日韩欧美免费精品| 国产精品电影一区二区三区| 他把我摸到了高潮在线观看| 人妻久久中文字幕网| 亚洲精品一区av在线观看| 两性夫妻黄色片| 亚洲欧美精品综合一区二区三区| 日本免费a在线| 99国产精品99久久久久| 深夜精品福利| 久久午夜综合久久蜜桃| av在线播放免费不卡| 国产成人欧美在线观看| 乱人伦中国视频| 亚洲人成电影免费在线| 欧美日韩黄片免| 欧美成狂野欧美在线观看| 在线观看www视频免费| 日本a在线网址| 一级a爱视频在线免费观看| 91国产中文字幕| 超碰成人久久| 亚洲精品中文字幕一二三四区| 免费无遮挡裸体视频| 麻豆av在线久日| 日韩中文字幕欧美一区二区| 亚洲五月天丁香| 岛国在线观看网站| 免费观看人在逋| 午夜福利影视在线免费观看| 叶爱在线成人免费视频播放| 咕卡用的链子| 日本免费a在线| 两性午夜刺激爽爽歪歪视频在线观看 | 一二三四社区在线视频社区8| 日本三级黄在线观看| 欧美日本亚洲视频在线播放| 亚洲人成电影免费在线| 精品国产亚洲在线| 色精品久久人妻99蜜桃| 欧美精品啪啪一区二区三区| 成在线人永久免费视频| 免费一级毛片在线播放高清视频 | 一卡2卡三卡四卡精品乱码亚洲| 自线自在国产av| 69精品国产乱码久久久| 亚洲成人免费电影在线观看| 人成视频在线观看免费观看| 熟妇人妻久久中文字幕3abv| 久久精品人人爽人人爽视色| 精品国内亚洲2022精品成人| 欧美激情 高清一区二区三区| 日本在线视频免费播放| 国产一区在线观看成人免费| av片东京热男人的天堂| 天天躁夜夜躁狠狠躁躁| 男人操女人黄网站| 一边摸一边做爽爽视频免费| 一级片免费观看大全| 久久久久国产精品人妻aⅴ院| 久久人妻福利社区极品人妻图片| 男男h啪啪无遮挡| 亚洲精品国产精品久久久不卡| 高潮久久久久久久久久久不卡| 在线永久观看黄色视频| 一卡2卡三卡四卡精品乱码亚洲| 日韩成人在线观看一区二区三区| 国产成人av教育| 大型黄色视频在线免费观看| 国产单亲对白刺激| 999久久久国产精品视频| 午夜精品久久久久久毛片777| 亚洲精品中文字幕在线视频| 桃红色精品国产亚洲av| 一级黄色大片毛片| 母亲3免费完整高清在线观看| 美女国产高潮福利片在线看| а√天堂www在线а√下载| 又黄又粗又硬又大视频| 久久中文字幕一级| 不卡av一区二区三区| 久久天躁狠狠躁夜夜2o2o| 精品不卡国产一区二区三区| 欧美日韩福利视频一区二区| 成人三级黄色视频| 国产片内射在线| 国产欧美日韩精品亚洲av| 日日爽夜夜爽网站| 久久久水蜜桃国产精品网| 国产激情欧美一区二区| 日韩三级视频一区二区三区| 国产99久久九九免费精品| 99在线人妻在线中文字幕| 国产精品电影一区二区三区| 日韩三级视频一区二区三区| 日日夜夜操网爽| 美女高潮到喷水免费观看| 中文字幕av电影在线播放| 国产亚洲欧美在线一区二区| 亚洲中文日韩欧美视频| 精品久久久久久久毛片微露脸| 99国产精品99久久久久| 国产精品爽爽va在线观看网站 | 在线观看午夜福利视频| 一边摸一边做爽爽视频免费| 亚洲人成电影观看| 精品国产亚洲在线| 一区在线观看完整版| 日韩免费av在线播放|