共查询到20条相似文献，搜索用时 187 毫秒
1.
单粒组密砂剪切带的直剪试验离散元数值分析 总被引：7，自引：1，他引：6
为分析砂土的剪切力学特性,采用离散元商业软件PFC2D对单粒组密砂在直剪试验中出现的剪切带进行了数值分析。改进传统的分条带观察方式,采用矩形格分割法观察试样内部不同区域的变形特性。对离散元商业软件进行了二次开发,实现了大小主应力及其应力主方向角的可视化,分析了试样内部的应力偏转情况。同时,以纯转动率和颗粒速度等微观变量为中心,观察了试样内部颗粒的运动状态,解析了直剪试验中砂土剪切带形成的微观机制。研究表明,直剪试验中应变局部化区域集中在剪切面附近的一个条带内。通过转动场和速度场的分析可知,剪切带宽度约为10～15倍的砂土平均粒径,带内颗粒转动明显,带边缘出现大的速度和转动变化梯度。对组构和接触力分布的分析可知,剪切过程中粒间接触点和接触力的主轴方向发生了相一致的偏转,偏转后的主轴方向为60°左右。 相似文献
2.
颗粒形状是影响粗粒土密实度、力学与渗流等特性的主要因素之一。为了分析颗粒形状对粗粒土剪切特性的影响,采用离散元法生成4种不同形状的颗粒组,进行粗粒土直剪试验模拟与剪切宏细观响应研究,得出了颗粒形状对剪应力-剪位移、体应变-剪位移的影响,分析了粗粒土剪切应力、应变特性与剪胀特性。通过分析剪切带厚度、颗粒旋转量值、平均接触数、孔隙率及接触力系等宏细观参量的演化规律,研究颗粒形状在宏细观尺度上对粗粒土的影响。研究表明:异形颗粒间的咬合自锁作用大于纯圆颗粒,试样的抗剪强度有随形状系数的减小而增大的趋势。试样颗粒在外荷载作用下发生运动,应变主要表现在颗粒运动剧烈、剪胀幅度较大的剪切带内。颗粒形状系数F减小,试样的初始平均接触数增加,内摩擦角φ增大,剪切带内孔隙率增量越大,剪胀幅度越大。剪切过程中强力链聚集于剪切带内并起骨架作用,随着形状系数的减小,力链长度在0～5所占百分比呈增大趋势;剪切带内强力链的数目随着形状系数的减小而增加,峰值含量在30%～35%之间。 相似文献
3.
填埋场衬垫系统中,土与土工膜界面剪切强度较低,易造成失稳破坏。目前国内外学者主要采用室内试验对土与土工膜界面的宏观剪切特性进行研究,而对界面剪切特性的细观研究较少。为了从细观角度研究土与土工膜界面的剪切特性,本文采用EsyS-particle程序对土工膜与土界面直剪试验进行了离散元数值模拟分析。采用摩擦接触模型模拟砂土;采用黏结模型颗粒模拟土工膜,通过紧密排列土工膜颗粒以模拟土工膜的光滑表面。通过室内拟合试验,选取和校准材料的细观参数。分析结果表明,离散元模型能较好的模拟界面应力-应变关系;剪切带的厚度约为两倍平均土颗粒直径;剪切带中的土颗粒发生较大位移,孔隙比增大,而剪切带之外的土颗粒位移和孔隙比变化较小;随着剪切位移的增加,颗粒间接触力逐渐向左端集中,力链方向由垂直逐渐倾斜。 相似文献
4.
《岩土力学》2017,(5):1472-1480
对于岩土类的颗粒材料,在特定的应变加载路径下会发生非局部化的失稳现象,此时应力状态处于Mohr-Coulomb屈服面内,试样整体急剧失稳。采用颗粒离散元方法,研究抗转动特性对颗粒材料在等比例应变加载路径下宏、细观力学特性的影响。模拟发现,较为松散的试样更易发生分散性失稳,此时颗粒集合体的应力-应变状态满足Hill材料失稳准则。采用考虑颗粒转动的接触模型进行离散元模拟,通过改变颗粒间接触的转动摩擦系数,从宏观和细观层面分析等比例应变加载路径中颗粒材料的稳定性。颗粒抗转动能力的增强可以降低材料发生分散性失稳的可能性,随着转动摩擦系数的增加,应力路径由应变软化逐渐转为应变硬化,原本会发生分散性失稳的松散颗粒集合体表现出与密实颗粒集合体相似的宏观力学特性;颗粒集合体的内部结构表现出相应的细观作用机制,转动摩擦系数的增加有效地抑制了颗粒转动,虽然降低了颗粒体系的配位数,但增加了颗粒之间的接触力,增强了颗粒体系力链结构的稳定性和各向异性,形成稳定的结构持续抵抗外荷载的施加,从而试样整体不会形成松散的接触状态而失去稳定性。 相似文献
5.
接触模型的宏?细观参数标定是成功使用离散元方法的关键。在离散元的接触模型中线性接触模型与抗转线性接触模型均可用于模拟砂性土的力学行为,其中抗转线性接触模型在模拟密砂的剪胀性方面具备优势。采用抗转线性接触模型对室内密实砂土三轴试验进行了离散元模拟,验证了抗转线性接触模型的可靠性;进而系统分析了颗粒间摩擦系数、刚度比和抗转动系数等细观参数与砂土峰值内摩擦角、残余内摩擦角、峰值剪胀角等宏观参数的相关关系并进行了验证;揭示了偏应力作用下,细观参数对密实砂土试样内部剪切带宽度与倾角变化的影响规律,提出了考虑剪胀角的剪切带倾角经验公式。通过研究建立了抗转线性接触模型宏?细观参数的量化关系并给出了标定参数的具体流程图,提出了快速标定宏观参数的方法并应用实例进行了验证,为采用抗转线性接触模型精准模拟密实砂土的力学特性提供依据。 相似文献
6.
为研究颗粒形态对砂土抗剪强度的影响，制备了5组不同粒径的球形石英砂和颗粒形态不规则的实验室标准砂，并进行一系列的直剪试验；离散元法能有效地模拟砂土的离散性和不均匀性，基于高性能离散元分析软件MatDEM，通过二次开发，建立直剪试验三维模拟器，对直剪试验过程进行了数值模拟，并对比分析了试验和数值模拟结果。数据表明:(1)石英砂试样的自然休止角与其抗剪强度有着正相关性；(2)球形石英砂较实验室标准砂能更快达到峰值强度，但其值相对较低；标准砂存在着一定的黏聚力；(3)在直剪过程中，直剪盒中部颗粒的位移较周围颗粒不同，形成了明显的剪切带，并且球形石英砂试样的剪切带内颗粒的位移更小；(4)砂土直剪试验中配位数与剪应力有着良好的正相关性。本文结果表明，在相同的荷载条件下不规则颗粒形态的石英砂的抗剪强度要明显大于球形石英砂。 相似文献
7.
三维变形体离散元法能够自动检索接触关系,并对具有不规则形状的粗粒土颗粒和结构物进行有限差分网格离散,因此,具有模拟离散-连续耦合问题的先天优势。采用随机模拟技术生成粗粒土三维数值试样,基于变形体离散元进行粗粒土与结构接触面特性的数值试验,研究了不同接触面粗糙程度的接触面力学特性,对比了粗粒土与结构物在单剪和直剪状态下的接触面力学特性,从宏观和细观两个层面分析了数值试验结果。结果表明,数值试验能较好地反映粗粒土与结构接触面的力学特性,其剪应力-相对剪切位移曲线与试验结果规律相似;接触面粗糙程度对接触面的强度和变形特性影响较大,其机制在于剪切对试样的扰动程度不同;直剪和单剪状态下试样剪应力-相对剪切位移均为双曲线,单剪试验的初始剪切刚度低于直剪试验,两种试验得到接触面抗剪强度指标比较接近。 相似文献
8.
采用空心圆柱扭剪仪对干燥TJ-1模拟月壤试样进行了应力主轴固定的定向剪切试验及不同偏应力比时主应力方向的纯旋转试验。从应力-应变关系角度出发,研究主应力方向、偏应力比对其各向异性的影响,并探讨了上述因素对TJ-1模拟月壤非共轴性的影响。试验结果表明：主应力方向和偏应力比对TJ-1模拟月壤的各向异性均有显著影响;应力主轴旋转引起的非共轴现象比定向剪切时明显,且偏应力比较小时非共轴角随应力主轴旋转呈先减小后增加的趋势,偏应力比较大时非共轴角一直减小直至试样破坏时非共轴现象消失。上述成果可弥补干砂试样非共轴领域的研究空白,并可望为将来月球上基础设施的修建提供技术支持。 相似文献
9.
针对冻土的各向异性特性,基于线性黏结接触模型,建立了能够反映冻土各向异性特性的修正线性黏结接触模型,并通过C++语言生成供颗粒流程序PFC3D调用的离散元本构子程序DLL。首先对单一接触进行了拉伸、直接剪切测试,通过对比数值与理论结果,验证了考虑各向异性影响的冻土修正线性黏结接触模型的正确性。此外,模拟了不同温度条件下的冻土三轴压缩试验,并与试验得到的应力-应变曲线进行对比,结果表明,所提出的修正线性黏结接触模型对冻土具有较好的适用性。基于标定后的模型细观参数,开展了一系列的三轴压缩离散元数值模拟,利用模拟结果探讨了虚拟弱面法向倾角对冻土的应力-应变曲线特征、强度及抗剪强度指标的影响,并分析了有效配位数、细观组构量的演化规律。研究结果可为冻土各向异性宏-细观力学特性提供数值基础。 相似文献
10.
砂土力学性质的细观模拟 总被引：23，自引：9，他引：14
基于颗粒流理论,对砂土的室内双轴试验和砂土剪切带形成与发展进行了数值模拟,基本再现了砂土试样应力-应变关系。主要研究了颗粒粒径、颗粒摩擦系数等细观参数变化时试样宏观性质的变化情况。对比分析了室内试验和颗粒流数值模拟试验的砂土剪切带特征。同时对比研究了模型试样颗粒粒径、颗粒刚度和摩擦系数等细观参数的变化,对剪切带的形成与发展过程的影响。结果表明：通过颗粒流数值模型试验可以有效模拟砂土剪切带的形成与发展机理。 相似文献
11.
3D analysis of kinematic behavior of granular materials in triaxial testing using DEM with flexible membrane boundary 总被引：2，自引：2，他引：0
This paper describes the constitutive behavior and particle-scale kinematics of granular materials in three-dimensional (3D) axisymmetric triaxial testing using discrete element method (DEM). PFC3D code was used to run the DEM simulations using a flexible membrane boundary model consisting of spherical particles linked through flexible contact bonds. The overall deformation behavior of the specimen was then compared with the specimen with rigid boundary and experimental measurements. Computed tomography was used to track the evolution of particle translation and rotation within a laboratory triaxial specimen in 3D. The DEM model of the flexible membrane specimen successfully predicted the stress–strain behavior when compared with laboratory experiment results at different confining pressures. The DEM results showed that the rigid specimen applies a uniform deformation and leads to non-uniformities in the confining stress along the particle-boundary interface in the lateral direction. In contrast, the flexible specimen better replicates the uniformly applied confining stress of a laboratory triaxial experiment. The 3D DEM simulations of the specimen with flexible membrane overpredicted particle translation and rotation in all directions when compared to a laboratory triaxial specimen. The difference between the particle translation and rotation distributions of DEM specimens with rigid and flexible membrane is almost negligible. The DEM specimen with flexible membrane produces a better prediction of the macroscopic stress–strain behavior and deformation characteristics of granular materials in 3D DEM simulations when compared to a specimen with rigid membrane. Comparing macroscale response and particle-scale kinematics between triaxial simulation results of rigid versus flexible membrane demonstrated the significant influence of boundary effects on the constitutive behavior of granular materials. 相似文献
12.
This paper presents a numerical investigation into mechanical behavior and strain localization in methane hydrate (MH) bearing sediments using the distinct element method (DEM). Based on the results of a series of laboratory tests on the bonded granules idealized by two glued aluminum rods and the available experimental data of methane hydrate samples, a pressure and temperature dependent bond contact model was proposed and implemented into a two-dimensional (2D) DEM code. This 2D DEM code was then used to numerically carry out a series of biaxial compression tests on the MH samples with different methane hydrate saturations, whose results were then compared with the experimental data obtained by Masui et al. [9]. In addition, stress, strain, void ratio and velocity fields, the distributions of bond breakage and averaged pure rotation rate (APR) as well as the evolution of strain localization were examined to investigate the relationships between micromechanical variables and macromechanical responses in the DEM MH samples. The numerical results show that: (1) the shear strength increases as methane hydrate saturation SMH increases, which is in good agreement with the experimental observation; (2) the strain localization in all the DEM MH samples develops with onset of inhomogeneity of void ratio, velocity, strain, APR, and distortion of stress fields and contact force chains; and (3) the methane hydrate saturation affects the type of strain localization, with one shear band developed in the case of 40.9% and 67.8% methane saturation samples, and two shear bands formed for 50.1% methane saturation sample. 相似文献
13.
In this paper, the Discrete Element Method (DEM) is employed to numerically explore the response of hollow cylinder specimens of granular soils under complex stress paths. Two series of numerical tests are conducted to clarify the effects of the principal stress direction α and the intermediate principal stress through the b-value on the mechanical response of granular materials. The effects of α and b-value on the non-coaxiality of the principal stress and the principal plastic strain increment directions are investigated. It is observed that b-value and α significantly affect the non-coaxial behavior of granular materials. Finally, the results are discussed and compared with those obtained from physical laboratory tests. 相似文献
14.
The normal and shear strains obtained in torsion shear tests may be interpreted in two different ways to gain insight into (1) the type of plastic potential to be employed in hardening plasticity stress–strain models, and (2) the coincidence in physical space of the plastic strain increment direction with the stress direction during principal stress rotation. Thirty-four drained torsion shear tests were performed on hollow cylinder specimens of Santa Monica Beach sand deposited by dry pluviation. Twenty-six tests were performed on tall specimens with height of 40 cm, and eight tests were performed on short specimens with height of 25 cm to investigate the effect of the specimen height on the soil behavior in hollow cylinder specimens. Each test was conducted with the same, constant inside and outside confining pressure, σr, thus tying the value of b = (σ2–σ3)/(σ1– σ3) to the inclination, β, of the major principal stress. The directions of strain increment vectors at failure are compared with the directions corresponding to associated and non-associated flow. The relation between the directions of major principal strain increment and major principal stress during rotation of principal stress axes in physical space are investigated. 相似文献
15.
We investigate the quasi‐static simple shear flow of a two‐dimensional assembly of cohesionless particles using discrete element method (DEM) simulations. We focus on the unsteady flow regime where the solid would experience significant evolution of stresses, mobilised shear strength and dilation. We construct the DEM model using a discretised‐wall confined granular cell where the apparent boundary is allowed to dilate or contract synchronously with the confined solid. A rather uniform simple shear field is achieved across the whole assembly, which benefits rheological studies in generalising constitutive laws for continuum methods. We examine two aspects of the simple shear behaviour: macroscopic stress and strain rate evolution, particularly the non‐coaxiality between the principal directions of the two; and micromechanics such as evolution of fabric. For an initially anisotropic specimen sheared under constant normal pressure, the direction of principal stress rotates towards that of the principal strain rate, gradually reducing the degree of non‐coaxiality from about 45° to fluctuating around 0°. The rate in approaching coaxiality is slower in samples with larger initial porosity, stress ratio and mean stress. Generally, a faster rate in approaching coaxiality in simple shear is observed in a more dilatant sample, which often shows a larger degree of mobilised fabric anisotropy, suggesting the possible important role of instantaneous internal friction angle. The evolution of principal fabric direction resembles that of the principal stress direction. © 2013 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons, Ltd. 相似文献
16.
实际交通荷载作用下,路基土单元内的竖向应力和水平应力大小不断发生变化,剪应力幅值和方向也不断变化,从而导致土体中的应力路径呈现出主应力轴连续旋转的现象。通过GDS空心圆柱扭剪仪模拟类似交通荷载作用下的应力路径,开展不同围压和不同循环应力比下的主应力轴连续旋转试验,旨在研究在交通荷载类轴向纯压缩条件下主应力轴方向连续旋转时循环应力比与围压对原状软黏土的强度、累积应变、回弹应变、软化等因素的影响。试验结果表明：随着孔压的不断累积,原状饱和软黏土试样逐渐软化,轴向模量和剪切模量均随着循环应力比和围压的增加而逐渐降低,并在主应力轴旋转一定的循环次数后达到稳定。当循环应力比较小时,轴向和剪切应力-应变滞回曲线均呈线性,不同主应力轴循环旋转次数下的轴向和剪切应力-应变滞回曲线近乎重合。随着主应力轴循环旋转次数的增加,轴向和剪切应力-应变滞回曲线越来越表现出明显的非线性,不同循环次数下试样的轴向和剪切应力-应变滞回圈不再重合且滞回圈逐渐向x轴倾斜。随着循环应力比的增加,在主应力轴连续旋转初期,轴向模量和剪切模量迅速衰减,且随着循环次数的增加而达到稳定,并且试样的轴向模量和剪切模量达到稳定时的主应力轴连续旋转的循环次数随循环应力比和围压的增大而不断增大。 相似文献
17.
利用配备高精度应变传感器的空心圆柱扭剪仪对重塑黄土开展了一系列定向剪切试验,剪切过程中保持平均主应力p、大主应力方向角? 以及中主应力系数b不变。共完成了包括3组不同大主应力方向角和1组不同中主应力系数的4组试验,共21个试样。对定向剪切条件下重塑粉土的应力-应变-强度以及小应变（<0.1%）刚度特性进行了研究,分析了不同大主应力方向角及中主应力系数的影响。试验研究结果表明,重塑黄土存在显著的强度各向异性特性,沿? =0°方向的强度最大,而最小强度方向因中主应力系数的不同而不同。? =0°时,b=0.25的强度最大,此时试样接近平面应变条件。小应变条件下重塑黄土的剪切刚度也因剪切时?和b的不同而不同,剪应变小于0.02%和大于2%时,刚度基本保持不变,初始剪切刚度很大,当剪应变大于0.02%后,随着剪应变的增加而迅速衰减,衰减曲线可以采用指数函数进行较好的拟合。 相似文献
18.
J. Antonio H. Carraro 《Geomechanics and Geoengineering》2017,12(3):169-180
The high cost of offshore infrastructure provides continuous encouragement for optimisation of design practices. Development of a more rational method to interpret results from simple shear tests with cell pressure confinement can reduce costs and improve reliability of offshore infrastructure. This paper addresses a commonly overlooked issue affecting design parameter selection: specimen shape varies from right cylinder to oblique cylinder after loading along a single shearing direction. Thus, horizontal stresses are not always equal to the cell pressure and their magnitude varies throughout the specimen’s lateral surface. An analysis is proposed that accounts for changing specimen geometry and lateral surface area during shearing and for the actual effect of cell pressure during testing. The analysis also describes how the intermediate principal stress can be assessed. Test results for medium dense silica sand are interpreted following de Josselin de Jong’s alternative shearing mechanism hypothesis. Conventional interpretation methods yield conservative design parameters for this soil. Failure states develop when the intermediate principal effective stress is halfway between major and minor principal effective stresses. Typical results for the soil tested show equipment performance meets standard direct simple shear requirements for shear strain rate, vertical stress and specimen height control. 相似文献
19.
采用基于PFC2D接触黏结模型的离散元数值模型,探讨含石量变化对土石混合体剪切特性的影响。利用室内大型直剪试验对表征土石混合体细观力学性质的模型参数进行标定,模拟分析了4种不同含石量土石混合体在4种不同围压作用下的剪切特性。试验结果表明:相同法向应力作用下,含石量越高,峰值剪应力越大,达到峰值剪应力时的剪应变越大,软化后强度越高,剪胀特性越强,剪切"跳跃"现象越明显。土石混合体摩擦角随含石量增高而增大;黏聚力随含石量增高而减小,含石量超过50%后基本不变。相同含石量的土石混合体,法向应力越高,峰值剪应力越大,应变软化特性越强。土石混合体剪切过程中以克服摩擦能和应变能为主,动能变化几乎为零。 相似文献
20.
Torsional tests are frequently used in material testing to determine shear stress–strain relationships. In the tests, torque is applied to a solid or hollow cylindrical specimen, whose twist angle is recorded as the response. While the raw test results are readily available in terms of torque–twist sets, converting them into intrinsic stress–strain relationships is not always straightforward, owing to the nonlinear stress–strain behavior and the non‐uniform radial distribution of the stress and strain fields. It follows that the approaches currently used to interpret torsional data are essentially empirical based. The present paper addresses this data interpretation issue from a more rigorous perspective. The paper includes two parts: the first part concerns hollow cylindrical specimens under torsion, where an analytical‐based recursive procedure is developed; and the second part deals with solid torsional specimens, for which an analytical solution given by Nadai is directly employed. Examples pertaining to these two cases are examined and the resulting shear stress–strain curves are compared with those based on a number of popular empirical methods. The proposed analytical approaches can serve as a benchmark for checking the accuracy and validity of those empirical propositions. It is worth noting that, for brittle materials exhibiting softening behavior, only the analytical method gives the correct interpretation. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献