The stability and responses of dynamic systems under parametric excitation are often encountered in many fields of science and engineering, such as slender columns and thin plates under axial loadings. This paper proposes a novel numerical simulation method to simultaneously construct stability diagrams and predict the responses of multiple degree-of-freedom (DOF) dynamic systems under arbitrary parametric loadings. The method divides an arbitrary load into discrete segments to approximate the variable excitation function by a step function and then accumulates the system responses of each segment using a matrix method. Dynamic stability and response analysis of undamped and damped 1DOF, 2DOF, 3DOF, and multiple DOF systems are conducted, and numerical examples are compared to conventional methods to show the accuracy and advantage of the proposed numerical simulation method. The instability diagrams are also substantiated by vibration response curves obtained from the same method. The method is applied to calibrate the validity and ranges of applicability of the Hill infinite determinants in the Floquet theory of 1DOF Mathieu-Hill equations. The second order approximation of the first instability region from the Hill determinant is acceptable. However, the fourth order approximation of the second instability region must be used to obtain a relatively accurate result. For 2DOF, 3DOF, and multi-parameter linear periodic systems, the method can simultaneously render four types of instability regions in a single procedure: no resonances, simple parametric resonances, combination additive resonances, and combination differential resonances. These results could provide a better understanding of stability and responses of dynamic systems under parametric excitation.
Read full abstract