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

    LabVIEW的瞬態(tài)熱線法導(dǎo)熱系數(shù)測(cè)試系統(tǒng)的開發(fā)

    2017-01-12 10:08:14趙曉東
    關(guān)鍵詞:萬(wàn)用表熱線瞬態(tài)

    王 巖,潘 江,趙曉東

    (中國(guó)計(jì)量大學(xué) 計(jì)量測(cè)試工程學(xué)院,浙江 杭州 310018)

    LabVIEW的瞬態(tài)熱線法導(dǎo)熱系數(shù)測(cè)試系統(tǒng)的開發(fā)

    王 巖,潘 江,趙曉東

    (中國(guó)計(jì)量大學(xué) 計(jì)量測(cè)試工程學(xué)院,浙江 杭州 310018)

    研制了一種基于LabVIEW平臺(tái)的瞬態(tài)熱線法測(cè)量導(dǎo)熱系數(shù)的測(cè)試系統(tǒng),從測(cè)量原理、系統(tǒng)硬件、系統(tǒng)軟件、數(shù)據(jù)處理等方面說(shuō)明了該系統(tǒng)的設(shè)計(jì)過(guò)程.采用有限容積法對(duì)瞬態(tài)熱線法導(dǎo)熱系數(shù)的測(cè)量過(guò)程進(jìn)行數(shù)值模擬,并建立迭代算法處理實(shí)驗(yàn)數(shù)據(jù).通過(guò)LabVIEW調(diào)用動(dòng)態(tài)鏈接庫(kù)混合編程實(shí)現(xiàn)了LabVIEW平臺(tái)對(duì)Fortran語(yǔ)言的調(diào)用,集成了數(shù)據(jù)采集和處理模塊,實(shí)現(xiàn)了測(cè)試系統(tǒng)的自動(dòng)化.通過(guò)測(cè)量氮?dú)獾膶?dǎo)熱系數(shù)驗(yàn)證了該測(cè)試系統(tǒng)的可靠性.

    導(dǎo)熱系數(shù);LabVIEW平臺(tái);混合編程;數(shù)據(jù)處理;自動(dòng)化

    導(dǎo)熱系數(shù)是表征物質(zhì)傳導(dǎo)熱能力的熱物理參數(shù),在日常生活、生產(chǎn)和科學(xué)活動(dòng)中有著廣泛的應(yīng)用.瞬態(tài)熱線法導(dǎo)熱系數(shù)測(cè)量方法測(cè)試速度快,可以完全避開自然對(duì)流的影響,被公認(rèn)為是最好的液體或氣體的導(dǎo)熱系數(shù)測(cè)量方法[1].然而現(xiàn)有的瞬態(tài)熱線法實(shí)驗(yàn)數(shù)據(jù)一般通過(guò)人為選取有效數(shù)據(jù)段,進(jìn)行線性擬合處理,使測(cè)試系統(tǒng)數(shù)據(jù)采集和數(shù)據(jù)處理模塊完全分開,不利于測(cè)試系統(tǒng)的自動(dòng)化.

    本文在測(cè)試裝置和LabVIEW軟件平臺(tái)的基礎(chǔ)上,開發(fā)了基于LabVIEW的數(shù)據(jù)采集程序以控制數(shù)據(jù)采集過(guò)程,開發(fā)了數(shù)據(jù)處理程序進(jìn)行實(shí)驗(yàn)數(shù)據(jù)處理和導(dǎo)熱系數(shù)計(jì)算.使用混合編程的方式將二者合二為一,從而在測(cè)試系統(tǒng)中實(shí)現(xiàn)了數(shù)據(jù)采集、轉(zhuǎn)化、求解、顯示等功能,簡(jiǎn)化了導(dǎo)熱系數(shù)的測(cè)試過(guò)程.

    1 測(cè)試原理

    瞬態(tài)熱線法是利用測(cè)量熱絲的溫度變化引起的阻值變化來(lái)求得物質(zhì)的導(dǎo)熱系數(shù).其理想模型為[2]:在無(wú)限大的各項(xiàng)同性的流體中垂直置入一個(gè)半徑為r0,長(zhǎng)度L無(wú)限長(zhǎng)的熱絲,熱絲的熱導(dǎo)率無(wú)限大,熱容量為零.在加熱之前,流體和熱絲處于熱平衡狀態(tài),溫度都為T0.對(duì)熱源施加階躍熱流,這時(shí)熱絲及周圍流體溫就會(huì)產(chǎn)生溫升,從而求得導(dǎo)熱系數(shù).當(dāng)αt/r02很大時(shí),熱絲的溫升為ΔT(r0,t)[3]

    (1)

    其中,q為單位長(zhǎng)度熱線加熱功率,λ為導(dǎo)熱系數(shù),α為擴(kuò)散系數(shù),t為加熱開始后的時(shí)間,C為exp(λ)=1.781….

    由上式可知,只要繪制出ΔTid~lnt曲線,就可以得到導(dǎo)熱系數(shù).

    2 測(cè)試系統(tǒng)

    2.1 硬件部分

    測(cè)試系統(tǒng)硬件部分分為測(cè)試裝置和測(cè)試電路.

    瞬態(tài)熱線法導(dǎo)熱系數(shù)測(cè)試裝置主要由兩部分組成:用于存儲(chǔ)待測(cè)流體的測(cè)試環(huán)境,以及熱線傳感器.該測(cè)試裝置如圖1.

    腔體材料為具有良好耐腐蝕性的不銹鋼(0Cr17Ni12Mo2),通過(guò)法蘭擠壓聚四氟乙烯O型圈密封;底部管路用于進(jìn)液,頂部管路用于出液、抽真空以及壓力測(cè)量;引出導(dǎo)線通過(guò)特殊設(shè)計(jì)的密封組件進(jìn)行密封.測(cè)試腔體能承受壓力30 MPa,溫度253~373 K.熱線傳感器是實(shí)驗(yàn)裝置的核心部分,它既是加熱絲,又可作為電阻溫度計(jì),通常由作為熱線的極細(xì)金屬絲、張緊熱線的支撐結(jié)構(gòu)和測(cè)試導(dǎo)線三部分構(gòu)成.系統(tǒng)熱線采用英國(guó)Goodfellow生產(chǎn)的鉭絲,直徑為25 μm,純度為99.9%;直徑2 mm的鉭棒作為主要支架.因?yàn)橹毋g棒和熱線是同種材料,具有相同的熱膨脹系數(shù),所以在溫度變化時(shí)能夠保證熱線始終處于相同的應(yīng)力狀態(tài)[4].

    1—密封組件,2—法蘭,3—連接螺栓,4—導(dǎo)線,5—腔體,6—鉭棒,7—無(wú)銹鋼支撐塊,8—進(jìn)液管路,9—長(zhǎng)、短絲,10—密封圈,11—出液管路,12—手動(dòng)閥門 圖1 測(cè)試裝置結(jié)構(gòu)示意圖Figure 1 Schematic diagram of test device structure

    由于存在軸向熱傳導(dǎo),熱線端部的溫度低于中間部分,因此采取雙熱線法用長(zhǎng)絲減去短絲以消除端部效應(yīng),中間部分為工作熱線.加熱后,長(zhǎng)短絲電阻差值即為工作熱線電阻值.為了實(shí)現(xiàn)這一目的,系統(tǒng)采用惠斯通電橋測(cè)量電路,如圖2.

    圖2 數(shù)據(jù)采集系統(tǒng)硬件構(gòu)成原理圖Figure 2 Hardware structure diagram of data acquisition system

    圖中R1、R2、R3、R4為五個(gè)相同的直流電阻箱,Rst為型號(hào)為1 Ω標(biāo)準(zhǔn)電阻,RL、RS分別代表長(zhǎng)、短熱線的電阻,RL1、RS1分別代表與熱線相連接的導(dǎo)線電阻,RM為熱線中部引出的測(cè)量導(dǎo)線的電阻.

    2.2 軟件部分

    2.2.1 采集模塊

    瞬態(tài)熱線法測(cè)量是在10秒的量級(jí),要在短時(shí)間內(nèi)采集足夠多的數(shù)據(jù)點(diǎn)完成數(shù)據(jù)擬合,需要測(cè)量?jī)x器有很高的采樣速率;此外,測(cè)得電壓的變化范圍只有幾十毫伏,而準(zhǔn)確的測(cè)量電壓就需要很高的測(cè)量精度.在滿足要求的前提下,研究者們通常采用高速數(shù)據(jù)采集卡或積分式數(shù)字萬(wàn)用表來(lái)協(xié)調(diào)二者的矛盾[5].系統(tǒng)采用三塊數(shù)字萬(wàn)用表作為數(shù)字采集儀器和NI公司生產(chǎn)的數(shù)據(jù)采集卡實(shí)施數(shù)據(jù)采集,如圖2.圖中Vs為Aglient E3617A可調(diào)精密直流電源,1點(diǎn)為常閉觸點(diǎn),開關(guān)S為繼電器受控于NI PCI6251數(shù)據(jù)采集卡.

    系統(tǒng)使用LabVIEW開發(fā)數(shù)據(jù)采集程序控制整個(gè)數(shù)據(jù)采集過(guò)程,實(shí)現(xiàn)采集數(shù)據(jù)的顯示并將數(shù)據(jù)和顯示圖形以文件的形式存儲(chǔ).在測(cè)量程序的邏輯控制中先讓Agilent 34972A、34410A、34410A進(jìn)行預(yù)先數(shù)據(jù)采集,然后數(shù)據(jù)采集卡NI PCI6251觸發(fā)繼電器動(dòng)作,惠斯通(Wheatstone)電橋測(cè)量電路開始工作[6].

    計(jì)算機(jī)與數(shù)字萬(wàn)用表通過(guò)GPIB/LAN總線連接后,通過(guò)萬(wàn)用表的特定編程指令,就可對(duì)萬(wàn)用表進(jìn)行操作.程序中計(jì)算機(jī)與數(shù)字萬(wàn)用表的通信都是用LabVIEW的VISA函數(shù)來(lái)完成.在對(duì)萬(wàn)用表進(jìn)行測(cè)量配置時(shí),為了得到恒定的采樣時(shí)間間隔,需將自動(dòng)量程、自動(dòng)歸零、進(jìn)行運(yùn)算和更新顯示等會(huì)帶來(lái)測(cè)量間隔變動(dòng)的功能關(guān)閉,Agilent34410還應(yīng)采用采樣間隔代替觸發(fā)延遲[7].

    程序的前面板如圖3所示.在前面板中,可以設(shè)置采樣頻率、采樣時(shí)間和預(yù)采樣時(shí)間,以滿足不同功率下測(cè)量的需要.

    圖3 數(shù)據(jù)采集程序前面板圖Figure 3 Front panel of data acquisition program

    2.2.2 數(shù)據(jù)處理

    在瞬態(tài)熱線法導(dǎo)熱系數(shù)測(cè)量中,一般在實(shí)驗(yàn)數(shù)據(jù)的有效段截取實(shí)驗(yàn)數(shù)據(jù),然后進(jìn)行線性擬合,獲取斜率,進(jìn)而得到導(dǎo)熱系數(shù).該方法存在一些不足,如像氣體類流體則會(huì)因?qū)嶒?yàn)數(shù)據(jù)太少而導(dǎo)致擬合偏差過(guò)大等問(wèn)題.研究表明:利用有限容積法對(duì)瞬態(tài)熱線法導(dǎo)熱系數(shù)的測(cè)量過(guò)程數(shù)值模擬,并在此基礎(chǔ)上建立迭代算法對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行處理,可以優(yōu)化數(shù)據(jù)處理結(jié)果[5].

    根據(jù)這一思路,開發(fā)了基于有限容積法的數(shù)值模擬程序,并通過(guò)迭代算法使得各個(gè)時(shí)刻計(jì)算的溫升數(shù)據(jù)和實(shí)驗(yàn)的溫升數(shù)據(jù)的偏差和最小.開始時(shí)給定導(dǎo)熱系數(shù)和熱擴(kuò)散率的初始值,計(jì)算對(duì)應(yīng)時(shí)刻計(jì)算的溫升數(shù)據(jù)和實(shí)驗(yàn)溫升數(shù)據(jù)的偏差.根據(jù)偏差調(diào)整待測(cè)流體導(dǎo)熱系數(shù)和熱擴(kuò)散率的數(shù)值大小,重新計(jì)算,并重復(fù)此過(guò)程,直到計(jì)算曲線和實(shí)驗(yàn)曲線的偏差最小.此時(shí)得到的導(dǎo)熱系數(shù)即為最終的數(shù)據(jù).處理過(guò)程如圖4所示.本文采用Fortran語(yǔ)言編寫數(shù)值模擬及迭代程序.

    圖4 數(shù)據(jù)處理流程圖Figure 4 Flowing chart of data processing

    2.2.3 模塊調(diào)用

    LabVIEW平臺(tái)與Fortran語(yǔ)言語(yǔ)法兼容性不好,所以使用LabVIEW調(diào)用Fortran程序沒(méi)有調(diào)用C程序方便,但是動(dòng)態(tài)鏈接庫(kù)(dynamic link library,DLL)可以隱藏語(yǔ)言的實(shí)現(xiàn)細(xì)節(jié),而且DLL的編制與具體的編程語(yǔ)言及編譯器無(wú)關(guān)[8].只要通過(guò)調(diào)整編譯指令封裝好Fortran程序的接口,就可以實(shí)現(xiàn)像常規(guī)動(dòng)態(tài)鏈接庫(kù)的調(diào)用.

    多語(yǔ)言混合編程中,各種語(yǔ)言不能共享相同的頭文件和接口,這種調(diào)用約定的不一致是混合編程中最大的問(wèn)題.LabVIEW與Fortran語(yǔ)言混合編程也不例外,要想調(diào)用成功,也必須全面協(xié)調(diào)兩者所使用的調(diào)用約定.

    默認(rèn)情況下,F(xiàn)ortran采用引用方式傳遞所有參數(shù),在C或者STDCALL調(diào)用約定下,F(xiàn)ortran以值傳遞傳遞所有參數(shù)(數(shù)組參數(shù)除外).因此,需設(shè)置LabVIEW的CLN(call library function node,CLN)節(jié)點(diǎn)中參數(shù)傳遞屬性,F(xiàn)ortran中的值傳遞與CLN節(jié)點(diǎn)的值傳遞匹配,而Fortran中的引用傳遞則與CLN節(jié)點(diǎn)的指針傳遞相匹配.本文程序中的參數(shù)是單向傳遞的,并不需要修改,因此只需通過(guò)規(guī)定STDCALL的調(diào)用約定指定參數(shù)傳遞方式為值傳遞.具體設(shè)置如圖5.

    封裝完成后,在LabVIEW平臺(tái)上通過(guò)庫(kù)函數(shù)節(jié)點(diǎn)調(diào)用DLL并傳遞初始導(dǎo)熱系數(shù)和熱擴(kuò)散率等相關(guān)參數(shù),完成數(shù)據(jù)處理過(guò)程.

    圖5 調(diào)用過(guò)程參數(shù)設(shè)置Figure 5 Parameter setting of calling procedure

    3 功能測(cè)試

    為了測(cè)試系統(tǒng)的可靠性,現(xiàn)對(duì)溫度為276.05 K,壓力為5.057 6MPa的氮?dú)獾膶?dǎo)熱系數(shù)進(jìn)行測(cè)量及處理.實(shí)驗(yàn)中使用的加熱功率為0.109 43 W·m-1,迭代初始參數(shù)導(dǎo)熱系數(shù)為0.015 W·m-1·K-1,迭代曲線見圖6,第七次迭代曲線和實(shí)驗(yàn)數(shù)據(jù)吻合良好,迭代七次后最終導(dǎo)熱系數(shù)為0.026 87 W·m-1·K-1,實(shí)驗(yàn)數(shù)據(jù)用斜率擬合的方法得到的導(dǎo)熱系數(shù)為0.026 90W·m-1·K-1,二者的偏差為0.1%,利用REFPROP[9]計(jì)算得到的數(shù)據(jù)為0.026 773 W·m-1·K-1,與迭代計(jì)算的導(dǎo)熱系數(shù)值偏差為0.36%,與傳統(tǒng)實(shí)驗(yàn)斜率擬合得到的導(dǎo)熱系數(shù)值的偏差為0.47%,從而驗(yàn)證了測(cè)試系統(tǒng)的可行性.在計(jì)算時(shí),采用不同的導(dǎo)熱系數(shù)初始參數(shù)進(jìn)行計(jì)算,最終計(jì)算得到結(jié)果與初始參數(shù)的選擇無(wú)關(guān),與測(cè)得的實(shí)驗(yàn)數(shù)據(jù)有關(guān)[10].

    圖6 氮?dú)獾膶?shí)驗(yàn)數(shù)據(jù)結(jié)果及處理Figure 6 Experimental data and processing of nitrogen

    4 結(jié) 論

    根據(jù)瞬態(tài)熱線法的原理,研制了瞬態(tài)熱線法測(cè)量導(dǎo)熱系數(shù)的裝置.實(shí)驗(yàn)裝置采用雙熱線法來(lái)消除端部效應(yīng),并將雙熱線接入惠斯通電橋電路測(cè)量工作熱線電阻的阻值.開發(fā)了一套LabVIEW數(shù)據(jù)采集程序,通過(guò)GPIB總線和數(shù)據(jù)采集卡控制萬(wàn)用表實(shí)施數(shù)據(jù)采集.利用數(shù)值模擬技術(shù),通過(guò)不斷調(diào)整導(dǎo)熱系數(shù)的數(shù)值使得計(jì)算曲線和實(shí)驗(yàn)曲線不斷接近求得導(dǎo)熱系數(shù)值.通過(guò)把數(shù)據(jù)采集和處理結(jié)合一起,提高了系統(tǒng)測(cè)試的自動(dòng)化程度.

    [1] ASSAEL M J, ANTONIADIS K D, WAKEHAM W A. Historical evolution of the transient hot-wire technique[J]. International Journal of Thermophysics,2010,31(6):1051-1072.

    [2] HEALY J J, GROOT DE J, KESTIN J. The theory of the transient hot-wire method for measuring thermal conductivity[J]. Physica B+C,1976,82(2):392-408.

    [3] WAKEHAM W A, NAGASHIMA A, SENGERS J V. Measurement of the Transport Properties of Fluids[M]. London: Blackwell Scientific,1991:163-188.

    [4] ASSAEL M J, CHEN C F, METAXA I, et al. Thermal conductivity of suspensions of carbon nanotubes in water[J]. International Journal of Thermophysics,2004,25(4):971-985.

    [5] WOODFIELD P L, MOROE S, FUKAI J, et al. Techniques for accurate resistance measurement in the transient short-hot-wire method applied to high thermal-diffusivity gas[J]. International Journal of Thermophysics,2009,30(6):1478-1772.

    [6] 張洪武,潘江,王玉剛.瞬態(tài)熱線法液體導(dǎo)熱系數(shù)測(cè)試系統(tǒng)的研制[J].中國(guó)計(jì)量學(xué)院學(xué)報(bào),2010,21(4):283-288. ZHANG Hongwu, PAN Jiang, WANG Yugang. Development of experimental system for measurement of liquid thermal conductivity by the transient hot-wire technique[J]. Journal of China University of Metrology,2010,21(4):283-288.

    [7] Agilent Technologies Inc. Agilent 34410A Multimeter User’s Guide[Z]. Beijing: Agilent Technologies Inc,2000.

    [8] 楊樂(lè)平,李海濤,趙勇,等.LabVIEW高級(jí)程序設(shè)計(jì)[M].北京:清華大學(xué)出版社,2003:133-135.

    [9] LEMMON E W, HUBER M L, MCLINDEN M O. NIST reference fluid thermodynamic and transport properties-REFPROP[EB/OL].(2013-04-01)[2016-05-10]https://www.nist.gov/sites/default/files/documents/srd/REFPROP9.PDF.

    [10] 潘江,林娜,王玉剛.瞬態(tài)熱線法導(dǎo)熱系數(shù)測(cè)量實(shí)驗(yàn)數(shù)據(jù)處理方法的研究[J].計(jì)量學(xué)報(bào),2015,36(4):384-388. PAN Jiang, LIN Na, WANG Yugang. Data analysis for thermal conductivity measurement using transient hot-wire method[J]. Acta Metrologica Sinica,2015,36(4):384-388.

    Development of a test system for thermal conductivity with the transient hot wire method based on LabVIEW

    WANG Yan, PAN Jiang, ZHAO Xiaodong

    (College of Metrology and Measurement Engineering, China Jiliang University, Hangzhou 310018, China)

    A system for measurement of thermal conductivity based on LabVIEW platform was developed. The principle hardware and software of the system were described. The finite volume method was used to simulate the thermal conductivity of the transient hot wire method and an iterative algorithm was established to deal with the experimental data. The LabVIEW platform calling the Fortran language was achieved through the mixed-language programming of the LabVIEW calling the dynamic link library. The automation of the test system was achieved through the integration of the data acquisition and data processing modules. The performance of the test system was verified by measurement of the thermal conductivity of nitrogen.

    thermal conductivity; LabVIEW platforn; mixed-language programming; data processing; automation

    2096-2835(2016)04-0389-05

    10.3969/j.issn.2096-2835.2016.04.006

    2016-06-01 《中國(guó)計(jì)量大學(xué)學(xué)報(bào)》網(wǎng)址:zgjl.cbpt.cnki.net

    國(guó)家質(zhì)檢性行業(yè)科研專項(xiàng)(No.201410133).

    TP273

    A

    猜你喜歡
    萬(wàn)用表熱線瞬態(tài)
    96365健康熱線 有問(wèn)必答 有求必應(yīng)
    96365健康熱線 有問(wèn)必答 有求必應(yīng)
    96365健康熱線 有問(wèn)必答 有求必應(yīng)
    Passage 3 Harvest from learning to repair cars
    高壓感應(yīng)電動(dòng)機(jī)斷電重啟時(shí)的瞬態(tài)仿真
    MF14萬(wàn)用表在軌道電路測(cè)量中存在問(wèn)題的探討
    十億像素瞬態(tài)成像系統(tǒng)實(shí)時(shí)圖像拼接
    基于瞬態(tài)流場(chǎng)計(jì)算的滑動(dòng)軸承靜平衡位置求解
    DC/DC變換器中的瞬態(tài)特性分析
    福祿克推出萬(wàn)用表升級(jí)版本
    日韩一区二区视频免费看| 永久免费av网站大全| 免费人妻精品一区二区三区视频| 欧美精品国产亚洲| 成人无遮挡网站| 超碰97精品在线观看| 日韩精品免费视频一区二区三区 | √禁漫天堂资源中文www| 国产亚洲精品久久久com| 自线自在国产av| 久久久久久久久久成人| 制服诱惑二区| 丰满乱子伦码专区| 麻豆精品久久久久久蜜桃| 91成人精品电影| 日韩熟女老妇一区二区性免费视频| 亚洲精品自拍成人| 国产精品女同一区二区软件| 成人二区视频| 伦精品一区二区三区| 2022亚洲国产成人精品| 欧美性感艳星| 在线免费观看不下载黄p国产| 日韩视频在线欧美| 日韩人妻精品一区2区三区| 亚洲精品色激情综合| 大片免费播放器 马上看| 香蕉精品网在线| 午夜福利,免费看| 亚洲av欧美aⅴ国产| 自拍欧美九色日韩亚洲蝌蚪91| 咕卡用的链子| 在线观看一区二区三区激情| 一级毛片黄色毛片免费观看视频| 精品国产一区二区久久| 一区二区三区乱码不卡18| 亚洲av在线观看美女高潮| 日韩熟女老妇一区二区性免费视频| 少妇被粗大猛烈的视频| 99久久综合免费| 纯流量卡能插随身wifi吗| 精品福利永久在线观看| √禁漫天堂资源中文www| 国产精品久久久久久精品电影小说| 日本欧美视频一区| 日韩av在线免费看完整版不卡| 日韩中文字幕视频在线看片| 黄色视频在线播放观看不卡| 18禁国产床啪视频网站| 国产熟女午夜一区二区三区| 亚洲av电影在线观看一区二区三区| 51国产日韩欧美| 精品亚洲成国产av| 永久免费av网站大全| 国产亚洲欧美精品永久| 亚洲精品第二区| 两个人免费观看高清视频| 天堂中文最新版在线下载| 一级毛片 在线播放| 日韩欧美精品免费久久| 精品国产一区二区三区四区第35| 午夜日本视频在线| 久久人人爽av亚洲精品天堂| 中文字幕免费在线视频6| 精品人妻熟女毛片av久久网站| av播播在线观看一区| 免费大片18禁| 久久人人97超碰香蕉20202| www.av在线官网国产| av福利片在线| 国产在视频线精品| 久久久a久久爽久久v久久| 香蕉丝袜av| 中国三级夫妇交换| 九色亚洲精品在线播放| 国产又色又爽无遮挡免| 一级,二级,三级黄色视频| 久久ye,这里只有精品| 日韩一区二区三区影片| www.av在线官网国产| 日韩不卡一区二区三区视频在线| 午夜福利视频精品| 十八禁高潮呻吟视频| 最近中文字幕2019免费版| 99久久中文字幕三级久久日本| 亚洲精品国产av成人精品| 在线观看人妻少妇| 欧美日韩一区二区视频在线观看视频在线| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩av久久| 久久精品aⅴ一区二区三区四区 | 人成视频在线观看免费观看| 国产精品欧美亚洲77777| 亚洲五月色婷婷综合| av视频免费观看在线观看| 亚洲精品日韩在线中文字幕| 国产成人免费无遮挡视频| 国产精品久久久久久久久免| 精品国产一区二区三区久久久樱花| 成人国语在线视频| 五月天丁香电影| 精品一品国产午夜福利视频| 制服诱惑二区| 欧美人与性动交α欧美软件 | 久久ye,这里只有精品| 国产亚洲一区二区精品| 免费av不卡在线播放| 国产精品国产三级国产av玫瑰| 欧美国产精品va在线观看不卡| 国产成人精品一,二区| 午夜福利网站1000一区二区三区| 精品少妇久久久久久888优播| 亚洲国产欧美日韩在线播放| 精品熟女少妇av免费看| 国产高清国产精品国产三级| 国产精品久久久久久精品古装| 91精品三级在线观看| 99热国产这里只有精品6| 热re99久久国产66热| 乱码一卡2卡4卡精品| av.在线天堂| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产极品天堂在线| 午夜福利视频精品| 少妇的逼水好多| 亚洲精品自拍成人| 日韩欧美精品免费久久| 熟女av电影| 婷婷色综合大香蕉| 成年人午夜在线观看视频| 人妻 亚洲 视频| 超色免费av| 亚洲精品久久午夜乱码| 精品人妻一区二区三区麻豆| 亚洲精品aⅴ在线观看| 欧美另类一区| 蜜桃国产av成人99| 天堂中文最新版在线下载| 一级毛片电影观看| 97在线人人人人妻| 国产精品偷伦视频观看了| 免费观看无遮挡的男女| 日韩电影二区| 五月伊人婷婷丁香| 美女国产高潮福利片在线看| 国产一区二区三区综合在线观看 | 成人综合一区亚洲| 波多野结衣一区麻豆| 热re99久久国产66热| 国产免费又黄又爽又色| 成人黄色视频免费在线看| 热re99久久国产66热| 极品少妇高潮喷水抽搐| freevideosex欧美| 日韩欧美精品免费久久| 日韩成人伦理影院| 精品久久久久久电影网| 亚洲精品成人av观看孕妇| 久久久久国产网址| 亚洲内射少妇av| 在线观看www视频免费| 色94色欧美一区二区| 久久久久视频综合| 丝袜在线中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 国产综合精华液| 18禁观看日本| 满18在线观看网站| 免费观看无遮挡的男女| 久久久欧美国产精品| 免费在线观看黄色视频的| 欧美国产精品va在线观看不卡| 日韩在线高清观看一区二区三区| 一本—道久久a久久精品蜜桃钙片| 水蜜桃什么品种好| 最近中文字幕2019免费版| 久久久精品94久久精品| videos熟女内射| 美女大奶头黄色视频| 看免费av毛片| 婷婷成人精品国产| 香蕉丝袜av| 制服人妻中文乱码| tube8黄色片| 最新的欧美精品一区二区| 国产成人免费无遮挡视频| 日韩视频在线欧美| 日韩制服骚丝袜av| 亚洲精品成人av观看孕妇| 日韩熟女老妇一区二区性免费视频| 亚洲av福利一区| 妹子高潮喷水视频| 日本-黄色视频高清免费观看| av天堂久久9| 十八禁高潮呻吟视频| 1024视频免费在线观看| 婷婷色综合www| 亚洲,欧美精品.| 女人精品久久久久毛片| 国产精品蜜桃在线观看| 狠狠婷婷综合久久久久久88av| 中国国产av一级| av黄色大香蕉| 亚洲天堂av无毛| 国产精品久久久久成人av| 18禁国产床啪视频网站| 一级片'在线观看视频| 激情视频va一区二区三区| 精品国产一区二区三区久久久樱花| 国产亚洲最大av| 久久久久久久国产电影| 日韩一本色道免费dvd| 卡戴珊不雅视频在线播放| 久久狼人影院| 国产成人免费观看mmmm| 美女视频免费永久观看网站| 91午夜精品亚洲一区二区三区| 男男h啪啪无遮挡| 日韩欧美一区视频在线观看| 亚洲美女视频黄频| 超碰97精品在线观看| 男女免费视频国产| 国产精品久久久久久av不卡| 国产伦理片在线播放av一区| 高清av免费在线| 女人被躁到高潮嗷嗷叫费观| 免费人妻精品一区二区三区视频| 下体分泌物呈黄色| 欧美日韩综合久久久久久| 伦理电影免费视频| 青春草视频在线免费观看| 精品国产乱码久久久久久小说| 熟女人妻精品中文字幕| 欧美xxⅹ黑人| 最后的刺客免费高清国语| 国产毛片在线视频| 色哟哟·www| 亚洲美女搞黄在线观看| 香蕉国产在线看| www日本在线高清视频| 免费看光身美女| 国产精品免费大片| 日本91视频免费播放| 国产精品成人在线| 少妇 在线观看| 99视频精品全部免费 在线| 国产综合精华液| 久久狼人影院| 黑人高潮一二区| 香蕉精品网在线| 日韩av在线免费看完整版不卡| 美国免费a级毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产男女内射视频| 一个人免费看片子| 菩萨蛮人人尽说江南好唐韦庄| 建设人人有责人人尽责人人享有的| 中文字幕另类日韩欧美亚洲嫩草| 大片免费播放器 马上看| 欧美激情 高清一区二区三区| av又黄又爽大尺度在线免费看| 18禁动态无遮挡网站| 午夜日本视频在线| 97人妻天天添夜夜摸| 国产高清不卡午夜福利| 国产1区2区3区精品| 午夜视频国产福利| 久久人人爽人人爽人人片va| 少妇精品久久久久久久| 草草在线视频免费看| 性色av一级| 男女下面插进去视频免费观看 | 欧美精品人与动牲交sv欧美| 亚洲四区av| 女性生殖器流出的白浆| 啦啦啦中文免费视频观看日本| 男女无遮挡免费网站观看| 国产麻豆69| 丰满饥渴人妻一区二区三| 交换朋友夫妻互换小说| 大香蕉97超碰在线| videosex国产| 日日摸夜夜添夜夜爱| 日韩在线高清观看一区二区三区| 男女高潮啪啪啪动态图| 秋霞伦理黄片| 成人国产av品久久久| 国产亚洲精品久久久com| 五月开心婷婷网| 国产日韩一区二区三区精品不卡| 国产午夜精品一二区理论片| 在线观看三级黄色| 又黄又爽又刺激的免费视频.| 少妇被粗大猛烈的视频| 欧美成人精品欧美一级黄| 亚洲av中文av极速乱| 色婷婷久久久亚洲欧美| 狠狠婷婷综合久久久久久88av| 在线观看免费高清a一片| 免费在线观看完整版高清| 国产亚洲午夜精品一区二区久久| 精品一品国产午夜福利视频| 日韩三级伦理在线观看| videos熟女内射| 国产一区二区在线观看日韩| 久久久久久人人人人人| 国产极品天堂在线| 日韩一区二区三区影片| 亚洲av欧美aⅴ国产| 久久国内精品自在自线图片| 亚洲第一区二区三区不卡| 亚洲成人手机| 国产白丝娇喘喷水9色精品| 黄色怎么调成土黄色| 国产又色又爽无遮挡免| av天堂久久9| 考比视频在线观看| 精品一区在线观看国产| 久久久精品免费免费高清| 爱豆传媒免费全集在线观看| 一级爰片在线观看| 国产男人的电影天堂91| 亚洲国产看品久久| 色哟哟·www| 精品国产一区二区三区四区第35| 国产亚洲欧美精品永久| h视频一区二区三区| 日韩中文字幕视频在线看片| 巨乳人妻的诱惑在线观看| 久久久久久伊人网av| 高清毛片免费看| 欧美bdsm另类| 丝瓜视频免费看黄片| 五月伊人婷婷丁香| 久热久热在线精品观看| 五月伊人婷婷丁香| 18禁观看日本| 国产麻豆69| 久久久久久人妻| 亚洲天堂av无毛| 另类精品久久| 午夜av观看不卡| 国产日韩欧美视频二区| 日产精品乱码卡一卡2卡三| 97人妻天天添夜夜摸| 人成视频在线观看免费观看| 日韩 亚洲 欧美在线| 高清黄色对白视频在线免费看| 国产亚洲精品久久久com| 香蕉丝袜av| 国语对白做爰xxxⅹ性视频网站| 免费高清在线观看日韩| 免费看不卡的av| 精品久久国产蜜桃| 欧美变态另类bdsm刘玥| 成人18禁高潮啪啪吃奶动态图| 国产成人午夜福利电影在线观看| 国产视频首页在线观看| 蜜桃在线观看..| 亚洲欧洲精品一区二区精品久久久 | 国产激情久久老熟女| 国产成人精品一,二区| 日韩三级伦理在线观看| 亚洲国产看品久久| 亚洲内射少妇av| h视频一区二区三区| 午夜久久久在线观看| 日韩精品免费视频一区二区三区 | 18禁观看日本| 国产精品 国内视频| 精品少妇内射三级| 久久国产精品男人的天堂亚洲 | 夜夜骑夜夜射夜夜干| 亚洲精品久久成人aⅴ小说| 黄片无遮挡物在线观看| 嫩草影院入口| 亚洲,欧美精品.| xxx大片免费视频| 国产精品久久久久久久久免| 高清黄色对白视频在线免费看| 国产av国产精品国产| 午夜免费观看性视频| 又大又黄又爽视频免费| 一边摸一边做爽爽视频免费| 你懂的网址亚洲精品在线观看| 啦啦啦中文免费视频观看日本| 成人午夜精彩视频在线观看| 最近最新中文字幕免费大全7| 婷婷色av中文字幕| 成人毛片a级毛片在线播放| 9色porny在线观看| 国产色婷婷99| 极品人妻少妇av视频| 999精品在线视频| 久久99一区二区三区| 好男人视频免费观看在线| 免费在线观看完整版高清| 一边亲一边摸免费视频| 国产深夜福利视频在线观看| 国产一区二区在线观看日韩| 亚洲成人一二三区av| 一区二区三区精品91| 中文乱码字字幕精品一区二区三区| 女性生殖器流出的白浆| 国产高清国产精品国产三级| 日本91视频免费播放| 99香蕉大伊视频| 9热在线视频观看99| 少妇高潮的动态图| 欧美性感艳星| 日韩精品免费视频一区二区三区 | 丰满迷人的少妇在线观看| 乱码一卡2卡4卡精品| 欧美日韩精品成人综合77777| 亚洲人与动物交配视频| 免费女性裸体啪啪无遮挡网站| 国产免费视频播放在线视频| 熟女av电影| 国产精品免费大片| 老女人水多毛片| 国产熟女午夜一区二区三区| 免费黄色在线免费观看| 国产精品久久久久久久电影| 热re99久久国产66热| 免费看光身美女| 搡老乐熟女国产| 亚洲精品久久成人aⅴ小说| 色吧在线观看| 亚洲美女搞黄在线观看| 亚洲熟女精品中文字幕| 99热全是精品| 男女啪啪激烈高潮av片| 色婷婷久久久亚洲欧美| 国产av精品麻豆| 成人无遮挡网站| 寂寞人妻少妇视频99o| 黄色视频在线播放观看不卡| 亚洲美女黄色视频免费看| 在线亚洲精品国产二区图片欧美| 一本—道久久a久久精品蜜桃钙片| 男人添女人高潮全过程视频| 亚洲av电影在线观看一区二区三区| 国产一区二区三区综合在线观看 | 母亲3免费完整高清在线观看 | av线在线观看网站| 又粗又硬又长又爽又黄的视频| 国产又色又爽无遮挡免| 黄色毛片三级朝国网站| 成年美女黄网站色视频大全免费| 秋霞在线观看毛片| 免费播放大片免费观看视频在线观看| 久久久a久久爽久久v久久| 国产精品免费大片| 久久精品国产亚洲av涩爱| 好男人视频免费观看在线| 高清黄色对白视频在线免费看| 一级黄片播放器| 国产一区二区在线观看日韩| 国产男女内射视频| a级片在线免费高清观看视频| av国产精品久久久久影院| 久久女婷五月综合色啪小说| 下体分泌物呈黄色| 免费观看a级毛片全部| 亚洲国产精品国产精品| 免费久久久久久久精品成人欧美视频 | 人成视频在线观看免费观看| 国产无遮挡羞羞视频在线观看| 亚洲国产av新网站| 日韩精品免费视频一区二区三区 | 日韩欧美一区视频在线观看| 啦啦啦中文免费视频观看日本| 99国产综合亚洲精品| 久久久久精品性色| 一区在线观看完整版| 乱码一卡2卡4卡精品| 一级毛片电影观看| 女人久久www免费人成看片| 2021少妇久久久久久久久久久| 一级片免费观看大全| 国产成人欧美| 18禁裸乳无遮挡动漫免费视频| 另类亚洲欧美激情| 久久97久久精品| 成年美女黄网站色视频大全免费| 国产精品一二三区在线看| 亚洲欧美成人综合另类久久久| 天天影视国产精品| 久久综合国产亚洲精品| www.av在线官网国产| 大香蕉97超碰在线| 欧美精品一区二区免费开放| 午夜福利视频在线观看免费| 成人18禁高潮啪啪吃奶动态图| 90打野战视频偷拍视频| 国产成人aa在线观看| 最黄视频免费看| 久久久久久久亚洲中文字幕| 美女xxoo啪啪120秒动态图| 久久久国产精品麻豆| 中文乱码字字幕精品一区二区三区| 亚洲美女黄色视频免费看| 国产日韩欧美亚洲二区| 国产精品一区二区在线不卡| 波野结衣二区三区在线| 国产黄色视频一区二区在线观看| 欧美bdsm另类| 久久精品熟女亚洲av麻豆精品| xxxhd国产人妻xxx| 男女免费视频国产| 香蕉精品网在线| 国产精品人妻久久久影院| 国产伦理片在线播放av一区| 人妻人人澡人人爽人人| 国产极品天堂在线| 涩涩av久久男人的天堂| 国产日韩一区二区三区精品不卡| 国产av精品麻豆| 亚洲精品成人av观看孕妇| 精品一区二区三区四区五区乱码 | 自拍欧美九色日韩亚洲蝌蚪91| 亚洲,欧美,日韩| av电影中文网址| av在线观看视频网站免费| 97在线人人人人妻| 成人18禁高潮啪啪吃奶动态图| 免费高清在线观看日韩| 欧美成人午夜精品| 亚洲国产精品专区欧美| 国产又色又爽无遮挡免| 日韩欧美精品免费久久| 国产精品偷伦视频观看了| 国产亚洲av片在线观看秒播厂| 成人18禁高潮啪啪吃奶动态图| 99精国产麻豆久久婷婷| 国产亚洲一区二区精品| 亚洲精品日本国产第一区| 欧美人与性动交α欧美精品济南到 | 女人被躁到高潮嗷嗷叫费观| 一级片'在线观看视频| 成年人免费黄色播放视频| 久久久久网色| 国产极品粉嫩免费观看在线| 两个人免费观看高清视频| 国产精品久久久久久av不卡| 日韩av在线免费看完整版不卡| 久久精品国产亚洲av天美| 97在线视频观看| 色94色欧美一区二区| 久久久久国产精品人妻一区二区| www.色视频.com| 欧美激情 高清一区二区三区| 国产在视频线精品| 欧美日韩视频精品一区| 国产淫语在线视频| 毛片一级片免费看久久久久| 这个男人来自地球电影免费观看 | 在线亚洲精品国产二区图片欧美| 一级a做视频免费观看| 久久精品人人爽人人爽视色| 一级黄片播放器| 成人国产av品久久久| videos熟女内射| 青青草视频在线视频观看| 国产免费福利视频在线观看| 国产精品蜜桃在线观看| 国产亚洲午夜精品一区二区久久| 免费av不卡在线播放| 99热这里只有是精品在线观看| 少妇的丰满在线观看| 亚洲国产日韩一区二区| 九色亚洲精品在线播放| av片东京热男人的天堂| 国产毛片在线视频| 精品一区二区三区视频在线| 亚洲美女视频黄频| 色5月婷婷丁香| 亚洲精品日本国产第一区| 国产国语露脸激情在线看| 18+在线观看网站| 人妻人人澡人人爽人人| 欧美人与性动交α欧美精品济南到 | 丝瓜视频免费看黄片| 久久久久人妻精品一区果冻| 永久网站在线| 精品99又大又爽又粗少妇毛片| 欧美另类一区| 中文字幕另类日韩欧美亚洲嫩草| 欧美成人精品欧美一级黄| 黑人猛操日本美女一级片| 99国产综合亚洲精品| 黄色 视频免费看| 久久久精品94久久精品| 欧美日韩一区二区视频在线观看视频在线| 欧美成人精品欧美一级黄| 男女高潮啪啪啪动态图| av.在线天堂| 18禁观看日本| 成人漫画全彩无遮挡| 国产成人精品无人区| 亚洲国产精品国产精品| 天堂8中文在线网| 精品人妻偷拍中文字幕| 老女人水多毛片| 久久 成人 亚洲| av在线观看视频网站免费| 男人添女人高潮全过程视频| 国产免费一区二区三区四区乱码| 亚洲美女黄色视频免费看| 日韩一区二区视频免费看| av电影中文网址| 国产极品粉嫩免费观看在线| 久久av网站|