Processing math: 100%
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

考虑砂土初始各向异性的单剪试验模拟分析

吴则祥, 陈佳莹, 尹振宇

吴则祥, 陈佳莹, 尹振宇. 考虑砂土初始各向异性的单剪试验模拟分析[J]. 岩土工程学报, 2021, 43(6): 1157-1165. DOI: 10.11779/CJGE202106020
引用本文: 吴则祥, 陈佳莹, 尹振宇. 考虑砂土初始各向异性的单剪试验模拟分析[J]. 岩土工程学报, 2021, 43(6): 1157-1165. DOI: 10.11779/CJGE202106020
WU Ze-xiang, CHEN Jia-ying, YIN Zhen-yu. Finite element simulation of simple shear tests considering inherent anisotropy[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(6): 1157-1165. DOI: 10.11779/CJGE202106020
Citation: WU Ze-xiang, CHEN Jia-ying, YIN Zhen-yu. Finite element simulation of simple shear tests considering inherent anisotropy[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(6): 1157-1165. DOI: 10.11779/CJGE202106020

考虑砂土初始各向异性的单剪试验模拟分析  English Version

基金项目: 

港研究资助局基金项目 R5037-18F

详细信息
    作者简介:

    吴则祥(1987—),男,讲师,博士,主要从事岩土试验及数值等方面的教学和科研工作。E-mail: zexiang.wu@wzu.edu.cn

    通讯作者:

    尹振宇, E-mail: zhenyu.yin@polyu.edu.hk

  • 中图分类号: TU43

Finite element simulation of simple shear tests considering inherent anisotropy

  • 摘要: 在单剪条件下,由于侧边界上不存在互补的剪切应力会导致应力不均匀,使得室内单剪试验研究变得复杂。提出一个考虑砂土初始各向异性的单剪试验数值分析方法。主要是通过引入初始各向异性修正的临界状态本构模型,并将模型编入有限元大变形计算平台。利用三维有限元分析实际的GDS型单剪试验条件,以分析在剪切过程中试样的不均匀特征。此项研究可以提高对单剪试验过程中边界效应的理解和认识,并且提供一种分析试样的不均匀性的计算方法。
    Abstract: The absence of complementary shear stress on the side boundary of a specimen in the simple shear tests will cause stress inhomogeneity. An enhanced critical state-based constitutive model is proposed by incorporating the inherent anisotropy fabric, and also implemented into the finite element code for the numerical simulation. In addition, a three-dimensional finite element analysis with the same size as the GDS-type simple shear apparatus is performed to illustrate the inhomogeneous features of the specimen. Above all, this study can improve the understanding and knowledge of boundary effects for the simple shear tests, and provide a calculation method for analyzing the inhomogeneous ness of the specimen.
  • 单剪试验被广泛应用于评价土体的剪切强度,和模拟现场加载条件下的力学响应(如坡度,打桩,滑坡,地震,顶管隧道)。然而,由于单剪过程中某些力学特性无法直接测量,比如确定试样剪切带位置、侧向边界的径向应力、主应力/应变旋转程度等[1-2],均是单剪试验研究所亟需解决的问题。此外,与理想纯剪试验相比,单剪试验侧向边界不存在互补的剪切应力,这将导致圆柱形试样的应力不均匀性[3],并使试验研究变的复杂。因此,建议提出一种有效的方法来评估单剪过程中试样的不均匀性。

    有限元法(FEM)是研究土样应力/应变分布的理想分析工具。已经有研究者利用二维模拟再现了单剪试验中的应力不均匀特性[4-5]。然而,由于二维平面应变假设只适用于立方体试样,因此只能提供有限的结果和关于边界效应的讨论。最近,Doherty等[6]使用修正剑桥模型提出了一种更先进的三维有限元建模方法。但他的是采用高岭土力学参数,建模中忽略了土体的各向异性。

    本文利用考虑初始各向异性的改进的SIMSAND模型和ABAQUS大变形计算模块,提出了一种模拟单剪试验的数值方法。首先将改进的SIMSAND模拟写入ABAQUS/explicit自定义用户子程序vumat来实现有限元的计算,然后基于GDS型单剪试验,建立相应等尺寸的三维计算模型模拟一系列枫丹白露砂试验,最后通过模拟结果研究试样不均匀性的边界效应。

    本文基于一种简单的临界状态砂土模型[7-11](SIMSAND模型),其基本本构方程如表1所示,综合考虑非线性弹性、非线性塑性硬化法则、修正剪胀剪缩特征、三维强度和临界状态等特征,来有效地描述砂土特征。

    表  1  SIMSAND模型的基本本构方程
    Table  1.  Basic constitutive equations of SIMSAND
    组成部分本构方程
    弹性准则ε˙ije=1+υ3K(12υ)σijυ3K(12υ)σkkδij K=K0pat(2.97e)2(1+e)(ppat)n 
    屈服面f=qpH 
    塑性势面gp=Ad(Mptqp), gq=1 
    硬化参数H=Mpεdpkp+εdp 
    临界状态参数ec=erefλ(ppat)ξ tanϕp=(ece)nptanϕμ   tanϕpt=(ece)ndtanϕμ 
    三维强度标准Mp=6sinϕp3sinϕp[2c141+c14+(1c14)sin3θ]14, c1=3sinϕp3+sinϕpMpt=6sinϕpt3sinϕpt[2c241+c24+(1c24)sin3θ]14 , c2=3sinϕpt3+sinϕpt 
    下载: 导出CSV 
    | 显示表格

    组构各向异性已被广泛证实会影响颗粒材料的强度[12-14]。在单剪条件下,引入组构各向异性联合不变量,可将原始各向同性强度扩展到砂的各向异性强度。在这种方法中,利用由Oda等[15]提出的组构张量Fij来描述组构各向异性的程度:

    Fij=[Fz000Fx000Fy]=13+Δ[1Δ0001+Δ0001+Δ]=13[100010001]+2Δ9+3Δ[200010001], (1)

    式中,Δ为描述颗粒材料组构各向异性分布的参数,0<Δ<1为典型砂土的典型范围。

    根据张量理论[16]采用联合不变量A,将组构张量的偏差部分与偏应力张量联系起来:

    A=sijdijsmnsmndmndmn, (2)

    式中,sij为偏应力张量(sij = σij-pδij),dij为组构张量的偏差部分(dij= Fij-Fkkδij/3),δij为Kronecker delta。

    结合各向异性强度的方式类似于Li等[17]、Gao等[14]的方法,是将各向异性的影响函数引入到强度准则中。各向异性插值函数由联合不变量A和组构各向异性参数Δ构成。利用各向异性修正函数修正原始各向同性峰值强度Mp和相变强度Mpt,

    Mp=Mpexp[g(A,Δ)] ,Mpt=Mptexp[g(A,Δ)] ,} (3)

    式中,g(A,Δ)为各向异性修正函数,该函数采用与Pietruszczak等[18-19]的研究相同的形式,

    g(A,Δ)=c1Δ(1+A)=c˜1(1+A), (4)

    其中,c1为材料参数。对于各向同性的砂样,Δ= 0,g(A,Δ) ≡ 0。对于0<Δ<1的组构各向同性砂样,g(A, Δ)随AΔ的改变而改变。在大多数方法中,组构各向异性的程度可能不便于测量。为了解决这个困难,组合项c1被认为是一个单一的材料参数,因此只需要校准而不是确切的值。

    另外,Vaid等[20]的三轴试验和单剪试验结果的比较表明,剪切刚度退化的主要原因是由于单剪阶段,主应力旋转过程中产生的应力各向异性。Yang等[21-23]通过引入非共轴塑性应变模拟了这种退化。参考各向异性影响函数g(A,Δ)修正强度准则的方法[14],将各向异性影响函数g(A,Δ)式(4))引入到塑性硬化方程中,即对控制q/pεdp曲线斜率的参数kp进行修正,来实现单剪过程中由主应力轴偏转引起的剪切刚度弱化现象,如各向异性校正函数如下所示:

    kp=kpexp[c˜2(1+A)] (5)

    式中 对于各向同性的砂样,Δ=0,g(A,Δ)c˜2(1+A)≡0,kp=kpexp[c˜2(1+A)]=0;对于组构各向同性砂样,0<Δ<1,c˜2(1+A)≠0,其中c˜2是材料参数,其数值大小可控制剪切刚度kp的变化程度。

    根据ABAQUS/Explicit所提供的用户材料子程序接口VUMAT[24],将SIMSAND模型写入有限元计算程序,其基本流程如图1所示。

    图  1  ABAQUS/Explicit分析的基本流程
    Figure  1.  Flow chart of explicit finite element analysis based on ABAQUS/Explicit

    ABAQUS/Explicit计算程序中,首先进行节点的位移计算,通过内力矢量I和外力矢量P的动力平衡方程,求解节点的在t时刻下的加速度u¨(t)

    Mu¨(t) =P(t)I(t) (6)

    节点的位移u可以通过时间中心积分法求得

    u˙(t+Δt2)=u˙(tΔt2)+(Δt(t+Δt)+Δt(Δt))2u¨(t) ,u(t+Δt)=u(t)+u˙(t+Δt2)Δt(t+Δt) ,} (7)

    式中,Δt为时间增量,其中单元的应变增量Δε可通过节点的位移增量Δu=Δt(t+Δt)u˙(t+Δt2)算得。

    然后通过用户材料子程序VUMAT(SIMSAND本构关系),采用切面算法[8, 25-26]计算单元的应力增量Δσ,如图2所示,其中塑性算子为

    图  2  切面法塑性修正迭代计算
    Figure  2.  Schematic diagram of a general cutting plane algorithm
    dλ=f(σn+dσ,κ*)fσDgσ+fκ*κ*εpgσ*, (8)

    式中,f为屈服函数,g为塑性势函数,D为弹性矩阵,κ*为硬化参数,εp为属性应变。

    假设应变增量Δε完全为弹性变形,弹性试算后的应力点为σn。当试算点σ¯n+1在屈服面外时f >0,此时进入塑性修正,此过程需要跟新应力点σ和硬化参数κ*,

    σ=σDdλgσ  ,κ*=κ*+κ*εpdλgσ*  ,} (9)

    式中,εp为屈服面向外发展过程中得到的塑性应变增量。若此时应力点第一次塑性修正后依然没与屈服面重合,需要进行再次塑性修正,跟新应力点和硬化参数,屈服面向外扩展。这个过程不断重复直至修正后的应力点和屈服面重合。

    为了使用改进模型验证有限元方法,使用商用GDS单剪设备进行了一系列单剪试验,其位移传感器,量程为50 mm,分辨率可达0.001 mm。可选的荷载范围为5 kN或10 kN,其荷重传感器精度可达到0.1%的传感器量程。试验材料是枫丹白露标准砂(NE34),为法国岩土工程应用的参考材料,其颗粒的平均粒径为0.21 mm,级配如图3所示,不均匀性系数为1.53,最小孔隙比为0.510,最大孔隙比为为0.882,密度为2.65 g/cm3。伺服控制系统能够在水平和垂直两个方向上传导应力或应变以控制加载路径。图3总结了该装置的基本特征:①圆柱形试样的直径D0 = 70 mm,高度H0 = 25 mm(高宽比H0/D0 = 0.36);②将试样密封在厚度为0.2 mm的软丁基薄膜中,其自身置于具有70 mm内径的刚性圆形聚四氟乙烯涂覆的薄环(每个1 mm)上,保持恒定的横截面积但允许单剪变形;③在试样的顶部和底部,砂与粗烧结不锈钢板接触,以防止基座和试样之间的任何滑动。

    图  3  GDS单剪试验
    Figure  3.  Simple shear tests

    表2所示,进行了一系列单调加载恒定体积和恒定正应力单剪试验。在K0固结后,通过将试样的高度限制在恒定值来确保体积恒定条件(相当于“不排水条件”)。值得注意的是,所有试验均在干燥条件下完成。同时,恒定法向应力条件(相当于“排水条件”)通过固定K0固结后试样顶部的法向应力来实现。所有的试验结果都被用于有限元分析。

    表  2  枫丹白露砂单调加载的单剪试验列表
    Table  2.  Summary of monotonic simple shear tests on fontainebleau sand
    编号加载方式e0(初始)e′(K0固结后)σn0 /kPa
    S-1法向应力恒定0.700.691104
    S-2法向应力恒定0.700.688208
    S-3法向应力恒定0.700.678416
    S-4法向高度恒定0.680.666208
    S-5法向高度恒定0.680.654416
    下载: 导出CSV 
    | 显示表格

    图4为单剪试验的三维模型。即该模型的尺寸与图4所示GDS单剪设备的尺寸相同。由于对称性,只建立一半的圆柱体试样模型。如图4所示,在ABAQUS中通过C3D8R单元对圆柱体试样进行划分,共有11000个单元。金属环和顶/底板是圆柱柱试样的边界,用刚体模拟。金属环的内径为70 mm,厚度为1 mm,而顶板和底板的直径均为70 mm,厚度为2.5 mm。总共有29个刚性环堆叠用以限制试样的横向位移,并且使用两个刚性(顶部和底部)平板来限制其竖向位移。在圆柱体试样和刚性边界(环和摩擦基座)之间施加无摩擦的接触。试样的顶部和底部表面分别连接到顶部和底部刚性板,以防止它们和试样之间发生任何滑动。

    图  4  单剪试验的三维有限元建模
    Figure  4.  Three-dimensional FEM model for simple shear tests

    单剪试验的模拟过程分为K0固结和剪切两个过程。在第一步(K0固结)中,分别在顶板上施加等于0.2,0.4,0.8 kN的力以产生分别对应于104,208,416 kPa的初始法向应力σn0。然后,通过对底板施加0.05 mm/s的位移来模拟剪切过程,同时将顶板在x,y方向上的位移固定为零。整个模拟过程花了100 s的时间产生了20%的剪应变。模拟的边界条件被认为与基于物理实验室的GDS简单剪切试验类似。恒定法向应力简单剪切试验通过保持顶板上的恒定法向力来模拟;通过在剪切过程中将试样保持在恒定高度来模拟恒定体积简单剪切试验。

    位移和力的测量方法与实际室内试验相似。有效法向应力σn和剪应力τ分别通过校准底板上测得的垂直和水平力来获得。体积应变εvH/H0,通过测量圆柱试样高度的变化获得(其中ΔH为试样高度的改变量,H0为剪切开始时的试样高度)。剪应变γ= Δd/H0(其中Δd是水平位移)。分别监测x轴(Path-1),y轴(Path-2)和25 mm半径圆(Path-3)的应力状态,来评估剪切过程中的应力不均匀性,如图4(c),(d)所示。

    为验证SINSAND模型的可应用性,选用了法国常用的枫丹白露标准砂的系列三轴试验进行模拟[27-28]。砂土的模型参数:K0=100,ν=0.25,n=0.51,ϕc=33.2,eref=0.811,λ=0.055,ξ=0.46,kp =0.0022,Ad =0.39,np =1.9,nd =4,c˜1=0.18,c˜1=2.0。其中弹性相关参数通过拟合各向同性压缩试验得到,临界状态和塑性颗粒咬合效应参数均采用Jin等[29-31],提出的参数优化识别方法获得。具体标定如图5所示。

    图  5  用于参数标定的枫丹白露砂三轴试验
    Figure  5.  Simulated results of triaxial tests using determined parameters for Fontainebleau NE34 sand

    通过有限元模拟初始状态(σn0= 104 kPa,e0 = 0.7,S-1)为恒定正应力条件下的单剪试验来确定枫丹白露砂的各向异性参数(c˜1,c˜2)。此外,还进行了类似的没有考虑各向异性的模拟(即c˜1=0和c˜2= 0),以便强调在考虑参数c˜1,c˜2时模型的性能。图6为模拟结果与试验结果的比较,结果表明引入参数c˜1= 0.18和c˜2= 2.0的模拟结果与试验结果一致。因此,参数(c˜1= 0.18和c˜2= 2.0)可用于后续的单剪试验模拟。

    图  6  用于参数标定的枫丹白露砂单剪试验
    Figure  6.  Calibration of anisotropic parameters

    为了验证改进模型的有限元分析的性能,基于单剪试验的足尺三维模型(图4),分别模拟了表2中所示的4个附加的单剪试验。图7为模拟结果和试验结果的比较,其模拟曲线中的数值波动是由有限元显式计算的应力波传递所造成。恒定体积试验是通过固定试样高度来实现的,即ΔH = 0(试样高度变化为零),Δr = 0(试样径向变化为零)。在恒定体积剪切的初始阶段,颗粒重新排列出现压缩密实(即剪缩)过程,此时由于体积应变恒定,砂样首先出现法向应力变小的趋势,随后当应力比达到相变状态时,土体开始发生剪胀,法向应力逐渐增加并达到最终的极限状态。因此可以基于足尺三维单剪模型来获得恒定正应力和恒定体积条件下的单调加载特性(剪缩或剪胀)。由于室内直接试验研究困难[3],可以结合有限元模拟结果来分析试样的应力不均匀性。

    图  7  三维单剪试验模拟结果
    Figure  7.  Simulated results by three-dimensional simple shear tests

    图8(a)为在恒定正应力单剪试验(S-3)下试样顶部表面沿着三条路径的垂直应力分布的演变。这3条路径先前在图4中进行了描述,这些路径用于评估沿x轴(路径-1),y轴(路径-2)和半径为25 mm的圆(路径-3)的应力不均匀性。在剪切过程中,应力不均匀性随剪切应变水平的增加而增加,最低和最高垂直应力位于侧向(右侧和左侧)边界。结果还表明,垂直应力在试样的中心区域几乎可以保持恒定。

    图  8  垂直有效应力不均匀分布
    Figure  8.  Nonuniform distribution of vertical effective stress

    恒定体积试验的垂直应力分布(S-5)的结果也在图8(b)中给出。与恒定法向应力条件不同,由于恒定体积的边界条件,垂直应力首先会由于低剪切应变水平相对减小(剪缩),然后增加(剪胀)。与恒定法向应力试验相比,可以发现垂直分布的类似分布特征,最低和最高值也位于横向边界中并且相当均匀地位于中心区域中。为了研究侧向边界互补剪应力缺失引起的应力不均匀特性,根据网格化之后的有限元模型汇总了11000个高斯点。图9为对于恒定正应力和恒定体积单剪试验在剪切过程中基于有限元中的所有高斯点的竖向应力σz的正态分布:

    图  9  恒法向应力试验的垂直应力的正态分布(S-3:σn=416 kPa)
    Figure  9.  Normal distribution of vertical stress under constant normal stress tests
    f(x|μ,σ2)=12πσ2exp((xμ)22σ2), (10)

    式中,f(x)为可能的密度函数,μ为试样x分布的平均至;σ为标准差,其中试样x对应竖向应力σz

    对于恒定正应力单剪试验,平均值μ对应于平均垂直应力σz,其保持稳定并接近初始法向应力(即σn =416 kPa)。偏差σ对应于应力不均匀程度,随着剪切应变的增加,应力逐渐增大,如图9所示。另外,如图10所示,也可以在恒定体积单剪试验中找到类似的应力不均匀性演变(标准偏差随剪切应变γ增加)。

    图  10  恒体积试验的垂直应力的正态分布(S-5:σn=416 kPa)
    Figure  10.  Normal distribution of vertical stress under constant volum tests

    根据对恒定法向应力单剪试验(S-3)的模拟结果,在图11中绘制了剪切应力τ的演变(相应于有限元模拟中的σxz),塑性偏应变εdp和各向异性变量A。对于早期的剪切阶段(应变γ小于或等于1%),试样中的剪切应力、塑性偏应变和各向异性变量A相当均匀,特别是在中央或中部区域。对于应变大于5%的剪切阶段,剪应力峰值逐渐从左上角和右下角扩展到中心区。在试样中可以找到塑性偏应变的对角线区域。此外,由于侧边界和试样之间没有摩擦力,在剪切过程中,侧向边界处的剪切应力和偏应力几乎为零,这可能是导致应力/应变不均匀性根本的原因,从而最终导致逐渐破坏。这些模拟当中的剪应力和塑性应变的对角线累积区域的趋势与Budhu[32]的射线照片所记录的实际物理破坏区域相似。

    图  11  恒定法向应力为416 kPa试验的连续单剪过程
    Figure  11.  Profiles of successive simple shearing process under constant normal stress at σn= 416 kPa

    图12绘制了恒定体积单剪(S-5)试验后的剪切应力τ,塑性偏应变εdp和各向异性变量A。类似于通过恒定法向应力试验获得的图形(S-3),从左上角到右下角形成了累积塑性偏应变的对角区域。由于试样的高度恒定的边界条件,对角区塑性偏应变的积累比恒定正应力条件下更明显。如图11(c)图12(c)所示,各向异性变量A也沿试样的对角线区域积聚。

    图  12  法向应力为416 kPa的恒定体积试验的连续单剪过程
    Figure  12.  Profiles of successive simple shearing process under constant volume at σn = 416 kPa

    图1给出了有限元建模中所有高斯点的3个模拟状态变量(τ,εdpA)的不均匀性特性的概率分析。可以看到,正态分布的演变对于剪应力和塑性偏应变是相似的,因为偏差随应变水平而增大,对应于剪切应力或塑性应变的不均匀性在大大增加时。此外,对于如图13,14所示的各向异性变量A,分布的偏差在低剪切应变水平下变化很大,如γ=0.5%,这意味着对剪应变的增加不敏感(超过1%)。

    图  13  基于恒定法向应力试验过程中的概率分布(S-3:σn=416 kPa)
    Figure  13.  Probability distribution with increase of shear strain based on constant normal stress
    图  14  恒定体积试验过程中的概率分布(S-5:σn=416 kPa)
    Figure  14.  Probability distribution based on constant volume

    本文考虑砂土初始各向异性的影响,提出了一个模拟全尺寸单剪试验的数值模拟方法。

    该方法能引入横观各向异性联合不变量的参数,将模型从原始各向同性强度扩展到各向异性强度,由此获得了单剪条件下固有的各向异性特性。

    考虑组构各向异性的SIMSAND模型子程序写入有限元模型,通过建立等尺寸的三维有限元模型,通过建立等尺寸的三维有限元模型,可以很好的模拟枫丹白露砂的GDS单剪试验。

    本研究可以提高对单剪条件下试样剪切特性的理解,为分析试样剪切过程中应力应变不均匀性提供有效的计算工具。

  • 图  1   ABAQUS/Explicit分析的基本流程

    Figure  1.   Flow chart of explicit finite element analysis based on ABAQUS/Explicit

    图  2   切面法塑性修正迭代计算

    Figure  2.   Schematic diagram of a general cutting plane algorithm

    图  3   GDS单剪试验

    Figure  3.   Simple shear tests

    图  4   单剪试验的三维有限元建模

    Figure  4.   Three-dimensional FEM model for simple shear tests

    图  5   用于参数标定的枫丹白露砂三轴试验

    Figure  5.   Simulated results of triaxial tests using determined parameters for Fontainebleau NE34 sand

    图  6   用于参数标定的枫丹白露砂单剪试验

    Figure  6.   Calibration of anisotropic parameters

    图  7   三维单剪试验模拟结果

    Figure  7.   Simulated results by three-dimensional simple shear tests

    图  8   垂直有效应力不均匀分布

    Figure  8.   Nonuniform distribution of vertical effective stress

    图  9   恒法向应力试验的垂直应力的正态分布(S-3:σn=416 kPa)

    Figure  9.   Normal distribution of vertical stress under constant normal stress tests

    图  10   恒体积试验的垂直应力的正态分布(S-5:σn=416 kPa)

    Figure  10.   Normal distribution of vertical stress under constant volum tests

    图  11   恒定法向应力为416 kPa试验的连续单剪过程

    Figure  11.   Profiles of successive simple shearing process under constant normal stress at σn= 416 kPa

    图  12   法向应力为416 kPa的恒定体积试验的连续单剪过程

    Figure  12.   Profiles of successive simple shearing process under constant volume at σn = 416 kPa

    图  13   基于恒定法向应力试验过程中的概率分布(S-3:σn=416 kPa)

    Figure  13.   Probability distribution with increase of shear strain based on constant normal stress

    图  14   恒定体积试验过程中的概率分布(S-5:σn=416 kPa)

    Figure  14.   Probability distribution based on constant volume

    表  1   SIMSAND模型的基本本构方程

    Table  1   Basic constitutive equations of SIMSAND

    组成部分本构方程
    弹性准则˙εeij=1+υ3K(12υ)σijυ3K(12υ)σkkδij K=K0pat(2.97e)2(1+e)(ppat)n 
    屈服面f=qpH 
    塑性势面gp=Ad(Mptqp), gq=1 
    硬化参数H=Mpεpdkp+εpd 
    临界状态参数ec=erefλ(ppat)ξ tanϕp=(ece)nptanϕμ   tanϕpt=(ece)ndtanϕμ 
    三维强度标准Mp=6sinϕp3sinϕp[2c411+c41+(1c41)sin3θ]14, c1=3sinϕp3+sinϕpMpt=6sinϕpt3sinϕpt[2c421+c24+(1c42)sin3θ]14 , c2=3sinϕpt3+sinϕpt 
    下载: 导出CSV

    表  2   枫丹白露砂单调加载的单剪试验列表

    Table  2   Summary of monotonic simple shear tests on fontainebleau sand

    编号加载方式e0(初始)e′(K0固结后)σn0 /kPa
    S-1法向应力恒定0.700.691104
    S-2法向应力恒定0.700.688208
    S-3法向应力恒定0.700.678416
    S-4法向高度恒定0.680.666208
    S-5法向高度恒定0.680.654416
    下载: 导出CSV
  • [1]

    DABEET A, Discrete Element Modeling of Direct Simple Shear Response of Granular Soils and Model Validation Using Laboratory Tests[D]. Vancouver: University of British Columbia, 2014.

    [2]

    WIJEWICKREME D, SRISKANDAKUMAR S, BYRNE P. Cyclic loading response of loose air-pluviated Fraser River sand for validation of numerical models simulating centrifuge tests[J]. Canadian Geotechnical Journal, 2005, 42(2): 550-561. doi: 10.1139/t04-119

    [3]

    BUDHU M. Nonuniformities imposed by simple shear apparatus[J]. Canadian Geotechnical Journal, 1984, 21(1): 125-137. doi: 10.1139/t84-010

    [4]

    WANG B, POPESCU R, PREVOST J H. Effects of boundary conditions and partial drainage on cyclic simple shear test results—a numerical study[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2004, 28(10): 1057-1082. doi: 10.1002/nag.377

    [5]

    GROGNET M. The Boundary Conditions in Direct Simple Shear Tests: Developments for Peat Testing at Low Normal Stress[M]. Delft: Delft University of Technology, 2011.

    [6]

    DOHERTY J, FAHEY M. Three-dimensional finite element analysis of the direct simple shear test[J]. Computers and Geotechnics, 2011, 38(7): 917-924. doi: 10.1016/j.compgeo.2011.05.005

    [7] 程马遥, 金银富, 尹振宇, 等. 改进DE-TMCMC法及其在高级模型参数识别上的应用[J]. 岩土工程学报, 2019, 41(12): 2281-2289. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201912020.htm

    CHENG Ma-yao, JIN Yin-fu, YIN Zhen-yu, et al. Improved DE-TMCMC method and its application in high-level model parameter identification[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(12): 2281-2289. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201912020.htm

    [8] 吴则祥, 金银富, 季慧, 等. 易破碎砂土地基中“平底桩”贯入数值模拟分析[J]. 岩土力学, 2017, 38(增刊2): 330-336. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2017S2050.htm

    WU Ze-xiang, JIN Yin-fu, JI Hui, et al. Numerical simulation analysis of "flat-bottomed pile" penetration in easily broken sand foundation[J]. Rock and Soil Mechanics, 2017, 38(S2): 330-336. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2017S2050.htm

    [9]

    JIN Y F, WU Z X, YIN Z Y, et al. Estimation of critical state-related formula in advanced constitutive modeling of granular material[J]. Acta Geotechnica, 2017: 1-23.

    [10]

    WU Z X, YIN Z Y, JIN Y F, et al. A straightforward procedure of parameters determination for sand: a bridge from critical state based constitutive modelling to finite element analysis[J]. European Journal of Environmental and Civil Engineering, 2017: 1-23.

    [11]

    YIN Z Y, JIN Z, KOTRONIS P, et al. Novel SPH SIMSAND-based approach for modeling of granular collapse[J]. International Journal of Geomechanics, 2018, 18(11).

    [12]

    YAO Y P, KONG Y X. Extended UH model: Three- dimensional unified hardening model for anisotropic clays[J]. Journal of Engineering Mechanics, 2011, 138(7): 853-866.

    [13]

    GAO Z, ZHAO J. A non-coaxial critical-state model for sandaccounting for fabric anisotropy and fabric evolution[J]. International Journal of Solids and Structures, 2017, 106: 200-212.

    [14]

    GAO Z, ZHAO J. Efficient approach to characterize strength anisotropy in soils[J]. Journal of Engineering Mechanics, 2012, 138(12): 1447-1456. doi: 10.1061/(ASCE)EM.1943-7889.0000451

    [15]

    ODA M, NAKAYAMA H. Yield function for soil with anisotropic fabric[J]. Journal of Engineering Mechanics, 1989, 115(1): 89-104. doi: 10.1061/(ASCE)0733-9399(1989)115:1(89)

    [16]

    WANG C C. A new representation theorem for isotropic functions: an answer to Professor GF Smith's criticism of my papers on representations for isotropic functions[J]. Archive for Rational Mechanics and Analysis, 1970, 36(3): 166-197. doi: 10.1007/BF00272241

    [17]

    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. doi: 10.1061/(ASCE)1090-0241(2002)128:10(868)

    [18]

    PIETRUSZCZAK S, MROZ Z. Formulation of anisotropic failure criteria incorporating a microstructure tensor[J]. Computers and Geotechnics, 2000, 26(2): 105-112. doi: 10.1016/S0266-352X(99)00034-8

    [19]

    PIETRUSZCZAK S, MROZ Z. On failure criteria for anisotropic cohesive‐frictional materials[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2001, 25(5): 509-524. doi: 10.1002/nag.141

    [20]

    VAID Y, SIVATHAYALAN S. Static and cyclic liquefaction potential of Fraser Delta sand in simple shear and triaxial tests[J]. Canadian Geotechnical Journal, 1996, 33(2): 281-289. doi: 10.1139/t96-007

    [21]

    YANG Y, YU H. 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

    [22]

    YANG Y, YU H. Numerical simulations of simple shear with non-coaxial soil models[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2006, 30(1): 1-19. doi: 10.1002/nag.468

    [23]

    YANG Y, YU H-S. Numerical aspects of non-coaxial model implementations[J]. Computers and Geotechnics, 2010, 37(1): 93-102.

    [24]

    Hibbitt, Karlsson, Sorensen. ABAQUS/Explicit: User's Manual[R]. Vol. 1. Providence: Dassault Systemes Simulia Corp, 2001.

    [25] 李舰, 蔡国庆, 尹振宇. 适用于弹黏塑性本构模型的修正切面算法[J]. 岩土工程学报, 2020, 42(2): 253-259. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202002008.htm

    LI Jian, CAI Guo-qing, YIN Zhen-yu. Modified section algorithm for elasto-viscoplastic constitutive model[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(2): 253-259. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202002008.htm

    [26] 杨杰, 尹振宇, 黄宏伟, 等. 面向边界面模型的切面算法扩展[J]. 岩土力学, 2017, 38(12): 3436-3444. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201712006.htm

    YANG Jie, YIN Zhen-yu, HUANG Hong-wei, et al. Extension of tangent surface algorithm for boundary surface model[J]. Rock and Soil Mechanics, 2017, 38(12): 3436-3444. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201712006.htm

    [27]

    ANDRIA-NTOANINA I, CANOU J, DUPLA J. Caractérisation mécanique du sable de Fontainebleau NE34 à l’appareil triaxial sous cisaillement monotone[J]. Laboratoire Navier-Géotechnique. CERMES, ENPC/LCPC, 2010.

    [28]

    GAUDIN C, SCHNAID F, GARNIER J. Sand characterization by combined centrifuge and laboratory tests[J]. International Journal of Physical Modelling in Geotechnics, 2005, 5(1): 42-56. doi: 10.1680/ijpmg.2005.050104

    [29]

    JIN Y F, YIN Z Y, SHEN S L, et al. Selection of sand models and identification of parameters using an enhanced genetic algorithm[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2016, 40(8): 1219-1240.

    [30]

    JIN Y F, YIN Z Y, SHEN S L, et al. Investigation into MOGA for identifying parameters of a critical-state-based sand model and parameters correlation by factor analysis[J]. Acta Geotechnica, 2016, 11(5): 1131-1145.

    [31]

    JIN Y F, YIN Z Y, SHEN S L, et al. A new hybrid real-coded genetic algorithm and its application to parameters identification of soils[J]. Inverse Problems in Science and Engineering, 2016: 1-24.

    [32]

    BUDHU M. Failure state of a sand in simple shear[J]. Canadian Geotechnical Journal, 1988, 25(2): 395-400.

  • 期刊类型引用(4)

    1. 刘红,陈怡馨,刘汉龙,孙增春,肖杨. 土体热力学性质试验研究进展. 中国科学:技术科学. 2024(01): 1-14 . 百度学术
    2. 王彤,高北,蔡汉竹,吴则祥,莫品强,扈萍,庄培芝. 层状砂土地基中静压桩贯入过程挤土效应数值模拟. 济南大学学报(自然科学版). 2024(05): 609-618+633 . 百度学术
    3. 朱俊高,陈鸽,王涛,夏勇,彭文明,罗启迅. 一种适用于细粒土的新型叠片式单剪仪研制与应用. 岩土工程学报. 2024(12): 2668-2674 . 本站查看
    4. 王学涛,王立忠,洪义,高智伟. 砂土各向异性临界状态模型及砂质海床板锚承载特性评价. 岩土工程学报. 2023(11): 2346-2356 . 本站查看

    其他类型引用(0)

图(14)  /  表(2)
计量
  • 文章访问数:  292
  • HTML全文浏览量:  31
  • PDF下载量:  308
  • 被引次数: 4
出版历程
  • 收稿日期:  2020-09-03
  • 网络出版日期:  2022-12-02
  • 刊出日期:  2021-05-31

目录

/

返回文章
返回