• <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
    大型風力機整機氣動彈性響應計算
    在线观看国产h片| 国产欧美日韩一区二区三区在线| 亚洲精品一区蜜桃| 国产一区二区激情短视频 | 性少妇av在线| 99国产精品99久久久久| 亚洲精品在线美女| 久久亚洲精品不卡| 国产精品av久久久久免费| 一区在线观看完整版| 久久久精品国产亚洲av高清涩受| 欧美日韩福利视频一区二区| 亚洲av成人精品一二三区| 国产一区二区三区av在线| 日韩大码丰满熟妇| 久久精品久久久久久噜噜老黄| 日本a在线网址| 亚洲欧美精品自产自拍| 交换朋友夫妻互换小说| av网站免费在线观看视频| 大话2 男鬼变身卡| 天天躁夜夜躁狠狠久久av| 国产女主播在线喷水免费视频网站| 亚洲成色77777| 黄色一级大片看看| 久久人妻熟女aⅴ| 免费女性裸体啪啪无遮挡网站| 日韩熟女老妇一区二区性免费视频| 99re6热这里在线精品视频| 99国产精品免费福利视频| 亚洲一区二区三区欧美精品| 欧美少妇被猛烈插入视频| 国产老妇伦熟女老妇高清| 无遮挡黄片免费观看| 黑人巨大精品欧美一区二区蜜桃| 午夜免费男女啪啪视频观看| 免费高清在线观看视频在线观看| 免费av中文字幕在线| 老司机影院成人| 男女午夜视频在线观看| 精品国产乱码久久久久久男人| 操出白浆在线播放| 99九九在线精品视频| 丝瓜视频免费看黄片| 国产欧美亚洲国产| 国产一区二区在线观看av| 脱女人内裤的视频| 精品亚洲乱码少妇综合久久| 日本a在线网址| 一本—道久久a久久精品蜜桃钙片| 黄色毛片三级朝国网站| 人人妻,人人澡人人爽秒播 | 国产精品一区二区免费欧美 | 精品一品国产午夜福利视频| 欧美精品人与动牲交sv欧美| 日本猛色少妇xxxxx猛交久久| 一区二区三区乱码不卡18| 99久久人妻综合| 看十八女毛片水多多多| 十分钟在线观看高清视频www| 国产av国产精品国产| 亚洲精品在线美女| 成人亚洲精品一区在线观看| 日本wwww免费看| 91精品国产国语对白视频| 黄色 视频免费看| 欧美黑人欧美精品刺激| 精品国产乱码久久久久久小说| 国产精品一区二区免费欧美 | 国产91精品成人一区二区三区 | 亚洲午夜精品一区,二区,三区| 久久午夜综合久久蜜桃| 亚洲国产欧美网| 最近手机中文字幕大全| 校园人妻丝袜中文字幕| 亚洲国产欧美在线一区| 午夜激情久久久久久久| 一区二区三区四区激情视频| 亚洲午夜精品一区,二区,三区| av国产久精品久网站免费入址| 在线观看国产h片| 日日夜夜操网爽| 99久久99久久久精品蜜桃| 最近手机中文字幕大全| 国产日韩欧美视频二区| 丁香六月天网| 深夜精品福利| 国产伦人伦偷精品视频| 制服诱惑二区| 啦啦啦在线观看免费高清www| 久久久久国产一级毛片高清牌| 亚洲成国产人片在线观看| 国产精品免费大片| 脱女人内裤的视频| 亚洲欧洲日产国产| 国产真人三级小视频在线观看| 狠狠精品人妻久久久久久综合| 999久久久国产精品视频| 国产又爽黄色视频| www.999成人在线观看| 欧美精品一区二区大全| 成人手机av| 最新在线观看一区二区三区 | 一级黄色大片毛片| av线在线观看网站| 午夜免费男女啪啪视频观看| 成在线人永久免费视频| 国产精品偷伦视频观看了| 国产亚洲欧美精品永久| 亚洲精品一卡2卡三卡4卡5卡 | av欧美777| 亚洲国产日韩一区二区| 亚洲精品美女久久久久99蜜臀 | 欧美另类一区| 国产成人a∨麻豆精品| 中文乱码字字幕精品一区二区三区| 欧美精品亚洲一区二区| 三上悠亚av全集在线观看| 在线观看免费视频网站a站| 欧美激情极品国产一区二区三区| 久久人妻熟女aⅴ| 中文字幕人妻熟女乱码| 欧美国产精品一级二级三级| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美一区二区三区国产| 精品少妇黑人巨大在线播放| 日本vs欧美在线观看视频| 欧美日韩国产mv在线观看视频| 久久精品国产a三级三级三级| 蜜桃在线观看..| 日韩精品免费视频一区二区三区| 中文精品一卡2卡3卡4更新| 18禁黄网站禁片午夜丰满| 亚洲国产av新网站| 国产高清videossex| 咕卡用的链子| 悠悠久久av| 欧美黑人精品巨大| 久久久国产精品麻豆| 国产成人影院久久av| 狂野欧美激情性xxxx| 亚洲 欧美一区二区三区| 国产精品久久久久久精品电影小说| 国产高清视频在线播放一区 | 久久久欧美国产精品| 久久性视频一级片| 又黄又粗又硬又大视频| 欧美老熟妇乱子伦牲交| 亚洲欧洲国产日韩| a级毛片黄视频| 亚洲第一青青草原| 成人国产av品久久久| 国产精品免费大片| 国产欧美日韩一区二区三区在线| 你懂的网址亚洲精品在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美一区二区三区久久| 国产亚洲精品第一综合不卡| 人人妻人人爽人人添夜夜欢视频| 日本a在线网址| 久久99精品国语久久久| 日本猛色少妇xxxxx猛交久久| 我的亚洲天堂| 欧美成人午夜精品| avwww免费| 亚洲成av片中文字幕在线观看| 天天影视国产精品| 午夜激情av网站| 亚洲中文字幕日韩| 国产在线一区二区三区精| av欧美777| 亚洲国产日韩一区二区| 国产福利在线免费观看视频| 国产色视频综合| 各种免费的搞黄视频| 两个人看的免费小视频| 欧美日韩国产mv在线观看视频| 下体分泌物呈黄色| 亚洲伊人色综图| 亚洲国产精品一区二区三区在线| 成人三级做爰电影| 欧美人与性动交α欧美软件| 女人高潮潮喷娇喘18禁视频| 久久久久精品人妻al黑| 这个男人来自地球电影免费观看| 制服诱惑二区| 制服人妻中文乱码| 国产黄色视频一区二区在线观看| 极品少妇高潮喷水抽搐| 大香蕉久久网| 91精品国产国语对白视频| 欧美亚洲日本最大视频资源| 免费在线观看完整版高清| 这个男人来自地球电影免费观看| 777米奇影视久久| 女人爽到高潮嗷嗷叫在线视频| 91国产中文字幕| 亚洲激情五月婷婷啪啪| 丰满人妻熟妇乱又伦精品不卡| 精品福利观看| 男女床上黄色一级片免费看| 婷婷色综合大香蕉| 午夜福利视频精品| 国产深夜福利视频在线观看| 亚洲av电影在线进入| 日本wwww免费看| 91九色精品人成在线观看| videos熟女内射| 成人亚洲欧美一区二区av| 日韩中文字幕视频在线看片| 免费在线观看视频国产中文字幕亚洲 | 极品人妻少妇av视频| av电影中文网址| 免费少妇av软件| 欧美精品人与动牲交sv欧美| 9色porny在线观看| 建设人人有责人人尽责人人享有的| 日韩欧美一区视频在线观看| 欧美在线一区亚洲| 欧美另类一区| 亚洲中文字幕日韩| 久久久欧美国产精品| 只有这里有精品99| 男女高潮啪啪啪动态图| 久久久久国产精品人妻一区二区| 自线自在国产av| 视频区欧美日本亚洲| 日韩制服骚丝袜av| 人妻人人澡人人爽人人| 男人添女人高潮全过程视频| 国产精品国产三级专区第一集| 成年女人毛片免费观看观看9 | 国产成人a∨麻豆精品| 性色av一级| kizo精华| 满18在线观看网站| 男女国产视频网站| 亚洲欧美日韩另类电影网站| 热99久久久久精品小说推荐| 亚洲精品日韩在线中文字幕| 黑人猛操日本美女一级片| 天堂俺去俺来也www色官网| 久久精品熟女亚洲av麻豆精品| 国产真人三级小视频在线观看| 国产av精品麻豆| 成人国产av品久久久| 国产精品麻豆人妻色哟哟久久| av有码第一页| 操美女的视频在线观看| 亚洲精品国产av成人精品| 欧美日韩亚洲高清精品| 十八禁网站网址无遮挡| 亚洲精品日韩在线中文字幕| 国产激情久久老熟女| 热re99久久国产66热| 亚洲av成人不卡在线观看播放网 | 久久久欧美国产精品| 天堂中文最新版在线下载| 最黄视频免费看| av一本久久久久| 久久久久久免费高清国产稀缺| 亚洲男人天堂网一区| 亚洲天堂av无毛| 久久人人爽av亚洲精品天堂| 亚洲精品在线美女| 免费在线观看完整版高清| 中文欧美无线码| 亚洲情色 制服丝袜| 性色av一级| 国产成人一区二区在线| 成人黄色视频免费在线看| 人人妻人人添人人爽欧美一区卜| 免费日韩欧美在线观看| 精品亚洲成a人片在线观看| 久久亚洲国产成人精品v| 欧美精品人与动牲交sv欧美| 亚洲免费av在线视频| 黄色毛片三级朝国网站| 欧美 日韩 精品 国产| 一本一本久久a久久精品综合妖精| 90打野战视频偷拍视频| 99香蕉大伊视频| 亚洲,一卡二卡三卡| 国产深夜福利视频在线观看| 中文字幕人妻丝袜制服| 久久久精品国产亚洲av高清涩受| av国产精品久久久久影院| 亚洲伊人色综图| 亚洲七黄色美女视频| 中国国产av一级| 女人被躁到高潮嗷嗷叫费观| 久久久精品免费免费高清| 纵有疾风起免费观看全集完整版| 免费不卡黄色视频| 99久久综合免费| 午夜福利影视在线免费观看| 久久久久久久国产电影| 亚洲成人免费av在线播放| 欧美激情 高清一区二区三区| 在线看a的网站| 各种免费的搞黄视频| 美女福利国产在线| 丝袜在线中文字幕| 日本黄色日本黄色录像| 中文欧美无线码| 国产片内射在线| 我要看黄色一级片免费的| 涩涩av久久男人的天堂| 美女脱内裤让男人舔精品视频| 一区二区av电影网| 日本欧美国产在线视频| 久久精品熟女亚洲av麻豆精品| 亚洲一码二码三码区别大吗| 国产视频首页在线观看| 国产精品久久久久久精品古装| 国产91精品成人一区二区三区 | 天堂俺去俺来也www色官网| 蜜桃在线观看..| 别揉我奶头~嗯~啊~动态视频 | 欧美乱码精品一区二区三区| 国产又色又爽无遮挡免| 69精品国产乱码久久久| av国产精品久久久久影院| 精品亚洲成a人片在线观看| 日本a在线网址| 国产一区二区 视频在线| 99久久人妻综合| 男女边摸边吃奶| 欧美另类一区| 亚洲一区中文字幕在线| 午夜影院在线不卡| 亚洲精品一二三| 国产女主播在线喷水免费视频网站| 一级毛片 在线播放| 国产精品熟女久久久久浪| 18禁国产床啪视频网站| 欧美精品一区二区免费开放| 国产1区2区3区精品| 嫩草影视91久久| 99精国产麻豆久久婷婷| 18禁裸乳无遮挡动漫免费视频| 国产女主播在线喷水免费视频网站| 免费女性裸体啪啪无遮挡网站| 日本vs欧美在线观看视频| 日韩 亚洲 欧美在线| 久久亚洲国产成人精品v| 18禁黄网站禁片午夜丰满| 97人妻天天添夜夜摸| 九草在线视频观看| 2018国产大陆天天弄谢| 午夜精品国产一区二区电影| 久久99一区二区三区| svipshipincom国产片| 欧美成人精品欧美一级黄| 色94色欧美一区二区| 成年女人毛片免费观看观看9 | 午夜福利在线免费观看网站| 国产精品秋霞免费鲁丝片| 丰满少妇做爰视频| 久久久精品免费免费高清| 最新的欧美精品一区二区| 丝袜美足系列| 日本91视频免费播放| 久久精品成人免费网站| 黄片播放在线免费| 国产野战对白在线观看| 国产精品 国内视频| 亚洲 欧美一区二区三区| 国产xxxxx性猛交| 制服诱惑二区| 在现免费观看毛片| 亚洲精品美女久久久久99蜜臀 | 中文字幕亚洲精品专区| 香蕉丝袜av| 天堂8中文在线网| 欧美日韩亚洲高清精品| 亚洲七黄色美女视频| 黄色怎么调成土黄色| 亚洲国产毛片av蜜桃av| 波多野结衣一区麻豆| 超碰成人久久| 日本黄色日本黄色录像| 亚洲欧洲精品一区二区精品久久久| 99久久人妻综合| 中文欧美无线码| 亚洲精品国产区一区二| 18禁黄网站禁片午夜丰满| 嫩草影视91久久| 日韩一卡2卡3卡4卡2021年| 十八禁高潮呻吟视频| 女人精品久久久久毛片| 久久狼人影院| 制服人妻中文乱码| 亚洲国产日韩一区二区| netflix在线观看网站| 美女脱内裤让男人舔精品视频| 大片免费播放器 马上看| 国产男女超爽视频在线观看| 亚洲国产最新在线播放| 久久99一区二区三区| 亚洲欧美中文字幕日韩二区| 日韩大码丰满熟妇| 女人被躁到高潮嗷嗷叫费观| 亚洲伊人久久精品综合| 欧美精品高潮呻吟av久久| 国产精品人妻久久久影院| 精品视频人人做人人爽| 国产日韩一区二区三区精品不卡| 日韩欧美一区视频在线观看| 精品第一国产精品| 精品少妇黑人巨大在线播放| 男男h啪啪无遮挡| 精品国产一区二区三区久久久樱花| 久久国产亚洲av麻豆专区| 久久女婷五月综合色啪小说| 一区二区三区乱码不卡18| 香蕉丝袜av| 免费日韩欧美在线观看| 国产黄色免费在线视频| 亚洲精品第二区| 性色av乱码一区二区三区2| netflix在线观看网站| 欧美激情高清一区二区三区| 国产成人精品久久二区二区免费| 欧美老熟妇乱子伦牲交| 久久精品久久精品一区二区三区| 五月开心婷婷网| 精品久久久精品久久久| 精品福利永久在线观看| 激情视频va一区二区三区| 两性夫妻黄色片| 午夜激情av网站| 亚洲欧美精品自产自拍| 国产免费视频播放在线视频| 国产欧美日韩一区二区三 | 成年美女黄网站色视频大全免费| 无限看片的www在线观看| 国产91精品成人一区二区三区 | 欧美黑人欧美精品刺激| √禁漫天堂资源中文www| 国产午夜精品一二区理论片| 久久久久久久精品精品| 美女脱内裤让男人舔精品视频| 咕卡用的链子| 电影成人av| 精品国产超薄肉色丝袜足j| 国产极品粉嫩免费观看在线| 免费在线观看影片大全网站 | 亚洲精品第二区| 亚洲美女黄色视频免费看| 91字幕亚洲| 一本—道久久a久久精品蜜桃钙片| 亚洲伊人色综图| 女警被强在线播放| 亚洲国产成人一精品久久久| 日韩av在线免费看完整版不卡| 少妇猛男粗大的猛烈进出视频| 看免费成人av毛片| av电影中文网址| 美女大奶头黄色视频| 啦啦啦啦在线视频资源| 性少妇av在线| 精品人妻在线不人妻| 久久热在线av| 欧美黄色片欧美黄色片| 搡老乐熟女国产| 免费看不卡的av| 亚洲国产中文字幕在线视频| www.自偷自拍.com| 69精品国产乱码久久久| 啦啦啦啦在线视频资源| 亚洲av国产av综合av卡| 丝袜脚勾引网站| 亚洲欧美一区二区三区黑人| 美国免费a级毛片| 99精国产麻豆久久婷婷| 九色亚洲精品在线播放| 精品熟女少妇八av免费久了| 免费少妇av软件| 亚洲精品国产色婷婷电影| 中文字幕人妻丝袜一区二区| 国产成人精品久久二区二区免费| 欧美xxⅹ黑人| 天堂俺去俺来也www色官网| 国产成人一区二区三区免费视频网站 | 天天影视国产精品| 亚洲欧洲精品一区二区精品久久久| 一本一本久久a久久精品综合妖精| 亚洲 欧美一区二区三区| 欧美激情高清一区二区三区| 国产精品国产三级国产专区5o| 免费高清在线观看视频在线观看| 亚洲国产成人一精品久久久| av又黄又爽大尺度在线免费看| 亚洲国产看品久久| 黄色怎么调成土黄色| 精品国产乱码久久久久久男人| 女性被躁到高潮视频| 中文欧美无线码| 97在线人人人人妻| 亚洲欧洲精品一区二区精品久久久| 亚洲第一青青草原| 午夜免费鲁丝| www.999成人在线观看| 久久精品久久精品一区二区三区| 久久毛片免费看一区二区三区| 久久人人爽av亚洲精品天堂| 亚洲,欧美精品.| 精品一区二区三区av网在线观看 | 一本综合久久免费| 欧美日韩成人在线一区二区| 男女午夜视频在线观看| 欧美大码av| 视频区欧美日本亚洲| a 毛片基地| 亚洲人成电影免费在线| 国产淫语在线视频| 国产av国产精品国产| 亚洲av日韩在线播放| 999精品在线视频| 国产成人欧美在线观看 | 亚洲中文日韩欧美视频| 日韩人妻精品一区2区三区| 久久久久久久国产电影| 免费久久久久久久精品成人欧美视频| 高清视频免费观看一区二区| 天天躁日日躁夜夜躁夜夜| 久久精品成人免费网站| av天堂在线播放| 欧美精品av麻豆av| 国产av精品麻豆| 国产日韩欧美视频二区| 国产97色在线日韩免费| 免费在线观看完整版高清| 人妻人人澡人人爽人人| 亚洲av欧美aⅴ国产| 在线观看www视频免费| 一本—道久久a久久精品蜜桃钙片| 日韩欧美一区视频在线观看| 国产成人a∨麻豆精品| 久久久久久久久久久久大奶| 日韩制服丝袜自拍偷拍| 伦理电影免费视频| 人体艺术视频欧美日本| 久久久精品区二区三区| 精品久久蜜臀av无| 久久鲁丝午夜福利片| 老司机亚洲免费影院| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人系列免费观看| 黄色怎么调成土黄色| 国产成人精品久久二区二区91| 热re99久久精品国产66热6| 中文字幕色久视频| 亚洲国产日韩一区二区| 精品福利观看| 菩萨蛮人人尽说江南好唐韦庄| 一区福利在线观看| 免费一级毛片在线播放高清视频 | 99国产综合亚洲精品| 久久人人爽人人片av| 七月丁香在线播放| 欧美黑人欧美精品刺激| 久久久久网色| 亚洲欧美成人综合另类久久久| 久久国产精品男人的天堂亚洲| 91老司机精品| 99国产精品一区二区蜜桃av | 中文字幕色久视频| 亚洲人成网站在线观看播放| 亚洲情色 制服丝袜| 久久久国产精品麻豆| www日本在线高清视频| 视频在线观看一区二区三区| 久久久精品区二区三区| av福利片在线| 丁香六月天网| 国产成人啪精品午夜网站| www.av在线官网国产| 久久久久国产精品人妻一区二区| 精品国产乱码久久久久久小说| 巨乳人妻的诱惑在线观看| 啦啦啦视频在线资源免费观看| 成年美女黄网站色视频大全免费| 久久久久久免费高清国产稀缺| 性少妇av在线| 精品少妇黑人巨大在线播放| 亚洲欧美中文字幕日韩二区| 亚洲国产av影院在线观看| 亚洲精品一二三| 免费日韩欧美在线观看| 亚洲欧洲国产日韩| 久久av网站| 男人爽女人下面视频在线观看| 男人操女人黄网站| 午夜免费观看性视频| 亚洲国产最新在线播放| 男人操女人黄网站| 午夜福利一区二区在线看| 久久精品国产综合久久久| 国产精品一国产av| 久久久精品免费免费高清| 日韩免费高清中文字幕av| 中文字幕人妻熟女乱码| 蜜桃在线观看..| 欧美精品啪啪一区二区三区 | 国产精品人妻久久久影院| 女性生殖器流出的白浆| 日韩中文字幕欧美一区二区 | 国产不卡av网站在线观看| 无限看片的www在线观看| 精品国产乱码久久久久久小说| 久久精品aⅴ一区二区三区四区| 精品高清国产在线一区| 久久99精品国语久久久|