• <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
    大型風力機整機氣動彈性響應計算
    大香蕉97超碰在线| 一二三四在线观看免费中文在 | 黄网站色视频无遮挡免费观看| 午夜av观看不卡| 免费高清在线观看日韩| 亚洲精华国产精华液的使用体验| 久久精品久久久久久久性| 卡戴珊不雅视频在线播放| 日韩中文字幕视频在线看片| 国产 一区精品| 亚洲av欧美aⅴ国产| 亚洲久久久国产精品| 欧美97在线视频| 日韩 亚洲 欧美在线| 亚洲内射少妇av| 大香蕉久久成人网| 99国产综合亚洲精品| 精品人妻熟女毛片av久久网站| 欧美亚洲日本最大视频资源| 咕卡用的链子| 日日爽夜夜爽网站| 黄色怎么调成土黄色| √禁漫天堂资源中文www| 菩萨蛮人人尽说江南好唐韦庄| 国产精品国产三级专区第一集| 国产女主播在线喷水免费视频网站| 国产深夜福利视频在线观看| 欧美3d第一页| 久久97久久精品| 成人手机av| 国产白丝娇喘喷水9色精品| 亚洲欧美色中文字幕在线| 卡戴珊不雅视频在线播放| 少妇人妻 视频| 日本与韩国留学比较| 精品人妻在线不人妻| 亚洲av免费高清在线观看| 日本av手机在线免费观看| 亚洲av成人精品一二三区| 成人免费观看视频高清| 另类亚洲欧美激情| 蜜桃国产av成人99| 制服人妻中文乱码| 成人毛片a级毛片在线播放| 国产av一区二区精品久久| 99久久人妻综合| 久久国产亚洲av麻豆专区| 交换朋友夫妻互换小说| 成年动漫av网址| 蜜桃国产av成人99| 嫩草影院入口| 亚洲欧美成人综合另类久久久| 成人免费观看视频高清| 国产欧美亚洲国产| 亚洲精品中文字幕在线视频| 久久久久人妻精品一区果冻| 汤姆久久久久久久影院中文字幕| 欧美xxⅹ黑人| 巨乳人妻的诱惑在线观看| www日本在线高清视频| 大片电影免费在线观看免费| 欧美日韩视频精品一区| 伦理电影大哥的女人| 丝袜美足系列| 日韩大片免费观看网站| 这个男人来自地球电影免费观看 | 18禁裸乳无遮挡动漫免费视频| 人人妻人人澡人人爽人人夜夜| 丝袜喷水一区| 久久人人97超碰香蕉20202| 丝袜美足系列| 日韩中文字幕视频在线看片| 久久99蜜桃精品久久| 中文字幕另类日韩欧美亚洲嫩草| 极品少妇高潮喷水抽搐| 亚洲四区av| 在线观看美女被高潮喷水网站| 色吧在线观看| av一本久久久久| 久久影院123| 91国产中文字幕| 欧美激情国产日韩精品一区| 精品亚洲成a人片在线观看| 蜜桃国产av成人99| 啦啦啦视频在线资源免费观看| 日韩大片免费观看网站| 欧美bdsm另类| videosex国产| 又黄又爽又刺激的免费视频.| 人妻少妇偷人精品九色| 黄色毛片三级朝国网站| 草草在线视频免费看| 免费观看在线日韩| 国产在线视频一区二区| 夜夜爽夜夜爽视频| 国产精品99久久99久久久不卡 | 91成人精品电影| 蜜桃国产av成人99| 大片免费播放器 马上看| 精品久久蜜臀av无| 久久久久精品久久久久真实原创| 中文天堂在线官网| 欧美激情 高清一区二区三区| 伊人久久国产一区二区| 亚洲国产欧美日韩在线播放| 国产免费视频播放在线视频| 各种免费的搞黄视频| 搡女人真爽免费视频火全软件| 免费女性裸体啪啪无遮挡网站| 久久久精品免费免费高清| 欧美成人精品欧美一级黄| 欧美国产精品va在线观看不卡| 一级毛片 在线播放| 亚洲精品aⅴ在线观看| 五月开心婷婷网| 蜜桃在线观看..| 成人二区视频| 精品卡一卡二卡四卡免费| av在线播放精品| 国产麻豆69| 边亲边吃奶的免费视频| 亚洲精品视频女| 国产一区二区激情短视频 | 国产片内射在线| 日本午夜av视频| 寂寞人妻少妇视频99o| 日本-黄色视频高清免费观看| 免费看不卡的av| 久久精品国产亚洲av天美| 亚洲国产成人一精品久久久| 又黄又粗又硬又大视频| 久久热在线av| 中文字幕av电影在线播放| 少妇的丰满在线观看| 视频在线观看一区二区三区| 天天操日日干夜夜撸| 精品少妇黑人巨大在线播放| 国产一区有黄有色的免费视频| 在线观看一区二区三区激情| 一级黄片播放器| 在线观看免费高清a一片| 亚洲精品av麻豆狂野| 国产亚洲午夜精品一区二区久久| 免费观看a级毛片全部| 国产成人欧美| 天天躁夜夜躁狠狠久久av| 深夜精品福利| 日韩,欧美,国产一区二区三区| 日韩成人av中文字幕在线观看| 观看av在线不卡| 另类亚洲欧美激情| 精品人妻一区二区三区麻豆| 免费黄频网站在线观看国产| 一级黄片播放器| 高清欧美精品videossex| 亚洲四区av| 日本av手机在线免费观看| 欧美精品国产亚洲| 午夜激情久久久久久久| 欧美人与性动交α欧美精品济南到 | 国产精品不卡视频一区二区| 久久精品熟女亚洲av麻豆精品| 免费少妇av软件| 女的被弄到高潮叫床怎么办| 久久久久国产精品人妻一区二区| 一级黄片播放器| 校园人妻丝袜中文字幕| 国产无遮挡羞羞视频在线观看| 91国产中文字幕| 人成视频在线观看免费观看| 婷婷色麻豆天堂久久| 国产精品一区二区在线观看99| kizo精华| 欧美激情 高清一区二区三区| 精品亚洲成国产av| 九九爱精品视频在线观看| 久久人人97超碰香蕉20202| 久久久久国产精品人妻一区二区| 麻豆精品久久久久久蜜桃| 99久久精品国产国产毛片| 黄色毛片三级朝国网站| 国内精品宾馆在线| 久久 成人 亚洲| 久久综合国产亚洲精品| 啦啦啦中文免费视频观看日本| 丰满乱子伦码专区| 啦啦啦在线观看免费高清www| 成人毛片60女人毛片免费| 国产熟女欧美一区二区| 男女啪啪激烈高潮av片| 成人亚洲欧美一区二区av| 边亲边吃奶的免费视频| 欧美最新免费一区二区三区| www.av在线官网国产| 久久国产精品大桥未久av| 91精品三级在线观看| 高清视频免费观看一区二区| 久久亚洲国产成人精品v| 亚洲国产看品久久| 九草在线视频观看| 国产伦理片在线播放av一区| 国产免费又黄又爽又色| 日韩欧美精品免费久久| 精品国产露脸久久av麻豆| 国产精品国产三级专区第一集| 久久99一区二区三区| 国产精品一国产av| 日韩一本色道免费dvd| 母亲3免费完整高清在线观看 | 热re99久久精品国产66热6| 少妇 在线观看| 亚洲精品456在线播放app| 99热全是精品| 亚洲av综合色区一区| 狂野欧美激情性bbbbbb| 美女国产高潮福利片在线看| 国产69精品久久久久777片| 日韩 亚洲 欧美在线| a级片在线免费高清观看视频| 日本-黄色视频高清免费观看| 黄色 视频免费看| 国产国拍精品亚洲av在线观看| 午夜福利网站1000一区二区三区| 精品第一国产精品| 国产av精品麻豆| tube8黄色片| 女性被躁到高潮视频| 精品久久久久久电影网| 在线天堂最新版资源| 国产淫语在线视频| 国产日韩欧美视频二区| 日本av手机在线免费观看| 制服诱惑二区| 免费大片黄手机在线观看| 韩国精品一区二区三区 | 大香蕉久久成人网| 国产免费又黄又爽又色| 久久久久国产精品人妻一区二区| 9色porny在线观看| 亚洲天堂av无毛| 夫妻性生交免费视频一级片| 天天影视国产精品| 国产精品99久久99久久久不卡 | 免费看av在线观看网站| 久久国产亚洲av麻豆专区| 精品卡一卡二卡四卡免费| 97精品久久久久久久久久精品| 男女无遮挡免费网站观看| 深夜精品福利| 一级爰片在线观看| 亚洲成人一二三区av| 一本久久精品| 亚洲欧洲精品一区二区精品久久久 | 少妇猛男粗大的猛烈进出视频| 一本—道久久a久久精品蜜桃钙片| 久久久久久久久久久久大奶| 国产在线视频一区二区| 丝袜喷水一区| 三级国产精品片| 咕卡用的链子| 搡老乐熟女国产| 极品人妻少妇av视频| 边亲边吃奶的免费视频| 不卡视频在线观看欧美| 国产亚洲精品第一综合不卡 | 在线看a的网站| 国产激情久久老熟女| 精品久久国产蜜桃| 精品国产乱码久久久久久小说| 制服丝袜香蕉在线| 久久久久精品人妻al黑| 日本欧美视频一区| 久热久热在线精品观看| 欧美人与性动交α欧美软件 | 成年av动漫网址| 亚洲欧美日韩另类电影网站| 男男h啪啪无遮挡| 丰满迷人的少妇在线观看| 亚洲精品乱久久久久久| 久久精品国产亚洲av涩爱| 在线观看人妻少妇| 免费黄色在线免费观看| 色吧在线观看| 国产精品久久久久久久电影| 少妇精品久久久久久久| 国产成人欧美| 熟女人妻精品中文字幕| 18禁在线无遮挡免费观看视频| 亚洲精品国产色婷婷电影| 精品久久久久久电影网| 中文欧美无线码| 高清黄色对白视频在线免费看| 黄色配什么色好看| 欧美老熟妇乱子伦牲交| 免费av不卡在线播放| 看非洲黑人一级黄片| 欧美国产精品一级二级三级| 卡戴珊不雅视频在线播放| 曰老女人黄片| 久久久久久久国产电影| 午夜福利乱码中文字幕| 亚洲av电影在线进入| 少妇高潮的动态图| 久久鲁丝午夜福利片| 久久av网站| 免费看不卡的av| av视频免费观看在线观看| 国产色婷婷99| 国产精品久久久av美女十八| 日日啪夜夜爽| 欧美成人午夜免费资源| 香蕉精品网在线| 热re99久久精品国产66热6| 97精品久久久久久久久久精品| a级毛片黄视频| 如何舔出高潮| 91精品三级在线观看| 最近最新中文字幕大全免费视频 | 在线观看免费日韩欧美大片| 日韩中文字幕视频在线看片| 成人无遮挡网站| 美女内射精品一级片tv| 狠狠婷婷综合久久久久久88av| 成人毛片60女人毛片免费| 亚洲欧美精品自产自拍| 狠狠婷婷综合久久久久久88av| 久久人妻熟女aⅴ| 欧美xxⅹ黑人| 天堂中文最新版在线下载| 你懂的网址亚洲精品在线观看| 99国产综合亚洲精品| 人体艺术视频欧美日本| 99国产综合亚洲精品| 你懂的网址亚洲精品在线观看| 欧美日韩一区二区视频在线观看视频在线| 夫妻性生交免费视频一级片| 插逼视频在线观看| 国产精品久久久久久久电影| 亚洲第一区二区三区不卡| 久久99一区二区三区| 22中文网久久字幕| 国产福利在线免费观看视频| 午夜福利网站1000一区二区三区| 久久久久久久久久久久大奶| 亚洲激情五月婷婷啪啪| 婷婷色麻豆天堂久久| 最新的欧美精品一区二区| 这个男人来自地球电影免费观看 | 国产精品女同一区二区软件| 伦理电影免费视频| 黑人欧美特级aaaaaa片| 日本av手机在线免费观看| 久久国产亚洲av麻豆专区| 成人无遮挡网站| 免费日韩欧美在线观看| 又黄又粗又硬又大视频| 水蜜桃什么品种好| 狂野欧美激情性bbbbbb| 啦啦啦啦在线视频资源| 丰满饥渴人妻一区二区三| 中文字幕亚洲精品专区| 亚洲成人一二三区av| 国产深夜福利视频在线观看| 亚洲第一av免费看| 亚洲一码二码三码区别大吗| 成人漫画全彩无遮挡| 自线自在国产av| freevideosex欧美| 亚洲精品视频女| 午夜免费鲁丝| 天天影视国产精品| 婷婷色综合大香蕉| 亚洲精品国产色婷婷电影| 免费观看性生交大片5| 亚洲婷婷狠狠爱综合网| 99久国产av精品国产电影| 久久人人爽人人片av| 亚洲精品国产av蜜桃| 国产亚洲精品第一综合不卡 | 日本色播在线视频| 日本vs欧美在线观看视频| 国产成人精品婷婷| 中文字幕最新亚洲高清| 卡戴珊不雅视频在线播放| 五月天丁香电影| 韩国高清视频一区二区三区| 国产爽快片一区二区三区| 国产亚洲av片在线观看秒播厂| 午夜福利乱码中文字幕| 免费大片黄手机在线观看| 国产成人精品久久久久久| 精品久久蜜臀av无| 三上悠亚av全集在线观看| 精品人妻偷拍中文字幕| 大片免费播放器 马上看| 人人妻人人爽人人添夜夜欢视频| 国产一区二区在线观看av| 色视频在线一区二区三区| 高清欧美精品videossex| 久久久a久久爽久久v久久| 欧美bdsm另类| 在线观看免费高清a一片| 老司机影院成人| 国产69精品久久久久777片| 精品亚洲成国产av| 国产精品不卡视频一区二区| 最近中文字幕2019免费版| 日韩av在线免费看完整版不卡| 日韩不卡一区二区三区视频在线| 多毛熟女@视频| 国产深夜福利视频在线观看| 亚洲精品,欧美精品| 精品一区二区三区视频在线| 夜夜爽夜夜爽视频| 狠狠精品人妻久久久久久综合| 亚洲成人av在线免费| 哪个播放器可以免费观看大片| 亚洲精品国产av成人精品| 国产在线视频一区二区| 精品第一国产精品| 国产av一区二区精品久久| av网站免费在线观看视频| videossex国产| 精品人妻熟女毛片av久久网站| 丰满乱子伦码专区| 高清毛片免费看| 免费高清在线观看视频在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲国产精品999| 午夜激情久久久久久久| 婷婷色av中文字幕| 亚洲欧美一区二区三区黑人 | 日本午夜av视频| 丰满少妇做爰视频| 国产成人午夜福利电影在线观看| 国产精品成人在线| 午夜福利视频精品| 日韩精品免费视频一区二区三区 | 国产在线一区二区三区精| 满18在线观看网站| 一边亲一边摸免费视频| 中文字幕亚洲精品专区| 制服诱惑二区| 欧美精品亚洲一区二区| 大陆偷拍与自拍| 黄色毛片三级朝国网站| 色5月婷婷丁香| 国产男女内射视频| 久久精品国产亚洲av天美| 欧美日韩成人在线一区二区| 高清欧美精品videossex| 国产精品.久久久| 国产欧美日韩综合在线一区二区| 中文字幕人妻丝袜制服| 久久精品国产自在天天线| 久久久久人妻精品一区果冻| a 毛片基地| 成人综合一区亚洲| 欧美bdsm另类| 男女国产视频网站| 亚洲精品日本国产第一区| 一个人免费看片子| 99国产综合亚洲精品| 草草在线视频免费看| 免费观看无遮挡的男女| 国产成人欧美| 国产成人一区二区在线| 1024视频免费在线观看| 欧美少妇被猛烈插入视频| 国产精品一区www在线观看| 精品国产一区二区三区四区第35| 香蕉丝袜av| 日日撸夜夜添| 国产免费视频播放在线视频| 999精品在线视频| 国产一区二区三区综合在线观看 | 男女高潮啪啪啪动态图| 欧美xxxx性猛交bbbb| h视频一区二区三区| 久久久亚洲精品成人影院| 欧美国产精品一级二级三级| 国产成人aa在线观看| h视频一区二区三区| 国产永久视频网站| 亚洲精品国产色婷婷电影| 人妻系列 视频| 欧美老熟妇乱子伦牲交| av播播在线观看一区| 日本av手机在线免费观看| 亚洲精华国产精华液的使用体验| 国产亚洲欧美精品永久| 精品国产一区二区三区四区第35| 久久久久精品久久久久真实原创| 狂野欧美激情性bbbbbb| 免费av不卡在线播放| 久久影院123| 欧美精品高潮呻吟av久久| 亚洲色图 男人天堂 中文字幕 | 精品视频人人做人人爽| 两性夫妻黄色片 | 丝袜脚勾引网站| 考比视频在线观看| 欧美日韩视频精品一区| 国产免费视频播放在线视频| 国产免费现黄频在线看| 午夜老司机福利剧场| 午夜免费观看性视频| 最新的欧美精品一区二区| 校园人妻丝袜中文字幕| 99久久人妻综合| 日本免费在线观看一区| 亚洲成国产人片在线观看| 老司机亚洲免费影院| 亚洲成人一二三区av| 欧美97在线视频| 国产一区二区三区综合在线观看 | 亚洲一级一片aⅴ在线观看| 在线观看美女被高潮喷水网站| 午夜老司机福利剧场| 日韩在线高清观看一区二区三区| 亚洲欧美成人综合另类久久久| 在线观看免费高清a一片| 蜜桃在线观看..| 亚洲成国产人片在线观看| 精品99又大又爽又粗少妇毛片| 亚洲精品自拍成人| 国产1区2区3区精品| 国产又爽黄色视频| 精品福利永久在线观看| 青春草亚洲视频在线观看| 亚洲成人一二三区av| 国产成人精品福利久久| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久精品久久久| 在线观看一区二区三区激情| 日本黄大片高清| 男女高潮啪啪啪动态图| 亚洲国产成人一精品久久久| 国产精品一区二区在线观看99| 黄片无遮挡物在线观看| 男女国产视频网站| 成年人午夜在线观看视频| 亚洲美女搞黄在线观看| 亚洲精品美女久久久久99蜜臀 | av在线播放精品| 一个人免费看片子| 超色免费av| 人妻人人澡人人爽人人| 国产成人午夜福利电影在线观看| 国产极品粉嫩免费观看在线| 亚洲av免费高清在线观看| 欧美人与善性xxx| 七月丁香在线播放| 三级国产精品片| 插逼视频在线观看| 久久久久久久久久成人| 亚洲精品乱久久久久久| 国产无遮挡羞羞视频在线观看| 亚洲成色77777| 久久ye,这里只有精品| 韩国精品一区二区三区 | 久久韩国三级中文字幕| 日韩视频在线欧美| 少妇被粗大猛烈的视频| 黄色视频在线播放观看不卡| 亚洲精品美女久久av网站| 免费不卡的大黄色大毛片视频在线观看| 国产精品国产三级国产专区5o| 精品国产一区二区三区久久久樱花| 91久久精品国产一区二区三区| 午夜福利乱码中文字幕| 黑人欧美特级aaaaaa片| 国产一区二区激情短视频 | 精品午夜福利在线看| 中文字幕免费在线视频6| 久久免费观看电影| 人妻人人澡人人爽人人| 久久久久久久久久成人| 大片电影免费在线观看免费| 中文字幕免费在线视频6| 视频在线观看一区二区三区| 亚洲综合精品二区| 欧美日本中文国产一区发布| 多毛熟女@视频| videos熟女内射| 免费观看a级毛片全部| 亚洲国产色片| 亚洲成色77777| 欧美日韩视频高清一区二区三区二| 国产色婷婷99| 国产成人91sexporn| 婷婷成人精品国产| av网站免费在线观看视频| 人妻人人澡人人爽人人| 岛国毛片在线播放| 国产视频首页在线观看| 天天影视国产精品| 菩萨蛮人人尽说江南好唐韦庄| 欧美3d第一页| 午夜日本视频在线| 在线 av 中文字幕| 人体艺术视频欧美日本| 啦啦啦视频在线资源免费观看| 亚洲欧美成人精品一区二区| 2018国产大陆天天弄谢| 97在线人人人人妻| 日韩伦理黄色片| a级毛色黄片| 久久久亚洲精品成人影院| 又黄又粗又硬又大视频| 亚洲经典国产精华液单| 久久精品久久精品一区二区三区| 日韩在线高清观看一区二区三区| 国产国拍精品亚洲av在线观看| 日本91视频免费播放|