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

    基于日本311大地震詳細(xì)破壞數(shù)據(jù)的結(jié)構(gòu)易損性分析

    2022-07-21 08:24:14仝文博王自法趙登科
    地震工程與工程振動 2022年3期
    關(guān)鍵詞:插值法易損性震動

    仝文博,王自法,趙登科,張 昕

    (1.河南大學(xué)土木建筑學(xué)院,河南開封 475004;2.中國地震局工程力學(xué)研究所,黑龍江哈爾濱 150080;3.中震科建(廣東)防災(zāi)減災(zāi)研究院,廣東韶關(guān) 512026)

    引言

    建筑物地震易損性是指建筑物在不同地震作用下所遭受的破壞程度,從工程的角度,可以把它看成是地震動輸入和結(jié)構(gòu)破壞之間的關(guān)系。結(jié)構(gòu)易損性關(guān)系的研究對震前建筑物的抗震性能和震后建筑損傷的快速準(zhǔn)確評估具有重要意義。影響易損性的因素主要有2 個:地震動和建筑物結(jié)構(gòu)信息。地震動的預(yù)測一般采用衰減模型來預(yù)測[1],由于地震動衰減模型有著相當(dāng)大的不確定性[2],通過衰減模型得到的地震易損性曲線同樣具有很大的不確定性。對于已經(jīng)發(fā)生的地震,建筑物損失的預(yù)測結(jié)果與實際的震害數(shù)據(jù)經(jīng)常有著較大的差別,雖然有時對這些預(yù)測誤差可能有合理的解釋,例如Wang 等[3]對2008年汶川大地震損失的預(yù)估,但是減少預(yù)測的誤差一直是地震損失研究的重點課題之一。為了改善地震動的預(yù)測精度,基于震后觀測到的大量臺站地震動實測數(shù)據(jù)插值得到近似實測地震動,研究基于實測地震動的地震易損性,并通過對比基于衰減模型的地震易損性研究,對有效地提高震后地震損失的預(yù)測精度就顯得尤為重要[4]。

    事實上,基于實際的震害數(shù)據(jù)以統(tǒng)計的方式進行易損性分析的方法差不多是最早提出的易損性分析方法,它是由Whitman 等[5]提出的。該方法的實質(zhì)就是基于震害的群體損失給出地震動(通常是烈度)和損失之間的矩陣關(guān)系。因為大地震不常發(fā)生,并且大地震發(fā)生后也難以收集詳細(xì)的建筑結(jié)構(gòu)損傷資料,所以基于實際詳細(xì)震害的易損性分析一直進展不大,但是,如果能夠得到大量詳細(xì)的建筑結(jié)構(gòu)損傷數(shù)據(jù),基于實際的震害數(shù)據(jù)得到的易損性關(guān)系最為準(zhǔn)確[4,6]。Lin 等[7]利用新西蘭Canterbury 地震中獲取的震害房屋數(shù)據(jù),得到了結(jié)構(gòu)易損性曲線,并將其應(yīng)用于建筑震后易損性評估?;趯嶋H震害的單體地震易損性研究總的來說存在幾個方面的問題:(1)缺少大量翔實的地震震害資料來驗證結(jié)果的正確性。目前國內(nèi)外很少有大量詳細(xì)的建筑物單體破壞的資料用來驗證地震易損性的精度和適用性。(2)衰減模型對于具體地震動預(yù)測的精度有待驗證,特別是需要利用單個地震中觀測到的密集臺網(wǎng)記錄來檢驗地震動衰減模型的預(yù)測結(jié)果。(3)預(yù)測結(jié)果的離散性沒有能夠得到很好地量化。關(guān)于地震損失離散性的分布規(guī)律,文中的相關(guān)作者等利用新西蘭詳細(xì)的地震損失數(shù)據(jù)提出了考慮邊界條件的混合概率分布模型[8]。

    為了解決上述3 個問題,最好的途徑就是通過大量詳細(xì)的實際震害資料和密集臺網(wǎng)實測的地震動進行單體建筑物易損性的研究。2011年3月11日發(fā)生的日本311地震,雖然造成了極大的人員傷亡和經(jīng)濟損失,但是也留下了寶貴的資料,包括詳細(xì)的地震破壞和大量的地震動觀測記錄。文中計劃利用近230 萬棟建筑的資料,其中包括近60 萬棟震后損傷建筑物的詳細(xì)破壞記錄,通過2 種方法獲取建筑場地的地震動信息:(1)采用日本最新的地震動衰減模型獲得建筑物場地的地震動信息;(2)基于臺站實測的地震動觀測數(shù)據(jù),通過空間插值的方法得到建筑場地的地震動信息。利用實際震害得出2 種地震動作用下的易損性關(guān)系,比較兩者結(jié)果的差異,為更加準(zhǔn)確地估算震后建筑物損傷提供更好的方法。

    1 地震數(shù)據(jù)分析

    1.1 地震事件信息

    日本處于亞歐板塊和太平洋板塊交界地帶,是地震多發(fā)的國家。2011年3月11日,在日本東北太平洋沿岸發(fā)生了里氏9 級特大地震,以及隨后一系列的余震并在日本其他地方也發(fā)生了2 次6 級以上的地震,這些地震導(dǎo)致日本東北等許多地區(qū)的房屋遭到嚴(yán)重破壞。Validus Research 收集了震后近60 萬棟建筑的詳細(xì)破壞數(shù)據(jù)(建筑物分布如圖1),其中顏色的深淺代表建筑損傷程度,顏色越深表示建筑的受損程度越嚴(yán)重,顏色越淺表示建筑的受損程度越輕微。對建筑損傷空間分布的分析發(fā)現(xiàn)這些建筑物主要受到了3次地震的影響,分別是1次俯沖帶板緣地震和2次淺地殼地震,3次地震的詳細(xì)信息見表1。

    圖1 震后損傷建筑分布圖Fig.1 Distribution of damaged buildings after the earthquakes

    表1 影響房屋損失的3次地震Table 1 Three earthquakes that caused loss to buildings

    1.2 觀測臺站信息

    地震動觀測臺網(wǎng)(臺站分布如圖2)數(shù)據(jù)來源于日本國家防災(zāi)科學(xué)技術(shù)研究所(NIED)的Kyoshin Network(K-NET)和Kiban Kyoshin Network(KIK-NET)數(shù)據(jù)庫,共987個臺站的數(shù)據(jù)信息。臺站記錄的數(shù)據(jù)信息包含三分向加速度、震源的經(jīng)緯度、震源深度和臺站的經(jīng)緯度等,臺站間隔一般不超過20 km[9-10]。

    圖2日本臺站分布圖Fig.2 Distribution of Japanese stations

    1.3 建筑物數(shù)據(jù)

    建筑年代和建筑結(jié)構(gòu)類別等是結(jié)構(gòu)易損性的重要影響因素[11],本研究根據(jù)Validus Research 所統(tǒng)計的建筑數(shù)據(jù)信息,將全部建筑和損傷建筑分別進行分類。該數(shù)據(jù)庫中建筑的建造年代大致按照明治(1867~1911年)、大正(1912~1925年)、昭和(1926~1988年)和平成(1988+)統(tǒng)計歸類,由于建造年代處于明治和大正2組的建筑數(shù)量極少,文中將明治、大正及昭和年代間的建筑數(shù)據(jù)合并后統(tǒng)稱為“建造年代為1988年之前的建筑”。建筑按建造年代分類時,僅分為1988年之前和1988年之后2組建筑數(shù)據(jù)(見表2和表3)。另一種是將建筑按照結(jié)構(gòu)類別分類:分為木結(jié)構(gòu)、鋼結(jié)構(gòu)和磚結(jié)構(gòu)共3組建筑數(shù)據(jù)。

    表2 全部建筑中不同類別建筑的數(shù)量和占比Table 2 The number and proportion of different types of buildings in all buildings

    表3 損傷建筑中不同類別建筑的數(shù)量和占比Table 3 The number and proportion of different types of buildings in damaged buildings

    2 建筑物場地的地震動

    常用的地震動參數(shù)指標(biāo)有2 種:一種是峰值加速度(PGA);另一種阻尼比為5%的加速度反應(yīng)譜值Sa。PGA 相對比較直觀[12-13],但建筑振動分析需要考慮建筑物自振周期的影響,Sa理論上可以更準(zhǔn)確地描述地震動對結(jié)構(gòu)的作用[14-16],因此,文中選取Sa作為地震動參數(shù)進行結(jié)構(gòu)易損性分析,建筑損傷用建筑結(jié)構(gòu)破壞程度的百分比表示[4,17-18]。

    關(guān)于建筑物場地Sa的估算,文中采用2 種方法。一種是基于趙興權(quán)等[1,19]提出的衰減模型,得到每棟建筑場地的衰減Sa;另一種以趙興權(quán)等提出的衰減模型為空間插值依據(jù),根據(jù)建筑物附近臺站實際觀測的地震動數(shù)據(jù)進行空間插值的方法得到每棟建筑場地的插值Sa。

    2.1 基于衰減模型的地震動計算

    趙興權(quán)等2016年提出的俯沖型地震衰減模型,所采用的數(shù)據(jù)庫中的一部分地震數(shù)據(jù)是自1968年以來一共76 次日本俯沖帶地震事件記錄;而淺地殼地震衰減模型數(shù)據(jù)來源是日本76 次淺地殼地震和46 次上地幔地震的6 482條強震記錄。日本2011年的311地震后收集到的損傷數(shù)據(jù)主要受到俯沖型地震和淺地殼地震2種類型地震影響,并由于主要的3次地震震源深度都小于25 km,因此文中選用趙興權(quán)等提出的2種淺源衰減模型,它們的表達式見式(1)和式(2)。

    由于這次收集到的損傷建筑可能受到來自3 次地震事件的影響,首先確定對建筑物場地地震動造成最大影響的地震事件,然后根據(jù)對建筑場地造成影響最大的地震事件信息(震源經(jīng)緯度、震源深度、地震類型等)和建筑的位置及場地等信息,代入到相應(yīng)的衰減模型得到建筑物場地的衰減Sa。

    淺源俯沖型地震的衰減模型為:

    淺源地殼地震的衰減模型為:

    上述2 種衰減模型中的距離參數(shù)Xi,j為斷層距或震源距[1]。式(1)中,fmintS是震級和震源深度的相關(guān)項;gintloge(ri,j)、gintSLloge(xi,j+ 200)和eintSxi.j為震源距或斷層距的相關(guān)項為地震波傳播路徑過程中,穿越火山帶距離的相關(guān)項;loge(A) 為場地類別的相關(guān)項;式(2)中,fmcr是震級和震源深度的相關(guān)項;gcrloge(ri,j)、gcrLloge(xi,j+ 200)、gN(ci,j)和ecrxi.j為震源距或斷層距的相關(guān)項為地震波傳播路徑過程中,穿越火山帶距離的相關(guān)項;loge(Acr)為場地類別的相關(guān)項;式(1)和式(2)中的距離參數(shù)Xi,j為斷層距或震源距,ξi,j和ηi為誤差的相關(guān)項[1]。對于311 地震,其斷層模型大小選取510 km×210 km[21],淺源俯沖型地震衰減模型(式1)的距離參數(shù)Xi,j選取為斷層距。由于兩次淺地殼地震類型缺少斷層模型,淺源地殼地震衰減模型(式2)的距離參數(shù)Xi,j選取為震源距。其它詳細(xì)參數(shù)見趙興權(quán)等[1,19-20]提出的衰減模型。

    首先利用311地震實測地震動檢驗一下所選取衰減模型的地震動預(yù)測精度。因為俯沖型地震的衰減模型斷層距適用范圍是大于25 km,所以選取離斷層距25 km 以外的臺站計算地震動,并將計算值和實測的值進行比較,從而檢驗衰減模型的精度。由于本次數(shù)據(jù)中的房屋主要是1 層和2 層居多,所以選取自振周期0.1 s時的Sa和自振周期0.2 s時的Sa進行比較研究,對比結(jié)果的散點圖如圖3所示。

    圖3 衰減模型的地震動和臺站實測地震動的對比Fig.3 Comparison of ground motion between the attenuation model and the observed data at the stations

    圖3 中橫軸表示通過衰減模型計算得到的臺站地震動數(shù)據(jù),縱軸表示相同臺站實測的地震動數(shù)據(jù)。兩者比較可見,雖然臺站點的實際地震動在總體上略小于根據(jù)衰減模型計算的地震動,在自振周期為0.1 s 時相差12%,在0.2 s時相差4%,但是離散度比較大,在0.1 s和0.2 s時的方差分別是0.089和0.164,特別是實測地震動在不少臺站點要遠遠大于基于衰減模型的地震動??傮w來看,在斷層距25 km之外時,趙興權(quán)等提出的衰減模型對于預(yù)測地震動有著不錯的精度。

    2.2 基于臺站實測值的空間插值地震動

    根據(jù)建筑物附近臺站的實測地震動插值得到建筑場地的地震動需要大量的地震動觀測記錄,同時也要考慮地震動的空間衰減特征[22]。因為臺站的密度有限,本研究計劃利用衰減關(guān)系模型且考慮空間衰減的影響,通過空間插值方法得到建筑場地實測地震動的近似值[23-24]。

    下面具體說明本研究提出的空間插值方法:

    第1步:建筑周邊臺站的選取。對于一棟建筑,先尋找距離建筑5 km之內(nèi)是否存在臺站,若5 km之內(nèi)存在臺站,則根據(jù)此范圍內(nèi)所有臺站記錄的地震動數(shù)據(jù)進行空間插值。若5 km 之內(nèi)沒有臺站時,則將搜索范圍依次擴展到10,15 km,直到20 km。因為日本觀測臺網(wǎng)比較密集,20 km之內(nèi)至少有一個臺站。

    第2步:選取對建筑影響最大的地震事件。建筑周邊僅有一個臺站時,該臺站記錄的地震事件作為最大地震事件。建筑周邊存在多個的臺站時,所有臺站的Sa記錄中,Sa最大值對應(yīng)的地震事件為最大地震事件。

    第3步:計算每個臺站場地的衰減Sa?;谶x取的最大地震事件和臺站信息,代入衰減模型得到每個臺站場地的衰減Sa。若地震事件為9 級地震,則選取淺源俯沖型地震衰減模型(式(1))。若地震事件為6.4 和6.7級地震,則選取淺源地殼地震衰減模型(式(2))。

    第4步:計算臺站場地與建筑場地類型相同時,臺站近似實測的Sa。僅將每個臺站的場地類型替換為建筑的場地類型,再代入衰減模型,得到和建筑場地類型相同時,每個臺站場地對應(yīng)的衰減Sa。再通過最大地震事件作用時,每個臺站實測Sa和臺站場地衰減Sa比值的修正,得到臺站場地與建筑場地類型相同時,每個臺站近似實測的Sa。當(dāng)建筑周邊僅一個臺站時,僅執(zhí)行第5 步計算;建筑周邊存在多個臺站時,則直接執(zhí)行第6,7步計算。

    第5 步:建筑周邊僅一個臺站時,即計算建筑場地的插值Sa。建筑周邊僅一個臺站時,將建筑信息代入衰減模型,得到建筑場地的衰減Sa,再通過臺站場地與建筑場地的類型相同時,該臺站近似實測的Sa和對應(yīng)的衰減Sa比值的修正。即可得到建筑場地的插值Sa。

    第6 步:建筑周邊存在多個臺站時,計算臺站區(qū)域形心點處近似實測的Sa。存在2 個或2 個以上的臺站時,每個臺站場地的地震動對建筑場地影響的大小不同,按照后述方法引入權(quán)重[25]。僅保留衰減模型中與距離參數(shù)Xi,j的相關(guān)項,即可得到每個臺站對應(yīng)權(quán)重,結(jié)合第4步中每個臺站近似實測的Sa,加權(quán)計算臺站形心點處近似實測的Sa。

    第7 步:建筑周邊存在多個臺站時,計算建筑場地的插值Sa。將形心處場地類型替換為建筑的場地類型,代入衰減模型,得到和建筑場地類型相同時,形心處對應(yīng)的衰減Sa。已知建筑場地處衰減Sa(第5步),再通過形心處場地類型與建筑場地類型相同時,形心處近似實測的Sa值和對應(yīng)衰減Sa比值的修正,即可計算建筑場地處的插值Sa。

    3 2種方法下的結(jié)構(gòu)易損性

    3.1 所有建筑的結(jié)構(gòu)易損性

    文中首先利用日本311地震后所有建筑的資料,分別采用衰減模型和插值法得到每棟建筑場地的Sa,再結(jié)合建筑損傷,得到2種方法下的建筑結(jié)構(gòu)易損性曲線。通過2種方法下的結(jié)構(gòu)易損性對比研究,為更加準(zhǔn)確地估算震后建筑物損傷提供更好的方法。

    文中的建筑破壞程度信息由Validus Research 公司進行統(tǒng)計計算,采用百分比表示,其量化依據(jù)是根據(jù)建筑損失價值與建筑總價值的比值得到。由于建筑損傷在0~5%時,建筑修復(fù)費用為自付額,故建筑損傷在5%以下時,均被記錄為0%。此外,日本311 大地震引發(fā)的海嘯也導(dǎo)致大量建筑破壞,進行建筑結(jié)構(gòu)易損性分析前,已將由海嘯造成破壞的建筑數(shù)據(jù)剔除。由此得到2種地震動下的易損性曲線如圖4所示,其中圖4(a)所示的是基于衰減模型的地震動,而圖4(b)所示的則是利用前述插值方法得到的地震動。圖中的箱型線的上端為損傷建筑數(shù)量超越75%時的分位數(shù),下端為損傷建筑物數(shù)量超越25%時的分位數(shù),箱型中的橫線位置所示的是損傷建筑物中位數(shù)位置,黑色圓圈描述建筑損傷的異常值。大量存在的異常值說明地震損失的離散性極強,即使在相同地震動水平下,結(jié)構(gòu)也會出現(xiàn)不同程度的破壞,且大量存在沒有破壞和全部破壞的建筑物。

    圖4 所有建筑物的結(jié)構(gòu)易損性曲線Fig.4 Structural vulnerability curves for all buildings

    圖4(a)中基于衰減模型的地震動強度最大可以達到1 500 Gal,建筑結(jié)構(gòu)在衰減Sa達到250 Gal 之后逐漸出現(xiàn)損傷,但Sa增加到700 Gal 之后隨著地震動強度的增加,建筑破壞不會出現(xiàn)明顯的增加,這和傳統(tǒng)的認(rèn)識不同,我們稱之為“地震破壞飽和現(xiàn)象”。圖4(b)所示的插值法下的結(jié)構(gòu)易損性曲線,地震動強度最大高達3 000 Gal,由此可見衰減模型對大震地震動的低估現(xiàn)象十分顯著,在震后損傷快速評估的時候有必要利用插值方法對地震動進行修正。和圖4(a)類似,在地震動達到250 Gal 后,建筑同樣開始逐漸出現(xiàn)損傷。當(dāng)?shù)卣饎訌姸瘸^1 000 Gal 后,隨著地震動強度的增加,地震損失不再明顯增加。通過圖4(a),(b)對比可知,相比于插值法下的結(jié)構(gòu)易損性曲線,衰減法下的結(jié)構(gòu)易損性曲線更加陡峭,同一地震動強度下,衰減法下建筑損傷略高于插值法下的建筑損傷。

    3.2 驗證“地震破壞飽和現(xiàn)象”

    2010年9月4日,新西蘭第二大城市的基督城以西40 km 發(fā)生了MW7.1 級地震,在次年的2月份MW6.2級地震伴隨而來,2 次地震均造成了廣泛的破壞,新西蘭國家地震委員會EQC(Earthquake Commission)收集到了150余萬棟房屋在2次地震中的詳細(xì)破壞信息。文中利用新西蘭的震害數(shù)據(jù)對“地震破壞飽和現(xiàn)象”的適用性進行有效地驗證。新西蘭和日本均位于太平洋板塊的邊緣,且都是狹長的島國,因此同樣可以采用趙興權(quán)等提出的衰減模型預(yù)測目標(biāo)場地處的地震動信息。為了更好地指導(dǎo)未來的抗震設(shè)計與工程建設(shè),文中依據(jù)新西蘭于2004年修訂的最新抗震規(guī)范(NZS1170),選取建造年代為2004年之后的建筑進行結(jié)構(gòu)易損性研究,并繪制衰減地震動下的結(jié)構(gòu)易損性曲線(如圖5)。

    圖5 衰減地震動下建造年代為2004年之后的建筑在2次地震影響下的結(jié)構(gòu)易損性曲線Fig.5 Structural vulnerability curves of buildings constructed after 2004 under attenuated ground motions under the influence of twice earthquakes

    圖5(a)中基于衰減模型的地震動強度最大可以達到900 Gal,建筑結(jié)構(gòu)在衰減Sa達到125 Gal 之后逐漸出現(xiàn)損傷,但Sa增加到500 Gal之后隨著地震動強度的增加,建筑破壞不會出現(xiàn)明顯的增加,出現(xiàn)明顯的“地震破壞飽和現(xiàn)象”。圖5(b)和圖5(a)類似,在地震動達到125 Gal后,建筑同樣開始逐漸出現(xiàn)損傷。當(dāng)?shù)卣饎訌姸冗_過700 Gal后,同樣出現(xiàn)明顯的“地震破壞飽和現(xiàn)象”。通過圖5可知“地震飽和現(xiàn)象”有良好的適用性。

    3.3 不同建造年代的結(jié)構(gòu)易損性

    文中將全部建筑按照不同的建造時間,分為1988年之前和1988年之后2 組。由于2 組數(shù)據(jù)中,同一結(jié)構(gòu)類別的建筑占比相差很?。ㄒ姳?),木結(jié)構(gòu)建筑均占80%以上。文中從統(tǒng)計學(xué)角度,選取全部建筑數(shù)據(jù)研究不同建造年代建筑對結(jié)構(gòu)易損性的影響,利用衰減和插值獲取每棟建筑場地的地震動并結(jié)合建筑損傷,繪制衰減和插值2種方法下的不同建造年代的建筑結(jié)構(gòu)易損性曲線(如圖6)。

    表4 不同建造年代建筑中不同結(jié)構(gòu)類別建筑的數(shù)量和占比Table 4 The number and proportion of buildings of different structural types in buildings of different construction years

    圖6 2種地震動對應(yīng)的不同建造年代的結(jié)構(gòu)易損性曲線Fig.6 Building vulnerability curves corresponding to two ground motions of different construction years

    圖6(a),(b)對比2 種方法下,不同建造年代的結(jié)構(gòu)易損性曲線。在衰減和插值法下,不同建造年代的結(jié)構(gòu)易損性均和全部建筑物的易損性關(guān)系大致相同。為研究不同建造年代建筑的易損性相對關(guān)系,文中以建造年代為1988年之前的建筑損傷為基準(zhǔn),計算相同地震動強度下,1988年之后建筑與1988年之前建筑損傷的比值,即“相對損傷比”。衰減方法與插值方法的具體結(jié)果如圖7所示。

    由圖7 可知,無論衰減法還是插值法,相對于1988年之前的建筑,1988年之后建筑的相對損傷比在不同地震動強度下都低于100%,說明在相同方法且同一地震動強度時,1988年之前建造的建筑相比于與1988年之后建造的建筑更容易出現(xiàn)損傷,隨著抗震規(guī)范水平的提高,建筑物的易損性逐步降低。從2 種方法下相對損傷比的均值來看,總體上同一地震動強度下,衰減法時1988年之后的建筑結(jié)構(gòu)損傷約是1988年之前的0.75 倍,而插值法時1988年之后的建筑結(jié)構(gòu)損傷約是1988年之前的0.66倍。

    圖7 2種方法下建造年代為1988年之后的相對損傷比例圖Fig.7 The relative damage ratio of construction with two methods after 1988

    3.4 不同結(jié)構(gòu)類型的結(jié)構(gòu)易損性對比

    為了研究不同結(jié)構(gòu)類型建筑對結(jié)構(gòu)易損性的影響,所有建筑物按照不同的結(jié)構(gòu)類別,分為木結(jié)構(gòu)、鋼結(jié)構(gòu)和磚結(jié)構(gòu)共3組。各種類型建筑物在兩種地震動下對應(yīng)的結(jié)構(gòu)易損性曲線如圖8所示。

    圖8 2種方法下不同結(jié)構(gòu)類型的結(jié)構(gòu)易損性曲線Fig.8 Building vulnerability curves of different structure types with two methods

    通過圖8(a)看出,木結(jié)構(gòu)的結(jié)構(gòu)易損性和包括所有建筑物的易損性關(guān)系比較接近,因為木結(jié)構(gòu)在所有建筑物中占比最大。從圖8(b)可知,衰減法下磚結(jié)構(gòu)的結(jié)構(gòu)易損性也略高于插值得到的結(jié)構(gòu)易損性,且隨著地震動強度的逐漸增加,結(jié)構(gòu)易損性曲線會出現(xiàn)明顯的振蕩現(xiàn)象,主要原因是磚結(jié)構(gòu)的建筑數(shù)量相對較少。此外可以明顯看出,磚結(jié)構(gòu)在地震動強度相對較小時,也易產(chǎn)生損傷,而且隨著地震動強度的增加損傷程度增加較快。表明磚結(jié)構(gòu)的結(jié)構(gòu)易損性較差,抗震性能不好。通過圖8(c)發(fā)現(xiàn),鋼結(jié)構(gòu)的結(jié)構(gòu)易損性和包括所有建筑物的易損性也比較接近。當(dāng)插值法下的地震動強度較大時,建筑損傷會隨著地震動強度的增加產(chǎn)生輕微下降。這可能是由于實際地震害資料中,此時鋼結(jié)構(gòu)的建筑數(shù)量相對較少。

    為研究不同結(jié)構(gòu)類型建筑的易損性相對關(guān)系,文中以木結(jié)構(gòu)建筑的結(jié)構(gòu)損傷為基準(zhǔn),計算相同地震動強度下,鋼結(jié)構(gòu)和磚結(jié)構(gòu)建筑的相對損傷比。衰減方法與插值方法的具體結(jié)果如圖9所示。

    由圖9 可知,無論衰減法還是插值法,隨著地震動強度的增加,鋼結(jié)構(gòu)的相對損傷比全部小于100%,而磚結(jié)構(gòu)的相對損傷比均大于100%,說明在相同方法且同一地震動強度時,鋼結(jié)構(gòu)相比于與木結(jié)構(gòu)的結(jié)構(gòu)易損性較好,抗震性能好,相反磚結(jié)構(gòu)相比于與木結(jié)構(gòu)的結(jié)構(gòu)易損性較差,抗震性能不好。從均值來看,總體上同一地震動強度下,衰減法時鋼結(jié)構(gòu)的建筑損傷約是木結(jié)構(gòu)的0.48 倍,磚結(jié)構(gòu)的建筑損傷約是木結(jié)構(gòu)的2.60倍;插值法時鋼結(jié)構(gòu)的建筑損傷約是木結(jié)構(gòu)的0.43倍,磚結(jié)構(gòu)的建筑損傷約是木結(jié)構(gòu)的1.70倍。

    圖9 2種方法下不同結(jié)構(gòu)類型的相對損傷比例圖Fig.9 The relative damage ratio of different structure types with two methods

    4 結(jié)論與討論

    建筑結(jié)構(gòu)易損性關(guān)系的研究對震前建筑物的抗震性能和震后建筑損傷的快速準(zhǔn)確評估具有重要意義。本研究基于日本近230 萬棟建筑資料和密集臺陣的地震動數(shù)據(jù),提出了基于臺站實測地震動且考慮空間衰減影響的插值方法。并對比衰減模型和插值法2種地震動作用下的結(jié)構(gòu)易損性關(guān)系,得到以下結(jié)論:

    (1)基于衰減模型得到的地震動整體偏小,通過臺站插值得到的地震動更大,衰減模型對地震動存在明顯的低估。同一地震動強度下,衰減法下的建筑損傷比插值法的整體偏高,說明臺站插值方法對更好估算地震損失的重要性。

    (2)無論是衰減法還是插值法,隨著地震動強度達到一定值后,地震損失不再隨著地震動強度的增加明顯增加,其中插值法更為明顯,文中稱這種現(xiàn)象為“震害飽和”現(xiàn)象。本研究組利用新西建造年代為2004年之后的震害建筑數(shù)據(jù)有效驗證了“地震破壞飽和現(xiàn)象”的適用性,說明經(jīng)過抗震設(shè)防的建筑物在遭遇超越設(shè)防水平的地震動時仍然具有良好的抗震能力,這可能是因為2個原因:(1)根據(jù)新西蘭震害資料的詳細(xì)分析,建筑物的損失主要是非結(jié)構(gòu)構(gòu)件的破壞,而結(jié)構(gòu)部分的破壞很有限,因此總損失也有限;(2)極高地震動對建筑的破壞能力并不隨地震動幅值的增加而增加。這個“震害飽和”現(xiàn)象的發(fā)現(xiàn),對于未來巨大地震的震后損傷評估具有重要意義,需要進一步的研究從理論上找出根據(jù)。

    (3)按建造年代分類時,無論是衰減法還是插值法,在相同地震動強度下1988年之后的建筑相比1988年之前的結(jié)構(gòu)損傷明顯降低,說明隨著抗震規(guī)范水平的提高,建筑物的易損性逐步降低。

    (4)按建筑結(jié)構(gòu)分類時,鋼結(jié)構(gòu)的易損性最好,木結(jié)構(gòu)的其次,而磚結(jié)構(gòu)的最差,結(jié)構(gòu)類型對于建筑物易損性有著重要影響。并從側(cè)面驗證了震害資料數(shù)量對于結(jié)構(gòu)易損性準(zhǔn)確性的重要性,只有通過大量詳細(xì)的建筑物單體破壞的資料得到的結(jié)構(gòu)易損性才有良好的精度和適用性。

    (5)得到了不同建造年代和不同結(jié)構(gòu)類別建筑之間的結(jié)構(gòu)易損性相對關(guān)系,能夠為未來不同類別建筑易損性關(guān)系的研究提供一定參考。

    (6)因為衰減模型的不確定性,利用大量真實記錄通過插值法得到的易損性穩(wěn)定性更好,對于特大地震的震后損傷準(zhǔn)確評估有著重要意義。

    猜你喜歡
    插值法易損性震動
    震動減脂儀可以減肥?
    基于IDA的預(yù)應(yīng)力混凝土連續(xù)梁橋易損性分析
    《計算方法》關(guān)于插值法的教學(xué)方法研討
    智富時代(2019年7期)2019-08-16 06:56:54
    水電工程場地地震動確定方法
    振動攪拌 震動創(chuàng)新
    中國公路(2017年18期)2018-01-23 03:00:38
    基于PSDM和IDA法的深水隔震橋梁地震易損性分析比較
    基于二次插值法的布谷鳥搜索算法研究
    Newton插值法在光伏發(fā)電最大功率跟蹤中的應(yīng)用
    基于性能的FRP加固RC框架結(jié)構(gòu)地震易損性分析
    人工合成最不利地震動
    肉色欧美久久久久久久蜜桃| 最黄视频免费看| 大码成人一级视频| 日韩免费高清中文字幕av| 妹子高潮喷水视频| 午夜免费男女啪啪视频观看| 91久久精品电影网| av又黄又爽大尺度在线免费看| 91精品伊人久久大香线蕉| 亚洲精品第二区| 青春草视频在线免费观看| 国产精品一区二区在线观看99| 国产精品成人在线| 久久久久久久精品精品| 国产视频内射| 欧美人与善性xxx| av在线观看视频网站免费| 午夜久久久在线观看| 大香蕉97超碰在线| 纯流量卡能插随身wifi吗| 尾随美女入室| freevideosex欧美| 国产日韩一区二区三区精品不卡 | 午夜影院在线不卡| 你懂的网址亚洲精品在线观看| 亚洲av成人精品一二三区| xxxhd国产人妻xxx| 免费看av在线观看网站| 久久鲁丝午夜福利片| 亚洲激情五月婷婷啪啪| 国产高清国产精品国产三级| 天堂俺去俺来也www色官网| 亚洲精品久久成人aⅴ小说 | 国产精品久久久久久av不卡| 十八禁高潮呻吟视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91成人精品电影| 成人亚洲欧美一区二区av| 日韩伦理黄色片| 精品少妇久久久久久888优播| 三级国产精品片| 国产精品国产av在线观看| 成人毛片a级毛片在线播放| 久久久久视频综合| 亚洲av免费高清在线观看| 极品少妇高潮喷水抽搐| 天堂中文最新版在线下载| 日韩三级伦理在线观看| 国产亚洲精品久久久com| 精品人妻一区二区三区麻豆| 亚洲成人av在线免费| 最新中文字幕久久久久| 91在线精品国自产拍蜜月| 亚洲成人一二三区av| 尾随美女入室| 亚洲怡红院男人天堂| 久久国产亚洲av麻豆专区| 如何舔出高潮| 老司机影院毛片| 草草在线视频免费看| 91久久精品国产一区二区成人| 国产精品久久久久久久久免| 亚洲国产日韩一区二区| 最黄视频免费看| 少妇猛男粗大的猛烈进出视频| 蜜桃久久精品国产亚洲av| 大香蕉久久网| 成人国产麻豆网| 一二三四中文在线观看免费高清| 中文字幕人妻熟人妻熟丝袜美| 老司机影院毛片| 我的女老师完整版在线观看| 国产精品久久久久成人av| 一个人免费看片子| 国产永久视频网站| 精品久久国产蜜桃| 久久久精品区二区三区| 女性生殖器流出的白浆| 各种免费的搞黄视频| 最近手机中文字幕大全| 日韩欧美一区视频在线观看| 欧美日韩成人在线一区二区| 久久 成人 亚洲| 国产欧美日韩一区二区三区在线 | 国产精品久久久久成人av| 夜夜爽夜夜爽视频| 国产精品国产三级专区第一集| 成人影院久久| 国产视频内射| 色94色欧美一区二区| 亚洲欧美成人综合另类久久久| 国产免费视频播放在线视频| videos熟女内射| 久久热精品热| xxxhd国产人妻xxx| 亚洲国产毛片av蜜桃av| 日韩人妻高清精品专区| 丰满迷人的少妇在线观看| 国产69精品久久久久777片| 人人妻人人爽人人添夜夜欢视频| 日本午夜av视频| 中文字幕av电影在线播放| 久久久久久人妻| 国产精品一区二区在线观看99| 中国三级夫妇交换| 亚洲国产精品成人久久小说| 久久精品熟女亚洲av麻豆精品| 国产免费一级a男人的天堂| 国产精品免费大片| 美女视频免费永久观看网站| 国产av精品麻豆| 成人国语在线视频| 简卡轻食公司| 成年人免费黄色播放视频| 日韩精品有码人妻一区| 在线 av 中文字幕| 日日摸夜夜添夜夜添av毛片| 亚洲美女视频黄频| 午夜激情av网站| 午夜福利,免费看| 久久久久精品久久久久真实原创| 一本一本综合久久| 91精品一卡2卡3卡4卡| 国产午夜精品久久久久久一区二区三区| 色视频在线一区二区三区| 中文字幕制服av| 亚洲人成77777在线视频| 成人国产av品久久久| 夜夜看夜夜爽夜夜摸| 在线观看三级黄色| 建设人人有责人人尽责人人享有的| 亚洲国产日韩一区二区| 99久久中文字幕三级久久日本| 在现免费观看毛片| 午夜福利,免费看| 亚洲av成人精品一二三区| 99久久精品国产国产毛片| 成人午夜精彩视频在线观看| 最黄视频免费看| 国产不卡av网站在线观看| 久久99精品国语久久久| 亚洲综合精品二区| 日韩一区二区视频免费看| 免费高清在线观看视频在线观看| av卡一久久| 亚洲国产毛片av蜜桃av| 免费观看a级毛片全部| 国产精品久久久久久av不卡| 秋霞伦理黄片| 国产极品粉嫩免费观看在线 | 国产欧美亚洲国产| 99久久中文字幕三级久久日本| 九色亚洲精品在线播放| av有码第一页| 高清av免费在线| 99热全是精品| 成人亚洲欧美一区二区av| 九色成人免费人妻av| 国产视频内射| 国产精品久久久久久久电影| 精品一区在线观看国产| 久久精品国产a三级三级三级| 精品熟女少妇av免费看| 日本av手机在线免费观看| 亚洲不卡免费看| 亚洲欧美一区二区三区黑人 | 国产亚洲欧美精品永久| 国产无遮挡羞羞视频在线观看| 久久久久久久久久久免费av| 国国产精品蜜臀av免费| 美女cb高潮喷水在线观看| 国产日韩欧美亚洲二区| 午夜免费鲁丝| av卡一久久| 午夜老司机福利剧场| 国产在线一区二区三区精| 18+在线观看网站| 亚洲欧美一区二区三区黑人 | 国国产精品蜜臀av免费| 久久97久久精品| 热99国产精品久久久久久7| 欧美3d第一页| 考比视频在线观看| 大香蕉久久网| 日本与韩国留学比较| 夫妻午夜视频| 黄色视频在线播放观看不卡| 狠狠精品人妻久久久久久综合| 精品国产露脸久久av麻豆| 国产熟女午夜一区二区三区 | 久久精品国产亚洲网站| 亚洲国产精品国产精品| 国产毛片在线视频| 人成视频在线观看免费观看| 亚洲国产精品专区欧美| 一本一本综合久久| kizo精华| 欧美日韩视频精品一区| 日韩三级伦理在线观看| 亚洲五月色婷婷综合| 亚洲三级黄色毛片| 日本免费在线观看一区| 精品久久蜜臀av无| 下体分泌物呈黄色| 老熟女久久久| 国产成人精品福利久久| 日韩中字成人| 久久国内精品自在自线图片| 亚洲av成人精品一区久久| 成年美女黄网站色视频大全免费 | 亚洲丝袜综合中文字幕| 日韩av在线免费看完整版不卡| 国产综合精华液| 久久精品国产a三级三级三级| 人妻夜夜爽99麻豆av| 777米奇影视久久| 秋霞在线观看毛片| 久热久热在线精品观看| 国产不卡av网站在线观看| videosex国产| 久久久久久伊人网av| 国产精品一国产av| 日本vs欧美在线观看视频| 久久久久网色| 亚洲精品一二三| 久久热精品热| 亚洲欧洲日产国产| 午夜激情福利司机影院| 国产成人freesex在线| 欧美人与性动交α欧美精品济南到 | 精品午夜福利在线看| av卡一久久| 内地一区二区视频在线| 成人免费观看视频高清| 极品人妻少妇av视频| 亚洲国产最新在线播放| 欧美 日韩 精品 国产| 一个人看视频在线观看www免费| 国产精品秋霞免费鲁丝片| 精品国产国语对白av| 久久婷婷青草| av福利片在线| 日日摸夜夜添夜夜爱| 免费黄网站久久成人精品| 人成视频在线观看免费观看| 国产国拍精品亚洲av在线观看| 99热全是精品| 日韩,欧美,国产一区二区三区| 亚洲美女黄色视频免费看| 亚洲国产精品一区三区| 亚洲少妇的诱惑av| 欧美xxⅹ黑人| 精品午夜福利在线看| 少妇人妻精品综合一区二区| 九草在线视频观看| 在线观看三级黄色| 能在线免费看毛片的网站| 丝袜在线中文字幕| 午夜福利在线观看免费完整高清在| 看非洲黑人一级黄片| 99久国产av精品国产电影| 一区二区日韩欧美中文字幕 | 亚洲中文av在线| 亚洲精品乱码久久久久久按摩| 国产亚洲午夜精品一区二区久久| 热99国产精品久久久久久7| 欧美日韩亚洲高清精品| 久久av网站| 国产午夜精品久久久久久一区二区三区| 亚洲国产精品一区三区| 欧美3d第一页| 久久精品国产亚洲av涩爱| 欧美人与性动交α欧美精品济南到 | 一级毛片我不卡| videos熟女内射| 久久精品人人爽人人爽视色| 男人操女人黄网站| 日韩中文字幕视频在线看片| 欧美成人精品欧美一级黄| 交换朋友夫妻互换小说| 毛片一级片免费看久久久久| 91精品一卡2卡3卡4卡| 久久人人爽人人爽人人片va| 精品久久久噜噜| 日韩精品有码人妻一区| 在线亚洲精品国产二区图片欧美 | a级毛色黄片| 日本黄大片高清| 男男h啪啪无遮挡| 久久毛片免费看一区二区三区| 国产一区亚洲一区在线观看| 乱码一卡2卡4卡精品| 亚洲天堂av无毛| 51国产日韩欧美| 欧美日韩一区二区视频在线观看视频在线| 91精品国产国语对白视频| 国产精品国产av在线观看| 在线观看三级黄色| 伊人久久国产一区二区| 亚洲伊人久久精品综合| 亚洲不卡免费看| 日本色播在线视频| 妹子高潮喷水视频| 母亲3免费完整高清在线观看 | 日本午夜av视频| 国产乱来视频区| 两个人的视频大全免费| 美女国产高潮福利片在线看| 天堂中文最新版在线下载| 99久久综合免费| 日日摸夜夜添夜夜添av毛片| 国产日韩一区二区三区精品不卡 | 免费不卡的大黄色大毛片视频在线观看| 如何舔出高潮| 日韩熟女老妇一区二区性免费视频| 观看美女的网站| 91精品一卡2卡3卡4卡| 男女边吃奶边做爰视频| 久久久精品免费免费高清| 亚洲精品乱久久久久久| 久久久久久久亚洲中文字幕| 色婷婷av一区二区三区视频| 精品亚洲成a人片在线观看| 久久精品国产自在天天线| 国产亚洲一区二区精品| 国产精品无大码| 最近的中文字幕免费完整| 交换朋友夫妻互换小说| 水蜜桃什么品种好| 亚洲精品乱码久久久久久按摩| 丝袜脚勾引网站| 精品熟女少妇av免费看| 在现免费观看毛片| 国产女主播在线喷水免费视频网站| 我要看黄色一级片免费的| 精品国产露脸久久av麻豆| 18禁在线无遮挡免费观看视频| 亚洲中文av在线| 亚洲内射少妇av| 天天影视国产精品| 十分钟在线观看高清视频www| 国产在线一区二区三区精| 18禁黄网站禁片午夜丰满| 女人被躁到高潮嗷嗷叫费观| 国产极品粉嫩免费观看在线| 老熟女久久久| 最新在线观看一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲免费av在线视频| 老熟妇仑乱视频hdxx| 嫁个100分男人电影在线观看| 侵犯人妻中文字幕一二三四区| 久久免费观看电影| 757午夜福利合集在线观看| 亚洲九九香蕉| 免费日韩欧美在线观看| 国产在线视频一区二区| www.精华液| 97在线人人人人妻| 久久久国产欧美日韩av| 日韩免费av在线播放| 女人久久www免费人成看片| av一本久久久久| 美女国产高潮福利片在线看| 黄色视频不卡| 2018国产大陆天天弄谢| 美女视频免费永久观看网站| 国产亚洲欧美在线一区二区| 最近最新中文字幕大全电影3 | 国产一区二区三区视频了| 亚洲视频免费观看视频| 亚洲欧洲日产国产| 精品一区二区三卡| 久久久久视频综合| 黑人巨大精品欧美一区二区蜜桃| 日本av免费视频播放| 美女福利国产在线| 日韩精品免费视频一区二区三区| 日本vs欧美在线观看视频| 亚洲午夜精品一区,二区,三区| 成人免费观看视频高清| 国产在线精品亚洲第一网站| 涩涩av久久男人的天堂| 久久香蕉激情| 久久久久网色| 久久久国产精品麻豆| 久久人妻av系列| 午夜精品国产一区二区电影| 黄片大片在线免费观看| 99国产精品一区二区三区| 女人久久www免费人成看片| 久久人人爽av亚洲精品天堂| 国产真人三级小视频在线观看| 精品久久久久久电影网| 97人妻天天添夜夜摸| 色94色欧美一区二区| 人人妻人人添人人爽欧美一区卜| 丁香六月天网| 国产熟女午夜一区二区三区| 女同久久另类99精品国产91| 亚洲国产精品一区二区三区在线| av国产精品久久久久影院| 黄片小视频在线播放| 老司机福利观看| 久久精品国产a三级三级三级| 亚洲色图 男人天堂 中文字幕| 国产精品熟女久久久久浪| 国产亚洲欧美在线一区二区| www.熟女人妻精品国产| 国产精品香港三级国产av潘金莲| bbb黄色大片| 飞空精品影院首页| 少妇裸体淫交视频免费看高清 | 精品亚洲乱码少妇综合久久| 久久精品亚洲av国产电影网| 侵犯人妻中文字幕一二三四区| 香蕉丝袜av| 久久99热这里只频精品6学生| 国产精品麻豆人妻色哟哟久久| 69精品国产乱码久久久| 9热在线视频观看99| 久久中文字幕人妻熟女| 国产av国产精品国产| 亚洲美女黄片视频| 香蕉国产在线看| 少妇 在线观看| 中文字幕另类日韩欧美亚洲嫩草| xxxhd国产人妻xxx| 国产精品香港三级国产av潘金莲| 国产97色在线日韩免费| 丝瓜视频免费看黄片| 国产黄频视频在线观看| 午夜福利影视在线免费观看| 19禁男女啪啪无遮挡网站| 久久精品国产99精品国产亚洲性色 | 另类亚洲欧美激情| 亚洲国产av新网站| 一级黄色大片毛片| 黄频高清免费视频| 丰满少妇做爰视频| 久久精品熟女亚洲av麻豆精品| 欧美精品人与动牲交sv欧美| 一本久久精品| 肉色欧美久久久久久久蜜桃| 日日夜夜操网爽| 日韩精品免费视频一区二区三区| 欧美精品人与动牲交sv欧美| 男女高潮啪啪啪动态图| 极品少妇高潮喷水抽搐| 在线观看免费午夜福利视频| 国产精品亚洲一级av第二区| 久久免费观看电影| 美女国产高潮福利片在线看| 99国产综合亚洲精品| 青青草视频在线视频观看| 搡老岳熟女国产| 日韩大片免费观看网站| 久久人人97超碰香蕉20202| a在线观看视频网站| 欧美 亚洲 国产 日韩一| 19禁男女啪啪无遮挡网站| 最近最新中文字幕大全电影3 | 中文字幕av电影在线播放| 色在线成人网| 极品人妻少妇av视频| 超碰成人久久| 黄色视频,在线免费观看| 丰满迷人的少妇在线观看| 精品一区二区三卡| 久久久久视频综合| 亚洲少妇的诱惑av| 国产熟女午夜一区二区三区| 纯流量卡能插随身wifi吗| 久久久精品国产亚洲av高清涩受| 嫩草影视91久久| 极品人妻少妇av视频| 97在线人人人人妻| 十分钟在线观看高清视频www| 女人被躁到高潮嗷嗷叫费观| 69精品国产乱码久久久| 老熟妇仑乱视频hdxx| 一级,二级,三级黄色视频| 五月天丁香电影| 90打野战视频偷拍视频| 久久99热这里只频精品6学生| 捣出白浆h1v1| 美女扒开内裤让男人捅视频| 国产淫语在线视频| 午夜成年电影在线免费观看| 两个人免费观看高清视频| 午夜成年电影在线免费观看| h视频一区二区三区| 午夜激情久久久久久久| 无遮挡黄片免费观看| 夜夜夜夜夜久久久久| 欧美+亚洲+日韩+国产| 黄色片一级片一级黄色片| 国产精品久久久久久精品古装| 后天国语完整版免费观看| 亚洲精品成人av观看孕妇| 欧美精品一区二区大全| 超碰成人久久| 国产精品一区二区在线观看99| 99久久99久久久精品蜜桃| 丝袜美腿诱惑在线| 免费日韩欧美在线观看| 老司机影院毛片| 午夜福利视频精品| 老汉色av国产亚洲站长工具| av一本久久久久| 久久精品aⅴ一区二区三区四区| 法律面前人人平等表现在哪些方面| 国产精品秋霞免费鲁丝片| 精品福利观看| 亚洲男人天堂网一区| av一本久久久久| 欧美精品高潮呻吟av久久| 日韩中文字幕视频在线看片| tube8黄色片| 热99re8久久精品国产| 亚洲成人免费av在线播放| 精品熟女少妇八av免费久了| 亚洲国产欧美日韩在线播放| 国产欧美日韩综合在线一区二区| 在线观看66精品国产| 91国产中文字幕| 亚洲精品国产区一区二| 欧美日韩福利视频一区二区| www.自偷自拍.com| 无人区码免费观看不卡 | 一边摸一边抽搐一进一小说 | 久久久久久久大尺度免费视频| 最黄视频免费看| 亚洲 国产 在线| 操出白浆在线播放| 日韩免费av在线播放| bbb黄色大片| 国产男靠女视频免费网站| 日韩一卡2卡3卡4卡2021年| 日本五十路高清| 免费不卡黄色视频| 99久久人妻综合| 老司机午夜十八禁免费视频| 日韩免费av在线播放| 女警被强在线播放| 亚洲国产成人一精品久久久| 日本av免费视频播放| 少妇粗大呻吟视频| 成人手机av| 午夜精品国产一区二区电影| 国产麻豆69| videos熟女内射| 大码成人一级视频| 精品一品国产午夜福利视频| 2018国产大陆天天弄谢| 久久精品aⅴ一区二区三区四区| 亚洲avbb在线观看| 欧美黑人精品巨大| svipshipincom国产片| 新久久久久国产一级毛片| 国产欧美日韩综合在线一区二区| 精品亚洲成a人片在线观看| 国产精品秋霞免费鲁丝片| 久久婷婷成人综合色麻豆| 一进一出抽搐动态| 亚洲专区中文字幕在线| 国产区一区二久久| 美女午夜性视频免费| 久久人妻av系列| 亚洲,欧美精品.| 伦理电影免费视频| 日韩三级视频一区二区三区| 一区二区三区乱码不卡18| 欧美精品人与动牲交sv欧美| 国产日韩欧美视频二区| 国产欧美日韩综合在线一区二区| 十八禁高潮呻吟视频| 成年女人毛片免费观看观看9 | 老熟妇仑乱视频hdxx| 亚洲色图 男人天堂 中文字幕| 亚洲精品在线观看二区| 午夜激情久久久久久久| 久久天堂一区二区三区四区| 亚洲精品久久午夜乱码| 国产欧美日韩精品亚洲av| 丁香六月欧美| 亚洲欧美一区二区三区黑人| 中文亚洲av片在线观看爽 | 国产精品秋霞免费鲁丝片| 在线观看一区二区三区激情| 亚洲avbb在线观看| 国产欧美日韩一区二区三| 精品第一国产精品| 大香蕉久久网| 高清欧美精品videossex| 最近最新免费中文字幕在线| 999精品在线视频| 精品国产一区二区久久| 亚洲欧美色中文字幕在线| 高清av免费在线| 99re6热这里在线精品视频| netflix在线观看网站| 亚洲av第一区精品v没综合| 午夜两性在线视频| 亚洲国产av新网站| 国产在线一区二区三区精| 80岁老熟妇乱子伦牲交| av天堂久久9| 一二三四社区在线视频社区8| 看免费av毛片| 亚洲性夜色夜夜综合| 超碰成人久久| 国产av又大| 天天影视国产精品| 国产精品一区二区在线不卡|