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

    火箭發(fā)動(dòng)機(jī)噴管擺動(dòng)對(duì)側(cè)向載荷的影響

    2016-11-03 05:45:03劉華偉葉正寅
    固體火箭技術(shù) 2016年5期
    關(guān)鍵詞:激波側(cè)向壁面

    劉華偉,葉正寅,葉 坤

    (西北工業(yè)大學(xué) 翼型葉柵空氣動(dòng)力學(xué)國(guó)防科技重點(diǎn)實(shí)驗(yàn)室,西安 710072)

    ?

    火箭發(fā)動(dòng)機(jī)噴管擺動(dòng)對(duì)側(cè)向載荷的影響

    劉華偉,葉正寅,葉坤

    (西北工業(yè)大學(xué) 翼型葉柵空氣動(dòng)力學(xué)國(guó)防科技重點(diǎn)實(shí)驗(yàn)室,西安710072)

    針對(duì)實(shí)際環(huán)境下噴管存在不可避免的擺動(dòng)情況,以VOLVO S1噴管為研究對(duì)象,基于動(dòng)網(wǎng)格技術(shù),采用CFD數(shù)值方法模擬了噴管擺動(dòng)的非定常過(guò)程,在NPR=16的受限激波分離(RSS)情況下,進(jìn)行了多方位計(jì)算研究,分析了噴管不同擺動(dòng)狀況下的內(nèi)部流動(dòng)分離與側(cè)向載荷。結(jié)果表明,噴管擺動(dòng)會(huì)顯著影響側(cè)向載荷,主要表現(xiàn)在側(cè)向載荷隨著擺動(dòng)頻率先增大、后減小,當(dāng)噴管擺動(dòng)頻率接近噴管流動(dòng)分離流場(chǎng)的特征頻率時(shí),側(cè)向載荷達(dá)到最大值;同時(shí),擺動(dòng)條件下的側(cè)向載荷隨擺動(dòng)振幅增加基本呈線性增大。此外,對(duì)不同擺動(dòng)工況下側(cè)向載荷的分析表明,側(cè)向載荷主要來(lái)自收縮段。

    噴管;非定常流場(chǎng);側(cè)向載荷;擺動(dòng);頻率;振幅

    0 引言

    隨著航天推進(jìn)技術(shù)的快速發(fā)展,現(xiàn)代火箭發(fā)動(dòng)機(jī)愈來(lái)愈多采用大面積噴管獲得高推進(jìn)效率和大比沖來(lái)提高整個(gè)火箭工作性能。對(duì)于火箭發(fā)動(dòng)機(jī)而言,噴管內(nèi)部流動(dòng)狀況直接影響著整個(gè)火箭的性能。然而,火箭發(fā)動(dòng)機(jī)噴管的工作環(huán)境會(huì)跨越很寬的高度范圍,為了在整個(gè)飛行過(guò)程中獲得可接受的性能,噴管都是按中間環(huán)境壓力設(shè)計(jì)的,而大部分情況下,大面積比噴管工作在非設(shè)計(jì)條件下,當(dāng)背壓高于設(shè)計(jì)值一定程度時(shí),容易發(fā)生流動(dòng)分離。由于分離流的不穩(wěn)定性和非定常性,噴管內(nèi)部流動(dòng)分離往往呈現(xiàn)復(fù)雜的非軸對(duì)稱(chēng)性,引起嚴(yán)重的側(cè)向載荷。嚴(yán)重的側(cè)向載荷可能會(huì)導(dǎo)致噴管擺動(dòng)、發(fā)動(dòng)機(jī)壽命縮短、噴管結(jié)構(gòu)破壞等后果。例如,美國(guó)的J-2S發(fā)動(dòng)機(jī)、航天飛機(jī)主發(fā)動(dòng)機(jī)SSME、Fastrac發(fā)動(dòng)機(jī)、日本的LE-7A和歐洲的火神發(fā)動(dòng)機(jī)都曾遇到側(cè)向載荷的問(wèn)題[1-2]。

    目前,對(duì)噴管的研究一方面集中在對(duì)噴管流動(dòng)分離的模擬和側(cè)向載荷的仿真,文獻(xiàn)[3-5]對(duì)噴管不同分離模式下,進(jìn)行了數(shù)值仿真以及與實(shí)驗(yàn)的對(duì)比,文獻(xiàn)[6-7]對(duì)噴管側(cè)向載荷進(jìn)行數(shù)值分析;另一方面,主要著眼于噴管動(dòng)載荷方面的研究,文獻(xiàn)[8]分析了分離模式下噴管的氣動(dòng)響應(yīng),文獻(xiàn)[9-11]利用數(shù)值模擬方法,對(duì)噴管進(jìn)行流固耦合的研究,文獻(xiàn)[12]對(duì)火箭發(fā)動(dòng)機(jī)噴管進(jìn)行了有關(guān)側(cè)向載荷和氣動(dòng)彈性穩(wěn)定性的研究。然而,實(shí)際中化學(xué)反應(yīng)、流動(dòng)分離以及噴管柔性的因素,勢(shì)必會(huì)導(dǎo)致噴管振動(dòng),振動(dòng)對(duì)噴管的內(nèi)部流動(dòng)和側(cè)向載荷均有一定的影響。美國(guó)火箭預(yù)燃試驗(yàn)中,測(cè)得發(fā)動(dòng)機(jī)擺動(dòng)引起的載荷高達(dá)190個(gè)過(guò)載[13],在國(guó)內(nèi)航天相關(guān)部門(mén)的地面測(cè)試過(guò)程中,也出現(xiàn)過(guò)擺動(dòng)側(cè)向載荷導(dǎo)致試驗(yàn)結(jié)構(gòu)破壞的事件[14]。目前,關(guān)于噴管的擺動(dòng)對(duì)側(cè)向載荷的影響還沒(méi)有具體研究。因此,詳細(xì)掌握噴管擺動(dòng)對(duì)側(cè)向載荷的影響,對(duì)發(fā)動(dòng)機(jī)設(shè)計(jì)和實(shí)際工作性能分析具有重要的研究?jī)r(jià)值。

    本文將以VOLVO S1噴管為研究對(duì)象,模擬壓比NPR(Nozzle Pressure Ratio,即燃燒室總壓和環(huán)境壓強(qiáng)的比值)為16下的受限激波分離,分析無(wú)擺動(dòng)下受限激波分離模式下的側(cè)向載荷與分布,對(duì)不同擺動(dòng)狀態(tài)進(jìn)行數(shù)值計(jì)算,研究擺動(dòng)頻率和振幅對(duì)側(cè)向載荷的影響規(guī)律,分析不同部位側(cè)向載荷的分布情況。

    1 控制方程與數(shù)值方法

    任意拉格朗日歐拉(ALE,Arbitrary Lagrangian Eulerian Method)是將拉格朗日方法和歐拉方法統(tǒng)一起來(lái)。當(dāng)網(wǎng)格速度為零時(shí)方程變成歐拉方法;當(dāng)網(wǎng)格速度等于流體運(yùn)動(dòng)速度時(shí),則變成拉格朗日方法。ALE描述下的守恒型非定常Navier-Stokes方程積分形式為[15]

    其中

    空間離散采用AUSM+格式,時(shí)間推進(jìn)采用LU-SGS格式;同時(shí),本文引入了kw-SST湍流模型,噴管內(nèi)壁面和外壁面采用無(wú)滑移絕熱壁條件。

    2 研究對(duì)象

    2.1幾何建模

    本文為沃爾沃宇航公司的VOLVO S1噴管為研究對(duì)象[5]。VOLVO S1噴管是依據(jù)Vulcain噴管幾何外形設(shè)計(jì)的一個(gè)縮比模型,是用來(lái)研究Vulcain系列噴管的流動(dòng)分離和側(cè)向載荷的TOP(Thrust Optimised Parabolic)噴管。

    根據(jù)VOLVO S1基本參數(shù),可建立TOP噴管軸對(duì)稱(chēng)截面幾何模型。Rao提出TOP噴管擴(kuò)張段NE采用傾斜拋物線,并滿足Rao-Shmyglevsky 方程(2)。由L、Dt、De、rtd/rt、θN和θE6個(gè)參數(shù)確定變量b、c、d和e,從而確定擴(kuò)張段的幾何外形。

    (2)

    噴管基本構(gòu)型簡(jiǎn)圖見(jiàn)圖1,VOLVO S1噴管變量參數(shù)見(jiàn)表1。

    圖1 噴管基本構(gòu)型簡(jiǎn)圖

    參數(shù)變量參數(shù)值膨脹比∑20噴管長(zhǎng)度L/mm350喉道直徑Dt/mm67.08噴管出口直徑De/mm300無(wú)量綱的擴(kuò)張段入口半徑rtd/rt0.5入口角θN/(°)35.025擴(kuò)張角θE/(°)4.0燃燒室總壓p0/MPa5.0燃?xì)饪倻豑0/K450氣體空氣

    2.2網(wǎng)格劃分

    軸向以噴管長(zhǎng)度L為參考長(zhǎng)度,徑向以噴管出口半徑re為參考長(zhǎng)度,圓柱外場(chǎng)計(jì)算域軸向與徑向長(zhǎng)度分別取對(duì)應(yīng)參考長(zhǎng)度的20倍;噴管內(nèi)部壁面網(wǎng)格局部分塊加密,附面層網(wǎng)格的第一層高度為1×10-4mm,壁面y+接近于1,從而保證準(zhǔn)確捕捉壁面的流動(dòng)細(xì)節(jié)。為了準(zhǔn)確捕捉噴管內(nèi)非對(duì)稱(chēng)側(cè)向載荷,減小網(wǎng)格因素對(duì)流動(dòng)非對(duì)稱(chēng)的影響,計(jì)算域由半噴管二維網(wǎng)格以X軸(噴管的軸線方向)旋轉(zhuǎn)得到三維網(wǎng)格,整個(gè)計(jì)算域采用六面體結(jié)構(gòu)網(wǎng)格(圖2),未做非對(duì)稱(chēng)假設(shè)。

    噴管內(nèi)部流動(dòng)分離較復(fù)雜,對(duì)網(wǎng)格依賴性較高,為了校驗(yàn)網(wǎng)格收斂性,本文采用了三套粗細(xì)不同網(wǎng)格,如表2所示。

    圖2 噴管計(jì)算域網(wǎng)格

    網(wǎng)格類(lèi)別噴管內(nèi)部軸向節(jié)點(diǎn)徑向節(jié)點(diǎn)附面層計(jì)算域總網(wǎng)格數(shù)Ⅰ(粗)14712025308484Ⅱ(中)22112625581952Ⅲ(密)25813025900054

    3 數(shù)值計(jì)算

    3.1實(shí)驗(yàn)結(jié)果對(duì)比與網(wǎng)格無(wú)關(guān)性分析

    VOLVO曾對(duì)S1噴管進(jìn)行冷流實(shí)驗(yàn),?stlund[5]提供了NPR=16條件下的受限激波分離的實(shí)驗(yàn)結(jié)果。為了校驗(yàn)本文所用數(shù)值方法計(jì)算結(jié)果的準(zhǔn)確性,同時(shí)考慮到網(wǎng)格密度對(duì)流動(dòng)分離的影響,對(duì)網(wǎng)格Ⅰ、網(wǎng)格Ⅱ、網(wǎng)格Ⅲ進(jìn)行了計(jì)算校驗(yàn),計(jì)算結(jié)果見(jiàn)圖3。

    網(wǎng)格Ⅱ計(jì)算得到的壁面壓力分布,與經(jīng)過(guò)加密的網(wǎng)格Ⅲ的結(jié)果相比,變化不大,說(shuō)明此時(shí)計(jì)算得到的穩(wěn)態(tài)收斂解與網(wǎng)格無(wú)關(guān),即NPR=16下,計(jì)算所采用的網(wǎng)格Ⅱ在計(jì)算收斂時(shí),基本可準(zhǔn)確模擬真實(shí)模擬噴管的內(nèi)部流動(dòng)。

    圖3 NPR=16不同網(wǎng)格類(lèi)型下的壁面壓力分布線

    3.2NPR=16無(wú)擺動(dòng)條件下的側(cè)向載荷

    燃燒室總壓為1.6 MPa,總溫為450 K,環(huán)境壓強(qiáng)為0.1 MPa,環(huán)境溫度為288 K,無(wú)擺動(dòng)的噴管進(jìn)行流場(chǎng)數(shù)值模擬NPR=16工況,計(jì)算得到受限激波分離下的噴管壁面壓力及馬赫云圖,見(jiàn)圖4。

    圖4中的壁面壓力分布圖顯示了受限激波分離模式下的特征點(diǎn)的定義。壁面壓力分布分析:

    (1)在X0~Xi,沿軸向噴管壁面壓強(qiáng)逐漸減小,流動(dòng)較穩(wěn)定。

    (2)在Xi~Xp激波附面層相互作用,流動(dòng)在分離區(qū)初始作用點(diǎn)Xi壁面壓力pi開(kāi)始偏離真空條件下的噴管壁面壓力值pw,在分離點(diǎn)Xs處流動(dòng)發(fā)生分離,在Xp達(dá)到一個(gè)相對(duì)短暫穩(wěn)定的壓強(qiáng);逆壓梯度使得Xp~Xr間壁面出現(xiàn)流動(dòng)分離區(qū),沿著軸向壓力值增大,在Xr處的壓強(qiáng)pr達(dá)到最大值,高于環(huán)境壓力pa。

    (3)Xr~Xe流動(dòng)再附著區(qū):由于pr高于pa,Xr之后的流動(dòng)發(fā)生再次附著壁面,壁面壓力降至pe,誘導(dǎo)壁面出現(xiàn)膨脹波,噴管內(nèi)出現(xiàn)強(qiáng)烈的非定常壓力擺動(dòng)。

    圖4 NPR=16下受限激波分離(RSS)壁面壓力分布與馬赫云圖的對(duì)比

    圖4中的馬赫云圖顯示了受限激波分離模式下的噴管內(nèi)部馬赫云圖。從激波相互作用分析:噴管內(nèi)內(nèi)激波與正激波相互作用產(chǎn)生反馬赫反射,偏離中軸線與壁面分離激波發(fā)生干涉形成錐狀激波,噴管內(nèi)形成特殊“帽”狀激波形態(tài),正激波與錐狀激波間形成亞音速的駐渦,足夠強(qiáng)的錐狀激波,使得超聲速射流經(jīng)過(guò)較短的分離區(qū),再次附著于壁面,故名受限激波分離(Restricted Shock Separation,RSS)。

    為研究噴管各部分側(cè)向載荷的分布,本文將噴管分為兩部分—收縮段(Convergent)和擴(kuò)張段(Divergent)。其中,擴(kuò)張段根據(jù)馬赫數(shù)云圖分布進(jìn)一步細(xì)化分為流動(dòng)穩(wěn)定區(qū)A(擴(kuò)張段_A)、分離區(qū)B(擴(kuò)張段_B)和再附著區(qū)C(擴(kuò)張段_C)(如圖4所示)。

    無(wú)擺動(dòng)噴管在NPR=16的工況下,計(jì)算達(dá)到穩(wěn)定收斂時(shí),噴管軸向力為12 081 N,Y方向的最大側(cè)向載荷僅為0.061 6 N。與軸向力相比而言,Y方向的側(cè)向載荷在數(shù)值上是一個(gè)微小量,且其基本由擴(kuò)張段貢獻(xiàn)。結(jié)果表明,NPR=16在無(wú)擺動(dòng)工況下,噴管的側(cè)向載荷極小(圖5)。

    選取噴管Y方向側(cè)向載荷隨時(shí)間的變化歷程為研究對(duì)象,對(duì)其進(jìn)行快速傅里葉變換(Fast Fourier Transform, FFT),得到噴管側(cè)向載荷FFT結(jié)果如圖6所示。結(jié)果表明,峰值脈動(dòng)頻率為2 395.21 Hz, 其次還有2個(gè)占比重較高的分支,其頻率分別為1 596.81 Hz 和2 994.01 Hz。計(jì)算時(shí)間步長(zhǎng)為1×10-5s,就高頻2 994.01 Hz而言,也足以保證在1個(gè)周期內(nèi)不低于30個(gè)時(shí)間步,從而較準(zhǔn)確地捕捉各個(gè)頻率脈動(dòng)。同時(shí),噴管側(cè)向載荷脈動(dòng)頻率范圍完全可能與周?chē)Y(jié)構(gòu)頻率接近,作為激勵(lì)載荷, 可能引起結(jié)構(gòu)振動(dòng)問(wèn)題。

    圖5 NPR=16無(wú)擺動(dòng)噴管的側(cè)向載荷

    圖6 噴管側(cè)向載荷FFT結(jié)果

    3.3NPR=16擺動(dòng)噴管的側(cè)向載荷

    為模擬擺動(dòng)噴管下的內(nèi)部流動(dòng)情況,分析側(cè)向載荷分布情況,在設(shè)置邊界條件時(shí),定義噴管內(nèi)壁面為WALL,定義噴管外壁面為WALL_FAR,借助動(dòng)網(wǎng)格技術(shù), WALL和WALL_FAR分別在Y方向上采用正弦擺動(dòng)函數(shù),針對(duì)采用不同的擺動(dòng)頻率和振幅,分別進(jìn)行比較分析。

    3.3.1側(cè)向載荷隨擺動(dòng)頻率的變化

    在噴管擺動(dòng)振幅為喉道直徑的1%(即0.670 8 mm)工況下,分別對(duì)頻率為10、50、100、200、400、500、800、1 000、1 500、2 000、2 500、3 000、3 500、4 000、4 500、5 000等不同頻率的狀態(tài)進(jìn)行了數(shù)值模擬仿真。為了能準(zhǔn)確模擬側(cè)向載荷隨時(shí)間的變化,將根據(jù)不同的頻率設(shè)置時(shí)間步長(zhǎng),確保不同頻率下的一個(gè)周期中,有50個(gè)時(shí)間步來(lái)捕捉側(cè)向載荷,從而計(jì)算得到側(cè)向載荷隨時(shí)間的變化結(jié)果。

    結(jié)合圖7和圖8,隨著頻率增加側(cè)向載荷先增大、后減小,在2 500 Hz下,Y向側(cè)向載荷達(dá)到峰值9 472.8 N。同時(shí)可見(jiàn),收縮段提供側(cè)向載荷起著主導(dǎo)作用,隨著頻率變化,收縮段側(cè)向載荷變化非常顯著。

    其中,f=2 500 Hz工況下的結(jié)果如圖9和圖10所示。

    圖7 噴管Y方向側(cè)向力隨頻率的變化曲線

    圖8 噴管不同段的側(cè)向載荷分布

    圖9 f=2 500 Hz側(cè)向載荷隨時(shí)間步長(zhǎng)關(guān)系

    圖10 f=2 500 Hz噴管側(cè)向載荷FFT結(jié)果

    結(jié)果分析表明,在正弦擺動(dòng)運(yùn)動(dòng)下,側(cè)向載荷隨時(shí)間也呈正弦函數(shù)變化,沿軸向不同段的側(cè)向載荷有一定的滯后效應(yīng)。對(duì)f=2 500 Hz側(cè)向力進(jìn)行FFT分析,發(fā)現(xiàn)Y側(cè)向載荷主頻為2 495.01 Hz。還有對(duì)應(yīng)頻率為7 485.0 Hz較小權(quán)重分支,高階頻率側(cè)向力幅值很低。與圖6對(duì)比,說(shuō)明f=2 395.21 Hz的頻率分支是產(chǎn)生大的側(cè)向力的主要因素。值得注意的是f=2 500 Hz下的擺動(dòng)噴管流場(chǎng)結(jié)構(gòu)與無(wú)擺動(dòng)噴管流場(chǎng)結(jié)構(gòu)相似,由于此工況下噴管擺動(dòng)頻率與無(wú)擺動(dòng)下噴管流動(dòng)分離流場(chǎng)的特征頻率十分接近,引起共振,從而使得側(cè)向載荷達(dá)到最大值。

    3.3.2側(cè)向載荷隨擺動(dòng)振幅的變化

    噴管擺動(dòng)頻率為500 Hz工況下,取擺動(dòng)振幅為1%Dt、2%Dt、3%Dt、4%Dt、5 %Dt不同振幅下工況進(jìn)行數(shù)值模擬仿真,計(jì)算結(jié)果如圖11所示。

    圖11 f=500 Hz噴管的側(cè)向載荷隨振幅的關(guān)系

    噴管在擺動(dòng)頻率為500 Hz工況下,計(jì)算結(jié)果表明,隨著振幅增大,整個(gè)噴管Y方向的側(cè)向載荷幾乎線性增大(圖12)。

    通過(guò)監(jiān)測(cè)不同區(qū)段的側(cè)向載荷可看出,隨著振幅增大,收縮段側(cè)向載荷在整個(gè)側(cè)向載荷中貢獻(xiàn)較大,均達(dá)到整個(gè)噴管的側(cè)向載荷50%左右;其次,再附著區(qū)C(擴(kuò)張段_C)與分離區(qū)B(擴(kuò)張段_B)側(cè)向載荷相當(dāng),流動(dòng)穩(wěn)定區(qū)A(擴(kuò)張段_A)由于流動(dòng)受擺動(dòng)影響相對(duì)較小,側(cè)向載荷變化最小。

    圖12 f=500 Hz噴管各段側(cè)向載荷隨振幅的變化

    從擺動(dòng)噴管不同振幅下流線圖(圖13)分析可知,無(wú)擺動(dòng)下噴管內(nèi)部流動(dòng)較穩(wěn)定,超聲速流上下基本呈對(duì)稱(chēng)狀態(tài)。

    (a) 無(wú)擺動(dòng)噴管   (b) 振幅為1%Dt 的擺動(dòng)噴管

    (c)振幅為2%Dt 的擺動(dòng)噴管 (d) 振幅為3%Dt的擺動(dòng)噴管

    (e) 振幅為4%Dt的擺動(dòng)噴管 (f) 振幅為5%Dt的擺動(dòng)噴管

    隨著振幅的增大,噴管內(nèi)部流動(dòng)由于擺動(dòng)的影響,不對(duì)稱(chēng)性變得愈來(lái)愈明顯,高壓區(qū)收縮段小擺動(dòng)引起的不對(duì)稱(chēng)性會(huì)激發(fā)較大的側(cè)向載荷。擴(kuò)張段內(nèi)帽狀激波后的出現(xiàn)不同程度大小的渦,上下分離點(diǎn)的位置發(fā)生偏離,促使分離區(qū)B(擴(kuò)張段_B)側(cè)向載荷的變化,再附著區(qū)C(擴(kuò)張段_C)段側(cè)向載荷的產(chǎn)生主要是由于超聲速流在噴管內(nèi)外均出現(xiàn)擺動(dòng)現(xiàn)象,使得壁面壓力發(fā)生劇烈變化,形成較大的壓力差。此外,噴管外的流動(dòng)隨著振幅增大,影響越明顯,影響范圍越大。

    4 結(jié)論

    (1)當(dāng)噴管在無(wú)擺動(dòng)條件下,側(cè)向載荷很小,基本由擴(kuò)張段提供。當(dāng)噴管在不同擺動(dòng)運(yùn)動(dòng)情況下,側(cè)向載荷會(huì)顯著增大,側(cè)向載荷則主要來(lái)自收縮段。

    (2)噴管做正弦擺動(dòng),且擺動(dòng)振幅相同時(shí),隨著擺動(dòng)頻率的增加,側(cè)向載荷先增大、后減小。當(dāng)擺動(dòng)頻率增大到f=2 500 Hz,接近噴管流動(dòng)分離流場(chǎng)的特征頻率時(shí),側(cè)向載荷取得最大值。此時(shí),收縮段提供側(cè)向載荷占總側(cè)向載荷的90%。

    (3)當(dāng)擺動(dòng)頻率恒定時(shí),隨擺動(dòng)振幅的增大,側(cè)向載荷基本呈線性增加。其中,收縮段側(cè)向載荷的貢獻(xiàn)突出,其占整個(gè)噴管側(cè)向載荷的50%左右,再附著區(qū)C(擴(kuò)張段_C)由于壁面附近超聲速流的非定常效應(yīng)和分離區(qū)B(擴(kuò)張段_B) 分離點(diǎn)位置的偏離產(chǎn)生的側(cè)向載荷均相當(dāng)明顯。

    [1]Hadjadj A, Onofri M. Nozzle flow separation[J]. Shock Waves, 2009, 19(3): 163-169.

    [2]Hagemann G, Frey M. Shock pattern in the plume of rocket nozzles: needs for design consideration[J]. Shock Waves, 2008, 17(8): 387-395.

    [3]Martelli E, Nasuti F, Onofri M. Numerical calculation of FSS/RSS transition in highly overexpanded rocket nozzle flows[J]. Shock Waves, 2010, 20(2): 139-46.

    [4]Lijo V, Kim H D, Setoguchi T, et al. Numerical simulation of transient flows in a rocket propulsion nozzle[J]. International Journal of Heat and Fluid Flow, 2010, 31: 409-417.

    [5]?stlund J. Flow processes in rocket engine nozzles with focus on flow separation and side-loads[R]. TRITA-MEK Technical Report, 2002: 9.

    [6]Smalley K B, Brown A M, Ruf J, et al. Flow separation side loads excitation of rocket nozzle FEM[C]//48th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Con., Honolulu, Hawaii, USA, 2007.

    [7]劉亞冰, 王長(zhǎng)輝, 許曉勇.噴管分離流動(dòng)與側(cè)向載荷定常數(shù)值模擬[J]. 航空動(dòng)力學(xué)報(bào), 2008, 23(11): 2115-2118.

    [8]Brown A M, Ruf J, Reed D, et al. Characterization of side load phenomena using measurement of fluid/structure interaction[C]//38th Joint Propulsion Conference and Exhibit, AIAA Paper 2002-3999.

    [9]Lefrancois E. Numerical validation of a stability model for a flexible over-expanded rocket nozzle[J]. International Journal for Numerical Methods in Fuids, 2005, 49(4): 349-369.

    [10]Zhang S J, Fuchiwaki T, Zhao X. Aeroelastic coupling and side loads in rocket nozzles[R]. AIAA 2008-4064.

    [11]Garelli L, Paz R R, Storti M A. Fluid-structure interaction study of the start-up of a rocket engine nozzle[J]. Computers & Fluids, 2010, 39(7): 1208-1218.

    [12]Lu Guang-liang, Ye Zheng-yin. Investigation on the mechanism of aeroelastic hazard during ground test of rochet nozzle[J]. Science China Technological Sciences, 2012, 55(9): 2462-2473.

    [13]Flandro G A, Majdalani J, Sims J D. Nonlinear longitudinal mode instability in liquid propellant rocket engine preburners[R]. AIAA 2004-4162.

    [14]葉正寅, 呂廣亮. 火箭發(fā)動(dòng)機(jī)非定常側(cè)向力和流固耦合研究發(fā)展[J]. 航空工程進(jìn)展, 2015, 6(1): 1-12.

    [15]許和勇. 基于動(dòng)態(tài)嵌套網(wǎng)格的非定常流動(dòng)數(shù)值模擬[D].西安:西北工業(yè)大學(xué), 2009.

    (編輯:崔賢彬)

    Research on the influence of swing on side-load in rocket motor

    LIU Hua-wei, YE Zheng-yin, YE Kun

    (National Key Laboratory of Aerodynamic Design and Research,Northwestern Polytechnical University, Xi'an710072, China)

    Side-loads in nozzle can cause serious structural damage. Owing to the presence of the inevitable swing nozzle in the actual environmental situation, the unsteady flow field of the swinging nozzle, which is based on dynamic mesh technology, was simulated by taking VOLVO S1 nozzle as the object of study and adopting CFD numerical simulation method. Meanwhile, a multi-faceted computing research was conducted on the restricted shock separation (RSS) atNPR=16. The internal nozzle flow separation and side-loads were analyzed under different swing states .The result shows that nozzle swing has notable effect on side-loads in nozzle, which is mainly embodied as the followings: Side-load enhances first,then decreases later with swing frequency. The loads attain the maximum while nozzle swing frequency approaches the frequency of flow separation. Side-load increases linearly with swing amplitude. Besides, the research on the analysis of side-loads in different swinging working conditions shows that side-load mainly comes from the convergent part.

    nozzle;unsteady flow field;side-loads;swing;frequency;amplitude

    2015-07-17;

    2015-09-21。

    國(guó)家自然科學(xué)基金(11272262)。

    劉華偉(1990—),男,碩士生,研究方向?yàn)榛鸺l(fā)動(dòng)機(jī)側(cè)向載荷研究。E-mail:huaweiliu@nwpu.edu.cn

    V435

    A

    1006-2793(2016)05-0619-06

    10.7673/j.issn.1006-2793.2016.05.004

    猜你喜歡
    激波側(cè)向壁面
    二維有限長(zhǎng)度柔性壁面上T-S波演化的數(shù)值研究
    軍航無(wú)人機(jī)與民航航班側(cè)向碰撞風(fēng)險(xiǎn)評(píng)估
    一種基于聚類(lèi)分析的二維激波模式識(shí)別算法
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    斜激波入射V形鈍前緣溢流口激波干擾研究
    適于可壓縮多尺度流動(dòng)的緊致型激波捕捉格式
    彎月薄鏡的側(cè)向支撐
    壁面溫度對(duì)微型內(nèi)燃機(jī)燃燒特性的影響
    側(cè)向風(fēng)場(chǎng)中無(wú)人機(jī)的飛行研究
    基于 L1自適應(yīng)控制的無(wú)人機(jī)橫側(cè)向控制
    一本色道久久久久久精品综合| 免费观看无遮挡的男女| 久久这里有精品视频免费| 黑人高潮一二区| 免费黄网站久久成人精品| 久久99蜜桃精品久久| 最近的中文字幕免费完整| 在线看a的网站| 国产亚洲av嫩草精品影院| 国产黄片美女视频| 国产av国产精品国产| a级一级毛片免费在线观看| 黑人高潮一二区| av在线老鸭窝| 亚洲av.av天堂| av黄色大香蕉| 亚洲av男天堂| 极品少妇高潮喷水抽搐| www.色视频.com| 亚洲精品中文字幕在线视频 | 少妇丰满av| 精品酒店卫生间| 18禁动态无遮挡网站| 亚洲高清免费不卡视频| h日本视频在线播放| videossex国产| 夜夜看夜夜爽夜夜摸| 老司机影院成人| 天堂中文最新版在线下载 | 在现免费观看毛片| 两个人的视频大全免费| 欧美xxxx性猛交bbbb| 九色成人免费人妻av| 亚洲伊人久久精品综合| 永久网站在线| 亚洲欧美成人精品一区二区| 国产精品久久久久久精品古装| 久久韩国三级中文字幕| 精品国产露脸久久av麻豆| 亚洲av免费高清在线观看| 欧美精品人与动牲交sv欧美| 国产乱来视频区| 精品人妻视频免费看| 18禁动态无遮挡网站| 99热网站在线观看| 国产成人a区在线观看| 高清av免费在线| 午夜视频国产福利| 欧美日韩视频高清一区二区三区二| 熟女av电影| 欧美zozozo另类| 精品国产露脸久久av麻豆| 好男人在线观看高清免费视频| 国产爱豆传媒在线观看| 少妇 在线观看| 久久国产乱子免费精品| 99热这里只有精品一区| 久久99热这里只有精品18| 国产精品麻豆人妻色哟哟久久| 一级片'在线观看视频| 看黄色毛片网站| 观看美女的网站| 夫妻性生交免费视频一级片| 国产亚洲精品久久久com| av网站免费在线观看视频| 建设人人有责人人尽责人人享有的 | 亚洲激情五月婷婷啪啪| 天堂中文最新版在线下载 | 国产精品久久久久久精品古装| 欧美bdsm另类| 亚洲av成人精品一区久久| 2018国产大陆天天弄谢| 亚洲精华国产精华液的使用体验| 日本av手机在线免费观看| 免费av不卡在线播放| 2021少妇久久久久久久久久久| 国产女主播在线喷水免费视频网站| 免费在线观看成人毛片| 黄片wwwwww| 亚洲人成网站高清观看| 亚洲欧洲日产国产| 久久鲁丝午夜福利片| 99久国产av精品国产电影| 免费高清在线观看视频在线观看| 99热这里只有精品一区| 一级片'在线观看视频| 一个人看视频在线观看www免费| 亚洲高清免费不卡视频| 欧美精品一区二区大全| 啦啦啦啦在线视频资源| 高清日韩中文字幕在线| 亚洲精品一二三| 午夜福利视频精品| 日韩人妻高清精品专区| 亚洲精品456在线播放app| 各种免费的搞黄视频| 人人妻人人爽人人添夜夜欢视频 | av在线播放精品| 身体一侧抽搐| 男人狂女人下面高潮的视频| 七月丁香在线播放| 亚洲无线观看免费| 最后的刺客免费高清国语| 在线观看av片永久免费下载| 男人狂女人下面高潮的视频| 日韩三级伦理在线观看| 亚洲性久久影院| 国产精品久久久久久精品电影小说 | 狂野欧美白嫩少妇大欣赏| 26uuu在线亚洲综合色| 国产精品av视频在线免费观看| 亚洲av中文av极速乱| 欧美区成人在线视频| 欧美+日韩+精品| 少妇的逼水好多| 在线 av 中文字幕| 精品国产三级普通话版| 亚洲欧美成人精品一区二区| 亚洲精品日本国产第一区| 久久久久精品久久久久真实原创| 看黄色毛片网站| 国产精品熟女久久久久浪| 看非洲黑人一级黄片| 小蜜桃在线观看免费完整版高清| 精品久久久久久久末码| 成人无遮挡网站| 日韩免费高清中文字幕av| 91在线精品国自产拍蜜月| 成人特级av手机在线观看| 一级毛片 在线播放| 日韩中字成人| 又大又黄又爽视频免费| 亚洲性久久影院| 日韩一区二区三区影片| 国产精品三级大全| 我的老师免费观看完整版| 日日啪夜夜爽| 六月丁香七月| 只有这里有精品99| 成年女人在线观看亚洲视频 | av在线天堂中文字幕| 男的添女的下面高潮视频| 真实男女啪啪啪动态图| 久久韩国三级中文字幕| 97人妻精品一区二区三区麻豆| 人人妻人人看人人澡| 日本欧美国产在线视频| 精品一区二区免费观看| 尾随美女入室| 久久久久精品性色| 亚州av有码| 波野结衣二区三区在线| 特大巨黑吊av在线直播| 我要看日韩黄色一级片| 国产高潮美女av| 国产精品麻豆人妻色哟哟久久| 大片免费播放器 马上看| 亚洲精品成人久久久久久| 97超碰精品成人国产| av女优亚洲男人天堂| 亚洲av.av天堂| 久久久久久久久久人人人人人人| 亚洲av不卡在线观看| 一级毛片久久久久久久久女| 久久久亚洲精品成人影院| 国产色婷婷99| 美女cb高潮喷水在线观看| 久久久成人免费电影| 波野结衣二区三区在线| 亚洲久久久久久中文字幕| 校园人妻丝袜中文字幕| 免费看a级黄色片| 久久久欧美国产精品| 人妻少妇偷人精品九色| 精品午夜福利在线看| 免费观看在线日韩| 日日啪夜夜爽| 嫩草影院入口| 亚洲国产av新网站| 老司机影院毛片| 国产久久久一区二区三区| 亚洲色图综合在线观看| 婷婷色综合大香蕉| 不卡视频在线观看欧美| 国产精品久久久久久久久免| 天美传媒精品一区二区| 91aial.com中文字幕在线观看| 人妻一区二区av| 韩国av在线不卡| 久久影院123| av播播在线观看一区| 国产一区二区三区综合在线观看 | 久久久久久久亚洲中文字幕| 久久国内精品自在自线图片| 亚洲自偷自拍三级| 永久免费av网站大全| 亚洲欧洲国产日韩| 一级毛片电影观看| av免费观看日本| 黄色一级大片看看| 婷婷色av中文字幕| 欧美xxxx黑人xx丫x性爽| 亚洲精品乱久久久久久| 亚洲精品亚洲一区二区| 免费看光身美女| videos熟女内射| 乱系列少妇在线播放| 美女xxoo啪啪120秒动态图| 久久久久精品久久久久真实原创| 你懂的网址亚洲精品在线观看| 免费人成在线观看视频色| 国产精品一二三区在线看| av国产精品久久久久影院| 日日撸夜夜添| 日韩av在线免费看完整版不卡| 亚洲欧美日韩东京热| 男女那种视频在线观看| 国产黄色视频一区二区在线观看| 晚上一个人看的免费电影| 一级av片app| 日韩亚洲欧美综合| 欧美xxxx黑人xx丫x性爽| 亚洲高清免费不卡视频| 一区二区三区四区激情视频| 久久精品国产自在天天线| 国产伦在线观看视频一区| 热re99久久精品国产66热6| 黄片无遮挡物在线观看| 亚洲最大成人中文| 极品教师在线视频| 久久6这里有精品| 精品国产三级普通话版| 岛国毛片在线播放| 高清毛片免费看| 中文字幕亚洲精品专区| 免费观看的影片在线观看| 久久午夜福利片| 亚洲精品影视一区二区三区av| 卡戴珊不雅视频在线播放| 国产亚洲91精品色在线| 中文字幕av成人在线电影| 人妻一区二区av| 亚洲精品中文字幕在线视频 | .国产精品久久| 中文在线观看免费www的网站| 亚洲最大成人手机在线| 视频区图区小说| 人人妻人人澡人人爽人人夜夜| 中文字幕人妻熟人妻熟丝袜美| 午夜福利视频1000在线观看| 国产一区亚洲一区在线观看| 精品人妻偷拍中文字幕| 亚洲无线观看免费| 99热这里只有精品一区| 国产精品女同一区二区软件| 最新中文字幕久久久久| 人妻一区二区av| 国产v大片淫在线免费观看| 国精品久久久久久国模美| 亚洲av.av天堂| 国内少妇人妻偷人精品xxx网站| 午夜视频国产福利| 国产视频内射| 在线观看免费高清a一片| 七月丁香在线播放| 大话2 男鬼变身卡| 80岁老熟妇乱子伦牲交| 国产精品一二三区在线看| 国产免费视频播放在线视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费看av在线观看网站| 六月丁香七月| 大片电影免费在线观看免费| 亚洲精品亚洲一区二区| 精品少妇黑人巨大在线播放| 亚洲国产成人一精品久久久| 在线 av 中文字幕| 18禁在线播放成人免费| 99久久精品热视频| 三级国产精品欧美在线观看| 国产毛片在线视频| 别揉我奶头 嗯啊视频| 在现免费观看毛片| 青春草亚洲视频在线观看| av免费观看日本| 日本免费在线观看一区| 久久亚洲国产成人精品v| 亚洲久久久久久中文字幕| 国产亚洲精品久久久com| 国产亚洲最大av| 少妇被粗大猛烈的视频| 熟女电影av网| 91精品一卡2卡3卡4卡| 在线亚洲精品国产二区图片欧美 | av在线亚洲专区| 久久精品国产a三级三级三级| 久久韩国三级中文字幕| 亚洲av福利一区| 久久热精品热| 国产成人一区二区在线| 亚洲国产成人一精品久久久| 日本欧美国产在线视频| 99九九线精品视频在线观看视频| 久久久成人免费电影| 日韩一本色道免费dvd| 一个人看视频在线观看www免费| 国产精品麻豆人妻色哟哟久久| 亚洲人与动物交配视频| 天天躁日日操中文字幕| 2018国产大陆天天弄谢| 婷婷色综合大香蕉| 欧美精品一区二区大全| 网址你懂的国产日韩在线| 夜夜爽夜夜爽视频| 亚洲久久久久久中文字幕| 国产精品一区www在线观看| 日本猛色少妇xxxxx猛交久久| 少妇人妻一区二区三区视频| 成人毛片60女人毛片免费| 国产精品一区二区在线观看99| av专区在线播放| 别揉我奶头 嗯啊视频| a级毛色黄片| 毛片女人毛片| 国精品久久久久久国模美| 日韩av免费高清视频| 草草在线视频免费看| 午夜日本视频在线| 亚洲国产av新网站| 在线观看一区二区三区| 国产精品人妻久久久影院| 欧美性猛交╳xxx乱大交人| 好男人视频免费观看在线| 男人爽女人下面视频在线观看| 国产免费视频播放在线视频| 欧美性感艳星| 亚洲av免费在线观看| 精品国产三级普通话版| 日本与韩国留学比较| 亚洲精品色激情综合| 国产精品偷伦视频观看了| 日韩av不卡免费在线播放| 看免费成人av毛片| 欧美激情国产日韩精品一区| 人妻 亚洲 视频| 中文字幕久久专区| 免费观看av网站的网址| 色吧在线观看| 久久久久九九精品影院| 久久精品国产亚洲av天美| 国产高清不卡午夜福利| 亚洲av免费在线观看| 亚洲欧美日韩另类电影网站 | 亚洲精品一二三| 日本与韩国留学比较| 国产精品女同一区二区软件| 人妻系列 视频| 色网站视频免费| 我的女老师完整版在线观看| 亚洲欧美中文字幕日韩二区| 在线精品无人区一区二区三 | 精品久久久久久久人妻蜜臀av| 成人黄色视频免费在线看| 2021少妇久久久久久久久久久| 久热这里只有精品99| 亚洲精品,欧美精品| 久久精品久久精品一区二区三区| 国产乱人偷精品视频| 永久网站在线| 国产一区亚洲一区在线观看| 麻豆乱淫一区二区| 国产在线一区二区三区精| 男的添女的下面高潮视频| 精品少妇久久久久久888优播| 久久精品国产鲁丝片午夜精品| 中文天堂在线官网| 王馨瑶露胸无遮挡在线观看| 日韩电影二区| 少妇人妻一区二区三区视频| 好男人视频免费观看在线| 亚洲国产成人一精品久久久| 亚洲,欧美,日韩| 久久精品人妻少妇| 久久ye,这里只有精品| 国产v大片淫在线免费观看| 成人毛片60女人毛片免费| 综合色av麻豆| 青春草国产在线视频| 看十八女毛片水多多多| 久久久久精品久久久久真实原创| 最近手机中文字幕大全| 成人一区二区视频在线观看| 久久99精品国语久久久| 97热精品久久久久久| 国产欧美日韩精品一区二区| 91aial.com中文字幕在线观看| 2021少妇久久久久久久久久久| 日韩视频在线欧美| 制服丝袜香蕉在线| 久久久国产一区二区| 51国产日韩欧美| 三级经典国产精品| 精品一区二区免费观看| 在线免费十八禁| 午夜免费观看性视频| av国产免费在线观看| 性色avwww在线观看| 欧美潮喷喷水| 寂寞人妻少妇视频99o| 成人漫画全彩无遮挡| 极品教师在线视频| 交换朋友夫妻互换小说| 一级二级三级毛片免费看| 蜜桃亚洲精品一区二区三区| 美女脱内裤让男人舔精品视频| 亚洲精品国产成人久久av| 成年av动漫网址| 成年女人看的毛片在线观看| 色哟哟·www| 成年免费大片在线观看| 免费av观看视频| 亚洲自拍偷在线| 男男h啪啪无遮挡| 老师上课跳d突然被开到最大视频| 日产精品乱码卡一卡2卡三| 亚洲精品久久午夜乱码| 别揉我奶头 嗯啊视频| av天堂中文字幕网| 欧美性猛交╳xxx乱大交人| 熟妇人妻不卡中文字幕| 日韩一本色道免费dvd| 黄色怎么调成土黄色| 亚洲av成人精品一区久久| 国产成人福利小说| 成人午夜精彩视频在线观看| 日本av手机在线免费观看| 国产伦理片在线播放av一区| 少妇猛男粗大的猛烈进出视频 | 久久久久久久午夜电影| 午夜免费观看性视频| 深爱激情五月婷婷| 国产精品99久久久久久久久| 久久韩国三级中文字幕| 99久久九九国产精品国产免费| 免费观看无遮挡的男女| 国产欧美亚洲国产| 性色avwww在线观看| 在线a可以看的网站| 久久久久网色| 日韩欧美 国产精品| 麻豆乱淫一区二区| 自拍偷自拍亚洲精品老妇| 亚洲人成网站在线播| 欧美bdsm另类| 国产精品麻豆人妻色哟哟久久| 最近最新中文字幕免费大全7| 精品一区二区三卡| 97在线视频观看| 亚洲在久久综合| 国产综合精华液| 久久女婷五月综合色啪小说 | 大又大粗又爽又黄少妇毛片口| 最近中文字幕2019免费版| 看免费成人av毛片| 99热网站在线观看| 久久99热6这里只有精品| 国产成人freesex在线| 婷婷色综合大香蕉| 国产亚洲一区二区精品| 国产美女午夜福利| 新久久久久国产一级毛片| 欧美激情久久久久久爽电影| 777米奇影视久久| 91精品国产九色| 九九在线视频观看精品| 国产视频首页在线观看| 精品亚洲乱码少妇综合久久| 男人舔奶头视频| 中文乱码字字幕精品一区二区三区| 国产成人freesex在线| 91久久精品国产一区二区三区| 国产色婷婷99| 男女下面进入的视频免费午夜| 水蜜桃什么品种好| 欧美最新免费一区二区三区| 亚洲,一卡二卡三卡| 80岁老熟妇乱子伦牲交| 亚洲国产精品专区欧美| 国产成人aa在线观看| 三级国产精品欧美在线观看| 欧美激情在线99| 亚洲无线观看免费| 亚洲精品国产色婷婷电影| 久久精品国产a三级三级三级| 国产精品无大码| 天堂俺去俺来也www色官网| 亚洲av日韩在线播放| 国产免费一区二区三区四区乱码| 亚洲色图av天堂| 国产亚洲最大av| 日韩 亚洲 欧美在线| 特大巨黑吊av在线直播| 中文字幕av成人在线电影| 一区二区三区乱码不卡18| 看十八女毛片水多多多| 欧美性感艳星| 人人妻人人看人人澡| 搡老乐熟女国产| 一级爰片在线观看| 日本色播在线视频| 91在线精品国自产拍蜜月| 看黄色毛片网站| 久久久精品94久久精品| 联通29元200g的流量卡| 99久久精品热视频| 成人鲁丝片一二三区免费| 男女边摸边吃奶| 高清欧美精品videossex| 国产精品一区二区三区四区免费观看| 国产国拍精品亚洲av在线观看| 极品少妇高潮喷水抽搐| 日本wwww免费看| 久久久久久久久久久丰满| 免费电影在线观看免费观看| 极品教师在线视频| 日日啪夜夜撸| 蜜桃亚洲精品一区二区三区| 黄片无遮挡物在线观看| 天堂俺去俺来也www色官网| 免费看a级黄色片| 久久人人爽人人片av| 亚洲精品久久久久久婷婷小说| 人妻少妇偷人精品九色| 亚洲av欧美aⅴ国产| 免费黄色在线免费观看| 亚洲精品日韩在线中文字幕| 青春草国产在线视频| 男的添女的下面高潮视频| 久久99热这里只频精品6学生| 成人免费观看视频高清| 国产黄片视频在线免费观看| 80岁老熟妇乱子伦牲交| 欧美性感艳星| 亚洲精品一二三| 高清av免费在线| 成年女人在线观看亚洲视频 | 色视频在线一区二区三区| 国产欧美亚洲国产| 最近中文字幕高清免费大全6| 久久亚洲国产成人精品v| 在线观看国产h片| 久久久精品免费免费高清| av国产免费在线观看| 嫩草影院新地址| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 大香蕉久久网| 欧美少妇被猛烈插入视频| 日韩电影二区| 综合色av麻豆| 国产精品久久久久久久久免| 国产 一区精品| 久久6这里有精品| h日本视频在线播放| 好男人视频免费观看在线| 免费观看在线日韩| 毛片女人毛片| 神马国产精品三级电影在线观看| 日韩一区二区三区影片| freevideosex欧美| 香蕉精品网在线| 精品国产三级普通话版| 免费少妇av软件| 精品久久国产蜜桃| 欧美+日韩+精品| 嫩草影院新地址| 成人综合一区亚洲| 中文精品一卡2卡3卡4更新| 真实男女啪啪啪动态图| 男女国产视频网站| 国产午夜精品一二区理论片| 狂野欧美激情性xxxx在线观看| 黄色欧美视频在线观看| 99久久中文字幕三级久久日本| 久久久久国产精品人妻一区二区| 日韩人妻高清精品专区| 一本一本综合久久| 精品国产露脸久久av麻豆| 亚洲四区av| 一级毛片黄色毛片免费观看视频| 日日摸夜夜添夜夜爱| 一个人看的www免费观看视频| 亚洲精品久久午夜乱码| 欧美高清性xxxxhd video| 亚洲综合色惰| 久久精品国产a三级三级三级| 一本色道久久久久久精品综合| 99久久人妻综合| 国产精品一及| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产精品专区欧美| 99热这里只有是精品50| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 中文字幕亚洲精品专区| 日韩人妻高清精品专区| 国产69精品久久久久777片| 久久久久精品久久久久真实原创| 亚洲欧美成人精品一区二区| 国精品久久久久久国模美| 2022亚洲国产成人精品| 午夜精品国产一区二区电影 | 韩国av在线不卡| 视频区图区小说| 亚洲国产精品999| 六月丁香七月|