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

    RC框架結(jié)構(gòu)基于構(gòu)件損傷的抗震性能評(píng)估研究*1

    2016-08-25 08:16:54張耀庭杜曉菊
    關(guān)鍵詞:易損性層間試件

    張耀庭,杜曉菊,楊 力

    (1.華中科技大學(xué) 土木工程與力學(xué)學(xué)院,湖北 武漢 430074;2.成都基準(zhǔn)方中建筑設(shè)計(jì)有限公司,四川 成都 610021 )

    ?

    RC框架結(jié)構(gòu)基于構(gòu)件損傷的抗震性能評(píng)估研究*1

    張耀庭?,杜曉菊1,2,楊力1

    (1.華中科技大學(xué) 土木工程與力學(xué)學(xué)院,湖北 武漢430074;2.成都基準(zhǔn)方中建筑設(shè)計(jì)有限公司,四川 成都610021 )

    首先根據(jù)鋼筋混凝土柱的擬靜力試驗(yàn)結(jié)果,探討王東升的修正Park-Ang模型的適用性;然后基于構(gòu)件層次的結(jié)構(gòu)損傷指標(biāo)為地震需求參數(shù),對(duì)鋼筋混凝土框架結(jié)構(gòu)進(jìn)行了基于IDA的易損性分析,并結(jié)合基于最大層間位移角的分析結(jié)果,探討了結(jié)構(gòu)基于損傷的抗震性能評(píng)估方法的可行性.結(jié)果表明:王東升提出的修正Park-Ang模型考慮了加載路徑的影響,能較準(zhǔn)確地反映首次超越破壞后累積損傷的發(fā)展過(guò)程,且總體上判別試件損傷狀態(tài)的準(zhǔn)確性相對(duì)較高;基于損傷的IDA能較好地反映整體結(jié)構(gòu)及局部的損傷發(fā)展過(guò)程、結(jié)構(gòu)的失效破壞機(jī)制,能準(zhǔn)確地判別結(jié)構(gòu)的薄弱環(huán)節(jié),但結(jié)構(gòu)損傷指標(biāo)會(huì)產(chǎn)生超出1.0的情況,不能很好地體現(xiàn)損傷指標(biāo)原始定義的基本含義;相較結(jié)構(gòu)基于最大層間位移角的抗震性能評(píng)估結(jié)果,基于結(jié)構(gòu)損傷的抗震性能評(píng)估方法綜合考慮了結(jié)構(gòu)的響應(yīng)與自身的能力,更全面地評(píng)估了結(jié)構(gòu)的性能水準(zhǔn),能預(yù)測(cè)結(jié)構(gòu)在不同地震強(qiáng)度下各性能狀態(tài)的失效概率.

    鋼筋混凝土框架;修正Park-Ang損傷模型;增量動(dòng)力分析;易損性分析;柱擬靜力試驗(yàn)

    隨著基于性能的抗震設(shè)計(jì)理論與損傷模型研究的發(fā)展,通過(guò)損傷模型來(lái)量化結(jié)構(gòu)性能指標(biāo),實(shí)現(xiàn)基于損傷的抗震性能評(píng)估方法將成為必然.目前,基于鋼筋混凝土結(jié)構(gòu)在地震作用下的兩種破壞模式,國(guó)內(nèi)外學(xué)者采用不同的響應(yīng)參數(shù)建立了各種各樣的損傷模型,包括基于變形、退化、能量、低周疲勞及其組合形式的損傷模型,但沒(méi)有一個(gè)統(tǒng)一、普遍適用的損傷模型.其中,基于變形與能量的Park-Ang損傷模型在地震工程領(lǐng)域中得到較多的應(yīng)用和研究,根據(jù)該模型的缺陷與不足,Kunnath[1]、Kumar[2]、Chai[3]、歐進(jìn)萍[4]、牛荻濤[5]和呂大剛[6]等學(xué)者從不同角度對(duì)其進(jìn)行了修正和改進(jìn).本文首先根據(jù)鋼筋混凝土柱的擬靜力試驗(yàn)結(jié)果,探討王東升的修正Park-Ang模型[7]的適用性.然后,以結(jié)構(gòu)損傷指標(biāo)為地震需求參數(shù),對(duì)鋼筋混凝土框架結(jié)構(gòu)進(jìn)行了基于IDA方法的易損性分析,并結(jié)合基于最大層間位移角的分析結(jié)果,探討了結(jié)構(gòu)基于損傷的抗震性能評(píng)估方法的可行性.

    1 鋼筋混凝土柱的擬靜力試驗(yàn)

    1.1鋼筋混凝土柱的擬靜力試驗(yàn)

    1.1.1試件設(shè)計(jì)

    按照我國(guó)現(xiàn)行規(guī)范設(shè)計(jì)了2根相同的鋼筋混凝土柱試件,分別為C1,C2.柱凈高為1 050 mm,保護(hù)層厚度為20 mm,剪跨比為3.5,設(shè)計(jì)軸壓比為0.15,配筋率為1.7%,配箍率為1.16%,試件詳細(xì)尺寸與配筋信息如圖1所示.混凝土強(qiáng)度等級(jí)為C40,縱筋、箍筋強(qiáng)度等級(jí)分別為HRB400,HPB300,混凝土立方體抗壓強(qiáng)度的實(shí)測(cè)值為48.20 MPa,軸心抗壓強(qiáng)度為36.63 MPa,試驗(yàn)軸壓比為0.076,鋼筋的實(shí)測(cè)性能指標(biāo)見(jiàn)表1.

    1.1.2試驗(yàn)加載裝置及加載制度

    試件采用懸臂自由端加載,如圖2,圖3所示.豎向液壓千斤頂通過(guò)滾軸裝置與反力架橫梁相連,以確保與試件一起平動(dòng),水平推拉液壓千斤頂固定在反力墻上,通過(guò)連接板作用于試件頂部加載中心線上.首先通過(guò)豎向液壓千斤頂一次性施加豎向荷載至預(yù)定值(250 kN),并使其在試驗(yàn)過(guò)程中保持穩(wěn)定,然后通過(guò)水平推拉液壓千斤頂施加水平反復(fù)荷載,并采用位移控制加載.試件C1的位移幅值分別為柱凈高的0.5%,1.0%,1.5%,2.0%,2.5%,3.0%,3.5%和4.0%,每級(jí)荷載循環(huán)兩次;試件C2分三級(jí)等幅值加載,位移幅值分別為柱凈高的2.0%,3.0%和3.5%,每級(jí)荷載循環(huán)10次.當(dāng)試件的承載力下降到最大承載力的85%時(shí),停止試驗(yàn).

    圖1 試件尺寸與配筋Fig.1 Dimension and reinforcement assembly of specimens表1 鋼筋的力學(xué)性能指標(biāo)Tab.1 Mechanical performance indexes of steels

    鋼筋等級(jí)屈服強(qiáng)度/MPa極限強(qiáng)度/MPa彈性模量/105MPaHRB400451.97638.901.98HPB300342.88532.962.08

    圖2 試驗(yàn)加載裝置示意圖Fig.2 Loading equipment diagram of the test

    圖3 試驗(yàn)實(shí)際加載裝置Fig.3 Actual loading equipment of the test

    1.1.3試驗(yàn)現(xiàn)象與結(jié)果

    隨著加載的進(jìn)行,柱底產(chǎn)生肉眼可見(jiàn)的裂紋,并逐漸向上擴(kuò)展、延伸.試件C1加載至第3周時(shí),柱底出現(xiàn)水平貫通的裂紋,柱身側(cè)面產(chǎn)生斜裂紋,試件產(chǎn)生輕微破壞;加載至第8周時(shí),柱底出現(xiàn)起皮掉皮現(xiàn)象,側(cè)面產(chǎn)生多條交叉斜裂縫,裂縫持續(xù)擴(kuò)展,試件產(chǎn)生中等破壞;加載至第14周時(shí),柱底出現(xiàn)多條貫通裂縫且裂縫顯著變寬,柱底保護(hù)層起拱,柱腳混凝土大塊脫落,試件產(chǎn)生嚴(yán)重破壞;加載至第17周時(shí),柱底混凝土大面積壓潰崩落,箍筋外露,縱筋壓屈向外鼓出,試件產(chǎn)生倒塌破壞,試件C1的破壞形態(tài)如圖4所示.

    試件C2加載至第2周時(shí),試件下部多處開(kāi)裂,柱底產(chǎn)生水平貫通的裂紋,側(cè)面產(chǎn)生交叉斜裂紋,反向加載時(shí)裂紋閉合,試件產(chǎn)生輕微破壞;第3~10周循環(huán)加載中,原有裂縫緩慢開(kāi)展,無(wú)其他明顯的破壞現(xiàn)象產(chǎn)生;加載至第11周時(shí),柱底產(chǎn)生豎向裂縫并出現(xiàn)起皮掉皮現(xiàn)象,側(cè)面斜裂縫延伸至柱底,試件產(chǎn)生中等破壞;加載至第15周時(shí),柱底保護(hù)層混凝土脫落,側(cè)面交叉裂縫寬度顯著增大,試件產(chǎn)生嚴(yán)重破壞;加載至第21周時(shí),試件嚴(yán)重傾斜,柱下部混凝土嚴(yán)重脫落,縱筋屈曲,試件產(chǎn)生倒塌破壞,第22~23周加載過(guò)程中,試件強(qiáng)度急劇下降,失去承載能力.試件C2的破壞形態(tài)如圖5所示.試件C1,C2的水平荷載與頂點(diǎn)位移(P-δ)滯回曲線如圖6所示.

    圖4 試件C1的破壞形態(tài)Fig.4 Failure mode of specimen C1

    圖5 試件C2的破壞形態(tài)Fig.5 Failure mode of specimen C2

    δ/mm (a) C1

    δ/mm (b) C2圖6 試件C1,C2的滯回曲線Fig.6 Hysteretic curves of specimens C1 and C2

    1.2Kunnath試驗(yàn)簡(jiǎn)介

    1996年,Kunnath等對(duì)2組共12根1/4比例的圓截面鋼筋混凝土橋墩試件分別進(jìn)行了擬靜力試驗(yàn)和振動(dòng)臺(tái)試驗(yàn)[8].由于本文試驗(yàn)數(shù)量有限,選用PEER數(shù)據(jù)庫(kù)提供的6根試件(A1~A6)的擬靜力試驗(yàn)結(jié)果作為補(bǔ)充.各試件的有效高度均為1 370 mm,截面直徑均為305 mm;縱筋配置21φ9.5 mm,屈服強(qiáng)度為448.0 MPa,配筋率為2.04%;箍筋配置φ4@19 mm,屈服強(qiáng)度為434.0 MPa,配箍率為0.96%,其他信息如表2所示.

    表2 試件A1~A6的基本信息Tab.2 Basic information of specimens A1 to A6

    試件A1~A6的加載制度分別為:

    A1:?jiǎn)握{(diào)加載至失效,其極限位移為150 mm;

    A2:位移幅值分別為1.0%,1.5%,2.0%,2.5%,3.0%,4.0%,5.0%和6.0%,每級(jí)荷載循環(huán)3周,不同位移幅值之間施加1周0.5%的循環(huán)荷載;

    A3~A6:均為等幅值加載,位移幅值分別為2.0%,4.0%,5.5%和7.0%.

    本文僅對(duì)試件A2~A6進(jìn)行分析,其力-位移(P-δ)滯回曲線如圖7所示.

    δ/mm (a) A2

    δ/mm (b) A3

    δ/mm (c) A4

    δ/mm (d) A5

    δ/mm (e) A6圖7 試件A2~A6的滯回曲線Fig.7 Hysteretic curves of specimens A2 to A6

    2 修正Park-Ang損傷模型研究

    2005年,王東升等通過(guò)鋼筋混凝土構(gòu)件的擬靜力試驗(yàn)結(jié)果指出當(dāng)位移延性系數(shù)較小時(shí),Park-Ang損傷模型不能正確反映構(gòu)件極限滯回耗能與位移延性系數(shù)的關(guān)系,不能有效地評(píng)估構(gòu)件的損傷狀態(tài)[7].其在能量項(xiàng)引入加權(quán)因子βi來(lái)考慮加載路徑的影響,修正了該損傷模型,如式(1),(2)所示,各損傷狀態(tài)對(duì)應(yīng)的損傷指標(biāo)范圍如表3所示.

    (1)

    (2)

    式中:Ei為第i個(gè)滯回環(huán)的塑性耗能;βi為加載路徑有關(guān)的加權(quán)因子;β為組合因子;γE為能量等效系數(shù),取γE=0.1;μ0為臨界延性系數(shù),可取為2~3,取μ0=2.5;μi為當(dāng)前構(gòu)件最大位移對(duì)應(yīng)的延性系數(shù);μp為單調(diào)加載下極限位移對(duì)應(yīng)的延性系數(shù).

    表3 不同損傷狀態(tài)對(duì)應(yīng)的損傷指標(biāo)范圍Tab.3 Damage index range of different damage states

    采用該損傷模型評(píng)估試件A2~A6,C1,C2在試驗(yàn)過(guò)程的損傷程度,損傷發(fā)展曲線(T-D)如圖8所示,各試件在加載過(guò)程中的破壞現(xiàn)象及對(duì)應(yīng)損傷指標(biāo)如表4所示.

    表4 試件的破壞現(xiàn)象及損傷指標(biāo)Tab.4 Damage phenomenon and damage index of specimens

    T(a) A2~A6試件

    T(b) C1,C2試件圖8 試件的損傷發(fā)展曲線Fig.8 Damage curve of specimens

    通過(guò)圖8及表3,表4,對(duì)比試驗(yàn)現(xiàn)象及相應(yīng)的損傷指標(biāo),可以看出:

    對(duì)于試件A2,損傷前期(第1周~第15周),損傷增長(zhǎng)緩慢,后期增長(zhǎng)較快,該損傷模型能體現(xiàn)出試件在后期所經(jīng)歷的混凝土明顯脫落、縱筋屈曲、箍筋斷裂等損傷過(guò)程,但其損傷值偏小,試件失效時(shí)其損傷值僅為0.644,低估了試件的破壞程度,與試驗(yàn)現(xiàn)象不吻合;試件A3加載至150周時(shí)僅產(chǎn)生中等破壞,該模型能反映構(gòu)件損傷發(fā)展緩慢的過(guò)程,但其損傷值偏小,低估了構(gòu)件的破壞程度;對(duì)于試件A4~A6,損傷值與試驗(yàn)現(xiàn)象能大致吻合,根據(jù)該損傷模型總體上能判別試件的損傷狀態(tài).

    對(duì)于試件C1,該損傷模型體現(xiàn)出試件在加載過(guò)程中隨位移幅值增加而損傷加重的發(fā)展過(guò)程,能較好地反映構(gòu)件的破壞情況;對(duì)于試件C2,該模型反映了試件在加載前期的損傷演化緩慢的現(xiàn)象,但稍微高估了試件的損傷程度,后期則稍微低估了構(gòu)件的損傷.

    通過(guò)以上分析可知:該損傷模型能較準(zhǔn)確地反映首次超越破壞后,累積損傷的發(fā)展過(guò)程;其考慮了加載路徑的影響,總體上判別試件損傷狀態(tài)的準(zhǔn)確性相對(duì)較高.

    3 RC框架結(jié)構(gòu)基于損傷的易損性分析

    3.1鋼筋混凝土框架結(jié)構(gòu)設(shè)計(jì)

    3.1.1結(jié)構(gòu)基本信息

    按照《建筑抗震設(shè)計(jì)規(guī)范》(GB 5001-2010)和《混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》(GB 50010—2010),采用PKPM設(shè)計(jì)某五層三跨的鋼筋混凝土框架結(jié)構(gòu)辦公樓,底層層高為4.2m,其他各層層高均為3.6m,結(jié)構(gòu)總高度為18.6m,其平面布置圖如圖9所示.抗震設(shè)防烈度為8度(0.2 g),環(huán)境類(lèi)別為Ⅰ類(lèi),場(chǎng)地類(lèi)別為Ⅱ類(lèi),設(shè)計(jì)地震分組為第二組.地表粗糙度為B類(lèi),基本雪壓為0.5 kN/m2,基本風(fēng)壓為0.45 kN/m2.混凝土強(qiáng)度等級(jí)為C40,縱筋強(qiáng)度等級(jí)為HRB400,箍筋強(qiáng)度等級(jí)為HPB300.選取一榀框架進(jìn)行分析,其立面布置圖、配筋圖如圖10所示.

    圖9 框架結(jié)構(gòu)平面布置圖(mm)Fig.9  Frame structural plan configuration(mm)

    圖10 框架立面布置圖及構(gòu)件配筋圖(mm)Fig.10 Frame structural plan configuration and reinforcement assembly of components(mm)

    3.1.2結(jié)構(gòu)有限元模型

    本文采用Opensees軟件對(duì)該鋼筋混凝土框架結(jié)構(gòu)進(jìn)行數(shù)值分析,結(jié)構(gòu)有限元模型采用桿系模型,以梁、柱構(gòu)件為基本單元,結(jié)構(gòu)質(zhì)量集中在節(jié)點(diǎn)上,質(zhì)量根據(jù)重力荷載代表值(恒載+0.5活載)計(jì)算.本文采用纖維截面模型,梁、柱的截面纖維劃分示意圖如圖11示.混凝土本構(gòu)模型采用Concrete01模型,鋼筋本構(gòu)模型采用Steel02模型,在非線性動(dòng)力分析中,各強(qiáng)度參數(shù)取其平均強(qiáng)度值,梁、柱構(gòu)件的混凝土本構(gòu)模型參數(shù)取值見(jiàn)表5.HRB400鋼筋的強(qiáng)度平均值為463.61 MPa,鋼筋的彈性模量為2.0×105MPa,鋼筋從彈性階段過(guò)渡至彈塑性階段的控制參數(shù)取Opensees的建議值,即R0=18.5,cR1=0.93,cR2=0.15,鋼筋的硬化率取0.01.梁、柱構(gòu)件均采用基于力插值函數(shù)的非線性梁柱(Nonlinear Beam Column)單元模型,一根構(gòu)件用一個(gè)單元來(lái)模擬.在進(jìn)行動(dòng)力時(shí)程分析時(shí),采用瑞雷阻尼模型定義結(jié)構(gòu)的阻尼矩陣.

    (a) 柱截面

    (b) 梁截面圖11 框架結(jié)構(gòu)各構(gòu)件截面纖維劃分示意圖Fig.11 Fiber division diagram of component cross section in the frame structure表5 梁、柱構(gòu)件混凝土本構(gòu)模型的參數(shù)取值Tab.5 Parameter values in concrete constitutive model of beam and column

    構(gòu)件抗壓強(qiáng)度/MPa峰值應(yīng)變極限強(qiáng)度/MPa極限應(yīng)變非約束混凝土36.050.00200.0043梁核心混凝土40.040.002219.550.016柱核心混凝土41.130.002326.840.019

    結(jié)構(gòu)有限元模型建立后,通過(guò)模態(tài)分析得到基本周期為0.955 s,與PKPM的結(jié)果0.949 s接近,表明框架結(jié)構(gòu)有限元模型是比較準(zhǔn)確與合理的.

    3.2RC框架結(jié)構(gòu)的損傷指標(biāo)計(jì)算

    對(duì)于鋼筋混凝土框架結(jié)構(gòu),構(gòu)件的損傷主要只集中在塑性鉸區(qū),因此,本文通過(guò)塑性鉸區(qū)截面彎矩-曲率來(lái)計(jì)算截面的損傷指標(biāo),并采用最大損傷指標(biāo)來(lái)表征構(gòu)件的損傷.文獻(xiàn)[7]研究表明該損傷模型的計(jì)算結(jié)果對(duì)β,γE和μo的取值不敏感.因此,對(duì)于梁、柱構(gòu)件,本文統(tǒng)一取μ0=2.5,γE=0.1,β=0.1;梁、柱構(gòu)件的屈服曲率與極限曲率通過(guò)對(duì)結(jié)構(gòu)進(jìn)行Pushover分析獲得.

    本文對(duì)結(jié)構(gòu)進(jìn)行Pushove分析時(shí),采用倒三角加載方式,即水平荷載分布與結(jié)構(gòu)各樓層的水平地震荷載分布一致.因該框架結(jié)構(gòu)梁、柱構(gòu)件分別統(tǒng)一配筋,對(duì)于梁構(gòu)件,其截面屈服曲率與極限曲率均近似相同;對(duì)于柱構(gòu)件,底層柱的軸壓比較大,且地震是往復(fù)作用,底層各邊柱之間、各中柱之間的軸壓比均不同,為方便后續(xù)構(gòu)件損傷指標(biāo)的計(jì)算,區(qū)分為邊柱、中柱:取各邊柱的屈服曲率與極限曲率均與底層邊柱的平均值相同,各中柱的屈服曲率與極限曲率均與底層中柱的平均值相同,這對(duì)于上層柱失效破壞時(shí)的損傷評(píng)估是偏安全的.

    Pushover分析時(shí),為使構(gòu)件達(dá)到極限狀態(tài),取較大的控制位移(最大層間位移角達(dá)到6%),結(jié)構(gòu)基底剪力-最大層間位移角曲線(P-θmax)如圖12所示,各構(gòu)件截面的屈服彎矩、屈服曲率及極限曲率如表6所示.

    θmax/rad圖12 框架結(jié)構(gòu)的Pushover曲線Fig.12 Pushover curve of the frame structure

    3.3增量動(dòng)力分析(IDA)

    影響結(jié)構(gòu)抗震性能的不確定因素主要包括地震

    動(dòng)的不確定性與結(jié)構(gòu)本身的隨機(jī)特性,相對(duì)于地震動(dòng)的不確定性,結(jié)構(gòu)的隨機(jī)性對(duì)結(jié)構(gòu)抗震性能的影響較小.對(duì)于中高層建筑結(jié)構(gòu),當(dāng)選取相對(duì)合理有效的地震強(qiáng)度指標(biāo)時(shí),10~20條地震動(dòng)記錄通常能為結(jié)構(gòu)的抗震性能評(píng)估提供足夠的精度[9].選波長(zhǎng)時(shí)綜合考慮地震動(dòng)的振幅、頻譜特征與持時(shí)要求,地震動(dòng)持續(xù)時(shí)間對(duì)結(jié)構(gòu)的低周疲勞與累積破壞效應(yīng)有重要的影響,一般要求包含地震動(dòng)最強(qiáng)烈部分在內(nèi)的地震動(dòng)持時(shí)為結(jié)構(gòu)基本周期的5~10倍且應(yīng)大于10 s,同時(shí)近場(chǎng)地震動(dòng)與遠(yuǎn)場(chǎng)地震動(dòng)中的響應(yīng)有所不同,為減小結(jié)構(gòu)響應(yīng)的離散性,美國(guó)ATC-63(2008)報(bào)告中通過(guò)震中距應(yīng)大于10 km來(lái)區(qū)分近、遠(yuǎn)場(chǎng)地震動(dòng)記錄.

    表6 梁、柱構(gòu)件截面彎矩-曲率特征值Tab.6 Moment -curvature characteristic value of beam and column cross section

    基于此,本文僅考慮地震動(dòng)的隨機(jī)性,根據(jù)框架結(jié)構(gòu)的場(chǎng)地類(lèi)別,在PEER Strong Motion Database 中選取18條滿足震中距、持時(shí)等要求的地震波對(duì)結(jié)構(gòu)進(jìn)行IDA分析,為結(jié)構(gòu)的地震易損性分析提供研究數(shù)據(jù),具體數(shù)據(jù)如表7所示.各地震波的反應(yīng)譜及均值反應(yīng)譜如圖13所示.

    表7 18條地震波記錄Tab.7 18 seismic records

    T/s圖13  各地震記錄的反應(yīng)譜Fig.13 Response spectrum of seismic records

    基于最大層間位移角θmax的設(shè)計(jì)方法是實(shí)現(xiàn)結(jié)構(gòu)性能化設(shè)計(jì)的主要手段,θmax是運(yùn)用比較廣泛地地震需求參數(shù).因此,本文分別選取基于構(gòu)件層次的結(jié)構(gòu)損傷指數(shù)DE,θmax作為地震需求參數(shù),通過(guò)對(duì)比分析來(lái)探討結(jié)構(gòu)基于損傷的抗震性能評(píng)估方法的可行性.

    《建筑抗震設(shè)計(jì)規(guī)范》(GB 0011-2010)劃分出結(jié)構(gòu)的5個(gè)性能狀態(tài)并給出了相應(yīng)的最大層間位移角θmax參考值,本文結(jié)合規(guī)范參考值及文獻(xiàn)[11],[12]等中的研究成果,確定出各性能狀態(tài)的θmax限值.各性能狀態(tài)下不同地震需求參數(shù)的取值范圍如表8所示.

    表8 各地震需求參數(shù)在不同性能狀態(tài)下取值范圍Tab.8 Value range of seismic demand parameters in different performance states

    在Opensees中可直接得到結(jié)構(gòu)的層間位移角響應(yīng)值,為了減少計(jì)算量以及在基于不同地震需求參數(shù)的IDA曲線上體現(xiàn)出各個(gè)性能點(diǎn),本文取IDA分析的極限狀態(tài)為結(jié)構(gòu)達(dá)到動(dòng)態(tài)不穩(wěn)定、非線性時(shí)程分析不收斂的點(diǎn)或θmax達(dá)到10%與相應(yīng)IDA曲線中20%的初始斜率所對(duì)應(yīng)的地震強(qiáng)度較小點(diǎn).此外,地震強(qiáng)度過(guò)大并無(wú)實(shí)際意義,本文取Sa(T1,5%)的最大值為2.5 gIDA分析時(shí)首次取Sa(T1,5%)=0.05g,此后以0.1 g為增量,依次在[0.1 g,2.5 g] 區(qū)間取值.

    1)基于構(gòu)件損傷的IDA結(jié)果

    采用構(gòu)件塑性鉸區(qū)截面彎矩-曲率計(jì)算整個(gè)結(jié)構(gòu)的損傷指標(biāo)DE,可得到相應(yīng)的IDA曲線(DE-Sa(T1,5%)).以1#地震記錄為例,其IDA曲線如圖14所示.

    DE圖14 1#地震波的IDA曲線Fig.14 IDA curve of 1# seismic record

    從圖14可以看出:當(dāng)Sa(T1,5%)較小時(shí),結(jié)構(gòu)損傷指數(shù)約為0,表明結(jié)構(gòu)處于基本完好狀態(tài);隨著Sa(T1,5%)的增加,IDA曲線逐漸變得平緩,表明結(jié)構(gòu)破壞程度逐步加深;最后結(jié)構(gòu)的損傷指數(shù)遠(yuǎn)遠(yuǎn)大于0.8,表明結(jié)構(gòu)已經(jīng)倒塌失效.

    根據(jù)表8,結(jié)構(gòu)的各性能點(diǎn)如圖15所示,當(dāng)DE=0.8,Sa(T1,5%)=0.508g時(shí),結(jié)構(gòu)達(dá)到倒塌破壞的臨界點(diǎn),此后,損傷指標(biāo)隨地震強(qiáng)度快速增長(zhǎng),能較好地表征結(jié)構(gòu)的失效狀態(tài).

    DE圖15 1#地震波 IDA曲線的性能點(diǎn)Fig.5  Performance point on the IDA curve of 1# seismic record

    1#地震記錄調(diào)幅至不同強(qiáng)度時(shí),結(jié)構(gòu)構(gòu)件塑性鉸區(qū)截面的損傷發(fā)展及各樓層的損傷(DET)分布情況分別如圖16,17所示.由圖16,17可知:隨著地震強(qiáng)度的增加,結(jié)構(gòu)、樓層、構(gòu)件的損傷逐步發(fā)展;樓層損傷由下至上逐漸減小,底層為結(jié)構(gòu)的薄弱層,其損傷發(fā)展較快,柱底破壞嚴(yán)重,Sa(T1,5%)=0.6g時(shí),柱底截面全部失效破壞;結(jié)構(gòu)破壞時(shí)未產(chǎn)生理想的“梁鉸機(jī)制”,而是“梁柱混合機(jī)制”;梁柱節(jié)點(diǎn)處,未完全實(shí)現(xiàn)“強(qiáng)柱弱梁”的設(shè)計(jì)原則.

    注:未標(biāo)注符號(hào)的構(gòu)件表示其處于基本完好狀態(tài)圖16 結(jié)構(gòu)的損傷發(fā)展過(guò)程Fig.16  Damage development process of the structure

    DET圖17 結(jié)構(gòu)各樓層損傷指標(biāo)隨地震強(qiáng)度的變化Fig.17  Damage index variation along with the change of seismic intensity of the structural storeys

    多條地震記錄基于構(gòu)件損傷的IDA曲線(DE-Sa(T1,5%))如圖18所示.不同地震記錄下的IDA曲線與1#地震波的IDA曲線發(fā)展規(guī)律相似,隨著地震強(qiáng)度的增加,IDA曲線逐漸變得平緩;結(jié)構(gòu)達(dá)到IDA分析的極限狀態(tài)時(shí),其損傷值遠(yuǎn)大于其倒塌破壞性能點(diǎn)的限值0.8,結(jié)構(gòu)已經(jīng)失效破壞.

    DE圖18 18條地震波的IDA曲線Fig.18 IDA curve of 18 seismic records

    2)基于最大層間位移角的IDA結(jié)果

    根據(jù)結(jié)構(gòu)的層間位移角響應(yīng),可得到基于最大層間位移角的IDA曲線(θmax-Sa(T1,5%)).以1#地震記錄為例,其IDA曲線如圖19所示.

    從圖19可以看出:當(dāng)Sa(T1,5%)較小時(shí),IDA曲線初始階段近似呈線性增長(zhǎng),表明結(jié)構(gòu)處于彈性狀態(tài);隨著Sa(T1,5%)的增加,曲線斜率逐漸減小,表明結(jié)構(gòu)處于彈塑性發(fā)展?fàn)顟B(tài);最后曲線變得平緩,地震強(qiáng)度僅稍微增大,結(jié)構(gòu)產(chǎn)生較大的變形,表明結(jié)構(gòu)處于倒塌破壞狀態(tài);IDA曲線較好地反映出結(jié)構(gòu)由完好至破壞的發(fā)展過(guò)程.

    θmax/rad圖19 1#地震波的IDA曲線Fig.19 IDA curve of 1# seismic record

    根據(jù)表8,在θmax-Sa(T1,5%)曲線上定義各個(gè)性能點(diǎn),如圖20所示,當(dāng)θmax=1/50,Sa(T1,5%)=0.349g時(shí),結(jié)構(gòu)達(dá)到嚴(yán)重破壞狀態(tài)的極限點(diǎn),這比前文定義的IDA分析的極限狀態(tài)點(diǎn)要保守很多,此時(shí)IDA曲線仍具有較大斜率,結(jié)構(gòu)還可承受更大的地震強(qiáng)度.

    θmax/rad圖20  1#地震波IDA曲線的性能點(diǎn)Fig.20 Performance point on the IDA curve of 1# seismic record

    1#地震波調(diào)幅至不同強(qiáng)度時(shí),結(jié)構(gòu)層間最大位移角的分布情況如圖21所示,結(jié)構(gòu)各層最大層間位移角均隨著地震強(qiáng)度的增大而增大;層間變形隨著樓層增加而減小,在不同地震強(qiáng)度下的最大層間位移角均發(fā)生在第1層,第1層為該結(jié)構(gòu)的薄弱層.

    θmax/rad圖21 結(jié)構(gòu)各層最大層間位移角隨地震強(qiáng)度的變化Fig.21 The maximum inter-story drift ratio variation

    多條地震記錄基于最大層間位移角的 IDA曲線(θmax-Sa(T1,5%))如圖22所示,當(dāng)結(jié)構(gòu)達(dá)到表8中倒塌破壞的性能點(diǎn),即θmax=1/50時(shí),結(jié)構(gòu)仍具有較大的承載力與變形能力.

    θmax/rad圖22 18條地震波的IDA曲線Fig.22 IDA curve of 18 seismic records along with the change of seismic intensity of the structural storeys

    對(duì)比基于不同地震需求參數(shù)的IDA分析結(jié)果可知:

    1)基于結(jié)構(gòu)損傷指標(biāo)的IDA曲線與基于最大層間位移角的IDA曲線發(fā)展規(guī)律相似,基于損傷指標(biāo)的IDA分析能較好地反映結(jié)構(gòu)構(gòu)件、樓層、整體結(jié)構(gòu)的損傷發(fā)展過(guò)程、結(jié)構(gòu)的失效破壞機(jī)制,能判別結(jié)構(gòu)的薄弱構(gòu)件與樓層,從而可以把握結(jié)構(gòu)的局部薄弱環(huán)節(jié),為結(jié)構(gòu)抗震設(shè)計(jì)與分析提供參考依據(jù).基于結(jié)構(gòu)最大層間位移角的IDA分析不能反映結(jié)構(gòu)局部構(gòu)件的損傷情況,僅能粗略地判斷結(jié)構(gòu)的薄弱層;

    2)在同一地震記錄同一強(qiáng)度下,基于不同地震需求參數(shù)所判別的結(jié)構(gòu)性能狀態(tài)可能不同,且結(jié)構(gòu)損傷指標(biāo)會(huì)產(chǎn)生超出1.0的情況,不能很好地體現(xiàn)損傷指標(biāo)原始定義的基本含義,參考規(guī)范給出的最大層間位移角倒塌性能點(diǎn)限值偏保守.

    3.4RC框架結(jié)構(gòu)的地震易損性分析

    3.4.1地震易損性分析模型

    結(jié)構(gòu)的地震易損性分析是指在不同強(qiáng)度地震作用下結(jié)構(gòu)的地震需求超過(guò)某一性能狀態(tài)或某一確定限值C的條件概率.文獻(xiàn)[13]中的研究表明,地震強(qiáng)度IM與結(jié)構(gòu)需求參數(shù)ID之間滿足如下關(guān)系:

    (3)

    其對(duì)數(shù)表達(dá)式為:

    (4)

    式中:α,a,b為回歸系數(shù).

    結(jié)構(gòu)的地震概率需求函數(shù)可用對(duì)數(shù)正態(tài)分布函數(shù)表示,即

    (5)

    (6)

    其對(duì)數(shù)形式為:

    (7)

    (8)

    (9)

    3.4.2基于構(gòu)件損傷的地震易損性分析

    βDE=0.530 2 .

    (10)

    βθmax=0.367 9.

    (11)

    ln(Sa(T1,5%)) (a) 計(jì)算模型ln(PGA)-ln(Sa)回歸曲線

    ln(Sa(T1,5%)) (b) 計(jì)算模型ln(θmax-ln(Sa)回歸曲線圖23 結(jié)構(gòu)概率地震需求模型Fig.23 Probabilistic seismic demand model of the structure

    根據(jù)表8中各個(gè)性能狀態(tài)下DE,θmax的取值范圍,由式(9)可得結(jié)構(gòu)超越各性能水平的概率:

    (12)

    (13)

    相應(yīng)的結(jié)構(gòu)地震易損性曲線如圖24,25所示.

    Sa(T1,5%)/g圖24 結(jié)構(gòu)基于損傷的易損性曲線Fig.24 Structural fragility curve based on damage

    Sa(T1,5%)/g圖25 結(jié)構(gòu)基于最大層間位移角的易損性曲線Fig.25 Structural fragility curve based on the maximum inter-story drift ratio

    由基于不同地震需求參數(shù)的結(jié)構(gòu)易損性分析結(jié)果可以看出:

    1) 各地震記錄對(duì)應(yīng)的ln(Sa(T1,5%)) -ln(DE),ln(Sa(T1,5%))-ln(θmax)近似呈線性關(guān)系,較好地滿足線性要求.

    2) 基于損傷指標(biāo)的易損性曲線與基于最大層間位移角的易損性曲線發(fā)展規(guī)律相似,其能評(píng)估結(jié)構(gòu)在不同地震強(qiáng)度下各性能狀態(tài)的失效概率;結(jié)構(gòu)由基本完好發(fā)展至倒塌性能狀態(tài),基于損傷的易損性曲線逐漸變得平緩,相同地震強(qiáng)度下,其失效概率逐漸減小,曲線發(fā)展趨勢(shì)與結(jié)構(gòu)設(shè)計(jì)原則相吻合.

    3.4.3地震易損性分析結(jié)果的對(duì)比研究

    在上述易損性分析的基礎(chǔ)上,對(duì)比各性能狀態(tài)下基于不同地震需求參數(shù)的分析結(jié)果,如圖26所示.

    Sa(T1,5%)/g (a)基本完好的超越概率

    Sa(T1,5%)/g (b)輕微破壞的超越概率

    Sa(T1,5%)/g (c)中等破壞的超越概率

    Sa(T1,5%)/g (d)嚴(yán)重破壞的超越概率圖26  結(jié)構(gòu)的超越概率曲線Fig.26 Failure probability curve of the structure

    由圖26可知:

    1)地震強(qiáng)度較小時(shí),結(jié)構(gòu)近似完好或損傷程度較小,在不同性能狀態(tài)下基于不同地震需求參數(shù)的易損性曲線基本重合;地震強(qiáng)度較大時(shí),結(jié)構(gòu)彈塑性變形充分發(fā)展,在同一性能狀態(tài)下,結(jié)構(gòu)基于不同地震需求參數(shù)的超越概率不同.

    2)地震強(qiáng)度較大時(shí),在不同性能狀態(tài)下,各易損性曲線均表現(xiàn)出結(jié)構(gòu)基于最大層間位移角的超越概率較大,結(jié)構(gòu)在較小地震強(qiáng)度下倒塌概率已達(dá)到100%,基于結(jié)構(gòu)最大層間位移角的抗震性能評(píng)估結(jié)果偏保守,結(jié)構(gòu)在倒塌性能狀態(tài)下仍有較強(qiáng)的承載能力與變形能力,結(jié)構(gòu)是偏安全的;基于結(jié)構(gòu)損傷的易損性分析結(jié)果則充分考慮了構(gòu)件的耗能能力與非線性變形能力,結(jié)構(gòu)在較大地震強(qiáng)度下才能達(dá)到倒塌性能點(diǎn).

    3) 基于結(jié)構(gòu)變形的單一性能指標(biāo)難以經(jīng)濟(jì)而合理地評(píng)估結(jié)構(gòu)的抗震性能,基于結(jié)構(gòu)損傷的抗震性能評(píng)估方法綜合考慮了結(jié)構(gòu)的響應(yīng)與自身的能力,更全面的評(píng)估了結(jié)構(gòu)的性能,具有良好的可行性.

    4 結(jié) 論

    1)王東升提出的修正Park-Ang模型考慮了加載路徑的影響,能較準(zhǔn)確地反映首次超越破壞后,累積損傷的發(fā)展過(guò)程,且總體上判別試件損傷狀態(tài)的準(zhǔn)確性相對(duì)較高.

    2)基于結(jié)構(gòu)損傷指標(biāo)的IDA曲線與基于最大層間位移角的IDA曲線發(fā)展規(guī)律相似,基于結(jié)構(gòu)損傷的IDA能較好地反映構(gòu)件、樓層、整體結(jié)構(gòu)的損傷發(fā)展過(guò)程及結(jié)構(gòu)的失效破壞機(jī)制,能準(zhǔn)確地判別結(jié)構(gòu)在地震下的局部薄弱環(huán)節(jié),為結(jié)構(gòu)抗震設(shè)計(jì)與分析提供依據(jù),但結(jié)構(gòu)損傷指標(biāo)會(huì)出現(xiàn)超出1.0的情況,不能很好地體現(xiàn)損傷指標(biāo)原始定義的基本含義.

    3)基于結(jié)構(gòu)最大層間位移角的抗震性能評(píng)估結(jié)果偏保守,難以全面、準(zhǔn)確地評(píng)估結(jié)構(gòu)的抗震性能;基于結(jié)構(gòu)損傷的抗震性能評(píng)估方法綜合考慮了結(jié)構(gòu)的響應(yīng)與自身的能力,更全面地評(píng)估了結(jié)構(gòu)的性能水準(zhǔn),能預(yù)測(cè)結(jié)構(gòu)在不同地震強(qiáng)度下各性能狀態(tài)的失效概率,其具有良好的可行性.

    [1]KUNNATH S K, REINHOM A M, PARK Y J. Analytical modeling of inelastic seismic response of R/C structures[J]. Journal of Structural Engineering, 1990, 116(4): 996-1017.

    [2]KUMAR S, USAMI T. Damage evaluation in steel box columns by cyclic loading tests[J]. Journal of Structural Engineering, 1996, 122(6): 626-634.

    [3]CHAI Y H, ROMSTAD K M, BIRD S M. Energy-based linear damage model for high-intensity seismic loading[J]. Journal of Structural Engineering,1995,121(5): 857-864.

    [4]歐進(jìn)萍,牛荻濤,王光遠(yuǎn). 多層非線性抗震鋼結(jié)構(gòu)的模糊動(dòng)力可靠性分析與設(shè)計(jì)[J]. 地震工程與工程振動(dòng),1990, 10(4): 27-37.

    OU Jin-ping, NIU Di-tao, WANG Guang-yuan. Fuzzy dynamical reliability analysis and design of multi-storey nonlinear aseismic steel structures[J].Earthquake Engineering and Engineering Vibration,1990,10(4): 27 -37.(In Chiaese)

    [5]牛荻濤,任利杰. 改進(jìn)的鋼筋混凝土結(jié)構(gòu)雙參數(shù)地震破壞模型[J]. 地震工程與工程振動(dòng),1996, 16(4): 44-54.

    NIU Di-tao, REN Li-jie. A modified seismic damage model with double variables for reinforced concrete structures [J]. Earthquake Engineering and Engineering Vibration, 1996, 16(4): 44-54. (In Chiaese)

    [6]呂大剛,王光遠(yuǎn). 基于損傷性能的抗震結(jié)構(gòu)最優(yōu)設(shè)防水準(zhǔn)的決策方法[J]. 土木工程學(xué)報(bào),2001, 34(1): 44-49.

    LV Dag-ang, WANG Guang-yuan. Decision- making method of optimal fortification level for aseismic structures based ondamage performance[J].China Civil Engineering Journal,2001,34( 1) : 44-49.(In Chiaese)

    [7]王東升,馮啟民,王國(guó)新. 考慮低周疲勞壽命的改進(jìn)Park-Ang地震損傷模型[J]. 土木工程學(xué)報(bào),2004, 37(11): 41-49.

    WANG Dong-sheng, FENG Qi-min, WANG Guo-xin. A modified Park-Ang seismic damage model considering low-cycle fatigue life [J]. China Civil Engineering Journal, 2004, 37(11): 41-49. (In Chiaese)

    [8]KUNNATH S K, El-BAHY A, TAYLOR A W,etal. Cumulative seismic damage of reinforced concrete bridge piers[R]. Technical Report NCEER, US National Center for Earthquake Engineering Research, Buffalo, 1997-0006.

    [9]陳昉健,易偉建. 近場(chǎng)地震作用下銹蝕鋼筋混凝土橋墩的IDA分析[J]. 湖南大學(xué)學(xué)報(bào):自然科學(xué)版, 2015,42(3):1-8.

    CHEN Fang-jian, YE Wei-jian. Incremental dynamic analysis of corroded reinforced concrete bridge columns subjected to near-field earthquake [J]. Journal of Hunan University :Natural Sciences, 2015,42(3):1-8.(In Chiaese)

    [10]-葉列平,馬千里,繆志偉. 結(jié)構(gòu)抗震分析用地震動(dòng)強(qiáng)度指標(biāo)的研究[J]. 地震工程與工程振動(dòng),2009,29(4): 9-22.

    YE Lie-ping,MA Qian-li,MIAO Zhi-wei. Study on earthquake intensities for seismic analysis of structures[J].Earthquake Engineering and Engineering Vibration,2009,29,(4): 9-22. (In Chiaese)

    [11]靳鑫.基于地震動(dòng)參數(shù)的鋼筋混凝土結(jié)構(gòu)易損性分析[D].太原:太原理工大學(xué)建筑與土木工程學(xué)院,2011:47-79.

    JIN Xin. The vulnerability analysis of the reinforced concrete structure based on the ground motion parameters[D]. Taiyuan:College of Architecture and Civil Engineering, Taiyuan University of Technology, 2011:47-79.(In Chiaese)

    [12]黃悠越. 基于構(gòu)件性能的RC框架結(jié)構(gòu)層間位移角性能指標(biāo)限值研究[D].廣州:華南理工大學(xué)土木與交通學(xué)院,2012: 47-64.

    HUANG You-yue. Research on the story drift limits of RC frame structures based on the components performance [D]. Guangzhou:School of Civil Engineering and Transportation, South China University of Technology,2012: 47-64.(In Chiaese)

    [13]龔思禮.建筑結(jié)構(gòu)設(shè)計(jì)系列手冊(cè)-建筑抗震設(shè)計(jì)手冊(cè)[M].2版.北京:中國(guó)建筑工業(yè)出版社,2003:175-187.

    GONG Si-li. Building structure design series manual- building seismic design manual [M]. 2nd ed. Beijing: China Building Industry Press, 2003:175-187.(In Chiaese)

    Research on Seismic Performance Assessment Based on Component Damage for RC Frame Structure

    ZHANG Yao-ting1?, DU Xiao-ju1,2, YANG Li1

    (1.Huazhong Univ of Science and Technology,School of Civil Engineering and Mechanics, Wuhan,Hubei430074, China; 2.Chengdu JZFZ Architecture Design Co Ltd, Chengdu,Sichuan610021,China)

    This paper study discussed the applicability of Park-Ang model modified by Wang based on the pseudo-static test results of reinforced concrete columns. The structural damage index based on the component level was selected, and both the incremental dynamic analysis and fragility analysis of the reinforced concrete frame were conducted. Additionally, the prediction by the modified method to evaluate the seismic performance in terms of the damage index of the RC frame was compared with the maximum inter-story drift ratio of the fragility analysis result. The analysis results show that the modified Park-Ang model that considers the loading path effect predicts well the cumulative damage process after the first damage, which discriminates the specimen damage status accurately, the damage development process, and the overall and local failure mechanism. Further, the weak links can be accurately distinguished. However, the damage index higher than 1.0 may be caused by using these methods, which cannot accurately reflect the basic meaning of the original definition. Compared with the maximum inter-story drift ratio, the seismic performance evaluation method using the damage index considers the structural response and properties. The modified method can be used to evaluate the structural performance more comprehensively and predict the failure probability of the structures under different earthquake loads.

    reinforced concrete frame; fixed Park-Ang damage model; incremental dynamic analysis; fragility analysis; pseudo-static test of columns

    2015-07-03

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51278218),National Natural Science Foundation of China(51278218)

    張耀庭(1965-),男,湖北紅安人,華中科技大學(xué)教授?通訊聯(lián)系人,E-mail:zyt1965@mail.hust.edu

    1674-2974(2016)05-0009-13

    TU375

    A

    猜你喜歡
    易損性層間試件
    復(fù)材管纖維纏繞角度對(duì)約束混凝土軸壓性能的影響研究
    廣東建材(2022年5期)2022-06-10 08:48:20
    基于超聲檢測(cè)的構(gòu)件層間粘接缺陷識(shí)別方法
    基于IDA的預(yù)應(yīng)力混凝土連續(xù)梁橋易損性分析
    不同參數(shù)對(duì)開(kāi)圓洞鋼板剪力墻抗震性能的影響
    基于層間接觸的鋼橋面鋪裝力學(xué)分析
    上海公路(2018年3期)2018-03-21 05:55:40
    基于PSDM和IDA法的深水隔震橋梁地震易損性分析比較
    結(jié)合多分辨率修正曲率配準(zhǔn)的層間插值
    國(guó)際標(biāo)準(zhǔn)攻堅(jiān)克難“S試件”美玉漸成
    大段合采油井層間干擾主控因素研究
    基于性能的FRP加固RC框架結(jié)構(gòu)地震易損性分析
    亚洲中文日韩欧美视频| 成人av在线播放网站| www日本在线高清视频| 国产成人aa在线观看| 国产黄色小视频在线观看| 国产精品,欧美在线| 国产精品香港三级国产av潘金莲| 午夜视频精品福利| 黄色成人免费大全| 黑人操中国人逼视频| 国内精品久久久久精免费| 欧美大码av| 色综合欧美亚洲国产小说| av片东京热男人的天堂| 国产成人av教育| 久久人妻av系列| 一本综合久久免费| 男人舔女人的私密视频| 欧美日韩中文字幕国产精品一区二区三区| 麻豆av在线久日| 国产单亲对白刺激| 亚洲av中文字字幕乱码综合| 成人国产综合亚洲| 婷婷精品国产亚洲av| 国产伦一二天堂av在线观看| 91在线观看av| 亚洲av成人一区二区三| 国产爱豆传媒在线观看| 99热只有精品国产| 久久久色成人| 日日夜夜操网爽| 国产精品久久久久久久电影 | 黄片小视频在线播放| av视频在线观看入口| xxxwww97欧美| 国产成+人综合+亚洲专区| 日本免费一区二区三区高清不卡| 久久久久久国产a免费观看| 男人的好看免费观看在线视频| 99热这里只有精品一区 | 天堂动漫精品| 成年人黄色毛片网站| 欧美国产日韩亚洲一区| 欧美日韩综合久久久久久 | 成人av在线播放网站| 一本精品99久久精品77| 久久精品91蜜桃| 久久这里只有精品中国| 亚洲精品久久国产高清桃花| 九色国产91popny在线| 亚洲精品一卡2卡三卡4卡5卡| 国产一区在线观看成人免费| 母亲3免费完整高清在线观看| 精品久久蜜臀av无| 曰老女人黄片| 丰满的人妻完整版| 亚洲黑人精品在线| 国产高清videossex| 老熟妇乱子伦视频在线观看| 国产精品久久久人人做人人爽| 最近最新中文字幕大全免费视频| 亚洲avbb在线观看| 久久国产精品影院| 欧美黑人欧美精品刺激| 后天国语完整版免费观看| 亚洲欧美日韩高清在线视频| ponron亚洲| 亚洲精品粉嫩美女一区| 首页视频小说图片口味搜索| 熟妇人妻久久中文字幕3abv| www.自偷自拍.com| 国产黄a三级三级三级人| 精品欧美国产一区二区三| 一本精品99久久精品77| 日韩欧美 国产精品| 亚洲18禁久久av| 欧美成狂野欧美在线观看| 99久久成人亚洲精品观看| av天堂中文字幕网| 美女高潮喷水抽搐中文字幕| 后天国语完整版免费观看| 看免费av毛片| 人人妻人人澡欧美一区二区| 亚洲精品中文字幕一二三四区| 久久草成人影院| 不卡av一区二区三区| 国产精品久久久av美女十八| 欧美黑人巨大hd| 午夜久久久久精精品| 中文亚洲av片在线观看爽| 国产私拍福利视频在线观看| 成人av在线播放网站| 中文字幕av在线有码专区| 男人舔奶头视频| 午夜福利免费观看在线| 久久久久国产一级毛片高清牌| 日韩国内少妇激情av| 美女免费视频网站| 国模一区二区三区四区视频 | 亚洲国产欧美人成| 亚洲片人在线观看| 免费人成视频x8x8入口观看| 国产精品一区二区免费欧美| 嫩草影视91久久| 舔av片在线| 91麻豆av在线| 黑人巨大精品欧美一区二区mp4| 1024香蕉在线观看| 国产美女午夜福利| 国产精品国产高清国产av| 亚洲成人中文字幕在线播放| 麻豆av在线久日| 久久久久久久久免费视频了| 人人妻,人人澡人人爽秒播| 欧美性猛交╳xxx乱大交人| 成年女人毛片免费观看观看9| 国内精品一区二区在线观看| 久久精品国产亚洲av香蕉五月| 成人无遮挡网站| 国产aⅴ精品一区二区三区波| 欧美三级亚洲精品| 琪琪午夜伦伦电影理论片6080| 久久这里只有精品19| 欧美乱码精品一区二区三区| 99久久精品热视频| av女优亚洲男人天堂 | 欧美黑人巨大hd| 九色国产91popny在线| 欧美一区二区精品小视频在线| 欧美一区二区国产精品久久精品| 国产单亲对白刺激| 在线观看66精品国产| 亚洲,欧美精品.| 国产成人aa在线观看| 国产精品99久久久久久久久| 一级a爱片免费观看的视频| 啪啪无遮挡十八禁网站| 啦啦啦韩国在线观看视频| 日韩精品青青久久久久久| 久久中文字幕一级| 狂野欧美激情性xxxx| 国产精品乱码一区二三区的特点| 免费电影在线观看免费观看| 99久久国产精品久久久| 成年版毛片免费区| 国产黄a三级三级三级人| 国产亚洲欧美98| h日本视频在线播放| 午夜免费观看网址| 欧美中文日本在线观看视频| 亚洲电影在线观看av| 中文字幕人妻丝袜一区二区| 美女高潮的动态| 曰老女人黄片| 亚洲av免费在线观看| 操出白浆在线播放| 国产成+人综合+亚洲专区| av欧美777| 久久久久国内视频| 国产激情欧美一区二区| 久久亚洲精品不卡| 亚洲人成伊人成综合网2020| 在线a可以看的网站| 国产成年人精品一区二区| 国产精品美女特级片免费视频播放器 | 黄色 视频免费看| 噜噜噜噜噜久久久久久91| 色综合欧美亚洲国产小说| 欧美日韩综合久久久久久 | 国产久久久一区二区三区| 日本黄色视频三级网站网址| 亚洲精品美女久久久久99蜜臀| 国产美女午夜福利| 一二三四在线观看免费中文在| 亚洲熟妇中文字幕五十中出| 成人鲁丝片一二三区免费| 亚洲av五月六月丁香网| 欧美日韩中文字幕国产精品一区二区三区| 久久久久精品国产欧美久久久| 日本黄色视频三级网站网址| 性色av乱码一区二区三区2| 国产精品永久免费网站| 久久欧美精品欧美久久欧美| 一本久久中文字幕| 精品久久久久久久久久免费视频| 老熟妇乱子伦视频在线观看| 欧美日本亚洲视频在线播放| 99久久精品热视频| 精品国产乱码久久久久久男人| 精品久久久久久久末码| 99热6这里只有精品| 日日夜夜操网爽| 人妻夜夜爽99麻豆av| 久久久成人免费电影| 观看美女的网站| 两人在一起打扑克的视频| 亚洲欧美日韩高清在线视频| 久久香蕉精品热| av欧美777| 成人午夜高清在线视频| 狂野欧美激情性xxxx| 熟妇人妻久久中文字幕3abv| 精品久久久久久,| 舔av片在线| 99热这里只有精品一区 | 99久久无色码亚洲精品果冻| 一级毛片女人18水好多| 美女午夜性视频免费| 欧美乱妇无乱码| 亚洲自偷自拍图片 自拍| 一进一出好大好爽视频| 久久精品国产99精品国产亚洲性色| 国产亚洲av高清不卡| 色综合站精品国产| 成熟少妇高潮喷水视频| 国产精品影院久久| 男插女下体视频免费在线播放| 9191精品国产免费久久| 国产精品久久久久久精品电影| 国产美女午夜福利| 精品国产亚洲在线| 露出奶头的视频| 午夜久久久久精精品| 国产精品影院久久| 俺也久久电影网| 91在线观看av| 国产av不卡久久| 一卡2卡三卡四卡精品乱码亚洲| av在线蜜桃| 69av精品久久久久久| 免费在线观看日本一区| 成人三级黄色视频| 午夜激情欧美在线| 全区人妻精品视频| 欧美在线一区亚洲| 国产成人精品久久二区二区免费| 婷婷六月久久综合丁香| 中文字幕精品亚洲无线码一区| 曰老女人黄片| 国产伦人伦偷精品视频| 99久久99久久久精品蜜桃| 91麻豆精品激情在线观看国产| 九色成人免费人妻av| 可以在线观看的亚洲视频| 亚洲午夜精品一区,二区,三区| 长腿黑丝高跟| 天堂影院成人在线观看| 中文字幕人成人乱码亚洲影| 欧美日韩福利视频一区二区| 午夜两性在线视频| 久久久久久久精品吃奶| 亚洲成av人片在线播放无| 97超视频在线观看视频| 99久久综合精品五月天人人| 国产精品久久久久久亚洲av鲁大| 在线观看日韩欧美| 在线观看免费视频日本深夜| 18禁美女被吸乳视频| 一级毛片女人18水好多| 岛国视频午夜一区免费看| 成人亚洲精品av一区二区| 老汉色∧v一级毛片| 国产成人精品无人区| 一级毛片精品| 黄片小视频在线播放| 91麻豆av在线| 亚洲欧美日韩无卡精品| 久久九九热精品免费| 欧美日韩瑟瑟在线播放| 视频区欧美日本亚洲| 丰满人妻熟妇乱又伦精品不卡| 久久国产乱子伦精品免费另类| 好男人在线观看高清免费视频| 国产精品98久久久久久宅男小说| 91av网站免费观看| 亚洲av片天天在线观看| 国产精品日韩av在线免费观看| 亚洲欧洲精品一区二区精品久久久| 午夜影院日韩av| 成人性生交大片免费视频hd| 黄色 视频免费看| 久久久久国内视频| 久久国产精品影院| 国产精品久久久久久人妻精品电影| 国产成人精品久久二区二区免费| 99在线人妻在线中文字幕| 丰满人妻熟妇乱又伦精品不卡| 村上凉子中文字幕在线| 非洲黑人性xxxx精品又粗又长| 日本黄色视频三级网站网址| 久久久久国内视频| 免费看光身美女| or卡值多少钱| 窝窝影院91人妻| 精品国产美女av久久久久小说| 97超级碰碰碰精品色视频在线观看| 一区二区三区激情视频| 久久天躁狠狠躁夜夜2o2o| 久久精品国产综合久久久| 免费电影在线观看免费观看| 亚洲天堂国产精品一区在线| 婷婷六月久久综合丁香| 亚洲欧美精品综合久久99| 国产三级中文精品| 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区色噜噜| 欧美在线一区亚洲| 曰老女人黄片| 国产成人精品无人区| 在线观看舔阴道视频| 久久精品影院6| 国产人伦9x9x在线观看| 久久精品夜夜夜夜夜久久蜜豆| 欧美色欧美亚洲另类二区| 悠悠久久av| 18禁裸乳无遮挡免费网站照片| 搞女人的毛片| 俄罗斯特黄特色一大片| 真人做人爱边吃奶动态| 亚洲精品美女久久久久99蜜臀| 99久久成人亚洲精品观看| 国产不卡一卡二| 最新在线观看一区二区三区| 免费看a级黄色片| 操出白浆在线播放| 伊人久久大香线蕉亚洲五| 一本一本综合久久| 12—13女人毛片做爰片一| 在线观看午夜福利视频| 亚洲人成网站在线播放欧美日韩| 久久久久久久精品吃奶| 欧美日韩精品网址| 成人午夜高清在线视频| 成年版毛片免费区| 免费人成视频x8x8入口观看| 白带黄色成豆腐渣| 99久久无色码亚洲精品果冻| 国产乱人伦免费视频| 亚洲av第一区精品v没综合| 中文字幕av在线有码专区| 国产午夜精品久久久久久| 久久热在线av| 宅男免费午夜| 最新美女视频免费是黄的| 国产精品久久久久久精品电影| 精品午夜福利视频在线观看一区| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产成人精品二区| 日韩欧美国产一区二区入口| 18禁黄网站禁片午夜丰满| 免费在线观看影片大全网站| 一本一本综合久久| 老汉色∧v一级毛片| 欧美乱妇无乱码| 国产午夜福利久久久久久| 又大又爽又粗| 欧美日韩综合久久久久久 | 久久精品亚洲精品国产色婷小说| 亚洲九九香蕉| 国产亚洲精品一区二区www| 久久精品综合一区二区三区| 黑人操中国人逼视频| 综合色av麻豆| 亚洲午夜理论影院| 一个人免费在线观看电影 | 又大又爽又粗| 日本黄色视频三级网站网址| 精品一区二区三区视频在线 | 88av欧美| 99热这里只有是精品50| 亚洲熟妇中文字幕五十中出| 国产激情偷乱视频一区二区| 中文字幕熟女人妻在线| 特级一级黄色大片| 欧美日韩乱码在线| 岛国在线观看网站| av视频在线观看入口| 国产精品一区二区三区四区免费观看 | 久久人人精品亚洲av| 午夜福利在线观看吧| 亚洲欧美日韩高清在线视频| 国产午夜精品久久久久久| 国产精品野战在线观看| 听说在线观看完整版免费高清| 亚洲av片天天在线观看| 日本免费a在线| 国产亚洲精品av在线| 亚洲狠狠婷婷综合久久图片| 超碰成人久久| 亚洲五月天丁香| 国产精品98久久久久久宅男小说| 美女cb高潮喷水在线观看 | 色视频www国产| 无遮挡黄片免费观看| 欧美3d第一页| 午夜福利欧美成人| 丰满人妻熟妇乱又伦精品不卡| 香蕉av资源在线| 一本精品99久久精品77| 亚洲熟妇熟女久久| 人人妻人人看人人澡| 90打野战视频偷拍视频| 老司机深夜福利视频在线观看| 日本免费a在线| 久久久久国产精品人妻aⅴ院| 男女视频在线观看网站免费| 精品不卡国产一区二区三区| 免费观看的影片在线观看| 日韩中文字幕欧美一区二区| 在线观看日韩欧美| 国产野战对白在线观看| 国产精品亚洲av一区麻豆| 日韩欧美精品v在线| 少妇的丰满在线观看| 一级a爱片免费观看的视频| 国产精品久久电影中文字幕| 久久久久久久久久黄片| 麻豆一二三区av精品| 淫秽高清视频在线观看| 亚洲最大成人中文| 亚洲av熟女| 国产乱人伦免费视频| 一本久久中文字幕| 午夜视频精品福利| 美女高潮的动态| 亚洲激情在线av| 午夜精品一区二区三区免费看| 国产精品九九99| 国产成人精品久久二区二区免费| 搡老妇女老女人老熟妇| 国产精品1区2区在线观看.| 亚洲五月婷婷丁香| 国产高清视频在线播放一区| 国产亚洲欧美98| 搡老岳熟女国产| 99热这里只有精品一区 | 中出人妻视频一区二区| 亚洲真实伦在线观看| 久久久久性生活片| 亚洲自拍偷在线| 少妇的丰满在线观看| 国产精品av视频在线免费观看| 亚洲欧洲精品一区二区精品久久久| 国产精品一区二区三区四区久久| 日日摸夜夜添夜夜添小说| av在线天堂中文字幕| 日韩国内少妇激情av| cao死你这个sao货| 亚洲中文字幕一区二区三区有码在线看 | 国产伦在线观看视频一区| 免费大片18禁| h日本视频在线播放| 俺也久久电影网| ponron亚洲| 亚洲成a人片在线一区二区| 一级a爱片免费观看的视频| 亚洲精品乱码久久久v下载方式 | 久9热在线精品视频| 人妻夜夜爽99麻豆av| 欧美日本视频| 亚洲一区二区三区不卡视频| 五月玫瑰六月丁香| 精品日产1卡2卡| 日本五十路高清| 又大又爽又粗| 国产精品98久久久久久宅男小说| 性色av乱码一区二区三区2| 亚洲av中文字字幕乱码综合| 国产高清视频在线观看网站| 99久久成人亚洲精品观看| 国产精品久久久av美女十八| 18禁黄网站禁片午夜丰满| 久久香蕉精品热| 91久久精品国产一区二区成人 | 久久九九热精品免费| 看黄色毛片网站| 黄片小视频在线播放| 性欧美人与动物交配| 亚洲精品一区av在线观看| av女优亚洲男人天堂 | 久久久久久国产a免费观看| 精品久久久久久久久久久久久| 日韩欧美国产一区二区入口| 免费看十八禁软件| 亚洲午夜精品一区,二区,三区| 国产黄色小视频在线观看| 中文字幕av在线有码专区| 国产伦人伦偷精品视频| 日本黄色视频三级网站网址| 热99re8久久精品国产| 搡老岳熟女国产| 日韩欧美国产在线观看| 琪琪午夜伦伦电影理论片6080| 特级一级黄色大片| 男插女下体视频免费在线播放| 亚洲专区中文字幕在线| 美女扒开内裤让男人捅视频| 国产精品免费一区二区三区在线| 天堂影院成人在线观看| 国产91精品成人一区二区三区| 欧美国产日韩亚洲一区| 日本免费一区二区三区高清不卡| 亚洲五月婷婷丁香| 色老头精品视频在线观看| 国产一区二区在线av高清观看| 伦理电影免费视频| www日本黄色视频网| 亚洲人成电影免费在线| 午夜福利在线观看免费完整高清在 | 欧美av亚洲av综合av国产av| 成人性生交大片免费视频hd| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产精品成人综合色| 亚洲av五月六月丁香网| 国产精品av久久久久免费| 国产精品野战在线观看| 一二三四社区在线视频社区8| 观看免费一级毛片| 美女高潮喷水抽搐中文字幕| 91老司机精品| 久久中文字幕一级| 精品国产超薄肉色丝袜足j| 免费在线观看亚洲国产| 免费av毛片视频| 精品久久久久久久人妻蜜臀av| 午夜亚洲福利在线播放| 亚洲中文字幕日韩| 日本撒尿小便嘘嘘汇集6| 热99在线观看视频| 一级毛片高清免费大全| 亚洲国产精品成人综合色| 日韩欧美 国产精品| 久久久久国产一级毛片高清牌| 黄片小视频在线播放| 亚洲一区二区三区不卡视频| 俺也久久电影网| 人人妻人人澡欧美一区二区| 国产精品久久电影中文字幕| 19禁男女啪啪无遮挡网站| 国产精品乱码一区二三区的特点| 99在线视频只有这里精品首页| 欧美在线黄色| 午夜精品一区二区三区免费看| 免费大片18禁| 身体一侧抽搐| 久久草成人影院| 国产欧美日韩精品亚洲av| 国产亚洲精品一区二区www| www.www免费av| 亚洲 国产 在线| 在线观看日韩欧美| 国产乱人视频| 在线看三级毛片| www日本黄色视频网| 精品久久久久久成人av| 九色成人免费人妻av| 国内精品一区二区在线观看| 成人av在线播放网站| 欧美在线黄色| 真人做人爱边吃奶动态| 母亲3免费完整高清在线观看| 亚洲av成人不卡在线观看播放网| av黄色大香蕉| 99国产综合亚洲精品| 一个人看的www免费观看视频| 真人一进一出gif抽搐免费| 91av网一区二区| 精品国内亚洲2022精品成人| 亚洲国产欧美一区二区综合| 男人舔奶头视频| 中文字幕最新亚洲高清| 最近最新中文字幕大全电影3| 99在线人妻在线中文字幕| 亚洲午夜精品一区,二区,三区| 国产成人av教育| 欧美午夜高清在线| 欧美绝顶高潮抽搐喷水| 亚洲成av人片免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 99精品在免费线老司机午夜| 亚洲五月天丁香| 国产av在哪里看| 天堂网av新在线| 亚洲黑人精品在线| 久久久久亚洲av毛片大全| 亚洲天堂国产精品一区在线| 亚洲熟妇熟女久久| 国产精品亚洲美女久久久| 69av精品久久久久久| 亚洲专区国产一区二区| 18美女黄网站色大片免费观看| 黄色成人免费大全| 高清毛片免费观看视频网站| 日本成人三级电影网站| 亚洲色图 男人天堂 中文字幕| 国产一级毛片七仙女欲春2| 欧美乱妇无乱码| 久久欧美精品欧美久久欧美| 美女高潮的动态| av欧美777| 精品国内亚洲2022精品成人| 九色成人免费人妻av| 国产av麻豆久久久久久久| 成在线人永久免费视频| 9191精品国产免费久久| 搡老妇女老女人老熟妇| 欧美日韩精品网址| 禁无遮挡网站| 国产高潮美女av| 首页视频小说图片口味搜索| 日本 av在线| 国产成人av激情在线播放| 露出奶头的视频| 国产精品亚洲美女久久久| 嫩草影院入口|