The conservation of
linear momentum of a
linear elastic material can be written as: :\boldsymbol{\nabla}\cdot(\mathsf{C}:\boldsymbol{\nabla}\mathbf{u}) = \rho~\ddot{\mathbf{u}} where \mathbf{u} is the
displacement vector and \mathsf{C} is the
stiffness tensor. Love waves are a special solution (\mathbf{u}) that satisfy this system of equations. We typically use a Cartesian coordinate system (x,y,z) to describe Love waves. Consider an isotropic linear elastic medium in which the elastic properties are functions of only the z coordinate, i.e., the
Lamé parameters and the mass
density can be expressed as \lambda(z), \mu(z), \rho(z). Displacements (u,v,w) produced by Love waves as a function of time (t) have the form : u(x,y,z,t) = 0 ~,~~ v(x,y,z,t) = \hat{v}(x,z,t) ~,~~ w(x,y,z,t) = 0 \,. These are therefore
antiplane shear waves perpendicular to the (x,z) plane. The function \hat{v}(x,z,t) can be expressed as the superposition of
harmonic waves with varying
wave numbers (k) and
frequencies (\omega). Consider a single harmonic wave, i.e., : \hat{v}(x,z,t) = V(k, z, \omega)\,\exp[i(k x - \omega t)] where i is the
imaginary unit, i.e. i^2 = -1. The
stresses caused by these displacements are : \sigma_{xx} = 0 ~,~~ \sigma_{yy} = 0 ~,~~ \sigma_{zz} = 0 ~, ~~ \tau_{zx} = 0 ~,~~ \tau_{yz} = \mu(z)\,\frac{dV}{dz}\,\exp[i(k x - \omega t)] ~,~~ \tau_{xy} = i k \mu(z) V(k, z, \omega) \,\exp[i(k x - \omega t)] \,. If we substitute the assumed displacements into the equations for the conservation of momentum, we get a simplified equation : \frac{d}{dz}\left[\mu(z)\,\frac{dV}{dz}\right] = [k^2\,\mu(z) - \omega^2\,\rho(z)]\,V(k,z,\omega) \,. The boundary conditions for a Love wave are that the
surface tractions at the free surface (z = 0) must be zero. Another requirement is that the stress component \tau_{yz} in a layer medium must be continuous at the interfaces of the layers. To convert the second order
differential equation in V into two first order equations, we express this stress component in the form : \tau_{yz} = T(k, z, \omega)\,\exp[i(k x - \omega t)] to get the first order conservation of momentum equations : \frac{d}{dz}\begin{bmatrix} V \\ T \end{bmatrix} = \begin{bmatrix} 0 & 1/\mu(z) \\ k^2\,\mu(z) - \omega^2\,\rho(z) & 0 \end{bmatrix} \begin{bmatrix} V \\ T \end{bmatrix} \,. The above equations describe an
eigenvalue problem whose solution
eigenfunctions can be found by a number of
numerical methods. Another common, and powerful, approach is the
propagator matrix method (also called the
matricant approach). ==See also==