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

    求解熱傳導(dǎo)系數(shù)反問(wèn)題的量子行為粒子群算法

    2014-08-05 02:41:28朱龍超
    關(guān)鍵詞:熱傳導(dǎo)正則量子

    田 娜,朱龍超

    1.江南大學(xué) 教育技術(shù)系,江蘇 無(wú)錫 214122

    2.中國(guó)船舶科學(xué)研究中心 信息技術(shù)室,江蘇 無(wú)錫 214082

    求解熱傳導(dǎo)系數(shù)反問(wèn)題的量子行為粒子群算法

    田 娜1,朱龍超2

    1.江南大學(xué) 教育技術(shù)系,江蘇 無(wú)錫 214122

    2.中國(guó)船舶科學(xué)研究中心 信息技術(shù)室,江蘇 無(wú)錫 214082

    1 引言

    熱力學(xué)中的熱傳導(dǎo)系數(shù),通常是通過(guò)對(duì)流或流體和固體之間的相變,來(lái)計(jì)算傳熱過(guò)程。對(duì)板表面的熱傳導(dǎo)系數(shù)精確地了解,在許多工程應(yīng)用中都有很重要的作用,包括連鑄坯和電子芯片的冷卻[1]。

    在過(guò)去的幾十年里,許多工作已經(jīng)致力于對(duì)隨時(shí)間變化(時(shí)變)的熱傳導(dǎo)系數(shù)的研究[1]。Su和Hewitt[2]使用Alifanov迭代正則化方法來(lái)估計(jì)強(qiáng)對(duì)流流體在加熱管外表面上沸騰時(shí)的熱傳導(dǎo)系數(shù)。在文獻(xiàn)[3]中,三種版本的共軛梯度法用來(lái)估計(jì)平板表面的熱傳導(dǎo)系數(shù),平板通過(guò)與降溫液體的對(duì)流而散熱。Chen和Wu[4]應(yīng)用了一種包含拉普拉斯變換,有限差分和最小二乘法的混合策略,連同時(shí)間順序的概念,三次樣條和溫度測(cè)量值一起,用于預(yù)測(cè)邊界表面的熱傳導(dǎo)系數(shù)分布。Slodicka[5]采用邊界元法和Tikhonov正則化來(lái)建構(gòu)隨時(shí)間變化的熱傳導(dǎo)系數(shù)。Chantasiriwan[6]采用順序函數(shù)估計(jì)法連同線(xiàn)性基函數(shù)和線(xiàn)性變化的未來(lái)邊界熱通量或溫度假設(shè)來(lái)估計(jì)隨時(shí)間變化的Biot數(shù)字。在本文中,采用QPSO算法和Tikhonov正則化方法[7-9]用來(lái)估計(jì)時(shí)變熱傳導(dǎo)系數(shù)。QPSO算法是由Sun[10-12]在PSO的基礎(chǔ)上提出的。PSO算法最初是由Kennedy和Eberhart提出的,是模擬魚(yú)群和鳥(niǎo)群的社會(huì)行為的過(guò)程。與GA比較,PSO和QPSO的優(yōu)點(diǎn)在于:有更少的參數(shù)需要控制,運(yùn)行起來(lái)更簡(jiǎn)單,并且,算法只需要簡(jiǎn)單的算術(shù)運(yùn)算符,而不需要像GA中選擇,交叉,變異之類(lèi)的操作。最后,結(jié)果與CGM估計(jì)得到的結(jié)果進(jìn)行了比較。

    2 問(wèn)題描述

    液體在平板的上面以一個(gè)恒定的溫度流動(dòng)[1],如圖1所示。

    圖1 電加熱平板

    如果平板突然被一個(gè)電加熱器加熱,板的溫度會(huì)上升。假設(shè)平板的材料沿y軸均勻,對(duì)流邊界條件在x=0和x=L是指定的。數(shù)學(xué)模型為:

    其中T(x,t)是在位置x和時(shí)間t的溫度分布,g(t)是熱源在x=xs的強(qiáng)度,T∞是周?chē)h(huán)境的溫度,T0是初始溫度分布。為了描述簡(jiǎn)單,物理屬性被設(shè)為K=ρC=L=1,這與采用無(wú)量綱的數(shù)據(jù)是一樣的。在這里,h(t)是未知的待確定的熱傳導(dǎo)系數(shù)。對(duì)方程(1)采用隱式有限差分法:

    在邊界x=L,對(duì)流邊界條件用二階離散得到:

    平均誤差herror,用來(lái)評(píng)估估計(jì)得到的熱傳導(dǎo)系數(shù),可以定義為:

    其中Nt是時(shí)間步數(shù),hj是估計(jì)得到的熱傳導(dǎo)系數(shù),h~j是準(zhǔn)確的熱傳導(dǎo)系數(shù)。

    最小二乘模型用來(lái)求解此類(lèi)反問(wèn)題:

    這種方法求解反問(wèn)題的過(guò)程是不適定性的,不可避免的測(cè)量和計(jì)算誤差常常導(dǎo)致不穩(wěn)定和不精確的結(jié)果。所以,采用正則化技術(shù)來(lái)得到穩(wěn)定解[7],它可以將目標(biāo)函數(shù)式(8)改寫(xiě)成如下適定形式:

    其中λ是正則參數(shù),公式(9)右邊第一項(xiàng)是差異項(xiàng),第二項(xiàng)是正則項(xiàng)。正則項(xiàng)通??梢员硎緸槿缦滦问剑?/p>

    其中n是正則階數(shù)。常用的有零階和一階正則項(xiàng)。當(dāng)n=0時(shí),為零階正則項(xiàng):

    3 量子行為粒子群優(yōu)化算法

    3.1 粒子群優(yōu)化算法

    粒子群優(yōu)化算法(PSO)是基于群體智能的優(yōu)化技術(shù),最初由Kennedy和Eberhart在1995提出的[13]。該算法的概念來(lái)源于鳥(niǎo)群或魚(yú)群的社會(huì)行為。該系統(tǒng)有一個(gè)粒子群體,其中每個(gè)粒子代表一個(gè)優(yōu)化問(wèn)題可能的解。已經(jīng)證明PSO算法與遺傳算法(GA)有相當(dāng)?shù)男阅躘14]。

    在有M個(gè)粒子,D維空間的原始PSO算法中,第i個(gè)粒子在第k次迭代步的位置向量和速度向量表示為:Xi(k)=(Xi1(k),Xi2(k),…,XiD(k)),Vi(k)=(Vi1(k),Vi2(k),…,ViD(k))。粒子根據(jù)以下迭代公式更新速度和位置:

    其中i=1,2,…,M ,j=1,2,…,D,c1和c2是加速系數(shù)。r1和r2是均勻分布在(0,1)中的隨機(jī)數(shù)。向量Pi= (Pi1,Pi2,…,PiD)是第i個(gè)粒子的歷史最佳位置。向量Pg=(Pg1,Pg2,…,PgD)是整個(gè)粒子群的歷史最佳位置。

    3.2 量子行為粒子群優(yōu)化算法

    原始PSO算法的缺點(diǎn)是不能確保收斂到全局最優(yōu)解[15],為了克服這個(gè)缺點(diǎn),Sun在2004年提出了量子行為粒子群優(yōu)化算法的概念,文獻(xiàn)[10-11]中的軌跡分析說(shuō)明,當(dāng)每個(gè)粒子收斂到 pi=(pi1,pi2,…,piD)時(shí),PSO算法實(shí)現(xiàn)收斂:

    其中φ∈(0,1)。

    在量子世界中,粒子的速度是沒(méi)有意義的。所以在QPSO中,位置是唯一描述粒子狀態(tài)的變量,更新公式如:

    其中mbest(k)稱(chēng)作平均最優(yōu)位置,定義為所有粒子歷史最優(yōu)位置的平均值:

    其中M是粒子的個(gè)數(shù)。公式(16)中的參數(shù)α是收縮系數(shù),通過(guò)調(diào)整可以控制收斂速度。與PSO不同的是,QPSO算法不需要算法向量,并且只有一個(gè)參數(shù)需要控制,使得算法更容易執(zhí)行。標(biāo)準(zhǔn)測(cè)試函數(shù)證明,QPSO算法的性能要優(yōu)于PSO算法[10-12]。

    在本文中,QPSO算法用來(lái)估計(jì)隨時(shí)間變化的熱傳導(dǎo)系數(shù)h(t),其中,每個(gè)粒子X(jué)i(k)作為h(t)的一個(gè)候選解:

    其中D=Nt。

    QPSO與Tikhonov正則化方法估計(jì)h(t)的過(guò)程可以描述為:

    選擇N個(gè)正則參數(shù)值{λ1,λ2,…,λN-1,λN}以及與之對(duì)應(yīng)的兩個(gè)向量Residual和Norm;對(duì)于每一個(gè)正則參數(shù) λj(j=1,2,…,N)

    Step1初始化:粒子位置

    X(0)={X1(0),X2(0),…,Xi(0),…,XM(0)}

    粒子歷史最優(yōu)位置

    P(0)={P1(0),P2(0),…,Pi(0),…,PM(0)}

    全局最優(yōu)位置Pg,縮放系數(shù)α=1.0,k=0。

    Step2 while(k<kmax)

    根據(jù)公式(17)計(jì)算mbest;根據(jù)公式(15)計(jì)算 p;根據(jù)公式(16)更新每個(gè)粒子的位置;根據(jù)公式(12)評(píng)估每個(gè)粒子在λj下的適應(yīng)值;更新P和Pg;更新α;

    Step3將向量Residual和Norm的點(diǎn)畫(huà)在坐標(biāo)系中,以獲得最優(yōu)正則系數(shù)λopt。Go to Step1。

    Step4輸出估計(jì)得到的最優(yōu)熱傳導(dǎo)系數(shù)hopt。

    4 數(shù)值測(cè)試

    在隨后的數(shù)值實(shí)驗(yàn)中,采用方波函數(shù)[2]:

    這是最難估計(jì)的一種情況。問(wèn)題的設(shè)置為:

    為了探討傳感器的數(shù)量和位置對(duì)結(jié)果的影響,表1列出了一系列不同的案例。圖2顯示了分別采用1,3,5個(gè)傳感器的結(jié)果。從中可以注意到,用3個(gè)和5個(gè)傳感器估計(jì)的結(jié)果并不優(yōu)于只用一個(gè)傳感器得到的結(jié)果。因此,這表明,一個(gè)傳感器足以獲得令人滿(mǎn)意的估計(jì)問(wèn)題。

    表1 傳感器的數(shù)量和位置對(duì)結(jié)果的影響

    圖2 QPSO用不同數(shù)量的傳感器估計(jì)的h(t)

    因此,在隨后的測(cè)試中只用一個(gè)傳感器來(lái)進(jìn)行溫度測(cè)量。圖3所示顯示了傳感器的位置對(duì)熱傳導(dǎo)系數(shù)的影響。表2列出了四種傳感器的位置和其相應(yīng)的平均誤差和目標(biāo)函數(shù)值。從中可以注意到,傳感器越接近對(duì)流邊界,結(jié)果越好,如圖3所示。

    圖3 QPSO采用不同位置的傳感器估計(jì)的h(t)

    表2 傳感器位置對(duì)結(jié)果的影響

    用模擬實(shí)驗(yàn)溫度來(lái)研究不同噪音水平對(duì)結(jié)果的影響。測(cè)試中采用了三個(gè)不同的噪聲水平。表3列出了QPSO和CGM算法的平均誤差和目標(biāo)函數(shù)值。圖4比較了二者的結(jié)果。從中注意到,由QPSO估計(jì)得到的出色結(jié)果,特別是對(duì)t=tf附近值的估計(jì)。圖5顯示了QPSO和CGM在估計(jì)熱傳導(dǎo)系數(shù)時(shí)的收斂曲線(xiàn)。注意到,CGM收斂速度很快,但是在一些迭代步之后,收斂就出現(xiàn)了停滯。雖然隨機(jī)初始化的QPSO算法的收斂速度相對(duì)緩慢,但是全局最優(yōu)解得到了保障。

    表3 在不同噪聲下QPSO和CGM的比較

    圖4 QPSO和CGM用不同水平的噪聲估計(jì)的h(t)

    圖5 CGM和QPSO估計(jì)的h(t)收斂曲線(xiàn)(ε=0)

    為了測(cè)試QPSO算法的魯棒性,再對(duì)兩個(gè)函數(shù)進(jìn)行測(cè)試,一個(gè)是持續(xù)很短時(shí)間的方波函數(shù):

    由QPSO和CGM估計(jì)的結(jié)果如圖6和圖7所示??梢钥闯?,由QPSO估計(jì)的時(shí)變熱傳導(dǎo)系數(shù)與準(zhǔn)確的熱傳導(dǎo)系數(shù)非常逼近。

    圖6 QPSO和CGM估計(jì)較小波函數(shù)的h(t)

    圖7 QPSO和CGM估計(jì)波形h(t)(ε=0)

    5 結(jié)論

    在本文中,QPSO算法用來(lái)對(duì)平板表面的時(shí)變熱傳導(dǎo)系數(shù)進(jìn)行估計(jì)。采用模擬測(cè)量溫度值得到的結(jié)果證明了QPSO算法的可行性和穩(wěn)定性。分析研究了傳感器的數(shù)量和位置對(duì)結(jié)果精確度的影響。通過(guò)對(duì)QPSO和CGM的比較,說(shuō)明了QPSO的優(yōu)越性。

    [1]Beck J V,Blackwell B,St-Clair C R.Inverse heat conduction:ill-posed problems[M].New York:Wiley-InterScience,1985.

    [2]Su J,Hewitt G.Inverse heat conduction problem of estimating time-varying heat transfer coefficient[J].Numerical Heat Transfer Part A:Applications,2004,45(8):777-789.

    [3]Colaco M J,Orlande H R B.Comparison of different versions of the conjugate gradient method of function estimation[J].Numerical Heat Transfer Part A:Applications,1999,36(2):229-249.

    [4]Chen H T,Wu X Y.Investigation of heat transfer coefficient in two-dimensional transient inverse heat conduction problems using the hybrid inverse scheme[J].International Journal for Numerical Methods in Engineering,2008,73(1):107-122.

    [5]Slodicka M,Lesnic D,Onyango T T M.Determination of a time-dependent heat transfer coefficient in a nonlinear inverse heat conduction problem[J].Inverse Problems in Science and Engineering,2010,18(1):65-81.

    [6]Chantasiriwan S.Inverse heat conduction problem of determining time-dependent heat transfer coefficient[J].International Journal of Heat Mass Transfer,2000,42(23):4275-4285.

    [7]Tikhonov A N,Arsenin V Y.Solution of ill-posed problems[M].Winston:Washington DC,1977.

    [8]Morozov V A.Methods for incorrectly posed problems[M]. [S.l.]:Springer,1984.

    [9]Hansen P C.The L-curve and its use in the numerical treatment of inverse problems[M]//Computational Inverse Problem in Electrocardiography,Advances in Computational Bioengineering.Holland:WIT Press,2001.

    [10]Sun J,F(xiàn)eng B,Xu W B.Particle swarm optimization with particles having quantum behaviour[C]//IEEE Congress on Evolutionary Computation,Portland,USA,2004.

    [11]Sun J,Xu W B,Liu J.Parameter selection of quantumbehaved particle swarm optimization[C]//Lecture Notes in Computer Science,2005,3612.

    [12]Sun J,Xu W B,F(xiàn)eng B.A global search strategy of quantum-behaved particle swarm optimization[C]//IEEE Conf on Cybernetics and Intelligent Systems,Singapore,2004.

    [13]Kennedy J,Eberhart R C.Particle swarm optimization[C]// IEEE Int Conf on Neural Networks,Perth,Australia,1995.

    [14]Eberhart R C,Shi Y.Comparison between genetic algorithm and particle swarm optimization[C]//Lecture Notes in Computer Science:Evolutionary Programming VII. Berlin,Heidelberg:Springer,1998,1447:611-616.

    [15]van den Bergh F.An analysis of particle swarm optimizers[D].South Africa:University of Pretoria,2001.

    [16]Clerc M,Kennedy J.The particle swarm:explosion,stability,and convergence in a multi-dimensional complex space[J].IEEE Transactions on Evolutionary Computing,2002,6(1):58-73.

    TIAN Na1,ZHU Longchao2

    1.Department of Educational Technology,Jiangnan University,Wuxi,Jiangsu 214122,China
    2.Department of Information Technology,China Ship Science Research Centre,Wuxi,Jiangsu 214082,China

    In this paper,the Quantum-behaved Particle Swarm Optimization(QPSO)with Tikhonov regularization is used to solve the inverse heat conduction problem of estimating the time dependent heat transfer coefficient of a flat plate.The prior information about the functional form of the unknown is unavailable.The estimation is based on transient temperature measurements taken by the sensors imbedded in the plate,which are used in the least square model,minimized by QPSO. The detail of choosing the best regularization parameter by L-curve method is presented.Numerical experiments are performed to test the proposed method.Effects of the location and number of sensors are also investigated.Comparison with conjugate gradient method is given as well.

    heat transfer coefficient;Quantum-behaved Particle Swarm Optimization(QPSO);Tikhonov regularization; conjugate gradient method;L-curve

    量子行為粒子群優(yōu)化算法(QPSO)和Tikhonov正則化方法用來(lái)求解熱傳導(dǎo)反問(wèn)題,近似估計(jì)平板隨時(shí)間變化的熱傳導(dǎo)系數(shù)。由于熱傳導(dǎo)系數(shù)的函數(shù)形式是未知的,所以問(wèn)題可以歸結(jié)為函數(shù)估計(jì)問(wèn)題。求解過(guò)程是基于最小二乘模型的,采用的是嵌在平板中的傳感器所測(cè)量得到的溫度,優(yōu)化過(guò)程由QPSO算法來(lái)求解。給出了由L曲線(xiàn)方法選擇正則參數(shù)的詳細(xì)過(guò)程。提出算法的有效性經(jīng)過(guò)了數(shù)值實(shí)驗(yàn)的驗(yàn)證。傳感器的位置和數(shù)量對(duì)結(jié)果的影響也做了研究。給出了與共軛梯度法的比較。

    熱傳導(dǎo)系數(shù);量子行為粒子群優(yōu)化算法;Tikhonov正則化;共軛梯度法;L曲線(xiàn)

    A

    TP391

    10.3778/j.issn.1002-8331.1305-0033

    TIAN Na,ZHU Longchao.Estimation of heat transfer coefficient using Quantum-behaved Particle Swarm Optimization.Computer Engineering and Applications,2014,50(24):266-270.

    江南大學(xué)自主科研項(xiàng)目基金(No.1245210382130120,No.1242050205142810)。

    田娜(1983—),女,博士,副教授,主要研究方向?yàn)橹悄苡?jì)算,模式識(shí)別,偏微分方程反問(wèn)題。

    2013-05-08

    2013-09-11

    1002-8331(2014)24-0266-05

    CNKI網(wǎng)絡(luò)優(yōu)先出版:2014-05-05,http∶//www.cnki.net/kcms/doi/10.3778/j.issn.1002-8331.1305-0033.html

    猜你喜歡
    熱傳導(dǎo)正則量子
    2022年諾貝爾物理學(xué)獎(jiǎng) 從量子糾纏到量子通信
    一類(lèi)三維逆時(shí)熱傳導(dǎo)問(wèn)題的數(shù)值求解
    決定未來(lái)的量子計(jì)算
    新量子通信線(xiàn)路保障網(wǎng)絡(luò)安全
    剩余有限Minimax可解群的4階正則自同構(gòu)
    類(lèi)似于VNL環(huán)的環(huán)
    熱傳導(dǎo)方程解的部分Schauder估計(jì)
    一類(lèi)非線(xiàn)性反向熱傳導(dǎo)問(wèn)題的Fourier正則化方法
    一種簡(jiǎn)便的超聲分散法制備碳量子點(diǎn)及表征
    有限秩的可解群的正則自同構(gòu)
    久久精品国产99精品国产亚洲性色| 亚洲专区中文字幕在线| 天堂影院成人在线观看| 欧美色视频一区免费| 国产精品精品国产色婷婷| 在线观看一区二区三区| 日本黄大片高清| 少妇高潮的动态图| 成人一区二区视频在线观看| 国产单亲对白刺激| 欧美黑人巨大hd| www.999成人在线观看| 一区福利在线观看| 老熟妇仑乱视频hdxx| 怎么达到女性高潮| 中文字幕精品亚洲无线码一区| 国产成人aa在线观看| 性欧美人与动物交配| 久久久久九九精品影院| 国产精品三级大全| 一区福利在线观看| www.999成人在线观看| 免费电影在线观看免费观看| 成人国产综合亚洲| 成人鲁丝片一二三区免费| www.999成人在线观看| 麻豆久久精品国产亚洲av| 韩国av一区二区三区四区| 久久99热6这里只有精品| 性色av乱码一区二区三区2| 日本黄大片高清| 如何舔出高潮| 99在线人妻在线中文字幕| 夜夜看夜夜爽夜夜摸| 狂野欧美白嫩少妇大欣赏| 男人舔女人下体高潮全视频| 91九色精品人成在线观看| 亚洲第一电影网av| 亚洲,欧美精品.| 国产精品国产高清国产av| 国产精品永久免费网站| 一个人观看的视频www高清免费观看| 色噜噜av男人的天堂激情| 中文资源天堂在线| 91麻豆av在线| 国产成人影院久久av| 精品国产三级普通话版| 日本黄色视频三级网站网址| 亚洲男人的天堂狠狠| 女人十人毛片免费观看3o分钟| 亚洲精品成人久久久久久| 又黄又爽又免费观看的视频| 麻豆av噜噜一区二区三区| 午夜亚洲福利在线播放| 亚洲美女视频黄频| 精品福利观看| 亚洲在线自拍视频| 精品熟女少妇八av免费久了| 亚洲精品粉嫩美女一区| 国产极品精品免费视频能看的| 亚洲在线自拍视频| 日韩 亚洲 欧美在线| 极品教师在线免费播放| 少妇丰满av| .国产精品久久| 国产久久久一区二区三区| 最近最新免费中文字幕在线| 国产黄a三级三级三级人| 91午夜精品亚洲一区二区三区 | 99热这里只有精品一区| 啦啦啦韩国在线观看视频| 国产伦精品一区二区三区四那| 在线a可以看的网站| 亚洲精品在线观看二区| 午夜a级毛片| 国产蜜桃级精品一区二区三区| 日本 av在线| 久久久久久国产a免费观看| 午夜免费激情av| 又爽又黄无遮挡网站| 国产 一区 欧美 日韩| 国产高清三级在线| 色综合欧美亚洲国产小说| 免费在线观看成人毛片| 欧美黑人巨大hd| 少妇裸体淫交视频免费看高清| 亚洲av.av天堂| 久久精品国产亚洲av香蕉五月| av国产免费在线观看| 日韩精品青青久久久久久| 91字幕亚洲| 国产主播在线观看一区二区| 久久久成人免费电影| 久久精品综合一区二区三区| 我的女老师完整版在线观看| 色噜噜av男人的天堂激情| 久久精品久久久久久噜噜老黄 | 国产一区二区在线av高清观看| 亚洲av.av天堂| 国产精品一及| 国产精品免费一区二区三区在线| 一区福利在线观看| 国产欧美日韩一区二区三| 久久精品国产清高在天天线| 欧美午夜高清在线| 国产欧美日韩一区二区三| 88av欧美| www.色视频.com| av中文乱码字幕在线| 日韩欧美 国产精品| 日本成人三级电影网站| 亚洲专区国产一区二区| 国产精品一区二区三区四区免费观看 | 在线免费观看不下载黄p国产 | 国产人妻一区二区三区在| 简卡轻食公司| 国产 一区 欧美 日韩| 天天一区二区日本电影三级| 国产一级毛片七仙女欲春2| 午夜免费成人在线视频| 国产精品,欧美在线| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩卡通动漫| 国产欧美日韩一区二区精品| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品久久男人天堂| 可以在线观看毛片的网站| 国产伦在线观看视频一区| 中文亚洲av片在线观看爽| 波多野结衣高清无吗| 婷婷六月久久综合丁香| 女生性感内裤真人,穿戴方法视频| 精品免费久久久久久久清纯| 国产成人影院久久av| 亚洲欧美日韩高清专用| 性色avwww在线观看| av女优亚洲男人天堂| 国内精品久久久久久久电影| 老司机午夜十八禁免费视频| 欧美日韩亚洲国产一区二区在线观看| 中文字幕人妻熟人妻熟丝袜美| 日日摸夜夜添夜夜添小说| 少妇熟女aⅴ在线视频| 国产精品久久久久久久电影| 精品久久久久久久久亚洲 | 国产高潮美女av| 成人三级黄色视频| 日韩欧美一区二区三区在线观看| 乱人视频在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美成狂野欧美在线观看| 首页视频小说图片口味搜索| 国产色爽女视频免费观看| 午夜免费男女啪啪视频观看 | 一个人免费在线观看电影| 在线观看午夜福利视频| 日韩 亚洲 欧美在线| 中文字幕精品亚洲无线码一区| 久久久国产成人免费| 国产精品美女特级片免费视频播放器| 精品人妻熟女av久视频| 在线看三级毛片| 美女 人体艺术 gogo| 国产成人欧美在线观看| 亚洲av电影不卡..在线观看| 内地一区二区视频在线| av在线老鸭窝| 免费在线观看日本一区| 脱女人内裤的视频| 真人做人爱边吃奶动态| 国产激情偷乱视频一区二区| 丁香六月欧美| 美女被艹到高潮喷水动态| 91在线观看av| 国内精品一区二区在线观看| 婷婷亚洲欧美| 网址你懂的国产日韩在线| 九九在线视频观看精品| 亚洲第一区二区三区不卡| 欧美一区二区精品小视频在线| 午夜福利在线在线| 中文字幕av成人在线电影| 成人国产一区最新在线观看| av在线观看视频网站免费| 少妇被粗大猛烈的视频| 一a级毛片在线观看| 丰满人妻一区二区三区视频av| 日日干狠狠操夜夜爽| 黄色视频,在线免费观看| 亚洲激情在线av| 99热只有精品国产| 高清日韩中文字幕在线| 美女高潮喷水抽搐中文字幕| 免费高清视频大片| 亚洲激情在线av| 国产私拍福利视频在线观看| 亚洲自偷自拍三级| 一a级毛片在线观看| 欧美午夜高清在线| 欧美黑人欧美精品刺激| 国产av不卡久久| 国产精品三级大全| 51午夜福利影视在线观看| 搡老妇女老女人老熟妇| 91久久精品国产一区二区成人| 亚洲精品久久国产高清桃花| 久99久视频精品免费| 亚洲人成伊人成综合网2020| 亚洲av电影不卡..在线观看| 久久国产精品人妻蜜桃| 精品一区二区三区av网在线观看| 日日夜夜操网爽| 亚洲国产日韩欧美精品在线观看| 在线播放无遮挡| 国产精品98久久久久久宅男小说| 国产精品一区二区性色av| 国产午夜福利久久久久久| 国产亚洲欧美98| 一区二区三区免费毛片| 久久久色成人| 国内精品久久久久久久电影| 亚洲国产精品久久男人天堂| 国产色婷婷99| 欧美激情国产日韩精品一区| 欧美最新免费一区二区三区 | 2021天堂中文幕一二区在线观| 国内揄拍国产精品人妻在线| 十八禁网站免费在线| 最近中文字幕高清免费大全6 | 在线观看一区二区三区| 欧美精品国产亚洲| 欧美性感艳星| 国产精品久久视频播放| 亚洲人与动物交配视频| 精品99又大又爽又粗少妇毛片 | 久久精品国产亚洲av天美| 精品久久久久久久久久免费视频| 91狼人影院| 成人特级黄色片久久久久久久| 国产乱人视频| 午夜a级毛片| 有码 亚洲区| 男女视频在线观看网站免费| 精品熟女少妇八av免费久了| 亚洲人成网站在线播| 欧美潮喷喷水| 国产精品一区二区三区四区免费观看 | 亚洲美女黄片视频| 国产精品免费一区二区三区在线| 午夜福利免费观看在线| 国模一区二区三区四区视频| 日本撒尿小便嘘嘘汇集6| 亚洲自偷自拍三级| 我要搜黄色片| 亚洲国产色片| 蜜桃久久精品国产亚洲av| 少妇的逼水好多| 九九在线视频观看精品| 深夜a级毛片| 亚洲自偷自拍三级| 欧美激情在线99| 别揉我奶头~嗯~啊~动态视频| 欧美另类亚洲清纯唯美| 天堂网av新在线| 免费无遮挡裸体视频| 午夜免费成人在线视频| 亚洲av免费在线观看| 99国产综合亚洲精品| 亚洲欧美清纯卡通| 国产单亲对白刺激| 看免费av毛片| 激情在线观看视频在线高清| 久久热精品热| 久久久成人免费电影| 国内毛片毛片毛片毛片毛片| 给我免费播放毛片高清在线观看| 免费看日本二区| 天堂av国产一区二区熟女人妻| 欧美日韩综合久久久久久 | 美女被艹到高潮喷水动态| 又紧又爽又黄一区二区| 97人妻精品一区二区三区麻豆| 久久人人爽人人爽人人片va | 日本三级黄在线观看| 国产探花极品一区二区| av在线观看视频网站免费| 国产一区二区在线av高清观看| 午夜福利在线观看吧| 国产亚洲欧美98| 国产精品国产高清国产av| 久久久久久久久大av| 男插女下体视频免费在线播放| 亚洲三级黄色毛片| 亚洲av二区三区四区| 久久久久久九九精品二区国产| 久久精品夜夜夜夜夜久久蜜豆| 综合色av麻豆| 一级作爱视频免费观看| 黄色丝袜av网址大全| 又黄又爽又免费观看的视频| 精品人妻1区二区| 麻豆国产97在线/欧美| 俺也久久电影网| 亚洲男人的天堂狠狠| 成人永久免费在线观看视频| 国产精品综合久久久久久久免费| 国内揄拍国产精品人妻在线| 亚洲av第一区精品v没综合| 亚洲欧美日韩高清专用| 特大巨黑吊av在线直播| 亚洲一区二区三区不卡视频| 国产精品影院久久| 少妇的逼好多水| 欧美成人a在线观看| 成人鲁丝片一二三区免费| 久久久久久久久中文| 国产一区二区三区视频了| 亚洲激情在线av| 欧美成人a在线观看| 人妻久久中文字幕网| 久久久久免费精品人妻一区二区| 我的老师免费观看完整版| 国产精品,欧美在线| 亚洲五月天丁香| 成熟少妇高潮喷水视频| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久久久亚洲av鲁大| 久久精品影院6| 在现免费观看毛片| 亚洲精品乱码久久久v下载方式| 老司机深夜福利视频在线观看| 国产三级中文精品| 女生性感内裤真人,穿戴方法视频| 成人av在线播放网站| 精品久久久久久久久久免费视频| 国产老妇女一区| 他把我摸到了高潮在线观看| 欧美zozozo另类| 精品日产1卡2卡| 国产免费男女视频| 自拍偷自拍亚洲精品老妇| 亚洲最大成人av| 能在线免费观看的黄片| 国产精品野战在线观看| 久久热精品热| 亚洲国产欧美人成| 一级作爱视频免费观看| 人妻丰满熟妇av一区二区三区| 小说图片视频综合网站| 中文字幕人成人乱码亚洲影| 精品免费久久久久久久清纯| 男人狂女人下面高潮的视频| 亚洲av中文字字幕乱码综合| 91久久精品电影网| 精品久久国产蜜桃| 亚洲国产色片| 国产精品野战在线观看| 午夜福利高清视频| 丰满乱子伦码专区| 国内精品久久久久久久电影| 91九色精品人成在线观看| 亚洲成人精品中文字幕电影| 深爱激情五月婷婷| 99热这里只有是精品在线观看 | 久久久久久久久久黄片| 成年人黄色毛片网站| 观看免费一级毛片| 亚洲精品在线美女| 国产三级黄色录像| 日韩 亚洲 欧美在线| 国产亚洲精品综合一区在线观看| 99国产极品粉嫩在线观看| 18禁在线播放成人免费| 中文字幕精品亚洲无线码一区| 国产亚洲精品av在线| 亚洲精品一区av在线观看| 两个人视频免费观看高清| 久久精品国产亚洲av天美| 亚洲第一区二区三区不卡| 变态另类丝袜制服| 国产真实乱freesex| 美女被艹到高潮喷水动态| 大型黄色视频在线免费观看| 国产一区二区三区在线臀色熟女| 99热6这里只有精品| 一本综合久久免费| 精品久久久久久,| 啪啪无遮挡十八禁网站| 中文字幕人妻熟人妻熟丝袜美| 欧美三级亚洲精品| 国产一区二区在线av高清观看| 首页视频小说图片口味搜索| 婷婷亚洲欧美| 香蕉av资源在线| 久久久久九九精品影院| 看免费av毛片| 亚洲精品亚洲一区二区| 国内毛片毛片毛片毛片毛片| 日本在线视频免费播放| 精品国产亚洲在线| 麻豆国产97在线/欧美| 久久99热6这里只有精品| 免费无遮挡裸体视频| 日本精品一区二区三区蜜桃| 精品久久久久久成人av| 午夜福利在线在线| 欧美成人性av电影在线观看| 99久国产av精品| 人妻制服诱惑在线中文字幕| 最好的美女福利视频网| 日韩欧美免费精品| 亚洲精品456在线播放app | 高潮久久久久久久久久久不卡| 国产一级毛片七仙女欲春2| 欧美黄色片欧美黄色片| 男人和女人高潮做爰伦理| 精品国内亚洲2022精品成人| 老熟妇仑乱视频hdxx| 久久久久久久午夜电影| 午夜激情福利司机影院| 欧美激情久久久久久爽电影| 男女下面进入的视频免费午夜| 亚洲国产精品成人综合色| 精品久久久久久成人av| av福利片在线观看| 五月玫瑰六月丁香| 亚洲人成电影免费在线| 99久久99久久久精品蜜桃| 色在线成人网| 亚洲av一区综合| 免费看日本二区| 亚洲av二区三区四区| 在线免费观看的www视频| 麻豆成人午夜福利视频| 亚洲中文日韩欧美视频| 一区二区三区高清视频在线| 亚洲无线在线观看| 性欧美人与动物交配| 成人国产综合亚洲| 美女黄网站色视频| 最近在线观看免费完整版| 赤兔流量卡办理| 最好的美女福利视频网| 女人十人毛片免费观看3o分钟| 亚洲美女搞黄在线观看 | 淫秽高清视频在线观看| 国产视频内射| 久久久久国产精品人妻aⅴ院| 亚洲精品影视一区二区三区av| 中国美女看黄片| 啦啦啦观看免费观看视频高清| 亚洲精品亚洲一区二区| 天美传媒精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 亚洲第一区二区三区不卡| 国产黄a三级三级三级人| 两人在一起打扑克的视频| 特级一级黄色大片| www日本黄色视频网| 性色avwww在线观看| 可以在线观看毛片的网站| 国内毛片毛片毛片毛片毛片| 日韩欧美在线二视频| 欧美日韩黄片免| 国产一区二区亚洲精品在线观看| 国产精品亚洲av一区麻豆| 黄色日韩在线| 88av欧美| 国产亚洲欧美98| 亚洲专区国产一区二区| 精品久久久久久久久av| 国产成人欧美在线观看| 国产老妇女一区| 成人特级av手机在线观看| 久久久国产成人精品二区| 小说图片视频综合网站| 久久久久久国产a免费观看| 国产精品久久久久久久久免 | 日本成人三级电影网站| 床上黄色一级片| 国产 一区 欧美 日韩| 宅男免费午夜| 中文在线观看免费www的网站| 高清日韩中文字幕在线| 色av中文字幕| 午夜免费成人在线视频| 97热精品久久久久久| 老熟妇乱子伦视频在线观看| 制服丝袜大香蕉在线| 欧美日本视频| a级毛片a级免费在线| 色av中文字幕| 亚洲精品成人久久久久久| 能在线免费观看的黄片| 久久精品人妻少妇| 午夜激情欧美在线| 我的老师免费观看完整版| 最近最新免费中文字幕在线| 国产精品亚洲美女久久久| 哪里可以看免费的av片| 亚洲国产精品999在线| 国产精品久久久久久人妻精品电影| 观看免费一级毛片| 深夜a级毛片| 国产成人av教育| 亚洲国产高清在线一区二区三| 亚洲av熟女| 久久久久久久久中文| 天堂动漫精品| av黄色大香蕉| 一a级毛片在线观看| 99久久精品国产亚洲精品| 亚洲电影在线观看av| а√天堂www在线а√下载| 美女被艹到高潮喷水动态| 久久婷婷人人爽人人干人人爱| 桃色一区二区三区在线观看| 日韩免费av在线播放| 午夜福利视频1000在线观看| 99久久九九国产精品国产免费| 制服丝袜大香蕉在线| 中文字幕av在线有码专区| 亚洲18禁久久av| 国产av一区在线观看免费| 欧美在线黄色| 色视频www国产| 日韩欧美国产一区二区入口| 在现免费观看毛片| 欧美最新免费一区二区三区 | 色播亚洲综合网| 中文字幕高清在线视频| 精品乱码久久久久久99久播| 免费黄网站久久成人精品 | 亚洲国产色片| 精品国产亚洲在线| 欧美黑人巨大hd| 老司机深夜福利视频在线观看| 精品久久国产蜜桃| 99热这里只有是精品50| 亚洲激情在线av| 91麻豆av在线| 男人舔奶头视频| 精品一区二区三区视频在线| 国产精品亚洲一级av第二区| or卡值多少钱| 亚洲精品久久国产高清桃花| 伦理电影大哥的女人| 超碰av人人做人人爽久久| 午夜视频国产福利| 少妇的逼好多水| 99热这里只有是精品50| 香蕉av资源在线| 亚洲av成人精品一区久久| 变态另类丝袜制服| 国产在视频线在精品| 精品人妻偷拍中文字幕| 午夜福利在线观看吧| 一个人看视频在线观看www免费| 亚洲av电影在线进入| 午夜免费男女啪啪视频观看 | 欧美高清性xxxxhd video| 丰满的人妻完整版| 男人狂女人下面高潮的视频| 97人妻精品一区二区三区麻豆| 亚洲第一欧美日韩一区二区三区| 18禁黄网站禁片免费观看直播| 精品午夜福利视频在线观看一区| 搡老岳熟女国产| 欧美另类亚洲清纯唯美| 草草在线视频免费看| 国产 一区 欧美 日韩| 欧美精品国产亚洲| 天堂av国产一区二区熟女人妻| a级毛片a级免费在线| 欧美日韩国产亚洲二区| 国产高清有码在线观看视频| 国内精品一区二区在线观看| 露出奶头的视频| 久久久久免费精品人妻一区二区| 嫩草影院精品99| 又爽又黄a免费视频| 五月玫瑰六月丁香| a级一级毛片免费在线观看| 露出奶头的视频| 熟妇人妻久久中文字幕3abv| 他把我摸到了高潮在线观看| 夜夜躁狠狠躁天天躁| 亚洲av日韩精品久久久久久密| 亚洲国产日韩欧美精品在线观看| 国产一区二区三区在线臀色熟女| 一本综合久久免费| 一区二区三区四区激情视频 | 全区人妻精品视频| 久久九九热精品免费| 一进一出抽搐gif免费好疼| 午夜激情福利司机影院| 午夜免费成人在线视频| 成年女人毛片免费观看观看9| 人人妻,人人澡人人爽秒播| x7x7x7水蜜桃| 18禁黄网站禁片免费观看直播| 亚洲人成网站在线播| 国产精品免费一区二区三区在线| av专区在线播放| 亚洲性夜色夜夜综合| 人妻制服诱惑在线中文字幕| 男女视频在线观看网站免费| 亚洲成人免费电影在线观看| 久99久视频精品免费| 色在线成人网| 丰满人妻熟妇乱又伦精品不卡| 深夜精品福利| 9191精品国产免费久久| 69人妻影院| 国产精品久久久久久久电影|