



Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Community
Ask the community for help and clear up your study doubts
Discover the best universities in your country according to Docsity users
Free resources
Download our free guides on studying techniques, anxiety management strategies, and thesis advice from Docsity tutors
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
1 / 6
This page cannot be seen from the preview
Don't miss anything!
∂ ui^2
∂ t^2
Generalize Hooke’s Law / Relationship between stress tensor & strain tensor
ui
∂ uk
∂ t ∂ xl ↑ ↑ ↑ f i =0 homogeneous isotropic
Using ∇ 2 u = ∇( ∇ ⋅ u ) − ∇×∇× u the above can be rewritten in a vector form:
(3) u =
& & (^) ⎜⎜ (^) ⎟⎟∇( ∇ ⋅ u ) −
∇ ×( ∇ × u )
change, equation (3) can be rewritten as:
Our equations for P -wave and S -wave velocities:
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:
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:
u = − ∂ z x , uy =^ , u^ z =
u
u
u
∂ z ∂ z x Æ^ pure shear y Æ^ pure shear z Æ^ pure^ P -wave
Recall: k
α =
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)
Even though it is not physical wave elastic theory can derive it.
Poisson’s medium ( λ=μ )
“modern papers”
Solutions to the Wave Equation
1. d’Alembert’s solution f(± x ± ct),
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 ( x − ct )) = A (cos( x − ct ) + i sin( x − ct ))
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.
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:
The vector potential for S-wave satisfying equation (7) is: