• <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
    飛走的氦氣球能飛多高?
    神奇的浮力
    久久久久国产精品人妻aⅴ院 | 日本vs欧美在线观看视频| 99精品欧美一区二区三区四区| 下体分泌物呈黄色| 多毛熟女@视频| 最近最新中文字幕大全免费视频| 国产日韩欧美亚洲二区| 欧美日韩av久久| 黄色a级毛片大全视频| 别揉我奶头~嗯~啊~动态视频| 色94色欧美一区二区| 亚洲伊人色综图| 欧美午夜高清在线| 欧美成人午夜精品| 国产男女内射视频| 成人手机av| 人人妻人人爽人人添夜夜欢视频| ponron亚洲| 国精品久久久久久国模美| 欧美丝袜亚洲另类 | 老司机影院毛片| 亚洲国产欧美网| 欧美亚洲 丝袜 人妻 在线| 一级a爱片免费观看的视频| 成人影院久久| 国产亚洲精品久久久久5区| 成人特级黄色片久久久久久久| svipshipincom国产片| 满18在线观看网站| 极品少妇高潮喷水抽搐| 久久青草综合色| 如日韩欧美国产精品一区二区三区| 免费黄频网站在线观看国产| 国产在线精品亚洲第一网站| 又大又爽又粗| 欧美老熟妇乱子伦牲交| 国产av一区二区精品久久| 国产精品久久久久久精品古装| 亚洲免费av在线视频| 最近最新中文字幕大全免费视频| 波多野结衣av一区二区av| 韩国av一区二区三区四区| av片东京热男人的天堂| 久久精品成人免费网站| 成年人免费黄色播放视频| 新久久久久国产一级毛片| 精品卡一卡二卡四卡免费| 夜夜夜夜夜久久久久| 欧美激情高清一区二区三区| 巨乳人妻的诱惑在线观看| 在线观看www视频免费| 中文字幕av电影在线播放| 欧美日韩视频精品一区| 亚洲欧美一区二区三区黑人| 国产xxxxx性猛交| 少妇 在线观看| 怎么达到女性高潮| 在线av久久热| 视频在线观看一区二区三区| 国产一区二区三区在线臀色熟女 | 免费在线观看日本一区| 国产成人免费观看mmmm| 成在线人永久免费视频| 国产又色又爽无遮挡免费看| 精品电影一区二区在线| 香蕉久久夜色| 最近最新中文字幕大全电影3 | 国产日韩欧美亚洲二区| 国产亚洲精品久久久久久毛片 | 国产精品久久久av美女十八| 成人特级黄色片久久久久久久| 老司机深夜福利视频在线观看| 最新美女视频免费是黄的| 久久香蕉激情| 大香蕉久久成人网| 满18在线观看网站| 老司机影院毛片| 亚洲男人天堂网一区| 精品一区二区三卡| 亚洲少妇的诱惑av| 激情视频va一区二区三区| 热99国产精品久久久久久7| 成人18禁高潮啪啪吃奶动态图| aaaaa片日本免费| 成人黄色视频免费在线看| 亚洲黑人精品在线| 久久精品人人爽人人爽视色| av福利片在线| 亚洲精品久久成人aⅴ小说| 久久狼人影院| 国产欧美日韩一区二区三区在线| 欧美激情久久久久久爽电影 | 美女国产高潮福利片在线看| 国产精品欧美亚洲77777| 1024视频免费在线观看| 99在线人妻在线中文字幕 | 国产av精品麻豆| 欧美黑人精品巨大| 亚洲伊人色综图| 丝袜在线中文字幕| 亚洲自偷自拍图片 自拍| 最新美女视频免费是黄的| 在线观看免费午夜福利视频| 欧美在线一区亚洲| 亚洲精品国产精品久久久不卡| 日韩 欧美 亚洲 中文字幕| 悠悠久久av| 美国免费a级毛片| 制服人妻中文乱码| 极品人妻少妇av视频| 1024视频免费在线观看| 1024视频免费在线观看| 女性生殖器流出的白浆| 午夜福利欧美成人| 涩涩av久久男人的天堂| 欧美 日韩 精品 国产| 一级a爱视频在线免费观看| 国产精品久久久久久精品古装| 亚洲国产精品合色在线| 欧美精品av麻豆av| 人妻久久中文字幕网| 操美女的视频在线观看| 性少妇av在线| av福利片在线| 国产精品成人在线| 啪啪无遮挡十八禁网站| 性色av乱码一区二区三区2| 亚洲av日韩精品久久久久久密| 岛国毛片在线播放| 精品人妻熟女毛片av久久网站| 午夜福利,免费看| 宅男免费午夜| 搡老岳熟女国产| 色在线成人网| av天堂在线播放| 中文亚洲av片在线观看爽 | 国产91精品成人一区二区三区| 日韩欧美一区二区三区在线观看 | 欧美日韩中文字幕国产精品一区二区三区 | 嫁个100分男人电影在线观看| 亚洲欧美日韩另类电影网站| 亚洲男人天堂网一区| 动漫黄色视频在线观看| 欧美精品啪啪一区二区三区| 国产成人精品久久二区二区免费| 欧美日韩亚洲国产一区二区在线观看 | 精品一品国产午夜福利视频| 12—13女人毛片做爰片一| 老司机在亚洲福利影院| 咕卡用的链子| 亚洲 欧美一区二区三区| 少妇的丰满在线观看| 无人区码免费观看不卡| 青草久久国产| 男女免费视频国产| 成人亚洲精品一区在线观看| 中文字幕人妻丝袜一区二区| 国产av一区二区精品久久| 国产又色又爽无遮挡免费看| 国产xxxxx性猛交| 丁香欧美五月| 免费少妇av软件| tube8黄色片| 村上凉子中文字幕在线| 99久久人妻综合| 午夜免费成人在线视频| 精品一区二区三卡| 国产精品香港三级国产av潘金莲| 精品一品国产午夜福利视频| 看片在线看免费视频| 老鸭窝网址在线观看| 日韩欧美免费精品| 日韩成人在线观看一区二区三区| 男女午夜视频在线观看| 搡老岳熟女国产| 国产一区有黄有色的免费视频| 国产xxxxx性猛交| 日韩欧美一区视频在线观看| 黄色女人牲交| 精品人妻在线不人妻| 一区福利在线观看| 美女午夜性视频免费| 日韩欧美一区视频在线观看| 国产xxxxx性猛交| 80岁老熟妇乱子伦牲交| 在线观看免费视频网站a站| 日韩精品免费视频一区二区三区| 巨乳人妻的诱惑在线观看| 午夜福利乱码中文字幕| 亚洲精品一二三| 美女高潮喷水抽搐中文字幕| 天堂动漫精品| 日本黄色日本黄色录像| 国产精品自产拍在线观看55亚洲 | 亚洲精品国产色婷婷电影| 亚洲在线自拍视频| 黄色怎么调成土黄色| 美女福利国产在线| 亚洲av成人一区二区三| 久久精品国产综合久久久| 欧美激情久久久久久爽电影 | 亚洲精品粉嫩美女一区| 视频区欧美日本亚洲| 日本a在线网址| 欧美精品一区二区免费开放| 在线观看日韩欧美| 视频区图区小说| 在线av久久热| 久久久久视频综合| svipshipincom国产片| 丰满迷人的少妇在线观看| 亚洲精品av麻豆狂野| 在线观看舔阴道视频| av电影中文网址| 成年人午夜在线观看视频| 亚洲精品成人av观看孕妇| 丝袜在线中文字幕| 午夜影院日韩av| 久久亚洲精品不卡| 91精品三级在线观看| 国产男女超爽视频在线观看| 欧美日韩亚洲综合一区二区三区_| 99国产精品一区二区蜜桃av | 好看av亚洲va欧美ⅴa在| 国产精品1区2区在线观看. | 免费少妇av软件| 亚洲男人天堂网一区| 91成年电影在线观看| 亚洲av日韩在线播放| a在线观看视频网站| 热99re8久久精品国产| 久久久久国产一级毛片高清牌| 久久影院123| 9191精品国产免费久久| a级毛片黄视频| 国产亚洲一区二区精品| 国产精品影院久久| 亚洲国产看品久久| 日本精品一区二区三区蜜桃| 19禁男女啪啪无遮挡网站| 中文字幕人妻丝袜制服| 日韩人妻精品一区2区三区| 别揉我奶头~嗯~啊~动态视频| 日本撒尿小便嘘嘘汇集6| 国产精品久久电影中文字幕 | www日本在线高清视频| 黄色成人免费大全| 国产精品综合久久久久久久免费 | 久久久精品免费免费高清| 国产91精品成人一区二区三区| 99在线人妻在线中文字幕 | 人人澡人人妻人| 国产精品电影一区二区三区 | 日本黄色日本黄色录像| 国产伦人伦偷精品视频| 这个男人来自地球电影免费观看| 中文字幕人妻熟女乱码| 老司机午夜福利在线观看视频| 9色porny在线观看| 国产亚洲精品一区二区www | 女性生殖器流出的白浆| 18禁美女被吸乳视频| 久久精品成人免费网站| 无人区码免费观看不卡| 亚洲国产中文字幕在线视频| 日日夜夜操网爽| 人妻久久中文字幕网| 久久久国产一区二区| 窝窝影院91人妻| 岛国在线观看网站| 天堂俺去俺来也www色官网| 精品亚洲成国产av| 国产精品二区激情视频| 麻豆国产av国片精品| 中文字幕色久视频| 我的亚洲天堂| 不卡av一区二区三区| 亚洲精品国产区一区二| 国产在线观看jvid| 久久久水蜜桃国产精品网| 天天操日日干夜夜撸| 波多野结衣一区麻豆| 亚洲精品国产一区二区精华液| 一级毛片精品| 成年版毛片免费区| 亚洲精品在线观看二区| 午夜福利,免费看| 18禁美女被吸乳视频| 99热国产这里只有精品6| 成人18禁在线播放| 9191精品国产免费久久| 午夜两性在线视频| 在线观看一区二区三区激情| 日本黄色日本黄色录像| 免费人成视频x8x8入口观看| 精品久久久久久电影网| 国产高清国产精品国产三级| 亚洲av成人一区二区三| 国产99白浆流出| 国产xxxxx性猛交| 亚洲精品国产精品久久久不卡| 母亲3免费完整高清在线观看| 午夜老司机福利片| 国内久久婷婷六月综合欲色啪| 91九色精品人成在线观看| 激情视频va一区二区三区| 欧美在线黄色| 国产av精品麻豆| 69精品国产乱码久久久| 91在线观看av| 亚洲欧美激情在线| 久久国产精品男人的天堂亚洲| 欧美黄色片欧美黄色片| 亚洲七黄色美女视频| 视频区图区小说| 伦理电影免费视频| 99国产精品一区二区三区| 欧美+亚洲+日韩+国产| 欧美老熟妇乱子伦牲交| 日韩欧美免费精品| 亚洲人成电影免费在线| 免费av中文字幕在线| 国产成+人综合+亚洲专区| 欧美+亚洲+日韩+国产| a级片在线免费高清观看视频| 男女之事视频高清在线观看| 欧美日韩一级在线毛片| 国产精品一区二区免费欧美| 一级a爱视频在线免费观看| 美女午夜性视频免费| 国产精品自产拍在线观看55亚洲 | 好看av亚洲va欧美ⅴa在| 国产精品98久久久久久宅男小说| 亚洲伊人色综图| 搡老岳熟女国产| 午夜福利,免费看| 亚洲国产欧美一区二区综合| 两人在一起打扑克的视频| 丝袜美足系列| 一进一出抽搐gif免费好疼 | 成人国产一区最新在线观看| 深夜精品福利| 久久精品91无色码中文字幕| 校园春色视频在线观看| 亚洲午夜精品一区,二区,三区| 国产深夜福利视频在线观看| 国产欧美亚洲国产| 成人黄色视频免费在线看| www.熟女人妻精品国产| a级毛片黄视频| 天天影视国产精品| 国产亚洲精品久久久久5区| 国产精品亚洲av一区麻豆| 欧美日韩国产mv在线观看视频| 成人av一区二区三区在线看| 日本黄色日本黄色录像| 午夜精品久久久久久毛片777| 免费观看精品视频网站| ponron亚洲| 黄色视频不卡| 人人妻人人澡人人看| 热re99久久国产66热| 无遮挡黄片免费观看| 99在线人妻在线中文字幕 | 咕卡用的链子| 日韩免费av在线播放| 久久中文看片网| 欧美日韩成人在线一区二区| 国产激情久久老熟女| 亚洲色图av天堂| 搡老岳熟女国产| 满18在线观看网站| 老司机福利观看| 美国免费a级毛片| 国产精品影院久久| 岛国毛片在线播放| 日本一区二区免费在线视频| 久久国产亚洲av麻豆专区| 女人被躁到高潮嗷嗷叫费观| 国产欧美日韩一区二区三| 日韩人妻精品一区2区三区| 国产亚洲欧美精品永久| 在线国产一区二区在线| 日韩免费av在线播放| 在线天堂中文资源库| 黄片大片在线免费观看| 91成人精品电影| 一夜夜www| 50天的宝宝边吃奶边哭怎么回事| 国产在线观看jvid| 久99久视频精品免费| 美女国产高潮福利片在线看| 成年人免费黄色播放视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久国产精品影院| 精品久久久久久久毛片微露脸| 在线观看免费午夜福利视频| 丝袜美足系列| 国产精品久久久人人做人人爽| 无人区码免费观看不卡| 国产乱人伦免费视频| 日本欧美视频一区| 一边摸一边抽搐一进一小说 | 日本欧美视频一区| a级毛片在线看网站| 精品国产国语对白av| 欧美午夜高清在线| 成人国语在线视频| 国产99白浆流出| 国产成人精品久久二区二区91| 免费日韩欧美在线观看| 交换朋友夫妻互换小说| 侵犯人妻中文字幕一二三四区| 久久亚洲真实| 国产亚洲精品久久久久久毛片 | 天天躁夜夜躁狠狠躁躁| 性少妇av在线| 国产精品久久久av美女十八| 亚洲avbb在线观看| 欧美午夜高清在线| 黄色a级毛片大全视频| 国产欧美日韩精品亚洲av| av国产精品久久久久影院| 国产欧美亚洲国产| 9191精品国产免费久久| 亚洲精品中文字幕一二三四区| 色婷婷久久久亚洲欧美| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美日韩另类电影网站| 男人的好看免费观看在线视频 | avwww免费| 成人18禁在线播放| 99精品欧美一区二区三区四区| 看免费av毛片| 天天添夜夜摸| av不卡在线播放| 天天操日日干夜夜撸| 欧美不卡视频在线免费观看 | av网站在线播放免费| 精品国产一区二区三区四区第35| 国产片内射在线| 国产一卡二卡三卡精品| 欧美最黄视频在线播放免费 | 国产精品免费一区二区三区在线 | 亚洲片人在线观看| 伦理电影免费视频| 久久久久精品人妻al黑| 99久久99久久久精品蜜桃| 国产亚洲精品久久久久5区| 老司机亚洲免费影院| 日韩人妻精品一区2区三区| 亚洲欧美激情在线| 久99久视频精品免费| 精品人妻1区二区| 大码成人一级视频| 久久精品国产清高在天天线| 国产精品电影一区二区三区 | 97人妻天天添夜夜摸| 欧美亚洲 丝袜 人妻 在线| 免费在线观看影片大全网站| 色婷婷av一区二区三区视频| 国产淫语在线视频| 最近最新中文字幕大全电影3 | av网站在线播放免费| 热99久久久久精品小说推荐| 国产99久久九九免费精品| 成人亚洲精品一区在线观看| 精品人妻1区二区| 在线观看免费日韩欧美大片| 不卡一级毛片| 日本撒尿小便嘘嘘汇集6| 巨乳人妻的诱惑在线观看| 亚洲国产中文字幕在线视频| 国产精品二区激情视频| 91在线观看av| 久久久久久久久久久久大奶| 成人特级黄色片久久久久久久| 99精品久久久久人妻精品| 男女午夜视频在线观看| 黑人猛操日本美女一级片| 中文字幕人妻丝袜一区二区| 91在线观看av| 亚洲五月天丁香| 身体一侧抽搐| 女人精品久久久久毛片| 黄色女人牲交| 午夜免费鲁丝| 国产有黄有色有爽视频| 久久国产精品影院| 免费在线观看黄色视频的| 男女午夜视频在线观看| 黑人猛操日本美女一级片| 精品人妻熟女毛片av久久网站| 少妇猛男粗大的猛烈进出视频| 性色av乱码一区二区三区2| 少妇 在线观看| av福利片在线| 国产xxxxx性猛交| 国产片内射在线| 欧美激情极品国产一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 俄罗斯特黄特色一大片| 黄色成人免费大全| 99热国产这里只有精品6| tube8黄色片| 成人影院久久| 亚洲av电影在线进入| 曰老女人黄片| 国产av又大| 精品亚洲成a人片在线观看| svipshipincom国产片| 99在线人妻在线中文字幕 | 这个男人来自地球电影免费观看| 亚洲专区中文字幕在线| 一级片'在线观看视频| 欧美黑人欧美精品刺激| 日韩免费高清中文字幕av| 这个男人来自地球电影免费观看| 欧美+亚洲+日韩+国产| 人妻久久中文字幕网| 99精品欧美一区二区三区四区| 黑人操中国人逼视频| 欧美在线一区亚洲| 午夜视频精品福利| 狂野欧美激情性xxxx| 国产欧美日韩一区二区三| 成人精品一区二区免费| 最近最新中文字幕大全免费视频| 999精品在线视频| 久久久国产成人精品二区 | 亚洲精品一二三| 制服人妻中文乱码| 亚洲成av片中文字幕在线观看| 亚洲情色 制服丝袜| 国产精品二区激情视频| 亚洲精品久久午夜乱码| 19禁男女啪啪无遮挡网站| 一级毛片精品| 亚洲人成77777在线视频| 91九色精品人成在线观看| 国产有黄有色有爽视频| 久久精品亚洲熟妇少妇任你| 免费人成视频x8x8入口观看| 国产精品免费视频内射| 1024视频免费在线观看| 一区二区三区激情视频| 成人黄色视频免费在线看| 女人被躁到高潮嗷嗷叫费观| 国产成人av教育| 新久久久久国产一级毛片| 国产精品亚洲av一区麻豆| 美女午夜性视频免费| 日韩免费高清中文字幕av| 欧美日韩乱码在线| 欧美不卡视频在线免费观看 | 人人妻人人爽人人添夜夜欢视频| 久久国产精品人妻蜜桃| 国产免费现黄频在线看| 变态另类成人亚洲欧美熟女 | 中文欧美无线码| av天堂在线播放| 欧美大码av| 精品国产乱码久久久久久男人| 丰满饥渴人妻一区二区三| 国产成人欧美在线观看 | 黄色毛片三级朝国网站| 久久久久久久精品吃奶| 18禁观看日本| 国产片内射在线| 俄罗斯特黄特色一大片| 精品国产超薄肉色丝袜足j| 老司机午夜福利在线观看视频| 老汉色av国产亚洲站长工具| 欧美精品啪啪一区二区三区| 在线观看免费午夜福利视频| 国产av一区二区精品久久| netflix在线观看网站| 一边摸一边抽搐一进一小说 | 正在播放国产对白刺激| 极品教师在线免费播放| 久久久久久免费高清国产稀缺| 国产成人欧美在线观看 | 欧美 亚洲 国产 日韩一| 青草久久国产| 在线视频色国产色| 怎么达到女性高潮| 啪啪无遮挡十八禁网站| 亚洲全国av大片| 成人影院久久| 精品免费久久久久久久清纯 | 精品午夜福利视频在线观看一区| 国产欧美日韩精品亚洲av| 欧美日韩av久久| 免费在线观看亚洲国产| 亚洲男人天堂网一区| 又黄又爽又免费观看的视频| 无人区码免费观看不卡| a级片在线免费高清观看视频| 国产成+人综合+亚洲专区| 欧美人与性动交α欧美精品济南到| 在线观看免费视频网站a站| 91在线观看av| 亚洲色图 男人天堂 中文字幕| 成年动漫av网址| 久久精品国产综合久久久| 黄片大片在线免费观看| 美女午夜性视频免费| 在线国产一区二区在线| 在线观看免费日韩欧美大片| 亚洲国产看品久久| 亚洲国产欧美网| 欧美成狂野欧美在线观看| 亚洲片人在线观看| 欧美日韩国产mv在线观看视频| 后天国语完整版免费观看|