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

    動力翼傘系統(tǒng)擬坐標(biāo)形式的多體動力學(xué)建模

    2019-10-24 05:20:14張青斌高峰郭銳豐志偉葛健全
    兵工學(xué)報(bào) 2019年9期
    關(guān)鍵詞:傘衣吊帶氣動力

    張青斌, 高峰, 郭銳, 豐志偉, 葛健全

    (1.國防科技大學(xué) 空天科學(xué)學(xué)院, 湖南 長沙 410073; 2.96796部隊(duì), 吉林 白山 135200)

    0 引言

    翼傘系統(tǒng)不僅具有折疊體積小、質(zhì)量輕的優(yōu)點(diǎn),而且擁有優(yōu)良的操穩(wěn)和滑翔性能,可以在一定初始偏差和風(fēng)場等干擾條件下實(shí)現(xiàn)定點(diǎn)無損著陸,廣泛應(yīng)用于高空投放、遠(yuǎn)距離滲透和精確著陸等任務(wù)中,具有極大的軍事和民用價(jià)值。目前國外研發(fā)了一系列頗具代表性的空投系統(tǒng),例如美國Atair Aerospace 公司研制的自主導(dǎo)航空投系統(tǒng)、美國Strong Enterprises公司研制的Screamer空投系統(tǒng)、美國Alabama大學(xué)和海軍研究院聯(lián)合研發(fā)的Snowflake空投系統(tǒng),以及加拿大米斯特機(jī)動綜合系統(tǒng)技術(shù)公司研制的Sherpa系統(tǒng)和德國Cassidian公司研制的Paralander空投系統(tǒng)等[1]。上述空投系統(tǒng)能完成從數(shù)千克到數(shù)十噸的貨物精確空投,落點(diǎn)偏差甚至在幾十米以內(nèi)。翼傘空投系統(tǒng)已成為世界各軍事強(qiáng)國空投技術(shù)的發(fā)展方向,但是由于理論研究難度大且技術(shù)復(fù)雜,目前世界上掌握精確空投核心技術(shù)并有產(chǎn)品正式列裝的國家較少。

    翼傘系統(tǒng)飛行動力學(xué)是開展翼傘系統(tǒng)設(shè)計(jì)研制的理論基礎(chǔ)。文獻(xiàn)[2-3]建立了多種自由度的翼傘多體動力學(xué)模型,分析了展弦比、后緣下拉方式、安裝角及系統(tǒng)長度對雀降性能的影響。Goodrick[4]提出并建立了翼傘和回收物剛性連接的6自由度模型,研究了翼傘系統(tǒng)縱向運(yùn)動的穩(wěn)定性和飛行性能。Zhang等[5]在建立翼傘6自由度模型基礎(chǔ)上,通過實(shí)驗(yàn)驗(yàn)證了6自由度模型的可行性,用該模型仿真得到了翼傘各項(xiàng)參數(shù)變化對翼傘飛行狀態(tài)的影響。Slegers[6-7]建立了9自由度多體動力學(xué)模型,分析了安裝角、傘衣彎曲度及后緣下偏量對翼傘側(cè)向運(yùn)動的影響。蔣萬松等[8]采用Lagrange乘子法建立了兩體8自由度動力學(xué)仿真模型,并通過空投試驗(yàn)驗(yàn)證了仿真模型的有效性。Meirovitch[9]利用分布系統(tǒng)的擴(kuò)展哈密頓原理,導(dǎo)出混合系統(tǒng)的標(biāo)準(zhǔn)Lagrange方程,隨后以擬坐標(biāo)表示剛體運(yùn)動,將混合系統(tǒng)的2階Lagrange方程轉(zhuǎn)化為一組適用于控制的狀態(tài)方程。陳建平等[10]通過引入Lagrange乘子,獲得了大型翼傘系統(tǒng)封閉的飛行動力學(xué)方程組,對某大型翼傘系統(tǒng)飛行動力學(xué)過程進(jìn)行了數(shù)值仿真。楊海鵬等[11]假定傘繩和吊帶分別通過鉸鏈與翼傘和載荷相連,考慮了傘繩和吊帶的彈性以及傘繩相交點(diǎn)與傘衣的相對運(yùn)動,通過建模仿真得到了吊帶、傘繩交點(diǎn)運(yùn)動的情況。

    但是在目前文獻(xiàn)中,還很少有關(guān)于載荷連接模式和穩(wěn)定性的研究成果。為此,本文借鑒航天器多柔體動力學(xué)建模,以某動力的彈載翼傘系統(tǒng)為工程背景(見圖1),考慮傘繩、載荷吊帶的連接方式,針對翼傘穩(wěn)定飛行的運(yùn)動特性,將傘衣和載荷簡化為剛體,傘繩和吊帶交匯點(diǎn)處理為質(zhì)點(diǎn);分析翼傘系統(tǒng)各體的受力特性,給出了翼傘多項(xiàng)式形式的氣動表達(dá)式;將各體平動速度、角速度作為18自由度擬坐標(biāo),采用Lagrange擬速度形式的建模原理[12],建立了18自由度的動力學(xué)模型。利用空投試驗(yàn)證明了模型的有效性,仿真分析下拉后緣、突風(fēng)擾動和不同吊掛模式下的相對運(yùn)動特性,對比了6自由度單剛體和18自由度多體動力學(xué)模型的仿真結(jié)果。

    圖1 某小型翼彈載系統(tǒng)Fig.1 A small wing projectile-borne parafoil system

    1 模型描述

    如圖2所示,翼傘系統(tǒng)4體模型主要包括傘衣、傘繩、吊帶和載荷部分。本文將穩(wěn)定飛行狀態(tài)下的傘衣、載荷以及2個(gè)連接點(diǎn)作為4個(gè)剛體,傘衣和載荷分別記為剛體A和剛體B;傘繩和吊帶的交匯點(diǎn)簡化為質(zhì)點(diǎn),記為L和R;傘繩簡化為4根繩索LA1、RA2、RA3和LA4,吊帶簡化為4根繩索LB1、RB2、RB3和LB4. 上述多剛體模型中剛體A和B分別有6個(gè)自由度,質(zhì)點(diǎn)L和質(zhì)點(diǎn)R分別有3個(gè)自由度,總共18個(gè)自由度。該模型可以較完整地描述傘衣與載荷的相對運(yùn)動,反映了實(shí)際的空投工況和載荷吊掛模式。

    圖2 翼傘系統(tǒng)4體模型Fig.2 4-body model of parafoil system

    圖3 翼傘系統(tǒng)4體模型坐標(biāo)系示意圖Fig.3 Schematic diagram of 4-body model coordinate system of parafoil system

    圖4 側(cè)面示意圖Fig.4 Side view of parafoil system

    2 翼傘系統(tǒng)受力分析

    當(dāng)翼傘系統(tǒng)在其穩(wěn)定狀態(tài)飛行時(shí),主要受力包括重力、氣動力和螺旋槳推力等。傘衣具有柔性大的變形特性,而且傘繩和載荷的連接吊帶均為柔性繩索,很難獲得其翼傘系統(tǒng)精確的氣動特性。為此,本文不考慮翼傘下拉后緣過程中的遲滯效應(yīng),并忽略各部件的氣動耦合作用,采用工程簡化模型表示各體的氣動力。

    2.1 傘衣受力分析

    關(guān)于傘衣氣動模型的建立,劉媛等[13]通過有限體積元法求解Navier-Stokes湍流模型的控制方程,并對美國國家航空航天局蘭利研究中心的三維翼傘進(jìn)行了氣動性能分析,驗(yàn)證了該方法的可行性;進(jìn)而利用該方法分析得到了不同迎角來流條件以及傘衣后緣不同程度下偏的氣動力數(shù)據(jù)。汪龍芳等[14]對傘衣透氣性進(jìn)行了二維和三維數(shù)值模擬,得到了透氣性和迎角對翼傘升阻力系數(shù)變化的影響。為了便于后續(xù)歸航控制設(shè)計(jì),本文在計(jì)算流體力學(xué)基礎(chǔ)上,采用多項(xiàng)式表示傘衣質(zhì)心處的氣動力和氣動力矩[1],即

    (1)

    (2)

    式中:ρ為大氣密度;v為傘衣質(zhì)心相對于氣流的速度大??;SP為傘衣面積;δS和δa為對稱舵偏和反對稱舵偏;p、q和r分別為傘衣滾轉(zhuǎn)角速度、俯仰角速度和偏航角速度;b為傘衣展長;c為傘衣弦長;其他參數(shù)k、CD0、CL、CDδS、CDδa、CL0、CLα、CLδS、CLδa、CYβ、CYp、CYr、CYδa、Clβ、Clp、Clr、Cl、Clδa、Cm0、Cmα、Cmq、CmδS、Cmδa、Cnβ、Cnp、Cnr、Cnδa均為氣動力常系數(shù),由試驗(yàn)得到。

    由于傘衣充氣張滿后,其平均密度與空氣密度相當(dāng),必須考慮附加質(zhì)量的影響。本文基于文獻(xiàn)[15-16],將傘衣附加質(zhì)量等價(jià)地處理為傘衣質(zhì)心的作用力和力矩,并在傘衣坐標(biāo)系下表示為

    FPM=-ωP×(MavP),

    (3)

    MPM=-ωP×(JaωP),

    (4)

    式中:ωP和vP為傘衣角速度和質(zhì)心相對氣流速度;Ma和Ja分別為附加質(zhì)量相對傘衣質(zhì)心的質(zhì)量矩陣和慣量矩陣,即

    (5)

    (6)

    ma,ii和Ia,ii(i=1,2,3)分別為附加質(zhì)量相對傘衣質(zhì)心的平動和轉(zhuǎn)動分量。

    作用于剛體A質(zhì)心的合力,可在剛體A體坐標(biāo)系下表示為

    (7)

    作用于剛體A的合力矩可在剛體A體坐標(biāo)系下表示為

    (8)

    式中:MrtA和MrqA分別為各傘繩張力和氣動力對剛體A作用的力矩矢量和。

    2.2 載荷受力分析

    考慮到坐標(biāo)變換,作用于剛體B質(zhì)心處的合外力在其體坐標(biāo)系下表示為

    (9)

    式中:Ft為作用于載荷的螺旋槳推力;GB為剛體B所受重力;FrtB和FrqB分別為與載荷相連的各吊帶張力總和與氣動力矢量和;FqB為載荷所受氣動力,忽略螺旋槳對載荷氣動力的影響,則載荷所受氣動力為

    (10)

    式中:vB和Sd分別為載荷質(zhì)心速度和阻力面積。

    作用于剛體B的合力矩在體坐標(biāo)系下表示為

    (11)

    式中:MrtB和MrqB分別為與載荷相連的各吊帶張力和氣動力對載荷的作用力矩矢量和,在剛體B體坐標(biāo)系下分別表示為

    (12)

    (13)

    式中:rOBi為載荷體坐標(biāo)系下吊帶端點(diǎn)Bi的位置矢量,i=1,2,3,4;FrtLB1、FrtLB4、FrtRB2、FrtRB3分別為吊帶LB1、LB4、RB2、RB3的張力,方向從載荷接觸點(diǎn)指向吊帶交匯點(diǎn);FrqLB1、FrqLB4、FrqRB2、FrqRB3分別為吊帶LB1、LB4、RB2、RB3受到的氣動力在B1、B4、B2、B3上的集中力。

    2.3 傘繩和吊帶交匯點(diǎn)的受力分析

    作用于連接點(diǎn)L的合力可表示為

    FL=GL+FrtL+FrqL,

    (14)

    式中:GL為質(zhì)點(diǎn)L所受重力;FrtL和FrqL分別為與質(zhì)點(diǎn)L相連的各傘繩(吊帶)張力和氣動力。類似地,作用于連接點(diǎn)R的合力在慣性坐標(biāo)系下可表示為

    FR=GR+FrtR+FrqR,

    (15)

    式中:GR為質(zhì)點(diǎn)R所受重力;FrtR和FrqR分別為與質(zhì)點(diǎn)R相連的各傘繩(吊帶)張力和氣動力。另外,在具體仿真計(jì)算中將作用在繩索上的氣動力等效為在兩端點(diǎn)的集中力。

    3 基于擬坐標(biāo)形式的翼傘動力學(xué)方程

    (16)

    本文采用擬坐標(biāo)形式Lagrange方程進(jìn)行動力學(xué)建模,經(jīng)過復(fù)雜推導(dǎo)可得到如(17)式的動力學(xué)方程:

    (17)

    式中:相關(guān)狀態(tài)變換矩陣

    (18)

    (19)

    廣義質(zhì)量矩陣

    (20)

    I3×3為3階單位矩陣,03×3為3階0矩陣;廣義力的相關(guān)表達(dá)式為

    (21)

    (22)

    (23)

    (24)

    (25)

    (26)

    (27)

    4 仿真分析

    某翼傘系統(tǒng)的傘衣展長b=3.404 m,弦長c=0.920 m,傘衣和載荷安裝角均為0°,傘衣質(zhì)量mP為0.5 kg、載荷質(zhì)量mS為7.0 kg,初始高度和系統(tǒng)速度分別為1 000.0 m、v0=[7.5 0 3.8]Tm/s,其他狀態(tài)量的初始值均為0. 下面分別基于單剛體和多剛體模型進(jìn)行飛行特性分析。需要說明的是,本文基于文獻(xiàn)[3]中的研究成果,不考慮載荷的姿態(tài)變化和吊掛模式,將R點(diǎn)和L點(diǎn)的載荷等效為質(zhì)點(diǎn),建立了翼傘系統(tǒng)6自由度單剛體動力學(xué)模型。

    1)單剛體和多剛體動力學(xué)模型的對比分析。圖5和圖6所示分別為無控狀態(tài)下(無推力、無后緣下拉偏量)翼傘姿態(tài)和傘繩連接點(diǎn)位移的變化曲線。由圖5和圖6可以看出,因?yàn)槌跏紶顟B(tài)值為非穩(wěn)定狀態(tài),所以仿真初始時(shí)刻出現(xiàn)一定周期的震蕩才能達(dá)到穩(wěn)定飛行狀態(tài)。仿真結(jié)果表明,盡管兩種模型仿真的軌跡曲線幾乎一致、衰減趨勢基本一致,但是傘衣姿態(tài)角度稍有差別。在單剛體動力學(xué)模型中,由于將載荷簡化為質(zhì)點(diǎn)、忽略了載荷在物- 傘運(yùn)動過程中的相對運(yùn)動,其穩(wěn)定過程中主要依靠傘衣來消耗能量,因此傘衣擺動幅度相對結(jié)果略大。

    圖5 無控時(shí)傘衣質(zhì)心的位移變化曲線Fig.5 Changing curves of center of mass of canopy without control

    圖6 無控時(shí)傘衣姿態(tài)變化曲線Fig.6 Changing curves of attitude of canopy without control

    2)轉(zhuǎn)彎特性分析。圖7和圖8所示為基于18自由度多體模型、單側(cè)下偏量100%的載荷質(zhì)心軌跡和角速度變化曲線。從圖7和圖8中可以看出:在恒定下偏量作用下翼傘系統(tǒng)具有穩(wěn)定的轉(zhuǎn)彎半徑(本例最小轉(zhuǎn)彎半徑約52.0 m);同時(shí)在翼傘系統(tǒng)穩(wěn)定轉(zhuǎn)彎下滑過程中,載荷也具有穩(wěn)定的角速度。

    圖7 單側(cè)下拉時(shí)載荷質(zhì)心軌跡Fig.7 Centroid locus of payload when pulling down on the side

    圖8 轉(zhuǎn)彎時(shí)載荷角速度變化曲線Fig.8 Changing curves of angular velocity of payload during turning

    3) 帶動力飛行時(shí)的運(yùn)動特性分析。圖9和圖10所示為30 s時(shí)在載荷上施加50 N螺旋槳推力情況下載荷的軌跡和姿態(tài)角變化曲線。從圖9和圖10中可以看出,在一定推力作用下,翼傘系統(tǒng)可以實(shí)現(xiàn)平飛,但是傘衣和載荷都具有較大的俯仰角度。

    圖9 平飛時(shí)的載荷質(zhì)心位移Fig.9 Centroid displacement of payload during level flight

    圖10 平飛時(shí)的傘衣和載荷俯仰角Fig.10 Attitude angles of canopy and payload during level flight

    4)突風(fēng)作用的影響分析。圖11所示為在突風(fēng)2 m/s情況下,設(shè)計(jì)狀態(tài)(工況1)、傘繩長度增加20%(工況2)、吊帶增加30%(工況3)3種工況下載荷的角速度變化曲線。由圖12可以看出,在突風(fēng)影響下,適當(dāng)增加吊帶長度有利于載荷姿態(tài)的穩(wěn)定。

    圖11 突風(fēng)作用下的載荷姿態(tài)變化曲線Fig.11 Changing curves of payload attitude under wind gust

    5)空投試驗(yàn)分析。采用直升機(jī)進(jìn)行了空投試驗(yàn)(天氣良好),測得了某次軌跡曲線(單側(cè)下拉后緣50%)。圖12~圖14為載荷位置隨時(shí)間的變化曲線。由圖12~圖14可見,理論仿真數(shù)據(jù)與試驗(yàn)數(shù)據(jù)基本一致,證明了模型的有效性。二者的差異性主要是受高空風(fēng)場的作用。

    圖12 載荷xB軸方向位移隨時(shí)間變化曲線Fig.12 Payload displacement versus time in xB direction

    圖13 載荷yB軸方向位移隨時(shí)間變化曲線Fig.13 Payload displacement versus time in yB direction

    圖14 載荷zB軸方向位移隨時(shí)間變化曲線Fig.14 Payload displacement versus time in zB direction

    5 結(jié)論

    本文針對翼傘穩(wěn)定飛行時(shí)的運(yùn)動狀態(tài),綜合分析各體受力特點(diǎn)和運(yùn)動特性,利用擬坐標(biāo)形式的Lagrange原理建立了反映傘繩構(gòu)型和載荷吊帶連接方式的4體18自由度動力學(xué)模型。以某小型動力翼傘系統(tǒng)為例,對比研究單剛體和多體動力學(xué)模型的仿真特點(diǎn),分析了多體模型下單側(cè)下拉、突風(fēng)作用和各種吊掛模式下的傘衣和載荷運(yùn)動特性,以及外加推力后的傘衣姿態(tài)角和載荷位移變化情況,并與試驗(yàn)飛行數(shù)據(jù)做比較。主要得出如下結(jié)論:

    1)相比于6自由度模型,18自由度模型能較真實(shí)地反映翼傘系統(tǒng)內(nèi)部的相對運(yùn)動,更接近真實(shí)情況。

    2)在轉(zhuǎn)彎過程中,傘衣和載荷均有穩(wěn)定的角速度。

    3)載荷施加推力作用,傘衣將具有較大的俯仰角。

    4)適當(dāng)增加吊帶的長度,可增大系統(tǒng)受突風(fēng)影響的魯棒性。

    5)建立的模型仿真結(jié)果與試驗(yàn)數(shù)據(jù)比較,擬合度較好,證明了該模型的準(zhǔn)確性。

    猜你喜歡
    傘衣吊帶氣動力
    小型降落傘充氣和傘衣塌陷過程
    飛行載荷外部氣動力的二次規(guī)劃等效映射方法
    設(shè)計(jì)參數(shù)及大氣參數(shù)對降落傘充氣性能的影響
    盆底超聲檢測評估無張力性尿道中段懸吊帶術(shù)后吊帶相關(guān)并發(fā)癥的意義
    吊帶背心
    鳳凰生活(2019年10期)2019-10-14 03:28:35
    你的偶像
    夏日吊帶風(fēng)
    側(cè)風(fēng)對拍動翅氣動力的影響
    傘衣透氣性對翼傘氣動特性的影響
    高速鐵路接觸線覆冰后氣動力特性的風(fēng)洞試驗(yàn)研究
    男女啪啪激烈高潮av片| 90打野战视频偷拍视频| 女人被躁到高潮嗷嗷叫费观| 大片免费播放器 马上看| 男人操女人黄网站| 国产男女超爽视频在线观看| 国产高清三级在线| 精品国产国语对白av| 久久久久久伊人网av| 色网站视频免费| av线在线观看网站| av国产久精品久网站免费入址| 波多野结衣一区麻豆| 久久狼人影院| 国产精品99久久99久久久不卡 | 精品久久久久久电影网| 免费看av在线观看网站| 在线观看三级黄色| 久久99蜜桃精品久久| 国产精品麻豆人妻色哟哟久久| 一二三四中文在线观看免费高清| av免费观看日本| 亚洲经典国产精华液单| 十八禁高潮呻吟视频| 久久99热这里只频精品6学生| 成人国产麻豆网| 你懂的网址亚洲精品在线观看| 精品视频人人做人人爽| 欧美丝袜亚洲另类| 大片免费播放器 马上看| 2018国产大陆天天弄谢| 免费高清在线观看视频在线观看| 捣出白浆h1v1| 黑人巨大精品欧美一区二区蜜桃 | 免费少妇av软件| 天美传媒精品一区二区| 日韩av不卡免费在线播放| 国产精品人妻久久久久久| 纯流量卡能插随身wifi吗| 国内精品宾馆在线| 成人18禁高潮啪啪吃奶动态图| 久久久亚洲精品成人影院| 久久亚洲国产成人精品v| 国产无遮挡羞羞视频在线观看| 久久韩国三级中文字幕| 国产精品国产三级国产av玫瑰| 欧美精品一区二区大全| 人妻一区二区av| 极品少妇高潮喷水抽搐| 国产黄色视频一区二区在线观看| 日日爽夜夜爽网站| xxx大片免费视频| 久久99蜜桃精品久久| 插逼视频在线观看| 97超碰精品成人国产| 女性被躁到高潮视频| 免费观看无遮挡的男女| 亚洲精品国产av成人精品| 一级片'在线观看视频| 国产精品偷伦视频观看了| 夫妻性生交免费视频一级片| 亚洲欧美色中文字幕在线| 男女免费视频国产| 免费人成在线观看视频色| 久久99热6这里只有精品| 日日爽夜夜爽网站| 国产成人精品婷婷| 免费人妻精品一区二区三区视频| 成人免费观看视频高清| 考比视频在线观看| 高清在线视频一区二区三区| 久久精品夜色国产| 亚洲国产最新在线播放| 美女xxoo啪啪120秒动态图| 日韩三级伦理在线观看| 中文字幕人妻熟女乱码| 视频在线观看一区二区三区| 国产av码专区亚洲av| 久久久久久久久久成人| 青春草视频在线免费观看| 国产一级毛片在线| 美女国产高潮福利片在线看| 在线观看免费高清a一片| 午夜久久久在线观看| 国产一区二区激情短视频 | 18禁裸乳无遮挡动漫免费视频| 国产精品 国内视频| 尾随美女入室| 亚洲av在线观看美女高潮| 国产老妇伦熟女老妇高清| 伦理电影免费视频| av福利片在线| 中文天堂在线官网| 欧美激情 高清一区二区三区| 精品亚洲成国产av| 亚洲国产成人一精品久久久| 韩国av在线不卡| 最近手机中文字幕大全| 啦啦啦视频在线资源免费观看| 久久午夜综合久久蜜桃| av一本久久久久| 男女边吃奶边做爰视频| 一边亲一边摸免费视频| 全区人妻精品视频| 欧美激情 高清一区二区三区| tube8黄色片| 欧美成人午夜免费资源| 欧美成人午夜精品| 日韩av免费高清视频| 亚洲精品乱码久久久久久按摩| av.在线天堂| 夫妻性生交免费视频一级片| 午夜福利,免费看| 久久婷婷青草| 国产精品99久久99久久久不卡 | 性色avwww在线观看| 国产av国产精品国产| 中文精品一卡2卡3卡4更新| xxx大片免费视频| 日韩三级伦理在线观看| 亚洲国产av影院在线观看| 欧美日韩视频精品一区| 男女边摸边吃奶| 久久人妻熟女aⅴ| 日本欧美视频一区| 国产深夜福利视频在线观看| 亚洲五月色婷婷综合| 最近手机中文字幕大全| 制服人妻中文乱码| 中国三级夫妇交换| 日本av免费视频播放| 国产精品 国内视频| 肉色欧美久久久久久久蜜桃| 亚洲精华国产精华液的使用体验| 丁香六月天网| 三级国产精品片| 国内精品宾馆在线| 日日爽夜夜爽网站| 五月玫瑰六月丁香| 免费高清在线观看视频在线观看| 久久久久久久久久人人人人人人| 秋霞伦理黄片| 国产69精品久久久久777片| 我要看黄色一级片免费的| 日韩成人av中文字幕在线观看| 高清av免费在线| 欧美xxⅹ黑人| 亚洲av国产av综合av卡| 欧美精品一区二区免费开放| videosex国产| 久久久久久久久久人人人人人人| 熟女电影av网| 国产福利在线免费观看视频| 午夜影院在线不卡| 成年人午夜在线观看视频| 日韩免费高清中文字幕av| 精品少妇内射三级| 久久午夜综合久久蜜桃| 女人精品久久久久毛片| av免费观看日本| 亚洲一区二区三区欧美精品| 国产免费现黄频在线看| 高清在线视频一区二区三区| 蜜桃国产av成人99| 欧美成人午夜免费资源| 亚洲av在线观看美女高潮| 久久精品夜色国产| 午夜福利影视在线免费观看| 中国国产av一级| 2018国产大陆天天弄谢| 色吧在线观看| 香蕉国产在线看| 成人毛片60女人毛片免费| 看十八女毛片水多多多| 亚洲精华国产精华液的使用体验| 新久久久久国产一级毛片| 久久99热6这里只有精品| 免费日韩欧美在线观看| 国产成人av激情在线播放| 国产精品人妻久久久影院| 中文字幕人妻熟女乱码| 在线观看免费日韩欧美大片| 亚洲av日韩在线播放| 女人久久www免费人成看片| 国产精品不卡视频一区二区| 国产一区有黄有色的免费视频| 18禁动态无遮挡网站| 日韩不卡一区二区三区视频在线| 香蕉国产在线看| 一级,二级,三级黄色视频| 男女免费视频国产| 美女脱内裤让男人舔精品视频| 不卡视频在线观看欧美| 十八禁高潮呻吟视频| 在现免费观看毛片| 欧美少妇被猛烈插入视频| 亚洲精品色激情综合| 国产国语露脸激情在线看| 亚洲精品,欧美精品| 欧美人与善性xxx| 人妻一区二区av| 宅男免费午夜| 秋霞伦理黄片| 国产又爽黄色视频| 久久精品国产鲁丝片午夜精品| 色婷婷久久久亚洲欧美| 丰满少妇做爰视频| 国产激情久久老熟女| 色婷婷av一区二区三区视频| 婷婷色麻豆天堂久久| 美女脱内裤让男人舔精品视频| 久久久久国产网址| 亚洲性久久影院| 精品视频人人做人人爽| 一个人免费看片子| 国产无遮挡羞羞视频在线观看| 亚洲五月色婷婷综合| 久久久久国产精品人妻一区二区| 99久国产av精品国产电影| 成年人午夜在线观看视频| 丝袜美足系列| 黄色怎么调成土黄色| av天堂久久9| 有码 亚洲区| 国产精品免费大片| 久久这里只有精品19| 国产麻豆69| 国产一区二区在线观看日韩| 99香蕉大伊视频| 欧美bdsm另类| 一区二区三区精品91| freevideosex欧美| 久久免费观看电影| 全区人妻精品视频| 日韩av不卡免费在线播放| 纵有疾风起免费观看全集完整版| 亚洲精品视频女| 日本与韩国留学比较| 精品久久久久久电影网| 国产精品 国内视频| 成人手机av| 日韩不卡一区二区三区视频在线| 久久这里只有精品19| 99国产综合亚洲精品| 亚洲av电影在线进入| 亚洲成色77777| 人人妻人人澡人人爽人人夜夜| 国产成人一区二区在线| 你懂的网址亚洲精品在线观看| 久久人妻熟女aⅴ| 亚洲第一区二区三区不卡| 日本-黄色视频高清免费观看| 亚洲国产日韩一区二区| 久久鲁丝午夜福利片| 99国产精品免费福利视频| 亚洲欧美成人精品一区二区| 人人妻人人爽人人添夜夜欢视频| 伦精品一区二区三区| 亚洲精品美女久久av网站| av网站免费在线观看视频| av福利片在线| 高清av免费在线| 一区二区三区乱码不卡18| 久久人人爽人人爽人人片va| 国产日韩一区二区三区精品不卡| 春色校园在线视频观看| 亚洲欧美中文字幕日韩二区| 欧美亚洲日本最大视频资源| 丝袜喷水一区| 极品少妇高潮喷水抽搐| 国产精品人妻久久久久久| 波野结衣二区三区在线| 2022亚洲国产成人精品| 欧美精品一区二区免费开放| 精品一区二区免费观看| 搡老乐熟女国产| 丝袜美足系列| 卡戴珊不雅视频在线播放| 咕卡用的链子| 国产乱来视频区| 新久久久久国产一级毛片| 国产色婷婷99| 国产又色又爽无遮挡免| 国产亚洲av片在线观看秒播厂| 久久久国产欧美日韩av| 久久国产精品男人的天堂亚洲 | 国产一区二区三区综合在线观看 | 日韩伦理黄色片| 视频区图区小说| 看非洲黑人一级黄片| 超色免费av| 美女主播在线视频| av免费观看日本| 精品一区二区免费观看| 午夜激情久久久久久久| 欧美国产精品一级二级三级| 看免费av毛片| av国产久精品久网站免费入址| 欧美人与性动交α欧美软件 | 亚洲精华国产精华液的使用体验| 伦精品一区二区三区| av有码第一页| 精品人妻熟女毛片av久久网站| 国产在线免费精品| 国产不卡av网站在线观看| 99视频精品全部免费 在线| 最新的欧美精品一区二区| 国产精品 国内视频| 街头女战士在线观看网站| 少妇被粗大猛烈的视频| 成人国产av品久久久| 熟女av电影| 不卡视频在线观看欧美| 精品国产一区二区三区久久久樱花| 大话2 男鬼变身卡| 天堂中文最新版在线下载| 国产极品天堂在线| 搡女人真爽免费视频火全软件| 十分钟在线观看高清视频www| 色吧在线观看| 精品第一国产精品| 亚洲精品色激情综合| 黑人猛操日本美女一级片| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美日韩卡通动漫| 中文字幕av电影在线播放| 热99国产精品久久久久久7| 99热全是精品| 国产在线一区二区三区精| 久久久久久人妻| av视频免费观看在线观看| 国产精品一区二区在线观看99| 欧美xxⅹ黑人| 成年人午夜在线观看视频| 精品国产国语对白av| 欧美日韩亚洲高清精品| av视频免费观看在线观看| 你懂的网址亚洲精品在线观看| 如日韩欧美国产精品一区二区三区| av免费在线看不卡| 日韩中字成人| 男女边吃奶边做爰视频| 9191精品国产免费久久| 国产精品久久久久久精品电影小说| 婷婷色av中文字幕| 高清视频免费观看一区二区| 亚洲精品久久久久久婷婷小说| 国产亚洲午夜精品一区二区久久| 99久久精品国产国产毛片| 国产精品嫩草影院av在线观看| 亚洲伊人色综图| 国产成人精品福利久久| 久久精品人人爽人人爽视色| 久久99精品国语久久久| 五月伊人婷婷丁香| 亚洲精品色激情综合| 9色porny在线观看| 午夜91福利影院| 欧美亚洲日本最大视频资源| 人妻 亚洲 视频| 国产熟女欧美一区二区| 日产精品乱码卡一卡2卡三| 免费不卡的大黄色大毛片视频在线观看| 黄片播放在线免费| 国产精品99久久99久久久不卡 | av卡一久久| 亚洲欧美日韩另类电影网站| 一边摸一边做爽爽视频免费| tube8黄色片| 国产伦理片在线播放av一区| 日本黄大片高清| 欧美日韩综合久久久久久| 日韩伦理黄色片| 亚洲欧洲日产国产| 午夜福利视频精品| 日韩 亚洲 欧美在线| 亚洲国产精品国产精品| 尾随美女入室| 亚洲第一区二区三区不卡| 一区在线观看完整版| 一级毛片电影观看| 久久韩国三级中文字幕| 国产又爽黄色视频| 成人毛片a级毛片在线播放| 亚洲内射少妇av| 日韩成人av中文字幕在线观看| 性色avwww在线观看| 亚洲国产成人一精品久久久| 国产成人精品福利久久| 中国美白少妇内射xxxbb| 欧美精品人与动牲交sv欧美| 国产精品欧美亚洲77777| 欧美精品亚洲一区二区| 巨乳人妻的诱惑在线观看| 午夜福利网站1000一区二区三区| 街头女战士在线观看网站| 中国国产av一级| 亚洲欧美中文字幕日韩二区| 91成人精品电影| 人妻 亚洲 视频| 国产成人精品在线电影| 日本-黄色视频高清免费观看| 男人操女人黄网站| 色吧在线观看| 一区二区三区四区激情视频| 一区二区三区乱码不卡18| 国产精品蜜桃在线观看| 久久久久久久亚洲中文字幕| 久久久久国产网址| 久久女婷五月综合色啪小说| 午夜av观看不卡| 国产极品粉嫩免费观看在线| 成人18禁高潮啪啪吃奶动态图| 99热6这里只有精品| 久久精品人人爽人人爽视色| 中国三级夫妇交换| 下体分泌物呈黄色| 亚洲av在线观看美女高潮| 日韩一区二区视频免费看| 久久人妻熟女aⅴ| 亚洲激情五月婷婷啪啪| 国产亚洲欧美精品永久| 中国国产av一级| 国产老妇伦熟女老妇高清| 亚洲伊人色综图| 国产一区二区激情短视频 | 亚洲欧美色中文字幕在线| 国产毛片在线视频| 大话2 男鬼变身卡| 久久婷婷青草| av网站免费在线观看视频| √禁漫天堂资源中文www| 丰满少妇做爰视频| 韩国精品一区二区三区 | 成人手机av| 91精品伊人久久大香线蕉| 日日啪夜夜爽| 99视频精品全部免费 在线| 国产激情久久老熟女| www.熟女人妻精品国产 | 免费在线观看完整版高清| 黄片无遮挡物在线观看| 日本猛色少妇xxxxx猛交久久| 激情五月婷婷亚洲| 美女xxoo啪啪120秒动态图| 日韩 亚洲 欧美在线| 亚洲精品一二三| 亚洲美女黄色视频免费看| 精品少妇久久久久久888优播| 亚洲成人一二三区av| 捣出白浆h1v1| 国产国语露脸激情在线看| 国语对白做爰xxxⅹ性视频网站| 尾随美女入室| 国产一区二区在线观看av| 夜夜骑夜夜射夜夜干| 精品99又大又爽又粗少妇毛片| 国产在线视频一区二区| 观看av在线不卡| 在线观看免费高清a一片| 丁香六月天网| 亚洲国产精品一区二区三区在线| 九色成人免费人妻av| 午夜激情av网站| 少妇被粗大的猛进出69影院 | 久久久欧美国产精品| 99视频精品全部免费 在线| 欧美日韩精品成人综合77777| 啦啦啦中文免费视频观看日本| 精品亚洲成a人片在线观看| 母亲3免费完整高清在线观看 | 亚洲欧美精品自产自拍| 国产成人av激情在线播放| 丝袜喷水一区| 国产成人精品久久久久久| 欧美另类一区| 黄色 视频免费看| 国产亚洲最大av| 熟女人妻精品中文字幕| 伊人久久国产一区二区| av.在线天堂| 日本wwww免费看| 中文字幕精品免费在线观看视频 | 看十八女毛片水多多多| 色哟哟·www| 国产精品久久久久久精品电影小说| 日韩在线高清观看一区二区三区| 国产极品粉嫩免费观看在线| 看免费av毛片| 国产黄色视频一区二区在线观看| 欧美激情 高清一区二区三区| 成人手机av| 久久午夜福利片| 久久精品熟女亚洲av麻豆精品| 久久久亚洲精品成人影院| 亚洲情色 制服丝袜| 99热全是精品| 日韩熟女老妇一区二区性免费视频| 最黄视频免费看| 欧美人与性动交α欧美软件 | 亚洲精品久久久久久婷婷小说| av国产久精品久网站免费入址| 宅男免费午夜| 亚洲精品色激情综合| 美女视频免费永久观看网站| 欧美激情 高清一区二区三区| 亚洲精品视频女| 91在线精品国自产拍蜜月| 老女人水多毛片| 男女边摸边吃奶| 日日摸夜夜添夜夜爱| 亚洲图色成人| a级片在线免费高清观看视频| 亚洲美女黄色视频免费看| 亚洲,一卡二卡三卡| 精品99又大又爽又粗少妇毛片| 建设人人有责人人尽责人人享有的| 一级a做视频免费观看| 欧美丝袜亚洲另类| 99re6热这里在线精品视频| 久久99一区二区三区| 啦啦啦视频在线资源免费观看| 97精品久久久久久久久久精品| 黑丝袜美女国产一区| 2022亚洲国产成人精品| freevideosex欧美| 国产av国产精品国产| 午夜福利视频在线观看免费| 日韩,欧美,国产一区二区三区| 又粗又硬又长又爽又黄的视频| 婷婷色综合www| 国产不卡av网站在线观看| 日韩精品有码人妻一区| 看免费av毛片| 一个人免费看片子| 亚洲精品久久久久久婷婷小说| 久久韩国三级中文字幕| 国产淫语在线视频| 波多野结衣一区麻豆| 国产亚洲精品久久久com| 国产又色又爽无遮挡免| 哪个播放器可以免费观看大片| 人妻人人澡人人爽人人| 久热这里只有精品99| 男人操女人黄网站| 国产精品熟女久久久久浪| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产老妇伦熟女老妇高清| 黄网站色视频无遮挡免费观看| 亚洲三级黄色毛片| 亚洲精品久久久久久婷婷小说| 亚洲欧洲精品一区二区精品久久久 | 亚洲在久久综合| av国产精品久久久久影院| 午夜福利在线观看免费完整高清在| 国产精品一二三区在线看| 少妇人妻 视频| 国产亚洲av片在线观看秒播厂| 国产亚洲最大av| 欧美日本中文国产一区发布| 亚洲 欧美一区二区三区| 国产成人精品久久久久久| 国产无遮挡羞羞视频在线观看| 亚洲伊人久久精品综合| 免费看不卡的av| 好男人视频免费观看在线| 国产免费一级a男人的天堂| 飞空精品影院首页| 日韩中文字幕视频在线看片| 99国产精品免费福利视频| 亚洲av免费高清在线观看| 超色免费av| tube8黄色片| www.熟女人妻精品国产 | 老司机亚洲免费影院| 一区在线观看完整版| 久久精品熟女亚洲av麻豆精品| 一区二区三区四区激情视频| 精品福利永久在线观看| 天天躁夜夜躁狠狠躁躁| 国产视频首页在线观看| 精品国产一区二区三区久久久樱花| 激情视频va一区二区三区| 亚洲国产av新网站| 一区二区日韩欧美中文字幕 | 国产精品人妻久久久久久| av电影中文网址| 欧美日韩视频精品一区| 亚洲三级黄色毛片| 丝袜美足系列| 18在线观看网站| tube8黄色片| 中文字幕av电影在线播放| 日韩人妻精品一区2区三区| 精品人妻一区二区三区麻豆| 在线观看美女被高潮喷水网站| 制服丝袜香蕉在线| 自线自在国产av| 久久午夜综合久久蜜桃| 狠狠婷婷综合久久久久久88av| 晚上一个人看的免费电影| 国产亚洲午夜精品一区二区久久| 久久精品人人爽人人爽视色| 少妇 在线观看| av片东京热男人的天堂| 久久精品国产亚洲av涩爱| av女优亚洲男人天堂| 亚洲国产av新网站| 欧美激情国产日韩精品一区| 精品一区二区三区视频在线| 日本wwww免费看| 尾随美女入室| 久久久久精品性色| 精品久久久久久电影网| 色5月婷婷丁香|