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

多操縱面飛行器控制分配技術的發展及應用

2012-11-06 03:49:23陳勇董新民薛建平王小平劉勤
飛行力學 2012年4期
關鍵詞:分配研究

陳勇, 董新民, 薛建平, 王小平, 劉勤

(空軍工程大學 工程學院, 陜西 西安 710038)

多操縱面飛行器控制分配技術的發展及應用

陳勇, 董新民, 薛建平, 王小平, 劉勤

(空軍工程大學 工程學院, 陜西 西安 710038)

控制分配是先進多操縱面飛行器控制理論及應用的一項關鍵技術。首先介紹了多操縱面布局的典型操縱面,然后從靜態、動態和非線性三個方面分別論述了控制分配技術的最新研究成果,概述了航空航天領域及其他工業領域的應用情況。最后總結了控制分配的幾個關鍵問題,討論了今后可能的研究方向。

飛行控制系統; 多操縱面布局; 控制分配; 級聯控制結構

引言

隨著對飛機機動性、安全性和可靠性要求的不斷提高,現代先進的氣動布局常采用過驅動配置。通過各軸向設計冗余的執行機構,為飛行控制律提供多種組合方式,極大地提高了飛機故障條件下的重構能力和戰場條件下的生存能力。但操縱面的增多使控制耦合度增強,執行器的非線性動態特性更為明顯,如何實現指令的有效分配已經成為多操縱面飛行控制系統設計中必須首先解決的問題。

近年來,控制分配作為向冗余執行器分配控制律指令的關鍵技術,已受到國內外專家的廣泛關注[1]。控制分配器有兩個突出優點[2]:能夠考慮執行器物理約束生成控制指令;操縱面故障時無需重新設計飛行控制律即可實現控制重構。

目前,控制分配技術除了在先進戰斗機、客機、導彈、飛碟、飛艇、再入飛行器等航空航天領域成功應用外,還在艦船、機動車、發動機等工業領域也得到了快速發展。本文在分析典型多操縱面布局的基礎上,構造了帶控制分配器的級聯飛行控制結構,從理論研究和工程應用兩個方面論述了控制分配技術的最新成果,指出了存在的問題和潛在的研究方向。

1 多操縱面布局及分類

隨著氣動布局的不斷發展,多操縱面飛行器除了配置方向舵、副翼、平尾等常規舵面外,還引入多種創新型廣義氣動操縱面,如升降副翼、全動翼梢、嵌入式舵面、前緣襟翼、翼梢小翼、開裂式方向舵、后緣襟副翼、擾流板、推力矢量等。國外典型先進飛行器的操縱面配置如圖1所示。

圖1 國外典型飛行器的操縱面配置

研究表明[3],創新型操縱面的氣動控制效率較常規舵面高,可提供足夠的各向控制力矩,但存在操縱面偏轉造成不利耦合力矩的問題。因此,消除或避免不利耦合力矩,開展綜合控制研究是開發使用新型高效氣動操縱面的關鍵。

2 多操縱面控制結構框架

圖2所示為帶控制分配器的多操縱面飛控系統級聯控制結構[2]。飛行控制律模塊基于虛擬指令設計,控制分配器則通過控制效能映射函數將虛擬指令進行合理分配,經執行機構實現對飛機的控制。

圖2 多操縱面飛控系統級聯控制結構

考慮多操縱面飛行器的仿射非線性模型為:

(1)

式中,x(t)∈Rn為系統的狀態;u(t),δ(t)∈Rm分別為執行器指令的期望值和實際值,通常δ(t)≠u(t);f(·)和fδ(·)分別為飛機和執行器的非線性動態特性;gδ(·)為控制輸入函數。在級聯控制結構中,通常引入飛機的控制力矩或偏轉角加速度作為虛擬控制v∈Rk,直接控制飛機三軸方向的運動,滿足:

v(t)=gu(x,u)u(t)

(2)

gv(x,v)v(t)=gδ(x,δ)δ(t)

(3)

式中,gv(·)為虛擬輸入函數;gu(·)為v(t)與u(t)之間的映射函數。

控制分配問題就是按照預定的優化目標將期望控制v(t)最優地分配到冗余受限執行器u(t)上[4]。對v(t)線性化可將線性控制分配問題描述為:

v(t)=Bu(t),u(t)∈Ω

(4)

式中,B∈Rk×m為行滿秩控制效率矩陣;Ω為執行器位置約束和速率約束構成的凸集:

Ω={δ(t)|δmin≤δ(t)≤δmax,

(5)

(6)

3 控制分配技術的理論研究現狀

控制分配技術經歷了從靜態到動態、從單目標到多目標、從線性到非線性的發展過程,已經延伸出多種分配策略。

3.1 靜態控制分配方法

靜態控制分配主要包括直接分配、廣義逆分配、串接鏈分配等經典方法。

3.1.1直接分配法

直接分配是1992年佛吉尼亞理工大學W C Durham教授基于轉矩可達集提出來的[5],存在計算量大的問題。為改進直接分配算法的解算速度,相繼提出了次優分配法和對邊搜索法,但仍可能得不到最優解。文獻[6]提出了具有共面控制的球面快速搜索算法。文獻[7]通過矩陣空間轉換構造了直接分配的線性規劃形式,保證最優的同時避免了高維控制量的空間計算。李衛琪等[8]通過步進式搜索可達集空間的所有表面提出了相鄰面搜索算法。

直接分配法的優點在于幾何意義直觀、虛擬指令不變向、能完全實現可達集指令,缺點是實時性難以保證。

3.1.2廣義逆分配法

廣義逆分配主要包括偽逆、加權偽逆、再分配偽逆、級聯廣義逆、自適應廣義逆等方法。

偽逆法是早期發展的一種以偏轉量最小為優化目標的控制分配方法。通過F-18仿真表明[9],偽逆法能獲取可達集中13.7%的虛擬指令,即使最優的偽逆映射函數也僅能實現可達集的42.7%。

按照不同的使用側重點,可對各控制面分別賦權值,形成了加權偽逆法。對于控制效能較高的操縱面,可分階段選擇不同的權值系數以降低過早進入飽和狀態的概率,于是出現了變權值自適應廣義逆[10]。文獻[11]考慮操縱面位置和速率約束,基于LMI分別提出了線性變權值和非線性變權值控制分配方案。為了改善分配效率,出現了再分配偽逆方案。

3.1.3串接鏈分配法

串接鏈采用層級思想,利用下一級控制機構產生前一級因飽和未實現的控制指令,如圖3所示。

圖3 串接鏈控制分配方法流程圖

圖3中,Bi為第i組操縱面控制效率矩陣,滿足BiPi=I。但使得BiNiPi≠I,可能導致閉環回路性能惡化,引起操縱面振蕩。在對串接鏈的研究中,文獻[12]分析了作動器速率限制對串接鏈分配的影響。文獻[13]根據Lyapunov定理給出了閉環系統零動態漸近穩定的充分條件。張曙光等[14]結合動態逆研究了推力矢量飛機超機動飛行操縱指令的協調分配問題,驗證了串接鏈的有效性。

串接鏈方法工程易實現、分配方式靈活,在先進推力矢量戰斗機中具有重要的應用價值。

3.2 動態控制分配法

靜態控制分配方法以忽略執行器動態為前提,即執行器頻帶寬度遠高于剛體飛機頻帶寬度。實際上操縱面具有不同的動態特性,操縱效率與理想情況也不盡相同,忽略控制分配和執行器的動態鉸鏈將對控制系統性能產生嚴重的影響[15]。研究表明[16],考慮執行器動態過程可有效提高控制分配的精度,從而提升整個控制系統的性能。

3.2.1頻率加權法

針對操縱面動態性能差異,文獻[17]將操縱面的頻域性能引入優化目標,提出了一種穩定的混合優化動態控制分配器設計方法,采用合理的控制加權,利用鴨翼的動態特性抵消了座艙過載,減小了配平阻力。文獻[18]采用方向保持法對特定控制量進行頻率加權,提出了二次規劃優化動態控制分配方法。

3.2.2動態補償法

在執行器動態響應過程不可忽略的情況下,受執行器速率約束的限制,控制指令常低于且滯后于控制分配器產生的指令。為補償執行器動態特性造成的控制指令衰減和滯后。文獻[15]研究了三種典型執行器動態模型的控制指令補償問題。文獻[16]基于LMI求解最優動態補償策略,表明補償后的控制分配策略能夠抑制導彈控制系統的閉環不穩定。文獻[19]則利用所有執行器的冗余控制效能共同補償舵機的動態影響。

3.2.3動態預測法

模型預測作為處理控制約束和狀態約束的有效方法,能夠直接考慮執行器動態特性建立穩定的控制分配律。文獻[20]研究了再入飛行器的模型預測動態控制分配問題,將執行器理想化為一階線性動態模型,通過預測控制求解包含執行器動態的控制分配指令。文獻[21]分析了模型預測動態控制分配系統的穩定性問題,缺點是對執行機構動力學模型依賴程度較高[1]。文獻[22]提出了一種子空間預測控制分配方法,實現了B747飛機多組關鍵操縱面故障時的容錯控制。馬建軍等[1]針對動態控制分配的不確定性建模,結合子空間預測控制理論實現了ADMIRE無人機的魯棒動態控制分配。

3.3 非線性控制分配法

大多數工業控制系統都是控制量耦合的非線性系統,且執行器本身存在非線性動態,直接針對真實系統進行控制分配,具有重要的研究價值。

3.3.1分段線性規劃法

分段線性規劃方法假設執行器位置與生成的虛擬控制量之間具有分段線性關系。文獻[23]對操縱面控制效率曲線分段線性化,從而將非線性控制分配問題轉化為整數線性規劃。文獻[24]通過利用仿射函數逼近操縱面的非線性控制效能,實現了非線性控制分配。王鵬等[25]運用分段線性函數處理飛翼布局的非線性控制分配問題,提高了控制響應精度。

3.3.2非線性規劃法

非線性規劃方法中最常見的有二次規劃和多目標規劃。文獻[26]針對結構奇異導致非線性控制分配非凸的問題,將消除奇異引入優化控制目標提出一種序列二次規劃分配方法。文獻[27]利用高階多項式擬合氣動數據,提出了一種非線性規劃控制分配方案。文獻[28]綜合權衡多個控制目標提出了一種非線性規劃控制分配方案,給出了多目標非線性控制分配解的評價指標及方法,實現了不同飛行任務的多種目標優化。

3.3.3非線性自適應法

非線性自適應方法根據漸近最優和準最優控制條件,將優化問題轉化為最優集的收斂問題,設計穩定且最優的非線性控制分配律。文獻[29]根據拉格朗日定理設計了非線性最優控制分配方案,并進一步研究了包含執行器動態的自適應控制分配問題[30],可推廣至操縱面故障或損傷時的自動控制重構。文獻[31]提出了一種有限時間收斂的非線性控制分配方案,結合模型參考動態逆控制律證明了閉環控制系統的穩定性。

4 控制分配技術的工程應用現狀

4.1 航空航天領域

化學示蹤劑以無機鹽類、熒光染料類、鹵代烴類為主,主要用于研究儲層物性、油層連通性和注水對地層的影響。由于化學示蹤劑存在相對用量大、需要井口作業、成本高和環境污染等問題,所以該類示蹤劑呈淘汰趨勢。

控制分配技術在航空航天領域的應用涉及飛機、導彈、飛碟、再入飛行器、衛星、飛艇等。

在飛機控制方面的研究成果最多。2003年,NASA和美國空軍在B747和C-17上研究了運輸類飛機關于損傷自適應性的控制再分配策略,改善了飛行品質。Wright實驗室進行了F-16推力矢量飛機的全包線非線性仿真,驗證了偽逆法的有效性。文獻[32]在新型控制面技術驗證機LMICE上運用偽逆分配驗證了飛行控制系統的有效性。文獻[33]將控制分配技術用于無尾飛機縱向控制系統中。

在導彈控制方面也有少量應用案例。文獻[34]分別通過偽逆法和串接鏈來協調導彈的推力矢量和傳統控制面的偏轉。文獻[35]研究了帶尾翼和燃氣舵的導彈動態控制分配問題,驗證了約束和無約束控制分配方法均可實現導彈的最優控制。

航天方面,文獻[36]采用整數規劃來綜合分配氣動操縱面和反作用控制系統,并運用Lyapunov理論證明了高度控制系統的穩定性。文獻[37]研究了下一代再入飛行器的控制分配問題。文獻[38]將控制分配技術應用于飛艇的縱向和側向控制中。

4.2 其他工程領域

控制分配技術在航空航天領域的成功應用,引導其在艦船、機動車等工程領域得到了快速發展。

在機動車控制方面,控制分配技術根據橫擺角速度、側偏角與參考值的誤差,考慮執行器約束和機輪工作狀態將控制力矩分配到機動車各車輪。文獻[41]運用非線性自適應控制分配策略實現了電動車的穩定控制。文獻[42]基于加權偽逆和滑模控制實現了汽車的動態協調控制,提出了一種自適應加速收斂的定點控制分配方法。

5 存在問題和研究方向

盡管多操縱面控制分配技術在理論研究和工程應用方面已取得了諸多成果,但仍然存在許多待完善的地方,需要在發展多操縱面飛行器的過程中展開進一步的深入研究,主要包括:

(1)現代先進飛行器普遍包含多個操縱面,尤其是大型軍用運輸機和民用客機,操縱面甚至在30副以上。由于各種操縱面具有不同的物理屬性,對飛機三軸將產生不同的控制效能,不合理的控制賦權容易出現操縱面間效能抵消的現象,導致控制分配對可達集空間分配效率降低,直接影響控制性能。如何考慮執行器物理特性以確定控制權值是提高控制分配效能的重要研究方向。

(2)關于包含執行器動態特性的控制分配方法,目前國內外公開發表的文獻不多,絕大部分都是集中研究典型執行器精確模型的動態補償問題。實際上執行機構作為影響飛行控制系統可靠性和安全性的關鍵部件之一,大量參數存在非線性和不確定性,常規方法難以建立其精確的數學模型。當執行器出現損傷等意外故障時,模型的動態不確定性將更為突出。考慮不確定執行器物理特性設計動態控制分配器是當前重要研究方向之一,越來越得到研究人員的高度重視。

(3)與常規飛機不同的是,多操縱面飛機包含多組性能參數各異的約束控制面,其控制效能在擴展的飛行包線內隨高度、馬赫數等時變參數不斷變化,直接進行控制器設計必然增加參數調節的難度。針對多操縱面飛行器隨狀態參數攝動的線性參變數學模型,結合操縱面氣動特性設計大包線范圍內穩定的飛行控制律具有十分重要的現實意義。

(4)受外界環境和自身參數影響,飛行器模型不可避免地存在未建模動態,基于精確模型設計的飛行控制律和控制分配器將具有一定的保守性。結合魯棒控制理論設計閉環穩定的級聯飛行控制器,以抑制模型攝動的不利影響,是提高飛控系統控制性能的重要手段之一。

(5)隨著智能控制技術的發展,啟發式優化算法在飛行控制領域已得到了初步應用。如何將遺傳算法、蟻群算法、粒子群優化等新興方法融入到控制分配優化問題中,把智能控制與控制分配結合起來,是未來控制分配研究的另一個熱點。

(6)許多控制分配方案都需解決計算量大、設計復雜、閉環不穩定的問題,工程應用有待進一步開拓。需要指出,盡管目前已有若干工程應用成功的例子,但更多的方案仍停留在數值驗證階段。

6 結束語

多操縱面控制分配技術是當前先進飛行控制理論及應用研究方面的前沿課題,是面向冗余控制系統設計的一項重要的關鍵技術,在航空、航天和其它工業領域具有廣泛的應用前景。本文在分析典型多操縱面布局的基礎上,概述了近年來控制分配技術的發展現狀。當然,控制分配理論和應用的研究尚未成熟,還存在一系列尚待解決的問題,值得未來在發展多操縱面飛行器的過程中做進一步探索。

[1] 馬建軍,鄭志強,胡德文.包含執行器動力學的子空間預測動態控制分配方法[J].自動化學報,2010,36(1):130-138.

[2] 陳勇,董新民,薛建平,等.兩類混合優化控制分配算法的靈敏度分析[J].飛行力學,2011,29(2):49-53.

[3] 王元元,張彬乾,沈冬.W型無尾布局復合式氣動舵面設計研究[J].西北工業大學學報,2008,26(6):698-702.

[4] 馬建軍,李文強,李鵬,等.飛行器控制分配技術研究現狀與展望[J].飛行力學,2009,27(3):1-5.

[5] Durham W C.Constrained control allocation three moment problem [J].Journal of Guidance,Control,and Dynamics,1994,17(2):330-336.

[6] Petersen J A,Bodson M.Fast implementation of direct allocation with extension to coplanar controls[J].Journal of Guidance,Control,and Dynamics,2002,25(3):464-473.

[7] Bodson M.Evaluation of optimization methods for control allocation[J].Journal of Guidance,Control,and Dynamics,2002,25(4):703-711.

[8] 李衛琪,魏晨,陳宗基.受限控制直接分配新算法[J].北京航空航天大學學報,2005,31(11):1177-1180.

[9] Bordignon K A,Durham W C.Closed-form solutions to constrained control allocation problem[J].Journal of Guidance,Control,and Dynamics,1995,18(5):1000-1007.

[10] 陳勇,董新民,薛建平,等.多操縱面飛控系統約束自適應控制分配策略[J].系統工程與電子技術,2011,33(5):1118-1123.

[12] Berg J,Hammett K,Schwartz C,et al.An analysis of the destabilizing effect of daisy chained rate-limited actuators[J].IEEE Transactions on Control Systems Technology,1996,4(2):171-176.

[13] Buffington J M,Enns D F.Lyapunov stability analysis of daisy chain control allocation[J].Journal of Guidance,Control,and Dynamics,1996,19(6):1226-1230.

[14] 張曙光,高浩.推力矢量飛機操縱指令的協調分配[J].西北工業大學學報,1996,14(3):361-365.

[15] Oppenheimer M,Doman D.A method for compensation of interactions between second-order actuators and control allocators[C]//2005 IEEE Aerospace Conference.Malta,2005:1-8.

[16] Kishore W C,Sen S,Ray G.Dynamic control allocation for tracking time-varying control demand[J].Journal of Guidance,Control,and Dynamics,2008,31(4):1150-1157.

[17] H?rkeg?rd O.Dynamic control allocation using constrained quadratic programming [J].Journal of Guidance,Control,and Dynamics,2004,27(6):1028-1034.

[18] Simmons A T,Hodel A S.Control allocation for the X-33 using existing and novel quadratic programming techniques[C]//Proceedings of the 2004 American Control Conference.Boston,2004:1701-1706.

[19] 陳勇,董新民,薛建平,等.舵機動態系統耦合補償控制分配策略[J].儀器儀表學報,2010,31(7):1495-1500.

[20] Yu Luo,Doman D.Model predictive dynamic control allocation with actuator dynamics[C]//Proceedings of the 2004 American Control Conference.Boston,2004:1695-1700.

[21] Vermillion C,Jing S,Butts K.Model predictive control allocation for overactuated systems stability and performance[C]//IEEE Conference on Decision and Control.New Orleans,2007:1251-1256.

[22] Hallouzi R,Verhaegen M.Reconfigurable fault tolerant control of a Boeing 747 using subspace predictive control[C]//Guidance,Navigation and Control Conference and Exhibit.Carolina,2007:1-18.

[23] Bolender M,Doman D.Non-linear control allocation using piecewise linear function[J].Journal of Guidance,Control,and Dynamics,2004,27(6):1017-1027.

[24] Doman D,Oppenheimer M. Improving control allocation accuracy for nonlinear aircraft dynamics[R].AIAA-2002-4667,2002.

[25] 王鵬,周洲,王睿.基于分段線性的飛翼布局飛機控制分配方法[J].西北工業大學學報,2009,27(3):321-325.

[26] Johansen T,Fossen T,Berge S.Constrained nonlinear control allocation with singularity avoidance using sequential quadratic programming [J].IEEE Transactions on Control Systems Technology,2004,12(1):211-216.

[27] Poonamallee V,Yurkovich S,Serrani A,et al.A nonlinear programming approach for control allocation[C]//Proceedings of the 2004 American Control Conference.Boston,2004:1689-1694.

[28] 楊恩泉,高金源,李衛琪.多目標非線性控制分配方法研究[J].航空學報,2008,29(4):995-1001.

[29] Tj?nn?s J,Johansen T A.Adaptive optimizing nonlinear control allocation[C]//The 16th IFAC World Congress.Prague,2005:1-6.

[30] Tj?nn?s J,Johansen T A.Optimizing adaptive control allocation with actuator dynamics[C]//Proceedings of the 46th IEEE Conference on Decision and Control.New Orleans,2007:3780-3785.

[31] Liao F,Lum K Y,Wagn J L,et al.Constrained nonlinear finite-time control allocation[C]//Proceedings of the 2007 American Control Conference.New York,2007:3801-3806.

[32] Raney D L,Montgomery R C.Flight control using distributed shape-change effecter arrays[R].AIAA-2000-1560,2000.

[33] 陳懷民,徐奎,馬松輝,等.控制分配技術在無尾飛機縱向控制系統中的應用研究[J].西北工業大學學報,2007,25 (2):199-203.

[34] Lee H,Choi J,Song T,et al.Agile missile autopilot design via time-varying eigenvalue assignment[C]//International Conference on Control,Automation,Robotics and Vision.Kunming,2004:1832-1837.

[35] Ridgely D B,Drake D,Triplett L,et al.Dynamic control allocation of a missile with tails and reaction jets[R].AIAA-2007-6671,2007.

[36] Doman D B,Gamble B J,Ngo A D.Quantized control allocation of reaction control jets and aerodynamic control surfaces[J].Journal of Guidance,Control,and Dynamics,2009,32(1):13-24.

[37] Shertzer R H.Control allocation for the next generation of entry vehicles[R].AIAA-2002-4849,2002.

[38] Battipede M,Gili P A,Lando M.Control allocation system for an innovative remotely-piloted airship[C]//AIAA Atmospheric Flight Mechanics Conference and Exhibit.Rhode Island,2004:1-13.

[39] Fossen T I,Johansen T A.A survey of control allocation methods for ships and underwater vehicles[C]//14th IEEE Mediterranean Conference on Control and Automation.Ancona,Italy,2006:1-6.

[40] Lindegaard K,Fossen T.Fuel-efficient rudder and propeller control allocation for marine craft:experiments with a model ship[J].IEEE Transaction on Control Systems Technology,2003,11(6):850-862.

[41] Tj?nn?s J,Johansen T A.Stabilization of automotive vehicles using active steering and adaptive brake control allocation[J].IEEE Transactions on Control Systems Technology,2010,18(3):545-558.

[42] Junmin Wang,Solis J M,Longoria R G.On the control allocation for coordinated ground vehicle dynamics control systems[C]//Proceeings of the 2007 American Control Conference.New York,2007:5724-5729.

(編輯:方春玲)

Developmentandapplicationofcontrolallocationformulti-effectorvehicle

CHEN Yong, DONG Xin-min, XUE Jian-ping, WANG Xiao-ping, LIU Qin

(Engineering Institute, Air Force Engineering University, Xi’an 710038, China)

Control allocation is one of the key technologies in the control theory and application of the advanced multi-effector vehicle. Firstly, the typical effectors for multi-effector configuration are presented in this overview. Then, the latest research results of control allocation are discussed respectively from three aspects of static, dynamic and nonlinear control allocation., which mainly concentrate on the applications in the aeronautic, astronautic and other industries. Finally, several key issues in control allocation are summarized, and the future potential research directions are discussed.

flight control system; multi-effector configuration; control allocation; cascaded control structure

V249.1

A

1002-0853(2012)04-0289-06

2011-10-08;

2012-02-22

陳勇(1984-),男,四川德陽人,博士研究生,研究方向為控制分配、故障診斷及控制重構。

猜你喜歡
分配研究
基于可行方向法的水下機器人推力分配
FMS與YBT相關性的實證研究
2020年國內翻譯研究述評
遼代千人邑研究述論
應答器THR和TFFR分配及SIL等級探討
視錯覺在平面設計中的應用與研究
科技傳播(2019年22期)2020-01-14 03:06:54
遺產的分配
EMA伺服控制系統研究
一種分配十分不均的財富
績效考核分配的實踐與思考
主站蜘蛛池模板: 综合色88| 无码人妻免费| 国产精品真实对白精彩久久| 久久96热在精品国产高清| 亚洲欧美成aⅴ人在线观看| 欧美精品在线观看视频| 91精品小视频| 国产精品免费p区| а∨天堂一区中文字幕| 亚洲综合第一页| 人妻夜夜爽天天爽| 福利小视频在线播放| 国产欧美视频在线观看| 丰满少妇αⅴ无码区| 国产一区二区视频在线| 一级毛片免费观看久| 亚洲精品不卡午夜精品| 国产成人成人一区二区| 亚洲日韩高清无码| 99成人在线观看| 欧美精品一区二区三区中文字幕| 一级毛片在线播放免费| 国产精品爽爽va在线无码观看| 国产乱子伦精品视频| 日韩毛片在线视频| 国产无码精品在线播放 | 国产精品福利尤物youwu| AV天堂资源福利在线观看| 国产成人精品高清不卡在线| 美女一区二区在线观看| 国产精品综合久久久 | 热这里只有精品国产热门精品| 色男人的天堂久久综合| 98超碰在线观看| 亚洲免费福利视频| 国产jizzjizz视频| аv天堂最新中文在线| 91麻豆国产视频| 欧美视频免费一区二区三区| 精品国产欧美精品v| 波多野结衣一二三| 久草青青在线视频| 午夜日b视频| 天天色天天综合网| 中文字幕人妻无码系列第三区| 国产精品美女自慰喷水| 欧美日本在线观看| 亚洲毛片网站| 久久网欧美| 亚洲国产亚综合在线区| 欧美成在线视频| 最近最新中文字幕在线第一页| 成年人国产视频| 欧美a级完整在线观看| 国产色婷婷| 91色在线观看| 午夜丁香婷婷| 欧美在线三级| 免费xxxxx在线观看网站| 91精品国产情侣高潮露脸| 欧美不卡视频一区发布| 久久精品免费看一| 国产成人91精品免费网址在线| 欧美成人怡春院在线激情| 三级国产在线观看| 亚洲综合精品香蕉久久网| 欧洲欧美人成免费全部视频| 欧美成a人片在线观看| 最新国产成人剧情在线播放| 欧美亚洲国产精品第一页| 亚洲成a∧人片在线观看无码| 日本午夜影院| 久久国产精品国产自线拍| 亚洲熟妇AV日韩熟妇在线| 国产真实乱了在线播放| 亚洲国产综合自在线另类| 日本一本正道综合久久dvd| 日本一区高清| 中文字幕乱妇无码AV在线| 日韩精品高清自在线| 高清无码手机在线观看| 亚洲电影天堂在线国语对白|