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

    三絲型質(zhì)點(diǎn)振速傳感器的模型分析

    2022-09-23 01:49:36王綏昌郭軒瑞朱林輝
    測(cè)控技術(shù) 2022年9期
    關(guān)鍵詞:質(zhì)點(diǎn)聲波穩(wěn)態(tài)

    王綏昌,郭軒瑞,朱林輝,張 穎,郭 航*

    (1.廈門大學(xué) 薩本棟微米納米科學(xué)技術(shù)研究院,福建 廈門 361005; 2.廈門大學(xué)馬來(lái)西亞分校 新能源科學(xué)與工程系,馬來(lái)西亞 雪邦 43900;3.廈門大學(xué)馬來(lái)西亞分校 電子信息工程系,馬來(lái)西亞 雪邦 43900)

    聲學(xué)測(cè)量系統(tǒng)在環(huán)境監(jiān)測(cè)、汽車工業(yè)、軍事國(guó)防等領(lǐng)域都發(fā)揮著重要作用。要完整地描述一個(gè)聲場(chǎng),需要知道標(biāo)量聲壓和矢量質(zhì)點(diǎn)振速兩個(gè)參數(shù)。1994年荷蘭de Bree提出了一種采用微機(jī)電系統(tǒng)即MEMS技術(shù)來(lái)加工制造的質(zhì)點(diǎn)振速傳感器(Particle Velocity Sensor,PVS)[1],它以氣團(tuán)運(yùn)動(dòng)代替質(zhì)量塊振動(dòng),通過(guò)測(cè)量熱源對(duì)稱位置的溫度差變化,從而確定引起器件微通道內(nèi)氣體流動(dòng)的聲場(chǎng)的質(zhì)點(diǎn)振速。在隨后的十多年間,以此技術(shù)為基礎(chǔ)的荷蘭Microflown Technology 公司所開發(fā)的由聲壓計(jì)與PVS所構(gòu)成的聲學(xué)探頭,即矢量傳聲器,在軍用與民用領(lǐng)域都得到了廣泛應(yīng)用[2-3],目前仍在深入研究當(dāng)中。國(guó)內(nèi)對(duì)PVS的研究起步較晚,2014年中國(guó)電子科技集團(tuán)公司第三研究所對(duì)矢量傳聲器技術(shù)進(jìn)行了研究探索,微加工制備了兩絲型PVS[4],對(duì)其性能進(jìn)行了測(cè)試分析與優(yōu)化;2019年北京大學(xué)微加工制備了三絲型PVS,進(jìn)行了測(cè)試研究[5-6]。國(guó)內(nèi)其他單位如中國(guó)計(jì)量學(xué)院[7]、廈門大學(xué)[8]等也做了深入研究。

    PVS是矢量傳聲器的核心部件,目前主要有兩種結(jié)構(gòu)形式:一種是兩絲型,包含兩根平行排列、具有一定間距的熱阻絲結(jié)構(gòu),如圖1所示,這兩根絲既是加熱絲,又被用作傳感器的敏感絲;另一種結(jié)構(gòu)是三絲型,如圖2所示,是包含三根平行排列具有一定間距的熱阻絲結(jié)構(gòu),兩邊的是敏感絲,中間的是加熱絲。雖然兩絲型結(jié)構(gòu)靈敏度比三絲型較高,但不是特別顯著,并且三絲型與兩絲型相比將熱源和傳感器隔離開來(lái),可有效避免熱噪聲的影響。

    圖1 兩絲型PVS橫截面及俯視示意圖

    圖2 三絲型PVS橫截面及俯視示意圖

    PVS是目前唯一的一種能夠直接測(cè)量空氣質(zhì)點(diǎn)振速的傳感器,其靈敏度、噪聲、功耗、工作頻段等直接決定矢量傳聲器的性能。因此,其分析設(shè)計(jì)至關(guān)重要。

    目前已有對(duì)兩絲型PVS的分析模型的研究,對(duì)于三絲型則較少。本文對(duì)三絲型PVS的模型進(jìn)行研究,從穩(wěn)態(tài)與瞬態(tài)的熱傳導(dǎo)分析出發(fā),建立起三絲型PVS的解析模型,并與數(shù)值模型進(jìn)行對(duì)比驗(yàn)證,從而為三絲型PVS的優(yōu)化設(shè)計(jì)與微加工制造提供理論基礎(chǔ)。

    1 三絲型PVS的解析分析

    如圖2所示,在三絲型PVS中,兩根熱阻絲在中間的加熱絲的左右兩邊對(duì)稱平行放置,作為敏感單元。整個(gè)傳感器的腔壁為熱絕緣且不可滲透的邊界,使通道內(nèi)部盡可能小地受外界環(huán)境影響。當(dāng)沒(méi)有聲場(chǎng)且加熱絲所受恒定電壓不變時(shí),微通道內(nèi)部溫度場(chǎng)處于穩(wěn)態(tài),兩根敏感絲的溫度差為0,即ΔT=0。當(dāng)外界施加一個(gè)聲場(chǎng)時(shí),離聲源較近一側(cè)的粒子會(huì)產(chǎn)生振動(dòng)[9],從而在靠近聲源一側(cè)的微通道入口處引起通道內(nèi)介質(zhì)發(fā)生受迫對(duì)流運(yùn)動(dòng),形成一個(gè)熱擾動(dòng)場(chǎng),使兩根敏感絲溫度不相等,即ΔT≠0。設(shè)中間加熱絲的長(zhǎng)度、寬度、厚度分別為l、w、h;左右兩根細(xì)絲僅作為敏感器件長(zhǎng)度與加熱絲相同,與加熱絲的距離均為d。

    兩根敏感絲的溫度差可以通過(guò)熱電阻效應(yīng)轉(zhuǎn)變?yōu)殡娮柚抵g的差異,經(jīng)過(guò)惠斯通電橋電路分別得到兩者的阻值。質(zhì)點(diǎn)振速與敏感絲溫度差的關(guān)系為

    (1)

    式中:Ub為施加在敏感絲上的電壓值;r為敏感絲材料的溫度系數(shù);v為腔體內(nèi)質(zhì)點(diǎn)振速;ΔT(v)為敏感絲的溫度差,是關(guān)于質(zhì)點(diǎn)振速v的函數(shù)。經(jīng)過(guò)相關(guān)計(jì)算與比較,后續(xù)的討論中將微通道內(nèi)總溫度場(chǎng)拆分為穩(wěn)態(tài)溫度場(chǎng)和熱擾動(dòng)場(chǎng)的耦合。

    1.1 穩(wěn)態(tài)溫度場(chǎng)模型

    如前所述,溫度場(chǎng)是由中間的加熱絲產(chǎn)生的。通道中氣體的低速運(yùn)動(dòng)可以近似為層流,其粒子沿平行于通道的直線平滑移動(dòng)。因此,從加熱絲到空氣的傳熱模式是基于熱傳導(dǎo)方式。有熱源的熱方程[9]可以用來(lái)描述這種導(dǎo)熱模式:

    (2)

    式中:U為腔體內(nèi)各位置的溫度分布;α為熱擴(kuò)散系數(shù);F(x,y,z,t)為與時(shí)間相關(guān)的腔體內(nèi)的熱量分布函數(shù);K為空氣的導(dǎo)熱系數(shù);p為空氣的密度;s為空氣的比熱;Q為腔體內(nèi)每單位體積的熱能。通過(guò)變量分離,可得:

    U(x,y,z,t)=T(x,y,z)L(t)

    (3)

    式中:T(x,y,z)為與時(shí)間無(wú)關(guān)的穩(wěn)態(tài)溫度分布;L(t)為溫度場(chǎng)與時(shí)間相關(guān)的變化幅度。假設(shè)溫度場(chǎng)在距離加熱絲l/2點(diǎn)處的衰減為0,因此,穩(wěn)態(tài)熱傳導(dǎo)方程可以簡(jiǎn)化為

    -k2T=Q

    (4)

    假設(shè)加熱絲在y方向的長(zhǎng)度為無(wú)限長(zhǎng)(l?w,l?h)且在x和z方向上無(wú)限薄。然后,可以近似認(rèn)為溫度場(chǎng)沿y軸沒(méi)有發(fā)生變化。因此,可以化簡(jiǎn)并省略三維解析模型中的y方向溫度場(chǎng)的分析,只考慮溫度在x和z方向的分布,從而將三維模型轉(zhuǎn)化為二維模型。假設(shè)施加在加熱絲的電功率為P,則單位長(zhǎng)度的功率為P/l,那么加熱線每體積產(chǎn)生的熱量為

    (5)

    式中:δ(x)和δ(z)為沖激函數(shù),表示加熱絲被近似認(rèn)為沿y方向沒(méi)有厚度和寬度。將式(5)代入式(4)可得:

    (6)

    將式(6)轉(zhuǎn)化為極坐標(biāo)并進(jìn)行積分運(yùn)算,則其通解為

    (7)

    (8)

    從而將三維解析模型化簡(jiǎn)為一維解析模型。然而,由于在上述分析中忽略了加熱絲的寬度,當(dāng)x接近零時(shí),T(x)不能完全描述從-w/2到w/2的溫度分布??梢酝ㄟ^(guò)將寬度視為連續(xù)點(diǎn)源從-w/2到w/2的連續(xù)對(duì)稱分布來(lái)計(jì)算加熱絲的平均溫度。因此,式(8)可通過(guò)設(shè)置x=0來(lái)計(jì)算平均值:

    (9)

    PVS幾何參數(shù)與物理常數(shù)如表1所示。將表1中參數(shù)代入式(9)可得加熱絲上溫度約為621.8 K,再經(jīng)式(8)計(jì)算,可得三絲型PVS微通道內(nèi)穩(wěn)態(tài)溫度場(chǎng)分布,如圖3所示。從圖3中可以看到圖像左右對(duì)稱,表明在沒(méi)有聲源情況下,兩根敏感絲的溫度差為0。

    表1 PVS幾何參數(shù)與物理常數(shù)

    圖3 PVS穩(wěn)態(tài)溫度場(chǎng)

    1.2 動(dòng)態(tài)溫度場(chǎng)模型

    當(dāng)聲波進(jìn)入傳感器的微通道時(shí),靠近聲源的粒子將開始振動(dòng),導(dǎo)致氣體在通道中移動(dòng)。

    對(duì)于粒子的這種運(yùn)動(dòng),假設(shè)總溫度場(chǎng)是穩(wěn)態(tài)溫度場(chǎng)和熱擾動(dòng)場(chǎng)的疊加。溫度可以從熱方程中求出

    Ps(?tTe+vTe)-K2Te=Q

    (10)

    式中:Te為當(dāng)聲波進(jìn)入微通道后的總溫度場(chǎng);v為質(zhì)點(diǎn)振動(dòng)的速度;vTe為由進(jìn)入微通道聲波所帶來(lái)的強(qiáng)迫對(duì)流運(yùn)動(dòng);Q為由式(5)定義的熱功率密度。

    當(dāng)比較熱擴(kuò)散過(guò)程和強(qiáng)迫對(duì)流運(yùn)動(dòng)時(shí),由于一個(gè)粒子從一根絲擴(kuò)散到另一根絲的時(shí)間為l2/D,而強(qiáng)迫對(duì)流運(yùn)動(dòng)的時(shí)間為l/v。比較這兩次運(yùn)動(dòng)的時(shí)間,可得

    (11)

    式中:熱擴(kuò)散系數(shù)D=k/ps≈1.9×10-5m2/s(空氣);強(qiáng)迫對(duì)流的速度v=4.4 mm/s(對(duì)應(yīng)100 dB下的較高的聲壓)。由于熱擴(kuò)散速度D/l≈0.2 m/s遠(yuǎn)大于質(zhì)點(diǎn)振速v。假設(shè)溫度為T+δT,其中T是穩(wěn)態(tài)情況下的溫度分布,δT是對(duì)流引起的校正。因此,用T+δT替換Te,將公式(9)轉(zhuǎn)換為

    (12)

    從前面的分析可知,三絲型PVS的物理模型已簡(jiǎn)化為一維解析模型,即僅為x方向的溫度分布??芍?tT=0,vT=v?xT,vδT=v?xδT,-K2T=Q。根據(jù)Mailly等[10]的實(shí)驗(yàn),x方向的溫度分布梯度遠(yuǎn)大于擾動(dòng)梯度。因此,位于x方向的項(xiàng)?xδT,在方程(12)中可以忽略,從而得到:

    (13)

    在聲學(xué)分析中,聲波的振動(dòng)是周期性變化的,聲波引起的質(zhì)點(diǎn)振速可定義為

    v=v0ei(2πft+φ0)

    (14)

    式中:f為聲波的頻率;φ0為聲波的初始相位,則溫度校正可以被表示為

    (15)

    將式(7)、式(14)和式(15)代入式(13),式(13)可寫為

    (16)

    當(dāng)x和z趨近于0時(shí),式(16)無(wú)法直接求解。因此,可以通過(guò)z坐標(biāo)上對(duì)該式進(jìn)行傅里葉變換來(lái)求解

    (17)

    式中:q是沿z方向的溫度梯度,通過(guò)將式(17)代入式(16),表達(dá)式可以轉(zhuǎn)換為

    (18)

    (19)

    式中:q應(yīng)解釋為絕對(duì)值|q|。此時(shí),式(18)的解在x=0時(shí)是連續(xù)的,并一直衰減,其表達(dá)式為

    (20)

    通過(guò)對(duì)式(20)進(jìn)行傅里葉逆變換并用貝塞爾函數(shù)求解,可以得到簡(jiǎn)化的熱擾動(dòng)場(chǎng)

    (21)

    當(dāng)聲場(chǎng)振動(dòng)頻率為 300 Hz、聲壓級(jí)為 94 dB、質(zhì)點(diǎn)振速為v0=2 mm/s時(shí),在式(21)中代入表1中數(shù)據(jù),可得到沿微通道中心方向上的熱擾動(dòng)場(chǎng)分布,如圖4所示,聲波從微通道一側(cè)(x負(fù)方向)進(jìn)入腔體內(nèi),引起的熱擾動(dòng)場(chǎng)在靠近通道入口處腔體一側(cè)為負(fù),靠近出口一側(cè)為正,從而使兩根敏感絲溫度值不相等。

    圖4 PVS熱擾動(dòng)場(chǎng)

    總溫度場(chǎng)是穩(wěn)態(tài)溫度場(chǎng)和熱擾動(dòng)場(chǎng)的疊加,可以表示為

    (22)

    上述方程為三絲型PVS的基本解析模型,傳感器墻體內(nèi)穩(wěn)態(tài)溫度場(chǎng)于聲波傳來(lái)時(shí)的總溫度場(chǎng)分別如圖5藍(lán)色實(shí)線與紅色虛線所示,證明了當(dāng)聲波引起空氣中質(zhì)點(diǎn)發(fā)生振動(dòng)使靜態(tài)溫度場(chǎng)與熱擾動(dòng)場(chǎng)耦合時(shí),使兩根敏感絲溫度不再相等,產(chǎn)生溫度差。

    圖5 PVS總溫度場(chǎng)

    2 三絲型PVS的有限元分析

    為了驗(yàn)證上述解析模型的準(zhǔn)確性和探討解析模型與實(shí)際器件性能上的差別,本文采用數(shù)值模擬的方法對(duì)矩形通道內(nèi)單熱源自然對(duì)流和強(qiáng)制對(duì)流進(jìn)行了仿真以貼近實(shí)際傳感器腔內(nèi)換熱情況。該模型采用有限元的方法(Finite Element Method,FEM)[11]將PVS的腔體、流體介質(zhì)和熱阻絲等沿坐標(biāo)軸方向劃分為有限數(shù)量控制單元。本文利用有限元分析軟件ANSYS中的Fluent模塊對(duì)上述問(wèn)題進(jìn)行求解。

    傳感器簡(jiǎn)化三維模型的微通道如圖6所示,其長(zhǎng)、寬、高分別為0.6 mm、0.6 mm、0.3 mm,兩根敏感絲與中間加熱絲的距離均為0.08 mm。當(dāng)聲波引起微通道入口處質(zhì)點(diǎn)發(fā)生振動(dòng)后,會(huì)使腔體內(nèi)部空氣發(fā)生受迫對(duì)流運(yùn)動(dòng),并在出口處流出。之后,對(duì)該模型進(jìn)行網(wǎng)格劃分如圖7所示,在三根熱阻絲附近采用了更高精度的網(wǎng)格劃分以提高仿真的精度,本次共選用15303個(gè)網(wǎng)格,兼具了簡(jiǎn)潔與準(zhǔn)確性。

    圖6 PVS數(shù)值仿真模型

    圖7 PVS數(shù)值仿真網(wǎng)格節(jié)點(diǎn)

    加熱絲和感應(yīng)絲以鉑為材料且密度和比熱為21450 kg/m3和130 J/kg·K。對(duì)于密度和比熱,空氣在流體域的邊界條件分別為1.225 kg/m3和1006.43 J/kg·K。以上述條件為基礎(chǔ),ANSYS Fluent進(jìn)行了如圖8所示的穩(wěn)態(tài)溫度場(chǎng)仿真。其結(jié)果表明在沒(méi)有熱擾動(dòng)的情況下,加熱絲產(chǎn)生的熱量均勻且逐漸地分布到兩側(cè),使兩敏感絲具有相同的溫度。

    圖8 PVS數(shù)值仿真穩(wěn)態(tài)溫度場(chǎng)

    為了驗(yàn)證這兩個(gè)部分的準(zhǔn)確性以及可靠性,數(shù)值仿真與解析模型進(jìn)行對(duì)比驗(yàn)證如圖9所示。在穩(wěn)態(tài)溫度場(chǎng)中,數(shù)值仿真與解析模型中間絲的最高溫度分別為621.0 K和621.8 K,這顯示了解析模型體現(xiàn)了與理論值極相似的最高溫度值,兩者間只有0.13%的誤差。在有溫度擾動(dòng)的情況下,入口速度為 2 mm/s時(shí),以模擬緩慢穩(wěn)定的流動(dòng),數(shù)值仿真與解析模型有少許變動(dòng)如圖 10 所示,兩感應(yīng)絲之間會(huì)有溫度偏差。模型與理論值,中間絲的最高溫度分別為620.2 K和621.4 K。與穩(wěn)態(tài)溫度場(chǎng)相比,中間絲的最高溫度,在數(shù)值仿真與解析模型中分別減少了0.8 K和0.4 K。數(shù)值仿真與解析模型具有相似的溫度趨勢(shì),只是當(dāng)中解析模型的整體溫度比數(shù)值仿真的整體溫度較低。

    圖9 穩(wěn)態(tài)溫度場(chǎng)PVS數(shù)值仿真與解析模型驗(yàn)證

    圖10 熱擾動(dòng)場(chǎng)PVS數(shù)值仿真與解析模型驗(yàn)證

    3 結(jié)束語(yǔ)

    本文從熱傳導(dǎo)方程出發(fā),初步建立了三絲型PVS的穩(wěn)態(tài)溫度場(chǎng)模型和熱擾動(dòng)場(chǎng)模型,利用ANSYS 軟件建立三維模型,對(duì)三絲型PVS進(jìn)行了熱傳導(dǎo)與層流運(yùn)動(dòng)耦合的數(shù)值仿真,并驗(yàn)證解析模型的準(zhǔn)確性,結(jié)果表明數(shù)值模擬和解析模型均具有很好的可靠性。

    猜你喜歡
    質(zhì)點(diǎn)聲波穩(wěn)態(tài)
    可變速抽水蓄能機(jī)組穩(wěn)態(tài)運(yùn)行特性研究
    碳化硅復(fù)合包殼穩(wěn)態(tài)應(yīng)力與失效概率分析
    巧用“搬運(yùn)法”解決連續(xù)質(zhì)點(diǎn)模型的做功問(wèn)題
    電廠熱力系統(tǒng)穩(wěn)態(tài)仿真軟件開發(fā)
    煤氣與熱力(2021年4期)2021-06-09 06:16:54
    元中期歷史劇對(duì)社會(huì)穩(wěn)態(tài)的皈依與維護(hù)
    中華戲曲(2020年1期)2020-02-12 02:28:18
    愛(ài)的聲波 將愛(ài)留在她身邊
    聲波殺手
    質(zhì)點(diǎn)的直線運(yùn)動(dòng)
    質(zhì)點(diǎn)的直線運(yùn)動(dòng)
    自適應(yīng)BPSK在井下鉆柱聲波傳輸中的應(yīng)用
    成人免费观看视频高清| 久久精品亚洲av国产电影网| 五月伊人婷婷丁香| 日韩精品有码人妻一区| av免费在线看不卡| 男女高潮啪啪啪动态图| 色视频在线一区二区三区| 亚洲欧美一区二区三区久久| 韩国av在线不卡| 午夜精品国产一区二区电影| 久久热在线av| 亚洲精品日本国产第一区| 亚洲欧美日韩另类电影网站| 国产爽快片一区二区三区| 久久精品久久久久久久性| 国产精品熟女久久久久浪| 国产成人精品福利久久| 日韩精品免费视频一区二区三区| 久久99精品国语久久久| 日韩中文字幕视频在线看片| 国产精品三级大全| 免费观看无遮挡的男女| 午夜福利影视在线免费观看| 中文字幕精品免费在线观看视频| 免费在线观看视频国产中文字幕亚洲 | 女人久久www免费人成看片| 一级片免费观看大全| 亚洲精品乱久久久久久| 久久午夜福利片| 欧美成人精品欧美一级黄| 国产又色又爽无遮挡免| 日韩人妻精品一区2区三区| 午夜激情久久久久久久| 一本久久精品| 久久国产精品男人的天堂亚洲| 女性生殖器流出的白浆| 美女主播在线视频| 黄片播放在线免费| 人妻系列 视频| 精品人妻在线不人妻| 高清黄色对白视频在线免费看| 久久精品国产自在天天线| 我的亚洲天堂| 亚洲精品国产av蜜桃| 十八禁网站网址无遮挡| 免费观看av网站的网址| 亚洲成人av在线免费| 一本久久精品| 久久久久网色| 如何舔出高潮| 国产又爽黄色视频| 亚洲在久久综合| 在线看a的网站| 男女边吃奶边做爰视频| 亚洲av中文av极速乱| 天堂中文最新版在线下载| 日韩一本色道免费dvd| 一区在线观看完整版| 久久这里只有精品19| 中文精品一卡2卡3卡4更新| 制服人妻中文乱码| 美女福利国产在线| 人妻人人澡人人爽人人| 国产黄色视频一区二区在线观看| 亚洲欧美清纯卡通| 韩国高清视频一区二区三区| 成年动漫av网址| 亚洲av日韩在线播放| 成年女人毛片免费观看观看9 | 狠狠婷婷综合久久久久久88av| 99热国产这里只有精品6| 在线天堂中文资源库| 一区在线观看完整版| 亚洲,一卡二卡三卡| 欧美成人精品欧美一级黄| 精品人妻熟女毛片av久久网站| 2022亚洲国产成人精品| 精品一区在线观看国产| 免费黄色在线免费观看| 日韩在线高清观看一区二区三区| 国产一区二区三区av在线| xxxhd国产人妻xxx| 亚洲国产欧美日韩在线播放| 国产精品一区二区在线观看99| 高清黄色对白视频在线免费看| 99久久中文字幕三级久久日本| 99久久综合免费| 高清在线视频一区二区三区| 午夜福利网站1000一区二区三区| 2018国产大陆天天弄谢| 日日爽夜夜爽网站| 日本欧美国产在线视频| 九九爱精品视频在线观看| 美女中出高潮动态图| 亚洲精品aⅴ在线观看| 色吧在线观看| 国产黄色视频一区二区在线观看| 亚洲男人天堂网一区| 日本午夜av视频| 一区二区三区四区激情视频| 国产日韩欧美亚洲二区| 肉色欧美久久久久久久蜜桃| 在线观看国产h片| 国产乱来视频区| 亚洲精品一二三| 黄片无遮挡物在线观看| 日韩成人av中文字幕在线观看| 一二三四在线观看免费中文在| 男人添女人高潮全过程视频| 久久精品亚洲av国产电影网| 热re99久久国产66热| 国产亚洲精品第一综合不卡| 亚洲av国产av综合av卡| 两个人看的免费小视频| 亚洲精品中文字幕在线视频| 精品国产乱码久久久久久男人| 精品久久久久久电影网| av电影中文网址| 欧美在线黄色| 日本色播在线视频| 免费观看无遮挡的男女| 亚洲美女搞黄在线观看| 久久久久久久久久人人人人人人| 亚洲av福利一区| 欧美日韩av久久| 熟女av电影| 国产精品成人在线| 国产精品一区二区在线不卡| 国产亚洲最大av| 最近最新中文字幕大全免费视频 | 日韩成人av中文字幕在线观看| 免费在线观看视频国产中文字幕亚洲 | 91成人精品电影| h视频一区二区三区| 亚洲av电影在线进入| 黄片无遮挡物在线观看| 视频在线观看一区二区三区| 国产一区二区在线观看av| 天堂8中文在线网| 在线观看免费高清a一片| 麻豆精品久久久久久蜜桃| 老司机影院成人| 妹子高潮喷水视频| 少妇被粗大的猛进出69影院| 看非洲黑人一级黄片| 国语对白做爰xxxⅹ性视频网站| 国产成人av激情在线播放| 999久久久国产精品视频| 久久久久久久久久人人人人人人| 十分钟在线观看高清视频www| 欧美人与性动交α欧美软件| 日韩在线高清观看一区二区三区| 美女主播在线视频| videossex国产| 男女无遮挡免费网站观看| 超色免费av| 国产成人免费无遮挡视频| videos熟女内射| 最近最新中文字幕大全免费视频 | 99香蕉大伊视频| 午夜福利影视在线免费观看| 国产探花极品一区二区| 性少妇av在线| 国产精品国产三级国产专区5o| 国产精品不卡视频一区二区| 久久这里只有精品19| 精品卡一卡二卡四卡免费| 久久久久久久大尺度免费视频| 亚洲在久久综合| 亚洲精品久久久久久婷婷小说| 国产成人91sexporn| av女优亚洲男人天堂| 亚洲在久久综合| 日韩大片免费观看网站| 欧美精品亚洲一区二区| 久久婷婷青草| 美女国产视频在线观看| 少妇人妻精品综合一区二区| 老司机影院毛片| 美女国产视频在线观看| 日韩 亚洲 欧美在线| 99国产精品免费福利视频| 亚洲一级一片aⅴ在线观看| 久久精品夜色国产| 亚洲精品美女久久av网站| 一级毛片 在线播放| 亚洲精品成人av观看孕妇| 性高湖久久久久久久久免费观看| 少妇人妻久久综合中文| 黄色视频在线播放观看不卡| 国产精品av久久久久免费| 亚洲精品国产色婷婷电影| 狠狠精品人妻久久久久久综合| 亚洲国产精品999| 成人国产麻豆网| 亚洲av免费高清在线观看| 久久亚洲国产成人精品v| 亚洲综合色惰| 女人被躁到高潮嗷嗷叫费观| 成人毛片a级毛片在线播放| 最近最新中文字幕免费大全7| 性色avwww在线观看| h视频一区二区三区| 国产熟女午夜一区二区三区| 97精品久久久久久久久久精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 丰满少妇做爰视频| 国产午夜精品一二区理论片| 日日撸夜夜添| 欧美 日韩 精品 国产| 男女下面插进去视频免费观看| 中文字幕最新亚洲高清| av线在线观看网站| 男女啪啪激烈高潮av片| 男男h啪啪无遮挡| 另类亚洲欧美激情| 亚洲国产av影院在线观看| 欧美精品亚洲一区二区| 人人妻人人澡人人爽人人夜夜| 国产一区亚洲一区在线观看| 久久久久久人人人人人| 美女大奶头黄色视频| 亚洲,一卡二卡三卡| h视频一区二区三区| 国产熟女午夜一区二区三区| 一边摸一边做爽爽视频免费| 亚洲成人av在线免费| 国产日韩一区二区三区精品不卡| 日日啪夜夜爽| 国产黄频视频在线观看| 美国免费a级毛片| 午夜福利,免费看| 国产精品人妻久久久影院| 欧美亚洲日本最大视频资源| 18在线观看网站| 亚洲国产最新在线播放| 深夜精品福利| 免费av中文字幕在线| 日日啪夜夜爽| 汤姆久久久久久久影院中文字幕| 综合色丁香网| av视频免费观看在线观看| 国产精品秋霞免费鲁丝片| 一边摸一边做爽爽视频免费| 伊人亚洲综合成人网| 毛片一级片免费看久久久久| 高清视频免费观看一区二区| 欧美国产精品一级二级三级| 国产精品久久久av美女十八| 日本wwww免费看| 人妻少妇偷人精品九色| 国产在线视频一区二区| 欧美少妇被猛烈插入视频| 男人舔女人的私密视频| 欧美日韩成人在线一区二区| 亚洲国产色片| 国产精品av久久久久免费| 黄色怎么调成土黄色| 日韩电影二区| 一本—道久久a久久精品蜜桃钙片| 国产精品熟女久久久久浪| 国产精品.久久久| 欧美国产精品一级二级三级| 亚洲国产看品久久| 视频区图区小说| 久久精品夜色国产| 熟女电影av网| 亚洲成人手机| 看十八女毛片水多多多| 国产综合精华液| h视频一区二区三区| 性少妇av在线| 亚洲av在线观看美女高潮| 香蕉国产在线看| 亚洲国产色片| 国产成人精品无人区| 9191精品国产免费久久| 熟女av电影| 亚洲色图 男人天堂 中文字幕| 女人久久www免费人成看片| 日韩欧美一区视频在线观看| 国产白丝娇喘喷水9色精品| 国产精品一国产av| 午夜老司机福利剧场| 精品一品国产午夜福利视频| 青草久久国产| 午夜91福利影院| 精品一品国产午夜福利视频| 亚洲av日韩在线播放| 国产精品一区二区在线不卡| 丁香六月天网| 久久久久久久久久久久大奶| 久久久久久久精品精品| 人妻系列 视频| 国产白丝娇喘喷水9色精品| 久久久国产一区二区| 女人被躁到高潮嗷嗷叫费观| 成年人免费黄色播放视频| 亚洲av福利一区| 女的被弄到高潮叫床怎么办| 满18在线观看网站| 久久午夜综合久久蜜桃| 欧美av亚洲av综合av国产av | 色播在线永久视频| 国产xxxxx性猛交| 婷婷色综合大香蕉| 性色avwww在线观看| 两个人免费观看高清视频| 丝袜人妻中文字幕| 久久久久久久精品精品| tube8黄色片| 高清黄色对白视频在线免费看| 91午夜精品亚洲一区二区三区| 国产精品 欧美亚洲| 欧美精品高潮呻吟av久久| 精品一区二区三区四区五区乱码 | 精品久久蜜臀av无| 久久精品国产a三级三级三级| 成年人免费黄色播放视频| 国产精品女同一区二区软件| 1024香蕉在线观看| 有码 亚洲区| av在线老鸭窝| 爱豆传媒免费全集在线观看| 国产深夜福利视频在线观看| 一区二区av电影网| 男人舔女人的私密视频| 26uuu在线亚洲综合色| 国产无遮挡羞羞视频在线观看| 亚洲一码二码三码区别大吗| 欧美国产精品va在线观看不卡| videosex国产| 好男人视频免费观看在线| 欧美激情极品国产一区二区三区| 久久午夜综合久久蜜桃| 久久毛片免费看一区二区三区| 99久久中文字幕三级久久日本| 欧美亚洲日本最大视频资源| 中文字幕人妻熟女乱码| 蜜桃国产av成人99| 只有这里有精品99| 日韩一区二区三区影片| 汤姆久久久久久久影院中文字幕| 97精品久久久久久久久久精品| 天堂中文最新版在线下载| 一本久久精品| 久久久久久久久久人人人人人人| 精品少妇内射三级| 成人亚洲精品一区在线观看| 亚洲国产欧美网| 国产精品.久久久| 亚洲色图综合在线观看| 亚洲精品美女久久av网站| 精品少妇久久久久久888优播| 一区福利在线观看| www.精华液| 亚洲精品久久成人aⅴ小说| 可以免费在线观看a视频的电影网站 | 久久鲁丝午夜福利片| 搡老乐熟女国产| 精品人妻在线不人妻| www.精华液| 日韩电影二区| 欧美中文综合在线视频| 成年女人在线观看亚洲视频| 国产亚洲欧美精品永久| 美女中出高潮动态图| 欧美日韩亚洲国产一区二区在线观看 | 精品一区二区免费观看| 亚洲精品国产av蜜桃| 欧美日韩成人在线一区二区| 蜜桃国产av成人99| 五月伊人婷婷丁香| 三上悠亚av全集在线观看| 免费在线观看黄色视频的| 三上悠亚av全集在线观看| 又大又黄又爽视频免费| 老熟女久久久| 亚洲精品乱久久久久久| 天天影视国产精品| 久久久久视频综合| 在线观看一区二区三区激情| 久久久久国产一级毛片高清牌| 女的被弄到高潮叫床怎么办| 亚洲av成人精品一二三区| 亚洲图色成人| 国产 一区精品| 另类精品久久| 99香蕉大伊视频| 熟妇人妻不卡中文字幕| 一区二区av电影网| 国产爽快片一区二区三区| 久久久久精品久久久久真实原创| 老汉色∧v一级毛片| 一区在线观看完整版| 国产欧美日韩一区二区三区在线| 国产欧美日韩综合在线一区二区| 国产爽快片一区二区三区| 少妇人妻 视频| 九九爱精品视频在线观看| 国产麻豆69| 亚洲精品一区蜜桃| 少妇的丰满在线观看| 国产色婷婷99| 中文字幕最新亚洲高清| 亚洲在久久综合| 日日爽夜夜爽网站| 国产高清不卡午夜福利| 99久久精品国产国产毛片| 在线看a的网站| 欧美精品亚洲一区二区| 免费人妻精品一区二区三区视频| 亚洲成人一二三区av| 亚洲熟女精品中文字幕| 欧美成人午夜免费资源| 中国三级夫妇交换| 日韩视频在线欧美| 午夜久久久在线观看| 亚洲精华国产精华液的使用体验| 亚洲情色 制服丝袜| 欧美人与善性xxx| 日韩成人av中文字幕在线观看| 99国产精品免费福利视频| 久久 成人 亚洲| 国产av精品麻豆| av不卡在线播放| 久久久久久久久久久久大奶| 超碰97精品在线观看| 欧美人与性动交α欧美精品济南到 | 亚洲男人天堂网一区| 一区二区三区四区激情视频| 精品国产一区二区三区久久久樱花| 啦啦啦视频在线资源免费观看| 精品酒店卫生间| 一级毛片我不卡| 成年美女黄网站色视频大全免费| av不卡在线播放| 久久久久久久久久久久大奶| 亚洲成人一二三区av| 人人妻人人澡人人看| 久久久国产欧美日韩av| 夫妻午夜视频| 午夜福利,免费看| 国产成人免费无遮挡视频| 美女脱内裤让男人舔精品视频| 亚洲综合精品二区| 久久免费观看电影| 丰满饥渴人妻一区二区三| 日韩伦理黄色片| 极品少妇高潮喷水抽搐| 在线观看三级黄色| 熟女av电影| 国产免费视频播放在线视频| 99精国产麻豆久久婷婷| 丝袜脚勾引网站| 亚洲精品久久久久久婷婷小说| 久久人人爽人人片av| 精品福利永久在线观看| 老司机亚洲免费影院| 黄频高清免费视频| 久久人人爽人人片av| 男女边摸边吃奶| 又黄又粗又硬又大视频| 国产免费福利视频在线观看| 男女边摸边吃奶| 国产日韩欧美视频二区| 看非洲黑人一级黄片| 免费黄频网站在线观看国产| 国产精品欧美亚洲77777| 久久精品国产亚洲av涩爱| 亚洲国产最新在线播放| 男女边摸边吃奶| 飞空精品影院首页| 国产97色在线日韩免费| 国产成人午夜福利电影在线观看| 欧美另类一区| av有码第一页| 最近手机中文字幕大全| freevideosex欧美| av线在线观看网站| 国产亚洲欧美精品永久| 不卡视频在线观看欧美| 一区二区三区乱码不卡18| 日韩一本色道免费dvd| 欧美 日韩 精品 国产| 中国三级夫妇交换| 黄片小视频在线播放| 国产免费视频播放在线视频| 午夜福利网站1000一区二区三区| 97人妻天天添夜夜摸| 欧美精品一区二区免费开放| 亚洲国产精品国产精品| 国产男女内射视频| 欧美激情极品国产一区二区三区| 国产在视频线精品| 超色免费av| 精品99又大又爽又粗少妇毛片| av福利片在线| 精品视频人人做人人爽| 啦啦啦中文免费视频观看日本| 性高湖久久久久久久久免费观看| av电影中文网址| 久久久亚洲精品成人影院| 亚洲精品中文字幕在线视频| 中国三级夫妇交换| 亚洲精品一二三| 王馨瑶露胸无遮挡在线观看| 精品人妻偷拍中文字幕| 免费黄网站久久成人精品| 一级片'在线观看视频| 免费在线观看完整版高清| 亚洲国产最新在线播放| 国产精品国产三级专区第一集| 精品国产超薄肉色丝袜足j| 国产成人免费观看mmmm| 日韩在线高清观看一区二区三区| 看免费av毛片| 国产成人精品久久二区二区91 | 成年av动漫网址| 久久久国产一区二区| 亚洲精品美女久久av网站| 亚洲av免费高清在线观看| 国产精品无大码| 久久久久久久亚洲中文字幕| 日韩成人av中文字幕在线观看| 大陆偷拍与自拍| 国产日韩欧美视频二区| 99久久中文字幕三级久久日本| 一区二区日韩欧美中文字幕| 男男h啪啪无遮挡| 久久久亚洲精品成人影院| 久久久久久人妻| 夫妻午夜视频| 国产精品一区二区在线观看99| 中文乱码字字幕精品一区二区三区| 亚洲精品日韩在线中文字幕| 97在线人人人人妻| 成人18禁高潮啪啪吃奶动态图| 女性生殖器流出的白浆| 久久精品夜色国产| 日韩一卡2卡3卡4卡2021年| 国产不卡av网站在线观看| av国产久精品久网站免费入址| 天天操日日干夜夜撸| 精品99又大又爽又粗少妇毛片| 久久久久国产一级毛片高清牌| 18禁国产床啪视频网站| 一本大道久久a久久精品| 欧美日韩成人在线一区二区| 波多野结衣av一区二区av| 亚洲情色 制服丝袜| 成年美女黄网站色视频大全免费| 国产精品一二三区在线看| 99久久中文字幕三级久久日本| 考比视频在线观看| 精品久久久精品久久久| 亚洲美女黄色视频免费看| 少妇精品久久久久久久| 看免费av毛片| 母亲3免费完整高清在线观看 | 国产精品国产三级专区第一集| 亚洲美女搞黄在线观看| 国产成人91sexporn| 亚洲精品一二三| 又黄又粗又硬又大视频| 日本午夜av视频| 亚洲精品久久成人aⅴ小说| 99久久综合免费| 久久99精品国语久久久| 亚洲第一av免费看| 久久国产亚洲av麻豆专区| 美女大奶头黄色视频| 国产乱人偷精品视频| 久久人人爽av亚洲精品天堂| 看十八女毛片水多多多| 不卡av一区二区三区| 亚洲一区二区三区欧美精品| 久久久久久免费高清国产稀缺| 日韩av在线免费看完整版不卡| www日本在线高清视频| 美女大奶头黄色视频| 久久久久国产精品人妻一区二区| 一级片免费观看大全| 少妇熟女欧美另类| 超色免费av| 晚上一个人看的免费电影| 日韩一卡2卡3卡4卡2021年| 亚洲第一av免费看| 免费av中文字幕在线| 国产探花极品一区二区| videossex国产| 亚洲一码二码三码区别大吗| 麻豆乱淫一区二区| 男人操女人黄网站| 新久久久久国产一级毛片| 啦啦啦在线观看免费高清www| 亚洲精品在线美女| 久久精品国产自在天天线| 一个人免费看片子| 国产亚洲最大av| 亚洲中文av在线| 亚洲av免费高清在线观看| 汤姆久久久久久久影院中文字幕| 一区二区三区激情视频| 久久精品国产自在天天线| 在线观看一区二区三区激情| av天堂久久9| 色94色欧美一区二区| 老司机影院成人| 亚洲精品,欧美精品| av线在线观看网站| 色哟哟·www| 欧美激情高清一区二区三区 | 欧美变态另类bdsm刘玥| 亚洲精品国产av成人精品|