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

    翼面熱靜氣動彈性的流固熱交錯迭代耦合分析*

    2022-11-04 08:37:20常斌黃杰姚衛(wèi)星
    振動、測試與診斷 2022年5期
    關(guān)鍵詞:氣動彈性翼面氣動力

    常斌,黃杰,2,姚衛(wèi)星,2

    (1.南京航空航天大學(xué)機(jī)械結(jié)構(gòu)力學(xué)及控制國家重點(diǎn)實(shí)驗(yàn)室 南京,210016)

    (2.南京航空航天大學(xué)飛行器先進(jìn)設(shè)計(jì)技術(shù)國防重點(diǎn)學(xué)科實(shí)驗(yàn)室 南京,210016)

    引言

    翼面的靜氣動彈性是飛行器設(shè)計(jì)必須考慮的問題,它涉及到氣動力與結(jié)構(gòu)彈性變形之間的耦合。隨著飛行器速度越來越快,氣動加熱引起翼面結(jié)構(gòu)溫度升高[1-3],導(dǎo)致翼面結(jié)構(gòu)剛度發(fā)生變化[4-5],并由此提出了熱靜氣動彈性的問題,熱靜氣動彈性的研究對高超聲速飛行器的設(shè)計(jì)至關(guān)重要。

    最早的高超聲速翼面熱靜氣動彈性分析方法假設(shè)翼面結(jié)構(gòu)具有均勻的溫度場分布,分析不同溫度下的翼面熱剛度,并在此熱剛度下進(jìn)行翼面的靜氣動彈性分析。隨著氣動熱分析技術(shù)研究的深入,一些學(xué)者采用以Eckert參考焓法[6]為代表的工程近似方法分析翼面氣動熱問題,將分析獲得的翼面熱流作為邊界條件,并通過有限差分法等求解熱傳導(dǎo)方程獲得翼面溫度場分布,再研究翼面的熱剛度和熱靜氣動彈性[7-8]。以上傳統(tǒng)翼面熱靜氣動彈性分析方法的計(jì)算精度較差,已無法滿足現(xiàn)代及未來高超聲速飛行器熱靜氣動彈性的分析要求。傳統(tǒng)方法最主要的缺點(diǎn)是未考慮結(jié)構(gòu)溫度場對氣動熱的反饋效應(yīng)。實(shí)際上氣動熱與結(jié)構(gòu)溫度場之間存在著強(qiáng)烈的耦合效應(yīng),但傳統(tǒng)分析方法僅考慮了氣動熱引起的結(jié)構(gòu)溫度場升高,而未考慮翼面溫度升高后對壁面熱流的影響。

    隨著計(jì)算流體力學(xué)(computational fluid dynamics,簡稱CFD)、數(shù)值傳熱學(xué)(numerical heat transfer,簡稱NHT)及計(jì)算結(jié)構(gòu)動力學(xué)(computational structural dynamics,簡稱CSD)的發(fā)展,近年來一些學(xué)者通過CFD和NHT之間的耦合分析求解結(jié)構(gòu)溫度場,研究熱環(huán)境對結(jié)構(gòu)熱剛度的影響,最后通過CFD和CSD之間的耦合分析求解熱靜氣動彈性響應(yīng)[9-10]。該方法分2步計(jì)算,能獲得精確的壁面熱流和結(jié)構(gòu)溫度場,并準(zhǔn)確評價(jià)翼面結(jié)構(gòu)熱剛度,具有計(jì)算精度高的特點(diǎn)。由于CFD與NHT的耦合以及CFD與CSD的耦合需要進(jìn)行2次CFD計(jì)算,故該方法還存在計(jì)算效率較低的缺點(diǎn)。此外該方法也無法考慮翼面變形對氣動熱帶來的影響,因此需要一種兼顧計(jì)算精度和計(jì)算效率的耦合分析方法。

    筆者提出了一種針對高超聲速翼面熱靜氣動彈性的流固熱交錯迭代耦合方法,其主要特點(diǎn)是能考慮翼面變形對氣動熱的影響,且不需要像傳統(tǒng)熱靜氣動彈性耦合方法分2步求解,并且利用該方法進(jìn)行了典型高超聲速翼面熱靜氣動彈性問題的分析。

    1 流體和結(jié)構(gòu)傳熱控制方程

    流體動力學(xué)的控制方程為Navier-Stokes方程,

    其在直角坐標(biāo)系下的積分形式為

    其中:W為守恒向量;Fc為 對 流 通量;Fv為黏性通量;dS為控制體V的邊界面;n為邊界面dS的 外法線單位向量。

    Navier-Stokes方程中對流通量Fc的空間離散采用AUSM+格式[11],黏性通量Fv的空間離散采用中心差分格 式,湍 流模擬采用Menter′s SSTk-ω兩方程模型[12],時間推進(jìn)采用LU-SGS格式[13],且采用聚合多重網(wǎng)格法[14]加速收斂。

    翼面結(jié)構(gòu)傳熱分析的控制方程為

    其中:kx,ky和kz分別為材料在x,y和z方向的導(dǎo)熱系數(shù);ρ和c分別為結(jié)構(gòu)材料的密度和比熱容。

    翼面結(jié)構(gòu)傳熱分析的邊界條件包括壁面熱流Qw和壁面輻射熱流Qr,其表達(dá)式為

    其中:ε為表面發(fā)射率;σ為玻爾茲曼常數(shù);Tw為壁面溫度。

    筆者采用有限元法進(jìn)行翼面結(jié)構(gòu)傳熱及變形的分析。

    2 交錯迭代耦合分析方法

    高超聲速飛行器做定常飛行時,翼面所受的氣動力和氣動熱均穩(wěn)定,其可通過數(shù)值方法計(jì)算同時得到,因此統(tǒng)稱為氣動環(huán)境。翼面熱靜氣動彈性分析中主要存在2種耦合效應(yīng):①氣動環(huán)境中的氣動熱與結(jié)構(gòu)傳熱之間的耦合效應(yīng),當(dāng)氣動熱作用于翼面時,翼面結(jié)構(gòu)溫度升高,翼面結(jié)構(gòu)溫度的升高導(dǎo)致邊界層內(nèi)氣體與壁面溫度的梯度減小,即氣動熱效應(yīng)將減弱;②氣動環(huán)境中的氣動力與結(jié)構(gòu)變形之間的耦合效應(yīng),當(dāng)氣動力作用于翼面時,翼面產(chǎn)生變形,變形會造成氣動力重分布。此外,翼面變形也會造成壁面熱流的重分布,且熱環(huán)境還會影響翼面結(jié)構(gòu)剛度(材料剛度和熱應(yīng)力引起的附加幾何剛度),進(jìn)而影響結(jié)構(gòu)變形。筆者考慮以上所述的耦合效應(yīng),而忽略其他次要因素,建立了如圖1所示的翼面熱靜氣動彈性分析耦合模型。

    圖1 耦合模型Fig.1 Coupled model

    根據(jù)圖1中的耦合模型,筆者采用流固熱交錯迭代耦合方法進(jìn)行高超聲速翼面熱靜氣動彈性的研究,交錯迭代耦合方法如圖2所示,圖中的i為耦合迭代步。采用CFD方法進(jìn)行翼面氣動力和氣動熱的計(jì)算,采用NHT方法進(jìn)行翼面結(jié)構(gòu)溫度場的計(jì)算,采用CSD方法進(jìn)行翼面結(jié)構(gòu)變形的計(jì)算,其中CFD方法基于有限體積法,而NHT和CSD方法均基于有限元法。流固熱交錯迭代耦合方法的基本假設(shè)和特點(diǎn)如下:

    圖2 交錯迭代耦合方法Fig.2 Staggered iterative coupled method

    1)在流場計(jì)算過程中(氣動力和氣動熱)翼面結(jié)構(gòu)參數(shù)(結(jié)構(gòu)變形和壁面溫度)保持不變,同理在結(jié)構(gòu)計(jì)算過程中(結(jié)構(gòu)變形和壁面溫度)翼面氣動環(huán)境(氣動力和氣動熱)保持不變;

    2)只需進(jìn)行一次耦合分析即可獲得翼面氣動力、壁面熱流、溫度場和結(jié)構(gòu)變形的全部穩(wěn)態(tài)結(jié)果。

    翼面熱靜氣動彈性的交錯迭代耦合方法分析流程如圖3所示,其主要步驟為:

    圖3 交錯迭代耦合分析流程Fig.3 Staggered iterative coupled analysis process

    1)首先建立CFD,NHT和CSD數(shù)值分析模型,通過CFD分析獲得初始?xì)鈩恿捅诿鏌崃鳎?/p>

    2)進(jìn)行第i步的CFD計(jì)算,將計(jì)算獲得的氣動力Fi和壁面熱流Qi分別傳遞給CSD模型和NHT模型;

    3)進(jìn)行第i步的NHT計(jì)算,將計(jì)算獲得的翼面結(jié)構(gòu)溫度場Ti傳遞給CSD模型,并更新翼面結(jié)構(gòu)剛度;

    4)進(jìn)行第i步的CSD計(jì)算,將獲得的翼面表面位移Uwi及表面溫度Twi傳遞給CFD模型,并進(jìn)行流體網(wǎng)格變形;

    5)若翼面結(jié)構(gòu)變形和溫度場均收斂,則結(jié)束迭代。否則重復(fù)步驟2~4,直至結(jié)果收斂,結(jié)束分析。

    由于流場和結(jié)構(gòu)的網(wǎng)格尺寸相差較大,耦合變量需要在耦合面上通過插值算法實(shí)現(xiàn)數(shù)據(jù)的傳遞。筆者采用文獻(xiàn)[15]中的基于控制面的雙向映射插值方法進(jìn)行壁面熱流、壁面溫度、氣動力以及翼面變形的數(shù)據(jù)傳遞。此外為了加速收斂,在流場計(jì)算中引入了子迭代。

    目前還沒有熱靜氣動彈性的標(biāo)準(zhǔn)試驗(yàn),且翼面熱靜氣動彈性分析的核心是精確計(jì)算氣動熱和結(jié)構(gòu)溫度場,故筆者采用NASA的高超聲速空心圓管模型[16]進(jìn)行氣動熱和結(jié)構(gòu)溫度場計(jì)算精度的驗(yàn)證。文獻(xiàn)[17]已經(jīng)進(jìn)行了空心圓管模型的耦合計(jì)算,外壁面熱流和溫度分布的計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合良好,驗(yàn)證了計(jì)算精度。

    3 翼面熱靜氣動彈性分析

    3.1 計(jì)算模型

    筆者采用小展弦比翼面進(jìn)行熱靜氣動彈性的分析,其平面和剖面如圖4所示。小展弦比翼面來源于F104戰(zhàn)斗機(jī)機(jī)翼,圖中參數(shù)C為翼根弦長。假設(shè)高超聲速自由來流馬赫數(shù)從6增加至9,飛行高度為10 km,攻角為0.5°,且整個翼面結(jié)構(gòu)的初始溫度為300 K。劃分了流場的結(jié)構(gòu)網(wǎng)格,網(wǎng)格量約為80萬,為了保證翼面熱流的計(jì)算精度,壁面第1層網(wǎng)格高度小于1×10-5m,圖5為CFD計(jì)算網(wǎng)格及結(jié)構(gòu)傳熱和變形分析的翼面結(jié)構(gòu)表面計(jì)算網(wǎng)格。筆者假設(shè)翼面前緣為耐高溫碳/碳復(fù)合材料結(jié)構(gòu),其導(dǎo)熱系數(shù)、彈性模量和熱膨脹系數(shù)分別為42 W/(m·K),95 GPa和4×10-6K-1。翼面其他部位采用TA7耐高溫鈦合金結(jié)構(gòu),材料密度為4.4×103kg/m3,但材料導(dǎo)熱系數(shù)、彈性模量和熱膨脹系數(shù)均與溫度相關(guān),可參考文獻(xiàn)[15]。文獻(xiàn)[15]進(jìn)行了翼面的CFD和NHT耦合計(jì)算,獲得了結(jié)構(gòu)溫度場,并研究了在熱環(huán)境下彈性模量和熱應(yīng)力對模態(tài)的影響。筆者考慮了氣動熱、氣動力、結(jié)構(gòu)傳熱和結(jié)構(gòu)變形,是典型的熱靜氣動彈性問題。翼面外表面有高輻射率的涂層,其發(fā)射率為0.85。此外流場計(jì)算中引入了子迭代,子迭代步數(shù)設(shè)置為30。

    圖4 小展弦比翼面的平面和剖面Fig.4 Platform and cross-sectional views of the low aspect ratio wing

    圖5 計(jì)算網(wǎng)格Fig.5 Computational grids

    3.2 翼面熱環(huán)境及熱靜氣動彈性變形

    筆者通過流固熱交錯迭代耦合方法獲得了馬赫數(shù)從6增加至9的計(jì)算結(jié)果。翼面結(jié)構(gòu)溫度場和翼面變形均在25步以內(nèi)達(dá)到收斂狀態(tài)。圖6給出了其表面溫度云圖,結(jié)果表明在不同馬赫數(shù)下計(jì)算獲得的翼面結(jié)構(gòu)溫度分布情況相似。由于翼面前緣當(dāng)量半徑很小,故該部位的壁面熱流和溫度最高,且往下游翼面結(jié)構(gòu)溫度逐漸降低。此外馬赫數(shù)從6增加到9,翼面結(jié)構(gòu)最高溫度從1 525.4 K上升到2 163.3 K,翼面結(jié)構(gòu)最低溫度從824.9 K上升到996.2 K。

    圖6 不同馬赫數(shù)下翼面結(jié)構(gòu)溫度場Fig.6 Temperature field of wing for different Mach number

    熱靜氣動彈性與常規(guī)靜氣動彈性的主要區(qū)別是熱環(huán)境下翼面結(jié)構(gòu)剛度會發(fā)生變化。首先,高溫會直接影響材料的彈性模量,針對本研究采用的TA7耐高溫鈦合金材料,溫度越高材料的彈性模量越低;其次,高溫下翼面結(jié)構(gòu)會膨脹,結(jié)構(gòu)受到約束及內(nèi)部溫度分布不均勻時均會產(chǎn)生熱應(yīng)力。熱應(yīng)力的本質(zhì)是預(yù)應(yīng)力,預(yù)應(yīng)力會造成幾何非線性現(xiàn)象,即翼面結(jié)構(gòu)會產(chǎn)生附加的幾何剛度。因此,必須綜合考慮熱環(huán)境造成的材料彈性模量的降低及熱應(yīng)力引起的附加幾何剛度才能準(zhǔn)確計(jì)算翼面變形。

    筆者通過流固熱交錯迭代耦合方法分析獲得了馬赫數(shù)從6增加至9的翼面法向(垂直于翼面方向)的熱靜氣動彈性變形,其熱環(huán)境下翼面法向位移如圖7所示,結(jié)果表明在不同馬赫數(shù)下翼面熱靜氣動彈性變形情況相似。圖8為300 K和熱環(huán)境下翼面最大厚度處的法向位移沿展向的分布情況,圖9為300 K和熱環(huán)境下翼面法向最大位移隨馬赫數(shù)的變化情況。

    圖7 熱環(huán)境下翼面法向位移Fig.7 Normal displacements of wing under thermal environments

    圖8 翼面法向位移沿展向分布情況Fig.8 Normal displacements of wing along the spanwise direction

    圖9 翼面最大法向位移Fig.9 Maximum normal displacements of wing

    由以上分析結(jié)果可知,隨著馬赫數(shù)的增加,翼面變形越大。此外在相同馬赫數(shù)下,考慮熱效應(yīng)時翼面法向位移明顯大于300 K時的分析結(jié)果。馬赫數(shù)為6,7,8和9時,熱環(huán)境下翼面法向最大位移分別比300 K時大41.05%,51.10%,57.64%和63.28%,即馬赫數(shù)越大,考慮熱環(huán)境和不考慮熱環(huán)境獲得的翼面變形量之差越大。這是由于熱環(huán)境下翼面結(jié)構(gòu)剛度降低了,且馬赫數(shù)越大翼面結(jié)構(gòu)剛度下降越快,造成變形量迅速增加。

    綜上所述,進(jìn)行高超聲速翼面靜氣動彈性分析時必須考慮其熱效應(yīng),熱環(huán)境導(dǎo)致翼面結(jié)構(gòu)剛度降低,造成翼面變形高于常溫下的結(jié)果,應(yīng)給予足夠重視。

    4 結(jié)論

    1)提出了一種針對高超聲速翼面熱靜氣動彈性的流固熱交錯迭代數(shù)值耦合方法。該方法充分考慮了氣動環(huán)境(氣動力和氣動熱)與結(jié)構(gòu)變形之間的耦合、氣動熱與結(jié)構(gòu)溫度場之間的耦合以及溫度場對結(jié)構(gòu)剛度的影響,能克服傳統(tǒng)熱靜氣動彈性耦合方法中未考慮翼面變形對氣動熱影響的缺點(diǎn)。

    2)進(jìn)行了翼面熱靜氣動彈性分析,針對本研究的翼面結(jié)構(gòu)形式和來流狀態(tài),熱環(huán)境造成了翼面結(jié)構(gòu)剛度下降,從而導(dǎo)致熱環(huán)境下翼面法向最大位移比300 K時增加了40%以上,且馬赫數(shù)越大,兩者之差越大。因此,熱環(huán)境對高超聲速翼面的靜氣動彈性變形影響明顯,在分析時必須考慮其影響才能獲得準(zhǔn)確的分析結(jié)果。

    猜你喜歡
    氣動彈性翼面氣動力
    基于拔銷器鎖定的飛行器氣動控制面解鎖控制方法
    飛行載荷外部氣動力的二次規(guī)劃等效映射方法
    固定翼二維彈道修正引信升力翼面位置的影響
    基于Aerobook平臺的復(fù)合材料翼面結(jié)構(gòu)設(shè)計(jì)流程
    側(cè)風(fēng)對拍動翅氣動力的影響
    飛翼無人機(jī)嗡鳴氣動彈性響應(yīng)分析
    模態(tài)選取對靜氣動彈性分析的影響
    直升機(jī)的氣動彈性問題
    大型風(fēng)力機(jī)整機(jī)氣動彈性響應(yīng)計(jì)算
    高速鐵路接觸線覆冰后氣動力特性的風(fēng)洞試驗(yàn)研究
    搡女人真爽免费视频火全软件 | 欧美精品啪啪一区二区三区| netflix在线观看网站| 日韩在线高清观看一区二区三区 | 久久久久国产精品人妻aⅴ院| 少妇高潮的动态图| 波多野结衣巨乳人妻| 人妻丰满熟妇av一区二区三区| 中文字幕精品亚洲无线码一区| 中国美女看黄片| a级毛片a级免费在线| 美女黄网站色视频| 国产精品人妻久久久影院| 国产午夜精品久久久久久一区二区三区 | 欧美一区二区亚洲| 天堂动漫精品| 高清日韩中文字幕在线| 美女黄网站色视频| 在线免费观看的www视频| 他把我摸到了高潮在线观看| 国产黄a三级三级三级人| 我要看日韩黄色一级片| 日韩强制内射视频| 久久精品国产鲁丝片午夜精品 | avwww免费| 五月玫瑰六月丁香| 搡老妇女老女人老熟妇| 国产亚洲精品久久久com| 国产日本99.免费观看| 成人亚洲精品av一区二区| 久久国产精品人妻蜜桃| 免费人成在线观看视频色| 中文字幕av在线有码专区| 国产aⅴ精品一区二区三区波| 乱码一卡2卡4卡精品| 国产精品一区二区免费欧美| 亚洲三级黄色毛片| 国产精品野战在线观看| 国产 一区 欧美 日韩| 婷婷色综合大香蕉| 一个人看的www免费观看视频| 亚洲人成网站在线播放欧美日韩| 精品国产三级普通话版| 久久6这里有精品| 成年女人毛片免费观看观看9| 韩国av一区二区三区四区| 欧美又色又爽又黄视频| 国产精品精品国产色婷婷| 亚洲人成伊人成综合网2020| 人人妻人人看人人澡| 亚洲 国产 在线| 久久久久久国产a免费观看| 亚洲国产色片| 中国美白少妇内射xxxbb| 22中文网久久字幕| 欧美日韩中文字幕国产精品一区二区三区| 99热这里只有是精品在线观看| 日韩在线高清观看一区二区三区 | 国产亚洲av嫩草精品影院| 亚洲精品乱码久久久v下载方式| 成年女人毛片免费观看观看9| 国产真实伦视频高清在线观看 | 免费高清视频大片| 有码 亚洲区| 国产成人福利小说| 欧美精品国产亚洲| 亚洲无线在线观看| 欧美极品一区二区三区四区| 中文字幕免费在线视频6| 在线播放国产精品三级| 99久久成人亚洲精品观看| 亚洲在线自拍视频| 波多野结衣巨乳人妻| 色5月婷婷丁香| 精品久久久久久,| 校园人妻丝袜中文字幕| 欧美日韩亚洲国产一区二区在线观看| 久久香蕉精品热| 国产私拍福利视频在线观看| 男人的好看免费观看在线视频| 国产综合懂色| 久久热精品热| 少妇裸体淫交视频免费看高清| 男人狂女人下面高潮的视频| 久久久久国产精品人妻aⅴ院| 成人国产麻豆网| 99久久无色码亚洲精品果冻| 国内久久婷婷六月综合欲色啪| 99视频精品全部免费 在线| 一进一出好大好爽视频| 日韩亚洲欧美综合| 亚洲中文日韩欧美视频| 国产午夜精品久久久久久一区二区三区 | 97超视频在线观看视频| 好男人在线观看高清免费视频| 免费人成在线观看视频色| 大型黄色视频在线免费观看| 国产精品精品国产色婷婷| 国产精品伦人一区二区| 99久久九九国产精品国产免费| 搡女人真爽免费视频火全软件 | 非洲黑人性xxxx精品又粗又长| 欧美高清成人免费视频www| 亚洲真实伦在线观看| 亚洲av一区综合| 亚洲久久久久久中文字幕| 日本成人三级电影网站| videossex国产| 欧美极品一区二区三区四区| 国产视频一区二区在线看| 精品午夜福利在线看| 久久久国产成人免费| 赤兔流量卡办理| 国产精品无大码| 免费人成视频x8x8入口观看| 国产av一区在线观看免费| 亚洲电影在线观看av| 国产一区二区在线观看日韩| 中文字幕熟女人妻在线| 午夜福利成人在线免费观看| 国产精品亚洲一级av第二区| 国产av不卡久久| 日本撒尿小便嘘嘘汇集6| 伊人久久精品亚洲午夜| 日日撸夜夜添| 色视频www国产| 99久久九九国产精品国产免费| 久久久久久伊人网av| 国产精品嫩草影院av在线观看 | 国产精品一区二区三区四区久久| 热99re8久久精品国产| 特级一级黄色大片| 国产精品久久久久久精品电影| 国产精华一区二区三区| 亚洲午夜理论影院| 久久亚洲精品不卡| 日韩亚洲欧美综合| 久久久久久久久久成人| 亚洲综合色惰| 久久中文看片网| 国产精品久久视频播放| 99精品在免费线老司机午夜| av在线观看视频网站免费| 最近中文字幕高清免费大全6 | 日韩欧美精品v在线| 最新中文字幕久久久久| 在现免费观看毛片| 18+在线观看网站| 日韩中字成人| 成人av在线播放网站| 欧美日韩国产亚洲二区| 国产91精品成人一区二区三区| 麻豆国产av国片精品| 男女那种视频在线观看| 欧美日韩黄片免| 国产精品爽爽va在线观看网站| 国产精品伦人一区二区| 日本成人三级电影网站| 国产一区二区三区在线臀色熟女| 午夜福利成人在线免费观看| 夜夜看夜夜爽夜夜摸| 亚洲国产欧美人成| 亚洲欧美精品综合久久99| 99视频精品全部免费 在线| 久久人人爽人人爽人人片va| 夜夜夜夜夜久久久久| 婷婷精品国产亚洲av| 国产欧美日韩精品亚洲av| 两人在一起打扑克的视频| 免费无遮挡裸体视频| 88av欧美| 久久精品国产99精品国产亚洲性色| 三级男女做爰猛烈吃奶摸视频| 久久精品影院6| 麻豆av噜噜一区二区三区| 亚洲色图av天堂| 国产欧美日韩精品亚洲av| 免费看av在线观看网站| 99视频精品全部免费 在线| 国国产精品蜜臀av免费| 成熟少妇高潮喷水视频| 国产精品一区二区性色av| 欧美性感艳星| 精品国内亚洲2022精品成人| 成人三级黄色视频| 真人一进一出gif抽搐免费| 色5月婷婷丁香| 高清日韩中文字幕在线| 欧美3d第一页| 国产老妇女一区| a级毛片免费高清观看在线播放| 中文字幕免费在线视频6| 我的老师免费观看完整版| 免费观看人在逋| 国产又黄又爽又无遮挡在线| 国产三级在线视频| 美女cb高潮喷水在线观看| 日韩在线高清观看一区二区三区 | 成年女人毛片免费观看观看9| 夜夜看夜夜爽夜夜摸| 18禁在线播放成人免费| 午夜久久久久精精品| av天堂中文字幕网| 3wmmmm亚洲av在线观看| 久久久久久国产a免费观看| 国产中年淑女户外野战色| 国产在线男女| 一进一出抽搐gif免费好疼| 免费电影在线观看免费观看| 久久99热这里只有精品18| 久久久久久久久大av| 啪啪无遮挡十八禁网站| 国产在视频线在精品| 中文字幕av成人在线电影| 免费搜索国产男女视频| 日韩在线高清观看一区二区三区 | 91精品国产九色| 国产一区二区三区视频了| 观看美女的网站| 淫妇啪啪啪对白视频| 美女大奶头视频| 最近在线观看免费完整版| 免费看光身美女| 亚洲av成人av| 国内久久婷婷六月综合欲色啪| 日韩精品青青久久久久久| 日本 欧美在线| 亚洲精华国产精华液的使用体验 | 日韩中字成人| 久久99热6这里只有精品| 久久香蕉精品热| 少妇裸体淫交视频免费看高清| 精品久久久噜噜| 亚洲精品日韩av片在线观看| 99精品久久久久人妻精品| 成年版毛片免费区| 18禁在线播放成人免费| 可以在线观看毛片的网站| 少妇丰满av| 午夜激情福利司机影院| 全区人妻精品视频| 亚洲av熟女| 99精品在免费线老司机午夜| 国产免费av片在线观看野外av| xxxwww97欧美| 亚洲精品亚洲一区二区| 男插女下体视频免费在线播放| 国产伦一二天堂av在线观看| 88av欧美| 乱码一卡2卡4卡精品| 亚洲欧美日韩无卡精品| 精品人妻熟女av久视频| 夜夜看夜夜爽夜夜摸| 亚洲七黄色美女视频| 国产伦精品一区二区三区四那| 国产蜜桃级精品一区二区三区| 最新在线观看一区二区三区| 精品一区二区三区人妻视频| 久久亚洲真实| 欧美激情久久久久久爽电影| 91久久精品国产一区二区三区| 一本精品99久久精品77| 久久精品久久久久久噜噜老黄 | 久久精品国产亚洲av香蕉五月| 两人在一起打扑克的视频| 国产一级毛片七仙女欲春2| 狠狠狠狠99中文字幕| 国产激情偷乱视频一区二区| 一进一出好大好爽视频| 久久精品国产99精品国产亚洲性色| 老熟妇仑乱视频hdxx| 日韩欧美精品免费久久| 欧美不卡视频在线免费观看| 超碰av人人做人人爽久久| 国产欧美日韩精品一区二区| 久久精品国产清高在天天线| 女的被弄到高潮叫床怎么办 | 日本在线视频免费播放| 亚洲成人中文字幕在线播放| 三级毛片av免费| 成人特级黄色片久久久久久久| 亚洲av熟女| 色5月婷婷丁香| 国产成人福利小说| 黄色欧美视频在线观看| 国产精品一及| 中文字幕人妻熟人妻熟丝袜美| 深夜a级毛片| 中文字幕精品亚洲无线码一区| 午夜福利在线观看免费完整高清在 | 久久国产乱子免费精品| 一进一出抽搐动态| 一区二区三区高清视频在线| 高清毛片免费观看视频网站| 国内精品一区二区在线观看| 少妇裸体淫交视频免费看高清| 看片在线看免费视频| 亚洲自偷自拍三级| 精品人妻一区二区三区麻豆 | 亚洲18禁久久av| 别揉我奶头 嗯啊视频| 无遮挡黄片免费观看| 狂野欧美白嫩少妇大欣赏| 日本黄大片高清| 99riav亚洲国产免费| 国产精品国产高清国产av| 美女xxoo啪啪120秒动态图| 成人综合一区亚洲| 国产男人的电影天堂91| 久久久久精品国产欧美久久久| 欧美bdsm另类| 精品日产1卡2卡| 国产av在哪里看| 国产精品一区www在线观看 | 亚洲一区高清亚洲精品| 嫩草影视91久久| 午夜精品久久久久久毛片777| 全区人妻精品视频| 久久久久九九精品影院| 九九在线视频观看精品| 久久午夜福利片| 又黄又爽又免费观看的视频| 精品免费久久久久久久清纯| 色5月婷婷丁香| 国产精品爽爽va在线观看网站| 窝窝影院91人妻| 亚洲av免费在线观看| 成人特级av手机在线观看| 蜜桃亚洲精品一区二区三区| 欧美日本亚洲视频在线播放| aaaaa片日本免费| 国产精品综合久久久久久久免费| 在现免费观看毛片| 免费在线观看影片大全网站| 毛片女人毛片| 国产免费一级a男人的天堂| 美女cb高潮喷水在线观看| 小说图片视频综合网站| videossex国产| 国语自产精品视频在线第100页| 搡老妇女老女人老熟妇| 欧美日韩瑟瑟在线播放| 精品免费久久久久久久清纯| 内射极品少妇av片p| 俄罗斯特黄特色一大片| 国产精品99久久久久久久久| 国内毛片毛片毛片毛片毛片| 国产精品不卡视频一区二区| 淫秽高清视频在线观看| 亚洲国产日韩欧美精品在线观看| 狂野欧美激情性xxxx在线观看| 成人欧美大片| 久久久成人免费电影| 日韩亚洲欧美综合| 欧美日韩瑟瑟在线播放| 久久人人爽人人爽人人片va| 两个人视频免费观看高清| 国产精品久久久久久久久免| 成年女人看的毛片在线观看| 男人的好看免费观看在线视频| 国产老妇女一区| 中国美女看黄片| 又紧又爽又黄一区二区| 俺也久久电影网| 在线观看免费视频日本深夜| 嫁个100分男人电影在线观看| av天堂中文字幕网| 禁无遮挡网站| 色哟哟哟哟哟哟| 国产欧美日韩一区二区精品| 国产精品综合久久久久久久免费| 国产成人aa在线观看| 欧美一区二区精品小视频在线| 亚洲国产色片| 神马国产精品三级电影在线观看| 精品人妻熟女av久视频| 日日啪夜夜撸| 国产精品嫩草影院av在线观看 | 亚洲精品一区av在线观看| 一级毛片久久久久久久久女| 免费一级毛片在线播放高清视频| 亚洲av.av天堂| 在线免费观看不下载黄p国产 | 他把我摸到了高潮在线观看| 久久99热6这里只有精品| 国产精品人妻久久久久久| 国内揄拍国产精品人妻在线| 老司机福利观看| 夜夜看夜夜爽夜夜摸| 国产精品三级大全| 蜜桃久久精品国产亚洲av| 国产午夜福利久久久久久| 久久久久久久午夜电影| 国产aⅴ精品一区二区三区波| 少妇人妻一区二区三区视频| 国产精品自产拍在线观看55亚洲| 我要搜黄色片| 亚洲图色成人| 国产一区二区三区av在线 | 91久久精品电影网| 亚洲aⅴ乱码一区二区在线播放| 男女那种视频在线观看| 神马国产精品三级电影在线观看| 久久天躁狠狠躁夜夜2o2o| 小蜜桃在线观看免费完整版高清| 黄片wwwwww| bbb黄色大片| 内射极品少妇av片p| 国产精品人妻久久久影院| 欧美成人免费av一区二区三区| 麻豆久久精品国产亚洲av| 国产精品精品国产色婷婷| 成人美女网站在线观看视频| 精品久久久久久,| 亚洲欧美清纯卡通| 黄色配什么色好看| 久久久久久大精品| 国产在线精品亚洲第一网站| av国产免费在线观看| 啦啦啦啦在线视频资源| 久久人人精品亚洲av| 久久精品国产亚洲网站| 久久久久久九九精品二区国产| 亚洲精品久久国产高清桃花| 欧美精品啪啪一区二区三区| 噜噜噜噜噜久久久久久91| 国产精品永久免费网站| 日韩人妻高清精品专区| 97热精品久久久久久| 内射极品少妇av片p| 久久精品综合一区二区三区| 亚洲天堂国产精品一区在线| 精品久久久久久久久亚洲 | 久久久午夜欧美精品| 禁无遮挡网站| 又黄又爽又免费观看的视频| 亚洲乱码一区二区免费版| 男女做爰动态图高潮gif福利片| 国产成人影院久久av| 不卡视频在线观看欧美| 国产一区二区在线av高清观看| 最近最新中文字幕大全电影3| 精品久久久久久久久亚洲 | 人妻夜夜爽99麻豆av| 欧美黑人巨大hd| 熟妇人妻久久中文字幕3abv| 91麻豆精品激情在线观看国产| 亚洲美女视频黄频| 如何舔出高潮| 91av网一区二区| 男人狂女人下面高潮的视频| 国产真实伦视频高清在线观看 | 动漫黄色视频在线观看| a在线观看视频网站| 变态另类成人亚洲欧美熟女| 亚洲av第一区精品v没综合| a级一级毛片免费在线观看| 国产成人福利小说| 欧美一级a爱片免费观看看| 色吧在线观看| 又爽又黄a免费视频| 日本 av在线| 亚洲avbb在线观看| 成人三级黄色视频| 国产69精品久久久久777片| 日日啪夜夜撸| 99热6这里只有精品| 91久久精品电影网| 欧美色视频一区免费| 永久网站在线| 亚洲欧美清纯卡通| 欧美区成人在线视频| 天天躁日日操中文字幕| 99久久九九国产精品国产免费| 国产高清有码在线观看视频| 成人高潮视频无遮挡免费网站| 国产精品av视频在线免费观看| 99国产精品一区二区蜜桃av| 91狼人影院| 亚洲最大成人av| a级一级毛片免费在线观看| 亚洲av美国av| 欧美在线一区亚洲| 国产中年淑女户外野战色| 国产激情偷乱视频一区二区| 十八禁网站免费在线| 观看免费一级毛片| 国产伦精品一区二区三区视频9| 人妻久久中文字幕网| 国产 一区精品| 欧美丝袜亚洲另类 | 91在线精品国自产拍蜜月| 欧美绝顶高潮抽搐喷水| 免费观看人在逋| 成人永久免费在线观看视频| 国产一区二区三区在线臀色熟女| 免费人成在线观看视频色| 91麻豆精品激情在线观看国产| 亚洲人成网站在线播| 午夜福利欧美成人| 欧美成人免费av一区二区三区| 搡老岳熟女国产| 91精品国产九色| www日本黄色视频网| 成熟少妇高潮喷水视频| 久久精品国产自在天天线| 欧美极品一区二区三区四区| 黄片wwwwww| 一卡2卡三卡四卡精品乱码亚洲| 精华霜和精华液先用哪个| 免费看日本二区| 婷婷色综合大香蕉| 一卡2卡三卡四卡精品乱码亚洲| 男女下面进入的视频免费午夜| 欧美高清性xxxxhd video| 18禁黄网站禁片免费观看直播| 中出人妻视频一区二区| www日本黄色视频网| 蜜桃久久精品国产亚洲av| АⅤ资源中文在线天堂| 欧美日韩瑟瑟在线播放| 99精品久久久久人妻精品| 亚洲精品成人久久久久久| 午夜精品一区二区三区免费看| 美女xxoo啪啪120秒动态图| 久久香蕉精品热| 亚洲成a人片在线一区二区| 香蕉av资源在线| 日韩在线高清观看一区二区三区 | 老司机深夜福利视频在线观看| 日韩高清综合在线| 午夜免费成人在线视频| 中文字幕人妻熟人妻熟丝袜美| 久久久精品欧美日韩精品| 欧美区成人在线视频| 俺也久久电影网| 日本-黄色视频高清免费观看| 九九久久精品国产亚洲av麻豆| 一区二区三区免费毛片| 女同久久另类99精品国产91| 亚洲av免费高清在线观看| 欧美一区二区精品小视频在线| 精品欧美国产一区二区三| 亚洲国产高清在线一区二区三| 日本爱情动作片www.在线观看 | 免费黄网站久久成人精品| 在线观看一区二区三区| 亚洲男人的天堂狠狠| 一进一出好大好爽视频| 国产色爽女视频免费观看| 国产av不卡久久| 看黄色毛片网站| 一级毛片久久久久久久久女| 黄色配什么色好看| 一区二区三区免费毛片| 精品一区二区三区视频在线观看免费| 一本一本综合久久| 成人国产一区最新在线观看| 女人被狂操c到高潮| 最新中文字幕久久久久| 级片在线观看| 亚洲av成人av| 国产精品人妻久久久久久| 一个人免费在线观看电影| 久久久久国产精品人妻aⅴ院| 嫩草影院入口| 国产精品久久久久久久电影| 国产单亲对白刺激| 亚洲精品粉嫩美女一区| 亚洲经典国产精华液单| 亚洲人成伊人成综合网2020| 97超视频在线观看视频| 人妻制服诱惑在线中文字幕| 亚洲成人久久爱视频| 丰满的人妻完整版| 久久久久久久久久成人| 午夜精品久久久久久毛片777| 搡女人真爽免费视频火全软件 | 一级黄色大片毛片| 亚洲一区高清亚洲精品| 九色成人免费人妻av| 国产精品1区2区在线观看.| 国产视频内射| 免费高清视频大片| 老熟妇乱子伦视频在线观看| 国内精品久久久久久久电影| 美女高潮喷水抽搐中文字幕| 国产精品国产高清国产av| 中文字幕精品亚洲无线码一区| 色吧在线观看| 黄色配什么色好看| 亚洲人成伊人成综合网2020| 亚洲自偷自拍三级| 精品福利观看| 久久午夜亚洲精品久久| 五月伊人婷婷丁香| 中文字幕高清在线视频| 18禁裸乳无遮挡免费网站照片| 97人妻精品一区二区三区麻豆| 亚洲av中文字字幕乱码综合| 一本久久中文字幕| 国产不卡一卡二| 18+在线观看网站| 国内精品宾馆在线| 男女下面进入的视频免费午夜| 91在线观看av| 男人和女人高潮做爰伦理| 午夜a级毛片| 久久久久免费精品人妻一区二区| 真人一进一出gif抽搐免费| 久久精品国产亚洲av香蕉五月| 成人永久免费在线观看视频| 久99久视频精品免费| 国产精品免费一区二区三区在线| 日韩精品中文字幕看吧|