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

    基巖與進(jìn)水塔的相互作用對(duì)塔體的地震損傷影響研究

    2022-06-21 10:23:04張?jiān)茲?/span>段錫志郭書(shū)亮鄭曉東
    水力發(fā)電 2022年4期
    關(guān)鍵詞:進(jìn)水塔塔體張開(kāi)

    張?jiān)茲?,段錫志,郭書(shū)亮,鄭曉東

    (1.河北省水利工程局集團(tuán)有限公司,河北 石家莊 050021;2.河北工程大學(xué)水利水電學(xué)院,河北 邯鄲 056038;3.重慶工貿(mào)職業(yè)技術(shù)學(xué)院,重慶 408000)

    0 引 言

    進(jìn)水塔是引水、泄水系統(tǒng)的控制性水工建筑物,它的安全性對(duì)整個(gè)引水、泄水系統(tǒng)甚至大壩的安全都有著舉足輕重的作用[1-3]。對(duì)于如此重要且條件復(fù)雜的高聳結(jié)構(gòu),很多學(xué)者在其穩(wěn)定性方面開(kāi)展了大量的研究[4-8]。曹偉[9]對(duì)于進(jìn)水塔在強(qiáng)震作用下局部混凝土進(jìn)入塑性階段產(chǎn)生破壞,建立混凝土拉壓彈塑性本構(gòu)關(guān)系,考慮混凝土材料的非線性對(duì)進(jìn)水塔進(jìn)行非線性時(shí)程分析。陳震等[10]模擬在強(qiáng)震作用下進(jìn)水塔內(nèi)部累計(jì)損傷導(dǎo)致的塔體變形,引入損傷彈塑性模型,對(duì)混凝土材料非線性進(jìn)行研究分析。唐青山等[11]采用有限差分法分析了回填過(guò)程混凝土和邊坡的穩(wěn)定性,主要考慮了回填混凝土與邊坡接觸面(巖-混凝土接觸面)的相互作用。李騫等[1]針對(duì)地震作用下應(yīng)力釋放,考慮塔體線彈性材料,引入動(dòng)態(tài)接觸方法和彈性連接研究了局部接觸非線性問(wèn)題,說(shuō)明了考慮局部接觸非線性的必要性。但是,專(zhuān)家們分析高聳結(jié)構(gòu)與地基巖體之間的相互作用時(shí),該結(jié)構(gòu)均按彈性材料進(jìn)行分析。

    塔體與巖體基礎(chǔ)之間的動(dòng)力相互作用和材料非線性是影響進(jìn)水塔抗震穩(wěn)定性的2個(gè)重要因素。對(duì)于進(jìn)水塔結(jié)構(gòu)抗震問(wèn)題,大多數(shù)研究都單一的考慮材料非線性因素的影響或者接觸非線性因素的影響,這樣對(duì)實(shí)際進(jìn)水塔的結(jié)構(gòu)和受力情況都不能夠準(zhǔn)確來(lái)描述。因此,為了更加接近進(jìn)水塔的真實(shí)實(shí)際情況,本文在考慮進(jìn)水塔材料非線性的基礎(chǔ)上,研究羊曲進(jìn)水塔在地震荷載作用下塔體與巖體基礎(chǔ)之間的動(dòng)力相互作用,利用面-面接觸單元模擬進(jìn)水塔塔體與基巖各個(gè)面之間的接觸非線性特性,分析地震動(dòng)持續(xù)作用下塔體頂部相對(duì)位移、塔體結(jié)構(gòu)的損傷過(guò)程以及塔體與基巖之間的接觸狀態(tài),從而為進(jìn)水塔的設(shè)計(jì)提供科學(xué)依據(jù)。

    1 工程實(shí)例與有限元模型

    1.1 工程概況及參數(shù)的選取

    本文運(yùn)用ABAQUS建立了羊曲進(jìn)水塔的塔體-地基有限元模型(如圖1所示),鋼筋本構(gòu)采用隨動(dòng)強(qiáng)化模型。進(jìn)水塔中的鋼筋采用三維線單元,鑲嵌到混凝土中,鋼筋與混凝土之間沒(méi)有滑動(dòng)。地基的寬度在塔的左岸、右岸各取28.0 m;塔前后各取40.0 m;深度取接近于塔體高度的1倍,即80.0 m;圍巖按照地基寬度取28.0 m。

    圖1 進(jìn)水塔塔體-地基整體模型示意

    1.1.1 塔體混凝土材料

    羊曲進(jìn)水塔塔座(2 663.00 m高程以下)混凝土強(qiáng)度等級(jí)為C30;塔筒(2 663.00 m高程以上)混凝土強(qiáng)度等級(jí)為C25;塔后回填混凝土強(qiáng)度等級(jí)為C20,見(jiàn)圖1a,混凝土材料動(dòng)態(tài)彈性模量值按照GB 51247—2018《水工建筑物抗震設(shè)計(jì)規(guī)范》確定,相應(yīng)的材料系數(shù)較靜態(tài)標(biāo)準(zhǔn)值提高50%,動(dòng)態(tài)抗拉強(qiáng)度值較靜力參數(shù)增加20%,混凝土材料具體參數(shù)見(jiàn)表1。

    表1 混凝土有關(guān)材料參數(shù)

    1.1.2 基巖材料參數(shù)

    羊曲進(jìn)水塔地基巖體以Ⅱ類(lèi)中厚層砂巖為主,堅(jiān)固系數(shù)為10。地基巖體的密度為2 720 kg/m3,地基變形模量為15 GPa,泊松比為0.167,摩擦系數(shù)為1.07,黏聚力取值為2 MPa,斷裂能為82 N/m。

    1.1.3 塑性損傷模型

    郭勝山[12]在Lee and Fenves基礎(chǔ)上建立了考慮損傷引起的不可恢復(fù)殘余變量的塑性損傷模型,塑性損傷模型本構(gòu)關(guān)系如圖2所示。

    圖2 塑性損傷模型本構(gòu)關(guān)系

    圖2中的塑性損傷模型本構(gòu)關(guān)系為

    ε=εe+εp

    (1)

    εe=εel+εenl

    (2)

    εd=εp+εenl

    (3)

    式中,ε為總應(yīng)變量;εe為彈性應(yīng)變;εP為不可恢復(fù)殘余應(yīng)變;εel為初始彈性應(yīng)變;εenl為剛度減小引起的非彈性應(yīng)變;εd為損傷應(yīng)變。

    塑性損傷模型應(yīng)力應(yīng)變關(guān)系

    本構(gòu)關(guān)系σ=(1-d)D(ε-εP)

    (4)

    (5)

    損傷演化關(guān)系為

    (6)

    1.2 荷載

    主要考慮的荷載有靜水壓力,動(dòng)水壓力,揚(yáng)壓力,塔體自重和地震荷載等。動(dòng)水壓力采用GB 51247—2018《水工建筑物抗震設(shè)計(jì)規(guī)范》中附加質(zhì)量模型進(jìn)行計(jì)算。對(duì)進(jìn)水塔考慮3個(gè)方向的地震荷載作用,本文通過(guò)極限抗震能力計(jì)算,確定出2倍的設(shè)計(jì)地震是羊曲進(jìn)水塔“破壞”前的最大地震荷載[13],即地震峰值加速度為0.608g。調(diào)整后的人工波作為設(shè)計(jì)人工地震波如圖3所示。

    圖3 設(shè)計(jì)人工地震波

    1.3 接觸力學(xué)模型

    接觸非線性屬于強(qiáng)非線性問(wèn)題,接觸狀態(tài)分為張開(kāi)、粘結(jié)、滑移等[14]。本文運(yùn)用面-面接觸方法來(lái)模擬進(jìn)水塔塔底板面與基巖、以及進(jìn)水塔塔后回填混凝土部分與基巖之間的接觸非線性。在進(jìn)水塔結(jié)構(gòu)和基巖受力面的法向方向設(shè)置為硬接觸,受力面的切向方向設(shè)置庫(kù)倫摩擦系數(shù)。

    1.3.1 接觸面約束條件

    (1)法向接觸條件[9]。①法向的不可慣入性,假設(shè)接觸面上點(diǎn)A、B的坐標(biāo)分別為uA、uB,定義點(diǎn)A指向點(diǎn)B的法向矢量為n,定義兩點(diǎn)方向距離為dn,則dn=(uA-uB)n,若dn>0,表示分離;若dn<0,表示相互侵入;若dn=0,表示接觸狀態(tài)。②法向不貫穿條件為dn=(uA-uB)n≥0。

    (2)切向接觸條件[9]。定義兩點(diǎn)接觸力中切向接觸力為λs,以庫(kù)倫模型考慮摩擦,定義界面摩擦系數(shù)為μ。當(dāng)兩點(diǎn)處于接觸狀態(tài)時(shí),若|λs|=μλn表示兩點(diǎn)將發(fā)生滑移;若|λs|<μλn處于粘結(jié)狀態(tài)。即

    粘結(jié)狀態(tài)

    (7)

    滑移狀態(tài)

    (8)

    1.3.2 罰函數(shù)接觸模型

    罰函數(shù)法是基于罰函數(shù)定義切向行為的庫(kù)倫摩擦準(zhǔn)則。罰函數(shù)可滿(mǎn)足法向不嵌入條件,罰函數(shù)摩擦允許接觸面有彈性滑移,能夠有效處理過(guò)約束問(wèn)題,適用于大多數(shù)的接觸問(wèn)題,本文選用罰函數(shù)法。通過(guò)式(9)來(lái)定義臨界剪應(yīng)力與接觸壓力之間的關(guān)系,即

    τcrit=μp

    (9)

    式中,τcrit為材料的剪切屈服應(yīng)力;p為接觸面間接觸壓力;μ為摩擦系數(shù)。接觸單元目標(biāo)單元幾何形狀為矩形,初始接觸狀態(tài)為閉合狀態(tài)。

    2 計(jì)算分析

    2.1 進(jìn)水塔結(jié)構(gòu)自振特性

    通過(guò)對(duì)進(jìn)水塔結(jié)構(gòu)的自振特性進(jìn)行分析求解水塔結(jié)構(gòu)的動(dòng)力特性,運(yùn)用ABAQUS有限元分析軟件進(jìn)行模態(tài)分析可以得到結(jié)構(gòu)的自振周期和自振頻率。羊曲進(jìn)水塔的前5階的自振頻率如表2所示。

    表2 進(jìn)水塔自振頻率

    進(jìn)水塔結(jié)構(gòu)阻尼采用瑞利阻尼,取進(jìn)水塔前5階振型頻率,阻尼比取7%,計(jì)算得到的滿(mǎn)庫(kù)情況的瑞利阻尼參數(shù)α=0.152 4;β=0.032。

    2.2 接觸狀態(tài)分析

    在地震荷載作用下,對(duì)進(jìn)水塔與基巖之間的接觸狀態(tài)進(jìn)行分析,塔體底板面、塔體后面、回填混凝土部分和塔體左右面與基巖接觸面狀態(tài)云圖如圖4、5所示。

    圖4 塔體底板面、塔體后面和回填混凝土與基巖接觸面狀態(tài)云圖

    由圖4中塔體底板面與基巖接觸面狀態(tài)云圖可知,在地震荷載作用下,8.52 s時(shí)進(jìn)水塔塔體底板面與巖基之間接觸面首次被拉裂。在塔體x向位移最大時(shí)刻(17.23 s),最大接觸壓力出現(xiàn)在底板底面靠上游處,塔體底板面與巖基之間后部分接觸面處于張開(kāi)狀態(tài),張開(kāi)值隨著位移值的逐漸增加從塔體與基巖接觸面的左側(cè)到右側(cè)呈現(xiàn)逐漸遞增。在y向最大位移最大時(shí)刻(23.15 s),最大接觸壓力出現(xiàn)在底板底面靠上游處,塔體底板面與巖基之間接觸的薄弱面處于張開(kāi)狀態(tài)。張開(kāi)值從塔體與基巖接觸面的左側(cè)到右側(cè)呈現(xiàn)遞減。

    由圖4中回填混凝土與基巖接觸面狀態(tài)云圖可知,在地震荷載的持續(xù)作用下,在8.52 s和17.23 s,塔后回填混凝土與基巖之間處于完全張開(kāi)狀態(tài)。在23.15 s和30.0 s,塔體被完全貫穿破壞等時(shí)間點(diǎn),其接觸面右側(cè)有小部分處于閉合接觸狀態(tài)。

    由圖4中塔體后面基巖接觸面狀態(tài)云圖可知,在8.52、23.15、30.0 s,塔體被完全貫穿破壞等時(shí)間點(diǎn)塔體后面基巖接觸面張開(kāi)狀態(tài)變化幾乎相同,在塔體x向位移最大時(shí)刻(17.23 s),塔體后面與基巖接觸面張開(kāi)狀態(tài)增大。造成以上兩種現(xiàn)象是由于在考慮塔體-基巖體系和塔后回填混凝土-基巖之間存在接觸非線性,他們之間通過(guò)摩擦系數(shù)和黏聚力約束,在靜、動(dòng)荷載的持續(xù)作用下,塔后回填混凝土與基巖之間處于完全張開(kāi)狀態(tài)。隨著x向地震強(qiáng)度的逐漸增加,使塔體在x向有較大的前傾。

    由圖5可知,塔體在x向最大位移時(shí)刻(17.23 s)和y向最大位移最大時(shí)刻(23.15 s)左面處于張開(kāi)狀態(tài)的面積大于右面,左側(cè)面在17.23~23.15 s之間處于張開(kāi)狀態(tài)的面積呈現(xiàn)逐漸遞增。張開(kāi)狀態(tài)的位置主要處于接觸面的上半部分的大部分面積。

    圖5 塔體左右面與基巖接觸面狀態(tài)云圖

    2.3 塔頂位移響應(yīng)分析

    考慮塔體-基巖體系之間有可靠的連接,不考慮接觸非線性和考慮接觸非線性?xún)煞N情況,分別對(duì)進(jìn)水塔結(jié)構(gòu)進(jìn)行時(shí)程分析,同時(shí)輸入3個(gè)方向的地震荷載,分析塔體在地震荷載作用下塔頂相對(duì)位移隨時(shí)間變化的規(guī)律,得到塔頂相對(duì)位移時(shí)程曲線,如圖6所示。

    圖6 塔體頂部相對(duì)位移時(shí)程曲線

    由圖6a可知,塔體x向相對(duì)位移時(shí)程曲線在14.0 s左右出現(xiàn)明顯變化。不考慮接觸非線性情況下塔體頂部發(fā)生的相對(duì)位移幅值為65.11 mm??紤]接觸非線性情況下,塔體x向的相對(duì)位移幅值為85.37 mm,與不考慮接觸非線性相比,考慮接觸非線性時(shí)塔頂最大位移幅值增加31.12%。這是由于塔體與基巖之間不考慮接觸非線性情況時(shí),將塔體、塔后回填混凝土與周?chē)鷰r體作為一個(gè)整體進(jìn)行研究,塔體和塔后回填混凝土相接觸部分容易產(chǎn)生應(yīng)力集中現(xiàn)象,導(dǎo)致塔體產(chǎn)生損傷開(kāi)裂,使塔體x向產(chǎn)生較小的位移幅值。考慮塔體、塔后回填混凝土與基巖之間存在接觸非線性,接觸面間主要是通過(guò)摩擦力來(lái)減小塔體順?biāo)鞣较虻奈灰?,?dāng)動(dòng)靜荷載的疊加作用力大于摩擦阻力的時(shí),塔后回填混凝土與基巖之間的接觸面局部就會(huì)處于張開(kāi)狀態(tài),使得塔體在順?biāo)鞣较虍a(chǎn)生較大塔頂相對(duì)位移。

    由圖6b可知,考慮塔體與基巖間存在接觸非線性和不存在接觸非線性?xún)煞N情況,塔體頂部相對(duì)位移變化趨勢(shì)幾乎同步,這是由于在兩種情況下,垂直水流y向上塔體與巖基之間都受法向應(yīng)力的約束作用,塔體與基巖之間的接觸面處于閉合狀態(tài),因此,使塔體頂部發(fā)生相對(duì)位移變化趨勢(shì)幾乎同步。

    由圖6c可知,考慮塔體與基巖之間接觸非線性和不考慮接觸非線性?xún)煞N情況下,塔體在地震荷載加載初期,兩種情況下的塔頂位移變化基本是一致的,從14.0 s前后塔頂相對(duì)位移出現(xiàn)逐漸遞增趨勢(shì),考慮接觸非線性時(shí)相對(duì)位移遞增速率更快,25.0 s左右塔頂相對(duì)位移達(dá)到最大值。考慮接觸非線性情況下塔體頂部最大相對(duì)位移值約為38.00 mm,不考慮接觸非線性情況下塔體頂部最大相對(duì)位移值為17.0 mm,考慮接觸非線性情況下相對(duì)于不考慮接觸非線性情況下最大位移值增長(zhǎng)123.52%。這是由于地震荷載加載初期,荷載強(qiáng)度較小,塔體處于彈性階段,所以?xún)煞N情況下初期塔頂相對(duì)位移的變化基本是一致的。隨著地震作用的逐漸加強(qiáng),在不考慮接觸非線性情況下,將塔體-基巖體系看做一個(gè)整體,塔體z向相對(duì)位移小,而考慮接觸非線性情況下,塔體與基巖之間的接觸面局部處于張開(kāi)狀態(tài),增大了塔體z向相對(duì)位移。

    2.4 進(jìn)水塔塔體拉損傷過(guò)程分析

    假設(shè)塔體與基巖之間有可靠的連接,不考慮非線性接觸情況和考慮塔體與基巖之間存在非線性接觸兩種情況計(jì)算得到的塔體腰部損傷云圖如圖7所示。

    由圖7可知,在不考慮接觸非線性情況下,進(jìn)水塔塔體在高程2 663.0 m出現(xiàn)貫穿損傷區(qū)域,該區(qū)域?yàn)椴煌炷敛牧辖唤用嫣帲鋽U(kuò)展方向從上游到下游。在考慮接觸非線性情況下,塔體損傷區(qū)域包括高程2 663.0 m混凝土材料交界面處和引水隧洞上表面幾何突變處。在考慮接觸非線性條件下?lián)p傷突變比不考慮接觸非線性損傷突變損傷破壞單元要多,塔體腰部損傷在在某一時(shí)間段內(nèi)呈現(xiàn)突發(fā)性的增長(zhǎng),總損傷程度在某一時(shí)間段逐漸累積。造成這一現(xiàn)象的主要原因是進(jìn)水塔塔體與基巖之間的相互作用,擴(kuò)大了地震能量,增加了塔體的實(shí)際損傷開(kāi)裂。

    圖7 進(jìn)水塔拉損傷發(fā)展過(guò)程演化

    圖8為進(jìn)水塔損傷耗散能時(shí)程曲線。由圖8可知,進(jìn)水塔-巖基系統(tǒng)在地震荷載作用下被完全貫穿損傷破壞,考慮接觸非線性模型計(jì)算的損傷耗散能為171 758.0 J,不考慮接觸非線性模型計(jì)算的損傷耗散能為113 251.0 J,考慮接觸非線性情況相對(duì)于不考慮接觸非線性情況下的損傷耗散能增長(zhǎng)了34.06%。

    圖8 進(jìn)水塔損傷耗散能時(shí)程曲線

    綜上所述,進(jìn)水塔-基巖體系在地震荷載作用下,不考慮接觸非線性情況,塔體損傷破壞比較慢,損傷單元較少;考慮接觸非線性情況,擴(kuò)大了地震能量,塔體損傷破壞比較快,損傷單元較多。因此,考慮接觸非線性情況下,進(jìn)水塔損傷模型偏不安全,但是進(jìn)水塔模型更加接近真實(shí)實(shí)際,可以為進(jìn)水塔及相近高聳混凝土建筑的設(shè)計(jì)提供更加合理的建議。

    3 結(jié) 論

    本文在考慮塔體材料非線性的基礎(chǔ)上,通過(guò)對(duì)比分析進(jìn)水塔與基巖之間考慮接觸非線性和不考慮接觸非線性2種方案,得到如下結(jié)論:

    (1)在不考慮接觸非線性情況下,進(jìn)水塔在塔后回填混凝土部分會(huì)產(chǎn)生應(yīng)力集中現(xiàn)象,塔體-基巖體系在應(yīng)力集中的影響下?lián)p傷單元主要集中于不同混凝土材料交界面處??紤]塔體與基巖之間存在接觸非線性,無(wú)應(yīng)力集中現(xiàn)象產(chǎn)生,塔體損傷主要集中于不同混凝土材料交界面處和引水隧洞上表面幾何突變處。

    (2)通過(guò)對(duì)塔體位移、損傷、能量等方面的分析,相對(duì)于不考慮接觸非線性情況下,考慮接觸非線性情況的塔頂最大相對(duì)位移幅值增加了31.12%,塔體的損傷區(qū)域擴(kuò)大,損傷耗散能增長(zhǎng)了34.06%,這是由于考慮進(jìn)水塔塔體與基巖之間的相互作用,擴(kuò)大了地震能量,增加了塔體的實(shí)際損傷開(kāi)裂。

    (3)考慮進(jìn)水塔與基巖接觸非線性的情況下,進(jìn)水塔模型偏不安全,但是該情況下進(jìn)水塔模型更接近于真實(shí)情況,可為進(jìn)水塔設(shè)計(jì)提供更加合理的建議。隨著地震強(qiáng)度逐漸增大,進(jìn)水塔與基巖基礎(chǔ)面會(huì)產(chǎn)生局部瞬態(tài)拉裂,為提高建筑物的穩(wěn)定性,可在開(kāi)裂部位進(jìn)行適當(dāng)配筋。該研究成果對(duì)實(shí)際工程設(shè)計(jì)具有一定的指導(dǎo)意義。

    猜你喜歡
    進(jìn)水塔塔體張開(kāi)
    外載荷作用下塔器開(kāi)孔補(bǔ)強(qiáng)計(jì)算探討
    含啟閉機(jī)房進(jìn)水塔結(jié)構(gòu)靜動(dòng)力分析
    某水電站進(jìn)水塔靜力穩(wěn)定性分析
    Preliminary Design of a Submerged Support Structure for Floating Wind Turbines
    復(fù)雜施工環(huán)境下的進(jìn)水塔安裝
    基于振型分解反應(yīng)譜法的深孔泄洪洞進(jìn)水塔動(dòng)力分析
    開(kāi)花
    詩(shī)潮(2019年10期)2019-11-19 13:58:55
    冷卻塔爆破拆除傾倒解體及振動(dòng)研究
    爆炸與沖擊(2019年2期)2019-02-27 02:25:00
    塔體現(xiàn)場(chǎng)改造技術(shù)
    午夜福利视频1000在线观看| 老熟妇乱子伦视频在线观看| 亚洲片人在线观看| 午夜视频国产福利| 69人妻影院| 日韩欧美精品v在线| .国产精品久久| 十八禁国产超污无遮挡网站| bbb黄色大片| 一区二区三区四区激情视频 | 精品午夜福利在线看| 亚洲成人精品中文字幕电影| 国产乱人伦免费视频| 啦啦啦观看免费观看视频高清| 国产淫片久久久久久久久 | 自拍偷自拍亚洲精品老妇| 国产精品久久久久久久电影| 91久久精品国产一区二区成人| www.色视频.com| 日本免费一区二区三区高清不卡| 搞女人的毛片| 国产又黄又爽又无遮挡在线| 国产亚洲欧美98| 免费大片18禁| 国产极品精品免费视频能看的| 免费在线观看日本一区| 免费在线观看亚洲国产| 丝袜美腿在线中文| 国产亚洲欧美在线一区二区| 国产激情偷乱视频一区二区| 欧美丝袜亚洲另类 | 97超视频在线观看视频| 宅男免费午夜| 18禁黄网站禁片免费观看直播| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕免费在线视频6| 99热6这里只有精品| 国模一区二区三区四区视频| 床上黄色一级片| 少妇人妻一区二区三区视频| 小说图片视频综合网站| 午夜免费男女啪啪视频观看 | 99热这里只有精品一区| 免费看美女性在线毛片视频| 国内少妇人妻偷人精品xxx网站| 成熟少妇高潮喷水视频| 国产精华一区二区三区| 日本 av在线| 91九色精品人成在线观看| 国产成人欧美在线观看| 免费人成视频x8x8入口观看| 成年女人永久免费观看视频| 欧美在线黄色| 91久久精品电影网| 人人妻人人看人人澡| 如何舔出高潮| 99国产综合亚洲精品| 久久精品国产亚洲av香蕉五月| 特级一级黄色大片| 国产在线精品亚洲第一网站| 久久久国产成人免费| 亚洲狠狠婷婷综合久久图片| eeuss影院久久| 国产视频一区二区在线看| 亚洲国产日韩欧美精品在线观看| 欧美日韩亚洲国产一区二区在线观看| 十八禁人妻一区二区| 乱人视频在线观看| 亚洲18禁久久av| 亚洲精品一卡2卡三卡4卡5卡| 久久人妻av系列| 成年免费大片在线观看| 中亚洲国语对白在线视频| 国产老妇女一区| 亚洲av中文字字幕乱码综合| 成人亚洲精品av一区二区| 久久九九热精品免费| a级毛片a级免费在线| 午夜影院日韩av| 久久热精品热| 变态另类成人亚洲欧美熟女| 97超级碰碰碰精品色视频在线观看| 欧美激情在线99| 国产高潮美女av| 搞女人的毛片| 中文字幕高清在线视频| 真实男女啪啪啪动态图| 国产麻豆成人av免费视频| 久久久久九九精品影院| 超碰av人人做人人爽久久| 男女床上黄色一级片免费看| 亚洲自拍偷在线| 夜夜看夜夜爽夜夜摸| 黄色一级大片看看| 日日摸夜夜添夜夜添av毛片 | 婷婷六月久久综合丁香| 亚洲精品亚洲一区二区| 日韩中文字幕欧美一区二区| 少妇高潮的动态图| 9191精品国产免费久久| 欧美最黄视频在线播放免费| 国产色爽女视频免费观看| 美女高潮的动态| 99久久无色码亚洲精品果冻| 嫩草影院入口| 69人妻影院| 乱人视频在线观看| 国产国拍精品亚洲av在线观看| 亚洲精品亚洲一区二区| 在线天堂最新版资源| 婷婷精品国产亚洲av| 日韩国内少妇激情av| 黄色配什么色好看| 成年女人看的毛片在线观看| 亚洲第一区二区三区不卡| 三级国产精品欧美在线观看| 欧美日韩综合久久久久久 | a在线观看视频网站| 90打野战视频偷拍视频| 深夜a级毛片| 亚洲乱码一区二区免费版| 色哟哟哟哟哟哟| 日本成人三级电影网站| 国产黄片美女视频| 熟妇人妻久久中文字幕3abv| 麻豆一二三区av精品| 国产在线精品亚洲第一网站| 校园春色视频在线观看| 久久久国产成人精品二区| 午夜激情福利司机影院| 免费观看的影片在线观看| 欧美xxxx黑人xx丫x性爽| 51午夜福利影视在线观看| 亚洲av日韩精品久久久久久密| 精品国产三级普通话版| 精品人妻视频免费看| 波多野结衣高清无吗| 亚洲专区中文字幕在线| 夜夜看夜夜爽夜夜摸| 麻豆成人av在线观看| 免费无遮挡裸体视频| 国产精品野战在线观看| 国内精品久久久久精免费| 成人美女网站在线观看视频| 国产精品永久免费网站| 国产精品永久免费网站| 国产精品久久电影中文字幕| 久久草成人影院| 免费黄网站久久成人精品 | 国产在线男女| 蜜桃久久精品国产亚洲av| 国产 一区 欧美 日韩| 久久精品国产亚洲av天美| 最后的刺客免费高清国语| 中亚洲国语对白在线视频| 岛国在线免费视频观看| 国产精品女同一区二区软件 | 日韩欧美精品v在线| av中文乱码字幕在线| 大型黄色视频在线免费观看| 精品人妻偷拍中文字幕| 精品人妻偷拍中文字幕| 麻豆久久精品国产亚洲av| av在线蜜桃| 一区二区三区高清视频在线| h日本视频在线播放| 白带黄色成豆腐渣| 亚洲欧美日韩东京热| 久久久久久久久久成人| 亚洲国产日韩欧美精品在线观看| 99精品在免费线老司机午夜| 赤兔流量卡办理| 欧美乱色亚洲激情| 成人精品一区二区免费| 亚洲三级黄色毛片| 欧美最黄视频在线播放免费| 又黄又爽又免费观看的视频| 看十八女毛片水多多多| 夜夜看夜夜爽夜夜摸| 精品一区二区三区人妻视频| 1000部很黄的大片| 久久热精品热| 男人舔女人下体高潮全视频| 亚洲精品日韩av片在线观看| 白带黄色成豆腐渣| 欧美区成人在线视频| 亚洲第一区二区三区不卡| 人妻夜夜爽99麻豆av| 久久久久久九九精品二区国产| 变态另类成人亚洲欧美熟女| av福利片在线观看| 日韩欧美 国产精品| 97碰自拍视频| 丰满的人妻完整版| 亚洲av.av天堂| 五月玫瑰六月丁香| 成年女人永久免费观看视频| or卡值多少钱| 国产黄色小视频在线观看| 亚洲欧美日韩卡通动漫| 国产精品久久视频播放| 高清毛片免费观看视频网站| 色综合站精品国产| 国产成人福利小说| 亚洲成人精品中文字幕电影| 婷婷精品国产亚洲av在线| 男人舔奶头视频| 人人妻人人澡欧美一区二区| 麻豆一二三区av精品| 精品久久久久久久人妻蜜臀av| 99久久精品一区二区三区| 精品久久久久久久人妻蜜臀av| 日韩欧美在线乱码| 亚洲无线观看免费| 日本熟妇午夜| 国产白丝娇喘喷水9色精品| 午夜两性在线视频| 精品人妻偷拍中文字幕| 美女大奶头视频| 久久亚洲精品不卡| 精品国产三级普通话版| 国产一区二区在线观看日韩| 丰满的人妻完整版| 草草在线视频免费看| 黄片小视频在线播放| 十八禁网站免费在线| 国产三级黄色录像| 国产三级黄色录像| 一个人免费在线观看的高清视频| 69人妻影院| 女人十人毛片免费观看3o分钟| 亚洲aⅴ乱码一区二区在线播放| 嫩草影院新地址| 欧美成人a在线观看| 成年人黄色毛片网站| 免费搜索国产男女视频| 国产亚洲av嫩草精品影院| 国产欧美日韩一区二区精品| 美女 人体艺术 gogo| 成人亚洲精品av一区二区| 欧美高清成人免费视频www| 亚洲最大成人中文| 国产午夜福利久久久久久| 成人毛片a级毛片在线播放| 成人三级黄色视频| 亚洲激情在线av| 一夜夜www| 一边摸一边抽搐一进一小说| 国产人妻一区二区三区在| 亚洲成a人片在线一区二区| 一级av片app| av在线老鸭窝| 赤兔流量卡办理| 国产精品国产高清国产av| 女同久久另类99精品国产91| 国产男靠女视频免费网站| 国产不卡一卡二| 日韩欧美一区二区三区在线观看| 亚洲熟妇中文字幕五十中出| 亚洲狠狠婷婷综合久久图片| 久久久久久久精品吃奶| 欧美成狂野欧美在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲人成网站在线播| 可以在线观看毛片的网站| 国产黄片美女视频| av福利片在线观看| 一区二区三区高清视频在线| 丁香欧美五月| 亚洲精品色激情综合| 99国产综合亚洲精品| 两个人视频免费观看高清| 可以在线观看的亚洲视频| 校园春色视频在线观看| 男女那种视频在线观看| 日韩高清综合在线| 一进一出抽搐动态| 亚洲avbb在线观看| 特大巨黑吊av在线直播| 国产69精品久久久久777片| 又黄又爽又刺激的免费视频.| 午夜精品一区二区三区免费看| 免费在线观看日本一区| 99热精品在线国产| 亚洲国产精品久久男人天堂| 欧美绝顶高潮抽搐喷水| 999久久久精品免费观看国产| 一级a爱片免费观看的视频| 亚州av有码| 少妇裸体淫交视频免费看高清| 欧美成人一区二区免费高清观看| 一进一出抽搐动态| 久久人妻av系列| 国产精品99久久久久久久久| 亚洲国产高清在线一区二区三| 12—13女人毛片做爰片一| 国产一区二区在线观看日韩| 九色成人免费人妻av| 免费人成视频x8x8入口观看| 免费在线观看成人毛片| 日韩欧美国产在线观看| 色5月婷婷丁香| 久久精品91蜜桃| 一级a爱片免费观看的视频| 日本精品一区二区三区蜜桃| 我要看日韩黄色一级片| 成年人黄色毛片网站| 日本黄色片子视频| 久久久精品大字幕| 午夜两性在线视频| 日韩成人在线观看一区二区三区| www.色视频.com| 国产精品女同一区二区软件 | 色哟哟哟哟哟哟| 欧美黑人欧美精品刺激| 国产精品精品国产色婷婷| 欧美丝袜亚洲另类 | 国产成人av教育| 欧洲精品卡2卡3卡4卡5卡区| 国产av麻豆久久久久久久| 波野结衣二区三区在线| 大型黄色视频在线免费观看| 日本 av在线| 国产精品女同一区二区软件 | 亚洲自偷自拍三级| 岛国在线免费视频观看| 久99久视频精品免费| 日韩欧美国产在线观看| 国产伦在线观看视频一区| 亚洲一区二区三区色噜噜| 国产不卡一卡二| 久久精品综合一区二区三区| 男插女下体视频免费在线播放| 校园春色视频在线观看| 在线播放无遮挡| 久久香蕉精品热| 亚洲成人精品中文字幕电影| 波多野结衣高清作品| av视频在线观看入口| 亚洲人成伊人成综合网2020| 亚洲欧美日韩高清在线视频| 黄色日韩在线| 自拍偷自拍亚洲精品老妇| 成年女人永久免费观看视频| 日韩有码中文字幕| 久久久久精品国产欧美久久久| 人妻制服诱惑在线中文字幕| 嫩草影视91久久| 国产精品一及| 两个人视频免费观看高清| 成熟少妇高潮喷水视频| 18+在线观看网站| 精品国产亚洲在线| 久久人人精品亚洲av| 国产精品日韩av在线免费观看| 一边摸一边抽搐一进一小说| 亚洲精品456在线播放app | 可以在线观看毛片的网站| 日韩精品中文字幕看吧| 国产精品久久久久久亚洲av鲁大| 精品午夜福利在线看| 欧美黑人巨大hd| 性欧美人与动物交配| 天堂av国产一区二区熟女人妻| 日韩欧美在线乱码| av中文乱码字幕在线| 国产亚洲欧美在线一区二区| 久久久久国内视频| 老司机午夜福利在线观看视频| 18+在线观看网站| 999久久久精品免费观看国产| 亚洲成人中文字幕在线播放| 我要看日韩黄色一级片| 国产真实乱freesex| 国内久久婷婷六月综合欲色啪| 国产亚洲精品综合一区在线观看| 欧美日韩综合久久久久久 | 日韩精品青青久久久久久| 久久中文看片网| 亚洲av.av天堂| 国产av不卡久久| 淫秽高清视频在线观看| 精品久久久久久久末码| 成人亚洲精品av一区二区| 看片在线看免费视频| 日韩欧美 国产精品| 亚州av有码| 久久久色成人| 九色国产91popny在线| 日本黄色视频三级网站网址| 久久亚洲精品不卡| 一个人免费在线观看的高清视频| 在现免费观看毛片| 国产在线精品亚洲第一网站| 日本成人三级电影网站| 国产不卡一卡二| 欧美在线一区亚洲| 一区二区三区四区激情视频 | 国产在线男女| 老熟妇仑乱视频hdxx| 国产精品一区二区三区四区久久| av视频在线观看入口| 色综合欧美亚洲国产小说| 特级一级黄色大片| 国产精品久久久久久久电影| 草草在线视频免费看| 精品国产三级普通话版| 国产精品一区二区三区四区免费观看 | 亚洲国产精品999在线| 制服丝袜大香蕉在线| 成年女人毛片免费观看观看9| 熟妇人妻久久中文字幕3abv| 成人三级黄色视频| 制服丝袜大香蕉在线| 色视频www国产| 人人妻人人看人人澡| 91字幕亚洲| 在线免费观看的www视频| 日韩亚洲欧美综合| av专区在线播放| 大型黄色视频在线免费观看| 国产亚洲欧美在线一区二区| 午夜精品久久久久久毛片777| 久久久精品大字幕| 老司机福利观看| 欧美精品啪啪一区二区三区| 97超视频在线观看视频| 日本一二三区视频观看| 老女人水多毛片| 麻豆久久精品国产亚洲av| 久久人人爽人人爽人人片va | 国产精华一区二区三区| 亚洲在线观看片| 一进一出抽搐gif免费好疼| 女生性感内裤真人,穿戴方法视频| 国产精品美女特级片免费视频播放器| 成人欧美大片| 91久久精品电影网| 国产成人福利小说| 色视频www国产| 欧美乱色亚洲激情| 美女免费视频网站| 日日干狠狠操夜夜爽| 人妻久久中文字幕网| 欧美最新免费一区二区三区 | 日韩欧美一区二区三区在线观看| 欧美日韩福利视频一区二区| 人人妻人人看人人澡| 欧美日本亚洲视频在线播放| 日本免费a在线| 成人亚洲精品av一区二区| 嫁个100分男人电影在线观看| 国产成人啪精品午夜网站| 成人永久免费在线观看视频| 国产午夜福利久久久久久| 久久精品影院6| 久久99热6这里只有精品| 精品国内亚洲2022精品成人| 日韩欧美国产一区二区入口| 女生性感内裤真人,穿戴方法视频| 午夜精品久久久久久毛片777| 国产亚洲精品综合一区在线观看| 草草在线视频免费看| 无人区码免费观看不卡| 小说图片视频综合网站| 色哟哟·www| 欧美激情在线99| 国产黄色小视频在线观看| 18禁在线播放成人免费| 最新在线观看一区二区三区| 亚洲国产日韩欧美精品在线观看| 亚洲国产精品成人综合色| 51午夜福利影视在线观看| 男女视频在线观看网站免费| 国产亚洲精品综合一区在线观看| 国内精品美女久久久久久| 免费av毛片视频| 久久99热6这里只有精品| 久久精品91蜜桃| 嫩草影院新地址| 真人做人爱边吃奶动态| 国产欧美日韩一区二区三| 美女高潮的动态| 亚洲经典国产精华液单 | 久久久国产成人免费| 国语自产精品视频在线第100页| 亚洲精品粉嫩美女一区| 国产在视频线在精品| 日韩欧美精品v在线| 国产真实乱freesex| av福利片在线观看| 国产黄色小视频在线观看| 免费在线观看日本一区| 国产不卡一卡二| 无人区码免费观看不卡| 欧美三级亚洲精品| 最近最新免费中文字幕在线| 日韩大尺度精品在线看网址| 久久精品久久久久久噜噜老黄 | 午夜两性在线视频| 长腿黑丝高跟| 一进一出抽搐动态| 亚洲一区二区三区不卡视频| 欧美在线黄色| а√天堂www在线а√下载| 成熟少妇高潮喷水视频| 国产精品一区二区三区四区免费观看 | 欧美最黄视频在线播放免费| 18+在线观看网站| 久久国产乱子免费精品| 日韩av在线大香蕉| 能在线免费观看的黄片| 免费一级毛片在线播放高清视频| 日本成人三级电影网站| 亚洲熟妇中文字幕五十中出| 最后的刺客免费高清国语| 级片在线观看| 亚洲最大成人av| 精品一区二区三区av网在线观看| 深夜精品福利| 中文字幕免费在线视频6| 搡老岳熟女国产| 久久精品夜夜夜夜夜久久蜜豆| 听说在线观看完整版免费高清| 天天躁日日操中文字幕| 综合色av麻豆| 欧美色视频一区免费| 中出人妻视频一区二区| 成人精品一区二区免费| 久久精品人妻少妇| 极品教师在线视频| a级一级毛片免费在线观看| 午夜影院日韩av| 欧美日韩乱码在线| 亚洲不卡免费看| 免费无遮挡裸体视频| 中文字幕人成人乱码亚洲影| 九色成人免费人妻av| 亚洲黑人精品在线| 国产大屁股一区二区在线视频| 神马国产精品三级电影在线观看| 天堂动漫精品| 夜夜爽天天搞| 99久国产av精品| 日韩中文字幕欧美一区二区| 亚洲欧美日韩无卡精品| 精品久久久久久久末码| 日韩大尺度精品在线看网址| 日韩欧美精品v在线| 俄罗斯特黄特色一大片| 99精品久久久久人妻精品| 夜夜躁狠狠躁天天躁| 婷婷精品国产亚洲av| 97超视频在线观看视频| 久久国产精品人妻蜜桃| 欧美区成人在线视频| 亚洲一区二区三区色噜噜| 在线观看av片永久免费下载| 午夜福利免费观看在线| 国产一区二区在线av高清观看| 久久6这里有精品| xxxwww97欧美| h日本视频在线播放| 亚洲国产欧洲综合997久久,| 热99re8久久精品国产| 黄色配什么色好看| 国产精品三级大全| 蜜桃久久精品国产亚洲av| 久久热精品热| 久久午夜福利片| 欧美一级a爱片免费观看看| 婷婷精品国产亚洲av在线| 国产成人影院久久av| 天堂av国产一区二区熟女人妻| 色哟哟哟哟哟哟| 两个人视频免费观看高清| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 俺也久久电影网| 国产 一区 欧美 日韩| 国产精品人妻久久久久久| 少妇被粗大猛烈的视频| 国产伦精品一区二区三区视频9| 久久天躁狠狠躁夜夜2o2o| 人妻丰满熟妇av一区二区三区| av中文乱码字幕在线| 久久久久国内视频| 国产成+人综合+亚洲专区| 日韩欧美三级三区| 成人亚洲精品av一区二区| 午夜两性在线视频| 狠狠狠狠99中文字幕| 99热精品在线国产| 亚洲美女黄片视频| 男人的好看免费观看在线视频| 亚洲欧美精品综合久久99| 国产亚洲精品av在线| 国产免费av片在线观看野外av| 青草久久国产| 亚洲在线自拍视频| 久久精品人妻少妇| 亚洲av免费高清在线观看| 久久久久免费精品人妻一区二区| 日韩欧美三级三区| 蜜桃亚洲精品一区二区三区| 久久热精品热| 亚洲狠狠婷婷综合久久图片| 乱码一卡2卡4卡精品| 2021天堂中文幕一二区在线观| 青草久久国产| 激情在线观看视频在线高清| 噜噜噜噜噜久久久久久91| 亚洲欧美日韩高清在线视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲av第一区精品v没综合| 高潮久久久久久久久久久不卡|