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

    基于后緣小翼的旋翼翼型動態(tài)失速控制分析

    2017-11-20 01:45:01馬奕揚招啟軍趙國慶
    航空學(xué)報 2017年3期
    關(guān)鍵詞:小翼后緣升力

    馬奕揚, 招啟軍,*, 趙國慶

    1.南京航空航天大學(xué) 直升機旋翼動力學(xué)國家級重點實驗室, 南京 210016 2.中航工業(yè)第一飛機設(shè)計研究院 總體氣動所, 西安 710089

    基于后緣小翼的旋翼翼型動態(tài)失速控制分析

    馬奕揚1, 招啟軍1,*, 趙國慶2

    1.南京航空航天大學(xué) 直升機旋翼動力學(xué)國家級重點實驗室, 南京 210016 2.中航工業(yè)第一飛機設(shè)計研究院 總體氣動所, 西安 710089

    針對后緣小翼(TEF)的典型運動參數(shù)對旋翼翼型動態(tài)失速特性的控制進行了研究。發(fā)展了一套適用于帶有后緣小翼控制的旋翼翼型非定常流動特性模擬的高效、高精度CFD方法。通過求解Poisson方程生成圍繞旋翼翼型的黏性貼體和正交網(wǎng)格,為保證后緣小翼附近的網(wǎng)格生成質(zhì)量,建立了基于翼型點重構(gòu)的方法來描述后緣小翼的偏轉(zhuǎn)運動;為克服變形網(wǎng)格方法可能導(dǎo)致網(wǎng)格畸變的不足,發(fā)展了一套適用于帶有后緣小翼控制的旋翼翼型運動嵌套網(wǎng)格方法。基于非定常雷諾平均Navier-Stokes(URANS)方程、雙時間法、Spalart-Allmaras(S -A)湍流模型和Roe-Monotone Upwind-centered Scheme for Conservation Laws(Roe-MUSCL)插值格式,發(fā)展了旋翼翼型非定常氣動特性分析的高精度數(shù)值方法,并采用Lower-Upper Symmetric Gauss-Seidel(LU-SGS)隱式時間推進方法及并行技術(shù)提高計算效率。以有試驗結(jié)果驗證的HH-02翼型和SC1095翼型為算例,精確捕捉了動態(tài)失速狀態(tài)下的氣動力遲滯效應(yīng),驗證了本文方法的有效性。著重針對SC1095旋翼翼型的動態(tài)失速狀態(tài)開展后緣小翼的控制分析,提出了可以體現(xiàn)翼型升力、阻力及力矩綜合特性的關(guān)系式Po和Pc,揭示了后緣小翼振蕩頻率、相位差和偏轉(zhuǎn)幅值對動態(tài)失速特性影響的規(guī)律。研究結(jié)果表明:當后緣小翼偏轉(zhuǎn)的相對運動頻率為1.0,且小翼運動規(guī)律與翼型振蕩規(guī)律之間的相位差為0°時,后緣小翼能夠更好地抑制翼型動態(tài)失速現(xiàn)象;在此狀態(tài)下,當偏轉(zhuǎn)幅值為10°時,SC1095翼型最大阻力系數(shù)和最大力矩系數(shù)可以分別降低19%和27%。

    旋翼; 翼型; 動態(tài)失速; 后緣小翼; 參數(shù)分析; 非定常雷諾平均Navier-Stokes方程

    旋翼工作在嚴重非對稱、非定常的渦流場中,容易產(chǎn)生氣流分離進而在槳葉后行側(cè)出現(xiàn)復(fù)雜的非定常動態(tài)失速現(xiàn)象[1-2]。該現(xiàn)象會造成槳葉局部升力的突降,同時阻力和力矩會有一個明顯的峰值,翼型氣動力中心不再穩(wěn)定,會引起旋翼顫振和額外的振動載荷,進而制約了直升機氣動性能和飛行速度的提高[3]。因此,研究旋翼(翼型)動態(tài)失速的發(fā)生機理,抑制動態(tài)失速引起的阻力和力矩的發(fā)散具有重要意義。

    通過主動流動控制技術(shù)來抑制旋翼動態(tài)失速是當前研究的一個新方向[4]。由于翼型前緣形成脫體渦和產(chǎn)生的氣流分離現(xiàn)象是動態(tài)失速的原因,動態(tài)失速主動控制的研究大都集中于翼型前緣,包括翼型固定前緣條[5]、前緣布置等離子作動器[6]、翼型前緣下垂裝置[7-8]以及可旋轉(zhuǎn)前緣[9]等。但是由于槳葉屬細長梁結(jié)構(gòu)、旋轉(zhuǎn)氣動面處在非定常氣動環(huán)境中,翼型前緣的變化往往會引起較大的重心位置變化和載荷變化,因此這些方法在直升機槳葉上實現(xiàn)都比較困難[10]。相比于前緣,槳葉后緣部分的結(jié)構(gòu)較為簡單,有較大的操作空間。因此,旋翼(翼型)后緣附加小翼(Trailing-Edge Flap,TEF)是一種很有潛力的主動流動控制形式。

    2004年,F(xiàn)eszty等[11]通過離散渦方法,對NACA0012翼型進行了后緣小翼對俯仰力矩影響研究,結(jié)果初步表明當翼型的迎角較大時,通過脈沖形式的小翼控制可以降低動態(tài)失速過程中附加的低頭力矩峰值。2006年,Krzysiak和Narkiewicz[12]對帶有后緣小翼的NACA0012翼型進行了低減縮頻率狀態(tài)下的風(fēng)洞試驗,研究后緣小翼的不同運動頻率和相位差對翼型非定常氣動載荷的影響。2007年,Gerontakos和Lee[13]對帶有脈沖形式后緣小翼運動的NACA0015翼型在動態(tài)失速情況下的氣動力進行了測量,后緣小翼的運動規(guī)律與Feszty的研究結(jié)果相似,發(fā)現(xiàn)脈沖形式的后緣小翼相對于翼型振蕩規(guī)律作用越晚,控制力矩系數(shù)效果越好。2011年,Lee和Su[14]進一步發(fā)展Gerontakos的工作,對低速狀態(tài)下帶后緣小翼的NACA0015翼型進行了風(fēng)洞試驗,后緣小翼按正弦規(guī)律變化。試驗結(jié)果進一步說明了后緣小翼對動態(tài)載荷遲滯回線有很大的影響,且情況較為復(fù)雜,尚未有定性的結(jié)論;但是通過試驗研究,他們發(fā)現(xiàn)正弦運動形式的后緣小翼與脈沖形式的后緣小翼均可以有效地抑制阻力和力矩系數(shù)發(fā)散,并且與后者相比,控制參數(shù)大大減少,因此更具發(fā)展?jié)摿Α?/p>

    國內(nèi)通過后緣小翼控制旋翼動態(tài)失速的研究較少。王進等[15]利用Fluent計算平臺對在小角度變化條件下,后緣襟翼(小翼)參數(shù)以及馬赫數(shù)等因素對翼型氣動特性進行了仿真,獲得了一些計算結(jié)果,但尚未對動態(tài)失速控制進行研究。劉洋和向錦武[16]研究了脈沖形式的后緣襟翼(小翼)激勵幅值和時長等對升力和俯仰力矩系數(shù)的影響。王榮和夏品奇[17]采用經(jīng)驗?zāi)P瓦M行了后緣小翼用于控制槳葉動態(tài)失速的研究,表明后緣小翼的合理偏轉(zhuǎn)可延遲動態(tài)失速的發(fā)生,但由于模型精度和適用范圍的限制,尚未開展系統(tǒng)性的參數(shù)影響分析。

    盡管先期的試驗和數(shù)值研究已表明后緣小翼在控制動態(tài)失速方面有很大的應(yīng)用潛力,但目前關(guān)于按正弦規(guī)律運動的后緣小翼控制旋翼翼型動態(tài)失速特性的系統(tǒng)性參數(shù)分析依然缺乏。此外,面對動態(tài)失速環(huán)境下的氣流分離、再附等非定常黏性流動現(xiàn)象,需要采用高精度的CFD方法來進行有效模擬,這在目前仍有很大的難度。一方面,翼型及小翼空間位置的不斷變化對翼型貼體網(wǎng)格的生成提出了很高的要求,單純采用變形網(wǎng)格方法很難滿足網(wǎng)格正交性和貼體性的控制要求,且當運動幅度較大時,變形網(wǎng)格方法可能導(dǎo)致網(wǎng)格畸變;另一方面,由于旋翼復(fù)雜的運動特征和氣動環(huán)境,后緣小翼控制下的旋翼翼型非定常渦流動特性的捕捉對CFD方法提出了很高的精度要求。

    針對這些問題,本文采用發(fā)展的Top-Map法生成圍繞帶有后緣小翼的旋翼翼型運動嵌套網(wǎng)格;為了提高求解的精度,將Roe-MUSCL(Roe-Monotone Upwind-centered Scheme for Conservation Laws)格式,隱式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)方法相結(jié)合,以基于S-A(Spalart-Allmaras)湍流模型的非定常雷諾平均Naiver-Stokes(URANS)方程模擬大范圍氣流分離現(xiàn)象以及翼型前緣脫體渦的流動過程。并通過對HH-02翼型和SC1095翼型動態(tài)失速狀態(tài)的計算,驗證了該數(shù)值方法的有效性。最后著重開展了對帶有后緣小翼的旋翼翼型進行系統(tǒng)性的參數(shù)分析,通過提出的綜合體現(xiàn)翼型升力、阻力及力矩特性的關(guān)系表達式Po和Pc,研究了相對運動頻率、運動規(guī)律間的相位差、偏轉(zhuǎn)幅值等參數(shù)對氣動力特性的影響,獲得了一些有意義的結(jié)論。

    1 網(wǎng)格生成方法

    1.1 翼型網(wǎng)格點重構(gòu)

    在模擬旋翼翼型振蕩的動態(tài)失速特性過程中,本文只考慮平均迎角α0和迎角的一階正弦變化量αm,所以翼型的振蕩規(guī)律可以寫為迎角α=α0+αmsin(2kt),k為減縮頻率。后緣小翼的運動規(guī)律為

    δ=δ0+δmsin(2k*kt-φ0)

    (1)

    式中:δ為翼型后緣小翼瞬時偏轉(zhuǎn)角,當δ為正數(shù)時,表示后緣小翼順時針偏轉(zhuǎn);δ0和δm分別為基準偏轉(zhuǎn)角和偏轉(zhuǎn)幅值,k*為相對于k的運動頻率;φ0為后緣小翼運動規(guī)律與翼型振蕩規(guī)律之間的相位差。

    在對翼型網(wǎng)格修正之前,需要對翼型后緣附近的網(wǎng)格點進行重構(gòu)。本文采用B樣條插值法,對后緣小翼偏轉(zhuǎn)后的翼型進行網(wǎng)格點的重構(gòu),避免小翼上下表面網(wǎng)格點分布過于集中或稀疏對后續(xù)網(wǎng)格生成質(zhì)量的不利影響。圖1給出了翼型振蕩中心和后緣小翼及網(wǎng)格點重構(gòu)的示意圖,從圖中可以分別看出后緣小翼的偏轉(zhuǎn)點、翼型的振蕩中心、后緣小翼的偏轉(zhuǎn)角和翼型的迎角。同時給出了初始翼型和帶有后緣小翼翼型的網(wǎng)格點分布。需要指出的是帶有后緣小翼的旋翼翼型瞬時迎角的定義與初始翼型瞬時迎角的定義相同。

    1.2 運動嵌套網(wǎng)格方法

    考慮到旋翼翼型的振蕩以及后緣小翼的周期性偏轉(zhuǎn)運動,本文發(fā)展了一套適用于帶有后緣小翼的旋翼翼型非定常流動特性模擬的運動嵌套網(wǎng)格方法,并以此研究后緣小翼對旋翼翼型動態(tài)失速的控制效果。主要步驟有:

    1) 圍繞帶有后緣小翼的旋翼翼型,通過求解Poisson方程的方法快速生成高質(zhì)量的貼體網(wǎng)格。計算平面的Poisson控制方程為

    (2)

    式中:ξ和η為網(wǎng)格點坐標,分別代表弦向和法向;γ1、γ2和γ3為坐標變換引入的參數(shù);φ和ψ為用來控制網(wǎng)格分布質(zhì)量的源項。圖2展示了圍繞初始翼型以及帶有后緣小翼翼型生成的貼體網(wǎng)格。

    2) 生成固定不動的笛卡兒背景網(wǎng)格;然后,采用發(fā)展的Top-Map[18]方法確定背景網(wǎng)格在翼型網(wǎng)格上的洞邊界,在此基礎(chǔ)上,通過最小距離法[19]進行背景網(wǎng)格人工內(nèi)邊界的貢獻單元搜索,從而獲得翼型與背景網(wǎng)格的運動嵌套關(guān)系。

    翼型網(wǎng)格邊界單元的貢獻單元由2個一維搜索得到,嵌套網(wǎng)格之間的信息傳遞通過雙線性插值完成。為減小插值過程中引入的數(shù)值誤差,在生成背景網(wǎng)格洞邊界時,使其遠離流場梯度較大的翼型表面。圖3給出了翼型網(wǎng)格運動、網(wǎng)格嵌套關(guān)系及貢獻單元搜索結(jié)果示意圖。

    2 數(shù)值方法

    2.1 非定常流場求解方法

    綜合考慮帶有后緣小翼的振蕩翼型非定常流場的復(fù)雜性,本文采用URANS方程模擬整個流場

    (3)

    式中:W為守恒變量;Fc和Fv分別為修正對流通量和黏性通量;τ和t分別為虛擬時間和物理時間;Ω為單元體的體積;S為網(wǎng)格單元的邊界。

    基于控制體表面的逆變速度Vt和初始的對流通量Fc,0,修正后的對流通量

    Fc=Fc,0-VtW

    (4)

    本文采用格心有限體積法在空間方向上進行離散,對流通量采用Roe-MUSCL格式進行計算,黏性通量采用Jameson二階中心差分方法求解。湍流黏性系數(shù)計算采用S-A湍流模型,因為S-A模型中沒有轉(zhuǎn)捩點,所以對流動做了全湍流假設(shè)。

    為模擬帶有后緣小翼的翼型流場非定常特性,采用雙時間方法進行時間推進。在進行虛擬時間推進時,考慮流場的收斂性,顯式方法的時間步長要求取得較小,這降低了數(shù)值模擬的效率。針對這一問題,虛擬時間步采用隱式LU-SGS格式的方法進行推進,從而有效加大時間步長,提高計算效率。此外,由于網(wǎng)格的位置、形狀均隨時間不斷變化,為了避免由于網(wǎng)格運動引入的額外誤差,幾何守恒律必須得到嚴格的滿足。

    2.2 數(shù)值方法的驗證

    為驗證本文CFD方法的有效性,分別對HH-02翼型在42212狀態(tài)[20]和SC1095翼型在37107狀態(tài)[20]下的動態(tài)失速特性進行了數(shù)值模擬,圖4給出了2種翼型的升力系數(shù)CL、阻力系數(shù)CD和力矩系數(shù)Cm隨迎角α的變化。由圖4可知,計算獲得的翼型升力、阻力和力矩系數(shù)與試驗值吻合較好,特別是對于力矩和阻力系數(shù)的峰值以及再附過程中的升力系數(shù)等的模擬,表明本文的數(shù)值方法能有效地模擬旋翼翼型動態(tài)失速狀態(tài)下的非定常氣動特性。

    3 動態(tài)失速的控制機理及參數(shù)分析

    采用建立的CFD方法開展后緣小翼運動參數(shù)對旋翼翼型動態(tài)失速特性的控制機理研究,參數(shù)包括后緣小翼的相對運動頻率k*、后緣小翼運動規(guī)律與翼型振蕩規(guī)律之間的相位差φ0以及偏轉(zhuǎn)幅值δm。

    為合理地表述旋翼翼型在動態(tài)失速過程中的氣動特性,并能直觀地反映參數(shù)分析中計算結(jié)果的差異性,本文提出了體現(xiàn)翼型升力、阻力及力矩特性的關(guān)系表達式Po和Pc,具體表示為

    (5)

    式中:Nα為一個周期的計算時間步數(shù);CLmax為一個翼型振蕩周期內(nèi)的最大升力系數(shù);KCL為升力系數(shù)線性增長段的斜率。參數(shù)a1、a2、b1、b2為非負權(quán)重系數(shù),并且需要滿足a1+a2=1.0和b1+b2=1.0;σ為一正值常數(shù),用來保證Pc中升力系數(shù)最大值的倒數(shù)項與前一項有相當?shù)闹怠?/p>

    可以看出Po表征了翼型的阻力和力矩特性,由于阻力系數(shù)和力矩系數(shù)在量級上相近,本文中的Po關(guān)系式將阻力系數(shù)和力矩系數(shù)進行了組合,較小的Po代表翼型動態(tài)振蕩過程中有較小的阻力系數(shù)和力矩系數(shù);Pc在一定程度上反映了翼型升力遲滯回線偏離線性增長段的程度及最大升力系數(shù),且翼型升力遲滯回線偏離線性增長段的程度能在較大程度上體現(xiàn)翼型的動態(tài)失速特性,所以Pc較小時表示翼型動態(tài)失速程度較弱。

    Po和Pc關(guān)系表達式中的非負權(quán)重系數(shù)a1、a2、b1、b2,反映了關(guān)系式中各成分所占的比重。對于Po,a1的大小反映了阻力系數(shù)所占的比重,a2的大小反映了力矩系數(shù)所占的比重;對于Pc,b1的大小反映了翼型升力遲滯回線偏離線性增長段的程度(類似于遲滯回線圍成的面積)所占的比重,b2的大小反映了最大升力系數(shù)所占的比重??紤]到翼型的力矩特性和升力遲滯回線偏離線性增長段的程度相對重要,因此所占比重較大。為統(tǒng)一起見,這里表達式中的權(quán)重系數(shù)設(shè)定為a1=0.4、a2=0.6、b1=0.8、b2=0.2和σ=10。

    以SC1095翼型動態(tài)失速的后緣小翼控制作為算例,翼型的振蕩規(guī)律為α=10°+10°sin(2kt),減縮頻率k=0.1,來流馬赫數(shù)為Ma∞=0.3。

    3.1 相對運動頻率的影響

    以后緣小翼的相對運動頻率k*作為研究對象,按照1.1節(jié)給出的翼型和后緣的運動規(guī)律,分別在一系列相對運動頻率下分析后緣小翼對SC1095翼型動態(tài)失速特性控制效果的影響。

    圖5給出了Po和Pc隨后緣小翼相對運動頻率的變化情況。從圖中可以明顯看出,相對運動頻率k*對動態(tài)失速過程中翼型特性的影響并非線性關(guān)系。整體來看,當k*=1時,體現(xiàn)翼型阻力和力矩綜合特性的Po有一個最小值,表明當后緣小翼的偏轉(zhuǎn)頻率與翼型的振蕩頻率相同時,后緣小翼對翼型阻力和力矩發(fā)散的控制效果最好。當k*>4時,隨著k*的增大Po值的變化幅度很小。與Po值的變化情況不同,代表動態(tài)失速過程中翼型升力特性的Pc值的變化情況更為復(fù)雜,Pc值在k*=1時也有一個最小值。從整體趨勢來看,當k*<5時,Pc值會隨著相對運動頻率的增大而增大,這表明翼型升力系數(shù)偏離線性段的幅度也越大,即動態(tài)失速現(xiàn)象越明顯。當k*>5時,Pc值略有下降并趨于平緩。

    圖6給出了在不同相對運動頻率的后緣小翼控制下的翼型阻力和力矩系數(shù)遲滯回線的對比,可以看出,當k*<5時,隨著相對運動頻率的增大,后緣小翼對翼型阻力和力矩系數(shù)發(fā)散的抑制效果逐漸減弱,與前文通過Po和Pc值分析得到的結(jié)果一致,也說明了后緣小翼對翼型阻力和力矩系數(shù)發(fā)散的抑制效果可以很好地通過Po的變化來反映。

    3.2 相位差的影響

    以后緣小翼運動規(guī)律中的φ0作為研究對象,φ0為后緣小翼運動規(guī)律與翼型振蕩規(guī)律之間的相位差。設(shè)置動態(tài)前緣的相對運動頻率k*=1.0,并且在分析時令δ0=0°、δm=10°。分別在一系列相位差φ0下分析后緣小翼對SC1095翼型動態(tài)失速特性控制效果的影響。

    圖7給出了Po和Pc隨后緣小翼相位差的變化情況。從圖中可以明顯看出,相位差φ0對動態(tài)失速過程中翼型特性的影響并非線性關(guān)系。整體來看,Po隨著相位差的增大,呈現(xiàn)先增大后減小的規(guī)律。當φ0=180° 時,體現(xiàn)翼型阻力和力矩綜合特性的Po有一個峰值,表明當后緣小翼的偏轉(zhuǎn)規(guī)律與翼型的振蕩規(guī)律相位差為180° 時,后緣小翼對翼型阻力和力矩發(fā)散的控制效果較差,這是因為當φ0=180° 時,后緣運動規(guī)律與翼型振蕩規(guī)律完全相反,使翼型的彎度變化更為劇烈。當φ0=0°,即兩者運動規(guī)律不存在相位差時,后緣小翼對翼型阻力和力矩發(fā)散的控制效果最好。與Po值的變化情況不同,Pc值受相位差φ0的影響較大,這與Krzysiak和Narkiewicz[12]的試驗結(jié)論:“動態(tài)氣動載荷回線取決于翼型運動與后緣運動的相位差”相一致。從圖中還可以看出,Pc值在φ0=90° 時有一個最大值。從整體趨勢來看,當φ0<90° 時,Pc值會隨著相位差的增大而快速增大,這表明翼型升力系數(shù)偏離線性段的幅度也越大,即動態(tài)失速現(xiàn)象越明顯。當φ0>90° 時,Pc值會隨著相位差的增大而快速減小,直到當φ0>180° 后,Pc值基本保持不變。

    圖8給出了在不同相位差的后緣小翼控制下的翼型阻力和力矩系數(shù)遲滯回線的對比,可以看出,隨著相位差的增大,后緣小翼對翼型阻力和力矩系數(shù)發(fā)散的抑制效果先變小、再變大。當兩者運動規(guī)律不存在相位差時,后緣小翼對翼型阻力和力矩發(fā)散的控制效果最好,這與前文通過Po和Pc值分析得到的結(jié)果相一致。

    3.3 偏轉(zhuǎn)幅值的影響

    進一步開展了后緣小翼偏轉(zhuǎn)幅值δm對翼型動態(tài)失速控制效果影響的數(shù)值分析。設(shè)置后緣小翼的相對運動頻率和相位差分別為k*=1和φ0=0°,并且在分析時令δ0=0°。需要指出的是,負的幅值δm表示翼型后緣順時針偏轉(zhuǎn)。

    圖9給出了Po和Pc的值隨后緣小翼偏轉(zhuǎn)幅值的變化情況。由圖可以看出,當后緣小翼幅值為負時,后緣小翼可能會加劇翼型動態(tài)失速的現(xiàn)象。隨著后緣小翼順時針偏轉(zhuǎn)幅值的增大,翼型動態(tài)失速下的氣動特性會迅速惡化。這主要是因為在較大的迎角條件下,后緣小翼的下偏會增大整個翼型的彎度,從而增大了翼型的阻力和力矩。而當δm為正值時,后緣小翼能夠抑制動態(tài)失速過程中翼型阻力和力矩系數(shù)的發(fā)散,并能夠減少升力系數(shù)在低頭過程中的損失。隨著后緣小翼偏轉(zhuǎn)幅值的增大,Po的值會很快減小然后呈現(xiàn)逐漸平穩(wěn)的趨勢,在δm>5° 時,Po的值趨于穩(wěn)定。Pc的值也是隨著δm的增大而減小,并逐漸趨于平穩(wěn)??傮w而言,后緣小翼對翼型動態(tài)失速的整體控制在δm>0° 時有較優(yōu)的效果。

    為進一步揭示后緣小翼對翼型動態(tài)失速特性的影響規(guī)律,圖10給出了在后緣小翼不同偏轉(zhuǎn)幅值控制下SC1095翼型在一個周期內(nèi)的平均升力系數(shù)和最大阻力系數(shù)、最大力矩系數(shù)隨后緣小翼偏轉(zhuǎn)幅值的變化情況。隨著后緣小翼偏轉(zhuǎn)幅值δm從-10° 增大到10°,翼型在一個周期內(nèi)的平均升力系數(shù)增大,而最大阻力系數(shù)和最大力矩系數(shù)均在降低。

    圖11給出了在后緣小翼不同偏轉(zhuǎn)幅值控制下SC1095翼型阻力系數(shù)和力矩系數(shù)的對比圖。隨著后緣小翼偏轉(zhuǎn)幅值δm從1° 增大到10°,阻力系數(shù)和力矩系數(shù)的峰值降低,但由于在小迎角條件下,翼型的彎度增大,翼型的低頭力矩系數(shù)有所增加;而在大迎角條件下,對翼型的力矩系數(shù)抑制較好。

    3.4 動態(tài)失速的綜合控制分析

    上述不同參數(shù)的影響研究表明,后緣小翼的相對運動減縮頻率k*=1及相位差φ0=0° 時有較好的控制效果,并且當偏轉(zhuǎn)角δm>5° 時,控制效果趨于最佳,為此后緣小翼的偏轉(zhuǎn)規(guī)律選擇了δ=10°sin(2k*kt)。下面給出上述參數(shù)組合影響下的動態(tài)失速控制特性。

    圖12給出了SC1095翼型在有無后緣小翼控制時的翼型阻力、力矩系數(shù)的計算值??梢钥闯?,后緣小翼可以有效地抑制阻力和低頭力矩的發(fā)散,阻力系數(shù)和低頭力矩系數(shù)峰值分別減小了19%和27%。

    為了進一步觀察翼型流場附近的流動細節(jié),分析后緣小翼對翼型周圍流場的影響,圖13進一步給出了有無后緣小翼控制時翼型表面附近的渦量云圖。從圖中可以看出,當翼型迎角超過動態(tài)失速臨界值時,前緣產(chǎn)生脫體渦,并進一步導(dǎo)致翼型前緣的氣流分離,脫體渦沿翼型表面對流直至脫落。通過有無后緣小翼控制時的對比,可以得出后緣小翼并不直接影響脫體渦的形成和再附過程;結(jié)合圖12的結(jié)果進行分析,得出后緣小翼可以降低動態(tài)失速過程中的阻力和力矩系數(shù)峰值。

    4 結(jié) 論

    本文建立了一套適用于帶有后緣小翼控制的旋翼翼型的非定常流場求解方法,可以有效的模擬后緣小翼對旋翼翼型動態(tài)失速的控制作用。并在此基礎(chǔ)上,開展了后緣小翼典型運動參數(shù)對旋翼翼型動態(tài)失速控制效果影響的分析。得到以下結(jié)論。

    1) 建立的包括翼型重構(gòu)、貼體網(wǎng)格生成方法、運動嵌套網(wǎng)格方法和基于URANS方程的非定常流場求解方法在內(nèi)的數(shù)值模擬方法可以有效地模擬帶有后緣小翼控制的旋翼翼型的非定常氣動特性。

    2) 當后緣小翼偏轉(zhuǎn)的相對運動頻率k*=1時,后緣小翼對翼型阻力和力矩系數(shù)發(fā)散的控制效果最好,隨著k*的增大,后緣小翼的控制效果減弱。

    3) 當后緣小翼運動規(guī)律與翼型振蕩規(guī)律同相位時,后緣小翼能夠較好地抑制旋翼翼型動態(tài)失速現(xiàn)象。當后緣小翼順時針偏轉(zhuǎn)時,后緣小翼會加劇翼型動態(tài)失速的現(xiàn)象;當后緣小翼逆時針偏轉(zhuǎn)時,可以較好地抑制翼型動態(tài)失速現(xiàn)象。

    4) 正弦運動形式的后緣小翼雖然不影響前緣渦的形成和再附,但是可以降低阻力和力矩系數(shù)的峰值;同時能夠保證翼型的升力系數(shù)維持在較高的水平,提高了整個周期內(nèi)的平均升力系數(shù)。

    [1] CONLISK A T. Modern helicopter rotor aerodynamics[J]. Progress in Aerospace Sciences, 2002, 37(5): 417-476.

    [2] VIEIRA B A O, MAUGHMER M D. An evaluation of dynamic stall onset prediction methods for rotorcraft airfoil design: AIAA-2013-1093[R]. Reston: AIAA, 2013.

    [3] GEISSLER W, RAFFEL M, DIETZ G, et al. Helicopter aerodynamics with emphasis placed on dynamic stall[M]. Heidelberg: Springer, 2007: 199-204.

    [4] CHOPRA I. Review of state of art of smart structures and integrated systems[J]. AIAA Journal, 2002, 40(11): 2145-2187.

    [5] CARR L. The effect of a leading-edge slat on the dynamic stall of an oscillating airfoil: AIAA-1983-2533[R]. Reston: AIAA, 1983.

    [6] POST M L, CORKE T C. Separation control using plasma actuators: Dynamic stall vortex control on oscillating airfoil[J]. AIAA Journal, 2006, 44(12): 3125-3135.

    [7] ZHAO G Q, ZHAO Q J. Dynamic stall control optimization of rotor airfoil via variable droop leading-edge[J]. Aerospace Science & Technology, 2015, 43(6): 406-414.

    [8] JAWORSKI J W. Thrust and aerodynamic forces from an oscillating leading-edge flap[J]. AIAA Journal, 2012, 50(12): 2928-2931.

    [9] PECHAN T. Design and analysis of an active leading edge wing: AIAA-2013-01222[R]. Reston: AIAA, 2013.

    [10] CHANDRASEKHARA M S. A review of compressible dynamic stall control principles and methods[C]//Proceedings of the Tenth Asian Congress of Fluid Mechanics. ACFM, 2004: 1-6.

    [11] FESZTY D, GILLIES E A, VEZZA M. Alleviation of airfoil dynamic stall moments via trailing-edge flap flow control[J]. AIAA Journal, 2004, 42(1): 17-25.

    [12] KRZYSIAK A, NARKIEWICZ J. Aerodynamic loads on airfoil with trailing-edge flap pitching with different frequencies[J]. Journal of Aircraft, 2006, 43(2): 407-418.

    [13] GERONTAKOS P, LEE T. Trailing-edge flap control of dynamic pitching moment[J]. AIAA Journal, 2012, 45(7): 1688-1694.

    [14] LEE T, SU Y Y. Unsteady airfoil with a harmonically deflected trailing-edge flap[J]. Journal of Fluids & Structures, 2011, 27(8): 1411-1424.

    [15] 王進, 楊茂, 陳鳳明. 帶后緣襟翼翼型的非定常氣動特性數(shù)值仿真[J]. 計算機仿真, 2011, 28(2): 88-92. WANG J, YANG M, CHEN F M. CFD simulation of unsteady aerodynamic of airfoil with trailing-edge flap[J]. Computer Simulation, 2011, 28(2): 88-92 (in Chinese).

    [16] 劉洋, 向錦武. 后緣襟翼對直升機旋翼翼型動態(tài)失速特性的影響[J]. 航空學(xué)報, 2013, 34(5): 1028-1035. LIU Y, XIANG J W. Effect of the trailing edge flap on dynamic stall performance of helicopter rotor airfoil[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(5): 1028-1035 (in Chinese).

    [17] 王榮, 夏品奇. 多片后緣小翼對直升機旋翼槳葉動態(tài)失速及槳轂振動載荷的控制[J]. 航空學(xué)報, 2013, 34(5): 1083-1091. WANG R, XIA P Q. Control of helicopter rotor blade dynamic stall and hub vibration loads by using multiple trailing edge flaps[J]. Acta Aeronautica et Astronautica Sinica, 2013,34(5): 1083-1091 (in Chinese).

    [18] 王博, 招啟軍, 徐廣, 等. 一種適合于旋翼前飛非定常流場計算的新型運動嵌套網(wǎng)格方法[J]. 空氣動力學(xué)學(xué)報, 2012, 30(1): 14-21. WANG B, ZHAO Q J, XU G, et al. A new moving-embedded grid method for numerical simulation of unsteady flow-field of the helicopter rotor in forward flight[J]. Acta Aerodynamica Sinica, 2012, 30(1): 14-21 (in Chinese).

    [19] 趙國慶, 招啟軍, 王清. 旋翼翼型非定常動態(tài)失速特性的CFD模擬及參數(shù)分析[J]. 空氣動力學(xué)學(xué)報, 2015, 33(1): 72-81. ZHAO G Q, ZHAO Q J, WANG Q. Simulations and parametric analyses on unsteady dynamic stall characteristics of rotor airfoil based on CFD method[J]. Acta Aerodynamica Sinica, 2015, 33(1): 72-81 (in Chinese).

    [20] MCALISTER K W, PUCCI S L, MCCROSKEY W J, et al. An experimental study of dynamic stall on advanced airfoil section. Volume 2: Pressure and force data: NASA TM-84245[R]. Washington, D.C.: NASA, 1982.

    (責(zé)任編輯:鮑亞平, 張晗)

    *Corresponding author. E-mail: zhaoqijun@nuaa.edu.cn

    Dynamic stall control of rotor airfoil via trailing-edge flap

    MA Yiyang1, ZHAO Qijun1,*, ZHAO Guoqing2

    1.NationalKeyLaboratoryofScienceandTechnologyonRotorcraftAeromechanics,NanjingUniversityofAeronauticsandAstronautics,Nanjing210016,China2.GeneralandAerodynamicInstitute,TheFirstAircraftInstituteofAVIC,Xi’an710089,China

    Control effects of typical motion parameters of trailing-edge flap (TEF) on the dynamic stall characteristics of rotor airfoil are investigated. A high-efficiency and high-precision CFD method for predicting the unsteady flow characteristics of rotor airfoil with TEF control is developed. The viscous and orthogonal body-fitted grids around the oscillatory rotor airfoil are regenerated by solving Poisson equations. To ensure the quality of the grids around TEF, a reconstruction of grid points on airfoil is conducted to describe the motion of TEF. Aiming at overcoming the shortcoming of deformable grid approach, which may result in the distortion of grid, a moving-embedded grid method is developed to predict the flowfield of the oscillatory airfoil with TEF control. Based on unsteady Reynolds averaged Navier-Stokes (URANS) equations, dual-time method, Spalart-Allmaras (S-A) turbulence model, and Roe-Monotone Upwind-centered Scheme for Conservation Laws (Roe-MUSCL) scheme, a high-precision CFD method for predicting the flowfield around airfoil is developed, and implicit Lower-Upper Symmetric Gauss-Seidel (LU-SGS) scheme and parallel techniques are adopted to improve computational efficiency. The dynamic stall cases of HH-02 and SC1095 rotor airfoils are calculated using the proposed method. It is demonstrated that the hysteresis effects are well captured, verifying the effectiveness of numerical simulation method. Focusing on the dynamic stall cases of SC1095 rotor airfoil, analyses on dynamic stall control via TEF are carried out. The functionPoandPcwhich can reflect the lift, drag and moment characteristics of airfoil are proposed. The effects of the non-dimensional frequency, the phase difference and the angular amplitude of the trailing-edge flap are revealed. The results indicate that dynamic stall phenomenon of an oscillatory airfoil could be significantly suppressed when relative motion frequency of the trailing-edge flap is 1.0 and the phase difference is about 0°. At this state, the maximum drag and negative moment coefficients of SC1095 airfoil can be reduced by about 19% and 27% respectively via TEF control when the angular amplitude is 10°.

    rotor; airfoil; dynamic stall; trailing-edge flap; parametric analysis; unsteady Reynolds averaged Navier-Stokes equation

    2016-04-11; Revised:2016-06-12; Accepted:2016-07-16; Published online:2016-08-02 09:49

    URL:www.cnki.net/kcms/detail/11.1929.V.20160802.0949.002.html

    s:National Natural Science Foundation of China (11272150, 11572156); Jiangsu Innovation Program for Graduate Education (KYLX15_0244); Priority Academic Program Development of Jiangsu Higher Education Institutions

    http://hkxb.buaa.edu.cn hkxb@buaa.edu.cn

    10.7527/S1000-6893.2016.0220

    2016-04-11; 退修日期:2016-06-12; 錄用日期:2016-07-16; 網(wǎng)絡(luò)出版時間:2016-08-02 09:49

    www.cnki.net/kcms/detail/11.1929.V.20160802.0949.002.html

    國家自然科學(xué)基金 (11272150, 11572156); 江蘇省普通高校研究生科研創(chuàng)新計劃項目 (KYLX15_0244); 江蘇高校優(yōu)勢學(xué)科建設(shè)工程基金

    *通訊作者.E-mail: zhaoqijun@nuaa.edu.cn

    馬奕揚, 招啟軍, 趙國慶. 基于后緣小翼的旋翼翼型動態(tài)失速控制分析[J]. 航空學(xué)報, 2017, 38(3): 120312. MA Y Y, ZHAO Q J, ZHAO G Q. Dynamic stall control of rotor airfoil via trailing-edge flap[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(3): 120312.

    V211.3, V211.52

    A

    1000-6893(2017)03-120312-11

    猜你喜歡
    小翼后緣升力
    漢字獵人(一)
    我家養(yǎng)了一只紙精靈(二)
    我家養(yǎng)了一只紙精靈(四)
    高速列車車頂–升力翼組合體氣動特性
    無人機升力測試裝置設(shè)計及誤差因素分析
    基于自適應(yīng)偽譜法的升力式飛行器火星進入段快速軌跡優(yōu)化
    我是霸王龍
    機翼后緣連續(xù)變彎度對客機氣動特性影響
    柔性后緣可變形機翼氣動特性分析
    升力式再入飛行器體襟翼姿態(tài)控制方法
    欧美性长视频在线观看| 国产成+人综合+亚洲专区| 九色国产91popny在线| 日韩大尺度精品在线看网址| 精品不卡国产一区二区三区| 男人舔女人下体高潮全视频| 91成年电影在线观看| 午夜精品在线福利| 欧美久久黑人一区二区| 亚洲欧美精品综合一区二区三区| 97人妻精品一区二区三区麻豆 | 十八禁人妻一区二区| 中出人妻视频一区二区| 免费在线观看视频国产中文字幕亚洲| 免费看a级黄色片| 国产精品日韩av在线免费观看| 给我免费播放毛片高清在线观看| 日本一本二区三区精品| 欧美三级亚洲精品| 免费av毛片视频| 亚洲国产精品成人综合色| 久久久久久亚洲精品国产蜜桃av| av电影中文网址| 婷婷精品国产亚洲av在线| 在线免费观看的www视频| 99久久无色码亚洲精品果冻| 国产精品永久免费网站| 久久香蕉精品热| 在线永久观看黄色视频| 啪啪无遮挡十八禁网站| 久久久水蜜桃国产精品网| 国产激情偷乱视频一区二区| 久久伊人香网站| 99热这里只有精品一区 | 日本撒尿小便嘘嘘汇集6| 69av精品久久久久久| 亚洲,欧美精品.| 日本精品一区二区三区蜜桃| 亚洲中文av在线| 精品国产乱码久久久久久男人| 1024手机看黄色片| 亚洲 欧美一区二区三区| 国产亚洲精品久久久久久毛片| 免费av毛片视频| 中国美女看黄片| 日本在线视频免费播放| 亚洲精品在线美女| 亚洲精品国产精品久久久不卡| 十八禁人妻一区二区| 午夜久久久久精精品| 精品少妇一区二区三区视频日本电影| 国产91精品成人一区二区三区| 久久天躁狠狠躁夜夜2o2o| 亚洲国产中文字幕在线视频| 在线观看免费视频日本深夜| 一进一出好大好爽视频| 大香蕉久久成人网| 黄色毛片三级朝国网站| 国产野战对白在线观看| 少妇裸体淫交视频免费看高清 | 法律面前人人平等表现在哪些方面| 亚洲va日本ⅴa欧美va伊人久久| 久久久国产欧美日韩av| 人成视频在线观看免费观看| 制服人妻中文乱码| 老熟妇仑乱视频hdxx| av有码第一页| 久久久国产成人精品二区| 久久久久免费精品人妻一区二区 | 午夜a级毛片| 一夜夜www| 婷婷丁香在线五月| 亚洲欧美精品综合一区二区三区| 国产精品一区二区三区四区久久 | 日本成人三级电影网站| 亚洲成人国产一区在线观看| 在线观看日韩欧美| 亚洲一区二区三区不卡视频| 亚洲第一av免费看| 99久久无色码亚洲精品果冻| 国内精品久久久久精免费| 亚洲电影在线观看av| 免费在线观看完整版高清| 在线观看免费日韩欧美大片| 精品久久久久久久人妻蜜臀av| 午夜福利高清视频| 搡老妇女老女人老熟妇| 在线播放国产精品三级| 午夜福利高清视频| 久久久精品欧美日韩精品| 18禁黄网站禁片午夜丰满| 亚洲中文日韩欧美视频| 精华霜和精华液先用哪个| 亚洲七黄色美女视频| 色综合欧美亚洲国产小说| 激情在线观看视频在线高清| 久久久久久久久久黄片| 国产麻豆成人av免费视频| 在线国产一区二区在线| 国内揄拍国产精品人妻在线 | 不卡av一区二区三区| 久久狼人影院| 成人手机av| 亚洲精品av麻豆狂野| 日韩精品青青久久久久久| 国产精品av久久久久免费| 少妇的丰满在线观看| 欧美不卡视频在线免费观看 | 黄片大片在线免费观看| 一级毛片高清免费大全| 欧美在线黄色| 真人一进一出gif抽搐免费| 欧美最黄视频在线播放免费| 91成人精品电影| 大型av网站在线播放| 国内久久婷婷六月综合欲色啪| 在线观看免费视频日本深夜| 欧美性猛交黑人性爽| 日本撒尿小便嘘嘘汇集6| 两性午夜刺激爽爽歪歪视频在线观看 | 国产野战对白在线观看| 亚洲精品久久国产高清桃花| 可以在线观看毛片的网站| 午夜免费成人在线视频| 久久精品成人免费网站| 亚洲欧美日韩高清在线视频| 2021天堂中文幕一二区在线观 | 亚洲一区二区三区色噜噜| 国产精品日韩av在线免费观看| 国产亚洲欧美98| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品色激情综合| 久9热在线精品视频| 亚洲av片天天在线观看| 男女做爰动态图高潮gif福利片| 桃红色精品国产亚洲av| 国产成+人综合+亚洲专区| 日本 av在线| 大香蕉久久成人网| 中文字幕精品免费在线观看视频| 亚洲国产看品久久| 欧美乱妇无乱码| 色播在线永久视频| 日本 av在线| 制服诱惑二区| 亚洲最大成人中文| 国产午夜精品久久久久久| 欧美色欧美亚洲另类二区| 亚洲一区中文字幕在线| 哪里可以看免费的av片| 国产亚洲欧美精品永久| 亚洲 欧美 日韩 在线 免费| 男女之事视频高清在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产一区在线观看成人免费| 日韩有码中文字幕| 激情在线观看视频在线高清| 欧美性猛交黑人性爽| 婷婷丁香在线五月| 亚洲精品中文字幕在线视频| 亚洲av第一区精品v没综合| 国产精品永久免费网站| 最好的美女福利视频网| 久久性视频一级片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产黄a三级三级三级人| 夜夜夜夜夜久久久久| 激情在线观看视频在线高清| 欧美性猛交黑人性爽| 国产区一区二久久| 亚洲在线自拍视频| 中文字幕精品免费在线观看视频| 亚洲精品国产区一区二| av在线播放免费不卡| 19禁男女啪啪无遮挡网站| 午夜福利视频1000在线观看| 免费在线观看成人毛片| 成人欧美大片| 国产一级毛片七仙女欲春2 | 午夜精品在线福利| 国产1区2区3区精品| 国产亚洲精品第一综合不卡| 制服诱惑二区| 欧美精品亚洲一区二区| 哪里可以看免费的av片| 国产av一区在线观看免费| 午夜视频精品福利| 免费搜索国产男女视频| 91九色精品人成在线观看| 婷婷六月久久综合丁香| 国产日本99.免费观看| 国产精品一区二区精品视频观看| 国产一区二区三区在线臀色熟女| 午夜久久久在线观看| 欧美亚洲日本最大视频资源| 欧美在线一区亚洲| 久久精品国产综合久久久| 少妇被粗大的猛进出69影院| 国产成人欧美在线观看| 12—13女人毛片做爰片一| 黄频高清免费视频| 看片在线看免费视频| 一区福利在线观看| 在线av久久热| ponron亚洲| 国产高清有码在线观看视频 | 窝窝影院91人妻| 韩国精品一区二区三区| 在线观看免费日韩欧美大片| av电影中文网址| 精品乱码久久久久久99久播| 亚洲欧美精品综合久久99| 不卡一级毛片| 18美女黄网站色大片免费观看| 一区二区三区国产精品乱码| 黄片小视频在线播放| av片东京热男人的天堂| 久久婷婷成人综合色麻豆| 精品免费久久久久久久清纯| 麻豆一二三区av精品| 巨乳人妻的诱惑在线观看| 精品欧美国产一区二区三| 欧美性长视频在线观看| 国产一区二区三区在线臀色熟女| 啦啦啦观看免费观看视频高清| 国产精品久久久av美女十八| 可以在线观看的亚洲视频| 国产精品99久久99久久久不卡| 欧美激情极品国产一区二区三区| 国内久久婷婷六月综合欲色啪| 大型黄色视频在线免费观看| 禁无遮挡网站| 国语自产精品视频在线第100页| 午夜福利免费观看在线| 色尼玛亚洲综合影院| 波多野结衣高清无吗| 丝袜人妻中文字幕| 亚洲国产精品合色在线| 久久人妻福利社区极品人妻图片| 午夜免费鲁丝| 精品国产美女av久久久久小说| 精品一区二区三区四区五区乱码| 不卡一级毛片| 俺也久久电影网| 好看av亚洲va欧美ⅴa在| 久久中文看片网| 久久草成人影院| 午夜福利高清视频| 制服诱惑二区| 成人特级黄色片久久久久久久| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av中文字字幕乱码综合 | 国产精品野战在线观看| 日本精品一区二区三区蜜桃| 国产人伦9x9x在线观看| 在线观看一区二区三区| 少妇的丰满在线观看| 国产高清视频在线播放一区| 欧美乱妇无乱码| 曰老女人黄片| 午夜福利免费观看在线| 欧美一级a爱片免费观看看 | tocl精华| 午夜两性在线视频| 久久精品91无色码中文字幕| 日本熟妇午夜| 久久中文字幕人妻熟女| 亚洲精品av麻豆狂野| 久久天堂一区二区三区四区| 日本三级黄在线观看| 女人高潮潮喷娇喘18禁视频| 热99re8久久精品国产| 日本黄色视频三级网站网址| 亚洲五月天丁香| 又黄又粗又硬又大视频| 国产精品精品国产色婷婷| 我的亚洲天堂| 亚洲第一青青草原| 国产人伦9x9x在线观看| 久久婷婷人人爽人人干人人爱| 十分钟在线观看高清视频www| 精品一区二区三区视频在线观看免费| 国产成人影院久久av| 熟妇人妻久久中文字幕3abv| 久久中文字幕人妻熟女| 精品卡一卡二卡四卡免费| 日韩欧美三级三区| 成人18禁高潮啪啪吃奶动态图| 精品国产一区二区三区四区第35| 国产单亲对白刺激| 色在线成人网| 欧美黄色淫秽网站| 午夜影院日韩av| 一进一出抽搐动态| 亚洲专区中文字幕在线| 免费在线观看黄色视频的| 熟女电影av网| 又大又爽又粗| 国产乱人伦免费视频| 亚洲九九香蕉| 在线观看舔阴道视频| 啦啦啦 在线观看视频| 国产av一区在线观看免费| 18禁黄网站禁片免费观看直播| 国产成人av教育| 亚洲成人久久性| 99riav亚洲国产免费| 美女国产高潮福利片在线看| 日日干狠狠操夜夜爽| 中文资源天堂在线| 嫩草影院精品99| 久久 成人 亚洲| 88av欧美| 亚洲欧美日韩高清在线视频| www.999成人在线观看| 国产av一区在线观看免费| www日本黄色视频网| 国产日本99.免费观看| 亚洲专区中文字幕在线| 一进一出抽搐动态| 日韩精品青青久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人欧美| 麻豆国产av国片精品| 国产亚洲av高清不卡| 国产成人精品久久二区二区91| 亚洲一区二区三区不卡视频| 少妇裸体淫交视频免费看高清 | 国产精品98久久久久久宅男小说| 久久久久久久精品吃奶| 免费无遮挡裸体视频| 最近最新中文字幕大全免费视频| 欧美 亚洲 国产 日韩一| 听说在线观看完整版免费高清| 亚洲av中文字字幕乱码综合 | 超碰成人久久| 搡老妇女老女人老熟妇| 色婷婷久久久亚洲欧美| 大型av网站在线播放| 搞女人的毛片| 亚洲av熟女| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久久精品吃奶| 白带黄色成豆腐渣| 午夜福利在线观看吧| 午夜激情av网站| 亚洲欧美一区二区三区黑人| 超碰成人久久| 大型av网站在线播放| 欧美性长视频在线观看| 69av精品久久久久久| 757午夜福利合集在线观看| 男女之事视频高清在线观看| 久久久久国产精品人妻aⅴ院| 国内毛片毛片毛片毛片毛片| 级片在线观看| 婷婷精品国产亚洲av| 波多野结衣av一区二区av| 叶爱在线成人免费视频播放| 国产精品永久免费网站| 丝袜在线中文字幕| 亚洲黑人精品在线| 中文字幕久久专区| 国产野战对白在线观看| 久久精品国产亚洲av高清一级| 午夜免费鲁丝| 伊人久久大香线蕉亚洲五| videosex国产| 成年免费大片在线观看| 日韩欧美 国产精品| 欧美人与性动交α欧美精品济南到| 不卡一级毛片| 亚洲成人久久爱视频| 18禁观看日本| 老司机午夜十八禁免费视频| 日韩成人在线观看一区二区三区| 成人三级做爰电影| 草草在线视频免费看| 欧洲精品卡2卡3卡4卡5卡区| 色综合亚洲欧美另类图片| 热re99久久国产66热| 亚洲精品国产精品久久久不卡| 一级毛片精品| 亚洲无线在线观看| 无限看片的www在线观看| 精品国产乱码久久久久久男人| 无限看片的www在线观看| 免费女性裸体啪啪无遮挡网站| 日韩国内少妇激情av| 18禁美女被吸乳视频| 亚洲国产精品合色在线| 欧美黑人巨大hd| 一级黄色大片毛片| 两性夫妻黄色片| 亚洲成人久久性| bbb黄色大片| 午夜福利在线在线| 精品欧美一区二区三区在线| 欧美黑人精品巨大| 国产视频内射| 成人av一区二区三区在线看| 十八禁人妻一区二区| 国产精品免费视频内射| 黄频高清免费视频| 99国产精品99久久久久| 黄色视频不卡| 热re99久久国产66热| 听说在线观看完整版免费高清| 一本综合久久免费| 最近最新中文字幕大全免费视频| 侵犯人妻中文字幕一二三四区| 丝袜人妻中文字幕| 一区二区三区高清视频在线| 国产日本99.免费观看| 欧美av亚洲av综合av国产av| 国产成人啪精品午夜网站| 国产精品久久久久久精品电影 | a级毛片a级免费在线| 成人精品一区二区免费| 黄频高清免费视频| 又紧又爽又黄一区二区| 亚洲色图 男人天堂 中文字幕| 日韩欧美 国产精品| 亚洲人成伊人成综合网2020| 国产真人三级小视频在线观看| 亚洲欧美精品综合久久99| 婷婷六月久久综合丁香| 人妻丰满熟妇av一区二区三区| 久久久久久久午夜电影| 不卡av一区二区三区| 夜夜爽天天搞| 欧美精品啪啪一区二区三区| a在线观看视频网站| 精品一区二区三区视频在线观看免费| 久久久国产成人精品二区| 男人舔女人的私密视频| 国产熟女午夜一区二区三区| 久久久久久免费高清国产稀缺| 在线观看免费午夜福利视频| 无遮挡黄片免费观看| 欧美另类亚洲清纯唯美| 在线免费观看的www视频| 人妻丰满熟妇av一区二区三区| 亚洲av日韩精品久久久久久密| 9191精品国产免费久久| 琪琪午夜伦伦电影理论片6080| 观看免费一级毛片| www日本黄色视频网| 老司机福利观看| 欧美在线黄色| 97超级碰碰碰精品色视频在线观看| 午夜视频精品福利| 国产精华一区二区三区| av免费在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 黑人操中国人逼视频| 欧美黑人欧美精品刺激| 黄色丝袜av网址大全| 男女下面进入的视频免费午夜 | 日日摸夜夜添夜夜添小说| 久久亚洲精品不卡| aaaaa片日本免费| 日韩欧美 国产精品| 俺也久久电影网| 亚洲av五月六月丁香网| av在线天堂中文字幕| 国产激情欧美一区二区| 午夜福利在线在线| 欧美日韩福利视频一区二区| 美女国产高潮福利片在线看| 黑丝袜美女国产一区| 亚洲第一电影网av| 国产在线观看jvid| 久久久精品国产亚洲av高清涩受| 亚洲欧美精品综合一区二区三区| 男女视频在线观看网站免费 | 美女高潮喷水抽搐中文字幕| 久久精品国产综合久久久| 免费人成视频x8x8入口观看| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利欧美成人| 可以免费在线观看a视频的电影网站| 一a级毛片在线观看| 制服人妻中文乱码| 国产精品自产拍在线观看55亚洲| 黄色成人免费大全| 免费在线观看黄色视频的| 国产在线观看jvid| 精品久久久久久久久久免费视频| 亚洲九九香蕉| 亚洲三区欧美一区| 亚洲成a人片在线一区二区| 男女下面进入的视频免费午夜 | 国产蜜桃级精品一区二区三区| 日本一本二区三区精品| 日韩欧美三级三区| 中文在线观看免费www的网站 | 亚洲精品色激情综合| 婷婷丁香在线五月| 韩国精品一区二区三区| 久久久国产成人免费| 亚洲精品国产区一区二| 免费搜索国产男女视频| 欧美日本视频| 午夜福利免费观看在线| 夜夜躁狠狠躁天天躁| 日韩欧美一区视频在线观看| 老司机午夜十八禁免费视频| 亚洲精品中文字幕在线视频| 免费看十八禁软件| 日韩精品中文字幕看吧| 欧美日韩乱码在线| netflix在线观看网站| 91成年电影在线观看| 精品国产乱子伦一区二区三区| 19禁男女啪啪无遮挡网站| 女性被躁到高潮视频| 男人操女人黄网站| 国产精品影院久久| 级片在线观看| 极品教师在线免费播放| 亚洲五月色婷婷综合| 精品无人区乱码1区二区| 久久婷婷人人爽人人干人人爱| 国产激情偷乱视频一区二区| 亚洲国产欧美一区二区综合| 亚洲专区中文字幕在线| 久久伊人香网站| 成人亚洲精品一区在线观看| 亚洲专区国产一区二区| aaaaa片日本免费| 午夜免费观看网址| 91在线观看av| a在线观看视频网站| 最好的美女福利视频网| 精品第一国产精品| 黄频高清免费视频| 亚洲国产欧美日韩在线播放| 欧美成人免费av一区二区三区| 亚洲成人国产一区在线观看| 夜夜躁狠狠躁天天躁| 一二三四社区在线视频社区8| xxx96com| 在线观看一区二区三区| 免费看日本二区| 久久国产精品男人的天堂亚洲| 国产精品久久电影中文字幕| av在线播放免费不卡| 欧美一级毛片孕妇| 国产欧美日韩一区二区三| 亚洲人成网站高清观看| 黄色视频不卡| netflix在线观看网站| 免费观看人在逋| www日本在线高清视频| 99国产精品一区二区蜜桃av| 国产免费男女视频| 老熟妇乱子伦视频在线观看| 一级片免费观看大全| 在线观看66精品国产| 青草久久国产| 亚洲 国产 在线| e午夜精品久久久久久久| 久久久国产欧美日韩av| 国产成人欧美在线观看| 久久狼人影院| 国产aⅴ精品一区二区三区波| 亚洲熟妇中文字幕五十中出| 亚洲av中文字字幕乱码综合 | 国产99久久九九免费精品| 精品久久久久久久毛片微露脸| 国产熟女xx| 日本a在线网址| 欧美精品亚洲一区二区| 午夜福利视频1000在线观看| 成人av一区二区三区在线看| 中文字幕精品亚洲无线码一区 | 亚洲精品中文字幕一二三四区| 少妇粗大呻吟视频| 在线播放国产精品三级| 少妇裸体淫交视频免费看高清 | 精华霜和精华液先用哪个| 国产99久久九九免费精品| 日本 欧美在线| 男女午夜视频在线观看| 中文字幕久久专区| 99久久无色码亚洲精品果冻| 韩国精品一区二区三区| 黑人欧美特级aaaaaa片| 亚洲一码二码三码区别大吗| 此物有八面人人有两片| 亚洲人成网站高清观看| 亚洲国产精品成人综合色| 99国产极品粉嫩在线观看| 黄色a级毛片大全视频| 99久久无色码亚洲精品果冻| 人人妻人人澡欧美一区二区| 免费在线观看日本一区| 午夜福利在线观看吧| 日韩成人在线观看一区二区三区| 亚洲精品av麻豆狂野| 这个男人来自地球电影免费观看| 国产精品久久久av美女十八| 麻豆一二三区av精品| 大香蕉久久成人网| 国产精品香港三级国产av潘金莲| 9191精品国产免费久久| 真人一进一出gif抽搐免费| 久久精品国产亚洲av香蕉五月| 久久精品亚洲精品国产色婷小说| 欧美日韩一级在线毛片| 国产91精品成人一区二区三区| 欧美国产精品va在线观看不卡| 中文字幕精品亚洲无线码一区 | 正在播放国产对白刺激|