首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 390 毫秒
1.
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.  相似文献   

2.
It is normally accepted that materials inside the shear band undergo severe rotation of the principal stress direction, which causes non‐coaxiality between the principal stress and principal plastic strain rate. However, classical plasticity flow theory implicitly assumes that the principal stress and the principal plastic strain rate are coaxial; thus, it may not correctly predict the onset of the shear band. In addition, classical continuum does not contain any internal length scales; as a result, it cannot provide a reasonable shear band thickness. In this study, the original vertex non‐coaxial plastic model based on the classical continuum is extended to the Cosserat continuum. The corresponding codes are implemented via the interface of the user defined element subroutine in ABAQUS. Through a simple shear test, the effectiveness of the user's codes is verified. Through a uniaxial compression test, the influence of non‐coaxiality on the onset, the orientation, and the thickness of the shear band is investigated. Results show that the onset of the shear localization is delayed, and the thickness of the shear band is widened when the non‐coaxial degree increases, while the orientation of the shear band is little affected by the non‐coaxial degree. In addition, it is found that the non‐coaxiality can weaken the micro‐polar effect to some extent; nonetheless, the Cosserat non‐coaxial model still has its advantage over the classical non‐coaxial model in capturing the pre‐bifurcation as well as the post‐bifurcation behaviors of strain localization. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

3.
In this article, we evaluate geomechanics of fluid injection from a fully penetrating vertical well into an unconsolidated formation confined with stiff seal rocks. The coupled behavior of an isotropic, homogeneous sand layer is studied under injection pressures that are high enough to induce plasticity yet not fracturing. Propagation of the significant influence zone surrounding the injection borehole, quantified by the extent of the plastic domain in the elasto‐plastic model, is examined for the first time. First, a new fully coupled axisymmetric numerical model is developed. A comprehensive assessment is performed on pore pressures, stresses/strains, and failure planes during the entire transient period of an injection cycle. Results anticipate existence of five distinctive zones in terms of plasticity state: liquefaction at wellbore; two inner plastic domains surrounding the wellbore, where failure occurs along two planes and major principal stress is in vertical direction; remaining of the plastic domain, where formation fails along one plane and major principal stress is in radial direction; and a non‐plastic region. Failure mechanism at the wellbore is found to be shear followed by liquefaction. Next, a novel methodology is proposed based on which new weakly coupled poro‐elasto‐plastic analytical solutions are derived for all three stress/strain components. Unlike previous studies, extension of the plastic zone is obtained as a function of injection pressure, incorporating plasticity effects on the subsequent elastic domain. Solutions, proven to be a good approximation of numerical simulations, offer a huge advantage as the run time of coupled numerical simulations is considerably long. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

4.
This paper presents an elasto‐plastic model for non‐linear analyses of cement‐treated sand. Various laboratory tests were systematically carried out to investigate the pre‐peak and post‐peak behaviours of a cement‐treated sand. On the basis of these experimental results, the new model was built within the framework of a relatively simple elasto‐plastic theory. Two failure criteria are employed to express tensile and shear failure characteristics observed in the experimental results of the cement‐treated sand. The proposed model can describe strain‐hardening and strain‐softening responses under both failure modes. In the strain‐softening rules, the smeared crack concept is used, and a characteristic length is considered to avoid the issue of mesh‐size dependency. Since the failure criterion and strain‐hardening/softening rules are based on the experimental evidences, the model is relatively easy to understand and the parameters used in the model have clear physical meaning. The proposed model was applied to simulate the behaviour of cement‐treated sand in various laboratory tests, allowing for a reasonable comprehensive evaluation. It was demonstrated that the proposed model is suitable for describing both the tensile and shear failure behaviours of cement‐treated sand. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

5.
The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

6.
A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

7.
Numerous constitutive models built on coaxial theory and validated under axi‐symmetric condition often describe the stress–stain relationships and predict the inceptions of shear banding in sands inaccurately under true triaxial condition. By adopting an elaborated Mohr–Coulomb yield function and using non‐coaxial non‐associated flow rule, a 3D non‐coaxial elasto‐plasticity model is proposed and validated by a series of true triaxial tests on loose sands. The bifurcation analysis of true triaxial tests on dense sands predicts the influence of the intermediate principal stress ratio on the onset of shear band accurately. The failure of soils is shown to be related to the formation of shear band under most intermediate principal stress ratio conditions except for those which are close to the axi‐symmetric compression condition. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

8.
In this paper, the non‐coaxial relation between the principal plastic strain increments and the principal stresses, which results from the internal friction in geomaterials, is analyzed, and the phenomenon of the unbalanced development of plastic flow in two conjugate directions is discussed. A non‐coaxial, unbalanced plastic flow model for Coulomb frictional materials is developed and used to determine the orientation of shear band in geomaterials. It is shown that the unbalanced index r of plastic flow has important effect on the orientation of the shear band, and the orientation determined by the conventional plastic flow theory is only a special case of the proposed model when r=0. This result soundly explains the reason that the geomaterials with the same internal friction angle and dilatancy angle can have very different shear band orientations. In addition, the difference between the intrinsic and apparent dilatancy angles is analyzed, and it is emphasized that the dilatancy angle commonly used in practice is indeed the apparent dilatancy angle. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

9.
The shear behavior at the interface between the soil and a structure is investigated at the macroscale and particle‐scale levels using a 3‐dimensional discrete element method (DEM). The macroscopic mechanical properties and microscopic quantities affected by the normalized interface roughness and the loading parameters are analyzed. The macro‐response shows that the shear strength of the interface increases as the normalized roughness of the interface increases, and stress softening and dilatancy of the soil material are observed in the tests that feature rough interfaces. The particle‐scale analysis illustrates that a localized band characterized by intense shear deformation emerges from the contact plane and gradually expands as shearing progresses before stabilizing at the residual stress state. The thickness of the localized band is affected by the normalized roughness of the interface and the normal stress, which ranges between 4 and 5 times that of the median grain diameter. A thicker localized band is formed when the soil has a rough shearing interface. After the localized band appears, the granular material structuralizes into 2 regions: the interface zone and the upper zone. The mechanical behavior in the interface zone is representative of the interface according to the local average stress analysis. Certain microscopic quantities in the interface zone are analyzed, including the coordination number and the material fabric. Shear at the interface creates an anisotropic material fabric and leads to the rotation of the major principal stress.  相似文献   

10.
3D finite strain analyses and kinematic vorticity measurements were carried out on the Loghon Anticline within the HP‐LT Sanandaj–Sirjan metamorphic belt (Neyriz area, SW Iran). Rƒ/φ and Fry methods were used on the strain markers (e.g. deformed fossils) to interpret geometric relationships between the fold axis, strain ellipsoid axes and shear zone boundaries. The results indicate the predominance of prolate strain in the anticline. Quantitative kinematic analyses show that the Wk parameter is 0. 67 ± 0. 06 (i.e. pure‐shear dominated non‐coaxial flow). This study quantitatively supports the establishment of a dextral transpressive system, which is responsible for the development of the large‐scale right‐lateral shear zones that strike sub‐parallel to the major folds. Flexural shear combined with regional dextral‐shear is suggested to be the most common mechanism of folding in this area. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

11.
  相似文献   

12.
The Maximum Effective Moment (MEM) criterion predicts that the initial orientation of ductile shear zones and shear bands is ~55° relative to the maximum principal stress axis (σ1) and that the kinematic vorticity number (Wk) is ~0.94. These preferred orientations should be reflected in the pattern of quartz -fabrics in shear zones and shear bands. Common quartz -fabrics in plane strain can be divided into low-temperature (L) and high-temperature (H) fabrics, with each group showing three patterns. A steady flow with a constant value of Wk≈0.94 gives rise to L-1 and H-1 patterns, which are commonly characterized by a single axis girdle normal to the shear zone and a single -point maximum parallel to the shear zone.Once the conjugate set develops, L-1 and H-1 have opening angles of ~70° and ~110°, respectively. L-2 and H-2 are asymmetric patterns associated with variable deformation partitioning and vorticity values of 0< Wk<0.94. In contrast, L-3 and H-3 are symmetric patterns associated with 100% deformation partitioning and Wk=0. The opening angle in quartz -fabrics is implicitly linked to the temperature during deformation. The opening angle is ~70° at low temperature and ~110° at high temperature. However, a linear correction between the opening angle and the temperature cannot be established. During deformation partitioning, synthetic shear bands form earlier than antithetic bands and are more easily developed. This may result in opening angles of <70° for low-temperature fabrics and of >110° for high-temperature fabrics. The following criteria can be used to recognize reworked shear zones that have experienced multiple orogenic phases and changes in the stress state: 1) the initial Wk is larger or smaller than ~0.94; 2) the change in Wk is abrupt, rather than progressive; 3) inconsistent shear senses are inferred for the different phases of deformation; and 4) a negative value of Wk is found in reworked shear zones.  相似文献   

13.
The calcite mylonites in the Xar Moron-Changchun shear zone show a significance dextral shearing characteristics. The asymmetric(σ-structure) calcite/quartz grains or aggregates, asymmetry of calcite c-axes fabric diagrams and the oblique foliation of recrystallized calcite grains correspond to a top-to-E shearing. Mineral deformation behaviors, twin morphology, C-axis EBSD fabrics, and quartz grain size-frequency diagrams demonstrate that the ductile shear zone was developed under conditions of greenschist facies, with the range of deformation temperatures from 200 to 300°C. These subgrains of host grains and surrounding recrystallized grains, strong undulose extinction, and slightly curved grain boundaries are probably results of intracrystalline deformation and dynamic recrystallization implying that the deformation took place within the dislocation-creep regime at shallow crustal levels. The calculated paleo-strain rates are between 10~(–7.87)s~(–1) and 10~(–11.49)s~(–1) with differential stresses of 32.63–63.94 MPa lying at the higher bound of typical strain rates in shear zones at crustal levels, and may indicate a relatively rapid deformation. The S-L-calcite tectonites have undergone a component of uplift which led to subhorizontal lifting in an already non-coaxial compressional deformation regime with a bulk pure shear-dominated general shear. This E-W large-scale dextral strike-slip movement is a consequence of the eastward extrusion of the Xing'an-Mongolian Orogenic Belt, and results from far-field forces associated with Late Triassic convergence domains after the final closure of the Paleo-Asian Ocean.  相似文献   

14.
Reaction and deformation microfabrics provide key information to understand the thermodynamic and kinetic controls of tectono‐metamorphic processes, however, they are usually analysed in two dimensions, omitting important information regarding the third spatial dimension. We applied synchrotron‐based X‐ray microtomography to document the evolution of a pristine olivine gabbro into a deformed omphacite–garnet eclogite in four dimensions, where the 4th dimension is represented by the degree of strain. In the investigated samples, which cover a strain gradient into a shear zone from the Western Gneiss Region (Norway), we focused on the spatial transformation of garnet coronas into elongated garnet clusters with increasing strain. The microtomographic data allowed quantification of garnet volume, shape and spatial arrangement evolution with increasing strain. The microtomographic observations were combined with light microscope and backscatter electron images as well as electron microprobe (EMPA) and electron backscatter diffraction (EBSD) analysis to correlate mineral composition and orientation data with the X‐ray absorption signal of the same mineral grains. With increasing deformation, the garnet volume almost triples. In the low‐strain domain, garnet grains form a well interconnected large garnet aggregate that develops throughout the entire sample. We also observed that garnet coronas in the gabbros never completely encapsulate olivine grains. In the most highly deformed eclogites, the oblate shapes of garnet clusters reflect a deformational origin of the microfabrics. We interpret the aligned garnet aggregates to direct synkinematic fluid flow, and consequently influence the transport of dissolved chemical components. EBSD analyses reveal that garnet shows a near‐random crystal preferred orientation that testifies no evidence for crystal plasticity. There is, however evidence for minor fracturing, neo‐nucleation and overgrowth. Microprobe chemical analysis revealed that garnet compositions progressively equilibrate to eclogite facies, becoming more almandine‐rich. We interpret these observations as pointing to a mechanical disintegration of the garnet coronas during strain localization, and their rearrangement into individual garnet clusters through a combination of garnet coalescence and overgrowth while the rock was deforming.  相似文献   

15.
A lower amphibolite Alpine shear zone from the Fibbia metagranite (Gotthard Massif, Central Alps) has been studied to better understand the parameters controlling strain localization in granitic rocks. The strain gradient on the metre‐scale shows an evolution from a weakly deformed metagranite (QtzI–KfsI–AbI–BtI ± PlII–ZoI–PhgI–Grt) to a fine banded ultramylonite (QtzII–KfsII–AbII–PlII–BtII–PhgII ± Grt–ZoII). Strain localization is coeval with dynamic recrystallization of the quartzofeldspathic matrix and a modal increase in mica, at the expense of K‐feldspar. The continuous recrystallization of plagioclase during deformation into a very fine‐grained assemblage forming anastomosed ribbons is interpreted as the dominant process in the shear zone initiation and development. The shear zone initiated under closed‐system conditions with the destabilization of metastable AbI–ZoI porphyroclasts into fine‐grained (20–50 μm sized) AbII–PlII aggregates, and with minor crystallization of phengite at the expense of K‐feldspar. The development of the shear zone requires a change in state of the system, which becomes open to externally derived fluids and mass transfer. Indeed, mass balance calculations and thermodynamic modelling show that the ultramylonite is characterized by gains in CaO, FeO and H2O. The progressive input of externally derived CaO drives the continuous metamorphic recrystallization of the fine‐grained AbII–PlII aggregate into a more PlII‐rich and finer aggregate. Input of water favours the crystallization of phengite at the expense of K‐feldspar to form an interconnected network of weak phases. Thus, recrystallization of 50% of the bulk rock volume would induce a decrease of the strength of the rock that might contribute to the development of the shear zone. This study emphasizes the major role of metamorphic reactions and more particularly plagioclase on strain localization process. Plagioclase represents at least one‐third of the bulk rock volume in granitic systems and forms a stress‐supporting framework that controls the rock rheology. Therefore, recrystallization of plagioclase due to changes in P–T conditions and/or bulk composition must be taken into account, together with quartz and K‐feldspar, in order to understand strain localization processes in granites.  相似文献   

16.
The microscopic and macroscopic behaviors of assemblages of monodisperse ellipsoids with different particle shapes were studied using the discrete element method. Four samples were created with 1170 identical prolate ellipsoids. The samples were compressed isotropically to 100 kPa. Then triaxial compression tests were carried out to very large strains until the ultimate state was reached. This paper presents typical macroscopic result including stress–strain relationship and volumetric behavior. In addition, the fabric of the samples was examined at the initial state, at the peak shear strength state, and at the ultimate state. We studied the evolution of three vector‐typed micromechanical arguments with strain including the particle orientation, branch vector, and normal contact force. The normal contact force (micromechanical argument) was found to have a direct relationship with the principal stress ratio (macroscopic parameter). The angles between these vectors were also investigated. The maximum angle between vectors is related to particle shape. The results indicate that the distributions and the maximum values of these angles do not change with loading. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

17.
Unsaturated soils are highly heterogeneous 3‐phase porous media. Variations of temperature, the degree of saturation, and density have dramatic impacts on the hydro‐mechanical behavior of unsaturated soils. To model all these features, we present a thermo‐hydro‐plastic model in which the hydro‐mechanical hardening and thermal softening are incorporated in a hierarchical fashion for unsaturated soils. This novel constitutive model can capture heterogeneities in density, suction, the degree of saturation, and temperature. Specifically, this constitutive model has 2 ingredients: (1) it has a “mesoscale” mechanical state variable—porosity and 3 environmental state variables—suction, the degree of saturation, and temperature; (2) both temperature and mechanical effects on water retention properties are taken into account. The return mapping algorithm is applied to implement this model at Gauss point assuming an infinitesimal strain. At each time step, the return mapping is conducted only in principal elastic strain space, assuming no return mapping in suction and temperature. The numerical results obtained by this constitutive model are compared with the experimental results. It shows that the proposed model can simulate the thermo‐hydro‐mechanical behavior of unsaturated soils with satisfaction. We also conduct shear band analysis of an unsaturated soil specimen under plane strain condition to demonstrate the impact of temperature variation on shear banding triggered by initial material heterogeneities.  相似文献   

18.
A literature review has shown that there exist adequate techniques to obtain ground reaction curves for tunnels excavated in elastic‐brittle and perfectly plastic materials. However, for strain‐softening materials it seems that the problem has not been sufficiently analysed. In this paper, a one‐dimensional numerical solution to obtain the ground reaction curve (GRC) for circular tunnels excavated in strain‐softening materials is presented. The problem is formulated in a very general form and leads to a system of ordinary differential equations. By adequately defining a fictitious ‘time’ variable and re‐scaling some variables the problem is converted into an initial value one, which can be solved numerically by a Runge–Kutta–Fehlberg method, which is implemented in MATLAB environment. The method has been developed for various common particular behaviour models including Tresca, Mohr–Coulomb and Hoek–Brown failure criteria, in all cases with non‐associative flow rules and two‐segment piecewise linear functions related to a principal strain‐dependent plastic parameter to model the transition between peak and residual failure criteria. Some particular examples for the different failure criteria have been run, which agree well with closed‐form solutions—if existing—or with FDM‐based code results. Parametric studies and specific charts are created to highlight the influence of different parameters. The proposed methodology intends to be a wider and general numerical basis where standard and newly featured behaviour modes focusing on obtaining GRC for tunnels excavated in strain‐softening materials can be implemented. This way of solving such problems has proved to be more efficient and less time consuming than using FEM‐ or FDM‐based numerical 2D codes. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

19.
The paper investigates the incipience of shear band with an incrementally non‐linear constitutive equation. Necessary conditions for the emergence of shear band are derived. The lower bound solution is obtained by taking the strain rate inside and outside the shear band into consideration. Numerical results of localized bifurcation for general stress and strain are presented and compared with experiments. In the principal stress space, the stresses at the onset of shear band form a surface, which is partially enclosed by the failure surface for homogeneous straining. The significance of the analysis for identification of the material parameters and verification of the constitutive model against experiments is discussed. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

20.
This paper investigates the possibility of interpreting progressive shear failure in hard soils and soft rocks as the result of shear propagation of a pre‐existing natural defect. This is done through the application of the principles of fracture mechanics, a slip‐weakening model (SWM) being used to simulate the non‐linear zone at the tips of the discontinuity. A numerical implementation of the SWM in a computation method based on the boundary element technique of the displacement discontinuity method (DDM) is presented. The crack and the non‐linear zone at the advancing tip are represented through a set of elements, where the displacement discontinuity (DD) in the tangential direction is determined on the basis of a friction law. A residual friction angle is assumed on the crack elements. Shear resistance decreases on elements in the non‐linear zone from a peak value at the tip, which is characteristic of intact material, to the residual value. The simulation of a uniaxial compressive test in plane strain conditions is carried out to exemplify the numerical methodology. The results emphasize the role played by the critical DD on the mechanical behaviour of the specimen. A validation of the model is shown through the back analysis of some experimental observations. The results of this back analysis show that a non‐linear fracture mechanics approach seems very promising to simulate experimental results, in particular with regards to the shear band evolution pattern. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号