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

    橋頭路堤-橋臺-地基作用體系的地震反應(yīng)分析

    2020-10-15 06:44:42程志明魏紅衛(wèi)
    關(guān)鍵詞:橋頭橋臺路堤

    程志明,魏紅衛(wèi)

    橋頭路堤-橋臺-地基作用體系的地震反應(yīng)分析

    程志明1, 2,魏紅衛(wèi)1

    (1. 中南大學(xué) 土木工程學(xué)院,湖南 長沙 410075;2. 武漢市政工程設(shè)計研究院有限責(zé)任公司,湖北 武漢 430023)

    考慮橋頭路堤、橋臺、地基相互作用,建立路橋過渡段三維動力分析模型,輸入汶川地震波,進行非線性地震反應(yīng)時程分析。將數(shù)值計算結(jié)果,與經(jīng)典理論和規(guī)范方法計算結(jié)果、離心試驗結(jié)果,以及震害調(diào)查資料對比分析,證實計算模型和計算方法的合理性。著重計算橋頭路堤、橋臺、地基的地震位移反應(yīng),橋臺的地震加速度、土壓力、彎矩以及基底剪應(yīng)力反應(yīng)等。計算結(jié)果表明:震后橋臺整體背離路堤滑移,并伴隨下沉和外傾,地震位移包含平移與轉(zhuǎn)動2部分,橋臺頂部加速度存在放大效應(yīng),承受的地震土壓力以及彎矩內(nèi)力較震前明顯增大,相比而言,公路規(guī)范與鐵路規(guī)范方法的計算值偏??;震后橋頭路堤向橋臺方向推移,存在拉裂趨勢和不均勻沉降破壞;地基與橋臺基礎(chǔ)相鄰部位,峰值水平位移與殘余水平位移以及最終沉降量相對較大,基底剪力大幅增加,橋臺與地基出現(xiàn)相對滑移。研究結(jié)果對橋梁抗震加固設(shè)計具有借鑒意義。

    橋頭路堤;橋臺;地基;地震反應(yīng);時程分析

    我國地處世界兩大地震帶之間,是一個多地震國家,四千二百余年的地震文獻記載表明,除浙江、江西兩省外,我國絕大部分地區(qū)都發(fā)生過震級較大的破壞性地震[1]。近年來的公路橋梁地震災(zāi)害調(diào)查表明,公路路基與橋梁銜接部位的地震破壞非常多見,是地震中的薄弱環(huán)節(jié)[2?5]。由于地震引起的橋頭路堤坍滑、沉陷,側(cè)向位移,橋臺地基不均勻沉降等,通常會引發(fā)橋梁結(jié)構(gòu)地震動放大,從而加劇結(jié)構(gòu)的地震破壞。因此,研究橋頭路堤、橋臺、地基的地震反應(yīng)及其相互作用特性,對公路橋梁的抗震加固設(shè)計很有實際意義。國內(nèi)不少外學(xué)者都很重視相關(guān)方面的研究,分別采用不同方法,從不同側(cè)重面進行了富有成效的工作,得到了很多有益結(jié)論。為了研究橋臺振動對橋梁結(jié)構(gòu)的影響,LI等[6]采用單跨梁橋模型,沿橋梁縱向激振進行振動臺試驗,研究橋臺固定方式對橋梁沖擊的影響。Anne等[7]選用美國加州典型的臺背回填材料粉砂進行橫向循環(huán)加載試驗,研究橋臺可犧牲背墻的土壓力特性,認為背墻位移為0.03(墻高)時,土壓力達到最大值,相應(yīng)被動土壓力系數(shù)為16.3;在深度處有一個對數(shù)螺旋主破壞面,被動楔塊中還有幾個相對較小的剪切面,可假定對數(shù)螺旋滑移破裂面,采用條分法估算土壓力;墻?土界面摩擦作用約為土體摩擦角的1/3到1/2。何度心[8]通過振動臺試驗研究橋頭路堤和橋臺的地震響應(yīng),發(fā)現(xiàn)路堤中孔隙水壓力增大,橋臺與臺后土體滑移,臺后填土裂縫,路面塌落,臺身滑移、傾斜、沉陷等。王建等[9]選取現(xiàn)場路堤填料,采用長3.5 m、寬1.5 m、高2 m的模型箱,輸入人工合成波進行試驗,結(jié)果表明路堤呈現(xiàn)上部拉裂、下部鼓脹的破壞性狀,印證了震害調(diào)查結(jié)果。為了估算橋臺的位移,蘇謙等[10]通過簡化橋臺滑動破壞模式,考慮土體強度衰減效應(yīng)對橋臺穩(wěn)定性的影響,提出了橋臺滑移計算公式。唐紅梅等[11]利用畢肖普法計算橋臺的地震位移,結(jié)合工程實例驗證方法的可靠性,并分析了黏聚力、振動孔隙水壓力和內(nèi)摩擦角等因素對位移的影響。Anoosh等[12]采用對數(shù)螺旋破壞面耦合修正的雙曲線土體應(yīng)力?應(yīng)變模型,用極限平衡法估算橋臺非線性位移,根據(jù)8個典型結(jié)構(gòu)回填土的現(xiàn)場試驗結(jié)果,驗證了計算方法的適應(yīng)性,并利用最小二乘法和試驗數(shù)據(jù),建立簡化的雙曲線力?位移方程,認為其計算結(jié)果用于橋臺抗震設(shè)計符合現(xiàn)有加州交通部門的要求。Al-Homoud等[13]進行了橋臺和填土非線性地震反應(yīng)分析,并考慮橋臺轉(zhuǎn)動和滑移條件,提出橋臺抗震設(shè)計方法,但計算分析中沒有考慮梁對橋臺的作用,并將其視為平面問題。孫治國等[14]利用ANSAYS建立了高原大橋三維有限元模型,研究橋臺破壞機理,并對抗震措施進行分析,認為胸墻、前墻以及其與翼墻交界處是地震中薄弱部位。這些研究從不同側(cè)面豐富和完善了關(guān)于橋頭路堤、橋臺、地基相互作用特性的內(nèi)涵,很有實際意義。由于橋頭段地震特性不但決定于橋臺、路堤、地基和錐坡自身振動特性,還與其相互作用特性密切相關(guān),地震反應(yīng)計算涉及到不同材料和不同結(jié)構(gòu)形式的組合,以及橋臺與路堤、地基以及錐坡的多重界面處理,考慮因素多,計算分析難度大,耗時長,有諸多方面還需進一步研究。本文基于FLAC3D計算平臺,考慮橋臺、路堤、地基和錐坡的相互作用,建立三維動力計算分析模型,研究橋臺?路堤?地基作用體系的地震反應(yīng)規(guī)律,旨在為相應(yīng)的抗震設(shè)計加固提供一定的理論借鑒。

    1 計算方法

    1.1 計算模型

    實際橋臺?路堤?地基作用體系包括橋臺、梁、地基、路堤與錐坡五部分,為了簡化計算,本文假定梁對橋臺僅傳遞豎向荷載,將梁與橋臺的作用簡化成豎向力,模型僅包括橋臺、地基、橋頭路堤和錐坡。

    幾何模型見圖1,其中路堤寬12 m,高6 m,路堤兩側(cè)邊坡坡度1:2。埋置式橋臺高8 m、厚2 m、寬12 m,擴大基礎(chǔ)埋深2 m,厚4 m,寬12 m。橋臺左右兩側(cè)與前側(cè)設(shè)置錐坡,坡度1:2,尺寸設(shè)置與公路工程技術(shù)標準典型斷面尺寸一致??紤]到模型邊界對計算精度的影響,經(jīng)試算,計算幅度確定為地基長80 m,寬48 m,厚20 m。為了波在網(wǎng)格中的合理傳播,網(wǎng)格單元尺寸均小于最小波長的1/10,鄰近橋臺的路堤和地基土體網(wǎng)格加密[1, 15-19],網(wǎng)格劃分見圖2。

    圖1 橋臺-路堤-地基作用體系示意

    圖2 計算模型

    1.2 材料參數(shù)

    相對土體,橋臺的剛度與強度大得多,為簡化計算,忽略橋臺的塑性變形,視為線彈性結(jié)構(gòu),土體為摩爾?庫倫材料,參數(shù)見表1。

    1.3 接觸面處理

    橋臺與地基、路堤以及錐坡的材料和結(jié)構(gòu)特性相差較大,其間存在滑移、分離,本文通過在橋臺與土體之間設(shè)置接觸面單元來考慮土?結(jié)構(gòu)界面作用,當(dāng)兩側(cè)的模量差異較大時,接觸面的法向剛度與切向剛度可取“較軟”材料等效剛度的10倍,按式1計算,由于接觸面有交叉,需要對其參數(shù)進行折減,接觸面摩擦角一般要小于填土的內(nèi)摩擦角,本模型取為10°,接觸面的具體參數(shù)見表2[16, 19]。

    表1 模型材料參數(shù)

    表2 接觸面參數(shù)

    1.4 邊界條件與阻尼設(shè)置

    靜力計算時,底部采用固定邊界,四周采用滑移支座約束水平位移,動力計算時,考慮波的反射,底部采用黏滯邊界,四周采用自由場邊界。

    阻尼通??刹捎萌鹄枘帷笞枘崤c局部阻尼等,Manica等[20]對比分析了3種阻尼的計算效果,表明瑞利阻尼得出的結(jié)果與實際的最吻合,因此本文采用瑞利阻尼。巖土工程中,土體達到塑性變形后會消耗一部分能量,因此阻尼比可取較小值,一般為2%~5%,本文取2%。

    1.5 地震荷載輸入

    研究時分別選取了Kobe波、Loma Prieta波、汶川波作為輸入波,為節(jié)約篇幅本文僅以汶川波為例,其他相應(yīng)分析詳見文獻[21]。

    由于三維動力分析模型計算耗時太長,為提高計算效率,選取汶川地震中臥龍臺站所記錄的東西向水平地震波,截取其中一個含波峰的,持時20 s的波段作為輸入波,并進行基線校正處理,以便線彈性結(jié)構(gòu)在動力計算結(jié)束后速度與位移為0[19]。處理后的汶川地震加速度時程與積分得到的位移時程見圖3。

    1.6 監(jiān)測點布置

    監(jiān)測點布置見圖4,從橋臺頂部向下每隔1 m布置1個監(jiān)測點,直至基礎(chǔ)底部,擴大基礎(chǔ)4個角另外各設(shè)1個監(jiān)測點,為闡述方便,將橋臺底部4個角編號為A,B,C和D。沿路線走向,遠離橋臺每1 m布置上、中、下3個監(jiān)控點,至距橋臺10 m處,每2 m布置監(jiān)測點,至距橋臺20 m處。基礎(chǔ)中心正下方的地基中,向下每2 m布置1個監(jiān) 測點。

    1.7 計算過程

    計算分為靜力計算與動力計算兩部分。先進行靜力計算,采用小變形模式進行靜力平衡,在施加梁對橋臺的豎向力時,為了避免產(chǎn)生過大的應(yīng)力集中,分5步進行加載。然后開啟大變形模式進行動力計算,將靜力計算后的位移場、速度場清零,保留應(yīng)力場,在模型底部施加水平方向地震荷載,采用動態(tài)時步模式。

    圖4 監(jiān)測點布置示意

    2 計算結(jié)果分析

    2.1 靜力反應(yīng)

    靜力平衡后,模型豎向應(yīng)力云圖見圖5,由圖可見,模型豎向應(yīng)力呈現(xiàn)明顯分層現(xiàn)象,隨深度增加豎向應(yīng)力逐漸增大,底部達到426.77 kPa。

    圖5 靜力計算豎向應(yīng)力云圖

    路堤對橋臺的土壓力分布見圖6,隨著深度增加,路堤土壓力非線性增大,呈內(nèi)凹拋物線形分布,路堤總土壓力165.313 kN。根據(jù)朗肯理論計算,路堤主動土壓力95.836 kN,被動土壓力876.266 kN,本文方法計算的靜止土壓力介于朗肯主動土壓力和被動土壓力之間。由于靜力平衡時,橋臺頂部位移為?0.624 mm,向路堤側(cè)偏移,推動回填土,因此,計算土壓力大于主動土壓力[16]。同時,從路堤土壓力分布看,其內(nèi)凹拋物線形與Springman等[22]的離心機試驗結(jié)果得到的分布規(guī)律吻合,由此說明本文計算模型和計算方法是合理的。

    圖6 震前土壓力分布

    橋臺基礎(chǔ)中心下方地基沉降沿深度變化見圖7,靠近橋臺底部的地基沉降較大,深處沉降相對較小。本文計算方法所得沉降與按《建筑地基基礎(chǔ)設(shè)計規(guī)范》(GB 50007—2011)[23]計算沉降的大小接近,最大差值僅為2.022 mm,沿深度的變化趨勢基本一致,都呈外凹拋物線形,進一步驗證了本文計算模型和計算方法的合理性。

    2.2 橋臺地震反應(yīng)

    輸入汶川波計算,橋臺頂部水平位移時程見圖8。

    為節(jié)約篇幅,橋臺其余各監(jiān)測點的位移時程參見文獻[21]。由圖可見,橋臺頂部位移呈波動變化,位移動力響應(yīng)很明顯,在第14.922 74 s時達到峰值位移312.526 mm后,隨著地震波主能量段的過去,橋臺發(fā)生部分回彈,出現(xiàn)67.762 mm的殘余水平位移。由于橋臺為線彈性,殘余水平位移主要由周圍彈塑性土體引起。賴杰等[24]采用汶川波激振,利用振動臺試驗研究抗滑樁的抗震性能,本文計算得到的橋臺位移時程曲線與其試驗監(jiān)測到的抗滑樁位移時程曲線特征類似。

    圖7 地基沉降

    圖8 橋臺頂部水平位移時程

    橋臺各監(jiān)測點峰值水平位移與殘余水平位移見圖9,橋臺殘余水平位移均為正值,表明橋臺震后整體向遠離路堤方向移動。同時,無論橋臺的峰值水平位移還是殘余位移,頂部與底部均不相同,說明橋臺發(fā)生了傾斜,表明地震過程中橋臺位移包含平移與轉(zhuǎn)動2種成分。橋臺殘余位移斜線與豎軸的夾角大于峰值位移斜線的夾角,說明地震結(jié)束時橋臺的傾斜角度大于峰值時刻,橋臺傾斜角度隨著震動持時而累積。

    橋臺擴大基礎(chǔ)底部4個角的豎向位移時程見圖10,A點與D點,B點與C點的豎向位移時程曲線基本一致,這是由于它們分別關(guān)于道路中線對稱的原因。橋臺沉降整體呈波動增大趨勢,A點在18.236 s時達到最大值14.957 mm,地震結(jié)束時沉降為14.604 mm,B點在17.328 s時達到最大值20.153 mm,最終沉降為19.493 mm,A,B相差4.889 mm,說明橋臺發(fā)生旋轉(zhuǎn),B點最終沉降大于A點,表明橋臺向外傾斜,不均勻沉降隨著地震持時而增大,即傾斜角度隨時間逐漸增大。

    圖9 橋臺峰值水平位移與殘余水平位移

    圖10 橋臺豎向位移時程

    橋臺峰值加速度分布見圖11,由圖11可見,橋臺峰值加速度與加速度放大系數(shù)沿橋臺高度整體呈增大趨勢,橋臺底部的峰值加速度大于距橋臺頂部5~7 m處的峰值加速度,這是因為該處是橋臺與地基交界面,受地震作用較強,故對應(yīng)峰值加速度較大。峰值加速度增大幅度隨著高度的增大而增大,由底部的5.948 m/s2增加至頂部的10.520 m/s2,整體上,加速度沿橋臺豎向上存在放大效應(yīng),而且剛開始增大緩慢,到達一定高度時,峰值加速度急劇增大,加速度放大系數(shù)達到2.630。

    圖11 橋臺加速度峰值

    路堤土壓力沿橋臺高度變化見圖12,由圖可見,地震土壓力最大值出現(xiàn)在距橋臺頂部4.5 m處;本文計算的路堤地震峰值土壓力與物部-岡部法(M-O法)以及規(guī)范方法的計算結(jié)果對比見表3,路堤地震峰值土壓力為353.588 kN,M-O法計算值為363.573 kN,公路規(guī)范為295.262 kN,鐵路規(guī)范為268.239 kN,本文計算路堤峰值土壓力大于公路規(guī)范與鐵路規(guī)范,與M-O法計算結(jié)果接近。

    圖12 土壓力分布

    表3 地震土壓力

    2.3 橋頭路堤地震反應(yīng)

    與臺背相鄰處路堤表面的水平位移時程見圖13,其余各點的位移時程曲線可參見文獻[21]。由圖可見,水平位移最大值為312.485 mm,殘余位移為67.727 mm;相同高度處橋臺的峰值水平位移為312.526 mm,殘余水平位移為67.762 mm,橋臺與路堤相鄰處位移不同,存在分離現(xiàn)象,這與震后現(xiàn)場橋頭路堤與橋臺常出現(xiàn)裂縫現(xiàn)象相吻合。

    圖13 臺背相鄰處路堤表面水平位移時程

    路堤表面各監(jiān)測點峰值水平位移與殘余水平位移見圖14,由圖14可見,距臺背0,5,10和20 m處路堤表面峰值水平位移分別為312.495,306.691,305.321和303.625 mm,殘余水平位移分別為67.727,60.404,58.787和56.948 mm,路堤表面的峰值水平位移與殘余水平位移隨遠離臺背而減小,可見路堤有拉裂趨勢。

    距臺背0,5和10 m處路堤表面的豎向位移時程曲線見圖15,由圖15可見,隨地震作用持續(xù),各點豎向位移逐漸增大,但距臺背0,10和20 m處各點豎向位移不同。

    路堤表面最終沉降見圖16,由圖16可見,最終沉降隨著距臺背距離的增加呈減小趨勢。從路堤的水平位移與豎向位移可發(fā)現(xiàn),路堤有拉裂的趨勢,橋頭段路堤變形較大,呈塌滑趨勢,是地震中的薄弱點,而且橋頭路堤已經(jīng)與橋臺發(fā)生分離,這些變形和破壞特征與震后調(diào)查結(jié)果一致[2?4]。

    圖15 路堤表面豎向位移時程

    圖16 路堤表面最終沉降

    2.4 地基地震反應(yīng)

    與橋臺底部相鄰處地基的水平位移時程見圖17,水平位移在14.899 59 s達到峰值307.279 mm,殘余水平位移為59.472 mm,而橋臺底部峰值水平位移為307.186 mm,殘余水平位移為60.132 mm,可見,在地震過程中,橋臺底部與地基產(chǎn)生相對滑移。地基峰值水平位移與殘余水平位移見圖18,地基峰值水平位移與殘余水平位移隨著深度的增加而減小,地基底部殘余水平位移為43.247 mm,為頂部的72.718%。

    圖17 橋臺相鄰處地基土水平位移時程

    圖18 地基峰值水平位移與殘余水平位移

    距橋臺底部0,10和18 m處的地基豎向位移時程見圖19,距橋臺底部0 m處峰值豎向位移?17.216 mm,殘余位移?16.976 mm;距橋臺底部10 m處峰值豎向位移?8.189 mm,殘余位移?7.689 mm;距橋臺底部18 m處峰值豎向位移為?1.592 mm,殘余位移為?1.068 mm。

    圖19 地基豎向位移時程

    圖20 基底剪力時程

    基底剪力時程見圖20,基底剪力在地震中呈波動增加趨勢,在第8.678 s達到最大值295.046 kN,相對震前剪力59.505 kN,增大幅度達206.6%,地震結(jié)束時的殘余剪力為106.174 kN。從前文已知存在相對滑移,抗震設(shè)計時應(yīng)進行抗滑移驗算。

    3 結(jié)論

    1) 本文計算臺背土壓力處于朗肯主動土壓力與被動土壓力之間,分布特征與離心機試驗結(jié)果相吻合;計算的地基沉降與規(guī)范方法計算結(jié)果接近,變化趨勢一致;計算的地震峰值土壓力與M-O法計算結(jié)果接近;計算得到的路堤和橋臺地震反應(yīng)特征和破壞趨勢與震后調(diào)查資料吻合良好,說明本文數(shù)值分析模型和計算方法較合理。

    2) 震后橋臺整體遠離路堤滑動,并呈整體下沉和外傾狀態(tài),地震位移包含平移與轉(zhuǎn)動;橋臺加速度在高度方向上存在放大效應(yīng),承受的地震峰值土壓力與最大彎矩較震前明顯增加,公路規(guī)范與鐵路規(guī)范計算值偏小。

    3) 路堤震后整體向橋臺方向推移,路堤表面峰值以及殘余水平和豎向位移反應(yīng)特征表明,橋頭路堤在地震中易發(fā)生拉裂和不均勻沉降破壞,反映了橋頭段的主要地震破壞特征。

    4) 毗鄰橋臺的地基土體的峰值水平位移與殘余水平位移以及最終沉降隨較大;地震過程中基底剪力較震前有較大幅度增加,存在相對滑移。

    [1] 陳國興. 巖土地震工程學(xué)[M]. 北京: 科學(xué)出版社, 2007. CHEN Guoxing. Geotechnical and seismic engineering [M]. Beijing: Science Press, 2007.

    [2] 王東升, 郭迅, 孫治國, 等. 汶川大地震公路橋梁震害初步調(diào)查[J]. 地震工程與工程振動, 2009, 29(3): 84?94. WANG Dongsheng, GUO Xun, SUN Zhiguo, et al. Damage to highway bridges during Wenchuan earth -quake[J]. Journal of Earthquake Engineering and Engineering Vibration, 2009, 29(3): 84?94.

    [3] 莊衛(wèi)林, 陳樂生. 汶川地震公路震害分析(橋梁與隧道) [M]. 北京: 人民交通出版社, 2013. ZHUANG Weilin, CHEN Lesheng. Analysis of highway earthquake damage in Wenchuan earthquake (Bridges and tunnels)[M]. Beijing: Communication Press, 2013.

    [4] 杜修力, 韓強, 李忠獻, 等. 汶川地震中山區(qū)公路橋梁震害及啟示[J]. 北京工業(yè)大學(xué)學(xué)報, 2008, 34(12): 1270?1279. DU Xiuli, HAN Qiang, LI Zhongxian, et al. Earthquake damage of highway bridges in mountainous areas and its enlightenment in 512 Wenchuan earthquake[J]. Journal of Beijing University of Technology, 2008, 34(12): 1270? 1279.

    [5] 莊衛(wèi)林, 陳樂生. 汶川地震公路震害分析(地質(zhì)災(zāi)害與路基)[M]. 北京: 人民交通出版社, 2013. ZHUANG Weilin, CHEN Lesheng. Analysis of highway earthquake damage in Wenchuan earthquake (Geological hazards and roadbed)[M]. Beijing: Communication Press, 2013.

    [6] LI B, BI K, Chouw N, et al. Effect of abutment excitation on bridge pounding[J]. Engineering Structures, 2013, 54(3): 57?68.

    [7] Anne Lemnitzer S, Eric R, Ahlberg A, et al. Lateral performance of full-scale bridge abutment wall with granular backfill[J]. J Geotech Geoenviron Eng, 2009, 135(4): 506?514.

    [8] 何度心. 橋臺抗震設(shè)計的理論與方法[J]. 世界地震工程, 1991, 50(3): 1?15. HE Duxin. Theory and method of seismic design of bridge abutment[J]. World Earthquake Engineering, 1991, 50(3): 1?15.

    [9] 王建, 姚令侃, 吳偉. 路堤震害模式及路堤動力特性研究[J]. 巖土力學(xué), 2010, 31(12): 3801?3808. WANG Jian, YAO Lingkan, WU Wei. Research on seismic damage mode and dynamic characteristics of road embankment[J]. Rock and Soil Mechanics, 2010, 31(12): 3801?3808.

    [10] 蘇謙, 鐘彪, 白皓, 等. 基于滑坡破壞模式的橋臺位移耦合計算[J]. 土木工程學(xué)報, 2010, 43(增1): 567?572. SU Qian, ZHONG Biao, BAI Hao, et al. Coupled calculation of abutment displacement based on landslide failure mode[J]. China Civil Engineering Journal, 2010, 43(Suppl 1): 567?572.

    [11] 唐紅梅, 劉榮燦, 陳洪凱. 地震作用下橋臺滑移計算方法及應(yīng)用[J]. 公路, 2010, 20(3): 42?47. TANG Hongmei, LIU Rongcan, CHEN Hongkai. Calculation method and application of bridge abutment slip under earthquake action[J]. Highway, 2010, 20(3): 42?47.

    [12] Anoosh Shamsabadi, Kyle M Rollins, Mike kapuskar. Nonlinear soil-abutment–bridge structure interaction for seismic performance-based design[J]. J Geotech Geoenviron Eng, 2007, 133(6): 707?720.

    [13] Al-Homoud A, Whitman R. Seismic analysis and design of rigid bridge abutments considering rotation and sliding incorporating non-linear soil behavior[J]. Soil Dynamics and Earthquake Engineering, 1999, 18(4): 247?277.

    [14] 孫治國, 王東升, 張蓓, 等. 高原大橋橋臺地震破壞機理與抗震措施分析[J]. 地震工程與工程振動, 2012, 32(4): 79?87. SUN Zhiguo, WANG Dongsheng, ZHANG Bei, et al. Analysis of seismic failure mechanism and earthquake resistant measures for the abutments of Gaoyuan[J]. Journal of Earthquake Engineering and Engineering Vibration, 2012, 32(4): 79?87.

    [15] 祁生文, 伍法權(quán),嚴福章, 等. 巖質(zhì)邊坡土動力反應(yīng)分析[M]. 北京: 科學(xué)出版社, 2007. QI Shengwen, WU Faquan,YAN Fuzhang, et al. Rock slope dynamic response analysis[M]. Beijing: Science Press, 2007.

    [16] 錢家歡, 殷宗澤. 土工原理與計算[M]. 北京: 中國水利水電出版社, 1996. QIAN Jiahuan, YIN Zongze. Geotechnical principles and computation[M]. Beijing: China Water Resources and Hydropower Press, 1996.

    [17] 吳世明. 土動力學(xué)[M]. 北京: 中國建筑工業(yè)出版社, 2000. WU Shiming. Soil dynamics[M]. Beijing: China Construction Industry Press, 2000.

    [18] Francisco M, Fernando R. Spurious reflections and accuracy of finite element models for unbounded wave propagation problems: the one-dimensional case[J]. Engineering Computations, 1987, 4(2):139?148.

    [19] 陳育民, 徐鼎平. FLAC/FLAC3D基礎(chǔ)與工程實例[M]. 北京: 中國水利水電出版社, 2013.

    CHEN Yumin, XU Dingping. FLAC/FLAC3Dfoundation and engineering example[M]. Beijing: China Water Resources and Hydropower Press, 2013

    [20] Manica M, Ovando E. Assessment of damping models in FLAC[J]. Computers & Geotechnics, 2014, 59(6): 12?20.

    [21] 程志明. 橋臺?路堤?地基作用體系的地震反應(yīng)分析[D]. 長沙: 中南大學(xué), 2019. CHENG Zhiming. Analysis of abutment-embankment- foundation interaction system subjected to seismic load [D]. Changsha: Central South University, 2019.

    [22] Springman S, Ng C, Ellis E. Centrifuge and analytical studies of full height bridge abutment on piled foundation subjected to lateral loading[R]. CUED/D-SOILS/TR278. 1994: 1?120.

    [23] GB50007—2011, 建筑地基基礎(chǔ)設(shè)計規(guī)范[S]. GB50007—2011, Code for design of building foundation [S].

    [24] 賴杰, 鄭穎人, 劉云, 等. 地震作用下雙排抗滑樁支護邊坡振動臺試驗研究[J]. 巖土工程學(xué)報, 2014, 36(4): 680?687. LAI Jie, ZHENG Yingren, LIU Yun, et al. Shaking table tests on double-row anti-slide piles of slopes under earthquakes[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(4): 680?687.

    Analysis of abutment-embankment-foundation interaction system subjected to seismic load

    CHENG Zhiming1, 2, WEI Hongwei1

    (1. School of Civil Engineering, Central South University, Changsha 410075, China;2. Wuhan Municipal Engineering Design and Research Institute Co., Ltd., Wuhan 430023, China)

    A three-dimensional non-linear numerical model of approach embankment was established considering abutment- embankment-foundation interaction, and the Wenchuan earthquake wave was used as seismic load to carry out the non-linear time history analyses of seismic response. Comparing the results of numerical calculation with those of classical theory and standard method, centrifugal test results and seismic damage investigation data, the rationality of the calculation model and method was verified. The main analyses included the seismic displacement responses of embankment, abutment and foundation, seismic accelerations, earth pressure, bending moments and base shear stress response of abutment. The results show that the abutment deviated from the embankment, accompanied with sinking and extroversion after earthquake, and the seismic displacements included translation and rotation. The acceleration at the top of the abutment appeared amplification effect, and the seismic earth pressure and bending moments were obviously increased compared with those before earthquake. Comparatively, the calculation results by highway code and railway code method appeared to be smaller. After earthquake, the main seismic damage behaviors in approach embankment were crack and differential settlement; The peak and residual horizontal displacement, and final settlement of foundation soil adjacent to abutment were relatively large, and the shear force between the foundation and abutment increased greatly in earthquake which led to a relative slip. The results in this paper will be beneficial to seismic reinforcement design of bridges.

    embankment; abutment; foundation; seismic response; time history analysis

    TU416.1

    A

    1672 ? 7029(2020)09 ? 2225 ? 10

    10.19713/j.cnki.43?1423/u. T20191094

    2019?12?06

    國家自然科學(xué)基金面上項目(51778639,51678573)

    魏紅衛(wèi)(1966?),男,河南汝州人,教授,博士,從事結(jié)構(gòu)工程與巖土工程的教學(xué)和研究工作;E?mail:whw.edu@163.com

    (編輯 涂鵬)

    猜你喜歡
    橋頭橋臺路堤
    土工合成材料加筋土橋臺技術(shù)及其應(yīng)用
    路堤下CFG樁復(fù)合地基穩(wěn)定分析方法探討
    解析橋梁下部結(jié)構(gòu)樁柱橋臺施工
    江西建材(2018年1期)2018-04-04 05:26:20
    橋頭跳車研究
    中國公路(2017年15期)2017-10-16 01:32:08
    橋頭跳車現(xiàn)象防治措施
    中國公路(2017年12期)2017-02-06 03:07:37
    多年凍土區(qū)鐵路路堤臨界高度研究
    橋頭牧歌
    仁存溝大橋高邊坡橋臺病害治理
    煤矸石浸水路堤不均勻沉降研究
    無處理軟基橋臺加固及驗算
    在现免费观看毛片| 最近最新中文字幕大全电影3| 在线免费观看不下载黄p国产 | 免费av不卡在线播放| 国产精品综合久久久久久久免费| avwww免费| a级毛片a级免费在线| 一本一本综合久久| 天堂网av新在线| 欧美最新免费一区二区三区| 他把我摸到了高潮在线观看| 亚洲精品影视一区二区三区av| 性欧美人与动物交配| 窝窝影院91人妻| 欧美三级亚洲精品| 99久久中文字幕三级久久日本| 一区二区三区四区激情视频 | 春色校园在线视频观看| 国产aⅴ精品一区二区三区波| 成人国产综合亚洲| 久久这里只有精品中国| 欧美成人a在线观看| 久久久国产成人精品二区| 亚洲中文日韩欧美视频| 日日啪夜夜撸| 深爱激情五月婷婷| 成人二区视频| 亚洲国产精品合色在线| avwww免费| 91在线精品国自产拍蜜月| 在线看三级毛片| 日韩一区二区视频免费看| 国产精品久久久久久久电影| 在线免费观看的www视频| 一区二区三区免费毛片| 久久人妻av系列| 一级黄色大片毛片| 我的女老师完整版在线观看| 高清毛片免费观看视频网站| 国产色爽女视频免费观看| 欧美日韩综合久久久久久 | 两人在一起打扑克的视频| 国产激情偷乱视频一区二区| 国产精品综合久久久久久久免费| 亚洲精品国产成人久久av| av女优亚洲男人天堂| 久久精品国产亚洲av天美| 国内精品一区二区在线观看| 国产主播在线观看一区二区| 亚洲va在线va天堂va国产| 亚洲不卡免费看| 国产探花极品一区二区| 久久精品久久久久久噜噜老黄 | 动漫黄色视频在线观看| 99热只有精品国产| 成人av在线播放网站| 一本精品99久久精品77| 悠悠久久av| 国产av不卡久久| 欧美色视频一区免费| 国产精品一区二区免费欧美| 成人av一区二区三区在线看| 久久久精品欧美日韩精品| 中文字幕精品亚洲无线码一区| 美女被艹到高潮喷水动态| 12—13女人毛片做爰片一| 十八禁网站免费在线| 桃红色精品国产亚洲av| 欧美潮喷喷水| 亚洲精品一区av在线观看| 舔av片在线| 99热这里只有是精品在线观看| 免费av观看视频| 老熟妇仑乱视频hdxx| 精品久久国产蜜桃| 国产亚洲欧美98| 国产成人一区二区在线| 久久精品国产亚洲av香蕉五月| 免费搜索国产男女视频| 日本一本二区三区精品| 99久久久亚洲精品蜜臀av| 亚洲熟妇熟女久久| 夜夜夜夜夜久久久久| 日日摸夜夜添夜夜添av毛片 | 中文字幕免费在线视频6| 国产免费男女视频| 国产高清激情床上av| 亚洲五月天丁香| 久久久久久久久久久丰满 | 午夜亚洲福利在线播放| 久久精品国产99精品国产亚洲性色| 亚洲成av人片在线播放无| 精品人妻1区二区| 最近视频中文字幕2019在线8| 韩国av在线不卡| 亚洲国产日韩欧美精品在线观看| 少妇丰满av| 日本 欧美在线| av黄色大香蕉| 99久久久亚洲精品蜜臀av| 天堂√8在线中文| 亚洲成人久久爱视频| 午夜免费男女啪啪视频观看 | 国产成年人精品一区二区| 久久久久久久午夜电影| 我要搜黄色片| 丝袜美腿在线中文| 别揉我奶头~嗯~啊~动态视频| 麻豆久久精品国产亚洲av| av国产免费在线观看| av中文乱码字幕在线| 国产精品人妻久久久影院| 午夜免费激情av| 久久人妻av系列| 免费在线观看影片大全网站| 我要搜黄色片| 国产aⅴ精品一区二区三区波| av在线蜜桃| 亚洲精品粉嫩美女一区| 91精品国产九色| 看十八女毛片水多多多| 波多野结衣高清作品| 国产成人aa在线观看| 波多野结衣巨乳人妻| 最好的美女福利视频网| 久久久久国内视频| 看片在线看免费视频| 久久久久久九九精品二区国产| 亚洲精品国产成人久久av| 男人狂女人下面高潮的视频| 国产视频一区二区在线看| 亚洲av免费在线观看| 亚洲四区av| 真人做人爱边吃奶动态| 精品国产三级普通话版| 免费观看的影片在线观看| 99热精品在线国产| 亚洲成人中文字幕在线播放| 国内精品宾馆在线| 韩国av一区二区三区四区| www日本黄色视频网| 夜夜爽天天搞| 国产精品一区二区免费欧美| 欧美一区二区精品小视频在线| 丰满乱子伦码专区| 亚洲av电影不卡..在线观看| 久久香蕉精品热| 欧美精品国产亚洲| 在线免费观看不下载黄p国产 | www.www免费av| 精品国内亚洲2022精品成人| 精华霜和精华液先用哪个| 一本久久中文字幕| 不卡视频在线观看欧美| 国产黄a三级三级三级人| 很黄的视频免费| 黄色一级大片看看| 亚洲av日韩精品久久久久久密| 性欧美人与动物交配| 精品不卡国产一区二区三区| 欧美性感艳星| 天堂网av新在线| 久久人人精品亚洲av| 欧美成人a在线观看| 免费在线观看日本一区| 午夜福利欧美成人| 日韩 亚洲 欧美在线| 亚洲国产精品合色在线| 欧美性猛交黑人性爽| 欧美成人免费av一区二区三区| av.在线天堂| or卡值多少钱| 久久国产精品人妻蜜桃| 亚洲精品日韩av片在线观看| 性色avwww在线观看| 麻豆一二三区av精品| 久久精品国产自在天天线| 欧美xxxx黑人xx丫x性爽| 精品国内亚洲2022精品成人| 黄色女人牲交| xxxwww97欧美| 日本黄大片高清| 制服丝袜大香蕉在线| 天堂影院成人在线观看| 少妇人妻一区二区三区视频| 国产真实乱freesex| 欧美+日韩+精品| 精品乱码久久久久久99久播| 国产精品久久久久久久久免| 又粗又爽又猛毛片免费看| 亚洲aⅴ乱码一区二区在线播放| 亚洲人成网站高清观看| 深夜a级毛片| 免费观看的影片在线观看| 国产成人一区二区在线| 床上黄色一级片| 成年女人毛片免费观看观看9| 熟女电影av网| 91久久精品国产一区二区成人| 久久久久久久亚洲中文字幕| 久久人人精品亚洲av| 高清毛片免费观看视频网站| 久久久久国产精品人妻aⅴ院| 国产成年人精品一区二区| 成人鲁丝片一二三区免费| 国产精品久久久久久久电影| 免费人成视频x8x8入口观看| 亚洲成av人片在线播放无| 一进一出好大好爽视频| 国产v大片淫在线免费观看| 999久久久精品免费观看国产| 深爱激情五月婷婷| 久久国产精品人妻蜜桃| 亚洲国产精品sss在线观看| 在现免费观看毛片| 精品不卡国产一区二区三区| 女人十人毛片免费观看3o分钟| 大又大粗又爽又黄少妇毛片口| 久久国产乱子免费精品| 韩国av在线不卡| 熟女电影av网| 丝袜美腿在线中文| 国产精品亚洲美女久久久| 天堂动漫精品| 日本在线视频免费播放| 国产 一区精品| 国内精品宾馆在线| 国产老妇女一区| 亚洲性夜色夜夜综合| 中文字幕人妻熟人妻熟丝袜美| 又粗又爽又猛毛片免费看| 欧美xxxx黑人xx丫x性爽| 中文资源天堂在线| 男女边吃奶边做爰视频| 欧美一级a爱片免费观看看| 亚洲欧美日韩无卡精品| 热99re8久久精品国产| 三级男女做爰猛烈吃奶摸视频| 淫秽高清视频在线观看| 成人一区二区视频在线观看| 我的女老师完整版在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久国产乱子免费精品| 亚洲成a人片在线一区二区| 亚洲电影在线观看av| 亚洲av第一区精品v没综合| 国产av麻豆久久久久久久| 黄色女人牲交| 精品99又大又爽又粗少妇毛片 | 亚洲不卡免费看| 国产男靠女视频免费网站| 亚洲黑人精品在线| 特大巨黑吊av在线直播| 国产一区二区三区在线臀色熟女| 中国美女看黄片| 一级黄片播放器| 中国美白少妇内射xxxbb| 国产一区二区三区av在线 | 欧美日韩瑟瑟在线播放| 少妇人妻一区二区三区视频| 欧美日韩中文字幕国产精品一区二区三区| 国产伦人伦偷精品视频| 三级毛片av免费| 18+在线观看网站| 欧美又色又爽又黄视频| 欧美一区二区亚洲| 91精品国产九色| h日本视频在线播放| 丰满的人妻完整版| 久久午夜亚洲精品久久| 99精品在免费线老司机午夜| 久久久久久九九精品二区国产| 噜噜噜噜噜久久久久久91| 久久精品国产99精品国产亚洲性色| 又黄又爽又免费观看的视频| 国产精品av视频在线免费观看| 久久久久九九精品影院| 99精品在免费线老司机午夜| 2021天堂中文幕一二区在线观| 久久久久久久精品吃奶| 国产成人a区在线观看| 午夜免费成人在线视频| 亚洲黑人精品在线| 成人美女网站在线观看视频| 精品一区二区免费观看| 人妻久久中文字幕网| 亚洲真实伦在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美日韩综合久久久久久 | 午夜日韩欧美国产| 美女黄网站色视频| 国产中年淑女户外野战色| 欧美高清性xxxxhd video| 可以在线观看的亚洲视频| 免费一级毛片在线播放高清视频| 亚洲真实伦在线观看| 精品免费久久久久久久清纯| 日韩欧美在线二视频| 亚洲第一区二区三区不卡| 国产成人福利小说| 女的被弄到高潮叫床怎么办 | 国产一区二区三区在线臀色熟女| 亚洲av不卡在线观看| 国产熟女欧美一区二区| 亚洲内射少妇av| 啦啦啦啦在线视频资源| 麻豆成人av在线观看| 九色成人免费人妻av| 一级av片app| 成人特级黄色片久久久久久久| or卡值多少钱| 免费av观看视频| 男女边吃奶边做爰视频| 国产在线男女| a级毛片a级免费在线| 99精品在免费线老司机午夜| 麻豆国产97在线/欧美| 精品久久久噜噜| 黄色女人牲交| 国产在视频线在精品| 男女下面进入的视频免费午夜| 久久精品国产99精品国产亚洲性色| 欧美又色又爽又黄视频| 久久久久性生活片| 日日啪夜夜撸| 夜夜夜夜夜久久久久| 此物有八面人人有两片| 乱系列少妇在线播放| 亚洲第一区二区三区不卡| 18禁裸乳无遮挡免费网站照片| 婷婷精品国产亚洲av| 日韩av在线大香蕉| 亚洲av五月六月丁香网| 国产精品久久久久久av不卡| 99久久九九国产精品国产免费| 免费看美女性在线毛片视频| 色5月婷婷丁香| 舔av片在线| 在线国产一区二区在线| 99九九线精品视频在线观看视频| 色5月婷婷丁香| 乱人视频在线观看| 熟妇人妻久久中文字幕3abv| 国产伦人伦偷精品视频| 欧美激情久久久久久爽电影| 久久亚洲真实| 精品不卡国产一区二区三区| 免费观看的影片在线观看| 国产熟女欧美一区二区| 亚洲av免费高清在线观看| 久久精品91蜜桃| 免费观看的影片在线观看| 成人特级黄色片久久久久久久| a在线观看视频网站| 精品不卡国产一区二区三区| 男插女下体视频免费在线播放| 99视频精品全部免费 在线| 天天躁日日操中文字幕| 免费av不卡在线播放| 中亚洲国语对白在线视频| 国产精品1区2区在线观看.| 少妇猛男粗大的猛烈进出视频 | 天天一区二区日本电影三级| 69av精品久久久久久| 欧美激情久久久久久爽电影| 久久天躁狠狠躁夜夜2o2o| 国产白丝娇喘喷水9色精品| 狠狠狠狠99中文字幕| 内射极品少妇av片p| 一个人免费在线观看电影| 亚洲性久久影院| 免费观看人在逋| 国产探花极品一区二区| 男女视频在线观看网站免费| 亚洲专区国产一区二区| 国内精品宾馆在线| 日韩一本色道免费dvd| 乱系列少妇在线播放| 在线观看av片永久免费下载| 精品无人区乱码1区二区| 精品人妻1区二区| 欧美丝袜亚洲另类 | 九九爱精品视频在线观看| 精品国产三级普通话版| 老女人水多毛片| 国产91精品成人一区二区三区| 精品久久久久久久久久久久久| 国产精品一及| 亚洲三级黄色毛片| 狂野欧美激情性xxxx在线观看| 在线国产一区二区在线| 偷拍熟女少妇极品色| 国产探花极品一区二区| 亚洲五月天丁香| 99视频精品全部免费 在线| 午夜福利视频1000在线观看| 午夜福利高清视频| 久久久久性生活片| 精品久久久久久久久亚洲 | 亚洲国产欧洲综合997久久,| 赤兔流量卡办理| 乱人视频在线观看| 国内精品一区二区在线观看| 午夜老司机福利剧场| 精品一区二区三区视频在线观看免费| 久99久视频精品免费| 一进一出抽搐gif免费好疼| 三级毛片av免费| 午夜福利18| 色av中文字幕| 久久精品国产清高在天天线| 中文字幕久久专区| 我要看日韩黄色一级片| 国产毛片a区久久久久| 欧美成人a在线观看| 99视频精品全部免费 在线| 亚洲av免费高清在线观看| 国产爱豆传媒在线观看| 欧美3d第一页| 搡女人真爽免费视频火全软件 | 欧美高清成人免费视频www| 超碰av人人做人人爽久久| 两性午夜刺激爽爽歪歪视频在线观看| 一个人观看的视频www高清免费观看| 麻豆成人av在线观看| 久久午夜亚洲精品久久| 此物有八面人人有两片| 婷婷丁香在线五月| 亚洲av.av天堂| 国内精品久久久久精免费| 深夜精品福利| 在线播放无遮挡| 两个人视频免费观看高清| 亚洲专区国产一区二区| 国产在线精品亚洲第一网站| 国产一区二区亚洲精品在线观看| 91狼人影院| 亚洲七黄色美女视频| 波野结衣二区三区在线| 1024手机看黄色片| 日本与韩国留学比较| 色av中文字幕| 88av欧美| 韩国av一区二区三区四区| 久久6这里有精品| 国产精品亚洲一级av第二区| 国产精品一区www在线观看 | 国产又黄又爽又无遮挡在线| 有码 亚洲区| 伦精品一区二区三区| 免费搜索国产男女视频| av天堂中文字幕网| 1024手机看黄色片| 亚洲乱码一区二区免费版| 国产精品亚洲一级av第二区| 给我免费播放毛片高清在线观看| 神马国产精品三级电影在线观看| 精品午夜福利视频在线观看一区| 久久久久久伊人网av| 天堂网av新在线| 国产精品自产拍在线观看55亚洲| 国产黄片美女视频| 99精品久久久久人妻精品| 国产三级在线视频| 尤物成人国产欧美一区二区三区| 成年女人永久免费观看视频| 欧洲精品卡2卡3卡4卡5卡区| 3wmmmm亚洲av在线观看| 成熟少妇高潮喷水视频| 搡老岳熟女国产| 国产单亲对白刺激| 波多野结衣高清无吗| 简卡轻食公司| 三级国产精品欧美在线观看| 免费观看的影片在线观看| 日本黄色视频三级网站网址| 久久久久国内视频| .国产精品久久| 中亚洲国语对白在线视频| 国产黄片美女视频| 夜夜爽天天搞| 亚洲经典国产精华液单| 亚洲国产精品成人综合色| 欧美人与善性xxx| 免费大片18禁| 午夜免费男女啪啪视频观看 | 亚洲av熟女| 日韩,欧美,国产一区二区三区 | 欧美日韩亚洲国产一区二区在线观看| 国产精品无大码| 美女高潮的动态| 亚洲国产精品sss在线观看| 成熟少妇高潮喷水视频| 欧美精品国产亚洲| 免费看光身美女| 久久久久九九精品影院| 亚洲精品一区av在线观看| 国内久久婷婷六月综合欲色啪| 国产一区二区在线av高清观看| 老司机深夜福利视频在线观看| 一边摸一边抽搐一进一小说| 亚洲 国产 在线| 精品久久久久久成人av| 精品久久久久久久久av| 亚洲欧美日韩东京热| 春色校园在线视频观看| 亚洲精品一区av在线观看| 99久久成人亚洲精品观看| 欧美一区二区国产精品久久精品| 少妇高潮的动态图| 69人妻影院| 动漫黄色视频在线观看| 性欧美人与动物交配| 免费不卡的大黄色大毛片视频在线观看 | 国产 一区 欧美 日韩| 成人特级av手机在线观看| 国产毛片a区久久久久| 精品一区二区三区视频在线| 国产一区二区在线观看日韩| 成人欧美大片| 91久久精品国产一区二区成人| 99久久久亚洲精品蜜臀av| 亚洲av不卡在线观看| 国产91精品成人一区二区三区| 亚洲av成人精品一区久久| 色哟哟·www| 欧美日本视频| 露出奶头的视频| 91在线观看av| 91久久精品国产一区二区三区| 波多野结衣巨乳人妻| 欧美性猛交╳xxx乱大交人| 久久久久久久久久成人| 日本a在线网址| 精品福利观看| 亚洲精华国产精华精| 日韩中文字幕欧美一区二区| 成年人黄色毛片网站| 欧美一区二区精品小视频在线| 在线免费观看不下载黄p国产 | 国产69精品久久久久777片| av天堂在线播放| 啦啦啦啦在线视频资源| 久久精品久久久久久噜噜老黄 | 在线天堂最新版资源| 国产欧美日韩一区二区精品| 亚洲av中文av极速乱 | 美女黄网站色视频| 国产免费男女视频| 精品一区二区三区视频在线| 午夜影院日韩av| 欧美高清性xxxxhd video| 国产精品国产三级国产av玫瑰| 精品国内亚洲2022精品成人| 成年女人永久免费观看视频| 在线国产一区二区在线| 亚洲av一区综合| 欧美日韩瑟瑟在线播放| 中亚洲国语对白在线视频| 日韩中文字幕欧美一区二区| 麻豆成人午夜福利视频| 亚洲精品在线观看二区| 国产精品伦人一区二区| 热99在线观看视频| 乱码一卡2卡4卡精品| 麻豆国产97在线/欧美| 18+在线观看网站| 99国产精品一区二区蜜桃av| 国产视频一区二区在线看| 久9热在线精品视频| 97碰自拍视频| 99久国产av精品| 国产精品久久久久久精品电影| 天堂网av新在线| 国产精品自产拍在线观看55亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 免费搜索国产男女视频| 亚洲国产高清在线一区二区三| 99精品久久久久人妻精品| 人人妻人人澡欧美一区二区| 一卡2卡三卡四卡精品乱码亚洲| 久久99热这里只有精品18| 亚洲欧美日韩高清在线视频| 九九热线精品视视频播放| 国产精品亚洲美女久久久| 免费搜索国产男女视频| 欧美三级亚洲精品| 无人区码免费观看不卡| av国产免费在线观看| 又黄又爽又刺激的免费视频.| a级毛片免费高清观看在线播放| 麻豆成人av在线观看| 免费看日本二区| 国产精品日韩av在线免费观看| 亚洲色图av天堂| 日本黄色片子视频| 精品久久久久久久久av| 国产精品亚洲一级av第二区| 蜜桃亚洲精品一区二区三区| 成年女人永久免费观看视频| 在线天堂最新版资源| 成人亚洲精品av一区二区| 国产一区二区三区在线臀色熟女| 日本黄大片高清| 色尼玛亚洲综合影院| 午夜亚洲福利在线播放| 大型黄色视频在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 日本a在线网址| 97碰自拍视频| 午夜精品一区二区三区免费看| 亚洲精华国产精华精| 亚洲精品日韩av片在线观看| 国产在线男女| 午夜精品在线福利|