data_05_d_000.csv - Curve of the RMS values of voltage induced on piezoelectric electrodes for p = 0.183 and randomly selected initial conditions δ = 0.0, κ = 0.5. data_05_d_015.csv - Curve of the RMS values of voltage induced on piezoelectric electrodes for p = 0.183 and randomly selected initial conditions δ = 0.15, κ = 0.5. data_05_d_030.csv - Curve of the RMS values of voltage induced on piezoelectric electrodes for p = 0.183 and randomly selected initial conditions δ = 0.3, κ = 0.5. data_05_d_060.csv - Curve of the RMS values of voltage induced on piezoelectric electrodes for p = 0.183 and randomly selected initial conditions δ = 0.6, κ = 0.5. data_06r_o_19.csv - The orbits of the periodic solutions presented in Fig. 6(a) ω = 1.9, Poincaré points = 2. data_06b_o_19.csv - The orbits of the periodic solutions presented in Fig. 6(a) ω = 1.9, Poincaré points = 3. data_06r_o_21.csv - The orbits of the periodic solutions presented in Fig. 6(b) ω = 2.1, Poincaré points = 2. data_06g_o_21.csv - The orbits of the periodic solutions presented in Fig. 6(b) ω = 2.1, Poincaré points = 3. data_06b_o_21.csv - The orbits of the periodic solutions presented in Fig. 6(b) ω = 2.1, Poincaré points = 9. data_06r_o_26.csv - The orbits of the periodic solutions presented in Fig. 6(c) ω = 2.6, Poincaré points = 2. data_06b_o_26.csv - The orbits of the periodic solutions presented in Fig. 6(c) ω = 2.6, Poincaré points = 3. data_06g_o_26.csv - The orbits of the periodic solutions presented in Fig. 6(c) ω = 2.6, Poincaré points = 6. data_06r_o_38.csv - The orbits of the periodic solutions presented in Fig. 6(d) ω = 3.8, Poincaré points = 3. data_06b_o_38.csv - The orbits of the periodic solutions presented in Fig. 6(d) ω = 3.8, Poincaré points = 4. data_06g_o_38.csv - The orbits of the periodic solutions presented in Fig. 6(d) ω = 3.8, Poincaré points = 5. data_06lb_o_38.csv - The orbits of the periodic solutions presented in Fig. 6(d) ω = 3.8, Poincaré points = 7. data_06p_o_38.csv - The orbits of the periodic solutions presented in Fig. 6(d) ω = 3.8, Poincaré points = 9. data_07b_d_015.csv - Numerical results showing the influence of potential asymmetry on the probability of occurrence of particular solutions for δ = 0.15. data_07g_d_030.csv - Numerical results showing the influence of potential asymmetry on the probability of occurrence of particular solutions for δ = 0.30. data_07lb_d_06.csv - Numerical results showing the influence of potential asymmetry on the probability of occurrence of particular solutions for δ = 0.60. data_07r_d_000.csv - Numerical results showing the influence of potential asymmetry on the probability of occurrence of particular solutions for δ = 0.0.

This repository contains the results of numerical simulations of a nonlinear bistable system for harvesting energy from ambient vibrating mechanical sources. Detailed model tests were carried out on an inertial energy harvesting system consisting of a piezoelectric beam with additional springs attached. The mathematical model was derived using the bond graph approach. Depending on the spring selection, the shape of the bistable potential wells was modified including the removal of wells’ degeneration. Consequently, the broken mirror symmetry between the potential wells led to additional solutions with corresponding voltage responses. The probability of occurrence for different high voltage/large orbit solutions with changes in potential symmetry was investigated. In particular, the periodicity of different solutions with respect to the harmonic excitation period were studied and compared in terms of the voltage output. The results showed that a large orbit period-6 subharmonic solution could be stabilized while some higher subharmonic solutions disappeared with the increasing asymmetry of potential wells. Changes in frequency ranges were also observed for chaotic solutions.