Most of the optimisation studies of Vibration Energy Harvesters (VEHs) account for a single output quantity, e.g. frequency bandwidth or maximum power output, but this approach does not necessarily maximise the system efficiency. In those applications where VEHs are suitable sources of energy, to achieve optimal design it is important to consider all these performance indexes simultaneously. This paper proposes a robust and straightforward multi-objective optimisation framework for Vibration Piezoelectric Energy Harvesters (VPEHs), considering simultaneously the most crucial performance indexes, i.e., the maximum power output, efficiency, and frequency bandwidth. For the first time, a rigorous formulation of efficiency for Multi-Degree of Freedom (MDOF) VPEHs is here proposed, representing an extension of previous definitions. This formulation lends itself to the optimisation of FE and MDOF harvesters models. The optimisation procedure is demonstrated using a planar-shape harvester and validated against numerical results. The effects of changing some structural parameters on the harvester performance are investigated via sensitivity analysis. The results show that the proposed methodology can effectively optimise the global performance of the harvester, although this does not correspond to an improvement of every single index. Furthermore, the optimisation of each performance index individually results in a variety of design configurations that greatly differs from one another. It is here demonstrated that the design obtained with the multi-objective function here proposed is similar to the design obtained when optimising the efficiency.

In the industry field, the increasingly stringent requirements of lightweight structures are exposing the ultimately nonlinear nature of mechanical systems. This is extremely true for systems with moving parts and loose fixtures which show piecewise stiffness behaviours. Nevertheless, the numerical solution of systems with ideal piecewise mathematical characteristics is associated with time-consuming procedures and a high computational burden. Smoothing functions can conveniently simplify the mathematical form of such systems, but little research has been carried out to evaluate their effect on the mechanical response of multi-degree-of-freedom systems. To investigate this problem, a slightly damped mechanical two-degree-of-freedom system with soft piecewise constraints is studied via numerical continuation and numerical integration procedures. Sigmoid functions are adopted to approximate the constraints, and the effect of such approximation is explored by comparing the results of the approximate system with the ones of the ideal piecewise counter-part. The numerical results show that the sigmoid functions can correctly catch the very complex dynamics of the proposed system when both the above-mentioned techniques are adopted. Moreover, a reduction in the computational burden, as well as an increase in numerical robustness, is observed in the approximate case.