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

    一維滲透力與浮力1)

    2017-11-11 01:55:28丁洲祥
    力學學報 2017年5期
    關鍵詞:滲透力土力學浮力

    丁洲祥

    ?(北京交通大學土木建筑工程學院,北京100044)

    ?(陜西省公路橋梁與隧道重點試驗室,西安710064)

    一維滲透力與浮力1)

    丁洲祥?,?,2)

    ?(北京交通大學土木建筑工程學院,北京100044)

    ?(陜西省公路橋梁與隧道重點試驗室,西安710064)

    為簡化分析,針對一維滲流問題,研究提出了土力學中滲透力和浮力的兩種推導方法,作為傳統(tǒng)的宏觀尺度孔隙水隔離體法的有益補充.彈性力學平衡微分方程、土力學Terzaghi有效應力方程和流體力學簡化Bernoulli方程構成本文分析滲透力和浮力的3個基本方程.在基本方程基礎上,容易推導出相應的骨架和孔隙水兩種隔離體的平衡微分方程,從而在靜力平衡范疇內揭示滲透力和浮力的內涵.單位體積飽和土體的滲透力,源于總水頭壓力的梯度,而浮力則源于位置水頭壓力在豎向的梯度,這兩者統(tǒng)一于骨架或孔隙水的平衡微分方程.實際工程關注的有效應力計算問題,一般可以直接應用3個基本方程來確定;只有在簡化條件下可按滲透力和浮力計算土體中有效應力分布規(guī)律.還討論了若干研究熱點問題,重點探討了當前一種滲透力新定義j=nγwi在形式上的合理性以及在實際應用中可能存在的風險,并驗證了一維滲透力的一種經典精細化表述結果中考慮滲流速度時間導數(shù)的嚴謹性.在土力學滲透力和浮力問題研究中應重視和正確應用Terzaghi有效應力方程.

    滲透力,浮力,有效應力原理,孔隙水壓力,平衡微分方程

    引言

    經典土力學[13]中滲透力和浮力的概念并不簡單,尤其是前者.過去一年多來,文獻[4-7]從機理和公式推導等角度對這些既有概念加以研究,但尚無終論,如蔣中明等[6]認為“滲透力概念定義的不統(tǒng)一性以及推導過程的多樣性、復雜性,是導致人們對滲透力的理解和應用產生誤解的根源”.滲透力和浮力概念再度成為研究熱點,這并非人們有意將簡單問題復雜化,而是因為大量的既有相關認識仍有不足乃至相互矛盾之處.這輪研究和反思的一個主要特點是,注重研究方法的一般性、系統(tǒng)性和推導過程的嚴謹性.自2013年來,筆者對滲透力和浮力概念格外加以關注,并做了一些探討性研究工作;現(xiàn)在意識到,可能存在一種觀點——有時在某種意義上,土力學應用中使用的滲透力和浮力概念既不宜絕對肯定地當成唯一正確的方法,也不宜無視其合理性及意義.

    見仁見智的學術研究,一般強調言貴有據(jù).傳統(tǒng)孔隙水隔離體法[19]將孔隙水假想地從土水整體中隔離開來,然后考慮浮力和滲透力兩者的反力,據(jù)此分析孔隙水隔離體平衡條件,其中滲透力和浮力貌似不可或缺,容易被絕對肯定化.此外,滲透力和浮力概念在傳統(tǒng)土力學中沿用多年,在評價滲透破壞等方面也有相關應用,這進一步強化了上述看法.然而,歷史上也不乏相反的學術觀點,對滲透力和浮力的研究意義基本上持否定的態(tài)度.

    目前少有文獻辯證地分析滲透力和浮力的上述絕對肯定與否定的觀點,這就需要我們謹慎地求證開展辯證分析的力學依據(jù).這樣做的意義是,從新視角探求滲透力和浮力含義與適用條件的同時,也發(fā)現(xiàn)了推導土力學滲透力的一種新方法,該方法有別于目前仍占主流的孔隙水隔離體法.為避免將問題復雜化,本文僅針對一維問題進行分析:首先,提出滲透力和浮力的3個基本方程;然后,導出骨架和孔隙水兩種隔離體平衡微分方程;其次,面向應用討論有效應力計算中的等價性問題;最后,就同行最新研究成果進行了對比探討.所得結論有助于深化對滲透力和浮力概念的認識.

    1 一維模型和3個基本方程

    飽和土體一維滲流模型如圖 1所示.其中,Z坐標軸向上為正,坐標原點位于同質土層的底面z=0;土層厚度為H;土層底面、頂面和坐標z的水頭分別為 h(0)、h(H)和 h(z).在圖示條件下,h(0)>h(H),故滲流方向豎直向上,但模型也適用于h(0)≤h(H)的情況.坐標增量(z,z+dz)微元段dz內土體受力情況:上下端面對應應力增量(σ,σ+dσ),體力為fzdz.

    圖1 滲透力和浮力的一維飽和土模型Fig.1 One-dimensional saturated soil model for seepage force and buoyance

    1.1 總應力平衡微分方程

    在土水整體總應力平衡微分方程方面,土力學和彈性力學是相通的——如果將土體固相在空間占據(jù)的區(qū)域明確定義為土體連續(xù)介質構形[10].這樣可沿用彈性力學中一維靜力平衡微分方程[11]

    式中,σ是飽和土體總應力;fz是單位體積土體的體力,fz=?γsat;σ,fz均服從彈性力學符號約定.

    1.2 Terzaghi有效應力方程

    在固液兩相應力分配關系上,土力學很有特色.經典土力學經常采用Terzaghi有效應力方程[3]

    式中,σ′是土體有效應力;u孔隙水壓力;σ,σ′和u按土力學習慣,常以受壓為正.

    1.3 孔隙水壓力方程

    在考慮滲流勢能時,土力學借鑒了流體力學的部分經典成果.當忽略速度水頭時,流體力學Bernoulli方程[12]更加簡化,由此可得到

    式中,h=h(z)是總水頭;z是土中考察點的位置,z=0是零勢參考面;γw是孔隙水單位體積的重度.

    本文認為,上述3個基本方程是研究傳統(tǒng)滲透力和浮力的聯(lián)合控制方程.

    2 骨架平衡微分方程與滲透力和浮力

    2.1 骨架平衡微分方程與浮力

    將式(2)和式(3)代入式(1)的過程中,注意應力符號約定的學科差異,不妨以壓為正,則

    式中,等號左端第3項γw其實就是單位體積飽和土體受到的浮力,它源于位置水頭壓力在豎向的梯度,在本文中以沿坐標軸正向為正(如果只看其模,就是俗稱的孔隙水重度),而左端第2項則表示滲流的影響.

    具體分析還涉及土力學中兩個相關概念.

    概念1:有1種水力梯度定義

    在數(shù)學上表示總水頭h的梯度對應的負矢量,反映驅動孔隙水滲流的一種作用.

    概念2:根據(jù)三相關系,飽和土體有效重度或浮重度(e ff ective or buoyant unit weight)[13]

    按上述兩概念,改寫式(4),得到式(1)的等價形式

    式中,若γwi和?γ′數(shù)值為正,則表示沿坐標軸正向.

    可見,單位體積的浮力與總體力 fz都是矢量,兩者合并成單位體積有效體力,記.也就是說,飽和土的浮力被整合到有效體力或通常的浮重度概念中,這較易理解.

    2.2 骨架平衡微分方程與滲透力

    接下來的問題是:滲透力在哪里?

    需要先回到土力學界對滲透力的常規(guī)理解.不妨以國外近期和較早的著作為例.Das[14]認為:“Flow ofwaterthroughasoilmassresultsinsomeforcebeingexerted on the soil itself”.Budhu[13]認為:“As water fl ows through soil,it exerts a frictional drag on the soil particles,resulting in head losses.The frictional drag is called seepage force in soil mechanics”.Lambe 等[15]認為:“The seepage force is applied by the moving water to the soil skeleton through frictional drag”.可見,滲透力定義至少有兩個要求:孔隙水滲流;對土體或骨架能產生某種作用.

    分析式(7)可知,該式左端第2項,剛好能滿足這兩個條件,所以有理由將其定義為一種滲透力

    它反映了滲流場對有效應力場的耦合作用.

    事實上,式(8)就是傳統(tǒng)土力學最常見的滲透力表述或推導結果.本文從骨架受力平衡角度,提供了考察滲透力和浮力的不同于傳統(tǒng)孔隙水隔離體法的新視角,也可作為推導滲透力方法的新補充.這種方法不失簡潔和意義深刻.

    利用式(8),可使式(7)的數(shù)學表述更簡潔

    式中,單位體積土體浮力γw與飽和體力 fz合并成?γ′,即

    若土體孔隙中無水,則有 j= 0,γw= 0,γsat= γ′= γd(干重度),此時式 (9)退化為彈性力學一維平衡微分方程 (按彈力符號約定)式 (1),即?σ/?z+fz=0,此時 σ = σ′,fz= ?γd,這表示干土的靜力平衡.

    3 孔隙水平衡微分方程與滲透力和浮力

    3.1 孔隙水平衡微分方程與浮力

    經典土力學中的孔隙水隔離體法[19]主要是借助一定宏觀尺寸的力學模型進行分析,而很少采用連續(xù)介質力學中局部化形式的平衡微分方程法.

    采用土力學應力符號約定,總應力平衡微分方程式(1)改寫為

    將總應力平衡微分方程(10)減去骨架平衡微分方程式(9),得到

    將有效應力方程式(2)代入式(11),得到

    等價簡化式(12),得到孔隙水平衡微分方程為

    式中,等號左端第3項和第2項分別表示土體受到的浮力反力以及滲透力反力作用.?γw和?j的數(shù)值為負值,表示反力方向沿坐標軸負向.

    當土體中無滲流時,孔隙水成為靜水,一般應有j=0.此時,式(13)就簡化為流體力學中一維的歐拉方程,即均勻重力場中靜止流體的一維平衡方程.

    若在一維歐拉方程中進一步忽略浮力即不計重力場對孔隙水的作用,則式(13)退化為?u/?z=0,此時孔隙水壓力u等于與坐標無關的某個常數(shù).這意味著孔壓在流體中每一點都相同;對于密閉流體而言,此即帕斯卡的靜壓傳遞原理.

    3.2 孔隙水平衡微分方程與滲透力

    現(xiàn)假定滲透力 j未知,則根據(jù)孔隙水平衡方程式(13),可推導出 j的具體表達式

    再依據(jù)孔隙水壓力方程式(3),得到滲透力為

    若同樣采用式(5)的水力梯度定義,則式(15)就完全等價于式(8),即 j=γwi.

    這樣就按主流的孔隙水隔離體法思想實現(xiàn)了對滲透力的新推導.

    重點分析以下3點:(1)結果對比表明,無論采用骨架平衡式(7)還是孔隙水平衡方程式(13),都可以非常一致地推出滲透力 j的表達式(8);而傳統(tǒng)土力學研究中通常忽略了基于骨架平衡的分析方法,這似乎是不夠理性和全面的.(2)對比式(13)與式(9),在骨架和孔隙水隔離體平衡微分方程中,浮力γw和滲透力 j前面的符號是相反的:這是相間的作用力和反作用力屬性的體現(xiàn).(3)如果將式(9)和式(13)兩端相加,則得到兩相整體總應力平衡方程式(10)或式(1)(按彈力應力符號約定)——兩式中均不出現(xiàn)滲透力和浮力,因為兩者都是固液兩相間成對出現(xiàn)的內部作用力.按經典力學觀點,如果對系統(tǒng)中所有質點受到的力求和,則內力將相互抵消而僅存外力.對本文的土水整體,外力只包含重力,因而在兩相整體平衡方程中外力僅有飽和體力fz(fz=?γsat).

    4 應用中的數(shù)學等價性

    在實際應用中,滲透力和浮力對骨架或固相的作用主要體現(xiàn)在對有效應力的影響方面.現(xiàn)針對一維有效應力計算,采用3種不同算法,簡要對比分析在采用/棄用滲透力與浮力概念情況下的計算結果.

    4.1 基于滲透力和浮力的有效應力算法

    當滲透力 j和浮重力γ′為常數(shù)時,式(9)的不定積分結果為

    式中,C是積分常數(shù),取決于地表載荷條件.如果設地表 z=H為排水邊界、地表邊界處有效應力σ′(H)=q以及孔壓為u(H)=uH,則式(9)的定積分結果即有效應力分布為

    有些土力學著作[16]給出的不同方向滲流作用下有效應力的計算公式,其實都可以統(tǒng)一成式(17).

    然而,當滲透力 j隨坐標z變化,式(16)和式(17)將不再適用.這體現(xiàn)了傳統(tǒng)滲透力和浮力算法的適用范圍.

    4.2 基于總水頭和浮力的有效應力算法

    考慮 j為變數(shù)而γ′仍為常數(shù)的情況,此時無法直接對j按常量積分,但可根據(jù)式(4)求出被積函數(shù)?σ′/?z的不定積分

    式中,D是積分常數(shù),依賴于地表載荷情況.引入上述邊界條件后的定積分結果為

    有效應力分布式(18)和式(19)具有更普遍的意義.當滲透力j為常數(shù)時,兩者分別能夠退化為基于j的有效應力函數(shù)式(16)和式(17).此時兩種算法在數(shù)學上是等價變換.然而,式(18)和式(19)根本不涉及滲透力,而主要取決于總水頭函數(shù)h(z)和有效體力 ?γ′.

    這初步表明,在進行有效應力計算時,滲透力并非不可或缺的中間量;個別情況下,滲透力算法的式(16)和式(17)反而會縮小問題的求解范圍.

    4.3 基于3個基本方程直接積分的有效應力算法

    對同一定解問題和邊界條件,此處按式(1)~式(3)進行求解,而不用從滲透力和浮力概念出發(fā).

    對總應力平衡方程式(1)進行不定積分,并按土力學應力習慣用法有

    式中,E為積分常數(shù).引入同樣的邊界條件,得到總應力的定積分

    根據(jù)有效應力方程式(2),得到有效應力

    將式(21)和孔壓方程式(3)代入式(22),同樣能得到有效應力沿埋深的分布函數(shù)σ′(z).該結果與式(19)完全相同(具體過程從略).

    鑒于上述論證,不妨據(jù)此推論:作為控制土體強度和變形的有效應力,其實不必根據(jù)額外考慮滲透力和浮力的特殊方程來推導.

    然而,滲透力和浮力在理論與工程中畢竟得到了一定的普及和應用.在實際應用中,偏執(zhí)于對滲透力和浮力概念絕對肯定或否定的任一極端觀點可能都未必明智.應意識到使用滲透力和浮力尤其是前者必然帶來的局限性.

    5 討論

    5.1 傳統(tǒng)土力學對滲透力的機理解釋準確嗎?

    李廣信[4]認為滲透力是滲透水流對土骨架中的顆粒的推動力和拖曳力,并從土顆粒某種受力分析出發(fā),推導出滲透力的計算公式,即本文式(8).錢家歡等[8]僅將滲透水流施于單位土體內土粒上的拖曳力稱為滲透力,這與Lambe等[15]的觀點基本一致,同樣得到本文式(8).

    對機理的深入研究是有益的,但基于不同機理解釋得到同樣的滲透力表達結果,卻很難理解:如果滲透力結果與機理解釋的差別無關,那么已有的某些機理解釋是否準確和合理?值得思考.

    5.2 傳統(tǒng)土力學的滲透力需要修正嗎?

    2004年,Li和Ming[17]提出滲透力的另一不同表達

    式中,n是土體孔隙率.

    2014年,邵龍?zhí)兜萚7]也提出了與式(23)基本一樣的研究結果,并給出了更詳細的論證過程,認為傳統(tǒng)土力學得到的滲透力式(8)是錯誤的,需要按式(23)進行修正.2016年,李廣信在滲透力有關的文章[4]和著作中均沒有提到和評價文獻[7,17]的這些不同觀點與建議.

    因此,傳統(tǒng)土力學中滲透力式(8)是否需要修正,如果需要修正,又該如何完善?

    5.3 既有推導方法充分利用有效應力方程了嗎?

    2015年,陳愈炯[18]報道了自 1954年以來,Terzaghi有效應力原理在我國的應用歷史和部分實例,并“鼓勵同仁們重視該原理的應用,在應用過程中提高自己的洞察力和處理工程問題的能力,從而取得意想不到的好效果”.

    縱觀國內外主要的土力學著作,筆者尚未發(fā)現(xiàn)有明確利用Terzaghi有效應力原理或方程來完整推導滲透力的相關論述.在論文中,蔣中明等[6]使用Terzaghi有效應力方程,按張量方法進行了多維滲透力推導,但未就最終結果進行深入討論.

    本文將Terzaghi有效應力方程作為基本方程之一進行滲透力的推導,是應用有效應力原理或方程的有益嘗試.使人“意想不到”的是,本文方法有助于避開滲透力具體物理意義或機理解釋方面的矛盾和爭議,將傳統(tǒng)土力學滲透力的本質闡釋主要歸于Terzaghi有效應力原理與方程內在的基礎理論前提.

    有趣的是,1996年,Terzaghi等[19]認為,“although no theoretical basis for Eq.15.2(即本文式(2).——筆者注)has been found,its empirical basis is so well established that a quantitative knowledge for the interparticle reactions is not needed”.按該觀點,理論基礎的缺失并不影響Terzaghi有效應力方程的實用性,因為經驗性基礎可以取代對粒間反力的定量知識.在該意義上,傳統(tǒng)滲透力式(8)也應視為經驗性的Terzaghi有效應力方程聯(lián)合其他基本方程的一個具體應用.

    以下對比蔣中明等[6]對有效應力方程的應用情況.按其推導思想,在總應力平衡方程中,利用σ=σ+χu?χu(χ為任一實數(shù))數(shù)學恒等式以及飽和土物理指標關系:γsat= γd+nγw(γd是土體干重度),同時考慮本文其他兩個基本方程,可以推導出以下一維平衡微分方程

    當χ=n時,式(24)退化為

    式(25)對應于文獻[6]在豎直方向滲流的一維結果(符號約定稍有不同).式(25)等號左端第4項就是式(23)給出的滲透力新定義形式.但這是滲透力嗎?

    按文獻[6],式(25)左端第3項被“規(guī)定定義”為浮力,理由是“重力方向的壓力差或壓力梯度就是浮力”;第4項nγwi被規(guī)定定義為“代表拖曳力”;而滲透力被詮釋為由式(25)等號左端最后兩項共同構成,容易得到,式(25)左端最后兩項的和

    顯然,式(26)的結果并不僅受滲流的影響,將“浮力”與“拖曳力”之和定義為“滲透力”值得商榷.

    此外,文獻[6]的“浮力”定義對于靜水壓力狀態(tài)才成立;因為此時重力勢和孔壓勢梯度之和為零,所以也可以用重力勢梯度來表述.但對于一般的滲流狀態(tài),將總孔壓梯度作為浮力是不合適的,因為此時?u/?z受滲流影響是變數(shù).這意味著單位體積的“浮力”不是常數(shù).所以,這種分析方法似乎沒有充分闡明和區(qū)分浮力和滲透力的準確內涵,也不能作為式(23)具有合理性的證明.

    如果將變換后平衡方程中出現(xiàn)的任一γwi有關項,都稱為滲透力,如式(24)中的χγwi,似乎也是不合適的,因為χ是任一實數(shù).在骨架平衡方程中,當且僅當γwi有關項全部反映滲流的影響時,將其歸為滲透力才是合理的.而骨架平衡方程式(24)和(25)中?u/?z其實仍然與滲流相關.

    容易看出,式(25)與式(9)或式(7)實際上完全等價,而且都采用了Terzaghi有效應力方程.

    5.4 滲透力新定義式(23)還有其他解讀方法嗎?

    這里將詳細論述一種證明式(23)的可能方法.

    前述3個基本方程中,用Biot有效應力方程[20]

    替代Terzaghi有效應力方程式(2),按同樣思路和推導即分別得到骨架和孔隙水平衡方程的另一形式

    當孔隙率n隨埋深不變或其梯度為零時,式(28)和式(29)分別簡化為

    對比式(30)與骨架平衡微分方程式(9)可知,只有式(30)等號左端第2項與滲流有關.所以,形式上nγwi也同樣符合滲透力定義的兩種條件,即式(23)可稱為一種滲透力,這在形式上可以得到論證.

    盡管文獻[7,17]并沒有闡明采用Biot有效應力方程式(27),但根據(jù)上述分析,異于傳統(tǒng)的滲透力式(23)隱含了Biot有效應力方程,而其他有效應力方程無法在形式上得到式(23)的結果.

    5.5 滲透力新定義式(23)的應用會帶來哪些變化?

    滲透力的傳統(tǒng)表述式(8)與新定義式(23)明顯不同.新滲透力式(23)若要得到巖土工程界的普遍接受,除了最終需要工程實踐檢驗之外,還應事先在理論上體現(xiàn)其不同以往的特點或優(yōu)點.

    因此就有必要對本文基本方程之一的有效應力方程加以討論.如前所述,兩種滲透力結果差異的主要原因是分別采用了Terzaghi有效應力方程式(2)和Biot有效應力方程式(27).對同一問題而言,總應力σ和孔隙水壓u相對客觀,而有效應力無法直接測量,因而可以按無量綱方法分析不同有效應力方程對應的無量綱有效應力的差異性規(guī)律,見圖2.其中,式 (2)對應于 σ′/σ =1?u/σ,而式 (27)對應于σ′/σ =1?nu/σ.

    圖2 Biot和Terzaghi兩種有效應力方程的無量綱形式Fig.2 Dimensionless forms of e ff ective stress equations given by Biot and Terzaghi,respectively

    分析圖2可知,只有在無量綱孔壓u/σ為零時,兩種無量綱有效應力σ′/σ結果才相同.隨著u/σ的增大,兩者差別逐漸顯著,其中 σ′/σ的 Biot理論值總是大于u/σ相同情況下的Terzaghi理論值.當u/σ =1時,按 Terzaghi理論得到的 σ′/σ 為零,這符合通常的理解和認識,而且常用作發(fā)生流土或砂沸的臨界判別條件.但此時若按Biot理論,則有σ′/σ=1?n>0,不會發(fā)生流土.

    而當基于 Biot有效應力方程的無量綱有效應力 σ′/σ =0時,則要求對應的無量綱孔壓 u/σ =1/n>1.此時相應的Terzaghi有效應力是負值,即σ′/σ=(n?1)/n<0.這對經典土力學觀點而言,將出現(xiàn)比較奇怪的現(xiàn)象:對于不能承受拉應力的土體,此時早已發(fā)生破壞而可能不存在該狀態(tài).所以,在相同的總應力和孔壓條件下,相比Terzaghi有效應力方程,在判斷流土現(xiàn)象時,基于Biot有效應力方程的結果在理論上偏于不安全.

    進一步看,出現(xiàn)流土的臨界力學條件一般為σ′=0.按本文式(17),在均質土層的表面載荷q=0條件下,當滲透力 j=γwi=γ′時,土層中任意點σ′(z)≡0,由此推導出通常的臨界水力梯度

    因此,有

    在實際工程中,土體平均孔隙率有一定的變化范圍[21],對常見砂土 n一般為 0.29~0.44,而對通常黏性土n為0.29~0.52.再根據(jù)傳統(tǒng)的臨界水力梯度icr,T值[14]大多介于0.85和1.1之間的經驗,以及式(34),就可以得到icr,B對應的大致變化范圍為2.6~6.2.

    因此,icr,B/icr,T在2.4~7.3之間,相對于傳統(tǒng)土力學中的經驗結果,基于常見土體參數(shù)統(tǒng)計規(guī)律的臨界水力梯度經驗值將偏大很多.可見,滲透力新定義式(23)及其對應的icr,B用于流土(砂)可能性的破壞評價時,是非常不安全的.

    再結合張邦芾[22]近兩年的基坑滲流試驗結果予以驗證.在該試驗中,4組試驗細砂土樣平均孔隙率 n為 0.48~ 0.49,飽和重度 γsat為 18.21~18.05kN/m3,傳統(tǒng)的臨界水力梯度值icr,T為0.858~0.842.4種工況的坑內止水帷幕插入深度范圍內平均水力梯度i為1.013~1.153,大于icr,T;試驗最終觀察到各工況均發(fā)生流土破壞.按傳統(tǒng)滲透力定義式(8)和對應臨界水力梯度式(32)可以較好地評估此類滲透破壞.

    而如果按滲透力新定義式 (23)及對應的臨界水力梯度式(33),發(fā)生流土的臨界水力梯度icr,B為2.867~2.755.因為i<icr,B,所以試樣理論上應處于穩(wěn)定狀態(tài),但實際試驗中卻發(fā)生了流土破壞.這表明滲透力新定義及對應臨界水力梯度在這種情況下是不適用的,將其普遍用于指導工程實踐的安全風險可能很高.

    5.6 其他

    此外,當?shù)乇磔d荷q≠0時,傳統(tǒng)流土臨界條件j=γ′將不再適用,因為按式(17)有σ′(z)≡q≠0.這也反映了傳統(tǒng)滲透力概念在應用中的局限性,此時判別流土臨界狀態(tài),也可以不采用滲透力概念而直接按有效應力來處理.

    筆者還注意到,1958年,蘇聯(lián)學者格爾西萬諾夫和包利興[2]在研究地下水層流上升時的現(xiàn)象時,采用了如下方程

    該方程與式(9)有些相似,但文獻[2]既沒有闡明該式的由來,也沒有提到本文的3個基本方程,而且還將式(35)等號左端最后1項闡明為干容重(干重度).按式(9)來看,這可能是錯誤的,應該為浮重度.

    在推導滲透力(動水壓力)時,文獻[2]并不是從式(35)出發(fā),而是根據(jù)前述宏觀尺寸的孔隙水隔離體法,在考慮孔隙水慣性力的影響下得到滲透力表達式

    式中,q是滲流速度,在文獻[2]中寫為g可能是刊誤.式(36)等號右端第2項反映孔隙水慣性力影響.

    筆者基于已有的固結理論方面的成果[10],在忽略骨架速度和加速度的前提下[23],探討并提出一種能夠考慮液相速度和加速度的廣義多維滲透力

    如果在式(37)中引入以下簡化假設:(1)Terzaghi有效應力方程成立ηij=δij;(2)一維問題假定i=j=z;(3)孔隙流體為孔隙水ρf=ρw,則式(37)嚴格退化為

    不難發(fā)現(xiàn),化簡后的多維廣義滲透力式(38)與式(36)等價,即這樣,初步證明了格爾西萬諾夫和包利興提出的一維滲透力式(36)在理論上依然是嚴謹?shù)?,同時也揭示了其滲透力推導方法中同樣隱含使用了Terzaghi有效應力方程,而其他有效應力方程無法導出式(36).

    值得關注的還有,國內學者雷國輝等[2425]還研究了滲透力在瑞典條分法邊坡穩(wěn)定分析中應用的注意事項與局限性,并提出作用在骨架上的滲透力并非體積力的觀點等.但是,Pan等[26]在三維邊坡穩(wěn)定分析中,卻將滲透力視為作用于離散單元上的體積力進行數(shù)值計算.另外,Al-Madhhachi等[27]在黏性土侵蝕率模型中考慮了滲透力概念以改進經典Wilson模型的預測精度.這些研究中采用的滲透力主要是式(8)的傳統(tǒng)形式,而沒有使用更精確的式(36)或式(38).這表明,在工程實踐中,亟待加強滲透力的理論和應用研究.

    6 結論

    (1)針對一維工況,明確了分析飽和土滲透力和浮力的3個基本方程,由此分別得到骨架平衡微分方程式(9)和孔隙水隔離體平衡微分方程式(13).這兩種方法和思路都能夠準確而不失簡潔地推導出滲透力和浮力.傳統(tǒng)和主流的宏觀尺度孔隙水隔離體法不再是分析滲透力和浮力的唯一方法.

    (2)土中無水時,骨架平衡微分方程式(9)可以退化為彈性力學平衡微分方程(干土).不考慮滲透力(和浮力)時,孔隙水平衡微分方程(13)可以退化為流體力學歐拉方程(和帕斯卡靜壓傳遞方程).

    (3)一維滲流問題中的有效應力計算,在簡單情況下可以按滲透力和浮力分析即式(17).但一般情況下可以直接按本文3個基本方程求解即式(19),此時無需利用滲透力等概念.

    (4)對最新滲透力研究成果的討論工作表明:①證明傳統(tǒng)滲透力式(8)存在合理性的關鍵因素是充分和正確地應用3個基本方程中的Terzaghi有效應力方程;②從Biot有效應力方程及本文其他2個基本方程來看,其他學者提出的滲透力新定義j=nγwi在形式上是合理的,但從有限的工程實踐經驗和試驗成果來分析,該新定義若用于流土破壞評估等工程問題則可能是很不安全的,值得注意和警惕.

    (5)格爾西萬諾夫和包利興提出的滲透力定義式(36),在Terzaghi有效應力方程和一維問題范疇內是嚴謹成立的,仍不過時.

    致謝感謝匿名審稿專家給本文提出的寶貴指導意見;科學出版社工程技術分社張曉娟編輯和北京交通大學土建學院李艦博士協(xié)助進行本文標點符號和個別語言的修改完善;北京交通大學陳曦教授對摘要的英文翻譯加以潤色和改進;研究生王思遠將圖1和圖2改繪為AUTOCAD格式.謹此致謝.

    1弗洛林BA.土力學原理(第1卷).同濟大學土力學及地基基礎教研室,譯.北京:中國工業(yè)出版社,1965(Φлoрин BA.Fundamentals of Sol Mechanics(Vol.1).The Faculty of Soil Mechanics and Foundation Engineering of Tongji Univerisity,transl.Beijing:China Industry Press,1965(in Chinese))

    2格爾西萬諾夫HM,包利興ДE.土力學的理論原理及其實際應用.魏汝龍,譯.北京:建筑工程出版社出版,1958(Γeрсeвaнoв HM,Πoльшин ДE.Theoretical Principles of Soil Mechanics with Its Practical Applications.Wei Rulong,transl.Beijing:China Architecture&Building Press,1958(in Chinese))

    3 Terzaghi,K.Theoretical Soil Mechanics.New York:John Wiley and Sons,Inc.,1943

    4李廣信.論土骨架與滲透力.巖土工程學報,2016,38(8):1522-1528(Li Guangxin.On soil skeleton and seepage force.Chinese Journal of Geotechnical Engineering,2016,38(8):1522-1528(in Chinese))

    5祝恩陽.考慮孔隙水浮力的滲透力公式推導.力學與實踐,2015,37(6):761-762(Zhu Enyang.Derivation of the seepage force.Mechanics in Engineering,2015,37(6):761-762(in Chinese))

    6蔣中明,龍芳,熊小虎等.邊坡穩(wěn)定性分析中的滲透力計算方法考證.巖土力學,2015,36(9):2478-2486,2493(Jiang Zhongming,Long Fang,Xiong Xiaohu,et al.Study of calculation methods of acting force of seepage in slope stability analysis.Rock and Soil Mechanics,2015,36(9):2478-2486,2493(in Chinese))

    7邵龍?zhí)?郭曉霞.有效應力新解.北京:中國水利水電出版社,2014(Shao Longtan,Guo Xiaoxia.Interpretation of E ff ective Stress.Beijing:China Water&Power Press,2014(in Chinese))

    8錢家歡.土力學.南京:河海大學出版社,1988(Qian Jiaohuan.Soil Mechanics.Nanjing:Hehai University Press,1988(in Chinese))

    9陳仲頤,周景星,王洪瑾.土力學.北京:清華大學出版社,1994(Chen Zhongyi,Zhou Jingxing,Wang Hongjin.Soil Mechanics.Beijing:Tsinghua University Press,1994(in Chinese))

    10丁洲祥,龔曉南,謝永利.歐拉描述的大變形固結理論.力學學報,2005,37(1):92-99(Ding Zhouxiang,Gong Xiaonan,Xie Yongli.Finite strain consolidation theory with eulerian description.Acta Mechanica Sinica,2005,37(1):92-99(in Chinese))

    11 Timoshenko S,Goodier JN.Theory of Elasticity. New York:McGraw-Hill Book Company,Inc.,1951

    12 Gerhart PM,Gerhart AL,Hochstein JI.Munson,Young and Okiishi’s Fundamentals of Fluid Mechanics,8th ed.New York:John Wiley&Sons,2016

    13 Budhu M.Soil Mechanics Fundamentals(Metric Version).Chichester,West Sussex,U.K.:John Wiley&Sons,Ltd.,2015

    14 Das BM.Advanced Soil Mechanics,4th ed.Boca Ration:CRC Press,2013

    15 Lambe TW,Whitman RV.Soil Mechanics.New York:John Wiley&Sons,1969

    16趙成剛,白冰.土力學原理(修訂本).北京:清華大學出版社,北京交通大學出版社,2009(Zhao Chenggang,Bai Bing.Principles of Soil Mechanics,Modi fi ed Edition.Beijing:Tsinghua University Press,2009(in Chinese))

    17 Li X,Ming H.Seepage driving e ff ect on deformations of San Fernando dams.Soil Dynamics and Earthquake Engineering,2004,24:979-992

    18陳愈炯.有效應力原理在中國的應用實例.巖土工程學報,2015,37(9):1674-1677(Chen Yujong.Examples of application of effective stress principle in China.Chinese Journal of Geotechnical Engineering,2015,37(9):1674-1677(in Chinese))

    19 Terzaghi,K,Peck RB,Mesri G.Soil Mechanics in Engineering Practice,3rd ed.New York:John Wiley&Sons,Inc.,1996

    20 Biot,MA.Theory of elasticity and consolidation for a porous anisotropic solid.Journal of Applied Physics,1955,26(2):182-185

    21《工程地質手冊》編委會.工程地質手冊,第4版.北京:中國建筑工業(yè)出版社,2007(Editorial Board of Handbook of Engineering Geology.Handbook of Engineering Geology,4th ed.Beijing:China Architecture&Building Press,2007(in Chinese))

    22張邦芾.基坑工程地下水滲流場特性研究.[碩士論文].北京:中國建筑科學研究院,2014(Zhang Bangfu.Study on the characteristics of groundwater seepage fi eld in foundation pit engineering.[Master Thesis].Beijing:China Academy of Building Research,2014(in Chinese))

    23丁洲祥.滲透力概念的力學分析及廣義化探討.巖土工程學報,2017,待刊 (Ding Zhouxiang.Mechanical fundamentals of seepage force concept and its generalization.Chinese Journal of Geotechnical Engineering,2017,in Press(in Chinese))

    24 Lei GH,Chiu ACF,Ng CWW.Is Terzaghi’s e ff ective stress a stress variable under seepage conditions.Journal of Central South University,2015,22:2316-2321

    25雷國輝,鄭強.瑞典條分法剖析引發(fā)的有效應力和滲流力概念問題.巖土工程學報,2012,34(4):667-676(Lei Guohui,Zheng Qiang.Issues on concepts of e ff ective stress and seepage force arising from anatomizing Swedish slice method.Chinese Journal of Geotechnical Engineering,2012,34(4):667-676(in Chinese))

    26 Pan Q,Xu J,Dias D.Three-dimensional stability of a slope subjected to seepage force.International Journal of Geomechanics,2017,17(8):04017035,1-9

    27 Al-Madhhachi AT,Fox GA,Hanson GJ,et al.Mechanistic detachment rate model to predict soil erodibility due to fl uvial and seepage force.Journal of Hydraulic Engineering,2014,140(5):04014010,1-12

    ONE-DIMENSIONAL SEEPAGE FORCE AND BUOYANCY1)

    Ding Zhouxiang?,?,2)?(School of Civil Engineering,Beijing Jiaotong University,Beijing 100044,China)
    ?(Shaanxi Provincial Major Laboratory for Highway Bridge&Tunnel,Xi’an 710064,China)

    Forthesakeofsimplicity,twomethodsofderivingseepageforceandbuoyancyassociatedwithone-dimensional fl uid fl ow in saturated soil were developed to give a uselful supplement to conventional macro-scale method of isolated pore fl uid.To study the seepage force and buoyancy,three fundamental governing equations,namely,equilibrium di ff erential equation in elasticity,Terzaghi’s e ff ective stress equation and simpli fi ed Bernoulli’s equation in fl uid mechanics,are utilized for the derivation of seepage force and buoyancy.Based on the fundamental governing equations,it is readily to derive the di ff erential equations of equilibrium corresponding to the isolated soil skeleton and the isolated pore fl uid,and hence,the connotation and nature of seepage force and buoyancy can be disclosed.The seepage force per unit volume of saturated soil arises from the gradient of fl uid pressure in terms of total head.The buoyancy originates in the vertical component of gradient of fl uid pressure in terms of position head.Both the seepage force and the buoyancy can be embodied in equilibirium di ff erential equation with respect to skeleton or pore fl uid of satured soil.In geotechnical practice,the three governing equations can be directly applied to the calculation of e ff ective stress.Only in some simpli fi ed cases,the distribution of e ff ective stress can be obtained using the seepage force and buoyancy.Additionally,some hot topics are also discussed in this study,and emphasis shall be put on the discussion on the rationality and the potential application risk of the new de fi nition of seepage force,i.e.,j=nγwi.The strictness of one classic formulation of seepage force in history was thoroughly validated when considering the derivative of seepage velocity with respect to time.It is noteworthy that in the study of seepage force and buoyance from the perspective of soil mechanics,the Terzaghi’s e ff ective stress equation shall be reasonably implemented with practical signi fi cance.

    seepage force,buoyancy,e ff ective stress principle,pore water pressure,di ff erential equilibrium equation.

    TU431

    A

    10.6052/0459-1879-17-001

    2017–01–01收稿,2017–06–08 錄用,2017–06–22 網(wǎng)絡版發(fā)表.

    1)國家自然科學基金(51278028,41172221)和中央高?;究蒲袠I(yè)務費專項資金(310821161119,2014JBM087)資助項目.

    2)丁洲祥,博士,副教授,主要研究方向:固結理論及其在巖土和隧道工程中的應用.E-mail:dingzhouxiang@163.com

    丁洲祥.一維滲透力與浮力.力學學報,2017,49(5):1154-1162

    Ding Zhouxiang.One-dimensional seepage force and bouyancy.Chinese Journal of Theoretical and Applied Mechanics,2017,49(5):1154-1162

    猜你喜歡
    滲透力土力學浮力
    “浮力”知識鞏固
    我們一起來“制服”浮力
    浮力大小由誰定
    垃圾土力學特性及其路基填筑應用研究
    土力學課程思政元素的探索與實踐
    土力學中土的密度和重度概念的分析與討論
    中華建設(2020年5期)2020-07-24 08:55:42
    提升高校思政教育“滲透力”
    當代貴州(2017年1期)2017-03-06 19:03:13
    神奇的浮力
    以企業(yè)發(fā)展為基礎如何做好職工思想政治教育工作
    夏惠明教授推拿治療腰椎間盤突出癥急性期經驗介紹
    精品久久久久久久久av| 亚洲av日韩在线播放| 99热全是精品| 精品一区二区免费观看| 色5月婷婷丁香| 麻豆乱淫一区二区| 狂野欧美激情性xxxx在线观看| 久久久久免费精品人妻一区二区| 久99久视频精品免费| 日日撸夜夜添| 色综合色国产| 51国产日韩欧美| 国产男人的电影天堂91| 白带黄色成豆腐渣| 最近2019中文字幕mv第一页| 高清av免费在线| 国产有黄有色有爽视频| 麻豆乱淫一区二区| av一本久久久久| 国产精品无大码| 啦啦啦韩国在线观看视频| 99久久精品国产国产毛片| 精品久久久久久电影网| 精品熟女少妇av免费看| 日韩欧美国产在线观看| 69av精品久久久久久| 精品国产一区二区三区久久久樱花 | 建设人人有责人人尽责人人享有的 | 日本一本二区三区精品| 日韩欧美精品免费久久| 国产人妻一区二区三区在| 国产精品国产三级国产专区5o| 18禁动态无遮挡网站| 两个人视频免费观看高清| 男女啪啪激烈高潮av片| 成人一区二区视频在线观看| 国产片特级美女逼逼视频| 亚洲内射少妇av| 国产精品久久久久久久电影| 日韩伦理黄色片| 精品久久久久久久末码| 国产精品福利在线免费观看| 成人漫画全彩无遮挡| 三级经典国产精品| 精品久久久久久久久久久久久| 午夜福利视频1000在线观看| 一边亲一边摸免费视频| 国产av在哪里看| 国产有黄有色有爽视频| 少妇的逼水好多| 国产精品久久久久久久久免| 在线观看人妻少妇| 综合色丁香网| 日韩大片免费观看网站| 免费高清在线观看视频在线观看| 国内精品一区二区在线观看| 欧美97在线视频| 国产单亲对白刺激| 午夜免费激情av| 精品一区在线观看国产| 欧美日韩精品成人综合77777| 波野结衣二区三区在线| 日本一本二区三区精品| 人妻系列 视频| 久久97久久精品| 国产 亚洲一区二区三区 | 天堂网av新在线| 免费黄频网站在线观看国产| 国产在线一区二区三区精| 韩国高清视频一区二区三区| 好男人视频免费观看在线| 少妇裸体淫交视频免费看高清| 成人特级av手机在线观看| 爱豆传媒免费全集在线观看| 婷婷六月久久综合丁香| av播播在线观看一区| 久久热精品热| 波野结衣二区三区在线| 伊人久久精品亚洲午夜| 欧美高清成人免费视频www| 日韩欧美国产在线观看| 极品教师在线视频| 国产精品一区www在线观看| 国语对白做爰xxxⅹ性视频网站| 综合色丁香网| 精品一区在线观看国产| 国产亚洲午夜精品一区二区久久 | av播播在线观看一区| 亚洲激情五月婷婷啪啪| 人妻夜夜爽99麻豆av| 国产有黄有色有爽视频| 精品不卡国产一区二区三区| 九草在线视频观看| 亚洲性久久影院| 成人午夜精彩视频在线观看| 免费黄频网站在线观看国产| 老师上课跳d突然被开到最大视频| 久久久久久久大尺度免费视频| 最近的中文字幕免费完整| 国产精品久久久久久av不卡| 国精品久久久久久国模美| or卡值多少钱| 男女国产视频网站| 九色成人免费人妻av| 亚洲精品乱久久久久久| 国产精品久久久久久精品电影| 日韩伦理黄色片| 亚洲成人精品中文字幕电影| 日本与韩国留学比较| 久久精品人妻少妇| 欧美精品一区二区大全| 国产午夜精品一二区理论片| 直男gayav资源| 久久精品国产鲁丝片午夜精品| 99热全是精品| 久久午夜福利片| 男女国产视频网站| 精品一区二区免费观看| 久久久久久伊人网av| 国产 亚洲一区二区三区 | 久久精品国产自在天天线| 国产 一区精品| 亚洲成人精品中文字幕电影| 亚洲欧美日韩无卡精品| 97人妻精品一区二区三区麻豆| 亚洲精品成人久久久久久| 自拍偷自拍亚洲精品老妇| 亚洲精品,欧美精品| 国产精品福利在线免费观看| freevideosex欧美| 日韩一本色道免费dvd| 在线观看av片永久免费下载| 日韩一本色道免费dvd| 欧美精品一区二区大全| 91久久精品电影网| 边亲边吃奶的免费视频| 欧美区成人在线视频| 大陆偷拍与自拍| 老司机影院毛片| 一个人看的www免费观看视频| 99久久九九国产精品国产免费| 欧美人与善性xxx| 十八禁网站网址无遮挡 | 久久久久九九精品影院| 97热精品久久久久久| 性色avwww在线观看| 在线观看免费高清a一片| 日本av手机在线免费观看| 成人午夜高清在线视频| 一级二级三级毛片免费看| 亚洲自偷自拍三级| 男女下面进入的视频免费午夜| 亚洲精品自拍成人| 亚洲精品自拍成人| 男的添女的下面高潮视频| 麻豆久久精品国产亚洲av| 1000部很黄的大片| 身体一侧抽搐| 美女脱内裤让男人舔精品视频| 日韩欧美一区视频在线观看 | 国产成人免费观看mmmm| 国产精品综合久久久久久久免费| 大陆偷拍与自拍| 亚洲最大成人中文| 秋霞伦理黄片| 十八禁国产超污无遮挡网站| 熟妇人妻不卡中文字幕| 人妻少妇偷人精品九色| 性插视频无遮挡在线免费观看| 青春草亚洲视频在线观看| 国产视频首页在线观看| 99久久精品国产国产毛片| 亚洲人成网站在线播| 亚洲人成网站高清观看| 十八禁国产超污无遮挡网站| 国产v大片淫在线免费观看| 一本一本综合久久| av国产久精品久网站免费入址| videossex国产| 亚洲av免费高清在线观看| 91精品伊人久久大香线蕉| 久久久精品94久久精品| 777米奇影视久久| 久久国产乱子免费精品| 大香蕉久久网| 亚洲在久久综合| 精品少妇黑人巨大在线播放| 成人综合一区亚洲| 老女人水多毛片| 午夜久久久久精精品| 成人av在线播放网站| 91久久精品国产一区二区成人| 成人特级av手机在线观看| 日本免费a在线| 男女国产视频网站| 欧美一区二区亚洲| av在线天堂中文字幕| av卡一久久| 国产亚洲av嫩草精品影院| 欧美xxxx黑人xx丫x性爽| 男女视频在线观看网站免费| 亚洲va在线va天堂va国产| 国产亚洲精品久久久com| 欧美xxxx性猛交bbbb| 国产麻豆成人av免费视频| 久99久视频精品免费| 嘟嘟电影网在线观看| 乱码一卡2卡4卡精品| 极品少妇高潮喷水抽搐| 99久久人妻综合| 最近最新中文字幕大全电影3| 免费黄网站久久成人精品| 国产综合懂色| 精品久久久噜噜| 午夜福利高清视频| 网址你懂的国产日韩在线| 你懂的网址亚洲精品在线观看| 亚洲精品一二三| 日本熟妇午夜| 日本免费在线观看一区| 亚洲精品aⅴ在线观看| 三级国产精品欧美在线观看| 亚洲精品亚洲一区二区| 国产精品无大码| 国产精品国产三级国产av玫瑰| 亚洲av成人av| 久久久精品欧美日韩精品| 一区二区三区乱码不卡18| 狂野欧美白嫩少妇大欣赏| 亚洲成人一二三区av| 观看免费一级毛片| 80岁老熟妇乱子伦牲交| 亚洲色图av天堂| 纵有疾风起免费观看全集完整版 | 色播亚洲综合网| 在线免费观看的www视频| 国国产精品蜜臀av免费| 偷拍熟女少妇极品色| 国产精品久久久久久av不卡| 亚洲欧美中文字幕日韩二区| 国产亚洲一区二区精品| 亚洲图色成人| 中文字幕av在线有码专区| 免费看光身美女| 亚洲欧美精品自产自拍| 午夜福利高清视频| 亚洲电影在线观看av| 国产高清有码在线观看视频| 日日摸夜夜添夜夜添av毛片| 免费看av在线观看网站| 亚洲欧洲国产日韩| 亚洲成人一二三区av| 精品人妻熟女av久视频| 久久久色成人| 国产69精品久久久久777片| 99热6这里只有精品| 女的被弄到高潮叫床怎么办| 可以在线观看毛片的网站| 精品酒店卫生间| 精品人妻偷拍中文字幕| 久久精品国产自在天天线| 久久久久久国产a免费观看| 天天躁日日操中文字幕| 国产黄色视频一区二区在线观看| 精品久久久久久成人av| 国产日韩欧美在线精品| 麻豆成人av视频| 欧美三级亚洲精品| 精品人妻视频免费看| 高清日韩中文字幕在线| 日韩强制内射视频| 国产精品久久久久久av不卡| 亚洲成人久久爱视频| 欧美 日韩 精品 国产| av黄色大香蕉| av在线亚洲专区| 免费观看在线日韩| 亚洲熟女精品中文字幕| 丰满人妻一区二区三区视频av| 激情五月婷婷亚洲| 亚洲婷婷狠狠爱综合网| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩视频高清一区二区三区二| videossex国产| 18禁动态无遮挡网站| 99久久精品国产国产毛片| 国产69精品久久久久777片| 国产精品嫩草影院av在线观看| 少妇熟女欧美另类| 精品久久久噜噜| av天堂中文字幕网| 蜜臀久久99精品久久宅男| 97超视频在线观看视频| 偷拍熟女少妇极品色| 最近2019中文字幕mv第一页| 午夜精品一区二区三区免费看| 国产精品.久久久| 2022亚洲国产成人精品| 国产精品一及| 丝袜美腿在线中文| 久久久久性生活片| 国产在线男女| 亚洲国产高清在线一区二区三| 国产一区二区在线观看日韩| 听说在线观看完整版免费高清| 国产精品久久久久久精品电影小说 | 国产久久久一区二区三区| 日韩一区二区视频免费看| 国产精品一区二区三区四区免费观看| 婷婷色麻豆天堂久久| 日韩三级伦理在线观看| 免费观看的影片在线观看| 午夜精品国产一区二区电影 | 一区二区三区四区激情视频| 免费观看的影片在线观看| 亚洲美女搞黄在线观看| 亚洲无线观看免费| av.在线天堂| 国产极品天堂在线| 国产精品久久久久久精品电影| 久久99热这里只频精品6学生| 国产美女午夜福利| 国产综合懂色| 国产白丝娇喘喷水9色精品| 久久99蜜桃精品久久| 一边亲一边摸免费视频| av国产久精品久网站免费入址| 超碰97精品在线观看| 日本一本二区三区精品| 国产老妇女一区| 日本午夜av视频| 精品一区在线观看国产| 高清欧美精品videossex| 干丝袜人妻中文字幕| 人人妻人人澡人人爽人人夜夜 | 久久99蜜桃精品久久| 精品久久国产蜜桃| 日韩欧美精品免费久久| 中文字幕免费在线视频6| 精品少妇黑人巨大在线播放| 久久精品人妻少妇| 久久综合国产亚洲精品| 少妇熟女欧美另类| 美女被艹到高潮喷水动态| 国产男女超爽视频在线观看| 日韩中字成人| av在线天堂中文字幕| 人妻少妇偷人精品九色| 99热这里只有是精品50| 免费大片黄手机在线观看| 国产精品国产三级专区第一集| 日韩强制内射视频| 天堂影院成人在线观看| 中文字幕av成人在线电影| 久久久精品欧美日韩精品| 国产69精品久久久久777片| 成人av在线播放网站| 久久精品国产亚洲av涩爱| 国产成人91sexporn| 成年免费大片在线观看| 国产伦在线观看视频一区| 美女黄网站色视频| 青春草国产在线视频| 亚洲第一区二区三区不卡| 国产精品av视频在线免费观看| 欧美zozozo另类| 成人午夜精彩视频在线观看| 一级毛片 在线播放| 免费看av在线观看网站| 亚洲精品亚洲一区二区| 成人无遮挡网站| 少妇熟女aⅴ在线视频| 久久精品国产自在天天线| 色播亚洲综合网| 午夜免费男女啪啪视频观看| av免费在线看不卡| 亚洲av免费高清在线观看| 日韩,欧美,国产一区二区三区| 高清av免费在线| 国产亚洲av嫩草精品影院| 91久久精品国产一区二区成人| 国产精品美女特级片免费视频播放器| 亚洲人成网站在线播| 亚洲色图av天堂| 国产不卡一卡二| 一级二级三级毛片免费看| 亚洲一区高清亚洲精品| 久久久久免费精品人妻一区二区| 日韩精品青青久久久久久| 五月玫瑰六月丁香| 国产精品美女特级片免费视频播放器| 伊人久久国产一区二区| 国产男女超爽视频在线观看| 七月丁香在线播放| 亚洲综合精品二区| 亚洲人成网站在线播| 三级男女做爰猛烈吃奶摸视频| 日韩三级伦理在线观看| 男人爽女人下面视频在线观看| 亚洲av中文字字幕乱码综合| 汤姆久久久久久久影院中文字幕 | 久久久久久久午夜电影| 亚洲四区av| 国产精品不卡视频一区二区| 男女边摸边吃奶| 综合色av麻豆| 久久久久久久大尺度免费视频| 三级经典国产精品| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品影视一区二区三区av| 精品人妻一区二区三区麻豆| 大话2 男鬼变身卡| 国产视频首页在线观看| 久久精品久久久久久久性| 精品少妇黑人巨大在线播放| 亚洲欧美成人综合另类久久久| 大香蕉97超碰在线| 看黄色毛片网站| 婷婷色综合大香蕉| 欧美一区二区亚洲| 色尼玛亚洲综合影院| 亚洲欧美精品自产自拍| 91久久精品国产一区二区三区| 久久99热6这里只有精品| 国内精品一区二区在线观看| www.色视频.com| 成人av在线播放网站| 不卡视频在线观看欧美| 精品久久久久久电影网| av女优亚洲男人天堂| 天堂av国产一区二区熟女人妻| 只有这里有精品99| 一夜夜www| 免费少妇av软件| 国产 一区精品| 99久久精品国产国产毛片| 亚洲在线自拍视频| 亚洲精品乱码久久久久久按摩| 欧美zozozo另类| 免费观看精品视频网站| 日本与韩国留学比较| 18禁裸乳无遮挡免费网站照片| 联通29元200g的流量卡| 亚洲成人精品中文字幕电影| 91久久精品国产一区二区成人| 日本一本二区三区精品| 亚洲精品自拍成人| 丝袜美腿在线中文| 老师上课跳d突然被开到最大视频| 国产永久视频网站| 免费观看性生交大片5| 99热这里只有是精品50| 日产精品乱码卡一卡2卡三| 国产黄色小视频在线观看| 99热6这里只有精品| 黄色配什么色好看| videos熟女内射| 国产一级毛片七仙女欲春2| 亚洲最大成人av| 国产在线一区二区三区精| 国产乱人视频| 91精品伊人久久大香线蕉| 美女被艹到高潮喷水动态| 看非洲黑人一级黄片| 99久国产av精品国产电影| 大香蕉久久网| 一级毛片 在线播放| 日本黄色片子视频| 精品熟女少妇av免费看| 国产黄色小视频在线观看| 99热这里只有是精品50| 精品人妻视频免费看| 男人和女人高潮做爰伦理| 日本熟妇午夜| 又粗又硬又长又爽又黄的视频| 中文字幕制服av| 国产爱豆传媒在线观看| 日日摸夜夜添夜夜爱| 精品久久久久久久人妻蜜臀av| 日韩不卡一区二区三区视频在线| 久99久视频精品免费| .国产精品久久| 中国国产av一级| 乱码一卡2卡4卡精品| av国产免费在线观看| 国产麻豆成人av免费视频| 老司机影院毛片| 在线观看美女被高潮喷水网站| 亚洲人成网站在线观看播放| 少妇丰满av| 一级av片app| 两个人的视频大全免费| 欧美zozozo另类| 18+在线观看网站| 国产精品人妻久久久影院| 中文字幕亚洲精品专区| 永久免费av网站大全| 卡戴珊不雅视频在线播放| 亚洲国产日韩欧美精品在线观看| 亚洲av男天堂| 国产精品久久久久久久久免| 一区二区三区高清视频在线| 欧美xxxx黑人xx丫x性爽| 国产精品99久久久久久久久| 2018国产大陆天天弄谢| 午夜精品国产一区二区电影 | 又爽又黄a免费视频| 少妇人妻一区二区三区视频| 久久人人爽人人片av| 伦理电影大哥的女人| 国产视频首页在线观看| 又黄又爽又刺激的免费视频.| 看非洲黑人一级黄片| 日韩一区二区三区影片| 国产午夜福利久久久久久| 大片免费播放器 马上看| 色5月婷婷丁香| 日韩av免费高清视频| av黄色大香蕉| 日本色播在线视频| 国产中年淑女户外野战色| 我的老师免费观看完整版| 嫩草影院入口| 久久久久久久久久成人| 综合色av麻豆| 久久精品熟女亚洲av麻豆精品 | 一区二区三区四区激情视频| 国产又色又爽无遮挡免| 色综合亚洲欧美另类图片| 99久国产av精品| 国产高潮美女av| 春色校园在线视频观看| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 一级毛片aaaaaa免费看小| 免费av毛片视频| 日本黄色片子视频| 国产精品综合久久久久久久免费| 2021少妇久久久久久久久久久| 国产黄片美女视频| 国产黄色小视频在线观看| 免费看不卡的av| 国产国拍精品亚洲av在线观看| 国产乱人偷精品视频| 久久精品国产自在天天线| 久久精品久久久久久久性| 日韩精品青青久久久久久| 精品欧美国产一区二区三| 国产精品一区二区在线观看99 | 免费黄色在线免费观看| 成人毛片a级毛片在线播放| 美女cb高潮喷水在线观看| 国产精品.久久久| 国产av国产精品国产| 久久这里有精品视频免费| 国产一区二区三区av在线| 18禁裸乳无遮挡免费网站照片| 18+在线观看网站| 亚洲国产精品专区欧美| 国产片特级美女逼逼视频| 天堂影院成人在线观看| av在线蜜桃| 欧美zozozo另类| 99久久人妻综合| 午夜福利在线观看免费完整高清在| 中文精品一卡2卡3卡4更新| 水蜜桃什么品种好| 午夜福利网站1000一区二区三区| 好男人在线观看高清免费视频| 你懂的网址亚洲精品在线观看| 欧美成人一区二区免费高清观看| 亚洲色图av天堂| 69av精品久久久久久| www.av在线官网国产| 日日摸夜夜添夜夜爱| 少妇高潮的动态图| 日本av手机在线免费观看| 日日摸夜夜添夜夜添av毛片| 成人漫画全彩无遮挡| 免费看不卡的av| a级毛片免费高清观看在线播放| 国产伦精品一区二区三区四那| 2018国产大陆天天弄谢| 啦啦啦韩国在线观看视频| 亚洲自偷自拍三级| 男女那种视频在线观看| 日韩欧美三级三区| 天堂影院成人在线观看| 一级毛片 在线播放| 国产 一区 欧美 日韩| 亚洲精品乱久久久久久| av在线老鸭窝| 中文精品一卡2卡3卡4更新| 午夜精品一区二区三区免费看| 午夜激情久久久久久久| 欧美成人精品欧美一级黄| 国产不卡一卡二| 欧美成人午夜免费资源| 亚洲综合色惰| 国产亚洲av片在线观看秒播厂 | 男女国产视频网站| 久久久久久久久中文| 乱人视频在线观看| 亚洲无线观看免费| 免费观看性生交大片5| 色综合色国产| 丰满少妇做爰视频| 在线a可以看的网站| 精品久久久精品久久久| 亚洲欧美精品自产自拍| 中文乱码字字幕精品一区二区三区 | 亚洲精品日韩av片在线观看| 99热这里只有是精品在线观看| 亚洲欧美日韩无卡精品| 精品久久久久久成人av|