Three-dimensional non-coaxial constitutive model for sand and its application in bearing capacity of foundation
-
摘要: 传统弹塑性模型隐含了塑性应变率与应力方向共轴的假设,无法模拟砂土在非比例加载条件下的非共轴特性。对传统模型隐含的共轴假设进行了剖析,将Hashiguchi等的广义非共轴塑性流动理论引入Li等建立的砂土状态相关剪胀模型中,构建了一个砂土的非共轴模型。Leighton-Buzzard砂的单剪试验模拟表明,非共轴模型能够更加合理地描试验结果。开发了相应的有限元程序,利用所建模型对不同孔隙比条件下的砂土地基承载力问题进行了模拟。探讨了非共轴参数取值对计算结果的影响,总体上预测结果展现出与单剪试验模拟相一致的规律。Abstract: Due to the coaxial assumption that the direction of plastic strain rate is the same as that of stress, the non-coaxial behavior of sand under non-proportional loading conditions can not be captured by the conventional elastic-plastic model. In this study, firstly the cause of coaxial hypothesis contained in the traditional model is analyzed. Then, a non-coaxial model for sand is constructed by introducing the generalized non-coaxial plastic flow theory developed by Hashiguchi et al into the state-dependent dilatancy model for sand established by Li et al. The simulation of simple shear tests on Leighton-Buzzard sand shows that the non-coaxial model can describe the experimental results more reasonably. The corresponding finite element program is developed, the bearing capacity of sand foundation under different void ratios is analyzed by using the proposed model. The influence of non-coaxial parameters on the calculation results is discussed, and the predicted results demonstrate similar characteristics with those of simple shear test simulation on the whole.
-
Keywords:
- sand /
- non-coaxiality /
- simple shear test /
- bearing capacity of foundation
-
0. 引言
岩土材料的本构模型在有限元仿真当中扮演着核心角色,其合理性与否直接决定着有限元计算结果的精度。在传统岩土弹塑性模型中,塑性势函数通常表示为主应力或应力不变量的函数并结合正交流动法则确定塑性应变率的方向,这一过程隐含了塑性应变率与应力的共轴性假设。该假设与含主应力轴旋转加载条件下岩土材料的变形规律并不吻合[1]。
主应力轴方向的改变会显著影响土体的变形特性,这一点早已被土力学的学者所认识到[2-4]。众多研究者进行了大量试验,研究含主应力旋转条件下土体的变形规律,例如Miura等[5]、Ishihara等[6]的空心圆柱扭剪试验表明,主应力轴旋转条件下产生的塑性应变增量方向与应力方向明显不同。Gutierrez等[7-8]对Toyoura砂进行的应力探测试验表明:即便在应力幅值保持不变的条件下,单纯的主应力方向改变也会引起砂土产生塑性变形,并且应力的方向通常滞后于塑性应变率的方向,但二者之间的非共轴程度随着应力比的增加而逐渐变弱。Li等[9]开展的细观试验也证明了主应力旋转条件下颗粒材料的变形存在显著的非共轴性。传统的共轴弹塑性模型并不能描述这一现象。针对这一不足,国内外学者提出了多种新模型来模拟应力与塑性应变率之间的非共轴特点。其中一大类是建立在由Rudnicki等[10]提出的屈服角点理论基础上。该理论认为非共轴塑性应变率是由屈服函数切线方向的应力率所引发的。基于该理论,Yang等[11-13]分别结合Drucker-Prager模型和临界状态土力学模型CASM模拟了单剪试验中砂土的非共轴变形行为;罗强等[14]在应变软化的Drucker-Prager模型中引入带角点结构的塑性流动理论,对密砂的单剪试验进行了模拟分析。但由于该理论早期定义的切线应力率是基于二维应力空间,并不适用于三维化的模型,国内学者钱建固等[15]应用Gram-Schimidt张量正交化理论,发展了三维的非共轴塑性流动理论;李学丰等[16]利用该方法建立了一个砂土的非共轴模型,但没有将其应用于工程分析;基于钱建固等的研究成果,陈洲泉等[17-18]提出了一个能够计算定向剪切条件下非共轴变形的砂土模型,然而Dafalias等[19]指出在该加载模式下是否应当考虑非共轴性仍是一个值得商榷的问题。Hashiguchi等[20-22]发展了广义的非共轴塑性流动理论,建立了适用于黏土且能够考虑非共轴塑性变形的次加载面模型。此外,一些学者包括Gao等[23]、Tian等[24]从各向异性角度对砂土非共轴性进行了研究。
虽然目前为止已经建立了一些砂土的非共轴模型,但大多数文献中并未对传统弹塑性模型所隐含的共轴假设给出详细的论证,此外将建立的非共轴模型应用于工程分析的情况也并不多见。因此,笔者首先对传统弹塑性模型隐含的共轴性假设进行剖析,然后将Hashiguchi等[20-22]发展的广义非共轴塑性理论引入Li等[25]提出的砂土状态相关剪胀模型,建立砂土的非共轴本构模型;通过已有的砂土单剪试验数据检验模型的合理性,并利用所建模型对砂土地基的承载力问题进行分析,探究考虑非共轴塑性变形对计算结果的影响。
1. 共轴性假设剖析
1928年Mises通过类比弹性势的概念提出了塑性势的概念[1],随后建立的德鲁克公设表明,屈服函数与塑性势函数等价;塑性应变率的方向必定沿着屈服面的法线方向,即正交流动法则,其数学表达式为
˙εp=L∂f∂σ, (1) 式中,
˙εp 为塑性应变率张量,σ 为应力张量,f 为屈服函数,L 为加载指标,由一致性条件确定。传统的岩土弹塑性模型继承了这一正交流动法则,并且通常视土体为各向同性材料,将屈服函数与塑性势函数表示为应力不变量的函数。这种表示方式的便捷之处是只需在某一固定坐标系下建立土体的屈服规律和塑性流动规律,当坐标系发生改变时,建立的准则仍然适用,但却隐性地规定了塑性应变增量方向必须与应力方向共轴。为了对此进行更加清楚的说明,假设塑性势函数的一般形式为
g=g(I1,J2,J3,H)。 (2) 式中 3个应力不变量分别为
I1=tr(σ) ,J2=tr(s2)/2 ,J3=tr(s3)/3 ;H 为硬化变量;s 为偏应力张量。由正交流动法则,塑性应变率张量˙εp 可表示为˙εp=L∂g∂σ=L(∂g∂I1∂I1∂σ+∂g∂J2∂J2∂σ+∂g∂J3∂J3∂σ) =L[(∂g∂I1−∂g∂J323J2)1+∂g∂J2s+∂g∂J3s2]=L(a11+a2σ+a3σ2), (3) 式中,
1 为二阶单位张量,3个系数分别为a1=∂g∂I1− ∂g∂J2I13+∂g∂J3(I219−23J2) ,a2=∂g∂J2−∂g∂J323I1 ,a3=∂g∂J3 。设
n1 ,n2 ,n3 分别为应力张量3个主方向上的单位基矢量,则必然存在某一单位正交张量Q 可以将应力张量σ 转化为主应力张量的形式,即Q⋅σ⋅QT=σ1n1⊗n1+σ2n2⊗n2+σ3n3⊗n3, (4) 式中,
σ1 ,σ2 ,σ3 分别为应力张量的3个主值。将塑性应变率张量˙εp 两侧点积张量Q ,并利用Q 的正交性质QT⋅Q=1 ,则可得Q⋅˙εp⋅QT=LQ⋅(a11+a2σ+a3σ2)⋅QT =L(a1Q⋅1⋅QT+a2Q⋅σ⋅QT+ a3Q⋅σ⋅QT⋅Q⋅σ⋅QT)=˙εp1n1⊗n1+˙εp2n2⊗n2+˙εp3n3⊗n3, (5) 式中,
˙εpi=L(a1+a2σi+a3σ2i) ,(i=1,2,3) 为塑性应变率张量的3个主值。式(5)表明,应力张量σ 与塑性应变率张量˙εp 在同一组基矢量中化为主张量的形式,即二者有相同主方向,因此按照传统位势理论计算的塑性应变率张量˙εp 总是与应力张量σ 共轴。2. 砂土的三维非共轴本构模型
由于在非比例加载条件下会涉及到三维应力状态,需要将模型三维化[19,26-27]。本文选择在偏平面上插入角隅函数的方法实现这一目的。
2.1 弹性部分
模型的剪切模量
G 、体积模量K 都与砂土密度和平均应力p 有关,其表达式为G=G0(2.97−e)21+e√ppa, (6) K=G2(1+ν)3(1−2ν), (7) 式中,
G0 为材料常数,e 为当前孔隙比,pa 为大气压力,ν 为泊松比。2.2 状态变量
模型的状态变量采用Been等[28]提出的形式,其定义为相同压力下当前孔隙比与临界孔隙比的差值,即
ψ=e−ec, (8) 式中,
ec 为临界孔隙比。对于砂土,Li等[25]提出了一个适用范围较广的形式,其数学形式为ec=eΓ−λc(p/pa)ξ, (9) 式中,
eΓ ,λc 和ξ 为3个材料常数。2.3 屈服函数与塑性模量
屈服函数采用广义Mohr-Coulomb形式,即
f=q−Mg(θ)p=0。 (10) 式中
p=I1/3 ;q=√3J2 ;M 为硬化内变量;应力罗德角θ=sin−1(3√3J3/2J32)/3 ;g(θ) 为引入的角隅函数,这里采用Li等[29]建议的形式:g(θ)=√(1+c2)2+4(1−c2)sin3θ−(1+c2)2(1−c)sin3θ, (11) 式中,
c=Me/Mc ,Me ,Mc 分别为三轴拉伸与压缩条件下的应力比比值。模型中定义了一个状态相关的塑性模量
Kp ,可以表示为Kp=hGenψ[Mcsg(θ)e−nψη−1]。 (12) 式中
Mcs 为临界状态应力比;应力比η=q/q ;n 为材料常数;Li等[29]建议h 取为孔隙比的函数h=h1−h2e ,h1 ,h2 为2个材料常数。2.4 剪胀方程
剪胀方程表示为塑性体应变率
˙εpv 与广义塑性剪应变率˙εpq 之比,具体形式为d=˙εpv˙εpq=˙εpvL=d0Mcsg(θ)[Mcsg(θ)emψ−η], (13) 式中,
d0 ,m 为2个材料常数。上述剪胀方程的等效势函数为g=q−d0Mcsg(θ)pd0−Mcsg(θ)[1−(pp0)d0Mcsg(θ)−1], (14) 式中,
p0 为塑性势面与p 轴的交点。2.5 硬化规律
根据式(10),由一致性条件可得
˙f=∂f∂σ:˙σ+∂f∂Μ∂Μ∂εpq˙εpq=0, (15) 式中,
˙σ 为应力率张量。根据塑性理论[30]˙εpq=L=1Kp∂f∂σ:˙σ。 (16) 结合式(12),(15),(16),可得
∂Μ∂εpq=hGenψg(θ)p[Mcsg(θ)e−nψη−1]。 (17) 2.6 非共轴塑性应变率
与传统塑性理论不同,Rudnicki等[10]提出的屈服角点理论将总塑性应变率分解为共轴塑性应变率和沿屈服函数切线方向的塑性应变率之和,则总应变率张量
˙ε 为˙ε=˙εe+˙εpc+˙εpt。 (18) 式中
˙εe 为弹性应变率张量;˙εpc 为共轴塑性应变率张量,由传统塑性位势理论确定;˙εpt 为非共轴塑性应变率张量,Rudnicki等[10]认为˙εpt 与静水压力部分无关,其方向沿着屈服函数切线方向并与s 正交,˙εpt 的表达式为˙εpt=1Ht(˙s−s:˙ss:ss), (19) 式中,
Ht 为非共轴塑性模量,˙s 为偏应力率张量。然而式(19)由于没有考虑第三应力不变量的影响,并不适用于一般化的三维模型。Hashiguchi等[20-22]发展了适用于三维应力空间的广义非共轴塑性流动理论,并假设非共轴塑性应变率张量与偏切线应力率张量之间为线性关系,
˙εpt 可表示为˙εpt=TEe−1:˙st, (20) 式中,
T 为一个无量纲的标量系数,Ee 为弹性刚度张量,˙st 为偏切线应力率张量,如图1所示。Gutierrez 等[7-8]的研究表明,当应力比接近临界应力比时,非共轴现象逐渐减弱,因此通常情况下非共轴系数并不是一个常量。Hashiguchi等[20-22]将非共轴系数表达为相似比的指数函数;Yang等[11]建议将T 表示为广义塑性剪应变的指数函数,并指出这种表示形式可显著提高非共轴模型的收敛性。本文中,T 取如下形式:T=T0exp(−χεpq), (21) 式中,
T0 为T 的初始值,χ 为材料常数。如图1所示,偏切线应力率张量
˙st 可以由下式计算:˙st=¯Ι′t:˙σ=˙s−˙sn, (22a) ¯Ι′t=ˉI′−n⊗n, (22b) n=(ˉn−13tr(ˉn)1)/‖ˉn−13tr(ˉn)1‖, (22c) ˉn=(∂f/∂σ)/‖(∂f/∂σ)‖。 (22d) 式中
ˉn 为屈服面的外法线方向;n 为偏法线方向;˙sn 为偏法线应力率张量;ˉI′ 为偏投影张量,其分量形式为ˉI′ijkl=(δikδjl+δilδjk)/2−δijδkl/3 ,δij 为Kronecker符号;¯Ι′t 为偏切线投影张量。结合式(16),(18),(20),(22a),总应变率张量可以表示为˙ε==(Ee−1+1Kp∂g∂σ⊗∂f∂σ+TEe−1:ˉΙ′t):˙σ。 (23) 对式(23)求逆可得非共轴弹塑性本构关系:
˙σ=Encep:˙ε=[Eep−T1+TˉΙ′t:Eep]:˙ε, (24a) Eep=Ee−Ee:∂g∂σ⊗∂f∂σ:EeKp+∂f∂σ:Ee:∂g∂σ, (24b) 式中,
Eep 为共轴弹塑性刚度张量,Encep 为非共轴弹塑性刚度张量。3. 模型的验证
为展示考虑非共轴塑性应变的影响并验证模型的合理性,采用上述建立的砂土非共轴模型对Leighton-Buzzard砂的单剪试验进行模拟。单剪试验结果来自文献[31],与Li等[25]模型相关的参数选择蔡正银等在文献[32]中给出的参数值,
χ 的值参考Yang等[11]的研究取为20,由于单剪试验中,一般很难获得塑性应变增量的大小,因此T0 的值通常需根据试验曲线拟合确定,具体参数见表1。蔡正银等[32]用于确定试验参数的Leighton-Buzzard砂为粉细砂,其粒径范围为0.09~0.15 mm,颗粒相对密度为2.65,最大和最小干密度分别为1.59,1.32 g/cm3,对应的最大和最小孔隙比分别为emax=1.0 和emin=0.67 ,而文献[31]的单剪试验中,Leighton-Buzzard砂的初始孔隙比e0 分为0.53,0.64,0.75共3种,该砂的最大、最小孔隙比显然与文献[32]中的砂并不相同,因此对部分模型参数进行了调整,将eΓ=0.998 ,n=1.45 ,Mcs=1.185 分别调整为eΓ=0.908 ,n=1.15 ,Mcs= 1.105 。调整该3个参数的主要原因是文献[31]中e0 为0.53,0.64,0.75时,分别对应于密砂、中密砂、松砂,而按照文献[32]中提供的原始模型参数,这3种孔隙比都对应于密砂情况。通过调整eΓ 可使砂在这3种初始孔隙比情况下的应力应变曲线与典型的密砂、中密砂、松砂的应力应变曲线相对应;由于峰值应力比对应于Kp=0 ,而Kp 对于参数n 是较为敏感的,因此调节参数n ,Mcs 以使预测的峰值应力比和临界应力比与试验数据基本一致。弹性参数 临界状态参数 剪胀参数 非共轴参数 G0=150 Mcs=1.105 d0=0.52 T0=10 ν=0.35 eΓ=0.908 m=2.42 χ=20 λ=0.0198 n=1.15 ξ=0.68 h=0.62 试验共分为3组,对应Leighton-Buzzard砂的初始孔隙比
e0 分别为0.53,0.64,0.75,静止侧压力系数K0=0.5 ,竖向应力σy=100 kPa 。试验过程中,水平应变εx=0 ,竖向应力σy 保持不变。加载过程中研究主应力方向偏转角α ,主塑性应变率方向偏转角β ,如图2所示,其中α 与β 的定义分别为tan2α=2τxyσy−σx, (25a) tan2β=2˙εpxy˙εpy−˙εpx, (25b) 式中,
σx 为水平向应力,τxy 为剪应力,˙εpx ,˙εpy 分别为水平与竖直向塑性应变率,˙εpxy 为塑性剪应变率。图3为模拟结果与试验结果的对比。由图3(a),(c),(e)可以看出,在加载前期,非共轴项的影响较为显著,非共轴模型预测的归一化应力比小于共轴模型的预测值,随着剪应变的增加,非共轴项的影响逐渐减弱,当剪应变
γxy>0.05 后,非共轴模型与共轴模型预测的应力–应变曲线基本重合。由式(24)可知,由于非共轴模型比共轴模型多了非共轴塑性应变项,在相同的应力增量条件下,非共轴模型将计算出更大的应变增量,因此预测的应力–应变曲线初始斜率减小,表现出刚度降低,更接近试验结果。图3(a),(c),(e)中,共轴模型预测的峰值应力比分别为0.85,0.72,0.615,非共轴模型预测结果分别为0.84,0.71,0.61。考虑非共轴影响,峰值应力比会有所下降,减弱模型的软化效应,但是程度很小,因此考虑非共轴塑性应变基本不影响峰值应力比,也不改变最终的临界应力比,但会延缓峰值应力的出现。对不同密度砂土而言,非共轴塑性应变的影响并不存在明显的差异。由图3(b),(d),(f)可以发现,主应力方向与主塑性应变率方向从开始加载就存在偏离。在加载早期阶段,应力的方向落后于塑性应变率的方向,随着剪应变的增加,应力比接近临界应力比时,主应力方向与塑性应变率方向之间的差别消失。非共轴模型能够反映主应力与主塑性应变率之间的非共轴情况,而共轴模型不能模拟这一现象。4. 砂土地基承载力分析
为了分析非共轴塑性变形对于地基承载力的影响,笔者基于MATLAB软件开发了相关的有限元计算程序。隐式的回退Euler算法不适用于角点类的非共轴模型,程序中单元高斯积分点上的应力更新采用了显式向前的Euler积分算法。由于非共轴项的引入以及程序采用非对称的刚度矩阵进行计算,导致有限元分析不易收敛,因此计算过程将总位移荷载分为2000个较小的加载步施加到地基上。此外,在每一加载步中的整体平衡迭代部分,程序也采用了适当放松的收敛准则:当最大节点不平衡力与等效节点荷载的比值小于1%时即认为已经平衡,程序进入下一加载步的计算。
地基模型如图4所示,地基宽为60 m,深为20 m,顶部条形基础宽为10 m,假设为完全粗糙的刚性基础。模型共划分为1200个单元,单元类型为四边形等参单元。加载共分为3步,首先生成初始地应力场,然后在模型顶部施加一个20 kPa的均布压力,最后对刚性基础施加2 m向下的位移荷载。加载过程中固定底部边界位移及两侧边界的水平位移。地基材料为标准Toyoura砂,其平均粒径为0.24 mm,颗粒相对密度为2.64,最大和最小孔隙比分别为
emax=0.97 和emin=0.63 。为简化计算,假设地基中的孔隙比均匀分布。将初始孔隙比分为e0=0.8 与e0=0.9 两种情况,对应的重度分别为14.6,13.8kN/m3 。与Li等[25]模型相关的参数取自文献[25],χ 取与单剪试验模拟相同值,T0 分别取0.0,2.5,5.0三个值以研究其大小对计算结果的影响。具体材料参数见表2。弹性参数 临界状态参数 剪胀参数 硬化参数 非共轴参数 G0=125 Mcs=1.25 d0=0.88 n=1.1 T0=(0.0, 2.5, 5.0) ν=0.05 eΓ=0.934 m=3.5 h1=3.15 χ=20 λ=0.019 h2=3.05 ξ=0.7 图5,6分别为孔隙比
e0=0.8 时计算得到的位移–压力曲线与塑性剪应变云图;图7,8分别为孔隙比e0=0.9 时计算的位移–压力曲线与塑性剪应变云图。由图5,7可知,在加载的前、中期阶段,非共轴模型预测的基础压力小于共轴模型的预测结果,并且T0 越大时,这种差异越明显,而在加载后期,差异逐渐减小,当竖向位移达到2 m时,差异基本消失。如图5所示,对于e0=0.8 的情况,共轴模型与非共轴模型预测的基础压力都出现了峰值,分别为3034 kPa(T0= 0.0 ),2990 kPa(T0=2.5 ),2970 kPa(T0=5.0 ),峰值压力基本一致,因此,非共轴塑性应变基本不影响峰值基础压力的大小,但会延缓峰值基础压力的出现,这与前述单剪试验模拟所展示的规律是相似的。如图5中的A′ ,B′ 两点所示,非共轴模型(T0=5.0 )与共轴模型预测的基础压力差在竖向位移为0.63 m时最大,相差为620 kPa,B′ 点压力值为A′ 点的74%;对于e0=0.9 的情况,如图7中的C′ ,D′ 两点所示,非共轴模型(T0=5.0 )与共轴模型预测值在竖向位移为0.626 m时相差最大,为305 kPa,D′ 点压力值为C′ 点的76%,这一比例与e0=0.8 的情况非常接近,此外,两种不同密度情况下,非共轴塑性应变的影响都存在先增大后减小的趋势,因此,对不同密度的砂土地基,非共轴项的影响并不存在明显差异,这与不同初始孔隙比的砂土单剪试验模拟结果也是一致的。由图6(a)可知,基础竖向位移达到2 m时,共轴模型预测显示地基内已出现连续的滑动面,滑动边界附近土体全部处于塑性状态,说明地基稳定性已经丧失;图6(b),(c)与图6(a)的对比表明考虑非共轴塑性应变会延缓极限状态的到来,但图6(b)中的塑性区略大于图6(c),暗示出图6(b)会更早出现连续滑动面,表明T0 越大,延缓趋势越明显。分析图8可知,随着T0 的增大,模型预测的塑性区域存在减小趋势,但差异并不如图6所示明显。整体来看,地基承载力问题的预测结果表现出单剪试验模拟结果相一致的规律,随着非共轴系数
T0 的增大,非共轴项对预测结果的影响也更加显著。5. 结论
(1)剖析了传统弹塑性理论的局限性,指出了传统位势理论所隐含共轴假设的原因。
(2)通过在状态相关剪胀模型中引入Hashiguchi等发展的广义非共轴塑性流动理论,建立了砂土的非共轴模型。对砂土单剪试验的模拟表明,建立的模型能够合理预测主应力轴旋转过程中的塑性应变率与应力之间的非共轴现象;引入非共轴塑性应变项会降低模型的刚度,延缓峰值应力的出现,但不改变峰值应力比与临界应力比的大小;对不同初始孔隙比的砂土,非共轴塑性应变的影响并不存在明显的差异。
(3)基于MATLAB软件开发了相应的有限元程序,通过对砂土地基承载力的分析表明,考虑非共轴塑性应变后,地基整体刚度减小,预测的位移-压力曲线更加平缓,峰值基础压力的大小基本不受影响,但非共轴项会延缓峰值基础压力的出现。此外,非共轴项的影响规律也与砂的密度关系不大。总体来看,预测结果所展示的规律与单剪试验模拟结果是相似的。
-
表 1 Leighton-Buzzard砂的模型参数
Table 1 Model parameters of Leighton-Buzzard sand
弹性参数 临界状态参数 剪胀参数 非共轴参数 G0=150 Mcs=1.105 d0=0.52 T0=10 ν=0.35 eΓ=0.908 m=2.42 χ=20 λ=0.0198 n=1.15 ξ=0.68 h=0.62 表 2 Toyoura砂的模型参数
Table 2 Model parameters of Toyoura sand
弹性参数 临界状态参数 剪胀参数 硬化参数 非共轴参数 G0=125 Mcs=1.25 d0=0.88 n=1.1 T0=(0.0, 2.5, 5.0) ν=0.05 eΓ=0.934 m=3.5 h1=3.15 χ=20 λ=0.019 h2=3.05 ξ=0.7 -
[1] 郑颖人, 孔亮. 岩土塑性力学[M]. 北京: 中国建筑工业出版社, 2010. ZHENG Ying-ren, KONG Liang. Geotechnical Plastic Mechanics[M]. Beijing: China Architecture and Building Press, 2010. (in Chinese)
[2] ARTHUR J R F, RODRIGUEZ DEL C J I, DUNSTAN T. Principal stress rotation: a missing parameter[J]. Journal of the Geotechnical Engineering Division, 1980, 106(4): 419-433. doi: 10.1061/AJGEB6.0000946
[3] ROSCOE K H. The influence of strains in soil mechanics[J]. Géotechnique, 1970, 20(2): 129-170. doi: 10.1680/geot.1970.20.2.129
[4] ODA M, KONISHI J. Microscopic deformation mechanism of granular material in simple shear[J]. Soils and Foundations, 1974, 14(4): 25-38. doi: 10.3208/sandf1972.14.4_25
[5] MIURA K, MIURA S, TOKI S. Deformation behaviour of sand under principal axes rotation[J]. Soils and Foundations, 1986, 26(1): 36-52. doi: 10.3208/sandf1972.26.36
[6] ISHIHARA K, TOWHATA I. Sand response to cyclic rotation of principal stress direction as induced by wave loads[J]. Soils and Foundations, 1983, 23: 11-16. doi: 10.3208/sandf1972.23.4_11
[7] GUTIERREZ M, ISHIHARA K, TOWHATA I. Flow theory for sand during rotation of principal stress direction[J]. Soils and Foundations, 1991, 31(4): 121-132. doi: 10.3208/sandf1972.31.4_121
[8] GUTIERREZ M, ISHIHARA K. Non-coaxiality and energy dissipation in granular materials[J]. Soils and Foundations, 2000, 40(2): 49-59. doi: 10.3208/sandf.40.2_49
[9] LI X, YANG D, YU H S. Macro deformation and micro structure of 3D granular assemblies subjected to rotation of principal stress axes[J]. Granular Matter, 2016, 18(3): 53. doi: 10.1007/s10035-016-0632-2
[10] RUDNICKI J W, RICE J R. Conditions for the localization of deformation in pressure-sensitive dilatant materials[J]. Journal of the Mechanics and Physics of Solids, 1975, 23(6): 371-394. doi: 10.1016/0022-5096(75)90001-0
[11] YANG Y, YU H S. Numerical aspects of non-coaxial model implementations[J]. Computers and Geotechnics, 2010, 37(1/2): 93-102.
[12] YANG Y, YU H S. Application of a non-coaxial soil model in shallow foundations[J]. Geomechanics and Geoengineering: An International Journal, 2006, 1(2): 139-150. doi: 10.1080/17486020600777101
[13] YANG Y, YU H S. A non-coaxial critical state soil model and its application to simple shear simulations[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2006, 30(13): 1369-1390. doi: 10.1002/nag.531
[14] 罗强, 郑伟花. 密砂单剪试验的非共轴本构数值分析[J]. 长江科学院院报, 2015, 32(5): 89-94. https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201505020.htm LUO Qiang, ZHENG Wei-hua. Numerical analysis of simple shear test on dense sand with non-coaxial constitutive model[J]. Journal of Yangtze River Scientific Research Institute, 2015, 32(5): 89-94. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201505020.htm
[15] 钱建固, 黄茂松, 杨峻. 真三维应力状态下土体应变局部化的非共轴[J]. 岩土工程学报, 2006, 28(4): 511-515. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200604015.htm QIAN Jian-gu, HUANG Mao-song, YANG Jun. Effect of non-coaxial plasticity on onset strain localization in soils under 3D stress condition[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(4): 511-515. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200604015.htm
[16] 李学丰, 黄茂松, 钱建固. 宏–细观结合的砂土单剪试验非共轴特性分析[J]. 岩土力学, 2013, 34(12): 3418-3424. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201312012.htm LI Xue-feng, HUANG Mao-song, QIAN Jian-gu. Analysis of non-coaxial characters of sand for simple shear test with the method of macro-meso-incorporation[J]. Rock and Soil Mechanics, 2013, 34(12): 3418-3424. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201312012.htm
[17] 陈洲泉, 黄茂松. 砂土各向异性与非共轴特性的本构模拟[J]. 岩土工程学报, 2018, 40(2): 244-251. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201802005.htm CHEN Zhou-quan, HUANG Mao-song. Constitutive modeling of anisotropic and non-coaxial behaviors of sand[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(2): 244-251. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201802005.htm
[18] 陈洲泉, 黄茂松. 基于状态相关本构模型的砂土非共轴特性模拟[J]. 岩土力学, 2017, 38(7): 1960-1966. doi: 10.16285/j.rsm.2017.07.015 CHEN Zhou-quan, HUANG Mao-song. Simulation of non-coaxial characteristics of sandy soil based on state-dependent constitutive model[J]. Rock and soil mechanics, 2017, 38(7): 1960-1966. (in Chinese) doi: 10.16285/j.rsm.2017.07.015
[19] DAFALIAS Y F, PAPADIMITRIOU A G, LI X S. Sand plasticity model accounting for inherent fabric anisotropy[J]. Journal of Engineering Mechanics, 2004, 130(11): 1319-1333. doi: 10.1061/(ASCE)0733-9399(2004)130:11(1319)
[20] HASHIGUCHI K. The tangential plasticity[J]. Metals and Materials, 1998, 4(4): 652-656. doi: 10.1007/BF03026374
[21] HASHIGUCHI K, TSUTSUMI S. Elastoplastic constitutive equation with tangential stress rate effect[J]. International Journal of Plasticity, 2001, 17(1): 117-145. doi: 10.1016/S0749-6419(00)00021-8
[22] HASHIGUCHI K, TSUTSUMI S. Shear band formation analysis in soils by the subloading surface model with tangential stress rate effect[J]. International Journal of Plasticity, 2003, 19(10): 1651-1677. doi: 10.1016/S0749-6419(02)00113-4
[23] GAO Z W, ZHAO J D. Strain localization and fabric evolution in sand[J]. International Journal of Solids and Structures, 2013, 50(22/23): 3634-3648.
[24] TIAN Y, YAO Y P. Constitutive modeling of principal stress rotation by considering inherent and induced anisotropy of soils[J]. Acta Geotechnica, 2018, 13(6): 1299-1311.
[25] LI X S, DAFALIAS Y F. Dilatancy for cohesionless soils[J]. Géotechnique, 2000, 50(4): 449-460.
[26] 姚仰平, 路德春, 周安楠. 岩土类材料的变换应力空间及其应用[J]. 岩土工程学报, 2005, 27(1): 24-29. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200501002.htm YAO Yang-ping, LU De-chun, ZHOU An-nan. Transformed stress space for geomaterials and its application[J]. Chinese Journal of Geotechnical Engineering, 2005, 27(1): 24-29. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200501002.htm
[27] 姚仰平. UH模型系列研究[J]. 岩土工程学报, 2015, 37(2): 193-217. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201502002.htm YAO Yang-ping. Advanced UH models for soils[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(2): 193-217. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201502002.htm
[28] BEEN K, JEFFERIES M G. A state parameter for sands[J]. Géotechnique, 1985, 35(2): 99-112.
[29] LI X S, DAFALIAS Y F. Constitutive modeling of inherently anisotropic sand behavior[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2002, 128(10): 868-880.
[30] DAFALIAS Y F. An anisotropic critical state soil plasticity model[J]. Mechanics Research Communications, 1986, 13(6): 341-347.
[31] COLE E R L. The Behavior of Soils in the Simple Shear Apparatus[D]. Cambridge: Cambridge University, 1967.
[32] 蔡正银, 李相菘. 砂土的剪胀理论及其本构模型的发展[J]. 岩土工程学报, 2007, 29(8): 1122-1128. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200708004.htm CAI Zheng-yin, LI Xiang-song. Development of dilatancy theory and constitutive model of sand[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(8): 1122-1128. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200708004.htm