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

靜不安定飛機縮比模型跨聲速顫振試驗技術

2017-03-01 07:57:30冉玉國李秋彥楊興華
關鍵詞:模態優化模型

冉玉國, 李秋彥, 楊興華

(1.中航工業成都飛機設計研究所, 成都610091;2.中國空氣動力研究與發展中心高速所, 四川綿陽621000)

靜不安定飛機縮比模型跨聲速顫振試驗技術

冉玉國1, 李秋彥1, 楊興華2

(1.中航工業成都飛機設計研究所, 成都610091;2.中國空氣動力研究與發展中心高速所, 四川綿陽621000)

全機縮比模型跨聲速顫振試驗是飛行器氣彈設計驗證工作的重要環節,借助該試驗可研究跨聲速區空氣壓縮性效應,并有效獲得飛行器的跨聲速顫振邊界。通過總結、梳理多次部件級跨聲速顫振試驗研究經驗,并基于某型翼身模態高度耦合的靜不安定飛機全模試驗研究,較全面地闡述了目前全模跨聲速顫振模型設計、地面試驗及風洞試驗技術,可供工程應用借鑒。

跨聲速顫振試驗;全機相似模型;模態耦合;模態試驗;靜不安定飛機

引言

飛行器顫振是機體結構動力特性與非定常氣動力的耦合問題,對飛行器設計起著至關重要的作用。隨著計算機技術與數值分析方法的廣泛發展,在結構方面,有限元方法日趨完善,結構動力學分析技術日漸成熟;在空氣動力方面,CFD技術得到廣泛發展,非定常氣動力數值模擬精度也越來越高、計算速度越來越快。CFD/CSD耦合分析技術的應用[1-4]大大促進了氣彈分析技術的發展,但跨聲速顫振問題仍舊是目前國內氣彈數值計算的難點,尚無在飛機型號設計中經過試驗驗證并得到廣泛認可的成熟分析手段,跨聲速顫振風洞試驗依舊是現階段型號飛機研制中廣泛采用的研究手段[5-7]。

由于跨聲速區域存在激波與混合流場等,跨聲速區非定常氣動力呈現高度的非線性,使得大多數飛機在該區域內的顫振邊界出現“凹坑”現象,因此,對飛機在該區域的非定常氣動力及顫振特性進行準確分析并開展針對性設計,是型號研制中進行低空大表速顫振試飛前的必要環節與前提條件。

國外從20世紀60年代即開始進行飛機部件及全機跨聲速試驗,尤以美俄兩國起步早、研究深入、設備先進以及技術領先為代表,其中,美國NASA Langley TDT風洞[8]先后承擔了B-58、C-141、F-111、C-5、B747、DC-10、F14、F15、B1、F16、B767、X29、C-17、B777、F18E/F等一系列軍民用飛機跨聲速氣彈試驗,包括數十次部件、半模與全模試驗;俄羅斯中央空氣流體動力研究院(TsAGI)T106/T109/T128風洞則承擔了蘇霍伊、米格等設計局多型飛機跨聲速顫振試驗項目。

國內跨聲速顫振試驗研究工作起步較晚,初期借助對外合作項目進行摸索,近年來發展迅速,多個科研院所已成功開展了多次軍民用飛機部件級(平尾、鴨翼、垂尾與機翼等)跨聲速顫振試驗,逐步建立起了完整的模型設計、制造、地面振動試驗及風洞試驗流程,形成了可靠高效的高精度設計與試驗方法。全機級(全模)跨聲速顫振風洞試驗方面,國內科研院所近年來也進行了摸索研究,取得了顯著的研究成果。全模顫振試驗技術包括全模顫振模型結構設計、模型防護系統設計、全模自由支持系統設計、地面試驗、風洞試驗等技術領域,該試驗技術難度大、安全風險高、執行成本高。本文基于某型先進布局飛機全模跨聲速顫振試驗項目實踐,較為全面地闡述了全模顫振試驗各項關鍵技術要求及常用方法,具有較強的工程實用價值,可供參考與借鑒。

1 模型設計技術

跨聲速顫振試驗模型以相似律為準則開展設計,須滿足與模擬對象的外形相似、質量特性相似、剛度分布相似,其中質量特性與剛度分布相似綜合表現為模態相似,包含模態頻率相似與模態振型相似兩方面。為了保證模型在風洞試驗過程中的安全,顫振模型設計需同時考慮其自身強度性能及試驗防護系統設計,對模型實現被動與主動雙重防護。

全機跨聲速顫振模型設計技術包括比例尺選取、結構設計、結構參數優化、質量配置、防護系統設計等一系列技術內容,其目的在于使得模型氣動外形、結構布局、質量剛度和強度特性以及模型支持防護系統實現最佳融合,以最大程度滿足對全尺寸飛機的相似模擬。

1.1 比例尺選取技術

跨聲速顫振模型比例尺選取主要包含尺寸、速壓及密度等三個基本比例尺。

為了提高顫振模型的模擬精度,往往希望將模型尺寸設計得較大,因此需要選擇較大截面尺寸的風洞進行試驗,而國內目前可供型號跨聲速顫振試驗使用的風洞很少,常用的為中國空氣動力研究與發展中心高速所FL26風洞,其試驗段截面尺寸為2.4 m×2.4 m,試驗馬赫數0.3~1.18,可開展固定馬赫數的增壓試驗。尺寸比例尺選取時需遵循截面尺寸比、風洞堵塞度、模型位置等三個方面的要求:翼展一般不超過0.6倍風洞試驗段寬度;堵塞度一般不大于0.02;模型頭部需在流場加速區后,位于均勻流場中,考慮支持裝置對流場的影響。

速壓比例尺選取受設計點馬赫數、試驗風洞速壓帶范圍、飛機顫振特性等三個因素的影響。設計速壓點選擇應綜合考慮試驗成本、試驗安全性、試驗工況(多部件、變參數)等各方面因素,使得模型速壓設計點處于風洞速壓帶的合適位置。

由于型號常選用的FL26風洞為暫沖式半回流風洞,其流場速壓帶窄、速壓下邊界高、載荷大,模型在風洞中將承受較大的氣動載荷且始終處于結構振動環境中,對模型頻率、強度(靜強度、動強度)性能提出了很高的要求,成為可能引起模型產生超重(設計完成后的模型實際重量大于相似律要求的目標重量)的要素之一;而模型的質量比(或密度比)效應對風洞試驗結果存在影響,過大的模型超重(即密度比例尺過大),可能導致試驗結果的不真實性,所以在模型設計中又必須通過優化速壓設計點、模型結構形式、材料和工藝等措施對模型超重進行嚴格控制。多種因素的相互制約增加了跨聲速顫振模型設計的難度與技術風險。

尺寸、速壓與密度等三個基本比例尺確定后,可根據相似推導出的誘導比例尺獲得誘導比例尺,如頻率比例尺、質量比例尺、轉動慣量比例尺、力比例尺、剛度與柔度比例尺等。

1.2 結構設計技術

全機跨聲速顫振模型(簡稱“全模”)結構設計大體上分為三種方法,即動力相似方法、結構相似方法和準結構相似方法。其中,動力相似設計方法通常采用梁架、變厚度板等模擬剛度分布,通過調整結構的尺寸使其滿足振動形態、頻率等要求;結構相似設計方法將真實飛機的每一個受力構件(梁、肋與蒙皮等)按相似比統一縮比,完全模擬結構布置,其優點是能精確反映真實結構的動力和顫振特性,且不易超重,但其設計、制造復雜,成本太高;準結構相似設計方法將受力構件進行適當合并簡化,通過調整結構尺寸使其滿足模態要求。

全模結構設計包含機身、翼面及其支持系統等三個方面。

機身結構常按動力相似要求開展設計,需滿足機身自身動力特性要求、各翼面支持及傳載等功能要求,一種典型結構形式如圖1 所示,由復合材料/金屬主梁、框板、蒙皮及填充硬泡等組成。前機身通常采用簡化的工程梁模擬,截面形式一般為“十”字型,通過調整梁元參數(高度、寬度等)實現機身垂直彎曲、水平彎曲與扭轉剛度的模擬;后機身則常采用變厚度板或框架形式模擬。機身結構設計時需設計簡單、高效的翼面安裝支持結構并開展減重優化,以實現對平尾、垂尾等翼面的有效支持,避免安裝結構與翼面自身產生模態耦合。翼面結構通常按準結構相似方法設計,一般設計為梁架式閉盒結構,包含金屬/復合材料梁架、蒙皮及填充硬泡等構件[9],典型結構形式如圖2 所示。翼身連接支持剛度則通過理論簡化后布置適當的連接元件進行模擬,如機翼根部布置連接彈簧片模擬翼身支持剛度;垂尾或平尾等全動翼面根部布置轉軸模擬機身彎曲支持剛度、彎曲彈簧片模擬機身旋轉支持剛度等。最后,根據設計需要對機身、翼面結構布局形式、參數尺寸與根部支持元件(連接彈簧片)進行優化,結合蒙皮氣彈剪裁,實現剛度模擬、氣動載荷傳遞等功能的同時有效減小模型超重,提高模擬精度。

圖1 機身典型結構形式圖

圖2 翼面典型結構形式圖

全模結構選用比剛度、比強度高的輕質材料,如玻璃纖維、碳纖維、高強度鋁合金等,選材時需考慮制備工藝、成本等各方面的約束,采用技術成熟、工藝可行、性能穩定、性價比高的材料。優化選取模型材料是降低模型超重系數的重要措施之一。

1.3 結構參數優化

采用動力相似方法及準結構相似方法設計的全模結構,按等剛度折算結果設置結構件初始參數后,其固有振動特性(包括頻率與振動形態兩部分)與要求值往往存在一定差異,此時需采用動力學優化方法對結構參數進行優化;此外,跨聲速顫振模型頻率高、在風洞中承受的載荷大,模型需要具有較高的強度,因此,為了保證模型的強度性能,在進行動力學優化的同時需聯合靜強度性能進行綜合優化。

采用復合材料的全模優化參數通常選擇梁架凸緣寬度與厚度、蒙皮鋪層厚度與方向等,可利用“全局-局部”綜合優化算法[10],發揮全局優化對初值的低依賴性及局部優化搜索的快速性等優勢,以固有頻率、固有振型與強度特性為約束,實現全模顫振模型結構參數的快速、有效優化。本文開展的靜不安定全模跨聲速試驗研究中,由于升力體機身及高度翼身融合布局的應用,使得機身模態與各翼面模態耦合嚴重,全機模態密集、復雜,對全模動力優化提出了嚴峻的考驗,應用圖3 所示的典型流程(動力優化部分)進行優化后,全機關鍵模態結果模態頻率與要求值誤差小于5%,模態振型與要求振型吻合良好。

圖3 典型“全局-局部”綜合優化算法

全模經過一輪動力優化設計完成后,其固有模態、剛度分布等滿足設計要求,需對其進行試驗載荷工況下的強度校核優化,若動力優化結果滿足強度要求,則該結果可作為全模系統的設計參數值;若動力優化結果不滿足強度要求,則需改變強度薄弱環節的參數值進行重新優化。作為強度校核優化的輸入,試驗載荷的準確預估是優化環節的關鍵之一。工程中常用的預估方法包括高階面元法、CFD數值模擬法等,為使模型具有合適的安全余量(安全且不過于超重),通常設置載荷計算速壓為風洞試驗最大速壓的1.5倍。本文為了驗證面元法與CFD數值分析方法的準確度,在全模顫振試驗過程中利用機身及翼面上設置的應變片對其所受的非定常氣動力載荷進行了標定與測試,獲得了良好效果。

1.4 質量配置技術

試驗模型完成動力/靜力有限元模型優化后,為了將其中的質量模型由僅具有數學意義的集中質量模型轉化為具有物理意義且可實現的實際結構質量模型,需要對結構進行質量配置設計,利用有限的配重塊調節全模各部件與全機的質量特性,使其滿足設計要求。根據配重設計的情況優化模型的超重比。本文開展的全模試驗中發展了基于數模的顫振模型質量配置技術,可綜合考慮模型各細節、連接件重量及膠接膠量等,設計出的模型質量與實物模型差異1.7%。

1.5 模型防護系統設計技術

跨聲速顫振試驗過程中模型(尤其是全動翼面模型)所受的載荷大,破壞風險高,因而有必要針對全動翼面模型設計可靠、有效的模型應急防護系統對其實施防護,保證試驗順利進行,減小損傷試驗設施的可能性。全動模型應急防護系統目前常采用兩種形式:第一種為主動防護系統,即使用氣壓作動筒推動制動滑槽,對模型根部運動進行卡止限位,增加模型根部旋轉支持剛度,特別使用于顫振速度隨旋轉支持剛度增加而升高的顫振模型,當模型發生顫振時,風洞監測控制系統在發出關車指令的同時,發出電信號激活防護系統運行,高度適用于顫振速度隨旋轉支持剛度增加而升高的顫振模型;第二種為被動防護系統,即采用限幅鋼索在模型展向尖部對模型實施保護,當模型發生顫振時自動限制其振動幅值,配合風洞監測控制系統發出的關車指令,適用于顫振速度隨彎曲支持剛度增加而升高的顫振模型。主動防護系統響應速度快、系統工作可靠、防護效果好,但需通過設計較為復雜的機構實施;被動防護系統實施成本低,但響應速度較慢,限幅的量值需要進行不斷調試優化,防護效果有時達不到預想效果。

2 地面試驗技術

全模地面試驗的目的是通過試驗方法對制造出的實際模型進行檢驗,以確定設計的有效性、準確性,并確定實際模型是否能用于后續風洞試驗。地面試驗主要分為質量特性測量試驗、剛度試驗及模態試驗三部分。質量特性測量試驗經典方法主要包括稱重法測量全模各部件質量;吊線法測定全模重心;擺線法、頻率法等測定全模轉動慣量。全模剛度試驗可通過在設置的加載點上加載、在設置的響應點上測量其位移而得到全模剛度分布特性。

全模模態試驗是整個地面試驗的難點與重點,其結果直接影響后續風洞試驗的有效性、準確性。顫振模型在制造過程中通常難免出現質量特性與剛度特性的允許偏差,因而在模態試驗過程中須通過調節各部件與機身上的連接構件(翼身連接彈簧片、全動翼面根部轉軸軸徑與彎曲彈簧片厚度等)參數對這些偏差進行修正,保證全模各部件的模態特性滿足要求,進而使得模型可用于后續風洞試驗。由于全模中各部件可能存在嚴重、復雜的模態耦合,全機密集模態的識別與調整成為模態試驗的核心,可首先對無耦合或弱耦合的部件模態進行試驗,確定其設計參數,再對耦合嚴重的部件模態進行試驗,對試驗結果進行分析,從而確定關鍵構件參數的優化方案。全模模態試驗通常采用隨機激勵法與純模態法相結合,對易于識別的模態利用隨機激勵法可有效提高試驗速度,而對難于識別的耦合模態,可通過調整激振器位置與數量、激振力幅值與相位等措施開展純模態試驗,提高模態試驗精度與準確度。

尤其值得注意的是在地面試驗之前,需對全模中全動翼面的間隙進行測量,采取必要措施減小間隙量,這往往關系到跨聲速試驗的成敗。全模模態試驗需要在風洞試驗采用的支持系統上進行,以確保模態試驗結果的正確性。

地面試驗完成后,需利用試驗結果對模型的密度、頻率及速壓等比例尺進行修正,用于將后續風洞試驗結果轉換到目標飛機上。

3 全模支持控制系統設計

全機跨聲速顫振模型支持系統需滿足支持頻率、靜/動穩定性及強度特性等要求。為了模擬飛機的自由飛行狀態,全模支持系統常常設計為柔性支持系統,該系統首先需具備足夠的強度特性以保證對全模支持的可靠、穩定與安全。利用該系統對全模實施支持后,模型的剛體模態頻率(沉浮、側擺、俯仰、滾轉與偏航)需滿足低于1/3全模最低彈性模態頻率的要求;模型與支持系統組成的整體系統在風洞試驗過程中既不能發生靜不穩定性也不能發生剛體模態耦合顫振等動不穩定性。由此可知,支持系統的剛度特性直接決定了試驗的安全性與成敗。此外,由于顫振模型在試驗過程中受非定常氣動力作用,模型的姿態(攻角、側滑角、滾轉角等)處于時變狀態,為了保證模型受載較小(攻角、滾轉角等不過大),支持系統需具有根據模型姿態主動調節其攻角、滾轉角的控制功能。

國內目前僅利用中國空氣動力研究與發展中心高速所FL26風洞開展了型號全模顫振試驗摸索研究,借助這些研究初步掌握了全模跨聲速顫振試驗技術、建立了相關試驗流程,在此過程中,設計完善了全模顫振支持控制系統(FSS),有限的試驗顯示該系統支持特性、控制效果均滿足設計要求。

針對本文設計的靜不安定飛機全模特點,分兩階段采取措施提升系統穩定性,即試驗前利用全模CFD計算結果提供較為準確的模型氣動系數及動導數,作為FL26風洞FSS系統(圖4 )初始控制參數;試驗過程中,在優化FSS系統張力與控制參數的同時對全模自身參數進行了調整優化。試驗情況顯示FSS工作有效,模型姿態變化幅值控制效果良好,如圖5 所示。

圖4 FL26風洞FSS系統

圖5 FSS伺服系統典型控制曲線

4 風洞試驗技術

全模經過地面試驗確定設計參數值及性能驗證后即可用于開展風洞試驗。跨聲速顫振試驗方法通常為在一個試驗點(馬赫數)上,穩定風洞流場馬赫數,通過增加總壓的方式增加流場的動壓,直至得到模型的顫振邊界。然而,在實際試驗中,為了保證模型不發生破壞,避免對風洞設備產生損壞,試驗往往僅增加至模型亞臨界顫振狀態,再采用亞臨界預測方法進行顫振邊界預測。為滿足全模顫振試驗運行、測試、防護需求,主要試驗設備包括風洞運行控制系統、電視攝像監控系統、模型響應采集測量系統、亞臨界顫振邊界預測分析系統及風洞應急關車系統等分系統。

全模顫振試驗過程為:首先由風洞運行控制系統按試驗要求調整風洞試驗M數和各增壓段的總壓,同時測量風洞試驗段的流場參數(馬赫數M、速壓q、密度ρ、總溫T),FSS伺服系統對全模姿態進行調節以減小模型攻角變化量等參數;通過安裝在模型部件上的應變電橋、加速度傳感器測量出模型振動響應的時間歷程,然后由亞臨界顫振邊界預測分析系統進行功率譜分析和顫振邊界預測,得到顫振頻率和速壓;電視攝像監控系統可以記錄和觀察試驗過程中模型的振動情況,如果模型發生顫振,風洞應急關車系統可以自動和手動關車,并同步啟動模型制動裝置保護模型。典型風洞試驗流場參數及模型響應如圖6 所示,亞臨界預測曲線如圖7 所示。

顫振模型實際發生顫振的速壓值及氣流密度與設計模型時所用的預設值往往存在偏差,在無法設計多個密度比模型開展試驗的情況下,通常采用密度修正方法對試驗結果進行修正[11],再通過修正后的比例尺將結果換算為飛機顫振速壓與頻率,得到飛機的顫振邊界。

圖6 典型試驗流場參數及模型響應

圖7 顫振試驗亞臨界預測曲線

5 結束語

本文從模型設計、地面試驗、模型支持系統設計及風洞試驗等多方面入手,較全面地介紹了目前全模跨聲速顫振試驗涉及到的各種技術,給出了目前通用的相關技術要求、解決方案等,可作為飛行器部件及全模跨聲速顫振試驗研究的參考與借鑒。

[1]KAPANIA R K,BHARDWAJ M K,REICHENBACH E,etal.Aeroelastic analysis of modern complex wings//Proceeding of 6th Symposium on Multidisciplinary Analysis and Optimization,Bellevue,WA,September 4-6,1996:258-265.

[2]安效民,徐敏,陳士櫓.多場耦合求解非線性氣動彈性的研究綜述.力學進展,2009,39(3):284-298.

[3]馬鐵林,馬東立,張華.大展弦比柔性機翼氣動特性分析.北京航空航天大學學報,2007,33(7):781-784.

[4]盧學成,葉正寅,張陳安.基于ANSYS/CFX耦合的機翼顫振分析.計算機仿真,2010,27(9):88-91.

[5]Stroud W J, Krishnamurthy T.Probabilistic design of a plate-like wing to meet flutter and strength requirements//Proceeding of 43rd AIAA/ASME/ASCE/AHS/ASC Structures,StructuralDynamics,and Materials Conference,Danver,Colorado,April 22-25,2002:1464-1470.

[6]羅時偉,楊鉅生.跨聲速顫振復合材料機翼模型設計研究//第七屆全國氣動彈性學術交流會論文集,浙江湖州,2001:125-130.

[7]錢衛,吳江鵬,趙鐵銘,等.顫振風洞試驗中結構相似動力學模型的設計//首屆全國航空航天領域中的力學問題學術研討會論文集(上冊),成都,2004:153-156.

[8]RIVERA J A,FLORANCE J R.Contributions of transonic dynamic tunnel testing to airplane flutter clearance//Proceeding of 41st Structures, Structural Dynamics, and Materials Conference and Exhibit,Atlanta,GA,April 5-6,2000:1768-1775.

[9]冉玉國,李秋彥,程勇.低超重高強度全復材跨聲速顫振模型設計制造與試驗.應用數學與力學,2014,35(S1):146-150.

[10]冉玉國,李秋彥.一種快速優化方法在顫振試驗模型設計中的應用.振動工程學報,2012,25(S1):339-342.

[11]管德.氣動彈性試驗.北京:北京航空學院出版社,1986.

Review of Transonic Flutter Test Techniques for Statically Unstable Aircraft Scaled Model

RANYuguo1,LIQiuyan1,YANGXinghua2

(1.Chengdu Aircraft Design & Research Institute, Chengdu 610091, China;2.China Aerodynamics Research and Development Center, Mianyang 621000, China)

The transonic flutter wind tunnel test of whole aircraft scaled model plays an important role in engineering certification of AC aeroelastic design, by which the aero compressibility and flutter boundary in transonic region can be investigated. The design experiences of transonic flutter tests for several components are summarized first in this paper, and the test techniques for whole aircraft scaled model are introduced in detail subsequently, which are composed of flutter scaled model design technique, ground tests technique and wind tunnel test technique. This paper is also abstracted based on the engineering project of scaled model design for a statically unstable aircraft, which has massive coupling modes between fuselage and other components. The contents introduced in this paper are supposed to be valuable for engineering application.

transonic flutter test; similar modelof whole aircraft; mode coupling; modal test; statically unstable aircraft

2016-04-06

冉玉國(1981-),男,四川成都人,高級工程師,碩士,主要從事飛行器氣彈設計與試驗方面的研究,(E-mail)rygnuaa@126.com

1673-1549(2017)01-0049-06

10.11863/j.suse.2017.01.09

TB115

A

猜你喜歡
模態優化模型
一半模型
超限高層建筑結構設計與優化思考
房地產導刊(2022年5期)2022-06-01 06:20:14
民用建筑防煙排煙設計優化探討
關于優化消防安全告知承諾的一些思考
一道優化題的幾何解法
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
3D打印中的模型分割與打包
國內多模態教學研究回顧與展望
基于HHT和Prony算法的電力系統低頻振蕩模態識別
主站蜘蛛池模板: 无码aaa视频| 91年精品国产福利线观看久久| 国产91视频免费| 国产一级在线播放| 再看日本中文字幕在线观看| 国产福利在线观看精品| 亚洲色图欧美视频| 亚洲香蕉伊综合在人在线| 亚洲欧洲天堂色AV| 欧美日韩福利| 亚洲最大看欧美片网站地址| AV熟女乱| 国产十八禁在线观看免费| 久久久久中文字幕精品视频| 久久影院一区二区h| 国产青榴视频在线观看网站| 国产激情无码一区二区免费| 亚洲中字无码AV电影在线观看| 亚洲精品国产精品乱码不卞| 激情无码视频在线看| 小蝌蚪亚洲精品国产| 国产97视频在线观看| 精品久久久久久中文字幕女| 精品视频第一页| 国产三级毛片| 色哟哟国产成人精品| 国产经典三级在线| 欧美成人综合视频| 日韩精品一区二区三区视频免费看| 在线观看无码a∨| 欧美日本一区二区三区免费| 无码aaa视频| 在线毛片免费| 99精品这里只有精品高清视频| 亚洲国产高清精品线久久| 2021国产乱人伦在线播放| 色国产视频| 最新国产精品第1页| 国产精品流白浆在线观看| 亚洲欧美激情小说另类| 亚洲乱亚洲乱妇24p| 被公侵犯人妻少妇一区二区三区| 亚洲综合二区| 欧美福利在线观看| 亚洲精品手机在线| 久99久热只有精品国产15| 一级毛片免费不卡在线视频| 亚洲视频免费在线| 国产又粗又爽视频| 亚洲最大看欧美片网站地址| 国产成熟女人性满足视频| 青青青视频免费一区二区| 99re这里只有国产中文精品国产精品 | 婷婷伊人五月| 全部无卡免费的毛片在线看| 国产乱码精品一区二区三区中文| a亚洲天堂| 在线网站18禁| 亚洲国产成人无码AV在线影院L| 亚洲天堂免费观看| 亚洲精品久综合蜜| 国产午夜一级毛片| 国产精品视频系列专区| 97成人在线观看| 无码电影在线观看| 亚洲国产欧美国产综合久久 | 特级欧美视频aaaaaa| 伊伊人成亚洲综合人网7777| 日本在线视频免费| 国产精品福利在线观看无码卡| 无码久看视频| 日本欧美在线观看| 日本一区中文字幕最新在线| 国产免费a级片| 精品三级在线| 午夜三级在线| 亚洲嫩模喷白浆| 久久精品国产91久久综合麻豆自制| 国产男女免费视频| 国产欧美在线观看视频| 亚洲无码视频一区二区三区| 在线人成精品免费视频|