• <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ū)鐵路路堤臨界高度研究
    橋頭牧歌
    仁存溝大橋高邊坡橋臺病害治理
    煤矸石浸水路堤不均勻沉降研究
    無處理軟基橋臺加固及驗算
    av视频在线观看入口| 真实男女啪啪啪动态图| 日本一二三区视频观看| 午夜亚洲福利在线播放| 我要搜黄色片| 69av精品久久久久久| 嫩草影院入口| 俺也久久电影网| 男插女下体视频免费在线播放| 69av精品久久久久久| 美女内射精品一级片tv| 国产精品久久久久久精品电影| 变态另类成人亚洲欧美熟女| 久久久a久久爽久久v久久| 麻豆av噜噜一区二区三区| 在线免费十八禁| 人人妻人人看人人澡| 九九爱精品视频在线观看| 国产麻豆成人av免费视频| 国产亚洲精品久久久com| 俺也久久电影网| 乱系列少妇在线播放| 国产精品美女特级片免费视频播放器| 亚洲人与动物交配视频| 一进一出抽搐动态| 老司机影院成人| 美女 人体艺术 gogo| 免费在线观看成人毛片| 免费看av在线观看网站| 亚洲人成网站在线观看播放| 久久精品人妻少妇| 一级毛片电影观看 | 大香蕉久久网| 欧美日韩精品成人综合77777| 最近视频中文字幕2019在线8| 欧美xxxx性猛交bbbb| 国产精品一区二区免费欧美| av.在线天堂| 高清日韩中文字幕在线| 最近手机中文字幕大全| 丝袜喷水一区| av.在线天堂| 国产精品免费一区二区三区在线| 嫩草影院入口| 中文字幕人妻熟人妻熟丝袜美| 国产黄片美女视频| 久久这里只有精品中国| 欧美成人一区二区免费高清观看| 亚洲国产精品国产精品| 白带黄色成豆腐渣| 日韩精品有码人妻一区| 女人十人毛片免费观看3o分钟| 久久精品国产99精品国产亚洲性色| av在线亚洲专区| 欧美+亚洲+日韩+国产| 丰满乱子伦码专区| 在线免费观看不下载黄p国产| 看片在线看免费视频| 日本免费一区二区三区高清不卡| 色哟哟·www| 精品人妻熟女av久视频| 在线播放无遮挡| 亚洲精品一区av在线观看| 噜噜噜噜噜久久久久久91| 热99re8久久精品国产| 日韩一区二区视频免费看| 亚洲成人久久性| 国产欧美日韩精品亚洲av| 久久久色成人| 免费一级毛片在线播放高清视频| 又粗又爽又猛毛片免费看| 色综合站精品国产| a级毛色黄片| 日韩成人伦理影院| 深夜a级毛片| 亚洲七黄色美女视频| 欧美成人精品欧美一级黄| 日韩国内少妇激情av| 搡老岳熟女国产| 狠狠狠狠99中文字幕| 人妻久久中文字幕网| 成人高潮视频无遮挡免费网站| av国产免费在线观看| 亚洲av熟女| 精品99又大又爽又粗少妇毛片| 欧美不卡视频在线免费观看| 国产蜜桃级精品一区二区三区| 丝袜美腿在线中文| 乱系列少妇在线播放| 一进一出抽搐动态| 乱系列少妇在线播放| 国产精品久久电影中文字幕| 少妇人妻精品综合一区二区 | 亚洲精品在线观看二区| 熟女电影av网| 欧美日韩一区二区视频在线观看视频在线 | 精品久久久久久久久av| av天堂在线播放| 国产午夜福利久久久久久| 亚洲色图av天堂| 联通29元200g的流量卡| 国产综合懂色| 一级av片app| 久久久久性生活片| 日韩高清综合在线| 97热精品久久久久久| 嫩草影院新地址| 色哟哟·www| 国产伦精品一区二区三区视频9| 久久精品人妻少妇| 久久久成人免费电影| 国产爱豆传媒在线观看| 老熟妇仑乱视频hdxx| 三级经典国产精品| 特级一级黄色大片| 亚洲性久久影院| 午夜视频国产福利| 麻豆一二三区av精品| 真实男女啪啪啪动态图| 乱码一卡2卡4卡精品| 国产又黄又爽又无遮挡在线| 蜜桃久久精品国产亚洲av| 成年女人看的毛片在线观看| 成人永久免费在线观看视频| 国产精品不卡视频一区二区| 亚洲中文日韩欧美视频| 精品少妇黑人巨大在线播放 | 全区人妻精品视频| 看免费成人av毛片| 99久久久亚洲精品蜜臀av| 日本黄色视频三级网站网址| 亚洲欧美精品自产自拍| 99久久久亚洲精品蜜臀av| 亚洲人成网站在线播放欧美日韩| 久久精品91蜜桃| 三级男女做爰猛烈吃奶摸视频| 日本熟妇午夜| 在线播放无遮挡| a级一级毛片免费在线观看| 国产高清不卡午夜福利| 亚洲人成网站在线播放欧美日韩| 国产老妇女一区| 变态另类成人亚洲欧美熟女| 日本熟妇午夜| 精品一区二区三区人妻视频| 日韩一本色道免费dvd| 永久网站在线| 又粗又爽又猛毛片免费看| 免费观看在线日韩| 日本黄色视频三级网站网址| 在线观看66精品国产| 久久精品国产亚洲av涩爱 | 美女高潮的动态| 成人鲁丝片一二三区免费| 久久久久性生活片| 可以在线观看的亚洲视频| 乱系列少妇在线播放| 麻豆成人午夜福利视频| 国产精品嫩草影院av在线观看| 久久国内精品自在自线图片| 欧美成人一区二区免费高清观看| 偷拍熟女少妇极品色| 国产高清有码在线观看视频| 男女之事视频高清在线观看| 久久久精品欧美日韩精品| av在线播放精品| 国产高清有码在线观看视频| 嫩草影院新地址| 一区二区三区四区激情视频 | av在线天堂中文字幕| 天美传媒精品一区二区| 春色校园在线视频观看| 97超级碰碰碰精品色视频在线观看| 国产精品一区二区免费欧美| 亚洲经典国产精华液单| 国产一区二区在线观看日韩| 免费人成视频x8x8入口观看| 极品教师在线视频| 少妇的逼水好多| 精品一区二区三区人妻视频| 香蕉av资源在线| 99久久无色码亚洲精品果冻| 非洲黑人性xxxx精品又粗又长| 日韩大尺度精品在线看网址| 国产免费一级a男人的天堂| 日日撸夜夜添| 97超级碰碰碰精品色视频在线观看| 婷婷精品国产亚洲av| 国产一区二区在线av高清观看| 特大巨黑吊av在线直播| 麻豆久久精品国产亚洲av| 18禁在线播放成人免费| 国产视频内射| 搞女人的毛片| 欧美精品国产亚洲| 亚洲av免费高清在线观看| 国产探花在线观看一区二区| 亚洲精品成人久久久久久| 日本一二三区视频观看| 搞女人的毛片| 变态另类丝袜制服| 精品日产1卡2卡| 亚洲成人久久性| 两个人视频免费观看高清| 成人综合一区亚洲| 色噜噜av男人的天堂激情| 国产片特级美女逼逼视频| 不卡视频在线观看欧美| 亚洲精品粉嫩美女一区| 亚洲图色成人| 国产男靠女视频免费网站| 亚洲婷婷狠狠爱综合网| 欧美高清成人免费视频www| 久久久久久伊人网av| 麻豆精品久久久久久蜜桃| 亚洲专区国产一区二区| 最近视频中文字幕2019在线8| 久久久久精品国产欧美久久久| 男人舔奶头视频| videossex国产| 久久久久久久久中文| 免费av毛片视频| 一级黄片播放器| 小说图片视频综合网站| 综合色丁香网| 免费一级毛片在线播放高清视频| 丰满的人妻完整版| 国产亚洲欧美98| 国产精品久久久久久久电影| 岛国在线免费视频观看| 亚洲人成网站在线观看播放| 免费不卡的大黄色大毛片视频在线观看 | 又黄又爽又免费观看的视频| 免费不卡的大黄色大毛片视频在线观看 | 亚洲欧美清纯卡通| 91久久精品国产一区二区成人| 国产单亲对白刺激| 婷婷精品国产亚洲av| 在线观看一区二区三区| 男女视频在线观看网站免费| 国产大屁股一区二区在线视频| 午夜福利在线观看免费完整高清在 | 欧美三级亚洲精品| 国产蜜桃级精品一区二区三区| 男人的好看免费观看在线视频| 在线播放国产精品三级| 精品不卡国产一区二区三区| 国产精品免费一区二区三区在线| 精品久久久久久久末码| 精品久久久久久久久久久久久| 亚洲精品一卡2卡三卡4卡5卡| av国产免费在线观看| 久久韩国三级中文字幕| 国产91av在线免费观看| 免费无遮挡裸体视频| 热99在线观看视频| 欧美xxxx黑人xx丫x性爽| 久久精品国产亚洲av天美| 国产白丝娇喘喷水9色精品| eeuss影院久久| videossex国产| 国内精品久久久久精免费| 国产乱人视频| 我要看日韩黄色一级片| 老司机福利观看| 日韩av不卡免费在线播放| 插阴视频在线观看视频| 国产精品一区二区三区四区免费观看 | 亚洲精品乱码久久久v下载方式| 日韩精品青青久久久久久| 春色校园在线视频观看| 看十八女毛片水多多多| 精品免费久久久久久久清纯| 欧美日本视频| 国产成人福利小说| 国产 一区 欧美 日韩| 亚洲精品成人久久久久久| 最近的中文字幕免费完整| 成年女人看的毛片在线观看| 97热精品久久久久久| 日本熟妇午夜| 久久综合国产亚洲精品| 亚洲人成网站在线观看播放| 床上黄色一级片| 国产乱人偷精品视频| 无遮挡黄片免费观看| 禁无遮挡网站| 女同久久另类99精品国产91| 69av精品久久久久久| 亚洲无线观看免费| 日韩强制内射视频| 最近最新中文字幕大全电影3| 搡老岳熟女国产| 99热全是精品| 国产三级中文精品| 别揉我奶头~嗯~啊~动态视频| 国产视频一区二区在线看| 最近中文字幕高清免费大全6| 国产成人a区在线观看| 国产精品av视频在线免费观看| 又爽又黄无遮挡网站| 日韩精品青青久久久久久| 欧美色视频一区免费| 免费av不卡在线播放| 欧美日韩乱码在线| 国产精品亚洲美女久久久| 有码 亚洲区| 精品久久久久久久久久免费视频| 亚洲综合色惰| 亚洲不卡免费看| 午夜福利18| 少妇人妻精品综合一区二区 | 日韩欧美精品v在线| 美女 人体艺术 gogo| 12—13女人毛片做爰片一| 亚洲欧美清纯卡通| 高清午夜精品一区二区三区 | 日韩欧美国产在线观看| 九九久久精品国产亚洲av麻豆| 女的被弄到高潮叫床怎么办| 18禁黄网站禁片免费观看直播| 久久久久性生活片| 婷婷精品国产亚洲av在线| 人妻少妇偷人精品九色| 国产爱豆传媒在线观看| 乱系列少妇在线播放| 国内精品美女久久久久久| 日韩高清综合在线| 插逼视频在线观看| 亚洲精品国产av成人精品 | 在线免费十八禁| 精品一区二区三区视频在线观看免费| 俺也久久电影网| av天堂中文字幕网| 午夜老司机福利剧场| 黄色日韩在线| 欧美xxxx黑人xx丫x性爽| 欧美不卡视频在线免费观看| 久久韩国三级中文字幕| 婷婷精品国产亚洲av| 欧美激情在线99| 欧美绝顶高潮抽搐喷水| 亚洲七黄色美女视频| 搞女人的毛片| 岛国在线免费视频观看| 色综合色国产| 99国产精品一区二区蜜桃av| 国产黄片美女视频| 97热精品久久久久久| av福利片在线观看| 亚洲va在线va天堂va国产| 黄色日韩在线| 18禁裸乳无遮挡免费网站照片| 一级黄色大片毛片| 亚洲精品国产av成人精品 | 99热全是精品| 亚洲美女黄片视频| 麻豆av噜噜一区二区三区| 蜜臀久久99精品久久宅男| 午夜激情福利司机影院| 欧美日本亚洲视频在线播放| 亚洲精品456在线播放app| 色在线成人网| 蜜桃亚洲精品一区二区三区| 久久国产乱子免费精品| 日本黄色片子视频| 成人毛片a级毛片在线播放| 亚洲欧美日韩卡通动漫| 淫妇啪啪啪对白视频| 此物有八面人人有两片| 最近中文字幕高清免费大全6| 精品一区二区三区av网在线观看| 国产精品福利在线免费观看| 欧美成人一区二区免费高清观看| 大又大粗又爽又黄少妇毛片口| 婷婷精品国产亚洲av在线| 亚洲欧美清纯卡通| 亚洲,欧美,日韩| 亚洲国产欧洲综合997久久,| 校园春色视频在线观看| 嫩草影视91久久| 中国美女看黄片| 色5月婷婷丁香| 亚洲自偷自拍三级| 成熟少妇高潮喷水视频| 22中文网久久字幕| 久久久久国内视频| 亚洲,欧美,日韩| 变态另类丝袜制服| 国产久久久一区二区三区| 久久人人精品亚洲av| 亚洲国产精品成人久久小说 | 男女边吃奶边做爰视频| 欧美日本视频| 岛国在线免费视频观看| 国产高清不卡午夜福利| 国产精品国产三级国产av玫瑰| 亚洲国产精品国产精品| 国产高清视频在线观看网站| 小蜜桃在线观看免费完整版高清| 一级毛片久久久久久久久女| 亚洲av五月六月丁香网| 18禁在线无遮挡免费观看视频 | 国产麻豆成人av免费视频| 久久久午夜欧美精品| 成年免费大片在线观看| 少妇被粗大猛烈的视频| av黄色大香蕉| 亚洲国产欧洲综合997久久,| 午夜免费男女啪啪视频观看 | av天堂中文字幕网| 一本一本综合久久| 91午夜精品亚洲一区二区三区| 麻豆乱淫一区二区| eeuss影院久久| 国产精品永久免费网站| 波多野结衣高清作品| 国产探花在线观看一区二区| 噜噜噜噜噜久久久久久91| 欧美+日韩+精品| 亚洲aⅴ乱码一区二区在线播放| av福利片在线观看| 1024手机看黄色片| 男人舔奶头视频| 免费看av在线观看网站| 久久精品国产亚洲av香蕉五月| 国产综合懂色| 国产精品日韩av在线免费观看| 国产精品不卡视频一区二区| 丝袜美腿在线中文| 成年女人看的毛片在线观看| 在线免费观看不下载黄p国产| 午夜福利18| 一卡2卡三卡四卡精品乱码亚洲| 久久久成人免费电影| 成人美女网站在线观看视频| 噜噜噜噜噜久久久久久91| 精品少妇黑人巨大在线播放 | 日本免费a在线| 网址你懂的国产日韩在线| 午夜老司机福利剧场| 三级国产精品欧美在线观看| 欧美色欧美亚洲另类二区| 日韩欧美精品v在线| 九九爱精品视频在线观看| 在线免费观看的www视频| 中文字幕久久专区| 国产精品女同一区二区软件| 亚洲专区国产一区二区| 别揉我奶头~嗯~啊~动态视频| 亚洲精品一区av在线观看| 成人亚洲精品av一区二区| 日韩欧美一区二区三区在线观看| 国产精品电影一区二区三区| 日本免费a在线| 波野结衣二区三区在线| 成人av在线播放网站| 国产男靠女视频免费网站| 亚洲av电影不卡..在线观看| 我要看日韩黄色一级片| 国产美女午夜福利| 日本黄大片高清| 俺也久久电影网| 一级黄色大片毛片| 久久韩国三级中文字幕| 熟妇人妻久久中文字幕3abv| 欧美性猛交╳xxx乱大交人| 赤兔流量卡办理| 日本精品一区二区三区蜜桃| 国内久久婷婷六月综合欲色啪| 日本色播在线视频| 夜夜夜夜夜久久久久| 日日啪夜夜撸| 国产成人aa在线观看| 国产av在哪里看| 九九爱精品视频在线观看| 亚洲国产欧洲综合997久久,| 久久久久免费精品人妻一区二区| 免费人成视频x8x8入口观看| 三级国产精品欧美在线观看| 美女高潮的动态| 菩萨蛮人人尽说江南好唐韦庄 | 国产人妻一区二区三区在| 亚洲美女黄片视频| 又爽又黄a免费视频| 哪里可以看免费的av片| 欧美绝顶高潮抽搐喷水| 午夜免费激情av| 国产乱人视频| 午夜精品国产一区二区电影 | 一个人免费在线观看电影| 黄色日韩在线| 熟女人妻精品中文字幕| 久久精品国产自在天天线| 国产成人a区在线观看| 最近视频中文字幕2019在线8| 亚洲经典国产精华液单| 精品人妻熟女av久视频| 欧美bdsm另类| 久久久精品大字幕| 午夜亚洲福利在线播放| 亚洲国产高清在线一区二区三| 免费在线观看成人毛片| 久久99热这里只有精品18| 国产av麻豆久久久久久久| 91在线精品国自产拍蜜月| 精品一区二区三区视频在线观看免费| 日本黄色片子视频| 国产精品,欧美在线| 国产亚洲精品综合一区在线观看| 天天躁日日操中文字幕| 午夜福利高清视频| 欧美成人精品欧美一级黄| 亚洲精品色激情综合| 国产成人一区二区在线| 亚洲va在线va天堂va国产| 99国产极品粉嫩在线观看| 直男gayav资源| 亚洲无线在线观看| 免费av不卡在线播放| 狂野欧美白嫩少妇大欣赏| 精品国产三级普通话版| 真人做人爱边吃奶动态| 久久韩国三级中文字幕| 中文字幕av成人在线电影| 成年av动漫网址| 久久精品影院6| 精品少妇黑人巨大在线播放 | 一区二区三区免费毛片| 久久欧美精品欧美久久欧美| 亚洲国产精品国产精品| 久久久精品欧美日韩精品| 国产高清不卡午夜福利| 欧美色欧美亚洲另类二区| 日日撸夜夜添| 国产成人一区二区在线| 亚洲av免费高清在线观看| 美女大奶头视频| 日本一本二区三区精品| 午夜亚洲福利在线播放| 精品一区二区三区人妻视频| 女的被弄到高潮叫床怎么办| 亚洲精品粉嫩美女一区| 亚洲色图av天堂| 51国产日韩欧美| 热99在线观看视频| 日韩欧美在线乱码| 午夜免费激情av| 亚洲欧美日韩无卡精品| 亚洲在线观看片| 一区二区三区四区激情视频 | 91久久精品国产一区二区三区| 网址你懂的国产日韩在线| 三级毛片av免费| 日韩中字成人| 级片在线观看| 波多野结衣高清作品| 国产黄色视频一区二区在线观看 | 国产精品三级大全| 中国美白少妇内射xxxbb| 日日摸夜夜添夜夜爱| 免费观看在线日韩| 给我免费播放毛片高清在线观看| 我的老师免费观看完整版| 乱系列少妇在线播放| 国产精品无大码| 又爽又黄a免费视频| 女生性感内裤真人,穿戴方法视频| 欧美日本视频| 欧美成人一区二区免费高清观看| 国产高潮美女av| 久久99热6这里只有精品| 菩萨蛮人人尽说江南好唐韦庄 | 精品午夜福利在线看| av专区在线播放| 免费看美女性在线毛片视频| 久久精品国产亚洲av涩爱 | 久久综合国产亚洲精品| 日韩在线高清观看一区二区三区| 国产高清三级在线| 久久精品国产亚洲av香蕉五月| 国产av一区在线观看免费| 亚洲一级一片aⅴ在线观看| 国产精品电影一区二区三区| 亚洲最大成人av| 搞女人的毛片| 亚洲高清免费不卡视频| 精品不卡国产一区二区三区| 免费看av在线观看网站| 亚洲av成人av| 一个人看的www免费观看视频| aaaaa片日本免费| 黑人高潮一二区| 免费无遮挡裸体视频| 一级黄片播放器| 亚洲18禁久久av| h日本视频在线播放| 桃色一区二区三区在线观看| 久久精品影院6| 在线观看66精品国产| 欧美中文日本在线观看视频| 波多野结衣高清作品| 网址你懂的国产日韩在线| 最新在线观看一区二区三区| 国产黄a三级三级三级人| 波野结衣二区三区在线| 干丝袜人妻中文字幕| 精品久久久久久久久av| 99视频精品全部免费 在线| 女的被弄到高潮叫床怎么办| 国产精品爽爽va在线观看网站| 亚洲中文字幕一区二区三区有码在线看| 成年av动漫网址| 免费大片18禁| 中文字幕av成人在线电影|