Modeling shear localization along granular soil–structure interfaces using elasto-plastic Cosserat continuum

The current study presents finite element simulations of shear localization along the interface between cohesionless granular soil and bounding structure under large shearing movement. Micro-polar (Cosserat) continuum approach is applied in the framework of elasto-plasticity in order to overcome the numerical problems of localization modeling seen in the conventional continuum mechanics. The effects of different micro-polar kinematic boundary conditions, along the interface, on the evolution and location of shear band are shown by the numerical results. Furthermore, shear band thickness is also investigated for its dependence on the initial void ratio, vertical pressure and mean grain size. Here, the distribution and evolution of static and kinematic quantities are the main focuses regarding infinite layer of micro-polar material during plane shearing, especially with advanced large movement of bounding structure. The influence of such movement has not been investigated yet in the literature. Based on the results obtained from this study, shear localization appears parallel to the direction of shearing. It occurs either in the middle of granular layer or near boundaries, regarding the assumed micro-polar kinematic boundary conditions at the bottom and top surfaces of granular soil layer. Narrower shear band is observed in lower rotation resistance of soil particles along the interface. It is emphasized that the displacement magnitude of bounding structure has significant effect on the distribution and evolution of state variables and polar quantities in the granular soil layer. However, continuous displacement has no meaningful effect on the thickness of shear band. Here, smooth distributions of void ratio and shear stress components are obtained within the shear band, what the other previous numerical investigations did not receive. Despite indirect linking of Lade’s model to the critical state soil mechanics, state variables tend towards asymptotical stationary condition in large shear deformation.     

Strain gradient crystal plasticity effects on flow localization

In metal grains one of the most important failure mechanisms involves shear band localization. As the band width is small, the deformations are affected by material length scales. To study localization in single grains a rate-dependent crystal plasticity formulation for finite strains is presented for metals described by the reformulated Fleck–Hutchinson strain gradient plasticity theory. The theory is implemented numerically within a finite element framework using slip rate increments and displacement increments as state variables. The formulation reduces to the classical crystal plasticity theory in the absence of strain gradients. The model is used to study the effect of an internal material length scale on the localization of plastic flow in shear bands in a single crystal under plane strain tension. It is shown that the mesh sensitivity is removed when using the nonlocal material model considered. Furthermore, it is illustrated how different hardening functions affect the formation of shear bands.     

微孔洞生长对板材拉伸过程中变形局部化的影响

本文在作者提出的含孔洞材料下限本构方程的基础上，采用了初始缺陷带模型对微孔洞生长及分布对板材拉伸过程中变形局部影响进行了，分析着重研究了细观损演化规律对变形局部化模式及临界应变的影响，并成功预测了ＡＩＳＩ４３４０钢板材拉伸试件变形局部化换稳为及失稳方向。     

Micromechanics based second gradient continuum theory for shear band modeling in cohesive granular materials following damage elasticity

Gradient theories, as a regularized continuum mechanics approach, have found wide applications for modeling strain localization failure process. This paper presents a second gradient stress–strain damage elasticity theory based upon the method of virtual power. The theory considers the strain gradient and its conjugated double stresses. Instead of introducing an intrinsic material length scale into the constitutive law in an ad hoc fashion, a microstructural granular mechanics approach is applied to derive the higher-order constitutive coefficients such that the internal length scale parameter reflects the natural granularity of the underlying material microstructure. The derivations of the required damage constitutive relationships, the strong form governing equations as well as its weak form for the second gradient model are described. The recently popularized Element-Free Galerkin (EFG) method is then employed to discretize the weak form equilibrium equation for accommodating the resultant higher-order continuity requirements and further handling the mesh sensitivity problem. Numerical examples for shear band simulations show that the proposed second gradient continuum model can produce stable, accurate as well as mesh-size independent solutions without a priori assumption of the shear band path.     

Shear band formation in thermo-viscoplastic solids under plane strain compression

The plane strain compression of a rectangular block is numerically investigated for the study of dynamic shear band development in thermo-elasto-viscoplastic materials from an internal inhomogeneity. As expected, it plays an important role in triggering the onset of shear, localization as well as thermal softening. And the competition between the strain, strain-rate hardening and thermal softening exists throughout the process. It is found that shear band develops at a 45-degree angle to the compression axis. In the light of given patterns of deformation and temperature, shear band evolution accelerated by thermal softening is retarded by the inertial effects. Interestingly, a similar temperature band is also formed along the trajectory of the localized deformation band. The calculations also show the energy evolution during the coupled thermo-mechanical process of shear band propagation. Finally, the mesh effect is discussed in terms of the numerical results from two different meshes. The project is supported by the National Natural Sciences Foundation of China.     

Shear bands in dense metallic granular materials

In this paper, the phenomenon of strain localization, i.e., shear banding, in densely distributed metallic assemblies has been studied. A discrete element methodology for analyzing metallic granular materials has been put forward. In this numerical model, elastoplastic contact, as well as friction, rolling resistance and cohesion between spheres, are explicitly taken into account. The calculations reveal that the shear banding mechanism in dense assemblies can be thought as an instability triggered by initial imperfections. Within a band, the motion, deformation and rearrangement of spheres soften the resistance of the aggregate, as these mechanisms create additional geometric imperfections. Additionally, the simulations showed that the shear-band width does not change conclusively with the friction, rolling resistance and plasticity parameters. However, cohesive strength, even in small amounts, drastically increased the shear-band width. Finally, the shear-band thickness and inclination angles are strongly dependent on the degree of initial imperfection. Whereas for a perfect assembly the shear band inclinations were consistently around 60°, more heterogeneous assemblies lead to shear band angles closer to the continuum mechanics solution, which is 45°. This was found to be in agreement with recent experimental observations.     

Plasticity size effects in voided crystals

The shear and equi-biaxial straining responses of periodic voided single crystals are analysed using discrete dislocation plasticity and a continuum strain gradient crystal plasticity theory. In the discrete dislocation formulation, the dislocations are all of edge character and are modelled as line singularities in an elastic material. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and annihilation are incorporated through a set of constitutive rules. Over the range of length scales investigated, both the discrete dislocation and strain gradient plasticity formulations predict a negligible size effect under shear loading. By contrast, under equi-biaxial loading both plasticity formulations predict a strong size dependence with the flow strength approximately scaling inversely with the void spacing. Excellent agreement is obtained between predictions of the two formulations for all crystal types and void volume fractions considered when the material length scale in the non-local plasticity model is chosen to be (about 10 times the slip plane spacing in the discrete dislocation models).     

Critical state shear behavior of the soil-structure interface determined by discrete element modeling

The interface between soil and structure can be referred to as a soil-structure system, and its behavior plays an important role in many geotechnical engineering practices. In this study, results are presented from a series of monotonic direct shear tests performed on a sand-structure interface under constant normal stiffness using the discrete element method (DEM). Strain localization and dilatancy behavior of the interface is carefully examined at both macroscopic and microscopic scales. The effects of soil initial relative density and normal stress on the interface shear behavior are also investigated. The results show that a shear band progressively develops along the structural surface as shear displacement increases. At large shear displacement a unique relationship between stress ratio and void ratio is reached in the shear band for a certain normal stress, indicating that a critical state exists in the shear band. It is also found that the thickness and void ratio of the shear band at the critical state decreases with increasing normal stress. Comparison of the DEM simulation results with experimental results provides insight into the shear behavior of a sand-structure interface and offers a means for quantitative modeling of such interfaces based on the critical state soil mechanics.     

Prediction of the initial thickness of shear band localization based on a reduced strain gradient theory

Plastic flow localization in ductile materials subjected to pure shear loading and uniaxial tension is investigated respectively in this paper using a reduced strain gradient theory, which consists of the couple-stress (CS) strain gradient theory proposed by Fleck and Hutchinson (1993) and the strain gradient hardening (softening) law (C–W) proposed by Chen and Wang (2000). Unlike the classical plasticity framework, the initial thickness of the shear band and the strain rate distribution in both cases are predicted analytically using a bifurcation analysis. It shows that the strain rate is obviously non-uniform inside the shear band and reaches a maximum at the center of the shear band. The initial thickness of the shear band depends on not only the material intrinsic length l_cs but also the material constants, such as the yield strength, ultimate tension strength, the linear hardening and softening shear moduli. Specially, in the uniaxial tension case, the most possible tilt angle of shear band localization is consistent qualitatively with the existing experimental observations. The results in this paper should be useful for engineers to predict the details of material failures due to plastic flow localization.     

Shear banding analysis of geomaterials by strain gradient enhanced damage model

Shear band localization is investigated by a strain-gradient-enhanced damage model for quasi-brittle geomaterials. This model introduces the strain gradients and their higher-order conjugate stresses into the framework of continuum damage mechanics. The influence of the strain gradients on the constitutive behaviour is taken into account through a generalized damage evolutionary law. A weak-form variational principle is employed to address the additional boundary conditions introduced by the incorporation of the strain gradients and the conjugate higher-order stresses. Damage localization under simple shear condition is analytically investigated by using the theory of discontinuous bifurcation and the concept of the second-order characteristic surface. Analytical solutions for the distributions of strain rates and strain gradient rates, as well as the band width of localised damage are found. Numerical analysis demonstrates the shear band width is proportionally related to the internal length scale through a coefficient function of Poisson’s ratio and a parameter representing the shape of uniaxial stress–strain curve. It is also shown that the obtained distributions of strains and strain gradients are well in accordance with the underlying assumptions for the second-order discontinuous shear band boundary and the weak discontinuous bifurcation theory.     

Effects of void band orientation and crystallographic anisotropy on void growth and coalescence

The effects of void band orientation and crystallographic anisotropy on void growth and linkage have been investigated. 2D model materials were fabricated by laser drilling a band of holes into the gage section of sheet tensile samples using various orientation angles with respect to the tensile axis normal. Both copper and magnesium sheets have been studied in order to examine the role of crystallographic anisotropy on the void growth and linkage processes. The samples were pulled in uniaxial tension inside the chamber of an SEM, enabling a quantitative assessment of the growth and linkage processes. The void band orientation angle has a significant impact on the growth and linkage of the holes in copper. As the void band orientation angle is increased from 0° to 45°, the processes of coalescence and linkage are delayed to higher strain values. Furthermore, the mechanism of linkage changes from internal necking to one dominated by shear localization. In contrast, the void band orientation does not have a significant impact on the void growth and linkage processes in magnesium. Void growth in these materials occurs non-uniformly due to interactions between the holes and the microstructure. The heterogeneous nature of deformation in magnesium makes it difficult to apply a coalescence criterion based on the void dimensions. Furthermore, the strain at failure does not show a relationship with the void band orientation angle. Failure associated with twin and grain boundaries interrupts the plastic growth of the holes and causes rapid fracture. Therefore, the impact of the local microstructure outweighs the effects of the void band orientation angle in this material.     

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.     

Discrete element analysis of the mechanical properties of deep-sea methane hydrate-bearing soils considering interparticle bond thickness

Due to increasing global energy demands, research is being conducted on the mechanical properties of methane hydrate-bearing soils (MHBSs), from which methane hydrate (MH) will be explored. This paper presents a numerical approach to study the mechanical properties of MHBSs. The relationship between the level of MH saturation and the interparticle bond thickness is first obtained by analyzing the scanning electron microscope images of MHBS samples, in which is the bridge connecting the micromechanical behavior captured by the DEM with the macroscopic properties of MHBSs. A simplified thermal-hydromechanical (THM) bond model that considers the different bond thicknesses is then proposed to describe the contact behavior between the soil particles and those incorporated into the discrete element method (DEM). Finally, a series of biaxial compression tests are carried out with different MH saturations under different effective confining pressures to analyze the mechanical properties of deep-sea MHBSs. The results of the DEM numerical simulation are also compared with the findings from triaxial compression tests. The results show that the macromechanical properties of deep-sea MHBSs can be qualitatively captured by the proposed DEM. The shear strength, cohesion, and volumetric contraction of deep-sea MHBSs increase with increasing MH saturation, although its influence on the internal friction angle is obscure. The shear strength and volumetric contraction increase with increasing effective confining pressure. The peak shear strength and the dilation of MHBSs increase as the critical bond thickness increases, while the residual deviator stress largely remains the same at a larger axial strain. With increasing the axial strain, the percentage of broken bonds increases, along with the expansion of the shear band.     

A second gradient elastoplastic cohesive-frictional model for geomaterials

Starting from some experimental observations on shear strength of stiff clays [3,2], a isotropic, geometrically non-linear second gradient elastoplastic model is proposed for pressure dependent, brittle geomaterials. The development of the model follows the general theory presented in [1]. Due to the internal length scale provided by the microstructure, the model is ideally suited for the analysis of failure problems in which strain localization into shear band occurs, see, e.g., [4,5].     

TOWARDS A NEW MODEL-FREE SIMULATION OF HIGH-REYNOLDS-FLOWS: LOCAL AVERAGE DIRECT NUMERICAL SIMULATION

Studies on the numerical simulation of high-Reynolds-number flows encounter difficulties due to the wide range of characteristic length and time scales existing in the flow field. These are often much smaller than the computational grid size. A new approach based on a ‘model-free ’ local average direct numerical simulation is presented which incorporates a strategy to filter the non-resolvable scales by means of an integration over the domain and to recover the contribution of the subgrid scales by using an integral formulation developed for them. The resulting weak formulation allows us to define a numerical flux that, thanks to the filtering operation, is highly accurate. Several computation test-cases concerning theoretical accuracy and the Navier–Stokes equations at high Reynolds number are carried out without using any turbulence model. The obtained accuracy for all computations confirms that this approach can be considered a valid contribution in the field of direct numerical simulation.     

A first-order strain gradient damage model for simulating quasi-brittle failure in porous elastic solids

In order to simulate quasi-brittle failure in porous elastic solids, a continuum damage model has been developed within the framework of strain gradient elasticity. An essential ingredient of the continuum damage model is the local strain energy density for pure elastic response as a function of the void volume fraction, the local strains and the strain gradients, respectively. The model adopts Griffith’s approach, widely used in linear elastic fracture mechanics, for predicting the onset and the evolution of damage due to evolving micro-cracks. The effect of those micro-cracks on the local material stiffness is taken into account by defining an effective void volume fraction. Thermodynamic considerations are used to specify the evolution of the latter. The principal features of the model are demonstrated by means of a one-dimensional example. Key aspects are discussed using analytical results and numerical simulations. Contrary to other continuum damage models with similar objectives, the model proposed here includes the effect of the internal length parameter on the onset of damage evolution. Furthermore, it is able to account for boundary layer effects.     

Cosserat连续体模型中本构参数对应变局部化模拟结果影响的数值分析

基于所发展的压力相关弹塑性Cosserat连续体模型及相应的数值方法,以一维剪切层及二维平板压缩问题为例,数值分析了Cosserat连续体模型中的本构参数Cosserat剪模、软化模量及内部长度参数对应变局部化数值模拟结果的影响.结果表明在一定取值范围内,Cosserat剪模对数值模拟结果几乎没有影响,并给出了具体数值计算时的取值范围;软化模量绝对值越大,后破坏段的荷载-位移曲线越陡,计算得到的剪切带宽度越窄;内部长度参数越大,后破坏段的荷载-位移曲线越平缓,计算得到的剪切带越宽.     

Rupture mechanisms in combined tension and shear—Micromechanics

A micromechanics model based on the theoretical framework of plastic localization into a band introduced by Rice is developed. The model consists of a planar band with a square array of equally sized cells, with a spherical void located in the centre of each cell. The periodic arrangement of the cells allows the study of a single unit cell for which fully periodic boundary conditions are applied. The micromechanics model is applied to analyze failure by ductile rupture in experiments on double notched tube specimens subjected to combined tension and torsion carried out by the present authors. The stress state is characterized in terms of the stress triaxiality and the Lode parameter. Two rupture mechanisms can be identified, void coalescence by internal necking at high triaxiality and void coalescence by internal shearing at low triaxiality. For the internal necking mechanism, failure is assumed to occur when the deformation localizes into a planar band and is closely associated with extensive void growth until impingement of voids. For the internal shearing mechanism, a simple criterion based on the attainment of a critical value of shear deformation is utilized. The two failure criteria capture the transition between the two rupture mechanisms successfully and are in good agreement with the experimental result.