Dynamic Modeling and Structural Optimization of a Bistable Electromagnetic Vibration Energy Harvester

: A novel bistable electromagnetic vibration energy harvester (BEMH) is constructed and optimized in this study, based on a nonlinear system consisting mainly of a ﬂexible membrane and a magnetic spring. A large-amplitude transverse vibration equation of the system is established with the general nonlinear geometry and magnetic force. Firstly, the mathematical model, considering the higher-order nonlinearities given by nonlinear Galerkin method, is applied to a membrane with a co-axial magnet mass and magnetic spring. Secondly, the steady vibration response of the membrane subjected to a harmonic base motion is obtained, and then the output power considering electromagnetic e ﬀ ect is analytically derived. On this basis, a parametric study in a broad frequency domain has been achieved for the BEMH with di ﬀ erent radius ratios and membrane thicknesses. It is demonstrated that model predictions are both in close agreement with results from the ﬁnite element simulation and experiment data. Finally, the proposed e ﬃ cient solution method is used to obtain an optimizing strategy for the design of multi-stable energy harvesters with the similar ﬂexible structure.


Introduction
With the rapid development of society, the issue of energy shortage and environmental pollution is becoming more and more concerned, which highlights the importance of energy harvesting technologies. Energy harvesting technology is to collect and convert environmental energy, such as vibration mechanical energy, magnetic energy, friction energy, temperature difference energy, wind energy, ocean energy, solar energy and so on [1,2], which are contained in the external environment or human body, into electricity through power generation devices to meet the energy supply of micro-electronic devices [3]. In the new micro-energy devices, the vibration energy harvester does not need the external power supply, and the energy required for working is fully self-sufficient. It has the advantages of long service life, no frequent replacement and high reliability [4], which provides a solution to the energy supply problem of micro-electronic systems [5][6][7]. The subject has been extensively studied by scholars from different countries [8,9]. According to different working principles, the current vibration energy harvesters are mainly divided into piezoelectric, electromagnetic, capacitive, and magnetostrictive [10][11][12][13][14][15][16]. The electromagnetic vibration energy harvester is composed of a permanent magnet, inductive coil winding and elastic elements. Under external excitation, the relative motion between a permanent magnet and the coil results in the change of magnetic flux through the coil, which generates inductive electromotive force [17]. Due to its low-frequency, simple structure, low-cost and no power-supply, the electromagnetic vibration energy harvester has wide application prospects, and attracts high research attention [18,19].
In the development of vibration energy harvester, different scholars have various research emphases, but their goal is always to expand working bandwidth and improve output power [20][21][22][23]. Traditional vibration energy harvesters are generally monostable systems with the shortcomings of the single resonance mode, narrow bandwidth and insufficient power output. Nonlinear energy harvesters with mono-stable, bi-stable, tri-stable, quad-stable and even penta-stable characteristics [24][25][26][27] have been extensively studied such as to improve the harvesting performance in natural environments. For example, the bistable energy harvesting system can achieve large-amplitude, broadband periodic and even chaotic vibrations under non-resonant excitation compared with monostable system so that the output power of the bistable energy harvesting system can be significantly improved [28]. Dhakar L et al. [29] studied the piezoelectric energy harvesting device with composite cantilever beam, which shows that the nonlinear effect can improve the energy harvesting efficiency of the device at a certain frequency. The influence of nonlinear effect on the piezoelectric and electromagnetic structure under random excitation is studied by De Paula A S et al. [30], who find that the nonlinear effect strengthen the performance of bistable energy harvesting system. Zhou S et al. [31] designed a tristable energy harvesting device with shallow potential well, and the result shows that the tristable energy harvester has a wide working bandwidth and collects energy well. The nonlinear electromagnetic force is a substitute for the traditional spring force to research a nonlinear electromagnetic vibration energy harvester [32]. A nonlinear energy harvester is designed with permanent magnets for the sake of covering a wider frequency range, and the simulation results are verified by experiments [33]. The bistable energy harvester with an additional linear oscillator is studied by Harne R L et al. [34], who also find that the device can increase the response amplitude to improve the harvesting power, and the differences of mass ratio and frequency modulation ratio have a great influence on the response characteristics of the system. Although the theoretical and experimental objects of the mentioned bistable vibration energy harvesting system are different, the simplified mechanical models are basically single-degree-of-freedom systems based on spring-mass model. These models can simplify the vibration equation of the system, and provide general technique opinions [35][36][37]. However, they may cause error in case of large-amplitude vibration of continuum structure, in which the influence of higher-order modes on structural vibration should not be neglected [38]. In related researches on the vibration of membranes, many studies are the small-amplitude vibration of membranes [39,40]. For the large deflection of membranes, the general researches consider the deformation of membranes under different loads, which belongs to the static problem [41][42][43]. However, there has been little research considering the large-amplitude vibration response of membranes with co-axial mass. In that case, higher-order components become necessary for computation. If the higher-order nonlinearities are not taken into account under large-amplitude vibration, the complex dynamic behaviors of the system may be evaluated imprecisely and the corresponding subsequent calculation results may also lead to larger deviations [44].
In order to study the nonlinear dynamic characteristics and reliable power optimization criteria for the energy harvester, this paper combines a flexible membrane with a magnetic spring to establish a distributed parameter model for the bistable electromagnetic vibration energy harvester (BEMH). A reduced-order model of the peripheral fixed-supported membrane with a co-axial magnet is obtained using the nonlinear Galerkin method to take into account the effect of higher-order vibration mode. It is verified that the vibration characteristics and power output of this proposed model are more consistent with that of the finite element simulation results than traditional Galerkin model. Along with the high computational efficiency of the model allows the variation characteristics of output power are investigated under different variants of main structural parameters: the radius ratio and membrane thickness, resulting in a general rule for power output optimization.

110
The working principle of the BEMH is that the co-axial permanent magnet driven by the 111 membrane vibrates up and down in the cavity wall, which changes the magnetic flux of the 112 conductive coil windings around the upper and lower end covers to generate inductive 113 electromotive force. If the polyimide membrane is replaced by piezoelectric membrane, the energy 114 harvester will become a hybrid vibration energy harvester.

117
The BEMH studied in this paper can be simplified to the mechanical model shown in Error!

118
Reference source not found., which is composed of a membrane surrounding with a co-axial 119 magnet. The membrane radius and thickness are b and h , respectively. The radius and mass of 120 the magnet are a and c m , respectively.

121
If the magnetic force is large enough, the effect of co-axis mass gravity on the system can be 125 where 0 Z and Ω are respectively the amplitude and angular frequency of excitation acceleration.

126
To represent the large-amplitude axisymmetric vibrations of the plate, the vibration governing The working principle of the BEMH is that the co-axial permanent magnet driven by the membrane vibrates up and down in the cavity wall, which changes the magnetic flux of the conductive coil windings around the upper and lower end covers to generate inductive electromotive force. If the polyimide membrane is replaced by piezoelectric membrane, the energy harvester will become a hybrid vibration energy harvester.
The BEMH studied in this paper can be simplified to the mechanical model shown in Figure 2, which is composed of a membrane surrounding with a co-axial magnet. The membrane radius and thickness are b and h, respectively. The radius and mass of the magnet are a and m c , respectively. The configuration of BEMH is outlined in Figure 1, which is mainly composed of a membrane, a 100 rigid magnet, two iron coil windings and a cavity wall. The polyimide membrane is bonded to the 101 clamp by epoxy resin, and the NdFeB magnet is adsorbed on both sides of the polyimide membrane.

110
The working principle of the BEMH is that the co-axial permanent magnet driven by the 111 membrane vibrates up and down in the cavity wall, which changes the magnetic flux of the 112 conductive coil windings around the upper and lower end covers to generate inductive 113 electromotive force. If the polyimide membrane is replaced by piezoelectric membrane, the energy 114 harvester will become a hybrid vibration energy harvester.

115
. 116  If the magnetic force is large enough, the effect of co-axis mass gravity on the system can be neglected. The energy harvesting efficiency of the BEMH can be investigated by setting a cosine excitation as a simple representative of the harmonic base motion. Mechanical vibration input for the system is achieved through: ..
where Z 0 and Ω are respectively the amplitude and angular frequency of excitation acceleration.
To represent the large-amplitude axisymmetric vibrations of the plate, the vibration governing equations are written in the form of Berger equations [45]: where D = Eh 3 12(1−ν 2 ) is bending stiffness of the circular plate. From the stress-strain relations, the time-dependent function may be expressed as α 2 = N r +N t (1+ν)D , where N r and N t represent radial and tangential stress, respectively.
Different from the generally employed Berger equations, in this article, it takes the transverse inertia term, foundation excitation, the mass inertia of the co-axial magnet, the mechanical damping, and the action of the nonlinear magnetic spring into consideration. Considering that the thickness of the membrane is small, the bending stiffness is neglected, i.e., D = 0. To describe the forced lateral axisymmetric vibration of a circular membrane with a center magnet, the governing equations can be written in the form like: ..
is the Kronecker delta, w(r, t) is the transverse displacement, u(r, t) is the radial displacement, E is the elasticity modulus, ν is Poisson's ratio, T is the tension of the membrane, ρ is the density of the membrane, and c m is the damping coefficient.
When the rigid magnet sustains an up-down symmetric transverse vibration, a vertical force is given in Figure 3. The Equations (4) and (5) are subjected to the following boundary conditions: 2πaT ∂w(r,t) ∂r r=a = m c ∂ 2 w(r,t) ∂t 2 , w(r, t) r=b = 0, u(r, t) r=a = u(r, t) r=b = 0, w(r, t) t=0 = u(r, t) t=0 = 0 (6) Based on the cubic nonlinear magnetic force, the magnetic force between magnets and iron 146 cores is expressed as follows [46]: where 1 β and 2 β represent the linear stiffness coefficient and the nonlinear stiffness coefficient 148 provided by the magnetic spring, respectively.

149
Equation Error! Reference source not found. can be rewritten in the equivalent form: Based on the cubic nonlinear magnetic force, the magnetic force between magnets and iron cores is expressed as follows [46]: where β 1 and β 2 represent the linear stiffness coefficient and the nonlinear stiffness coefficient provided by the magnetic spring, respectively. Equation (5) can be rewritten in the equivalent form: Because k is only a time-dependent function, one may thus multiply Equation (8) by rdr and integrate from the radius a to the outside radius b: b a ∂(ur) ∂r When the membrane is clamped, the radial displacement u(r, t) vanishes at the boundaries r = a and r = b. The left side of Equation (9) may be integrated by parts and expressed in the alternate form: The time-dependent function k can be obtained by Equation (10). Based on the standard method of separation of variables, the relative vibration response of the membrane can be separated to the spatial and temporal functions as the following series form: where q i (t) is the time-dependent response, and φ i (r) is the i-th mode function of vibration of the membrane.
To investigate the dynamic response of the BEMH system, the natural frequencies of the membrane and associated mode eigenfunction are obtained by the undamped free vibration Equation of the membrane, as follows [47]: Substituting Equation (11) into Equation (12), the mode function φ i (r) is given by: The associated boundary conditions can be formulated as: where ω i is the undamped natural frequency of the i−th vibration mode, and the eigenvalue of i−th vibration mode is given by: Using the differential eigenvalue problem, the eigenfunction is expressed as [48]: where the coefficient σ i is given by: Regarding the eigenfunction and the boundary conditions, the eigenequation of the differential eigenvalue problem can be obtained as follows: It should be noted that, according to the orthogonality conditions, the eigenfunction can be mass normalized and stiffness normalized respectively as followed [49]: where δ ij is the Kronecker delta, and the coefficient C i can be calculated by using the above orthogonality relations. Then, substituting Equations (1), (7) and (10) into Equation (4), the distributed parameter model of transverse vibration for the BEMH can be given by: Substituting Equation (11) into Equation (21), multiplying both sides by φ i (r)r and integrating over the radius of the membrane, and considering the orthogonality relations of Equations (19) and (20), the following ordinary differential equation is obtained: ..
where the parameters in Equation (22) are: and the mechanical damping ratio ζ j of the j−th mode of vibration can be obtained by experimental measurement and m e is the equivalent mass of the membrane and the magnet. For large deflection vibration of the membrane, the influence of higher-order nonlinear components on dynamic characteristics of the system should not be neglected. However, the traditional Galerkin discretization method usually ignores the influence of higher-order components, and obtains the lower-order model by its reduction is inaccurate in the case of large-amplitude vibration of the membrane, which will lead to greater deviations in power optimization analysis based on reduced-order model. Therefore, the nonlinear Galerkin method is chosen to reduce the order of the model, considering higher-order nonlinearities when the membrane undergoes large-amplitude vibration.

Reduced-Order Model
Many practical structures are infinite-dimensional systems, which usually have the characteristics of dense and continuous modes. In order to discuss the vibration of the structures, order reduction methods are often used to obtain a lower-order model but sufficient precision.
Convergence problem is the primary concern when dealing with continuum model by dimension reduction method. The first three-order vibration modes are intercepted for analysis. As shown in Figure 4, the contribution of the first order component is 84%, the second order component is 14%, and the third order component is only 0.9%. However, the computational complexity is increasing rapidly. To guarantee the accuracy and convergence rate of calculation of the system, the second-order truncation is chosen to calculate and analyze the model.

221
The traditional Galerkin method usually truncates Equation

222
Error! Reference source not found. directly to obtain the reduced-order model without considering 223 the influence of higher-order components. For example, when = = n j 1, 1 , the reduced-order 224 single-degree-of-freedom (SDOF) system is given by: In traditional Galerkin method, the low-order model is often based on experience interception directly, without considering the influence of higher-order components. However, there is a drawback of the direct truncation method, ignoring the higher-order components influences on structure, and getting the imprecise conclusions for a nonlinear system when a large-amplitude vibration occurs. Therefore, Marion and Temam [50] proposed the nonlinear Galerkin method based on the traditional Galerkin method to refine the reduced-model. Referring to the construction of the approximate inertial manifold, the key issue of the nonlinear Galerkin method is to describe the original system with the low-order model and consider the influence of higher-order components in the same order equation. Accordingly, it has higher accuracy and better convergence speed than the traditional method. The major difference between the two methods can be found in the flowchart of Figure 5. Convergence problem is the primary concern when dealing with continuum model by 201 dimension reduction method. The first three-order vibration modes are intercepted for analysis. As

221
The traditional Galerkin method usually truncates Equation

222
Error! Reference source not found. directly to obtain the reduced-order model without considering 223 the influence of higher-order components. For example, when = = n j 1, 1 , the reduced-order 224 single-degree-of-freedom (SDOF) system is given by: 225 Figure 5. The difference between two methods of dimension reduction.
The traditional Galerkin method usually truncates Equation (22) directly to obtain the reduced-order model without considering the influence of higher-order components. For example, when n = 1, j = 1, the reduced-order single-degree-of-freedom (SDOF) system is given by: ..
In contrast, the reduced-order SDOF equation obtained by nonlinear Galerkin method includes not only lower-order components but also higher-order components of the system. In particular, the higher-order components are taken as second-order modes, i.e., n = 2. For j = 1, 2, the two-degree-of-freedom differential equations are as follows: .. ..
The approximate inertial manifold reflecting the relation between q 1 (t) and q 2 (t) is solved by means of the fixed point method. Because q 2 (t) is a small quantity, considering that the derivative of high frequency components for time is very small, let .. q 2 (t) = 0 and . q 2 (t) = 0 , ignoring the second-order and above infinitesimals of q 2 (t), the result is given by: Substituting Equation (30) into Equation (25), the reduced-order model solved by nonlinear Galerkin method is obtained. Because this SODF equation is very large and complex, it is not listed separately here.

The Output Displacement Response and Power of the BEMH
The transverse vibration displacement response of the system, that is, the approximate solution of Equation (21) can be obtained by numerical solution of the reduced-order SODF equation as follows: Involving the electromagnetic damping, the mechanical energy of the vibration system can be converted into electric energy. The instantaneous output power and average output power over one vibration period, transferred from the vibration source to the electromagnetic damping, can be expressed respectively as [51]: w(a, t) 2 (32) where .
w(a, t) is the vibration velocity of the magnet, and the electromagnetic damping coefficient can be expressed as: where R coil is coil resistance, R load is external resistance, and Θ is the commonly assumed as the electromagnetic coupling constant [52] between the moving magnet and coil windings, which reflects the strength of electromagnetic interaction. Electric power is usually considered to be the output power, which is upon the assumption: R load R coil . In the structure of the BEMH, the distance between the co-axial magnet and the upper and lower iron cores are variable to adjust the magnetic force on the co-axial magnet by changing the relative position between the upper and lower end covers. The energy harvester is still a monostable system with small magnetic force, and turns to a bistable model as the magnetic force increases to a certain range. Adjusting the stiffness of magnetic spring to keep the system in bistable state and fixing the upper and lower iron cores, the magnetic force of the co-axial levitated magnet is obtained using Ansoft Maxwell software for finite element analysis of the electromagnetic field in this paper. According to the vibration range of the magnet, the magnetic simulation and polynomial fit results are shown in Figure 6.

238
The transverse vibration displacement response of the system, that is, the approximate solution where & w a t ( , ) is the vibration velocity of the magnet, and the electromagnetic damping coefficient 246 can be expressed as: where coil R is coil resistance, load R is external resistance, and Θ is the commonly assumed as the 248 electromagnetic coupling constant [52] between the moving magnet and coil windings, which 249 reflects the strength of electromagnetic interaction. Electric power is usually considered to be the 250 output power, which is upon the assumption: ?
load coil R R .

251
In the structure of the BEMH, the distance between the co-axial magnet and the upper and  In order to analyze the vibration response characteristics and the generated electric power across the resistive load of the EBMH under different harmonic excitations, the numerical comparison is carried out for the reduced-order system via the traditional Galerkin method and nonlinear Galerkin method. Then the dynamic characteristics obtained by the mathematical models are comparatively analyzed with the finite element simulation results. In this system, the physical parameters in numerical solution and simulation analysis are designated consulting the entity structure as shown in Table 1.

Finite Element Simulation Validations and Numerical Results
The bifurcation diagram of amplitude-frequency response relation is obtained, calculating the dynamic displacement response of the reduced-order model based on the nonlinear Galerkin method at different frequencies, as shown in Figure 7. Through the dynamic analysis of the theoretical model, we know that under the excitation amplitude of 5 m/s 2 , the system cannot obtain enough energy at a small excitation frequency. In this stage, the system makes small-amplitude vibration near one of its equilibrium points as shown in Figure 8a. With the increase of excitation frequency, the system gains enough energy from the excitation and begins to transfer from the bistable movement to large-amplitude steady vibration. Among them, the critical frequency value of the bistable movement is 6 Hz as shown in Figure 8b. Then, as the excitation frequency continues to increase, the system begins to vibrate substantially and steadily as shown in Figure 8c. That results verify the low frequency character of the membrane harvester, and the advantage of bistable movement in improving the output performance. Accordingly, the architecture model in this article is very suitable for the common engineering applications. In order to verify the theoretical model, the finite-element simulation model of the membranemagnet is established by using the finite-element software ANSYS Workbench. The transient dynamic analysis of the model under different excitation conditions is carried out. The results of vibration response are shown in Figure 9. As the excitation frequency reaches 6 Hz, the system undergoes the bistable vibration. Then, the system generates large-amplitude stable vibration with the increase of excitation frequency.  For the comparison between the results in Figures 8 and 9, the excitation frequency calculated by the mathematical model is basically coherent with the finite-element calculation result, at which the BEMH begins to be the bistable vibration under the same excitation amplitude. When Z 0 = 6 m/s 2 and Ω/2π = 13 Hz, the amplitude of large periodic motion of the system is also similar to the finite-element result. Therefore, the obtained computational dynamics model is effective to analyze the dynamic behavior of the BEMH.
The dynamic characteristics of large-amplitude vibration and power output of the membrane and magnetic spring system are numerically studied by nonlinear Galerkin methods and traditional Galerkin method under same excitation amplitude of 5 m/s 2 and frequency of 13 Hz. The BEMH generates large-amplitude and steady vibration under this condition, and the vibration displacement and instantaneous output power of the two models are compared as shown in Figure 10. It is found that the displacement responses and instantaneous output powers of the two models are quite different, and the deviations of instantaneous powers are non-negligible.
finite-element calculation result, at which the BEMH begins to be the bistable vibration under the 313 same excitation amplitude. When 0 Z =6 m/s 2 and Ω π / 2 =13 Hz, the amplitude of large periodic 314 motion of the system is also similar to the finite-element result. Therefore, the obtained 315 computational dynamics model is effective to analyze the dynamic behavior of the BEMH.

316
The dynamic characteristics of large-amplitude vibration and power output of the membrane 317 and magnetic spring system are numerically studied by nonlinear Galerkin methods and traditional   Within the excitation frequency range of the large-amplitude and stable vibration of the system, the vibration displacement and average output power are calculated by different methods. Figure 11 shows that the amplitude results obtained by the model of traditional Galerkin method are within 18.5% error of the finite-element simulation results, and the results analyzed by the model of nonlinear Galerkin method is within 4.6% error. Therefore, the calculation model using nonlinear Galerkin method is closer to the finite-element model.   As is shown in Figure 10b, the instantaneous output powers of the BEMH obtained by the model of traditional Galerkin method and nonlinear Galerkin method also have major differences. The average power generated by the BEMH in one vibration period is calculated utilizing three models as shown in Figure 12, in which the maximum calculation errors gained by the models of traditional Galerkin method and nonlinear Galerkin method are 31.1% and 5.4% respectively, compared with the results of the finite-element analysis. In conclusion, the dynamic model reduced via nonlinear Galerkin method is more accurate and available for the dynamic analysis of the large-amplitude vibrations on membrane structures.

Experimental Results
The main objective of the experiment is to validate the theoretical model and investigate the both quantitative comparison between the theoretical model and qualitative agreement with numerical simulation predictions. The laboratory-scale experimental devices and the prototype of energy harvester as shown in Figure 13 are designed and fabricated, whose geometry property is listed in Table 1

360
In the signal generator, the input frequency is 6 Hz, and the voltage is 3 V. According to the 361 previous numerical calculation, the system can produce bistable motion only in a certain interval.  In the signal generator, the input frequency is 6 Hz, and the voltage is 3 V. According to the previous numerical calculation, the system can produce bistable motion only in a certain interval. Therefore, the gain knob of the signal amplifier makes the amplitude of the sinusoidal motion of the prototype increase from zero until bistable vibration occurs. The amplitude and the frequency of the excitation signal is recorded under bistable vibration. The bistable signal generated by the system is shown in Figure 14. In that case, the power converted from the measured voltage is very close to the theoretical value. previous numerical calculation, the system can produce bistable motion only in a certain interval.

369
Under different excitation frequencies, the excitation amplitudes are different required for the 370 system to achieve bistable vibration. As can be seen from   Under different excitation frequencies, the excitation amplitudes are different required for the system to achieve bistable vibration. As can be seen from Table 2, the model obtained by this nonlinear Galerkin method can be used to solve the bifurcation parameters of the system in which the bistable vibration occurs. Within the excitation frequency range of the large-amplitude and stable vibration of the system. Figure 15 shows measured output power of the system at two different accelerations, i.e., 3 m/s 2 and 4 m/s 2 , respectively. It also compares the higher-order nonlinear model predictions of output power against measured data. In general, at both accelerations, 3 m/s 2 and 4 m/s 2 , experiment and mathematical model show similar trends.

Parametric Study of the BEMH
With the validated model of the nonlinear Galerkin method, further numerical computations are performed to reveal the vibration response and average output power of the BEMH under periodic excitations. Amongst other factors, structure radius ratio and membrane thickness affect more significantly on the dynamic response of the system in the vibration of a circular membrane with a co-axis mass [53]. This paper mainly concerns the changing rules of the corresponding amplitude and the output power with the variation of these two parameters. When the parameters of the center magnet are determined, the radius ratio a/b of the system can be adjusted by changing the radius of the membrane structure. Figures 16 and 17 show that under the same excitation condition, the amplitude and average output power of the BEMH decrease with the increase of the radius ratio. Especially, when radius ratios are 0.05 to 0.1, the amplitude and average output power of the BEMH change significantly. Then the change is relatively slower when the radius ratio is greater than 0.1. In addition, with the radius ratio of 0.05, the system has a maximum error of 30.7% in amplitude and 52.9% in average power through comparing the calculated results of the two models. Furthermore, the resulted error decreases with the increase of the radius ratio.     Additionally, the thickness of the membrane also affects the vibration response and the average output power of the BEMH. As shown in Figures 18 and 19, both the vibration amplitude and the average output power will decrease with the increase of membrane thickness. In the change of radius ratios from 0.05 to 0.3, the analysis results have the largest decline with the membrane thickness of 0.04 mm that the amplitude falls by up to 75.3% and the average power reduces by 93.3%. In the change of membrane thickness from 0.04 mm to 0.12 mm, the analysis results have the largest decline with the radius ratio of 0.05 that the amplitude falls by up to 49.8% and the average power reduces by 74.3%. When the radius ratio is 0.3, the influence of membrane thickness on the output power of energy harvesting system is negligible. To sum up, the BEMH model acquired by nonlinear Galerkin method leads to more accurate results. Moreover, the output power of the energy harvester can be improved by decreasing the radius ratio of the structure and the membrane thickness, which provides a theoretical reference for optimizing the structure parameters.

Conclusions
In this paper, a novel architecture of BEMH is proposed to collect the mechanical energy characterized by low frequency and small amplitude environmental vibration. The nonlinear distributed parameter model is established for the large deflection transverse vibration membrane attached with a co-axial magnet mass. The low-order vibration equation is obtained via the nonlinear Galerkin method to cope up with the higher-order geometrical and magnetic nonlinearities. It provides an efficient foundation for considering the influence of higher-order nonlinearity, such as the nonlinearity caused by the geometric deformation of large-amplitude vibration or nonlinear magnetic force. Accordingly, the vibration characteristics and output power of the system under different excitation conditions can be accurately analyzed, thus providing an efficient theoretical basis to optimize the power output performance of the BEMH. The conclusions are as follows: 1. The proposed BEMH can harvest low-frequency and small-amplitude vibration energy in a broad bandwidth. It can carry out bistable vibration, thereby increasing the vibration displacement and improving the efficiency of energy acquisition.
2. By taking the higher-order components as second-order modes via the nonlinear Galerkin method, the model order is reduced, which allows the model to be solved efficiently. Moreover, the nonlinear Galerkin method gives a solution both close agreeable with that of FE result and experiment data.
3. With this efficient and accurate approach, it has been verified that both the radius ratio and the membrane thickness show a significant improvement on the output power with a lower value. The influence of the membrane thickness on output power is negligible when the radius ratio reaches a certain value. This finding provides a theoretical basis for designing the structure in order to obtain a higher power output.