- Art Gallery -

# .

In classical mechanics, moment of inertia, also called mass moment of inertia, rotational inertia, polar moment of inertia of mass, or the angular mass, (SI units kg·m²) is a measure of an object's resistance to changes to its rotation. It is the inertia of a rotating body with respect to its rotation. The moment of inertia plays much the same role in rotational dynamics as mass does in linear dynamics, describing the relationship between angular momentum and angular velocity, torque and angular acceleration, and several other quantities. The symbols I and sometimes J are usually used to refer to the moment of inertia or polar moment of inertia.

While a simple scalar treatment of the moment of inertia suffices for many situations, a more advanced tensor treatment allows the analysis of such complicated systems as spinning tops and gyroscopic motion.

The concept was introduced by Leonhard Euler in his book Theoria motus corporum solidorum seu rigidorum in 1765.[1] In this book, he discussed the moment of inertia and many related concepts, such as the principal axis of inertia.

Overview
Spinning figure skater pulling in her arms, reducing her moment of inertia, to rotate faster

The moment of inertia of an object about a given axis describes how difficult it is to change its angular motion about that axis. Therefore, it encompasses not just how much mass the object has overall, but how far each bit of mass is from the axis. The further out the object's mass is, the more rotational inertia the object has, and the more rotational force (torque, the force multiplied by its distance from the axis of rotation) is required to change its rotation rate.

For example, consider two wheels suspended so they can turn freely, a large bicycle wheel and a small baby buggy wheel. Suppose they have the same weight, which means they have equal mass. It requires more effort (torque) to spin (accelerate) the bicycle wheel to a given angular velocity than the baby buggy wheel. This is because the rim of the bicycle wheel is further from its axis than the rim of the baby buggy wheel. Even though they have the same amount of mass, most of the mass of the bicycle wheel is located further from the axis than the mass in the baby buggy wheel, so it must move faster for a given rotation rate. So the bicycle wheel has a larger moment of inertia than the baby buggy wheel.

The moment of inertia of an object can change if its shape changes. Figure skaters who begin a spin with arms outstretched provide a striking example. By pulling in their arms, they reduce their moment of inertia, causing them to spin faster by the conservation of angular momentum.

The moment of inertia has two forms, a scalar form, I, (used when the axis of rotation is specified) and a more general tensor form that does not require the axis of rotation to be specified. The scalar moment of inertia, I, (often called simply the "moment of inertia") allows a succinct analysis of many simple problems in rotational dynamics, such as objects rolling down inclines and the behavior of pulleys. For instance, while a block of any shape will slide down a frictionless decline at the same rate, rolling objects may descend at different rates, depending on their moments of inertia. A hoop will descend more slowly than a solid disk of equal mass and radius because more of its mass is located far from the axis of rotation. However, for (more complicated) problems in which the axis of rotation can change, the scalar treatment is inadequate, and the tensor treatment must be used (although shortcuts are possible in special situations). Examples requiring such a treatment include gyroscopes, tops, and even satellites, all objects whose alignment can change.
A flywheel is a wheel with a large moment of inertia used to smooth out motion in machines. This example is in a Russian factory.

The moment of inertia is also called the mass moment of inertia (especially by mechanical engineers) to avoid confusion with the second moment of area, which is sometimes called the area moment of inertia (especially by structural engineers). The easiest way to differentiate these quantities is through their units (kg·m² as opposed to m4). In addition, moment of inertia should not be confused with polar moment of inertia (more specifically, polar moment of inertia of area), which is a measure of an object's ability to resist torsion (twisting) only, although, mathematically, they are similar: if the solid for which the moment of inertia is being calculated has uniform thickness in the direction of the rotating axis, and also has uniform mass density, the difference between the two types of moments of inertia is a factor of mass per unit area.
Scalar moment of inertia for single body

Consider a massless rigid rod of length l with a point mass m at one end and rotating about the other end. Suppose the rod rotates at a constant rate so that the mass moves at speed v. Then the kinetic energy T of the mass is:

$$T = \frac12\,mv^2$$

Using v = lω, where ω is the angular velocity, one obtains:

$$T = \frac12\,m (\omega l)^2$$

which can be rearranged to give:

$$T = \frac12\, ml^2 \omega^2$$

This equation resembles the original expression for the kinetic energy, but in the place of the linear velocity v is the angular velocity ω, and instead of the mass m is ml2. The quantity ml2 can therefore be seen as an analogue of mass for rotational motion; in other words, it is a measure of rotational inertia.

Scalar moment of inertia for many bodies

Consider a rigid body rotating with angular velocity ω around a certain axis. The body consists of N point masses mi whose distances to the axis of rotation are denoted ri. Each point mass will have the speed vi = ωri, so that the total kinetic energy T of the body can be calculated as

$$T = \sum_{i=1}^N \frac12\,m_i v_i^2 = \sum_{i=1}^N \frac12\,m_i (\omega r_i)^2 = \frac12\, \omega^2 \Big( \textstyle \sum_{i=1}^N m_i r_i^2 \Big).$$

In this expression the quantity in parentheses is called the moment of inertia of the body (with respect to the specified axis of rotation). It is a purely geometric characteristic of the object, as it depends only on its shape and the position of the rotation axis. The moment of inertia is usually denoted with the capital letter I:

$$I = \sum_{i=1}^N m_i r_i^2\ .$$

It is worth emphasizing that ri here is the distance from a point to the axis of rotation, not to the origin. As such, the moment of inertia will be different when considering rotations about different axes.

Similarly, the moment of inertia of a continuous solid body rotating about a known axis can be calculated by replacing the summation with the integral:

$$I = \int_V \rho(\mathbf{r})\,d(\mathbf{r})^2 \, \mathrm{d}V\!(\mathbf{r}),$$

where r is the radius vector of a point within the body, ρ(r) is the mass density at point r, and d(r) is the distance from point r to the axis of rotation. The integration is evaluated over the volume V of the body.
Moment of inertia about a point

In discussions of the virial theorem, a different definition of moment of inertia is introduced, often qualified as moment of inertia about a point (as opposed to an axis). The defining equations look the same as the ones above, but $$r_i$$ and $$d(\mathbf{r})$$ in these cases are understood to be the distance to the origin.
Moment of inertia theorems

Calculations for the moment of inertia (MOI) of a body are not easy in general. The process can be simplified in the following ways:

Choosing axes to take advantage of geometric symmetries
Physical homeogeneity (i.e. uniform mass distribution) making the density function ρ(r) become constant (elementary calculations, or generally an approximation)
Use of the following theorems[2], see below for details.

Theorem Nomenclature Equation
Superposition Principle for
Moment Of Inertia (MOI) about any chosen Axis
Inet = Resultant MOI (about any one axis) $$I_\mathrm{net} = \sum_j I_j\,\!$$
Parallel axis theorem M = Total mass of body

d = Perpendicular distance from an axis
through the Center Of Mass (COM) to another parallel axis
ICOM = MOI about the axis through
the COM
Id = MOI about the parallel axis

$$I_d = I_\mathrm{com} + Md^2 \,\!$$
Perpendicular axis theorem i, j, k refer to MOI about any three mutually

perpendicular axes.

The sum of MOI about any two is greater than
or equal to the third.

$$I_k \leq I_i + I_j \,\!$$

Properties

Superposition

The moment of inertia of the body is additive. That is, if a body can be decomposed (either physically or conceptually) into several constituent parts, then the moment of inertia of the whole body about a given axis is equal to the sum of moments of inertia of each part around the same axis.[3]

Basing just on the dimensional analysis, the moment of inertia must take the form I = c·M·L2, where M is the mass, L is the “size” of the body in the direction perpendicular to the axis of rotation, and c is a dimensionless inertial constant. Additionally, the length $$\scriptstyle R\,=\,\sqrt{I/M}$$ is called the radius of gyration of the body.

Perpendicular Axes

If Ix, Iy, Iz are moments of inertia around three perpendicular axes passing through the body’s center of mass, then each of them cannot be greater than the sum of two others: for example IzIx + Iy. Here the equality holds only if the body is flat and located in the Oxy coordinate plane.

Parallel Axes

If the object’s moment of inertia ICOM around a certain axis passing through the center of mass is known, then the parallel axis theorem or Huygens–Steiner theorem provides a convenient formula to compute the moment of inertia Id of the same body around a different axis, which is parallel to the original and located at a distance d from it. The formula is only suitable when the initial and final axes are parallel. In order to compute the moment of inertia about an arbitrary axis, one has to use the object’s moment of inertia tensor.
Energy, angular momentum, torque

The rotational kinetic energy of a rigid body with angular velocity ω (in radians per second) is expressed in terms of the object’s moment of inertia:

$$E_r = \frac12 \, I \omega^2 .$$

This formula is similar to the translational kinetic energy Et = ½mv 2. Thus, the moment of inertia I plays the role of mass in rotational dynamics. One key difference between the mass and the (scalar) moment of inertia is that the latter depends on the axis of rotation and so is not truly invariant. The invariant characteristic of the body in rotational motion is the tensor of moment of inertia I, defined later.

The angular momentum of the body rotating around one of its principal axes is also proportional to the moment of inertia:

$$L = I\omega\, .$$

This expression is parallel to the formula for translational momentum p = mv, where the moment of inertia I plays the role of mass m, and the angular velocity ω stands for the velocity v. The scalar formula is valid for rotations around one of the principal axes of the body or for rotation about any fixed axis. The equivalent formula involving the tensor moment of inertia is always correct and must be used in cases of free rotation about a non-principal axis.

Also when the body rotates around one of its principal axes, and the direction of the axis of rotation remains constant, one can relate the torque on an object and its angular acceleration in a similar equation:

$$\tau = I\alpha \,,$$

where τ is the torque and α is the angular acceleration.
Examples
Main article: List of moments of inertia
Sodium-hydride-monomer-CRC-MW-3D-balls.png

Diatomic molecule, with atoms m1 and m2 at a distance d from each other, rotating around the axis which passes through the molecule’s center of mass and is perpendicular to the direction of the molecule.

The easiest way to calculate this molecule’s moment of inertia is to use the parallel axis theorem. If we consider rotation around the axis passing through the atom m1, then the moment of inertia will be I1 = m1·0 + m2·d 2 = m2d 2. On the other hand, by the parallel axis theorem this moment is equal to I1 = I + (m1 + m2a 2, where I is the moment of inertia around the axis passing through the center of mass, and a is the distance between the center of mass and the first atom. By the center of mass formula, this distance is equal to a = m2d / (m1 + m2). Thus,

$$I = I_1 - (m_1+m_2)a^2 = m_2d^2 - \frac{m_2^2d^2}{m_1+m_2} = \frac{m_1m_2}{m_1+m_2}\,d^2 .$$

Moment of inertia rod center.png

Thin rod of mass m and length , rotating around the axis which passes through its center and is perpendicular to the rod.

Let Oz be the axis of rotation, and Ox the axis along the rod. If ρ is the density, and s the cross-section of the rod (so that m = ρℓs), then the volume element for the integral formula will be equal to dV = s·dx, where x changes from −½ to ½. The moment of inertia can be found by computing the integral:

$$I = \int_{-\ell/2}^{\ell/2} \rho\,x^2 s \mathrm{d}x = \rho s \,\frac{x^3}{3}\bigg|_{-\ell/2}^{\ell/2} = \frac{m}{s\ell} \cdot s \cdot 2\frac{\ell^3/8}{3} = \frac{1}{12}\, m\ell^2 .$$

Solid ball of mass m and radius R, rotating around an axis which passes through the center.

Suppose Oz is the axis of rotation. The distance from point r = (x, y, z) towards the axis Oz is equal to d(r)2 = x 2 + y 2. Thus, in order to compute the moment of inertia Iz, we need to evaluate the integral ∭(x 2 + y 2) dV. The calculation considerably simplifies if we notice that by symmetry of the problem, the moments of inertia around all axes are equal: Ix = Iy = Iz. Then

$$I = \frac13(I_x + I_y + I_z) = \frac13 \iiint \rho\cdot(y^2+z^2 + x^2+z^2 + x^2+y^2)\,\mathrm{d}V = \frac23\,\rho \int r^2 \,\mathrm{d}V ,$$

where r2 = x 2 + y 2 + z2 is the distance from point r to the origin. This integral is easy to evaluate in the spherical coordinates, the volume element will be equal to dV = 4πr 2 dr, where r goes from 0 to R. Thus,

$$I = \frac23\,\rho \int_0^R \!\!4\pi r^4\,\mathrm{d}r = \frac23\,\rho\cdot4\pi\frac{R^5}{5} = \frac{m}{\frac43\pi R^3}\cdot\frac{8\pi R^5}{15} = \frac25\,mR^2 .$$

In mechanics of machines, when designing rotary parts like gears, pulleys, shafts, couplings etc., which are used to transmit torques, the moment of inertia has to be considered. The moment of inertia is given about an axis and it depends on the shape, density of a rotating element.

When considering mechanisms like gear trains, worm and wheel, where there are more than one rotating element, more than one axis of rotation, an equivalent moment of inertia for the system should be found. Practically when a geared system is enclosed, equivalent moment of inertia can be measured by measuring the angular acceleration for a known torque or theoretically it can be estimated when the masses and dimensions of the rotating elements and shafts are known. In this practical the equivalent moment of inertia of a worm and wheel system is measured using above mention methods.
Moment of inertia tensor

In three dimensions, if the axis of rotation is not given, we need to be able to generalize the scalar moment of inertia to a quantity that allows us to compute a moment of inertia about arbitrary axes. This quantity is known as the moment of inertia tensor and can be represented as a symmetric positive semi-definite matrix, I. This representation elegantly generalizes the scalar case: The angular momentum vector is related to the rotation velocity vector ω by

$$L=\mathbf{I}\omega$$

and the kinetic energy is given by

$$\frac{1}{2} \omega^\top \mathbf{I} \omega$$

as compared with

$$\frac{1}{2} I\omega^2$$

in the scalar case.

Like the scalar moment of inertia, the moment of inertia tensor may be calculated with respect to any point in space, but for practical purposes, the center of mass is almost always used. In general, its components are time dependent.
Definition

For a rigid object of N point masses m_{k}, the moment of inertia tensor (with respect to the origin) has components given by

\mathbf{I} = \begin{bmatrix} I_{11} & I_{12} & I_{13} \\ I_{21} & I_{22} & I_{23} \\ I_{31} & I_{32} & I_{33} \end{bmatrix} ,\)

where

$$I_{11} = I_{xx} \ \stackrel{\mathrm{def}}{=}\ \sum_{k=1}^{N} m_{k} (y_{k}^{2}+z_{k}^{2}),\,\!$$
$$I_{22} = I_{yy} \ \stackrel{\mathrm{def}}{=}\ \sum_{k=1}^{N} m_{k} (x_{k}^{2}+z_{k}^{2}),\,\!$$
$$I_{33} = I_{zz} \ \stackrel{\mathrm{def}}{=}\ \sum_{k=1}^{N} m_{k} (x_{k}^{2}+y_{k}^{2}),\,\!$$
$$I_{12} = I_{xy} \ \stackrel{\mathrm{def}}{=}\ -\sum_{k=1}^{N} m_{k} x_{k} y_{k},\,\!$$
$$I_{13} = I_{xz} \ \stackrel{\mathrm{def}}{=}\ -\sum_{k=1}^{N} m_{k} x_{k} z_{k},\,\!$$
$$I_{23} = I_{yz} \ \stackrel{\mathrm{def}}{=}\ -\sum_{k=1}^{N} m_{k} y_{k} z_{k},\,\!$$

and $$I_{12}=I_{21}, I_{13}=I_{31}, and I_{23}=I_{32}$$ . (Thus I is a symmetric tensor.) Note that the scalars I_{ij} with i\ne j are called the products of inertia.

Here $$I_{xx}$$ denotes the moment of inertia around the x-axis when the objects are rotated around the x-axis, $$I_{xy}$$ denotes the moment of inertia around the y-axis when the objects are rotated around the x-axis, and so on.

These quantities can be generalized to an object with distributed mass, described by a mass density function, in a similar fashion to the scalar moment of inertia. One then has

$$\mathbf{I}=\iiint_V \rho(x,y,z)\left( (x^2+y^2+z^2) \mathbf{E}_{3} - \begin{bmatrix}xx & xy & xz \\ yx & yy & yz \\ zx & zy & zz\end{bmatrix}\right)\, dx\,dy\,dz$$

In more-concise notation:

$$\mathbf{I}=\int_V \rho(\mathbf{r}) \left( \left( \mathbf{r} \cdot \mathbf{r} \right) \mathbf{E}_{3} - \mathbf{r}\otimes \mathbf{r}\right)\, dV$$

where $$\mathbf{r}$$ is the vector from the center of mass to a point in the volume, V and $$\mathbf{r}\otimes \mathbf{r}$$ is their outer product, E3 is the 3 × 3 identity matrix, and V is a region of space completely containing the object. Alternatively, the above can be described in Einstein notation as:

{{I}_{ij}} = \int_{V}{\rho(r) \left( r_l^2 \delta_{ij} - r_i r_j \right)\,dV} \)

where $$\delta_{ij}$$ is the Kronecker delta.

The diagonal elements of I are called the principal moments of inertia.

As is clear from the definition, I and its components are in general time-dependent. They are constant only if calculated in a reference frame moving rigidly with the rigid body. On the other hand, the eigenvalues I are invariant, that is independent of the reference frame used to calculate I.
Derivation of the tensor components

The distance r of a particle at $$\mathbf{x}$$ from the axis of rotation passing through the origin in the $$\mathbf{\hat{n}}$$ direction is $$|\mathbf{x}-(\mathbf{x} \cdot \mathbf{\hat{n}}) \mathbf{\hat{n}}|$$ . By using the formula $$I=mr^2$$ (and some simple vector algebra) it can be seen that the moment of inertia of this particle (about the axis of rotation passing through the origin in the \mathbf{\hat{n}} direction) is $$I=m(|\mathbf{x}|^2 (\mathbf{\hat{n}} \cdot \mathbf{\hat{n}})-(\mathbf{x} \cdot \mathbf{\hat{n}})^2)$$ This is a quadratic form in $$\mathbf{\hat{n}}$$ and, after a bit more algebra, this leads to a tensor formula for the moment of inertia

$${I} = m [n_1\,n_2\,n_3]\begin{bmatrix} y^2+z^2 & -xy & -xz \\ -y x & x^2+z^2 & -yz \\ -zx & -zy & x^2+y^2 \end{bmatrix} \begin{bmatrix} n_1 \\ n_2\\ n_3 \end{bmatrix} . This is exactly the formula given below for the moment of inertia in the case of a single particle. For multiple particles we need only recall that the moment of inertia is additive in order to see that this formula is correct. Reduction to scalar For any axis \( \hat{\mathrm{n}}$$ , represented as a column vector with elements ni, the scalar form I can be calculated from the tensor form I as

$$I = \mathbf{\hat{n}^\top} \mathbf{I}\, \mathbf{\hat{n}} = \sum_{j=1}^{3} \sum_{k=1}^{3} n_{j} I_{jk} n_{k}$$

The range of both summations correspond to the three Cartesian coordinates.

The following equivalent expression avoids the use of transposed vectors which are not supported in maths libraries because internally vectors and their transpose are stored as the same linear array,

$$I = \mathbf{{I}^\top} \mathbf{\hat{n}} \cdot \mathbf{\hat{n}}$$

However it should be noted that although this equation is mathematically equivalent to the equation above for any matrix, inertia tensors are symmetrical. This means that it can be further simplified to:

$$I = \mathbf{{I}} \mathbf{\hat{n}} \cdot \mathbf{\hat{n}}$$

Principal axes of inertia

By the spectral theorem, since the moment of inertia tensor is real and symmetric, there exists a Cartesian coordinate system in which it is diagonal, having the form

$$\mathbf{I} = \begin{bmatrix} I_{1} & 0 & 0 \\ 0 & I_{2} & 0 \\ 0 & 0 & I_{3} \end{bmatrix}$$

where the coordinate axes are called the principal axes and the constants $$I_{1}, I_{2}$$ and $$I_{3}$$ are called the principal moments of inertia.

This result was first shown by J. J. Sylvester (1852), and is a form of Sylvester's law of inertia. The principal axis with the highest moment of inertia is sometimes called the figure axis or axis of figure.

When all principal moments of inertia are distinct, the principal axes through center of mass are uniquely specified. If two principal moments are the same, the rigid body is called a symmetrical top and there is no unique choice for the two corresponding principal axes. If all three principal moments are the same, the rigid body is called a spherical top (although it need not be spherical) and any axis can be considered a principal axis, meaning that the moment of inertia is the same about any axis.

The principal axes are often aligned with the object's symmetry axes. If a rigid body has an axis of symmetry of order m, meaning it is symmetrical under rotations of 360°/m about the given axis, that axis is a principal axis. When m>2, the rigid body is a symmetrical top. If a rigid body has at least two symmetry axes that are not parallel or perpendicular to each other, it is a spherical top, for example, a cube or any other Platonic solid.

The motion of vehicles is often described about these axes with the rotations called yaw, pitch, and roll.

A practical example of this mathematical phenomenon is the routine automotive task of balancing a tire, which basically means adjusting the distribution of mass of a car wheel such that its principal axis of inertia is aligned with the axle so the wheel does not wobble.
Parallel axis theorem
Main article: Parallel axis theorem

Once the moment of inertia tensor has been calculated for rotations about the center of mass of the rigid body, there is a useful labor-saving method to compute the tensor for rotations offset from the center of mass.

If the axis of rotation is displaced by a vector R from the center of mass, the new moment of inertia tensor equals

$$\mathbf{I}^{\mathrm{displaced}} = \mathbf{I}^{\mathrm{center}} + m \left[\left(\mathbf{R} \cdot \mathbf{R}\right) \mathbf{E}_{3} - \mathbf{R} \otimes \mathbf{R} \right]$$

where m is the total mass of the rigid body, E3 is the 3 × 3 identity matrix, and \otimes is the outer product.
Rotational symmetry

Using the above equation to express all moments of inertia in terms of integrals of variables either along or perpendicular to the axis of symmetry usually simplifies the calculation of these moments considerably.
Comparison with covariance matrix
Main article: Moment (mathematics)

The moment of inertia tensor about the center of mass of a 3-dimensional rigid body is related to the covariance matrix of a trivariate random vector whose probability density function is proportional to the pointwise density of the rigid body by:

$$I=n (\mathbf{1}_{3\times 3} \operatorname{tr}(\Sigma) - \Sigma)$$

where n is the number of points.

The structure of the moment-of-inertia tensor comes from the fact that it is to be used as a bilinear form on rotation vectors in the form

$$\frac{1}{2}\omega^\top I \omega\,\!$$

Each element of mass has a kinetic energy of

$$\frac{1}{2} m |v|^2 \,\!$$

The velocity of each element of mass is ω × r where r is a vector from the center of rotation to that element of mass. The cross product can be converted to matrix multiplication so that

$$\omega\times r = [r]_\times^\top \omega\,\!$$

and similarly

$$(\omega\times r)^\top = ([r]_\times \omega)^\top = \omega^\top [r]_\times^\top \,\!$$

Thus,

$$|v|^2 = (\omega\times r)^\top(\omega\times r)=\omega^\top [r]_\times^\top [r]_\times \omega \,\!$$

plugging in the definition of $$[\cdot]_\times the [r]_\times^\top [r]_\times$$ term leads directly to the structure of the moment tensor.

List of moments of inertia
List of moment of inertia tensors
Rotational energy
Parallel axis theorem
Perpendicular axis theorem
Stretch rule
Tire balance
Poinsot's ellipsoid
Instant centre of rotation

Notes

^ Euler, Leonhard (1765-01-01) (in Latin), Theoria motus corporum solidorum seu rigidorum: Ex primis nostrae cognitionis principiis stabilita et ad omnes motus, qui in huiusmodi corpora cadere possunt, accommodata, Cornell University Library, ISBN 978-1-4297-4281-8
^ 3000 Solved Problems in Physics, Schaum Series, A. Halpern, Mc Graw Hill, 1988, ISBN 978-0-07-025734-4

References

Goldstein, H. (1980). Classical Mechanics (2nd ed.). Addison-Wesley. ISBN 0-201-02918-9..
Landau, LD; Lifshitz (1976). Mechanics (3rd ed.). Pergamon Press. ISBN 0-08-021022-8 (hardcover).; ISBN 0-08-029141-4 (softcover).
Marion, JB; Thornton, ST. (1995). Classical Dynamics of Systems and Particles (4th ed.). Thomson. ISBN 0-03-097302-3..
Sylvester, J J (1852), "A demonstration of the theorem that every homogeneous quadratic polynomial is reducible by real orthogonal substitutions to the form of a sum of positive and negative squares", Philosophical Magazine IV: 138–142, retrieved 2008-06-27
Symon, KR (1971). Mechanics (3rd ed.). Addison-Wesley. ISBN 0-201-07392-7..
Tenenbaum, RA (2004). Fundamentals of Applied Dynamics. Springer. ISBN 0-387-00887-X..