• <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)定窗
    寂寞人妻少妇视频99o| 欧美 亚洲 国产 日韩一| 91成人精品电影| 97精品久久久久久久久久精品| 在线 av 中文字幕| 美女cb高潮喷水在线观看| 国语对白做爰xxxⅹ性视频网站| 日本av手机在线免费观看| 午夜久久久在线观看| 亚洲av.av天堂| 国产一区二区在线观看av| 五月伊人婷婷丁香| 伊人亚洲综合成人网| 少妇被粗大的猛进出69影院 | 97在线人人人人妻| 极品少妇高潮喷水抽搐| 伊人久久国产一区二区| 久久久国产精品麻豆| 成人亚洲精品一区在线观看| 成人亚洲欧美一区二区av| 99精国产麻豆久久婷婷| 看十八女毛片水多多多| 久久精品熟女亚洲av麻豆精品| 免费黄色在线免费观看| 2022亚洲国产成人精品| 天天躁夜夜躁狠狠久久av| 日韩人妻高清精品专区| 国产精品无大码| 看非洲黑人一级黄片| 国产精品免费大片| 午夜老司机福利剧场| 2021少妇久久久久久久久久久| 欧美精品亚洲一区二区| 久热久热在线精品观看| 99热这里只有是精品50| 看免费成人av毛片| a级毛片在线看网站| 日韩制服骚丝袜av| av一本久久久久| 国产av国产精品国产| 在线观看美女被高潮喷水网站| 噜噜噜噜噜久久久久久91| 日韩av免费高清视频| 日韩强制内射视频| 9色porny在线观看| 久久久久久久久久久丰满| 五月天丁香电影| av卡一久久| 免费观看性生交大片5| 一区二区三区免费毛片| 日韩强制内射视频| 国产日韩一区二区三区精品不卡 | 一级a做视频免费观看| 99热网站在线观看| 久久国产精品男人的天堂亚洲 | 久久久精品免费免费高清| 日韩,欧美,国产一区二区三区| 国产亚洲5aaaaa淫片| 黑人猛操日本美女一级片| 男人和女人高潮做爰伦理| 观看免费一级毛片| 午夜福利网站1000一区二区三区| 日本猛色少妇xxxxx猛交久久| 99久久精品国产国产毛片| 国产乱人偷精品视频| h日本视频在线播放| 如日韩欧美国产精品一区二区三区 | 午夜福利,免费看| 91久久精品电影网| 亚洲av电影在线观看一区二区三区| 熟女人妻精品中文字幕| 大又大粗又爽又黄少妇毛片口| 亚洲欧美一区二区三区黑人 | 日韩精品免费视频一区二区三区 | 久久国产精品男人的天堂亚洲 | 女性生殖器流出的白浆| 少妇丰满av| 黄色毛片三级朝国网站 | 日日撸夜夜添| 精品卡一卡二卡四卡免费| 高清毛片免费看| 国产在线视频一区二区| 日本-黄色视频高清免费观看| 国产成人精品福利久久| 国产免费一区二区三区四区乱码| av一本久久久久| 国产欧美日韩综合在线一区二区 | 国产精品成人在线| 国产精品福利在线免费观看| 十分钟在线观看高清视频www | 国产又色又爽无遮挡免| 一区二区三区精品91| 久久午夜综合久久蜜桃| 国产一级毛片在线| 黑丝袜美女国产一区| 色婷婷av一区二区三区视频| 女性被躁到高潮视频| 久久精品熟女亚洲av麻豆精品| 永久免费av网站大全| 又粗又硬又长又爽又黄的视频| 国产精品熟女久久久久浪| av在线老鸭窝| 久久久亚洲精品成人影院| 少妇的逼好多水| 国产 精品1| 国产伦精品一区二区三区四那| 久久久久久久久久久久大奶| 久久午夜综合久久蜜桃| 一级av片app| 国产精品不卡视频一区二区| 久久久国产精品麻豆| 大又大粗又爽又黄少妇毛片口| 日日摸夜夜添夜夜爱| 亚洲四区av| 欧美xxxx性猛交bbbb| 国产精品久久久久久精品古装| 人妻 亚洲 视频| 岛国毛片在线播放| 精品人妻偷拍中文字幕| av又黄又爽大尺度在线免费看| 男女边吃奶边做爰视频| 国产亚洲最大av| 国模一区二区三区四区视频| 狠狠精品人妻久久久久久综合| 伊人久久精品亚洲午夜| 在线观看免费日韩欧美大片 | 国产精品偷伦视频观看了| 菩萨蛮人人尽说江南好唐韦庄| 菩萨蛮人人尽说江南好唐韦庄| 国产精品免费大片| 麻豆成人午夜福利视频| 国产免费又黄又爽又色| 大香蕉久久网| 免费人妻精品一区二区三区视频| 亚洲av成人精品一二三区| 国产视频首页在线观看| 男人爽女人下面视频在线观看| tube8黄色片| 看十八女毛片水多多多| 欧美精品一区二区大全| 成人国产麻豆网| 一本—道久久a久久精品蜜桃钙片| 有码 亚洲区| 搡老乐熟女国产| 各种免费的搞黄视频| 久久久久视频综合| 久久97久久精品| 国产精品一区www在线观看| 国产av一区二区精品久久| 免费高清在线观看视频在线观看| 嫩草影院入口| 精品人妻一区二区三区麻豆| 国产在线视频一区二区| 内射极品少妇av片p| 性色avwww在线观看| 中文字幕免费在线视频6| 五月开心婷婷网| 街头女战士在线观看网站| 亚洲一级一片aⅴ在线观看| 亚洲高清免费不卡视频| 嘟嘟电影网在线观看| 80岁老熟妇乱子伦牲交| 日本黄大片高清| 欧美bdsm另类| 亚洲av男天堂| 亚洲av二区三区四区| 伊人亚洲综合成人网| av免费在线看不卡| 偷拍熟女少妇极品色| 中文乱码字字幕精品一区二区三区| 久久精品久久精品一区二区三区| 女性被躁到高潮视频| 97精品久久久久久久久久精品| 国产一区二区在线观看日韩| 亚洲欧美精品专区久久| 一本色道久久久久久精品综合| 丰满乱子伦码专区| 一级毛片电影观看| 18禁在线无遮挡免费观看视频| 国产亚洲精品久久久com| 丝袜在线中文字幕| 久久精品久久久久久久性| 爱豆传媒免费全集在线观看| 六月丁香七月| av黄色大香蕉| 中文字幕人妻丝袜制服| 亚洲综合色惰| 欧美日韩国产mv在线观看视频| 九九爱精品视频在线观看| 亚洲中文av在线| 99视频精品全部免费 在线| a级毛片在线看网站| 在线观看av片永久免费下载| 亚洲美女视频黄频| 久久久国产一区二区| 18禁在线无遮挡免费观看视频| 日日啪夜夜爽| 深夜a级毛片| 日韩欧美精品免费久久| 国产毛片在线视频| 各种免费的搞黄视频| 日韩一区二区三区影片| 插逼视频在线观看| 偷拍熟女少妇极品色| 国产在线一区二区三区精| 人妻夜夜爽99麻豆av| 色视频在线一区二区三区| 十八禁高潮呻吟视频 | 一二三四中文在线观看免费高清| 亚洲丝袜综合中文字幕| 成人国产麻豆网| 简卡轻食公司| 国产白丝娇喘喷水9色精品| 777米奇影视久久| 久久99热6这里只有精品| 欧美最新免费一区二区三区| 国产综合精华液| 乱码一卡2卡4卡精品| 91午夜精品亚洲一区二区三区| 一二三四中文在线观看免费高清| 两个人免费观看高清视频 | 国产亚洲av片在线观看秒播厂| freevideosex欧美| 一区二区三区免费毛片| 亚洲人成网站在线观看播放| 成人毛片60女人毛片免费| 69精品国产乱码久久久| 久久久久人妻精品一区果冻| 亚洲国产精品成人久久小说| 狠狠精品人妻久久久久久综合| 国产无遮挡羞羞视频在线观看| 亚洲欧美成人精品一区二区| 男女免费视频国产| 久久久久久伊人网av| 久久国产亚洲av麻豆专区| freevideosex欧美| 亚洲av中文av极速乱| 色婷婷av一区二区三区视频| 久久 成人 亚洲| 热re99久久精品国产66热6| 又爽又黄a免费视频| 99久久精品一区二区三区| 国产熟女欧美一区二区| 一区在线观看完整版| 高清av免费在线| 香蕉精品网在线| 亚洲欧美精品自产自拍| 久久韩国三级中文字幕| 蜜桃在线观看..| 我的女老师完整版在线观看| 国产亚洲一区二区精品| 纵有疾风起免费观看全集完整版| 卡戴珊不雅视频在线播放| 观看免费一级毛片| 五月伊人婷婷丁香| 国产视频内射| 永久免费av网站大全| 91精品国产国语对白视频| 午夜视频国产福利| 久久久久视频综合| 久久精品国产亚洲av天美| 亚洲精品视频女| 成人无遮挡网站| 中文资源天堂在线| 97在线视频观看| 精品亚洲乱码少妇综合久久| 日本色播在线视频| 人妻系列 视频| 欧美xxⅹ黑人| 一区二区三区乱码不卡18| 欧美人与善性xxx| 一级片'在线观看视频| 国产精品蜜桃在线观看| 久久99精品国语久久久| 日韩 亚洲 欧美在线| h日本视频在线播放| 搡老乐熟女国产| tube8黄色片| 久久精品国产自在天天线| 国产成人aa在线观看| 精品国产一区二区三区久久久樱花| 在线观看国产h片| 99热这里只有是精品50| 精品久久久久久电影网| 国产亚洲一区二区精品| 亚洲av中文av极速乱| 亚洲四区av| 性高湖久久久久久久久免费观看| 一级毛片 在线播放| 18禁在线播放成人免费| 久久精品久久久久久噜噜老黄| 在线观看免费日韩欧美大片 | 91午夜精品亚洲一区二区三区| 日韩av免费高清视频| 亚洲精品日本国产第一区| 国产精品蜜桃在线观看| 熟女av电影| 国产av精品麻豆| 久久久久久久久大av| 亚洲国产最新在线播放| 久久久久网色| 日韩一本色道免费dvd| 一区在线观看完整版| 在线观看三级黄色| 亚洲美女黄色视频免费看| 搡老乐熟女国产| 亚洲精品成人av观看孕妇| 亚洲无线观看免费| 蜜臀久久99精品久久宅男| 国产高清国产精品国产三级| 人妻系列 视频| 精品少妇黑人巨大在线播放| 99热国产这里只有精品6| videossex国产| 精品人妻一区二区三区麻豆| 国产成人freesex在线| av网站免费在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品不卡视频一区二区| av免费观看日本| 能在线免费看毛片的网站| 在线观看一区二区三区激情| 精品国产一区二区三区久久久樱花| 欧美激情国产日韩精品一区| 最黄视频免费看| 男男h啪啪无遮挡| 亚洲精品一区蜜桃| 日韩一区二区视频免费看| 成年美女黄网站色视频大全免费 | 伊人亚洲综合成人网| 高清不卡的av网站| 欧美激情国产日韩精品一区| 视频中文字幕在线观看| 秋霞伦理黄片| 日韩av免费高清视频| 欧美精品高潮呻吟av久久| 国产一区二区在线观看av| 纯流量卡能插随身wifi吗| 免费人妻精品一区二区三区视频| 一级黄片播放器| 美女大奶头黄色视频| 一本一本综合久久| 亚洲欧美一区二区三区黑人 | 中文天堂在线官网| 在现免费观看毛片| 人妻夜夜爽99麻豆av| 中文字幕精品免费在线观看视频 | 一个人看视频在线观看www免费| 搡老乐熟女国产| 国产av一区二区精品久久| 国产成人精品婷婷| 亚洲三级黄色毛片| 有码 亚洲区| 成人综合一区亚洲| 少妇的逼好多水| 国产免费又黄又爽又色| 国产亚洲欧美精品永久| 中文字幕久久专区| 午夜久久久在线观看| 国产色婷婷99| 九色成人免费人妻av| 在线观看人妻少妇| 久久久久视频综合| 精品一区二区三卡| 有码 亚洲区| 精品视频人人做人人爽| 亚洲内射少妇av| 中文字幕免费在线视频6| 永久免费av网站大全| 狂野欧美激情性xxxx在线观看| 91精品国产九色| 大陆偷拍与自拍| 国国产精品蜜臀av免费| 国产无遮挡羞羞视频在线观看| 女人久久www免费人成看片| 最近最新中文字幕免费大全7| 男人爽女人下面视频在线观看| 久久免费观看电影| 一区二区三区四区激情视频| 亚洲精品一区蜜桃| 91午夜精品亚洲一区二区三区| 美女国产视频在线观看| 欧美日韩国产mv在线观看视频| 国产91av在线免费观看| 插逼视频在线观看| 国产国拍精品亚洲av在线观看| 十八禁高潮呻吟视频 | 十八禁高潮呻吟视频 | 熟女人妻精品中文字幕| 黄色配什么色好看| 亚洲综合色惰| 国产中年淑女户外野战色| 成人亚洲精品一区在线观看| 九色成人免费人妻av| 纯流量卡能插随身wifi吗| 久久人妻熟女aⅴ| 大片电影免费在线观看免费| 国产精品嫩草影院av在线观看| 国产成人午夜福利电影在线观看| 久久国产乱子免费精品| 亚洲精品国产av成人精品| 久久国产精品大桥未久av | 伊人亚洲综合成人网| 黄色视频在线播放观看不卡| 日日啪夜夜撸| 亚洲精品456在线播放app| 五月伊人婷婷丁香| 久久精品国产a三级三级三级| 国产一区二区在线观看av| av黄色大香蕉| 黑人高潮一二区| 又黄又爽又刺激的免费视频.| 爱豆传媒免费全集在线观看| 在线观看免费视频网站a站| 亚洲内射少妇av| 午夜视频国产福利| 国产综合精华液| 日本黄色日本黄色录像| 美女内射精品一级片tv| 免费观看av网站的网址| 丝袜脚勾引网站| 久久免费观看电影| a级毛色黄片| 91aial.com中文字幕在线观看| 亚洲欧美一区二区三区国产| 欧美日韩一区二区视频在线观看视频在线| 青春草视频在线免费观看| av网站免费在线观看视频| 成人影院久久| 亚洲美女搞黄在线观看| 久久免费观看电影| 亚洲高清免费不卡视频| 久久久久久久亚洲中文字幕| 天堂俺去俺来也www色官网| 国产毛片在线视频| 国产亚洲91精品色在线| 在线观看www视频免费| 亚洲精品aⅴ在线观看| 精品卡一卡二卡四卡免费| 免费看不卡的av| 成人国产av品久久久| 免费观看a级毛片全部| 中文字幕av电影在线播放| 黄色欧美视频在线观看| 26uuu在线亚洲综合色| 蜜桃久久精品国产亚洲av| 亚洲av.av天堂| 欧美一级a爱片免费观看看| 午夜久久久在线观看| 97超碰精品成人国产| 狂野欧美激情性xxxx在线观看| 亚洲久久久国产精品| 久久6这里有精品| 三级经典国产精品| 色婷婷久久久亚洲欧美| 国产精品国产三级专区第一集| 热re99久久精品国产66热6| 久久精品久久久久久噜噜老黄| 搡女人真爽免费视频火全软件| 天堂8中文在线网| 三上悠亚av全集在线观看 | 精品久久久精品久久久| 精品少妇黑人巨大在线播放| 大片电影免费在线观看免费| 黑人巨大精品欧美一区二区蜜桃 | 日本91视频免费播放| 欧美精品一区二区免费开放| 国产美女午夜福利| 国产av码专区亚洲av| 久久6这里有精品| 亚洲欧美成人精品一区二区| 国产亚洲91精品色在线| 91成人精品电影| 国产极品粉嫩免费观看在线 | 久久影院123| 欧美97在线视频| 精品久久久久久电影网| 精品一区二区三区视频在线| 久热这里只有精品99| 免费不卡的大黄色大毛片视频在线观看| 99热全是精品| 天天操日日干夜夜撸| 性色avwww在线观看| 高清午夜精品一区二区三区| 少妇人妻精品综合一区二区| 免费观看a级毛片全部| 精品亚洲成a人片在线观看| 一级a做视频免费观看| 亚洲久久久国产精品| av国产久精品久网站免费入址| 国产成人精品无人区| 欧美丝袜亚洲另类| 99热网站在线观看| 赤兔流量卡办理| 久久97久久精品| 欧美日韩一区二区视频在线观看视频在线| 少妇人妻精品综合一区二区| 最黄视频免费看| 国产一区有黄有色的免费视频| 日韩强制内射视频| 日韩熟女老妇一区二区性免费视频| 大香蕉97超碰在线| 国产精品成人在线| 乱码一卡2卡4卡精品| a级毛片免费高清观看在线播放| 国产欧美亚洲国产| 老司机影院成人| 欧美精品一区二区免费开放| 中文在线观看免费www的网站| 免费av中文字幕在线| 亚洲国产日韩一区二区| 99久久精品国产国产毛片| 最新的欧美精品一区二区| 国产色爽女视频免费观看| 少妇人妻久久综合中文| 激情五月婷婷亚洲| 性色avwww在线观看| 国产成人精品久久久久久| 久久精品久久精品一区二区三区| 麻豆成人av视频| 日本av手机在线免费观看| 美女xxoo啪啪120秒动态图| 人妻 亚洲 视频| 国产视频首页在线观看| 18禁动态无遮挡网站| 精品一区二区免费观看| 新久久久久国产一级毛片| 又粗又硬又长又爽又黄的视频| 午夜免费男女啪啪视频观看| 国产 一区精品| 高清毛片免费看| 亚洲av男天堂| 午夜影院在线不卡| 一级片'在线观看视频| 日韩一区二区视频免费看| 一级毛片电影观看| 亚洲精华国产精华液的使用体验| 国产精品国产三级国产av玫瑰| 久久久久久久久久久丰满| 人人妻人人澡人人看| 在线亚洲精品国产二区图片欧美 | 国产免费一级a男人的天堂| 色婷婷久久久亚洲欧美| 永久免费av网站大全| 丝袜脚勾引网站| 中文欧美无线码| 欧美日韩亚洲高清精品| 亚洲精品日韩在线中文字幕| 曰老女人黄片| 国产精品成人在线| 午夜影院在线不卡| 国产高清有码在线观看视频| 日韩欧美 国产精品| 中文字幕免费在线视频6| 一级a做视频免费观看| 日本欧美视频一区| 有码 亚洲区| 欧美高清成人免费视频www| 成人影院久久| 色5月婷婷丁香| 黄色怎么调成土黄色| 欧美激情极品国产一区二区三区 | 亚洲成人av在线免费| av国产久精品久网站免费入址| 夫妻性生交免费视频一级片| 亚洲高清免费不卡视频| av天堂中文字幕网| 能在线免费看毛片的网站| 国产有黄有色有爽视频| 一本大道久久a久久精品| 欧美日韩一区二区视频在线观看视频在线| 日韩欧美精品免费久久| 一级毛片黄色毛片免费观看视频| 99久久精品国产国产毛片| 天天操日日干夜夜撸| 免费在线观看成人毛片| 一区在线观看完整版| 久久久久久久久久久丰满| 亚洲一级一片aⅴ在线观看| 在线播放无遮挡| 国产精品久久久久成人av| av免费观看日本| 男女无遮挡免费网站观看| 男男h啪啪无遮挡| 日产精品乱码卡一卡2卡三| 少妇精品久久久久久久| 国产av国产精品国产| 日韩av不卡免费在线播放| 国产在线一区二区三区精| 久久ye,这里只有精品| 女人久久www免费人成看片| 在线 av 中文字幕| 亚洲欧洲国产日韩| 日韩一区二区视频免费看| 嫩草影院入口| 街头女战士在线观看网站| 男人添女人高潮全过程视频| 热re99久久国产66热| 国产黄频视频在线观看| 99热这里只有是精品50| 国产免费视频播放在线视频| 最新的欧美精品一区二区| 丁香六月天网| 人人妻人人添人人爽欧美一区卜| 男的添女的下面高潮视频| 十八禁高潮呻吟视频 | 五月开心婷婷网| 国产乱来视频区| 亚洲一级一片aⅴ在线观看| 黄色配什么色好看| 少妇人妻精品综合一区二区| 七月丁香在线播放| a级一级毛片免费在线观看| 91久久精品国产一区二区成人| 在线观看av片永久免费下载|