Approximate Solution to Nonlinear Dynamics of a Piezoelectric Energy Harvesting Device Subject to Mechanical Impact and Winkler–Pasternak Foundation

• Published in: Materials Vol. 18; no. 7; p. 1502
• Main Authors: Marinca, Vasile, Herisanu, Nicolae, Marinca, Bogdan
• Format: Journal Article
• Language: English
• Published: Switzerland MDPI AG 27.03.2025
• Subjects: Beam theory (structures); Bifurcations; Convergence; Damping; Dynamic models; Dynamical systems; Electric power production; Energy; Energy harvesting; Euler-Bernoulli beams; Finite element analysis; Force and energy; Foundations; Investigations; Liapunov functions; Magnetic fields; Nonlinear dynamics; Nonlinearity; Parameters; Piezoelectricity; Pontryagin principle; Routh-Hurwitz criterion; Stability; Vibration; Romania; energy harvesting; piezoelectricity; OAFM; Routh–Hurwitz; Lyapunov
• DOI: 10.3390/ma18071502
• ISSN: 1996-1944

To explore the nonlinear dynamics of a piezoelectric energy harvesting device, we consider the simultaneous parametric and external excitations. Based on Bernoulli–Euler beam theory, a new dynamic model is proposed taking into account the curvature of the beam, geometric and electro-mechanical coupling nonlinearities, and damping nonlinearity, with inextensible deformation. The system is discretized by using the Galerkin–Bubnov procedure and then is investigated by the optimal auxiliary functions method. Explicit analytical expressions of the approximate solutions are presented for a complex problem near the primary resonance. The main novelty of our approach relies on the presence of different auxiliary functions, the involvement of a few convergence-control parameters, the construction of the initial and first iteration, and much freedom in selecting the procedure for obtaining the optimal values of the convergence-control parameters. Our procedure proves to be very efficient, simple, easy to implement, and very accurate to solve a complicated nonlinear dynamical system. To study the stability of equilibrium points, the Routh–Hurwitz criterion is adopted. The Hopf and saddle node bifurcations are studied. Global stability is analyzed by the Lyapunov function, La Salle’s invariance principle, and Pontryagin’s principle with respect to the control variables.