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

    隨機裂隙巖體滲流-應(yīng)力耦合分析的塊體單元法

    2011-09-20 06:18:18殷德勝汪衛(wèi)明陳勝宏
    巖土力學(xué) 2011年9期
    關(guān)鍵詞:塊體滲流裂隙

    殷德勝,汪衛(wèi)明,陳勝宏

    (武漢大學(xué) 水資源與水電工程科學(xué)國家重點實驗室,武漢 430072)

    1 引 言

    1959年法國Malpasset拱壩的潰壩事故,導(dǎo)致裂隙巖體滲流與應(yīng)力的耦合問題成為巖石水力學(xué)的熱門課題。由于該問題的邊界條件和材料性質(zhì)的復(fù)雜性,很難進行物理模型試驗,人們普遍采用數(shù)值解法。

    有限單元法是連續(xù)介質(zhì)力學(xué)的代表,它有著完備的理論體系,在巖體滲流-應(yīng)力耦合研究中的應(yīng)用最為廣泛。Noorishad等[1]利用土體固結(jié)分析理論,并運用有限元法中的節(jié)理單元,研究了非連續(xù)介質(zhì)中的固液兩相介質(zhì)的耦合間題;Oda[2]將裂隙巖體滲透張量與彈性張量統(tǒng)一用裂隙網(wǎng)絡(luò)幾何張量表示,建立了較為嚴(yán)密的裂隙巖體滲流特性與應(yīng)力特性之間的耦合關(guān)系;在國內(nèi),陳平和張有天[3]考慮了裂隙變形的非線性特性,以裂隙平均隙寬為計算參數(shù),對立方定理作了修正,提出了滲流-應(yīng)力耦合計算方法;王媛等[4]則建立的同時以滲流水壓力和位移為未知值的4自由度全耦合法,這降低了耦合分析的難度。由于已有多種成熟的商業(yè)軟件如ABAQUS、ADINA等能夠模擬多場耦合問題,這將進一步加強有限元法在巖土工程滲流-應(yīng)力耦合這一領(lǐng)域中的應(yīng)用。

    對于復(fù)雜的工程,已有的研究成果大多基于二維情況,由于三維裂隙巖體的幾何復(fù)雜性,目前尚無一個強大的前處理器能夠輕松地對其進行網(wǎng)格離散,這是限制有限元法在該領(lǐng)域進一步應(yīng)用的主要原因。學(xué)者們只能退而求其次,大多采用等效連續(xù)介質(zhì)模型進行耦合研究。然而,裂隙中的水主要通過裂隙等結(jié)構(gòu)面流動,其滲流與傳統(tǒng)的多孔介質(zhì)滲流存在著本質(zhì)上的差別。對裂隙巖體進行滲流分析,宜采用裂隙網(wǎng)絡(luò)模型或雙重介質(zhì)模型[5]。

    雙重介質(zhì)模型除裂隙網(wǎng)格外,將巖塊視為滲透系數(shù)較小的連續(xù)滲透介質(zhì),且能考慮巖塊孔隙與巖體裂隙之間水量的交換,是更接近于實際的模型。柴軍瑞等[6]將各種裂隙和孔隙按規(guī)模和滲透性分為4級,建立起巖體滲流場與應(yīng)力場耦合分析的多重裂隙網(wǎng)絡(luò)模型。綜合來說,同連續(xù)介質(zhì)方法相比,雙重介質(zhì)模型優(yōu)勢明顯,但由于繼承了網(wǎng)格剖分的困難,實現(xiàn)難度相對更大。綜合起來,不連續(xù)介質(zhì)力學(xué)分析方法結(jié)合裂隙網(wǎng)絡(luò)滲流模型則是個相對合理可行的思路。

    不連續(xù)變形分析方法(discontinuous deformation analysis method,DDA)屬于不連續(xù)介質(zhì)力學(xué)范疇,十分適合處理裂隙巖體滲流、應(yīng)力的計算問題。張國新等[7]應(yīng)用DDA法,對日本某事故隧道進行了二維滲流-應(yīng)力耦合分析,得到了邊坡失穩(wěn)坍塌的機制。應(yīng)用DDA法進行三維滲流-應(yīng)力的耦合分析,在計算中尚存在接觸難以準(zhǔn)確確定、罰函數(shù)如何施加等諸多困難,目前的研究成果很少。

    塊體單元法也屬于不連續(xù)介質(zhì)力學(xué)范疇,是一種能對節(jié)理裂隙巖體進行滲流、應(yīng)力分析的數(shù)值方法[8]。在塊體單元法的基礎(chǔ)上,本文考慮了節(jié)理裂隙的剛度和導(dǎo)水系數(shù)隨開度的變化關(guān)系,并將常規(guī)的塊體單元法結(jié)構(gòu)計算模塊和滲流分析模塊整合成一個系統(tǒng),應(yīng)用兩場交叉迭代技術(shù),完成了隨機裂隙巖體的滲流-應(yīng)力耦合研究。該方法前、后處理簡單,計算過程也不復(fù)雜,能有效地應(yīng)用于三維隨機裂隙巖體的耦合分析中。

    2 塊體單元法原理簡介

    巖體的變形主要都集中于節(jié)理、裂隙等結(jié)構(gòu)面上,相比而言,塊體的變形可以忽略不計。在此假設(shè)條件下,塊體單元法認為塊體為剛體,只考慮結(jié)構(gòu)面的變形和強度特性,以塊體單元形心位移為基本未知量,通過建立塊體單元的平衡方程、結(jié)構(gòu)面變形和塊體單元位移的幾何相容方程、結(jié)構(gòu)面的本構(gòu)方程,綜合推導(dǎo)出塊體系統(tǒng)的整體平衡方程。

    同時,斷層、夾層、節(jié)理和裂隙等結(jié)構(gòu)面不僅是巖體中的軟弱面,且由于其滲透系數(shù)通常遠大于完整巖石,故它們也是主要的滲流通道。假定巖塊不透水,即在巖塊中滲透系數(shù)k =0,滲流只發(fā)生在裂隙面上,根據(jù)裂隙網(wǎng)絡(luò)交叉處的流量守恒和分片變分原理,建立了與塊體單元法相配套的隨機裂隙巖體的滲流分析方法[9]。

    有關(guān)塊體單元法的應(yīng)力、滲流分析的詳細過程請參閱文獻[8]。

    3 耦合機制

    巖體的變形大部分由結(jié)構(gòu)面產(chǎn)生,結(jié)構(gòu)面隙寬的微小改變將引起滲流量的重大變化。同時,滲流場的變化會引起水荷載的變化,從而導(dǎo)致應(yīng)力場的重分布。因而,裂隙巖體的結(jié)構(gòu)場和滲流場之間存在著重要的耦合關(guān)系。

    然而,在應(yīng)力場和滲流場相互影響的定量描述上,目前尚無定論。由于裂隙面上的法向應(yīng)力引起裂隙變形相對顯著,為抓住主要矛盾,本文僅研究法向應(yīng)力對裂隙變形特性、滲透特性的影響。

    3.1 結(jié)構(gòu)面變形對裂隙剛度系數(shù)的影響

    結(jié)構(gòu)面的相對變形影響著剛度系數(shù),從而間接地影響應(yīng)力場的分布。陳勝宏等[10]曾提出有關(guān)節(jié)理面變形和應(yīng)力關(guān)系的“充填模型”,將節(jié)理裂隙視為一種均勻介質(zhì),并據(jù)此推導(dǎo)了有關(guān)節(jié)理面法向、切向剛度系數(shù)與應(yīng)力之間簡單而實用的關(guān)系式。

    以單裂隙巖體說明裂隙滲流與應(yīng)力的耦合關(guān)系。設(shè)圖1中的裂隙面受法向應(yīng)力σz和切向應(yīng)力τzx,τzy作用,由于裂隙面的開度較小,故應(yīng)變滿足:

    式中: εx、 εy分別為x、y向的正應(yīng)變, γxy、 γyx為xy平面上的剪應(yīng)變。

    圖1 單裂隙巖體Fig.1 Rock mass containing one fracture

    “充填模型”假設(shè)裂隙面上,“充填介質(zhì)”是各向同性體。據(jù)此,可得到裂隙面局部坐標(biāo)下的彈性本構(gòu)關(guān)系為

    式中:λ、G為裂隙內(nèi)充填物的拉梅常數(shù)。記e0為裂隙面應(yīng)力為0時的開度,ux、 uy、 uz分別為x、y、z三向的相對位移,則裂隙開度滿足:

    于是式(2)可以變成:

    式中: kn、 ks分別為裂隙面的法向、切向的剛度系數(shù),可由充填物的彈性參數(shù)和開度確定:

    由于開度e與法向變形uz相關(guān),而uz又與應(yīng)力有關(guān)聯(lián),故 kn、 ks都是應(yīng)力的函數(shù)。

    從式(4)中取第1行,可以得到:

    結(jié)合公式(3),上式可寫成:

    將式(7)積分后得到:

    把式(9)代入式(5),即可得到法向、切向剛度系數(shù)與法向應(yīng)力的關(guān)系:

    記 kn0= λ+ 2G /e0,ks0= G/e0,用它來表示初始剛度,即裂隙在無應(yīng)力狀態(tài)時的剛度系數(shù),則式(10)可寫成:

    式(11)表明了裂隙剛度系數(shù)與法向應(yīng)力之間存在著指數(shù)關(guān)系。

    3.2 結(jié)構(gòu)面變形對導(dǎo)水系數(shù)的影響

    立方定律是研究裂隙巖體滲流的基礎(chǔ)。然而,天然節(jié)理并非完全光滑平直,也不是純粹的孔隙,因而立方定律的適用條件有限。近幾十年來,許多考慮了節(jié)理面的粗糙性和凸起物高度分布情況的修正立方定律[11-12]相繼被提出,其主要思想就是用等效水力隙寬代替幾何平均隙寬。然而,等效水力隙寬和幾何平均隙寬是完全不同的兩個概念[13],若將其直接用于耦合研究,不但理論說不通,且流量與隙寬的3次方成正比,會導(dǎo)致耦合分析中計算過于敏感和收斂困難。故立方定律在耦合研究中直接應(yīng)用較少。

    “充填模型” 同樣適用于結(jié)構(gòu)面的滲流特性的研究。該模型中假設(shè)裂隙中填充物的滲透系數(shù)不變,而導(dǎo)水系數(shù)隨著開度的變化而變化,并結(jié)合式(9),得到:

    式中:kd為導(dǎo)水系數(shù);kd0為初始導(dǎo)水系數(shù),即裂隙面無應(yīng)力時的導(dǎo)水系數(shù),可通過試驗獲得,其他同式(10)。

    式(11)、(12)能較好地模擬裂隙剛度系數(shù)和導(dǎo)水系數(shù)與應(yīng)力之間的關(guān)系。本文運用這套關(guān)系式,結(jié)合塊體單元法進行了裂隙巖體滲流-應(yīng)力的耦合研究。

    3.3 耦合步驟

    根據(jù)結(jié)構(gòu)面變形對應(yīng)力場、滲流場的影響關(guān)系式,使用Fortran語言,編制了塊體元滲流-應(yīng)力的耦合分析程序。該程序中,充分考慮塊體元應(yīng)力、滲流分析的共同點,將2個分析模塊整合在一個程序中,顯著減少了公共計算部分的計算量,且相比傳統(tǒng)的兩場迭代,不需要2個程序相互調(diào)用時讀寫0時文件的時間開銷,且降低了文件傳遞導(dǎo)致的累積誤差,耦合計算的效率和精度均較高。

    程序中,耦合分析的具體實施步驟如下:

    (1)由初始剛度系數(shù)和初始導(dǎo)水系數(shù),結(jié)合邊界條件,計算得到初始應(yīng)力場 {σ}(0)和初始滲流場{φ}(0);(2)結(jié)合式(11)、(12),根據(jù){σ}(0)計算裂隙面新的法向、切向剛度系數(shù)和導(dǎo)水系數(shù),并重新生成結(jié)構(gòu)計算的剛度矩陣和滲流計算的導(dǎo)水矩陣;(3)由新的導(dǎo)水矩陣和滲流邊界條件,重新計算滲流場{φ}(i);(4)由新的滲流場 {φ}(i)計算出滲透荷載,并加上其它荷載,重新計算應(yīng)力場 {σ}(i);(5)計算應(yīng)力、水頭的增量: {Δ σ} ={σ}(i)-{σ}(i-1)和{Δφ}={φ}(i)-{φ}(i-1),并取{Δσ}和{Δφ}的2范數(shù)。若和小于給定的誤差值,則迭代結(jié)束,否則重新執(zhí)行步驟(2),直至滲流場和應(yīng)力場都滿足收斂精度為止。

    4 算例驗證

    薛孌鸞[14]曾用基于雙重介質(zhì)模型的復(fù)合單元法研究過裂隙巖體的滲流-應(yīng)力耦合問題。本文采用她文中的算例,從規(guī)律上對本文方法進行驗證。如圖2所示的長方體巖塊試件,坐標(biāo)系在左下角,試件的x方向長為1.0 m,y方向?qū)挒?.1 m,z方向高為0.2 m。在z =0.1 m處有一條初始開度為1 mm的水平貫穿性裂隙。由于塊體單元法假設(shè)塊體為剛體,本文采用卓家壽等提出的界面元法[15-16],用自定義的虛擬結(jié)構(gòu)面將試件等分成若干塊體,將變形累積在結(jié)構(gòu)面上,以此來近似模擬。

    對于均質(zhì)材料,根據(jù)文[16]的建議,虛擬結(jié)構(gòu)面的法向剛度、切向剛度按式(13)選?。?/p>

    式中: h1、 h2為塊體形心到結(jié)構(gòu)面的距離;E為彈性模量;μ為泊松比。如圖2所示,總共引入9條虛擬結(jié)構(gòu)面,將試件分成20個塊體。另外,虛擬結(jié)構(gòu)面的導(dǎo)水系數(shù)的取值為實際裂隙導(dǎo)水系數(shù)的1%。所有結(jié)構(gòu)面的參數(shù)見表1:

    取巖塊試件的左側(cè)水位為 10 m,右側(cè)水位為5 m,其它外表面為不透水邊界??紤]水壓力并在上表面正中間的 23、53兩點處各施加集中力 F =0.1 MN,不考慮試件的自重。在除上表面外的其余5個側(cè)面均施加法向約束。

    圖2 算例塊體系統(tǒng)Fig.2 Block element system of the example

    表1 裂隙面的參數(shù)Table 1 Parameters of fractures

    整理耦合和不耦合兩種工況下的成果,圖3為結(jié)構(gòu)面上的水頭分布,圖4為不耦合和耦合兩種工況下y =0.05 m斷面上的流速矢量圖。

    圖3 結(jié)構(gòu)面上的水頭分布Fig.3 Hydraulic distribution on the structure

    圖4 流速矢量圖Fig.4 Hydraulic velocity vector diagram

    從圖3中可以看出,不考慮耦合作用時,結(jié)構(gòu)面上的水頭呈線性分布;考慮耦合作用后,集中力作用點的左側(cè)部位,結(jié)構(gòu)面的水頭趨于上游水位,右側(cè)部位的水頭則趨于下游水位。對比圖4可以看出:當(dāng)不考慮耦合時,整個裂隙面上的流速是相同的;考慮耦合后,集中力作用點底部壓應(yīng)力較大區(qū)域的導(dǎo)水系數(shù)減小,滲透坡降增加,流速增加,其他區(qū)域的流速減小。結(jié)果顯示了本文耦合機制和程序的合理性。

    5 工程算例

    圖5是貴州省烏江沙沱水電站非溢流壩段的典型剖面的示意圖。其巖基上含有F88斷層,j4、j5、j6三條夾層,并有兩組傾角約為45°的正交裂隙。沿壩軸線方向取寬度為1 m的壩段參與計算,對于斷層和夾層,采用確定性方式模擬;對于裂隙,采用蒙特卡羅法隨機生成,并用塊體識別方法[9]找出所有塊體。塊體系統(tǒng)如圖6所示。

    圖5 重力壩剖面Fig.5 Profile diagram of gravity dam

    圖6 塊體系統(tǒng)Fig.6 Block element system

    對于地基部分,斷層、夾層、裂隙等結(jié)構(gòu)面的參數(shù)可以通過實驗的方式獲得;對于壩體部分,仍按卓家壽的方法,按5 m的間距獲得虛擬結(jié)構(gòu)面的參數(shù)。帷幕采用文獻[17]的方法,即將帷幕塊體周邊結(jié)構(gòu)面視為相對不透水面,其導(dǎo)水系數(shù)取為較小值,本文取值為1×10-8m/s。

    斷層、夾層的材料參數(shù)見表 2,裂隙參數(shù)見表3,壩體和基巖的參數(shù)見表4。

    表2 結(jié)構(gòu)面參數(shù)Table 2 Parameters of discontinuities

    表3 裂隙參數(shù)Table 3 Parameters of fractures

    表4 基巖和壩體混凝土參數(shù)Table 4 Parameters of rock foundation and dam concrete

    同樣,計算也分為不耦合和耦合兩種工況。上下游水位見圖5,考慮的荷載有壩體和基巖的自重、靜水壓力和滲透荷載作用。

    整理了兩種工況的計算成果對比。圖7為裂隙網(wǎng)絡(luò)的水頭等勢線,圖8為水力坡降矢量圖。

    圖7 裂隙網(wǎng)絡(luò)等水頭線(單位:m)Fig.7 Hydraulic potential isolines of fracture network (unit: m)

    圖8 水力坡降矢量圖Fig.8 Hydraulic gradient vector diagrms

    從圖 7、8中可以看出:考慮耦合作用以后,壩前地基的水力坡降逐漸減小,水頭普遍增大;壩后基礎(chǔ)的變化規(guī)律與壩前類似;壩體下方地基的水平向主導(dǎo)的裂隙上的壓應(yīng)力增加,導(dǎo)水系數(shù)減小,坡降、流速均有所增加,進而導(dǎo)致建基面的揚壓力增加,對壩體的穩(wěn)定不利。

    6 結(jié) 論

    (1)塊體單元法本身即是為解決節(jié)理、裂隙巖體的應(yīng)力、滲流問題而提出,將其應(yīng)用于裂隙巖體的滲流-應(yīng)力耦合中,有著先天的優(yōu)勢。而且能考慮三維巖體中含有的復(fù)雜隨機裂隙網(wǎng)絡(luò),這是其他模型所不能考慮或很難考慮的;

    (2)法向應(yīng)力對裂隙的滲透特性影響很大。在高壓應(yīng)力的區(qū)域,導(dǎo)水系數(shù)減小,滲透坡降增加,流速增大;其他區(qū)域流速減小。從算例中可以看出,重力壩分析中若不考慮耦合作用,則可能導(dǎo)致?lián)P壓力的計算值偏小。對實際工程會產(chǎn)生不利影響。

    然而,就目前的計算機發(fā)展水平,采用塊體單元法進行滲流-應(yīng)力耦合分析的計算量仍顯偏大,需要進一步優(yōu)化數(shù)據(jù)結(jié)構(gòu)、改進積分算法,甚至采用分布式并行計算,使該方法真正意義上用于工程的三維分析中。

    [1]NOORISHAD J, AYATOLLAHI M S, WITHERSPOON P A. A finite-element method for coupled stress and fluid flow analysis in fractured rock masses[J]. International Journal of Rock Mechanics and Mining Sciences, 1982:185-193.

    [2]ODA M. An equivalent continuum model for coupled stress and fluid flow analysis in jointed rock masses[J].Water Resources Research, 1986, 22(13): 1845-1856.

    [3]陳平, 張有天. 裂隙巖體滲流與應(yīng)力耦合分析[J]. 巖石力學(xué)與工程學(xué)報, 1994, 13(4): 299-308.CHEN Ping, ZHANG You-tian. Coupling analysis of seepage/stress for joint rock[J]. Chinese Journal of Rock Mechanics and Engineering, 1994, 13(4): 299-308.

    [4]王媛, 徐志英, 速寶玉. 裂隙巖體滲流與應(yīng)力耦合分析的四自由度全耦合法[J]. 水利學(xué)報, 1998, (7): 55-59.WANG Yuan, XU Zhi-ying, SU Bao-yu. Four-freedom complete method for the seepage-stress coupled analysis in fissured rock masses[J]. Journal of Hydraulic Engineering, 1998, (7): 55-59.

    [5]張有天. 巖石水力學(xué)與工程[M]. 北京: 中國水利水電出版社, 2005.

    [6]柴軍瑞, 仵彥卿. 巖體滲流場與應(yīng)力場耦合分析的多重裂隙網(wǎng)絡(luò)模型[J]. 巖石力學(xué)與工程學(xué)報, 2000, 19(6):712-717.CHAI Jun-rui, WU Yan-qing. Research on multiple-level fracture network model for coupled seepage and stress fields in rock mass[J]. Chinese Journal of Rock Mechanics and Engineering, 2000, 19(6): 712-717.

    [7]張國新, 武曉峰. 裂隙滲流對巖石邊坡穩(wěn)定的影響——滲流、變形耦合作用的DDA法[J]. 巖石力學(xué)與工程學(xué)報, 2003, 22(8): 1269-1275.ZHANG Guo-xin, WU Xiao-feng. Influence of seepage on the stability of rock slope coupling of seepage and deformation by DDA method[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(8): 1269-1275.

    [8]陳勝宏. 計算巖體力學(xué)與工程[M]. 中國水利水電出版社, 2006.

    [9]殷德勝, 汪衛(wèi)明, 陳勝宏. 三維隨機裂隙巖體滲流分析的塊體單元法[J]. 巖土力學(xué), 2009, 30(8): 1269-1275.YIN De-sheng, WANG Wei-ming, CHEN Sheng-hong.Block element method for seepage analysis in three dimensional random fracture network[J]. Rock and Soil Mechanics, 2009, 30(8): 1269-1275.

    [10]陳勝宏, 王鴻儒, 熊文林. 節(jié)理面滲流性質(zhì)的探討[J].武漢水利電力學(xué)院學(xué)報, 1989, 22(1): 53-60.CHEN Sheng-hong, WANG Hong-ru, XIONG Wen-lin.Study on the seepage characteristics of joint surface[J].Journal of Wuhan University of Hydraulic and Electric Engineering, 1989, 22(1): 53-60.

    [11]LOMIZE G M. Flow in fractured rocks[M]. Moscow:Gesenergoizdat, 1951.

    [12]LOUIS C. Rock mechanics[M]. New York: Elsevier Science, 1974.

    [13]王媛. 單裂隙面滲流與應(yīng)力的耦合特性[J]. 巖石力學(xué)與工程學(xué)報, 2002, 21(1): 83-87.WANG Yuan. Coupling characteristic of stress and fluid flow within a single fracture[J]. Chinese Journal of Rock Mechanics and Engineering, 2002, 21(1): 83-87.

    [14]薛孌鸞, 陳勝宏. 巖石裂隙滲流與法向應(yīng)力耦合的復(fù)合單元模型[J]. 巖石力學(xué)與工程學(xué)報, 2007, 26(增刊):2613-2619.XUE Luan-luan, CHEN Sheng-hong. Composite element model of seepage-normal stress coupling for rock fractures[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(Supp.): 2613-2619.

    [15]卓家壽, 趙寧. 不連續(xù)介質(zhì)靜、動力分析的剛體-彈簧元法[J]. 河海大學(xué)學(xué)報, 1993, 21(5): 34-43.ZHUO Jia-shou, ZHAO Ning. Piecewise rigid body-interface spring method for problems of discontinuous medium[J]. Journal of Hohai University,1993, 21(5): 34-43.

    [16]方義琳, 卓家壽, 章青. 具有任意形狀單元離散模型的界面元法[J]. 工程力學(xué), 1998, 15(2): 27-37.FANG Yi-lin, ZHUO Jia-shou, ZHANG Qing. Interface stress element method for distinct model with elements of arbitrary shape[J]. Engineering Mechanics, 1998, 15(2):27-37.

    [17]汪衛(wèi)明, 徐明毅, 陳勝宏. 復(fù)雜邊界條件下的巖體網(wǎng)絡(luò)滲流分析[J]. 巖石力學(xué)與工程學(xué)報, 2001, 20(4): 473-476.WANG Wei-ming, XU Ming-yi, CHEN Sheng-hong.Analysis method of network seepage for blocky rock mass with complicated boundary condition[J]. Chinese Journal of Rock Mechanics and Engineering, 2001,20(4): 473-476.

    猜你喜歡
    塊體滲流裂隙
    裂隙腦室綜合征的診斷治療新進展
    一種新型單層人工塊體Crablock 的工程應(yīng)用
    水運工程(2020年11期)2020-11-27 06:36:58
    裂隙燈檢查的個性化應(yīng)用(下)
    《老炮兒》:在時代裂隙中揚棄焦慮
    電影新作(2016年1期)2016-02-27 09:16:48
    一種Zr 基塊體金屬玻璃的納米壓入蠕變行為研究
    上海金屬(2015年3期)2015-11-29 01:09:58
    塊體非晶合金及其應(yīng)用
    波浪作用下斜坡上護面塊體斷裂破壞的數(shù)值模擬
    水道港口(2015年1期)2015-02-06 01:25:35
    簡述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    非完全充填裂隙的應(yīng)力應(yīng)變特性研究
    河南科技(2014年11期)2014-02-27 14:09:42
    日本91视频免费播放| 日本一区二区免费在线视频| 婷婷成人精品国产| 丝瓜视频免费看黄片| 成人影院久久| 欧美久久黑人一区二区| 国精品久久久久久国模美| 免费日韩欧美在线观看| 日韩一区二区三区影片| 自线自在国产av| 欧美97在线视频| 久久天堂一区二区三区四区| 国产精品二区激情视频| 少妇裸体淫交视频免费看高清 | 建设人人有责人人尽责人人享有的| 又紧又爽又黄一区二区| 悠悠久久av| 新久久久久国产一级毛片| 久久久精品免费免费高清| 国产欧美亚洲国产| 永久免费av网站大全| 制服诱惑二区| 欧美一级毛片孕妇| 50天的宝宝边吃奶边哭怎么回事| 少妇猛男粗大的猛烈进出视频| 国产成人a∨麻豆精品| 久久久久网色| 丰满饥渴人妻一区二区三| 免费黄频网站在线观看国产| 亚洲视频免费观看视频| 国产有黄有色有爽视频| 1024视频免费在线观看| 在线观看免费午夜福利视频| 一本—道久久a久久精品蜜桃钙片| 国产男女内射视频| 啦啦啦在线免费观看视频4| 色婷婷久久久亚洲欧美| 亚洲国产成人一精品久久久| 精品国产一区二区三区四区第35| 亚洲av美国av| 亚洲国产精品一区二区三区在线| 中国美女看黄片| 99热网站在线观看| www.熟女人妻精品国产| 精品卡一卡二卡四卡免费| 一二三四在线观看免费中文在| 国产无遮挡羞羞视频在线观看| 亚洲精品一区蜜桃| 久久久欧美国产精品| 91成年电影在线观看| 久久天躁狠狠躁夜夜2o2o| 中文字幕人妻熟女乱码| 亚洲欧洲日产国产| 黑人操中国人逼视频| 精品少妇一区二区三区视频日本电影| 飞空精品影院首页| h视频一区二区三区| 久久久国产精品麻豆| 亚洲av电影在线观看一区二区三区| 精品免费久久久久久久清纯 | 亚洲专区国产一区二区| 蜜桃国产av成人99| 日韩精品免费视频一区二区三区| 亚洲成av片中文字幕在线观看| 国产av精品麻豆| 日韩大码丰满熟妇| 久久精品久久久久久噜噜老黄| 成年女人毛片免费观看观看9 | 美女高潮喷水抽搐中文字幕| 69av精品久久久久久 | 老司机午夜十八禁免费视频| 成人影院久久| 99国产精品一区二区三区| 午夜视频精品福利| 欧美黑人精品巨大| 亚洲人成77777在线视频| 成年人黄色毛片网站| 国产男人的电影天堂91| 日日摸夜夜添夜夜添小说| 亚洲精品第二区| 国产精品影院久久| 日本精品一区二区三区蜜桃| 欧美日韩亚洲综合一区二区三区_| 交换朋友夫妻互换小说| 18禁黄网站禁片午夜丰满| www.熟女人妻精品国产| 又黄又粗又硬又大视频| 久久久久久人人人人人| 在线观看免费视频网站a站| 亚洲人成77777在线视频| 国产av又大| 国产成人啪精品午夜网站| 国产亚洲精品一区二区www | 午夜成年电影在线免费观看| 丝袜美足系列| 啦啦啦啦在线视频资源| 国产主播在线观看一区二区| 黄色视频,在线免费观看| 亚洲色图 男人天堂 中文字幕| 最新在线观看一区二区三区| av国产精品久久久久影院| 水蜜桃什么品种好| 免费少妇av软件| 国产成人精品在线电影| 国产色视频综合| 国产成人啪精品午夜网站| 亚洲专区字幕在线| 国产一级毛片在线| 欧美激情 高清一区二区三区| 久久人妻福利社区极品人妻图片| 欧美亚洲 丝袜 人妻 在线| 国产男女内射视频| 免费在线观看视频国产中文字幕亚洲 | 欧美 日韩 精品 国产| 亚洲国产看品久久| 免费高清在线观看日韩| 在线观看免费日韩欧美大片| 亚洲精品国产一区二区精华液| 亚洲精品av麻豆狂野| 亚洲熟女毛片儿| 美女中出高潮动态图| 日韩欧美免费精品| 日韩有码中文字幕| 淫妇啪啪啪对白视频 | 亚洲精品久久成人aⅴ小说| 久久国产精品大桥未久av| 波多野结衣一区麻豆| 夫妻午夜视频| 9191精品国产免费久久| 久久久久久久大尺度免费视频| 欧美av亚洲av综合av国产av| www.av在线官网国产| 国产精品成人在线| 多毛熟女@视频| 亚洲自偷自拍图片 自拍| 黑人巨大精品欧美一区二区蜜桃| 亚洲av美国av| 午夜福利在线观看吧| 精品欧美一区二区三区在线| 一本久久精品| 韩国精品一区二区三区| 日本黄色日本黄色录像| 久久精品亚洲av国产电影网| 中文字幕色久视频| 久久精品国产综合久久久| 伊人久久大香线蕉亚洲五| 日本黄色日本黄色录像| 中文欧美无线码| 亚洲情色 制服丝袜| 国产日韩欧美视频二区| 成年动漫av网址| 亚洲七黄色美女视频| 性少妇av在线| 丝袜美足系列| www.自偷自拍.com| 欧美日本中文国产一区发布| 伊人亚洲综合成人网| 亚洲精品乱久久久久久| 9191精品国产免费久久| www.自偷自拍.com| 亚洲av电影在线进入| 黄色片一级片一级黄色片| 脱女人内裤的视频| 欧美亚洲日本最大视频资源| 91国产中文字幕| 老司机福利观看| 国产在视频线精品| 丝瓜视频免费看黄片| 成人免费观看视频高清| 如日韩欧美国产精品一区二区三区| 久久久久国产一级毛片高清牌| 无限看片的www在线观看| www.自偷自拍.com| 女性生殖器流出的白浆| 99热国产这里只有精品6| xxxhd国产人妻xxx| 国产成人欧美| 美女高潮喷水抽搐中文字幕| 在线观看免费日韩欧美大片| 久久午夜综合久久蜜桃| 秋霞在线观看毛片| 两个人免费观看高清视频| 99国产精品一区二区蜜桃av | 亚洲av日韩精品久久久久久密| 两性夫妻黄色片| 亚洲情色 制服丝袜| 老司机靠b影院| www.av在线官网国产| 精品熟女少妇八av免费久了| xxxhd国产人妻xxx| 免费在线观看影片大全网站| 日本黄色日本黄色录像| 在线观看免费午夜福利视频| 国产免费视频播放在线视频| 欧美亚洲日本最大视频资源| 老汉色av国产亚洲站长工具| 91国产中文字幕| 三上悠亚av全集在线观看| 国产亚洲欧美在线一区二区| 女警被强在线播放| 丁香六月欧美| 99国产极品粉嫩在线观看| 国产成人精品在线电影| 免费不卡黄色视频| 99久久99久久久精品蜜桃| 亚洲少妇的诱惑av| 亚洲性夜色夜夜综合| 99久久精品国产亚洲精品| netflix在线观看网站| 一二三四在线观看免费中文在| 啦啦啦免费观看视频1| 久久久久久久久久久久大奶| 丰满人妻熟妇乱又伦精品不卡| 在线十欧美十亚洲十日本专区| 男女午夜视频在线观看| 欧美日韩亚洲综合一区二区三区_| 日本vs欧美在线观看视频| 午夜91福利影院| 免费黄频网站在线观看国产| 国内毛片毛片毛片毛片毛片| 女人精品久久久久毛片| 国产成+人综合+亚洲专区| 搡老岳熟女国产| 香蕉丝袜av| 在线av久久热| 99久久综合免费| 亚洲欧美色中文字幕在线| 欧美精品人与动牲交sv欧美| 午夜久久久在线观看| 欧美日韩av久久| 国产欧美日韩一区二区三 | 久久久久久久久免费视频了| 欧美亚洲日本最大视频资源| 成人亚洲精品一区在线观看| 国产97色在线日韩免费| 色婷婷久久久亚洲欧美| 一区福利在线观看| 亚洲人成电影观看| 黄色怎么调成土黄色| 成人国产一区最新在线观看| 91av网站免费观看| 久久久久国产一级毛片高清牌| 欧美国产精品va在线观看不卡| 国产野战对白在线观看| 国产成人精品无人区| 久久ye,这里只有精品| 啦啦啦啦在线视频资源| 色94色欧美一区二区| 成人国产av品久久久| 王馨瑶露胸无遮挡在线观看| 新久久久久国产一级毛片| 亚洲七黄色美女视频| 黑人猛操日本美女一级片| 啦啦啦视频在线资源免费观看| 五月天丁香电影| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久精品人妻al黑| 国产精品二区激情视频| 亚洲天堂av无毛| 在线观看免费高清a一片| 美女中出高潮动态图| 黄片大片在线免费观看| 看免费av毛片| 国产深夜福利视频在线观看| 51午夜福利影视在线观看| 欧美国产精品一级二级三级| 性色av一级| 国产成人精品在线电影| 久久久久国产一级毛片高清牌| 国产精品99久久99久久久不卡| 国产欧美日韩一区二区三 | 美女福利国产在线| 涩涩av久久男人的天堂| av网站免费在线观看视频| 成人18禁高潮啪啪吃奶动态图| 亚洲中文字幕日韩| 国产一区二区三区在线臀色熟女 | a 毛片基地| 蜜桃在线观看..| 亚洲人成电影免费在线| 丰满少妇做爰视频| 精品一品国产午夜福利视频| 欧美在线一区亚洲| 国产精品秋霞免费鲁丝片| 亚洲专区字幕在线| 少妇被粗大的猛进出69影院| 亚洲色图 男人天堂 中文字幕| 日日爽夜夜爽网站| 在线观看一区二区三区激情| 永久免费av网站大全| 淫妇啪啪啪对白视频 | 香蕉丝袜av| 正在播放国产对白刺激| 亚洲 国产 在线| 亚洲精品中文字幕一二三四区 | 精品少妇久久久久久888优播| 动漫黄色视频在线观看| 日韩欧美免费精品| 国产高清国产精品国产三级| 十八禁高潮呻吟视频| a级毛片在线看网站| 啦啦啦在线免费观看视频4| 午夜影院在线不卡| 亚洲精品日韩在线中文字幕| 免费在线观看日本一区| 成人国产一区最新在线观看| 国产福利在线免费观看视频| xxxhd国产人妻xxx| 国产成人系列免费观看| 欧美成狂野欧美在线观看| 搡老熟女国产l中国老女人| 欧美黑人精品巨大| 最黄视频免费看| 美女国产高潮福利片在线看| 99国产精品免费福利视频| 国产成人免费观看mmmm| 国产精品免费视频内射| 在线天堂中文资源库| 色综合欧美亚洲国产小说| 亚洲三区欧美一区| 老司机影院毛片| 国产精品欧美亚洲77777| 在线亚洲精品国产二区图片欧美| 亚洲av电影在线进入| 午夜久久久在线观看| 欧美 日韩 精品 国产| 男女之事视频高清在线观看| 中文字幕最新亚洲高清| 人人妻,人人澡人人爽秒播| 1024视频免费在线观看| 91精品国产国语对白视频| av欧美777| 欧美日韩一级在线毛片| 色老头精品视频在线观看| 精品少妇黑人巨大在线播放| 操美女的视频在线观看| 久热爱精品视频在线9| 久久久国产精品麻豆| 欧美国产精品一级二级三级| 狠狠精品人妻久久久久久综合| 美女脱内裤让男人舔精品视频| 亚洲综合色网址| 免费在线观看影片大全网站| 日本a在线网址| 国产99久久九九免费精品| 巨乳人妻的诱惑在线观看| 成人手机av| 色婷婷久久久亚洲欧美| 中文字幕人妻丝袜一区二区| 成人国语在线视频| 日韩中文字幕欧美一区二区| 日韩欧美免费精品| 高潮久久久久久久久久久不卡| 国产深夜福利视频在线观看| 日韩人妻精品一区2区三区| 一级毛片精品| 在线观看免费视频网站a站| bbb黄色大片| 久久国产精品影院| 亚洲av电影在线进入| 91大片在线观看| 97在线人人人人妻| 视频在线观看一区二区三区| 午夜福利影视在线免费观看| 99精品欧美一区二区三区四区| 日韩一卡2卡3卡4卡2021年| 亚洲精品成人av观看孕妇| 国产免费现黄频在线看| svipshipincom国产片| 国产精品影院久久| 成人手机av| 久久久久网色| 亚洲av日韩精品久久久久久密| 国产男女内射视频| 日韩制服丝袜自拍偷拍| 久久精品成人免费网站| 国产亚洲精品一区二区www | 亚洲专区国产一区二区| 天天添夜夜摸| 久久综合国产亚洲精品| 国产xxxxx性猛交| 人人妻人人澡人人看| 国产视频一区二区在线看| 国产有黄有色有爽视频| 亚洲欧美日韩高清在线视频 | 亚洲,欧美精品.| 精品久久蜜臀av无| 亚洲 国产 在线| 国产精品成人在线| 在线 av 中文字幕| 免费高清在线观看日韩| 亚洲情色 制服丝袜| 日韩 欧美 亚洲 中文字幕| 一级片免费观看大全| 伦理电影免费视频| av在线老鸭窝| 久久99一区二区三区| 亚洲精品国产色婷婷电影| 少妇精品久久久久久久| 中文字幕人妻熟女乱码| 亚洲专区国产一区二区| 久久人妻熟女aⅴ| 少妇猛男粗大的猛烈进出视频| 久热这里只有精品99| 日韩制服骚丝袜av| 国产精品一二三区在线看| 男女下面插进去视频免费观看| 一区二区三区四区激情视频| 久久毛片免费看一区二区三区| 国产精品秋霞免费鲁丝片| 一区二区三区激情视频| 亚洲av欧美aⅴ国产| 一本久久精品| 国产三级黄色录像| 成人18禁高潮啪啪吃奶动态图| 无遮挡黄片免费观看| 亚洲欧美成人综合另类久久久| 一级a爱视频在线免费观看| 欧美xxⅹ黑人| 国产成人a∨麻豆精品| 亚洲久久久国产精品| 欧美少妇被猛烈插入视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产欧美日韩综合在线一区二区| 亚洲av成人一区二区三| 久久久久精品人妻al黑| 久久午夜综合久久蜜桃| 日韩欧美一区视频在线观看| 成人国产av品久久久| 久久久久久人人人人人| 18在线观看网站| 男人爽女人下面视频在线观看| 欧美激情极品国产一区二区三区| 国产精品欧美亚洲77777| 亚洲国产精品一区二区三区在线| 母亲3免费完整高清在线观看| 啦啦啦在线免费观看视频4| 18禁国产床啪视频网站| 国产激情久久老熟女| 精品少妇久久久久久888优播| 欧美97在线视频| 国产区一区二久久| 成人国产一区最新在线观看| 国产av又大| 国产亚洲av高清不卡| 国产福利在线免费观看视频| 亚洲欧美清纯卡通| 亚洲性夜色夜夜综合| 精品少妇内射三级| 国产欧美亚洲国产| 人人澡人人妻人| 日韩欧美一区二区三区在线观看 | 国产日韩欧美亚洲二区| 欧美日韩一级在线毛片| 老汉色∧v一级毛片| 黄片大片在线免费观看| 久久国产精品影院| 黄色视频,在线免费观看| 精品少妇内射三级| 性高湖久久久久久久久免费观看| 日本av手机在线免费观看| 亚洲精品一区蜜桃| 女警被强在线播放| 午夜视频精品福利| 成人黄色视频免费在线看| 午夜两性在线视频| 日韩人妻精品一区2区三区| 亚洲午夜精品一区,二区,三区| 老司机影院毛片| 精品高清国产在线一区| 成人免费观看视频高清| 黄片播放在线免费| 精品乱码久久久久久99久播| 91大片在线观看| 1024香蕉在线观看| 欧美精品av麻豆av| 久久精品国产亚洲av香蕉五月 | 欧美日韩中文字幕国产精品一区二区三区 | 亚洲人成电影观看| 午夜激情av网站| 成年美女黄网站色视频大全免费| 正在播放国产对白刺激| 欧美在线黄色| 日韩大片免费观看网站| 一区二区三区乱码不卡18| 91字幕亚洲| 亚洲精华国产精华精| 91av网站免费观看| 91大片在线观看| 久久久久久久国产电影| 最近最新中文字幕大全免费视频| 欧美 亚洲 国产 日韩一| 妹子高潮喷水视频| 美女大奶头黄色视频| 热99国产精品久久久久久7| 国产av精品麻豆| 50天的宝宝边吃奶边哭怎么回事| 老汉色∧v一级毛片| 两性夫妻黄色片| 免费在线观看视频国产中文字幕亚洲 | 免费黄频网站在线观看国产| a级毛片在线看网站| 日日摸夜夜添夜夜添小说| 亚洲精品一二三| 欧美精品av麻豆av| 啦啦啦免费观看视频1| 日韩制服骚丝袜av| 亚洲精品国产区一区二| 久久久久久久久久久久大奶| 又黄又粗又硬又大视频| 国产免费现黄频在线看| 男女高潮啪啪啪动态图| 波多野结衣av一区二区av| 久久久精品国产亚洲av高清涩受| 日韩 亚洲 欧美在线| 日本a在线网址| 国精品久久久久久国模美| 日韩欧美一区视频在线观看| 亚洲成国产人片在线观看| 午夜精品国产一区二区电影| 男女下面插进去视频免费观看| 我要看黄色一级片免费的| 亚洲精品日韩在线中文字幕| 纯流量卡能插随身wifi吗| 精品欧美一区二区三区在线| 亚洲国产精品999| 亚洲av成人一区二区三| av超薄肉色丝袜交足视频| 丝袜在线中文字幕| 老司机午夜福利在线观看视频 | 精品人妻在线不人妻| 一级毛片女人18水好多| 男女之事视频高清在线观看| 日韩中文字幕欧美一区二区| 色婷婷av一区二区三区视频| 在线永久观看黄色视频| 少妇猛男粗大的猛烈进出视频| 免费高清在线观看日韩| 久热爱精品视频在线9| 黄色片一级片一级黄色片| 久久狼人影院| 亚洲欧美日韩高清在线视频 | 亚洲人成电影免费在线| 精品一区二区三卡| 人人妻人人添人人爽欧美一区卜| 精品久久久久久久毛片微露脸 | 免费高清在线观看日韩| 亚洲精品成人av观看孕妇| 人人妻人人澡人人看| 老汉色∧v一级毛片| 91国产中文字幕| 亚洲中文日韩欧美视频| 亚洲精品中文字幕在线视频| 国产成人精品无人区| 亚洲午夜精品一区,二区,三区| 亚洲精品美女久久久久99蜜臀| av片东京热男人的天堂| 丝袜喷水一区| 色综合欧美亚洲国产小说| 午夜福利一区二区在线看| 高清在线国产一区| 日韩电影二区| 亚洲国产av影院在线观看| av免费在线观看网站| 久久99一区二区三区| 我要看黄色一级片免费的| 欧美激情高清一区二区三区| 搡老熟女国产l中国老女人| 成在线人永久免费视频| 中文字幕高清在线视频| 视频区欧美日本亚洲| 精品高清国产在线一区| 国产免费现黄频在线看| 精品卡一卡二卡四卡免费| 90打野战视频偷拍视频| 色播在线永久视频| 黄频高清免费视频| 免费日韩欧美在线观看| 久久天躁狠狠躁夜夜2o2o| 麻豆乱淫一区二区| 在线观看免费日韩欧美大片| 999久久久精品免费观看国产| 亚洲av国产av综合av卡| 亚洲五月色婷婷综合| 国产黄频视频在线观看| 亚洲精品中文字幕在线视频| 国产一区二区在线观看av| 日本五十路高清| 操美女的视频在线观看| 黄片小视频在线播放| 一本大道久久a久久精品| 香蕉国产在线看| 久久国产精品男人的天堂亚洲| 国产精品一区二区在线不卡| 51午夜福利影视在线观看| 飞空精品影院首页| 免费久久久久久久精品成人欧美视频| 99精品欧美一区二区三区四区| 亚洲av成人一区二区三| 如日韩欧美国产精品一区二区三区| 新久久久久国产一级毛片| 高清在线国产一区| 成年美女黄网站色视频大全免费| 秋霞在线观看毛片| 91麻豆av在线| 亚洲精品美女久久久久99蜜臀| 色婷婷久久久亚洲欧美| 欧美日韩一级在线毛片| 午夜福利一区二区在线看| 国产成+人综合+亚洲专区| 麻豆国产av国片精品| 在线永久观看黄色视频| 欧美性长视频在线观看| 免费在线观看黄色视频的| 久久中文看片网|