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

    考慮殘余變形的震損框架結(jié)構(gòu)性能評(píng)估及維修優(yōu)化策略

    2022-04-04 16:33:29何浩祥程時(shí)濤廖李燦
    振動(dòng)工程學(xué)報(bào) 2022年1期

    何浩祥 程時(shí)濤 廖李燦

    摘要:針對(duì)實(shí)際震損結(jié)構(gòu)可測(cè)信息有限、修復(fù)加固過程難以模擬以及修復(fù)加固后結(jié)構(gòu)性能不易評(píng)估的現(xiàn)狀,建立梁柱構(gòu)件端部曲率與頂部位移的關(guān)系,厘清殘余位移、位移、曲率、損傷指數(shù)之問的關(guān)聯(lián),提出在有限元中模擬構(gòu)件修復(fù)加固的方法和流程以及根據(jù)不同損傷物理量對(duì)構(gòu)件損傷進(jìn)行量化的方法。根據(jù)不同的損傷物理量以及修復(fù)策略對(duì)鋼筋混凝土震損結(jié)構(gòu)進(jìn)行快速修復(fù),對(duì)修復(fù)后的結(jié)構(gòu)抗震性能進(jìn)行評(píng)估并進(jìn)行地震易損性分析。結(jié)果表明:對(duì)震損結(jié)構(gòu)進(jìn)行快速修復(fù)能夠有效降低結(jié)構(gòu)再遇地震時(shí)發(fā)生嚴(yán)重破壞和倒塌的風(fēng)險(xiǎn);在不同的損傷物理量中根據(jù)殘余位移角量化構(gòu)件損傷,并指導(dǎo)震損結(jié)構(gòu)快速修復(fù)能夠最大限度提高結(jié)構(gòu)抗震性能;相對(duì)于擇件修復(fù)策略,擇層修復(fù)具有更好的效果。

    關(guān)鍵詞:震損結(jié)構(gòu);損傷評(píng)估;修復(fù)加固;殘余變形;易損性

    中圖分類號(hào):TU312+.3;TU375.4

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):10044523(2022)01-0023-11

    DOI: 10.1638 5/j .cnki.issn.10044523.2022.01.003

    引 言

    在地震災(zāi)害中,建筑結(jié)構(gòu)的破壞和倒塌是造成人員傷亡和財(cái)產(chǎn)損失的直接原因,因此結(jié)構(gòu)的抗震設(shè)計(jì)、鑒定和性能評(píng)估具有重大意義[1]。震損結(jié)構(gòu)是指在地震中發(fā)生損傷導(dǎo)致性能部分失效甚至倒塌的建筑結(jié)構(gòu)[2]。在震后救援和結(jié)構(gòu)功能恢復(fù)過程中,對(duì)震損結(jié)構(gòu)的力學(xué)性能和損傷程度做出準(zhǔn)確評(píng)估,并據(jù)此進(jìn)行快速修復(fù),能夠有效提升結(jié)構(gòu)的抗震性能,避免結(jié)構(gòu)在余震中發(fā)生倒塌而造成更嚴(yán)重的人員傷亡和經(jīng)濟(jì)損失。此外,震損結(jié)構(gòu)的快速評(píng)估和修復(fù)優(yōu)化方案也是迅速恢復(fù)災(zāi)區(qū)結(jié)構(gòu)使用功能,對(duì)災(zāi)民進(jìn)行及時(shí)妥善安置的必要技術(shù),成為區(qū)域抗震能力恢復(fù)中至關(guān)重要的環(huán)節(jié)。

    良好的震損結(jié)構(gòu)損傷性能評(píng)估方法不僅可以對(duì)構(gòu)件和結(jié)構(gòu)的損傷程度進(jìn)行量化分析,還可以為結(jié)構(gòu)修復(fù)加固方案和性能提升評(píng)定提供有效的指導(dǎo)和優(yōu)化策略。目前對(duì)構(gòu)件或結(jié)構(gòu)的損傷程度進(jìn)行評(píng)估除采用常用的層間位移角、頂點(diǎn)位移角、剛度變化系數(shù)、構(gòu)件端部曲率和頻率變化率等指標(biāo)引之外,最具有代表性的是Park等[4]在對(duì)大量試驗(yàn)結(jié)果進(jìn)行分析之后提出的由最大位移及滯回耗能線性組合而成的損傷評(píng)價(jià)指標(biāo),但其不能正確反映構(gòu)件極限耗能隨幅值變化的情況,同時(shí)存在閾值范圍不嚴(yán)格、無法直接應(yīng)用于結(jié)構(gòu)動(dòng)力時(shí)程分析等不足。雖然不同研究者提出了Park-Ang修正模型[5-6],但仍不能完全改變其局限性。提出并應(yīng)用準(zhǔn)確而有效的構(gòu)件和結(jié)構(gòu)損傷指數(shù)是提升結(jié)構(gòu)震害性能評(píng)價(jià)的關(guān)鍵。值得指出的是,上述損傷指標(biāo)中的大部分通常是基于結(jié)構(gòu)有限元模擬或試驗(yàn)來獲取的,具有瞬時(shí)性,測(cè)量難度大,對(duì)于實(shí)際震損結(jié)構(gòu)而言,如不能在地震過程中進(jìn)行監(jiān)測(cè),并不能被準(zhǔn)確有效地獲取,這嚴(yán)重阻礙了實(shí)際震損結(jié)構(gòu)性能評(píng)估的實(shí)現(xiàn)和普及。受制于目前震后檢測(cè)的能力和手段,能夠獲得的震損結(jié)構(gòu)性能指標(biāo)一般僅包括裂縫尺寸和數(shù)量、構(gòu)件和整體殘余變形、構(gòu)件破壞類型和數(shù)量等有限信息。能否及如何利用上述指標(biāo)進(jìn)行損傷評(píng)估,并結(jié)合有限元分析和傳統(tǒng)評(píng)估方法提出切實(shí)可行的加固維修策略和優(yōu)化方案是亟需解決的工程問題。

    另一方面,對(duì)震損結(jié)構(gòu)修復(fù)加固之后的性能進(jìn)行準(zhǔn)確模擬及分析,可為制定維修策略和結(jié)構(gòu)易損性評(píng)估提供有效支持。由于震損結(jié)構(gòu)中的部分構(gòu)件已經(jīng)發(fā)生不同程度的損傷甚至嚴(yán)重破壞,在加固維修模擬過程中,必須對(duì)各個(gè)構(gòu)件的既有損傷信息進(jìn)行有效保留,并對(duì)維修部分的構(gòu)件實(shí)現(xiàn)材料和單元模型的有效更新才能準(zhǔn)確計(jì)算其修復(fù)過程的真實(shí)性能,因此對(duì)震損修復(fù)結(jié)構(gòu)進(jìn)行合理的有限元建模和分析具有十分重要的研究意義。然而,在傳統(tǒng)的結(jié)構(gòu)有限元分析中一般只針對(duì)完好結(jié)構(gòu)進(jìn)行分析,或僅選取單根柱等典型構(gòu)件進(jìn)行局部修復(fù)模擬,極少涉及震損結(jié)構(gòu)整體的彈塑性分析、修復(fù)加固模擬及綜合評(píng)估,這極大阻礙了震損結(jié)構(gòu)的性能模擬和加固維修優(yōu)化策略研究,亟需發(fā)展相關(guān)的有限元分析方法。

    此外,目前的震損結(jié)構(gòu)維修加固研究主要考慮構(gòu)件級(jí)別的新型修復(fù)加固技術(shù)和整體減震維修技術(shù)[7],缺乏從力學(xué)機(jī)制、優(yōu)化方法和經(jīng)濟(jì)效益等方面著手的系統(tǒng)深入的研究方法和成套技術(shù)[8]?;趯?shí)際工程能力和需求,在僅已知諸如構(gòu)件和整體的殘余變形及裂縫特征等可測(cè)的損傷性能參數(shù)之后,如果能夠根據(jù)理論分析方法和工程經(jīng)驗(yàn)進(jìn)一步確定構(gòu)件的最優(yōu)維修順序和技術(shù)方案,則可在有限時(shí)間內(nèi)及時(shí)有效地提升結(jié)構(gòu)抗震性能,為結(jié)構(gòu)和區(qū)域的工程韌性提供有利支持。因此,研究震損結(jié)構(gòu)的維修加固理論依據(jù)、擇優(yōu)策略及技術(shù)實(shí)現(xiàn)方法具有重要的理論和工程意義。

    有鑒于此,本文從實(shí)際震損結(jié)構(gòu)維修策略研究的迫切性出發(fā),通過理論分析厘清了構(gòu)件端部殘余變形與位移、曲率和損傷指數(shù)之間的關(guān)系,提出以殘余變形為主要參數(shù)的震損結(jié)構(gòu)損傷評(píng)定思路和分析方法。針對(duì)目前有限元分析方法的不足,提出基于單元替換和重啟動(dòng)分析來實(shí)現(xiàn)震損結(jié)構(gòu)修復(fù)加固的有限元模擬方法。在此基礎(chǔ)上,考慮地震動(dòng)主余震效應(yīng),按不同的損傷評(píng)價(jià)指標(biāo)對(duì)受時(shí)間成本約束的震損結(jié)構(gòu)快速修復(fù)方案進(jìn)行全面細(xì)致的模擬,并對(duì)修復(fù)結(jié)構(gòu)抗震性能進(jìn)行分析,從而實(shí)現(xiàn)對(duì)不同快速修復(fù)方案的比較和擇優(yōu)。通過地震易損性分析獲得了采用不同修復(fù)策略的震損結(jié)構(gòu)加固前后產(chǎn)生損傷的風(fēng)險(xiǎn),為實(shí)際震損結(jié)構(gòu)的修復(fù)和加固提供建議和支撐。

    1 框架結(jié)構(gòu)變形參量與損傷指數(shù)關(guān)聯(lián)性

    1.1 構(gòu)件端部位移、曲率和損傷指數(shù)的關(guān)系

    對(duì)于鋼筋混凝土框架結(jié)構(gòu),其抗震性能主要由鋼筋混凝土柱和梁的性能決定[9]。對(duì)于震損結(jié)構(gòu),震后的柱和梁端部的最大變形和殘余變形是評(píng)價(jià)其損傷性能的主要指標(biāo)。作為構(gòu)件和截面層面的變形量,位移和曲率均能夠直接反映構(gòu)件的材料特性、力學(xué)性能和損傷特征,因此通常采用力一位移或彎矩一曲率兩種力學(xué)參量體系表征鋼筋混凝土柱或梁的抗震性能和滯回能力。雖然構(gòu)件端部的位移或變形較之曲率更直觀且更易測(cè)量,但在理論和有限元分析中也經(jīng)常使用截面曲率來表征損傷,因此需要從理論上建立二者之間的聯(lián)系以便為震損評(píng)定提供統(tǒng)一的標(biāo)準(zhǔn)。在彈塑性狀態(tài)下,構(gòu)件的變形可以近似看作構(gòu)件端部塑性鉸區(qū)域的轉(zhuǎn)動(dòng),通過建立框架結(jié)構(gòu)梁柱構(gòu)件頂部位移與端部曲率的關(guān)系,便可根據(jù)構(gòu)件頂部位移計(jì)算得到相對(duì)應(yīng)的端部曲率。gzslib202204041633

    鋼筋混凝土懸臂柱在軸向力N和水平力P的共同作用下產(chǎn)生變形,如圖1所示。當(dāng)構(gòu)件進(jìn)入彈塑性狀態(tài)時(shí)將產(chǎn)生損傷,柱子以塑性鉸區(qū)域?yàn)橹行陌l(fā)生轉(zhuǎn)動(dòng),其頂部位移“為屈服位移uy和塑性轉(zhuǎn)動(dòng)位移up之和,同時(shí)RC柱的端部曲率φ為屈服曲率φy和塑性曲率φp之和,且可近似認(rèn)為φy沿柱高H線性分布而φp存在于塑性鉸區(qū)域并均勻分布。因此,端部曲率φ和頂部位移u的關(guān)系可表示為:

    此外,Priestley等[11]通過大量試驗(yàn)認(rèn)為塑性鉸長度與構(gòu)件高度H以及縱筋直徑d有關(guān),并提出等效塑性鉸長度表達(dá)式為:

    lp=0.08H+6d

    (4)

    綜上,通過聯(lián)立式(1)~(4),鋼筋混凝土懸臂柱在彈塑性狀態(tài)下的端部曲率φ與頂部位移u之間的關(guān)系為:

    一般情況下,塑性鉸長度lp遠(yuǎn)小于柱高H,因此式(5)可以簡化為:

    框架結(jié)構(gòu)中梁柱構(gòu)件的邊界條件與懸臂柱邊界條件存在一定差別,可利用反彎點(diǎn)法將梁或普通層的柱拆分為兩個(gè)相同的懸臂柱,并對(duì)式(6)進(jìn)行修正,則其端部曲率與頂部位移之間的關(guān)系為:

    底層柱的反彎點(diǎn)一般位于2/3柱高處,可以按同樣方法進(jìn)行修正,不再贅述。由以上結(jié)果可知梁或柱的頂部位移和端部曲率具有明確的線性關(guān)系,力一位移或彎矩一曲率兩種力學(xué)參量體系是相通和內(nèi)洽的,可以互相轉(zhuǎn)換。

    構(gòu)件的耗能能力和抗震性能指數(shù)一般也通過以上兩種力學(xué)參量體系來表征。例如,構(gòu)件延性系數(shù)也可以用曲率延性系數(shù)和位移延性系數(shù)描述,對(duì)于普通鋼混柱,二者近似關(guān)系如下:

    在構(gòu)件地震損傷評(píng)價(jià)中,可采用Park等[4]提出的考慮變形和滯回耗能的雙參數(shù)模型(即Park-Ang模型)進(jìn)行分析,其基本公式如下:式中 D表示構(gòu)件的損傷值;dE表示構(gòu)件的塑性耗能增量;um表示構(gòu)件在地震作用下的實(shí)際位移;us表示構(gòu)件在地震作用下的極限位移;Fy表示構(gòu)件屈服力;β表示權(quán)重系數(shù),其取值范圍一般在0.10 --0.85之間。

    Kunnath等[12]對(duì)Park-Ang模型進(jìn)行改進(jìn),提出基于彎矩一曲率的損傷評(píng)價(jià)模型,公式如下:式中 φm表示構(gòu)件在地震作用下的實(shí)際曲率;φs表示在地震作用下的極限曲率;φy表示在地震作用下的屈服曲率;My表示屈服彎矩;β的含義與Park-Ang模型中相同,表示權(quán)重系數(shù)。

    由上兩式可以看出,基于Park-Ang模型的損傷指數(shù)受位移比和能量比兩項(xiàng)因素的影響。對(duì)大部分構(gòu)件而言第一項(xiàng)的影響較大,損傷指數(shù)與位移之間呈近似線性關(guān)系,但受第二項(xiàng)能量比的影響,最終損傷指數(shù)隨位移變化呈上凹型冪函數(shù)關(guān)系(階次大于1)。

    此外,何浩祥等”J基于彈塑性耗能差的概念認(rèn)為結(jié)構(gòu)的損傷程度可以由結(jié)構(gòu)理想彈性狀態(tài)耗能和實(shí)際耗能的差值與理想彈性耗能的比值來確定。如圖2所示,假定構(gòu)件或結(jié)構(gòu)在彈性狀態(tài)下的剛度為ke,其屈服荷載和屈服位移分別為Fy和uy,當(dāng)進(jìn)入塑性狀態(tài)后,其屈服剛度系數(shù)為a。當(dāng)位移為um時(shí),構(gòu)件或結(jié)構(gòu)所受荷載為Fm,割線剛度為km。此時(shí),構(gòu)件或結(jié)構(gòu)在理想彈性狀態(tài)下和實(shí)際彈塑性狀態(tài)下的耗能差值為三角形ABC的面積。由此,提出單調(diào)力下構(gòu)件或結(jié)構(gòu)的損傷指數(shù)公式如下:

    文獻(xiàn)[14]通過對(duì)梁柱節(jié)點(diǎn)進(jìn)行擬靜力試驗(yàn)驗(yàn)證了基于彈塑性耗能差的損傷評(píng)價(jià)模型的準(zhǔn)確性。由理論和試驗(yàn)結(jié)果可知,基于彈塑性耗能差的損傷指數(shù)與位移之間存在上凸型冪函數(shù)關(guān)系(階次小于1)。此外,Kratzig損傷模型[15-16]和等效黏滯阻尼比的演變趨勢(shì)14均與基于彈塑性耗能差的損傷模型的類似??梢姡M管損傷指數(shù)和結(jié)構(gòu)位移之間具有強(qiáng)相關(guān)性,但由于各損傷模型機(jī)理和內(nèi)涵不同,各損傷指數(shù)與位移的關(guān)系函數(shù)存在一定區(qū)別,在震損評(píng)估時(shí)應(yīng)選擇更嚴(yán)謹(jǐn)準(zhǔn)確的損傷模型。

    1.2 構(gòu)件端部殘余變形與位移的關(guān)系

    值得指出的是,曲率、位移等參數(shù)可以在有限元分析中直接獲得,滯回曲線、延性系數(shù)和損傷指數(shù)等也可通過計(jì)算公式得到,并為結(jié)構(gòu)損傷評(píng)估提供重要依據(jù)。然而對(duì)于真實(shí)震損結(jié)構(gòu)而言,如果不進(jìn)行性能監(jiān)測(cè),在震后僅能得到構(gòu)件裂縫及殘余位移等有限信息,與頂部位移和端部曲率相關(guān)的損傷評(píng)價(jià)指標(biāo)在實(shí)際震損結(jié)構(gòu)性能評(píng)價(jià)中受到了制約??紤]到裂縫產(chǎn)生機(jī)理的復(fù)雜性及隨機(jī)性,建立震后殘余位移和殘余曲率與位移和曲率過程量及損傷指數(shù)的關(guān)系更可行,可為震損結(jié)構(gòu)的損傷程度和性質(zhì)提供直接依據(jù),相應(yīng)的理論和工程意義也更突出。

    結(jié)構(gòu)或構(gòu)件在經(jīng)歷地震作用或往復(fù)加載后,由于剛度退化的緣故在卸載時(shí)會(huì)有不同程度的殘余位移。朱伯龍[17]在三線型滯回曲線模型的基礎(chǔ)上分析了構(gòu)件在近屈服點(diǎn)和近極限點(diǎn)處的滯回環(huán)剛度退化規(guī)律,得到如下式所示的卸載剛度kT的表達(dá)式,并認(rèn)為從屈服滯回環(huán)剛度ky到極限滯回環(huán)剛度ku之間的剛度退化隨位移u的絕對(duì)值增加而發(fā)展。式中 u,uy和uu 分別表示結(jié)構(gòu)或構(gòu)件的位移、屈服位移以及極限位移。

    通過卸載剛度kT,位移u以及與位移u相應(yīng)的荷載P,即可求得從位移u處開始卸載后的殘余位移ut,其計(jì)算公式如下:

    ur=u-P/kT

    (13)

    日本的《公路橋梁抗震規(guī)范》[18]引根據(jù)Takeda模型以及多次試驗(yàn)得到的彎矩一曲率曲線,給出了鋼筋混凝土構(gòu)件在任意時(shí)刻t處的卸載剛度kT,其具體計(jì)算公式如下式所示:

    同時(shí)在阪神地震之后,將對(duì)殘余位移的限制要求首次寫入抗震規(guī)范中,并給出如下式所示的橋墩殘余位移簡化計(jì)算方法:。

    上述兩組公式分別從力一位移、彎矩一曲率兩種力學(xué)體系提出計(jì)算殘余位移的方法,并能夠通過殘余位移反推得到相應(yīng)的位移和曲率。由上文可知構(gòu)件端部曲率與頂部位移具有線性函數(shù)關(guān)系,且可互相換算。因此可根據(jù)殘余位移計(jì)算結(jié)構(gòu)或構(gòu)件位移和曲率,并通過位移獲得損傷指數(shù),從而量化結(jié)構(gòu)損傷程度。

    2 結(jié)構(gòu)震損動(dòng)態(tài)評(píng)估方法

    如前文論述,對(duì)于未采用監(jiān)測(cè)手段的實(shí)際震損結(jié)構(gòu)來說,其損傷評(píng)估一般只能通過觀測(cè)震后表觀裂縫和殘余位移并結(jié)合相關(guān)鑒定和檢測(cè)規(guī)范進(jìn)行定性或定量判斷來實(shí)現(xiàn),如果在震前已獲取結(jié)構(gòu)無損頻率,也可以通過地震前后結(jié)構(gòu)的頻率變化來確定結(jié)構(gòu)的損傷程度。除以上數(shù)據(jù)之外,很難獲取更豐富的信息。通過結(jié)構(gòu)彈塑性有限元模擬對(duì)結(jié)構(gòu)在地震下的損傷演變和特性進(jìn)行仿真既可以預(yù)測(cè)結(jié)構(gòu)的損傷位置和程度,也可以為結(jié)構(gòu)的維修加固提供有效的指導(dǎo)。由于在結(jié)構(gòu)有限元時(shí)程分析中可以提取全面的結(jié)構(gòu)響應(yīng)信息,因此可以運(yùn)用傳統(tǒng)的損傷指數(shù)和分析方法來完成基本的損傷評(píng)估。然而,傳統(tǒng)的損傷指數(shù)存在諸多不足:以最大層間位移角或最大頂點(diǎn)位移角為代表的極值型指標(biāo)并不能反映結(jié)構(gòu)的損傷演變過程,也不能表征具體構(gòu)件的損傷情況;以Park-Ang損傷指數(shù)[4]和Kratzig損傷指數(shù)[15]為代表的組合累積型指標(biāo)限于周期滯回耗能的假設(shè)在動(dòng)力時(shí)程分析中無法有效應(yīng)用,且存在閾值不嚴(yán)格等局限。建立準(zhǔn)確有效的損傷指數(shù)是對(duì)震損結(jié)構(gòu)及修復(fù)結(jié)構(gòu)進(jìn)行損傷評(píng)估的必要手段。gzslib202204041633

    針對(duì)目前動(dòng)態(tài)損傷評(píng)價(jià)方法存在的缺陷,何浩祥等[13]提出基于彈塑性耗能差的結(jié)構(gòu)動(dòng)態(tài)損傷指數(shù),計(jì)算公式如下:

    試驗(yàn)和數(shù)值模擬得到的結(jié)果表明基于彈塑性耗能差的結(jié)構(gòu)損傷指數(shù)具有準(zhǔn)確度高、閾值嚴(yán)格、動(dòng)態(tài)表征、適用于多種構(gòu)件和結(jié)構(gòu)整體等優(yōu)點(diǎn)。根據(jù)崔濟(jì)東等[19]建議的結(jié)構(gòu)不同損傷程度對(duì)應(yīng)的位移角限值,將與位移角限值對(duì)應(yīng)的耗能差損傷指數(shù)劃分到同一損傷等級(jí),由此得到如表1所示的結(jié)構(gòu)不同損傷程度所對(duì)應(yīng)的耗能差損傷指數(shù)取值范圍。

    上述方法能夠反映具體地震動(dòng)下的結(jié)構(gòu)損傷演變過程。結(jié)構(gòu)地震易損性分析[20-21]則能夠展現(xiàn)在一定地震動(dòng)強(qiáng)度和隨機(jī)性下結(jié)構(gòu)地震反應(yīng)超越某種極限狀態(tài)的概率,從而反映在不同地震動(dòng)強(qiáng)度下結(jié)構(gòu)的破壞程度。綜上,對(duì)結(jié)構(gòu)有限元模型進(jìn)行大量時(shí)程分析,并將損傷計(jì)算方法和地震易損性分析結(jié)合才能實(shí)現(xiàn)精準(zhǔn)的結(jié)構(gòu)損傷評(píng)估和預(yù)測(cè)。某一確定地震動(dòng)強(qiáng)度指標(biāo)值下結(jié)構(gòu)性能的超越概率可表達(dá)為:式中 DMlim表示該狀態(tài)下結(jié)構(gòu)響應(yīng)或損傷指數(shù)的限值;DMmd為結(jié)構(gòu)響應(yīng)或損傷指數(shù)中位值;β表示易損性函數(shù)的對(duì)數(shù)標(biāo)準(zhǔn)差。

    綜上,本文建議將基于彈塑性耗能差的結(jié)構(gòu)損傷指數(shù)作為評(píng)估結(jié)構(gòu)震損破壞程度的依據(jù),并應(yīng)用于結(jié)構(gòu)地震易損性分析,以定量反映結(jié)構(gòu)發(fā)生不同程度破壞的超越概率。同時(shí),通過比較震損結(jié)構(gòu)修復(fù)前后以及不同修復(fù)方案下的地震易損性曲線特征來評(píng)價(jià)不同修復(fù)方案的優(yōu)劣,為震損結(jié)構(gòu)的實(shí)際修復(fù)提供指導(dǎo)。

    3 基于有限元模擬震損結(jié)構(gòu)修復(fù)加固流程

    受限于經(jīng)濟(jì)成本,大批量制作框架結(jié)構(gòu)試件并進(jìn)行震損后的加固性能試驗(yàn)較難實(shí)現(xiàn)。為了在有限元中準(zhǔn)確模擬震損結(jié)構(gòu)經(jīng)過部分修復(fù)加固后的抗震性能,首先要在有限元中對(duì)結(jié)構(gòu)先施加地震作用使結(jié)構(gòu)產(chǎn)生損傷成為震損結(jié)構(gòu),隨后對(duì)震損結(jié)構(gòu)中所有構(gòu)件的損傷程度進(jìn)行量化,最后根據(jù)量化結(jié)果對(duì)損傷嚴(yán)重的梁柱構(gòu)件進(jìn)行性能恢復(fù)以模擬修復(fù)加固效果。然而,大部分傳統(tǒng)有限元軟件并沒有模擬該過程的功能。目前利用有限元軟件模擬加固后結(jié)構(gòu)抗震性能的研究較薄弱且主要是針對(duì)構(gòu)件進(jìn)行的,對(duì)于結(jié)構(gòu)整體抗震性能模擬分析的研究較少。針對(duì)這些問題,提出一種在有限元中切實(shí)可行的修復(fù)加固模擬方法。

    本文利用開源有限元軟件OpenSEES來實(shí)現(xiàn)模擬震損結(jié)構(gòu)修復(fù)加固的過程。在OpenSEES軟件中常通過刪除構(gòu)件單元這一命令模擬框架結(jié)構(gòu)中梁、柱構(gòu)件損傷失效,引起結(jié)構(gòu)內(nèi)力重新分布,導(dǎo)致結(jié)構(gòu)連續(xù)倒塌的動(dòng)態(tài)過程。而實(shí)際震損結(jié)構(gòu)中損傷構(gòu)件的修復(fù)過程從一定角度可以理解為是原損傷構(gòu)件拆除和新構(gòu)件建立的過程。因此,為在有限元軟件中實(shí)現(xiàn)震損結(jié)構(gòu)修復(fù)這一過程,在OpenSEES中將刪除單元和建立單元這兩個(gè)功能相結(jié)合,通過單元替換的方式實(shí)現(xiàn)對(duì)震損結(jié)構(gòu)中損傷構(gòu)件的修復(fù)與加固過程,其具體流程如圖3所示。

    精準(zhǔn)的震損結(jié)構(gòu)損傷評(píng)估應(yīng)該對(duì)所有構(gòu)件的損傷程度進(jìn)行綜合評(píng)價(jià),在修復(fù)加固時(shí)也宜從構(gòu)件層面著手,通過對(duì)損傷構(gòu)件進(jìn)行精細(xì)修復(fù)使結(jié)構(gòu)的使用功能得以恢復(fù)。因此,需要根據(jù)構(gòu)件損傷指數(shù)對(duì)構(gòu)件損傷程度進(jìn)行量化,判斷是否需要修復(fù)并選擇合理的修復(fù)加固方式。隨后,通過刪除單元命令將損傷構(gòu)件單元?jiǎng)h除,在原損傷構(gòu)件單元的節(jié)點(diǎn)上新建單元,并根據(jù)修復(fù)時(shí)所采取的材料設(shè)置新建單元的材料屬性。另外,也可通過增加防屈曲支撐等耗能構(gòu)件單元提高修復(fù)結(jié)構(gòu)的抗震能力。

    4 損傷表征指標(biāo)及震損結(jié)構(gòu)快速修復(fù)策略

    對(duì)構(gòu)件損傷進(jìn)行評(píng)估量化的常用物理量指標(biāo)包括最大位移角、殘余位移(角)、端部曲率、剛度變化系數(shù)、頻率變化率等。大量研究者圍繞不同的表征損傷物理量指標(biāo)給出了代表不同損傷程度的規(guī)定限值。如崔濟(jì)東等[19]通過大量試驗(yàn)求得了RC梁、柱及剪力墻等在不同破壞程度時(shí)所對(duì)應(yīng)最大位移角限值,為根據(jù)最大位移角對(duì)構(gòu)件損傷進(jìn)行量化提供理論依據(jù)。

    受荷載作用方向以及配筋等因素的影響,構(gòu)件在兩個(gè)主軸方向上的震損程度并不相同。本文建議按以下公式對(duì)各損傷表征指標(biāo)進(jìn)行綜合處理后再進(jìn)行損傷評(píng)估。

    等效最大位移角θm的計(jì)算公式為:

    裂縫作為構(gòu)件損傷最直接的反映,其數(shù)量和寬度能直接表征損傷程度。為對(duì)震損結(jié)構(gòu)中所有梁柱構(gòu)件損傷程度實(shí)現(xiàn)更全面充分的量化,建議將裂縫與表征損傷的物理量指標(biāo)結(jié)合使用,計(jì)算公式如下:式中 d ir表示所有梁柱構(gòu)件中第i個(gè)構(gòu)件第7-條裂縫的寬度;Si表示根據(jù)表征損傷物理量指標(biāo)對(duì)第i個(gè)構(gòu)件的損傷量化值;Smax表示所有梁柱構(gòu)件根據(jù)表征損傷物理量指標(biāo)進(jìn)行損傷量化中的最大值。

    在基于有限元模型進(jìn)行結(jié)構(gòu)整體損傷評(píng)估時(shí),也可采用剛度變化系數(shù)22和頻率變化率[23]等指標(biāo)來實(shí)現(xiàn),其特點(diǎn)是考慮了完好結(jié)構(gòu)與震損結(jié)構(gòu)的性能差異。剛度變化系數(shù)的計(jì)算公式如下:

    由于結(jié)構(gòu)的頻率與剛度相關(guān),因此頻率變化率是結(jié)構(gòu)損傷的直接表征。何浩祥等[24]提出基于基本頻率變化的時(shí)變損傷指標(biāo)及相關(guān)有限元分析方法,以反映結(jié)構(gòu)在地震作用下任意瞬時(shí)的損傷程度及演變過程,其計(jì)算式如下:

    D(t)=1- [ To/T.(t)]2=1 = [f0(t)/f0]2(23)式中 T0和f0分別表示結(jié)構(gòu)在彈性狀態(tài)下的周期和頻率;Te(t)和fe(t)分別表示結(jié)構(gòu)的塑彈性時(shí)變周期和彈塑性時(shí)變頻率。對(duì)于結(jié)構(gòu)有限元模型,可通過提取結(jié)構(gòu)質(zhì)量矩陣和瞬時(shí)剛度矩陣并進(jìn)行特征值分解來獲得瞬時(shí)頻率;對(duì)于有監(jiān)測(cè)信號(hào)的實(shí)際結(jié)構(gòu),可以通過小波包分解或短時(shí)傅里葉變換等信號(hào)處理技術(shù)來提取結(jié)構(gòu)的瞬時(shí)頻率。

    綜上,在無監(jiān)測(cè)的實(shí)際結(jié)構(gòu)震損評(píng)估中,可以將殘余位移(角)和裂縫特征作為損傷表征指標(biāo),對(duì)有監(jiān)測(cè)的實(shí)際結(jié)構(gòu),可將最大位移角、端部曲率、剛度變化系數(shù)、頻率變化率等作為補(bǔ)充指標(biāo)。在有限元分析中,可以將基于彈塑性耗能差的損傷指數(shù)作為構(gòu)件或結(jié)構(gòu)級(jí)別的損傷評(píng)估指標(biāo),而將頻率變化率或剛度變化系數(shù)作為整體損傷評(píng)估指標(biāo)。

    在地震發(fā)生后,為保障抗震救援的順利開展,亟需在較短時(shí)間內(nèi)對(duì)諸如醫(yī)院等關(guān)鍵震損結(jié)構(gòu)進(jìn)行快速修復(fù)加固以使其恢復(fù)部分或全部使用功能,從而保證結(jié)構(gòu)在遭遇余震時(shí)具有一定的抗震能力。由于機(jī)理有所區(qū)別,采用不同的損傷表征指標(biāo)對(duì)震損結(jié)構(gòu)梁柱構(gòu)件損傷進(jìn)行量化時(shí),得到的構(gòu)件損傷程度排序并不一致。因此,需要比對(duì)分析并篩選出能夠最大限度提高震損結(jié)構(gòu)抗震性能的損傷表征指標(biāo),并以此為依據(jù)指導(dǎo)震損結(jié)構(gòu)的快速修復(fù)加固,進(jìn)而形成優(yōu)化策略。震損結(jié)構(gòu)需要實(shí)現(xiàn)的優(yōu)化目標(biāo)或策略為:在僅修復(fù)加固有限個(gè)損傷梁柱構(gòu)件的前提下能夠最大限度地提升震損結(jié)構(gòu)抗震能力。gzslib202204041633

    針對(duì)震損結(jié)構(gòu)的快速修復(fù)加固,通常存在兩種不同的基本策略:一種是直接對(duì)結(jié)構(gòu)整體中損傷最嚴(yán)重的樓層或構(gòu)件群進(jìn)行修復(fù)以保證其不會(huì)喪失使用功能(即擇層修復(fù));另一種認(rèn)為需要對(duì)每一層中損傷嚴(yán)重構(gòu)件均進(jìn)行修復(fù)(即擇件修復(fù)),以防止剛度突變發(fā)生薄弱層轉(zhuǎn)移的現(xiàn)象。兩種修復(fù)加固策略均具有相應(yīng)的理論基礎(chǔ),在指導(dǎo)震損結(jié)構(gòu)快速修復(fù)時(shí)的效果優(yōu)劣尚需要進(jìn)一步研究和探討。

    5 實(shí)例分析

    為驗(yàn)證在OpenSEES中通過單元替換模擬結(jié)構(gòu)修復(fù)加固過程的可行性和準(zhǔn)確性,選擇一個(gè)二維彈性混凝土門式框架結(jié)構(gòu)進(jìn)行驗(yàn)證。結(jié)構(gòu)模型如圖4所示,梁、柱截面尺寸分別為0.25 m×0.50 m和0.50 m×0.50 m,彈性模量為3.0×104MPa。對(duì)結(jié)構(gòu)進(jìn)行地震時(shí)程分析,隨后對(duì)單元1進(jìn)行單元替換并再次輸入地震波進(jìn)行時(shí)程分析,分別比較在不替換單元、替換單元但不修改彈性模量、替換單元且將彈性模量設(shè)為原來2倍或0.5倍這四種T況。主要分析頂點(diǎn)3處的位移響應(yīng),結(jié)果如圖5所示。

    由節(jié)點(diǎn)3的結(jié)果可以看出:在不替換單元和替換單元但不修改彈性模量兩種T況下得到的節(jié)點(diǎn)3處位移時(shí)程完全一致;在替換單元并將彈性模量分別設(shè)為原來0.5倍和2倍后,位移變化則分別大于、小于不替換單元時(shí)節(jié)點(diǎn)3處的位移變化。由此可見,通過單元替換可以準(zhǔn)確模擬損傷構(gòu)件修復(fù)后的效果,這為分析震損結(jié)構(gòu)修復(fù)加固后的性能提供了技術(shù)支持。

    下面對(duì)比分析按不同的損傷表征指標(biāo)指導(dǎo)修復(fù)有限個(gè)構(gòu)件后對(duì)震損結(jié)構(gòu)性能的提升程度,從而為實(shí)際震損結(jié)構(gòu)快速修復(fù)加固提供理論指導(dǎo)。以7度區(qū)某5層鋼筋混凝土框架結(jié)構(gòu)作為研究對(duì)象。結(jié)構(gòu)的平面布置如圖6所示,結(jié)構(gòu)層高3.3 m,混凝土強(qiáng)度C30,梁、柱截面尺寸分別為:250 mm×500 mm,500 mm×500 mm,梁、柱配筋率分別為:1.96%,1.57%。在OpenSEES中建立該混凝土框架結(jié)構(gòu)模型,其中混凝土和鋼筋的本構(gòu)模型分別選用Con-crete01和Steel01,單元類型選用纖維單元。選用最大峰值加速度(PGA)為0.3g的El Centro波作為地震輸入,使框架結(jié)構(gòu)產(chǎn)生損傷成為震損結(jié)構(gòu)。

    取最大位移角、殘余位移角、曲率和剛度變化系數(shù)作為損傷表征指標(biāo),并根據(jù)式(18)-(23)對(duì)所有構(gòu)件的損傷程度進(jìn)行量化,結(jié)果如圖7所示。可以看出按不同損傷表征指標(biāo)對(duì)構(gòu)件損傷進(jìn)行量化時(shí)得到嚴(yán)重次序并不一致,因而待修復(fù)的構(gòu)件也不相同。

    根據(jù)損傷量化結(jié)果,以單元替換的方法分別對(duì)震損結(jié)構(gòu)中損傷最嚴(yán)重的前16,24和32個(gè)柱構(gòu)件進(jìn)行修復(fù)。隨后輸入PGA為0.3g的El Centro波,根據(jù)基于彈塑性耗能差的損傷評(píng)估方法對(duì)震損結(jié)構(gòu)修復(fù)加固后再遇地震時(shí)的整體損傷程度進(jìn)行量化,以對(duì)比不同損傷表征指標(biāo)指導(dǎo)震損結(jié)構(gòu)修復(fù)的效果優(yōu)劣,得到的結(jié)果如圖8所示。此外,在地震作用下,不同修復(fù)方案下用頻率變化率表征損傷演變曲線如圖9所示。

    由圖8可以看出,以殘余位移角和曲率作為損傷表征指標(biāo)指導(dǎo)震損結(jié)構(gòu)快速修復(fù)后結(jié)構(gòu)整體損傷指數(shù)最小,抗震性能最優(yōu),而其他指標(biāo)對(duì)應(yīng)的結(jié)果并不理想。同時(shí),隨著修復(fù)損傷構(gòu)件數(shù)量的逐漸增多,震損結(jié)構(gòu)在余震下的整體損傷指數(shù)逐漸降低,抗震性能得到逐步提高。由圖9也可以看出按殘余位移角和曲率指導(dǎo)震損結(jié)構(gòu)修復(fù)后,余震下的損傷程度最小。

    將曲率和殘余位移角作為損傷表征指標(biāo),分別按擇層修復(fù)和擇件修復(fù)兩種策略修復(fù)16,24,32個(gè)柱構(gòu)件,并輸入PGA為0.3g的El Centro波,分析不同修復(fù)策略對(duì)其抗震性能的影響。得到震損結(jié)構(gòu)修復(fù)后再歷地震時(shí)的整體損害指數(shù)如圖10所示。

    可以看出盡管根據(jù)兩種策略對(duì)震損結(jié)構(gòu)進(jìn)行修復(fù)都能極大提高結(jié)構(gòu)的抗震性能,但按擇層修復(fù)策略修復(fù)后表現(xiàn)出更好的效果。因此,在對(duì)震損結(jié)構(gòu)進(jìn)行快速修復(fù)時(shí),宜直接對(duì)結(jié)構(gòu)整體中損傷最嚴(yán)重的樓層或構(gòu)件群進(jìn)行修復(fù)使震損結(jié)構(gòu)修復(fù)后具有更好的抗震性能。

    為進(jìn)一步分析按不同損傷表征指標(biāo)指導(dǎo)震損結(jié)構(gòu)修復(fù)后的抗震性能,同時(shí)分析震損結(jié)構(gòu)修復(fù)后遭遇地震時(shí)產(chǎn)生不同程度損傷的概率,以上文五層鋼筋混凝土框架為例,取12組地震波記錄作為地震輸入進(jìn)行結(jié)構(gòu)地震易損性分析。將PGA和基于彈塑性耗能差的損傷指數(shù)分別作為地震動(dòng)強(qiáng)度指數(shù)和結(jié)構(gòu)損傷評(píng)估指數(shù)??紤]修復(fù)24個(gè)柱構(gòu)件的情況,按式(16)分別計(jì)算PGA從0.1g提高至1.0g時(shí)結(jié)構(gòu)的損傷指數(shù),并按式(17)計(jì)算完好結(jié)構(gòu)、震損結(jié)構(gòu)、震損結(jié)構(gòu)按曲率修復(fù)構(gòu)件、按殘余位移角修復(fù)構(gòu)件以及按殘余位移角和曲率加權(quán)統(tǒng)一修復(fù)構(gòu)件后的地震易損性。其中,加權(quán)統(tǒng)一修復(fù)是指對(duì)于遭受不同地震作用的12個(gè)工況,對(duì)相同構(gòu)件的曲率、殘余位移角相加求和,并根據(jù)結(jié)構(gòu)統(tǒng)一修復(fù)最嚴(yán)重的24個(gè)構(gòu)件。其中,按曲率加權(quán)后統(tǒng)一修復(fù)的24個(gè)損傷嚴(yán)重的構(gòu)件與按殘余位移角加權(quán)后修復(fù)的24個(gè)損傷嚴(yán)重的構(gòu)件相同,可將其表示為同一種。分析得到不同T況下結(jié)構(gòu)在不同地震波下?lián)p傷指數(shù)與PGA的關(guān)系以及相應(yīng)的結(jié)構(gòu)地震易損性曲線分別如圖11和圖12所示。

    由圖11可以看出,震損結(jié)構(gòu)在遭遇余震后的損傷程度遠(yuǎn)大于完好結(jié)構(gòu)遭遇地震時(shí)的損傷程度,因此在地震后需要迅速采取合理的措施對(duì)震損結(jié)構(gòu)進(jìn)行修復(fù)或加固,提高結(jié)構(gòu)的抗震能力,降低結(jié)構(gòu)發(fā)生嚴(yán)重破壞甚至倒塌的風(fēng)險(xiǎn)。僅修復(fù)有限個(gè)損傷構(gòu)件的震損結(jié)構(gòu)在遭遇余地震后的損傷指數(shù)雖然得到降低,且不同修復(fù)方案的結(jié)果不同,但還是大于完好結(jié)構(gòu)在主震下的損傷指數(shù)。在三種修復(fù)方案中,根據(jù)殘余位移角對(duì)震損結(jié)構(gòu)進(jìn)行修復(fù)后,結(jié)構(gòu)發(fā)生嚴(yán)重破壞以及倒塌的概率最小,修復(fù)后結(jié)構(gòu)的抗震性能提高最顯著。同時(shí),由于實(shí)際震損結(jié)構(gòu)的殘余位移角恰好可以快速測(cè)量,因而宜將殘余位移角作為首選表征損傷物理量進(jìn)行震損評(píng)估。綜上,在地震災(zāi)害后可根據(jù)實(shí)際震損結(jié)構(gòu)中所有構(gòu)件的殘余位移角以及裂縫尺寸和數(shù)量對(duì)損傷程度進(jìn)行量化,并根據(jù)量化結(jié)果和“擇層修復(fù)”的策略去指導(dǎo)震損結(jié)構(gòu)修復(fù),從而能從最快速度和最大限度兩方面提高結(jié)構(gòu)的抗震能力,以保障人民的生命財(cái)產(chǎn)安全。gzslib202204041634

    6 結(jié) 論

    根據(jù)合理的損傷表征指標(biāo)對(duì)震損結(jié)構(gòu)中構(gòu)件的損傷進(jìn)行量化,并用于指導(dǎo)震損結(jié)構(gòu)的快速修復(fù)加固,能夠有效地提高震損結(jié)構(gòu)抗震性能。本文針對(duì)目前真實(shí)震損結(jié)構(gòu)中構(gòu)件的損傷不易量化且修復(fù)后結(jié)構(gòu)的抗震性能難以評(píng)估的現(xiàn)狀,厘清了位移、曲率、殘余位移以及損傷指數(shù)之間的關(guān)系,提出了在有限元軟件中模擬損傷構(gòu)件修復(fù)加固的方法和流程以及按不同物理量指標(biāo)對(duì)構(gòu)件損傷進(jìn)行量化的方法。同時(shí),在OpenSEES中分析了震損結(jié)構(gòu)按不同的損傷物理量指標(biāo)對(duì)構(gòu)件進(jìn)行快速修復(fù)后的抗震性能。主要結(jié)論如下:

    (1)為克服真實(shí)震損結(jié)構(gòu)可測(cè)信息有限、損傷難以量化的局限,在塑性鉸理論的基礎(chǔ)上,建立梁柱構(gòu)件曲率與位移的關(guān)系,并進(jìn)一步證明梁柱構(gòu)件頂部位移和端部曲率之間具有明確線性關(guān)系、損傷指數(shù)與位移之間存在冪函數(shù)關(guān)系,為根據(jù)殘余位移量化損傷程度提供途徑。

    (2)針對(duì)在有限元軟件中震損結(jié)構(gòu)的修復(fù)加固難以模擬,修復(fù)后結(jié)構(gòu)的抗震性能難以預(yù)估的問題,提出一種通過單元替換模擬震損結(jié)構(gòu)中損傷構(gòu)件進(jìn)行修復(fù)的方法,并經(jīng)過實(shí)例驗(yàn)證可知該方法具有合理性和準(zhǔn)確性。

    (3)建立了根據(jù)不同損傷物理量指標(biāo)對(duì)震損結(jié)構(gòu)中構(gòu)件的損傷進(jìn)行量化的方法,并通過分析震損結(jié)構(gòu)按損傷表征指標(biāo)進(jìn)行快速修復(fù)后再歷地震時(shí)的損傷程度可知,根據(jù)曲率和殘余位移角對(duì)構(gòu)件損傷進(jìn)行量化并用于指導(dǎo)快速修復(fù)時(shí),結(jié)構(gòu)的損傷程度最低,抗震能力提升最明顯。由于在震后檢測(cè)中可快速準(zhǔn)確地獲取殘余位移(角)數(shù)值,因此宜將其作為震損結(jié)構(gòu)最關(guān)鍵的損傷表征指標(biāo)。此外,震損結(jié)構(gòu)的修復(fù)加固宜按照擇層修復(fù)策略進(jìn)行開展。

    (4)由震損結(jié)構(gòu)以及震損修復(fù)結(jié)構(gòu)的地震易損性曲線可知,震損結(jié)構(gòu)再歷地震時(shí)有很大的風(fēng)險(xiǎn)出現(xiàn)嚴(yán)重破壞甚至倒塌的風(fēng)險(xiǎn),而根據(jù)殘余位移角對(duì)震損結(jié)構(gòu)進(jìn)行修復(fù),能夠有效提高結(jié)構(gòu)的抗震能力,降低結(jié)構(gòu)發(fā)生嚴(yán)重破壞和倒塌的風(fēng)險(xiǎn)。這證實(shí)了對(duì)震損結(jié)構(gòu)進(jìn)行修復(fù)加固具有十分重要的意義,為實(shí)際震損結(jié)構(gòu)修復(fù)加固提供有力的理論支撐。

    參考文獻(xiàn):

    [1]李宏男,肖詩云,霍林生.汶川地震震害調(diào)查與啟示[J].建筑結(jié)構(gòu)學(xué)報(bào),2008,29(4):10-19.

    LI Hongnan, XIAO Shiyun, HUO Linsheng. Damageinvestigation and analysis of engineering structures inthe Wenchuan earthquake[J].Journal of Building Struc-tures, 2008, 29(4): 10-19.

    [2] 張超,周云,翁大根,等.震損RC框架結(jié)構(gòu)減震加固的可修性研究[J].土木工程學(xué)報(bào),2014,47(S1):207-213.

    ZHANG Chao, ZHOU Yun, WENG Dagen, et al.Reparability of earthquakedamaged RC frame structurefor seismic retrofit with energy dissipation devices[Jl.Journal of Civil Engineering, 2014, 47( S1): 207-213.

    [3]季靜,喬升訪,韓小雷,等.基于構(gòu)件變形的RC框架結(jié)構(gòu)可修復(fù)能力評(píng)估方法研究[J].建筑結(jié)構(gòu)學(xué)報(bào),2018,39(S1):205211.

    JI Jing,QIAO Shengfang,HAN Xia01ei,et al.RepairIability of RC frame using component performance basedmethod[J].Joumal of Building Structures,2018,39(S1):205 211.

    [4] Park Y,Ang A H.Mechanistic seismic damage modelfor reinforced concrete[J].Joumal of Structural Engi—neenng,1985,111(4):722 739.

    [5]杜修力,歐進(jìn)萍.建筑結(jié)構(gòu)地震破壞評(píng)估模型[J].世界地震工程,1991(3):52 58.

    DU Xiuli,OU Jinping.Earthquake damage evaluanonmodel of building structure[J].World Earthquake Engi—neering,1991(3):52 58.

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

    WANG Dongsheng,F(xiàn)ENG Qimin,WANG Guoxin.A modified Park Ang seismic damage model considerIing low cvcle fatigue life[J].Joumal of Civil Engineer_ing,2004,17(11):41 49.

    [7] 周云,劉季.新型耗能(阻尼)減震器的開發(fā)與研究[J].地震工程與工程振動(dòng),1998,18(1):71-79.

    ZHOU Yun,LIU Ji Development and study of new en—ergy dissipaters(dampers)[J].Earthquake Engineeringand Engineering Vibration,1998,18(1):71-79.

    [8] 中華人民共和國住房和城鄉(xiāng)建設(shè)部.地震災(zāi)后建筑鑒定與加固技術(shù)指南[M],北京:中國建筑工業(yè)出版社,2008.gzslib202204041634

    Ministry of Housing and Urban Rural Development ofthe People s Republic of China.Technical Guide forPost Earthquake SafetV Evaluation and Retrofit of Build—ings[M].Beijing:China Architecture Publishing&Me—dia Co.Ltd..2008.

    [9]歐進(jìn)萍,何政,吳斌,等.鋼筋混凝土結(jié)構(gòu)基于地震損傷性能的設(shè)計(jì)[J].地震工程與工程振動(dòng),1999(1): 21-30.

    OU Jinping,HE Zheng,WU Bin,et al.Seismic dam—age performance based design of reinforced concretestructures[J].Earthquake Engineering and EngineeringVibration,1999(1):21-30.

    [10]張寧.鋼筋混凝土柱變形能力與水平位移計(jì)算研究[D].西安:長安大學(xué),2018.

    ZHANG Ning.Research on deformation ability and hor—izontal displacement calculation of reinforced concretecolumns[D].Xi‘a(chǎn)n:Chan9 an University,2018.

    [11]Priestley M,Park R.Strength and ductility of concretebridge columns under seismic loading[J].ACI Stmctural Joumal,1987,84(1):61 76.

    [12]Kunnath S K,Reinhom A M,Lobo R F.IDARC Version 3.0:A program for the inelastic damage analysis ofreinforced concrete structures[R].Technical ReportNCEER 92 0022,National Center for Earthquake En—gineering ResearCh, State University of New Y0rk,Buffalo,NY.USA.

    [13]何浩祥,陳奎,范少勇.基于彈塑性耗能差率的地震損傷評(píng)估模型及分析方法[J].振動(dòng)工程學(xué)報(bào),2018,31(3):382 390.HE Haoxiang, CHEN Kui, FAN Shaoyong. Seismicdamage model based on differential ratio of elastic plastic dissipated energy and application[Jl. Journal of Vibration Engineering, 2018, 31(3):382-390.

    [14]何浩祥,陳奎,王小兵.梁端填充低屈服點(diǎn)鋼材的梁柱連接減震性能試驗(yàn)與損傷分析[J].建筑結(jié)構(gòu)學(xué)報(bào),2017. 38(5):1-10.

    HE Haoxiang, CHEN Kui. WAN(J Xiaobing. Experiment research and damage analysis on vibration absorption performance of filled steel connection[J].Joumal ofBuilding Structures. 2017, 38(5):1-10.

    [15] KRATZIG W B. MEYER I F. MESKOURIS K.Damage evolution in reinforced concrete members undercyclic loading[Cl. Proceedings of the 5th InternationalConference on Structural Safety and Reliability, SanFrancisco. 1989: 795802.

    [16]何利,葉獻(xiàn)國.Kratzig及ParkAng損傷指數(shù)模型比較研究[J].土木工程學(xué)報(bào),2010,43( 12):1-6.

    HE Li. YE Xianguo. Comparative of Kratzig and ParkAng damage index models[Jl. China Civil EngineeringJournal. 2010, 43( 12): 1-6.

    [17]朱伯龍.結(jié)構(gòu)抗震試驗(yàn)[M].北京:地震出版社,1989.ZHU Bolong. Structural Seismic Test[Ml. Beijing:Earthquake Press, 1989.

    [18]日本道路協(xié)會(huì).道路橘示方耆·同解說:V耐震設(shè)計(jì)編[S].束京:日本道路協(xié)會(huì),2002.

    Association of Highway of Japan. Design specificationsof highway bridges. Part V:seismic design[S].Tokyo:Association of Highway of Japan, 2002.

    [19]崔濟(jì)東.RC梁、柱及剪力墻變形性能指標(biāo)限值研究與試驗(yàn)驗(yàn)證[D].廣州:華南理工大學(xué),2017.

    CUI Jidong.Research and experimental verification ofdefOrmation index limits of RC beams. columns andshear wa11s[D].Guangzhou:South China University ofTechn0109y,2017.gzslib202204041634

    [20]于曉輝.鋼筋混凝土框架結(jié)構(gòu)的概率地震易損性與風(fēng)險(xiǎn)分析[D].哈爾濱:哈爾濱工業(yè)大學(xué),2012.

    YU Xiaohui.Probabilistic seismic fragility and risk anal-ysis of reinforced concrete frame structures[D].Har_bin:Harbin Institute of Techn0109y,2012.

    [21]Mitropoulou C C,Papadrakakis M.Devloping fragilitycurves based on neural network IDA prediction[J].Englneenng Structures,2011,33(12):735 752.

    [22]劉晶波,王文暉,趙冬冬,等.循環(huán)往復(fù)加載的地下結(jié)構(gòu)Pushover分析方法及其在地震損傷分析中的應(yīng)用[J].地震工程學(xué)報(bào),2013,35(1):21 28.

    LIU Jingbo,WANG Wenhui,ZHA()Dongd.ng,eta1.Pushover analysis method of underground structures under reversa1 1oad and its application in seismic damageanalysis[J].Journal of Earthquake Engineering,2013,35(1):21 28.

    [23]Gu X L,Shen Z Y.Damage analysis of reinfOrced concrete structures under earthquake series[C].Proceedings of The 7th International Conference on Computingin Civil and Building Engineering,Seoul,Korea,1997:1019-1024。

    [24]何浩祥,陳奎,閆維明.基于小波包變換和時(shí)變頻率的結(jié)構(gòu)地震損傷評(píng)估[J].振動(dòng)與沖擊,2016,35(7):23 30.

    老司机午夜福利在线观看视频| 成人高潮视频无遮挡免费网站| 亚洲欧美成人精品一区二区| 亚洲成人久久性| 亚洲国产精品sss在线观看| 久久人妻av系列| 一本精品99久久精品77| 欧美日本视频| 成年女人毛片免费观看观看9| 国产综合懂色| 国产精品乱码一区二三区的特点| 国产午夜精品久久久久久一区二区三区 | 久久人妻av系列| 日韩精品中文字幕看吧| 蜜臀久久99精品久久宅男| 国产高清激情床上av| 美女xxoo啪啪120秒动态图| 久久午夜亚洲精品久久| 人妻久久中文字幕网| 日韩高清综合在线| ponron亚洲| 毛片女人毛片| 我要搜黄色片| 国语自产精品视频在线第100页| 欧美极品一区二区三区四区| 精品久久久久久久久av| 两个人的视频大全免费| 成年女人毛片免费观看观看9| 亚洲av成人av| 91av网一区二区| 国产极品精品免费视频能看的| 亚洲一区二区三区色噜噜| 黄色配什么色好看| 亚洲av二区三区四区| 亚洲精品一卡2卡三卡4卡5卡| 全区人妻精品视频| 成人三级黄色视频| 亚洲综合色惰| 色综合站精品国产| 男女之事视频高清在线观看| 国产欧美日韩精品亚洲av| 亚洲av免费高清在线观看| 欧美日韩一区二区视频在线观看视频在线 | av视频在线观看入口| 国产成人a∨麻豆精品| 免费搜索国产男女视频| 最近视频中文字幕2019在线8| 十八禁网站免费在线| 色吧在线观看| 久久精品久久久久久噜噜老黄 | 搡老岳熟女国产| 国产一区二区在线观看日韩| 日韩欧美 国产精品| 成年版毛片免费区| 女同久久另类99精品国产91| 麻豆久久精品国产亚洲av| 日本五十路高清| 亚洲国产精品sss在线观看| 九九久久精品国产亚洲av麻豆| 99久久中文字幕三级久久日本| 欧美日韩在线观看h| 色噜噜av男人的天堂激情| 精品国产三级普通话版| 香蕉av资源在线| 国产男人的电影天堂91| 床上黄色一级片| 一级毛片久久久久久久久女| 97热精品久久久久久| 69av精品久久久久久| av卡一久久| 久久久国产成人免费| 日本撒尿小便嘘嘘汇集6| 亚洲成a人片在线一区二区| 精品午夜福利视频在线观看一区| 欧美另类亚洲清纯唯美| 国产精品伦人一区二区| 国产亚洲欧美98| 桃色一区二区三区在线观看| 一级黄片播放器| 狂野欧美白嫩少妇大欣赏| 亚洲av免费高清在线观看| 美女内射精品一级片tv| 深夜a级毛片| 国产精品一区二区性色av| 看黄色毛片网站| 久久久久久久久中文| 欧美潮喷喷水| 男女那种视频在线观看| 老女人水多毛片| 久久久精品大字幕| 久久久久精品国产欧美久久久| 国产91av在线免费观看| 国产精品久久电影中文字幕| 直男gayav资源| 亚洲av不卡在线观看| 免费看av在线观看网站| 久久精品国产亚洲网站| 午夜精品在线福利| 一个人观看的视频www高清免费观看| 少妇裸体淫交视频免费看高清| 欧美性猛交黑人性爽| 狂野欧美激情性xxxx在线观看| 国产精品日韩av在线免费观看| 日本精品一区二区三区蜜桃| 99热只有精品国产| 91麻豆精品激情在线观看国产| 久久久久久久久大av| 丰满的人妻完整版| 日韩制服骚丝袜av| 久久久久久久久大av| 五月玫瑰六月丁香| 欧美精品国产亚洲| 麻豆国产97在线/欧美| 长腿黑丝高跟| 最近视频中文字幕2019在线8| 久久精品影院6| 国产一区二区激情短视频| 国产精品一区二区三区四区久久| 国产精品人妻久久久影院| 男女边吃奶边做爰视频| 亚洲国产欧洲综合997久久,| 久久精品影院6| 十八禁网站免费在线| 色哟哟·www| 高清日韩中文字幕在线| 色哟哟·www| 亚洲国产欧美人成| 露出奶头的视频| 69av精品久久久久久| 午夜免费激情av| 日韩欧美一区二区三区在线观看| 少妇被粗大猛烈的视频| 久久精品国产亚洲av香蕉五月| 国产蜜桃级精品一区二区三区| 国产av不卡久久| 一级毛片久久久久久久久女| 免费人成视频x8x8入口观看| a级毛片a级免费在线| 床上黄色一级片| 男女之事视频高清在线观看| 亚洲国产精品成人综合色| 精品欧美国产一区二区三| 能在线免费观看的黄片| 国产中年淑女户外野战色| 欧美3d第一页| 欧美国产日韩亚洲一区| 亚洲在线观看片| 免费人成在线观看视频色| 亚洲丝袜综合中文字幕| 国产伦在线观看视频一区| 黄色配什么色好看| 亚洲成人久久性| 国产精品无大码| 99久久中文字幕三级久久日本| 舔av片在线| 一级毛片久久久久久久久女| 一区二区三区免费毛片| 乱系列少妇在线播放| 亚洲第一电影网av| 真实男女啪啪啪动态图| 亚洲精华国产精华液的使用体验 | 亚洲丝袜综合中文字幕| 久久亚洲国产成人精品v| 日韩国内少妇激情av| 九九在线视频观看精品| 国产免费一级a男人的天堂| 91午夜精品亚洲一区二区三区| 精品国产三级普通话版| 久久久精品94久久精品| 免费看美女性在线毛片视频| 岛国在线免费视频观看| 高清毛片免费观看视频网站| 国产亚洲精品久久久com| 亚洲国产精品成人久久小说 | 日日摸夜夜添夜夜爱| 18禁裸乳无遮挡免费网站照片| 成年免费大片在线观看| 国产精品伦人一区二区| 高清毛片免费看| 国产精品野战在线观看| 18+在线观看网站| 国产视频一区二区在线看| 国内少妇人妻偷人精品xxx网站| 日本黄大片高清| 男女啪啪激烈高潮av片| 3wmmmm亚洲av在线观看| 国产精品爽爽va在线观看网站| 乱系列少妇在线播放| 老女人水多毛片| 黄色配什么色好看| 婷婷精品国产亚洲av在线| 99久久中文字幕三级久久日本| 亚洲七黄色美女视频| 久久久久性生活片| 看片在线看免费视频| 人妻夜夜爽99麻豆av| 国产成人91sexporn| 久久天躁狠狠躁夜夜2o2o| 亚洲不卡免费看| 久久鲁丝午夜福利片| www日本黄色视频网| 国产一区二区在线av高清观看| 99久国产av精品国产电影| 亚洲人成网站在线播| 日韩,欧美,国产一区二区三区 | 欧美极品一区二区三区四区| 91狼人影院| 最好的美女福利视频网| 久久人人爽人人爽人人片va| 在线观看午夜福利视频| 精品久久久久久成人av| 国内精品一区二区在线观看| eeuss影院久久| 丰满乱子伦码专区| 午夜精品在线福利| 国产高清三级在线| 小说图片视频综合网站| 五月伊人婷婷丁香| 欧美性猛交黑人性爽| 久久精品夜色国产| 亚洲在线自拍视频| 白带黄色成豆腐渣| 成人三级黄色视频| 最近视频中文字幕2019在线8| 人妻少妇偷人精品九色| 啦啦啦观看免费观看视频高清| 国产精品1区2区在线观看.| 美女免费视频网站| 天美传媒精品一区二区| 亚洲精品456在线播放app| 亚洲,欧美,日韩| 人人妻,人人澡人人爽秒播| 欧美性猛交黑人性爽| 热99re8久久精品国产| 尤物成人国产欧美一区二区三区| 成人高潮视频无遮挡免费网站| 精华霜和精华液先用哪个| 国产精品精品国产色婷婷| 99在线视频只有这里精品首页| 成人美女网站在线观看视频| 色综合色国产| 最近视频中文字幕2019在线8| 国产精品精品国产色婷婷| 午夜激情欧美在线| 国产综合懂色| 中文亚洲av片在线观看爽| 久久婷婷人人爽人人干人人爱| 夜夜爽天天搞| 极品教师在线视频| 欧美一级a爱片免费观看看| 在线观看免费视频日本深夜| 国产午夜精品久久久久久一区二区三区 | 婷婷精品国产亚洲av在线| 欧美性猛交黑人性爽| 欧美xxxx性猛交bbbb| a级一级毛片免费在线观看| 99久久久亚洲精品蜜臀av| 精品午夜福利在线看| 最近手机中文字幕大全| 国产精品久久电影中文字幕| 人妻制服诱惑在线中文字幕| 3wmmmm亚洲av在线观看| 日日摸夜夜添夜夜爱| 别揉我奶头 嗯啊视频| 搞女人的毛片| 一级a爱片免费观看的视频| 91久久精品国产一区二区三区| 亚洲欧美中文字幕日韩二区| 青春草视频在线免费观看| 寂寞人妻少妇视频99o| 天堂√8在线中文| 乱系列少妇在线播放| 午夜日韩欧美国产| 欧美日本亚洲视频在线播放| 搡老妇女老女人老熟妇| 国产精品久久久久久久久免| 精品少妇黑人巨大在线播放 | 久久精品国产鲁丝片午夜精品| 成人亚洲欧美一区二区av| 免费电影在线观看免费观看| 午夜免费激情av| 99精品在免费线老司机午夜| 日本黄大片高清| 国产一区二区激情短视频| 秋霞在线观看毛片| 看黄色毛片网站| 久久国产乱子免费精品| 欧美在线一区亚洲| 亚洲人成网站在线播| 亚洲专区国产一区二区| 日本精品一区二区三区蜜桃| 国产片特级美女逼逼视频| 欧美一级a爱片免费观看看| 国内精品美女久久久久久| 国产精品人妻久久久影院| 国产麻豆成人av免费视频| 亚洲av中文字字幕乱码综合| 日产精品乱码卡一卡2卡三| 狂野欧美激情性xxxx在线观看| 夜夜夜夜夜久久久久| 国产成人一区二区在线| 一级a爱片免费观看的视频| 毛片一级片免费看久久久久| 国产综合懂色| 国产精品综合久久久久久久免费| 久久久精品欧美日韩精品| 老司机午夜福利在线观看视频| 伊人久久精品亚洲午夜| 男女下面进入的视频免费午夜| 一区二区三区免费毛片| 观看免费一级毛片| 日韩人妻高清精品专区| 一个人免费在线观看电影| 欧美丝袜亚洲另类| 久久久久免费精品人妻一区二区| 男人舔女人下体高潮全视频| 六月丁香七月| 精品人妻一区二区三区麻豆 | 亚洲综合色惰| 午夜爱爱视频在线播放| 欧美日韩国产亚洲二区| 色在线成人网| 97热精品久久久久久| 天堂网av新在线| 亚洲精品一区av在线观看| av在线播放精品| 亚洲人成网站高清观看| 午夜福利18| 99riav亚洲国产免费| 人人妻人人看人人澡| 国产高清激情床上av| 熟女人妻精品中文字幕| 内射极品少妇av片p| 免费不卡的大黄色大毛片视频在线观看 | 99九九线精品视频在线观看视频| 老司机福利观看| 女同久久另类99精品国产91| 国产精品日韩av在线免费观看| 欧美成人精品欧美一级黄| 亚洲av二区三区四区| 麻豆国产97在线/欧美| 欧美成人精品欧美一级黄| 精品不卡国产一区二区三区| 天堂√8在线中文| a级毛色黄片| 少妇高潮的动态图| 久久婷婷人人爽人人干人人爱| 禁无遮挡网站| 亚洲人成网站在线观看播放| 亚洲成人精品中文字幕电影| 国产精品久久视频播放| 真实男女啪啪啪动态图| 又粗又爽又猛毛片免费看| 欧美一区二区亚洲| 国产精品久久电影中文字幕| 色哟哟哟哟哟哟| 一级毛片久久久久久久久女| 国产av一区在线观看免费| 成人av一区二区三区在线看| 少妇裸体淫交视频免费看高清| 干丝袜人妻中文字幕| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧美人成| 欧美成人精品欧美一级黄| 亚洲欧美精品自产自拍| 噜噜噜噜噜久久久久久91| 成人二区视频| 欧美日韩一区二区视频在线观看视频在线 | 女人十人毛片免费观看3o分钟| 久久人人爽人人爽人人片va| 国产精品美女特级片免费视频播放器| 色5月婷婷丁香| 九九久久精品国产亚洲av麻豆| 精品欧美国产一区二区三| 大又大粗又爽又黄少妇毛片口| 日本五十路高清| 一区二区三区四区激情视频 | 老司机福利观看| 成人特级黄色片久久久久久久| 国产精品免费一区二区三区在线| 精品人妻熟女av久视频| 亚洲性久久影院| 亚洲精品色激情综合| 欧美丝袜亚洲另类| 国产亚洲精品av在线| 在线观看免费视频日本深夜| 亚洲精华国产精华液的使用体验 | 午夜福利在线在线| av在线老鸭窝| 女人被狂操c到高潮| 99热这里只有是精品50| 大型黄色视频在线免费观看| 免费人成视频x8x8入口观看| 99热精品在线国产| 精品国产三级普通话版| 美女 人体艺术 gogo| 精品福利观看| 长腿黑丝高跟| 欧美成人免费av一区二区三区| 久久久久久九九精品二区国产| 欧美激情在线99| 国产色婷婷99| 亚洲精品日韩av片在线观看| 最近手机中文字幕大全| 日韩国内少妇激情av| 亚洲18禁久久av| 又爽又黄无遮挡网站| 在线天堂最新版资源| 一个人免费在线观看电影| 成人美女网站在线观看视频| 亚洲精品日韩av片在线观看| 99在线视频只有这里精品首页| 99久久九九国产精品国产免费| 国产国拍精品亚洲av在线观看| 91精品国产九色| 日韩av不卡免费在线播放| 久久久久九九精品影院| 在线播放无遮挡| 亚洲av成人av| 日本在线视频免费播放| 国产三级中文精品| 免费观看的影片在线观看| 内射极品少妇av片p| 97在线视频观看| 国产乱人偷精品视频| 午夜影院日韩av| 在现免费观看毛片| 少妇高潮的动态图| 欧美性猛交黑人性爽| 又爽又黄a免费视频| 日本与韩国留学比较| 一本一本综合久久| 色综合色国产| 淫秽高清视频在线观看| 亚洲国产精品国产精品| 色在线成人网| 最近在线观看免费完整版| 国产精品一区二区三区四区免费观看 | 淫秽高清视频在线观看| 国产高清三级在线| 精品久久久久久久久亚洲| 成人性生交大片免费视频hd| 国产精品1区2区在线观看.| 国产精品日韩av在线免费观看| 精品无人区乱码1区二区| 黄色配什么色好看| 国内精品宾馆在线| 一级黄片播放器| 身体一侧抽搐| 免费看美女性在线毛片视频| 亚洲av成人精品一区久久| 91av网一区二区| 国产美女午夜福利| 国产白丝娇喘喷水9色精品| 嫩草影院精品99| 最近最新中文字幕大全电影3| 国产真实伦视频高清在线观看| 99久久精品国产国产毛片| 欧美激情久久久久久爽电影| 精品久久国产蜜桃| 最近最新中文字幕大全电影3| 美女免费视频网站| 在线看三级毛片| 国产精品综合久久久久久久免费| 国产精品久久视频播放| 欧美日韩国产亚洲二区| 麻豆国产97在线/欧美| 久久久a久久爽久久v久久| 最新中文字幕久久久久| 亚洲最大成人手机在线| 一个人看的www免费观看视频| 免费观看人在逋| 午夜免费男女啪啪视频观看 | 又爽又黄a免费视频| 国产精品无大码| 国产美女午夜福利| 国产一级毛片七仙女欲春2| 国产成人a∨麻豆精品| 欧美+日韩+精品| 欧美性猛交黑人性爽| 亚洲在线观看片| 看非洲黑人一级黄片| 久久99热6这里只有精品| 欧美一区二区精品小视频在线| 99热这里只有精品一区| 国产爱豆传媒在线观看| 禁无遮挡网站| 黄色日韩在线| 男女之事视频高清在线观看| 麻豆一二三区av精品| 亚洲一级一片aⅴ在线观看| 51国产日韩欧美| 精品一区二区三区视频在线观看免费| 久久天躁狠狠躁夜夜2o2o| 在线国产一区二区在线| 麻豆国产97在线/欧美| 亚洲久久久久久中文字幕| 日日啪夜夜撸| 亚洲最大成人手机在线| 波野结衣二区三区在线| 激情 狠狠 欧美| 国产在线男女| 亚洲精品粉嫩美女一区| 国产成人a区在线观看| 免费av毛片视频| 欧美高清性xxxxhd video| 精华霜和精华液先用哪个| 国产片特级美女逼逼视频| 国产亚洲精品久久久久久毛片| 色播亚洲综合网| 搞女人的毛片| av在线天堂中文字幕| 精品国内亚洲2022精品成人| 久久国产乱子免费精品| aaaaa片日本免费| 国产三级中文精品| 熟女人妻精品中文字幕| 1000部很黄的大片| 99久国产av精品国产电影| 精品不卡国产一区二区三区| 嫩草影院入口| 免费无遮挡裸体视频| 18禁在线播放成人免费| 午夜视频国产福利| 国产精品日韩av在线免费观看| 联通29元200g的流量卡| 精品午夜福利在线看| av在线老鸭窝| 日韩欧美三级三区| 天堂影院成人在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av五月六月丁香网| 两个人视频免费观看高清| 特级一级黄色大片| 淫秽高清视频在线观看| 看片在线看免费视频| 欧美性猛交黑人性爽| 少妇的逼水好多| 麻豆国产97在线/欧美| 国产精品电影一区二区三区| 十八禁网站免费在线| 99热精品在线国产| 变态另类丝袜制服| 久久久精品大字幕| 国产精品一区二区三区四区久久| 国产高清三级在线| 国产精品,欧美在线| 丝袜喷水一区| 国产精品人妻久久久影院| 国产色婷婷99| 国产黄a三级三级三级人| 日日摸夜夜添夜夜添小说| 久久久久久久午夜电影| 赤兔流量卡办理| 九九热线精品视视频播放| 精品久久久久久久人妻蜜臀av| 色综合站精品国产| 日韩一区二区视频免费看| 国产国拍精品亚洲av在线观看| 久久久久久久亚洲中文字幕| 日本精品一区二区三区蜜桃| 亚洲欧美日韩高清在线视频| 天堂√8在线中文| 国产激情偷乱视频一区二区| 淫秽高清视频在线观看| 91狼人影院| 成人午夜高清在线视频| 亚洲欧美日韩东京热| 色播亚洲综合网| 精品不卡国产一区二区三区| 男插女下体视频免费在线播放| 在线a可以看的网站| 亚洲高清免费不卡视频| 久久精品国产亚洲av涩爱 | 国产伦精品一区二区三区视频9| 美女高潮的动态| 老熟妇仑乱视频hdxx| 久久精品国产亚洲网站| 美女免费视频网站| 欧美高清性xxxxhd video| 亚洲av.av天堂| 国内精品久久久久精免费| 国产一区二区在线av高清观看| 小蜜桃在线观看免费完整版高清| 激情 狠狠 欧美| 国产伦在线观看视频一区| 精品免费久久久久久久清纯| 人妻制服诱惑在线中文字幕| 欧美xxxx黑人xx丫x性爽| 国产亚洲精品久久久久久毛片| 久久久久九九精品影院| 国产免费一级a男人的天堂| 国产精品一区www在线观看| 久久久午夜欧美精品| 最新在线观看一区二区三区| 日韩中字成人| 乱码一卡2卡4卡精品| 哪里可以看免费的av片| 亚洲精华国产精华液的使用体验 | 偷拍熟女少妇极品色| 日韩国内少妇激情av| 亚洲精品影视一区二区三区av| 插逼视频在线观看| 1024手机看黄色片| 亚洲精品日韩在线中文字幕 | 校园春色视频在线观看| 99热网站在线观看| 国产欧美日韩精品一区二区| 中文字幕熟女人妻在线| 成人永久免费在线观看视频| 一夜夜www| h日本视频在线播放| 男人的好看免费观看在线视频| 久久久久久九九精品二区国产| 色视频www国产|