Approximate Methods for Multiple Degree of Freedom Systems: Dunkerley’s Formula
Dunkerley’s Formula is another method of estimating the lowest (fundamental) natural frequency of a system without having to solve an eigenvalue problem. Rather than using the stiffness matrix, Dunkerley’s Method makes use of the flexibility matrix which is the inverse of the stiffness matrix.
Starting with the equations of motion for a general degree of freedom system
and assuming simple simultaneous harmonic motion of the form
Premultiplication of (9.9) by results in
As we have seen previously, for a non-trivial solution to exist it is required that
For simplicity we will assume that the mass matrix is diagonal. In that case, the determinant in (9.10) has the form
Expanding (9.11) produces the characteristic equation which has the form
where and are constants (which will depend on and ) which we do not need to consider here.
Let the roots of (9.12) be given by
Since these values are the roots, the characteristic equation could also be written in the form
To this point no approximations have been made. Now however, we note that for many systems the higher natural frequencies are significantly larger than the fundamental natural frequency so that
As a result, ignoring all but the first term on the LHS of (9.15) results in
or more generally
Equation (9.16) is known as Dunkerley’s Formula and provides an estimate of the fundamental natural frequency of a system. Due to the terms that have been neglected, the natural frequency obtained using Equation (9.16) will be lower than the actual fundamental natural frequency. Dunkerley’s Formula therefore provides a lower bound on the lowest natural frequency of the system (in contrast to Rayleigh’s Quotient which provided an upper bound on the lowest natural frequency).
Note that as long as the mass matrix is diagonal, Equation (9.16)
can be written as
where is defined as the natural frequency of the single degree of freedom system that results when mass is the only mass in the system (i.e. all other masses in the system are assumed to be zero). This is because of the fact that if is the only mass in the system then represents the deflection of the mass (in the positive i coordinate direction) when a unit load is applied to the mass (again in the positive i coordinate direction). The effective stiffness of the (single degree of freedom) system is then
as defined above. This interpretation can be useful in applying Dunkerley’s Formula in some instances.
Estimate the fundamental natural frequency for the system shown below using Dunkerley’s Method.
Consider the following cases:
Note that the stiffness matrix for this problem is
so that is (Note: show full screen view)
Use Dunkerley’s method to estimate the fundamental natural frequency for a light simply supported beam supporting three equally spaced identical masses.
Note that for this system
Dunkerley’s Formula to Estimate The Highest Natural Frequency
The ideas behind Dunkerley’s Formula, with a slightly different formulation, can also be used to estimate the highest natural frequency in a system. To see this start with (9.9)
Premultiplication by produces
This again requires, for non-trivial solutions,
If we again assume that the mass matrix is diagonal, then the inverse is simply
As a result, the determinant in (9.17) has the form
Expanding (9.18) produces a characteristic equation with the form
where and are again constants which depend on and . The roots of (9.19) correspond to the natural frequencies of the system given by
Once again, since these values are the roots of the characteristic equation, we could equivalently write the characteristic equation as
which expanded gives
Since is the largest natural frequency in the system, we ignore all of the other terms on the LHS of (9.22) to get
or more generally
Due to the terms that have been ignored in (9.22), the frequency obtained using (9.23) will be larger than the highest natural frequency in the system. As a result, this version of Dunkerley’s Formula provides an upper bound on the largest natural frequency of the system.
For the system shown below:
(a) Use Rayleigh’s quotient to find an upper bound to the lowest natural frequency. Use a trial vector of
(b) Use Dunkerley’s method to find a lower bound for the lowest natural frequency.
(c) Use Dunkerley’s method to find an upper bound for the highest natural frequency.
Note that for this problem,