• <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é)用多功能支架
    窝窝影院91人妻| 日本免费一区二区三区高清不卡| 91字幕亚洲| 每晚都被弄得嗷嗷叫到高潮| 亚洲乱码一区二区免费版| 白带黄色成豆腐渣| 国产免费一级a男人的天堂| www.色视频.com| 美女黄网站色视频| 在线免费观看的www视频| 亚洲真实伦在线观看| 熟妇人妻久久中文字幕3abv| 欧美大码av| 91麻豆av在线| 免费在线观看成人毛片| 男人和女人高潮做爰伦理| 亚洲五月婷婷丁香| 91在线精品国自产拍蜜月 | 女警被强在线播放| 欧美zozozo另类| 极品教师在线免费播放| 亚洲内射少妇av| 日韩成人在线观看一区二区三区| 亚洲精品粉嫩美女一区| 欧美成狂野欧美在线观看| 天堂影院成人在线观看| 一个人看视频在线观看www免费 | tocl精华| svipshipincom国产片| 亚洲欧美日韩东京热| 身体一侧抽搐| 久久精品国产清高在天天线| 一边摸一边抽搐一进一小说| 久久久国产成人精品二区| 久久精品国产综合久久久| 一个人看的www免费观看视频| 日韩欧美 国产精品| 国内精品美女久久久久久| 色吧在线观看| 女人十人毛片免费观看3o分钟| 午夜精品久久久久久毛片777| 精品午夜福利视频在线观看一区| 9191精品国产免费久久| 日韩中文字幕欧美一区二区| 天堂√8在线中文| 69av精品久久久久久| 99精品久久久久人妻精品| bbb黄色大片| 在线观看午夜福利视频| 99riav亚洲国产免费| 久久久久九九精品影院| 久久午夜亚洲精品久久| 午夜精品在线福利| 18美女黄网站色大片免费观看| 国产精品女同一区二区软件 | 蜜桃亚洲精品一区二区三区| aaaaa片日本免费| 三级毛片av免费| 99在线视频只有这里精品首页| 国产亚洲精品久久久com| a在线观看视频网站| 最近最新免费中文字幕在线| 欧美zozozo另类| 国产精品日韩av在线免费观看| 亚洲天堂国产精品一区在线| 亚洲人成伊人成综合网2020| 久久精品91无色码中文字幕| 色综合亚洲欧美另类图片| 十八禁网站免费在线| 久久久久亚洲av毛片大全| 国产成人a区在线观看| 成人鲁丝片一二三区免费| 母亲3免费完整高清在线观看| 日韩欧美三级三区| 人人妻,人人澡人人爽秒播| 欧美绝顶高潮抽搐喷水| 精品久久久久久久久久久久久| 首页视频小说图片口味搜索| 久久久久久久精品吃奶| 亚洲av一区综合| 深夜精品福利| 日韩精品青青久久久久久| 搡老熟女国产l中国老女人| 亚洲精品456在线播放app | 51国产日韩欧美| 岛国在线免费视频观看| 99精品欧美一区二区三区四区| 国产精品av视频在线免费观看| 在线观看免费视频日本深夜| 亚洲午夜理论影院| 久久午夜亚洲精品久久| 免费搜索国产男女视频| 一区福利在线观看| 手机成人av网站| 欧美最黄视频在线播放免费| 免费人成视频x8x8入口观看| 99热只有精品国产| av天堂中文字幕网| 久久久精品大字幕| 99精品久久久久人妻精品| 最近最新中文字幕大全免费视频| 一个人免费在线观看电影| 91字幕亚洲| 夜夜躁狠狠躁天天躁| 免费人成视频x8x8入口观看| av在线蜜桃| 午夜日韩欧美国产| 九九在线视频观看精品| 亚洲 国产 在线| 亚洲av成人不卡在线观看播放网| 身体一侧抽搐| 人人妻,人人澡人人爽秒播| 久久国产乱子伦精品免费另类| 亚洲精品一区av在线观看| 3wmmmm亚洲av在线观看| 亚洲精品成人久久久久久| 最好的美女福利视频网| 国产男靠女视频免费网站| 亚洲精品成人久久久久久| 国产真实伦视频高清在线观看 | 99久久无色码亚洲精品果冻| 精品日产1卡2卡| 一级毛片高清免费大全| 国产高潮美女av| 国产精品三级大全| 黄色片一级片一级黄色片| 国产aⅴ精品一区二区三区波| 国产毛片a区久久久久| 最好的美女福利视频网| 噜噜噜噜噜久久久久久91| 好男人在线观看高清免费视频| 亚洲在线观看片| 久久国产精品影院| 三级国产精品欧美在线观看| 色综合亚洲欧美另类图片| 男女视频在线观看网站免费| 国产v大片淫在线免费观看| 一个人看视频在线观看www免费 | 老司机福利观看| 国产亚洲精品一区二区www| 大型黄色视频在线免费观看| 亚洲精品国产精品久久久不卡| 国产成人影院久久av| 人妻久久中文字幕网| 欧美绝顶高潮抽搐喷水| 最新中文字幕久久久久| 久久天躁狠狠躁夜夜2o2o| 色av中文字幕| 免费av不卡在线播放| 长腿黑丝高跟| 日韩高清综合在线| 亚洲天堂国产精品一区在线| 成人av一区二区三区在线看| 在线免费观看的www视频| 精品熟女少妇八av免费久了| 欧美激情久久久久久爽电影| 久久精品影院6| e午夜精品久久久久久久| 少妇高潮的动态图| 19禁男女啪啪无遮挡网站| 日本撒尿小便嘘嘘汇集6| 成人国产综合亚洲| 国产一区二区亚洲精品在线观看| 色尼玛亚洲综合影院| 国产99白浆流出| 黄色视频,在线免费观看| 无人区码免费观看不卡| 精品福利观看| 亚洲激情在线av| 国产成人影院久久av| 国产精品一区二区免费欧美| 黄色成人免费大全| 色综合站精品国产| 少妇人妻精品综合一区二区 | 女人十人毛片免费观看3o分钟| 成年女人毛片免费观看观看9| 国产探花极品一区二区| 国产91精品成人一区二区三区| 久久天躁狠狠躁夜夜2o2o| 国产高清videossex| 激情在线观看视频在线高清| 久久国产乱子伦精品免费另类| 性色av乱码一区二区三区2| 欧美一级毛片孕妇| 欧美最新免费一区二区三区 | 亚洲中文字幕一区二区三区有码在线看| 18禁国产床啪视频网站| 成年人黄色毛片网站| 丁香六月欧美| av福利片在线观看| 国产黄a三级三级三级人| 日韩高清综合在线| 国产探花在线观看一区二区| 18禁黄网站禁片免费观看直播| 久久久久久久精品吃奶| 午夜精品一区二区三区免费看| 亚洲成人久久性| 色综合站精品国产| 97超级碰碰碰精品色视频在线观看| 国产午夜精品论理片| 亚洲精品日韩av片在线观看 | 亚洲成a人片在线一区二区| 欧美在线黄色| 色综合亚洲欧美另类图片| 久久久久国内视频| 精品久久久久久久久久免费视频| 91久久精品电影网| 欧美成人性av电影在线观看| 精品久久久久久久久久久久久| 亚洲av二区三区四区| 99精品欧美一区二区三区四区| 99在线视频只有这里精品首页| 一级毛片高清免费大全| 精华霜和精华液先用哪个| 亚洲欧美日韩卡通动漫| 最新美女视频免费是黄的| 又紧又爽又黄一区二区| 美女高潮喷水抽搐中文字幕| 欧美另类亚洲清纯唯美| 麻豆成人av在线观看| 国产91精品成人一区二区三区| tocl精华| 久久精品国产亚洲av涩爱 | 欧美日本亚洲视频在线播放| 国产精品 国内视频| 精品人妻一区二区三区麻豆 | 欧美bdsm另类| 国产成人影院久久av| 亚洲在线观看片| 免费在线观看影片大全网站| 99国产综合亚洲精品| 欧美一级毛片孕妇| 国产一区二区三区视频了| 国产一区二区在线av高清观看| 亚洲最大成人手机在线| 一区二区三区国产精品乱码| 国产成人系列免费观看| 亚洲av日韩精品久久久久久密| 国产黄色小视频在线观看| 免费在线观看日本一区| 在线天堂最新版资源| 国产91精品成人一区二区三区| 一区二区三区免费毛片| 欧美成人免费av一区二区三区| 国产精品99久久久久久久久| ponron亚洲| 99热只有精品国产| 性色av乱码一区二区三区2| 国产麻豆成人av免费视频| 免费在线观看亚洲国产| 亚洲精品乱码久久久v下载方式 | 我的老师免费观看完整版| 国产真实伦视频高清在线观看 | 天堂网av新在线| 真人做人爱边吃奶动态| 小蜜桃在线观看免费完整版高清| 国产精品亚洲一级av第二区| 国产伦精品一区二区三区视频9 | 99热精品在线国产| 国产不卡一卡二| 中国美女看黄片| 又黄又粗又硬又大视频| 久久中文看片网| 可以在线观看的亚洲视频| 国产伦在线观看视频一区| 热99在线观看视频| 别揉我奶头~嗯~啊~动态视频| 真人做人爱边吃奶动态| 精华霜和精华液先用哪个| 男插女下体视频免费在线播放| 可以在线观看毛片的网站| 国产真实伦视频高清在线观看 | 国产高潮美女av| 国产精品 国内视频| 国产精品亚洲美女久久久| 亚洲精品粉嫩美女一区| 国产综合懂色| 在线观看一区二区三区| 免费看日本二区| 亚洲av成人不卡在线观看播放网| 久久天躁狠狠躁夜夜2o2o| 中文字幕久久专区| 国产私拍福利视频在线观看| 一级黄片播放器| 少妇高潮的动态图| 久久精品国产清高在天天线| 99热6这里只有精品| 国产乱人伦免费视频| 波多野结衣高清作品| 黄色片一级片一级黄色片| 99国产精品一区二区三区| 19禁男女啪啪无遮挡网站| 91麻豆av在线| 又粗又爽又猛毛片免费看| 丝袜美腿在线中文| 99热精品在线国产| 国产高清视频在线观看网站| 51国产日韩欧美| 国产探花极品一区二区| 亚洲精品粉嫩美女一区| 法律面前人人平等表现在哪些方面| 亚洲成人精品中文字幕电影| 久久久色成人| 国产伦在线观看视频一区| 高清日韩中文字幕在线| 熟女人妻精品中文字幕| 亚洲精品在线美女| 久久6这里有精品| 精品国内亚洲2022精品成人| eeuss影院久久| 久久久久性生活片| 国语自产精品视频在线第100页| 国产av不卡久久| 欧美极品一区二区三区四区| 欧美又色又爽又黄视频| 高清毛片免费观看视频网站| 一区二区三区国产精品乱码| 丰满人妻一区二区三区视频av | 男人和女人高潮做爰伦理| 高清在线国产一区| 日韩成人在线观看一区二区三区| 成人av在线播放网站| 一个人看视频在线观看www免费 | 黄色丝袜av网址大全| 午夜免费成人在线视频| 好男人在线观看高清免费视频| 99久久精品一区二区三区| 日韩欧美一区二区三区在线观看| 99久久精品国产亚洲精品| 国产老妇女一区| 偷拍熟女少妇极品色| 制服人妻中文乱码| 亚洲精品日韩av片在线观看 | 日韩欧美 国产精品| 精品免费久久久久久久清纯| 日日夜夜操网爽| 午夜福利欧美成人| 一进一出抽搐动态| 在线天堂最新版资源| 我要搜黄色片| 人妻久久中文字幕网| 国产午夜福利久久久久久| 欧美日韩乱码在线| 99精品在免费线老司机午夜| 91麻豆精品激情在线观看国产| aaaaa片日本免费| 久久精品亚洲精品国产色婷小说| 一二三四社区在线视频社区8| 内地一区二区视频在线| 色综合亚洲欧美另类图片| 国产精品 欧美亚洲| av欧美777| 成年女人毛片免费观看观看9| 午夜福利18| 欧美丝袜亚洲另类 | 欧美日韩综合久久久久久 | 精品乱码久久久久久99久播| 欧美大码av| 国产私拍福利视频在线观看| 精品久久久久久久久久免费视频| 亚洲熟妇熟女久久| 日日摸夜夜添夜夜添小说| 19禁男女啪啪无遮挡网站| 国产激情欧美一区二区| 琪琪午夜伦伦电影理论片6080| 欧美大码av| 欧美在线黄色| 欧美3d第一页| 久久精品人妻少妇| 欧美性感艳星| netflix在线观看网站| 最近视频中文字幕2019在线8| 国产在线精品亚洲第一网站| 国产一区二区在线观看日韩 | 内地一区二区视频在线| 男女做爰动态图高潮gif福利片| aaaaa片日本免费| 免费无遮挡裸体视频| 色综合站精品国产| 色综合婷婷激情| 男女床上黄色一级片免费看| 国产在视频线在精品| 亚洲乱码一区二区免费版| 亚洲五月天丁香| av天堂中文字幕网| 亚洲av第一区精品v没综合| av专区在线播放| 丰满人妻一区二区三区视频av | 午夜日韩欧美国产| 18禁在线播放成人免费| 窝窝影院91人妻| 国产91精品成人一区二区三区| 欧美成人免费av一区二区三区| 国产乱人视频| 热99re8久久精品国产| 全区人妻精品视频| 亚洲精品影视一区二区三区av| 特级一级黄色大片| 搡老妇女老女人老熟妇| 国产成+人综合+亚洲专区| 国产精品日韩av在线免费观看| 特大巨黑吊av在线直播| 精品国产三级普通话版| 一卡2卡三卡四卡精品乱码亚洲| 97碰自拍视频| 可以在线观看的亚洲视频| 久久精品人妻少妇| 欧美日韩一级在线毛片| 国内精品美女久久久久久| 18美女黄网站色大片免费观看| 伊人久久精品亚洲午夜| 欧美bdsm另类| 男人的好看免费观看在线视频| 午夜精品一区二区三区免费看| 男女做爰动态图高潮gif福利片| 狂野欧美激情性xxxx| 悠悠久久av| 成人性生交大片免费视频hd| 亚洲精品在线美女| 日本一二三区视频观看| 一进一出抽搐gif免费好疼| 美女cb高潮喷水在线观看| 亚洲精品一卡2卡三卡4卡5卡| 制服丝袜大香蕉在线| 免费av观看视频| 国产一区二区在线av高清观看| 国产成人a区在线观看| 国产欧美日韩精品亚洲av| 日本一本二区三区精品| 免费av毛片视频| 精品熟女少妇八av免费久了| 2021天堂中文幕一二区在线观| 99精品久久久久人妻精品| 精品一区二区三区视频在线观看免费| 淫秽高清视频在线观看| 国产精品 欧美亚洲| 中文字幕人成人乱码亚洲影| www国产在线视频色| 国产三级在线视频| 久久精品国产99精品国产亚洲性色| 免费在线观看影片大全网站| 日韩亚洲欧美综合| 夜夜看夜夜爽夜夜摸| 欧美色视频一区免费| 精品一区二区三区人妻视频| 嫁个100分男人电影在线观看| 内射极品少妇av片p| 少妇的逼好多水| 日本一本二区三区精品| 97超视频在线观看视频| 国产精品永久免费网站| www日本黄色视频网| 精品无人区乱码1区二区| 中文字幕av成人在线电影| 国产淫片久久久久久久久 | 午夜影院日韩av| 99久久综合精品五月天人人| 国产精品,欧美在线| 一级毛片女人18水好多| 听说在线观看完整版免费高清| 日韩精品中文字幕看吧| 欧美成人性av电影在线观看| 在线观看66精品国产| 亚洲av一区综合| 波野结衣二区三区在线 | 免费看十八禁软件| 成人欧美大片| 美女大奶头视频| 禁无遮挡网站| 一区二区三区国产精品乱码| 男人舔女人下体高潮全视频| 一级作爱视频免费观看| 亚洲国产精品合色在线| 床上黄色一级片| 欧美国产日韩亚洲一区| 欧美一区二区精品小视频在线| 欧美成人性av电影在线观看| 最新在线观看一区二区三区| 婷婷精品国产亚洲av| 午夜福利免费观看在线| 高清在线国产一区| 久久欧美精品欧美久久欧美| 一个人观看的视频www高清免费观看| 免费观看精品视频网站| 久久精品国产自在天天线| 最近在线观看免费完整版| 国产成+人综合+亚洲专区| 老汉色av国产亚洲站长工具| 亚洲人成伊人成综合网2020| 激情在线观看视频在线高清| 91麻豆av在线| 动漫黄色视频在线观看| 欧美在线一区亚洲| 一级毛片女人18水好多| www日本黄色视频网| 一进一出抽搐动态| 三级男女做爰猛烈吃奶摸视频| 国产一区二区亚洲精品在线观看| 日本在线视频免费播放| 99国产极品粉嫩在线观看| 免费观看的影片在线观看| 韩国av一区二区三区四区| 久久精品国产清高在天天线| 老司机午夜福利在线观看视频| 亚洲成人久久性| 日本 欧美在线| 国产亚洲精品综合一区在线观看| 高清日韩中文字幕在线| 久久草成人影院| 人妻夜夜爽99麻豆av| 中亚洲国语对白在线视频| 午夜日韩欧美国产| 欧美性感艳星| 麻豆一二三区av精品| 欧美午夜高清在线| 久久这里只有精品中国| 亚洲狠狠婷婷综合久久图片| av黄色大香蕉| 亚洲国产精品sss在线观看| av片东京热男人的天堂| 国产真实伦视频高清在线观看 | 黄色女人牲交| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 老熟妇仑乱视频hdxx| 亚洲欧美日韩高清专用| 亚洲专区中文字幕在线| 免费在线观看亚洲国产| 国产真实伦视频高清在线观看 | 99久久精品热视频| 亚洲男人的天堂狠狠| 国产精品一区二区三区四区久久| 少妇的逼水好多| 欧美性感艳星| 伊人久久精品亚洲午夜| 午夜精品在线福利| 国产精品影院久久| 首页视频小说图片口味搜索| www日本黄色视频网| 国产探花在线观看一区二区| 九九热线精品视视频播放| 一本综合久久免费| 两个人视频免费观看高清| 12—13女人毛片做爰片一| 中文字幕精品亚洲无线码一区| 久久久久亚洲av毛片大全| 日本撒尿小便嘘嘘汇集6| 国产三级中文精品| 蜜桃亚洲精品一区二区三区| 久久这里只有精品中国| 99热这里只有是精品50| 国产精品久久久人人做人人爽| 精品久久久久久久人妻蜜臀av| 在线观看美女被高潮喷水网站 | 日韩大尺度精品在线看网址| 国产亚洲欧美在线一区二区| 久久精品国产亚洲av涩爱 | 久久国产精品影院| 精品无人区乱码1区二区| 一a级毛片在线观看| av黄色大香蕉| 日韩欧美三级三区| 99精品在免费线老司机午夜| 黄色视频,在线免费观看| 亚洲国产精品合色在线| 亚洲精品久久国产高清桃花| 国产一级毛片七仙女欲春2| 18禁黄网站禁片午夜丰满| 亚洲成人精品中文字幕电影| 国产成人啪精品午夜网站| 国产精品av视频在线免费观看| 亚洲国产欧美网| АⅤ资源中文在线天堂| 久99久视频精品免费| 在线看三级毛片| 国产99白浆流出| 免费观看的影片在线观看| 国产高清视频在线播放一区| 少妇人妻精品综合一区二区 | 99久久综合精品五月天人人| eeuss影院久久| 欧美日韩亚洲国产一区二区在线观看| 欧美午夜高清在线| 国产激情欧美一区二区| 欧美一区二区亚洲| 少妇高潮的动态图| 久久性视频一级片| 亚洲 国产 在线| 国产v大片淫在线免费观看| 日韩欧美一区二区三区在线观看| 免费观看精品视频网站| 国产国拍精品亚洲av在线观看 | 国产精品亚洲av一区麻豆| 在线观看一区二区三区| 国产在视频线在精品| 国产成年人精品一区二区| 中文字幕av在线有码专区| 一个人免费在线观看电影| 波野结衣二区三区在线 | 午夜视频国产福利| 午夜亚洲福利在线播放| 小说图片视频综合网站| av女优亚洲男人天堂| 人人妻,人人澡人人爽秒播| 亚洲av电影不卡..在线观看| 免费电影在线观看免费观看| 看片在线看免费视频| 日本a在线网址| 国产视频一区二区在线看| 国产美女午夜福利| 3wmmmm亚洲av在线观看| 99热这里只有是精品50| 亚洲性夜色夜夜综合| 婷婷精品国产亚洲av| 精品一区二区三区视频在线 | 国产成人欧美在线观看|