Energy Bounds in Linear Elasticity
|
|
- Lionel Chester Harper
- 6 years ago
- Views:
Transcription
1 Chapter 1 Energy Bounds in Linear Elasticity Focus of this Recitation: Thefocusofthisrecitationistoputtoworktheenergyboundsoflinearelasticity for the derivation of upper bound and lower bound energy solutions of the form: ) E com σ, ) T E com σ, T S.A. K.A. < = < E }{{ } ) E pot ε, ) pot ε, ξ ) 1.1) ξ }{{ } ) Ψ σ ) W T Ψε ) W ξ The main idea of these energy bounds is that the actual energy change that occurs in a material structural system can be approximated by means of the complementary energy E com σ, T ) and of the potential energy E pot ε, ξ ), which both realize a minimum for the solution: The complementary energy realizes a minimum among all possible statically admissible stress fieldss.a.): ) ) E com σ, T = min E σ com σ, T S.A. ) <E com σ, T 1.) S.A. The potential energy realizes a minimum among all possible kinematically admissible displacement fieldsk.a.): E pot ε, ) ξ = min E pot ε, ξ ) K.A. <E pot ε, ξ ) 1.3) ξ K.A. 18
2 1.1. REVIEW Review The two energy bounds defined by 1.1) provide a means toestimate the actual energy change in a material structural system Target For this actual exact solution, the solution or target energy can be evaluated from the sole knowledge of the external works, which is known as Clapeyron s formula: E com σ, T ) = 1 ) )) ) W T W ξ = E pot ξ,ε 1.4) ) wherew T istheworkbyprescribeddisplacements ξ d : ) W T = d ξ Tda 1.5) Ω ξ d ) andw ξ istheworkbyprescribedforcesρ g dω; T d da): ) W ξ = Ω ξ ρ g dω+ Ω T d ξ T d da 1.6) We will refer to this actual energy solution as the target solution Complementary Energy Approach = Stress Approach The complementary energy approach is a "stress approach", in which the complementaryenergye com σ, T ) is evaluated with statically admissible stress fieldss.a.): E com σ, T ) ) =Ψ σ ) W T 1.7)
3 130 CHAPTER 1. ENERGY BOUNDS IN LINEAR ELASTICITY whereψ σ )isthecomplementaryfreeenergy: ) Isotropic: Ψ σ 1 σ )= m Ω K +s dω G with: 1.8) σ m= 1 3 trσ ; s = 1 σ :σ 3σ m) ) andw T istheworkbyprescribeddisplacements ξ d definedby1.5). This stress approach involves the following steps: 1. Evaluate the target complementary energy solution from Clapeyron s formula1.4)[careful with sign!].. Chooseastressfieldσ andcheckthatitisstaticallyadmissible;thatis: on Ω T d : ons : T d =σ n T n)+ T n)=0 1.9) inω : divσ +ρ g =0 3. Determine the complementary energy from Eqs. 1.7),1.8) and1.5). 4. Minimize this complementary energy expression w.r.t. any"stress-degreeof-freedom" to obtain the lowest possible complementary energy for the chosenstressfieldσ. 5. Compare with Target : E com σ, T ) S.A. <E com σ, T ) 1.10) Potential Energy Approach = Displacement Approach The potential energy approach is a"displacement approach", in which the potentialenergye pot ε, ξ ) isevaluatedwithkinematicallyadmissibledisplacementfields ξ : E pot ε, ξ ) ) =Ψε ) W ξ 1.11)
4 1.. EXAMPLE 1: FROM LECTURE#8) 131 whereψε )isthehelmholtz)freeenergy: Isotropic: Ψε 1 )= Kǫ Ω v +Gǫ d) dω with: 1.1) ǫ v=trε ; ǫ d = ε :ε 1 ) 3 ǫ v andw ξ ) istheworkbyprescribedforcesdefinedby1.6). This displacement approach involves the following steps: 1. Evaluate the target potential energy solution from Clapeyron s formula 1.4).. Choose a displacement field ξ x) and check that it is kinematically admissible; that is on Ω ξ d: ξ d = ξ x) 1.13) Calculate the strain tensor: ε = 1 grad ξ +grad ) ξ ) T 1.14) 3. Determine the potential energy from Eqs. 1.11),1.1) and1.6). 4. Minimize this potential energy expression w.r.t. any"displacement-degreeof-freedom"inordertoobtainforthechosen ξ x)thesmallestpossible value of the potential energy. 5. Compare with Target : E pot ε, ξ ) K.A. <E pot ε, ξ ) 1.15) 1. Example 1: from lecture#8) 1..1 Problem Formulation We consider a homogeneous isotropic elastic soil layer of height h, subjected at itssurfacez=0toasurfacepressurep. Thefocusofthisexerciseistoestimate the surface settlement s. Using the energy approach, determine the following:
5 13 CHAPTER 1. ENERGY BOUNDS IN LINEAR ELASTICITY z p Figure 1.1: Soild layer subjected to a surface pressure. i) Target : Determine the potential energy and the complementary energy of the solution. ii) Lower Bound Approach: Determine an upper bound of the surface settlement by using the complementary energy approach, with a diagonal stress fieldoftheform: σ =a e x e x + e y e y )+b e z e z 1.16) wherea, b) are two stress constants to be determined. iii) Upper Bound Approach: Determine a lower bound of the displacement by using the potential energy approach, with a displacement field of the form: ξ =s 1 h) x e z 1.17) iv) Conclude this exercise by comparing the two displacement bounds so obtained. 1.. Target For the Target, we use Clapeyron s formula1.4) to evaluate the potential energy from the external work. In this problem, the work provided
6 1.. EXAMPLE 1: FROM LECTURE#8) 133 ) frotheoutsideisonlyduetothesurfacepressureatz=0;thusw T =0 and: ) E com σ, T = 1 ) W ξ = 1 ) psa = E pot ξ,ε whereaisthesurfacearea. 1.18) 1..3 Lower Bound Energy Approach = Upper Displacement Bound Thestressfieldσ mustbestaticallyadmissible. Fora,bconstant,theconditionσ needstosatisfyisthestressboundarycondition: z=0: T d =p e z σ e z 1.19) b= p The second stress parameter a remains undetermined. It is a"stress degree of freedom", and will be determined from minimizing the complementary energy. To calculate the complementary energy, we determine the stress invariants: σ m= 1 3 trσ = a p 1.0) 3 s = 1 σ :σ 3σ 1 m) = trσ σ ) 3σ m) = 1 ) ) a p a +p 3 3 = 3 ap+1 3 a p 1.1) [Note that σ m and s are constant in this problem, which facilitates the evaluation of subsequent integrals].
7 134 CHAPTER 1. ENERGY BOUNDS IN LINEAR ELASTICITY Thus, noting that W = 0 no work by prescribed displacements), we obtain from1.7) and1.8): E com σ, T ) ) = Ψ σ 1 σ )= m Ω K +s dω G = Ah 1 a 1 p) ap+ ) a p 1.) K G InordertoobtainthepossiblysmallestvalueofE com σ, T ) we minimize this function w.r.t. the"stress-degree-of-freedom", a: E com a) a = Ah 9GK a4g+3k) pg 3K))=0 1.3) a= G 3K 4G+3K p Resubstitution in the complementary energy expression yields: min a E com a)= 3 Ah p 4G+3K Finally, a comparison with the target solution yields: E com σ, T )= 1 psa< 3 Ah p 4G+3K =E com σ, T ) 1.4) 1.5) p s<3h 4G+3K Thus, the lower energy bound approach based on statically admissible stress fields yields an upper bound for the displacement Upper Bound Energy Approach = Lower Displacement Bound Thedisplacementfield ξ mustbekinematicallyadmissible,whichisthecase here, since: z=h: ξ d =0 s e z 1 z ) 1.6) h z=h
8 1.. EXAMPLE 1: FROM LECTURE#8) 135 In this displacement approach, s is a "displacement degree of freedom", which will be determined from minimizing the potential energy. To determine the potential energy1.11), we note that the external work by prescribed forces1.6) is: W ) ξ = ξ T d da=ap ξ e z =Aps 1.7) Ω T d To determine the free energy, we determine the strain invariants: ε = ε zz e z e z = s e z e z 1.8) h ǫ v = trε = s 1.9) h ǫ d = ε :ε 1 ) 3 ǫ v = 4 ) s 1.30) 3 h Hence,notingthatǫ v andǫ d areconstants: Ψε )=Ah 1 ) Kǫ v +Gǫ 1 d =Ah s h ) K+G 4 ) 3 The potential energy thus is: E pot ε, ξ ) ) = Ψε ) W ξ [ = A h 1 ) s ] K+G 4 ) ps h ) 1.3) InordertoobtainthesmallestvalueofE pot ε, ξ ) weminimizetheexpression w.r.t. the"displacement-degree-of-freedom": E pot s ) s =A s h ) s = h ) K+G 4 ) ] p = ) p ) K+G 4 3
9 136 CHAPTER 1. ENERGY BOUNDS IN LINEAR ELASTICITY The smallest potential energy for the chosen displacement field is thus: mine pot s )= 3 p ha s 4G+3K Finally, a comparison with the target solution yields: 1.34) ) E pot ξ,ε = 1 psa< 3 p ha 4G+3K =min E pot s ) s 1.35) s s ph = ) K+G Comparison: Lower Upper Bound Energy and Displacements Letusbringthedifferentenergyboundstogetherinthesenseofthefundamental relation1.1): mine com a)= 3 a Ah p 4G+3K ) E com σ, T S.A. < = 1 psa = E pot ε, ) ξ K.A. < min s E pot s )= 3 p ha 4G+3K S.A. < K.A. < 1.36) From a comparison of the energies, we readily realize that the lower energy bound min a E com a) is equaltothe upper bound K.A. < min s E pot s ). This means that the chosen stress field and displacement field are actual the solution of the problem. This is also found by comparing the displacement bounds. Since the upper energy bound provides lower settlement boundi.e. it overestimates the rigidity), and since the lower energy bound provides an upper settlement bound, the very fact that the settlements obtained with the two different approaches coincide, means that this is the actual solution of the problem: p 3h 4G+3K S.A. s K.A. s = ph K+G 4 3 ) 1.37)
10 1.3. EXAMPLE : DEEP TUNNELING THROUGH HETEROGENEOUS ROCK Example : Deep Tunneling Through Heterogeneous Rock What s the Meaning of Heterogeneous XYZ The previous example was restricted to homogeneous materials, meaning thattheelasticpropertiesdonotvaryfromonepointtotheother. However, the same approach also applied to heterogeneous materials and structures, for which the elastic properties vary from one point to another: K=K x);g=g x) 1.38) The problem is quite complex, and there are often no closed form solutions possible for such heterogeneous problems. The energy bounds then provide a unique opportunity to develop upper and lower energy bounds. In this case, the spatial variability of the elastic properties must be considered in the evaluation of the free energy integral1.1) and the complementary free energy integral1.8): Ψε ) = Ψ σ ) = Ω Ω Problem Formulation K x) ǫ v +G x) ǫ σ m K x) + s G x) d) dω 1.39) ) dω 1.40) WeconsiderthetunnelingofacylinderofradiusRisathighdepthH R) throughamountainthat is composedof two rocks separatedalong a plane Γ Fig. 1.). Giventhehighdepth,thetworockscanbemodeledastwosemiinfinite domains which are perfectly connected along the separation plane Γ. The tunnel is driven in the z-direction, normal to the separation plane Γ, which separates the tunnel section into two half-circles. Throughout this exercise, we willusecylindercoordinatesr,θ,z),withzbeingthetunnelaxis. We have seen a similar problem in a previous recitation. However, the main difference here is that the surrounding rock is not homogeneousconstant elasticity properties), but composed of two different rocks having different elastic properties: { K1,G Kθ),Gθ))= 1 ); 0<θ<π 1.41) K,G ); π<θ<π
11 138 CHAPTER 1. ENERGY BOUNDS IN LINEAR ELASTICITY H >> R r r K,G 1 1 θ Γ K 1,G 1 R θ Γ K,G K,G σ = 01 0 p δ Figure 1.: Deep tunneling through heterogeneous rock: two infinite halfspaces separated along a plane Γ.
12 1.3. EXAMPLE : DEEP TUNNELING THROUGH HETEROGENEOUS ROCK139 In contrast to the homogeneous deep tunneling problem, there is no closed form solution of this elasticity problem. By putting the energy bounds to work, we want to provide a back-of-the-enveloppe estimate of the on-average elastic radial displacementpositive inward) that occurs at the tunnel wall; that is an upper and a lower bound of the displacementfig. 1.): where n= e r. δ= 1 π π 0 ξr=r,θ) ndθ 1.4) 1. Target : By using Clapeyron s formula, show that the energy of the actual solution is equal to: E pot ξ,ε ) = = p 0 Rπδ 1.43) wherep 0 =ρ 0 gh isthehydrostaticgeological)pressureprevailinginthe rock before tunneling.. Upper Bound Approach: Consider a radial displacement field of the form: ξ =ur) er = C r e r;i=1, 1.44) where C is a"displacement-degree-of-freedom". Using the potential energyapproach,determinealowerboundofthedisplacementδ. Whyis ξ not the displacement solution of the heterogeneous problem? 3. Lower Bound Approach: Consider a radial stress field of the form: σ σ 0=p 0 R r) e r e r + e θ e θ ) 1.45) whereσ 0= p 0 1isthegeostaticstressstate. Usingthecomplementary energy approach, determine an upper bound of the displacement δ. Why isσ notthesolutionoftheheterogeneousproblem? [The solution is developed in the lecture notes, pp ]
1.033/1.57 Q#2: Elasticity Bounds Conical Indentation Test
1.033/1.57 Q#: Elasticity Bounds Conical Indentation Test November 14, 003 MIT 1.033/1.57 Fall 003 Instructor: Franz-Josef UM Instrumented nano-indentation is a new technique in materials science and engineering
More information1 Hooke s law, stiffness, and compliance
Non-quilibrium Continuum Physics TA session #5 TA: Yohai Bar Sinai 3.04.206 Linear elasticity I This TA session is the first of three at least, maybe more) in which we ll dive deep deep into linear elasticity
More informationWe briefly discuss two examples for solving wave propagation type problems with finite differences, the acoustic and the seismic problem.
Excerpt from GEOL557 Numerical Modeling of Earth Systems by Becker and Kaus 2016 1 Wave propagation Figure 1: Finite difference discretization of the 2D acoustic problem. We briefly discuss two examples
More informationChapter 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 informationLASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE
LASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE H. M. Al-Qahtani and S. K. Datta University of Colorado Boulder CO 839-7 ABSTRACT. An analysis of the propagation of thermoelastic waves
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 informationHomework 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 informationMechanics 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 informationFinite Elements for Thermal Analysis
Chapter 4 Finite Elements for Thermal Analysis 4.1 Governing Equations for the Thermal Problem CONSTITUTIVE EQUATIONS Introducing the temperature field T(x i,t ), the heat flux vector q(x i,t) q n (x i,t)
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 informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationSome 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 information1.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 informationRocking behaviour of a rigid foundation with an arbitrary embedment
Rocking behaviour of a rigid foundation with an arbitrary embedment H. Moghaddasi Islamic Azad University, Qazvin, Iran. M. Kalantari University of Tehran, Tehran, Iran N. Chouw The University of Auckland,
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationLecture Notes #10. The "paradox" of finite strain and preferred orientation pure vs. simple shear
12.520 Lecture Notes #10 Finite Strain The "paradox" of finite strain and preferred orientation pure vs. simple shear Suppose: 1) Anisotropy develops as mineral grains grow such that they are preferentially
More information1.033/1.57 Q#1: Stress & Strength Conical Indentation Tests
1.033/1.57 Q#1: Stress & Strength Conical Indentation Tests October 8, 003 MIT 1.033/1.57 Fall 003 Instructor: Franz-Josef ULM Instrumented nano-indentation is a new technique in materials science and
More informationFig. 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 informationContinuum 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 informationDynamics 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 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 informationChapter 2 Examples of Optimization of Discrete Parameter Systems
Chapter Examples of Optimization of Discrete Parameter Systems The following chapter gives some examples of the general optimization problem (SO) introduced in the previous chapter. They all concern the
More informationModule 9 : Foundation on rocks
LECTURE 32 9.3 BEARING CAPCITY contd... 9.3.2 Safe bearing pressure There are different methods are available to determine the safe bearing pressure on rocks. However, the applicability of different methods
More informationDistributed: Wednesday, March 17, 2004
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 019.00 MECHANICS AND MATERIALS II QUIZ I SOLUTIONS Distributed: Wednesday, March 17, 004 This quiz consists
More informationPoint Load Generated Surface Waves in a Transversely Isotropic Half-Space using Elastodynamic Reciprocity
Advances in Copyright 4 Tech Science Press 5 Point Load Generated Surface Waves in a Transversely Isotropic Half-Space using Elastodynamic eciprocity A. C. Wijeyewickrema 1, M. Hattori 1, S. Leungvichcharoen
More informationProceedings 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 information2 Basic Equations in Generalized Plane Strain
Boundary integral equations for plane orthotropic bodies and exterior regions G. Szeidl and J. Dudra University of Miskolc, Department of Mechanics 3515 Miskolc-Egyetemváros, Hungary Abstract Assuming
More informationModeling of Fiber-Reinforced Membrane Materials Daniel Balzani. (Acknowledgement: Anna Zahn) Tasks Week 2 Winter term 2014
Institute of echanics and Shell Structures Faculty Civil Engineering Chair of echanics odeling of Fiber-Reinforced embrane aterials OOC@TU9 Daniel Balani (Acknowledgement: Anna Zahn Tasks Week 2 Winter
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 informationIntroduction to the J-integral
Introduction to the J-integral Instructor: Ramsharan Rangarajan February 24, 2016 The purpose of this lecture is to briefly introduce the J-integral, which is widely used in fracture mechanics. To the
More informationMechanics of materials Lecture 4 Strain and deformation
Mechanics of materials Lecture 4 Strain and deformation Reijo Kouhia Tampere University of Technology Department of Mechanical Engineering and Industrial Design Wednesday 3 rd February, 206 of a continuum
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationBasic 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 informationInfluence of Irregularity and Rigidity on the Propagation of Torsional Wave
Applied Mathematical Sciences, Vol. 4, 00, no. 7, 805-86 Influence of Irregularity and Rigidity on the Propagation of Torsional Wave S. Gupta, A. Chattopadhyay and S. Kundu Department of Applied Mathematics
More informationFundamentals 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 informationVariational 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. 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 informationparticle p = m v F ext = d P = M d v cm dt
Lecture 11: Momentum and Collisions; Introduction to Rotation 1 REVIEW: (Chapter 8) LINEAR MOMENTUM and COLLISIONS The first new physical quantity introduced in Chapter 8 is Linear Momentum Linear Momentum
More information1.5 STRESS-PATH METHOD OF SETTLEMENT CALCULATION 1.5 STRESS-PATH METHOD OF SETTLEMENT CALCULATION
Module 6 Lecture 40 Evaluation of Soil Settlement - 6 Topics 1.5 STRESS-PATH METHOD OF SETTLEMENT CALCULATION 1.5.1 Definition of Stress Path 1.5. Stress and Strain Path for Consolidated Undrained Undrained
More informationWater is sloshing back and forth between two infinite vertical walls separated by a distance L: h(x,t) Water L
ME9a. SOLUTIONS. Nov., 29. Due Nov. 7 PROBLEM 2 Water is sloshing back and forth between two infinite vertical walls separated by a distance L: y Surface Water L h(x,t x Tank The flow is assumed to be
More informationChapter 0. Preliminaries. 0.1 Things you should already know
Chapter 0 Preliminaries These notes cover the course MATH45061 (Continuum Mechanics) and are intended to supplement the lectures. The course does not follow any particular text, so you do not need to buy
More information2017 Soil Mechanics II and Exercises Final Exam. 2017/7/26 (Wed) 10:00-12:00 Kyotsu 4 Lecture room
2017 Soil Mechanics II and Exercises Final Exam 2017/7/26 (Wed) 10:00-12:00 Kyotsu 4 Lecture room Attention: The exam consists of five questions for which you are provided with five answer sheets. Write
More informationS.P. Timoshenko Institute of Mechanics, National Academy of Sciences, Department of Fracture Mechanics, Kyiv, Ukraine
CALCULATION OF THE PREFRACTURE ZONE AT THE CRACK TIP ON THE INTERFACE OF MEDIA A. Kaminsky a L. Kipnis b M. Dudik b G. Khazin b and A. Bykovtscev c a S.P. Timoshenko Institute of Mechanics National Academy
More informationThe Finite Element Method
The Finite Element Method 3D Problems Heat Transfer and Elasticity Read: Chapter 14 CONTENTS Finite element models of 3-D Heat Transfer Finite element model of 3-D Elasticity Typical 3-D Finite Elements
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 informationReference 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 informationConstitutive Relations
Constitutive Relations Andri Andriyana, Ph.D. Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of
More informationPHYS 771, Quantum Mechanics, Final Exam, Fall 2011 Instructor: Dr. A. G. Petukhov. Solutions
PHYS 771, Quantum Mechanics, Final Exam, Fall 11 Instructor: Dr. A. G. Petukhov Solutions 1. Apply WKB approximation to a particle moving in a potential 1 V x) = mω x x > otherwise Find eigenfunctions,
More informationVariation Principle in Mechanics
Section 2 Variation Principle in Mechanics Hamilton s Principle: Every mechanical system is characterized by a Lagrangian, L(q i, q i, t) or L(q, q, t) in brief, and the motion of he system is such that
More informationRutgers University Department of Physics & Astronomy. 01:750:271 Honors Physics I Fall Lecture 19. Home Page. Title Page. Page 1 of 36.
Rutgers University Department of Physics & Astronomy 01:750:271 Honors Physics I Fall 2015 Lecture 19 Page 1 of 36 12. Equilibrium and Elasticity How do objects behave under applied external forces? Under
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 informationLecture II: Rigid-Body Physics
Rigid-Body Motion Previously: Point dimensionless objects moving through a trajectory. Today: Objects with dimensions, moving as one piece. 2 Rigid-Body Kinematics Objects as sets of points. Relative distances
More informationElements 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 informationSurface force on a volume element.
STRESS and STRAIN Reading: Section. of Stein and Wysession. In this section, we will see how Newton s second law and Generalized Hooke s law can be used to characterize the response of continuous medium
More informationChapter 6: Momentum Analysis
6-1 Introduction 6-2Newton s Law and Conservation of Momentum 6-3 Choosing a Control Volume 6-4 Forces Acting on a Control Volume 6-5Linear Momentum Equation 6-6 Angular Momentum 6-7 The Second Law of
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More informationBACKGROUNDS. 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 informationSPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES
2010/2 PAGES 1 8 RECEIVED 21. 9. 2009 ACCEPTED 20. 1. 2010 Y. KOLEKOVÁ, M. PETRONIJEVIĆ, G. SCHMID SPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES ABSTRACT
More information7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim
69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the
More informationMicro-mechanics in Geotechnical Engineering
Micro-mechanics in Geotechnical Engineering Chung R. Song Department of Civil Engineering The University of Mississippi University, MS 38677 Fundamental Concepts Macro-behavior of a material is the average
More informationLinear Cosserat elasticity, conformal curvature and bounded stiffness
1 Linear Cosserat elasticity, conformal curvature and bounded stiffness Patrizio Neff, Jena Jeong Chair of Nonlinear Analysis & Modelling, Uni Dui.-Essen Ecole Speciale des Travaux Publics, Cachan, Paris
More informationConstitutive Relations
Constitutive Relations Dr. Andri Andriyana Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of field
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 informationA non-linear elastic/perfectly plastic analysis for plane strain undrained expansion tests
Bolton, M. D. & Whittle, R. W. (999). GeÂotechnique 49, No., 33±4 TECHNICAL NOTE A non-linear elastic/perfectly plastic analysis for plane strain undrained expansion tests M. D. BOLTON and R. W. WHITTLE{
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationChapter 6: Momentum Analysis of Flow Systems
Chapter 6: Momentum Analysis of Flow Systems Introduction Fluid flow problems can be analyzed using one of three basic approaches: differential, experimental, and integral (or control volume). In Chap.
More informationLecture 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 informationTopology Optimization of Compliant Mechanism with Geometrical Advantage
610 Topology Optimization of Compliant Mechanism with Geometrical Advantage Seungjae MIN and Younggi KIM A compliant mechanism is a mechanism that produces its motion by the flexibility of some or all
More informationStresses in Curved Beam
Stresses in Curved Beam Consider a curved beam subjected to bending moment M b as shown in the figure. The distribution of stress in curved flexural member is determined by using the following assumptions:
More informationProf. B V S Viswanadham, Department of Civil Engineering, IIT Bombay
29 Module 7: Lecture - 4 on Geotechnical Physical Modelling Variation of centrifugal acceleration in models of a 10 m soil layer [N =50] 1.0 m Merit of large centrifuge 4.0 m 18.4 47.4g 45g 50g 0.2 m 4.4
More informationFigure 2-1: Stresses under axisymmetric circular loading
. Stresses in Pavements.1. Stresses in Fleible Pavements.1.1. Stresses in Homogeneous Mass Boussinesq formulated models for the stresses inside an elastic half-space due to a concentrated load applied
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 informationLecture 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 informationDynamics 12e. Copyright 2010 Pearson Education South Asia Pte Ltd. Chapter 20 3D Kinematics of a Rigid Body
Engineering Mechanics: Dynamics 12e Chapter 20 3D Kinematics of a Rigid Body Chapter Objectives Kinematics of a body subjected to rotation about a fixed axis and general plane motion. Relative-motion analysis
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 informationElastic Scattering. R = m 1r 1 + m 2 r 2 m 1 + m 2. is the center of mass which is known to move with a constant velocity (see previous lectures):
Elastic Scattering In this section we will consider a problem of scattering of two particles obeying Newtonian mechanics. The problem of scattering can be viewed as a truncated version of dynamic problem
More informationUNIT IV FLEXIBILTY AND STIFFNESS METHOD
SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : SA-II (13A01505) Year & Sem: III-B.Tech & I-Sem Course & Branch: B.Tech
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 13 http://acousticalsociety.org/ ICA 13 Montreal Montreal, Canada - 7 June 13 Structural Acoustics and Vibration Session 4aSA: Applications in Structural
More informationMechanical Design in Optical Engineering
Torsion Torsion: Torsion refers to the twisting of a structural member that is loaded by couples (torque) that produce rotation about the member s longitudinal axis. In other words, the member is loaded
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 informationÜbungen zu RT2 SS (4) Show that (any) contraction of a (p, q) - tensor results in a (p 1, q 1) - tensor.
Übungen zu RT2 SS 2010 (1) Show that the tensor field g µν (x) = η µν is invariant under Poincaré transformations, i.e. x µ x µ = L µ νx ν + c µ, where L µ ν is a constant matrix subject to L µ ρl ν ση
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 informationThe Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density
Applied Mathematics & Information Sciences 23 2008, 237 257 An International Journal c 2008 Dixie W Publishing Corporation, U. S. A. The Rotating Inhomogeneous Elastic Cylinders of Variable-Thickness and
More information9 Kinetics of 3D rigid bodies - rotating frames
9 Kinetics of 3D rigid bodies - rotating frames 9. Consider the two gears depicted in the figure. The gear B of radius R B is fixed to the ground, while the gear A of mass m A and radius R A turns freely
More informationMATERIAL ELASTIC HERRMANN INCOMPRESSIBLE command.
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
More informationFundamental principles
Dynamics and control of mechanical systems Date Day 1 (03/05) - 05/05 Day (07/05) Day 3 (09/05) Day 4 (11/05) Day 5 (14/05) Day 6 (16/05) Content Review of the basics of mechanics. Kinematics of rigid
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 informationTransduction Based on Changes in the Energy Stored in an Electrical Field
Lecture 6-1 Transduction Based on Changes in the Energy Stored in an Electrical Field Electric Field and Forces Suppose a charged fixed q 1 in a space, an exploring charge q is moving toward the fixed
More informationMATH45061: SOLUTION SHEET 1 V
1 MATH4561: SOLUTION SHEET 1 V 1.) a.) The faces of the cube remain aligned with the same coordinate planes. We assign Cartesian coordinates aligned with the original cube (x, y, z), where x, y, z 1. The
More informationChapter 3. Load and Stress Analysis
Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3
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 informationMidterm Examination. Please initial the statement below to show that you have read it
EN75: Advanced Mechanics of Solids Midterm Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use two pages
More informationClusters in granular flows : elements of a non-local rheology
CEA-Saclay/SPEC; Group Instabilities and Turbulence Clusters in granular flows : elements of a non-local rheology Olivier Dauchot together with many contributors: D. Bonamy, E. Bertin, S. Deboeuf, B. Andreotti,
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationQuestion: Find displacement field ξ(x,t) and stress field σ(x,t) complying with the system of equations
4 An example of solution to a poblem of an elasto-plastic stuctue Small petubations : ε = 2 ( ξ + T ξ, Equilibium : div ( σ + F = 0, Constitutive laws: f (σ 0, ε = S : σ + ε p, ε p = λ f σ (σ, λ f (σ =
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 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 informationPLAXIS. 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 informationCHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.
1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy
More informationMacroscopic 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