Docsity
Docsity

Prepare for your exams
Prepare for your exams

Study with the several resources on Docsity


Earn points to download
Earn points to download

Earn points by helping other students or get them with a premium plan


Guidelines and tips
Guidelines and tips

Generalize Hooke’s Law - Seismology - Lecture Notes, Study notes of Geology

These Lecture Notes cover the following aspects of Seismology : Generalize Hooke’S Law, Relationship, Hooke’S Law, Stress Tensor, Strain Tensor, Homogeneous, Isotropic, Rewritten, Vector Form, Shear

Typology: Study notes

2012/2013

Uploaded on 07/19/2013

shantii
shantii 🇮🇳

4.4

(14)

98 documents

1 / 6

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
(1) Σ F = am i
ρ
ui
2
=
σ
+ fi
i ij j
t
2
Generalize Hooke’s Law / Relationship between stress tensor & strain tensor
2
(2)
ρ
ui = j (cijklekl ) = cijkl
j(e ) = (
λ
+
µ
) uk +
µ
2uikl i
2
t l
x
fi=0 homogeneous isotropic
media media (i.e., )( jkiljlikklij
cijkl
δ
δ
µ
δ
δ
δ
λδ
+
+
=
)
Using
2u = ( u) ×× u the above can be rewritten in a vector form:
(3) u =
λ
+ 2
µ
&&
( u)
µ
× ( ×
u)
ρ
ρ
We set
θ
= u , which is the cubic dilatation; and ω
=
×
u , which is shear/rotation
change, equation (3) can be rewritten as:
2
&& 2
(4) u =
α
θ
β
× ω
Our equations for P-wave and S-wave velocities:
k + 4/ 3
µ
(5a)
α
=
λ
+ 2
µ
, or with k =
λ
+
2 / 3
µ
, we can get
α
=
ρρ
(5b)
β
=
µ
ρ
N.B. Equations (3) and (4) begin to look like wave equations, but we refer to them as
equations of motion. Solving them for
u is not easy, but can be done numerically.
We first show that we can obtain wave equations for P and S waves. Then the solution
is used to go back and solve the equation of motion. There are two ways of doing
thi
1. (This is not done in Stein and Wysession’s book) Isolate the part that only results in
volume change and the part that only results in shear.
To isolate the volume component take the divergence:
2
2(
2
α
θ
β
× ω) = ( u) =
θ
t
2 t
2
Because × ω)
=
0 , this becomes:
(
2 2 &&
(6)
α
θ
=
θ
1
2
pf3
pf4
pf5

Partial preview of the text

Download Generalize Hooke’s Law - Seismology - Lecture Notes and more Study notes Geology in PDF only on Docsity!

(1) Σ F = m ai → ρ

ui^2

i ⋅^ =σ^ ij ∂^ j + fi

t^2

Generalize Hooke’s Law / Relationship between stress tensor & strain tensor

ui

= ∂ j ( cijklekl ) = cijkl ∂ j ( e ) = ( λ +μ )∂

uk

2 kl^ i +μ∇^2 ui

txl ↑ ↑ ↑ f i =0 homogeneous isotropic

media media ( i.e. , c ijkl = λδ ij δ kl + μ( δ ik δ jl + δ il δ jk ))

Using ∇ 2 u = ∇( ∇ ⋅ u ) − ∇×∇× u the above can be rewritten in a vector form:

(3) u =

& & (^) ⎜⎜ (^) ⎟⎟∇( ∇ ⋅ u ) −

∇ ×( ∇ × u )

⎝^ ρ^ ⎠ ⎝

We set θ = ∇ ⋅ u , which is the cubic dilatation; and ω = ∇ × u , which is shear/rotation

change, equation (3) can be rewritten as:

(4) u & & =α ∇^2 θ −β 2 ∇ × ω

Our equations for P -wave and S -wave velocities:

k + 4 / 3 μ

(5a) α =

, or with k = λ + 2 / 3 μ , we can get α =

(5b) β =

N.B. Equations (3) and (4) begin to look like wave equations, but we refer to them as equations of motion. Solving them for u is not easy, but can be done numerically. We first show that we can obtain wave equations for P and S waves. Then the solution is used to go back and solve the equation of motion. There are two ways of doing thi

1. (This is not done in Stein and Wysession’s book) Isolate the part that only results in volume change and the part that only results in shear. To isolate the volume component take the divergence:

2 ∂^2 (

α ∇⋅∇ θ − β ∇ ⋅∇× ω ) = ( ∇ ⋅ u ) = θ ∂ t^2 ∂ t^2 Because ∇ ⋅(∇× ω ) = 0 , this becomes:

(6) α ∇^2 2 θ =^ θ&&

2

To separate the shear part, take the curl:

∂ 2 α 2 ∇×∇ θ −β 2 ∇ ×( ∇ × ω ) =

( ∇ × u ) = ωt^2 ∂ t^2 Because ∇ ×(∇ θ ) = 0 , this becomes:

(7) β ∇^2 2 ω = ω &&

This is a simple way of isolating change in volume and change in rotation. A more elegant/powerful method is our second method:

2. Helmholtz decompositic of a vector field (☼) (8) ☼= ∇ A + ∇× B , with ∇ ⋅ B = 0 B is the vector potential, and A is the scalar potential. There are 4 components: 1 gradient and 3 curl. NEED the constant that ∇ ⋅ B = 0.

Examples: Gravity g = −∇ U (^) grav rotation free because gravity is a conservative force field Magnetic Field B = −∇ Vmag + ∇ × f e.g. Lorentz force → Maxwell’s equations

Here we can write for displacement:

(9) u = ∇ φ + ∇× ψ , with ∇ ⋅ ψ = 0

This is the Helmholtz representation with displacement potentials φ and ψ.

This can be substituted into the equations of motion to get:

⎫ (10) ∇⎨⎜⎜ &&^

⎨⎜⎜^

⎟⎟∇^

+^ ∇×

⎟⎟∇^ ψ^ −^ ψ ⎬ =^0

⎩⎝^ ρ^ ⎠^ ⎭ ⎩⎝^ ρ^ ⎠^ ⎭

One way to solve this equation is to say that both parts are zero at the same time.

The P -wave potential equation:

(11) φ&^ &^ =^ λ^ +^2 μ^ ∇ 2 φ ρ The S -wave potential equation:

(12) ψ &^ & =

μ ⎞ ⎟⎟∇^

(^2) ψ ⎝^ ρ^ ⎠ These are not the only solutions to (10). There could be a coupling between the P -wave and S -wave parts at low frequencies.

Propagation direction

In 1D, just vertical propagation:

∂ψ y ∂ψ x ∂φ

u = − ∂ z x , uy =^ , u^ z =

u

u

u

zz x Æ^ pure shear y Æ^ pure shear z Æ^ pure^ P -wave

Recall: k

α =

λ + 2 μ

k + 4 / 3 μ ρ ρ Z (5a, b)

β =

μ ρ

where α is the P -wave velocity and β is the S -wave velocity. The bulk modulus ( k ) can be measured in the lab, hence, it is sometimes useful to write the velocities in terms of this constant. The acoustic sound wave speed can be written as:

k (15)

Bulk sound speed: V φ^2 =α 2 −^4 β 2 = k

Even though it is not physical wave elastic theory can derive it.

Poisson’s medium ( λ=μ )

α = 3 β This will work for most areas of the mantle…this is not used generally in

“modern papers”

Solutions to the Wave Equation

φ =&& α ∇^2 2 φ

Ψ&&^2

SV =β ∇^ Ψ SV

Ψ&&^2

SH =β ∇^ Ψ SH

1. d’Alembert’s solution f(± x ± ct),

φ ( x , t ) = f ( x − ct ) + g ( x + ct ) , f and g are arbitrary functions as long as they are at least

twice differentiable at time and space. (x-ct) propagation in the positive x-direction (x+ct) propagation in the negative x-direction c is the wave speed or phase velocity

c = x 1 − x 0 looks at how long the peak propagates through the system t 1 − t 0

f ( ) f^ ( )

Example of a plane wave propagation in positive x-axis; function f ( ) remains same if argument (the phase!) remains same; this happens if x increases when t increases → motion in positive x axis.

Plane Waves: Surfaces of constant phase change Æ wave front

Wave front ≡ surface contains points of equal phase.

Ray k ˆ

One example of a general solution is the harmonic function: A exp( i ( xct )) = A (cos( xct ) + i sin( xct ))

x − ct = phase =ϕ

so ... Aei^ ϕ

The real part of the solution gives the displacement; the complex part can be used to describe amplitude decay; makes the notation elegant.

Angular frequency: ω = kc , dispersion relationship k = ω / c

Wave length: λ = cT = 2

Wave number: k =^2 π^ λ, high wave number corresponding to short wave length

k

ψ=(kx-ωt) a.k.a. Fourier duals will transform from time to frequency domain in 3D ( k · x - ωt)

P -wave and S -wave potential The scalar potential for a harmonic plane P -wave satisfying motions of equation (6) is:

φ ( x , t ) = A exp( i ( k ⋅ x − ω t ))

The vector potential for S-wave satisfying equation (7) is: