Computational models of diamond anvil cell compression

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Computational models of diamond anvil cell compression"

Transcription

1 UDC Computational models of diamond anvil cell compression A. I. Kondrat yev Independent Researcher, 5944 St. Alban Road, Pensacola, Florida 32503, USA Abstract. Diamond anvil cells (DAC) are extensively utilized in the study of material properties at extreme conditions of pressures and temperatures. The pressures in excess of 400 GPa and temperatures in excess of 4000 K have been reported in DAC devices. Most studies on the optimization of diamond geometry and gasket materials used in diamond anvil cell devices have been carried out by trial and error using experimental high pressure data. This article addresses design and optimization issues in DAC using finite element modeling (FEM) and computational analysis. The computational approach is centered around the recent advances in the growth of isotopically enriched carbon-13 layers on diamond anvil and their use as pressure sensor in diamond anvil cell devices. In particular, the carbon-13 diamond layer of 6 microns in thickness grown on top of existing anvil has been demonstrated to serve as a universal pressure sensor to 156 GPa. If a thin enough pressure sensor can be fabricated then the calibration of this sensor is independent of the stress or strain distribution in DAC. Some questions about the geometry and other properties of DAC remain open as well the theory of a universal pressure sensor. Keywords: Diamond anvil cell device, Finite-element modeling, Diamond-coated rhenium gasket, Diamond anvil failure conditions. 1. Introduction The object of our research is DAC device main component of which is diamond anvil. In DAC configuration a sample chamber is placed between the polished culets of diamond anvils. The metal gasket is preindented between the two anvils before drilling a hole that serves as a pressure chamber. All components are thoroughly assembled in a specially sealed cylinder. Then using tools we can apply pressure by turning special nuts in the assembly. Object of our investigation is the DAC device. We study separately DAC device components: diamond anvils, gasket, and sample material. In those components we study radial and axial stresses as well as pressure distribution with respect to geometry of anvils, sample and gasket material properties, and radial and axial coordinates. We have complete description of geometry and material properties of all DAC components. Because of DAC special design and diamond transparency we have X-ray diffraction experiments data for sample materials obtained within DAC compression. In some cases, diamonds were compressed to failure to test the upper limit of the pressures that can be created in diamond anvil cell

2 device. In order to use results of diffraction experiments we have to use equations of state for diamond, sample material, and a gasket. As a result, we obtain dependencies of pressure distribution vs. radial distance. The pressure distribution in DAC can be experimentally measured only at the diamond/gasket interface and diamond/sample interface or at z = 0. Obtaining experimental data of pressure vs. axial distance is more complicated problem. At this stage of DAC device technology we are not able to measure axial and radial stresses. In order to calculate this information we have to use finite element analysis and modeling (FEAM) and computer simulation as well as analytical modeling and solution to main problems. Further development of finite element modeling (FEM) give us the full computer simulation model of DAC compression process. We develop and extend diamond anvil capabilities in order to use it as a high pressure sensor. When a thin film of 13 C layer was grown on a top of regular diamond anvil this thin film maybe used as a high pressure sensor using Raman spectroscopy. Our last study is by using analytical description and analysis of stresses and strains in an anvil and in a layer and also by using FEAM and computer simulation to analyze pressure distribution in a layer by describing the behavior of optical transverse phonon mode in an anvil and in a layer. Our working approach is based on use FEM and variety of analytical models with conjunction of experimental data. We validated these models as well computer simulation models on the results of other researchers and verify them on experimental data. The use of analytical model is essential because in some cases Nike2D supported simulation models do not provide this type of information, e.g. equations of state, dependence on initial and boundary conditions, solution stability, internal and external strain, stress, loads and phonoelastic tensor properties. Phonoelastic analytical models allow connect strain measurements with experimental data on Raman data under high pressure. 2. Computational Models. Equations of Motion and State in Nike2D Let ρ be the density, u i are the displacements Ω is the continuum domain, b i is body force per unit volume, τ ij is the Cauchy s stress tensor. As always u i = u i t, ü i = 2 u i t 2. FE equations of motion in Ω: ρ ü i = b i + τ ij,j. Let continuum domain Ω has a boundary Γ u, where spatial displacements u i are defined, and also a boundary Γ τ, where stresses τ are defined. u Boundary conditions are i = ũ i Γ u, τ ij n j = τ i Γ τ.

3 Initial conditions are u i(0) = u 0i, u i (0) = u 0i. Rate of deformation tensor is d ij = 1 2 ( u i,j + u j,i ). The Cauchy tensor τ in general is a function of d ij and a set of history variables H and temperature T (we consider temperature to be constant in our investigation) τ ij = τ ij ( d ij, H, T ) = τ ij ( d ij, H). In our work we use two Nike2D material models: elastic and elastoplastic. Nike2D uses axes r and z. Equations of state in the way they usually defined in solid mechanics are not inbuilt in Nike2D. Sets of parameters are used for elastic and elastoplastic materials. These parameters are related through some definitive equations of general form which we present later. This leads us in our research to a rather strong commitment: try to use well known equations of state for DAC components as well as different types of analytical models valid for FEM DAC (e.g. Timoshenko-Goodier s model). We use Lagrangian and Eulerian theories in order to implement

4 some of Nielsen s stress-strain models into FEAM. In order to validate analytical models we compare stress-strain curves obtained by Nike2D FEM with the ones obtained from Timoshenko-Goodier s, Nielsen s, and Hanfland et al. s models. We also develop adjusting procedures with use of some controlling parameters of these models in order to inbuilt them into FEAM. Notice that use of analytical models in our FEAM research is perfectly justified because these models are not inbuilt in FEM supported by Nike2D. Process of DAC compression is described by components, geometry model, material model, motion model and equation of state. Maze language commands allow to describe material models for each of DAC components. Diamond anvil and pusher are considered to be elastic and each of them is described by density, Young s modulus, and Poisson s ratio. This model describes isotropic, linear elastic material behavior. Gasket and sample material are considered to be kinematic/isotropic elastic-plastic material and described by density, Young s modulus, Poisson s ratio, yield stress, hardening modulus, hardening parameter, number of points in stress-effective plastic strain curve, effective plastic strain, and effective stress. This model includes linear or nonlinear strain hardening. Kinematic and isotropic hardening models which describe sample and gasket material yield identical behavior under monotonic loading. We consider only double-beveled diamond anvils. Using Nike2D Maze input language and using axial and radial symmetry of DAC we describe all four components of DAC. Each component is described by its own geometry and its own mesh. Geometry includes a series of line definitions. Points are introduced in the line definitions. Using these definitions element and nodal topologies are generated by Maze. Each DAC FE is defined by four boundary, nodal points which are in nondeformed state connected by line segments. The material number describing material properties of the element is also included in the definition. From lines parts and regions are constructed. Mesh assembly is defined from all described parts. Process of DAC compression is described by components, geometry model, material model, motion model and equation of state. Sample of DAC geometry and mesh for all components is shown. Maze language

5 commands allow to describe material models for each of DAC components. Diamond anvil and pusher are considered to be elastic and each of them is described by density, Young s modulus, and Poisson s ratio. This model describes isotropic, linear elastic material behavior. Gasket and sample material are considered to be kinematic/isotropic elasticplastic material and described by density, Young s modulus, Poisson s ratio, yield stress, hardening modulus, hardening parameter, number of points in stress-effective plastic strain curve, effective plastic strain, and effective stress. This model includes linear or nonlinear strain hardening. Kinematic and isotropic hardening models which describe sample and gasket material yield identical behavior under monotonic loading. Parameters describing elastic and elastoplastic material behavior usually are considered as linear or nonlinear functions of pressure. 3. Results and Conclusions The following results were obtained: 1. The validation of the modeling approaches. Pressure distribution results obtained on DAC compression model were verified using existing literature data; 2. DAC compression model in Nike2D computer code was investigated for different types of boundary and initial conditions, verified on several practical applications and theoretical results including equations of state, Timoshenko-Goodier s, Nielsen s and Hanfland et al. s models; 3. The experimentally measured pressure distribution in DAC to a peak pressure of 213 GPa by X-ray diffraction methods using a diamond coated rhenium gasket was compared with the finite element modeling results and a good agreement was obtained; 4. The role of a thin film of 13 C layer and its use as a universal sensor in DAC to 156 GPa and corresponding mechanisms involved were analyzed; 5. Finally, the role of diamond geometry and gasket materials were investigated to get a realistic estimate of the ultra high pressure conditions that can be generated in diamond anvil cells. The radial and axial stresses as well as shear stresses were examined and a failure criterion for diamond anvils in high pressure devices was developed. Acknowledgments This material is based upon work supported by the Department of Energy (DOE) National Nuclear Security Administration (NNSA) under Grant No DE-FG52-06NA26168.

6 References 1. Kondrat yev A. I., Vohra Y. K. Finite-element modeling of stresses and strains in a diamond anvil cell device: case of a diamond-coated rhenium gasket // High Pres. Res Vol. 27, no. 3. P Kondrat yev A. I. Mathematical modeling and computer simulation of compression process. VDM Verlag, 2009.

Example-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium

Example-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 information

Nonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su. ME 563: Nonlinear Finite Element Analysis.

Nonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su. ME 563: Nonlinear Finite Element Analysis. ME 563: Nonlinear Finite Element Analysis Spring 2016 Nonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su Department of Mechanical and Nuclear Engineering,

More information

Constitutive models: Incremental plasticity Drücker s postulate

Constitutive 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 information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

More information

Enhancing Prediction Accuracy In Sift Theory

Enhancing 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 information

Onedimensional ratchetting Twodimensional ratchetting Response of unified models Response of 2M1C model Response of polycrystalline models

Onedimensional ratchetting Twodimensional ratchetting Response of unified models Response of 2M1C model Response of polycrystalline models Ratchetting Onedimensional ratchetting Twodimensional ratchetting Response of unified models Response of 2M1C model Response of polycrystalline models Still an open problem Ratchetting effects Primary

More information

Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study

Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study Biplab Chatterjee, Prasanta Sahoo 1 Department of Mechanical Engineering, Jadavpur University

More information

Plane Strain Test for Metal Sheet Characterization

Plane Strain Test for Metal Sheet Characterization Plane Strain Test for Metal Sheet Characterization Paulo Flores 1, Felix Bonnet 2 and Anne-Marie Habraken 3 1 DIM, University of Concepción, Edmundo Larenas 270, Concepción, Chile 2 ENS - Cachan, Avenue

More information

Theory of Plasticity. Lecture Notes

Theory of Plasticity. Lecture Notes Theory of Plasticity Lecture Notes Spring 2012 Contents I Theory of Plasticity 1 1 Mechanical Theory of Plasticity 2 1.1 Field Equations for A Mechanical Theory.................... 2 1.1.1 Strain-displacement

More information

ENGN 2290: Plasticity Computational plasticity in Abaqus

ENGN 2290: Plasticity Computational plasticity in Abaqus ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element

More information

ME 582 Advanced Materials Science. Chapter 2 Macromechanical Analysis of a Lamina (Part 2)

ME 582 Advanced Materials Science. Chapter 2 Macromechanical Analysis of a Lamina (Part 2) ME 582 Advanced Materials Science Chapter 2 Macromechanical Analysis of a Lamina (Part 2) Laboratory for Composite Materials Research Department of Mechanical Engineering University of South Alabama, Mobile,

More information

Application of nanoindentation technique to extract properties of thin films through experimental and numerical analysis

Application of nanoindentation technique to extract properties of thin films through experimental and numerical analysis Materials Science-Poland, Vol. 28, No. 3, 2010 Application of nanoindentation technique to extract properties of thin films through experimental and numerical analysis A. WYMYSŁOWSKI 1*, Ł. DOWHAŃ 1, O.

More information

Introduction to Engineering Materials ENGR2000. Dr. Coates

Introduction to Engineering Materials ENGR2000. Dr. Coates Introduction to Engineering Materials ENGR2 Chapter 6: Mechanical Properties of Metals Dr. Coates 6.2 Concepts of Stress and Strain tension compression shear torsion Tension Tests The specimen is deformed

More information

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics

More information

Stress analysis of a stepped bar

Stress analysis of a stepped bar Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.

More information

Concept Question Comment on the general features of the stress-strain response under this loading condition for both types of materials

Concept Question Comment on the general features of the stress-strain response under this loading condition for both types of materials Module 5 Material failure Learning Objectives review the basic characteristics of the uni-axial stress-strain curves of ductile and brittle materials understand the need to develop failure criteria for

More information

MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS

MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS 1 MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS Version 2011-01-14 Stress tensor Definition of traction vector (1) Cauchy theorem (2) Equilibrium (3) Invariants (4) (5) (6) or, written in terms of principal

More information

Continuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms

Continuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive

More information

Chapter 6: Mechanical Properties of Metals. Dr. Feras Fraige

Chapter 6: Mechanical Properties of Metals. Dr. Feras Fraige Chapter 6: Mechanical Properties of Metals Dr. Feras Fraige Stress and Strain Tension Compression Shear Torsion Elastic deformation Plastic Deformation Yield Strength Tensile Strength Ductility Toughness

More information

Contact Modeling of Rough Surfaces. Robert L. Jackson Mechanical Engineering Department Auburn University

Contact Modeling of Rough Surfaces. Robert L. Jackson Mechanical Engineering Department Auburn University Contact Modeling of Rough Surfaces Robert L. Jackson Mechanical Engineering Department Auburn University Background The modeling of surface asperities on the micro-scale is of great interest to those interested

More information

UNIT I SIMPLE STRESSES AND STRAINS

UNIT I SIMPLE STRESSES AND STRAINS Subject with Code : SM-1(15A01303) Year & Sem: II-B.Tech & I-Sem SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) UNIT I SIMPLE STRESSES

More information

1 Introduction. Abstract

1 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 information

Module 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains

Module 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains Introduction In this lecture we are going to introduce a new micromechanics model to determine the fibrous composite effective properties in terms of properties of its individual phases. In this model

More information

Bulk Metal Forming II

Bulk Metal Forming II Bulk Metal Forming II Simulation Techniques in Manufacturing Technology Lecture 2 Laboratory for Machine Tools and Production Engineering Chair of Manufacturing Technology Prof. Dr.-Ing. Dr.-Ing. E.h.

More information

A PAPER ON DESIGN AND ANALYSIS OF PRESSURE VESSEL

A PAPER ON DESIGN AND ANALYSIS OF PRESSURE VESSEL A PAPER ON DESIGN AND ANALYSIS OF PRESSURE VESSEL P.Palanivelu 1, R.Siva Prasad 2, 1 PG Scholar, Department of Mechanical Engineering, Gojan School of Business and Technology, Redhills, Chennai, India.

More information

Esben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer

Esben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics

More information

Mechanics of Solids. Mechanics Of Solids. Suraj kr. Ray Department of Civil Engineering

Mechanics of Solids. Mechanics Of Solids. Suraj kr. Ray Department of Civil Engineering Mechanics Of Solids Suraj kr. Ray (surajjj2445@gmail.com) Department of Civil Engineering 1 Mechanics of Solids is a branch of applied mechanics that deals with the behaviour of solid bodies subjected

More information

BACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)

BACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM) BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple

More information

Load Cell Design Using COMSOL Multiphysics

Load Cell Design Using COMSOL Multiphysics Load Cell Design Using COMSOL Multiphysics Andrei Marchidan, Tarah N. Sullivan and Joseph L. Palladino Department of Engineering, Trinity College, Hartford, CT 06106, USA joseph.palladino@trincoll.edu

More information

4.MECHANICAL PROPERTIES OF MATERIALS

4.MECHANICAL PROPERTIES OF MATERIALS 4.MECHANICAL PROPERTIES OF MATERIALS The diagram representing the relation between stress and strain in a given material is an important characteristic of the material. To obtain the stress-strain diagram

More information

HSNV140 - Thermoplasticity with restoration of work hardening: test of blocked dilatometry

HSNV140 - Thermoplasticity with restoration of work hardening: test of blocked dilatometry Titre : HSNV140 - Thermo-plasticité avec restauration d'éc[...] Date : 29/06/2015 Page : 1/11 HSNV140 - Thermoplasticity with restoration of work hardening: test of blocked dilatometry Summary: This test

More information

CIVL4332 L1 Introduction to Finite Element Method

CIVL4332 L1 Introduction to Finite Element Method CIVL L Introduction to Finite Element Method CIVL L Introduction to Finite Element Method by Joe Gattas, Faris Albermani Introduction The FEM is a numerical technique for solving physical problems such

More information

ANSYS Mechanical Basic Structural Nonlinearities

ANSYS 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 information

Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method

Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method V.M.F. Nascimento Departameto de ngenharia Mecânica TM, UFF, Rio de Janeiro

More information

Nonlinear Finite Element Analysis of Gasketed Flange Joints under Combined Internal Pressure and Different Thermal Loading Conditions

Nonlinear Finite Element Analysis of Gasketed Flange Joints under Combined Internal Pressure and Different Thermal Loading Conditions Failure of Engineering Materials & Structures Code 50 UET TAXILA MECHNICAL ENGINEERING DEPARTMENT Nonlinear Finite Element Analysis of Gasketed Flange Joints under Combined Internal Pressure and Different

More information

ELASTICITY (MDM 10203)

ELASTICITY (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 information

Lectures on. Constitutive Modelling of Arteries. Ray Ogden

Lectures on. Constitutive Modelling of Arteries. Ray Ogden Lectures on Constitutive Modelling of Arteries Ray Ogden University of Aberdeen Xi an Jiaotong University April 2011 Overview of the Ingredients of Continuum Mechanics needed in Soft Tissue Biomechanics

More information

SSNV221 Hydrostatic test with a behavior DRUCK_PRAGER linear and parabolic

SSNV221 Hydrostatic test with a behavior DRUCK_PRAGER linear and parabolic Titre : SSNV221 - Essai hydrostatique avec un comportement[...] Date : 08/08/2011 Page : 1/10 SSNV221 Hydrostatic test with a behavior DRUCK_PRAGER linear and parabolic Summary: The case test proposes

More information

Curvature of a Cantilever Beam Subjected to an Equi-Biaxial Bending Moment. P. Krulevitch G. C. Johnson

Curvature of a Cantilever Beam Subjected to an Equi-Biaxial Bending Moment. P. Krulevitch G. C. Johnson UCRL-JC-30440 PREPRINT Curvature of a Cantilever Beam Subjected to an Equi-Biaxial Bending Moment P. Krulevitch G. C. Johnson This paper was prepared for submittal to the Materials Research Society Spring

More information

Chapter 7. Highlights:

Chapter 7. Highlights: Chapter 7 Highlights: 1. Understand the basic concepts of engineering stress and strain, yield strength, tensile strength, Young's(elastic) modulus, ductility, toughness, resilience, true stress and true

More information

Back Analysis of Measured Displacements of Tunnels

Back Analysis of Measured Displacements of Tunnels Rock Mechanics and Rock Engineering 16, 173--180 (1983) Rock Mechanics and Rock Engineering 9 by Springer-Verlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe

More information

INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS

INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS SHEN R. WU and LEI GU WILEY A JOHN WILEY & SONS, INC., PUBLICATION ! PREFACE xv PARTI FUNDAMENTALS 1 1 INTRODUCTION 3

More information

Chapter 5 Torsion STRUCTURAL MECHANICS: CE203. Notes are based on Mechanics of Materials: by R. C. Hibbeler, 7th Edition, Pearson

Chapter 5 Torsion STRUCTURAL MECHANICS: CE203. Notes are based on Mechanics of Materials: by R. C. Hibbeler, 7th Edition, Pearson STRUCTURAL MECHANICS: CE203 Chapter 5 Torsion Notes are based on Mechanics of Materials: by R. C. Hibbeler, 7th Edition, Pearson Dr B. Achour & Dr Eng. K. El-kashif Civil Engineering Department, University

More information

Strength of Material. Shear Strain. Dr. Attaullah Shah

Strength of Material. Shear Strain. Dr. Attaullah Shah Strength of Material Shear Strain Dr. Attaullah Shah Shear Strain TRIAXIAL DEFORMATION Poisson's Ratio Relationship Between E, G, and ν BIAXIAL DEFORMATION Bulk Modulus of Elasticity or Modulus of Volume

More information

Implementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis

Implementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty

More information

Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas

Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas San Miguel de Tucuman, Argentina September 14 th, 2011 Seminary on Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas Antonio Caggiano 1, Guillermo Etse 2, Enzo

More information

A Critical Plane-energy Model for Multiaxial Fatigue Life Prediction. of Homogeneous and Heterogeneous Materials. Haoyang Wei

A Critical Plane-energy Model for Multiaxial Fatigue Life Prediction. of Homogeneous and Heterogeneous Materials. Haoyang Wei A Critical Plane-energy Model for Multiaxial Fatigue Life Prediction of Homogeneous and Heterogeneous Materials by Haoyang Wei A Thesis Presented in Partial Fulfillment of the Requirements for the Degree

More information

MECHANICS OF MATERIALS

MECHANICS OF MATERIALS Third E CHAPTER 2 Stress MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Lecture Notes: J. Walt Oler Texas Tech University and Strain Axial Loading Contents Stress & Strain:

More information

Stress, Strain, Mohr s Circle

Stress, 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 information

SIMPLIFIED MODELING OF THIN-WALLED TUBES WITH OCTAGONAL CROSS SECTION AXIAL CRUSHING. Authors and Correspondance: Abstract:

SIMPLIFIED MODELING OF THIN-WALLED TUBES WITH OCTAGONAL CROSS SECTION AXIAL CRUSHING. Authors and Correspondance: Abstract: SIMPLIFIED MODELING OF THIN-WALLED TUBES WITH OCTAGONAL CROSS SECTION AXIAL CRUSHING Authors and Correspondance: Yucheng Liu, Michael L. Day Department of Mechanical Engineering University of Louisville

More information

Basic concepts to start Mechanics of Materials

Basic concepts to start Mechanics of Materials Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen

More information

6.730 Physics for Solid State Applications

6.730 Physics for Solid State Applications 6.730 Physics for Solid State Applications Lecture 5: Specific Heat of Lattice Waves Outline Review Lecture 4 3-D Elastic Continuum 3-D Lattice Waves Lattice Density of Modes Specific Heat of Lattice Specific

More information

AN IMPORTANT PITFALL OF PSEUDO-STATIC FINITE ELEMENT ANALYSIS

AN IMPORTANT PITFALL OF PSEUDO-STATIC FINITE ELEMENT ANALYSIS AN IMPORTANT PITFALL OF PSEUDO-STATIC FINITE ELEMENT ANALYSIS S. Kontoe, L. Pelecanos & D.M. Potts ABSTRACT: Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified

More information

MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD

MODELING 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 information

3.032 Problem Set 2 Solutions Fall 2007 Due: Start of Lecture,

3.032 Problem Set 2 Solutions Fall 2007 Due: Start of Lecture, 3.032 Problem Set 2 Solutions Fall 2007 Due: Start of Lecture, 09.21.07 1. In the beam considered in PS1, steel beams carried the distributed weight of the rooms above. To reduce stress on the beam, it

More information

Experimental and theoretical characterization of Li 2 TiO 3 and Li 4 SiO 4 pebbles

Experimental and theoretical characterization of Li 2 TiO 3 and Li 4 SiO 4 pebbles Experimental and theoretical characterization of Li 2 TiO 3 and Li 4 SiO 4 s D. Aquaro 1 N. Zaccari ABSTRACT Dipartimento di Ingegneria Meccanica Nucleare e della Produzione University of Pisa (Italy)

More information

Generic Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials

Generic Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Engineering Mechanics Dissertations & Theses Mechanical & Materials Engineering, Department of Winter 12-9-2011 Generic

More information

Figure 1 Lifting Lug Geometry with Weld

Figure 1 Lifting Lug Geometry with Weld Should you Perform Nonlinear Stress Analysis? Many of our clients inquire whether nonlinearity should be considered in their analyses. The answer to that question is not simple. Sometimes, as in certain

More information

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3 MECE 3321 MECHANICS OF SOLIDS CHAPTER 3 Samantha Ramirez TENSION AND COMPRESSION TESTS Tension and compression tests are used primarily to determine the relationship between σ avg and ε avg in any material.

More information

Experiment Two (2) Torsional testing of Circular Shafts

Experiment Two (2) Torsional testing of Circular Shafts Experiment Two (2) Torsional testing of Circular Shafts Introduction: Torsion occurs when any shaft is subjected to a torque. This is true whether the shaft is rotating (such as drive shafts on engines,

More information

EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES

EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES OF WOOD OR PLYWOOD CONSIDERED AS ORTHOTROPIC MATERIALS Information Revied and Reaffirmed March 1956 No. 1510 EFFECT OF ELLIPTIC

More information

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS H. Kumazawa and T. Takatoya Airframes and Structures Group, Japan Aerospace Exploration Agency 6-13-1, Ohsawa, Mitaka,

More information

An alternative design method for the double-layer combined die using autofrettage theory

An alternative design method for the double-layer combined die using autofrettage theory https://doi.org/10.5194/ms-8-267-2017 Author(s) 2017. This work is distributed under the Creative Commons Attribution.0 License. An alternative design method for the double-layer combined die using autofrettage

More information

DEVELOPMENT 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 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 information

Transactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN

Transactions on Engineering Sciences vol 14, 1997 WIT Press,  ISSN On the Computation of Elastic Elastic Rolling Contact using Adaptive Finite Element Techniques B. Zastrau^, U. Nackenhorst*,J. Jarewski^ ^Institute of Mechanics and Informatics, Technical University Dresden,

More information

*MAT_PAPER - a new orthotropic elastoplastic model for paper materials

*MAT_PAPER - a new orthotropic elastoplastic model for paper materials *MAT_PAPER - a new orthotropic elastoplastic model for paper materials Jesper Karlsson, Dynamore Nordic Mikael Schill, Dynamore Nordic Johan Tryding, Tetra Pak *MAT_PAPER (*MAT_274) A new orthotropic elastoplastic

More information

Temperature-dependent Thermal Stress Determination for Through-Silicon-Vias (TSVs) by Combining Bending Beam Technique with Finite Element Analysis

Temperature-dependent Thermal Stress Determination for Through-Silicon-Vias (TSVs) by Combining Bending Beam Technique with Finite Element Analysis Temperature-dependent Thermal Stress Determination for Through-Silicon-Vias (TSVs) by Combining Bending Beam Technique with Finite Element Analysis Kuan H. Lu, Suk-Kyu Ryu*, Qiu Zhao, Klaus Hummler**,

More information

Nonlinear Buckling Prediction in ANSYS. August 2009

Nonlinear Buckling Prediction in ANSYS. August 2009 Nonlinear Buckling Prediction in ANSYS August 2009 Buckling Overview Prediction of buckling of engineering structures is a challenging problem for several reasons: A real structure contains imperfections

More information

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ. Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,

More information

Modeling and Simulations of Aircraft Structures Stiffness, Damage, and Failure Prediction for Laminated Composites

Modeling and Simulations of Aircraft Structures Stiffness, Damage, and Failure Prediction for Laminated Composites Modeling and Simulations of Aircraft Structures Stiffness, Damage, and Failure Prediction for Laminated Composites H.E.Pettermann 1, C.Schuecker 1, D.H.Pahr 2, F.G.Rammerstorfer 2 1 Austrian Aeronautics

More information

MATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f

MATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown

More information

Continuum Mechanics. Continuum Mechanics and Constitutive Equations

Continuum Mechanics. Continuum Mechanics and Constitutive Equations Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform

More information

ANALYSIS AND NUMERICAL MODELLING OF CERAMIC PIEZOELECTRIC BEAM BEHAVIOR UNDER THE EFFECT OF EXTERNAL SOLICITATIONS

ANALYSIS AND NUMERICAL MODELLING OF CERAMIC PIEZOELECTRIC BEAM BEHAVIOR UNDER THE EFFECT OF EXTERNAL SOLICITATIONS Third International Conference on Energy, Materials, Applied Energetics and Pollution. ICEMAEP016, October 30-31, 016, Constantine,Algeria. ANALYSIS AND NUMERICAL MODELLING OF CERAMIC PIEZOELECTRIC BEAM

More information

Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models

Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models I. Rhee, K.J. Willam, B.P. Shing, University of Colorado at Boulder ABSTRACT: This paper examines the global

More information

6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa ( psi) and

6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa ( psi) and 6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa (15.5 10 6 psi) and an original diameter of 3.8 mm (0.15 in.) will experience only elastic deformation when a tensile

More information

TRESS - STRAIN RELATIONS

TRESS - 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 information

1.103 CIVIL ENGINEERING MATERIALS LABORATORY (1-2-3) Dr. J.T. Germaine Spring 2004 LABORATORY ASSIGNMENT NUMBER 6

1.103 CIVIL ENGINEERING MATERIALS LABORATORY (1-2-3) Dr. J.T. Germaine Spring 2004 LABORATORY ASSIGNMENT NUMBER 6 1.103 CIVIL ENGINEERING MATERIALS LABORATORY (1-2-3) Dr. J.T. Germaine MIT Spring 2004 LABORATORY ASSIGNMENT NUMBER 6 COMPRESSION TESTING AND ANISOTROPY OF WOOD Purpose: Reading: During this laboratory

More information

HONGJUN LI Department of Mechanical Engineering University of Strathclyde Glasgow, Scotland, UK

HONGJUN LI Department of Mechanical Engineering University of Strathclyde Glasgow, Scotland, UK HONGJUN LI Department of Mechanical Engineering University of Strathclyde Glasgow, Scotland, UK Introduction FEA is an established analysis method used in Pressure Vessel Design Most DBA is still based

More information

CONTACT MODEL FOR A ROUGH SURFACE

CONTACT MODEL FOR A ROUGH SURFACE 23 Paper presented at Bucharest, Romania CONTACT MODEL FOR A ROUGH SURFACE Sorin CĂNĂNĂU Polytechnic University of Bucharest, Dep. of Machine Elements & Tribology, ROMANIA s_cananau@yahoo.com ABSTRACT

More information

LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS

LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS XXII. LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS Introduction The lamination theory for the elastic stiffness of fiber composite materials is the backbone of the entire field, it holds

More information

1. A pure shear deformation is shown. The volume is unchanged. What is the strain tensor.

1. A pure shear deformation is shown. The volume is unchanged. What is the strain tensor. Elasticity Homework Problems 2014 Section 1. The Strain Tensor. 1. A pure shear deformation is shown. The volume is unchanged. What is the strain tensor. 2. Given a steel bar compressed with a deformation

More information

3-D Finite Element Analysis of Bolted Flange Joint of Pressure Vessel

3-D Finite Element Analysis of Bolted Flange Joint of Pressure Vessel ISSN No. 2230 7699 MIT Publications 35 3-D Finite Element Analysis of Bolted Flange Joint of Pressure Vessel Nomesh Kumar, P.V.G. Brahamanandam and B.V. Papa Rao Advanced Systems Laboratory, Defence Research

More information

DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION

DYNAMIC 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 information

Size Effects In the Crushing of Honeycomb Structures

Size Effects In the Crushing of Honeycomb Structures 45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference 19-22 April 2004, Palm Springs, California AIAA 2004-1640 Size Effects In the Crushing of Honeycomb Structures Erik C.

More information

Finite Element Modeling of an Aluminum Tricycle Frame

Finite Element Modeling of an Aluminum Tricycle Frame Finite Element Modeling of an Aluminum Tricycle Frame A. Rodríguez, B. Chiné*, and J. A. Ramírez Costa Rica Institute of Technology, School of Materials Science and Engineering, Cartago, Costa Rica *Corresponding

More information

Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE

Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE 1 Chapter 3 Load and Stress Analysis 2 Chapter Outline Equilibrium & Free-Body Diagrams Shear Force and Bending Moments in Beams Singularity Functions Stress Cartesian Stress Components Mohr s Circle for

More information

CHAPTER 3 THE EFFECTS OF FORCES ON MATERIALS

CHAPTER 3 THE EFFECTS OF FORCES ON MATERIALS CHAPTER THE EFFECTS OF FORCES ON MATERIALS EXERCISE 1, Page 50 1. A rectangular bar having a cross-sectional area of 80 mm has a tensile force of 0 kn applied to it. Determine the stress in the bar. Stress

More information

Proceedings of the ASME th International Conference on Ocean, Offshore and Arctic Engineering OMAE2016 June 19-24, 2016, Busan, South Korea

Proceedings of the ASME th International Conference on Ocean, Offshore and Arctic Engineering OMAE2016 June 19-24, 2016, Busan, South Korea Proceedings of the ASME 26 35th International Conference on Ocean, Offshore and Arctic Engineering OMAE26 June 9-24, 26, Busan, South Korea OMAE26-54554 LOCAL STRAIN AND STRESS CALCULATION METHODS OF IRREGULAR

More information

Finite element simulation of residual stresses in laser heating

Finite element simulation of residual stresses in laser heating IAS-2008-66-546ST Finite element simulation of residual stresses in laser heating G. H. Farrahi 1, M. Sistaninia 2, H. Moeinoddini 3 1,2-School of Mechanical Engineering, Sharif University of Technology,

More information

Unit I - Properties of Matter

Unit I - Properties of Matter Unit I - Properties of Matter Elasticity: Elastic and plastic materials Hooke s law elastic behavior of a material stress - strain diagram factors affecting elasticity. Three moduli of elasticity Poisson

More information

PLAXIS. Scientific Manual

PLAXIS. Scientific Manual PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3

More information

Mechanics of Irregular Honeycomb Structures

Mechanics of Irregular Honeycomb Structures Mechanics of Irregular Honeycomb Structures S. Adhikari 1, T. Mukhopadhyay 1 Chair of Aerospace Engineering, College of Engineering, Swansea University, Singleton Park, Swansea SA2 8PP, UK Sixth International

More information

After lecture 16 you should be able to

After lecture 16 you should be able to Lecture 16: Design of paper and board packaging Advanced concepts: FEM, Fracture Mechanics After lecture 16 you should be able to describe the finite element method and its use for paper- based industry

More information

Basic Principles of Weak Galerkin Finite Element Methods for PDEs

Basic Principles of Weak Galerkin Finite Element Methods for PDEs Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element

More information

GLOBAL AND LOCAL LINEAR BUCKLING BEHAVIOR OF A CHIRAL CELLULAR STRUCTURE

GLOBAL AND LOCAL LINEAR BUCKLING BEHAVIOR OF A CHIRAL CELLULAR STRUCTURE GLOBAL AND LOCAL LINEAR BUCKLING BEHAVIOR OF A CHIRAL CELLULAR STRUCTURE Alessandro Spadoni, Massimo Ruzzene School of Aerospace Engineering Georgia Institute of Technology Atlanta, GA 30332 Fabrizio Scarpa

More information

Coupled CFD-FE-Analysis for the Exhaust Manifold of a Diesel Engine

Coupled CFD-FE-Analysis for the Exhaust Manifold of a Diesel Engine Coupled CFD-FE-Analysis for the Exhaust Manifold of a Diesel Engine Yasar Deger*, Burkhard Simperl*, Luis P. Jimenez** *Sulzer Innotec, Sulzer Markets and Technology Ltd, Winterthur, Switzerland **Guascor

More information

Elasticity, the fourth-rank tensor defining the strain of crystalline

Elasticity, the fourth-rank tensor defining the strain of crystalline Elasticity of MgO and a primary pressure scale to 55 GPa Chang-Sheng Zha*, Ho-kwang Mao, and Russell J. Hemley Geophysical Laboratory and Center for High Pressure Research, Carnegie Institution of Washington,

More information

Module 2: Thermal Stresses in a 1D Beam Fixed at Both Ends

Module 2: Thermal Stresses in a 1D Beam Fixed at Both Ends Module 2: Thermal Stresses in a 1D Beam Fixed at Both Ends Table of Contents Problem Description 2 Theory 2 Preprocessor 3 Scalar Parameters 3 Real Constants and Material Properties 4 Geometry 6 Meshing

More information

. D CR Nomenclature D 1

. D CR Nomenclature D 1 . D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the

More information

Reference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity",

Reference material Reference books: Y.C. Fung, Foundations of Solid Mechanics, Prentice Hall R. Hill, The mathematical theory of plasticity, Reference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity", Oxford University Press, Oxford. J. Lubliner, "Plasticity

More information