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

    兩種快速測(cè)定非飽和土 水力學(xué)參數(shù)方法的對(duì)比分析

    2012-01-08 07:12:30伊盼盼韋昌富魏厚振曹華峰徐文強(qiáng)
    巖土力學(xué) 2012年7期
    關(guān)鍵詞:水力學(xué)非飽和吸力

    伊盼盼 ,韋昌富,魏厚振,曹華峰,徐文強(qiáng)

    (1.中國科學(xué)院武漢巖土力學(xué)研究所 巖土力學(xué)與工程國家重點(diǎn)實(shí)驗(yàn)室,武漢 430071;2.武昌理工學(xué)院,武漢 430223)

    1 引 言

    非飽和土水力學(xué)參數(shù)包括土-水特征曲線和滲透函數(shù)。土-水特征曲線(SWCC)表示非飽和土的基質(zhì)吸力與含水率(如飽和度或體積含水率)之間的關(guān)系。滲透函數(shù)是滲透系數(shù)與含水率之間的關(guān)系。這兩大水力學(xué)參數(shù)在分析非飽和土變形與強(qiáng)度、降雨入滲與邊坡穩(wěn)定、污染物地下遷移等問題中扮演著重要角色[1-2]。

    傳統(tǒng)的非飽和土水力學(xué)參數(shù)的測(cè)定方法均是在試樣處于平衡態(tài)下測(cè)得的,耗時(shí)較長(zhǎng)。另外,由于加載時(shí)間較長(zhǎng),高壓氣體會(huì)通過水的流動(dòng)在陶土板背面析出,影響測(cè)試精度。因此,如何能夠快速、有效地確定非飽和土水力學(xué)參數(shù)成為需要解決的關(guān)鍵問題。

    通過流動(dòng)試驗(yàn)并利用數(shù)值反算方法對(duì)結(jié)果進(jìn)行分析能夠快速得到非飽和土的水力學(xué)參數(shù)。Gardner[3]首次提出流動(dòng)試驗(yàn)法,該方法通過求解非線性流動(dòng)方程得到非飽和土水力學(xué)參數(shù)。Kool 等[4]提出一步流動(dòng)試驗(yàn)方法,結(jié)合數(shù)值反算方法得到非飽和土水力學(xué)參數(shù)。盡管一步流動(dòng)試驗(yàn)?zāi)軐?duì)試驗(yàn)結(jié)果進(jìn)行很好地?cái)M合,但容易出現(xiàn)參數(shù)不惟一的現(xiàn)象。van Dam[5-6]等指出若能同時(shí)測(cè)得試樣中的孔隙水壓力,則可以降低參數(shù)不惟一的問題。另外,Eching等[7-9]同時(shí)開展了一步流動(dòng)和多步流動(dòng)試驗(yàn),結(jié)果顯示多步流動(dòng)試驗(yàn)?zāi)芨玫慕鉀Q參數(shù)不惟一的問題。

    文中提出能夠快速測(cè)定非飽和土水力學(xué)參數(shù)的一步流動(dòng)和動(dòng)態(tài)多步流動(dòng)方法,并對(duì)這兩種方法進(jìn)行了對(duì)比分析。其中,一步流動(dòng)方法是利用HYDRUS-1D 一維水分運(yùn)移模型擬合試樣一步脫濕流動(dòng)試驗(yàn)下溢出水量隨時(shí)間變化的關(guān)系曲線,進(jìn)而反算得出非飽和土的土-水特征曲線和滲透函數(shù)。該方法能夠降低參數(shù)的不惟一性,同時(shí)考慮了陶土板材料參數(shù)對(duì)試驗(yàn)溢出水隨時(shí)間變化規(guī)律的影響。動(dòng)態(tài)多步流動(dòng)方法與上述多步流動(dòng)試驗(yàn)不同,它是基于Wei 和Muraleetharan[10]、Wei 和Dewoolkar[11]提出的能夠描述靜態(tài)土-水特征曲線和動(dòng)態(tài)土-水特征曲線間關(guān)系的動(dòng)態(tài)模型,開展多步流動(dòng)試驗(yàn)測(cè)定非飽和土的水力學(xué)參數(shù)。在此動(dòng)態(tài)多步流動(dòng)試驗(yàn)過程中不需要土樣處于水力學(xué)平衡狀態(tài),因而節(jié)省了試驗(yàn)所需時(shí)間。該試驗(yàn)中試樣在每一吸力步下都達(dá)到平衡。

    2 兩種方法的原理

    2.1 動(dòng)態(tài)多步流動(dòng)方法

    2.1.1 動(dòng)態(tài)(非平衡態(tài))土-水特征關(guān)系

    考慮一圓柱形非飽和土樣,其兩端分別與氣源和水源相連接。在初始時(shí),土樣處于靜力平衡狀態(tài),其飽和度為,吸力為=,這里( Sr)是平衡態(tài)的土-水特征函數(shù)。此時(shí),若改變?cè)嚇觾啥说膲毫Σ?,試樣中的水?huì)從一端流向另外一端,故試樣的飽和度也會(huì)隨之發(fā)生變化。Wei 和 Muralee- tharan[10]與Wei 和Dewoolkar[11]提出基質(zhì)吸力的變化與飽和度之間存在以下關(guān)系:

    式中:基質(zhì)吸力 pc= pa- pw,這里 pa為孔隙氣壓力, pw為孔隙水壓力;μ 為材料參數(shù);θ 為體積含水率。

    如果不考慮試樣的體積變化,式(1)可以轉(zhuǎn)化為

    式中: nμ μ′= ,n 為孔隙率。

    若吸力增加一小量 cpδ ,則飽和度的變化 rSδ與 cpδ 的關(guān)系為

    式(3)類似于線性、黏彈性固體本構(gòu)關(guān)系Kelvin模型,可以用圖1 來表示。由此可見,當(dāng)吸力改變時(shí),飽和度不是突然變化的,而是隨時(shí)間慢慢地改變,這一過程與實(shí)際是相吻合的。

    圖1 線性黏彈性Kelvin模型 Fig.1 Linear viscoelastic model of Kelvin

    假設(shè)擾動(dòng)(吸力)是突然增加的,故有

    式中: ( )H t 為Heaviside 函數(shù)。

    那么求解式(3),得

    如果吸力是多步施加的,即

    同理,求解式(3)可以得到多步流動(dòng)下試樣飽和度隨時(shí)間的演化方程為

    式中:kC 和分別為吸力為時(shí)的穩(wěn)態(tài)土-水特征曲線斜率(容水率)和擾動(dòng)后孔隙水流動(dòng)的特征時(shí)間。

    2.1.2 平衡條件下土-水特征曲線和滲透系數(shù)的確定

    由式(7)可以看出,式中主要由C 和cτ 兩大參數(shù)。如果這兩個(gè)參數(shù)值已知,那么試樣的飽和度隨時(shí)間的演化曲線也就明確了,因此,確定這兩大參數(shù)值非常重要。這里列出兩種求解的方法。

    方法1 為兩點(diǎn)法,考慮某時(shí)間段 ti-1< t < ti中的兩個(gè)時(shí)刻和,由式(7)可以得到一個(gè)非線性方程組,采用迭代法可以求出每級(jí)吸力步對(duì)應(yīng)的C 和τc值。

    方法2 為最小二乘法。根據(jù)實(shí)測(cè)的飽和度隨時(shí)間的演化曲線,在上面取足夠多具有代表性的點(diǎn),采用最小二乘法對(duì)實(shí)測(cè)數(shù)據(jù)進(jìn)行擬合,求出每級(jí)吸力步對(duì)應(yīng)的C 和cτ 值。

    確定了各級(jí)吸力作用下的容水率C 和特征時(shí)間τc,就可以得到在基質(zhì)吸力作用下試樣處于平衡態(tài)時(shí)對(duì)應(yīng)的飽和度:

    對(duì)于一定尺寸的土樣,孔隙水壓力耗散時(shí)間(特征時(shí)間)cτ 與土樣的滲透性有關(guān),通過對(duì)孔隙水的滲流過程分析,可以建立特征時(shí)間cτ 與滲透系數(shù)的關(guān)系。其過程如圖2 所示。由圖可知,初始時(shí)土樣處于平衡狀態(tài),在某個(gè)時(shí)刻,土樣底部的孔隙水壓力降低了 wpδ ,而其他條件不變。為了描述這一過程,假定在整個(gè)過程中孔隙氣體壓力ap 保持不變,同時(shí)假定土顆粒和水是不可壓縮的(即質(zhì)量密度sρ 和wρ 為常數(shù)),且不考慮變形的影響(孔隙率為常數(shù))。

    圖2 土中Darcy流動(dòng)示意圖 Fig.2 Darcy flow diagram of the soil

    在這些前提下,孔隙水質(zhì)量平衡方程為

    式中:w

    n 為水相孔隙率。

    式(9)可變形為

    孔隙水的流動(dòng)由以下Darcy 類方程控制:

    利用土-水特征關(guān)系rc( )S p 及C 的定義,將式(11)代入式(10)可得:

    亦可以寫成:

    這就是在上述試驗(yàn)條件下土樣孔隙水壓力耗散過程的控制方程,式中vc 為孔隙水壓力耗散系數(shù)。

    設(shè)試樣的長(zhǎng)度為L(zhǎng),那么由式(13)得出Darcy流動(dòng)的特征時(shí)間cτ 為

    2.2 一步流動(dòng)方法

    一步流動(dòng)方法的原理為:通過對(duì)飽和試樣施加一步吸力,測(cè)得試樣溢出水量隨時(shí)間的變化曲線。利用一維水分運(yùn)移過程分析軟件HYDRUS-1D 對(duì)溢出水隨時(shí)間流量曲線進(jìn)行擬合,同時(shí)分析反算出一些難以實(shí)測(cè)的水力學(xué)參數(shù),進(jìn)而得到非飽和土-水特征曲線和滲透函數(shù)曲線。

    土樣中水分變化可由Richards 方程描述,即

    式中:h 為壓力水頭; ( )K h 為滲透函數(shù);t 為時(shí)間;z 為土樣深度。

    試樣的初邊值條件表示為

    對(duì)試樣施加一步吸力P 之前,先施加一小吸力0P (略大于試樣進(jìn)氣值),使試樣由飽和狀態(tài)變?yōu)榉秋柡蜖顟B(tài)。這里00/h P gρ= , 0z= 是試樣上部邊界,z L= 為試樣下部邊界, ( )h t 為試樣下部邊界的水頭。

    土-水特征曲線模型選用比較常用的van Genu- chten[12]模型: )

    將van Genuchten 模型與Mualem[13]模型相結(jié)合可得到滲透函數(shù)的表達(dá)式為

    式中:sθ 為飽和含水率;rθ 為殘余含水率;α 、β 、γ 為經(jīng)驗(yàn)參數(shù);sK 為飽和滲透系數(shù);eS 為有效含水率;l 為孔隙結(jié)構(gòu)參數(shù)。

    3 試驗(yàn)設(shè)備和試驗(yàn)方法

    3.1 試驗(yàn)設(shè)備

    為了開展前面提出的多步流動(dòng)試驗(yàn),采用聯(lián)合測(cè)試系統(tǒng),其示意圖如圖3 所示。

    圖3 聯(lián)合測(cè)試系統(tǒng)示意圖 Fig.3 Sketch of combined testing system

    該系統(tǒng)主要由壓力控制系統(tǒng)、壓力室、儲(chǔ)水容器、稱量系統(tǒng)、數(shù)據(jù)采集系統(tǒng)等主要部分構(gòu)成。

    壓力控制系統(tǒng)由高壓氮?dú)馄?、減壓閥、壓力調(diào)節(jié)泵等部分構(gòu)成。減壓閥把氮?dú)馄坷锏母邏簹怏w變成低壓氣體,起粗調(diào)作用。壓力調(diào)節(jié)泵能精確控制給壓力室輸入的氣壓力值,起到細(xì)調(diào)作用。

    壓力室主要由裝樣容器和上、下蓋板等3 部分組成,其中下蓋板裝有高進(jìn)氣值陶土板。

    儲(chǔ)水容器的用途有兩個(gè),一是頂部和真空泵相連,可以用來對(duì)試樣進(jìn)行飽和;二是用于對(duì)系統(tǒng)管路的飽和。

    稱量系統(tǒng)主要由天平及上方開有一小孔的盛水杯組成,開一小孔有利于防止水分揮發(fā)。

    數(shù)據(jù)采集系統(tǒng)是通過LabVIEW 軟件來完成數(shù)據(jù)采集功能,它有兩個(gè)數(shù)據(jù)通道分別與天平和壓力調(diào)節(jié)泵相連,可以采集試樣的排出水量和對(duì)試樣施加的吸力隨時(shí)間變化的數(shù)據(jù)。采集的時(shí)間間隔可以從幾秒到幾小時(shí)不等,可根據(jù)試驗(yàn)需要任意設(shè)置。

    3.2 試驗(yàn)方法

    試驗(yàn)所用土樣取自黃河三角洲,其基本物理特性指標(biāo)如表1 所示。

    表1 土樣基本物理性質(zhì) Table 1 Physical properties of the soil sample

    取該土樣以初始含水率為9%,孔隙率以42%為控制目標(biāo),制備一個(gè)直徑為5.4 cm,高為3 cm 的圓柱形試樣。隨后,將制備好試樣采用抽氣飽和法進(jìn)行飽和。先對(duì)飽和試樣開展多步流動(dòng)試驗(yàn),具體過程如下:

    在 t0時(shí)刻,對(duì)試樣施加第1 級(jí)吸力增量,然后在吸力保持不變的條件下,直至 t1時(shí)刻(此時(shí)試樣并沒有達(dá)到平衡狀態(tài))。接著,在 t1時(shí)刻,施加第2 級(jí)吸力增量,并保持吸力不變至 t2時(shí)刻。以此類推,施加第N 級(jí)吸力增量,直至 tn時(shí)刻。根據(jù)這組多步流動(dòng)試驗(yàn)數(shù)據(jù),利用前面的理論模型及求解辦法,可以推算出試樣脫濕條件下平衡態(tài)的土-水特征曲線。

    由于文中采用的試樣在試驗(yàn)過程中其體變可忽略不計(jì),且經(jīng)歷了脫濕、吸濕之后試樣的結(jié)構(gòu)基本趨于穩(wěn)定,故在多步流動(dòng)試驗(yàn)后把試樣重新飽和,然后對(duì)其開展一步流動(dòng)試驗(yàn),初始飽和度的影響可以忽略,其過程如下:對(duì)該飽和試樣施加一小吸力步(略大于試樣的進(jìn)氣值),使之由飽和狀態(tài)變?yōu)榉秋柡蜖顟B(tài)。當(dāng)流動(dòng)停滯(穩(wěn)定狀態(tài))時(shí),對(duì)試樣施加一較大吸力步,直至溢出水量穩(wěn)定。這樣就可以得到脫濕段溢出水量隨時(shí)間變化的關(guān)系曲線。利用一維水分運(yùn)移過程分析軟件HYDRUS-1D 分別對(duì)脫濕段和吸濕段溢出水隨時(shí)間的流量曲線進(jìn)行擬合,反算得出一些難以實(shí)測(cè)的水力學(xué)參數(shù),進(jìn)而得到非飽和土干、濕土-水特征曲線和滲透函數(shù)曲線。

    4 試驗(yàn)結(jié)果分析

    4.1 動(dòng)態(tài)多步流動(dòng)試驗(yàn)結(jié)果

    由于每級(jí)吸力下加載的時(shí)間不同,測(cè)得的試樣飽和度隨時(shí)間的演化關(guān)系也不同,本文分別進(jìn)行了兩組不同時(shí)間步的流動(dòng)試驗(yàn)。每組流動(dòng)試驗(yàn)施加的吸力步及各吸力步的加載時(shí)間如表2 所示。為了進(jìn)行對(duì)比分析,試樣在各吸力步下達(dá)到平衡的時(shí)間也在表2 中列出。

    由表可以看出,這兩組流動(dòng)試驗(yàn)在各吸力下的加載時(shí)間遠(yuǎn)小于試樣達(dá)到平衡所需的時(shí)間,第2 組流動(dòng)試驗(yàn)的加載時(shí)間比第1 組稍長(zhǎng)。另外,這兩組流動(dòng)試驗(yàn)在初級(jí)吸力下都達(dá)到了平衡,原因是試樣由飽和區(qū)向非飽和區(qū)過渡時(shí),土-水特征曲線的斜率即C 值變化較大,故先對(duì)試樣施加一小吸力(略大于進(jìn)氣值),待達(dá)到平衡時(shí),試樣處于非飽和狀態(tài),以該狀態(tài)為初始狀態(tài)開展多步流動(dòng)試驗(yàn)?zāi)苡行У臏p少試驗(yàn)誤差。

    這兩組流動(dòng)試驗(yàn)得到的試樣飽和度隨時(shí)間的演化曲線如圖4 中數(shù)據(jù)點(diǎn)所示,利用前面提出的動(dòng)態(tài)模型,并采用求解方法2,分別對(duì)這兩組試驗(yàn)曲線進(jìn)行最小二乘擬合,其擬合結(jié)果如圖4 中實(shí)線所示,最小二乘擬合的精度如圖5 所示。

    從以上圖4、5 可以看出,這兩組流動(dòng)試驗(yàn)的最小二乘擬合精度都比較高,進(jìn)而說明了由最小二乘擬合得出的各級(jí)吸力步下C、τ 值的可靠性。

    4.2 一步流動(dòng)試驗(yàn)結(jié)果

    先對(duì)飽和試樣施加一較小吸力1 kPa,有少量水溢出,這時(shí)試樣從飽和狀態(tài)進(jìn)入非飽和狀態(tài),然后增加吸力至80 kPa,進(jìn)行脫濕試驗(yàn)。待溢出水量恒定后,根據(jù)試樣初始飽和質(zhì)量和溢出水的質(zhì)量可以得到體積含水率隨時(shí)間變化的關(guān)系。

    根據(jù)該關(guān)系曲線、初邊值條件、陶土板材料參數(shù)運(yùn)用HYDRUS-1D 一維水分運(yùn)移模型進(jìn)行反算,擬合出試樣脫濕段溢出水量隨時(shí)間的演化曲線。實(shí)測(cè)數(shù)據(jù)曲線和擬合曲線如圖6 所示。

    圖4 兩組流動(dòng)試驗(yàn)中實(shí)測(cè) S r-t 曲線與擬合結(jié)果的對(duì)比 Fig.4 Comparison between measured and fitted curves of S r-t during the two outflow tests

    圖5 兩組流動(dòng)試驗(yàn)中最小二乘的擬合精度 Fig.5 Accuracy of least square fitting during the two outflow tests

    圖6 脫濕段θ -t 曲線模型預(yù)測(cè)值與實(shí)測(cè)值比較 Fig.6 Comparison between estimated and measured values of θ -t curves during drying

    從圖中可以看出,HYDRUS-1D 內(nèi)部模型能夠很好地?cái)M合實(shí)測(cè)數(shù)據(jù),其擬合參數(shù)如表3 所示。

    表3 脫濕條件下的擬合參數(shù) Table 3 Fitting parameters for drying

    4.3 試驗(yàn)結(jié)果對(duì)比分析

    對(duì)于多步流動(dòng)試驗(yàn),由于各級(jí)吸力步下的C、τ 值已知,根據(jù)前面的理論模型,可得到平衡狀態(tài)下試樣的土-水特征曲線和滲透函數(shù)曲線;對(duì)于一步流動(dòng)試驗(yàn),根據(jù)擬合參數(shù)可反算得出平衡狀態(tài)下試樣的土-水特征曲線和滲透函數(shù)曲線。由這兩種試驗(yàn)方法得到的試樣在平衡態(tài)下的土-水特征曲線和實(shí)測(cè)值的對(duì)比結(jié)果如圖7 所示。

    從圖可以看出,一步流動(dòng)試驗(yàn)得到的試樣平衡態(tài)土-水特征曲線與實(shí)測(cè)值有些偏差,而多步流動(dòng)試驗(yàn)結(jié)果則與實(shí)測(cè)值接近,第2 組流動(dòng)試驗(yàn)的加載時(shí)間較長(zhǎng),得到的土-水特征曲線與實(shí)測(cè)值基本吻合。

    由這兩種試驗(yàn)方法得到的試樣滲透函數(shù)曲線,如圖8 所示。

    圖7 兩種流動(dòng)試驗(yàn)及實(shí)測(cè)得到的SWCC對(duì)比 Fig.7 Comparison of SWCC between the measured and the two outflow tests

    圖8 兩種流動(dòng)試驗(yàn)得到的滲透函數(shù)曲線對(duì)比 Fig.8 Comparison curves of hydraulic permeability function between the two outflow tests

    由圖可以看出,這兩種試驗(yàn)方法得到的滲透函數(shù)接近。對(duì)比多步流動(dòng)結(jié)果,一步流動(dòng)試驗(yàn)反算得到的滲透系數(shù)在試樣體積含水率高時(shí)偏高,在體積含水率低時(shí)偏低。這是由于一步流動(dòng)方法是根據(jù)溢出水量隨時(shí)間的變化關(guān)系反算得到,受邊界條件和陶土板材料參數(shù)的制約。

    5 結(jié) 論

    傳統(tǒng)的測(cè)試非飽和土力學(xué)參數(shù)的方法,需要很長(zhǎng)時(shí)間,并且難以克服由于加載時(shí)間較長(zhǎng)高壓氣體在陶土板背面析出的問題。因此,快速測(cè)定非飽和水力學(xué)參數(shù)的方法研究顯得尤為重要,通過研究可得如下結(jié)論:

    (1)一步流動(dòng)方法和動(dòng)態(tài)多步流動(dòng)方法均能在短時(shí)間內(nèi)測(cè)定非飽和土的水力學(xué)參數(shù),克服了傳統(tǒng)測(cè)試方法耗時(shí)較長(zhǎng)的問題,并且有效地避免了高壓氣體通過水的流動(dòng)在陶土板背面析出的問題,提高了測(cè)定非飽和土水力學(xué)參數(shù)的速度和準(zhǔn)確性。

    (2)動(dòng)態(tài)多步流動(dòng)方法得到試樣平衡態(tài)的土-水特征曲線與實(shí)測(cè)值基本吻合,無需等到試樣達(dá)到平衡狀態(tài)后再進(jìn)行測(cè)定,大大縮短了測(cè)試時(shí)間;一步流動(dòng)方法僅對(duì)試樣施加一步吸力,比動(dòng)態(tài)多步流動(dòng)方法所用的時(shí)間短,但得出的試樣平衡態(tài)土-水特征曲線和實(shí)測(cè)數(shù)值有些偏差。因此,通過對(duì)比分析,動(dòng)態(tài)多步流動(dòng)方法的結(jié)果要優(yōu)于一步流動(dòng)試驗(yàn)方法。

    [1] FREDLUND D G, XING A Q, FREDLUND M D, et al. The relationship of the unsaturated soil shear strength function to the soil-water characteristic curve[J]. Cana- dian Geotechnical Journal, 1996, 32(3): 440-448.

    [2] VANAPALLI S K, FREDLUND D G, PUFAHL D E, et al. Model for the prediction of shear strength with respect to soil suction[J]. Canadian Geotechnical Journal, 1996, 33(3): 379-392.

    [3] GARDNER W R. Calculation of capillary conductivity from pressure plate outflow data[J]. Soil Science Society of America Journal, 1956, 20: 317-320.

    [4] KOOL J B, PARKER J C, VAN GENUCHTEN M TH. Determining soil hydraulic properties from one-step outflow experiments by parameter estimation: I. Theory and numerical studies[J]. Soil Science Society of America Journal, 1985, 49: 1348-1354.

    [5] VAN DAM J C, STRICKER J N M, DROOGERS P. Evaluation of the inverse method for determining soil hydraulic functions from one-step outflow experiments[J]. Soil Science Society of America Journal, 1992, 56: 1042-1050.

    [6] TOORMAN A F, WIERENGA P J, HILLS R G. Parameter estimation of hydraulic properties from one- step outflow data[J]. Water Resources Research, 1992, 28(11): 3021-3028.

    [7] ECHING S O, HOPMANS J W. Optimization of hydraulic functions from transient outflow and soil-water pressure data[J]. Soil Science Society of America Journal, 1993, 57(5): 1167-1175.

    [8] ECHING S O, HOPMANS J W, WENDROTH O. Unsaturated hydraulic conductivity from transient multistep outflow and soil-water pressure data[J]. Soil Science Society of America Journal, 1994, 58: 687-695.

    [9] VAN DAM J C, STRICKER J N M, DROOGERS P. Inverse method to determine soil hydraulic functions from multistep outflow experiments[J]. Soil Science Society of America Journal, 1994, 58: 647-652.

    [10] WEI C F, DEWOLKAR M M. Formulation of capillary hysteresis with internal state variables[J]. Water Resources Research, 2006, 42(7): 1-16.

    [11] WEI C F, MURALEETHARAN K K. Linear viscoelastic behavior of porous media with non-uniform saturation[J]. International Journal of Engineering Science, 2007, 45(2-8): 698-715.

    [12] VAN GENUCHTEN M T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal, 1980, 44: 892-898.

    [13] MUALEM Y. A new model for predicting the hydraulic conductivity of unsaturated porous media[J]. Water Resources Research, 1976, 12(3): 513-522.

    猜你喜歡
    水力學(xué)非飽和吸力
    深水大型吸力錨測(cè)試技術(shù)
    ROV在海上吸力樁安裝場(chǎng)景的應(yīng)用及安裝精度和風(fēng)險(xiǎn)控制
    化工管理(2022年11期)2022-06-03 07:08:24
    飽和紫色土初始態(tài)和穩(wěn)定態(tài)細(xì)溝水力學(xué)特征研究*
    深水吸力樁施工技術(shù)研究
    非飽和原狀黃土結(jié)構(gòu)強(qiáng)度的試驗(yàn)研究
    二維水力學(xué)模型在紅光大橋洪水影響評(píng)價(jià)中的應(yīng)用
    非飽和多孔介質(zhì)應(yīng)力滲流耦合分析研究
    基于管網(wǎng)理論的人口遷移動(dòng)力學(xué)模型構(gòu)建
    科技視界(2016年27期)2017-03-14 23:09:34
    非飽和土基坑剛性擋墻抗傾覆設(shè)計(jì)與參數(shù)分析
    自排式沉沙池上游渠道水力學(xué)特性研究
    亚洲精品国产一区二区精华液| 亚洲全国av大片| 精品国内亚洲2022精品成人 | 国产在线免费精品| 久久久精品区二区三区| 9色porny在线观看| 久久精品成人免费网站| 亚洲精品美女久久久久99蜜臀| 999精品在线视频| 窝窝影院91人妻| 最新在线观看一区二区三区| 伊人亚洲综合成人网| 80岁老熟妇乱子伦牲交| 亚洲熟女毛片儿| 亚洲av片天天在线观看| svipshipincom国产片| av网站免费在线观看视频| 亚洲,欧美精品.| 少妇粗大呻吟视频| 日本av手机在线免费观看| av福利片在线| 免费高清在线观看日韩| 国产亚洲av片在线观看秒播厂| 国产日韩欧美视频二区| 热re99久久国产66热| 两人在一起打扑克的视频| 国产成人精品在线电影| 成人国产一区最新在线观看| 一区二区三区四区激情视频| 99久久国产精品久久久| 午夜福利免费观看在线| 亚洲av成人不卡在线观看播放网 | bbb黄色大片| 精品免费久久久久久久清纯 | 美女高潮到喷水免费观看| 黄片大片在线免费观看| 欧美久久黑人一区二区| 啦啦啦啦在线视频资源| 久久久国产精品麻豆| 99热国产这里只有精品6| 黄色怎么调成土黄色| 动漫黄色视频在线观看| 久久久国产精品麻豆| 国产欧美日韩综合在线一区二区| 国产成人av激情在线播放| 国产视频一区二区在线看| 久热这里只有精品99| 亚洲黑人精品在线| 成人影院久久| 国产亚洲一区二区精品| 一本一本久久a久久精品综合妖精| av超薄肉色丝袜交足视频| 久久久国产一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久毛片免费看一区二区三区| 满18在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 一级,二级,三级黄色视频| 亚洲国产中文字幕在线视频| 国产精品国产三级国产专区5o| 久久毛片免费看一区二区三区| 麻豆av在线久日| 国产亚洲午夜精品一区二区久久| 国产精品一二三区在线看| 午夜福利,免费看| 欧美精品一区二区免费开放| 美女福利国产在线| videos熟女内射| 国产一区二区三区在线臀色熟女 | 夜夜骑夜夜射夜夜干| a级毛片黄视频| 国产真人三级小视频在线观看| 高清黄色对白视频在线免费看| 一级片'在线观看视频| 曰老女人黄片| 国产精品二区激情视频| 热99re8久久精品国产| 91精品国产国语对白视频| 三上悠亚av全集在线观看| 19禁男女啪啪无遮挡网站| 在线观看免费午夜福利视频| 一区二区三区乱码不卡18| 少妇精品久久久久久久| 狠狠婷婷综合久久久久久88av| 国产熟女午夜一区二区三区| 亚洲av片天天在线观看| avwww免费| 一级片免费观看大全| 亚洲欧洲精品一区二区精品久久久| 亚洲专区国产一区二区| 国产在线免费精品| 日韩大片免费观看网站| 日韩 亚洲 欧美在线| 国产一卡二卡三卡精品| 精品国产一区二区久久| 一区二区三区乱码不卡18| 国产日韩欧美在线精品| 色精品久久人妻99蜜桃| 日日爽夜夜爽网站| 波多野结衣av一区二区av| 欧美激情高清一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 女性被躁到高潮视频| 久久精品熟女亚洲av麻豆精品| 交换朋友夫妻互换小说| 淫妇啪啪啪对白视频 | 黄片大片在线免费观看| 9191精品国产免费久久| 免费黄频网站在线观看国产| 伊人亚洲综合成人网| 中国国产av一级| 久久久久久久大尺度免费视频| 五月开心婷婷网| 99久久精品国产亚洲精品| 女人爽到高潮嗷嗷叫在线视频| 在线亚洲精品国产二区图片欧美| kizo精华| 日本撒尿小便嘘嘘汇集6| 无限看片的www在线观看| 国产淫语在线视频| 波多野结衣一区麻豆| 91精品国产国语对白视频| 精品亚洲成a人片在线观看| 男女之事视频高清在线观看| 免费久久久久久久精品成人欧美视频| 成人av一区二区三区在线看 | 亚洲国产欧美网| 午夜福利视频在线观看免费| 欧美亚洲 丝袜 人妻 在线| 9色porny在线观看| 嫩草影视91久久| 免费观看人在逋| 多毛熟女@视频| 啦啦啦啦在线视频资源| 免费观看av网站的网址| 青青草视频在线视频观看| 女人被躁到高潮嗷嗷叫费观| 交换朋友夫妻互换小说| 午夜福利视频精品| 日韩欧美一区二区三区在线观看 | av有码第一页| 最新在线观看一区二区三区| 97人妻天天添夜夜摸| 另类亚洲欧美激情| 视频在线观看一区二区三区| 91精品国产国语对白视频| 嫩草影视91久久| 亚洲国产精品一区二区三区在线| 国产人伦9x9x在线观看| 亚洲少妇的诱惑av| 丝袜人妻中文字幕| 99国产精品99久久久久| 亚洲第一av免费看| 一区二区三区精品91| 亚洲中文日韩欧美视频| 午夜福利一区二区在线看| 欧美日韩视频精品一区| 99国产精品一区二区三区| 久久天堂一区二区三区四区| 99热国产这里只有精品6| 美国免费a级毛片| 一区二区日韩欧美中文字幕| 国产成人av教育| 欧美激情久久久久久爽电影 | 亚洲天堂av无毛| 精品久久久精品久久久| 精品乱码久久久久久99久播| 最新在线观看一区二区三区| 免费在线观看影片大全网站| 人成视频在线观看免费观看| 久久国产亚洲av麻豆专区| 天天影视国产精品| 在线观看免费日韩欧美大片| 一本综合久久免费| 欧美黑人欧美精品刺激| 久久久久国产一级毛片高清牌| 黄网站色视频无遮挡免费观看| 国产免费现黄频在线看| 国产高清videossex| 嫩草影视91久久| 精品福利永久在线观看| 9191精品国产免费久久| 中文字幕人妻丝袜一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 丰满饥渴人妻一区二区三| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看免费日韩欧美大片| 在线十欧美十亚洲十日本专区| 精品国产一区二区三区四区第35| 欧美精品一区二区免费开放| 女性被躁到高潮视频| 91精品三级在线观看| 精品一区二区三区四区五区乱码| 水蜜桃什么品种好| 午夜福利在线免费观看网站| 一级毛片电影观看| 日本91视频免费播放| 麻豆av在线久日| 国产1区2区3区精品| 女人高潮潮喷娇喘18禁视频| 国产黄色免费在线视频| 99久久99久久久精品蜜桃| 美女中出高潮动态图| 如日韩欧美国产精品一区二区三区| 精品国产国语对白av| 日韩欧美一区视频在线观看| 精品久久久久久电影网| 动漫黄色视频在线观看| 美女大奶头黄色视频| 久久99热这里只频精品6学生| 精品亚洲乱码少妇综合久久| 中国美女看黄片| 久久 成人 亚洲| 久热这里只有精品99| 热99re8久久精品国产| 无限看片的www在线观看| 欧美97在线视频| 成人手机av| 日韩免费高清中文字幕av| 91麻豆精品激情在线观看国产 | av在线老鸭窝| 亚洲人成电影免费在线| 久久人人爽人人片av| 叶爱在线成人免费视频播放| 啦啦啦视频在线资源免费观看| 久久亚洲精品不卡| 桃花免费在线播放| 美女视频免费永久观看网站| 一级,二级,三级黄色视频| 国精品久久久久久国模美| 精品国产一区二区三区四区第35| 两人在一起打扑克的视频| 日韩中文字幕欧美一区二区| 亚洲一码二码三码区别大吗| 精品卡一卡二卡四卡免费| 18禁国产床啪视频网站| 亚洲国产欧美在线一区| 女人高潮潮喷娇喘18禁视频| 亚洲精品久久久久久婷婷小说| 一本久久精品| 国产三级黄色录像| 欧美大码av| 国产野战对白在线观看| 精品国产乱码久久久久久小说| 一二三四在线观看免费中文在| 大陆偷拍与自拍| 成人18禁高潮啪啪吃奶动态图| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产成人免费| 下体分泌物呈黄色| 亚洲精品中文字幕一二三四区 | 亚洲全国av大片| 中国美女看黄片| 一本综合久久免费| av在线天堂中文字幕| 黄色a级毛片大全视频| 午夜a级毛片| 午夜老司机福利片| 麻豆久久精品国产亚洲av| 国产精品 国内视频| 在线国产一区二区在线| 国产av不卡久久| 久久久久精品国产欧美久久久| av中文乱码字幕在线| 国产精华一区二区三区| 在线a可以看的网站| 在线观看美女被高潮喷水网站 | 日韩欧美三级三区| 亚洲国产精品久久男人天堂| av国产免费在线观看| 国产成人av教育| 精品国产亚洲在线| a级毛片在线看网站| 亚洲av成人av| 天堂影院成人在线观看| 黄色女人牲交| 精品午夜福利视频在线观看一区| 亚洲国产精品sss在线观看| 亚洲人与动物交配视频| 我的老师免费观看完整版| 1024视频免费在线观看| 国产精品98久久久久久宅男小说| 亚洲黑人精品在线| 日韩精品中文字幕看吧| а√天堂www在线а√下载| 亚洲国产欧美网| 天堂动漫精品| 黄色a级毛片大全视频| 久久精品国产亚洲av香蕉五月| 亚洲成a人片在线一区二区| 免费在线观看影片大全网站| 嫁个100分男人电影在线观看| 少妇熟女aⅴ在线视频| 久久久国产精品麻豆| 怎么达到女性高潮| 少妇被粗大的猛进出69影院| 亚洲精品国产一区二区精华液| 亚洲乱码一区二区免费版| 国产精品美女特级片免费视频播放器 | 国产高清视频在线播放一区| 两性夫妻黄色片| 国产乱人伦免费视频| 丝袜美腿诱惑在线| 国产1区2区3区精品| 欧美黑人欧美精品刺激| 亚洲av成人精品一区久久| 久久精品夜夜夜夜夜久久蜜豆 | 日本免费一区二区三区高清不卡| 后天国语完整版免费观看| 亚洲精品一卡2卡三卡4卡5卡| 日本a在线网址| 在线国产一区二区在线| 亚洲狠狠婷婷综合久久图片| 9191精品国产免费久久| 国产精品一区二区三区四区久久| 久久人妻av系列| 99久久精品热视频| 高清在线国产一区| 国产精品av视频在线免费观看| 欧美一级毛片孕妇| 村上凉子中文字幕在线| 国产成人aa在线观看| 欧美日韩黄片免| 一本久久中文字幕| 精品国产亚洲在线| 午夜精品在线福利| 人妻夜夜爽99麻豆av| 午夜福利视频1000在线观看| 全区人妻精品视频| 高清毛片免费观看视频网站| 精品国产亚洲在线| 好男人在线观看高清免费视频| 亚洲欧美精品综合一区二区三区| 日韩av在线大香蕉| 国产精品98久久久久久宅男小说| 99久久无色码亚洲精品果冻| 国产精品一区二区精品视频观看| 亚洲精华国产精华精| 日本熟妇午夜| 51午夜福利影视在线观看| 亚洲成av人片在线播放无| 亚洲av日韩精品久久久久久密| 亚洲欧美日韩高清在线视频| 天堂√8在线中文| 50天的宝宝边吃奶边哭怎么回事| АⅤ资源中文在线天堂| 久久 成人 亚洲| 黄色视频不卡| 成人欧美大片| or卡值多少钱| 亚洲国产日韩欧美精品在线观看 | 99精品欧美一区二区三区四区| 国产激情久久老熟女| 精华霜和精华液先用哪个| 国产97色在线日韩免费| 国产精品久久久av美女十八| 成年女人毛片免费观看观看9| 亚洲精品一区av在线观看| 亚洲人成伊人成综合网2020| 亚洲美女黄片视频| 久久精品国产99精品国产亚洲性色| 最好的美女福利视频网| 首页视频小说图片口味搜索| 少妇人妻一区二区三区视频| 国产v大片淫在线免费观看| 亚洲乱码一区二区免费版| 两人在一起打扑克的视频| 女人高潮潮喷娇喘18禁视频| 狂野欧美激情性xxxx| 亚洲国产精品sss在线观看| 欧美3d第一页| 男女下面进入的视频免费午夜| 亚洲avbb在线观看| 亚洲精品久久成人aⅴ小说| 国产成人啪精品午夜网站| 国产主播在线观看一区二区| 久久伊人香网站| 一a级毛片在线观看| 欧美日韩亚洲国产一区二区在线观看| 99久久精品热视频| 精品久久久久久久久久免费视频| 日韩三级视频一区二区三区| 看黄色毛片网站| 舔av片在线| netflix在线观看网站| 舔av片在线| 国产精品亚洲美女久久久| 黄频高清免费视频| 免费无遮挡裸体视频| 亚洲aⅴ乱码一区二区在线播放 | 香蕉国产在线看| 国产一级毛片七仙女欲春2| 好男人在线观看高清免费视频| 99久久久亚洲精品蜜臀av| 精品熟女少妇八av免费久了| 97碰自拍视频| 亚洲精华国产精华精| 黑人操中国人逼视频| 美女黄网站色视频| 国产69精品久久久久777片 | 波多野结衣高清作品| 亚洲欧美日韩高清在线视频| 九九热线精品视视频播放| 亚洲成av人片在线播放无| 免费看美女性在线毛片视频| 国产精品一区二区三区四区免费观看 | 少妇被粗大的猛进出69影院| 淫妇啪啪啪对白视频| 久久久水蜜桃国产精品网| 日日夜夜操网爽| a级毛片a级免费在线| 成人永久免费在线观看视频| 夜夜爽天天搞| 精品国产超薄肉色丝袜足j| 久9热在线精品视频| 18禁观看日本| www国产在线视频色| 午夜a级毛片| 天堂影院成人在线观看| 日韩三级视频一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 男女之事视频高清在线观看| 精品福利观看| 91国产中文字幕| 国产高清视频在线播放一区| 亚洲免费av在线视频| 一级a爱片免费观看的视频| 99热这里只有精品一区 | 亚洲中文av在线| 亚洲国产日韩欧美精品在线观看 | 亚洲无线在线观看| 日本熟妇午夜| 高清毛片免费观看视频网站| 精品电影一区二区在线| 可以免费在线观看a视频的电影网站| 欧美乱色亚洲激情| 19禁男女啪啪无遮挡网站| 中文字幕高清在线视频| 国产精品精品国产色婷婷| 国产亚洲精品av在线| 久久久久久国产a免费观看| 国产1区2区3区精品| 免费看美女性在线毛片视频| 亚洲精品一卡2卡三卡4卡5卡| 久久性视频一级片| 久久久久亚洲av毛片大全| 成人国语在线视频| 国产又色又爽无遮挡免费看| 国产精品一区二区免费欧美| www.999成人在线观看| 亚洲av成人一区二区三| 免费人成视频x8x8入口观看| 三级男女做爰猛烈吃奶摸视频| tocl精华| 欧美日本亚洲视频在线播放| 国产精品永久免费网站| 亚洲精品粉嫩美女一区| 亚洲一区二区三区不卡视频| 久久天躁狠狠躁夜夜2o2o| 日韩中文字幕欧美一区二区| 中文字幕精品亚洲无线码一区| 午夜影院日韩av| 国产成人啪精品午夜网站| 亚洲精品色激情综合| 日日摸夜夜添夜夜添小说| 免费观看人在逋| 久热爱精品视频在线9| 欧美色视频一区免费| 女人高潮潮喷娇喘18禁视频| 久久伊人香网站| 深夜精品福利| 女人爽到高潮嗷嗷叫在线视频| 男女视频在线观看网站免费 | 97超级碰碰碰精品色视频在线观看| 国产在线观看jvid| 两个人看的免费小视频| 欧美性猛交黑人性爽| 在线十欧美十亚洲十日本专区| 久久天躁狠狠躁夜夜2o2o| 成人av一区二区三区在线看| 可以免费在线观看a视频的电影网站| 国产精品久久久久久精品电影| 黄色a级毛片大全视频| 久久久精品欧美日韩精品| www日本在线高清视频| 嫁个100分男人电影在线观看| 好男人电影高清在线观看| 特大巨黑吊av在线直播| 男女床上黄色一级片免费看| 黑人操中国人逼视频| 欧美黑人巨大hd| 少妇粗大呻吟视频| 亚洲欧美精品综合久久99| 日韩欧美三级三区| 搡老熟女国产l中国老女人| 国产精品久久视频播放| 成人国产一区最新在线观看| 国产精品久久电影中文字幕| 国产伦在线观看视频一区| 日韩av在线大香蕉| 午夜视频精品福利| 成熟少妇高潮喷水视频| 老汉色av国产亚洲站长工具| 久久精品亚洲精品国产色婷小说| 嫁个100分男人电影在线观看| 久久欧美精品欧美久久欧美| www.精华液| 国产高清视频在线观看网站| 搞女人的毛片| 欧美乱妇无乱码| 男女午夜视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 老汉色av国产亚洲站长工具| 国产成+人综合+亚洲专区| 免费搜索国产男女视频| 日本黄大片高清| 亚洲成人久久爱视频| 丁香欧美五月| 国内毛片毛片毛片毛片毛片| 国产黄片美女视频| 男插女下体视频免费在线播放| 日韩精品青青久久久久久| 国产成人啪精品午夜网站| 久久香蕉国产精品| 精品福利观看| 1024手机看黄色片| 国产真人三级小视频在线观看| 高潮久久久久久久久久久不卡| 又粗又爽又猛毛片免费看| 蜜桃久久精品国产亚洲av| 国产午夜精品久久久久久| 又黄又爽又免费观看的视频| 欧美日韩一级在线毛片| 一区福利在线观看| 一级毛片高清免费大全| 精品久久久久久久末码| 99国产极品粉嫩在线观看| 久久伊人香网站| 欧美三级亚洲精品| 成人国语在线视频| 天堂动漫精品| 国内精品一区二区在线观看| 亚洲人成电影免费在线| 久久精品aⅴ一区二区三区四区| 天堂av国产一区二区熟女人妻 | 亚洲国产日韩欧美精品在线观看 | 免费人成视频x8x8入口观看| 亚洲全国av大片| 国产精品永久免费网站| 天天添夜夜摸| 亚洲电影在线观看av| 久久久国产成人精品二区| 五月玫瑰六月丁香| 亚洲国产精品合色在线| 国产成人aa在线观看| 国产蜜桃级精品一区二区三区| 欧美高清成人免费视频www| 日韩大码丰满熟妇| 午夜福利18| 亚洲欧美日韩高清在线视频| 久久久久久久精品吃奶| 嫩草影视91久久| 日韩欧美国产在线观看| 少妇裸体淫交视频免费看高清 | 在线观看午夜福利视频| 亚洲精品在线观看二区| 国产精品一区二区免费欧美| 久久精品国产亚洲av高清一级| 老熟妇仑乱视频hdxx| 久久伊人香网站| 亚洲性夜色夜夜综合| 欧美三级亚洲精品| 人妻久久中文字幕网| 亚洲国产欧美人成| 精品不卡国产一区二区三区| 国产精品亚洲av一区麻豆| 天堂√8在线中文| 777久久人妻少妇嫩草av网站| 亚洲人成伊人成综合网2020| 99热只有精品国产| 99久久99久久久精品蜜桃| 精品国产乱码久久久久久男人| 国产精品乱码一区二三区的特点| 最近最新免费中文字幕在线| 免费看十八禁软件| 真人做人爱边吃奶动态| 国产一区二区三区在线臀色熟女| 九色国产91popny在线| 亚洲中文av在线| 久久人妻av系列| 免费高清视频大片| 国产午夜精品论理片| 夜夜爽天天搞| 波多野结衣高清无吗| 熟女少妇亚洲综合色aaa.| av欧美777| 亚洲第一电影网av| 亚洲国产欧洲综合997久久,| 精品国内亚洲2022精品成人| 亚洲狠狠婷婷综合久久图片| 国产精品98久久久久久宅男小说| 精品熟女少妇八av免费久了| 国产一区在线观看成人免费| 黄色a级毛片大全视频| 久久久久久久午夜电影| 好看av亚洲va欧美ⅴa在| 一级片免费观看大全| 免费看a级黄色片| 黄色a级毛片大全视频| 日韩三级视频一区二区三区| 在线国产一区二区在线| 国产av又大| 色综合亚洲欧美另类图片| 欧美性长视频在线观看| 日本一二三区视频观看|