where M and C stand for the mass and damping matrices, u x and u y are the displacements, f x and f y are the excitations and q x and q y are the restoring forces acting in two orthogonal (perpendicular) directions, which are given by