Influence of Hysteresis on the Vibration Control of a Smart Beam with a Piezoelectric Actuator by the Bouc–Wen Model

.e hysteresis property in a smart structure has attracted much attention from researchers for several decades. Hysteresis not only affects the response precision of the smart structure but also threatens the stability of the system. .is paper focuses on how the hysteresis property influences the control effect of vibration suppression for a smart beam. Furthermore, the Bouc–Wen model is adopted to describe the hysteresis property of a smart beam and the hysteresis parameters of the hysteresis model are identified with a genetic algorithm. Based on the identification results, the hysteresis model is validated to represent the hysteresis property of the smart beam. Based on the hysteresis model, model reference adaptive control is designed to explore the influence of hysteresis on the vibration control of the smart beam. With some simulations and experiments, it is found that the vibration control effect is influenced when the hysteresis item changes. .e vibration control effect will be improved when the hysteresis coefficient in the Bouc–Wen model, as the expected objective model of the adaptive reference model, is within a proper numerical range where the control system is stable. Furthermore, when the time delay is considered in the closed-loop control system, the principle of the hysteresis influence is different. .e results indicate that the hysteresis property affects not only the control effect but also the stability of the control system for a smart cantilever beam.


Introduction
Hysteresis is a common nonlinearity of piezoelectric materials [1][2][3]. Hysteresis is a lagging phenomenon where a physical effect on the smart material lags behind its cause. When a smart material acts as an actuator or a sensor of space structure, the response of the smart structure lags behind the input excitation of the system [4,5]. e nonlinearity not only affects the control precision of the smart structure with a control law but also threatens the stability of the control system [6]. erefore, it is an urgent challenge to address how the hysteresis property exactly influences the control precision and the stability of the control system. Some researchers have not considered the hysteresis property in their studies when performing active control investigations of smart structures [7]. Furthermore, some linear control laws are commonly designed based on a linear model for the smart structure. For example, PID control was implemented to damp both the free and forced vibrations of a thin-walled structure and compared with LQR optimal control [8], a new intelligent methodology was introduced to mitigate the vibration response of flexible cantilever plates [9], and a minimax-linear quadratic Gaussian controller was designed based on an uncertain system model for the positioning of a piezoelectric tube scanner used in an atomic force microscope [10]. ese linear control approaches are limited in terms of their applications to linear system models. Once the controlled system is requested to have high precision and a fast response, these linear control methods cannot meet the requirements. erefore, modeling the hysteresis nonlinearities of smart structures is important for designing a controller to explore the influence of hysteresis on the control of the smart structures.
e Preisach model and Prandtl-Ishlinskii model are rate-independent hysteresis models and are also phenomenological hysteresis models [16]. Based on the two types of models, the main control laws were adopted with only their inverse hysteresis model as open-loop control to compensate for the hysteresis property of the smart structure [17,18]. It is difficult to design closedloop control strategies on account of the two models. However, the Bouc-Wen model is a rate-dependent model and a physics-based hysteresis model. e hysteresis parameters of the Bouc-Wen model are easily identified [19][20][21]. In addition, feedback or feedforward control may be easily realized with its inverse model and a certain control law. e control laws applied to suppress vibration mainly include PID control [22], adaptive control [6,23], robust control [24], and intelligent control [25,26]. ese control laws have been used to compensate for the hysteresis property or track the expected trajectory. However, few studies have analyzed how to affect the control result by the hysteresis of a smart structure. Furthermore, compared with other control methods, the adaptive control has strongly robust, adaptive, and anti-interference properties. Because of these advantages, the model reference adaptive control [23,27] is adopted to explore the influence of hysteresis on the control effect of vibration suppression of smart structures with the Bouc-Wen model in this paper. e rest of this paper is organized as follows. A dynamic model for a smart flexible beam, dynamic model validations, and stability analyses are presented in Section 2. Model reference adaptive control with the Bouc-Wen model as an expected reference model is introduced in Section 3. Simulations and experiments on vibration suppression of the smart beam are discussed in Section 4. Finally, conclusions are presented in Section 5.

Dynamic Modeling of a Smart
Cantilever Beam

A Model of a Smart Beam with the Bouc-Wen Equation.
A flexible beam bonded with a piezoelectric actuator is shown in Figure 1. e smart beam is regarded as an Euler-Bernoulli beam, and its bending vibration is assumed to be ω(x, t). us, the dynamic partial differential equation of the flexible beam is given as where m is the linear mass, c is the damping coefficient, E is the elasticity modulus, I is the inertia moment, and Q is a line distributed force. e constitutive equation of the piezoelectric material is described as where s 1 and S 1 are the strain and stress, respectively, E 3 is the electric field, and D 3 is the electric displacement.
Moreover, s E 11 , d 31 , and ε X 33 are the elastic compliance, piezoelectric strain, and dielectric constants, respectively.
When there is no external stress and the electric field is E 3 � U/t p (U is the external applied voltage and t p is the thickness of the piezoelectric actuator), the strain of the piezoelectric actuator is derived as e stress is obtained as where E p is the elasticity modulus of the piezoelectric actuator. erefore, the moment M of force is expressed as where A p is the cross-sectional area of the piezoelectric actuator. Substituting equations (3) and (4) into equation (5), the moment is obtained to be as follows: where t b and h p are the thickness of the beam and the height of the piezoelectric actuator, respectively, and Furthermore, the input voltage U is expressed as where V(t) is a time function and H(x) is the Heaviside function. e moment of force and line distributed force are related as follows: Substituting equations (6) and (8) into equation (1) yields e bending vibration ω(x, t) is obtained by a superposition of vibration modes as where ϕ i (x) is a mode and q i (t) is the generalized coordinate: 2 Shock and Vibration where β 1 l b � 1.875, β 2 l b � 4.694, β 3 l b � 7.855,. . ., and l b is the length of the beam. Substituting equation (10) into equation (9) yields Because of the mode shape orthogonality, It is assumed that . e mode ϕ j (x) is multiplied by both sides of equation (12) ; then, equation (12) is integrated with respect to x as follows: e natural frequency is obtained as where m b and m p are the masses of the cantilever beam and the piezoelectric patch, respectively, I b and I p are the moments of inertia of the cantilever beam, and E p is the elastic modulus of the cantilever beam. en, assuming Both sides of equation (17) are multiplied by where ζ j is the j-th order damping coefficient.
Some corresponding physical parameters of the smart cantilever beam are given in Table 1. However, the smart structure has a nonlinear property, namely, the hysteresis phenomenon.
erefore, when the first order of the vibration modes of the beam with a piezoelectric actuator is considered with the Bouc-Wen equation, a dynamic model with hysteresis of the smart beam is derived from equation (16) as where z(t) denotes the hysteresis item and the parameters A z , n, α, β, and c are used to obtain the form of the hysteresis property.

Model Identification by a Genetic Algorithm.
For the above hysteresis parameters, A z , α, β, and c, a genetic identification algorithm [28,29] is adopted to identify the hysteresis phenomenon of the smart beam. e parameters A z , α, β, and c are estimated by the input voltage and output response (bending deflection ω(x 0 , t) at x 0 � 0.115 m) of the smart beam in the experiment. ese estimated parameters are substituted in the original model to identify the output response. erefore, an objective function is first defined as where f s is the standard deviation, X k is the k-th error between the identified data and experimental response, k is the sampling point, X is the average error, and N is the number of identified data points. Before the genetic algorithm is implemented, some initial conditions are set. e maximum number of generations is MAXGEN � 100, the number of individuals is NIND � 100, the crossover probability and mutation probability are 0.7 and 0.02, respectively, and the generation gap is GGAP � 0.8. e default ranges of the identified parameters in Table 2 are very important for the convergence rate of the iterations for the genetic algorithm. Moreover, the experimental response is collected with a laser displacement sensor when a sinusoidal voltage V � 150 × sin(2π × 2.5t) is applied to the piezoelectric actuator. erefore, with the input voltage and output response data from the experiment, the identified parameters A z , α, β, and c are obtained in the default ranges, as shown in Table 2. Figure 2 shows the convergence of the genetic algorithm when the hysteresis parameters A z , α, β, and c are identified. e minimum standard deviation converges to a stable state after a generation number of approximately 20. Moreover, when the identified parameters in the third column of Table 2 at the last generation number 100 are substituted into equation (18), the identified displacement is obtained and compared with the experimental displacement in Figure 3(a), and a part of Figure 3(a) from 26 s to 28 s is magnified in Figure 3(b). en, to display the detailed hysteresis property, a period of experimental and identified data is excerpted from Figure 3(a) and is shown Figure 4. e red identified discrete data are consistent with the blue experimental solid curve. As a result, the dynamic model equation (18) is valid for describing the hysteresis property of the smart cantilever beam with the genetic identification algorithm.

Stability Analysis for the Bouc-Wen Model.
For the constructed dynamic model, the stability of the identified system is analyzed. Based on the above identified parameters, equation (18) is first transferred as where en, if n � 1, a further equation simplification is made and equation (20) is expressed as where δ 1 � 1 when u 2 > 0, δ 1 � − 1 and δ 2 � 1 when u 3 > 0, and Next, equation (21) is given in vector form as where u � u 1 u 2 u 3 T , If P(u s ) � 0, u s is the equilibrium point in the phase plane. By the nonlinear function P(u), u s � u s1 u s2 0 T and u s2 � − (Ω 2 1 /A z )u s1 . erefore, the stability analyses of equation (22) are explored near the equilibrium point u s . In addition, equation (22) without external excitation voltage is written as where and O(·) is the minimum value of approximately ‖u − u s ‖ 2 . erefore, equation (23) is simplified as where Furthermore, the eigenvalues λ are calculated by the following equation: where I is the unit matrix. By equation (28), the characteristic equation is derived as If the system is stabilized, the real parts of its eigenvalues are negative. As a result, the third item plus the fourth item is greater than zero and is given as Finally, for the stabilized model of smart beam, the relation between the identified parameters A z , α, β, and c is presented as If u s2 � 0, then According to the identified parameter values in Table 2, assuming A z � 0.110a z , the smart beam system is stable at a z > − 4.5044 from equation (32).

Adaptive Control Law.
In this paper, the vibration suppression problem of a smart beam considering the hysteresis property is explored by using a model reference adaptive control (MRAC) [22,25]. For adaptive control, a reference model is assumed as an expected objective model. According to the error between the output states of the smart beam and the expected objective hysteresis model, the adaptive control law is designed online with gains K r and K u , as shown in Figure 5. First, the expected objective hysteresis model is given as where ζ eo is set as 0.707, a z is a hysteresis coefficient, and a z A z satisfies equation (32).
Substituting equation (34) into equation (18), the identified model is derived as en, the error between the outputs of the expected objective model and the smart beam is expressed as e � u eo − u. (36)

Furthermore, equation (36) is transformed into a statespace equation with equations (33) and (37), which is displayed as
where If the second item and the third item on the right of equation (37) are equal to zero, the error e of equation (36) is convergent and lim t⟶∞ e(t) � 0. e fourth item of equation (36) is only considered when the influence of hysteresis on the control effect is studied. erefore, the control gains K r (t ⟶ ∞) � K r∞ and K u (t ⟶ ∞) � K u∞ must satisfy the following equations: or

Moreover, substituting equation (40) into equation (37) yields
where e � e _ e T , K u � K u∞ − K u , and K r � K r∞ − K r . Furthermore, for the stability analyses of the control system, a Lyapunov function is constructed as where P, P u , and P r are positive definite matrices. In addition, trace (·) calculates the trace of a matrix. For equation (42), the derivative is obtained as Substituting equation (41) into equation (43) yields If the first item of equation (44) is negative, namely, A T eo P + PA eo � − Q, Q is a positive definite matrix. erefore, the second and third items are equal to zero. With K u � K u∞ − K u and K r � K r∞ − K r , the control gain matrices K u and K r are expressed as where R 1 � P u (K − 1 u∞ ) T and R 2 � P r (K − 1 u∞ ) T . If the actual control gains are realized in discrete form, equation (45) is also represented as where k is the sampling point and h is the sampling interval.   Shock and Vibration Figure 6 shows the control result with no hysteresis in the expected objective hysteresis model of equation (33). Figure 6(a) shows the control voltage. Moreover, Figure 6(b) displays the free vibration of the cantilever beam with adaptive control at a z � 0 compared with the result with no control. If the hysteresis coefficient a z changes, the control result for the free vibration of the smart beam is different. When the hysteresis coefficient a z is set to − 4, − 3,. . ., 8,. . ., it is found that the control effect at a z � 8 is the best. Figure 7 shows the control result considering the hysteresis property in the expected objective hysteresis model at a z � 8. Figure 7(a) shows the control voltage and Figure 7(b) presents the free vibration of the smart beam at a z � 8 compared with that at a z � 0. Obviously, the vibration amplitude at a z � 8 is less than that at a z � 0. Figure 8 displays the frequency response of the curves in Figure 7(b) obtained by a fast Fourier transform (FFT). e amplitude peak value of the frequency response at a z � 8 is 0.03083 mm, the amplitude peak value of the frequency response at a z � 0 is 0.05875 mm, and the amplitude peak value of the frequency response with no control is 0.2567 mm. e control effect considering the hysteresis property in the expected objective hysteresis model at a z � 8 is superior to that with a linear expected objective model at a z � 0.

Simulation and Experimental Verifications of Model Reference Adaptive Control
When the hysteresis coefficient a z changes from − 4 to 15, the amplitude peak values of the frequency response are given in Figure 9(a). If the practical time delay of the sensor and the processing time spent in the computer are not considered in the simulation, the amplitude peak value of the frequency response tends to be constant under a limited control voltage applied by the power equipment with hysteresis item a z gradually increasing. Moreover, when the hysteresis coefficient a z is less than − 4, the control system of the smart beam is not stable. However, when the time delay    Shock and Vibration 7 of the sensing data and the processing time of the PC are added in the closed-loop control in the simulation, the control system of the smart beam is not stable at a z < 5 and a z > 12. e amplitude peak value of the frequency response at approximately a z � − 2 is the minimum of the dotted curve in Figure 9(a). Figure 9(b) shows peak of amplitudefrequency characteristic curves with different time delays (0.02 s, 0.05 s, 0.056 s, 0.058 s, 0.059 s, and 0.06 s), and the limited peak value is the minimum at a z � − 5, − 6, and − 7, when the time delay are set to 0.06 s, 0.059 s, and 0.058 s, respectively. Moreover, different time delays will affect the stability of the control system. erefore, adding the time delay in the closed-loop control system will produce a different control result.

Experimental
Verifications. An experimental platform is set up to verify the vibration control effect of the smart beam by the adaptive model reference control with a Bouc-Wen model as the expected objective model. e experimental equipment is shown in Figure 10. e adaptive control law is realized in LabVIEW. e control voltage is transported to the data acquisition unit. A power supply is used to amplify the control voltage from the data acquisition unit. en, the control voltage is applied to the piezoelectric patch to suppress the vibration of the cantilever beam when the beam suffers from an initial disturbance. A laser displacement sensor is used to sense the bending deflection of the smart beam. Figure 11 shows the vibration control result at a z � 0 for the smart beam. Figure 11(a) shows the control voltage, while Figure 11(b) displays the free vibration of the beam when the initial displacement at the free end of the beam is approximately 2 cm. It is observed easily that the vibration amplitude with control at a z � 0 is 0.15 mm at 0.45 s. Furthermore, Figure 12 shows the control result at a z � − 6. Figure 12(a) shows the control voltage, and Figure 12(b) displays the free vibration of the cantilever beam with control at a z � − 6. Compared with the vibration amplitude at a z � 0, the vibration amplitude a z � − 6 is 0.09 mm at 0.4 s. Notably, the control effect at a z � − 6 is better than that at a z � 0.
Similarly, the frequency response of the three curves in Figure 12(b) is shown in Figure 13. e amplitude peak of the frequency response at a z � − 6 is 0.01157 mm at a frequency of 5.2490 Hz, the amplitude peak of the frequency response at a z � 0 is 0.02431 mm at a frequency of 8.8060 Hz, and the amplitude peak of the frequency response with no control is 0.2063 mm at a frequency of 2.5740 Hz. It is clearly indicated that the vibration control effect a z � − 6 for the cantilever beam is better than that at a z � 0.
When the hysteresis coefficient a z changes, the peak value in the experiment exhibits the same changing trend in contrast to the result of the simulation considering the time delay in closed-loop control system, as shown in Figure 14.
e solid curve is a polyfit of the curve from the experiment. e dashed curve is the amplitude peak value of the frequency response from the simulation considering the time delay. en, both the solid curve and the dashed curve have their own minimum peak values. However, the hysteresis coefficients a z of the two minimum peak values are different. One value is − 2 and the other value is − 6. e deviation between the experiment and simulation is caused by the time delay. In the experimental system, the time delay includes not only the sensing delay of the sensor but also the processing time of the PC. e sensing delay of the sensor is  Shock and Vibration obtained from the user manual of the sensor. However, the processing time of the PC in the simulation is not calculated accurately. When too many application programs are running on the computer, the processing velocity of the PC slows down and the processing time of PC is large, and vice versa. rough the simulation and experimental verifications, adaptive model reference control with a Bouc-Wen model as an expected objective model is proven to be feasible for suppressing the vibration of the cantilever beam when it suffers from an initial displacement at the free end. It is illustrated that the hysteresis property will affect the stability and control effect of the control system when considering the time delay.

Conclusion
In this paper, a Bouc-Wen model is constructed with a genetic identification algorithm to describe the hysteresis property of a smart beam with a piezoelectric actuator. With the identification result, the identified hysteresis loop can be  Shock and Vibration 9 consistent with the experimental hysteresis loop. Based on the hysteresis model as an expected objective model, a stability analysis of the smart system is presented, and adaptive model reference control is designed to explore the influence of the hysteresis characteristic on the control effect. By the simulation and experimental verifications, adaptive model reference control is demonstrated to be feasible for suppressing the vibration of the cantilever beam and the hysteresis item in the expected objective model will affect the stability and control effect of adaptive control for the smart beam with a time delay.

Nomenclature
Bouc-Wen model: A hysteresis model PID: Proportional, integral, and differential control LQR: Linear quadratic regulator Preisach model: A hysteresis model PI model: A hysteresis model MAXGEN: e maximum number of generations NIND: e number of individuals GGAP: e generation gap MRAC: Model reference adaptive control ODE45: A numerical solution method for an ordinary differential equation FFT: Fast Fourier transform PC: Personal computer.

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

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.