Abstract

Acoustic metamaterial, which can prohibit effectively the elastic wave propagation in the bandgap frequency range, has broad application prospects in the vibration and noise reduction areas. The Lindstedt–Poincaré method was utilized to analyze the dispersion curves of nonlinear metamaterial with a pure Duffing oscillator. The first-order perturbation solutions of acoustic and optical branches were obtained. Both the starting and cutoff frequencies of the bandgap are determined consequently. It was found that the soft/hard characteristics of pure Duffing oscillators could lead to the lower/upper movement of the starting and cutoff frequencies of the bandgap. By further researching the degraded linear system, the conclusion that actual nonlinear metamaterial bandgap region is wider than effective negative mass region is drawn and that both mass and stiffness ratio effect on the starting frequency is obtained. Effective positive mass can also lead to the vibration attenuation in bandgap. For nonlinear metamaterial, the translation effect of the external excitation amplitude on the bandgap range and the zero mass at the nonlinear bandgap cutoff frequency were discussed, and all above conclusion are identified by numerical analysis.

1. Introduction

In the last ten years, a significant attention is directed toward so called metamaterial, a kind of artificial structures, which is composed of small substructure that behaves like a continuous material. The most remarkable property of metamaterials is the bandgap in which the acoustic/elastic waves propagate is prohibited [1]. With this unique characteristic, metamaterials can be applied to the wave filtering [24], vibration attenuation [5], acoustic isolation [6, 7], and sonic transmission [8, 9]. In civil engineering, acoustic metamaterials with the band characteristics were also created to replace the traditional seismic designs for isolating seismic waves from buildings. Shi and his coworkers [1014] proposed a novel concept of foundation, periodic foundation, to isolate civil structures from seismic wave with frequencies that might resonate with the structures, which has been investigated in a number of comprehensive theoretical, numerical, and experimental studies. In particular, Brûlé et al. [15] carried out the large-scale test of seismic metamaterials [15] such as seismic metawedge capable [16] capable of creating seismic waves and studied the attenuation effect of seismic waves by using finite periodic lattices of mass-in-mass barriers [17].

Effective negative mass metamaterials come from mass effective ratio, an artificial parameter, which is calculated analogically as the parameter of electromagnetic metamaterials [18, 19] and has not only positive but also negative values. Based on the idea of locally resonant microstructure, Liu et al. [20] designed the first effective negative mass metamaterial in certain frequency region. Subsequently, Liu et al. [21], Milton and Willis [22], Willis [23], and Huang et al. [24, 25] proposed various analytic models to analyze the related effective negative mass mechanisms. Yao et al. [26] realized effective negative mass and zero-mass phenomena in a 1D mass-spring system experimentally. Other different types of acoustic metamaterials with negative effective mass have been proposed [2732]. From the existing literatures, the starting frequency bandgap of metamaterials is defined by local resonance frequency where the effective negative mass just began. However, the bandgap starting frequency will change with outer spring stiffness when the local resonance frequency remains unchanged [26], and actually the starting frequency is lower than resonance frequency in experimental research.

Nonlinear acoustic metamaterials also deserve special attention. Vakakis and his coauthors [33, 34] utilized the multiple scales perturbation method to analyze nonlinear chains subjected to external forcing and ground springs and developed nonlinear dispersion relationships that exhibit amplitude dependence. Chakraborty and Mallik [35] studied the cubic chain and the effect of nonlinearities on the propagation constant and natural frequencies. Then Lazarov and Jensen [36] considered a linear chain with attached nonlinear damped oscillators by balance approach. Marathe and Chatterjee [37] looked at a damped nonlinearity and used harmonic balance and multiple scales to uncover the decay rate in the propagation zone. Narisetti et al. [38] developed a Lindstedt–Poincaré perturbation technique and analyzed monoatomic cubic chains to capture dispersion and bandgap shifts. However, their models are not based on effective negative mass, and Duffing oscillator is often a damped one. The damp may have great influence on the starting and cutoff frequencies of bandgap [39]. As a result, in order to investigate how effect of the nonlinearity on the effective negative mass nonlinear metamaterial dispersion curves, the pure Duffing oscillator, which is no damped pure Duffing oscillator is adaptive

It is worth mentioning that, there is no literature discussion whether positive effective mass phenomenon exists in the nonlinear metamaterial bandgap based on effective negative mass.

In this paper, effective negative mass nonlinear metamaterial mathematical model is given firstly, and Lindstedt–Poincaré perturbation method is utilized to calculate metamaterial dispersion curves. Then, the precise expression of both starting and cutoff frequencies is deduced consequently. After that, the result is obtained that positive effective mass phenomenon exists in the nonlinear bandgap. Finally, the numerical calculation verification results agree with our theoretical analysis.

2. The First-Order Solution for the Dispersive Curve of Acoustic Metamaterial of Effective Negative Mass with Pure Cubic Oscillator

Consider an acoustic metamaterial containing a pure Duffing oscillator, illustrated as Figure 1. The shell mass is M, the rigidity coefficient of the spring connecting the two shells is K, and the mass of the internal oscillator is m. The internal oscillator is a pure Duffing oscillator whose force is proportional to displacement and the cubic of displacement, and the linear rigidity coefficient of the spring connecting the oscillator and the shell is k; the coefficient of the cubic term is Γ.

In the case of period j, for instance, and are the displacements of the shell and the oscillator in period j; and are the displacements of the shell and the oscillator in period j−1; and and are the displacements of the shell and the oscillator in period j+ 1. Next, we are going to discuss kinematic modeling of this system.

From the load on the shell in period , its differential equation of motion iswhere overdots denote derivative with respect to time t.

According to the load on the oscillator in period , its differential equation of motion is

Rewrite Equations (1) and (2) into a matrix

Define the external/internal spring rigidity ratio as , the shell/oscillator mass ratio as , and the linear natural frequency of oscillator as by introducing dimensionless time and dimensionless nonlinear coefficient . This way, we can nondimensionalize Equation (3):here is a dimensionless frequency. Now we are going to find the first-order perturbation solution of by Lindstedt–Poincaré method.

Using the asymptotic expansion below:where the superscript and subscript 0 and 1 are the linear and first-order asymptotic expansion alternatively.

Substitute Equation (5) into Equation (4):

Separate the and coefficients and order them to be zero. Then, we get,

coefficient:

coefficient:

2.1. Linear Dispersive Curve Equation

The steady-state solution of Equation (7) is expressed aswhere is the periodic unit size, is the periodic phase factor, and are the steady-state amplitude of shell and oscillator, respectively. stands for the conjugation of the equations above.

Substitute Equations (9) and (10) into Equation (7):

In order for Equation (11) to have nonzero solutions , the determinant of its coefficient matrix must be zero. Hence,and and have the following relation:

By finding in Equation (12), we get the dispersive curve equation of the two frequency branches: where and stand for the acoustic and optical frequency branches. An acoustic frequency branch is a low-frequency dispersive curve while an optical frequency branch is a high-frequency dispersive curve. On this basis, we can obtain linear dispersive curve Equations (14) and (15).

2.2. Nonlinear Dispersive Curve Equation

Next, let us consider the first-order perturbation solution which, for Equation (8), is expressed as

Substitute Equations (16) and (17) into Equation (8):where .

To remove the secular term, considering only the equality between the coefficients on both sides of Equation (18). Then we have

Order , and consider Equation (12). Then, we have

Through linear transform of the augmented matrix of Equation (19), we get

In order for Equation (19) to have nonzero solutions , the rank of its coefficient matrix must be the same as that of its augmented rank. Hence,

We get

So, the perturbed solution is expressed as

Considering the relation between the and shown in Equation (13), Equation (23) can be written as

By substituting and shown in Equations (14) and (15) into the equation above, we can obtain the dispersive curve equations of the first-order asymptotically expanded acoustic and optical frequency branches:

These are the dispersive curve equations of the acoustic and optical frequency branches that include for nonlinear factors. We can see that these dispersive curve equations are relevant to nonlinear small perturbation , coefficient before nonlinear term , and steady-state vibration amplitude of the pure Duffing oscillator .

At , the influence of coefficient before the nonlinear term , on the dispersive curve is analyzed as presented in Figure 2. Here, represents a hard Duffing oscillator, represents a soft Duffing oscillator, and represents linearity. From this diagram, when the local resonance oscillator is a hard Duffing oscillator, its two dispersive curves will translate toward the high band relative to the linear oscillator; when the local resonance oscillator is a soft Duffing oscillator, on the contrary, its two dispersive curves will translate toward the low band relative to the linear oscillator. Also, under the same conditions, nonlinear factors make greater differences to the dispersive equation of the optical frequency branch than to that of the acoustic frequency branch.

2.3. Rediscussion on the Starting and Cutoff Frequencies of a Corresponding Linear System

For acoustic frequency branch in Equation (14), when , the bandgap starting frequency is

It is easy to prove that, in Equation (27), . Hence, the bandgap dimensionless starting frequency is smaller than 1. Namely, the bandgap starting frequency is smaller than the resonance frequency . Furthermore, is the start area of effective negative mass. In the literature [26], is used to define the bandgap starting frequency. Comparing with the value, we can see that it is higher than the actual starting frequency. It is relatively conservative for bandgap design, but does not reflect the influence of spring stiffness of the shell.

Figure 3 shows how the starting frequency changes when and are defined differently. From Figure 3, the variation of the nondimensionalized starting frequency is not only relevant to the shell/oscillator rigidity ratio , but also to their mass ratio . Hence it is inappropriate to simply define the bandgap starting frequency as a resonance frequency. Literature [26] offers a good example supporting our observation with resonance frequency of 6.35 Hz and the starting frequency at 5.8 Hz.

For in Equation (15), when , the bandgap cutoff frequency is

From Equation (28), the bandgap cutoff frequency is only relevant to the shell/oscillator mass ratio and decreases when the mass ratio increases. And the cutoff frequency is corresponding to the effective zero mass.

Considering the dispersive curve equations of the first-order asymptotically expanded acoustic and optical frequency branches shown in Equations (25) and (26), for nonlinear acoustic metamaterial with weak nonlinearity (ε is small), nonlinear system characteristics depend on the corresponding degraded linear systems. As a result, in linear system the phenomenon that both effective and negative mass in bandgap and the phenomenon that cutoff frequency is effective zero mass should also exist in the nonlinear system However, none of the above has been reported. Then, verify them by numerical methods as mentioned in the following.

3. Numerical Simulation and Validation

The bandgap characteristics of infinite acoustic metamaterial can be assumed by the finite one (the periodic number is less than 5) [40]. In order to validate our theoretical analysis above, we are going to build a model of an effective negative mass system with pure Duffing oscillator in five periods as shown in Figure 4. By adding a simple harmonic exciting force on the shell of period 1 and picking up the response in period 5, the steady-state displacements of the shell and the oscillator are calculated at each of the excitation frequencies. This way, we get the vibration transmission characteristics of this finite periodic structure to analyze the bandgap and vibration transmission pattern of this structure. When the coefficient before the nonlinear term is , it represents the corresponding linear system.

3.1. Influence of Nonlinear Terms on the Bandgap Starting and Cutoff Frequencies

To validate the theoretical analysis result above, we may compare it effectively against experimental parameters provided in the literature [26] and establish a periodic structure vibration system of five periods with shell mass  = 0.1011 kg, internal oscillator mass  = 0.04647 kg, internal nonlinear spring linear rigidity  = 37 N/m, external linear spring damping coefficient c = 0.05 NS/m, and rigidity coefficient of spring linking the external oscillators K = 117 N/m. If the nonlinear perturbation term is taken as ε = 0.01, to facilitate expression, we define

When is positive, Duffing is hard rigid; when n is negative, Duffing is soft rigid; when is zero, the system has been degraded into a linear system. The vibration transmission characteristic curves of the corresponding system are calculated at  = 0, 1, −1, 3, −3, 10, and −10 as presented in Figure 5. The vibration transmission can be confirmed by the ratio of the fifth shell steady-state amplitude to the first shell steady-state amplitude.

As illustrated in Figure 5, existence of a nonlinear term does not make much difference to the bandgap, mainly because the Duffing oscillator we selected is weakly nonlinear. By enlarging the points of the bandgap starting and cutoff frequencies in the chart, we will see how a nonlinearity factor affects the bandgap.

When the longitudinal coordinate (i.e., transmission rate) of the transmission characteristic curve turns from positive to negative, the vibration in period 1 cannot be transmitted to the shell of period 5. Then, a bandgap will appear. Hence, we have to find the point where the transmission curve runs across the horizontal axis, i.e., the bandgap starting frequency, as presented in Figure 6.

From Figure 6(a), when is positive, the corresponding internal Duffing oscillator is a hardened Duffing oscillator. Hence, we can observe that, as increases, the bandgap starting frequency gradually translates toward the high-frequency region. When is negative, the corresponding internal Duffing oscillator is a softened Duffing oscillator. Hence, as decreases, the bandgap starting frequency gradually translates toward the low-frequency region.

The same method also applies to the influence of nonlinear factor on the bandgap cutoff frequency. As shown in Figure 6(b), when is positive, the corresponding internal Duffing oscillator is a hardened Duffing oscillator. Hence, we can observe that, as increases, the bandgap cutoff frequency gradually translates toward the high-frequency region. When is negative, the corresponding internal Duffing oscillator is a softened Duffing oscillator. As decreases, the bandgap cutoff frequency gradually translates toward the low-frequency region.

As discussed above, existence of nonlinear factor affects the bandgap starting and cutoff frequencies of the structure. The smaller the value is, the smaller the corresponding starting or cutoff frequency is; the larger the value is, the larger the starting or cutoff frequency is. Under the same conditions, nonlinear factor makes greater difference to cutoff frequency than it does to starting frequency.

3.2. Influence of the Shell Rigidity on the Bandgap Frequency

As analyzed above, existence of a nonlinear term makes a little difference to the bandgap starting frequency of the corresponding linear system; starting frequency is relevant to shell/oscillator rigidity ratio and mass ratio , and cutoff frequency is only relevant to mass ratio . As decreases, the starting frequency decreases, the cutoff frequency remains the same, and the bandgap width increases. Now, let us draw the vibration transmission profile of the corresponding system by only changing the rigidity coefficient of the spring between the oscillator and the shell, K, from K = 117 N/m to 37 N/m and taking the n value as n = 0, 1, −1, 3, −3, 10, and −10 as presented in Figure 7.

When the shell connecting rigidity K = 117 N/m is changed to 37 N/m, the masses of the Duffing oscillator and the shell remain the same, and according to analysis results in Equation (27) and Figure 3, after the shell rigidity is changed, the bandgap cutoff frequency will remain the same, but the starting frequency will decrease. That is, the bandwidth of the bandgap will increase, and the bandgap will expand toward the low-frequency region. As can be observed from Figure 5, the bandgap cutoff frequency remains the same at 7.76 Hz while the bandgap starting frequency has decreased to 3.76 Hz, proving our analysis result. Hence, if we want a wider bandgap width, when other conditions remain unchanged, the rigidity coefficient of the linear spring connecting the shell must be as small as practically possible.

3.3. Responses of the Shell and the Oscillator on the Bandgap

As can be summarized from the discussions on bandgap starting frequency above, existence of weakly nonlinear factor results in some minor changes near the frequency value of its degraded linear system, and that both the Duffing oscillator/shell linear rigidity ratio and mass ratio make a difference to the starting frequency and that the starting frequency is lower than the resonance frequency of a linear oscillator In other words, generation of bandgap is not always resulted from effective negative mass. Within the spectrum of bandgap, it must be effective positive mass below the linear resonance frequency constituted by the Duffing oscillator and the shell in each period, and effective negative mass does not appear except above the resonance frequency and below the cutoff frequency. Next, the responses of the shells and the oscillators of the system in the effective negative mass area and the effective positive mass area within the spectrum of bandgap are examined with rigidity ratio between springs connecting external large oscillators K = 37 N/m while keeping all the other parameters unchanged.

3.3.1. Effective Negative Mass Phenomenon in Bandgap

From Figure 7, with bandgap between 3.7 Hz and 7.67 Hz and resonance frequency of linear oscillator at 6.35 Hz, analysis is carried out under effective negative mass taking the external excitation frequency  = 7 Hz. Figure 8 shows the vibration time history curves and phase diagrams of the shells in five periods. Figure 9 shows the oscillators in five periods.

From Figures 8 and 9, some time later, the vibration amplitude of each shell and oscillator no longer attenuates but has stabilized; the displacement-velocity phase diagrams of the shell and the oscillator both consist of a circle, and all the shells and oscillators of this periodic structure have achieved steady-state vibration.

Now, let us compare the shell vibration in each period. From period 1 through period 5, the shell vibration gradually decreases. This vibration is almost zero in the last period. The internal Duffing oscillator shows very similar vibration characteristics to the shell. This indicates that vibration loaded on the shell in period 1 cannot be propagated in this periodic structure but is suppressed, which concurs with the characteristics of bandgap.

Figure 10 shows the vibration time histories of the shell and the Duffing oscillator in period 1. Here, the dotted lines are the time histories of the shell, and the solid lines are the time histories of the Duffing oscillator. As can be observed, when the motion displacement of the shell is positive amplitude, the motion of the Duffing oscillator, m, is a negative amplitude; when the motion displacement of is a negative amplitude, the motion displacement of m is also positive amplitude. This means that and m always have an opposite motion displacement, and the vibrations of the shell and the oscillator are reverse phased. Here, the black curves are the time displacement curves of m vibration of the internal mass bock. This is because has a larger mass than m, and under the same external excitation amplitude, objects with a smaller mass will find it easier to change their velocity, resulting in a larger displacement. Further analysis revealed that the Duffing oscillator has a greater momentum than the shell. Now, the system comprising the shell and the Duffing oscillator in individual periods has a negative momentum. As the effective velocity is positive, the system possesses negative mass, and consequently the system helps control vibration.

3.3.2. Effective Positive Mass Phenomenon in Bandgap

Based on the above analysis, bandgap in 6.35 Hz–7.67 Hz, as shown in Figure 7, where the effective negative mass phenomenon occurs, meanwhile bandgap in 3.7 Hz–6.35 Hz should be the effective positive mass region. Analysis is carried out under effective positive mass taking the external excitation frequency  = 5 Hz. Figure 11 shows the vibration time history curves and phase diagrams of the shells in five periods, and Figure 12 shows the oscillators in five periods.

From Figures 11 and 12, all the shells and oscillators of this periodic structure have achieved steady-state vibration. The vibrations loaded on the shell and oscillators in period 1 cannot be propagated in this periodic structure which concurs with the characteristics of bandgap.

Figure 13 shows the vibration time histories of the shell and the Duffing oscillator in period 1. As can be observed, when the motion displacement of the shell is positive amplitude, the motion of the Duffing oscillator, , is also of the same positive amplitude. and always have the same motion displacement, and the vibrations of the shell and the oscillator in the same period are cophasal. Now, the system in every period must have a positive momentum, and the system helps control vibration just as effective as negative mass systems does. That is, effective positive mass phenomenon also occurs in the bandgap of nonlinear metamaterial.

3.3.3. Discussion on Effective Zero Mass

As analyzed above, when the effective mass gradually changes from negative to positive through a frequency point at which the effective mass is zero, the place corresponding to zero mass is where the bandgap cutoff frequency stands. Existence of nonlinear factor, however, makes a minor difference to the frequency point at zero mass, and it would be possible that the actual cutoff frequency of bandgap cannot be found. To solve this problem, we can closely sample the external excitation frequency so as to bind the bandgap cutoff frequency of the nonlinear system, i.e., the point of zero-mass frequency. Finally, with a constant nonlinear factor n = 10, shell mass  = 0.1011 kg, internal oscillator mass  = 0.04647 kg, internal nonlinear spring linear rigidity  = 37 N/m, and external linear spring damping coefficient c = 0.1 NS/m, we found the frequency point at zero mass to be 8.2205 Hz. Figures 14 and 15 show the time histories of the shells and the Duffing oscillators using this frequency point. From these curves, the shells and Duffing oscillators in all the periods are moving. From the phase diagrams, both the shells and the Duffing oscillators have reached a steady state, and the shells in all the periods have the same maximum displacements as the Duffing oscillators in all the periods.

From these charts in Figures 14 and 15, the shells of all periodic units move in the same phase and with the same displacement pattern, keeping the spring K linking the periodic units from deformation. If we ignore the motion of the internal Duffing oscillators from outside, the entire system will appear to be a rigid rod. The entire system maintains an in-phase translation while all the internal Duffing oscillators move in the same phase but in the opposite direction to the shells, and there is no phase difference among the internal Duffing oscillators. This is because a zero mass unit signifies that its inertia force is zero, making it appear that the entire system is composed of massless rigid bodies. It can be proved that the springs connecting the units will not deform, and the entire system looks as if one massless rigid rod is translating. The lattice kymatology of the periodic system can also provide an explanation. This phenomenon similar to zero mass is how electromagnetic waves propagate through a matching metamaterial with zero refraction index in electromagnetics. This material also features zero dielectric constant and zero magnetic permeability and no phase difference in the electromagnetic field.

3.4. Influence of the Excition Amplitude on the Start and Stop Frequencies of the Bandgap

As is known to all, for a linear system, the system response is irrelevant to the external excitation amplitude. For a nonlinear system, on the contrary, the system response is usually sensitive to and can vary with the external excitation amplitude. Hence, we are going to analyze how bandgap changes with the external excitation amplitude.

With a constant nonlinear factor n, we change external excitation amplitude to examine how it affects the periodic system again with shell mass  = 0.1011 kg, internal oscillator mass  = 0.04647 kg, internal nonlinear spring linear rigidity  = 37 N/m, external linear spring damping coefficient c = 0.05 NS/m, and rigid coefficient between springs connecting external large oscillators  = 117 N/m. If the nonlinear perturbation term is taken as , , f is taken as 1, 2, 3, 5, and 10 to calculate the vibration transmission characteristics. The result is presented in Figure 16. When the external excitation amplitude changes, so do the bandgap start and stop frequencies and the bandgap bandwidth. Figure 17 shows enlarged areas of the starting and cutoff frequencies.

As can be observed from Figure 17(a), the external excitation amplitude makes quite a big difference to the bandgap starting frequency. As the external excitation amplitude f increases, the bandgap start frequency moves toward the high band; when f decreases, the bandgap starting frequency moves toward the low band.

Figure 17(b) shows an enlarged view near the stop frequency. From this, we can see that, as f increases, the bandgap stop frequency moves toward the high band; as f decreases, the bandgap stop frequency moves toward the low band. As the transmission characteristic chart is much more separated than the previous two subsections, changing the f value makes a greater difference to the starting and cutoff frequencies.

Through comparison above, we can see that external excitation amplitude makes a great difference to bandgap. As the external excitation amplitude increases, so does the bandgap starting frequency of the metamaterial. This will further reduce the bandgap width, which is negative for the application of vibration control. Hence, when a metamaterial of effective negative mass with Duffing oscillator is used for vibration control, the external excitation amplitude must be kept within a reasonable limit.

4. Conclusions

In this paper, the Lindstedt–Poincaré method is utilized to calculate the first-order perturbation solution of acoustic and optical branches for effective negative mass nonlinear acoustic metamaterial with pure Duffing oscillator. The starting and cutoff frequencies of bandgap are deduced. Then, the influence of mass ratio and stiffness ratio and nonlinearity influence on starting and cutoff frequency of bandgap are discussed. The main conclusions can be drawn as follows:(1)The nonlinear part of the dispersion curves solution has relationships with coefficient of nonlinear item, steady-state amplitude of Duffing oscillator, and . Both soft and hard Duffing oscillators have effect on the starting and cutoff frequencies. The soft/hard characteristic of Duffing oscillator leads to the trend that both starting and cutoff frequencies move down/up relative to the degraded linear system, and the stronger the nonlinearity is the more obvious the trend is.(2)Both starting and cutoff frequencies of this nonlinear metamaterial bandgap depend on and . The starting frequency is lower than the resonant frequency, and it is monotonically increasing with and decreasing with . The cutoff frequency is only decreasing with mass ratio.(3)The actual nonlinear metamaterial bandgap region is wider than effective negative mass region. The frequency range from the resonant frequency to cutoff frequency is corresponding to effective negative mass. The frequency range from starting frequency to the resonant frequency is corresponding to effective positive mass. The vibration attenuation of both effective negative and effective positive mass occurs in bandgap.(4)Effective zero-mass phenomenon also exists in this nonlinear metamaterial. This is corresponding to the cutoff frequency. On this frequency, there is no relative motion for every shell of this nonlinear metamaterial, and every Duffing oscillator also has the same characteristics. The phases of shell and Duffing oscillator in every period are opposite, and the aggregate momentum of every periodic unit is close to zero.(5)The external excitation amplitude has influence on bandgap. As the amplitude increases, the starting and cutoff frequencies shift to the high-frequency band.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

Gao Ming now works in College of Mechanical and Electronic Engineering, Shandong Agricultural University in Taian, China.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge the State Key Research Development Program of China (Grant no. 2016YFC0600708), the National Natural Science Foundation of China (Grant nos. 11172198 and 51605264), the National Basic Research Program of China (Grant nos. 2013CB035402 and 2014CB046800), the Shandong Province Higher Educational Science and Technology Program (Grant no. J15LH04), the Natural Science Foundation of Shandong Province (Grant nos. ZR2018 and MEE001), and the State Key Laboratory of Open Funds (Grant no. MDPC201601).