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

    跨聲速風(fēng)洞中使用短軸探管測量試驗段核心流馬赫數(shù)影響研究

    2022-04-27 01:46:10鄧海均羅新福洪少尊王偉仲
    關(guān)鍵詞:探管馬赫數(shù)試驗段

    鄧海均,熊 波,羅新福,洪少尊,李 強,王偉仲,劉 俊

    (中國空氣動力研究與發(fā)展中心高速空氣動力研究所,綿陽 621000)

    一直以來,風(fēng)洞試驗都是對飛行器進行空氣動力學(xué)研究最為有效的手段之一。跨聲速風(fēng)洞試驗馬赫數(shù)范圍通常為0.4~1.4,是現(xiàn)代戰(zhàn)機進行空戰(zhàn)的主要速度范圍,其他如火箭、導(dǎo)彈等雖然速度遠大于聲速,但在其加速過程中仍然必須經(jīng)過跨聲速范圍[1]。所有這些飛行器必須要在跨聲速風(fēng)洞中開展大量的試驗研究,從而獲取準確的氣動特性數(shù)據(jù),為研發(fā)設(shè)計提供數(shù)據(jù)支持。風(fēng)洞優(yōu)良的流場品質(zhì)是飛行器獲得高質(zhì)量試驗數(shù)據(jù)的前提,更是評判其性能的重要指標。新建跨聲速風(fēng)洞在投入試驗運行前應(yīng)按要求完成速度場、方向場、跨聲速通氣壁消波特性、洞壁邊界層、氣流噪聲、湍流度等項目的流場校測內(nèi)容[2?3],當(dāng)風(fēng)洞進行洞體大修、測控系統(tǒng)改造、測量設(shè)備有重大更改、發(fā)現(xiàn)模型試驗質(zhì)量有問題且懷疑與流場品質(zhì)有關(guān)時,均應(yīng)進行相應(yīng)的流場校測,其中速度場是最為核心的項目。

    常規(guī)速度場校測通常采用軸探管測得的靜壓與風(fēng)洞總壓計算得到試驗段中心線(核心流)上的馬赫數(shù),最終獲得核心流馬赫數(shù)軸向分布特性和風(fēng)洞試驗馬赫數(shù)與駐室馬赫數(shù)的修正關(guān)系。為減小頭錐對測壓點的干擾,通常將軸探管前伸至收縮段,為確保試驗安全,保證軸探管中心線與風(fēng)洞軸線重合以及軸探管的強度和剛度,需要在噴管或試驗段多個位置用鋼繩將軸探管拉緊固定(圖1),設(shè)備安裝工序復(fù)雜且耗時較長,對試驗段的流場也存在一定的干擾[4?7]。同時,張線和固定工裝需在噴管段和試驗段開孔或開槽,對洞體結(jié)構(gòu)也存在一定不利的影響。短軸探管作為一種流場檢測的有效工具,在國外的大型跨聲速風(fēng)洞流場測試中受到較高的重視,日本JAXA 中心連續(xù)式跨聲速風(fēng)洞JT?WT、美國波音跨聲速風(fēng)洞BTWT 和AMES 研究中心跨聲速風(fēng)洞UPWT 的研究人員均利用短軸探管開展了相關(guān)研究性試驗[8?10]。為進一步探索短軸探管在跨聲速試驗段中的應(yīng)用效果,文中設(shè)計了一種頭部位于試驗段內(nèi)的新型短軸探管,無需張線和固定工裝,安裝相對方便、簡單,測量區(qū)域可覆蓋整個試驗段模型區(qū)核心流馬赫數(shù)范圍。同時,依據(jù)研究計劃開展了短軸探管數(shù)值仿真計算工作,并在CARDC 0.6 m連續(xù)式跨聲速風(fēng)洞中進行了相關(guān)試驗。

    圖1 某跨聲速風(fēng)洞中的長軸探管Fig.1 Long centerline probe in a transonic wind tunnel

    1 數(shù)值仿真

    1.1 計算模型建立

    為研究短軸探管頭部在不同馬赫數(shù)下對流場的影響規(guī)律,針對馬赫數(shù)0.95,1.0,1.2 三個典型馬赫數(shù)狀態(tài)開展了數(shù)值模擬研究。

    數(shù)值計算所采用的短軸探管直徑為30 mm,柱段長度為1 200 mm,頭部為與8°尖錐型頭部等長的改型圓弧頭部。采用Gridgen 軟件生成六面體的結(jié)構(gòu)化網(wǎng)格,如圖2 所示。為提高計算精度,在短軸探管頭部最尖點處截取直徑為0.4 mm 的平臺進行處理,以避免極性軸網(wǎng)格的出現(xiàn)。短軸探管柱體長度一直延伸至遠場邊界,靠近壁面的區(qū)域網(wǎng)格進行加密處理,以更好地捕捉近壁區(qū)域的邊界層。計算域的入口和出口條件均設(shè)為壓力遠場邊界條件。

    圖2 短軸探管計算網(wǎng)格Fig.2 Calculation grid of short centerline probe

    1.2 控制方程及湍流模型

    在求解N?S 方程時,F(xiàn)luent 使用有限體積法:首先將計算區(qū)域劃分為網(wǎng)格,使各個網(wǎng)格點四周的控制體積保持不重復(fù)的狀態(tài),然后將待求解的微分方程對每個控制體積積分求解。

    湍流模型的作用是使方程組可以采用封閉的計算方法。湍流模擬采用k?ωSST 湍流模型,該模型使用混合函數(shù)從壁面附近的標準k?ω模型逐漸過渡到邊界層外部高雷諾數(shù)k?ω模型,模型中包含修正的湍流黏性公式用以解決湍流剪應(yīng)力引起的輸運效果,并且湍流黏度考慮了湍流剪應(yīng)力的傳播。

    1.3 計算求解

    計算模型建立后,利用Fluent 求解器進行求解[11],采用基于密度的耦合隱式求解法,壁面條件為無滑移條件,所有其他標量采用不可滲透壁面條件,動量選擇二階迎風(fēng)格式,湍流動能及湍流耗散率選用一階迎風(fēng)格式,數(shù)值仿真計算中按照Y+=1的原則模擬壁面第一層網(wǎng)格高度為8.6e-3,分別以68 萬個、300 萬個、680 萬個網(wǎng)格數(shù)進行了網(wǎng)格無關(guān)性驗證,判斷依據(jù)為參考測點壓力系數(shù)相差值小于0.01,驗證結(jié)果表明68 萬個網(wǎng)格數(shù)足夠滿足仿真要求。

    2 風(fēng)洞試驗

    2.1 試驗平臺

    試驗風(fēng)洞是一座由AV90?3 型軸流式壓縮機驅(qū)動的連續(xù)式跨聲速風(fēng)洞,水平布置全鋼結(jié)構(gòu)。風(fēng)洞氣動輪廓如圖3 所示。風(fēng)洞本體主要包括風(fēng)洞主回路和風(fēng)洞輔助系統(tǒng)。風(fēng)洞主回路由穩(wěn)定段、收縮段、噴管段、試驗段、模型支架段、二喉道段、再導(dǎo)入段、動力段、冷卻器以及第一、二擴散段和4 個拐角段等部段組成[12]。

    2.2 性能參數(shù)

    0.6 m 連續(xù)式跨聲速風(fēng)洞主要性能參數(shù)為[13?15]:

    (1)試驗段尺寸:0.6 m(寬)×0.6 m(高)×1.85 m(長);

    (2)試驗段馬赫數(shù):0.2~1.6;

    (3)馬赫數(shù)控制精度:≤0.002;

    (4)穩(wěn)定段總壓:(0.15~2.5)×105Pa;

    (5)總壓控制精度:≤0.2%;

    (6)氣流總溫:273~323 K;

    (7)試驗雷諾數(shù):Re·c=(0.1~2.25)×106(c=0.06 m);

    (8)風(fēng)洞試驗時間:連續(xù)運行不小于2 h;

    (9)試驗段模型區(qū)截面氣流溫度均勻性:|ΔT0|≤1~2 K;

    (10)試驗段氣流溫度穩(wěn)定性:單條極曲級|ΔT0|≤1 K ;

    (11)風(fēng)洞壓縮機軸功率:主壓縮機軸功率Nm≤3.8 MW;

    (12)輔壓縮機軸功率Na≤1.5 MW。

    2.3 試驗條件及測試設(shè)備

    試驗全部采用常壓方式運行,穩(wěn)定段總壓控制在100 kPa。分別選用孔壁試驗段和槽壁試驗段開展試驗,如圖4、5 所示??妆谠囼灦蔚拈_孔方式為:距離試驗段入口120~520 mm,6%;距離試驗段入口520~1 600 mm,4%;距離試驗段入口1 600~1 850 mm,6%。槽壁試驗段的開槽方式為:6 條槽,開閉比為6%。

    圖5 槽壁試驗段Fig.5 Test section of groove wall

    軸探管各測點的壓力信號通過掃描閥進行測量,為了保證壓力測量的同步性,將穩(wěn)定段總壓引至掃描閥模塊進行采集。常壓工況時,掃描閥模塊量程為±15 psi;掃描閥測量精度為0.05%。

    試驗采用代號C?T 的長軸探管,其主要由測壓段、延伸段以及頭錐組成,頭錐位于收縮段低速來流中,對試驗結(jié)果基本無影響。直徑d=40 mm,在風(fēng)洞中堵塞度為0.35%,試驗結(jié)果可作為短軸探管的參照標準。長軸探管尾部固定安裝在流場校測專用中部支架上,在收縮段入口設(shè)置有上下兩根鋼絲繩以提高軸探管的剛度并實現(xiàn)管體沿風(fēng)洞軸線方向的調(diào)節(jié)功能。如圖6、7 所示,C?T 在風(fēng)洞中安裝完成后全長5 350 mm,測點覆蓋距離試驗段入口520~1 850 mm 的核心流區(qū)域,其中模型區(qū)前后(距離試驗段入口640~1 690 mm 區(qū)域)進行了加密處理,測點間距為25 mm,其他區(qū)域測點間距為40 mm,測壓孔徑0.5 mm,測點50 個,采用上下交錯開孔形式,模型區(qū)定義為距離試驗段入口960~1 560 mm 區(qū)域。

    圖6 長軸探管C-T 尺寸示意圖Fig.6 Dimension diagram of long centerline probe C-T

    圖7 風(fēng)洞長軸探管C-TFig.7 Long centerline probe C-T in wind tunnel

    為方便對比,代號DYH?T 的短軸探管采用同一測壓段,頭部以8°尖錐形頭部長度為基準設(shè)計了等長的圓弧型頭部(圖8),實現(xiàn)與柱段光滑過渡,可在一定程度上減小對試驗段高速流場的擾動。短軸探管安裝完成后總長1 881 mm,在風(fēng)洞中如圖9 所示。

    圖8 短軸探管DYH-T 尺寸示意圖Fig.8 Dimension diagram of short centerline probe DYH-T

    圖9 風(fēng)洞短軸探管DYH-TFig.9 Short centerline probe DYH-T in wind tunnel

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

    在亞跨聲速,軸探管測點和駐室測點的靜壓通過靜壓管接到掃描閥進行測量,根據(jù)穩(wěn)定段總壓及各測點靜壓,依據(jù)馬赫數(shù)計算公式得出相應(yīng)的馬赫數(shù)

    由式(8)可知,當(dāng)均值為0 時表示兩次測值曲線完全重合。

    3 數(shù)值計算結(jié)果分析

    對馬赫數(shù)0.95,1.0,1.2 三個典型馬赫數(shù)狀態(tài)開展了數(shù)值計算,得到了遠場工況下短軸探管靜壓及馬赫數(shù)分布情況,如圖10~15 所示,統(tǒng)一選取X?Y截面進行分析。

    圖10 馬赫數(shù)0.95 的短軸探管靜壓分布Fig.10 Static pressure distribution of short centerline probe with Ma=0.95

    圖11 馬赫數(shù)0.95 的短軸探管馬赫數(shù)分布Fig.11 Mach number distribution of short centerline probe with Ma=0.95

    圖12 馬赫數(shù)1.0 的短軸探管靜壓分布Fig.12 Static pressure distribution of short centerline probe with Ma=1.0

    圖13 馬赫數(shù)1.0 的短軸探管馬赫數(shù)分布Fig.13 Mach number distribution of short centerline probe with Ma=1.0

    圖14 馬赫數(shù)1.2 的短軸探管靜壓分布Fig.14 Static pressure distribution of short centerline probe with Ma=1.2

    圖15 馬赫數(shù)1.2 的短軸探管馬赫數(shù)分布Fig.15 Mach number distribution of short centerline probe with Ma=1.2

    馬赫數(shù)0.95時來流總壓和靜壓分別為100 000 Pa和55 946 Pa。從數(shù)值計算結(jié)果中可以看出:在圓弧形頭部區(qū)靜壓值升高,馬赫數(shù)降低,自肩部以后靜壓逐漸降低然后升高,馬赫數(shù)逐步升高然后降低,頭部對靜壓和馬赫數(shù)的影響從前往后遞減,最終很快恢復(fù)至來流馬赫數(shù)水平。柱段馬赫數(shù)與來流馬赫數(shù)差量<0.001 5 時靜壓值為55 852 Pa,位置距離肩部3.4 倍管徑左右。

    馬赫數(shù)1.0時來流總壓和靜壓分別為100 000 Pa和52 828 Pa。從數(shù)值計算結(jié)果中可以看出:該馬赫數(shù)下頭部產(chǎn)生與來流幾乎垂直的激波和膨脹波,影響區(qū)域略擴大。柱段馬赫數(shù)與來流馬赫數(shù)差量<0.001 5 時靜壓值為52 736 Pa,位置距離肩部3.5 倍管徑左右。

    馬赫數(shù)1.2時來流總壓和靜壓分別為100 000 Pa和41 238 Pa。從數(shù)值計算結(jié)果中可以看出:該馬赫數(shù)下頭部產(chǎn)生了較強的斜激波和膨脹波,頭部對馬赫數(shù)的影響從前往后遞減。柱段馬赫數(shù)與來流馬赫數(shù)差量<0.001 5 時靜壓值為41 158 Pa,位置距離肩部8.5 倍管徑左右。

    數(shù)值計算結(jié)果表明,在馬赫數(shù)1.0 以下時,短軸探管頭部在流場中引起的干擾相對較小,馬赫數(shù)很快就能恢復(fù)到來流馬赫數(shù)水平,短軸探管頭部在流場中不會產(chǎn)生激波和膨脹波。當(dāng)馬赫數(shù)等于1.0 時,頭部產(chǎn)生的激波和膨脹波強度相對較弱且波系方向接近垂直于來流方向。當(dāng)馬赫數(shù)大于1.0 時,短軸探管頭部在低超聲速來流中會產(chǎn)生相對較強的激波和膨脹波,頭部干擾區(qū)域相對于馬赫數(shù)1.0 及以下狀態(tài)急劇增加,干擾區(qū)域隨著馬赫數(shù)增大而增大,從而引起流場產(chǎn)生較大變化,馬赫數(shù)恢復(fù)到與來流馬赫數(shù)相當(dāng)?shù)乃叫枰木嚯x也越長。

    由于孔壁/槽壁試驗段條件的跨聲速流場難以精確模擬,文中僅考慮遠場邊界條件進行數(shù)值仿真,仿真定量分析結(jié)果與風(fēng)洞試驗結(jié)果存在一定差異,可為風(fēng)洞試驗數(shù)據(jù)分析提供一定參考。

    4 風(fēng)洞試驗結(jié)果分析

    在孔壁試驗段和槽壁試驗段中分別使用長軸探管和短軸探管開展風(fēng)洞試驗,試驗馬赫數(shù)范圍為0.8~1.4,長軸探管試驗結(jié)果作為參照值。表1給出了短軸探管與長軸探管在模型區(qū)各對應(yīng)測點差量的平均值數(shù)據(jù)。試驗結(jié)果表明,在孔壁試驗段和槽壁試驗段中長、短軸探管的測值差異均隨著馬赫數(shù)上升而增大,且在馬赫數(shù)1.0 以上時,其差異幅度迅速擴大。相對而言,孔壁試驗段條件下的平均值差異幅度整體上明顯低于槽壁試驗段結(jié)果。

    Table 1 模型區(qū)馬赫數(shù)測點差異均值結(jié)果Table 1 Mean difference of Mach number in model region

    為進一步具體分析二者的差異特性,圖16~19 給出了馬赫數(shù)0.95 和1.1 條件下短軸探管與長軸探管的全區(qū)域測點分布差異曲線。

    圖16 孔壁試驗段核心流馬赫數(shù)分布曲線差異(Ma=0.95)Fig.16 Difference of Mach number distribution curves of core flow in test section of hole wall(Ma=0.95)

    圖17 槽壁試驗段核心流馬赫數(shù)分布曲線差異(Ma=0.95)Fig.17 Difference of Mach number distribution curves of core flow in test section of groove wall(Ma =0.95)

    在Ma=0.95 時,長、短軸探管的測值呈現(xiàn)的流場波動規(guī)律一致性較好,孔壁試驗段及槽壁試驗段條件下模型區(qū)馬赫數(shù)測點差異均值較小,從測值曲線上看均無明顯變化。說明在該馬赫數(shù)下短軸探管頭部引起的干擾對試驗段模型區(qū)流場影響較小,核心流馬赫數(shù)測量結(jié)果準度和精度無明顯差異。

    圖18 孔壁試驗段核心流馬赫數(shù)分布曲線差異(Ma =1.1)Fig.18 Difference of Mach number distribution curves of core flow in test section of hole wall(Ma =1.1)

    圖19 槽壁試驗段核心流馬赫數(shù)分布曲線差異(Ma=1.1)Fig.19 Difference of Mach number distribution curves of core flow in test section of groove wall(Ma=1.1)

    在Ma=1.1 時,孔壁試驗段條件下長、短軸探管的全部測值呈現(xiàn)的流場波動規(guī)律仍具有較好的一致性,但局部流場的細節(jié)存在一定差異,模型區(qū)馬赫數(shù)測點差異均值為0.003 4。槽壁試驗段條件下短軸探管的測值曲線與長軸探管差異明顯,模型區(qū)馬赫數(shù)測點差異均值為0.004 8,顯著高于孔壁試驗段結(jié)果。

    圖20、21 給出了孔壁試驗段和槽壁試驗段條件下長、短軸探管模型區(qū)內(nèi)馬赫數(shù)均方根偏差對比結(jié)果。從模型區(qū)均方根偏差指標來看,長、短軸探管在孔壁試驗段和槽壁試驗段中的測試結(jié)果均隨馬赫數(shù)增加呈上升趨勢,在Ma≤0.95 時,長、短軸探管測值結(jié)果中均方根偏差指標差異≤0.000 2,整體上無明顯差異;在Ma=1.0 時,孔壁試驗段和槽壁試驗段中長、短軸探管的均方根偏差均存在較小差異,相對于亞音速狀態(tài)小幅增加,主要原因是Ma=1.0 時主要產(chǎn)生正激波和膨脹波,由于激波角激波呈90°使得在較短的軸向距離內(nèi),激波就可以形成多次壁面反射,而通氣壁具有一定的消波能力,因此經(jīng)過一定軸向距離后,通氣壁較好地消除了模型頭部和肩部產(chǎn)生的大部分激波和膨脹波;而當(dāng)1.0<Ma≤1.4 時,隨著馬赫數(shù)的增加,激波角的減小及激波強度的增加,經(jīng)過頭部和肩部所產(chǎn)生的激波和膨脹波沿風(fēng)洞軸線方向傳遞的距離更遠,這些波系經(jīng)過壁面反射的現(xiàn)象也是如此,但隨著軸線方向通氣壁對波系的消波作用,使得頭部產(chǎn)生的激波和膨脹波對風(fēng)洞速度場的影響呈現(xiàn)先強后弱的規(guī)律。要想削弱這些波系對流場的擾動,需要風(fēng)洞試驗段具有更長的消波區(qū)域,且馬赫數(shù)越高,長度越長。而對于風(fēng)洞來說,模型旋轉(zhuǎn)中心是固定的,即用于風(fēng)洞試驗的模型區(qū)是固定的。因此相同長度的通氣壁,馬赫數(shù)越高,消波能力越低,模型頭部產(chǎn)生的激波和膨脹波對模型區(qū)速度場的擾動增加,相對于長軸探管的測量結(jié)果差異就越大,速度場的均勻性越差。在低超聲速時,相同長度孔壁試驗段的消波特性要明顯優(yōu)于槽壁試驗段,因此風(fēng)洞試驗結(jié)果中孔壁試驗段長、短軸探管的模型區(qū)各測點均方根偏差量差異明顯低于槽壁試驗段。

    圖20 孔壁試驗段模型區(qū)均方根偏差對比Fig.20 Comparison of root mean square deviation of Mach number in model area of hole wall test section

    圖21 槽壁試驗段模型區(qū)均方根偏差對比Fig.21 Comparison of root mean square deviation of Mach number in model area of slot wall test section

    5 結(jié)論

    (1)Ma≤0.95 時,仿真結(jié)果表明短軸探管頭部在流場中不會產(chǎn)生激波和膨脹波,頭部在流場中引起的干擾相對較小,從風(fēng)洞試驗結(jié)果來看孔壁試驗段和槽壁試驗段條件下長、短軸探管的測值結(jié)果均未發(fā)生明顯變化,試驗數(shù)據(jù)結(jié)論與數(shù)值仿真結(jié)論一致,該條件下可以采用短軸探管來進行流場校測試驗。

    (2)1.0≤Ma≤1.4 時,仿真結(jié)果中短軸探管頭部在低超聲速來流中會產(chǎn)生逐漸增強的激波和膨脹波,影響區(qū)域和干擾強度明顯大于亞聲速條件下的結(jié)果。風(fēng)洞試驗中隨著馬赫數(shù)的增加,激波角的減小和激波強度的增加,以及洞壁反射波的影響,短軸探管頭部產(chǎn)生的頭波對流場擾動進一步增大,使得孔壁試驗段和槽壁試驗段條件下長、短軸探管的測值結(jié)果差異較亞音速狀態(tài)均明顯擴大。

    (3)1.0≤Ma≤1.4 時,由于孔壁試驗段的消波特性明顯優(yōu)于槽壁試驗段,因此孔壁試驗段條件下長、短軸探管的模型區(qū)各測點馬赫數(shù)差異量均值和均方根偏差量差異明顯低于槽壁試驗段的試驗結(jié)果。

    猜你喜歡
    探管馬赫數(shù)試驗段
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    小口徑γ 能譜深測井探管耐高溫技術(shù)研究
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    跨聲速風(fēng)洞槽壁試驗段流場品質(zhì)提升措施研究
    載荷分布對可控擴散葉型性能的影響
    合肥軌道交通1號線試驗段Ⅰ標基坑監(jiān)測及結(jié)果分析研究
    鄰井隨鉆電磁測距防碰工具模擬試驗研究
    煤田物探測井探管的技術(shù)改進探討
    魅力中國(2017年34期)2017-09-11 10:49:34
    HL-2M真空室試驗段制造工藝技術(shù)
    中國核電(2017年2期)2017-08-11 08:00:54
    能源監(jiān)控管理系統(tǒng)在滬寧高速公路照明試驗段的應(yīng)用
    av在线播放精品| 高清视频免费观看一区二区| 91狼人影院| 亚洲最大成人中文| 中文乱码字字幕精品一区二区三区| 91在线精品国自产拍蜜月| av国产免费在线观看| .国产精品久久| 国产乱人偷精品视频| 晚上一个人看的免费电影| 黄色视频在线播放观看不卡| 日韩 亚洲 欧美在线| 视频区图区小说| 婷婷色综合大香蕉| 草草在线视频免费看| 五月伊人婷婷丁香| 国内揄拍国产精品人妻在线| 欧美精品国产亚洲| av又黄又爽大尺度在线免费看| 欧美日韩视频高清一区二区三区二| 热99国产精品久久久久久7| 欧美极品一区二区三区四区| 亚洲色图av天堂| 免费看a级黄色片| 国产一区二区亚洲精品在线观看| 少妇人妻精品综合一区二区| 国产在线男女| 亚洲精品自拍成人| av又黄又爽大尺度在线免费看| 五月玫瑰六月丁香| 大香蕉97超碰在线| 精品人妻熟女av久视频| 国产亚洲av嫩草精品影院| 日日撸夜夜添| freevideosex欧美| 伦精品一区二区三区| 蜜桃亚洲精品一区二区三区| 亚洲精品影视一区二区三区av| 直男gayav资源| 久久久久久久久久久免费av| 人妻制服诱惑在线中文字幕| 99精国产麻豆久久婷婷| 日韩欧美 国产精品| 久久热精品热| 欧美性感艳星| 69av精品久久久久久| 日本欧美国产在线视频| 天堂俺去俺来也www色官网| 免费黄频网站在线观看国产| 久久久欧美国产精品| 国产白丝娇喘喷水9色精品| 亚洲第一区二区三区不卡| 亚洲真实伦在线观看| 各种免费的搞黄视频| 直男gayav资源| 免费电影在线观看免费观看| 一级毛片久久久久久久久女| 午夜免费鲁丝| 99热网站在线观看| 中文精品一卡2卡3卡4更新| 少妇人妻精品综合一区二区| 一本久久精品| 少妇 在线观看| 97精品久久久久久久久久精品| 色婷婷久久久亚洲欧美| 丝袜美腿在线中文| 免费观看性生交大片5| 美女xxoo啪啪120秒动态图| 嫩草影院精品99| 国产极品天堂在线| 交换朋友夫妻互换小说| freevideosex欧美| 亚洲欧美成人综合另类久久久| 亚州av有码| 色哟哟·www| 国产黄片视频在线免费观看| 一本久久精品| 少妇熟女欧美另类| 午夜精品国产一区二区电影 | 日本一二三区视频观看| 国产黄色视频一区二区在线观看| 欧美xxxx性猛交bbbb| 国产一区有黄有色的免费视频| 日本免费在线观看一区| 久久精品综合一区二区三区| 最近的中文字幕免费完整| 国产高清三级在线| av福利片在线观看| 日本色播在线视频| 国模一区二区三区四区视频| 少妇 在线观看| 欧美97在线视频| 国产女主播在线喷水免费视频网站| 少妇的逼好多水| 精品一区二区三区视频在线| 久久久欧美国产精品| 国产精品偷伦视频观看了| 偷拍熟女少妇极品色| 婷婷色av中文字幕| 国产老妇女一区| 丰满人妻一区二区三区视频av| 亚洲高清免费不卡视频| 亚洲va在线va天堂va国产| 成人特级av手机在线观看| 日韩制服骚丝袜av| 80岁老熟妇乱子伦牲交| 精品视频人人做人人爽| 日韩在线高清观看一区二区三区| 国产毛片a区久久久久| 欧美精品国产亚洲| 日韩欧美 国产精品| 欧美+日韩+精品| 国国产精品蜜臀av免费| 夫妻午夜视频| 久久6这里有精品| 亚洲欧美一区二区三区国产| 国产成人a区在线观看| 黄色配什么色好看| 欧美日韩综合久久久久久| 国产免费福利视频在线观看| 国产一区二区在线观看日韩| 国产亚洲91精品色在线| 国产视频首页在线观看| 婷婷色综合www| 亚洲国产成人一精品久久久| 国产精品麻豆人妻色哟哟久久| 一个人看视频在线观看www免费| 三级男女做爰猛烈吃奶摸视频| 精品人妻偷拍中文字幕| 日韩精品有码人妻一区| 成人特级av手机在线观看| 精品久久久精品久久久| 国产精品久久久久久精品古装| 免费看不卡的av| av.在线天堂| 黑人高潮一二区| 18禁裸乳无遮挡免费网站照片| 午夜精品国产一区二区电影 | 91精品一卡2卡3卡4卡| 寂寞人妻少妇视频99o| 美女主播在线视频| 在线观看三级黄色| 91在线精品国自产拍蜜月| 成年免费大片在线观看| 久久亚洲国产成人精品v| 精品国产露脸久久av麻豆| 26uuu在线亚洲综合色| 国产男女内射视频| 91在线精品国自产拍蜜月| 精品人妻熟女av久视频| 成人一区二区视频在线观看| 国产成人精品福利久久| av国产免费在线观看| 中文精品一卡2卡3卡4更新| 久久精品久久精品一区二区三区| 精品人妻熟女av久视频| 国产成人精品一,二区| 日日摸夜夜添夜夜添av毛片| 伊人久久精品亚洲午夜| 看黄色毛片网站| 色5月婷婷丁香| 欧美+日韩+精品| 高清欧美精品videossex| 久久久色成人| 国产一级毛片在线| 少妇丰满av| 最近最新中文字幕免费大全7| 日本wwww免费看| 国产成人福利小说| 亚洲国产最新在线播放| 一级av片app| 久久这里有精品视频免费| 三级男女做爰猛烈吃奶摸视频| 97人妻精品一区二区三区麻豆| 久久ye,这里只有精品| 国产精品久久久久久精品电影小说 | 久久99蜜桃精品久久| 亚洲精品国产成人久久av| 在线免费十八禁| 老女人水多毛片| 我的老师免费观看完整版| 我要看日韩黄色一级片| 免费黄频网站在线观看国产| 日韩中字成人| 禁无遮挡网站| 日韩,欧美,国产一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本熟妇午夜| 一边亲一边摸免费视频| 成人鲁丝片一二三区免费| 大陆偷拍与自拍| 国产一区二区亚洲精品在线观看| 老师上课跳d突然被开到最大视频| 亚洲av国产av综合av卡| 男插女下体视频免费在线播放| 国产精品人妻久久久久久| 亚洲av成人精品一二三区| 简卡轻食公司| 看免费成人av毛片| 男女国产视频网站| 中文字幕人妻熟人妻熟丝袜美| 日韩av不卡免费在线播放| 亚洲,一卡二卡三卡| 乱码一卡2卡4卡精品| 亚洲精品自拍成人| 禁无遮挡网站| 伦精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 国产精品三级大全| 一本一本综合久久| 丝瓜视频免费看黄片| 亚洲性久久影院| www.av在线官网国产| 国产高潮美女av| 深夜a级毛片| 国产精品爽爽va在线观看网站| 毛片女人毛片| 中文字幕亚洲精品专区| av专区在线播放| 搞女人的毛片| 久久99精品国语久久久| 亚洲一级一片aⅴ在线观看| 免费观看性生交大片5| 亚洲国产色片| 内地一区二区视频在线| 国产伦理片在线播放av一区| av网站免费在线观看视频| 精品久久国产蜜桃| 人妻少妇偷人精品九色| 国产免费又黄又爽又色| 在线看a的网站| 777米奇影视久久| 乱系列少妇在线播放| 99久久精品一区二区三区| 亚洲精品日韩在线中文字幕| 久久久久久久午夜电影| 韩国av在线不卡| 特级一级黄色大片| 深夜a级毛片| 午夜福利视频精品| 真实男女啪啪啪动态图| 三级国产精品欧美在线观看| 色哟哟·www| 日韩一区二区三区影片| 边亲边吃奶的免费视频| 久久久久九九精品影院| 亚洲av二区三区四区| a级毛片免费高清观看在线播放| 亚洲国产精品国产精品| 亚洲成人精品中文字幕电影| 亚洲精品成人av观看孕妇| 成人亚洲欧美一区二区av| 我要看日韩黄色一级片| 成人午夜精彩视频在线观看| 各种免费的搞黄视频| 特大巨黑吊av在线直播| 蜜臀久久99精品久久宅男| 大香蕉97超碰在线| 精品国产露脸久久av麻豆| 天美传媒精品一区二区| 日日摸夜夜添夜夜添av毛片| 777米奇影视久久| 亚洲av不卡在线观看| 99精国产麻豆久久婷婷| 免费观看无遮挡的男女| 国产成人一区二区在线| 一级毛片电影观看| 成人一区二区视频在线观看| 久久久久久伊人网av| 国产精品伦人一区二区| 国产精品国产三级专区第一集| 99久国产av精品国产电影| 男人添女人高潮全过程视频| 成人美女网站在线观看视频| 国产成人免费无遮挡视频| 亚洲四区av| 久久精品久久久久久久性| 黄色配什么色好看| 国产成人免费观看mmmm| 一本一本综合久久| 少妇的逼水好多| 五月伊人婷婷丁香| 亚洲欧美日韩卡通动漫| 免费黄频网站在线观看国产| 99久久精品国产国产毛片| 内地一区二区视频在线| 高清毛片免费看| av福利片在线观看| 一级a做视频免费观看| 国产中年淑女户外野战色| 特级一级黄色大片| 精品国产乱码久久久久久小说| 国产精品久久久久久av不卡| 干丝袜人妻中文字幕| 伦理电影大哥的女人| 免费电影在线观看免费观看| 精品人妻熟女av久视频| 亚洲欧美日韩东京热| 18禁裸乳无遮挡动漫免费视频 | 久久久久国产网址| 亚洲伊人久久精品综合| 色婷婷久久久亚洲欧美| kizo精华| 少妇人妻 视频| 午夜福利视频1000在线观看| 免费看不卡的av| 久久99热这里只有精品18| 国产美女午夜福利| 日本熟妇午夜| 国产精品国产三级国产av玫瑰| 免费在线观看成人毛片| 男的添女的下面高潮视频| 婷婷色综合www| 亚洲成色77777| 久久精品久久精品一区二区三区| 不卡视频在线观看欧美| 在线观看美女被高潮喷水网站| 男人添女人高潮全过程视频| 高清av免费在线| av在线老鸭窝| av免费在线看不卡| 新久久久久国产一级毛片| 久久精品国产自在天天线| 性色av一级| 婷婷色麻豆天堂久久| 中文乱码字字幕精品一区二区三区| 国产av码专区亚洲av| 91精品伊人久久大香线蕉| 一个人看视频在线观看www免费| 青春草亚洲视频在线观看| 成人午夜精彩视频在线观看| 综合色av麻豆| 久久久欧美国产精品| 夫妻午夜视频| 人人妻人人看人人澡| 久久人人爽人人爽人人片va| 成年女人在线观看亚洲视频 | 国产伦理片在线播放av一区| 欧美日韩视频高清一区二区三区二| kizo精华| 国产综合精华液| 人妻一区二区av| 舔av片在线| 免费看不卡的av| a级毛片免费高清观看在线播放| 美女被艹到高潮喷水动态| 一级片'在线观看视频| 亚洲av成人精品一区久久| 精品一区二区免费观看| 国产 一区 欧美 日韩| 九九久久精品国产亚洲av麻豆| 欧美日韩亚洲高清精品| 精品国产三级普通话版| 欧美成人午夜免费资源| 亚洲国产精品专区欧美| 老女人水多毛片| 日韩大片免费观看网站| a级毛片免费高清观看在线播放| 久久精品国产亚洲av天美| 日韩中字成人| 国产久久久一区二区三区| 26uuu在线亚洲综合色| 在线亚洲精品国产二区图片欧美 | 伊人久久国产一区二区| 我的女老师完整版在线观看| 在线观看一区二区三区激情| xxx大片免费视频| 天堂网av新在线| 97超视频在线观看视频| 亚洲在线观看片| 大码成人一级视频| 青春草视频在线免费观看| 大码成人一级视频| 久久这里有精品视频免费| 亚洲欧洲日产国产| 99久久精品一区二区三区| xxx大片免费视频| 亚洲三级黄色毛片| 国产精品一及| 亚洲三级黄色毛片| 一级片'在线观看视频| 新久久久久国产一级毛片| 国产成人免费观看mmmm| 夜夜爽夜夜爽视频| 激情五月婷婷亚洲| 欧美激情国产日韩精品一区| 69人妻影院| 99久久九九国产精品国产免费| 色播亚洲综合网| 国产精品一区二区三区四区免费观看| 男男h啪啪无遮挡| 一个人看的www免费观看视频| 亚洲精品影视一区二区三区av| av天堂中文字幕网| 热99国产精品久久久久久7| 成人毛片60女人毛片免费| 在线 av 中文字幕| freevideosex欧美| 极品教师在线视频| 日韩不卡一区二区三区视频在线| 最近最新中文字幕大全电影3| 色吧在线观看| 欧美97在线视频| 色哟哟·www| 免费av观看视频| 欧美性感艳星| 久久久久久久久久久免费av| 亚洲欧美清纯卡通| 在线看a的网站| 国产成人aa在线观看| 汤姆久久久久久久影院中文字幕| 国产精品久久久久久久久免| 亚洲欧洲国产日韩| 最近最新中文字幕大全电影3| 在线观看三级黄色| 国产在线男女| 免费大片18禁| 五月开心婷婷网| 亚洲aⅴ乱码一区二区在线播放| 黄色日韩在线| 亚洲四区av| 国产乱人视频| 日本黄色片子视频| 2022亚洲国产成人精品| 欧美丝袜亚洲另类| 深夜a级毛片| 熟妇人妻不卡中文字幕| 久久99热这里只有精品18| 日本三级黄在线观看| 国产成人免费观看mmmm| 一级毛片aaaaaa免费看小| 国国产精品蜜臀av免费| 国产毛片在线视频| 国产免费视频播放在线视频| 成年女人看的毛片在线观看| 精品一区二区三区视频在线| 国产精品成人在线| 亚洲av在线观看美女高潮| 男女下面进入的视频免费午夜| 国产亚洲av片在线观看秒播厂| 最近最新中文字幕免费大全7| 精品亚洲乱码少妇综合久久| 国产亚洲91精品色在线| 国产伦理片在线播放av一区| 69av精品久久久久久| 最近中文字幕2019免费版| av国产久精品久网站免费入址| 大话2 男鬼变身卡| 一级黄片播放器| 精品人妻视频免费看| 别揉我奶头 嗯啊视频| 成人国产av品久久久| 亚洲av在线观看美女高潮| av天堂中文字幕网| 有码 亚洲区| 日韩精品有码人妻一区| 老司机影院成人| 亚洲国产精品999| 亚洲一级一片aⅴ在线观看| 国产精品蜜桃在线观看| 亚洲精品,欧美精品| 黄色视频在线播放观看不卡| 寂寞人妻少妇视频99o| 亚洲va在线va天堂va国产| 人人妻人人澡人人爽人人夜夜| 69人妻影院| 伦精品一区二区三区| 嫩草影院精品99| 波多野结衣巨乳人妻| 2018国产大陆天天弄谢| 最近最新中文字幕大全电影3| 亚洲国产精品专区欧美| 熟妇人妻不卡中文字幕| 日本午夜av视频| 尾随美女入室| 欧美日韩综合久久久久久| 自拍偷自拍亚洲精品老妇| 亚洲精品国产av成人精品| 麻豆成人av视频| 在线观看三级黄色| 在线观看一区二区三区激情| 亚洲国产欧美在线一区| 国产一区有黄有色的免费视频| 九九久久精品国产亚洲av麻豆| 精品99又大又爽又粗少妇毛片| 哪个播放器可以免费观看大片| 人妻少妇偷人精品九色| 成人亚洲欧美一区二区av| 性插视频无遮挡在线免费观看| 日日摸夜夜添夜夜添av毛片| 欧美另类一区| 成人欧美大片| 国产亚洲午夜精品一区二区久久 | 2022亚洲国产成人精品| 午夜爱爱视频在线播放| 老司机影院毛片| 国产av国产精品国产| 91狼人影院| 亚洲最大成人av| 久久精品熟女亚洲av麻豆精品| 亚洲自偷自拍三级| 国产精品久久久久久久久免| 久久精品人妻少妇| 国产淫片久久久久久久久| h日本视频在线播放| 久久久久久久精品精品| 九九爱精品视频在线观看| 亚洲一区二区三区欧美精品 | 欧美日韩综合久久久久久| 卡戴珊不雅视频在线播放| 国产精品国产av在线观看| 日韩三级伦理在线观看| 国产亚洲一区二区精品| 神马国产精品三级电影在线观看| 春色校园在线视频观看| 成人高潮视频无遮挡免费网站| 精品午夜福利在线看| 久久久a久久爽久久v久久| 日韩伦理黄色片| 色视频www国产| 亚洲激情五月婷婷啪啪| 极品教师在线视频| 国产欧美另类精品又又久久亚洲欧美| 午夜免费男女啪啪视频观看| 18禁在线无遮挡免费观看视频| 国产成人精品久久久久久| 成年免费大片在线观看| 七月丁香在线播放| 午夜免费观看性视频| 国产免费视频播放在线视频| 亚洲欧美精品自产自拍| 极品教师在线视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人国产麻豆网| 成人美女网站在线观看视频| 欧美日韩视频高清一区二区三区二| 欧美日韩一区二区视频在线观看视频在线 | 国产大屁股一区二区在线视频| 一级毛片电影观看| 网址你懂的国产日韩在线| 亚洲欧美中文字幕日韩二区| 日韩中字成人| 99热这里只有是精品50| 视频中文字幕在线观看| 人体艺术视频欧美日本| 人妻系列 视频| 亚洲av成人精品一二三区| 小蜜桃在线观看免费完整版高清| 亚洲色图综合在线观看| 午夜福利高清视频| 99热这里只有是精品在线观看| 午夜福利在线在线| 久热久热在线精品观看| 亚洲国产精品专区欧美| 日韩成人伦理影院| 搡老乐熟女国产| 特级一级黄色大片| 能在线免费看毛片的网站| 交换朋友夫妻互换小说| 亚洲图色成人| 97在线视频观看| 久久久久久久久久久免费av| 亚洲欧美中文字幕日韩二区| 在线观看一区二区三区| 日韩一区二区视频免费看| 少妇猛男粗大的猛烈进出视频 | 晚上一个人看的免费电影| 国产成人精品婷婷| 人妻一区二区av| 日韩亚洲欧美综合| 国产精品偷伦视频观看了| 国产 一区 欧美 日韩| 一级二级三级毛片免费看| 你懂的网址亚洲精品在线观看| 国产一级毛片在线| 一级片'在线观看视频| 成人亚洲精品一区在线观看 | 内地一区二区视频在线| 成人特级av手机在线观看| 看十八女毛片水多多多| 国产乱来视频区| 亚洲av男天堂| 免费看光身美女| 男人和女人高潮做爰伦理| 欧美bdsm另类| 国产欧美日韩精品一区二区| 欧美成人a在线观看| 午夜福利在线在线| 国产成人a区在线观看| 免费观看性生交大片5| 久久久亚洲精品成人影院| 九九爱精品视频在线观看| 国产亚洲最大av| 在线天堂最新版资源| 国产亚洲一区二区精品| 亚洲自偷自拍三级| 人人妻人人爽人人添夜夜欢视频 | 麻豆国产97在线/欧美| 亚洲色图av天堂| 99热这里只有是精品在线观看| 黑人高潮一二区| 看十八女毛片水多多多| 久久人人爽人人片av| 国产精品国产三级国产av玫瑰| 日韩国内少妇激情av| 欧美xxxx性猛交bbbb| 国产精品一及| 黄色配什么色好看| 狂野欧美白嫩少妇大欣赏| 七月丁香在线播放| 一个人看视频在线观看www免费| av女优亚洲男人天堂| 日韩精品有码人妻一区| kizo精华| 国产精品熟女久久久久浪|