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

    剛性航天器的預(yù)定義時間滑??刂?/h1>
    2022-01-24 04:51:18賽華陽徐振邦張恩陽
    光學精密工程 2021年12期
    關(guān)鍵詞:角速度航天器滑模

    賽華陽,徐振邦,賀 帥,張恩陽,秦 超

    (1.中科院長春光學精密機械與物理研究所 中科院空間光學系統(tǒng)在軌制造與集成重點實驗室,吉林 長春 130033;2.中國科學院大學,北京 100049;3.中國科學院大學材料與光電研究中心,北京 100049)

    1 引言

    近年來,航天器的姿態(tài)跟蹤問題由于在空間應(yīng)用中的重要作用而得到了廣泛的研究[1-2]。航天器的姿態(tài)調(diào)整是一個復雜的控制問題,因為其存在相互作用的非線性運動學和動力學模型,以及空間環(huán)境中不可預(yù)測的外界擾動[3-4]。

    滑??刂疲⊿liding Mode Control,SMC)是處理具有不確定性和外界擾動的非線性系統(tǒng)的最有效的方法之一[5]。由于其對系統(tǒng)不確定性和外界干擾的強魯棒性,SMC 已被廣泛應(yīng)用到航天器的姿態(tài)控制中。但傳統(tǒng)的SMC 算法只能獲得漸近穩(wěn)定的結(jié)果,這意味著航天器的姿態(tài)跟蹤誤差需要很長的時間才能收斂到平衡點,而許多任務(wù)要求航天器實現(xiàn)快速的姿態(tài)調(diào)整[6]。為了滿足這一要求,有限時間控制的概念被提出,它可以提高系統(tǒng)誤差的收斂速度,使其在有限的時間內(nèi)收斂。終端滑??刂凭褪且环N典型的有限時間控制方法,其中,一種新型的非奇異終端滑模控制(Non-singular Terminal Sliding Mode Con?trol,NTSMC)已經(jīng)被應(yīng)用到航天器的姿態(tài)調(diào)整中[7]。但有限時間控制方案的收斂時間依賴于系統(tǒng)的初始狀態(tài),因此難以提前獲得航天器準確的姿態(tài)穩(wěn)定時間界限。與有限時間控制方案不同的是,固定時間控制可以保證系統(tǒng)的收斂時間獨立于系統(tǒng)的初始值,即系統(tǒng)收斂時間的上界是一個不依賴于系統(tǒng)初始狀態(tài)的常數(shù)[8]。目前,一些固定時間的SMC 也已經(jīng)被成功地應(yīng)用于航天器的姿態(tài)跟蹤控制中[1,9]。

    雖然與有限時間控制相比,固定時間在系統(tǒng)穩(wěn)定時間的建立上有顯著的優(yōu)勢,但固定時間控制器的控制參數(shù)與穩(wěn)定時間之間的關(guān)系往往是復雜的。為了進一步解決這個問題,一種更先進的控制概念,即預(yù)定義時間控制被提出[10]。預(yù)定義時間控制為被控系統(tǒng)提供了一種先進的穩(wěn)定性特征,即通過調(diào)整預(yù)定義時間參數(shù)可以很容易地確定系統(tǒng)的穩(wěn)定性上界,從而為系統(tǒng)的行為提供高度的確定性。在文獻[10-12]中,預(yù)定義時間控制被用來和SMC 相結(jié)合以提高一階控制系統(tǒng)的魯棒性,但這些方法僅能保證系統(tǒng)在SMC中到達階段的預(yù)定義時間穩(wěn)定性。在文獻[13-15]中,預(yù)定義時間滑模控制的相關(guān)理論被進一步提出,并被拓展到二階控制系統(tǒng)。此外,在文獻[3]中,預(yù)定義時間控制被應(yīng)用于航天器的姿態(tài)跟蹤控制中。但上述預(yù)定義時間控制方案的穩(wěn)定時間上界往往過于保守,這意味著控制器會遠提前于所設(shè)定的穩(wěn)定時間完成航天器的姿態(tài)調(diào)整。對于航天器姿態(tài)控制而言,過快的姿態(tài)調(diào)節(jié)會導致系統(tǒng)能量的浪費,甚至使控制力矩超過系統(tǒng)所允許的范圍。

    在本文中,我們提出了一種應(yīng)用于剛性航天器姿態(tài)跟蹤控制的新型預(yù)定義時間滑模控制(Predefined-Time Sliding Mode Control,PTSMC)方案用于解決現(xiàn)有預(yù)定義時間控制方案的穩(wěn)定時間上界過于保守的問題。首先,我們設(shè)計了一種新型滑模變量和預(yù)定義時間滑模面,然后提出了一種新型PTSMC 方案,并采用邊界層技術(shù)降低了系統(tǒng)抖振。接著,通過Lyapunov函數(shù)證明了所提出的控制器的預(yù)定義時間穩(wěn)定性和系統(tǒng)收斂上界的非保守性。最后通過仿真實驗,分別與現(xiàn)有的預(yù)定義時間控制方案,NTSMC 方案和比例微分(Proportional-Deriva?tive,PD)控制進行了對比,并利用3 自由度氣浮平臺進行了實驗驗證,證明了所提出控制方案的有效性和優(yōu)越性。

    2 剛性航天器模型和性質(zhì)

    如圖1 所示,在建立航天器的姿態(tài)控制系統(tǒng)時常建立3 個參考坐標系,包括慣性坐標系FI,參考坐標系FO和航天器本體坐標系FB。在本文中,F(xiàn)I被選擇為以地球為中心的慣性框架,F(xiàn)B的原點位于被控航天器的質(zhì)心,并圍繞FI以ω0∈R+的軌道速率運行。FO與FB有相同的原點,其滾轉(zhuǎn)軸XO指向飛行器的飛行方向,俯仰軸YO垂直于飛行器的運行軌道,偏航軸ZO指向地心。

    圖1 航天器的參考坐標系Fig.1 Coordinate reference frame system of spacecraft

    剛性航天器的姿態(tài)可以通過航天器本體坐標系FB相對于慣性坐標系FI的姿態(tài)四元數(shù)進行表示,即(qv,q4)∈R3×R。剛性航天器的運動學方程可以表示為:

    其中:ω∈R3表示航天器相對于本體坐標系FB的角速度,I3表示單位矩陣,(·)×表示為由向量元素構(gòu)建的偏斜對稱矩陣,即:

    航天器的角速度誤差可以通過下式給出:

    其中,C表示從FO到FB的坐標變換矩陣,并且矩陣C可以通過下式進行計算:

    此外,航天器的姿態(tài)動力學方程可以表示為:

    其中:J∈R3×3是航天器關(guān)于本體坐標系FB的慣性矩陣,u∈R3和d∈R3分別表示系統(tǒng)的控制力矩和未知但有界的擾動力矩。

    結(jié)合式(3),(5)和(7),航天器姿態(tài)的誤差運動學和動力學方程可以被寫為:

    考慮到航天器模型的不確定性,可以假設(shè):

    其中,J0和ΔJ分別表示已知的系統(tǒng)的名義慣性矩陣和慣性矩陣J的不確定部分。式(8)中給出的動力學方程可以被寫為:

    其中,系統(tǒng)耦合不確定部分ρ∈R3被定義為如下形式:

    參考文獻[1],系統(tǒng)耦合不確定部分的上界可以通過下式確定:

    本文要解決的問題如下:通過設(shè)計控制輸入u,使包含系統(tǒng)不確定性和外界干擾的剛性航天器的姿態(tài)誤差ev和角速度跟蹤誤差ωe在給定的時間內(nèi)收斂到零。

    3 預(yù)定義時間穩(wěn)定性的基本理論

    首先介紹關(guān)于有限時間、固定時間和預(yù)定時間穩(wěn)定性的一些的定義和引理以用于控制器預(yù)定義時間穩(wěn)定性的證明。

    考慮一個非線性系統(tǒng):

    其中:x∈Rn表示系統(tǒng)狀態(tài),向量ρ∈Rb代表系統(tǒng)(13)的系統(tǒng)參數(shù)。函數(shù)f:Rn→Rn可以被認為是非線性和連續(xù)的,并且它的原點被假設(shè)為系統(tǒng)(13)的平衡點,即f(0;ρ)=0。

    定義1.全局有限時間穩(wěn)定性[16]:如果系統(tǒng)(13)的原點是全局漸近穩(wěn)定的,且系統(tǒng)(13)的任何解x(t,x0)在某個有限時刻達到平衡點,即?t≥T(x0)=0,則系統(tǒng)(13)的平衡點被認為是全局有限時間穩(wěn)定的,其中T:Rn→R+?{0}被稱為穩(wěn)定時間函數(shù)。

    定義2.固定時間穩(wěn)定性[8]:如果系統(tǒng)(13)的原點是全局有限時間穩(wěn)定的且穩(wěn)定時間函數(shù)有界,即?Tmax>0:?x0∈Rn,T(x0)≤Tmax,那么系統(tǒng)(13)的平衡點是固定時間穩(wěn)定的。

    定義3.預(yù)定義時間穩(wěn)定性[17]:對于系統(tǒng)(13)的系統(tǒng)參數(shù)ρ和常數(shù)Tc=Tc(ρ)>0,如果系統(tǒng)(13)的平衡點是固定時間穩(wěn)定的,且穩(wěn)定時間函數(shù)T:Rn→R 滿足以下條件:

    則系統(tǒng)(13)的平衡點被認為是預(yù)定義時間穩(wěn)定的,且Tc叫做預(yù)定義時間。此外,如果穩(wěn)定時間函數(shù)滿足,那 么Tc被稱作 強預(yù)定義時間,否則,Tc被稱作弱預(yù)定義時間。

    理論1.[13]如果存在滿足以下條件的連續(xù)正定徑向無界函數(shù)V:Rn→R:

    其中:x∈Rn{0};α,β,p,q,k為滿足kp<1,kq>1 的正常數(shù)。那么,系統(tǒng)(13)的平衡點是預(yù)定義時間穩(wěn)定的,預(yù)定義時間為Tc。如果設(shè)參數(shù)向量,則預(yù)定義時間函數(shù)滿足Tf=γ(ρ),其中γ(ρ)可以通過下式計算得出:

    4 預(yù)定義時間滑??刂破髟O(shè)計與穩(wěn)定性分析

    4.1 預(yù)定義時間滑??刂破髟O(shè)計

    首先,設(shè)計如下滑模變量:

    其中:ωe和ev分別是航天器的角速度跟蹤誤差和姿態(tài)跟蹤誤差;α,β,p,q分別是定義的正常數(shù)且滿足01;γ,mp,mq分別是滿足γ=關(guān)系的正常數(shù);Tc1是定義的滑動階段的預(yù)定義時間常數(shù)

    接著,定義如下的預(yù)定義時間滑模面:

    根據(jù)上述的滑模變量(17)和滑模面(18),所提出的PTSMC 可以表示如下:

    其中:Tc2表示系統(tǒng)到達階段的預(yù)定義時間常數(shù);b1,b1,b2是定義的正常數(shù)并滿足不等式(12),對于向量x∈Rn,|x|表示分別對向量的每個元素取絕對值。

    圖2 給出了(17)~(19)中所提出的PTSMC的控制信號流程圖,其中u1是使系統(tǒng)在預(yù)定義時間內(nèi)收斂的滑模控制器,u2是用于抵抗系統(tǒng)不確定性和外界擾動的魯棒控制器。

    圖2 PTSMC 方案的控制信號流程圖Fig.2 Block diagram showing the flow of the control sig?nals for the PTSMC scheme.

    4.2 穩(wěn)定性分析

    對于所提出的PTSMC 的穩(wěn)定性分析可以分為到達階段和滑動階段進行分析。同時,為了簡化證明過程,我們將式(10)重寫為以下形式:

    4.2.1 系統(tǒng)到達階段的穩(wěn)定性分析

    所提出的滑模面(18)關(guān)于時間的導數(shù)為:

    帶入式(20),可以得到:

    將式(22)帶入所提出的控制輸入(19),有:

    對式(23)進一步化簡可以得到:

    考慮一個Lyapunov 候選函數(shù)為V1(s)=|s|,則它的一階導數(shù)為:

    根據(jù)理論1,我們可以得到,系統(tǒng)在到達階段是預(yù)定義時間穩(wěn)定的,且預(yù)定義時間為Tc2。

    4.2.2 系統(tǒng)滑動階段的穩(wěn)定性分析

    一旦系統(tǒng)的跟蹤軌跡被約束到滑模面上,即當t>Tc2時s=0,系統(tǒng)進入滑動階段,根據(jù)式(17)和(18),可以得到:

    設(shè)滑動階段的Lyapunov 候選函數(shù)為V2(ev)=|ev|,于是可以得到:

    同樣地,根據(jù)理論1,系統(tǒng)在預(yù)定義時間Tc2內(nèi)沿滑模面收斂到平衡點。

    根據(jù)上述證明,我們可以得到,所提出的控制器可以在系統(tǒng)不確定和存在外界擾動的情況下實現(xiàn)剛性航天器的預(yù)定義時間收斂,其收斂時間為Tc=Tc1+Tc2。

    備注1.由于存在系統(tǒng)的耦合不確定性ρ,即使系統(tǒng)的初始誤差ev→∞,系統(tǒng)誤差也會在預(yù)定義時間Tc之前收斂到原點,即Tc為弱預(yù)定義時間。同時,我們可以考慮一種特殊情況,當k0=成立時,式(25)和(26)中的等號成立,此時Tc為系統(tǒng)的最小收斂時間上界,根據(jù)定義3,Tc為強預(yù)定義時間。相較于現(xiàn)有的預(yù)定義時間控制方法,如文獻[5-9]中所提出的控制方案,本文中的Lyapunov 函數(shù)的導數(shù)僅在式(25)考慮耦合不確定時進行了縮放,而其他控制方案則在推導過程中進行了多次縮放,因此本文中所提出收斂時間上界是更不保守的。

    備注2.考慮到符號函數(shù)sign(·)會導致系統(tǒng)控制輸入的抖振,為了改善符號函數(shù)導致的系統(tǒng)抖振,將所提出的控制器(19)中的sign(s)修改為同時,將式(19)中的函數(shù)a修改為如下形式:

    其中,ε∈R 是一個很小的正常數(shù)。

    5 實驗驗證與分析

    在本章節(jié)中,為了證明所提出控制方案在系統(tǒng)收斂時間的保守性和跟蹤誤差精度等方面的優(yōu)越性,分別與現(xiàn)有的預(yù)定義時間控制方案,NTSMC 方案和傳統(tǒng)PD 控制方案進行仿真對比。所有的數(shù)值仿真實驗均基于Simulink/MATLAB 2020a 進行,基本采樣時間設(shè)為10-4s。

    與文獻[7]中給出的航天器系統(tǒng)的參數(shù)相同,假設(shè)名義慣性矩陣為J0=[201.20.9;1.217 1.4;0.91.415]kg·m2,慣性矩陣的不確定度設(shè)為ΔJ=diag(2,2,3)kg·m2。期望的角速度設(shè)為ωd=0.05[sin(πt/100)sin(2πt/100)sin(3πt/100)]rad/s,外界擾動設(shè)為d=[0.1sin(t)0.2 sin(1.2t)0.3sin(1.5t)]N·m。系統(tǒng)的初始狀態(tài)設(shè)為q0=[0.3-0.2-0.3 0.883 2]Trad 和ω0=[0.06-0.04 0.05]Trad/s。對于所提出的PTSMC 方案,其控制參數(shù)分別設(shè)為p=1,q=3,k=0.5,α=4,β=0.25,Tc1=Tc2=3,ε=10?4,b0=1.8,b1=17.2,b2=3。

    5.1 與現(xiàn)有的預(yù)定義時間控制方案的對比仿真實驗

    考慮文獻[14]中所提出的二階預(yù)定義時間滑??刂破鳎⊿econd Order Predefined-Time Slid?ing Mode Control,SOPTSMC):

    為了公平比較,修改控制器(30)中的k為本文中所提出的k0,sign(σ)修改為本文中在備注2中所設(shè)定的其他參數(shù)參考文獻[14]中給出的相應(yīng)參數(shù)值,即q1=1.2,p1=1,q2=0.5,p2=1,α2=10?3,β2=1 將其預(yù)定義參數(shù)設(shè)為Tc1=Tc2=3。

    對比仿真結(jié)果如圖3~5 所示。圖3 和圖4 分別表示與SOPTSMC 控制方案相比的姿態(tài)跟蹤誤差和角速度跟蹤誤差,從圖中可以看出,給定相同的預(yù)定義收斂時間Tc=Tc1+Tc2=6 s,所提出控制方案中的跟蹤誤差在4 s 后收斂到零,而SOPTSMC 方案的實際收斂時間小于3 s。一般而言,我們總是希望航天器系統(tǒng)跟蹤誤差的實際收斂時間更接近于所設(shè)定的預(yù)定義收斂時間,這將有利于系統(tǒng)收斂時間的合理設(shè)置,因為更快的收斂時間往往意味著更多的能量消耗。圖5 顯示了控制力矩的對比,從圖中可以看出,所提出控制方案的初始力矩約為50 N,遠小于SOPTSMC 方案中的初始輸出力矩。此外,所提出的控制方案的控制力矩更加平滑,抖振更小,這對于剛性航天器的姿態(tài)跟蹤控制是更有利的。

    圖3 PTSMC 與SOPTSMC 的姿態(tài)跟蹤誤差Fig.3 Attitude tracking errors of PTSMC and SOPTSMC

    圖4 PTSMC 與SOPTSMC 的角速度跟蹤誤差Fig.4 Angular velocity tracking errors of PTSMC and SOPTSMC

    圖5 PTSMC 與SOPTSMC 的控制輸入Fig.5 Control inputs of PTSMC and SOPTSMC

    5.2 與NTSMC 方案和傳統(tǒng)PD 控制方案的對比仿真實驗

    參考文獻[7]中的所提出的NTSMC 方案

    其中:控制參數(shù)滿足0

    參考文獻[18]中給出的PD 控制方案:

    其中,kp和kd為給定的控制參數(shù),分別定義kp=6I3,kd=8I3。

    圖6 和圖7 分別顯示了所提出的控制方案與PD 控制和NTSMC 方案的姿態(tài)跟蹤誤差和角速度誤差,從圖中可以看出,PD 控制中的姿態(tài)跟蹤誤差漸進地收斂到零點。同時,由于PD 控制對系統(tǒng)的耦合不確定性不具有魯棒性,因此它有較差的誤差跟蹤精度。相較于NTSMC 方案,所提出的控制方案中系統(tǒng)的跟蹤誤差精度要遠高于NTSMC 方案,其姿態(tài)跟蹤誤差精度高于1.5×10-6rad,角速度跟蹤誤差精度高于2×10-6rad/s。最重要的是,PD 控制的跟蹤誤差是漸進收斂的,它的收斂時間趨近于無窮大,NTSMC 是有限時間收斂的,但它的收斂時間取決于系統(tǒng)的初始狀態(tài)。而本文中所提出的PTSMC 方案的收斂時間可以通過控制器中的預(yù)定義參數(shù)直接進行調(diào)節(jié),對于航天器的姿態(tài)控制而言,這有利于它在規(guī)定的時間內(nèi)完成姿態(tài)調(diào)整。

    圖6 PTSMC,PD 和NTSMC 的姿態(tài)跟蹤誤差Fig.6 Attitude tracking errors of PTSMC,PD and NTSMC

    圖7 PTSMC,PD 和NTSMC 的角速度誤差Fig.7 Angular velocity errors of PTSMC,PD and NTSMC

    總之,仿真結(jié)果表明,本文所提出的PTSMC方案相較于目前存在的預(yù)定義時間控制方案,其收斂時間的上界更加不保守。相較于現(xiàn)有的NTSMC 和PD 控制方案,有更高的跟蹤精度和魯棒性,并具有預(yù)定義時間收斂的特性。同時,在所提出的控制器中只需要采集航天器的實際姿態(tài)與角速度,這對于航天器的姿態(tài)控制是可實現(xiàn)的,可以被應(yīng)用在航天器的實際姿態(tài)控制中。

    5.3 基于三自由度氣浮平臺的實驗驗證

    由于在地面難以開展空間狀態(tài)的飛行器姿態(tài)跟蹤實驗,為了進一步驗證所提出的控制方案在實際應(yīng)用中的有效性,本章節(jié)利用本課題組研制的三自由度自由飛行機器人氣浮式平臺[19]模擬飛行器姿態(tài)控制,進行了多次姿態(tài)跟蹤控制實驗,實驗裝置如圖8 所示。該裝置的運動是一個剛體的平面運動,包括平面x,y軸方向的移動,和一個繞自身的轉(zhuǎn)動量?。

    圖8 三自由度自由飛行機器人氣浮平臺Fig.8 3-DOF free-flying robot floating platform

    為了測試文中所提出的控制方案,分別進行了四次關(guān)于直線軌跡的跟蹤實驗,同時,在直線軌跡的跟蹤過程中伴隨平臺角度的跟蹤。四次實驗的軌跡跟蹤結(jié)果如圖9 所示,實驗均在氣浮平臺所在的大理石平臺上進行。可以看出,四次實驗均獲得了較高的位置跟蹤精度。對于軌跡1,在跟蹤軌跡末端位置誤差變大,這是由于軌跡1 中氣浮平臺運行時間較長,平臺后期供氣不足導致的。對于軌跡1 和軌跡2,氣浮平臺的初始位置與命令軌跡均有較大誤差,但氣浮平臺在文中所提出的控制方案的作用下,均能迅速地跟蹤到指令軌跡。

    圖9 三自由度自由飛行機器人氣浮平臺的姿態(tài)跟蹤實驗Fig.9 Attitude tracking experimental results of the float?ing platform of a 3-DOF free-flying robot

    進一步地,在四次跟蹤實驗中,氣浮平臺的關(guān)于x,y軸和轉(zhuǎn)動量?隨時間的指令變化量與實際變化量如圖10 所示。在每次實驗中,控制器的預(yù)定義時間被設(shè)為20 s,從圖中可以看出,各自由度的姿態(tài)跟蹤誤差大約均在10 s 內(nèi)收斂,這說明了所提出的預(yù)定義時間滑??刂破鞯挠行?。同時,從實驗結(jié)果可以看出,氣浮平臺具有優(yōu)異的姿態(tài)跟蹤性能,其中,角度跟蹤誤差小于0.1 rad,x和y軸方向的跟蹤誤差小于0.2 m。

    圖10 氣浮平臺的x,y 和? 隨時間的變化Fig.10 Changes of x,y and ? of the floating platform with time

    6 結(jié)論

    本文在保證性能的前提下,解決了航天器的預(yù)定義時間姿態(tài)跟蹤問題。通過Lyapunov 控制理論,我們嚴格證明了所提出的控制方案對于航天器的姿態(tài)跟蹤的調(diào)整時間是可以通過控制參數(shù)被指定的。同時,相較于現(xiàn)有的預(yù)定義時間控制策略,文中所提出的控制方案的收斂時間是更不保守的,這有利于航天器姿態(tài)調(diào)整時間的合理設(shè)置。最后,通過仿真實驗證明剛性航天器的姿態(tài)跟蹤誤差精度可達1.5×10-6rad,角速度跟蹤誤差精度可達2×10-6rad/s,同時也說明了本方案相較于現(xiàn)有控制方案的優(yōu)越性,包括更不保守的收斂時間上界,良好的抗干擾性和較高的跟蹤誤差精度。通過3 自由度氣浮平臺的姿態(tài)跟蹤實驗說明了控制方案的有效性,角度跟蹤誤差小于0.1 rad,位置跟蹤誤差小于0.2 m。

    猜你喜歡
    角速度航天器滑模
    2022 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2022年7期)2022-08-16 09:52:50
    2019 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2019年9期)2019-10-23 01:55:34
    基于組合滑??刂频慕^對重力儀兩級主動減振設(shè)計
    2018 年第三季度航天器發(fā)射統(tǒng)計
    國際太空(2018年12期)2019-01-28 12:53:20
    測控技術(shù)(2018年4期)2018-11-25 09:47:26
    并網(wǎng)逆變器逆系統(tǒng)自學習滑??箶_控制
    2018年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2018年9期)2018-10-18 08:51:32
    圓周運動角速度測量方法賞析
    半捷聯(lián)雷達導引頭視線角速度提取
    基于構(gòu)架點頭角速度的軌道垂向長波不平順在線檢測

    3wmmmm亚洲av在线观看| 精品久久久久久久人妻蜜臀av| 国产成+人综合+亚洲专区| 1000部很黄的大片| 天堂动漫精品| 18美女黄网站色大片免费观看| 亚洲一区高清亚洲精品| 天天一区二区日本电影三级| 美女大奶头视频| 国产精品日韩av在线免费观看| 成年人黄色毛片网站| 国内久久婷婷六月综合欲色啪| 国产亚洲精品综合一区在线观看| 好男人电影高清在线观看| 国产免费av片在线观看野外av| 床上黄色一级片| 国产亚洲精品久久久com| 亚洲avbb在线观看| 好看av亚洲va欧美ⅴa在| 亚洲专区中文字幕在线| АⅤ资源中文在线天堂| 精品熟女少妇八av免费久了| 日韩欧美在线乱码| 亚洲美女黄片视频| 久久6这里有精品| 国产97色在线日韩免费| 一本综合久久免费| 国产三级中文精品| a级毛片a级免费在线| 国产精华一区二区三区| 夜夜躁狠狠躁天天躁| 亚洲人成网站在线播放欧美日韩| 欧美黄色片欧美黄色片| 欧美色视频一区免费| 国产高清三级在线| 12—13女人毛片做爰片一| 岛国视频午夜一区免费看| 国产精品久久视频播放| 亚洲在线自拍视频| 午夜福利高清视频| 俄罗斯特黄特色一大片| 国产精品久久久久久久久免 | 婷婷丁香在线五月| 成年版毛片免费区| 欧美午夜高清在线| 亚洲精华国产精华精| 亚洲男人的天堂狠狠| 成人高潮视频无遮挡免费网站| 亚洲成人精品中文字幕电影| 偷拍熟女少妇极品色| 99热6这里只有精品| 国内精品久久久久精免费| 欧美三级亚洲精品| 一二三四社区在线视频社区8| 桃红色精品国产亚洲av| 一级毛片女人18水好多| 欧美又色又爽又黄视频| 91在线精品国自产拍蜜月 | 国产伦精品一区二区三区视频9 | 国内精品一区二区在线观看| 九九久久精品国产亚洲av麻豆| 久久久久久久精品吃奶| 人妻久久中文字幕网| 天天一区二区日本电影三级| 午夜免费观看网址| 99精品久久久久人妻精品| 国产在线精品亚洲第一网站| 日本免费一区二区三区高清不卡| 99国产综合亚洲精品| 日本撒尿小便嘘嘘汇集6| 9191精品国产免费久久| 伊人久久大香线蕉亚洲五| 久久精品综合一区二区三区| 欧美中文日本在线观看视频| 欧美最新免费一区二区三区 | 亚洲精品粉嫩美女一区| 在线观看日韩欧美| 无限看片的www在线观看| 日韩欧美精品v在线| 欧美日韩乱码在线| 色尼玛亚洲综合影院| 麻豆成人av在线观看| 久久婷婷人人爽人人干人人爱| 母亲3免费完整高清在线观看| 亚洲不卡免费看| av天堂中文字幕网| 国产国拍精品亚洲av在线观看 | 国产亚洲精品久久久久久毛片| 国产欧美日韩一区二区三| 性色avwww在线观看| 观看免费一级毛片| 午夜福利在线在线| 一夜夜www| 欧美另类亚洲清纯唯美| 亚洲成人免费电影在线观看| 色综合站精品国产| xxxwww97欧美| 国产精品99久久99久久久不卡| 婷婷六月久久综合丁香| 国产欧美日韩一区二区精品| 好看av亚洲va欧美ⅴa在| 五月伊人婷婷丁香| 男人的好看免费观看在线视频| 国产精品一区二区免费欧美| 欧美日韩黄片免| 亚洲精品国产精品久久久不卡| 国产午夜福利久久久久久| 男女床上黄色一级片免费看| 日本黄色片子视频| 一区二区三区激情视频| 亚洲性夜色夜夜综合| 欧美中文日本在线观看视频| 亚洲色图av天堂| 国产精品,欧美在线| 午夜免费激情av| 日本在线视频免费播放| 大型黄色视频在线免费观看| 亚洲欧美一区二区三区黑人| 国产久久久一区二区三区| 在线天堂最新版资源| 中文字幕人妻丝袜一区二区| 色视频www国产| 国产毛片a区久久久久| 在线观看一区二区三区| 好男人在线观看高清免费视频| 少妇的逼水好多| 美女大奶头视频| 欧美成人一区二区免费高清观看| 少妇的逼水好多| 少妇人妻一区二区三区视频| 国产精品综合久久久久久久免费| 波多野结衣高清作品| 午夜激情欧美在线| 日韩精品中文字幕看吧| 国产欧美日韩精品亚洲av| 老司机午夜十八禁免费视频| 九九在线视频观看精品| 国产精品爽爽va在线观看网站| 最新在线观看一区二区三区| 国产精品98久久久久久宅男小说| 国产精品,欧美在线| 在线观看免费午夜福利视频| 成人国产综合亚洲| 床上黄色一级片| 欧美丝袜亚洲另类 | 国产精品电影一区二区三区| 美女大奶头视频| 在线观看午夜福利视频| 欧美日本视频| 岛国视频午夜一区免费看| 无人区码免费观看不卡| 欧美日韩精品网址| 90打野战视频偷拍视频| 丰满人妻熟妇乱又伦精品不卡| 两个人的视频大全免费| 欧美激情在线99| 制服人妻中文乱码| 国产三级在线视频| 国产亚洲av嫩草精品影院| 国产激情偷乱视频一区二区| 日韩欧美精品v在线| 免费看十八禁软件| 看片在线看免费视频| 国产野战对白在线观看| 国产精华一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国产伦精品一区二区三区视频9 | 日韩中文字幕欧美一区二区| 免费看美女性在线毛片视频| 淫妇啪啪啪对白视频| 日韩欧美国产一区二区入口| 俺也久久电影网| 天堂av国产一区二区熟女人妻| 美女大奶头视频| 桃红色精品国产亚洲av| 成人高潮视频无遮挡免费网站| 日韩欧美三级三区| 老司机午夜十八禁免费视频| 夜夜看夜夜爽夜夜摸| 欧美性猛交黑人性爽| 欧美极品一区二区三区四区| 日韩免费av在线播放| 性色av乱码一区二区三区2| 丰满人妻熟妇乱又伦精品不卡| 午夜福利高清视频| 99riav亚洲国产免费| 午夜日韩欧美国产| 五月玫瑰六月丁香| 日本三级黄在线观看| 1000部很黄的大片| 中国美女看黄片| 久久久久久久久中文| 欧美bdsm另类| 欧美色欧美亚洲另类二区| 变态另类成人亚洲欧美熟女| 观看美女的网站| 亚洲天堂国产精品一区在线| 国产一区在线观看成人免费| 美女 人体艺术 gogo| 一个人看的www免费观看视频| 长腿黑丝高跟| 日本三级黄在线观看| 日本免费一区二区三区高清不卡| 高清毛片免费观看视频网站| 亚洲最大成人中文| 久久久国产精品麻豆| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人精品中文字幕电影| 操出白浆在线播放| 久久精品91蜜桃| 国内毛片毛片毛片毛片毛片| 人妻丰满熟妇av一区二区三区| 97人妻精品一区二区三区麻豆| 国产精品嫩草影院av在线观看 | 欧美极品一区二区三区四区| 精品人妻1区二区| 性色avwww在线观看| 91九色精品人成在线观看| 久久久久久久久大av| aaaaa片日本免费| 欧美区成人在线视频| 国产精品精品国产色婷婷| 观看免费一级毛片| 国内揄拍国产精品人妻在线| 99久久精品国产亚洲精品| 小说图片视频综合网站| 精品免费久久久久久久清纯| 午夜免费男女啪啪视频观看 | 精品久久久久久久人妻蜜臀av| 热99在线观看视频| 一级黄片播放器| 国产视频一区二区在线看| 偷拍熟女少妇极品色| 淫妇啪啪啪对白视频| 3wmmmm亚洲av在线观看| 国产亚洲欧美98| 一区二区三区激情视频| 天天躁日日操中文字幕| 伊人久久精品亚洲午夜| 欧美性猛交╳xxx乱大交人| 国产精品一区二区免费欧美| 99久久精品热视频| 亚洲第一欧美日韩一区二区三区| 人人妻,人人澡人人爽秒播| 黄片大片在线免费观看| 亚洲色图av天堂| 嫩草影院精品99| 国产精品一区二区免费欧美| 免费av毛片视频| 欧美中文日本在线观看视频| 99国产极品粉嫩在线观看| 黄片大片在线免费观看| 成人国产综合亚洲| 亚洲av二区三区四区| 乱人视频在线观看| 99精品在免费线老司机午夜| 18禁国产床啪视频网站| 日韩欧美在线二视频| 欧美又色又爽又黄视频| 亚洲人成电影免费在线| 精品国产超薄肉色丝袜足j| av视频在线观看入口| 国产伦精品一区二区三区四那| 婷婷六月久久综合丁香| 成人亚洲精品av一区二区| 精品电影一区二区在线| 久久久久久久精品吃奶| 日韩欧美在线乱码| 亚洲精品国产精品久久久不卡| 黄色丝袜av网址大全| 亚洲精品一区av在线观看| 最近在线观看免费完整版| 国产精品1区2区在线观看.| 亚洲av电影在线进入| 九色国产91popny在线| 丁香欧美五月| 色视频www国产| 欧美一区二区精品小视频在线| 一个人免费在线观看电影| 国产伦精品一区二区三区四那| 国产91精品成人一区二区三区| 婷婷丁香在线五月| 97人妻精品一区二区三区麻豆| 国产成人av教育| 国产亚洲精品久久久久久毛片| 日日干狠狠操夜夜爽| 小蜜桃在线观看免费完整版高清| 长腿黑丝高跟| 午夜久久久久精精品| 丁香欧美五月| 亚洲精品456在线播放app | 三级国产精品欧美在线观看| 午夜老司机福利剧场| 国产精品乱码一区二三区的特点| 免费观看人在逋| 中文字幕高清在线视频| 国产一区二区在线av高清观看| 国产 一区 欧美 日韩| 国产一区在线观看成人免费| 99久久综合精品五月天人人| 91av网一区二区| 国产日本99.免费观看| 91久久精品电影网| 精品熟女少妇八av免费久了| 国产精品综合久久久久久久免费| 麻豆成人午夜福利视频| 免费在线观看影片大全网站| 毛片女人毛片| 国产亚洲精品久久久久久毛片| 免费搜索国产男女视频| 欧美日韩精品网址| 亚洲精品影视一区二区三区av| 757午夜福利合集在线观看| 亚洲av电影在线进入| 成人一区二区视频在线观看| 国产高清视频在线观看网站| 国产高清三级在线| 国产成年人精品一区二区| 亚洲国产精品久久男人天堂| 搡女人真爽免费视频火全软件 | 男女午夜视频在线观看| 成年免费大片在线观看| 国产av在哪里看| 精品熟女少妇八av免费久了| 精品电影一区二区在线| 无人区码免费观看不卡| 欧美又色又爽又黄视频| xxxwww97欧美| 久久人人精品亚洲av| 国产美女午夜福利| 日韩免费av在线播放| 真人一进一出gif抽搐免费| 国产精品亚洲一级av第二区| 亚洲五月婷婷丁香| 亚洲中文字幕一区二区三区有码在线看| 一边摸一边抽搐一进一小说| 女同久久另类99精品国产91| 男人的好看免费观看在线视频| 在线a可以看的网站| av视频在线观看入口| 久久精品国产亚洲av涩爱 | 精品午夜福利视频在线观看一区| 少妇人妻一区二区三区视频| 成年女人毛片免费观看观看9| 婷婷精品国产亚洲av在线| 欧美日韩福利视频一区二区| 欧美+亚洲+日韩+国产| 精品一区二区三区视频在线 | 成年免费大片在线观看| 神马国产精品三级电影在线观看| 亚洲成a人片在线一区二区| 久久久久亚洲av毛片大全| 亚洲人成电影免费在线| 97碰自拍视频| 美女高潮的动态| 国产男靠女视频免费网站| 亚洲人成伊人成综合网2020| 国产av在哪里看| 97超级碰碰碰精品色视频在线观看| 亚洲专区中文字幕在线| 国内精品久久久久久久电影| 亚洲自拍偷在线| 欧美丝袜亚洲另类 | 国内少妇人妻偷人精品xxx网站| 久9热在线精品视频| 午夜福利在线观看免费完整高清在 | 亚洲精品日韩av片在线观看 | 亚洲精品亚洲一区二区| 在线看三级毛片| 此物有八面人人有两片| 国产97色在线日韩免费| 精品久久久久久,| 老司机午夜十八禁免费视频| 女人高潮潮喷娇喘18禁视频| 欧美丝袜亚洲另类 | 亚洲精品亚洲一区二区| 97超级碰碰碰精品色视频在线观看| 婷婷精品国产亚洲av在线| 日本撒尿小便嘘嘘汇集6| 2021天堂中文幕一二区在线观| e午夜精品久久久久久久| 国产野战对白在线观看| 亚洲av不卡在线观看| av国产免费在线观看| 特大巨黑吊av在线直播| 精品国产三级普通话版| 精品人妻一区二区三区麻豆 | 内地一区二区视频在线| 久久精品国产清高在天天线| 日本一本二区三区精品| 国产一区二区激情短视频| 国产伦精品一区二区三区四那| 国产精品日韩av在线免费观看| 久久久久久久久中文| 舔av片在线| 欧美日本视频| 女人十人毛片免费观看3o分钟| 欧美不卡视频在线免费观看| www.www免费av| 麻豆成人av在线观看| 男人的好看免费观看在线视频| 国产极品精品免费视频能看的| tocl精华| 午夜福利在线在线| 特大巨黑吊av在线直播| 亚洲激情在线av| 最近最新中文字幕大全电影3| 免费电影在线观看免费观看| 成人欧美大片| 噜噜噜噜噜久久久久久91| 亚洲成av人片在线播放无| 国产精品女同一区二区软件 | 欧美日本视频| 亚洲欧美日韩卡通动漫| 午夜两性在线视频| 熟女人妻精品中文字幕| 日韩欧美精品v在线| 色av中文字幕| 亚洲,欧美精品.| 亚洲欧美日韩高清在线视频| 丰满乱子伦码专区| 又黄又粗又硬又大视频| 国产精品 国内视频| 夜夜看夜夜爽夜夜摸| 亚洲国产色片| 人人妻人人澡欧美一区二区| 久久久久久久久中文| 变态另类丝袜制服| 一本综合久久免费| 亚洲精品影视一区二区三区av| 又紧又爽又黄一区二区| 午夜精品在线福利| 婷婷精品国产亚洲av在线| 免费看日本二区| 草草在线视频免费看| 国产高清videossex| 国产欧美日韩精品亚洲av| 日日摸夜夜添夜夜添小说| 欧美性猛交╳xxx乱大交人| 国产午夜精品论理片| 免费观看人在逋| 国产乱人伦免费视频| a在线观看视频网站| 999久久久精品免费观看国产| 在线观看免费午夜福利视频| 99国产精品一区二区三区| 18禁美女被吸乳视频| 国语自产精品视频在线第100页| 亚洲av熟女| 久久国产精品影院| 成人性生交大片免费视频hd| 日韩成人在线观看一区二区三区| 国产精品永久免费网站| 一二三四社区在线视频社区8| 国产激情偷乱视频一区二区| 亚洲av第一区精品v没综合| 亚洲专区国产一区二区| 精品熟女少妇八av免费久了| 村上凉子中文字幕在线| 99久久成人亚洲精品观看| 99精品欧美一区二区三区四区| 看黄色毛片网站| 制服丝袜大香蕉在线| 免费无遮挡裸体视频| 一边摸一边抽搐一进一小说| 国产乱人视频| 人妻久久中文字幕网| 久久中文看片网| 国产一区二区在线观看日韩 | 色在线成人网| 婷婷丁香在线五月| 国产成年人精品一区二区| 看黄色毛片网站| 99国产精品一区二区蜜桃av| 亚洲国产高清在线一区二区三| 国产麻豆成人av免费视频| 亚洲成人中文字幕在线播放| 国产在线精品亚洲第一网站| av在线天堂中文字幕| 国产毛片a区久久久久| 久久性视频一级片| 国产三级在线视频| 级片在线观看| a级毛片a级免费在线| 欧美激情久久久久久爽电影| 青草久久国产| 国产精品影院久久| 法律面前人人平等表现在哪些方面| 国产高清视频在线观看网站| 国产伦一二天堂av在线观看| a在线观看视频网站| 欧美+亚洲+日韩+国产| 人妻久久中文字幕网| 99精品欧美一区二区三区四区| 国产亚洲精品av在线| 老鸭窝网址在线观看| 每晚都被弄得嗷嗷叫到高潮| 99久久精品国产亚洲精品| 国产精品av视频在线免费观看| 国产探花在线观看一区二区| 天美传媒精品一区二区| 夜夜看夜夜爽夜夜摸| 免费在线观看成人毛片| 国产高清激情床上av| 午夜精品一区二区三区免费看| 国内毛片毛片毛片毛片毛片| 亚洲中文字幕日韩| 色视频www国产| 久久久久久久久大av| 亚洲中文日韩欧美视频| 天天躁日日操中文字幕| 成人国产一区最新在线观看| 亚洲不卡免费看| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩东京热| 亚洲av免费高清在线观看| 亚洲自拍偷在线| 国产高清videossex| 淫秽高清视频在线观看| 国语自产精品视频在线第100页| 精品熟女少妇八av免费久了| avwww免费| 日日摸夜夜添夜夜添小说| 精品免费久久久久久久清纯| 人妻夜夜爽99麻豆av| 欧美一级毛片孕妇| 天堂√8在线中文| 亚洲欧美日韩东京热| 真人一进一出gif抽搐免费| 成人三级黄色视频| 中文字幕久久专区| 欧美在线一区亚洲| 中文字幕av成人在线电影| 法律面前人人平等表现在哪些方面| 国产一区二区在线av高清观看| 成人国产综合亚洲| 一进一出抽搐动态| 国产精品久久久人人做人人爽| 午夜日韩欧美国产| 18美女黄网站色大片免费观看| 波多野结衣高清无吗| av黄色大香蕉| 偷拍熟女少妇极品色| 亚洲 欧美 日韩 在线 免费| 亚洲精品456在线播放app | 免费电影在线观看免费观看| 欧美日韩国产亚洲二区| 午夜精品在线福利| 日韩 欧美 亚洲 中文字幕| 成人欧美大片| 亚洲五月天丁香| 午夜福利在线在线| 久久人妻av系列| 日本免费a在线| 亚洲中文字幕一区二区三区有码在线看| 美女高潮的动态| 亚洲国产欧洲综合997久久,| 久久久久国产精品人妻aⅴ院| 在线观看日韩欧美| 99久久精品热视频| 亚洲无线观看免费| 毛片女人毛片| 久久精品国产自在天天线| 亚洲av中文字字幕乱码综合| 一边摸一边抽搐一进一小说| 亚洲乱码一区二区免费版| 天堂网av新在线| 日韩av在线大香蕉| 国产精品久久久久久亚洲av鲁大| 一本精品99久久精品77| 日本免费a在线| 国产高潮美女av| 久久伊人香网站| av在线天堂中文字幕| 日本黄大片高清| 亚洲人成网站在线播| 麻豆国产97在线/欧美| 免费av不卡在线播放| 国产亚洲欧美98| 长腿黑丝高跟| 国产精品1区2区在线观看.| av片东京热男人的天堂| www.999成人在线观看| 白带黄色成豆腐渣| xxxwww97欧美| 免费高清视频大片| 99久久九九国产精品国产免费| 午夜激情欧美在线| 亚洲av免费在线观看| 在线十欧美十亚洲十日本专区| bbb黄色大片| 在线免费观看的www视频| 国内精品一区二区在线观看| 午夜激情欧美在线| 88av欧美| 成年版毛片免费区| 色视频www国产| 变态另类丝袜制服| 欧美日韩综合久久久久久 | 久久亚洲真实| 午夜福利在线在线| 国内精品久久久久精免费| 亚洲五月天丁香| 亚洲欧美日韩无卡精品| 国产精品野战在线观看| 高清毛片免费观看视频网站| eeuss影院久久| 在线播放无遮挡| 精品一区二区三区视频在线 | 中文字幕av在线有码专区| 美女cb高潮喷水在线观看| 男人舔女人下体高潮全视频| 色综合婷婷激情| 精品欧美国产一区二区三| 国产免费av片在线观看野外av|