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

    空間多體繩系編隊動力學(xué)與控制研究綜述

    2023-08-22 04:49:20王長青孫佩杰李愛軍
    宇航學(xué)報 2023年7期
    關(guān)鍵詞:繩系系繩構(gòu)型

    王長青,楊 航,孫佩杰,李愛軍

    (西北工業(yè)大學(xué)自動化學(xué)院,西安 710072)

    0 引 言

    空間多體繩系編隊是由空間雙體繩系系統(tǒng)演變而來,通常由多個衛(wèi)星或航天器通過特殊材質(zhì)系繩連接組成的某種特定構(gòu)形的編隊系統(tǒng)[1-3],如圖1所示??臻g多體繩系編隊系統(tǒng)在進入預(yù)定空間軌道后,會根據(jù)任務(wù)要求展開形成相應(yīng)的編隊構(gòu)型,在編隊系統(tǒng)中,一些大型航天器可由多個結(jié)構(gòu)簡單的衛(wèi)星代替,衛(wèi)星之間可以通過通信來進行信息共享。此外,空間多體繩系編隊系統(tǒng)還可以根據(jù)任務(wù)需求及變化,通過改變系繩長度等方式完成編隊構(gòu)型重組,以實現(xiàn)不同類別任務(wù)的目標。相比于無繩系編隊系統(tǒng),多體繩系編隊系統(tǒng)的主要優(yōu)勢一方面在于這種系統(tǒng)可通過調(diào)整系繩張力來保持或改變隊形,而傳統(tǒng)的編隊系統(tǒng)在隊形控制過程中需要依靠推進器產(chǎn)生推力,會消耗大量燃料;另一方面繩系編隊衛(wèi)星系統(tǒng)可獲得較長的測量基線,并通過收放系繩來改變測量基線。相對于傳統(tǒng)的自由編隊方式,空間多體繩系編隊系統(tǒng)的研究雖然起步較晚,但其巨大的應(yīng)用前景已經(jīng)得到學(xué)術(shù)界以及工程界的廣泛認可,被認為在空間目標捕獲[4]、在軌服務(wù)[5-6]、空間干涉儀[7]、電子監(jiān)視等[8]多方面具有很高的應(yīng)用潛力和價值。例如近些年來世界各國開展的大量在軌繩系衛(wèi)星試驗,包括NASA的SEDS系列、俄羅斯和歐洲聯(lián)合進行的YES2空間試驗以及美國ISAS與日本JAXA合作開展的T-REX等太空試驗[9-11],均驗證了空間繩系系統(tǒng)太空任務(wù)的可行性和應(yīng)用前景。

    圖1 空間多體繩系編隊系統(tǒng)

    空間多體繩系編隊運行于特殊的太空環(huán)境,不僅構(gòu)形種類繁多,還存在較為復(fù)雜的動力學(xué)特性,因此需要對編隊構(gòu)型、動力學(xué)模型、編隊構(gòu)型形成、編隊構(gòu)型穩(wěn)定保持以及星體姿態(tài)機動等方面展開研究。本文介紹了常見的繩系編隊構(gòu)型及其特點,對常見的繩系系統(tǒng)動力學(xué)模型進行了總結(jié);梳理了繩系編隊任務(wù)全流程控制研究現(xiàn)狀,包括繩系編隊構(gòu)型形成控制、構(gòu)型穩(wěn)定保持研究以及剛體姿態(tài)機動控制;分析總結(jié)了現(xiàn)有的地面試驗裝置及方法,展望了未來的發(fā)展趨勢。

    1 空間多體繩系編隊構(gòu)型及動力學(xué)模型

    1.1 繩系編隊構(gòu)型

    空間多體繩系編隊系統(tǒng)因航天器個數(shù)的不同而存在多種構(gòu)型組成方式,在進入到預(yù)定空間軌道之后,多體繩系編隊系統(tǒng)會根據(jù)不同的空間任務(wù)要求產(chǎn)生不同的構(gòu)型。繩系編隊的構(gòu)型一般按照末端體在空間的分布形式或編隊穩(wěn)定飛行方式進行分類。按照空間分布形式可分為三大類:一維構(gòu)型、二維構(gòu)型和三維構(gòu)型。其中一維構(gòu)型即衛(wèi)星均位于一條直線上,屬于此類構(gòu)型的繩系編隊系統(tǒng)稱為直鏈式構(gòu)型[12-13],如圖2所示;二維構(gòu)型即在空間中衛(wèi)星位于同一平面內(nèi),常見的此類構(gòu)型有閉環(huán)三角形構(gòu)型[14]、開環(huán)-輪輻式構(gòu)型[15]、閉環(huán)-輪輻構(gòu)型[16]以及TetraStar構(gòu)型[17]等,如圖3所示;三維構(gòu)型即系統(tǒng)中的衛(wèi)星處于不同平面的構(gòu)型,研究較多的此類構(gòu)型包括四面體構(gòu)型[18]、雙四面體構(gòu)型[19]、金字塔構(gòu)型[20]以及雙金字塔構(gòu)型[21]等,如圖4所示。

    圖2 一維構(gòu)型

    圖3 二維構(gòu)型

    圖4 三維構(gòu)型

    按照繩系編隊的穩(wěn)定飛行形式,系統(tǒng)還可以分為靜態(tài)構(gòu)型和動態(tài)構(gòu)型兩大類[22]。在靜態(tài)構(gòu)型編隊中,整個編隊系統(tǒng)與軌道坐標系保持相對靜止。靜態(tài)構(gòu)型編隊主要以重力梯度、阻力以及電磁力等方式來保持構(gòu)型的穩(wěn)定,例如圖2所示的一維重力梯度穩(wěn)定構(gòu)型、一維阻力穩(wěn)定構(gòu)形。對于動態(tài)構(gòu)型,編隊系統(tǒng)會相對于軌道坐標系進行旋轉(zhuǎn)運動,從而使得系繩拉緊來實現(xiàn)整個編隊系統(tǒng)的穩(wěn)定。例如二維離心力穩(wěn)定構(gòu)型、三維離心力-重力穩(wěn)定構(gòu)型,如圖5所示。

    圖5 多維動態(tài)構(gòu)型

    動態(tài)構(gòu)型編隊為了使系繩在編隊中保持張力,系統(tǒng)必須旋轉(zhuǎn)。在張力的作用下,整個系統(tǒng)會產(chǎn)生剛性作用,即系統(tǒng)可以被視作剛體。相反,如果系繩沒有張力,或者當(dāng)系統(tǒng)在拉緊和松弛狀態(tài)之間變換時,各星體之間會發(fā)生碰撞或纏繞等事故。如圖6所示,對于在圓軌道上軸對稱旋轉(zhuǎn)的繩系編隊系統(tǒng)有三種相對平衡:圓柱形平衡、雙曲線平衡和圓錐形平衡,在這種平衡下,整個編隊系統(tǒng)自旋軸以與軌道運動相同的速度旋轉(zhuǎn)[23]。在這三種平衡狀態(tài)中,圓柱形平衡態(tài)自旋軸垂直于軌道平面,雙曲線平衡態(tài)自旋軸垂直于當(dāng)?shù)卮咕€,編隊運動在地球上的投影是一條直線,因此這兩種平衡態(tài)不適用于繩系編隊對地球進行觀測。而對于圓錐平衡,其自旋軸垂直于軌道切線方向,其運動投影到地球表面將呈橢圓形,可以增大對地球的觀測范圍。

    圖6 相對平衡狀態(tài)

    不同構(gòu)型的繩系編隊具有不同的結(jié)構(gòu)特性,例如系統(tǒng)穩(wěn)定方式、適用任務(wù)場景等。余本嵩等在文獻[24]中對不同構(gòu)型下的繩系編隊系統(tǒng)穩(wěn)定性進行了總結(jié),如表1所示。本文查閱了從2000年至2022年期間的多篇空間多體繩系編隊研究文獻,參考了102篇較為權(quán)威的文獻,給出了不同構(gòu)型繩系編隊的文獻研究數(shù)量對比,其中一維構(gòu)型繩系編隊研究占比38%,二維構(gòu)型研究占比52%,三維構(gòu)型占比10%。在二維構(gòu)型繩系編隊中,三角形構(gòu)型和輪輻式構(gòu)型是研究熱點,占比分別達到了49%和34%。可以看出,目前對于空間多體繩系編隊系統(tǒng)的研究主要集中于一維構(gòu)型和平面式構(gòu)型,而空間構(gòu)型由于結(jié)構(gòu)比較復(fù)雜,模型的復(fù)雜程度較高,因此對于該類構(gòu)型研究較少。

    表1 不同構(gòu)型下的繩系編隊系統(tǒng)穩(wěn)定方式[24]

    1.2 系繩模型

    系繩模型建立是空間繩系系統(tǒng)研究的基礎(chǔ),目前應(yīng)用最廣泛的系繩模型主要包括三種,剛性桿模型、彈性桿模型以及珠點模型。在早期的研究中,研究人員通常將系繩看作剛體,不考慮系繩的彈性和柔性,這種簡易的模型主要用來設(shè)計系統(tǒng)的標稱展開軌跡。隨著研究的不斷深入,學(xué)者們又提出了彈性桿模型,這種模型將系繩看作一種只能在受拉時產(chǎn)生拉力而無法在受壓時產(chǎn)生彈力的“彈簧”,相當(dāng)于只考慮了系繩的縱向振蕩。但實際的系繩除了在縱向會產(chǎn)生振蕩,在橫向也會有振蕩。針對這種情況,研究人員建立了系繩的珠點模型,將系繩看作一個個珠點,珠點之間通過彈性桿連接,這種模型也是現(xiàn)在研究中常用的相對精確的一種系繩模型。

    1.2.1剛性桿模型

    剛性桿模型是一種較為簡單的系繩模型,如圖7所示,該模型忽略了系繩的彈性和柔性,將其視為一根不可拉壓的剛性桿,即系繩剛度趨于無窮大,兩端分別連接著系統(tǒng)中的不同航天器。

    圖7 剛性桿模型

    剛性桿模型在90年代被眾多學(xué)者提出,主要用于雙體繩系系統(tǒng)的標稱展開軌跡設(shè)計[25-26],之后諸多學(xué)者將剛性桿模型拓展至多體繩系編隊模型的建立,通過簡化系繩模型,對繩系編隊飛行的可行性以及動力學(xué)進行分析。Tragesser等將剛性桿模型分別用于任意數(shù)量衛(wèi)星的環(huán)形繩系編隊和雙金字塔型繩系編隊標稱軌跡的設(shè)計[27]。之后Kumar等忽略了系繩質(zhì)量、阻尼效應(yīng)以及橫向振動,通過所提出的約束系統(tǒng)非線性運動控制方程的數(shù)值模擬結(jié)果和平衡分析,驗證了在軌道平面上實現(xiàn)三角形編隊飛行的可行性[28]。以上研究均未考慮平臺姿態(tài),Menon等將系繩建模為剛性桿模型,考慮了星體姿態(tài),研究了編隊中各平臺姿態(tài)的動力學(xué)和控制[29]。由于平動點在深空探測中有巨大的應(yīng)用前景,眾多學(xué)者對平動點軌道上的繩系編隊動力學(xué)特性進行了深入研究。Sanjurjo-Rivo等利用剛性桿模型,研究了繩系編隊系統(tǒng)在日-地平動點處的動力學(xué)分析,并給出了幾個穩(wěn)定系統(tǒng)的控制律,為繩系系統(tǒng)在太空探索提供了一定的理論支撐[30]。Cai等[31]考慮系繩為不可伸展且忽略質(zhì)量的剛性桿,分析了一種輪輻式繩系衛(wèi)星編隊在平動點附近的動力學(xué),并基于希爾近似,提出了一種新的非旋轉(zhuǎn)或旋轉(zhuǎn)繩系編隊的動力學(xué)公式。在此基礎(chǔ)上,Cai等[32-33]又對閉環(huán)三角形繩系編隊動力學(xué)穩(wěn)定性進行了深入研究,分析了軌道振幅、旋轉(zhuǎn)速率以及系繩長度等參數(shù)對系統(tǒng)穩(wěn)定性的影響。為了研究日地平動點繩系衛(wèi)星編隊的重定向控制,Luo等[34]利用剛性桿模型,設(shè)計了用于系統(tǒng)狀態(tài)和控制輸入的標稱軌跡。Huang等[35]研究了一種輪輻式繩系機器人編隊系統(tǒng),忽略了系繩質(zhì)量和彈性等影響,提出了一種系繩張力和主星扭矩的協(xié)同控制策略。之后,又在另一種輪輻式繩系編隊建模中,采用剛性桿假設(shè),研究了系繩的展開和回收問題[36]。Zhai等[37]針對眾多研究中未考慮星體和系繩之間的耦合問題,基于剛性桿模型,考慮了星體自轉(zhuǎn)與系繩展開之間的耦合,研究了軌道平面內(nèi)輪輻式繩系編隊的展開控制問題。Chen等[38]將編隊系統(tǒng)建模為質(zhì)點-剛性桿系統(tǒng),對近地軌道上三角形繩系納米衛(wèi)星系統(tǒng)的展開問題進行了研究。Wang等[39]假設(shè)系繩為剛性桿,針對線性繩系衛(wèi)星編隊系統(tǒng)設(shè)計了分階段展開的方案,利用Lagrange方程建立了編隊系統(tǒng)的標稱模型,通過參數(shù)優(yōu)化得到了繩系衛(wèi)星編隊系統(tǒng)的標稱軌跡。Zhang等[40]基于簡化的剛性桿模型,分析了直鏈式繩系編隊在構(gòu)型轉(zhuǎn)換過程中的非線性動力學(xué)特性。

    1.2.2彈性桿模型

    在實際情況下,系繩會隨著所受外力的變化而發(fā)生形變,進而產(chǎn)生內(nèi)應(yīng)力。因此學(xué)者們又引入了較剛性桿模型更為貼合實際繩系系統(tǒng)的彈性桿模型。如圖8所示,彈性桿模型考慮了系繩的彈性而忽略了其柔性,將系繩視為一根不會彎曲的直桿。

    圖8 彈性桿模型

    Pizarro等[41]考慮系繩彈性并將星體近似為質(zhì)點,研究了開環(huán)和閉環(huán)兩種輪輻式結(jié)構(gòu)繩系編隊的動力學(xué),使用Lagrange方法推導(dǎo)了兩個系統(tǒng)的運動方程。Williams等[42]基于彈性桿模型,提出了一種考慮軌道運動、系繩彈性以及重力梯度效應(yīng)的三角形繩系編隊動力學(xué)模型,該模型被眾多學(xué)者在研究三角形繩系編隊時采用[43-45]。耿凌璐[46]考慮系繩彈性桿模型,利用Newton第二定律建立了三角形繩系編隊在非慣性參考系下的動力學(xué)方程,研究了J2攝動、熱效應(yīng)等因素對編隊穩(wěn)定性的影響。Qi等[47]考慮系繩彈性,利用Lagrange方程得到了一種庫侖雙金字塔繩系編隊的運動方程,探討了該編隊應(yīng)用于地球同步衛(wèi)星配置的可行性。Zhang等[48]也提出了一種金字塔型繩系庫侖編隊,考慮彈性桿模型,并利用Lagrange方程推導(dǎo)了編隊的動力學(xué)模型。Huang等[49]在Williams的研究基礎(chǔ)上,考慮系繩彈性,利用Lagrange方程建立了三角形繩系編隊動力學(xué)模型,并提出了一種基于濾波器的狀態(tài)估計方案,解決了繩系編隊狀態(tài)測量的問題?;谠撃P?Su等[50]進行了一系列三角形編隊動力學(xué)特性以及穩(wěn)定展開控制研究,計算出了編隊自旋穩(wěn)定時的自轉(zhuǎn)角速度范圍,為后續(xù)的編隊控制提供了很好的理論支撐。

    1.2.3珠點模型

    珠點模型將系繩離散為一系列由無質(zhì)量彈簧連接的質(zhì)點,兩端分別與主星和子星相連接,如圖9所示。隨著所選取的離散點的數(shù)量越多,模型與真實系繩系統(tǒng)越接近,但同時也使系統(tǒng)變得復(fù)雜,同時增大了理論研究的難度與仿真運算的時間。

    圖9 珠點模型

    Yan等[51]首次將系繩珠點模型應(yīng)用在直鏈式構(gòu)型繩系編隊中。之后,Tragesser等將這種模型用在了閉環(huán)自旋多體繩系衛(wèi)星編隊和雙四面體型繩系編隊上面,用于實現(xiàn)地球定向觀測任務(wù)[27]。Nakaya等[52-53]也將珠點模型用于三角形繩系衛(wèi)星編隊系統(tǒng)上,并提出了釋放新珠點時克服系統(tǒng)不連續(xù)性的方法。Williams等在文獻[54]中也利用珠點模型對三角形繩系編隊的動力學(xué)進行了分析,并將其拓展至雙金字塔型繩系編隊。為了更加清晰的描述系統(tǒng)的動力學(xué)響應(yīng)進而分析控制系統(tǒng)的性能,在基于該模型對繩系系統(tǒng)進行研究時,研究者們一般采用Galerkin模態(tài)法[55]和有限元法[56]完成方程的離散化。劉壯壯等[57]考慮了系繩的密度非線性、非線性彈力以及空間攝動因素的影響改進了珠點模型,并定量分析了執(zhí)行任務(wù)時系繩的橫縱向振動。余本嵩等[58]在前人研究的基礎(chǔ)上提出了一種時變的珠點模型,隨著系繩的不斷展開,系繩上的珠點也不斷的增加。為了突出系繩質(zhì)量對系統(tǒng)動力學(xué)的影響,Avanzini等[59]利用珠點模型建模獲得的結(jié)果與簡單無質(zhì)量系繩模型進行比較,結(jié)果表明系繩質(zhì)量對封閉構(gòu)型的編隊動力學(xué)產(chǎn)生較大影響。之后,Avanzini等又在文獻[60]中利用珠點模型分析了軌道偏心率對閉環(huán)繩系編隊的影響。Luo等[13]將系繩建模為珠點模型,基于希爾近似,開發(fā)了一種新的柔性多體繩系編隊動力學(xué)模型,該模型可應(yīng)用于任意構(gòu)型的N體繩系衛(wèi)星編隊系統(tǒng)?;谠撃P?Luo等[34]又研究了日地平動點直鏈式繩系編隊的重定向控制。

    通過前人的研究可以看出,雖然剛性桿模型和彈性桿模型無法反映真實情況下的繩系系統(tǒng)狀態(tài),但有利于降低控制器設(shè)計難度并提高計算效率,便于對系統(tǒng)進行一些簡單直接的分析,因此被廣泛應(yīng)用于繩系系統(tǒng)控制方法驗證、標稱軌跡設(shè)計等研究。當(dāng)需要提高模型精度時,則需要同時考慮系繩質(zhì)量、彈性、形變等因素,從而建立更加精確的珠點模型。雖然珠點模型的精確度較高,更加貼切實際任務(wù)場景中的系繩,但是該模型的復(fù)雜度很高,會導(dǎo)致計算耗時很長,不便于對系統(tǒng)進行分析。所以大部分研究者會先使用系繩較為簡單的剛性桿模型、彈性模型來驗證猜想,然后在珠點模型上施加實際控制器得到驗證。

    2 空間多體繩系編隊全流程控制

    空間多體繩系編隊在進入太空后,根據(jù)任務(wù)控制需求可分為三個部分,包括繩系編隊構(gòu)型形成控制、構(gòu)型穩(wěn)定保持以及星體姿態(tài)機動控制。

    2.1 空間多體繩系編隊構(gòu)型形成控制

    編隊構(gòu)型形成是空間多體繩系編隊一切應(yīng)用技術(shù)的基礎(chǔ),包括實現(xiàn)繩系編隊系統(tǒng)的快速穩(wěn)定展開、回收。此外,針對不同的太空任務(wù),繩系編隊還可能通過編隊重構(gòu)形成其他編隊構(gòu)型。為使繩系編隊按照合理的方式完成構(gòu)型形成,達到系繩與航天器之間不發(fā)生纏繞、盡可能的節(jié)省燃料等實際目標,需要對系繩展開方案進行具體的設(shè)計。目前已經(jīng)有許多關(guān)于繩系編隊展開及回收的研究,控制方法主要有張力控制、長度控制以及速度控制等。Nakaya等[52]基于虛擬結(jié)構(gòu)方法,使用預(yù)先確定的角動量分布以及系繩張力分布來計算自旋角速度,從而實現(xiàn)三角形繩系編隊的展開。Williams等[61-62]利用張力控制確定了不同自旋條件下的最佳展開和回收軌跡,并對系統(tǒng)旋轉(zhuǎn)速率和機動時間的影響進行了參數(shù)研究,數(shù)值結(jié)果說明了系繩展開和回收的對稱性。Liu等[63]針對輪輻式構(gòu)型繩系衛(wèi)星編隊開發(fā)了一種純張力控制策略,以抑制展開過程中的系繩振動,而不向子衛(wèi)星施加推力,最終系繩能夠完全穩(wěn)定地展開到所需的輪輻式結(jié)構(gòu)。劉麗麗等[64]研究了輪輻式三體繩系衛(wèi)星編隊面內(nèi)飛行時的最優(yōu)回收控制問題,針對不同回收初值及受擾情況,計算了繩系衛(wèi)星編隊飛行的最優(yōu)控制張力和飛行軌跡??紤]到張力控制的局限性,Chen等[38]提出了一種系繩張力和外部低推力結(jié)合的控制律,在實現(xiàn)穩(wěn)定展開的同時保證了控制律的簡潔性和有效性。Kumar等[65]設(shè)計了一套繩長控制律,使三體直鏈式繩系衛(wèi)星編隊可以沿水平方向逐一穩(wěn)定展開/回收。文獻[66-67]針對日地平動點附近的輪輻式旋轉(zhuǎn)繩系衛(wèi)星編隊,利用數(shù)值仿真研究了編隊重構(gòu)階段的動力學(xué)穩(wěn)定性。蘇飛等[68]針對圓軌道上的開環(huán)輪輻式繩系衛(wèi)星編隊,分別研究了在重力梯度力補償和無補償兩種情況下的自旋展開控制策略,通過速度控制實現(xiàn)了繩系編隊的有效展開,仿真結(jié)果表明兩種展開控制策略能夠保證編隊穩(wěn)定展開。以上針對系繩的展開控制研究均未考慮星體姿態(tài)運動對系繩的影響,于是,Zhai等[37]針對一種開環(huán)輪輻式繩系編隊,考慮到了母星的自旋運動與系繩展開之間的耦合,研究了軌道平面內(nèi)自旋穩(wěn)定繩系編隊的展開控制問題,并且開發(fā)了具有重力梯度補償?shù)目刂撇呗浴|S攀峰等針對三角形繩系編隊展開問題進行了深入研究,在文獻[69]中定量分析了編隊自旋穩(wěn)定時的自轉(zhuǎn)角速度范圍,為之后的穩(wěn)定展開控制提供了理論依據(jù)。由于針對繩系編隊狀態(tài)估計的研究較少,Huang等在文獻[70]中研究了具有狀態(tài)觀測器的繩系編隊的穩(wěn)定展開,針對系統(tǒng)速度不可測的情況,給出了一種高增益觀測器來對速度變量進行估計。之后,Huang等[71-72]又針對在實際空間環(huán)境中繩系編隊可能存在的通信資源受限問題,基于事件觸發(fā)機制設(shè)計了一系列展開控制策略,不僅可以實現(xiàn)編隊的穩(wěn)定展開,還能減少對通信資源的消耗。

    2.2 空間多體繩系編隊構(gòu)型穩(wěn)定保持

    繩系衛(wèi)星編隊系統(tǒng)構(gòu)型穩(wěn)定性是繩系編隊領(lǐng)域的研究重點和熱點,由于太空中重力梯度力和低軌道高度下的空氣阻力的存在,一維線性構(gòu)型繩系編隊能夠滿足重力梯度力穩(wěn)定和阻力穩(wěn)定的條件,并保持線形的構(gòu)型穩(wěn)定。Misra等[73]研究了直鏈式構(gòu)型編隊系統(tǒng)在構(gòu)型形成過程中的四種平衡位置及其穩(wěn)定性,驗證了其中沿當(dāng)?shù)卮咕€的平衡位置是穩(wěn)定的。Avanzini等[59]從Misra的研究結(jié)果出發(fā),對一種由四顆衛(wèi)星連接得到的編隊系統(tǒng)進行了研究,得出除了四顆衛(wèi)星共線的平衡方式外,還有16種平衡解決方案的特殊情況。在線形構(gòu)型的基礎(chǔ)上,又有學(xué)者研究了多顆衛(wèi)星平面排布的平面構(gòu)型方法。Tragesser等在文獻[74]中基于錐形Likins-Pringle相對平衡的軌道設(shè)計,研究了具有任意數(shù)量衛(wèi)星的環(huán)形編隊系統(tǒng)穩(wěn)定性。Kumar等[28]對一種三角形繩系編隊穩(wěn)定性進行了研究,實驗結(jié)果證明要使系統(tǒng)穩(wěn)定則需要編隊旋轉(zhuǎn)速率與軌道角速率的比值不能小于0.58。Pizarro等[75]對Y字型輪輻式繩系編隊和閉環(huán)十字型繩系編隊穩(wěn)定性進行了研究,發(fā)現(xiàn)當(dāng)編隊在軌道平面上旋轉(zhuǎn)時,對于具有四顆星體的輪輻式構(gòu)型是穩(wěn)定的。當(dāng)星體數(shù)量超過4顆后,則需要外部系繩構(gòu)成閉環(huán)式編隊才能保持穩(wěn)定。陳志明等[76-77]根據(jù)力學(xué)方法證明了Thomson平衡下繩系衛(wèi)星編隊系統(tǒng)自旋角速率與軌道角速率的比值在0.61和2.61之間時,系統(tǒng)可以保持穩(wěn)定。Yu等[78]使用Floquet理論得到了使旋轉(zhuǎn)三角形編隊系統(tǒng)留位狀態(tài)穩(wěn)定的初始速度取值范圍。蔡志勤等[79]深入研究了位于日地平動點附近的旋轉(zhuǎn)三角形繩系編隊系統(tǒng)的動力學(xué)與控制問題,重點分析了編隊系統(tǒng)在構(gòu)型形成和保持時的穩(wěn)定性,研究結(jié)果均表明三角形繩系編隊可以通過旋轉(zhuǎn)保持穩(wěn)定。

    在平面構(gòu)型的基礎(chǔ)上可以構(gòu)建更為復(fù)雜空間編隊構(gòu)型,并通過自旋和重力梯度力等使其在繞地飛行中保持穩(wěn)定。Sabatini等[80]對一種雙金字塔構(gòu)型繩系編隊進行了穩(wěn)定性研究,認為該種構(gòu)型保持穩(wěn)定的條件極為苛刻,在真實條件下很難成功。Williams等[54]也對雙金字塔構(gòu)型繩系編隊系統(tǒng)進行了動力學(xué)仿真研究,從其結(jié)構(gòu)特點出發(fā),對在旋轉(zhuǎn)軸與地面垂直時的運動穩(wěn)定性進行了研究。Avanzini等[60]對一種由重力梯度力及其自旋離心力共同作用保持穩(wěn)定的雙四面體繩系編隊的穩(wěn)定性進行了研究,并對其構(gòu)型穩(wěn)定性的保持進行了控制?;陔p金字塔型穩(wěn)定性的研究,黃攀峰等[81]首次提出了雙金字塔繩系衛(wèi)星編隊的狀態(tài)估計方案研究。

    2.3 空間多體繩系編隊姿態(tài)機動控制

    繩系編隊運行在太空中,針對不同的任務(wù)需求,編隊星體往往需要進行一些姿態(tài)機動來達到理想的工作位置。由于編隊各星體通常具有較大的質(zhì)量和尺寸,因此星體姿態(tài)運動對系統(tǒng)產(chǎn)生的影響不可忽視。目前對空間繩系系統(tǒng)姿態(tài)動力學(xué)的研究主要聚焦于雙體繩系系統(tǒng),在很多利用雙體繩系系統(tǒng)進行碎片捕獲、軌道轉(zhuǎn)移等研究的文獻中,均對星體的姿態(tài)動力學(xué)進行了深入分析。然而多體繩系編隊由于結(jié)構(gòu)復(fù)雜,大部分研究均考慮星體為質(zhì)點,對于姿態(tài)動力學(xué)的研究還較少。Lemke等[82]首先探索了通過利用系繩張力提供的扭矩來控制繩系衛(wèi)星姿態(tài)的可能性。Modi等[83]研究了繩系衛(wèi)星編隊剛體平臺的姿態(tài)動力學(xué)和控制,并基于Lyapunov第二種方法提出了一種非常有效的偏移控制算法。Nakaya等[52]將衛(wèi)星視作剛體并考慮其姿態(tài),使用虛擬結(jié)構(gòu)法實現(xiàn)了三角形繩系衛(wèi)星編隊的展開控制。Chung等[84-85]利用分散控制討論了多體繩系編隊的自旋速率以及衛(wèi)星剛體姿態(tài)控制問題,這是第一次提出了用于繩系編隊的無推進劑欠驅(qū)動控制方法,并基于麻省理工學(xué)院SPHERES試驗臺進行了驗證。Chang等[86]研究了一種三體直鏈式繩系編隊,并且考慮了每一個星體的姿態(tài),采用狀態(tài)相關(guān)的Riccati方程(SDRE)控制器調(diào)節(jié)姿態(tài)誤差,并將集中式和分散式方法應(yīng)用于動態(tài)系統(tǒng),仿真驗證了SDRE控制器用于繩系衛(wèi)星姿態(tài)控制具有良好的控制性能。蘇飛等[68]和Zhai等[87]均考慮到母星姿態(tài)與系繩展開之間的耦合關(guān)系,對輻射開環(huán)形繩系編隊穩(wěn)定展開控制進行了研究。Sun等[88]提出了一種任意數(shù)量的直鏈式繩系衛(wèi)星編隊,考慮了母星姿態(tài)運動,利用絕對節(jié)點坐標公式(ANCF)和自然坐標公式(NCF)分別用于精確描述柔性系繩和剛性衛(wèi)星的大變形和大整體運動。Luo等[34]考慮衛(wèi)星的姿態(tài)運動,研究了日地L2平動點處柔性多繩系衛(wèi)星編隊重定向控制,所建立的靈活模型可用于描述具有任意結(jié)構(gòu)的多繩系衛(wèi)星編隊系統(tǒng)。考慮到系繩的擺動問題,賈程等[89]基于母星的姿態(tài)機動,通過調(diào)整母星推力方向?qū)崿F(xiàn)了繩系衛(wèi)星編隊系統(tǒng)軌道跟蹤。Chen等[90]考慮了母星旋轉(zhuǎn)與系繩展開之間的耦合效應(yīng),建立了軌道平面上輪輻式構(gòu)型繩系編隊的動力學(xué)模型。并基于Bellman動態(tài)規(guī)劃方法,考慮了展開過程中系繩張力的物理約束,提出了一種最優(yōu)控制器,通過數(shù)值仿真驗證了所提控制策略的有效性。

    3 繩系系統(tǒng)地面試驗

    由于太空任務(wù)復(fù)雜,成本高昂,因此在進行太空任務(wù)前需要進行充分的研究與試驗,而地面模擬試驗對在軌繩系系統(tǒng)的研究起到非常重要的作用。目前國內(nèi)外已經(jīng)開展了多種空間多體繩系編隊系統(tǒng)的地面試驗研究方案。

    3.1 地面氣浮平臺法

    氣浮平臺法具有建造周期短、成本低、易于實現(xiàn)及維護等優(yōu)點,通過設(shè)計平板止推軸承的大小能夠?qū)崿F(xiàn)高達幾噸的飛行器微重力實驗,且實驗的時間不受限制。通過更換接口部件即可實現(xiàn)重復(fù)利用,可靠性和魯棒性較高,對飛行器結(jié)構(gòu)沒有太多限制,且根據(jù)不同實驗?zāi)康目梢暂^為容易地進行二次設(shè)計和改造。地面氣浮臺作為一種結(jié)構(gòu)簡單、成本較低的試驗裝置,已經(jīng)被諸多學(xué)者應(yīng)用于空間繩系系統(tǒng)的地面試驗驗證,是目前繩系系統(tǒng)地面試驗應(yīng)用最廣泛的一種方法。

    Nakaya等[91]為了研究空間繩系衛(wèi)星編隊問題,提出了實現(xiàn)旋轉(zhuǎn)三角形繩系編隊構(gòu)型形成的方案,如圖10所示。由同時發(fā)射的航天器(步驟a)組成的初始系統(tǒng),通過系繩展開(步驟b)形成穩(wěn)定的編隊構(gòu)型(步驟c)。為了驗證所提出方案的可行性,Nakaya等開發(fā)了一套在地面情況下用于模擬繩系編隊的氣浮平臺系統(tǒng),該系統(tǒng)包含氣浮平臺、三個衛(wèi)星模擬器和一個CCD相機組成的位置測量子系統(tǒng),最終成功實現(xiàn)了三角形繩系衛(wèi)星編隊的穩(wěn)定展開。

    圖10 繩系衛(wèi)星編隊構(gòu)型形成方案[91]

    麻省理工學(xué)院空間系統(tǒng)實驗室(MIT-SSL)開發(fā)了名為SPHERES的多顆納米衛(wèi)星試驗臺,主要用于編隊飛行和自主對接、交會和重構(gòu)算法試驗驗證[92]。在該平臺基礎(chǔ)上,Chung等[93-94]成功進行了一系列直鏈式、三角形繩系編隊系統(tǒng)的旋轉(zhuǎn)展開動力學(xué)試驗,驗證了在未來的空間干涉儀任務(wù)中不依賴推進器的情況下控制繩系編隊的可行性。之后,Chung又在文獻[95]中使用多顆SPHERES微型衛(wèi)星對繩系編隊飛行的一些關(guān)鍵控制策略進行了試驗驗證?;谠撈脚_,Mantellato等[96]開發(fā)了一個用于模擬雙質(zhì)量繩系系統(tǒng)在軌動力學(xué)的軟件模型,并通過SPHERES平臺對軟件模型進行了驗證,為后續(xù)在太空微重力環(huán)境開展更精確的試驗進行了預(yù)研。Iki等[97]通過電動力繩系(Electro-Dynamic Tether,EDT)系統(tǒng)的數(shù)值模擬,利用氣浮臺研究了使用推進器從卷軸式卷筒上展開系繩。推進器用于確保展開長度為數(shù)公里的系繩,通過地面試驗評估了幾個關(guān)鍵的參數(shù),如推力和推進器激活周期,明確了展開所需推進器的要求,該裝置主要用于研究EDT技術(shù)進行軌道轉(zhuǎn)移的主動碎片清除系統(tǒng)。

    為了研究繩系編隊動力學(xué)問題,金棟平等[98]為空間繩系系統(tǒng)研究建立了氣浮平臺設(shè)施,該平臺可以模擬繩系系統(tǒng)所經(jīng)歷的微重力場,通過使用平面空氣軸承試驗臺和科里奧利力,在反饋控制下使用機載推力。之后,該團隊利用動態(tài)相似原理構(gòu)造了等效地面試驗系統(tǒng)[78],該系統(tǒng)由花崗巖試驗臺、三個由系繩連接的衛(wèi)星模擬器和一個立體視覺測量系統(tǒng)(SVMS)組成,如圖11所示?;谠撛囼炏到y(tǒng),該團隊研究了在近地軌道上飛行的三角形繩系衛(wèi)星編隊的旋轉(zhuǎn)穩(wěn)定性,成功實現(xiàn)了直鏈式構(gòu)型和三角形構(gòu)型繩系編隊的穩(wěn)定展開和保持。

    圖11 繩系衛(wèi)星編隊地面氣浮臺試驗[78]

    為了驗證繩系編隊動力學(xué)與初始狀態(tài)之間的關(guān)系,Yu等[99]又在該試驗臺上對長繩系編隊系統(tǒng)(LSTF)動力學(xué)進行了分析,討論了LTSF的三種典型面內(nèi)運動以及面外運動抑制,給出了三種運動形式的擺狀振蕩、旋轉(zhuǎn)運動和不規(guī)則運動。圖12為LSTF旋轉(zhuǎn)運動軌跡仿真,可以看出,地面試驗?zāi)M器軌跡與衛(wèi)星軌跡基本保持一致,驗證了地面等效試驗的準確性。

    圖12 LSTF旋轉(zhuǎn)運動軌跡[99]

    3.2 機構(gòu)測試法

    為了驗證繩系系統(tǒng)釋放、制動和回收等試驗機構(gòu)的可靠性,研究人員提出了機構(gòu)測試法,主要是通過一些特殊設(shè)計的機構(gòu)來測試繩系系統(tǒng)內(nèi)部裝置的特性。這些機構(gòu)相對整個系統(tǒng)而言質(zhì)量較小,受重力影響較低,因此在使用該方法進行相關(guān)模擬實驗和可靠性驗證時,對微重力環(huán)境的要求較低。

    易琳等[100]研究了用于空間捕獲的繩網(wǎng)收口方式,介紹了一種由雙轉(zhuǎn)子電機、雙卷筒和自由導(dǎo)繩套筒組成的自適應(yīng)系繩卷取機構(gòu)。該機構(gòu)能適應(yīng)負載不平衡變化,不會發(fā)生系繩纏繞現(xiàn)象,具有較強的自適應(yīng)性。Schultz等[101]為模擬Tips試驗[102]中系繩內(nèi)部張力變化的情況,采取控制系繩水平位移,使系繩只在豎直方向上移動的結(jié)構(gòu),利用擺錘來提供可變張力,如圖13所示。通過選取合適的映射關(guān)系,分析系繩在試驗中的內(nèi)部張力情況可以研究在軌任務(wù)中系繩的內(nèi)部張力情況。

    圖13 繩系衛(wèi)星水平結(jié)構(gòu)配置地面測試裝置示意圖[101]

    Nakamura等[103]開發(fā)了能夠展開和回收微型衛(wèi)星的微型系繩纏繞機構(gòu)(MTRM)。它由一個卷軸、一個驅(qū)動電機和一個平風(fēng)機構(gòu)組成,包括一根2 km長的系繩,可控制卷軸電機扭矩以避免系繩松弛。此外在該機構(gòu)基礎(chǔ)上搭建了一個名為“繩系動力學(xué)模擬器”(TDSM)的地面試驗系統(tǒng)。利用該模擬器對雙體繩系編隊進行了地面試驗,用于驗證該機構(gòu)對系繩展開和回收的可靠性和有效性。試驗最終成功將系繩展開至2 km。但在回收階段,振動角在最后100 m發(fā)散。

    3.3 失重飛機飛行法

    除了以上提到的試驗裝置外,另一種通過飛機拋物線飛行構(gòu)造微重力環(huán)境的方案也被一些學(xué)者應(yīng)用于繩系系統(tǒng)的試驗驗證。其原理是飛機通過先爬升再迅速向下俯沖進行飛行的過程可以營造時間約幾十秒10-2~10-3g左右的失重環(huán)境,可以用于精度較高的三維空間微重力實驗?zāi)M,但這種方法成本高昂,而且實驗周期較短,目前僅有少數(shù)學(xué)者進行了失重飛機環(huán)境試驗測試。

    Benvenuto等[104]設(shè)計了一種可裝載在失重飛機上的系繩展開試驗裝置,該裝置主要由兩顆立方體衛(wèi)星結(jié)構(gòu)組成,衛(wèi)星結(jié)構(gòu)內(nèi)部附有系繩展開機構(gòu)、高速攝像測量模塊和數(shù)據(jù)存儲模塊。在完成展開控制律的地面試驗驗證后,Benvenuto等將該裝置裝在飛機上,在失重環(huán)境下進行試驗。Zhu等[105]通過空氣軸承轉(zhuǎn)臺和微重力拋物線飛行試驗,對為微重力環(huán)境中的繩系立方體衛(wèi)星展開進行了試驗驗證。繩系立方體衛(wèi)星任務(wù)包括兩顆1U立方體星,100 m長的帶狀系繩。系繩通過拉動展開,兩顆立方體以期望的速度分離。該試驗進行了各種地面測試,以驗證設(shè)計的分離速度、估計內(nèi)部摩擦力和測量展開的系繩長度。最后,通過微重力拋物線飛行試驗驗證了繩系系統(tǒng)展開的工程模型。試驗過程如圖14所示。

    圖14 拋物線飛行下的系繩展開試驗[105]

    該試驗給出了兩個立方體衛(wèi)星之間相對距離的實測時間歷程,如圖15所示??梢钥闯?立方體衛(wèi)星在自由飛行狀態(tài)下的相對分離距離軌跡幾乎是一條線性直線。虛線表示自由飛行狀態(tài)下的相對分離速度,為1.42 m/s,高于地面試驗測得的1.23 m/s的速度。試驗結(jié)果證明了繩系系統(tǒng)在微重力環(huán)境中的內(nèi)摩擦力小于在重力環(huán)境中測量的內(nèi)摩擦力,驗證了所設(shè)計的繩系系統(tǒng)的有效性,為后續(xù)的試驗提供了很好的理論參考價值。

    圖15 衛(wèi)星相對距離和分離速度[105]

    地面試驗是在進行昂貴、耗時的航天試驗之前,驗證繩系系統(tǒng)動力學(xué)和控制理論研究的重要工具。表2比較了幾種地面模擬試驗方法特點??梢钥闯?在已經(jīng)進行的地面仿真試驗研究中,氣浮平臺法和機構(gòu)測試法由于造價低、設(shè)備維護簡單等優(yōu)點,已經(jīng)被諸多學(xué)者用于繩系系統(tǒng)地面試驗驗證,但氣浮平臺只能在系統(tǒng)層面通過極低的摩擦力來模擬失重環(huán)境,無法為系統(tǒng)中的每一個結(jié)構(gòu)進行失重模擬。通過氣浮平臺實現(xiàn)三維空間的復(fù)雜運動的實驗較為困難,主要用于二維平面的飛行器模擬試驗。而機構(gòu)測試法僅適用于一些簡單的內(nèi)部機構(gòu)可靠性驗證,無法進行高精度的模擬試驗。失重飛機試驗雖然可以實現(xiàn)對微重力環(huán)境的高精度模擬,但由于設(shè)備造價及維護成本高昂,試驗條件極為苛刻,且單次實驗的微重力實驗時間較短,對飛行器的各種設(shè)備性能指標無法全面考核,目前已開展的相關(guān)實驗較少。

    表2 各種地面模擬試驗方法對比

    4 結(jié)束語

    空間多體繩系編隊研究的快速發(fā)展將促進太空探索和開發(fā)技術(shù)的進步。本文首先簡要介紹了空間多體繩系編隊的構(gòu)型,并總結(jié)了國內(nèi)外研究學(xué)者對繩系編隊動力學(xué)建模與控制、地面試驗方面的研究成果。可以看出,作為一種新型的空間編隊技術(shù),繩系編隊已經(jīng)吸引了越來越多的關(guān)注,并且已經(jīng)取得了很多研究成果,但仍存在一些不足,未來應(yīng)從這幾個方面開展進一步的研究。

    1) 目前關(guān)于繩系編隊動力學(xué)建模方面,針對諸如J2攝動、熱效應(yīng)、大氣阻力以及太陽光壓等空間環(huán)境外部因素的擾動,已經(jīng)積累了很多研究。但在編隊星體姿態(tài)的研究上,大部分系統(tǒng)建模均將星體考慮為質(zhì)點,忽略了星體姿態(tài)的影響,僅有少量文獻考慮到母星的姿態(tài)運動,但仍忽略了子星姿態(tài)。另外關(guān)于進一步考慮系繩與星體運動之間的耦合、系統(tǒng)姿態(tài)與軌道之間耦合的研究則更加稀少。因此,考慮繩系編隊系統(tǒng)耦合關(guān)系的動力學(xué)建模分析將是未來研究的重點。

    2) 空間多體繩系編隊作為一種節(jié)省能量的編隊構(gòu)成方式,缺乏基于能量分析和對比方面的文獻研究,未來需要更多從能量角度考慮的研究,例如與無繩系衛(wèi)星編隊的能量消耗對比,說明繩系編隊在能量節(jié)省方面的可靠性。

    3) 在現(xiàn)有的空間多體繩系編隊星體姿態(tài)研究中,大部分文獻沒有考慮星體姿態(tài)機動的局限性,針對繩系編隊星體姿態(tài)運動的限制范圍分析,值得更加深入的研究。

    4) 在編隊構(gòu)型研究上,目前大多文獻僅針對某一種特定構(gòu)型繩系編隊進行控制,幾乎沒有考慮繩系編隊在不同構(gòu)型之間進行切換的研究。編隊重構(gòu)控制將是未來繩系編隊研究的難點與重點。

    5) 以往關(guān)于繩系編隊的地面試驗主要集中于結(jié)構(gòu)較為簡單的直鏈式或三角形繩系編隊,沒有針對空間構(gòu)型繩系編隊的地面試驗研究。為了對繩系編隊系統(tǒng)的可行性、可靠性等做更加深入的研究,復(fù)雜結(jié)構(gòu)繩系編隊的地面試驗將是未來的研究重點。

    6) 目前已開展的地面試驗裝置大部分還是以氣浮臺的方式進行研究,無法真正模擬太空失重環(huán)境。而失重飛機飛行法成本過高,維護周期太長。因此,未來可考慮建立更多微重力環(huán)境下的試驗,例如可利用水下環(huán)境或真空罐進行試驗裝置搭建,目前在繩系系統(tǒng)領(lǐng)域還沒有這方面的研究。

    猜你喜歡
    繩系系繩構(gòu)型
    航天員出艙可伸縮安全系繩設(shè)計與驗證
    基于模糊應(yīng)力-強度模型的空間系繩碰撞可靠性
    分子和離子立體構(gòu)型的判定
    電動力繩系離軌技術(shù)性能與任務(wù)適應(yīng)性分析
    空間繩系機器人目標抓捕魯棒自適應(yīng)控制器設(shè)計
    航天器受迫繞飛構(gòu)型設(shè)計與控制
    繩系衛(wèi)星被動釋放無超調(diào)脈寬脈頻調(diào)制控制
    空間繩系拖拽系統(tǒng)擺動特性與平穩(wěn)控制
    遙感衛(wèi)星平臺與載荷一體化構(gòu)型
    兩個具stp三維拓撲構(gòu)型的稀土配位聚合物{[Ln2(pda)3(H2O)2]·2H2O}n(Ln=Nd,La)
    国产精品秋霞免费鲁丝片| 国产欧美另类精品又又久久亚洲欧美| 91午夜精品亚洲一区二区三区| 亚洲在久久综合| 在线观看人妻少妇| 亚洲精品成人av观看孕妇| 黄色毛片三级朝国网站| 极品人妻少妇av视频| 五月开心婷婷网| 国产色爽女视频免费观看| 丰满饥渴人妻一区二区三| 女性被躁到高潮视频| 成人影院久久| 2018国产大陆天天弄谢| 免费播放大片免费观看视频在线观看| 亚洲av日韩在线播放| tube8黄色片| 久热久热在线精品观看| 亚洲综合精品二区| 国产一区亚洲一区在线观看| 超色免费av| 美女脱内裤让男人舔精品视频| 热99久久久久精品小说推荐| 欧美一级a爱片免费观看看| 国产一级毛片在线| 99久国产av精品国产电影| 亚洲精品av麻豆狂野| 色哟哟·www| 日本午夜av视频| 2022亚洲国产成人精品| 国产男女超爽视频在线观看| 午夜免费男女啪啪视频观看| 在线观看人妻少妇| 高清毛片免费看| 欧美丝袜亚洲另类| 亚洲欧美中文字幕日韩二区| 精品卡一卡二卡四卡免费| 最近中文字幕高清免费大全6| 各种免费的搞黄视频| 69精品国产乱码久久久| 制服诱惑二区| 日韩人妻高清精品专区| 免费高清在线观看日韩| 午夜福利,免费看| 精品人妻在线不人妻| av线在线观看网站| 欧美亚洲日本最大视频资源| 黄色视频在线播放观看不卡| 成年人午夜在线观看视频| 在线观看免费高清a一片| 久久久久人妻精品一区果冻| 18禁裸乳无遮挡动漫免费视频| 亚洲怡红院男人天堂| 国产精品免费大片| 在线观看一区二区三区激情| 久久青草综合色| 久久99热这里只频精品6学生| 男女国产视频网站| 91久久精品国产一区二区三区| 伦理电影免费视频| 久久精品国产自在天天线| 日韩人妻高清精品专区| 22中文网久久字幕| 夜夜看夜夜爽夜夜摸| 国产成人精品无人区| 高清午夜精品一区二区三区| 成人免费观看视频高清| 成年人免费黄色播放视频| 在线播放无遮挡| 97超碰精品成人国产| kizo精华| 亚洲无线观看免费| 国产精品不卡视频一区二区| 亚洲欧美成人精品一区二区| 80岁老熟妇乱子伦牲交| 久久久精品免费免费高清| 亚洲国产精品国产精品| 免费黄网站久久成人精品| 国产乱来视频区| 午夜激情av网站| 国产精品99久久99久久久不卡 | 岛国毛片在线播放| 午夜免费观看性视频| 中文字幕久久专区| 亚洲少妇的诱惑av| 免费黄频网站在线观看国产| 丁香六月天网| 亚洲av在线观看美女高潮| 久久精品国产a三级三级三级| 晚上一个人看的免费电影| 国产精品一二三区在线看| 999精品在线视频| 欧美3d第一页| 亚洲精品456在线播放app| 晚上一个人看的免费电影| 国产亚洲一区二区精品| 免费看av在线观看网站| 国产视频内射| 啦啦啦啦在线视频资源| 国产亚洲一区二区精品| 超碰97精品在线观看| 飞空精品影院首页| 亚洲精品第二区| 51国产日韩欧美| 丰满饥渴人妻一区二区三| 麻豆精品久久久久久蜜桃| 久久人人爽人人片av| 精品人妻熟女av久视频| 狂野欧美激情性bbbbbb| 国产一级毛片在线| 亚洲精品一区蜜桃| 在线观看人妻少妇| 欧美97在线视频| 考比视频在线观看| 午夜福利视频在线观看免费| 精品视频人人做人人爽| 少妇的逼水好多| 中国美白少妇内射xxxbb| 亚洲,欧美,日韩| 日韩,欧美,国产一区二区三区| 亚洲性久久影院| 欧美人与善性xxx| 久久久国产精品麻豆| 夫妻性生交免费视频一级片| av.在线天堂| 久久久午夜欧美精品| 免费观看a级毛片全部| 亚洲国产av新网站| 国产精品免费大片| 18禁在线播放成人免费| 日本av手机在线免费观看| 一区二区三区四区激情视频| 日本av手机在线免费观看| 18禁在线无遮挡免费观看视频| 午夜91福利影院| 日韩免费高清中文字幕av| 国产精品一区www在线观看| 欧美日韩av久久| 成年av动漫网址| 久久国产精品男人的天堂亚洲 | 一级黄片播放器| 啦啦啦啦在线视频资源| 久久久久久人妻| 丝袜喷水一区| 在线播放无遮挡| 少妇高潮的动态图| 国产高清不卡午夜福利| 成年美女黄网站色视频大全免费 | 青春草国产在线视频| 精品久久久噜噜| 高清不卡的av网站| 毛片一级片免费看久久久久| 久久国内精品自在自线图片| 国模一区二区三区四区视频| 久久久久精品性色| 久久精品久久久久久久性| 久久精品国产亚洲av天美| 美女中出高潮动态图| 国产精品99久久久久久久久| videosex国产| 满18在线观看网站| 男女边摸边吃奶| 亚洲熟女精品中文字幕| 亚洲一级一片aⅴ在线观看| 亚洲内射少妇av| 最近中文字幕高清免费大全6| 最近手机中文字幕大全| 久久精品国产亚洲网站| 99久久人妻综合| h视频一区二区三区| 欧美日本中文国产一区发布| 在线天堂最新版资源| 男女边摸边吃奶| av福利片在线| 中文精品一卡2卡3卡4更新| 亚洲国产精品一区三区| 亚洲成人手机| 三级国产精品片| 午夜福利视频在线观看免费| 在线精品无人区一区二区三| 国产av码专区亚洲av| 欧美日韩成人在线一区二区| av免费在线看不卡| 国产精品女同一区二区软件| 欧美 亚洲 国产 日韩一| 国产亚洲一区二区精品| 男女无遮挡免费网站观看| 国产一区二区三区综合在线观看 | 欧美日本中文国产一区发布| 纵有疾风起免费观看全集完整版| 日本wwww免费看| 国产亚洲午夜精品一区二区久久| 国产成人免费无遮挡视频| 久久人人爽人人片av| 精品视频人人做人人爽| 成人免费观看视频高清| 欧美精品一区二区大全| 少妇人妻 视频| 伦理电影免费视频| 大片免费播放器 马上看| 亚洲五月色婷婷综合| 丁香六月天网| 精品少妇黑人巨大在线播放| 欧美精品一区二区大全| 人人妻人人澡人人看| 在线观看国产h片| 黄色视频在线播放观看不卡| 亚洲少妇的诱惑av| 亚洲在久久综合| 天天操日日干夜夜撸| 久久久久视频综合| 久久青草综合色| 在线免费观看不下载黄p国产| 久久久久久久国产电影| 日韩精品免费视频一区二区三区 | 香蕉精品网在线| 亚洲丝袜综合中文字幕| 精品国产一区二区三区久久久樱花| 国产黄色免费在线视频| 精品久久蜜臀av无| 最后的刺客免费高清国语| 欧美3d第一页| 午夜免费男女啪啪视频观看| 日本黄大片高清| 国产一区二区在线观看av| 日韩精品有码人妻一区| 青春草视频在线免费观看| 18禁在线无遮挡免费观看视频| 国产女主播在线喷水免费视频网站| 成年av动漫网址| 国产精品熟女久久久久浪| 日本黄色片子视频| 午夜福利视频在线观看免费| 人妻一区二区av| 交换朋友夫妻互换小说| 搡女人真爽免费视频火全软件| 99九九线精品视频在线观看视频| 国产精品嫩草影院av在线观看| 国产在线一区二区三区精| 在线精品无人区一区二区三| 高清黄色对白视频在线免费看| 精品久久久久久电影网| 我要看黄色一级片免费的| 久久精品熟女亚洲av麻豆精品| 精品人妻偷拍中文字幕| 丰满迷人的少妇在线观看| 在线观看国产h片| 妹子高潮喷水视频| 丝袜喷水一区| 日韩,欧美,国产一区二区三区| 九九在线视频观看精品| 亚洲美女视频黄频| 精品一区在线观看国产| 美女国产高潮福利片在线看| 精品99又大又爽又粗少妇毛片| 亚洲欧美成人精品一区二区| 永久免费av网站大全| 国产精品秋霞免费鲁丝片| 免费观看性生交大片5| 久久久久国产精品人妻一区二区| 精品人妻一区二区三区麻豆| 亚洲在久久综合| 最新的欧美精品一区二区| 少妇的逼好多水| 中文字幕精品免费在线观看视频 | 一区在线观看完整版| 一级a做视频免费观看| a级毛片黄视频| 国产亚洲精品第一综合不卡 | 不卡视频在线观看欧美| 久久热精品热| 又粗又硬又长又爽又黄的视频| 精品久久久精品久久久| 久久精品国产鲁丝片午夜精品| 大香蕉久久网| 国产精品一国产av| 亚洲精品国产色婷婷电影| 2021少妇久久久久久久久久久| 国产爽快片一区二区三区| 国产一区二区三区综合在线观看 | 超色免费av| 国产亚洲最大av| 人妻少妇偷人精品九色| 国产亚洲av片在线观看秒播厂| 亚洲成人av在线免费| 亚洲国产av影院在线观看| 午夜免费鲁丝| 999精品在线视频| 成年女人在线观看亚洲视频| 亚洲欧美一区二区三区黑人 | 大又大粗又爽又黄少妇毛片口| 亚洲欧美一区二区三区黑人 | 亚洲精品一区蜜桃| 国产欧美亚洲国产| 亚洲欧洲国产日韩| 王馨瑶露胸无遮挡在线观看| 久久精品国产鲁丝片午夜精品| 大陆偷拍与自拍| 亚洲av.av天堂| 一本一本综合久久| 久热这里只有精品99| 永久免费av网站大全| 丰满乱子伦码专区| a级毛色黄片| 国产探花极品一区二区| 日本午夜av视频| 天堂8中文在线网| 在线观看免费高清a一片| 中文字幕亚洲精品专区| 人妻制服诱惑在线中文字幕| 搡女人真爽免费视频火全软件| 亚洲av男天堂| 卡戴珊不雅视频在线播放| 美女cb高潮喷水在线观看| 91精品国产九色| 精品少妇黑人巨大在线播放| 中文欧美无线码| 一级毛片电影观看| 少妇 在线观看| 80岁老熟妇乱子伦牲交| 国产免费一级a男人的天堂| 久久久久视频综合| 国产老妇伦熟女老妇高清| 久久婷婷青草| 久久久午夜欧美精品| 国产精品三级大全| 亚洲国产色片| 久久婷婷青草| 日韩一本色道免费dvd| 国产亚洲午夜精品一区二区久久| 在线 av 中文字幕| 免费少妇av软件| 久久久久久久久久久久大奶| 亚州av有码| 国产精品嫩草影院av在线观看| 国产色婷婷99| 高清黄色对白视频在线免费看| 国产亚洲欧美精品永久| 精品人妻熟女av久视频| 少妇被粗大的猛进出69影院 | 蜜桃久久精品国产亚洲av| 国产精品一区二区在线观看99| 视频区图区小说| 美女内射精品一级片tv| 97在线人人人人妻| 成年美女黄网站色视频大全免费 | 特大巨黑吊av在线直播| 在线精品无人区一区二区三| 精品久久久久久久久亚洲| 久久国产亚洲av麻豆专区| 老熟女久久久| 日本免费在线观看一区| 亚洲国产精品专区欧美| 人妻少妇偷人精品九色| 一二三四中文在线观看免费高清| 丰满少妇做爰视频| 老女人水多毛片| 三级国产精品欧美在线观看| 国产亚洲午夜精品一区二区久久| 黄色视频在线播放观看不卡| 国产又色又爽无遮挡免| 99热国产这里只有精品6| 91精品伊人久久大香线蕉| 新久久久久国产一级毛片| 老司机影院毛片| 九色成人免费人妻av| 99热6这里只有精品| 99国产精品免费福利视频| 一区在线观看完整版| 男男h啪啪无遮挡| 又黄又爽又刺激的免费视频.| 熟女电影av网| 亚洲内射少妇av| 在线免费观看不下载黄p国产| 黄色一级大片看看| 国产男女内射视频| 国产欧美日韩一区二区三区在线 | 男女边吃奶边做爰视频| 午夜激情av网站| 少妇被粗大的猛进出69影院 | tube8黄色片| 欧美精品一区二区大全| 狂野欧美激情性xxxx在线观看| 午夜激情福利司机影院| 亚洲精品美女久久av网站| 久久久久视频综合| 久久久久网色| 国产又色又爽无遮挡免| 精品久久久久久久久亚洲| 能在线免费看毛片的网站| 丝袜喷水一区| 午夜日本视频在线| 日韩成人av中文字幕在线观看| 亚洲不卡免费看| 国产免费福利视频在线观看| 9色porny在线观看| 国产欧美日韩一区二区三区在线 | 国产精品久久久久久av不卡| 亚洲国产毛片av蜜桃av| 边亲边吃奶的免费视频| 欧美少妇被猛烈插入视频| 嘟嘟电影网在线观看| 成人亚洲精品一区在线观看| 亚洲美女视频黄频| 久久久久精品性色| videosex国产| 99九九线精品视频在线观看视频| 极品少妇高潮喷水抽搐| 国产熟女欧美一区二区| 精品久久蜜臀av无| 国产日韩欧美视频二区| 国产精品.久久久| 亚洲精品色激情综合| 久久精品久久久久久久性| 国产精品一区二区三区四区免费观看| 97超视频在线观看视频| 日韩一区二区三区影片| 一个人看视频在线观看www免费| 国产成人91sexporn| 老司机亚洲免费影院| 欧美人与善性xxx| 一区二区av电影网| 高清欧美精品videossex| 亚洲性久久影院| 九九久久精品国产亚洲av麻豆| 好男人视频免费观看在线| 久久这里有精品视频免费| 3wmmmm亚洲av在线观看| 亚洲欧美清纯卡通| 丝袜喷水一区| 日本av免费视频播放| 一区二区三区免费毛片| 蜜桃在线观看..| 日本色播在线视频| av女优亚洲男人天堂| 不卡视频在线观看欧美| 777米奇影视久久| 久久久久久久大尺度免费视频| 成人毛片60女人毛片免费| 啦啦啦中文免费视频观看日本| 欧美一级a爱片免费观看看| 18在线观看网站| 我要看黄色一级片免费的| 日韩在线高清观看一区二区三区| 亚洲综合色惰| 亚洲欧美清纯卡通| 91成人精品电影| 街头女战士在线观看网站| videossex国产| 伦精品一区二区三区| 久久久国产一区二区| 欧美成人精品欧美一级黄| 成年美女黄网站色视频大全免费 | 国产又色又爽无遮挡免| 中国三级夫妇交换| 亚洲人成网站在线播| 欧美精品一区二区大全| 热re99久久国产66热| 国产精品国产av在线观看| 亚洲美女黄色视频免费看| 亚洲国产色片| 久久国产精品大桥未久av| 久久午夜综合久久蜜桃| 国产精品麻豆人妻色哟哟久久| 亚洲精品,欧美精品| 曰老女人黄片| 国产男女超爽视频在线观看| 极品人妻少妇av视频| 韩国高清视频一区二区三区| 欧美最新免费一区二区三区| 青春草亚洲视频在线观看| 国产日韩欧美在线精品| 最新中文字幕久久久久| 亚洲人与动物交配视频| 欧美xxⅹ黑人| 王馨瑶露胸无遮挡在线观看| 亚洲无线观看免费| 777米奇影视久久| 久久精品国产a三级三级三级| 纵有疾风起免费观看全集完整版| 国产在视频线精品| 久久鲁丝午夜福利片| 伦理电影免费视频| 一个人免费看片子| 久久精品国产自在天天线| 国产精品一二三区在线看| 男女啪啪激烈高潮av片| 夜夜骑夜夜射夜夜干| 亚洲经典国产精华液单| 国产片内射在线| 少妇精品久久久久久久| 人人妻人人澡人人看| 麻豆乱淫一区二区| 中文天堂在线官网| 成人无遮挡网站| 一个人看视频在线观看www免费| 国产国拍精品亚洲av在线观看| 黑人巨大精品欧美一区二区蜜桃 | 夜夜看夜夜爽夜夜摸| 日日爽夜夜爽网站| 国产午夜精品一二区理论片| 又粗又硬又长又爽又黄的视频| 极品少妇高潮喷水抽搐| 国产国语露脸激情在线看| 欧美3d第一页| freevideosex欧美| 日韩亚洲欧美综合| 夫妻性生交免费视频一级片| 黄色怎么调成土黄色| 国产av一区二区精品久久| 亚洲一级一片aⅴ在线观看| 中文字幕制服av| 国产精品无大码| 51国产日韩欧美| 欧美精品一区二区大全| 2021少妇久久久久久久久久久| 国产男女超爽视频在线观看| 午夜免费鲁丝| h视频一区二区三区| 大码成人一级视频| 日日啪夜夜爽| 欧美日韩精品成人综合77777| 亚洲人与动物交配视频| 日本-黄色视频高清免费观看| 美女内射精品一级片tv| 大香蕉久久网| 亚洲色图综合在线观看| 蜜桃在线观看..| 日韩强制内射视频| 精品视频人人做人人爽| 日韩中文字幕视频在线看片| 丁香六月天网| 久久久久久久久大av| 丝瓜视频免费看黄片| 国产极品粉嫩免费观看在线 | 久久久欧美国产精品| 亚洲内射少妇av| 成人综合一区亚洲| 免费不卡的大黄色大毛片视频在线观看| 自线自在国产av| 亚洲欧洲国产日韩| 狠狠精品人妻久久久久久综合| 欧美精品一区二区免费开放| 天堂俺去俺来也www色官网| 亚洲精品久久久久久婷婷小说| 最近2019中文字幕mv第一页| 日本-黄色视频高清免费观看| 亚洲精品国产av成人精品| 99国产综合亚洲精品| 国产成人a∨麻豆精品| 看免费成人av毛片| 国产高清有码在线观看视频| 伊人久久精品亚洲午夜| 久久人人爽av亚洲精品天堂| 熟女av电影| 中文字幕最新亚洲高清| 少妇的逼水好多| 午夜久久久在线观看| 一级毛片aaaaaa免费看小| 午夜激情久久久久久久| 美女中出高潮动态图| 亚州av有码| 国产高清三级在线| 国产成人免费观看mmmm| 丁香六月天网| 97超视频在线观看视频| 91成人精品电影| 美女国产视频在线观看| av黄色大香蕉| 女人精品久久久久毛片| 一个人免费看片子| 免费不卡的大黄色大毛片视频在线观看| 丁香六月天网| 一本—道久久a久久精品蜜桃钙片| 国产精品麻豆人妻色哟哟久久| 国产精品国产三级国产av玫瑰| 国产亚洲精品第一综合不卡 | 最黄视频免费看| 嘟嘟电影网在线观看| 哪个播放器可以免费观看大片| 纯流量卡能插随身wifi吗| 成人国产麻豆网| 黑人巨大精品欧美一区二区蜜桃 | 亚洲综合色惰| 极品少妇高潮喷水抽搐| 中国三级夫妇交换| 精品一区二区三卡| av又黄又爽大尺度在线免费看| 男女无遮挡免费网站观看| 日韩av免费高清视频| 少妇的逼水好多| 亚洲国产精品一区三区| 看免费成人av毛片| 日韩欧美精品免费久久| 熟女电影av网| a级毛片黄视频| 一级,二级,三级黄色视频| a级毛色黄片| av福利片在线| 欧美成人午夜免费资源| 女人精品久久久久毛片| 国产色婷婷99| 亚洲欧洲精品一区二区精品久久久 | 日韩一本色道免费dvd| 一级黄片播放器| av免费观看日本| 三上悠亚av全集在线观看| 亚洲综合色惰| 国产免费又黄又爽又色| 国产亚洲午夜精品一区二区久久| a级片在线免费高清观看视频| 日韩成人av中文字幕在线观看| 看十八女毛片水多多多| av电影中文网址| 我的老师免费观看完整版|