A new simultaneous identification of the harmonic excitations and nonlinear damping of forced damped nonlinear oscillations: a parametric approach.
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 life time. 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 . For example, Iourtchenko and Dimentberg  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 . Doyle  studied a wavelet-based deconvolution method to recover impact forces using experimentally measured responses. Zen et al.  proposed a method to obtain contact forces from experimental measurement of a single impact. Jang et al.  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 straight forward. 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 ill-posed in the sense of stability [3, 4, 8-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.
2. Forced Damped Nonlinear Oscillation
We consider a nonlinear damped dynamical system, subject to an external harmonic load [GAMMA] x cos [omega]t, of an amplitude [GAMMA] > 0 and frequency [OMEGA] > 0. The equation of motion is a second-order nonlinear ordinary differential equation for the response q
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1)
Here, the constant m > 0 represents the mass of a particle of an oscillator. The term R(q) is a general nonlinear restoring force and the term D([??]), appearing in (1), is the nonlinear damping of a system, which is assumed to have a quadratic nonlinear term
D([??]) = [c.sub.1][??] + [c.sub.2][??][absolute value of [??]], (2)
where [c.sub.1] and [c.sub.2] are the linear and nonlinear damping coefficients, respectively. The functional form, (2), of nonlinearity is often encountered in the area of engineering field .
The sum of the inertia force m[??] and restoring force R(q) and the derivative of the sum will play an important role in identifying both the excitation r and nonlinear damping D([??]) in (1). We define the sum as
J(q, [??]) = m[??] + R(q). (3)
The derivative of / with respect to time then becomes
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)
Noting that (1) is equivalent to
J = [GAMMA] x cos [omega]t - D([??]), (5)
we can express j in terms of the excitation and damping as follows:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)
by differentiating (5) withrespect to time.
In what follows, we assume that (1) has two dynamic responses of T(= 2[pi]/[omega]) -periodic limit cycle solutions [q.sub.i](t), which are forced by [[GAMMA].sub.i] for i = 1, 2
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (7)
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.
3. 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 J for i = 1, 2, denoted by [([t.sup.ze.sub.J]).sub.i], satisfying
[J.sub.i][[([t.sup.ze.sub.J]).sub.i]] = 0 for i = 1, 2. (8)
Here, [J.sub.i] stands for the sum (3) which corresponds to [q.sub.i], that is, a function from (0, [infinity]) to R
[J.sub.i] [equivalent to] m[[??].sub.i] + R([q.sub.1]). (9)
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 [([t.sup.ze.sub.j]).sub.i], i = 1, 2, is a zero-crossing time in terms of J for the system (7). Then the following holds:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)
Proof. A simple substitution of (9) into (7) yields
[J.sub.i] = [[GAMMA].sub.i] x cos [omega]t - D([[??].sub.i]). (11)
Since (11) remains valid for 0 < t < [infinity], it should also hold for the zero-crossing time of [([t.sup.ze.sub.J]).sub.i] > 0. However, the left hand side of (11) vanishes for the zero-crossing time [([t.sup.ze.sub.J]).sub.i] due to the definition (8). That is,
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)
which immediately proves (10).
Next, we introduce a zero-crossing time in terms of acceleration, ([t.sup.ze.sub.acc]) for i = 1, 2, such that
[[??].sub.i][[([t.sup.ze.sub.acc]).sub.i]] = 0 for i = 1, 2, (13)
where [[??].sub.i] 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.sub.i] defined by
[j.sub.i](t) = dJ/dt [[q.sub.i](t), [[??].sub.i](t)], for i = 1, 2, (14)
[j.sub.i] [[([t.sup.ze.sub.acc]).sub.i]] = -[[GAMMA].sub.i][omega] x sin [[[omega]([t.sup.ze.sub.acc]).sub.i]], for i = 1, 2. (15)
Here [([t.sup.ze.sub.acc]).sub.i] 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, t = [([t.sup.ze.sub.acc]).sub.i], for i = 1, 2. However, (14) is rewritten as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (16)
in which the second term of the right hand side becomes zero because the acceleration vanishes when t = [([t.sup.ze.sub.acc]).sub.i] by the definition (13). It follows that the above equation becomes (15), which completes the proof.
4. 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 [[GAMMA].sub.1], [[GAMMA].sub.2], [[GAMMA].sub.1], and [[GAMMA].sub.2] in Section 2. This results in that the developed inverse formalism also involves the parameters of [[GAMMA].sub.1], [[GAMMA].sub.2], [[GAMMA].sub.1], and [[GAMMA].sub.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 [[GAMMA].sub.1], [[GAMMA].sub.2], [[GAMMA].sub.1], and [[GAMMA].sub.2]. These can be written in a compact way as follows:
Ax = b. (17)
Here, A is an operator from [R.sup.4] to [R.sup.4] whose matrix is
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18)
and x, b [member of] [R.sup.4] are represented as (one will call them a solution vector and data vector, resp., in this paper)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)
Proof. Equations (10) and (15) in Lemmas 1 and 2 can be rewritten as the following:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (20)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (21)
-[[GAMMA].sub.1][omega] x sin [[omega][([t.sup.ze.sub.acc]).sub.1]] = [j.sub.1] [[([t.sup.ze.sub.acc]).sub.1]], (22)
-[[GAMMA.sub.]2][omega] x sin [[omega][([t.sup.ze.sub.acc]).sub.2]] = [j.sub.2][[([t.sup.ze.sub.acc]).sub.2]] (23)
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.
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
x = [A.sup.-1] b, (24)
where the notation A-1 denotes the inverse operator for the linear operator A 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.
Theorem 4 (excitations). The excitations T for i=1, 2 in (7) are identified as an analytical form
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)
provided that sin [[omega][([t.sup.ze.sub.acc]).sub.i]] = 0 for i = 1, 2.
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 -[omega] x sin [[omega][([t.sup.ze.sub.acc]).sub.1]] to obtain [[GAMMA].sub.1] as follows:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (26)
Similarly, we have
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (27)
from (23). This completes the proof.
The formulas of the damping coefficients [c.sub.1] and [c.sub.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 [c.sub.1] and the nonlinear damping coefficient [c.sub.2] in (2) are determined as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (28)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (29)
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 [c.sub.1] of (28) and then [c.sub.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 A in (18) should satisfy det A [not equal to] 0, which is to be discussed in the following theorem.
Theorem 7 (uniqueness). Let one denote by [DELTA],
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (30)
Then, the excitations and nonlinear damping in (7) can be identified in a unique way whenever [DELTA] [not equal to] 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 [DELTA] [not equal 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.
6. 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 R(q) in (1) involves the cubic nonlinearity, that is,
R(q) = [k.sub.1]q + [k.sub.2][q.sup.3], (31)
where the values for [k.sub.1] and [k.sub.2] are chosen as 40 and 10, respectively. Then, the responses for the system of (1) can be observed as following subsection.
6.1. System Responses. With the initial conditions imposed as
[q.sub.i](0) = 0, [[??].sub.i](0) = 0 for i = 1, 2, (32)
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 [c.sub.1] and [c.sub.2] are set to be 4 and 2, respectively. The excitation frequency co is taken as 4 with [[GAMMA].sub.1] = 1 and [[GAMMA].sub.2] = 2. The calculated responses of displacement and velocity of the system are depicted in Figure 1 and their phase diagrams in Figure 2.
6.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, [([t.sup.ze.sub.J]).sub.i] and [([t.sup.ze.sub.acc]).sub.i] 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 equation solver, the results are calculated as [([t.sup.ze.sub.J]).sub.1] = 8.401 sec, [([t.sup.ze.sub.J]).sub.2] = 8.408 sec, [(t.sup.ze.sub.acc).sub.1] = 9.182 sec, and [([t.sup.ze.sub.acc]).sub.2] = 9.186 sec, each of which is indicated in Figure 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, time-simulated 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 listed in Table 2. Finally (third case), we choose [([t.sup.ze.sub.J]).sub.i] from the transient region and [([t.sup.ze.sub.acc]).sub.i] from the steady-state region for each i = 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 A in (18) are of O([10.sup.2]), meaning that our system equation (17) is well conditioned.
7. 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.
This research is supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (Grant no. 2011-0010090). In addition, this research was also supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (Grant no. K20903002030-11E0100-04610). Finally, the author would like to extend his special thanks to the students, Miss. Hyeree Bae and Mr. Jinsoo Park in Naval Architecture and Ocean Engineering Department, Pusan National University for their help in the paper.
 M. Kulisiewicz, "A nonparametric method of identification of vibration damping in non-linear dynamic systems, " International Journal of Solids and Structures, vol. 19, pp. 601-609, 1983.
 D. V. Iourtchenko and M. F. Dimentberg, "In-service Identification of non-linear damping from measured random vibration, " Journal ofSound and Vibration, vol.255, pp.549-554, 2002.
 T. S. Jang, H. S. Choi, and S. L. Han, "A new method for detecting non-linear damping and restoring forces in nonlinear oscillation systems from transient data, " International Journal of Non-Linear Mechanics, vol.44, pp.801-808, 2009.
 T. S. Jang, "Non-parametric simultaneous identification of both the nonlinear damping and restoring characteristics of nonlinear systems whose dampings depend on velocity alone, " Mechanical Systems and Signal Processing, vol. 25, pp. 1159-1173, 2011.
 S. F. Masri and T. K. Caughey, "A nonparametric identification technique for nonlinear dynamic problems, " ASME Journal of Applied Mechanics, vol. 46, pp. 433-447, 1979.
 J. Doyle, "A wavelet deconvolution method for impact force identification, " Experimental Mechanics, vol.37, pp.403-408, 1997.
 Q. H. Zen, M. Chapin, and D. B. Bogy, "A force identification method for slider/disk contact force measurement, " IEEE Transactions on Magnetics, vol. 36, pp. 2667-2670, 2000.
 T. S. Jang, H. S. Baek, H. S. Choi, and S. G. Lee, "Anew method for measuring nonharmonic periodic excitation forces in nonlinear damped systems, " Mechanical Systems and Signal Processing, vol. 256, pp. 2219-2228, 2011.
 T. S. Jang, H. S. Back, S. L. Han, and T. Kinoshita, "Indirect measurement of the impulsive load to a nonlinear system from dynamic responses: inverse problem formulation, " Mechanical Systems and Signal Processing, vol.24, pp.1665-1681, 2010.
 T. S. Jang, S. L. Han, and T. Kinoshita, "An inverse measurement of the sudden underwater movement of the sea-floor by using the time-history record of the water-wave elevation, " Wave Motion, vol.47, no.3, pp.146-155, 2010.
 M. Taylan, "The effect of nonlinear damping and restoring in ship rolling, " Ocean Engineering, vol.27, pp.921-932, 2000.
T. S. Jang
Department of Naval Architecture and Ocean Engineering, Pusan National University, Busan 609-735, Republic of Korea
Correspondence should be addressed to T. S. Jang; firstname.lastname@example.org
Received 27 September 2012; Revised 21 December 2012; Accepted 23 December 2012
Academic Editor: George Jaiani
TABLE 1: Identified coefficients compared with the exact ones: steady-state region. [c.sub.1] [c.sub.2] [[GAMMA]. [[GAMMA]. sub.1] sub.2] Exact 4 2 1 2 Identified 3.9761 2.0195 0.9988 2.0003 Error (%) 0.6 1.0 0.1 0.02 TABLE 2: Identified coefficients compared with the exact ones: transient region. [c.sub.1] [c.sub.2] [[GAMMA]. [[GAMMA]. sub.1] sub.2] Exact 4 2 1 2 Identified 3.9120 2.4773 0.9984 2.0227 Error (%) 2.2 23.87 0.16 1.14 TABLE 3: Identified coefficients compared with the exact ones: [([t.sup.ze.sub.J]).sub.i]--the transient region, [([t.sup.ze.sub.acc]).sub.i]--the steady-state region (i = 1,2). [c.sub.1] [c.sub.2] [[GAMMA]. [[GAMMA]. sub.1] sub.2] Exact 4 2 1 2 Identified 3.9702 2.0948 0.9988 2.0003 Error (%) 0.75 4.74 0.1 0.02
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Publication:||Journal of Applied Mathematics|
|Date:||Jan 1, 2013|
|Previous Article:||Electroencephalography signal grouping and feature classification using harmony search for BCI.|
|Next Article:||Cost-sensitive feature selection of numeric data with measurement errors.|