Enhanced Adaptive Filtering Algorithm Based on Sliding Mode Control for Active Vibration Rejection of Smart Beam Structures

This article investigates vibration rejection for a continuous smart structure using piezoelectric bimorph patches. Such a structure has inherent nonlinearities, such as hysteresis and creep, and the whole system may experience unexpected disturbances, uncertainties, and noise from external sources. Thus, it is very important to design the active control scheme carefully with adaptive filtering systems to deal with these conditions. An advanced adaptive filtering algorithm was developed based on the conventional least mean squares (LMS) method and sliding mode control for the active vibration rejection system. The sliding mode controller is applied to the standard LMS algorithm to overcome problems with misadjustment and excess error in an optimal manner. A numerical analysis and laboratory experiment show that the technique can significantly attenuate the vibration of the smart structure at different levels and broadband frequency spectra. In addition, unidentified impedance is chosen to change the distribution of the mass, and the robustness and the adaptivity of the proposed approach are verified. The experimental results show that the method can isolate impulse-type vibrations of at least 2.8 dB, even with the adjusted mass arrangement.


Introduction
The primary structural components that experience external forces are usually modeled as continuous cantilevered beams in engineering structures such as vehicles, aircrafts, buildings, and bridges.Such structures commonly have dynamic behavior, which cannot be modeled perfectly and occur from unexpected external sources.These influences could induce undesirable vibrations, resulting in fatigue and likely severe failure.Smart structures with active control methods are one possible solution for isolating these detrimental vibrations, but the control system itself and mistuned control parameters can cause instabilities.Therefore, it is important to design control methods that are inherently insensitive to parameter changes to attenuate structural vibrations.
Lead zirconate titanate (PZT) patches have a number of real-world applications for active control for reducing noise, vibration, and harshness (NVH).Some automotive parts have thin walls and can transmit unexpected noise and vibration from many sources.Thus, there have been many efforts to apply PZT patches to automotive parts such as the body panels [1], muffler system [2], and powertrain oil pan [3].PZT patches are also used for mechanical and aerodynamic NVH reduction applications, such as aircraft fin-tips [4], rotorcraft blades [5], and aeroelastic flutter [6].PZT patches are also used in active fault detections and energy harvesting devices.Piezoelectric shunt circuits can control vibrations from the actuator arms of hard disk drives [7], and they can be used for detecting or fixing cracked beams [8,9].A hybrid-type switching control method has been used in order to isolate unwanted vibration as well as to improve the energy efficiency [10].Active vibration isolation techniques have been implemented with piezoelectric actuators and sensors for the functionally graded material beams [11].A smart structure with shunted piezoelectric patches was constructed, and topology optimization was performed for maximizing the electromechanical coupling factor in terms of the damping efficiency, piezo patch geometries, and patch locations [12].A visualization technique was proposed to identify and reduce the vibrational energy of a 2D plate using a constrained damping layer [13].
To achieve satisfactory vibration attenuation, the structures and vibrational responses should be analyzed.An analytical solution was introduced for a piezoelectric bimorph cantilever beam with suitable boundary conditions [14].The Airy stress function can also be used to study other different piezoelectric cantilevers, including multilayer piezoelectric cantilevers under different loads.The active vibration control of parametrically excited systems with resonance was examined to avoid the parametric instability, and an active control system was introduced, which includes a beam subjected to an axial harmonic load [15].An enhanced method was proposed to stabilize the system via velocity feedback and pole placement by relocating the transition curves.Finally, generally layered smart beams were observed with a one-dimensional finite element method [16].This state space model is based on first-order shear deformation beam theory with Hermite shape functions and the effect of external electric loads.
Appropriate control algorithms must be designed for individual systems to achieve desirable performance.Piezoelectric shunt circuits have been applied extensively for the semi-active control of vibrating structures, but a single piezoelectric shunt circuit can deal with only one vibration mode at a time [17].Many methods have been explored to manage more than one mode, such as multiple piezoelectric shunt circuits [18], synchronized switch voltage damping [19], and periodic piezoelectric shunts [20].However, these techniques are still passive methods, for which adjustment is needed in advance.Thus, a number of active control schemes have been examined.These schemes include positive position feedback [21,22], acceleration feedback [23,24], full-state feedback [25], dynamic hysteresis compensators [26], pole placement control [27], and a hybrid of variable structure control and Lyapunov-based control [28].Decoupled models were proposed to manage many degrees of freedom (DOF), and a backstepping method was introduced to design a robust controller with sliding mode and adaptive strategies [29].To avoid the high computational burden of full-order controllers, a direct minimization method was introduced to design reduced-order H-infinity controllers and genetic algorithms [30].
The conventional approaches for active vibration control are often too complicated and thus difficult to implement.Furthermore, the control targets have to be determined, which limits the tracking, convergence, and stability performance.If a correct mathematical model of flexible structures can be derived, the structures can be regulated actively to isolate the structural vibration.However, several control schemes show performance degradation when the model is not precise enough and the structural dynamics varies with time due to structural modifications.To solve these issues, a simple and robust control algorithm should be developed without any correlation to the original structural dynamics and with low computational burden.
This study uses an enhanced active vibration control method attenuate the vibration of a nonlinear elastic smart structure.A new control scheme called as the sliding mode least mean squares (SM-LMS) algorithm is proposed and validated both numerically and experimentally [31].The control algorithm combines the optimal control characteristics of the least mean squares algorithm with the robustness of sliding mode control (SMC).It has been shown that SM-LMS can be successfully used for the displacement control of piezoelectric actuators [32].However, no study has applied SM-LMS to the active vibration control of flexible structures.To validate the scheme, the controller was implemented on a smart structure constructed with an aluminum beam with two PZT actuators attached to form a bimorph structure and one PZT sensor.For the actual tests, we would not know the exact model of the system due to unknown changes in the operating environment, and the plant may vary over time.This changes the control parameters needed to achieve the best tracking performance.Thus, experiments of two different cases were prepared: (1) control for a baseline structure; and (2) control for a structure with additional dynamics that cannot be predicted.Additional impedance is added to the baseline structure to simulate a sudden change in the plant.The scope is limited to a single-input single-output system, and the first four peaks of the disturbance plant (the natural modes) are considered.Specific objectives of this research are described below.
(1) Introduce the characteristics of the novel nonlinear control methodology (sliding mode least mean squares) and examine the performance on active vibration control through numerical simulations.
(2) Identify the dynamic model of the active beam in the frequency domain with an impact hammer test.Estimate parameters of the linear beam model through a frequency domain identification method.(3) Conduct numerical validations with a conventional approach for attenuating vibration level of the structure under impulsive excitations and observe limitations.(4) Conduct experimental validations and examine the performance of the controller for the baseline structure and the modified structures in order to show the robustness of the new control technique.

Adaptive Filtering Systems for Active Control
Conventionally, digital filters have been applied for modeling dynamic responses to any type of given inputs in engineering structures and acoustic platforms.They apply numerical processes in the discrete time domain for a certain signal to obtain a desired signal shape, and they are broadly used in electronics, computer science, and mechanical engineering.The finite impulse response (FIR) filter is extensively used to represent system responses because of its inherent simplicity, stability, and ease of implementation.In addition, digital filters notably applicable in adaptive structural design and control.Adaptive digital signal processing is applied to update digital filter coefficients automatically in order to create a desired signal from a reference signal and error based on recursive algorithms.Digital filters are thus commonly used in active noise and vibration isolation, system identification, and inverse modeling.An adaptive digital filtering system was designed for active control of the unwanted vibration of continuous structures with a new adaptive filtering algorithm.

Least Mean Squares Method
The least mean squares (LMS) algorithm is the most widely used feedforward control algorithm for adaptive digital filtering systems, since it is relatively simple and easy to implement.The stochastic gradient of the squared error signal is the core of this algorithm [33], and it iteratively regulates the coefficients of transversal filters.These filters use the following two values for the updating formulation: (1) a reference input u[k] containing frequency components to be dealt with; and (2) the measured error value e[k] at the current instant, which is the difference between the filter output y[k] and the disturbance (desired value) d[k] to be tracked.The algorithm adjusts the adaptive digital filter to determine a filter output y[k] that is as close as possible to the disturbance value d [k].A formulation of the LMS algorithm for updating these coefficients can be derived by minimizing the mean squared error with a pre-defined cost function.
The most attractive feature of the LMS method is that it does not require average correlation functions from gradient estimates and matrix inversions [34].It starts from the steepest descent method, and the output of the FIR filter is calculated by multiplying a reference vector u by a vector of filter coefficients or weights w.In Equation (1), u[k − j] is the reference input signal at time index k − j, and w j [k] is the j th filter coefficient at time index k.The desired signal d[k] is the output of the disturbance plant to be tracked and reduced, and the error e[k] is the deviation between the filter output y[k] and the desired signal d A cost function is chosen to minimize the error.The mean squared error ξ is typically selected as a cost function since it is a quadratic function and thus guarantees a unique minimum, which is a very important advantage for the digital filtering system.
The cost function is expanded by substituting Equations ( 1) and (2) into Equation ( 3): P is a cross correlation vector between the reference input and the desired signal, and R is the autocorrelation matrix of the reference signal.The steepest descent concept is used to develop an algorithm that adjusts the vector of filter weights w.The reason why this method is applied here is that it is very straightforward and has been validated in many applications.
Next, Equation ( 7) is partially differentiated with respect to w[k], and then the gradient of the cost function ∇ξ can be derived: From Equations ( 7) and ( 8), the steepest descent algorithm is obtained: where µ is a parameter that regulates the convergence speed and the stability of the adaptation procedure.Figure 1 shows a block diagram of a typical adaptive digital filtering system.
A cost function is chosen to minimize the error.The mean squared error  is typically selected as a cost function since it is a quadratic function and thus guarantees a unique minimum, which is a very important advantage for the digital filtering system. 2 The cost function is expanded by substituting Equations ( 1) and (2) into Equation (3): P is a cross correlation vector between the reference input and the desired signal, and R is the autocorrelation matrix of the reference signal.The steepest descent concept is used to develop an algorithm that adjusts the vector of filter weights w .The reason why this method is applied here is that it is very straightforward and has been validated in many applications.
Next, Equation ( 7) is partially differentiated with respect to [ ] k w , and then the gradient of the cost function   can be derived: From Equations ( 7) and ( 8), the steepest descent algorithm is obtained: where  is a parameter that regulates the convergence speed and the stability of the adaptation procedure.Figure 1 shows a block diagram of a typical adaptive digital filtering system.The gradient of the cost function   is calculated in every iteration with an average of squared errors at recent sampling times.However, it is usually impossible to predict the cross correlation The gradient of the cost function ∇ξ is calculated in every iteration with an average of squared errors at recent sampling times.However, it is usually impossible to predict the cross correlation vector P and the correlation matrix R at every sample time.Thus, to simplify the updating equation for the filter coefficients, the LMS algorithm uses instant estimations for P, R, and the cost function ξ [35].Therefore, e[k] 2 is chosen as an estimate of ξ, and the correlation matrix and cross correlation vector are assumed to be: Next, the estimation of the gradient ∇ξ of the cost function is obtained from Equations ( 8) and ( 10)-( 12): Finally, substituting Equation ( 13) into Equation ( 7), the final form of the LMS algorithm is derived.
Notably, the filter coefficients ŵ[k] cannot trace the actual steepest descent surface because the changes of the weights are based on estimated or simplified gradients, which produce noisy tracking.

Filtered-X Least Mean Squares Method
In most practical applications, there are usually effects from secondary path dynamics due to AD/DA converters, filters, actuators, sensors, amplifiers, and circuit noise, among others.However, we assume that the whole control system includes no such dynamics and that the algorithm will not behave properly with such effects.The most common types of secondary path dynamics are phase shifts or time delays.The whole control structure become unstable or converges to an inappropriate solution under these effects, which changes the filter output f [k] and the signal y[k] subtracted from the desired signal d[k].
The filtered-X LMS (FX-LMS) algorithm is a modified version of the traditional LMS algorithm that manages the secondary path dynamics, which are identified in advance [36].This algorithm requires a filtered reference signal u[k] with a numerical model of the secondary path dynamics Ps (z).The signal is applied in the original LMS algorithm to obtain the adaptive filter output y[k].FX-LMS is more stable than the LMS method when the secondary path dynamics cannot be disregarded.Figure 2 shows a schematic of the FX-LMS algorithm.
vector P and the correlation matrix R at every sample time.Thus, to simplify the updating equation for the filter coefficients, the LMS algorithm uses instant estimations for P , R , and the cost function  [35].Therefore,

2
[ ] e k is chosen as an estimate of  , and the correlation matrix and cross correlation vector are assumed to be: Next, the estimation of the gradient ˆ  of the cost function is obtained from Equations ( 8) and ( 10)-( 12): Finally, substituting Equation ( 13) into Equation ( 7), the final form of the LMS algorithm is derived.
Notably, the filter coefficients [ ] ˆk w cannot trace the actual steepest descent surface because the changes of the weights are based on estimated or simplified gradients, which produce noisy tracking.

Filtered-X Least Mean Squares Method
In most practical applications, there are usually effects from secondary path dynamics due to AD/DA converters, filters, actuators, sensors, amplifiers, and circuit noise, among others.However, we assume that the whole control system includes no such dynamics and that the algorithm will not behave properly with such effects.The most common types of secondary path dynamics are phase shifts or time delays.The whole control structure become unstable or converges to an inappropriate solution under these effects, which changes the filter output [ ] f k and the signal [ ] y k subtracted from the desired signal [ ] d k .
The filtered-X LMS (FX-LMS) algorithm is a modified version of the traditional LMS algorithm that manages the secondary path dynamics, which are identified in advance [36].This algorithm requires a filtered reference signal [ ] k u  with a numerical model of the secondary path dynamics ˆ( ) s P z .The signal is applied in the original LMS algorithm to obtain the adaptive filter output .
FX-LMS is more stable than the LMS method when the secondary path dynamics cannot be disregarded.Figure 2 shows a schematic of the FX-LMS algorithm.[ ] y k The formulation for the FX-LMS algorithm is almost the same as that of the LMS method except that the filtered reference input u[k] is used in the updating equation of the filter weights: The most important property of the filter Ps (z) is that its impulse response contains the same amount of time lag as the secondary path dynamics P s (z).However, a significant defect of FX-LMS is that an offline system identification procedure is essential to obtain Ps (z) beforehand.Although online system identification might be possible, it makes the overall controller complex and difficult to implement.

Limitations
Many studies have modified the LMS algorithm to address defects such as misadjustment and excessive error with respect to the reference inputs.Hybrid structures were proposed by combining the LMS method with feedback control [37,38].The vibration isolation was improved at main frequencies and the speed of the convergence was increased, but the error over a broadband frequency range could not be managed.An alternative feedback structure was suggested with internal model control.This method achieves a better transient response and faster convergence, but an incorrect secondary path model seems to degrade the total control performance [39].An LMS algorithm with adaptive step size was investigated to increase the convergence rate [40,41].Similarly, a variable-step LMS scheme was applied along with a creep model to accomplish rapid convergence and stability and to identify the hysteresis parameters [42].
A variable convergence coefficient corresponding with an input signal can improve the performance, but it may result in noisy step sizes with unwanted disturbances and thus influence the mitigation at other frequencies [43].Another LMS scheme prevents instability with impulsive noise, and the reference is pre-defined to assure stability and convergence, particularly when the magnitude of the signal is larger than a certain threshold [44].An anti-saturation LMS algorithm was introduced to overcome the controller output with large disturbance, which showed improved performance.A different cost function is needed to obtain an enhanced recursive formulation [45].
An extra filter has been used to suppress disturbances in the primary path resulting from uncorrelated noise in the secondary path.This method accomplishes faster convergence, but additional filters increase the complexity and noise in the control system [46].A new LMS algorithm with phase compensation was introduced to deal with signals that include DC components.The traditional LMS may not be effective for such signals, and when a high pass filter is used to eliminate the DC components, it can produce a phase variation [47].

Sliding Mode Control
SMC is one of the most popular nonlinear control methodologies in use today.The strength of this controller is that it is robust, insensitive to plant dynamics, and able to decouple higher-order systems, regardless of their nonlinearity and uncertainty [48].Once a sliding surface is defined, the state trajectory moves between both sides of the surface by high frequency switching control with discontinuous inputs, and it eventually arrives on the pre-defined sliding surface.The motion of a system is not associated with its original dynamics, and it depends on only the sliding mode equation.The main advantage of SMC lies in its robustness and insensitivity to noise, since the system prescribes the dynamics associated with desirable performance characteristics and the system response is insensitive to disturbances or uncertainties in the system parameters [49].
The SMC design procedure includes the following two parts: (1) a discontinuity surface called the sliding surface is chosen to have dynamics related to desirable performance properties; and (2) a high-frequency switching control input is designed to drive the initial system states to the sliding surface.Because the discontinuous input leads the system states towards the sliding surface, the states eventually reach the surface and move along it, and this condition is referred to as "sliding mode."The sliding surface is defined for linear systems: where x is a state vector, and G s is a row vector that decides the dynamics of the sliding mode.The discontinuous inputs u + s (x, t) and u − s (x, t) are then chosen to fulfill the requirements for the existence and the reachability of the sliding mode.A common model of the control input for linear system is: where M is a positive gain that is fairly large to insure s .s < 0 so that the state trajectory will reach the sliding surface.
Figure 3 shows conceptual figures of an ideal sliding mode and a real sliding mode.The system trajectory in the state space (x, . x) arrives at the sliding surface defined by s = 0 and slides towards the origin along the surface.To accomplish ideal sliding mode control, discontinuous control inputs should be able to switch at a fairly high frequency (theoretically infinite), as shown in the figure.However, most mechanical actuators cannot realistically replicate this condition, which results in chattering problems that could possibly make the overall system unstable when combined with unmodeled dynamics.states eventually reach the surface and move along it, and this condition is referred to as "sliding mode."The sliding surface is defined for linear systems: where x is a state vector, and s G is a row vector that decides the dynamics of the sliding mode.
The discontinuous inputs where M is a positive gain that is fairly large to insure 0 ss   so that the state trajectory will reach the sliding surface.
Figure 3 shows conceptual figures of an ideal sliding mode and a real sliding mode.The system trajectory in the state space ( , ) x x  arrives at the sliding surface defined by 0 s  and slides towards the origin along the surface.To accomplish ideal sliding mode control, discontinuous control inputs should be able to switch at a fairly high frequency (theoretically infinite), as shown in the figure.However, most mechanical actuators cannot realistically replicate this condition, which results in chattering problems that could possibly make the overall system unstable when combined with unmodeled dynamics.The SMC should also be considered in the discrete time domain since controllers are usually implemented on modern digital computers.For traditional discrete time sliding mode control (DSMC), the sliding mode is forced in the very next sampling moment [50], the point of which may cause controller saturation.The control system will then not work as expected and can even destabilize the target structure.Therefore, there is a need for a novel control technique that can solve these issues.

Sliding Mode LMS (SM-LMS)
A different structure of the adaptive filtering system is introduced using SMC for active vibration isolation.This new system prevents misadjustment from improper initial states and excessive error problems with respect to the amplitude of the desired values for the sake of the SMC.The sliding mode control is improved by the LMS in terms of smooth convergence and stability.With these properties, it is expected that the technique will enable multi-spectral control of disturbances to manage noise signals and a number of frequency elements simultaneously with optimized tracking schemes.
The SMC technique is combined with the LMS to track or control ( ) x t due to its intrinsic benefits of insensitivity to the plant dynamics, robustness, and ability to decouple systems with high dimension [49].The proposed SM-LMS technique is illustrated in Figure 4.The SMC should also be considered in the discrete time domain since controllers are usually implemented on modern digital computers.For traditional discrete time sliding mode control (DSMC), the sliding mode is forced in the very next sampling moment [50], the point of which may cause controller saturation.The control system will then not work as expected and can even destabilize the target structure.Therefore, there is a need for a novel control technique that can solve these issues.

Sliding Mode LMS (SM-LMS)
A different structure of the adaptive filtering system is introduced using SMC for active vibration isolation.This new system prevents misadjustment from improper initial states and excessive error problems with respect to the amplitude of the desired values for the sake of the SMC.The sliding mode control is improved by the LMS in terms of smooth convergence and stability.With these properties, it is expected that the technique will enable multi-spectral control of disturbances to manage noise signals and a number of frequency elements simultaneously with optimized tracking schemes.
The SMC technique is combined with the LMS to track or control x(t) due to its intrinsic benefits of insensitivity to the plant dynamics, robustness, and ability to decouple systems with high dimension [49].The proposed SM-LMS technique is illustrated in Figure 4. 2 is defined as follows: Appl.Sci.2017, 7, 750 8 of 23 [ ] f k in Figure 2 is defined as follows: Next, the sliding mode [ ] s k in the error state plane is designated as follows.Here, α determines the slope of the sliding surface. [ Using a similar procedure to the derivation of the LMS algorithm, the cost function  for the SM-LMS is now defined as Equation (21), where E means the expected value.Note that the cost function is applied in the LMS algorithm.
A unique minimum value is guaranteed since is a quadratic function.By expanding using Equations ( 4) and ( 20), the following results can be obtained: The vector P indicates the expected cross correlation between the reference and unwanted signals, and the matrix R is the expected autocorrelation of the reference:   where w j [k] indicates the j th filter coefficient at sampling index [k].e[k] is defined with the assumption that P s (z) can be disregarded: Next, the sliding mode s[k] in the error state plane is designated as follows.Here, α determines the slope of the sliding surface.
Using a similar procedure to the derivation of the LMS algorithm, the cost function ξ for the SM-LMS is now defined as Equation (21), where E means the expected value.Note that the cost A unique minimum value is guaranteed since ξ is a quadratic function.By expanding ξ using Equations ( 4) and ( 20), the following results can be obtained: Appl.Sci.2017, 7, 750 9 of 23 The vector P indicates the expected cross correlation between the reference and unwanted signals, and the matrix R is the expected autocorrelation of the reference: The next step is using the steepest descent scheme to obtain an algorithm for regulating the weights w Equation ( 21) is partially differentiated with respect to w[k] in order to derive the gradient of the cost function ∇ξ: For this control algorithm, the gradient ∇ξ along with P and R, have to be determined in every time instant k, as mentioned in the previous section.This is a clumsy process, and to simplify the calculating operations, instantaneous estimates of ξ, P, and R should be obtained at a specific sampling instant in the LMS algorithm [34].Thus estimations of these values indicated by the hat symbol ("ˆ") are defined for time step k: After substituting Equations ( 34) and ( 35) for P and R in Equation (30), a new estimate for ∇ ξ is derived: Upon combining Equations ( 29) and ( 36), w[k] is updated as follows: As shown in Figure 5, the inessential maneuvers of error states are decreased by following a pre-defined sliding surface.Therefore, less input effort or energy would be needed, and faster convergence can be achieved with the introduction of sliding mode control.would accomplish the optimal solution * ŵ after a specific number of iterations when the gradient  becomes zero ideally.Thus, * ŵ can be derived from Equation ( 38): Using the expectation operator on both sides of Equation (37) gives the following equation: The principal axis coordinate should be changed to solve this equation.Hence, the translated variable is defined as ˆ  * V w w , and Equation ( 40) becomes: The coordinates are rotated again with '  V QV to produce the following. where Equation ( 42) can be solved by extending it from time steps 1 k  to 1 , and the higher-power terms of ( ) From Equation (45), it is confirmed that the convergence of the SM-LMS is guaranteed only if the following rules are satisfied:  E[ ŵ[k]] would accomplish the optimal solution ŵ * after a specific number of iterations when the gradient ∇ becomes zero ideally.Thus, ŵ * can be derived from Equation (38): Using the expectation operator on both sides of Equation (37) gives the following equation: The principal axis coordinate should be changed to solve this equation.Hence, the translated variable is defined as V = ŵ − ŵ * , and Equation (40) becomes: The coordinates are rotated again with V = QV to produce the following. where Equation ( 42) can be solved by extending it from time steps k + 1 to 1, and the higher-power terms of (I − µα 2 Λ 1 ) and (−µαΛ 2 ) are ignored.
From Equation (45), it is confirmed that the convergence of the SM-LMS is guaranteed only if the following rules are satisfied: where λ 1,max and λ 2,max indicate the largest eigenvalues of R 1 and R 3 , respectively.Relatively fast convergence could be achieved with the introduction of the sliding mode control for the LMS algorithm.Nonetheless, the possible range of the convergence parameter µ is reduced compared with that of the LMS algorithm, which could compromise the benefits of using the sliding mode control.

Numerical Evaluation of the SM-LMS
A simple numerical validation was carried out to examine the signal tracking performance of the proposed algorithm under the following conditions: (1) a second-order linear time-invariant (LTI) system with a primary frequency of ω n = 30 Hz and damping ratio of ζ = 0.2; (2) a finite impulse response (FIR) filter with three coefficients; and (3) a single-frequency sine wave u(t) = A sin(ωt + φ) with amplitude A = 1, excitation frequency ω = 120 Hz, and no phase (φ = 0).The control parameters are selected as α = 15 and µ = 0.3.Figure 6a shows the time-domain tracking results with the two algorithms.When the signal has sudden changes in the direction and amount, the conventional LMS algorithm has worse performance than the SM-LMS algorithm, and higher frequency fluctuations are observed compared to the estimation error for both cases, as shown in Figure 6b.It is also important to ensure that the filter weights at steady state approach the same or close values for both algorithms, as shown in Figure 7a.The results show that the sliding mode control only supplements the original LMS algorithm to achieve better optimized convergence without changing its primary characteristics.Figure 5b shows the trajectory of the coefficient plane (the filter weights), and the path to the desired steady state is relatively straight and fast for the SM-LMS.Therefore, changes of the filter weights in the SM-LMS case show superior control without inappropriate system movement.
where 1,max  and 2,max  indicate the largest eigenvalues of 1 R and 3 R , respectively.Relatively fast convergence could be achieved with the introduction of the sliding mode control for the LMS algorithm.Nonetheless, the possible range of the convergence parameter  is reduced compared with that of the LMS algorithm, which could compromise the benefits of using the sliding mode control.

Numerical Evaluation of the SM-LMS
A simple numerical validation was carried out to examine the signal tracking performance of the proposed algorithm under the following conditions: (1) a second-order linear time-invariant (LTI) system with a primary frequency of  6a shows the time-domain tracking results with the two algorithms.When the signal has sudden changes in the direction and amount, the conventional LMS algorithm has worse performance than the SM-LMS algorithm, and higher frequency fluctuations are observed compared to the estimation error for both cases, as shown in Figure 6b.It is also important to ensure that the filter weights at steady state approach the same or close values for both algorithms, as shown in Figure 7a.The results show that the sliding mode control only supplements the original LMS algorithm to achieve better optimized convergence without changing its primary characteristics.Figure 5b shows the trajectory of the coefficient plane (the filter weights), and the path to the desired steady state is relatively straight and fast for the SM-LMS.Therefore, changes of the filter weights in the SM-LMS case show superior control without inappropriate system movement., with least mean squares; , with sliding mode least mean squares).
where 1,max  and 2,max  indicate the largest eigenvalues of 1 R and 3 R , respectively.Relatively fast convergence could be achieved with the introduction of the sliding mode control for the LMS algorithm.Nonetheless, the possible range of the convergence parameter  is reduced compared with that of the LMS algorithm, which could compromise the benefits of using the sliding mode control.

Numerical Evaluation of the SM-LMS
A simple numerical validation was carried out to examine the signal tracking performance of the proposed algorithm under the following conditions: (1) a second-order linear time-invariant (LTI) system with a primary frequency of  6a shows the time-domain tracking results with the two algorithms.When the signal has sudden changes in the direction and amount, the conventional LMS algorithm has worse performance than the SM-LMS algorithm, and higher frequency fluctuations are observed compared to the estimation error for both cases, as shown in Figure 6b.It is also important to ensure that the filter weights at steady state approach the same or close values for both algorithms, as shown in Figure 7a.The results show that the sliding mode control only supplements the original LMS algorithm to achieve better optimized convergence without changing its primary characteristics.Figure 5b shows the trajectory of the coefficient plane (the filter weights), and the path to the desired steady state is relatively straight and fast for the SM-LMS.Therefore, changes of the filter weights in the SM-LMS case show superior control without inappropriate system movement., with least mean squares; , with sliding mode least mean squares).
, original signal; where 1,max  and 2,max  indicate the largest eigenvalues of 1 R and 3 R , respectively.Relatively fast convergence could be achieved with the introduction of the sliding mode control for the LMS algorithm.Nonetheless, the possible range of the convergence parameter  is reduced compared with that of the LMS algorithm, which could compromise the benefits of using the sliding mode control.

Numerical Evaluation of the SM-LMS
A simple numerical validation was carried out to examine the signal tracking performance of the proposed algorithm under the following conditions: (1) a second-order linear time-invariant (LTI) system with a primary frequency of  6a shows the time-domain tracking results with the two algorithms.When the signal has sudden changes in the direction and amount, the conventional LMS algorithm has worse performance than the SM-LMS algorithm, and higher frequency fluctuations are observed compared to the estimation error for both cases, as shown in Figure 6b.It is also important to ensure that the filter weights at steady state approach the same or close values for both algorithms, as shown in Figure 7a.The results show that the sliding mode control only supplements the original LMS algorithm to achieve better optimized convergence without changing its primary characteristics.Figure 5b shows the trajectory of the coefficient plane (the filter weights), and the path to the desired steady state is relatively straight and fast for the SM-LMS.Therefore, changes of the filter weights in the SM-LMS case show superior control without inappropriate system movement.
(a) (b) where 1,max  and 2,max  indicate the largest eigenvalues of 1 R and 3 R , respectively.Relatively fast convergence could be achieved with the introduction of the sliding mode control for the LMS algorithm.Nonetheless, the possible range of the convergence parameter  is reduced compared with that of the LMS algorithm, which could compromise the benefits of using the sliding mode control.

Numerical Evaluation of the SM-LMS
A simple numerical validation was carried out to examine the signal tracking performance of the proposed algorithm under the following conditions: (1) a second-order linear time-invariant (LTI) system with a primary frequency of  6a shows the time-domain tracking results with the two algorithms.When the signal has sudden changes in the direction and amount, the conventional LMS algorithm has worse performance than the SM-LMS algorithm, and higher frequency fluctuations are observed compared to the estimation error for both cases, as shown in Figure 6b.It is also important to ensure that the filter weights at steady state approach the same or close values for both algorithms, as shown in Figure 7a.The results show that the sliding mode control only supplements the original LMS algorithm to achieve better optimized convergence without changing its primary characteristics.Figure 5b shows the trajectory of the coefficient plane (the filter weights), and the path to the desired steady state is relatively straight and fast for the SM-LMS.Therefore, changes of the filter weights in the SM-LMS case show superior control without inappropriate system movement.
(a) (b) , with least mean squares; , with sliding mode least mean squares).
, with least mean squares; Appl.Sci.2017, 7, 750 11 of 23 where 1,max  and 2,max  indicate the largest eigenvalues of 1 R and 3 R , respectively.Relatively fast convergence could be achieved with the introduction of the sliding mode control for the LMS algorithm.Nonetheless, the possible range of the convergence parameter  is reduced compared with that of the LMS algorithm, which could compromise the benefits of using the sliding mode control.

Numerical Evaluation of the SM-LMS
A simple numerical validation was carried out to examine the signal tracking performance of the proposed algorithm under the following conditions: (1) a second-order linear time-invariant (LTI) system with a primary frequency of  6a shows the time-domain tracking results with the two algorithms.When the signal has sudden changes in the direction and amount, the conventional LMS algorithm has worse performance than the SM-LMS algorithm, and higher frequency fluctuations are observed compared to the estimation error for both cases, as shown in Figure 6b.It is also important to ensure that the filter weights at steady state approach the same or close values for both algorithms, as shown in Figure 7a.The results show that the sliding mode control only supplements the original LMS algorithm to achieve better optimized convergence without changing its primary characteristics.Figure 5b shows the trajectory of the coefficient plane (the filter weights), and the path to the desired steady state is relatively straight and fast for the SM-LMS.Therefore, changes of the filter weights in the SM-LMS case show superior control without inappropriate system movement.
(a) (b) , with least mean squares; , with sliding mode least mean squares).
, with sliding mode least mean squares).

Problem Formulation
A one-dimensional continuous structure was prepared for the validation of the active control strategy.The structure is a thin, aluminum, cantilevered elastic beam with piezoelectric ceramic patches (PZT) attached to both sides.The specific type of PZT material used is PZT-5A and "31" mode has been used for the actuator.It is better for actuator applications for dynamic operating conditions, while it is still good for utilized as sensor applications.The beam can be bent when the ceramic has deformations with applied voltage, which is called a "bimorph".Figure 8 shows the concept of different types of bimorph structures.The red arrows show the polarization direction, and green arrows indicate the strain direction.Table 1 shows the material properties of the PZT patch., least mean squares; , sliding mode least mean squares).

Problem Formulation
A one-dimensional continuous structure was prepared for the validation of the active control strategy.The structure is a thin, aluminum, cantilevered elastic beam with piezoelectric ceramic patches (PZT) attached to both sides.The specific type of PZT material used is PZT-5A and "31" mode has been used for the actuator.It is better for actuator applications for dynamic operating conditions, while it is still good for utilized as sensor applications.The beam can be bent when the ceramic has deformations with applied voltage, which is called a "bimorph".Figure 8 shows the concept of different types of bimorph structures.The red arrows show the polarization direction, and green arrows indicate the strain direction.Table 1 shows the material properties of the PZT patch., least mean squares; , sliding mode least mean squares).

Problem Formulation
A one-dimensional continuous structure was prepared for the validation of the active control strategy.The structure is a thin, aluminum, cantilevered elastic beam with piezoelectric ceramic patches (PZT) attached to both sides.The specific type of PZT material used is PZT-5A and "31" mode has been used for the actuator.It is better for actuator applications for dynamic operating conditions, while it is still good for utilized as sensor applications.The beam can be bent when the ceramic has deformations with applied voltage, which is called a "bimorph".Figure 8 shows the concept of different types of bimorph structures.The red arrows show the polarization direction, and green arrows indicate the strain direction.Table 1 shows the material properties of the PZT patch., least mean squares; , sliding mode least mean squares).

Problem Formulation
A one-dimensional continuous structure was prepared for the validation of the active control strategy.The structure is a thin, aluminum, cantilevered elastic beam with piezoelectric ceramic patches (PZT) attached to both sides.The specific type of PZT material used is PZT-5A and "31" mode has been used for the actuator.It is better for actuator applications for dynamic operating conditions, while it is still good for utilized as sensor applications.The beam can be bent when the ceramic has deformations with applied voltage, which is called a "bimorph".Figure 8 shows the concept of different types of bimorph structures.The red arrows show the polarization direction, and green arrows indicate the strain direction.Table 1 shows the material properties of the PZT patch.control only supplements the original LMS algorithm to achieve better optimized convergence without changing its primary characteristics.Figure 5b shows the trajectory of the coefficient plane (the filter weights), and the path to the desired steady state is relatively straight and fast for the SM-LMS.Therefore, changes of the filter weights in the SM-LMS case show superior control without inappropriate system movement., with least mean squares; , with sliding mode least mean squares).
, least mean squares; ontrol only supplements the original LMS algorithm to achieve better optimized convergence ithout changing its primary characteristics.Figure 5b shows the trajectory of the coefficient plane he filter weights), and the path to the desired steady state is relatively straight and fast for the SM-MS.Therefore, changes of the filter weights in the SM-LMS case show superior control without appropriate system movement., with least mean squares; , with sliding mode least mean squares).

Problem Formulation
A one-dimensional continuous structure was prepared for the validation of the active control strategy.The structure is a thin, aluminum, cantilevered elastic beam with piezoelectric ceramic patches (PZT) attached to both sides.The specific type of PZT material used is PZT-5A and "31" mode has been used for the actuator.It is better for actuator applications for dynamic operating conditions, while it is still good for utilized as sensor applications.The beam can be bent when the ceramic has deformations with applied voltage, which is called a "bimorph".Figure 8 shows the concept of different types of bimorph structures.The red arrows show the polarization direction, and green arrows indicate the strain direction.Table 1 shows the material properties of the PZT patch., least mean squares; , sliding mode least mean squares).

Problem Formulation
A one-dimensional continuous structure was prepared for the validation of the active control strategy.The structure is a thin, aluminum, cantilevered elastic beam with piezoelectric ceramic patches (PZT) attached to both sides.The specific type of PZT material used is PZT-5A and "31" mode has been used for the actuator.It is better for actuator applications for dynamic operating conditions, while it is still good for utilized as sensor applications.The beam can be bent when the ceramic has deformations with applied voltage, which is called a "bimorph".Figure 8 shows the concept of different types of bimorph structures.The red arrows show the polarization direction, and green arrows indicate the strain direction.Table 1 shows the material properties of the PZT patch.Figure 9 shows a schematic of the experimental setup.Serial bimorph PZT actuators are attached to the surface of the aluminum beam, and an additional small piezoelectric ceramic is attached to one side as a strain sensor for the feedback control system.PZT sensors are selected to check their feasibility for the active vibration control applications.An accelerometer is attached at the tip of the beam as a performance indicator.A bench vise and an electrodynamic shaker are installed on a vibration isolation table, and the whole beam structure is gripped with the vise, resulting clamped-free boundary conditions.The electrodynamic shaker provides impulsive force to the vise, which results in lateral vibration of the beam that has to be mitigated.Additional mass can be attached to the right side of the beam, which works as an impedance change of the structure.
Appl.Sci.2017, 7, 750 13 of 23 Figure 9 shows a schematic of the experimental setup.Serial bimorph PZT actuators are attached to the surface of the aluminum beam, and an additional small piezoelectric ceramic is attached to one side as a strain sensor for the feedback control system.PZT sensors are selected to check their feasibility for the active vibration control applications.An accelerometer is attached at the tip of the beam as a performance indicator.A bench vise and an electrodynamic shaker are installed on a vibration isolation table, and the whole beam structure is gripped with the vise, resulting clampedfree boundary conditions.The electrodynamic shaker provides impulsive force to the vise, which results in lateral vibration of the beam that has to be mitigated.Additional mass can be attached to the right side of the beam, which works as an impedance change of the structure.A data acquisition and control system (dSPACE 1104) is used to implement real-time control, and it is connected to a desktop computer with MATLAB/Simulink installed.The designed controllers are uploaded to the dSPACE module after compiling the MATLAB files, and the control system is prepared as follows: (1) signals from the PZT sensor and the accelerometer are acquired by ADC channels and used for the control; and (2) the control input calculated in the dSPACE module is produced from a DAC channel and provided to the PZT actuators after amplification. Figure 10 shows the main components of the experimental setup.The following frequency response function was obtained from the impact test, and the relationship between the input and output (the transfer function) can be predicted from the results.The objective is to investigate the controller performance with this baseline structure and a modified structure with additional mass, which dramatically changes the structural dynamics and transfer function.Figure 11 compares the frequency response functions of the baseline and the modified structures, and obvious changes in natural frequencies can be observed.A data acquisition and control system (dSPACE 1104) is used to implement real-time control, and it is connected to a desktop computer with MATLAB/Simulink installed.The designed controllers are uploaded to the dSPACE module after compiling the MATLAB files, and the control system is prepared as follows: (1) signals from the PZT sensor and the accelerometer are acquired by ADC channels and used for the control; and (2) the control input calculated in the dSPACE module is produced from a DAC channel and provided to the PZT actuators after amplification. Figure 10 shows the main components of the experimental setup.
Appl.Sci.2017, 7, 750 13 of 23 Figure 9 shows a schematic of the experimental setup.Serial bimorph PZT actuators are attached to the surface of the aluminum beam, and an additional small piezoelectric ceramic is attached to one side as a strain sensor for the feedback control system.PZT sensors are selected to check their feasibility for the active vibration control applications.An accelerometer is attached at the tip of the beam as a performance indicator.A bench vise and an electrodynamic shaker are installed on a vibration isolation table, and the whole beam structure is gripped with the vise, resulting clampedfree boundary conditions.The electrodynamic shaker provides impulsive force to the vise, which results in lateral vibration of the beam that has to be mitigated.Additional mass can be attached to the right side of the beam, which works as an impedance change of the structure.A data acquisition and control system (dSPACE 1104) is used to implement real-time control, and it is connected to a desktop computer with MATLAB/Simulink installed.The designed controllers are uploaded to the dSPACE module after compiling the MATLAB files, and the control system is prepared as follows: (1) signals from the PZT sensor and the accelerometer are acquired by ADC channels and used for the control; and (2) the control input calculated in the dSPACE module is produced from a DAC channel and provided to the PZT actuators after amplification. Figure 10 shows the main components of the experimental setup.The following frequency response function was obtained from the impact test, and the relationship between the input and output (the transfer function) can be predicted from the results.The objective is to investigate the controller performance with this baseline structure and a modified structure with additional mass, which dramatically changes the structural dynamics and transfer function.Figure 11 compares the frequency response functions of the baseline and the modified structures, and obvious changes in natural frequencies can be observed.The following frequency response function was obtained from the impact test, and the relationship between the input and output (the transfer function) can be predicted from the results.The objective is to investigate the controller performance with this baseline structure and a modified structure with additional mass, which dramatically changes the structural dynamics and transfer function.Figure 11 compares the frequency response functions of the baseline and the modified structures, and obvious changes in natural frequencies can be observed.

System Identification and Notch Filter Design
To design an appropriate controller, the structural dynamics of the system should be fully characterized.A frequency-domain system identification method was applied to derive the model of the system.Figure 12 shows the frequency response function of the structural model to be identified.The damping ratio of each peak can be estimated by the quadrature peak picking method.It is a classical method for determining the damping at resonance, and half power (−3 dB) points of each resonance magnitude should be identified.For a particular mode, the damping ratio m  can be found from Equation (49). Figure 13 illustrates the concept of this method.
The accuracy of this method depends on the frequency resolution used for the measurement, since it determines how accurately the peak magnitude can be measured.For lightly damped structures, high-resolution analysis is required to measure the peak exactly.Consequently, repeated measurements at each resonance frequency are normally required to achieve sufficient accuracy.The obtained resonances are shown in Table 2.

System Identification and Notch Filter Design
To design an appropriate controller, the structural dynamics of the system should be fully characterized.A frequency-domain system identification method was applied to derive the model of the system.Figure 12 shows the frequency response function of the structural model to be identified.The damping ratio of each peak can be estimated by the quadrature peak picking method.It is a classical method for determining the damping at resonance, and half power (−3 dB) points of each resonance magnitude should be identified.

System Identification and Notch Filter Design
To design an appropriate controller, the structural dynamics of the system should be fully characterized.A frequency-domain system identification method was applied to derive the model of the system.Figure 12 shows the frequency response function of the structural model to be identified.The damping ratio of each peak can be estimated by the quadrature peak picking method.It is a classical method for determining the damping at resonance, and half power (−3 dB) points of each resonance magnitude should be identified.For a particular mode, the damping ratio m  can be found from Equation (49). Figure 13 illustrates the concept of this method.
The accuracy of this method depends on the frequency resolution used for the measurement, since it determines how accurately the peak magnitude can be measured.For lightly damped structures, high-resolution analysis is required to measure the peak exactly.Consequently, repeated measurements at each resonance frequency are normally required to achieve sufficient accuracy.The obtained resonances are shown in Table 2.For a particular mode, the damping ratio ζ m can be found from Equation (49). Figure 13 illustrates the concept of this method.
The accuracy of this method depends on the frequency resolution used for the measurement, since it determines how accurately the peak magnitude can be measured.For lightly damped structures, high-resolution analysis is required to measure the peak exactly.Consequently, repeated measurements at each resonance frequency are normally required to achieve sufficient accuracy.The obtained resonances are shown in Table 2.A transfer function can be identified for the object structure based on the frequency response data.The least square estimation of the linear system parameters is presented as a complex-domain formulation [51].Figure 14 shows a general model of the frequency-domain identification of linear systems.The system is represented by its transfer function ( ) H s has a rational form that is extended by a delay term The excitation signal is represented by a Fourier series, with complex amplitudes i U at angular frequencies i  .The response of the system is . Both the measured input and output complex amplitudes are corrupted by noise terms u n and y n , respectively (errors in variables model), which are usually assumed to be Gaussian, uncorrelated between the input and output, and uncorrelated between different frequency points.Based on the above assumptions, a maximumlikelihood estimate of the system parameters is developed in real terms.A transfer function can be identified for the object structure based on the frequency response data.The least square estimation of the linear system parameters is presented as a complex-domain formulation [51].Figure 14 shows a general model of the frequency-domain identification of linear systems.The system is represented by its transfer function H(s), where s = jω in the Laplace domain.H(s) has a rational form that is extended by a delay term e −jωt d : The excitation signal is represented by a Fourier series, with complex amplitudes U i at angular frequencies ω i .The response of the system is Y i = H(ω i )U i .Both the measured input and output complex amplitudes are corrupted by noise terms n u and n y , respectively (errors in variables model), which are usually assumed to be Gaussian, uncorrelated between the input and output, and uncorrelated between different frequency points.Based on the above assumptions, a maximum-likelihood estimate of the system parameters is developed in real terms.A transfer function can be identified for the object structure based on the frequency response data.The least square estimation of the linear system parameters is presented as a complex-domain formulation [51].Figure 14 shows a general model of the frequency-domain identification of linear systems.The system is represented by its transfer function The excitation signal is represented by a Fourier series, with complex amplitudes i U at angular frequencies i  .The response of the system is . Both the measured input and output complex amplitudes are corrupted by noise terms u n and y n , respectively (errors in variables model), which are usually assumed to be Gaussian, uncorrelated between the input and output, and uncorrelated between different frequency points.Based on the above assumptions, a maximumlikelihood estimate of the system parameters is developed in real terms.The time-domain measurements are transformed into the frequency domain, and the calculated complex input and output amplitudes at selected angular frequencies ω i , i = 1, • • • , G are U mi and Y mi , respectively.The unknown parameters are the unknown coefficients of the numerator and the denominator of the transfer function (Equation ( 49)), which are collected in the vector C. Furthermore, the exact complex input and output amplitudes U and Y are also unknown.The basic equations can be written as follows for the exact values (Equation ( 50)) and for the noisy observations (Equation ( 51)): The maximum-likelihood estimate of the unknown quantities is equivalent to the minimization of the following cost function: which is subject to the constraints The upper bar means the complex conjugate, and σ ui and σ yi are the standard deviations of the input and output, respectively.Using this above methodology, a transfer function with 9 zeros and 11 poles has been identified, which has four major peaks up to around 500 Hz.These four peaks should be attenuated with proper control methods.Figure 15 compares the frequency response functions from the experiment and the system identification, and Equation (53) shows the formulation of the transfer function.This result seems to fit with the experimental results for up to 500 Hz. Figure 15 also shows the root locus for checking which poles are related to each peak.Some parts of the locus go over the right half plane, which indicates instability.
Appl.Sci.2017, 7, 750 16 of 23 Furthermore, the exact complex input and output amplitudes U and Y are also unknown.The basic equations can be written as follows for the exact values (Equation ( 50)) and for the noisy observations (Equation ( 51)): ( , ) ( , ) ( ) The maximum-likelihood estimate of the unknown quantities is equivalent to the minimization of the following cost function: which is subject to the constraints ( , ) The upper bar means the complex conjugate, and ui  and yi  are the standard deviations of the input and output, respectively.Using this above methodology, a transfer function with 9 zeros and 11 poles has been identified, which has four major peaks up to around 500 Hz.These four peaks should be attenuated with proper control methods.Figure 15 compares the frequency response functions from the experiment and the system identification, and Equation (53) shows the formulation of the transfer function.This result seems to fit with the experimental results for up to 500 Hz.
With the identified poles and zeros of the system, notch filters are designed to reduce the amplitude of each peak.A stop band filter (BSF) or a notch filter refers to a band-elimination or blocking filter that exhibits sharp attenuation characteristics in only a specific frequency band.Peak amplitudes can be reduced by adding zeros near poles to maintain the DC gain.Since two zeros are added for each peak, two poles should also be added, which is similar to replacing the damping ratio at each resonance frequency.These are the resonance frequencies and damping ratios of each peak.The filters are derived by maintaining n  and changing the damping ratio  to 0.7.This is done to confirm that an actual solution exists and that enough control authority can be obtained from the actuators for the actual system.Table 3 shows the transfer functions of the notch filters based on Equation (54), in which the original damping ratio z  replaces the higher damping ratio p  at each resonance frequency so that the amplitudes at each peak can be attenuated.G(s) = With the identified poles and zeros of the system, notch filters are designed to reduce the amplitude of each peak.A stop band filter (BSF) or a notch filter refers to a band-elimination or blocking filter that exhibits sharp attenuation characteristics in only a specific frequency band.Peak amplitudes can be reduced by adding zeros near poles to maintain the DC gain.Since two zeros are added for each peak, two poles should also be added, which is similar to replacing the damping ratio at each resonance frequency.These are the resonance frequencies and damping ratios of each peak.The filters are derived by maintaining ω n and changing the damping ratio ζ to 0.7.This is done to confirm that an actual solution exists and that enough control authority can be obtained from the actuators for the actual system.Table 3 shows the transfer functions of the notch filters based on Equation (54), in which the original damping ratio ζ z replaces the higher damping ratio ζ p at each resonance frequency so that the amplitudes at each peak can be attenuated.As shown in Figure 16, the amplitudes of the peaks are significantly reduced compared to the original system, and the root locus of the filtered system has been changed.The part in the right-hand plane is automatically removed by these filters with the added zeros near poles.As shown in Figure 16, the amplitudes of the peaks are significantly reduced compared to the original system, and the root locus of the filtered system has been changed.The part in the right-hand plane is automatically removed by these filters with the added zeros near poles.The control results are also described in the time domain.Numerical simulations were conducted with impulse and step inputs to the system.Figures 17 and 18 show the impulse response and the step response of the structure.The notch filters show dramatically reduced vibration amplitudes compared to the case with no control.However, these filters will work on only the baseline structure, assuming that the identified model is accurate.The control results are also described in the time domain.Numerical simulations were conducted with impulse and step inputs to the system.Figures 17 and 18 show the impulse response and the step response of the structure.The notch filters show dramatically reduced vibration amplitudes compared to the case with no control.However, these filters will work on only the baseline structure, assuming that the identified model is accurate.In a real situation, the exact model of the system cannot be known due to changes in the operational environment, and the plant may vary over time.To incorporate the fact that the plant will change and to highlight the efficacy of the controller, it is assumed that an additional mass (about a 20% increase) is attached to the original system, as mentioned in Section 4.
Since an added mass implies decreased resonance frequencies, a modified transfer function would exist with reduced resonance frequencies for each peak.In the next section, the SM-LMS algorithm's target-tracking performance is observed for vibration attenuation using the baseline and modified beam structures.

Experiments
Several tests were performed with the cantilevered flexible beam.The tests are divided into two cases of vibration reduction control: one for the baseline structure of the beam, and one for a beam with variable mass.The LMS and SM-LMS algorithms are applied to compare the vibration reduction performance.The control parameters were tuned to obtain the best results.Figure 19 shows the control results for the baseline beam structure.The plot on the left indicates the strain change of the aluminum beam in the time domain from the PZT sensor outputs.Note that the vibration stops in 0.2 to 0.3 s with LMS and SM-LMS, whereas it takes longer than 1 s to stop without a controller.The plot on the right compares the frequency response plots between open-loop and closed-loop cases, and some amplitude reductions can be seen for each peak frequency.In a real situation, the exact model of the system cannot be known due to changes in the operational environment, and the plant may vary over time.To incorporate the fact that the plant will change and to highlight the efficacy of the controller, it is assumed that an additional mass (about a 20% increase) is attached to the original system, as mentioned in Section 4.
Since an added mass implies decreased resonance frequencies, a modified transfer function would exist with reduced resonance frequencies for each peak.In the next section, the SM-LMS algorithm's target-tracking performance is observed for vibration attenuation using the baseline and modified beam structures.

Experiments
Several tests were performed with the cantilevered flexible beam.The tests are divided into two cases of vibration reduction control: one for the baseline structure of the beam, and one for a beam with variable mass.The LMS and SM-LMS algorithms are applied to compare the vibration reduction performance.The control parameters were tuned to obtain the best results.Figure 19 shows the control results for the baseline beam structure.The plot on the left indicates the strain change of the aluminum beam in the time domain from the PZT sensor outputs.Note that the vibration stops in 0.2 to 0.3 s with LMS and SM-LMS, whereas it takes longer than 1 s to stop without a controller.The plot on the right compares the frequency response plots between open-loop and closed-loop cases, and some amplitude reductions can be seen for each peak frequency.In a real situation, the exact model of the system cannot be known due to changes in the operational environment, and the plant may vary over time.To incorporate the fact that the plant will change and to highlight the efficacy of the controller, it is assumed that an additional mass (about a 20% increase) is attached to the original system, as mentioned in Section 4.
Since an added mass implies decreased resonance frequencies, a modified transfer function would exist with reduced resonance frequencies for each peak.In the next section, the SM-LMS algorithm's target-tracking performance is observed for vibration attenuation using the baseline and modified beam structures.

Experiments
Several tests were performed with the cantilevered flexible beam.The tests are divided into two cases of vibration reduction control: one for the baseline structure of the beam, and one for a beam with variable mass.The LMS and SM-LMS algorithms are applied to compare the vibration reduction performance.The control parameters were tuned to obtain the best results.Figure 19 shows the control results for the baseline beam structure.The plot on the left indicates the strain change of the aluminum beam in the time domain from the PZT sensor outputs.Note that the vibration stops in 0.2 to 0.3 s with LMS and SM-LMS, whereas it takes longer than 1 s to stop without a controller.The plot on the right compares the frequency response plots between open-loop and closed-loop cases, and some amplitude reductions can be seen for each peak frequency.After 0.3 s, the vibration amplitudes with both control algorithms are very similar, but LMS shows worse vibration isolation performance than SM-LMS.The frequency response shows a more obvious performance difference between the two methods, as shown in Table 4.There is excellent attenuation of mode 1 and slight attenuation of mode 2 with LMS, while higher frequencies show slight amplifications.SM-LMS shows excellent attenuation of modes 1 and 2 and slight attenuation of mode 4, while mode 3 shows slight amplification.At every peak, the difference in attenuated amplitude between the two methods is about 3 to 5 dBs.The amplitudes of 1st, 2nd, and 4th modal frequencies are reduced, and the 3rd modal frequency has some spill-over with SM-LMS.However, with LMS, the 3rd and 4th peaks have spillovers, and the amounts are higher than in the case with SM-LMS.When observing the vibration level in broadband, the SM-LMS shows better vibration isolation performance, thus reducing the overall RMS value.
Next, to incorporate the fact that the plant will change and to highlight the efficacy of the controller, an additional mass is attached to the original system (increasing the total weight by about a 20%).A bracket is fastened near the tip of the beam, as shown in Figure 11.Since an added mass implies decreased resonance frequencies, a modified frequency response with reduced resonance frequencies for each peak is obtained to measure the performance of the controller in these cases.Note that the SM-LMS is not related to the original structural dynamics and can make the system have desired dynamics.Moreover, the parameters of the SM-LMS used for the baseline structures are maintained for this case of control in order to observe the robustness of the controller.
Again, the time-domain responses caused by impulsive disturbances show significant reduction visible, as shown in Figure 20.Furthermore, the amplitudes of each resonance peak and other frequencies are decreased, regardless of the significant change in the plant dynamics.Similarly, the vibration stops in 0.4 to 0.5 s with LMS and SM-LMS, whereas it takes longer than 1 s without a controller.The plot on the right compares the frequency response between the open-loop and closedloop cases, in which some amplitude reductions can be seen for each peak frequency.After 0.3 s, the vibration amplitudes with both control algorithms are very similar, but LMS shows worse vibration isolation performance than SM-LMS.The frequency response shows a more obvious performance difference between the two methods, as shown in Table 4.There is excellent attenuation of mode 1 and slight attenuation of mode 2 with LMS, while higher frequencies show slight amplifications.SM-LMS shows excellent attenuation of modes 1 and 2 and slight attenuation of mode 4, while mode 3 shows slight amplification.At every peak, the difference in attenuated amplitude between the two methods is about 3 to 5 dBs.The amplitudes of 1st, 2nd, and 4th modal frequencies are reduced, and the 3rd modal frequency has some spill-over with SM-LMS.However, with LMS, the 3rd and 4th peaks have spillovers, and the amounts are higher than in the case with SM-LMS.When observing the vibration level in broadband, the SM-LMS shows better vibration isolation performance, thus reducing the overall RMS value.
Next, to incorporate the fact that the plant will change and to highlight the efficacy of the controller, an additional mass is attached to the original system (increasing the total weight by about a 20%).A bracket is fastened near the tip of the beam, as shown in Figure 11.Since an added mass implies decreased resonance frequencies, a modified frequency response with reduced resonance frequencies for each peak is obtained to measure the performance of the controller in these cases.Note that the SM-LMS is not related to the original structural dynamics and can make the system have desired dynamics.Moreover, the parameters of the SM-LMS used for the baseline structures are maintained for this case of control in order to observe the robustness of the controller.
Again, the time-domain responses caused by impulsive disturbances show significant reduction visible, as shown in Figure 20.Furthermore, the amplitudes of each resonance peak and other frequencies are decreased, regardless of the significant change in the plant dynamics.Similarly, the vibration stops in 0.4 to 0.5 s with LMS and SM-LMS, whereas it takes longer than 1 s without a controller.The plot on the right compares the frequency response between the open-loop and closed-loop cases, in which some amplitude reductions can be seen for each peak frequency.After 0.5 s, the vibration amplitudes for both control algorithms are very similar, but LMS shows worse vibration isolation performance than SM-LMS.The frequency response shows a more obvious performance difference between the two methods, as indicated in Table 5. LMS shows excellent attenuation of mode 1 but has very little effect on higher modes.SM-LMS shows excellent attenuation of mode 1 and slight attenuation of modes 2 and 3, while mode 4 shows very slight amplification.At every peak, the difference in attenuated amplitude between the two methods is about 2 to 3 dBs.The amplitudes of the 1st, 2nd, and 3rd modal frequencies are reduced, and the 4th modal frequency has some spill-over with SM-LMS.However, with LMS, the 2nd and 4th peaks have spillovers, and the amounts are higher than in the case with SM-LMS.When observing the vibration level in broadband, the SM-LMS shows better vibration isolation performance, thus showing a reduced overall RMS although it is 1.1 dB lower than the baseline case.As shown in the simulation results, the time-domain responses caused by impulsive disturbances are significantly reduced, and the amplitude of each resonance peak is nearly flattened.These results verify that the SM-LMS method can reduce the vibration of the cantilevered beam considerably, even though the structural dynamics is dramatically changed with the additional impedance.

Conclusions and Future Work
An enhanced LMS algorithm has been proposed and applied to an adaptive digital filtering system combined with sliding mode control.The method benefits from the nonlinear control technique inside the control system structure for better adaptability and tracking performance.This adaptive filtering system was used for the vibration reduction and control of a continuous smart structure with piezoelectric patches attached to it.The SM-LMS method allows for significant and robust vibration reduction for not only the original structure, but also for a beam with modified dynamics.Thus, it can be a good solution for reducing the vibration of structures with unknown dynamics.After 0.5 s, the vibration amplitudes for both control algorithms are very similar, but LMS shows worse vibration isolation performance than SM-LMS.The frequency response shows a more obvious performance difference between the two methods, as indicated in Table 5. LMS shows excellent attenuation of mode 1 but has very little effect on higher modes.SM-LMS shows excellent attenuation of mode 1 and slight attenuation of modes 2 and 3, while mode 4 shows very slight amplification.At every peak, the difference in attenuated amplitude between the two methods is about 2 to 3 dBs.The amplitudes of the 1st, 2nd, and 3rd modal frequencies are reduced, and the 4th modal frequency has some spill-over with SM-LMS.However, with LMS, the 2nd and 4th peaks have spillovers, and the amounts are higher than in the case with SM-LMS.When observing the vibration level in broadband, the SM-LMS shows better vibration isolation performance, thus showing a reduced overall RMS value, although it is 1.1 dB lower than the baseline case.As shown in the simulation results, the time-domain responses caused by impulsive disturbances are significantly reduced, and the amplitude of each resonance peak is nearly flattened.These results verify that the SM-LMS method can reduce the vibration of the cantilevered beam considerably, even though the structural dynamics is dramatically changed with the additional impedance.

Conclusions and Future Work
An enhanced LMS algorithm has been proposed and applied to an adaptive digital filtering system combined with sliding mode control.The method benefits from the nonlinear control technique inside the control system structure for better adaptability and tracking performance.This adaptive filtering system was used for the vibration reduction and control of a continuous smart structure with piezoelectric patches attached to it.The SM-LMS method allows for significant and robust vibration reduction for not only the original structure, but also for a beam with modified dynamics.Thus, it can be a good solution for reducing the vibration of structures with unknown dynamics.
The contributions of this research are summarized as: First, the proposed control system has improved stability, robustness, and capability to manage signals with multiple spectra with increasing complexity and frequency of components to be controlled in the output of the disturbance plant.Numerical analyses demonstrated that the method has adequate tracking performance.Second, a smart structure with piezoelectric patches was constructed and its dynamic model has been identified with frequency domain identification method; third, a conventional control method has been validated numerically in the case of active vibration control and showed limitations.Third, the method was validated for the original structure and a modified structure.The controller had acceptable tracking performance in both cases.
This work could be validated with several real-world applications in active vibration control, such as vehicle vibration reduction and vibration attenuation for engineering structures.In addition, a hybrid controller could be proposed to control two or more directions simultaneously for realistic applications by effectively managing an extensive range of multi-spectral signals.Self-sensing actuator could be implemented with one PZT material, which acts as sensor and actuator simultaneously.Intelligent control methods such as neural networks could also be more adaptive and robust to changeable system dynamics.

Figure 3 .
Figure 3. Phase trajectory in sliding mode: (a) ideal sliding mode; and (b) real sliding mode.

Figure 3 .
Figure 3. Phase trajectory in sliding mode: (a) ideal sliding mode; and (b) real sliding mode.

Figure 5 . 2 [
Figure 5. Error surface illustration for: (a) the conventional LMS method; and (b) SM-LMS.We now discuss the convergence of ˆ[ ] k w with the optimized solutions of SM-LMS when the cost function

Figure 6 .
Figure 6.Control of a single sinusoid at 120 Hz in the time domain: (a) predicted with least mean squares and sliding mode least mean squares (Key: , original signal; , filter output); and (b) predicted estimation error (Key:, with least mean squares; , with sliding mode least mean squares).

Figure 6 .
Figure 6.Control of a single sinusoid at 120 Hz in the time domain: (a) predicted with least mean squares and sliding mode least mean squares (Key:

Figure 6 .
Figure 6.Control of a single sinusoid at 120 Hz in the time domain: (a) predicted with least mean squares and sliding mode least mean squares (Key: , original signal; , filter output); and (b) predicted estimation error (Key:, with least mean squares; , with sliding mode least mean squares).

Figure 6 .
Figure 6.Control of a single sinusoid at 120 Hz in the time domain: (a) predicted with least mean squares and sliding mode least mean squares (Key: , original signal; , filter output); and (b) predicted estimation error (Key:, with least mean squares; , with sliding mode least mean squares).

Figure 6 .
Figure 6.Control of a single sinusoid at 120 Hz in the time domain: (a) predicted with least mean squares and sliding mode least mean squares (Key: , original signal; , filter output); and (b) predicted estimation error (Key:, with least mean squares; , with sliding mode least mean squares).

Figure 6 .
Figure 6.Control of a single sinusoid at 120 Hz in the time domain: (a) predicted with least mean squares and sliding mode least mean squares (Key: , original signal; , filter output); and (b) predicted estimation error (Key:, with least mean squares; , with sliding mode least mean squares).

Figure 7 .
Figure 7.Comparison of least mean squares and sliding mode least mean squares: (a) change in filter coefficients (Key: , w0; , w1; , w2); and (b) comparison of trajectories in the plane of filter coefficients (Key:, least mean squares; , sliding mode least mean squares).

Figure 7 .Figure 7 .
Figure 7.Comparison of least mean squares and sliding mode least mean squares: (a) change in filter coefficients (Key:

Figure 7 .
Figure 7.Comparison of least mean squares and sliding mode least mean squares: (a) change in filter coefficients (Key: , w0; , w1; , w2); and (b) comparison of trajectories in the plane of filter coefficients (Key:, least mean squares; , sliding mode least mean squares).

Figure 7 .
Figure 7.Comparison of least mean squares and sliding mode least mean squares: (a) change in filter coefficients (Key: , w0; , w1; , w2); and (b) comparison of trajectories in the plane of filter coefficients (Key:, least mean squares; , sliding mode least mean squares).
; and (b) comparison of trajectories in the plane of filter coefficients (Key:

Figure 6 .
Figure 6.Control of a single sinusoid at 120 Hz in the time domain: (a) predicted with least mean squares and sliding mode least mean squares (Key: , original signal; , filter output); and (b) predicted estimation error (Key:, with least mean squares; , with sliding mode least mean squares).

Figure 6 .
Figure 6.Control of a single sinusoid at 120 Hz in the time domain: (a) predicted with least mean squares and sliding mode least mean squares (Key: , original signal; , filter output); and (b) predicted estimation error (Key:, with least mean squares; , with sliding mode least mean squares).

Figure 7 .
Figure 7.Comparison of least mean squares and sliding mode least mean squares: (a) change in filter coefficients (Key: , w0; , w1; , w2); and (b) comparison of trajectories in the plane of filter coefficients (Key:, least mean squares; , sliding mode least mean squares).

Figure 11 .
Figure 11.Frequency response function of beam structure: (a) two different configurations; and (b) comparison of frequency response function.

Figure 12 .
Figure 12.Frequency response function of the structural model to be identified.

Figure 11 .
Figure 11.Frequency response function of beam structure: (a) two different configurations; and (b) comparison of frequency response function.

Figure 11 .
Figure 11.Frequency response function of beam structure: (a) two different configurations; and (b) comparison of frequency response function.

Figure 12 .
Figure 12.Frequency response function of the structural model to be identified.

Figure 12 .
Figure 12.Frequency response function of the structural model to be identified.

Figure 14 .
Figure 14.Frequency-domain identification of linear systems.The time-domain measurements are transformed into the frequency domain, and the calculated complex input and output amplitudes at selected angular frequencies , 1, , i i G 
a rational form that is extended by a delay term

Figure 14 .Y
Figure 14.Frequency-domain identification of linear systems.The time-domain measurements are transformed into the frequency domain, and the calculated complex input and output amplitudes at selected angular frequencies , 1, , i i G 

Figure 14 .
Figure 14.Frequency-domain identification of linear systems.

Figure 16 .
Figure 16.Notch filter control result: (a) frequency response functions; and (b) changed root locus.

Figure 16 .
Figure 16.Notch filter control result: (a) frequency response functions; and (b) changed root locus.

Figure 17 .
Figure 17.Notch filter control result: (a) tracking error of impulse response; and (b) zoomed plot of (a).

Figure 18 .
Figure 18.Notch filter control result: (a) tracking error of step response; and (b) zoomed plot of (a).

Figure 19 .
Figure 19.Experimental result for the baseline structure: (a) time domain; and (b) frequency domain.

Figure 19 .
Figure 19.Experimental result for the baseline structure: (a) time domain; and (b) frequency domain.

Figure 20 .
Figure 20.Experimental result for the modified structure: (a) time domain; and (b) frequency domain.

Table 1 .
Material properties of the piezoelectric ceramic patches (PZT) patch.

Table 1 .
Material properties of the piezoelectric ceramic patches (PZT) patch.

Table 1 .
Material properties of the piezoelectric ceramic patches (PZT) patch.

Table 1 .
Material properties of the piezoelectric ceramic patches (PZT) patch.

Table 1 .
Material properties of the piezoelectric ceramic patches (PZT) patch.

Table 1 .
Material properties of the piezoelectric ceramic patches (PZT) patch.

Table 2 .
Natural frequency and damping ratio of each peak.

Table 2 .
Natural frequency and damping ratio of each peak.

Table 2 .
Natural frequency and damping ratio of each peak.

Table 4 .
Vibrational amplitudes of each case for the original beam structure.

Table 4 .
Vibrational amplitudes of each case for the original beam structure.

Table 5 .
Vibrational amplitudes of each case for the original beam structure.