• <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)
    视频中文字幕在线观看| 亚洲欧美清纯卡通| 一二三四在线观看免费中文在 | 中文字幕制服av| 999精品在线视频| 午夜精品国产一区二区电影| 久热这里只有精品99| 91精品三级在线观看| 国产免费现黄频在线看| 一级,二级,三级黄色视频| 日韩中文字幕视频在线看片| 99久久综合免费| 中文字幕精品免费在线观看视频 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av在线观看美女高潮| 日本91视频免费播放| videos熟女内射| 蜜桃在线观看..| 成人毛片a级毛片在线播放| 国产精品一国产av| 国产亚洲午夜精品一区二区久久| 国产精品久久久久成人av| 少妇熟女欧美另类| 国产免费又黄又爽又色| 天天操日日干夜夜撸| 99热国产这里只有精品6| 女人精品久久久久毛片| 2021少妇久久久久久久久久久| 亚洲欧美中文字幕日韩二区| 国产亚洲av片在线观看秒播厂| a级片在线免费高清观看视频| 国产成人精品福利久久| 在线看a的网站| 26uuu在线亚洲综合色| 观看美女的网站| 国产精品嫩草影院av在线观看| 99热这里只有是精品在线观看| 熟女人妻精品中文字幕| 色婷婷久久久亚洲欧美| 国产日韩欧美视频二区| 亚洲人与动物交配视频| 日韩欧美精品免费久久| 亚洲情色 制服丝袜| 亚洲av福利一区| av一本久久久久| 日本猛色少妇xxxxx猛交久久| 国产精品三级大全| 又粗又硬又长又爽又黄的视频| 最近最新中文字幕免费大全7| 欧美精品一区二区免费开放| 成人黄色视频免费在线看| 午夜日本视频在线| 永久免费av网站大全| 亚洲第一区二区三区不卡| 热99久久久久精品小说推荐| av线在线观看网站| 51国产日韩欧美| 日韩 亚洲 欧美在线| 国产一级毛片在线| 久久鲁丝午夜福利片| 高清黄色对白视频在线免费看| 观看av在线不卡| 日本爱情动作片www.在线观看| av不卡在线播放| 欧美成人精品欧美一级黄| 国产一区二区三区综合在线观看 | 又黄又粗又硬又大视频| 亚洲成人一二三区av| 日本欧美视频一区| 久久久久国产网址| av免费在线看不卡| 精品少妇内射三级| 亚洲人与动物交配视频| 午夜福利乱码中文字幕| 人体艺术视频欧美日本| 色吧在线观看| 少妇猛男粗大的猛烈进出视频| 永久网站在线| 亚洲综合色惰| 2021少妇久久久久久久久久久| 亚洲欧美色中文字幕在线| 国产一区二区三区综合在线观看 | 国产精品久久久久成人av| tube8黄色片| 欧美人与善性xxx| 成人毛片60女人毛片免费| 色吧在线观看| 夜夜骑夜夜射夜夜干| 免费看av在线观看网站| 久久久国产精品麻豆| 国产精品偷伦视频观看了| 亚洲内射少妇av| 国产 精品1| 久久久a久久爽久久v久久| 免费看光身美女| 亚洲色图 男人天堂 中文字幕 | 欧美xxxx性猛交bbbb| 精品视频人人做人人爽| 久久久久久久久久人人人人人人| 国产极品粉嫩免费观看在线| 午夜老司机福利剧场| 国产无遮挡羞羞视频在线观看| 免费人成在线观看视频色| 精品人妻偷拍中文字幕| 一级毛片 在线播放| 久久国产亚洲av麻豆专区| 免费黄色在线免费观看| 女人久久www免费人成看片| av播播在线观看一区| 国产激情久久老熟女| 精品一区二区三区视频在线| 免费黄频网站在线观看国产| 国产精品一区二区在线观看99| 亚洲av.av天堂| 一区二区三区四区激情视频| 国产成人欧美| 日韩,欧美,国产一区二区三区| 亚洲精品久久久久久婷婷小说| 国产精品女同一区二区软件| 亚洲欧美一区二区三区国产| av黄色大香蕉| 2022亚洲国产成人精品| 成人影院久久| 午夜精品国产一区二区电影| 人人妻人人澡人人爽人人夜夜| 色吧在线观看| 性色avwww在线观看| 国产乱来视频区| 99国产精品免费福利视频| 国产精品 国内视频| 99热网站在线观看| 少妇 在线观看| 又黄又爽又刺激的免费视频.| 国产女主播在线喷水免费视频网站| av片东京热男人的天堂| 人妻人人澡人人爽人人| 国产av一区二区精品久久| 99精国产麻豆久久婷婷| 亚洲av.av天堂| 在线 av 中文字幕| 午夜免费观看性视频| 不卡视频在线观看欧美| 大片电影免费在线观看免费| 色婷婷久久久亚洲欧美| 亚洲图色成人| 香蕉丝袜av| 一级a做视频免费观看| 国产免费视频播放在线视频| 成人二区视频| 国产有黄有色有爽视频| 国产精品人妻久久久影院| 又黄又粗又硬又大视频| 亚洲,一卡二卡三卡| 国产成人精品在线电影| 国产成人一区二区在线| 香蕉精品网在线| 麻豆乱淫一区二区| 精品亚洲成a人片在线观看| 亚洲欧美精品自产自拍| 欧美另类一区| 蜜桃国产av成人99| 国产女主播在线喷水免费视频网站| 爱豆传媒免费全集在线观看| av女优亚洲男人天堂| 午夜福利,免费看| freevideosex欧美| 国产精品女同一区二区软件| 欧美97在线视频| 中国美白少妇内射xxxbb| 亚洲精华国产精华液的使用体验| 欧美日韩视频精品一区| 国产国拍精品亚洲av在线观看| 国产av国产精品国产| 日本av免费视频播放| 亚洲国产欧美日韩在线播放| a 毛片基地| 久久久久精品性色| 久久精品aⅴ一区二区三区四区 | 欧美丝袜亚洲另类| 肉色欧美久久久久久久蜜桃| 日韩av免费高清视频| 国产精品久久久久久精品电影小说| 国产精品不卡视频一区二区| 亚洲av国产av综合av卡| 国产精品一区二区在线不卡| 国产一区有黄有色的免费视频| 久久久久久久久久久免费av| 桃花免费在线播放| 国产黄色视频一区二区在线观看| 尾随美女入室| 久久久久视频综合| 九色亚洲精品在线播放| 99精国产麻豆久久婷婷| 日本猛色少妇xxxxx猛交久久| 99久久精品国产国产毛片| 国产成人午夜福利电影在线观看| 观看av在线不卡| 成人影院久久| 亚洲少妇的诱惑av| 天堂中文最新版在线下载| 国产高清国产精品国产三级| 超碰97精品在线观看| 国产毛片在线视频| 国产熟女欧美一区二区| 永久免费av网站大全| 精品久久国产蜜桃| 精品99又大又爽又粗少妇毛片| 亚洲经典国产精华液单| 日韩伦理黄色片| 亚洲欧美日韩卡通动漫| 高清在线视频一区二区三区| 91午夜精品亚洲一区二区三区| 国产片内射在线| 日韩 亚洲 欧美在线| 久久久欧美国产精品| 国产极品天堂在线| av在线播放精品| 亚洲国产看品久久| 国产精品成人在线| 国产老妇伦熟女老妇高清| 午夜视频国产福利| 18在线观看网站| 亚洲成人手机| 久久国产精品男人的天堂亚洲 | 日本爱情动作片www.在线观看| 狂野欧美激情性bbbbbb| 九九在线视频观看精品| 男女边吃奶边做爰视频| 午夜福利视频在线观看免费| 男女边摸边吃奶| 少妇猛男粗大的猛烈进出视频| 久久精品国产a三级三级三级| 久久久久久久久久久久大奶| 国产成人免费观看mmmm| 久久韩国三级中文字幕| videos熟女内射| 亚洲成av片中文字幕在线观看 | 亚洲av欧美aⅴ国产| 国产亚洲午夜精品一区二区久久| 人妻系列 视频| 亚洲成人av在线免费| av免费在线看不卡| 9热在线视频观看99| 精品人妻熟女毛片av久久网站| 搡女人真爽免费视频火全软件| 狂野欧美激情性bbbbbb| 精品国产一区二区三区久久久樱花| 成人国语在线视频| 妹子高潮喷水视频| 精品国产一区二区三区久久久樱花| 色婷婷av一区二区三区视频| 免费日韩欧美在线观看| 中国美白少妇内射xxxbb| 777米奇影视久久| 亚洲欧美一区二区三区国产| 国产在线视频一区二区| 国产精品一区二区在线不卡| 精品亚洲成国产av| 看非洲黑人一级黄片| 99久久综合免费| 国产高清三级在线| 亚洲色图 男人天堂 中文字幕 | 久久99一区二区三区| 久久人人爽av亚洲精品天堂| 成人午夜精彩视频在线观看| 欧美日韩精品成人综合77777| 亚洲国产色片| 亚洲第一区二区三区不卡| 啦啦啦中文免费视频观看日本| av在线app专区| 91国产中文字幕| 久久 成人 亚洲| 99热全是精品| 亚洲精品视频女| 最近中文字幕高清免费大全6| 国产xxxxx性猛交| freevideosex欧美| 丰满少妇做爰视频| 黄色一级大片看看| 伦理电影免费视频| 国产成人精品一,二区| 久久久久久久久久久久大奶| 在现免费观看毛片| 精品酒店卫生间| 国产欧美日韩综合在线一区二区| www日本在线高清视频| 少妇的逼水好多| 两性夫妻黄色片 | 乱人伦中国视频| 狠狠精品人妻久久久久久综合| 国产精品麻豆人妻色哟哟久久| 永久网站在线| 人人澡人人妻人| 国产在线一区二区三区精| 街头女战士在线观看网站| av国产精品久久久久影院| 亚洲精品一区蜜桃| 69精品国产乱码久久久| 中文乱码字字幕精品一区二区三区| 男人添女人高潮全过程视频| 国产黄色免费在线视频| 精品卡一卡二卡四卡免费| 久久亚洲国产成人精品v| 国产极品粉嫩免费观看在线| 国产老妇伦熟女老妇高清| 亚洲中文av在线| 观看av在线不卡| 2018国产大陆天天弄谢| 91在线精品国自产拍蜜月| 国产成人精品无人区| 蜜臀久久99精品久久宅男| 人人妻人人澡人人看| 欧美人与性动交α欧美软件 | 久久国产精品大桥未久av| 久久人人爽av亚洲精品天堂| videos熟女内射| 久久人妻熟女aⅴ| 国产成人免费无遮挡视频| 国产成人精品一,二区| 日韩精品有码人妻一区| 狂野欧美激情性bbbbbb| 欧美日韩一区二区视频在线观看视频在线| 久久影院123| 午夜免费观看性视频| 欧美人与性动交α欧美精品济南到 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 99视频精品全部免费 在线| av片东京热男人的天堂| 亚洲欧美中文字幕日韩二区| 日本黄色日本黄色录像| 大片电影免费在线观看免费| 涩涩av久久男人的天堂| 亚洲天堂av无毛| 黄色毛片三级朝国网站| 亚洲精品一二三| 国产国语露脸激情在线看| 少妇熟女欧美另类| 熟女电影av网| 亚洲精品国产av成人精品| 五月开心婷婷网| 久久久国产一区二区| 丁香六月天网| 蜜臀久久99精品久久宅男| 亚洲婷婷狠狠爱综合网| 久久久久久久久久久免费av| 亚洲精品国产色婷婷电影| 亚洲国产色片| 伦理电影大哥的女人| 一边亲一边摸免费视频| 少妇精品久久久久久久| 69精品国产乱码久久久| 国产精品国产三级国产专区5o| 国产精品一区二区在线观看99| 菩萨蛮人人尽说江南好唐韦庄| 大香蕉97超碰在线| 亚洲精品456在线播放app| 国产精品人妻久久久久久| 交换朋友夫妻互换小说| 国产精品女同一区二区软件| 久久精品国产亚洲av天美| 成人免费观看视频高清| av国产精品久久久久影院| 日本欧美视频一区| 精品午夜福利在线看| 男人舔女人的私密视频| 视频区图区小说| 少妇人妻精品综合一区二区| 国产成人精品福利久久| 啦啦啦视频在线资源免费观看| 国产精品久久久av美女十八| 久久久久视频综合| 99热全是精品| 欧美精品一区二区免费开放| 曰老女人黄片| 亚洲精品乱码久久久久久按摩| 啦啦啦中文免费视频观看日本| 亚洲精品美女久久久久99蜜臀 | www日本在线高清视频| 91精品三级在线观看| xxx大片免费视频| 午夜福利视频在线观看免费| 日本午夜av视频| 亚洲精品456在线播放app| 大话2 男鬼变身卡| 卡戴珊不雅视频在线播放| 女性生殖器流出的白浆| 啦啦啦中文免费视频观看日本| 伦精品一区二区三区| av又黄又爽大尺度在线免费看| a 毛片基地| 纵有疾风起免费观看全集完整版| 国产免费一区二区三区四区乱码| 色网站视频免费| 久久人人爽av亚洲精品天堂| 国产亚洲精品第一综合不卡 | 另类精品久久| 成人黄色视频免费在线看| 人妻少妇偷人精品九色| 大话2 男鬼变身卡| 一区二区三区乱码不卡18| 久久久久国产网址| 亚洲情色 制服丝袜| 国产无遮挡羞羞视频在线观看| 国产精品无大码| 国产精品嫩草影院av在线观看| 免费观看在线日韩| 国产成人免费无遮挡视频| h视频一区二区三区| 午夜激情久久久久久久| 精品人妻一区二区三区麻豆| 亚洲av成人精品一二三区| 久久精品夜色国产| 亚洲经典国产精华液单| 春色校园在线视频观看| 亚洲中文av在线| 国产成人a∨麻豆精品| 亚洲精品一二三| 免费av中文字幕在线| 国产一区有黄有色的免费视频| 黄片无遮挡物在线观看| 久久人人爽人人片av| 久久99一区二区三区| 少妇猛男粗大的猛烈进出视频| 性高湖久久久久久久久免费观看| 亚洲精品久久午夜乱码| 国产精品.久久久| 欧美精品av麻豆av| 蜜桃在线观看..| 亚洲av福利一区| 国产黄色视频一区二区在线观看| 国产精品国产三级专区第一集| 欧美日韩国产mv在线观看视频| 伦理电影大哥的女人| 国产日韩一区二区三区精品不卡| 精品视频人人做人人爽| 亚洲精品国产色婷婷电影| 久久国产精品男人的天堂亚洲 | 纯流量卡能插随身wifi吗| 26uuu在线亚洲综合色| 国产精品99久久99久久久不卡 | 久久精品国产a三级三级三级| 男女啪啪激烈高潮av片| 99热这里只有是精品在线观看| 精品人妻一区二区三区麻豆| 美女主播在线视频| 涩涩av久久男人的天堂| 极品人妻少妇av视频| 亚洲精品成人av观看孕妇| 大香蕉久久成人网| 狂野欧美激情性bbbbbb| 成人手机av| 欧美精品亚洲一区二区| 视频在线观看一区二区三区| 欧美日韩av久久| 久久久亚洲精品成人影院| www.熟女人妻精品国产 | 精品亚洲乱码少妇综合久久| 天天操日日干夜夜撸| 青春草亚洲视频在线观看| 婷婷色av中文字幕| 国产在线免费精品| 久久午夜综合久久蜜桃| 丝瓜视频免费看黄片| 国产福利在线免费观看视频| 天堂8中文在线网| 精品亚洲乱码少妇综合久久| 人成视频在线观看免费观看| 精品第一国产精品| 精品一区二区三卡| 一区在线观看完整版| 狂野欧美激情性bbbbbb| 黄片无遮挡物在线观看| 黄网站色视频无遮挡免费观看| 波多野结衣一区麻豆| 亚洲精品456在线播放app| 久久久亚洲精品成人影院| a级毛片黄视频| 寂寞人妻少妇视频99o| 午夜老司机福利剧场| 亚洲中文av在线| 欧美人与性动交α欧美软件 | 国产有黄有色有爽视频| 中文乱码字字幕精品一区二区三区| 欧美精品人与动牲交sv欧美| 国产乱来视频区| 日产精品乱码卡一卡2卡三| 婷婷色av中文字幕| 亚洲精品,欧美精品| 最后的刺客免费高清国语| 国产一区有黄有色的免费视频| 免费观看a级毛片全部| 最近手机中文字幕大全| 中文字幕av电影在线播放| 午夜91福利影院| 亚洲三级黄色毛片| 最黄视频免费看| 久久久久网色| 男女下面插进去视频免费观看 | 男女午夜视频在线观看 | 婷婷成人精品国产| 精品福利永久在线观看| 精品久久久精品久久久| 深夜精品福利| 桃花免费在线播放| 欧美国产精品va在线观看不卡| 国产xxxxx性猛交| 国产极品天堂在线| 国产黄色视频一区二区在线观看| 国产免费福利视频在线观看| 亚洲精品成人av观看孕妇| 国产精品 国内视频| 亚洲成人一二三区av| freevideosex欧美| 51国产日韩欧美| 久久国产精品大桥未久av| 久久久精品94久久精品| 高清毛片免费看| 亚洲一码二码三码区别大吗| 亚洲欧美清纯卡通| 久久久久精品人妻al黑| 最后的刺客免费高清国语| 香蕉丝袜av| 少妇人妻精品综合一区二区| 久久久亚洲精品成人影院| 中文字幕亚洲精品专区| 亚洲av成人精品一二三区| 日韩精品免费视频一区二区三区 | 婷婷色麻豆天堂久久| 七月丁香在线播放| 中文字幕亚洲精品专区| 国产永久视频网站| 少妇熟女欧美另类| 亚洲精品,欧美精品| 99久久精品国产国产毛片| 久久毛片免费看一区二区三区| 伦理电影免费视频| av黄色大香蕉| 1024视频免费在线观看| 999精品在线视频| 国产精品成人在线| 国产 精品1| 国产一区二区激情短视频 | 丝袜美足系列| 老女人水多毛片| 最新的欧美精品一区二区| 久久青草综合色| 久久国产精品男人的天堂亚洲 | 欧美日韩国产mv在线观看视频| 久久久国产欧美日韩av| 日韩av在线免费看完整版不卡| 久久久国产一区二区| 久久久久久久大尺度免费视频| 天堂中文最新版在线下载| 男女国产视频网站| 十八禁高潮呻吟视频| 欧美日韩视频精品一区| 国产欧美日韩一区二区三区在线| 久久久久视频综合| 国产欧美日韩一区二区三区在线| 国产男女超爽视频在线观看| 久久午夜福利片| h视频一区二区三区| 亚洲美女视频黄频| 乱码一卡2卡4卡精品| av免费在线看不卡| 18禁国产床啪视频网站| 色视频在线一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 99久国产av精品国产电影| 国产精品一国产av| 自线自在国产av| 丰满乱子伦码专区| 少妇被粗大的猛进出69影院 | 久久精品久久久久久久性| 亚洲激情五月婷婷啪啪| 中文字幕免费在线视频6| 国产一区二区在线观看av| 80岁老熟妇乱子伦牲交| 在线观看国产h片| 国产成人午夜福利电影在线观看| 91久久精品国产一区二区三区| 精品亚洲成a人片在线观看| 亚洲情色 制服丝袜| 在线观看人妻少妇| 高清av免费在线| 18禁观看日本| 男女下面插进去视频免费观看 | 国产一区有黄有色的免费视频| 国产极品天堂在线| 国语对白做爰xxxⅹ性视频网站| 这个男人来自地球电影免费观看 | 久久免费观看电影| 精品亚洲乱码少妇综合久久| 国产男人的电影天堂91| 久久鲁丝午夜福利片| 大片免费播放器 马上看| 日韩av在线免费看完整版不卡| 亚洲精品中文字幕在线视频| 日韩不卡一区二区三区视频在线| 欧美3d第一页| 1024视频免费在线观看| 国产精品熟女久久久久浪| 亚洲av电影在线进入| 亚洲成人手机| 亚洲久久久国产精品| 2018国产大陆天天弄谢| 久久亚洲国产成人精品v| 色哟哟·www| 久久久久久久亚洲中文字幕| 欧美性感艳星| 精品酒店卫生间| 亚洲人与动物交配视频| 成人黄色视频免费在线看| 丁香六月天网| 一区二区av电影网|