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

    對(duì)流擴(kuò)散方程的多尺度解耦小波算法研究

    2013-12-07 05:57:38汪友明沈建冬
    關(guān)鍵詞:對(duì)流小波尺度

    汪友明,吳 青,沈建冬

    (1.西安郵電大學(xué)自動(dòng)化學(xué)院, 陜西西安 710121;2.大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧大連 116023)

    對(duì)流擴(kuò)散方程的多尺度解耦小波算法研究

    汪友明1,2,吳 青1,沈建冬1

    (1.西安郵電大學(xué)自動(dòng)化學(xué)院, 陜西西安 710121;2.大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧大連 116023)

    針對(duì)傳統(tǒng)有限元方法在求解對(duì)流擴(kuò)散問(wèn)題時(shí)常會(huì)出現(xiàn)的數(shù)值震蕩和數(shù)值耗散等缺點(diǎn),提出一種對(duì)流擴(kuò)散方程的尺度解耦小波求解方法。介紹第二代小波多分辨分析,推導(dǎo)有限元多分辨空間的兩尺度關(guān)系,提出對(duì)流擴(kuò)散方程的多尺度計(jì)算框架。推導(dǎo)對(duì)流擴(kuò)散方程的解耦條件,并利用提升方案構(gòu)造多尺度解耦小波。提出多尺度解耦小波算法,該方法通過(guò)向求解域添加解耦小波,逐步逼近問(wèn)題精確解。數(shù)值算例證明,解耦小波是一種求解對(duì)流擴(kuò)散方程性能優(yōu)良的小波基。

    對(duì)流擴(kuò)散方程;多尺度計(jì)算;解耦小波

    對(duì)流擴(kuò)散問(wèn)題是電子、流體力學(xué)、傳熱學(xué)等學(xué)科以及水利工程、環(huán)境工程、化工、冶金、航空等應(yīng)用領(lǐng)域中經(jīng)常遇到的基礎(chǔ)問(wèn)題之一。傳統(tǒng)有限元方法在求解對(duì)流擴(kuò)散問(wèn)題時(shí)常會(huì)出現(xiàn)數(shù)值震蕩和數(shù)值耗散等缺點(diǎn)[1-2]。小波數(shù)值方法是近年來(lái)出現(xiàn)的新興數(shù)值算法,其優(yōu)點(diǎn)是多尺度、多分辨、緊支撐等特性,可根據(jù)求解域的變化梯度或奇異性任意改變小波的分析尺度或提高逼近階,而小波的緊支性、正則性使得小波基函數(shù)聚焦到待求解問(wèn)題的任意細(xì)節(jié),從而實(shí)現(xiàn)工程問(wèn)題的高效求解[3-6]。近年來(lái)將小波方法與傳統(tǒng)的數(shù)值方法相結(jié)合出現(xiàn)了很多性能優(yōu)越的新算法,包括小波Galerkin法[7-8],小波配置法[9],小波Petrov-Galerkin法[10],小波最小二乘法[11]和小波有限元法[12-13]等。

    以伸縮、平移為基礎(chǔ)的第一代小波無(wú)法在有限、不規(guī)則、非結(jié)構(gòu)化的網(wǎng)格上構(gòu)造,且在求解工程問(wèn)題前所選定小波的特性已無(wú)法改變,這些缺點(diǎn)常難以滿足各種工程問(wèn)題的多尺度計(jì)算要求,使得工程多尺度計(jì)算存在難以插值、強(qiáng)耦合及收斂性差等缺點(diǎn)[14]。第二代小波理論的提出有效地解決了這些問(wèn)題,它不再依賴于伸縮、平移變換來(lái)構(gòu)造小波,而是利用提升系數(shù)靈活地構(gòu)造具有期望特性的小波基,例如緊支撐、對(duì)稱、消失矩等特性,給工程問(wèn)題提供了非常豐富的小波基[15-16]。目前,第二代小波基在逐步應(yīng)用于工程數(shù)值計(jì)算方面。2004年,PINHO等應(yīng)用第二代小波基求解Maxwell方程,獲得較高的求解效率[17]。2006年,WANG等利用第二代小波構(gòu)造動(dòng)態(tài)計(jì)算網(wǎng)格求解波動(dòng)方程[18]。2007年,HE等提出基于預(yù)測(cè)算子和更新算子的第二代小波求解算法[19]。2009年,MEHRA等提出自適應(yīng)多尺度第二代小波方法,有效求解了橢球面方程[20]。

    如何根據(jù)對(duì)流擴(kuò)散問(wèn)題的需要靈活、有效構(gòu)造具有期望特性的第二代小波基成為擺在數(shù)學(xué)家、工程研究人員面前的難題。為解決這些問(wèn)題,本文提出基于提升格式的解耦小波構(gòu)造方法,根據(jù)問(wèn)題來(lái)構(gòu)造解耦小波基,以求達(dá)到高效求解問(wèn)題的目的。

    1 第二代小波

    通過(guò)提升方案構(gòu)造的小波即為第二代小波,第二代小波多分辨分析如下。

    L2(R)空間上的多分辨分析是一族閉子空間序列R={Vj?L2(R)|j∈Z},滿足下述條件1)—條件3)[11]:

    1)Vj?Vj+1;

    2)所有Vj的并集在L2(R)中是稠密的,即

    3)對(duì)于每一個(gè)j∈Z,Vj具有一個(gè)由尺度函數(shù){φj,k|k∈K(j)}給定的Riesz基,K(j)為j尺度k取值的集合。設(shè)Wj是Vj在Vj+1中的互補(bǔ)空間,Wj由小波函數(shù)ψj,m(x)張成,對(duì)于每個(gè)點(diǎn)m∈M(j),有M(j)=K(j+1)K(j),即M(j)是差集。

    (1)

    (2)

    式中hj,k,l和gj,m,l分別為低通濾波器系數(shù)和高通濾波器系數(shù)。根據(jù)文獻(xiàn)[21]定理:如果L2(R)有限元空間中的插值函數(shù)集{φj,k}k∈K(j)是稠密完備的,那么有限元空間V可構(gòu)成一個(gè)多分辨空間Vj。令有限元插值函數(shù)為小波的尺度函數(shù),那么對(duì)應(yīng)的小波函數(shù)就是尺度函數(shù)的補(bǔ)空間中的細(xì)節(jié)插值函數(shù)。有限元多分辨空間中的尺度函數(shù)滿足式(3)的細(xì)化關(guān)系:

    圖1 二次插值函數(shù)的兩尺度關(guān)系Fig.1 Two-scale relation for the quadratic interpolation function

    (3)

    其中g(shù)j,k,m為小波的細(xì)化系數(shù),是通過(guò)聯(lián)立相鄰兩尺度上的尺度函數(shù)和小波函數(shù)節(jié)點(diǎn)方程組來(lái)求解的。第二代小波多分辨空間中小波函數(shù)和尺度函數(shù)滿足關(guān)系ψj,m=φj+1,m。圖1表示的是單元端點(diǎn)尺度函數(shù)和小波函數(shù)及其兩尺度關(guān)系。

    二次Lagrange多分辨有限元空間的二尺度關(guān)系式為

    (4)

    2 對(duì)流擴(kuò)散方程的多尺度算法

    2.1多尺度計(jì)算框架

    對(duì)流-擴(kuò)散方程的雙線性弱解形式為

    a(u,v)=f(v), ?v∈V,

    (5)

    利用小波數(shù)值求解方法對(duì)解空間V構(gòu)造任意平方可積實(shí)數(shù)空間L2(R)上的多分辨分析{Vj,j∈Z}, 對(duì)式(5)中的試探函數(shù)u和容許函數(shù)v分別選取為小波函數(shù)和尺度函數(shù),那么可表示為式(6)形式:

    a(φj+1,ψj+1)=f(φj+1),?vj∈Vj;

    (6)

    若采用多分辨分析來(lái)逼近真解u(x),則j+1尺度的小波解可以由j尺度的尺度函數(shù)和小波函數(shù)表示如下:

    (7)

    對(duì)應(yīng)于式(7)的j尺度的小波離散方程組為

    Kj+1uj+1=fj+1,

    (8)

    其中

    (9)

    (10)

    Kjuj=fj,

    (11)

    Ke,jrj=pj,

    (12)

    Ke,j+1rj+1=pj+1,…,

    (13)

    具有解耦特性的提升小波可為解耦小波,其優(yōu)點(diǎn)在于系統(tǒng)剛度矩陣具有沿對(duì)角線的強(qiáng)稀疏性,實(shí)現(xiàn)工程問(wèn)題在每個(gè)尺度上的獨(dú)立、快速求解,使得系統(tǒng)方程的求解效率得到迅速提高。

    2.2解耦條件

    (14)

    式中sj,k,m為提升系數(shù)。提升方案提供了一種簡(jiǎn)單、靈活的工具,通過(guò)選擇合適的提升系數(shù)sj,k,m來(lái)獲得用戶定義的消失矩特性,可設(shè)計(jì)各種滿足問(wèn)題需要的解耦小波。解耦小波通過(guò)式(15)求解提升系數(shù)sj,k,m:

    (15)

    為保證多尺度剛度矩陣的解耦特性,必須根據(jù)小波的消失矩來(lái)保證尺度函數(shù)和小波函數(shù)之間的算子正交性。在多分辨有限元空間中,如果小波在積分邊界上為零,那么小波的導(dǎo)數(shù)將會(huì)有一階消失矩。一個(gè)n階小波必須提升到n階消失矩來(lái)保證算子正交性。具體推導(dǎo)形式如下:

    (16)

    (17)

    2.3多尺度解耦小波算法

    多尺度解耦小波算法是通過(guò)提升尺度,即在全域添加解耦小波,逐步逼近問(wèn)題精確解。具體步驟如下:

    1) 給定誤差閾值γ,生成初始尺度空間Vj的有限元網(wǎng)格,并求解初始尺度的方程K0u0=f0,計(jì)算初始尺度的相對(duì)誤差η0;

    2)比較相對(duì)誤差ηj(j=0)和誤差閾值γ,如果ηj≤γ,停止計(jì)算,輸出結(jié)果;

    3) 提高逼近空間的尺度,j=j+1,并利用提升方法構(gòu)造j尺度的解耦小波;

    4) 添加小波空間Wj+1的解耦小波到原始尺度的多尺度剛陣,回到步驟2)。

    (18)

    3 算 例

    設(shè)對(duì)流擴(kuò)散方程為-εu″+bu′=f,設(shè)b=1,ε=0.01 ,并給定邊界條件

    u(0)=0,u(1)=1 , 右端項(xiàng)為f(x)=-1,

    其解析解為

    該問(wèn)題在靠近邊界1處存在奇異,通過(guò)式(8)可推導(dǎo)該問(wèn)題的耦合矩陣為

    (19)

    分別采用二次和三次Lagrange解耦小波求解此對(duì)流擴(kuò)散方程,構(gòu)造解耦小波如圖2所示。給定閾值γ=0.02%,采用二次及三次Lagrange解耦小波多尺度算法求解該對(duì)流擴(kuò)散方程問(wèn)題。圖3表示的是對(duì)應(yīng)于三次Lagrange解耦小波的多尺度剛度矩陣(nz表示非零元素的個(gè)數(shù)),該剛陣是沿對(duì)角高度稀疏的完全解耦矩陣,適合高效的多尺度計(jì)算。為了分析該對(duì)流擴(kuò)散方程多尺度求解的相對(duì)誤差,給出了解耦小波解的細(xì)節(jié)圖及相對(duì)誤差圖。圖4和圖5表示的是多尺度解耦小波解的相對(duì)誤差與尺度和自由度之間的關(guān)系曲線圖。表1列出的是二次和三次解耦小波解。表2列出的是二次和三次解耦小波算法的求解時(shí)間比較。由圖可知,該問(wèn)題的三次解耦小波解比二次解耦小波解具有更快的收斂速度。因此,對(duì)流擴(kuò)散方程的解耦小波解法具有較高的求解精度。從表2分析可得出,多尺度解耦小波算法的計(jì)算時(shí)間較傳統(tǒng)自適應(yīng)有限元法少,其原因是利用解耦小波所構(gòu)造的多尺度剛度矩陣為沿對(duì)角的稀疏矩陣,提高了剛度矩陣的稀疏性,避免了重新劃分網(wǎng)格或提高插值函數(shù)的階次所帶來(lái)的大量有限元前處理工作和時(shí)間,從而提高了對(duì)流擴(kuò)散問(wèn)題的求解效率。

    圖2 解耦小波 Fig.2 Decoupling wavelets

    圖3 剛度矩陣稀疏度Fig.3 Sparsity of stiffness matrix

    圖4 解耦小波解的誤差-尺度Fig.4 Estimation of decoupled wavelet solution at each scale

    圖5 解耦小波解的誤差-自由度Fig.5 Estimation of decoupled wavelet solution with degrees of freedom

    表1 對(duì)流擴(kuò)散方程的多尺度解耦小波解

    表2 各種自適應(yīng)方法的運(yùn)行時(shí)間

    4 結(jié) 語(yǔ)

    提升方案的優(yōu)點(diǎn)在于根據(jù)對(duì)流擴(kuò)散問(wèn)題的需要靈活設(shè)計(jì)各種具有緊支撐、高消失矩的解耦小波。本文提出的方法可根據(jù)對(duì)流擴(kuò)散方程求解的需要構(gòu)造不同階次的解耦小波。解耦小波的解耦特性使得系統(tǒng)剛度矩陣具有沿對(duì)角線的強(qiáng)稀疏性,實(shí)現(xiàn)解耦小波在每個(gè)尺度上獨(dú)立求解,使得系統(tǒng)方程的求解效率得到迅速提高。與傳統(tǒng)有限元相比,多尺度解耦小波算法避免了重新劃分網(wǎng)格或提高插值函數(shù)的階次所帶來(lái)的大量有限元前處理工作和時(shí)間,提高了對(duì)流擴(kuò)散問(wèn)題的分析效率。計(jì)算結(jié)果表明,通過(guò)提高解耦小波的階次或尺度,可實(shí)現(xiàn)對(duì)流擴(kuò)散方程的高效求解。

    /

    [1] 張引娣,封建湖,余劍生. 二維非線性對(duì)流擴(kuò)散方程特征有限元方法[J]. 長(zhǎng)安大學(xué)學(xué)報(bào),2010,30(3): 105-110.

    ZHANG Yindi,FENG Jianhu,YU Jiansheng. Characteristic finite element method for 2-D nonlinear convection-diffusion problems[J]. Journal of Chang′an University,2010,30(3): 105-110.

    [2] 蔣建國(guó), 吳吉春. 格子行走方法求解一維對(duì)流擴(kuò)散方程[J]. 南京大學(xué)學(xué)報(bào), 2011, 47(3): 276-280.

    JIANG Jianguo,WU Jichun. Lattice-walk method for the convection-diffusion equation[J].Journal of Nanjing University, 2011, 47(3): 276-280.

    [3] 紀(jì)冬梅, 軒福貞, 涂善東. 矩形板屈曲問(wèn)題的一個(gè)小波解[J]. 應(yīng)用力學(xué)學(xué)報(bào),2011,28(1):95-99.

    JI Dongmei, XUAN Fuzhen, TU Shandong.A wavelet solution of high order buckling problem for clamped rectangular plates[J].Chinese Journal of Applied Mechanics, 2011,28(1):95-99.

    [4] 王保國(guó), 吳俊宏, 朱俊強(qiáng). 基于小波奇異分析的流場(chǎng)計(jì)算方法及應(yīng)用[J].航空動(dòng)力學(xué)報(bào), 2010,25(12): 2 747-2 748.

    WANG Baoguo,WU Junhong,ZHU Junqiang.Method based on wavelet sigularity analysis for complicated flow and its application[J].Journal of Aerospace Power, 2010,25(12): 2 747-2 748.

    [5] 張 悅, 杜守軍, 張麗梅.小波奇異性在鋼結(jié)構(gòu)損傷檢測(cè)中的應(yīng)用[J].河北科技大學(xué)學(xué)報(bào),2010,31(2):151-157.

    ZHANG Yue,DU Shoujun,ZHANG Limei.Application of wavelet singularity to steel structural damage detection[J].Journal of Hebei University of Science and Technology,2010,31(2):151-157.

    [6] 陳軍霞,孫幸成, 基于小波變換的神經(jīng)網(wǎng)絡(luò)容錯(cuò)控制研究[J].河北科技大學(xué)學(xué)報(bào),2011,32(6): 567-570.

    CHEN Junxia,SUN Xingcheng. Study on fault-tolerant controller based on wavelet transform compensatory network[J].Journal of Hebei University of Science and Technology,2011,32(6): 567-570.

    [7] HASHISH H, BEHIRY S H,ELSAID A. Solving the 2-D heat equations using wavelet-Galerkin method with variable time step [J]. Applied Mathematics and Computation, 2009, 213(1): 209-215.

    [8] DOU Fangfang, FU Chuli, MA Yunjie. A wavelet-Galerkin method for high order numerical differentiation[J]. Applied Mathematics and Computation, 2010, 215(10): 3 702-3 712.

    [9] BERTOLUZZA S, NALDI G. A wavelet collocation method for the numerical solution of partial differential equations[J]. Applied and Computational Harmonic Analysis, 1996, 3(1): 1-9.

    [10] KANEKO H, NOREN R D, NOVAPRATEEP R. Wavelet applications to the Petrov-Galerkin method for Hammerstein equations [J]. Applied Numerical Mathematics, 2003, 45(2/3): 255-273.

    [11] DAHMEN W, KUNOTH A, SCHNEIDER R. Wavelet least squares methods for boundary value problems[J]. SIAM Journal on Numerical Analysis, 2002, 39(6): 1 985-2 013.

    [12] DONG Hongbo, CHEN Xuefeng, LI Bing, et al. Rotor crack detection based on high-precision modal parameter identification method and wavelet finite element model [J]. Mechanical Systems and Signal Processing, 2009, 23(3): 869-883.

    [13] CHEN Xuefeng, XIANG Jiawei, LI Bing, et al. A study of multiscale wavelet-based elements for adaptive finite element analysis[J]. Advances in Engineering Software, 2010, 41(2): 196-205.

    [14] 何正嘉,陳雪峰,李 兵, 等. 小波有限元理論及其工程應(yīng)用[M]. 北京:科學(xué)出版社,2006.

    HE Zhengjia, CHEN Xuefeng, LI Bing, et al. Wavelet Finite Element Theory and Its Engineering Applications [M].Beijing: Science Press, 2006.

    [15] SWELDENS W. The lifting scheme:A construction of second generation wavelets[J]. SIAM Journal of Mathematical Analysis, 1997, 29 (2): 511-546.

    [16] SWELDENS W. The lifting scheme: A custom-design construction of biorthogonal wavelets [J]. Applied and Computational Harmonic Analysis, 1996, 3 (2): 186-200.

    [17] PINHO P, FERREIRA P, PEREIRA J R. Multiresolution analysis using biorthogonal and interpolating wavelets [J]. IEEE Antennas and Propagation Society Symposium, 2004, 2: 1 483-1 486.

    [18] WANG Yibo, YANG Huizhu. Second generation wavelet based on adaptive solution of wave equation[J]. International Journal of Nonlinear Sciences and Numerical Simulation, 2006,7(4): 435-438.

    [19] HE Yumin, CHEN Xuefeng, XIANG Jiawei, et al. Adaptive multiresolution finite element method based on second generation wavelets[J]. Finite Elements in Analysis and Design, 2007, 43: 566-579.

    [20] MEHRA M, KEVLAHAN N. An adaptive multilevel wavelet solver for elliptic equations on an optimal spherical geodesic grid [J]. SIAM Journal on Scientific Computing, 2009, 30(6): 3 073-3 086.

    [21] BATHE K J. Finite Element Procedures [M]. Prentice-Hall: Upper Saddle River, 1996.

    A multi-scale decoupling wavelet algorithm for convection-diffusion equations

    WANG Youming1,2, WU Qing1, SHEN Jiandong1

    (1. School of Automation, Xi′an University of Posts and Telecommunications, Xi′an Shaanxi 710121, China;2. State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian Liaoning 116023, China)

    In order to solve the convection-diffusion problems of numerical oscillation and dissipation in traditional finite element method, a scaled decoupled wavelet solution method is proposed. Firstly, the multi-resolution analysis of second generation wavelets is introduced, and the finite element two multi-resolution spacial scaling relation is derived, then the multi-scale computational framework is presented for solving convection-diffusion equations. The decoupling conditions of convection-diffusion equations are developed, and the multi-scale decoupling wavelets are constructed by the lifting scheme. A multi-level decoupled wavelets algorithm is proposed for approximating the exact solution by adding the decoupled wavelets into the solving domain. Numerical example shows that the decoupled wavelets have good computational performance in solving convection-diffusion equations.

    convection-diffusion equations; multi-scale computation; decoupled wavelets.

    1008-1542(2013)04-0269-06

    10.7535/hbkd.2013yx04023

    O241.82;O242.21

    A

    2013-03-24;

    2013-05-28;責(zé)任編輯:李 穆

    國(guó)家自然科學(xué)基金(51205309, 61100165);工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室開放課題基金(GZ1209); 陜西省教育廳自然科學(xué)專項(xiàng)(2013JK0992)

    汪友明(1981-),男,湖北黃岡人,講師,博士,主要從事小波數(shù)值分析與工程應(yīng)用方面的研究。

    E-mail: xautroland@163.com

    猜你喜歡
    對(duì)流小波尺度
    齊口裂腹魚集群行為對(duì)流態(tài)的響應(yīng)
    構(gòu)造Daubechies小波的一些注記
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    基于改進(jìn)的G-SVS LMS 與冗余提升小波的滾動(dòng)軸承故障診斷
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    基于ANSYS的自然對(duì)流換熱系數(shù)計(jì)算方法研究
    二元驅(qū)油水界面Marangoni對(duì)流啟動(dòng)殘余油機(jī)理
    9
    基于FPGA小波變換核的設(shè)計(jì)
    男人舔女人的私密视频| 亚洲av国产av综合av卡| 国产成人91sexporn| 国产精品国产三级国产专区5o| 波多野结衣一区麻豆| 亚洲国产欧美日韩在线播放| 插逼视频在线观看| 黑人欧美特级aaaaaa片| 午夜福利乱码中文字幕| 人人妻人人澡人人看| 国产在线一区二区三区精| 日韩伦理黄色片| 精品少妇久久久久久888优播| 国产黄色视频一区二区在线观看| 免费看光身美女| 老司机亚洲免费影院| 亚洲av日韩在线播放| 美女中出高潮动态图| 一二三四中文在线观看免费高清| av在线观看视频网站免费| 岛国毛片在线播放| 欧美 日韩 精品 国产| 9色porny在线观看| 国产精品 国内视频| 赤兔流量卡办理| 国产免费一区二区三区四区乱码| 视频区图区小说| 九草在线视频观看| 在线观看免费日韩欧美大片| 精品国产乱码久久久久久小说| 日韩一区二区视频免费看| 一区二区三区四区激情视频| 欧美日韩精品成人综合77777| 国产精品99久久99久久久不卡 | 国产精品 国内视频| 咕卡用的链子| 啦啦啦视频在线资源免费观看| 人妻 亚洲 视频| 全区人妻精品视频| 黄色毛片三级朝国网站| 最近最新中文字幕大全免费视频 | 少妇高潮的动态图| 久久久久国产网址| 欧美成人午夜免费资源| 内地一区二区视频在线| 成人国语在线视频| 精品第一国产精品| 人妻 亚洲 视频| 91在线精品国自产拍蜜月| 亚洲经典国产精华液单| 永久网站在线| 成人影院久久| 香蕉国产在线看| 久久精品aⅴ一区二区三区四区 | 久久影院123| 在线天堂中文资源库| 午夜福利影视在线免费观看| 久久99蜜桃精品久久| 精品一区二区三卡| 日韩 亚洲 欧美在线| 香蕉丝袜av| 精品福利永久在线观看| 欧美国产精品一级二级三级| 在线天堂中文资源库| 2021少妇久久久久久久久久久| 日本vs欧美在线观看视频| 亚洲,欧美,日韩| 久久精品夜色国产| 日韩欧美精品免费久久| 99热全是精品| 亚洲激情五月婷婷啪啪| 国产一区有黄有色的免费视频| 日日爽夜夜爽网站| 成人黄色视频免费在线看| 狂野欧美激情性bbbbbb| 美女国产视频在线观看| 亚洲人成77777在线视频| 在线观看www视频免费| 成人漫画全彩无遮挡| 亚洲成人手机| 91午夜精品亚洲一区二区三区| 国产爽快片一区二区三区| 久久久国产精品麻豆| 2022亚洲国产成人精品| 涩涩av久久男人的天堂| 在线观看美女被高潮喷水网站| 90打野战视频偷拍视频| 亚洲综合精品二区| 26uuu在线亚洲综合色| 看十八女毛片水多多多| 三级国产精品片| 日本vs欧美在线观看视频| 亚洲性久久影院| 最近的中文字幕免费完整| 国产男女内射视频| 男人添女人高潮全过程视频| 两性夫妻黄色片 | 男人舔女人的私密视频| 日韩制服骚丝袜av| 亚洲人成77777在线视频| 尾随美女入室| 成年女人在线观看亚洲视频| 极品少妇高潮喷水抽搐| a级毛色黄片| 国产免费视频播放在线视频| av福利片在线| 日韩一本色道免费dvd| 国产精品久久久久久久电影| 国产色婷婷99| 日本av手机在线免费观看| 少妇人妻久久综合中文| 美女大奶头黄色视频| 国产精品久久久久久av不卡| 欧美日韩一区二区视频在线观看视频在线| 精品国产一区二区久久| 精品人妻在线不人妻| a级毛片黄视频| av又黄又爽大尺度在线免费看| 精品人妻偷拍中文字幕| 18+在线观看网站| 精品少妇内射三级| 欧美成人午夜精品| av天堂久久9| 亚洲精品自拍成人| 男人舔女人的私密视频| 天美传媒精品一区二区| 亚洲第一av免费看| 欧美老熟妇乱子伦牲交| 色婷婷av一区二区三区视频| 精品酒店卫生间| 欧美 日韩 精品 国产| 国产精品免费大片| 成人午夜精彩视频在线观看| 久久青草综合色| 美国免费a级毛片| 国产在线一区二区三区精| 2018国产大陆天天弄谢| 这个男人来自地球电影免费观看 | 高清不卡的av网站| 最近的中文字幕免费完整| 国产欧美日韩综合在线一区二区| 热99国产精品久久久久久7| 久久99蜜桃精品久久| 日韩电影二区| 最新的欧美精品一区二区| 国产精品久久久久久久久免| 精品少妇久久久久久888优播| 精品国产露脸久久av麻豆| 国产成人精品无人区| 夫妻午夜视频| 国产又爽黄色视频| 午夜激情久久久久久久| 午夜福利,免费看| 国产免费一区二区三区四区乱码| 成人免费观看视频高清| 最新的欧美精品一区二区| 91精品伊人久久大香线蕉| 国产又色又爽无遮挡免| 精品国产国语对白av| xxx大片免费视频| 国产av精品麻豆| 精品第一国产精品| 国产精品国产三级国产av玫瑰| 看十八女毛片水多多多| 看免费av毛片| a级毛片在线看网站| 伦精品一区二区三区| 一级毛片我不卡| 亚洲精品久久久久久婷婷小说| 一本久久精品| 国产精品三级大全| 国产色婷婷99| 免费观看av网站的网址| 青春草视频在线免费观看| 美女国产高潮福利片在线看| 午夜激情久久久久久久| 黄色一级大片看看| 久久精品夜色国产| 亚洲国产精品成人久久小说| 国产男女内射视频| 免费女性裸体啪啪无遮挡网站| 国产69精品久久久久777片| 国产精品久久久久久久久免| 一区二区三区精品91| 日本欧美视频一区| 国产精品久久久久久精品古装| 久久热在线av| 欧美亚洲 丝袜 人妻 在线| 亚洲成人手机| 久久99蜜桃精品久久| 波野结衣二区三区在线| 人妻 亚洲 视频| 国产精品久久久久久久久免| 51国产日韩欧美| 一边摸一边做爽爽视频免费| 熟女电影av网| 亚洲熟女精品中文字幕| 晚上一个人看的免费电影| 草草在线视频免费看| 69精品国产乱码久久久| 中文字幕制服av| 精品少妇久久久久久888优播| 青春草亚洲视频在线观看| 男女下面插进去视频免费观看 | 精品亚洲乱码少妇综合久久| 波多野结衣一区麻豆| 精品少妇黑人巨大在线播放| 亚洲精品第二区| 国国产精品蜜臀av免费| 精品国产一区二区三区四区第35| 日韩免费高清中文字幕av| 久久青草综合色| 久久久精品94久久精品| 国产精品一二三区在线看| 大香蕉97超碰在线| 日本欧美国产在线视频| 日韩成人av中文字幕在线观看| 青春草视频在线免费观看| 丝袜脚勾引网站| 国产成人免费无遮挡视频| 国产国拍精品亚洲av在线观看| 永久免费av网站大全| 一区二区日韩欧美中文字幕 | 丰满少妇做爰视频| 国产激情久久老熟女| 香蕉精品网在线| 两个人免费观看高清视频| 一本—道久久a久久精品蜜桃钙片| 国产成人一区二区在线| 蜜桃在线观看..| 国产精品一二三区在线看| 中国三级夫妇交换| 免费女性裸体啪啪无遮挡网站| 久久精品国产亚洲av天美| 免费黄色在线免费观看| 国产免费视频播放在线视频| 在线观看美女被高潮喷水网站| 多毛熟女@视频| 中国三级夫妇交换| 伦精品一区二区三区| 免费少妇av软件| 亚洲国产精品999| 九九在线视频观看精品| 男女免费视频国产| 大话2 男鬼变身卡| 久久久久久久久久人人人人人人| 国产xxxxx性猛交| 新久久久久国产一级毛片| 国产亚洲精品第一综合不卡 | 人成视频在线观看免费观看| 蜜桃在线观看..| 曰老女人黄片| 免费av不卡在线播放| 男女午夜视频在线观看 | 黄片播放在线免费| 国产麻豆69| 69精品国产乱码久久久| 亚洲天堂av无毛| 亚洲精品美女久久久久99蜜臀 | 亚洲欧洲日产国产| 国产免费又黄又爽又色| 国产69精品久久久久777片| 90打野战视频偷拍视频| 国产老妇伦熟女老妇高清| 一区在线观看完整版| 色网站视频免费| 欧美日韩综合久久久久久| 蜜桃在线观看..| 欧美人与善性xxx| 亚洲国产日韩一区二区| 伊人亚洲综合成人网| 美女国产视频在线观看| 精品第一国产精品| 日韩制服骚丝袜av| 国产老妇伦熟女老妇高清| 男女下面插进去视频免费观看 | 黑人猛操日本美女一级片| 久久久久视频综合| 久久久久精品久久久久真实原创| 黑人高潮一二区| 国产男女内射视频| 午夜av观看不卡| 一区二区av电影网| 美女国产视频在线观看| 久久久国产一区二区| 国产日韩一区二区三区精品不卡| 国产国拍精品亚洲av在线观看| 一级a做视频免费观看| 国产免费又黄又爽又色| 亚洲欧美成人精品一区二区| 深夜精品福利| 亚洲国产日韩一区二区| 美女脱内裤让男人舔精品视频| 成年美女黄网站色视频大全免费| 人体艺术视频欧美日本| 天天操日日干夜夜撸| 亚洲色图 男人天堂 中文字幕 | 中文天堂在线官网| 在线观看免费日韩欧美大片| 草草在线视频免费看| 亚洲综合色网址| 国产亚洲欧美精品永久| 青青草视频在线视频观看| 久久 成人 亚洲| 侵犯人妻中文字幕一二三四区| 美女视频免费永久观看网站| 精品亚洲成a人片在线观看| 高清视频免费观看一区二区| 国产欧美亚洲国产| 永久免费av网站大全| 免费观看a级毛片全部| 中文欧美无线码| 日本午夜av视频| 草草在线视频免费看| 嫩草影院入口| 女的被弄到高潮叫床怎么办| 午夜影院在线不卡| 久久女婷五月综合色啪小说| 日日撸夜夜添| 亚洲精品av麻豆狂野| 麻豆精品久久久久久蜜桃| 午夜免费观看性视频| 欧美精品人与动牲交sv欧美| 日本猛色少妇xxxxx猛交久久| 美女内射精品一级片tv| 超碰97精品在线观看| 午夜影院在线不卡| 成年美女黄网站色视频大全免费| 一本色道久久久久久精品综合| 天美传媒精品一区二区| 五月天丁香电影| 欧美日韩亚洲高清精品| 在线观看人妻少妇| 久久精品aⅴ一区二区三区四区 | 亚洲国产av影院在线观看| 国产精品一区二区在线不卡| 9191精品国产免费久久| 国产精品久久久久久精品古装| 亚洲丝袜综合中文字幕| 精品视频人人做人人爽| 中文字幕另类日韩欧美亚洲嫩草| 人人妻人人爽人人添夜夜欢视频| 一级片免费观看大全| 69精品国产乱码久久久| 亚洲婷婷狠狠爱综合网| 国产精品欧美亚洲77777| 日韩成人av中文字幕在线观看| 天堂8中文在线网| 日韩欧美精品免费久久| 国产亚洲精品久久久com| 欧美亚洲日本最大视频资源| 久久国产精品男人的天堂亚洲 | 亚洲欧美精品自产自拍| 亚洲av成人精品一二三区| 免费av中文字幕在线| 咕卡用的链子| 亚洲精品456在线播放app| 日韩视频在线欧美| 天天影视国产精品| 久久久久久久久久久免费av| av网站免费在线观看视频| 各种免费的搞黄视频| 高清不卡的av网站| 精品卡一卡二卡四卡免费| 熟妇人妻不卡中文字幕| 久久精品国产亚洲av天美| 三上悠亚av全集在线观看| 99久久人妻综合| av有码第一页| 久久精品国产自在天天线| 女人精品久久久久毛片| 久久热在线av| 久久久久久人妻| 日韩熟女老妇一区二区性免费视频| 国产精品成人在线| 午夜久久久在线观看| 丝袜脚勾引网站| 又黄又爽又刺激的免费视频.| 人体艺术视频欧美日本| 亚洲欧美成人综合另类久久久| 日日摸夜夜添夜夜爱| 最近2019中文字幕mv第一页| 中文字幕免费在线视频6| 国产成人午夜福利电影在线观看| 男女下面插进去视频免费观看 | 日韩大片免费观看网站| 最近中文字幕高清免费大全6| 免费看av在线观看网站| 捣出白浆h1v1| 日韩中字成人| 中文字幕最新亚洲高清| 国产成人aa在线观看| 国产亚洲av片在线观看秒播厂| 满18在线观看网站| 欧美国产精品一级二级三级| 欧美 日韩 精品 国产| 日本免费在线观看一区| 国产一区亚洲一区在线观看| 黄网站色视频无遮挡免费观看| 国产高清不卡午夜福利| 一边亲一边摸免费视频| 婷婷成人精品国产| 亚洲成人av在线免费| 欧美国产精品一级二级三级| 97在线人人人人妻| 午夜日本视频在线| 久久综合国产亚洲精品| 天天躁夜夜躁狠狠躁躁| 亚洲国产精品999| 狠狠婷婷综合久久久久久88av| 大片电影免费在线观看免费| 黄色视频在线播放观看不卡| 免费人成在线观看视频色| 肉色欧美久久久久久久蜜桃| 极品人妻少妇av视频| 搡老乐熟女国产| 三级国产精品片| 人人妻人人爽人人添夜夜欢视频| 99久久中文字幕三级久久日本| videossex国产| 欧美成人午夜免费资源| 午夜福利,免费看| 波多野结衣一区麻豆| 青春草亚洲视频在线观看| 亚洲精品第二区| 亚洲精品久久成人aⅴ小说| 国产高清不卡午夜福利| 亚洲综合色惰| 国产亚洲精品久久久com| 乱码一卡2卡4卡精品| 中文字幕另类日韩欧美亚洲嫩草| 十分钟在线观看高清视频www| 天美传媒精品一区二区| 51国产日韩欧美| 日韩精品免费视频一区二区三区 | 天天躁夜夜躁狠狠久久av| 国产一区二区在线观看av| av电影中文网址| 久久人人97超碰香蕉20202| 免费女性裸体啪啪无遮挡网站| 色吧在线观看| 日本欧美国产在线视频| 国产又色又爽无遮挡免| 亚洲国产精品一区二区三区在线| 国产麻豆69| av电影中文网址| 成人综合一区亚洲| av在线老鸭窝| 免费在线观看完整版高清| 久久99蜜桃精品久久| 久久影院123| 我要看黄色一级片免费的| 久久国产精品男人的天堂亚洲 | 国产xxxxx性猛交| 国产福利在线免费观看视频| 中国三级夫妇交换| 国产高清国产精品国产三级| 国产精品嫩草影院av在线观看| 精品视频人人做人人爽| 国产黄频视频在线观看| 搡老乐熟女国产| 亚洲精品国产av蜜桃| 午夜福利,免费看| 亚洲精品乱码久久久久久按摩| 制服诱惑二区| 久久久精品94久久精品| 久久影院123| 国产精品国产三级国产专区5o| 久久精品aⅴ一区二区三区四区 | 成年人免费黄色播放视频| 欧美最新免费一区二区三区| 亚洲av国产av综合av卡| 欧美精品av麻豆av| 日韩制服丝袜自拍偷拍| 制服丝袜香蕉在线| 在线观看美女被高潮喷水网站| 精品亚洲成a人片在线观看| 亚洲精品视频女| 久久青草综合色| 国产亚洲av片在线观看秒播厂| 亚洲av欧美aⅴ国产| 亚洲美女视频黄频| 水蜜桃什么品种好| 欧美老熟妇乱子伦牲交| 大香蕉久久成人网| 精品一区二区免费观看| 国产男女超爽视频在线观看| 丰满乱子伦码专区| 久久精品国产a三级三级三级| 婷婷色麻豆天堂久久| 亚洲综合精品二区| 亚洲在久久综合| 一级a做视频免费观看| 国产爽快片一区二区三区| 欧美精品一区二区大全| 久久韩国三级中文字幕| 精品一品国产午夜福利视频| 五月伊人婷婷丁香| av女优亚洲男人天堂| 波野结衣二区三区在线| 黄色怎么调成土黄色| 国产日韩欧美亚洲二区| 中文欧美无线码| 亚洲国产精品成人久久小说| 亚洲中文av在线| 精品久久久久久电影网| 国产精品国产三级国产专区5o| 国产成人精品久久久久久| 啦啦啦中文免费视频观看日本| 精品一区二区三卡| 欧美人与性动交α欧美精品济南到 | 久久久久视频综合| 性色av一级| 交换朋友夫妻互换小说| 嫩草影院入口| videos熟女内射| 人体艺术视频欧美日本| 男女国产视频网站| 蜜桃国产av成人99| 色94色欧美一区二区| 啦啦啦在线观看免费高清www| 91在线精品国自产拍蜜月| 又黄又爽又刺激的免费视频.| 香蕉国产在线看| 久久久久网色| 精品国产露脸久久av麻豆| 亚洲成人手机| 国产片特级美女逼逼视频| 欧美xxⅹ黑人| 国精品久久久久久国模美| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产精品成人久久小说| 69精品国产乱码久久久| 高清在线视频一区二区三区| 99久久中文字幕三级久久日本| 精品一品国产午夜福利视频| 免费观看性生交大片5| 久久国产精品大桥未久av| 最近手机中文字幕大全| 亚洲精品视频女| 亚洲国产成人一精品久久久| 日本wwww免费看| av线在线观看网站| 国产精品三级大全| 亚洲精品美女久久av网站| 欧美性感艳星| 免费大片黄手机在线观看| 一本大道久久a久久精品| 蜜臀久久99精品久久宅男| 狠狠精品人妻久久久久久综合| 国产亚洲欧美精品永久| 欧美日本中文国产一区发布| 国产老妇伦熟女老妇高清| 哪个播放器可以免费观看大片| 在线观看一区二区三区激情| 免费大片18禁| 人体艺术视频欧美日本| 国产精品一二三区在线看| 久久精品久久久久久噜噜老黄| 日韩制服丝袜自拍偷拍| 国产亚洲欧美精品永久| 捣出白浆h1v1| 一本久久精品| a级毛片在线看网站| 亚洲成人一二三区av| 中国三级夫妇交换| 熟女av电影| 亚洲四区av| 久久99精品国语久久久| 国产又爽黄色视频| 久久久久久久精品精品| 人成视频在线观看免费观看| 中文字幕制服av| 男女国产视频网站| 免费看光身美女| 国产精品熟女久久久久浪| 亚洲美女黄色视频免费看| 天天影视国产精品| 午夜免费鲁丝| 欧美日韩一区二区视频在线观看视频在线| 热re99久久精品国产66热6| 中国三级夫妇交换| 十八禁网站网址无遮挡| 高清欧美精品videossex| 亚洲色图 男人天堂 中文字幕 | 亚洲精品美女久久av网站| 精品国产一区二区三区四区第35| 多毛熟女@视频| 视频中文字幕在线观看| 中国美白少妇内射xxxbb| 久久影院123| a级片在线免费高清观看视频| 成人漫画全彩无遮挡| 欧美+日韩+精品| 久久青草综合色| 黄色配什么色好看| 色婷婷av一区二区三区视频| 美女福利国产在线| 亚洲成av片中文字幕在线观看 | videosex国产| 久久热在线av| 日韩视频在线欧美| 日韩三级伦理在线观看| av黄色大香蕉| 精品国产一区二区久久| 亚洲精品国产av蜜桃| 成人国产麻豆网| 国产无遮挡羞羞视频在线观看| 久久免费观看电影| 亚洲精品日本国产第一区| 日本免费在线观看一区| 春色校园在线视频观看| 一级毛片黄色毛片免费观看视频| 午夜老司机福利剧场| 亚洲精品美女久久av网站|