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

    改進(jìn)規(guī)則下自適應(yīng)神經(jīng)網(wǎng)絡(luò)的UUV 水平面路徑跟隨控制

    2014-04-13 04:18:40周佳加嚴(yán)浙平賈鶴鳴侯恕萍
    關(guān)鍵詞:水平面航向制導(dǎo)

    周佳加,嚴(yán)浙平,賈鶴鳴,侯恕萍

    (1. 哈爾濱工程大學(xué) 自動(dòng)化學(xué)院,黑龍江 哈爾濱,150001;2. 東北林業(yè)大學(xué) 機(jī)電工程學(xué)院,黑龍江 哈爾濱,150040;3. 哈爾濱工程大學(xué) 機(jī)電工程學(xué)院,黑龍江 哈爾濱,150001)

    無(wú)人水下航行器(unmanned underwater vehicle,UUV)的水平面路徑跟隨能力是實(shí)現(xiàn)其水下管道檢測(cè)、海底地形勘測(cè)和打撈施工的重要技術(shù)前提[1-4]。實(shí)際中,待檢測(cè)海底管道、電纜及海底地形等靜態(tài)目標(biāo)常常認(rèn)為可用一系列直線段的集合來(lái)近似,此時(shí),相鄰直線段的交角可能為銳角、直角或者鈍角。這意味著當(dāng)UUV 在相鄰直線段運(yùn)動(dòng)時(shí),跟隨目標(biāo)將發(fā)生跳變,系統(tǒng)狀態(tài)也極有可能朝著發(fā)散的方向發(fā)展。為了保證檢測(cè)的可靠性,要求攜帶觀測(cè)設(shè)備的載體在指定高度、準(zhǔn)確的水平面位置航行,因此,減速轉(zhuǎn)彎是實(shí)現(xiàn)準(zhǔn)確路徑跟隨的措施之一,其可有效地減小轉(zhuǎn)彎時(shí)航向的超調(diào)。但是,當(dāng)UUV 浮力與重力不均衡時(shí),UUV 一旦減速,隨著舵效的降低,其深度將會(huì)有明顯變化,即深度控制能力降低。另外,在不同的航行環(huán)境中[5-6],UUV 的水動(dòng)力參數(shù)是變化的,很難精確地估計(jì)水動(dòng)力參數(shù)的先驗(yàn)值,加之外界海流干擾的影響,使得UUV在航道改變過(guò)程中既滿足深度或高度控制精度,又降低轉(zhuǎn)彎引起的水平面位置誤差成為亟需解決的工程難題。目前,在解決UUV 的運(yùn)動(dòng)控制中,大部分方法集中在控制器設(shè)計(jì)方面去求解。畢鳳陽(yáng)等[7]為了實(shí)現(xiàn)具有時(shí)變水動(dòng)力系數(shù)、未建模動(dòng)態(tài)和時(shí)滯等不確定性的欠驅(qū)動(dòng)無(wú)人水下航行器的魯棒控制,以變結(jié)構(gòu)控制的切換函數(shù)及其變化率為模糊控制器的輸入,以變結(jié)構(gòu)控制律的變化率為模糊控制器的輸出設(shè)計(jì)了一個(gè)模糊變結(jié)構(gòu)控制器。但是,該控制器在目標(biāo)跟蹤狀態(tài)發(fā)生變化時(shí)有抖振現(xiàn)象。施淑偉等[8]針對(duì)欠驅(qū)動(dòng)AUV(autonomous underwater vehicle)在常值海流影響下的水平面路徑跟蹤控制問(wèn)題,建立了常值海流作用下的AUV 運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)模型,并采用以自由參考點(diǎn)為原點(diǎn)的Serret-Frenet 坐標(biāo)系描述路徑跟蹤誤差及其動(dòng)態(tài)模型,綜合應(yīng)用李亞普諾夫方法和反步法設(shè)計(jì)了AUV 路徑跟蹤的運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)控制律,水平面跟蹤控制效果較為理想,但是,從其速度變化曲線可發(fā)現(xiàn)在軌跡跟蹤開(kāi)始后的18 s 左右載體的橫向速度大于縱向運(yùn)動(dòng)速度,實(shí)際中這將影響欠驅(qū)動(dòng)無(wú)人水下航行器的水平面運(yùn)動(dòng)穩(wěn)定性,載體容易發(fā)生側(cè)滑而引起航行器失控[8]。上述研究單純從控制器角度設(shè)計(jì)UUV水平面運(yùn)動(dòng)控制器,并未考慮制導(dǎo)規(guī)則在運(yùn)動(dòng)控制問(wèn)題中的重要性。對(duì)UUV 在水平面的運(yùn)動(dòng)分析知,該機(jī)動(dòng)問(wèn)題可分為2 個(gè)子問(wèn)題:一方面為幾何學(xué)問(wèn)題,即利用設(shè)計(jì)的制導(dǎo)律和控制器使UUV 的運(yùn)動(dòng)到預(yù)先規(guī)劃的路徑上;另一方面為動(dòng)力學(xué)問(wèn)題,即使UUV的前向速度保持在設(shè)計(jì)航速。對(duì)于前者,本文利用Line-of-Sight(LOS)系統(tǒng)[9-11]將期望位置(xd, yd)映射到指令航向角,并通過(guò)設(shè)計(jì)的自適應(yīng)神經(jīng)網(wǎng)絡(luò)控制器[12-15](adaptive NN controller,即ANNC)解算控制量實(shí)現(xiàn)載體的水平面路徑跟隨控制,執(zhí)行機(jī)構(gòu)為方向舵。ANNC 采取RBF 神經(jīng)網(wǎng)絡(luò)估計(jì)不確定項(xiàng),選用合適的自適應(yīng)律來(lái)保證神經(jīng)網(wǎng)絡(luò)權(quán)值的最優(yōu)估計(jì)使得估計(jì)更加精確。同時(shí),為了便于設(shè)計(jì)ANCC 控制器,將UUV的水平面動(dòng)力學(xué)特性簡(jiǎn)化為一個(gè)非仿射純反饋系統(tǒng),并通過(guò)恰當(dāng)?shù)倪x擇設(shè)計(jì)參數(shù)保證閉環(huán)系統(tǒng)的控制性能。基于李亞普諾夫穩(wěn)定性定理,該控制器可以保證路徑跟隨控制系統(tǒng)中的誤差狀態(tài)是漸近穩(wěn)定的,使得系統(tǒng)輸出收斂到臨近設(shè)定路徑的很小范圍內(nèi)。對(duì)于后者,利用主推進(jìn)器實(shí)現(xiàn)載體航速控制。

    1 UUV 水平面運(yùn)動(dòng)模型

    由剛體動(dòng)力學(xué)角度分析,UUV 在水平面運(yùn)動(dòng)控制問(wèn)題就是解決其在該平面內(nèi)的三自由度路徑跟隨控制。本文研究的UUV 在水平面機(jī)動(dòng)時(shí)利用主推進(jìn)器來(lái)提供縱向推力,方向舵用于實(shí)現(xiàn)航向控制,即控制輸入維數(shù)少于系統(tǒng)維數(shù),屬于欠驅(qū)動(dòng)UUV。另外,當(dāng)全驅(qū)動(dòng)水下無(wú)人航行器前向運(yùn)動(dòng)速度較高時(shí),橫向推進(jìn)器推力對(duì)載體的作用被弱化,即在高速時(shí)其表現(xiàn)出欠驅(qū)動(dòng)特性。因此,研究欠驅(qū)動(dòng)條件約束下的UUV在水平面的路徑跟隨控制有重要意義[4,9]。

    本文所研究的UUV 長(zhǎng)×寬×高為6.77 m×1.98 m×1.72 m,如圖1 所示。UUV 重心與船體坐標(biāo)系的原點(diǎn)一致。該UUV 在oxy 平面上沿縱軸對(duì)稱(chēng)地裝有2個(gè)主推進(jìn)器;利用艉部十字舵不同控制面進(jìn)行機(jī)動(dòng)航行,水平面機(jī)動(dòng)時(shí)方向舵舵角的轉(zhuǎn)動(dòng)約束為±30°。設(shè)定UUV 的最大對(duì)流速度是6 節(jié)(約3.08 m/s)。

    圖1 UUV 坐標(biāo)系Fig.1 Earth-fixed and body-fixed coordinate systems

    結(jié)合UUV 載體的結(jié)構(gòu)特性,描述其水平面路徑跟隨控制的簡(jiǎn)化的運(yùn)動(dòng)方程可如式(1)~(4)所示。

    載體動(dòng)力學(xué):

    載體運(yùn)動(dòng)學(xué):

    2 水平面LOS 航行制導(dǎo)律

    在載體縱向運(yùn)動(dòng)航速可控的情況下,欲實(shí)現(xiàn)UUV在水平面的機(jī)動(dòng)即為了解決幾何學(xué)問(wèn)題,下面將從LOS 航向制導(dǎo)律和自適應(yīng)神經(jīng)網(wǎng)絡(luò)路徑跟隨控制器2方面對(duì)其進(jìn)行設(shè)計(jì),最終得到帶約束的自適應(yīng)神經(jīng)網(wǎng)絡(luò)水平面路徑跟隨控制器。

    2.1 直線段路徑的LOS 航向制導(dǎo)律

    假設(shè)在水平面UUV 的航行速度可表示為

    此時(shí),與航行機(jī)動(dòng)相關(guān)的航跡角為

    已知約束為平面內(nèi)的預(yù)規(guī)劃的路徑,要使UUV跟隨到指定位置,則需要研究與路徑相關(guān)的航行制導(dǎo)律以及控制載體的運(yùn)動(dòng)速度。若規(guī)劃的路徑可通過(guò)一系 列 點(diǎn) 連 接 表 示pk,k=0,1, …,n , 其 中pk=[ξk,ηk]T∈R2,則取跟隨路徑上2 點(diǎn)pk和pk+1,并以pk為原點(diǎn),沿路徑方向建立路徑坐標(biāo)系{p},如圖2 所示。

    圖2 LOS 航向制導(dǎo)律原理圖Fig.2 Principle of LOS guidance for straight-lines in horizontal plane

    由圖2 可知:沿路徑方向建立的坐標(biāo)系{p}的x軸與大地坐標(biāo)系北向的夾角為

    此時(shí),可在路徑坐標(biāo)系{p}中描述運(yùn)動(dòng)載體的坐標(biāo):

    2.2 圓弧段路徑的LOS 航向制導(dǎo)律

    參考直線段路徑的LOS 航向制導(dǎo)律的思想設(shè)計(jì)

    根據(jù)圖3 可知:圓弧段路徑的LOS 航向制導(dǎo)律為:

    該制導(dǎo)律在形式上與直線段路徑的LOS 航向制導(dǎo)律完全一致。

    圖3 圓弧段路徑的LOS 航向制導(dǎo)律的原理圖Fig.3 Principle of LOS guidance for circular arcs in horizontal plane

    3 自適應(yīng)神經(jīng)網(wǎng)絡(luò)路徑跟隨控制器設(shè)計(jì)

    為方便控制器設(shè)計(jì),定義

    由UUV 水平面機(jī)動(dòng)模型知,欲控制UUV 航向需引入航向運(yùn)動(dòng)差分方程,將水平面路徑跟隨控制的運(yùn)動(dòng)方程改寫(xiě)為

    定義航向跟蹤誤差:

    取K ∈R,使得式(9)描述的系統(tǒng)穩(wěn)定。

    那么,航向跟蹤誤差的微分形式可表示為

    其中:A= -K ,B=1??烧业絇>0,Q≥0 使得

    如果 f ( υ) 和b 均已知,并且系統(tǒng)干擾d=0,那么利用極點(diǎn)配置方法可得如下的線性運(yùn)動(dòng)控制器:

    結(jié)合航向控制運(yùn)動(dòng)方程,可知

    由于K 保證式(9)描述的系統(tǒng)穩(wěn)定,因此,

    實(shí)際中, f ( υ) 和b 有可能未知,同時(shí)系統(tǒng)干擾d也存在。為了解決非線性不確定項(xiàng)以及未知擾動(dòng)給控制器設(shè)計(jì)帶來(lái)的困難, 設(shè)計(jì)一個(gè)徑向基神經(jīng)網(wǎng)絡(luò)(RBF NN) 估計(jì) f ( υ) 及d:

    由于ε 是RBF NN 的估計(jì)誤差,故

    同時(shí),定義該RBF NN 的權(quán)值誤差為

    結(jié)合公式(12)知:

    在式(19)中,為補(bǔ)償神經(jīng)網(wǎng)絡(luò)的估計(jì)誤差及系統(tǒng)干擾項(xiàng)引起的控制偏差,設(shè)計(jì)1 個(gè)虛擬控制量

    因此,式(19)可改寫(xiě)為

    引入自適應(yīng)律:

    為保證設(shè)計(jì)的自適應(yīng)神經(jīng)網(wǎng)絡(luò)控制器給出的控制量可使與路徑跟隨相關(guān)的狀態(tài)均有界,選擇1 個(gè)Lyapunov 函數(shù):

    式(23)的微分形式為

    將式(21)代入式(24)得

    因?yàn)?/p>

    其中:ρ >0。

    基于Lyapunov 穩(wěn)定性理論可知:該路徑跟隨控制系統(tǒng)的狀態(tài)都是一致最終有界的,通過(guò)合適的選擇控制參數(shù)可以使得航向控制誤差漸近收斂至零。

    最終可得自適應(yīng)神經(jīng)網(wǎng)絡(luò)路徑跟隨控制器的具體表達(dá)式:

    4 仿真結(jié)果與分析

    在相同的初始條件,相同的控制目標(biāo)下進(jìn)行水平面路徑跟隨試驗(yàn),2 種路徑方案包括:直線段跟隨路徑和圓弧連接的直線段跟隨路徑。結(jié)合直線段LOS 航向制導(dǎo)律和圓弧路徑的圓弧LOS 航向制導(dǎo)律,利用自適應(yīng)神經(jīng)網(wǎng)絡(luò)控制器進(jìn)行水平面路徑跟隨控制對(duì)比試驗(yàn),初始值設(shè)定如下。

    表1 UUV 水平面無(wú)因次水動(dòng)力系數(shù)Table 1 UUV’s non-dimensional hydrodynamic coefficients

    潛器初始位置:(ξ ,η ,ζ ) = (0, -35,0.5) m ;

    初始姿態(tài)角:(φ ,θ ,ψ ) =(0°,0°,0°);

    指令速度:ud=1.8 m/s。

    另外,自適應(yīng)神經(jīng)網(wǎng)絡(luò)控制器的增益及控制參數(shù)設(shè)定為K=2.8,λ=0.1,γ1=2和γ2=1.5,RBF 神經(jīng)網(wǎng)絡(luò)的中心值為2,寬度為5,網(wǎng)絡(luò)初始權(quán)值取為0。

    環(huán)境中存在時(shí)不變海流,海流方向朝西,速度約為0.25 m/s。為盡可能接近實(shí)際情況,在直線段跟隨路徑中路徑為4 個(gè)直線段順序連接,其中包括3 種轉(zhuǎn)角:直角,鈍角和銳角,分別為90°,135°和45°。2種路徑方案下,仿真結(jié)果如圖4~8 所示。

    圖4 ANNC 下直線段跟隨路徑曲線Fig.4 Path following of straight-lines with adaptive NN controller

    圖4 和圖5 所示為2 種路徑下UUV 在北東坐標(biāo)系下的位置變化曲線,在速度為0.25 m/s,方向?yàn)橛蓶|向西的海流干擾下,顯然兩者的最大不同在航行段的交界處。在相同的初始條件,相同的控制目標(biāo)下進(jìn)行水平面路徑跟隨試驗(yàn),由圓弧連接的直線段路徑跟隨方案,當(dāng)航行段通過(guò)圓弧連接后,相比由直線段連接的直線段路徑跟隨方案,在轉(zhuǎn)角處的位置跟隨誤差大大減小。

    圖5 ANNC 下圓弧連接的直線段跟隨路徑曲線Fig.5 Path following of straight-lines and arcs with adaptive NN controller

    圖6 圓弧連接的直線段跟隨路徑時(shí)的速度曲線Fig.6 Speed response of straight-lines and arcs with adaptive NN controller

    圖7 2 種路徑方案下航行制導(dǎo)律的變化曲線Fig.7 Response of heading guidance under the two path profiles

    圖8 2 種路徑方案下橫向誤差曲線Fig.8 Cross-track errors of the two path profiles

    圖6(a)和(b)分別為圓弧連接的直線段跟隨路徑時(shí)的縱向運(yùn)動(dòng)速度曲線和橫向運(yùn)動(dòng)速度曲線。由圖6(a)可知:載體縱向運(yùn)動(dòng)航速可控,且航行速度在直線段航行過(guò)程中一直穩(wěn)定地保持在設(shè)定航速,僅在轉(zhuǎn)彎處有小幅度變化,這可消除減速轉(zhuǎn)彎對(duì)航行深度帶來(lái)的影響。根據(jù)圖7 可知:當(dāng)通過(guò)圓弧連接相鄰的航行段后,跟隨路徑的航向隨位置連續(xù)變化不再引起跳躍變化,則航向制導(dǎo)律得到的航向指令不會(huì)產(chǎn)生大幅度跳躍,此時(shí)在轉(zhuǎn)角處的位置誤差自然可減小。

    5 結(jié)論

    1) 針對(duì)UUV 水動(dòng)力參數(shù)變化和不確定項(xiàng)干擾下的精確水平面路徑跟隨問(wèn)題,基于李雅普諾夫穩(wěn)定性理論設(shè)計(jì)了自適應(yīng)神經(jīng)網(wǎng)絡(luò)控制器,可保證水平面路徑跟隨控制系統(tǒng)中所有的誤差狀態(tài)漸近穩(wěn)定。

    2) 在每個(gè)拐點(diǎn)處設(shè)計(jì)合適半徑的弧線連接兩個(gè)路徑部分,使得水平面機(jī)動(dòng)時(shí)拐角處的跟隨誤差限制在設(shè)計(jì)范圍內(nèi)。

    3) 仿真實(shí)驗(yàn)表明該方法是有效的,與單獨(dú)直線段路徑跟隨控制相比,改進(jìn)規(guī)則下自適應(yīng)神經(jīng)網(wǎng)絡(luò)的路徑跟隨控制器可使得拐點(diǎn)處的跟蹤誤差較大幅度減小,具有較好的工程使用價(jià)值。

    [1] 賈鶴鳴, 張利軍, 邊信黔, 等. 基于非線性迭代滑模的欠驅(qū)動(dòng)UUV 三維航跡跟蹤控制[J]. 自動(dòng)化學(xué)報(bào), 2012, 38(2):308-314.JIA Heming, ZHANG Lijun, BIAN Xinqian, et al. Three dimensional path following control for an underactuated UUV based on nonlinear iterative sliding mode[J]. Acta Automatica Sinica, 2012, 38(2): 308-314.

    [2] Wang B, Su Y M, Wan L, et al. Adaptive PID control system for an autonomous underwater vehicle[J]. High Technology Letters,2011, 17(1): 7-12.

    [3] Li J H, Lee P M. Path tracking in dive plane for a class of torpedotype underactuated AUVs[C]// Proceedings of the 7th Asian Control Conference. Piscataway, NJ, USA: IEEE, 2009, 8:360-365.

    [4] 張利軍, 賈鶴鳴, 邊信黔, 等. 基于L2干擾抑制的水下機(jī)器人三維航跡跟蹤控制[J]. 控制理論與應(yīng)用, 2011, 28(5): 645-651.ZHANG Lijun, JIA Heming, BIAN Xinqian, et al.Three-dimensional path tracking control for an autonomous underwater vehicle based on L-two disturbance attenuation method[J]. Control Theory & Applications, 2011, 28(5):645-651.

    [5] 唐旭東, 龐永杰, 李曄, 等. 基于混沌過(guò)程神經(jīng)元的水下機(jī)器人運(yùn)動(dòng)控制方法[J]. 控制與決策, 2010, 25(2): 213-217.TANG Xudong, PANG Yongjie, LI Ye, et al. Chaotic process neuron control for AUVs[J]. Control and Decision, 2010, 25(2):213-217.

    [6] 邊信黔, 程相勤, 賈鶴鳴, 等. 基于迭代滑模增量反饋的欠驅(qū)動(dòng) AUV 的地形跟蹤控制[J]. 控制與決策, 2011, 26(2):289-292.BIAN Xinqian, CHENG Xiangqin, JIA Heming, et al. A bottom-following controller for underactuated AUV based on iterative sliding and increment feedback[J]. Control and Decision,2011, 26(2): 289-292.

    [7] 畢鳳陽(yáng), 張嘉鐘, 魏英杰, 等. 模型不確定時(shí)滯欠驅(qū)動(dòng)AUV的模糊變結(jié)構(gòu)控制[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2010, 42(3):358-363.BI Fengyang, ZHANG Jiazhong, WEI Yingjie, et al. A fuzzy sliding mode controller for underactuated autonomous underwater vehicles with uncertain time delay[J]. Journal of Harbin Institute of Technology, 2010, 42(3): 358-363.

    [8] 施淑偉, 嚴(yán)衛(wèi)生, 高劍, 等. 常值海流作用下的AUV 水平面路徑跟蹤控制[J]. 兵工學(xué)報(bào), 2010, 31(3): 375-379.SHI Shuwei, YAN Weisheng, GAO Jian, et al. Path following control of an AUV in the horizontal plane with constant ocean currents[J]. Acta Armamentar II, 2010, 31(3): 375-379.

    [9] Fossen T I. Handbook of marine craft hydrodynamics and motion control[M]. Norway: A John Wiley & Sons, Ltd., 2011:243-246.

    [10] Oh S R, SUN Jing. Path following of underactuated marine surface vessels using line-of-sight based model predictive control[J]. Ocean Engineering, 2010, 37: 289-295.

    [11] B?rhaug E, Pavlov A, Pettersen K Y. Integral LOS control for path following of underactuated marine surface vessels in the presence of constant ocean currents[C]// Proceedings of the 47th IEEE Conference on Decision and Control. Cancun, Mexico,2008: 4984-4991.

    [12] BIAN Xinqian, ZHOU Jiajia, JIA Heming, et al. Adaptive neural network control system of bottom following for an underactuated AUV[C]// OCEANS'10 MTS/IEEE. Seattle, 2010: 1-6.

    [13] 賈鶴鳴, 張利軍, 齊雪, 等. 基于神經(jīng)網(wǎng)絡(luò)的水下機(jī)器人三維航跡跟蹤控制[J]. 控制理論與應(yīng)用, 2012, 29(7): 877-883.JIA Heming, ZHANG Lijun, QI Xue, et al. Three-dimensional path tracking control for autonomous underwater vehicle based on neural network[J]. Control Theory & Applications, 2012,29(7): 877-883.

    [14] Antonelli G. On the use of adaptive/integral actions for six-degrees-of-freedom control of autonomous underwater vehicles[J]. Ocean Engineering, 2007, 32(2): 300-312.

    [15] 于金霞, 蔡自興, 鄒小兵, 等. 基于神經(jīng)網(wǎng)絡(luò)辨識(shí)的移動(dòng)機(jī)器人航向誤差校準(zhǔn)方法[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2005,36(5): 745-750.YU Jinxia, CAI Zixing, ZOU Xiaobing, et al. Calibration method for heading error of mobile robot based on neural networks identification[J]. Journal of Central South University(Science and Technology), 2005, 36(5): 745-750.

    猜你喜歡
    水平面航向制導(dǎo)
    知坐標(biāo),明航向
    考慮幾何限制的航向道模式設(shè)計(jì)
    坡角多大,圓柱體在水平面滾得最遠(yuǎn)
    基于干擾觀測(cè)器的船舶系統(tǒng)航向Backstepping 控制
    電子制作(2017年24期)2017-02-02 07:14:16
    基于MPSC和CPN制導(dǎo)方法的協(xié)同制導(dǎo)律
    基于在線軌跡迭代的自適應(yīng)再入制導(dǎo)
    水平面內(nèi)勻速圓周運(yùn)動(dòng)的四種模型解讀
    帶有攻擊角約束的無(wú)抖振滑模制導(dǎo)律設(shè)計(jì)
    使命:引領(lǐng)航向與保持穩(wěn)定
    法大研究生(2015年2期)2015-02-27 10:13:55
    水平面上圓周運(yùn)動(dòng)中臨界問(wèn)題的分析和解題策略
    亚洲成人久久爱视频| av.在线天堂| 亚洲真实伦在线观看| 亚洲综合色惰| 日本爱情动作片www.在线观看 | 黄色丝袜av网址大全| 成人高潮视频无遮挡免费网站| 日韩精品有码人妻一区| 国产精品国产三级国产av玫瑰| 日日干狠狠操夜夜爽| 国产亚洲欧美98| 中文字幕高清在线视频| 人人妻,人人澡人人爽秒播| 国产欧美日韩精品亚洲av| 岛国在线免费视频观看| 免费黄网站久久成人精品| 亚洲,欧美,日韩| 男女啪啪激烈高潮av片| 久久久精品欧美日韩精品| 午夜免费激情av| 又爽又黄无遮挡网站| 看黄色毛片网站| 免费在线观看影片大全网站| 国产成人福利小说| 1024手机看黄色片| 美女被艹到高潮喷水动态| 永久网站在线| 久久精品综合一区二区三区| x7x7x7水蜜桃| 亚洲在线观看片| 亚洲精华国产精华液的使用体验 | 国产亚洲精品久久久com| 美女 人体艺术 gogo| 亚洲欧美精品综合久久99| 亚洲在线自拍视频| 看片在线看免费视频| a级一级毛片免费在线观看| 中文字幕高清在线视频| a级一级毛片免费在线观看| 尾随美女入室| 小说图片视频综合网站| 日本免费一区二区三区高清不卡| 国产一区二区三区av在线 | 啦啦啦韩国在线观看视频| 日韩国内少妇激情av| 夜夜看夜夜爽夜夜摸| 99国产极品粉嫩在线观看| 热99在线观看视频| 国产一区二区在线av高清观看| 日日干狠狠操夜夜爽| 日韩欧美精品免费久久| 成年人黄色毛片网站| 成人av在线播放网站| 最新在线观看一区二区三区| 色视频www国产| 欧美高清性xxxxhd video| 99久久成人亚洲精品观看| a级一级毛片免费在线观看| 人妻制服诱惑在线中文字幕| 国产一区二区三区视频了| 国产男靠女视频免费网站| 黄色女人牲交| 欧美一区二区国产精品久久精品| 中亚洲国语对白在线视频| 国产真实乱freesex| 亚洲欧美精品综合久久99| 国产欧美日韩精品一区二区| 网址你懂的国产日韩在线| 神马国产精品三级电影在线观看| 免费av毛片视频| 91久久精品国产一区二区三区| 观看美女的网站| 日韩欧美精品免费久久| 国产真实乱freesex| 色视频www国产| 国产精品永久免费网站| 亚洲经典国产精华液单| 亚洲国产欧洲综合997久久,| 免费搜索国产男女视频| 好男人在线观看高清免费视频| 国产亚洲欧美98| 又黄又爽又刺激的免费视频.| 成人永久免费在线观看视频| 国产精品福利在线免费观看| 国产一区二区在线av高清观看| 久久久久九九精品影院| 欧美丝袜亚洲另类 | 精品不卡国产一区二区三区| 在线观看av片永久免费下载| 国产一区二区三区在线臀色熟女| 欧美色欧美亚洲另类二区| 免费av观看视频| 自拍偷自拍亚洲精品老妇| 免费大片18禁| 国产精品伦人一区二区| h日本视频在线播放| 黄色女人牲交| bbb黄色大片| 在线观看舔阴道视频| 亚洲乱码一区二区免费版| 天美传媒精品一区二区| 黄色丝袜av网址大全| av.在线天堂| 少妇的逼好多水| 亚洲成人久久爱视频| 欧美日本视频| 又粗又爽又猛毛片免费看| 高清日韩中文字幕在线| 两个人视频免费观看高清| 亚洲av第一区精品v没综合| 九九久久精品国产亚洲av麻豆| 欧美成人性av电影在线观看| 欧美中文日本在线观看视频| 亚洲成av人片在线播放无| 一区二区三区高清视频在线| 免费高清视频大片| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播| 免费人成在线观看视频色| 99视频精品全部免费 在线| 搡女人真爽免费视频火全软件 | 淫秽高清视频在线观看| 在线免费十八禁| 日本三级黄在线观看| netflix在线观看网站| 国产午夜福利久久久久久| 国产精品免费一区二区三区在线| 亚洲国产精品成人综合色| 亚洲欧美精品综合久久99| 欧美高清成人免费视频www| 亚洲精品一卡2卡三卡4卡5卡| 免费看a级黄色片| 亚洲一级一片aⅴ在线观看| 免费人成视频x8x8入口观看| 在线免费观看不下载黄p国产 | 亚洲综合色惰| 91久久精品国产一区二区成人| 性欧美人与动物交配| 99精品在免费线老司机午夜| 九九在线视频观看精品| 婷婷色综合大香蕉| 亚洲无线在线观看| 国产精品嫩草影院av在线观看 | 国产精品久久久久久久电影| 国产 一区精品| 美女xxoo啪啪120秒动态图| 91麻豆av在线| 成人av一区二区三区在线看| 69av精品久久久久久| 老女人水多毛片| 22中文网久久字幕| av视频在线观看入口| 乱码一卡2卡4卡精品| 中文字幕av成人在线电影| 免费大片18禁| 啦啦啦观看免费观看视频高清| 久久欧美精品欧美久久欧美| 嫁个100分男人电影在线观看| 午夜视频国产福利| 91麻豆精品激情在线观看国产| 亚洲成人久久爱视频| 久久精品国产亚洲av涩爱 | 婷婷色综合大香蕉| 一进一出抽搐动态| 女同久久另类99精品国产91| 制服丝袜大香蕉在线| 国产精品永久免费网站| 久久人人精品亚洲av| 嫩草影院入口| 在线播放无遮挡| 少妇丰满av| 给我免费播放毛片高清在线观看| 国产亚洲精品av在线| 欧美bdsm另类| 简卡轻食公司| 女同久久另类99精品国产91| 免费观看人在逋| 国产综合懂色| 国产一区二区在线观看日韩| 在线免费观看不下载黄p国产 | 在线免费十八禁| 两性午夜刺激爽爽歪歪视频在线观看| 免费高清视频大片| 精品日产1卡2卡| 成年女人看的毛片在线观看| 日日干狠狠操夜夜爽| 99热这里只有是精品50| 精品久久久久久,| a在线观看视频网站| 免费一级毛片在线播放高清视频| 婷婷六月久久综合丁香| 在线看三级毛片| 看十八女毛片水多多多| 成人特级av手机在线观看| 精品午夜福利视频在线观看一区| 在线观看av片永久免费下载| 中国美女看黄片| 精品国内亚洲2022精品成人| 日日啪夜夜撸| 亚洲最大成人手机在线| 国语自产精品视频在线第100页| 亚洲内射少妇av| 91精品国产九色| 免费观看在线日韩| 欧美zozozo另类| 联通29元200g的流量卡| 久久久久久久久大av| 51国产日韩欧美| 白带黄色成豆腐渣| 免费观看人在逋| 精品一区二区三区av网在线观看| 日本 av在线| 国产成人av教育| 看十八女毛片水多多多| 舔av片在线| 99久久成人亚洲精品观看| 久久欧美精品欧美久久欧美| 日韩欧美精品免费久久| 一个人看视频在线观看www免费| 97超视频在线观看视频| 亚洲不卡免费看| 男女啪啪激烈高潮av片| 在线观看午夜福利视频| 国产毛片a区久久久久| 亚洲内射少妇av| 蜜桃久久精品国产亚洲av| 成人毛片a级毛片在线播放| 两个人视频免费观看高清| 窝窝影院91人妻| 日本 av在线| 国内精品久久久久久久电影| 九九在线视频观看精品| 99久久无色码亚洲精品果冻| 久久精品国产亚洲av涩爱 | 丰满的人妻完整版| 欧美成人免费av一区二区三区| 老熟妇仑乱视频hdxx| 日本 欧美在线| 在线国产一区二区在线| 国产男人的电影天堂91| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利18| 久9热在线精品视频| 伦理电影大哥的女人| 男女边吃奶边做爰视频| 男女那种视频在线观看| 黄片wwwwww| bbb黄色大片| 亚洲人成网站在线播放欧美日韩| 少妇猛男粗大的猛烈进出视频 | 日韩一区二区视频免费看| 久久亚洲精品不卡| 久久草成人影院| 成人av一区二区三区在线看| 久久九九热精品免费| 色综合色国产| 午夜亚洲福利在线播放| 久久久久久久久中文| 99热只有精品国产| 国产av一区在线观看免费| 亚洲乱码一区二区免费版| 麻豆成人午夜福利视频| 熟女电影av网| 午夜福利在线观看吧| 国产午夜精品论理片| 88av欧美| 精品人妻一区二区三区麻豆 | 一区二区三区高清视频在线| 美女cb高潮喷水在线观看| 内射极品少妇av片p| 九色成人免费人妻av| 九九爱精品视频在线观看| av中文乱码字幕在线| 国产精品嫩草影院av在线观看 | 午夜亚洲福利在线播放| 韩国av一区二区三区四区| 精品午夜福利在线看| 校园春色视频在线观看| 3wmmmm亚洲av在线观看| 91精品国产九色| 男人舔奶头视频| 男人的好看免费观看在线视频| 亚洲最大成人中文| 最近视频中文字幕2019在线8| 精品久久久久久久久久免费视频| 久久99热这里只有精品18| 久久99热6这里只有精品| 麻豆久久精品国产亚洲av| 国模一区二区三区四区视频| 免费av不卡在线播放| 欧美一区二区亚洲| 尾随美女入室| 乱人视频在线观看| 久久久久久久久中文| 国产男人的电影天堂91| 久久久久久久亚洲中文字幕| 国产主播在线观看一区二区| 久久久久性生活片| 久久精品人妻少妇| 国产人妻一区二区三区在| 夜夜看夜夜爽夜夜摸| 欧美成人a在线观看| 久久久久久久久久成人| 少妇裸体淫交视频免费看高清| 日韩一区二区视频免费看| 久久中文看片网| 国产久久久一区二区三区| 少妇人妻一区二区三区视频| 久久久国产成人精品二区| 日本黄大片高清| videossex国产| 成人国产综合亚洲| 免费搜索国产男女视频| 亚洲色图av天堂| 超碰av人人做人人爽久久| 久久久久久久久中文| 床上黄色一级片| 最近最新中文字幕大全电影3| 在线观看一区二区三区| 国产精品久久电影中文字幕| 亚洲美女视频黄频| 欧美高清性xxxxhd video| 一级a爱片免费观看的视频| 国产日本99.免费观看| 亚洲 国产 在线| 1024手机看黄色片| 99久久久亚洲精品蜜臀av| 桃红色精品国产亚洲av| 动漫黄色视频在线观看| 真人一进一出gif抽搐免费| 亚洲成人免费电影在线观看| 午夜福利在线观看吧| 成年版毛片免费区| 啪啪无遮挡十八禁网站| av福利片在线观看| 日韩欧美精品免费久久| 中文亚洲av片在线观看爽| 欧美黑人巨大hd| 欧美精品啪啪一区二区三区| 色播亚洲综合网| 不卡视频在线观看欧美| 女人十人毛片免费观看3o分钟| 久久精品国产亚洲av涩爱 | 成年女人永久免费观看视频| 日日摸夜夜添夜夜添小说| 欧美黑人欧美精品刺激| 99热这里只有是精品50| 天堂√8在线中文| 淫妇啪啪啪对白视频| 97超视频在线观看视频| 国内精品久久久久精免费| 国产精品免费一区二区三区在线| 精品一区二区三区av网在线观看| 99在线人妻在线中文字幕| 日韩大尺度精品在线看网址| 久久精品影院6| 亚洲在线观看片| 精品午夜福利在线看| 男女那种视频在线观看| 哪里可以看免费的av片| 看黄色毛片网站| 亚洲无线在线观看| 色5月婷婷丁香| av在线天堂中文字幕| netflix在线观看网站| 亚洲狠狠婷婷综合久久图片| 中国美白少妇内射xxxbb| 国产精品自产拍在线观看55亚洲| 久久久久久久久久成人| 女人十人毛片免费观看3o分钟| 看免费成人av毛片| 亚洲五月天丁香| 在线观看66精品国产| 欧美xxxx黑人xx丫x性爽| 国语自产精品视频在线第100页| 国产精品久久久久久av不卡| 亚洲人成网站在线播放欧美日韩| 在现免费观看毛片| 日韩人妻高清精品专区| 桃色一区二区三区在线观看| 日本色播在线视频| 亚洲avbb在线观看| 久久久久久伊人网av| 黄色丝袜av网址大全| 久久精品国产清高在天天线| 成人美女网站在线观看视频| 国产一区二区亚洲精品在线观看| 韩国av一区二区三区四区| 欧美性猛交╳xxx乱大交人| 18+在线观看网站| 极品教师在线免费播放| 看免费成人av毛片| 一区二区三区免费毛片| 精品一区二区三区视频在线观看免费| 在线免费十八禁| 男女啪啪激烈高潮av片| 欧美丝袜亚洲另类 | 欧美日韩亚洲国产一区二区在线观看| 久久婷婷人人爽人人干人人爱| 亚洲第一电影网av| 88av欧美| 小蜜桃在线观看免费完整版高清| 我要搜黄色片| 伦理电影大哥的女人| 日本 欧美在线| 精品不卡国产一区二区三区| 日韩亚洲欧美综合| 波多野结衣高清无吗| 国产亚洲欧美98| 一级a爱片免费观看的视频| 性欧美人与动物交配| 亚洲成人中文字幕在线播放| 欧美国产日韩亚洲一区| 免费不卡的大黄色大毛片视频在线观看 | 日本 欧美在线| 女同久久另类99精品国产91| 最近在线观看免费完整版| 国内少妇人妻偷人精品xxx网站| 热99re8久久精品国产| av天堂在线播放| 久久草成人影院| 国产主播在线观看一区二区| 精品一区二区三区av网在线观看| 淫秽高清视频在线观看| 欧美日韩瑟瑟在线播放| 免费搜索国产男女视频| 国产一级毛片七仙女欲春2| 午夜精品久久久久久毛片777| 亚洲精品亚洲一区二区| 又爽又黄无遮挡网站| 成人国产麻豆网| 精品人妻一区二区三区麻豆 | 亚洲精品亚洲一区二区| 18禁在线播放成人免费| 97热精品久久久久久| 中国美女看黄片| 国产亚洲精品久久久久久毛片| 国产成人影院久久av| 看免费成人av毛片| 一夜夜www| 色5月婷婷丁香| 国产成人福利小说| 中文字幕久久专区| 久久精品影院6| 精品久久久噜噜| 亚州av有码| 国产高清三级在线| 久久久久久久久大av| 亚洲欧美精品综合久久99| 亚洲专区中文字幕在线| 少妇猛男粗大的猛烈进出视频 | 国产成人影院久久av| 亚洲专区国产一区二区| 国产高清有码在线观看视频| 欧美高清成人免费视频www| 亚洲国产高清在线一区二区三| 嫩草影院入口| 特级一级黄色大片| 舔av片在线| 联通29元200g的流量卡| www.色视频.com| 欧洲精品卡2卡3卡4卡5卡区| 日韩中字成人| 成年免费大片在线观看| 欧美+亚洲+日韩+国产| 91久久精品国产一区二区成人| 欧美日本视频| 国产蜜桃级精品一区二区三区| 午夜久久久久精精品| 久久中文看片网| 国产视频内射| 12—13女人毛片做爰片一| 成年人黄色毛片网站| 欧美区成人在线视频| 日本免费a在线| 天堂av国产一区二区熟女人妻| 性插视频无遮挡在线免费观看| 波多野结衣巨乳人妻| 99在线人妻在线中文字幕| 欧美一级a爱片免费观看看| 亚洲av.av天堂| 高清在线国产一区| 国产高清激情床上av| 精品久久久久久成人av| 久久精品国产自在天天线| 日本黄色片子视频| 热99在线观看视频| 国产精品电影一区二区三区| 天美传媒精品一区二区| av.在线天堂| 国产日本99.免费观看| 国产高清激情床上av| 精品国内亚洲2022精品成人| 91在线精品国自产拍蜜月| 日本黄色视频三级网站网址| aaaaa片日本免费| 日韩精品有码人妻一区| 99精品在免费线老司机午夜| 91在线精品国自产拍蜜月| 免费观看在线日韩| 国产日本99.免费观看| 国产伦精品一区二区三区视频9| 久久久精品大字幕| 在线观看美女被高潮喷水网站| 在线观看免费视频日本深夜| 亚洲欧美激情综合另类| 他把我摸到了高潮在线观看| 国产在线男女| 中亚洲国语对白在线视频| 久久热精品热| 成人亚洲精品av一区二区| 国产精品人妻久久久影院| 欧美日韩国产亚洲二区| 美女cb高潮喷水在线观看| 亚洲真实伦在线观看| 久久久久精品国产欧美久久久| 色综合站精品国产| 成年版毛片免费区| 乱系列少妇在线播放| 欧美bdsm另类| 国产av麻豆久久久久久久| 欧美zozozo另类| 悠悠久久av| 极品教师在线视频| 欧美精品国产亚洲| 久久精品国产亚洲av涩爱 | 久久6这里有精品| 国产一区二区三区视频了| 成人二区视频| 国产亚洲91精品色在线| 成人高潮视频无遮挡免费网站| 日韩亚洲欧美综合| 精品久久久久久久久亚洲 | 日本a在线网址| 动漫黄色视频在线观看| 久久草成人影院| 免费高清视频大片| 国产精品国产高清国产av| 亚洲国产色片| 搡女人真爽免费视频火全软件 | 亚洲经典国产精华液单| 无人区码免费观看不卡| 他把我摸到了高潮在线观看| 韩国av一区二区三区四区| 一本久久中文字幕| 精品人妻视频免费看| 九色国产91popny在线| 18禁在线播放成人免费| 999久久久精品免费观看国产| 久久草成人影院| 国产国拍精品亚洲av在线观看| 国产毛片a区久久久久| 能在线免费观看的黄片| 99久久精品热视频| 国产精品不卡视频一区二区| 春色校园在线视频观看| 亚洲精品日韩av片在线观看| 99精品在免费线老司机午夜| 国产精品久久久久久av不卡| 日日啪夜夜撸| 又黄又爽又刺激的免费视频.| 国产黄色小视频在线观看| 我的老师免费观看完整版| 99久国产av精品| 亚洲精品粉嫩美女一区| 欧美最新免费一区二区三区| 一级黄片播放器| 精品乱码久久久久久99久播| 日本熟妇午夜| 国产精品美女特级片免费视频播放器| 久久99热这里只有精品18| 韩国av一区二区三区四区| 成年人黄色毛片网站| 午夜福利欧美成人| 中文字幕av在线有码专区| 欧美成人一区二区免费高清观看| 午夜福利在线在线| 国产又黄又爽又无遮挡在线| 亚洲av一区综合| 99久久久亚洲精品蜜臀av| 亚洲不卡免费看| 国产精品精品国产色婷婷| 国产一级毛片七仙女欲春2| 校园春色视频在线观看| a级毛片a级免费在线| 亚洲成人免费电影在线观看| 极品教师在线视频| 丰满乱子伦码专区| 国产视频一区二区在线看| 免费看av在线观看网站| 国产成人福利小说| 小说图片视频综合网站| 成年女人毛片免费观看观看9| 国产av不卡久久| 一个人看的www免费观看视频| 亚洲成a人片在线一区二区| 亚洲av成人av| 国产白丝娇喘喷水9色精品| 国产中年淑女户外野战色| 亚洲精品色激情综合| 又黄又爽又免费观看的视频| 一级黄片播放器| 精品午夜福利在线看| 亚洲经典国产精华液单| 麻豆av噜噜一区二区三区| 床上黄色一级片| 国产视频一区二区在线看| 如何舔出高潮| 欧美一区二区国产精品久久精品| 亚洲中文字幕日韩| 99久国产av精品| 国产一区二区三区在线臀色熟女| 国内精品美女久久久久久| 直男gayav资源| 69人妻影院|