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

    地震作用下水-軸對(duì)稱柱體相互作用的子結(jié)構(gòu)分析方法

    2021-03-02 09:59:12王曉靜王丕光杜修力程星磊
    工程力學(xué) 2021年2期
    關(guān)鍵詞:有限元結(jié)構(gòu)

    王曉靜,趙 密,王丕光,杜修力,程星磊

    (北京工業(yè)大學(xué)城市與工程安全減災(zāi)教育部重點(diǎn)實(shí)驗(yàn)室,北京100124)

    隨著經(jīng)濟(jì)的發(fā)展,近海工程結(jié)構(gòu)的發(fā)展也越來(lái)越受到關(guān)注,如海上風(fēng)電、人工島、跨海橋梁等[1-3]。地震作用下水-結(jié)構(gòu)的相互作用會(huì)在結(jié)構(gòu)的表面產(chǎn)生動(dòng)水壓力作用,該動(dòng)水壓力會(huì)影響結(jié)構(gòu)的動(dòng)力響應(yīng)和動(dòng)力特性[4-6]。因此研究水-柱體結(jié)構(gòu)相互作用對(duì)近海工程結(jié)構(gòu)的抗震設(shè)計(jì)具有重要的意義。

    地震作用下水-結(jié)構(gòu)的相互作用的研究可以追溯到1933年Westergaard 對(duì)豎直壩面上的地震動(dòng)水壓力的求解。地震作用下水-柱體相互作用的研究則始于1965年橢圓截面剛性柱體上動(dòng)水壓力的推導(dǎo)[7]。之后,Chopra 在討論了水體壓縮性和表面波對(duì)圓柱結(jié)構(gòu)地震響應(yīng)的影響[8]。水-柱體結(jié)構(gòu)相互作用的分析方法主要有解析法、數(shù)值法、和物理實(shí)驗(yàn)法[9]等。解析法是通過(guò)首先求解水-柱體結(jié)構(gòu)相互作用的定解問(wèn)題的動(dòng)水壓力解析式,進(jìn)而將動(dòng)水壓力公式代入結(jié)構(gòu)的運(yùn)動(dòng)方程對(duì)結(jié)構(gòu)進(jìn)行動(dòng)力響應(yīng)分析[10-15]。解析法一般只能分析幾何形狀簡(jiǎn)單的等截面柱體結(jié)構(gòu),如圓柱[10-13]和橢圓柱[14-15]。數(shù)值法則是基于有限元[16-17]、邊界元[18-19]、無(wú)限單元[20]等數(shù)值方法分析水-柱體結(jié)構(gòu)相互作用問(wèn)題的方法,數(shù)值方法對(duì)結(jié)構(gòu)的幾何形狀沒有那么高的要求,可以求解幾何形狀復(fù)雜的結(jié)構(gòu)或者傾斜結(jié)構(gòu)。

    軸對(duì)稱柱體結(jié)構(gòu)是近海工程中常見的結(jié)構(gòu)形式,如圓錐柱體、復(fù)合桶形基礎(chǔ)等[17,21-22]。數(shù)值方法可以用于求解地震作用下水-軸對(duì)稱柱體的相互作用問(wèn)題。近海工程中一般將水體介質(zhì)視為無(wú)限域,數(shù)值模型中可通過(guò)引入人工邊界條件將水體分為近場(chǎng)有限域和遠(yuǎn)場(chǎng)無(wú)限域。近場(chǎng)有限域可以用動(dòng)力有限元法進(jìn)行模擬,遠(yuǎn)場(chǎng)無(wú)限域則用人工邊界條件模擬[23-26]。王丕光等[24]提出了一種三維圓柱形高精度時(shí)域人工邊界條件求解三維復(fù)雜結(jié)構(gòu)與水體的動(dòng)力相互作用問(wèn)題。本文則是針對(duì)地震作用下水-軸對(duì)稱柱體的相互作用問(wèn)題的特點(diǎn),采用分離變量法將三維問(wèn)題轉(zhuǎn)變?yōu)橐环N環(huán)向解析、豎向和徑向數(shù)值的二維分析模型;二維分析模型中,基于比例邊界有限元法(SBFEM)[27-29]推導(dǎo)了一種模擬遠(yuǎn)場(chǎng)無(wú)限域的高精度人工邊界。比例邊界有限元法能夠降低問(wèn)題的維度,從而減小計(jì)算量、提高計(jì)算效率,并且能自動(dòng)滿足無(wú)限遠(yuǎn)處的輻射邊界條件。

    1 問(wèn)題描述

    地震作用下水-軸對(duì)稱柱體的動(dòng)力相互作用問(wèn)題如圖1 所示,圖中h表示水深,H為柱體高度;在直角坐標(biāo)系和柱坐標(biāo)系中,z軸為沿著柱體軸線向上,坐標(biāo)原點(diǎn)位于柱體底部;a(z)表示z處截面半徑。假定地基為剛性,地面運(yùn)動(dòng)沿x方向運(yùn)動(dòng),相應(yīng)的位移時(shí)程為ug;結(jié)構(gòu)-水體系統(tǒng)初始靜止;水體假定為小擾動(dòng)不可壓縮流體,并忽略表面重力波的影響,水體的密度為 ρw=1000 kg/m3。柱坐標(biāo)系(r,θ,z)下水體的控制方程為:

    式中,p表示水體的動(dòng)水壓力。

    圖1 軸對(duì)稱結(jié)構(gòu)-水相互作用示意圖Fig.1 Axisymmetric cylinder-water interaction

    水底、水面和水與結(jié)構(gòu)交界面的邊界條件以及無(wú)窮遠(yuǎn)處的輻射條件分別為:

    動(dòng)力壓力p的可以通過(guò)分離變量的方式并結(jié)合邊界條件表示為p=p1(r,z,t)cosθ,則相應(yīng)的軸對(duì)稱結(jié)構(gòu)-水相互作用問(wèn)題的控制方程和邊界條件可以用動(dòng)水壓力的環(huán)向分量p1表示為:

    為方便敘述在下文的推導(dǎo)過(guò)程中p1稱為動(dòng)水壓力。

    在r=R0處引入人工邊界條件,則無(wú)限域水體可以分為近場(chǎng)有限域 ΩI和遠(yuǎn)場(chǎng)無(wú)限域 ΩII。近場(chǎng)有限域可以用軸對(duì)稱有限元離散。

    2 軸對(duì)稱有限元推導(dǎo)

    采用與有限元相同的四節(jié)點(diǎn)四邊形單元等參單元變換,再利用高斯積分即可得到剛度矩陣的值。

    3 高精度人工邊界條件

    3.1 比例邊界坐標(biāo)變換

    基于比例邊界中心線的思想,本文將oz軸選取為比例邊界中心線,如圖2所示,圖中比例邊界坐標(biāo) ξ 和 η,為局部坐標(biāo)。徑向坐標(biāo) ξ可以看作是比例系數(shù),因此域內(nèi)任一點(diǎn)的坐標(biāo)(r,z)可以用局部比例邊界坐標(biāo) ξ 和 η表示為:

    式中:N=N(η)=[(1-η)/2 (1+η)/2]為形函數(shù);zb=[z1z2]T為邊上的節(jié)點(diǎn)坐標(biāo)。

    圖2 SBFE 坐標(biāo)變換Fig.2 Thetransformation of SBFE coordinates

    3.2 比例邊界有限元方程

    對(duì)控制方程式(6)運(yùn)用加權(quán)余量法可得:

    將相應(yīng)的坐標(biāo)變換關(guān)系式(21)、式(22)和式(23)代入式(24),并對(duì)其中的兩項(xiàng)進(jìn)行分部積分得:

    對(duì)動(dòng)水壓力和權(quán)函數(shù)采用相同的形函數(shù)插值:

    式中:p1=p1(ξ)為單元節(jié)點(diǎn)動(dòng)水壓力向量;w=w(ξ)為權(quán)函數(shù)向量。定義B1=b1N,B2=b2N,η,下標(biāo)撇號(hào)表示導(dǎo)數(shù),如N,η即表示形函數(shù)N 對(duì)坐標(biāo)η的導(dǎo)數(shù)。

    將式(26)、式(27)代入式(25),并消去任意權(quán)函數(shù)w和對(duì) ξ的積分可得:

    其中:

    將式(28)沿人工邊界處進(jìn)行裝配,消去系數(shù)為零的項(xiàng),并利用上下表面的邊界條件可得:

    4 水-軸對(duì)稱柱體系統(tǒng)的時(shí)域方程

    4.1 水體有限元方程

    整個(gè)水域的動(dòng)水壓力可以表示為:

    4.2 水-軸對(duì)稱柱體系統(tǒng)的時(shí)域方程

    柱體結(jié)構(gòu)可以用梁?jiǎn)卧M(jìn)行模擬,單元?jiǎng)偠群唾|(zhì)量矩陣可以用一下公式表示[30]:

    5 數(shù)值算例

    本文中算例的結(jié)構(gòu)的彈性模量、密度、泊松比和阻尼比分別為:21 GPa,7800 kg/m3,0.3和0.02,結(jié)構(gòu)的厚度為0.06 m。圖3為本文中選取的作用在剛性地基上的地震動(dòng)時(shí)程曲線,圖4為對(duì)應(yīng)的地震動(dòng)的傅里葉幅值曲線。

    圖3 地震動(dòng)時(shí)程Fig.3 Time history of earthquakes

    5.1 方法驗(yàn)證

    首先,對(duì)軸對(duì)稱結(jié)構(gòu)的等效梁?jiǎn)卧M(jìn)行驗(yàn)證。表1為本文軸對(duì)稱結(jié)構(gòu)的等效梁?jiǎn)卧虯BAQUS中的實(shí)體梁?jiǎn)卧?jì)算的自振頻率的比較,可以看出相對(duì)誤差在5%以內(nèi),因此本文選用的等效變換方法是符合精度要求的。

    其次,對(duì)二維分析模型中的人工邊界的驗(yàn)證。分別選取等截面圓柱和圓錐柱體兩個(gè)模型:模型1等截面圓柱的半徑和高度分別為3 m、80 m,模型2圓錐形柱體的半徑、傾角和高度分別為a=3 m 、θ0=atan(80/1)和80 m。模型1的參考解為解析解,模型2的參考解為20 000×80的擴(kuò)展有限元解,網(wǎng)格大小為5×5。人工邊界的設(shè)置在r=50m 處。圖5為剛性軸對(duì)稱柱體上的動(dòng)水力分布本文解和參考解的對(duì)比。由圖5可以看出,采用本文提出的人工邊界條件得到的本文解與參考解極好地吻合。

    最后,對(duì)本文提出的水-軸對(duì)稱相互作用的子結(jié)構(gòu)分析方法的驗(yàn)證。等截面圓柱的參考解為將柱體結(jié)構(gòu)與動(dòng)水壓力的解析解耦合求得的結(jié)構(gòu)動(dòng)力反應(yīng);圓錐柱體模型的參考解為柱體結(jié)構(gòu)有限元方程與動(dòng)水力的擴(kuò)展有限元方程與柱體結(jié)構(gòu)有限元方程耦合解得的結(jié)構(gòu)動(dòng)力反應(yīng)。在模型1、模型2施加圖3所示的編號(hào)為No.1地震動(dòng)荷載,將結(jié)構(gòu)的頂部位移作為參考量。圖6為結(jié)構(gòu)頂部位移響應(yīng)本文解和參考解的對(duì)比。由圖6可以看出本文方解與參考解吻合較好。

    圖4 地震動(dòng)的傅里葉幅值曲線Fig.4 Furrier amplitude spectrum of ground motions

    表1 等效變換梁?jiǎn)卧蛯?shí)體單元的自振頻率Table 1 The natural frequency of vibration comparison of solid element and equivalent beam

    圖5 剛性柱體上無(wú)量綱動(dòng)水力沿高度分布Fig.5 The distribution of dimensionless hydrodynamic pressure

    5.2 應(yīng)用分析

    計(jì)算了不同水深時(shí)水-軸對(duì)稱結(jié)構(gòu)相互作用系統(tǒng)的地震響應(yīng),分析水深對(duì)水-軸對(duì)稱結(jié)構(gòu)相互作用。引入無(wú)量綱系數(shù)Ru=uim/u0m,其中uim表示不同水深時(shí)水-軸對(duì)稱結(jié)構(gòu)相互作用時(shí)位移響應(yīng)的最大值,u0m表示軸對(duì)稱結(jié)構(gòu)在空氣中時(shí)位移響應(yīng)的最大值。

    以模型2為例,分別計(jì)算了水深為20 m、40 m、60 m 和80 m 時(shí),荷載分別為圖3所示的7條地震動(dòng)時(shí)的結(jié)構(gòu)的動(dòng)力響應(yīng);圖7為Ru隨水深的變化。由圖7可以看出,地震動(dòng)水壓力對(duì)結(jié)構(gòu)地震反應(yīng)的影響隨水深相對(duì)結(jié)構(gòu)高度(h/H)的增加呈增大趨勢(shì),這主要是由于作用在結(jié)構(gòu)上的地震動(dòng)水壓力隨水深的增加而增大。需要指出的是,由于不同地震動(dòng)的頻譜特性不同,因此不同地震動(dòng)作用下動(dòng)水壓力對(duì)結(jié)構(gòu)地震響應(yīng)的影響也不相同。

    定義自振頻率降低率為:

    圖6 柱體頂部位移響應(yīng)Fig.6 Top displacement of the axisymmetric structure

    圖7 R u 與水深的關(guān)系Fig.7 R u vary with different earthquakes

    式中,ωws和ωs分別是指水-結(jié)構(gòu)相互作用系統(tǒng)的自振頻率和結(jié)構(gòu)無(wú)水時(shí)的自振頻率。

    表2為不同水深時(shí)水-軸對(duì)稱相互作用系統(tǒng)的自振頻率。表3為不同水深時(shí)自振頻率的相對(duì)值。由表2和3可以看出當(dāng)水深相對(duì)結(jié)構(gòu)高度較小時(shí)(h/H≤0.25),系統(tǒng)的基頻和結(jié)構(gòu)無(wú)水時(shí)的基頻相差不大,由圖7也可以看出此時(shí)結(jié)構(gòu)的地震響應(yīng)與無(wú)水時(shí)沒有明顯差異;這表明水深相對(duì)結(jié)構(gòu)高度較淺時(shí),地震動(dòng)水壓力對(duì)柱體結(jié)構(gòu)地震響應(yīng)的影響可忽略。隨著水深相對(duì)結(jié)構(gòu)的高度的增大(0.25<h/H≤1),系統(tǒng)的頻率降低率逐漸減小,自振頻率的減低可能會(huì)使得系統(tǒng)的二階自振頻率出現(xiàn)在地震動(dòng)的卓越頻率范圍內(nèi),圖7中也可以看出此時(shí)結(jié)構(gòu)的地震響應(yīng)與無(wú)水時(shí)相差較大;因此水深相對(duì)結(jié)構(gòu)高度較深時(shí),地震動(dòng)水力對(duì)柱體結(jié)構(gòu)的動(dòng)力響應(yīng)和自振特性的影響不可忽略。

    表2 水-軸對(duì)稱結(jié)構(gòu)相互作用系統(tǒng)的自振頻率Table2 Thenatural frequency of vibration of wateraxisymmetric interaction system

    表 3自振頻率降低率Table 3 The reduction ratio of natural vibration frequency

    6 結(jié)論

    基于不可壓縮水體,本文提出了一種地震作用下水-軸對(duì)稱柱體相互作用的動(dòng)力響應(yīng)的子結(jié)構(gòu)分析方法。首先將水體三維模型轉(zhuǎn)化為環(huán)向解析、豎向和徑向離散的二維模型,其次用高精度人工邊界條件模擬無(wú)限域水體的作用,在保證了精度的同時(shí)大大提高了計(jì)算效率。

    數(shù)值算例中分別驗(yàn)證了結(jié)構(gòu)的等效梁?jiǎn)卧?、?dòng)水壓力和水-結(jié)構(gòu)動(dòng)力相互作用耦合模型的精確性。以圓錐柱體為例分析了地震作用下水深對(duì)水-軸對(duì)稱結(jié)構(gòu)相互作用的地震響應(yīng)的影響,數(shù)值算例表明隨著水深的增加,動(dòng)水力對(duì)結(jié)構(gòu)的自振頻率和動(dòng)力響應(yīng)的影響呈增大的趨勢(shì)。

    由于本文提出的子結(jié)構(gòu)可以將水-軸對(duì)稱結(jié)構(gòu)的相互作用通過(guò)附加質(zhì)量的模擬,易于在商業(yè)有限元軟件中實(shí)現(xiàn)和方便于工程應(yīng)用。

    猜你喜歡
    有限元結(jié)構(gòu)
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    新型有機(jī)玻璃在站臺(tái)門的應(yīng)用及有限元分析
    基于有限元的深孔鏜削仿真及分析
    論結(jié)構(gòu)
    基于有限元模型對(duì)踝模擬扭傷機(jī)制的探討
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    論《日出》的結(jié)構(gòu)
    創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長(zhǎng)
    磨削淬硬殘余應(yīng)力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    精品一区二区三区视频在线| 女人久久www免费人成看片| av网站免费在线观看视频| 26uuu在线亚洲综合色| 久久久精品欧美日韩精品| 国产亚洲精品久久久com| 色视频www国产| 亚洲在线观看片| 国产老妇女一区| 久久久久精品久久久久真实原创| 听说在线观看完整版免费高清| 麻豆成人午夜福利视频| 久久精品熟女亚洲av麻豆精品| 欧美日韩一区二区视频在线观看视频在线 | 国产精品成人在线| 美女xxoo啪啪120秒动态图| 在现免费观看毛片| av在线老鸭窝| 大陆偷拍与自拍| 天堂中文最新版在线下载 | 69av精品久久久久久| 狂野欧美白嫩少妇大欣赏| 国产 一区 欧美 日韩| 亚洲成色77777| 国产一区二区在线观看日韩| 精品久久久久久久末码| 亚洲欧美日韩另类电影网站 | 久久99热这里只频精品6学生| 美女cb高潮喷水在线观看| 国产精品久久久久久精品电影小说 | 欧美少妇被猛烈插入视频| 美女xxoo啪啪120秒动态图| 久久午夜福利片| 国产精品久久久久久精品电影小说 | 在线观看国产h片| 国精品久久久久久国模美| 亚洲av免费高清在线观看| 色视频www国产| 亚洲不卡免费看| 免费人成在线观看视频色| 特级一级黄色大片| av福利片在线观看| 国产乱人视频| 少妇被粗大猛烈的视频| 又粗又硬又长又爽又黄的视频| 国产高清国产精品国产三级 | 在线观看三级黄色| 国产亚洲91精品色在线| 日韩av免费高清视频| 久久精品综合一区二区三区| 午夜免费鲁丝| 一区二区三区乱码不卡18| 少妇 在线观看| 91久久精品电影网| 岛国毛片在线播放| 日本黄大片高清| 精品人妻视频免费看| 各种免费的搞黄视频| 亚洲av一区综合| 亚洲国产最新在线播放| 久久精品国产鲁丝片午夜精品| 亚洲经典国产精华液单| 国产成年人精品一区二区| 高清欧美精品videossex| tube8黄色片| 国产老妇女一区| 99久久中文字幕三级久久日本| 成人无遮挡网站| 国产乱人偷精品视频| 国产白丝娇喘喷水9色精品| 免费播放大片免费观看视频在线观看| 禁无遮挡网站| 国产午夜精品一二区理论片| 国产真实伦视频高清在线观看| 久久综合国产亚洲精品| 中文字幕亚洲精品专区| 亚洲电影在线观看av| 97超视频在线观看视频| 最近2019中文字幕mv第一页| 国产精品一区二区性色av| 欧美亚洲 丝袜 人妻 在线| 精品一区在线观看国产| 香蕉精品网在线| 欧美日本视频| 欧美精品一区二区大全| 偷拍熟女少妇极品色| 国产精品一二三区在线看| 人人妻人人爽人人添夜夜欢视频 | 狠狠精品人妻久久久久久综合| 久久鲁丝午夜福利片| 3wmmmm亚洲av在线观看| 久久韩国三级中文字幕| 精品久久久噜噜| 国产爱豆传媒在线观看| 综合色丁香网| 真实男女啪啪啪动态图| 婷婷色麻豆天堂久久| 偷拍熟女少妇极品色| 如何舔出高潮| 国产av码专区亚洲av| 我的女老师完整版在线观看| 国内精品宾馆在线| 久久久a久久爽久久v久久| 欧美三级亚洲精品| 夫妻午夜视频| 久热久热在线精品观看| 大陆偷拍与自拍| 爱豆传媒免费全集在线观看| 免费在线观看成人毛片| 久久久久国产网址| 国产精品人妻久久久影院| 亚洲国产精品成人久久小说| 亚洲,欧美,日韩| av国产免费在线观看| 六月丁香七月| 又黄又爽又刺激的免费视频.| 中文天堂在线官网| 色视频在线一区二区三区| 内地一区二区视频在线| 国产高潮美女av| 精品一区在线观看国产| 国产视频内射| 久久久久久久久久人人人人人人| 国产一区亚洲一区在线观看| 欧美3d第一页| 色网站视频免费| 欧美精品国产亚洲| 热re99久久精品国产66热6| 精品酒店卫生间| 成人国产av品久久久| 成年版毛片免费区| 中文精品一卡2卡3卡4更新| 国产成人aa在线观看| 麻豆成人av视频| 欧美日韩视频精品一区| 午夜福利视频1000在线观看| 中国美白少妇内射xxxbb| 欧美日韩精品成人综合77777| 欧美成人一区二区免费高清观看| 国产淫语在线视频| 精品久久久久久久久亚洲| 一区二区三区免费毛片| 久久久欧美国产精品| 日本色播在线视频| 精品一区二区三区视频在线| 黄色一级大片看看| 爱豆传媒免费全集在线观看| 建设人人有责人人尽责人人享有的 | 男女无遮挡免费网站观看| 久久久久久九九精品二区国产| 欧美另类一区| 狠狠精品人妻久久久久久综合| 国产精品久久久久久久电影| 免费大片18禁| 在线观看一区二区三区激情| 成人黄色视频免费在线看| 久久午夜福利片| 男人添女人高潮全过程视频| av播播在线观看一区| 在线观看免费高清a一片| 亚洲av成人精品一区久久| 亚洲av中文av极速乱| 丝袜美腿在线中文| 国产精品国产三级国产专区5o| 看黄色毛片网站| 美女主播在线视频| 久久久久精品性色| 久久久久国产网址| 国产精品成人在线| 毛片女人毛片| 国产精品一区二区三区四区免费观看| 国产成人福利小说| 久久精品国产亚洲av天美| 国内精品美女久久久久久| 激情 狠狠 欧美| 蜜桃久久精品国产亚洲av| 亚洲欧美精品自产自拍| 校园人妻丝袜中文字幕| 男女无遮挡免费网站观看| 亚洲熟女精品中文字幕| 色婷婷久久久亚洲欧美| 网址你懂的国产日韩在线| 国产国拍精品亚洲av在线观看| 超碰97精品在线观看| 色网站视频免费| 成人特级av手机在线观看| 成人亚洲精品av一区二区| 寂寞人妻少妇视频99o| 成人漫画全彩无遮挡| 亚洲国产精品国产精品| 亚洲精品,欧美精品| 久久午夜福利片| 少妇的逼水好多| 一区二区av电影网| 啦啦啦啦在线视频资源| 欧美精品一区二区大全| 一区二区三区乱码不卡18| kizo精华| 亚洲av福利一区| 婷婷色av中文字幕| 91午夜精品亚洲一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 亚洲在线观看片| 精品人妻视频免费看| 婷婷色综合www| 国产伦在线观看视频一区| av免费在线看不卡| 国产69精品久久久久777片| 精品国产一区二区三区久久久樱花 | 精品人妻视频免费看| 大香蕉97超碰在线| 国产成人精品福利久久| 国产高清国产精品国产三级 | 熟女电影av网| 深爱激情五月婷婷| 日韩电影二区| 国产精品一及| 精品99又大又爽又粗少妇毛片| 在线看a的网站| 亚洲精品久久午夜乱码| av在线亚洲专区| 丰满人妻一区二区三区视频av| 麻豆乱淫一区二区| 天堂中文最新版在线下载 | 国产精品国产av在线观看| 亚洲欧美成人综合另类久久久| 国产精品一区www在线观看| 黄色一级大片看看| 久久综合国产亚洲精品| a级一级毛片免费在线观看| 久久精品久久久久久噜噜老黄| 熟女人妻精品中文字幕| 夜夜看夜夜爽夜夜摸| 午夜免费鲁丝| 亚洲精华国产精华液的使用体验| 五月玫瑰六月丁香| 中文字幕免费在线视频6| av网站免费在线观看视频| 毛片一级片免费看久久久久| 夜夜爽夜夜爽视频| 草草在线视频免费看| 国产美女午夜福利| 三级国产精品欧美在线观看| 嫩草影院入口| 日韩,欧美,国产一区二区三区| 一级片'在线观看视频| 久久久久久久国产电影| 欧美日韩视频精品一区| 青春草国产在线视频| 一级毛片黄色毛片免费观看视频| 性插视频无遮挡在线免费观看| 亚洲人成网站在线观看播放| 日韩av免费高清视频| 99久久人妻综合| 街头女战士在线观看网站| 国产片特级美女逼逼视频| 亚洲色图av天堂| 久久影院123| 午夜精品一区二区三区免费看| av免费观看日本| 狂野欧美白嫩少妇大欣赏| 全区人妻精品视频| 精品久久久久久久末码| 99re6热这里在线精品视频| 大片电影免费在线观看免费| 久久久久久久亚洲中文字幕| 99热国产这里只有精品6| 99热6这里只有精品| 九九在线视频观看精品| 亚洲高清免费不卡视频| 人妻制服诱惑在线中文字幕| 国产综合精华液| 国产欧美日韩精品一区二区| 精品午夜福利在线看| 丝瓜视频免费看黄片| 黄片无遮挡物在线观看| 亚洲精品第二区| 日本免费在线观看一区| 一区二区三区免费毛片| 女人久久www免费人成看片| 我要看日韩黄色一级片| 国产午夜精品一二区理论片| 久久久久久国产a免费观看| 97在线人人人人妻| 婷婷色综合大香蕉| 成人鲁丝片一二三区免费| 欧美精品一区二区大全| 成年人午夜在线观看视频| 亚洲成人av在线免费| 久久97久久精品| 久热这里只有精品99| 国产精品一及| 涩涩av久久男人的天堂| 亚洲欧美成人综合另类久久久| 美女高潮的动态| 久久久亚洲精品成人影院| 婷婷色综合www| h日本视频在线播放| 亚洲av电影在线观看一区二区三区 | 3wmmmm亚洲av在线观看| 香蕉精品网在线| 2022亚洲国产成人精品| 97精品久久久久久久久久精品| 免费观看在线日韩| 国语对白做爰xxxⅹ性视频网站| 黄色配什么色好看| 亚洲av一区综合| 午夜免费男女啪啪视频观看| 国国产精品蜜臀av免费| 久久久精品免费免费高清| 精品午夜福利在线看| av黄色大香蕉| 菩萨蛮人人尽说江南好唐韦庄| 一边亲一边摸免费视频| 国产精品麻豆人妻色哟哟久久| 亚洲精品视频女| av女优亚洲男人天堂| 人体艺术视频欧美日本| 日韩欧美精品免费久久| 国产精品一区二区三区四区免费观看| 欧美zozozo另类| 欧美人与善性xxx| 欧美另类一区| 国产欧美日韩精品一区二区| 国产亚洲午夜精品一区二区久久 | 国产亚洲最大av| 大又大粗又爽又黄少妇毛片口| 国产精品一及| 国产色爽女视频免费观看| 国产伦精品一区二区三区视频9| 亚洲精品中文字幕在线视频 | 亚洲无线观看免费| 18禁在线播放成人免费| 欧美人与善性xxx| 国产精品精品国产色婷婷| 看黄色毛片网站| 男女边吃奶边做爰视频| 赤兔流量卡办理| 一级a做视频免费观看| 日韩视频在线欧美| 一级a做视频免费观看| 中文字幕久久专区| 国产国拍精品亚洲av在线观看| 国产日韩欧美在线精品| 成人国产av品久久久| 亚洲av成人精品一二三区| 久久久久久久大尺度免费视频| 欧美日韩视频精品一区| 18禁动态无遮挡网站| 哪个播放器可以免费观看大片| 天堂俺去俺来也www色官网| 亚洲精品视频女| 国产成人aa在线观看| 久久久久久九九精品二区国产| 精品一区二区三区视频在线| 午夜福利网站1000一区二区三区| 丝袜脚勾引网站| 两个人的视频大全免费| 黄色视频在线播放观看不卡| 国产免费一级a男人的天堂| 亚洲最大成人手机在线| 身体一侧抽搐| 黄色视频在线播放观看不卡| 婷婷色麻豆天堂久久| 亚洲精品第二区| 亚洲一区二区三区欧美精品 | 亚洲精品成人久久久久久| 男人狂女人下面高潮的视频| 26uuu在线亚洲综合色| 亚洲精品国产av蜜桃| 男女下面进入的视频免费午夜| 18禁动态无遮挡网站| eeuss影院久久| 午夜精品国产一区二区电影 | 国产亚洲一区二区精品| 成人黄色视频免费在线看| 免费看日本二区| 哪个播放器可以免费观看大片| 国产v大片淫在线免费观看| 亚洲精品成人久久久久久| 18禁动态无遮挡网站| 国产一区二区在线观看日韩| 亚洲精品国产av蜜桃| 亚洲欧洲国产日韩| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久精品电影| 新久久久久国产一级毛片| 久久久久久伊人网av| 亚洲成人av在线免费| 精品久久久精品久久久| 亚洲欧美精品自产自拍| 亚洲欧美一区二区三区国产| 一区二区三区精品91| 成人二区视频| 日韩不卡一区二区三区视频在线| 精品人妻一区二区三区麻豆| 欧美精品国产亚洲| 国产又色又爽无遮挡免| 视频中文字幕在线观看| 白带黄色成豆腐渣| 春色校园在线视频观看| 22中文网久久字幕| 性色av一级| 日本色播在线视频| 久久久久久久久久久丰满| 美女国产视频在线观看| 男人添女人高潮全过程视频| 精品一区在线观看国产| 不卡视频在线观看欧美| 一区二区三区免费毛片| 亚洲av日韩在线播放| 亚洲成人久久爱视频| 男人舔奶头视频| 国产乱来视频区| 国产 一区精品| 如何舔出高潮| 七月丁香在线播放| 国产欧美日韩精品一区二区| 偷拍熟女少妇极品色| 日韩,欧美,国产一区二区三区| 国产精品女同一区二区软件| 国产精品.久久久| 成人无遮挡网站| 国产午夜福利久久久久久| 亚洲欧美中文字幕日韩二区| 国产精品偷伦视频观看了| 香蕉精品网在线| 国产精品秋霞免费鲁丝片| 最近最新中文字幕大全电影3| 久久人人爽人人爽人人片va| 国产在线男女| 天美传媒精品一区二区| 一个人看的www免费观看视频| 精品久久久精品久久久| 国产爱豆传媒在线观看| 免费看av在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 久久久a久久爽久久v久久| 少妇的逼好多水| 在线播放无遮挡| 免费观看性生交大片5| 国产黄片美女视频| 一本一本综合久久| 国产亚洲5aaaaa淫片| 特大巨黑吊av在线直播| 亚洲欧洲国产日韩| 夜夜爽夜夜爽视频| 国产精品精品国产色婷婷| 人妻制服诱惑在线中文字幕| 男插女下体视频免费在线播放| 91久久精品电影网| 国产精品久久久久久精品古装| 久久久久精品久久久久真实原创| 在线观看一区二区三区激情| 我要看日韩黄色一级片| 三级经典国产精品| 久久久久网色| 91aial.com中文字幕在线观看| 国产爽快片一区二区三区| 国产精品国产三级国产av玫瑰| 成人欧美大片| 欧美日韩视频高清一区二区三区二| 三级经典国产精品| av女优亚洲男人天堂| 2021天堂中文幕一二区在线观| 国产91av在线免费观看| 91狼人影院| 亚洲精品久久午夜乱码| 尾随美女入室| 亚洲精品影视一区二区三区av| 日本黄大片高清| 久久久久性生活片| 亚洲怡红院男人天堂| 国产色婷婷99| 99九九线精品视频在线观看视频| 欧美最新免费一区二区三区| 久久久久精品性色| 哪个播放器可以免费观看大片| 精品人妻一区二区三区麻豆| 国国产精品蜜臀av免费| 国产欧美亚洲国产| 亚洲四区av| 蜜臀久久99精品久久宅男| 午夜精品国产一区二区电影 | 99热这里只有精品一区| 自拍欧美九色日韩亚洲蝌蚪91 | 男女无遮挡免费网站观看| 欧美性猛交╳xxx乱大交人| 三级经典国产精品| 午夜爱爱视频在线播放| 亚洲精品成人久久久久久| 国产 一区 欧美 日韩| 丰满乱子伦码专区| 久久久久久久久久久免费av| 亚洲欧美中文字幕日韩二区| 亚洲成人中文字幕在线播放| 一级爰片在线观看| 偷拍熟女少妇极品色| 成人亚洲欧美一区二区av| 亚洲国产色片| 3wmmmm亚洲av在线观看| 亚洲三级黄色毛片| 精品少妇久久久久久888优播| 国产有黄有色有爽视频| 麻豆乱淫一区二区| 亚洲伊人久久精品综合| 波多野结衣巨乳人妻| 国产高清有码在线观看视频| 身体一侧抽搐| 如何舔出高潮| 啦啦啦在线观看免费高清www| 日日摸夜夜添夜夜添av毛片| 成人无遮挡网站| 亚洲精品视频女| 欧美成人午夜免费资源| 国国产精品蜜臀av免费| 看十八女毛片水多多多| 下体分泌物呈黄色| 在线观看一区二区三区| 欧美潮喷喷水| 欧美激情在线99| 在线看a的网站| 成人二区视频| 自拍欧美九色日韩亚洲蝌蚪91 | 精品久久久久久久人妻蜜臀av| 日日摸夜夜添夜夜爱| 国产精品一区二区性色av| 嘟嘟电影网在线观看| 国产精品人妻久久久影院| 伦精品一区二区三区| 日本一二三区视频观看| 国产精品成人在线| 国产精品麻豆人妻色哟哟久久| 超碰av人人做人人爽久久| 亚洲国产高清在线一区二区三| 日本黄色片子视频| 爱豆传媒免费全集在线观看| 内射极品少妇av片p| 亚洲av免费高清在线观看| 少妇人妻一区二区三区视频| 国产一区有黄有色的免费视频| 视频中文字幕在线观看| 三级国产精品欧美在线观看| 欧美成人精品欧美一级黄| 国产男人的电影天堂91| 久久精品国产亚洲av涩爱| 在现免费观看毛片| 少妇被粗大猛烈的视频| 久热这里只有精品99| 性插视频无遮挡在线免费观看| 男女边摸边吃奶| 久久精品国产亚洲av天美| 激情五月婷婷亚洲| 99精国产麻豆久久婷婷| 国产69精品久久久久777片| 亚州av有码| 搞女人的毛片| 国产老妇伦熟女老妇高清| 五月伊人婷婷丁香| 秋霞伦理黄片| 国产免费一区二区三区四区乱码| 午夜福利在线在线| 国产高清三级在线| 91在线精品国自产拍蜜月| 91久久精品电影网| 亚洲图色成人| 欧美性感艳星| 亚洲欧美一区二区三区黑人 | 国产黄色视频一区二区在线观看| 国产精品av视频在线免费观看| 80岁老熟妇乱子伦牲交| 成人无遮挡网站| 国产午夜精品久久久久久一区二区三区| 久久久久久久亚洲中文字幕| 国产av国产精品国产| 日韩一本色道免费dvd| 久久精品国产鲁丝片午夜精品| 久久久精品欧美日韩精品| 免费观看在线日韩| 中文字幕制服av| 精品一区在线观看国产| 国产精品久久久久久久电影| 国产v大片淫在线免费观看| 日韩av免费高清视频| 免费av不卡在线播放| 天堂俺去俺来也www色官网| 在线观看三级黄色| 国产成人精品福利久久| 2021少妇久久久久久久久久久| 亚洲激情五月婷婷啪啪| 如何舔出高潮| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 色网站视频免费| 中文乱码字字幕精品一区二区三区| 亚洲精品自拍成人| 国产精品久久久久久久久免| 菩萨蛮人人尽说江南好唐韦庄| 欧美一级a爱片免费观看看| 国产视频首页在线观看| 国产欧美亚洲国产| 日本一二三区视频观看| 日韩人妻高清精品专区| 一级av片app| 美女国产视频在线观看| 美女被艹到高潮喷水动态| 18禁在线播放成人免费| 久久精品久久精品一区二区三区| 国产毛片a区久久久久| 午夜老司机福利剧场| 精品人妻偷拍中文字幕| 91精品国产九色| 中国三级夫妇交换| 国产又色又爽无遮挡免| 久久久久久久精品精品| 国产老妇女一区|