Lyapunov exponent calculation of a two-degree-of-freedom vibro-impact system with symmetrical rigid stops

A two-degree-of-freedom vibro-impact system having symmetrical rigid stops and subjected to periodic excitation is investigated in this paper. By introducing local maps between different stages of motion in the whole impact process, the Poincar'e map of the system is constructed. Using the Poincar'e map and the Gram-Schmidt orthonormalization, a method of calculating the spectrum of Lyapunov exponents of the above vibro-impact system is presented. Then the phase portraits of periodic and chaotic attractors for the system and the corresponding convergence diagrams of the spectrum of Lyapunov exponents are given out through the numerical simulations. To further identify the validity of the aforementioned computation method, the bifurcation diagram of the system with respect to the bifurcation parameter and the corresponding largest Lyapunov exponents are shown.