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

    基于ABAQUS軟件的覆蓋層中塑性混凝土防滲墻的應(yīng)力變形分析

    2013-12-23 17:09:59李龍舟
    科技創(chuàng)新與應(yīng)用 2014年2期
    關(guān)鍵詞:土石壩覆蓋層防滲墻

    李龍舟

    摘 要:塑性混凝土防滲墻在土石壩加固工程中得到了廣泛的應(yīng)用,防滲墻的設(shè)計(jì)和施工對(duì)土石壩的安全運(yùn)行至關(guān)重要。文章基于有限元分析軟件ABAQUS,實(shí)現(xiàn)Goodman非線性接觸模型的二次開(kāi)發(fā),并應(yīng)用于土體與結(jié)構(gòu)材料接觸面的數(shù)值模擬。結(jié)合某土石壩工程實(shí)例,對(duì)軟弱覆蓋層上塑性混凝土防滲墻進(jìn)行數(shù)值計(jì)算,分析防滲墻在蓄水前和蓄水后等不同階段的工作特性,同時(shí)考慮混凝土彈模的改變和泥皮參的變化對(duì)防滲墻應(yīng)力變形的影響,為塑性混凝土防滲墻在土石壩中的應(yīng)用提供一定參考依據(jù)。

    關(guān)鍵詞:土石壩;覆蓋層;塑性混凝土;防滲墻;接觸模型;應(yīng)力變形

    引言

    在深厚覆蓋層上修建土石壩,常采用澆筑混凝土防滲墻、帷幕灌漿等技術(shù),來(lái)截?cái)嗟叵聺B流通道,保證壩基的抗滑穩(wěn)定性[1]。大壩壩體在施工和蓄水以后,由于受到大壩自重和水荷載的作用,大壩和壩基會(huì)產(chǎn)生較大的變形?;炷练罎B墻作為大壩防滲體,由于其彈性模量與周?chē)馏w的巨大差異,變形的不協(xié)調(diào)性,會(huì)使防滲墻承受很高的拉應(yīng)力和壓應(yīng)力而致使混凝土防滲墻遭受破壞。20世紀(jì)60年代國(guó)外通過(guò)對(duì)常規(guī)混凝土的改性,摻加一定數(shù)量的膨潤(rùn)土或黏土來(lái)替代常規(guī)混凝土中的大部分水泥,而出現(xiàn)了塑性混凝土材料(其彈模一般小于2000MP,抗壓強(qiáng)度低于5MP),并開(kāi)始應(yīng)用于土石壩壩體的防滲處理,引起了國(guó)際大壩工程界的高度關(guān)注。20世紀(jì)80年代中后期,我國(guó)開(kāi)始研究塑性混凝土,并在一些臨時(shí)圍堰中得以應(yīng)用,同時(shí)少量試用于大壩[2]。

    在防滲墻的施工過(guò)程中,墻體混凝土在開(kāi)挖成槽時(shí),多采用泥漿護(hù)壁的方法來(lái)保證槽孔穩(wěn)定,泥漿滲入地層會(huì)在槽孔孔壁表面形成泥皮,在混凝土澆筑完成后,泥皮依然會(huì)粘結(jié)在槽壁上。因此,墻體上下游面與覆蓋層土體之間都會(huì)有泥漿滲透形成的泥皮[2]。很多學(xué)者對(duì)土與結(jié)構(gòu)接觸面之間的力學(xué)特性及其影響因素進(jìn)行了研究,證實(shí)泥皮的存在會(huì)對(duì)接觸面的力學(xué)性質(zhì)產(chǎn)生顯著影響。張 嘎,張建民進(jìn)行了粗粒土與結(jié)構(gòu)之間夾有泥皮和無(wú)泥皮兩種狀態(tài)的接觸面單調(diào)和循環(huán)剪切試驗(yàn),表明泥皮對(duì)接觸面力學(xué)特性的影響程度不僅與泥皮本身的特性及其厚度有關(guān),還與構(gòu)成接觸面的結(jié)構(gòu)面材料與粗粒土的特性有關(guān)[3]。張治軍、饒錫保、王志軍等采用大型疊環(huán)單剪儀,對(duì)砂礫石料與夾泥皮的結(jié)構(gòu)物所形成接觸面的力學(xué)特性進(jìn)行了相關(guān)試驗(yàn)分析,分析了泥皮厚度對(duì)接觸面特性的影響并提出了接觸面在不同泥皮厚度下所適用的數(shù)學(xué)模型[4]。楊春鳴、邵生俊、劉鑫通過(guò)實(shí)驗(yàn)確定了夾有泥皮時(shí)不同粗粒土級(jí)配與結(jié)構(gòu)接觸面的強(qiáng)度和摩擦系數(shù)[5]。

    在塑性混凝土防滲墻數(shù)值模擬中往往忽略了泥皮對(duì)防滲墻應(yīng)力變形的影響。本文利用有限元分析軟件ABAQUS的用戶(hù)子程序,實(shí)現(xiàn)Goodman模型的二次開(kāi)發(fā),對(duì)在覆蓋層中采用塑性混凝土防滲墻防滲的某土石壩進(jìn)行了數(shù)值計(jì)算,分析防滲墻在蓄水前和蓄水后等不同階段的工作狀態(tài),同時(shí)考慮塑性混凝土彈模的改變和泥皮參的變化對(duì)防滲墻應(yīng)力變形的影響。

    1 計(jì)算原理

    1.1 Goodman模型

    本文利用Goodman接觸模型來(lái)模擬塑性混凝土防滲墻和周?chē)馏w之間的接觸作用,能夠更好地考慮墻體和土體之間的承載力分配問(wèn)題。

    Goodman單元不考慮兩個(gè)方向的相互耦合,接觸面是在受力之間完全吻合,即單元沒(méi)有厚度只有長(zhǎng)度,法向應(yīng)力只與法向相對(duì)位移有關(guān),剪應(yīng)力只與切向相對(duì)位移有關(guān)[6]。本構(gòu)關(guān)系為

    (1)

    式中: ,ks,kn分別為切向剛度系數(shù)和法向剛度系數(shù)。

    1.2 Goodman模型在ABAQUS中的二次開(kāi)發(fā)

    ABAQUS是由美國(guó)HKS公司開(kāi)發(fā)的非線性有限元分析軟件系統(tǒng)。它在材料、幾何和接觸非線性方面的分析能力居世界領(lǐng)先水平,以高求解效率和高計(jì)算精度在工程界和學(xué)術(shù)界贏得了聲譽(yù)和信賴(lài),是強(qiáng)大的非線性有限元計(jì)算分析工具。ABAQUS為廣大的用戶(hù)提供了大量的單元庫(kù)和求解模型,用戶(hù)可以利用這些模型處理大多數(shù)的問(wèn)題。并能利用其提供的用戶(hù)子程序接口,增加單元類(lèi)型、本構(gòu)關(guān)系、接觸模型等。在ABAQUS所提供的通用平臺(tái)上有可能實(shí)現(xiàn)特定問(wèn)題的有限元數(shù)值仿真計(jì)算,以滿(mǎn)足用戶(hù)的個(gè)性化分析計(jì)算的需求[7]。

    ABAQUS提供用戶(hù)自定義接觸面摩擦模型的子程序是FRIC。ABAQUS的提供的面與面接觸和Goodman模型相同的是沒(méi)有厚度只有長(zhǎng)度,區(qū)別僅在于摩擦模型不同。因此,在FRIC子程序中只要考慮接觸面上摩擦接觸特性就可以,而把判斷接觸面是否脫開(kāi)的任務(wù)交給ABAQUS進(jìn)行。在FRIC程序中,我們只需給出兩個(gè)方向的ks??死蚝袜嚳险J(rèn)為剪應(yīng)力和相對(duì)剪切位移之間符合雙曲線關(guān)系[8],則ks1、ks2表示為

    (2)

    (3)

    式中K1、K2、Rf、n為非線性指標(biāo),由試驗(yàn)確定;?啄為接觸面摩擦角。

    2 工程算例

    2.1 工程概況

    某水庫(kù)土石壩加固對(duì)舊壩體與覆蓋層采用塑性混凝土防滲墻、強(qiáng)風(fēng)化基巖層采用帷幕灌漿的綜合防滲措施進(jìn)行處理。壩頂高程1136.90m,壩高36.7m,壩頂寬度6m。防滲墻樁號(hào)范圍為:0-070~0+936.8,全長(zhǎng)1006.8m,設(shè)計(jì)墻厚0.8m,最大墻深59.92m。大壩加固后樁號(hào)0+562.5斷面圖如圖1所示,壩體填土為中等壓縮低液限黏土,壩基從上到下依次為覆蓋層、強(qiáng)風(fēng)化巖、若風(fēng)化巖。防滲墻厚布置在壩軸線上有1.8m處,底部深入弱風(fēng)化巖5m。

    圖1 0+562.5斷面大壩剖面圖

    2.2 模型概況

    模型計(jì)算按平面應(yīng)變問(wèn)題處理,計(jì)算整個(gè)范圍自壩址和壩踵分別向上游延伸100m,模型中壩基兩端采用滑動(dòng)支座,約束水平方向位移,壩基底端采用固定支座,約束一切變形。加載過(guò)程按土石壩施工和蓄水過(guò)程模擬。壩體、覆蓋層和塑性混凝土防滲墻材料均按照非線性材料考慮,計(jì)算模型采用鄧肯·張E-B模型。土與結(jié)構(gòu)的接觸面采用基于ABAQUS二次開(kāi)發(fā)的Goodman單元模擬。計(jì)算材料參數(shù)見(jiàn)表1。由于防滲墻頂端的高塑性土和防滲墻底部的灌漿帷幕與壩址處的排水棱體的面積相對(duì)較小,且缺乏相對(duì)可靠的土工參數(shù),故計(jì)算中忽略其影響。

    2.3 計(jì)算結(jié)果

    圖2~圖5給出了墻體(K=3000)在蓄水前和蓄水后的應(yīng)力變形沿高程的分布規(guī)律。

    (1)蓄水前墻體的應(yīng)力變形分析。墻體豎向應(yīng)力均為壓應(yīng)力,墻體上下游面應(yīng)力分布情況基本一致,對(duì)于墻體頂部和下部,都出現(xiàn)了上游面豎向應(yīng)力略大于下游面,對(duì)于墻體中部,上下游面豎向應(yīng)力差別不大,受力較均勻。墻體的大主應(yīng)力的分布規(guī)律與豎向應(yīng)力基本相同,小主應(yīng)力隨墻體高程的增加而逐漸減小,上下游面小主應(yīng)力差別不大,受力較均勻。

    墻體有向上游傾斜的趨勢(shì),最大水平變形為23.29cm,發(fā)生在墻體頂處;墻體最大垂直變形為37.69cm,豎直向下,同樣發(fā)生在墻體頂部。

    圖2 防滲墻的水平位移和豎向位移沿高程分布圖

    圖3 墻體大主應(yīng)力沿高程分布圖

    (a)蓄水前;(b)蓄水后

    圖4 墻體小主應(yīng)力沿高程分布圖

    (a)蓄水前;(b)蓄水后

    (2)蓄水后墻體的應(yīng)力變形分析。蓄水到正常蓄水位,墻體的中下部在上游面受到較大的水壓力,可以看到墻體的水平變形向下游有一定的回復(fù),但仍然偏向上游,這就使得墻體的豎向應(yīng)力或者大主應(yīng)力有所減小。同時(shí),在水壓力的作用下,其圍壓提高,小主應(yīng)力明顯增大。從圖4、圖5可以看出,在墻體下部覆蓋層與基巖交界面處出現(xiàn)了一定程度的應(yīng)力集中,最大壓應(yīng)力達(dá)到了1.54MP;墻體頂部新填壩體和覆蓋層交界處也出現(xiàn)了一定程度的應(yīng)力集中。蓄水對(duì)墻體的垂直變形影響不大。

    (3)塑性混凝土模量變化對(duì)墻體的影響。為了考慮塑性混凝土模量的變化對(duì)墻體應(yīng)力變形的影響,取模量K為3000、5000、8000這三種方案進(jìn)行數(shù)值計(jì)算,分析比較應(yīng)力變形的變化。計(jì)算結(jié)果如表2所示。

    表2 蓄水后各方案墻體應(yīng)力變形分布特征值

    隨著塑性混凝土模量的增加,墻體變形減小,但變化不大;墻體的大主應(yīng)力和小主應(yīng)力整體呈增大趨勢(shì),最大壓應(yīng)力普遍增高,墻體的受拉程度明顯增加。在彈模達(dá)到8000時(shí),墻體下游面出現(xiàn)了受拉區(qū)域,小主應(yīng)力隨高程的分布出現(xiàn)了劇烈的振蕩,對(duì)防滲墻的受力狀態(tài)產(chǎn)生不利影響。模量較小對(duì)防滲墻的受力狀態(tài)不會(huì)產(chǎn)生較大的影響,而防滲墻的強(qiáng)度則隨著模量的降低顯著的減小;模量增大,防滲墻墻體強(qiáng)度大大的提高了,但防滲墻也出現(xiàn)了明顯的受拉區(qū),對(duì)防滲墻不利。因此,在設(shè)計(jì)時(shí)應(yīng)根據(jù)工程地質(zhì)條件、水位及壩型等綜合選取防滲墻的彈性模量,盡量降低塑性混凝土的模量,同時(shí)提高墻體的強(qiáng)度。

    (4)泥皮的參數(shù)變化對(duì)墻體的影響分析。為了考慮泥皮的作用對(duì)防滲墻應(yīng)力變形的影響,分別對(duì)接觸面摩擦系數(shù)為0.2,0.3和0.4的情況進(jìn)行模擬,對(duì)防滲墻在蓄水前和蓄水后的應(yīng)力變形情況進(jìn)行分析。結(jié)果見(jiàn)表3所示。從表中可以看出,無(wú)論是在蓄水前還是蓄水后,摩擦系數(shù)越大,防滲墻的大主應(yīng)力、小主應(yīng)力和豎向應(yīng)力都會(huì)增加,而墻體的最大水平位移減小。摩擦系數(shù)越大,墻體和周?chē)馏w之間的變形協(xié)調(diào)能力略微下降。接觸面參數(shù)的變化對(duì)防滲墻的應(yīng)力變形影響是有規(guī)律的,但變化量相對(duì)較小。

    表2 墻體在各個(gè)階段的最大應(yīng)力和位移

    3 結(jié)束語(yǔ)

    本文基于有限元分析軟件ABAQUS,實(shí)現(xiàn)Goodman模型的二次開(kāi)發(fā),對(duì)覆蓋層中采用塑性混凝土防滲墻的某土石壩進(jìn)行了數(shù)值計(jì)算,同時(shí)分析比較了彈性模量和接觸面參數(shù)對(duì)壩體塑性混凝土防滲墻的應(yīng)力變形影響。通過(guò)數(shù)值計(jì)算分析,可得到如下結(jié)論:

    (1)防滲墻在蓄水前和蓄水后均處于受壓狀態(tài),無(wú)拉應(yīng)力。在基巖到覆蓋層和覆蓋層到填土的過(guò)渡段都出現(xiàn)了一定程度的應(yīng)力集中,最大壓應(yīng)力不超過(guò)2.0MP,所以從抗壓強(qiáng)度上考慮,墻體是安全可靠的。

    (2)墻體的最大應(yīng)力值受混凝土模量的影響較為敏感。隨著混凝土模量的提高,塑性混凝土的強(qiáng)度越高,同時(shí)混凝土的最大應(yīng)力值也越高。但墻體的變形對(duì)混凝土模量的影響不敏感。盡量降低塑性混凝土的模量,對(duì)降低墻體的應(yīng)力是有好處的。

    (3)防滲墻的墻體與覆蓋土體接觸參數(shù)的變化會(huì)對(duì)墻體應(yīng)力變形產(chǎn)生影響較小。

    參考文獻(xiàn)

    [1]高鐘璞.大壩基礎(chǔ)防滲墻[M].北京:中國(guó)電力出版社,2000.

    [2]王清友,孫萬(wàn)功,熊歡.塑性混凝土防滲墻[M].北京:中國(guó)水利水電出版社,2008.

    [3]張嘎,張建民.夾有泥皮粗粒土與結(jié)構(gòu)接觸面力學(xué)特性試驗(yàn)研究[J].巖土力學(xué),2005,26(9):1374-1378.

    [4]張治軍,饒錫保,王志軍.泥皮厚度對(duì)結(jié)構(gòu)接觸面力學(xué)特性影響的試驗(yàn)研究[J].巖土力學(xué),2008,29(9):2433-2438.

    [5]楊春鳴,邵生俊,劉鑫.考慮泥皮作用的深厚覆蓋層壩基防滲墻應(yīng)力變形特性分析[J].西北農(nóng)林科技大學(xué)學(xué)報(bào),2013,41(8):1~10.

    [6]Goodman R F,Taylor R L,Brekke T L. A Model for the Mechanics of Jointed Rock[J]. Journal of the Soil Mechanics and Foundation Div. ASCE,1968,94(3):637-660.

    [7]費(fèi)康,張建偉.ABAQUS在巖土工程中的應(yīng)用[M].北京:中國(guó)水利水電出版社,2010.

    [8]Clough,G.W.Duncan,J. M. Finite Element Analysis of Retaining Wall Behavior[J]. Journal of the Soil Mechanics and Foundation Div. ,ASCE,1971,97(12):1657-1674.

    2.3 計(jì)算結(jié)果

    圖2~圖5給出了墻體(K=3000)在蓄水前和蓄水后的應(yīng)力變形沿高程的分布規(guī)律。

    (1)蓄水前墻體的應(yīng)力變形分析。墻體豎向應(yīng)力均為壓應(yīng)力,墻體上下游面應(yīng)力分布情況基本一致,對(duì)于墻體頂部和下部,都出現(xiàn)了上游面豎向應(yīng)力略大于下游面,對(duì)于墻體中部,上下游面豎向應(yīng)力差別不大,受力較均勻。墻體的大主應(yīng)力的分布規(guī)律與豎向應(yīng)力基本相同,小主應(yīng)力隨墻體高程的增加而逐漸減小,上下游面小主應(yīng)力差別不大,受力較均勻。

    墻體有向上游傾斜的趨勢(shì),最大水平變形為23.29cm,發(fā)生在墻體頂處;墻體最大垂直變形為37.69cm,豎直向下,同樣發(fā)生在墻體頂部。

    圖2 防滲墻的水平位移和豎向位移沿高程分布圖

    圖3 墻體大主應(yīng)力沿高程分布圖

    (a)蓄水前;(b)蓄水后

    圖4 墻體小主應(yīng)力沿高程分布圖

    (a)蓄水前;(b)蓄水后

    (2)蓄水后墻體的應(yīng)力變形分析。蓄水到正常蓄水位,墻體的中下部在上游面受到較大的水壓力,可以看到墻體的水平變形向下游有一定的回復(fù),但仍然偏向上游,這就使得墻體的豎向應(yīng)力或者大主應(yīng)力有所減小。同時(shí),在水壓力的作用下,其圍壓提高,小主應(yīng)力明顯增大。從圖4、圖5可以看出,在墻體下部覆蓋層與基巖交界面處出現(xiàn)了一定程度的應(yīng)力集中,最大壓應(yīng)力達(dá)到了1.54MP;墻體頂部新填壩體和覆蓋層交界處也出現(xiàn)了一定程度的應(yīng)力集中。蓄水對(duì)墻體的垂直變形影響不大。

    (3)塑性混凝土模量變化對(duì)墻體的影響。為了考慮塑性混凝土模量的變化對(duì)墻體應(yīng)力變形的影響,取模量K為3000、5000、8000這三種方案進(jìn)行數(shù)值計(jì)算,分析比較應(yīng)力變形的變化。計(jì)算結(jié)果如表2所示。

    表2 蓄水后各方案墻體應(yīng)力變形分布特征值

    隨著塑性混凝土模量的增加,墻體變形減小,但變化不大;墻體的大主應(yīng)力和小主應(yīng)力整體呈增大趨勢(shì),最大壓應(yīng)力普遍增高,墻體的受拉程度明顯增加。在彈模達(dá)到8000時(shí),墻體下游面出現(xiàn)了受拉區(qū)域,小主應(yīng)力隨高程的分布出現(xiàn)了劇烈的振蕩,對(duì)防滲墻的受力狀態(tài)產(chǎn)生不利影響。模量較小對(duì)防滲墻的受力狀態(tài)不會(huì)產(chǎn)生較大的影響,而防滲墻的強(qiáng)度則隨著模量的降低顯著的減??;模量增大,防滲墻墻體強(qiáng)度大大的提高了,但防滲墻也出現(xiàn)了明顯的受拉區(qū),對(duì)防滲墻不利。因此,在設(shè)計(jì)時(shí)應(yīng)根據(jù)工程地質(zhì)條件、水位及壩型等綜合選取防滲墻的彈性模量,盡量降低塑性混凝土的模量,同時(shí)提高墻體的強(qiáng)度。

    (4)泥皮的參數(shù)變化對(duì)墻體的影響分析。為了考慮泥皮的作用對(duì)防滲墻應(yīng)力變形的影響,分別對(duì)接觸面摩擦系數(shù)為0.2,0.3和0.4的情況進(jìn)行模擬,對(duì)防滲墻在蓄水前和蓄水后的應(yīng)力變形情況進(jìn)行分析。結(jié)果見(jiàn)表3所示。從表中可以看出,無(wú)論是在蓄水前還是蓄水后,摩擦系數(shù)越大,防滲墻的大主應(yīng)力、小主應(yīng)力和豎向應(yīng)力都會(huì)增加,而墻體的最大水平位移減小。摩擦系數(shù)越大,墻體和周?chē)馏w之間的變形協(xié)調(diào)能力略微下降。接觸面參數(shù)的變化對(duì)防滲墻的應(yīng)力變形影響是有規(guī)律的,但變化量相對(duì)較小。

    表2 墻體在各個(gè)階段的最大應(yīng)力和位移

    3 結(jié)束語(yǔ)

    本文基于有限元分析軟件ABAQUS,實(shí)現(xiàn)Goodman模型的二次開(kāi)發(fā),對(duì)覆蓋層中采用塑性混凝土防滲墻的某土石壩進(jìn)行了數(shù)值計(jì)算,同時(shí)分析比較了彈性模量和接觸面參數(shù)對(duì)壩體塑性混凝土防滲墻的應(yīng)力變形影響。通過(guò)數(shù)值計(jì)算分析,可得到如下結(jié)論:

    (1)防滲墻在蓄水前和蓄水后均處于受壓狀態(tài),無(wú)拉應(yīng)力。在基巖到覆蓋層和覆蓋層到填土的過(guò)渡段都出現(xiàn)了一定程度的應(yīng)力集中,最大壓應(yīng)力不超過(guò)2.0MP,所以從抗壓強(qiáng)度上考慮,墻體是安全可靠的。

    (2)墻體的最大應(yīng)力值受混凝土模量的影響較為敏感。隨著混凝土模量的提高,塑性混凝土的強(qiáng)度越高,同時(shí)混凝土的最大應(yīng)力值也越高。但墻體的變形對(duì)混凝土模量的影響不敏感。盡量降低塑性混凝土的模量,對(duì)降低墻體的應(yīng)力是有好處的。

    (3)防滲墻的墻體與覆蓋土體接觸參數(shù)的變化會(huì)對(duì)墻體應(yīng)力變形產(chǎn)生影響較小。

    參考文獻(xiàn)

    [1]高鐘璞.大壩基礎(chǔ)防滲墻[M].北京:中國(guó)電力出版社,2000.

    [2]王清友,孫萬(wàn)功,熊歡.塑性混凝土防滲墻[M].北京:中國(guó)水利水電出版社,2008.

    [3]張嘎,張建民.夾有泥皮粗粒土與結(jié)構(gòu)接觸面力學(xué)特性試驗(yàn)研究[J].巖土力學(xué),2005,26(9):1374-1378.

    [4]張治軍,饒錫保,王志軍.泥皮厚度對(duì)結(jié)構(gòu)接觸面力學(xué)特性影響的試驗(yàn)研究[J].巖土力學(xué),2008,29(9):2433-2438.

    [5]楊春鳴,邵生俊,劉鑫.考慮泥皮作用的深厚覆蓋層壩基防滲墻應(yīng)力變形特性分析[J].西北農(nóng)林科技大學(xué)學(xué)報(bào),2013,41(8):1~10.

    [6]Goodman R F,Taylor R L,Brekke T L. A Model for the Mechanics of Jointed Rock[J]. Journal of the Soil Mechanics and Foundation Div. ASCE,1968,94(3):637-660.

    [7]費(fèi)康,張建偉.ABAQUS在巖土工程中的應(yīng)用[M].北京:中國(guó)水利水電出版社,2010.

    [8]Clough,G.W.Duncan,J. M. Finite Element Analysis of Retaining Wall Behavior[J]. Journal of the Soil Mechanics and Foundation Div. ,ASCE,1971,97(12):1657-1674.

    2.3 計(jì)算結(jié)果

    圖2~圖5給出了墻體(K=3000)在蓄水前和蓄水后的應(yīng)力變形沿高程的分布規(guī)律。

    (1)蓄水前墻體的應(yīng)力變形分析。墻體豎向應(yīng)力均為壓應(yīng)力,墻體上下游面應(yīng)力分布情況基本一致,對(duì)于墻體頂部和下部,都出現(xiàn)了上游面豎向應(yīng)力略大于下游面,對(duì)于墻體中部,上下游面豎向應(yīng)力差別不大,受力較均勻。墻體的大主應(yīng)力的分布規(guī)律與豎向應(yīng)力基本相同,小主應(yīng)力隨墻體高程的增加而逐漸減小,上下游面小主應(yīng)力差別不大,受力較均勻。

    墻體有向上游傾斜的趨勢(shì),最大水平變形為23.29cm,發(fā)生在墻體頂處;墻體最大垂直變形為37.69cm,豎直向下,同樣發(fā)生在墻體頂部。

    圖2 防滲墻的水平位移和豎向位移沿高程分布圖

    圖3 墻體大主應(yīng)力沿高程分布圖

    (a)蓄水前;(b)蓄水后

    圖4 墻體小主應(yīng)力沿高程分布圖

    (a)蓄水前;(b)蓄水后

    (2)蓄水后墻體的應(yīng)力變形分析。蓄水到正常蓄水位,墻體的中下部在上游面受到較大的水壓力,可以看到墻體的水平變形向下游有一定的回復(fù),但仍然偏向上游,這就使得墻體的豎向應(yīng)力或者大主應(yīng)力有所減小。同時(shí),在水壓力的作用下,其圍壓提高,小主應(yīng)力明顯增大。從圖4、圖5可以看出,在墻體下部覆蓋層與基巖交界面處出現(xiàn)了一定程度的應(yīng)力集中,最大壓應(yīng)力達(dá)到了1.54MP;墻體頂部新填壩體和覆蓋層交界處也出現(xiàn)了一定程度的應(yīng)力集中。蓄水對(duì)墻體的垂直變形影響不大。

    (3)塑性混凝土模量變化對(duì)墻體的影響。為了考慮塑性混凝土模量的變化對(duì)墻體應(yīng)力變形的影響,取模量K為3000、5000、8000這三種方案進(jìn)行數(shù)值計(jì)算,分析比較應(yīng)力變形的變化。計(jì)算結(jié)果如表2所示。

    表2 蓄水后各方案墻體應(yīng)力變形分布特征值

    隨著塑性混凝土模量的增加,墻體變形減小,但變化不大;墻體的大主應(yīng)力和小主應(yīng)力整體呈增大趨勢(shì),最大壓應(yīng)力普遍增高,墻體的受拉程度明顯增加。在彈模達(dá)到8000時(shí),墻體下游面出現(xiàn)了受拉區(qū)域,小主應(yīng)力隨高程的分布出現(xiàn)了劇烈的振蕩,對(duì)防滲墻的受力狀態(tài)產(chǎn)生不利影響。模量較小對(duì)防滲墻的受力狀態(tài)不會(huì)產(chǎn)生較大的影響,而防滲墻的強(qiáng)度則隨著模量的降低顯著的減??;模量增大,防滲墻墻體強(qiáng)度大大的提高了,但防滲墻也出現(xiàn)了明顯的受拉區(qū),對(duì)防滲墻不利。因此,在設(shè)計(jì)時(shí)應(yīng)根據(jù)工程地質(zhì)條件、水位及壩型等綜合選取防滲墻的彈性模量,盡量降低塑性混凝土的模量,同時(shí)提高墻體的強(qiáng)度。

    (4)泥皮的參數(shù)變化對(duì)墻體的影響分析。為了考慮泥皮的作用對(duì)防滲墻應(yīng)力變形的影響,分別對(duì)接觸面摩擦系數(shù)為0.2,0.3和0.4的情況進(jìn)行模擬,對(duì)防滲墻在蓄水前和蓄水后的應(yīng)力變形情況進(jìn)行分析。結(jié)果見(jiàn)表3所示。從表中可以看出,無(wú)論是在蓄水前還是蓄水后,摩擦系數(shù)越大,防滲墻的大主應(yīng)力、小主應(yīng)力和豎向應(yīng)力都會(huì)增加,而墻體的最大水平位移減小。摩擦系數(shù)越大,墻體和周?chē)馏w之間的變形協(xié)調(diào)能力略微下降。接觸面參數(shù)的變化對(duì)防滲墻的應(yīng)力變形影響是有規(guī)律的,但變化量相對(duì)較小。

    表2 墻體在各個(gè)階段的最大應(yīng)力和位移

    3 結(jié)束語(yǔ)

    本文基于有限元分析軟件ABAQUS,實(shí)現(xiàn)Goodman模型的二次開(kāi)發(fā),對(duì)覆蓋層中采用塑性混凝土防滲墻的某土石壩進(jìn)行了數(shù)值計(jì)算,同時(shí)分析比較了彈性模量和接觸面參數(shù)對(duì)壩體塑性混凝土防滲墻的應(yīng)力變形影響。通過(guò)數(shù)值計(jì)算分析,可得到如下結(jié)論:

    (1)防滲墻在蓄水前和蓄水后均處于受壓狀態(tài),無(wú)拉應(yīng)力。在基巖到覆蓋層和覆蓋層到填土的過(guò)渡段都出現(xiàn)了一定程度的應(yīng)力集中,最大壓應(yīng)力不超過(guò)2.0MP,所以從抗壓強(qiáng)度上考慮,墻體是安全可靠的。

    (2)墻體的最大應(yīng)力值受混凝土模量的影響較為敏感。隨著混凝土模量的提高,塑性混凝土的強(qiáng)度越高,同時(shí)混凝土的最大應(yīng)力值也越高。但墻體的變形對(duì)混凝土模量的影響不敏感。盡量降低塑性混凝土的模量,對(duì)降低墻體的應(yīng)力是有好處的。

    (3)防滲墻的墻體與覆蓋土體接觸參數(shù)的變化會(huì)對(duì)墻體應(yīng)力變形產(chǎn)生影響較小。

    參考文獻(xiàn)

    [1]高鐘璞.大壩基礎(chǔ)防滲墻[M].北京:中國(guó)電力出版社,2000.

    [2]王清友,孫萬(wàn)功,熊歡.塑性混凝土防滲墻[M].北京:中國(guó)水利水電出版社,2008.

    [3]張嘎,張建民.夾有泥皮粗粒土與結(jié)構(gòu)接觸面力學(xué)特性試驗(yàn)研究[J].巖土力學(xué),2005,26(9):1374-1378.

    [4]張治軍,饒錫保,王志軍.泥皮厚度對(duì)結(jié)構(gòu)接觸面力學(xué)特性影響的試驗(yàn)研究[J].巖土力學(xué),2008,29(9):2433-2438.

    [5]楊春鳴,邵生俊,劉鑫.考慮泥皮作用的深厚覆蓋層壩基防滲墻應(yīng)力變形特性分析[J].西北農(nóng)林科技大學(xué)學(xué)報(bào),2013,41(8):1~10.

    [6]Goodman R F,Taylor R L,Brekke T L. A Model for the Mechanics of Jointed Rock[J]. Journal of the Soil Mechanics and Foundation Div. ASCE,1968,94(3):637-660.

    [7]費(fèi)康,張建偉.ABAQUS在巖土工程中的應(yīng)用[M].北京:中國(guó)水利水電出版社,2010.

    [8]Clough,G.W.Duncan,J. M. Finite Element Analysis of Retaining Wall Behavior[J]. Journal of the Soil Mechanics and Foundation Div. ,ASCE,1971,97(12):1657-1674.

    猜你喜歡
    土石壩覆蓋層防滲墻
    深水淺覆蓋層傾斜巖面河床圍堰設(shè)計(jì)及應(yīng)用
    聲子晶體覆蓋層吸聲機(jī)理研究
    無(wú)限元法在深覆蓋層土石壩動(dòng)力分析中的應(yīng)用
    淺談新舊防滲墻接頭連接處理施工技術(shù)
    淺薄覆蓋層傾斜巖面大直徑鋼護(hù)筒施工方案比選及應(yīng)用
    試述中小型水庫(kù)的土石壩除險(xiǎn)加固施工技術(shù)應(yīng)用
    水利施工中土石壩施工技術(shù)的應(yīng)用實(shí)踐
    關(guān)于水利水電工程中土石壩施工技術(shù)的探討
    水電站圍堰防滲墻中鉆劈法施工應(yīng)用探究
    面向水庫(kù)土石壩的除險(xiǎn)加固技術(shù)研究
    科技資訊(2015年21期)2015-11-14 19:27:16
    国产成人一区二区三区免费视频网站| 日韩欧美在线二视频| 50天的宝宝边吃奶边哭怎么回事| 国产精品一区二区精品视频观看| 色在线成人网| 免费高清在线观看日韩| xxxwww97欧美| 色综合婷婷激情| 色哟哟哟哟哟哟| 国产欧美日韩一区二区三| 最近最新中文字幕大全免费视频| 国产v大片淫在线免费观看| 亚洲av片天天在线观看| 亚洲av成人不卡在线观看播放网| 亚洲人成77777在线视频| 1024手机看黄色片| 日韩欧美在线二视频| 国产在线精品亚洲第一网站| 中文在线观看免费www的网站 | 亚洲成av人片免费观看| 中国美女看黄片| 午夜福利视频1000在线观看| 精品熟女少妇八av免费久了| 99热6这里只有精品| 特大巨黑吊av在线直播 | 国产真人三级小视频在线观看| 一区二区日韩欧美中文字幕| 亚洲男人天堂网一区| 国产精品日韩av在线免费观看| 日本免费a在线| 18禁黄网站禁片午夜丰满| 欧美在线一区亚洲| 久久久久久久精品吃奶| 亚洲中文av在线| 美女午夜性视频免费| 欧美又色又爽又黄视频| 成人三级黄色视频| 少妇被粗大的猛进出69影院| 国产成年人精品一区二区| 色综合站精品国产| 国产成人精品无人区| 给我免费播放毛片高清在线观看| 美女 人体艺术 gogo| 99久久综合精品五月天人人| 亚洲天堂国产精品一区在线| 精品一区二区三区视频在线观看免费| 女生性感内裤真人,穿戴方法视频| 99在线视频只有这里精品首页| 观看免费一级毛片| 麻豆久久精品国产亚洲av| 久久天堂一区二区三区四区| 欧美成狂野欧美在线观看| 高清在线国产一区| 亚洲熟女毛片儿| 在线永久观看黄色视频| a在线观看视频网站| 午夜久久久在线观看| 757午夜福利合集在线观看| 男女午夜视频在线观看| 亚洲一码二码三码区别大吗| 99久久综合精品五月天人人| 日韩精品中文字幕看吧| 国产精品香港三级国产av潘金莲| 两个人视频免费观看高清| 精品无人区乱码1区二区| 麻豆一二三区av精品| 一区二区三区高清视频在线| 国产精品日韩av在线免费观看| 黄色片一级片一级黄色片| 日日爽夜夜爽网站| 国产高清视频在线播放一区| 亚洲人成77777在线视频| 免费在线观看亚洲国产| 色婷婷久久久亚洲欧美| 在线永久观看黄色视频| 一区二区三区精品91| 午夜成年电影在线免费观看| 国产午夜福利久久久久久| 99热只有精品国产| 一进一出抽搐动态| 亚洲专区字幕在线| 一区福利在线观看| 成年人黄色毛片网站| 日本一本二区三区精品| 亚洲成人免费电影在线观看| 久久久久久国产a免费观看| 可以在线观看的亚洲视频| 日本一本二区三区精品| 国产成人欧美| 久久久水蜜桃国产精品网| 亚洲人成电影免费在线| 亚洲一区二区三区不卡视频| av欧美777| 中国美女看黄片| 最近最新中文字幕大全电影3 | 日韩欧美免费精品| www.精华液| 亚洲专区国产一区二区| 亚洲国产欧洲综合997久久, | 亚洲国产毛片av蜜桃av| 日本 欧美在线| 高清在线国产一区| 少妇粗大呻吟视频| 精品国产乱码久久久久久男人| 精品午夜福利视频在线观看一区| 又紧又爽又黄一区二区| 久久久水蜜桃国产精品网| 中文字幕精品亚洲无线码一区 | 99在线人妻在线中文字幕| 欧美日韩一级在线毛片| 国产成人精品久久二区二区免费| 亚洲精品在线观看二区| 一夜夜www| 欧美另类亚洲清纯唯美| 午夜视频精品福利| 性色av乱码一区二区三区2| 老司机福利观看| 国产精品99久久99久久久不卡| 免费av毛片视频| 在线观看免费日韩欧美大片| 人人妻人人看人人澡| 久久国产乱子伦精品免费另类| 18禁美女被吸乳视频| svipshipincom国产片| 在线观看免费视频日本深夜| 久久精品影院6| 国产精品美女特级片免费视频播放器 | 久久精品成人免费网站| 国产精品野战在线观看| 啦啦啦免费观看视频1| 一个人观看的视频www高清免费观看 | 日韩欧美 国产精品| 视频区欧美日本亚洲| 国产视频内射| 国产精品免费视频内射| 91成年电影在线观看| 成人国语在线视频| 精品少妇一区二区三区视频日本电影| 国产av不卡久久| 中文字幕精品亚洲无线码一区 | 露出奶头的视频| 国产亚洲欧美98| 色婷婷久久久亚洲欧美| 国产精品,欧美在线| 国产欧美日韩一区二区精品| 免费在线观看视频国产中文字幕亚洲| 黄片小视频在线播放| 别揉我奶头~嗯~啊~动态视频| 国产亚洲精品综合一区在线观看 | 亚洲七黄色美女视频| 国产精品 欧美亚洲| 亚洲精品国产区一区二| 国内毛片毛片毛片毛片毛片| 黄片播放在线免费| 波多野结衣av一区二区av| 99久久无色码亚洲精品果冻| 一级毛片精品| 色在线成人网| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品av在线| 久久国产精品男人的天堂亚洲| 麻豆av在线久日| 国产精品电影一区二区三区| 搡老岳熟女国产| 一边摸一边抽搐一进一小说| 亚洲一区二区三区不卡视频| 免费高清视频大片| av片东京热男人的天堂| 在线观看免费视频日本深夜| 欧美日韩乱码在线| 高清在线国产一区| 精品人妻1区二区| 欧美zozozo另类| 久久午夜综合久久蜜桃| 亚洲美女黄片视频| 一本久久中文字幕| 国产伦在线观看视频一区| 久久亚洲精品不卡| 可以在线观看的亚洲视频| 欧美精品啪啪一区二区三区| 女性生殖器流出的白浆| 色尼玛亚洲综合影院| 啪啪无遮挡十八禁网站| 亚洲欧美日韩无卡精品| 搞女人的毛片| 欧美 亚洲 国产 日韩一| 国产激情偷乱视频一区二区| 老司机深夜福利视频在线观看| 久久久久久人人人人人| 99国产精品一区二区蜜桃av| 精品第一国产精品| 欧美又色又爽又黄视频| 成人午夜高清在线视频 | 熟妇人妻久久中文字幕3abv| 久久中文看片网| 很黄的视频免费| 日韩国内少妇激情av| 亚洲激情在线av| 国产极品粉嫩免费观看在线| 亚洲av五月六月丁香网| 一个人免费在线观看的高清视频| 高清在线国产一区| 国产亚洲精品第一综合不卡| 亚洲色图 男人天堂 中文字幕| 亚洲自偷自拍图片 自拍| 1024香蕉在线观看| 欧美成人免费av一区二区三区| 欧美亚洲日本最大视频资源| 久久久久久久精品吃奶| 国产1区2区3区精品| 久久国产乱子伦精品免费另类| 欧美黑人巨大hd| 国产激情欧美一区二区| 亚洲精品久久成人aⅴ小说| 欧美色视频一区免费| 精品卡一卡二卡四卡免费| 午夜福利一区二区在线看| 午夜老司机福利片| 国产亚洲av高清不卡| 欧美性猛交╳xxx乱大交人| 欧美大码av| 国产精品久久久久久亚洲av鲁大| 亚洲色图 男人天堂 中文字幕| 在线十欧美十亚洲十日本专区| 天堂√8在线中文| 国产成人欧美在线观看| 亚洲精品国产一区二区精华液| ponron亚洲| 中亚洲国语对白在线视频| 国产成人精品久久二区二区免费| 亚洲熟妇中文字幕五十中出| 成人18禁在线播放| 巨乳人妻的诱惑在线观看| 禁无遮挡网站| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文av在线| 欧美另类亚洲清纯唯美| 村上凉子中文字幕在线| 中文字幕人妻熟女乱码| 中文亚洲av片在线观看爽| 日韩视频一区二区在线观看| 啦啦啦观看免费观看视频高清| 国产真人三级小视频在线观看| 最近在线观看免费完整版| 岛国在线观看网站| 91av网站免费观看| 精品国产国语对白av| 亚洲第一青青草原| 精品久久久久久成人av| 亚洲精品中文字幕在线视频| 亚洲精品久久国产高清桃花| 超碰成人久久| 欧美性猛交╳xxx乱大交人| 亚洲在线自拍视频| 午夜免费观看网址| 欧美日韩亚洲国产一区二区在线观看| 色在线成人网| 在线看三级毛片| 精品一区二区三区视频在线观看免费| 1024视频免费在线观看| 亚洲av五月六月丁香网| 日本a在线网址| 精品久久久久久久人妻蜜臀av| 久久天堂一区二区三区四区| 日韩欧美在线二视频| 日本五十路高清| 国产精华一区二区三区| 欧美日韩精品网址| 满18在线观看网站| 亚洲av成人一区二区三| 琪琪午夜伦伦电影理论片6080| 91大片在线观看| 国产真人三级小视频在线观看| 国产成人一区二区三区免费视频网站| 色哟哟哟哟哟哟| 99国产精品一区二区蜜桃av| 欧美性长视频在线观看| 91av网站免费观看| 国产亚洲精品一区二区www| 99精品欧美一区二区三区四区| 色综合欧美亚洲国产小说| 天天躁夜夜躁狠狠躁躁| 后天国语完整版免费观看| 国产精品98久久久久久宅男小说| cao死你这个sao货| 亚洲国产精品合色在线| 91av网站免费观看| 美国免费a级毛片| 国产欧美日韩精品亚洲av| 色尼玛亚洲综合影院| 国产精品久久久av美女十八| 国产日本99.免费观看| 亚洲中文字幕一区二区三区有码在线看 | 免费无遮挡裸体视频| 黑丝袜美女国产一区| 国产精品久久电影中文字幕| 亚洲va日本ⅴa欧美va伊人久久| e午夜精品久久久久久久| 两个人看的免费小视频| 夜夜看夜夜爽夜夜摸| 久久精品国产亚洲av香蕉五月| 亚洲一卡2卡3卡4卡5卡精品中文| 国产一区二区三区视频了| 国产精品 国内视频| 国产伦在线观看视频一区| 天堂动漫精品| 亚洲欧美日韩高清在线视频| 中文字幕人妻丝袜一区二区| 淫妇啪啪啪对白视频| 两个人免费观看高清视频| 性欧美人与动物交配| 欧美又色又爽又黄视频| 免费在线观看成人毛片| 日本五十路高清| 免费在线观看黄色视频的| 亚洲精品中文字幕在线视频| 男人操女人黄网站| 悠悠久久av| 青草久久国产| 亚洲精品国产区一区二| 欧美色视频一区免费| 亚洲 欧美 日韩 在线 免费| 日韩大码丰满熟妇| 日韩中文字幕欧美一区二区| 成人手机av| 婷婷精品国产亚洲av| 男女之事视频高清在线观看| 少妇熟女aⅴ在线视频| 久久久久久大精品| 9191精品国产免费久久| 亚洲av电影在线进入| 久久国产亚洲av麻豆专区| 天堂√8在线中文| 黄片播放在线免费| 亚洲精品色激情综合| 在线播放国产精品三级| 国产精品久久久久久精品电影 | 丁香欧美五月| 精品午夜福利视频在线观看一区| 精品国产一区二区三区四区第35| 亚洲avbb在线观看| 午夜a级毛片| 欧美成人免费av一区二区三区| 久久国产乱子伦精品免费另类| 国产日本99.免费观看| 日韩欧美国产一区二区入口| 非洲黑人性xxxx精品又粗又长| 日韩免费av在线播放| 嫁个100分男人电影在线观看| 亚洲专区国产一区二区| 1024视频免费在线观看| 夜夜看夜夜爽夜夜摸| 十八禁人妻一区二区| 一二三四在线观看免费中文在| 国产伦人伦偷精品视频| 又黄又爽又免费观看的视频| 一级片免费观看大全| 老司机在亚洲福利影院| 99国产综合亚洲精品| 国产人伦9x9x在线观看| 国产精品久久久久久精品电影 | 一区福利在线观看| 人人妻人人澡人人看| 伊人久久大香线蕉亚洲五| 给我免费播放毛片高清在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成电影免费在线| 亚洲 国产 在线| 99久久国产精品久久久| 午夜福利在线在线| 久久精品国产99精品国产亚洲性色| 人人妻,人人澡人人爽秒播| bbb黄色大片| 日韩欧美国产一区二区入口| 日本撒尿小便嘘嘘汇集6| 最好的美女福利视频网| 老司机午夜十八禁免费视频| 亚洲国产毛片av蜜桃av| 成人亚洲精品av一区二区| 久久精品夜夜夜夜夜久久蜜豆 | 午夜激情av网站| 久久久久国产精品人妻aⅴ院| 久久婷婷人人爽人人干人人爱| 亚洲色图 男人天堂 中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久国产一级毛片高清牌| 亚洲五月色婷婷综合| 日韩视频一区二区在线观看| 国产精品久久电影中文字幕| 亚洲熟女毛片儿| 亚洲男人天堂网一区| 欧美色视频一区免费| 久久青草综合色| 亚洲专区中文字幕在线| 国产成人精品无人区| 国产99白浆流出| 满18在线观看网站| 一级毛片高清免费大全| 宅男免费午夜| 欧美黄色淫秽网站| 久久久久久人人人人人| 精品久久久久久久末码| netflix在线观看网站| 深夜精品福利| e午夜精品久久久久久久| 亚洲精品国产区一区二| 亚洲国产欧美一区二区综合| 精品高清国产在线一区| 国产精品免费一区二区三区在线| 免费看日本二区| www.自偷自拍.com| 成年人黄色毛片网站| av免费在线观看网站| 99国产综合亚洲精品| 两个人看的免费小视频| 国产免费av片在线观看野外av| www.熟女人妻精品国产| 欧美黑人欧美精品刺激| 免费一级毛片在线播放高清视频| 日韩视频一区二区在线观看| 成人国产综合亚洲| 久久精品成人免费网站| 亚洲精品国产区一区二| 老司机午夜十八禁免费视频| 国产不卡一卡二| 国产精品 欧美亚洲| 亚洲激情在线av| 高潮久久久久久久久久久不卡| 女生性感内裤真人,穿戴方法视频| 真人做人爱边吃奶动态| 欧美大码av| 88av欧美| 亚洲午夜理论影院| 狠狠狠狠99中文字幕| 老熟妇乱子伦视频在线观看| 十八禁人妻一区二区| 12—13女人毛片做爰片一| 日本一区二区免费在线视频| 国产精品久久久av美女十八| 国产片内射在线| 99热只有精品国产| 中文字幕av电影在线播放| 国产精品爽爽va在线观看网站 | 亚洲男人天堂网一区| 国产精品野战在线观看| 免费看美女性在线毛片视频| 欧美精品亚洲一区二区| 香蕉国产在线看| 波多野结衣高清无吗| 成人国语在线视频| 久久婷婷人人爽人人干人人爱| 人人澡人人妻人| 色综合站精品国产| 看片在线看免费视频| 黄色视频不卡| 熟女少妇亚洲综合色aaa.| 国产又黄又爽又无遮挡在线| 国产av一区二区精品久久| 亚洲avbb在线观看| 亚洲美女黄片视频| 可以在线观看的亚洲视频| 在线观看免费视频日本深夜| 亚洲精品中文字幕一二三四区| 精品久久久久久久久久免费视频| 国产主播在线观看一区二区| 欧美一级毛片孕妇| 久久久久久九九精品二区国产 | 亚洲av熟女| 成人永久免费在线观看视频| 久久久久久久午夜电影| 国产精品野战在线观看| 日本精品一区二区三区蜜桃| 狂野欧美激情性xxxx| 91字幕亚洲| 99精品久久久久人妻精品| 午夜视频精品福利| 国产熟女午夜一区二区三区| 日韩中文字幕欧美一区二区| 午夜福利在线观看吧| 热re99久久国产66热| 99在线人妻在线中文字幕| 国产又黄又爽又无遮挡在线| 两个人免费观看高清视频| 亚洲欧美精品综合久久99| 国产三级在线视频| 麻豆成人午夜福利视频| 亚洲熟女毛片儿| 视频在线观看一区二区三区| 色婷婷久久久亚洲欧美| 搡老妇女老女人老熟妇| 久久精品亚洲精品国产色婷小说| 极品教师在线免费播放| 亚洲国产日韩欧美精品在线观看 | 亚洲色图 男人天堂 中文字幕| 亚洲三区欧美一区| 十八禁人妻一区二区| 一边摸一边做爽爽视频免费| 日本 欧美在线| 黄色成人免费大全| 久久亚洲真实| 少妇 在线观看| 亚洲精品久久国产高清桃花| 国产视频一区二区在线看| 久久精品91无色码中文字幕| 亚洲午夜精品一区,二区,三区| 日韩大码丰满熟妇| 男人舔奶头视频| 国产久久久一区二区三区| 久久香蕉国产精品| 在线十欧美十亚洲十日本专区| 青草久久国产| 日韩成人在线观看一区二区三区| 日韩欧美国产在线观看| 成在线人永久免费视频| 欧美成人午夜精品| 亚洲人成伊人成综合网2020| 一进一出好大好爽视频| 亚洲中文日韩欧美视频| 麻豆av在线久日| 日韩欧美国产一区二区入口| 国产av一区在线观看免费| 国产一区二区在线av高清观看| 欧美 亚洲 国产 日韩一| 成人手机av| 无遮挡黄片免费观看| 男女那种视频在线观看| 99国产综合亚洲精品| 熟女电影av网| 老司机午夜十八禁免费视频| 午夜免费观看网址| 香蕉久久夜色| 久久九九热精品免费| 在线视频色国产色| 中文字幕久久专区| 国产成人精品无人区| 一区二区三区高清视频在线| 色综合欧美亚洲国产小说| 99在线人妻在线中文字幕| 色老头精品视频在线观看| 99精品久久久久人妻精品| 男人舔女人下体高潮全视频| 久久国产亚洲av麻豆专区| 亚洲国产欧美一区二区综合| 神马国产精品三级电影在线观看 | 午夜两性在线视频| 亚洲一区高清亚洲精品| 亚洲,欧美精品.| 91九色精品人成在线观看| 亚洲精品久久国产高清桃花| 国产亚洲精品第一综合不卡| 淫妇啪啪啪对白视频| 狠狠狠狠99中文字幕| ponron亚洲| 国产精品 欧美亚洲| 露出奶头的视频| 一级a爱视频在线免费观看| 很黄的视频免费| 国产精品,欧美在线| 中文资源天堂在线| 亚洲天堂国产精品一区在线| 听说在线观看完整版免费高清| 午夜两性在线视频| www国产在线视频色| 99在线人妻在线中文字幕| 亚洲av五月六月丁香网| 很黄的视频免费| 国产精品免费一区二区三区在线| 熟妇人妻久久中文字幕3abv| 欧美人与性动交α欧美精品济南到| 男女下面进入的视频免费午夜 | 国产精品一区二区三区四区久久 | 天天添夜夜摸| 在线观看www视频免费| 久久久国产成人精品二区| 日韩欧美三级三区| 少妇熟女aⅴ在线视频| 操出白浆在线播放| 精品国产美女av久久久久小说| 99久久99久久久精品蜜桃| 在线观看免费日韩欧美大片| 无遮挡黄片免费观看| 脱女人内裤的视频| 亚洲五月天丁香| 91成人精品电影| 久久青草综合色| 中亚洲国语对白在线视频| 2021天堂中文幕一二区在线观 | 99国产综合亚洲精品| 午夜免费观看网址| www.熟女人妻精品国产| 欧美激情高清一区二区三区| 亚洲精品在线美女| 亚洲avbb在线观看| 国产精品精品国产色婷婷| 亚洲精品中文字幕在线视频| 日韩欧美在线二视频| 久久中文字幕一级| 一区二区三区高清视频在线| 国产av一区在线观看免费| 免费女性裸体啪啪无遮挡网站| 欧美一区二区精品小视频在线| 精品久久久久久久久久久久久 | 久久 成人 亚洲| 国产成人影院久久av| www日本黄色视频网| 午夜日韩欧美国产| 香蕉国产在线看| 久热这里只有精品99| 无限看片的www在线观看| 免费在线观看黄色视频的| 免费电影在线观看免费观看| 搞女人的毛片| 国产成人av教育| 黄片小视频在线播放| 757午夜福利合集在线观看|