Mode computation of immersed multilayer plates by solving an eigenvalue problem