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

    小天體著陸動(dòng)力學(xué)參數(shù)不確定性影響分析

    2014-03-06 05:08:45袁旭朱圣英喬棟崔平遠(yuǎn)
    深空探測學(xué)報(bào) 2014年2期

    袁旭,朱圣英,喬棟,崔平遠(yuǎn)

    (1.北京理工大學(xué)小天體探測與防御實(shí)驗(yàn)室,北京100081;2.飛行器動(dòng)力學(xué)與控制教育部重點(diǎn)實(shí)驗(yàn)室,北京100081)

    小天體著陸動(dòng)力學(xué)參數(shù)不確定性影響分析

    袁旭1,2,朱圣英1,2,喬棟1,2,崔平遠(yuǎn)1,2

    (1.北京理工大學(xué)小天體探測與防御實(shí)驗(yàn)室,北京100081;2.飛行器動(dòng)力學(xué)與控制教育部重點(diǎn)實(shí)驗(yàn)室,北京100081)

    針對小天體不規(guī)則程度高、引力場復(fù)雜,且物理參數(shù)存在較高不確定性的問題,基于小天體著陸動(dòng)力學(xué)方程線性化近似解析解,對各動(dòng)力學(xué)參數(shù)不確定性的影響進(jìn)行了分析??紤]動(dòng)力學(xué)方程線性化帶來的誤差,引入線性化誤差補(bǔ)償校正方法,建立了探測器軌跡對動(dòng)力學(xué)參數(shù)不確定性的敏感度方程。以小行星Eros 433為例,重點(diǎn)分析了目標(biāo)小天體質(zhì)量、自轉(zhuǎn)角速度、引力勢函數(shù)系數(shù),以及探測器初始狀態(tài)、推力加速度等動(dòng)力學(xué)參數(shù)不確定性對探測器著陸軌跡的影響。數(shù)學(xué)仿真分析表明,針對本文選取的目標(biāo)小天體,推力加速度擾動(dòng)為主要影響因素,探測器初始狀態(tài)的不確定性為次要影響因素,其他參數(shù)擾動(dòng)的影響較小。

    小天體;著陸;參數(shù)不確定性;敏感性分析

    0 引言

    對小天體的探測研究,有利于人們認(rèn)識和研究太陽系的起源與演化。因此,小天體探測是21世紀(jì)深空探測活動(dòng)的重要內(nèi)容[12]。其中,著陸探測由于能夠進(jìn)行現(xiàn)場科學(xué)探測,因而引起學(xué)者的廣泛興趣[34]。迄今,各國的小天體探測任務(wù)中,有2次成功完成了著陸任務(wù),分別為美國航空航天局(NASA)發(fā)射的NEAR探測器于2001年2月著陸于Eros 433小行星[5],和日本宇宙航空研究開發(fā)機(jī)構(gòu)(JAXA)發(fā)射的Hayabusa探測器于2005年11月著陸于Itokawa小行星并采樣返回[6]。此外,歐洲太空局(ESA)發(fā)射的Rosetta探測器預(yù)計(jì)將于2014年11月對Churyumov-Gerasimenko彗星進(jìn)行著陸探測[7]。

    形狀的不規(guī)則性是小天體顯著的物理特征。由于形狀極度不規(guī)則,小天體引力場環(huán)境較為復(fù)雜,小天體的各物理參數(shù)也存在較大的不確定性[89]。這種參數(shù)不確定性會(huì)對小天體著陸過程產(chǎn)生影響,給小天體著陸任務(wù)帶來挑戰(zhàn)。同時(shí),發(fā)動(dòng)機(jī)推力誤差等探測器自有誤差也對小天體著陸精度造成影響。因此,有必要對小天體著陸過程中參數(shù)不確定性的影響進(jìn)行分析與研究,通過分析各參數(shù)不確定性對探測器狀態(tài)的影響程度,可重點(diǎn)關(guān)注影響較大的參數(shù)并設(shè)法減小相應(yīng)的誤差。Broschart等對小天體固聯(lián)坐標(biāo)系下探測器的平移與下降運(yùn)動(dòng)特點(diǎn)進(jìn)行了大量研究,并進(jìn)行了敏感性分析,取得了顯著的研究成果[10]。

    本文針對小天體下降與著陸過程,首先給出動(dòng)力學(xué)方程并對其線性化,利用線性化方程對探測器狀態(tài)及其他變量求解。基于這些方程和敏感性分析方法,研究了探測器狀態(tài)對參數(shù)不確定性的敏感度,對相同條件下處于不同區(qū)域時(shí)探測器狀態(tài)對各參數(shù)的敏感度進(jìn)行了分析與比較,并以Eros 433小行星參數(shù)[9]為例給出了仿真分析結(jié)果。

    1 系統(tǒng)動(dòng)力學(xué)模型

    1.1 坐標(biāo)系

    本文采用小天體固聯(lián)坐標(biāo)系Σb:原點(diǎn)ob位于小天體質(zhì)心,zb軸與小天體最大慣量軸即自轉(zhuǎn)軸重合,xb與yb軸分別與最小和中間慣量軸重合,xb, yb,zb三軸滿足右手法則,如圖1所示。

    圖1 小天體固聯(lián)坐標(biāo)系Fig.1 Small celestial body fixed frame

    1.2 動(dòng)力學(xué)方程

    假設(shè)小天體密度均勻并繞其最大慣量主軸勻速轉(zhuǎn)動(dòng),忽略其他干擾力,在小天體固聯(lián)坐標(biāo)系下,選取探測器位置、速度為狀態(tài)變量,即

    則探測器的動(dòng)力學(xué)方程可表示為

    其中,ω為小天體自轉(zhuǎn)角速度;Tx,Ty,Tz為探測器控制力形成的三軸加速度分量;V為小天體引力勢函數(shù);為三軸的小天體引力加速度分量。

    小天體引力勢函數(shù)采用球諧函數(shù)形式,即

    其中,Pnm為締結(jié)勒讓德多項(xiàng)式;Cnm和Snm為小天體引力場對應(yīng)的球諧函數(shù)各階系數(shù);a為小天體名義半徑;θ,φ,r分別為探測器所處的經(jīng)度、緯度和半徑,滿足

    1.3 動(dòng)力學(xué)方程線性化

    由于探測器受到的小天體引力加速度與探測器位置有關(guān),因此式(1)所示的動(dòng)力學(xué)方程為非線性方程,且無法求出精確解析解。在一些情況下,需要解析形式的近似解。此時(shí)可將式(1)進(jìn)行線性化處理,得到近似的線性化方程

    其中

    則可求出此線性方程的解析解

    雖然探測器的運(yùn)動(dòng)存在很大的非線性,但線性化假設(shè)可在一定精度下簡化分析過程,并能獲得解析形式的近似解,便于進(jìn)一步的分析,在初期的理論推導(dǎo)和任務(wù)設(shè)計(jì)階段是非常方便和重要的分析手段。

    由于矩陣A+δ可逆,若在研究的時(shí)間范圍內(nèi),推力加速度T為常矢量,則此解可進(jìn)一步化簡為

    式(5)即為常推力加速度下動(dòng)力學(xué)方程式(3)的線性近似解析解。

    2 線性化誤差校正

    2.1 利用線性化方程求解未知變量

    式(5)中含有6個(gè)方程,包括始末狀態(tài)、推力加速度、作用時(shí)間共計(jì)16個(gè)變量。若已知其中10個(gè)變量,則可求解出6個(gè)未知變量。

    例如,若已知初始狀態(tài)X(0),欲在常推力加速度T作用下經(jīng)時(shí)間tf到達(dá)末端位置rf,則可求得T和到達(dá)末端位置時(shí)的速度vf分別為:

    2.2 線性化誤差及誤差校正

    線性化方程求解出的探測器狀態(tài)為近似解,與真實(shí)狀態(tài)存在誤差。探測器的真實(shí)狀態(tài)X*(t)可表示為

    則線性化誤差可表示為

    假設(shè)在適當(dāng)時(shí)間范圍內(nèi),線性化解(5)的軌跡接近真實(shí)軌跡,且此區(qū)域內(nèi)也很接近,則可認(rèn)為在此區(qū)域內(nèi)進(jìn)行線性化帶來的誤差總是近似相等的。這樣,可對式(6)所求得的推力加速度T進(jìn)行校正[10]:按照非線性方程式(1)進(jìn)行數(shù)值求解后,可求出末時(shí)刻的線性化誤差E(tf)?;赬*(tf)附近不同軌跡對應(yīng)的E(tf)近似相等的假設(shè),可設(shè)一個(gè)用于校正的,滿足

    圖2 線性化誤差修正示意圖Fig.2 Scheme of linearization error calibration

    圖3 線性化誤差修正局部放大圖Fig.3 Calibration scheme partially amplified

    3 探測器狀態(tài)對參數(shù)不確定性的敏感度

    從動(dòng)力學(xué)方程的解式(5)及引力勢函數(shù)式(2)可以看出,探測器t時(shí)刻的狀態(tài)變量X(t)與其初始狀態(tài)、發(fā)動(dòng)機(jī)推力加速度、小天體自轉(zhuǎn)角速度、質(zhì)量和球諧函數(shù)各階系數(shù)有關(guān),即

    于是,式(10)中各參數(shù)的誤差或不確定性都會(huì)對探測器的狀態(tài)與著陸軌跡產(chǎn)生影響。

    在tf時(shí)刻,以αi,i=1,2,...,k表示參數(shù)ω, M,Cnm,Snm,X(0),T,則X(tf)的協(xié)方差矩陣可表示為:

    由式(11)可以看出,X(tf)的不確定度與各參數(shù)的不確定度和X(tf)對該參數(shù)的敏感度矩陣有關(guān),即由系統(tǒng)狀態(tài)方程和各參數(shù)不確定度共同決定。探測器狀態(tài)對某參數(shù)的敏感度矩陣可由以下過程求得[11]

    結(jié)合動(dòng)力學(xué)方程

    對某一參數(shù)αi,其敏感度矩陣滿足

    通過解上式中的微分方程即可求出對某參數(shù)的敏感度矩陣,再結(jié)合該參數(shù)的不確定度,即可得到某時(shí)刻探測器狀態(tài)X(tf)的不確定度。

    由式(11)可見,已知各參數(shù)不確定的情況下,探測器狀態(tài)對各參數(shù)的敏感度矩陣隨時(shí)間變化,因此由某一參數(shù)擾動(dòng)引起的探測器狀態(tài)不確定度也隨時(shí)間變化。可研究相同條件下,同一時(shí)刻各參數(shù)不確定性對探測器狀態(tài)的影響,來比較某一任務(wù)中不同參數(shù)擾動(dòng)的作用大小;對某一參數(shù),可研究探測器處于不同區(qū)域(初始狀態(tài))時(shí),其狀態(tài)受該參數(shù)擾動(dòng)影響的差異。這里,取探測器位置向量協(xié)方差陣最大特征值的平方根作為探測器位置不確定度的度量指標(biāo)

    4 仿真與分析

    4.1 線性化誤差校正的仿真分析

    以Eros 433小行星為例進(jìn)行仿真分析。假設(shè)探測器在常推力作用下沿徑向向小行星表面下降(即保持經(jīng)緯度不變)。仿真中選取探測器初始位置半徑為20 km,經(jīng)緯度選取多組數(shù)據(jù)以比較,初始速度大小為3.0 m/s,方向指向小天體質(zhì)心即坐標(biāo)系原點(diǎn),末端位置半徑為17.8 km,運(yùn)動(dòng)時(shí)間為1 000 s。在某一初始位置,可由式(6)計(jì)算所需推力加速度并通過式(9)校正。仿真末端位置與目標(biāo)末端位置的偏差(距離)隨初始位置的經(jīng)緯度不同而變化,可繪制等高線圖表示。仿真所用參數(shù)參照NEAR任務(wù)獲得的數(shù)據(jù),如表1所示,未經(jīng)誤差校正和校正后末端位置偏差的區(qū)域分布如圖4和圖5所示。

    表1 仿真參數(shù)Table 1 Simulation parameters

    圖4 未校正時(shí)末端位置與目標(biāo)位置偏差等高線圖(單位:m)Fig.4 Contour map of terminal position offset from target without calibration(unit:m)

    由圖4和圖5可見,在此仿真條件下:

    1)校正前的末端位置偏差一般處于米級,校正后減小至厘米級,精度提高了2個(gè)數(shù)量級,表明此校正方法非常有效。

    圖5 校正后末端位置與目標(biāo)位置偏差等高線圖(單位:cm)Fig.5 Contour map of terminal position offset from target after calibration(unit:cm)

    2)緯度愈靠近0°,經(jīng)度愈靠近0°或±180°,即愈靠近xb軸,則偏差愈大。這表明愈靠近長軸,小天體形狀的不規(guī)則度愈大,引力加速度的非線性愈強(qiáng),線性化引入的誤差愈大。

    4.2 探測器狀態(tài)對參數(shù)敏感度的仿真分析

    仿真參數(shù)仍如表1所示,各參數(shù)的1—σ不確定度參照NEAR任務(wù)中獲得的數(shù)據(jù)以及文獻(xiàn)[10],如表2所示。

    表2 仿真參數(shù)的1-σ不確定度Table 2 1-σuncertainty of the simulation parameters

    其中,對球諧函數(shù)系數(shù)的不確定度考察到3階, Tmag,?T,θT分別代表推力加速度大小和推力方向角。

    探測器末端位置不確定度以式(14)中σ表達(dá),得到圖6所示的等高線圖。由各參數(shù)擾動(dòng)引起的不確定度數(shù)量級如表3所示,其中推力加速度擾動(dòng)引起的末端位置不確定度分布如圖7所示。

    在此仿真條件下:

    1)將圖6與圖7進(jìn)行比較發(fā)現(xiàn),二者等高線圖非常相似,數(shù)值也非常接近,表明末端位置的不確定度主要源于推力加速度的不確定性,即末端位置不確定度受發(fā)動(dòng)機(jī)推力不確定性影響最大,發(fā)動(dòng)機(jī)推力偏差是影響探測器著陸軌跡不確定度的主要因素。

    表3 各參數(shù)擾動(dòng)引起的末端位置不確定度的數(shù)量級Table 3 Order of magnitude of terminal position uncertainty due to uncertainty of each parameter alone

    圖6 末端位置不確定度等高線圖(單位:m)Fig.6 Contour map of terminal position uncertainty(unit:m)

    圖7 參數(shù)T擾動(dòng)引起的末端位置不確定度等高線圖(單位:m)Fig.7 Contour map of terminal position uncertainty due to uncertainty of T(unit:m)

    2)末端位置的不確定度為10 m級,主要分布于30~50 m范圍。推力加速度偏差引起的末端位置不確定度為10 m級,主要分布于25~45 m范圍,為主要影響因素;初始狀態(tài)偏差引起的末端位置不確定度為米級,主要分布于5~6 m范圍,為次要影響因素;小天體質(zhì)量、引力勢函數(shù)系數(shù)、自轉(zhuǎn)角速度不確定性引起的末端位置不確定度均小于米級??梢娫诖朔抡鏃l件下,探測器軌跡的不確定度主要受推力偏差影響,其次受初始狀態(tài)誤差影響,受其他參數(shù)擾動(dòng)的影響較小。

    5 結(jié)論

    本文對小天體著陸過程中探測器狀態(tài)對小天體及探測器參數(shù)的敏感性進(jìn)行了簡要的分析。首先,建立非線性動(dòng)力學(xué)方程并進(jìn)行了線性化,得出線性化方程的近似解析解?;诰€性化方程,得出常推力下求解未知變量的方法,并引入一種線性化誤差校正方法。探測器狀態(tài)的不確定度由各參數(shù)不確定度和探測器狀態(tài)對各參數(shù)的敏感度共同決定。通過敏感性分析,比較了探測器著陸軌跡在相同條件下受各參數(shù)擾動(dòng)影響的區(qū)別和處于不同區(qū)域時(shí)受某一參數(shù)擾動(dòng)影響的差異。在本文的仿真條件下,推力加速度的擾動(dòng)為主要影響因素,初始狀態(tài)的不確定性為次要影響因素,其他參數(shù)擾動(dòng)的影響較小。

    [1] 崔平遠(yuǎn),喬棟.小天體附近軌道動(dòng)力學(xué)與控制研究現(xiàn)狀與展望[J].力學(xué)進(jìn)展,2013,43(5):526- 539.[Cui P Y,Qiao D. Research progress and prospect of orbital dynamics and control near small bodies[J].Advances in Mechanics,2013, 43(5):526- 539.]

    [2] Furfaro R,Cersosimo D,Wibben D R.Asteroid precision landing via multiple sliding surfaces guidance techniques[J]. Journal of Guidance,Control,and Dynamics,2013:1-18.

    [3] Barucci M,Cheng A F,Michel P,et al.MarcoPolo-R nearearth asteroid sample return mission[J].Experimental Astronomy,2012,33(2- 3):645- 684.

    [4] Michel P,Barucci M A,Cheng A F,et al.MarcoPolo-R: near-earth asteroid sample return mission selected for the assessment study phase of the ESA program cosmic vision [J].Acta Astronautica,2014,93(0):530- 538.

    [5] Dunham D W,Farquhar R W,Mcadams J V,et al. Implementation of the first asteroid landing[J].Icarus, 2002,159(2):433- 438.

    [6] Uo M,Shirakawa K,Hashimoto T,et al.Hayabusa's touching-down to Itokawa-autonomous guidance and navigation[J].The Journal of Space Technology and Science, 2006,22(1):41.

    [7] Ulamec S,Biele J,Bousquet P,et al.Landing on small bodies:from the Rosetta lander to{MASCOT}and beyond [J].Acta Astronautica,2014,93(0):460- 466.

    [8] 崔平遠(yuǎn),朱圣英,崔祜濤.小天體自主軟著陸脈沖機(jī)動(dòng)控制方法研究[J].宇航學(xué)報(bào),2008,29(2):511- 516.[Cui P Y,Zhu S H,Cui H T.Autonomous impulse maneuver control method for soft landing on small bodies[J].Journal of Astronautics,2008,29(2):511- 516.]

    [9] Miller J K,Konopliv A S,Antreasian P G,et al. Determination of shape,gravity,and rotational state of asteroid 433 Eros[J].Icarus,2002,155(1):3- 17.

    [10] Broschart S B,Scheeres D J.Spacecraft descent and translation in the small-body fixed frame[C]//Proc.AIAA/ AAS Astrodynamics Specialists Meeting.Providence,RI, USA:[s.n.],2004.

    [11] Montenbruck O,Gill E.Satellite orbits:models,methods, and applications[M].Berlin Heidelberg:Springer,2000.

    電話:(010)68918910

    E-mail:yuanxu1020@126.com

    [責(zé)任編輯:宋宏]

    Impact Analysis of Dynamic Parameters Uncertainty on Small Celestial Body Landing

    YUAN Xu1,2,ZHU Shengying1,2,QIAO Dong1,2,CUI Pingyuan1,2
    (1.Laboratory of Asteroid Exploration and Defense,Beijing Institute of Technology,Beijing 100081,China;2.Key Laboratory of Dynamics and Control of Flight Vehicle,Ministry of Education,Beijing 100081,China)

    To the challenge of high degree irregularity,complex gravity field and relatively high degree of physical parameters uncertainty of small celestial bodies,the impact of dynamic parameters uncertainty is analyzed, based on the approximate linearized analytical solution of the small celestial body landing dynamic equations.A calibrating method is used to compensate the errors brought in by linearization.And the sensitivity equations of spacecraft trajectories to dynamic parameters are established.Taking the asteroid Eros 433 as an example,the impact of the uncertainty of the mass,rotation rate and gravitational potential coefficients of the target body on the spacecraft landing trajectory,the initial state of the spacecraft and the thrust acceleration are analyzed elaborately. Mathematical simulation shows that for this selected target body,the thrust acceleration disturbance is the primary factor of such impact,and the initial state of the spacecraft the secondary one.And the impact of other parameters are relatively small.

    small celestial body;landing;parameters uncertainty;sensitivity analysis

    V448.21

    :A

    :2095-7777(2014)02-0134-06

    袁旭(1986—),男,博士研究生,主要研究方向:深空探測器自主導(dǎo)航與控制

    2013-11-01

    2013-12-30

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)(2012CB720000);國家自然科學(xué)基金(61304248);基礎(chǔ)科研計(jì)劃(B2220110012);民用航天預(yù)先研究計(jì)劃

    国产伦在线观看视频一区| 校园人妻丝袜中文字幕| 久久亚洲国产成人精品v| 99久国产av精品| 国产精品久久久久久av不卡| 五月玫瑰六月丁香| 亚洲av中文av极速乱| 黄色一级大片看看| 亚洲第一电影网av| 男人的好看免费观看在线视频| 亚洲精品国产成人久久av| 黄片wwwwww| 九九在线视频观看精品| 久久这里有精品视频免费| 国产91av在线免费观看| 1024手机看黄色片| 久久久精品欧美日韩精品| 国产精品国产三级国产av玫瑰| 黄片wwwwww| 两个人的视频大全免费| 久久久欧美国产精品| 男女那种视频在线观看| 麻豆精品久久久久久蜜桃| 免费观看的影片在线观看| 偷拍熟女少妇极品色| 有码 亚洲区| 我的女老师完整版在线观看| av在线蜜桃| 国产伦一二天堂av在线观看| 天堂影院成人在线观看| 非洲黑人性xxxx精品又粗又长| 天天躁日日操中文字幕| 国国产精品蜜臀av免费| 搡女人真爽免费视频火全软件| 人妻夜夜爽99麻豆av| 黄色日韩在线| 成人午夜精彩视频在线观看| 男女那种视频在线观看| 国产av一区在线观看免费| 日本撒尿小便嘘嘘汇集6| 国产片特级美女逼逼视频| av在线播放精品| 综合色av麻豆| 成年女人永久免费观看视频| 色哟哟·www| 99热只有精品国产| 亚洲成人中文字幕在线播放| 欧美又色又爽又黄视频| 99久久久亚洲精品蜜臀av| 白带黄色成豆腐渣| 看黄色毛片网站| 欧美潮喷喷水| 国产v大片淫在线免费观看| 中出人妻视频一区二区| 精品人妻熟女av久视频| 精品久久久久久久末码| 噜噜噜噜噜久久久久久91| 黑人高潮一二区| 中文欧美无线码| 白带黄色成豆腐渣| 99久久精品一区二区三区| 亚洲人成网站在线播| 51国产日韩欧美| 久久九九热精品免费| 18禁在线无遮挡免费观看视频| 免费观看人在逋| 国产一区二区亚洲精品在线观看| 亚洲国产欧美在线一区| 免费不卡的大黄色大毛片视频在线观看 | 国产av麻豆久久久久久久| 亚洲一区高清亚洲精品| 麻豆成人午夜福利视频| 美女cb高潮喷水在线观看| 女人十人毛片免费观看3o分钟| 激情 狠狠 欧美| 舔av片在线| a级毛片免费高清观看在线播放| 我的女老师完整版在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产午夜精品久久久久久一区二区三区| 悠悠久久av| 天堂中文最新版在线下载 | 51国产日韩欧美| 春色校园在线视频观看| 日韩av在线大香蕉| 免费看美女性在线毛片视频| 又粗又硬又长又爽又黄的视频 | 久久精品综合一区二区三区| 国产亚洲91精品色在线| 久久亚洲精品不卡| 少妇被粗大猛烈的视频| 欧美日韩精品成人综合77777| 久久欧美精品欧美久久欧美| 久久精品91蜜桃| 看免费成人av毛片| 尾随美女入室| 亚洲三级黄色毛片| 日本-黄色视频高清免费观看| 亚洲,欧美,日韩| 欧美最黄视频在线播放免费| 日本三级黄在线观看| 女人被狂操c到高潮| 韩国av在线不卡| av女优亚洲男人天堂| 久久午夜亚洲精品久久| 18禁在线无遮挡免费观看视频| 欧美zozozo另类| 岛国在线免费视频观看| 激情 狠狠 欧美| 欧美一区二区国产精品久久精品| АⅤ资源中文在线天堂| 99热网站在线观看| 国产蜜桃级精品一区二区三区| 久久久久久久久久黄片| 黑人高潮一二区| 国产精品人妻久久久影院| 国产伦精品一区二区三区视频9| av免费观看日本| av国产免费在线观看| 欧美bdsm另类| 精华霜和精华液先用哪个| 久久国产乱子免费精品| 男人舔奶头视频| 毛片女人毛片| 18禁在线播放成人免费| 亚洲国产欧洲综合997久久,| 97超视频在线观看视频| 亚洲国产欧美在线一区| av又黄又爽大尺度在线免费看 | 一级毛片aaaaaa免费看小| 欧美最黄视频在线播放免费| 日韩国内少妇激情av| 国产午夜精品论理片| 国产成人精品一,二区 | 少妇丰满av| 亚洲成人中文字幕在线播放| 五月伊人婷婷丁香| 国产一区二区在线观看日韩| 国产蜜桃级精品一区二区三区| 亚洲av男天堂| 搡女人真爽免费视频火全软件| 国产私拍福利视频在线观看| 亚洲av男天堂| 国产高清视频在线观看网站| 九色成人免费人妻av| 久久婷婷人人爽人人干人人爱| 精品一区二区免费观看| 亚洲一级一片aⅴ在线观看| 免费电影在线观看免费观看| 亚州av有码| 一级毛片电影观看 | 波野结衣二区三区在线| 欧美不卡视频在线免费观看| 小说图片视频综合网站| 麻豆久久精品国产亚洲av| 成人毛片60女人毛片免费| 日本黄大片高清| 亚洲精品自拍成人| av天堂在线播放| 一区二区三区高清视频在线| 亚洲欧美精品自产自拍| 蜜臀久久99精品久久宅男| 2022亚洲国产成人精品| 久久久午夜欧美精品| 亚洲真实伦在线观看| 精品人妻视频免费看| 亚洲精品成人久久久久久| 亚洲欧美精品专区久久| 午夜a级毛片| 久久国产乱子免费精品| 精品久久久久久久久亚洲| eeuss影院久久| eeuss影院久久| 久久亚洲精品不卡| 国产成人福利小说| 一区二区三区四区激情视频 | 久久精品国产亚洲av涩爱 | videossex国产| 美女 人体艺术 gogo| 国产在视频线在精品| 国产精品国产高清国产av| 不卡一级毛片| 久久99热这里只有精品18| 国产视频首页在线观看| 精品欧美国产一区二区三| 成年免费大片在线观看| 别揉我奶头 嗯啊视频| 免费一级毛片在线播放高清视频| 欧美日本视频| 中国美白少妇内射xxxbb| 婷婷精品国产亚洲av| av在线观看视频网站免费| 高清日韩中文字幕在线| 在线观看美女被高潮喷水网站| 国产成年人精品一区二区| 九九热线精品视视频播放| 亚洲av.av天堂| 国产高清有码在线观看视频| 男插女下体视频免费在线播放| 午夜精品在线福利| 久久草成人影院| 国产综合懂色| 免费观看精品视频网站| kizo精华| 国产精品国产高清国产av| 亚洲中文字幕日韩| 国产爱豆传媒在线观看| 亚洲自拍偷在线| 久久久久久伊人网av| 国产av不卡久久| 亚洲国产欧美人成| 国产黄色小视频在线观看| 欧美xxxx黑人xx丫x性爽| 国产老妇伦熟女老妇高清| 国产中年淑女户外野战色| 精品99又大又爽又粗少妇毛片| 精品久久久久久久久久免费视频| 国产在线男女| 白带黄色成豆腐渣| 国产一级毛片七仙女欲春2| 黄片wwwwww| 久久精品国产亚洲av香蕉五月| 国产一区二区亚洲精品在线观看| av在线播放精品| 欧美日韩在线观看h| 久久精品久久久久久噜噜老黄 | 欧美日韩在线观看h| 久久99热这里只有精品18| 久久久久久国产a免费观看| 18禁裸乳无遮挡免费网站照片| 91精品国产九色| 国产久久久一区二区三区| 成人午夜精彩视频在线观看| 国产亚洲精品av在线| 国产精品.久久久| 人妻夜夜爽99麻豆av| 日日啪夜夜撸| 日日撸夜夜添| 听说在线观看完整版免费高清| 午夜精品在线福利| 国产单亲对白刺激| 老女人水多毛片| 两性午夜刺激爽爽歪歪视频在线观看| 你懂的网址亚洲精品在线观看 | 国内精品美女久久久久久| 国产精品久久久久久亚洲av鲁大| 日日撸夜夜添| 欧美成人a在线观看| 国产亚洲欧美98| 国产亚洲av片在线观看秒播厂 | 成人毛片60女人毛片免费| 中文欧美无线码| 亚洲内射少妇av| 中国美女看黄片| 最近手机中文字幕大全| 晚上一个人看的免费电影| 亚洲欧美日韩高清在线视频| 国产成人一区二区在线| 精品不卡国产一区二区三区| 亚洲精品乱码久久久久久按摩| 久久亚洲精品不卡| 国产成人精品婷婷| 99九九线精品视频在线观看视频| 性插视频无遮挡在线免费观看| 麻豆成人午夜福利视频| 午夜精品国产一区二区电影 | 免费搜索国产男女视频| 好男人在线观看高清免费视频| 男人狂女人下面高潮的视频| 精品久久久久久久久久免费视频| 日韩,欧美,国产一区二区三区 | 日韩,欧美,国产一区二区三区 | 人妻制服诱惑在线中文字幕| 夫妻性生交免费视频一级片| 日本五十路高清| 亚洲一级一片aⅴ在线观看| 亚洲av第一区精品v没综合| 午夜老司机福利剧场| 亚洲精品456在线播放app| 啦啦啦韩国在线观看视频| 特级一级黄色大片| 亚洲成人久久性| 99久久中文字幕三级久久日本| 亚洲精品国产成人久久av| 给我免费播放毛片高清在线观看| 内地一区二区视频在线| 成人午夜精彩视频在线观看| 国产成人精品久久久久久| 26uuu在线亚洲综合色| av在线天堂中文字幕| 18禁在线无遮挡免费观看视频| 欧美3d第一页| 日韩大尺度精品在线看网址| 黄色日韩在线| 哪个播放器可以免费观看大片| 久久久久久久久久黄片| 美女xxoo啪啪120秒动态图| 亚洲av熟女| 亚洲国产精品成人综合色| 成人性生交大片免费视频hd| 国产片特级美女逼逼视频| 99久久成人亚洲精品观看| 日韩国内少妇激情av| 一进一出抽搐gif免费好疼| 男女下面进入的视频免费午夜| 亚洲欧洲日产国产| 久久午夜亚洲精品久久| 人妻久久中文字幕网| 久久这里有精品视频免费| 国产精品久久久久久久电影| 美女大奶头视频| 国产一区二区在线av高清观看| 给我免费播放毛片高清在线观看| 亚洲不卡免费看| 人人妻人人看人人澡| av在线蜜桃| 亚洲av免费在线观看| 欧美xxxx黑人xx丫x性爽| 九九爱精品视频在线观看| 精品少妇黑人巨大在线播放 | 国产精品久久久久久久电影| 日本色播在线视频| 一夜夜www| 麻豆av噜噜一区二区三区| 国产精品人妻久久久久久| 久久久久久九九精品二区国产| 亚洲自偷自拍三级| 人妻少妇偷人精品九色| 啦啦啦韩国在线观看视频| 22中文网久久字幕| 白带黄色成豆腐渣| 亚洲自拍偷在线| 国产精品爽爽va在线观看网站| 精品一区二区三区视频在线| 久久久a久久爽久久v久久| 精品不卡国产一区二区三区| 少妇被粗大猛烈的视频| 我要搜黄色片| 日韩 亚洲 欧美在线| 亚洲第一电影网av| 国产精品一区二区在线观看99 | 观看美女的网站| 韩国av在线不卡| 亚洲人成网站在线观看播放| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久视频播放| 亚洲美女搞黄在线观看| 国产免费男女视频| 男女啪啪激烈高潮av片| 久久精品影院6| 亚洲精品久久国产高清桃花| 欧美成人免费av一区二区三区| 亚洲在久久综合| 精华霜和精华液先用哪个| 中国国产av一级| 久久久久性生活片| 联通29元200g的流量卡| 欧美日韩一区二区视频在线观看视频在线 | 少妇人妻精品综合一区二区 | 国产精品一区二区三区四区免费观看| 2021天堂中文幕一二区在线观| 超碰av人人做人人爽久久| 99热全是精品| 午夜免费男女啪啪视频观看| 国产av在哪里看| 舔av片在线| 秋霞在线观看毛片| 人人妻人人看人人澡| 99国产极品粉嫩在线观看| a级毛片a级免费在线| 麻豆精品久久久久久蜜桃| 国产午夜精品论理片| 亚洲欧美日韩无卡精品| 69人妻影院| 级片在线观看| 一进一出抽搐gif免费好疼| 欧美性猛交黑人性爽| 国产精品久久电影中文字幕| 欧美区成人在线视频| 国产一区亚洲一区在线观看| 亚洲四区av| 国产伦精品一区二区三区四那| 欧美又色又爽又黄视频| 国产黄色小视频在线观看| 久久精品国产亚洲av天美| 免费观看在线日韩| 一边亲一边摸免费视频| 亚洲欧美日韩高清专用| 欧美变态另类bdsm刘玥| 人人妻人人澡欧美一区二区| 18禁在线播放成人免费| 91精品国产九色| 看十八女毛片水多多多| 爱豆传媒免费全集在线观看| 国产美女午夜福利| 国产午夜精品论理片| 3wmmmm亚洲av在线观看| 一级毛片我不卡| 国产日韩欧美在线精品| 波多野结衣高清无吗| 亚洲av一区综合| 一级av片app| 一级黄色大片毛片| 色综合亚洲欧美另类图片| 久久久久九九精品影院| 欧美最黄视频在线播放免费| 久久99蜜桃精品久久| 亚洲精品影视一区二区三区av| 欧美一级a爱片免费观看看| 此物有八面人人有两片| 中文字幕熟女人妻在线| 亚洲av不卡在线观看| 欧美日韩综合久久久久久| 国产av不卡久久| 蜜桃亚洲精品一区二区三区| 美女脱内裤让男人舔精品视频 | 精品免费久久久久久久清纯| 禁无遮挡网站| 国产精品麻豆人妻色哟哟久久 | 国产精品一及| 在线观看午夜福利视频| 中文字幕久久专区| 亚洲成人精品中文字幕电影| 成人午夜高清在线视频| 精品久久久久久久久久免费视频| av在线观看视频网站免费| 91久久精品国产一区二区成人| 日本一本二区三区精品| 波多野结衣高清无吗| 亚洲五月天丁香| 99热精品在线国产| 99久久精品热视频| 国产美女午夜福利| 国产精品一区二区三区四区免费观看| 色吧在线观看| 麻豆av噜噜一区二区三区| 变态另类成人亚洲欧美熟女| av.在线天堂| 久久人人爽人人爽人人片va| 久久99热6这里只有精品| 国产高清不卡午夜福利| 一边亲一边摸免费视频| 久久草成人影院| 日本黄大片高清| 淫秽高清视频在线观看| 蜜桃久久精品国产亚洲av| 欧美一区二区精品小视频在线| 最近中文字幕高清免费大全6| 精品久久久久久久久久久久久| 欧美+亚洲+日韩+国产| 岛国在线免费视频观看| 国产成人精品一,二区 | 夜夜爽天天搞| 美女被艹到高潮喷水动态| 精品国产三级普通话版| 精品人妻熟女av久视频| 欧美激情久久久久久爽电影| 欧美性猛交黑人性爽| 男女那种视频在线观看| 亚洲综合色惰| 麻豆av噜噜一区二区三区| 国内精品美女久久久久久| 69人妻影院| 99热精品在线国产| 九九久久精品国产亚洲av麻豆| 又爽又黄无遮挡网站| 国产高潮美女av| 男女视频在线观看网站免费| 成人欧美大片| 久久婷婷人人爽人人干人人爱| 国产中年淑女户外野战色| 亚洲中文字幕日韩| 成人午夜高清在线视频| 国产精品野战在线观看| 国产精品av视频在线免费观看| 秋霞在线观看毛片| 全区人妻精品视频| 亚洲七黄色美女视频| 成人一区二区视频在线观看| 亚洲美女视频黄频| 久久人人爽人人片av| 久久精品91蜜桃| av在线播放精品| 欧美变态另类bdsm刘玥| 国产一区亚洲一区在线观看| 欧美丝袜亚洲另类| av天堂中文字幕网| 国产精品女同一区二区软件| 欧美一级a爱片免费观看看| 欧美bdsm另类| 青春草国产在线视频 | 国产真实伦视频高清在线观看| 国产片特级美女逼逼视频| 国产黄片美女视频| 听说在线观看完整版免费高清| 亚洲在线自拍视频| 久久精品久久久久久噜噜老黄 | 乱系列少妇在线播放| 国产高清不卡午夜福利| 最近中文字幕高清免费大全6| 小蜜桃在线观看免费完整版高清| 我的女老师完整版在线观看| 欧美成人免费av一区二区三区| 男人舔奶头视频| 久久久国产成人免费| av天堂中文字幕网| 欧美日韩乱码在线| 国产一区二区在线av高清观看| 亚洲欧美精品综合久久99| 亚洲精品久久久久久婷婷小说 | avwww免费| 干丝袜人妻中文字幕| 日本免费一区二区三区高清不卡| 男女视频在线观看网站免费| 国产伦在线观看视频一区| 国产av不卡久久| 精品久久久久久久久久久久久| 免费观看a级毛片全部| 麻豆乱淫一区二区| 国产女主播在线喷水免费视频网站 | 1024手机看黄色片| 小说图片视频综合网站| 性插视频无遮挡在线免费观看| 人妻制服诱惑在线中文字幕| 舔av片在线| 成人欧美大片| 精品午夜福利在线看| 亚洲久久久久久中文字幕| 中文字幕av成人在线电影| 久久久色成人| 美女脱内裤让男人舔精品视频 | 麻豆精品久久久久久蜜桃| 两个人的视频大全免费| 日韩欧美精品免费久久| 久久久久久久久久久免费av| 卡戴珊不雅视频在线播放| 国语自产精品视频在线第100页| 国产成人a∨麻豆精品| 九九在线视频观看精品| 毛片一级片免费看久久久久| 九草在线视频观看| 欧美日韩乱码在线| 中文字幕av在线有码专区| 一边亲一边摸免费视频| 最近最新中文字幕大全电影3| 菩萨蛮人人尽说江南好唐韦庄 | 深爱激情五月婷婷| 看黄色毛片网站| 亚洲精品粉嫩美女一区| 精品久久久噜噜| 国产精品一区二区三区四区久久| 中文在线观看免费www的网站| 久久亚洲精品不卡| 91精品国产九色| 黄片无遮挡物在线观看| 桃色一区二区三区在线观看| 此物有八面人人有两片| 一进一出抽搐动态| 我的老师免费观看完整版| 国内精品宾馆在线| 成年女人看的毛片在线观看| 日韩欧美精品免费久久| 亚洲成人av在线免费| 亚洲高清免费不卡视频| 久久久久久久久久久免费av| 久久6这里有精品| 看免费成人av毛片| 亚洲av熟女| 岛国毛片在线播放| 99国产精品一区二区蜜桃av| 成人高潮视频无遮挡免费网站| 午夜亚洲福利在线播放| 欧美日韩国产亚洲二区| 超碰av人人做人人爽久久| 国产午夜精品论理片| 成人漫画全彩无遮挡| 亚洲成av人片在线播放无| 日韩强制内射视频| 国产精品人妻久久久久久| 综合色av麻豆| 久久亚洲国产成人精品v| 超碰av人人做人人爽久久| 欧美+日韩+精品| 少妇裸体淫交视频免费看高清| 国产伦理片在线播放av一区 | 一本一本综合久久| 永久网站在线| 亚洲国产色片| 91午夜精品亚洲一区二区三区| 可以在线观看的亚洲视频| 日韩人妻高清精品专区| 最近2019中文字幕mv第一页| 国内少妇人妻偷人精品xxx网站| 能在线免费观看的黄片| 国产爱豆传媒在线观看| 国产美女午夜福利| 亚洲精品自拍成人| 日韩一区二区三区影片| 国产一区二区三区av在线 | 欧美+亚洲+日韩+国产| 三级毛片av免费| 最好的美女福利视频网| 亚洲一区高清亚洲精品| av国产免费在线观看| 欧美另类亚洲清纯唯美| 成人无遮挡网站| 人体艺术视频欧美日本| 麻豆成人午夜福利视频| 国产精品久久久久久亚洲av鲁大| 99热这里只有是精品在线观看| 亚洲自拍偷在线| 亚洲欧美精品综合久久99| 黄片wwwwww| 尤物成人国产欧美一区二区三区| 深夜精品福利|