Abstract

Nonviscously damped structural system has been raised in many engineering fields, in which the damping forces depend on the past time history of velocities via convolution integrals over some kernel functions. This paper investigates stochastic stability of coupled viscoelastic system with nonviscously damping driven by white noise through moment Lyapunov exponents and Lyapunov exponents. Using the coordinate transformation, the coupled Itô stochastic differential equations of the norm of the response and angles process are obtained. Then the problem of the moment Lyapunov exponent is transformed to the eigenvalue problem, and then the second-perturbation method is used to derive the moment Lyapunov exponent of coupled stochastic system. Lyapunov exponent also can be obtained according to the relationship between moment Lyapunov exponent and Lyapunov exponent. Finally, the effects of various physical quantities of stochastic coupled system on the stochastic stability are discussed in detail. These results are validated by the direct Monte Carlo simulation technique.

1. Introduction

Viscoelastic materials are widely used in aerospace, construction, textile, and other industries, because they have a series of excellent properties, including light weight, high strength, wide source, and good shock absorption. These materials’ stress depends on the past time history of stain; that is, the stress will increase correspondingly with the increase of time and the bucking will be more likely to occur. Therefore, the research of dynamical behavior of viscoelastic system has received a lot of interests in recent years [14].

To better understand the dynamical behavior of the viscoelastic system, many methods had been put forward. McTavish et al. [5] presented GHM method to analyze the linear multiple degree of freedom viscoelastic damping systems. Li et al. [6] studied the decay rate of energy functional of nonlinear viscoelastic full Marguerre-von Kármán shallow shell system. Later, Li [7] also showed the existence and uniqueness of weak solution of this system by applying the Galerkin finite element method. Nieto and Ahmad [8] used the generalized quasi-linearization method to obtain the explicit approximate solutions of an initial and terminal value problem for the forced Duffing equation with nonviscous damping. Recently, Li and Du [9] studied general energy decay of a degenerate viscoelastic Petrovsky-type plate equation with boundary feedback through a priori estimate and analysis of Lyapunov-like functional.

The dynamic stability as an important field of nonlinear dynamics has attached increasing attention in recent years. Guo et al. [10] derived a new simple sufficient condition of the global asymptotic stability of the integrodifferential systems with delay through constructing suitable Lyapunov functionals combined with the analytical technique. And then Meng et al. [11, 12] obtained a sufficient condition of uniform stability for nonlinear integrodifferential system. However, when we consider the dynamical behavior of the system in a real environment, the stochastic excitation cannot be ruled out, such as wind gusts, earthquakes, and ocean waves [1315]. Stochastic stability of the system under the stochastic excitation has attracted more and more attention [1620]. The Lyapunov exponent specially has been researched by many scholars as an important indication for judging the stability of stochastic system. For example, Potapov [21] derived the sufficient condition of the almost-sure asymptotic stability of elastic and viscoelastic systems by the Lyapunov’s direct method. He [22] also studied the stability of elastic and viscoelastic systems under the non-Gaussian excitation. For stochastic stability of high dimensional coupled system, Pavlović et al. [23] proposed the direct Lyapunov method to investigate the almost-sure stochastic stability of a viscoelastic double-beam system under parametric excitation. Then they used the same method to analyze the instability of coupled nanobeam systems subject to compressive axial loading [24]. These works only indicate that the system is almost-surely stable.

According to the theory of large deviation, which was first proposed by Baxendale and Stroock [25], the almost-sure stability of stochastic system does not mean the moment stability. Therefore, it is necessary to study the moment Lyapunov exponent of a stochastic system. The moment Lyapunov exponent can judge not only almost-sure stability but also the moment stability [26]. Nevertheless, analytic solutions of the moment Lyapunov exponent are very difficult to derive. To overcome this, many scholars developed different approximation methods to study the moment Lyapunov exponents of the stochastic system. For example, Sri Namachchivaya and Van Roessel [27] studied the moment Lyapunov exponents of two-degrees-of-freedom coupled elastic oscillators under real noise excitation by combining an asymptotic approximation for the moment Lyapunov exponents. Kozić and his associates [2830] developed the first-order perturbation approach to obtain weak noise expansion of moment Lyapunov exponents and Lyapunov exponents for a stochastically coupled double-beam system and Timoshenko beam system, respectively. Subsequently, Stojanović and Petković [31] used a perturbation method to study the moment Lyapunov exponents and the Lyapunov exponents of the three elastically connected Euler beams. More recently, Deng [32] applied the averaging method to establish the moment Lyapunov exponent of coupled gyroscopic stochastic system under bounded noise excitation. But it is rare to see studies of the moment Lyapunov exponent of a viscoelastic coupled system with nonviscous damping. Deng [33] investigated stochastic stability of two-degrees-of-freedom coupled viscoelastic system under white noise through moment Lyapunov exponent, but its damping term is viscous.

The purpose of this paper is to study stochastic stability of a viscoelastic coupled system with nonviscous damping subject to Gaussian white noise excitation through moment Lyapunov exponents and Lyapunov exponents, in which the nonviscous damping term is expressed by Boltzmanns superposition integral with a hereditary relaxation kernel. The paper is organized as follows. Section 2 derives the governing equations of motion of a stochastic viscoelastic system with a nonviscous damping. And then solving the problem of the moment Lyapunov exponent is changed into the eigenvalue and eigenfunction problem through a suitable transformation. In Section 3, the zeroth-order and the first-order and second-order perturbation solutions of the moment Lyapunov exponents are derived, respectively, based on the second-order perturbation method. The effects of different physical quantities on the stochastic stability of the coupled system are discussed under the given parameters. Correspondingly, the results are verified by means of the direct Monte Carlo simulation in Section 4. Finally, conclusions will be drawn in Section 5.

2. Problem Formulation

Considering a coupled viscoelastic system with nonviscous damped structure driven by white noise excitation, the governing equations can be expressed aswhere and are generalized displacements. and are natural frequencies. is a small parameter. and are nonviscous damping coefficients. and are constants. is a Gaussian white noise with zero mean and noise intensity . is a relaxation function of nonviscous damping, which can be expressed by in which is the relaxation parameter. Note that the limit case , the nonviscous damping force reduces to classical viscous damping.

Let , , , ; system (1) can be represented aswhere

Using a coordinate transformation,where represents the norm of the response. and are the angles process of the coupled stochastic system. denotes the coupling or exchange of energy between the coupled system. Equation (3) can be rewritten asTherefore, the solution forms of can be expressed as from (6), where is quickly varying process relative to the process due to . Therefore, the following approximations will be obtained:Substituting (7) into the nonviscous damping, one obtainsSimilarly, the following results also can be derived:whereSubstituting (8) and (9) into (6) and letting , (6) can be written as the following Itô stochastic differential equations:where is the standard Wiener process; the drift and diffusion coefficients can be rewritten as

To obtain the eigenvalue problem for the th moment Lyapunov exponents of a four-dimensional linear Itô stochastic system, a linear stochastic transformation [34] is introducedwhere is a new norm process depending on the transformation function . The transformation function is defined on the independent stationary phase processes , , and in the ranges , , . Applying Itô formula, one can obtain

If the transformation function is bounded and nonsingular, both processes and have the same stability behavior. Therefore, transformation function is chosen so that the drift term of the Itô differential equation (14) does not depend on the phase processes , , and , so one can obtain

Comparing (14) and (15), such transformation function can be written by the following equation:where and are the following first-order and second-order differential operators:in which

In the following investigation, the perturbation theory will be used to obtain the solution of (16).

3. Moment Lyapunov Exponents

The method of deriving the eigenvalue problem for the moment Lyapunov exponents of a two-dimensional linear Itô stochastic system was first applied by Wedig [34]. The eigenvalue problem for a differential operator of three independent variables , , and will be identified from (16), in which is the eigenvalue and is the associated eigenfunction. Meanwhile, the eigenvalue is seen to be the Lyapunov exponent of the th moment of system (1) from (15). Applying the perturbation theory, both the moment Lyapunov exponent and the eigenfunction are expanded in power series of ; that is,

Substituting (19) into (16) and equating terms of the equal powers of lead to the following equations:where must be positive and periodic in the range , .

3.1. The Zeroth-Order Perturbation

Substituting (17) into (20), the zeroth-order perturbation equation isThrough , thus, , will be obtained. Then (21) does not contain , and the eigenvalue is independent of . Therefore, will lead to for arbitrary . At this point, (21) will be reduced toUsing the method of separation of variables, so we can express aswhere , , and are functions about , , and , respectively.

Substituting (23) into (22), (22) will become Solving the above differential equation (24), one can easily obtain In addition, since function is periodic about and , the following boundary conditions can be obtained:Substituting (23) and (25) into (26), we can determine the constants . Therefore, is a function about only to be determined; that is, .

The adjoint equation of (22) is Applying the method of separation of variables to solve (27), one obtainswhere is an arbitrary function.

3.2. The First-Order Perturbation

From (20), the first-order perturbation equation isBased on the analysis results in Section 3.1, one obtainsThen the solvability condition of (30) isApplying the sense of the arbitrary of function , (31) will becomewhere then (32) is simplified as the following equation:Then the boundary conditions of (34) are determined by considering the adjoint equation for the case of ,where is the Fokker-Planck operator and and are the entrance boundaries [27]. The eigenfunction satisfies zero Neumann boundary condition, and is the largest eigenvalue of (34) with zeros Neumann boundary. Therefore, the solution of (34) can be calculated from an orthogonal expansion [34]. Under zeros Neumann boundary conditions, can be expressed by a Fourier cosine series [27]; that is,Substituting (36) into (34), multiplying by in both sides and integrating for , we can obtainwhere In order to guarantee the existence of the nontrivial solution for each , the coefficient matrix must equal zero. Therefore, the problem of calculating is translated into calculating the leading eigenvalue of . The sequence of approximations will be constructed by the eigenvalues of a sequence of the following submatrices:

Obviously, the set of approximate eigenvalues obtained by this method converges to the associated true eigenvalues as . In general, we can obtain the approximate eigenvalues by truncating , which can be written as and

The comparison results of approximate analytical moment Lyapunov exponents obtained by truncating and the direct Monte Carlo simulation results are shown in Figure 1. It is easily found that the approximate results agree well with the simulation results. Therefore, the first-order perturbation will be reduced toWithout loss of generality, there is a relationship between two frequencies of the form , in which and are integers. The second frequency can be rewritten as . Then the function in (41) can be rewritten aswhere function is the periodic function on and given asNote that the general solution of (41) can not be obtained except in some special cases. Considering the nature of the coefficients of (41), the series expansion of function can be presented in the following form:

Substituting (44) into (41) and equating terms of the equal trigonometry function to give a set of partial differential equations, one obtains where the function can be written asin which are arbitrary functions of two variables, and we make assumptions as follows:Here, the unknown constants , ,    will be determined by using the following conditions:That is,

3.3. The Second-Order Perturbation

To further investigate the effects of parameter on stochastic stability of system, the second-order perturbation equation of (20) will be used; that is,Based on the solvability condition, one obtainsThe solution has the following form:where the values , , are given in Appendix A.

Substituting (40) and (52) into (19), the second-order approximate solution of the moment Lyapunov exponent is obtained asThe corresponding Lyapunov exponents can be expressed aswhere the values , , and are given in Appendix B.

4. Stochastic Stability Analysis

Through the above-mentioned analysis, the analytic expressions of moment Lyapunov exponents and Lyapunov exponents have been derived. Next, we use two indices to go on stochastic stability analysis. The th moment Lyapunov exponent is defined aswhere denotes the expected value and is a solution process of a random system. If , as , and it means that the th moment is stable. Lyapunov exponent is the slope of at the origin; that is,If , it indicates that the system is almost-surely stable.

The th moment and almost-sure stability boundary can be determined from and , respectively. By comparing the second-order perturbation with the first-order perturbation, we find that the moment Lyapunov exponents and the Lyapunov exponents under the second-order perturbation contain all the parameter information of the coupled stochastic system. Therefore, we only consider the second-order perturbation in the subsequent research.

The boundaries of the th moment and the almost-sure stability with respect to nonviscous damping parameters and are shown in Figure 2. We find that the moment stability boundaries are more conservative than the almost-sure stability boundaries. This means that the almost-sure stability cannot assure the moment stability. And the boundaries of stability for th moment would increase with the increase of the order .

The effects of relaxation parameter on moment Lyapunov exponents are shown in Figure 3. It describes that the stability index and the root of will decrease and the system becomes more unstable. Correspondingly, the analytical results and the direct Monte Carlo results are shown in Figure 3(b) for the given . And Lyapunov exponent is negative; it means the system is almost-surely stable.

The moment Lyapunov exponents for different nonviscous damping parameter are discussed in Figure 4. It is clear that the bigger values of can enhance the system's stability for . Correspondingly, the analytical results and the direct Monte Carlo results are shown in Figure 4(b) for the given different . It illustrates that the system also is almost-surely stable. However, when becomes sufficiently large, , the system becomes more unstable.

Next, the effects of on the moment Lyapunov exponents are studied in Figure 5. One easily finds the system is moment stability when , and the stability of the system increases with the increase of from 0.01 to 0.1. Afterward, the moment Lyapunov exponents are changed from negative to positive; the system is varied from moment stability to moment instability with the increase of . Numerical results are verified by the direct Monte Carlo in Figure 5(b).

Then, in order to discuss the influences of the noise intensity on stochastic stability of coupled stochastic system, the moment Lyapunov exponents are shown in Figure 6 with change of values of . They describe that the coupled stochastic system could be almost-surely stable. Meanwhile, the moment stability of the coupled system would reduce with the increase of the noise intensity .

The variations of the moment Lyapunov exponents for different are also analyzed. From Figure 7, we find that the coupled stochastic system may be almost-surely stable due to the value of slope less than 0 at the origin. However, the moment stability of system is not assured for or is sufficiently large due to the moment Lyapunov exponents greater than 0. Besides, one can easily see that the moment stability of the system enhances with the increase of .

Finally, the effects of on the moment Lyapunov exponents are investigated and shown in Figure 8. From those figures, the system may be almost-surely stable in terms of the Lyapunov exponents. It is also indicated that the system has the characteristic of the moment stability. And the moment stability of the coupled system will increase with the increase of under given parameters. The approximate analytical results agree well with numerical simulation results by the direct Monte Carlo.

5. Conclusions

In this paper, the moment stability and almost-sure stability of coupled viscoelastic system with nonviscous damping subject to Gaussian white noise excitation are investigated. The nonviscously damped structure is assumed to follow the exponential integral relation, which is possible to model the rate of energy dissipation. Motion equations are transformed into the coupled Itô stochastic differential equations by means of the method of stochastic averaging. Then the linear transformation is introduced to derive the eigenvalue problem. In order to obtain the analytical expressions of eigenvalue, the second-order perturbation method is used. Finally, the effects of Gaussian noise, the relaxation parameter, the nonviscous damping coefficient, and two physical quantities, that is , , on the moment stability are discussed in detail. Then the corresponding results are verified through the direct Monte Carlo simulation technique. Based on the results above, the appropriate parameters will help to strengthen the stability of the coupled viscoelastic system with nonviscous damping.

Appendix

A. Values of , , and

B. Values of , , and

Data Availability

All the numerical calculated data used to support the findings of this study can be obtained by calculating the equations in the paper, and the white noise is generated randomly. The codes used in this paper are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work is supported by the National Nature Science Foundation of China (no. 11402139, no. 11572247) and the Outstanding Innovation Programs of Graduate Students in Shanxi (no. 2017BY121). The writers are grateful to the stochastic dynamics team of the Northwestern Polytechnical University for discussions and valuable suggestions during the course of this study.