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

    基于位移的格子法模擬彈性波傳播與應(yīng)用*

    2010-12-01 03:58:18高廣運(yùn)馮世進(jìn)歐湘萍
    關(guān)鍵詞:差分法衰減系數(shù)邊界條件

    熊 浩 高廣運(yùn) 馮世進(jìn) 歐湘萍

    (臺(tái)州學(xué)院建筑工程學(xué)院1) 臺(tái)州 318000) (同濟(jì)大學(xué)地下建筑與工程系2) 上海 200092)

    (同濟(jì)大學(xué)巖土及地下工程教育部重點(diǎn)實(shí)驗(yàn)室3) 上海 200092) (武漢理工大學(xué)交通學(xué)院4) 武漢 430063)

    0 引 言

    彈性波傳播的數(shù)值模擬在鐵道及城市軌道交通環(huán)境振動(dòng)分析和隔振設(shè)計(jì)、地震勘探等領(lǐng)域有廣泛的應(yīng)用.有限差分法、有限元法和偽譜法是三種主要的數(shù)值模擬方法.Virieux等[1]提出了基于應(yīng)力-速度混合變量彈性波方程的交錯(cuò)網(wǎng)格差分法,提高了差分法的計(jì)算精度和穩(wěn)定性.裴正林[2]將交錯(cuò)網(wǎng)格計(jì)算技術(shù)與高精度差分計(jì)算格式結(jié)合起來,用于三維波動(dòng)模擬.Chen[3]則進(jìn)一步提出了旋轉(zhuǎn)交錯(cuò)網(wǎng)格計(jì)算技術(shù),增加了差分法模擬非均勻介質(zhì)中波動(dòng)問題的能力.有限元法雖理論成熟,適應(yīng)性強(qiáng),但計(jì)算所需存儲(chǔ)空間較大,計(jì)算效率稍差.偽譜法[4]則采用快速傅立葉變換計(jì)算空間導(dǎo)數(shù),用差分法計(jì)算時(shí)間導(dǎo)數(shù).通過引入映射函數(shù)雖可以使其在計(jì)算中考慮地表形狀的變化,但如何合理選擇映射函數(shù)有一定困難.

    近年來,張劍鋒[5-6]給出了基于應(yīng)力-速度混合變量彈性波方程及任意四邊形網(wǎng)格差分算子的數(shù)值方法——格子法,它融合了有限元可以適應(yīng)復(fù)雜形狀邊界條件、有限差分法計(jì)算簡(jiǎn)便及無須計(jì)算剛度矩陣的優(yōu)點(diǎn),可以方便地引入任意形狀自由表面邊界條件.

    本文基于以位移為變量彈性波方程,給出了二維格子法的計(jì)算方法,推導(dǎo)適合該法的人工邊界條件,并通過經(jīng)典Lam b問題算例驗(yàn)證了該方法的正確性.與傳統(tǒng)的基于混合變量的格子法相比,該方法采用單點(diǎn)高斯積分法計(jì)算線積分,提高了積分的準(zhǔn)確性;文中所給出的人工邊界條件計(jì)算簡(jiǎn)便,透射外行波效果良好.應(yīng)用上述數(shù)值計(jì)算方法,本文分析了在不同填充溝深度條件下的隔振效果,即地表豎向振幅衰減系數(shù)變化.

    1 計(jì)算方法

    1.1 基本計(jì)算公式

    在笛卡爾坐標(biāo)系中,不計(jì)體力,均質(zhì)彈性介質(zhì)中彈性波傳播控制方程為

    式中:ρ為彈性介質(zhì)的密度,kg/m3;μ和λ為介質(zhì)的Lamé系數(shù),Pa;σx,σy,σz和 τxy為應(yīng)力張量,Pa;ux和uy為位移,m.

    將計(jì)算區(qū)域以三角形或四邊形網(wǎng)格進(jìn)行離散,局部網(wǎng)格如圖1所示.圖中的實(shí)線是離散時(shí)所生成網(wǎng)格線,實(shí)心點(diǎn)j,k,l,…是網(wǎng)格線頂點(diǎn);圖中的虛線是各段實(shí)線中點(diǎn)與網(wǎng)格中心點(diǎn)的連線,四邊形網(wǎng)格中空心點(diǎn)(如1,2,3,…)是各段虛線的中點(diǎn),而三角形網(wǎng)格中的空心點(diǎn)(如點(diǎn)7)是三角形的中心點(diǎn).如果在實(shí)心點(diǎn)處定義位移分量ui(i=x,y),那么對(duì)于四邊形網(wǎng)格而言,可以利用任意四邊形網(wǎng)格差分算子[6]計(jì)算位移在域內(nèi)任意點(diǎn)處的空間一階導(dǎo)數(shù)值.

    圖1 局部網(wǎng)格示意圖

    對(duì)于三角形網(wǎng)格,以圖1中三角形 fm j為例,其中心點(diǎn)處的位移空間一階導(dǎo)數(shù)計(jì)算公式為[7]

    式中:A為三角形 fm j的面積;(ui)f,(ui)m和(ui)j為三個(gè)頂點(diǎn)f,m,j處的位移值;bf=(yjym)/2,cf=(xj-xm)/2,(f→m→j→t),通過下標(biāo)輪換可以類似地計(jì)算出 bm,cm等;(xj,yj),(xm,ym)分別為j,m點(diǎn)處的坐標(biāo).

    從圖1可見,各個(gè)實(shí)心點(diǎn)周圍的虛線都圍成一個(gè)閉合區(qū)域.假定結(jié)點(diǎn) f周圍虛線所圍的區(qū)域?yàn)?Ω,其邊界1-2-3-4-5-6-7記為 Γ.在 Ω上對(duì)式(1)積分可得

    式中:α,β分別為邊界 的外法線方向余弦.采用集中質(zhì)量模型,式(4)~(5)中的等號(hào)左端的面積積分可簡(jiǎn)化為僅與 f點(diǎn)加速度有關(guān),即分別為.式中:M f是區(qū)域 Ω的質(zhì)量.對(duì)于式(3)右端的線積分,在邊界Γ每一段上采用單點(diǎn)高斯積分公式進(jìn)行數(shù)值計(jì)算.此時(shí),每段上的積分點(diǎn)正是邊界 Γ上的各個(gè)空心點(diǎn).設(shè)式(4)和式(5)線積分值分別為(Ax)f和(Ay)f,則

    對(duì)式(6)中二階導(dǎo)數(shù)應(yīng)用中心差分公式計(jì)算,可得

    式中:上標(biāo) n為nΔt時(shí)刻,類似地n+1代表(n+1)Δt時(shí)刻,n-1為(n-1)Δt時(shí)刻;Δt為計(jì)算時(shí)間步長(zhǎng).根據(jù)式(7),可求得

    1.2 計(jì)算步驟

    對(duì)節(jié)點(diǎn) f而言,應(yīng)用單點(diǎn)高斯積分公式,在邊界Γ上計(jì)算式(3)中線積分(Ax)f和(Ay)f,并把它們代入到式(7)中,計(jì)算出(n+1)Δt時(shí)刻的位移(un+1x)f,(un+1y)f.這樣,就完成了 f點(diǎn)從nΔt時(shí)刻到(n+1)Δt時(shí)刻的位移分量遞推計(jì)算.對(duì)于其他實(shí)心點(diǎn),可以類似于 f點(diǎn)進(jìn)行處理.從而整個(gè)實(shí)線網(wǎng)格上的全部節(jié)點(diǎn)在(n+1)Δt時(shí)刻位移可以完全確定.類似地,可以計(jì)算出(n+1)Δt時(shí)刻的位移場(chǎng),直到設(shè)定的計(jì)算時(shí)間終點(diǎn)為止.

    為了保證計(jì)算的穩(wěn)定,時(shí)間步長(zhǎng)Δt與空間最小網(wǎng)格邊長(zhǎng)h應(yīng)滿足穩(wěn)定條件:Δt≤h/vp,式中:v p為縱波的波速.

    1.3 自由表面邊界條件

    假定在圖1中h,l,k是自由表面邊界上的節(jié)點(diǎn).在 l點(diǎn)處,在圍線 l-Ma-9-3-2-8-Mb所圍區(qū)域內(nèi)類似于式(4)~(5)建立平衡方程

    式中:M l為上述圍線所圍區(qū)域的質(zhì)量.

    式(9)中線積分分為兩部分,其中Ma-9-3-2-8-Mb線段上的積分計(jì)算方法與前述式(4)~(5)中的線積分計(jì)算辦法完全相同.而Mbl-Ma線段上的積分事實(shí)上就等于作用于表面邊界上的已知外荷載.因此,在自由邊界上邊界條件可以方便地引入到計(jì)算中.

    2 人工邊界條件

    考慮到實(shí)際上波的傳播是在半無限域中進(jìn)行的,而計(jì)算模型只能是有限的,因此必須在計(jì)算區(qū)域的邊界上引入人工邊界條件,以準(zhǔn)確地模擬波的無限傳播.本文給出一種簡(jiǎn)單實(shí)用的人工邊界條件.

    將式(2)代入到式(1)中,可得

    假定xoy坐標(biāo)系旋轉(zhuǎn)θ角得到新坐標(biāo)系xoy′,其坐標(biāo)軸x′與平面上某處波的傳播方向相一致,忽略另一個(gè)軸方向波的傳播,亦即假定?/?y′=0.兩個(gè)坐標(biāo)系之間有如下關(guān)系式

    由此,可得

    把式(13)代入到式(10~11)中,可得以矩陣表示的方程

    對(duì)Q作特征值分解,Q=P-1ΛP,于是式(14)可以寫成

    把xoy系中PU表達(dá)式代入上式中,可得

    式中 :ux(vs)=ux(x-vsΔt cosθ,y+vsΔt sin θ,t-Δt),其余類似.由式(18)不難得到ux,yy的計(jì)算表達(dá),可以發(fā)現(xiàn)邊界節(jié)點(diǎn)位移由前一時(shí)刻的鄰近節(jié)點(diǎn)的位移計(jì)算確定.在不同的邊界上,應(yīng)選擇不同的θ,例如在邊界法線與 x正方向平行的邊界上,應(yīng)取θ=0°計(jì)算,而在邊界法線方向與y軸的負(fù)方向上一致的邊界上,則應(yīng)取θ=90°計(jì)算.

    3 計(jì)算方法的驗(yàn)證

    為了本文計(jì)算方法的正確性,這里計(jì)算了Lamb問題,并將計(jì)算結(jié)果與解析解[8]進(jìn)行了對(duì)比.假定半無限地基密度 ρ=2 250 kg/m3,其Lamé系數(shù) μ=1.048×103MPa,λ=1.863×103MPa;計(jì)算區(qū)域?yàn)? 000 m×1 000 m的矩陣區(qū)域,采用邊長(zhǎng)為5 m的正方形網(wǎng)格,時(shí)間步長(zhǎng)取為0.002 s;振源為一垂直作用在地表中點(diǎn)的集中荷載,其表達(dá)式為 f(t)=exp(-200(t-0.178)2);計(jì)算采用上述人工邊界條件.圖2同時(shí)給出了地表上距垂直力源距離為500 m的點(diǎn)水平位移的解析解和本文數(shù)值計(jì)算結(jié)果,并給出了無人工邊界時(shí)的數(shù)值解結(jié)果.由圖可知,采用本文所述的計(jì)算方法及人工邊界條件所得的數(shù)值結(jié)果與解析解答基本吻合,由此可見本文方法的正確性;對(duì)比無人工邊界條件計(jì)算結(jié)果與有人工邊界條件計(jì)算結(jié)果,可見本文所述的人工邊界可以較好地透射出外行波,是一種有效的人工邊界條件.

    圖2 水平位移分量的解析解與數(shù)值解的比較

    4 填充溝隔振分析

    振動(dòng)對(duì)人們的生活與工作環(huán)境造成了不可忽視的影響,振動(dòng)已被列七大公害之一.治理環(huán)境振動(dòng)的有效途徑之一是設(shè)置屏障隔振,填充溝就是屏障隔振措施之一.利用上述的計(jì)算方法,本文研究了采用填充溝進(jìn)行隔振時(shí),不同的填充溝深度條件下的豎向振幅衰減變化.

    圖3是計(jì)算模型.均勻彈性半空間地基的密度為ρs=2 000 kg/m3;剪切模量為G s=1.923×107Pa;泊松比v s=0.3.填充溝內(nèi)填充物為橡膠,其密度ρ=1 200 kg/m3;剪切模量為G=9.655×105Pa;泊松比υ=0.48.振源是頻率為5 H z的正弦力,表達(dá)式為 f(t)=sin(10πt).假定填充溝的寬度不變,為2m;深度h變化,分別取10,12,15,18和20m;填充溝距振源為L(zhǎng)1=50m;在填充溝后的地表,選取與填充溝距離為L(zhǎng)2的若干個(gè)觀測(cè)點(diǎn).

    圖3 填充溝計(jì)算示意圖

    Woods[9]研究屏障隔振時(shí),提出在頻域分析中,對(duì)于隔振效果可采用振幅衰減系數(shù) AR來衡量.類似地,本文定義豎向振幅衰減系數(shù)

    用于評(píng)價(jià)隔振效果.根據(jù)這一定義,如果ˉAR值越小,則表明隔振的效果越好,反之則表明隔振效果不佳.

    圖4是計(jì)算結(jié)果.為便于分析,圖中長(zhǎng)度采用量綱一的量的長(zhǎng)度:L*2=L2/λR,h*=h/λR,λR 是R 波波長(zhǎng),λR=18.188 m.由圖可知,對(duì)于所設(shè)置的不同填充溝深度,都可以取得一定的隔振效果.在計(jì)算模型中,隨著溝深度的逐步增加,各點(diǎn)隔振效果先有不同程度地降低,但當(dāng)h*≈λR時(shí),各點(diǎn)隔振效果卻達(dá)到最佳水平.這表明,增加填充溝深度最終是有助于增加隔振效果的.對(duì)于某一填充溝深度而言,填充溝后地表各觀測(cè)點(diǎn)的豎向振幅衰減系數(shù)均在 L*2=1.10λR~3.10λR,L*2=3.50 λR~4.75λR范圍內(nèi)相對(duì)偏大,而在 L*2=3.30λR和L*2=4.95λR處相對(duì)較小.這說明:一方面,填充溝的深度達(dá)一定值后其變化對(duì)于豎向振幅衰減系數(shù)ˉA R影響不大;另一方面,填充溝的隔振效果是隨距離L*2而不斷變化的,在不同的距離處隔振效果相差較大,且存在著某個(gè)點(diǎn),對(duì)應(yīng)著最小豎向振幅衰減系數(shù).所以在設(shè)計(jì)填充溝時(shí),應(yīng)注意選擇填充溝的位置,確保隔振保護(hù)區(qū)域?yàn)樽钚∝Q向振幅衰減系數(shù),以達(dá)到最佳的隔振效果.

    圖4 不同填充溝深度h*時(shí)豎向振幅衰減系數(shù)隨距離L*2變化曲線

    5 結(jié) 論

    1)本文基于以位移為變量彈性波方程,給出了二維格子法的計(jì)算方法,并推導(dǎo)適合該法的人工邊界條件.

    2)經(jīng)典的Lam b算例表明了本文計(jì)算方法的正確性和邊界條件的適應(yīng)性.該方法能很好地適應(yīng)不規(guī)則的邊界條件,且計(jì)算簡(jiǎn)便、準(zhǔn)確,可用于彈性波傳播及散射問題研究.

    3)填充溝屏障隔振分析表明,增加填充溝深度有助于增加隔振效果,但填充溝深度并不影響最小豎向振幅衰減系數(shù)點(diǎn)的位置.

    [1]V irieux J.P-SV w ave p ropagation in heterogeneous media:velocity-stress finite-difference method[J].Geophysics,1986,52(4):888-901.

    [2]裴正林.三維各向同性介質(zhì)彈性波方程交錯(cuò)網(wǎng)格高階有限差分法模擬[J].石油物探,2005,44(4):308-316.

    [3]Chen H,W ang X M,Zhao H B.A rotated staggered grid finite-difference with the absorbing boundary condition of a perfec tly matched layer[J].Chinese Science Bulletin,2006,51(19):2304-2314.

    [4]Tessmer E,Kosloff D,Behle A.Elastic wave propagation simu lation in the p resence of surface topography[J].Geophysical Journal International.1992,108:621-632.

    [5]Zhang JF,Liu T.P-SV-wave propagation in heterogeneousmedia:G ridmethod[J].Geophysical Journal International.1999,136,431-438.

    [6]張劍鋒.彈性波數(shù)值模擬的非規(guī)則網(wǎng)格差分法[J].地球物理學(xué)報(bào),1998,41(增刊):357-366.

    [7]Cook R D.Concepts and applications of finite element analysis[M].New York:John Wiley&Sons,1974.

    [8] Lamb H.On the propagation of tremors over the surface of an elastic so lid[J].Philosophical Transactions o f the Royal Society of London,1904,203:1-44.

    [9]Woods R D,Barnet N E,Sangesser R H.A new tool for soil dynamics[J].Journal of Geotechnical Engineering,1974,100(11):1231-1247.

    猜你喜歡
    差分法衰減系數(shù)邊界條件
    二維粘彈性棒和板問題ADI有限差分法
    一類帶有Stieltjes積分邊界條件的分?jǐn)?shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動(dòng)邊值問題的漸近解
    復(fù)合材料孔隙率的超聲檢測(cè)衰減系數(shù)影響因素
    近岸及內(nèi)陸二類水體漫衰減系數(shù)的遙感反演研究進(jìn)展
    對(duì)《電磁波衰減系數(shù)特性分析》結(jié)果的猜想
    HT250材料超聲探傷中的衰減性探究
    基于SQMR方法的三維CSAMT有限差分法數(shù)值模擬
    帶Robin邊界條件的2維隨機(jī)Ginzburg-Landau方程的吸引子
    有限差分法模擬電梯懸掛系統(tǒng)橫向受迫振動(dòng)
    亚洲国产欧美网| 深夜精品福利| 人妻 亚洲 视频| 午夜精品国产一区二区电影| 亚洲欧美激情在线| a级毛片黄视频| 亚洲熟女毛片儿| 国产日韩欧美视频二区| 亚洲欧美日韩另类电影网站| 中文字幕av电影在线播放| 中国国产av一级| 亚洲情色 制服丝袜| 色视频在线一区二区三区| 丰满乱子伦码专区| 美女扒开内裤让男人捅视频| 亚洲精品国产一区二区精华液| 亚洲国产中文字幕在线视频| 中文字幕最新亚洲高清| 久久久久久久精品精品| 丰满迷人的少妇在线观看| 一级毛片电影观看| 丰满迷人的少妇在线观看| 好男人视频免费观看在线| 久久午夜综合久久蜜桃| 久久久久精品国产欧美久久久 | 女性被躁到高潮视频| a级毛片黄视频| 一级黄片播放器| 国产无遮挡羞羞视频在线观看| 最新在线观看一区二区三区 | 国产精品久久久久久精品古装| 国产 精品1| 久久久久国产精品人妻一区二区| 国产日韩欧美视频二区| 精品国产超薄肉色丝袜足j| 激情视频va一区二区三区| 99热全是精品| 亚洲欧美激情在线| 婷婷色av中文字幕| 中文字幕制服av| 国产日韩欧美亚洲二区| 搡老岳熟女国产| 国产精品麻豆人妻色哟哟久久| tube8黄色片| 亚洲国产毛片av蜜桃av| 午夜福利一区二区在线看| 国产精品久久久av美女十八| 日韩制服丝袜自拍偷拍| 黄色 视频免费看| 亚洲美女搞黄在线观看| 涩涩av久久男人的天堂| 国产免费又黄又爽又色| 黄片小视频在线播放| 老汉色av国产亚洲站长工具| 国产成人精品福利久久| 国产男女超爽视频在线观看| 国产成人精品无人区| 超色免费av| 国产精品成人在线| 中文天堂在线官网| 久久精品亚洲av国产电影网| 看非洲黑人一级黄片| 国产免费福利视频在线观看| 成人毛片60女人毛片免费| 日本黄色日本黄色录像| 99精国产麻豆久久婷婷| 精品免费久久久久久久清纯 | 中文字幕高清在线视频| 成人18禁高潮啪啪吃奶动态图| 免费观看性生交大片5| 色婷婷久久久亚洲欧美| 久久久久网色| 建设人人有责人人尽责人人享有的| 日韩制服丝袜自拍偷拍| 少妇精品久久久久久久| 亚洲国产欧美在线一区| 精品久久久久久电影网| 国产精品人妻久久久影院| 亚洲国产精品一区三区| 久久97久久精品| 99久国产av精品国产电影| 亚洲精品国产av成人精品| 久久久久国产精品人妻一区二区| 纯流量卡能插随身wifi吗| 熟妇人妻不卡中文字幕| 伊人久久国产一区二区| 久久狼人影院| 久久亚洲国产成人精品v| 成年人免费黄色播放视频| 高清黄色对白视频在线免费看| 一区在线观看完整版| e午夜精品久久久久久久| 成人三级做爰电影| 亚洲精品国产色婷婷电影| 99久久综合免费| 国产片特级美女逼逼视频| 欧美国产精品va在线观看不卡| 狂野欧美激情性bbbbbb| 免费黄网站久久成人精品| 午夜免费观看性视频| 女人久久www免费人成看片| 肉色欧美久久久久久久蜜桃| 亚洲精品视频女| 少妇人妻久久综合中文| 国产黄色免费在线视频| 久久精品亚洲av国产电影网| 最近最新中文字幕免费大全7| 亚洲欧洲精品一区二区精品久久久 | 国产欧美亚洲国产| 成人免费观看视频高清| 亚洲国产精品一区二区三区在线| 制服诱惑二区| 十分钟在线观看高清视频www| www.熟女人妻精品国产| 人人妻人人添人人爽欧美一区卜| 99九九在线精品视频| 色婷婷久久久亚洲欧美| 久久久久久久大尺度免费视频| 黑人巨大精品欧美一区二区蜜桃| 日韩精品免费视频一区二区三区| 久久av网站| 精品免费久久久久久久清纯 | 亚洲精品美女久久久久99蜜臀 | 国产av码专区亚洲av| 亚洲精品日韩在线中文字幕| 色吧在线观看| 成人影院久久| 精品久久久精品久久久| 麻豆乱淫一区二区| 飞空精品影院首页| 久久精品aⅴ一区二区三区四区| a级片在线免费高清观看视频| 伦理电影大哥的女人| 亚洲av日韩在线播放| 制服诱惑二区| www.精华液| 亚洲国产最新在线播放| 日韩一区二区三区影片| 波多野结衣av一区二区av| 精品少妇黑人巨大在线播放| 男女免费视频国产| 国产精品无大码| 高清av免费在线| 热re99久久国产66热| 丰满饥渴人妻一区二区三| 涩涩av久久男人的天堂| 亚洲第一青青草原| 国产成人精品久久二区二区91 | 天美传媒精品一区二区| 又黄又粗又硬又大视频| 青春草国产在线视频| 亚洲成人免费av在线播放| 亚洲国产看品久久| 成人免费观看视频高清| 亚洲一区二区三区欧美精品| 亚洲欧美成人综合另类久久久| 一级片'在线观看视频| 久久韩国三级中文字幕| 青春草国产在线视频| 久久国产精品男人的天堂亚洲| www.av在线官网国产| 日韩av不卡免费在线播放| 日韩大片免费观看网站| 亚洲 欧美一区二区三区| 国产 精品1| 男女午夜视频在线观看| 亚洲第一青青草原| 国产成人精品久久二区二区91 | 国语对白做爰xxxⅹ性视频网站| 久久久久久免费高清国产稀缺| 两个人免费观看高清视频| 中文字幕另类日韩欧美亚洲嫩草| 又大又黄又爽视频免费| 91老司机精品| av女优亚洲男人天堂| 国产亚洲av片在线观看秒播厂| 国产探花极品一区二区| 亚洲激情五月婷婷啪啪| 午夜老司机福利片| 亚洲情色 制服丝袜| 九色亚洲精品在线播放| 免费在线观看完整版高清| 欧美在线一区亚洲| 成人国语在线视频| 咕卡用的链子| 高清在线视频一区二区三区| 91精品三级在线观看| 国产精品一区二区在线观看99| 国产一区二区三区av在线| 夫妻午夜视频| 国产片特级美女逼逼视频| 国产亚洲av片在线观看秒播厂| 精品少妇黑人巨大在线播放| 日韩大码丰满熟妇| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产欧美日韩在线播放| 久久婷婷青草| 人体艺术视频欧美日本| 最黄视频免费看| 青青草视频在线视频观看| 欧美精品亚洲一区二区| 国产亚洲av片在线观看秒播厂| 国产在线免费精品| 只有这里有精品99| 午夜久久久在线观看| netflix在线观看网站| 亚洲专区中文字幕在线 | 久久影院123| 久久人人爽人人片av| 考比视频在线观看| 这个男人来自地球电影免费观看 | 国产无遮挡羞羞视频在线观看| 国产av一区二区精品久久| 丝袜美腿诱惑在线| 免费在线观看黄色视频的| 观看美女的网站| 国产亚洲一区二区精品| 久久久精品免费免费高清| 色精品久久人妻99蜜桃| 高清av免费在线| 青春草亚洲视频在线观看| 免费日韩欧美在线观看| 国产成人系列免费观看| 国产人伦9x9x在线观看| 激情视频va一区二区三区| 日日爽夜夜爽网站| 久久久久视频综合| 最黄视频免费看| 亚洲国产av影院在线观看| 七月丁香在线播放| 一区二区三区精品91| 国产精品偷伦视频观看了| 制服诱惑二区| 亚洲天堂av无毛| 免费在线观看黄色视频的| 亚洲国产最新在线播放| 尾随美女入室| 日韩av免费高清视频| 亚洲精品国产区一区二| 国产女主播在线喷水免费视频网站| 久久久久久久精品精品| 国产又爽黄色视频| 久久精品久久久久久噜噜老黄| 国产老妇伦熟女老妇高清| 国产在线视频一区二区| av一本久久久久| 国产在线一区二区三区精| 高清黄色对白视频在线免费看| 国产精品女同一区二区软件| 久久狼人影院| 99热网站在线观看| 国产精品蜜桃在线观看| 亚洲国产欧美一区二区综合| 午夜福利影视在线免费观看| 大香蕉久久网| 制服人妻中文乱码| 男人添女人高潮全过程视频| 少妇人妻久久综合中文| 欧美亚洲日本最大视频资源| 麻豆av在线久日| 亚洲婷婷狠狠爱综合网| 国产探花极品一区二区| www.熟女人妻精品国产| 巨乳人妻的诱惑在线观看| 亚洲熟女精品中文字幕| 在线免费观看不下载黄p国产| 老司机影院毛片| 久久午夜综合久久蜜桃| 一区福利在线观看| 亚洲精品乱久久久久久| 老熟女久久久| 哪个播放器可以免费观看大片| 老司机深夜福利视频在线观看 | 看免费成人av毛片| 中文字幕最新亚洲高清| 精品视频人人做人人爽| 成人国语在线视频| 亚洲精品在线美女| 久久久久久久久久久久大奶| 欧美人与性动交α欧美软件| 女性生殖器流出的白浆| 精品免费久久久久久久清纯 | 黄色视频在线播放观看不卡| 国产成人欧美| 日韩精品免费视频一区二区三区| 精品少妇黑人巨大在线播放| 狠狠婷婷综合久久久久久88av| 十八禁人妻一区二区| 日本vs欧美在线观看视频| 久久影院123| 在线 av 中文字幕| 自线自在国产av| 国产又色又爽无遮挡免| 最近最新中文字幕大全免费视频 | 电影成人av| 搡老乐熟女国产| 亚洲欧美色中文字幕在线| 欧美日韩一级在线毛片| 亚洲精品在线美女| 亚洲精品,欧美精品| 女性被躁到高潮视频| 亚洲av电影在线进入| 天美传媒精品一区二区| 欧美老熟妇乱子伦牲交| 国产精品 欧美亚洲| 最近最新中文字幕大全免费视频 | 国产无遮挡羞羞视频在线观看| 人成视频在线观看免费观看| 中文字幕精品免费在线观看视频| 夫妻性生交免费视频一级片| 久久国产亚洲av麻豆专区| 国产毛片在线视频| 久久影院123| 一区二区日韩欧美中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 国产人伦9x9x在线观看| av视频免费观看在线观看| 一本—道久久a久久精品蜜桃钙片| 久久久久久人妻| 美女大奶头黄色视频| 午夜福利影视在线免费观看| 大陆偷拍与自拍| 国产99久久九九免费精品| 蜜桃在线观看..| 精品午夜福利在线看| 人人妻,人人澡人人爽秒播 | 多毛熟女@视频| 新久久久久国产一级毛片| 免费在线观看黄色视频的| 妹子高潮喷水视频| 只有这里有精品99| 麻豆精品久久久久久蜜桃| 伊人亚洲综合成人网| 中文字幕精品免费在线观看视频| 在线天堂最新版资源| 午夜福利视频精品| 亚洲av日韩精品久久久久久密 | 午夜福利影视在线免费观看| 国产有黄有色有爽视频| av电影中文网址| 亚洲av电影在线进入| 一区二区三区四区激情视频| 国产一区二区三区综合在线观看| 观看美女的网站| 午夜免费男女啪啪视频观看| 久久久久精品国产欧美久久久 | 少妇 在线观看| 国产成人精品久久二区二区91 | 性高湖久久久久久久久免费观看| 亚洲色图综合在线观看| 国产一区二区激情短视频 | 国产片内射在线| 日韩电影二区| √禁漫天堂资源中文www| 免费少妇av软件| 日韩大片免费观看网站| 涩涩av久久男人的天堂| 最新在线观看一区二区三区 | 日韩欧美一区视频在线观看| 欧美国产精品一级二级三级| 国产精品久久久人人做人人爽| 美女扒开内裤让男人捅视频| 两个人免费观看高清视频| 2018国产大陆天天弄谢| 日韩一区二区视频免费看| 一级片'在线观看视频| 一级爰片在线观看| 一级片'在线观看视频| 日韩视频在线欧美| 国产男女内射视频| 亚洲人成网站在线观看播放| av天堂久久9| a级片在线免费高清观看视频| 丰满乱子伦码专区| 午夜福利,免费看| 两个人看的免费小视频| 99国产精品免费福利视频| 青青草视频在线视频观看| 高清视频免费观看一区二区| 精品少妇久久久久久888优播| 十八禁网站网址无遮挡| 国产99久久九九免费精品| 19禁男女啪啪无遮挡网站| av国产精品久久久久影院| 黄色视频不卡| 亚洲国产欧美网| 波野结衣二区三区在线| 亚洲精品国产av成人精品| 国产伦人伦偷精品视频| 中文乱码字字幕精品一区二区三区| 男女之事视频高清在线观看 | 国产精品一二三区在线看| 又大又爽又粗| 亚洲精品成人av观看孕妇| 丝袜人妻中文字幕| 亚洲av综合色区一区| 美女主播在线视频| 韩国精品一区二区三区| 少妇的丰满在线观看| 亚洲欧美色中文字幕在线| 中文精品一卡2卡3卡4更新| 日韩电影二区| 青春草亚洲视频在线观看| 欧美老熟妇乱子伦牲交| 国产av国产精品国产| av有码第一页| 欧美日韩国产mv在线观看视频| 成人亚洲欧美一区二区av| 一区二区三区乱码不卡18| 免费观看人在逋| 2021少妇久久久久久久久久久| 免费观看av网站的网址| av女优亚洲男人天堂| 国产精品 国内视频| 国产精品免费视频内射| 国产精品成人在线| 18禁动态无遮挡网站| 久久久久精品久久久久真实原创| 人人妻,人人澡人人爽秒播 | 精品少妇久久久久久888优播| 国产成人精品福利久久| 男男h啪啪无遮挡| 黑丝袜美女国产一区| 美国免费a级毛片| 亚洲欧美一区二区三区黑人| 亚洲欧美精品综合一区二区三区| 国产精品成人在线| 嫩草影视91久久| 亚洲一级一片aⅴ在线观看| 亚洲第一av免费看| 久久久精品94久久精品| 捣出白浆h1v1| 美国免费a级毛片| 久久99一区二区三区| 水蜜桃什么品种好| 国产精品无大码| 九九爱精品视频在线观看| 日本vs欧美在线观看视频| 搡老岳熟女国产| 国产激情久久老熟女| 国产av一区二区精品久久| 久久影院123| 国产成人精品福利久久| 欧美 亚洲 国产 日韩一| 男女下面插进去视频免费观看| 一二三四在线观看免费中文在| 欧美日韩成人在线一区二区| 亚洲精品国产av蜜桃| 久久精品久久久久久噜噜老黄| 涩涩av久久男人的天堂| 欧美少妇被猛烈插入视频| 国产精品99久久99久久久不卡 | 欧美成人精品欧美一级黄| 国产麻豆69| 男女午夜视频在线观看| www.精华液| 久久亚洲国产成人精品v| 黄色毛片三级朝国网站| 亚洲欧美精品自产自拍| 青春草国产在线视频| av福利片在线| 国产又爽黄色视频| 久久综合国产亚洲精品| 伦理电影大哥的女人| 不卡av一区二区三区| 视频在线观看一区二区三区| 国产精品蜜桃在线观看| 国产男人的电影天堂91| 人人妻人人爽人人添夜夜欢视频| 丝袜美足系列| 亚洲欧美一区二区三区黑人| 天堂俺去俺来也www色官网| bbb黄色大片| 日本欧美视频一区| 亚洲欧洲精品一区二区精品久久久 | 一区二区三区精品91| 自线自在国产av| 亚洲综合色网址| 免费在线观看完整版高清| 久久国产亚洲av麻豆专区| 亚洲欧美一区二区三区黑人| 啦啦啦啦在线视频资源| 亚洲国产av影院在线观看| 精品视频人人做人人爽| 国产99久久九九免费精品| 免费观看人在逋| 亚洲欧美日韩另类电影网站| 欧美日韩一级在线毛片| 日韩 欧美 亚洲 中文字幕| av天堂久久9| 国产免费现黄频在线看| 久久久久精品性色| 99香蕉大伊视频| 国产精品久久久av美女十八| 国产精品女同一区二区软件| 丝袜美腿诱惑在线| 少妇 在线观看| 国产国语露脸激情在线看| av在线老鸭窝| svipshipincom国产片| 岛国毛片在线播放| 男人爽女人下面视频在线观看| 久久99精品国语久久久| 涩涩av久久男人的天堂| 在线观看国产h片| 欧美最新免费一区二区三区| 侵犯人妻中文字幕一二三四区| 亚洲一卡2卡3卡4卡5卡精品中文| av在线老鸭窝| 国产成人av激情在线播放| 黑丝袜美女国产一区| 少妇人妻精品综合一区二区| 人人妻人人澡人人看| 精品视频人人做人人爽| 高清不卡的av网站| 永久免费av网站大全| 999精品在线视频| 一本—道久久a久久精品蜜桃钙片| 国产一卡二卡三卡精品 | 亚洲成av片中文字幕在线观看| 午夜91福利影院| 久久免费观看电影| 老熟女久久久| 欧美精品av麻豆av| 人体艺术视频欧美日本| 日韩一区二区三区影片| 日本vs欧美在线观看视频| 国产又色又爽无遮挡免| 亚洲欧美一区二区三区久久| 免费高清在线观看视频在线观看| 97在线人人人人妻| av天堂久久9| 99精国产麻豆久久婷婷| 天堂俺去俺来也www色官网| 色吧在线观看| 熟女av电影| 国产在视频线精品| 成人手机av| 在线观看人妻少妇| 色视频在线一区二区三区| 精品免费久久久久久久清纯 | 久久精品人人爽人人爽视色| 水蜜桃什么品种好| 国产一区二区在线观看av| 国产精品国产av在线观看| 国产免费一区二区三区四区乱码| 男女边摸边吃奶| 高清在线视频一区二区三区| 成人漫画全彩无遮挡| 叶爱在线成人免费视频播放| 亚洲熟女毛片儿| 亚洲精品乱久久久久久| 日本黄色日本黄色录像| av女优亚洲男人天堂| 人人妻人人澡人人看| 午夜免费男女啪啪视频观看| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩成人在线一区二区| 国产成人系列免费观看| 一区二区三区四区激情视频| 亚洲国产欧美一区二区综合| 国产乱人偷精品视频| 成人18禁高潮啪啪吃奶动态图| 肉色欧美久久久久久久蜜桃| 亚洲精品在线美女| 午夜激情久久久久久久| 五月开心婷婷网| 九色亚洲精品在线播放| 国产一区二区在线观看av| 亚洲,欧美精品.| 国产精品熟女久久久久浪| 在线观看免费午夜福利视频| 久热这里只有精品99| 午夜91福利影院| 国产一区二区三区av在线| 精品视频人人做人人爽| 色播在线永久视频| 欧美日韩成人在线一区二区| 免费少妇av软件| 久久国产亚洲av麻豆专区| 亚洲一区二区三区欧美精品| 极品少妇高潮喷水抽搐| 精品少妇久久久久久888优播| 国产精品久久久久久人妻精品电影 | 丝袜人妻中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久午夜乱码| a级毛片在线看网站| 国产野战对白在线观看| 久久99热这里只频精品6学生| 老鸭窝网址在线观看| av线在线观看网站| 精品一区二区三区av网在线观看 | 超碰成人久久| 国产欧美日韩综合在线一区二区| 另类亚洲欧美激情| 国产精品久久久久成人av| 日本午夜av视频| 久久毛片免费看一区二区三区| 女人爽到高潮嗷嗷叫在线视频| av免费观看日本| 午夜免费观看性视频| 午夜91福利影院| av免费观看日本| 久久久欧美国产精品| 2021少妇久久久久久久久久久| 久久久久精品久久久久真实原创| 18禁国产床啪视频网站| av不卡在线播放| 操出白浆在线播放| 久久天堂一区二区三区四区| 我要看黄色一级片免费的| 日韩伦理黄色片| 亚洲成人国产一区在线观看 | 少妇被粗大的猛进出69影院|