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

    亞聲速條件下總壓探針臨壁效應(yīng)的數(shù)值研究

    2020-05-06 06:25:28鄭光華王肖郭林惠呈程呂文召
    航空工程進(jìn)展 2020年2期
    關(guān)鍵詞:支桿總壓測量誤差

    鄭光華,王肖,郭林,惠呈程,呂文召

    (1.西北工業(yè)大學(xué) 動(dòng)力與能源學(xué)院,西安 710129) (2.中國航天空氣動(dòng)力技術(shù)研究院 特種飛行器總體技術(shù)設(shè)計(jì)部,北京 100074)

    0 引 言

    探針因其結(jié)構(gòu)簡單、可靠性高被廣泛應(yīng)用于流場測量,但在實(shí)際測量時(shí),被測流場中近壁效應(yīng)產(chǎn)生的影響難以在探針標(biāo)定時(shí)消除[1]。研究發(fā)現(xiàn),探針在測量過程中,不可避免地引起局部流場的擾動(dòng),不僅使壓力的測量值與真實(shí)值產(chǎn)生偏差[2],還會(huì)造成流動(dòng)損失,降低機(jī)械效率[3]。當(dāng)測量通道太窄時(shí),為了避免臨壁效應(yīng)對(duì)測量準(zhǔn)確性的影響,同時(shí)減小流場內(nèi)的流動(dòng)損失,確定總壓探針在流場中的影響規(guī)律就顯得尤為重要。

    J.Lepicovsky[4]將探針安裝在轉(zhuǎn)子上,研究其產(chǎn)生的尾流和流動(dòng)阻塞對(duì)流場的影響,發(fā)現(xiàn)將探針安裝在轉(zhuǎn)子中會(huì)導(dǎo)致轉(zhuǎn)子出口流動(dòng)惡化,靜壓系數(shù)降低,且葉柵通道下部阻塞增加使得通道上部未阻擋部分的軸向速度分量增加,勢必會(huì)影響壓氣機(jī)的穩(wěn)定裕度;向宏輝等[5]采用三維模擬方法分析了壓氣機(jī)在安裝探針后的性能及其內(nèi)部流場的變化,發(fā)現(xiàn)級(jí)間局部靜葉安裝探針后,壓氣機(jī)下堵點(diǎn)流量減小0.1%,最高效率降低0.4%,失速點(diǎn)流量增大0.15%,且探針的擾流加劇還可能導(dǎo)致壓氣機(jī)氣動(dòng)失穩(wěn);孫志強(qiáng)等[6]針對(duì)不同型號(hào)的皮托管開展試驗(yàn)以研究皮托管的支桿和安裝角對(duì)測量結(jié)果的影響,結(jié)果表明,在低速測量時(shí),皮托管支桿直徑對(duì)總壓的測量影響較小,安裝角在±10°以內(nèi)速度測量誤差小于±2%;林其勛等[7]、馬宏偉等[8]對(duì)總壓探針在測量中臨壁效應(yīng)的影響進(jìn)行了試驗(yàn)研究,結(jié)果表明,可以用疊加點(diǎn)渦的方法來描述其臨壁效應(yīng),且點(diǎn)渦強(qiáng)度隨著近壁距離的減小而增強(qiáng),近壁效應(yīng)一般對(duì)靜壓測量影響較大,對(duì)總壓測量影響相對(duì)較?。籋.Bubeck等[9-10]針對(duì)楔形四孔探針和球頭五孔探針,研究探針和壁面的相互作用對(duì)測量誤差的影響,發(fā)現(xiàn)近壁測量時(shí)探針的俯仰角測量值存在誤差,前者偏離10°,后者偏離5°。

    上述文獻(xiàn)多是研究探針對(duì)流場或來流參數(shù)對(duì)單一探針測量結(jié)果的影響。而本文考慮流道壁面與探針的相互作用可能改變探針測量的準(zhǔn)確性,且對(duì)比分析四種不同頭部結(jié)構(gòu)的總壓探針(A型、B型、C型、D型)在不同近壁距離的條件下對(duì)機(jī)匣近壁流場及自身總壓測量值造成的影響差異,并利用文獻(xiàn)實(shí)驗(yàn)結(jié)果驗(yàn)證數(shù)值模擬結(jié)果。

    1 計(jì)算模型和方法

    1.1 計(jì)算模型

    計(jì)算域選擇圓柱形,為了減小來流擾動(dòng)的影響,計(jì)算域入口到總壓探針前端的距離為40 mm,總壓探針末端到計(jì)算域出口的距離設(shè)定為100 mm。考慮中小型發(fā)動(dòng)機(jī)實(shí)際尺寸[11],計(jì)算域的外徑為150 mm,總壓探針計(jì)算模型直徑最大為4 mm,堵塞比遠(yuǎn)遠(yuǎn)小于0.5%,堵塞效應(yīng)造成的影響可忽略不計(jì),后經(jīng)數(shù)值模擬驗(yàn)證,計(jì)算域尺寸選擇合理。計(jì)算域坐標(biāo)統(tǒng)一定義為:沿y軸正方向?yàn)闅饬鞣较?,z軸為計(jì)算域豎直高度方向,x軸為計(jì)算域水平徑向方向,總壓探針入口所在截面中心位置為y=0處。計(jì)算模型如圖1所示。

    圖1 計(jì)算模型Fig.1 Calculation model

    研究近機(jī)匣臨壁效應(yīng)的總壓探針計(jì)算模型主要有四種:A型、B型、C型、D型,如圖2所示。A型為具有圓柱形外形及內(nèi)錐角內(nèi)腔的單點(diǎn)總壓探針,內(nèi)錐角β=25°,外徑d=2 mm,內(nèi)徑d′=0.5 mm;B型、C型總壓探針無內(nèi)錐角且為平頭圓柱形總壓探針,外徑d都為2 mm,內(nèi)徑d′分別為0.3和0.5 mm;D型總壓探針具有半球形頭部結(jié)構(gòu),外徑d=2 mm,內(nèi)徑d′=0.5 mm。

    (a) A型總壓探針

    (b) B/C型總壓探針

    (c) D型總壓探針圖2 總壓探針結(jié)構(gòu)Fig.2 Structures of total pressure probe

    1.2 網(wǎng)格劃分

    利用ANSYS ICEM對(duì)計(jì)算模型進(jìn)行進(jìn)一步處理并劃分網(wǎng)格,整體采用非結(jié)構(gòu)體網(wǎng)格,總壓探針近壁面及計(jì)算域壁面處采用邊界層網(wǎng)格,保證近壁端湍流模型所需的y+為1左右,并且在近壁區(qū)內(nèi)設(shè)定10層邊界層網(wǎng)格。網(wǎng)格數(shù)量達(dá)到400萬以后,總壓測量值基本穩(wěn)定,總壓測量相對(duì)誤差遠(yuǎn)小于1%,可認(rèn)為網(wǎng)格數(shù)量對(duì)計(jì)算結(jié)果的影響可忽略,故將數(shù)量大于400萬范圍的網(wǎng)格作為本文計(jì)算網(wǎng)格。網(wǎng)格劃分及邊界層處理如圖3所示。

    (a) 整體計(jì)算域網(wǎng)格

    (b) 邊界層處理

    (c) 邊界層處網(wǎng)格放大圖3 計(jì)算網(wǎng)格Fig.3 Calculation grid

    1.3 邊界條件

    計(jì)算域入口邊界條件設(shè)定為壓力入口:設(shè)置總溫、總壓、速度方向;環(huán)境溫度、總溫、總靜壓與馬赫數(shù)的關(guān)系如式(1)~式(2)所示[12],其中環(huán)境溫度為293.15 K、靜壓為101 325 Pa。

    (1)

    (2)

    式中:Tt為總溫設(shè)定值;Pt為總壓設(shè)定值;Tfree為自由來流環(huán)境溫度;Pfree為自由來流靜壓;Mafree為自由來流馬赫數(shù);k為氣體的比熱比。

    計(jì)算域出口邊界條件設(shè)定為壓力出口:出口靜壓為101 325 Pa,計(jì)算域周向壁面與總壓探針壁面給定無滑移壁面,采用FLUENT求解,湍流模型為standardk-ε模型;馬赫數(shù)為0.4,即為亞聲速流動(dòng)。

    1.4 參數(shù)定義

    總壓測量誤差:總壓測量值與理論總壓值的差值,再與理論總壓值的比值。其計(jì)算公式為

    (3)

    總壓損失系數(shù):不同徑向高度下當(dāng)?shù)乜倝褐蹬c來流總壓之差占動(dòng)壓頭的百分比。其計(jì)算公式為

    (4)

    速度均勻度:基于面積加權(quán)平均速度和質(zhì)量加權(quán)平均速度的流動(dòng)均勻性指標(biāo)[13-14]。計(jì)算公式為

    (5)

    (6)

    (7)

    式中:Vn為面積加權(quán)平均速度;Vm為質(zhì)量加權(quán)平均速度;A為總表面積;Aj為第j個(gè)單元面上的面積矢量;vj為第j個(gè)單元面上的速度矢量;ρj為第j個(gè)單元面上的流體密度;n為同流截面被劃分成的單元面數(shù)。

    2 結(jié)果與討論

    臨壁效應(yīng)不同于總壓探針進(jìn)入附面層時(shí)由于總壓探針空間分辨率不夠而引起的誤差,二者概念不同。相關(guān)實(shí)驗(yàn)[7]表明,當(dāng)總壓探針的布置位置未達(dá)到附面層厚度時(shí),臨壁效應(yīng)的影響就已經(jīng)較大。當(dāng)總壓探針進(jìn)入附面層后,臨壁效應(yīng)與附面層影響同時(shí)起作用,使得總壓探針測量誤差更大。

    已知,從物體表面速度為0處開始,沿法線方向至速度等于來流速度的99%位置處的距離為附面層厚度,Ma=0.4時(shí)主流速度約為137.6 m/s,在z=74.1 mm處速度達(dá)到主流速度,根據(jù)差值法,得出99%速度處所對(duì)應(yīng)的位置為z=73.7 mm,而速度為0處壁面位置為z=75 mm,則附面層厚度為1.3 mm。

    為了保證臨壁效應(yīng)為獨(dú)立的影響因素,應(yīng)使得總壓探針完全處于附面層之外,因此,當(dāng)近壁距離c(計(jì)算域壁面到總壓探針中心軸線的距離)大于2.3 mm,即c/d>1.15時(shí),總壓探針完全處于附面層之外,使臨壁效應(yīng)在很廣的c范圍內(nèi)成為獨(dú)立的因素。

    為了驗(yàn)證數(shù)值計(jì)算的可靠性,本文針對(duì)平頭圓柱形單孔總壓探針(B型)開展不同進(jìn)氣偏轉(zhuǎn)角(φ)的數(shù)值計(jì)算,并與試驗(yàn)結(jié)果[15]進(jìn)行對(duì)比,如圖4所示。

    圖4 結(jié)果對(duì)比Fig.4 Comparison of results

    從圖4可以看出:數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果在小角度內(nèi)重合,大角度稍有偏離,總體上基本符合。本文研究內(nèi)容都為0°進(jìn)氣攻角,可以保證數(shù)值計(jì)算的可靠性。

    2.1 臨壁效應(yīng)對(duì)測量結(jié)果的影響

    總壓測量誤差隨相對(duì)近壁距離c/d的變化曲線如圖5所示,近壁距離最小取2.5 mm,即c/d=1.25時(shí),總壓探針相當(dāng)靠近壁面但還未達(dá)到附面層邊界。

    圖5 總壓測量誤差隨相對(duì)近壁距離的變化曲線Fig.5 Curve of total pressure measurement error with relative near wall distance

    從圖5可以看出:四種總壓探針的總壓測量誤差隨著相對(duì)近壁距離的減小整體向負(fù)半軸偏轉(zhuǎn)且斜率逐漸增大,總壓測量誤差逐漸增大;當(dāng)總壓探針臨壁時(shí)使得氣流方向改變,向總壓探針下側(cè)面偏轉(zhuǎn),致使總壓探針產(chǎn)生總壓測量誤差;當(dāng)相對(duì)近壁距離c/d=1.25時(shí),最大總壓測量誤差小于1%;隨著近壁距離的增大,總壓測量誤差越小,當(dāng)c/d>3時(shí),總壓測量誤差曲線趨近于0,說明此時(shí)的臨壁效應(yīng)已消失。

    圖5中四種總壓探針的變化曲線整體規(guī)律一致,但是曲線間存在差異,主要是由于總壓探針頭部形狀不同引起的。從圖5可以看出:A型總壓探針整體總壓測量誤差最小,D型總壓探針整體總壓測量誤差最大,B型和C型總壓探針整體總壓測量誤差大小介于最大值和最小值之間,且B型總壓探針總壓測量誤差大于C型總壓探針,表明孔口帶倒角的總壓探針受到臨壁效應(yīng)的影響最小,半球形頭部結(jié)構(gòu)的總壓探針受到臨壁效應(yīng)的影響最大,兩者總壓測量誤差最大差值為0.45%。

    2.2 臨壁效應(yīng)對(duì)總壓探針附近流場繞流變化的影響

    總壓探針與壁面的距離越小,其附近流場的氣流流線受到的干擾將逐漸增大,現(xiàn)以外徑2 mm、內(nèi)徑0.5 mm的C型平頭圓柱形總壓探針為例進(jìn)行分析。近壁距離c分別為3、10 mm時(shí)總壓探針附近流場繞流情況及速度分布如圖6所示,速度范圍為0~152 m/s。

    從速度分布可以看出:總壓探針入口處速度逐漸降低并在測壓通道內(nèi)滯止,速度為0;總壓探針前端外壁面兩側(cè)產(chǎn)生低速區(qū);總壓探針支桿背部速度明顯低于主流速度,且越靠近機(jī)匣壁面速度越低,并且隨著總壓探針近壁距離越大,支桿越長,支桿背部低速區(qū)越寬。

    (a) 近壁距離為3 mm

    (b) 近壁距離為10 mm圖6 總壓探針附近繞流情況及速度分布Fig.6 Flow around the total pressure probe and velocity distribution

    從總壓探針附近流線分布可以看出:主流流線在遇到總壓探針頭部時(shí)向其下表面發(fā)生偏轉(zhuǎn)并向下游流去,在到達(dá)總壓探針支桿背部處流線向機(jī)匣壁面發(fā)生偏轉(zhuǎn),近壁距離越大支桿背部越容易產(chǎn)生較大旋渦;總壓探針前端壁面由于氣流發(fā)生偏轉(zhuǎn)在其表面產(chǎn)生邊界層分離并形成旋渦,當(dāng)靠近機(jī)匣壁面較近時(shí),總壓探針前端上下表面產(chǎn)生不對(duì)稱旋渦,近壁距離c=3 mm時(shí)總壓探針下表面邊界層厚度明顯大于上表面,上側(cè)由于氣流受到機(jī)匣壁面影響邊界層較?。贿h(yuǎn)離壁面時(shí),近壁距離c=10 mm時(shí)臨壁效應(yīng)不起作用,總壓探針前端兩側(cè)產(chǎn)生近似對(duì)稱旋渦。

    分析可知,臨壁效應(yīng)主要影響的是近壁面區(qū)域流場分布以及近壁面的測量值,在實(shí)際測量中不排除近壁測量的極端情況,故對(duì)臨壁效應(yīng)的研究十分有意義。

    2.3 臨壁效應(yīng)對(duì)被測流場總壓分布的影響

    近壁距離c=3 mm和c=10 mm時(shí)x=0截面的C型總壓探針附近總壓分布云圖如圖7所示。

    (a) 近壁距離為3 mm

    (b) 近壁距離為10 mm圖7 總壓探針附近總壓分布Fig.7 Total pressure distribution near the total pressure probe

    從圖7可以看出:總壓探針上游總壓分布較均勻,其對(duì)上游的總壓分布影響較小,總壓探針表面及下游分別出現(xiàn)低壓區(qū)及尾跡區(qū),其對(duì)下游的總壓分布影響較大;隨著近壁距離增大總壓探針插入深度增大,尾跡沿z方向的寬度也相應(yīng)增大。

    為進(jìn)一步對(duì)比分析不同近壁距離時(shí)的總壓探針尾跡對(duì)下游的干擾情況,給出距離總壓探針末端30 mm位置處的尾跡區(qū)總壓損失系數(shù)沿計(jì)算域z軸(高度)方向的分布曲線,如圖8所示,橫坐標(biāo)表示總壓損失系數(shù),縱坐標(biāo)表示z軸高度,圖例表示相對(duì)近壁距離c/d只取1.5~5.5范圍。

    圖8 沿高度方向總壓損失系數(shù)分布Fig.8 Distribution of total pressure loss coefficient along the height direction

    從圖8可以看出:z<55 mm時(shí),總壓損失系數(shù)曲線幾乎與0重合,說明總壓探針的尾跡對(duì)計(jì)算域中心區(qū)域的總壓分布幾乎無影響;當(dāng)z>55 mm

    時(shí),隨著總壓探針近壁距離的不同,總壓損失系數(shù)曲線呈現(xiàn)一定規(guī)律的變化趨勢,相對(duì)近壁距離c/d越大總壓損失曲線與0軸的分離點(diǎn)(圖中畫圈標(biāo)明處)越早,說明總壓探針與機(jī)匣近壁距離越大,支桿插入深度越大,沿z軸方向的尾跡寬度就越大,低壓區(qū)造成的總壓損失范圍越大;隨著近壁距離的增大,支桿長度越長,支桿與機(jī)匣壁面接觸處的總壓損失強(qiáng)度越大,c/d為1.5~5.5時(shí),總壓損失強(qiáng)度從50%增大到70%;隨著近壁距離增大,在總壓損失系數(shù)達(dá)到10%處曲線出現(xiàn)轉(zhuǎn)折,主要原因是總壓探針支桿長度增加,背部低壓區(qū)分布較均勻,總壓損失分布也相應(yīng)較均勻,總壓損失系數(shù)較接近。

    2.4 臨壁效應(yīng)對(duì)流場流動(dòng)均勻性的影響

    計(jì)算域各截面示意圖如圖9所示,總壓探針上游及下游各截取兩個(gè)截面,總壓探針測量段及支桿處各截取一個(gè)截面,根據(jù)截面所在位置分別用i1、i2、i3、i4、i5、i6表示,如表1所示。

    圖9 計(jì)算域截面示意圖Fig.9 Schematic diagram of the section in calculation domain 表1 各截面坐標(biāo)定義 Table 1 Definition of each section coordinate

    截面名稱截面位置/mm截面名稱截面位置/mmi1y=-20i4y=13i2y=-5i5y=20i3y=5i6y=60

    被測流場的流動(dòng)均勻性變化如圖10所示,橫坐標(biāo)表示各截面位置,縱坐標(biāo)表示速度均勻度。

    圖10 被測流場流動(dòng)均勻性曲線Fig.10 Flow field uniformity curve of the measured flow field

    從圖10可以看出:從計(jì)算域進(jìn)口至出口流動(dòng)均勻性整體下降,說明被測流場內(nèi)速度存在擾動(dòng);由于速度在到達(dá)總壓探針入口截面時(shí)逐漸降低,從計(jì)算域進(jìn)口至總壓探針前端速度均勻度曲線呈下降趨勢;總壓探針附近出現(xiàn)低壓區(qū)且靠近總壓探針支桿后方處出現(xiàn)速度最小值,總壓探針附近截面(如i4~i5)速度均勻度最低;隨著氣流向下游流動(dòng),尾跡損失強(qiáng)度逐漸減小,直到沿程損失占主導(dǎo)作用,故下游速度均勻度曲線先增大后減小;隨著總壓探針與壁面相對(duì)距離由5.5到1.5,兩壁面間相互干擾的程度越大,總壓探針附近處截面的速度均勻度減小0.035%。

    3 結(jié) 論

    (1) 為了研究臨壁效應(yīng)單獨(dú)對(duì)總壓探針總壓測量誤差的影響,近壁距離須大于或近似等于附面層厚度,使得總壓探針不受附面層效應(yīng)的影響。已知附面層厚度為1.3 mm,此時(shí)近壁距離越小,氣流偏轉(zhuǎn)角越大,臨壁效應(yīng)造成的總壓測量誤差越大,最小相對(duì)近壁距離c/d=1.25時(shí)總壓測量誤差最大,最大值小于1%。為使單點(diǎn)總壓探針獲得較準(zhǔn)確的總壓測量值,布置時(shí)近壁距離應(yīng)大于等于3倍支桿外徑。

    (2) 當(dāng)總壓探針頭部形狀不同時(shí),臨壁效應(yīng)造成的總壓測量誤差不同,對(duì)半球形頭部結(jié)構(gòu)的總壓探針影響最大,對(duì)帶半錐角型總壓探針影響最小,且隨著總壓探針支桿直徑d越大影響越大。

    (3) 當(dāng)c/d減小時(shí),總壓探針前端氣流偏轉(zhuǎn)角產(chǎn)生誤差,總壓探針上表面邊界層厚度減小、下表面邊界層厚度增大;c/d為1.5~5.5,隨著支桿插入深度越深,沿高度方向的尾跡區(qū)寬度增加,總壓探針末端30 mm處尾跡總壓損失系數(shù)最大值從50%增大到70%。

    (4) 被測流場流動(dòng)均勻性整體減小,靠近總壓探針后方i5截面處速度曲線下降幅度最大,c/d從5.5到1.5時(shí),速度均勻度減小0.035%。

    猜你喜歡
    支桿總壓測量誤差
    總壓探針性能結(jié)構(gòu)敏感性分析
    單支桿腹撐支架干擾測力及PIV試驗(yàn)研究
    密度測量誤差分析
    航空發(fā)動(dòng)機(jī)進(jìn)氣總壓畸變地面試驗(yàn)數(shù)據(jù)處理方法綜述
    可調(diào)式總壓耙設(shè)計(jì)及應(yīng)用
    縱向數(shù)據(jù)下變系數(shù)測量誤差模型的漸近估計(jì)
    2 m超聲速風(fēng)洞流場變速壓控制方法研究
    新型移動(dòng)支撐架的制作與應(yīng)用
    教學(xué)用多功能支架
    教學(xué)用多功能支架
    国产片内射在线| 无限看片的www在线观看| 久热爱精品视频在线9| 日本黄色视频三级网站网址| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av美国av| 97人妻精品一区二区三区麻豆 | 欧美+亚洲+日韩+国产| 美女大奶头视频| 国产精品一区二区三区四区久久 | 欧美绝顶高潮抽搐喷水| 国产一区二区三区视频了| 国产精品久久久人人做人人爽| 69av精品久久久久久| 免费看日本二区| 国产亚洲精品av在线| 老熟妇仑乱视频hdxx| 国产欧美日韩一区二区三| 啪啪无遮挡十八禁网站| 精品久久蜜臀av无| 国产亚洲精品综合一区在线观看 | 亚洲熟妇熟女久久| 久久九九热精品免费| 国产精品永久免费网站| 亚洲午夜精品一区,二区,三区| 亚洲欧美日韩无卡精品| 听说在线观看完整版免费高清| 99精品久久久久人妻精品| 中文在线观看免费www的网站 | 国产欧美日韩一区二区三| a级毛片在线看网站| 色综合欧美亚洲国产小说| 女人爽到高潮嗷嗷叫在线视频| 波多野结衣高清无吗| 国产成人精品无人区| 欧美日韩中文字幕国产精品一区二区三区| 欧美黄色淫秽网站| 99久久国产精品久久久| 最新在线观看一区二区三区| 国产久久久一区二区三区| 亚洲精品粉嫩美女一区| 久久人人精品亚洲av| 12—13女人毛片做爰片一| 国产男靠女视频免费网站| 少妇被粗大的猛进出69影院| 美女高潮到喷水免费观看| 88av欧美| 三级毛片av免费| 波多野结衣巨乳人妻| 99国产精品一区二区三区| 看黄色毛片网站| 99精品欧美一区二区三区四区| 亚洲国产欧洲综合997久久, | 午夜激情福利司机影院| 午夜久久久在线观看| 国产99久久九九免费精品| 两性午夜刺激爽爽歪歪视频在线观看 | 看免费av毛片| 亚洲精品一区av在线观看| av在线天堂中文字幕| 免费在线观看亚洲国产| 色播在线永久视频| 国产亚洲精品第一综合不卡| 久久久久久久精品吃奶| 午夜成年电影在线免费观看| 国产一级毛片七仙女欲春2 | 国产精品99久久99久久久不卡| 国产又爽黄色视频| av有码第一页| 免费看十八禁软件| 在线av久久热| 国产激情偷乱视频一区二区| 亚洲五月色婷婷综合| 天天躁夜夜躁狠狠躁躁| av欧美777| 久久久久久久午夜电影| 久久久久国内视频| 精品欧美一区二区三区在线| 中文亚洲av片在线观看爽| 男女下面进入的视频免费午夜 | 两个人看的免费小视频| 国产97色在线日韩免费| √禁漫天堂资源中文www| 波多野结衣巨乳人妻| 午夜两性在线视频| 欧美黑人精品巨大| 国产成人影院久久av| 久久人妻av系列| 一级a爱视频在线免费观看| 午夜福利在线在线| 婷婷六月久久综合丁香| 老汉色av国产亚洲站长工具| 欧美性猛交╳xxx乱大交人| a在线观看视频网站| 狂野欧美激情性xxxx| 黄片播放在线免费| 身体一侧抽搐| 女生性感内裤真人,穿戴方法视频| 亚洲成人久久爱视频| 看黄色毛片网站| 无遮挡黄片免费观看| 亚洲av第一区精品v没综合| 亚洲av电影不卡..在线观看| 黄网站色视频无遮挡免费观看| 欧美黑人欧美精品刺激| 女警被强在线播放| 在线观看舔阴道视频| 亚洲欧美日韩无卡精品| 久久亚洲真实| 日日夜夜操网爽| 亚洲专区中文字幕在线| 岛国视频午夜一区免费看| 欧美日韩中文字幕国产精品一区二区三区| 欧美绝顶高潮抽搐喷水| 岛国在线观看网站| 亚洲一区二区三区色噜噜| 国产精品久久久久久亚洲av鲁大| 老鸭窝网址在线观看| 制服人妻中文乱码| ponron亚洲| 久久久精品欧美日韩精品| 精品欧美一区二区三区在线| 精品少妇一区二区三区视频日本电影| 女性生殖器流出的白浆| 99热这里只有精品一区 | 国产日本99.免费观看| 中文字幕精品免费在线观看视频| 久热这里只有精品99| 久久午夜亚洲精品久久| 亚洲av中文字字幕乱码综合 | 最新美女视频免费是黄的| 国产一区二区三区视频了| 亚洲无线在线观看| 精品午夜福利视频在线观看一区| 亚洲成人久久性| 黑丝袜美女国产一区| 成熟少妇高潮喷水视频| 日韩中文字幕欧美一区二区| 俄罗斯特黄特色一大片| 亚洲精品国产一区二区精华液| 欧美三级亚洲精品| 91成年电影在线观看| www国产在线视频色| 黄色成人免费大全| 男女床上黄色一级片免费看| 亚洲欧美精品综合久久99| 又紧又爽又黄一区二区| 哪里可以看免费的av片| 欧美黄色淫秽网站| 女同久久另类99精品国产91| 老熟妇仑乱视频hdxx| 国产不卡一卡二| 亚洲精品国产精品久久久不卡| 99国产精品一区二区蜜桃av| 欧美一级毛片孕妇| 国产亚洲精品久久久久久毛片| 黄色视频不卡| 久久久国产成人精品二区| 国产精品99久久99久久久不卡| 波多野结衣高清作品| 免费在线观看黄色视频的| 日韩国内少妇激情av| 欧美黑人精品巨大| av福利片在线| 超碰成人久久| 国产av不卡久久| 亚洲免费av在线视频| 亚洲专区国产一区二区| 久久香蕉国产精品| 久久草成人影院| 亚洲成人免费电影在线观看| 天天添夜夜摸| 91老司机精品| 亚洲成人精品中文字幕电影| 午夜亚洲福利在线播放| 亚洲美女黄片视频| 成人av一区二区三区在线看| 亚洲激情在线av| 久久中文字幕人妻熟女| aaaaa片日本免费| 午夜激情av网站| 黄色视频不卡| 久久午夜综合久久蜜桃| 精品人妻1区二区| 夜夜爽天天搞| 在线观看舔阴道视频| 国产蜜桃级精品一区二区三区| 两个人视频免费观看高清| 18禁黄网站禁片午夜丰满| 国产av一区在线观看免费| 精品福利观看| 在线视频色国产色| 久9热在线精品视频| 中文字幕精品免费在线观看视频| 啦啦啦 在线观看视频| 精品久久蜜臀av无| 国产一区二区三区视频了| 中文字幕另类日韩欧美亚洲嫩草| 国内精品久久久久精免费| 国产男靠女视频免费网站| 又大又爽又粗| 美国免费a级毛片| 在线免费观看的www视频| 国产区一区二久久| 亚洲国产高清在线一区二区三 | 精品一区二区三区四区五区乱码| 在线观看www视频免费| 亚洲九九香蕉| 丁香六月欧美| 岛国视频午夜一区免费看| 女人被狂操c到高潮| 亚洲精品在线美女| 老鸭窝网址在线观看| 一区二区三区国产精品乱码| 99国产极品粉嫩在线观看| 午夜a级毛片| 色播在线永久视频| 亚洲最大成人中文| 天堂动漫精品| 正在播放国产对白刺激| 欧美日本亚洲视频在线播放| av福利片在线| 精品一区二区三区视频在线观看免费| 大型av网站在线播放| 国产精品电影一区二区三区| 精品国产美女av久久久久小说| 国内精品久久久久久久电影| 久久久久久人人人人人| 99久久无色码亚洲精品果冻| 韩国av一区二区三区四区| 国产91精品成人一区二区三区| 成人国语在线视频| 亚洲精品中文字幕在线视频| 午夜福利在线观看吧| 成人精品一区二区免费| 宅男免费午夜| 欧美中文综合在线视频| 午夜福利在线在线| 国产亚洲av嫩草精品影院| 国产精品乱码一区二三区的特点| a级毛片a级免费在线| 级片在线观看| 国产高清激情床上av| 国产日本99.免费观看| 精品国产亚洲在线| 亚洲国产精品成人综合色| 国产激情偷乱视频一区二区| 亚洲国产欧洲综合997久久, | 巨乳人妻的诱惑在线观看| 一级作爱视频免费观看| 亚洲欧美精品综合一区二区三区| 两个人视频免费观看高清| 久久国产亚洲av麻豆专区| 国产精品久久久久久精品电影 | 一进一出好大好爽视频| 欧美黄色淫秽网站| 亚洲专区中文字幕在线| 国产精华一区二区三区| 手机成人av网站| 亚洲人成77777在线视频| 最近最新中文字幕大全电影3 | 99久久国产精品久久久| 夜夜看夜夜爽夜夜摸| 18禁美女被吸乳视频| 中文字幕av电影在线播放| 久久天堂一区二区三区四区| 国产激情久久老熟女| 亚洲精品在线观看二区| 亚洲色图 男人天堂 中文字幕| 成年版毛片免费区| 少妇裸体淫交视频免费看高清 | 久久久久国产精品人妻aⅴ院| 欧美性猛交黑人性爽| 日本一本二区三区精品| 久久精品夜夜夜夜夜久久蜜豆 | 正在播放国产对白刺激| 美女午夜性视频免费| 美女 人体艺术 gogo| 黄网站色视频无遮挡免费观看| 99国产精品一区二区蜜桃av| 99国产精品一区二区三区| 一个人免费在线观看的高清视频| 麻豆一二三区av精品| 亚洲av美国av| 两个人看的免费小视频| 两个人免费观看高清视频| 狂野欧美激情性xxxx| 亚洲av成人不卡在线观看播放网| 久99久视频精品免费| 欧美成人午夜精品| 免费搜索国产男女视频| 国产精品亚洲美女久久久| 国产三级黄色录像| 亚洲在线自拍视频| 香蕉av资源在线| 婷婷丁香在线五月| 久久久国产欧美日韩av| 国内精品久久久久久久电影| 欧美又色又爽又黄视频| 久99久视频精品免费| 久久久久久免费高清国产稀缺| 精品久久久久久久末码| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利免费观看在线| 国产激情欧美一区二区| 桃色一区二区三区在线观看| 黄色a级毛片大全视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产av一区在线观看免费| 日韩欧美 国产精品| 亚洲无线在线观看| 亚洲天堂国产精品一区在线| 国产免费av片在线观看野外av| 欧美成人免费av一区二区三区| 波多野结衣高清无吗| 精品熟女少妇八av免费久了| 一区二区三区高清视频在线| 日韩有码中文字幕| 丝袜在线中文字幕| 国产色视频综合| av免费在线观看网站| 性色av乱码一区二区三区2| 9191精品国产免费久久| 久久久久久国产a免费观看| 我的亚洲天堂| 国产精品一区二区三区四区久久 | 精品国产超薄肉色丝袜足j| 久久精品成人免费网站| 午夜激情av网站| av视频在线观看入口| 国产伦一二天堂av在线观看| 99国产综合亚洲精品| 免费一级毛片在线播放高清视频| 亚洲av日韩精品久久久久久密| 免费观看精品视频网站| 波多野结衣高清作品| 国产在线精品亚洲第一网站| 亚洲成av片中文字幕在线观看| 国产高清videossex| 久久亚洲精品不卡| 国产精品久久电影中文字幕| 草草在线视频免费看| 特大巨黑吊av在线直播 | xxx96com| 亚洲av成人av| 精品国内亚洲2022精品成人| 久久久久久国产a免费观看| 可以在线观看毛片的网站| 国产av不卡久久| 中文字幕久久专区| 日韩三级视频一区二区三区| 久久 成人 亚洲| 日韩三级视频一区二区三区| 成人精品一区二区免费| 亚洲无线在线观看| 亚洲午夜理论影院| 久久香蕉国产精品| 一区二区三区国产精品乱码| 欧美性猛交黑人性爽| 一区二区三区国产精品乱码| 极品教师在线免费播放| 成人永久免费在线观看视频| а√天堂www在线а√下载| 老司机深夜福利视频在线观看| 欧美激情高清一区二区三区| 女性被躁到高潮视频| 日本a在线网址| 天天一区二区日本电影三级| 亚洲精品中文字幕一二三四区| 亚洲自偷自拍图片 自拍| 精品国产乱码久久久久久男人| 国产精品亚洲美女久久久| 欧美激情高清一区二区三区| 久久久久国产精品人妻aⅴ院| 国产伦人伦偷精品视频| 在线十欧美十亚洲十日本专区| 免费女性裸体啪啪无遮挡网站| 搡老熟女国产l中国老女人| 免费高清在线观看日韩| 女人爽到高潮嗷嗷叫在线视频| 国产精品av久久久久免费| 成人精品一区二区免费| 黑丝袜美女国产一区| 在线观看舔阴道视频| 亚洲精品美女久久av网站| 精品日产1卡2卡| 亚洲av中文字字幕乱码综合 | 久久欧美精品欧美久久欧美| 欧美zozozo另类| 精品第一国产精品| 免费在线观看完整版高清| 久久久久久久久免费视频了| 19禁男女啪啪无遮挡网站| 可以在线观看毛片的网站| 色综合站精品国产| 国产成人精品久久二区二区免费| 久久人妻福利社区极品人妻图片| 欧美激情久久久久久爽电影| av视频在线观看入口| 女性被躁到高潮视频| 啪啪无遮挡十八禁网站| 国产黄片美女视频| 亚洲精华国产精华精| 女同久久另类99精品国产91| 久9热在线精品视频| 不卡一级毛片| 一进一出抽搐动态| 午夜影院日韩av| 亚洲欧美精品综合久久99| 日韩 欧美 亚洲 中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 欧美黑人精品巨大| 最近最新中文字幕大全电影3 | 免费看日本二区| av在线播放免费不卡| 黄色丝袜av网址大全| 国产成人系列免费观看| www国产在线视频色| 一a级毛片在线观看| 久久久久精品国产欧美久久久| 亚洲男人天堂网一区| 怎么达到女性高潮| 嫩草影院精品99| 熟妇人妻久久中文字幕3abv| 国产精品久久视频播放| 欧美三级亚洲精品| 可以在线观看的亚洲视频| 欧美又色又爽又黄视频| 欧美中文日本在线观看视频| 亚洲国产毛片av蜜桃av| 两人在一起打扑克的视频| 成人av一区二区三区在线看| 日本 欧美在线| 搞女人的毛片| 最近最新中文字幕大全免费视频| 国产视频一区二区在线看| 国产欧美日韩一区二区精品| 免费看a级黄色片| 亚洲欧美精品综合久久99| 国产精品久久久人人做人人爽| 精品久久久久久久久久免费视频| 国产一区二区在线av高清观看| 亚洲中文日韩欧美视频| 亚洲九九香蕉| 婷婷精品国产亚洲av| 俺也久久电影网| 午夜福利高清视频| 成年女人毛片免费观看观看9| 一区福利在线观看| 丝袜人妻中文字幕| 亚洲成人久久性| 日韩一卡2卡3卡4卡2021年| 国产精品久久久av美女十八| 中文在线观看免费www的网站 | 欧美乱妇无乱码| 97碰自拍视频| 亚洲男人的天堂狠狠| 18禁国产床啪视频网站| 成人精品一区二区免费| 欧美zozozo另类| 亚洲 欧美 日韩 在线 免费| 人人妻人人澡人人看| 天天躁狠狠躁夜夜躁狠狠躁| 国产片内射在线| 啦啦啦 在线观看视频| av电影中文网址| 色播亚洲综合网| 欧美日韩黄片免| x7x7x7水蜜桃| 午夜a级毛片| 黄色a级毛片大全视频| 亚洲人成77777在线视频| 女性被躁到高潮视频| 黄片大片在线免费观看| 人人妻,人人澡人人爽秒播| 特大巨黑吊av在线直播 | 韩国av一区二区三区四区| 色播亚洲综合网| 老汉色∧v一级毛片| 日韩一卡2卡3卡4卡2021年| av欧美777| 啦啦啦观看免费观看视频高清| 黄色视频,在线免费观看| 国产黄片美女视频| 可以在线观看的亚洲视频| 搡老岳熟女国产| 欧美国产日韩亚洲一区| 国产亚洲欧美在线一区二区| 成人国产综合亚洲| tocl精华| 国产在线观看jvid| 淫妇啪啪啪对白视频| 亚洲成人久久爱视频| 色在线成人网| 两个人看的免费小视频| 久久中文字幕人妻熟女| 国产精品亚洲一级av第二区| 人人妻人人澡人人看| av在线播放免费不卡| 精品国产乱子伦一区二区三区| 欧美性猛交╳xxx乱大交人| 国产精品免费视频内射| 国产精品美女特级片免费视频播放器 | 搡老熟女国产l中国老女人| 中亚洲国语对白在线视频| 999久久久国产精品视频| 午夜福利成人在线免费观看| 老汉色av国产亚洲站长工具| 欧美精品啪啪一区二区三区| tocl精华| 18禁黄网站禁片免费观看直播| 久久午夜综合久久蜜桃| 国产精品精品国产色婷婷| 少妇 在线观看| bbb黄色大片| 午夜福利欧美成人| 啦啦啦韩国在线观看视频| 狠狠狠狠99中文字幕| 精品国产亚洲在线| 久久精品国产99精品国产亚洲性色| 欧美午夜高清在线| 亚洲国产日韩欧美精品在线观看 | 国内少妇人妻偷人精品xxx网站 | 亚洲免费av在线视频| 国产91精品成人一区二区三区| 亚洲国产欧美网| 亚洲国产欧美一区二区综合| 自线自在国产av| 日本在线视频免费播放| 露出奶头的视频| 亚洲专区国产一区二区| 免费在线观看亚洲国产| 日韩欧美一区二区三区在线观看| 亚洲av中文字字幕乱码综合 | 一进一出抽搐gif免费好疼| e午夜精品久久久久久久| 色综合亚洲欧美另类图片| 国产亚洲精品综合一区在线观看 | 九色国产91popny在线| 丝袜美腿诱惑在线| 夜夜看夜夜爽夜夜摸| 成人国语在线视频| 亚洲精品美女久久久久99蜜臀| 国产伦一二天堂av在线观看| 欧美乱妇无乱码| 18禁黄网站禁片免费观看直播| 精品国产乱子伦一区二区三区| 狠狠狠狠99中文字幕| 99久久久亚洲精品蜜臀av| 成人亚洲精品av一区二区| 亚洲av成人av| 手机成人av网站| 午夜老司机福利片| 国产一卡二卡三卡精品| 一级毛片高清免费大全| av超薄肉色丝袜交足视频| 神马国产精品三级电影在线观看 | 少妇 在线观看| 99久久精品国产亚洲精品| 亚洲av美国av| 日韩精品中文字幕看吧| 欧美日韩一级在线毛片| 一a级毛片在线观看| 成人精品一区二区免费| 日韩大码丰满熟妇| 亚洲国产精品999在线| 国内精品久久久久精免费| 久久久久久久精品吃奶| 日韩精品青青久久久久久| 正在播放国产对白刺激| 免费观看人在逋| 免费在线观看影片大全网站| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av日韩精品久久久久久密| av免费在线观看网站| e午夜精品久久久久久久| 亚洲成av片中文字幕在线观看| 又黄又爽又免费观看的视频| 国产成人一区二区三区免费视频网站| 黄色女人牲交| 看黄色毛片网站| 在线观看免费日韩欧美大片| 亚洲国产精品sss在线观看| 久久青草综合色| 免费人成视频x8x8入口观看| 91字幕亚洲| 欧美中文综合在线视频| 听说在线观看完整版免费高清| 亚洲中文字幕一区二区三区有码在线看 | 国产三级黄色录像| 女同久久另类99精品国产91| 亚洲第一欧美日韩一区二区三区| 19禁男女啪啪无遮挡网站| 少妇 在线观看| 午夜激情福利司机影院| 久久精品国产清高在天天线| 亚洲va日本ⅴa欧美va伊人久久| 国产精品美女特级片免费视频播放器 | 亚洲专区国产一区二区| 欧美丝袜亚洲另类 | 国产精品久久久人人做人人爽| 一本综合久久免费| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇被粗大的猛进出69影院| 搞女人的毛片| 无人区码免费观看不卡| 美女高潮喷水抽搐中文字幕| 免费无遮挡裸体视频| 亚洲成a人片在线一区二区| 天堂√8在线中文| 亚洲欧美精品综合一区二区三区| 亚洲一区高清亚洲精品| 国产精品九九99| 曰老女人黄片| 国产视频内射|