Advances in Acoustics and Vibration

Volume 2010 (2010), Article ID 284187, 9 pages

http://dx.doi.org/10.1155/2010/284187

## Study of the Response Prediction of a Parametric Vibration Based on Modulation Feedback

Department of Mechanical Engineering and Automation, Shanghai University, 149 Yanchang Road, Shanghai 200072, China

Received 29 September 2010; Accepted 31 December 2010

Academic Editor: Jorge Arenas

Copyright © 2010 Dishan Huang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The problem of response prediction is investigated for parametric vibration in terms of a new concept. The response solution is presented in the special form of Fourier series for signal degree freedom of parametric vibration based on modulation feedback. By applying harmonic balance and limitation operation, all coefficients of a harmonic component are fully determined with a set of series. Meanwhile, some important dynamic behaviors are exposed through mathematical deduction, and an instability phenomenon can be discussed through given frequency factors. The investigation result shows that the new approach has an advantage in the complete response expression, and it is very significant for the theoretical research and engineering application concerning parametric vibration.

#### 1. Introduction

The problem of parametric vibration arises in many branches of physics and engineering, and investigation of stability and response prediction are the two most significant dynamic problems in the parametric vibration system. In the past, several methods have been used to study the stability of systems with periodic coefficients. These include Hill’s Method [1], the perturbation method [2], the averaging approach, Floquet theory with numerical integration [3], and Sinha’s numerical scheme with the shifted Chebyshev polynomial [4, 5].

Hill’s infinite determinant was used to compute the transition curve or the instability boundary. However, it does not yield the response vector at an arbitrary point in the parametric space. The perturbation method supplies a method for obtaining the solution vector and instability boundary. The solution vector is expanded as a power series in terms of the small parametric and substituted into the original equation. The method is applicable only if the parametric multiplying the period term is small [6]. Authors [7–9] have tried to use the approximate approach to determine the stability and response from an approximate system of equations, which is usually obtained by replacing the elements of the periodic coefficients matrix with piecewise constants or linear functions. Sinha and Wu presented a numerical scheme for computing the Floquet transition matrix associated with a class of periodic systems. In the authors’ approach, the solution vector is expanded in terms of the shifted Chebyshev polynomials. This technique reduces the original differential system to a system of linear algebraic equation from which the solutions in the interval of one period can be easily obtained.

In science and engineering practice, it will be very useful to express the response solution in the form of a Fourier series, which is convenient in an online vibration monitoring system. Especially in a mechanical fault diagnosis, such expression is important and not replaceable. One typical example is to diagnose the mechanical fault in a rotational rotor with a cross crack by analyzing a characteristic of the parametric system [10]. However, so far, no related response solution has directly been considered in all of the above approaches.

In this paper, the concept of modulation feedback is introduced to set a response form that is similar to a Fourier series. Using harmonic balance in the equation and limitation operation, a full response expression is conducted. Though this concept comes from the content of a modulation system, it reflects the intrinsic physical character of a parametric system. Moreover, the presented solution conducts the response prediction that could be seen as an expansion of the Perturbation solution.

#### 2. Modulation Feedback Conception

Consider a classical parametric vibration equation: and rewrite it as Based on (2), the problem of parametric vibration can be schematically described as a system drawn in Figure 1, and it is a special feedback system that contains a second-order linear unit and a frequency modulation. The output of the system is the general solution of parametric vibration.

In this system, in addition to enhancing the input and feedback signal, the second-order linear unit will inherently generate the harmonic component after the moment . However, the generated component will be multiplied by , and this operation is called the frequency modulation, which conducts frequency fission in the system.

The phenomenon of frequency fission can be stated as the following: the component is fissile into the combination of harmonic components and by a modulation operation within the first . Then, the modulation result directly feeds back to the second-order linear unit as an input of the new harmonic components. Thus, within the second , three components , and will be the output of the system and be involved in another fission operation, a sequential combination of five harmonic components ,,, , and will become the modulation results and will feed back to the second-order linear unit, and the iterate fission operation continues. The whole physical procedure of the frequency fission is schematically described in Figure 2.

Assuming the frequency fission procedure continues, except for the case of the instable state, the frequency fission will attain a state of balance. Therefore, there are many linear combinations of harmonic components related to and in the output of the system as the frequency fission balance, and they can be mathematically expressed as follows: Because the energy of the harmonic components is concentrated on the narrow band range in the frequency domain, the coefficient while .

Therefore, the response solution problem of the parametric vibration equation is simplified to the determination of coefficient values of in (3).

#### 3. Fission Coordinate Equation

After substituting the solution (3) into (1), we obtain the following infinite set of linear, algebraic, homogeneous equations for the by equating the coefficients of each harmonic on both sides of (1): Here, we introduce notations: First, we consider finite equations, and put all coefficient relations together and form a linear algebra equation called a fission coordinate equation (which is similar to Hill’s form).

#### 4. Determination of Coefficient

To get a full response expression, the determination of coefficient is carried out by the following operation.

(a) Consider the lower half of (6): Let frequency factors be Then One of the solutions to (9): As mentioned in the above section, while . Let ; thus the coefficient can be expanded as where For the same reason, can be determined. They are where

(b) Consider the upper half of (6): where the frequency factors The solution of (17) is where From (11), (13), (19), and (20), we obtain The final solution form of (4) is

#### 5. Special Case

While ( is an integral number) and , the response solution (3) is changed to If is an even number, one gets If is an odd number, one gets (a)There is a direct component in the response while is an even number in (29). This is a kind of particular characteristic exposed by the given approach in the response study of parametric vibration.(b)Based on (30) and (31), the theoretical spectrum of the response of the parametric vibration can be easily predicted.(c)While in (30) or in (31), if the coefficients the response will be the case of the instability state.

#### 6. Complete Form of the Solution

##### 6.1. Complete Form of the Solution

According to the theory of second-order differential equation, the strict form of a general solution for the given parametric vibration equation should be Meanwhile, it needs to meet the sufficient and necessary condition that the Wronskian determinant is equal to the constant.

In fact, has no relation with time , where Because is a constant, (33) is the complete form of the general solution of the given parametric vibration equation, where and are the arbitrary constants in the solution in (33) and depend on the initial condition.

##### 6.2. Determination of Coefficient

The coefficient is determined by using the harmonic balance as the determination of the coefficient value of : Thus, in the special case of , if is an even number, the complete form of the solution will be If is an odd number, the complete form of the solution will be

##### 6.3. Initial Condition

Equation (33) is the parametric vibration solution after the system attains a state of balance through frequency fission. However, when , equation (1) degenerates to a linear equation. Therefore, the coefficients and in (33) can be determined by

##### 6.4. Comparison with Another Method

If in (33), we obtain This is a form of the parametric vibration solution conducted by singular Perturbation called the two variables expansion. Thus, the presented method conducts the response solution that is the expansion of the solution obtained by Perturbation.

In fact, equation (6) is similar to Hill’s form, which was used in the study of the instability boundary of parametric vibration during that time [11]. However, this paper determines the coefficients of the harmonic series by using a special property of in (6).

#### 7. Instability

There are infinite terms of frequency factor in the expression of the obtained solution, which causes instability in the system while . The frequency factors are related to and , and the roots of distribute around the instability areas, which are summarized as follows.

##### 7.1. Principal Instability Area

While and are selected, the frequency factor . As a result, the values of and in (23) will tend to be infinite while . This means the harmonic response corresponding to tends to be infinite in the principal instability area.

##### 7.2. High-Order Instability Boundary Curve

Based on the roots of (), the group of harmonic resonance (higher-order instability) boundary curves are drawn in the parameter plan of and , as shown in Figures 3(a)–3(d) (generated by commercial software Maple 9.5), which have the following features.

*Left Bias*

Due to the term having influence on frequency factor , each group of multiple harmonic resonance curves is under bias toward (, and ). In other words, each harmonic resonance boundary curve is located at the left side of the point that is the high-order instability point while is close to 0.

*Instability Region Overlap*

There are different order harmonic resonance boundary curves in each instability area around (). For example, as shown in Figure 4, the multiple harmonic resonance regions exit near the third-order instability area, and they are based on , ; but for , almost overlap in the same region while . Wherein the curve is the left side of the instability region, and the curve (, in theory) is the right side of the instability region. Therefore, the harmonic resonance boundary curves () are composed of the boundaries of the instability region.

*Cross Direction*

As known from (26), the harmonic resonance strength is different between and . Thus, the harmonic resonance strength is sensitive to the cross direction (from left to right or from right to left) to the high-order instability area.

#### 8. Example

Let , , and in (1), the initial conditions are and .

From (37), the response solution can be expressed as On the initial condition, obtain .

The calculation result of coefficient values and the corresponding frequency component are summarized in Table 1 in which the theoretical spectrum of the response solution of parametric vibration can be predicted. The spectrum is schematically shown in Figure 5.

The response solution is very sensitive to while is closed to 0.5. If is a little less than 0.5, the high-order instability will be reflected in the solution (33). For example, in the case of in the above example, we get . The result can be seen as a case of infinity, and it is governed by the left bias of higher-order instability. Simultaneously, the direct component cannot be neglected. The detailed results are listed in Table 2, and the theoretical response spectrum is drawn in Figure 6.

#### 9. Conclusion

The response prediction of the single freedom of parametric vibration can be determined with a set of combinations of harmonic components related to and . Once the coefficients and are given, whether the general case or the special case of , the harmonic components coefficients and can be fully determined by a set of series. On the coefficients and , the corresponding theoretical spectrum can be drawn; it will be very useful in identifying and predicting the characteristics in the mechanical system on a model of parametric vibration, especially in the case of .

The instability region can be determined when the frequency factor is equal to zero. The investigation result matches the case of Hill’s instability boundary curve. The harmonic resonance in the high-order instability region has three features, which are left bias, region overlap, and cross direction.

The presented concept of modulation feedback reflects the physical inherency in the parametric vibration system, though it comes from the scope of a modulation system. This can become an available tool for us to analyze the response in the parametric vibration. Therefore, it is expected to obtain a response prediction of parametric vibration with an external force excitation and of a multiple freedom of parametric vibration, and so on.

#### Acknowledgment

The authors wish to acknowledge the support of the 211 Leading Academic Discipline Project of Shanghai University.

#### References

- V. A. Yakubovitch and V. M. Starzhinskii,
*Linear Differential Equation with Periodic Coefficients*, John Wiley & Sons, New York, NY, USA, 1975. - A. H. Nayfeh,
*Problems in Perturbation*, John Wiley & Sons, New York, NY, USA, 1985. - G. H. Gaonkar, D. S. Simha Prasad, and D. Sastry, “On computing Floquet transition matrices of rotorcraft,”
*Journal of the American Helicopter Society*, vol. 26, no. 3, pp. 56–61, 1981. View at Google Scholar · View at Scopus - S. C. Sinha and D. H. Wu, “An efficient computational scheme for the analysis of periodic systems,”
*Journal of Sound and Vibration*, vol. 151, no. 1, pp. 91–117, 1991. View at Google Scholar · View at Scopus - S. C. Sinha, D. H. Wu, V. Juneja, and P. Joseph, “Analysis of dynamic systems with periodically varying parameters via Chebyshev polynomials,”
*Transaction of the ASME, Journal of Vibration and Acoustics*, vol. 115, no. 1, pp. 96–102, 1993. View at Google Scholar · View at Scopus - P. P. Friedmann, “Numerical methods for determining the stability and response of periodic systems with applications to helicopter rotor dynamics and aeroelasticity,”
*Computers and Mathematics with Applications*, vol. 12, no. 1, pp. 131–148, 1986. View at Google Scholar · View at Scopus - H. D’Angelo,
*Linear Time-Varying System: Analysis and Synthesis*, Allyn and Bacon, Boston, Mass, USA, 1970. - C. S. Hsu, “On approximating a general linear periodic system,”
*Journal of Mathematical Analysis and Applications*, vol. 45, no. 1, pp. 234–251, 1974. View at Google Scholar · View at Scopus - J. A. Richards,
*Analysis of Periodically Time-Varying System*, Springer, Berlin, Germany, 1983. - R. Gasch, “Dynamic behavior of a simple rotor with cross-sectional crack,”
*IMechE*, vol. 170, pp. 123–128, 1976. View at Google Scholar - A. H. Nayfeh and D. T. Mook,
*Non-Linear Oscillations*, John Wiley & Sons, New York, NY, USA, 1979.