Micro-scale behavior of granular materials during cyclic loading

This study presents the micro-scale behavior of granular materials under biaxial cyclic loading for differ- ent confining pressures using the two-dimensional （2D） discrete element method （DEM）. Initially, 8450 ovals were generated in a rectangular frame without any overlap. Four dense samples having confining pressures of 15, 25, 50, and 100 kPa were prepared from the initially generated sparse sample. Numeri- cal simulations were performed under biaxial cyclic loading using these isotropically compressed dense samples. The numerical results depict stress-strain-dilatancy behavior that was similar to that observed in experimental studies. The relationship between the stress ratio and dilatancy rate is almost indepen- dent of confining pressures during loading but significantly dependent on the confining pressures during unloading. The evolution of the coordination number, effective coordination number and slip coordina- tion number depends on both the confining pressures and cyclic loading. The cyclic loading significantly affects the microtopology of the granular assembly. The contact fabric and the fabric-related anisotropy are reported, as well. A strong correlation between the stress ratio and the fabric related to contact normals is observed during cyclic loading, irrespective of confining pressures.     

DEM investigation on the effect of sample preparation on the shear behavior of granular soil

The effect of initial fabric anisotropy produced by sample preparation on the shear behavior of granular soil is investigated by performing discrete element method (DEM) simulations of fourteen biaxial tests in drained conditions. Numerical test specimens are prepared by three means: gravitational deposition, multi-layer compression, and isotropic compression, such that different initial inherent soil fabrics are created. The DEM simulation results show that initial fabric anisotropy exerts a considerable effect on the shear behavior of granular soil, and that the peak stress ratio and peak dilatancy increase with an increase in the fabric index a_n that is estimated from the contact orientations. The stress–dilatancy relationship is found to be independent of the initial fabric anisotropy. The anisotropy related to the contact orientation and contact normal force accounts for the main contribution to the mobilized friction angle. Also, the occurrence of contractive shear response in an initial shearing stage is accompanied by the most intense particle rearrangement and microstructural reorganization, regardless of the sample preparation method. Furthermore, the uniqueness of the critical state line in e–log p′ and q–p′ plots is observed, suggesting that the influence of initial fabric anisotropy is erased at large shear strains.     

基于细观拓扑结构演化的颗粒材料剪胀性分析

剪胀性是包括岩土材料在内的摩擦性颗粒材料的重要特征之一,其形成机制与颗粒体系内部拓扑结构的演化有关.基于颗粒体系细观数据,可对颗粒体系内部的拓扑结构特征及演化进行分析,进而建立拓扑演化与宏观剪胀变形之间的联系.采用离散单元法,根据密实、中密和松散摩擦性颗粒材料双轴试验的宏微观数据,从拓扑参量演化及接触网络拓扑变化所引起...     

饱和砂土的循环边界面本构模

饱和砂土在循环加卸载过程中,土体发生显著的组构变化和塑性变形的累积。试验现象分析表明,循环过程中,砂土的首次与再循环塑性模量既有区别又有联系。因此,论文引入考虑组构变化的剪胀内变量,并提出循环塑性模量的关系式,真实描述上述两种现象;进而基于已有本构模型的基础上,建立饱和砂土的循环边界面本构模型。最后,将饱和砂土的模型预测结果与三轴试验结果进行验证对比,得到较好地吻合,这表明该模型能够合理反映饱和砂土循环加卸载的变形行为。     

Temperature dependent bulge test for elastomers

The bulge test is a particularly convenient testing method for characterizing elastomers under biaxial loading. In addition, it is convenient to utilize this test for validating material models in simulation due to the heterogeneous strain field induced during inflation. During the bulge test the strain field for elastomers covers uniaxial tension at the border to pure shear and equibiaxial tension at the pole. Elastomeric materials exhibit a hyperelastic material behavior, with a dependency on temperature and loading rate. The temperature effect on the mechanical behavior during biaxial loading is considered in the present study. A bulge test setup combined with a temperature chamber is developed in order to characterize this effect, and an exemplary temperature dependent characterization of a poly(norbornene) elastomer is performed with this setup. The equibiaxial stress–strain curves measured at 60 °C, 20 °C and −20 °C are presented.     

含抗转能力散粒体的宏微观力学特性数值分析

颗粒间抗转因素对散粒体宏微观特性的影响十分显著。本文采用离散单元法,将蒋明镜（2005年）等提出的抗转动接触模型植入离散元分析软件PFC2D中,模拟了大量刚性边界下的砂土双轴试验,研究了压缩过程中试样的宏观力学性质及部分微观参量的变化情况,用以定性分析颗粒间抗转作用对散粒体的宏微观力学特性的影响。研究表明：颗粒间的抗转...     

DEM investigation of strain behaviour and force chain evolution of gravel–sand mixtures subjected to cyclic loading

The strain characteristic and load transmission of mixed granular matter are different from those of homogeneous granular matter. Cyclic loading renders the mechanical behaviours of mixed granular matter more complex. To investigate the dynamic responses of gravel–sand mixtures, the discrete element method (DEM) was used to simulate the cyclic loading of gravel–sand mixtures with low fines contents. Macroscopically, the evolution of the axial strain and volumetric strain was investigated. Mesoscopically, the coordination number and contact force anisotropy were studied, and the evolution of strong and weak contacts was explored from two dimensions of loading time and local space. The simulation results show that increasing fines content can accelerate the development of the axial strain and volumetric strain but has little effect on the evolution of contact forces. Strong contacts tend to develop along the loading boundary, presenting the spatial difference. Weak contacts are firstly controlled by confining pressure and then controlled by axial stress, while strong contacts are mainly controlled by axial stress throughout the whole cyclic loading. Once compression failure occurs, the release of axial stress causes the reduction of strong contact proportion in all local regions. These findings are helpful to understand the dynamic responses of gravel–sand mixtures, especially in deformation behaviours and the Spatio-temporal evolution of contact forces.     

A micromechanics-based elastoplastic damage model for granular materials at low confining pressure

This paper is devoted to the formulation of a micromechanics-based constitutive model for granular materials under relatively low confining pressure. The constitutive formulation is performed within the general framework of homogenization for granular materials. However, new rigorous stress localization laws are proposed. Some local constitutive relations are established under the consideration of irreversible thermodynamics. Macroscopic plastic deformation is obtained by considering local plastic sliding in a limit number of families of contact planes. The plastic sliding at each contact plane is described by a non-associated plastic flow rule, taking into account pressure sensitivity and normal dilatancy. Nonlinear elastic deformation related to progressive compaction of contacts is also taken into account. Material softening is described by involving damage process related to degradation of microstructure fabric. The proposed model is applied to some typical granular materials (sands). The numerical predictions are compared with experimental data.     

Lattice orientation effects on void growth and coalescence in fcc single crystals

Void growth and coalescence in fcc single crystals were studied using crystal plasticity under uniaxial and biaxial loading conditions and various orientations of the crystalline lattice. A 2D plane strain unit cell with one and two cylindrical voids was employed using three-dimensional 12 potentially active slip systems. The results were compared to five representative orientations of the tensile axis on the stereographic triangle. For uniaxial tension conditions, the void volume fraction increase under the applied load is strongly dependent on the crystallographic orientation with respect to the tensile axis. For some orientations of the tensile axis, such as [1 0 0] or [1 1 0], the voids exhibited a growth rate twice as fast compared with other orientations ([1 0 0], [2 1 1]). Void growth and coalescence simulations under uniaxial loading indicated that during deformation along some orientations with asymmetry of the slip systems, the voids experienced rotation and shape distortion, due mainly to lattice reorientation. Coalescence effects are shown to diminish the influence of lattice orientation on the void volume fraction increase, but noteworthy differences are still present. Under biaxial loading conditions, practically all differences in the void volume fraction for different orientations of the tensile axes during void growth vanish. These results lead to the conclusion that at microstructural length scales in regions under intense biaxiality/triaxiality conditions, such as crack tip or notched regions, the plastic anisotropy due to the initial lattice orientation has only a minor role in influencing the void growth rate. In such situations, void growth and coalescence are mainly determined by the stress triaxiality, the magnitude of accumulated strain, and the spatial localization of such plastic strains.     

Experimental investigation on strength and failure behavior of pre-cracked marble under conventional triaxial compression

Fractures in natural rocks have an important effect on the strength and failure behavior of rock mass, which are often evaluated in rock engineering practice. The theoretical evaluation of mechanical behavior of fractured rock mass has no satisfactory answer due to the role of confining pressure and crack geometry. Therefore, in this paper, conventional triaxial compression experiments were carried out to study the strength and failure behavior of marble samples with two pre-existing closed cracks in non-overlapping geometry. Based on the experimental results of a number of triaxial compression tests, the effect of crack coalescence on the axial supporting capacity and deformation property were investigated with different confining pressures. The results show that intact samples and flawed samples (marble with pre-existing cracks) have different deformation properties after peak stress, which change from brittleness to plasticity and ductility with the increase of confining pressure. The peak strength and failure mode are found depending not only on the geometry of flaw, but also on the confining pressure. The strength of flawed samples shows distinct non-linear behavior, which is in a better agreement with non-linear Hoek–Brown criterion than linear Mohr–Coulomb criterion. For a kind of rock that has been evaluated as a Hoek–Brown material, a new evaluation criterion is put forward by adopting optimal approximation polynomial theory, which can be used to confirm more precisely the strength parameters (cohesion and internal friction angle) of flawed samples. For intact samples, the marble leads to typical shear failure mode with a single fracture surface under different confining pressures, while for flawed samples, under uniaxial compression and a lower confining pressure (σ₃ = 10 MPa), tests for coarse and medium marble (the coarse and medium refer to the grain size) exhibit three basic failure modes, i.e., tensile mode, shear mode, and mixed mode (tensile and shear). Shear mode is associated with lower strength behavior. However, under higher confining pressures (σ₃ = 30 MPa), for coarse marble, the axial supporting capacity is not related to the geometry of flaw. The friction among crystal grains determines the strength behavior of coarse marble. For medium marble, the failure mode and deformation behavior are dependent on the crack coalescence in the sample. The present research provides increased understanding of the fundamental nature of rock failure under conventional triaxial compression.     

剪胀性砂土边界面模型的研究

剪胀性对于砂土,尤其是中密以及密实砂土,是一个非常显著的特性。相变线是剪胀性砂土的特征曲线,能够反映砂土的围压以及初时孔隙比对变形特性的影响。本文在边界面塑性理论的框架内,把相变状态参量引入到剪胀方程以及塑性硬化模量中,建立了一个能够描述砂土剪胀性以及循环特性的本构模型。本模型采用一套参量可以模拟不同初时孔隙比、不同围压、排水(或不排水)条件下单调(或循环)加载的应力-应变特性。验证表明本模型数值计算与试验结果相吻合。     

Incorporating the effects of fabric in the dilatant double shearing model for planar deformation of granular materials

The objective of this paper is to incorporate the effects of fabric and its evolution into the Dilatant Double Shearing Model [Mehrabadi, M.M., Cowin, S.C., 1978. Initial planar deformation of dilatant granular materials. J. Mech. Phys. Solids 26, 269–284] for granular materials in order to capture the anisotropic behavior and the complex response of granular materials in cyclic shear loading. An important consequence of considering the fabric is that one can have unequal shearing rates along the two slip directions. This property leads to the non-coaxiality of the principal axes of stress and strain rate, which is more appropriate for a material that exhibits initial and induced anisotropy. In addition, we employ a fabric-dependent elasticity tensor with orthotropic symmetry. The model developed in this paper also predicts one of the experimentally observed characteristics of granular materials: the gradual concentration of the contact normals towards the maximum principal stress direction.We implement the constitutive equations into ABAQUS/Explicit by writing a user material subroutine in order to predict the strength anisotropy of granular materials in a plane strain biaxial compression test and investigate the mechanical behavior of granular materials under the cyclic shear loading conditions. The predictions from this model show good quantitative agreement with the experiments of [Park, C.S., 1990. Anisotropy in deformation and strength properties of sands in plane strain compression, Masters Thesis, University of Tokyo; Park, C.S., Tatsuoka, F., 1994. Anisotropic strength and deformation of sands in plane strain compression. In: XIII ICSMFE, New Delhi, India; Okada, N., 1992. Energy dissipation in inelastic flow of cohesionless granular media. Ph.D. Thesis, University of California, San Diego].     

Multi-cycle deformation of semicrystalline polymers: Observations and constitutive modeling

Experimental data are reported on isotactic polypropylene in multi-cycle uniaxial tensile tests where a specimen is stretched up to some maximum strain and retracted down to the zero minimum stress, while maximum strains increase with number of cycles. Fading memory of deformation history is observed: when two samples are subjected to loading programs that differ along the first n − 1 cycles only, their stress–strain diagrams coincide starting from the nth cycle. Constitutive equations are developed in cyclic viscoelasticity and viscoplasticity of semicrystalline polymers, and adjustable parameters in the stress–strain relations are found by fitting the experimental data. Results of numerical simulation demonstrate that the model predicts the fading memory effect quantitatively. To confirm that the observed phenomenon is typical of semicrystalline polymers, experimental data are presented in tensile cyclic tests with large strains on low density polyethylene and compressive cyclic tests on poly(oxymethylene).     

Structural stability and jamming of self-organized cluster conformations in dense granular materials

We examine emergent, self-organized particle cluster conformations in quasistatically deforming dense granular materials from the perspective of structural stability. A structural mechanics approach is employed, first, to devise a new stability measure for such conformations in equilibrium and, second, to use this measure to explore the evolving stability of jammed states of specific cluster conformations, i.e. particles forming force chains and minimal contact cycles. Knowledge gained on (a) the spatial and temporal evolution of stability of individual jammed conformations and (b) their relative stability levels, offer valuable clues on the rheology and, in particular, self-assembly of granular materials. This study is undertaken using data from assemblies of nonuniformly sized circular particles undergoing 2D deformation in two biaxial compression tests: a discrete element simulation of monotonic loading under constant confining pressure, and cyclic loading of a photoelastic disk assembly under constant volume. Our results suggest that the process of self-assembly in these systems is realized at multiple length scales, and that jammed force chains and minimal cycles form the basic building blocks of this process. In particular, 3-cycles are stabilizing agents that act as granular trusses to the load-bearing force chain columns. This co-evolutionary synergy between force chains and 3-cycles proved common to the different materials under different loading conditions. Indeed, the remarkable similarities in the evolution of stability, prevalence and persistence of minimal cycles and force chains in these systems suggest that these structures and their co-evolution together form a generic feature of dense granular systems under quasistatic loading.     

Soil compaction and tires for harvesting and transporting sugarcane

Four tire types (A, block-shape tread; B, rib-shape tread; C, low-lug tread; D, high-lug tread) used to harvest and transport sugarcane were compared regarding the compaction induced to the soil. Tires were tested at three inflation pressures (207, 276, 345 kPa) and six loads ranging from 20 to 60 kN/tire. Track impressions were traced, and 576 areas were measured to find equations relating inflation pressure, load, contact surface and pressure. Contact surface increased with increasing load and decreasing inflation pressure; however, the contact pressure presented no defined pattern of variation, with tire types A and B generating lower contact pressure. The vertical stresses under the tires were measured and simulated with sensors and software developed at the Colombian Sugarcane Research Center (Cenicaña). Sensors were placed at 10, 30, 50 and 70 cm depth. Tire types A and B registered vertical stresses below 250 kPa at the surface. These two tires were better options to reduce soil compaction. The equations characterizing the tires were introduced into a program to simulate the vertical stress. Simulated and measured stresses were adjusted in an 87–92% range. Results indicate a good correlation between the tire equations, the vertical stress simulation and the vertical stress measurement.     

A thermomechanical constitutive model for phase transformations in silicon under pressure and contact loading conditions

Most of the technologically relevant abrasive machining techniques for silicon (Si) such as lapping, sawing and grinding are based on the interaction of the silicon surface with a hard particle or asperity. It has been long established that the governing deformation mechanism for Si under such contact loading conditions is stress induced phase transformation. The present work introduces a novel phenomenological constitutive model for phase transformations of silicon set up in a thermomechanical framework of broad applicability. Taking into account experimental observations as well as first principle and molecular dynamics calculations, it captures both the cd-Si → β-Si transition upon compression and the β-Si → a-Si transition upon rapid decompression, which are most relevant for indenter loading. The model was numerically implemented in analogy to incremental plasticity and successfully applied for finite-element (FE) simulations of nanoindentation.     

颗粒材料三维应力路径下的接触组构特性

颗粒材料的宏观应力变形特征与其微观接触力、组构等紧密相关.一般而言,强接触系统属于颗粒内部体系的传力结构,其对应的组构张量是影响宏观应力性质的重要因素.细观数值方法(如离散单元法)能够反映物理试验的基本规律,并且可以方便地提取宏微观数据来研究颗粒体系的应力变形机制.采用离散单元法(discrete element method,DEM)进行一系列等$p$等$b$应力路径下颗粒材料的真三轴试验,在此基础上研究了三维应力路径下颗粒材料的宏微观力学参数的演化过程、三维组构张量与应力张量多重联系以及强接触体系反映的宏观应力特征.研究表明:颗粒体系偏应力峰值状态和临界状态均存在与加载路径无关的宏微观特征;三维应力路径下组构张量与应力张量存在非共轴性,但其联合不变量演化过程表现出加载路径无关的特征;与弱接触系统的组构张量相比,强接触系统的组构张量更能反映宏观应力张量的特征;强弱接触体系的组构张量对颗粒体系宏观响应的贡献不同,其分界点存在一定取值范围,但采用平均接触力较为简单合理.      

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

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