31.08.2018 Views

Eduardo Kausel-Fundamental solutions in elastodynamics_ a compendium-Cambridge University Press (2006)

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

1.3 Coord<strong>in</strong>ate systems and differential operators 9<br />

Laplacian<br />

∂<br />

∂r + 1 ∂ 2<br />

r 2 ∂θ + ∂2<br />

2 ∂z 2<br />

∇ 2 =∇· ∇= ∂2<br />

∂r + 1 2 r<br />

[ ∂<br />

∇ 2 2 u r<br />

u =∇· ∇u = ˆr<br />

∂r + 1 ∂u r<br />

2 r<br />

[ ∂<br />

+ ˆt<br />

2 u θ<br />

∂r + 1 ∂u θ<br />

2 r<br />

∂r + 1 ( ∂ 2 u r<br />

r 2 ∂θ 2<br />

∂r + 1 ( ∂ 2 u θ<br />

r 2 ∂θ 2<br />

[ ∂<br />

+ ˆk<br />

2 u z<br />

∂r + 1 ∂u z<br />

2 r ∂r + 1<br />

∂ 2 u z<br />

r 2 ∂θ 2<br />

− u r − 2 ∂u ) ]<br />

θ<br />

+ ∂2 u r<br />

∂θ ∂z 2<br />

− u θ + 2 ∂u r<br />

∂θ<br />

]<br />

+<br />

∂2 u z<br />

∂z 2<br />

) ]<br />

+ ∂2 u θ<br />

∂z 2<br />

(1.27)<br />

(Note: ∂ 2 /∂θ 2 <strong>in</strong> ∇ 2 acts both on the components of u and on the basis vectors ˆr, ˆt.)<br />

Waveequation<br />

(λ + µ)∇∇ · u + µ∇ · ∇u + b = ρü (1.28)<br />

Expansion of avector <strong>in</strong> Fourier series<strong>in</strong>the azimuth<br />

∞∑<br />

( ) cos nθ<br />

u = u n (1.29a)<br />

s<strong>in</strong> nθ<br />

n=0<br />

∞∑<br />

( ) −s<strong>in</strong> nθ<br />

v = v n (1.29b)<br />

cos nθ<br />

n=0<br />

∞∑<br />

( ) cos nθ<br />

w = w n (1.29c)<br />

s<strong>in</strong> nθ<br />

n=0<br />

<strong>in</strong> which u ≡ u r , v ≡ u θ , w ≡ u z , and either the lower or the upper element <strong>in</strong> the parentheses<br />

must be used, as may be necessary. Also, u n , v n , w n are the coefficients of the Fourier<br />

series, which do not depend on θ, but only on r and z, that is, u n = u n (r, z), and so forth.<br />

1.3.3 Spherical coord<strong>in</strong>ates<br />

Source–receiver distance R = √ x 2 + y 2 + z 2 (1.30a)<br />

Range r = √ x 2 + y 2 (1.30b)<br />

Azimuth tan θ = y/x, 0 ≤ θ ≤ 2π (1.30c)<br />

Polar angle φ = arccos (z/R), 0 ≤ φ ≤ π (1.30d)<br />

Direction cos<strong>in</strong>es γ 1 = s<strong>in</strong> φ cos θ ⇒ cos θ =<br />

γ 1<br />

√<br />

1 − γ 2 3<br />

(1.30e)<br />

γ 2 = s<strong>in</strong> φ s<strong>in</strong> θ ⇒ s<strong>in</strong> θ = √<br />

(1.30f)<br />

1 − γ3<br />

2<br />

√<br />

γ 3 = cos φ ⇒ s<strong>in</strong> φ = 1 − γ3 2 (1.30g)<br />

γ 2

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!