Vibration Annihilation of Sandwiched Beam with MROF DTSMC

In the present paper, an analytical model of a flexible beam fixed at an end with embedded shear sensors and actuators is developed. The smart cantilever beam model is evolved using a piezoelectric sandwich beam element, which accommodates sensor and actuator embedded at distinct locations and a regular sandwiched beam element, having rigid foam at the core. A FE model of a piezoelectric sandwich beam is evolved using laminated beam theory in MATLAB®. Each layer behaves as a Timoshenko beam and the cross-section of the beam remains plane and rotates about the neutral axis of the beam, but it does not remain normal to the deformed longitudinal axis. Keeping the sensor and actuator location fixed in a MIMO system, state space models of the smart cantilever beam is obtained. The proper selection of control strategy is very crucial in order to obtain the better control. In this paper a DSM controller designed to control the first three modes of vibration of the smart cantilever beam and their performances are represented on the basis of control signal input, sensor output and sliding functions. It is found that DSM controller provides superior control than other conventional controllers and also MROF DSM controller is much better than SISO DSM controller.


Introduction
Smart structures [1] are systems having particular functions viz.sensing, processing, actuation and making them suitable for structural health conditioning, vibration suppression of structures.Piezoelectric materials are found most suitable to be used as active components in smart structures [2].The apposite-V.Rathi, A. A. Khan DOI: 10.4236/eng.2017.99046 756 Engineering ness of piezoelectric materials as sensors and actuators has gained the focus in health monitoring of structures like beams, plates, and shells [3]- [15].Krommer [16], Rao and Sunar [17] have shown the implementation of piezoelectric materials as both for sensing and actuation.Active control through bonded piezo components was studied by Moita et al. [18].An optimal linear quadratic generator control strategy to control the structures is advised by Ulrich et al.Young et al. [20] presented a finite element simulation of flexible structures with output feedback controllers.Aldraihem et al. [21] developed the model of the laminated beam based on EBT and TBT.Abramovich [22] has obtained an analytical formulation and closed form results of the laminated beam based on TBT with piezoelectric sensors and actuators.Chandrashekhara and Vardarajan [23] acquainted a finite element model of the laminated beam to evolve deflection in beams with various end conditions.Sun and Zhang [24] have suggested the basis of shear mode to produce transverse deflection in embedded structures.Aldraihem and Khdeir [25] expounded the analytical model and exact solution of Timoshenko beam with shear and extensional piezoelectric actuators.Zhang and Sun [26] have presented an analytical model of surface mounted beam with shear piezo actuators at the core.The top and bottom layers obey EBT and core obeys TBT.Donthireddy and Chandrashekhara [27] have proposed a model with embedded piezoelectric components.Rathi and Khan [28] have modeled a smart cantilever beam with surface mounted and embedded shear sensors and actuators on the basis of TBT and justify that embedded components of flexible structures provide better control than surface mounted arrangement and also emphasized on optimal location of sensors and actuators in embedded beam.Chammas and Leondes [29] [30] have presented the pole assignment by piecewise constant output feedback for LTI systems while Werner and Furuta [31] [32] focussed on fast output sampling for LTI system.Janardhan et al. [33] designed a controller based on MROF using the samples of control input and sensor output at different sampling rates.Bandyopadhyay et al. [34] adduced a DTSM control that has the use of switching function in control results in QSMC.
A numerous types of control policies for the SISO and MIMO state space presentation of the active structures using the Multiple Rate Output Feedback (MROF) dependent Discrete Sliding Mode Control (DSMC) approach is depicted in this monograph.The key objective instigating this control technique is to constrain and damp out the flexural or transverse vibrations of active beam when they are subjected to external annoyance.The control technique used on the basis of Bartoszewicz law and does not need to use switching in control function and hence eradicate chattering.This method does not need the reconnaissance of the system states for feedback being using solely the output samples for designing the controller.The schematic espousal is more viable and may be easy to accomplish in true life applications.out using a switching function in control and has the property of finite time convergence to the QSM band.A discrete output feedback sliding mode control algorithm in [33] based Bartoszewicz's control law [35] and MROF [36] is used for the vibration suppression of flexible structures.
In the present situation, the disturbance is the external force input ( ) r t in form of impulsive force applied to the free end of the beam and hence producing the vibration.DSM controllers with multirate output feedback plan evolved and applied to the system with the plant to attenuate the vibrations earliest.The methodology is described as follows: Consider a CT SISO system sampled with an interval α seconds and given as where, α ∆Θ is the uncertainty in the state, ( ) Let the desired sliding manifold be governed by the parameter vector T p such that T 0 p α ϒ ≠ and resulting quasi-sliding motion is stable and assume that disturbance be bounded such that Which satisfies the inequality ( ) where, ζ − and 1 ζ + are lower and upper bounds on the disturbance.We take, ( ) ( ) The switching surface is given by The QS mode is the motion such that ( ) ≤ , where the positive constant η is termed as quasi-sliding mode bandwidth.A significant reduction of con- trol effort and better quality of quasi-sliding mode control is found.A reaching law advised by Bartoszewicz [35] is as follows where, ( ) The value of the positive integer * n is chosen by Engineer so as to have a compromise in between rapid convergence an amplitude of control input ( ) By controlling the decay rate ( ) * n , the convergence of ( ) , can be described as The control law satisfying the reaching law (Equation ( 7)) and get sliding mode for the system as given in Equation ( 4) can be computed to be When control input given in Equation ( 11) substituted into the system, it would sure for any * , n n > the switching function would satisfy the expression ( ) ( ) ( ) Thus, system states adjudicate within a QSM band having less than half bandwidth as given in [37].From [33] MROF based algorithm using an advanced reaching law can be attained.Let the advanced reaching law be [35] given as A new variable ( ) o n is incorporated here.The control input generated can be given by using algorithm in [33], ( ) ( ) Here, ( ) The bounds on ( ) o n are given as ( ) sen to be greater the observability index j of the system defined as for system given by triplet ( ) , , A B C is the minimum integer value of j such that Thus the control input can be estimated by using the past output samples and immediate past input.But at 0 n = , there are no past outputs for use in control, here ( ) 0 u is obtained by neglecting ( ) and 0 o (no disturbance before 0 n = to affect the system), so we have, When control input obtained from eq. ( 14) is used in system obeys reaching law ( ) ( ) ( ) ( ) When ( ) So we have, ( ) ( ) ( ) This can be written as It can be emphasized that this algorithm does not need the assessment of system states for the creation of control input.This control technique is used to design a multi-rate output feedback based DSM control to attenuate the transverse disturbance in a flexible structure which is modeled on the basis of Timoshenko beam theory for 3 vibratory modes.

Finite Element Modeling of an Embedded Beam
An embedded beam consists of three layers having a piezoelectric patch with the obdurate foam in between two thick steel beams shown in Figure 1.The lead zirconate titanate (PZT) layer acts as both actuator and sensor in thickness shear actuation mode.The foam and PZT together behave like a core element to obtain embedded beam model [28].
The presumption is that the mid layer is perfectly bonded to the rest of the structure and thickness of binder is neglected (hence preventing shear-lag, slip or layer delamination during vibration) resulting a strong blend between parent structure and piezoelectric patches.The binder used between the layers have been assumed no added mass or stiffness to sensor or actuator.Engineering For the parts having no PZT patch, the auxiliary space is being filled full with a material like obdurate foam.Again, there is a strong blend between foam and parent structure.Thus, embedded beam consists of slabs and a light weight core are effectively good in producing bending and shear action.
In analysis of embedded beam, the poling orientation of piezo patch in the axial direction.The displacement domain is based on first-order shear deformation theory (FSDT).The element is considered to have invariable elasticity modulus, moment of inertia, mass density, and length.The wiring capacitance is neglected between the sensor and signal conditioning device.The gain is assumed to be 100 for signal conditioning device.
Consider a beam having an element with two nodes.The longitudinal axis of the embedded beam element stays along x-axis and beam vibrates along x-z plane.The beam element has three degrees of freedom these are, axial displacement of the node u, transverse displacement of the node w, and bending rotation θ .An auxiliary degree of freedom in the form of sensor voltage occurs.As sen- sor voltage is invariable through the element, the number of electrical degree of freedom is one.At each node, a bending moment and a transverse shear force act.The slope of the beam ( ) possesses two parts first one is the bending and the second one is shear deformation angle ( )

Equations of Motion
The displacement of the beam is written as; Strains are; The constitutive equations of the beam element are where, xx N is in-plane force resultant in longitudinal direction, xx M is in-plane moment resultant in transverse direction and xz F is shear force resul- tant in transverse direction and they are given as Here, b is beam width, z is depth of material direction from beam reference plane, h is the height of beam and piezoelectric patch.A are extensional, bending-extensional, bending and transverse shear stiff- nesses and expounded by using lamination theory ( ) ( ) ( ) ( ) ( ) ( ) ) where, n z is the distance of th n lamina from longitudinal axis, N is the total number of laminas, κ is shear correction factor and The angle ψ is the angle between the fiber direction and x-axis of beam.
Various material constants are obtained individually for steel, PZT and foam by relations listed in appendix., , ( ) ( ) ( ) , when PZT layer is oriented along longitudinal direction, ( ) where , U K δ δ and W δ correspond to virtual strain energy, virtual kinetic energy and virtual work done by external forces respectively and are given as where, 0 q is transverse load (equals to external force applied at the free end of beam).
( ) are mass inertias of beam cross-section and are defined as ( ) For case of static loading with invariable beam properties.We have simplified form of equation of motion as For the solution of unknowns, the degree of polynomial used for axial displacement, u and bending rotation, θ must be one order lower than that used for transverse displacement, w to satisfy the compatibility.Here we used V. Rathi , q q x q x θ = + + (51) The boundary conditions are Using the boundary conditions in Equations ( 50)-( 52), the unknown coefficients , i j p q and j r are determined.Substituting the evaluated unknowns into Equations ( 50)-( 52) and arranging them into matrix form, we obtain expressions for , w u and θ in terms of shape functions and nodal displacements.
The mass matrix of beam element is given as where, [ ] I is the inertia vector and given as [ ] V. Rathi  is a symmetric local mass matrix of size 6 6 × for a beam element, its coefficients are given in the appendix.
The stiffness matrix of beam element is given as where A is the area of beam cross-section and Finally, the stiffness matrix of the beam element is given as  is a symmetric local stiffness matrix of size 6 6 × for a beam ele- ment, its coefficients are given in appendix.
The mass matrix and stiffness matrix of the regular beam are obtained by placing foam core in between two laminas of steel.Similarly, a piezoelectric patch is used in place of foam between two laminas to obtain piezoelectric beam element.

Equation of Sensing Component
Sensor works on direct piezoelectric effect, which is used to evaluate the output charge developed due to straining of the structure.The electric displacement produced by the sensor is directly proportional to strain rate.The charge ( ) appeared on piezoelectric sensor surface is given by Gauss law as ( ) where, z D is electric displacement in thickness direction and PZ A is the area of shear PZT patch.If poling is done along the thickness direction having electrodes on top and bottom surfaces, the electric displacement is given as Here,  The current induced by the sensor is where, d  is the strain rate and T b is a constant vector of ( ) The actuator input voltage is

Equation of Actuating Component
Actuator works on converse piezoelectric effect which is used to evaluate the straining effect caused due to actuator.The strain produced in PZT patch is directly proportional to electric potential applied to lamina and is given as Or we may write Here, 2 is the distance between neutral axis of beam and neutral axis of PZT patch.
The control force ctrl f produced by actuator and applied to beam is obtained from Equation (84) as ( ) Or simplified as where, g is a constant vector of size ( ) × for a double node beam element and depends on actuator location and type.The total force vector in existence of any external force vector is Engineering where, ext f representing the external disturbance vector.

Dynamic Equation of Smart Structure for A MIMO Model
The dynamic equation for the smart structure is obtained by using both the regular beam element and the piezoelectric beam element.The equation of motion of the smart structure and sensor output is The mass and stiffness matrices of the beam in the equation of system (64) can be changed by varying the location of PZT patch on beam and by varying the number of regular and piezoelectric beam elements.The generalized coordinates are introduced into Equation (64) with a transformation = d Ta , in order to reduce it further so the resulting equation showing the dynamics of desired vibratory modes.Here T is a modal matrix containing the eigenvectors showing the desired number of modes of vibration.This procedure is applied to derive the uncoupled equations governing forced vibration in terms of principal coordinates by inducing linear transformation between generalized coordinates d and principal coordinates a and hence decoupling into equations related to each mode.Using the transformation = d Ta .Equation (88), Equation (89) be- come Premultiplying Equation (90) by T T , we get Which may be written as where, T are called the generalized mass and stiffness matrices.
The generalized external force vector ( ) where, ( ) r t is external force input to beam.
V. The generalized control force vector ( ) ( ) where, ( ) r t is control force input to beam.
The structural modal damping matrix C by using Rayleigh proportional damping as where, , α β are frictional damping constants and structural damping constant in C .The dynamic equation and sensor output of smart structure is finally

State Space Formulation for A MIMO Model
Here in the present case of actively controlled cantilever beam only first three vibratory modes are controlled since more energy is stored in lower order modes as similar to lower order Fourier components are larger in magnitude and the higher frequency components are smaller as the harmonics increase in number [38].The state space model for first three vibratory modes can be obtained as, let = a x where, and Which can be further simplified as In state space form ( ) ( ) ( ) where, A is state matrix, B is input matrix, C is output matrix, D is transmission matrix, E is external load matrix coupling the disturbance to the system and all are in continuous time model of LTI system.The size of matrices , , , × with D being a null matrix.

Simulation for Controllers for Smart Beams with MIMO Using Embedded Piezo
A cantilever beam of proposed parameters as given in Table 1, the piezo material properties in Table 2 and material constants in Table 3.The beam is divided into 10 finite elements and shear piezo are embedded into parent structure as sensors and actuators as presented in Figure 1.The actuators are placed in between two thick steel beams at FE position 2 and 5, while the sensors are placed at FE location 6 and 10, hence developing a single MIMO system with 2 inputs and 2 outputs.The MIMO model achieved by using the TBT, piezoelectric coupling, FE modeling and state space approach by taking first three vibratory modes into consideration.An external impulsive force ext f of 10 N is employed for 60 ms at the free end of the cantilever beam.There are three inputs to the system, the first one is the external force ext f responsible for the disturbance.Other inputs are the control inputs ( ) to actuators by the controller.The control strategy presented in this monograph is implemented to design a multi-rate output feedback based discrete sliding mode controller to attenuate the first three modes of vibration of a cantilever beam by using smart structure approach of smart embedded beam with MIMO.
The performance of the model with multiple inputs and multiple outputs for active vibration attenuation by performing simulations in MATLAB® and analyzing different responses.The discrete sliding mode controller is designed and implemented to MIMO system.The responses viz.control inputs, sensor outputs, and switching are demonstrated in included   It is obvious from previous works that the control is more persuasive at the root with the sensor output voltage is significantly more due to the substantial dispensation of the bending moment near the firm end for the rudimentary mode, thus provoking a stupendous strain rate and the susceptibility of the sensor/actuator duo rely on its placement in the beam and the vibration attributes of the system precarious on collocation of the piezo pair and also on some other numerous facet viz. the gain of amplifier employed, the mode number and the placement of piezo patches at the nodal points from fixed end [28].Modelling a smart structure inclusive of sensor/actuator mass and stiffness and by altering its orientation in the beam from the free end to the fixed end acquaint an ample modification in the system's structural response attributes.Sensor voltage is lower when the piezo patch is imposed at the free end due to the exiguous strain rate and hence demand more control endeavor.MIMO control is superior over SISO control due to its multifarious interactions of input and output and all-inclusive control endeavor needed by MIMO controller is less than SISO controller and also placing the piezo at two distinct FE locations on the beam establishing the significant modification in the system structural traits than placing it lonely at a location.[39].

Elucidation of Results
The multirate output feedback dependent DSMC strategy are more harmonious as compared to the other control approaches viz.periodic output feedback (POF) and fast output sampling (FOS) controllers [40].The multirate output feedback based DSMC policy are more episodic as compared to the other control techniques.In discrete quasi sliding mode control (DQSMC) with output samples, there is a necessity of switching function for control and hence engendering some chattering phenomenon [41], while control strategy presented in the present article is the MROF based DSMC technique obtained from Bartoszewicz's law does not demand any use of switching function and provides control input directly in form of past control data and past samples.The system responds well in closed loop and does not manifest inexpedient chattering phenomenon.MROF based DSMC employ the signum function in the control input and the control is computed from the immediate past control value and the past control output samples.The fractious system takes an extended time to damp out the oscillations in contrast to the system with the designed sliding mode control input means without control the transient response was preeminent and with control, the vibrations are quashed.
From simulation results, it can be inferred that sensor output at FE 6 is more than sensor output at FE 10 by approximately 10 times due to its high strain rate at FE 6 as compared at FE 10 and also the control input are approximately 10 times smaller in case of MROF as compared to SISO case [41].In case of MROF Hence, it can be concluded that the multivariable control is best among all the models due to its multilevel interactions on both input and output.A MIMO model furnishes excellent energy distribution and even good administration of actuation forces and minimal requirement of control forces as compared to SISO model for the case of smart cantilever beam with embedded sensors and actuators.

1 o + and 1 o
− are the upper and lower bounds of ( ) o n .The variable ( ) o n represents the disturbance effect on sampled output

Figure 1 .
Figure 1.A three layered embedded beam (stacking sequence: steel/PZT or foam/steel) with MIMO. 55 from Equations (40)-(42) in to Equation(39), we get equation of motion for general, unsymmetric piezoelectric laminated beam as per FSDT with shear deformation and rotary inertia as, ( )

N
θ are modal shape functions due to , w u and θ which are given as shape functions in matrix form, the relations between inertial forces vector  and nodal displacement vector d as

.
element which depends on sensor type, its properties and its location in embedded structure and described as The mass and stiffness matrices of the regular beam and piezoelectric beam element are known as local mass and stiffness matrices and give only the mass and stiffness matrices of only one finite element.The mass and stiffness matrices of entire beam i.e. divided into 10 finite elements are obtained by assembling the local matrices by applying finite element technique and resulting matrices are called global mass matrix of smart structure have both sensor and actuator masses and stiffnesses.

Figure 2 .Figure 3 .
Figure 2. Plot for control input for embedded smart cantilever beam.

Figure 4 .
Figure 4. Plot for Sliding functions for embedded smart cantilever beam.
DSM controllers have designed to control first three modes of vibration of a flexible cantilever beam modelled based on Timoshenko beam theory.New active vibration control scheme to suppress the vibrations of MIMO model has developed.The actuators are located at 2 nd and 5 th FE positions while sensors are set at 6 th and 10 th FE position to form the embedded smart cantilever beam with V. Rathi, A. A. Khan DOI: 10.4236/eng.2017.99046772 Engineering 10 finite elements.The piezo crystals are located in the central core at prescribed locations, rest core filled with rigid foam, and this central core is sandwiched between two regular steel beams as shown in Figure 1.
) Mass Matrix Coefficients for Embedded Beam Element ( ) t is applied voltage to th n actuator having thickness ( ) The mass matrix for beam element is finally given as , A. A. Khan DOI: 10.4236/eng.2017.99046764 Engineering

Table 1 .
Properties of the steel cantilever Timoshenko beam.

Table 2 .
Properties of the piezoelectric shear sensor and actuator when the beam is divided into 10 finite elements.

Table 3 .
Material properties and constants.
better enduringness, brisk convergence and exalted steady state authenticity of the system.The technique used is more feasible as the output being used rather than states.
technique, the states of the system are needed neither for switching function assessment nor for the feedback denotation.DSMC algorithm are computationally V.Rathi, A. A. Khan DOI: 10.4236/eng.2017.99046773 Engineering unpretentious, ensures