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

    高負荷擴壓葉柵吹風試驗流場二維性控制技術(shù)研究

    2021-05-18 02:28:06高麗敏黎浩學
    實驗流體力學 2021年2期
    關(guān)鍵詞:附面層葉柵馬赫數(shù)

    高麗敏, 劉 哲, 蔡 明, 黎浩學

    西北工業(yè)大學 動力與能源學院, 翼型、葉柵空氣動力學國家級重點實驗室, 西安 710129

    0 引 言

    現(xiàn)代航空發(fā)動機推重比的不斷提高要求壓氣機在減少級數(shù)的同時實現(xiàn)高壓比、高效率,設(shè)計使用高負荷擴壓葉片是最直接的解決方法,而平面葉柵吹風試驗是先進葉型設(shè)計不可或缺的試驗環(huán)節(jié)[1-6]。

    平面葉柵試驗中葉柵流場需沿展向滿足一定的二維性[1, 7-10],才能獲得較為可靠的二維葉型性能數(shù)據(jù)。然而,隨著擴壓葉柵負荷的不斷提升,流道內(nèi)逆壓梯度與附面層干涉會加劇附面層累積,導(dǎo)致流道收縮,測量截面(50%葉高)出口的軸向速度密流增加,從而使得流場二維性惡化,造成測量結(jié)果失真[7-13]。一般試驗中通過增大展弦比來保證測量截面的二維性[1,8,13],但流道收縮程度隨負荷升高而加劇,所需展弦比更大,受限于風洞尺寸和成本,此方法不具備普適性,且變工況性能較差。目前試驗中多采用端壁抽吸的主動控制方法保證流場二維性[8-10,12-19]。

    美國NACA于20世紀20年代中期發(fā)展的“多孔壁”(Porous wall)技術(shù)[14]將多孔材料作為葉柵端壁以吸除端壁附面層,Briggs[15]、Herrig等[16]通過平面葉柵試驗證明了此技術(shù)的可靠性。1967年P(guān)ollard等[13]發(fā)現(xiàn)“多孔壁”技術(shù)雖然能獲得較好的二維結(jié)果,但無法完全吸除角區(qū)低能流體。1975年德國Starken等[17]在大負荷擴壓葉柵吹風預(yù)試驗中,發(fā)現(xiàn)不同位置抽吸槽對氣源能力的利用率以及對來流的干擾不同。國內(nèi)的相關(guān)研究起步較晚,1995年中國燃氣渦輪研究院姜正禮[18]在試驗中采取尾部開槽的方式改善二維性,研究了二維性對葉柵性能的影響規(guī)律。2014年西北工業(yè)大學鄧熙等[7]研究了出氣角的二維修正方法,但發(fā)現(xiàn)全弦長柵狀槽僅可實現(xiàn)局部控制。2020年沈陽發(fā)動機研究所王東等[19]通過試驗研究發(fā)現(xiàn)流道端壁抽吸對減小密流比起決定作用,柵前抽吸作用不大。綜上所述,抽吸對改善葉柵流場二維性的作用顯而易見。但目前抽吸控制葉柵流場二維性的研究大多聚焦于控制效果,對抽吸控制方法研究較少,如抽吸位置的選擇缺乏理論依據(jù)、缺乏新的抽吸方法等等。

    近年來,對吸附式葉柵的研究[20-23]表明,不同葉柵的抽吸布局、最佳抽吸位置及吸氣量并不相同,組合抽吸可以進一步優(yōu)化控制效果。因此,本文借鑒吸附式葉柵的研究成果,針對高負荷擴壓葉柵,采用數(shù)值方法分析影響葉柵流場二維性的因素,探索不同端壁抽吸位置對流場二維性的控制機理,驗證分布式抽吸方法對流場二維性的控制效果,旨在為高負荷擴壓葉柵吹風試驗流場二維性的控制技術(shù)優(yōu)化提供參考。

    1 研究對象及研究方案

    1.1 研究對象

    本文研究的二維平面葉柵幾何定義如圖1所示,主要幾何尺寸參數(shù)見表1。其擴散因子為0.63,稠度為2.77,彎角為54.12°,為典型的高負荷擴壓葉柵。

    圖1 平面葉柵幾何定義Fig.1 Linear cascade configuration

    表1 平面葉柵幾何參數(shù)Table 1 Linear cascade geometric parameter

    按照《航空發(fā)動機設(shè)計手冊》中的葉柵試驗件標準,參考西北工業(yè)大學高亞聲速平面葉柵風洞試驗段尺寸,將二維葉柵等比模化,常規(guī)葉柵試驗件的計算模型展向高度h=100 mm,弦長b=65 mm。

    依據(jù)本文研究目的,計算模型除常規(guī)試驗件外,另在常規(guī)試驗件基礎(chǔ)上,于柵板開設(shè)抽吸槽,進行端壁抽吸以控制吹風試驗中葉柵流場二維性。葉柵流場二維性的變化與端壁附面層發(fā)展直接相關(guān)。端壁附面層沿軸向發(fā)展經(jīng)歷了附面層累積階段、分離起始階段、充分發(fā)展階段,即端壁附面層形態(tài)沿著軸向位置而不同。因此,本文提出了3種軸向位置的柵板抽吸方案:Front:(0~33%)b,Middle:(33%~66%)b,Tail:(66%~99%)b,控制單一變量為軸向位置,使抽吸槽沿周向覆蓋整個端壁。參考現(xiàn)有研究的建模思路,為了兼顧柵板強度和抽吸能力[22-23],葉柵端壁開設(shè)的抽吸槽寬度為1.5 mm,周向相鄰抽吸槽間距2 mm,抽吸槽布置如圖2所示。

    圖2 抽吸控制葉柵開槽方案Fig.2 Slot configuration scheme of controlled cascade

    由于研究內(nèi)容是葉柵吹風試驗的流場二維性,因此忽略周期性的影響,分別對二維平面葉柵、常規(guī)試驗葉柵和抽吸控制葉柵進行單通道計算,其中抽吸控制葉柵包括3種不同軸向位置的柵板抽吸方案,共計5種方案。為便于表示,將二維平面葉柵命名為Ideal,常規(guī)試驗葉柵命名為No suc, 前部、中部、尾部抽吸控制葉柵分別命名為Front、Middle、Tail。計算工況為進口馬赫數(shù)Ma1=0.40,攻角i=-8°~8°。

    1.2 數(shù)值方法及有效性驗證

    本文5種幾何模型的主流道網(wǎng)格相同,均為常規(guī)試驗件的流道。葉柵主流道網(wǎng)格劃分使用軟件NUMECA中AutoGrid5模塊生成O4H型結(jié)構(gòu)化網(wǎng)格拓撲,對于抽吸槽網(wǎng)格則使用IGG模塊手動生成H型結(jié)構(gòu)化網(wǎng)格。抽吸槽網(wǎng)格與葉柵通道網(wǎng)格之間采用完全非匹配連接方式實現(xiàn)界面間的數(shù)值傳遞,葉柵主流道網(wǎng)格及抽吸槽局部放大網(wǎng)格如圖3所示。

    圖3 葉柵計算域網(wǎng)格及抽吸槽局部放大網(wǎng)格Fig.3 Computational grid of cascade and local enlargement of slot grid

    對主流道網(wǎng)格進行網(wǎng)格無關(guān)性驗證,如圖4所示,當網(wǎng)格量達到700 000時,總壓損失系數(shù)ω幾乎不隨網(wǎng)格量增加而增加,因此主流道網(wǎng)格量選擇700 000。抽吸控制葉柵的抽吸槽網(wǎng)格與相鄰主流道端壁網(wǎng)格尺度保持一致,同時對壁面網(wǎng)格進行加密,保證壁面y+值小于1。

    圖4 網(wǎng)格無關(guān)性驗證Fig.4 Grid independence verification

    數(shù)值方法采用Fine/Turbo模塊求解定常雷諾平均N-S方程,湍流模型選用SA模型??臻g離散為中心差分格式,時間項采用4階Runge-Kutta方法迭代求解,CFL數(shù)為3。進口邊界條件給定總參數(shù)及進氣方向,出口給定背壓,葉片表面給定絕熱無滑移壁面,二維平面葉柵端壁給定mirror模擬無限長葉片,常規(guī)葉柵試驗件端壁、抽吸控制葉柵端壁及其抽吸槽壁面給定絕熱無滑移壁面,抽吸控制葉柵抽吸槽出口給定流量出口。

    因本文研究對象的葉柵試驗數(shù)據(jù)缺乏,為驗證本文數(shù)值方法的正確性,故選取試驗數(shù)據(jù)庫完整、可信度高的NPU-A1葉型模化后的三維常規(guī)試驗葉柵進行對比工作[7,9,22]。工況為i=0°、Ma1=0.55。

    試驗前期準備:1) 檢測空風洞流場。其核心區(qū)分別占展向、周向范圍的90%、80%,馬赫數(shù)0.55±0.005,展向、周向氣流偏角±0.5°,空風洞流場均勻。2) 檢測葉柵流場。葉柵中部2個葉柵通道滿足Ma1=0.55±0.005,來流均勻,尾跡深度、寬度一致性較高,滿足試驗要求。

    圖5為測量截面的葉片表面等熵馬赫數(shù)分布對比[22],試驗時由于葉片厚度的原因,近前緣和尾緣處無法開設(shè)靜壓孔,因此有部分數(shù)據(jù)缺失,但兩者等熵馬赫數(shù)峰值點位置相同,分布趨勢一致,數(shù)值結(jié)果能夠較好地對應(yīng)試驗結(jié)果。

    圖5 表面等熵馬赫數(shù)分布[22]Fig.5 Isentropic Mach number distribution[22]

    圖6為總壓損失的展向分布對比??梢钥闯觯瑪?shù)值方法結(jié)果與試驗結(jié)果的總壓損失沿展向具有一致的分布及變化趨勢,數(shù)值雖然存在差異,但處于可接受范圍,可以認為本文的數(shù)值方法正確可靠。

    圖6 總壓損失系數(shù)展向分布[22]Fig.6 Span distribution of ω[22]

    1.3 數(shù)據(jù)處理方法

    平面葉柵試驗測量截面選取50%葉高處,進口參數(shù)通過距離前緣1倍弦長處的上游測量站測量均勻來流得到,出口參數(shù)通過尾緣后1倍弦長處的下游測量站(保證主流與尾跡充分摻混)測量得到。通過質(zhì)量平均的方法處理1個柵距內(nèi)的流動參數(shù),得到表征葉柵性能的平均參數(shù),如式(1)所示。

    (1)

    依據(jù)試驗情況,使用離散格式進行處理:

    (2)

    其中,C為平均參數(shù),如壓力p、總壓損失系數(shù)ω等。t為葉柵標號,dy為相鄰離散點間隔,Ci、ρi、vai分別為下游測量站的離散點處的物理參數(shù)值、密度和軸向速度,離散點為n個,Δl為相鄰離散點間隔。

    2 端壁附面層對流場二維性影響

    本節(jié)對比二維平面葉柵與常規(guī)葉柵流場,旨在分析常規(guī)試驗葉柵吹風試驗中流場二維性的影響因素。

    由于葉柵吹風試驗的流場二維性受端壁附面層發(fā)展造成的流道收縮影響,因此進出口流道收縮比可以作為流場二維性的衡量標準。進出口流道收縮比可用出口與進口的軸向速度密流比(Axial Velocity Density Ratio,AVDR)表示:

    (3)

    其中,A、β、ρ、va分別為流道面積、氣流角、密度、軸向速度,下標1、2分別表示進口、出口。

    圖7為攻角i=0°、進口馬赫數(shù)Ma1=0.40時二維平面葉柵與常規(guī)葉柵50%葉高截面的馬赫數(shù)云圖及流線分布對比。由圖7可知,受擴壓葉柵內(nèi)逆壓梯度的影響,二維平面葉柵通道中部流動發(fā)生分離,該分離區(qū)向下游發(fā)展過程中并未附著,最終形成開式分離區(qū)。與之相比,常規(guī)試驗葉柵主流流線并未分離,持續(xù)附著,近尾緣處的馬赫數(shù)比二維平面葉柵高0.07,擴壓能力下降,流場已經(jīng)失真。

    圖7 馬赫數(shù)云圖及流線Fig.7 Mach number contour and streamline

    圖8的葉片體流線表明,由于逆壓梯度的存在,常規(guī)試驗葉柵端壁附面層在近前緣處即發(fā)生了分離,進而產(chǎn)生了大范圍的回流區(qū),向中部擠壓主流,馬赫數(shù)云圖沿軸向出現(xiàn)“凸起”,說明主流在擠壓下加速,抗逆壓能力增強,導(dǎo)致常規(guī)葉柵50%葉高截面流場與二維平面葉柵不同,造成流場失真。

    圖8 體流線Fig.8 Volume line

    常規(guī)試驗葉柵的軸向速度密流比的攻角特性曲線及常規(guī)試驗葉柵與二維平面葉柵總壓損失的攻角特性曲線如圖9和10所示。圖9表明,常規(guī)試驗葉柵軸向速度密流比均大于1.0,隨攻角線性增加,葉柵流場二維性逐漸惡化。結(jié)合圖10,軸向速度密流比增大會使失真程度增加,總壓損失偏差增大,其中0°攻角下?lián)p失與真實情況相比減小50%左右,整體最小偏差達23%,偏差過大。由此可見,葉柵試驗結(jié)果的正確性依賴于葉柵流場二維性,調(diào)控手段是必要的。

    圖9 軸向速度密流比攻角特性Fig.9 AVDR of incidence characteristics

    圖10 總壓損失攻角特性Fig.10 ω of incidence characteristics

    3 控制結(jié)果分析

    改善平面葉柵流場二維性的最佳抽吸結(jié)構(gòu)設(shè)計應(yīng)是:1) 抽吸結(jié)構(gòu)對來流擾動小。2) 抽吸功耗小,最大化利用真空泵能力。3) 改善效果明顯,包括葉柵性能參數(shù)、葉片負荷、展向二維區(qū)域?qū)挾鹊雀纳啤?/p>

    實際平面葉柵試驗中流場達到準二維狀態(tài)即可,相關(guān)學者建議[8],下游測量站的軸向速度密流比臨界值為1.1,因此本文在此條件下對比不同抽吸位置的優(yōu)劣。常規(guī)葉柵在攻角i=0°、進口馬赫數(shù)Ma1=0.40下軸向速度密流比可達1.4,改善效果對比明顯,選擇其為對比工況。

    3.1 抽吸對來流影響及抽吸功耗對比

    對來流的擾動主要表現(xiàn)為對來流方向的擾動,圖11為不同抽吸位置下進氣角β1隨軸向速度密流比的變化曲線。

    圖11 進氣角隨軸向速度密流比變化Fig.11 Variation of β1 with AVDR

    整體來看,抽吸槽距進口越遠,對進氣角的擾動越小。相比之下,中部抽吸在軸向速度密流比為1.1時進氣角最接近理想二維狀態(tài),前部抽吸對進氣角擾動最大,尾部抽吸時的進氣角則始終維持在常規(guī)試驗葉柵的進氣角附近。整體上,3種位置抽吸的進氣角與設(shè)計進氣角偏差不超過0.25°,滿足試驗要求。

    抽吸功耗則由相對抽吸流量qmsi(抽吸流量占主流流量的百分比)、抽吸槽入口相對主流進口的靜壓系數(shù)cp表征:

    (4)

    其中,p為靜壓,下標1、si分別表示葉柵進口、抽吸槽進口。

    對比下游測量站軸向速度密流比為1.1時的抽吸流量百分比qmsi及抽吸槽入口靜壓系數(shù)cp,即可對比3種位置所需的抽吸功耗。

    表2 抽吸功耗對比Table 2 Suction power consumption

    前部抽吸的抽吸流量和靜壓系數(shù)均最小,相比之下,中部抽吸時分別增加67%、15%,尾部抽吸時則分別增加了3倍和13倍。這是因為擴壓葉柵內(nèi)沿軸向端壁附面層會持續(xù)累積并在角區(qū)位置形成大范圍低能區(qū),所需要的抽吸流量隨之增加;越靠近葉柵尾部,靜壓越高。兩者共同作用導(dǎo)致尾部抽吸功耗最大,因此在真空泵抽吸能力一定時,應(yīng)避免在尾部抽吸。

    3.2 改善效果及機理分析

    以二維平面葉柵為基準,通過對比不同抽吸位置的葉柵性能參數(shù)及其偏差衡量改善效果。分別給出不同抽吸位置控制下的靜壓升p2/p1和總壓損失系數(shù)ω隨軸向速度密流比的變化曲線,同時給出常規(guī)試驗葉柵和二維平面葉柵的相關(guān)性能參數(shù),如圖12和13所示。

    圖12 靜壓升隨軸向速度密流比變化Fig.12 Variation of p2/p1 with AVDR

    圖13 總壓損失系數(shù)隨軸向速度密流比變化Fig.13 Variation of ω with AVDR

    相對于二維平面葉柵,常規(guī)試驗葉柵由于流道收縮、主流加速、擴壓性下降,靜壓升下降了0.02,占比2%。主流加速造成開式分離消失,會減小葉型損失,總壓損失系數(shù)下降了0.04,占比達50%。對于控制葉柵,隨軸向速度密流比下降,葉柵擴壓性明顯升高,損失也隨之增加,逐漸向理想二維狀態(tài)靠攏。軸向速度密流比達到1.1時,前部、中部抽吸的靜壓升與二維平面葉柵一致,尾部抽吸則表現(xiàn)出了過改善;但3種抽吸位置對損失改善均有不足,前部、中部、尾部抽吸與二維平面葉柵的偏差分別為12.5%、7.5%和25%。相比之下,尾部抽吸對靜壓升p2/p1改善過度,對總壓損失系數(shù)ω則表現(xiàn)出改善不足。

    葉片表面等熵馬赫數(shù)分布表征擴壓葉柵的擴壓程度以及表面附面層流動狀況,圖14和15分別為葉片表面等熵馬赫數(shù)Mais分布、與二維平面葉柵的等熵馬赫數(shù)偏差ΔMais分布。

    圖14 等熵馬赫數(shù)Mais對比Fig.14 Comparison of isentropic Mach number Mais

    常規(guī)試驗葉柵受流道收縮的影響,對主流有加速效應(yīng),因此等熵馬赫數(shù)整體最大,其偏差ΔMais整體超過0.05,尾緣等熵馬赫數(shù)升高會導(dǎo)致擴壓性下降;抽吸控制葉柵的表面等熵馬赫數(shù)分布則均有不同程度的改善。下文結(jié)合不同軸向位置的軸向速度密流比變化進行分析。

    由于葉柵通道狹窄,試驗中無法測得通道內(nèi)部的軸向速度密流比,因此現(xiàn)有研究中軸向速度密流比數(shù)據(jù)大都基于下游測量站數(shù)據(jù)[13],但其僅能反映端壁附面層對流場二維性的影響在下游測量站的總和,不利于精細化調(diào)控,因此,有必要了解流道內(nèi)軸向速度密流比變化。由于葉片的厚度分布及葉片表面附面層的發(fā)展都會引起實際流道收縮,因此,為單獨研究端壁附面層發(fā)展對流道收縮的影響,采用與二維平面葉柵流道內(nèi)相同軸向位置的軸向速度密流比差值(ΔAVDR)代表端壁附面層的作用,如圖16所示。

    圖16 ΔAVDR分布對比Fig.16 Comparison of ΔAVDR Distribution

    常規(guī)試驗葉柵在進口前端壁附面層逐漸累積,ΔAVDR增長較緩。進入葉柵流道內(nèi),受逆壓梯度影響,端壁附面層累積加劇,發(fā)生分離,流道收縮劇烈,ΔAVDR劇烈增大。0.04 m后ΔAVDR略微下降,結(jié)合第2節(jié)的分析可知,這是開式分離消失與端壁附面層發(fā)展的耦合作用所致。開式分離消失使通流面積增加,端壁附面層發(fā)展使通流面積減小。開式分離消失本質(zhì)上是附面層發(fā)展對主流的加速作用所致,實際上流場已經(jīng)失真。葉柵出口端壁附面層向下游發(fā)展與主流摻混會引起ΔAVDR繼續(xù)增大,較遠下游處,其影響逐漸減弱,最終趨于平緩。

    前部、中部抽吸分別作用于前部附面層、附面層分離起始區(qū),能夠在改善當?shù)囟S性的同時對下游產(chǎn)生正效果,使ΔAVDR整體維持在較低的水平,Mais分布整體接近理想二維值。相比之下,前部抽吸會對流道前部產(chǎn)生過抽吸,使當?shù)氐褥伛R赫數(shù)低于理想二維值,但其偏差ΔMais不大(見圖15),中部抽吸則對全流道端壁附面層的控制效果最佳。尾部抽吸直接作用于附面層完全分離區(qū)域,ΔAVDR于葉片流道中部開始驟降,僅局部改善葉片負荷,前部并無改善甚至惡化,偏差ΔMais最大達0.08。

    圖15 等熵馬赫數(shù)偏差ΔMais分布對比Fig.15 Comparison of isentropic Mach number deviation ΔMais

    為直觀上理解不同位置抽吸控制的作用機理,觀察通道內(nèi)端壁附面層的軸向發(fā)展。圖17為S3流面的總壓損失系數(shù)分布,在葉柵通道內(nèi)垂直于軸向截取間距相等的7個截面,自前緣向尾緣分別標記為面1~7。

    圖17 端壁附面層發(fā)展Fig.17 End wall boundary layer development

    常規(guī)試驗葉柵自面1至面3附面層累積,面3之后端壁附面層在逆壓梯度作用下開始發(fā)生分離并逐漸發(fā)展,在面7近端壁處形成大范圍低能區(qū)域,主流流道明顯收縮。端壁附面層經(jīng)過前部抽吸始終保持較薄的狀態(tài),未見明顯高損失區(qū),角區(qū)低能區(qū)強度整體減弱,周向影響范圍明顯減小,因此對全葉展均有改善效果,但其展向影響范圍擴大,會影響展向二維區(qū)域?qū)挾?。中部抽吸與前部抽吸相比,面3之前的附面層有所累積,但經(jīng)過抽吸后,對面3之后的端壁附面層仍表現(xiàn)出較強的控制能力,出口處低能流體的范圍及損失強度同樣有明顯改善。尾部抽吸時,常規(guī)試驗葉柵面1~5的附面層并未削減,僅在近尾緣處損失強度驟減,低能區(qū)域向兩側(cè)移動,因此其對葉片負荷僅有局部改善。

    平面葉柵試驗中,葉柵的中部應(yīng)保持一定的展向二維區(qū),以避免端壁低能區(qū)的影響,確保下游測量站所測參數(shù)的真實性。圖18為不同抽吸位置下出口展向軸向速度密流比分布。

    圖18表明,常規(guī)試驗葉柵由端壁向中心軸向速度密流比急劇升高,50%葉高處軸向速度密流比最高,可達1.4,說明兩側(cè)存在大范圍低能區(qū),中部主流加速明顯。抽吸控制葉柵50%葉高處軸向速度密流比均為1.1。相比之下,前部、中部抽吸的全葉高二維性整體獲得改善,但會對下游造成擾動,使出口展向二維區(qū)偏窄。尾部抽吸的展向二維區(qū)域達展向高度60%左右,低速區(qū)完全分布在兩側(cè)端壁。

    圖18 出口展向軸向速度密流比分布Fig.18 Span distribution of AVDR

    試驗中,希望在全工況下都能得到二維數(shù)據(jù),因此不同抽吸位置的控制能力也是我們所關(guān)心的。圖19為出口軸向速度密流比為1.1時不同抽吸槽的相對抽吸流量(qmsi)的攻角特性線。

    圖19 不同抽吸位置流量攻角特性Fig.19 Suction flow of incidence at different positions

    相比之下,前部抽吸主要作用于附面層,流量隨攻角變化很小,控制效果最佳,全攻角特性很好。尾部抽吸流量始終最大,隨攻角線性增加, 其作用于附面層完全分離區(qū),對出口軸向速度密流比值有著直接影響,變化趨勢與軸向速度密流比的發(fā)展趨勢一致(對比圖9)。中部抽吸的負攻角特性較好,當攻角i>0°控制效果明顯下降,隨攻角增加,附面層分離提前,中部位置逐漸成為附面層完全分離區(qū),因此正攻角特性與尾部抽吸類似。

    3.3 分布式抽吸探索

    前部、中部抽吸對葉柵二維性的整體控制更好,但出口展向二維區(qū)偏??;尾部抽吸的展向二維區(qū)更寬,但僅能局部改善。為實現(xiàn)較好的控制效果且保證出口展向二維區(qū)寬度,探索了分布式(Distributed)抽吸方法。

    對前部、中部和尾部同時抽吸并分別控制流量,在i=0°、Ma1=0.40時進行控制,軸向速度密流比控制在1.1左右。圖20為分布式抽吸控制后50%葉高截面云圖與二維流場對比(含流線)。

    圖20 馬赫數(shù)云圖及流線Fig.20 Mach number contour and streamline

    可以看出,分布式抽吸后,50%葉高截面出現(xiàn)二維平面葉柵流場的開式分離特征。相比常規(guī)試驗葉柵,50%葉高截面近尾緣處與理想二維狀態(tài)的馬赫數(shù)偏差由0.07降為0.01(對比圖7),擴壓性上升,高度還原了二維流場。

    分布式抽吸下出口展向二維區(qū)域分布如圖21所示。圖21表明,分布式抽吸能夠明顯改善前部、中部抽吸出口二維區(qū)偏窄的問題。圖22為通道內(nèi)端壁附面層的軸向變化,端壁附面層自面3分離后強度持續(xù)減弱,并無單一尾部抽吸的低能區(qū)域驟減的現(xiàn)象,因此對通道整體流場均有改善。同時在尾部抽吸作用下,角區(qū)低能區(qū)遠離中部區(qū)域,因此流道中部二維區(qū)域明顯加寬。

    圖21 出口展向軸向速度密流比分布Fig.21 Span distribution of AVDR

    圖22 端壁附面層發(fā)展Fig.22 End wall boundary layer development

    因此,針對本文高負荷擴壓葉柵,分布式抽吸能夠擴展出口二維區(qū)域并實現(xiàn)流場二維性的整體控制,值得探索與應(yīng)用。

    4 結(jié) 論

    本文研究了高負荷擴壓葉柵吹風試驗流場二維性的影響因素及不同軸向位置端壁抽吸的主動控制技術(shù),得到以下結(jié)論:

    1) 常規(guī)試驗葉柵端壁附面層發(fā)生了分離,低能區(qū)擠壓主流并使其加速,相比二維平面葉柵擴壓性下降,造成流場失真。軸向速度密流比隨攻角線性增加,損失與理想二維狀態(tài)相比最小偏差達23%??刂剖侄问潜匾摹?/p>

    2) 對比3種位置抽吸控制方法,抽吸槽離進口越遠,對來流影響越小,但整體上偏差控制在±0.25°,滿足試驗要求。前部抽吸的抽吸流量和抽吸槽靜壓系數(shù)均最小,相比之下,中部抽吸時分別增加67%、15%,尾部抽吸時則分別增加了3倍和13倍。

    3) 前部、中部抽吸可整體控制葉柵二維性,但出口展向二維區(qū)偏窄。尾部抽吸出口展向二維區(qū)較寬,但僅局部改善近尾緣處二維性。前部抽吸在全攻角下控制良好,中部抽吸的負攻角特性較好,尾部抽吸流量則隨攻角呈線性變化。

    4) 分布式抽吸能夠擴展出口二維區(qū)域并實現(xiàn)流場二維性的整體控制,值得探索與應(yīng)用。

    展望:本文針對附面層沿軸向的發(fā)展階段進行抽吸,未考慮附面層沿周向的發(fā)展,后期將針對周向不同位置抽吸進行研究。

    猜你喜歡
    附面層葉柵馬赫數(shù)
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    基于網(wǎng)格框架的非結(jié)構(gòu)附面層網(wǎng)格生成技術(shù)
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    亞聲速壓氣機平面葉柵及其改型的吹風試驗
    基于數(shù)值模擬的流場附面層邊緣識別方法
    載荷分布對可控擴散葉型性能的影響
    超聲壓氣機葉柵流場的數(shù)值模擬與試驗驗證
    超音速透平葉柵造型設(shè)計及分析
    針對軸流壓氣機的非軸對稱端壁造型優(yōu)化設(shè)計
    發(fā)動機進口附面層測量試驗與數(shù)值模擬
    国产成人精品无人区| 亚洲欧美成人精品一区二区| 精品一区二区三卡| 精品亚洲成国产av| 免费黄频网站在线观看国产| 久久国产精品男人的天堂亚洲| 精品一区二区三卡| 国产一卡二卡三卡精品 | 国产亚洲午夜精品一区二区久久| 免费在线观看黄色视频的| 男人操女人黄网站| 国产精品无大码| 别揉我奶头~嗯~啊~动态视频 | 高清视频免费观看一区二区| 国产免费现黄频在线看| 在线天堂最新版资源| 老鸭窝网址在线观看| 只有这里有精品99| 99热国产这里只有精品6| 日本vs欧美在线观看视频| 啦啦啦在线观看免费高清www| 2018国产大陆天天弄谢| 国产精品欧美亚洲77777| 欧美精品av麻豆av| 午夜免费鲁丝| av卡一久久| 又黄又粗又硬又大视频| 国产老妇伦熟女老妇高清| 在线免费观看不下载黄p国产| 美女视频免费永久观看网站| 十八禁高潮呻吟视频| 日韩 欧美 亚洲 中文字幕| 尾随美女入室| 国产精品一区二区在线观看99| 久久99一区二区三区| 丝袜人妻中文字幕| 看非洲黑人一级黄片| 色吧在线观看| 久久久精品94久久精品| 精品少妇久久久久久888优播| 99久久99久久久精品蜜桃| 亚洲欧美清纯卡通| 99热全是精品| 久久久久久久国产电影| 亚洲国产看品久久| 男女床上黄色一级片免费看| 国产在线免费精品| 人成视频在线观看免费观看| 久久精品亚洲av国产电影网| 观看av在线不卡| 日韩一卡2卡3卡4卡2021年| 久久久国产一区二区| 黑人欧美特级aaaaaa片| 热re99久久国产66热| 777久久人妻少妇嫩草av网站| 欧美 日韩 精品 国产| 免费少妇av软件| 国产毛片在线视频| 老司机影院毛片| 午夜日本视频在线| 日韩精品有码人妻一区| 国产一区亚洲一区在线观看| 丰满迷人的少妇在线观看| 好男人视频免费观看在线| 狂野欧美激情性bbbbbb| 三上悠亚av全集在线观看| 一区福利在线观看| 国产免费现黄频在线看| 日韩av不卡免费在线播放| 伊人亚洲综合成人网| 成人国产麻豆网| 成人影院久久| 午夜91福利影院| 亚洲精品成人av观看孕妇| 国产日韩欧美亚洲二区| 国产精品二区激情视频| 黄频高清免费视频| 一区二区三区精品91| 成人国产av品久久久| 欧美中文综合在线视频| 人人妻人人添人人爽欧美一区卜| 王馨瑶露胸无遮挡在线观看| 一级黄片播放器| 国产野战对白在线观看| 国产伦理片在线播放av一区| 欧美日韩综合久久久久久| 电影成人av| 一级毛片 在线播放| 久久久久久久久免费视频了| 在线 av 中文字幕| 中文字幕制服av| 国产精品久久久人人做人人爽| 国产精品.久久久| 国产精品 国内视频| 国产一区二区在线观看av| 黑人猛操日本美女一级片| 亚洲国产最新在线播放| 女人被躁到高潮嗷嗷叫费观| 在线观看免费视频网站a站| 成人午夜精彩视频在线观看| 看十八女毛片水多多多| 日韩大码丰满熟妇| 少妇的丰满在线观看| 两个人看的免费小视频| 9191精品国产免费久久| 少妇人妻精品综合一区二区| 亚洲伊人色综图| 国产成人一区二区在线| 亚洲成人av在线免费| 国产黄频视频在线观看| 777米奇影视久久| 好男人视频免费观看在线| 超碰97精品在线观看| 亚洲精品久久成人aⅴ小说| 丝袜美腿诱惑在线| 十八禁人妻一区二区| 欧美精品一区二区免费开放| e午夜精品久久久久久久| 日本爱情动作片www.在线观看| 少妇 在线观看| 美女午夜性视频免费| 国产免费现黄频在线看| 国产精品一国产av| 久久久久人妻精品一区果冻| 成年人午夜在线观看视频| 乱人伦中国视频| 免费在线观看完整版高清| 欧美在线一区亚洲| 99热网站在线观看| 日日爽夜夜爽网站| 久久久精品区二区三区| 秋霞伦理黄片| 亚洲第一区二区三区不卡| 操出白浆在线播放| 国产日韩一区二区三区精品不卡| 不卡av一区二区三区| 午夜av观看不卡| 成人手机av| 国产精品三级大全| 夜夜骑夜夜射夜夜干| 国产爽快片一区二区三区| 亚洲,欧美精品.| 亚洲国产精品一区二区三区在线| 又大又黄又爽视频免费| 久久综合国产亚洲精品| 女性被躁到高潮视频| av在线app专区| 色94色欧美一区二区| 亚洲熟女精品中文字幕| www.熟女人妻精品国产| 无限看片的www在线观看| 亚洲av成人不卡在线观看播放网 | 亚洲精品日本国产第一区| 97在线人人人人妻| 黄色一级大片看看| 美女大奶头黄色视频| av又黄又爽大尺度在线免费看| 日日爽夜夜爽网站| 亚洲五月色婷婷综合| 欧美97在线视频| 亚洲成人av在线免费| av免费观看日本| 久久ye,这里只有精品| 水蜜桃什么品种好| 黑人巨大精品欧美一区二区蜜桃| 国产成人一区二区在线| 亚洲人成电影观看| 成人三级做爰电影| 日韩免费高清中文字幕av| 熟妇人妻不卡中文字幕| 一本一本久久a久久精品综合妖精| 一级,二级,三级黄色视频| 国产精品亚洲av一区麻豆 | 亚洲第一区二区三区不卡| 久久这里只有精品19| 久久国产亚洲av麻豆专区| 蜜桃国产av成人99| 国产成人午夜福利电影在线观看| 欧美在线一区亚洲| 国产成人欧美在线观看 | 欧美日韩福利视频一区二区| 国产精品女同一区二区软件| 一边摸一边抽搐一进一出视频| 亚洲精品中文字幕在线视频| 国产xxxxx性猛交| 欧美xxⅹ黑人| 菩萨蛮人人尽说江南好唐韦庄| 久久韩国三级中文字幕| 国产日韩欧美视频二区| 蜜桃在线观看..| 午夜日本视频在线| av一本久久久久| 欧美精品亚洲一区二区| 日韩电影二区| 精品福利永久在线观看| xxx大片免费视频| 午夜福利影视在线免费观看| 伊人久久国产一区二区| 久久久久久久久久久免费av| 久久久久精品久久久久真实原创| xxxhd国产人妻xxx| 亚洲三区欧美一区| 国产精品麻豆人妻色哟哟久久| 你懂的网址亚洲精品在线观看| 操美女的视频在线观看| av在线播放精品| 亚洲精品中文字幕在线视频| 欧美黑人欧美精品刺激| 香蕉国产在线看| 青青草视频在线视频观看| 亚洲自偷自拍图片 自拍| 亚洲少妇的诱惑av| av有码第一页| 久久久久精品久久久久真实原创| 日日爽夜夜爽网站| 午夜精品国产一区二区电影| 免费在线观看完整版高清| 又大又爽又粗| 国产一区二区 视频在线| 青春草国产在线视频| 尾随美女入室| 免费高清在线观看视频在线观看| 制服人妻中文乱码| 国产精品国产三级国产专区5o| 亚洲视频免费观看视频| 亚洲精品自拍成人| 久久午夜综合久久蜜桃| 日韩中文字幕欧美一区二区 | 黑丝袜美女国产一区| 极品人妻少妇av视频| 国产免费一区二区三区四区乱码| 少妇人妻久久综合中文| 黄色怎么调成土黄色| 老汉色∧v一级毛片| 黄色一级大片看看| 久久这里只有精品19| 国语对白做爰xxxⅹ性视频网站| 免费黄频网站在线观看国产| 国产成人精品久久久久久| 18禁裸乳无遮挡动漫免费视频| 欧美 亚洲 国产 日韩一| 国产又爽黄色视频| 日韩人妻精品一区2区三区| 久久久精品94久久精品| 亚洲欧美激情在线| 熟女av电影| 大香蕉久久成人网| 色视频在线一区二区三区| 国产成人欧美| 国产一区二区激情短视频 | 蜜桃在线观看..| 亚洲美女黄色视频免费看| 久久97久久精品| 亚洲av成人不卡在线观看播放网 | 丝袜在线中文字幕| 国产精品久久久人人做人人爽| 久久久亚洲精品成人影院| 搡老乐熟女国产| 亚洲精品国产区一区二| 久久久久精品人妻al黑| 少妇人妻 视频| 亚洲综合精品二区| 日本欧美视频一区| 母亲3免费完整高清在线观看| 制服人妻中文乱码| 亚洲国产成人一精品久久久| 黄色一级大片看看| 色吧在线观看| 日韩 亚洲 欧美在线| 一本一本久久a久久精品综合妖精| 国产熟女午夜一区二区三区| www.精华液| 超碰97精品在线观看| 亚洲成av片中文字幕在线观看| 日本午夜av视频| 国产成人免费无遮挡视频| 97在线人人人人妻| 成年美女黄网站色视频大全免费| 2018国产大陆天天弄谢| 侵犯人妻中文字幕一二三四区| 亚洲欧美成人精品一区二区| 欧美日韩一级在线毛片| 亚洲婷婷狠狠爱综合网| 国产xxxxx性猛交| 亚洲人成电影观看| 国产99久久九九免费精品| 人人妻,人人澡人人爽秒播 | 大香蕉久久网| 成人国产av品久久久| 九色亚洲精品在线播放| 18禁国产床啪视频网站| 18在线观看网站| 久久韩国三级中文字幕| 国产精品麻豆人妻色哟哟久久| 欧美在线一区亚洲| 日本爱情动作片www.在线观看| 国产爽快片一区二区三区| 国产精品av久久久久免费| 亚洲成国产人片在线观看| 丝袜脚勾引网站| 尾随美女入室| 国产欧美日韩综合在线一区二区| 国产在视频线精品| 成年动漫av网址| 一区在线观看完整版| 亚洲综合精品二区| 亚洲人成网站在线观看播放| 日本av免费视频播放| 色播在线永久视频| 秋霞伦理黄片| 中文欧美无线码| 久久 成人 亚洲| 日本欧美视频一区| 在线观看免费日韩欧美大片| 久久青草综合色| 女的被弄到高潮叫床怎么办| 一级毛片黄色毛片免费观看视频| www.熟女人妻精品国产| 悠悠久久av| 中文字幕人妻丝袜制服| 精品一区在线观看国产| 在线观看www视频免费| 啦啦啦 在线观看视频| 在线观看三级黄色| 久久久久久免费高清国产稀缺| 久久精品人人爽人人爽视色| 国产精品人妻久久久影院| 麻豆精品久久久久久蜜桃| 亚洲av日韩精品久久久久久密 | 久久天躁狠狠躁夜夜2o2o | 男人爽女人下面视频在线观看| 亚洲视频免费观看视频| 亚洲伊人久久精品综合| 天天躁日日躁夜夜躁夜夜| 99国产综合亚洲精品| 纵有疾风起免费观看全集完整版| 黑人欧美特级aaaaaa片| 亚洲精品一区蜜桃| 中文字幕人妻丝袜制服| 激情五月婷婷亚洲| 欧美亚洲日本最大视频资源| 人体艺术视频欧美日本| 亚洲av福利一区| 亚洲国产欧美日韩在线播放| 久久人妻熟女aⅴ| 久久 成人 亚洲| 久久精品aⅴ一区二区三区四区| 亚洲熟女精品中文字幕| 日日爽夜夜爽网站| 日韩欧美一区视频在线观看| 亚洲色图 男人天堂 中文字幕| 又黄又粗又硬又大视频| 精品人妻熟女毛片av久久网站| 菩萨蛮人人尽说江南好唐韦庄| 国产老妇伦熟女老妇高清| 欧美日本中文国产一区发布| 少妇的丰满在线观看| 婷婷色av中文字幕| 亚洲第一区二区三区不卡| 久久综合国产亚洲精品| 涩涩av久久男人的天堂| 国产精品偷伦视频观看了| 不卡av一区二区三区| 人人妻人人爽人人添夜夜欢视频| 亚洲国产欧美在线一区| 欧美 日韩 精品 国产| 欧美精品人与动牲交sv欧美| 日本爱情动作片www.在线观看| 青春草视频在线免费观看| 日韩精品有码人妻一区| 老司机影院成人| 日韩,欧美,国产一区二区三区| 一级片'在线观看视频| 亚洲国产成人一精品久久久| 欧美在线一区亚洲| 毛片一级片免费看久久久久| 午夜免费男女啪啪视频观看| 日韩一卡2卡3卡4卡2021年| 天堂俺去俺来也www色官网| 国产欧美亚洲国产| 多毛熟女@视频| 精品亚洲乱码少妇综合久久| 一级,二级,三级黄色视频| 亚洲国产最新在线播放| 午夜福利一区二区在线看| 99精品久久久久人妻精品| 免费在线观看视频国产中文字幕亚洲 | 久久av网站| 日本色播在线视频| 亚洲欧美一区二区三区国产| 天天添夜夜摸| 天堂8中文在线网| 波多野结衣av一区二区av| 婷婷色综合大香蕉| 欧美日韩av久久| 夫妻性生交免费视频一级片| 制服人妻中文乱码| 日本91视频免费播放| 一区在线观看完整版| 国产精品99久久99久久久不卡 | 日韩欧美精品免费久久| 免费在线观看完整版高清| 国产亚洲欧美精品永久| 啦啦啦啦在线视频资源| 亚洲欧美一区二区三区黑人| 国产免费视频播放在线视频| 宅男免费午夜| 爱豆传媒免费全集在线观看| 18禁观看日本| 亚洲精品日韩在线中文字幕| 热99国产精品久久久久久7| 美女高潮到喷水免费观看| 久久ye,这里只有精品| 欧美精品高潮呻吟av久久| 又大又黄又爽视频免费| 韩国高清视频一区二区三区| 纵有疾风起免费观看全集完整版| 国产精品香港三级国产av潘金莲 | 最新在线观看一区二区三区 | 九色亚洲精品在线播放| 天天添夜夜摸| 最近中文字幕高清免费大全6| 1024香蕉在线观看| 日本猛色少妇xxxxx猛交久久| 久久久国产欧美日韩av| 69精品国产乱码久久久| 777久久人妻少妇嫩草av网站| 国产精品av久久久久免费| 国产免费视频播放在线视频| 亚洲在久久综合| 国产亚洲午夜精品一区二区久久| 亚洲少妇的诱惑av| 美女主播在线视频| 国产亚洲一区二区精品| 亚洲第一区二区三区不卡| 午夜日韩欧美国产| 九草在线视频观看| 国产一区二区 视频在线| 亚洲成人免费av在线播放| 亚洲国产av新网站| 晚上一个人看的免费电影| 欧美日韩成人在线一区二区| 国产黄色视频一区二区在线观看| 人体艺术视频欧美日本| 国产成人免费无遮挡视频| 高清视频免费观看一区二区| 日韩伦理黄色片| 99久久综合免费| 欧美日本中文国产一区发布| www.熟女人妻精品国产| 美女午夜性视频免费| 国产精品成人在线| 中文字幕av电影在线播放| 亚洲成人国产一区在线观看 | 99香蕉大伊视频| 超碰97精品在线观看| 丝袜人妻中文字幕| 国产成人a∨麻豆精品| 欧美激情极品国产一区二区三区| 亚洲av电影在线观看一区二区三区| 日韩 欧美 亚洲 中文字幕| 亚洲精品乱久久久久久| 亚洲少妇的诱惑av| 岛国毛片在线播放| 90打野战视频偷拍视频| 欧美中文综合在线视频| 少妇精品久久久久久久| 91aial.com中文字幕在线观看| 国产精品免费视频内射| 成人毛片60女人毛片免费| 人妻人人澡人人爽人人| 欧美日韩综合久久久久久| 人人妻人人添人人爽欧美一区卜| 国产av一区二区精品久久| 午夜激情久久久久久久| tube8黄色片| 视频区图区小说| 三上悠亚av全集在线观看| 国产 一区精品| 韩国av在线不卡| 亚洲人成电影观看| 最近2019中文字幕mv第一页| 最近最新中文字幕大全免费视频 | 国产乱来视频区| 日韩av在线免费看完整版不卡| 国产精品99久久99久久久不卡 | 黄色视频在线播放观看不卡| 在线看a的网站| 国产精品久久久久成人av| 国产有黄有色有爽视频| 亚洲综合精品二区| 99精国产麻豆久久婷婷| 美国免费a级毛片| 中文字幕人妻丝袜制服| 精品久久久精品久久久| 久久国产精品大桥未久av| 9191精品国产免费久久| 国产成人精品在线电影| 亚洲人成电影观看| 久久久久网色| 亚洲第一区二区三区不卡| 欧美人与性动交α欧美软件| 亚洲第一区二区三区不卡| 亚洲四区av| 亚洲欧美精品自产自拍| 久久精品亚洲熟妇少妇任你| 你懂的网址亚洲精品在线观看| 国产精品久久久久成人av| 制服丝袜香蕉在线| 曰老女人黄片| 纵有疾风起免费观看全集完整版| 老司机在亚洲福利影院| 美女大奶头黄色视频| a级毛片黄视频| 日韩人妻精品一区2区三区| a 毛片基地| 精品福利永久在线观看| 亚洲精品日本国产第一区| 国产福利在线免费观看视频| av女优亚洲男人天堂| 日韩av不卡免费在线播放| 精品少妇一区二区三区视频日本电影 | 色网站视频免费| 一本久久精品| 精品免费久久久久久久清纯 | 综合色丁香网| 亚洲男人天堂网一区| 看免费av毛片| 亚洲av日韩精品久久久久久密 | 人人澡人人妻人| 不卡视频在线观看欧美| 一区二区三区四区激情视频| 大香蕉久久网| 99久久99久久久精品蜜桃| 叶爱在线成人免费视频播放| 国语对白做爰xxxⅹ性视频网站| 欧美日韩一区二区视频在线观看视频在线| 又黄又粗又硬又大视频| 国产成人精品无人区| 少妇人妻精品综合一区二区| 亚洲精品,欧美精品| 亚洲国产欧美日韩在线播放| 少妇被粗大猛烈的视频| 十八禁高潮呻吟视频| 狂野欧美激情性xxxx| 日韩精品免费视频一区二区三区| 精品一区在线观看国产| 宅男免费午夜| 亚洲色图综合在线观看| 亚洲av电影在线进入| 国产在视频线精品| 永久免费av网站大全| 亚洲精品成人av观看孕妇| 亚洲美女黄色视频免费看| 男人舔女人的私密视频| 最近最新中文字幕免费大全7| 国产亚洲精品第一综合不卡| 午夜激情久久久久久久| 热99久久久久精品小说推荐| 久久人人爽av亚洲精品天堂| 久久精品人人爽人人爽视色| 无限看片的www在线观看| a级毛片在线看网站| 另类精品久久| 看非洲黑人一级黄片| 国产亚洲精品第一综合不卡| 国产亚洲欧美精品永久| 亚洲视频免费观看视频| 黄色怎么调成土黄色| 亚洲成色77777| 色婷婷久久久亚洲欧美| 国产福利在线免费观看视频| 韩国av在线不卡| 看非洲黑人一级黄片| 在线观看人妻少妇| 国产亚洲欧美精品永久| 亚洲精品一二三| 日本av手机在线免费观看| 九草在线视频观看| 国产精品久久久久成人av| 亚洲精品自拍成人| 久久久久久免费高清国产稀缺| 久久这里只有精品19| 午夜影院在线不卡| 国产精品一区二区在线不卡| 深夜精品福利| 亚洲成人手机| 大话2 男鬼变身卡| 麻豆乱淫一区二区| 99精品久久久久人妻精品| 丝瓜视频免费看黄片| 各种免费的搞黄视频| 少妇被粗大的猛进出69影院| 久久综合国产亚洲精品| 免费高清在线观看日韩| 午夜激情av网站| 男人操女人黄网站| av线在线观看网站| 亚洲综合色网址| 欧美激情 高清一区二区三区| 中文字幕av电影在线播放| 男女之事视频高清在线观看 | 高清欧美精品videossex| 日韩av在线免费看完整版不卡| 18禁动态无遮挡网站| 久久久久国产一级毛片高清牌| 久久亚洲国产成人精品v| 黄色毛片三级朝国网站| 在线看a的网站| 国产一级毛片在线| 国产成人欧美在线观看 | 国产深夜福利视频在线观看| 日韩av免费高清视频| 男女床上黄色一级片免费看| 久久精品国产综合久久久|