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

    運(yùn)用Excel計(jì)算純流體飽和熱力學(xué)性質(zhì)

    2016-12-22 07:40:26陳漓莫小梅孫美娟
    河池學(xué)院學(xué)報(bào) 2016年5期
    關(guān)鍵詞:蒸汽壓狀態(tài)方程單元格

    陳漓 莫小梅 孫美娟

    (百色學(xué)院 a.材料科學(xué)與工程學(xué)院; b.化學(xué)與環(huán)境工程學(xué)院, 廣西 百色 533000)

    ?

    運(yùn)用Excel計(jì)算純流體飽和熱力學(xué)性質(zhì)

    陳漓a莫小梅b孫美娟a

    (百色學(xué)院 a.材料科學(xué)與工程學(xué)院; b.化學(xué)與環(huán)境工程學(xué)院, 廣西 百色 533000)

    利用Excel的數(shù)組公式和相關(guān)數(shù)據(jù)處理功能,通過(guò)簡(jiǎn)單的設(shè)計(jì),解決立方型狀態(tài)方程計(jì)算流體飽和熱力學(xué)性質(zhì)所遇到的非線性高次方程、多次迭代問(wèn)題,有效避開手工計(jì)算難度大和計(jì)算機(jī)計(jì)算需編程問(wèn)題等,而且計(jì)算過(guò)程直觀、可視,利于掌握,計(jì)算結(jié)果準(zhǔn)確、客觀。該方法可為研究者、學(xué)習(xí)者和教學(xué)者提供參考。

    Excel;飽和蒸汽壓;改進(jìn)PR狀態(tài)方程;熱力學(xué)性質(zhì)

    0 引言

    計(jì)算純流體飽和熱力學(xué)性質(zhì),常常需要狀態(tài)方程,由于立方型狀態(tài)方程參數(shù)少且有較好的計(jì)算精度,因此得到廣泛的運(yùn)用。在眾多立方形方程中PR狀態(tài)方程應(yīng)用最為普遍。PR方程在預(yù)測(cè)液體體積上較之之前方程有所改善,但要進(jìn)一步外推就比較困難。DAKEWU和SHULINCHEN修正PR方程,使方程的計(jì)算精度有顯著改進(jìn),在計(jì)算純流體的飽和蒸氣壓、飽和液相密度方面有良好的準(zhǔn)確度。在此我們以改進(jìn)PR狀態(tài)方程為模型結(jié)合Excel電子表格計(jì)算純流體的飽和熱力學(xué)性質(zhì)。

    純流體處于飽和狀態(tài)時(shí),一般用4個(gè)狀態(tài)參量描述,溫度T、壓強(qiáng)p、汽相摩爾體積Vν和液相摩爾體積Vl。然而系統(tǒng)平衡時(shí)只有一個(gè)獨(dú)立變量,通常取溫度T或壓強(qiáng)p,故有兩類計(jì)算過(guò)程:第一類是蒸汽壓計(jì)算;第二類是沸點(diǎn)計(jì)算。

    以第一類的蒸氣壓計(jì)算為例說(shuō)明,并以改進(jìn)PR狀態(tài)方程為模型[1]

    (1)

    純流體在飽和狀態(tài)時(shí)需滿足汽液平衡準(zhǔn)則式:φν=φl(shuí)

    (2)

    其中φν為汽相逸度,φl(shuí)為液相逸度。

    結(jié)合狀態(tài)方程,其逸度系數(shù)表達(dá)式可寫為:

    (3)

    純流體的汽液平衡的計(jì)算,首先要根據(jù)方程輸入臨界參數(shù)和偏心因子,計(jì)算給定溫度T下的改進(jìn)PR方程常數(shù)a,b。假設(shè)p的初值,求狀態(tài)方程得到汽相摩爾體積Vν和液相摩爾體積Vl,由此判別方程式(2)是否滿足收斂條件,若不滿足,通過(guò)調(diào)節(jié)p,直到方程式(2)收斂,此時(shí)的p、Vν和Vl就是方程組式(1)和式(2)的解。由于在計(jì)算中涉及到高次方程運(yùn)算,需借助計(jì)算機(jī)進(jìn)行編程計(jì)算,操作過(guò)程較為繁瑣。如果能用辦公軟件Excel來(lái)解決,將給計(jì)算帶來(lái)很大的便利。

    Excel是微軟辦公軟件的一個(gè)重要部分,它可以進(jìn)行各種數(shù)據(jù)處理、統(tǒng)計(jì)分析等操作。對(duì)于熱力學(xué)數(shù)據(jù)的處理和運(yùn)算也具有良好的交互界面,可簡(jiǎn)單明了地展示計(jì)算原理與過(guò)程,有利于解決繁瑣的熱力學(xué)系統(tǒng)計(jì)算問(wèn)題[2-5]。

    1 純流體汽液平衡的理論基礎(chǔ)

    1.1 計(jì)算飽和蒸汽壓

    飽和蒸汽壓p進(jìn)行首輪迭代計(jì)算時(shí),要預(yù)先估計(jì)蒸汽壓的初值。運(yùn)用Ednister偏心因子估算式ω=[3Tbr/7(1-Tbr)]lgpc-1能得到臨界參數(shù)和偏心因子估計(jì)蒸汽壓的初值的計(jì)算式[6]:

    (4)

    對(duì)于飽和蒸汽壓p的迭代式,可從式(2)利用Newton迭代法得到

    p(k+1)=p(k)-[RT(lnφν-lnφ1)/(Vν-Vl)](k)

    (5)

    1.2 計(jì)算汽相、液相摩爾體積

    對(duì)于立方形狀態(tài)方程摩爾體積的計(jì)算,當(dāng)T

    Vk+1=Vk-f(Vk)/f′(Vk)

    (6)

    為了方便迭代,將改進(jìn)PR方程用摩爾體積的三次展開式來(lái)表示,即

    (7)

    f ′(V)為f(V)的一階導(dǎo)。

    汽相摩爾體積Vν的求取通常以理想氣體體積V=RT/p為初值,而液相摩爾體積以V=b為初值,代入式(6),得到Vl值后再代入等式右邊,一直迭代下去,直到滿足Vl<ε。利用Excel的計(jì)算功能,可方便同時(shí)迭代出Vν和Vl,這樣就極大提高計(jì)算的效率。

    1.3 計(jì)算偏離焓、偏離熵等熱力學(xué)性質(zhì)

    流體的各種熱力學(xué)性質(zhì),比如焓、熵值和熱容等,可以利用偏離函數(shù)并結(jié)合狀態(tài)方程計(jì)算出來(lái)[7]。

    (8)

    (9)

    (10)

    偏離定壓摩爾熱容:

    (11)

    2 在Excel中進(jìn)行純流體飽和熱力學(xué)性質(zhì)的計(jì)算

    為計(jì)算純流體汽液平衡下的飽和熱力學(xué)性質(zhì),在Excel中進(jìn)行如下操作:

    (1)以乙烷為例,分別在B3到E3單元格輸入臨界溫度、臨界壓強(qiáng)、氣體常數(shù)和偏心因子等數(shù)值。A6到F6單元格為方程常數(shù)a,b的值以及相應(yīng)一階導(dǎo)數(shù)、二階導(dǎo)數(shù)的數(shù)值。

    (2)在B10單元格輸入某一溫度“220”如圖1所示,即假設(shè)此時(shí)沸點(diǎn)溫度為220 K,運(yùn)用式(4)計(jì)算飽和蒸氣壓的初值,即在A12單元格輸入“=C3*10^(7*(1+E3)/3*(1-B3/B10))”,把計(jì)算結(jié)果以粘貼數(shù)值的方式粘貼到A10單元格里。

    (3)在單元格C9和D9分別輸入Vν的初值(V=RT/p)“=D3*B10/(A10*10^6)”和Vl的初值(V=b)“=B6”,單元格C10-C16和D10-D16分別輸入Newton迭代式(6),如C10-C16單元格輸入如下數(shù)組公式“{=C9:C15-(C9:C15^3-(D3*B10/(A10*10^6)-0.645*B6)*C9:C15^2+1/(A10*10^6)*(D6-1.645*B6*D3*B10-2.29*A10*10^6*B6^2)*C9:C15-D6*B6/(A10*10^6)+0.645*D3*B10*B6^2/(A10*10^6)+0.645*B6^3)/(3*C9:C15^2-2*(D3*B10/(A10*10^6)-0.645*B6)*C9:C15+1/(A10*10^6)*(D6-1.645*B6*D3*B10-2.29*A10*10^6*B6^2))}”。同樣D10-D16單元格中式子與C10-C16內(nèi)容相類似,所不同之處僅僅是把式中“C9:C15”改為“D9:D15”即可。一般來(lái)說(shuō)大多數(shù)情況下迭代5-6次即可得到滿意結(jié)果,如圖1所示。為確保迭代的精度我們進(jìn)行了7次迭代,單元格C16和D16得到的數(shù)值分別為汽相摩爾體積Vν和液相摩爾體積Vl,參見圖1。單元格F9和G9分別為式(7)f(V)的Vν和Vl收斂情況,f(V)應(yīng)趨于零。

    (4)接下來(lái)進(jìn)一步計(jì)算逸度系數(shù)。在單元格F11和G11分別輸入計(jì)算汽相和液相逸度的關(guān)系式(3),在單元格F11輸入“=A10*10^6*F10/(D3*B10)-1-LN(A10*10^6*(F10-B6)/(D3*B10))-D6/(2.299135707*B6*D3*B10)*LN((F10+1.9720678535*B6)/(F10-0.3270678535*B6))”。同理在單元格G11輸入液相逸度的關(guān)系式,計(jì)算后如果不滿足|lnφν-lnφl(shuí)|≤ε則需要調(diào)整蒸汽壓p進(jìn)行迭代計(jì)算,迭代關(guān)系式參見式(5),即在單元格A14輸入式(5),即“=(A10*10^6*(1-(F11-G11)/(A10*10^6*F10/(D3*B10)-A10*10^6*G10/(D3*B10))))/10^6”。把A14計(jì)算結(jié)果以粘貼數(shù)值的方式到A10單元格里,重復(fù)粘貼數(shù)值,Excel將對(duì)整個(gè)過(guò)程進(jìn)行自動(dòng)迭代計(jì)算,直至滿足|lnφν-lnφl(shuí)|≤ε,一般迭代3-4次即能得到滿意的結(jié)果,此時(shí)得到p、Vν和Vl分別為飽和蒸汽壓、汽相摩爾體積和液相摩爾體積的數(shù)值。

    (5)計(jì)算其它的熱力學(xué)性質(zhì)。在單元格F13到F16和G13到G16分別輸入式(8)、式(9)、式(10)和式(11),這樣我們同時(shí)得到乙烷汽相和液相偏離焓、偏離熵、偏離定容摩爾熱容和偏離定壓摩爾熱容等數(shù)值。例如,在單元格F13輸入“=A10*10^6*F10-D3*B10-1/(2.299135707*B6)*(D6-B10*E6)*LN((F10+1.9720678535*B6)/(F10-0.3270678535*B6))”,計(jì)算得到的數(shù)值為汽相偏離焓。

    3 應(yīng)用分析

    為了檢驗(yàn)Excel在計(jì)算純流體飽和性質(zhì)的可靠性,我們以乙烷為研究對(duì)象,對(duì)該氣體飽和蒸汽壓等性質(zhì)進(jìn)行計(jì)算。在單元格B4、B5分別輸入乙烷的臨界溫度Tc和臨界壓強(qiáng)pc,單元格B7輸入偏心因子ω的數(shù)值參見圖1,計(jì)算乙烷溫度為160~305 K的飽和蒸汽壓、摩爾體積和偏離性質(zhì)。計(jì)算的結(jié)果參見表1。

    圖1 乙烷飽和性質(zhì)的計(jì)算

    溫度T/K飽和汽壓p/MPa摩爾體積汽相Vν/m3·mol-1液相Vl/m3·mol-1偏離焓汽相[HR]ν/J·mol-1液相[HR]l/J·mol-1偏離熵汽相[SR]ν/J·(mol·K)-1液相[SR]l/J·(mol·K)-1160002126225×10-25022×10-5-289508-156550863-01076-977709180007781875×10-25238×10-5-925438-150629562-03116-834806200021567300×10-35511×10-5-2285760-144285276-07085-717083220049043363×10-35867×10-5-4753996-137265643-13765-616090240096791728×10-36356×10-5-8817424-129190432-24176-525730260172119451×10-47081×10-5-15216885-119398562-40066-440765280282805238×10-48317×10-5-25472798-106428417-65514-354641300437522585×10-41153×10-4-45660183-84250988-118712-247348305484071903×10-41503×10-4-57951513-69336618-154363-191691

    從表1計(jì)算的數(shù)據(jù)可以看出隨著溫度的增大,蒸發(fā)加劇密度增大飽和蒸汽壓升高,當(dāng)溫度升高到305 K時(shí),飽和汽壓逐步接近臨界壓強(qiáng)。同時(shí)溫度升高汽相摩爾體積Vν將隨著飽和汽壓升高而減小,液相摩爾體積Vl則隨著溫度升高而增大,在臨界點(diǎn)附近它們數(shù)值趨于一致。同樣偏離理想氣體性質(zhì)焓和熵,汽相隨飽和汽壓升高,越遠(yuǎn)離理想氣體狀態(tài),液相隨溫度升高,偏離程度減小,在臨界點(diǎn)附近它們數(shù)值趨于一致。

    把表1計(jì)算的飽和汽壓、摩爾體積、偏離焓和偏離熵結(jié)果與國(guó)內(nèi)外計(jì)算熱力學(xué)性質(zhì)軟件得到的數(shù)值是一致的,這也說(shuō)明該方程能夠較好的描述純流體的熱力學(xué)性質(zhì)。

    上面的飽和熱力學(xué)性質(zhì)的計(jì)算是以溫度T為獨(dú)立變量進(jìn)行的,同樣也可以以蒸汽壓p為獨(dú)立變量來(lái)進(jìn)行飽和熱力學(xué)性質(zhì)的計(jì)算,只需在單元格B14輸入如下的迭代式即可

    在單元格B14輸入“=B10*(1+(F11-G11)/((F13-G13)/(D3*B10)))”,多次迭代直到滿足|lnφν-lnφl(shuí)|≤ε,此時(shí)溫度T為飽和蒸汽壓沸點(diǎn)的數(shù)值。

    4 結(jié)論

    通過(guò)上面的例子可以看出,運(yùn)用Excel進(jìn)行純流體汽液兩相飽和性質(zhì)的計(jì)算,無(wú)需編程,只需在表格里輸入純流體的臨界性質(zhì)和偏心因子的數(shù)值,通過(guò)簡(jiǎn)單的迭代設(shè)計(jì),Excel即可計(jì)算各種純流體的飽和熱力學(xué)性質(zhì)。整個(gè)運(yùn)算過(guò)程有良好的交互界面,并能有效地提高計(jì)算結(jié)果的準(zhǔn)確度和可靠性。

    Exce1具有強(qiáng)大的運(yùn)算功能,但以往多用于簡(jiǎn)單計(jì)算。實(shí)際上,充分利用Exce1提供的運(yùn)算功能,只需輸入相應(yīng)的關(guān)系式,就能完成復(fù)雜的熱力學(xué)工程計(jì)算,計(jì)算過(guò)程簡(jiǎn)單、直觀,學(xué)生非常容易掌握。由于Exce1具有較強(qiáng)的數(shù)組公式和公式的復(fù)制功能,所以使用Exce1計(jì)算純流體熱力學(xué)性質(zhì)具有無(wú)可比擬的優(yōu)越性。

    [1]DAKEWU,SHULINCHEN.AmodifiedPeng-Robinsonequationonstate[J].ChemicalEngineeringCommunications, 1997, 156(1):215-225.

    [2]陳方.用EXCEL2000解決方程求根問(wèn)題[J].湖南工業(yè)職業(yè)技術(shù)學(xué)院學(xué)報(bào),2003,4(1):14-16.

    [3]班玉鳳,常圣泉,朱海峰,等.EXCEL在非理想系泡露點(diǎn)計(jì)算中的應(yīng)用[J].計(jì)算機(jī)應(yīng)用與軟件,2011,28(10):275-277.

    [4]于志家,陳傳棋,李香琴,等.應(yīng)用Excel進(jìn)行泡點(diǎn)與露點(diǎn)計(jì)算[J].化工高等教育,2012,126(4):73-76.

    [5]王智娟,胡粉娥,楊曉麗,等.Excel在化工熱力學(xué)計(jì)算中的應(yīng)用[J].廣州化工,2013,41(3):167-169.

    [6]陳鐘秀,顧飛燕,胡望明,等.化工熱力學(xué)[M].3版.北京:化學(xué)工業(yè)出版社,2012:287-288.

    [7]陳則韶.高等工程熱力學(xué)[M].2版.合肥:中國(guó)科學(xué)技術(shù)大學(xué)出版社,2014:120-125.

    [責(zé)任編輯 劉景平]

    Using Excel to Calculate the Thermodynamic Properties of Pure Fluid Saturation

    CHEN Lia, MO Xiaomeib, SUN MeiJuana

    (a.School of Materials Science and Engineering;b.School of Chemistry and Environmental Engineering, Baise University, Baise, Guangxi 533000, China)

    The study uses Excel array formulas and data processing functions, through a simple design, to solve the problems of nonlinear equation of higher degree and multiple iterations faced when calculating thermodynamic properties of fluid saturation. It can effectively avoid the problems of manual calculation being difficult and computer calculation needing programming; the calculation process is intuitive, visual, conducive to grasp, and the calculation results are accurate and objective. The method provides a reference for researchers, learners and educators.

    Excel; saturated vapor pressure; modified PR equation of state; thermodynamic properties

    O642

    A

    1672-9021(2016)05-0113-05

    陳漓(1962-),男,廣西桂林人,百色學(xué)院材料科學(xué)與工程學(xué)院副教授,主要研究方向:熱力學(xué)。

    廣西高??茖W(xué)技術(shù)研究基金資助項(xiàng)目(2013YB244)。

    2016-09-12

    猜你喜歡
    蒸汽壓狀態(tài)方程單元格
    LKP狀態(tài)方程在天然氣熱物性參數(shù)計(jì)算的應(yīng)用
    煤氣與熱力(2021年6期)2021-07-28 07:21:30
    普通玉米、糯玉米和蒸汽壓片玉米對(duì)生長(zhǎng)豬能量和營(yíng)養(yǎng)物質(zhì)消化率的影響
    玩轉(zhuǎn)方格
    玩轉(zhuǎn)方格
    蒸汽壓片玉米加工工藝及其在肉牛生產(chǎn)中應(yīng)用的研究進(jìn)展
    淺談Excel中常見統(tǒng)計(jì)個(gè)數(shù)函數(shù)的用法
    西部皮革(2018年6期)2018-05-07 06:41:07
    基于隨機(jī)與區(qū)間分析的狀態(tài)方程不確定性比較
    用狀態(tài)方程模擬氨基酸水溶液的熱力學(xué)性質(zhì)
    頁(yè)巖中甲烷虛擬飽和蒸汽壓的計(jì)算方法研究
    混合星物質(zhì)的狀態(tài)方程和奇異夸克物質(zhì)的穩(wěn)定窗
    两性午夜刺激爽爽歪歪视频在线观看 | 亚洲aⅴ乱码一区二区在线播放 | 波多野结衣一区麻豆| 午夜福利成人在线免费观看| 老司机深夜福利视频在线观看| 伊人久久大香线蕉亚洲五| 欧美激情 高清一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 欧美一级毛片孕妇| 曰老女人黄片| 久久国产精品男人的天堂亚洲| 国产精品精品国产色婷婷| 操出白浆在线播放| 91老司机精品| 女人被躁到高潮嗷嗷叫费观| 国产三级黄色录像| 免费在线观看亚洲国产| 9191精品国产免费久久| 亚洲少妇的诱惑av| 亚洲中文av在线| 久久久国产成人免费| 国产日韩一区二区三区精品不卡| 身体一侧抽搐| 琪琪午夜伦伦电影理论片6080| 好男人在线观看高清免费视频 | 午夜两性在线视频| 免费av毛片视频| 69av精品久久久久久| 可以免费在线观看a视频的电影网站| 9191精品国产免费久久| 日本五十路高清| 久久久国产精品麻豆| 午夜精品久久久久久毛片777| 久久久久久久久免费视频了| 夜夜看夜夜爽夜夜摸| 日韩欧美一区二区三区在线观看| 久久久久九九精品影院| 久久久久久久久免费视频了| 亚洲精品在线美女| 亚洲av成人av| 久久久精品国产亚洲av高清涩受| 国产精品亚洲美女久久久| 久久久久亚洲av毛片大全| 久久久久久国产a免费观看| 久久久国产成人免费| 免费观看精品视频网站| 好男人电影高清在线观看| 国产午夜精品久久久久久| 久久香蕉激情| 黑人操中国人逼视频| 国产精品久久久久久亚洲av鲁大| 制服丝袜大香蕉在线| 亚洲一区高清亚洲精品| 一本大道久久a久久精品| 亚洲美女黄片视频| 亚洲国产精品成人综合色| 色综合欧美亚洲国产小说| 亚洲va日本ⅴa欧美va伊人久久| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久久精品吃奶| 久久热在线av| 亚洲专区中文字幕在线| 免费看美女性在线毛片视频| 99久久综合精品五月天人人| 变态另类成人亚洲欧美熟女 | 国产1区2区3区精品| 国产欧美日韩一区二区三区在线| 大陆偷拍与自拍| 欧美丝袜亚洲另类 | 中国美女看黄片| 级片在线观看| 午夜精品国产一区二区电影| 天堂动漫精品| 免费不卡黄色视频| 欧美日韩精品网址| 丰满的人妻完整版| 久9热在线精品视频| 免费观看精品视频网站| 亚洲少妇的诱惑av| 中亚洲国语对白在线视频| 国产精品永久免费网站| 丝袜在线中文字幕| 叶爱在线成人免费视频播放| 色尼玛亚洲综合影院| 一级黄色大片毛片| 欧美激情 高清一区二区三区| 午夜精品在线福利| 黑人操中国人逼视频| 最好的美女福利视频网| 亚洲国产欧美一区二区综合| 人妻丰满熟妇av一区二区三区| 国产高清视频在线播放一区| 搡老妇女老女人老熟妇| 正在播放国产对白刺激| 日本黄色视频三级网站网址| 亚洲精品久久成人aⅴ小说| 19禁男女啪啪无遮挡网站| av超薄肉色丝袜交足视频| 久久草成人影院| 一边摸一边抽搐一进一出视频| 久久精品国产亚洲av香蕉五月| 1024视频免费在线观看| 窝窝影院91人妻| 日韩欧美国产一区二区入口| 欧美日韩亚洲综合一区二区三区_| 首页视频小说图片口味搜索| 在线观看免费视频日本深夜| 在线观看免费视频网站a站| 美女高潮到喷水免费观看| 久久性视频一级片| 国产精品一区二区免费欧美| 亚洲 欧美 日韩 在线 免费| 国产一卡二卡三卡精品| 老鸭窝网址在线观看| 老司机午夜福利在线观看视频| 在线国产一区二区在线| 夜夜看夜夜爽夜夜摸| 国产91精品成人一区二区三区| 视频在线观看一区二区三区| 97超级碰碰碰精品色视频在线观看| 成熟少妇高潮喷水视频| 搡老熟女国产l中国老女人| 国产成人欧美| 亚洲自拍偷在线| 国产一区二区三区综合在线观看| 久久久久久久久久久久大奶| 夜夜夜夜夜久久久久| 日韩大尺度精品在线看网址 | 黄频高清免费视频| 国产成人av教育| 国产精品一区二区免费欧美| 亚洲色图av天堂| 成人国产一区最新在线观看| 日本精品一区二区三区蜜桃| 夜夜爽天天搞| 亚洲成国产人片在线观看| 国产1区2区3区精品| 久久婷婷成人综合色麻豆| videosex国产| 欧美成人免费av一区二区三区| 精品少妇一区二区三区视频日本电影| 亚洲精品美女久久久久99蜜臀| 国产成人av教育| 欧美国产精品va在线观看不卡| 国产精品久久久人人做人人爽| 午夜福利视频1000在线观看 | 宅男免费午夜| av网站免费在线观看视频| 欧美另类亚洲清纯唯美| 久久精品国产综合久久久| 乱人伦中国视频| 黄色视频不卡| 如日韩欧美国产精品一区二区三区| 久久影院123| 一边摸一边抽搐一进一出视频| 最好的美女福利视频网| 中文字幕人妻丝袜一区二区| 狂野欧美激情性xxxx| 亚洲专区国产一区二区| 黄色视频,在线免费观看| 国产精品久久电影中文字幕| 韩国精品一区二区三区| 国产精品一区二区三区四区久久 | 久久天躁狠狠躁夜夜2o2o| 国产精品一区二区免费欧美| 国产又爽黄色视频| av福利片在线| 国产精品98久久久久久宅男小说| 亚洲av电影在线进入| 老司机福利观看| 999久久久精品免费观看国产| 一个人免费在线观看的高清视频| 久久国产精品人妻蜜桃| 亚洲少妇的诱惑av| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品综合一区在线观看 | 国产亚洲精品第一综合不卡| 成人手机av| 亚洲国产日韩欧美精品在线观看 | 丝袜人妻中文字幕| 国产精品亚洲美女久久久| av视频免费观看在线观看| 免费在线观看黄色视频的| 欧美 亚洲 国产 日韩一| 亚洲,欧美精品.| 国产精品秋霞免费鲁丝片| www日本在线高清视频| 99热只有精品国产| 久久草成人影院| 精品国产超薄肉色丝袜足j| 亚洲人成77777在线视频| 午夜成年电影在线免费观看| 免费看a级黄色片| 国产精品久久视频播放| 中文字幕精品免费在线观看视频| 黄色视频不卡| 久久影院123| 免费看a级黄色片| 丝袜人妻中文字幕| www日本在线高清视频| 十八禁网站免费在线| www国产在线视频色| 亚洲国产高清在线一区二区三 | 中亚洲国语对白在线视频| 国产1区2区3区精品| av欧美777| 99精品久久久久人妻精品| 少妇熟女aⅴ在线视频| 在线免费观看的www视频| 免费不卡黄色视频| 一进一出抽搐gif免费好疼| 51午夜福利影视在线观看| 黄色a级毛片大全视频| 日韩成人在线观看一区二区三区| 在线观看免费午夜福利视频| 黑人巨大精品欧美一区二区mp4| 99在线人妻在线中文字幕| www国产在线视频色| 欧美日韩黄片免| 男女之事视频高清在线观看| 热re99久久国产66热| 久热这里只有精品99| 精品无人区乱码1区二区| 精品不卡国产一区二区三区| 亚洲av熟女| 可以在线观看毛片的网站| 在线观看一区二区三区| 欧美日韩福利视频一区二区| 亚洲三区欧美一区| 啪啪无遮挡十八禁网站| 久久精品亚洲精品国产色婷小说| av有码第一页| 亚洲在线自拍视频| av视频免费观看在线观看| aaaaa片日本免费| 一边摸一边做爽爽视频免费| 香蕉丝袜av| 91老司机精品| 午夜福利一区二区在线看| 操美女的视频在线观看| 丝袜在线中文字幕| 亚洲五月婷婷丁香| 97人妻天天添夜夜摸| 90打野战视频偷拍视频| 中文字幕久久专区| 午夜免费观看网址| 欧美一级毛片孕妇| 欧美日本亚洲视频在线播放| 日本一区二区免费在线视频| 999久久久精品免费观看国产| 国产成人精品无人区| 国产1区2区3区精品| 精品国产美女av久久久久小说| 首页视频小说图片口味搜索| 久久欧美精品欧美久久欧美| 成年人黄色毛片网站| 亚洲午夜精品一区,二区,三区| 亚洲无线在线观看| 黑人欧美特级aaaaaa片| 两个人免费观看高清视频| 国内精品久久久久久久电影| 老熟妇仑乱视频hdxx| 欧美成人一区二区免费高清观看 | 亚洲五月婷婷丁香| 在线观看66精品国产| 日韩欧美免费精品| 午夜日韩欧美国产| 啦啦啦 在线观看视频| 伦理电影免费视频| 日韩精品免费视频一区二区三区| 欧美中文日本在线观看视频| 日本五十路高清| 欧美色视频一区免费| 777久久人妻少妇嫩草av网站| 成在线人永久免费视频| 热re99久久国产66热| 深夜精品福利| www.熟女人妻精品国产| avwww免费| av片东京热男人的天堂| 97碰自拍视频| 国产麻豆成人av免费视频| 久久精品91蜜桃| 久久亚洲真实| 亚洲午夜精品一区,二区,三区| 国产精品精品国产色婷婷| 久久青草综合色| 国产免费av片在线观看野外av| 久久久国产成人精品二区| 999久久久国产精品视频| 国产亚洲av嫩草精品影院| 日韩精品中文字幕看吧| 黄频高清免费视频| 中文字幕精品免费在线观看视频| 成人亚洲精品一区在线观看| 国产成人免费无遮挡视频| 岛国视频午夜一区免费看| 熟女少妇亚洲综合色aaa.| 亚洲美女黄片视频| 最新在线观看一区二区三区| 欧美日本亚洲视频在线播放| 两人在一起打扑克的视频| 亚洲一区二区三区色噜噜| 久久精品91无色码中文字幕| 欧美日本视频| 午夜a级毛片| 久久性视频一级片| 国产成人av激情在线播放| 国产av在哪里看| 久久亚洲真实| 日日摸夜夜添夜夜添小说| 久久狼人影院| 99久久国产精品久久久| 侵犯人妻中文字幕一二三四区| 国产av精品麻豆| 天堂影院成人在线观看| 精品人妻1区二区| 日本 欧美在线| 色综合婷婷激情| 亚洲成人国产一区在线观看| 欧美丝袜亚洲另类 | 男人舔女人的私密视频| 人成视频在线观看免费观看| 脱女人内裤的视频| 国产不卡一卡二| 亚洲一卡2卡3卡4卡5卡精品中文| av电影中文网址| 国产精品,欧美在线| 免费在线观看黄色视频的| 欧美国产精品va在线观看不卡| 免费在线观看完整版高清| 99久久精品国产亚洲精品| 在线观看66精品国产| videosex国产| 法律面前人人平等表现在哪些方面| 久久婷婷成人综合色麻豆| 免费观看人在逋| xxx96com| 欧美av亚洲av综合av国产av| 一本久久中文字幕| 亚洲国产欧美网| 丰满的人妻完整版| 久久这里只有精品19| 一级a爱视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 精品久久久久久成人av| www.自偷自拍.com| x7x7x7水蜜桃| 人人妻人人澡欧美一区二区 | av欧美777| 成熟少妇高潮喷水视频| 一区二区日韩欧美中文字幕| 在线播放国产精品三级| 国产一卡二卡三卡精品| 亚洲在线自拍视频| 在线天堂中文资源库| 中文字幕久久专区| 女性生殖器流出的白浆| 亚洲午夜理论影院| 欧美激情高清一区二区三区| 757午夜福利合集在线观看| 黄频高清免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲专区中文字幕在线| 精品国产一区二区久久| 欧美日本亚洲视频在线播放| 久久久久国产一级毛片高清牌| 久久国产乱子伦精品免费另类| 国产97色在线日韩免费| 亚洲第一av免费看| 久久伊人香网站| 黄色女人牲交| 精品第一国产精品| 变态另类丝袜制服| 精品高清国产在线一区| 久久人妻福利社区极品人妻图片| www日本在线高清视频| 亚洲七黄色美女视频| 欧美性长视频在线观看| 久久亚洲真实| 日韩中文字幕欧美一区二区| 麻豆国产av国片精品| 黄色 视频免费看| 一区在线观看完整版| 国产精品综合久久久久久久免费 | 精品国产亚洲在线| 麻豆国产av国片精品| 欧美一区二区精品小视频在线| 欧美日韩福利视频一区二区| 亚洲少妇的诱惑av| 亚洲欧美精品综合久久99| av在线播放免费不卡| 亚洲国产精品久久男人天堂| 欧美黄色片欧美黄色片| 亚洲在线自拍视频| 亚洲av日韩精品久久久久久密| 色播在线永久视频| 叶爱在线成人免费视频播放| 欧美激情久久久久久爽电影 | 欧美精品啪啪一区二区三区| 淫秽高清视频在线观看| 亚洲专区字幕在线| 叶爱在线成人免费视频播放| 性欧美人与动物交配| 老汉色∧v一级毛片| a在线观看视频网站| 国产精品1区2区在线观看.| 国产伦人伦偷精品视频| 国产成人精品久久二区二区免费| 久久久水蜜桃国产精品网| 久久久久精品国产欧美久久久| 国产又色又爽无遮挡免费看| 夜夜夜夜夜久久久久| 身体一侧抽搐| 亚洲欧美精品综合久久99| 国产精品久久久久久亚洲av鲁大| 午夜免费激情av| 国产精品乱码一区二三区的特点 | 国产99久久九九免费精品| av视频免费观看在线观看| 在线av久久热| 每晚都被弄得嗷嗷叫到高潮| 国产成人影院久久av| 欧美 亚洲 国产 日韩一| 久久久久久久午夜电影| 欧美大码av| 中文字幕人成人乱码亚洲影| 性少妇av在线| 少妇粗大呻吟视频| 免费看十八禁软件| 国产精品乱码一区二三区的特点 | 亚洲国产日韩欧美精品在线观看 | 亚洲国产精品999在线| 国产成人一区二区三区免费视频网站| 精品国产超薄肉色丝袜足j| 后天国语完整版免费观看| 动漫黄色视频在线观看| 久久久精品国产亚洲av高清涩受| 巨乳人妻的诱惑在线观看| 成人18禁在线播放| 亚洲av成人av| 亚洲一区二区三区不卡视频| 一个人观看的视频www高清免费观看 | 免费在线观看完整版高清| 精品国内亚洲2022精品成人| 久久中文字幕一级| 欧美日韩精品网址| 男人操女人黄网站| 久久人人精品亚洲av| 欧美成人免费av一区二区三区| 色综合亚洲欧美另类图片| 人妻丰满熟妇av一区二区三区| 91av网站免费观看| 变态另类成人亚洲欧美熟女 | 欧美av亚洲av综合av国产av| 亚洲精品国产色婷婷电影| 久久精品影院6| 亚洲av成人一区二区三| 色播亚洲综合网| 亚洲熟女毛片儿| 可以在线观看的亚洲视频| 国产野战对白在线观看| 亚洲国产中文字幕在线视频| АⅤ资源中文在线天堂| av电影中文网址| av超薄肉色丝袜交足视频| 999精品在线视频| 免费不卡黄色视频| 久久婷婷人人爽人人干人人爱 | 国产伦人伦偷精品视频| 女同久久另类99精品国产91| 18禁黄网站禁片午夜丰满| 亚洲国产精品sss在线观看| www.熟女人妻精品国产| 国产三级黄色录像| 亚洲国产精品合色在线| tocl精华| 午夜影院日韩av| 亚洲第一电影网av| 国产亚洲欧美在线一区二区| 久久久久国产精品人妻aⅴ院| a在线观看视频网站| 自线自在国产av| 午夜福利成人在线免费观看| 免费高清视频大片| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产欧美一区二区综合| 欧美一区二区精品小视频在线| www国产在线视频色| 久久精品国产清高在天天线| 亚洲av片天天在线观看| 日本 欧美在线| 久久久久久人人人人人| 少妇的丰满在线观看| 宅男免费午夜| 91在线观看av| 久久久久国内视频| 亚洲av成人不卡在线观看播放网| 午夜成年电影在线免费观看| 美女高潮到喷水免费观看| 少妇熟女aⅴ在线视频| 午夜福利视频1000在线观看 | 成人18禁在线播放| 久久久久久久午夜电影| 午夜激情av网站| 精品熟女少妇八av免费久了| 一a级毛片在线观看| 日韩大尺度精品在线看网址 | 激情视频va一区二区三区| 性少妇av在线| 久久精品国产清高在天天线| 91在线观看av| a在线观看视频网站| 亚洲精品在线观看二区| www.熟女人妻精品国产| 一个人免费在线观看的高清视频| 久热爱精品视频在线9| 国产欧美日韩综合在线一区二区| 日韩精品免费视频一区二区三区| 可以在线观看毛片的网站| 两个人视频免费观看高清| 999精品在线视频| 中文字幕色久视频| 88av欧美| 不卡av一区二区三区| 国产亚洲av嫩草精品影院| 久久久国产欧美日韩av| 最新美女视频免费是黄的| 精品久久久久久久毛片微露脸| 国产精品一区二区精品视频观看| 国产国语露脸激情在线看| 欧美黑人精品巨大| 波多野结衣巨乳人妻| 日本 欧美在线| 嫩草影视91久久| 丰满人妻熟妇乱又伦精品不卡| 亚洲av日韩精品久久久久久密| 国产精品99久久99久久久不卡| 日韩有码中文字幕| 欧美黑人精品巨大| 中文字幕高清在线视频| 免费看十八禁软件| 涩涩av久久男人的天堂| 精品第一国产精品| 男女之事视频高清在线观看| 免费在线观看影片大全网站| 成熟少妇高潮喷水视频| 久久香蕉国产精品| 午夜激情av网站| 久久久精品国产亚洲av高清涩受| 少妇被粗大的猛进出69影院| 在线国产一区二区在线| 十八禁人妻一区二区| 国产亚洲精品第一综合不卡| 淫妇啪啪啪对白视频| 日日干狠狠操夜夜爽| 99久久久亚洲精品蜜臀av| 熟女少妇亚洲综合色aaa.| 色综合站精品国产| 国产视频一区二区在线看| 国产区一区二久久| 亚洲熟妇中文字幕五十中出| 激情在线观看视频在线高清| 热99re8久久精品国产| 国产免费av片在线观看野外av| 精品卡一卡二卡四卡免费| 亚洲人成电影观看| 日韩成人在线观看一区二区三区| 中文字幕久久专区| av超薄肉色丝袜交足视频| 精品人妻1区二区| 黄片播放在线免费| 99国产综合亚洲精品| 欧美不卡视频在线免费观看 | 亚洲成人国产一区在线观看| 久久青草综合色| 久久久久久大精品| 中文字幕精品免费在线观看视频| 这个男人来自地球电影免费观看| 精品少妇一区二区三区视频日本电影| 人成视频在线观看免费观看| 女性被躁到高潮视频| www.www免费av| 日韩中文字幕欧美一区二区| a级毛片在线看网站| 亚洲国产精品合色在线| 久久精品成人免费网站| 激情视频va一区二区三区| 在线观看日韩欧美| 精品人妻1区二区| 熟女少妇亚洲综合色aaa.| 美女高潮喷水抽搐中文字幕| 亚洲中文日韩欧美视频| 免费观看人在逋| 侵犯人妻中文字幕一二三四区| 69精品国产乱码久久久| 免费看a级黄色片| 欧美中文综合在线视频| 十八禁网站免费在线| 在线国产一区二区在线| 国产成人精品无人区| 精品一区二区三区视频在线观看免费| 99久久国产精品久久久| 国产真人三级小视频在线观看| 亚洲情色 制服丝袜| 久久久国产精品麻豆| 岛国在线观看网站| 国产精华一区二区三区| 91字幕亚洲| 色综合欧美亚洲国产小说| 免费高清视频大片| 变态另类成人亚洲欧美熟女 | 久久天堂一区二区三区四区| 亚洲va日本ⅴa欧美va伊人久久| 午夜老司机福利片|