0
Research Papers

Stationary Response of Multidegree-of-Freedom Strongly Nonlinear Systems to Fractional Gaussian Noise

[+] Author and Article Information
Qiang Feng Lü, Mao Lin Deng

Department of Mechanics,
Zhejiang University,
Hangzhou 310027, China;
State Key Laboratory of Fluid Power and
Mechatronic Systems,
Zhejiang University,
Hangzhou 310027, China;
Key Laboratory of Soft Machines and
Smart Devices of Zhejiang Province,
Zhejiang University,
Hangzhou 310027, China

Wei Qiu Zhu

Department of Mechanics,
Zhejiang University,
Hangzhou 310027, China;
State Key Laboratory of Fluid Power and
Mechatronic Systems,
Zhejiang University,
Hangzhou 310027, China;
Key Laboratory of Soft Machines and
Smart Devices of Zhejiang Province,
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 22, 2017; final manuscript received July 25, 2017; published online August 18, 2017. Editor: Yonggang Huang.

J. Appl. Mech 84(10), 101001 (Aug 18, 2017) (14 pages) Paper No: JAM-17-1213; doi: 10.1115/1.4037409 History: Received April 22, 2017; Revised July 25, 2017

The stationary response of multidegree-of-freedom (MDOF) strongly nonlinear system to fractional Gaussian noise (FGN) with Hurst index 1/2 < H < 1 is studied. First, the system is modeled as FGN-excited and -dissipated Hamiltonian system. Based on the integrability and resonance of the associated Hamiltonian system, the system is divided into five classes: partially integrable and resonant, partially integrable and nonresonant, completely integrable and resonant, completely integrable and nonresonant, and nonintegrable. Then, the averaged fractional stochastic differential equations (SDEs) for five classes of quasi-Hamiltonian systems with lower dimension and involving only slowly varying processes are derived. Finally, the approximate stationary probability densities and other statistics of two example systems are obtained by numerical simulation of the averaged fractional SDEs to illustrate the application and compared with those from original systems to show the advantages of the proposed procedure.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.

References

Zhu, W. Q. , 2006, “ Nonlinear Stochastic Dynamics and Control in Hamiltonian Formulation,” ASME Appl. Mech. Rev., 59(4), pp. 230–248. [CrossRef]
Zeng, Y. , and Zhu, W. Q. , 2011, “ Stochastic Averaging of Quasi Non-Integrable Hamiltonian Systems Under Poisson White Noise Excitation,” ASME J. Appl. Mech., 78(2), p. 021002. [CrossRef]
Liu, W. Y. , and Zhu, W. Q. , 2015, “ Lyapunov Function Method for Analyzing Stability of Quasi-Hamiltonian Systems Under Combined Gaussian and Poisson White Noise Excitations,” Nonlinear Dyn., 81(4), pp. 1879–1893. [CrossRef]
Jia, W. T. , Zhu, W. Q. , Xu, Y. , and Liu, W. Y. , 2013, “ Stochastic Averaging of Quasi-Integrable and Resonant Hamiltonian Systems Under Combined Gaussian and Poisson White Noise Excitations,” ASME J. Appl. Mech., 81(4), p. 041009. [CrossRef]
Xu, Y. , Guo, R. , Liu, D. , Zhang, H. Q. , and Duan, J. Q. , 2014, “ Stochastic Averaging Principle for Dynamical Systems With Fractional Brownian Motion,” AIMS Discrete Continuous Dyn. Syst., Ser. B, 19(4), pp. 1197–1212. [CrossRef]
Xu, Y. , Duan, J. Q. , and Xu, W. , 2011, “ An Averaging Principle for Stochastic Dynamical Systems With Lévy Noise,” Physica D, 240(17), pp. 1395–1401. [CrossRef]
Calif, R. , and Schmitt, F. G. , 2012, “ Modeling of Atmospheric Wind Speed Sequence Using a Lognormal Continuous Stochastic Equation,” J. Wind Eng. Ind. Aerodyn., 109, pp. 1–8. [CrossRef]
Guo, W. J. , Wang, Y. X. , Xie, M. X. , and Cui, Y. J. , 2009, “ Modeling Oil Spill Trajectory in Coastal Waters Based on Fractional Brownian Motion,” Mar. Pollut. Bull., 58(9), pp. 1339–1346. [CrossRef] [PubMed]
Sliusarenko, O. Y. , Gonchar, V. Y. , Chechkin, A. V. , Sokolov, I. M. , and Metzler, R. , 2010, “ Kramers-Like Escape Driven by Fractional Gaussian Noise,” Phys. Rev. E., 81(4), p. 041119. [CrossRef]
Ai, B. , He, Y. , and Zhong, W. , 2010, “ Transport in Periodic Potentials Induced by Fractional Gaussian Noise,” Phys. Rev. E., 82(6), p. 061102. [CrossRef]
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]
Diop, M. A. , and Garrido-Atienza, M. J. , 2014, “ Retarded Evolution Systems Driven by Fractional Brownian Motion With Hurst Parameter H > 1/2,” Nonlinear Anal. Theor., 97, pp. 15–29. [CrossRef]
Deng, M. L. , and Zhu, W. Q. , 2015, “ Responses of Linear and Nonlinear Oscillators to Fractional Gaussian Noise With Hurst Index Between 1/2 and 1,” ASME J. Appl. Mech., 82(10), p. 101008. [CrossRef]
Lin, S. J. , 1995, “ Stochastic Analysis of Fractional Brownian Motions,” Stochastics, 55(1–2), pp. 121–140.
Decreusefond, L. , 1999, “ Stochastic Analysis of the Fractional Brownian Motion,” Potential Anal., 10(2), pp. 177–214. [CrossRef]
Zähle, M. , 1998, “ Integration With Respect to Fractal Functions and Stochastic Calculus—I,” Probab. Theory Relat. Fields, 111(3), pp. 333–374. [CrossRef]
Duncan, T. E. , Hu, Y. , and Pasik-Duncan, B. , 2000, “ Stochastic Calculus for Fractional Brownian Motion I,” Theory SIAM J. Optim., 38(2), pp. 582–612. [CrossRef]
Biagini, F. , Hu, Y. Z. , Øksendal, B. , and Zhang, T. S. , 2008, Stochastic Calculus for Fractional Brownian Motion and Applications, Springer-Verlag, London. [CrossRef]
Zhu, W. Q. , Huang, Z. L. , and Yang, Y. Q. , 1997, “ Stochastic Averaging of Quasi-Integrable Hamiltonian Systems,” ASME J. Appl. Mech., 64(4), pp. 975–984. [CrossRef]
Zhu, W. Q. , Huang, Z. L. , and Suzuki, Y. , 2002, “ Stochastic Averaging and Lyapunov Exponent of Quasi Partially Integrable Hamiltonian Systems,” Int. J. Nonlinear Mech., 37(3), pp. 419–437. [CrossRef]
Zhu, W. Q. , and Yang, Y. Q. , 1997, “ Stochastic Averaging of Quasi-Nonintegrable-Hamiltonian Systems,” ASME J. Appl. Mech., 64(1), pp. 157–164. [CrossRef]
Baker, C. J. , 1995, “ The Development of a Theoretical Model for the Windthrow of Plants,” J. Theor. Biol., 175(3), pp. 355–372. [CrossRef]
Peltola, H. , and Kellomäki, S. , 1993, “ A Mechanistic Model for Calculating Windthrow and Stem Breakage of Scots Pines at Stand Age,” Silva Fennica, 27(2), pp. 99–111.
Letchford, C. W. , 2001, “ Wind Loads on Rectangular Signboards and Hoardings,” J. Wind Eng. Ind. Aerodyn., 89(2), pp. 135–151. [CrossRef]
Simiu, E. , and Scanlan, R. H. , 1996, Wind Effects on Structures: Fundamentals and Application to Design, Wiley, New York.
Kaimal, J. C. , Wyngaard, J. C. , Izumi, Y. , and Cote, O. R. , 1972, “ Spectral Characteristics of Surface-Layer Turbulence,” Q. J. R. Meteorol. Soc., 98(417), pp. 563–589. [CrossRef]
Davenport, A. G. , 1961, “ The Spectrum of Horizontal Gustiness Near the Ground in High Winds,” Q. J. R. Meteorol. Soc., 87(372), pp. 194–211. [CrossRef]
Von Karman, T. , 1948, “ Progress in the Statistical Theory of Turbulence,” Proc. Natl. Acad. Sci. U.S.A., 34(11), pp. 530–539. [CrossRef] [PubMed]
Vepa, R. , 2011, “ Nonlinear Optimal Control of a Wind Turbine Generator,” IEEE Trans. Energy Convers., 26(2), pp. 468–478. [CrossRef]
van der Male, P. , van Dalen, K. N. , and Metrikine, A. V. , 2016, “ The Effect of the Nonlinear Velocity and History Dependencies of the Aerodynamic Force on the Dynamic Response of a Rotating Wind Turbine Blade,” J. Sound. Vib, 383, pp. 191–209. [CrossRef]
Mishura, Y. S. , 2008, Stochastic Calculus for Fractional Brownian Motion and Related Processes, Springer-Verlag, Berlin. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Samples of processes I1,I2,H3,Φ, Θ2,Q3,Q4,P4 of system (43)

Grahic Jump Location
Fig. 7

Contour of stationary probability densities p(q1,q2),p(q1,q3) of system (43): (a), (c) simulated from fractional SDEs (47); (b), (d) simulated from original system (43)

Grahic Jump Location
Fig. 8

Stationary probability densities p(q1), p(q2), p(q3), p(q4) of system (43)

Grahic Jump Location
Fig. 9

Stationary probability density of phase difference p(ϕ) of system (43)

Grahic Jump Location
Fig. 10

Mean-square values E[Q12],E[Q22],E[Q32],E[Q42] of system (43) as a function of Hurst index H

Grahic Jump Location
Fig. 11

The autocorrelation function R(τ) of FGN WH(t) with Hurst index 0.5, 0.6, 0.8, and 0.98

Grahic Jump Location
Fig. 12

Samples of FGN with different Hurst index 0.5, 0.6, 0.8, and 0.98

Grahic Jump Location
Fig. 13

Autocorrelation function RQ3(τ) of response Q3(t) of system (43) under FGN excitation with different Hurst index 0.5, 0.9, and 0.98. Part of (a) (i.e., 30<τ<100) are zoomed in (b).

Grahic Jump Location
Fig. 14

Sample functions of processes I1,I2,H3,Φ of system (43) under FGN excitation with different Hurst index 0.5, 0.8, 0.98, and 1

Grahic Jump Location
Fig. 15

The mechanical model of system (58)

Grahic Jump Location
Fig. 2

stationary probability densities p(I1), p(H3) of system (43) under FGN or GWN excitation. Solid line is obtained from the method in Ref. [20]. Other curves are obtained from the method proposed in this paper.

Grahic Jump Location
Fig. 3

Contour of stationary probability densities p(I1,I2), p(I1,H3) of system (43): (a), (c) simulated from fractional SDEs (47); (b), (d) simulated from original system (43)

Grahic Jump Location
Fig. 4

Stationary probability densities p(I1),p(I2),p(H3) of system (43)

Grahic Jump Location
Fig. 5

Mean-values E[I1],E[I2],E[H3] of system (43) as a function of Hurst index H

Grahic Jump Location
Fig. 6

Mean-square values E[I12],E[I22],E[H32] of system (43) as a function of Hurst index H

Grahic Jump Location
Fig. 16

The Kaimal spectrum and the spectral density of FGN with parameters: v*2=0.0264, z=10 m, v¯(z)=40 m/s, and Hurst=0.973

Grahic Jump Location
Fig. 17

Stationary probability densities p(H) of system (58)

Grahic Jump Location
Fig. 18

Stationary probability densities p(q1), p(q2) of system (58)

Grahic Jump Location
Fig. 19

Sample of X(t) generated from Eqs. (A5) and (A6), respectively, where α(X,t)=−γX+D/ω2X, β(X,t)=D/ωX with parameter γ=0.03,D=0.1,ω=1

Tables

Errata

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