A New Simultaneous Identification of the Harmonic Excitations and Nonlinear Damping of Forced Damped Nonlinear Oscillations : A Parametric Approach

This paper presents a novel, general method that is aimed at identifying both the harmonic excitations and nonlinear damping of forced nonlinear oscillation systems. An inverse formalism is developed for the identification method, which is shown to be well posed.That is, the inverse formalism has a unique solution and it depends continuously on the data.The unique solution is derived in the closed form. In addition to the mathematical analysis, a numerical example involving a highly nonlinear system is examined for demonstrating the workability of the proposed method for simultaneously identifying both the two aforementioned excitations and nonlinear damping.


Introduction
Over the past decades, a considerable number of studies have been made on forced oscillations arising in many diverse fields of science and engineering.This is mainly because the studies are closely related in the investigation of system performance such as the effective control and the estimation of the structure lifetime.Thus, as a forward problem, many researchers investigated nonlinear responses and stability analysis of the second-order nonlinear differential equations of a general form which have a periodic forcing term.However, these studies strongly depend on how accurately we can estimate the characteristics of forced loading condition or/and nonlinear damping.This is regarded as an inverse problem of the system identification, which is not only essential but crucial for the forward problem.
There are a variety of researches regarding the system identification devoted to the study about damping [1].For example, Iourtchenko and Dimentberg [2] suggested a procedure for in-service identification of the damping based on the stochastic averaging method from a measured stationary response.Jang et al. [3,4] proposed a new nonparametric method to recover the functional form of both nonlinear damping and nonlinear restoring forces of the dynamic system using the transient data.As for the force identification, there are also a lot of articles [5].Doyle [6] studied a waveletbased deconvolution method to recover impact forces using experimentally measured responses.Zen et al. [7] proposed a method to obtain contact forces from experimental measurement of a single impact.Jang et al. [8] introduced a new technique for the identification of external loading of long time duration, which is nonharmonic but periodic, acting on a nonlinear dynamic system.
As mentioned above, although a large number of studies have been carried out on the system identification of nonlinear damping or external force, little is known about their combinations, that is, the simultaneous identification of nonlinear damping as well as external force.In fact, many researchers in engineering fields have wondered if the problem of the simultaneous identification is mathematically solvable.Motivated by this, an inverse formalism is, in this present work, mathematically developed for identifying both the nonlinear damping and external forces in a forced nonlinear system.
In this study, we propose a method for the simultaneous identification, which is a new procedure with remarkable improvements.First, based on the concept of a zero-crossing time, we develop an inverse formalism which is linear, even if the original dynamical equation of motion is nonlinear.Second, we prove that the inverse formalism has a unique solution and its closed form solution is derived, which makes the proposed method quite simple and straightforward.Third, the inverse formalism developed is shown to have stability properties of solution.This means that the method by the inverse formalism work without the use of any regularization, so that it can have a wide range of physical applications.This contrasts the usual inverse problems illposed in the sense of stability [3,4,[8][9][10], which require an additional process of regularization.
Finally, given the model nonlinear differential equation concerning the quadratic nonlinear damping and cubic nonlinear restoring characteristics, the proposed method's workability for simultaneously identifying both excitations and nonlinear damping is demonstrated through a numerical example.

Forced Damped Nonlinear Oscillation
We consider a nonlinear damped dynamical system, subject to an external harmonic load Γ ⋅ cos , of an amplitude Γ > 0 and frequency  > 0. The equation of motion is a second-order nonlinear ordinary differential equation for the response   q +  ( q ) +  () = Γ ⋅ cos ,  > 0.
Here, the constant  > 0 represents the mass of a particle of an oscillator.The term () is a general nonlinear restoring force and the term ( q ), appearing in (1), is the nonlinear damping of a system, which is assumed to have a quadratic nonlinear term where  1 and  2 are the linear and nonlinear damping coefficients, respectively.The functional form, (2), of nonlinearity is often encountered in the area of engineering field [11].
The sum of the inertia force  q and restoring force () and the derivative of the sum will play an important role in identifying both the excitation Γ and nonlinear damping ( q ) in (1).We define the sum as  (, q ) ≡  q +  () .
The derivative of  with respect to time then becomes Noting that (1) is equivalent to we can express J in terms of the excitation and damping as follows: by differentiating (5) with respect to time.
In what follows, we assume that (1) has two dynamic responses of (= 2/)-periodic limit cycle solutions   (), which are forced by The response analysis of the nonlinear damped dynamical system considered can be regarded as a forward problem, whose solution procedures are well established theoretically as well as numerically.However, an inverse problem can now be formulated in the form of a question: given a motion response of the nonlinear dynamical system, is it possible to identify the excitations and nonlinear damping of the system at the same time?The question may be interesting mathematically as well as practically.In order to answer the question, we will begin with introducing a zero-crossing time in the next section.

Some Lemmas on Zero-Crossing Time
This section is concerned with a zero-crossing time associated with motion responses of a system.Thereby, we derive some lemmas which are essential in the present study.We start by defining a zero-crossing time in terms of  for  = 1, 2, denoted by ( ze  )  , satisfying Here,   stands for the sum (3) which corresponds to   , that is, a function from (0, ∞) to R Physically, the zero-crossing time defined in ( 8) is a useful concept for finding the force balancing between the excitation and the damping in a forced damped oscillatory motion of (7) as discussed in the following lemma.

Lemma 1. Suppose that (𝑡 ze
)  ,  = 1, 2, is a zero-crossing time in terms of  for the system (7).Then the following holds: Proof.A simple substitution of ( 9) into (7) yields Since (11) remains valid for 0 <  < ∞, it should also hold for the zero-crossing time of ( ze  )  > 0. However, the left hand side of (11) vanishes for the zero-crossing time ( ze  )  due to the definition (8).That is, which immediately proves (10).
Journal of Applied Mathematics 3 Next, we introduce a zero-crossing time in terms of acceleration, ( ze acc where q  denotes an acceleration response of (7).We notice that the time rate of change of the damping is always zero at a zero-crossing time in terms of acceleration, which shall be examined in Lemma 2.
Lemma 2. For the function J defined by Here ( ze acc )  refers to a zero-crossing time in terms of acceleration.
Proof.In the same way as Lemma 1, (14) should hold for the zero-crossing times in terms of acceleration, that is,  = ( ze acc )  , for  = 1, 2. However, ( 14) is rewritten as in which the second term of the right hand side becomes zero because the acceleration vanishes when  = ( ze acc )  by the definition (13).It follows that the above equation becomes (15), which completes the proof.

Inverse Formalism
It is a high time that we develop an inverse formalism for the inverse problem raised by Section 2, based on the lemmas in the previous section.We recall that the excitations and nonlinear damping are parameterized by Γ 1 , Γ 2 ,  1 , and  2 in Section 2. This results in that the developed inverse formalism also involves the parameters of Γ 1 , Γ 2 ,  1 , and  2 as illustrated in Lemma 3.
Lemma 3 (inverse formalism).Equations (10) and (15) in Lemmas 1 and 2 constitute simultaneous linear algebraic equations for the four unknowns Γ 1 , Γ 2 ,  1 , and  2 .These can be written in a compact way as follows: Here, Α is an operator from R 4 to R 4 whose matrix is and x, b ∈ R 4 are represented as (one will call them a solution vector and data vector, resp., in this paper) Proof.Equations ( 10) and (15) in Lemmas 1 and 2 can be rewritten as the following: It is straightforward to show that the above four equations are equivalent to the operator equation ( 17) from the definition of the operator A, the vectors x and b of ( 18) and ( 19).This completes the proof.

Method
In the previous section, we developed an inverse formalism (17), whose solution procedure is to be established in this section.The formal expression for the solution may be symbolically written as where the notation Α −1 denotes the inverse operator for the linear operator Α of (18).According to linear algebra, it is well known that the inverse operator is bounded (or continuous), because it is defined on finite dimensional spaces.Thus, the solution x in (24) depends continuously on the data vector b: that is, the solution x is stable.The construction of x is examined through the following theorems.
Proof.It is easy to notice that the formula of (25) immediately follows from the relations of ( 22) and (23).We first divide both sides of (22) by − ⋅ sin[( ze acc ) 1 ] to obtain Γ 1 as follows: Similarly, we have from ( 23).This completes the proof.
The formulas of the damping coefficients  1 and  2 in (2) can also be obtained in an essentially simple way, as to be examined in Theorem 5.

Theorem 5 (nonlinear damping).
The linear damping coefficient  1 and the nonlinear damping coefficient  2 in (2) are determined as Proof.It suffices to multiply (20) and ( 21) by the corresponding factors, then subtract them from each other.As a result of taking (25) into account, we find  1 of (28) and then  2 of (29), which completes the proof.
Remark 6.Because of the stability properties of solution x, our solution procedure does not concern any kind of regularization, which is necessary for the usual inverse problems arising in science and engineering.
As to the solvability and uniqueness condition, the operator Α in (18) should satisfy det A ̸ = 0, which is to be discussed in the following theorem.

Theorem 7 (uniqueness). Let one denote by
Then, the excitations and nonlinear damping in (7) can be identified in a unique way whenever Δ ̸ = 0.
Proof.It suffices to notice that the formulas (25), (28), and (29) are valid if and only if the corresponding denominators are different from zero.After writing these restrictions, we obtain a condition equivalent to Δ ̸ = 0, which completes the proof.Remark 8. Theorem 7 tells us that the stable solution x is unique.This implies that the inverse formalism developed in the present study is a well-posed problem.
Our method for the simultaneous identification consists of the closed form solutions for the inverse formalism as in Theorems 4 and 5.However, the solutions are still expressed through the unknown zero crossing times.We thus are required to determine the unknown times for realizing the solutions, which is discussed in more detail in the next section.

Numerical Example
In this section, we will examine the workability of the proposed method to identify the excitations and nonlinear damping through numerical experiments.For that, let us assume that the restoring term () in (1) involves the cubic nonlinearity, that is, where the values for  1 and  2 are chosen as 40 and 10, respectively.Then, the responses for the system of (1) can be observed as following subsection.

System Responses.
With the initial conditions imposed as we integrate ( 7) by employing Runge-Kutta schemes to seek the system responses of (7).Here, the mass is normalized as unit and the coefficients of  1 and  2 are set to be 4 and 2, respectively.The excitation frequency  is taken as 4 with Γ 1 = 1 and Γ 2 = 2.The calculated responses of displacement and velocity of the system are depicted in Figure 1 and their phase diagrams in Figure 2.

Recovering Solutions by Zero-Crossing Time.
As mentioned in Section 5, even if we derived the closed form solutions for the present identification, they depend on the unknown two zero-crossing times, namely, ( ze  )  and ( ze acc )  in ( 8) and (13), respectively.
As a first case, we attempt to find the zero-crossing times from the steady-state region.By a proper numerical   3. We then substitute the zero-crossing times into (25), (28), and (29) for realizing the closed form solutions of identification.Table 1 presents the numerical values for the solutions by the proposed method, which is compared with the exact ones.In addition, timesimulated oscillatory motions using the identified excitations and damping by the proposed method are depicted in Figure 6, where the real (exact) time-simulated oscillatory motions also appear for comparison.
In a similar manner, as a second case, Figure 4 shows the zero-crossing times selected from the transient region, which are substituted again into (25), (28), and (29) for the identification.The results, compared with the exact ones, are  )  from the steady-state region for each  = 1, 2 as shown in Figure 5. Their corresponding results are listed in Table 3.For all three cases, the values for the condition numbers of Α in (18) are of (10 2 ), meaning that our system equation ( 17) is well conditioned.

Concluding Remarks
As discussed in Section 1 of introduction, many of engineers have doubted if it is possible to identify not only the external forcing but also nonlinear damping characteristics by only using system responses.This is a basic idea underlying this study.We have introduced the concept of zero-crossings with   regard to the system responses, which becomes a crucial key to answer this question successfully.Thereby, we present a mathematical analysis for a systematic process of identification as well as a general method for the nonlinear simultaneous identification.Even though the numerical example provided herein is a highly nonlinear system, the proposed method yields a good approximation of the exact expression in identifying the two nonlinear physical quantities in a stable and accurate manner.

𝑞 2 and q 2 Figure 2 :
Figure 2: Phase diagrams for the system responses.

2 Figure 3 :
Figure 3: Illustration of zero-crossing times in the region of steady state.

2 Figure 4 :
Figure 4: Illustration of zero-crossing times from the transient region.

Table 1 :
Identified coefficients compared with the exact ones: steady-state region.

Table 2 :
Identified coefficients compared with the exact ones: transient region.

Table 2 .
Finally (third case), we choose ( ze  )  from the transient region and ( ze acc