Research on Vibration Energy Harvester Based on Two-Dimensional Acoustic Black Hole

The wave energy focus effect of an acoustic black hole (ABH) is used for broadband vibration energy harvesting and boosts the harvested power. A vibration energy harvester based on two-dimensional ABH is proposed in this study, which consists of a rectangle plate with 2-D ABH and PZT film attached. The structure of ABH was designed and analyzed based on numerical simulation. The optimal parameters of the ABH were obtained, such as the power index, truncation thickness, cross-sectional length, and round table diameter, which were 3, 0.4 mm, 40 mm, and 24 mm, respectively. The quadratic velocity of the plate surface with ABH is up to 22.33 times that of a flat plate, and PZT film adheres to the corresponding positions of the ABH structure and plate structure, respectively. In the same condition, the average output power of a PZT with an ABH structure is higher than that of a flat plate under the same excitation-vibration condition.


Introduction
With the continuous development of Internet of Things (IoT) technology, low-power electronic products, such as wireless sensor networks (WSN) and portable wearable devices, have been widely used in recent years. Traditional battery methods are usually applied to low-power devices' power supplies. The limited service life of batteries makes lowpower electronic devices hard to maintain, especially when deployed in a remote area or harsh environment.
Various forms of energy sources widely exist in the natural environment, such as solar, wind, vibration, tidal, etc. Mechanical vibration energy is the most common energy source, which exists in production workshops, transportation such as vehicles and ships, construction such as tunnels and bridges, and human movement. A vibration energy harvester (VEH) [1][2][3] is an alternative battery power supply method, and it converts vibration energy into electric energy. VEH can be divided into three types according to their operation principle, such as piezoelectric, electrostatic, or electromagnetic energy harvesters. The piezoelectric VEH is most commonly used due to its simple structure and high electromechanical conversion efficiency. However, traditional piezoelectric VEHs with linear cantilever structures have narrow operating frequency bands and low harvesting efficiency.
Nonlinear methods can be applied in VEHs for frequency broadening and harvesting efficiency improvement [4][5][6], which consist of nonlinear structures and external nonlinear forces. A. J. Sneller [7] realized the transition from a monostable system to a bistable system by introducing an axial prestress beam in the vibration energy harvester. When the center mass increases from 0 to 6 µm, the output power of the energy harvester increases from 0.439 µW to 70.1 µW. M. Darakhshani [8] developed a bistable buckled structure by connecting the main beam and two cantilever beams with a torsion beam. The

ABH Theory
In a one-dimensional uniform dielectric beam, an acoustic black hole (ABH)'s wave resistance and energy concentration effects appear when the bending wave passes through a section that cuts the structure's thickness according to a decreasing power-law relationship. Assume that the power-law expression of the variable thickness section is: where H is the thickness of the beam, X is the distance to the tip of the power-law curve, ε is a constant coefficient, m is the power exponent under a positive rational number, and m ≥ 2. The thickness of truncation (h 1 ) cannot be zero during the actual manufacturing. The expression of the truncation thickness h 1 at the center of the actual 2-D ABH structure is: Assuming that the thickness of the beam is constant, which is H, and the length of the ABH section is L, then L and ε are: Substituting the coefficient ε into Equation (2):  The center of the actual two-dimensional ABH structure will form a circular platform with a truncated thickness of h1. The so-called two-dimensional ABH is the area formed by the one-dimensional ABH section with x = 0 as the circle's center and one rotation. Therefore, the section still follows the power-law expression. For the thin plate with changing section thickness, its bending wave equation is expressed as: When the moment of inertia and shear effect of the thin plate is not considered, the wave number k is obtained, and then the phase velocity c = ω/k of the bending wave is obtained. When the thickness of the section strictly follows the change law of the power function and the material is uniform, that is, Under ideal conditions, h1 = 0 and when m ≥ 2, as the thickness of the thin plate section decreases, the wave velocity of the flexural wave gradually decreases to 0, the wave amplitude increases slowly, and the flexural wave is concentrated at the center of the ABH, resulting in an energy focus effect.

Structural Design of 2-D ABH
A rectangular plate is designed, and the material of the plate is aluminum. As shown in Figure 2, a 2-D ABH is integrated at the position of x = 250 mm on the rectangle plate. The left boundary of the plate is fixed, and a 2 N point load with sinusoidal excitation is applied to the free right side of the plate. The center of the actual two-dimensional ABH structure will form a circular platform with a truncated thickness of h 1 . The so-called two-dimensional ABH is the area formed by the one-dimensional ABH section with x = 0 as the circle's center and one rotation. Therefore, the section still follows the power-law expression. For the thin plate with changing section thickness, its bending wave equation is expressed as: When the moment of inertia and shear effect of the thin plate is not considered, the wave number k is obtained, and then the phase velocity c = ω/k of the bending wave is obtained. When the thickness of the section strictly follows the change law of the power function and the material is uniform, that is, Under ideal conditions, h 1 = 0 and when m ≥ 2, as the thickness of the thin plate section decreases, the wave velocity of the flexural wave gradually decreases to 0, the wave amplitude increases slowly, and the flexural wave is concentrated at the center of the ABH, resulting in an energy focus effect.

Structural Design of 2-D ABH
A rectangular plate is designed, and the material of the plate is aluminum. As shown in Figure 2, a 2-D ABH is integrated at the position of x = 250 mm on the rectangle plate. The left boundary of the plate is fixed, and a 2 N point load with sinusoidal excitation is applied to the free right side of the plate.
The specific material properties and geometric dimensions are shown in Table 1. The thickness of the ABH changes according to the power function law. Since it is difficult to reduce the thickness to 0 due to the processing conditions, the truncated thickness is taken as 0.4 mm, so the central area of the ABH is a cylinder with a diameter of d and a height of 0.4 mm. Referring to Formula (5), the ABH structure parameters are shown in Table 2. The structure diagram is shown in Figure 3. The specific material properties and geometric dimensions are shown in Table 1. The thickness of the ABH changes according to the power function law. Since it is difficult to reduce the thickness to 0 due to the processing conditions, the truncated thickness is taken as 0.4 mm, so the central area of the ABH is a cylinder with a diameter of d and a height of 0.4 mm. Referring to Formula (5), the ABH structure parameters are shown in Table 2. The structure diagram is shown in Figure 3.     The specific material properties and geometric dimensions are shown in Table 1. The thickness of the ABH changes according to the power function law. Since it is difficult to reduce the thickness to 0 due to the processing conditions, the truncated thickness is taken as 0.4 mm, so the central area of the ABH is a cylinder with a diameter of d and a height of 0.4 mm. Referring to Formula (5), the ABH structure parameters are shown in Table 2. The structure diagram is shown in Figure 3.   As shown in Figure 4, the PZT is pasted on the back of the central part of the ABH structure. PZT-5H can be selected as the PZT material, and the shape is round and leads out the circuit to connect the load on the positive and negative electrodes of the PZT. This is the vibration energy harvester based on the two-dimensional acoustic black hole designed in this paper. Figure 4, the PZT is pasted on the back of the central part of the A structure. PZT-5H can be selected as the PZT material, and the shape is round and le out the circuit to connect the load on the positive and negative electrodes of the PZT. T is the vibration energy harvester based on the two-dimensional acoustic black hole signed in this paper.

The Energy Concentration Effect of 2-D ABH
Based on the COMSOL modal analysis for flat plate and plate with ABH struct (ABH plate) to obtain the characteristic frequencies. Some frequencies are shown in Ta 3. It can be seen from Table 3 that the mode of the ABH plate has changed greatly a the characteristic frequency has decreased. It is calculated that the number of characte tic frequencies of the plate within 5000 Hz is 31, while the ABH plate is 36, which sho that the number of modes of the ABH structure has increased, as shown in Figure 5. It be seen that the energy concentration effect of the ABH structure is not obvious under mode of lower characteristic frequency, but with an increase in characteristic frequen the ABH structure shows a very strong energy accumulation effect.

The Energy Concentration Effect of 2-D ABH
Based on the COMSOL modal analysis for flat plate and plate with ABH structure (ABH plate) to obtain the characteristic frequencies. Some frequencies are shown in Table 3. It can be seen from Table 3 that the mode of the ABH plate has changed greatly and the characteristic frequency has decreased. It is calculated that the number of characteristic frequencies of the plate within 5000 Hz is 31, while the ABH plate is 36, which shows that the number of modes of the ABH structure has increased, as shown in Figure 5. It can be seen that the energy concentration effect of the ABH structure is not obvious under the mode of lower characteristic frequency, but with an increase in characteristic frequency, the ABH structure shows a very strong energy accumulation effect. The frequency domain analysis is carried out in the range of 10-5000 Hz, and the quadratic velocity of the upper surface of the ABH region is compared with that of the corresponding region of the flat plate, as shown in Figure 6. It can be seen that the characteristic frequency of an ABH plate is reduced, and additional characteristic frequencies have been generated near 2800 Hz and 4000 Hz. At the same time, it can be seen that the advantages of the quadratic velocity of the ABH structure are not obvious compared with that of the flat plate before 1000 Hz, but after that, the quadratic velocity is greater. In the 5000 Hz frequency band, the average value of the quadratic velocity of the ABH plate is 7.24 mm/s, and the flat plate is 3.77 mm/s, increasing by 1.9 times.  The frequency domain analysis is carried out in the range of 10-5000 Hz, and the quadratic velocity of the upper surface of the ABH region is compared with that of the corresponding region of the flat plate, as shown in Figure 6. It can be seen that the characteristic frequency of an ABH plate is reduced, and additional characteristic frequencies have been generated near 2800 Hz and 4000 Hz. At the same time, it can be seen that the advantages of the quadratic velocity of the ABH structure are not obvious compared with that of the flat plate before 1000 Hz, but after that, the quadratic velocity is greater. In the 5000 Hz frequency band, the average value of the quadratic velocity of the ABH plate is 7.24 mm/s, and the flat plate is 3.77 mm/s, increasing by 1.9 times. The frequency domain analysis is carried out in the range of 10-5000 Hz, and the quadratic velocity of the upper surface of the ABH region is compared with that of the corresponding region of the flat plate, as shown in Figure 6. It can be seen that the characteristic frequency of an ABH plate is reduced, and additional characteristic frequencies have been generated near 2800 Hz and 4000 Hz. At the same time, it can be seen that the advantages of the quadratic velocity of the ABH structure are not obvious compared with that of the flat plate before 1000 Hz, but after that, the quadratic velocity is greater. In the 5000 Hz frequency band, the average value of the quadratic velocity of the ABH plate is 7.24 mm/s, and the flat plate is 3.77 mm/s, increasing by 1.9 times.

Influence of the ABH Power Index on Vibration
When the ABH length L and the truncation thickness h 1 are fixed, the power index m affects the ABH section depth, as shown in Figure 7. When the ABH length L and the truncation thickness h1 are fixed, the power index m affects the ABH section depth, as shown in Figure 7. When m = 0, it is a flat plate. With the gradual increase in m, the local area of the section will have a smaller curvature radius. Here, m is taken as 2-20, and the step size is taken as 3. Table 4 shows the corresponding characteristic frequencies of high-order modes under different power exponents. It can be seen from Table 4 that the corresponding characteristic frequency of each mode shows a decreasing trend with the increase in power exponent m under the higher order mode, and this trend is gradual, that is, the characteristic frequency is closer to the increase in m value. Therefore, the vibration mode has a large deviation when the power index is low. The transient simulation of the plate with ABH is carried out with a power exponent m, under a 2.4 kHz excitation frequency. The change in displacement at the center point of the circular plane of ABH over time is shown in Figure 8. The results show that, with the increase in power exponent m, the propagation velocity of the wave lags and decreases, so the corresponding wave amplitude increases. In addition, when m is larger, the increase in the focus effect of ABH does not change much. When m = 0, it is a flat plate. With the gradual increase in m, the local area of the section will have a smaller curvature radius. Here, m is taken as 2-20, and the step size is taken as 3. Table 4 shows the corresponding characteristic frequencies of high-order modes under different power exponents. It can be seen from Table 4 that the corresponding characteristic frequency of each mode shows a decreasing trend with the increase in power exponent m under the higher order mode, and this trend is gradual, that is, the characteristic frequency is closer to the increase in m value. Therefore, the vibration mode has a large deviation when the power index is low. The transient simulation of the plate with ABH is carried out with a power exponent m, under a 2.4 kHz excitation frequency. The change in displacement at the center point of the circular plane of ABH over time is shown in Figure 8. The results show that, with the increase in power exponent m, the propagation velocity of the wave lags and decreases, so the corresponding wave amplitude increases. In addition, when m is larger, the increase in the focus effect of ABH does not change much.

Influence of ABH Truncation Thickness on Vibration
The frequency response of the plate with ABH is numerically simulated with different truncation thicknesses h1, under an excitation frequency ranging from 10 Hz to 5000

Influence of ABH Truncation Thickness on Vibration
The frequency response of the plate with ABH is numerically simulated with different truncation thicknesses h 1 , under an excitation frequency ranging from 10 Hz to 5000 Hz. To evaluate the energy focus effect of ABH, the logarithm of the ratio of the quadratic velocities of the ABH area and the flat area is recorded as Ra, which calls the focus ratio: The focus ratio of the plate with ABH integrated under different truncation thicknesses h 1 was carried out through numerical simulation, which is shown in Figure 9. The influence of truncation thickness is not different in the lower frequency region. However, the difference in the focus effect of ABH becomes significant when the excitation frequency is higher than 1500 Hz. When the plate with 2-D ABH is at the resonance frequency of 3150 Hz, the focus effect is increased with the smaller truncation thickness. The maximum Ra can reach one at the minimum truncation thickness, h 1 = 0.4 mm.

Influence of ABH Truncation Thickness on Vibration
The frequency response of the plate with ABH is numerically simulated with different truncation thicknesses h1, under an excitation frequency ranging from 10 Hz to 5000 Hz. To evaluate the energy focus effect of ABH, the logarithm of the ratio of the quadratic velocities of the ABH area and the flat area is recorded as Ra, which calls the focus ratio: The focus ratio of the plate with ABH integrated under different truncation thicknesses h1 was carried out through numerical simulation, which is shown in Figure 9. The influence of truncation thickness is not different in the lower frequency region. However, the difference in the focus effect of ABH becomes significant when the excitation frequency is higher than 1500 Hz. When the plate with 2-D ABH is at the resonance frequency of 3150 Hz, the focus effect is increased with the smaller truncation thickness. The maximum Ra can reach one at the minimum truncation thickness, h1 = 0.4 mm.

Influence of ABH Cross-Sectional Length on Vibration
The frequency response of the plate with ABH is numerically simulated with different cross-sectional lengths L, under excitation frequencies ranging from 10 Hz to 5000 Hz. Figure 10 shows that the focus ratio Ra increases with the ABH cross-sectional length L, which means a better energy focus effect. However, for each cross-sectional length at a specific frequency, Ra will have a peak. For example, when the excitation frequency is 4000 Hz and L = 30 mm, the Ra can reach one. This may be due to the influence of the resonant frequency. The frequency response of the plate with ABH is numerically simulated with different cross-sectional lengths L, under excitation frequencies ranging from 10 Hz to 5000 Hz. Figure 10 shows that the focus ratio Ra increases with the ABH cross-sectional length L, which means a better energy focus effect. However, for each cross-sectional length at a specific frequency, Ra will have a peak. For example, when the excitation frequency is 4000 Hz and L = 30 mm, the Ra can reach one. This may be due to the influence of the resonant frequency.

Optimization and Analysis of ABH Size
To achieve the best energy focus effect of ABH, four structural parameters are opti-

Optimization and Analysis of ABH Size
To achieve the best energy focus effect of ABH, four structural parameters are optimized in this paper, including the section length L of ABH, the power exponent m, the truncation thickness h 1 , and the diameter of the central circular platform d. Orthogonal test table L25 (56) is used for studying the effects of the influencing factors, as shown in Table 5. According to the orthogonal test table, 25 experiments were completed. The kinetic energy density of the upper and lower round surfaces (diameter 15 mm) of the ABH structure was integrated. Twenty-five groups of experimental results, according to the integral results of Experiment 1, were normalized. The experimental results are shown in Table 6. It can be seen from the table that the order of the range values of each factor is: truncation thickness > cross-sectional length > round table diameter > power index. The order of the error square sum of each factor is the same. Therefore, the order of influence of ABH structure factors on the simulation results is the same. Table 6. Visual analysis of ABH structural parameters.

Factor
Power Index Truncation Thickness

Cross-Sectional Length
Round In Table 6, the energy focus effect of Experiment 10 is the best, so the parameters of Experiment 10 are taken as the basis for the optimization of the structure parameters. Based on the parameters of Experiment 10, the parametric sweep is carried out on the section's length, central cone diameter, and power exponent, respectively. The results were normalized according to Experiment 1, as shown in Figure 11. To sum up, when the power index, truncation thickness, cross-sectional length, and round table diameter of the ABH are 3, 40 mm, 0.4 mm, and 24 mm, respectively, the ABH has the best energy focus effect. This group of structural parameters is taken as the parameters of Experiment 26.

Frequency Domain Analysis of Different Experiments
The simulation results of Experiment 10, Experiment 3, and the flat plate are compared to verify the above conclusions in the frequency domain analysis. Apply a simple harmonic excitation with an amplitude of 10 N at the right end of the cantilever plate, and the frequency range is 10-5000 Hz. Divide the frequency band into low frequency (10-1400 Hz) region, medium frequency (1400-3600 Hz) region, and high frequency (3600-5000 Hz) region. Compare the effects of the three groups of experiments with the quadratic velocities of the upper and lower surfaces of the ABH central platform (the circular area of the flat plate at the same position). The quadratic velocity at different frequencies is shown in Figure 12. The average value of the quadratic velocity at each frequency region is shown in Figure 13. To sum up, when the power index, truncation thickness, cross-sectional length, and round table diameter of the ABH are 3, 40 mm, 0.4 mm, and 24 mm, respectively, the ABH has the best energy focus effect. This group of structural parameters is taken as the parameters of Experiment 26.

Frequency Domain Analysis of Different Experiments
The simulation results of Experiment 10, Experiment 3, and the flat plate are compared to verify the above conclusions in the frequency domain analysis. Apply a simple harmonic excitation with an amplitude of 10 N at the right end of the cantilever plate, and the frequency range is 10-5000 Hz. Divide the frequency band into low frequency (10-1400 Hz) region, medium frequency (1400-3600 Hz) region, and high frequency (3600-5000 Hz) region. Compare the effects of the three groups of experiments with the quadratic velocities of the upper and lower surfaces of the ABH central platform (the circular area of the flat plate at the same position). The quadratic velocity at different frequencies is shown in Figure 12. The average value of the quadratic velocity at each frequency region is shown in Figure 13.  As can be seen from Figure 12, when the frequency is below 600 Hz, the experimental results of the three groups of parameters are close, and the energy focus effect of the ABH has almost disappeared. When the excitation frequency is higher than 600 Hz, the quadratic velocity of the plate with ABH is higher than that of the flat plate, which means the energy focus effect of the ABH has appeared. The best performance of the energy focus effect is Experiment 26. The quadratic velocity of the plate decreases rapidly when the excitation frequency is far away from the characteristic frequency of the plate, however, the reduction in speed of the ABH plate is smaller than that of the flat plate. It shows that the energy focus effect of ABH can be affected over a wider excitation frequency range.
The same conclusion can be drawn from Figure 13. In the low-frequency region, the average quadratic velocity of Experiment 26 is 1.79 times and 1.33 times higher than that of the flat plate and Experiment 3, respectively. In the medium-frequency region, the multiples are 7.6 and 3.05, respectively. In the medium-frequency region, the multiples are 7.6 and 3.05, respectively. In the high-frequency region, the multiples are 22.33 and 2.02, respectively. It can be seen that the ABH structure has obvious advantages over the flat plate, especially in the high-frequency region. In the frequency region of 10-5000 Hz, the  As can be seen from Figure 12, when the frequency is below 600 Hz, the experimental results of the three groups of parameters are close, and the energy focus effect of the ABH has almost disappeared. When the excitation frequency is higher than 600 Hz, the quadratic velocity of the plate with ABH is higher than that of the flat plate, which means the energy focus effect of the ABH has appeared. The best performance of the energy focus effect is Experiment 26. The quadratic velocity of the plate decreases rapidly when the excitation frequency is far away from the characteristic frequency of the plate, however, the reduction in speed of the ABH plate is smaller than that of the flat plate. It shows that the energy focus effect of ABH can be affected over a wider excitation frequency range.
The same conclusion can be drawn from Figure 13. In the low-frequency region, the average quadratic velocity of Experiment 26 is 1.79 times and 1.33 times higher than that of the flat plate and Experiment 3, respectively. In the medium-frequency region, the multiples are 7.6 and 3.05, respectively. In the medium-frequency region, the multiples are 7.6 and 3.05, respectively. In the high-frequency region, the multiples are 22.33 and 2.02, respectively. It can be seen that the ABH structure has obvious advantages over the flat plate, especially in the high-frequency region. In the frequency region of 10-5000 Hz, the As can be seen from Figure 12, when the frequency is below 600 Hz, the experimental results of the three groups of parameters are close, and the energy focus effect of the ABH has almost disappeared. When the excitation frequency is higher than 600 Hz, the quadratic velocity of the plate with ABH is higher than that of the flat plate, which means the energy focus effect of the ABH has appeared. The best performance of the energy focus effect is Experiment 26. The quadratic velocity of the plate decreases rapidly when the excitation frequency is far away from the characteristic frequency of the plate, however, the reduction in speed of the ABH plate is smaller than that of the flat plate. It shows that the energy focus effect of ABH can be affected over a wider excitation frequency range.
The same conclusion can be drawn from Figure 13. In the low-frequency region, the average quadratic velocity of Experiment 26 is 1.79 times and 1.33 times higher than that of the flat plate and Experiment 3, respectively. In the medium-frequency region, the multiples are 7.6 and 3.05, respectively. In the medium-frequency region, the multiples are 7.6 and 3.05, respectively. In the high-frequency region, the multiples are 22.33 and 2.02, respectively. It can be seen that the ABH structure has obvious advantages over the flat plate, especially in the high-frequency region. In the frequency region of 10-5000 Hz, the average value of the quadratic velocity of Experiment 26 is increased by 5.03 times and 1.81 times compared with the flat plate and Experiment 3, respectively, showing the apparent advantages of the energy focus effect after the optimization of the ABH structure parameters. It shows that the energy focus effect is more significant with optimized ABH structure parameters.

Output Performance Simulation of VEH
Based on the above research, to further verify the advantages of the VEH composed of ABH. The frequency domain output performance of VEH is simulated. The circular PZT (PZT-5H) is pasted on the back of the central part of the ABH structure (corresponding to the position of the flat plate). An external resistance (1000 Ω) of the circuit on the positive and negative electrodes of the PZT. This is the VEH based on the two-dimensional acoustic black hole designed in this paper. The output performance of Experiment 26 VEH and the flat plate VEH is analyzed in the frequency domain by using simple harmonic excitation with an amplitude of 5 N. The simulation results are shown in Figure 14.
Micromachines 2023, 14, x FOR PEER REVIEW 14 of 15 average value of the quadratic velocity of Experiment 26 is increased by 5.03 times and 1.81 times compared with the flat plate and Experiment 3, respectively, showing the apparent advantages of the energy focus effect after the optimization of the ABH structure parameters. It shows that the energy focus effect is more significant with optimized ABH structure parameters.

Output Performance Simulation of VEH
Based on the above research, to further verify the advantages of the VEH composed of ABH. The frequency domain output performance of VEH is simulated. The circular PZT (PZT-5H) is pasted on the back of the central part of the ABH structure (corresponding to the position of the flat plate). An external resistance (1000 Ω) of the circuit on the positive and negative electrodes of the PZT. This is the VEH based on the two-dimensional acoustic black hole designed in this paper. The output performance of Experiment 26 VEH and the flat plate VEH is analyzed in the frequency domain by using simple harmonic excitation with an amplitude of 5 N. The simulation results are shown in Figure 14. From Figure 14, in the low-frequency region, the average power of VEH in Experiment 26 is two times that of the flat plate. In the medium-frequency region, the multiple is 1.26. In the high-frequency region, the multiple is 4.25. It can be seen that the ABH structure has obvious advantages over the flat plate, especially in the high-frequency region. This is consistent with the result of the reaction in Figure 13, which further shows that the optimized ABH structural parameter VEH has a more excellent output performance.

Conclusions
The vibration energy harvester based on a two-dimensional ABH is studied in this article. Firstly, based on the working principle of ABH, the ABH effect and the influence of ABH structural parameters on the vibration are studied by simulation. Secondly, based on the orthogonal test method, the optimal ABH structure parameters are determined in the time domain analysis, namely, when the power index is 3, the truncated thickness is 0.4 mm, the cross-sectional length is 40 mm, and the round table diameter is 24 mm. Thirdly, the optimal ABH structure parameters are verified by frequency domain analysis. The results show that ABH structure can play a significant role in broadband, and the energy focalization effect is more significant after optimizing the ABH structure parameters. Finally, the output voltage contrast experiment shows that the output performance of VEH in Experiment 26 is significantly improved compared with that of the flat plate. From Figure 14, in the low-frequency region, the average power of VEH in Experiment 26 is two times that of the flat plate. In the medium-frequency region, the multiple is 1.26. In the high-frequency region, the multiple is 4.25. It can be seen that the ABH structure has obvious advantages over the flat plate, especially in the high-frequency region. This is consistent with the result of the reaction in Figure 13, which further shows that the optimized ABH structural parameter VEH has a more excellent output performance.

Conclusions
The vibration energy harvester based on a two-dimensional ABH is studied in this article. Firstly, based on the working principle of ABH, the ABH effect and the influence of ABH structural parameters on the vibration are studied by simulation. Secondly, based on the orthogonal test method, the optimal ABH structure parameters are determined in the time domain analysis, namely, when the power index is 3, the truncated thickness is 0.4 mm, the cross-sectional length is 40 mm, and the round table diameter is 24 mm. Thirdly, the optimal ABH structure parameters are verified by frequency domain analysis. The results show that ABH structure can play a significant role in broadband, and the energy focalization effect is more significant after optimizing the ABH structure parameters. Finally, the output voltage contrast experiment shows that the output performance of VEH in Experiment 26 is significantly improved compared with that of the flat plate.