Closed Form Solutions for Nonlinear Oscillators Under Discontinuous and Impulsive Periodic Excitations

: Periodic responses of linear and nonlinear systems under discontinuous and impulsive excitations are analyzed with non-smooth temporal transformations incorporating temporal symmetries of periodic processes. The related analytical manipulations are illustrated on a strongly nonlinear oscillator whose free vibrations admit an exact description in terms of elementary functions. As a result, closed form analytical solutions for the non-autonomous strongly nonlinear case are obtained. Conditions of existence for such solutions are represented as a family of period-amplitude curves. The family is represented by different couples of solutions associated with different numbers of vibration half cycles between any two consecutive pulses. Poincaré sections showed that the oscillator can respond quite chaotically when shifting from the period-amplitude curves.

Impact loads applied to physical systems are often modeled with Dirac's generalized functions (distributions) by ignoring the short-term load profiles and using their integral effects instead.However, the presence of distributions in a differential equation of motion requires the entire equation to be considered within the theory of distributions that leads to theoretical complications especially in nonlinear cases [1,2].In the present work, we illustrate non-smooth temporal substitutions incorporating the pulse effects into a periodic time folding argument.Substituting such an argument into the differential equation of motion generates singularities of derivatives that can be effectively used to balance the loading Dirac's delta impulses or discontinuities from the differential equations.In this introductory section, we outline necessary analytical manipulations based on general time symmetries of periodic processes regardless their temporal mode shapes.Namely, any periodic process admits time shifts with reflections.Such properties can be easily incorporated into the corresponding differential equations at the preliminary stage of study by means of the following standard pair of non-smooth periodic functions As follows from the graphical illustration in Figure 1, the couple of functions τ(t) and e(t) can be interpreted as a coordinate and velocity of a small particle oscillating between two perfectly stiff barriers without energy loss.since ( ) e t is a rectangular cosine-wave with step-wise discontinuities, its first generalized derivative is given by the periodic sequence of Dirac's δ -functions as [ ] ) Distribution (2) associates with the impulsive force applied from the barriers to the particle whenever the particle strikes the barriers.In addition, the following differential and algebraic properties hold Although functions (1) are quite specific, they create a general basis for considering a broader class of oscillating processes.This is based on the statement that any periodic process ( ) x t of the period 4 T a = admits representation in the form of 'hyperbolic number' [3,4] Here, the terms ( ) X τ and ( ) Y τ depict the sin-and cos-type symmetries, and can be viewed as "real" and "imaginary" parts of the hyperbolic number due to the relationship 2 1 e = (3).Note that the term 'hyperbolic number' represents one of the multiple versions typically applied to the abstract algebraic combination [5], In particular, if a periodic function ( ) x t is even with respect to the quarter of period t a = , then ( ) 0 Y τ ≡ and identity (4) has only the "real component" ( ) x X τ = .The corresponding geometrical interpretation is shown in Figure 2. Note that the derivative .τ(t) remains undefined at discontinuity points t : τ(t) = ±1 .Also, since e(t) is a rectangular cosine-wave with step-wise discontinuities, its first generalized derivative is given by the periodic sequence of Dirac's δfunctions as .
Distribution ( 2) associates with the impulsive force applied from the barriers to the particle whenever the particle strikes the barriers.In addition, the following differential and algebraic properties hold Although functions (1) are quite specific, they create a general basis for considering a broader class of oscillating processes.This is based on the statement that any periodic process x(t) of the period T = 4a admits representation in the form of 'hyperbolic number' [3,4] x = X(τ) + Y(τ)e where τ = τ(t/a), e = e(t/a), and Here, the terms X(τ) and Y(τ) depict the sin-and cos-type symmetries, and can be viewed as "real" and "imaginary" parts of the hyperbolic number due to the relationship e 2 = 1 (3).Note that the term 'hyperbolic number' represents one of the multiple versions typically applied to the abstract algebraic combination [5], α + βj ( j 2 = 1).In our case, the unipotent j is represented by the periodic step-wise discontinuous function e(t/a) = τ (t/a) describing the velocity of impact oscillator.Identity (4) can be easily proved by substituting (5) in ( 4) and considering the result on one period −a ≤ t ≤ 3a; see Appendix A. The following example reveals temporal symmetries of the both components in (4) In particular, if a periodic function x(t) is even with respect to the quarter of period t = a, then Y(τ) ≡ 0 and identity (4) has only the "real component" x = X(τ).The corresponding geometrical interpretation is shown in Figure 2.
Geometrical interpretation of the particular case with a sine-wave temporal symmetry: observing the coordinate x does not reveal which of the two temporal variables, τ or t, is 'in play'.
Relationship (6)  Let us outline some algebraic and differential properties of representation (4).
• Algebraic properties.Taking into account that e 2 = 1, e 3 = e, e 4 = 1, . . ., gives, for instance, where Relationship ( 6) is proved by setting e = ±1 and then solving the two equations for R and I.

•
Differential properties.Despite of the presence of non-smooth and discontinuous functions on the right-hand side of identity (4), its class of smoothness coincides with a class of smoothness of the function x(t).Namely, the behavior of the right-hand side of (4) in the neighborhoods of singularities Λ = t : τ(t/a) = ±1 is determined by the functions X(τ) and Y(τ).For instance, if However, if the function x(t) has inherent step-wise discontinuities at points Λ then quantities Y(±1) 0 will represent the steps of such discontinuities.Now assuming x(t) ∈ C(R) and taking the formal time derivative of (4) gives where e (t/a) is a generalized derivative with δ-type singularities at Λ to be eliminated by imposing condition Therefore, due to condition (8), derivative (7) still preserves the structure of hyperbolic number.If x(t) ∈ C n (R) then the procedure ( 7)-( 8) can be reiterated n times by sequentially imposing boundary conditions on the functions X(τ), Y(τ) and their derivatives at τ = ±1.
Further, we use the identity (4) as a representation for solutions of differential equations by considering X(τ) and Y(τ) as a couple of new unknown functions of the temporal argument τ ∈ [−1, 1].In order to illustrate the corresponding manipulations, let us consider the equation where the input function is periodic step-wise discontinuous rectangle cosine of the period T = 4. Let us seek continuous periodic solutions of Equation ( 9) of the same period in the form of representation (4).Setting a = 1, substituting (4) in ( 9), and taking into account condition (8) gives It is seen that the left-hand side of Equation ( 10) is an element of 'hyperbolic algebra'.Setting its both components separately to zero and including conditions (8) gives the autonomous boundary value problem with no discontinuities, which can be treated in a classical way, Namely, eliminating X from the second equation gives the second-order linear differential equation Then satisfying the boundary conditions, determining X = −Y /λ, and substituting both X and Y in (4) gives periodic solution of Equation ( 9) in the final form Although Equation ( 9) is linear and can be solved by other means, such as Laplace transforms, Fourier expansions, or matching different solutions between the impulses at pulse times, solution (12) has the closed form and represents an element of the hyperbolic algebra.These properties can essentially facilitate further algebraic, differential, and integral manipulations whenever such types of solutions serve as a basis of different asymptotic procedures [4].

Nonlinear Oscillators under Periodic Impulses
Now let us consider an oscillator with a general form of the restoring force characteristic under the periodic impulsive excitation ..
where the upper dot indicates time derivative, the impulses are assumed to be two-directional, and the restoring force function, f (x, . x, t), to be continuous with respect to all the arguments.We seek a steady state continuous periodic solution of Equation (13) in the form of representation (4).Assuming that condition (8) holds and taking first derivative of (4) gives Then, substituting ( 4) and ( 14) in (13) and using the algebraic rules, as outlined in the previous section, one obtains Equation ( 15) gives two differential equations for the components X and Y with the boundary condition for X .Then the resulting boundary value problem takes the form Equations ( 17) do not include any singularities but appear to be coupled.However, under certain symmetry conditions, the number of equations can be reduced by one.For instance, one can consider the boundary value problem only for X provided that I f ≡ 0 as Y ≡ 0. The latter condition always holds, for instance in the conservative case, when the restoring force characteristic depends only upon the coordinate, f = f (x), and thus (17) is reduced to Let us assume now that the restoring force characteristic has the form Figure 3a shows dependence (20) whereas Figure 3b gives its possible physical interpretation, where the effect of combined springs and amplitude limiters is approximated with the smooth function (20) within the interval |x|< π/2 .The reason for using function (20) as a fit is that the free conservative oscillator ..
x + tan x + tan 3 x = 0 is known to have the exact analytical solution in terms of elementary functions [6,7] x = arcsin[sin A sin(t/ cos A)].Note that, despite of the formal similarity of this function with the triangle wave τ(t) (1), transition to the non-smooth limit A → π/2 is impossible due to the singularity caused by the term 1/ cos A. Typically, the presence of external loads invalidates solutions obtained for the corresponding free oscillators.However, in the present case of the impulsive load, the form of solution obviously holds in between any two consecutive pulses.From the physical standpoint, the impulsive load can practically occur in the situation illustrated by Figure 3b.The moving platform periodically strikes the perfectly stiff walls creating the impulsive inertia load on the oscillator inside the moving noninertial frame.A one parameter family of solutions can be represented as where A is an arbitrary constant, which is sufficient for satisfying both equations in (19) due to the oddness of solution (21), X(−τ) ≡ −X(τ).
( ) arcsin[sin sin( / cos )] where A is an arbitrary constant, which is sufficient for satisfying both equations in  19) and conducting analytical manipulations with trigonometric functions gives four sets of solutions for the quarter of period.Such solutions can be combined as and As follows from expression (21), |A|< π/2 and X(τ) ∼ A sin(aτ) as A → 0 .Therefore A can be viewed as a parameter characterizing the amplitude of and thus restricted to be positive.
Substituting (21) in ( 19) and conducting analytical manipulations with trigonometric functions gives four sets of solutions for the quarter of period.Such solutions can be combined as and where T = 4a is the period of response, and α = cos A. Note that solution (21) represents a one parameter set of particular solutions imposing the initial conditions on the original variables as x(0) = 0 and .
x(0) > 0. In ( 22) and ( 23), such type of the initial conditions associates with even k, whereas the odd numbers k reveal another subset of the periodic solutions with the negative initial velocities.Both of the subsets are covered by the following modification of (21) It can be shown that solutions ( 22) and ( 23) exist in the interval Figure 4 illustrates the sequence of branches of solutions ( 22) and (23) at different numbers k under the fixed parameter of the strength of pulses, F. The diagram reveals such combinations of the period and amplitude at which the oscillator has a periodic response with the loading period, T = 4a.The upper and lower branches of each loop correspond to positive and negative signs in expression (22), respectively.In particular, relationship (23) corresponds to the number k = 0 and has only the upper branch.For the selected value of F = 2.0, the minimal amplitude is A min = 1.10715, corresponding to the left boundary of the amplitude-period loops in Figure 4.A slight shift to the right from this boundary leads to bifurcations of solutions as seen from the time histories in Figure 5.The series of graphs illustrates only the first three couples of new solutions from the infinite set of solutions.Further evolution of the temporal shapes due to the amplitude increase is shown in Figure 6.One can observe quite a significant difference in both periods and temporal shapes of vibrations comparing the left and right fragments of Figure 6.This is explained by the high sensitivity of solutions near the left boundary with respect to the amplitude variation as follows from Figure 4. Finally, when the amplitude becomes close to its maximum A max = π/2, the restoring force of the oscillator itself generates quite strong pulses compared to the external loading as seen from the temporal mode shapes in Figure 7.As a result, the effect of external pulses is effectively vanishing.The oscillator' response can be viewed rather as a free vibration whose amplitude and period are dictated by the loading period., corresponding to the left boundary of the amplitude-period loops in Figure 4.A slight shift to the right from this boundary leads to bifurcations of solutions as seen from the time histories in Figure 5.The series of graphs illustrates only the first three couples of new solutions from the infinite set of solutions.Further evolution of the temporal shapes due to the amplitude increase is shown in Figure 6.One can observe quite a significant difference in both periods and temporal shapes of vibrations comparing the left and right fragments of Figure 6.This is explained by the high sensitivity of solutions near the left boundary with respect to the amplitude variation as follows from Figure 4. Finally, when the amplitude becomes close to its maximum max = / 2 A π , the restoring force of the oscillator itself generates quite strong pulses compared to the external loading as seen from the temporal mode shapes in Figure 7.As a result, the effect of external pulses is effectively vanishing.The oscillator' response can be viewed rather as a free vibration whose amplitude and period are dictated by the loading period.
(a) (b)   Analyzing the link between Figures 4 and 5 shows that solutions corresponding to different loops of the period-amplitude dependence differ by the number of half cycles between any two consecutive pulses applied to the oscillator.These graphs were reproduced by solving the original equation numerically with the Mathematica ® NDSolve package using the WhenEvent option for conditioning the velocity steps ±2F when passing the impulses under the initial conditions imposed by solution (24); see Appendix B for the code., corresponding to the left boundary of the amplitude-period loops in Figure 4.A slight shift to the right from this boundary leads to bifurcations of solutions as seen from the time histories in Figure 5.The series of graphs illustrates only the first three couples of new solutions from the infinite set of solutions.Further evolution of the temporal shapes due to the amplitude increase is shown in Figure 6.One can observe quite a significant difference in both periods and temporal shapes of vibrations comparing the left and right fragments of Figure 6.This is explained by the high sensitivity of solutions near the left boundary with respect to the amplitude variation as follows from Figure 4. Finally, when the amplitude becomes close to its maximum max = / 2 A π , the restoring force of the oscillator itself generates quite strong pulses compared to the external loading as seen from the temporal mode shapes in Figure 7.As a result, the effect of external pulses is effectively vanishing.The oscillator' response can be viewed rather as a free vibration whose amplitude and period are dictated by the loading period.
(a) (b) Analyzing the link between Figures 4 and 5 shows that solutions corresponding to different loops of the period-amplitude dependence differ by the number of half cycles between any two consecutive pulses applied to the oscillator.These graphs were reproduced by solving the original equation numerically with the Mathematica ® NDSolve package using the WhenEvent option for conditioning the velocity steps 2F ± when passing the impulses under the initial conditions imposed by solution (24); see Appendix B for the code.Note that, away from the curves of Figure 4, the response may appear to be quite complicated.This is confirmed by the Poincaré sections in Figures 8 and 9 below.In particular, Figure 8 shows the Poincaré sections obtained for the case of short loading periods under a set of randomly chosen initial conditions.Then Figure 8a-d represents a weakly damped case, in which the energy inflow is sufficient for keeping the oscillator amplitudes large enough for a "resonance" interaction between the oscillator and the impulsive load.With the increase of dissipation, the initial energies of the oscillator dissipate leading to the decrease of its own frequencies.As a result, the high frequency impulsive loading becomes effectively quasi dynamic as reflected by the drastic change in Poincaré sections; see Figure 8e,f   Note that, away from the curves of Figure 4, the response may appear to be quite complicated.This is confirmed by the Poincaré sections in Figures 8 and 9 below.In particular, Figure 8 shows the Poincaré sections obtained for the case of short loading periods under a set of randomly chosen initial conditions.Then Figure 8a-d represents a weakly damped case, in which the energy inflow is sufficient for keeping the oscillator amplitudes large enough for a "resonance" interaction between the oscillator and the impulsive load.With the increase of dissipation, the initial energies of the oscillator dissipate leading to the decrease of its own frequencies.As a result, the high frequency impulsive loading becomes effectively quasi dynamic as reflected by the drastic change in Poincaré sections; see Figure 8e,f.Note that, away from the curves of Figure 4, the response may appear to be quite complicated.This is confirmed by the Poincaré sections in Figures 8 and 9 below.In particular, Figure 8 shows the Poincaré sections obtained for the case of short loading periods under a set of randomly chosen initial conditions.Then Figure 8a-d represents a weakly damped case, in which the energy inflow is sufficient for keeping the oscillator amplitudes large enough for a "resonance" interaction between the oscillator and the impulsive load.With the increase of dissipation, the initial energies of the oscillator dissipate leading to the decrease of its own frequencies.As a result, the high frequency impulsive loading becomes effectively quasi dynamic as reflected by the drastic change in Poincaré sections; see Figure 8e,f  Finally, Figure 9 illustrates the case of a long-period impulsive loading under a relatively large damping, compared to the cases of Figure 8a-d.The difference between this case and the case represented by Figure 8e,f is that the distance between pulses is now long enough compared to the period of linearized oscillator whose period is the longest possible.As a result, the quasi dynamic branching in the Poincaré maps disappears.(e) (f) Finally, Figure 9 illustrates the case of a long-period impulsive loading under a relatively large damping, compared to the cases of Figure 8a-d.The difference between this case and the case represented by Figure 8e,f is that the distance between pulses is now long enough compared to the period of linearized oscillator whose period is the longest possible.As a result, the quasi dynamic branching in the Poincaré maps disappears.Finally, Figure 9 illustrates the case of a long-period impulsive loading under a relatively large damping, compared to the cases of Figure 8a-d.The difference between this case and the case represented by Figure 8e,f is that the distance between pulses is now long enough compared to the period of linearized oscillator whose period is the longest possible.As a result, the quasi dynamic branching in the Poincaré maps disappears.

Conclusions
Non-smooth temporal transformations are adapted for linear and nonlinear dynamical systems under discontinuous and impulsive periodic excitations.In particular, closed form analytical solutions for the strongly nonlinear exactly solvable oscillator are obtained.The solutions are represented by a family of period-amplitude dependencies with the corresponding time histories.Nonstationary responses under the presence of damping are illustrated by Poincaré sections at different magnitudes and periods of impulses.Practical importance of such a strongly nonlinear oscillator is due to the fact that the discontinuities of the restoring force characteristic can play the role of amplitude limiters without typical nonholonomic constraint conditions.

Figure 3 . 0 A
Figure 3. (a) Possible smooth fitting curve (20) approximating the restoring force of combined springs including the amplitude limiters, and (b) moving platform periodically strikes the perfectly stiff walls creating the impulsive inertia load on the oscillator in the moving noninertial frame; it is assumed that both spatial and temporal variables are appropriately scaled.As follows from expression (21), | | / 2 A π < and ( ) sin( ) X A a τ τ 

where = 4 Figure 3 .
Figure 3. (a) Possible smooth fitting curve (20) approximating the restoring force of combined springs including the amplitude limiters, and (b) moving platform periodically strikes the perfectly stiff walls creating the impulsive inertia load on the oscillator in the moving noninertial frame; it is assumed that both spatial and temporal variables are appropriately scaled.
The upper and lower branches of each loop correspond to positive and negative signs in expression (22), respectively.In particular, relationship (23) corresponds to the number = 0 k and has only the upper branch.For the selected value of 2.0 F = , the minimal amplitude is min = 1.10715A
The upper and lower branches of each loop correspond to positive and negative signs in expression (22), respectively.In particular, relationship (23) corresponds to the number = 0 k and has only the upper branch.For the selected value of 2.0 F = , the minimal amplitude is min = 1.10715A

Figure 5 .
Figure 5.Time histories of the oscillator' response near the minimum amplitude parameter A = 1.12: k = 1 (a,b); k = 2 (c,d); k = 3 (e,f); the left and right columns correspond to the lower and upper branches of loops shown in Figure 4.

Figure 5 .
Figure 5.Time histories of the oscillator' response near the minimum amplitude parameter A = 1.12: k = 1 (a,b); k = 2 (c,d); k = 3 (e,f); the left and right columns correspond to the lower and upper branches of loops shown in Figure 4.

Figure 7 .Figure 6 .Figure 7 .Figure 8 .
Figure 7. Case k = 1, larger amplitude parameter A = 1.4708:(a) corresponds to the lower part of the loop, and (b) corresponds the upper part of the loop shown in Figure 4 for k = 1.
In our case, the unipotent j is represented by