This paper introduces higher-order solutions of the quadratic nonlinear stochastic oscillatory equation. Solutions with different orders and different number of corrections are obtained with the WHEP technique which uses the Wiener Hermite expansion and perturbation technique. The equivalent deterministic equations are derived for each order and correction. The solution ensemble average and variance are estimated and compared for different orders, different number of corrections and different strengths of the nonlinearity. The solutions are simulated using symbolic computa tion software such as Mathematica. The comparisons between different orders and different number of corrections show the importance of higher-order and higher corrected WHEP solutions for the nonlinear stochastic differential equations.
Analysis of the response of linear and nonlinear systems subjected to random excitations is of considerable interest to the fields of mechanical and structural engineering [
According to [
The application of the WHE [4-10] aims at finding a truncated series solution to the solution process of a stochastic differential equation. The truncated series composes of two major parts; the first is the Gaussian part which consists of the first two terms, while the rest of the series constitute the non-Gaussian part. In non-linear cases, there exists always difficulties of solving the resultant set of deterministic integro-differential equations got from the applications of a set of comprehensive averages on the stochastic integro-differential equation obtained after the direct application of WHE. Many authors introduced different methods to face these obstacles. Among them, the WHEP technique [
The WHE was originally started and developed by Norbert Wiener in 1938 and 1958 [
As in [
As in [5,6], quadrate oscillation arises through many applied models in applied sciences and engineering when studying oscillatory systems [
In [
In the current work the generalized WHEP technique developed in [
This paper is organized as follows. The formulation of the quadratic nonlinear stochastic oscillatory equation is outlined in Section 2. The WHEP technique is reviewed in Section 3. The equivalent deterministic system is derived in Section 4. In Section 5, the mean and variance of the solution is simulated with different order, different number of corrections and different values of the nonlinearity strengths.
In this paper, the nonlinear oscillatory equation:
is considered under stochastic excitation and with the proper set of initial conditions which is assumed to be deterministic. The operator L is a general linear operator and in the case of the oscillatory equation it will be:
where is the undamped angular frequency of the oscillator and is the damping ratio. The nonlinearity is introduced as losses of degree strengthened by a deterministic small parameter. The uncertainty is introduced through white noise scaled by a deterministic envelope function. The white noise is considered here as a function of time but it can be generalized in time and space as it was declared in [
Theorem (1): The solution of Equation (1), if exists, is a power series in, i.e.
The theorem can be proved using the mathematical induction with the Pickard iterative technique [
The WHEP technique will be used in this work to determine the equivalent deterministic set of equations. The deterministic equations are then solved to obtain the solution kernels and hence the mean and variance of the response.
As a consequence of the completeness of the WienerHermite set [
The solution function can be expanded in terms of Wiener-Hermite functionals as [
Or after eliminating the parameters, for the sake of brevity, we get:
where and is a k-dimensional integral over the variables. The first term in the expansion (3) is the non-random part or ensemble mean of the function. The first two terms represent the normally distributed (Gaussian) part of the solution. Higher terms in the expansion depart more and more from the Gaussian form. The Gaussian approximation is usually a bad approximation for nonlinear problems, especially when high order statistics are concerned [
The components are called the (deterministic) kernels of the WHE for. The variable w is a random output of a triple probability space, where is a sample space, is a - algebra associated with and P is a probability measure. For simplicity, will be dropped later on.
The functional is the order Wiener-Hermite time-independent functional. The Wiener-H functionals form a complete set with and: the white noise. By construction, the Wiener-Hermite functions are symmetric in their arguments and are statistically orthonormal, i.e.
The average of almost all Wiener-Hermite functionals vanishes, particularly,
The expectation and variance of the solution will be:
The WHE method can be elementary used in solving stochastic differential equations by expanding the solution as well as the stochastic input processes via the WHE. The resultant equation is more complex than the original one due to being a stochastic integro-differential equation [
The WHEP technique for general nonlinear exponent, general order and general number of corrections follow the steps [
1. Truncate the expansion (3) to contain only
kernels, i.e.
and then 2. Substitute into the stochastic partial differential Equation (1);
3. Use the multinomial theorem to expand the nonlinear term in (1);
4. Multiply by and then apply the ensemble average. This will lead to equations in the kernels;
5. For each kernel, apply the perturbation technique up to corrections, i.e.
;
6. Equating the coefficients of in both sides to get equations for each kernel .
This will lead to the following equations [
where
And the expectations are computed as:
It was explained in [
Equations (5) and (6) can always be solved using the proper sequence. The first Equations (5) are solved independently to get then they are used to compute the other components in (6). For, the component is obtained by solving
with the original initial conditions which are assumed deterministic. For, the component
is obtained by solving with zero initial conditions. The other components will be zeros due to zero right hand side and zero initial conditions. Equations (6) specify the solution sequence to be followed. The component is evaluated in terms of the previously computed components
. This means that the 1st corrections for all kernels, are solved firstly then solving the 2nd corrections for all kernels, ,
up to the corrections for all kernels
.
These results are consistent with the known results obtained using WHE. In WHE, higher order kernels are driven by lower order kernels, and at the bottom, the Gaussian kernels are driven by the random forcing directly. So, the lower order kernels are usually dominant in magnitude [
The statistical properties of the solution will now be calculated as:
If, then it will be convergent if [
for. This means that should obey an upper bound condition after which divergence is obtained.
The formulation given in Equations (5) and (6) are quite general and could be used for analysis of the response of an any linear operator with degree nonlinearity and subjected to an arbitrary, stationary, or nonstationary Gaussian or non-Gaussian random excitation.
Consider the quadratic nonlinear oscillatory equation with excitation function:
With the initial conditions
In case of zero initial conditions, the exact solution can be obtained using different methods such as the theory of linear differential equations or the Laplace transform, and it will be the convolution:
where withwhich is the angular frequency of the underdamped harmonic oscillator.
For, the solution will results in:
The solution (11) of the model Equation (10) can be used as a model solution that is used in all kernels after considering the proper right hand side for the kernel equation.
Applying the above mentioned WHEP algorithm to get the following systems of equations of the quadratic (n = 2) nonlinear oscillatory equation and first order (m = 1) Gaussian approximation and different number of corrections. The initial conditions are assumed deterministic and hence only the zero-order and zerocorrection kernel equation will has the initial conditions. Other kernels equations will have zero initial conditions.
:
: The above equations in addition to:
: The above equations in addition to:
: The above equations in addition to:
: The above equations in addition to:
In case of zero initial conditions and zero deterministic excitation [i.e.], we shall have:
Which means the all of and are become zeros.
The second order equations will be:
:
: The above equations in addition to:
: The above equations in addition to:
: The above equations in addition to:
The third order equations are:
:
: The above equations in addition to:
: The above equations in addition to:
: The above equations in addition to:
The fourth order equations are:
:
: The above equations in addition to:
The following output is simulated using Mathematica. The solution (11) of the model Equation (10) is used to get all kernels with the proper right hand side. The mean response and the response variance are then calculated from the kernels using Equation (9):
Figures 1 and 2 show the response mean and variance, respectively, for the case of zero initial conditions, zero deterministic exciatation and unit envelope function multiplied by the white noise. The angular frequency and the damping ratio. The nonlinearity strength is changed to study its effect on the response mean and variance. As it is shown in the figures, the nonlinearity strength greatly affect the amplitudes of the mean and variance. It should not be increased after a certain value to obtain convergent solution. This value depends on the different parameters of the problem.
Also, we can notice that higher correction solutions are required with longer time intervals.
In Figures 3 and 4, the envelope function is taken as. This will attenuate the effect of the white noise as the time increases. We can notice the attenuation effect on the response variance as the time increases. In this case, the variance vanishes with time and the solution becomes nearly deterministic.
In Figures 5 and 6, non-zero initial condition is considered; and. The third correction mean and variance differs from the first and second corrections noticeably. This ensures the need for higher corrected WHEP solutions especially with larger values for the nonlinearity strength.
In Figures 7 and 8, the second order solution is simulated for different nonlinearity strengths with the first and second corrections. Higher corrected solutions will be time consuming even with multi-core machines. Mathematica automatically runs in parallel when multiple cores are available. The estimation of the response variance is more difficult than the response mean for higher corrections.
In Figures 9 and 10, the third order solution is simulated. Also, higher corrected solutions will be time con-
suming especially when estimating the response variance.
The fourth order first correction solution is also obtained but it is the same as the third and the second order with first corrections. Higher corrections will be more and more expensive. There is a need for an alternative method such as the numerical estimations to overcome the difficulties in using symbolic packages.
In the present paper, we investigate the mean response of the quadratic nonlinear oscillatory system subjected to nonstationary random excitation using WHEP technique. The equivalent deterministic equations have been derived up to fourth order. The corrections are considered up to the fifth correction for lower orders and second corrections for higher orders. The mean and variance are simulated using Mathematica. It is observed that the WEHP technique can be applied to study the non-Gaussian response of any random systems. Moreover, the higherorder terms of the WHE should be adopted to obtain the responses of nonlinear random systems, and the integrodifferential equations for the kernels should be calculated. There is a need for numerical estimations when higher order and higher corrections are required.
We should thank Prof. Magdy El-Tawil (has passed away) for motivating us to cooperate and do this work.