Utilizing the time delayed PPF controller to suppress vibrations of a nonlinear system containing real power exponents in damping and restoring forces

The time delayed Positive Position Feedback (PPF) controller is utilized to suppress the primary resonance of vibrations of an excited base oscillator by real power exponents of the restoring and damping forces. Multiple scales method is conducted to get the frequency response equations. The stability of the system is studied by using the Lyapunov first method. The influences of system parameters and time delay on the system response are investigated to avoid the jump phenomenon for better system performance. Time margin is deduced for most possible values of controller gain. Analytic results are verified by numerical integration of the original system equations.


Introduction
The restoring force related to the classical model of linear springs usually obeys Hook's law i.e., the force needed to deform the spring is proportional to its deflection.Practically, the restoring force is a nonlinear function of the deflection because it decreases or increases more rapidly than the deflection.Usually, the function that relates the restoring force to the deflection is a polynomial with a linear term and another one with integer powers higher than unity (mainly quadratic and cubic).Many structural materials of the springs can cause this type of nonlinearity.Also, the origin of this type of nonlinearity may be geometrical or physical configurations, as in the case of a pendulum performing small vibrations, beams, cables, electrical circuits.This type of restoring force was widely studied in many scientific contributions, as was reported [1].However, many materials whose elastic deformations cannot be exactly, described as a polynomial [2,3].In these materials, the force-deflection relationship includes non-integer exponents of a fractional order [2][3][4][5][6][7][8][9].Specifically, the nonlinear materials of Ramberg-Osgood type, which have been applied to depict the elastic properties of the aircraft materials such as aluminum and titanium [4].While, The Ramberg-Osgood equation was created to describe the nonlinear relationship between stress and strain in some metals.
Vibration isolators having nonlinearity in both stiffness and damping terms were analyzed under harmonic excitation [5] using the method of harmonic balance.Oscillators with a non-negative real-power restoring force and quadratic damping are considered [6].The expressions for the energy displacement function as well as the closed form exact solutions for the relationship between subsequent amplitudes were derived in terms of incomplete Gamma functions.The relative and absolute displacement transmissibility of a base excited vibration isolation system were analyzed, and the performance characteristics were expressed in terms of the damping parameters [7].Kovacic [8] studied the steady-state response of forced damped nonlinear oscillators.In the considered oscillators, the restoring force has a nonlinear term in the form of a non-integer exponent with a fractional order and a linear term.The damping term was assumed in a power form.The method of multiple scales was improved to cover the cases in which the nonlinearity is not necessarily small.Huang et al. [9] utilized the time-delayed cubic velocity feedback to improve the performance of a vibration isolation system whose restoring and damping forces are in a real-power form.They also studied the primary resonance, dynamic stability and transmissibility for this system under base excitation using the method of multiple scales.In the last decades, active control techniques have been widely used in mechanical and structural vibrations for sensing structural vibrations and delivering appropriate control forces to counteract structural motions.However, active control techniques may lead control systems to instability due to presence of time delays.The reason is that the time delays are unavoidable during the process of sensing and calculating the necessary control forces.Researchers paid a great attention to the stability analysis of delayed controlled systems [10][11][12][13][14][15][16][17][18][19].Abdelhafez and coauthor applied positive position feedback (PPF) controller to reduce the vibrations of a forced and self-excited nonlinear beam [19].They deduced that the time margin of the system depends on the sum of all time delays.To our knowledge, previous works didn't apply PPF controller on systems with real-power exponents of the restoring and damping forces.
In this work, the PPF controller is utilized to suppress the primary resonance vibrations of the base excited oscillator with real-power exponents in the restoring and damping forces given [9].This oscillator is harmonically excited on its base.The method of multiple scales is applied to get the possible resonance conditions and the amplitude-phase modulating equations.Also, the frequency response equations and the stability conditions are derived.The effects of some effective parameters on the dynamical behavior of the system will be explored.The optimal values of system parameters for avoiding the jump-phenomenon are deduced.Effects of time delays on the system response are studied extensively.It is known that, time margin is the amount of time delay that system can bear without being unstable.Time margin of the system is deduced for most possible values of controller gain and the optimal conditions for system operation are derived.Some analytical results are verified by numerical integration of the original closed-loop system equations.Finally, main results are summarized.

Mathematical model
The system considered here comprises of a mass  isolated from a harmonically excited base through a linear spring, a nonlinear spring composed of nonlinear materials "Ramberg-Osgood type", and a nonlinear damper similar to that reported [7].The equation of motion of the base is  = cos(Ω ) and the equation of motion of this primary system is given in [9] as follows: where  =  −  is the relative displacement of the mass  with respect to the vibrating base, dots indicate derivative with respect to time ,  is an arbitrary positive real number, and the coefficient  symbolizes the corresponding damping coefficient. is the coefficient of the linear stiffness. is an arbitrary non-negative real number, and the coefficient  stands for the corresponding stiffness coefficient.sgn() and sgn( ) are the sign functions.Introducing the dimensionless parameters  = /,  =  , Ω = Ω / , 2 =   / ,  =  / ,  =   / and  = / yields the following dimensionless equation of motion of the system:  () + 2 sgn  ()  ()| +  () + sgn () ()| = cos(Ω). ( Utilizing the positive position feedback controller to isolate the mass  from the base vibrations yields the following system model: where dots denote derivatives with respect to ,  is the controller mode,  and  are the controller gains.The block diagram of this closed loop system is presented in Fig. 1.

Mathematical analysis
Applying the method of multiple time scales [20] yields a first-order approximate solution to Eqs. ( 3) and ( 4) by seeking the solution in the form: where  is a small dimensionless perturbation parameter,  = ,  =  are the two time scales that represent the fast and slow times.The time derivatives are expressed in terms of  and  as follows: For the necessity of applying the multiple time-scales method, we may scale system parameters as follows: Substituting Eqs. ( 5) to (8) into Eqs.( 3) and ( 4) and equating coefficients of the same power of  give the following set of differential equations: The general solutions of Eqs. ( 9) and ( 10) can be expressed as follows: where ( ),  ( ), ( ) and  ( ) are functions of the slow time scale  .| | sgn( ) and |  | sgn(  ) can be expanded by using Fourier series as follows: |  | sgn(  ) = | ( )| (− sin(  +  ( )) +  sin(3  + 3 ( ))), (16) where coefficients  ,  ,  and  are given in the Appendix.Also,  and  can be defined as follows: ( ,  ) =  ( )cos( ( −  ) +  ( )), (17)  ( ,  ) =  ( )cos( ( −  ) +  ( )). ( Expanding  and  by using Taylor series yields: where, prime denotes derivative with respect to  .Substituting Eqs. ( 13) to (20) into Eqs.( 11) and ( 12), we get: where  stands for complex conjugate.In this work, the primary resonance case is treated in presence of internal resonance between the system and the controller, so we have: where  and  are detuning parameters used to describe the closeness of the considered resonance cases.After inserting Eq. (23) into Eqs.( 21) and ( 22), the conditions for a bounded solution of Eqs. ( 21) and ( 22) are: where  =   −  =   −  ,  =   −  +  =   −  +  .Separating real and imaginary parts of Eqs. ( 24) and (25) and rescaling parameters to their original form yield the following system of autonomous differential equations: (26)

Equilibrium solution and stability analysis
At steady-state oscillations, the derivatives of system states are zeros ( =  =  =  = 0).Inserting this condition into Eq.( 26), we can get the following closed form equations: =    . ( We have utilized Lyapunov first method to study the stability of the equilibrium solution.So, we need to linearize Eq. ( 26) about the equilibrium solution  ,  ,  , and  by using Taylor expansion to first-order.Each system state consists of its value at equilibrium plus a small perturbation as follows:  =  +  ,  =  +  ,  =  +  ,  =  +  ,  =  ,  =  ,  =  ,  =  , where  ,  ,  , and  are perturbations which are small with respect to  ,  ,  ,  .Linearizing Eq. ( 26) about the equilibrium solution as illustrated in Eq. (29) yields the following linear system of ordinary differential equations: where coefficients  ; ,  = 1, 2, 3, 4 are given in the appendix.The solution is stable, if and only if, all eigenvalues of the Jacobian matrix of Eq. (30) have negative real parts.
Fig. 2 presents a comparison between the FRC of the system controlled by positive position feedback (PPF) controller versus the uncontrolled system.This figure shows that the PPF Controller is very effective in reducing vibrations of the primary system under control.The solid small circles in the figure represent a numerical validation of the analytical results.Also, from the figure it could be observed that, the numerical solutions coincide with the analytical solution when | | tends to zero so that, the accuracy of the analytical solution increases when the system approaches the studied resonance case.However, both of the analytical and numerical solutions are generally of same trend.Fig. 3(a) presents the frequency response curve (FRC) of the system and the controller under different values of .From those figures it can be observed that, increasing  leads to an increase in the peak amplitudes as indicated by the dashed arrow in them.Also, the minimum amplitude " " at  = 0 increases slightly when  increases.The behaviors of the peak amplitudes  and the minimum amplitude  are indicated in Fig. 3(b).When the exponent of damping force increases, the minimum amplitude of the FRC of the system increases monotonically until the value  = 1, the left peak amplitude increases monotonically until  = 4, and the right peak amplitude increases monotonically until  ≃ 9. So, the FRC of the system nearly saturates after  = 9 as seen in Figs.3(a  From Eq. (31) we obtain four solutions of  as funntions of .These four solutions are combined to plot the FRC.The maximum peaks are obtained by the equation  () = 0, with substitution for all parameters as constants except .Substituting  () = 0 in Eq. (31) yields the minimum amplitudes.In this way can study the effect of any parameter on the peaks of the amplitudes of the system FRC.
From Fig. 3 we deduced that, increasing  does not cause multi-valued solution regions or frequency-island phenomenon if PPF controller is implemented.However, variation of the restoring force exponent  can lead to multi-valued solution regions in the system response.Fig. 3(a) depicts the possible cases of the system behavior for any value of .
In Fig. 4(a), it is clear that there are three regions and each of these regions has a possible different case of system behavior.Fig. 4(b) presents the system behavior at  = 0.05 as a sample of the region R1.In this region, the system response is distorted which means that there is no real solution for the frequency response equations when  → 0. Fig. 4(c) depicts the FRC of the system when  belongs to the region R2.In region R2, jump phenomenon occurs as shown by arrows.The FRC of the system at  = 1.5 as a sample of region R3 is presented in Fig. 4(d).For  > 0.28, there is no jump and there is a good vibration suppression bandwidth.So, R3 is the best α region for system operation when controlled by the PPF controller.The equivalent natural frequency of the system can be estimated approximately from the following equation:  =  +   .
(32) Fig. 5 shows the plots of the equivalent natural frequency  as a function of  and the system amplitude .It can be observed that  increases largely for small values of  or  as indicated in the figure.This is the reason behind the bad results given by PPF controller in Fig. 4(b) (region R3).The effect of  on the system performance was also tested and gave nearly the same effect of  .The system amplitude at  = 0 decreases as  increases as seen in Fig. 7(c).The analytical solution is validated numerically using small circles as shown in the figure .The time margin "the amount of time delay that the system can bear without being unstable" of the system decreases as  increases.The relation between the time margin of the system and controller gain  is explained later on.Eq. ( 28) shows that the system behavior depends on the quantity   +   , but  = 1 and the PPF controller is perfectly tuned i.e.,  = 1.So, the system response and its stability are dependent on the quantity  =  +  .Hence,  represents the overall time-delay of the system.The following data in Table 1 of the overall time-delay  against the controller gain  can be deduced using stability criteria.The values of  in Table 1 represent the maximum overall delay that the system can bear without being unstable at the considered values of  .Plotting the data in Table 1 yields Fig. 8(a).From Fig. 8(a) it could be observed that  "Time margin" is a decreasing function of  .Utilizing fitting methods yields the following formula that links  and  , such that: (0.0007cos(0.063 ) + 4.686sin(0.063)) Eq. ( 33) is used to get the time margin of the system corresponding to any value of  .When  increases the time margin of the system decreases.The stable solutions region of the system for all possible values of  and  is indicated by the yellow region in Fig. 8(a).Any point on the curve in Fig. 7(a) represents the time margin of the system at the corresponding value of  .Other values of  and  were used to obtain Fig. 8(a) and it is observed that, the same results are nearly obtained when those values increases, but the time margin and the area of stable solutions region decrease when  increases.Fig. 8(b) describes the behavior of the system eigenvalues at  = 0.428 when  changes from 0 to 0.87.A pair of eigenvalues move towards the imaginary axis when  increases until intersecting it at  = 0.5.So, Hopf-bifurcation occurs at the point P1 in Fig. 8(a), the system vibrates by a bounded quasi-periodic motion at this point as seen in Fig. 8(e), and the system is critical stable.All points on the curve of Fig. 8(a) have the same response as the point P1.At the point P1,  = 0.428 and  = 0.5, however  =  +  so  = 0.428 and  = 0.5.
The red line in Fig. 8(c) is drawing at  +  = 0.428.So,  = 0.428 at any point on the line plotted in Fig. 8(c).Then this line symbolizes all possible combinations of  and  that represent the time margin of the system at  = 0.5.The stable solution regions indicated in Fig. 8(a) and 8(c) are the regions at which all eigenvalues have negative real part.The points P2, P3, and P4 in Fig. 8(c) are studied in Figs.8(d) and 8(f) to explain the possible different responses of the system.The time response depicted in Fig. 8(d) shows that the system vibrates by a stable periodic motion when  belongs to the stable solutions region.At  equals the time margin, the system response becomes critical stable and vibrates by a bounded quasi-periodic motion as shown in Fig. 8(e).The overall delay  must not exceed the time margin to avoid occurrence of system instability indicated in Fig. 8(f).Fig. 9(a) depicts the system amplitude as a function of the external detuning parameter  and the magnitude of the external force .From this figure it can be seen that, the peak amplitudes and the minimum amplitude of the system are increasing when the value of  is increased.Also, the vibration suppression bandwidth is not affected largely by increasing .The force response curves of the system before and after control are presented and validated numerically in Fig. 9(b) at  = 0.The system amplitudes after controller are very small compared to them before control as seen in Fig. 9

Conclusions
Utilizing PPF controller Guarantees wider vibration suppression bandwidth around  = 0.When the exponent of damping force  increases, the smallest amplitude of the FRC of the system increases monotonically until  = 1, the left peak amplitude increases monotonically until  = 4, and the right peak amplitude increases monotonically until  ≃ 9. So, the frequency response curve of the system nearly saturates after  = 9.The PPF controller reduces vibrations of the system effectively around  = 0 for all possible values of  as seen in Fig. 3. Increasing  does not cause multi-valued solution regions or frequency-island phenomenon.Variation of the restoring force exponent  yields three possible cases of system response as shown in Fig. 4.
The vibration suppression bandwidth increases as  and/or  increases.The minimum amplitude of the system decreases when  and/or  increases.The effects of time-delays on the system stability were studied with the deduction that, the time margin of the system depends on the sum  +  .When the overall delay of the system is small than the time margin, the system vibrates by a stable periodic motion.When the overall delay of the system equals the time margin, the system response becomes critical stable, the system vibrates by a bounded quasi-periodic, and Hopf-bifurcation occurs.The system becomes unstable when the overall delay  exceeds the time margin.Increasing , within the stable solutions region, boosts the right peak amplitude, suppress the left peak amplitude and has no effect on the vibration suppression bandwidth.The system vibrations amplitudes after controller are very small if compared to their values before control.

Fig. 3 .Fig. 4 .
a) FRC of the system, b) FRC of the controller, c) Maximum and minimum values of the system amplitude  against , d) System amplitude as a function of the detuning parameter  and  a) Jump and no-jump regions in the system response under different values , b) FRC of the system at  = 0.05, c) FRC of the system at  = 0.24, d) FRC of the system at  = 1.5

Fig. 5 .
Fig. 5. Equivalent natural frequency  as a function of  and the system amplitude  Fig. 6(a) gives us the key features of the FRC of the system when  belongs to region R3.If  4 then the right peak amplitude is larger than the left peak amplitude.The right peak amplitude becomes equal to the left peak for the values of  > 4. The minimum amplitude of the system increases monotonically as  increases until  ≃ 1 and becomes constant for  > 1.Hence, Fig. 6(b) verifies the general description presented in Fig. 6(a).

Fig. 8 .
a) The overall delay  against controller gain  at  = 0, b) Imaginary part of eigenvalues versus their real parts when  = 0.428, c) The stable solutions region w.r.t. and  at the point P1 in a), d) Time response of the system at the point P2 in c), e) Time response of the system at the point P3 in c), f) Time response of the system at the point P4 in (b)