Abstract

A procedure for analyzing stationary responses of lightly nonlinear vibroimpact system with inelastic impact subjected to external Poisson white noise excitation is proposed. First, the original vibroimpact system is transformed to a new system without velocity jump in terms of the Zhuravlev nonsmooth coordinate transformation and the Dirac delta function. Second, the averaged generalized Fokker-Planck-Kolmogorov (FPK) equation for transformed system under parametric excitation of Poisson white noise is derived by stochastic averaging method. Third, the averaged generalized FPK equation is solved by using the perturbation technique and inverse transformation of the Zhuravlev nonsmooth coordinate transformation to obtain the approximately stationary solutions for response probability density functions of original vibroimpact system. Last, analytical and numerical results for two typical lightly nonlinear vibroimpact systems are presented to assess the effectiveness of the proposed method. It is found that they are in good agreement and the proposed method is quite effective.

1. Introduction

Vibroimpact system, as a class of typical nonsmooth system [1], has attracted much attention in the past few decades due to its existence in various engineering applications [2, 3]. Certain applications or problems involved impacts include vibratory pile drivers, heat exchanger tube wear in nuclear power stations, large gaps in labyrinth seals for rotating shafts of high-pressure pumps, and a moored body with an inextensible mooring line [47]. There are two commonly used methods to model impacts. One is the classical or hard impact model [1], and the other is the Hertz contact model or soft impact model [8, 9]. The classical impact model consists of two parts: the ordinary motion described by a differential equation and the impact condition. In the Hertz contact model, the contact force between two elastic bodies is assumed to be proportional to the 3/2 (sometimes with different values) power of the relative displacement between them. By means of the two models, many interesting dynamical behaviors were observed such as chaos motion and bifurcations [10, 11], Hopf bifurcation [12], grazing bifurcation [13, 14], torus bifurcation [15], chattering bifurcation [16, 17], sticking motion [18, 19], and border collision bifurcation [20].

As is known to all, random excitations widely exist in vibroimpact systems. Their effects usually can not be neglected and sometimes are even quite large [21]. Therefore, studies on the prediction of stochastic responses of vibroimpact systems are naturally important but difficult due to the presence of random excitations and nonsmooth factors. In recent years, some methods [2227] have been developed to explore the response probability density function (PDF) which is an important characteristic of a stochastic system. For the classical impact model, Dimentberg et al. [28, 29] studied stochastic response of linear vibroimpact system under Gaussian white noise. Feng et al. [30, 31] explored the response PDFs of Duffing-Van der Pol vibroimpact system under independent and correlated Gaussian white noises, respectively. Zhu [32] obtained the response PDFs of a vibroimpact Duffing system with a randomly varying damping term. For the Hertz contact model, Huang et al. [33] investigated stationary responses of multi-degree-of-freedom vibroimpact system under white noise excitations. Xu et al. [34] analyzed random vibration problems of vibroimpact system with inelastic impact. Although stochastic responses of vibroimpact system have been reported, the above studies are limited to responses of vibroimpact systems subject to Gaussian white noise. In practice, many random excitations are essentially non-Gaussian, for example, those due to sea wave [35] and wind [36]. The Poisson white noise is the most used model for non-Gaussian random excitations [3740]. However, the research on stochastic responses of vibroimpact system subject to Poisson white noise has been rarely addressed [41] and is far from enough. This paper is devoted to presenting a procedure to predict the stochastic responses of lightly nonlinear vibroimpact system subject to external Poisson white noise excitation. In the following parts, the Zhuravlev nonsmooth coordinate transformation [22], the stochastic averaging method [42], and the perturbation technique [43] are applied in succession to obtain the stationary solutions of response PDFs. Two examples are also presented to assess the effectiveness of the procedure. Comparison between results obtained from the proposed method and those from Monte-Carlo simulation shows that they are in good agreement. Therefore, the proposed procedure is quite effective.

2. Problem Statement

Without loss of generality, a simple mass-spring-damper model (Figure 1) for single-degree-of-freedom lightly nonlinear vibroimpact system with inelastic impact subjected to external Poisson white noise excitation is considered, and the corresponding equation of motion can be expressed aswhere is a small parameter; denotes lightly linear or nonlinear damping force; is the stiffness coefficient; is Poisson white noise; and denote the rebound and impact velocity, respectively; the subscripts “+” and “−” refer to values of velocity just before and after impact; is the restitution factor whose value is among zero and one. Particularly, when equals one, the system reduces to an elastic impact one. denotes that the rigid barrier is located at the system’s equilibrium position. Equation (1a) governs the system’s motion between impacts as . Equation (1b) describes the impact law, which means that the oscillator’s velocity has an instantaneous change when it impacts with the barrier at .

The Poisson white noise can be treated as the formal derivative of the compound Poisson process ; that is,in whichwhere is the Poisson counting process with a constant impulse arrival rate which also means the expected number of event occurrences per second. is a random variable representing the random magnitude of th impulse arriving at time . is a unit step function. Assume that are independent identically distributed random variables with zero mean which are independent of the impulse arrival time . is a process with independent increments, and the th correlation function of whose increment process is as follows [44]:

By the relation between the correlation function and the moment function of random processes, omitting the small quantity that higher than , the kth moment function can be derived as

Given that the nonsmooth characteristic of vibroimpact system described by (1a) and (1b), the following coordinate transformation introduced by Zhuravlev [22] is used asThe relation between the original variables , and the new variables , can be simply represented by Figures 2 and 3.

Substituting (6) into (1a) and (1b), the transformation maps the domain of original phase plane onto the whole phase plane (see Figure 3). Since when , this implies also that . When , the impact/rebound condition for the transformed state variables is found to be according to Dimentberg and Menyailov [28]. Note that the identity , and the original equation of motion can be replaced bywhere represents the time instant of impact corresponding to , that is, Also, the velocity jump becomes proportional to instead of for original one. Using the Dirac delta function, an additional impulsive term can be introduced as follows:According to [28], it is reasonable to treat (7b) as an additional impulse of (7a) as Therefore, (7a), (7b), and (8) can be combined into one new equation asThe above transformed equation may permit rigorous analytical study since there contains no velocity jump and the impact condition is taken place by an additional impulsive term.

3. Stochastic Responses

3.1. Stochastic Averaging

When the values of approaches to zero, the energy loss of impact is very small, and the system becomes a quasiconservative one as long as the damping coefficients and random excitations are weak. Therefore, the stochastic averaging method can be applied to study the response of the transformed system analytically. For simplicity, is assumed of order . The transformed equation (9) is equivalent to a pair of first-order equations:whereThe total energy of system isSubstituting (12) into (10) yields the following Stratonovich stochastic differential equations about displacement and energy with parametric excitation of Poisson white noise:The corresponding Itô stochastic differential equations are derived using the converting rule for systems with Poisson white noise excitation proposed by Di Paola and Falsone [37]:Give thatEquation (14) can be simplified asThe generalized FPK equation associated with (16) iswhere .

Particularly, if the probability distribution of is symmetric, all odd order moments are zero. Thus, we haveThe conditional probability density given by Stratonovich [45] isBy means of the formulation and integrating (18) with respect to , the following averaged generalized FPK equation for probability density is obtained:whereIt is worth noting that the impulsive damping term should be averaged over a half period since there have two impacts in each period (see Figure 3).

3.2. Approximately Stationary Solution

The averaged generalized FPK (20) is too complicated to be solved analytically. Herein, we seek the approximately stationary solutions. Then, the left-hand side of (20) vanishes and the original equation can be rewritten aswhereIt is easy to know from (23) that is a constant. In another way, the probability current function satisfies . Thus, we only need to solve the following equation to obtain the approximately stationary solutions:The above equation can be solved by a perturbation technique [43]. The perturbation solution is assumed to be of the formwhere satisfies the boundary condition:Substituting (26) into (25) yields the following equations:Generally speaking, the approximately stationary solution of can be obtained analytically and works well with the numerical simulation results up to order of . The higher terms can be neglected since their effects are relatively small.

The joint PDF for the variables and can be further obtained asUsing the inverse transformation of formula (6), the joint PDF for original displacement and velocity is calculated asThus, the marginal PDFs for displacement and velocity are easy to be derived, respectively:

4. Illustrative Examples

In this section, two examples are shown to illustrate the proposed procedure. In order to compare the analytical results and numerical results, the impulses magnitudes of Poisson white noise are assumed to be Gaussian distributed with zero mean. Thus, an interesting result can be used to simplify the practical computation.

4.1. Rayleigh Oscillator

In this example, a typical Rayleigh vibroimpact system under external Poisson white noise excitation is considered:In this case, . The following transformed equation without velocity jump is obtained by means of the Zhuravlev transformation and Dirac delta function.The corresponding Itô stochastic differential equation is derived asAccording to (20) and (21), is calculated as follows in the averaged generalized FPK equation for probability density :Using the perturbation technique, the approximately stationary solution for energy PDF is obtained:whereFurthermore, the joint PDF for original displacement and velocity can be calculated by substituting (37) and (38) into (29) and (30) as

To verify the effectiveness of the proposed procedure, analytical results and Monte-Carlo simulation results are shown in Figure 4 for PDFs of system energy with , , , and and two different values of small parameter ( and ). The solid line is the analytical results obtained by the proposed method while • represents the Monte-Carlo simulation results based on the original system (33a) and (33b). It is obvious that the two results agree well with each other. In order to illustrate the effects of average arrival rate of external Poisson white noise on the accuracy of the approximately stationary solution, Figure 5 shows PDFs of system energy analytically and numerically with different values of ( and ) when is fixed as 0.1. It is seen that the proposed method works well.

To investigate the effects of restitution coefficient on the stationary response of vibroimpact system, numerical and analytical results are shown in Figures 68 with the parameters , , , and and different values of . It is seen from Figure 6 that the stationary PDFs of total energy are shown for = 0.97, 0.98, 0.99, and . The PDFs of total energy present lower peaks with the increase of restitution coefficient. It is found that the values of peak are at the system’s static equilibrium position when = 0.97, 0.98 while they do not stay there when = 0.99, 0.995. This may mean that the steady state has a qualitative change with the increase of restitution coefficient. To show it more clearly, Figures 7 and 8 provide analytical and numerical results of joint PDFs of the displacement and velocity for and , respectively. It is seen that the analytical results agree well with the numerical results. In Figure 7, the joint PDFs of the displacement and velocity have a singular peak which appears at the static equilibrium position. When increases to , the joint PDFs present the shape of crater in Figure 8. This illuminates the occurrence of P-bifurcation according to the concept of stochastic bifurcation. In other words, the restitution coefficient can induce the occurrence of P-bifurcation.

4.2. Van der Pol Oscillator

As an example, a typical Van der Pol vibroimpact system under external Poisson white noise excitation is considered:Based on (9), the transformed system without velocity is obtained asThe Itô stochastic differential equation about displacement and energy is derived asIn the averaged generalized FPK (20), here is calculated asThus, is obtained analytically by solving the corresponding averaged generalized FPK equation aswhereFurthermore, the PDFs , , and can be calculated following (29)–(31), respectively.

To demonstrate the effectiveness of the proposed method, numerical and analytical results of response PDFs are provided in Figures 913. The effects of damping coefficient on the accuracy of the energy PDFs are shown in Figure 9 with , , , and and two different values of ( and ). Stationary PDFs of total energy obtained numerically and analytically for different values of restitution coefficient are presented in Figure 10. It is clear from Figures 9 and 10 that the analytical results give powerful validation with Monte-Carlo simulation results. The joint PDFs of the displacement and velocity for different impulse intensities are shown numerically and analytically in Figures 11 and 12. The corresponding marginal PDFs and are provided in Figures 13(a) and 13(b), respectively. It is seen that the results obtained by using the proposed method and those from Monte-Carlo simulation are in good agreement. Additionally, the increase of impulse intensity only effects the values of peaks of stationary response PDFs, and there is no occurrence of P-bifurcation.

5. Conclusions

The paper is devoted to presenting a solution procedure for predicting stationary responses of single-degree-of-freedom lightly nonlinear vibroimpact systems with inelastic impact subjected to external Poisson white noise excitations. The classical model with instantaneous impacts is used and the rigid barrier is located at the static equilibrium of the system. The solution procedure consists of three steps. First, the Zhuravlev nonsmooth coordinate transformation and the Dirac delta function are adopted to convert the original vibroimpact system to a new system without velocity jump. Second, the Itô stochastic differential equations with correction terms are obtained and the averaged generalized Fokker-Planck-Kolmogorov (FPK) equation is derived by stochastic averaging method. Third, the perturbation method is utilized to solve the averaged generalized FPK equation and the inverse transformation of the Zhuravlev nonsmooth coordinate transformation is used to establish the PDFs for original vibroimpact system. Furthermore, two examples are presented to show the effectiveness of the proposed method. It is found that the analytical results obtained by the proposed method and the Monte-Carlo simulation results agree well. Therefore, the proposed procedure is quite promising and can be extended to more vibroimpact systems subject to Poisson white noise excitation.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 11702200, 11702214, and 11702201) and the Fundamental Research Funds for the Central Universities (Grant no. JB170705).