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

    中軌道Walker導航星座在軌備份方案優(yōu)化設計

    2021-03-16 09:16:02王許煜胡敏趙玉龍張學陽李玖陽
    中國空間科學技術 2021年1期
    關鍵詞:軌位增量星座

    王許煜,胡敏,趙玉龍,張學陽,李玖陽

    航天工程大學,北京 101400

    導航星座為了滿足系統(tǒng)對可用性、連續(xù)性和完整性的嚴格要求,通常會在空間部署若干顆備份衛(wèi)星。若星座中有衛(wèi)星失效,則通過在軌備份衛(wèi)星變軌對其進行替換。如GPS(global positioning system)星座將備份衛(wèi)星部署在故障概率最大的衛(wèi)星附近形成“衛(wèi)星對”[1],以便其中的一顆衛(wèi)星發(fā)生故障時,備份衛(wèi)星能夠在較短時間內(nèi)通過軌道機動來實現(xiàn)對故障衛(wèi)星的快速替換[2],從而降低對用戶的影響。在軌備份方案的優(yōu)化設計關系到整個導航星座運行期間的服務性能,同時還涉及到空間軌道資源的維護利用等問題[3],因此必須對在軌備份方案進行合理設計。

    目前,與導航星座在軌備份方案優(yōu)化設計相關的研究較少。文獻[3]基于系統(tǒng)服務性能提升,從PDOP(position dilution of precision)值、地面可見星數(shù)和地面對新增衛(wèi)星可視情況3個方面對比分析了北斗衛(wèi)星導航(區(qū)域)系統(tǒng)中GEO和IGSO備份星在不同軌位時,對現(xiàn)有系統(tǒng)服務性能的影響。文獻[4]基于5GEO/5IGSO/4MEO區(qū)域?qū)Ш叫l(wèi)星系統(tǒng),在僅有2顆備份衛(wèi)星的條件下利用星座的幾何關系給出了同軌道面?zhèn)浞莺彤愜壍烂鎮(zhèn)浞輹r,備份星IGSO的最優(yōu)相位分布,從而有效地保證了系統(tǒng)的連續(xù)性和可用性。上述文獻只對GEO和IGSO備份星軌位進行了分析,而沒有考慮MEO備份星軌位的設計,同時由于GEO和IGSO軌道面內(nèi)衛(wèi)星數(shù)量較少,且備份星軌位相對固定,因此也都未涉及到在軌備份星替換的優(yōu)化設計。

    本文以中軌道Walker導航星座為研究對象,首先針對星座運行期間在軌備份星與工作衛(wèi)星存在共同提供服務的情況,通過建立軌位優(yōu)化模型,以PDOP值和可見衛(wèi)星數(shù)為目標函數(shù),利用NSGA-Ⅱ(non-dominated sorting genetic algorithm)算法分析在軌備份星在不同軌位下對系統(tǒng)服務性能的影響。其次,基于在軌備份星軌位的優(yōu)化結(jié)果,建立在軌備份星替換的軌道機動模型,分別以替換所需的速度增量最少和替換時間最少為優(yōu)化目標進行在軌備份星替換方案的優(yōu)化設計,并分析比較兩種替換方案。

    1 在軌備份星軌位優(yōu)化設計

    1.1 軌位優(yōu)化模型

    在軌備份星的軌位直接決定了備份星對系統(tǒng)服務性能的增強效果,與星座優(yōu)化設計一樣,在軌備份星軌位的設計同樣屬于多目標優(yōu)化問題。本文以中軌道Walker導航星座為分析對象,該星座由24顆衛(wèi)星組成,星座參數(shù)為24/3/1,軌道高度為21 528 km,傾角為55°,同時,模型中每個軌道面上分別部署2顆在軌備份星。

    (1) 優(yōu)化變量

    在軌備份星軌位優(yōu)化設計中,模型優(yōu)化變量為每個軌道面上在軌備份星的軌位fi,j,i∈(1,2,3),j∈(1,2),其中i為軌道面編號,j為在軌備份星編號。如f1,2表示第1個軌道面上的第2個在軌備份星的軌位。

    (2) 目標函數(shù)

    在導航星座中,定位精度是其性能評估的重要指標,而該指標除了受到各偽距測量值的影響,還與星座的幾何構(gòu)型有關[5],通過計算星座位置精度衰減因子(PDOP)值可以對構(gòu)型進行量化評估[6]。因此,本文選取PDOP值作為備份星軌位優(yōu)化的評價指標之一,其計算如下[7]。

    假設用戶坐標為(X0,Y0,Z0),此時在用戶本地坐標系中,滿足最小觀測仰角α的N顆衛(wèi)星坐標可表示為:

    ri=[Xi,Yi,Zi],i=1,2,…,N

    (1)

    則相應的系數(shù)矩陣H為:

    (2)

    記矩陣Q為:

    (3)

    最終用戶坐標的PDOP值為:

    (4)

    同時,導航星座除了定位精度要求外,還有其它的性能要求,本文選取可見衛(wèi)星數(shù)作為另一個備份星軌位優(yōu)化的評價指標??梢娦l(wèi)星數(shù)是指在一定的仰角范圍內(nèi),用戶能接收到衛(wèi)星導航信號的衛(wèi)星數(shù)目,其值的大小不僅決定精度因子計算結(jié)果的準確性,還與星座對地面的覆蓋特性密切相關,因此同樣是衡量導航星座性能優(yōu)劣的重要指標。

    為了更全面地評價導航星座服務區(qū)域的性能,本文采取網(wǎng)格分析法,對導航星座的服務區(qū)域進行網(wǎng)格劃分,并統(tǒng)計在每一時刻所有網(wǎng)格點的PDOP值和可見衛(wèi)星數(shù),從而獲得仿真時間內(nèi)每個網(wǎng)格點的平均PDOP值和平均可見衛(wèi)星數(shù),最終對所有網(wǎng)格點的數(shù)值求平均,得到服務區(qū)域內(nèi)PDOP值和可見衛(wèi)星數(shù)的平均值,分別用FP(X)和FM(X)表示。綜上可得,導航星座在軌備份星軌位優(yōu)化設計的數(shù)學模型為:

    (5)

    (3) 約束條件

    在軌備份星軌位優(yōu)化設計中,模型的約束條件為在軌備份星的軌位fi,j∈[0°,360°],同時,由于工作衛(wèi)星的存在,在軌備份星的軌位不能取值為同一軌道面中已有工作衛(wèi)星的軌位。

    1.2 基于NSGA-Ⅱ算法的多目標優(yōu)化設計

    (1) 軌位優(yōu)化變量的編碼方式

    在軌備份星軌位優(yōu)化變量屬于連續(xù)型變量,要利用NSGA-Ⅱ算法對其進行求解,必須對優(yōu)化變量進行編碼。變量編碼有多種編碼方式,其中應用最廣泛的為二進編碼方式,然而當模型優(yōu)化變量取值范圍較大時,為了保證優(yōu)化變量具有較高的精度,需要增大二進制字符串位數(shù),這會導致染色體基因過長,從而降低搜索效率,影響算法的收斂速度。因此,本文采用浮點數(shù)編碼方法,染色體中每個基因代表一個變量,該編碼方法能夠有效降低染色體編碼長度,提高運算效率,同時能確保變量具有較高的精度。對每個變量進行編碼后,將它們串聯(lián)成一個染色體,從而完成對一個個體的編碼,如圖1所示。

    圖1 染色體編碼示意Fig.1 Chromosome coding diagram

    (2) 基于Pareto占優(yōu)的選擇機制

    與單目標優(yōu)化問題不同的是,由于在求解過程中難以獲得一個能使所有目標函數(shù)都達到最優(yōu)值的理想解,當進一步優(yōu)化一個目標時通常會引起其他目標的劣化,因此多目標優(yōu)化問題的最優(yōu)解是一組最優(yōu)解集。已有大量學者對多目標優(yōu)化問題的求解展開了研究并提出相應的求解思路,其中應用較為廣泛的是基于Pareto 占優(yōu)思想的優(yōu)化方法[8]。

    假設多目標優(yōu)化是求解目標函數(shù)的最小值,若可行解x1所對應的任意一個目標函數(shù)值fk(x1)都小于等于可行解x2對應的目標函數(shù)值fk(x2),并且存在一個目標函數(shù)值fl(x1)小于x2對應的目標函數(shù)值fl(x2),則稱可行解x1相比于x2占優(yōu)或x1支配x2。

    (6)

    如果可行解x不被其他任一解所支配,則該可行解稱為非支配解,經(jīng)過算法迭代進化,可以得到所有非支配解,并組成非支配解集,稱為Pareto前沿[9]。在NSGA-Ⅱ算法中采取Pareto占優(yōu)的選擇機制來獲得Pareto前沿,其公式為[10]:

    (7)

    式中:μi,g為試驗向量,Xi,g為目標向量。

    (3) NSGA-Ⅱ算法的優(yōu)化流程

    非支配排序遺傳算法(NSGA)是Srinivas和Deb于1994年提出的基于Pareto占優(yōu)選擇的多目標優(yōu)化算法[11]。在此基礎上,Deb等人于2002年進一步提出了NSGA的改進算法NSGA-Ⅱ[12],該算法較上一代算法的優(yōu)勢主要體現(xiàn)在NSGA-Ⅱ算法采取快速非支配排序法,從而降低了搜索非支配解的復雜度,并定義擁擠距離來計算各層中點之間的擁擠度,代替了NSGA算法中共享函數(shù)的使用,避免了人為確定共享參數(shù)對解空間分布的影響,保持了種群的多樣性,同時還引入精英策略,通過對解空間進行分層確定非支配序列,保留了種群中解的優(yōu)越性,也擴大了采樣空間[13]。NSGA-Ⅱ算法主要包括初始化、選擇、交叉和變異等步驟,其具體流程如圖2所示。

    圖2 NSGA-Ⅱ算法流程框圖Fig.2 NSGA-Ⅱ algorithm flow chart

    1.3 仿真分析

    依據(jù)上述模型及優(yōu)化算法對在軌備份星軌位進行優(yōu)化設計,中軌道Walker導航星座的服務區(qū)域為全球區(qū)域,并按5°×5°的經(jīng)緯線對其進行網(wǎng)格劃分,仿真時間為一個星座回歸周期,用戶最小觀測仰角α為5°,數(shù)據(jù)統(tǒng)計步長為600 s。優(yōu)化算法的初始參數(shù)為:種群個數(shù)N為50,最大進化代數(shù)G為50,交叉因子為1,變異因子為0.16。在軌備份星軌位優(yōu)化結(jié)果如圖3所示。

    圖3 在軌備份星軌位優(yōu)化結(jié)果Fig.3 Optimization results of in-orbit backup satellite orbital position

    從圖3可得,算法迭代獲得的優(yōu)化解集中可見衛(wèi)星數(shù)對應的區(qū)間為[10.330 3,10.330 6],PDOP值對應的區(qū)間為[1.665,1.695]。同時,對于優(yōu)化解集中不同的解,其可見衛(wèi)星數(shù)的變化范圍較小,而且在這小范圍內(nèi),不同解所對應的PDOP值呈現(xiàn)出隨著可見衛(wèi)星數(shù)增大而增大的趨勢。最終,從優(yōu)化解集中篩選出非支配解集,即圖3中的Pareto等級1,并按照PDOP值進行升序排列,如表1所示。

    表1 在軌備份星軌位的非支配解集

    從表1中的非支配解集可得,可見衛(wèi)星數(shù)變化范圍較小,因此本文將選擇非支配解集中PDOP值最小的解作為最優(yōu)解,即f1,1=202.9°,f1,2=247.9°,f2,1=263.4°,f2,2=306.6°,f3,1=142.5°,f3,2=97.5°,并對其進行優(yōu)化結(jié)果分析。最優(yōu)解的在軌備份星軌位如圖4所示,圖4(b)中M11表示第1個軌道面上編號為1的衛(wèi)星。

    圖4 在軌備份星軌位示意Fig.4 Schematic diagram of in-orbit backup satellite orbital position

    針對最優(yōu)解的在軌備份星軌位,對全球區(qū)域內(nèi)中軌道Walker導航星座與加入在軌備份星后的星座PDOP值和可見衛(wèi)星數(shù)進行對比分析,圖5和圖6所示的是在每個統(tǒng)計時刻下,全球區(qū)域內(nèi)所有網(wǎng)格點的PDOP值和可見衛(wèi)星數(shù)的平均值??梢钥闯?,在星座回歸周期內(nèi),中軌道Walker導航星座在不同時刻下PDOP平均值的最大值為1.888 9,最小值為1.875 9,平均值為1.882 4。加入在軌備份星后星座PDOP平均值的最大值降至1.6797,降低幅度為11.1%;最小值降至1.661 8,降低幅度為11.4%;平均值降至1.669 7,降低幅度為11.3%。同理,中軌道Walker導航星座在不同時刻下可見衛(wèi)星數(shù)平均值的最大值為8.284 5,最小值為8.247 6,平均值為8.264 2。加入在軌備份星后可見衛(wèi)星數(shù)平均值的最大值升至10.352 9,增大幅度為24.9%;最小值升至10.307 5,增大幅度為24.9%;平均值升至10.330 4,增大幅度為25%。由此可得,按照優(yōu)化所得的軌位進行在軌備份星的部署,可以顯著提高導航星座的服務性能。

    圖5 星座PDOP值比較Fig.5 Comparison of PDOP values of constellation

    圖 6 星座可見衛(wèi)星數(shù)比較Fig.6 Comparison of the number of visible satellites in constellation

    2 在軌備份星替換方案優(yōu)化設計

    在軌備份星替換方案優(yōu)化設計問題不僅關系到衛(wèi)星故障后導航星座服務性能的恢復,還涉及到在軌備份星替換故障衛(wèi)星過程中的軌道機動問題。目前,針對軌道機動問題通常以軌道機動所需的能量和機動所持續(xù)的時間為評價指標[14]。由于衛(wèi)星燃料的消耗會縮短衛(wèi)星的工作壽命,而單顆衛(wèi)星的壽命會極大地影響整個星座的運行周期,因此,為了最大限度延長衛(wèi)星的工作壽命,保證星座的預期運行周期,應盡量減少在軌備份星燃料的消耗。同時,考慮到衛(wèi)星發(fā)生故障后會導致導航星座性能的下降,從而造成相應的損失,所以,應使在軌備份星替換時間盡可能短,以確保星座服務性能盡快恢復。

    本文基于在軌備份星軌位優(yōu)化結(jié)果,考慮在軌備份星替換故障衛(wèi)星時軌道機動所需的速度增量以及替換時間,對在軌備份星替換方案進行優(yōu)化設計。

    2.1 軌道機動模型

    根據(jù)衛(wèi)星相位調(diào)整過程中是否涉及軌道面的變化,衛(wèi)星相位的調(diào)整方式可以分為同軌道面相位調(diào)整和異軌道面相位調(diào)整[15]。由于異面軌道調(diào)整所消耗的燃料較大,因此在軌備份星通常只對同一軌道面內(nèi)的故障衛(wèi)星進行替換。根據(jù)故障衛(wèi)星相位和在軌備份星相位的關系,備份星同軌道面的替換又可以分為相位超前和相位滯后兩種情況。

    (1) 相位超前

    相位超前時,在軌備份星沿著運行方向到故障衛(wèi)星之間的地心角θ∈[180°,360°)。此時,在軌備份星可以通過抬升軌道高度,在初始位置施加一次沖量進入高軌過渡橢圓軌道,該過渡軌道的半長軸a須滿足:

    a0

    (8)

    式中:a0為原始軌道半長軸。

    當在軌備份星在過渡軌道上運行若干圈后,故障衛(wèi)星位置正好運行到在軌備份星的初始位置,此時,再對在軌備份星施加一次沖量,使其從過渡軌道進入原始軌道,即可完成對故障衛(wèi)星的替換,如圖7所示。

    (2) 相位滯后

    相位滯后時,在軌備份星沿著運行方向到故障衛(wèi)星之間的地心角θ∈(0°,180°)。此時,在軌備份星可以通過降低軌道高度,在初始位置施加一次沖量進入低軌過渡橢圓軌道,該過渡軌道的半長軸a須滿足:

    (9)

    式中:R為地球半徑。

    圖7 相位超前替換示意Fig.7 Schematic diagram of phase ahead replacement

    當在軌備份星在過渡軌道上運行若干圈后,故障衛(wèi)星正好運行到在軌備份星的初始位置,此時,再對在軌備份星施加一次沖量,使其從過渡軌道進入原始軌道,即可完成對故障衛(wèi)星的替換,如圖8所示。

    圖 8 相位滯后替換示意Fig.8 Schematic diagram of phase lag replacement

    依據(jù)上述的軌道機動模型,可以獲得在軌備份星替換故障衛(wèi)星的數(shù)學模型,在軌備份星在原始軌道上的速度為v1,在過渡軌道上變軌點處的速度為v2,由于經(jīng)過2次變軌,因此替換所需的總能量ΔV為2倍的速度增量。

    (10)

    (11)

    ΔV=2|v2-v1|

    (12)

    式中:μ為地球引力常數(shù)。

    假設在軌備份星替換過程中,在軌備份星在過渡軌道上運行的圈數(shù)為n,故障衛(wèi)星在原始軌道上運行的圈數(shù)為m,若不到1圈則為0。在軌備份星完成替換所需的時間為ΔT,則有:

    (13)

    (14)

    2.2 以速度增量最少為優(yōu)化目標

    以在軌備份星替換故障衛(wèi)星時所需的速度增量最少為優(yōu)化目標,分析各在軌備份星對同一軌道面上的8顆工作衛(wèi)星的替換結(jié)果。由于中軌道Walker導航星座中每個軌道面分別部署了2顆在軌備份星,因此,比較2顆在軌備份星的替換結(jié)果,可以得到每個軌道面內(nèi)替換不同故障衛(wèi)星時所需的最少速度增量。模型目標函數(shù)為:

    F=min(ΔVx,y,z),x∈(1,2,3),

    y∈(1,2),z∈(1,2,…8)

    (15)

    式中:ΔVx,y,z為第x個軌道面中編號為y的在軌備份星替換編號為z的故障衛(wèi)星所需的速度增量。

    模型約束條件為:在軌備份星相位超前時采取高軌變相,相位滯后時采取低軌變相,同時,針對實際情況,假設在軌備份星在過渡軌道上運行的圈數(shù)n和故障衛(wèi)星位置在原始軌道上運行的圈數(shù)m都小于等于10。

    最終以速度增量最少為優(yōu)化目標的優(yōu)化設計結(jié)果如圖9所示,圖9(a)為每個軌道面上替換不同衛(wèi)星時所需的最少速度增量,圖9(b)為該替換方案所對應的替換時間。從圖中可以看出,速度增量的最大值為0.106 2 km/s,最小值為0.015 km/s。雖然替換所需的速度增量較小,但是替換所需的時間較長,最短的替換時間需要126.4 h,最長的替換時間則需要134.5 h。

    圖9 以速度增量最少為優(yōu)化目標的優(yōu)化設計結(jié)果Fig.9 Optimization results with the minimum increment of velocity as the optimization objective

    2.3 以替換時間最少為優(yōu)化目標

    以在軌備份星替換故障衛(wèi)星時所需的替換時間最少為優(yōu)化目標,再次對在軌備份星替換軌道面內(nèi)的8顆工作衛(wèi)星的優(yōu)化結(jié)果進行分析。根據(jù)比較每個軌道面中2顆在軌備份星的替換結(jié)果,從而得到每個軌道面內(nèi)替換不同故障衛(wèi)星時所需的最少替換時間。模型目標函數(shù)為:

    F=min(ΔTx,y,z),x∈(1,2,3),

    y∈(1,2),z∈(1,2,…8)

    (16)

    式中:ΔTx,y,z為第x個軌道面中編號為y的在軌備份星替換編號為z的故障衛(wèi)星所需的替換時間。

    模型約束條件為:在軌備份星相位超前時采取高軌變相,相位滯后時采取低軌變相,同時,在軌備份星在過渡軌道上運行的圈數(shù)n和故障衛(wèi)星在原始軌道上運行的圈數(shù)m都小于等于10。

    最終以替換時間最少為優(yōu)化目標的優(yōu)化設計結(jié)果如圖10所示,圖10(a)為每個軌道面上替換不同衛(wèi)星時所需的最少替換時間,圖10(b)為該替換方案所對應的替換速度增量。從圖中可以看出,以替換時間最少為優(yōu)化目標時,替換衛(wèi)星所需的最短時間為8.8 h,最長替換時間為18.5 h,極大地縮短了替換所需的時間。同時,替換所需的速度增量雖然比以速度增量最小為優(yōu)化目標時的大,但速度增量最大值為1.16 km/s,最小值為0.14 km/s,仍處于較低水平。

    圖10 以替換時間最少為優(yōu)化目標的優(yōu)化設計結(jié)果Fig.10 Optimization results with the minimum replacement time as the optimization objective

    2.4 結(jié)果分析

    基于上述對兩種替換方案的分析可得,以速度增量最少為優(yōu)化目標的軌道機動方案,雖然所需的能量較少,但是替換所需的時間較長,適用于主要考慮節(jié)省能量的星座;而對于需要快速恢復服務性能的導航星座而言,該替換方案并不十分合適。而以替換時間最少為優(yōu)化目標的軌道機動方案,其最長的替換時間和最短的替換時間分別降低了86.25%和93.04%,因此該方案能夠在很大程度上減少替換時間;同時,速度增量較以速度增量最小為優(yōu)化目標的替換方案有所增加,備份星單次機動的速度增量最大值為0.58 km/s,最小值為0.07 km/s,而對于衛(wèi)星軌道機動而言,該值仍保持在較低水平。綜合兩種替換方案,本文選用以替換時間最少為優(yōu)化目標的替換方案作為在軌備份星替換的最終方案。

    3 結(jié)束語

    本文對中軌道Walker導航星座在軌備份方案進行了優(yōu)化設計,得到以下結(jié)論:

    1)NSGA-Ⅱ多目標優(yōu)化算法能夠有效實現(xiàn)在軌備份星軌位最優(yōu)解的快速求解,最終在軌備份星軌位的最優(yōu)解為:f1,1=202.9°,f1,2=247.9°,f2,1=263.4°,f2,2=306.6°,f3,1=142.5°,f3,2=97.5°。

    2)按照本文所得的軌位設計方案進行在軌備份星的部署,可以顯著提高導航星座的服務性能,在回歸周期內(nèi),導航星座在全球服務區(qū)域上的PDOP值和可見衛(wèi)星數(shù)的平均值分別降低了11.3%和提高了25%。

    3)通過分析比較以速度增量最少和以替換時間最少為優(yōu)化目標的替換方案,可得以替換時間最少為優(yōu)化目標的替換方案能夠在短時間內(nèi)以較低的能量完成衛(wèi)星的替換,更符合導航星座的需求。

    猜你喜歡
    軌位增量星座
    提質(zhì)和增量之間的“辯證”
    當代陜西(2022年6期)2022-04-19 12:12:22
    “價增量減”型應用題點撥
    ABS要出售15個軌位
    太空探索(2018年8期)2018-11-29 14:50:57
    ABS要出售15個軌位自己用不完
    星座
    健康女性(2017年3期)2017-04-27 21:34:31
    12星座之我愛洗澡
    星座
    都市麗人(2016年12期)2016-12-14 09:54:15
    星座
    都市麗人(2016年7期)2016-08-03 10:06:05
    基于均衡增量近鄰查詢的位置隱私保護方法
    電信科學(2016年9期)2016-06-15 20:27:25
    德州儀器(TI)發(fā)布了一對32位增量-累加模數(shù)轉(zhuǎn)換器(ADC):ADS1262和ADS126
    2018国产大陆天天弄谢| 最近的中文字幕免费完整| 国产精品无大码| 婷婷色综合大香蕉| 人妻夜夜爽99麻豆av| 边亲边吃奶的免费视频| 五月开心婷婷网| 制服丝袜香蕉在线| 这个男人来自地球电影免费观看 | 亚洲中文av在线| 国产精品久久久久久久久免| 干丝袜人妻中文字幕| 在线免费观看不下载黄p国产| 黑人高潮一二区| 99久久中文字幕三级久久日本| 国产淫语在线视频| 亚洲av成人精品一区久久| 亚洲图色成人| 亚洲图色成人| 久久久精品免费免费高清| 国产欧美日韩精品一区二区| 国产 一区 欧美 日韩| 免费人成在线观看视频色| 熟妇人妻不卡中文字幕| 中文字幕人妻熟人妻熟丝袜美| 建设人人有责人人尽责人人享有的 | 中国国产av一级| 在线观看国产h片| 国产男人的电影天堂91| 国产精品无大码| 日日撸夜夜添| 国产精品一区二区在线不卡| 啦啦啦中文免费视频观看日本| 午夜老司机福利剧场| 久久久久久久精品精品| 黄色怎么调成土黄色| av又黄又爽大尺度在线免费看| 男女国产视频网站| 日本色播在线视频| 18禁裸乳无遮挡动漫免费视频| 黄色日韩在线| 永久免费av网站大全| 老司机影院毛片| 身体一侧抽搐| h日本视频在线播放| 肉色欧美久久久久久久蜜桃| av天堂中文字幕网| 午夜日本视频在线| 日本av免费视频播放| 男女免费视频国产| 99视频精品全部免费 在线| 一区二区三区精品91| 婷婷色综合大香蕉| 亚洲第一区二区三区不卡| 国产免费一级a男人的天堂| 亚洲,欧美,日韩| 国产熟女欧美一区二区| 中文在线观看免费www的网站| 免费人成在线观看视频色| 99国产精品免费福利视频| 九草在线视频观看| 亚洲中文av在线| 婷婷色麻豆天堂久久| 男人舔奶头视频| 久久精品国产自在天天线| 熟女电影av网| 国产精品成人在线| 亚洲精品日韩av片在线观看| 黄片无遮挡物在线观看| 三级国产精品欧美在线观看| 欧美日韩视频精品一区| 视频区图区小说| 欧美日韩在线观看h| 精品国产三级普通话版| 精品久久久久久久久亚洲| 久久人人爽人人爽人人片va| 国内揄拍国产精品人妻在线| 亚洲欧美一区二区三区国产| 国产视频首页在线观看| 亚洲国产精品一区三区| 在线观看三级黄色| 97在线视频观看| 精品午夜福利在线看| 黄色怎么调成土黄色| 国产男女内射视频| 黑人猛操日本美女一级片| 伦理电影免费视频| 爱豆传媒免费全集在线观看| 色哟哟·www| 大片电影免费在线观看免费| 看非洲黑人一级黄片| 高清av免费在线| 精品一区二区免费观看| 精品午夜福利在线看| 一边亲一边摸免费视频| 久久国产亚洲av麻豆专区| 六月丁香七月| 国产精品三级大全| 美女福利国产在线 | 婷婷色综合www| 国产老妇伦熟女老妇高清| 偷拍熟女少妇极品色| 99久国产av精品国产电影| 97精品久久久久久久久久精品| 欧美精品亚洲一区二区| 97精品久久久久久久久久精品| 99热6这里只有精品| 一个人免费看片子| 久热这里只有精品99| 国产在线视频一区二区| 亚洲无线观看免费| 交换朋友夫妻互换小说| 哪个播放器可以免费观看大片| 日韩中字成人| 激情五月婷婷亚洲| 菩萨蛮人人尽说江南好唐韦庄| 久久人人爽人人爽人人片va| 国产综合精华液| 国产成人aa在线观看| 少妇人妻精品综合一区二区| 中文字幕制服av| 美女视频免费永久观看网站| 国产伦精品一区二区三区四那| 亚洲久久久国产精品| 国产精品人妻久久久久久| 国产白丝娇喘喷水9色精品| 欧美激情国产日韩精品一区| 日韩一本色道免费dvd| 色吧在线观看| 国产精品三级大全| 天天躁日日操中文字幕| 国产视频内射| 久久这里有精品视频免费| 国产大屁股一区二区在线视频| 99re6热这里在线精品视频| 午夜福利在线观看免费完整高清在| 亚洲真实伦在线观看| 国产精品偷伦视频观看了| 啦啦啦在线观看免费高清www| 97在线视频观看| 深爱激情五月婷婷| 精品一区在线观看国产| 噜噜噜噜噜久久久久久91| 99久久精品国产国产毛片| 乱系列少妇在线播放| 欧美成人精品欧美一级黄| 99热全是精品| 青春草视频在线免费观看| 成人午夜精彩视频在线观看| 色综合色国产| 一个人看视频在线观看www免费| 午夜福利在线观看免费完整高清在| 美女脱内裤让男人舔精品视频| 精品久久久久久久久亚洲| 18禁裸乳无遮挡动漫免费视频| 大片免费播放器 马上看| 日韩免费高清中文字幕av| 成年av动漫网址| 久久精品夜色国产| 又爽又黄a免费视频| 人妻一区二区av| 人人妻人人看人人澡| tube8黄色片| 久久久久久久亚洲中文字幕| 亚洲av成人精品一区久久| 欧美变态另类bdsm刘玥| 你懂的网址亚洲精品在线观看| 黄片无遮挡物在线观看| 欧美精品国产亚洲| 最近中文字幕高清免费大全6| 一级片'在线观看视频| 婷婷色综合www| 九九爱精品视频在线观看| 国产日韩欧美在线精品| 精品亚洲成国产av| 久久久久久久久久久免费av| 日韩精品有码人妻一区| 高清日韩中文字幕在线| 91久久精品国产一区二区成人| 亚洲欧美一区二区三区国产| 亚洲婷婷狠狠爱综合网| 久久精品国产亚洲av天美| 欧美激情极品国产一区二区三区 | 成人18禁高潮啪啪吃奶动态图 | 国产精品国产av在线观看| 国产欧美亚洲国产| 3wmmmm亚洲av在线观看| 久久久久久久久大av| 黄色怎么调成土黄色| 性色avwww在线观看| 中文字幕人妻熟人妻熟丝袜美| 人人妻人人澡人人爽人人夜夜| 亚洲精品视频女| 成人特级av手机在线观看| 在线观看国产h片| 下体分泌物呈黄色| 免费黄网站久久成人精品| 精品久久久噜噜| 大话2 男鬼变身卡| 亚洲一级一片aⅴ在线观看| 久久久久久久久大av| 美女内射精品一级片tv| 观看美女的网站| 免费在线观看成人毛片| 国产成人午夜福利电影在线观看| 高清黄色对白视频在线免费看 | 欧美激情国产日韩精品一区| 黄片无遮挡物在线观看| 亚洲精品国产色婷婷电影| av国产精品久久久久影院| 久久午夜福利片| 欧美高清性xxxxhd video| 熟女av电影| 亚洲欧美日韩无卡精品| 国产黄片视频在线免费观看| 99久国产av精品国产电影| 成人美女网站在线观看视频| 高清视频免费观看一区二区| 久久久久久久久久人人人人人人| 美女国产视频在线观看| 欧美少妇被猛烈插入视频| 免费黄色在线免费观看| 免费少妇av软件| 成人免费观看视频高清| 免费观看a级毛片全部| 丰满乱子伦码专区| 韩国av在线不卡| 五月天丁香电影| 中文在线观看免费www的网站| 精品熟女少妇av免费看| 久久久久久九九精品二区国产| 97超视频在线观看视频| 国模一区二区三区四区视频| 蜜桃久久精品国产亚洲av| 乱码一卡2卡4卡精品| 99视频精品全部免费 在线| 国产精品国产三级专区第一集| 97在线视频观看| 能在线免费看毛片的网站| 99久久人妻综合| 精品少妇黑人巨大在线播放| 最近中文字幕2019免费版| 成人国产av品久久久| 久久毛片免费看一区二区三区| 国产伦理片在线播放av一区| 黄片无遮挡物在线观看| 97超视频在线观看视频| 久久 成人 亚洲| 欧美高清性xxxxhd video| 国产免费又黄又爽又色| 国产一区二区三区综合在线观看 | 18+在线观看网站| 国产免费一级a男人的天堂| 菩萨蛮人人尽说江南好唐韦庄| 久久韩国三级中文字幕| 欧美成人午夜免费资源| 亚洲国产日韩一区二区| 久久久久久九九精品二区国产| 日日撸夜夜添| 久久97久久精品| 一区二区三区免费毛片| 欧美丝袜亚洲另类| 亚洲三级黄色毛片| 三级经典国产精品| 日韩中字成人| 亚洲色图综合在线观看| 国产色爽女视频免费观看| 黄色配什么色好看| 免费观看a级毛片全部| 成人影院久久| 国产免费视频播放在线视频| 欧美极品一区二区三区四区| 亚洲经典国产精华液单| 亚洲精品日韩av片在线观看| 亚洲av电影在线观看一区二区三区| 建设人人有责人人尽责人人享有的 | 极品少妇高潮喷水抽搐| 99视频精品全部免费 在线| 日韩中文字幕视频在线看片 | 亚洲精品久久久久久婷婷小说| 日韩成人av中文字幕在线观看| 极品教师在线视频| 岛国毛片在线播放| 精品久久久久久电影网| 久久人人爽av亚洲精品天堂 | 欧美亚洲 丝袜 人妻 在线| 久久人人爽av亚洲精品天堂 | 亚洲av国产av综合av卡| 国产亚洲最大av| 青春草国产在线视频| 国产精品不卡视频一区二区| 精品国产乱码久久久久久小说| 久久亚洲国产成人精品v| 99热国产这里只有精品6| 简卡轻食公司| www.av在线官网国产| 一级a做视频免费观看| 欧美极品一区二区三区四区| 精品国产露脸久久av麻豆| 免费不卡的大黄色大毛片视频在线观看| 欧美3d第一页| 国产精品一区二区三区四区免费观看| 免费播放大片免费观看视频在线观看| 春色校园在线视频观看| 国产伦精品一区二区三区四那| 少妇人妻一区二区三区视频| 男男h啪啪无遮挡| 中文在线观看免费www的网站| 大香蕉久久网| 亚洲av男天堂| 亚洲精品国产av成人精品| 免费黄频网站在线观看国产| 久久人妻熟女aⅴ| 国产精品爽爽va在线观看网站| 夜夜骑夜夜射夜夜干| 亚洲av成人精品一区久久| 在线观看免费日韩欧美大片 | 日本与韩国留学比较| 国产高清有码在线观看视频| 久久精品国产a三级三级三级| h日本视频在线播放| 我的女老师完整版在线观看| 亚洲精品成人av观看孕妇| 亚洲精华国产精华液的使用体验| 三级国产精品欧美在线观看| 欧美性感艳星| 亚洲人成网站高清观看| 全区人妻精品视频| 国产亚洲91精品色在线| 久久人人爽人人片av| 久久亚洲国产成人精品v| 蜜臀久久99精品久久宅男| 国产成人a区在线观看| 午夜免费鲁丝| 91久久精品国产一区二区成人| 99国产精品免费福利视频| 久久久久久久亚洲中文字幕| 国产在线男女| 欧美人与善性xxx| 亚洲成人av在线免费| 亚洲精品日韩在线中文字幕| 国产黄频视频在线观看| 亚州av有码| 有码 亚洲区| 18禁在线播放成人免费| 六月丁香七月| 亚洲av不卡在线观看| 国产成人精品福利久久| 国产精品福利在线免费观看| 久久影院123| 在线免费观看不下载黄p国产| 久久久欧美国产精品| 美女脱内裤让男人舔精品视频| 国产真实伦视频高清在线观看| 青青草视频在线视频观看| 啦啦啦啦在线视频资源| 91狼人影院| av国产免费在线观看| 欧美一区二区亚洲| 成年免费大片在线观看| 亚洲精品一二三| 亚洲av中文av极速乱| 成人18禁高潮啪啪吃奶动态图 | 日本黄色日本黄色录像| 观看av在线不卡| 91精品伊人久久大香线蕉| 国产大屁股一区二区在线视频| 日本-黄色视频高清免费观看| 精品午夜福利在线看| 我要看黄色一级片免费的| 观看免费一级毛片| 最近中文字幕2019免费版| 日韩中文字幕视频在线看片 | 久久国产精品男人的天堂亚洲 | 黄色欧美视频在线观看| 国产精品.久久久| 18禁裸乳无遮挡动漫免费视频| 一二三四中文在线观看免费高清| 国产亚洲一区二区精品| 久久av网站| 免费人妻精品一区二区三区视频| 在线观看国产h片| 日韩 亚洲 欧美在线| 国产精品国产av在线观看| av在线观看视频网站免费| 日韩强制内射视频| 91精品国产九色| 日日摸夜夜添夜夜添av毛片| 国产亚洲5aaaaa淫片| 综合色丁香网| 久久久久国产精品人妻一区二区| 国产欧美日韩精品一区二区| 大话2 男鬼变身卡| 一级毛片 在线播放| 毛片一级片免费看久久久久| 欧美 日韩 精品 国产| 看十八女毛片水多多多| 日本与韩国留学比较| 国产综合精华液| 国产乱来视频区| 精品一品国产午夜福利视频| 91狼人影院| 国产 一区精品| 亚洲国产欧美人成| 97精品久久久久久久久久精品| 国产伦在线观看视频一区| 日韩三级伦理在线观看| 男人爽女人下面视频在线观看| kizo精华| 在线观看一区二区三区| 97精品久久久久久久久久精品| 国产黄片美女视频| 97在线视频观看| 最近2019中文字幕mv第一页| 丝袜脚勾引网站| 午夜精品国产一区二区电影| 高清在线视频一区二区三区| 黄色一级大片看看| 最黄视频免费看| 丰满迷人的少妇在线观看| 少妇人妻精品综合一区二区| 在线观看免费视频网站a站| 欧美精品人与动牲交sv欧美| 久久久久久伊人网av| 亚洲国产最新在线播放| 亚洲图色成人| 国产爱豆传媒在线观看| 免费大片18禁| 麻豆乱淫一区二区| 亚洲最大成人中文| 亚洲av.av天堂| 国产成人aa在线观看| 97超碰精品成人国产| 最黄视频免费看| 天堂8中文在线网| 亚洲精品国产av蜜桃| 成人免费观看视频高清| 日日撸夜夜添| 欧美精品亚洲一区二区| 人人妻人人看人人澡| 日本av免费视频播放| 毛片一级片免费看久久久久| 中文在线观看免费www的网站| 高清黄色对白视频在线免费看 | 亚洲精品日韩在线中文字幕| 在线观看美女被高潮喷水网站| 香蕉精品网在线| 亚洲成人av在线免费| 各种免费的搞黄视频| 免费黄色在线免费观看| 亚洲欧美日韩东京热| 精品一区二区免费观看| 国语对白做爰xxxⅹ性视频网站| 国产精品国产av在线观看| 国产免费又黄又爽又色| 亚洲性久久影院| 国产精品成人在线| 少妇裸体淫交视频免费看高清| 极品教师在线视频| 亚洲电影在线观看av| 欧美少妇被猛烈插入视频| 国产无遮挡羞羞视频在线观看| 成人免费观看视频高清| 啦啦啦在线观看免费高清www| 美女视频免费永久观看网站| 亚洲av欧美aⅴ国产| 我要看黄色一级片免费的| 欧美日韩一区二区视频在线观看视频在线| 一二三四中文在线观看免费高清| 尾随美女入室| 国产男女内射视频| av在线app专区| 亚洲欧洲国产日韩| 国产无遮挡羞羞视频在线观看| 在线观看免费日韩欧美大片 | 免费人妻精品一区二区三区视频| 日韩 亚洲 欧美在线| 午夜福利在线观看免费完整高清在| 91精品国产国语对白视频| 国产欧美亚洲国产| 18禁在线播放成人免费| 久久av网站| 国产高清国产精品国产三级 | 国国产精品蜜臀av免费| 大又大粗又爽又黄少妇毛片口| 一本色道久久久久久精品综合| 97在线人人人人妻| 国产精品国产三级专区第一集| 水蜜桃什么品种好| 夜夜看夜夜爽夜夜摸| 能在线免费看毛片的网站| 国产成人freesex在线| 91在线精品国自产拍蜜月| 国产又色又爽无遮挡免| 一边亲一边摸免费视频| 久久久a久久爽久久v久久| 日韩精品有码人妻一区| 成人二区视频| 国产永久视频网站| 高清欧美精品videossex| 国产黄片美女视频| 日韩免费高清中文字幕av| 人人妻人人澡人人爽人人夜夜| 男男h啪啪无遮挡| h日本视频在线播放| 一区二区三区四区激情视频| 国产精品99久久久久久久久| 青春草亚洲视频在线观看| 97超碰精品成人国产| 噜噜噜噜噜久久久久久91| 亚洲丝袜综合中文字幕| 大码成人一级视频| 久久久a久久爽久久v久久| 免费播放大片免费观看视频在线观看| 99热全是精品| 国产爱豆传媒在线观看| 婷婷色综合www| 精品熟女少妇av免费看| 亚洲成人手机| 网址你懂的国产日韩在线| 久久久精品94久久精品| 老女人水多毛片| 五月伊人婷婷丁香| 日韩亚洲欧美综合| 亚洲经典国产精华液单| 18禁在线播放成人免费| kizo精华| 777米奇影视久久| 高清毛片免费看| 日韩伦理黄色片| 亚洲av综合色区一区| 久久久亚洲精品成人影院| 夜夜骑夜夜射夜夜干| 日本黄色日本黄色录像| 欧美一区二区亚洲| 欧美变态另类bdsm刘玥| 午夜激情久久久久久久| 国产v大片淫在线免费观看| 自拍偷自拍亚洲精品老妇| 精品一区二区三区视频在线| h日本视频在线播放| 免费观看无遮挡的男女| 国产精品人妻久久久久久| 久久99蜜桃精品久久| 免费观看的影片在线观看| a级一级毛片免费在线观看| 久久精品久久久久久噜噜老黄| 国产人妻一区二区三区在| 高清午夜精品一区二区三区| 免费观看性生交大片5| 男女无遮挡免费网站观看| 九九久久精品国产亚洲av麻豆| 欧美成人精品欧美一级黄| 亚洲av日韩在线播放| 国产免费视频播放在线视频| 亚洲欧美成人精品一区二区| 欧美日韩视频精品一区| 国产乱来视频区| 精品人妻偷拍中文字幕| 国产成人a∨麻豆精品| 久久av网站| 免费av中文字幕在线| 免费少妇av软件| 一区二区av电影网| 亚洲精品456在线播放app| 国内少妇人妻偷人精品xxx网站| 日韩欧美精品免费久久| 国产成人一区二区在线| 日本猛色少妇xxxxx猛交久久| 在线观看人妻少妇| 亚洲国产av新网站| 国产亚洲一区二区精品| 男人舔奶头视频| 男人和女人高潮做爰伦理| 国产高清不卡午夜福利| 嫩草影院新地址| 人妻 亚洲 视频| 久久久久久伊人网av| 又爽又黄a免费视频| 最新中文字幕久久久久| 欧美日韩亚洲高清精品| 亚洲一区二区三区欧美精品| 天天躁日日操中文字幕| 欧美+日韩+精品| 国产精品精品国产色婷婷| 中国美白少妇内射xxxbb| 在线观看三级黄色| 大又大粗又爽又黄少妇毛片口| 一级毛片 在线播放| 亚洲av成人精品一区久久| 亚洲av.av天堂| 国产乱人偷精品视频| 精品久久久久久久久亚洲| 美女高潮的动态| 免费观看a级毛片全部| 晚上一个人看的免费电影| 亚洲欧洲国产日韩| 蜜桃亚洲精品一区二区三区| 观看av在线不卡| 极品教师在线视频| 少妇人妻久久综合中文| 国产精品一二三区在线看| 亚洲精品日韩av片在线观看| 亚洲国产精品国产精品| 日韩大片免费观看网站| 观看免费一级毛片| 日韩中文字幕视频在线看片 | 午夜日本视频在线| 久久久久国产精品人妻一区二区| 黄片无遮挡物在线观看| 99热这里只有是精品在线观看| 日本色播在线视频| 熟妇人妻不卡中文字幕| 日本午夜av视频| av在线app专区|