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

    基于傅立葉平均法下的連續(xù)小推力動力學分析①

    2012-07-09 09:11:50馮維明
    固體火箭技術 2012年3期
    關鍵詞:傅立葉根數(shù)級數(shù)

    馮維明,李 源,苗 楠

    (山東大學 工程力學系,濟南 250062)

    基于傅立葉平均法下的連續(xù)小推力動力學分析①

    馮維明,李 源,苗 楠

    (山東大學 工程力學系,濟南 250062)

    通過將小推力展開為偏近點角的傅立葉級數(shù),并對高斯攝動方程在一個軌道周期上的平均,將原方程的推力轉化為僅由14個傅立葉系數(shù)表示的控制變量。仿真計算表明,平均化后的高斯方程使計算量與牛頓積分相比顯著減少,且對小推力而言有足夠的精度。對利用平均化后的高斯方程計算軌道根數(shù)時產(chǎn)生誤差的原因進行了研究,并進一步分析小推力的范圍和小推力近似表達式對上述誤差的影響,為今后小推力下非開普勒軌道動力學分析提供了理論依據(jù)和參數(shù)。

    傅立葉級數(shù);連續(xù)小推力;平均法;高斯方程;動力學分析;非開普勒軌道

    0 引言

    小推力推進系統(tǒng)為許多星際航行和地球軌道任務提供了一種高效的新選擇,但是這種系統(tǒng)卻對最優(yōu)化控制提出了新的挑戰(zhàn)[1]對于某些特殊的小推力作用的情況下,軌道轉移的最優(yōu)控制問題已經(jīng)出現(xiàn)了解析法和數(shù)值近似解的方法[2-3],但是對于一般的小推力問題,需要對每一個初始條件和推力值進行完整的數(shù)值積分,軌道參數(shù)的確定往往對這些變量很敏感。因此,普通的最優(yōu)控制法則對于確定螺旋軌道的數(shù)十甚至上百圈軌道是十分困難的[4]。

    在一些特殊情形的軌道轉移中,數(shù)值解法已經(jīng)得到了挖掘。目前較多的解法應用了變分法或直接優(yōu)化法[5-7],以確定特定約束條件下的小推力的最優(yōu)控制律。此外,文獻[8]運用李亞普諾夫反饋控制的方法,用以解決開環(huán)軌道最優(yōu)時間控制和最優(yōu)軌道問題。然而實踐證明,把平均法和其他方法結合的思路在解決對初始軌道和推力變量預測敏感的問題上是很有效的[9-10],但所有這些方法都局限于特定的推力和軌道參數(shù)范圍。

    本文討論了一種有效解決航天器在小推力下軌道動力學問題的新方法。推力加速度的各分量以偏近點角展開為傅立葉級數(shù),而高斯變分方程則在一個軌道周期上進行平均并且通過正交條件進行簡化,從而定義出了一組久期方程。這組久期方程是一個含有14個推力傅立葉系數(shù)的函數(shù)(無論初始傅立葉級數(shù)的階數(shù)是多少)。因此,該方法將普遍形式的推力分布情況簡化為僅有14個參數(shù)的形式,這使計算量與牛頓積分相比大大減少。并以此分析了影響傅立葉系數(shù)表示的平均久期方程求解精度的主要因素。

    1 數(shù)學建模

    考慮一個在軌航天器,其質量對于其環(huán)繞的中心體可忽略不計。將航天器認為是一個質點。它受一個大小和方向可能隨時間變化的連續(xù)推力加速度的作用。其軌道可用牛頓運動方程來表述:

    力加速度矢量F可以沿著徑向FR、法向FW和周向FS分解,即

    式中r和w分別為徑向和法向的單位矢量。

    牛頓方程可分解為拉格朗日行星方程,高斯形式的拉格朗日行星方程可表示為

    式中 a為半長軸;e為偏心率;i為軌道傾角;Ω為升交點赤經(jīng);ω為近地點幅角;f為真近點角;E為偏近點角,而ε1+∫ndt=l為平黃經(jīng)。

    平近點角是平黃經(jīng)和近地點幅角的差值:

    根據(jù)傅立葉理論,對于任何一個在(0,L)上只存在有限個跳躍間斷點的分段光滑的函數(shù)f(θ),都可以表示為一系列周期性延伸的傅立葉級數(shù),這個級數(shù)無限逼近于函數(shù)本身,當跳躍間斷點存在時,傅立葉級數(shù)將逼近于左右極限的平均值。因此,這種表示方法可以應用于幾乎任何一種普通的小推力航天器的控制律中。對于給定的任意加速度矢量F,其各個分量都可以展開為任意時間間隔上的傅立葉級數(shù)。實際上,傅立葉級數(shù)可按時間展開,也可按隨時間變化的軌道參數(shù)展開,如真近點角、偏近點角或平近點角??紤]一個軌道周期L=2π,以θ來表示這一任意的軌道參數(shù):

    現(xiàn)在開始一階平均化分析,先假設一個加速度矢量,它可在一個軌道周期(L=2π)中表示,這個加速度矢量的數(shù)值足夠小,以使軌道的大小和形狀在一圈中的改變不是很明顯。因此,可將高斯方程相對于平近點角在一個軌道周期上平均來得到平均軌道根數(shù)。

    式中 γ代表任意的軌道根數(shù)。

    在對推力加速度矢量的各分量進行傅立葉展開時,所使用的軌道參數(shù)的選擇是很重要的,有的參數(shù)會使由此得出的久期方程將變得冗長而復雜(如按真近點角展開成傅立葉級數(shù)的話),而有的參數(shù)會使久期方程會被簡化。經(jīng)比較,加速度矢量的每一個分量都是按偏近點角進行傅立葉級數(shù)展開是最佳選擇,并且各軌道參數(shù)的平均化也是按偏近點角進行,由dM=(1-ecosE)dE,式(13)可寫為

    將式(11)代入式(4)~式(9)中,并由式(14)可得到關于偏近點角的平均化高斯方程。由傅立葉級數(shù)的正交性,可消掉推力加速度傅立葉展開級數(shù)中第二階以上及個別第二階傅立葉系數(shù)。因此,無論初始推力加速度傅立葉級數(shù)的階數(shù)如何,軌道根數(shù)a、e、i、Ω、ω和ε1的平均變化率只與這14個傅立葉系數(shù)有關,即(k=0,1,2)(k=0,1,2)(k=0,1,2),(k=1,2)(1,2)則高斯平均方程為

    2 計算分析

    為在計算時數(shù)據(jù)的簡便起見,將方程做歸一化處理。在此,將重力參數(shù)μ歸一化常數(shù),μ=3.986×105km3/s2=3.986×1014m3/s2,因此相當于各變量的長度(以m為單位)除以。初始的軌道根數(shù)由表1給出。

    表1 初始軌道根數(shù)Table 1 initial orbit elements

    設控制率為周向加速度FS=1×10-8×(1-sinE)、徑向加速度FR=1×10-8×sinE、法向加速度FW=1×10-8×(1+cosE),上述推力的加速度分量為歸一化后的值,如航天器為1 000 kg,推力的最大值幅值為2.208 N,可算作小推力范圍。用龍格-庫塔法進行積分來估算軌道根數(shù)的變化情況。該情況下由牛頓方程和平均變分方程求解出的軌道根數(shù)分別用實線和虛線在圖1中給出。

    圖1 航天器受到3個方向推力加速度作用時軌道根數(shù)吻合情況Fig.1 Osculating orbital elements of spacecraft subject to acceleration in three directions

    由圖1可看到,所有的軌道根數(shù)均給出了非常吻合的結果。即用傅立葉系數(shù)表示的平均久期方程的求解結果和直接對拉格朗日方程進行積分得到的結果具有緊密的一致性。

    再考慮一個簡單的周向分步加速度的問題,F(xiàn)R=FW=0,F(xiàn)S按如圖2中所示變化。每個周期中包含2次點火和2個滑行弧,初始軌道同表1所示。切向分步推力加速度作用時軌道根數(shù)吻合情況見圖3。

    圖2 周向分步加速度Fig.2 Step tangential acceleration

    圖3 切向分步推力加速度作用時軌道根數(shù)吻合情況Fig.3 Osculating orbital elements of spacecraft subject to step tangential acceleration

    由圖3知,平均久期方程和牛頓方程給出的結果具有完全相同的趨勢,數(shù)值上也相差不大。

    另外,基于傅立葉平均化后的高斯方程最顯著優(yōu)勢之一就是計算量大規(guī)模減少,尤其在進行非開普勒軌道優(yōu)化計算時,計算精度和計算耗時是困擾著眾多學者的問題,考慮在表2所列初始條件下,歸一化后的周向推力為FS=5×10-9(實際推力加速度為7.36×10-5m/s2),分別用牛頓方程和平均變分方程求解出的軌道根數(shù)及誤差列入表2。

    推力表達式為如下形式:

    表2 連續(xù)小推力作用下牛頓方程和平均方程10圈后的計算結果比較Table 2 Comparison of the calculated results determined by the two methods on low-thrust continuous controls after running 10 laps

    航天器運行時間為65 973.4 s,約運行10圈。由表2計算結果看出,由平均化后的高斯方程計算得到的10圈后的軌道的6個根與牛頓方程精確積分所得到的結果相比,除近地點幅角ω相對誤差略大些,其他相對誤差非常之小。利用牛頓法進行精確積分計算耗時(CPU耗時)為13.556 407 s,而利用基于傅立葉平均化后的高斯方程計算耗時為0.450 808 s,后者比前者快30倍。

    3 影響計算精度的主要因素

    平均久期方程能夠精確有效解決小推力螺旋軌道問題,且與通常的牛頓問題的積分方法相比能大大減少計算量。但每一種近似計算方法都有其局限性,對于平均久期方程也是如此,這也是很多文獻所未曾涉及到的問題,因此很有必要對利用該方法進行軌道計算可能導致的精度進行分析。

    3.1 推力幅值大小的影響

    由該方法的限制條件可知,推力幅值過大,會導致周期內軌道變形太大,從而使平均法的結果產(chǎn)生累計誤差,于是,關心在保證精度條件下,推力的極限值是多少。為討論推力增大對計算結果帶來的影響,研究推力加速度大小在切向分步加速度形式(如圖2示)的控制率中對計算結果的影響。令周向加速度FS的幅值分別增大為 1 ×10-7(歸一化后)的 1.0、1.5、2.0和2.5倍,限于篇幅僅討論對軌道形狀影響最大的根數(shù)解半長軸a和偏心率e,圖4給出其吻合情況(幅值為1×10-7見圖3)。

    圖4 不同切向分步推力加速度作用時軌道根數(shù)吻合情況Fig.4 Osculating orbital elements of spacecraft subject to different step tangential acceleration

    為更清晰地表現(xiàn)牛頓方程與平均方程的計算結果誤差,現(xiàn)將不同幅值的推力作用下航天器運行20圈2種方法計算得到的半長軸a、偏心率e及相對誤差列入表3。

    由圖4和表3可看到,隨著加速度幅值的增大,a和e的誤差都在增大,a的誤差增幅較小,而e的誤差顯著增大。在FS達到2.5×10-7時,偏心率的2種求解結果已有顯著偏離,此時對應質量為1 000 kg的航天器,推力大小約為18.4 N,約20圈后,偏心率的結果已經(jīng)完全不可靠了。另外,由圖4也可看到,即使推力較大時,在前12圈之內,平均久期方程和牛頓方程給出的結果仍然保持著良好的一致性。由表3可見,半長軸變化很大,因此,若在一定圈數(shù)內完成變軌任務,平均方程仍有足夠的可靠度,即小推力的上限應根據(jù)變軌任務確定。

    表3 推力加速度大小對a、e計算誤差的影響Table 3 Influence of acceleration on the accuracy of a and e

    3.2 近似推力表達式的影響

    由式(15)~式(20)可知,平均化后的高斯方程僅與14個傅立葉系數(shù)有關,推力的傅立葉級數(shù)高階項會因平均過程而被消去,平均化后非零的傅立葉系數(shù)是否能精確的描述原推力形式將成為關鍵。在本文第1個算例中,由傅立葉系數(shù)表達式(12)不難得到各推力加速度矢量分量的傅立葉表達式中的系數(shù)為=1 ×10-8=1 ×10-8=-1 ×10-8=1 ×10-8和=1×10-8,而其他系數(shù)均為零。將上述非零系數(shù)帶入式(11)可得

    高斯方程經(jīng)平均化后得到的推力加速度分量傅立葉級數(shù)近似表達式與原初始推力表達式完全一致,因此傅立葉級數(shù)近似表達式對計算精度沒有影響。下面再來考慮3個推力分量皆為分步加速度的形式,即

    在圖5中可表示加速度分量近似表達式與原表達式的差異。圖5中用點表示的線是高斯方程平均后3個加速度分量變化曲線,注意到此時徑向推力FR已變成常數(shù),為分段推力的平均值;而周向推力FS和法向推力FW近似為正弦函數(shù)形式,從圖形上看比較接近原推力形式。這也就是在前例中(僅有周向推力FS)計算精度較高的原因之一。平均方程與解析解得到的計算結果如圖6所示。

    圖5 分段切向加速度的傅立葉級數(shù)表示Fig.5 Fourier series for step circumferential acceleration

    圖6 三向分步推力加速度作用時軌道根數(shù)吻合情況Fig.6 Osculating orbital elements of spacecraft subject to step circumferential acceleration in three directions

    圖6表示的是航天器運行了10圈的根數(shù)時程圖。長半軸a和偏心率e平均方程的解與牛頓積分法的解吻合的較好,軌道傾角i,升交點赤經(jīng)Ω,近地點幅角ω在10圈后誤差開始增大,較明顯的為Ω的值,平均結果開始偏離高斯方程積分結果。

    4 結論

    (1)將推力加速度分量展開為偏近點角表示的傅立葉級數(shù),而高斯變分方程則在一個軌道周期上進行平均并且通過正交條件進行簡化。該方法將整個連續(xù)控制問題的參數(shù)減少為14個(無論初始傅立葉級數(shù)的階數(shù)是多少),這使計算量與牛頓積分相比大大減少。

    (2)通過平面變軌和空間變軌的計算結果分析,用該方法計算連續(xù)和非連續(xù)小推力下的軌道轉移是正確的和準確的。

    (3)對造成計算誤差的主要因素進行了定性和定量分析,推力超出“小推力”的范圍和傅立葉級數(shù)下推力的近似表達式(高斯方程平均后)是影響精度的主要因素,即便如此,在一定范圍內平均法計算結果仍與牛頓方程的計算結果吻合良好。

    [1]Gao Y.Advances in low-thrust trajectory optimization and flight mechanics ,dissertation thesis[R].University of Missouri-Columbia,2003,64-12(B):6178.

    [2]Craig A Kluever.Direct approach for computing near-optimal low-thrust earth-orbit transfers[J].Journal of Spacecraft and Rockets,1998,29(1):45-61.

    [3]Akella M R,Broucke R A.Anatomy of the constant radial thrust problem[J].Journal of Guidance,Control,and Dy namics,2002,25(3):563-570.

    [4]李俊峰,龔勝平.非開普勒軌道動力學與控制[J].宇航學報,2009,30(1):47-53.

    [5]Kluever C A.Optimal low-thrust interplanetary trajectories by direct method techniques[J].Journal of the Astronautical Sciences,1997,45(3):162-247.

    [6]Petropoulos A E.Some analytic integrals of the averaged variational equations for a thrusting spacecraft[R].Interplanetary Network Progress Rept.42-150,Jet Propulsion Lab.,California Inst.of Technology,Pasadena,CA,Aug.2002:1-29.

    [7]John T Betts.Very low-thrust trajectory optimization using a direct SQP method[J].Journal of Computational and Applied Mathematics,2000,120:27-40.

    [8]Gurfill P.Nonlinear feedback control of low-thrust orbital transfer in a central graviational field[J].Acta Astronautica,2007,60(8-9):631-648.

    [9]Jennifer S Hudson,Daniel J Scheeres.Reduction of lowthrust continuous controls for trajectory dynamics[J].Journal of Guidance,Control,and Dynamics,2009,32(3):780-787.

    [10]尚海濱,崔平遠,欒恩杰.基于平均法的小推力轉移軌道優(yōu)化研究[C]//25屆中國控制論文集.2006.

    Dynamic analysis of continuous low-thrust based on fourier average method

    FENG Wei-ming,LI Yuan,MIAO Nan
    (Department of Engineering Mechanics,Shandong University,Jinan 250062,China)

    Each component of the thrust vector was expanded as Fourier series in eccentric anomaly and Gauss variational equations were averaged over one orbit period,then the thrust vector was translated to a variable controlled by fourteen Fourier's parameters.Simulation results show that these secular equations are sufficient to accurately determine a low-thrust spiral trajectory with significantly reduced computation as compared with integration of the full Newtonian problem.In addition,error causes of orbit elements witch were calculated by the averaged Gauss equations were studied,and influence of low-thrust range and approximate expressions on the error was further analyzed,providing theoretical basis and parameters for dynamic analysis of the Non-Keplerian orbits of lowthrust.

    Fourier series;continuous low-thrust;average method;Gauss equations;dynamic analysis;non-Keplerian orbits

    V412

    A

    1006-2793(2012)03-0285-05

    2011-08-18;;

    2011-10-10。

    國家863項目。

    馮維明(1957—),男,教授,主要研究方向為非線性動力學和軌道動力學。E-mail:fwm@sdu.edu.cn

    (編輯:呂耀輝)

    猜你喜歡
    傅立葉根數(shù)級數(shù)
    更正
    中國棉花(2022年2期)2022-11-23 20:02:04
    尋找規(guī)律巧算根數(shù)
    不同坐標系下傅立葉變換性質
    三角函數(shù)的傅立葉變換推導公式
    電力電纜非金屬外護套材料的傅立葉紅外光譜鑒別法
    電線電纜(2018年6期)2019-01-08 02:45:36
    Dirichlet級數(shù)及其Dirichlet-Hadamard乘積的增長性
    基于傅立葉變換的CT系統(tǒng)參數(shù)標定成像方法探究
    智富時代(2018年7期)2018-09-03 03:47:26
    玉米的胡須
    幾個常數(shù)項級數(shù)的和
    基于改進型號第二婁無廳點根數(shù)的北斗CEO衛(wèi)星廣播星歷擬合算法及實現(xiàn)
    成人高潮视频无遮挡免费网站| 国产探花极品一区二区| 国产精品一二三区在线看| 最近中文字幕高清免费大全6| 汤姆久久久久久久影院中文字幕 | 内地一区二区视频在线| 校园人妻丝袜中文字幕| 精华霜和精华液先用哪个| av专区在线播放| 精品久久久久久久末码| 亚洲真实伦在线观看| 免费黄频网站在线观看国产| 久久国内精品自在自线图片| 亚洲18禁久久av| 成人毛片a级毛片在线播放| 精品国内亚洲2022精品成人| 老司机影院成人| 汤姆久久久久久久影院中文字幕 | 国产精品麻豆人妻色哟哟久久 | 国产视频内射| 久久精品综合一区二区三区| 九九在线视频观看精品| 成人鲁丝片一二三区免费| 一边亲一边摸免费视频| 91久久精品电影网| 乱码一卡2卡4卡精品| 日韩欧美精品免费久久| 久久99精品国语久久久| 国产伦在线观看视频一区| 国产精品美女特级片免费视频播放器| 精品久久久久久久久av| 一区二区三区高清视频在线| 国产乱人偷精品视频| 18禁裸乳无遮挡免费网站照片| 男女下面进入的视频免费午夜| 大香蕉久久网| 成人鲁丝片一二三区免费| 欧美性感艳星| av.在线天堂| 精品人妻一区二区三区麻豆| 春色校园在线视频观看| 久久人人爽人人爽人人片va| 最近中文字幕2019免费版| 亚洲不卡免费看| videossex国产| 免费播放大片免费观看视频在线观看| 亚洲欧美成人精品一区二区| 青春草国产在线视频| 亚洲精品国产成人久久av| freevideosex欧美| 精品酒店卫生间| 欧美性感艳星| 精品酒店卫生间| 国产乱人偷精品视频| 亚洲精品影视一区二区三区av| 黄片wwwwww| 最近中文字幕2019免费版| 成人欧美大片| 91精品伊人久久大香线蕉| 亚洲欧洲国产日韩| 亚洲av二区三区四区| 51国产日韩欧美| 欧美另类一区| 久久久午夜欧美精品| 大又大粗又爽又黄少妇毛片口| 亚洲不卡免费看| 插阴视频在线观看视频| 国产免费一级a男人的天堂| 成人美女网站在线观看视频| 亚洲av免费高清在线观看| 色综合站精品国产| 1000部很黄的大片| 免费看不卡的av| 亚洲四区av| 我的女老师完整版在线观看| 国产女主播在线喷水免费视频网站 | 亚洲精品影视一区二区三区av| 国产成人免费观看mmmm| 人体艺术视频欧美日本| 99热这里只有是精品50| 久久久色成人| 一级毛片久久久久久久久女| 少妇被粗大猛烈的视频| 色网站视频免费| 国产在视频线在精品| 熟女人妻精品中文字幕| 免费看美女性在线毛片视频| 亚洲成人精品中文字幕电影| h日本视频在线播放| 国产大屁股一区二区在线视频| 男插女下体视频免费在线播放| 天堂√8在线中文| 91aial.com中文字幕在线观看| 国产成人精品福利久久| 国产在视频线在精品| 久久国内精品自在自线图片| 亚洲欧美精品专区久久| 亚洲欧美一区二区三区黑人 | 色尼玛亚洲综合影院| 卡戴珊不雅视频在线播放| 亚洲内射少妇av| 日本与韩国留学比较| 九草在线视频观看| 亚洲三级黄色毛片| 熟女人妻精品中文字幕| 视频中文字幕在线观看| av免费在线看不卡| 国产免费福利视频在线观看| 日韩一区二区视频免费看| 亚洲精品国产av蜜桃| 日本一二三区视频观看| av一本久久久久| 哪个播放器可以免费观看大片| 精品国产一区二区三区久久久樱花 | 一边亲一边摸免费视频| 国产午夜精品久久久久久一区二区三区| 天天一区二区日本电影三级| 午夜久久久久精精品| 少妇人妻精品综合一区二区| 男女视频在线观看网站免费| 精品一区在线观看国产| 深夜a级毛片| 精品不卡国产一区二区三区| 国产午夜精品一二区理论片| av线在线观看网站| 少妇人妻精品综合一区二区| av在线播放精品| 久久久欧美国产精品| 亚洲精品日本国产第一区| 18禁在线无遮挡免费观看视频| 99九九线精品视频在线观看视频| 高清av免费在线| 欧美日韩在线观看h| 大话2 男鬼变身卡| 在线a可以看的网站| 久久鲁丝午夜福利片| 国产又色又爽无遮挡免| 国产精品一区www在线观看| 一夜夜www| 九九久久精品国产亚洲av麻豆| 国内精品美女久久久久久| 观看免费一级毛片| 亚洲一区高清亚洲精品| 国产精品熟女久久久久浪| 18禁在线播放成人免费| 国产视频首页在线观看| 一级毛片久久久久久久久女| 美女主播在线视频| 日韩一本色道免费dvd| 国产成人精品福利久久| 天堂av国产一区二区熟女人妻| 亚洲成人一二三区av| 久久久久久伊人网av| 少妇被粗大猛烈的视频| 精华霜和精华液先用哪个| 国产精品福利在线免费观看| 男女啪啪激烈高潮av片| 三级国产精品欧美在线观看| 成人亚洲精品一区在线观看 | 国产高清三级在线| 国产在视频线在精品| 国产日韩欧美在线精品| 日韩视频在线欧美| 一级毛片aaaaaa免费看小| 极品少妇高潮喷水抽搐| 亚洲精品国产av蜜桃| 高清在线视频一区二区三区| 久久久久性生活片| 大片免费播放器 马上看| 激情五月婷婷亚洲| 免费av毛片视频| 97人妻精品一区二区三区麻豆| 纵有疾风起免费观看全集完整版 | 三级男女做爰猛烈吃奶摸视频| 美女大奶头视频| 国内少妇人妻偷人精品xxx网站| 国产欧美日韩精品一区二区| 91精品一卡2卡3卡4卡| kizo精华| 精品酒店卫生间| 国产精品1区2区在线观看.| 人妻制服诱惑在线中文字幕| 成人亚洲精品一区在线观看 | 久久精品夜色国产| 成人av在线播放网站| 大香蕉97超碰在线| 国产免费视频播放在线视频 | 大陆偷拍与自拍| 国产午夜精品论理片| 精品国内亚洲2022精品成人| 久久久欧美国产精品| av国产免费在线观看| 欧美bdsm另类| 久久国内精品自在自线图片| 一级a做视频免费观看| 国产精品无大码| 纵有疾风起免费观看全集完整版 | 欧美丝袜亚洲另类| 99久久精品国产国产毛片| 久久97久久精品| 欧美精品国产亚洲| 国产伦精品一区二区三区视频9| 午夜福利在线在线| 2022亚洲国产成人精品| 国产精品国产三级专区第一集| 成年女人看的毛片在线观看| 欧美不卡视频在线免费观看| 日日摸夜夜添夜夜添av毛片| 看黄色毛片网站| 夜夜爽夜夜爽视频| 国产免费又黄又爽又色| 在线 av 中文字幕| 男人爽女人下面视频在线观看| 欧美成人a在线观看| 美女xxoo啪啪120秒动态图| www.色视频.com| 久久久久久伊人网av| 亚洲自偷自拍三级| eeuss影院久久| 99久久精品热视频| 亚洲怡红院男人天堂| 亚洲乱码一区二区免费版| 国产精品国产三级专区第一集| 国产高清不卡午夜福利| 91久久精品电影网| 精品国产一区二区三区久久久樱花 | 狂野欧美白嫩少妇大欣赏| 成人二区视频| 久久久精品免费免费高清| 伦理电影大哥的女人| 国产精品综合久久久久久久免费| 日韩中字成人| 精华霜和精华液先用哪个| 日韩大片免费观看网站| 久久午夜福利片| 日韩电影二区| 1000部很黄的大片| 丰满人妻一区二区三区视频av| 亚洲婷婷狠狠爱综合网| 三级国产精品片| 中文字幕亚洲精品专区| av福利片在线观看| 亚洲伊人久久精品综合| 一二三四中文在线观看免费高清| 如何舔出高潮| 看免费成人av毛片| 久久精品久久久久久噜噜老黄| 欧美成人a在线观看| 久久97久久精品| 一级a做视频免费观看| 丰满少妇做爰视频| 国产久久久一区二区三区| 老师上课跳d突然被开到最大视频| 欧美不卡视频在线免费观看| 亚洲熟女精品中文字幕| 国产精品人妻久久久影院| 久久99热这里只频精品6学生| 中文字幕制服av| 欧美变态另类bdsm刘玥| 亚洲欧美成人综合另类久久久| 亚洲av免费高清在线观看| 精品久久久久久久人妻蜜臀av| 五月伊人婷婷丁香| 中文乱码字字幕精品一区二区三区 | 亚洲av在线观看美女高潮| 自拍偷自拍亚洲精品老妇| 性色avwww在线观看| 赤兔流量卡办理| 26uuu在线亚洲综合色| 久久久久九九精品影院| 国产探花极品一区二区| 国产成人免费观看mmmm| 狂野欧美白嫩少妇大欣赏| 中文资源天堂在线| 色哟哟·www| av线在线观看网站| 日韩一本色道免费dvd| 免费观看无遮挡的男女| 亚洲性久久影院| 高清日韩中文字幕在线| 天堂俺去俺来也www色官网 | 国产不卡一卡二| 好男人视频免费观看在线| 中文字幕免费在线视频6| 综合色av麻豆| 最近中文字幕高清免费大全6| 亚洲精品第二区| 禁无遮挡网站| 又爽又黄a免费视频| 国产精品1区2区在线观看.| 亚洲自偷自拍三级| 乱系列少妇在线播放| 欧美变态另类bdsm刘玥| 寂寞人妻少妇视频99o| 22中文网久久字幕| 日韩av免费高清视频| 国产高清有码在线观看视频| 五月天丁香电影| 国产综合懂色| kizo精华| 人妻制服诱惑在线中文字幕| 亚洲第一区二区三区不卡| 天天躁夜夜躁狠狠久久av| 丰满人妻一区二区三区视频av| 久久97久久精品| 啦啦啦啦在线视频资源| av专区在线播放| 欧美最新免费一区二区三区| av在线观看视频网站免费| 国产单亲对白刺激| 成人欧美大片| 久久精品国产自在天天线| 中文乱码字字幕精品一区二区三区 | 亚洲国产欧美人成| 国产真实伦视频高清在线观看| 亚洲精品乱码久久久久久按摩| 日韩三级伦理在线观看| 久久韩国三级中文字幕| 干丝袜人妻中文字幕| av线在线观看网站| 久久这里有精品视频免费| 国产高清不卡午夜福利| 国产成人免费观看mmmm| 两个人视频免费观看高清| 国产精品久久视频播放| 插阴视频在线观看视频| 极品教师在线视频| 麻豆成人av视频| 男的添女的下面高潮视频| 日韩人妻高清精品专区| 在线观看av片永久免费下载| 日韩三级伦理在线观看| 精品久久久久久电影网| 久久99热这里只有精品18| 精品一区在线观看国产| 亚洲国产精品成人久久小说| 中文字幕人妻熟人妻熟丝袜美| 欧美不卡视频在线免费观看| 国产永久视频网站| 麻豆国产97在线/欧美| 一区二区三区四区激情视频| 大香蕉97超碰在线| 亚洲久久久久久中文字幕| 日日摸夜夜添夜夜添av毛片| 日韩一本色道免费dvd| 精品午夜福利在线看| 最后的刺客免费高清国语| 别揉我奶头 嗯啊视频| 成人午夜精彩视频在线观看| 久久精品人妻少妇| 午夜福利视频1000在线观看| 日日摸夜夜添夜夜添av毛片| 美女主播在线视频| av女优亚洲男人天堂| 国产成人一区二区在线| 亚洲av福利一区| 成年av动漫网址| 亚洲精品视频女| 能在线免费看毛片的网站| 国产伦一二天堂av在线观看| 日本爱情动作片www.在线观看| 亚洲欧美日韩卡通动漫| 美女黄网站色视频| 在线 av 中文字幕| 国产伦精品一区二区三区四那| 国产精品一区二区性色av| 国产 亚洲一区二区三区 | 99久久中文字幕三级久久日本| 国语对白做爰xxxⅹ性视频网站| 日日啪夜夜爽| 成人亚洲欧美一区二区av| 婷婷六月久久综合丁香| 午夜福利视频精品| 成年av动漫网址| 网址你懂的国产日韩在线| 日韩不卡一区二区三区视频在线| 少妇猛男粗大的猛烈进出视频 | 欧美日韩国产mv在线观看视频 | 神马国产精品三级电影在线观看| 夫妻午夜视频| 水蜜桃什么品种好| 国产伦精品一区二区三区四那| 国产精品综合久久久久久久免费| 2018国产大陆天天弄谢| 国产精品久久久久久精品电影小说 | 青春草视频在线免费观看| 精品一区二区三区人妻视频| 欧美精品国产亚洲| 亚洲国产精品成人久久小说| 成人一区二区视频在线观看| 国产永久视频网站| 国产高清国产精品国产三级 | 精品国产一区二区三区久久久樱花 | 亚洲第一区二区三区不卡| 国产精品一区二区三区四区久久| 毛片女人毛片| 精品一区二区三卡| 亚洲精品乱码久久久久久按摩| 我要看日韩黄色一级片| 亚洲精品第二区| 美女xxoo啪啪120秒动态图| 18禁在线无遮挡免费观看视频| 久久久精品94久久精品| 少妇人妻精品综合一区二区| 免费av观看视频| 国产精品福利在线免费观看| 亚洲av一区综合| 国产精品久久久久久久久免| 亚洲欧美日韩东京热| 亚洲精品成人久久久久久| 69av精品久久久久久| 别揉我奶头 嗯啊视频| 欧美最新免费一区二区三区| 国产精品国产三级国产av玫瑰| 在现免费观看毛片| 免费人成在线观看视频色| ponron亚洲| 我的女老师完整版在线观看| 亚洲精品国产av成人精品| 真实男女啪啪啪动态图| 国产成人精品久久久久久| 99久久精品热视频| 久久国内精品自在自线图片| 大又大粗又爽又黄少妇毛片口| 国产永久视频网站| 我要看日韩黄色一级片| 国产日韩欧美在线精品| 少妇被粗大猛烈的视频| 欧美成人一区二区免费高清观看| 国产午夜精品论理片| 午夜激情久久久久久久| 晚上一个人看的免费电影| 美女cb高潮喷水在线观看| 亚洲,欧美,日韩| 高清日韩中文字幕在线| 国产精品久久久久久久久免| 日韩一本色道免费dvd| 毛片女人毛片| 精品国产三级普通话版| 免费无遮挡裸体视频| 婷婷六月久久综合丁香| 亚洲欧美清纯卡通| 国产淫语在线视频| 久久久亚洲精品成人影院| 久久综合国产亚洲精品| 建设人人有责人人尽责人人享有的 | 91精品伊人久久大香线蕉| 又爽又黄a免费视频| 三级国产精品欧美在线观看| av福利片在线观看| 国产一区二区在线观看日韩| 色网站视频免费| 黑人高潮一二区| 久久人人爽人人片av| 国内精品一区二区在线观看| 国产 亚洲一区二区三区 | 国语对白做爰xxxⅹ性视频网站| 日韩视频在线欧美| av在线观看视频网站免费| 丝袜美腿在线中文| 99热这里只有是精品50| 亚洲精品久久久久久婷婷小说| 日日摸夜夜添夜夜爱| 一级爰片在线观看| 日产精品乱码卡一卡2卡三| 国产精品久久久久久av不卡| av专区在线播放| 3wmmmm亚洲av在线观看| 91精品国产九色| 欧美日韩在线观看h| 女人被狂操c到高潮| 亚洲熟女精品中文字幕| 国产精品一二三区在线看| 久久久久久九九精品二区国产| 欧美成人一区二区免费高清观看| 精品久久久久久久人妻蜜臀av| 日韩不卡一区二区三区视频在线| 国产老妇伦熟女老妇高清| 久久久久久久久久成人| 最近最新中文字幕免费大全7| 少妇人妻精品综合一区二区| 精品酒店卫生间| 久久久国产一区二区| 一级av片app| 日本-黄色视频高清免费观看| 嫩草影院精品99| 91精品一卡2卡3卡4卡| 亚洲国产精品sss在线观看| 最近手机中文字幕大全| 国产亚洲精品av在线| 一级毛片 在线播放| 麻豆精品久久久久久蜜桃| 老司机影院毛片| av一本久久久久| 国产亚洲av片在线观看秒播厂 | 婷婷六月久久综合丁香| 毛片女人毛片| 我的老师免费观看完整版| 亚洲av.av天堂| 麻豆精品久久久久久蜜桃| 夫妻性生交免费视频一级片| 国产一区二区亚洲精品在线观看| 女的被弄到高潮叫床怎么办| 国产高潮美女av| 亚洲国产精品成人综合色| 亚洲成人久久爱视频| 免费看光身美女| 一边亲一边摸免费视频| 蜜桃久久精品国产亚洲av| 男女边摸边吃奶| 毛片一级片免费看久久久久| 中文字幕av在线有码专区| 一边亲一边摸免费视频| 久久久色成人| 大香蕉久久网| av一本久久久久| 久久久a久久爽久久v久久| 亚洲国产精品sss在线观看| 久久久久网色| 久久久久久久久久久免费av| 久久久精品94久久精品| 亚洲精品一区蜜桃| 2022亚洲国产成人精品| 国产亚洲精品久久久com| 看免费成人av毛片| 国产一区二区亚洲精品在线观看| 久久久久久久久久久丰满| 亚洲欧美精品专区久久| 久久久久精品性色| 精品国产一区二区三区久久久樱花 | av网站免费在线观看视频 | 狂野欧美白嫩少妇大欣赏| 国产午夜福利久久久久久| 中文天堂在线官网| 精品久久久久久久久亚洲| 国产精品嫩草影院av在线观看| 色综合色国产| 最近视频中文字幕2019在线8| 大又大粗又爽又黄少妇毛片口| 97人妻精品一区二区三区麻豆| 能在线免费观看的黄片| 寂寞人妻少妇视频99o| 老司机影院成人| 亚洲精品视频女| 干丝袜人妻中文字幕| 神马国产精品三级电影在线观看| 欧美丝袜亚洲另类| 午夜福利高清视频| 久99久视频精品免费| 永久网站在线| 一区二区三区免费毛片| 日本熟妇午夜| 国产高清国产精品国产三级 | 亚洲av中文av极速乱| 精品欧美国产一区二区三| 国产老妇女一区| 丝袜美腿在线中文| 人妻一区二区av| 亚洲aⅴ乱码一区二区在线播放| 国产在视频线在精品| 亚洲精品,欧美精品| 国产精品熟女久久久久浪| 国产精品.久久久| 免费看美女性在线毛片视频| 亚洲av中文字字幕乱码综合| 亚洲自偷自拍三级| 极品少妇高潮喷水抽搐| 又大又黄又爽视频免费| 国产v大片淫在线免费观看| 亚洲真实伦在线观看| 能在线免费观看的黄片| 日本-黄色视频高清免费观看| 午夜福利在线观看吧| 色综合色国产| 91在线精品国自产拍蜜月| 国产亚洲精品av在线| 久久久久久国产a免费观看| 一级爰片在线观看| 亚洲av男天堂| 亚洲不卡免费看| 哪个播放器可以免费观看大片| 国产69精品久久久久777片| 亚洲av中文av极速乱| 成人av在线播放网站| 日产精品乱码卡一卡2卡三| 国产v大片淫在线免费观看| 日日撸夜夜添| 午夜免费观看性视频| 一夜夜www| 国产精品久久久久久久电影| 在线免费观看不下载黄p国产| 久久精品夜夜夜夜夜久久蜜豆| 天堂网av新在线| 亚洲国产精品成人综合色| 国产精品蜜桃在线观看| 国产av在哪里看| 麻豆精品久久久久久蜜桃| 国产av在哪里看| 色视频www国产| 亚洲成人精品中文字幕电影| 日本黄色片子视频| 欧美日韩在线观看h| 国产一级毛片在线| 欧美日韩精品成人综合77777| 亚洲av不卡在线观看| 久久精品久久久久久噜噜老黄| 成年av动漫网址| 色综合亚洲欧美另类图片| 天堂影院成人在线观看| 丰满人妻一区二区三区视频av| 日本一本二区三区精品| 国产一区亚洲一区在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲内射少妇av| 欧美变态另类bdsm刘玥|