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

    黃土地區(qū)地鐵車站地震反應(yīng)的頻域分析及空間效應(yīng)

    2016-12-15 10:42:27權(quán)登州王毅紅馬蓬渤井彥林
    振動與沖擊 2016年21期
    關(guān)鍵詞:傅里葉震動峰值

    權(quán)登州 , 王毅紅 , 馬蓬渤 , 井彥林 , 陳 蘇

    (1.長安大學(xué) 建筑工程學(xué)院,西安 710061; 2.中國地震局地球物理研究所,北京 100081)

    ?

    黃土地區(qū)地鐵車站地震反應(yīng)的頻域分析及空間效應(yīng)

    權(quán)登州1, 王毅紅1, 馬蓬渤1, 井彥林1, 陳 蘇2

    (1.長安大學(xué) 建筑工程學(xué)院,西安 710061; 2.中國地震局地球物理研究所,北京 100081)

    在黃土場地與地鐵地下結(jié)構(gòu)動力相互作用振動臺試驗(yàn)中,測得模型地基及地鐵車站的加速度反應(yīng)及應(yīng)變反應(yīng)。基于實(shí)測數(shù)據(jù)對模型地基、地鐵車站的加速度反應(yīng)進(jìn)行了頻域分析,對比不同觀測斷面內(nèi)地鐵車站的地震反應(yīng),分析地鐵車站加速度反應(yīng)的空間效應(yīng)。結(jié)果表明:隨輸入峰值加速度增加,模型體系加速度反應(yīng)傅里葉譜基頻減小,相應(yīng)譜值增大;隨埋深減小地基加速度反應(yīng)傅里葉譜低頻成分增大,高頻成分減小,且增大部分頻帶寬度隨輸入峰值加速度增加而縮小;隨輸入峰值加速度增加,結(jié)構(gòu)加速度反應(yīng)傅里葉譜主頻集中、變窄,高頻成分呈現(xiàn)增強(qiáng)效應(yīng);模型體系加速度反應(yīng)傅里葉譜主頻帶寬呈現(xiàn)出松潘波作用時(shí)較寬,西安人工波作用時(shí)居中,Taft波作用時(shí)較窄;受結(jié)構(gòu)端頭效應(yīng)及土結(jié)相互作用中傾斜與扭轉(zhuǎn)影響,地鐵車站地震反應(yīng)具有顯著的空間效應(yīng)。研究結(jié)論可為黃土地區(qū)地鐵地下結(jié)構(gòu)的抗震設(shè)計(jì)及相關(guān)理論研究提供重要參考。

    黃土地區(qū);地鐵車站;振動臺試驗(yàn);加速度反應(yīng);頻域分析;空間效應(yīng)

    近年,隨著我國西部地區(qū)城市地下空間開發(fā)利用的加快發(fā)展,黃土地區(qū)地鐵工程日益增多。如西安地鐵已建成通車50.3 km,遠(yuǎn)景規(guī)劃總里程達(dá)660 km,蘭州地鐵一號線已開工建設(shè),計(jì)劃2016年建成通車,其遠(yuǎn)景規(guī)劃里程約130 km[1]。由于黃土具有柱狀節(jié)理、大孔隙與弱膠結(jié)的特殊結(jié)構(gòu)性及對水的特殊敏感性[2],使得黃土地區(qū)震害相對其它土類要嚴(yán)重得多。動力、靜力及浸水作用分別引起黃土震陷、壓密及濕陷變形,黃土震陷變形大小與動力荷載相對于靜力荷載和浸水作用的先后有密切關(guān)系。因此,黃土場地對地震作用十分敏感,對黃土場地條件下地鐵地下結(jié)構(gòu)地震動力反應(yīng)特征進(jìn)行研究具有重要意義。

    1995年日本阪神大地震發(fā)生以來,國內(nèi)外學(xué)者對地下結(jié)構(gòu)地震反應(yīng)特性開展了廣泛深入的研究,JAFARZADEH等[3]進(jìn)行了系列振動臺試驗(yàn),研究了埋地管線在均質(zhì)和非均質(zhì)場地中的地震反應(yīng),AMOROSI等[4]對圓形隧道的地震響應(yīng)進(jìn)行了數(shù)值分析,AZADI等[5]基于FLAC 3D軟件研究了可液化地基隧道的地震反應(yīng),CHIAN[6]研究了振動臺試驗(yàn)中不同邊界條件對模型結(jié)構(gòu)上浮特性的影響,楊林德等[7]進(jìn)行了軟土地層中地鐵車站結(jié)構(gòu)及地鐵車站與區(qū)間隧道接頭結(jié)構(gòu)的振動臺模型試驗(yàn);陳國興等[8-10]在軟弱場地、可液化土層等不同地質(zhì)條件下進(jìn)行了土-地鐵地下結(jié)構(gòu)動力相互作用大型振動臺試驗(yàn)和數(shù)值模擬研究;劉晶波等[11]對砂土地基中地下結(jié)構(gòu)進(jìn)行了離心機(jī)振動臺試驗(yàn)研究;門玉明等[12]對地裂縫場地中地鐵區(qū)間隧道的地震反應(yīng)進(jìn)行了研究。除此之外,其他學(xué)者也對地下結(jié)構(gòu)地震反應(yīng)進(jìn)行了有意義的探索[13-14]。然而,針對黃土場地中地鐵地下結(jié)構(gòu)地震動力反應(yīng)的研究較少。

    我國黃土主要分布在西北、華北等地區(qū),西安地鐵是我國首次在黃土地區(qū)修建地鐵工程。本文以西安地鐵工程為背景,進(jìn)行了黃土與地鐵地下車站結(jié)構(gòu)動力相互作用地震模擬振動臺試驗(yàn),測得了模型地基及地鐵車站的加速度及應(yīng)變反應(yīng)?;趯?shí)測數(shù)據(jù)對模型地基、地鐵車站的加速度反應(yīng)傅里葉譜進(jìn)行了分析,從頻域分析角度對模型體系地震動力反應(yīng)特性進(jìn)行了研究;同時(shí),對比不同觀測斷面內(nèi)地鐵車站的地震反應(yīng),分析了地鐵車站地震反應(yīng)的空間效應(yīng),研究結(jié)論可為黃土地區(qū)地鐵地下車站、區(qū)間隧道及地下商業(yè)街等地下結(jié)構(gòu)抗震設(shè)計(jì)及相關(guān)理論研究提供可靠資料。

    1 試驗(yàn)簡介

    1.1 模型相似設(shè)計(jì)

    本試驗(yàn)主要研究黃土場地中地鐵車站的抗震機(jī)理及土-結(jié)構(gòu)動力相互作用機(jī)理,由于地鐵車站模型尺寸較小,通過附加全配重的人工質(zhì)量模型消除重力失真效應(yīng)很難實(shí)現(xiàn),因此采用欠人工質(zhì)量模型。根據(jù)Bockingham π定理,選取長度、彈性模量及加速度為基本物理量,并充分考慮試驗(yàn)系統(tǒng)的臺面尺寸、動力性能、承載噸位及其它配套設(shè)備性能等制約,確定模型結(jié)構(gòu)的相似關(guān)系見表1,試驗(yàn)中采用在模型結(jié)構(gòu)中粘貼鉛塊的方法實(shí)現(xiàn)附加人工質(zhì)量。

    1.2 模型制作

    以西安地鐵四號線飛天路車站為原型結(jié)構(gòu),該車站位于黃土臺塬區(qū),為典型的兩層雙跨矩形斷面結(jié)構(gòu)。采用微粒混凝土和鍍鋅鋼絲制作地鐵車站模型,根據(jù)相似設(shè)計(jì)確定模型結(jié)構(gòu)尺寸,按照原型與模型結(jié)構(gòu)構(gòu)件彎矩、剪力等效及施工可操作性原則,確定模型結(jié)構(gòu)中鋼筋布置。制作模型結(jié)構(gòu)時(shí),使用木板制作澆筑結(jié)構(gòu)的外圍邊界,內(nèi)部構(gòu)件采用在聚苯乙烯塑料泡沫實(shí)體中挖出構(gòu)件澆筑空間的方法制作,制作完成的車站模型如圖1所示。試驗(yàn)使用的疊層剪切模型土箱凈尺寸為3.5 m(振動方向)×2.0 m(縱向)×1.7 m(豎向),采用15層疊層方鋼管框架并輔之以雙側(cè)面鋼板約束,在水平振動方向形成可自由滑動的邊界,如圖2所示。該疊層剪切模型土箱能較好地消除邊界上地震波的反射或散射效應(yīng)[15]。試驗(yàn)中黃土取自西安地鐵四號線飛天路車站施工基坑,制備模型地基時(shí),將黃土分層裝入模型土箱,按原型場地土天然含水量和天然密度控制模型地基含水量及密度。

    表1 模型相似常數(shù)

    圖1 制作完成的模型結(jié)構(gòu)Fig.1 Model structure

    圖2 疊層剪切模型土箱Fig.2 Laminar shear container

    1.3 傳感器布置

    基于ABAQUS大型非線性有限元分析平臺,采用有限元-無限元耦合的建模方法,對黃土場地中地鐵車站的地震響應(yīng)進(jìn)行數(shù)值模擬;根據(jù)數(shù)值模擬結(jié)果、模型尺寸及試驗(yàn)條件,確定試驗(yàn)傳感器的布置。本試驗(yàn)中加速度傳感器及布置方案如圖3、圖4所示。圖4中A為加速度傳感器,共26個(gè),其中A20~A22設(shè)于2號次觀測斷面中與A3~A5相對應(yīng)的位置,A23設(shè)于3號次觀測斷面中與A5相對應(yīng)的位置;S為應(yīng)變片,共44個(gè)。

    圖3 加速度傳感器Fig.3 Accelerometer in the test

    圖4 黃土地區(qū)地鐵車站振動臺試驗(yàn)傳感器布置Fig.4 Arrangement plan of sensors for shaking table test

    1.4 試驗(yàn)加載方案

    本試驗(yàn)采用的激振系統(tǒng)為美國MTS公司生產(chǎn)的水平單向高性能地震模擬振動臺,臺面尺寸為3.36 m× 4.86 m,最大載重量25T,最大加速度±1.0 g。為了研究不同特性地震動作用下黃土場地地鐵車站的地震反應(yīng),分別選擇Taft波、松潘波和西安人工波作為輸入地震動。Taft地震波距震中約43.5 km,記錄長達(dá)54 s,峰值加速度175.9 cm/s2;松潘波距震中約122 km,峰值加速度為40.2 cm/s2,強(qiáng)震部分持續(xù)約120 s;西安人工波是根據(jù)西安地鐵飛天路車站黃土場地特性,按概率方法人工合成的地震波,重現(xiàn)期為475年的基巖地震動峰值加速度為200.6 cm/s2,強(qiáng)震持時(shí)為60 s。本次振動臺試驗(yàn)輸入的地震動如圖5所示,加載工況見表2。

    圖5 輸入地震動的加速度時(shí)程及傅里葉譜Fig.5 Time-histories and Fourier spectra of input accelerations

    編號工況地震動加速度峰值/g編號工況地震動加速度峰值/g編號工況地震動加速度峰值/g1B0白噪聲0.0510S3松潘波0.219T5Taft波0.62S1松潘波0.0511T3Taft波0.220X5西安人工波0.63T1Taft波0.0512X3西安人工波0.221B5白噪聲0.054X1西安人工波0.0513B3白噪聲0.0522S6松潘波0.85B1白噪聲0.0514S4松潘波0.423T6Taft波0.86S2松潘波0.115T4Taft波0.424X6西安人工波0.87T2Taft波0.116X4西安人工波0.425X7西安人工波1.28X2西安人工波0.117B4白噪聲0.0526B7白噪聲0.059B2白噪聲0.0518S5松潘波0.6---

    2 模型地基加速度反應(yīng)的峰值及頻域分析

    2.1 模型地基峰值加速度

    在模型地基中由下而上分別布置加速度計(jì)A11、A12、A13和A14以測試模型地基的地震動力反應(yīng)特性。試驗(yàn)記錄的加速度數(shù)據(jù)表明,模型地基內(nèi)各測點(diǎn)的加速度反應(yīng)時(shí)程曲線形狀特征均與輸入地震波相似,強(qiáng)震持時(shí)基本相同,峰值加速度反應(yīng)呈現(xiàn)出一定的規(guī)律。輸入不同峰值加速度時(shí),松潘波、Taft波和西安人工波作用下模型地基中不同深度處峰值加速度見表3。由表可知,相同類型地震動作用下,模型地基內(nèi)各測點(diǎn)的峰值加速度均隨輸入峰值加速度的增加而增大;輸入峰值加速度相同時(shí),模型地基內(nèi)各測點(diǎn)表現(xiàn)為西安人工波作用下峰值加速度顯著大于松潘波和Taft波;除個(gè)別工況外,模型地基內(nèi)各測點(diǎn)的峰值加速度均隨著埋深的減小而增大。

    2.2 模型地基加速度反應(yīng)的頻域分析

    輸入不同峰值加速度的Taft波時(shí),模型地基加速度反應(yīng)的傅里葉譜如圖6所示(圖中字母和數(shù)字組合代表工況或測點(diǎn),下同)。由圖可知,地表測點(diǎn)A14加速度反應(yīng)傅里葉譜發(fā)生了明顯變化:隨輸入峰值加速度的增加,地表加速度反應(yīng)傅里葉譜的基頻逐漸減小,基頻對應(yīng)的傅里葉譜值逐漸增大。在0.1 g地震動作用下,A14加速度反應(yīng)傅里葉譜基頻為7.15 Hz,相應(yīng)譜值為0.734 cm/s;0.4 g地震動作用下,其基頻和譜值分別為4.81 Hz和3.891 cm/s;0.8 g地震動作用下,其基頻和譜值分別為2.39 Hz和7.787 cm/s。這是由于隨輸入峰值加速度的增加,地基加速度反應(yīng)逐漸增大,土體非線性殘余變形增大,剛度不斷降低,從而使其基頻逐漸減?。惠斎敕逯导铀俣仍黾?,地震動中蘊(yùn)含的能量增大,因此基頻處的傅里葉譜值逐漸增大。同時(shí),隨輸入峰值加速度的增加,地表測點(diǎn)A14加速度反應(yīng)傅里葉譜在0~10 Hz低頻部分逐漸增大,在10~30 Hz高頻部分逐漸較小。這是由于土體基頻較低,輸入地震動的低頻成分在共振作用下被放大,而高頻成分被濾掉,隨著輸入峰值加速度的增加,土體發(fā)生非線性變形而軟化,剛度降低,地基低頻放大和高頻過濾效應(yīng)更顯著。模型地基內(nèi)不同深度處加速度反應(yīng)傅里葉譜的變化特征與地表測點(diǎn)相似;松潘波和西安人工波作用下,模型地基內(nèi)同深度處加速度反應(yīng)傅里葉譜的變化與輸入Taft波時(shí)的規(guī)律相似。

    表3 模型地基不同深度處的峰值加速度

    西安人工波作用下,模型地基中不相同深度處加速度反應(yīng)傅里葉譜如圖7所示。由圖可知,輸入相同峰值加速度時(shí),模型地基不同深度處加速度反應(yīng)傅里葉譜具有不同的特征。在0.1 g地震動作用下,隨著土層深度減小模型地基加速度反應(yīng)傅里葉譜在0~15 Hz低頻部分逐漸增大,在15~30 Hz高頻部分先減小,到地表時(shí)又增大;在0.8 g地震動作用下,隨著土層深度減小,模型地基加速度反應(yīng)傅里葉譜在0~5 Hz低頻部分逐漸增大,在5~30 Hz高頻部分逐漸減小。由于隨著土層深度減小,土體加速度反應(yīng)增大,層間變形增大,殘余變形累計(jì)效應(yīng)增強(qiáng),模型地基的剛度退化,對輸入地震波的低頻放大和高頻過濾效應(yīng)更顯著,從而加速度反應(yīng)傅里葉譜低頻部分逐漸增大,高頻部分逐漸減小。輸入0.1 g地震動時(shí),模型地基加速度反應(yīng)傅里葉譜的高頻成分到達(dá)地基表面時(shí)又表現(xiàn)出一定程度的增強(qiáng),這可能與地震波傳播至地基表面時(shí)發(fā)生了反射和散射效應(yīng)有關(guān)。當(dāng)輸入峰值加速度較小時(shí),模型地基加速度反應(yīng)較小,地表處地震波的反射和散射效應(yīng)對加速度反應(yīng)高頻成分增強(qiáng)作用較顯著;當(dāng)輸入峰值加速度較大時(shí),模型地基加速度反應(yīng)增大,地表處地震波的反射和散射效應(yīng)對加速度反應(yīng)影響程度有限, 因此地基加速度反應(yīng)傅里葉譜高頻成分到達(dá)地表時(shí)未出現(xiàn)增大現(xiàn)象。同時(shí),隨著輸入峰值加速度增加,土體大范圍軟化整體剛度逐漸降低,模型地基加速度反應(yīng)傅里葉譜中增強(qiáng)的低頻成分頻帶寬度逐漸縮小。松潘波和Taft波作用下,隨著土層深度的減小模型地基加速度反應(yīng)傅里葉譜變化特征與西安人工波作用時(shí)相似。此外,在0.4 g西安人工波作用下,模型地基不同深度處各測點(diǎn)的加速度反應(yīng)傅里葉譜均在4 Hz位置出現(xiàn)了集中、放大現(xiàn)象。這是由于強(qiáng)震作用下,土體內(nèi)發(fā)生了較大幅度的塑性變形,模型地基剛度衰減顯著,加速度反應(yīng)傅里葉譜出現(xiàn)單一頻率的集中與放大現(xiàn)象。

    圖6 不同峰值加速度的Taft波作用下地基加速度反應(yīng)的傅里葉譜Fig.6 Fourier spectrum of accelerations in soil under Taft wave of different PGA

    圖7 西安人工波作用下地基不同深度處加速度反應(yīng)的傅里葉譜Fig.7 Fourier spectrum of accelerations in different soil depth under Xi’an artificial wave

    輸入0.8 g峰值加速度時(shí),不同頻譜特性地震動作用下模型地基加速度反應(yīng)的傅里葉譜如圖8所示。由圖可知,模型地基中不同深度處加速度反應(yīng)傅里葉譜主頻帶寬具有顯著差異。由于松潘波的主頻帶較寬,西安人工波主頻帶寬居中,Taft波主頻帶較窄,因此,模型地基中不同深度處加速度反應(yīng)傅里葉譜的主頻帶寬也呈現(xiàn)出松潘波作用時(shí)較寬,西安人工波作用時(shí)居中,Taft波作用時(shí)較窄的現(xiàn)象。不同地震動峰值加速度作用下,模型地基中各深度處加速度反應(yīng)傅里葉譜在輸入不同頻譜特性地震動時(shí)具有相似的規(guī)律。同時(shí),在0.8 g松潘波作用下,模型地基中A11測點(diǎn)的加速度反應(yīng)傅里葉譜呈現(xiàn)出雙峰現(xiàn)象,這可能與模型地基中土層在強(qiáng)震持時(shí)內(nèi)產(chǎn)生了較大幅度的塑性殘余變形有關(guān)。

    圖8 不同頻譜特性地震動作用下地基加速度反應(yīng)的傅里葉譜Fig.8 Fourier spectrum of accelerations in soil under different types of ground motions

    3 地鐵車站加速度反應(yīng)的峰值及頻域分析

    3.1 地鐵車站峰值加速度

    在地鐵車站底板、中板及頂板高度處分別布置加速度傳感器A20、A21和A22,以測試結(jié)構(gòu)不同部位的加速度反應(yīng)。振動臺試驗(yàn)表明,地鐵車站內(nèi)各測點(diǎn)的加速度反應(yīng)時(shí)程曲線形狀均與輸入地震波相似,強(qiáng)震持時(shí)基本相同,峰值加速度反應(yīng)呈現(xiàn)出一定的規(guī)律。各工況下地鐵車站不同高度處峰值加速度見表4。由表可知,相同類型地震動作用下,地鐵車站內(nèi)各測點(diǎn)的峰值加速度隨輸入峰值加速度的增加而增大;輸入峰值加速度相同時(shí),西安人工波作用下地鐵車站的峰值加速度大于松潘波和Taft波;隨著輸入峰值加速度及結(jié)構(gòu)高度的增加,不同頻譜特性地震動作用下地鐵車站峰值加速度的差異更加顯著;當(dāng)輸入峰值加速度小于0.4 g時(shí),地鐵車站內(nèi)各測點(diǎn)的峰值加速度隨高度的增加變化幅度較小,輸入峰值加速度大于0.4 g時(shí),地鐵車站內(nèi)峰值加速度均呈現(xiàn)出頂部和底部較大,中部相對較小,且結(jié)構(gòu)頂部峰值加速度略微小于底部。地鐵車站不同高度處峰值加速度分布與土結(jié)動力相互作用及結(jié)構(gòu)自身質(zhì)量、剛度分布有關(guān)。

    表4 地鐵車站結(jié)構(gòu)高度處的峰值加速度

    3.2 地鐵車站加速度反應(yīng)頻域分析

    輸入不同峰值加速度的Taft波時(shí),地鐵車站內(nèi)不同測點(diǎn)加速度反應(yīng)傅里葉譜如圖9所示。由圖可知,隨著輸入峰值加速度的增加,結(jié)構(gòu)底部A20加速度反應(yīng)傅里葉譜的基頻逐漸減小,基頻對應(yīng)譜值逐漸增大。在0.1 g地震動作用下,A20加速度反應(yīng)傅里葉譜基頻為7.59 Hz,相應(yīng)譜值為0.675 cm/s;0.4 g地震動作用下,其基頻和譜值分別為4.81 Hz和3.460 cm/s;0.8 g地震動作用下,其基頻和譜值分別為2.84 Hz和7.604 cm/s。隨著輸入峰值加速度的增加,土體非線性殘余變形增大,地基剛度降低,同時(shí)結(jié)構(gòu)內(nèi)部產(chǎn)生微裂縫并逐漸發(fā)展,損傷變形不斷積累,其固有頻率降低,從而結(jié)構(gòu)加速度反應(yīng)傅里葉譜的基頻逐漸減?。浑S著輸入峰值加速度的增加,地震動中蘊(yùn)含的能量增大,因此基頻對應(yīng)的傅里葉譜值逐漸增大。此外,隨著輸入峰值加速度的增加,結(jié)構(gòu)內(nèi)各測點(diǎn)的加速度反應(yīng)傅里葉譜主頻帶寬度逐漸變窄,高頻部分呈現(xiàn)出一定的增強(qiáng)效果,且隨高度增加,傅里葉譜高頻部分增強(qiáng)效果更加顯著。這是由于隨著輸入峰值加速度的增加,結(jié)構(gòu)產(chǎn)生了損傷并不斷累積,固有頻率降低,對地震動主頻帶中與其固有頻率相近成分產(chǎn)生了放大效應(yīng),對固有頻率兩側(cè)的頻率成分具有一定的抑制作用,因此結(jié)構(gòu)加速度反應(yīng)傅里葉譜的主頻集中,且?guī)捵冋煌瑫r(shí),隨著輸入峰值加速度的增加,土結(jié)構(gòu)動力相互作用逐漸劇烈,地鐵車站與周圍土體之間出現(xiàn)偶發(fā)性的脫離與碰撞現(xiàn)象,土體撞擊地鐵車站時(shí)產(chǎn)生了高頻振動效應(yīng),因此結(jié)構(gòu)加速度反應(yīng)傅里葉譜高頻部分呈現(xiàn)出一定的增強(qiáng)效果;隨著高度增加,土結(jié)動力相互作用增強(qiáng),結(jié)構(gòu)受土體撞擊產(chǎn)生的高頻振動效應(yīng)更明顯,因此結(jié)構(gòu)加速度反應(yīng)傅里葉譜高頻部分增強(qiáng)效果更加顯著。結(jié)構(gòu)內(nèi)不同高度處加速度傅里葉譜的變化特征與底部測點(diǎn)類似;松潘波和西安人工波作用時(shí),結(jié)構(gòu)加速度傅里葉譜的變化規(guī)律與輸入Taft波時(shí)相似。

    圖9 不同峰值加速度的Taft波作用下地鐵車站加速度反應(yīng)的傅里葉譜Fig.9 Fourier spectrum of accelerations in station under Taft wave of different PGA

    西安人工波作用下,地鐵車站內(nèi)不同高度處加速度反應(yīng)傅里葉譜如圖10所示。由圖可知,輸入相同峰值加速度時(shí),結(jié)構(gòu)內(nèi)不同高度處加速度反應(yīng)傅里葉譜具有不同的特征。地震動較小時(shí),隨著高度增加結(jié)構(gòu)加速度反應(yīng)傅里葉譜呈現(xiàn)出低頻成分先減小后增大,高頻成分變化規(guī)律不明顯;地震動較大時(shí),隨著高度增加結(jié)構(gòu)加速度反應(yīng)傅里葉譜呈現(xiàn)出低頻成分先減小后增大,而高頻成分顯著增大。在0.1 g地震動作用下,隨著高度增加,結(jié)構(gòu)加速度反應(yīng)傅里葉譜在0~8 Hz低頻部分表現(xiàn)為先減小后增大,在8~30 Hz高頻部分變化規(guī)律不明顯;在0.8 g 地震動作用下,隨著高度增加,結(jié)構(gòu)加速度反應(yīng)傅里葉譜在0~8 Hz低頻部分亦呈現(xiàn)先減小后增大,而在8~30 Hz高頻部分則顯著增大。隨著高度增大,結(jié)構(gòu)內(nèi)各測點(diǎn)加速度反應(yīng)傅里葉譜低頻成分均為先減小后增大,與地鐵車站峰值加速度反應(yīng)變化規(guī)律一致,這可能與結(jié)構(gòu)自身質(zhì)量及剛度分布特征有關(guān)。地震動較小時(shí),由于土結(jié)構(gòu)動力相互作用較弱,結(jié)構(gòu)內(nèi)不同高度加速度反應(yīng)高頻成分變化規(guī)律不明顯;隨著輸入峰值加速度的增大,土與結(jié)構(gòu)動力相互作用逐漸增強(qiáng);地震動較大時(shí),土結(jié)構(gòu)動力相互作用劇烈,結(jié)構(gòu)受周圍土體脫離與撞擊作用增強(qiáng),從而產(chǎn)生明顯的高頻振動效應(yīng),且隨著高度增加高頻振動效應(yīng)更明顯,因此,結(jié)構(gòu)加速度反應(yīng)高頻成分增大效應(yīng)隨著高度增加更加顯著。松潘波和Taft波作用下,隨著高度增加結(jié)構(gòu)加速度反應(yīng)傅里葉譜的變化規(guī)律與西安人工波作用時(shí)相似。此外,在0.4 g西安人工波作用下,結(jié)構(gòu)不同高度處加速度反應(yīng)傅里葉譜均在4 Hz位置出現(xiàn)了集中、放大現(xiàn)象,與該工況下模型地基中加速度反應(yīng)傅里葉譜具有相似的特征。這是由于較強(qiáng)地震動作用下,結(jié)構(gòu)周圍土體內(nèi)發(fā)生了較大幅度的塑性變形,剛度衰減顯著,從而使結(jié)構(gòu)與模型地基加速度反應(yīng)傅里葉譜均出現(xiàn)了單一頻率的集中與放大現(xiàn)象。這也說明周圍土體對地鐵車站地震動力反應(yīng)有重要影響。

    圖10 西安人工波作用下地鐵車站不同高度處加速度反應(yīng)的傅里葉譜Fig.10 Fourier spectrum of accelerations at different height of station under Xi’an artificial wave

    輸入0.8 g峰值加速度時(shí),不同頻譜特性地震動作用下地鐵車站加速度反應(yīng)傅里葉譜如圖11所示。由圖可知,結(jié)構(gòu)內(nèi)加速度反應(yīng)傅里葉譜主頻帶寬具有顯著差異。由于松潘波主頻帶較寬,西安人工波主頻帶寬居中,Taft波主頻帶較窄,因此,結(jié)構(gòu)內(nèi)各測點(diǎn)加速度反應(yīng)傅里葉譜主頻帶寬也呈現(xiàn)出松潘波作用時(shí)較寬,西安人工波作用時(shí)居中,Taft波作用時(shí)較窄的現(xiàn)象。不同峰值加速度作用下,結(jié)構(gòu)加速度傅里葉譜隨輸入地震動頻譜特性的變化規(guī)律相似。同時(shí),在0.8 g松潘波作用下,地鐵車站頂部A22加速度反應(yīng)傅里葉譜呈現(xiàn)出雙峰現(xiàn)象,這可能與地鐵車站結(jié)構(gòu)在該工況的強(qiáng)震持時(shí)內(nèi)產(chǎn)生了較大的損傷破壞有關(guān)。

    圖11 不同頻譜特性地震動作用下地鐵車站加速度反應(yīng)的傅里葉譜Fig.11 Fourier spectrum of accelerations in station under different types of ground motions

    4 地鐵車站地震反應(yīng)的空間效應(yīng)

    4.1 地鐵車站加速度反應(yīng)的空間效應(yīng)

    地鐵車站結(jié)構(gòu)不同觀察斷面的峰值加速度如圖12所示。對比主觀測斷面和2號次觀測斷面的加速度反應(yīng)發(fā)現(xiàn),A3與A20測點(diǎn)的峰值加速度反應(yīng)相差甚微,當(dāng)輸入峰值加速度小于0.4 g時(shí),各工況下其峰值加速度差值不大于0.005 g,當(dāng)輸入峰值加速度大于0.4 g時(shí),其峰值加速度差值不大于0.05 g;而A4與A21的峰值加速度反應(yīng)表現(xiàn)出一定的差異,A4測點(diǎn)峰值加速度始終大于A20測點(diǎn)的峰值加速度,但各工況下其峰值加速度的差值不大于0.05 g。由于主觀測斷面和2號次觀測斷距離結(jié)構(gòu)端部距離相同,因此底部測點(diǎn)A3與A20的峰值加速度反應(yīng)較接近;同時(shí),受模型地基內(nèi)局部土體震密及不均勻的非線性塑性變形等影響,結(jié)構(gòu)在地震動作用下出現(xiàn)傾斜與扭轉(zhuǎn),且隨著高度增加, 結(jié)構(gòu)傾斜與扭轉(zhuǎn)作用對不同觀測斷面內(nèi)峰值加速度的影響更加顯著,從而使A4的峰值加速度大于A20的峰值加速度。對比2號和3號次觀測斷面加速度反應(yīng)發(fā)現(xiàn),在各種工況下A23測點(diǎn)峰值加速度始終大于A22測點(diǎn)峰值加速度;在相同峰值加速度作用下,輸入松潘波時(shí)A23與A22的峰值加速度的差值最大,輸入Taft波和西安人工波時(shí)A23與A22的峰值加速度的差值相對較小;在相同頻譜特性地震動作用下,隨著輸入峰值加速度的增大,A23與A22的峰值加速度的差值逐漸增大,當(dāng)輸入峰值加速度不大于0.4 g時(shí),其差值均小于0.1 g,當(dāng)輸入峰值加速度大于0.4 g時(shí),其差值相對較大。與中間段相比,結(jié)構(gòu)縱向端頭處截面剛度較小,同時(shí)受縱向端頭截?cái)嗵庍吔鐥l件影響,地鐵車站縱向端頭處的地震動力反應(yīng)較劇烈,稱為端頭效應(yīng)。由于A23位于3號次觀測斷面,距離地鐵車站端頭較近,受端頭效應(yīng)影響較強(qiáng),且隨著地震動峰值加速度的增大,3號次觀測斷面受端頭效應(yīng)影響更加顯著,因此,A23測點(diǎn)的峰值加速度始終大于A22測點(diǎn)的峰值加速度,且隨著輸入峰值加速度的增加其峰值加速度的差值逐漸增大。綜上所述,受模型結(jié)構(gòu)縱向端頭效應(yīng)及土結(jié)相互作用中結(jié)構(gòu)傾斜與扭轉(zhuǎn)變形影響,地鐵車站的加速度反應(yīng)具有顯著的空間效應(yīng)。

    圖12 地鐵車站峰值加速度反應(yīng)的空間分布Fig.12 Spatial distribution of PGA values of station under different ground motions

    4.2 地鐵車站應(yīng)變反應(yīng)的空間效應(yīng)

    輸入不同峰值加速度時(shí),Taft波作用下主觀測斷面及1號次觀測斷面的應(yīng)變幅值見表5~6。由表5可知,除少數(shù)工況外,隨著輸入峰值加速度增大,地鐵車站各測點(diǎn)應(yīng)變幅值逐漸增大;中柱應(yīng)變幅值較大,側(cè)墻的應(yīng)變幅值居中,頂板及底板的應(yīng)變幅值較小。由于地鐵車站與黃土介質(zhì)動力相互作用過程中結(jié)構(gòu)不同高度處動力反應(yīng)的差異及構(gòu)件自身質(zhì)量與剛度的不同,

    地鐵車站各構(gòu)件的應(yīng)變幅值呈現(xiàn)出不同的分布特征。輸入峰值加速度小于0.4 g時(shí),中柱應(yīng)變幅值表現(xiàn)為上層柱頂大于上層柱底,下層柱底大于下層柱頂,且上層柱頂大于下層柱底;側(cè)墻內(nèi)與頂板連接處應(yīng)變幅值較大,與底板連接處居中,與中板連接處較小。在較大地震動作用下,土與結(jié)構(gòu)動力相互作用劇烈,由于作用在結(jié)構(gòu)上動土壓力的不均勻性及結(jié)構(gòu)局部損傷破壞等影響,結(jié)構(gòu)中各構(gòu)件應(yīng)變幅值有所不同。

    對比表5及表6中不同觀測斷面上層中柱頂部S4與S32及下層中柱底部S1與S29應(yīng)變幅值發(fā)現(xiàn),輸入峰值加速度不大于0.6 g時(shí),由于1號次觀測斷面距端部較近,其截面剛度小于遠(yuǎn)離結(jié)構(gòu)端部的主觀測斷面,整體變形較大,因此其應(yīng)變幅值均大于主觀測斷面;隨著輸入峰值加速度的增加,由于1號次觀測斷面應(yīng)變較大,在較大地震動作用下,1號次觀測斷面上層柱頂及下層柱底先出現(xiàn)受拉破壞,結(jié)構(gòu)內(nèi)力重分布,從而出現(xiàn)1號次觀測斷面拉應(yīng)變幅值減小。對比分析不同觀測斷面內(nèi)側(cè)墻頂部S12與S36及側(cè)墻底部S9與S33的應(yīng)變反應(yīng)表明,在較小峰值加速度作用下,側(cè)墻應(yīng)變幅值較小,受應(yīng)變片測量精度的影響,應(yīng)變反應(yīng)規(guī)律不明顯。在較大峰值加速度作用下,側(cè)墻頂部主觀測斷面的拉應(yīng)變幅值大于1號次觀測斷面,壓應(yīng)變幅值小于1號次觀測斷面,而側(cè)墻底部呈現(xiàn)出相反的規(guī)律,這說明在土與結(jié)構(gòu)動力相互作用過程中1號次觀測斷面相對主觀測斷面發(fā)生了逆時(shí)針扭轉(zhuǎn),從而使拉應(yīng)變與壓應(yīng)變呈現(xiàn)出不同的特征。綜上所述,受結(jié)構(gòu)縱向端部效應(yīng)及土結(jié)相互作用中結(jié)構(gòu)傾斜與扭轉(zhuǎn)變形影響,使得結(jié)構(gòu)應(yīng)變反應(yīng)具有顯著的空間效應(yīng)。

    表5 Taft波作用下主觀測斷面結(jié)構(gòu)應(yīng)變幅值(με)

    表6 Taft波作用下1號次觀測斷面結(jié)構(gòu)應(yīng)變幅值(με)

    5 地鐵車站地震破壞

    試驗(yàn)結(jié)束后模型結(jié)構(gòu)發(fā)生了嚴(yán)重破壞。主觀測斷面中柱出現(xiàn)典型的剪壓破壞,產(chǎn)生豎向裂縫,1號次觀測斷面中柱產(chǎn)生斜裂縫,混凝土大量剝落,縱向受力鋼筋外露,如圖13(a)所示;中柱與頂、底板連接處破壞最為嚴(yán)重,如圖13(b)所示;1號次觀測斷面處側(cè)墻與底板、頂板連接處產(chǎn)生裂縫,局部腋角鋼筋從混凝土基材中拉出如圖13(c)所示。由于框架式地鐵車站結(jié)構(gòu)為中柱承重,水平地震動作用下,中柱承受豎向壓力的同時(shí)承受較大的水平向剪力;且與側(cè)墻、板構(gòu)件相比,中柱截面小,整體剛度小,因此,地震時(shí),中柱在雙向動力作用下發(fā)生嚴(yán)重的剪壓破壞。同時(shí),由于中柱內(nèi)與頂、底板連接處附加彎矩較大,因此破壞最為嚴(yán)重;1號次觀測面破壞較主觀測面嚴(yán)重,這與1號次觀測面距離模型結(jié)構(gòu)縱向端頭較近有關(guān)。

    綜上所述,地鐵車站中柱頂、底部及側(cè)墻與板構(gòu)件連接處地震破壞嚴(yán)重;距結(jié)構(gòu)縱向端頭較近的1號次觀測斷面破壞較主觀測斷面更為嚴(yán)重,模型結(jié)構(gòu)的地震破壞呈現(xiàn)出顯著的空間效應(yīng)。

    圖13 地鐵車站的地震破壞Fig.13 Seismic damage of subway station

    6 結(jié) 論

    本文進(jìn)行了黃土場地與地鐵車站結(jié)構(gòu)動力相互作用大型地震模擬振動臺試驗(yàn),根據(jù)實(shí)測數(shù)據(jù)對模型地基及地鐵車站的加速度反應(yīng)進(jìn)行了頻域分析;同時(shí),對比不同觀測斷面內(nèi)地鐵車站的地震反應(yīng),分析了地鐵車站地震反應(yīng)的空間效應(yīng)。研究結(jié)論可為黃土地區(qū)地鐵地下結(jié)構(gòu)的抗震設(shè)計(jì)及相關(guān)理論研究提供重要參考,主要結(jié)論如下:

    (1) 由于地基對地震動有低頻放大和高頻過濾效應(yīng),且隨著輸入峰值加速度增大而更加顯著,因此,地震動較大時(shí),隨土層深度減小地基加速度反應(yīng)傅里葉譜的低頻成分逐漸增大,高頻成分逐漸減?。煌瑫r(shí),隨著輸入峰值加速度增加,地基加速度反應(yīng)傅里葉譜中增強(qiáng)的低頻成分頻帶寬度逐漸縮小。

    (2) 隨著輸入峰值加速度增加,由于結(jié)構(gòu)損傷累計(jì),固有頻率降低,地鐵車站加速度反應(yīng)傅里葉譜主頻集中,帶寬變窄;同時(shí),隨著輸入峰值加速度增加,土與結(jié)構(gòu)相互作用逐漸劇烈,地鐵車站與周圍土體之間出現(xiàn)偶發(fā)性的脫離與碰撞現(xiàn)象,土體撞擊結(jié)構(gòu)產(chǎn)生高頻振動,從而結(jié)構(gòu)加速度反應(yīng)傅里葉譜的高頻成分呈現(xiàn)出一定的增強(qiáng)效應(yīng)。

    (3) 隨著輸入峰值加速度增加,地基與結(jié)構(gòu)中加速度反應(yīng)傅里葉譜的基頻逐漸減小,相應(yīng)傅里葉譜值逐漸增大;由于松潘波的主頻帶較寬,西安人工波主頻帶寬居中,Taft波主頻帶較窄,因此,地基與結(jié)構(gòu)中加速度反應(yīng)傅里葉譜主頻帶寬也呈現(xiàn)出松潘波作用時(shí)較寬,西安人工波作用時(shí)居中,Taft波作用時(shí)較窄。

    (4) 地基內(nèi)各測點(diǎn)的峰值加速度均隨土層埋深的減小而增大;地震動較大時(shí),結(jié)構(gòu)的峰值加速度均呈現(xiàn)出頂部和底部較大,中部相對較小,且頂部峰值加速度略小于底部。因此,在黃土地區(qū)地鐵地下結(jié)構(gòu)的抗震設(shè)計(jì)中計(jì)算模型慣性力的施加應(yīng)考慮構(gòu)件加速度反應(yīng)的差異。

    (5) 結(jié)構(gòu)應(yīng)變反應(yīng)及震害觀測均表明中柱頂、底部及側(cè)墻與板構(gòu)件連接處為地鐵車站抗震薄弱部位,因此,在黃土地區(qū)地鐵地下結(jié)構(gòu)抗震設(shè)計(jì)中應(yīng)加強(qiáng)中柱頂、底部及側(cè)墻與板構(gòu)件連接處的延性,以提高其變形性能。

    (6) 受結(jié)構(gòu)縱向端頭效應(yīng)及土結(jié)相互作用中結(jié)構(gòu)傾斜與扭轉(zhuǎn)變形影響,地鐵車站結(jié)構(gòu)地震反應(yīng)具有顯著的空間效應(yīng)。在振動臺模型試驗(yàn)中,可調(diào)整模型結(jié)構(gòu)縱向端頭處剛度及控制模型地基密實(shí)均勻性,以減小模型結(jié)構(gòu)地震反應(yīng)的空間效應(yīng),提高黃土地區(qū)地鐵地下結(jié)構(gòu)振動臺試驗(yàn)的可靠性。

    [1] 權(quán)登州, 王毅紅, 井彥林,等.黃土地區(qū)地鐵車站數(shù)值模型及測試位置研究[J].震災(zāi)防御技術(shù), 2015, 10(1): 108-115. QUAN Dengzhou, WANG Yihong, JING Yanlin, et al. Numerical model and measuring position for subway stations in loess area [J]. Technology for Earthquake Disaster Prevention, 2015, 10(1): 108-115.

    [2] 王志杰, 駱亞生, 楊永俊. 不同地區(qū)非飽和黃土動力結(jié)構(gòu)性研究[J]. 巖土力學(xué), 2010, 31(8): 2459-2464. WANG Zhijie,LUO Yasheng,YANG Yongjun. Study of dynamic structural characteristics of unsaturated loesses in different regions [J]. Rock and Soil Mechanics, 2010, 31(8): 2459-2464.

    [3] JAFARZDEH F, FARAHI JAHROMI H, ABZAARI TORGHABEH E. Investigating dynamic response of a buried pipeline in sandy soil layer by 1G shaking table test [J]. International Journal of Civil Engineering, 2010,8(2): 107-124.

    [4] AMOROSI A, BOLDINI D. Numerical modelling of the transverse dynamic behavior of circular tunnels in clayey soils [J]. Soil Dynamics and Earthquake Engineering, 2009,29(6): 1059-1072.

    [5] AZADI M, MIR MOHAMMAD HOSSEINI S M. Analyses of the effect of seismic behavior of shallow tunnels in liquefiable grounds [J]. Tunnelling and Underground Space Technology, 2010, 25:543-552.

    [6] CHIAN S C. Effect of buried depth and diameter on uplift of underground structures in liquefied soils[J] Soil Dynamics and Earthquake Engineering, 2012(41):181-190.

    [7] 楊林德, 季倩倩, 鄭永來, 等. 軟土地鐵車站結(jié)構(gòu)的振動臺試驗(yàn)[J]. 現(xiàn)代隧道技術(shù), 2003, 40(1): 7-11. YANG Linde, JI Qianqian, ZHENG Yonglai, et al. Shaking table test on metro station structures in soft soil[J]. Modern Tunnelling Technolog, 2003, 40(1): 7-11.

    [8] CHEN Guoxing, WANG Zhihua, ZUO Xi, et al. Shaking table test on the seismic failure characteristics of a subway station structure on liquefiable ground [J]. Earthquake Engineering & Structural Dynamics,2013; 42:1489-1507.

    [9] CHEN Guoxing, CHEN Su, ZUO Xi, et al. Shaking-table tests and numerical simulations on a subway structure in soft soil[J]. Soil dynamic and earthquake engineering, 2015,76:13-28.

    [10] 陳國興, 左熹, 王志華, 等. 近遠(yuǎn)場地震作用下液化地基上地鐵車站結(jié)構(gòu)動力損傷特性的振動臺試驗(yàn) [J]. 土木工程學(xué)報(bào),2010,43(12):120-126. CHEN Guoxing, ZUO Xi, WANG Zhihua, et al. Large scale shaking table test study of the dynamic damage behavior of subway station structures in liquefiable foundation under near-fault and far-field ground motions[J]. China Civil Engineering Journal, 2010,43(12):120-126.

    [11] 劉晶波,劉祥慶,王宗綱,等土-結(jié)構(gòu)動力相互作用系統(tǒng)離心機(jī)振動臺模型試驗(yàn)[J]. 土木工程學(xué)報(bào). 2010,43(11): 114-121. LIU Jingbo, LIU Xiangqing, WANG Zonggang, et al. Dynamic centrifuge model test of a soil-structure interaction system[J]. China Civil Engineering Journal, 2010,43(11): 114-121.

    [12] 劉妮娜,門玉明. 地裂縫環(huán)境下的地鐵隧道-地層地震動力相互作用研究[D]. 西安:長安大學(xué),2010: 4-9.

    [13] VARGHESE R M, LATHA G M. Shaking table tests to investigate the influence of various factors on the liquefaction resistance of sands [J]. Nat Hazards, 2014, 73:1337-1351.

    [14] LANZANO G, BILOTTA E, RUSSO G, et al. Centrifuge modeling of seismic loading on tunnels in sand [J]. Geotechnical Testing Journal, 2012,35:1-16.

    [15] 陳國興,王志華,左熹,等. 振動臺試驗(yàn)疊層剪切型土箱的研制[J]. 巖土工程學(xué)報(bào),2010,32 (1): 89-97. CHEN Guoxing, WANG Zhihua, ZUO Xi, et al. Development of laminar shear soil container for shaking table tests [J]. Chinese Journal of Geotechnical Engineering, 2010,32 (1): 89-97.

    Spatial effects and frequency domain analysis for seismic responses of subway station in loess area

    QUAN Dengzhou1, WANG Yihong1, MA Pengbo1, JING Yanlin1, CHEN Su2

    (1. School of Civil Engineering, Chang’an University, Xi’an 710061, China;2. Institute of Geophysics, China Seismological Bureau, Beijing 100081, China)

    The shaking table tests of underground subway station in loess area were performed, the acceleration and strain responses of soil and the structure were measured. The frequency domain analysis for acceleration responses of the model system was conducted, and the spatial effects on seismic responses of the station were analyzed by comparing data of different observation sections. The results showed that the fundamental frequencies of the model system’s acceleration Fourier spectra decrease with increase in PGA, meanwhile, the corresponding amplitudes of Fourier spectra increase with increase in PGA; the Fourier spectra values of low-frequency components of accelerations increase and their high-frequency components decrease, and the bandwidth of low-frequency components decreases with increase in PGA; the Fourier spectra values of high-frequency components of accelerations of the station increase and the main frequencies of the spectra are concentrated and become narrower with increase in PGA; the main frequency band width of Fourier spectra of accelerations of the model system is wider under the action of Songpan wave, it is medium under the action of Xi’an artificial wave and narrower under the action of Taft wave; meanwhile, the spatial effects on seismic responses of the station are obvious because of end effects of the structure and torsion in interaction between soil and structure. The results provided a reference for aseismic design and theoretical study of underground structures in loess area.

    loess area; subway station; shaking table test; acceleration response; frequency domain analysis; spatial effect

    國家自然科學(xué)基金(41472267);中央級公益性科研專項(xiàng)基金項(xiàng)目(DQJB14B50);西安市地下鐵道有限責(zé)任公司科研基金(D4-YJ-042014048)

    2015-09-09 修改稿收到日期:2015-10-13

    權(quán)登州 男,博士,講師,1983年8月生

    P315.97

    A

    10.13465/j.cnki.jvs.2016.21.016

    猜你喜歡
    傅里葉震動峰值
    “四單”聯(lián)動打造適齡兒童隊(duì)前教育峰值體驗(yàn)
    震動減脂儀可以減肥?
    雙線性傅里葉乘子算子的量化加權(quán)估計(jì)
    基于小波降噪的稀疏傅里葉變換時(shí)延估計(jì)
    水電工程場地地震動確定方法
    振動攪拌 震動創(chuàng)新
    中國公路(2017年18期)2018-01-23 03:00:38
    寬占空比峰值電流型準(zhǔn)PWM/PFM混合控制
    基于峰值反饋的電流型PFM控制方法
    基于傅里葉變換的快速TAMVDR算法
    快速離散傅里葉變換算法研究與FPGA實(shí)現(xiàn)
    電測與儀表(2015年5期)2015-04-09 11:30:44
    深夜a级毛片| 美女 人体艺术 gogo| 欧美zozozo另类| 淫秽高清视频在线观看| 国产精品,欧美在线| 男女下面进入的视频免费午夜| 在线观看一区二区三区| 国产乱人视频| 亚洲精品一卡2卡三卡4卡5卡| 赤兔流量卡办理| 久久久精品大字幕| 又爽又黄a免费视频| 精品久久久久久久久久久久久| 亚洲第一电影网av| 色尼玛亚洲综合影院| 亚洲欧美精品综合久久99| 国产v大片淫在线免费观看| 成年女人毛片免费观看观看9| 国产成人av教育| 在线免费观看不下载黄p国产 | 久久国产精品人妻蜜桃| 亚洲人成网站在线播| 精品99又大又爽又粗少妇毛片 | 欧美激情国产日韩精品一区| 国产精品人妻久久久久久| 中文字幕高清在线视频| 尾随美女入室| 91狼人影院| 熟妇人妻久久中文字幕3abv| 不卡视频在线观看欧美| 黄色一级大片看看| 在线a可以看的网站| 女同久久另类99精品国产91| 精品福利观看| 国内精品一区二区在线观看| 日本熟妇午夜| 黄色配什么色好看| 黄色配什么色好看| 国产白丝娇喘喷水9色精品| 一个人看的www免费观看视频| 亚洲欧美日韩东京热| 99久久九九国产精品国产免费| 成年版毛片免费区| .国产精品久久| 中文亚洲av片在线观看爽| 一区二区三区高清视频在线| 91午夜精品亚洲一区二区三区 | 春色校园在线视频观看| x7x7x7水蜜桃| 亚洲美女搞黄在线观看 | 国产在视频线在精品| 高清在线国产一区| 淫妇啪啪啪对白视频| 日韩欧美在线乱码| 国产精品99久久久久久久久| 国内精品美女久久久久久| 国产高清视频在线播放一区| 婷婷精品国产亚洲av在线| 国产午夜精品久久久久久一区二区三区 | 日韩精品青青久久久久久| 99国产极品粉嫩在线观看| 亚洲美女视频黄频| 国产精品免费一区二区三区在线| 亚洲熟妇熟女久久| 露出奶头的视频| a级毛片a级免费在线| 听说在线观看完整版免费高清| 日本免费a在线| 久久久久国内视频| 精品欧美国产一区二区三| x7x7x7水蜜桃| 女的被弄到高潮叫床怎么办 | 国产免费一级a男人的天堂| av福利片在线观看| 在线观看66精品国产| 人妻制服诱惑在线中文字幕| 特级一级黄色大片| 婷婷亚洲欧美| 国语自产精品视频在线第100页| 国产精品亚洲一级av第二区| 国产一级毛片七仙女欲春2| 91午夜精品亚洲一区二区三区 | 成年女人毛片免费观看观看9| 乱码一卡2卡4卡精品| 欧美最黄视频在线播放免费| 天堂av国产一区二区熟女人妻| 少妇的逼水好多| 国产免费男女视频| 午夜福利高清视频| 成年女人看的毛片在线观看| 日韩欧美 国产精品| 亚洲欧美日韩无卡精品| 欧美日韩精品成人综合77777| 少妇丰满av| 国产精品免费一区二区三区在线| 91麻豆精品激情在线观看国产| 变态另类丝袜制服| 欧美性猛交黑人性爽| 黄色视频,在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久视频播放| 18禁黄网站禁片午夜丰满| 人妻久久中文字幕网| 日韩欧美国产在线观看| 国内精品久久久久精免费| 免费一级毛片在线播放高清视频| 中文字幕av在线有码专区| 麻豆成人av在线观看| 日韩欧美国产一区二区入口| 18禁黄网站禁片免费观看直播| 五月玫瑰六月丁香| 久久久久久国产a免费观看| 亚洲精品粉嫩美女一区| 亚洲国产欧洲综合997久久,| 国产欧美日韩精品亚洲av| 国产亚洲91精品色在线| 88av欧美| 丰满的人妻完整版| 亚洲av日韩精品久久久久久密| 精品久久国产蜜桃| 亚洲黑人精品在线| www.www免费av| 国产成人影院久久av| 成年免费大片在线观看| 九九爱精品视频在线观看| 久久人妻av系列| 啦啦啦韩国在线观看视频| 日本熟妇午夜| 999久久久精品免费观看国产| 日本色播在线视频| 国产精品,欧美在线| 国产一区二区激情短视频| 日本a在线网址| 国产欧美日韩一区二区精品| 最新在线观看一区二区三区| 中文字幕免费在线视频6| 啦啦啦观看免费观看视频高清| 国产精品99久久久久久久久| 久久久久久久久久成人| 又紧又爽又黄一区二区| 日韩一本色道免费dvd| 精品久久久久久,| 国产精品一区二区三区四区久久| 免费大片18禁| 毛片女人毛片| 婷婷丁香在线五月| 一卡2卡三卡四卡精品乱码亚洲| 成人欧美大片| 欧美日韩瑟瑟在线播放| 欧美日韩综合久久久久久 | 亚洲精品亚洲一区二区| 欧美日韩精品成人综合77777| 天堂动漫精品| 偷拍熟女少妇极品色| 日韩高清综合在线| 欧美色欧美亚洲另类二区| 国产在线精品亚洲第一网站| 国产精品人妻久久久久久| 夜夜看夜夜爽夜夜摸| 97人妻精品一区二区三区麻豆| 九九久久精品国产亚洲av麻豆| 成年女人永久免费观看视频| 偷拍熟女少妇极品色| 久久久国产成人免费| 成年免费大片在线观看| 大又大粗又爽又黄少妇毛片口| 国产精品1区2区在线观看.| 亚洲精品一区av在线观看| 亚洲av免费在线观看| 国产精品精品国产色婷婷| 春色校园在线视频观看| 国内精品久久久久久久电影| 久久久久久久久久久丰满 | 网址你懂的国产日韩在线| 欧美一区二区精品小视频在线| 日本 av在线| 亚洲狠狠婷婷综合久久图片| 69av精品久久久久久| 亚洲av免费高清在线观看| 精品午夜福利视频在线观看一区| 欧美高清成人免费视频www| 亚洲中文字幕一区二区三区有码在线看| 久久精品人妻少妇| 99热只有精品国产| 丝袜美腿在线中文| 男人狂女人下面高潮的视频| 啦啦啦观看免费观看视频高清| 在线a可以看的网站| 国产极品精品免费视频能看的| 精华霜和精华液先用哪个| 国产精品人妻久久久影院| 亚洲aⅴ乱码一区二区在线播放| 亚洲av免费高清在线观看| 亚洲美女搞黄在线观看 | 在线播放国产精品三级| 一a级毛片在线观看| 波多野结衣高清作品| 免费搜索国产男女视频| 麻豆精品久久久久久蜜桃| 亚洲天堂国产精品一区在线| 特级一级黄色大片| 精品人妻偷拍中文字幕| 一进一出好大好爽视频| 偷拍熟女少妇极品色| 成人av在线播放网站| av黄色大香蕉| 国产伦精品一区二区三区四那| 欧美日韩综合久久久久久 | 免费观看人在逋| 亚洲成人久久爱视频| 亚洲在线自拍视频| 真人做人爱边吃奶动态| 国产激情偷乱视频一区二区| 女同久久另类99精品国产91| 村上凉子中文字幕在线| 九九热线精品视视频播放| 午夜爱爱视频在线播放| 99久国产av精品| 精品一区二区三区av网在线观看| 婷婷亚洲欧美| 欧美bdsm另类| av在线天堂中文字幕| 1024手机看黄色片| 亚洲人成网站高清观看| 国产视频一区二区在线看| 香蕉av资源在线| 午夜免费激情av| 日韩欧美在线二视频| 老熟妇乱子伦视频在线观看| 国产不卡一卡二| 日日干狠狠操夜夜爽| 亚洲av不卡在线观看| 国产毛片a区久久久久| 少妇丰满av| 女人十人毛片免费观看3o分钟| 啦啦啦啦在线视频资源| 午夜免费男女啪啪视频观看 | 一进一出抽搐gif免费好疼| 欧美精品国产亚洲| 国产精品美女特级片免费视频播放器| 国内精品久久久久精免费| 成年女人看的毛片在线观看| 长腿黑丝高跟| 精品一区二区三区视频在线| av专区在线播放| 1000部很黄的大片| 亚洲自偷自拍三级| 日韩亚洲欧美综合| 99精品在免费线老司机午夜| 亚洲一区二区三区色噜噜| 亚洲精品色激情综合| 国模一区二区三区四区视频| 日本一本二区三区精品| 12—13女人毛片做爰片一| 国产精品三级大全| 桃色一区二区三区在线观看| 天堂av国产一区二区熟女人妻| 久久久久久九九精品二区国产| 国产精品免费一区二区三区在线| 男人和女人高潮做爰伦理| 国产精品电影一区二区三区| 欧美又色又爽又黄视频| 亚洲av日韩精品久久久久久密| 国产中年淑女户外野战色| 欧美日韩国产亚洲二区| 一级a爱片免费观看的视频| 听说在线观看完整版免费高清| 欧美3d第一页| 久久精品国产亚洲av涩爱 | 国产大屁股一区二区在线视频| 精品不卡国产一区二区三区| 国产精品人妻久久久影院| 乱系列少妇在线播放| 国产欧美日韩精品一区二区| 午夜免费成人在线视频| а√天堂www在线а√下载| 国产精品久久电影中文字幕| 波野结衣二区三区在线| 久久久久久久久大av| 美女xxoo啪啪120秒动态图| 亚洲欧美精品综合久久99| 婷婷精品国产亚洲av| 高清日韩中文字幕在线| 欧美色视频一区免费| 在线国产一区二区在线| 欧美又色又爽又黄视频| 一个人免费在线观看电影| 午夜激情福利司机影院| 天堂av国产一区二区熟女人妻| 国产av在哪里看| 欧美一区二区亚洲| 特大巨黑吊av在线直播| 国产精品98久久久久久宅男小说| 亚州av有码| 人妻少妇偷人精品九色| 亚洲男人的天堂狠狠| 三级男女做爰猛烈吃奶摸视频| 十八禁网站免费在线| 桃色一区二区三区在线观看| 亚洲最大成人中文| 亚洲熟妇中文字幕五十中出| 午夜激情福利司机影院| 免费看光身美女| 丰满的人妻完整版| 日本欧美国产在线视频| 国内精品久久久久精免费| 老熟妇仑乱视频hdxx| 极品教师在线免费播放| 精品午夜福利视频在线观看一区| 亚洲av二区三区四区| 色吧在线观看| 可以在线观看的亚洲视频| 亚洲最大成人中文| 99热网站在线观看| 俺也久久电影网| 国产色婷婷99| 国产精品福利在线免费观看| 91久久精品国产一区二区成人| 精品人妻一区二区三区麻豆 | 麻豆一二三区av精品| 日本色播在线视频| 99久久无色码亚洲精品果冻| 亚洲四区av| 亚洲av成人精品一区久久| 免费在线观看影片大全网站| 91久久精品国产一区二区三区| 成人二区视频| 一个人看视频在线观看www免费| 欧美bdsm另类| 成年人黄色毛片网站| 韩国av一区二区三区四区| 午夜福利欧美成人| 久久久久九九精品影院| 尤物成人国产欧美一区二区三区| 日韩欧美 国产精品| 我要看日韩黄色一级片| 亚洲精品乱码久久久v下载方式| 小说图片视频综合网站| 欧美日韩中文字幕国产精品一区二区三区| 国产日本99.免费观看| 黄色配什么色好看| 精品乱码久久久久久99久播| 久99久视频精品免费| 联通29元200g的流量卡| 无遮挡黄片免费观看| 欧美激情在线99| 亚洲欧美日韩高清在线视频| 亚洲精品亚洲一区二区| 成人午夜高清在线视频| 精品久久久久久久久久久久久| 国产免费一级a男人的天堂| 最近视频中文字幕2019在线8| 亚洲精品456在线播放app | 中国美女看黄片| 人妻丰满熟妇av一区二区三区| 国内少妇人妻偷人精品xxx网站| 久久久久久大精品| 免费黄网站久久成人精品| 12—13女人毛片做爰片一| 亚洲狠狠婷婷综合久久图片| 亚洲欧美激情综合另类| 亚洲av五月六月丁香网| 免费av毛片视频| 成人性生交大片免费视频hd| 久久久久久久久久成人| 丰满人妻一区二区三区视频av| 三级毛片av免费| 99久久九九国产精品国产免费| 丰满人妻一区二区三区视频av| 99热6这里只有精品| 免费av不卡在线播放| 有码 亚洲区| 啦啦啦啦在线视频资源| 一级黄色大片毛片| 久久久久久久午夜电影| 99热这里只有精品一区| 久久热精品热| 日韩av在线大香蕉| 国内精品久久久久精免费| 搡老岳熟女国产| 亚洲四区av| 不卡视频在线观看欧美| 五月玫瑰六月丁香| 亚洲欧美精品综合久久99| 精品一区二区三区人妻视频| av视频在线观看入口| 国产探花在线观看一区二区| 精品人妻视频免费看| 欧美色视频一区免费| 亚洲欧美日韩东京热| 波多野结衣高清作品| 人妻久久中文字幕网| 中文字幕av成人在线电影| 性插视频无遮挡在线免费观看| 小蜜桃在线观看免费完整版高清| 久久精品夜夜夜夜夜久久蜜豆| 午夜爱爱视频在线播放| 看免费成人av毛片| 国产男人的电影天堂91| 极品教师在线视频| 日本-黄色视频高清免费观看| 免费看美女性在线毛片视频| 又爽又黄a免费视频| 国产午夜精品久久久久久一区二区三区 | 99热这里只有是精品50| 2021天堂中文幕一二区在线观| 高清毛片免费观看视频网站| 中文字幕高清在线视频| 久久久成人免费电影| 麻豆精品久久久久久蜜桃| h日本视频在线播放| 变态另类丝袜制服| 人妻久久中文字幕网| av.在线天堂| 两个人视频免费观看高清| 小蜜桃在线观看免费完整版高清| 男女做爰动态图高潮gif福利片| 久久久久性生活片| av福利片在线观看| 久久精品影院6| 91精品国产九色| 国产一级毛片七仙女欲春2| 我要搜黄色片| 国产在线精品亚洲第一网站| 精品久久久久久,| 国产精品嫩草影院av在线观看 | 少妇丰满av| 给我免费播放毛片高清在线观看| 天堂av国产一区二区熟女人妻| 亚洲精品成人久久久久久| 欧美色视频一区免费| 不卡视频在线观看欧美| 午夜a级毛片| 桃色一区二区三区在线观看| 日本黄色视频三级网站网址| 制服丝袜大香蕉在线| 国产精品电影一区二区三区| 国产主播在线观看一区二区| 麻豆国产av国片精品| 99久久九九国产精品国产免费| 日本色播在线视频| av中文乱码字幕在线| 亚洲国产精品成人综合色| 大又大粗又爽又黄少妇毛片口| 国产高清激情床上av| 啦啦啦啦在线视频资源| 亚洲中文字幕日韩| 亚洲成a人片在线一区二区| 欧美区成人在线视频| 久久久久久久精品吃奶| 真人做人爱边吃奶动态| 97超视频在线观看视频| 成人亚洲精品av一区二区| 亚洲真实伦在线观看| 成人av在线播放网站| 亚洲色图av天堂| 99九九线精品视频在线观看视频| 日韩大尺度精品在线看网址| 日韩精品有码人妻一区| 欧美色欧美亚洲另类二区| 国产精品一区二区性色av| 色精品久久人妻99蜜桃| 不卡视频在线观看欧美| 国产精品精品国产色婷婷| 精品久久久噜噜| 久久精品国产亚洲av香蕉五月| 欧美色欧美亚洲另类二区| 午夜福利高清视频| 亚洲精品久久国产高清桃花| 欧美日韩瑟瑟在线播放| 搡老妇女老女人老熟妇| 99久久成人亚洲精品观看| 波多野结衣高清无吗| 久久久色成人| 国产视频一区二区在线看| 欧美成人性av电影在线观看| 亚洲av中文字字幕乱码综合| 美女高潮喷水抽搐中文字幕| 国产久久久一区二区三区| 91av网一区二区| 免费在线观看成人毛片| 国产熟女欧美一区二区| 色播亚洲综合网| 色综合婷婷激情| 亚洲国产欧美人成| 亚洲av成人精品一区久久| 日韩欧美一区二区三区在线观看| 99久国产av精品| 日本与韩国留学比较| 亚洲成人久久爱视频| 国产精品久久久久久久电影| 天天躁日日操中文字幕| 国产伦在线观看视频一区| 变态另类丝袜制服| av在线老鸭窝| 色综合亚洲欧美另类图片| 亚洲人成网站在线播| 精华霜和精华液先用哪个| 欧美色欧美亚洲另类二区| 老熟妇仑乱视频hdxx| 99在线视频只有这里精品首页| 欧美日韩综合久久久久久 | 三级毛片av免费| 亚洲成人久久性| 国产精品嫩草影院av在线观看 | 五月玫瑰六月丁香| 黄色丝袜av网址大全| 精品国产三级普通话版| 亚洲av电影不卡..在线观看| 日韩欧美精品v在线| 热99re8久久精品国产| 国内少妇人妻偷人精品xxx网站| 成年人黄色毛片网站| 亚洲国产高清在线一区二区三| 日韩强制内射视频| 日日夜夜操网爽| 有码 亚洲区| 亚洲av中文字字幕乱码综合| 有码 亚洲区| 国产精品野战在线观看| 亚洲精品粉嫩美女一区| 日韩大尺度精品在线看网址| 国产亚洲91精品色在线| 国产高清三级在线| 亚洲自偷自拍三级| 亚洲欧美日韩高清在线视频| 国产精品一及| 国产久久久一区二区三区| 久久亚洲精品不卡| 亚洲avbb在线观看| 少妇丰满av| 啦啦啦啦在线视频资源| 五月伊人婷婷丁香| 日本黄色片子视频| 日日啪夜夜撸| 色视频www国产| 国产精品久久久久久亚洲av鲁大| 精品人妻熟女av久视频| 日韩人妻高清精品专区| 非洲黑人性xxxx精品又粗又长| 91午夜精品亚洲一区二区三区 | 日韩一本色道免费dvd| 国产麻豆成人av免费视频| 国产极品精品免费视频能看的| 禁无遮挡网站| 欧美日韩综合久久久久久 | 精品99又大又爽又粗少妇毛片 | 最近最新中文字幕大全电影3| 国产精品人妻久久久影院| 成人美女网站在线观看视频| 国产亚洲91精品色在线| 窝窝影院91人妻| 国产精品免费一区二区三区在线| 99热这里只有是精品在线观看| 嫩草影视91久久| 亚洲成a人片在线一区二区| 国产毛片a区久久久久| 国产精品爽爽va在线观看网站| 小蜜桃在线观看免费完整版高清| 听说在线观看完整版免费高清| 国产探花极品一区二区| 成人午夜高清在线视频| 日本成人三级电影网站| 免费观看的影片在线观看| 亚洲电影在线观看av| 亚洲四区av| 我的老师免费观看完整版| 亚洲欧美日韩高清专用| 亚洲成人久久性| 热99re8久久精品国产| 深夜a级毛片| 国产精品久久久久久久久免| 国产精品福利在线免费观看| 少妇人妻精品综合一区二区 | 国产一区二区三区av在线 | 少妇丰满av| 国产久久久一区二区三区| 国产蜜桃级精品一区二区三区| 亚洲图色成人| 久久中文看片网| 免费人成在线观看视频色| 久久精品综合一区二区三区| 日韩一区二区视频免费看| 亚洲精品粉嫩美女一区| 午夜日韩欧美国产| 国产精品自产拍在线观看55亚洲| 真人做人爱边吃奶动态| 露出奶头的视频| 久久精品国产亚洲av涩爱 | 精品久久久久久成人av| 丰满的人妻完整版| 亚洲欧美日韩高清在线视频| 国国产精品蜜臀av免费| 亚洲午夜理论影院| 超碰av人人做人人爽久久| 国产精品乱码一区二三区的特点| 一个人观看的视频www高清免费观看| av天堂中文字幕网| 好男人在线观看高清免费视频| 亚洲av五月六月丁香网| 欧美在线一区亚洲| 变态另类成人亚洲欧美熟女| 久久久久久久精品吃奶| 小说图片视频综合网站| 亚洲av免费高清在线观看| 精品不卡国产一区二区三区| 亚洲最大成人手机在线| 日本免费一区二区三区高清不卡| 在线看三级毛片| 一区福利在线观看| 天堂影院成人在线观看| 亚洲自偷自拍三级| 日韩大尺度精品在线看网址| 色哟哟·www| 国产免费男女视频|