Stress Integration for the Drucker-Prager Material Model Without Hardening Using the Incremental Plasticity Theory
|
|
- Elfrieda Pierce
- 6 years ago
- Views:
Transcription
1 Journal of the Serbian Society for Computational Mechanics / Vol. / No., 008 / pp UDC: 59.74:004.0 Stress Integration for the Drucker-rager Material Model Without Hardening Using the Incremental lasticity heory D. Rakić, M. Živković, R. Slavković, M. Kojić 4,5,6 Faculty of Mechanical ngineering, University of Kragujevac drakic@kg.ac.yu zile@kg.ac.yu radovan@kg.ac.yu 4 Harvard School of ublic Health, Boston, USA mkojic@hsph.harvard.edu 5 Research and Development Center for Bioengineering BIOIRC, Kragujevac, Serbia 6 Department of Biomedical ngineering, University of exas Medical Center at Houston, USA Abstract his paper presents the formulation of the stress integration procedure for the Drucker-rager material model with associative yielding condition by using the incremental plasticity theory. he idea of this method is to reach the solution by calculating the plastic matrix according to the method of incremental plasticity (used for elastic matrix corrections), and with use of the total strain increment. he computational procedure is implemented within the program package AK. Results of this procedure were compared with the solutions obtained by the governing parameter method (Kojic and Bathe 005), as well as with the results obtained by the program package NX-Nastran. Keywords: Drucker-rager model, FM, Geomechanics, Incremental plasticity. Introduction Stress integration represents calculation of stress change during an incremental step, corresponding to strain increments in the step. It is in essence the incremental integration of inelastic constitutive relationships to trace the history of material deformation. he stress integration is an important ingredient in the overall finite element inelastic analysis of structures. It is important that the integration algorithm accurately reproduces the material behavior since the mechanical response of the entire structure is directly dependent on this accuracy. he algorithm should be also computationally efficient because the stress integration is performed at all integration points. For general applications, this computational procedure should be robust, providing reliable results under all possible loading (straining) conditions (see e.g. discussion of these issues in e.g. Kojic and Bathe 005). In this paper we present a formulation of the computational algorithm for the Drucker- rager material model without hardening using an incremental plasticity approach. he results using this approach were compared with the results obtained using the implicit integrations
2 Journal of the Serbian Society for Computational Mechanics / Vol. / No., scheme termed as the governing parameter method (Kojic, 996), as well as with results obtained by using other software packages. In the next section we present formulation the Drucker-rager material model, followed by the derivation of the elastic-plastic constitutive matrix in general associated plasticity. hen, the general relations are implemented on of the Drucker-rager material model.. Formulation of the Drucker-rager model his material model is defined by the yield surface called the Drucker-rager yield surface. In the space of main stresses, this surface represents a cone, whose axis matches the space diagonal of the coordinate system of principal stresses σ, σ, σ (Figure ). his material model is a special case of the Drucker-rager material model with a cap (Kojic and Bathe 005). Fig.. Drucker-rager failure surface. he model can be represented in the I J coordinate system, as shown in Fig., where Fig.. Drucker-rager model in space I J. is the first stress invariant, and I σ m σδ ij ij σ ii ()
3 8 D. Rakić et al. : Stress Integration for the Drucker-rager Material Model J SS () ij ij is the second invariant of the deviatoric stress. Deviatoric stress S ij is defined as: Sij σ ij σmδij () he elastic domain is bounded by the Drucker-rager surface and by the maximum tension line with the tension cutoff stress. his representation of the model is suitable for development of the computational procedure. If the stress point of a trial solution is above the Drucker-rager line, the stresses at the end of step have to be reduced to satisfy the following yield condition: f I J k (4) tδ t tδ t tδt D α 0 where α and k are material constants; and the upper left index t Δ t denotes the end of step. If the stress exceeds the allowed tension stress value, the stress reduction is performed by equalizing all three components of normal stress to one third of the maximum tension stress, i. e.: σ σ σ (5) while the shear stress components are set to zero: σ ij 0 i j (6). lastic-plastic matrix For an elastic-plastic material, the stress increment { dσ } can be expressed by the total strain increment { de } as { dσ} C { de} (7) where C is the elastic-plastic constitutive matrix defined below. In the case of small strain material deformation, the total strain increment can be divided into elastic and plastic part: from which follows { } { de de } { de } (8) { de } { de} { de } (9) he stress and elastic strain increments are related by the elastic constitutive law, { } dσ C { de } (0) where C is elastic constitutive matrix. Substituting (9) into (0), we obtain
4 Journal of the Serbian Society for Computational Mechanics / Vol. / No., { dσ} C ( { de} { de }) () We next include into this derivation a yield function F, which in general has a form: ( ) ( ) A differential of the yield the function can be expressed as F F σ F σ ij () F F F F F F df dσ dσ dσ dσ dσ dσ σ σ σ σ σ σ xx yy zz xy yz zx xx yy zz xy yz zx () or in matrix form, df { dσ} In the incremental theory of plasticity we have that the yield function during plastic deformation must be equal to zero (when F 0 we have elastic unloading), hence F 0 and consequently df 0 is satisfied over the load step. he last condition can be written in the following form: (4) df { dσ} 0 (5) where is whereas { dσ } is F F F F F F σxx σyy σzz σxy σyz σzx (6) { dσ} dσxx dσ yy dσzz dσxy dσ yz dσzx (7) Further, we assume a non-associative plasticity, therefore the increment of plastic strain is { de } G dλ where: G ( σ ) is the plastic potential function, and dλ is a positive scalar. (8) Substituting the plastic strain increment (8) into (), we obtain the stress increment: { σ} { } G d C de dλ C σ Now, we substitute this stress increment into (5): (9) or F df C de d C G σ { } λ 0 (0)
5 84 D. Rakić et al. : Stress Integration for the Drucker-rager Material Model hen, the parameter dλ is df C de d C F F G σ σ { } λ 0 () { de} C σ dλ F G C Finally, substituting dλ from () into (9), obtain the relationship between { dσ } and the total strain increment { de } : () { dσ} C { de} G C C de σ σ G C By comparing this and equation (7) it follows: { } () G C C σ σ C C G C his is a general expression for the elastic-plastic constitutive matrix. Also, by substituting () into (8) we can calculate the plastic strain increment{ de } : (4) { de } G F C { de} σ σ G C In the case of associative plasticity the yield function F is at the same time the plastic potential function. In our implementation of the above general expressions to the Drucker- rager model we will adopt the associative plasticity assumption. (5) 4. Implementation of the general associative plasticity relations to the Drucker-rager material model We further use the Drucker-rager yielding surface defined as: F α I J k (6) and assume the associative plasticity, hence we have that the plastic potential is defined by the same function:
6 Journal of the Serbian Society for Computational Mechanics / Vol. / No., G α I J k (7) dil he first invariant I is defined by equation (), while the second stress invariant J is defined by (): J ( σ σ ) ( σ σ ) ( σ σ ) σ σ σ where the one-index notation is used (as in the subsequent text). Also, we can express J as J I I (9) where the second stress invariant I is: (8) I σ σ σ σ σ σ σ σ σ (0) Most frequently used material constants of geological materials (which are experimentally evaluated) are the internal friction angle φ and the cohesion c. hey are the characteristic of the Mohr-Coulomb material model. here are direct relations between these parameters and parameters k and α of the Ducker-rager model: α sinφ ( φ ) ( sinφ ) sin k 6c cosφ () hese relations are obtained by coinciding the Drucker-rager cone with the outer apexes of the Mohr-Coulomb hexagon locus (compressive cone). he derivatives of the Drucker-rager yield function with respect to stresses can be calculated by using the chain rule: J F I F () I J σ Note that this expression is applicable to all models where the yield function is expressed in terms the stress invariants I and J (which is the case for a large number of plasticity models). Also we have the following relations from (6): F α () I and F J J he derivatives of the stress invariants with respect to stresses are: (4) and I [ 0 0 0] (5)
7 86 D. Rakić et al. : Stress Integration for the Drucker-rager Material Model J ( ) ( ) D σ σ σ σ σ σ ( σ σ σ) σ4 σ5 σ6 (6) σ Now, substituting () to (6) into () we obtain: σ σ σ α 6 J σ σ σ α 6 J σ σ σ α 6 J σ 4 J σ 5 J σ 6 J (7) he above relations for the stress integration over a load (time) step are summarized in able. he left upper indices t and t Δ t denote start and end of a time step. tδ Known quantities: { t e}, { t e }, { t t p σ }, { e } A. rial (elastic) results: tδt t { dσ} [ Ce]{ de } [ Ce] ({ e} { e} ) { t Δ t σ} { t σ} { dσ} I σ σ σ J ( σ σ ) ( σ σ ) ( σ σ ) σ σ σ F α I J k B. Yielding condition check: IF ( F < 0 ) trial result correct (go to ) IF ( F 0 ) elastic-plastic result (CONINU) J F I F I J σ
8 Journal of the Serbian Society for Computational Mechanics / Vol. / No., Associative yielding assumption: G [ Ce ]{ de} σ dλ F G [ Ce ] C. Stress correction (solution by bisection - bisection loop): { de } D. IF ( ABS ( F ) G dλ { de } { de} { de } { dσ } [ C ]{ e de } { t Δ t σ} { t σ} { dσ} Calculation of new invariants: I σ σ σ J ( σ σ ) ( σ σ ) ( σ σ ) σ σ σ F α I J k tδt t { e } { e } { de } tδ. ND: { t tδt p σ},{ e } OL ) go back to C with new dλ from bisection able. Algorithm for incremental stress integration 5. Verification of the computational procedure Verification of the proposed computational procedure was done by solving elastic-plastic deformation of a sand box specimen, as specified in the literature (Kojic and Bathe 005). he loading is modeled by prescribed displacement on the top, while only vertical sliding is allowed on the sides. he bottom model side is completely restricted. A half of the specimen is modeled with the use of the appropriate symmetry boundary conditions (Figure ). ight-node D finite elements were used. Displacement is applied at the nodes on the model upper side.
9 88 D. Rakić et al. : Stress Integration for the Drucker-rager Material Model [] step [] displacement u [] 0 [4] 0.0 [5] 0 [6] 0.0 [7] 0 [8] 0.05 Material data: Young s modulus, 00 oisson s ratio, ν 0. Yield function parameter, α 0.05 Yield function parameter, k 0. ension Cutoff limit, 0.0 Dimension, D.0 Dimension, H 0.5 Fig.. Sand specimen geometry, material data and load. he problem was solved using the program package AK (Kojic et al. 998), by employing the described method and the governing parameter method. Also the program NX NASRAN was used. Figures 4 and 5 show the results. e xx Nas tran AK σ xx Fig. 4. he calculated axial stress as a function of the axial strain.
10 Journal of the Serbian Society for Computational Mechanics / Vol. / No., Failure surface AK sqrt(j ) -5 I Fig. 5. he second invariant deviatoric stress as function of first stress invariant. It can be seen from the presented results that all three solutions are practically the same. 6. Conclusions he presented methodology gives very accurate results when compared with other methods, such as the governing parameter method. he advantage of the presented computational procedure is that it is formulated in general form so that it can be applied to various yield functions. Also, this procedure can be implemented in an explicit computationally efficient incremental scheme, where yield condition check D (able ) is not performed, but then very small load steps are necessary. References Bathe K. J. (98), Finite lement rocedures in ngineering Analysis, rentice-hall, Inc., nglewood Cliffs, New Jersey. Geo-Slope Office (00), Sigma W-for finite element stress and deformation analysis V5, User s guide. Kojic M. (996), he governing parameter method for implicit integration of viscoplastic constitutive relations for isotropic and orthotropic metals, Computational Mechanics, Vol. 9, No., pp Kojic M., Bathe K. J. (005), Inelastic Analysis of Solids and Structures, Springer Berlin Heidelberg New York. Kojic M., Slavkovic R., Zivkovic M., Grujovic N. (998), AK-finite element program for linear and nonlinear structural analysis and heat transfer, Faculty for Mechanical ngineering, Kragujevac, University of Kragujevac. Woods R., Rahim A. (007), CRIS Geotechnical Finite lement Analysis Software, echnical Reference Manual.
Finite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationMechanical Properties of Materials
Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of
More informationGEO E1050 Finite Element Method Mohr-Coulomb and other constitutive models. Wojciech Sołowski
GEO E050 Finite Element Method Mohr-Coulomb and other constitutive models Wojciech Sołowski To learn today. Reminder elasticity 2. Elastic perfectly plastic theory: concept 3. Specific elastic-perfectly
More informationDEFORMATION PATTERN IN ELASTIC CRUST
DEFORMATION PATTERN IN ELASTIC CRUST Stress and force in 2D Strain : normal and shear Elastic medium equations Vertical fault in elastic medium => arctangent General elastic dislocation (Okada s formulas)
More informationConstitutive models: Incremental plasticity Drücker s postulate
Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager
More informationPLAXIS 3D TUNNEL. Material Models Manual version 2
PLAXIS 3D TUNNEL Material Models Manual version 2 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction...1-1 1.1 On the use of different models...1-1 1.2 Limitations...1-2 2 Preliminaries on material modelling...2-1
More informationDEFORMATION PATTERN IN ELASTIC CRUST
DEFORMATION PATTERN IN ELASTIC CRUST Stress and force in 2D Strain : normal and shear Elastic medium equations Vertical fault in elastic medium => arctangent General elastic dislocation (Okada s formulas)
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More information1. Background. is usually significantly lower than it is in uniaxial tension
NOTES ON QUANTIFYING MODES OF A SECOND- ORDER TENSOR. The mechanical behavior of rocks and rock-like materials (concrete, ceramics, etc.) strongly depends on the loading mode, defined by the values and
More informationSOIL MODELS: SAFETY FACTORS AND SETTLEMENTS
PERIODICA POLYTECHNICA SER. CIV. ENG. VOL. 48, NO. 1 2, PP. 53 63 (2004) SOIL MODELS: SAFETY FACTORS AND SETTLEMENTS Gabriella VARGA and Zoltán CZAP Geotechnical Department Budapest University of Technology
More informationBone Tissue Mechanics
Bone Tissue Mechanics João Folgado Paulo R. Fernandes Instituto Superior Técnico, 2016 PART 1 and 2 Introduction The objective of this course is to study basic concepts on hard tissue mechanics. Hard tissue
More informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationVerification of the Hyperbolic Soil Model by Triaxial Test Simulations
1 Introduction Verification of the Hyperbolic Soil Model by Triaxial Test Simulations This example simulates a series of triaxial tests that can be used to verify that the Hyperbolic constitutive model
More informationClassical fracture and failure hypotheses
: Chapter 2 Classical fracture and failure hypotheses In this chapter, a brief outline on classical fracture and failure hypotheses for materials under static loading will be given. The word classical
More information1 Introduction. Abstract
Abstract This paper presents a three-dimensional numerical model for analysing via finite element method (FEM) the mechanized tunneling in urban areas. The numerical model is meant to represent the typical
More informationEffect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities
Effect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities Hany El Naggar, Ph.D., P. Eng. and M. Hesham El Naggar, Ph.D., P. Eng. Department of Civil Engineering
More informationELASTICITY (MDM 10203)
ELASTICITY () Lecture Module 3: Fundamental Stress and Strain University Tun Hussein Onn Malaysia Normal Stress inconstant stress distribution σ= dp da P = da A dimensional Area of σ and A σ A 3 dimensional
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationELASTICITY (MDM 10203)
LASTICITY (MDM 10203) Lecture Module 5: 3D Constitutive Relations Dr. Waluyo Adi Siswanto University Tun Hussein Onn Malaysia Generalised Hooke's Law In one dimensional system: = (basic Hooke's law) Considering
More information2 CONSTITUTIVE MODELS: THEORY AND IMPLEMENTATION
CONSTITUTIVE MODELS: THEORY AND IMPLEMENTATION 2-1 2 CONSTITUTIVE MODELS: THEORY AND IMPLEMENTATION 2.1 Introduction There are twelve basic constitutive models provided in, arranged into null, elastic
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationExample-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium
Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic
More informationNUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS
IGC 009, Guntur, INDIA NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS Mohammed Younus Ahmed Graduate Student, Earthquake Engineering Research Center, IIIT Hyderabad, Gachibowli, Hyderabad 3, India.
More informationCONSTITUTIVE RELATIONS FOR LINEAR ELASTIC SOLIDS
Chapter 9 CONSTITUTIV RLATIONS FOR LINAR LASTIC SOLIDS Figure 9.1: Hooke memorial window, St. Helen s, Bishopsgate, City of London 211 212 CHAPTR 9. CONSTITUTIV RLATIONS FOR LINAR LASTIC SOLIDS 9.1 Mechanical
More informationPLAXIS. Material Models Manual
PLAXIS Material Models Manual 2015 Build 7519 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 7 1.1 On the use of different models 7 1.2 Limitations 9 2 Preliminaries on material modelling 13 2.1 General
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationCard Variable MID RO E PR ECC QH0 FT FC. Type A8 F F F F F F F. Default none none none 0.2 AUTO 0.3 none none
Note: This is an extended description of MAT_273 input provided by Peter Grassl It contains additional guidance on the choice of input parameters beyond the description in the official LS-DYNA manual Last
More information1. Vladimir MILOVANOVIĆ, 2. Dragan RAKIĆ, 3. Miroslav ŽIVKOVIĆ, 4. Snežana VULOVIĆ, 5. Miroslav MILUTINOVIĆ
1. Vladimir MILOVANOVIĆ, 2. Dragan RAKIĆ, 3. Miroslav ŽIVKOVIĆ, 4. Snežana VULOVIĆ, 5. Miroslav MILUTINOVIĆ THERMO-MECHANIC ANALYSIS OF CEMENT TRANSPORT WAGON - IDENTIFICATION OF THE CAUSE OF CRACKS 1-3.
More informationDeformable Materials 2 Adrien Treuille. source: Müller, Stam, James, Thürey. Real-Time Physics Class Notes.
Deformable Materials 2 Adrien Treuille source: Müller, Stam, James, Thürey. Real-Time Physics Class Notes. Goal Overview Strain (Recap) Stress From Strain to Stress Discretization Simulation Overview Strain
More informationFEM for elastic-plastic problems
FEM for elastic-plastic problems Jerzy Pamin e-mail: JPamin@L5.pk.edu.pl With thanks to: P. Mika, A. Winnicki, A. Wosatko TNO DIANA http://www.tnodiana.com FEAP http://www.ce.berkeley.edu/feap Lecture
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationTIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL
TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL Qassun S. Mohammed Shafiqu and Maarib M. Ahmed Al-Sammaraey Department of Civil Engineering, Nahrain University, Iraq
More informationNonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface
Transaction A: Civil Engineering Vol. 16, No. 6, pp. 512{519 c Sharif University of Technology, December 2009 Research Note Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield
More informationA circular tunnel in a Mohr-Coulomb medium with an overlying fault
MAP3D VERIFICATION EXAMPLE 9 A circular tunnel in a Mohr-Coulomb medium with an overlying fault 1 Description This example involves calculating the stresses and displacements on a fault overlying a 5 m
More informationModified Cam-clay triaxial test simulations
1 Introduction Modified Cam-clay triaxial test simulations This example simulates a series of triaxial tests which can be used to verify that Modified Cam-Clay constitutive model is functioning properly.
More informationPrediction of torsion shear tests based on results from triaxial compression tests
Prediction of torsion shear tests based on results from triaxial compression tests P.L. Smith 1 and N. Jones *2 1 Catholic University of America, Washington, USA 2 Geo, Lyngby, Denmark * Corresponding
More informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More informationCONTENTS. Lecture 1 Introduction. Lecture 2 Physical Testing. Lecture 3 Constitutive Models
CONTENTS Lecture 1 Introduction Introduction.......................................... L1.2 Classical and Modern Design Approaches................... L1.3 Some Cases for Numerical (Finite Element) Analysis..........
More informationA simple elastoplastic model for soils and soft rocks
A simple elastoplastic model for soils and soft rocks A SIMPLE ELASTO-PLASTIC MODEL FOR SOILS AND SOFT ROCKS by Roberto Nova Milan University of Technology 1. MODEL HISTORY The model is the result of the
More informationANSYS Mechanical Basic Structural Nonlinearities
Lecture 4 Rate Independent Plasticity ANSYS Mechanical Basic Structural Nonlinearities 1 Chapter Overview The following will be covered in this Chapter: A. Background Elasticity/Plasticity B. Yield Criteria
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Uniaxial Model: Strain-Driven Format of Elastoplasticity
MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #3: Elastoplastic Concrete Models Uniaxial Model: Strain-Driven Format of Elastoplasticity Triaxial
More information2D Liquefaction Analysis for Bridge Abutment
D Liquefaction Analysis for Bridge Abutment Tutorial by Angel Francisco Martinez Integrated Solver Optimized for the next generation 64-bit platform Finite Element Solutions for Geotechnical Engineering
More information14- Hardening Soil Model with Small Strain Stiffness - PLAXIS
14- Hardening Soil Model with Small Strain Stiffness - PLAXIS This model is the Hardening Soil Model with Small Strain Stiffness as presented in PLAXIS. The model is developed using the user-defined material
More informationA return mapping algorithm for unified strength theory model
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2015; 104:749 766 Published online 10 June 2015 in Wiley Online Library (wileyonlinelibrary.com)..4937 A return mapping
More informationA Simple and Accurate Elastoplastic Model Dependent on the Third Invariant and Applied to a Wide Range of Stress Triaxiality
A Simple and Accurate Elastoplastic Model Dependent on the Third Invariant and Applied to a Wide Range of Stress Triaxiality Lucival Malcher Department of Mechanical Engineering Faculty of Tecnology, University
More informationUnified Constitutive Model for Engineering- Pavement Materials and Computer Applications. University of Illinois 12 February 2009
Unified Constitutive Model for Engineering- Pavement Materials and Computer Applications Chandrakant S. Desai Kent Distinguished i Lecture University of Illinois 12 February 2009 Participation in Pavements.
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More informationCombined Isotropic-Kinematic Hardening Laws with Anisotropic Back-stress Evolution for Orthotropic Fiber-Reinforced Composites
Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution for Orthotropic Fiber- Reinforced Composites Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution
More informationAuthor(s) Okajima, Kenji; Tanaka, Tadatsugu; Symposium on Backwards Problem in G.
Title Backwards Analysis for Retaining Wa based upon ateral Wall Displacemen Author(s) Okajima, Kenji; Tanaka, Tadatsugu; Proceeding of TC302 Symposium Osaka Citation Symposium on Backwards Problem in
More informationNumerical Modeling of Direct Shear Tests on Sandy Clay
Numerical Modeling of Direct Shear Tests on Sandy Clay R. Ziaie Moayed, S. Tamassoki, and E. Izadi Abstract Investigation of sandy clay behavior is important since urban development demands mean that sandy
More informationWellbore stability analysis in porous carbonate rocks using cap models
Wellbore stability analysis in porous carbonate rocks using cap models L. C. Coelho 1, A. C. Soares 2, N. F. F. Ebecken 1, J. L. D. Alves 1 & L. Landau 1 1 COPPE/Federal University of Rio de Janeiro, Brazil
More informationEffect of Plasticity on Residual Stresses Obtained by the Incremental Hole-drilling Method with 3D FEM Modelling
Effect of Plasticity on Residual Stresses Obtained by the Incremental Hole-drilling Method with 3D FEM Modelling Evy Van Puymbroeck 1, a *, Wim Nagy 1,b and Hans De Backer 1,c 1 Ghent University, Department
More informationEnhancing Prediction Accuracy In Sift Theory
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS Enhancing Prediction Accuracy In Sift Theory J. Wang 1 *, W. K. Chiu 1 Defence Science and Technology Organisation, Fishermans Bend, Australia, Department
More informationABSTRACT INTRODUCTION
Optimization of soil anchorages M. Teschner, C. Mattheck Kernforschungszentrum Karlsruhe GmbH, Institut fur Materialforschung II, W-7500 Karlsruhe 1, Postfach 3640, Germany ABSTRACT A new optimization
More informationNumerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour
13 th International Brick and Block Masonry Conference Amsterdam, July 4-7, 2004 Numerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour H. R.
More informationHypoplastic Cam-clay model
Hypoplastic model David Mašín Charles University in Prague corresponence address: David Mašíın Charles University in Prague Faculty of Science Albertov 6 12843 Prague 2, Czech Republic E-mail: masin@natur.cuni.cz
More informationUseful Formulae ( )
Appendix A Useful Formulae (985-989-993-) 34 Jeremić et al. A.. CHAPTER SUMMARY AND HIGHLIGHTS page: 35 of 536 A. Chapter Summary and Highlights A. Stress and Strain This section reviews small deformation
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationCOMPARISON OF CONSTITUTIVE SOIL MODELS FOR THE SIMULATION OF DYNAMIC ROLLER COMPACTION
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2012) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 10-14, 2012 COMPARISON OF CONSTITUTIVE SOIL MODELS
More informationApplication of Three Dimensional Failure Criteria on High-Porosity Chalk
, 5-6 May 00, Trondheim, Norway Nordic Energy Research Programme Norwegian U. of Science and Technology Application of Three Dimensional Failure Criteria on High-Porosity Chalk Roar Egil Flatebø and Rasmus
More informationAdvanced model for soft soils. Modified Cam-Clay (MCC)
Advanced model for soft soils. Modified Cam-Clay (MCC) c ZACE Services Ltd August 2011 1 / 62 2 / 62 MCC: Yield surface F (σ,p c ) = q 2 + M 2 c r 2 (θ) p (p p c ) = 0 Compression meridian Θ = +π/6 -σ
More informationDYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION
October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,
More information20. Rheology & Linear Elasticity
I Main Topics A Rheology: Macroscopic deformation behavior B Linear elasticity for homogeneous isotropic materials 10/29/18 GG303 1 Viscous (fluid) Behavior http://manoa.hawaii.edu/graduate/content/slide-lava
More informationEffect of the intermediate principal stress on fault strike and dip - theoretical analysis and experimental verification
Effect of the intermediate principal stress on fault strike and dip - theoretical analysis and experimental verification B. Haimson University of Wisconsin, USA J. Rudnicki Northwestern University, USA
More informationCONSOLIDATION BEHAVIOR OF PILES UNDER PURE LATERAL LOADINGS
VOL., NO., DECEMBER 8 ISSN 89-8 -8 Asian Research Publishing Network (ARPN). All rights reserved. CONSOLIDATION BEAVIOR OF PILES UNDER PURE LATERAL LOADINGS Qassun S. Mohammed Shafiqu Department of Civil
More information1.8 Unconfined Compression Test
1-49 1.8 Unconfined Compression Test - It gives a quick and simple measurement of the undrained strength of cohesive, undisturbed soil specimens. 1) Testing method i) Trimming a sample. Length-diameter
More informationChapter 6: Plastic Theory
OHP Mechanical Properties of Materials Chapter 6: Plastic Theory Prof. Wenjea J. Tseng 曾文甲 Department of Materials Engineering National Chung Hsing University wenjea@dragon.nchu.edu.tw Reference: W. F.
More informationThe Finite Element Method II
[ 1 The Finite Element Method II Non-Linear finite element Use of Constitutive Relations Xinghong LIU Phd student 02.11.2007 [ 2 Finite element equilibrium equations: kinematic variables Displacement Strain-displacement
More informationRock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth
Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of
More informationLocalization in Undrained Deformation
Localization in Undrained Deformation J. W. Rudnicki Dept. of Civil and Env. Engn. and Dept. of Mech. Engn. Northwestern University Evanston, IL 621-319 John.Rudnicki@gmail.com Fourth Biot Conference on
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationTRESS - STRAIN RELATIONS
TRESS - STRAIN RELATIONS Stress Strain Relations: Hook's law, states that within the elastic limits the stress is proportional to t is impossible to describe the entire stress strain curve with simple
More informationNumerical model comparison on deformation behavior of a TSF embankment subjected to earthquake loading
Numerical model comparison on deformation behavior of a TSF embankment subjected to earthquake loading Jorge Castillo, Yong-Beom Lee Ausenco, USA Aurelian C. Trandafir Fugro GeoConsulting Inc., USA ABSTRACT
More informationEffect Of The In-Situ Stress Field On Casing Failure *
Effect Of The In-Situ Stress Field On Casing Failure * Tang Bo Southwest Petroleum Institute, People's Republic of China Lian Zhanghua Southwest Petroleum Institute, People's Republic of China Abstract
More informationNumerical Assessment of the Influence of End Conditions on. Constitutive Behavior of Geomaterials
Numerical Assessment of the Influence of End Conditions on Constitutive Behavior of Geomaterials Boris Jeremić 1 and Zhaohui Yang 2 and Stein Sture 3 ABSTRACT In this paper we investigate the behavior
More informationDUCTILITY BEHAVIOR OF A STEEL PLATE SHEAR WALL BY EXPLICIT DYNAMIC ANALYZING
The 4 th World Conference on arthquake ngineering October -7, 008, Beijing, China ABSTRACT : DCTILITY BHAVIOR OF A STL PLAT SHAR WALL BY XPLICIT DYNAMIC ANALYZING P. Memarzadeh Faculty of Civil ngineering,
More informationFinite Element Solutions for Geotechnical Engineering
Release Notes Release Date: June, 2017 Product Ver.: GTSNX 2017(v1.1) Integrated Solver Otimized for the next generation 64-bit latform Finite Element Solutions for Geotechnical Engineering 1. Analysis
More information8.1. What is meant by the shear strength of soils? Solution 8.1 Shear strength of a soil is its internal resistance to shearing stresses.
8.1. What is meant by the shear strength of soils? Solution 8.1 Shear strength of a soil is its internal resistance to shearing stresses. 8.2. Some soils show a peak shear strength. Why and what type(s)
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Standard Solids and Fracture Fluids: Mechanical, Chemical Effects Effective Stress Dilatancy Hardening and Stability Mead, 1925
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Thermodynamics Derivation Hooke s Law: Anisotropic Elasticity
More informationModelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases
Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Georg Haasemann Conrad Kloß 1 AIMCAL Conference 2016 MOTIVATION Wrinkles in web handling system Loss of
More informationFinite Element Solutions for Geotechnical Engineering
Release Notes Release Date: July, 2015 Product Ver.: GTSNX 2015 (v2.1) Integrated Solver Optimized for the next generation 64-bit platform Finite Element Solutions for Geotechnical Engineering Enhancements
More informationA Constitutive Model for DYNEEMA UD composites
A Constitutive Model for DYNEEMA UD composites L Iannucci 1, D J Pope 2, M Dalzell 2 1 Imperial College, Department of Aeronautics London, SW7 2AZ l.iannucci@imperial.ac.uk 2 Dstl, Porton Down, Salisbury,
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationLimit analysis of brick masonry shear walls with openings under later loads by rigid block modeling
Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling F. Portioli, L. Cascini, R. Landolfo University of Naples Federico II, Italy P. Foraboschi IUAV University,
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More information3.2 Hooke s law anisotropic elasticity Robert Hooke ( ) Most general relationship
3.2 Hooke s law anisotropic elasticity Robert Hooke (1635-1703) Most general relationship σ = C ε + C ε + C ε + C γ + C γ + C γ 11 12 yy 13 zz 14 xy 15 xz 16 yz σ = C ε + C ε + C ε + C γ + C γ + C γ yy
More informationLateral Earth Pressure
1 of 11 6/2/2012 4:28 AM Lateral Earth Pressure The magnitude of lateral earth pressure depends on: 1. Shear strength characteristics of soil 2. Lateral strain condition 3. Pore water pressure 4. State
More informationPLAXIS LIQUEFACTION MODEL UBC3D-PLM
PLAXIS LIQUEFACTION MODEL UBC3D-PLM Alexandros Petalas Assistant Researcher, PLAXIS B.V Vahid Galavi Researcher, PLAXIS B.V May 30, 2012 Contents 1 Key Features of UBC3D 2 1.1 Yield Surfaces...........................
More informationConstitutive models. Part 2 Elastoplastic
Contitutive model art latoplatic latoplatic material model latoplatic material are aumed to behave elatically up to a certain tre limit after which combined elatic and platic behaviour occur. laticity
More informationMATERIAL PROPERTIES. Material Properties Must Be Evaluated By Laboratory or Field Tests 1.1 INTRODUCTION 1.2 ANISOTROPIC MATERIALS
. MARIAL PROPRIS Material Properties Must Be valuated By Laboratory or Field ests. INRODUCION he fundamental equations of structural mechanics can be placed in three categories[]. First, the stress-strain
More informationA Two Dimensional Finite Element Analysis of a Plane Tillage Tool in Soil Using a Non-linear Elasto-Plastic Model
American-Eurasian J. Agric. & Environ. Sci., 3 (3): 498-505, 2008 ISSN 88-6769 IDOSI Publications, 2008 A Two Dimensional Finite Element Analysis of a Plane Tillage Tool in Soil Using a Non-linear Elasto-Plastic
More informationThe Influence of Strain Amplitude, Temperature and Frequency on Complex Shear Moduli of Polymer Materials under Kinematic Harmonic Loading
Mechanics and Mechanical Engineering Vol. 21, No. 1 (2017) 157 170 c Lodz University of Technology The Influence of Strain Amplitude, Temperature and Frequency on Complex Shear Moduli of Polymer Materials
More informationNUMERICAL SIMULATIONS OF CORNERS IN RC FRAMES USING STRUT-AND-TIE METHOD AND CDP MODEL
Numerical simulations of corners in RC frames using Strut-and-Tie Method and CDP model XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J.
More informationFLAC3D analysis on soil moving through piles
University of Wollongong Research Online Faculty of Engineering - Papers (Archive) Faculty of Engineering and Information Sciences 211 FLAC3D analysis on soil moving through piles E H. Ghee Griffith University
More informationConstitutive modelling of fabric anisotropy in sand
Geomechanics from Micro to Macro Soga et al. (Eds) 2015 Taylor & Francis Group, London, ISBN 978-1-138-02707-7 Constitutive modelling of fabric anisotropy in sand Z.W. Gao School of Engineering, University
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More information