Example - Iterative Krylov method

Problem

The given cantilevered beam ($//<![CDATA[ \begin{array}{l}l=0.8\end{array} //]]>$$//<![CDATA[ \begin{array}{l}\mathrm{m}\end{array} //]]>$) has a quadratic cross section ($//<![CDATA[ \begin{array}{l}a=0.01\end{array} //]]>$$//<![CDATA[ \begin{array}{l}\mathrm{m}\end{array} //]]>$) with $//<![CDATA[ \begin{array}{l}E=2.1\times 10^{11}\end{array} //]]>$$//<![CDATA[ \begin{array}{l}\mathrm{N/m^2}\end{array} //]]>$, $//<![CDATA[ \begin{array}{l}\nu=0.3\end{array} //]]>$, $//<![CDATA[ \begin{array}{l}\rho=7850\end{array} //]]>$$//<![CDATA[ \begin{array}{l}\mathrm{kg/m^3}\end{array} //]]>$. It is discretized using the Finite Element Method and excited at the tip with $//<![CDATA[ \begin{array}{l}F=1\end{array} //]]>$$//<![CDATA[ \begin{array}{l}\mathrm{N}\end{array} //]]>$ and of varying circular frequency ($//<![CDATA[ \begin{array}{l}\Omega=\frac{s}{\mathrm{i}}\end{array} //]]>$).

We now include Rayleigh damping ($//<![CDATA[ \begin{array}{l}\mathbf{C}=\alpha\mathbf{M} + \beta\mathbf{K}\end{array} //]]>$) with coefficients $//<![CDATA[ \begin{array}{l}\alpha=1.64\end{array} //]]>$ $//<![CDATA[ \begin{array}{l}\mathrm{1/s}\end{array} //]]>$, $//<![CDATA[ \begin{array}{l}\beta=4.1 e-6\end{array} //]]>$ $//<![CDATA[ \begin{array}{l}\mathrm{s}\end{array} //]]>$.

Solution

We choose $//<![CDATA[ \begin{array}{l}s_{1,2}=1860i \frac{\text{rad}}{\text{s}},\ s_{3,4}=8000\frac{\text{rad}}{\text{s}}\end{array} //]]>$. The expansion points lie between two poles of the transfer function.

It can be seen how IRKA iterates towards the desired result.

Iteration 1

Iteration 2

Iteration 3

Iteration 4

Iteration 5

Iteration 6

Iteration 7

Iteration 8

Iteration 9

Iteration 10