• <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)
    欧美国产精品va在线观看不卡| 久久久水蜜桃国产精品网| 老熟妇仑乱视频hdxx| 999久久久精品免费观看国产| 精品一区二区三区视频在线观看免费 | 久久ye,这里只有精品| 国产av精品麻豆| 国产在线观看jvid| 久久久国产欧美日韩av| 亚洲第一av免费看| 久久狼人影院| 国产av一区二区精品久久| 人妻一区二区av| 日韩三级视频一区二区三区| 岛国毛片在线播放| 俄罗斯特黄特色一大片| 免费女性裸体啪啪无遮挡网站| 久久久久久久精品吃奶| 久久天躁狠狠躁夜夜2o2o| 国产成人影院久久av| 国产麻豆69| 欧美成人免费av一区二区三区 | 成人特级黄色片久久久久久久 | 99精品欧美一区二区三区四区| 亚洲熟妇熟女久久| 亚洲精品一二三| 久久婷婷成人综合色麻豆| 手机成人av网站| 又紧又爽又黄一区二区| 飞空精品影院首页| 精品亚洲成a人片在线观看| av网站免费在线观看视频| a在线观看视频网站| 久久香蕉激情| 精品欧美一区二区三区在线| 黄色怎么调成土黄色| 欧美国产精品va在线观看不卡| 久久国产精品影院| 精品少妇久久久久久888优播| 黄色 视频免费看| 精品少妇内射三级| 超碰成人久久| 久久国产精品影院| 亚洲专区中文字幕在线| 日韩视频在线欧美| 人妻一区二区av| 国产成人精品在线电影| 久久人妻熟女aⅴ| 中文字幕人妻丝袜制服| 99国产精品一区二区蜜桃av | 视频在线观看一区二区三区| 女人精品久久久久毛片| bbb黄色大片| 在线播放国产精品三级| 国产精品免费视频内射| 久久亚洲精品不卡| 嫁个100分男人电影在线观看| av一本久久久久| 国产成+人综合+亚洲专区| 热99国产精品久久久久久7| 国产91精品成人一区二区三区 | 丁香欧美五月| 巨乳人妻的诱惑在线观看| 一个人免费看片子| 日本五十路高清| 一区二区日韩欧美中文字幕| 大香蕉久久网| 欧美精品一区二区大全| 国产熟女午夜一区二区三区| 两人在一起打扑克的视频| 国产精品一区二区在线不卡| 国产精品免费大片| 黄色毛片三级朝国网站| 成年人午夜在线观看视频| 午夜视频精品福利| 精品视频人人做人人爽| av有码第一页| a级毛片在线看网站| 午夜老司机福利片| 高清av免费在线| 啦啦啦视频在线资源免费观看| 国产在视频线精品| 在线观看www视频免费| 免费黄频网站在线观看国产| 大香蕉久久网| 久久性视频一级片| 亚洲精品国产一区二区精华液| 黄色毛片三级朝国网站| 精品福利观看| 91字幕亚洲| 国产欧美日韩一区二区三区在线| av在线播放免费不卡| 91精品三级在线观看| 露出奶头的视频| 十八禁人妻一区二区| 久久久久久亚洲精品国产蜜桃av| 久久影院123| 日韩免费高清中文字幕av| 精品国产亚洲在线| 男女午夜视频在线观看| 精品国产一区二区三区四区第35| 一边摸一边抽搐一进一出视频| 久久人人97超碰香蕉20202| 大香蕉久久网| 夜夜爽天天搞| 女警被强在线播放| 欧美国产精品一级二级三级| 午夜日韩欧美国产| 热re99久久精品国产66热6| 国产又色又爽无遮挡免费看| 久久精品人人爽人人爽视色| 久久性视频一级片| 久久人妻av系列| 亚洲五月色婷婷综合| 欧美亚洲日本最大视频资源| 欧美变态另类bdsm刘玥| 一个人免费在线观看的高清视频| 叶爱在线成人免费视频播放| 欧美日韩视频精品一区| 成人国产av品久久久| 男男h啪啪无遮挡| 一级黄色大片毛片| 免费黄频网站在线观看国产| 亚洲精品美女久久av网站| 国产av精品麻豆| 欧美精品人与动牲交sv欧美| 香蕉丝袜av| 精品久久久久久久毛片微露脸| 母亲3免费完整高清在线观看| 一区二区av电影网| 日本av手机在线免费观看| 精品国产亚洲在线| 久久中文字幕一级| 久久人妻av系列| 午夜视频精品福利| 亚洲精品久久成人aⅴ小说| 国产精品二区激情视频| 99国产精品免费福利视频| 国产深夜福利视频在线观看| av线在线观看网站| 蜜桃在线观看..| 狠狠婷婷综合久久久久久88av| 亚洲色图 男人天堂 中文字幕| 精品欧美一区二区三区在线| 桃花免费在线播放| 国产黄色免费在线视频| 免费观看a级毛片全部| 夫妻午夜视频| 丝袜美足系列| 每晚都被弄得嗷嗷叫到高潮| 精品人妻在线不人妻| 一个人免费在线观看的高清视频| 在线播放国产精品三级| 老熟妇仑乱视频hdxx| 午夜福利,免费看| 久久久久久久精品吃奶| 亚洲天堂av无毛| 18禁黄网站禁片午夜丰满| 人妻久久中文字幕网| 午夜91福利影院| 精品一区二区三区视频在线观看免费 | 考比视频在线观看| 少妇裸体淫交视频免费看高清 | 高清毛片免费观看视频网站 | 国产成人啪精品午夜网站| 男男h啪啪无遮挡| 亚洲欧美激情在线| 欧美在线黄色| 日韩视频一区二区在线观看| 麻豆国产av国片精品| 欧美久久黑人一区二区| 国产精品1区2区在线观看. | 另类精品久久| 亚洲国产毛片av蜜桃av| 国产av一区二区精品久久| 国产精品一区二区在线观看99| 99在线人妻在线中文字幕 | 国产成人影院久久av| 可以免费在线观看a视频的电影网站| 性少妇av在线| 久久人妻av系列| 国产成人精品久久二区二区免费| 一本久久精品| 色播在线永久视频| 亚洲自偷自拍图片 自拍| 每晚都被弄得嗷嗷叫到高潮| 成年人午夜在线观看视频| 757午夜福利合集在线观看| 成人av一区二区三区在线看| 日韩精品免费视频一区二区三区| 国产免费现黄频在线看| 精品国产乱码久久久久久小说| 大片电影免费在线观看免费| 亚洲欧美日韩高清在线视频 | 一级a爱视频在线免费观看| 国产不卡av网站在线观看| 欧美日韩视频精品一区| 免费一级毛片在线播放高清视频 | 久久国产精品男人的天堂亚洲| 欧美日本中文国产一区发布| 亚洲,欧美精品.| 国产精品偷伦视频观看了| 狂野欧美激情性xxxx| 久久人人爽av亚洲精品天堂| 亚洲中文字幕日韩| 午夜福利乱码中文字幕| 最近最新免费中文字幕在线| 欧美另类亚洲清纯唯美| 成年人午夜在线观看视频| 嫩草影视91久久| 国产三级黄色录像| 99精品久久久久人妻精品| 精品一区二区三区av网在线观看 | kizo精华| 日韩大码丰满熟妇| 老司机午夜福利在线观看视频 | 一级片免费观看大全| 大香蕉久久网| 国产有黄有色有爽视频| 中国美女看黄片| 成人手机av| 成年女人毛片免费观看观看9 | 免费av中文字幕在线| svipshipincom国产片| 久久99热这里只频精品6学生| 久久99热这里只频精品6学生| 亚洲熟女毛片儿| 久久亚洲精品不卡| 国产亚洲精品久久久久5区| 国产欧美日韩一区二区精品| 欧美日韩亚洲国产一区二区在线观看 | 亚洲伊人久久精品综合| 国产成+人综合+亚洲专区| 大码成人一级视频| 亚洲中文日韩欧美视频| 精品人妻1区二区| 天天添夜夜摸| 在线av久久热| 亚洲国产欧美网| 老司机深夜福利视频在线观看| 国产高清视频在线播放一区| 国产真人三级小视频在线观看| 久久人妻福利社区极品人妻图片| 麻豆国产av国片精品| 欧美日韩国产mv在线观看视频| 男女无遮挡免费网站观看| 久久亚洲真实| 9色porny在线观看| 91字幕亚洲| 免费高清在线观看日韩| 热99re8久久精品国产| 人人妻人人添人人爽欧美一区卜| 一二三四在线观看免费中文在| 国产又爽黄色视频| 岛国毛片在线播放| 欧美日韩亚洲高清精品| 午夜成年电影在线免费观看| 亚洲精品国产区一区二| 99精品欧美一区二区三区四区| 国产精品美女特级片免费视频播放器 | 久久国产精品人妻蜜桃| 国产视频一区二区在线看| 亚洲精品一二三| 国产色视频综合| 岛国在线观看网站| 中文字幕人妻丝袜制服| 黑人操中国人逼视频| 亚洲av美国av| 久久久精品94久久精品| 午夜福利欧美成人| 日本wwww免费看| 麻豆乱淫一区二区| 免费观看av网站的网址| 亚洲国产欧美网| 亚洲成av片中文字幕在线观看| 欧美日韩亚洲国产一区二区在线观看 | 两性午夜刺激爽爽歪歪视频在线观看 | 九色亚洲精品在线播放| 欧美日韩亚洲综合一区二区三区_| 亚洲精品久久午夜乱码| 精品乱码久久久久久99久播| 色综合欧美亚洲国产小说| 狠狠精品人妻久久久久久综合| 日韩人妻精品一区2区三区| 久久精品国产综合久久久| 99国产精品免费福利视频| www.精华液| 免费看a级黄色片| 欧美乱妇无乱码| 一级a爱视频在线免费观看| 99国产综合亚洲精品| 国产99久久九九免费精品| 久久毛片免费看一区二区三区| 桃花免费在线播放| 老司机靠b影院| 这个男人来自地球电影免费观看| 午夜激情av网站| 国产精品av久久久久免费| 亚洲 欧美一区二区三区| 国产一卡二卡三卡精品| 国产高清videossex| 日本欧美视频一区| 久久久久久亚洲精品国产蜜桃av| 国产男靠女视频免费网站| 成人黄色视频免费在线看| 水蜜桃什么品种好| 国产福利在线免费观看视频| 超碰成人久久| 性色av乱码一区二区三区2| 啦啦啦视频在线资源免费观看| 欧美午夜高清在线| 多毛熟女@视频| 亚洲第一av免费看| av又黄又爽大尺度在线免费看| 99久久人妻综合| 亚洲色图 男人天堂 中文字幕| 一级黄色大片毛片| 亚洲一区二区三区欧美精品| 精品福利观看| 交换朋友夫妻互换小说| 下体分泌物呈黄色| 黄频高清免费视频| 我要看黄色一级片免费的| 国产精品九九99| 两个人看的免费小视频| 美女高潮喷水抽搐中文字幕| av天堂久久9| 日韩欧美一区视频在线观看| 国产亚洲午夜精品一区二区久久| 久热这里只有精品99| 狂野欧美激情性xxxx| 97在线人人人人妻| 夜夜夜夜夜久久久久| 日韩免费av在线播放| 波多野结衣一区麻豆| 欧美成人午夜精品| 国产黄色免费在线视频| 欧美成人午夜精品| 黄片大片在线免费观看| 日韩中文字幕欧美一区二区| 国产精品偷伦视频观看了| 久久国产精品影院| 十八禁人妻一区二区| av免费在线观看网站| 久久久久久久精品吃奶| 水蜜桃什么品种好| 一边摸一边抽搐一进一小说 | 在线观看免费视频网站a站| 国产单亲对白刺激| 日韩视频在线欧美| www.精华液| av欧美777| 欧美日韩亚洲国产一区二区在线观看 | 久久天堂一区二区三区四区| tube8黄色片| 99riav亚洲国产免费| 亚洲成人手机| 久久精品国产亚洲av香蕉五月 | 天天躁狠狠躁夜夜躁狠狠躁| 精品人妻熟女毛片av久久网站| 日本av免费视频播放| 午夜福利影视在线免费观看| 后天国语完整版免费观看| 99精品在免费线老司机午夜| 久久久国产欧美日韩av| av片东京热男人的天堂| 欧美av亚洲av综合av国产av| 成人三级做爰电影| h视频一区二区三区| 色综合欧美亚洲国产小说| 国产一区有黄有色的免费视频| 日韩大片免费观看网站| 巨乳人妻的诱惑在线观看| 亚洲国产欧美在线一区| 精品一品国产午夜福利视频| √禁漫天堂资源中文www| 成年动漫av网址| 国产激情久久老熟女| 精品一品国产午夜福利视频| 国产片内射在线| 久久毛片免费看一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 久久精品aⅴ一区二区三区四区| 国产精品免费大片| 日韩欧美一区二区三区在线观看 | 色视频在线一区二区三区| 中文字幕人妻丝袜一区二区| 看免费av毛片| 黄色视频不卡| 国产成人精品久久二区二区免费| 欧美日韩中文字幕国产精品一区二区三区 | 99国产精品99久久久久| 露出奶头的视频| 午夜免费鲁丝| 午夜福利一区二区在线看| 日韩中文字幕视频在线看片| 亚洲欧美色中文字幕在线| 亚洲免费av在线视频| 日韩成人在线观看一区二区三区| 麻豆国产av国片精品| 亚洲午夜精品一区,二区,三区| 国产精品av久久久久免费| 久久精品亚洲精品国产色婷小说| 国产午夜精品久久久久久| 日本vs欧美在线观看视频| 极品少妇高潮喷水抽搐| 亚洲综合色网址| 亚洲国产看品久久| 国产成人啪精品午夜网站| 国产免费av片在线观看野外av| 久久久国产一区二区| 天天影视国产精品| 国产男女内射视频| 亚洲精品一二三| 国产在线精品亚洲第一网站| 国产一区二区激情短视频| 中文字幕高清在线视频| 性色av乱码一区二区三区2| 女性被躁到高潮视频| 人成视频在线观看免费观看| 女同久久另类99精品国产91| 免费观看a级毛片全部| 999久久久国产精品视频| 国产精品秋霞免费鲁丝片| 新久久久久国产一级毛片| 操美女的视频在线观看| 美女午夜性视频免费| 午夜久久久在线观看| 成年动漫av网址| 正在播放国产对白刺激| 欧美日本中文国产一区发布| 咕卡用的链子| 精品国产乱码久久久久久男人| 99精品在免费线老司机午夜| 99精国产麻豆久久婷婷| 又黄又粗又硬又大视频| 欧美在线一区亚洲| 动漫黄色视频在线观看| 午夜老司机福利片| 巨乳人妻的诱惑在线观看| 午夜成年电影在线免费观看| 日本精品一区二区三区蜜桃| 亚洲欧美日韩高清在线视频 | 中文字幕精品免费在线观看视频| 亚洲国产欧美一区二区综合| 亚洲精品中文字幕在线视频| 啦啦啦免费观看视频1| 99热网站在线观看| 啦啦啦免费观看视频1| 一级毛片电影观看| 91精品三级在线观看| 一区二区三区国产精品乱码| 亚洲精品中文字幕一二三四区 | 国产熟女午夜一区二区三区| 极品少妇高潮喷水抽搐| 法律面前人人平等表现在哪些方面| 久久精品熟女亚洲av麻豆精品| 成人影院久久| 五月开心婷婷网| 日韩免费高清中文字幕av| 亚洲黑人精品在线| 欧美日韩黄片免| 如日韩欧美国产精品一区二区三区| 999久久久国产精品视频| 女人精品久久久久毛片| 考比视频在线观看| 啦啦啦中文免费视频观看日本| 黑丝袜美女国产一区| www日本在线高清视频| 国产黄色免费在线视频| 999久久久国产精品视频| 亚洲成国产人片在线观看| 国产高清激情床上av| 精品一区二区三区av网在线观看 | 国产精品一区二区免费欧美| 国产精品一区二区在线观看99| 我要看黄色一级片免费的| 母亲3免费完整高清在线观看| 在线观看www视频免费| 国产成人一区二区三区免费视频网站| 99re在线观看精品视频| 两个人看的免费小视频| 天天添夜夜摸| 国产精品1区2区在线观看. | 可以免费在线观看a视频的电影网站| 人妻 亚洲 视频| 日本av手机在线免费观看| 99精国产麻豆久久婷婷| 国产精品免费视频内射| 最黄视频免费看| 国产野战对白在线观看| 欧美精品av麻豆av| 成人亚洲精品一区在线观看| 手机成人av网站| 欧美精品一区二区免费开放| 国产日韩欧美亚洲二区| 久久久久精品国产欧美久久久| 精品少妇黑人巨大在线播放| 国产亚洲精品一区二区www | 国产片内射在线| 国产欧美亚洲国产| 啦啦啦中文免费视频观看日本| 男女床上黄色一级片免费看| 啦啦啦视频在线资源免费观看| 国产亚洲精品久久久久5区| 又紧又爽又黄一区二区| 欧美人与性动交α欧美精品济南到| 国产日韩欧美亚洲二区| 欧美成人免费av一区二区三区 | 女人高潮潮喷娇喘18禁视频| 成人国语在线视频| 一个人免费看片子| 自拍欧美九色日韩亚洲蝌蚪91| 国产伦理片在线播放av一区| 亚洲国产欧美网| 国产精品久久久久久人妻精品电影 | 国产精品av久久久久免费| 九色亚洲精品在线播放| 精品国产乱码久久久久久小说| 亚洲第一av免费看| 色视频在线一区二区三区| 国产精品av久久久久免费| 在线观看66精品国产| 精品国产乱码久久久久久小说| 国产又爽黄色视频| 我的亚洲天堂| 久久人妻熟女aⅴ| 黄网站色视频无遮挡免费观看| 亚洲国产av新网站| 国产免费av片在线观看野外av| 丁香六月欧美| 亚洲av日韩在线播放| 妹子高潮喷水视频| 女人精品久久久久毛片| 欧美老熟妇乱子伦牲交| 一区二区三区乱码不卡18| 亚洲一区二区三区欧美精品| av线在线观看网站| 黑人猛操日本美女一级片| 性高湖久久久久久久久免费观看| a级毛片黄视频| 99香蕉大伊视频| 99re在线观看精品视频| av国产精品久久久久影院| 狠狠狠狠99中文字幕| 婷婷成人精品国产| 国产黄频视频在线观看| 午夜福利视频在线观看免费| 国产在线精品亚洲第一网站| 91av网站免费观看| 超碰97精品在线观看| 18禁裸乳无遮挡动漫免费视频| 免费在线观看完整版高清| 国产在视频线精品| 欧美成狂野欧美在线观看| 免费看a级黄色片| 久久影院123| 亚洲九九香蕉| 精品一区二区三区视频在线观看免费 | 大型黄色视频在线免费观看| √禁漫天堂资源中文www| 久久精品人人爽人人爽视色| 亚洲专区国产一区二区| 国产精品美女特级片免费视频播放器 | 两个人免费观看高清视频| 亚洲综合色网址| 黑人巨大精品欧美一区二区蜜桃| 亚洲一码二码三码区别大吗| 欧美日韩国产mv在线观看视频| 精品久久久久久久毛片微露脸| 看免费av毛片| 俄罗斯特黄特色一大片| 777久久人妻少妇嫩草av网站| 中文欧美无线码| 亚洲五月色婷婷综合| 最黄视频免费看| 午夜免费成人在线视频| 99久久国产精品久久久| 久久婷婷成人综合色麻豆| 视频在线观看一区二区三区| 老熟妇仑乱视频hdxx| 国产精品亚洲av一区麻豆| 99re在线观看精品视频| 午夜老司机福利片| 午夜福利欧美成人| 日本wwww免费看| 国产日韩欧美亚洲二区| 黄色片一级片一级黄色片| 无遮挡黄片免费观看| 精品午夜福利视频在线观看一区 | 久久中文字幕人妻熟女| 免费日韩欧美在线观看| 一级,二级,三级黄色视频| 国产精品久久久久久精品古装| 一进一出抽搐动态| 亚洲av成人不卡在线观看播放网| 成人国语在线视频| 人妻 亚洲 视频| 91老司机精品| 精品少妇黑人巨大在线播放| 亚洲av片天天在线观看| 国产福利在线免费观看视频| 侵犯人妻中文字幕一二三四区| 美女主播在线视频| 一个人免费看片子| av一本久久久久| 国产午夜精品久久久久久| av电影中文网址| 最近最新免费中文字幕在线| 婷婷成人精品国产| 日本vs欧美在线观看视频| 久久精品熟女亚洲av麻豆精品| 久久人妻av系列| 午夜福利免费观看在线| 欧美激情久久久久久爽电影 |