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

    高空氣球發(fā)放過程中的流固耦合動(dòng)力學(xué)研究

    2023-09-04 14:36:02張向強(qiáng)
    計(jì)算機(jī)仿真 2023年7期
    關(guān)鍵詞:系纜氦氣浮力

    蘇 一,張向強(qiáng),劉 強(qiáng)

    (1. 中國科學(xué)院空天信息創(chuàng)新研究院,北京 100094;2. 中國科學(xué)院大學(xué),北京 100049)

    1 引言

    臨近空間(Near Space)通常被定義為距離地面18-100km高度的高空區(qū)域,其位于地球大氣中的平流層和中層,是一個(gè)待開拓的具有重要戰(zhàn)略意義的空間圈層。在這一空間圈層的中低層(18-50km)運(yùn)用最為廣泛的平臺(tái)是高空氣球,高空氣球根據(jù)球膜是否承受內(nèi)外壓差通常被劃分為零壓氣球和超壓氣球,從上世紀(jì)五六十年代開始,零壓氣球由于其成本低、載重大、發(fā)放簡易等優(yōu)勢,在空間科學(xué)領(lǐng)域方面發(fā)揮了重要的作用[5-7];1997年,NASA提出了ULDB計(jì)劃,其目的是完成超壓氣球載重1.6噸、駐空33.5km,續(xù)航100天的飛行任務(wù),目前已實(shí)現(xiàn)駐空54天的飛行[8]。由于臨近空間的大氣密度較小,海拔20km的大氣密度是地面的1/14,海拔35km的大氣密度是地面的1/150,即為維持浮力重力平衡,高空氣球的體積很大,同時(shí)材料面密度要小或厚度薄。目前國內(nèi)外常用的球膜材料一般為聚乙烯,體積達(dá)幾萬至幾十萬立方米,球膜厚度在20~80微米之間。若設(shè)計(jì)升限為40~50km,材料的厚度降為幾微米左右[13]。因此,基于氣球的形變和應(yīng)力分析的球形設(shè)計(jì)始終是高空氣球的難點(diǎn)[14]。

    高空氣球發(fā)放方式大體上可分為靜態(tài)發(fā)放、半動(dòng)態(tài)發(fā)放和動(dòng)態(tài)發(fā)放[2-4],不管采用哪種發(fā)放方式,解除約束后氣球受力和形態(tài)變化劇烈,典型的氣球發(fā)放過程如圖1所示。

    圖1 高空氣球發(fā)放過程

    2 發(fā)放過程系纜分析

    2.1 數(shù)學(xué)模型

    假設(shè)氣泡(球體充氣部分)剛性,結(jié)纜系統(tǒng)和氣泡以下的球體部分視為柔性系纜,簡化后的模型受力如下圖所示。

    考慮在X,Y坐標(biāo)下氣球某位置上球體微段的受力,它包括張力T和T′=T-dT,重量dW.g,風(fēng)對該線段的阻力dDr,在瞬時(shí)位置要求每一個(gè)微段受力平衡,在系纜沒有觸地的時(shí)候滿足以下平衡方程[1]:

    X方向有,

    T·cosθ+dDr·cosα=(T-dT)·cos(θ-dθ)

    (1)

    Y方向有,

    T·sinθ=dW·g+dDr·sinα+(T-dT)·sin(θ-dθ)

    (2)

    化簡得,

    dT=dDr·cosα·cosθ+(dW·g+dDr·sinα)·sinθ

    (3)

    另外有

    dx=cosθ·dl

    dy=sinθ·dl

    (4)

    (5)

    ODE45是一種自適應(yīng)步長的采用Runge-Kutta算法常微分方程數(shù)值解法,給定微分方程、積分范圍和參數(shù)初值可以得到積分范圍內(nèi)的解析解。積分的初始值為[T0θ0X0Y0],X0=0,Y0由高空氣球所處的位置決定。通過打靶法,并以二分法優(yōu)化初值的搜索過程,初值T0和θ0建立以下的關(guān)聯(lián)

    T0=F·sinθ0

    (6)

    其中,F為高空氣球產(chǎn)生的凈浮力,這樣T0的初始值由θ0決定,初始值變?yōu)閇F·sinθ0θ0X0Y0],變量只有一個(gè),確定θ0的取值范圍,通過打靶法找到確定的θ0,就可以完成曲線的擬合。

    假設(shè)有

    dDr=fDr·dl

    dW=ρ·dl

    (7)

    其中fDr系纜單位長度所受風(fēng)阻,ρ為線密度,根據(jù)(7)化簡(3),這樣積分時(shí)以dl為自變量,積分區(qū)間為系纜長度,根據(jù)公式,并考慮初值的正負(fù)性,則ODE45初值公式為

    =-fDr·cosα·cosθ+(ρ·g-fDr·sinα)·sinθ

    (8)

    其中α和θ之間有

    (9)

    當(dāng)系纜觸地后每一個(gè)微段會(huì)有

    T′+Ff=T

    (10)

    Ff=f·FN

    (11)

    根據(jù)圖2可知,觸地后系纜微段的微分方程為

    圖2 系纜受力圖

    (12)

    確定系纜的曲線后,考慮系纜是否觸地,將曲線做出后和地面的高度比較,判斷系纜曲線是否與地面直線相交,若不相交,則系纜曲線為所求;若相交,則分段打靶求曲線。第一段打靶根據(jù)初始條件找出滿足與地面相切的點(diǎn)C1,相切條件為目標(biāo)點(diǎn)C1(y=0,θ=0),然后根據(jù)觸地后的微分方程,根據(jù)區(qū)間設(shè)定點(diǎn)C2的范圍為(C2,O),作為第三段打靶的初始條件,打靶最后收斂于目標(biāo)點(diǎn)O,反過來確定C2的位置。

    高空氣球在發(fā)放過程中會(huì)受到來自空氣的阻力

    (13)

    其中,Fair為空氣阻力,ρa(bǔ)ir為空氣密度,V為高空氣球體積,Cd空氣阻力系數(shù),Sd為高空氣球受力面積。由上文可知,每次可以通過確定的位置得到當(dāng)時(shí)球體的受力情況,如圖3所示,確定系纜的曲線初值為[F·sinθ0θ0X0Y0],而只需要給定Y0就可以通過打靶法確定曲線,故只需分析Y方向的受力情況,根據(jù)基本的運(yùn)動(dòng)學(xué)公式

    圖3 高空氣球受力圖

    (14)

    Y方向受力情況為

    F=T·sinθ+Fair+(MHe+Mballoon)·a

    (15)

    本次數(shù)學(xué)模型采用VBalloon=7m3的小型高空氣球,充入30%的氦氣作為浮升氣體,經(jīng)過

    FHe=ρa(bǔ)ir·VHe·g

    GHe=ρHe·VHe·g

    (16)

    F=FHe-GHe-MBalloon·g

    (17)

    其中,VHe=0.3·VBalloon,經(jīng)過計(jì)算可知F=22.12N,因?yàn)槌錃饬看嬖谡`差并且球膜的質(zhì)量不可控,經(jīng)測量凈浮力為21.52N。

    通過MATLAB中ODE45函數(shù)求解后,高空氣球在高度Y方向的加速度、速度、位移和系纜末端O′點(diǎn)拉力隨時(shí)間的變化如圖4所示。

    圖4 高空氣球運(yùn)動(dòng)曲線

    由高空氣球運(yùn)動(dòng)軌跡可知,氣球在初始時(shí)加速度較大,速度在到達(dá)峰值后逐漸變小。

    2.2 有限元仿真

    通過MATLAB中的計(jì)算結(jié)果可知,在初始時(shí)刻發(fā)放車和高空球底部的直線距離為20.76m,在DYNA的前處理中將系纜未與高空氣球相連的一端通過添加初始位移的方法將距離調(diào)至相同,并通過添加重力達(dá)到發(fā)放的初始狀態(tài)的參數(shù)設(shè)置相同距離,除此之外,高空氣球高度、發(fā)放車高度、系纜的粗細(xì)、長度和線密度以及系纜所受的重力和高空氣球產(chǎn)生的浮力均與數(shù)學(xué)模型一致,在LS-DYNA中建立模型如圖5所示。

    圖5 高空氣球系纜隨時(shí)間變化(間隔1s)

    圖6 LS-DYNA模型初始狀態(tài)

    圖7 LS-DYNA仿真系纜形狀變化(間隔1s)

    高空氣球在上升中受浮力和重力相互影響,在LS-DYNA中取全局阻尼為3.5,計(jì)算結(jié)果如圖8~11所示。

    圖8 氣球加速度隨時(shí)間變化

    圖9 氣球速度隨時(shí)間變化

    圖10 氣球高度隨時(shí)間變化

    圖11 系纜拉力隨時(shí)間變化

    圖12 高空氣球受力分析

    圖13 高空氣球初始形狀

    通過LS-DYNA仿真可以看出,MATLAB和LS-DYNA計(jì)算結(jié)果說明高空氣球在發(fā)放過程中確實(shí)會(huì)受到系纜拉力作用,但LS-DYNA中無法做到根據(jù)用戶要求定義阻尼大小,故速度峰值會(huì)有所不同。

    3 發(fā)放仿真

    在LS-DYNA中進(jìn)行動(dòng)力學(xué)仿真,觀察球膜的形狀以及應(yīng)力變化,因此球膜需要轉(zhuǎn)化成柔性體,通過在球膜里面加入SPH粒子來模擬氦氣,通過SPH和球膜的相互作用,對發(fā)生沖擊時(shí)球膜和氣體進(jìn)行耦合分析。

    3.1 高空氣球球型模型

    氣球球膜一般由厚度為20μm左右的薄膜制成,球膜的厚度和球體尺度相比很小,符合無矩薄殼理論使用的條件[12]。

    由回轉(zhuǎn)薄殼無力矩理論可知:

    (18)

    其中

    (19)

    (20)

    將式(19)(20)代入(18),并化簡得

    (21)

    在截面圓上的總負(fù)載T=2πrσm在弧長方向的改變程度和周向應(yīng)力在弧長方向的分力以及膜重的和應(yīng)保持平衡,如式(22)所示。

    (22)

    化簡式(22)得

    (23)

    球形參數(shù)λ和∑的定義如式(24)所示,λ和Σ表示了氣球負(fù)載L及飛行高度對球形的影響程度

    (24)

    自然形高空氣球應(yīng)當(dāng)保證周向應(yīng)力σc=0,上式共同組成球形母線的求解微分方程組,在NASA高空氣球手冊中有Σ表,不同的Σ值對應(yīng)不同的母線形狀,本文采用的Σ=0.3,成型后的形狀圖14所示。

    圖14 高空氣球柔性球膜發(fā)放過程

    3.2 SPH方法

    SPH方法(Smoothed Particle Hydrodynamics)全稱為光滑粒子流體動(dòng)力學(xué)方法[10],是一種為求解流體的動(dòng)力學(xué)問題而提出的一種無網(wǎng)格方法。一般流體動(dòng)力學(xué)問題需要基于密度、速度、能量等變量場的偏微分方程組,而SPH方法是通過采用插值方法將連續(xù)的流體用相互作用的質(zhì)點(diǎn)組來描述,每個(gè)質(zhì)點(diǎn)上包含各種物理量的信息,這樣再過連續(xù)介質(zhì)力學(xué)的守恒定律可以將整個(gè)流體過程描述出來,從這個(gè)方面來說,只要采用的質(zhì)點(diǎn)粒子數(shù)夠多,就能準(zhǔn)確的描述流體的力學(xué)行為。

    在有限元計(jì)算中,SPH是一種純Lagrange方法,SPH由于不存在網(wǎng)格關(guān)系,所以不存在由于變形過大導(dǎo)致的網(wǎng)格畸變問題,因此適合求解高速碰撞過程中的大變形問題。

    3.3 接觸算法

    LS-dyna在處理接觸問題的過程中主要采用以下三種不同的算法:

    1)動(dòng)力約束法

    2)分配參數(shù)法

    3)對稱罰函數(shù)法。

    LS-DYNA的一般默認(rèn)采用對稱罰函數(shù)法,是計(jì)算中最常用的一種方法。

    其首先要確定接觸面的主從關(guān)系,一般從接觸面比主接觸面的網(wǎng)格更密,然后在計(jì)算過程中的每一個(gè)時(shí)間步都判斷從節(jié)點(diǎn)是否會(huì)穿透主表面,沒穿透則不操作;若穿透,則對從節(jié)點(diǎn)和接觸面之間引入一種限制穿透作用的類似法向彈簧的接觸力,接觸力的大小和穿透深度、接觸剛度有關(guān),具體如式(25)所示

    F=K·δ

    (25)

    其中,K為接觸界面剛度,δ為穿透量。這種接觸算法應(yīng)用簡單,不需要添加碰撞和釋放條件,而且計(jì)算過程中保持動(dòng)量守恒。在球膜與氦氣耦合過程中,氦氣為從節(jié)點(diǎn),球膜為主接觸面[9]。

    3.4 沖擊仿真

    3.4.1 材料參數(shù)

    表1 仿真參數(shù)表

    3.4.2 高空氣球發(fā)放沖擊模擬

    本文研究的高空氣球的發(fā)放過程是高空氣球在充氣完畢后在滾筒的約束下保持穩(wěn)定為發(fā)放零時(shí)刻,之后打開發(fā)放滾筒到氣球升起并在升到系纜的極限長度后在下方發(fā)放車的約束下停止并保持穩(wěn)定的一個(gè)過程。

    3.4.2.1 建立模型

    首先根據(jù)Σ=0.3做出球型的母線,將母線導(dǎo)入U(xiǎn)G建立模型,再用ICEM為模型畫網(wǎng)格,用LS-prepost進(jìn)行前處理,并加入SPH粒子,建立好的網(wǎng)格模型如下圖所示。其中高空氣球體積為7m3,沖入的氦氣為氣球體積的30%,然后對氦氣施加浮力,對球膜和系纜施加重力,獲得初始穩(wěn)定的形狀。

    利用文章2.2章節(jié)中所采用的仿真方法,如圖5所示,將球膜轉(zhuǎn)換成柔性體并沖入30%氦氣,會(huì)得到如圖15所示的球膜以及系纜曲線。

    圖15 高空氣球充氦氣后形狀與實(shí)際對比圖

    發(fā)放過程中會(huì)出現(xiàn)沖擊是在系纜拉直,球膜會(huì)向系纜方向收緊,和豎直向上所受沖擊模式類似,故可以采用豎直向上發(fā)放的方式研究受到系纜和氣體沖擊時(shí)球膜所受應(yīng)力變化。

    3.4.2.2 添加場力

    因?yàn)榻⒌哪P屯獠坎]有氣體,所以要讓氦氣產(chǎn)生向上的浮力氣體的密度并非氦氣的密度,而是通過浮力公式進(jìn)行換算,算出的密度應(yīng)當(dāng)將外部氣體和氦氣的密度都考慮進(jìn)去,根據(jù)式(16),有

    F凈浮力=FHe-GHe=ρa(bǔ)ir·VHe·g-ρHe·VHe·g

    =ρHe′·VHe·g

    (26)

    故有

    ρHe′=ρa(bǔ)ir-ρHe

    (27)

    這樣,預(yù)先充入的氣體在受到向上的加速度g時(shí),可以產(chǎn)生向上的且符合實(shí)際的浮力,球膜和系纜通過關(guān)鍵字*LOAD_GRAVITY_PART施加重力,獲得初始形狀如圖16所示。

    圖16 高空氣球沖擊受力圖

    高空氣球成型后的形狀是內(nèi)部流場和外部流場對球膜共同作用的結(jié)果,由于仿真后的高空氣球模型沒有加入外部流場,故來自于外部流場的壓強(qiáng)和風(fēng)場并沒有作用在球膜上,圖15中球膜的形狀是在內(nèi)部流場作用下穩(wěn)定的結(jié)果。

    3.4.2.3 后處理

    獲得初始穩(wěn)定形狀后,仿真模型在穩(wěn)定形狀的基礎(chǔ)上給其中的SPH粒子和球膜節(jié)點(diǎn)通過關(guān)鍵字*INITIAL_VELOSITY施加初速度4m/s,在仿真時(shí)不施加阻尼。 以下為沖擊過程中兩次沖擊間隔過程中球膜的應(yīng)力以及形狀變化。

    通過球膜的受力情況可知,開始時(shí)球頂中心有比較大的應(yīng)力,之后隨著氣體擴(kuò)散,球頂部受力比較均勻,但較大的應(yīng)力還是集中在球頂和氣球底部。當(dāng)氣體和球膜沖擊后會(huì)造成球膜有比較大形狀變化,并由于氣體運(yùn)動(dòng)的關(guān)系和球膜的接觸變少,所以此時(shí)球膜應(yīng)力并不大。在浮力的作用下,氣體再次上升,會(huì)和球膜再一次發(fā)生碰撞,這時(shí)球頂部和球底部應(yīng)力較大,通過應(yīng)力云圖發(fā)現(xiàn)球底部應(yīng)力較頂部更大。

    取高空氣球頂部中心的單元,做出其速度和應(yīng)力曲線如圖17所示。

    圖17 球頂?shù)乃俣扰c應(yīng)力隨時(shí)間變圖

    由圖中應(yīng)力和速度曲線可知,在速度下降的同時(shí)伴隨著球膜應(yīng)力的上升,說明氣體沖擊會(huì)造成球膜應(yīng)力增加,同時(shí)球膜會(huì)降低氣體的速度,說明球膜和氣體產(chǎn)生了耦合作用。

    4 結(jié)論

    本文以7m3高空氣球?yàn)槔?通過MATLAB和LS-DYNA對高空氣球的發(fā)放過程中系纜隨時(shí)間的變化以及受到?jīng)_擊時(shí)球膜形狀及應(yīng)力的變化的進(jìn)行了模擬。通過MATLAB的數(shù)值仿真可以得到系纜的拉力,并可以對高空氣球做受力分析,通過運(yùn)動(dòng)學(xué)方程得到每一時(shí)刻高空氣球的位置,并得到高空氣球每個(gè)位置的受力情況,通過LS-DYNA將簡化后的模型進(jìn)行了有限元仿真,并對受到?jīng)_擊時(shí)的球膜的形狀及應(yīng)力變化進(jìn)行了仿真,通過對比可以得到以下結(jié)論:

    1)高空氣球發(fā)放時(shí)會(huì)產(chǎn)生2~3g的加速度,這時(shí)候由于氦氣和球膜的初始狀態(tài)均是靜止的,相對速度不大,并且系纜并沒有拉緊,初始狀態(tài)更接近于高空氣球做自由上升運(yùn)動(dòng),系纜對球膜的拉力隨著高空氣球高度的增加逐漸增大,阻力也隨著球膜的速度的增大而增大,這樣高空氣球在阻力和系纜拉力的共同作用下速度上升不會(huì)太快,能有效減小到達(dá)沖擊點(diǎn)時(shí)的速度,達(dá)到減小沖擊的效果。

    2)通過LS-DYNA將發(fā)放時(shí)球膜所受沖擊進(jìn)行仿真,通過應(yīng)力云圖可知?dú)馇蝽敳亢蜌馇虻撞繛檩^大的應(yīng)力集中區(qū)域,而由于球頂部受力面積較大,應(yīng)力相對較小,而受力更加集中的底部反而是更容易遭到破壞的地方。在以后球頂部中心要加強(qiáng)材料強(qiáng)度,而球底部要加強(qiáng)對球膜的保護(hù),避免應(yīng)力集中。

    在以后計(jì)算高空氣球發(fā)放過程中,可以采用MATLAB數(shù)值仿真的方法,改變系纜長度,高空氣氣球的大小,和充氣量的不同,計(jì)算不同情況下高空氣球的速度峰值。通過LS-DYNA可以仿真小型高空氣球球膜應(yīng)力及形狀變化,但難以采用SPH方法對大型高空氣球進(jìn)行模擬,因?yàn)殡S著高空氣球的體積增大,網(wǎng)格數(shù)量和SPH粒子數(shù)量也在增加,導(dǎo)致計(jì)算量呈幾何級增長,需要找到更加合適的計(jì)算方法。

    猜你喜歡
    系纜氦氣浮力
    神奇的氦氣
    “浮力”知識(shí)鞏固
    我們一起來“制服”浮力
    浮力大小由誰定
    30萬t級油船靠泊過程實(shí)測系纜力統(tǒng)計(jì)分析
    水道港口(2022年1期)2022-05-18 13:01:20
    跟氣球上天
    廉政瞭望(2020年17期)2020-11-17 07:37:32
    坦桑尼亞發(fā)現(xiàn)巨型氦氣礦
    低溫與特氣(2018年1期)2018-04-16 13:19:36
    三峽船閘浮式系纜樁系泊安全性分析及實(shí)船試驗(yàn)
    船海工程(2018年1期)2018-03-01 02:47:04
    飛走的氦氣球能飛多高?
    神奇的浮力
    日本黄色日本黄色录像| 久久久久久久精品精品| 久久影院123| 老司机深夜福利视频在线观看 | 午夜福利一区二区在线看| 叶爱在线成人免费视频播放| 少妇人妻 视频| 久久国产精品人妻蜜桃| 黄网站色视频无遮挡免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲七黄色美女视频| 视频区图区小说| 国产日韩一区二区三区精品不卡| 激情视频va一区二区三区| 乱人伦中国视频| 在线亚洲精品国产二区图片欧美| 在线观看免费午夜福利视频| 亚洲精品自拍成人| 亚洲美女黄色视频免费看| 中文字幕人妻熟女乱码| 少妇人妻久久综合中文| 国产国语露脸激情在线看| 校园人妻丝袜中文字幕| 久久久久久久国产电影| 国产成人精品在线电影| 免费观看a级毛片全部| 男女高潮啪啪啪动态图| 一边摸一边做爽爽视频免费| 亚洲,欧美,日韩| 亚洲精品久久久久久婷婷小说| 国产亚洲午夜精品一区二区久久| 精品少妇一区二区三区视频日本电影| 欧美黑人精品巨大| 国产精品久久久久久人妻精品电影 | 一本久久精品| 亚洲伊人色综图| 日韩熟女老妇一区二区性免费视频| 69精品国产乱码久久久| xxxhd国产人妻xxx| 嫩草影视91久久| 黄色一级大片看看| 菩萨蛮人人尽说江南好唐韦庄| 久久 成人 亚洲| 国产无遮挡羞羞视频在线观看| 蜜桃在线观看..| 久久免费观看电影| 91老司机精品| 国产爽快片一区二区三区| 国产在线免费精品| 又大又黄又爽视频免费| 欧美日韩国产mv在线观看视频| 亚洲伊人色综图| 免费高清在线观看日韩| 国产男人的电影天堂91| 97精品久久久久久久久久精品| 两个人免费观看高清视频| 电影成人av| 一级毛片电影观看| 国产亚洲一区二区精品| 啦啦啦 在线观看视频| 欧美精品一区二区免费开放| 青青草视频在线视频观看| 国产亚洲一区二区精品| videosex国产| 精品久久蜜臀av无| 亚洲人成77777在线视频| 男女边摸边吃奶| 国产黄色视频一区二区在线观看| 黄色a级毛片大全视频| 国产一区亚洲一区在线观看| 操美女的视频在线观看| 性少妇av在线| 国产视频首页在线观看| 夫妻午夜视频| 精品第一国产精品| 欧美日韩国产mv在线观看视频| 97人妻天天添夜夜摸| 欧美 日韩 精品 国产| 午夜免费成人在线视频| 狠狠婷婷综合久久久久久88av| 女人精品久久久久毛片| 欧美精品啪啪一区二区三区 | 成年av动漫网址| 久久精品人人爽人人爽视色| 成年女人毛片免费观看观看9 | 一二三四社区在线视频社区8| 欧美变态另类bdsm刘玥| 欧美黄色淫秽网站| 如日韩欧美国产精品一区二区三区| 视频在线观看一区二区三区| 男人操女人黄网站| 大片电影免费在线观看免费| 高潮久久久久久久久久久不卡| 欧美国产精品va在线观看不卡| 亚洲精品成人av观看孕妇| 国产不卡av网站在线观看| 亚洲精品久久久久久婷婷小说| 性少妇av在线| 久久久久精品国产欧美久久久 | 成人亚洲欧美一区二区av| 永久免费av网站大全| av福利片在线| 黑人猛操日本美女一级片| 色婷婷久久久亚洲欧美| 欧美日韩黄片免| a级片在线免费高清观看视频| 国产色视频综合| av国产精品久久久久影院| 欧美日韩亚洲高清精品| 黄色一级大片看看| 爱豆传媒免费全集在线观看| 欧美成人午夜精品| 亚洲人成77777在线视频| 少妇的丰满在线观看| 亚洲国产精品一区三区| 99久久99久久久精品蜜桃| 老汉色av国产亚洲站长工具| 午夜福利免费观看在线| 人人妻人人爽人人添夜夜欢视频| 免费观看人在逋| 国产成人91sexporn| 午夜日韩欧美国产| 亚洲中文字幕日韩| 男人添女人高潮全过程视频| 十八禁网站网址无遮挡| 国产成人精品在线电影| 桃花免费在线播放| 国产激情久久老熟女| 丰满少妇做爰视频| 久久久久久人人人人人| 国产欧美日韩精品亚洲av| 国产黄频视频在线观看| 国产精品三级大全| 老司机亚洲免费影院| 免费在线观看影片大全网站 | 久久人人97超碰香蕉20202| 午夜免费成人在线视频| 亚洲国产精品999| 日韩av在线免费看完整版不卡| 菩萨蛮人人尽说江南好唐韦庄| 国产1区2区3区精品| 成年av动漫网址| 欧美日本中文国产一区发布| 伊人亚洲综合成人网| 国产精品免费大片| 久久久久网色| 国产亚洲精品第一综合不卡| 蜜桃在线观看..| 90打野战视频偷拍视频| 久久精品成人免费网站| 又粗又硬又长又爽又黄的视频| 精品高清国产在线一区| 日韩一区二区三区影片| av欧美777| 国产精品香港三级国产av潘金莲 | 日本av免费视频播放| 日韩伦理黄色片| 国产精品欧美亚洲77777| av电影中文网址| 国产亚洲av高清不卡| 国产野战对白在线观看| 日本a在线网址| 日本a在线网址| 成人手机av| 国产精品亚洲av一区麻豆| 天堂8中文在线网| 亚洲欧美一区二区三区国产| 性色av一级| 性色av一级| 欧美黄色片欧美黄色片| 久久国产精品大桥未久av| 日本a在线网址| 在线观看免费午夜福利视频| 亚洲精品乱久久久久久| 黄片播放在线免费| 亚洲国产精品成人久久小说| 国产高清视频在线播放一区 | 热99久久久久精品小说推荐| 免费女性裸体啪啪无遮挡网站| 日本一区二区免费在线视频| 国产精品一区二区在线不卡| 在线 av 中文字幕| 精品国产一区二区三区四区第35| www.av在线官网国产| 男的添女的下面高潮视频| 人人妻人人澡人人看| 在线观看免费午夜福利视频| 国产成人免费观看mmmm| 啦啦啦观看免费观看视频高清| 精品不卡国产一区二区三区| 神马国产精品三级电影在线观看 | 日韩 欧美 亚洲 中文字幕| 亚洲熟妇熟女久久| 一级作爱视频免费观看| 亚洲av成人一区二区三| 亚洲一区二区三区色噜噜| 亚洲欧美日韩无卡精品| 免费高清在线观看日韩| 亚洲成人久久爱视频| 国产99久久九九免费精品| 日韩中文字幕欧美一区二区| 欧美在线黄色| 99国产综合亚洲精品| 成年免费大片在线观看| 可以在线观看毛片的网站| 动漫黄色视频在线观看| 日韩欧美一区二区三区在线观看| 欧美性猛交黑人性爽| 欧美成狂野欧美在线观看| 亚洲人成电影免费在线| 日本三级黄在线观看| 欧美日韩精品网址| 搡老岳熟女国产| 亚洲av成人av| 男女之事视频高清在线观看| 午夜精品在线福利| 在线播放国产精品三级| 又黄又粗又硬又大视频| 久久久久精品国产欧美久久久| 一本综合久久免费| 亚洲精品中文字幕在线视频| 在线观看一区二区三区| 日本免费a在线| 老司机午夜十八禁免费视频| 成人一区二区视频在线观看| 亚洲国产欧美日韩在线播放| 亚洲成人国产一区在线观看| 国产精品久久久久久人妻精品电影| 久久香蕉激情| 精品久久久久久久久久免费视频| 不卡一级毛片| 国产成人系列免费观看| 亚洲人成网站高清观看| 国产又黄又爽又无遮挡在线| 免费在线观看视频国产中文字幕亚洲| 亚洲av成人不卡在线观看播放网| 听说在线观看完整版免费高清| 欧美激情高清一区二区三区| 久久国产精品男人的天堂亚洲| 亚洲九九香蕉| 国产极品粉嫩免费观看在线| 人人妻人人看人人澡| 少妇粗大呻吟视频| 深夜精品福利| 亚洲五月婷婷丁香| 欧洲精品卡2卡3卡4卡5卡区| 国产精品1区2区在线观看.| 无遮挡黄片免费观看| 听说在线观看完整版免费高清| 久久天躁狠狠躁夜夜2o2o| 欧美国产日韩亚洲一区| 精品高清国产在线一区| 大型黄色视频在线免费观看| 亚洲精品色激情综合| 可以在线观看的亚洲视频| 一本精品99久久精品77| 精品少妇一区二区三区视频日本电影| 丁香欧美五月| 亚洲精品中文字幕一二三四区| 后天国语完整版免费观看| 午夜成年电影在线免费观看| 男男h啪啪无遮挡| 在线观看www视频免费| 亚洲精品久久成人aⅴ小说| 啪啪无遮挡十八禁网站| 50天的宝宝边吃奶边哭怎么回事| 婷婷精品国产亚洲av| 国产爱豆传媒在线观看 | 国产成人系列免费观看| 天天添夜夜摸| 久久久精品国产亚洲av高清涩受| 看免费av毛片| 日韩欧美国产一区二区入口| 啦啦啦 在线观看视频| 色播亚洲综合网| 亚洲国产精品久久男人天堂| 亚洲中文字幕一区二区三区有码在线看 | 视频区欧美日本亚洲| 日本一区二区免费在线视频| 免费在线观看日本一区| 成人三级黄色视频| 99riav亚洲国产免费| 正在播放国产对白刺激| 欧美三级亚洲精品| 国产亚洲欧美在线一区二区| 国内精品久久久久久久电影| 精品久久久久久久人妻蜜臀av| 精品电影一区二区在线| 51午夜福利影视在线观看| 黄色丝袜av网址大全| 精品国产美女av久久久久小说| 巨乳人妻的诱惑在线观看| 国内毛片毛片毛片毛片毛片| 香蕉国产在线看| 中文字幕久久专区| 午夜免费成人在线视频| 99国产精品一区二区蜜桃av| xxx96com| 大型黄色视频在线免费观看| 性色av乱码一区二区三区2| 欧美中文日本在线观看视频| 大型黄色视频在线免费观看| 日本一区二区免费在线视频| 999久久久国产精品视频| 午夜免费激情av| 老熟妇乱子伦视频在线观看| 国产成人精品久久二区二区91| 亚洲午夜理论影院| 久久久国产精品麻豆| 我的亚洲天堂| 岛国视频午夜一区免费看| aaaaa片日本免费| 欧美激情久久久久久爽电影| 丝袜美腿诱惑在线| 搡老熟女国产l中国老女人| 日本精品一区二区三区蜜桃| 黑人欧美特级aaaaaa片| 好男人电影高清在线观看| 中文字幕另类日韩欧美亚洲嫩草| 看片在线看免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 给我免费播放毛片高清在线观看| 熟女电影av网| 十八禁人妻一区二区| 欧美精品亚洲一区二区| 精品久久久久久久末码| 午夜激情av网站| 满18在线观看网站| 亚洲最大成人中文| 亚洲av第一区精品v没综合| 韩国av一区二区三区四区| 女生性感内裤真人,穿戴方法视频| 国产在线精品亚洲第一网站| 亚洲成av片中文字幕在线观看| 中文字幕最新亚洲高清| 久久99热这里只有精品18| 特大巨黑吊av在线直播 | 一本一本综合久久| 欧美三级亚洲精品| 后天国语完整版免费观看| 日本黄色视频三级网站网址| 久久精品91无色码中文字幕| 欧美不卡视频在线免费观看 | 色综合欧美亚洲国产小说| 波多野结衣av一区二区av| 国产视频内射| 久久人人精品亚洲av| 亚洲欧美激情综合另类| 国产野战对白在线观看| 欧美中文日本在线观看视频| 国产真人三级小视频在线观看| 成年免费大片在线观看| 国产成人系列免费观看| 可以免费在线观看a视频的电影网站| 99久久久亚洲精品蜜臀av| 成年人黄色毛片网站| 国内精品久久久久久久电影| 午夜影院日韩av| 在线看三级毛片| 午夜福利在线在线| 亚洲一区高清亚洲精品| 丁香欧美五月| 国产精品久久电影中文字幕| 黄色女人牲交| 级片在线观看| 亚洲av美国av| 日日夜夜操网爽| 黄色视频不卡| 亚洲av成人一区二区三| 日韩欧美免费精品| 国产精品国产高清国产av| 一夜夜www| 欧美 亚洲 国产 日韩一| 精品国产乱码久久久久久男人| 好看av亚洲va欧美ⅴa在| 叶爱在线成人免费视频播放| 久久国产亚洲av麻豆专区| 国产黄a三级三级三级人| 我的亚洲天堂| 久久久久久久午夜电影| 精品国产国语对白av| 淫秽高清视频在线观看| 色播在线永久视频| 免费女性裸体啪啪无遮挡网站| 无限看片的www在线观看| 一夜夜www| 女性被躁到高潮视频| 免费人成视频x8x8入口观看| 在线观看免费视频日本深夜| 亚洲五月婷婷丁香| 亚洲成人久久性| 免费在线观看黄色视频的| ponron亚洲| 日韩有码中文字幕| 男女午夜视频在线观看| 色av中文字幕| 久久国产乱子伦精品免费另类| e午夜精品久久久久久久| 最好的美女福利视频网| 色综合欧美亚洲国产小说| 久久精品夜夜夜夜夜久久蜜豆 | 午夜免费鲁丝| 欧美日韩中文字幕国产精品一区二区三区| 中文资源天堂在线| 日日干狠狠操夜夜爽| 手机成人av网站| 日本 欧美在线| 侵犯人妻中文字幕一二三四区| 亚洲五月天丁香| 欧美成人午夜精品| 一级作爱视频免费观看| av中文乱码字幕在线| 两人在一起打扑克的视频| 欧美av亚洲av综合av国产av| 久久草成人影院| 国产在线观看jvid| 国产片内射在线| 久久精品国产清高在天天线| 午夜免费鲁丝| 黑丝袜美女国产一区| 老司机深夜福利视频在线观看| 欧美在线一区亚洲| 国产亚洲精品久久久久久毛片| 欧美成狂野欧美在线观看| 国产爱豆传媒在线观看 | 亚洲一区二区三区不卡视频| 精品免费久久久久久久清纯| 在线观看66精品国产| 久久性视频一级片| 国产一级毛片七仙女欲春2 | 男人舔女人下体高潮全视频| 亚洲欧洲精品一区二区精品久久久| 久热这里只有精品99| 91麻豆av在线| 午夜影院日韩av| 日本a在线网址| 精品免费久久久久久久清纯| 岛国视频午夜一区免费看| 色哟哟哟哟哟哟| 91av网站免费观看| 怎么达到女性高潮| 一级a爱片免费观看的视频| 久99久视频精品免费| 老司机在亚洲福利影院| 亚洲国产精品成人综合色| 欧美精品亚洲一区二区| 国产精品久久久久久人妻精品电影| 在线观看舔阴道视频| 国产精品爽爽va在线观看网站 | 99久久99久久久精品蜜桃| 精品久久久久久久毛片微露脸| 在线国产一区二区在线| 午夜福利免费观看在线| 亚洲精品av麻豆狂野| 欧美激情极品国产一区二区三区| 国产99久久九九免费精品| 1024香蕉在线观看| 国产精品久久电影中文字幕| 美女午夜性视频免费| 麻豆成人午夜福利视频| 亚洲电影在线观看av| 免费看十八禁软件| 亚洲国产精品合色在线| 亚洲七黄色美女视频| 少妇熟女aⅴ在线视频| 51午夜福利影视在线观看| 日本熟妇午夜| 欧美日韩乱码在线| 亚洲专区字幕在线| 国产一区二区在线av高清观看| 日本一本二区三区精品| 18禁国产床啪视频网站| 午夜久久久久精精品| 宅男免费午夜| 欧美大码av| 宅男免费午夜| 九色国产91popny在线| 99热只有精品国产| 淫妇啪啪啪对白视频| 99精品欧美一区二区三区四区| 在线播放国产精品三级| 久久精品aⅴ一区二区三区四区| 狂野欧美激情性xxxx| 97碰自拍视频| 久久这里只有精品19| 大型黄色视频在线免费观看| 中文字幕精品亚洲无线码一区 | 午夜福利免费观看在线| 一级毛片高清免费大全| 哪里可以看免费的av片| 亚洲av中文字字幕乱码综合 | 无限看片的www在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲av片天天在线观看| or卡值多少钱| 久久青草综合色| 日韩三级视频一区二区三区| 夜夜躁狠狠躁天天躁| 欧美激情 高清一区二区三区| av片东京热男人的天堂| 亚洲成人久久性| 在线观看午夜福利视频| 精品少妇一区二区三区视频日本电影| 一级片免费观看大全| 欧美精品亚洲一区二区| 免费看十八禁软件| 国产一级毛片七仙女欲春2 | 亚洲国产欧美网| 国产私拍福利视频在线观看| 黄片播放在线免费| 亚洲精品美女久久av网站| 99精品在免费线老司机午夜| 国产精品一区二区三区四区久久 | 亚洲全国av大片| bbb黄色大片| 日韩 欧美 亚洲 中文字幕| 久久午夜亚洲精品久久| 一级黄色大片毛片| 免费看十八禁软件| 欧美三级亚洲精品| 午夜精品在线福利| 999精品在线视频| 一本久久中文字幕| 黄网站色视频无遮挡免费观看| 一本久久中文字幕| 亚洲av五月六月丁香网| 国产精品久久久久久人妻精品电影| 欧美日韩精品网址| 男女视频在线观看网站免费 | 岛国在线观看网站| 亚洲五月婷婷丁香| 亚洲国产毛片av蜜桃av| 亚洲在线自拍视频| 久久久国产精品麻豆| 91在线观看av| 一夜夜www| 婷婷亚洲欧美| 丝袜人妻中文字幕| 老司机深夜福利视频在线观看| 欧美激情高清一区二区三区| 国产免费男女视频| 狠狠狠狠99中文字幕| 午夜激情av网站| 国产精品 国内视频| 黑人巨大精品欧美一区二区mp4| 免费看十八禁软件| 欧美三级亚洲精品| 在线观看免费午夜福利视频| 日本五十路高清| 最新在线观看一区二区三区| avwww免费| 校园春色视频在线观看| 精品久久久久久,| 99热这里只有精品一区 | 欧美一级毛片孕妇| 在线观看一区二区三区| 亚洲第一电影网av| 亚洲,欧美精品.| 男女那种视频在线观看| 天堂动漫精品| 国内精品久久久久精免费| 变态另类丝袜制服| 久久国产亚洲av麻豆专区| 久久中文字幕人妻熟女| 亚洲最大成人中文| 久久精品91无色码中文字幕| 国产午夜福利久久久久久| 三级毛片av免费| 中文字幕人成人乱码亚洲影| 高清在线国产一区| 免费在线观看完整版高清| 午夜免费鲁丝| 久久久久久国产a免费观看| 欧美 亚洲 国产 日韩一| 18美女黄网站色大片免费观看| 亚洲一区二区三区色噜噜| 97超级碰碰碰精品色视频在线观看| 美女免费视频网站| 日本三级黄在线观看| 欧美黑人欧美精品刺激| 日本三级黄在线观看| 久久狼人影院| 欧美久久黑人一区二区| 中文字幕人妻丝袜一区二区| 看片在线看免费视频| 久久中文字幕一级| 少妇熟女aⅴ在线视频| 亚洲国产欧美一区二区综合| 黑人巨大精品欧美一区二区mp4| 1024视频免费在线观看| 国产伦在线观看视频一区| 国产高清有码在线观看视频 | 亚洲av中文字字幕乱码综合 | 国产激情欧美一区二区| 亚洲九九香蕉| 激情在线观看视频在线高清| 麻豆av在线久日| 国产精品1区2区在线观看.| 嫁个100分男人电影在线观看| 亚洲精品中文字幕在线视频| av福利片在线| 悠悠久久av| 两人在一起打扑克的视频| 18禁黄网站禁片免费观看直播| 校园春色视频在线观看| 男女做爰动态图高潮gif福利片| 国产高清videossex| 一区二区日韩欧美中文字幕| 欧美日韩亚洲国产一区二区在线观看| 男女做爰动态图高潮gif福利片| 免费看日本二区| 久久国产精品男人的天堂亚洲| 国产高清videossex| 天堂√8在线中文| 亚洲自偷自拍图片 自拍|