Shock and Vibration 11 (2004) 685–692 685 IOS Press

Parameter identification of weakly nonlinear vibration system in frequency domain

Jiehua Penga,b,∗, Jiashi Tanga and Zili Chena,c aDepartment of Mechanics, University, , 410082, P.R. bDepartment of physics, University, Shaoyang, 422000, P.R. China cCollege of Architectural Engineering, Central South Forestry University, , 412008, P.R. China

Received 1 November 2001 Revised 11 December 2003

Abstract. A new method of identifying parameters of nonlinearly vibrating system in frequency domain is presented in this paper. The problems of parameter identification of the nonlinear dynamic system with nonlinear elastic force or nonlinear damping force are discussed. In the method, the mathematic model of parameter identification is frequency response function. Firstly, by means of perturbation method the frequency response function of weakly nonlinear vibration system is derived. Next, a parameter transformation is made and the frequency response function becomes a linear function of the new parameters. Then, based on this function and with the least square method, physical parameters of the system are identified. Finally, the applicability of the proposed technique is confirmed by numerical simulation.

Keywords: Parameter identification, weakly nonlinear vibration system, least square method, frequency domain

1. Introduction

The parameter identification is one of the important problems of the vibration research. There are many methods of identifying parameters of the linear system. But the parameter identification of nonlinear system is much more difficult than the linear system. Hence, this problem attracts much attention. A. H. Nayfeh has presented an excellent method of parametric identification of nonlinear system that exploits nonlinear resonances and comparisons of the system to be identified with those of known system [1], but the difficult usually appears when we find the known system. Tang et al. [2–4] have suggested time domain method and phase plane method for identifying parameters of nonlinear system which are effective for integrable nonlinear system. K. Yasuda et al. [5,6] have proposed an experiment method of parameter identification which is useful for geometrically nonlinear systems. R. Bachmayer et al. [7] have established adaptive parameter identification method for marine thruster which is an on-line technique for adaptive identification. There are some other methods of parameter identification, such as orthogonal function method [8], wavelet-based technique [9], Volterraseries method [10], cross-correlation analysis [11], power spectrum technique [12], and so on. A new method of identifying parameters of the weakly nonlinear vibration system in the frequency domain is presented in this paper. We discuss two kinds of nonlinear vibration systems with nonlinear elastic force or nonlinear damping force respectively. First, by using the multiple scale method, we obtain the frequency response function of the system which is regarded as a mathematical model of the curve fitting. Next, a new set of parameters are introduced to replace the parameters to be identified and the frequency response function is transformed into the linear function of the new parameters so that the least square method can be used to identify the parameters. This

∗Corresponding author. E-mail: [email protected].

ISSN 1070-9622/04/$17.00  2004 – IOS Press and the authors. All rights reserved 686 J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain is the key step of our method. Then, with the least square method the normal equations are established and they are linear equations of the new parameters. From the solution of normal equations we can get the values of the parameters to be identified. The method we present is a combination of the perturbation method and the least square method, hence, the method has the same precision with those of perturbation method and least square method. On the other hand, the perturbation method is an effective method to deal with weakly nonlinear system and the least square method is an excellent method of curve fitting such that this parameter identification method has a wide application in the study of weakly nonlinear vibration system. Although frequency response function is the first order approximation of the relationship between amplitude and frequency of nonlinear system, it is the most useful and important mathematical model to study the property of nonlinear vibration system [13]. It can be determined by experiment with good accuracy because in practice any complicated nonlinear system will finally move in a steady motion state if the excitation is harmonic and the system doesn’t come into chaotic motion. Compared with response function in time domain, frequency response function is independent of initial conditions and the test to determine this relationship possesses good repeatability and stability, such that this method of parameter identification can be easily realized by experiment.

2. Principle of the method

The governing equation of the vibration system with nonlinear damping force is as follows mx¨ + cx˙ + f(˙x)+kx = G cos Ωt (1) The governing equation of the vibration system with nonlinear elastic force is as follows mx¨ + cx˙ + kx + f(x)=G cos Ωt (2) where f(˙x) is a nonlinear function of x˙, f(x) is a nonlinear function of x. Rewrite Eqs (1) and (2) into the standard form 2 x¨ + ω x = ε[−2nx˙ − f1(˙x)+g cos Ωt] (3) and 2 x¨ + ω x = ε[−2nx˙ − f1(x)+g cos Ωt] (4) k c where ε being a small positive parameter, ω = m being the linear natural frequency of the system, n = 2m , f(˙x) f(x) G f1(˙x)= m , f1(x)= m and g = m . When Ω approaches to ω, we introduce a detuning parameter σ, letting Ω=ω + εσ (5) We can find the approximate frequency response function of the nonlinear system Eqs (3) and (4) by using approximate analytical method, such as the methods of iteration, multiple scales, KBM and so on [13–15]. The general approximate frequency response function of Eqs (3) and (4) can be written in the form F (ω,n,...,a,σ)=0 (6) Where ω,n,..., are physical parameters to be identified of the system Eqs (3) and (4), a is the amplitude, σ is the detuning parameter whose dimension is the same with frequency. Equation (6) plays key role in revealing the properties of nonlinear vibration system and here we take it as a mathematical model of the curve fitting. But F (ω,n,...,a,σ) is a nonlinear function of parameters (ω,n,...,) we must replace the parameters (ω,n,...,) with a new set of parameters (A,B,...,) before employing the least square method to identify the parameters, such that F (ω,n,...,a,σ) can be transformed into a linear function of (A,B,...,). That is, if we chose (ω,n,...,)=R(A,B,...,) (7) then J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain 687

F (ω,n,...,a,σ)=F1(A,B,...,a,σ)=0 (8)

Where F1(A,B,...,a,σ) is a linear function of parameters (A,B,...,). Transform function R(A,B,...,) is determined by the frequency response function of the systems. We use the test data (σ i,ai)(i =1, 2,...,m) for fitting the curve. Substituting σi,ai into Eq. (8) yields small errors

Ei = F1(A,B,...,ai,σi) i =1, 2,...,m The sum of the square of all errors is m m 2 2 E = Ei = [F1(A,B,...,ai,σi)] (9) i=1 i=1 It is obvious that the parameters ω,n...will fit to make the totle square errors E get the minimum. This means that we can use the method of the minimum squares to identify the parameters of the system. By determining the extreme value of E and solving a set of normal equations resulted from the extreme value of E, we can identify parameters ω,n,...of the nonlinear system.

3. The system with nonlinear damping

Let us consider a weakly nonlinear vibration system with a small cubic nonlinear damping force, its governing equation is as follows. 3 mx¨ + cx˙ + c1x˙ + kx = G cos Ωt (10) Rewrite the standard equation in the form mx¨ + ω2x = ε(−2nx˙ − αx˙ 3 + g cos Ωt) (11) where k g c G ω2 = , 2n = ,α= ,g= (12) m m m m ω,n and α are the physical parameters to be identified of the system Eq. (11). We begin by assuming an expansion of the solution having the form

x(t, ε)=x0(T0,T1)+εx1(T0,T1)+... (13) New independent variables are introduced according to n Tn = ε t for n =0, 1,... (14) Substituting Eq. (13) into Eq. (11) and equating the coefficients of ε i, (i =0, 1, 2,...) to zero, we obtain 2 2 D0x0 + ω x0 =0 (15)

2 2 3 D0x1 + ω x1 = −2D0D1x0 − 2nD0x0 − α(D0x0) + g cos(ωT0 + σT1) (16) The solution of the Eq. (15) is

x0 = a cos(ωT0 + ϕ) (17) Substituting Eq. (17) into Eq. (16) leads to 1 2 + 2 =2 sin( + )+2 cos( + )+2 sin( + )+ 3 3[3 sin( + ) D0x1 ω x1 ωa ωT0 ϕ ωaβ ωT0 ϕ nωa ωT0 ϕ 4αa ω ωT0 ϕ

− sin 3(ωT0 + ϕ)] + g cos(ωT0 + ϕ)cos(σT1 − ϕ) − g sin(ωT0 + ϕ)sin(σT1 − ϕ)

To eliminate secular terms from x1, we must put 688 J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain  3 3 3 2ωa +2naω + 4 αa ω −g sin(σT1 − ϕ)=0 (18)  2ωaϕ + g cos(σT1 − ϕ)=0 Letting

σT1 − ϕ = γ (19) Determining the derivatives with respect to T1 yields ϕ = σ − γ . Considering the steady-state response, that is, a =0,γ =0,wehave  2naω + 3 αa3ω3 − g sin γ =0 4 (20) 2ωaσ + g cos γ =0 The result by eliminating γ in Eq. (20) is 9 4 2 2 2 +3 4 4 + 2 6 6 +4 2 2 2 − 2 =0 n ω a αnω a 16α ω a ω σ a g (21) This is the first-order approximation frequency response function of the weakly nonlinear system Eq. (11). For identifying physical parameters ω, n and α of the system with test data {σ i,ai}(i =1, 2,...,m), Following transformation should be introduced 9 =4 2 2 =3 4 = 2 6 =4 2 A n ω ,B αnω ,C 16α ω ,D ω (22) therefore Eq. (21) becomes Aa2 + Ba4 + Ca6 + Dσ2a2 − g2 =0 (23) Equation (23) is a frequency response equation of the nonlinear system Eq. (11). We regard it as a mathematical model of the parameter identification. Setting the testing values of (σ, a) are {σ i,ai}(i =1, 2 ...,m) and substituting them into Eq. (23) generating errors as follow 2 4 6 2 2 2 Ei = ai A + ai B + ai C + σi ai D − g i =1, 2 ...,m The sum of the square of all errors is m 2 4 6 2 2 2 2 E = (ai A + ai B + ai C + σi ai D − g ) (24) i=1 The best way to determine A, B, C and D is choosing those values of them that make the sum of the square of all errors E be minimum. This means that the difference between the theoretical result and the experiment result is least. Noting that the error E is the function of A, B, C and D, the necessary conditions of E getting minimum are that the partial derivatives of E with respect to A, B, C and D equal zero, respectively, that is, ∂E ∂E ∂E ∂E =0, =0, =0, =0 ∂A ∂B ∂C ∂D namely, m 2 4 6 2 2 2 2 (ai A + ai B + ai C + σi ai D − g )ai =0 i=1 m 2 4 6 2 2 2 4 (ai A + ai B + ai C + σi ai D − g )ai =0 i=1 m 2 4 6 2 2 2 6 (ai A + ai B + ai C + σi ai D − g )ai =0 i=1 m 2 4 6 2 2 2 2 2 (ai A + ai B + ai C + σi ai D − g )σi ai =0 i=1 J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain 689

These equations are usually called normal equations in numerical analysis. They are linear algebraic equations of parameters A, B, C, D. According the theory of numerical analysis, the solution of the normal equations is sole and it makes the errors E get minimum definitively [16]. For the simplification we write the normal equations in the matrix form, Q y = s (25) where y =(ABCD)T   m m m m T 2 2 2 4 2 6 2 2 2 s = g ai g ai g ai g σi ai i=1 i=1 i=1 i=1 and   m m m m 4 6 8 2 4  ai ai ai σi ai     i=1 i=1 i=1 i=1  m m m m   6 8 10 2 6   ai ai ai σi ai    =  i=1 i=1 i=1 i=1  Q m m m m  (26)  8 10 12 2 8   ai ai ai σi ai   i=1 i=1 i=1 i=1  m m m m   2 2 2 6 2 6 4 4  σi ai σi ai σi ai σi ai i=1 i=1 i=1 i=1 As mentioned above, the solution of normal equation is sole and this means matrix Q is reversible. Solving matrix Eq. (25) yields y = Q−1 s. Substituting the values of A, B, C and D into Eq. (22) we obtain physical parameters of the nonlinear system in the form √ √ D A B 4 C ω = ,n= ,α= or α = (27) 2 D 2nω4 3ω3

4. The system with cubic elastic force

In the case of a nonlinear vibration system with a small cubic nonlinear elastic force and linear viscous damping, the governing equation of it is 3 mx¨ + cx˙ + kx + k1x = G cos Ωt (28) Equation (28) can be rewritten as x¨ + ω2x = ε(−2nx˙ − βx3 + g cos Ωt) (29)

k c k1 G ω2 = , 2n = ,β= ,g= . (30) m m m m Using the method of multiple scales or KBM, we can derive the frequency response equation of Eq. (29) as follows 9 4 2 2 2 +4 2 2 2 − 3 4 + 2 6 − 2 =0 n ω a ω σ a βωσa 16β a g (31) Similarly, A, B, C and D are introduced, Eq. (31) becomes Aa2 + Bσ2a2 + Cσa4 + Da6 − g2 =0 (32) Where 690 J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain 9 =4 2 2 =4 2 = −3 = 2 A n ω ,B ω ,C βω,D 16β . (33)

The sum of the square of all error arising in the procedure of substituting the testing data {σ i,ai}(i =1, 2,...,m) into Eq. (32) is m 2 2 2 4 6 2 2 E = (ai A + σi ai B + σiai C + ai D − g ) (34) i=1 The normal equation corresponding to the minimum of error function E can be obtained as follows       m m m m m 4 2 4 6 8 A 2 2  ai σi ai σiai ai    g ai         i=1 i=1 i=1 i=1     i=1  m m m m     m   2 4 4 4 3 6 2 8     2 2 2   σi ai σi ai σi ai σi ai   B  g σi ai         i=1 i=1 i=1 i=1    =  i=1  m m m m     m  (35)  6 3 6 2 8 10     2 4   σiai σi ai σi ai σiai   C  g σiai   i=1 i=1 i=1 i=1     i=1  m m m m     m   8 2 8 10 12     2 6  ai σi ai σiai ai g a i=1 i=1 i=1 i=1 D i=1 Solving Eq. (35) we obtain the values of A, B, C and D. Therefore physical parameters of the nonlinear system Eq. (29) are √ B A C 4√ ω = ,n= ,β= − , or,β= D (36) 2 B 3ω 3

5. Example and discussion

With the method of numerical simulation we can check up the precision of the method. We take Eq. (29) for example, letting ω =5,n=0.1 β =0.4 g =9.8 By substituting them into Eq. (29) and employing multiple scales method, we get the frequency response equation as follows a2 + 100σ2a2 − 6.0σa4 +0.09a6 − 96.04 = 0 From this equation we can get the data as in Table 1. Noting that in the multivalue section in Table 1 there are several amplitudes for each value of σ and it is consistent to theoretical results. The values in the blod correspond to the unstable state of the system as σ increases and mutation occurs at the point σ =0.75. Substituting {σi,ai} in Table 1 into Eq. (35), except the values in dark area, and solving them we get the values of A, B, C and D as follows A =1.0387 B =99.8408 C = −5.9802 D =0.0891 Substituting them into Eq. (36), we obtain the physical parameters of the nonlinear system as follows ω =4.996 n =0.102 β =0.399 These results of ω,n and β are very close to the values of them we set above, thus the accuracy of the method presented in this paper is excellent. This conclusion can be explained by the principle of the method. From Section II we see our method is the combination of perturbation techniques and least square method which are effective methods J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain 691

Table 1 Frequency response data simulated by computer σ a −1 0.49394 −0.5 0.93359 −0.4 1.11467 −0.3 1.35531 −0.2 1.66498 −0.1 2.03576 0 2.44104 0.1 2.85192 0.2 3.25015 0.3 3.62786 0.4 3.98285 2.96945 1.40922 0.5 4.31488 3.69048 1.04664 0.6 4.62311 4.2309 0.85208 0.7 4.90001 4.70466 0.72298 0.75 5.00000 0.67279 0.8 0.62938 0.9 0.55787 1 0.50126 to deal with the weakly nonlinear vibration system and curve fitting, respectively, and the frequency response function is obtained from the first-order approximation of multiple scale method or other approximation methods, so that the accuracy of the method is as well as the one of first-order approximation of perturbation methods or the least square method. Meanwhile, the mathematical model, frequency response function of the weakly nonlinear system, of parameter identification of present method is a most important model describing the properties of weakly nonlinear systems. Above discussion shows that this method is a useful means for the parameter identification of weakly nonlinear vibration system. Though the data {σi,ai}(i =1, 2 ...,m) in the example is obtained by numerical simulation, we can achieve these data in practice. Noting that, in essence, the detuning parameter σ is a small difference between the frequency of excitation and the linear natural frequency of the nonlinear system, thus, in practice test, we may change the frequency of excitation in a small range around the linear natural frequency, that is varying the detuning parameter σ. Recording down σ and the corresponding amplitude a of the system we obtain the test data {σ i,ai}(i =1, 2,...,m). This means that the present method of parameter identification can be realized by experiment. The study only involves one-degree-of –freedom nonlinear vibration system with cubic nonlinearity, but it is obvious that the proposed method of parameter identification can be applied to the system with higher nonlinearity or multi-degree-of-freedom nonlinear vibration systems.

References

[1] A.H. Nayfeh, Parametric identification of nonlinear systems, Computer and Structures 20 (1988), 487–493. [2] J. Tang and Q. Huo, Identifying Parameters of Nonlinear System in Time Domain by the Response of Free Vibration, Journal of Vibration Engineering 2 (1989), 62–66, in Chinese. [3] J. Tang and X. Li, Normal Perturbation Solution for Large Damped Autonomous Nonlinear Systems, Journal of Vibration and Shock 6(4) (1987). [4] J. Tang, Phase plane method for the parameter identification of strongly nonlinear systems, Shock and Vibration 8 (2001), 49–53. [5] K. Yasuda et al., Experimental identification technique of vibrating structures with geometrical nonlinear, ASME J. Appl. Mech. 64 (1997), 275–280. [6] K. Yasuda and K. Kamiya, Experimental identification technique of nonlinear beams in time domain, Nonlinear Dynamics 18 (1999), 185–202. [7] R. Bachmayer and L. Whitcomb, Adaptive parameter identification of an accurate nonlinear dynamical model for marine thrusters, J. Dynamic Systems, Measurement, and Control. 125 (2003), 489–491. [8] R. P. Pacheco and V. Steffen Jr., Using orthogonal functions for identification and sensitivity analysis of mechanical systems, J Vibration and Control 8 (2002), 993–1021. [9] R. Ghanem and F. Romeo, A wavelet-based approach for model and parameter identification of nonlinear systems, International J. Nonlinear Mechanics 36 (2001), 835-859. 692 J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain

[10] A. Chatterjee and N. Vyas, Nonlinear parameter estimation using Volterra series with Multi-Tone excitation, XXIMAC Proceedings of the 20th International Model Analysis Conference, Los Angeles CA, 2002, pp. 880–885. [11] D.Q. Nykamp and D.L. Ringach, Full identification of a linear-nonlinear system via cross-correlation analysis, Journal of Vision 2 (2002), 1–11. [12] H. Robert and K. Laszl´ o,´ Nonlinear System Identification Ð Input-Output Modeling Approach, Dordrecht: Kluwer Academic Publishers, 1999. [13] A.H. Nayfeh and D.T. Mook, Nonlinear oscillations, Wiley: New York, 1979. [14] A.H. Nayfeh, Introduction to perturbation techniques, Wiley: New York, 1988. [15] W.T. Thomson, Theory of vibration with applications, New Jersey: Prentice-hall Inc., 1972. [16] J. Deng, R. Ge and Z. Cheng, Calculation Method, Xi-an: Xi-an Jiaotong University Publishing Inc., 1985, in Chinese. International Journal of

Rotating Machinery

International Journal of Journal of The Scientific Journal of Distributed Engineering World Journal Sensors Sensor Networks Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014

Journal of Control Science and Engineering

Advances in Civil Engineering Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014

Submit your manuscripts at http://www.hindawi.com

Journal of Journal of Electrical and Computer Robotics Engineering Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014

VLSI Design Advances in OptoElectronics International Journal of Modelling & International Journal of Simulation Aerospace Navigation and in Engineering Observation Engineering

Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2010 http://www.hindawi.com Hindawi Publishing Corporation http://www.hindawi.com Volume 2014

International Journal of International Journal of Antennas and Active and Passive Advances in Chemical Engineering Propagation Electronic Components Shock and Vibration Acoustics and Vibration

Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014