基于PFC~(3D)的砂土直剪模拟及宏细观分析

为能更好地分析砂土在剪切过程中的体积变化及力学行为,通过PFC~(3D)颗粒流程序,按照室内实际级配建立数值试样,进行了不同垂直压力下的砂土直剪试验模拟,并将模拟结果与室内试验结果进行了对比;研究了砂土在剪切过程中的体积变化,并用两种方法分析了砂土剪切带的演变过程;从细观角度对试样颗粒的速度场及力链网络的发展变化进行了分析研究。结果表明,数值模拟结果与室内试验结果基本一致,砂土的体积变化表现为先剪缩后剪胀,剪胀量与垂直压力成反比;剪切带厚度约为11～12倍d_(50),在剪切带内颗粒的位移和欧拉角变化较大;试样内部强力链的演变较为明显,力链网络对外荷载变化的敏感性很高;剪切过程中上盒颗粒速度场方向的变化能够较好地解释砂土产生的剪胀现象。     

制样方法对初始各向异性砂土三轴剪切力学特性的影响研究

基于砂土颗粒的毛细效应特性和冻融试验原理,提出一种砂土试样制备方法,能够在0°~90°范围内任意选取沉积方向,且适用于不同颗粒级配的试样。在中围压条件下进行三轴固结排水剪切试验,研究结果表明,毛细效应、冻融相结合的制样方法能够减缓压缩剪切试验的应力–应变关系增长速度,降低试样的抗剪峰值强度,对残余强度的影响相对较弱。在其他应力路径条件下,例如拉伸、减压压缩、等应力压缩等试验中,可以不考虑装样方法的影响。建议的制样方法能够合理应用到砂土初始各向异性的研究工作中。     

INFLUENCE OF SAMPLE PREPARATION TECHNOLOGY ON SHEAR BEHAVIOR OF INHERENT ANISOTROPIC SAND IN TRIAXIAL TESTS 1)

基于砂土颗粒的毛细效应特性和冻融试验原理,提出一种砂土试样制备方法,能够在0°~90°范围内任意选取沉积方向,且适用于不同颗粒级配的试样。在中围压条件下进行三轴固结排水剪切试验,研究结果表明,毛细效应、冻融相结合的制样方法能够减缓压缩剪切试验的应力-应变关系增长速度,降低试样的抗剪峰值强度,对残余强度的影响相对较弱。在其他应力路径条件下,例如拉伸、减压压缩、等应力压缩等试验中,可以不考虑装样方法的影响。建议的制样方法能够合理应用到砂土初始各向异性的研究工作中。     

沉积方向和粒径组成对砂土力学特性的影响

运用砂土毛细效应所形成的假凝聚力制备具有不同沉积方向的试样,开展三轴固结排水剪切试验,研究粒径组成对初始各向异性砂土力学特性的影响,分析了粒径组成、沉积方向、强度特性和体积变化之间的相互关系。研究结论表明,试样剪切带的形成受到颗粒级配的影响,它是应力应变关系中陡降阶段的产生原因。沉积方向对粗砂和细砂的偏应力峰值强度均有影响,对残余强度的影响受制于颗粒级配。沉积方向与试样破坏面接近时,试样强度低、剪胀变形小且更容易破坏,粒径组成对砂土三轴剪切力学特性的影响更显著。     

循环荷载下硅藻土应变和孔压特性

针对嵊州硅藻土的高结构性及地基应力场的复杂多变，应用GDS 动三轴仪进行不同围压和不同循环应力比下的不排水循环加载试验，分析应力-应变滞回曲线演化、应变累积特性、回弹特性及孔压特性．研究表明：循环应力比及围压对原状硅藻土应变和孔压的发展规律影响明显，随着循环应力比的增加，滞回曲线由线性状态转变为非线性状态，且逐渐向x 轴倾斜，累计应变和回弹应变都随着循环应力比的增加而略微增加，但总应变基本小于2 %；残余孔压比随循环应力比增大而增大，范围为0.1~0.5．并且存在一个临界循环应力比，其值随有效围压的增大而减小，范围为0.8~2.2，当循环应力比大于该值时，试样内部产生一个破裂面迅速发生破坏．在相同循环应力比下，试样的应变随着有效围压的增大而增大．     

应力历史对颗粒材料强度和变形特征影响的离散元模拟

主要关注了颗粒材料前期所受的应力历史对其后期宏观力学响应的影响。该应力历史由一段等比例加载应力路径以及卸载垂直方向应力至与水平方向围压相同的卸载段描述。具体工作为:基于PFC2D双轴压缩数值实验,调查了应力历史对颗粒样本的强度、变形特征、细观参量如组构的影响,得出颗粒样本的偏应力-应变、体积应变曲线及其发生破坏时的名义应变云图。数值结果表明:高低围压下样本分别发生剪切破坏和弥散破坏,高围压下弹性阶段的刚度受应力历史影响较大,而低围压下样本在刚进入塑性至应力峰值点阶段的弹塑性刚度变化较为明显。随着应力历史中等比例加载系数的增大,剪胀加快,变形局部化范围有所不同;另一方面,颗粒形状的不规则性会增强颗粒材料的各向异性,导致样本强度更高。     

不同围压下C60混凝土三轴压缩过程能量分析

针对复杂应力状态下高强混凝土受压变形破坏过程中的能量演化机制问题,开展不同围压下C60高强混凝土试样三轴压缩试验,分析其在受力全过程中的变形与破坏特征。根据试验结果,探究了不同围压下高强混凝土三轴压缩过程能量演化机制。研究结果表明:围压越大,混凝土试样破坏时,其峰值应力与峰值应力对应的轴向应变越大,破坏形式由张拉破坏向剪切破坏过渡;峰值应力前,混凝土试样主要以弹性应变能储存为主,峰值应力对应的输入能密度和耗散能密度均随围压的增大而增大,且均与围压满足指数函数关系,其形状改变系数FX与轴向应力呈正比关系,体积改变系数FV与轴向应力呈反比关系;达到峰值应力时,体积改变系数FV小于形状改变系数FX;峰值应力后,主要以弹性应变能释放为主,并随着混凝土试样的破坏转化为各种形式的能量耗散。研究结果可为今后从能量角度研究高强混凝土本构关系提供有益参考。     

类岩石材料力学特性的试验及数值模拟研究

在进行多组不同配比类岩石材料单轴压缩试验和巴西试验的基础上,详细分析了石膏水泥比和石英砂含量对类岩石材料的单轴抗压强度、抗拉强度及弹性模量等力学参数的影响规律,力图找到适合模拟现场砂质泥岩的类岩石材料及配合比。利用颗粒流程序(PFC)模拟,进一步研究了高径比和围压对类岩石材料力学特性的影响。结果表明:随着石膏水泥比的增大,抗压强度和弹性模量均逐渐减小,而抗拉强度逐渐增大;随着石英粉含量的增大,抗压强度和弹性模量均先增大后减小,而抗拉强度则为先减小后增大。结合单轴压缩过程的声发射特征,揭示了裂纹扩展与声发射有密切的关系。PFC2D模拟获得的力学参数与室内试验相近,破裂模式也与实际情况相似。通过尺寸效应的研究可知试样的高径比在2.0~2.5较合理。随着围压的增大,试样的峰值强度、残余强度、峰值应变及弹性模量等力学参数均增大,且围压会改变试样的破裂模式。     

不同围压下含天然贯通节理岩石卸载破坏机理研究

二郎山特长高速公路隧道断裂带较多，受此影响，隧道开挖后岩爆灾害十分严重．从此隧道中取样制备成含天然贯通节理试件，进行常规单轴压缩室内试验，获取基本力学参数和破坏模式．本文以室内实验结果为依据，确定数值试验时所采用的岩石和节理的力学参数，利用RFPA软件开展数值试验，研究了不同围压同轴向荷载下卸掉围压诱发岩石失稳破坏机理．试验结果表明，五种不同围压(5、10、15、20、25 MPa)卸载后应力-应变曲线均呈“断崖式”跌落，试样均突然破坏，与现场岩爆情况吻合；受贯通节理影响，不同围压试样的具体破坏形态、声发射数量最大值出现时刻有所差异；在卸载阶段五种试样声发射累积数量都突然急剧增大．     

基于室内试验与数值模拟的膨胀土裂隙对强度影响规律研究

针对室内试验测试多裂隙发育膨胀土抗剪强度指标的困难,提出以室内试验与数值模拟相结合的研究方法;以不同产状单裂隙面膨胀土三轴试样在200kPa围压下的CD试验结果为样本,基于颗粒流原理,采用PFC3D软件进行仿真试验,获得了与三轴试验结果相一致的PFC3D球颗粒参数和裂隙参数;以此参数为基础,进行不同裂隙产状膨胀土的仿真试验,初步探寻了裂隙发育对广西宁明膨胀土的强度影响规律。研究发现:裂隙发育是膨胀土强度衰减的主要原因;裂隙的空间角度、不同裂隙的组合、裂隙发育程度均对膨胀土抗剪强度指标产生影响,且对黏聚力的影响较大;与无裂隙情况相比,裂隙率达14.06%时,C值降低71.2%。     

Strain localization and fabric evolution in sand

Strain localization is frequently observed in sand and is considered an important precursor related to major geohazards such as landslides, debris flow and failure of relevant geo-structures. This paper presents a numerical study on strain localization in sand, with a special emphasis on the influence of soil fabric and its evolution on the initiation and development of shear band. In particular, a critical state sand plasticity model accounting for the effect of fabric and its evolution is used in the finite element analysis of plane strain compression tests. It is found that the initiation of shear band is controlled by the initial fabric, while the development of shear band is governed by two competing physical mechanisms, namely, the structural constraint and the evolution of fabric. The evolution of fabric generally makes the sand response more coaxial with the applied load, while the structural constraint induced by the sample ends leads to more inhomogeneous deformation within the sand sample when the initial fabric is non-coaxial with the applied stress. In the case of smooth boundary condition, structural constraint dominates over the fabric evolution and leads to the formation of a single shear band. When the boundary condition is rough, the structural constraint may play a comparable role with fabric evolution, which leads to symmetric cross-shape shear bands. If the fabric is prohibited from evolving in the latter case, a cross-shape shear band pattern is found with the one initiated first by the structural constraint dominating over the second one. In all cases, significantly larger dilation and fabric evolution are observed inside the shear band than outside. The simulated shear band orientation coincides with the Roscoe’s angle for cases with high confining pressure and lies in between the Roscoe’s angle and Arthur’s angle for the low confining pressure cases.     

黄河三角洲粉质土的动模量和阻尼比试验研究

结合室内共振柱和动三轴实验,对黄河三角洲饱和原状粉质土体(粉土、粉砂、粉质粘土)动模量和阻尼比的影响因素和发展规律进行了详细的研究。研究表明,在粉粒和粘粒含量对动模量的共同影响中,粉粒含量起着举足起重的作用;侧限压力对归一化剪模比和阻尼比的影响均较显著,相比粘粒含量的影响不大。通过与Seed建议的砂土及饱和粘土的G/Gmax~γ曲线和λ~γ曲线进行对比,结果显示研究区的粉质土相比一般的砂土和饱和粘土而言,其动力变形特性更接近于砂土,但是与砂土也存在着非常明显的差异;其发展规律与其他地区沉积粉质土也较为不同,具有明显的区域性。采用修正了的Hard in-D rnevich模型和对数模型分别对G/Gmax~γ曲线和λ~γ曲线进行拟合,给出了三类粉质土的归一化动力变形G/Gmax~γ/γr关系曲线,对模型中有关参数的影响因素做了初步的探讨。     

A mesoscale investigation of strain rate effect on dynamic deformation of single-crystal copper

A combined finite element (FE) simulation and discrete dislocation dynamics (DD) approach has been developed in this paper to investigate the dynamic deformation of single-crystal copper at mesoscale. The DD code yields the plastic strain based on the slip of dislocations and serves as a substitute for the 3D constitutive form used in the usual FE computation, which is implemented into ABAQUS/Standard with a user-defined material subroutine. On the other hand, the FE code computes the displacement and stress field during the dynamic deformation. The loading rate effects on the yield stress and the deformation patterning of single-crystal copper are investigated. With the increasing of strain rate, the yield stress of single-crystal copper increases rapidly. A critical strain rate exists in each single-crystal copper block for the given size and dislocation sources, below which the yield stress is relatively insensitive to the strain rate. The dislocation patterning changes from non-uniform to uniform under high-strain-rate. The shear stresses in the bands are higher than that in the neighboring regions, which are formed shear bands in the crystal. The band width increases with the strain rate, which often take places where the damage occurs.     

基于邓肯-张的冻结粉质粘土本构模型试验研究

为研究冻结粉质粘土强度和变形特性,以沈阳地铁DK11+395联络通道处人工冻结粉质粘土为研究对象,通过冻土三轴剪切试验,研究了不同试验条件下冻结粉质粘土的强度和变形特性。试验结果表明：围压、温度和试件初始含水量是影响冻结粉质粘土强度和变形的主要因素。冻结粉质粘土偏应力-应变曲线呈应变硬化型,其破坏强度和切线弹性模量随围压和试件初始含水量的增加而增大,随温度的升高而减小,破坏偏应力比Rf取值在0.79~0.96之间。基于试验数据建立了以上述三因素为影响因子的冻结粉质粘土Duncan-Chang模型。通过回归分析,建立了模型参数a和b与围压、初始含水量和温度之间的线性回归公式。将依据模型计算的偏应力-应变曲线与对比试验曲线相比较,发现两者具有较好的吻合程度,说明建立的模型能够准确反映围压、初始含水量和温度等条件对冻结粉质粘土强度和变形的影响规律。上述研究成果为冻结粉质粘土强度和变形特性的研究提供了参考,为人工冻结法施工提供了指导,具有重要的理论意义和工程应用价值。     

金属材料平面应变拉伸变形局部化有限元分析

本文对于涉及韧性金属大变形中颈缩与剪切带断裂一类高度非线性变形局部化问题进行了弹塑性有限元数值模拟。采用改进的Ｊ２形变理论微分形式公式与交叉三角形四边形单元有限元网格，详细研究了应变硬化指数及初始表面不均匀特性的平面应变拉伸颈缩和剪切带形成的综合影响，给出此类问题的断裂机制图。     

Simulation and interpretation of diffuse mode bifurcation of elastoplastic solids

The diffuse mode bifurcation of elastoplastic solids at finite strain is investigated. The multiplicative decomposition of deformation gradient and the hyperelasto-plastic constitutive relationship are adapted to the numerical bifurcation analysis of the elastoplastic solids. First, bifurcation analyses of rectangular plane strain specimens subjected to uniaxial compression are conducted. The onset of the diffuse mode bifurcations from a homogeneous state is detected; moreover, the post-bifurcation states for these modes are traced to arrive at localization to narrow band zones, which look like shear bands. The occurrence of diffuse mode bifurcation, followed by localization, is advanced as a possible mechanism to create complex deformation and localization patterns, such as shear bands. These computational diffuse modes and localization zones are shown to be in good agreement with the associated experimental ones observed for sand specimens to ensure the validity of this mechanism. Next, the degradation of horizontal sway stiffness of a rectangular specimen due to plane strain uniaxial compression is pointed out as a cause of the bifurcation of the first antisymmetric diffuse mode, which triggers the tilting of the specimen. Last, circular and punching failures of a footing on a foundation are simulated.     

Discrete element modeling of the effect of particle size distribution on the small strain stiffness of granular soils

Discrete element modeling was used to investigate the effect of particle size distribution on the small strain shear stiffness of granular soils and explore the fundamental mechanism controlling this small strain shear stiffness at the particle level. The results indicate that the mean particle size has a negligible effect on the small strain shear modulus. The observed increase of the shear modulus with increasing particle size is caused by a scale effect. It is suggested that the ratio of sample size to the mean particle size should be larger than 11.5 to avoid this possible scale effect. At the same confining pressure and void ratio, the small strain shear modulus decreases as the coefficient of uniformity of the soil increases. The Poisson’s ratio decreases with decreasing void ratio and increasing confining pressure instead of being constant as is commonly assumed. Microscopic analyses indicate that the small strain shear stiffness and Poisson’s ratio depend uniquely on the soil’s coordination number.     

三轴试验土样剪切破坏过程的判别

本文在土工三轴试验中应用土样全表面变形数字图像测量系统,通过跟踪包裹土样的橡皮膜表面的标记点,得到加载过程中其全表面的变形和应变分布。根据土样表面的应变分布和实测的应力应变关系曲线,可以得到每一时刻土样表面的应力分布。于是,土样表面剪切破坏开始出现的时刻和剪切破坏带贯穿的时刻可分别依据应变和应力判定。依据应变,将橡皮膜上各点的轴向应变的拐点最晚出现的时刻作为剪切带完全贯穿的时刻;依据应力,将全部点中某一个或几个点最早达到峰值应力的时刻作为剪切带开始出现的时刻;土样各点相继达到应力峰值的过程就是剪切破坏的过程。以剪切破坏的起始点和剪切带贯穿的点为分界,土样剪切破坏过程可以划分成三个阶段:破坏前阶段、破坏阶段和破坏后阶段,在不同阶段土样具有不同的变形特点和变形机理。     

Numerical simulation for the simple shear test in amorphous glassy polymers

In this paper, a driving stress finite element method of elastic-plastic large deformation based on implicit time integrating algorithm and an eight-chain molecular network model is used for the numerical simulation of the simple shear test of polycarbonate (PC) materials. The simulated results are compared with experimental ones. The strain localization propagation for the shear band deformation for simple shear deformation is investigated numerically. The effects of microstructure parameters in the model on strain softening and orientation hardening of the PC are discussed in detail. Supported by the National Natural Science Foundation of China.