A Novel Computational Method to Calculate Nonlinear Normal Modes of Complex Structures