0
Research Papers

Responses of Linear and Nonlinear Oscillators to Fractional Gaussian Noise With Hurst Index Between 1/2 and 1

[+] Author and Article Information
Mao Lin Deng

Department of Mechanics,
State Key Laboratory of Fluid Power
Transmission and Control,
Zhejiang University,
Hangzhou 310027, China

Wei Qiu Zhu

Department of Mechanics,
State Key Laboratory of Fluid Power
Transmission and Control,
Zhejiang University,
Hangzhou 310027, China
e-mail: wqzhu@zju.edu.cn

1Corresponding author.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received April 2, 2015; final manuscript received July 2, 2015; published online July 22, 2015. Editor: Yonggang Huang.

J. Appl. Mech 82(10), 101008 (Jul 22, 2015) Paper No: JAM-15-1180; doi: 10.1115/1.4031009 History: Received April 02, 2015

The responses of linear and nonlinear oscillators to fractional Gaussian noise (fGn) are studied. First, some preliminary concepts and properties of fractional Brownian motion (fBm) and fGn with Hurst index 1/2<H<1 are introduced. Then, the exact sample solution, correlation function, spectral density, and mean-square value of the response of linear oscillator to fGn are obtained. Based on the sample solution, it is proved that the long-range correlation index of displacement response of linear oscillator is the same as that of excitation fGn, i.e., 2-2H, while the velocity response has no such long-range correlation. An interesting discovery is that the ratio of kinetic energy to total energy decreases as increasing Hurst index H. Finally, for the responses of one and two degrees-of-freedom (DOF) nonlinear oscillators to fGn, the equivalent linearization method is applied to obtain the sample functions, correlation functions and mean-square values of the responses. Plenty of digital simulation results are obtained to support these solutions. It is shown that the approximate solution is effective for weakly nonlinear oscillators and it is feasible to apply the equivalent linearization to study multi-DOF weakly nonlinear oscillators.

FIGURES IN THIS ARTICLE
<>
Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.

References

Roberts, J. B. , and Spanos, P. D. , 1990, Random Vibration and Statistical Linearization, Wiley, New York.
Soong, T. T. , and Grigoriu, M. , 1993, Random Vibration of Mechanical and Structural System, PTR Prentice-Hall, Englewood Cliffs, NJ.
Lin, Y. K. , and Cai, G. Q. , 1995, Probabilistic Structural Dynamics, Advanced Theory and Applications, McGraw-Hill, New York.
Zhu, W. Q. , 2006, “Nonlinear Stochastic Dynamics and Control in Hamiltonian Formulation,” ASME Appl. Mech. Rev., 59(4), pp. 230–248. [CrossRef]
Yaglom, A. M. , 1958, “Correlation Theory of Processes With Random Stationary n th Increments,” AMS Transl., 2(8), pp. 87–141.
Mandelbrot, B. B. , and VanNess, J. W. , 1968, “Fractional Brownian Motions, Fractional Noises and Applications,” SIAM Rev., 10(4), pp. 422–437. [CrossRef]
Biagini, F. , Hu, Y. , Øksendal, B. , and Zhang, T. , 2008, Stochastic Calculus for Fractional Brownian Motion and Applications, Springer-Verlag, London.
Mishura, Y. S. , 2008, Stochastic Calculus for Fractional Brownian Motion and Related Processes, Springer-Verlag, Berlin.
Mandelbrot, B. B. , 1997, Fractals and Scaling in Finance: Discontinuity, Concentration, Risk, Springer Verlag, New York.
Cheridito, P. , 2003, “Arbitrage in Fractional Brownian Motion Models,” Finance Stochastic, 7(4), pp. 533–553. [CrossRef]
Leland, W. E. , Taqqu, M. S. , Willinger, W. , and Wilson, D. V. , 1994, “On the Self-Similar Nature of Ethernet Traffic,” IEEE/ACM Trans. Networking, 2(1), pp. 1–15. [CrossRef]
Yang, H. , Luo, G. B. , Karnchanaphanurach, P. , Louie, T. M. , Rech, I. , Cova, S. , Xun, L. , and Xie, X. S. , 2003, “Protein Conformational Dynamics Probed by Single-Molecule Electron Transfer,” Science, 302(5643), pp. 262–266. [CrossRef] [PubMed]
Kou, S. C. , and Xie, X. S. , 2004, “Generalized Langevin Equation With Fractional Gaussian Noise: Subdiffusion Within a Single Protein Molecule,” Phys. Rev. Lett., 93(18), p. 180603. [CrossRef] [PubMed]
Cheridito, P. , Kawaguchi, H. , and Maejima, M. , 2003, “Fractional Ornstein–Uhlenbeck Processes,” Electron. J. Probab., 8(3), pp. 1–14.
Kaarakka, T. , and Salminen, P. , 2011, “On Fractional Ornstein–Uhlenbeck Processes,” Commun. Stochastic Anal., 5(1), pp. 121–133.
Dung, N. T. , 2012, “On Delayed Logistic Equation Driven by Fractional Brownian Motion,” ASME J. Comput. Nonlinear Dyn., 7(3), p. 031005. [CrossRef]
Mandelbrot, B. B. , 1982, The Fractal Geometry of Nature, W. H. Freeman, San Francisco, CA.
Papoulis, A. , 1991, Probability, Random Variables and Stochastic Processes, 3rd ed., McGraw-Hill, New York.
Goldberg, R. R. , 2008, Fourier Transforms, Cambridge University, Cambridge, MA.
Uchaikin, V. V. , 2012, Fractional Derivatives for Physicists and Engineers: Vol. I Background and Theory, Springer-Verlag, Berlin.
Wong, E. , and Zakai, M. , 1965, “On the Relation Between Ordinary and Stochastic Differential Equations,” Int. J. Eng. Sci., 3(2), pp. 213–219. [CrossRef]

Figures

Grahic Jump Location
Fig. 4

Autocorrelation function RX(τ) of response X(t) of oscillator (27). The system parameters are k=2, D=0.1, γ=0.5. Part of data in (a) (i.e., 0.1< τ < 600) are zoomed in and show in (b). The analytical results come from Eq. (39).

Grahic Jump Location
Fig. 1

Mean-square displacements E[X2] and mean-square velocity E[V2] of oscillator (27). The system parameters are k=2, D=0.1, γ=2ξk. (a) Analytical solution for mean-square E[X2] from Eq. (29); (b) analytical solution for mean-square E[V2] from Eq. (31).

Grahic Jump Location
Fig. 2

Ratio of kinetic energy to total energy rK (a) and ratio of potential energy to total energy rP (b) with different damping ratio ξ. The system parameters are the same as those in Fig. 1.

Grahic Jump Location
Fig. 3

Plot of ratio rK and rP against damping ratio ξ with different Hurst index H. The system parameters are the same as those in Fig. 2.

Grahic Jump Location
Fig. 5

Mean-square of displacements E[X2] and velocity E[V2] of oscillator (44). The system parameters are k=2, D=0.1, λ=0.75k. Analytical solutions come from Eq. (51).

Grahic Jump Location
Fig. 6

The influence of nonlinear parameter λ on mean-square displacement E[X2]. H=0.75. The other system parameters are the same as Fig. 5.

Grahic Jump Location
Fig. 7

Two sample functions calculated from Eq. (35)X(t)|k→ke and simulated from Eq. (27). The system parameters are k=2, D=0.1, γ=k, λ=3k, H=0.75, initial state [X(0),V(0)]T=[0,0]T.

Grahic Jump Location
Fig. 10

The influence of nonlinear parameter λ1 on mean-square displacement E[X12]. λ2=0, H=0.75. The other system parameters are the same as Fig. 9.

Grahic Jump Location
Fig. 11

Two sample functions calculated from Eq. (21)X1(t)|C→Ce,K→Ke and simulated from Eq. (52). k1=k2=2, ks=0.5, D=0.1, λ1=3k1, λ2=2k2, γ=k1, H=0.75, initial state [X1(0),V1(0),X2(0),V2(0)]T=0.

Grahic Jump Location
Fig. 8

Autocorrelation function RX(τ)|k→ke of response X(t) of nonlinear oscillator (44). Where RX(τ)|k→ke means RX(τ) is determined in Eq. (39) and k is replaced by ke in Eq. (50). k=2, D=0.1, γ=k, and H=0.75.

Grahic Jump Location
Fig. 9

Mean-square of displacements E[X12] and E[X22] of oscillator (52). The system parameters are k1=k2=2, ks=0.5, D=0.1, λ1=0.75k1, λ2=0.5k2. The analytical solutions come from Eq. (55).

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In