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

    柔性后緣可變形機翼氣動彈性分析方法研究

    2022-12-27 11:45:42趙仕偉闞梓李道春
    航空科學技術 2022年12期
    關鍵詞:氣動彈性陣風后緣

    趙仕偉,闞梓,李道春

    北京航空航天大學,北京 100191

    可變形機翼可以根據(jù)飛行環(huán)境和飛行任務的變化光滑、連續(xù)地改變自身形狀,使飛機在整個任務周期內(nèi)的總體性能達到最優(yōu)。然而其光滑連續(xù)變形的氣動外形特性導致非定常氣動力變化劇烈,需采用高精度氣動力分析方法。傳統(tǒng)計算流體力學(CFD)氣動分析方法[1]相比面元法[2]和降階方法計算精度高,但計算量大、耗時多,不便于開展氣動彈性分析。氣動力模型降階可建立一種精度高、適應性強的可變形機翼非定常氣動力模型,用于可變形機翼氣動彈性分析計算。

    目前,常見的有兩類降階非定常空氣動力學建模方法:流場特征分析方法和系統(tǒng)辨識方法。流場特征分析方法針對整個流場進行模型降階分析,得到流場的特征狀態(tài)參數(shù);包括正交分解(POD)方法[3-6]和諧波平衡(HB)方法[7-9]。基于流體計算模型的輸入響應獲得正交基,模型的狀態(tài)量可視為該組正交基的線性疊加,可以通過增加模型階數(shù)從而更準確地模擬實際流場運動特性。但該方法通常用于流場的分析,由于模型考慮整個流場狀態(tài)特性,模型較為復雜、階數(shù)較高,不便于開展氣動彈性分析。

    系統(tǒng)辨識方法對流場的某些特定參數(shù)(如升力系數(shù)、力矩系數(shù)等)進行系統(tǒng)特征識別。非定常氣動力系統(tǒng)辨識常見的有Volterra級數(shù)方法[10-11]、線性狀態(tài)空間模型[12-13]、自回歸移動平均(ARMA)模型[14]、基于人工神經(jīng)網(wǎng)絡的替代模型[15-16]、徑向基函數(shù)(RBF)插值和Kriging模型[17-20]、支持向量機(SVM)方法[21-22]、狀態(tài)觀測器[23]、面向塊的Wiener 模型[24-25]等。

    采用非定常氣動力降階模型,可以提高柔性后緣可變形機翼的氣動彈性響應計算效率,常見的模型降階方法有帶外輸入的自回歸模型(ARX)方法和Volterra 級數(shù)方法。ARX方法的非定常氣動力模型階數(shù)相比Volterra級數(shù)方法要低一個量級左右,計算效率比Volterra 級數(shù)高。因此,本文采用ARX方法建立三維非定常氣動力,與結(jié)構(gòu)狀態(tài)空間模型相耦合,得到柔性后緣可變形機翼氣動彈性系統(tǒng)的狀態(tài)空間模型,在此基礎上計算柔性后緣可變形機翼的顫振特性和陣風響應。

    1 氣動彈性狀態(tài)空間模型

    非定常氣動力ARX模型可以寫成如下表達式

    式中,y(k)表示第k步的廣義氣動力系數(shù),na、nb是系統(tǒng)廣義氣動力系數(shù)輸入和廣義位移輸入階數(shù),可以將該系統(tǒng)改寫成離散狀態(tài)空間的形式

    式中

    狀態(tài)矢量如下

    對于多自由度系統(tǒng),通過歸一化可以轉(zhuǎn)換為模態(tài)坐標運動,結(jié)構(gòu)模態(tài)坐標下的動力學方程可表示如下

    其中,M為系統(tǒng)質(zhì)量矩陣,歸一化處理后是單位矩陣;C為系統(tǒng)阻尼矩陣;K為該多自由度結(jié)構(gòu)系統(tǒng)剛度矩陣;x,x?,x?分別為該系統(tǒng)位移、速度、加速度。各階模態(tài)對應的廣義氣動力FA為動壓q、壓力系數(shù)Cp和結(jié)構(gòu)模態(tài)分析得到的各階振型矢量Si的乘積在各個氣動表面元上的積分

    該系數(shù)可以Fluent 數(shù)值仿真獲取。通過編寫UDF 程序,在每一個時間步輸出對應的各階模態(tài)對應的廣義氣動力,從而得到各階模態(tài)輸入下的各階模態(tài)對應的廣義氣動力。

    將上述模態(tài)坐標下的動力學方程寫成狀態(tài)空間表達形式,可以得到

    式中

    狀態(tài)矢量xs(t)為

    將上述連續(xù)時間系統(tǒng)狀態(tài)空間結(jié)構(gòu)模型進行離散化處理,得到對應的離散狀態(tài)空間模型如下

    其中:

    式中,T為時間步長,上述過程為系統(tǒng)的結(jié)構(gòu)動力學模型的獲取流程,廣義氣動力為系統(tǒng)輸入,結(jié)構(gòu)彈性變形為系統(tǒng)輸出。

    通過ARX 方法建立離散形式的氣動力狀態(tài)空間模型如下

    式中,u(k)為系統(tǒng)的廣義結(jié)構(gòu)位移;yA(k)為系統(tǒng)的廣義氣動力。

    將基于ARX 方法建立的氣動力狀態(tài)空間模型和結(jié)構(gòu)狀態(tài)空間模型相耦合,可以得到整個氣動彈性系統(tǒng)的狀態(tài)空間模型,可以用于分析系統(tǒng)的氣動彈性特性,計算系統(tǒng)的氣動彈性響應。氣動彈性系統(tǒng)狀態(tài)空間模型表達式如下

    上述狀態(tài)空間模型可以通過氣動彈性響應計算或者矩陣特征值來分析該氣動彈性系統(tǒng)的穩(wěn)定性。通過調(diào)整速度和調(diào)整動壓,從而獲得對應速度下的系統(tǒng)響應,當響應幅值隨著時間變化保持不變,該速度對應該馬赫數(shù)下的顫振速度。也可以求解狀態(tài)空間的特征值來分析該系統(tǒng)的穩(wěn)定性。對于離散系統(tǒng),當所有特征值模長均小于等于1時,該系統(tǒng)穩(wěn)定;若至少有一個特征值模長大于1,該系統(tǒng)不穩(wěn)定??梢詫㈦x散狀態(tài)空間模型轉(zhuǎn)化為連續(xù)狀態(tài)空間模型來分析系統(tǒng)穩(wěn)定性,當且僅當所有系統(tǒng)的特征值實部小于0,該系統(tǒng)穩(wěn)定。

    通過ARX 方法同樣可以建立離散形式的陣風氣動力狀態(tài)空間模型,如下所示

    式中,xg(k)為陣風輸入;u(k)是系統(tǒng)的廣義結(jié)構(gòu)位移;yg(k)是系統(tǒng)的陣風廣義氣動力。

    將基于ARX 方法建立的氣動彈性狀態(tài)空間模型和陣風響應氣動力狀態(tài)空間模型相耦合,可以得到考慮陣風響應氣動彈性系統(tǒng)的狀態(tài)空間模型,用于分析系統(tǒng)的陣風響應。氣動彈性系統(tǒng)狀態(tài)空間表達式如下

    2 氣動彈性特性分析方法驗證

    為驗證基于ARX 方法建立的氣動力降階模型的有效性,本文采用Agard445.6 機翼模型。Agard445.6 幾何形狀如圖1 所示,機翼根部弦長0.559m,半展長0.762m,展弦比為1.65,梢根比為0.66,后掠角為45°,翼型為NACA 65A004。設置結(jié)構(gòu)網(wǎng)格最大尺寸為10mm,對機翼根部所在的平面使用固支約束,對該結(jié)構(gòu)模型進行模態(tài)分析。圖2為Agard445.6 機翼對應的前四階模態(tài)振型,前四階模態(tài)的振動頻率與實驗結(jié)果對比見表1,最大誤差為3.15%,有限元與試驗結(jié)果[26]吻合較好。

    表1 前四階模態(tài)的振動頻率與試驗結(jié)果對比Table 1 Comparison between the vibration frequency of the first four modes and the experimental results

    圖1 Agard445.6幾何形狀Fig.1 Geometry of Agard445.6

    圖2 Agard445.6機翼前四階模態(tài)振型圖Fig.2 The first four modal shapes of Agard445.6 wing

    圖3 所示為Agard445.6 機翼非結(jié)構(gòu)網(wǎng)格并導入Fluent中計算。遠場邊界位于距機翼20倍翼根弦長處,并設置為壓力遠場條件。機翼表面設置為無滑動靜態(tài)壁面條件,網(wǎng)格在翼型附近被細化。采用動態(tài)網(wǎng)格結(jié)合用戶自定義函數(shù)(UDF)接口編程,采用彈簧方法進行網(wǎng)格重構(gòu)實現(xiàn)機翼的模態(tài)位移連續(xù)變形。在數(shù)值計算中,采用雙精度求解器進行數(shù)值模擬,采用Spalart-Allmaras湍流模型。

    圖3 Agard445.6機翼氣動網(wǎng)格Fig.3 Aerodynamic grid of Agard445.6 wing

    圖4所示為顫振馬赫數(shù)和參考文獻[24]的對比,可以看出,二者的數(shù)值和趨勢均具有良好的一直性,最大誤差小于8%。

    圖4 顫振馬赫數(shù)對比驗證Fig.4 Flutter Mach number verification

    圖5所示為動壓對應的馬赫數(shù)為0.29時Agard445.6機翼各階模態(tài)時域響應,從圖5中可以看出,各階模態(tài)位移隨著時間變化位移逐漸減少呈收斂趨勢,說明在該動壓下該機翼尚未達到顫振。通過分析離散狀態(tài)空間模型可以看出,最大特征值模長為0.994,小于1,該系統(tǒng)穩(wěn)定。通過分析連續(xù)狀態(tài)空間模型的特征值也可以看出,在該狀況下所有特征值的實部均小于0,該系統(tǒng)穩(wěn)定。圖6所示為動壓對應的馬赫數(shù)為0.3時Agard445.6機翼各階模態(tài)時域響應,從圖6 中可以看出,各階模態(tài)位移隨著時間變化位移基本保持不變,說明動壓對應的馬赫數(shù)為0.3 時非常接近顫振速度。圖7所示為動壓對應的馬赫數(shù)為0.31時Agard445.6機翼各階模態(tài)時域響應,從圖7中可以看出,各階模態(tài)位移隨著時間變化位移逐漸增大,說明該馬赫數(shù)已經(jīng)超過顫振速度了。通過分析離散狀態(tài)空間模型可以看出,最大特征值模長為1.0003,該系統(tǒng)發(fā)散。通過分析連續(xù)狀態(tài)空間模型的特征值也可以看出,在該狀況下所有特征值存在特征值0.2917±91.8737i實部大于0,該系統(tǒng)發(fā)散。

    圖5 Ma=0.29時Agard445.6機翼各階模態(tài)時域響應Fig.5 Time domain response of various modes of Agard445.6 wing at Ma=0.29

    圖6 Ma=0.3時Agard445.6機翼各階模態(tài)時域響應Fig.6 Time domain response of various modes of Agard445.6 wing at Ma=0.3

    圖7 Ma=0.31時Agard445.6機翼各階模態(tài)時域響應Fig.7 Time domain response of various modes of Agard445.6 wing at Ma=0.31

    3 柔性后緣可變形機翼氣動彈性分析計算

    3.1 顫振分析計算

    本節(jié)考慮柔性后緣機翼在不同偏角下的顫振速度計算,僅考慮氣動外形變化對氣動彈性特性的影響,柔性后緣可變形機翼的流體網(wǎng)格如圖8(a)所示。弦向變化范圍為翼型60%到后緣部分,這一段翼型中弧線變形采用作拋物線軌跡變彎,定義60%翼型處中弧線所在的位置和變形后的后緣點的連線與初始翼型中弧線的夾角為后緣偏角β,向下偏為正偏轉(zhuǎn)方向。

    為節(jié)約網(wǎng)格繪制時間,對不同柔性后緣偏角情況采用同一套流體網(wǎng)格進行CFD 仿真計算。在計算某一特定偏角時,可以基于UDF程序預先將該部分后緣偏轉(zhuǎn)至給定的角度,圖8(b)所示為機翼柔性后緣偏轉(zhuǎn)6°時流體網(wǎng)格的變化圖。等到氣動力響應基本穩(wěn)定后,輸入對應的各階位移激勵,從而得到對應偏角下的各階廣義氣動力系數(shù)響應。

    圖8 機翼柔性后緣偏轉(zhuǎn)前后網(wǎng)格變化Fig.8 Grid changes before and after flexible trailing edge deflection

    在訓練馬赫數(shù)為0.9 時,后緣偏角為0°、2°、4°、6°時的顫振馬赫數(shù)分別為0.202、0.203、0.211、0.224??梢钥闯觯S著后緣偏角的增大,顫振馬赫數(shù)相比無偏角情況有一定的提高。

    3.2 陣風響應分析計算

    取馬赫數(shù)為0.7 工況下計算1-cos 陣風響應,圖9 所示為陣風尺寸分別為60m、90m、140m 時柔性后緣偏角為6°下的機翼各階模態(tài)位移響應??梢钥闯?,在該工況下,飛機遭遇陣風響應隨著陣風尺寸的增大,模態(tài)位移響應呈下降趨勢。

    圖9 模態(tài)位移時域響應Fig.9 Time domain response of modal displacement

    4 結(jié)論

    本文提出一種基于氣動降階模型建立柔性后緣可變形機翼氣動彈性分析方法,應用于柔性后緣可變形機翼的顫振特性和陣風響應分析計算。通過結(jié)構(gòu)有限元方法得到各階模態(tài),將各階模態(tài)變形導入CFD 計算,得到對應的氣動力,基于ARX 方法得到氣動狀態(tài)空間模型,耦合結(jié)構(gòu)狀態(tài)空間模型從而建立可變形機翼的氣動彈性狀態(tài)空間模型,用于可變形機翼氣動彈性顫振特性分析和陣風響應計算。研究結(jié)果表明,隨著后緣偏角的增大,顫振馬赫數(shù)相比無偏角情況有一定的提高。在馬赫數(shù)為0.7的工況下,隨著陣風尺寸的增大,模態(tài)位移響應呈下降趨勢。本文尚未考慮柔性后緣結(jié)構(gòu)變形對氣動彈性特性的影響,可在后續(xù)工作中綜合考慮結(jié)構(gòu)氣動影響的氣動彈性特性。

    猜你喜歡
    氣動彈性陣風后緣
    陣風戰(zhàn)斗機
    軍事文摘(2021年19期)2021-10-10 13:29:36
    法國陣風戰(zhàn)斗機
    軍事文摘(2021年17期)2021-09-24 01:11:46
    陣風勁吹
    航空世界(2018年12期)2018-07-16 08:34:50
    機翼后緣連續(xù)變彎度對客機氣動特性影響
    柔性后緣可變形機翼氣動特性分析
    飛翼無人機嗡鳴氣動彈性響應分析
    TNF-α和PGP9.5在椎體后緣離斷癥軟骨終板的表達及意義
    模態(tài)選取對靜氣動彈性分析的影響
    航空學報(2015年4期)2015-05-07 06:43:37
    直升機的氣動彈性問題
    航空學報(2015年4期)2015-05-07 06:43:25
    大型風力機整機氣動彈性響應計算
    a级毛片在线看网站| 日日啪夜夜爽| 下体分泌物呈黄色| 能在线免费看毛片的网站| 日本爱情动作片www.在线观看| 成人国产麻豆网| 午夜福利网站1000一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 99热6这里只有精品| 久久久精品94久久精品| 久久国产乱子免费精品| 黑丝袜美女国产一区| 日韩成人伦理影院| 成年av动漫网址| 午夜福利影视在线免费观看| 国国产精品蜜臀av免费| 最新中文字幕久久久久| 亚洲第一区二区三区不卡| 成人综合一区亚洲| 涩涩av久久男人的天堂| 亚洲精品成人av观看孕妇| 美女cb高潮喷水在线观看| 热re99久久国产66热| 综合色丁香网| 男的添女的下面高潮视频| 亚洲久久久国产精品| av专区在线播放| 日本欧美视频一区| 亚洲内射少妇av| 国产男人的电影天堂91| 亚洲美女视频黄频| 国产 一区精品| 看非洲黑人一级黄片| 伊人亚洲综合成人网| 亚洲精品色激情综合| 蜜桃久久精品国产亚洲av| 国产片特级美女逼逼视频| 亚洲色图综合在线观看| 国产亚洲欧美精品永久| 少妇丰满av| 精品国产一区二区久久| 国产av码专区亚洲av| 最近中文字幕2019免费版| 日本免费在线观看一区| 国产亚洲5aaaaa淫片| 欧美精品一区二区免费开放| 久久精品国产自在天天线| 久久久久久伊人网av| 亚洲精品乱久久久久久| 国产国拍精品亚洲av在线观看| 亚洲av福利一区| 亚洲欧洲日产国产| 久久精品久久精品一区二区三区| 日韩强制内射视频| a级毛片免费高清观看在线播放| 中文精品一卡2卡3卡4更新| 欧美日韩精品成人综合77777| 久久热精品热| 亚洲精品自拍成人| 成人黄色视频免费在线看| 国产精品无大码| 亚洲精品日韩在线中文字幕| 国产女主播在线喷水免费视频网站| 在线观看av片永久免费下载| 天天躁夜夜躁狠狠久久av| 永久网站在线| 大陆偷拍与自拍| h日本视频在线播放| 成人亚洲欧美一区二区av| 日韩熟女老妇一区二区性免费视频| 插逼视频在线观看| 全区人妻精品视频| 久久人妻熟女aⅴ| 欧美最新免费一区二区三区| 亚洲欧美日韩东京热| 精品一品国产午夜福利视频| 国产午夜精品一二区理论片| 蜜桃在线观看..| 亚洲av二区三区四区| 成人亚洲欧美一区二区av| 成人毛片60女人毛片免费| 日韩av不卡免费在线播放| 亚洲精品久久久久久婷婷小说| 亚洲精品久久久久久婷婷小说| 汤姆久久久久久久影院中文字幕| 国产 一区精品| 日本wwww免费看| 午夜激情久久久久久久| 色视频在线一区二区三区| 人人妻人人添人人爽欧美一区卜| 日韩强制内射视频| 伦理电影免费视频| 亚洲av电影在线观看一区二区三区| √禁漫天堂资源中文www| 精品久久久久久久久av| 午夜免费鲁丝| 欧美 亚洲 国产 日韩一| 91精品国产九色| 亚洲人成网站在线播| 欧美97在线视频| 日韩成人伦理影院| 国产av码专区亚洲av| 最新的欧美精品一区二区| 亚洲欧洲国产日韩| 久久人人爽人人爽人人片va| 欧美少妇被猛烈插入视频| 久久国产亚洲av麻豆专区| 老司机亚洲免费影院| 国产精品久久久久久久久免| 中文字幕亚洲精品专区| 日本91视频免费播放| 18+在线观看网站| 成人亚洲精品一区在线观看| 成人亚洲欧美一区二区av| 在线观看av片永久免费下载| 一级,二级,三级黄色视频| 日产精品乱码卡一卡2卡三| 一区在线观看完整版| 日本午夜av视频| 国产一区二区在线观看日韩| 女性生殖器流出的白浆| 日韩一区二区视频免费看| 91成人精品电影| 一级爰片在线观看| 免费人成在线观看视频色| 精品一区二区三卡| 少妇精品久久久久久久| 十八禁高潮呻吟视频 | 人人妻人人澡人人爽人人夜夜| 美女国产视频在线观看| 国产成人精品久久久久久| 国产熟女欧美一区二区| 人体艺术视频欧美日本| 少妇高潮的动态图| 国产探花极品一区二区| 色5月婷婷丁香| 国产成人午夜福利电影在线观看| 少妇被粗大的猛进出69影院 | 欧美 亚洲 国产 日韩一| 大片免费播放器 马上看| 嫩草影院入口| 亚洲av免费高清在线观看| 美女内射精品一级片tv| 国产综合精华液| 免费大片18禁| 国产成人aa在线观看| 免费看av在线观看网站| 国产精品麻豆人妻色哟哟久久| 日日摸夜夜添夜夜添av毛片| 国产成人精品福利久久| 永久免费av网站大全| 日产精品乱码卡一卡2卡三| 国内揄拍国产精品人妻在线| 日韩av免费高清视频| 男男h啪啪无遮挡| 国产无遮挡羞羞视频在线观看| av视频免费观看在线观看| 国产精品一区二区三区四区免费观看| 99精国产麻豆久久婷婷| 精品人妻偷拍中文字幕| 少妇人妻 视频| 18禁裸乳无遮挡动漫免费视频| 多毛熟女@视频| 亚洲精品久久午夜乱码| 能在线免费看毛片的网站| 久久久久精品性色| 在线观看www视频免费| 黑人巨大精品欧美一区二区蜜桃 | a级毛片在线看网站| 成人国产av品久久久| 一区二区三区精品91| 蜜臀久久99精品久久宅男| 成人国产麻豆网| 成人国产麻豆网| 国产综合精华液| 午夜福利影视在线免费观看| 精品亚洲乱码少妇综合久久| 亚洲电影在线观看av| 亚洲精品色激情综合| 如日韩欧美国产精品一区二区三区 | 伦理电影免费视频| 国产精品久久久久成人av| 九九爱精品视频在线观看| 纯流量卡能插随身wifi吗| 五月开心婷婷网| 精品国产乱码久久久久久小说| 久久久久久久精品精品| 国产精品伦人一区二区| 十八禁网站网址无遮挡 | 天堂俺去俺来也www色官网| 亚洲精品乱码久久久v下载方式| 亚洲精品第二区| 国产精品麻豆人妻色哟哟久久| 国产熟女欧美一区二区| 中国三级夫妇交换| 美女cb高潮喷水在线观看| 亚洲精品456在线播放app| 国产成人一区二区在线| 国产午夜精品久久久久久一区二区三区| 中国国产av一级| a级毛色黄片| 久久久久久久亚洲中文字幕| 寂寞人妻少妇视频99o| 国国产精品蜜臀av免费| 久热久热在线精品观看| 高清午夜精品一区二区三区| 在线观看三级黄色| 精品一区二区三区视频在线| 国产精品人妻久久久久久| 亚洲精品,欧美精品| 亚洲高清免费不卡视频| videos熟女内射| 精品人妻偷拍中文字幕| 国产精品三级大全| 一级av片app| 国产成人一区二区在线| 寂寞人妻少妇视频99o| 国产成人精品婷婷| 国产日韩欧美在线精品| 欧美激情国产日韩精品一区| 日本色播在线视频| 王馨瑶露胸无遮挡在线观看| 香蕉精品网在线| 成人美女网站在线观看视频| 亚洲经典国产精华液单| 免费高清在线观看视频在线观看| 日韩成人伦理影院| 永久网站在线| 我的女老师完整版在线观看| 在线观看免费高清a一片| 国产一级毛片在线| 少妇被粗大的猛进出69影院 | 日产精品乱码卡一卡2卡三| 少妇的逼好多水| 2021少妇久久久久久久久久久| 久久久久久久亚洲中文字幕| 欧美日韩综合久久久久久| 免费观看a级毛片全部| 一级a做视频免费观看| 好男人视频免费观看在线| 美女视频免费永久观看网站| 色5月婷婷丁香| 亚洲av中文av极速乱| 欧美性感艳星| 亚洲av综合色区一区| 午夜福利,免费看| 最新的欧美精品一区二区| 欧美国产精品一级二级三级 | 插逼视频在线观看| 深夜a级毛片| 国产高清国产精品国产三级| 性色av一级| 亚洲,欧美,日韩| 婷婷色综合大香蕉| 色视频在线一区二区三区| 久久99热6这里只有精品| 亚洲人成网站在线观看播放| 天天操日日干夜夜撸| 亚洲精品乱久久久久久| 欧美日韩国产mv在线观看视频| 国产深夜福利视频在线观看| 国产探花极品一区二区| 成人黄色视频免费在线看| 最近中文字幕高清免费大全6| 国产视频内射| 人妻一区二区av| 国产黄频视频在线观看| 中国国产av一级| 性色av一级| .国产精品久久| 草草在线视频免费看| 天美传媒精品一区二区| 性色av一级| 色5月婷婷丁香| 亚洲真实伦在线观看| 久久久久人妻精品一区果冻| 99热这里只有是精品在线观看| 欧美高清成人免费视频www| 亚洲精品久久久久久婷婷小说| 啦啦啦在线观看免费高清www| 亚洲四区av| 美女中出高潮动态图| 另类亚洲欧美激情| 国产黄色免费在线视频| 免费观看无遮挡的男女| 国产深夜福利视频在线观看| 欧美区成人在线视频| 99热网站在线观看| 丰满少妇做爰视频| 日日啪夜夜爽| 国产有黄有色有爽视频| 国产欧美另类精品又又久久亚洲欧美| 大片免费播放器 马上看| 午夜精品国产一区二区电影| 色5月婷婷丁香| 日日啪夜夜撸| 国产精品一区二区在线不卡| 青春草国产在线视频| 国产伦理片在线播放av一区| 国产一区二区三区av在线| 我要看日韩黄色一级片| 国产视频内射| 内射极品少妇av片p| 亚洲美女黄色视频免费看| 日韩成人av中文字幕在线观看| 在线免费观看不下载黄p国产| 国产伦理片在线播放av一区| 精品国产一区二区久久| 免费不卡的大黄色大毛片视频在线观看| 九九久久精品国产亚洲av麻豆| av卡一久久| 九九爱精品视频在线观看| 中文字幕久久专区| 免费少妇av软件| 2022亚洲国产成人精品| 一级毛片久久久久久久久女| av线在线观看网站| 人妻夜夜爽99麻豆av| 午夜免费鲁丝| 91成人精品电影| 这个男人来自地球电影免费观看 | 高清视频免费观看一区二区| 免费播放大片免费观看视频在线观看| 女性被躁到高潮视频| 精品亚洲成a人片在线观看| 伊人亚洲综合成人网| 人人澡人人妻人| 高清毛片免费看| 欧美精品人与动牲交sv欧美| 人人澡人人妻人| 成人二区视频| av天堂中文字幕网| 精品一区在线观看国产| 日本黄色片子视频| 日本黄色日本黄色录像| 免费不卡的大黄色大毛片视频在线观看| 午夜91福利影院| 国产淫语在线视频| 午夜日本视频在线| 免费观看性生交大片5| 观看免费一级毛片| 成年美女黄网站色视频大全免费 | 免费久久久久久久精品成人欧美视频 | 久久97久久精品| 国产在线一区二区三区精| 纵有疾风起免费观看全集完整版| 国产亚洲最大av| 婷婷色av中文字幕| 肉色欧美久久久久久久蜜桃| 激情五月婷婷亚洲| 免费看日本二区| 免费av不卡在线播放| 大码成人一级视频| 欧美丝袜亚洲另类| 99精国产麻豆久久婷婷| 国产精品免费大片| 超碰97精品在线观看| 欧美日韩综合久久久久久| 久久人人爽人人爽人人片va| 免费黄色在线免费观看| .国产精品久久| 一级二级三级毛片免费看| 亚洲精品乱码久久久v下载方式| 日韩人妻高清精品专区| 国产av国产精品国产| 在线观看一区二区三区激情| 精品国产一区二区三区久久久樱花| av又黄又爽大尺度在线免费看| 国产视频内射| 精品国产露脸久久av麻豆| 搡老乐熟女国产| 久久精品久久久久久久性| 高清午夜精品一区二区三区| 免费久久久久久久精品成人欧美视频 | 国内精品宾馆在线| 欧美激情国产日韩精品一区| 精品一区二区免费观看| 国产精品久久久久成人av| 伦精品一区二区三区| 国产熟女午夜一区二区三区 | 大香蕉97超碰在线| 欧美+日韩+精品| 美女国产视频在线观看| 美女cb高潮喷水在线观看| 97超碰精品成人国产| 国产伦精品一区二区三区视频9| 午夜激情福利司机影院| 欧美最新免费一区二区三区| 大片电影免费在线观看免费| 纯流量卡能插随身wifi吗| 精品熟女少妇av免费看| 夫妻午夜视频| 九九爱精品视频在线观看| 男的添女的下面高潮视频| 久久 成人 亚洲| 久久国产精品男人的天堂亚洲 | 老司机影院成人| 久久精品国产亚洲av涩爱| 桃花免费在线播放| 亚洲人与动物交配视频| 少妇高潮的动态图| 欧美精品国产亚洲| 欧美 日韩 精品 国产| 国产女主播在线喷水免费视频网站| 久久久久久久精品精品| 日韩精品免费视频一区二区三区 | 国产在视频线精品| 国产成人一区二区在线| 女人精品久久久久毛片| 久久精品夜色国产| 亚洲成色77777| 国产视频首页在线观看| 自线自在国产av| 男女啪啪激烈高潮av片| 国产免费一区二区三区四区乱码| 国产探花极品一区二区| 久久人人爽人人爽人人片va| 久久久久久人妻| 两个人的视频大全免费| videossex国产| 国产精品一区二区三区四区免费观看| 久久人人爽av亚洲精品天堂| 精品亚洲乱码少妇综合久久| 国产精品国产三级专区第一集| 久久韩国三级中文字幕| 免费人成在线观看视频色| 深夜a级毛片| 老司机亚洲免费影院| 国产成人freesex在线| 亚洲精品色激情综合| 国产精品国产三级国产专区5o| 色网站视频免费| 噜噜噜噜噜久久久久久91| 97在线人人人人妻| 免费高清在线观看视频在线观看| videossex国产| 伦理电影大哥的女人| 9色porny在线观看| 又爽又黄a免费视频| 日韩中字成人| 麻豆精品久久久久久蜜桃| 亚洲欧美日韩东京热| av在线观看视频网站免费| 日韩欧美精品免费久久| 国产精品人妻久久久久久| 日日撸夜夜添| 亚洲av免费高清在线观看| 老司机影院毛片| 欧美另类一区| av在线app专区| 亚洲美女黄色视频免费看| 免费观看性生交大片5| 亚洲第一av免费看| 五月伊人婷婷丁香| 在现免费观看毛片| 日韩,欧美,国产一区二区三区| 日韩三级伦理在线观看| 亚洲av综合色区一区| 三级经典国产精品| 免费人成在线观看视频色| 亚洲天堂av无毛| 老司机影院成人| 黄色配什么色好看| 18禁在线播放成人免费| 十八禁高潮呻吟视频 | 欧美日韩国产mv在线观看视频| 青春草视频在线免费观看| 久久久亚洲精品成人影院| av女优亚洲男人天堂| h日本视频在线播放| 国产精品国产三级专区第一集| 国产精品久久久久久av不卡| 国产精品国产av在线观看| 日韩精品有码人妻一区| 亚洲精品第二区| 久久99一区二区三区| 国产精品久久久久久av不卡| 99re6热这里在线精品视频| √禁漫天堂资源中文www| 精品亚洲成国产av| 国产精品秋霞免费鲁丝片| 成人二区视频| 日韩av不卡免费在线播放| 插阴视频在线观看视频| 精品少妇黑人巨大在线播放| xxx大片免费视频| 伊人久久国产一区二区| 中文欧美无线码| 国产成人aa在线观看| 国产精品久久久久久久电影| 久久韩国三级中文字幕| 伊人久久精品亚洲午夜| 日韩强制内射视频| 欧美97在线视频| 亚洲一级一片aⅴ在线观看| 九九爱精品视频在线观看| 日本av免费视频播放| 韩国高清视频一区二区三区| 国产极品天堂在线| 亚洲精品成人av观看孕妇| 国产精品一区二区性色av| 嫩草影院入口| 欧美+日韩+精品| 激情五月婷婷亚洲| 日本欧美视频一区| 亚洲av男天堂| 久久久久久人妻| av免费观看日本| 国产精品一二三区在线看| 一级毛片黄色毛片免费观看视频| 国产成人a∨麻豆精品| 国产探花极品一区二区| 久久久久人妻精品一区果冻| 男人添女人高潮全过程视频| 日本午夜av视频| 另类精品久久| 如日韩欧美国产精品一区二区三区 | 久久精品国产a三级三级三级| 大话2 男鬼变身卡| 久久6这里有精品| 在线播放无遮挡| 一区二区av电影网| 午夜激情福利司机影院| 日本黄大片高清| 日韩免费高清中文字幕av| 国产黄片美女视频| 亚洲欧美日韩卡通动漫| 欧美+日韩+精品| 久久99热这里只频精品6学生| 99re6热这里在线精品视频| 9色porny在线观看| a级片在线免费高清观看视频| 欧美少妇被猛烈插入视频| 男人添女人高潮全过程视频| 丰满人妻一区二区三区视频av| 亚洲精品乱码久久久v下载方式| 日本午夜av视频| 国产亚洲一区二区精品| 日日啪夜夜撸| 亚洲人与动物交配视频| 纵有疾风起免费观看全集完整版| 下体分泌物呈黄色| 纵有疾风起免费观看全集完整版| 日本欧美国产在线视频| 你懂的网址亚洲精品在线观看| 亚洲高清免费不卡视频| 伦理电影大哥的女人| 国产成人a∨麻豆精品| 色网站视频免费| 亚洲综合精品二区| 久久99一区二区三区| 日韩亚洲欧美综合| 最近中文字幕2019免费版| 18禁在线无遮挡免费观看视频| 精品酒店卫生间| 一二三四中文在线观看免费高清| 久久精品国产亚洲av涩爱| 久久久久人妻精品一区果冻| 午夜免费观看性视频| 欧美另类一区| 亚洲欧美成人精品一区二区| 美女xxoo啪啪120秒动态图| 中文字幕久久专区| 2022亚洲国产成人精品| 久久久久国产精品人妻一区二区| av网站免费在线观看视频| 日韩亚洲欧美综合| 精品熟女少妇av免费看| 亚洲精品国产av成人精品| av天堂中文字幕网| videossex国产| 少妇高潮的动态图| 高清不卡的av网站| 久久鲁丝午夜福利片| 久久97久久精品| 国产高清不卡午夜福利| 91精品伊人久久大香线蕉| 波野结衣二区三区在线| 日韩欧美精品免费久久| 免费观看av网站的网址| 日韩亚洲欧美综合| 永久免费av网站大全| a 毛片基地| 91午夜精品亚洲一区二区三区| 99热6这里只有精品| 91aial.com中文字幕在线观看| 欧美3d第一页| 国产极品粉嫩免费观看在线 | 亚洲欧洲日产国产| 久久毛片免费看一区二区三区| 日本av手机在线免费观看| 亚洲成色77777| 夜夜骑夜夜射夜夜干| 中文天堂在线官网| 国产视频内射| 亚洲精品日韩av片在线观看| 欧美日韩亚洲高清精品| 又大又黄又爽视频免费| 亚洲精品一区蜜桃| 一级毛片 在线播放| tube8黄色片| 国产亚洲午夜精品一区二区久久| 国产精品人妻久久久久久| 欧美性感艳星| www.av在线官网国产| 六月丁香七月| av.在线天堂| 人妻少妇偷人精品九色| 男人舔奶头视频| 亚洲av综合色区一区| 日韩中字成人| 精品亚洲乱码少妇综合久久| 国产一区亚洲一区在线观看| 人人妻人人看人人澡| 一二三四中文在线观看免费高清|