999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

雙軸旋轉(zhuǎn)式激光捷聯(lián)慣導(dǎo)系統(tǒng)的轉(zhuǎn)位方案研究

2016-03-16 06:36:44陳家斌韓勇強(qiáng)宋春雷王明杰
導(dǎo)航定位與授時(shí) 2016年4期
關(guān)鍵詞:系統(tǒng)

秦 沖,陳家斌,韓勇強(qiáng),宋春雷,王明杰

(北京理工大學(xué) 自動(dòng)化學(xué)院,北京 100081)

雙軸旋轉(zhuǎn)式激光捷聯(lián)慣導(dǎo)系統(tǒng)的轉(zhuǎn)位方案研究

秦 沖,陳家斌,韓勇強(qiáng),宋春雷,王明杰

(北京理工大學(xué) 自動(dòng)化學(xué)院,北京 100081)

慣性導(dǎo)航系統(tǒng)的誤差隨時(shí)間累積,旋轉(zhuǎn)調(diào)制技術(shù)可以有效地提高慣導(dǎo)系統(tǒng)的長(zhǎng)航時(shí)精度,旋轉(zhuǎn)調(diào)制方案是決定旋轉(zhuǎn)式捷聯(lián)慣導(dǎo)系統(tǒng)導(dǎo)航精度的一個(gè)重要因素。針對(duì)雙軸旋轉(zhuǎn)慣導(dǎo)系統(tǒng),相較于16次序轉(zhuǎn)位方案,提出了一種新的32次序雙軸旋轉(zhuǎn)調(diào)制方案。根據(jù)捷聯(lián)慣導(dǎo)系統(tǒng)的誤差方程,推導(dǎo)出旋轉(zhuǎn)捷聯(lián)慣導(dǎo)的誤差方程,分析了誤差補(bǔ)償?shù)臋C(jī)理,研究了慣性器件常值偏置誤差、標(biāo)度因數(shù)誤差和安裝角誤差的傳播特性。仿真結(jié)果表明,32次序雙軸旋轉(zhuǎn)調(diào)制方案相對(duì)于16次序轉(zhuǎn)位方案有明顯的優(yōu)勢(shì),可以有效地降低姿態(tài)角誤差和經(jīng)緯度誤差。

捷聯(lián)慣導(dǎo)系統(tǒng);激光陀螺;旋轉(zhuǎn)調(diào)制方案;誤差補(bǔ)償

0 引言

旋轉(zhuǎn)調(diào)制技術(shù)是一種誤差自補(bǔ)償技術(shù),利用IMU周期性轉(zhuǎn)動(dòng)完成對(duì)慣性器件誤差的調(diào)制,從而提高導(dǎo)航系統(tǒng)的精度。目前,國(guó)外雙軸旋轉(zhuǎn)式激光捷聯(lián)慣性導(dǎo)航系統(tǒng)已經(jīng)大量裝備于海軍艦船,Sperry公司的MK49型雙軸旋轉(zhuǎn)式激光陀螺捷聯(lián)慣導(dǎo)系統(tǒng)經(jīng)過(guò)海試,被選為北約的標(biāo)準(zhǔn)船用慣導(dǎo)系統(tǒng),裝備了大量的艦艇。在MK39/49基礎(chǔ)上,Sperry公司發(fā)展了AN/WSN-7系列雙軸旋轉(zhuǎn)式激光陀螺捷聯(lián)式慣導(dǎo)系統(tǒng),重調(diào)周期可達(dá)14d,成為美國(guó)海軍艦艇的新一代標(biāo)準(zhǔn)慣性導(dǎo)航系統(tǒng),更高精度的激光陀螺旋轉(zhuǎn)調(diào)制系統(tǒng)正在研制中[1-2]。國(guó)內(nèi),在2000年左右誤差自補(bǔ)償技術(shù)的研究才開(kāi)始被重視起來(lái),在雙軸旋轉(zhuǎn)慣導(dǎo)系統(tǒng)調(diào)制機(jī)理、高精度IMU設(shè)計(jì)、轉(zhuǎn)位機(jī)構(gòu)設(shè)計(jì)和轉(zhuǎn)位方案設(shè)計(jì)方面也取得了一定的成果。

目前,公開(kāi)資料上很難找到MK-49、AN/WSN-7A等國(guó)外雙軸旋轉(zhuǎn)系統(tǒng)所使用的旋轉(zhuǎn)方案[2]。轉(zhuǎn)位方案是影響雙軸旋轉(zhuǎn)慣導(dǎo)系統(tǒng)導(dǎo)航精度的重要因素,不同的轉(zhuǎn)位方案將會(huì)產(chǎn)生不同的調(diào)制效果。文獻(xiàn)[4]在靜電陀螺平臺(tái)翻滾方案的基礎(chǔ)上,提出了旋轉(zhuǎn)式系統(tǒng)的8次序雙軸轉(zhuǎn)位方案,并進(jìn)一步提出了16次序雙軸轉(zhuǎn)位方案。文獻(xiàn)[5]研究了轉(zhuǎn)位機(jī)構(gòu)的調(diào)制角速度、時(shí)間參數(shù)設(shè)置問(wèn)題,指出調(diào)制角速度應(yīng)設(shè)置在15(°)/s以下,不同的調(diào)制方案,調(diào)制角速度選擇不盡相同,須進(jìn)行比較分析。文獻(xiàn)[7]提出了一種改進(jìn)的16次序雙軸交替旋轉(zhuǎn)調(diào)制方法,進(jìn)一步抑制了速度誤差積累所引起的位置誤差。文獻(xiàn)[8]在比較16次序轉(zhuǎn)位方案和64次序轉(zhuǎn)位方案基礎(chǔ)上,提出了改進(jìn)的64次序雙軸轉(zhuǎn)位方案,有效地提高了系統(tǒng)的導(dǎo)航精度。文獻(xiàn)[14]在16次序轉(zhuǎn)位方案基礎(chǔ)上,將每次180°的旋轉(zhuǎn)改為兩次90°的旋轉(zhuǎn),在每個(gè)位置停留相同時(shí)間來(lái)補(bǔ)償二次諧波誤差,從而構(gòu)成了一個(gè)32次序的雙軸轉(zhuǎn)位方案。

本文綜合考慮以上研究?jī)?nèi)容,以雙軸旋轉(zhuǎn)捷聯(lián)慣性導(dǎo)航系統(tǒng)為研究對(duì)象,從旋轉(zhuǎn)調(diào)制型激光捷聯(lián)慣性導(dǎo)航系統(tǒng)的基本原理出發(fā),提出了一種新的32次序雙軸轉(zhuǎn)位方案,可以對(duì)慣性器件的常值誤差、標(biāo)度因數(shù)誤差和安裝角誤差進(jìn)行有效的補(bǔ)償,可有效地提高導(dǎo)航精度。

1 旋轉(zhuǎn)調(diào)制原理

旋轉(zhuǎn)調(diào)制型激光陀螺捷聯(lián)慣性導(dǎo)航系統(tǒng)是在激光陀螺捷聯(lián)慣性導(dǎo)航系統(tǒng)的外面加上轉(zhuǎn)動(dòng)機(jī)構(gòu)和測(cè)角裝置,導(dǎo)航解算采用捷聯(lián)慣性導(dǎo)航算法,計(jì)算出IMU的姿態(tài),并根據(jù)IMU相對(duì)于載體的轉(zhuǎn)動(dòng)角度,從而轉(zhuǎn)化得到載體的姿態(tài)信息[5]。旋轉(zhuǎn)調(diào)制原理框圖如圖1所示。

以傳統(tǒng)正交安裝的IMU為研究對(duì)象,IMU固聯(lián)在旋轉(zhuǎn)機(jī)構(gòu)上,旋轉(zhuǎn)機(jī)構(gòu)繞相互正交的坐標(biāo)軸旋轉(zhuǎn),構(gòu)成一個(gè)新的旋轉(zhuǎn)坐標(biāo)系。定義:s系為旋轉(zhuǎn)坐標(biāo)系;b系為載體坐標(biāo)系;n系為導(dǎo)航坐標(biāo)系;i系為慣性坐標(biāo)系。初始時(shí)刻旋轉(zhuǎn)坐標(biāo)系與載體坐標(biāo)系重合。

圖1 旋轉(zhuǎn)調(diào)制原理框圖Fig.1 Principles of rotation strap-down algorithm

慣性元件的輸出描述如下:

(1)

(2)

由式(1)、式(2)展開(kāi),略去二階小量,可推出慣性元件的輸出誤差方程:

(3)

δfs=(δKa+δMa)fs+▽

(4)

根據(jù)捷聯(lián)慣性導(dǎo)航系統(tǒng)的誤差方程,結(jié)合旋轉(zhuǎn)捷聯(lián)解算框圖,推出旋轉(zhuǎn)慣性導(dǎo)航系統(tǒng)的誤差傳播方程如下:

(5)

(6)

(7)

2 旋轉(zhuǎn)調(diào)制方案設(shè)計(jì)與分析

單軸旋轉(zhuǎn)調(diào)制系統(tǒng)只能消除與旋轉(zhuǎn)軸垂直方向上的慣性器件偏差,旋轉(zhuǎn)軸方向上的器件誤差依然會(huì)引起慣性導(dǎo)航系統(tǒng)的誤差[6],而雙軸旋轉(zhuǎn)調(diào)制系統(tǒng)可以完全調(diào)制3個(gè)方向上的常值陀螺漂移。在旋轉(zhuǎn)過(guò)程中,IMU在每個(gè)位置駐停時(shí)間越短,慣性器件誤差的調(diào)制效果越好,因此可采用繞水平軸和天向軸的連續(xù)轉(zhuǎn)位方案,實(shí)現(xiàn)3個(gè)方向上慣性器件偏差的調(diào)制。

2.1 旋轉(zhuǎn)調(diào)制方案的設(shè)計(jì)

目前,采用較多的是16次序的雙軸轉(zhuǎn)位方案,此方案中的轉(zhuǎn)位機(jī)構(gòu)轉(zhuǎn)角可達(dá)720°,如表1所示。

表1 16次序雙軸轉(zhuǎn)位方案

為了改善16次序轉(zhuǎn)位方案中轉(zhuǎn)角過(guò)大問(wèn)題,提高系統(tǒng)的導(dǎo)航精度,從而設(shè)計(jì)了一種32次序的轉(zhuǎn)位方案。一個(gè)合理的雙軸旋轉(zhuǎn)方案設(shè)計(jì)原則應(yīng)滿足:繞水平軸和天向軸周期性的交互旋轉(zhuǎn),繞每個(gè)軸的旋轉(zhuǎn)均具有正反性和對(duì)稱(chēng)性;在一個(gè)轉(zhuǎn)動(dòng)周期內(nèi),慣性元件引起的累計(jì)角度誤差和速度誤差應(yīng)為0。文獻(xiàn)[9]指出,在一個(gè)小的旋轉(zhuǎn)周期內(nèi),每一對(duì)同軸旋轉(zhuǎn)的轉(zhuǎn)動(dòng)方向須相同,每一對(duì)異軸旋轉(zhuǎn)的轉(zhuǎn)動(dòng)方向須相反。

本文根據(jù)上述規(guī)則,通過(guò)調(diào)整轉(zhuǎn)位次序,設(shè)計(jì)了32次序的轉(zhuǎn)位方案,每8個(gè)次序?yàn)橐粋€(gè)小周期,每16次序奇對(duì)稱(chēng),每32次序偶對(duì)稱(chēng),調(diào)整后的32次序方案轉(zhuǎn)角在-180°~+180°,轉(zhuǎn)位機(jī)構(gòu)水平軸的轉(zhuǎn)動(dòng)過(guò)程中可以省掉滑環(huán),從而提高系統(tǒng)的穩(wěn)定性。轉(zhuǎn)位次序如表2所示。

表2 32次序雙軸轉(zhuǎn)位方案

續(xù)表

2.2 慣性器件常值誤差調(diào)制

這里,只討論在一個(gè)小周期8次序內(nèi)的情況,設(shè)轉(zhuǎn)動(dòng)過(guò)程的時(shí)間為T(mén)s。

次序1、3、6、8繞z軸轉(zhuǎn)動(dòng)過(guò)程中,陀螺漂移在載體坐標(biāo)系x軸的投影經(jīng)過(guò)積分后結(jié)果為:

(8)

(9)

(10)

(11)

同理,陀螺漂移在載體坐標(biāo)系y軸的投影經(jīng)過(guò)積分后的結(jié)果為:

(12)

(13)

(14)

(15)

(16)

(17)

次序2、4、5、7繞y軸轉(zhuǎn)動(dòng)過(guò)程中,x、z軸陀螺漂移在載體坐標(biāo)系中呈現(xiàn)正反各一周的變化規(guī)律,整周期內(nèi)積分產(chǎn)生的常值偏置為0,即:

(18)

(19)

32次序的轉(zhuǎn)動(dòng)方案由8次序小周期組成,每16次序奇對(duì)稱(chēng),每32次序偶對(duì)稱(chēng),因此32次序轉(zhuǎn)動(dòng)方案中,陀螺漂移完全調(diào)制。同理,加速度計(jì)的常值偏置誤差也可以完全調(diào)制。

2.3 慣性器件標(biāo)度因數(shù)誤差調(diào)制

32次序轉(zhuǎn)動(dòng)方案主要為繞z軸和y軸進(jìn)行旋轉(zhuǎn),這里只分析繞z軸和繞y軸的轉(zhuǎn)動(dòng)情況。

設(shè)初始時(shí)刻導(dǎo)航坐標(biāo)系與載體坐標(biāo)系重合,繞z軸正向旋轉(zhuǎn),導(dǎo)航坐標(biāo)系下由于標(biāo)度因數(shù)誤差存在導(dǎo)致陀螺儀輸出誤差為

(20)

由于陀螺儀標(biāo)度因數(shù)誤差引起的姿態(tài)角誤差:

(21)

(22)

(23)

設(shè)初始時(shí)刻導(dǎo)航坐標(biāo)系與載體坐標(biāo)系重合,繞y軸正向轉(zhuǎn)動(dòng),由于陀螺標(biāo)度因數(shù)誤差引起的姿態(tài)角誤差為

(24)

由式(21)~式(24)可以得到,標(biāo)度因數(shù)誤差引起的姿態(tài)角誤差在導(dǎo)航坐標(biāo)系下水平東向的誤差為0。在水平北向和天向方向上均產(chǎn)生了常值偏差,姿態(tài)角誤差大小與旋轉(zhuǎn)角速度和標(biāo)度因數(shù)大小有關(guān)。

2.4 慣性器件安裝誤差調(diào)制

設(shè)初始時(shí)刻導(dǎo)航坐標(biāo)系與載體坐標(biāo)系重合,繞z軸正向旋轉(zhuǎn),導(dǎo)航坐標(biāo)系下由于安裝誤差導(dǎo)致的陀螺儀輸出誤差為

(25)

(26)

(27)

由式(26)、式(27)得到,陀螺安裝誤差經(jīng)過(guò)調(diào)制不能消除姿態(tài)角誤差的積累,出現(xiàn)了安裝誤差與地球自轉(zhuǎn)角速度分量耦合。

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

對(duì)上述16次序轉(zhuǎn)位方案和32次序轉(zhuǎn)位方案進(jìn)行仿真驗(yàn)證。仿真條件設(shè)置如下:陀螺常值漂移0.01(°)/h,陀螺的刻度系數(shù)誤差為1×10-5,陀螺安裝誤差為10″,加速度計(jì)常值偏置為50μg,加速度計(jì)刻度系數(shù)誤差為5×10-5,加速度計(jì)安裝誤差角為10″,初始位置誤差、速度誤差、初始對(duì)準(zhǔn)誤差均為0,旋轉(zhuǎn)角速度設(shè)置為10(°)/s,忽略高度誤差。

在32次序和16次序完整周期內(nèi),只考慮陀螺常值漂移情況下的實(shí)驗(yàn)仿真結(jié)果如圖2和圖3所示。

圖2 16次序和32次序轉(zhuǎn)位方案導(dǎo)航姿態(tài)角誤差Fig.2 Navigation output attitude angle errors in different schemes

圖3 16次序和32次序轉(zhuǎn)位方案導(dǎo)航經(jīng)緯度誤差Fig.3 Longitude and latitude errors of navigation output position in different schemes

在32次序和16次序完整周期內(nèi),只考慮陀螺標(biāo)度因數(shù)誤差情況下的實(shí)驗(yàn)仿真結(jié)果如圖4和圖5所示。

圖4 16次序和32次序轉(zhuǎn)位方案導(dǎo)航姿態(tài)角誤差Fig.4 Navigation output attitude angle errors in different schemes

圖5 16次序和32次序轉(zhuǎn)位方案導(dǎo)航經(jīng)緯度誤差Fig.5 Longitude and latitude errors of navigation output position in different schemes

考慮陀螺儀和加速度計(jì)的所有誤差,并考慮長(zhǎng)航時(shí)情況下的導(dǎo)航誤差情況,設(shè)置了80h的仿真時(shí)間,仿真結(jié)果圖如圖6和圖7所示。

圖6 16次序和32次序轉(zhuǎn)位方案導(dǎo)航姿態(tài)角誤差Fig.6 Navigation output attitude angle errors in different schemes

圖7 16次序和32次序轉(zhuǎn)位方案導(dǎo)航經(jīng)緯度誤差Fig.7 Longitude and latitude errors of navigation output position in different schemes

仿真結(jié)果中,圖2和圖3表明16次序轉(zhuǎn)位方案和32次序轉(zhuǎn)位方案對(duì)姿態(tài)角誤差調(diào)制區(qū)別較小,雖然32次序東向姿態(tài)角誤差峰值稍大,但緯度誤差明顯減小。圖4和圖5中顯示32次序轉(zhuǎn)位方案中東向姿態(tài)角誤差和北向姿態(tài)角誤差明顯的減小,經(jīng)緯度誤差也有明顯的減小。圖6和圖7表明,考慮所有誤差情況并在長(zhǎng)航時(shí)狀態(tài)下,32次序轉(zhuǎn)位方案相較于16次序轉(zhuǎn)位方案,姿態(tài)角誤差和經(jīng)緯度誤差都有了較明顯的減小,導(dǎo)航精度有所提高。

4 結(jié)論

本文以雙軸旋轉(zhuǎn)慣性導(dǎo)航系統(tǒng)為研究對(duì)象,為了提高雙軸旋轉(zhuǎn)調(diào)制系統(tǒng)的精度,針對(duì)16次序轉(zhuǎn)位方案,提出了32次序的調(diào)制方案。文章分析研究了旋轉(zhuǎn)過(guò)程中慣性器件常值誤差、標(biāo)度因數(shù)誤差、安裝誤差的傳播特性,分析表明了繞水平軸和天向軸旋轉(zhuǎn)可以有效地調(diào)制慣性器件常值偏差,但是標(biāo)度因數(shù)誤差和安裝誤差不能被完全調(diào)制。通過(guò)仿真驗(yàn)證,32次序轉(zhuǎn)位方案可以有效地提高慣性導(dǎo)航系統(tǒng)的導(dǎo)航精度,有一定的工程應(yīng)用價(jià)值。

[1] Hibbard R, Wylie B, Levison E.Sperry marine MK-49, The world’s best ring laser gyro ship’s inertial navigation system[C].JSDE Proceedings, 1996: 35-46.

[2] Tucker T, Levinson E.The AN/WSN-7B marine gyrocompass navigator[C].Proceedings of the 2000 National Technical Meeting of the Institute of Navigation.London, 2000: 348-357.

[3] 孫騫.雙軸旋轉(zhuǎn)式光纖陀螺捷聯(lián)慣導(dǎo)系統(tǒng)技術(shù)[D].哈爾濱工程大學(xué),2013.

[4] 袁保倫.四頻激光陀螺旋轉(zhuǎn)式慣導(dǎo)系統(tǒng)研究[D].長(zhǎng)沙:國(guó)防科學(xué)技術(shù)大學(xué),2007.

[5] 王庭軍.旋轉(zhuǎn)式捷聯(lián)慣導(dǎo)系統(tǒng)關(guān)鍵技術(shù)研究[D].哈爾濱工程大學(xué),2013.

[6] 毛玉良.激光陀螺捷聯(lián)慣導(dǎo)系統(tǒng)誤差辨識(shí)與修正技術(shù)研究[D].北京理工大學(xué),2014.

[7] 劉潔瑜,蔚國(guó)強(qiáng),楊建業(yè).雙軸連續(xù)旋轉(zhuǎn)激光捷聯(lián)慣導(dǎo)誤差高精度補(bǔ)償方法[J].系統(tǒng)工程與電子技術(shù),2015,37(1):148-154.

[8] 魏國(guó).二頻機(jī)抖激光陀螺雙軸旋轉(zhuǎn)慣性導(dǎo)航系統(tǒng)若干關(guān)鍵技術(shù)研究[D].長(zhǎng)沙:國(guó)防科學(xué)技術(shù)大學(xué),2013.

[9] 于飛,阮雙雙.光纖捷聯(lián)慣導(dǎo)系統(tǒng)的雙軸旋轉(zhuǎn)調(diào)制方案[J].哈爾濱工程大學(xué)報(bào),2014,35(12):1536-1542.

[10] 孫偉,初婧,李瑞豹,張媛.旋轉(zhuǎn)慣導(dǎo)系統(tǒng)中IMU轉(zhuǎn)停時(shí)間分配技術(shù)研究[J].壓電與聲光,2014,36(2):225-229+233.

[11] Levinson E.Accuracy enhancement techniques to the marine ring laser inertial navigation[J].Navigation, 1987, 34 (1): 64-87.

[12] 尹洪亮,楊功流,宋凝芳,王麗芬.旋轉(zhuǎn)激光陀螺慣導(dǎo)系統(tǒng)誤差傳播特性分析[J].北京航空航天大學(xué)學(xué)報(bào),2012,38(3):345-350.

[13] 陸志東,王曉斌.系統(tǒng)級(jí)雙軸旋轉(zhuǎn)調(diào)制捷聯(lián)慣導(dǎo)誤差分析及標(biāo)校[J].中國(guó)慣性技術(shù)學(xué)報(bào),2010,18(2):135-141.

[14] 張鵬.激光陀螺捷聯(lián)慣導(dǎo)系統(tǒng)旋轉(zhuǎn)調(diào)制技術(shù)研究[D].哈爾濱工業(yè)大學(xué),2010.

Research on Rotating Scheme for Dual-Axis Rotation Laser SINS

QIN Chong, CHEN Jia-bin, HAN Yong-qiang, SONG Chun-lei, WANG Ming-jie

(School of Automation,Beijing Institute of Technology,Beijing 100081,China)

The error of inertial navigation system has the characteristics of accumulation over time.However, error rotating compensation method can effectively improve the long-time navigation accuracy.And the rotating scheme of the turntable is an important factor to determine the precision of rotary SINS.Based on the dual-axis rotating SINS, compared with the sixteen sequence dual-axis indexing scheme, a thirty-two sequence indexing scheme is put forward to modulate the inertial instrument errors better.By deriving the error model of strap-down inertial navigation system and analyzing the navigation propagating function, the mechanism of error rotating compensation is described in detail.And the propagation characteristics of inertial sensor errors, scale factor errors and installation errors are also studied.The simulation results show that the thirty-two sequence dual-axis indexing scheme has obvious advantages when compared with the sixteen sequence dual-axis indexing scheme.And it can effectively reduce the amplitude of attitude angle errors and longitude and latitude errors.

Strap-down inertial navigation system; Ring laser gyroscope; Rotation modulation scheme; Error compensation

10.19306/j.cnki.2095-8110.2016.04.004

2016-04-08;

2016-05-05。

國(guó)家國(guó)防基金(9140A09050313BQ01127);國(guó)家自然科學(xué)基金(91120010)

秦沖(1991-),男,碩士,主要從事慣性導(dǎo)航方面的研究。E-mail:qinchongbit@126.com

U666.1

A

2095-8110(2016)04-0019-06

猜你喜歡
系統(tǒng)
Smartflower POP 一體式光伏系統(tǒng)
WJ-700無(wú)人機(jī)系統(tǒng)
ZC系列無(wú)人機(jī)遙感系統(tǒng)
基于PowerPC+FPGA顯示系統(tǒng)
基于UG的發(fā)射箱自動(dòng)化虛擬裝配系統(tǒng)開(kāi)發(fā)
半沸制皂系統(tǒng)(下)
FAO系統(tǒng)特有功能分析及互聯(lián)互通探討
連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
一德系統(tǒng) 德行天下
PLC在多段調(diào)速系統(tǒng)中的應(yīng)用
主站蜘蛛池模板: 国产精品熟女亚洲AV麻豆| 午夜人性色福利无码视频在线观看| 日本精品影院| 中文字幕欧美日韩高清| 午夜国产在线观看| 久久这里只有精品8| 四虎精品国产永久在线观看| 亚洲av无码久久无遮挡| 国产亚洲欧美日韩在线一区| 日韩 欧美 小说 综合网 另类| 91精品福利自产拍在线观看| 久久国产精品娇妻素人| 美女无遮挡被啪啪到高潮免费| 免费不卡在线观看av| 她的性爱视频| 亚洲一级色| 亚欧美国产综合| 97精品国产高清久久久久蜜芽 | 97青青青国产在线播放| 日韩在线2020专区| 国产无码在线调教| 四虎精品黑人视频| 在线a网站| 亚洲美女久久| 亚洲色图综合在线| 国产成人精品无码一区二| 国产欧美网站| 狠狠色丁香婷婷综合| 亚洲另类色| 欧美成人二区| 欧洲极品无码一区二区三区| 国产丰满大乳无码免费播放| 黄色三级网站免费| 精品无码人妻一区二区| 国产jizzjizz视频| 亚洲黄网视频| 亚洲美女一级毛片| 欧美精品啪啪一区二区三区| 中文字幕第4页| 99热这里只有成人精品国产| 污网站免费在线观看| 狼友av永久网站免费观看| 国产成人成人一区二区| 亚洲va视频| 制服丝袜在线视频香蕉| 免费无码网站| 亚洲欧美一区二区三区蜜芽| 欧美成人午夜视频免看| 狠狠亚洲婷婷综合色香| 日韩高清成人| 午夜无码一区二区三区| 久久福利网| 日本黄色不卡视频| 久久精品66| www.精品国产| 一区二区无码在线视频| 国产欧美精品专区一区二区| 亚洲人成在线精品| 亚洲视频免费在线| 欧美精品黑人粗大| 欧美一级高清片欧美国产欧美| 少妇高潮惨叫久久久久久| 亚洲欧美日本国产专区一区| 四虎永久在线精品国产免费| 国产福利拍拍拍| 青青青国产视频| 亚洲成a人在线播放www| 在线视频97| 亚洲永久色| 欧美一级一级做性视频| 午夜国产理论| 日本影院一区| 国产精品欧美激情| 国产一区二区免费播放| 日本国产一区在线观看| 亚洲欧洲自拍拍偷午夜色| 亚洲日韩日本中文在线| 免费a级毛片视频| 精品久久国产综合精麻豆| 无码专区国产精品第一页| 欧美成人区| 亚洲一区色|