• <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)彈性問題
    又黄又粗又硬又大视频| 国产免费一级a男人的天堂| 国产亚洲精品久久久com| 亚洲av一区综合| 淫秽高清视频在线观看| 90打野战视频偷拍视频| 亚洲精华国产精华精| 嫁个100分男人电影在线观看| 最近在线观看免费完整版| 午夜激情欧美在线| 听说在线观看完整版免费高清| 18禁在线播放成人免费| 亚洲精品影视一区二区三区av| 性色avwww在线观看| 青草久久国产| 国产美女午夜福利| 国产av麻豆久久久久久久| 在线播放无遮挡| 国产探花极品一区二区| 12—13女人毛片做爰片一| 午夜福利在线在线| 久久精品国产清高在天天线| 久久精品综合一区二区三区| 亚洲av免费高清在线观看| 51国产日韩欧美| www国产在线视频色| 色av中文字幕| 国产欧美日韩精品一区二区| 黑人欧美特级aaaaaa片| 免费无遮挡裸体视频| 国语自产精品视频在线第100页| 夜夜看夜夜爽夜夜摸| 神马国产精品三级电影在线观看| 丰满的人妻完整版| 18美女黄网站色大片免费观看| 亚洲av免费高清在线观看| 一级作爱视频免费观看| 国产69精品久久久久777片| 女警被强在线播放| 成人鲁丝片一二三区免费| 可以在线观看的亚洲视频| 最近在线观看免费完整版| 首页视频小说图片口味搜索| 淫妇啪啪啪对白视频| netflix在线观看网站| 国产精品1区2区在线观看.| 国产爱豆传媒在线观看| 国产精品免费一区二区三区在线| 国产av不卡久久| 欧美乱妇无乱码| 两人在一起打扑克的视频| 18+在线观看网站| 黄色日韩在线| 亚洲欧美日韩无卡精品| 欧美zozozo另类| 中文字幕精品亚洲无线码一区| 法律面前人人平等表现在哪些方面| 99热这里只有是精品50| 日韩免费av在线播放| 最新中文字幕久久久久| 一夜夜www| 国产精品久久视频播放| 蜜桃亚洲精品一区二区三区| 中文字幕高清在线视频| 老汉色∧v一级毛片| 不卡一级毛片| 欧美黑人欧美精品刺激| av在线天堂中文字幕| 日日干狠狠操夜夜爽| 成人av在线播放网站| 日韩国内少妇激情av| 不卡一级毛片| 久久精品综合一区二区三区| 好男人在线观看高清免费视频| 91麻豆av在线| 成人欧美大片| 日本与韩国留学比较| 久久中文看片网| 深爱激情五月婷婷| 岛国在线免费视频观看| 老鸭窝网址在线观看| 老司机午夜十八禁免费视频| 人妻丰满熟妇av一区二区三区| 国产激情偷乱视频一区二区| 色在线成人网| 成人特级av手机在线观看| 18美女黄网站色大片免费观看| x7x7x7水蜜桃| 日本三级黄在线观看| 看黄色毛片网站| 亚洲成人免费电影在线观看| 一区福利在线观看| 极品教师在线免费播放| 久99久视频精品免费| 一本精品99久久精品77| 一本一本综合久久| 成人特级av手机在线观看| eeuss影院久久| 男人舔女人下体高潮全视频| 久久亚洲精品不卡| 久久久久久久久大av| 国产精品免费一区二区三区在线| 国产成年人精品一区二区| 欧美黄色片欧美黄色片| 免费电影在线观看免费观看| 久久中文看片网| 国产欧美日韩精品亚洲av| 亚洲激情在线av| 国产一区二区在线av高清观看| 老鸭窝网址在线观看| 欧美激情在线99| 亚洲欧美精品综合久久99| 亚洲成av人片在线播放无| 欧美三级亚洲精品| 欧美bdsm另类| 日日干狠狠操夜夜爽| 女警被强在线播放| 99久久九九国产精品国产免费| 男女午夜视频在线观看| 久久久国产精品麻豆| 一个人观看的视频www高清免费观看| tocl精华| 国产综合懂色| 免费av不卡在线播放| 日韩人妻高清精品专区| 88av欧美| 国产一区二区亚洲精品在线观看| 国产精品久久久久久久久免 | 天堂√8在线中文| 国内毛片毛片毛片毛片毛片| 国产单亲对白刺激| 亚洲av中文字字幕乱码综合| 麻豆久久精品国产亚洲av| 国产蜜桃级精品一区二区三区| av在线天堂中文字幕| 久久久久九九精品影院| 国产色婷婷99| 午夜两性在线视频| 深夜精品福利| 欧美黑人巨大hd| 国产高清videossex| 激情在线观看视频在线高清| 色噜噜av男人的天堂激情| 噜噜噜噜噜久久久久久91| 日韩高清综合在线| 亚洲成a人片在线一区二区| 亚洲欧美日韩高清在线视频| 女同久久另类99精品国产91| 国产高清视频在线观看网站| 久久天躁狠狠躁夜夜2o2o| 国产高清videossex| 欧美绝顶高潮抽搐喷水| 女人高潮潮喷娇喘18禁视频| 免费电影在线观看免费观看| 最近最新免费中文字幕在线| 国内久久婷婷六月综合欲色啪| 色综合婷婷激情| 欧美日本亚洲视频在线播放| 欧美一区二区亚洲| 美女高潮的动态| 国产欧美日韩精品一区二区| 亚洲中文字幕一区二区三区有码在线看| 日本 av在线| 51午夜福利影视在线观看| 狠狠狠狠99中文字幕| 日韩高清综合在线| svipshipincom国产片| 在线十欧美十亚洲十日本专区| 亚洲自拍偷在线| 久久精品影院6| 97超级碰碰碰精品色视频在线观看| 国内精品久久久久精免费| 国产黄片美女视频| 国产精品香港三级国产av潘金莲| 久久久成人免费电影| 欧美在线一区亚洲| 亚洲男人的天堂狠狠| 亚洲最大成人中文| 日韩 欧美 亚洲 中文字幕| 长腿黑丝高跟| 久久性视频一级片| 观看免费一级毛片| 国产精品综合久久久久久久免费| 又粗又爽又猛毛片免费看| 我要搜黄色片| 蜜桃亚洲精品一区二区三区| 欧美不卡视频在线免费观看| 国产精品久久久人人做人人爽| 18禁美女被吸乳视频| 国产午夜福利久久久久久| 久久这里只有精品中国| 免费观看人在逋| 九色国产91popny在线| 又紧又爽又黄一区二区| 日韩成人在线观看一区二区三区| 1024手机看黄色片| 一区二区三区国产精品乱码| 成年版毛片免费区| 亚洲乱码一区二区免费版| 香蕉丝袜av| 性欧美人与动物交配| 男女那种视频在线观看| 亚洲精品在线美女| 精品国产美女av久久久久小说| 精品一区二区三区视频在线 | avwww免费| 亚洲狠狠婷婷综合久久图片| 最新美女视频免费是黄的| 日韩中文字幕欧美一区二区| 国产欧美日韩精品亚洲av| 精品一区二区三区视频在线 | 此物有八面人人有两片| 丁香欧美五月| 人人妻人人澡欧美一区二区| 一进一出抽搐gif免费好疼| 12—13女人毛片做爰片一| 亚洲国产精品sss在线观看| 亚洲精品亚洲一区二区| 亚洲男人的天堂狠狠| aaaaa片日本免费| 亚洲成a人片在线一区二区| 一本一本综合久久| 日韩亚洲欧美综合| 一个人观看的视频www高清免费观看| 欧美zozozo另类| 国产精品日韩av在线免费观看| 亚洲,欧美精品.| 国产精品影院久久| 一区二区三区激情视频| 老司机深夜福利视频在线观看| 人妻久久中文字幕网| aaaaa片日本免费| 男人的好看免费观看在线视频| 成人午夜高清在线视频| 日韩欧美精品免费久久 | 国产亚洲精品一区二区www| 99久久综合精品五月天人人| 日韩有码中文字幕| 天天躁日日操中文字幕| 亚洲狠狠婷婷综合久久图片| 亚洲精品国产精品久久久不卡| www.色视频.com| 午夜福利在线观看免费完整高清在 | 久久久久亚洲av毛片大全| 每晚都被弄得嗷嗷叫到高潮| 一区二区三区激情视频| 手机成人av网站| 国内精品久久久久久久电影| 舔av片在线| 黄色女人牲交| 在线十欧美十亚洲十日本专区| 欧美中文综合在线视频| 99riav亚洲国产免费| 在线免费观看的www视频| 人妻夜夜爽99麻豆av| 韩国av一区二区三区四区| xxx96com| 日本三级黄在线观看| 尤物成人国产欧美一区二区三区| 欧美极品一区二区三区四区| 欧美国产日韩亚洲一区| 中文在线观看免费www的网站| 天堂影院成人在线观看| 国内精品久久久久精免费| 国产v大片淫在线免费观看| 亚洲欧美日韩高清专用| 久久精品国产亚洲av香蕉五月| 精品一区二区三区视频在线观看免费| 国产精品美女特级片免费视频播放器| 欧美日韩乱码在线| 亚洲内射少妇av| www.999成人在线观看| 久久精品国产清高在天天线| 久久久久久久久大av| www.色视频.com| 国产精品免费一区二区三区在线| 国产免费av片在线观看野外av| 中文字幕人成人乱码亚洲影| 老汉色∧v一级毛片| 国产毛片a区久久久久| 一级毛片女人18水好多| 精品福利观看| 亚洲一区高清亚洲精品| 久久久久九九精品影院| 日本精品一区二区三区蜜桃| 天堂动漫精品| 国产精品自产拍在线观看55亚洲| 观看美女的网站| 免费看十八禁软件| 男人的好看免费观看在线视频| 免费在线观看成人毛片| 色综合欧美亚洲国产小说| 国产精品 欧美亚洲| 国产成人欧美在线观看| 精品久久久久久久久久免费视频| 久久久久免费精品人妻一区二区| 国产精品久久久久久久电影 | 精品电影一区二区在线| 国产男靠女视频免费网站| 美女黄网站色视频| 一夜夜www| 亚洲av一区综合| 日本成人三级电影网站| 久久人妻av系列| 亚洲久久久久久中文字幕| 国产单亲对白刺激| 人妻久久中文字幕网| 免费在线观看成人毛片| 很黄的视频免费| 在线观看午夜福利视频| 亚洲黑人精品在线| 久9热在线精品视频| 中出人妻视频一区二区| 人妻丰满熟妇av一区二区三区| 国产高清三级在线| 校园春色视频在线观看| 欧美中文综合在线视频| 99久国产av精品| 性欧美人与动物交配| 一卡2卡三卡四卡精品乱码亚洲| 人妻久久中文字幕网| 观看免费一级毛片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产乱人视频| 久久香蕉精品热| 免费在线观看影片大全网站| 女同久久另类99精品国产91| 亚洲国产高清在线一区二区三| 99视频精品全部免费 在线| 免费在线观看影片大全网站| 99国产极品粉嫩在线观看| 久久久久久九九精品二区国产| 在线观看av片永久免费下载| 精品人妻偷拍中文字幕| 亚洲国产高清在线一区二区三| 三级国产精品欧美在线观看| 国产精品综合久久久久久久免费| 亚洲自拍偷在线| www日本黄色视频网| av欧美777| 久久久久久久久大av| 一二三四社区在线视频社区8| 国产一区二区三区在线臀色熟女| 最近最新免费中文字幕在线| 亚洲av第一区精品v没综合| 18禁美女被吸乳视频| 亚洲七黄色美女视频| 91在线精品国自产拍蜜月 | 99久久无色码亚洲精品果冻| 日韩欧美 国产精品| 一本综合久久免费| 757午夜福利合集在线观看| 麻豆成人av在线观看| 成年女人看的毛片在线观看| 国产精品影院久久| 搡老岳熟女国产| 久久久久久久久大av| 成人特级av手机在线观看| 日本a在线网址| 18美女黄网站色大片免费观看| 可以在线观看的亚洲视频| 午夜激情福利司机影院| 欧美色欧美亚洲另类二区| 亚洲国产精品久久男人天堂| 午夜精品一区二区三区免费看| 亚洲中文日韩欧美视频| 精品久久久久久久人妻蜜臀av| 国产精品乱码一区二三区的特点| 麻豆久久精品国产亚洲av| 午夜福利在线在线| 一a级毛片在线观看| 91麻豆精品激情在线观看国产| 天堂动漫精品| 狠狠狠狠99中文字幕| 欧美日韩一级在线毛片| 色视频www国产| 99久国产av精品| 国产成年人精品一区二区| 性色av乱码一区二区三区2| 青草久久国产| 亚洲欧美激情综合另类| 中文字幕人妻熟人妻熟丝袜美 | av在线蜜桃| 免费看a级黄色片| 国产精品,欧美在线| 亚洲精品国产精品久久久不卡| 日本a在线网址| 国产亚洲精品久久久com| 午夜福利成人在线免费观看| 尤物成人国产欧美一区二区三区| 国产亚洲精品av在线| 此物有八面人人有两片| 露出奶头的视频| 成人午夜高清在线视频| 两人在一起打扑克的视频| 国产伦精品一区二区三区四那| 国产乱人视频| 亚洲精品国产精品久久久不卡| a级毛片a级免费在线| 欧美黄色淫秽网站| 亚洲国产精品合色在线| 国产一区二区在线观看日韩 | 母亲3免费完整高清在线观看| 亚洲欧美激情综合另类| 51国产日韩欧美| 成人性生交大片免费视频hd| 国产淫片久久久久久久久 | 久久久久久久亚洲中文字幕 | 美女黄网站色视频| 免费人成视频x8x8入口观看| 久久久国产成人精品二区| 十八禁人妻一区二区| 少妇裸体淫交视频免费看高清| 一区二区三区免费毛片| 岛国在线免费视频观看| 99国产精品一区二区蜜桃av| 性色avwww在线观看| 亚洲精品456在线播放app | 久久人妻av系列| 怎么达到女性高潮| 欧美不卡视频在线免费观看| 99久久久亚洲精品蜜臀av| 小说图片视频综合网站| 国产男靠女视频免费网站| 搡女人真爽免费视频火全软件 | 操出白浆在线播放| 一本精品99久久精品77| 欧美精品啪啪一区二区三区| 日本五十路高清| 3wmmmm亚洲av在线观看| 九色成人免费人妻av| 禁无遮挡网站| 怎么达到女性高潮| 日日夜夜操网爽| 国产久久久一区二区三区| 一区二区三区国产精品乱码| 老司机在亚洲福利影院| 亚洲电影在线观看av| 免费看a级黄色片| 黄色女人牲交| 嫩草影院精品99| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久久末码| 国产成年人精品一区二区| 在线观看日韩欧美| 免费观看精品视频网站| 一区二区三区高清视频在线| 国产av不卡久久| 国产一区二区激情短视频| 亚洲成av人片免费观看| 国产亚洲av嫩草精品影院| 亚洲av不卡在线观看| 中文字幕av在线有码专区| 免费在线观看影片大全网站| 色吧在线观看| 欧美黑人巨大hd| 91麻豆av在线| 特级一级黄色大片| 国产精品影院久久| 久久久国产精品麻豆| 热99re8久久精品国产| 人人妻人人看人人澡| 手机成人av网站| 成人三级黄色视频| 免费观看人在逋| 国产单亲对白刺激| 琪琪午夜伦伦电影理论片6080| 欧美成狂野欧美在线观看| 中文资源天堂在线| 亚洲18禁久久av| 欧美一级毛片孕妇| 日韩中文字幕欧美一区二区| 久久久久久久久久黄片| 久久久色成人| 国产激情偷乱视频一区二区| 欧美日韩黄片免| 一进一出抽搐gif免费好疼| 黄片小视频在线播放| 午夜福利在线观看吧| 欧美黄色淫秽网站| 国产精品乱码一区二三区的特点| 日本五十路高清| 欧美性感艳星| 久99久视频精品免费| 淫妇啪啪啪对白视频| 老司机在亚洲福利影院| 嫩草影视91久久| 岛国视频午夜一区免费看| 日本一二三区视频观看| 麻豆一二三区av精品| 亚洲乱码一区二区免费版| av国产免费在线观看| 午夜亚洲福利在线播放| 欧美成人免费av一区二区三区| 99久久99久久久精品蜜桃| 女人高潮潮喷娇喘18禁视频| 少妇裸体淫交视频免费看高清| 一卡2卡三卡四卡精品乱码亚洲| 天堂av国产一区二区熟女人妻| 国产精品久久久久久久电影 | 国语自产精品视频在线第100页| 日本免费一区二区三区高清不卡| 男女床上黄色一级片免费看| 欧美成人a在线观看| 欧美极品一区二区三区四区| 国产伦在线观看视频一区| 国产午夜福利久久久久久| 久久久久久久久大av| a级毛片a级免费在线| 国产欧美日韩一区二区精品| 天堂av国产一区二区熟女人妻| 精品人妻1区二区| 十八禁网站免费在线| 亚洲成a人片在线一区二区| 人人妻人人看人人澡| 日韩精品青青久久久久久| 亚洲av美国av| 亚洲欧美激情综合另类| 搡老熟女国产l中国老女人| 女人高潮潮喷娇喘18禁视频| 少妇裸体淫交视频免费看高清| 欧美极品一区二区三区四区| 国语自产精品视频在线第100页| 久久99热这里只有精品18| 最近在线观看免费完整版| 两个人看的免费小视频| 99国产精品一区二区三区| 国产精品98久久久久久宅男小说| 男女做爰动态图高潮gif福利片| 日本一二三区视频观看| 亚洲一区高清亚洲精品| 国内精品一区二区在线观看| 国产一区二区亚洲精品在线观看| 69人妻影院| 欧美乱妇无乱码| 国产精品久久久久久久电影 | 男女那种视频在线观看| 亚洲精品美女久久久久99蜜臀| 精品一区二区三区av网在线观看| 国产不卡一卡二| 国产久久久一区二区三区| 国产主播在线观看一区二区| 国产精品免费一区二区三区在线| 成人av在线播放网站| 老汉色∧v一级毛片| 国产美女午夜福利| 欧美又色又爽又黄视频| 日韩av在线大香蕉| 嫁个100分男人电影在线观看| 91麻豆av在线| 久久久久久久亚洲中文字幕 | 国产亚洲精品一区二区www| 亚洲,欧美精品.| 久久久久性生活片| 白带黄色成豆腐渣| 国产一区二区亚洲精品在线观看| 国产又黄又爽又无遮挡在线| avwww免费| 亚洲男人的天堂狠狠| 国产一区二区在线观看日韩 | 在线播放无遮挡| 久久久精品欧美日韩精品| 午夜日韩欧美国产| 亚洲熟妇熟女久久| 日日夜夜操网爽| 老汉色av国产亚洲站长工具| 男女那种视频在线观看| 免费电影在线观看免费观看| 国产又黄又爽又无遮挡在线| 国产精品一及| 国产aⅴ精品一区二区三区波| 国产欧美日韩一区二区三| 成人18禁在线播放| 国产欧美日韩一区二区三| 午夜久久久久精精品| 91麻豆精品激情在线观看国产| 尤物成人国产欧美一区二区三区| 他把我摸到了高潮在线观看| 久久人妻av系列| 国产国拍精品亚洲av在线观看 | 成人高潮视频无遮挡免费网站| 99国产综合亚洲精品| 18+在线观看网站| 亚洲狠狠婷婷综合久久图片| 欧美乱码精品一区二区三区| 99国产综合亚洲精品| 又黄又粗又硬又大视频| 91在线精品国自产拍蜜月 | 黄色视频,在线免费观看| 国产亚洲精品av在线| 黄色视频,在线免费观看| 国产伦一二天堂av在线观看| 制服丝袜大香蕉在线| 搞女人的毛片| 久久6这里有精品| 国内精品一区二区在线观看| 99久久无色码亚洲精品果冻| 免费av不卡在线播放| 黄片大片在线免费观看| 精品乱码久久久久久99久播| 天堂影院成人在线观看| 岛国视频午夜一区免费看| 99精品在免费线老司机午夜| 综合色av麻豆| 国产伦在线观看视频一区| 久久久久久久精品吃奶| 麻豆成人av在线观看| 久久午夜亚洲精品久久| 91久久精品电影网| 美女 人体艺术 gogo| 午夜福利18| 内地一区二区视频在线| 中文字幕人成人乱码亚洲影| 老司机福利观看| 国产免费av片在线观看野外av|