土木建筑与环境工程  2018, Vol. 40 Issue (6): 39-45   PDF    
不同应力路径下旋转硬化方程的统一表示及验证
张浩1,2, 吴志鸿1, 张峰1, 戴自航2    
1. 中建海峡建设发展有限公司, 福州 350003;
2. 福州大学 土木工程学院, 福州 350116
收稿日期:2017-12-06
作者简介:张浩(1980-), 男, 博士, 高级工程师, 主要从事数值模拟及地下工程施工管理研究, E-mail:zhanghao08022@163.com
戴自航(通信作者), 男, 教授, 博士生导师, E-mail:1062573585@qq.com
摘要:在收集到的大量试验数据的基础上,结合旋转硬化表达式,提出一个适用于不同应力路径下的旋转硬化表达式,该公式中共有两个参数,通过理论分析得到这两个参数的确定方法。将该公式导入修正剑桥模型来模拟土体在不同加载路径下的应力应变行为,将模拟结果与试验结果进行对比,结果显示,采用该旋转硬化表达式的各向异性本构模型得到的结果与试验结果均吻合较好,从而验证了该旋转硬化表达式的合理性。
关键词旋转硬化    应力路径    统一表示    本构模型    
Unified expression and verification of rotational hardening equations under different stress paths
Zhang Hao1,2, Wu Zhihong1, Zhang Feng1, Dai Zihang2    
1. Strait Construction and Development Co., Ltd, Fuzhou 350003, P. R. China;
2. College of Civil Engineering, Fuzhou University, Fuzhou 350116, P. R. China
Received: 2017-12-06
Author brief: Zhang Hao(1980-), PhD, main research interests:numerical simulation and construction management of underground engineering, zhanghao08022@163.com.
Dai Zihang(corresponding author), professor, doctoral supervisor, E-mail: 1062573585@qq.com.
Abstract: Based on the collected series of experimental data and many well-known rotational hardening expressions, a rotational hardening expression for different stress paths is proposed, mainly including 2 parameters which can be determined in a direct way through analysis. The proposed rotational hardening expression has been implemented to the MCC model and adopted to represent the stress-strain behavior of clays under different loading paths. The simulations are compared with the experimental results. The comparisons show that the simulated results with this proposed method have good agreement with the experimental results, verifying the rationality of the proposed rotational hardening expression.
Keywords: rotational hardening    stress path    unified expression    constitutive model    

原状土体是在一定应力条件下形成的,土体颗粒结构产生定向排列,在宏观结构上形成一定的各向异性。同时,土体的变形与加载应力路径密切相关,采用不同的应力加载路径会得到不同的应变,引起初始各向异性的变化,形成土体颗粒之间的新的排列,这种新的各向异性称为诱发性各向异性。在屈服面模型中,这种各向异性通过屈服面的旋转来表示,这就是模型中采用的旋转硬化表达式。

研究人员对旋转硬化规律采用了不同的表达形式。Wheeler等[1]提出的旋转硬化表达式同时考虑了塑性体积应变和塑性剪切应变的影响,但旋转硬化参数不能直接确定,需要通过反分析,并且没有考虑洛德角的影响。Whittle等[2]提出的旋转硬化表达式采用临界状态线作为旋转的边界线,只考虑了塑性剪切应变对旋转硬化的影响,这与试验现象不吻合。Taiebat等[3]没有考虑塑性剪切应变的影响,同时,参数值不能直接确定,需要通过反分析, 屈服面的旋转存在边界线。Zhang等[4]依据交变加载试验结果提出了一个旋转硬化表达式,该表达式考虑了塑性剪切应变的作用,但是,参数需要通过反分析才能确定,给计算带来不便。Hashiguchi[5]认为屈服面旋转过程中塑性体积应变没有影响,采用K0线作为屈服面的倾斜线,这与已有的试验结果不符。Hueckel等[6]同时考虑了塑性体积应变和塑性剪切应变对旋转速度的影响,但是,旋转硬化参数的确定缺乏依据。Newson等[7]采用应力表示屈服面的旋转速度,屈服面的旋转只考虑了每一步的应力增量,没有考虑剪切应力的作用,同时,由于该本构模型采用非相关联流动法则,使得该旋转硬化表达式非常复杂。王立忠等[8]提出的表达式同时考虑了塑性体积应变和塑性剪切应变的作用,认为屈服面的旋转存在一定的范围,超过该范围时,对应屈服面的倾斜角度平衡值为0。对土体旋转硬化规律表达形式方面的研究成果还有很多[9-14],目前提出的旋转硬化表达式均是针对某一种具体的应力路径所提出的,当应力路径发生改变时,旋转硬化表达式的适用性有待验证。

1 旋转硬化方程的统一表示
1.1 旋转硬化方程的分类

研究人员针对旋转硬化规律采用了不同的表达形式,不同表达形式之间,塑性剪切应变和塑性体积应变对应的平衡值也不相同。依据加载应力比η和相对应屈服面的倾斜角度平衡值α之间的关系,这些旋转硬化表达式可划分为3种类别。

第1种:ηα之间为线性关系,不存在旋转极限线,方程式为

$ \alpha = x\eta $ (1)

式中:xηα之间的线性比例关系参数,其取值一般为1和0,这种关系式比较典型的旋转硬化规律包括Whittle等[2]、Taiebat等[3]及Zhang等[4]提出的旋转硬化表达式。

第2种:ηα之间为线性关系,但是存在旋转极限线,超出旋转范围极限线之后,屈服面对应的倾斜角度平衡值α为0,表达式为

$ \left\{ {\begin{array}{*{20}{l}} {\alpha = \eta }&{\left| \eta \right| \le \left| {M - \left| \alpha \right|} \right|}\\ {\alpha = 0}&{\left| \eta \right| > \left| {M - \left| \alpha \right|} \right|} \end{array}} \right. $ (2)

这种关系式比较典型的旋转硬化规律包括王立忠等[8]提出的旋转硬化表达式。

第3种:ηα之间为曲线关系,比较典型的包括Wheeler等[1]提出的旋转硬化规律表达式。

$ 3\left( {3\eta - 4\alpha } \right)\left( {{M^2} - {\eta ^2}} \right) = \pm 8\beta \left( {3\alpha - \eta } \right)\left( {\eta - \alpha } \right) $ (3)

为了分析这3种表达方式的合理性,收集了14种数据的应力比η/M和与应力比相对应的平衡值α/M之间的试验结果([1], [15-17])M为考虑了洛德角影响的临界状态线的斜率,试验结果与采用这3种表达式进行模拟的结果对比如图 1所示。

图 1 归一化的应力比η/M与平衡值α/M Fig. 1 Normalization of the stress ratio η/M and equilibrium valueα/M

图 1可以看出,在归一化应力比η/M较低的情况下,除了第1种表达式中x取0的情况外,依据这3种表达式得到的模拟结果与试验结果之间都吻合得比较好,当η/M的值达到一定值之后,塑性剪切应变的影响程度开始增加,随着归一化应力比η/M的增加,α/M反而开始下降。采用第2种关系式得到的曲线虽然也能反映出这种趋势,但是模拟曲线与试验结果之间存在明显差距。

1.2 旋转硬化方程的统一表示及参数确定
1.2.1 旋转硬化表达式

图 1可以看出,第3种表达式能较好地反映α/Mη/M的增加而下降的趋势,但试验结果与模拟结果之间存在一定误差,尤其在拉伸状态下,这种误差比较明显。同时,通过试验数据也可以看出,在应力比η较低的情况下,塑性体积应变占主导位置,ηα之间相等,当应力比η趋向于M时,α接近于0。基于此,提出旋转硬化规律的统一表达式为

$ {\rm{d}}\alpha {\rm{ = }}\mu \left[ {\left( {a\eta - \alpha } \right)\left[ {{\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}}} \right] - \beta \left( {b\eta - \alpha } \right)\left| {{\rm{d}}\varepsilon _{\rm{s}}^{\rm{p}}} \right|} \right] $ (4)

式中总共包含4个参数,分别是μβab,当这4个参数取不同的值时,就可以得到不同的旋转硬化规律表达式。同时,从图 1可以看出,在应力比η/M较低时,η/M近似等于α/M,当达到一定程度之后,α/M随着η/M的增大逐渐减小,最后趋向于0。为了反映这种趋势,本文提出参数ab分别取1和0,得到的旋转硬化规律表达式为

$ {\rm{d}}\alpha {\rm{ = }}\mu \left[ {\left( {\eta - \alpha } \right)\left[ {{\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}}} \right] + \beta \alpha \left| {{\rm{d}}\varepsilon _{\rm{s}}^{\rm{p}}} \right|} \right] $ (5)

为验证该旋转硬化方程的合理性,将该表达式嵌入修正剑桥模型,有关修正剑桥模型的具体介绍可参考文献[18]。建立了一个能描述各向异性的弹塑性本构模型,模型的屈服面方程表示为

$ f = \left( {q - \alpha p} \right) - \left( {{M^2} - {\alpha ^2}} \right)\left( {{p_{\rm{c}}} - p} \right)p=0 $ (6)
1.2.2 旋转硬化表达式参数的确定

旋转硬化表达式包含两个参数,分别是μβ,这里介绍对这两个参数的推导过程。若继续一维压缩过程,屈服面不发生旋转,即

$ {\rm{d}}\alpha = \mu \left[ {\left( {\eta - \alpha } \right)\left[ {{\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}}} \right] - \beta \alpha \left| {{\rm{d}}\varepsilon _{\rm{s}}^{\rm{p}}} \right|} \right] = 0 $ (7)

一维压缩过程中

$ \frac{{{\rm{d}}\varepsilon _{\rm{s}}^{\rm{p}}}}{{{\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}}}} = \frac{{\frac{2}{3}\left( {{\rm{d}}{\varepsilon _1} - {\rm{d}}{\varepsilon _3}} \right)}}{{\left( {{\rm{d}}{\varepsilon _1} + 2{\rm{d}}{\varepsilon _3}} \right)}} = \frac{2}{3} $ (8)

式中:dε1和dε3分别表示压缩过程中的单元体竖向应变和水平应变,联合式(7)和式(8),得到

$ \beta = \frac{{3\left( {{\eta _{{K_0}}} - {\alpha _{{K_0}}}} \right)}}{{2{\alpha _{{K_0}}}}} $ (9)

结合方程式(6),一维压缩中塑性体积应变和塑性剪切应变之间的关系表示为

$ \frac{{{\rm{d}}\varepsilon _{\rm{s}}^{\rm{p}}}}{{{\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}}}} = \frac{{\frac{{\partial f}}{{\partial q}}}}{{\frac{{\partial f}}{{\partial p}}}} = \frac{{2\left( {{\eta _{{K_0}}} - {\alpha _{{K_0}}}} \right)}}{{{M^2} - {\eta _{K_0^2}}}} = \frac{2}{3} $ (10)

得到

$ {\alpha _{{K_0}}} = \frac{{{\eta _{K_0^2}} + 3{\eta _{{K_0}}} - {M^2}}}{3} $ (11)

将式(11)代入式(9),可得到

$ \beta = \frac{{3\left( {{\eta _{{K_0}}} - {\alpha _{{K_0}}}} \right)}}{{2{\alpha _{{K_0}}}}} = \frac{{3\left( {{M^2} - \eta _{{K_0}}^2} \right)}}{{2\left( {\eta _{{K_0}}^2 + 3{\eta _{{K_0}}} - {M^2}} \right)}} $ (12)

如果将土体式样各向同性固结到先期压力的2~3倍左右,则土体的各向异性将基本消失[19],此时, 发生的塑性体积应变可表示为

$ \varepsilon _{\rm{v}}^{\rm{p}} = \frac{{\lambda - k}}{{1 + {e_0}}}\ln \left( {\frac{p}{{{p_0}}}} \right) = \frac{{\lambda - k}}{{1 + {e_0}}}\ln \left( {2 \sim 3} \right) \approx \frac{{\lambda - k}}{{1 + {e_0}}} $ (13)

各向同性固结时,应力比η=0, 塑性剪切应变和塑性体积应变之间的比例关系可表示为

$ \frac{{{\rm{d}}\varepsilon _{\rm{s}}^{\rm{p}}}}{{{\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}}}} = \frac{{\frac{{\partial f}}{{\partial q}}}}{{\frac{{\partial f}}{{\partial p}}}} = \frac{{2\left( {\eta - \alpha } \right)}}{{{M^2} - {\eta ^2}}} = \frac{{ - 2\alpha }}{{{M^2}}} $ (14)

将式(14)代入式(5), 得

$ \frac{1}{{\alpha \left( {2\alpha \beta - {M^2}} \right)}}{\rm{d}}\alpha {\rm{ = }}\frac{\mu }{{{M^2}}}{\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}} $ (15)

对式(15)左边积分,得

$ \int_{{\alpha _0}}^\alpha {\frac{1}{{\alpha \left( {2\alpha \beta - {M^2}} \right)}}{\rm{d}}\alpha {\rm{ = }}\frac{{2\beta }}{{{M^2}}}\left. {\left( {\ln \frac{{2\alpha \beta - {M^2}}}{{2\alpha \beta }}} \right)} \right|} _{{\alpha _0}}^\alpha $ (16)

假定屈服面累计旋转到初始值α0的0.1倍时,土体显示各向同性性质,这样联合式(13)、式(15)和式(16)就可求出参数μ的表达式。

$ \mu = \frac{{2\beta \left( {1 + {e_0}} \right)}}{{\lambda - \kappa }}\ln \frac{{10{M^2} - 2{\alpha _0}\beta }}{{{M^2} - 2{\alpha _0}\beta }} $ (17)

联合式(5)、式(12)和式(17),可建立α/M和η/M之间的关系曲线,如图 2所示。

图 2 应力比η/M与平衡值α/M Fig. 2 Stress ratio η/M and equilibrium value α/M

图 2可以看出,采用本文提出的旋转硬化表达式得到的曲线与试验结果之间吻合得较好,能很好地反映屈服面的倾斜角度先增后减,最后趋向于0的过程。将提出的旋转硬化表达式导入修正剑桥模型,模拟屈服面的旋转以及不同应力路径下的土体试样的应力-应变行为,模型需要输入的模型与修正剑桥模型需要输入的参数相同,总共包含6个,分别为e0λκMυpm0, 参数的物理意义可参考文献[18]。

2 不同应力路径试验的模拟

为了分析旋转硬化对不同加载路径下土体应力-应变关系的影响,选择3种不同的应力路径进行分析,分别是温州软土三轴排水试验[20]、BBC(Boston Blue Clay)土的三轴不排水试验[21]、Otaniemi土的循环加载试验[22],这3种土需要输入模型的参数见表 1

表 1 输入模型的参数 Table 1 Input parameter of the model

2.1 三轴排水试验

沈恺伦[20]对温州软土进行了三轴排水试验,试验过程先依照原状土的应力状态进行固结,然后再分别进行排水状态下的三轴压缩试验。图 3图 4表示的是应力路径为BCD50和BED33时,三轴压缩试验结果和采用不同旋转硬化表达式得到的模拟结果之间的对比。

图 3 应力路径BCD50的试验结果和模拟结果 Fig. 3 Simulated and experiment results of BCD50

图 4 应力路径BED33的试验结果和模拟结果 Fig. 4 Simulated and experiment results of BED33

图 3是应力路径为BCD50时试验结果和采用不同旋转硬化表达式时得到的模拟结果之间的对比,从图中可以看出,总体上模拟结果与试验结果之间吻合得比较好,采用不同旋转硬化表达式得到的模拟曲线之间非常接近,在模拟p-εv时,随着体积应变的发展,模拟曲线之间开始出现差距。当体积应变达到3%时,模拟曲线之间开始出现一定的差距,与试验结果之间的差距也慢慢增大,在模拟q-εs时,采用不同旋转硬化表达式得到的模拟结果与试验结果之间非常接近,都能较准确地模拟q-εs的发展趋势。

图 4是三轴拉伸试验BED-33时试验结果与模拟结果之间的对比,从图中可以看出,在模拟p-εv曲线时,采用不同旋转硬化表达式对模拟结果的影响不是很明显,采用4种旋转硬化表达式得到的模拟曲线之间非常接近,都能较准确地模拟p-εv之间的发展趋势,但在模拟q-εs曲线时,模拟结果与试验结果存在一定的误差,当剪切应变较小时,模拟曲线之间比较接近,随着剪切应变的发展,模拟曲线之间开始出现一定程度的误差,采用本文提出的旋转硬化表达式得到的模拟结果与试验结果之间吻合较好。

2.2 三轴不排水试验

Ladd等[21]对BBC(Bonston Blue Clay)土进行了三轴不排水压缩试验,所取土试样首先各向同性固结到150 kPa, 然后根据不同的超固结比值Roc,分别卸载到相应的应力值,这里选取的Roc分别为1和4,土样在整个加载过程中保持不排水条件,需要输入模型的土体参数见表 1。Yu等[23]也对该试验进行了数值研究。

图 5为模拟结果与试验结果之间的对比,从图中可以看出,在对应力路径p-q进行模拟时,采用不同旋转硬化表达式得到的模拟曲线与试验曲线之间比较接近,采用这4种旋转硬化表达式都能较准确地反映出p-q之间的变化趋势。在模拟q-εs时,当OCR为1,即当试样为正常固结时,模拟曲线位于试验曲线的上部,在出现相同的剪切应变εs时,模拟得到的剪切应力q要大于试验过程中的剪切应力,同时,采用不同旋转硬化表达式得到的模拟曲线之间也不是非常接近,相比较而言,采用本文提出的旋转硬化表达式得到的模拟曲线与试验曲线之间最为接近。当OCR为4时,在模拟应力路径p-q时,不同旋转硬化表达式得到的模拟结果与试验结果之间非常接近,都能准确地反映出应力路径p-q的变化趋势,在模拟q-εs时,模拟结果与试验结果之间存在一定的误差,但这种误差要小于OCR为1时试验结果与模拟结果之间的误差,总体而言,这4个旋转硬化表达式都能较准确地反映出变形趋势,采用本文提出的旋转硬化表达式得到的模拟结果与试验结果之间最为接近。

图 5 不同超固结比的试验结果与模拟结果 Fig. 5 Simulated and experiment results of samples with different overconsolidation rato

2.3 循环加载试验

循环加载试验选取Otaniemi重塑土进行了研究,这种土的相关参数及试验过程见Karstunen等[22]的研究成果,在三轴压缩之前,先将试样按照一定的应力路径η0进行固结,随后将试样进行卸载到一定值, 然后对试样进行加载卸载循环试验。本文选取了两组试样,分别是试样CAE3516R和试样CAE3519R,试样CAE3516R模拟结果如图 6所示。

图 6 CAE3516R试验结果与分析结果对比 Fig. 6 Simulated and experiment results of CAE3516R

图 6可以看出,对于CAE3516R试样,在模拟应力-应变曲线p-εv时,采用不同旋转硬化表达式对模拟结果有较大影响。初始加载时,采用Whittle等和王立忠等提出的旋转硬化表达式得到的模拟曲线位于试验曲线的上部,采用本文提出的旋转硬化表达式与Hueckel等提出的旋转硬化表达式得到的模拟曲线之间比较接近。卸载阶段,不同旋转硬化表达式得到的模拟结果之间的差距比加载阶段不同旋转硬化表达式得到的模拟结果之间的差距要明显一些,相比较而言,采用本文提出的旋转硬化表达式得到的模拟结果与试验结果之间更为接近。在模拟应力-应变关系q-εs时, 采用4种旋转硬化表达式都能反映出应力-应变关系的变化趋势,但是,相比较应力-应变关系p-εv而言,不同旋转硬化表达式得到的模拟曲线之间的差距更明显。

图 7为CAE3519R的试验结果与模拟结果之间的对比,从图中可以看出,在模拟应力-应变曲线p-εv时,加载阶段,不同旋转硬化表达式得到的模拟结果之间存在一定差距,卸载阶段,不同模拟结果之间则比较吻合。在模拟应力-应变关系q-εs时,不同旋转硬化表达式得到的模拟曲线之间差距比较明显,采用本文提出的旋转硬化表达式得到的模拟曲线与实测数据之间最为接近。

图 7 CAE3519R试验结果与分析结果对比 Fig. 7 Simulated and experiment results of CAE3519

3 结论

对提出的旋转硬化规律表达式进行总结,依据加载应力比η和相对应屈服面的倾斜角度平衡值α之间的关系,将提出的旋转硬化表达式划分为3种类型,分别是线性关系同时不存在旋转极限线、线性关系同时存在旋转极限线以及曲线关系。建立了这3种关系式的函数曲线,函数曲线与试验结果之间的对比结果显示,采用曲线表达形式与试验结果之间最为吻合。

在采用曲线型表达方程的基础上,提出了一个适用于不同应力路径下的旋转硬化表达式,同时考虑了塑性体积应变和塑性剪切应变对旋转速度的影响。表达式共包含两个参数,通过理论分析,得到这两个参数的确定方法。将所提出的旋转硬化表达式导入修正剑桥模型,建立一个能描述各向异性的本构模型,采用此模型模拟了不同应力路径下的应力-应变关系,模拟过程中,同时采用了其他3种旋转硬化表达式。将试验结果与模拟结果之间以及采用不同旋转硬化表达式得到的模拟结果之间进行对比分析,分析结果验证了提出的旋转硬化表达式的合理性。

参考文献
[1] WHEELER S J, NAATANEN A, KARSTUNEN M, et al. An anisotropic elastoplastic model for soft clay[J]. Canadian Geotechnical Journal, 2003, 40(2): 403–418. DOI:10.1139/t02-119
[2] WHITTLE A J, KAVVADAS M J. Formulation of MIT-E3 constitutive model for overconsolidated clays[J]. Journal of Geotechnical Engineering, 1994, 120(1): 173–198. DOI:10.1061/(ASCE)0733-9410(1994)120:1(173)
[3] TAIEBAT M, DAFALIAS Y F. Simple anisotropic sand plasticity model[J]. International Journal of Numerical and Analytical Methods in Geomechanics, 2008, 32(8): 915–948. DOI:10.1002/(ISSN)1096-9853
[4] ZHANG F, YE B, NODA T, et al. Explanation of cyclic mobility of soils approach by stress-induced anisotropy[J]. Soils and Foundations, 2007, 47(4): 635–648. DOI:10.3208/sandf.47.635
[5] HASHIGUCHI K. Elastoplastic constitutive equation of soils with the subloading surface and the rotational hardening[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1998, 12(8): 197–227.
[6] HUECKEL T, TUTUMLUER E. Modeling of elastic anisotropy due to one-dimensional plastic consolidation of clays[J]. Computers and Geotechnics, 1994, 16(4): 311–349. DOI:10.1016/0266-352X(94)90013-2
[7] NEWSON T A, DAVIES C R. A rotational hardening constitutive model for anisotropically consolidated clay[J]. Soils and Foundations, 1996, 36(3): 13–20. DOI:10.3208/sandf.36.3_13
[8] 王立忠, 沈恺伦. K0固结结构性软黏土的旋转硬化规律研究[J]. 岩土工程学报, 2008, 30(6): 863–872.
WANG L Z, SHEN K L. Rotational hardening law of K0 consolidated structured soft clays[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(6): 863–872. DOI:10.3321/j.issn:1000-4548.2008.06.014 (in Chinese)
[9] YIN Z Y, YIN J H, HUANG H W. Rate-dependent and long-term yield stress and strength of soft Wenzhou marine clay:experiments and modeling[J]. Marine Georesources and Geotechnology, 2015, 33(1): 79–91. DOI:10.1080/1064119X.2013.797060
[10] YIN Z Y, WANG J H. A one-dimensional strain-rate based model for soft structured clays[J]. Science in China Series, 2012, 55(1): 90–100.
[11] 姚仰平, 张丙印, 朱俊高. 土的基本特性、本构关系及数值模拟研究综述[J]. 土木工程学报, 2012, 45(3): 128–150.
YAO Y P, ZHANG B Y, ZHU J G. Behaviors, constitutive models and numerical simulation of soils[J]. China Civil Engineeing Journal, 2012, 45(3): 128–150. (in Chinese)
[12] 杨召焕, 王建华. 循环荷载作用下饱和软土的各向异性边界面模型[J]. 岩土力学, 2016, 31(Sup1): 63–70.
YANG Z H, WANG J H. Anisotropic bounding surface model for saturated soft clay under cyclic loading[J]. Rock and Soil Mechanics, 2016, 31(Sup1): 63–70. (in Chinese)
[13] 褚福永, 朱俊高, 赵颜辉, 等. 粗粒土初始各向异性弹塑性模型[J]. 中南大学学报, 2012, 43(5): 1914–1919.
ZHU F Y, ZHU J G, ZHAO Y H, et al. An elastoplastic model for granular soil considering initial anisotropy[J]. Journal of Central South University, 2012, 43(5): 1914–1919. (in Chinese)
[14] 刘正义.软黏土屈服特性及各向异性屈服面方程研究[D].杭州: 浙江大学, 2014.
LIU Z Y. Research on the yield behavior and anisotropic yield surface equation of soft clay[D]. Hangzhou: Zhejiang University, 2014. (in Chinese)
[15] KEVIN M. The stress-strain behavior of bothkennar clay[D]. Glasgow: University of Glasgow, 2006.
[16] 但汉波, 王立忠. 基于弹黏塑性本构模型的旋转硬化规律[J]. 岩石力学与工程学报, 2010, 29(1): 184–192.
DAN H B, WANG L Z. Rotational hardening law based on elastoviscoplastic constitutive model[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(1): 184–192. (in Chinese)
[17] 但汉波.天然软粘土的流变特性[D].杭州: 浙江大学, 2009.
DAN H B. Time dependent behavior of natural soft clays[D]. Hangzhou: Zhejiang University, 2009. (in Chinese)
[18] ROSCOE K H, BURLAND J B. On generalized stress strain behavior of wet clay in engineering plasticity[M]. London: Cambridge University Press, 1968.
[19] ANANDARAJAH A, KUGANENTHIRA N, ZHAO D. Variation of fabric anisotropy of kaolinite in triaxial loading[J]. Journal of Geotechnical Engineering, 1996, 122(8): 633–640. DOI:10.1061/(ASCE)0733-9410(1996)122:8(633)
[20] 沈恺伦.软粘土结构性、塑性各向异性及其演化[D].杭州: 浙江大学, 2006.
SHEN K L. Degradation of soil structure and anisotropic plasticity evolution of natural clay[D]. Hangzhou: Zhejiang University, 2006.(in Chinese)
[21] LADD C C, VARALLYAY J. The influence of the stress system on the behavior of saturated clays during undrained shear[R]. Research Rep. No. R65-11, Department of Civil Engineering, MIT, Cambridge, MA, 1965.
[22] KARSTUNEN M, KOSKINEN M. Plastic anisotropy of soft reconstituted clays[J]. Canadian Geotechnical Journal, 2008, 45(3): 314–328. DOI:10.1139/T07-073
[23] YU C, XU Q, YIN Z Y. Softening response under undrained compression following anisotropic consolidation[J]. Journal of Central South University, 2013, 20(6): 1703–1712. DOI:10.1007/s11771-013-1663-z