红层滑坡滑带土蠕变力学特性及经验模型研究

第47卷第3期

2021年3月

水力发电

贾逸12，魏良帅12，黄海峰3

(1.中国地质科学院探矿工艺研究所，四川成都611734 ；

2-中国地质调查局地质灾害防治技术中心，四川成都611734；

3-四川成绵苍巴高速公路有限责任公司成都分公司，四川成都641400)

摘要：土体蠕变本构经验模型由于拟合参数无明确物理意义难以推广#为解决经验模型适用性弱的难题，在传统Mesri经验模型的基础上，保留原有双曲线方程描述应力-应变关系，引入分数阶微积分理论反映土体应变-寸间关系，使经验模型参数物理意义更为明确，得到改进后的蠕变经验模型#利用该模型对西南地区某红层滑坡滑带土和相关文献中的黄土和高含盐细砂土蠕变数据进行辨识对比，证明本模型的可行性和适用性，可为土体蠕变经验模型构建及红层边坡长期稳定性分析提供参考#

关键词：红层；滑带土；蠕变；经验模型；双曲线；分数阶微积分

Study on Crrep Mechanics Characteristics and Empiricai Modei of Soil in the Slip Zone of Red Layer Landslite

JIA Yi1'2,WEI Liangshuai1'2,HUANG Haifeng3

(1.InstotuteofECpaoeatoon Technoaogy,ChoneseAcademyofGeoaogocaaScoences,Chengdu611734,Sochuan,Chona;

2.Technicat Center for Geologicat Hazard Prevention and Control,China Geelogicat Survey,Chengdu611734,Sichuan,China;

3.Chengdu Branch of Sichuan Chengmian Cangba Expressway Co.,Ltd.,Chengdu641400,Sichuan,China)

Abstract:It is difficult to generalize the empiricat creep ccnstitutive model because the fitting parameters have no clear physicat meaning.In order to solve the problem of weak applicability of the empiricat model,based on the traditionat Mesri empoeocaBmode,theoeogonaBhypeebooc equatoon oseetaoned todesceobethe steess-steaon eeatoonshop,and the feactoonaB caacuaustheoeyisinteod

uced toeefaectthesoiasteain time ee aation ship,so that the phy sica ameaning of the empi eicaamode a paeameteesiscaeaeee,and theimpeoeed ceeep empi eica amode ai sobtained.The mode ai su sed toidentifyand compaeetheceeep data of a red bed landslide sliding zone soit and the loess and high salt ccntent fine sand in related literature in Southwest China,which p eo e<sth<f a sibi aity and appaicabiaity of th<mod<a.Th<e<s<aech e<suatscan peoeid<e<fe<nc<s foeth< consteuction ofsoiaceeep empieicaamodeaand theaong-teem stabiaityanaaysisofeed bed saope.

Key Words:red bed;sliding zone soil;creep;empiricat model;hyperbola;fractional calculus

中图分类号：TU411.3文献标识码：A文章编号：0559-9342(2021)03-0025-06

0引言

红层在我国西南地区广泛分布，孕育了大量滑坡体，威胁人类生命财产安全［1-2］#滑带是指在滑坡发育、发生过程中经历物理化学作用，力学性质大幅衰减的下滑界面［3I O当滑带土蠕变累积至某一程度，易在降雨、地震等外界因素的诱发下发生滑坡。因此，研究滑带土蠕变力学特性对边坡长期稳定性分析及灾害防治具有十分重要的意义(4)。对于滑带土蠕变的研究已有较多成果，赖小玲等(51在传统Mesri蠕变经验模型的基础上，将基质吸力作为独立变量，建立考虑基质吸力的滑带土蠕变经验模型；Bhat等(61采用扭剪蠕变仪研究滑带土在残余状态下的蠕变行为；Maia等⑺开展直剪蠕变试验，研

究孔隙水含盐量对滑带土残余蠕变强度的影响；Wen 等(81研究发现含砾量与滑带土蠕变强度之间无直接关联；李梦姿(91构建了滑带土非线性Burgers模型；

收稿日期:2020-09-02

基金项目：中国地质调查局二级项目(DD20190325)

作者简介：贾逸(1987—)，男，四川成都人，工程师，硕士，研究方向为岩土体稳定性及工程环境效应；魏良帅(通讯作者)•

Watrr Power Vol.47No.

水力发电2021 年 3 月

黄海峰等（101通过Log-Modified 函数描述应变-时间关

系，建立新的滑带土蠕变经验模型；蒋树等1引入非线性黏塑性损伤元件，改进传统Cvisc 模型，对滑带土蠕变行为模拟取得良好成效；朱荣森1以黄土滑坡滑带土为研究对象，开展三轴压缩蠕变试验，

得到蠕变力学性质与含水率和围压之间的规律#

土体蠕变本构模型以针对性较强的经验模型和适用性较强的元件模型为主#经验模型由于拟合参数无明确物理意义难以推广#为解决经验模型适用

性弱的难题，本文在传统Mesri 经验模型的基础上，

保留原有双曲线方程描述应力-应变关系，引入分数阶微积分理论反映土体应变-时间关系，使经验模型

参数物理意义更为明确，得到改进后的蠕变经验模型#利用该模型对西南地区某红层滑坡滑带土和相 1.2试验方法

首先开展滑带土三轴不排水剪切试验，得到不排水剪切强度由于滑带在滑坡体内所处位置较

浅，实际围压较低，故将围压设置为100、200 kPa

和300 kPa ，应力水平；取为0.5、0.55、0. 6、 0.65、0. 7和0.75，则每级偏应力对应为0.5'、

0.55'、0.6'、0.65'、0. 7'和 0.75应力路径

示意见图1 （以围压100 kPa 为例）。蠕变试验加载方案见表2 #蠕变试验首先以0. 2 kPa/min 的速率将围

压加载至预定值，接着固结土样，固结时间设定为 4 h ，之后再以0. 05 mm/min 的剪切速率施加剪切荷

载，每0. 2 mm 的剪切位移记数1次#

关文献中的黄土和高含盐细砂土蠕变数据进行辨识

对比，证明本模型的可行性和适用性，可为土体蠕变经验模型构建及红层边坡长期稳定性分析提供参考#

1滑带土蠕变试验

-----应力水平

围压

100806040

1 - 1试样制备及试验仪器

试样取自西南地区某红层滑坡滑带，由于钻孔取样时不可避免有一定扰动，取原状土条件受限,

故本文以重塑样进行固结不排水三轴压缩蠕变试验#

土样密封运回试验室风干后将其碾散，过1mm 筛,

通过击实法并以天然含水率12%制备直径39. 1 mm 、高度80 mm 的圆柱体试样，最大允许差值在1%以

内，符合GB/T 50123&2019《土工试验方法标准》的

规定，用保鲜膜充分包裹，静置1 d 使试样水分均

匀# 土样物理力学参数见表1 #

表1 土样物理力学参数

干密度p/-3

g* cm 液限%/%塑限!p/%

黏聚力c/kPa 内摩擦角

）/（

。

1.6631.6

15.28.37

16.45

试验仪器为CSS-2901TS 土体三轴流变试验机# 该试验机由轴向加载系统、围压加载系统和变形测

量系统等组成，最大轴向荷载达10 kN ，围压及孔隙水压力变化范围皆为0〜2 MPa ，轴向位移量测范围为 0 ~ 20 mm #

时间/h 0.4————1 00

160240320

400

图1应力路径

2试验成果

2. 1滑带土蠕变曲线

通过图1所示应力路径开展固结不排水三轴蠕

变试验，得到分级加载蠕变曲线，见图2 ；通过 Boltzmann 线性叠加（切处理图2，绘制分别加载蠕变曲线，见图3 #图2、3中，将各级应力水平标识于

曲线上#从图2可看出，分级加载蠕变曲线为“台阶式”攀爬的曲线，每一级“台阶”高于前一级#高围

压下的应变量值总是大于低围压#从图3可看出，土样在加载瞬间首先表现出瞬弹性变形，紧接着是

历时较短的衰减蠕变阶段，该阶段曲线斜率逐渐降

低，蠕变速率递减，随后进入稳定蠕变阶段，此时蠕变速率基本保持平衡，应变量值基本保持恒定#

土样轴向应变由瞬时应变和蠕变应变组成，蠕变应变包含衰减、稳定蠕变阶段#围压300 kPa 时，第6

kPs

表2蠕变试验加载方案

每一级加载的偏应力

试样编号最小王应力$剪切强度 '

123456CU-2100692.25346

381

415

450484—

CU-3 200 851.42

426 468 511 553

596 638CU-4

300

1 045. 17 523

575

627 679

732

784

] Watem Power Vol. 47 No.

第47卷第3期贾逸，等：红层滑坡滑带土蠕变力学特性及经验模型研穷

级加载下稳定蠕变阶段的曲线斜率呈明显增大趋势，

说明蠕变变形现象更为显著#对于同一滑带土，随

着土体垂直压力的增加，围压不断增大，轴向应变

值递增#

图2分级加载蠕变

2.5

0.65

------▲

0.7

~a--------------a-----0.6•-----*------7x0(5

-------

0.5

0204060

时间Zh

a围压100kPa

时间Zh

图3分别加载蠕变

2.2滑带土的等时偏应力-应变特性

从图3中择取1、11、21、31、41、51、61h和71h共8个时间节点的偏应力和应变值，绘制等时偏应力-应变曲线，见图4#从图4可以看出，滑带土等时偏应力-应变曲线近乎为一束曲线簇，其形态近似为双曲线，具有明显的非线性特征#蠕变时间为1h的曲线与11、21、31、41、51、61h和71h 这7条曲线组成的曲线簇明显分离，这说明土样在初始荷载下存在明显的弹性变形，第1级加载情况下尤为显著#分析可知，滑带土总应变由瞬时弹性应变和蠕变应变组成#

煤

辱

c围压300kPa

图4等时偏应力-应变关系

3蠕变经验模型

3.1经验型蠕变模型的建立

滑带土总应变由瞬时应变和蠕变应变组成，暂不考虑瞬时应变，将滑带土蠕变应变+表示为

+==1(7)Z(t)(1%式中，a(7、a(r分别为土样的应力-应变关系函数、应变-时间关系函数；7、t分别为函数A(7)、

Water Power Vol.47No.3|

水力发电

2021年3月

A ( #的自变量#

3.2应力-应变关系

由于滑带土等时偏应力-应变曲线近似双曲线，参考传统MesC 经验模型(141引入Kondnar 双曲线型偏

应力-应变方程，即

式中，—(—为伽马函数，-为大于0的分数阶数；s 为用于拉普拉斯变换的某一自变量#函数/( #)的-

阶微分相应地定义为

d - [/(#] = £[ ' ( t - s ) ”--d t d t # J0 — ( n -—)

式中,$、$分别为最大和最小主应力；a 、6为双

曲线方程参数。

对式(2 )等号两侧同时求微分，得到初始切线模量E u 为

当+$s 时，对式(2 )取极限值即为最终偏应力

差($ - $ ) ull

!即

($ - $ ) ui = lim — = +

( 4 )

+$ s a + b 彳 b

实际上，土样应变不可能无穷大，当外界荷载达到排水剪切强度'时土样开始破坏#为了在

Kondnar 双曲线中体现破坏点[+ , ( $-$)『]，其

中,l 为破坏点对应的应变,($ - $力为破坏偏应力#引入破坏比R A 1 ，公式为

( $1 - $3 ) A ( $1 - $3 ) A

R = ( $-$) ui = 1/C

)

将式(3 )、(5 )代入式(2 )得到应力-应变本构关系式, 即

式中，； = ($-$)/( $-$ ) f o

3.3应变-时间关系

蠕变经验模型中，应变-时间关系多通过双曲线

函数、幂次函数和对数函数甚至生长函数等(10,151多

种形式进行描述，尽管能取得一定的拟合效果，但会导致应力-应变关系参数无实际物理意义，难以表

征蠕变发展过程的内在性质，且由于其具有较强的针对性，无法进行推广，故本文引入分数微积分理

论反映土样应变-时间关系。

分数阶微积分的核心在于其阶数为有理分数、无理数和复数，可表征与某一现象的整个发展历史

有关的性质[161#分数阶微积分物理意义清晰，结构简单，在非线性动力系统中优势显著#滑带土是一种非线性的地质材料，故引入应用较为广泛的

Riemann-Liouvilla 型分数阶微积分#依据该理论，函数 /( t )在可积区间[0 , 11 的(阶 Riemann -Liouvilla 积分定义为

D ：代 t ) :

(r r1 = -(-v '

( t - E (-f ( s ) d s ( 7 )

d r

— (-) J 。

式中，n = [ —1# —(—)定义为

—(—)=

t a -r -dt Ra ( —) > 0 (9)

在滑带土蠕变发展过程中，土样蠕变应变一般用牛顿体本构方程描述，即

$ = - l

( 10 )

式中，-为黏滞系数；l 为蠕变应变率；为便于描述，将偏应力($-$)记为$#

基于分数阶微积分理论, 可将该本构关系其转化为

式中，-为分数阶黏滞系数#

当应力$恒定时，依据R-iemann-Liouvilla 积分

理论, 对式(7) 进行 Romann-Looueo a 分数阶积分得

l =■$

_7，0 "—"】

(12 )

—厂(1 + —

式(12 )即为土样应变-时间关系#取$为200

kPa , -%为5 X105 kPa-h ，通过式(12 )绘制不同—值

的蠕变曲线，见图5 #从图5可看出，当—增大时, 曲线斜率不断增加，曲线逐渐远离时间轴，表现出明显的非线性特征；当—增至1时，表现出完全线

性关系#

图5不同a 值的蠕变

3.4经验型蠕变模型

将式(6)、 (12)同时代入式(1)可得

£ =

($ - $)/. —D . $—t — (13)

2 -

E u

1 - R ； - —(1 + —) 丿

瞬时弹性应变1服从Hooka 定律，i = $/E u ，由于式(13 )未考虑瞬时应变，故将式(13 )改为

Wat€% Power Vol. 47 No.

第47卷第3期贾逸，等：红层滑坡滑带土蠕变力学特性及经验模型研究

$($-$)/；$#

+=+++=D+瓦1'-,(1+()

(14)

由于；=($1-$3)/($1-$3)f，将其代入式(14)化简可得

+=-$+_______________1_____#____(15)

D,D(1-?；)-，(1+('丿式(15)即为本文所建蠕变经验模型#

4模型验证及参数求解

传统经验模型皆是分别求解应力-应变关系和应变-时间关系参数，实际土体蠕变发展过程中变形同时受到应力和时间的影响，通过应力-应变-时间模型整体同时求解参数更为合理。本文建立的蠕变经验模型有D,、？、-和(共4个参数，通过一般的非线性最小二乘算法求解即可。模型参数见表3。引用文献:14

］中传统Mesri经验模型与所建模型进行对比，理论值与试验值对比曲线见图6。

从图6可看出，引用经验模型对稳定蠕变阶段的拟合偏差较大，拟合精度相对较低，平均相关系数？仅有0.9476#而本文建立的模型对衰减蠕变阶段和稳定蠕变阶段的辨识更为精准，平均相关系数达到0.9839#传统Mesri经验模型属于传统经验模型，具有很强的针对性，仅针对某一种土体，难以推广。而本文建立的模型在传统经验模型的基础上，通过分数阶微积分描述应变-时间关系，从而得到新的分数阶经验模型。该模型参数具有一定的物理意义，模型内涵较为丰富，可从某种意义上表征土体蠕变过程的内在性质，具有一定的推广价值#

为验证本文所建模型描述土体蠕变力学性质的适用性，引用文献［17-18］中黄土、高含盐细砂土的蠕变数据，利用本文所建模型进行模拟辨识，得到的理论值与试验值对比曲线见图7。从图7可看出，本文模型对黄土、高含盐细砂土蠕变曲线具有较强的辨识能力，平均相关系数？分别为0.9915和0.9863，理论值和试验值吻合较好。由此证明，本文建立的经验模型能较好地描述土体蠕变力学行为，具有一定的适用性。

由于元件模型的参数多具有明确物理意义，辨识不同种类、区域的土体蠕变时每个参数表现出不同的规律，故土体元件模型在蠕变模拟中应用更为广泛。经验模型一般仅针对某种土体，本文在传统Mesri经验模型的基础上，保留原有双曲线方程描述应力-应变关系，引入分数阶微积分理论反映土体应变-时

间关系，使得改进后的经验模型具有一定的物理意义，突破传统经验模型拟合参数无意义的限制，从某种程度上能体现土体蠕变过程的内在性质，为蠕变经验模型的深入研究提供一种新思路。但是，经验模型如何借鉴元件模型的优点，通过一组或多组参数的差异反映土体性质的规律性差异，这是接下来研究的重点。

5结语

滑带土蠕变变形过程中，应变包含瞬时弹性应

表3模型参数

$3skPa($i-$3)/kPa D“/kPa—%/MPa-h?a相关系数？

100 200 3003460.5499.4268.74 1.540.080.9895 381

415

450

484

426

468

511

553

596

638

523

575

627

679

732

784

0.55

0.6

0.65

0.7

0.5

0.55

0.6

0.65

0.7

0.75

0.5

0.55

0.6

0.65

0.7

0.75

1214.92

1387.65

3641.53

3236.78

708.89

1231.58

1786.52

3950.03

4264.78

3652.84

816.79

1462.73

1985.96

2836.50

3872.54

3964.28

51.36

44.28

46.85

31.90

80.45

72.33

54.19

42.36

38.36

35.47

87.49

68.55

57.42

51.39

42.56

38.76

1.23

1.06

0.73

0.59

1.36

1.18

0.95

0.78

0.56

0.51

1.14

0.87

0.72

0.55

0.38

0.32

0.10

0.13

0.18

0.25

0.06

0.12

0.17

李梦个人资料简介

0.15

0.18

0.22

0.03

0.08

0.14

0.16

0.20

0..23

0.9876

0.9867

0.9838

0.9745

0.9843

0.9884

0.9822

0.9861

0.9853

0.9786

0.9818

0.9855

0.9795

0.9868

0.9848

0.9810

Wa#Power Vol47No.

红层滑坡滑带土蠕变力学特性及经验模型研究

发布评论取消回复

最近发表

热门文章

标签列表