Research on the Directional Adaptability of a Self-Adaptive Energy Harvester

With the continuous development of wireless sensor networks (WSNs), multi-directional energy harvesting technology has received widespread attention from scholars. In order to evaluate the performance of multi-directional energy harvesters, this paper uses a directional self-adaptive piezoelectric energy harvester (DSPEH) as an example, defines the direction of the excitation in three-dimensional space, and studies the influence of excitations on the key parameters of the DSPEH. The rolling angle and pitch angle are used to define complex excitations in three-dimensional space, and the dynamic response of the excitation changes in a single direction and multiple directions is discussed. It is noteworthy that this work presents the concept of “Energy Harvesting Workspace” to describe the working ability of a multi-directional energy harvesting system. The workspace is expressed by the excitation angle and voltage amplitude, and energy harvesting performance is evaluated by the volume-wrapping method and area-covering method. The DSPEH exhibits good directional adaptability in two-dimensional space (rolling direction); in particular, when the mass eccentricity coefficient is r = 0 mm, 100% of the workspace in two-dimensional space is obtained. The total workspace in three-dimensional space depends entirely on the energy output in the pitch direction.


Introduction
The wireless sensor network (WSN) is a revolution in the field of information sensing and harvesting, and it has introduced profound impacts on human production and life [1,2]. WSN consists of tens of thousands of sensor nodes, each of which includes a sensing module, an information processing module, a wireless communication module, a power management module, and a power supply module, where the power supply module supplies power to the other modules and plays a very important role in WSNs [3,4]. In order to provide clean and sustainable electrical energy for WSNs, researchers use various electromechanical conversion mechanisms, such as piezoelectric [5][6][7], electromagnetic [8][9][10], electrostatic [11][12][13], thermoelectric [14][15][16], and triboelectric [17][18][19] mechanisms to convert the vibration energy near the working environment of WSNs into electrical energy. Piezoelectric energy harvesting technology uses the direct piezoelectric effect of piezoelectric materials to realize the conversion of mechanical energy to electrical energy and has high electromechanical conversion efficiency, which has been widely researched by scholars [20,21].
Generally, the traditional piezoelectric energy harvester (PEH) is designed as a cantilever beam structure with a tip mass. When the external excitation direction is consistent with the normal direction of the cantilever beam, the vibration energy in the environment can be efficiently harvested; when the external excitation direction deviates from the normal direction of the cantilever beam, the harvesting efficiency of the traditional PEH will be Compared with the conventional PEH, the DSPEH has additional rotational DOFs, which can drive the rotation of the cantilever beam in order to track the direction of the external excitation. When there is an initial angle, θ0, between the normal direction of the cantilever beam and the excitation direction, the vibrating cantilever beam is in a state of force imbalance. Under the action of the eccentric gravity of the tip mass and vibrating centrifugal force, the cantilever beam begins to rotate to reconstruct the stable state. When the cantilever beam reaches a force equilibrium state, it begins to generate stable vibrations at a certain position without rotation. This adjustment process without any manual intervention or power supply is called a self-adaptive rotation mechanism, which enables the cantilever beam to achieve multi-directional energy harvesting.

Definition of the Three-Dimensional Excitation
As shown in Figure 2, in order to describe the direction of the external excitation in coordinate system x″-y″-z″, we use the RPY angle representation method to define three angles: θ, φ, and γ. Here, rolling angle θ represents the angle of the external excitation rotating along the z″ axis; pitch angle φ represents the angle of the external excitation rotating along y″ axis; yaw angle γ represents the angle of the external excitation rotating along x″ axis. The three-dimensional excitation direction is determined by (θ, φ, and γ). The direction of the external excitation can be described as follows: the external excitation first rotates θ angle around the z″ axis; then, it rotates φ angle around the y″ axis and finally rotates γ angle around the x″ axis, as shown in Figure 3. The corresponding rotation transformation formula is shown as follows:  Compared with the conventional PEH, the DSPEH has additional rotational DOFs, which can drive the rotation of the cantilever beam in order to track the direction of the external excitation. When there is an initial angle, θ 0 , between the normal direction of the cantilever beam and the excitation direction, the vibrating cantilever beam is in a state of force imbalance. Under the action of the eccentric gravity of the tip mass and vibrating centrifugal force, the cantilever beam begins to rotate to reconstruct the stable state. When the cantilever beam reaches a force equilibrium state, it begins to generate stable vibrations at a certain position without rotation. This adjustment process without any manual intervention or power supply is called a self-adaptive rotation mechanism, which enables the cantilever beam to achieve multi-directional energy harvesting.

Definition of the Three-Dimensional Excitation
As shown in Figure 2, in order to describe the direction of the external excitation in coordinate system x -y -z , we use the RPY angle representation method to define three angles: θ, ϕ, and γ. Here, rolling angle θ represents the angle of the external excitation rotating along the z axis; pitch angle ϕ represents the angle of the external excitation rotating along y axis; yaw angle γ represents the angle of the external excitation rotating along x axis. The three-dimensional excitation direction is determined by (θ, ϕ, and γ). Compared with the conventional PEH, the DSPEH has additional rotational DOFs, which can drive the rotation of the cantilever beam in order to track the direction of the external excitation. When there is an initial angle, θ0, between the normal direction of the cantilever beam and the excitation direction, the vibrating cantilever beam is in a state of force imbalance. Under the action of the eccentric gravity of the tip mass and vibrating centrifugal force, the cantilever beam begins to rotate to reconstruct the stable state. When the cantilever beam reaches a force equilibrium state, it begins to generate stable vibrations at a certain position without rotation. This adjustment process without any manual intervention or power supply is called a self-adaptive rotation mechanism, which enables the cantilever beam to achieve multi-directional energy harvesting.

Definition of the Three-Dimensional Excitation
As shown in Figure 2, in order to describe the direction of the external excitation in coordinate system x″-y″-z″, we use the RPY angle representation method to define three angles: θ, φ, and γ. Here, rolling angle θ represents the angle of the external excitation rotating along the z″ axis; pitch angle φ represents the angle of the external excitation rotating along y″ axis; yaw angle γ represents the angle of the external excitation rotating along x″ axis. The three-dimensional excitation direction is determined by (θ, φ, and γ). The direction of the external excitation can be described as follows: the external excitation first rotates θ angle around the z″ axis; then, it rotates φ angle around the y″ axis and finally rotates γ angle around the x″ axis, as shown in Figure 3. The corresponding rotation transformation formula is shown as follows:  The direction of the external excitation can be described as follows: the external excitation first rotates θ angle around the z axis; then, it rotates ϕ angle around the y axis and finally rotates γ angle around the x axis, as shown in Figure 3. The corresponding rotation transformation formula is shown as follows: where cϕ = cosϕ and sϕ = sinϕ. where cφ = cosφ and sφ = sinφ. Compared to the research in reference [49], which only discusses the rolling angle along the z" axis by assuming an initial angle, θ0, between the normal direction of the cantilever beam and the excitation direction, this paper adds two angles (φ and γ) to describe the direction of the external excitation, as shown in Figure 4. The included angles, θ and φ, are achieved by adjusting the installation angle of the exciter, and the included angle, γ, is achieved by adjusting the installation angle of the energy harvester. It should be noted that since the energy harvester is always installed in a certain position, the directional changes with respect to the external excitation in the x" axis will not affect the working state of the energy harvester. Therefore, only included angles θ and φ can be used to determine the excitation direction.

The Variation Law of System Parameters under Three-Dimensional Excitation
When the external excitation changes in the three-dimensional space, the external excitation and tip mass will exert an axial load on the cantilever beam under the action of the included angle, thereby affecting the natural frequency of the cantilever beam. Additionally, the eccentric gravity of the tip mass will also correspondingly increase the resistance moment on the rotation's DOF. Therefore, it is necessary to analyze the parameter changes with respect to the cantilever beam under complex excitations in order to correct the relevant parameters of the differential motion equation.
According to Equation (1), based on the coordinate transformation matrix, excitation force F acting on the cantilever beam can be decomposed into Compared to the research in reference [49], which only discusses the rolling angle along the z" axis by assuming an initial angle, θ 0 , between the normal direction of the cantilever beam and the excitation direction, this paper adds two angles (ϕ and γ) to describe the direction of the external excitation, as shown in Figure 4. The included angles, θ and ϕ, are achieved by adjusting the installation angle of the exciter, and the included angle, γ, is achieved by adjusting the installation angle of the energy harvester. It should be noted that since the energy harvester is always installed in a certain position, the directional changes with respect to the external excitation in the x" axis will not affect the working state of the energy harvester. Therefore, only included angles θ and ϕ can be used to determine the excitation direction. where cφ = cosφ and sφ = sinφ. Compared to the research in reference [49], which only discusses the rolling angle along the z" axis by assuming an initial angle, θ0, between the normal direction of the cantilever beam and the excitation direction, this paper adds two angles (φ and γ) to describe the direction of the external excitation, as shown in Figure 4. The included angles, θ and φ, are achieved by adjusting the installation angle of the exciter, and the included angle, γ, is achieved by adjusting the installation angle of the energy harvester. It should be noted that since the energy harvester is always installed in a certain position, the directional changes with respect to the external excitation in the x" axis will not affect the working state of the energy harvester. Therefore, only included angles θ and φ can be used to determine the excitation direction.

The Variation Law of System Parameters under Three-Dimensional Excitation
When the external excitation changes in the three-dimensional space, the external excitation and tip mass will exert an axial load on the cantilever beam under the action of the included angle, thereby affecting the natural frequency of the cantilever beam. Additionally, the eccentric gravity of the tip mass will also correspondingly increase the resistance moment on the rotation's DOF. Therefore, it is necessary to analyze the parameter changes with respect to the cantilever beam under complex excitations in order to correct the relevant parameters of the differential motion equation.
According to Equation (1), based on the coordinate transformation matrix, excitation force F acting on the cantilever beam can be decomposed into

The Variation Law of System Parameters under Three-Dimensional Excitation
When the external excitation changes in the three-dimensional space, the external excitation and tip mass will exert an axial load on the cantilever beam under the action of the included angle, thereby affecting the natural frequency of the cantilever beam. Additionally, the eccentric gravity of the tip mass will also correspondingly increase the resistance moment on the rotation's DOF. Therefore, it is necessary to analyze the parameter changes with respect to the cantilever beam under complex excitations in order to correct the relevant parameters of the differential motion equation.
According to Equation (1), based on the coordinate transformation matrix, excitation force F acting on the cantilever beam can be decomposed into where F x represents the component of the excitation in the normal direction of the cantilever beam, which can cause the vibration of the cantilever beam; F y represents the component of the excitation in the tangent direction of the cantilever beam, which can cause the rotation of the cantilever beam; F z represents the component of the excitation along the z axis of the cantilever beam, which generates a pre-tension or pre-pressure on the cantilever beam. The DSPEH in this paper adopts a cantilever beam structure with a tip mass. Although it has an additional rotation DOF, it still belongs to the ordinary cantilever beam's structure. Therefore, this paper adopts the same method as ordinary cantilever beams to establish a distributed parameter model, solve the changes in the natural frequency of the system when the excitation direction changes, and optimize the lumped parameter model under harmonic excitations in order to more accurately describe the motion state of the system.
Due to the large length-to-thickness ratio of the cantilever beam and the fact that the cantilever beam has rotational DOF in DSPEH, the influence of the shear deformation and rotational inertia of the cantilever beam on the system can be ignored. Therefore, this configuration conforms to the Euler-Bernoulli beam theory. Here, it is assumed that the length of the cantilever beam is L, the elastic modulus is E, the cross-section of the cantilever beam is equal, the moment of inertia of the cross-section is I, and the mass of the cantilever beam per unit length is ρ. The lateral vibration displacement of any point on the cantilever beam at time t is represented by function w(x, t). When the pre-tension or pre-pressure applied to the cantilever beam is P (P = F z ), any microelement on the cantilever beam can be taken and its stress is analyzed, as shown in Figure 5.
where Fx represents the component of the excitation in the normal direction of the cantilever beam, which can cause the vibration of the cantilever beam; Fy represents the component of the excitation in the tangent direction of the cantilever beam, which can cause the rotation of the cantilever beam; Fz represents the component of the excitation along the z″ axis of the cantilever beam, which generates a pre-tension or pre-pressure on the cantilever beam.

Influence of the Three-Dimensional Excitation on the Natural Frequency of the Cantilever Beam
The DSPEH in this paper adopts a cantilever beam structure with a tip mass. Although it has an additional rotation DOF, it still belongs to the ordinary cantilever beam's structure. Therefore, this paper adopts the same method as ordinary cantilever beams to establish a distributed parameter model, solve the changes in the natural frequency of the system when the excitation direction changes, and optimize the lumped parameter model under harmonic excitations in order to more accurately describe the motion state of the system.
Due to the large length-to-thickness ratio of the cantilever beam and the fact that the cantilever beam has rotational DOF in DSPEH, the influence of the shear deformation and rotational inertia of the cantilever beam on the system can be ignored. Therefore, this configuration conforms to the Euler-Bernoulli beam theory. Here, it is assumed that the length of the cantilever beam is L, the elastic modulus is E, the cross-section of the cantilever beam is equal, the moment of inertia of the cross-section is I, and the mass of the cantilever beam per unit length is ρ. The lateral vibration displacement of any point on the cantilever beam at time t is represented by function w(x, t). When the pre-tension or pre-pressure applied to the cantilever beam is P (P = Fz), any microelement on the cantilever beam can be taken and its stress is analyzed, as shown in Figure 5. According to the moment balance and force balance equations of the microelement, we can obtain Ignoring the second-order small quantity in the first equation of Equation (3), the partial derivative of x at both ends of the equation is obtained. According to the moment balance and force balance equations of the microelement, we can obtain Ignoring the second-order small quantity in the first equation of Equation (3), the partial derivative of x at both ends of the equation is obtained.
Accordingly, the second equation in Equations (3) and (4) can be changed into Regardless of geometric nonlinearity, it can be approximately considered that Then, Equation (6) is substituted into Equation (5).
The time and space functions in Equation (7) are separated and solved by the method of separating variables, assuming that Equation (8) is then substituted into Equation (7).
Equation (9) can be written as follows.
The left side of Equation (10) depends only on x, and the right side depends only on t. Since x and t are independent variables, the standard parameters in the method of separating variables stipulate that both sides of Equation (10) must be equal to the same constant υ.
Equation (11) can be decomposed into two independent equations.
The general solution of Equation (12) can be expressed as Substituting Equation (14) into Equation (12) yields The four solutions of the above equation are where The four solutions are substituted into Equation (14), and they are replaced with trigonometric and hyperbolic functions, yielding where A, B, C, and D; and α 1 and α 2 are determined by boundary conditions. Equation (13) is a single DOF vibration differential equation, and the general solution can be expressed as where H and ψ are determined by initial conditions. For a cantilever beam clamped at x = 0 and free at x = L with a tip mass, the boundary conditions of the vibration equation are where Equations (19) and (20) are the geometric boundary conditions at x = 0, and Equations (21) and (22) are the geometric boundary conditions at x = L. Here, m t is the mass of the tip mass, and I t is the mass inertia moment of the tip mass at x = L. Substituting Equations (8), (17), and (18) into Equations (19) and (20) yields Substituting Equations (8), (17), (18), and (23) into Equations (21) and (22) where In order to obtain the non-zero solution of the above equation system, the coefficient matrix of the equation must be singular; that is, the following is obtained.
According to the above equation, the relationship between the natural frequency of cantilever beam ω and axial force P can be obtained via MATLAB numerical simulations using the Newton Iteration, as shown in Figure 6. From the figure, it can be observed that when the axial force changes within the range of −1.0~+1.0 N, the natural frequency is approximately linearly related to the axial force. When the cantilever beam is stretched, the natural frequency increases with an increase in axial force; when the cantilever beam is compressed, the natural frequency decreases with an increase in axial force.

EI  
In order to obtain the non-zero solution of the above equation system, the coefficient matrix of the equation must be singular; that is, the following is obtained. 11 22 12 21 0 a a a a   (25) According to the above equation, the relationship between the natural frequency of cantilever beam ω and axial force P can be obtained via MATLAB numerical simulations using the Newton Iteration, as shown in Figure 6. From the figure, it can be observed that when the axial force changes within the range of −1.0~+1.0 N, the natural frequency is approximately linearly related to the axial force. When the cantilever beam is stretched, the natural frequency increases with an increase in axial force; when the cantilever beam is compressed, the natural frequency decreases with an increase in axial force. The natural frequency of the cantilever beam is related to the axial force. According to Equation (2), the axial force is related to the magnitude of external excitations and pitch angle φ, and it is independent of rolling angle θ. Figure 7 shows the variation of the natural frequency of the cantilever beam with respect to pitch angle φ when the excitation force is F = −0.05 N. From the figure, it can be observed that within the range of 0~+90°, the natural frequency of the cantilever beam decreases as pitch angle φ increases; within the range of −90~0°, the natural frequency increases as pitch angle φ increases. The natural frequency of the cantilever beam is related to the axial force. According to Equation (2), the axial force is related to the magnitude of external excitations and pitch angle ϕ, and it is independent of rolling angle θ. Figure 7 shows the variation of the natural frequency of the cantilever beam with respect to pitch angle ϕ when the excitation force is F = −0.05 N. From the figure, it can be observed that within the range of 0~+90 • , the natural frequency of the cantilever beam decreases as pitch angle ϕ increases; within the range of −90~0 • , the natural frequency increases as pitch angle ϕ increases.

Influence of the Three-Dimensional Excitation on the Rotational Resistance Moment
When external excitation changes in the rolling and pitch directions, the eccentric gravity of the tip mass will generate a resistance moment, TC, that acts on the z" axis, which will affect the steady-state vibration position and output voltage amplitude of the harvester. This resistance moment can also be obtained using the coordinate transformation matrix.

cos sin
The variation of resistance moment TC with respect to the excitation direction is shown in Figure 8. From the figure, it can be seen that the resistance moment increases with an increase in rolling angle θ and decreases with an increase in pitch angle φ within the range of −90~+90°. However, when rolling angle θ is zero, the change in pitch angle does not affect the change in resistance moment, so the key factor affecting the resistance moment is rolling angle θ.

Influence of the Three-Dimensional Excitation on the Rotational Resistance Moment
When external excitation changes in the rolling and pitch directions, the eccentric gravity of the tip mass will generate a resistance moment, T C , that acts on the z" axis, which will affect the steady-state vibration position and output voltage amplitude of the harvester. This resistance moment can also be obtained using the coordinate transformation matrix.
The variation of resistance moment T C with respect to the excitation direction is shown in Figure 8. From the figure, it can be seen that the resistance moment increases with an increase in rolling angle θ and decreases with an increase in pitch angle ϕ within the range of −90~+90 • . However, when rolling angle θ is zero, the change in pitch angle does not affect the change in resistance moment, so the key factor affecting the resistance moment is rolling angle θ.
vester. This resistance moment can also be obtained using the coordinate transformation matrix.
cos sin C T mgr    (26) The variation of resistance moment TC with respect to the excitation direction is shown in Figure 8. From the figure, it can be seen that the resistance moment increases with an increase in rolling angle θ and decreases with an increase in pitch angle φ within the range of −90~+90°. However, when rolling angle θ is zero, the change in pitch angle does not affect the change in resistance moment, so the key factor affecting the resistance moment is rolling angle θ. Based on the above analysis, the dynamic equation of the DSPEH can be expressed as Based on the above analysis, the dynamic equation of the DSPEH can be expressed as where m is the equivalent mass; c 1 is the viscous damping coefficient in the x direction; k represents the equivalent stiffness; g is the gravitational acceleration; r is the mass eccentricity coefficient, which represents the distance of the equivalent mass from the center of the cantilever beam, as shown in Figure 1; α represents the equivalent electromechanical coupling coefficient; V represents the output voltage of the piezoelectric patch; F represents the external excitation force; ω' represents the frequency of the external excitation; c 2 is the viscous damping coefficient in θ direction; C P represents the equivalent capacitance of the piezoelectric patch; R represents the equivalent resistance of the energy harvesting circuit.

Prototype and Experimental Setup
As shown in Figure 9, the experimental device includes a signal generator, a power amplifier, an electromagnetic exciter, an oscilloscope, an angle sensor, and a DC power supply. The signal generator is used to send out a sine signal, and the power amplifier amplifies and transmits the sine signal to the electromagnetic exciter. The electromagnetic exciter generates a three-dimensional force, F, acting on the DSPEH, which is installed on the electromagnetic exciter and outputs the voltage signals of the piezoelectric patch to the oscilloscope. The rotation of the cantilever beam is monitored by the angle sensor and transmitted to the oscilloscope, and the DC power supply is used to supply power to the angle sensor. During the experiment, when rolling angle θ or pitch angle ϕ changes separately, adjustment nut B is used on the exciter to adjust the parameters to the required angle, and assembly fixture C is used to ensure the horizontal installation position of the DSPEH. When rolling angle θ and pitch angle ϕ change simultaneously, the included angle is adjusted by adjusting fixture A and adjusting nut B, respectively; moreover, fixture C is used to ensure the horizontal installation position of the DSPEH. Rolling angle θ and pitch angle ϕ can be adjusted within the range of ±90 • . The parameters of the DSPEH are shown in Table 1.
arately, adjustment nut B is used on the exciter to adjust the parameters to the required angle, and assembly fixture C is used to ensure the horizontal installation position of the DSPEH. When rolling angle θ and pitch angle φ change simultaneously, the included angle is adjusted by adjusting fixture A and adjusting nut B, respectively; moreover, fixture C is used to ensure the horizontal installation position of the DSPEH. Rolling angle θ and pitch angle φ can be adjusted within the range of ±90°. The parameters of the DSPEH are shown in Table 1.   Piezoelectric materials have a critical impact on electromechanical energy conversion. The most common piezoelectric materials include PZT, PVDF, MFC, and so on. Although PZT is widely used for generating electrical energy, its inherent brittleness makes it safe to use only within a certain stress-strain range. PVDF is a piezoelectric polymer with better flexibility than PZT, and it can withstand greater strain and is conducive to generating more electrical energy. Although MFC has good comprehensive performance, the current market price of MFC is tens or hundreds of times higher than that of PZT and PVDF. Therefore, considering the large strain of the cantilever beam and the cost of the experimental system, PVDF is chosen as the piezoelectric material in this paper. Figure 10 shows the output voltage and cantilever beam's position versus time when the rolling angle is θ = 60 • . The red line represents the position of the cantilever beam, and the black line represents the output voltage of the piezoelectric patch. For the convenience of understanding, the position of the cantilever beam is represented by the orientation angle between its normal direction and the excitation direction; that is, when the normal direction of the cantilever beam is consistent with the excitation's direction, the orientation angle of the cantilever beam position is 0 • . It can be observed from the figure that the orientation angle of the cantilever beam position gradually decreases to the steady-state vibration position, and the output voltage amplitude of the cantilever beam increases with a decrease in the orientation angle, ultimately reaching a stable voltage amplitude. At the initial stage, the orientation angle of the cantilever beam position is 60 • , and at the final steady-state vibration, the orientation angle is approximately 35 • . The cantilever beam is not located at the ideal 0 • position during steady-state vibrations due to the effect of the resistance moment generated by the eccentric gravity of the tip mass on the energy harvester.  Figure 11 shows the steady-state vibration positions and voltage output amplitudes of the cantilever beam under different rolling angles when the mass eccentricity coefficient is r = 3 mm. It can be observed from the figure that different rolling angles correspond to different steady-state vibration positions, resulting in different voltage output amplitudes. On the one hand, eccentric gravity drives the cantilever beam to rotate under external excitation; on the other hand, eccentric gravity generates a resistance moment that prevents the cantilever beam from rotating. When the rolling angle of the excitation is small (such as θ = 10° or θ = 20°), the amplitude response of the cantilever beam is large, the rotational torque generated by the centrifugal force of vibrations can overcome the resistance moment of eccentric gravity, driving the cantilever beam to rotate to a position that is roughly consistent with the excitation's direction. When the rolling angle is large (such as θ = 60°), the amplitude response of the cantilever beam is very small and can only drive the cantilever beam to rotate to a certain position. Figure 10 shows that the cantilever beam produces balance near an orientation angle of 35°, and as the rolling angle further increases (such as θ = 80°), the equilibrium position hardly changes as the resistance moment and vibration's centrifugal force reach an equilibrium state.

The Output Characteristics of the DSPEH When the Excitation Changes in the Rolling Direction
According to Figure 11, it can be observed that due to the presence of eccentric gravity, when the rolling angle is large, the cantilever beam cannot rotate to the ideal 0° position for large-amplitude vibrations. The steady-state vibration position of the cantilever beam is related to the magnitude of mass eccentricity coefficient r. Figure 12 shows the steady-state vibration positions and voltage output amplitudes of the cantilever beam under different rolling angles when the mass eccentricity coefficient is r = 0 mm. It can be observed from the figure that due to the absence of an additional resistance moment, the cantilever beam can rotate to the 0° position for large-amplitude vibrations regardless of the change in external excitation direction. It should be noted that when the rolling angle is θ = 90°, theoretically, the cantilever beam cannot rotate to the 0° position, but in practice, due to unbalanced mass or other reasons caused by manufacturing and assembly processes, the cantilever beam can also rotate to the 0° position, so there is also a large voltage output from the cantilever beam in this excitation direction. Figures 13 and 14 show the steady-state output voltage amplitude and steady-state vibration position under different mass eccentricity coefficients: The steady-state output voltage amplitude is related to the position of the cantilever beam. When the orientation angle of the steady-state vibration position is small, the output voltage amplitude is larger, and vice versa, the output voltage  Figure 11 shows the steady-state vibration positions and voltage output amplitudes of the cantilever beam under different rolling angles when the mass eccentricity coefficient is r = 3 mm. It can be observed from the figure that different rolling angles correspond to different steady-state vibration positions, resulting in different voltage output amplitudes. On the one hand, eccentric gravity drives the cantilever beam to rotate under external excitation; on the other hand, eccentric gravity generates a resistance moment that prevents the cantilever beam from rotating. When the rolling angle of the excitation is small (such as θ = 10 • or θ = 20 • ), the amplitude response of the cantilever beam is large, the rotational torque generated by the centrifugal force of vibrations can overcome the resistance moment of eccentric gravity, driving the cantilever beam to rotate to a position that is roughly consistent with the excitation's direction. When the rolling angle is large (such as θ = 60 • ), the amplitude response of the cantilever beam is very small and can only drive the cantilever beam to rotate to a certain position. Figure 10 shows that the cantilever beam produces balance near an orientation angle of 35 • , and as the rolling angle further increases (such as θ = 80 • ), the equilibrium position hardly changes as the resistance moment and vibration's centrifugal force reach an equilibrium state.
According to Figure 11, it can be observed that due to the presence of eccentric gravity, when the rolling angle is large, the cantilever beam cannot rotate to the ideal 0 • position for large-amplitude vibrations. The steady-state vibration position of the cantilever beam is related to the magnitude of mass eccentricity coefficient r. Figure 12 shows the steady-state vibration positions and voltage output amplitudes of the cantilever beam under different rolling angles when the mass eccentricity coefficient is r = 0 mm. It can be observed from the figure that due to the absence of an additional resistance moment, the cantilever beam can rotate to the 0 • position for large-amplitude vibrations regardless of the change in external excitation direction. It should be noted that when the rolling angle is θ = 90 • , theoretically, the cantilever beam cannot rotate to the 0 • position, but in practice, due to unbalanced mass or other reasons caused by manufacturing and assembly processes, the cantilever beam can also rotate to the 0 • position, so there is also a large voltage output from the cantilever beam in this excitation direction. Figures 13 and 14 show the steady-state output voltage amplitude and steady-state vibration position under different mass eccentricity coefficients: The steady-state output voltage amplitude is related to the position of the cantilever beam. When the orientation angle of the steady-state vibration position is small, the output voltage amplitude is larger, and vice versa, the output voltage amplitude is smaller. From the figure, it can be observed that when the mass eccentricity coefficient changes from large to small, the output voltage amplitude of the cantilever beam increases, and the wrapped area within the excitation range of 0 •~+ 90 • increases. When the mass eccentricity coefficient is r = 0 mm, the maximum working area indicates that the energy harvester has the best directional adaptability in the two-dimensional space.   Figure 15 shows the sweep frequency curve for different rolling angles when the mass eccentricity coefficient is r = 3 mm. When the rolling angle changes, due to the directional self-adaptive characteristics, the cantilever beam can always rotate to a certain position where the resistance moment generated by the eccentric gravity and the rotation       Figure 15 shows the sweep frequency curve for different rolling angles when the mass eccentricity coefficient is r = 3 mm. When the rolling angle changes, due to the directional self-adaptive characteristics, the cantilever beam can always rotate to a certain position where the resistance moment generated by the eccentric gravity and the rotation moment generated by the vibration centrifugal force reach a balanced torque state. Different rolling angles correspond to different steady-state vibration positions, which correspond to different output voltage amplitudes. However, the natural frequency of the cantilever beam does not change, which corresponds to Equation (2). In summary, under actual installation angles, when the external excitation changes in the rolling direction, due to the self-adaptive feature of the cantilever beam, the energy harvesting system can automatically adjust to a certain position for stable vibration. Due to the existence of eccentric gravity, changes in the rolling angle will generate a rotational resistance moment, which affects the steady-state vibration position and voltage output amplitude of the system, but it does not affect the natural frequency of the cantilever beam.

The Output Characteristics of the System When the Excitation Changes in the Pitch Direction
This section analyzes the output characteristics of the energy harvesting system when the excitation direction changes in the pitch direction. Due to the pitch angle between the normal direction of the cantilever beam and the external excitation direction, the excitation acting on the cantilever beam is decomposed into two parts: one part is applied to the normal direction of the cantilever beam, driving the cantilever beam to vibrate; the other part is applied to the axis direction of the cantilever beam, generating an axial force that affects the natural frequency of the cantilever beam.   Figure 15 shows the sweep frequency curve for different rolling angles when the mass eccentricity coefficient is r = 3 mm. When the rolling angle changes, due to the directional self-adaptive characteristics, the cantilever beam can always rotate to a certain position where the resistance moment generated by the eccentric gravity and the rotation moment generated by the vibration centrifugal force reach a balanced torque state. Different rolling angles correspond to different steady-state vibration positions, which correspond to different output voltage amplitudes. However, the natural frequency of the cantilever beam does not change, which corresponds to Equation (2). In summary, under actual installation angles, when the external excitation changes in the rolling direction, due to the self-adaptive feature of the cantilever beam, the energy harvesting system can automatically adjust to a certain position for stable vibration. Due to the existence of eccentric gravity, changes in the rolling angle will generate a rotational resistance moment, which affects the steady-state vibration position and voltage output amplitude of the system, but it does not affect the natural frequency of the cantilever beam.

The Output Characteristics of the System When the Excitation Changes in the Pitch Direction
This section analyzes the output characteristics of the energy harvesting system when the excitation direction changes in the pitch direction. Due to the pitch angle between the normal direction of the cantilever beam and the external excitation direction, the excitation acting on the cantilever beam is decomposed into two parts: one part is applied to the normal direction of the cantilever beam, driving the cantilever beam to vibrate; the other part is applied to the axis direction of the cantilever beam, generating an axial force that affects the natural frequency of the cantilever beam. In summary, under actual installation angles, when the external excitation changes in the rolling direction, due to the self-adaptive feature of the cantilever beam, the energy harvesting system can automatically adjust to a certain position for stable vibration. Due to the existence of eccentric gravity, changes in the rolling angle will generate a rotational resistance moment, which affects the steady-state vibration position and voltage output amplitude of the system, but it does not affect the natural frequency of the cantilever beam.

The Output Characteristics of the System When the Excitation Changes in the Pitch Direction
This section analyzes the output characteristics of the energy harvesting system when the excitation direction changes in the pitch direction. Due to the pitch angle between the normal direction of the cantilever beam and the external excitation direction, the excitation acting on the cantilever beam is decomposed into two parts: one part is applied to the normal direction of the cantilever beam, driving the cantilever beam to vibrate; the other part is applied to the axis direction of the cantilever beam, generating an axial force that affects the natural frequency of the cantilever beam. Figure 16 shows the steady-state vibration position and output voltage of the cantilever beam at different pitch angles. The red curve in the figure represents the position of the cantilever beam in the rolling direction. From the figure, it can be observed that as the pitch angle increases, the steady-state position of the cantilever beam in the rolling direction remains unchanged, but the output voltage amplitude decreases as the pitch angle increases. According to Equation (2), it can be observed that as the pitch angle increases, the component force of the external excitation in the normal direction of the cantilever beam will gradually decrease in a manner exhibiting a trigonometric function; that is, the excitation intensity of the cantilever beam will decrease, so its output voltage will also decrease. Similar results are obtained during the experimental frequency sweep process, as shown in Figure 17. From this perspective, the adaptability of the DSPEH in the pitch direction is poor, which is the same as a conventional PEH.
the cantilever beam in the rolling direction. From the figure, it can be observed that as the pitch angle increases, the steady-state position of the cantilever beam in the rolling direction remains unchanged, but the output voltage amplitude decreases as the pitch angle increases. According to Equation (2), it can be observed that as the pitch angle increases, the component force of the external excitation in the normal direction of the cantilever beam will gradually decrease in a manner exhibiting a trigonometric function; that is, the excitation intensity of the cantilever beam will decrease, so its output voltage will also decrease. Similar results are obtained during the experimental frequency sweep process, as shown in Figure 17. From this perspective, the adaptability of the DSPEH in the pitch direction is poor, which is the same as a conventional PEH.  Figure 18 shows the sweep frequency curves of the energy harvesting system under different pitch angles. It can be observed from the figure that as the pitch angle changes, the natural frequency shifts to the right and the output voltage peak decreases. The natu-  Figure 18 shows the sweep frequency curves of the energy harvesting system under different pitch angles. It can be observed from the figure that as the pitch angle changes, the natural frequency shifts to the right and the output voltage peak decreases. The natural frequency shift in the experiment indicates that the additional axial load caused by the pitch angle has an impact on the cantilever beam. Figure 19 shows the natural frequency curves of the cantilever beam under different pitch angles. The solid line represents the experimental results, and the dashed line represents the simulation data. Due to the small excitation acceleration (0.7 g) during the experimental process, the additional axial force generated by the excitation component is also very small, so the natural frequency changes little. In addition, there are uncontrollable factors and measurement errors during the experimental process, and the natural frequencies of the experiment and simulation cannot fully correspond with one another. However, the trend of the two is the same: in the range of 0~+90 • , the natural frequency of the cantilever beam decreases with an increase in the pitch angle; within the range of −90~0 • , the natural frequency increases with an increase in pitch angle.  Figure 18 shows the sweep frequency curves of the energy harvesting system u different pitch angles. It can be observed from the figure that as the pitch angle cha the natural frequency shifts to the right and the output voltage peak decreases. The n ral frequency shift in the experiment indicates that the additional axial load caused b pitch angle has an impact on the cantilever beam. Figure 19 shows the natural frequ curves of the cantilever beam under different pitch angles. The solid line represent experimental results, and the dashed line represents the simulation data. Due to the s excitation acceleration (0.7 g) during the experimental process, the additional axial generated by the excitation component is also very small, so the natural frequency cha little. In addition, there are uncontrollable factors and measurement errors during th perimental process, and the natural frequencies of the experiment and simulation ca fully correspond with one another. However, the trend of the two is the same: in the r of 0~+90°, the natural frequency of the cantilever beam decreases with an increase i pitch angle; within the range of −90~0°, the natural frequency increases with an inc in pitch angle.   For different initial rolling angles, the DSPEH can track the external excitation direction by automatically calibrating the position of the cantilever beam, ultimately reaching the same position and stabilizing the vibration, indicating that the DSPEH also has the feature of self-adaptive energy harvesting in the rolling direction under different pitch directions. Although the final steady-state position for different pitch directions is the same, the output voltage is different, because the pitch angle causes different excitation components acting on the cantilever beam, which is the same as the analysis in the previous section. Figure 21 shows the steady-state output voltage amplitude of the cantilever beam at different rolling angles when the pitch angles are 0°, 30°, and 60°. From the figure, it can be observed that due to the self-adaptive feature of the energy harvester in the rolling direction, the output characteristics of the harvester are the same at different rolling angles, while the output voltage varies at different pitch angles. As the pitch angle increases, the area covered by the voltage amplitude in the graph decreases. For different initial rolling angles, the DSPEH can track the external excitation direction by automatically calibrating the position of the cantilever beam, ultimately reaching the same position and stabilizing the vibration, indicating that the DSPEH also has the feature of self-adaptive energy harvesting in the rolling direction under different pitch directions. Although the final steady-state position for different pitch directions is the same, the output voltage is different, because the pitch angle causes different excitation components acting on the cantilever beam, which is the same as the analysis in the previous section. Figure 21 shows the steady-state output voltage amplitude of the cantilever beam at different rolling angles when the pitch angles are 0 • , 30 • , and 60 • . From the figure, it can be observed that due to the self-adaptive feature of the energy harvester in the rolling direction, the output characteristics of the harvester are the same at different rolling angles, while the output voltage varies at different pitch angles. As the pitch angle increases, the area covered by the voltage amplitude in the graph decreases.  Figure 22 shows the steady-state output voltage amplitude of the cantilever be different pitch angles when the rolling angles are 0°, 30°, and 60°. As observed from figure, due to the self-adaptive feature of the energy harvester in the rolling directio output characteristics of the harvester are the same at different rolling angles. Acco to the above analysis, the main factor affecting the multi-directional performance o energy harvester is the pitch angle, and the directional self-adaptive feature allow harvester to fully adapt to the rolling angle.   Figure 22 shows the steady-state output voltage amplitude of the cantilever beam at different pitch angles when the rolling angles are 0 • , 30 • , and 60 • . As observed from the figure, due to the self-adaptive feature of the energy harvester in the rolling direction, the output characteristics of the harvester are the same at different rolling angles. According to the above analysis, the main factor affecting the multi-directional performance of the energy harvester is the pitch angle, and the directional self-adaptive feature allows the harvester to fully adapt to the rolling angle. different pitch angles when the rolling angles are 0°, 30°, and 60°. As observed from the figure, due to the self-adaptive feature of the energy harvester in the rolling direction, the output characteristics of the harvester are the same at different rolling angles. According to the above analysis, the main factor affecting the multi-directional performance of the energy harvester is the pitch angle, and the directional self-adaptive feature allows the harvester to fully adapt to the rolling angle.

Performance Analysis and Evaluation of a Multi-Directional Energy Harvester
In summary, when the excitation direction changes in three-dimensional space, the external excitation direction has a significant impact on the output characteristics of the energy harvesting system. In order to evaluate the output performance of energy harvest ers under complex excitations, the concept of "Energy Harvesting Workspace" is put for ward to describe the ability of energy harvesting in multiple directions. The energy har vesting workspace is represented by the angles of external excitation and the amplitude of output voltage. In the coordinate system, the horizontal axis represents the angles o the external excitation, and the vertical axis represents the amplitude of the output volt age. The energy harvesting workspaces of the DSPEH under different mass eccentricity coefficients are shown in Figures 23 and 24.
From Figure 23, it can be observed that when the mass eccentricity coefficient is r = 0 mm, the three-dimensional workspace of the DSPEH is arched, and its boundary is limited by the different voltage amplitudes generated by pitch angle φ. By converting angles to radians, the volume wrapped by the arch is 50.9 V, occupying 63% of the total three-di mensional space (80.9 V). When external excitation changes in the two-dimensional space the two-dimensional workspace of the DSPEH is equivalent to the two-dimensional cross section of the three-dimensional workspace. As shown on the right side of Figure 23, due to its direction self-adaptive feature, the DSPEH can always vibrate stably at a certain

Performance Analysis and Evaluation of a Multi-Directional Energy Harvester
In summary, when the excitation direction changes in three-dimensional space, the external excitation direction has a significant impact on the output characteristics of the energy harvesting system. In order to evaluate the output performance of energy harvesters under complex excitations, the concept of "Energy Harvesting Workspace" is put forward to describe the ability of energy harvesting in multiple directions. The energy harvesting workspace is represented by the angles of external excitation and the amplitude of output voltage. In the coordinate system, the horizontal axis represents the angles of the external excitation, and the vertical axis represents the amplitude of the output voltage. The energy harvesting workspaces of the DSPEH under different mass eccentricity coefficients are shown in Figures 23 and 24. position, and the output voltage amplitude in the rolling direction is consistent. When the pitch angle is φ = 0°, the output voltage's amplitude curve covers an area of 25.8 V within −90~+90°, occupying 100% of the total two-dimensional space (25.8 V). However, when the pitch angle is φ = 60 °, due to the limited excitation component acting on the cantilever beam, the area covered by the output voltage amplitude curve only accounts for 50% of the total area, which is significantly lower than that of φ = 0°. In addition, as shown on the left side of Figure 23, the area of the two-dimensional workspace obtained from the pitch angle cross-section accounts for 63% of the total area and remains unchanged, which is due to the self-adaptive feature in the rolling direction.  When the mass eccentricity coefficient is r = 3 mm, the workspace is an arched shape From Figure 23, it can be observed that when the mass eccentricity coefficient is r = 0 mm, the three-dimensional workspace of the DSPEH is arched, and its boundary is limited by the different voltage amplitudes generated by pitch angle ϕ. By converting angles to radians, the volume wrapped by the arch is 50.9 V, occupying 63% of the total three-dimensional space (80.9 V). When external excitation changes in the two-dimensional space, the two-dimensional workspace of the DSPEH is equivalent to the two-dimensional cross-section of the three-dimensional workspace. As shown on the right side of Figure 23, due to its direction self-adaptive feature, the DSPEH can always vibrate stably at a certain position, and the output voltage amplitude in the rolling direction is consistent. When the pitch angle is ϕ = 0 • , the output voltage's amplitude curve covers an area of 25.8 V within −90~+90 • , occupying 100% of the total two-dimensional space (25.8 V). However, when the pitch angle is ϕ = 60 • , due to the limited excitation component acting on the cantilever beam, the area covered by the output voltage amplitude curve only accounts for 50% of the total area, which is significantly lower than that of ϕ = 0 • . In addition, as shown on the left side of Figure 23, the area of the two-dimensional workspace obtained from the pitch angle cross-section accounts for 63% of the total area and remains unchanged, which is due to the self-adaptive feature in the rolling direction.  When the mass eccentricity coefficient is r = 3 mm, the workspace is an arched shape with a raised middle, and its boundary is affected by both pitch and rolling angles, as shown in Figure 24. The rolling angle affects the height of the arch: When the rolling angle changes from −90° to +90°, the height of the arch first increases and then decreases. The height of the arch at both ends is significantly smaller than that of r = 0 mm. This specially shaped arch is generated by mass eccentricity coefficient r, and its wrapping volume is 38.8 V, occupying 48% of the total three-dimensional space (80.9 V). This arch also correspondingly reduces the two-dimensional workspace; for example, when the pitch angle is φ = 0°, the two-dimensional space in the rolling direction accounts for 89% of the total space. When the rolling angle is θ = 60°, the two-dimensional space in the pitch direction accounts for 51% of the total space. From this perspective, the energy harvester proposed in this paper has good adaptability in the direction of the two-dimensional space (rolling direction). In particular, when When the mass eccentricity coefficient is r = 3 mm, the workspace is an arched shape with a raised middle, and its boundary is affected by both pitch and rolling angles, as shown in Figure 24. The rolling angle affects the height of the arch: When the rolling angle changes from −90 • to +90 • , the height of the arch first increases and then decreases. The height of the arch at both ends is significantly smaller than that of r = 0 mm. This specially shaped arch is generated by mass eccentricity coefficient r, and its wrapping volume is 38.8 V, occupying 48% of the total three-dimensional space (80.9 V). This arch also correspondingly reduces the two-dimensional workspace; for example, when the pitch angle is ϕ = 0 • , the two-dimensional space in the rolling direction accounts for 89% of the total space. When the rolling angle is θ = 60 • , the two-dimensional space in the pitch direction accounts for 51% of the total space. From this perspective, the energy harvester proposed in this paper has good adaptability in the direction of the two-dimensional space (rolling direction). In particular, when the mass eccentricity coefficient is r = 0 mm, it obtains all workspaces in the two-dimensional space. The total workspace in the three-dimensional direction depends entirely on the energy output of the second angle (pitch direction) as the pitch angle affects the size of the excitation component, which limits the ability of the energy harvester to work in multiple directions.

Conclusions
This paper discusses in detail the directional self-adaptation performance of the DSPEH. The rolling angle and the pitch angle are used to define complex excitations in three-dimensional space. A distributed parameter model of the energy harvesting system is established based on the moment balance and force balance equations of the microelement. The variation law of the natural frequency and rotational resistance moment of the cantilever beam when the external excitation direction changes are studied, as well as the dynamic response of the DSPEH in different excitation directions. The main conclusions are as follows. Firstly, the natural frequency of the cantilever beam is related to the pitch angle of the external excitation and is not related to the rolling angle. The resistance moment of the cantilever beam exhibits a trigonometric relationship relative to the pitch angle and the rolling angle, and it increases with an increase in the rolling angle. Secondly, the output characteristics of the DSPEH are studied when the rolling angle and pitch angle of the external excitation change individually or simultaneously. Due to the directional selfadaptive feature of the energy harvester in the rolling direction, the main factor affecting multi-directional performance is pitch angle. Thirdly, in order to describe the output performance of energy harvesting systems in multiple directions, this paper puts forward the concept of "Energy Harvesting Workspace" and establishes a performance evaluation method for multi-directional energy harvesting. The workspace is expressed using external excitation angles and output voltage amplitudes, and the output performance of energy harvesters is evaluated using volume-wrapping and area-covering methods. The DSPEH has excellent directional adaptability in the two-dimensional space (rolling direction), especially when the mass eccentricity coefficient is r = 0 mm, achieving 100% of the workspace in the two-dimensional space. The total workspace in the three-dimensional direction depends entirely on the energy harvesting efficiency of the pitch angle. The volume of the workspace package is 50.9 V, occupying 63% of the total space (80.9 V).