易長(zhǎng)圭

(湘潭大學(xué) 土木工程與力學(xué)學(xué)院, 湖南 湘潭, 411105)

砂土微宏觀性質(zhì)相依關(guān)系研究

易長(zhǎng)圭

(湘潭大學(xué) 土木工程與力學(xué)學(xué)院, 湖南 湘潭, 411105)

為了研究砂土微觀力學(xué)性質(zhì), 采用了離散單元法對(duì)砂土進(jìn)行三軸試驗(yàn)?zāi)M, 并引入抗轉(zhuǎn)動(dòng)模型來(lái)表征砂土顆粒的抗轉(zhuǎn)動(dòng)效應(yīng)。通過(guò)對(duì)選取的不同法向接觸剛度、剛度比、顆粒內(nèi)摩擦角、滾動(dòng)剛度系數(shù)以及塑性彎矩系數(shù)進(jìn)行三軸排水?dāng)?shù)值試驗(yàn), 得到了砂土不同微觀參數(shù)與其宏觀力學(xué)性質(zhì)的對(duì)應(yīng)關(guān)系。試驗(yàn)結(jié)果表明:法向接觸剛度對(duì)砂土彈性模量影響顯著, 并呈現(xiàn)正相關(guān)線性關(guān)系; 剛度比與砂土泊松比呈現(xiàn)較為顯著的負(fù)相關(guān)線性關(guān)系; 顆粒內(nèi)摩擦角對(duì)砂土宏觀內(nèi)摩擦角影響顯著, 滾動(dòng)剛度系數(shù)影響砂土殘余強(qiáng)度, 塑性彎矩系數(shù)主要影響峰值強(qiáng)度。

離散單元法; 砂土; 微觀參數(shù); 宏觀參數(shù)

離散單元法是上個(gè)世紀(jì)70年代由Cundall[1]提出的一種基于非連續(xù)介質(zhì)力學(xué)的數(shù)值計(jì)算方法。其基本思想是把所研究材料離散為單個(gè)單元, 各個(gè)單元運(yùn)動(dòng)由牛頓第二運(yùn)動(dòng)定律控制。由于離散單元法在顆粒尺度上對(duì)物理量變化進(jìn)行考量, 因此, 其得到了廣泛應(yīng)用。在巖土工程領(lǐng)域內(nèi), 離散單元法已經(jīng)成功應(yīng)用于砂土剪切[2]、結(jié)構(gòu)性土壓縮試驗(yàn)[3]、非飽和土吸力[4]等研究中。離散單元法計(jì)算原理雖然簡(jiǎn)單, 但由于接觸模型參數(shù)屬于顆粒尺度參數(shù), 這些細(xì)觀參數(shù)通常是未知的, 又無(wú)法通過(guò)試驗(yàn)直接獲取, 所以選取一個(gè)合理的計(jì)算參數(shù)十分困難, 而在進(jìn)行離散單元法模擬時(shí), 只有細(xì)觀計(jì)算參數(shù)正確的前提下, 才能獲得一個(gè)合理的計(jì)算結(jié)果[5]。目前, 離散元微觀參數(shù)的選取方式多種多樣, 目前大多數(shù)研究者采用試錯(cuò)法[6]。針對(duì)微觀參數(shù)選取的研究主要集中在細(xì)觀顆粒力學(xué)分析以及宏觀試驗(yàn)推導(dǎo)兩個(gè)方面。武力等[7]將遺傳神經(jīng)網(wǎng)絡(luò)與三軸試驗(yàn)離散元數(shù)值模擬結(jié)合起來(lái), 并進(jìn)行改性砂土顆粒離散元接觸模型參數(shù)反演;邢紀(jì)波等[8]根據(jù)應(yīng)力傳播條件給出了顆粒間接觸剛度的理論計(jì)算公式, 并給出了求解靜態(tài)問(wèn)題的自適應(yīng)阻尼; 王泳嘉等[9]依據(jù)兩彈性球受集中力作用情況, 得到了法向和切向剛度的理論公式; 此外,Mindlin[10]、Hentz[11]、張振南等[12]也通過(guò)相關(guān)研究得到了接觸剛度表達(dá)式; 徐小敏[13]等通過(guò)室內(nèi)三軸試驗(yàn)進(jìn)行回歸分析得到了宏微觀參數(shù)之間的經(jīng)驗(yàn)公式; 焦紅光等[14]通過(guò)離散元模擬指出接觸剛度與計(jì)算時(shí)步存在較大的關(guān)聯(lián); Alessio等[15]提出了巖土材料宏微觀結(jié)合的微觀本構(gòu)模型, 并采用該模型對(duì)凝灰?guī)r的本構(gòu)模型參數(shù)進(jìn)行了預(yù)測(cè)。

綜上所述, 目前針對(duì)離散元微觀參數(shù)和其宏觀力學(xué)特征直接相依關(guān)系的研究較為少見(jiàn), 本文運(yùn)用大型離散元程序SDEC進(jìn)行了大量的三軸試驗(yàn)?zāi)M, 并將適用于模擬砂土力學(xué)性質(zhì)的抗轉(zhuǎn)動(dòng)模型引入進(jìn)程序中。針對(duì)模型中關(guān)鍵微觀參數(shù)進(jìn)行分析, 研究了剛度比、法向接觸剛度、顆粒間摩擦角、滾動(dòng)剛度系數(shù)以及塑性彎矩系數(shù)對(duì)砂土宏觀力學(xué)性能的影響, 并詳細(xì)分析了各個(gè)微觀參數(shù)在宏觀上所對(duì)應(yīng)的力學(xué)響應(yīng), 得到了各個(gè)微觀參數(shù)與砂土宏觀參數(shù)的相依規(guī)律。

1 DEM抗轉(zhuǎn)動(dòng)模型

顆粒的不規(guī)則形狀是影響砂土力學(xué)性能的主要因素。為充分表征砂土的力學(xué)性能, 在線性接觸模型的基礎(chǔ)上, 將Blheine所提出的抗轉(zhuǎn)動(dòng)模型引入進(jìn)該數(shù)值計(jì)算中。如圖1所示, 顆粒之間的相互作用受到法向接觸、切向接觸以及轉(zhuǎn)動(dòng)接觸三個(gè)部分共同作用。

圖1 顆?？罐D(zhuǎn)動(dòng)接觸模型

顆粒間法向接觸力Fn和切向接觸力Ft分別由2個(gè)方向上的顆粒相對(duì)位移及接觸剛度控制, 計(jì)算公式分別為Fn= KnUn和ΔFt= KsdUs。式中: Un為法向位移; dUs為切向位移增量; Kn和Ks分別為顆粒間法向接觸剛度和切向接觸剛度。當(dāng)顆粒間接觸力Fn和Ft滿足關(guān)系Ft≥ Fntanμ時(shí), 顆粒之間不產(chǎn)生滑動(dòng), μ為顆粒間摩擦角。

圖2給出了接觸模型不同部分的力學(xué)響應(yīng)關(guān)系。從圖2(c)可以看出, 轉(zhuǎn)動(dòng)彎矩逐漸從零增長(zhǎng)到最大值, 其值與相對(duì)轉(zhuǎn)動(dòng)角呈現(xiàn)出一定的關(guān)系, 具體計(jì)算式為ΔMr= krΔθr。式中, kr= βKsr2為顆粒滾動(dòng)剛度,其中β則表示無(wú)量綱滾動(dòng)剛度系數(shù), r表示接觸顆粒的平均半徑, Δθr為兩接觸顆粒間單位時(shí)間內(nèi)的相對(duì)轉(zhuǎn)動(dòng)角。顆粒間最大轉(zhuǎn)動(dòng)彎矩可由方程Mp= ηr|Fn|計(jì)算, 式中, η為無(wú)量綱塑性彎矩系數(shù)。

圖2 接觸模型力學(xué)特征

2 模型驗(yàn)證

為驗(yàn)證該離散元模型正確性, 采用抗轉(zhuǎn)動(dòng)模型對(duì)Labenne砂土進(jìn)行三軸試驗(yàn)?zāi)M, 并與相關(guān)文獻(xiàn)的試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比分析。Labenne砂土物理性質(zhì)如表1所示。

三軸試驗(yàn)采用開(kāi)源離散元程序SDEC進(jìn)行模擬, 試樣為50 mm × 50 mm × 50 mm立方體, 試樣由8 120個(gè)顆粒組成(圖3(a)), 邊界為六面剛性墻體, 在固結(jié)過(guò)程中, 六面剛性墻體往試樣中心移動(dòng)直到達(dá)到指定圍壓。剪切過(guò)程中, 周邊墻體通過(guò)伺服控制其應(yīng)力不變, 以達(dá)到穩(wěn)定圍壓的狀態(tài), 而上下墻體同時(shí)以一定速率向中心移動(dòng)進(jìn)行剪切。為充分模擬Labenne砂土性質(zhì), 首先在100 kPa圍壓下對(duì)接觸模型微觀參數(shù)進(jìn)行標(biāo)定, 得到合理微觀參數(shù)如表2所示。然后采用所標(biāo)定參數(shù), 進(jìn)行200 kPa圍壓下三軸試驗(yàn)?zāi)M, 并與室內(nèi)試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比, 結(jié)果如圖3(b)、3(c)所示。

表1 Labenne砂土物理性質(zhì)

表2 接觸模型參數(shù)

由圖3(b)可知, 在不同圍壓下, 數(shù)值試驗(yàn)應(yīng)力—應(yīng)變曲線和室內(nèi)試驗(yàn)曲線吻合良好, 所采用抗轉(zhuǎn)動(dòng)模型能較為準(zhǔn)確地描述砂土應(yīng)力—應(yīng)變的非線性關(guān)系。從圖3(c)可知, 在固結(jié)排水狀態(tài)三軸試驗(yàn)中, 采用抗轉(zhuǎn)動(dòng)模型可反映試樣剪切過(guò)程中的先剪縮后剪脹現(xiàn)象, 能準(zhǔn)確描述剪切過(guò)程中試樣體積變化。

圖3 模型試驗(yàn)及試驗(yàn)結(jié)果對(duì)比

3 三軸試驗(yàn)設(shè)計(jì)

為研究砂土中微觀力學(xué)性能對(duì)其宏觀力學(xué)性質(zhì)的影響, 采用離散單元法進(jìn)行不同微觀參數(shù)下的三軸試驗(yàn)?zāi)M。模擬所采用圍壓均為200 kPa, 顆粒密度為2650 kg/m3, 時(shí)步長(zhǎng)為10-7s。所選取微觀參數(shù)分別為法向接觸剛度、剛度比、顆粒內(nèi)摩擦角、塑性彎矩系數(shù)和滾動(dòng)剛度系數(shù)。具體實(shí)施方案如表3所示。

表3 三軸數(shù)值試驗(yàn)方案

4 結(jié)果與分析

4.1 微觀彈性參數(shù)

圖4 微觀彈性參數(shù)宏觀響應(yīng)

維持其他微觀參數(shù)不變的情況下, 分別改變剛度比和法向接觸剛度進(jìn)行三軸剪切試驗(yàn), 得到了剛度比與泊松比以及法向接觸剛度與彈性模量之間的關(guān)系, 結(jié)果如圖4所示。

由圖4(a)可知, 隨著剛度比從0.1增長(zhǎng)至0.6, 砂土泊松比從0.23下降至0.13, 剛度比與砂土泊松比呈現(xiàn)出明顯的負(fù)相關(guān)線性關(guān)系。此外, 從圖4(b)可知, 隨著法向接觸剛度從1 × 108增長(zhǎng)到1 × 109N/m,試樣彈性模量從156 kPa增長(zhǎng)至907 kPa, 且顆粒間法向接觸剛度和試樣彈性模量呈現(xiàn)出良好的正相關(guān)線性關(guān)系。

4.2 不同微觀強(qiáng)度參數(shù)下應(yīng)力—應(yīng)變關(guān)系

通過(guò)整理試驗(yàn)過(guò)程中不同顆粒內(nèi)摩擦角、不同滾動(dòng)剛度系數(shù)以及不同塑性彎矩系數(shù)下的試樣強(qiáng)度變化曲線, 得到相關(guān)偏應(yīng)力—應(yīng)變關(guān)系如圖5所示。從圖5可以看出, 顆粒間內(nèi)摩擦角以及2個(gè)抗轉(zhuǎn)動(dòng)系數(shù)均對(duì)峰值強(qiáng)度有一定影響, 而僅滾動(dòng)剛度系數(shù)對(duì)試樣殘余強(qiáng)度產(chǎn)生影響。

由圖5(a)可知, 隨著顆粒間摩擦系數(shù)的增大, 峰值應(yīng)力逐漸增大, 但顆粒間摩擦角增長(zhǎng)至40°后, 隨著摩擦角增大,峰值應(yīng)力上升趨勢(shì)逐漸減慢。由圖5(b)可知, 隨著滾動(dòng)剛度系數(shù)增大, 峰值強(qiáng)度和殘余強(qiáng)度均呈現(xiàn)增大趨勢(shì), 但當(dāng)滾動(dòng)剛度系數(shù)大于1.0時(shí), 峰值強(qiáng)度以及殘余強(qiáng)度幾乎不再變化。由圖5(c)可知, 塑性彎矩系數(shù)從0.1增長(zhǎng)至0.3時(shí), 峰值強(qiáng)度從394 kPa增長(zhǎng)至480 kPa, 而隨著塑性彎矩系數(shù)繼續(xù)增大至1.0后, 峰值強(qiáng)度不再產(chǎn)生變化。

圖5 應(yīng)力應(yīng)變關(guān)系曲線

4.3 不同微觀強(qiáng)度參數(shù)下的體變

通過(guò)整理試驗(yàn)過(guò)程中不同顆粒內(nèi)摩擦角、不同滾動(dòng)剛度系數(shù)以及不同塑性彎矩系數(shù)下的試樣體積變化曲線, 得到如圖6所示結(jié)果。

由圖6(a)可知, 顆粒間摩擦角對(duì)體積變化影響顯著, 尤其體現(xiàn)在后期剪脹階段。摩擦角從20°增至40°時(shí), 試樣最大體變值從7.1%增至9.4%。而顆粒間摩擦角繼續(xù)增長(zhǎng)至50°時(shí), 體變值僅增加了0.6%。

圖6 體變關(guān)系曲線

從圖6(b)可以看出, 滾動(dòng)剛度系數(shù)對(duì)體積變化的影響較小, 僅在后期剪脹階段存在一定影響。隨著滾動(dòng)剛度系數(shù)從0.1增長(zhǎng)至1.0, 最大體變量從9.7%增長(zhǎng)至10.7%, 但當(dāng)滾動(dòng)剛度系數(shù)繼續(xù)增加至3.0時(shí), 最大體變量幾乎不再變化。由圖6(c)可知, 塑性彎矩系數(shù)幾乎不對(duì)試樣的體變產(chǎn)生影響。

5 結(jié)論

通過(guò)采用離散單元法對(duì)砂土性質(zhì)進(jìn)行表征, 并采用室內(nèi)試驗(yàn)與數(shù)值試驗(yàn)相結(jié)合的方法對(duì)所建立模型進(jìn)行驗(yàn)證, 建立了符合砂土力學(xué)性質(zhì)的離散元模型。通過(guò)一系列三軸排水試驗(yàn)?zāi)M, 確定模型中微觀參數(shù)與砂土宏觀力學(xué)性質(zhì)之間的對(duì)應(yīng)關(guān)系, 得到如下主要結(jié)論: (1) 抗轉(zhuǎn)動(dòng)模型能較好地描述砂土非線性應(yīng)力—應(yīng)變關(guān)系以及剪縮、剪脹等體積變化特征。(2) 剛度比對(duì)泊松比呈現(xiàn)較為明顯地影響, 并隨著剛度比的增大, 泊松比逐漸降低,兩者呈現(xiàn)明顯線性關(guān)系。(3) 法向接觸剛度對(duì)彈性模量影響較大, 隨著法向接觸剛度增大, 彈性模量逐漸增大, 且兩者呈現(xiàn)較好的線性關(guān)系。(4) 顆粒內(nèi)摩擦角對(duì)顆粒宏觀摩擦角影響較大; 滾動(dòng)剛度系數(shù)對(duì)殘余強(qiáng)度有較為明顯地影響; 塑性彎矩系數(shù)主要影響試件應(yīng)力峰值強(qiáng)度。

[1] Cundall P A. A Computer model for simulating progressive large scale movements in blocky rock system [C]// Proceedings of the Symposium of the International Society of Rock Mechanics. Nancy, France, 1971: 11-19.

[2] Thornton C. Numerical simulations of deviatoric shear deformation of granular media [J]. Geotechnique, 2003, 50(1):43-53.

[3] Jiang Ming-jing, Yu H S, Harris D. Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses [J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2016, 30(7): 723-761.

[4] Liu S H, Sun D A. Simulating the collapse of unsaturated soil by DEM [J]. Int J Numer Anal Meth Geomech, 2012, 26:633–46.

[5] 周健, 賈敏才. 土工細(xì)觀模型試驗(yàn)與數(shù)值模擬[M]. 北京: 科學(xué)出版社, 2008.

[6] Manuel J M M, Luis E M R. Discrete numerical model for analysis of earth pressure balance tunnel excavation [J]. Journal of geotechnical and geoenvironmental engineering, 2015, 131(10): 1 234–1 242.

[7] 武力, 屈福政, 李守巨. 土壓平衡盾構(gòu)改性砂土離散元模型參數(shù)反演方法研究[J]. 大連理工大學(xué)學(xué)報(bào), 2013, 50(6):860–866.

[8] 邢紀(jì)波, 俞良群, 張瑞豐, 等. 離散單元法的計(jì)算參數(shù)和求解方法選擇[J]. 計(jì)算力學(xué)學(xué)報(bào), 1999, 16(1): 47–51, 99.

[9] 王泳嘉, 邢紀(jì)波. 離散單元法及其在巖土力學(xué)中的作用[M]. 沈陽(yáng): 東北工學(xué)院出版社, 1991.

[10] Mindlin R D. Compliance of elastic bodies in contact [J]. Journal of Applied Mechanics, 1949, 71(3): 258–269.

[11] Hentz S, Daudevile L, Donze F V. Identification and validation of a discrete element model for concrete [J]. Journal of Engineering Mechanics, 2012, 130(6): 709–719.

[12] 張振南, 葛修潤(rùn), 李永和. 基于虛內(nèi)鍵理論的材料多尺度力學(xué)模型[J]. 固體力學(xué)學(xué)報(bào), 2006, 27(4): 325–329.

[13] 徐小敏, 凌道盛, 陳云敏, 等. 基于線性接觸模型的顆粒材料細(xì)–宏觀彈性常數(shù)相關(guān)關(guān)系研究[J]. 巖土工程學(xué)報(bào),2010, 32(7): 991–998.

[14] 焦紅光, 李靖如, 趙繼芬, 等. 關(guān)于離散元法計(jì)算參數(shù)的探討[J]. 河南理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2007, 26(1):88–93.

[15] Alession N, Bernhard S, Marek L. Application of artificial neural network for identification of parameters of a constitutive law for soils [J]. Lecture Notes in Artificial Intelligence, 2013, 2 718(2): 545–554.

(責(zé)任編校: 江河)

Research on relationships between microscopic and macroscopic properties of sand

Yi Changgui
(College of Civil Engineering and Mechanics, Xiangtan University, Xiangtan 411105, China)

To investigate the microscopic mechanical properties of sand, a series of consolidated drained triaxial test is performed via discrete element method. In addition, the rolling resistance contact model is added into simulations.Five different microscopic parameters are discussed, namely, the normal contact stiffness, the stiffness ratio, the internal friction angle, the coefficient of rolling stiffness, as well as the coefficient of plastic moment. Finally, the relationship between microscopic parameters and macroscopic properties of sand is obtained. The results show that the normal contact stiffness shows a positive linear correlation with the Young’s Modulus, and the stiffness ratio shows a negative linear correlation with the Possion’s ratio, and the internal friction angle has a significant influence on friction angle; the coefficient of rolling resistance is chosen to control the residual strength; the coefficient of plastic moment play a major role in peak stress.

discrete element method; sand; microscopic parameters; macroscopic parameters

TU 433

: A

1672–6146(2017)03–0064–05

10.3969/j.issn.1672–6146.2017.03.014

易長(zhǎng)圭, yichangguixtu@163.com。

: 2017–04–13