DOCSLIB.ORG

Application of the Mohr-Coulomb Yield Criterion for Rocks with Multiple Joint Sets Using Fast Lagrangian Analysis of Continua 2D (FLAC2D) Software

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Application of the Mohr-Coulomb Yield Criterion for Rocks with Multiple Joint Sets Using Fast Lagrangian Analysis of Continua 2D (FLAC2D) Software energies Article Application of the Mohr-Coulomb Yield Criterion for Rocks with Multiple Joint Sets Using Fast Lagrangian Analysis of Continua 2D (FLAC2D) Software Lifu Chang 1,2,* ID and Heinz Konietzky 1 1 Institut für Geotechnik, Technische Universität Bergakademie Freiberg, Freiberg 09599, Germany; [email protected] 2 School of Mechanics & Civil Engineering, China University of Mining & Technology, Beijing 100083, China * Correspondence: [email protected]; Tel.: +86-189-1102-2674 Received: 6 February 2018; Accepted: 6 March 2018; Published: 9 March 2018 Abstract: The anisotropic behavior of a rock mass with persistent and planar joint sets is mainly governed by the geometrical and mechanical characteristics of the joints. The aim of the study is to develop a continuum-based approach for simulation of multi jointed geomaterials. Within the continuum methods, the discontinuities are regarded as smeared cracks in an implicit manner and all the joint parameters are incorporated into the equivalent constitutive equations. A new equivalent continuum model, called multi-joint model, is developed for jointed rock masses which may contain up to three arbitrary persistent joint sets. The Mohr-Coulomb yield criterion is used to check failure of the intact rock and the joints. The proposed model has solved the issue of multiple plasticity surfaces involved in this approach combined with multiple failure mechanisms. The multi-joint model was implemented into Fast Lagrangian Analysis of Continua software (FLAC) and is verified against the strength anisotropy behavior of jointed rock. A case study considering a circular excavation under uniform and non-uniform in-situ stresses is used to illustrate the practical application. The multi-joint model is compared with the ubiquitous joint model. Keywords: jointed rock masses; multi-joint constitutive model; strength anisotropy; tunnel excavation; numerical simulation 1. Introduction The behavior of a jointed rock (within this paper this term includes also jointed rock masses) is anisotropic, non-linear and stress path dependent due to the presence of discontinuities. Compared to intact rock, jointed rock generally exhibit reduced stiffness, higher permeability and lower strength. Various analytical, numerical and empirical methods are suggested in order to take into account the influence of discontinuities on the mechanical behavior of jointed rock [1–6]. In the past several decades, a lot of lab tests with samples containing a single or multiple joint sets under confined and unconfined conditions have been implemented to investigate the joint orientation effect [7–9]. With the ongoing development of computer based simulation techniques, more sophisticated models have been developed in recent years and it is possible to incorporate more aspects into the simulations [10–13]. Two numerical techniques are common in rock mechanics: continuum-based methods and discontinuum based methods. Discontinuum methods consider the rock mass as an assemblage of rigid or deformable blocks connected along discontinuities [14], and investigate the microscopic mechanisms in granular material and crack development in rocks [3,15]. Boon et al. [16] performed 2D discrete element analyses to identify potential failure mechanisms around a circular unsupported tunnel in a rock mass intersected by three independent joint sets. In continuum methods, the discontinuities are considered as smeared planes of weakness (joints) incorporated in each zone. The equivalent Energies 2018, 11, 614; doi:10.3390/en11030614 www.mdpi.com/journal/energies Energies 2018, 11, x FOR PEER REVIEW 2 of 23 tunnelEnergies in2018 a, 11rock, 614 mass intersected by three independent joint sets. In continuum methods,2 ofthe 23 discontinuities are considered as smeared planes of weakness (joints) incorporated in each zone. The equivalentcontinuum continuum method contains method three contains parts: three the equivalent parts: the continuumequivalent compliancecontinuum partcompliance (incremental part (incrementalelastic law), theelastic yield law), criterion the yield and criterion the plastic an correctionsd the plastic part. corrections part. WillWill [ [1717]] implemented implemented a a multi multi-surface-surface plasticity plasticity model model into into a a finite finite element element code code with with internal internal iterationsiterations in in case case that that several several plasticity plasticity surfaces surfaces are are touched touched at at the the same same time. time. Zhang Zhang [ [1818]] developed developed furtherfurther the the Goodman’s Goodman’s joint joint element element model model and and derived derived the the equivalent equivalent method method for for rock rock masses masses with with multimulti-set-set joints joints in in a a global global coordinate coordinate system. system. Rihai Rihai [19 [19,20,20] ]proposed proposed a a method method which which takes takes into into accountaccount the the influence influence of of micro micro moments moments on on the the behavior behavior of of the the equivalent equivalent continuum continuum based based on on the the CosseratCosserat theory. HurleyHurley etet al.al. [ 21[21]] presented presented a numericala numerical method method which which implemented implemented the the method method into intoGEODYN-L GEODYN code.-L code. This This method method uses uses a series a series of stress of stress and strain-rateand strain- updaterate update algorithms algorithms to rule tojoints rule jointsclosure, closure, slip, and slip, stress and distributionstress distribution inside the inside computational the computational cells which cells contain which multiple contain embedded multiple embeddedjoints. A model joints. with A model multi-sets with ofmulti ubiquitous-sets of jointsubiquitous with anjoints extended with an Barton extended empirical Barton formula empirical was formulaproposed was by proposed Wang [22] by and Wang pre- and[22] post-peakand pre- and deformation post-peak characteristics deformation characteristics of the rock mass of werethe rock also massconsidered. were also Many considered. models haveMany been models formulated have been to simulateformulated the to effect simulate of joint the angle effect on of the joint strength angle onand the deformability strength and of deformability rock masses, butof rock only masses, a few consider but only several a few jointconsider sets several having differentjoint sets strengthhaving differentparameters strength and their parameters effects on and failure their characteristics.effects on failure characteristics. ThisThis paper paper presents presents a a new new equivalent equivalent continuum continuum constitutive constitutive model model named named multi multi-joint-joint model model toto simulate simulate the the behavior behavior of ofjointed jointed rock. rock. This This model model is suited is suited for jointed for jointed rock containing rock containing up to three up to arbitrarythree arbitrary joint sets joint with sets corresponding with corresponding spatial spatialdistribution. distribution. The equivalent The equivalent compliance compliance matrix matrixof the rockof the mass rock is mass establis is establishedhed and the and Mohr the Mohr-Coulomb-Coulomb yield yield criterion criterion is used is used to tocheck check the the failure failure characteristicscharacteristics of thethe intactintact rock rock and and the the joints. joints. A circular A circular tunnel tunnel subjected subjected to uniform to uniform and non-uniform and non- uniformin-situ stresses in-situ isstresses used to is illustrateused to illustrate the application the application of this model. of this model. 2.2. Methodo Methodologylogy for for Equivalent Equivalent Continuum Continuum Multi Multi-Joint-Joint Model Model 2.1.2.1. Ubiquitous Ubiquitous Joint Joint and and Related Related Approaches Approaches TheThe term term ‘ubiquitous ‘ubiquitous joint’ joint’ used used by by Goodman Goodman [ [2323]] implies implies that that the the joint joint sets sets may may occur occur at at any any pointpoint in in the the rock rock mass mass without without a a fixed fixed location. location. The The joints joints are are embedded embedded i inn a a Mohr Mohr-Coulomb-Coulomb solid solid (Figure(Figure 11aa).). FigureFigure 1. 1. SchematicSchematic for for isotropic isotropic and and anisotropic anisotropic models: models: (a) ubiquitous (a) ubiquitous joint model, joint model, (b) multi (b )joint multi model, joint (model,c) anisotropic (c) anisotropic model with model transversely with transversely isotropic isotropicrock matrix rock and matrix representative and representative joint orientation joint orientation.. Failure may occur in either the rock matrix or along the joints, or both, according to the stress Failure may occur in either the rock matrix or along the joints, or both, according to the stress state, the orientation of the joints and the parameters of matrix and joints [24]. More advanced models state, the orientation of the joints and the parameters of matrix and joints [24]. More advanced models based on ubiquitous joint or strain-harding/softening ubiquitous joint (subiquitous) were developed based on ubiquitous joint or strain-harding/softening ubiquitous joint (subiquitous) were developed in recent years [13,22,25–27] as shown in Figure 1b,c. These approaches consider aspects like stiffness Energies 2018, 11, 614 3 of 23 Energies 2018, 11, x FOR PEER REVIEW 3 of 23 in recent years [13,22,25–27]
Load more

Recommended publications
  • Mechanics of Materials Plasticity
    Provided for non-commercial research and educational use. Not for reproduction, distribution or commercial use. This article was originally published in the Reference Module in Materials Science and Materials Engineering, published by Elsevier, and the attached copy is provided by Elsevier for the author’s benefit and for the benefit of the author’s institution, for non-commercial research and educational use including without limitation use in instruction at your institution, sending it to specific colleagues who you know, and providing a copy to your institution’s administrator. All other uses, reproduction and distribution, including without limitation commercial reprints, selling or licensing copies or access, or posting on open internet sites, your personal or institution’s website or repository, are prohibited. For exceptions, permission may be sought for such use through Elsevier’s permissions site at: http://www.elsevier.com/locate/permissionusematerial Lubarda V.A., Mechanics of Materials: Plasticity. In: Saleem Hashmi (editor-in-chief), Reference Module in Materials Science and Materials Engineering. Oxford: Elsevier; 2016. pp. 1-24. ISBN: 978-0-12-803581-8 Copyright © 2016 Elsevier Inc. unless otherwise stated. All rights reserved. Author's personal copy Mechanics of Materials: Plasticity$ VA Lubarda, University of California, San Diego, CA, USA r 2016 Elsevier Inc. All rights reserved. 1 Yield Surface 1 1.1 Yield Surface in Strain Space 2 1.2 Yield Surface in Stress Space 2 2 Plasticity Postulates, Normality and Convexity
    [Show full text]
  • An Effective Stress Equation for Unsaturated Granular Media in Pendular Regime
    University of Calgary PRISM: University of Calgary's Digital Repository Graduate Studies The Vault: Electronic Theses and Dissertations 2014-04-28 An Effective Stress Equation for Unsaturated Granular Media in Pendular Regime Khosravani, Sarah Khosravani, S. (2014). An Effective Stress Equation for Unsaturated Granular Media in Pendular Regime (Unpublished master's thesis). University of Calgary, Calgary, AB. doi:10.11575/PRISM/24849 http://hdl.handle.net/11023/1443 master thesis University of Calgary graduate students retain copyright ownership and moral rights for their thesis. You may use this material in any way that is permitted by the Copyright Act or through licensing that has been assigned to the document. For uses that are not allowable under copyright legislation or licensing, you are required to seek permission. Downloaded from PRISM: https://prism.ucalgary.ca UNIVERSITY OF CALGARY An Effective Stress Equation for Unsaturated Granular Media in Pendular Regime by Sarah Khosravani A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE DEPARTMENT OF CIVIL ENGINEERING CALGARY, ALBERTA April, 2014 c Sarah Khosravani 2014 Abstract The mechanical behaviour of a wet granular material is investigated through a microme- chanical analysis of force transport between interacting particles with a given packing and distribution of capillary liquid bridges. A single effective stress tensor, characterizing the ten- sorial contribution of the matric suction and encapsulating evolving liquid bridges, packing, interfaces, and water saturation, is derived micromechanically. The physical significance of the effective stress parameter (χ) as originally introduced in Bishop’s equation is examined and it turns out that Bishop’s equation is incomplete.
    [Show full text]
  • V- 3 31 Finite Element Implementation of State Variable-Based
    https://ntrs.nasa.gov/search.jsp?R=19910023209 2020-03-17T14:54:26+00:00Z View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by NASA Technical Reports Server ,,/v-_3 31 / NASA Contractor Report 187212 Finite Element Implementation of State Variable-Based Viscoplasticity Models (UA_A-CR-I_7212) FINITE ELEMENT N91-32523 IMPLEMENTATIL3N OF STATE VA&IABLF-_ASED VISCOPLA_TICITY MnOELS Final Report (Akron Univ.) 41 p CSCL 2OK G3/39 I. Iskovitz, T.Y.P. Chang, and A.F. Saleeb University of Akron Akron, Ohio September 1991 Prepared for Lewis Research Center Under Grant NAG3-901 National Aeronautics and Space Administration Finite Element Implementation of State Variable-Based Viscoplasticity Models I. Iskovitz, T.Y.P. Chang and A.F. Saleeb Department of Civil Engineering University of Akron Akron, Ohio 44325 Abstract Implementation of state variable-based viscoplasticity models is made in a general purpose finite element code for structural appli- cations of metals deformed at elevated temperature. Two consti- tutive models, i.e. Walker's and Robinson's models, are studied in conjunction with two implicit integration methods: the trapezoidal rule with Newton-Raphson iterations and an asymptotic integra- tion algorithm. A comparison is made between the two integration methods, and the latter method appears to be computationally more appealing in terms of numerical accuracy and CPU time. However, in order to make the asymptotic algorithm robust, it is necessary to include a self adaptive scheme with subincremental step control and error checking of the Jacobian matrix at the integration points. Three examples are given to illustrate the numerical aspects of the integration methods tested.
    [Show full text]
  • A Plasticity Model with Yield Surface Distortion for Non Proportional Loading
    page 1 A plasticity model with yield surface distortion for non proportional loading. Int. J. Plast., 17(5), pages :703–718 (2001) M.L.M. FRANÇOIS LMT Cachan (ENS de Cachan /CNRS URA 860/Université Paris 6) ABSTRACT In order to enhance the modeling of metallic materials behavior in non proportional loadings, a modification of the classical elastic-plastic models including distortion of the yield surface is proposed. The new yield criterion uses the same norm as in the classical von Mises based criteria, and a "distorted stress" Sd replacing the usual stress deviator S. The obtained yield surface is then “egg-shaped” similar to those experimentally observed and depends on only one new material parameter. The theory is built in such a way as to recover the classical one for proportional loading. An identification procedure is proposed to obtain the material parameters. Simulations and experiments are compared for a 2024 T4 aluminum alloy for both proportional and nonproportional tension-torsion loading paths. KEYWORDS constitutive behavior; elastic-plastic material. SHORTENED TITLE A model for the distortion of yield surfaces. ADDRESS Marc François, LMT, 61 avenue du Président Wilson, 94235 Cachan CEDEX, France. E-MAIL [email protected] INTRODUCTION Most of the phenomenological plasticity models for metals use von Mises based yield criteria. This implies that the yield surface is a hypersphere whose center is the back stress X in the 5-dimensional space of stress deviator tensors. For proportional loadings, the exact shape of the yield surface has no importance since all the deviators governing plasticity remain colinear tensors.
    [Show full text]
  • Identification of Rock Mass Properties in Elasto-Plasticity D
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Archive Ouverte en Sciences de l'Information et de la Communication Identification of rock mass properties in elasto-plasticity D. Deng, D. Nguyen Minh To cite this version: D. Deng, D. Nguyen Minh. Identification of rock mass properties in elasto-plasticity. Computers and Geotechnics, Elsevier, 2003, 30, pp.27-40. 10.1016/S0266-352X(02)00033-2. hal-00111371 HAL Id: hal-00111371 https://hal.archives-ouvertes.fr/hal-00111371 Submitted on 26 Sep 2019 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Identification of rock mass properties in elasto-plasticity Desheng Deng*, Duc Nguyen-Minh Laboratoire de Me´canique des Solides, E´cole Polytechnique, 91128 Palaiseau, France Abstract A simple and effective back analysis method has been proposed on the basis of a new cri-terion of identification, the minimization of error on the virtual work principle. This method works for both linear elastic and nonlinear elasto-plastic problems. The elasto-plastic rock mass properties for different criteria of plasticity can be well identified based on field measurements.
    [Show full text]
  • Crystal Plasticity Model with Back Stress Evolution. Wei Huang Louisiana State University and Agricultural & Mechanical College
    Louisiana State University LSU Digital Commons LSU Historical Dissertations and Theses Graduate School 1996 Crystal Plasticity Model With Back Stress Evolution. Wei Huang Louisiana State University and Agricultural & Mechanical College Follow this and additional works at: https://digitalcommons.lsu.edu/gradschool_disstheses Recommended Citation Huang, Wei, "Crystal Plasticity Model With Back Stress Evolution." (1996). LSU Historical Dissertations and Theses. 6154. https://digitalcommons.lsu.edu/gradschool_disstheses/6154 This Dissertation is brought to you for free and open access by the Graduate School at LSU Digital Commons. It has been accepted for inclusion in LSU Historical Dissertations and Theses by an authorized administrator of LSU Digital Commons. For more information, please contact [email protected]. INFORMATION TO USERS This manuscript has been reproduced from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type o f computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, charts) are reproduced by sectioning the original, beginning at the upper left-hand comer and continuing from left to right in equal sections with small overlaps.
    [Show full text]
  • Yield Surface Effects on Stablity and Failure
    XXIV ICTAM, 21-26 August 2016, Montreal, Canada YIELD SURFACE EFFECTS ON STABLITY AND FAILURE William M. Scherzinger * Solid Mechanics, Sandia National Laboratories, Albuquerque, New Mexico, USA Summary This work looks at the effects of yield surface description on the stability and failure of metal structures. We consider a number of different yield surface descriptions, both isotropic and anisotropic, that are available for modeling material behavior. These yield descriptions can all reproduce the same nominal stress-strain response for a limited set of load paths. However, for more general load paths the models can give significantly different results that can have a large effect on stability and failure. Results are shown for the internal pressurization of a cylinder and the differences in the predictions are quantified. INTRODUCTION Stability and failure in structures is an important field of study that provides solutions to many practical problems in solid mechanics. The boundary value problems require careful consideration of the geometry, the boundary conditions, and the material description. In metal structures the material is usually described with an elastic-plastic constitutive model, which can account for many effects observed in plastic deformation. Crystal plasticity models can account for many microstructural effects, but these models are often too fine scaled for production modeling and simulation. In this work we limit the study to rate independent continuum plasticity models. When using a continuum plasticity model, the hardening behavior is important and much attention is focused on it. Equally important, but less well understood, is the effect of the yield surface description. For classical continuum metal plasticity models that assume associated flow the yield surface provides a definition for the elastic region of stress space along with the direction of plastic flow when the stress state is on the yield surface.
    [Show full text]
  • Durability of Adhesively Bonded Joints in Aerospace Structures
    Durability of adhesive bonded joints in aerospace structures Washington State University Date: 11/08/2017 AMTAS Fall meeting Seattle, WA Durability of Bonded Aircraft Structure • Principal Investigators & Researchers • Lloyd Smith • Preetam Mohapatra, Yi Chen, Trevor Charest • FAA Technical Monitor • Ahmet Oztekin • Other FAA Personnel Involved • Larry Ilcewicz • Industry Participation • Boeing: Will Grace, Peter VanVoast, Kay Blohowiak 2 Durability of adhesive bonded joints in aerospace structures 1. Influence of yield criteria Yield criterion 2. Biaxial tests (Arcan) Plasticity 3. Cyclic tests Hardening rule 4. FEA model Nonlinearity in Bonded Joints 1. Creep, non-linear response Tension (closed form) 2. Ratcheting, experiment Time Dependence 3. Creep, model development Shear (FEA) 4. Ratcheting, model application Approach: Plasticity vPrimary Research Aim: Modeling of adhesive plasticity to describe nonlinear stress-strain response. ØSub Task: 1. Identify the influence of yield criterion and hardening rule on bonded joints (complete) 2. Characterizing hardening rule (complete Dec 2017) 3. Characterizing yield criterion (complete Dec 2017) 4. Numerically combining hardening rule and yield criterion (complete Feb 2017) Approach: Plasticity 1. Identify yield criterion and hardening rule of bonded joints • Input properties: Bulk and Thin film adhesive properties in Tension. 60 Characterization of Tensile properties 55 50 45 40 35 30 25 Bulk Tension: Tough Adheisve 20 Tensile Stress [MPa] Bulk Tension: Standard Adhesive 15 10 Thin film Tension: Tough Adhesive 5 Thin film Tension: Standard Adhesive 0 0 0.05 0.1 0.15 0.2 0.25 0.3 Axial Strain Approach: Plasticity 1. Identify the influence of yield criterion and hardening rule on bonded joints • Background: Investigate yield criterion and hardening rule of bonded joints.
    [Show full text]
  • Tuning the Performance of Metallic Auxetic Metamaterials by Using Buckling and Plasticity
    materials Article Tuning the Performance of Metallic Auxetic Metamaterials by Using Buckling and Plasticity Arash Ghaedizadeh 1, Jianhu Shen 1, Xin Ren 1,2 and Yi Min Xie 1,* Received: 28 November 2015; Accepted: 8 January 2016; Published: 18 January 2016 Academic Editor: Geminiano Mancusi 1 Centre for Innovative Structures and Materials, School of Engineering, RMIT University, GPO Box 2476, Melbourne 3001, Australia; [email protected] (A.G.); [email protected] (J.S.); [email protected] (X.R.) 2 Key Laboratory of Traffic Safety on Track, School of Traffic & Transportation Engineering, Central South University, Changsha 410075, China * Correspondence: [email protected]; Tel.: +61-3-9925-3655 Abstract: Metallic auxetic metamaterials are of great potential to be used in many applications because of their superior mechanical performance to elastomer-based auxetic materials. Due to the limited knowledge on this new type of materials under large plastic deformation, the implementation of such materials in practical applications remains elusive. In contrast to the elastomer-based metamaterials, metallic ones possess new features as a result of the nonlinear deformation of their metallic microstructures under large deformation. The loss of auxetic behavior in metallic metamaterials led us to carry out a numerical and experimental study to investigate the mechanism of the observed phenomenon. A general approach was proposed to tune the performance of auxetic metallic metamaterials undergoing large plastic deformation using buckling behavior and the plasticity of base material. Both experiments and finite element simulations were used to verify the effectiveness of the developed approach. By employing this approach, a 2D auxetic metamaterial was derived from a regular square lattice.
    [Show full text]
  • Constitutive Behavior of Granitic Rock at the Brittle-Ductile Transition
    CONSTITUTIVE BEHAVIOR OF GRANITIC ROCK AT THE BRITTLE-DUCTILE TRANSITION Josie Nevitt, David Pollard, and Jessica Warren Department of Geological and Environmental Sciences, Stanford University, Stanford, CA 94305 e-mail: [email protected] fracture to ductile flow with increasing depth in the Abstract earth’s crust. This so-called “brittle-ductile transition” has important implications for a variety of geological Although significant geological and geophysical and geophysical phenomena, including the mechanics processes, including earthquake nucleation and of earthquake rupture nucleation and propagation propagation, occur in the brittle-ductile transition, (Hobbs et al., 1986; Li and Rice, 1987; Scholz, 1988; earth scientists struggle to identify appropriate Tse and Rice, 1986). Despite the importance of constitutive laws for brittle-ductile deformation. This deformation within this lithospheric interval, significant paper investigates outcrops from the Bear Creek field gaps exist in our ability to accurately characterize area that record deformation at approximately 4-15 km brittle-ductile deformation. Such gaps are strongly depth and 400-500ºC. We focus on the Seven Gables related to the uncertainty in choosing the constitutive outcrop, which contains a ~10cm thick aplite dike that law(s) that govern(s) the rheology of the brittle-ductile is displaced ~45 cm through a contractional step transition. between two sub-parallel left-lateral faults. Stretching Rheology describes the response of a material to an and rotation of the aplite dike, in addition to local applied force and is defined through a set of foliation development in the granodiorite, provides an constitutive equations. In addition to continuum excellent measure of the strain within the step.
    [Show full text]
  • Plasticity, Localization, and Friction in Porous Materials
    RICE UNIVERSITY Plasticity, Localization, and Friction in Porous Materials by Daniel Christopher Badders A T h e s is S u b m i t t e d in P a r t ia l F u l f il l m e n t o f t h e R equirements f o r t h e D e g r e e Doctor of Philosophy A p p r o v e d , T h e s is C o m m i t t e e : Michael M. Carroll, Chairman Dean of the George R. Brown School of Engineering and Burton J. and Ann M. McMurtry Professor of Engineering Mary F. Wbjeeler and Applied Mathematics Yves C. Angel Associate Professor of Mechanical Engineering and Materials Science Houston, Texas May, 1995 Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. INFORMATION TO USERS This manuscript has been reproduced from the microfilm master. UMI films the text directly from the original or copy submitted Ihus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted Also, if unauthorized copyright material had to be removed anote will indicate the deletion.
    [Show full text]
  • On the Path-Dependence of the J-Integral Near a Stationary Crack in an Elastic-Plastic Material
    On the Path-Dependence of the J-integral Near a Stationary Crack in an Elastic-Plastic Material Dorinamaria Carka and Chad M. Landis∗ The University of Texas at Austin, Department of Aerospace Engineering and Engineering Mechanics, 210 East 24th Street, C0600, Austin, TX 78712-0235 Abstract The path-dependence of the J-integral is investigated numerically, via the finite element method, for a range of loadings, Poisson's ratios, and hardening exponents within the context of J2-flow plasticity. Small-scale yielding assumptions are employed using Dirichlet-to-Neumann map boundary conditions on a circular boundary that encloses the plastic zone. This construct allows for a dense finite element mesh within the plastic zone and accurate far-field boundary conditions. Details of the crack tip field that have been computed previously by others, including the existence of an elastic sector in Mode I loading, are confirmed. The somewhat unexpected result is that J for a contour approaching zero radius around the crack tip is approximately 18% lower than the far-field value for Mode I loading for Poisson’s ratios characteristic of metals. In contrast, practically no path-dependence is found for Mode II. The applications of T or S stresses, whether applied proportionally with the K-field or prior to K, have only a modest effect on the path-dependence. Keywords: elasto-plastic fracture mechanics, small scale yielding, path-dependence of the J-integral, finite element methods 1. Introduction The J-integral as introduced by Eshelby [1,2] and Rice [3] is perhaps the most useful quantity for the analysis of the mechanical fields near crack tips in both linear elastic and non-linear elastic materials.
    [Show full text]