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

    基于振動(dòng)臺(tái)實(shí)驗(yàn)的結(jié)構(gòu)損傷識(shí)別研究①

    2016-04-07 07:44:11RezaShahneam
    地震工程學(xué)報(bào) 2016年1期

    王 飛, Reza Shahneam

    (1.中國地震局地球物理研究所,北京 100081; 2.北京市地震局,北京100080; 3.美國地質(zhì)調(diào)查局,美國 加州 94025)

    ?

    基于振動(dòng)臺(tái)實(shí)驗(yàn)的結(jié)構(gòu)損傷識(shí)別研究①

    王飛1,2, Reza Shahneam3

    (1.中國地震局地球物理研究所,北京 100081; 2.北京市地震局,北京100080; 3.美國地質(zhì)調(diào)查局,美國 加州 94025)

    摘要:選擇美國加州大學(xué)圣地亞哥分校7層鋼筋混凝土剪力墻足尺結(jié)構(gòu)振動(dòng)臺(tái)實(shí)驗(yàn),開展結(jié)構(gòu)損傷識(shí)別研究,實(shí)驗(yàn)采用白噪聲、環(huán)境振動(dòng)和不同強(qiáng)度的地震動(dòng)交替激發(fā),記錄地震動(dòng)激發(fā)實(shí)驗(yàn)前后的結(jié)構(gòu)反應(yīng)?;谠撚涗浻?jì)算和對(duì)比自振頻率和振型曲率的變化、剪切波走時(shí)及其變化和結(jié)構(gòu)層間位移角,分析發(fā)現(xiàn)一層和二層振型曲率較大,走時(shí)較長,走時(shí)變化也較大,現(xiàn)場(chǎng)檢查發(fā)現(xiàn)一層和二層的破壞也較為嚴(yán)重,這些參數(shù)可用于識(shí)別結(jié)構(gòu)損傷程度和定位損傷位置,而自振頻率和層間位移角變化僅可反映出結(jié)構(gòu)損傷程度,難以揭示結(jié)構(gòu)損傷位置。

    關(guān)鍵詞:振動(dòng)臺(tái)實(shí)驗(yàn); 損傷識(shí)別; 模態(tài)參數(shù); 等效剪切走時(shí); 層間位移角

    0引言

    近年來地震摧毀了許多不同類型的建筑結(jié)構(gòu),因此針對(duì)建筑結(jié)構(gòu)的地震反應(yīng)觀測(cè)研究逐漸引起地震工程學(xué)家的重視,建筑抗震設(shè)計(jì)規(guī)范[1]和許多地方法規(guī)[2]要求重大建筑結(jié)構(gòu)安裝結(jié)構(gòu)地震反應(yīng)觀測(cè)臺(tái)陣。結(jié)構(gòu)地震反應(yīng)觀測(cè)旨在了解和掌握結(jié)構(gòu)體系在強(qiáng)地震作用下的反應(yīng)性狀[3]。一旦獲得結(jié)構(gòu)強(qiáng)震反應(yīng)記錄,就相當(dāng)于對(duì)該結(jié)構(gòu)進(jìn)行了一次原型實(shí)驗(yàn),大量信息可用來揭示結(jié)構(gòu)的抗震性能,這有助于提高結(jié)構(gòu)抗震設(shè)計(jì)水平,改進(jìn)抗震設(shè)計(jì)標(biāo)準(zhǔn),幫助減輕未來的地震災(zāi)害[4]。結(jié)構(gòu)地震反應(yīng)記錄還能通過重建結(jié)構(gòu)地震反應(yīng)的方式來檢驗(yàn)結(jié)構(gòu)地震反應(yīng)分析中數(shù)值模型或?qū)嶒?yàn)?zāi)P偷目煽啃訹5]。健康建筑結(jié)構(gòu)能有效保障人們生命財(cái)產(chǎn)安全,這就需要對(duì)結(jié)構(gòu)開展定期安全性能評(píng)估,尤其是在遭受破壞性地震后,評(píng)估結(jié)果可作為結(jié)構(gòu)使用決策的重要依據(jù)[6]。然而當(dāng)前很多工作都集中在基于設(shè)計(jì)方法的結(jié)構(gòu)抗震能力的理論分析或數(shù)值模擬,很少開展基于振動(dòng)臺(tái)實(shí)驗(yàn)的整體結(jié)構(gòu)抗震性能評(píng)價(jià)。

    基于振動(dòng)的結(jié)構(gòu)健康監(jiān)測(cè)和損傷識(shí)別是以定量方式評(píng)估結(jié)構(gòu)安全性能的有效工具。結(jié)構(gòu)健康監(jiān)測(cè)通常采用現(xiàn)場(chǎng)無損探測(cè)和系統(tǒng)特性分析方法,通過分析時(shí)間域或頻率域中系統(tǒng)特性的變化來表明體系存在損傷或者抗震性能退化[7]。本研究調(diào)研了過去幾十年工程專家們發(fā)展的基于振動(dòng)的結(jié)構(gòu)健康監(jiān)測(cè)和損傷識(shí)別方法[8-10]。這些方法中多數(shù)僅通過某個(gè)結(jié)構(gòu)構(gòu)件或簡(jiǎn)單結(jié)構(gòu)體系的數(shù)值模擬構(gòu)建多個(gè)損傷情景進(jìn)行分析驗(yàn)證。由于結(jié)構(gòu)建模中的不確定性和測(cè)試結(jié)果的簡(jiǎn)化,這種單純基于數(shù)值模型的分析驗(yàn)證并不充分,需通過實(shí)驗(yàn)或?qū)嶋H記錄分析驗(yàn)證,如果可能,開展基于破壞性地震中的現(xiàn)場(chǎng)原型實(shí)驗(yàn)數(shù)據(jù)的分析驗(yàn)證將成為檢驗(yàn)損傷識(shí)別方法的有力工具。盡管當(dāng)前國內(nèi)外建成了一些結(jié)構(gòu)地震反應(yīng)觀測(cè)臺(tái)陣,但目前很難獲得破壞性地震作用下的數(shù)據(jù)記錄,尚不能用于開展損傷識(shí)別的系統(tǒng)研究。因此本工作選擇原型破壞性振動(dòng)臺(tái)實(shí)驗(yàn)來開展結(jié)構(gòu)損傷識(shí)別研究。

    1損傷識(shí)別依據(jù)

    關(guān)于彈性波理論的有關(guān)研究表明[11-12],系統(tǒng)的剪切模量G和系統(tǒng)內(nèi)剪切波速密切相關(guān),其統(tǒng)計(jì)關(guān)系為:

    其中:ρ為材料密度;vS為等效剪切波速。而彈性模量E和剪切模量G之間的關(guān)系為:

    式中:μ為泊松比。

    系統(tǒng)剛度K是彈性模量E或剪切模量G與質(zhì)量慣性矩的乘積,而在經(jīng)典動(dòng)力學(xué)中,質(zhì)量慣性矩I只跟剛體形狀、剛體自身質(zhì)量分布和慣性轉(zhuǎn)軸的位置有關(guān),與剛體的轉(zhuǎn)動(dòng)狀態(tài)無關(guān)。其表達(dá)式為:

    其中:mi為第i個(gè)質(zhì)量元的質(zhì)量;ri為第i個(gè)質(zhì)量元到慣性轉(zhuǎn)軸的垂直距離。因此體現(xiàn)系統(tǒng)運(yùn)動(dòng)狀態(tài)的剪切波速vS也直接決定著系統(tǒng)剛度。根據(jù)系統(tǒng)自振頻率ω:

    式中:K表示系統(tǒng)剛度;M表示系統(tǒng)質(zhì)量。在系統(tǒng)質(zhì)量保持不變的情況下,系統(tǒng)的等效剪切波速vS、系統(tǒng)剛度K和自振頻率ω等特性參數(shù)存在正相關(guān)性,即體系一旦受到損傷,層間位移角將增大,體系的剛度減小,等效剪切波速和自振頻率也相應(yīng)減小。隨著損傷逐漸加重,上述自振特性改變更加顯著。因此結(jié)構(gòu)模態(tài)參數(shù)、剪切波走時(shí)和層間位移角可作為結(jié)構(gòu)損傷識(shí)別參數(shù)用于結(jié)構(gòu)損傷識(shí)別。基于上述原理,開展根據(jù)振動(dòng)臺(tái)試驗(yàn)中的記錄數(shù)據(jù)分析出損傷識(shí)別參數(shù)的變化,以確定結(jié)構(gòu)損傷的程度和可能位置的深入分析。

    2振動(dòng)臺(tái)實(shí)驗(yàn)

    2.1 測(cè)試結(jié)構(gòu)概況

    建造某足尺的7層結(jié)構(gòu)的薄片進(jìn)行不同地震輸入下的振動(dòng)測(cè)試,該原型測(cè)試結(jié)構(gòu)如圖1。測(cè)試在加州大學(xué)圣地亞哥分校的高性能室外振動(dòng)臺(tái)上進(jìn)行。測(cè)試的目的在于驗(yàn)證鋼筋混凝土墻體的地震反應(yīng)。

    圖1 某足尺7層原型測(cè)試結(jié)構(gòu) Fig.1 The 7-story full-scale test structure

    該墻體中的橫向力由基于位移設(shè)計(jì)方法獲得[9]。承重腹墻是該測(cè)試墻體中的主要單元,墻體截面尺寸為3.65 m×0.20 m,可大大提高該測(cè)試結(jié)構(gòu)在振動(dòng)方向的橫向剛度。每一層設(shè)置矩形樓板,尺寸為8.15 m×3.65 m,腹墻單元和樓板四角處有鋼管立柱支撐樓板。一組翼墻和預(yù)應(yīng)力墻為該體系提供了扭轉(zhuǎn)和水平剛度,這兩組墻通過開槽和鋼支撐跟樓板相連,另外還有橫墻用于維持結(jié)構(gòu)在振動(dòng)實(shí)驗(yàn)中的穩(wěn)定性。該測(cè)試結(jié)構(gòu)的總高度約為19.20 m,總重量約為2 450 kN[13]。振動(dòng)臺(tái)和測(cè)試結(jié)構(gòu)及其周圍總共布設(shè)139個(gè)加速度計(jì),88個(gè)位移傳感器,314個(gè)應(yīng)變計(jì)和23個(gè)壓力傳感器。設(shè)置在每一層上的三個(gè)加速度計(jì)記錄結(jié)構(gòu)在振動(dòng)方向的橫向反應(yīng),固定在腹墻上的加速度計(jì)記錄每層中部的反應(yīng),樓板處正交于振動(dòng)方向的兩個(gè)加速度計(jì)記錄該方向的相應(yīng)振動(dòng),兩個(gè)放置在樓板上的加速度計(jì)用于記錄垂直振動(dòng)。所有加速度數(shù)據(jù)的采樣頻率為240 sps。

    2.2測(cè)試進(jìn)程

    比對(duì)兩組臨床免疫檢驗(yàn)的質(zhì)量以及真受控率、檢驗(yàn)環(huán)境等,采用SPSS 18.0統(tǒng)計(jì)學(xué)軟件對(duì)數(shù)據(jù)進(jìn)行處理,計(jì)量資料以“±s”表示,采用t檢驗(yàn),以P<0.05)為差異有統(tǒng)計(jì)學(xué)意義。

    實(shí)驗(yàn)中選擇一系列地震動(dòng)時(shí)程和白噪聲作為振動(dòng)方向的振動(dòng)輸入,其中天然地震動(dòng)時(shí)程為1971年圣費(fèi)南多和1999年北嶺的地震記錄。在實(shí)驗(yàn)過程中地震動(dòng)輸入峰值不斷增大,使得結(jié)構(gòu)逐漸產(chǎn)生破壞。兩次地震激發(fā)振動(dòng)以后,記錄結(jié)構(gòu)的環(huán)境振動(dòng)反應(yīng)。測(cè)試進(jìn)程的順序如表1所示。

    表 1 實(shí)驗(yàn)進(jìn)程的順序

    3測(cè)試結(jié)果分析

    基于上述的實(shí)驗(yàn)進(jìn)程,加速度計(jì)記錄了全實(shí)驗(yàn)過程中的結(jié)構(gòu)反應(yīng),這些記錄被用于開展結(jié)構(gòu)損傷識(shí)別方法分析。通過白噪聲實(shí)驗(yàn)計(jì)算出結(jié)構(gòu)的自振頻率和振型,并估算剪切波從該結(jié)構(gòu)基底傳播至頂部的時(shí)間,而層間位移角則根據(jù)四次地震激發(fā)實(shí)驗(yàn)的結(jié)構(gòu)反應(yīng)記錄數(shù)據(jù)得出。結(jié)構(gòu)損傷識(shí)別參數(shù)的變化反映出結(jié)構(gòu)的損傷程度,對(duì)于能夠用于定位損傷的參數(shù),其變化明顯的位置即為結(jié)構(gòu)損傷所在的樓層。

    3.1模態(tài)參數(shù)變化

    首先從第一次白噪聲測(cè)試數(shù)據(jù)識(shí)別出結(jié)構(gòu)的前三階模態(tài),它們只與無損傷的結(jié)構(gòu)有關(guān),為后續(xù)比較提供了一系列的基線值。地震激發(fā)實(shí)驗(yàn)以后,利用實(shí)驗(yàn)后的白噪聲數(shù)據(jù)重復(fù)計(jì)算結(jié)構(gòu)模態(tài)識(shí)別參數(shù)。規(guī)準(zhǔn)處理后的頻率識(shí)別結(jié)果列于表2中。隨著結(jié)構(gòu)損傷的不斷加重,每次地震激發(fā)實(shí)驗(yàn)后前三階模態(tài)頻率估算值逐漸減小。分析表明地震1激發(fā)后,結(jié)構(gòu)第一階模態(tài)頻率降低9%,結(jié)構(gòu)構(gòu)件的裂縫不斷發(fā)展使得結(jié)構(gòu)的有效截面剛度減小[13]。隨著裂縫更多更大,結(jié)構(gòu)的有效剛度進(jìn)一步減小,結(jié)構(gòu)自振特性變化非常明顯。地震2和3激發(fā)實(shí)驗(yàn)后,自振頻率規(guī)準(zhǔn)值減小為66%和74%,直到地震4激發(fā)實(shí)驗(yàn)后,結(jié)構(gòu)的自振頻率規(guī)準(zhǔn)值減小到49%。此時(shí)結(jié)構(gòu)的剩余剛度僅僅只有初始剛度的一半。

    圖2 第一階模態(tài)的振型曲率隨著地震激發(fā)實(shí)驗(yàn)    順序的變化Fig.2 Changes of the modal curvature with different    earthquake test procedure in the first mode

    白噪聲實(shí)驗(yàn)振型1估算頻率/Hz規(guī)準(zhǔn)頻率振型2估算頻率/Hz規(guī)準(zhǔn)頻率振型3估算頻率/Hz規(guī)準(zhǔn)頻率11.721.0010.171.0023.571.0031.560.919.840.9724.131.0251.270.748.660.8523.430.9971.140.667.910.7721.110.9090.850.495.100.5015.260.65

    模態(tài)分析結(jié)果表明自振頻率的變化可以作為識(shí)別結(jié)構(gòu)是否損傷的有效參數(shù),但不能用來定位結(jié)構(gòu)損傷。本研究采用結(jié)構(gòu)振型曲率來定位結(jié)構(gòu)損傷,以作為主導(dǎo)結(jié)構(gòu)振動(dòng)的第一階振型進(jìn)行振型曲率分析,如圖2所示。結(jié)果表明,振型曲率的變化率對(duì)損傷更敏感,且能夠準(zhǔn)確地定位損傷的位置。振型曲率的變化主要集中在一層和二層上,變化值詳見表3。變化結(jié)果與地震激發(fā)強(qiáng)度變化規(guī)律和現(xiàn)場(chǎng)觀測(cè)到的損傷程度一致。

    表 3 第一階振型在第二層處振型曲率變化

    3.2剪切波走時(shí)

    一維波動(dòng)模型可以用來分析建筑結(jié)構(gòu)的振動(dòng)特性[14]。通過對(duì)不同樓層記錄的反卷積分析,可以獲得剪切波在結(jié)構(gòu)的走時(shí)和平均等效剪切波速。剪切波在結(jié)構(gòu)某樓層中的走時(shí)主要取決于該樓層的質(zhì)量和剛度,因此根據(jù)結(jié)構(gòu)中的走時(shí)變化也可以識(shí)別結(jié)構(gòu)的損傷。結(jié)構(gòu)某位置上的剛度變化使得該樓層上的走時(shí)變化,說明結(jié)構(gòu)該位置上存在損傷[15-16]。同樣可利用每次地震激發(fā)實(shí)驗(yàn)后的白噪聲數(shù)據(jù)來計(jì)算剪切波走時(shí)。本文將所有樓層上的白噪聲數(shù)據(jù)對(duì)頂層數(shù)據(jù)進(jìn)行反卷積分析,通過12.5 s的移動(dòng)時(shí)間窗來計(jì)算走時(shí),如圖3(a)所示。

    根據(jù)反卷積干涉法構(gòu)建的波場(chǎng)如圖3(b), 波場(chǎng)中顯示了各層數(shù)據(jù)與頂層數(shù)據(jù)之間的脈沖響應(yīng)函數(shù)。地震波以脈沖方式從一層樓板進(jìn)入該結(jié)構(gòu)向上傳播至頂層,在頂層反射下行并傳播至一層。這種具有明顯的上行波和下行波的波場(chǎng)可用于計(jì)算等效剪切波速。分析表明,通過該時(shí)間窗計(jì)算出的走時(shí)與其他時(shí)間窗的計(jì)算結(jié)果基本一致,因此可基于該時(shí)間窗方法計(jì)算分層等效剪切波速及其變化。

    圖3 用于開展反卷積干涉分析和走時(shí)計(jì)算的白噪   聲速度時(shí)程及時(shí)間窗和反卷積干涉波場(chǎng)Fig.3 Velocity time history and time window of the white noise    for the deconvolution analysis and travel-times,calcula-   tion and the deconvolved interferometric waveforms

    圖4 原結(jié)構(gòu)與四次地震激發(fā)實(shí)驗(yàn)后的地震波絕    對(duì)走時(shí)及其走時(shí)變化Fig.4 The absolute travel-time of seismic wave and its change    before and after the four earthquake excitations

    白噪聲實(shí)驗(yàn)基礎(chǔ)輸入面剪切波走時(shí)絕對(duì)走時(shí)/s變化率(%)10.05726-3(地震1)0.058562.25(地震2)0.062168.67(地震3)0.0657014.79(地震4)0.0827644.5

    在結(jié)構(gòu)遭受地震作用后,不單絕對(duì)走時(shí)受到影響,走時(shí)變化也比較明顯[圖4(b)]。四次地震激發(fā)實(shí)驗(yàn)后,走時(shí)變化同樣表現(xiàn)在底部?jī)蓪由?,其中?次地震激發(fā)實(shí)驗(yàn)后基礎(chǔ)輸入面上的走時(shí)變化為0.001 3 s,而在第4次地震激發(fā)實(shí)驗(yàn)后走時(shí)變化為0.025 s,增大為前者的約20倍。由于地震作用水平的不斷增加,結(jié)構(gòu)損傷逐漸加重,走時(shí)變化同樣表現(xiàn)為逐漸增加的趨勢(shì)。剪切波走時(shí)和走時(shí)變化同樣反映出結(jié)構(gòu)在一二層上損傷較上部各層嚴(yán)重,因此可用走時(shí)和走時(shí)變化來識(shí)別結(jié)構(gòu)的損傷程度并確定結(jié)構(gòu)損傷的位置。

    3.3層間位移角

    層間位移角是結(jié)構(gòu)設(shè)計(jì)和抗震性能鑒定中的常用參數(shù)。美國聯(lián)邦法律規(guī)定了結(jié)構(gòu)在地震作用下的4種極限狀態(tài):可以使用、直接入住、生命安全和防止倒塌。根據(jù)FEMA 45014和ASCE 41-615規(guī)定,對(duì)于測(cè)試結(jié)構(gòu),其直接入住、生命安全和防止倒塌三種極限狀態(tài)下的瞬態(tài)層間位移角限值分別為0.5%、1%和2%,這些限值可直接用于判定結(jié)構(gòu)在地震作用下的損傷狀況。本研究中的層間位移角是根據(jù)實(shí)驗(yàn)中的加速度時(shí)程記錄,通過基線校正和二次積分獲得的層位移時(shí)程得到,將相連樓層的層位移時(shí)程相減除以對(duì)應(yīng)樓層的高度即可得出該層的層間位移角。取最大層間位移角進(jìn)行對(duì)比分析,4次地震激發(fā)后的結(jié)構(gòu)最大層間位移角如圖5所示。在第1次地震激發(fā)實(shí)驗(yàn)后,層間位移角最大值為0.39%,小于0.5%,結(jié)構(gòu)處于直接入住的極限狀態(tài)下。在第2和第3次地震作用下,層間位移角逐漸增大,除一層外,上部各層全部進(jìn)入生命安全的極限狀態(tài)階段,表明此時(shí)該結(jié)構(gòu)不能直接使用,需經(jīng)過一定的修補(bǔ)或加固后方可使用;結(jié)構(gòu)五層上的層間位移角最大值分別為0.91%和0.96%,接近1%,即將進(jìn)入防倒塌極限狀態(tài)。而第4個(gè)地震激發(fā)實(shí)驗(yàn)后,層間位移角全部超過了1%,完全進(jìn)入防止倒塌的極限狀態(tài),最大層間位移角為1.7%,并未超限值2%,結(jié)構(gòu)尚未倒塌?,F(xiàn)場(chǎng)檢查同樣表明第四次地震激發(fā)后測(cè)試結(jié)構(gòu)遭受了嚴(yán)重?fù)p傷,但由于其腹墻的存在使得結(jié)構(gòu)并沒有倒塌。分析表明層間位移角也可用于識(shí)別結(jié)構(gòu)的損傷程度。

    圖5 結(jié)構(gòu)不同極限狀態(tài)下的層間位移角限值和   地震激發(fā)實(shí)驗(yàn)后的結(jié)構(gòu)最大層間位移角Fig.5 Limit values inter-story drift ratio for various    limit states and the maximum ratios after 4    earthquake excitations

    4結(jié)論和展望

    本研究基于美國加州大學(xué)圣地亞哥分校的7層鋼筋混凝土剪力墻結(jié)構(gòu)的原型振動(dòng)實(shí)驗(yàn),檢驗(yàn)了三種方法的有效性。該實(shí)驗(yàn)為基于結(jié)構(gòu)地震反應(yīng)觀測(cè)臺(tái)陣記錄的結(jié)構(gòu)損傷識(shí)別研究提供了強(qiáng)大的數(shù)據(jù)庫支持和現(xiàn)場(chǎng)檢驗(yàn)對(duì)比平臺(tái)。利用結(jié)構(gòu)上安裝的結(jié)構(gòu)地震反應(yīng)觀測(cè)設(shè)備獲得結(jié)構(gòu)反應(yīng)記錄分析結(jié)構(gòu)的模態(tài)參數(shù)、剪切波走時(shí)和層間位移角等參數(shù),并利用其變化來綜合判定結(jié)構(gòu)損傷程度及其損傷位置,損傷識(shí)別的結(jié)果和實(shí)驗(yàn)數(shù)據(jù)的對(duì)比表明定性的結(jié)構(gòu)損傷描述與損傷識(shí)別方法給出的定量結(jié)果相同,分析表明上述參數(shù)可用于開展結(jié)構(gòu)損傷識(shí)別分析,并設(shè)定了對(duì)應(yīng)于不同損傷狀態(tài)下的損傷識(shí)別參數(shù)的閾值。下一步將挖掘更多的結(jié)構(gòu)損傷識(shí)別參數(shù)并檢驗(yàn)各閾值的合理性,確定更準(zhǔn)確的結(jié)構(gòu)損傷閾值,以此推進(jìn)結(jié)構(gòu)地震反應(yīng)觀測(cè)臺(tái)陣的工程應(yīng)用,為結(jié)構(gòu)抗震設(shè)計(jì)和建筑結(jié)構(gòu)的抗震性能鑒定等工作提供參考。

    致謝:感謝Ulusoy Hasan為本研究中使用的損傷識(shí)別方法提供的指導(dǎo)和幫助,同時(shí)也感謝Babak Moaveni為本研究提供的振動(dòng)臺(tái)實(shí)驗(yàn)數(shù)據(jù)。

    參考文獻(xiàn)(References)

    [1]GB50011-2010,建筑抗震設(shè)計(jì)規(guī)范[S].北京:中國建筑工業(yè)出版社, 2010.

    GB50011-2010,Code for Seismic Design of Buildings[S].Beijing: Building Industry Press of China, 2010. (in Chinese)

    [2]北京市地震局. 北京市實(shí)施《中華人民共和國防震減災(zāi)法》規(guī)定[S].北京, 2013.

    Earthquake Administration of Beijing Municipality.Provisions for Implementing Law of China on Protecting Against and Mitigating Earthquake Disasters[S].Beijing,2013. (in Chinese)

    [3]謝禮立,于雙久.強(qiáng)震觀測(cè)與分析原理[M].北京:地震出版社,1982.

    XIE Li-li,YU Shuang-jiu.Strong Motion Observation and Analyzing Theory[M].Beijing: Seismological Press,1982. (in Chinese)

    [4]李鴻晶,朱士云, Mehmet C.強(qiáng)震觀測(cè)建筑結(jié)構(gòu)的地震反應(yīng)分析[J].地震工程與工程振動(dòng),2003,23(6): 31-36.

    LI Hong-jing,ZHU Shi-yun, Mehmet C.Seismic Response Analysis of an Instrumented Building Structure[J].Earthquake Engineering and Engineering Vibration,2003,23(6):31-36.(in Chinese)

    [5]Mehmet C. Current Practice and Guidelines for USGS Instrumentation of Buildings Including Federal Buildings[R]//Prepared for 2001 Workshop on Structural Instrumentation for Consortium of Organizations for Strong-motion Observation Systems.2001.

    [6]鄭山鎖,田進(jìn),韓言召,等.考慮銹蝕的鋼結(jié)構(gòu)地震易損性分析[J].地震工程學(xué)報(bào),2014,36(1):1-6.

    ZHENG Shan-suo,TIAN Jin, HAN Yan-zhao,et al.Steel Structure Seismic Vulnerability Analysis Considering Steel Corrosion[J].China Earthquake Engineering Journal,2014,36(1):1-6. (in Chinese)

    [7]鄭山鎖,代曠宇,孫龍飛,等.鋼框架結(jié)構(gòu)的地震損傷研究[J].地震工程學(xué)報(bào),2015,37(2):290-297.

    ZHENG Shan-suo,DAI Kuang-yu,SUN Long-fei,et al.Research on the Seismic Damage of Steel Frame Structure[J].China Earthquake Engineering Journal,2015,37(2):290-297. (in Chinese)

    [8]Carden E P,Fanning P.Vibration Based Condition Monitoring: A Review[J].Structural Health Monitoring,2004,3(4):355-377.

    [9]Doebling S W,Farrar C R,Prime M B,et al.Damage Identification and Health Monitoring of Structural and Mechanical Systems from Changes in Their Vibration Characteristics: a Literature Review[R].Los Alamos National Laboratory,Los Alamos,NM,1996.

    [10]Sohn H,Farrar C R,Hemez F M,et al.A Review of Structural Health Monitoring Literature:1996-2001[R].Los Alamos National Laboratory,Los Alamos,NM,2003.

    [11]Schneider J A,Hoyos L,Mayne P W,et al.Field and laboratory Measurements of Dynamic Shear Modulus of Piedmont Residual Soils[C]//Behavioral Characteristics of Residual Soils,GSP92,American Society of Civil Engineers,Reston,VA.1999.

    [12]Suleyman K,Yalcin M.A Simple Soil-structure Interaction Model[J].Applied Mathematical Modelling,2000,24:607-635.

    [13]Panagiotou M,Restrepo J I,Conte J P.Shake-table Test of 7-Story Building Slice. Phase I:Rectangular Wall[J].Journal of Structural Engineering,2011,137(6):691-704.

    [14]Snieder R,Safak E.Extracting the Building Response Using Seismic Interferometry:Theory and Application to the Millikan Library in Pasadena, California[J].Bulletin of the Seismological Society of America,2006,96(2):586-598.

    [15]Todorovska M I,Trifunac M D.Impulse Response Analysis of the Van Nuys 7-story Hotel during 11 Earthquakes and Earthquake Damage Identification[J].Structural Control and Health Monitoring,2008,15(1):90-116.

    [16]Todorovska M I,Trifunac M D.Earthquake Damage Identification in the Imperial County Services Building-I:The Data and Time-frequency Analysis[J].Soil Dynamics and Earthquake Engineering,2007,27(6):564-576.

    Structural Damage Identification Based on Shaking Table Tests

    WANG Fei1, 2, REZA Shahneam3

    (1.InstituteofGeophysics,CEA,Beijing100081,China;2.EarthquakeAdministrationofBeijingMunicipalityBeijing100080,China; 3.U.S.GeologicalSurvey,California94025,USA)

    Abstract:In this study, we chose high-performance shaking-table tests of a full-scale seven-story reinforced-concrete shear wall structure at the University of California, San Diego to test this method with respect to damage identification. We alternately tested the structure under the excitations of white noise, the environment, and nine earthquakes. For each case, we scaled the amplitudes of the input ground motions to various levels. We recorded the acceleration responses before and after the earthquake excitations with seismometers located on the seven floors. We determined the vibration characteristics for each earthquake excitation by analyzing the acceleration responses mentioned above. These characteristics include the modal information, the shear-wave propagation characteristics, and the inter-story drift ratio. We estimated the modal frequencies of the first three modes from the recordings when white noise was first applied to the building, and considered these as the criterion. The subsequent modal frequencies were then normalized and compared with this criterion. The normalized frequencies diminished gradually with the load case tests and the normalized frequency reduced by 51 percent for the first mode. The reduction in the modal characteristics indicates that crevices develop as the amplitude of the input ground motions increase, and thereby decrease the rigidity. Lower rigidity suggests that damage throughout the building has been aggravated. However, changes in modal frequencies cannot be used to locate damage. Mode shape curvatures of the building were similarly applied to identify the building damage. Test results demonstrate that the mode shape curvatures increase significantly with the test process and the main changes were concentrated on the second floor. After the excitations of earthquakes 1, 2, 3, and 4, the curvature values were 0.214, 1.214, 7.101, and 9.641, respectively. Therefore, we conclude that the damage on the second floor was more severe. Subsequently, we used a one-dimensional shear-wave propagation model to form the virtual waveform by deconvolving the recordings on each floor with the signal on the seventh floor. This waveform has a wave equation that is identical with that of a physical waveform and reflects the propagation characteristics of the shear wave in the building. Upward traveling and downward traveling waves are recognized in the virtual waveform. The travel time of the shear waves is inferred from the upward and downward traveling waves. At the same time, we obtained the changes in the travel time. The travel time and its changes both increase with the amplitude of the input ground motions. The travel change after the earthquake-4 excitation rose by 44.5 percent on the first floor. The travel time and its changes suggest that the lower two floors were more damaged than the upper floors, and are appropriate for damage identification as well. Finally, we computed the inter-story drift ratio and compared the results with the response after the excitations of the four earthquakes. The inter-story drift ratio increases after the input ground motions and breaks through the limit values of immediate occupancy of 0.5 percent and life safety of 1.0 percent for a reinforced concrete building. After excitation by earthquake 4, the drift ratio approaches the limit value of collapse prevention of 2.0 percent. Our analysis indicates that the parameters described above are sufficient to identify the damage.

    Key words:shaking table test; damage identification; modal parameters; equivalent travel time of shear wave; inter-story drift ratio

    DOI:10.3969/j.issn.1000-0844.2016.01.0129

    中圖分類號(hào):TU317

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

    文章編號(hào):1000-0844(2016)01-0129-07

    作者簡(jiǎn)介:王飛(1979-),男,副研究員,碩士,主要從事結(jié)構(gòu)地震反應(yīng)觀測(cè)研究。E-mail: wangfei@bjseis.gov.cn。

    基金項(xiàng)目:中國地震局青年科技骨干人才項(xiàng)目(201204190030)

    收稿日期:①2015-01-05

    久久九九热精品免费| 十八禁网站网址无遮挡| 成人18禁高潮啪啪吃奶动态图| 欧美在线一区亚洲| 757午夜福利合集在线观看| 免费人妻精品一区二区三区视频| 久久久欧美国产精品| 建设人人有责人人尽责人人享有的| 天堂8中文在线网| 99久久精品国产亚洲精品| 国产福利在线免费观看视频| 亚洲精品久久成人aⅴ小说| 国产日韩欧美视频二区| 亚洲avbb在线观看| 久久久久久久久久久久大奶| 老汉色av国产亚洲站长工具| 另类精品久久| 国产aⅴ精品一区二区三区波| 国产精品av久久久久免费| 在线观看www视频免费| 午夜免费鲁丝| bbb黄色大片| 黑丝袜美女国产一区| 国产三级黄色录像| 少妇被粗大的猛进出69影院| 亚洲少妇的诱惑av| 12—13女人毛片做爰片一| 日韩免费高清中文字幕av| 久久人人爽av亚洲精品天堂| 国产高清国产精品国产三级| 黄色片一级片一级黄色片| 亚洲av国产av综合av卡| 中文字幕人妻丝袜一区二区| 精品亚洲成a人片在线观看| 黄色怎么调成土黄色| 日本vs欧美在线观看视频| 成人黄色视频免费在线看| 国产男女超爽视频在线观看| 日韩视频在线欧美| 男人舔女人的私密视频| 侵犯人妻中文字幕一二三四区| 国产亚洲一区二区精品| 欧美精品人与动牲交sv欧美| 巨乳人妻的诱惑在线观看| 国产99久久九九免费精品| 正在播放国产对白刺激| 亚洲中文日韩欧美视频| 亚洲 国产 在线| 熟女少妇亚洲综合色aaa.| 桃红色精品国产亚洲av| 搡老岳熟女国产| 亚洲熟女毛片儿| 热re99久久国产66热| 日韩大片免费观看网站| 亚洲熟女精品中文字幕| 中国美女看黄片| 一夜夜www| 亚洲精品粉嫩美女一区| 国产深夜福利视频在线观看| 亚洲五月色婷婷综合| 不卡av一区二区三区| 国产在线精品亚洲第一网站| 亚洲国产av影院在线观看| 成人黄色视频免费在线看| 一边摸一边抽搐一进一小说 | 啦啦啦视频在线资源免费观看| 精品人妻在线不人妻| 女人精品久久久久毛片| av天堂久久9| 久久精品国产99精品国产亚洲性色 | 成年人午夜在线观看视频| 免费av中文字幕在线| 亚洲国产av影院在线观看| 视频在线观看一区二区三区| 亚洲精品粉嫩美女一区| avwww免费| 999久久久国产精品视频| 一本综合久久免费| 久久影院123| 女人精品久久久久毛片| 别揉我奶头~嗯~啊~动态视频| 精品免费久久久久久久清纯 | 在线亚洲精品国产二区图片欧美| 亚洲成人手机| 免费一级毛片在线播放高清视频 | 国产亚洲一区二区精品| 久久99一区二区三区| 久久久久久久国产电影| 亚洲成国产人片在线观看| 久久久久精品人妻al黑| 亚洲熟女毛片儿| 国产精品欧美亚洲77777| 久久久国产精品麻豆| 久久久精品区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 国产激情久久老熟女| 一边摸一边抽搐一进一小说 | 久久婷婷成人综合色麻豆| 免费在线观看视频国产中文字幕亚洲| 欧美另类亚洲清纯唯美| 亚洲黑人精品在线| 久久久久精品人妻al黑| 中文欧美无线码| 国产一区有黄有色的免费视频| 亚洲av成人不卡在线观看播放网| 欧美亚洲日本最大视频资源| 亚洲专区中文字幕在线| 又黄又粗又硬又大视频| 国产精品久久久av美女十八| 国产成人av激情在线播放| 中文字幕最新亚洲高清| 熟女少妇亚洲综合色aaa.| 黄色片一级片一级黄色片| 久久亚洲真实| 狠狠婷婷综合久久久久久88av| 国产免费现黄频在线看| 亚洲成人国产一区在线观看| 欧美激情 高清一区二区三区| 欧美精品高潮呻吟av久久| 极品少妇高潮喷水抽搐| 少妇粗大呻吟视频| 老熟妇仑乱视频hdxx| 亚洲精品国产区一区二| 国产一区二区三区视频了| 丰满饥渴人妻一区二区三| 国产无遮挡羞羞视频在线观看| 午夜视频精品福利| 亚洲熟女毛片儿| 丰满饥渴人妻一区二区三| 午夜两性在线视频| 在线观看www视频免费| 亚洲中文字幕日韩| 建设人人有责人人尽责人人享有的| 成人特级黄色片久久久久久久 | 少妇粗大呻吟视频| 亚洲国产成人一精品久久久| 婷婷成人精品国产| 成人特级黄色片久久久久久久 | 免费少妇av软件| 国产男靠女视频免费网站| 自线自在国产av| 成人特级黄色片久久久久久久 | 国产野战对白在线观看| 亚洲全国av大片| 午夜福利视频在线观看免费| 亚洲伊人色综图| 12—13女人毛片做爰片一| 麻豆成人av在线观看| 在线观看66精品国产| 午夜激情av网站| 久久国产精品影院| 国产精品国产高清国产av | 男男h啪啪无遮挡| 欧美日韩精品网址| 国产亚洲av高清不卡| 美女视频免费永久观看网站| 不卡av一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 亚洲五月色婷婷综合| 在线观看www视频免费| 12—13女人毛片做爰片一| 在线观看www视频免费| 午夜福利视频在线观看免费| 久久天堂一区二区三区四区| 视频区欧美日本亚洲| videosex国产| 亚洲av第一区精品v没综合| 满18在线观看网站| 最近最新免费中文字幕在线| 亚洲一区二区三区欧美精品| 国产黄频视频在线观看| 别揉我奶头~嗯~啊~动态视频| bbb黄色大片| 夜夜夜夜夜久久久久| 男女高潮啪啪啪动态图| tocl精华| 精品国产乱码久久久久久小说| 日韩中文字幕视频在线看片| 黄色丝袜av网址大全| 国产成人啪精品午夜网站| 精品一区二区三区av网在线观看 | 超碰成人久久| 中文亚洲av片在线观看爽 | 老司机在亚洲福利影院| 欧美国产精品va在线观看不卡| 亚洲伊人久久精品综合| 久久久久久久久免费视频了| 性色av乱码一区二区三区2| 午夜免费鲁丝| 久久天躁狠狠躁夜夜2o2o| 国产不卡一卡二| 在线十欧美十亚洲十日本专区| av片东京热男人的天堂| 少妇裸体淫交视频免费看高清 | 亚洲精品一卡2卡三卡4卡5卡| 欧美国产精品一级二级三级| 久久精品91无色码中文字幕| 亚洲中文av在线| 亚洲第一av免费看| 久久久久久久大尺度免费视频| 99久久99久久久精品蜜桃| 大码成人一级视频| 国产成人精品久久二区二区91| 午夜福利一区二区在线看| 欧美人与性动交α欧美软件| 777米奇影视久久| 99精国产麻豆久久婷婷| 日本黄色视频三级网站网址 | 亚洲精品国产精品久久久不卡| 久久天躁狠狠躁夜夜2o2o| 蜜桃国产av成人99| 9色porny在线观看| 在线看a的网站| 亚洲成人国产一区在线观看| 女同久久另类99精品国产91| 成人三级做爰电影| 久久久久久久久久久久大奶| 日韩视频在线欧美| 亚洲专区国产一区二区| 成人18禁高潮啪啪吃奶动态图| 中文字幕人妻熟女乱码| 1024视频免费在线观看| 欧美精品一区二区大全| 叶爱在线成人免费视频播放| 国产麻豆69| 视频区图区小说| 欧美日韩视频精品一区| 国产黄色免费在线视频| 少妇裸体淫交视频免费看高清 | 如日韩欧美国产精品一区二区三区| 男女免费视频国产| 久久中文字幕人妻熟女| 蜜桃在线观看..| 国产精品久久久久成人av| 美女扒开内裤让男人捅视频| 亚洲成国产人片在线观看| 热99re8久久精品国产| 亚洲人成电影免费在线| 麻豆国产av国片精品| 久久精品国产亚洲av高清一级| 日韩熟女老妇一区二区性免费视频| 日韩 欧美 亚洲 中文字幕| 一二三四在线观看免费中文在| 热re99久久国产66热| 久久国产精品人妻蜜桃| av一本久久久久| 久久ye,这里只有精品| 伦理电影免费视频| 国产亚洲午夜精品一区二区久久| 国产区一区二久久| 精品一区二区三区四区五区乱码| 日本一区二区免费在线视频| 99热国产这里只有精品6| 国产成人av教育| 麻豆av在线久日| 超色免费av| 久久免费观看电影| 高清黄色对白视频在线免费看| 国产免费福利视频在线观看| 91av网站免费观看| 性色av乱码一区二区三区2| 国产伦人伦偷精品视频| 免费女性裸体啪啪无遮挡网站| 久久99一区二区三区| 国产福利在线免费观看视频| 精品免费久久久久久久清纯 | 黑人巨大精品欧美一区二区蜜桃| 精品一区二区三区视频在线观看免费 | 在线观看人妻少妇| 欧美精品高潮呻吟av久久| 欧美精品高潮呻吟av久久| 精品视频人人做人人爽| 久久久国产欧美日韩av| 欧美av亚洲av综合av国产av| 国产日韩欧美亚洲二区| 久久久久久久大尺度免费视频| 在线亚洲精品国产二区图片欧美| 999久久久国产精品视频| 亚洲精品中文字幕在线视频| 动漫黄色视频在线观看| 美女视频免费永久观看网站| 三上悠亚av全集在线观看| 亚洲精品国产区一区二| 久久久国产欧美日韩av| 国产又爽黄色视频| 国产av一区二区精品久久| 午夜福利,免费看| av有码第一页| 狠狠婷婷综合久久久久久88av| av超薄肉色丝袜交足视频| 黄片大片在线免费观看| 视频在线观看一区二区三区| 国产黄色免费在线视频| 又大又爽又粗| 日本一区二区免费在线视频| 久久精品国产99精品国产亚洲性色 | 久久久久视频综合| 精品午夜福利视频在线观看一区 | 黑人欧美特级aaaaaa片| 人人妻人人添人人爽欧美一区卜| 窝窝影院91人妻| 久久精品国产综合久久久| 国产精品1区2区在线观看. | 欧美亚洲 丝袜 人妻 在线| 大香蕉久久成人网| 一区在线观看完整版| 国产精品香港三级国产av潘金莲| 啦啦啦免费观看视频1| 精品少妇黑人巨大在线播放| videosex国产| 少妇裸体淫交视频免费看高清 | 可以免费在线观看a视频的电影网站| 首页视频小说图片口味搜索| 大片电影免费在线观看免费| 国产成人啪精品午夜网站| 人妻久久中文字幕网| √禁漫天堂资源中文www| 正在播放国产对白刺激| 日日爽夜夜爽网站| 咕卡用的链子| 九色亚洲精品在线播放| 飞空精品影院首页| 日本欧美视频一区| 久久国产精品男人的天堂亚洲| 精品少妇黑人巨大在线播放| 肉色欧美久久久久久久蜜桃| 高清毛片免费观看视频网站 | 这个男人来自地球电影免费观看| 十八禁人妻一区二区| 成年人黄色毛片网站| 国产麻豆69| 欧美国产精品一级二级三级| 一本大道久久a久久精品| 免费在线观看影片大全网站| 丝瓜视频免费看黄片| 久久国产精品男人的天堂亚洲| 天堂动漫精品| 高清黄色对白视频在线免费看| 麻豆国产av国片精品| 亚洲免费av在线视频| svipshipincom国产片| 黄色怎么调成土黄色| 国产激情久久老熟女| 精品一区二区三卡| 丝袜美腿诱惑在线| 99精国产麻豆久久婷婷| 久久香蕉激情| 久久久久久免费高清国产稀缺| 啪啪无遮挡十八禁网站| 精品乱码久久久久久99久播| 18禁美女被吸乳视频| √禁漫天堂资源中文www| 一边摸一边做爽爽视频免费| 久久天堂一区二区三区四区| 91国产中文字幕| 成人精品一区二区免费| 午夜福利视频在线观看免费| 欧美精品人与动牲交sv欧美| 久久人妻av系列| 精品国产超薄肉色丝袜足j| 美女高潮到喷水免费观看| 最新在线观看一区二区三区| 老司机在亚洲福利影院| 欧美黑人欧美精品刺激| 成年人黄色毛片网站| 啦啦啦在线免费观看视频4| 十八禁高潮呻吟视频| 欧美精品一区二区大全| 91老司机精品| 久久人人97超碰香蕉20202| 欧美日韩亚洲综合一区二区三区_| av又黄又爽大尺度在线免费看| 在线观看66精品国产| 黄片小视频在线播放| 桃花免费在线播放| 女同久久另类99精品国产91| 成年人免费黄色播放视频| 国产成人免费观看mmmm| 动漫黄色视频在线观看| 在线天堂中文资源库| 18禁黄网站禁片午夜丰满| 女性被躁到高潮视频| 久久精品国产亚洲av香蕉五月 | 丝袜美腿诱惑在线| 777米奇影视久久| 国产精品亚洲一级av第二区| av有码第一页| av一本久久久久| 丰满饥渴人妻一区二区三| 高清视频免费观看一区二区| 午夜老司机福利片| 精品国产国语对白av| 97在线人人人人妻| 亚洲av第一区精品v没综合| 久久久精品94久久精品| 国产精品99久久99久久久不卡| 91精品国产国语对白视频| 婷婷成人精品国产| 99精国产麻豆久久婷婷| bbb黄色大片| 91大片在线观看| 国产在线一区二区三区精| 国产免费现黄频在线看| 精品一区二区三区四区五区乱码| 午夜免费鲁丝| 飞空精品影院首页| 老汉色av国产亚洲站长工具| 久久久国产精品麻豆| av福利片在线| 免费在线观看视频国产中文字幕亚洲| 两个人免费观看高清视频| 亚洲天堂av无毛| 色播在线永久视频| 无遮挡黄片免费观看| 国产成人影院久久av| 露出奶头的视频| 国产av又大| 精品少妇内射三级| 十八禁高潮呻吟视频| 一区二区三区精品91| 亚洲午夜理论影院| 久久久精品免费免费高清| 久久久国产成人免费| 757午夜福利合集在线观看| 人妻一区二区av| 夜夜夜夜夜久久久久| 动漫黄色视频在线观看| 国产有黄有色有爽视频| 成人18禁高潮啪啪吃奶动态图| av天堂在线播放| 超色免费av| 色婷婷久久久亚洲欧美| 国产在线视频一区二区| 亚洲国产欧美在线一区| 国产无遮挡羞羞视频在线观看| 性色av乱码一区二区三区2| 国产精品免费大片| 国产无遮挡羞羞视频在线观看| 亚洲人成77777在线视频| tube8黄色片| 亚洲色图 男人天堂 中文字幕| 欧美人与性动交α欧美软件| www.999成人在线观看| 老司机午夜十八禁免费视频| 久热爱精品视频在线9| 露出奶头的视频| 青青草视频在线视频观看| 国产精品久久久av美女十八| 欧美日韩视频精品一区| 国产老妇伦熟女老妇高清| 国产一区二区三区综合在线观看| 国产亚洲欧美精品永久| 99久久人妻综合| 黄片大片在线免费观看| 免费在线观看视频国产中文字幕亚洲| 十八禁网站网址无遮挡| 多毛熟女@视频| 极品少妇高潮喷水抽搐| 久久久国产成人免费| 欧美精品亚洲一区二区| 久久人妻福利社区极品人妻图片| 最近最新中文字幕大全电影3 | 久久中文字幕一级| 91老司机精品| av在线播放免费不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美精品综合一区二区三区| 91av网站免费观看| 日本av手机在线免费观看| 亚洲第一av免费看| 美女主播在线视频| 一边摸一边做爽爽视频免费| 肉色欧美久久久久久久蜜桃| 免费在线观看日本一区| 男女无遮挡免费网站观看| 午夜老司机福利片| 超碰97精品在线观看| 欧美精品人与动牲交sv欧美| 欧美乱码精品一区二区三区| 日本精品一区二区三区蜜桃| 国产精品 欧美亚洲| 十八禁高潮呻吟视频| 亚洲天堂av无毛| 亚洲精品成人av观看孕妇| 18禁美女被吸乳视频| 黑人操中国人逼视频| 99国产精品免费福利视频| 免费av中文字幕在线| 国产日韩欧美在线精品| 久久久精品国产亚洲av高清涩受| 日本黄色视频三级网站网址 | 69av精品久久久久久 | 亚洲精品美女久久av网站| 一级毛片电影观看| 丁香欧美五月| 亚洲国产av影院在线观看| 中文字幕制服av| 嫁个100分男人电影在线观看| 午夜激情av网站| 免费女性裸体啪啪无遮挡网站| netflix在线观看网站| 91精品三级在线观看| 免费久久久久久久精品成人欧美视频| 又紧又爽又黄一区二区| 日本精品一区二区三区蜜桃| 99久久国产精品久久久| 欧美性长视频在线观看| 亚洲视频免费观看视频| 男人操女人黄网站| 亚洲av成人一区二区三| 欧美黑人精品巨大| 欧美精品一区二区大全| 精品第一国产精品| 日韩欧美国产一区二区入口| 狠狠婷婷综合久久久久久88av| 男女高潮啪啪啪动态图| 亚洲中文字幕日韩| 天天添夜夜摸| 国产精品 欧美亚洲| 妹子高潮喷水视频| 欧美日韩亚洲综合一区二区三区_| 法律面前人人平等表现在哪些方面| 一边摸一边抽搐一进一出视频| 男女床上黄色一级片免费看| av片东京热男人的天堂| 久热爱精品视频在线9| 19禁男女啪啪无遮挡网站| 日韩中文字幕视频在线看片| 中文亚洲av片在线观看爽 | 男女午夜视频在线观看| 久久这里只有精品19| 成人av一区二区三区在线看| 欧美成狂野欧美在线观看| 国产日韩欧美视频二区| 十八禁人妻一区二区| 亚洲欧美激情在线| 欧美精品一区二区大全| 国产精品久久久久久精品电影小说| 男女之事视频高清在线观看| 一区二区三区乱码不卡18| videosex国产| 大片免费播放器 马上看| 老汉色av国产亚洲站长工具| 色视频在线一区二区三区| 美女扒开内裤让男人捅视频| 巨乳人妻的诱惑在线观看| 亚洲avbb在线观看| 国产淫语在线视频| 国产激情久久老熟女| 人人妻,人人澡人人爽秒播| 99久久国产精品久久久| 日韩欧美国产一区二区入口| 视频区图区小说| 丝瓜视频免费看黄片| 色婷婷av一区二区三区视频| av天堂在线播放| e午夜精品久久久久久久| 亚洲成人免费av在线播放| 一本色道久久久久久精品综合| 波多野结衣av一区二区av| 国产97色在线日韩免费| 国产日韩欧美亚洲二区| 在线av久久热| 欧美在线黄色| 色尼玛亚洲综合影院| 国产精品久久久人人做人人爽| 成人影院久久| 99国产精品免费福利视频| 亚洲 欧美一区二区三区| 久久久久精品人妻al黑| 啦啦啦视频在线资源免费观看| 久久久国产精品麻豆| 搡老熟女国产l中国老女人| 亚洲欧美一区二区三区黑人| 757午夜福利合集在线观看| 国产成人精品无人区| 一级a爱视频在线免费观看| 桃红色精品国产亚洲av| 亚洲成人免费电影在线观看| 亚洲五月婷婷丁香| 国产高清视频在线播放一区| 久久这里只有精品19| videos熟女内射| 免费不卡黄色视频| 在线观看一区二区三区激情| 老司机午夜福利在线观看视频 | 欧美日韩亚洲综合一区二区三区_| 欧美人与性动交α欧美软件| 69av精品久久久久久 | 国产精品久久久久久人妻精品电影 | 麻豆成人av在线观看| 欧美精品一区二区免费开放| 亚洲精品国产区一区二| 欧美日韩一级在线毛片| 国产欧美日韩一区二区三区在线| 性少妇av在线| 国产又爽黄色视频| 免费av中文字幕在线| 一本一本久久a久久精品综合妖精| 精品少妇久久久久久888优播| 国产一区二区激情短视频| 亚洲精品av麻豆狂野| 叶爱在线成人免费视频播放| 两性午夜刺激爽爽歪歪视频在线观看 | cao死你这个sao货| 麻豆成人av在线观看| 麻豆国产av国片精品| 99国产极品粉嫩在线观看| 成人国语在线视频| 亚洲国产看品久久| 高清在线国产一区| 国产一区二区三区综合在线观看| 丝袜美足系列| 51午夜福利影视在线观看| 精品福利观看|