# Normal form for high-dimensional nonlinear system and its application to a viscoelastic moving belt.

1. Introduction

Bifurcation and stability analysis of nonlinear differential equations is one of the challenging problems of mathematicians and engineers. Normal form theory is one of the most important tools for such analysis [1-3]. The normal form theory for differential equations can be dated back to the pioneer works of the renowned mathematician Poincare [4]. He tried to use some change of variables to alter nonlinear systems into linear ones. The idea of the method is to simplify the system such that the topological behavior of the system in the vicinity of a singularity point remains unchanged.

Some recent developments of the theory of normal form can be found in [5-15]. Chen et al. [5] presented the renormalization group theory, which was used for the search of normal forms for large classes of finite-dimensional vector fields [6-8]. Strozyna and Zondek [9, 10] made use of the time rescaling to achieve the results on the orbital equivalence of vector fields. In addition, Dullin and Meiss [11] and Murdock [12] dealt with the problems for the normal forms of nilpotent systems, whose linear part about the origin is a nilpotent matrix. Benderesky and Churchill [13, 14] and Sanders [15] studied the spectral sequences for the normal forms of vectors. Such spectral sequences can provide valuable information on the normal forms.

On the other hand, Kuznetsov [16] considered the normal form theory in an application-oriented way for computation of high-dimensional nonlinear systems. Zhang and his coworkers [17,18] employed the adjoint operator method to obtain the higher-order normal forms of high-dimensional nonlinear dynamical systems and the associated nonlinear transformations. Zhang and Leung [19] considered a general four-dimensional normal form of a double Hopf bifurcation. Yu and his associates [20-22] developed efficient computing methods for parametric normal forms. They also applied the new method to consider controlling bifurcations of the nonlinear dynamical systems. Chen and Dora [23] were devoted to the development of effective methods for further reductions of the classical normal forms for complex dynamical systems. From a practical point of view, a better understanding and knowledge of the normal forms of various complex nonlinear systems will further promote the potential interest for the analysis of real engineering problems.

The normal form theory plays an important role in the study of bifurcation behavior of differential dynamical systems. Itovich and Moiola [24] made use of the frequency domain and the normal form methodologies to analyze the unfolding of a nonresonant double Hopf singularity. Zhang et al. [25] employed the center manifold reduction and normal form method to obtain the singular bifurcation of a ring of three coupled advertising oscillators with delay. Jiang and Yuan [26] studied the classical Van Der Pol equation with Bogdanov-Takens singularity and bifurcation. Gattulli et al. [27] analyzed the postcritical behavior of a single degree of freedom system equipped with a Tuned Mass Damper for double Hopf bifurcation in the neighbourhood of 1: 1 resonance. Using the normal form method and the center manifold theory, Li and his associates [28] investigated the double Hopf bifurcation of the trivial equilibrium for delay-coupled limit cycle oscillators. Buono and Belair [29] studied the normal form of a vector field, which is generated by a scalar delay-differential equations at nonresonant double Hopf bifurcation points.

Some of the above works use the Jordan canonical form of the leading matrix A. However, it is well known that handling the eigenvalues and Jordan canonical forms is very difficult in computer algebra system. In this paper, a new computation method by direct computation is developed to refine the normal forms for high-dimensional nonlinear systems. We do not need to compute the Jordan canonical form of A nor its eigenvalues. Our method is applicable in both the nilpotent and the nonnilpotent cases.

In this paper, we will develop an efficient method for computing the normal forms directly for general four-dimension systems and apply the method to consider controlling bifurcations. The approach is efficient since it does not require the computation of the Jordan canonical form of A or its eigenvalues. Besides, the proposed method is applied to investigate the nonlinear oscillations of a viscoelastic moving belt under parametric excitations. The rest of the paper is organized as follows. In Section 2, the essential idea behind the method in [23] is briefly introduced. The new computation method is described in detail by analyzing the four-dimensional nonlinear dynamical systems in Section 3. The applications to stability and bifurcation analysis on the viscoelastic moving belt are presented in Sections 4 and 5 to show the efficiency of the method. Finally, conclusions are drawn in Section 6.

2. Normal Forms for Nonlinear System

Consider a dynamical system described by the following differential equation:

[??] = F(x) = Ax + [[infinity].summation over (k=2)] [f.sup.k] (x), x [member of] [R.sup.n], (1)

where Ax represents the linear part, A is the Jacobian matrix, and [f.sup.k](x) denotes the kth-order vector homogeneous polynomials of x.

Without loss of generality, A is expressed in terms of the standard Jordan canonical form. Note that system (1) is assumed to have an equilibrium at the origin x = 0.

We take the coordinate transformation as follows:

x = y + [[phi].sup.k] (y). (2)

Substituting (2) into (1) gives

[??] = [(I + [[partial derivative].sub.y] < [[phi].sup.k]).sup.-1] F(y + < [[phi].sup.k] (y)), (3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)

where [[partial derivative].sub.y][[phi].sup.k] is the Jacobian matrix of [[phi].sup.k] with respect to y.

Then, (4) is substituted back into (3) to form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

Employing the normal form theory, we introduce a linear operator as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

Hence, we can write

[H.sup.n.sub.k] = [R.sup.k] + [C.sup.k], (7)

where [R.sup.k] represents the range of [L.sup.k.sub.A] and [C.sup.k] is the complementary space to [R.sup.k].

Hence, the kth-order terms can be simplified to

[g.sup.k](y) [member of] [C.sup.k]. (8)

The rationale for the classical normal form theory can be explained by the following theorem (see [30]).

Theorem 1. Let the notations be the same as above. Suppose that the decomposition (7) is given for k = 2, ..., N. Then, there exists a sequence of near identity changes of variables x = y + [[phi].sup.k](y), in which [[phi].sup.k](y) [member of] [H.sup.n.sub.k] k. Therefore, the dynamical system (1) is transformed into

[??] = Ay + [g.sup.2] (y) + *** + [g.sup.N] (y) + O([parallel]y[[parallel].sup.N+1]), (9)

where [g.sup.k](y) [member of] [C.sup.k] for k = 2, ..., N.

By applying the Takens normal form theory [30], one arrives at the kth-order normal form [g.sup.k](y), while those parts belonging to [R.sup.k] can be removed by properly choosing the coefficients of the nonlinear transformation [[phi].sup.k](y).

3. Computation of Normal Forms and Their Coefficients

Consider a four-dimensional generalized averaged system with [Z.sub.2][cross product][Z.sub.2]-symmetry governed by

[??] = F(x) = Ax + [f.sup.3] (x), x [member of] [R.sup.4], (10)

where [f.sup.3](x) [member of] [H.sup.3.sub.4]; that is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

and [absolute value of m] = [m.sub.1] + [m.sub.2] + [m.sub.3] + [m.sub.4].

At the same time, (6) becomes

[[partial derivative].sub.y] [[phi].sup.3] (y)Ay - A[[phi].sup.3] (y) = [f.sup.3] - [g.sup.3], (12)

The problem at hand is to determine [phi], so that [g.sup.3] contains the smallest possible number of monomials.

Consider the following three cases of the Jordan matrix A in four-dimensional nonlinear systems.

(i) The Jordan matrix A has two pairs of pure imaginary eigenvalues.

(ii) The Jordan matrix A has one nonsemisimple double zero and a pair of pure imaginary eigenvalues.

(iii) The Jordan matrix A has two nonsemisimple double zero eigenvalues.

The forms of the Jordan matrix A in these cases (i)-(iii) can be, respectively, represented by (13)-(15) as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (13)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (14)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

It is easy to see that the three-order polynomial solutions in four variables can be obtained from (12). To achieve this, we write

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (16)

with [a.sub.q], [b.sub.q], [c.sub.q], [d.sub.q], [[alpha]'.sub.q], [[beta]'.sub.q], [[gamma]'.sub.q] and [[eta]'.sub.q] to be determined later.

For the case (i), (12) can be expressed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17d)

where [e.sub.1] = (1, 0, 0, 0), [e.sub.2] = (0, 1, 0, 0), [e.sub.3] = (0, 0, 1, 0), and [e.sub.4] = (0, 0, 0, 1).

Then, we solve the following equations:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18d)

For the case (ii), we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (19a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (19b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (19c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (19d)

and we need to solve the following equations:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20d)

Finally, the case (iii) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21d)

Similarly, the following equations must be solved:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22d)

Let

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

Substituting (23) into (18a)-(18d), (20a)-(20d), and (22a)-(22d), separately, the following three sets of 3-order nonlinear algebraic equations are resulted.

(i) For the case of two pairs of pure imaginary eigenvalues, we arrive at

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

(ii) For the case of one nonsemisimple double zero and a pair of pure imaginary eigenvalues, we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)

(iii) For the case of two nonsemisimple double zero eigenvalues, we derive

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (26)

It should be clear that [M.sub.q] = 0if one of the components of q is completely negative. We choose the lexicographical order on the set {q = ([q.sub.1], [q.sub.2], [q.sub.3], [q.sub.4]) [member of] [N.sup.4]: [absolute value of q] = 3}.

Let [q.sub.1] = i, [q.sub.2] = j, [q.sub.3] = l, and [q.sub.4] = 3 - i - j - l. Then, we have

[Z.sub.i,j,l] = [M.sub.(i,j,l,3-i-j-l)], [[LAMBDA].sub.i,j,l] = [f.sup.3], [[LAMBDA]'.sub.i,j,l] = [g.sup.3]. (27)

Substituting (27) into (24)-(26) gives the following three conclusions.

(i) For the case of two pairs of pure imaginary eigenvalues, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (28)

Let [Z.sub.0,j,l] (0 [less than or equal to] j [less than or equal to] 3, 0 [less than or equal to] l [less than or equal to] 3) be given arbitrarily. We can determine [Z.sub.i+1,j-1,l] by the above formulae to achieve [[LAMBDA]'.sub.i,j,l] = 0 for 0 [less than or equal to] i [less than or equal to] k - 1, 1 [less than or equal to] j [less than or equal to] k, and 0 [less than or equal to] l [less than or equal to] k. In fact, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (29)

The remaining equation is (for j = 0)

[[LAMBDA].sub.i,0,l] = [[LAMBDA].sub.i,0,l] - [w.sub.1][Z.sub.i-1,1,l] + [w.sub.2](l + 1)[Z.sub.i,0,l+1] -[w.sub.2](4 - i - l)[Z.sub.i,0,l-1] + A[Z.sub.i,0,l]. (30)

(ii) For the case of one nonsemisimple double zero and a pair of pure imaginary eigenvalues, (25) becomes

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (31)

In view of (31), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (32)

Hence, the remaining equation is given by (for j = 0)

[[LAMBDA].sub.i,0,l] = [[LAMBDA].sub.i,0,l] + w (l + 1) [Z.sub.i,0,l+1] -w (4 - i - l) [Z.sub.i,0,l-1] + A[Z.sub.i,0,l]. (33)

(iii) For the case of two nonsemisimple double zero eigenvalues, we deduce

(i+1)[Z.sub.i+1,j-1,l] + (l + 1)[Z.sub.i,j,l+1] - A[Z.sub.i,j,l] = [[LAMBDA].sub.i,j,l] - [[LAMBDA]'.sub.i,j,l]. (34)

Similarly, making use of (34) leads to

(i + 1)[Z.sub.i+1,j-1,l] = [[LAMBDA].sub.i,j,l] - (l + 1)[Z.sub.i,j,l+1] + A[Z.sub.i,j,l]. (35)

From (34) and (35), the remaining equation can be written as (for j = 0)

[[LAMBDA]'.sub.i,0,l] = [[LAMBDA].sub.i,0,l] - (l + 1)[Z.sub.i,0,l+1] + A[Z.sub.i,0,l]. (36)

We determine all [Z.sub.i,j,l] for i [greater than or equal to] 1 as functions of [Z.sub.0,j,l], which may be solved in some cases by making some of the terms [[LAMBDA]'.sub.i,0,l] equal zeroes.

4. Application to a Viscoelastic Moving Belt Model

In this section, we apply the proposed method in Section 3 to a parametrically excited viscoelastic moving belt with the external damping. Consider the viscoelastic moving belt model [31] with cross-sectional area A, length L between two end supports, axial velocity v, and viscous damping coefficient ?? as shown in Figure 1. A Cartesian coordinate system (Oxyz) is adopted, which is located in the plane of the viscoelastic moving belt. Another coordinate system is a moving coordinate fixed on the belt. The u and w denote the displacements in the x and y directions, respectively. It is assumed that the tension T is characterized as a small periodic perturbation [T.sub.1] cos [OMEGA]t on the steady-state tension [T.sub.0]; that is,T = [T.sub.0] + [T.sub.1] cos [OMEGA]t. Since the belt tension is assumed to dominate the transverse stiffness, the bending stiffness of the viscoelastic moving belt is neglected. The equations of motion for the transverse vibration of the belt are based on an axially moving string model. The nondimensional nonlinear governing equation of motion for the viscoelastic moving belt under parametric excitations can be written as follows [31-33]:

[w,.sub.tt] + 2[gamma][w.sub.,tx] + ([[gamma].sup.2] - 1 - a cos [omega]t) [w.sub.,xx] + 2[mu][w.sub.,t] - N(w) = 0, (37)

where the comma subscript denotes the partial differentiation, and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (38)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (39)

The boundary conditions are imposed by

w(0, t) = 0, w(L, t) = 0. (40)

In the subsequent analysis, we use the method of multiple scales and Galerkin's approach in the partial differential governing equations of the viscoelastic moving belt. We introduce the mass, gyroscopic, and linear stiffness operators as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (41)

Substituting (41) into (37) leads to the standard symbolic form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (42)

To obtain a system that is suitable for the application of the method of multiple scales, we introduce the scale transformations

a [right arrow] [epsilon]a, N [right arrow] [epsilon]N, [mu] [right arrow] [epsilon][mu], (43)

where [epsilon] is the small perturbation parameter.

Substituting (43) into (42), we obtain the following dimensionless nonlinear system under parametric excitations

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (44)

The method of multiple scales can now be applied to search for the uniform solutions of (44) in the following form:

w(x, t, [epsilon]) = [w.sub.0] (x, [T.sub.0], [T.sub.1]) + [epsilon][w.sub.1]) + (x, [T.sub.0], [T.sub.1]) + ***, (45)

where [T.sub.0] = t, [T.sub.1] = [epsilon]t.

The differential operators of the method of multiple scales can be defined as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (46a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (46b)

where [D.sub.n] = [partial derivative]/[partial derivative][T.sub.n], n = 0, 1,....

In this paper, we investigate the case of primary parametric resonance for the nth and lth order modes of (44), and we introduce the following linear transformation:

[x.sub.1] [right arrow] [x.sub.2], [x.sub.2] [right arrow] [x.sub.1], [x.sub.3] [right arrow] [x.sub.4], [x.sub.4] [right arrow] [x.sub.3]. (47)

Then, the four-dimensional averaged equations in the Cartesian form can be obtained as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (48a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (48b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (48c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (48d)

where the coefficients presented in (48a)-(48d) are listed in the appendix.

The above algorithm applied to the system (48a)-(48d) leads to the following normal form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (49a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (49b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (49c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (49d)

Comparing the method developed here with other methods given in [17, 31], it is observed that normal form of the averaged (49a)-(49d) is simpler than normal forms obtained in [17, 31].

5. Stability and Bifurcation Analysis on the Viscoelastic Moving Belt

It is known that (49a)-(49d) has a trivial zero solution ([x.sub.1], [x.sub.2], [x.sub.3], [x.sub.4]) = (0, 0, 0, 0) in which the Jacobi matrix can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (50)

The characteristic equation corresponding to the trivial zero solution is

([[lambda].sup.2] + 2[mu][lambda] + [[mu].sup.2] + [[sigma].sup.2.sub.1] - [f.sup.2]) ([[lambda].sup.2] + 2[mu][lambda] + [[mu].sup.2] + [[sigma].sup.2.sub.2] - [D.sup.2]) = 0. (51)

The eigenvalues of the above equations are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (52)

Take a linear transformation of coordinate as follows:

{X} = [D]{Z}, (53)

where [D] is the matrix of the eigenvectors of the linear part in (49a)-(49d),

{X} = [[[x.sub.1], [x.sub.2], [x.sub.3], [x.sub.4]].sup.T], {Z} = [[[z.sub.1], [[bar.z].sub.1], [z.sub.2], [[bar.z].sub.2]].sup.T]. (54)

After substituting the above transformation into (49a)-(49d), one obtains the equations in new complex coordinateas follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (55a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (55b)

Let

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (56)

In polar coordinates ([r.sub.1], [r.sub.2], [[phi].sub.1], [[phi].sub.2]), systems (55a) and (55b) can be written as

[[??].sub.1] = [r.sub.1] (-u +4[a.sub.e2v1] [[sigma].sub.1]/[[sigma].sub.1] + f [r.sup.2.sub.1] + 4[b.sub.e2v1] [[sigma].sub.2]/[[sigma].sub.2] + D[r.sup.2.sub.2]), (57a)

[[??].sub.2] = [r.sub.2] (-u +4[B.sub.e2v1] [[sigma].sub.1]/[[sigma].sub.1] + f [r.sup.2.sub.1] + 4[A.sub.e2v1] [[sigma].sub.2]/[[sigma].sub.2] +D [r.sup.2.sub.2]), (57b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (57c)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (57d)

Since our analysis is local, we can truncate the higher order and consider the system

[[??].sub.1] = [r.sub.1] (-u +[4[a.sub.e2V1] [[sigma].sub.1]/[[sigma].sub.1] + f] [r.sup.2.sub.1] + 4[b.sub.e2V1] [[sigma].sub.2] [[sigma].sub.2] +D [r.sup.2.sub.2]), (58a)

[[??].sub.2] = [r.sub.2] (-u + 44[B.sub.e2V1] [[sigma].sub.1]/[[sigma].sub.1] + f] [r.sup.2.sub.1]+ 4[A.sub.e2V1] [[sigma].sub.2] [[sigma].sub.2] +D [r.sup.2.sub.2]), (58b)

[[??].sub.1] = [square root of ([[sigma].sup.2.sub.1] - [f.sup.2])], (58c)

[[??].sub.2] = [square root of ([[sigma].sup.2.sub.2] -[D.sup.2])]. (58d)

The first two equations are independent of the last two. The last two equations describe rotations in the planes [r.sub.2] = 0 and [r.sub.1] = 0 with angular velocities [square root of ([[sigma].sup.2.sub.1] - [f.sup.2])] and [square root of ([[sigma].sup.2.sub.2] -[D.sup.2])], respectively. This does not change the bifurcation diagrams. So, it is enough to study the first two equations. Therefore, we can study the original four-dimensional system (49a)-(49d) by analyzing the planar system

[[??].sub.1] = [r.sub.1] (-u + 4[a.sub.e2v1] [[sigma].sub.1]/[[sigma].sub.1] + f [r.sup.2.sub.1] + 4[b.sub.e2v1][[sigma].sub.2]/[[sigma].sub.2] +D [r.sup.2.sub.2]), (59a)

[[??].sub.2] = [r.sub.2] (-u + 4[B.sub.e2v1] [[sigma].sub.1]/[[sigma].sub.1] + f [r.sup.2.sub.1] + 4[A.sub.e2v1] [[sigma].sub.2]/[[sigma].sub.2] +D [r.sup.2.sub.2]). (59b)

This system is called the amplitude system. The trivial equilibrium, [r.sub.1] = 0, [r.sub.2] = 0, corresponds to the trivial equilibrium of the original system. The study of the amplitude system is simplified if we use squares [[rho].sub.1,2] of the amplitudes as follows:

[[rho].sub.k] = [r.sup.2.sub.k], k = 1,2. (60)

[[??].sub.1] = 2[[rho].sub.1] (-u + 4[a.sub.e2v1] [[sigma].sub.1]/[[sigma].sub.1] + f [[rho].sub.1] + 4[b.sub.e2v1][[sigma].sub.2]/[[sigma].sub.2] +D [[rho].sub.2]), (61a)

[[??].sub.2] = 2[[rho].sub.2] (-u + 4[B.sub.e2v1] [[sigma].sub.1]/[[sigma].sub.1] + f [[rho].sub.1]+ 4[A.sub.e2v1][[sigma].sub.2]/[[sigma].sub.2] + D [[rho].sub.2]). (61b)

The behavior of systems (61a) and (61b) depends on the coefficients of [[rho].sub.1] and [[rho].sub.2]. We start our analysis with the case when 4[a.sub.e2v1][[sigma].sub.1]/([[sigma].sub.1] + f) < 0 and 4[A.sub.e2v1] [[sigma].sub.2]/([[sigma].sub.2] + D) < 0. The case 4[a.sub.e2v1] [[sigma].sub.1]/([[sigma].sub.1] + f) > 0 and 4[A.sub.e2v1] [[sigma].sub.2]/([[sigma].sub.2] + D) > 0 can be reduced to the previous one by time reversal. The choices 4[a.sub.e2v1] [[sigma].sub.1]/([[sigma].sub.1] + f) < 0 and 4[A.sub.e2v1] [[sigma].sub.2]/([[sigma].sub.2]+ D) < 0 imply that both of the primary Hopf bifurcations are supercritical and stable.

For the case of 16[A.sub.e2v1] [a.sub.e2v1] [[sigma].sub.1] [[sigma].sub.2]/([[sigma].sub.1] + f)([[sigma].sub.2]+ D) < 0, First we can reduce the number of the coefficients in (61a) and (61b) by rescaling. Let

[[xi].sub.1] = 4[a.sub.e2v1][[sigma].sub.1]/[[sigma].sub.1] + f [[rho].sub.1], [[xi].sub.2] = 4[a.sub.e2v1][[sigma].sub.2]/[[sigma].sub.2] + D [[rho].sub.2], [tau] = 2t, (62)

We obtain the system

[[??].sub.1] = [[xi].sub.1] (-u + [[xi].sub.1] - [b.sub.e2v1]/[A.sub.e2v1] [[xi].sub.2]), (63a)

[[??].sub.2] = [[xi].sub.2] (-u + 4[B.sub.e2v1]/[a.sub.e2v1] [[xi].sub.1] - [[xi].sub.2]). (63b)

The trivial and the nontrivial equilibria of (63a) and(63b) have the representations

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (64)

The nontrivial equilibrium [E.sub.3] can bifurcate. Suppose that 4[a.sub.e2v1] [[sigma].sub.1]/([[sigma].sub.1] + f) > 0 and 4[A.sub.e2v1] [[sigma].sub.2]/([[sigma].sub.2] + D) < 0. The opposite case can be treated similarly. The Hopf bifurcation and consequent existence of cycles are only possible in following three cases: (1) [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

The numerical simulation result is given in Figure 2, which is the trajectory from the 4-dimensional space ([x.sub.1], [x.sub.2], [x.sub.3], [x.sub.4]) into the 3-dimensional space ([x.sub.1], [x.sub.2], [x.sub.3]). Figure 2 shows the trajectory of the systems (61a) and (61b) for [mu] = 0.0001, [[sigma].sub.1] = 2, f = -124, [a.sub.e2v1] = 56, [a.sub.e1v2] = 2.3, [b.sub.e2v1] = 15, [b.sub.e1V2] = 19.6, [[sigma].sub.2] = 4.5, [A.sub.e2V1] = -127, [A.sub.e1V2] = 12, 4[B.sub.e2v1] = 1.82, [B.sub.e1v2] = 12, and D = 190.

6. Conclusions

An efficient method for computing the normal form of high-dimensional nonlinear systems is presented in this paper. This computation method is applied to obtain the normal form of the averaged equation for the viscoelastic moving belt under parametric excitations. Based on the current studies, it is found that the newly developed computation method improves the classical normal form. Therefore, it is a further reduction of the classical normal form. Meanwhile, the normal form derived herein is used to explore the bifurcation and stability analysis of axially viscoelastic moving belts under parametric excitations.

In contrast to earlier works, this article argues that the normal forms of high-dimensional nonlinear systems may always be achieved without computing either the Jordan canonical form of A or its eigenvalues. This significantly reduces the difficulties which other computing methods may face in obtaining normal forms. Therefore, it is more convenient to utilize the approach developed here to compute the normal forms of averaged equations for different resonant cases. It is also found that the normal form of the averaged equation by using the new method is simpler than those obtained in [17, 31]. This is very useful for the complex behavior patterns analysis of the nonlinear dynamical systems in an abstract sense. Therefore, this has opened the area for more research works.

Appendix

The coefficients in the averaged system (48a)-(48d) are presented as follows,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.1)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.2)

where a, [E.sub.e], [E.sub.v], [[omega].sub.l] and [[omega].sub.n] are given in (38) and (48a)-(48d). The terms [m.sub.ij] (i = 0, 2, 3; j = l, n) are complex eigenfunctions of the displacement field, the detailed derivations of [m.sub.ij] (i = 0, 2, 3; j = l, n) in (A.1) and (A.2) can be referred to [25, 27].

http://dx.doi.org/10.1155/2014/879564

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this article.

Acknowledgments

The authors gratefully acknowledge the support of the National Natural Science Foundations of China (NNSFC), through Grant nos. 11302184 and 11202189, the Scientific Research Fund of Zhejiang Provincial Education Department, through Grant no. Y201121157, and the Scientific Research Foundation of Xiamen University of Technology, through Grant no. 90030631. We are also grateful for helpful comments which were provided by Wei Zhang at the Beijing University of Technology.

References

[1] J. Murdock, Normal Forms and Unfoldings for Local Dynamical Systems, Springer Monographs in Mathematics, Springer, New York, NY, USA, 2003.

[2] Y. A. Kuznetsov, Elements of Applied Bifurcation Theory, vol. 112 of Applied Mathematical Sciences, Springer, New York, NY, USA, 2004.

[3] A. Rincon, F. Angulo, and G. Olivar, "Control of an anaerobic digester through normal form of fold bifurcation," Journal of Process Control, vol. 19, no. 8, pp. 1355-1367, 2009.

[4] H. Poincare, Sur les proprietes des fonctions definies par des equations aux differences partielles these inaugural, Gauthier-Villars, Paris, France, 1879.

[5] L. Y. Chen, N. Goldenfeld, and Y. Oono, "Renormalization group theory for global asymptotic analysis," Physical Review Letters, vol. 73, no. 10, pp. 1311-1315, 1994.

[6] R. E. L. DeVille, A. Harkin, M. Holzer, K. Josic, and T. J. Kaper, "Analysis of a renormalization group method and normal form theory for perturbed ordinary differential equations," PhysicaD: Nonlinear Phenomena, vol. 237, no. 8, pp. 1029-1052, 2008.

[7] H. Chiba, "Simplified renormalization group method for ordinary differential equations," Journal of Differential Equations, vol. 246, no. 5, pp. 1991-2019, 2009.

[8] M. Holzer, Renormalization group methods for singularly perturbed systems, normal forms and stability of traveling waves in a reaction-diffusion-mechanics systems [Ph.D. thesis], Boston University, 2010.

[9] E. Strozyna and H. Zoladek, "The analytic and formal normal form for the nilpotent singularity," Journal of Differential Equations, vol. 179, no. 2, pp. 479-537, 2002.

[10] E. Strozyna and H. Zoladek, "Orbital formal normal forms for general Bogdanov-Takens singularity," Journal of Differential Equations, vol. 193, no. 1, pp. 239-259, 2003.

[11] H. R. Dullin and J. D. Meiss, "Nilpotent normal form for divergence-free vector fields and volume-preserving maps," Physica D: Nonlinear Phenomena, vol. 237, no. 2, pp. 156-166, 2008.

[12] J. Murdock, "On the structure of nilpotent normal form modules," Journal of Differential Equations, vol. 180, no. 1, pp. 198-237, 2002.

[13] M. Benderesky and R. Churchill, "A spectral sequence approach to normal forms," in Recent Developments in Algebraic Topology, vol. 407 of Contemporary Mathematics, pp. 27-81, 2006.

[14] M. Bendersky and R. C. Churchill, "Normal forms in a cyclically graded Lie algebra," Journal of Symbolic Computation, vol. 41, no. 6, pp. 633-662, 2006.

[15] J. A. Sanders, "Normal form theory and spectral sequences," Journal of Differential Equations, vol. 192, no. 2, pp. 536-552, 2003.

[16] Y. A. Kuznetsov, "Practical computation of normal forms on center manifolds at degenerate Bogdanov-Takens bifurcations," International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 15, no. 11, pp. 3535-3546, 2005.

[17] W. Zhang, F. X. Wang, and J. W. Zu, "Computation of normal forms for high dimensional non-linear systems and application to non-planar non-linear oscillations of a cantilever beam," Journal of Sound and Vibration, vol. 278, no. 4-5, pp. 949-974, 2004.

[18] W. Zhang, Y. Chen, and D. Cao, "Computation of normal forms for eight-dimensional nonlinear dynamical system and application to a viscoelastic moving belt," International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 35-58, 2006.

[19] Q. C. Zhang and A. Y. T. Leung, "Normal form of double Hopf bifurcation in forced oscillators," Journal of Sound and Vibration, vol. 231, no. 4, pp. 1057-1069, 2000.

[20] P. Yu and G. Chen, "The simplest parametrized normal forms of HOPf and generalized HOPf bifurcations," Nonlinear Dynamics, vol. 50, no. 1-2, pp. 297-313, 2007. Abstract and Applied Analysis 11

[21] P. Yu and A. Y. T. Leung, "Normal forms of vector fields with perturbation parameters and their application," Chaos, Solitons & Fractals, vol. 34, no. 2, pp. 564-579, 2007.

[22] M. Gazor and P. Yu, "Infinite order parametric normal form of Hopf singularity," International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 18, no. 11, pp. 3393-3408, 2008.

[23] G. Chen and J. D. Dora, "Further reductions of normal forms for dynamical systems," Journal of Differential Equations, vol. 166, no. 1, pp. 79-106, 2000.

[24] G. R. Itovich and J. L. Moiola, "Non-resonant double Hopf bifurcations: the complex case," Journal of Sound and Vibration, vol. 322, no. 1-2, pp. 358-380, 2009.

[25] C. Zhang, H. Yin, and H. Zheng, "Simple bifurcation of coupled advertising oscillators with delay," Applied Mathematics Letters, vol. 24, no. 11, pp. 1840-1844, 2011.

[26] W. Jiang and Y. Yuan, "Bogdanov-Takens singularity in van der Pol's oscillator with delayed feedback," Physica D: Nonlinear Phenomena, vol. 227, no. 2, pp. 149-161, 2007.

[27] V. Gattulli, F. Di Fabio, and A. Luongo, "One to one resonant double Hopf bifurcation in aeroelastic oscillators with tuned mass dampers," Journal of Sound and Vibration, vol. 262, no. 2, pp. 201-217, 2003.

[28] Y. Li, W. Jiang, and H. Wang, "Double Hopf bifurcation and quasi-periodic attractors in delay-coupled limit cycle oscillators," Journal of Mathematical Analysis and Applications, vol. 387, no. 2, pp. 1114-1126, 2012.

[29] P. Buono and J. Belair, "Restrictions and unfolding of double Hopf bifurcation in functional differential equations," Journal of Differential Equations, vol. 189, no. 1, pp. 234-266, 2003.

[30] F. Takens, "Singularities of vector fields," Publications Mathematiques de l'IHES, vol. 43, pp. 47-100, 1974.

[31] W. Zhang and M. H. Yao, "Multi-pulse orbits and chaotic dynamics in motion of parametrically excited viscoelastic moving belt," Chaos, Solitons and Fractals, vol. 28, no. 1, pp. 42-66, 2006.

[32] L. H. Chen, W. Zhang, and Y. Q. Liu, "Modeling of nonlinear oscillations for viscoelastic moving belt using generalized Hamilton's principle," Journal of Vibration and Acoustics, Transactions of the ASME, vol. 129, no. 1, pp. 128-132, 2007.

[33] C. Z. Song, Studies on non-linear dynamics of a multi-degree of freedom viscoelastic drive belt system [M.S. thesis], Beijing University of Technology, 2005.

S. P. Chen (1) and Y. H. Qian (2)

(1) College of Mathematics, Xiamen University of Technology, Xiamen 361024, China

(2) College of Mathematics, Physics and Information Engineering, Zhejiang Normal University, Jinhua, Zhejiang 321004, China

Correspondence should be addressed to S. P. Chen; shupingchen@126.com

Received 10 May 2014; Revised 24 July 2014; Accepted 11 August 2014; Published 27 August 2014