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

    信息熵控制的流場動態(tài)間距流線放置算法

    2021-12-08 00:20:46秦緒佳程宇軒左少華鄭紅波張美玉
    小型微型計算機系統(tǒng) 2021年12期
    關(guān)鍵詞:流線信息熵流場

    秦緒佳,程宇軒,左少華,鄭紅波,張美玉

    (浙江工業(yè)大學(xué) 計算機科學(xué)與技術(shù)學(xué)院,杭州 310023) E-mail:zhb@zjut.edu.cn

    1 引 言

    流場可視化是科學(xué)計算可視化領(lǐng)域的一個重要且最具挑戰(zhàn)性的研究熱點,涉及到對流場數(shù)據(jù)的轉(zhuǎn)換,繪制和操作,流場可視化領(lǐng)域應(yīng)用范圍廣泛,因此受到廣泛關(guān)注.流場存在于我們現(xiàn)實生活中的任何地方,例如水流,氣流等物理現(xiàn)象,我們無法直接觀察流場的內(nèi)部結(jié)構(gòu),因此需要借助科學(xué)的方法對流場進(jìn)行特殊處理,通過可視化技術(shù)將其復(fù)雜的數(shù)據(jù)結(jié)構(gòu)以一種簡潔直觀的方式展示,有效地獲得流場的內(nèi)部規(guī)律,這就是流場可視化作為可視化研究中重要領(lǐng)域的意義.流場可視化在許多學(xué)科領(lǐng)域例如空氣動力學(xué),流體力學(xué),氣象學(xué),海洋科學(xué)等有著重要的作用.

    流線可視化方法是使用流線來表達(dá)流場中的信息,可以準(zhǔn)確直觀地反映出流場的結(jié)構(gòu)特征,被廣泛應(yīng)用于許多工程領(lǐng)域.近些年來,隨著工程需求的不斷增強,許多研究人員對流線可視化方法進(jìn)行了大量研究,并在實踐中取得了很大的進(jìn)展與成果.Turk等[1]提出了一種基于圖像引導(dǎo)的流線放置方法,首先定義一個能量函數(shù),通過該函數(shù)對流線進(jìn)行一系列操作,然后反復(fù)迭代直到能量函數(shù)收斂,就能的得到分布均勻的流線,但該方法計算時間過長,效率較低,產(chǎn)生的流線中也存在較多無效值的短流線.Jobard等[2]提出了一種均勻放置算法,該方法先隨機生成一條初始流線,并在流線周圍放置種子點,然后計算種子點到已生成的流線的歐氏距離,判斷歐氏距離是否大于所給定的距離閾值,若大于則是有效種子點,可以生成新流線,若小于則拋棄該種子點,流線停止生成,該方法在所有流線的生成過程中不斷地檢測并判斷已放置的流線之間的距離,使流線能均勻分布,但計算時間較長,運行速度較慢,之后作者對該方法進(jìn)行了優(yōu)化,加快了運行速度并推向了非穩(wěn)定場[3].Mebarki等[4]提出了一種最遠(yuǎn)種子點放置算法,該方法每次生成流線時選擇距離其他流線最遠(yuǎn)的空白區(qū)域,結(jié)果是使得生成的流線較長,具有較好的連續(xù)性.Verma等[5]提出一種基于流引導(dǎo)的種子點放置算法,首先計算流場的特征點,然后選取適合的模板放置種子點,生成的流線能更好地表現(xiàn)流場的特征類型.Chen等[6]提出了一種相似度引導(dǎo)的流線放置方法,在計算流線及種子點之間的歐氏距離外又引入幾何體的相似性,通過結(jié)合這兩種方法度量不同流線間的相似度距離,能有效地判斷流場中冗余的流線并去除.郭雨蒙等[7]在Chen等提出的方法上增加了并行相似度引導(dǎo)的方法,提高了可視化效率.孔龍星等[8]提出了一種特征保持的視點相關(guān)三維矢量場流線簡化方法,通過對流線視覺效果進(jìn)行度量,有效地簡化了流線集,保持了流場特征且具有較好的視覺效果.魯大營等[9]提出三維流場的流線提取算法,采用了迭代最鄰近點(ICP)和K均值聚類,有效地簡化了三維流場的流線.巴振宇等[10]提出一種基于特征信息種子點選取的多層次流線可視化,在計算信息熵時對流場進(jìn)行區(qū)域劃分,加快了特征點的計算速度,并在放置種子點時選擇合適的特征模板.Han等[11]提出了一種基于深度學(xué)習(xí)的矢量場重構(gòu)方法,通過跟蹤原始矢量場中每個時間步的流線,并輸入到兩階段矢量場重建方法中,有效地簡化了矢量場.高天成等[12]提出了一種基于CVT的動態(tài)流場可視化方法,通過繼承不同時間步的生成元提高了計算效率,同時引入短流線解決了流線相交的問題,提高了可視化效果.

    本文主要采用了幾何可視化中基于流線的可視化方法來研究二維流場,傳統(tǒng)的二維流場可視化方法中使用固定間距放置流線,雖然能達(dá)到均勻布置流線的可視化效果,但卻依賴于固定間距的大小,間距過小會導(dǎo)致流線分布較密集,過大則容易丟失特征信息.本文提出了基于信息熵控制的流場動態(tài)間距流線放置算法,在傳統(tǒng)算法的基礎(chǔ)上引入信息熵來動態(tài)計算流線間隔,根據(jù)流場不同區(qū)域的不同特征合理布置種子點,通過該算法可以實現(xiàn)在流場的重要區(qū)域放置密集的流線,突出特征細(xì)節(jié),在其他非重要區(qū)域放置稀疏的流線,減少不必要的計算,提高可視化效率.

    2 流線生成方式

    流線是速度矢量場的場線,流線上任意一點的速度矢量都與該流線相切,流線方程定義如公式(1)所示:

    (1)

    其中μ(τ)代表點的位置,ν(μ(τ))表示點μ(τ)在流線上的切線方向,τ可以表示時間、弧長等參數(shù),則流線的求解公式如公式(2)所示:

    (2)

    流線生成過程如圖1所示,使用公式(1)定義流線方程,通過公式(2)求解,最終生成流線圖.

    圖1 流線生成示意圖Fig.1 Diagram of streamline generation

    流場數(shù)據(jù)集中的數(shù)據(jù)通常是離散的,在沒有解析表達(dá)式的情況下不能直接通過積分來計算流線,一般采用離散數(shù)值積分方法使流線逐步生長.

    常見的積分方法有一階歐拉法和四階龍格庫塔法,不同的方法生成的流線精度不同.其中歐拉法計算速度較快,但生成的流線精度較低;四階龍格庫塔計算量較大,速度較慢,但生成的流線精度較高.

    歐拉方法計算公式如公式(3)所示,其中Δt是積分步長,pi是流線上當(dāng)前點,pi+1是流線追蹤的下一點,ν(pi)是pi處的矢量值.

    pi+1=pi+Δt×ν(pi)

    (3)

    四階龍格庫塔方法計算公式如公式(4)所示,公式中使用了k1-k4這4個采樣點來計算,因此生成的流線精度更高.

    (4)

    3 信息熵控制的動態(tài)間隔流線放置算法

    3.1 流場信息熵定義

    (5)

    信息熵度量了一個系統(tǒng)的復(fù)雜程度,如果這個系統(tǒng)越混亂,它的復(fù)雜程度越高,不確定的信息也越多,信息熵就越大.反之這個系統(tǒng)越有序,復(fù)雜程度越低,不確定信息也越少,它的信息熵就越小.從香農(nóng)計算公式可以看出,當(dāng)變量x的每個可能值出現(xiàn)的概率p相等時,變量x的信息熵H(x)達(dá)到最大值,表示此時出現(xiàn)不同情況的種類最多,系統(tǒng)最復(fù)雜,不確定性最大.

    通過引入信息熵可以量化流場中矢量信息變化的程度,研究者可以借此有效地確定流場中變化程度不同的區(qū)域.Xu[13]等人將信息論中的信息熵概念引入到流場可視化中,并提出了一種流場信息熵的定義與信息熵計算的理論框架.首先將二維流場中的矢量方向區(qū)間(0°,360°)劃分成n個同等大小的子區(qū)間,記為xi,i∈1,2,…,n,然后計算流場中每個矢量的方向,并根據(jù)計算后的方向?qū)⒚總€矢量放入對應(yīng)的子區(qū)間中,通過統(tǒng)計每個子區(qū)間中矢量數(shù)量Si,i∈1,2,…,n,獲得一個矢量分布統(tǒng)計直方圖,則可以通過公式(6)計算每個子區(qū)間出現(xiàn)的概率pi,然后根據(jù)公式(5)可計算出流場的信息熵值,通過該值可以反映流場的穩(wěn)定程度.

    (6)

    迭代數(shù)據(jù)信息直到條件熵值收斂并迭代布置種子點,直到滿意的結(jié)果,這樣生成的流線具有較好的可視化效果,并且該框架能夠有效地可視化二維和三維流動數(shù)據(jù).Cheng等[14]在上述方法中加入了對流線的相似分組,優(yōu)化了遮擋問題.黃冬梅等[15]提出一種基于信息熵種子點選取的流線可視化方法,通過基于貪婪策略和蒙特卡洛的兩種種子點選取方法,優(yōu)化了流線生成的數(shù)量以及分布,解決了流線過多導(dǎo)致的遮擋與雜亂問題,具有良好的可視化效果與顯著的可用性.

    3.2 局部熵場計算

    牛嬋等[16]提出了一種基于局部最大熵值的特征點檢測方法,解決在提取數(shù)據(jù)時特征點被遺漏的問題,其次又提出了一種基于特征區(qū)域信息熵值變化規(guī)律的界定算法,提高流場特征信息提取的準(zhǔn)確度.通過計算流場的信息熵可以對流場整體的穩(wěn)定程度進(jìn)行量化度量,但為了突出流場中局部區(qū)域的變化程度,還需要將流場進(jìn)行區(qū)域劃分,分別計算局部位置的熵值.在流場中取多個網(wǎng)格點,計算以其為中心的一個L×L鄰域的信息熵作為該點處的熵值,通過計算所有網(wǎng)格點的熵值,可以得到一個與原流場數(shù)據(jù)規(guī)模一樣的熵場數(shù)據(jù).

    計算局部信息熵值需要給定鄰域范圍和區(qū)間數(shù)量.選定的鄰域范圍過小會導(dǎo)致該區(qū)域的穩(wěn)定程度較模糊,鄰域范圍過大會導(dǎo)致各區(qū)域相似度較高,網(wǎng)格點的熵值差距較小,計算量較大且影響計算結(jié)果的準(zhǔn)確性.對于區(qū)間數(shù)量,選的過小導(dǎo)致計算結(jié)果不準(zhǔn)確,相似度和偶然性較高,過多會增加熵場的計算量,使得效率下降.根據(jù)測試分析以及經(jīng)驗總結(jié),最合適的鄰域范圍L為13,區(qū)間數(shù)量為60,詳細(xì)計算過程如下:

    1)計算矢量方向角

    獲取二維流場數(shù)據(jù)的u,v分量并使用公式(7)計算其方向角度dir.

    (7)

    2)建立概率統(tǒng)計模型

    在流場中的局部鄰域建立統(tǒng)計模型,通過模型獲得概率密度函數(shù)計算該區(qū)域的信息熵.流場中矢量方向在(0°,360°)之間,將該范圍等分為60個區(qū)間,使用公式(8)計算dir所在區(qū)間I.

    (8)

    劃分完區(qū)間后需要根據(jù)每個區(qū)間內(nèi)的矢量個數(shù)統(tǒng)計局部鄰域內(nèi)的區(qū)間數(shù).設(shè)i為流場中的網(wǎng)格點,取以點i為中心的局部鄰域,根據(jù)統(tǒng)計直方圖獲取該鄰域中矢量數(shù)不為0的區(qū)間總數(shù),設(shè)為bin[i],i∈1,2,…,N,設(shè)bin[i][j],j∈1,2,…,bin[i]為統(tǒng)計直方圖中以點i為中心的局部鄰域范圍內(nèi)第j個區(qū)間中的矢量個數(shù),sumBin[i]表示以網(wǎng)格點i為中心的局部鄰域中矢量個數(shù),這里取為13×13,通過公式(9)可以計算出以點i為中心的局部鄰域中第j個區(qū)間的概率pj.

    (9)

    3)計算局部熵值生成熵場

    對每個網(wǎng)格點i通過上述方式計算概率pj,以此建立概率統(tǒng)計模型,則每個網(wǎng)格點所在的局部鄰域的信息熵可由公式(10)計算得出.

    (10)

    圖2是使用上述信息熵計算方法生成的熵場,并通過顏色映射直觀反映計算結(jié)果.圖2(a)、圖2(b)中左側(cè)為流場圖,右側(cè)為該流場數(shù)據(jù)生成的熵場圖,最右側(cè)為顏色映射條,越靠上的顏色代表的熵值越大,通過給不同熵值設(shè)定不同的顏色值來反映不同區(qū)域的熵值變化.通過對比可以看出在流場變化劇烈的區(qū)域熵值較大,如圖2中圓形框處的區(qū)域,對應(yīng)顏色映射中頂端的顏色;變化平緩的區(qū)域熵值較小,如圖2中矩形框處的區(qū)域,對應(yīng)顏色映射表中底端的顏色.

    圖2 熵場顏色映射結(jié)果Fig.2 Color mapping result of information entropy field

    3.3 流線間距控制

    在流場中生成過多的流線時,如果只是考慮到流線脫離網(wǎng)格邊界或速度減小到零就終止生長卻沒有對流線間的距離加以控制,其結(jié)果就是生成的流線相互之間較緊密,距離較小,視覺上表現(xiàn)得較擁擠.圖3是在無流線間距控制的情況下生成的流線效果,可以明顯地看到流線間距離較小,過于靠近,尤其是在流線的匯聚點,使得視覺效果較差.

    為了解決流線過于密集的問題,很多的流場都采用基于歐式距離的流線間距控制來生成流線,將流線間距作為流線的終止判定條件之一.給定一個常量ds為流線間距,設(shè)p為當(dāng)前流線生長過程中的下一個采樣點,計算p到其周圍某條已存在的流線的距離,因為流線是由一連串離散的采樣點組合而成,因此在計算距離時只需要計算p到流線上采樣點的距離,若p到其它流線上的某一采樣點距離小于ds,則認(rèn)定當(dāng)前正在生成的流線與已生成的流線間的距離過小,流線停止生長.

    然而隨著流場中生成的流線越來越多,這種方法在計算距離時所選取的現(xiàn)存流線上的采樣點沒有受到約束,會產(chǎn)生大量不必要的計算比較,導(dǎo)致算法效率較低.因此,需要在流場區(qū)域建立一個橫豎間隔為ds的虛擬網(wǎng)格,將已生成的流線的采樣點存儲在相應(yīng)的網(wǎng)格單元中,每次生成新流線的采樣點時,通過網(wǎng)格來建立約束條件.

    如圖4所示,上方與下方的方塊點所在的線條為已生成完畢的流線,中間的圓點所在的線條為正在生長的新流線,線上有已生成的采樣點,圓點p為當(dāng)前最新的采樣點,方塊點為已存在的流線上的采樣點,根據(jù)該區(qū)域建立的虛擬控制網(wǎng)格,在計算采樣點距離時,只需要對以p所在網(wǎng)格為中心,選取周圍8個網(wǎng)格中其它流線上的采樣點進(jìn)行距離計算,若存在某一個采樣點與p的距離小于給定的最小間距,則流線停止生長.

    圖3 無間距控制流線效果Fig.3 Streamlineeffectwithoutdistancecontrol圖4 距離控制網(wǎng)格Fig.4 Distancecontrolgrid圖5 間距控制流線效果Fig.5 Streamlineeffectwithdistancecontrol

    圖5是采用了間距控制后流線放置的效果,與圖3進(jìn)行對比,流線的種子點仍相同,并且可以明顯地看到原本流線過于密集的地方通過間距控制后,使得該處流線因距離較小而停止生長,優(yōu)化了流線的放置,在視覺上有較好的效果.

    3.4 信息熵控制的動態(tài)間距控制

    新生成的流線上的采樣點會通過流線間距控制算法計算與周圍已生成流線上的采樣點的距離,并與給定距離ds進(jìn)行判斷,從而控制新采樣點的生成,達(dá)到流線均勻放置的效果.但存在的問題是,這種給定的流線間距離ds是固定的,流線分布的位置完全依賴于流線間距ds的大小,若給定的ds較小,就會導(dǎo)致流線數(shù)量過大,間距過小,分布過于密集,使得算法運算效率低下;若ds較大,則流線分布較稀疏,容易丟失流場中特征區(qū)域的細(xì)節(jié),難以表現(xiàn)流場的所有特征,因此這種使用固定間距來控制流線生成的方法并不具有良好的適應(yīng)性.在前幾章中已經(jīng)詳細(xì)介紹了流場信息熵以及場熵如何計算,通過在流場中引入信息熵可以對流場穩(wěn)定程度進(jìn)行量化度量,可以突出流場中的特征細(xì)節(jié).本文使用信息熵對流線的間距進(jìn)行動態(tài)控制,通過計算每個區(qū)域的局部信息熵,并在計算流線間距時引入,使得流線間距不再是固定值而是根據(jù)流場特征進(jìn)行動態(tài)調(diào)整,根據(jù)流場信息熵的計算,流場中非重要區(qū)域的信息熵值較小,則該區(qū)域計算后的流線間距較大,生成的流線較稀疏,流向較穩(wěn)定;流場中重要特征區(qū)域的信息熵值較大,則該區(qū)域計算后的流線間距較小,生成流線較密集,流場中不同區(qū)域的流線密度不一致能夠使得特征區(qū)域更加突出,在視覺上有更好的效果.

    給定流線間距的最大值為dmax,最小值為dmin,在流場中建立大小為dmax×dmax的虛擬控制網(wǎng)格.在流線生成過程中,流線控制間距值ds根據(jù)采樣點所在區(qū)域的信息熵值來動態(tài)計算,計算如公式(11)所示.

    (11)

    其中emax為信息熵最大值,emin為信息熵最小值,ei為當(dāng)前點信息熵值,完整的信息熵控制的動態(tài)間距流線放置算法如表1所示.

    表1 信息熵控制的動態(tài)間距控制算法Table 1 Dynamic distance control algorithm based on information entropy

    4 實驗結(jié)果與分析

    實驗結(jié)果生成的不同流線間距的流場效果對比如圖6所示.圖6(a)-圖6(d)是使用相同數(shù)據(jù)生成的流場,其中圖6(a)、圖6(b)、圖6(c)采用了固定間距的流線生成算法,給定的間距依次增大.分析圖6(a)和圖6(b),可以看出在間距較小時,生成的流線排列過于緊密,較混亂,視覺效果上較差;如圖6(c)所示,間距逐漸增大后,可視化效果上得到了優(yōu)化,但丟失了流線特征區(qū)域的細(xì)節(jié),如流場中的渦流與鞍點處.圖6(d)則是采用了本文基于信息熵的動態(tài)流線間距放置生成的效果圖,在優(yōu)化可視化效果的基礎(chǔ)上,保持了流場中的特征細(xì)節(jié)與變化規(guī)律.

    圖6 不同流線間距的流場效果對比Fig.6 Comparison of flow field effects with different streamline spacing

    詳細(xì)的流線放置效果對比如圖7所示,圖7(a)是采用傳統(tǒng)的固定流線間距生成的效果圖,圖7(b)是采用本文算法生成的效果圖,圖7(b)中給定的流線間距最大值與圖7(a)中的固定間距一致.圖中圓形框處為流場的特征區(qū)域,通過對比可以明顯地看出圖7(b)中該區(qū)域生成的流線間距小于圖7(a),流線密度高于圖7(a),矩形框處為平穩(wěn)區(qū)域,在該區(qū)域圖7(b)的流線密度與圖7(a)基本保持一致.

    圖7 流線放置效果對比Fig.7 Comparison of streamline placement

    綜合實驗結(jié)果分析,在使用傳統(tǒng)的固定間距算法生成流線時,給定的固定間距較大則容易丟失重要區(qū)域細(xì)節(jié),固定間距較小則流線密度較大導(dǎo)致計算效率較低且可視化效果較差,通過本文基于信息熵的動態(tài)間距控制可以合理分配區(qū)域流線密度,在提高了流線生成的速率的同時又保持了流場重要特征.

    5 結(jié) 論

    流場可視化是可視化領(lǐng)域的重要課題,合理地放置流線能突出流場的特征與全局規(guī)律,流場中流線的間距是影響流線放置效果的重要因素.本文研究了二維流場的流線放置算法,傳統(tǒng)的均勻放置算法采用的固定的流線間隔,生成的流線具有較好地均勻分布效果,但采用的固定間隔較小時流線整體密度大,算法耗時;采用的固定間隔較大時流線整體密度低,容易丟失流場中特征區(qū)域的重要信息.本文提出了基于信息熵控制的動態(tài)間距流線放置算法,在流線生成過程中引入信息熵來度量流場區(qū)域的穩(wěn)定程度,并以此為依據(jù)動態(tài)地控制流線間距,使得流線分布具有較好的自適應(yīng)性,在變化劇烈的特征區(qū)域放置較多較密集的流線來完整地表達(dá)重要信息;在平穩(wěn)區(qū)域放置較少較稀疏的流線來提高可視化效率.

    猜你喜歡
    流線信息熵流場
    基于信息熵可信度的測試點選擇方法研究
    大型空冷汽輪發(fā)電機轉(zhuǎn)子三維流場計算
    幾何映射
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    基于信息熵的實驗教學(xué)量化研究
    電子測試(2017年12期)2017-12-18 06:35:48
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計分析
    一種基于信息熵的雷達(dá)動態(tài)自適應(yīng)選擇跟蹤方法
    基于信息熵的IITFN多屬性決策方法
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    密集導(dǎo)向流線算法研究
    男女床上黄色一级片免费看| 人人妻人人看人人澡| 成年免费大片在线观看| 日韩有码中文字幕| 51国产日韩欧美| 91在线精品国自产拍蜜月 | 校园春色视频在线观看| 国产精品久久久久久人妻精品电影| 亚洲激情在线av| 精品乱码久久久久久99久播| 色噜噜av男人的天堂激情| 波多野结衣高清无吗| 婷婷精品国产亚洲av在线| 国产免费一级a男人的天堂| 亚洲色图av天堂| 国产精品综合久久久久久久免费| 欧美激情久久久久久爽电影| 一个人看视频在线观看www免费 | 亚洲美女视频黄频| 国产欧美日韩一区二区三| 亚洲精品粉嫩美女一区| 成人特级av手机在线观看| 久久久久久大精品| 午夜亚洲福利在线播放| 99热只有精品国产| 午夜免费成人在线视频| 男人舔奶头视频| 亚洲精品在线观看二区| 麻豆成人av在线观看| 熟女电影av网| 日韩中文字幕欧美一区二区| 日日摸夜夜添夜夜添小说| 午夜老司机福利剧场| 最后的刺客免费高清国语| www.熟女人妻精品国产| 亚洲中文字幕日韩| 亚洲熟妇熟女久久| x7x7x7水蜜桃| 国产精品99久久久久久久久| 国产午夜精品久久久久久一区二区三区 | 亚洲国产欧洲综合997久久,| 精品一区二区三区视频在线 | 舔av片在线| 在线观看日韩欧美| 91九色精品人成在线观看| 波多野结衣高清无吗| 亚洲无线在线观看| 国产69精品久久久久777片| 久久久久久久久中文| 亚洲va日本ⅴa欧美va伊人久久| 国产爱豆传媒在线观看| 精品国产亚洲在线| 久久久久久九九精品二区国产| 久久久久国产精品人妻aⅴ院| 久久久久性生活片| 别揉我奶头~嗯~啊~动态视频| 美女黄网站色视频| 国产亚洲精品久久久com| 欧美成人一区二区免费高清观看| 国产精品日韩av在线免费观看| 国产亚洲av嫩草精品影院| 中文在线观看免费www的网站| 他把我摸到了高潮在线观看| 免费大片18禁| 日本黄色片子视频| 少妇高潮的动态图| 欧美黄色片欧美黄色片| 在线免费观看不下载黄p国产 | 免费在线观看成人毛片| 欧美精品啪啪一区二区三区| 国内揄拍国产精品人妻在线| 国产高清有码在线观看视频| 99久久成人亚洲精品观看| 在线观看66精品国产| 免费av毛片视频| www.www免费av| 中文字幕人妻熟人妻熟丝袜美 | 18禁黄网站禁片午夜丰满| 日日摸夜夜添夜夜添小说| 嫩草影视91久久| 久久久精品欧美日韩精品| 色视频www国产| 久久精品亚洲精品国产色婷小说| 日本 av在线| 亚洲精品久久国产高清桃花| 窝窝影院91人妻| 亚洲专区中文字幕在线| 18禁美女被吸乳视频| 老司机福利观看| 白带黄色成豆腐渣| 波多野结衣高清作品| 亚洲欧美日韩高清专用| 中文资源天堂在线| 国产真人三级小视频在线观看| 亚洲av熟女| 熟女少妇亚洲综合色aaa.| 三级男女做爰猛烈吃奶摸视频| 日韩 欧美 亚洲 中文字幕| 啪啪无遮挡十八禁网站| 久久久久久久亚洲中文字幕 | 亚洲av日韩精品久久久久久密| 精品无人区乱码1区二区| 国产伦一二天堂av在线观看| 国产蜜桃级精品一区二区三区| 男人的好看免费观看在线视频| 黑人欧美特级aaaaaa片| 此物有八面人人有两片| 日韩欧美免费精品| 蜜桃久久精品国产亚洲av| 在线观看免费视频日本深夜| 亚洲国产欧洲综合997久久,| 国产麻豆成人av免费视频| 欧美高清成人免费视频www| 亚洲自拍偷在线| 亚洲av熟女| 尤物成人国产欧美一区二区三区| 日本免费一区二区三区高清不卡| 嫩草影院精品99| 国产单亲对白刺激| 亚洲精品在线美女| 超碰av人人做人人爽久久 | 最好的美女福利视频网| 亚洲精华国产精华精| 黑人欧美特级aaaaaa片| 日韩成人在线观看一区二区三区| 久久精品91无色码中文字幕| 少妇的逼水好多| 亚洲自拍偷在线| 夜夜爽天天搞| 看黄色毛片网站| 搡老岳熟女国产| 中国美女看黄片| 欧美一级a爱片免费观看看| 精品人妻偷拍中文字幕| 天天躁日日操中文字幕| 精品久久久久久久久久免费视频| 国产午夜精品久久久久久一区二区三区 | aaaaa片日本免费| 欧美+日韩+精品| 国产伦精品一区二区三区视频9 | 国产午夜福利久久久久久| 19禁男女啪啪无遮挡网站| xxxwww97欧美| 全区人妻精品视频| 国产精品日韩av在线免费观看| 又爽又黄无遮挡网站| 黄色成人免费大全| 亚洲精品色激情综合| 日韩有码中文字幕| 成人高潮视频无遮挡免费网站| 欧美午夜高清在线| av片东京热男人的天堂| 亚洲国产高清在线一区二区三| 亚洲精品国产精品久久久不卡| 色视频www国产| 免费av不卡在线播放| 日本精品一区二区三区蜜桃| 免费看日本二区| 国产精品98久久久久久宅男小说| 欧美日韩中文字幕国产精品一区二区三区| 日本在线视频免费播放| 国产乱人视频| 亚洲精品久久国产高清桃花| 日韩欧美在线乱码| 国产精品久久久人人做人人爽| 99精品久久久久人妻精品| 特级一级黄色大片| 19禁男女啪啪无遮挡网站| 久久精品国产99精品国产亚洲性色| 麻豆成人av在线观看| 亚洲国产欧美人成| 制服人妻中文乱码| 久久国产乱子伦精品免费另类| 中国美女看黄片| 麻豆一二三区av精品| 五月伊人婷婷丁香| 18禁黄网站禁片午夜丰满| 又黄又爽又免费观看的视频| 岛国视频午夜一区免费看| 非洲黑人性xxxx精品又粗又长| 成人国产综合亚洲| 久久精品影院6| 免费看日本二区| 成人特级黄色片久久久久久久| 露出奶头的视频| 欧美另类亚洲清纯唯美| 99久久久亚洲精品蜜臀av| 欧美一区二区亚洲| 国产探花在线观看一区二区| 久久这里只有精品中国| 免费大片18禁| 色吧在线观看| 色哟哟哟哟哟哟| 亚洲国产欧洲综合997久久,| 十八禁人妻一区二区| 午夜久久久久精精品| 国产欧美日韩精品一区二区| 精品国内亚洲2022精品成人| 91九色精品人成在线观看| 国产精品98久久久久久宅男小说| 日本五十路高清| 琪琪午夜伦伦电影理论片6080| 亚洲国产色片| 亚洲在线自拍视频| 熟女电影av网| 99热只有精品国产| 高清毛片免费观看视频网站| 法律面前人人平等表现在哪些方面| 日本熟妇午夜| 成人三级黄色视频| 日日干狠狠操夜夜爽| 成熟少妇高潮喷水视频| 少妇高潮的动态图| 久久精品人妻少妇| АⅤ资源中文在线天堂| 欧洲精品卡2卡3卡4卡5卡区| 久久人人精品亚洲av| 久久久久九九精品影院| 757午夜福利合集在线观看| 久久精品影院6| 校园春色视频在线观看| 午夜免费成人在线视频| 成人av在线播放网站| 欧美日韩瑟瑟在线播放| 午夜福利免费观看在线| 国产精品一区二区免费欧美| 中文字幕人成人乱码亚洲影| 非洲黑人性xxxx精品又粗又长| 在线a可以看的网站| 露出奶头的视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 熟女电影av网| 色播亚洲综合网| 一区二区三区国产精品乱码| 亚洲av美国av| 国产av一区在线观看免费| 日韩欧美国产一区二区入口| 真人一进一出gif抽搐免费| 国产毛片a区久久久久| 女生性感内裤真人,穿戴方法视频| 法律面前人人平等表现在哪些方面| 欧美+亚洲+日韩+国产| 一区二区三区免费毛片| 欧美激情久久久久久爽电影| 偷拍熟女少妇极品色| 亚洲欧美日韩高清在线视频| 少妇丰满av| 日韩有码中文字幕| 国产一区二区三区视频了| 亚洲七黄色美女视频| 午夜视频国产福利| 国产av不卡久久| 国产精品久久电影中文字幕| 亚洲国产精品sss在线观看| 久久久国产成人免费| 网址你懂的国产日韩在线| 国产精品精品国产色婷婷| 99热这里只有是精品50| 丰满的人妻完整版| 亚洲国产中文字幕在线视频| 美女大奶头视频| 久99久视频精品免费| 久久国产精品影院| 久久亚洲真实| 成人国产综合亚洲| 亚洲七黄色美女视频| 一区二区三区国产精品乱码| 19禁男女啪啪无遮挡网站| 级片在线观看| 国产高清激情床上av| 久久久久性生活片| 国内少妇人妻偷人精品xxx网站| 国产精品久久视频播放| 久久精品国产综合久久久| 天堂√8在线中文| 高清在线国产一区| 国产精品av视频在线免费观看| 欧美色欧美亚洲另类二区| 人妻夜夜爽99麻豆av| 在线观看午夜福利视频| 成人国产综合亚洲| 最近最新中文字幕大全免费视频| 99久久99久久久精品蜜桃| www.999成人在线观看| 高清在线国产一区| 黄色女人牲交| 亚洲国产精品久久男人天堂| 男人舔女人下体高潮全视频| 国产精品亚洲av一区麻豆| 国产高清视频在线播放一区| 精品欧美国产一区二区三| av天堂中文字幕网| 久久天躁狠狠躁夜夜2o2o| 亚洲一区二区三区不卡视频| 亚洲性夜色夜夜综合| 成人一区二区视频在线观看| 欧美黑人欧美精品刺激| 国产aⅴ精品一区二区三区波| 搡老岳熟女国产| 99久久精品国产亚洲精品| 亚洲成人久久性| 日本一二三区视频观看| 成人特级黄色片久久久久久久| 男女视频在线观看网站免费| 激情在线观看视频在线高清| 亚洲av第一区精品v没综合| 午夜日韩欧美国产| 噜噜噜噜噜久久久久久91| 欧美最黄视频在线播放免费| 亚洲av美国av| 国产亚洲av嫩草精品影院| a级一级毛片免费在线观看| 高清日韩中文字幕在线| 精品99又大又爽又粗少妇毛片 | 99久久久亚洲精品蜜臀av| 一本久久中文字幕| 日韩中文字幕欧美一区二区| 美女cb高潮喷水在线观看| 亚洲色图av天堂| 成人午夜高清在线视频| 免费看光身美女| 精品久久久久久久毛片微露脸| 高潮久久久久久久久久久不卡| 高清日韩中文字幕在线| 国产精品 欧美亚洲| av国产免费在线观看| 亚洲欧美激情综合另类| 国产精品久久久久久亚洲av鲁大| 免费看十八禁软件| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 搡老熟女国产l中国老女人| 三级国产精品欧美在线观看| 九九热线精品视视频播放| 人妻久久中文字幕网| 亚洲精品一区av在线观看| 成人性生交大片免费视频hd| 一进一出抽搐gif免费好疼| 青草久久国产| 老汉色∧v一级毛片| a在线观看视频网站| av天堂中文字幕网| 18禁裸乳无遮挡免费网站照片| 亚洲精品一卡2卡三卡4卡5卡| 51午夜福利影视在线观看| 好看av亚洲va欧美ⅴa在| 噜噜噜噜噜久久久久久91| 亚洲中文日韩欧美视频| 黄色视频,在线免费观看| 午夜免费男女啪啪视频观看 | 国内精品一区二区在线观看| 亚洲av熟女| 国内精品一区二区在线观看| 99久久无色码亚洲精品果冻| 一区二区三区国产精品乱码| 五月伊人婷婷丁香| 国产精品久久视频播放| 天天躁日日操中文字幕| 精品国内亚洲2022精品成人| 午夜免费观看网址| 悠悠久久av| 九色成人免费人妻av| 老司机福利观看| 91字幕亚洲| 日本一二三区视频观看| 又黄又爽又免费观看的视频| 国产成人福利小说| 欧美中文综合在线视频| 最好的美女福利视频网| 久久久国产成人免费| 网址你懂的国产日韩在线| 九色成人免费人妻av| 精品午夜福利视频在线观看一区| 国内精品一区二区在线观看| 国产主播在线观看一区二区| 少妇人妻精品综合一区二区 | 在线十欧美十亚洲十日本专区| 五月伊人婷婷丁香| 久久久久九九精品影院| 欧美乱色亚洲激情| 亚洲av成人av| 99热这里只有精品一区| 欧美+亚洲+日韩+国产| 国产麻豆成人av免费视频| 欧美+亚洲+日韩+国产| 国产中年淑女户外野战色| 国产精品影院久久| 免费看日本二区| 五月伊人婷婷丁香| 成人18禁在线播放| 97超视频在线观看视频| 成人欧美大片| 深爱激情五月婷婷| 少妇高潮的动态图| 高清毛片免费观看视频网站| 悠悠久久av| 亚洲成人久久爱视频| 此物有八面人人有两片| 少妇的丰满在线观看| 人人妻人人澡欧美一区二区| 免费看日本二区| 午夜福利欧美成人| 在线播放国产精品三级| 亚洲熟妇中文字幕五十中出| 波多野结衣巨乳人妻| 天天添夜夜摸| 久久久久国产精品人妻aⅴ院| 啪啪无遮挡十八禁网站| 成人特级黄色片久久久久久久| 午夜精品在线福利| 日韩欧美国产一区二区入口| 久久久国产成人免费| 成人永久免费在线观看视频| 91九色精品人成在线观看| 国产高清激情床上av| 欧美色视频一区免费| 欧美绝顶高潮抽搐喷水| 亚洲精品国产精品久久久不卡| 日本三级黄在线观看| 免费在线观看影片大全网站| 十八禁人妻一区二区| 激情在线观看视频在线高清| 制服丝袜大香蕉在线| 一区二区三区国产精品乱码| 99在线人妻在线中文字幕| xxx96com| 一进一出好大好爽视频| 欧美成人免费av一区二区三区| 亚洲精品亚洲一区二区| 俺也久久电影网| 免费无遮挡裸体视频| 欧美一区二区国产精品久久精品| 国产色爽女视频免费观看| 久久婷婷人人爽人人干人人爱| 欧美3d第一页| 日韩欧美国产一区二区入口| 好男人在线观看高清免费视频| 三级男女做爰猛烈吃奶摸视频| 欧美激情在线99| 又黄又粗又硬又大视频| 色综合站精品国产| 亚洲在线自拍视频| 国产精品一及| 国产精品亚洲av一区麻豆| 精品一区二区三区视频在线观看免费| 日本 欧美在线| 欧美一区二区亚洲| 久久久久精品国产欧美久久久| 丁香欧美五月| 在线观看日韩欧美| 精品国产超薄肉色丝袜足j| 一区二区三区免费毛片| 成人高潮视频无遮挡免费网站| 99国产综合亚洲精品| 午夜视频国产福利| 哪里可以看免费的av片| 久久久久精品国产欧美久久久| 美女免费视频网站| 亚洲成a人片在线一区二区| 亚洲精品亚洲一区二区| 午夜精品久久久久久毛片777| 久久精品91蜜桃| 久久久久国内视频| 国产精品98久久久久久宅男小说| av视频在线观看入口| 欧美日韩瑟瑟在线播放| 精品免费久久久久久久清纯| 国产伦在线观看视频一区| 九九热线精品视视频播放| 窝窝影院91人妻| 久久精品影院6| 天堂网av新在线| 欧美黄色片欧美黄色片| 天堂影院成人在线观看| 黄色女人牲交| 亚洲国产精品久久男人天堂| 男女下面进入的视频免费午夜| 国产v大片淫在线免费观看| 日本五十路高清| 很黄的视频免费| 人人妻,人人澡人人爽秒播| 在线观看免费午夜福利视频| 欧美日韩黄片免| 高清毛片免费观看视频网站| 女同久久另类99精品国产91| 久久精品91无色码中文字幕| 波野结衣二区三区在线 | 青草久久国产| 动漫黄色视频在线观看| 中国美女看黄片| 欧美日韩综合久久久久久 | 又紧又爽又黄一区二区| 中亚洲国语对白在线视频| 成年版毛片免费区| 99久久无色码亚洲精品果冻| 国产精品一区二区三区四区免费观看 | 欧美zozozo另类| 日韩国内少妇激情av| 岛国在线免费视频观看| 国产精品乱码一区二三区的特点| 超碰av人人做人人爽久久 | 黄色女人牲交| 丰满人妻熟妇乱又伦精品不卡| 男女做爰动态图高潮gif福利片| 日韩欧美在线二视频| 黄色女人牲交| 午夜久久久久精精品| 97超视频在线观看视频| 亚洲在线观看片| 中文字幕人妻丝袜一区二区| 啦啦啦观看免费观看视频高清| 他把我摸到了高潮在线观看| 欧美极品一区二区三区四区| 婷婷精品国产亚洲av在线| 国产成人av教育| 日韩免费av在线播放| x7x7x7水蜜桃| 国产精品久久久人人做人人爽| 麻豆一二三区av精品| 色播亚洲综合网| 久久久久久人人人人人| 亚洲avbb在线观看| 午夜精品久久久久久毛片777| 欧美三级亚洲精品| 国产97色在线日韩免费| 日韩欧美三级三区| 久久久久久久久中文| 最近最新免费中文字幕在线| 欧美+亚洲+日韩+国产| 日本三级黄在线观看| 色综合站精品国产| 黄色成人免费大全| 精品无人区乱码1区二区| 久久99热这里只有精品18| 一区二区三区国产精品乱码| 午夜福利在线观看吧| 97碰自拍视频| 国产在视频线在精品| 90打野战视频偷拍视频| 精品不卡国产一区二区三区| 97超级碰碰碰精品色视频在线观看| 午夜免费观看网址| 国产午夜精品论理片| 午夜福利视频1000在线观看| 欧美av亚洲av综合av国产av| 亚洲av成人精品一区久久| 性色av乱码一区二区三区2| 国产一区二区三区视频了| 中出人妻视频一区二区| 男女床上黄色一级片免费看| 亚洲欧美日韩东京热| 欧美色视频一区免费| 一a级毛片在线观看| 天堂av国产一区二区熟女人妻| 国产精品久久久久久人妻精品电影| 日韩免费av在线播放| 国产久久久一区二区三区| 久久久久久大精品| 久久香蕉精品热| 欧美3d第一页| 亚洲精品美女久久久久99蜜臀| 日本熟妇午夜| 精品熟女少妇八av免费久了| 男女午夜视频在线观看| 校园春色视频在线观看| 国产精品,欧美在线| 亚洲人与动物交配视频| 丁香六月欧美| 叶爱在线成人免费视频播放| 免费av不卡在线播放| 18禁在线播放成人免费| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久久久人妻精品电影| h日本视频在线播放| 日韩av在线大香蕉| 精品久久久久久成人av| 欧美最新免费一区二区三区 | 91麻豆精品激情在线观看国产| 天堂网av新在线| 丰满乱子伦码专区| 老熟妇乱子伦视频在线观看| 精品一区二区三区视频在线 | 亚洲精品乱码久久久v下载方式 | 亚洲精品456在线播放app | 中文资源天堂在线| 天天躁日日操中文字幕| 色哟哟哟哟哟哟| 亚洲一区高清亚洲精品| 午夜免费观看网址| 在线国产一区二区在线| 搡老妇女老女人老熟妇| 在线观看免费视频日本深夜| 天天躁日日操中文字幕| 午夜福利高清视频| 99久久综合精品五月天人人| 免费电影在线观看免费观看| 亚洲中文日韩欧美视频| 亚洲国产精品成人综合色| 欧美三级亚洲精品| xxxwww97欧美| 精品国产美女av久久久久小说| 国产探花在线观看一区二区| 天美传媒精品一区二区| 亚洲中文字幕日韩| xxxwww97欧美| 少妇人妻一区二区三区视频| 天堂影院成人在线观看| 国产精品日韩av在线免费观看| 757午夜福利合集在线观看| 亚洲最大成人手机在线| 免费看a级黄色片| 国产在视频线在精品| 丰满乱子伦码专区| 黄色日韩在线| 亚洲国产色片| av欧美777|