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

    前緣鈍度對(duì)高速翼型顫振邊界的影響

    2020-03-05 04:25:20劉杰平呂計(jì)男王昕江劉子強(qiáng)
    宇航總體技術(shù) 2020年1期
    關(guān)鍵詞:氣動(dòng)彈性動(dòng)壓激波

    郭 力,劉杰平, 呂計(jì)男,王昕江,季 辰,劉子強(qiáng)

    (1.中國(guó)航天空氣動(dòng)力技術(shù)研究院, 北京 100074;2. 中國(guó)運(yùn)載火箭技術(shù)研究院,北京 100076)

    0 引言

    為了提高高速飛行器的裝載比例,其結(jié)構(gòu)的剛度與質(zhì)量下降,從而可能引起結(jié)構(gòu)的氣動(dòng)彈性振動(dòng)問題。針對(duì)高速飛行器的氣動(dòng)彈性問題,McNamara等[1]采用NASA Lanley研究中心的CFL3D程序計(jì)算了Mach數(shù)在2.5~28范圍內(nèi)高速飛行器的控制面與整機(jī)的氣動(dòng)彈性問題,并與活塞理論進(jìn)行了對(duì)比,分析了活塞理論與CFD/CSD耦合計(jì)算的誤差范圍。Lamorte等[2]針對(duì)高速飛行器的控制面進(jìn)行了氣動(dòng)彈性計(jì)算,考察了真實(shí)氣體效應(yīng)、湍流模型、轉(zhuǎn)捩位置對(duì)于顫振邊界的影響,發(fā)現(xiàn)上述效應(yīng)通過影響熱環(huán)境,影響氣動(dòng)彈性性能。Ilie等[3-5]采用大渦模擬(LES)方法,對(duì)平板及翼型的氣動(dòng)彈性響應(yīng)問題開展了算例研究。季辰等[6]發(fā)展了在高超聲速風(fēng)洞中對(duì)結(jié)構(gòu)進(jìn)行氣動(dòng)彈性實(shí)驗(yàn)的方法,研究了鈍前緣梯形翼[7]、舵面[8]的氣動(dòng)彈性特性,并獲得了梯形舵面的顫振邊界。竇怡彬等[9]采用分層求解技術(shù)研究了根部固支的高超聲速舵面顫振響應(yīng)和頻率。葉坤等[10]也采用分層求解技術(shù)分析了全動(dòng)舵面結(jié)構(gòu)的熱邊界影響。

    國(guó)內(nèi)外學(xué)者對(duì)于高超聲速環(huán)境下的舵面氣動(dòng)彈性問題進(jìn)行了數(shù)值、實(shí)驗(yàn)探索,所研究的對(duì)象具有工程實(shí)際意義,但所研究的外形變化較大,較難形成系統(tǒng)性的認(rèn)識(shí)。在高超聲速條件下,為了減小高超聲速時(shí)激波產(chǎn)生的阻力,高速飛行器的舵翼面大多采用尖前緣構(gòu)型,同時(shí)為了降低加熱對(duì)結(jié)構(gòu)的影響,需要對(duì)前緣進(jìn)行鈍化。本文以二維翼型為研究對(duì)象,通過變化翼型前緣的鈍化程度,探討單純的鈍度因素對(duì)翼型氣動(dòng)彈性性能的影響,為高超聲速下翼型前緣鈍度對(duì)顫振特性的影響研究提供參考。

    1 計(jì)算模型

    本文采用二維翼型作為模型考察高速下不同前緣鈍度對(duì)顫振邊界的影響。本研究采用了3個(gè)不同的前緣鈍度進(jìn)行了模擬,如圖1所示。第一種為尖前緣構(gòu)型(圖1(a)),此構(gòu)型在高速來流條件下會(huì)在前緣產(chǎn)生斜激波,激波的強(qiáng)度與正激波相比較小,波后的壓力增加與熱流較小。后兩種分別為1/60m半徑前緣(圖1(b))與1/30m半徑前緣構(gòu)型(圖1(c))。這兩種構(gòu)型會(huì)在前緣產(chǎn)生較強(qiáng)的正激波,產(chǎn)生的阻力較大。但是此類鈍前緣構(gòu)型,有利于結(jié)構(gòu)的防熱設(shè)計(jì),在氣動(dòng)加熱效應(yīng)比較明顯的高速速度范圍能夠較好地保持結(jié)構(gòu)的完整性。

    (a)

    (b)

    (c)圖1 計(jì)算采用的翼型截面形狀Fig.1 Illustration of wing sections

    計(jì)算來流Mach數(shù)為5,參考高度為10km。為了得到翼型結(jié)構(gòu)的顫振邊界,通過增加來流密度的方法增加動(dòng)壓。在高度10km時(shí),Mach為5,來流空氣的動(dòng)壓為7.6×105Pa,為了得到顫振邊界,又分別考慮了來流動(dòng)壓為15.1×105Pa、30.2×105Pa、60.5×105Pa這3種工況。通過對(duì)這4種情況的考慮,得到了各個(gè)來流動(dòng)壓下的時(shí)域響應(yīng)曲線,并根據(jù)這些曲線提取阻尼比,通過插值得到阻尼比為0時(shí)的動(dòng)壓即為顫振邊界。

    結(jié)構(gòu)模型如圖2所示。翼型結(jié)構(gòu)為剛體,具有沉浮與俯仰兩個(gè)自由度。翼型的轉(zhuǎn)動(dòng)中心在翼型的中心,質(zhì)心在轉(zhuǎn)動(dòng)中心前0.1m。質(zhì)量為10kg,轉(zhuǎn)動(dòng)慣量為14kg·m2,沉浮剛度為2×105N·m,扭轉(zhuǎn)剛度為1×106N·m2。在計(jì)算過程中,結(jié)構(gòu)與氣動(dòng)力耦合,在俯仰與沉浮兩個(gè)自由度上運(yùn)動(dòng)。

    圖2 結(jié)構(gòu)模型示意圖Fig.2 Illustration of structure model

    在實(shí)際計(jì)算過程中,由于考察的是系統(tǒng)的顫振邊界,只需要考慮結(jié)構(gòu)振動(dòng)的趨勢(shì),從中得到振動(dòng)的阻尼比,所以結(jié)構(gòu)運(yùn)動(dòng)的振幅不需要較大,只需要考慮小振幅就可以。根據(jù)此小振幅假設(shè),可以將振動(dòng)的特征向量進(jìn)行線性化,得到線性的振動(dòng)幅度,從而對(duì)轉(zhuǎn)軸的轉(zhuǎn)動(dòng)位移進(jìn)行線性組合運(yùn)算。翼型的質(zhì)心在15%弦長(zhǎng)處,翼型質(zhì)量與轉(zhuǎn)動(dòng)慣量均設(shè)為1,由此得到翼型的沉浮與俯仰兩階模態(tài)的頻率分別為ω1和ω2,如表1所示。

    表1 兩階模態(tài)的固有頻率

    2 數(shù)值方法

    為了方便描述動(dòng)網(wǎng)格上的流體運(yùn)動(dòng),流場(chǎng)采用積分形式的任意拉格朗日-歐拉(ALE)N-S方程[11]描述,其表達(dá)式如下

    (1)

    (2)

    (3)

    式(1)~式(3)中,ρ為密度,vi為各個(gè)方向的速度,vm,j為各個(gè)方向的控制體表面速度,p為壓力,E為總能量,T為溫度。計(jì)算過程中空間離散采用Roe格式,時(shí)間推進(jìn)采用雙時(shí)間步的LU-SGS方法。

    z=ψzψ

    (4)

    模態(tài)坐標(biāo)的方程為

    (5)

    3 顫振分析結(jié)果

    3.1 流場(chǎng)分析

    在Mach數(shù)為5.0時(shí),各個(gè)構(gòu)型的流場(chǎng)如圖3所示。對(duì)于尖前緣的翼型,在翼型周圍形成了斜激波,激波的傾角由尖前緣的前緣角決定。對(duì)于鈍前緣的兩個(gè)截面,在前緣前方形成了弓形激波。前緣鈍度的增加,使得前面的弓形激波強(qiáng)度增加,并且后面的斜激波的傾角增加。

    (a)

    (b)

    (c)圖3 不同翼型截面與對(duì)應(yīng)的定常流場(chǎng)Fig.3 Flow fields of the 3 wing section shapes

    3.2 不同來流動(dòng)壓下的時(shí)域數(shù)據(jù)

    計(jì)算初始時(shí)刻對(duì)翼型施加初始速度擾動(dòng),通過計(jì)算得到其時(shí)域響應(yīng)數(shù)據(jù),對(duì)得到的時(shí)域數(shù)據(jù)進(jìn)行分析可以得到系統(tǒng)的阻尼比和特征頻率。

    (1) 尖前緣外形變動(dòng)壓得到的時(shí)域數(shù)據(jù)

    從圖4中看到,隨著動(dòng)壓的增加,振動(dòng)幅度增大,并最終發(fā)散。

    (2) 1/60m半徑前緣外形變動(dòng)壓得到的時(shí)域數(shù)據(jù)

    從圖5中看到,隨著動(dòng)壓的增加,振動(dòng)幅度增大,并最終發(fā)散。1/60m半徑前緣的時(shí)域數(shù)據(jù)與尖前緣的時(shí)域數(shù)據(jù)相比,從定性角度看到發(fā)散較慢。

    (a)

    (b)

    (c)

    (d)

    (a)

    (b)

    (c)

    (d)

    (3) 1/30m半徑前緣外形變動(dòng)壓得到的時(shí)域數(shù)據(jù)

    通過對(duì)上面兩組計(jì)算從定性角度觀察,即使對(duì)于相同的動(dòng)壓增長(zhǎng)趨勢(shì),1/60m半徑前緣與尖前緣相比發(fā)散趨勢(shì)減小。反映到計(jì)算顫振邊界上,可以預(yù)計(jì),前緣鈍度的增加會(huì)使得顫振邊界增加。為了進(jìn)一步確認(rèn)此趨勢(shì),更進(jìn)一步增加了前緣的鈍度,采用圖1(c)中的翼型進(jìn)行了顫振計(jì)算。計(jì)算采用的Mach數(shù)與動(dòng)壓增加梯度與前面兩組計(jì)算相同,并通過固定在翼型前緣的點(diǎn)的位移得到了時(shí)域數(shù)據(jù)。

    1/30m半徑前緣的計(jì)算的時(shí)域數(shù)據(jù)如圖6所示,與尖前緣與1/60m半徑前緣的計(jì)算相比,從定性上看,隨著動(dòng)壓的增加,發(fā)散的趨勢(shì)更加不明顯,同時(shí)在最大動(dòng)壓點(diǎn)時(shí)域數(shù)據(jù)未顯示出發(fā)散的趨勢(shì),而是呈現(xiàn)出等幅振動(dòng)的趨勢(shì)。

    3.3 顫振頻率與阻尼比分析

    3.3.1 銳前緣翼型

    采用自回歸移動(dòng)平均方法[13](ARMA)對(duì)銳前緣翼型前緣的位移時(shí)域數(shù)據(jù)進(jìn)行了分析,分別得到了不同動(dòng)壓下的頻率與阻尼比,如表2所示。時(shí)域分析得到的頻譜圖如圖7所示。從頻率上可以看出,隨著動(dòng)壓的增加,最高頻率與最低頻率逐漸靠近。在顫振沒有發(fā)生時(shí),時(shí)域數(shù)據(jù)的頻譜分析具有兩個(gè)峰值,分別對(duì)應(yīng)兩個(gè)模態(tài)的振動(dòng)頻率。當(dāng)顫振發(fā)生時(shí),時(shí)域數(shù)據(jù)發(fā)散,此時(shí)只有一個(gè)對(duì)應(yīng)的頻率。

    將來流動(dòng)壓表示為阻尼比的函數(shù),則在阻尼比為0時(shí)對(duì)應(yīng)的動(dòng)壓為顫振動(dòng)壓。阻尼比與顫振動(dòng)壓的函數(shù)關(guān)系通過已知的動(dòng)壓-阻尼比數(shù)據(jù)插值得到。相應(yīng)的插值曲線與數(shù)據(jù)點(diǎn)如圖8所示。根據(jù)此函數(shù)關(guān)系,得到阻尼比為0時(shí)對(duì)應(yīng)的顫振動(dòng)壓為3.01×106Pa。

    (a)

    (b)

    (c)

    (d)

    表2 銳前緣翼型頻率與阻尼比

    (a)

    (b)

    (c)

    (d)

    圖8 銳前緣翼型阻尼比與顫振動(dòng)壓Fig.8 Sharp-edge wing damping ratio vs. incoming flow dynamic pressure

    3.3.2 1/60m半徑前緣翼型

    對(duì)1/60m半徑前緣翼型前緣的位移時(shí)域數(shù)據(jù)進(jìn)行分析,分別得到了不同動(dòng)壓下的頻率與阻尼比,如表3所示。時(shí)域分析得到的頻譜圖如圖9所示。從頻率上可以看出,隨著動(dòng)壓的增加,最高頻率與最低頻率逐漸靠近。在顫振沒有發(fā)生時(shí),時(shí)域數(shù)據(jù)的頻譜分析具有兩個(gè)峰值,分別對(duì)應(yīng)兩個(gè)模態(tài)的振動(dòng)頻率。當(dāng)顫振發(fā)生時(shí),時(shí)域數(shù)據(jù)發(fā)散,此時(shí)只有一個(gè)對(duì)應(yīng)的頻率。

    表3 1/60m半徑前緣翼型頻率與阻尼比

    (b)

    (c)

    (d)圖9 1/60m半徑前緣翼型頻率分析Fig.9 Frequency of 1/60m blunted-edge wing

    將來流動(dòng)壓表示為阻尼比的函數(shù),則在阻尼比為0時(shí)對(duì)應(yīng)的動(dòng)壓為顫振動(dòng)壓。阻尼比與顫振動(dòng)壓的函數(shù)關(guān)系通過已知的動(dòng)壓-阻尼比數(shù)據(jù)插值得到,相應(yīng)的插值曲線與數(shù)據(jù)點(diǎn)如圖10所示。根據(jù)此函數(shù)關(guān)系,得到阻尼比為0時(shí)對(duì)應(yīng)的顫振動(dòng)壓為3.25×106Pa。

    圖10 1/60m半徑前緣翼型阻尼比與顫振動(dòng)壓Fig.10 1/60m blunted-edge damping ratio vs. incoming flow dynamic pressure

    3.3.3 1/30m半徑前緣翼型

    同樣對(duì)1/30m半徑前緣翼型前緣的位移時(shí)域數(shù)據(jù)進(jìn)行了分析,分別得到了不同動(dòng)壓下的頻率與阻尼比,如表4所示。時(shí)域分析得到的頻譜圖如圖11所示。從頻率上可以看出,隨著動(dòng)壓的增加,最高頻率與最低頻率逐漸靠近。在顫振沒有發(fā)生時(shí),時(shí)域數(shù)據(jù)的頻譜分析具有兩個(gè)峰值,分別對(duì)應(yīng)兩個(gè)模態(tài)的振動(dòng)頻率。與銳前緣和1/60m半徑前緣算例不同的是,1/30m半徑前緣即使在8倍動(dòng)壓的情況下仍沒有發(fā)散。

    表4 1/30m半徑前緣翼型頻率與阻尼比

    將來流動(dòng)壓表示為阻尼比的函數(shù),相應(yīng)的數(shù)據(jù)點(diǎn)如圖12所示??梢钥闯?即使在最大動(dòng)壓狀態(tài)下,振動(dòng)仍沒有發(fā)散,因此1/30m半徑前緣的翼型顫振邊界高于前面兩種構(gòu)型。雖然在計(jì)算設(shè)置的動(dòng)壓情況下沒有得到結(jié)構(gòu)振動(dòng)發(fā)散,從而沒有得到具體的顫振邊界數(shù)據(jù),但是此結(jié)果與前面結(jié)果的趨勢(shì)一致,說明了隨著翼型前緣鈍度的增加,顫振邊界逐漸提高。

    (a)

    (b)

    (c)

    (d)圖11 1/30m半徑前緣翼型頻率分析Fig.11 Frequency of 1/30m blunted-edge wing

    圖12 1/30m半徑前緣翼型阻尼比與顫振動(dòng)壓Fig.12 1/30m blunted-edge damping ratio vs. incoming flow dynamic pressure

    4 結(jié)論

    根據(jù)從銳前緣到1/30m半徑前緣的氣動(dòng)彈性分析,發(fā)現(xiàn)隨著前緣鈍度的增加,即使翼型的剛度、質(zhì)量不發(fā)生改變,顫振邊界也會(huì)不斷提高。隨著翼型前緣鈍度的增加,在高速下,頭部激波結(jié)構(gòu)從斜激波轉(zhuǎn)變?yōu)槊擉w的弓形激波。激波結(jié)構(gòu)的變化使得翼型在運(yùn)動(dòng)情況下產(chǎn)生的使得翼型發(fā)生俯仰運(yùn)動(dòng)的力矩發(fā)生了變化。而只有翼型在平動(dòng)過程中產(chǎn)生力矩并與沉浮運(yùn)動(dòng)耦合,才會(huì)出現(xiàn)翼型結(jié)構(gòu)的自激振動(dòng),從而發(fā)散出現(xiàn)顫振現(xiàn)象。在銳前緣情況下,力矩較大,隨著前緣鈍度的增加,翼型上受到的力矩逐漸減小。一個(gè)極端的情況為球形,此時(shí)如果忽略黏性作用,球形受到的力矩為0。因此,增大前緣鈍度可以提高顫振邊界。

    猜你喜歡
    氣動(dòng)彈性動(dòng)壓激波
    國(guó)內(nèi)首個(gè)現(xiàn)代箔片氣體動(dòng)壓軸承技術(shù)培訓(xùn)班在長(zhǎng)沙成功舉辦
    一種基于聚類分析的二維激波模式識(shí)別算法
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    斜激波入射V形鈍前緣溢流口激波干擾研究
    適于可壓縮多尺度流動(dòng)的緊致型激波捕捉格式
    飛翼無人機(jī)嗡鳴氣動(dòng)彈性響應(yīng)分析
    南屯煤礦深部泵房硐室群動(dòng)壓失穩(wěn)機(jī)理及控制對(duì)策
    強(qiáng)烈動(dòng)壓巷道支護(hù)技術(shù)探討
    模態(tài)選取對(duì)靜氣動(dòng)彈性分析的影響
    直升機(jī)的氣動(dòng)彈性問題
    免费在线观看黄色视频的| 91字幕亚洲| 两性夫妻黄色片| 久久人人精品亚洲av| 91九色精品人成在线观看| x7x7x7水蜜桃| 性色av乱码一区二区三区2| 国产不卡一卡二| av电影中文网址| 18美女黄网站色大片免费观看| 久久国产精品男人的天堂亚洲| 亚洲av五月六月丁香网| 欧洲精品卡2卡3卡4卡5卡区| 两性夫妻黄色片| 久久久精品欧美日韩精品| 一级黄色大片毛片| 国产精品 欧美亚洲| 亚洲av熟女| 在线观看一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 欧美色欧美亚洲另类二区 | 久久精品成人免费网站| www.自偷自拍.com| 亚洲av成人一区二区三| 91麻豆av在线| 午夜免费激情av| 两个人免费观看高清视频| 淫妇啪啪啪对白视频| 两个人视频免费观看高清| 亚洲激情在线av| 国语自产精品视频在线第100页| 我的亚洲天堂| 亚洲五月天丁香| 美女高潮喷水抽搐中文字幕| 一级a爱视频在线免费观看| 久久青草综合色| 欧美成人免费av一区二区三区| 18美女黄网站色大片免费观看| 亚洲成人国产一区在线观看| 欧美日韩精品网址| 亚洲欧美一区二区三区黑人| 亚洲av成人不卡在线观看播放网| 咕卡用的链子| 国产麻豆69| 视频区欧美日本亚洲| 成人av一区二区三区在线看| 久久热在线av| 亚洲欧美激情综合另类| 长腿黑丝高跟| 欧美日韩瑟瑟在线播放| 亚洲国产中文字幕在线视频| 给我免费播放毛片高清在线观看| 日韩三级视频一区二区三区| 侵犯人妻中文字幕一二三四区| 自拍欧美九色日韩亚洲蝌蚪91| 1024视频免费在线观看| 99国产综合亚洲精品| 亚洲全国av大片| 亚洲va日本ⅴa欧美va伊人久久| 高清在线国产一区| 国产亚洲av嫩草精品影院| 美女高潮到喷水免费观看| 成人av一区二区三区在线看| 日韩大尺度精品在线看网址 | 青草久久国产| 美女高潮到喷水免费观看| 亚洲熟妇熟女久久| 国产精品久久久av美女十八| 国产精品1区2区在线观看.| 欧美乱码精品一区二区三区| 欧美大码av| 久久草成人影院| 高清黄色对白视频在线免费看| 亚洲国产精品久久男人天堂| 精品午夜福利视频在线观看一区| 琪琪午夜伦伦电影理论片6080| 国产高清激情床上av| 国产三级在线视频| 女同久久另类99精品国产91| √禁漫天堂资源中文www| 色播亚洲综合网| 午夜亚洲福利在线播放| 深夜精品福利| 国产欧美日韩一区二区三区在线| 女同久久另类99精品国产91| 搡老熟女国产l中国老女人| 国产成人啪精品午夜网站| 国产欧美日韩综合在线一区二区| 久久久久久久午夜电影| 久久婷婷人人爽人人干人人爱 | 亚洲专区中文字幕在线| 久久精品亚洲熟妇少妇任你| 国产片内射在线| 亚洲天堂国产精品一区在线| 18禁黄网站禁片午夜丰满| 麻豆av在线久日| 久久久国产精品麻豆| 亚洲精品久久成人aⅴ小说| 国产av在哪里看| 成年版毛片免费区| 成人欧美大片| 精品一区二区三区av网在线观看| 好男人在线观看高清免费视频 | 女同久久另类99精品国产91| av电影中文网址| √禁漫天堂资源中文www| 99久久精品国产亚洲精品| 国产成人影院久久av| 国产一卡二卡三卡精品| 波多野结衣av一区二区av| 久久香蕉精品热| 啦啦啦观看免费观看视频高清 | 一级毛片精品| bbb黄色大片| 女人精品久久久久毛片| 免费不卡黄色视频| 大码成人一级视频| 色尼玛亚洲综合影院| 亚洲 欧美一区二区三区| 国产激情久久老熟女| 成人亚洲精品av一区二区| 一级片免费观看大全| 国产午夜精品久久久久久| 女人高潮潮喷娇喘18禁视频| 午夜视频精品福利| 国产免费男女视频| 亚洲午夜理论影院| 多毛熟女@视频| 欧美黄色淫秽网站| 电影成人av| 亚洲狠狠婷婷综合久久图片| 99re在线观看精品视频| 国产99白浆流出| 女同久久另类99精品国产91| 一二三四在线观看免费中文在| 亚洲中文日韩欧美视频| 国内久久婷婷六月综合欲色啪| 亚洲欧美精品综合久久99| 中国美女看黄片| 男女床上黄色一级片免费看| 巨乳人妻的诱惑在线观看| 黄网站色视频无遮挡免费观看| 午夜福利影视在线免费观看| 可以在线观看毛片的网站| 青草久久国产| e午夜精品久久久久久久| 国产亚洲欧美精品永久| 精品一区二区三区av网在线观看| 看片在线看免费视频| 亚洲自偷自拍图片 自拍| 久久国产精品影院| 999久久久国产精品视频| 国产野战对白在线观看| 精品福利观看| 久久久久精品国产欧美久久久| 欧美一级毛片孕妇| 青草久久国产| 亚洲五月天丁香| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品中文字幕一二三四区| 久久久久久久午夜电影| 国产精品香港三级国产av潘金莲| 99re在线观看精品视频| 一区二区三区高清视频在线| 欧美日韩瑟瑟在线播放| 一本大道久久a久久精品| 国产伦一二天堂av在线观看| 男男h啪啪无遮挡| 精品午夜福利视频在线观看一区| 国产精品永久免费网站| 国产av又大| 日日干狠狠操夜夜爽| 亚洲av日韩精品久久久久久密| 国产亚洲欧美在线一区二区| 亚洲国产精品sss在线观看| 久久久久久九九精品二区国产| 国产亚洲精品久久久久久毛片| 一个人看的www免费观看视频| 久久国产乱子免费精品| 国产精品人妻久久久久久| 国产在视频线在精品| 夜夜夜夜夜久久久久| 一本久久中文字幕| 色播亚洲综合网| 欧美日韩乱码在线| 变态另类丝袜制服| 国产精品一及| 久久99热这里只有精品18| 国产精品福利在线免费观看| av.在线天堂| 黄色日韩在线| 搡老熟女国产l中国老女人| 人人妻人人看人人澡| 精品99又大又爽又粗少妇毛片 | 桃色一区二区三区在线观看| 国产欧美日韩一区二区精品| 日韩高清综合在线| 一夜夜www| 男女边吃奶边做爰视频| 毛片一级片免费看久久久久 | 在线看三级毛片| 国产精品久久久久久久久免| 国产久久久一区二区三区| 国产探花在线观看一区二区| 久久精品国产亚洲av涩爱 | 亚洲成人免费电影在线观看| 麻豆成人av在线观看| 两个人的视频大全免费| 97碰自拍视频| 精品久久久久久久久亚洲 | 日韩欧美国产一区二区入口| 国产精品久久久久久精品电影| 99热网站在线观看| 18禁黄网站禁片免费观看直播| 日韩欧美国产一区二区入口| 天天一区二区日本电影三级| 成人亚洲精品av一区二区| 成人美女网站在线观看视频| 在线观看66精品国产| 国产激情偷乱视频一区二区| av.在线天堂| 最新中文字幕久久久久| 91午夜精品亚洲一区二区三区 | 国产精品av视频在线免费观看| 2021天堂中文幕一二区在线观| 久久久久久久久久黄片| 国产高清有码在线观看视频| 欧美成人a在线观看| 在线免费观看不下载黄p国产 | 亚洲av二区三区四区| 亚洲成人免费电影在线观看| 中文字幕av成人在线电影| 我的老师免费观看完整版| 18禁黄网站禁片免费观看直播| 一a级毛片在线观看| www.色视频.com| 丰满的人妻完整版| 亚洲18禁久久av| 国产男人的电影天堂91| 少妇猛男粗大的猛烈进出视频 | 午夜亚洲福利在线播放| 男人的好看免费观看在线视频| 99久国产av精品| 欧美日本视频| 两个人的视频大全免费| 99热网站在线观看| 亚洲七黄色美女视频| 亚洲人与动物交配视频| 久久精品国产鲁丝片午夜精品 | 99久久精品热视频| 丰满乱子伦码专区| www.www免费av| 桃色一区二区三区在线观看| 联通29元200g的流量卡| 一级黄色大片毛片| 深夜a级毛片| 狂野欧美激情性xxxx在线观看| 免费观看的影片在线观看| 精品日产1卡2卡| 中文字幕高清在线视频| 1000部很黄的大片| 久久久久久久久久黄片| 免费大片18禁| 精品人妻熟女av久视频| 99热这里只有是精品50| 在线a可以看的网站| 亚洲精品成人久久久久久| 色哟哟·www| 一本一本综合久久| 啦啦啦啦在线视频资源| 国产精品永久免费网站| 日韩在线高清观看一区二区三区 | 午夜日韩欧美国产| ponron亚洲| 丰满的人妻完整版| 一个人观看的视频www高清免费观看| av在线天堂中文字幕| 成人国产一区最新在线观看| 真实男女啪啪啪动态图| 很黄的视频免费| 国产成人aa在线观看| 51国产日韩欧美| 啦啦啦啦在线视频资源| 亚洲欧美精品综合久久99| 欧美日本视频| 男女啪啪激烈高潮av片| 一进一出好大好爽视频| 国产一区二区在线av高清观看| 日日干狠狠操夜夜爽| 亚洲四区av| 国产欧美日韩一区二区精品| 非洲黑人性xxxx精品又粗又长| 国产精品伦人一区二区| 91精品国产九色| 真实男女啪啪啪动态图| 中文字幕高清在线视频| 在线观看免费视频日本深夜| 国产高潮美女av| 九九久久精品国产亚洲av麻豆| 欧美成人性av电影在线观看| 午夜爱爱视频在线播放| 色尼玛亚洲综合影院| 91狼人影院| 黄色配什么色好看| 麻豆久久精品国产亚洲av| eeuss影院久久| 国产精品一及| 俄罗斯特黄特色一大片| 啪啪无遮挡十八禁网站| 日本欧美国产在线视频| 最近在线观看免费完整版| 久久久久国内视频| 亚洲成人久久爱视频| av在线老鸭窝| 少妇高潮的动态图| 日韩欧美国产一区二区入口| www.色视频.com| 久久精品综合一区二区三区| 亚洲最大成人手机在线| 国产日本99.免费观看| 免费一级毛片在线播放高清视频| 久久久久免费精品人妻一区二区| 欧美+亚洲+日韩+国产| 国产免费av片在线观看野外av| 国产国拍精品亚洲av在线观看| 九九热线精品视视频播放| 亚洲国产精品成人综合色| 欧美色视频一区免费| 亚洲在线自拍视频| 欧美日韩国产亚洲二区| 此物有八面人人有两片| 少妇的逼好多水| 日本免费a在线| 精品无人区乱码1区二区| 亚洲国产色片| 午夜福利视频1000在线观看| 亚洲欧美日韩高清专用| 精品日产1卡2卡| 免费一级毛片在线播放高清视频| 草草在线视频免费看| 免费av毛片视频| 搡老妇女老女人老熟妇| 午夜激情欧美在线| 国产欧美日韩精品一区二区| 九九久久精品国产亚洲av麻豆| 在线播放国产精品三级| 99久久精品一区二区三区| 欧美+亚洲+日韩+国产| 国产精品三级大全| 午夜激情福利司机影院| 18+在线观看网站| 色噜噜av男人的天堂激情| 欧美性感艳星| 亚洲第一电影网av| 国产高清不卡午夜福利| 99热网站在线观看| a在线观看视频网站| 亚洲天堂国产精品一区在线| 免费av不卡在线播放| 精华霜和精华液先用哪个| 日本爱情动作片www.在线观看 | 少妇的逼好多水| 99riav亚洲国产免费| 精品久久久久久久久久免费视频| 黄色日韩在线| 国产精品久久久久久精品电影| 日日干狠狠操夜夜爽| 老司机福利观看| 真人做人爱边吃奶动态| 悠悠久久av| 亚洲精品色激情综合| 一本久久中文字幕| 国产久久久一区二区三区| 国产黄片美女视频| 一进一出好大好爽视频| 精品欧美国产一区二区三| 亚洲中文日韩欧美视频| 久久国产乱子免费精品| videossex国产| 国内久久婷婷六月综合欲色啪| 亚洲国产高清在线一区二区三| 国产欧美日韩精品亚洲av| 免费在线观看成人毛片| 国产伦在线观看视频一区| 国产 一区精品| 毛片一级片免费看久久久久 | 人妻久久中文字幕网| 亚洲人成网站在线播放欧美日韩| 国产精品国产三级国产av玫瑰| 简卡轻食公司| 亚洲成a人片在线一区二区| 日本在线视频免费播放| 国产亚洲精品久久久久久毛片| 麻豆国产97在线/欧美| 精品午夜福利视频在线观看一区| 免费观看在线日韩| 国产一区二区在线av高清观看| 1024手机看黄色片| 国产精品野战在线观看| 午夜福利高清视频| 一个人观看的视频www高清免费观看| 日韩欧美在线二视频| 久久精品人妻少妇| 午夜福利在线观看吧| 国内久久婷婷六月综合欲色啪| 欧美最新免费一区二区三区| 日日撸夜夜添| av女优亚洲男人天堂| av在线亚洲专区| 男人舔女人下体高潮全视频| 禁无遮挡网站| 国产精品人妻久久久久久| 亚洲欧美日韩无卡精品| 岛国在线免费视频观看| 亚洲天堂国产精品一区在线| 我的老师免费观看完整版| 国产精品自产拍在线观看55亚洲| 热99re8久久精品国产| 又黄又爽又免费观看的视频| 十八禁网站免费在线| 亚洲人成伊人成综合网2020| 日韩精品中文字幕看吧| 亚洲一级一片aⅴ在线观看| 桃色一区二区三区在线观看| 在线a可以看的网站| 亚洲性夜色夜夜综合| 精品午夜福利视频在线观看一区| 国产熟女欧美一区二区| 淫秽高清视频在线观看| 直男gayav资源| 男人舔奶头视频| 亚洲欧美日韩卡通动漫| 大又大粗又爽又黄少妇毛片口| 日韩欧美精品v在线| 亚洲精品久久国产高清桃花| 久久久久久久久中文| av国产免费在线观看| av女优亚洲男人天堂| 成年版毛片免费区| 哪里可以看免费的av片| 欧美中文日本在线观看视频| .国产精品久久| 久久久久久伊人网av| 欧美色欧美亚洲另类二区| 成人国产综合亚洲| 99热只有精品国产| 热99re8久久精品国产| 少妇熟女aⅴ在线视频| 我的女老师完整版在线观看| 一本一本综合久久| 国产精品不卡视频一区二区| 国产免费av片在线观看野外av| 日韩亚洲欧美综合| 韩国av一区二区三区四区| 麻豆av噜噜一区二区三区| 国产 一区 欧美 日韩| av在线蜜桃| 免费大片18禁| 国产 一区 欧美 日韩| 中文字幕av成人在线电影| 欧美成人性av电影在线观看| 久久久久久久亚洲中文字幕| 大型黄色视频在线免费观看| 亚洲国产精品成人综合色| 国内精品久久久久精免费| 2021天堂中文幕一二区在线观| 欧美xxxx黑人xx丫x性爽| 国产精品电影一区二区三区| 日本与韩国留学比较| 欧美黑人巨大hd| 99热这里只有是精品50| 精品人妻1区二区| 一个人免费在线观看电影| 中文在线观看免费www的网站| 12—13女人毛片做爰片一| 欧美性猛交黑人性爽| 亚洲人成网站在线播| 国产免费一级a男人的天堂| 91久久精品国产一区二区成人| 午夜精品一区二区三区免费看| 日韩精品有码人妻一区| 国产视频内射| 99久久久亚洲精品蜜臀av| 少妇猛男粗大的猛烈进出视频 | 在线播放无遮挡| 国产精品爽爽va在线观看网站| 色综合婷婷激情| 极品教师在线免费播放| 成人精品一区二区免费| 亚洲av日韩精品久久久久久密| 波多野结衣巨乳人妻| 淫秽高清视频在线观看| 91狼人影院| 女人十人毛片免费观看3o分钟| 一个人看视频在线观看www免费| 老熟妇仑乱视频hdxx| 日本免费a在线| 99热精品在线国产| 久久久成人免费电影| 无人区码免费观看不卡| 国产精品一区二区性色av| 精品欧美国产一区二区三| 日本成人三级电影网站| 午夜福利视频1000在线观看| 成人特级黄色片久久久久久久| 春色校园在线视频观看| 最近最新中文字幕大全电影3| 午夜福利18| 97超级碰碰碰精品色视频在线观看| 男人舔奶头视频| 51国产日韩欧美| 成人av一区二区三区在线看| 亚洲四区av| av在线蜜桃| 精品一区二区三区视频在线| av.在线天堂| 免费大片18禁| 两性午夜刺激爽爽歪歪视频在线观看| 免费大片18禁| 黄色配什么色好看| 99国产精品一区二区蜜桃av| 韩国av一区二区三区四区| 国产精品三级大全| 99久久精品一区二区三区| 日韩欧美精品v在线| 一进一出抽搐动态| 国产免费男女视频| 亚洲av美国av| 亚洲第一电影网av| 久久国产乱子免费精品| 午夜激情福利司机影院| 国产一区二区三区视频了| 午夜老司机福利剧场| 国产探花极品一区二区| 国产精品精品国产色婷婷| 欧美一级a爱片免费观看看| 一进一出好大好爽视频| 亚洲avbb在线观看| 日韩欧美精品v在线| 在线免费观看不下载黄p国产 | 制服丝袜大香蕉在线| 日日摸夜夜添夜夜添av毛片 | 麻豆一二三区av精品| 长腿黑丝高跟| 精品久久久噜噜| 99久久中文字幕三级久久日本| 国内久久婷婷六月综合欲色啪| 久久亚洲精品不卡| 亚洲一区高清亚洲精品| 久久欧美精品欧美久久欧美| 51国产日韩欧美| 国产高清视频在线播放一区| 麻豆国产av国片精品| 久久精品91蜜桃| 99riav亚洲国产免费| 99视频精品全部免费 在线| 欧美成人一区二区免费高清观看| 欧美日韩瑟瑟在线播放| 免费人成视频x8x8入口观看| av国产免费在线观看| 免费在线观看成人毛片| 成人av一区二区三区在线看| 桃红色精品国产亚洲av| 亚洲人成网站在线播放欧美日韩| 久久久久国内视频| 色哟哟哟哟哟哟| 在线免费观看不下载黄p国产 | 欧美色视频一区免费| 国产精品三级大全| 亚洲av中文字字幕乱码综合| 国产一区二区三区在线臀色熟女| 国产伦在线观看视频一区| 亚洲aⅴ乱码一区二区在线播放| 欧美不卡视频在线免费观看| 精品久久久久久久末码| 又紧又爽又黄一区二区| 国产成人影院久久av| 色哟哟哟哟哟哟| 亚洲成人久久爱视频| 午夜免费成人在线视频| 色视频www国产| 日韩中字成人| 亚洲avbb在线观看| 一进一出好大好爽视频| 欧美国产日韩亚洲一区| 国产精品免费一区二区三区在线| 精品久久久久久久久久久久久| 日日啪夜夜撸| 人妻少妇偷人精品九色| 亚洲精品在线观看二区| 老司机福利观看| 国产高清激情床上av| 国产精品一区二区免费欧美| 国产伦精品一区二区三区视频9| 最近最新中文字幕大全电影3| 精品久久久久久久久av| 日韩欧美在线乱码| 欧美zozozo另类| 国产熟女欧美一区二区| 欧美色视频一区免费| 日韩人妻高清精品专区| 免费在线观看成人毛片| 少妇丰满av| av在线老鸭窝| 国产成人福利小说| 中文资源天堂在线| 日韩av在线大香蕉| .国产精品久久| 日本欧美国产在线视频| 日韩一本色道免费dvd| 两人在一起打扑克的视频| 色噜噜av男人的天堂激情| 黄色丝袜av网址大全| 亚洲国产欧洲综合997久久,| 午夜视频国产福利|