MATERIAL ELASTIC HERRMANN INCOMPRESSIBLE command.

Size: px
Start display at page:

Download "MATERIAL ELASTIC HERRMANN INCOMPRESSIBLE command."

Transcription

1 MATERIAL ELASTIC HERRMANN INCOMPRESSIBLE command Synopsis The MATERIAL ELASTIC HERRMANN INCOMPRESSIBLE command is used to specify the parameters associated with an isotropic, linear elastic material idealization valid for incompressible as well as compressible materials The associated constitutive relation was first proposed in 1965 by L R Herrmann [1] The MATERIAL ELASTIC HERRMANN INCOMPRESSIBLE idealization is intended for use with mixed elements for the analysis of incompressible or nearly incompressible solids Consequently, this material idealization should be used in conjunction with mixed elements possessing discontinuous pressure approximations Presently APES contains the following elements of this type: Q4P1d quadrilateral elements Q9P3d quadrilateral elements T6P1d triangular elements T7P3d triangular elements H8P1d hexahedral elements H20P4d hexahedral elements H27P4d hexahedral elements Syntax The following syntax is used to describe a typical isotropic elastic Herrmann material idealization: MATerial ELAstic HERrmann incompressible NUMber # (DEScription string ) (DENsity ##) (LAMbda ## ) (SHEar modulus ##) 1 V N Kaliakin

2 Explanatory Notes The NUMBER keyword is used to specify the (global) number of the material associated with the incompressible isotropic elastic idealization The default material number is one (1) The default material DENSITY is zero (00) The optional alphanumeric string associated with the DESCRIPTION keyword must be enclosed in double quotes ( ) It is used solely to describe the material being idealized to the analyst The DESCRIPTION string is printed as part of the echo of the material information For mixed (as opposed to mixed/penalty) versions the Q4P1, Q9P3, T6P1 and H8P1 elements the keyword LAMBDA is used to input an appropriate value for the Lamé parameter λ Denoting the elastic modulus by E and Poisson s ratio by ν, λ = Eν (1 + ν)(1 2ν) The default value of LAMBDA is 250e+05 The inverse of λ is used in the element pressure submatrix; viz, ( ) 1 L (e) = λ (e) N(p)T N (p) dω Ω e where N (p) is a vector of pressure interpolation functions used to approximate the pressure in the element domain Ω e For the case of ν = 00 one cannot input λ = GAMMA = 00, for this would cause a divideby-zero runtime error In such cases, the analyst should simply input a small number (eg, 10) Past experience shows that the element will then yield results that are identical to those associated with the corresponding mixed elements with the incompressibility constraint exactly enforced (but with L (e) = 0 and thus a computationally more complex solution) As a material approaches the incompressible limit λ approaches infinity For slightly compressible materials, experience shows that the mixed versions of the Q4P1, Q9P3, T6P1 and H8P1 elements give accurate results If the mixed/penalty version of the Q4P1, Q9P3, T6P1 and H8P1 elements is used, the submatrix L (e) is replaced by (1/λ) I, where λ now represents a large but finite penalty number that approximately enforces the incompressibility constraint and I is an identity matrix of appropriate size In this case, LAMBDA should be taken equal to a value that is approximately one order of magnitude greater than the actual value of λ associated with the material idealization As the incompressible limit is approached, the check on the suitability of the value chosen for LAMBDA is the magnitude of the element volumetric strain, which should be zero (at least to the significant digits displayed in the result) 2 V N Kaliakin

3 The keyword SHEAR MODULUS is used to specify the elastic shear modulus G = µ, where E G = 2(1 + ν) The default SHEAR MODULUS is 1150e+05 Example of Command Usage Assuming that E = 20 x 10 5 and ν = 04999, and noting that and that λ = Eν (1 + ν)(1 2ν) = (20 x 10 5 )(04999) = 3333 x 109 ( )(1 2(04999)) G = E 2(1 + ν) = (20 x 105 ) = 6667 x 105 2( ) the following commands are sufficient to specify an incompressible Herrmann isotropic elastic material for these parameter values: mat elas Herrmann_incomp num 1 des "nearly incompressible idealization" & lambda 3333e+09 shear 6667e+05 3 V N Kaliakin

4 Theoretical Considerations The isotropic, linear elastic material idealization valid for incompressible as well as compressible materials was first proposed in 1965 by L R Herrmann [1] This formulation is based on the observation that the main difficulty in applying standard displacement-based formulations to incompressible or nearly incompressible media lies in the determination of the mean stress or pressure, which is related to the volumetric part of the strain For this reason, it is convenient to separate this pressure from the total stress field and to treat it as a primary dependent variable independent from the displacement Utilyzing a special case of the Hellinger-Reissner mixed variational principle in which displacements and the average hydrostatic pressure were selected as primary dependent variables, Herrmann [1] proposed a mixed formulation consisting of the constitutive relation and the following constraint equation for the additional unknown p: σ ij = pδ ij + 2µε ij (1) u i,i + p λ = 0 (2) In equations (1) and (2) compressive pressure is taken as being negative For the incompressible case, ν = 1/2 λ = Equation (2) thus correctly enforces the kinematic constraint of u i,i = 0 For compressible materials ν < 1/2 and λ is finite Solving equation (2) gives p = λu i,i Substituting this relation into equation (1), and noting that δ ij u k,k = ε 11 + ε 22 + ε 33 = ε kk gives σ ij = λδ ij ε kk + 2µε ij (3) which is the standard constitutive relation for isotropic linear elasticity Remark 1 The pressure p appearing in equations (1) and (2) does not necessarily equal the mean pressure From equation (3) it follows that σ ii = 3λε kk + 2µε ii = (3λ + 2µ) ε kk From the definition of the mean or hydrostatic pressure p, it follows that p = σ ii (λ 3 = + 23 ) µ ε kk = Kε kk Now from equation (2) p = λu i,i = λε kk Since K λ, it follows that for the compressible case, p does not represent the hydrostatic pressure Indeed, p λ = p K For the nearly incompressible case µ << λ, K λ, and thus p p Finally, for the incompressible case, the standard constitutive relations from isotropic elasticity do not apply, thus no 4 V N Kaliakin

5 meaningful relationship can be made However, for the mixed method, for λ =, equation (2) correctly gives u i,i = ε kk = 0 Thus, from equation (1) σ kk = 3p + 2µε kk = 3p showing that now p = σ kk /3 = p 5 V N Kaliakin

6 Bibliography [1] Herrmann, L R, Elasticity Equations for Nearly Incompressible Materials by a Variational Theorem, AIAA Journal, 3: (1965) 6

MATERIAL ELASTIC ANISOTROPIC command

MATERIAL ELASTIC ANISOTROPIC command MATERIAL ELASTIC ANISOTROPIC command.. Synopsis The MATERIAL ELASTIC ANISOTROPIC command is used to specify the parameters associated with an anisotropic linear elastic material idealization. Syntax The

More information

The following syntax is used to describe a typical irreducible continuum element:

The following syntax is used to describe a typical irreducible continuum element: ELEMENT IRREDUCIBLE T7P0 command.. Synopsis The ELEMENT IRREDUCIBLE T7P0 command is used to describe all irreducible 7-node enhanced quadratic triangular continuum elements that are to be used in mechanical

More information

Finite Element Method in Geotechnical Engineering

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 information

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations 6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =

More information

3D Elasticity Theory

3D 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 information

Lecture 4 Implementing material models: using usermat.f. Implementing User-Programmable Features (UPFs) in ANSYS ANSYS, Inc.

Lecture 4 Implementing material models: using usermat.f. Implementing User-Programmable Features (UPFs) in ANSYS ANSYS, Inc. Lecture 4 Implementing material models: using usermat.f Implementing User-Programmable Features (UPFs) in ANSYS 1 Lecture overview What is usermat.f used for? Stress, strain and material Jacobian matrix

More information

Lecture notes Models of Mechanics

Lecture notes Models of Mechanics Lecture notes Models of Mechanics Anders Klarbring Division of Mechanics, Linköping University, Sweden Lecture 7: Small deformation theories Klarbring (Mechanics, LiU) Lecture notes Linköping 2012 1 /

More information

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16. CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo

More information

Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature

Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte

More information

Basic Equations of Elasticity

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

Common pitfalls while using FEM

Common pitfalls while using FEM Common pitfalls while using FEM J. Pamin Instytut Technologii Informatycznych w Inżynierii Lądowej Wydział Inżynierii Lądowej, Politechnika Krakowska e-mail: JPamin@L5.pk.edu.pl With thanks to: R. de Borst

More information

3.2 Hooke s law anisotropic elasticity Robert Hooke ( ) Most general relationship

3.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 information

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

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

DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD

DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro

More information

Constitutive models. Constitutive model: determines P in terms of deformation

Constitutive models. Constitutive model: determines P in terms of deformation Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function

More information

Macroscopic theory Rock as 'elastic continuum'

Macroscopic theory Rock as 'elastic continuum' Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave

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

Stretching of a Prismatic Bar by its Own Weight

Stretching of a Prismatic Bar by its Own Weight 1 APES documentation (revision date: 12.03.10) Stretching of a Prismatic Bar by its Own Weight. This sample analysis is provided in order to illustrate the correct specification of the gravitational acceleration

More information

Mechanics of Earthquakes and Faulting

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

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004 Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic

More information

Module 2 Stresses in machine elements

Module 2 Stresses in machine elements Module 2 Stresses in machine elements Lesson 3 Strain analysis Instructional Objectives At the end of this lesson, the student should learn Normal and shear strains. 3-D strain matri. Constitutive equation;

More information

Fundamentals of Linear Elasticity

Fundamentals of Linear Elasticity Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy

More information

PREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO

PREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG

More information

A Stabilized Formulation for Incompressible Elasticity Using Linear Displacement and Pressure Interpolations

A Stabilized Formulation for Incompressible Elasticity Using Linear Displacement and Pressure Interpolations A Stabilized Formulation for Incompressible Elasticity Using Linear Displacement and Pressure Interpolations M. Chiumenti y, Q. Valverde Λ, C. Agelet de Saracibar y and M. Cervera y ydept. de Resistència

More information

20. Rheology & Linear Elasticity

20. 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 information

Lecture 8: Tissue Mechanics

Lecture 8: Tissue Mechanics Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials

More information

Resource reading (If you would like to brush up on stress, strain, and elasticity):

Resource reading (If you would like to brush up on stress, strain, and elasticity): 12.108 Struct. & Prop s. of Earth Mat. Lecture 9 1 Single-crystal elasticity Assigned Reading: Nye JF (1957) Physical Properties of Crystals. Oxford University Press, Oxford, UK (Chapters 5 and 6). Nye,

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

Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet

Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x

More information

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design MAE4700/5700 Finite Element Analsis for Mechanical and Aerospace Design Cornell Universit, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laborator Sible School of Mechanical and Aerospace

More information

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation

More information

Continuum Mechanics and the Finite Element Method

Continuum Mechanics and the Finite Element Method Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after

More information

Basic Concepts of Strain and Tilt. Evelyn Roeloffs, USGS June 2008

Basic Concepts of Strain and Tilt. Evelyn Roeloffs, USGS June 2008 Basic Concepts of Strain and Tilt Evelyn Roeloffs, USGS June 2008 1 Coordinates Right-handed coordinate system, with positions along the three axes specified by x,y,z. x,y will usually be horizontal, and

More information

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich

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

Mechanics of Biomaterials

Mechanics of Biomaterials Mechanics of Biomaterials Lecture 7 Presented by Andrian Sue AMME498/998 Semester, 206 The University of Sydney Slide Mechanics Models The University of Sydney Slide 2 Last Week Using motion to find forces

More information

Constitutive Equations

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

The Finite Element Method for Computational Structural Mechanics

The Finite Element Method for Computational Structural Mechanics The Finite Element Method for Computational Structural Mechanics Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) January 29, 2010 Martin Kronbichler (TDB) FEM for CSM January

More information

RHEOLOGY & LINEAR ELASTICITY. B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous isotropic materials

RHEOLOGY & LINEAR ELASTICITY. B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous isotropic materials GG303 Lecture 2 0 9/4/01 1 RHEOLOGY & LINEAR ELASTICITY I II Main Topics A Rheology: Macroscopic deformation behavior B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous

More information

MATERIAL PROPERTIES. Material Properties Must Be Evaluated By Laboratory or Field Tests 1.1 INTRODUCTION 1.2 ANISOTROPIC MATERIALS

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

1 Exercise: Linear, incompressible Stokes flow with FE

1 Exercise: Linear, incompressible Stokes flow with FE Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski

More information

Other state variables include the temperature, θ, and the entropy, S, which are defined below.

Other state variables include the temperature, θ, and the entropy, S, which are defined below. Chapter 3 Thermodynamics In order to complete the formulation we need to express the stress tensor T and the heat-flux vector q in terms of other variables. These expressions are known as constitutive

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139

MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 MASSACHUSETTS NSTTUTE OF TECHNOLOGY DEPARTMENT OF MATERALS SCENCE AND ENGNEERNG CAMBRDGE, MASSACHUSETTS 39 3. MECHANCAL PROPERTES OF MATERALS PROBLEM SET SOLUTONS Reading Ashby, M.F., 98, Tensors: Notes

More information

Mathematical Background

Mathematical Background CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous

More information

Hydrogeophysics - Seismics

Hydrogeophysics - Seismics Hydrogeophysics - Seismics Matthias Zillmer EOST-ULP p. 1 Table of contents SH polarized shear waves: Seismic source Case study: porosity of an aquifer Seismic velocities for porous media: The Frenkel-Biot-Gassmann

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

Some Remarks on the Reissner Mindlin Plate Model

Some Remarks on the Reissner Mindlin Plate Model Some Remarks on the Reissner Mindlin Plate Model Alexandre L. Madureira LNCC Brazil Talk at the Workshop of Numerical Analysis of PDEs LNCC February 10, 2003 1 Outline The 3D Problem and its Modeling Full

More information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals of Fluid Dynamics: Elementary Viscous Flow Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research

More information

RHEOLOGY & LINEAR ELASTICITY

RHEOLOGY & LINEAR ELASTICITY GG303 Lecture 20 10/25/09 1 RHEOLOGY & LINEAR ELASTICITY I Main Topics A Rheology: Macroscopic deformation behavior B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous

More information

INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA

INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA Problems in Solid Mechanics A Symposium in Honor of H.D. Bui Symi, Greece, July 3-8, 6 INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA M. HORI (Earthquake Research

More information

Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros

Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration

More information

The New Boundary Condition Effect on The Free Vibration Analysis of Micro-beams Based on The Modified Couple Stress Theory

The New Boundary Condition Effect on The Free Vibration Analysis of Micro-beams Based on The Modified Couple Stress Theory International Research Journal of Applied and Basic Sciences 2015 Available online at www.irjabs.com ISSN 2251-838X / Vol, 9 (3): 274-279 Science Explorer Publications The New Boundary Condition Effect

More information

Variational principles in mechanics

Variational principles in mechanics CHAPTER Variational principles in mechanics.1 Linear Elasticity n D Figure.1: A domain and its boundary = D [. Consider a domain Ω R 3 with its boundary = D [ of normal n (see Figure.1). The problem of

More information

Intrinsic finite element modeling of a linear membrane shell problem

Intrinsic finite element modeling of a linear membrane shell problem arxiv:3.39v [math.na] 5 Mar Intrinsic finite element modeling of a linear membrane shell problem Peter Hansbo Mats G. Larson Abstract A Galerkin finite element method for the membrane elasticity problem

More information

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

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Fig. 1. Circular fiber and interphase between the fiber and the matrix. Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In

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

Verification of the Hyperbolic Soil Model by Triaxial Test Simulations

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

Physics of Continuous media

Physics of Continuous media Physics of Continuous media Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 October 26, 2012 Deformations of continuous media If a body is deformed, we say that the point which originally had

More information

Mechanics PhD Preliminary Spring 2017

Mechanics PhD Preliminary Spring 2017 Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n

More information

Chapter 2. General concepts. 2.1 The Navier-Stokes equations

Chapter 2. General concepts. 2.1 The Navier-Stokes equations Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work

More information

A Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements

A Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements A Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment

More information

Cubzac-les-Ponts Experimental Embankments on Soft Clay

Cubzac-les-Ponts Experimental Embankments on Soft Clay Cubzac-les-Ponts Experimental Embankments on Soft Clay 1 Introduction In the 197 s, a series of test embankments were constructed on soft clay at Cubzac-les-Ponts in France. These full-scale field tests

More information

Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth

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

A super convergent hybridisable discontinuous Galerkin method for linear elasticity

A super convergent hybridisable discontinuous Galerkin method for linear elasticity A super convergent hybridisable discontinuous Galerkin method for linear elasticity arxiv:82.223v [math.na] 6 Feb 28 Ruben Sevilla Zienkiewicz Centre for Computational Engineering, College of Engineering,

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

MOOSE Workshop MOOSE-PF/MARMOT Training Mechanics Part Nashville, February 14th Daniel Schwen

MOOSE Workshop MOOSE-PF/MARMOT Training Mechanics Part Nashville, February 14th Daniel Schwen www.inl.gov MOOSE Workshop MOOSE-PF/MARMOT Training Mechanics Part Nashville, February 14th Daniel Schwen www.inl.gov III. Mechanics Mechanics Outline 1. Introduction 2. Tensor-based System 1. Compute

More information

Understand basic stress-strain response of engineering materials.

Understand basic stress-strain response of engineering materials. Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities

More information

Lecture 8. Stress Strain in Multi-dimension

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

1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity

1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity 1.050 Engineering Mechanics Lecture 22: Isotropic elasticity 1.050 Content overview I. Dimensional analysis 1. On monsters, mice and mushrooms 2. Similarity relations: Important engineering tools II. Stresses

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

Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity

Lecture #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 information

MICROMECHANICS AND HOMOGENIZATION. Materials Containing Coated and Uncoated Spherical Inhomogeneities

MICROMECHANICS AND HOMOGENIZATION. Materials Containing Coated and Uncoated Spherical Inhomogeneities MICROMECHANICS AND HOMOGENIZATION Materials Containing Coated and Uncoated Spherical Inhomogeneities Václav Nežerka, 2012 i Acknowledgement I would like to thank my supervisor, Jan Zeman, who supported

More information

GEOMETRIC NONLINEAR ANALYSIS

GEOMETRIC NONLINEAR ANALYSIS GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.

More information

- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)

- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum) 2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities

More information

pyoptfem Documentation

pyoptfem Documentation pyoptfem Documentation Release V0.0.6 F. Cuvelier November 09, 2013 CONTENTS 1 Presentation 3 1.1 Classical assembly algorithm (base version).............................. 6 1.2 Sparse matrix requirement........................................

More information

DEFORMATION PATTERN IN ELASTIC CRUST

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

Bending Load & Calibration Module

Bending Load & Calibration Module Bending Load & Calibration Module Objectives After completing this module, students shall be able to: 1) Conduct laboratory work to validate beam bending stress equations. 2) Develop an understanding of

More information

Poly3d Input File Format

Poly3d Input File Format Poly3d Input File Format Frantz Maerten, www.igeoss.com Version 1.0 2005, May 9th Contents 1 Conventions 2 1.1 Global coordinate system: global....................... 2 1.2 Element coordinate system: elocal......................

More information

SEISMOLOGY I. Laurea Magistralis in Physics of the Earth and of the Environment. Elasticity. Fabio ROMANELLI

SEISMOLOGY I. Laurea Magistralis in Physics of the Earth and of the Environment. Elasticity. Fabio ROMANELLI SEISMOLOGY I Laurea Magistralis in Physics of the Earth and of the Environment Elasticity Fabio ROMANELLI Dept. Earth Sciences Università degli studi di Trieste romanel@dst.units.it 1 Elasticity and Seismic

More information

ME FINITE ELEMENT ANALYSIS FORMULAS

ME FINITE ELEMENT ANALYSIS FORMULAS ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness

More information

Kirchhoff Plates: Field Equations

Kirchhoff Plates: Field Equations 20 Kirchhoff Plates: Field Equations AFEM Ch 20 Slide 1 Plate Structures A plate is a three dimensional bod characterized b Thinness: one of the plate dimensions, the thickness, is much smaller than the

More information

HÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling

HÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling HÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling TID: 2016-28-10, kl 14.00-19.00 Maximalt 60 poäng kan erhållas på denna tenta. För godkänt krävs 30 poäng. Tillåtet hjälpmedel: räknare

More information

Constitutive modelling of fabric anisotropy in sand

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

MIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires

MIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or

More information

Exercise: concepts from chapter 8

Exercise: 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

IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.

IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1. IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.978 PDF) http://web.mit.edu/mbuehler/www/teaching/iap2006/intro.htm

More information

Analysis of Planar Truss

Analysis of Planar Truss Analysis of Planar Truss Although the APES computer program is not a specific matrix structural code, it can none the less be used to analyze simple structures. In this example, the following statically

More information

Lecture 15 Strain and stress in beams

Lecture 15 Strain and stress in beams Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME

More information

Reduction of Finite Element Models of Complex Mechanical Components

Reduction of Finite Element Models of Complex Mechanical Components Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department

More information

Mechanical Properties of Materials

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

Deflection Analysis of Spur Polymer Gear Teeth

Deflection Analysis of Spur Polymer Gear Teeth Deflection Analysis of Spur Polymer Gear Teeth Boštjan Trobentar, Srečko Glodež Faculty of mechanical engineering University of Maribor Maribor, Slovenia bostjan.trobentar@student.um.si Boštjan Zafošnik

More information

Homework Problems. ( σ 11 + σ 22 ) 2. cos (θ /2), ( σ θθ σ rr ) 2. ( σ 22 σ 11 ) 2

Homework Problems. ( σ 11 + σ 22 ) 2. cos (θ /2), ( σ θθ σ rr ) 2. ( σ 22 σ 11 ) 2 Engineering Sciences 47: Fracture Mechanics J. R. Rice, 1991 Homework Problems 1) Assuming that the stress field near a crack tip in a linear elastic solid is singular in the form σ ij = rλ Σ ij (θ), it

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 Institute of Structural Engineering Method of Finite Elements II 1 Constitutive

More information

An orthotropic damage model for crash simulation of composites

An orthotropic damage model for crash simulation of composites High Performance Structures and Materials III 511 An orthotropic damage model for crash simulation of composites W. Wang 1, F. H. M. Swartjes 1 & M. D. Gan 1 BU Automotive Centre of Lightweight Structures

More information

Computational models of diamond anvil cell compression

Computational models of diamond anvil cell compression UDC 519.6 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

More information

Stress-Strain Behavior

Stress-Strain Behavior Stress-Strain Behavior 6.3 A specimen of aluminum having a rectangular cross section 10 mm 1.7 mm (0.4 in. 0.5 in.) is pulled in tension with 35,500 N (8000 lb f ) force, producing only elastic deformation.

More information

Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur

Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur-613 401 Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 9 Table of Contents 1. Plasticity:... 3 1.1 Plastic Deformation,

More information

Composite Structures. Indian Institute of Technology Kanpur

Composite Structures. Indian Institute of Technology Kanpur Mechanics of Laminated Composite Structures Nachiketa Tiwari Indian Institute of Technology Kanpur Lecture 23 Analysis of an Orthotropic Ply Lecture Overview Introduction Engineering constants for an 2

More information

Solving the torsion problem for isotropic matrial with a rectangular cross section using the FEM and FVM methods with triangular elements

Solving the torsion problem for isotropic matrial with a rectangular cross section using the FEM and FVM methods with triangular elements Solving the torsion problem for isotropic matrial with a rectangular cross section using the FEM and FVM methods with triangular elements Nasser M. Abbasi. June 0, 04 Contents Introduction. Problem setup...................................

More information