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

    地震作用下由落梁引起的曲線梁橋連續(xù)倒塌過程分析

    2016-04-21 06:35:13林心潔
    公路工程 2016年1期
    關鍵詞:落梁系梁梁段

    林心潔, 羅 韌, 孫 鵬

    (南京工業(yè)大學 交通學院, 江蘇 南京 210009)

    ?

    地震作用下由落梁引起的曲線梁橋連續(xù)倒塌過程分析

    林心潔, 羅韌, 孫鵬

    (南京工業(yè)大學 交通學院, 江蘇 南京210009)

    [摘要]近年來,大量連續(xù)梁橋在強烈地震中出現(xiàn)了構(gòu)件破壞甚至發(fā)生了整體倒塌的現(xiàn)象,其中橋梁上部結(jié)構(gòu)的落梁又是引起橋梁發(fā)生連續(xù)倒塌的一個關鍵因素。利用有限元軟件ABAQUS,采用合理的材料和結(jié)構(gòu)模型,對一座小半徑曲線梁橋建立模型并對其從發(fā)生落梁到連續(xù)倒塌的整個過程進行細致地模擬和分析。結(jié)果顯示,落梁的碰撞和橋梁的傾斜使得橋墩接連破壞,最終導致橋梁連續(xù)倒塌,橋墩底部和系梁處易發(fā)生塑性破壞,需要對其進行針對性的加固或抗震設計。

    [關鍵詞]連續(xù)倒塌; 落梁; 地震; 有限元; 曲線梁橋

    0前言

    橋梁的震害幾乎伴隨著每次地震的發(fā)生,近十年中我國地震頻發(fā),如2008年汶川地震和2010年玉樹地震中,公路、橋梁、隧道等交通設施大量破壞和倒塌[1]。橋梁連續(xù)倒塌從局部破壞開始并發(fā)生蔓延,最終導致橋梁大范圍破壞和坍塌。通常認為,由于局部構(gòu)件失效或者梁段相對位移過大從支座脫落導致的落梁是引起橋梁連續(xù)倒塌的最主要原因[2]。

    曲線梁橋?qū)儆诜且?guī)則空間結(jié)構(gòu),在我國西南、西北地區(qū)的強震帶,小半徑曲線梁橋被大量采用,目前國內(nèi)外橋梁抗震規(guī)范僅給出其抗震設計應遵循的一般原則,而對于小半徑曲線梁橋則并未明確給出相應的規(guī)定,橋梁抗震研究工作存在一定的局限性,缺乏系統(tǒng)研究[3]。針對曲線橋梁連續(xù)倒塌問題的研究,目前尚處于初始階段,其破壞過程和機理研究不充分,試驗也較少。本文采用數(shù)值模擬的方法,利用ABAQUS軟件,結(jié)合合理的材料本構(gòu)和結(jié)構(gòu)模型,對一座小半徑曲線梁橋在地震中連續(xù)倒塌的過程進行了模擬。

    1混凝土和鋼筋的材料模型

    在軟件ABAQUS中,有多種適用于不同情況的鋼筋和混凝土本構(gòu)模型[4],保證了數(shù)值結(jié)果的準確性。本文中鋼筋和混凝土采用下列本構(gòu)模型。

    1.1混凝土本構(gòu)模型

    受拉段:

    σ=E0ε;0<ε<εt0

    (1)

    (2)

    受壓段:

    (3)

    (4)

    σ=fcu;ε<εcu

    (5)

    式中:fc0為混凝土軸心抗壓強度,εc0為與混凝土軸心抗壓強度所對應的應變;εcu為混凝土極限壓應變,fcu為與混凝土極限壓應變所對應的混凝土強度。

    圖1 混凝土本構(gòu)圖Figure 1 Constitutive model of concrete

    通過引入一個虛擬參考點R來控制混凝土在卸載與再加載時的剛度退化。損傷參數(shù)dcu的大小決定了R的位置。以下表達式(6)即可求得與R點所對應的應變,其對應的應力為E0εR。dcu的取值范圍如式(7):

    (6)

    (7)

    R點確定后,受壓骨架線以內(nèi)的加卸載規(guī)則為: ①卸載點D與R點之間的連線與橫軸(應變軸)的交點為殘余應變εre,該連線的斜率為損傷剛度,記為dcE0,其中dc為受壓損傷系數(shù),如式(8)。 ②卸載時,先按照初始剛度E0卸載,當卸載至經(jīng)過殘余應變點且剛度為0.5dcE0的直線時,改為按0.5dcE0的卸載剛度卸載。 ③在卸載中途再加載時,按初始剛度E0加載,直至達到卸載點與R點的連線時改為按損傷剛度dcE0加載(圖1(b))。

    (8)

    式中:εcu為極限壓應變;εcmax為最大壓應變。

    1.2鋼筋本構(gòu)模型

    (9)

    (10)

    式中:α為屈服剛度系數(shù);Ei為滯回耗能,第i循環(huán);εi為拉壓應變最值,第i循環(huán);εf為單向加載鋼筋破壞時的應變。

    圖2 循環(huán)荷載下的鋼筋雙線性本構(gòu)Figure 2 Under cyclic loading of reinforced bilinear constitutive

    2支座模型與碰撞模型

    2.1支座模型及參數(shù)設置

    支座是能將橋梁上部結(jié)構(gòu)的反力和變形可靠地傳遞給橋梁下部結(jié)構(gòu)的重要部件,但在橋梁抗震中又是一個比較薄弱的環(huán)節(jié)。本文采用雙線型支座模型來模擬橋梁支座,雙線型模型是支座性能分析中最常用的一種非線性模型(見圖3)。其單向軸力-位移關系為:

    Fb=Ke·ub;ub

    Fb=α·Ke·ub+(1-α)Fy;ub>uy

    (11)

    如果用等效線性模型來對具有雙線型模型特性的支座進行模擬,則等效剛度和等效粘滯阻尼比為:

    (12)

    圖3 雙線型模型示意圖Figure 3 Bilinear model diagram

    在模擬支座時,3個方向的剛度k1、k2、k3采用連接單元。其基本屬性包括平移連接屬性(translational type)和旋轉(zhuǎn)屬性(rotational type)。平移連接屬性影響兩個節(jié)點的平移自由度,還可以影響第1個連接點的旋轉(zhuǎn)自由度;旋轉(zhuǎn)連接屬性只影響兩個連接點的旋轉(zhuǎn)自由度。該兩個基本屬性可以相互組合。如圖4、圖5為有限元中的支座模型,本文模型采用Proj-Cartesian與Align結(jié)合使用。已有的大量研究表明,支座水平方向上的彈簧剛度取值對橋梁主體結(jié)構(gòu)的地震響應影響較大,而豎向和轉(zhuǎn)動方向上的彈簧剛度取值影響則較小[7]。因此,3個轉(zhuǎn)動方向上的剛度可簡單得取為0,豎向的剛度取為k2=1.25×109N/m(近似于無窮大),水平方向上的剛度取值根據(jù)橋梁在荷載作用下的可活動性而定,本文中根據(jù)橋梁的實際荷載情況,取k1=k3=1×105N/m,支座的水平摩擦力系數(shù)取μ=0.3[8],本文按照《公路橋梁盆式支座》選用單向活動支座GPZ5DX±50進行建模,其支座最大豎向承壓力為Rck=5 MN,豎向最大應變?yōu)?.02ε。最大位移偏移為Δl=±50 mm,支座厚度20 cm?;瑒又ё瑒拥姆较騽偠葹?,滑動的摩擦系數(shù)取0.02。

    圖4 平動轉(zhuǎn)動約束圖Figure 4 Translational and rotational constraints

    2.2接觸模型及參數(shù)確定

    Kelvin模型的接觸單元是由一個彈簧單元和一個阻尼器并聯(lián)而成的,該模型考慮了碰撞過程中能量的損耗,如圖6、圖7。

    圖6 Kelvin模型Figure 6 Kelvin model

    圖7 Kelvin模型力-位移關系Figure 7 Force-displacement relations of Kelvin model

    兩相鄰跨橋梁之間所產(chǎn)生的碰撞力可表示為:

    (13)

    其中:kk為彈簧單元的剛度值;Ck為阻尼系數(shù);v為鄰梁的相對速度,u1、u2為相鄰橋跨結(jié)構(gòu)的位移,gp為初始間隙(伸縮縫間距)。

    本文在模擬梁間碰撞時采用基于Kelvin模型的接觸單元法,其中碰撞彈簧剛度是Kelvin模型中一個極為重要的參數(shù),參數(shù)取值的不同會導致不同的分析結(jié)果。在定義碰撞剛度時,根據(jù)夏琪[9]對碰撞剛度修正進行的分析,這里取碰撞剛度為:

    3橋梁結(jié)構(gòu)模型的建立

    本文選取汶川地震過程中發(fā)生倒塌的一座小半徑曲線連續(xù)梁橋(垮塌部分曲率半徑為66 m)作為實例,橋梁上部結(jié)構(gòu)共分為6聯(lián),跨徑組合:4×25 m+5×25 m+1×50 m+3×25 m+5×20 m+2×20 m(本文選取第4聯(lián)、發(fā)生倒塌的第5聯(lián)和第6聯(lián)建立有限元模型)。如圖8,為了方便描述落梁過程,將橋墩編號和主梁段編號,梁B10-B13在第4聯(lián),B13-B17在第5聯(lián),B18、B19在第6聯(lián)。10#和19#墩高12 m,11#~15#墩高30.3 m,16#~18#墩高分別為26、 22、 18 m。橋墩采用雙圓柱形式,墩中部設有系梁,16#墩頂部設有固定支座,其它墩頂部為滑動支座;主梁采用C50混凝土,橋墩為C30混凝土,混凝土配筋采用HRB335鋼筋。

    在建立模型時主要采用多尺度的建模方法,主梁為實體單元;橋墩劃分彈性區(qū)域和塑性區(qū)域,彈性區(qū)域采用梁單元,塑性區(qū)域采用實體單元,系梁采用實體單元(見圖8)。支座采用雙線型模型[10],碰撞接觸基于Kelvin[11]模型并結(jié)合ABAQUS自帶的搜索接觸算法處理。

    圖8 橋梁模型圖Figure 8 The model of bridge

    4結(jié)果分析

    由于縱橫向分量共同存在于地震中,斜彎橋會比直線橋產(chǎn)生更大的位移,往往一個任意方向的地震都會使曲線梁橋伸縮縫處產(chǎn)生縱橫向耦合的移位,本文模型中的小半徑曲線梁橋牛腿處搭接構(gòu)造擱置長度有限且未采取有效的限位裝置,發(fā)生落梁難以避免,進而觸發(fā)了橋梁的連續(xù)倒塌。為了便于分析,沿著橋梁軸線向右定義為橋梁縱向(x)的正方向;重力向上的方向為豎向(y)正方向;水平垂直于橋梁軸線指向紙外定義為橋梁橫向(z)的正方向。沿x、y、z三個方向輸入地震波EL-Centro波。

    4.1連續(xù)倒塌過程的各瞬間結(jié)構(gòu)形態(tài)

    數(shù)值模擬分析結(jié)果顯示:在地震力作用下,該曲線梁橋第五聯(lián)發(fā)生連續(xù)倒塌,從開始落梁到完全倒塌整個過程持續(xù)大約9 s的時間,見圖9,倒塌過程可分為4個階段。

    (a) t=33.5 s,B17牛腿處出現(xiàn)落梁

    (b) t=35 s,B17梁段墜落

    (c) t=35.5 s,B13梁段左端落梁,B17梁段撞擊17#墩底部

    (d) t=37.5 s,B13梁段墜落撞擊14#墩

    (e) t=40 s,第五聯(lián)向左側(cè)傾倒

    (f) t=42.5 s,第五聯(lián)完全倒塌

    ① 階段1:落梁前階段(t=0.0 s~33.5 s)。

    橋梁結(jié)構(gòu)整體與地震波保持著一致的振動趨勢,牛腿左端位移波動圖見圖10,圖中可明顯看出牛腿部位的縱向和橫向相對位移波動較大。

    ② 階段2:B17梁段落梁階段(t=33.5 s~35.5 s)。

    在t=33.5 s,B17牛腿段開始落梁(見圖9(a)),梁左端出現(xiàn)塑性鉸,梁段繞塑性鉸墜落并撞擊17#橋墩底部(t=35.5 s)。17#橋墩在B17梁段墜落過程中底部出現(xiàn)塑性鉸破壞,被撞擊后系梁部位出現(xiàn)塑性鉸,系梁以上橋墩倒向18#墩,系梁以下橋墩折向16#墩。此過程中橋梁其他部分梁段向右傾斜,如圖9(b)、圖9(c)。

    如圖10(b)所示,在t=33.5 s時,牛腿左端和牛腿右端的相對位移增大至0.56 m,大于0.4 m的牛腿搭接長度;同時,豎向相對位移急劇增大至2 m,表明B17梁牛腿段開始落梁。如圖10(c),在t=33~35.5 s這段時間內(nèi),B17梁段牛腿的縱向位移和豎向位移增大至極值,表明梁段B17從落梁到撞擊17#橋墩的過程。如圖10(d)、圖10(f),在B17梁段牛腿落梁階段(t=33.5~35.5 s),14#、15#、16#、17#橋墩和B13梁段左端的縱向位移都在正向增大,表明第五聯(lián)整體向著牛腿的方向傾斜。

    ③ 階段3:B13落梁階段(t=35.5~37.5 s)。

    在t=35.5 s,B13梁段左端縱向位移超過搭接長度,B13梁段開始落梁(圖9(c)、圖9(d))。B13梁段變成由14#墩支撐的懸臂梁,右端由于彎矩較大出現(xiàn)塑性鉸;B13梁段繞塑性鉸墜落并撞擊14#橋墩下部系梁(t=37.5 s),使14#橋墩系梁部分發(fā)生塑性破壞。從B13落梁開始,第五聯(lián)開始反向13#墩傾斜,與前階段傾斜方向相反,這種傾斜趨勢一直持續(xù)下去。

    如圖10(d)所示:在t=35.5 s時,B13梁段左端的縱向位移增大至1 m,大于0.8 m的支座搭接長度;同時,豎向相對位移急劇增大至2.3 m,表明B13梁段開始落梁。如圖10(e),在t=35.5~37.5這段時間內(nèi),B13梁段牛腿的縱向位移增大至極值,是梁段B13從落梁到撞擊14#橋墩的過程。如圖10(f),在B17梁段牛腿落梁(t=35.5 s)后,14#、15#、16#橋墩和17#墩的縱向位移出現(xiàn)反方向增大,表明第五聯(lián)整體反向向著13#墩的方向傾斜。出現(xiàn)這種傾斜方向突變的原因是第五聯(lián)的橋墩高度從13#橋墩到17#橋墩依次減小、質(zhì)量分布上靠近13#墩部分較大的緣故。

    圖10 位移Figure 10 Displacement

    ④ 階段4:第五聯(lián)完全倒塌階段(t=37.5~42.5 s)。在地震和落梁作用下,第五聯(lián)整體向13#橋墩傾倒,橋墩底部逐漸破壞,整體落梁倒塌,如圖9(e)、圖9(f)。

    4.2橋墩內(nèi)力分析

    在倒塌過程中,13#和18#橋墩的內(nèi)力較小,基本沒有發(fā)生破壞;14#~17#橋墩倒塌破壞,在此分析13#~17#橋墩的內(nèi)力變化情況。

    ① 階段1:(t=0~33.5 s)。

    如圖11(a),在落梁前,17#墩底軸力較大,橫向和縱向剪力較小,其他橋墩的內(nèi)力在此階段也較小。

    ② 階段2、階段3(t=33.5~37.5 s)。

    如圖11(a)~圖11(c),t=33.5 s,B17梁段落梁,17#橋墩的軸力和彎矩開始增大。t=35.5 s,17#橋墩遭到落梁撞擊,墩底軸力達到最大值,墩底和中部的彎矩也達到最大值,相應位置出現(xiàn)塑性鉸變形,而其后軸力迅速下降。由于17#橋墩的破壞,軸力不斷減小,導致14#、15#和16#橋墩軸力增大。t=37.5 s左右,14#橋墩受B13落梁撞擊,橋墩系梁部位和墩底彎矩達到最大值,因而出現(xiàn)塑性鉸變形。

    ③ 階段4(t=37.5~42.5 s)。

    在牛腿處發(fā)生落梁后,14#~17#橋墩的墩底和系梁部位軸力增大,這是由于第五聯(lián)主梁失去18#和13#橋墩的支撐,導致14#~17#橋墩的承載變大,底部容易出現(xiàn)受壓破壞。較大的軸力會導致系梁和墩連接處發(fā)生剪切破壞。

    4.3應力

    圖11 橋墩內(nèi)力Figure 11 Internal force of piers

    如圖12所示:這是橋墩各部分應力大小對比。從圖上可以看出橋墩墩底豎向應力最大,位于30~40 MPa之間,超過了橋墩所用C30混凝土的抗壓強度,說明橋墩底部發(fā)生受壓破壞。橋墩和系梁連接部位的豎向剪應力也在10~25 MPa之間,剪應力較大;墩頂應力較小,幾乎未破壞。

    圖12 橋墩應力對比Figure 12 Comparison of bridge piers stress

    4.4損傷

    受壓損傷系數(shù)是ABAQUS中混凝土塑性損傷模型中的一個參數(shù),數(shù)值越接近1,表明混凝土接近完全破壞[12]。圖13是百花大橋第5聯(lián)倒塌的14#至17#橋墩的受壓損傷系數(shù)的變化情況圖。

    圖13 橋墩受壓損傷Figure 13 Compressive damage of piers

    從圖13可以看出: ①在受壓損傷方面,系梁處首先受損,然后是橋墩底部,橋墩頂部在橋梁倒塌后才受損; ②系梁與橋墩連接處的豎向剪應力較大,發(fā)生受剪破壞; ③橋墩底部的豎向壓應力較大,發(fā)生受壓損壞; ④墩頂部在連續(xù)倒塌最后階段才出現(xiàn)損傷,這是由于完全倒塌后的碰撞造成。

    5結(jié)論

    a. 采用的考慮抗拉強度的混凝土塑性損傷模型和循環(huán)荷載下的鋼筋多折線模型能較好地模擬橋梁連續(xù)倒塌過程,與實際倒塌情況較好地吻合。利用雙線型支座模型、Kelvin碰撞模型、合理的多尺度建模方法,能較好地模擬橋梁的地震連續(xù)倒塌過程。

    b. 地震作用下,曲線梁橋的倒塌形式并非從初始落梁一端向另一端依次倒塌;而是主梁一端落梁后,相鄰橋墩被落梁撞擊而發(fā)生破壞,橋梁向落梁方向傾斜,引起另一端落梁后,橋梁逐漸倒塌,靠近中間的橋墩一般最后破壞。

    c. 地震作用下: ①伸縮縫處墩梁相對位移過大是引起上部結(jié)構(gòu)落梁的直接原因;橋梁發(fā)生落梁是導致橋梁發(fā)生地震連續(xù)倒塌的觸發(fā)原因,而橋墩破壞則是橋梁發(fā)生地震連續(xù)倒塌的直接原因; ②橋梁兩端的橋墩受損傷較小,未發(fā)生倒塌;其余中間橋墩發(fā)生倒塌。

    d. 橋墩各部分的損傷是不同步的:橋墩系梁部位和橋墩底部發(fā)生損傷較早,橋墩頂部損傷出現(xiàn)在連續(xù)倒塌最后階段。

    e. 在碰撞力方面: ①較早落梁的主梁梁段與橋墩撞擊更加劇烈,初始撞擊力較大且持續(xù)時間長; ②較晚落梁的梁段與相鄰橋墩撞擊力較小,撞擊頻率高且持續(xù)時間短。

    f. 橋梁在地震連續(xù)倒塌過程中的破壞形式如下: ①連續(xù)梁橋橋墩的底部受壓破壞; ②橋墩底部和系梁部位發(fā)生受彎塑性破壞; ③系梁與橋墩連接處發(fā)生受剪破壞。

    [參考文獻]

    [1]李鴻晶,路鳴,羅韌.汶川地震橋梁震害的特征[J].南京工業(yè)大學學報:自然科學版,2009,31(1):24-29.

    [2]屠冰冰,趙冬.連續(xù)倒塌基本問題研究綜述[J].空間結(jié)構(gòu),2013,19(3):20-25.

    [3]王東升,孫治國,李曉莉,等.汶川大地震曲線梁橋震害及破壞機理分析[J].防災減災工程學報,2010,30(5):572-579.

    [4]李廷秀,張研,蔣林華.混凝土單軸動態(tài)強度理論的探討[J].防災減災工程學報,2013,33(2):206-209.

    [5]薛志成,馬強,楊璐,等.混凝土受壓損傷本構(gòu)模型及其數(shù)值驗證[J].混凝土,2012,271(5):41-43.

    [6]汪訓流,陸新征,葉列平.往復荷載下鋼筋混凝土柱受力性能的數(shù)值模擬[J].工程力學,2007,24(12):76-81.

    [7]羅如登,葉梅新,莫朝慶.橋梁支座水平靜力約束方向抗震中的彈簧剛度取值方法對比研究[J].鐵道科學與工程學報,2008,5(2):23-28.

    [8]JTGD62-2004,公路鋼筋混凝土及預應力混凝土橋涵設計規(guī)范[S].

    [9]夏琪.地震中梁式橋的碰撞模擬分析[D].南京:南京工業(yè)大學,2013.

    [10]孫廣俊,李鴻晶,趙鵬飛,等.非規(guī)則多跨簡支梁橋縱向地震反應及參數(shù)影響分析[J].防災減災工程學報,2013,33(4):442-448.

    [11]石巖,王軍文,李建中,等.橫橋向地震作用下簡支梁橋偏心碰撞反應參數(shù)研究[J].振動與沖擊,2012,31(18):132-146.

    [12]莊茁,由小川,廖劍輝,等.基于ABAQUS的有限元分析和應用[M].北京:清華大學出版社,2009.

    The Analysis of the Process of Progressive CollaTPe of Curved Girder Bridge under the Seismic Action

    LIN Xinjie, LUO Ren, SUN Peng

    (College of Transportation Science and Engineering, Nanjing Tech University, Nanjing, Jiangsu 210009, China)

    [Abstract]The progressive collaTPe of bridges caused by falling beams often occurs in earthquake.With the finite element software ABAQUS,a curved continuous girder bridge model was built by using applicable constitutive model of material,bearing and collision.Then its procedure of the progressive collaTPe is simulated.It is caused by unseating and the lean of bridge.The damage of the bottom and the tie beam of piers is more serious,we should strengthen them in design or reinforce them.

    [Key words]progressive collaTPe; earthquake; falling beam; finite element; curved girder bridge

    [中圖分類號]TU 352.1; U 447

    [文獻標識碼]A

    [文章編號]1674-0610(2016)01-0038-08

    [作者簡介]林心潔(1989-),女,福建武平人,碩士研究生。(Linxinjie1206@126.com)

    [基金項目]國家自然科學基金(51278245);江蘇省屬高校自然科學研究重大項目(10KJA560012)。

    [收稿日期]2014-11-12

    猜你喜歡
    落梁系梁梁段
    中小跨徑橋梁柱式橋墩系梁設計的一些探討
    高速鐵路跨海特大橋預制梁段施工技術
    價值工程(2022年25期)2022-09-01 14:04:08
    主跨245 m連續(xù)組合梁落梁技術研究
    系梁對啞鈴型承臺基礎承載特性的影響研究
    偏心支撐耗能梁段內(nèi)力分析及構(gòu)造措施
    高強鋼組合K形偏心支撐框架抗震性能影響參數(shù)分析(Ⅱ)
    大永高速跨高速公路鋼箱梁頂推施工控制
    大跨度組合梁斜拉橋主梁安裝技術研究
    安徽建筑(2019年5期)2019-06-17 02:13:10
    水化熱和徐變收縮對V形墩系梁應力的影響
    上海公路(2018年1期)2018-06-26 08:37:34
    雙柱式橋墩柱間系梁對其抗震性能的影響
    国产精品熟女久久久久浪| 免费一级毛片在线播放高清视频 | 国产精品偷伦视频观看了| 成人午夜精彩视频在线观看| 啦啦啦啦在线视频资源| 99国产精品99久久久久| 777米奇影视久久| 久9热在线精品视频| 叶爱在线成人免费视频播放| 黄片小视频在线播放| 韩国精品一区二区三区| 欧美日韩精品网址| 91精品伊人久久大香线蕉| 亚洲欧洲日产国产| 亚洲人成电影免费在线| 男女下面插进去视频免费观看| 国产不卡av网站在线观看| 999久久久国产精品视频| av不卡在线播放| 久久国产精品人妻蜜桃| 两人在一起打扑克的视频| 最新的欧美精品一区二区| 国产av精品麻豆| 99国产精品免费福利视频| 国产又爽黄色视频| 老司机靠b影院| 午夜福利影视在线免费观看| 男女国产视频网站| 久久精品久久久久久噜噜老黄| 黄频高清免费视频| 91精品三级在线观看| 久久狼人影院| 晚上一个人看的免费电影| 狂野欧美激情性bbbbbb| 99国产精品99久久久久| 国产人伦9x9x在线观看| 亚洲av男天堂| 国产一区二区三区av在线| 国产亚洲精品久久久久5区| 高清视频免费观看一区二区| 精品国产一区二区三区四区第35| 一区二区三区激情视频| 日本欧美国产在线视频| 99九九在线精品视频| 久久国产精品人妻蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 天天操日日干夜夜撸| 欧美日韩综合久久久久久| a级片在线免费高清观看视频| 国产野战对白在线观看| 国产精品国产三级国产专区5o| 大片免费播放器 马上看| 亚洲情色 制服丝袜| 国产男女内射视频| 国产麻豆69| 国产色视频综合| 亚洲国产精品一区二区三区在线| 女性生殖器流出的白浆| 免费看不卡的av| 男人添女人高潮全过程视频| 午夜免费成人在线视频| 午夜两性在线视频| 亚洲精品日本国产第一区| 日日夜夜操网爽| 国产精品久久久av美女十八| 精品福利观看| 中文字幕精品免费在线观看视频| 日日爽夜夜爽网站| 中国国产av一级| 黑人猛操日本美女一级片| 精品久久久久久电影网| 在线观看免费日韩欧美大片| 黄频高清免费视频| 90打野战视频偷拍视频| 亚洲久久久国产精品| 亚洲成人国产一区在线观看 | 丝袜美腿诱惑在线| 丁香六月天网| 国产在视频线精品| 成人18禁高潮啪啪吃奶动态图| 亚洲欧洲日产国产| 午夜激情久久久久久久| 各种免费的搞黄视频| 男人操女人黄网站| 欧美人与性动交α欧美精品济南到| 丝袜人妻中文字幕| 亚洲av欧美aⅴ国产| 男女免费视频国产| 亚洲九九香蕉| 波野结衣二区三区在线| 国产精品.久久久| 中文字幕人妻丝袜一区二区| a级片在线免费高清观看视频| 两个人免费观看高清视频| 精品人妻一区二区三区麻豆| 中文字幕制服av| 亚洲三区欧美一区| 日韩 亚洲 欧美在线| 在现免费观看毛片| 国产免费一区二区三区四区乱码| 黄片播放在线免费| 欧美久久黑人一区二区| 男女高潮啪啪啪动态图| 亚洲国产欧美在线一区| 菩萨蛮人人尽说江南好唐韦庄| 操出白浆在线播放| 国产一区二区三区av在线| 男女床上黄色一级片免费看| 国产免费又黄又爽又色| 久久国产精品男人的天堂亚洲| 精品人妻一区二区三区麻豆| 十八禁高潮呻吟视频| 韩国精品一区二区三区| 国产深夜福利视频在线观看| 蜜桃国产av成人99| 女人被躁到高潮嗷嗷叫费观| 国产精品久久久久久精品古装| 亚洲av综合色区一区| 国产亚洲av高清不卡| 亚洲五月色婷婷综合| 久久国产精品影院| 一级黄色大片毛片| 一边摸一边做爽爽视频免费| 日韩 亚洲 欧美在线| a级毛片黄视频| 一本久久精品| 一级毛片我不卡| 老司机深夜福利视频在线观看 | 成人三级做爰电影| 十分钟在线观看高清视频www| 欧美精品一区二区免费开放| 欧美亚洲 丝袜 人妻 在线| 一区二区三区精品91| 亚洲熟女精品中文字幕| 国产精品偷伦视频观看了| 国产色视频综合| 波多野结衣一区麻豆| 操出白浆在线播放| 亚洲国产av影院在线观看| 午夜老司机福利片| 国产精品秋霞免费鲁丝片| 女警被强在线播放| 99热国产这里只有精品6| 69精品国产乱码久久久| 亚洲色图 男人天堂 中文字幕| 青草久久国产| 别揉我奶头~嗯~啊~动态视频 | 精品福利永久在线观看| 日韩大码丰满熟妇| 国产成人影院久久av| 国产黄频视频在线观看| 只有这里有精品99| 国产精品国产三级专区第一集| 欧美日韩成人在线一区二区| 亚洲五月婷婷丁香| 久9热在线精品视频| 国产国语露脸激情在线看| 性少妇av在线| 国产片内射在线| 精品熟女少妇八av免费久了| 国产色视频综合| 2018国产大陆天天弄谢| 精品久久蜜臀av无| 亚洲欧洲日产国产| 久久久久久久久免费视频了| 一区二区三区四区激情视频| 日日夜夜操网爽| 精品久久久精品久久久| 国产老妇伦熟女老妇高清| 午夜激情av网站| 亚洲欧美精品综合一区二区三区| 久久人人爽人人片av| 韩国高清视频一区二区三区| 成人国产一区最新在线观看 | 欧美精品一区二区大全| 亚洲精品久久久久久婷婷小说| 日韩大码丰满熟妇| 观看av在线不卡| 少妇人妻久久综合中文| 国产精品人妻久久久影院| 国产精品欧美亚洲77777| 亚洲av电影在线观看一区二区三区| 中文字幕av电影在线播放| 亚洲国产欧美网| www日本在线高清视频| 少妇 在线观看| 国产精品久久久久久精品古装| 精品人妻1区二区| 一区福利在线观看| 十八禁网站网址无遮挡| 久久国产精品大桥未久av| 欧美精品人与动牲交sv欧美| 成人黄色视频免费在线看| 日韩一卡2卡3卡4卡2021年| 亚洲精品乱久久久久久| 国产精品人妻久久久影院| 国产片内射在线| 99国产精品99久久久久| 免费日韩欧美在线观看| 日韩一区二区三区影片| 三上悠亚av全集在线观看| 国产又爽黄色视频| 久久精品国产亚洲av涩爱| 中文字幕另类日韩欧美亚洲嫩草| cao死你这个sao货| 汤姆久久久久久久影院中文字幕| 成人黄色视频免费在线看| 国产成人免费观看mmmm| 中国美女看黄片| 美女高潮到喷水免费观看| 午夜福利免费观看在线| 亚洲欧美精品自产自拍| 国产精品成人在线| 久久久精品国产亚洲av高清涩受| 久久久久久亚洲精品国产蜜桃av| av一本久久久久| netflix在线观看网站| 美女国产高潮福利片在线看| 黄网站色视频无遮挡免费观看| 三上悠亚av全集在线观看| av在线老鸭窝| 久久av网站| 波野结衣二区三区在线| 老鸭窝网址在线观看| 日本一区二区免费在线视频| 一本—道久久a久久精品蜜桃钙片| 多毛熟女@视频| 亚洲精品av麻豆狂野| 色婷婷久久久亚洲欧美| 性色av一级| 精品亚洲成a人片在线观看| 免费看不卡的av| 国产亚洲欧美在线一区二区| 在线看a的网站| 日本午夜av视频| av片东京热男人的天堂| 多毛熟女@视频| 9191精品国产免费久久| 黄网站色视频无遮挡免费观看| 日本五十路高清| 蜜桃在线观看..| 女性生殖器流出的白浆| 日韩人妻精品一区2区三区| 日韩伦理黄色片| 男女边吃奶边做爰视频| 国精品久久久久久国模美| 久久久久网色| 天天躁夜夜躁狠狠躁躁| 国产91精品成人一区二区三区 | 久久人人爽av亚洲精品天堂| 99久久综合免费| 日韩中文字幕欧美一区二区 | 制服诱惑二区| 成人三级做爰电影| 亚洲欧美一区二区三区久久| 日韩中文字幕欧美一区二区 | 99热国产这里只有精品6| 亚洲国产中文字幕在线视频| 亚洲人成网站在线观看播放| 亚洲精品一二三| 搡老乐熟女国产| 爱豆传媒免费全集在线观看| 精品久久久久久电影网| 欧美成狂野欧美在线观看| 久久天躁狠狠躁夜夜2o2o | 一级毛片我不卡| 我的亚洲天堂| 久久青草综合色| 99热国产这里只有精品6| 性少妇av在线| 777米奇影视久久| 亚洲欧洲日产国产| 欧美大码av| 中文字幕高清在线视频| 曰老女人黄片| 丝袜脚勾引网站| 久热爱精品视频在线9| 亚洲国产欧美日韩在线播放| 男人舔女人的私密视频| 极品少妇高潮喷水抽搐| 国产深夜福利视频在线观看| 国产亚洲精品第一综合不卡| 国产又爽黄色视频| 成人午夜精彩视频在线观看| 亚洲精品一二三| 国产亚洲欧美在线一区二区| 免费观看人在逋| 亚洲欧美精品自产自拍| 一级片免费观看大全| 久久99精品国语久久久| 最黄视频免费看| 多毛熟女@视频| 亚洲综合色网址| 久久 成人 亚洲| 午夜两性在线视频| 新久久久久国产一级毛片| 免费看不卡的av| 亚洲精品美女久久久久99蜜臀 | 看免费av毛片| 9热在线视频观看99| 三上悠亚av全集在线观看| 精品视频人人做人人爽| 九草在线视频观看| 天天躁夜夜躁狠狠久久av| 视频区图区小说| 欧美精品av麻豆av| 丰满少妇做爰视频| 999久久久国产精品视频| 国产又色又爽无遮挡免| 精品国产一区二区三区久久久樱花| 久久天躁狠狠躁夜夜2o2o | 午夜视频精品福利| 亚洲国产看品久久| 午夜福利在线免费观看网站| av福利片在线| 人妻人人澡人人爽人人| 精品亚洲成a人片在线观看| av欧美777| 波野结衣二区三区在线| 免费高清在线观看视频在线观看| 日韩大片免费观看网站| 9色porny在线观看| 亚洲国产欧美网| 日韩一区二区三区影片| 成年动漫av网址| 国产成人精品在线电影| 天天躁夜夜躁狠狠久久av| 日韩中文字幕视频在线看片| av国产久精品久网站免费入址| 欧美成人精品欧美一级黄| 男女无遮挡免费网站观看| 人人妻人人澡人人爽人人夜夜| 69精品国产乱码久久久| 搡老乐熟女国产| 精品福利观看| 国产亚洲精品第一综合不卡| 波野结衣二区三区在线| bbb黄色大片| 午夜免费鲁丝| 国产亚洲一区二区精品| 久久亚洲国产成人精品v| 日韩一本色道免费dvd| 国产成人免费观看mmmm| 热re99久久国产66热| 久久久久久久国产电影| 亚洲,欧美精品.| 日韩一本色道免费dvd| 免费久久久久久久精品成人欧美视频| 美国免费a级毛片| 每晚都被弄得嗷嗷叫到高潮| av欧美777| 国产成人精品无人区| av福利片在线| 日日爽夜夜爽网站| 久久精品久久久久久久性| 无遮挡黄片免费观看| 亚洲精品美女久久av网站| 亚洲激情五月婷婷啪啪| 大码成人一级视频| 国产免费一区二区三区四区乱码| 国产爽快片一区二区三区| 午夜福利,免费看| 久久人妻熟女aⅴ| 日韩 亚洲 欧美在线| 欧美日本中文国产一区发布| 欧美精品一区二区大全| 国产一区二区在线观看av| 国产成人一区二区在线| 女人久久www免费人成看片| 国产精品国产三级专区第一集| 激情五月婷婷亚洲| 久热爱精品视频在线9| 午夜免费男女啪啪视频观看| 亚洲精品自拍成人| 欧美少妇被猛烈插入视频| 欧美 亚洲 国产 日韩一| 午夜影院在线不卡| 成在线人永久免费视频| 99热国产这里只有精品6| 青青草视频在线视频观看| 亚洲欧美精品综合一区二区三区| 新久久久久国产一级毛片| 99国产综合亚洲精品| 少妇人妻 视频| 91精品国产国语对白视频| 国产一级毛片在线| 中文字幕人妻丝袜一区二区| 性色av乱码一区二区三区2| 首页视频小说图片口味搜索 | 国产淫语在线视频| 大型av网站在线播放| 色综合欧美亚洲国产小说| 久久久久网色| 国产主播在线观看一区二区 | 丁香六月欧美| 国产在线一区二区三区精| 亚洲av日韩精品久久久久久密 | 婷婷丁香在线五月| 99久久99久久久精品蜜桃| 成年av动漫网址| 久久国产精品大桥未久av| av福利片在线| 在现免费观看毛片| 美女午夜性视频免费| 老司机影院毛片| 免费不卡黄色视频| 国产亚洲午夜精品一区二区久久| 中文乱码字字幕精品一区二区三区| 午夜福利一区二区在线看| 男女边摸边吃奶| 9色porny在线观看| 狠狠精品人妻久久久久久综合| 国产免费现黄频在线看| 最新在线观看一区二区三区 | 少妇粗大呻吟视频| 黄片小视频在线播放| 精品国产乱码久久久久久小说| 啦啦啦视频在线资源免费观看| 乱人伦中国视频| 久久精品人人爽人人爽视色| 国产精品亚洲av一区麻豆| 不卡av一区二区三区| 免费少妇av软件| 天天躁夜夜躁狠狠久久av| 欧美大码av| 欧美亚洲日本最大视频资源| 十八禁网站网址无遮挡| 日本午夜av视频| 亚洲国产欧美网| 久久国产精品男人的天堂亚洲| 亚洲欧美一区二区三区黑人| 亚洲精品乱久久久久久| 国产精品熟女久久久久浪| 国产精品久久久人人做人人爽| 日韩大码丰满熟妇| 亚洲国产最新在线播放| 丁香六月天网| 99国产精品一区二区蜜桃av | 两人在一起打扑克的视频| 国产精品免费视频内射| 亚洲精品美女久久久久99蜜臀 | 亚洲国产成人一精品久久久| 免费高清在线观看日韩| av福利片在线| 欧美国产精品va在线观看不卡| 国产一卡二卡三卡精品| 国产精品国产三级专区第一集| 国产av国产精品国产| 亚洲图色成人| 天天躁夜夜躁狠狠久久av| 我的亚洲天堂| 国产野战对白在线观看| 男女边摸边吃奶| 亚洲中文日韩欧美视频| 男女午夜视频在线观看| 久久久亚洲精品成人影院| 人人澡人人妻人| 午夜影院在线不卡| 亚洲美女黄色视频免费看| 日韩免费高清中文字幕av| 两个人免费观看高清视频| 日本午夜av视频| 色婷婷久久久亚洲欧美| 91麻豆精品激情在线观看国产 | 999久久久国产精品视频| 国产精品九九99| xxx大片免费视频| 久久精品久久久久久久性| 国产高清不卡午夜福利| 一边亲一边摸免费视频| 国产精品一二三区在线看| 国产男女超爽视频在线观看| 夫妻午夜视频| 丰满饥渴人妻一区二区三| 精品国产乱码久久久久久男人| 日韩一卡2卡3卡4卡2021年| 一级毛片女人18水好多 | 麻豆国产av国片精品| 无限看片的www在线观看| 欧美黑人精品巨大| 日本午夜av视频| 亚洲欧美中文字幕日韩二区| 亚洲精品美女久久av网站| 91麻豆av在线| 婷婷成人精品国产| 亚洲欧美精品自产自拍| 久久精品久久精品一区二区三区| 一个人免费看片子| 久久精品成人免费网站| 在线观看免费高清a一片| 两个人免费观看高清视频| 精品一区二区三区四区五区乱码 | 欧美av亚洲av综合av国产av| 免费人妻精品一区二区三区视频| 午夜福利影视在线免费观看| 人妻一区二区av| 99久久精品国产亚洲精品| 天堂俺去俺来也www色官网| 天堂中文最新版在线下载| 青春草亚洲视频在线观看| 丁香六月天网| 观看av在线不卡| 在现免费观看毛片| 看免费成人av毛片| 一级毛片电影观看| 亚洲图色成人| 视频在线观看一区二区三区| 无限看片的www在线观看| 国产在线视频一区二区| 首页视频小说图片口味搜索 | 交换朋友夫妻互换小说| netflix在线观看网站| 黄片小视频在线播放| 午夜激情久久久久久久| 男人舔女人的私密视频| 男女边摸边吃奶| 精品国产国语对白av| 亚洲成人国产一区在线观看 | 乱人伦中国视频| 成人黄色视频免费在线看| 亚洲专区国产一区二区| 亚洲欧美一区二区三区黑人| 美女福利国产在线| 深夜精品福利| 国产精品三级大全| 丁香六月天网| 欧美亚洲 丝袜 人妻 在线| 不卡av一区二区三区| 国产精品一区二区在线观看99| 高潮久久久久久久久久久不卡| 国产淫语在线视频| 国产成人91sexporn| 久久久久久免费高清国产稀缺| 天天躁夜夜躁狠狠久久av| 2021少妇久久久久久久久久久| 亚洲av电影在线进入| 国产精品一区二区免费欧美 | 午夜免费成人在线视频| 欧美日韩视频精品一区| 一区二区三区激情视频| 视频在线观看一区二区三区| 无限看片的www在线观看| 国产精品一区二区在线观看99| 一区二区三区四区激情视频| 亚洲国产精品999| 国产激情久久老熟女| 男女午夜视频在线观看| 人人妻人人澡人人看| 悠悠久久av| 欧美精品一区二区免费开放| 一级黄色大片毛片| 51午夜福利影视在线观看| 啦啦啦在线免费观看视频4| 久久精品熟女亚洲av麻豆精品| 精品久久蜜臀av无| 一个人免费看片子| 精品久久久久久电影网| 嫩草影视91久久| 国产极品粉嫩免费观看在线| 老汉色∧v一级毛片| 免费高清在线观看视频在线观看| 少妇 在线观看| 久久久亚洲精品成人影院| 一级a爱视频在线免费观看| 亚洲av日韩在线播放| 亚洲精品日本国产第一区| 中文字幕人妻丝袜制服| 好男人视频免费观看在线| 成人亚洲精品一区在线观看| 久久 成人 亚洲| 黄色视频不卡| av国产精品久久久久影院| 黄片小视频在线播放| 久久久精品94久久精品| 少妇被粗大的猛进出69影院| 可以免费在线观看a视频的电影网站| 亚洲伊人色综图| 视频区图区小说| 少妇粗大呻吟视频| 欧美日韩成人在线一区二区| 少妇 在线观看| 久久天躁狠狠躁夜夜2o2o | 午夜福利免费观看在线| 亚洲,一卡二卡三卡| 国产一区二区激情短视频 | 午夜av观看不卡| 亚洲精品国产区一区二| 国产黄频视频在线观看| 另类精品久久| 日韩大码丰满熟妇| 日韩一区二区三区影片| 91国产中文字幕| 超碰97精品在线观看| 成年人黄色毛片网站| 妹子高潮喷水视频| 一区二区三区精品91| 精品国产超薄肉色丝袜足j| 国产精品国产三级专区第一集| 一二三四在线观看免费中文在| 日韩一本色道免费dvd| 久久久久国产精品人妻一区二区| a级毛片在线看网站| 精品福利观看| 久久天躁狠狠躁夜夜2o2o | 亚洲欧美精品自产自拍| 亚洲精品日本国产第一区| cao死你这个sao货| 女性被躁到高潮视频| 性色av乱码一区二区三区2| 欧美乱码精品一区二区三区| 国产在线免费精品| 9191精品国产免费久久| 国产伦人伦偷精品视频| 高清欧美精品videossex| cao死你这个sao货|