I would like to find the frequency response of a spring mass system of multiple degrees of freedom by using the Laplace transform. I think I know how to do this with one mass oscillating, however I find difficulty in extending the approach to multiple masses.
My approach for one mass
The equation of motion for one damped spring mass is as follows: $$m\ddot{x}(t) + c\dot{x}(t) + kx(t) = f(t) \tag{1}$$ To solve this in the frequency domain, I employ the Laplace transform on Equation (1) above: $$ms^2X(s) + csX(s) + kX(s) = F(s) \\ \Rightarrow H(s) = \frac{X(s)}{F(s)}=\frac{1}{s^2m + sc + k} \tag{2}$$ Substituting $j\omega$ in the place of $s$, and taking the modulus of the transfer function $H(s)$ above, we obtain the frequency response for the single mass for all $\omega$: $$|H(j\omega)|=\left| \frac{1}{-\omega^2m + j\omega c+ k} \right| \tag{3}$$
My approach for multiple masses
I pick up where I left of from Equation (2), however this time I replace the masss, damping and stiffness coefficients with their respective matrices, and make my substitution for $s$ as in Equation (3) above: $$[-\omega^2\mathbf{M} + j\omega\mathbf{C} + \mathbf{K}]\mathbf{X}=\mathbf{F}$$
However this is where I get stuck. Intuitively I can tell that my transfer function will be something like $\mathbf{H} =[-\omega^2\mathbf{M} + j\omega\mathbf{C} + \mathbf{K}]^{-1}$, however firstly I don't know how to prove it, and secondly I don't know how from this matrix, I will end up getting a vector for every value of $\omega$ that I insert in this equation, such that I get the frequency response for each degree of freedom
Best Answer
In frequency space a one dimensional linear system is characterized by the dynamical equation
$$ A(s)X(s) = F(s) $$
Which can typically be found from the equations of motion and a Laplace transform. Inverting this expression we find
$$ X(s) = A(s)^{-1}F(s) = H(s)F(s) $$
Where I've defined $H(s) = A(s)^{-1}$. $H(s)$ is the system transfer function or susceptibility. It tells us the response of the linear system (characterized by variable $X(s)$) due to excitation $F(s)$ at a particular frequency.
The system becomes more physically complicated when multiple degrees of freedom are included but the math is not much more complicated. The dynamical equation (from equations of motion and Laplace transform) is now
$$ \boldsymbol{A}(s) \boldsymbol{X}(s) = \boldsymbol{F}(s) $$
Here $\boldsymbol{A}$ is an $n\times n$ matrix To find the susceptibility matrix we invert
$$ \boldsymbol{X}(s) = \boldsymbol{A}(s)^{-1}\boldsymbol{F}(s) = \boldsymbol{H}(s)\boldsymbol{F}(s) $$
The susceptibility matrix $\boldsymbol{H}$ is an $n\times n$ matrix. The $(i,j)$ component of this matrix, $H_{ij}$ tells us the response of variable $X_i$ to "force" $F_j$.
If we want to "visualize" the susceptibility matrix we, in general, must plot all $n\times n$ components $H_{ij}(s)$ versus frequency $s$. This is the analog of the plot of $H(s)$ in the 1d case.