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

    葉片表面粗糙度對高負荷低壓渦輪的流動影響

    2019-06-03 03:39張宗辰喬渭陽白濤

    張宗辰 喬渭陽 白濤

    摘 要:為了得到不同工況下表面粗糙度對渦輪葉片葉型損失的影響規(guī)律,采用數(shù)值模擬的方法對某前加載葉型在不同攻角和不同雷諾數(shù)下的流動進行了詳細的分析。結(jié)果表明,當(dāng)攻角i=0°、10°時,葉片表面并無明顯的分離現(xiàn)象出現(xiàn),當(dāng)i=20°、25°、30°時,葉片表面都出現(xiàn)了不同程度的分離,且攻角越大分離越嚴重。當(dāng)攻角一定時,增大雷諾數(shù)對抑制分離泡的出現(xiàn)有促進作用;當(dāng)雷諾數(shù)也一定時,增大葉片表面粗糙度對抑制附面層的分離有明顯的效果,且雷諾數(shù)越大抑制分離所需的粗糙度值就越低。攻角為20°,雷諾數(shù)分別等于25000、50000、100000、150000、200000時,抑制分離所需的最佳粗糙度值依次為38、14、5.1、2.5、1.7mm;攻角為25°,相同雷諾數(shù)下抑制分離所需的最佳粗糙度值依次為230、50、11、4、2.2mm;攻角為30°,雷諾數(shù)分別等于50000、100000、150000、200000時,抑制分離所需的最佳粗糙度值依次為3200、800、120、29mm。最后,建立了一套不同攻角下抑制分離的最佳粗糙度-雷諾數(shù)關(guān)系模型,并編寫了相應(yīng)的C語言程序。通過該程序,只要得知葉片工作的攻角與雷諾數(shù)大小,便可直接算出抑制附面層分離的最佳粗糙度值。

    關(guān)鍵詞:高負荷低壓渦輪;表面粗糙度;前加載;附面層分離;葉型損失

    DOI:10.15938/j.jhust.2019.02.010

    中圖分類號: V231.3

    文獻標志碼: A

    文章編號: 1007-2683(2019)02-0059-14

    收稿日期: 2017-04-18

    基金項目: 國家自然科學(xué)基金(51476134).

    作者簡介:

    喬渭陽(1963—),男,博士,教授;

    白 濤(1988—),女,碩士.

    通信作者:

    張宗辰(1994—),男,碩士研究生,Email:744493714@qq.com.

    Abstract:In order to get the effects of surface roughness on turbine profile loss in different working conditions, the flow of a frontloaded blade at different incidences and different Reynolds number was analyzed in detail by numerical simulation. The results show that when the incidence i=0° and 10°, there are no obvious separation on the blade surface. When i=20°, 25° and 30°, the boundary layer on blade surface are separated and the incidence is more serious. When the incidence is certain, the increasing of Reynolds number on the suppression of separation bubble appears to have a role in promoting; when the Reynolds number also is certain, increase the blade surface roughness of suppression of boundary layer separation have obvious effect, and Reynolds number greater suppression required to separate the roughness value is lower. When i=20° and the Reynolds number is 25000,50000,100000,150000,200000, the optimal roughness value of suppressing separation is 38mm,14mm,5.1mm,2.5mm,1.7mm, When i=25°, the optimal roughness value of suppressing separation is 230mm,50mm,11mm, 4mm,2.2mm to the same Reynolds number, When i=30° and the Reynolds number is 50000, 100000, 150000, 200000, the optimal roughness value of suppressing separation is 3200mm,800mm,120mm,29mm. Finally, a modal of optimal roughness vs. Reynolds number was built for suppressing separation in different incidences, and a C language program was written by the model. By using the program, the optimal roughness value of suppressing boundary layer separation can be calculated directly from the certain incidence and Reynolds number.

    Keywords:highlift lowpressure turbine; surface roughness; frontloaded; boundary layer separation; profile loss

    0 引 言

    20世紀50年代末以來,渦扇發(fā)動機開始應(yīng)用于民用航空工業(yè),并逐漸占據(jù)主導(dǎo)地位。隨著民用航空工業(yè)對經(jīng)濟效益日益增長的追求,提高發(fā)動機的工作效率、降低其使用成本便成為各大科研單位、發(fā)動機制造廠商著手研究的工作,而高負荷低壓渦輪由于具有許多優(yōu)良的性能,很早就成為研究的焦點[1-2]。研究表明,影響渦輪效率的關(guān)鍵因素在于葉片表面的附面層是否發(fā)生分離[3],當(dāng)發(fā)生分離時,會造成嚴重的葉型損失,使渦輪的效率急劇降低。因此,如何控制葉片表面的氣流流動,使附面層分離得到最大程度的抑制就具有重要的意義。

    目前,國內(nèi)外關(guān)于抑制附面層分離的方法主要有3類:一是通過改變?nèi)~片負荷分布來降低局部逆壓力梯度;二是主動控制,控制手段包括抽吸氣、射流等;三是被動控制,控制手段包括添加凸臺、拌線、球渦、粗糙度等[4-5]。國內(nèi)方面,西北工業(yè)大學(xué)的楊林等[6]研究了低雷諾數(shù)狀態(tài)下前、后加載葉型對高負荷低壓渦輪葉型流動損失的影響規(guī)律;中科院工程熱物理研究所的杜強等[7]對不同雷諾數(shù)、不同自由來流湍流度下高負荷低壓渦輪葉型邊界層的流動進行了詳細的數(shù)值模擬,并探討了這兩個因素對該高負荷渦輪葉型表面邊界層流動的影響;中科院工程熱物理研究所的孫爽等[8]專門研究了不同粗糙高度、不同粗糙度布置位置對高負荷低壓渦輪邊界層的影響。國外方面,Benner M W等[9-10]用數(shù)值模擬的方法研究了葉型攻角對低壓渦輪葉型損失、二次流損失的影響; 2006年,Stephen K[11]研究了被動控制中的不同沙粒高度、沙粒外形、沙粒偏斜角度、沙粒間距對轉(zhuǎn)捩的影響;同年Francesco Montomoli延續(xù)了Ramesh O N的工作[12],研究了尾跡掃掠、粗糙度的耦合控制方法對轉(zhuǎn)捩的影響。

    可以看出,到目前為止,國內(nèi)外的各種研究大部分是著重于探究某種控制手段對抑制附面層分離所產(chǎn)生的效果,或者從該控制手段的影響機理角度分析是否能夠起到抑制作用。在工程實際當(dāng)中,為了便于判斷不同工況下葉片是否發(fā)生附面層分離,以及抑制分離所需施加控制手段的具體方法和數(shù)值,應(yīng)該建立起一套某控制手段在各個工況下對抑制附面層分離所產(chǎn)生效果的模型,通過該模型可以直接找出該控制手段在不同工況下抑制附面層分離的最佳點。而目前為止,國內(nèi)外的眾多文獻當(dāng)中還沒有發(fā)現(xiàn)能夠直接反映某控制手段對抑制附面層分離所產(chǎn)生效果的模型?;诖?,本文采用數(shù)值模擬的方法,通過對某前加載葉型在不同攻角、不同雷諾數(shù)的工況下進行計算,分析其表面是否發(fā)生附面層分離。對于發(fā)生分離的工況,采用被動控制方法,通過不斷增加葉片表面的粗糙度來抑制附面層的分離,并找出能夠抑制附面層分離的最佳粗糙度值。最后,將不同攻角、不同雷諾數(shù)下的最佳粗糙度值進行整理,建立一套不同攻角、不同雷諾數(shù)下抑制附面層分離的最佳粗糙度關(guān)系模型。由于國內(nèi)外已經(jīng)做過很多的實驗,證明葉片表面粗糙度只對吸力面的附面層分離有影響,對壓力面幾乎沒有影響[],因此本文不考慮壓力面發(fā)生分離即攻角為負值的情況。

    1 數(shù)值計算方法

    ,數(shù)值模擬時,同時給出

    2.1 實驗算例驗證

    為了確保數(shù)值模擬結(jié)果的準確性,本節(jié)首先進行算例驗證。采用德國慕尼黑大學(xué)的Marco montis對本葉型的實驗測量結(jié)果進行算例校驗,該實驗測量了多組表面粗糙度下渦輪葉柵的損失系數(shù)[18]。驗證時,主要給出損失系數(shù)的數(shù)值模擬結(jié)果和實驗結(jié)果的對比。

    圖3~5分別為表面粗糙度roughness=10mm、1mm、0.1mm時渦輪葉柵的損失系數(shù)隨雷諾數(shù)變化的數(shù)值模擬結(jié)果和實驗結(jié)果的對比。從圖中可以看出,當(dāng)roughness=10mm時,損失系數(shù)的計算值與實驗值比較接近,二者的變化趨勢相同,即雷諾數(shù)越高損失系數(shù)越小,且計算值越接近實驗值。隨著粗糙度的減小,損失系數(shù)的計算值越來越靠近實驗值。當(dāng)roughness=0.1mm時,計算值與實驗值幾乎完全重合,即使雷諾數(shù)較低,計算值與實驗值也非常接近。

    由此可知,當(dāng)表面粗糙度較小時,數(shù)值模擬的結(jié)果與實驗結(jié)果非常接近。隨著粗糙度的增大,數(shù)值模擬的結(jié)果開始逐漸小于實驗結(jié)果,這可能是由于實際工作環(huán)境中條件復(fù)雜,任何微小的擾動都會造成額外的損失所致。但是,對于寬廣的雷諾數(shù)的工況來說,數(shù)值模擬的結(jié)果整體上與實驗結(jié)果保持一致,即使粗糙度最大、雷諾數(shù)最低時,計算值與實驗值的偏差也不是很大。因此,對于不同數(shù)量級的表面粗糙度值,數(shù)值模擬的結(jié)果與實驗結(jié)果基本吻合,數(shù)值模擬可以作為本文的研究方法。

    2.2 工況的選定

    2.3 0°攻角下葉片表面的流動狀態(tài)分析

    2.4 10°攻角下葉片表面的流動狀態(tài)分析

    將氣流的攻角加至10°。同樣,先驗證10°攻角下葉片的工作情況。

    圖7為i=10°,Re=25000工況下葉片表面的流線圖和速度云圖。在流線圖中,氣流在葉片前緣處較i=0°時有較大的折轉(zhuǎn),除此之外,葉片吸力面并無明顯的分離;在速度云圖中,葉片吸力面前緣存在輕微的逆速度梯度,但在逆速度梯度區(qū)以外,速度梯度又恢復(fù)了正常。這與正攻角下氣流發(fā)生了折轉(zhuǎn)有關(guān)。

    由此可知,10°攻角下氣流在前緣處的折轉(zhuǎn)更大,由此造成的前緣損失也更大,但在葉片中部及后部附面層依然保持良好,沒有發(fā)生明顯的分離現(xiàn)象,因此,不用對10°攻角下的葉片施加粗糙度。

    2.5 20°攻角下表面粗糙度對流動的影響分析

    20°攻角時,攻角較大,附面層已經(jīng)發(fā)生了明顯的分離。因此,本節(jié)分別對Re=25000、50000、100000、150000、200000的工況進行分析。

    布圖、壁面剪切力圖。從流線圖中可以看出,20°攻角時葉片的吸力面發(fā)生了明顯的附面層分離,分離的起始點位于葉片前緣上部,結(jié)束點位于葉背中部,葉片前緣處的氣流也有更明顯的折轉(zhuǎn)。在速度云圖中,葉背處存在明顯的逆速度梯度區(qū),即氣流速度由外到內(nèi)先減小后增大。在壓力云圖中,吸力面上方的壓力分布由外到內(nèi)逐漸減小,這與壓力面上的壓力分布情況形成了鮮明的對比。在總壓云圖中,葉片吸力面上方存在嚴重的總壓損失區(qū),而壓力面上的總壓分布均勻,并未發(fā)生明顯的總壓損失。

    載荷分布圖表示葉片表面的載荷分布情況,它是根據(jù)式(8)得出的葉片表面各處的靜壓力系數(shù)所連成的曲線。若葉片屬于沒有發(fā)生分離的理想情況,則載荷分布圖呈標準的加速-擴壓結(jié)構(gòu),整個曲線在變化上過渡得十分平緩。如果發(fā)生分離,則曲線在分離段會出現(xiàn)一個明顯的“平臺”,且分離越大“平臺”越高,如圖9所示。當(dāng)分離過大時,“平臺”甚至還會轉(zhuǎn)化為“波浪”,曲線變得十分不規(guī)則。在圖8所示的載荷分布圖中,曲線并無明顯的加速段,而是從一開始就已進入擴壓段結(jié)構(gòu),且擴壓段上的“平臺區(qū)”跨度較長,曲線過渡得并不均勻。由此可以判斷,吸力面上方存在著較長的分離區(qū)。

    壁面剪切力圖表示葉片表面所受氣流摩擦力的大小和方向情況。若葉片沒有發(fā)生分離,則氣流對葉片沿Z軸的摩擦力永遠為正(向右),如果發(fā)生分離,則在分離泡內(nèi)氣流對葉片沿Z軸的摩擦力為負(向左)。在圖8所示的壁面剪切力圖中,吸力面的壁面剪切力從一開始就已變?yōu)樨撝?,直到葉片的靠近尾緣處才重新恢復(fù)正值。由此可知,吸力面沿Z軸的剪切力負值區(qū)就是氣流的分離區(qū)。

    通過以上6個角度的分析可知,在i=20°,Re=25000的工況下,葉片已出現(xiàn)嚴重的附面層分離,必須對其進行抑制。現(xiàn)對葉片表面逐步增加粗糙度,同時觀察其對分離的抑制情況。通常,在開始出現(xiàn)分離的臨界Re下,附面層的厚度可以達到渦輪正常工作的最大厚度。基于此,本文將臨界Re下葉片前緣點至附面層轉(zhuǎn)捩點之間的平均附面層厚度的1.1倍定義為葉片的基礎(chǔ)粗糙度高度[8]。圖10為基礎(chǔ)粗糙度高度的布置示意圖。根據(jù)高磊、孫爽等的結(jié)果[21-22],首次設(shè)置粗糙度時,將葉片吸力面的幾何表面粗糙度設(shè)置為基礎(chǔ)粗糙度高度的0.2倍,若該粗糙度下的分離現(xiàn)象沒有明顯的變化,則將粗糙度值增加一個數(shù)量級;若該粗糙度下的分離開始出現(xiàn)變化,則在該粗糙度的數(shù)量級內(nèi)逐步增加粗糙度值,直至找出抑制分離的最佳粗糙度值。

    失,壓力梯度恢復(fù)正常。在總壓云圖中,葉片周圍的總壓分布均勻,只有吸力面上方還存在輕微的摩擦損失。在載荷分布圖中,曲線重新恢復(fù)加速-擴壓結(jié)構(gòu),且擴壓段上的“平臺區(qū)”已完全消失,曲線過渡得非常均勻。在壁面剪切力圖中,吸力面上的剪切力負值區(qū)已明顯變小,甚至快要消失。

    圖13為根據(jù)式(7)得出的損失系數(shù)關(guān)系曲線。通過曲線可以看出,損失系數(shù)隨著粗糙度的增大而減小。

    2.6 25°攻角下表面粗糙度對流動的影響分析

    逐步增加表面粗糙度。由于本工況下的分離區(qū)很大,因此估計粗糙度值比i=20°時要大很多。當(dāng)粗糙度加至roughness=230mm時,葉片表面的流線圖、速度云圖、壓力云圖、總壓云圖、載荷分布圖、壁面剪切力圖如圖15所示。

    在流線圖中,分離泡已被完全抑制,但葉片前緣處的流線存在較大的折轉(zhuǎn),這種折轉(zhuǎn)會造成一定的前緣損失。在速度云圖中,吸力面已無明顯的逆速度梯度區(qū),速度分布比較均勻。在壓力云圖中,吸力面上方的壓力分布重新恢復(fù)正常結(jié)構(gòu),壓力梯度變化規(guī)律。在總壓云圖中,葉背上方的分離損失區(qū)已基本消失,但摩擦損失比較明顯,且摩擦損失區(qū)比i=20°時更寬。在載荷分布圖中,曲線重新恢復(fù)加速-擴壓結(jié)構(gòu),且擴壓段上的“平臺區(qū)”已完全消失,曲線過渡得非常均勻。在壁面剪切力圖中,除了葉片前緣處存在一小段負值區(qū)外,其余部分均沒有負值區(qū)出現(xiàn),這說明分離現(xiàn)象已被完全抑制,而前緣處的負值區(qū)是由于氣流發(fā)生折轉(zhuǎn)所致。

    2.7 30°攻角下表面粗糙度對流動的影響分析

    30°攻角時,由于攻角太大,葉片已經(jīng)失去了正常工作的能力。數(shù)值模擬結(jié)果表明,當(dāng) 時,抑制附面層分離所需的粗糙度值已經(jīng)達到了米的量級,而實際中顯然不可能把粗糙度做到這種程度。因此,本節(jié)所研究的內(nèi)容僅作為一種理論上的預(yù)測,它可用來與小攻角下抑制分離的模型做對比。

    由于本工況下攻角太大,且雷諾數(shù)較低,經(jīng)數(shù)值模擬后發(fā)現(xiàn)無論粗糙度加到多大,都不能完全抑制附面層的分離。因此,本工況下無抑制分離的最佳粗糙度值。

    當(dāng)攻角為30°,雷諾數(shù)分別等于50000、100000、150000、200000時,抑制附面層分離的最佳粗糙度值如表4所示:

    3 不同工況下抑制附面層分離所需粗糙度的模型

    前面的數(shù)值模擬結(jié)果表明,當(dāng)攻角較小時,附面層沒有發(fā)生明顯的分離,葉片工作性能良好。當(dāng)攻角i在10°~20°之間時,附面層開始發(fā)生分離,且攻角越大分離越嚴重。對于固定的攻角,當(dāng)雷諾數(shù)較小時分離較難控制,抑制分離往往需要較大的粗糙度值;當(dāng)雷諾數(shù)較大時,附面層接近湍流附面層,分離泡對粗糙度的變化較為敏感,很小的粗糙度就能完全抑制分離泡的出現(xiàn)。由此可知,攻角一定時,抑制附面層分離所需的粗糙度值與雷諾數(shù)成反比關(guān)系,雷諾數(shù)越小粗糙度越大,雷諾數(shù)越大粗糙度越小。如果能在不同的攻角下分別得到抑制分離的最佳粗糙度值與雷諾數(shù)的函數(shù)關(guān)系,則這個關(guān)系式在工程實際當(dāng)中就有很大的意義,只要給定某一攻角下流動雷諾數(shù)的大小,就可由關(guān)系式直接得出抑制分離的最佳粗糙度值而不用再通過實驗一組一組尋找,節(jié)省了大量的時間與人力、物力,這對于工程實際來說是十分重要的。

    3.1 數(shù)學(xué)模型分析

    根據(jù)數(shù)值模擬的結(jié)果可知,固定攻角下,雷諾數(shù)與抑制分離所需的粗糙度值成反比關(guān)系,雷諾數(shù)越小粗糙度越大,雷諾數(shù)越大粗糙度越小。如果用y表示粗糙度值,用x表示雷諾數(shù),則y與x的關(guān)系如圖18所示。

    同攻角下k與a的值不同。

    由式(9)可知,若想求得y與x的具體關(guān)系,只需求出k與a的值,給定任意兩組y與x的數(shù)據(jù),就可求出k與a。因此,在某一攻角下,只要將兩組雷諾數(shù)與粗糙度的數(shù)據(jù)代入,就可得到該攻角下y與x的函數(shù)關(guān)系。由此可知,式(9)即為本節(jié)所尋找的函數(shù)關(guān)系式。

    3.2 模型的建立

    關(guān)系曲線如圖19所示。

    為了保證建立出的函數(shù)模型能盡可能地符合所有的點,這里選擇第二組與第四組數(shù)據(jù),即 。將這兩組數(shù)據(jù)代入式(9),可得

    由2.6節(jié)的分析可知, 時抑制分離的粗糙度值與雷諾數(shù)的關(guān)系如表6所示。

    關(guān)系曲線如圖20所示。

    因此,粗糙度y與雷諾數(shù)x的關(guān)系為

    關(guān)系曲線如圖21所示。

    因此,粗糙度y與雷諾數(shù)x的關(guān)系為

    式(12)即為i=30°時抑制附面層分離的函數(shù)模型。

    3.2.4 誤差分析

    等[23]指出,當(dāng)測量值與真值的相對誤差在10%以內(nèi)時,可以認為測量值符合工程實際的需要,即

    式中δ為每次測量的相對誤差,ω為測量值,μ為真值。

    表(8)~(10)分別是由式(10)~(12)計算出的最佳粗糙度值與數(shù)值模擬結(jié)果的比較。在表(8)和表(9)中,不同雷諾數(shù)下計算值與真值的相對誤差都在10%以內(nèi),且能看出高雷諾數(shù)下的相對誤差要普遍小于低雷諾數(shù)下的,說明建立出的函數(shù)模型在高雷諾數(shù)下的計算結(jié)果更接近真實結(jié)果。

    在表(10)中,當(dāng)Re=100000、150000時,由于函數(shù)模型存在一定的精度誤差,導(dǎo)致計算結(jié)果并不完全等于真實結(jié)果。當(dāng)Re=200000時,計算值與真值的相對誤差δ=7.6%<10%,這符合上述要求。而當(dāng)Re=50000時,相對誤差δ高達539.6%,這說明低雷諾數(shù)下函數(shù)模型的計算值與真實值存在較大的偏差,且攻角越大偏差越明顯,這也說明當(dāng)攻角較大且雷諾數(shù)較低時,不宜采用本模型求解最佳粗糙度值。

    1

    3.2.5 結(jié)果分析

    由圖可知,在低雷諾數(shù)下,如果攻角較小,則抑制附面層分離很容易,很小的粗糙度值就能完全抑制分離泡的產(chǎn)生;當(dāng)攻角逐漸增大時,抑制分離的粗糙度值急劇增大,甚至達到了實際中根本不會出現(xiàn)的程度。隨著雷諾數(shù)的增大,3種攻角下抑制分離所需的粗糙度值都在減小,且攻角越大減小得越劇烈,3條曲線有逐漸匯合的趨勢。當(dāng)雷諾數(shù)達到最大時,3條曲線已接近重合,抑制分離的粗糙度值非常相近。由此可以推斷,隨著雷諾數(shù)的繼續(xù)增大,3種攻角下的粗糙度值會繼續(xù)減小,曲線之間的距離也越來越小,粗糙度值會更相近。

    圖23是根據(jù)式(10)~(12)編成的C語言程序,通過此程序可直接求出不同工況下抑制分離所需的粗糙度值。求解時,需根據(jù)提示輸入氣流的攻角i與雷諾數(shù)Re,當(dāng)攻角i=0°或10°時,程序會自動提示“不需增加粗糙度”,當(dāng) 時,程序會算出輸入工況下的最佳粗糙度。

    4 結(jié) 論

    本文采用數(shù)值模擬的方法分析了某前加載葉型在攻角 和雷諾數(shù) 。對于出現(xiàn)分離的工況,采用增加粗糙度的方法抑制其分離,得出了這些工況下抑制附面層分離的最佳粗糙度值,建立了一套在攻角 幾種工況下抑制附面層分離的粗糙度-雷諾數(shù)關(guān)系模型,通過給定攻角與雷諾數(shù)的大小便可得出抑制分離的最佳粗糙度值。本文得出的主要結(jié)論如下:

    1)附面層分離現(xiàn)象會造成嚴重的分離損失,且攻角越大分離越嚴重。對于本文所選用的葉型,當(dāng)攻角在10°~20°之間時,附面層開始發(fā)生分離。攻角達到30°時,分離泡會占據(jù)葉柵通道的大部分位置,葉片已不能正常工作。

    2)附面層分離現(xiàn)象與雷諾數(shù)有關(guān),雷諾數(shù)越小附面層越容易分離,雷諾數(shù)越大附面層越不易分離。在低雷諾數(shù)下,抑制分離往往需要較大的粗糙度值,且當(dāng)攻角過大時,粗糙度甚至大到實際中無法存在的程度。在高雷諾數(shù)下,分離泡對粗糙度的變化十分敏感,粗糙度的微小改變會使分離泡發(fā)生很大變化,這對抑制附面層的分離有促進作用。

    3)攻角變化引起的損失不僅只有分離損失。當(dāng)大攻角下的分離被抑制時,依然存在較大的前緣損失。

    4)固定攻角下,最佳粗糙度值與雷諾數(shù)呈負指數(shù)冪函數(shù)關(guān)系,當(dāng)雷諾數(shù)無限增大時,粗糙度趨近于0,當(dāng)雷諾數(shù)趨近于0時,粗糙度趨近于無窮大。20°攻角下,粗糙度與雷諾數(shù)的函數(shù)關(guān)系為 攻角下,粗糙度與雷諾數(shù)的函數(shù)關(guān)系為y=3.177×1012x2.299;30°攻角下,粗糙度與雷諾數(shù)的函數(shù)關(guān)系為 攻角越小且雷諾數(shù)越高,函數(shù)模型的計算值越接近真實值;攻角越大且雷諾數(shù)越低,函數(shù)模型的計算值與真實值偏差越大。

    參 考 文 獻:

    [1] WILSER D C. The Technical and Economic Relevance of Understanding Boundary Layer Transition in Gas Turbine Engines[R].NASA/CP-1998-206958,1998.

    [2] HOWELL R J. Wake Separation Bubble Interactions in Low Reynolds Number Turbomachinery [D].Cambridge:Cambridge University,1999.

    [3] DENTON J D. Loss Mechanisms in Turbomachines[J]. Journal of Turbomachinery, 1993, 115(4):621.

    [4] SCHULTE V, HODSON H P. Unsteady WakeInduced Boundary Layer Transition in High Lift LP Turbines[J]. Journal of Turbomachinery, 1998, 120(1):28.

    [5] ZHANG X F, VERA M, HODSON H P, et al. Separation and Transition Control on An Aftloaded Ultrahighlift LP Turbine Blade at Low Reynolds Numbers:LowSpeed Investigation[R]. ASME paper 2005-GT-68892,2005.

    [6] 楊林,喬渭陽,羅華玲,等.低雷諾數(shù)高負荷低壓渦輪葉型的氣動設(shè)計[J].航空動力學(xué)報,2013,28(5):1019.

    [7] 杜強,朱俊強,溫殿中.高負荷低壓渦輪邊界層轉(zhuǎn)捩預(yù)測及其機理分析[J].工程熱物理學(xué)報,2010,31(5):761.

    [8] 孫爽,雷志軍,朱俊強,等.粗糙度對超高負荷低壓渦輪邊界層影響[J].推進技術(shù),2014,35(3):347.

    [9] BENNER M W, SJOLANDER S A, MOUSTAPHA S H. Measurements of Secondary Flows Downstream of a Turbine Cascade at Offdesign Incidence[R].ASME GT 2004-53786,2004.

    [10]BENNER M W, SJOLANDER S A, MOUSTAPHA S H. Measurements of Secondary Flows in Turbine Cascade at Offdesign Incidence[R].ASME Turbo Expo, GT 97-382,1997.

    [11]ROBERTS S K, YARAS M I. Effects of SurfaceRoughness Geometry on SeparationBubble Transition[J]. Journal of Turbomachinery,2006, 128(2):349.

    [12]FRANCESCO Montomoli, HOWARD Hodson. Effect of Roughness and Unsteadiness on the Performance of a New Low Pressure Turbine Blade at Low Reynolds Numbers[J]. Journal of Turbomachinery, 2010, 132(3):031018.

    [13]白濤.葉片幾何偏差對渦輪性能的影響[D].北京:北京航空航天大學(xué),2013.

    [14]張儀,王曉東,胡昊,等.湍流模型對湍流射流CFD模擬的影響[J].推進技術(shù),2016,37(6):1049.

    [15]張波,李偉,杜強,等.U型槽對高負荷低壓渦輪葉型攻角特性影響[J].航空動力學(xué)報,2012,27(7):1503.

    [16]SCHILICHTING H. Experiment Investigation of the Problem of Surface Roughness[R].NACA Technical Memorandum,1937,823(4),1937.

    [17]NIKURADSE J. Law of Flow in Rough Pipe[M].Technical Memordum 1292 NACA,1950.

    [18]MARCO Montis, ANDREAS Fiala. Effect of Surface Roughness on Loss Behavior, Aerodynamic Loading and Boundary Layer Development of a Lowpressure Gas Turbine Airfoil[R].ASME GT2010-23317,2010.

    [19]孫爽,雷志軍,盧新根,等.來流條件對超高負荷低壓渦輪附面層非定常特性影響的實驗研究[J].推進技術(shù),2016,37(4):653.

    [20]HODSON H P, HOWELL R J. The Role of Transition in Highlift Lowpressure Turbines for Aero Engine[J]. Progress in Aerospace Science, 2005, 41:419.

    [21]高磊,王子楠,耿少娟,等.粗糙度對壓氣機葉柵損失特性影響的實驗研究[J].推進技術(shù),2016,37(7):1263.

    [22]孫爽,雷志軍,盧新根,等.基于表面粗糙度的超高負荷低壓渦輪葉片附面層控制[J].航空動力學(xué)報,2016, 31(4):836.

    [23]GEERS T L. Objective Error Measure for the Comparison of Calculated and Measured Transient Response Histories[J]. Shock and Vibration Bulletin, 1984,54:99.

    (編輯:溫澤宇)

    亚洲国产欧美人成| 午夜福利免费观看在线| 欧美xxxx黑人xx丫x性爽| 国产成人欧美在线观看| 中文资源天堂在线| 午夜福利在线在线| 国产精品爽爽va在线观看网站| 亚洲成av人片免费观看| 欧美极品一区二区三区四区| 黄色一级大片看看| 五月伊人婷婷丁香| 狠狠狠狠99中文字幕| 97超级碰碰碰精品色视频在线观看| 男女那种视频在线观看| 色哟哟哟哟哟哟| 免费av毛片视频| 久久久久久久精品吃奶| 亚洲人成网站在线播| 熟妇人妻久久中文字幕3abv| 午夜影院日韩av| 欧美一区二区国产精品久久精品| 成年人黄色毛片网站| 99久久成人亚洲精品观看| 久久久久国产精品人妻aⅴ院| 免费看a级黄色片| 国产亚洲欧美在线一区二区| aaaaa片日本免费| 天堂av国产一区二区熟女人妻| 亚洲人成网站在线播放欧美日韩| 高清日韩中文字幕在线| 免费观看的影片在线观看| 国内久久婷婷六月综合欲色啪| 成人高潮视频无遮挡免费网站| 村上凉子中文字幕在线| 欧美激情在线99| 国产精品伦人一区二区| 久久久久国产精品人妻aⅴ院| 久久国产乱子免费精品| а√天堂www在线а√下载| 桃色一区二区三区在线观看| 99国产精品一区二区三区| 人人妻人人澡欧美一区二区| 亚洲av日韩精品久久久久久密| 欧美乱色亚洲激情| 很黄的视频免费| 91久久精品国产一区二区成人| 超碰av人人做人人爽久久| 深夜a级毛片| 丁香欧美五月| 欧美又色又爽又黄视频| 成年版毛片免费区| 成年女人毛片免费观看观看9| 黄色配什么色好看| av在线老鸭窝| 脱女人内裤的视频| 窝窝影院91人妻| 成人永久免费在线观看视频| 亚洲av五月六月丁香网| 小说图片视频综合网站| 国产精品一区二区性色av| 亚洲国产高清在线一区二区三| 国产精品久久久久久精品电影| 免费人成视频x8x8入口观看| 观看免费一级毛片| 国产爱豆传媒在线观看| 欧美日韩中文字幕国产精品一区二区三区| 很黄的视频免费| 韩国av一区二区三区四区| 人人妻,人人澡人人爽秒播| 综合色av麻豆| 久9热在线精品视频| 亚洲欧美精品综合久久99| 色5月婷婷丁香| 亚洲专区中文字幕在线| 一个人免费在线观看的高清视频| 国产中年淑女户外野战色| 丁香六月欧美| 自拍偷自拍亚洲精品老妇| 欧美性猛交╳xxx乱大交人| 人妻夜夜爽99麻豆av| 久久婷婷人人爽人人干人人爱| 久久亚洲真实| 亚洲在线观看片| 国产一级毛片七仙女欲春2| 国产私拍福利视频在线观看| 一进一出抽搐动态| 国产熟女xx| 2021天堂中文幕一二区在线观| av在线天堂中文字幕| 精品无人区乱码1区二区| 久久久久免费精品人妻一区二区| 国内揄拍国产精品人妻在线| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久成人av| 国产极品精品免费视频能看的| 亚洲自偷自拍三级| 无人区码免费观看不卡| 国内精品一区二区在线观看| 久久久久久久久大av| 一边摸一边抽搐一进一小说| 中文字幕av在线有码专区| 99国产极品粉嫩在线观看| 亚洲性夜色夜夜综合| 美女高潮的动态| 国产精品亚洲美女久久久| 日本免费一区二区三区高清不卡| 国产精品精品国产色婷婷| 成年人黄色毛片网站| 窝窝影院91人妻| 婷婷精品国产亚洲av在线| 超碰av人人做人人爽久久| 午夜福利成人在线免费观看| 精品一区二区三区视频在线观看免费| 深爱激情五月婷婷| 一本精品99久久精品77| 国产乱人视频| 亚洲七黄色美女视频| 免费黄网站久久成人精品 | 日本 av在线| 欧美成人一区二区免费高清观看| 丁香六月欧美| 一区福利在线观看| 国产亚洲av嫩草精品影院| 国产精品电影一区二区三区| 午夜免费成人在线视频| 欧美在线一区亚洲| 怎么达到女性高潮| 美女免费视频网站| 亚洲精品久久国产高清桃花| 小蜜桃在线观看免费完整版高清| 精品久久久久久,| 丝袜美腿在线中文| 国产色婷婷99| 日韩成人在线观看一区二区三区| 最后的刺客免费高清国语| 看免费av毛片| 一级a爱片免费观看的视频| 亚洲黑人精品在线| 99热精品在线国产| 久久人人精品亚洲av| 久久6这里有精品| 91麻豆av在线| 国产一区二区三区视频了| a级毛片a级免费在线| 99久国产av精品| 欧美最新免费一区二区三区 | 中出人妻视频一区二区| 午夜两性在线视频| 成年女人永久免费观看视频| 人妻久久中文字幕网| 女同久久另类99精品国产91| 国产精品爽爽va在线观看网站| 如何舔出高潮| 1000部很黄的大片| 狠狠狠狠99中文字幕| 中文在线观看免费www的网站| 国产亚洲精品久久久久久毛片| 十八禁网站免费在线| 久久午夜亚洲精品久久| 国产毛片a区久久久久| 国产91精品成人一区二区三区| 99精品久久久久人妻精品| 在线观看舔阴道视频| eeuss影院久久| 欧美成人性av电影在线观看| 国产色婷婷99| 少妇丰满av| 毛片女人毛片| 亚洲欧美日韩东京热| 美女免费视频网站| 亚洲五月婷婷丁香| 亚洲一区二区三区不卡视频| 精品乱码久久久久久99久播| 亚洲成a人片在线一区二区| 伦理电影大哥的女人| 深夜精品福利| 欧美日韩亚洲国产一区二区在线观看| 国产高清视频在线播放一区| 欧美成人a在线观看| 男女那种视频在线观看| 亚洲国产精品999在线| 午夜老司机福利剧场| 可以在线观看毛片的网站| 亚洲精品日韩av片在线观看| 桃红色精品国产亚洲av| 内射极品少妇av片p| 亚洲av.av天堂| 熟女人妻精品中文字幕| 国产精品国产高清国产av| 亚州av有码| 日本撒尿小便嘘嘘汇集6| 成人国产综合亚洲| 亚洲无线在线观看| 亚洲国产精品久久男人天堂| 久久人人精品亚洲av| 高清毛片免费观看视频网站| 熟妇人妻久久中文字幕3abv| АⅤ资源中文在线天堂| 亚洲国产精品久久男人天堂| 亚洲自拍偷在线| 亚洲五月婷婷丁香| 伦理电影大哥的女人| 精品久久久久久久久av| 丰满人妻熟妇乱又伦精品不卡| 91午夜精品亚洲一区二区三区 | 亚洲男人的天堂狠狠| 三级国产精品欧美在线观看| 色播亚洲综合网| 99久久精品一区二区三区| 18禁黄网站禁片免费观看直播| 嫁个100分男人电影在线观看| 中文资源天堂在线| 少妇的逼好多水| 欧美成人免费av一区二区三区| 国产免费一级a男人的天堂| 狠狠狠狠99中文字幕| 啦啦啦韩国在线观看视频| av国产免费在线观看| 午夜精品在线福利| av福利片在线观看| 一个人免费在线观看电影| 国产精品一区二区免费欧美| 欧美bdsm另类| 亚洲人与动物交配视频| 亚洲国产高清在线一区二区三| 在线观看美女被高潮喷水网站 | 久久精品夜夜夜夜夜久久蜜豆| 久久久久免费精品人妻一区二区| 国产蜜桃级精品一区二区三区| 可以在线观看的亚洲视频| 又黄又爽又免费观看的视频| 国产一区二区三区在线臀色熟女| 午夜精品一区二区三区免费看| 欧美在线黄色| xxxwww97欧美| 国产精品久久久久久久电影| 69av精品久久久久久| 欧美最黄视频在线播放免费| 久久精品国产亚洲av香蕉五月| 无人区码免费观看不卡| 动漫黄色视频在线观看| 欧美zozozo另类| 国产欧美日韩精品一区二区| 久久这里只有精品中国| 欧美色欧美亚洲另类二区| 亚洲 欧美 日韩 在线 免费| 久久亚洲真实| 免费看光身美女| 亚洲精品粉嫩美女一区| 国产aⅴ精品一区二区三区波| 精品人妻视频免费看| 午夜精品久久久久久毛片777| 国产乱人伦免费视频| 成人特级av手机在线观看| 欧美国产日韩亚洲一区| 搡老岳熟女国产| 国产人妻一区二区三区在| 成年版毛片免费区| 狠狠狠狠99中文字幕| 国产在线男女| 亚洲综合色惰| 麻豆av噜噜一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 免费看美女性在线毛片视频| 精品久久久久久久人妻蜜臀av| 久久精品国产亚洲av香蕉五月| 三级毛片av免费| 国产亚洲精品综合一区在线观看| 小说图片视频综合网站| 中文亚洲av片在线观看爽| 天堂√8在线中文| 一边摸一边抽搐一进一小说| 在线观看舔阴道视频| www.色视频.com| 一本一本综合久久| 欧美绝顶高潮抽搐喷水| 精品乱码久久久久久99久播| 亚洲人成电影免费在线| 欧美日韩国产亚洲二区| 久久这里只有精品中国| 一本一本综合久久| 成人一区二区视频在线观看| www日本黄色视频网| 国产精品嫩草影院av在线观看 | 国产中年淑女户外野战色| 看十八女毛片水多多多| www.999成人在线观看| 成人鲁丝片一二三区免费| 成人国产综合亚洲| 有码 亚洲区| 成熟少妇高潮喷水视频| 99久久99久久久精品蜜桃| www.色视频.com| 国产成人av教育| 亚洲专区国产一区二区| 波多野结衣高清作品| 免费看日本二区| 在线天堂最新版资源| a级毛片免费高清观看在线播放| 如何舔出高潮| 成人毛片a级毛片在线播放| 国产精品久久久久久人妻精品电影| 亚洲在线观看片| 国产野战对白在线观看| 两个人视频免费观看高清| 国内久久婷婷六月综合欲色啪| 精品福利观看| 亚洲片人在线观看| 18禁裸乳无遮挡免费网站照片| 天堂网av新在线| 日韩人妻高清精品专区| 90打野战视频偷拍视频| 久久99热6这里只有精品| 成人永久免费在线观看视频| 日韩欧美精品v在线| 97热精品久久久久久| 亚洲人成网站高清观看| 好男人电影高清在线观看| 亚洲国产欧美人成| 国产美女午夜福利| 九色国产91popny在线| 我的老师免费观看完整版| 久久久久久国产a免费观看| 国内久久婷婷六月综合欲色啪| 少妇被粗大猛烈的视频| 欧美丝袜亚洲另类 | 国产欧美日韩一区二区三| 亚洲人成网站在线播| 亚洲精品日韩av片在线观看| 桃色一区二区三区在线观看| 村上凉子中文字幕在线| 欧美3d第一页| 日韩欧美精品v在线| 99久久精品国产亚洲精品| 欧美一区二区精品小视频在线| 久久久久九九精品影院| 国产精品亚洲一级av第二区| www.www免费av| 精品福利观看| 欧美成人一区二区免费高清观看| av福利片在线观看| 俺也久久电影网| 又紧又爽又黄一区二区| 久久久久九九精品影院| 亚洲中文字幕一区二区三区有码在线看| 性色avwww在线观看| 婷婷亚洲欧美| 天天躁日日操中文字幕| .国产精品久久| 亚洲国产精品成人综合色| 欧美日韩瑟瑟在线播放| 偷拍熟女少妇极品色| 免费av观看视频| a级毛片a级免费在线| 赤兔流量卡办理| 亚洲国产欧洲综合997久久,| 精品一区二区三区av网在线观看| 国产av不卡久久| 久久国产精品人妻蜜桃| 十八禁网站免费在线| 天堂网av新在线| 五月玫瑰六月丁香| 亚洲三级黄色毛片| 精品久久国产蜜桃| 欧美高清成人免费视频www| 成年女人永久免费观看视频| 欧美高清成人免费视频www| or卡值多少钱| 国产免费男女视频| 狠狠狠狠99中文字幕| 亚洲最大成人中文| 夜夜爽天天搞| 又紧又爽又黄一区二区| 午夜福利在线在线| 国产精品不卡视频一区二区 | 国产精品三级大全| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩高清在线视频| av中文乱码字幕在线| 亚洲aⅴ乱码一区二区在线播放| 成年免费大片在线观看| 国产69精品久久久久777片| 亚洲最大成人中文| 热99在线观看视频| 男人和女人高潮做爰伦理| 美女被艹到高潮喷水动态| 黄色女人牲交| 欧美日韩乱码在线| 国产成人aa在线观看| 精品一区二区三区av网在线观看| 亚洲 国产 在线| 能在线免费观看的黄片| 少妇的逼水好多| 18禁黄网站禁片免费观看直播| 中亚洲国语对白在线视频| 亚洲男人的天堂狠狠| 成年女人看的毛片在线观看| 亚洲av电影在线进入| 亚洲美女搞黄在线观看 | 久久久久性生活片| 免费在线观看影片大全网站| 成人av在线播放网站| 中文字幕av成人在线电影| 18禁黄网站禁片午夜丰满| 91字幕亚洲| 97超级碰碰碰精品色视频在线观看| 村上凉子中文字幕在线| 久久草成人影院| 亚洲精品成人久久久久久| 中文字幕免费在线视频6| 久久精品国产亚洲av香蕉五月| 中文在线观看免费www的网站| 深夜精品福利| 中文字幕av在线有码专区| 国产精品精品国产色婷婷| 搞女人的毛片| 欧美成人a在线观看| 中文亚洲av片在线观看爽| 精品欧美国产一区二区三| 欧美一区二区精品小视频在线| 中文字幕精品亚洲无线码一区| 2021天堂中文幕一二区在线观| 免费在线观看成人毛片| av国产免费在线观看| 久久伊人香网站| 国产精品一及| 免费人成在线观看视频色| 亚洲av日韩精品久久久久久密| 色5月婷婷丁香| 女人十人毛片免费观看3o分钟| 欧美最新免费一区二区三区 | 免费无遮挡裸体视频| 精品久久久久久久末码| 亚洲第一区二区三区不卡| 欧美日韩福利视频一区二区| 亚洲精品一区av在线观看| 婷婷丁香在线五月| 黄片小视频在线播放| 国内精品久久久久精免费| 成人特级av手机在线观看| 久久久久久久久久成人| 18美女黄网站色大片免费观看| 亚洲国产高清在线一区二区三| 国产人妻一区二区三区在| 日本五十路高清| 色综合欧美亚洲国产小说| 老司机福利观看| 亚洲av成人不卡在线观看播放网| 国产精品久久久久久久久免 | 如何舔出高潮| 亚洲激情在线av| 国产探花在线观看一区二区| 欧美一区二区亚洲| 亚洲av二区三区四区| 欧美成狂野欧美在线观看| 成人国产一区最新在线观看| 精品日产1卡2卡| 久久久久久久久大av| 免费人成在线观看视频色| av在线天堂中文字幕| 国产精品一区二区三区四区免费观看 | 国产av麻豆久久久久久久| 久久久精品欧美日韩精品| 国产毛片a区久久久久| 嫁个100分男人电影在线观看| 国内少妇人妻偷人精品xxx网站| 美女高潮的动态| 国产伦精品一区二区三区视频9| 精品一区二区三区av网在线观看| 国产精品自产拍在线观看55亚洲| 网址你懂的国产日韩在线| 男插女下体视频免费在线播放| 好男人电影高清在线观看| 国产精品伦人一区二区| 免费电影在线观看免费观看| 色综合婷婷激情| 亚洲欧美日韩卡通动漫| 国产一级毛片七仙女欲春2| 亚洲人成伊人成综合网2020| 搡女人真爽免费视频火全软件 | 禁无遮挡网站| 午夜免费成人在线视频| 直男gayav资源| www.999成人在线观看| 亚洲国产日韩欧美精品在线观看| 中亚洲国语对白在线视频| 一区二区三区高清视频在线| 亚洲欧美精品综合久久99| 丰满的人妻完整版| 真人做人爱边吃奶动态| 91久久精品国产一区二区成人| 亚洲第一区二区三区不卡| 在线看三级毛片| 婷婷六月久久综合丁香| 天堂av国产一区二区熟女人妻| 91在线精品国自产拍蜜月| 国产视频一区二区在线看| 精品人妻熟女av久视频| 久久精品久久久久久噜噜老黄 | 人妻夜夜爽99麻豆av| 久久久久性生活片| 丰满的人妻完整版| 午夜日韩欧美国产| 亚洲精华国产精华精| ponron亚洲| 美女高潮的动态| 高清日韩中文字幕在线| 午夜日韩欧美国产| 18美女黄网站色大片免费观看| 日韩人妻高清精品专区| 国产成人啪精品午夜网站| 亚洲内射少妇av| 欧美精品啪啪一区二区三区| 最近在线观看免费完整版| 日韩欧美精品v在线| 国产aⅴ精品一区二区三区波| or卡值多少钱| 在线a可以看的网站| 狂野欧美白嫩少妇大欣赏| 亚洲片人在线观看| 真人做人爱边吃奶动态| 日本 欧美在线| 91字幕亚洲| 国产成人啪精品午夜网站| 国产精品三级大全| 国产黄片美女视频| 少妇人妻精品综合一区二区 | 欧美黄色淫秽网站| 日韩欧美三级三区| 久久久久久久久久成人| 欧美乱妇无乱码| 午夜影院日韩av| 十八禁国产超污无遮挡网站| 国内毛片毛片毛片毛片毛片| 夜夜看夜夜爽夜夜摸| 美女高潮的动态| 亚州av有码| 欧美绝顶高潮抽搐喷水| 1024手机看黄色片| 少妇被粗大猛烈的视频| 精品国内亚洲2022精品成人| 日韩av在线大香蕉| 在线观看av片永久免费下载| 内地一区二区视频在线| 成年女人永久免费观看视频| 一个人看的www免费观看视频| 欧美xxxx性猛交bbbb| 久久99热这里只有精品18| 国产激情偷乱视频一区二区| 夜夜爽天天搞| 国产精品影院久久| 久久久久精品国产欧美久久久| 国产一区二区三区在线臀色熟女| 久久香蕉精品热| 99热这里只有是精品在线观看 | 在线观看舔阴道视频| 男插女下体视频免费在线播放| 国产一区二区在线av高清观看| 伊人久久精品亚洲午夜| 午夜福利在线观看吧| 欧美一区二区精品小视频在线| 成熟少妇高潮喷水视频| 国产一区二区亚洲精品在线观看| 伦理电影大哥的女人| 热99re8久久精品国产| 久久精品久久久久久噜噜老黄 | 精品人妻1区二区| av黄色大香蕉| 俄罗斯特黄特色一大片| 亚洲一区二区三区色噜噜| 成人鲁丝片一二三区免费| www.999成人在线观看| 亚洲最大成人中文| 国产美女午夜福利| 在线观看一区二区三区| 久久久久久久久中文| 日本 欧美在线| 一区二区三区免费毛片| 国产av在哪里看| 国产在视频线在精品| 噜噜噜噜噜久久久久久91| 欧美最新免费一区二区三区 | 1000部很黄的大片| 亚洲avbb在线观看| 嫩草影院新地址| 亚洲人成网站在线播放欧美日韩| 最新在线观看一区二区三区| 亚洲精品色激情综合| 99国产精品一区二区三区| 1024手机看黄色片| 国产野战对白在线观看| 国产色爽女视频免费观看| 搡老妇女老女人老熟妇| 精品人妻偷拍中文字幕| 高清在线国产一区| 国产伦人伦偷精品视频| 每晚都被弄得嗷嗷叫到高潮| 色av中文字幕| 国产精品三级大全| 久久国产精品影院| 欧美黄色淫秽网站| 亚洲av成人av| 三级男女做爰猛烈吃奶摸视频| x7x7x7水蜜桃| 欧美黄色片欧美黄色片| 欧美不卡视频在线免费观看| 少妇人妻精品综合一区二区 | 精品日产1卡2卡| 久久久久久国产a免费观看| 中文亚洲av片在线观看爽| 久久草成人影院| 亚洲国产欧美人成| 成人特级黄色片久久久久久久| 国产v大片淫在线免费观看| 亚洲熟妇熟女久久|