• <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)控制方法
    中文字幕精品免费在线观看视频| 丁香六月欧美| 级片在线观看| 日韩一卡2卡3卡4卡2021年| 日韩大码丰满熟妇| 亚洲九九香蕉| 99在线人妻在线中文字幕| 午夜福利一区二区在线看| 99国产极品粉嫩在线观看| 老汉色∧v一级毛片| 曰老女人黄片| 亚洲av片天天在线观看| 免费观看精品视频网站| 国产精品一区二区在线不卡| 熟女少妇亚洲综合色aaa.| 午夜精品久久久久久毛片777| 两性夫妻黄色片| 51午夜福利影视在线观看| 一个人免费在线观看的高清视频| 天堂动漫精品| 久久精品影院6| 国产私拍福利视频在线观看| 国产精品久久久人人做人人爽| av中文乱码字幕在线| 欧美+亚洲+日韩+国产| 国产亚洲av高清不卡| 午夜激情av网站| 99国产精品一区二区三区| 亚洲国产欧美网| 啦啦啦韩国在线观看视频| 国产精品永久免费网站| 日韩欧美在线二视频| 熟女少妇亚洲综合色aaa.| 99久久综合精品五月天人人| 国产区一区二久久| 精品久久蜜臀av无| 成年女人毛片免费观看观看9| 性少妇av在线| 一进一出抽搐动态| 久久九九热精品免费| 日韩高清综合在线| 中亚洲国语对白在线视频| 男女下面插进去视频免费观看| 久久人妻福利社区极品人妻图片| 日韩视频一区二区在线观看| 久久久久久国产a免费观看| 在线天堂中文资源库| 亚洲中文av在线| 国产一区二区在线av高清观看| 国产精品二区激情视频| 国产精品一区二区三区四区久久 | 搞女人的毛片| 又黄又爽又免费观看的视频| 一边摸一边抽搐一进一出视频| 中文字幕人妻丝袜一区二区| av视频免费观看在线观看| 国产不卡一卡二| 啦啦啦免费观看视频1| 亚洲国产看品久久| 国产一区在线观看成人免费| 亚洲国产精品sss在线观看| 亚洲欧美日韩无卡精品| 男女午夜视频在线观看| 亚洲全国av大片| 成人国产一区最新在线观看| www.自偷自拍.com| 亚洲少妇的诱惑av| 亚洲 国产 在线| 嫩草影视91久久| 日韩精品青青久久久久久| 久久久久亚洲av毛片大全| 久久久久亚洲av毛片大全| 国产三级在线视频| 99国产极品粉嫩在线观看| 无遮挡黄片免费观看| 久久人人精品亚洲av| 国产av在哪里看| 国产亚洲精品第一综合不卡| 村上凉子中文字幕在线| 日韩成人在线观看一区二区三区| 久久人人精品亚洲av| 丁香欧美五月| 一区福利在线观看| 99re在线观看精品视频| 一区福利在线观看| 精品国产美女av久久久久小说| 欧美午夜高清在线| 国产精品二区激情视频| 亚洲精品国产精品久久久不卡| av有码第一页| 国内毛片毛片毛片毛片毛片| 满18在线观看网站| av网站免费在线观看视频| 丝袜人妻中文字幕| 欧美成狂野欧美在线观看| 久久久国产成人免费| 人人妻人人澡人人看| 免费在线观看亚洲国产| 桃色一区二区三区在线观看| 久久久久国产一级毛片高清牌| 老司机福利观看| 多毛熟女@视频| 9色porny在线观看| 久久 成人 亚洲| 亚洲专区字幕在线| 日韩精品中文字幕看吧| 亚洲专区国产一区二区| 精品一区二区三区视频在线观看免费| 国产高清视频在线播放一区| 久久人人97超碰香蕉20202| 国产av在哪里看| 日本 av在线| 黄片小视频在线播放| 天堂动漫精品| 国内毛片毛片毛片毛片毛片| 一个人观看的视频www高清免费观看 | 国产高清videossex| 亚洲情色 制服丝袜| 丁香六月欧美| 国产伦一二天堂av在线观看| 亚洲一区中文字幕在线| 国产亚洲精品久久久久久毛片| 亚洲精品中文字幕在线视频| 好男人电影高清在线观看| 色在线成人网| 免费av毛片视频| 岛国视频午夜一区免费看| 国产亚洲欧美在线一区二区| 日韩视频一区二区在线观看| 满18在线观看网站| 一本久久中文字幕| 欧美人与性动交α欧美精品济南到| 欧美色欧美亚洲另类二区 | av天堂在线播放| 精品国产一区二区三区四区第35| 一进一出好大好爽视频| 长腿黑丝高跟| 亚洲人成77777在线视频| 黄频高清免费视频| 一边摸一边抽搐一进一小说| 午夜精品在线福利| 国产成人av教育| 91麻豆av在线| 久久久久精品国产欧美久久久| 成人18禁高潮啪啪吃奶动态图| 脱女人内裤的视频| 精品久久久精品久久久| 国产精品98久久久久久宅男小说| 久久婷婷成人综合色麻豆| 91成人精品电影| avwww免费| 美女高潮到喷水免费观看| 午夜亚洲福利在线播放| 亚洲五月婷婷丁香| 十分钟在线观看高清视频www| 亚洲男人的天堂狠狠| av片东京热男人的天堂| 亚洲第一欧美日韩一区二区三区| 级片在线观看| 一区二区三区国产精品乱码| 老熟妇乱子伦视频在线观看| 纯流量卡能插随身wifi吗| 巨乳人妻的诱惑在线观看| 宅男免费午夜| 国产精品爽爽va在线观看网站 | av中文乱码字幕在线| 日韩欧美一区视频在线观看| 亚洲av日韩精品久久久久久密| 欧美在线一区亚洲| 午夜免费激情av| 久久久久国产一级毛片高清牌| 高潮久久久久久久久久久不卡| 不卡av一区二区三区| 黑人巨大精品欧美一区二区mp4| 亚洲精品中文字幕在线视频| 国产真人三级小视频在线观看| 天堂动漫精品| 免费女性裸体啪啪无遮挡网站| 18禁裸乳无遮挡免费网站照片 | 视频区欧美日本亚洲| 亚洲精品中文字幕一二三四区| 高清毛片免费观看视频网站| 久久这里只有精品19| 亚洲av电影不卡..在线观看| 国产人伦9x9x在线观看| 欧美乱码精品一区二区三区| 侵犯人妻中文字幕一二三四区| 亚洲 国产 在线| 免费av毛片视频| 久久人妻熟女aⅴ| 男人操女人黄网站| 老司机福利观看| 一级毛片高清免费大全| 又黄又爽又免费观看的视频| 亚洲av成人不卡在线观看播放网| 国产主播在线观看一区二区| 亚洲无线在线观看| 午夜福利在线观看吧| 在线观看日韩欧美| 伦理电影免费视频| 精品久久久精品久久久| 国产99白浆流出| 欧美在线黄色| 国产av精品麻豆| 欧美成狂野欧美在线观看| 在线观看66精品国产| 亚洲中文字幕一区二区三区有码在线看 | 亚洲专区字幕在线| 亚洲人成网站在线播放欧美日韩| 精品人妻1区二区| 在线免费观看的www视频| 大陆偷拍与自拍| 午夜精品在线福利| www.精华液| 日日摸夜夜添夜夜添小说| 亚洲av五月六月丁香网| 999精品在线视频| 久久狼人影院| 国产精品精品国产色婷婷| 91国产中文字幕| 熟女少妇亚洲综合色aaa.| netflix在线观看网站| 精品一区二区三区视频在线观看免费| 午夜福利欧美成人| 色综合站精品国产| 最近最新免费中文字幕在线| 国产区一区二久久| 久久青草综合色| 女同久久另类99精品国产91| 美女免费视频网站| 变态另类丝袜制服| 一级片免费观看大全| 极品人妻少妇av视频| 国产av一区在线观看免费| 9191精品国产免费久久| 亚洲电影在线观看av| 搡老岳熟女国产| or卡值多少钱| 亚洲人成电影观看| 欧美日韩亚洲国产一区二区在线观看| 欧美最黄视频在线播放免费| 亚洲人成网站在线播放欧美日韩| 巨乳人妻的诱惑在线观看| 久久中文字幕一级| 女警被强在线播放| 成人18禁高潮啪啪吃奶动态图| 亚洲一区中文字幕在线| 精品人妻1区二区| 女生性感内裤真人,穿戴方法视频| 免费看十八禁软件| 午夜福利成人在线免费观看| 亚洲专区字幕在线| 国产精品久久久久久亚洲av鲁大| 99精品欧美一区二区三区四区| 久久久久久久精品吃奶| e午夜精品久久久久久久| 看免费av毛片| 亚洲三区欧美一区| 男女做爰动态图高潮gif福利片 | av福利片在线| 精品第一国产精品| 国语自产精品视频在线第100页| 午夜福利欧美成人| av天堂在线播放| 在线天堂中文资源库| 18禁裸乳无遮挡免费网站照片 | 亚洲国产精品成人综合色| 精品国产一区二区久久| 精品久久久久久久久久免费视频| 亚洲欧美日韩另类电影网站| 亚洲国产精品sss在线观看| 亚洲九九香蕉| 夜夜夜夜夜久久久久| 怎么达到女性高潮| 悠悠久久av| 久久久水蜜桃国产精品网| 波多野结衣一区麻豆| 人人澡人人妻人| 国产亚洲精品久久久久久毛片| 日韩欧美一区视频在线观看| 麻豆av在线久日| √禁漫天堂资源中文www| 成年版毛片免费区| 在线视频色国产色| 亚洲人成77777在线视频| 精品熟女少妇八av免费久了| 高潮久久久久久久久久久不卡| 亚洲成人免费电影在线观看| 在线观看免费视频网站a站| 欧美在线一区亚洲| 中文字幕色久视频| 欧美乱码精品一区二区三区| 亚洲伊人色综图| 国产成人精品在线电影| 老司机午夜十八禁免费视频| 精品国内亚洲2022精品成人| 国产成人欧美| 亚洲久久久国产精品| 99riav亚洲国产免费| 九色国产91popny在线| 精品不卡国产一区二区三区| 国产av一区二区精品久久| 久久精品国产亚洲av高清一级| www国产在线视频色| 免费在线观看视频国产中文字幕亚洲| 99国产极品粉嫩在线观看| 国产av精品麻豆| 成人国语在线视频| 色播在线永久视频| 好男人电影高清在线观看| 90打野战视频偷拍视频| 国产精品秋霞免费鲁丝片| 日日夜夜操网爽| 国产欧美日韩一区二区三区在线| 19禁男女啪啪无遮挡网站| 麻豆国产av国片精品| 亚洲色图综合在线观看| 999精品在线视频| 亚洲熟妇中文字幕五十中出| 大香蕉久久成人网| 欧美久久黑人一区二区| 女同久久另类99精品国产91| av视频免费观看在线观看| www.www免费av| 久久人人爽av亚洲精品天堂| 激情视频va一区二区三区| 欧美黄色片欧美黄色片| 99riav亚洲国产免费| 人人澡人人妻人| 国产精品综合久久久久久久免费 | 亚洲精品久久国产高清桃花| 午夜久久久久精精品| 在线观看www视频免费| 男女床上黄色一级片免费看| 午夜精品久久久久久毛片777| 久久 成人 亚洲| 欧美+亚洲+日韩+国产| 国产成人精品在线电影| av片东京热男人的天堂| 欧美性长视频在线观看| 日本 av在线| 国产精品 国内视频| 国产aⅴ精品一区二区三区波| 国产精品一区二区精品视频观看| 女人被躁到高潮嗷嗷叫费观| 99在线人妻在线中文字幕| 国产高清有码在线观看视频 | 热re99久久国产66热| 神马国产精品三级电影在线观看 | 国产成+人综合+亚洲专区| 午夜福利18| 亚洲精品中文字幕一二三四区| 欧美中文综合在线视频| 一级黄色大片毛片| 久久精品影院6| 老司机午夜十八禁免费视频| 欧美性长视频在线观看| 少妇被粗大的猛进出69影院| 色精品久久人妻99蜜桃| 两个人看的免费小视频| 午夜视频精品福利| www.999成人在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 热re99久久国产66热| 亚洲欧美日韩无卡精品| 大型av网站在线播放| 国产精品久久电影中文字幕| 亚洲人成网站在线播放欧美日韩| 国产精品二区激情视频| 老司机午夜十八禁免费视频| 亚洲最大成人中文| 午夜免费观看网址| 在线观看日韩欧美| 亚洲一区二区三区不卡视频| 日韩精品免费视频一区二区三区| 亚洲七黄色美女视频| av天堂久久9| 亚洲精品中文字幕在线视频| 国产精品爽爽va在线观看网站 | 香蕉丝袜av| 国产成人啪精品午夜网站| 欧美成人午夜精品| 国产一区二区三区综合在线观看| 99精品久久久久人妻精品| 黑人操中国人逼视频| 久久国产精品男人的天堂亚洲| 国产精品久久视频播放| 精品国产乱子伦一区二区三区| 韩国精品一区二区三区| 一进一出抽搐gif免费好疼| 久久香蕉国产精品| 夜夜夜夜夜久久久久| 国语自产精品视频在线第100页| 亚洲成av片中文字幕在线观看| 非洲黑人性xxxx精品又粗又长| 欧美 亚洲 国产 日韩一| 精品国产美女av久久久久小说| 久久久国产成人免费| 国产免费男女视频| 在线永久观看黄色视频| 亚洲五月婷婷丁香| 国产亚洲精品综合一区在线观看 | 成年女人毛片免费观看观看9| 亚洲欧美精品综合一区二区三区| 中文字幕人妻丝袜一区二区| 久久久精品国产亚洲av高清涩受| 亚洲va日本ⅴa欧美va伊人久久| 亚洲成人国产一区在线观看| 国产精品秋霞免费鲁丝片| 亚洲最大成人中文| 长腿黑丝高跟| 成年女人毛片免费观看观看9| 在线观看免费视频日本深夜| 久久香蕉国产精品| 大码成人一级视频| 亚洲欧美激情在线| 国产精品自产拍在线观看55亚洲| 美女扒开内裤让男人捅视频| 亚洲精品在线美女| 99国产综合亚洲精品| 日本三级黄在线观看| 国产精品综合久久久久久久免费 | 日韩一卡2卡3卡4卡2021年| 国产精品国产高清国产av| www.999成人在线观看| www国产在线视频色| 妹子高潮喷水视频| 久久精品人人爽人人爽视色| 亚洲男人的天堂狠狠| 国产高清视频在线播放一区| 国产成人精品久久二区二区91| 九色国产91popny在线| 乱人伦中国视频| 亚洲精华国产精华精| 成人欧美大片| 国产主播在线观看一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩国内少妇激情av| 日韩欧美三级三区| 性欧美人与动物交配| 精品国产超薄肉色丝袜足j| 国内毛片毛片毛片毛片毛片| 欧美久久黑人一区二区| 91成人精品电影| 大香蕉久久成人网| 国产免费男女视频| 88av欧美| 中出人妻视频一区二区| 少妇 在线观看| 久久久水蜜桃国产精品网| 亚洲av电影不卡..在线观看| 欧美日韩黄片免| 亚洲aⅴ乱码一区二区在线播放 | 欧美老熟妇乱子伦牲交| 桃色一区二区三区在线观看| 激情视频va一区二区三区| 叶爱在线成人免费视频播放| 欧美久久黑人一区二区| 亚洲国产精品999在线| 男人操女人黄网站| 亚洲av五月六月丁香网| 亚洲精品粉嫩美女一区| 午夜成年电影在线免费观看| 日韩欧美国产一区二区入口| 久久久久九九精品影院| 国产亚洲精品一区二区www| 91麻豆精品激情在线观看国产| 日韩一卡2卡3卡4卡2021年| 一区二区三区激情视频| 叶爱在线成人免费视频播放| 精品国产乱子伦一区二区三区| 日韩 欧美 亚洲 中文字幕| 黄色片一级片一级黄色片| 午夜激情av网站| av在线天堂中文字幕| 国产成人av激情在线播放| 99re在线观看精品视频| 黄色片一级片一级黄色片| 国语自产精品视频在线第100页| 国产一区二区三区在线臀色熟女| 国产日韩一区二区三区精品不卡| 人人妻人人澡欧美一区二区 | 欧美精品亚洲一区二区| 色综合站精品国产| 日韩欧美国产在线观看| 丁香欧美五月| 日本一区二区免费在线视频| 亚洲中文字幕日韩| 免费女性裸体啪啪无遮挡网站| 日日摸夜夜添夜夜添小说| 日韩欧美一区二区三区在线观看| 少妇的丰满在线观看| 自线自在国产av| 国产亚洲av嫩草精品影院| 国产乱人伦免费视频| 黄色视频,在线免费观看| 两个人看的免费小视频| 日本欧美视频一区| 性少妇av在线| 国产欧美日韩一区二区三区在线| 一边摸一边抽搐一进一小说| 在线十欧美十亚洲十日本专区| 国产一区二区激情短视频| 久久影院123| 精品一品国产午夜福利视频| 很黄的视频免费| 国产精品综合久久久久久久免费 | 欧美日韩亚洲国产一区二区在线观看| 神马国产精品三级电影在线观看 | 身体一侧抽搐| 久久精品91无色码中文字幕| 99国产精品一区二区三区| 淫秽高清视频在线观看| 精品久久久久久久久久免费视频| 欧美黑人精品巨大| 极品人妻少妇av视频| 日韩免费av在线播放| 久久性视频一级片| 亚洲色图av天堂| 欧美成人性av电影在线观看| 精品卡一卡二卡四卡免费| 日本精品一区二区三区蜜桃| 日本三级黄在线观看| 成人三级黄色视频| 亚洲人成伊人成综合网2020| 亚洲一区二区三区色噜噜| 亚洲五月婷婷丁香| 国产精品1区2区在线观看.| av在线天堂中文字幕| 亚洲午夜理论影院| 19禁男女啪啪无遮挡网站| 欧美大码av| 欧美黑人精品巨大| 91麻豆av在线| 日本五十路高清| 国产精品国产高清国产av| 久久中文字幕人妻熟女| 精品国内亚洲2022精品成人| 国产私拍福利视频在线观看| 亚洲国产中文字幕在线视频| 男人舔女人下体高潮全视频| 性欧美人与动物交配| 涩涩av久久男人的天堂| 日本撒尿小便嘘嘘汇集6| 一进一出抽搐gif免费好疼| 十八禁人妻一区二区| 国产精品一区二区精品视频观看| 狠狠狠狠99中文字幕| 一本综合久久免费| 亚洲无线在线观看| 一区二区三区国产精品乱码| 亚洲精品粉嫩美女一区| 99精品在免费线老司机午夜| АⅤ资源中文在线天堂| 欧美一级a爱片免费观看看 | 亚洲第一电影网av| 欧美黑人精品巨大| 欧美中文综合在线视频| 日韩 欧美 亚洲 中文字幕| 精品久久蜜臀av无| 久久久国产欧美日韩av| 国产国语露脸激情在线看| 国产精品综合久久久久久久免费 | 精品国产乱子伦一区二区三区| 国产精品二区激情视频| 欧美久久黑人一区二区| 18禁国产床啪视频网站| 久久伊人香网站| av视频免费观看在线观看| www.精华液| 日本 av在线| 国产精品 欧美亚洲| 在线永久观看黄色视频| 亚洲欧洲精品一区二区精品久久久| 男女做爰动态图高潮gif福利片 | 日韩国内少妇激情av| 亚洲欧美一区二区三区黑人| 中文字幕精品免费在线观看视频| 欧美激情高清一区二区三区| 久久精品国产综合久久久| 大陆偷拍与自拍| 熟妇人妻久久中文字幕3abv| 一夜夜www| 中文字幕精品免费在线观看视频| 叶爱在线成人免费视频播放| 亚洲午夜理论影院| 国产av精品麻豆| 熟妇人妻久久中文字幕3abv| 狠狠狠狠99中文字幕| 手机成人av网站| 电影成人av| 亚洲中文字幕日韩| tocl精华| 精品电影一区二区在线| 亚洲人成伊人成综合网2020| 国产三级在线视频| 精品欧美一区二区三区在线| 在线观看一区二区三区| www.熟女人妻精品国产| 1024香蕉在线观看| 欧美大码av| 深夜精品福利| 国产1区2区3区精品| 亚洲狠狠婷婷综合久久图片| 成人精品一区二区免费| 久久精品人人爽人人爽视色| 亚洲国产日韩欧美精品在线观看 | 香蕉久久夜色| 两个人视频免费观看高清| 国产精品99久久99久久久不卡| 1024视频免费在线观看| 亚洲精品中文字幕在线视频| 国产欧美日韩综合在线一区二区| www.自偷自拍.com| 狂野欧美激情性xxxx|