Misplaced Pages

This article has multiple issues. Please help improve it or discuss these issues on the talk page. (Learn how and when to remove these messages) This article may require cleanup to meet Misplaced Pages's quality standards. The specific problem is: Needs linked into Misplaced Pages. Please help improve this article if you can. (February 2017) (Learn how and when to remove this message) This article includes a list of general references, but it lacks sufficient corresponding inline citations. Please help to improve this article by introducing more precise citations. (February 2017) (Learn how and when to remove this message) This article needs additional citations for verification. Please help improve this article by adding citations to reliable sources. Unsourced material may be challenged and removed. Find sources: "Rayleigh's quotient in vibrations analysis" – news · newspapers · books · scholar · JSTOR (February 2017) (Learn how and when to remove this message) (Learn how and when to remove this message)

The Rayleigh's quotient represents a quick method to estimate the natural frequency of a multi-degree-of-freedom vibration system, in which the mass and the stiffness matrices are known.

The eigenvalue problem for a general system of the form $${\displaystyle M\,{\ddot {\textbf {q}}}(t)+C\,{\dot {\textbf {q}}}(t)+K\,{\textbf {q}}(t)={\textbf {Q}}(t)}$$ in absence of damping and external forces reduces to $${\displaystyle M\,{\ddot {\textbf {q}}}(t)+K\,{\textbf {q}}(t)=0}$$

The previous equation can be written also as the following: $${\displaystyle K\,{\textbf {u}}=\lambda \,M\,{\textbf {u}}}$$ where ${\displaystyle \lambda =\omega ^{2}}$, in which ${\displaystyle \omega }$ represents the natural frequency, M and K are the real positive symmetric mass and stiffness matrices respectively.

For an n-degree-of-freedom system the equation has n solutions ${\displaystyle \lambda _{m}}$, ${\displaystyle {\textbf {u}}_{m}}$ that satisfy the equation $${\displaystyle K\,{\textbf {u}}_{m}=\lambda _{m}\,M\,{\textbf {u}}_{m}}$$

By multiplying both sides of the equation by ${\displaystyle {\textbf {u}}_{m}^{T}}$ and dividing by the scalar ${\displaystyle {\textbf {u}}_{m}^{T}\,M\,{\textbf {u}}_{m}}$, it is possible to express the eigenvalue problem as follow: $${\displaystyle \lambda _{m}=\omega _{m}^{2}={\frac {{\textbf {u}}_{m}^{T}\,K\,{\textbf {u}}_{m}}{{\textbf {u}}_{m}^{T}\,M\,{\textbf {u}}_{m}}}}$$ for m = 1, 2, 3, ..., n.

In the previous equation it is also possible to observe that the numerator is proportional to the potential energy while the denominator depicts a measure of the kinetic energy. Moreover, the equation allow us to calculate the natural frequency only if the eigenvector (as well as any other displacement vector) ${\displaystyle {\textbf {u}}_{m}}$ is known. For academic interests, if the modal vectors are not known, we can repeat the foregoing process but with ${\displaystyle \lambda =\omega ^{2}}$ and ${\displaystyle {\textbf {u}}}$ taking the place of ${\displaystyle \lambda _{m}=\omega _{m}^{2}}$ and ${\displaystyle {\textbf {u}}_{m}}$, respectively. By doing so we obtain the scalar ${\displaystyle R({\textbf {u}})}$, also known as Rayleigh's quotient: $${\displaystyle R({\textbf {u}})=\lambda =\omega ^{2}={\frac {{\textbf {u}}^{T}\,K\,{\textbf {u}}}{{\textbf {u}}^{T}\,M\,{\textbf {u}}}}}$$

Therefore, the Rayleigh's quotient is a scalar whose value depends on the vector ${\displaystyle {\textbf {u}}}$ and it can be calculated with good approximation for any arbitrary vector ${\displaystyle {\textbf {u}}}$ as long as it lays reasonably far from the modal vectors ${\displaystyle {\textbf {u}}_{i}}$, i = 1,2,3,...,n.

Since, it is possible to state that the vector ${\displaystyle {\textbf {u}}}$ differs from the modal vector ${\displaystyle {\textbf {u}}_{m}}$ by a small quantity of first order, the correct result of the Rayleigh's quotient will differ not sensitively from the estimated one and that's what makes this method very useful. A good way to estimate the lowest modal vector ${\displaystyle (u_{1})}$, that generally works well for most structures (even though is not guaranteed), is to assume ${\displaystyle (u_{1})}$ equal to the static displacement from an applied force that has the same relative distribution of the diagonal mass matrix terms. The latter can be elucidated by the following 3-DOF example.

Example – 3DOF

As an example, we can consider a 3-degree-of-freedom system in which the mass and the stiffness matrices of them are known as follows: $${\displaystyle M={\begin{bmatrix}1&0&0\\0&1&0\\0&0&3\end{bmatrix}}\;,\quad K={\begin{bmatrix}3&-1&0\\-1&3&-2\\0&-2&2\end{bmatrix}}}$$

To get an estimation of the lowest natural frequency we choose a trial vector of static displacement obtained by loading the system with a force proportional to the masses: $${\displaystyle {\textbf {F}}=k{\begin{bmatrix}m_{1}\\m_{2}\\m_{3}\end{bmatrix}}=1{\begin{bmatrix}1\\1\\3\end{bmatrix}}}$$

Thus, the trial vector will become $${\displaystyle {\textbf {u}}=K^{-1}{\textbf {F}}={\begin{bmatrix}2.5\\6.5\\8\end{bmatrix}}}$$ that allow us to calculate the Rayleigh's quotient: $${\displaystyle R={\frac {{\textbf {u}}^{T}\,K\,{\textbf {u}}}{{\textbf {u}}^{T}\,M\,{\textbf {u}}}}=\cdots =0.137214}$$

Thus, the lowest natural frequency, calculated by means of Rayleigh's quotient is: $${\displaystyle w_{\text{Ray}}=0.370424}$$

Using a calculation tool is pretty fast to verify how much it differs from the "real" one. In this case, using MATLAB, it has been calculated that the lowest natural frequency is: ${\displaystyle w_{\text{real}}=0.369308}$ that has led to an error of ${\displaystyle 0.302315\%}$ using the Rayleigh's approximation, that is a remarkable result.

The example shows how the Rayleigh's quotient is capable of getting an accurate estimation of the lowest natural frequency. The practice of using the static displacement vector as a trial vector is valid as the static displacement vector tends to resemble the lowest vibration mode.

References

1. Meirovitch, Leonard (2003). Fundamentals of Vibration. McGraw-Hill Education. p. 806. ISBN 9780071219839.

• Abstract algebra
• Linear algebra
• Mathematical physics
• Mechanical vibrations