Introduction to elastic wave equation. Salam Alnabulsi University of Calgary Department of Mathematics and Statistics October 15,2012

Similar documents
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

In this section is given an overview of the common elasticity models.

EQUATION CHAPTER 1 SECTION 1STRAIN IN A CONTINUOUS MEDIUM

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

MEMBRANE ELEMENT WITH NORMAL ROTATIONS

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

Mathematical Preparations

Lecture Note 3. Eshelby s Inclusion II

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

C PLANE ELASTICITY PROBLEM FORMULATIONS

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Study Guide For Exam Two

Introduction to Turbulence Modelling

PHYS 705: Classical Mechanics. Newtonian Mechanics

Spin-rotation coupling of the angularly accelerated rigid body

AE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

One Dimensional Axial Deformations

BAR & TRUSS FINITE ELEMENT. Direct Stiffness Method

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

Spring 2002 Lecture #13

CLOSED-FORM CHARACTERIZATION OF THE CHANNEL CAPACITY OF MULTI-BRANCH MAXIMAL RATIO COMBINING OVER CORRELATED NAKAGAMI FADING CHANNELS

Indeterminate pin-jointed frames (trusses)

Week 9 Chapter 10 Section 1-5

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

Homework Notes Week 7

A Porous Media Approach for Complex Heat and Fluid Flow System

Generalized form of reflection coefficients in terms of impedance matrices of qp-qp and qp-qs waves in TI media

On the symmetric character of the thermal conductivity tensor

Iterative General Dynamic Model for Serial-Link Manipulators

Handout: Large Eddy Simulation I. Introduction to Subgrid-Scale (SGS) Models

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

Linear Momentum. Center of Mass.

Physics 181. Particle Systems

On Vectors and Tensors, Expressed in Cartesian Coordinates

Conservation of Angular Momentum = "Spin"

Lecture 6. P ω ω ε χ ω ω ω ω E ω E ω (2) χ ω ω ω χ ω ω ω χ ω ω ω (2) (2) (2) (,, ) (,, ) (,, ) (2) (2) (2)

Please initial the statement below to show that you have read it

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Tensor Smooth Length for SPH Modelling of High Speed Impact

Week 8: Chapter 9. Linear Momentum. Newton Law and Momentum. Linear Momentum, cont. Conservation of Linear Momentum. Conservation of Momentum, 2

Physics 207: Lecture 20. Today s Agenda Homework for Monday

PHYSICS 231 Review problems for midterm 2

EN40: Dynamics and Vibrations. Homework 4: Work, Energy and Linear Momentum Due Friday March 1 st

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Invariant deformation parameters from GPS permanent networks using stochastic interpolation

Solid Mechanics Z. Suo

Physics 106a, Caltech 11 October, Lecture 4: Constraints, Virtual Work, etc. Constraints

Navier Stokes Second Exact Transformation

Effect of loading frequency on the settlement of granular layer

Chapter 12 Equilibrium & Elasticity

C PLANE ELASTICITY PROBLEM FORMULATIONS

PHZ 6607 Lecture Notes

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Thermodynamics General

THEOREMS OF QUANTUM MECHANICS

How Differential Equations Arise. Newton s Second Law of Motion

Lecture 5.8 Flux Vector Splitting

FUZZY FINITE ELEMENT METHOD

PHYS 705: Classical Mechanics. Calculus of Variations II

Δ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions

Δ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions

Estimation of homogenized elastic coefficients of pre-impregnated composite materials

Lecture 14: Forces and Stresses

Turbulence. Lecture 21. Non-linear Dynamics. 30 s & 40 s Taylor s work on homogeneous turbulence Kolmogorov.

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Fracture analysis of FRP composites using a meshless finite point collocation method

The classical spin-rotation coupling

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics

Physics 207 Lecture 13. Lecture 13

COMSOL Conference 2008

Short historical notes

MAGNUM - A Fortran Library for the Calculation of Magnetic Configurations

10/23/2003 PHY Lecture 14R 1

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15

SIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

ESS 265 Spring Quarter 2005 Time Series Analysis: Error Analysis

An (almost) unbiased estimator for the S-Gini index

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

1 Matrix representations of canonical matrices

Chapter 3. r r. Position, Velocity, and Acceleration Revisited

Experimental Errors and Error Analysis

Turbulent Transport in Single-Phase Flow. Peter Bernard, University of Maryland

MECHANICS OF MATERIALS

11. Dynamics in Rotating Frames of Reference

A boundary element method with analytical integration for deformation of inhomogeneous elastic materials

3/31/ = 0. φi φi. Use 10 Linear elements to solve the equation. dx x dx THE PETROV-GALERKIN METHOD

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

An introduction to continuum mechanics and elastic wave propagation

Feb 14: Spatial analysis of data fields

Numerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Bending and Vibrations of a Thick Plate with Consideration of Bimoments

Center of Mass and Linear Momentum

Transcription:

Introdcton to elastc wave eqaton Salam Alnabls Unversty of Calgary Department of Mathematcs and Statstcs October 15,01

Otlne Motvaton Elastc wave eqaton Eqaton of moton, Defntons and The lnear Stress- Stran relatonshp The Sesmc Wave Eqaton n Isotropc Meda Sesmc wave eqaton n homogeneos meda Acostc wave eqaton Short Smmary

Motvaton Elastc wave eqaton has been wdely sed to descrbe wave propagaton n an elastc medm, sch as sesmc waves n Earth and ltrasonc waves n hman body. Sesmc waves are waves of energy that travel throgh the earth, and are a reslt of an earthqake, exploson, or a volcano.

Elastc wave eqaton The standard form for sesmc elastc wave eqaton n homogeneos meda s : ρ ( ). ρ :s the densty :s the dsplacement, : Lame parameters

Eqaton of Moton We wll depend on Newton s second law F=ma m : mass ρdx dx 1 a : acceleraton t The total force from stress feld: dx 3 F F F body F F body f x dx dx 1 1 dx dx dx 3 dx 3 dx dx 1 dx 3

Eqaton of Moton Combnng these nformaton together we get the Momentm eqaton (Eqaton of Moton) t f where s thedensty, s and s thestresstensor. the dsplacement,

Defntons Stress : A measre of the nternal forces actng wthn a deformable body. (The force actng on a sold to deform t) The stress at any pont n an obect, assmed to behave as a contnm, s completely defned by nne component stresses: three orthogonal normal stresses and sx orthogonal shear stresses.

11 1 31 1 3 13 3 33 Ths can be expressed as a second-order tensor known as the Cachy stress tensor.

Defntons Stran : A local measre of relatve change n the dsplacement feld, that s, the spatal gradents n the dsplacement feld. And t related to deformaton, or change n shape, of a materal rather than any change n poston. e 1 ( )

Some possble strans for two- dmensonal element

The lnear Stress-Stran Relatonshp Stress and Stran are lnked n elastc meda by Stress - Stran or constttve relatonshp. The most general lnear relatonshp between Stress and Stran s : where, C kl e kl C Stffness(or Elastc coeffcent) C kl s termed the elastc tensor.

The lnear Stress-Stran Relatonshp The elastc tensor C kl, s forth-order wth 81 components ( 1,,k,l 3 ). Becase of the symmetry of the stress and stran tensors and the thermodynamc consderatons, only 1 of these components are ndependent. The 1 components are necessary to specfy the stress-stran relatonshp for the most general form of an elastc sold.

The lnear Stress-Stran Relatonshp The materal s sotropc f the propertes of the sold are the same n all drectons. The materal s ansotropc f the propertes of the meda vary wth drecton.

The lnear Stress-Stran Relatonshp If we assme sotropy, the nmber of the ndependent parameters s redced to two : C ( where and are called the Lame parameters δ kl 0 for : A measre of the resstance of the materal xy e xy kl 1for, δ l : Has no smple physcal explanaton. k k l ) to shearng

The lnear Stress-Stran Relatonshp The stress-stran eqaton for an sotropc meda : [ kl ( l k k l )] e kl e kk e

The lnear Stress-Stran Relatonshp The lnear sotropc stress-stran relatonshp e kk e The stran tensor s defned as : 1 e ( ) (1) () Sbstttng for () n (1) we obtan : k k ( ) (3)

The Sesmc Wave Eqaton n Isotropc Meda Sbstttng (3) n the homogeneos eqaton of moton : ] ) ( [ k k t k k k k ) ( k k k k ) (

The Sesmc Wave Eqaton n Isotropc Meda Defnng we can wrte ths n vector form as t T ρ (. ).[ ( ) ] ( ). se the vector dentty weobtan : ρ (. ).[ ( ) T. ] ( ).

The Sesmc Wave Eqaton n Isotropc Meda Ths s one form of the sesmc wave eqaton ρ (. ).[ ( ) The frst two terms on the (r.h.s) nvolve gradent n the Lame parameters and are nonzero whenever the materal s nhomogeneos (.e. : contans velocty gradent) Incldng these factors makes the eqatons very complcated and dffclt to solve effcently. T ] ( ).

The Sesmc Wave Eqaton n Isotropc Meda If velocty s only a fncton of depth, then the materal can be modeled as a seres of homogeneos layers. Wthn each layer, there are no gradents n the Lames parameters and so these terms go to zero. The standard form for sesmc wave eqaton n homogeneos meda s : ρ ( ). Note : Here we neglected the gravty and velocty gradent terms and has assmed a lnear, sotropc Earth model

Sesmc Wave Eqaton n homogeneos meda If, and are constants,the waveeqaton s smplfed. as: thes- wherethe P - wave velocty wave velocty

Acostc Wave Eqaton If the Lame parameter µ = 0 (.e. No shearng ) then we get : 1 c t where c the speedof propagaton In ths case, the Elastc wave eqaton s redced to an acostc wave eqaton.

Short Smmary We ntrodced defntons of Stress and Stran and the relatonshp between them. We depend on Newton s nd law to get the eqaton of moton and from t we Derve the general form of Elastc wave eqaton. We smplfy t to the standard form by modelng the materal as seres of homogeneos layers. We dscssed two types of waves P-waves(Compressonal) S-waves(Shear) Fnally, f we assme no shearng then we redced t to an acostc wave eqaton.