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

    多排環(huán)境下轉(zhuǎn)子葉片氣動彈性穩(wěn)定性機理分析

    2017-02-15 00:48:02李振鵬
    振動與沖擊 2017年1期
    關(guān)鍵詞:氣動彈性靜子導葉

    楊 慧, 李振鵬

    (北京航空航天大學 能源與動力工程學院, 先進航空發(fā)動機協(xié)同創(chuàng)新中心,北京 100191)

    多排環(huán)境下轉(zhuǎn)子葉片氣動彈性穩(wěn)定性機理分析

    楊 慧, 李振鵬

    (北京航空航天大學 能源與動力工程學院, 先進航空發(fā)動機協(xié)同創(chuàng)新中心,北京 100191)

    傳統(tǒng)葉片顫振分析多是基于單轉(zhuǎn)子研究模型,發(fā)動機的緊湊性要求導致級間距減小,多排耦合作用對顫振的影響將不容忽視。采用自行開發(fā)的程序?qū)δ承?.5級高壓壓氣機進行了流固耦合數(shù)值模擬,分析上、下游葉排對轉(zhuǎn)子葉片顫振特性的影響。針對典型工況,分別進行了單轉(zhuǎn)子模型,導葉轉(zhuǎn)子模型,轉(zhuǎn)子靜子模型,導葉轉(zhuǎn)子靜子模型的葉片氣動彈性穩(wěn)定性分析。研究表明,激波振蕩對顫振特性影響顯著;多排環(huán)境下存在非定常壓力波的反射和疊加,明顯改變轉(zhuǎn)子葉片表面的非定常壓力幅值和相位,進而改變轉(zhuǎn)子葉片氣動彈性穩(wěn)定性。多排干涉作用提高了轉(zhuǎn)子葉片的氣動阻尼,尤其是上、下游葉排同時作用時阻尼提高了近732.7%。

    顫振;全環(huán)多排;流固耦合;氣動阻尼;壓力波

    顫振是葉輪機械葉片振動的一種,嚴重制約著發(fā)動機的結(jié)構(gòu)完整性和可靠性[1]。當今對高推重比的追求[2]與輕薄葉片的使用使得葉片顫振問題越發(fā)劇烈,并成為設計人員不得不攻克的技術(shù)難題。顫振是流體誘導葉片振動的一種,涉及復雜的流固耦合作用,因其多發(fā)性及災難性,吸引了眾多學者對其進行深入研究。

    準確分析葉片顫振問題需要求解非定常流場和結(jié)構(gòu)動力學方程,為了節(jié)省數(shù)值模擬時間,達到工程應用的要求,需要根據(jù)特定的問題,使用簡化的數(shù)值模型和縮減的空間模型以減少求解方程的時間。SRIVASTAVA等[3]使用單通道模型施加相移邊界條件,進行了正激波對葉片顫振影響的無粘分析,指出激波對葉片氣彈的穩(wěn)定作用或失穩(wěn)作用與激波位置、葉片間相位角(InterBlade Phase Angle, IBPA)和振動頻率有關(guān)。VASANTHAKUMAR[4]使用頻域線化法,在單通道模型上施加相移邊條,指出激波位置及運動,以及激波與IBPA的關(guān)系對葉片氣動阻尼影響顯著,振動模態(tài)是重要的影響因素之一。但使用小擾動假設的線化方法分析非定常流場,不能很好地解決激波、分離流等強非線性問題。

    ISOMURA等[5]采用準三維的粘性非定常CFD程序,在單葉片模型上證明了激波振蕩而非葉片失速是跨音速風扇顫振失穩(wěn)的主要原因。張小偉等[6-7]采用能量法發(fā)現(xiàn)葉片間相位角和振動模態(tài)對葉片顫振有關(guān)鍵性影響。這些基于能量法的氣動彈性模型,不能考慮顫振問題中流體對固體的作用,也無法真實地捕捉流場的非定常流動特征。IM等[8]通過流固耦合方法,在半環(huán)葉排模型上研究發(fā)現(xiàn)脫體激波和葉尖間隙泄漏渦及吸力面附面層的非定常流動干擾作用是高負荷跨音速風扇顫振失穩(wěn)的主要原因。這些基于單轉(zhuǎn)子模型的研究成果,主要考慮同排葉片間的氣動耦合作用。

    發(fā)動機緊湊型的設計要求導致葉輪機部件轉(zhuǎn)、靜子間距越來越小,轉(zhuǎn)靜干涉作用明顯增強。目前多排干涉作用的研究重點在于其對氣動性能、非定常效應以及時序效應的影響[9-11],而多排葉片間的氣動耦合作用是否會對葉片氣動彈性穩(wěn)定性產(chǎn)生明顯作用,其影響機理如何,以及如何有效的預測多排環(huán)境中的顫振問題,都需要深入研究。BUFFUM等[12-14]使用頻域法,在二維單通道模型上,采用耦合模態(tài)方法首先分析了多排間的耦合作用。研究表明在多排干涉作用下,氣動阻尼與孤立振蕩葉柵顯著不同,轉(zhuǎn)靜軸向間距是關(guān)鍵參數(shù)。LI等[15]采用頻域法,在三維單通道模型上研究了軸向間距對轉(zhuǎn)子葉片氣動阻尼的影響。在特定間距值下,靜子的干涉作用會使轉(zhuǎn)子葉片的氣動阻尼增長100%,且葉片氣動阻尼和間距的關(guān)系與靜子葉片數(shù)有關(guān),并存在一個可以顯著改善顫振和強迫響應問題的最佳間距值。HSU等[16]為了考察使用多排模型和復雜氣動彈性模型的必要性,對比了全環(huán)多排流固耦合法和傳統(tǒng)單排單通道能量法的顫振特性計算結(jié)果,發(fā)現(xiàn)單排單通道模型不能正確預測顫振發(fā)生,而全環(huán)多排模型的預測結(jié)果和實驗一致。

    綜上可知,轉(zhuǎn)靜干涉作用對多排環(huán)境中的葉片顫振特性產(chǎn)生明顯影響,其影響機理需要深入研究。只有采用全環(huán)多排模型才能精確預測多排效應對葉片顫振的影響以及非定常壓力波在各個方向上的變化;考慮流固耦合作用才能真實地捕捉流場的非定常流動特征,特別是跨音速工況下。本文采用全環(huán)多排模型的流固耦合算法進行葉片顫振特性分析。針對典型氣動工況,通過單轉(zhuǎn)子模型,導葉-轉(zhuǎn)子模型,轉(zhuǎn)子-靜子模型,導葉-轉(zhuǎn)子-靜子模型的顫振特性結(jié)果對比,研究多排干涉作用對轉(zhuǎn)子葉片氣彈穩(wěn)定性的影響,為更深入地理解多排環(huán)境中葉片顫振發(fā)作機理提供參考。

    1 數(shù)值方法

    本文采用自行開發(fā)的流固耦合分析程序HGAE,通過葉片在非定常流場中隨時間的位移變化歷程來判斷葉片的氣動彈性穩(wěn)定性。

    流體域求解動邊界和動態(tài)變形網(wǎng)格下的守恒型積分形式的三維非定??蓧嚎s雷諾平均Navier-Stokes方程,其方程為:

    (1)

    葉片的動力學方程通過廣義坐標解耦,得到N個線性無關(guān)的常微分方程組:

    (2)

    模態(tài)位移向量ζ與物理位移向量x的關(guān)系為

    x=Φζ

    (3)

    式中:Φ為質(zhì)量歸一化的系統(tǒng)振型矩陣,由單葉片振型組集[18]而來。通過迭代求解方程(1)和方程(2),得到每一時間步的流場參數(shù)及葉片位移,從而得到模態(tài)位移隨時間變化的歷程。

    采用指數(shù)函數(shù)曲線擬合模態(tài)位移曲線,計算各階模態(tài)的氣動阻尼,判定葉片的氣動彈性穩(wěn)定性,該流固耦合求解方法針對全環(huán)模型可同時考慮多個振動模態(tài)和所有葉片間相位角,詳細的數(shù)值方法介紹見參考文獻 [18-19] ,其在氣動彈性領(lǐng)域的有效性驗證見文獻[17,20]。

    2 計算說明

    本文采用某型高壓壓氣機前1.5級,算例的基本參數(shù)見表1。

    表1 基本參數(shù)

    圖1(左)為導葉-轉(zhuǎn)子-靜子模型的全環(huán)網(wǎng)格,其拓撲結(jié)構(gòu)均為O4H型。導葉的全環(huán)網(wǎng)格節(jié)點數(shù)約為780×104;轉(zhuǎn)子的全環(huán)網(wǎng)格節(jié)點數(shù)約為580×104,靜子的全環(huán)網(wǎng)格節(jié)點數(shù)約為523×104。第1層網(wǎng)格距壁面厚度均為0.01 mm。圖1(右)為轉(zhuǎn)子葉片的有限元網(wǎng)格,采用八節(jié)點六面體單元,根部固支,網(wǎng)格數(shù)為24×2×20。

    圖1 導葉-轉(zhuǎn)子-靜子模型網(wǎng)格(左)和葉片有限元網(wǎng)格(右)Fig.1 Grids of IGV-rotor-stator(left)and blade FEM grids (right)

    為確定葉片氣動彈性分析的工況點,通過1.5級模型完成特性計算。進口給定總壓101 325 Pa,總溫288.15 K,軸向進氣,出口給定輪轂處靜壓并結(jié)合徑向壓力平衡方程,轉(zhuǎn)靜交界面使用混合面[21]方法。非定常流動計算時,轉(zhuǎn)靜交界面使用雙線性插值的滑移面[22]方法。

    顫振計算工況選取近設計點,流場較為“干凈”,便于分析上、下游葉排的干涉作用對轉(zhuǎn)子葉片氣動彈性穩(wěn)定性的影響。分別針對單轉(zhuǎn)子模型,導葉-轉(zhuǎn)子模型,轉(zhuǎn)子-靜子模型,導葉-轉(zhuǎn)子-靜子模型進行顫振特性分析,各模型的邊界條件由特性計算時交界面的參數(shù)分布給出,從而保證各個模型下計算流量以及葉表的定常氣動負荷一致,進而實現(xiàn)計算工況一致。

    葉片顫振計算在全環(huán)定常流動結(jié)果上給定葉片微小的模態(tài)速度以激勵葉片初始振動。為了簡化分析,計算中僅考慮葉片的1階彎曲振動模態(tài),包含所有行波節(jié)徑(如28個葉片對應28個行波節(jié)徑,分別為0,1,-1,2,-2,…,13,-13,14)。計算物理時間步長取為葉片振動周期的1/100,共計算了20個周期。

    3 計算結(jié)果

    3.1 氣動阻尼計算結(jié)果

    圖2所示為四種模型下轉(zhuǎn)子葉片氣動阻尼隨節(jié)徑的變化。對于單轉(zhuǎn)子模型,氣動阻尼隨節(jié)徑正弦規(guī)律變化,節(jié)徑對氣動阻尼影響顯著。這與文獻[4]和文獻[20]的計算結(jié)果一致。氣動彈性最不穩(wěn)定狀態(tài)(最小值)對應前行波1節(jié)徑,氣動阻尼為負,葉片氣動彈性失穩(wěn)。

    導葉沒有使最危險節(jié)徑發(fā)生變化,下游靜子的作用使得最危險節(jié)徑由1變?yōu)?,上下游葉排同時存在時,危險節(jié)徑由1變?yōu)?。由圖可知,葉排干涉作用使得氣動阻尼隨節(jié)徑振蕩變化,在單轉(zhuǎn)子正弦規(guī)律下呈現(xiàn)多峰值現(xiàn)象,增加或降低了單轉(zhuǎn)子的氣動阻尼值。表明非定常壓力波沿葉排周向的傳播受節(jié)徑的影響,如果使用單通道模型必須遍歷所有的節(jié)徑,才能獲得正確的結(jié)果。

    圖2 不同模型下氣動阻尼對比Fig.2 Aerodynamic damping comparision under different models

    圖3所示為不同模型的最小氣動阻尼值對比,由圖可知,導葉使得轉(zhuǎn)子最危險節(jié)徑下的氣動阻尼增加了174.8%,下游靜子使其提高近327.8%,上下游葉排同時存在時,最小氣動阻尼提高近732.7%。在本文選取的級間距值下,多排環(huán)境顯著提高了最不穩(wěn)定狀態(tài)下轉(zhuǎn)子葉片的氣彈穩(wěn)定性,尤其是上下游葉排同時存在的情況下,并且下游靜子對氣彈穩(wěn)定性的影響大于上游導葉。

    圖3 最小氣動阻尼對比 Fig.3 Minimum aerodynamic damping comparison

    3.2 單轉(zhuǎn)子模型結(jié)果分析

    圖4為85%葉高和葉尖區(qū)域的流場Ma分布。由圖可知,在轉(zhuǎn)子通道中,存在一道較強的脫體激波,與吸力面相交于50%弦長處;在葉尖區(qū)域,泄漏渦與激波相互作用,引起局部氣流堵塞。

    圖4 流場Ma云圖(85%葉高)Fig.4 Ma contour of Flow field (85% span)

    圖5為非定常壓力的時空分布圖,非定常壓力定義為在葉片一個振動周期內(nèi),瞬時壓力和平均壓力之差。

    圖5 非定常壓力時空分布(單轉(zhuǎn)子,85%葉高)Fig.5 Time-space distribution of unsteady pressure(single rotor, 85% span)

    由圖5可知,吸力面(SS)非定常壓力強于壓力面(PS),在轉(zhuǎn)子葉片吸力面50%弦長處存在激波附面層干擾現(xiàn)象,激波振蕩是引起葉表非定常壓力的主要原因。

    葉尖位置葉片振動幅值最大,引起的能量轉(zhuǎn)換也最劇烈,對葉片氣彈穩(wěn)定性影響明顯,而本文脫體激波與葉尖泄漏渦在此處又存在干涉作用(如圖4(b))。分析葉尖區(qū)域A點(激波附近)、B點和C點(葉尖其他區(qū)域)的非定常壓力和葉片位移的相位關(guān)系(如圖6)可知,激波附近A點的壓力和葉片位移反相變化,表明激波對氣彈穩(wěn)定性起穩(wěn)定作用[23];而在葉尖其他區(qū)域(B點和C點)壓力和葉片位移相位近似同相變化,產(chǎn)生氣彈失穩(wěn)效應。由此推斷在本文算例中,激波振蕩通過改變轉(zhuǎn)子葉尖非定常壓力的相位,對氣動彈性起到一定穩(wěn)定作用

    圖6 轉(zhuǎn)子吸力面壓力云圖及點A、B、C的壓力和位移時間歷程Fig.6 Pressure contour on rotor suction surface and time histories of pressure and displacement at point A, B and C

    3.3 導葉-轉(zhuǎn)子模型結(jié)果分析

    圖7為上游導葉軸向力的頻譜分析,由圖可知,多排環(huán)境下導葉上的擾動頻率主要為葉片振動頻率和轉(zhuǎn)子通過頻率,在上游導葉中,振動頻率的作用甚至高于通過頻率??梢娤掠无D(zhuǎn)子葉片振動引起的非定常壓力波向上游傳播,作用在前排導葉上。圖8為轉(zhuǎn)子葉片吸力面非定常壓力幅值和相位,通過單轉(zhuǎn)子模型和導葉-轉(zhuǎn)子模型的葉表非定常壓力對比可知,導葉-轉(zhuǎn)子模型激波前非定常高壓力幅值區(qū)域沿徑向擴展到近50%葉高,且葉尖前緣高幅值范圍也變大;激波前的相位產(chǎn)生近似180°的變化,激波對導葉-轉(zhuǎn)子模型的顫振穩(wěn)定性作用明顯。結(jié)果表明導葉將轉(zhuǎn)子葉片振動產(chǎn)生的壓力擾動波反射回轉(zhuǎn)子表面,并與轉(zhuǎn)子葉片上的非定常壓力波疊加。文獻[23]研究表明,反射波與葉表壓力擾動波的疊加情況由兩者之間的相位差決定,而此相位差受轉(zhuǎn)子葉片的振動節(jié)徑影響。圖2中氣動阻尼值隨節(jié)徑的變化規(guī)律,驗證了文獻[23]的這一結(jié)論。

    圖7 導葉軸向力(Fx)頻譜分析Fig.7 Frequency spectrum of axial force (Fx) on IGV

    圖8 轉(zhuǎn)子葉片非定常壓力幅值及相位(吸力面,振動頻率)Fig.8Unsteady pressure amplitude and phase on rotor blade (suction surface, vibration frequency)

    為了分析葉片振動對流場的影響,針對導葉-轉(zhuǎn)子模型分別進行了非定常流動計算和葉片振動的顫振特性計算。圖9列出了轉(zhuǎn)子葉片85%葉高處葉表非定常壓力的時空分布圖。非定常流動計算時,勢干擾作用引起葉片前緣劇烈的壓力脈動,流場呈現(xiàn)周期性變化;當考慮葉片振動的氣動彈性穩(wěn)定性分析時,葉片的衰減振動導致流場勢干擾引起的周期性變化明顯減弱,表現(xiàn)為壓力脈動的幅值明顯降低。而在吸力面50%弦長處出現(xiàn)了激波振蕩引起的周期性變化。

    圖10為兩種計算條件下導葉和轉(zhuǎn)子葉片85%葉高處通過頻率下的非定常壓力幅值對比。可見葉片的振動作用顯著降低了上下游通過頻率下葉表非定常壓力幅值。多排環(huán)境中,葉片振動會對強迫響應產(chǎn)生影響。

    圖9 轉(zhuǎn)子葉表非定常壓力時空分布(導葉轉(zhuǎn)子,85%葉高)Fig.9 Time-space distribution of unsteady pressure on rotor surface (IGV-rotor, 85% span)

    圖10 兩種計算下非定常壓力幅值對比(85%葉高,通過頻率)Fig.10 Comparison of two calculated amplitudes of unsteady pressure (85% span, blade passing frequency)

    3.4 轉(zhuǎn)子-靜子模型結(jié)果分析

    圖11為轉(zhuǎn)子-靜子模型中葉片表面振動頻率下非定常壓力幅值和相位。與單轉(zhuǎn)子(圖8(a))結(jié)果對比,激波后高幅值范圍增加了約20%弦長,且在尾緣和葉尖前緣區(qū)域幅值均增加;葉尖尾緣大面積區(qū)域相位發(fā)生近120°變化,激波后50%葉高以下相位發(fā)生了180°轉(zhuǎn)變。這些變化使得轉(zhuǎn)子葉片最小氣動阻尼提高(圖3),且最危險節(jié)徑發(fā)生改變(圖2)。

    圖11 轉(zhuǎn)子葉片壓力幅值及相位(轉(zhuǎn)子-靜子模型,吸力面,振動頻率)Fig.11 Unsteady pressure amplitude and phase on rotor blade (rotor-stator, suc.sur, vibration freq.)

    如圖12所示,在靜子葉片尾緣處存在大面積分離,分離頻率約為1 516 Hz(如圖12(b))。周期性分離是流場主要的擾動源,由圖13中軸向力的頻譜分析可知,即使在轉(zhuǎn)子域中,流場分離仍是主要擾動源,振動頻率和通過頻率的作用十分微小。在靜子域中并沒有振動頻率的作用,表明流場的大分離影響了壓力擾動波的傳播。下游靜子大分離流場對壓力波反射的影響改變了多排環(huán)境下葉片氣彈穩(wěn)定性。

    圖12 靜子流場分離區(qū)(85%葉高)Fig.12 Separation flow in the flow field of stator

    圖13 軸向力(Fx)頻譜分析Fig.13 Frequency spectrum of axial force (Fx)

    3.5 導葉-轉(zhuǎn)子-靜子模型結(jié)果分析

    導葉-轉(zhuǎn)子-靜子模型的轉(zhuǎn)子葉表非定常壓力幅值和相位如圖14所示。與單轉(zhuǎn)子模型(圖8(a))對比,葉表激波位置,前緣和尾緣壓力幅值均增加明顯,多排的耦合作用劇烈;在激波振蕩位置及激波前大面積區(qū)域相位接近0°,使得原來葉尖前緣及激波處的相位發(fā)生了近180°變化,在激波后靠近尾緣大面積區(qū)域相位同樣發(fā)生了近180°變化。從整體上看,導葉-轉(zhuǎn)子-靜子模型的轉(zhuǎn)子葉表非定常壓力幅值和相位近似為導葉-轉(zhuǎn)子和轉(zhuǎn)子-靜子兩種模型的耦合效果,導致轉(zhuǎn)子葉片最小氣動阻尼較單轉(zhuǎn)子模型大幅提高(圖3)。

    圖14 轉(zhuǎn)子葉片壓力幅值及相位(導葉-轉(zhuǎn)子-靜子模型,吸力面,振動頻率)Fig.14 Unsteady pressure amplitude and phase on rotor blade (IGV-rotor-stator, suc.sur., vibration freq.)

    4 結(jié) 論

    本文采用自行開發(fā)的流固耦合程序詳細分析了多排耦合作用對轉(zhuǎn)子葉片顫振特性的影響。研究了上游導葉、下游靜子及上下游葉排同時存在時對轉(zhuǎn)子葉片一階彎曲模態(tài)氣動彈性穩(wěn)定性的影響機理,有如下結(jié)論:

    (1)對于單轉(zhuǎn)子模型,激波振蕩是葉表非定常壓力的主要原因;其通過改變?nèi)~尖非定常壓力的相位,對葉片氣動彈性起到一定穩(wěn)定作用。

    (2)葉片振動作用,將導致上下游葉排通過頻率下的非定常壓力幅值降低。

    (3)多葉排環(huán)境下,導葉、靜子作用使得轉(zhuǎn)子葉表非定常壓力波反射,顯著改變轉(zhuǎn)子葉表非定常壓力的幅值和相位,使得氣動阻尼隨節(jié)徑呈現(xiàn)多峰值振蕩現(xiàn)象。本文算例中,多排的干涉作用明顯改善了葉片一階彎曲振動最不穩(wěn)定狀態(tài)的氣彈穩(wěn)定性,下游靜子的作用大于上游導葉。

    (4)由本文研究結(jié)果可以推斷,使用單轉(zhuǎn)子模型分析轉(zhuǎn)子葉片的氣動彈性穩(wěn)定性將得到不準確的預測結(jié)果,特別是跨音速工況下,同時存在復雜的分離流動和葉尖間隙渦的情況。

    (5)本算例中,單轉(zhuǎn)子模型給出過保守的氣動彈性穩(wěn)定性預測。下游靜子的周期性分離流動會影響非定常壓力波的反射,其影響機制需要進一步深入研究。

    [ 1 ] SRINIVASAN A V. Flutter and resonant vibration characteristics of engine blades[J]. Journal of Engineering for Gas Turbines and Power, 1997, 119(4): 742-775.

    [ 2 ] 陳懋章. 風扇/壓氣機技術(shù)發(fā)展和對今后工作的建議[J]. 航空動力學報, 2002, 17(1): 1-15. CHEN Maozhang. Development of fan/compressor techniques and suggestions on further researches[J]. Journal of Aerospace Power,2002,17(1):1-15.

    [ 3 ] SRIVASTAVA R, KEITH T G. Shock induced flutter of turbomachinery blade row[C]//ASME Turbo Expo 2004: Power for Land, Sea, and Air.American Society of Mechanical Engineers, 2004:487-496.

    [ 4 ] VASANTHAKUMAR P. Computation of aerodynamic damping for flutter analysis of a transonic fan[C]//ASME 2011 Turbo Expo: Turbine Technical Conference and Exposition. American Society of Mechanical Engineers, 2011:1429-1437.

    [ 5 ] ISOMURA K, GILES M B. A numerical study of flutter in a transonic fan[J]. Journal of Turbomachinery,1998,120(3): 500-507.

    [ 6 ] 張小偉,王延榮,許可寧.葉輪機械葉片顫振的影響參數(shù)[J].航空動力學報, 2011,26(7):1557-1561. ZHANG Xiaowei,WANG Yanrong,XU Kening.Effects of parameters on blade flutter in turbomachinery [J].Journal of Aerospace Power, 2011, 26(7):1557-1561.

    [ 7 ] 張小偉,王延榮.葉片間相位角對葉片顫振的影響[J].航空動力學報, 2010,25(2):412-416. ZHANG Xiaowei,WANG Yanrong.Influence of interblade phase angle on the flutter of rotor blades[J].Journal of Aerospace Power, 2010, 25(2):412-416.

    [ 8 ] IM H S, ZHA G C.Flutter prediction of a transonic fan with travelling wave using fully coupled fluid/structure interaction[C]//ASME Turbo Expo 2013:Turbine Technical Conference and Exposition. American Society of Mechanical Engineers,2013.V07BT33A003.

    [ 9 ] SUN Y, REN Y X,FU S,et al. The influence of rotor-stator spacing on the loss in one-stage transonic compressor[C]//ASME Turbo Expo 2009: Power for land, Sea and Air . American Society of Mechanical Engineers, 2009:1707-1715.

    [10] NG W F, OBRIEN W F, OLSEN T L. Experimental investigation of unsteady fan flow interaction with downstream struts[J]. Journal of Propulsion and Power, 1987, 3(2): 157-163.

    [11] SAREN V E, SAVIN N M, DORNEY D J, et al. Experimental and numerical investigation of airfoil clocking and inner-blade-row gap effects on axial compressor performance [J]. International Journal of Turbo and Jet Engines, 1998(15): 235-252.

    [12] BUFFURN D H. Blade row interaction effects on flutter and forced response[J]. Journal of Propulsion and Power, 1995, 11(2): 205-212.

    [13] HALL K C, SILKOWSKI P D.The influence of neighboring blade rows on the unsteady aerodynamic response of cascades[J]. Journal of Turbomachinery, 1997(119):85-93.

    [14] SILKOWSKI P D, HALL K C. A coupled mode analysis of unsteady multistage flows in turbomachinery[J]. Journal of Turbomachinery,1998(120): 410-421.

    [15] LI H D,HE L.Blade aerodynamic damping variation with rotor-stator gap:A computational study using single-passage approach[J]. Journal of Turbomachinery, 2005,127(3): 573-579.

    [16] HSU K,HOYNIAK D,ANAND M S.Full-annulus multi-row flutter analyses[C]//ASME Turbo Expo 2012: Turbine Technical Conference and Exposition.American Society of Mechanical Engineers, 2012:1453-1462.

    [17] 鄭赟.基于非結(jié)構(gòu)網(wǎng)格的氣動彈性數(shù)值方法研究[J]. 航空動力學報, 2009,24(9):2069-2077. ZHENG Yun.Computational aeroelasticity with an unstructured grid method[J].Journal of Aerospace Power, 2009, 24(9): 2069-2077.

    [18] 張錦,劉曉平.葉輪機振動模態(tài)分析理論及數(shù)值方法[M].北京:國防工業(yè)出版社, 2001.

    [19] 鄭赟,楊慧.跨音速風扇全環(huán)葉片顫振特性的流固耦合分析[J].北京航空航天大學學報,2013,39(5):626-630. ZHENG Yun,YANG Hui.Full assembly fluid/structured flutter analysis of a transonic fan[J].Journal of Beijing University of Aeronautics and Astronautics, 2013,39(5):626-630.

    [20] ZHENG Yun, YANG Hui. Coupled fluid-structure flutter analysis of a transonic fan[J]. Chinese Journal of Aeronautics, 2011, 24(3): 258-264.

    [21] CHIMA R V.Calculation of multistage turbomachinery using steady characteristic boundary conditions[R].Reston,VA:36th Aerospace Sciences Meeting &Exhibit,1998.

    [22] MATHUR S R.Unsteady flow simulations using unstruck tured sliding meshes[R]. Reston, VA:25th AIAA Fluid Dynamics Conference, 1994.

    [23] VAHDATI M, SIMPSON G, IMTEGUN M.Mechanisms for wide-chord fan blade flutter[J]. Journal of Turbomachinery, 2011, 133(4): 1396-1402.

    Aeroelastic stability analysis of rotor blades under multi-row environment

    YANG Hui, LI Zhenpeng

    (Collaborative Innovation Center for Advanced Aero-Engine,School of Energy and Power Engineering, Beihang University, Beijing 100191, China)

    Conventional blade flutter analysis is normally based on an isolated blade row model, the influence of multi-row aerodynamic coupling on blade flutter characteristics can’t be ignored when rotor-stator gaps decrease due to aeroengine compact requirements. A fluid-structure coupled simulation for a 1.5stage HPC was conducted with a self-developed algorithm to analyze the influence of upstream and downstream blade rows on rotor blade flutter characteristics. Aiming at a typical operation condition, rotor blades’ aeroelastic stability analyses were performed with an isolated rotor model, an IGV-rotor model, a rotor-stator one and an IGV-rotor-stator one, respectively. The results showed that the shock wave vibration influences the flutter stability significantly; there are reflection and superposition of unsteady pressure waves under the multi-row environment, the amplitude and phase of unsteady pressures on the rotor blade surface are changed obviously and furthermore the blade aeroelastic stability is changed; multi-row interferences raise the aerodynamic damping of rotor blade, especially, when the upstream and downstream blade rows act simultaneously, the damping value increases by nearly 732.7%.

    blade flutter; full-annulus/multi-row; fluid-structure interaction; aerodynamic damping; pressure wave

    2015-07-10 修改稿收到日期:2015-12-29

    楊慧 女,博士,講師,1970年10月生

    V232.4

    A

    10.13465/j.cnki.jvs.2017.01.022

    猜你喜歡
    氣動彈性靜子導葉
    壓氣機緊湊S形過渡段內(nèi)周向彎靜子性能數(shù)值計算
    何必喧囂慰寂寥
    河北畫報(2020年22期)2020-07-12 12:51:08
    堅持了十年的書信
    山鄉(xiāng)一瞥
    壓氣機導葉用GTD222合金鑄造技術(shù)研究
    模具制造(2019年4期)2019-06-24 03:36:50
    飛翼無人機嗡鳴氣動彈性響應分析
    模態(tài)選取對靜氣動彈性分析的影響
    航空學報(2015年4期)2015-05-07 06:43:37
    直升機的氣動彈性問題
    航空學報(2015年4期)2015-05-07 06:43:25
    大型風力機整機氣動彈性響應計算
    混流式模型水輪機導葉水力矩試驗
    少妇 在线观看| 亚洲美女搞黄在线观看| 国产亚洲5aaaaa淫片| 亚洲人成网站在线观看播放| 欧美成人一区二区免费高清观看| h日本视频在线播放| 色播亚洲综合网| 日韩一本色道免费dvd| 老师上课跳d突然被开到最大视频| 色吧在线观看| 日韩 亚洲 欧美在线| 国产一区亚洲一区在线观看| 国产熟女欧美一区二区| 欧美97在线视频| 久久国内精品自在自线图片| 亚洲国产欧美人成| 国产黄频视频在线观看| 国产精品一区二区三区四区免费观看| 日韩视频在线欧美| 欧美日韩视频高清一区二区三区二| 边亲边吃奶的免费视频| 十八禁网站网址无遮挡 | a级毛片免费高清观看在线播放| 国产又色又爽无遮挡免| 日日摸夜夜添夜夜爱| 欧美日韩精品成人综合77777| 国产精品人妻久久久久久| 一级毛片久久久久久久久女| 亚洲精品日韩在线中文字幕| 亚洲丝袜综合中文字幕| 亚洲国产最新在线播放| 国产精品偷伦视频观看了| 小蜜桃在线观看免费完整版高清| 成年人午夜在线观看视频| 一边亲一边摸免费视频| av免费观看日本| 亚洲欧美日韩无卡精品| 久久精品国产亚洲av涩爱| 九九爱精品视频在线观看| 在现免费观看毛片| 3wmmmm亚洲av在线观看| 亚洲最大成人中文| 国产91av在线免费观看| 日本爱情动作片www.在线观看| av网站免费在线观看视频| 国产一级毛片在线| 久久久精品欧美日韩精品| 1000部很黄的大片| 国产精品一区www在线观看| 亚洲精品色激情综合| 熟女电影av网| 中文字幕久久专区| 麻豆成人午夜福利视频| 色哟哟·www| 插逼视频在线观看| 国产成人aa在线观看| 少妇人妻久久综合中文| 亚洲在线观看片| 我的女老师完整版在线观看| 亚洲人成网站在线播| 免费大片18禁| 国产一区二区在线观看日韩| 白带黄色成豆腐渣| 中文字幕免费在线视频6| 性色av一级| 亚洲内射少妇av| 菩萨蛮人人尽说江南好唐韦庄| 十八禁网站网址无遮挡 | a级一级毛片免费在线观看| 亚洲精品456在线播放app| 高清欧美精品videossex| 精品一区在线观看国产| 精品国产三级普通话版| 国产免费一级a男人的天堂| 久久精品国产亚洲av涩爱| 精品视频人人做人人爽| 狠狠精品人妻久久久久久综合| 精品人妻熟女av久视频| 国产成人免费无遮挡视频| 亚州av有码| 午夜福利在线在线| 国产片特级美女逼逼视频| 欧美三级亚洲精品| 欧美高清性xxxxhd video| 精品久久久久久久末码| 国产午夜精品久久久久久一区二区三区| 午夜免费男女啪啪视频观看| 别揉我奶头 嗯啊视频| 国产美女午夜福利| av女优亚洲男人天堂| 日韩伦理黄色片| 亚洲第一区二区三区不卡| 国产一区有黄有色的免费视频| 午夜免费鲁丝| 99久久精品热视频| 青春草视频在线免费观看| 毛片一级片免费看久久久久| 日本午夜av视频| 亚洲av中文字字幕乱码综合| 最近2019中文字幕mv第一页| 国产成人a∨麻豆精品| 舔av片在线| 国产男人的电影天堂91| 26uuu在线亚洲综合色| 国产亚洲5aaaaa淫片| 久久久久国产网址| 精品久久久噜噜| 久久久久久久国产电影| 日本黄色片子视频| 国产精品一区www在线观看| 亚洲精品国产av成人精品| av国产精品久久久久影院| 99久久中文字幕三级久久日本| 午夜福利在线在线| 成年女人看的毛片在线观看| 亚洲人成网站在线观看播放| 欧美日韩亚洲高清精品| 一级毛片我不卡| 2018国产大陆天天弄谢| 男女下面进入的视频免费午夜| 久久亚洲国产成人精品v| 国产探花在线观看一区二区| 免费黄色在线免费观看| 97热精品久久久久久| 精品久久久精品久久久| 亚洲精品久久午夜乱码| 国产精品秋霞免费鲁丝片| 日产精品乱码卡一卡2卡三| 寂寞人妻少妇视频99o| 国产免费又黄又爽又色| 久久久久久久大尺度免费视频| 六月丁香七月| 又粗又硬又长又爽又黄的视频| 2022亚洲国产成人精品| 色5月婷婷丁香| 中文字幕免费在线视频6| 最近最新中文字幕免费大全7| 大话2 男鬼变身卡| 香蕉精品网在线| 插阴视频在线观看视频| 丰满乱子伦码专区| 亚洲一区二区三区欧美精品 | 欧美精品一区二区大全| 网址你懂的国产日韩在线| 国产欧美亚洲国产| 国产成人freesex在线| 国产伦理片在线播放av一区| 97超碰精品成人国产| 日韩一区二区视频免费看| 伊人久久国产一区二区| 99久久人妻综合| 午夜精品国产一区二区电影 | 老女人水多毛片| 精品99又大又爽又粗少妇毛片| 国产精品国产三级国产专区5o| 亚洲综合精品二区| 亚洲,一卡二卡三卡| 午夜激情久久久久久久| 日本色播在线视频| 天堂中文最新版在线下载 | 18+在线观看网站| 三级国产精品片| 久久久久精品性色| 一级黄片播放器| 少妇人妻精品综合一区二区| 成年人午夜在线观看视频| 久久久精品欧美日韩精品| 噜噜噜噜噜久久久久久91| 久久精品久久精品一区二区三区| 欧美少妇被猛烈插入视频| 久久久久久久久大av| 欧美人与善性xxx| 熟女av电影| 成人一区二区视频在线观看| 特大巨黑吊av在线直播| 久久精品夜色国产| av在线播放精品| 一级毛片 在线播放| 国产成人福利小说| 久久精品人妻少妇| 日韩人妻高清精品专区| 深夜a级毛片| 久久国产乱子免费精品| 亚洲精品成人久久久久久| 亚洲av中文av极速乱| 日本wwww免费看| 神马国产精品三级电影在线观看| 99精国产麻豆久久婷婷| 国产熟女欧美一区二区| 汤姆久久久久久久影院中文字幕| 精品国产三级普通话版| 天堂中文最新版在线下载 | 亚洲精品中文字幕在线视频 | 韩国av在线不卡| 99九九线精品视频在线观看视频| 观看美女的网站| 夫妻午夜视频| 日本色播在线视频| 2022亚洲国产成人精品| 日韩欧美 国产精品| www.av在线官网国产| 大又大粗又爽又黄少妇毛片口| 亚洲图色成人| 亚洲人成网站高清观看| 国产高清三级在线| 国产毛片a区久久久久| 熟女电影av网| 亚洲精品日韩在线中文字幕| 看非洲黑人一级黄片| 又粗又硬又长又爽又黄的视频| 亚洲精品影视一区二区三区av| 91精品伊人久久大香线蕉| 国内精品美女久久久久久| 国产高潮美女av| 深夜a级毛片| 久久影院123| 99久久精品国产国产毛片| 国产成人aa在线观看| 久久精品国产亚洲av天美| 国产在线男女| 久久这里有精品视频免费| 日韩电影二区| 午夜福利在线观看免费完整高清在| 99热国产这里只有精品6| 国产欧美另类精品又又久久亚洲欧美| 大码成人一级视频| 好男人视频免费观看在线| 亚洲精品第二区| 成年免费大片在线观看| 久久精品国产鲁丝片午夜精品| 小蜜桃在线观看免费完整版高清| 99久久九九国产精品国产免费| 欧美亚洲 丝袜 人妻 在线| 欧美老熟妇乱子伦牲交| 麻豆成人av视频| 午夜免费观看性视频| 丝瓜视频免费看黄片| 国产v大片淫在线免费观看| 久久久久网色| 一二三四中文在线观看免费高清| 久久ye,这里只有精品| 搡老乐熟女国产| 国产亚洲最大av| 久久99热这里只有精品18| 亚洲av国产av综合av卡| 精品国产一区二区三区久久久樱花 | 观看美女的网站| 91久久精品电影网| 国产熟女欧美一区二区| 午夜免费观看性视频| 精品一区二区三卡| 人人妻人人爽人人添夜夜欢视频 | 在线看a的网站| 国产成人免费观看mmmm| 大片免费播放器 马上看| 啦啦啦啦在线视频资源| 水蜜桃什么品种好| 日本欧美国产在线视频| 免费看日本二区| 黄色配什么色好看| 2021天堂中文幕一二区在线观| 激情五月婷婷亚洲| 狂野欧美激情性bbbbbb| a级一级毛片免费在线观看| 午夜福利在线观看免费完整高清在| 18禁裸乳无遮挡动漫免费视频 | 免费看不卡的av| 国产av码专区亚洲av| 99热网站在线观看| 国产成人a∨麻豆精品| 少妇人妻久久综合中文| 国产亚洲91精品色在线| 国产亚洲av嫩草精品影院| 国产日韩欧美亚洲二区| 三级经典国产精品| 欧美bdsm另类| 国产成人精品婷婷| av国产久精品久网站免费入址| 成年av动漫网址| 亚洲色图av天堂| 蜜臀久久99精品久久宅男| 久久久久久久久大av| 免费观看无遮挡的男女| 国产精品熟女久久久久浪| 内地一区二区视频在线| 大又大粗又爽又黄少妇毛片口| 成人毛片60女人毛片免费| 日韩一本色道免费dvd| 久久久久久伊人网av| 中文精品一卡2卡3卡4更新| 视频中文字幕在线观看| 春色校园在线视频观看| 亚洲精品国产av成人精品| 亚洲av男天堂| av免费在线看不卡| 国内少妇人妻偷人精品xxx网站| 日韩av不卡免费在线播放| 99视频精品全部免费 在线| 久久久久久久久久成人| 18+在线观看网站| 久久精品国产亚洲av涩爱| 看黄色毛片网站| av在线app专区| 欧美97在线视频| 国产在线一区二区三区精| 99热6这里只有精品| 国产一级毛片在线| 午夜精品一区二区三区免费看| 日韩欧美一区视频在线观看 | 国产精品.久久久| 永久网站在线| 两个人的视频大全免费| av国产免费在线观看| 一本色道久久久久久精品综合| 精品久久久久久久久av| 一边亲一边摸免费视频| 最近中文字幕高清免费大全6| 亚洲va在线va天堂va国产| 韩国av在线不卡| 亚洲在线观看片| av在线蜜桃| 有码 亚洲区| 亚洲欧美一区二区三区国产| 欧美变态另类bdsm刘玥| 成年女人看的毛片在线观看| 久久久国产一区二区| 精品少妇久久久久久888优播| 亚洲va在线va天堂va国产| 亚洲精品成人av观看孕妇| 国产高清三级在线| 人人妻人人看人人澡| 久久久亚洲精品成人影院| 啦啦啦在线观看免费高清www| 国产免费一区二区三区四区乱码| 午夜精品国产一区二区电影 | 不卡视频在线观看欧美| 亚洲一区二区三区欧美精品 | 狂野欧美白嫩少妇大欣赏| 午夜免费男女啪啪视频观看| 日日啪夜夜撸| 久久精品久久精品一区二区三区| 中文字幕免费在线视频6| av福利片在线观看| 日韩大片免费观看网站| 国内精品宾馆在线| 免费黄频网站在线观看国产| 国产白丝娇喘喷水9色精品| 日本午夜av视频| 91久久精品国产一区二区三区| 一区二区三区免费毛片| 国产 精品1| 国产日韩欧美在线精品| 国产高清不卡午夜福利| 97人妻精品一区二区三区麻豆| 色视频在线一区二区三区| 午夜福利视频精品| 成年免费大片在线观看| 欧美日本视频| 午夜福利高清视频| 国产男女内射视频| 日本猛色少妇xxxxx猛交久久| 97热精品久久久久久| 秋霞在线观看毛片| 成人高潮视频无遮挡免费网站| 丝袜喷水一区| 国产精品一二三区在线看| 国产精品不卡视频一区二区| 国产爱豆传媒在线观看| 成人特级av手机在线观看| 亚洲va在线va天堂va国产| 国产亚洲av片在线观看秒播厂| 色网站视频免费| 欧美亚洲 丝袜 人妻 在线| 日韩大片免费观看网站| 搡老乐熟女国产| 三级国产精品片| 99久久精品一区二区三区| 久久久午夜欧美精品| av在线天堂中文字幕| 极品教师在线视频| 黄色配什么色好看| 高清日韩中文字幕在线| 国内精品宾馆在线| 午夜免费男女啪啪视频观看| 婷婷色综合大香蕉| 视频中文字幕在线观看| av在线app专区| 少妇高潮的动态图| 久久精品人妻少妇| 嫩草影院精品99| 青春草国产在线视频| kizo精华| 深夜a级毛片| 国产免费一区二区三区四区乱码| 欧美 日韩 精品 国产| 另类亚洲欧美激情| av在线观看视频网站免费| 69人妻影院| 乱码一卡2卡4卡精品| 亚洲第一区二区三区不卡| 麻豆国产97在线/欧美| 亚洲精品456在线播放app| 亚洲经典国产精华液单| 欧美另类一区| 蜜桃亚洲精品一区二区三区| 肉色欧美久久久久久久蜜桃 | 乱系列少妇在线播放| 99久久九九国产精品国产免费| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美一区二区三区国产| 天天一区二区日本电影三级| 黑人高潮一二区| 国产乱来视频区| 白带黄色成豆腐渣| 国国产精品蜜臀av免费| 国产日韩欧美亚洲二区| 一区二区三区精品91| 日韩免费高清中文字幕av| 国产毛片在线视频| 久久久久久久久大av| 97精品久久久久久久久久精品| 国产精品一区二区三区四区免费观看| 丰满人妻一区二区三区视频av| 精品久久久久久久末码| 国产探花在线观看一区二区| 久久久午夜欧美精品| 久久午夜福利片| 新久久久久国产一级毛片| 日本午夜av视频| 能在线免费看毛片的网站| 男女无遮挡免费网站观看| 国产精品偷伦视频观看了| av在线播放精品| 国产日韩欧美亚洲二区| 亚洲精品aⅴ在线观看| 欧美另类一区| 美女视频免费永久观看网站| 少妇高潮的动态图| 亚洲国产欧美在线一区| 亚洲精品456在线播放app| 成年版毛片免费区| 中国三级夫妇交换| 国产精品国产三级国产专区5o| 亚洲在久久综合| 国产精品一及| 久久人人爽人人片av| 黄色日韩在线| 一边亲一边摸免费视频| 在线免费观看不下载黄p国产| 伦精品一区二区三区| 国内揄拍国产精品人妻在线| 久久精品久久精品一区二区三区| 久久久久久久久久久丰满| 国产成人精品久久久久久| 亚洲国产欧美人成| 嘟嘟电影网在线观看| 欧美区成人在线视频| 欧美日韩视频精品一区| 国产精品三级大全| 久热久热在线精品观看| 波野结衣二区三区在线| 国产av不卡久久| 国产欧美日韩一区二区三区在线 | av专区在线播放| 91aial.com中文字幕在线观看| 欧美激情在线99| 在线观看一区二区三区激情| av网站免费在线观看视频| 黄色怎么调成土黄色| 新久久久久国产一级毛片| 欧美精品人与动牲交sv欧美| av国产精品久久久久影院| 精品99又大又爽又粗少妇毛片| 国产一区有黄有色的免费视频| 亚洲一区二区三区欧美精品 | 99久久九九国产精品国产免费| 日韩强制内射视频| 精品亚洲乱码少妇综合久久| 国产黄色免费在线视频| 一区二区三区乱码不卡18| av专区在线播放| 国产精品秋霞免费鲁丝片| 美女高潮的动态| 网址你懂的国产日韩在线| 69人妻影院| 久久久久性生活片| 精品一区二区三卡| 日日啪夜夜撸| 嫩草影院新地址| 亚洲美女搞黄在线观看| 在线天堂最新版资源| 欧美日韩视频高清一区二区三区二| 国产 一区 欧美 日韩| 国产精品成人在线| 毛片女人毛片| 免费大片黄手机在线观看| 草草在线视频免费看| 黄色视频在线播放观看不卡| av天堂中文字幕网| 卡戴珊不雅视频在线播放| 欧美高清成人免费视频www| 九色成人免费人妻av| av在线app专区| 日本午夜av视频| 欧美3d第一页| 亚洲精品日韩在线中文字幕| av国产精品久久久久影院| 联通29元200g的流量卡| 国产伦精品一区二区三区视频9| 午夜福利视频1000在线观看| 日韩成人伦理影院| 精品久久久久久电影网| 韩国高清视频一区二区三区| 亚洲人成网站在线播| 国产一区二区三区综合在线观看 | 禁无遮挡网站| 欧美成人精品欧美一级黄| 午夜福利在线观看免费完整高清在| 久久久久久国产a免费观看| 亚洲国产精品国产精品| 国产精品国产三级专区第一集| 久久人人爽av亚洲精品天堂 | 欧美高清性xxxxhd video| 亚洲欧美日韩无卡精品| 亚洲精品成人av观看孕妇| 亚洲色图av天堂| 啦啦啦中文免费视频观看日本| 国产伦精品一区二区三区四那| 中文欧美无线码| 在线免费观看不下载黄p国产| 在线观看av片永久免费下载| 男插女下体视频免费在线播放| 岛国毛片在线播放| 视频中文字幕在线观看| 日韩电影二区| 自拍欧美九色日韩亚洲蝌蚪91 | 精品人妻视频免费看| 男人和女人高潮做爰伦理| 网址你懂的国产日韩在线| 性插视频无遮挡在线免费观看| www.av在线官网国产| 免费观看无遮挡的男女| 国产精品99久久久久久久久| 91久久精品国产一区二区三区| 国产亚洲av片在线观看秒播厂| 成人国产麻豆网| 国产91av在线免费观看| 免费观看a级毛片全部| 成人二区视频| 街头女战士在线观看网站| 久久精品久久久久久久性| 麻豆久久精品国产亚洲av| 久久国内精品自在自线图片| .国产精品久久| 午夜亚洲福利在线播放| 可以在线观看毛片的网站| 久久人人爽人人爽人人片va| 成年av动漫网址| 日本午夜av视频| 晚上一个人看的免费电影| 国产真实伦视频高清在线观看| 人妻系列 视频| 在线a可以看的网站| www.色视频.com| 久久久久久久国产电影| 国产精品福利在线免费观看| 乱系列少妇在线播放| 三级男女做爰猛烈吃奶摸视频| 国产免费视频播放在线视频| 久久久欧美国产精品| 国产成人免费观看mmmm| 日日啪夜夜爽| 最新中文字幕久久久久| 免费在线观看成人毛片| 国产女主播在线喷水免费视频网站| 大香蕉97超碰在线| 黄色日韩在线| 97在线视频观看| a级一级毛片免费在线观看| 一边亲一边摸免费视频| 久热这里只有精品99| 国产精品三级大全| 国产亚洲av片在线观看秒播厂| 一级毛片 在线播放| 内地一区二区视频在线| 99热国产这里只有精品6| 亚洲精品中文字幕在线视频 | 精品久久久久久电影网| 免费看日本二区| 日本色播在线视频| 亚洲美女搞黄在线观看| 国产精品一及| 欧美潮喷喷水| 综合色av麻豆| 中文资源天堂在线| 综合色丁香网| 岛国毛片在线播放| 色视频在线一区二区三区| 天天一区二区日本电影三级| 综合色丁香网| 99re6热这里在线精品视频| 观看美女的网站| 在线播放无遮挡| 噜噜噜噜噜久久久久久91| 赤兔流量卡办理| 久久ye,这里只有精品| 18+在线观看网站| 噜噜噜噜噜久久久久久91| 看黄色毛片网站| 欧美日韩视频高清一区二区三区二| 小蜜桃在线观看免费完整版高清| 亚洲av电影在线观看一区二区三区 | 久久精品熟女亚洲av麻豆精品| 有码 亚洲区| 久久亚洲国产成人精品v|