Research Article | Open Access
Dynamical Performances of a Vibration Absorber for Continuous Structure considering Time-Delay Coupling
The nonlinear effect incurred by time delay in vibration control is investigated in this study via a vibration absorber coupled with a continuous beam structure. The stability of the vibration absorber coupled structure system with time-delay coupling is firstly studied, which provides a general guideline for the potential time delay to be introduced to the system. Then it is shown that there is a specific region for the time delay which can bring bifurcation modes to the dynamic response of the coupling system, and the vibration energy at low frequencies can be transferred or absorbed due to the bifurcation mode and the vibration in the corresponding frequency range is thus suppressed. The nonlinear mechanism of this vibration suppression incurred by the coupling time delay is discussed in detail, which provides a novel and alternative approach to the analysis, design, and control of vibration absorbers in engineering practice.
Vibration absorbers including passive and active dynamic vibration absorbers received increasing attention in structural vibration control in the past several years [1–5]. The mechanism of a passive vibration absorber is to distribute the vibration energy through coupling a slave structure to the master one , while an active vibration absorber produces forces to counteract the vibration through adjusting the slave structure according to the motion state of the master .
Considering the mechanism of various vibration absorbers developed in the literature, many different works have been done, which provide a fundamental basis for the understanding of nonlinear dynamics in vibration control. Wong and Cheung  employed the damping characteristics of a vibration absorber to control the master system by dissipating the vibration energy with the damper of the absorber. Kim et al.  considered optimization of the damping characteristics of vibration absorbers. The optimal damping ratio for the master and slave structure under study is therefore obtained so that the amplitude of the master structure can be kept in a small range. Natsiavas  investigated vibration absorbers with nonlinear springs since soft and hard springs have very different effect on system damping characteristics. However, it is noticed that passive vibration absorbers are often unsatisfactory for their uncontrolled reaction under complex working conditions in practice. Huang et al.  used a servo motor as actuator to drive a vibration absorber and achieve remarkable vibration suppression performance. Kobori et al.  proposed a concept of active variable stiffness in active vibration control and explored the potential benefits of nonlinear dynamics incurred by time delay in vibration control.
Noticeably, Xu et al. [13, 14] presented that time delay could be used as a simple but efficient “switch” to control motions of a system for different applications, which was proved by relative experiment of a linear time-delayed control absorber. A delayed-resonator vibration absorber is first proposed by Olgac et al. [15–17]. Vibration energy of a linear master structure is absorbed to suppress undesired oscillations by coupling a time delay to a master mass. To broaden a valid range in frequency for suppression of vibrations, Zhao and Xu  proposed a delayed-resonator vibration absorber with nonlinear master and slave structures, and then a linear time-delayed active absorber was realized in relative experiment [19, 20]. The slave structure in this study has geometrical nonlinearity, which can expand the frequency range so that vibration energy could be absorbed easily in engineering practice. It is shown that the amplitude of the master system is decreased by at least 10% as compared to that without the time-delayed feedback control, and the complex dynamics of the system induced by the time delay are also preliminarily studied. Sun et al. [14, 21] utilized linear and nonlinear velocity feedback to control a continuous beam for parametrical resonance. It is shown that the active control designed suppresses the parametric resonance for force excitation in the theoretical analysis. The significant effect of the vibration control has been proved and the bifurcation in parametric excitation system can be eliminated in the related experiment. In previous studies, the complex nonlinear dynamics induced by time delay are usually considered to be harmful for structural vibration control. However, it is more and more noticed that some predicable nonlinear dynamics could have beneficial influence in vibration control, due to the fact that multiple vibration modes could result in energy transmission or suppression under certain conditions. Although some preliminary results have shown that the nonlinear dynamics incurred by time delay in the coupling of slave and master structures could be advantageous in absorbing vibration energy, a systematic and theoretical study on the nonlinear mechanism in affecting the stability of continuous structures and absorbing vibration energy is still yet to be further developed.
In this paper, to investigate the nonlinear mechanism of vibration absorbers with time delays, a vibration absorber system (VAS) under harmonic excitation is studied. The VAS includes a continuous beam structure as the master structure and a continuous structure consisting of two beams as the slave system, which also acts as an absorber. After coupling the master structure and the slave system with a spring and a servo motor, there exists energy transmission between these two structures . Basically, the geometry of the slave structure and the coupling strength should be carefully designed so that the frequency of vibration modes bifurcated by the introduced time delay could be proportional to the excitation frequency. In this way, it is known that the coupled modes of nonlinear systems can realize energy transmission and possible suppression . The slave structure in this study has geometrical nonlinearity, which can expand the frequency range where vibration energy could be absorbed . In the present research, it is shown that there is a specific region for the time delay which can bring bifurcation modes to the dynamic response of the system, while the latter can absorb the vibration energy at low frequencies and thus suppress the vibration over the corresponding frequency range. The nonlinear mechanism of this vibration suppression incurred by the coupling time delay is discussed in detail. This provides a novel and alternative insight into the analysis, design, and control of vibration absorbers in engineering practice.
The paper is organized as follows. Firstly, the continuous structure system is modeled and the nonlinear dynamics of both master and slave structures are considered in Section 2. The relationship between the coupling spring stiffness and the amplitude of the master structure is derived. Secondly, the nonlinear effect induced by the time delay via the servo motor is studied. The stability region is determined through the eigenvalues of vibration equations in Section 3. Owing to the stability of fixed point, the master structure is under periodic motion which can be described by Harmonic Balance Method (HBM) [25, 26]. Complex dynamical behavior would appear because of the bifurcation induced by time delay. The beneficial region of the time delay for vibration suppression is developed in Section 4. The quasi-periodic solution is solved by Perturbation Incrementation Scheme (PIS) [27, 28] and Multiscale Method (MSM) . The nonlinear mechanism of the vibration absorber under different regions of the time delay is also discussed in Section 4. A conclusion is given in Section 5.
2. The Nonlinear Vibration and Its Modeling
The system under consideration is shown schematically in Figure 1. The master (Beam 1) and the slave structures (Beams 2 and 3) are coupled by two springs and a servo motor since the slave structure may perform a geometrically nonlinear nature , where and are the linear stiffness of coupling and controllable springs, respectively. The servo motor is used to provide the strength of and the coupling delay in demand so that and the delay are considered as two controllable parameters to suppress or absorb vibrations of the master structure. It should be noted that the equivalent mass of the slave structure is much smaller than the master structure for the convenience of the installation and realization. The geometric sizes of the three beams are indicated in Figure 1 and the lengths of Beams 1, 2, and 3 are m, m, and m, respectively. The material of the model is steel, and Young’s modulus of the material is GPa.
The state based actuator in the control spring by the servo motor is given by with time delay . The potential damping effect of the time delay would be focused for this control strategy. The dynamical equations of the continuous structure with the actuator are
All the parameters in the equations above such as , , , , , and are defined in Physical Parameters of the Structure. The deflections of each beam are , , and . The high order modes are ignored due to our goal to suppress or absorb the buckling vibration in the present research. Therefore, the dynamics of the structure are truncated to the first mode as , , and . The first mode of each beam in the structure is determined by boundary conditions. The master beam (Beam 1) is a simply supported beam at both ends and Beam 2 is fixed at both ends. Beam 3 is fixed at one end and is fixed to the midpoint of Beam 2 at another end, which provides a boundary condition by deformation compatibility with the midpoint of Beam 2. The first mode of each beam is given as follows:
Utilizing the compatible displacement , the coefficients of mode equations above can be computed as , , and , and the system is a two-degree-of-freedom vibration system because of the compatibility condition of Beam 2 and Beam 3. The vibration equations can then be obtained aswhere
3. Natural Frequency and Stability of the System
3.1. Without Time Delay
The dynamic response of the system without time delay is considered firstly. The natural frequency of the structure can be obtained through solving the eigenvalue problem of the linear part of (3). The natural frequencies of the structure are defined as and . Ignoring the damping effect, the frequencies can be solved bywhere the stiffness matrixand the system mass matrixThe natural frequencies can be obtained as
The relationship between the natural frequencies and of the structure and the stiffness of coupling spring is shown in Figure 2. The coupling spring actually increases the natural frequency of the structure. This results in the first natural frequency of the master structure being moved to a higher frequency. The vibration amplitude of the master structure under different stiffness of the coupling spring is given in Figure 3. The vibration amplitude of the master structure decreases after the slave structure is coupled through the coupling spring.
Figure 3 shows that the coupling spring helps to suppress the vibration of the master structure. However, as the stiffness of the coupling spring increases from zero to 4000 N·m−1, the amplitude of the master structure decreases from 1.6 mm to 1.1 mm. The vibration suppression due to the coupling spring actually is not significant to meet engineering requirements. For vibration suppression effect, a time-delayed active control is proposed in the following section via the servo motor, which makes the slave structure a much more powerful vibration absorber. The time delay can be adjusted and controlled to give the best vibration suppression.
3.2. With Time Delay
Because of the unsatisfactory vibration suppression effect of the coupling spring, the time delay is proposed to be introduced in this section as a key element for performance improvement. Owing to the fact that the amplitude of the external excitation could be small, the small external excitation in this paper cannot cause complex dynamics in that the magnitude of the external excitation is much smaller than other terms in the equation, but the nonlinearity and the time delay could induce bifurcation and even unstable dynamics, which should be systematically investigated.
Introducing the following dimensionless variables,then the vibration equations can be rewritten as
For the values of dimensionless parameters, please refer to Table 2.
After setting , the linear part of the dimensionless vibration equation can be written as
The characteristic equation can be obtained from (12) as
The imaginary part of eigenvalues is the natural frequency. Two eigenvalues of the characteristic equation (14) can be computed, whose imaginary parts are the smallest and real parts are the biggest, and then the first and second natural frequencies of the continuous structure can be obtained, which are shown in Figure 4.
From Figure 4, it can be seen that the second natural frequency decreases quickly as the time delay increases from 1.0 to 6.532. When time delay reaches 6.532, the first natural frequency can drop. The decrease of the first natural frequency makes it easier for the system to be resonant under low frequency periodic excitation. Therefore, the value of time delay should not be larger than 6.532.
The real parts of eigenvalues determine the stability of the corresponding mode. The complex dynamics would be caused if there is more than one eigenvalue having positive real parts. The real parts of the eigenvalues corresponding to the first and second frequencies above can be obtained by solving (14) and thus shown in Figure 5.
In Figure 5, the dynamic behavior of the structure is periodic subject to external excitation when where no eigenvalues have positive real parts. The structure has complex dynamic behavior as due to the appearance of eigenvalues with positive real parts.
From the discussions above, the nonlinear effects of the time delay on natural frequency and stability indicate that the time delay can greatly change the dynamic response of the structure through the variation of natural frequency and bifurcation behaviors induced. The detailed nonlinear mechanism and influence of the time delay on vibration suppression will be discussed further in the following section.
4. The Vibration Response of the Nonlinear Vibration Absorber with Time Delay
4.1. Periodic Vibration Response with Time Delay
In this section, the output spectrum of the vibration response of the master structure is derived with the Harmonic Balance Method for the time delay . To this aim, it is assumed that the amplitudes caused by external excitations are as follows:
Substituting the expressions into (10) and finding the coefficients proportional to and , the following nonlinear algebraic equations can be obtained:
The amplitude of the master structure is . Considering only the case when , Figure 6 presents the frequency response curves under different values of the time delay.
Figure 6 presents the frequency response amplitude of the master structure around the resonance frequency rad/s under , , , and . The time delay plays a significant role in vibration suppression when . It is noticed that, with the value of the time delay increasing, the peak response amplitude of the master structure decreases obviously. Also, the curve corresponding to jumps drastically when rad/s, which shows that the cubic coefficients and affect the dynamic response of the structure. Figure 7 is the time-series response of the master structure for and . It is shown that the vibration amplitude decreases by 50% with compared to that without time delay as rad/s. The mechanism of the vibration absorber with time delay is that the introduced time delay adjusts the output response of the servo motor so that the active force acting on the master structure could counteract the original motion.
4.2. Coupled Vibration Modes
It has been mentioned that the system characteristic equation will have eigenvalues with positive real parts for . The bifurcation behavior would be incurred by the time delay, which can be coupled with the nonlinear motion of the system induced by external excitation. When the natural frequency of the master structure is designed around half of resonance frequency of the slave structure, vibration energy transmission could happen between these coupled modes.
4.2.1. Derivation of the Solutions of Vibration Modes
To investigate the coupling vibration modes, the vibration solution can be derived, which can be expressed in two parts including motions induced by the time delay and external excitation, respectively. The frequency of the first part in the solution should equal the input frequency of external excitation while the frequency of the second part should equal the frequency at bifurcation point. From (14), the critical time delay for producing the bifurcation motion is and the dimensionless frequency of bifurcation mode is . The solution for can be expressed as follows:where is the motion by the external excitation and is the bifurcation motion caused by the time delay. As mentioned in (17), the frequency of equals the input frequency of external excitation, and the frequency of equals the frequency at bifurcation point.
The two parts of the solution are difficult to be derived because the two modes are coupled by nonlinear terms. But the linear part of the vibration equations can be separated to find the eigenvector of each mode. Substituting the expression of solutions (17) into (10), the linear part can be separated from the equations. The expression of the motion excited by the external excitation can therefore be simplified.
Two time scales are introduced and given by and ; then one has
Perturbing in asand expanding the vector , the time-delay items are as follows:
The first subscript is the number of degrees of freedom, the second is the number of modes, and the third is the order of perturbation. Substituting (18)–(20) into the dynamical equation (10), one has
The frequency of the external excitation is near the natural frequency of the structure while the bifurcation frequency is close to . Tuning parameters and are introduced to describe the difference between and and and , respectively. The following expressions for and are given:
For , the nontrivial solution to (10) can be written aswhere stands for the complex conjugation of the preceding terms and the eigenvector is determined by and is determined by
To obtain the amplitude and frequency of the motion, and are set as
After substituting the values of parameters, the critical value of time delay and , into the equations above, the normal form of (28) is given in (A.1). The solution of the two vibration modes can therefore be obtained by solving , , , and .
4.2.2. The Bifurcation Modes and Critical Points
Figure 8 shows the vibration amplitude of the master structure as the time delay varies. There are three bifurcation points in Figure 8(b), which are , , and . The resonance behavior among different modes and the transmission of energy are different at , , and due to the fact that the critical time delay and bifurcation frequency at these three bifurcation points are different. The amplitude drops continuously when and jumps from 1 mm to nearly 0.5 mm at where . The amplitude curves are discontinuous at , which means that the bifurcation undergoes a subcritical Hopf bifurcation induced by the time delay at point and makes the amplitude of the master structure jump.
As time delay , the vibration amplitude keeps at a low level. The time delay is helpful for vibration suppression in this range. The second critical value of the time delay for bifurcation is where . The frequency of bifurcation mode is and is close to two times the first natural frequency. The amplitude of the master structure increases continuously for the occurrence of supercritical bifurcation at point where the bifurcation mode and excitation mode are in weak resonance. The weak resonance helps the vibration suppression in that the bifurcation mode absorbs the vibration energy of the master structure. The vibration energy is transferred from to because of mode coupling by the nonlinearity of the system. The slave structure absorbs the vibration energy from the master structure.
When the time delay reaches the value , the third Hopf bifurcation occurs. From Figure 5, an eigenvalue crosses the real axis, which results in a subcritical Hopf bifurcation. The amplitude of the master structure jumps to a large value there in Figure 8. At this bifurcation point, the frequency of the bifurcation mode is , which is close to three times the first natural frequency as . The two modes are in strong resonance when time delay . The strong resonance cannot help the vibration absorbing because of the extra energy imported.
Figure 9 can clearly show the ratio of the frequencies of the motion at around the bifurcation points and . The frequency spectrum of the dynamic response of the master structure contains two frequencies subject to a single-tone excitation. The ratio of the frequencies of bifurcation mode and the first mode is 1/2 at around point while it is 1/3 at around . Due to the magnitude level, the former mode can be regarded as being in weak resonances while the latter can be regarded as being in strong resonance.
4.2.3. Vibration Performance in Different Bifurcation Modes
Figure 10 shows the dynamic behavior of the master structure as . Figure 10(a) is the time-series response of the master structure; Figure 10(b) is the phase diagram and Figure 10(c) is the Poincare mapping section. There is a cycle on the Poincare mapping, which indicates that the dynamic behavior of the system is quasi-periodic motion including two irreducible frequencies. As the bifurcation mode and excitation mode are in 1/2 weak resonance, the amplitude of the master structure can be controlled in a small range. The vibration suppression performance of the nonlinear vibration absorber is significant for owing to the fact that the energy could be transferred from the excitation mode to the bifurcation mode.
Figure 11 shows the motion of the master structure as . The system undergoes a subcritical bifurcation at where . From Figure 11, the dynamics of the structure contain two modes and when the vibration amplitude of the master structure is much larger than that when . The ratio of the bifurcation mode frequency and the first natural frequency is about 1/3 here; thus the two modes incurred by bifurcation and excitation are in strong resonance. The amplitude of the master structure jumps to 10 mm at for the strong energy transmission between the two modes. The amplitude of the master structure cannot be suppressed because the time delay cannot dissipate vibration energy in this case. Therefore, different energy transmission mechanism in weak and strong resonance results in very different effect in vibration suppression.
4.2.4. Energy Transfer Phenomena between Two Vibration Modes
It has been shown that the motion bifurcated by the time delay is helpful to vibration suppression for . The vibration absorber effect for is more significant than that in other ranges because the vibration absorber may utilize the energy transmission of coupling modes to realize the active vibration absorbing. The mode bifurcated by the time delay could absorb the vibration energy of external excitation.
To demonstrate this more clearly, Figure 12 indicates the energy distribution of the vibration response under different time delays. When , the vibration of the master structure contains huge energy at the low frequency. The frequency spectrum for shows that the vibration energy diminishes with the time delay increasing while it keeps the same frequency as . When the time delay increases to , the vibration energy at the first frequency reduces several times lower than that when . Also, there is vibration energy at the second frequency, which implies that the bifurcation induced by the time delay transfers some energy from a low frequency to a high one. When , the dynamic response of the system contains more than two frequencies and the vibration energy is higher than other cases. It can be seen that the phenomena observed here are consistent with these by the stability analysis in Section 3.
The most significant region of the time delay affecting vibration absorbing is . The value of the time delay in the analysis above is a dimensionless parameter. After carrying out an inverse procedure of dimensional transformation, the value of time delay in practical engineering can be obtained by (s). As a result, the best range of time delay in practical engineering is that (s), which can be realized in practice.
In this study, the time-delayed active control of vibration absorbers coupled with a continuous beam structure is studied. The vibration amplitude of the master structure decreases nearly by 50% when the time delay falls into a specific range; that is, . The time-delayed vibration absorber engenders opposite motion on the master structure, which can counteract part of its vibration energy. Moreover, it is shown that different values of time delay have different effects on vibration absorbing behavior of the structure. The bifurcation mode induced by the time delay helps the vibration absorbing for . The mode bifurcated by the time delay and the first mode of the system are in 1/2 weak resonance, which can transfer some vibration energy out of the mode at low frequency and therefore decrease the amplitude of the master structure at corresponding frequency range. However, at , the 1/3 strong resonance makes the amplitude of the master structure enlarged because the bifurcation caused by the time delay amplifies vibration energy of the system. The results obtained in this study could provide a significant insight into the analysis and design of active vibration absorbers with deliberately introduced time delay by exploiting nonlinear benefits for realizing superior vibration suppression.
The normal form at around for (28) is given by
Physical Parameters of the Structure
|:||Length of each beam (m)|
|:||Damping coefficient (Ns)|
|:||Stiffness of coupling spring (N·)|
|:||Stiffness of control spring (N·)|
|:||Elastic modulus of the material, (Pa)|
|:||Moment of inertia (m4)|
|:||Frequency of external force (rad·)|
|:||Amplitude of external force|
|:||Deflection of each beam|
|:||Cross-section area of each beam (m2)|
|:||Linear density of each beam, (kg·)|
|:||Time delay in control (s).|
The authors declare that they have no competing interests.
The authors would like to gratefully acknowledge the support from Shanghai Sailing Program (16YF1408000), Natural Science Foundation of Shanghai (16ZR1423600), and the National Nature Science Foundation of China (11572224).
- Y. Du, R. A. Burdisso, and E. Nikolaidis, “Control of internal resonances in vibration isolators using passive and hybrid dynamic vibration absorbers,” Journal of Sound and Vibration, vol. 286, no. 4-5, pp. 697–727, 2005.
- H. L. Sun, P. Q. Zhang, H. B. Chen, K. Zhang, and X. L. Gong, “Application of dynamic vibration absorbers in structural vibration control under multi-frequency harmonic excitations,” Applied Acoustics, vol. 69, no. 12, pp. 1361–1367, 2008.
- T. Mizuno, M. Moriya, and K. Araki, “Robust disturbance cancellation in an active dynamic vibration absorber system,” Control Engineering Practice, vol. 3, no. 6, pp. 773–781, 1995.
- S. Chatterjee, “Optimal active absorber with internal state feedback for controlling resonant and transient vibration,” Journal of Sound and Vibration, vol. 329, no. 26, pp. 5397–5414, 2010.
- S.-M. Kim, S. Y. Wang, and M. J. Brennan, “Optimal and robust modal control of a flexible structure using an active dynamic vibration absorber,” Smart Materials and Structures, vol. 20, no. 4, Article ID 045003, 2011.
- N. D. Anh, H. Matsuhisa, L. D. Viet, and M. Yasuda, “Vibration control of an inverted pendulum type structure by passive mass-spring-pendulum dynamic vibration absorber,” Journal of Sound and Vibration, vol. 307, no. 1-2, pp. 187–201, 2007.
- M. Sayed and M. Kamel, “1:2 and 1:3 internal resonance active absorber for non-linear vibrating system,” Applied Mathematical Modelling, vol. 36, no. 1, pp. 310–332, 2012.
- W. O. Wong and Y. L. Cheung, “Optimal design of a damped dynamic vibration absorber for vibration control of structure excited by ground motion,” Engineering Structures, vol. 30, no. 1, pp. 282–286, 2008.
- D. H. Kim, J. W. Park, G. S. Lee, and K. I. Lee, “Active impact control system design with a hydraulic damper,” Journal of Sound and Vibration, vol. 250, no. 3, pp. 485–501, 2002.
- S. Natsiavas, “Steady state oscillations and stability of non-linear dynamic vibration absorbers,” Journal of Sound and Vibration, vol. 156, no. 2, pp. 227–245, 1992.
- S.-J. Huang, K.-S. Huang, and K.-C. Chiou, “Development and application of a novel radial basis function sliding mode controller,” Mechatronics, vol. 13, no. 4, pp. 313–329, 2003.
- T. Kobori, H. Kanayama, and S. Kamagata, “Active seismic response control systems for nuclear power plant equipment facilities,” Nuclear Engineering and Design, vol. 111, no. 3, pp. 351–356, 1989.
- J. Xu and K. W. Chung, “Effects of time delayed position feedback on a van der Pol-Duffing oscillator,” Physica D, vol. 180, no. 1-2, pp. 17–39, 2003.
- Y. Sun and J. Xu, “Experiments and analysis for a controlled mechanical absorber considering delay effect,” Journal of Sound and Vibration, vol. 339, pp. 25–37, 2015.
- N. Olgac and B. T. Holm-Hansen, “A novel active vibration absorption technique: delayed resonator,” Journal of Sound and Vibration, vol. 176, no. 1, pp. 93–104, 1994.
- N. Olgae and B. Holm-Hansen, “Design considerations for delayed-resonator vibration absorbers,” Journal of Engineering Mechanics, vol. 121, no. 1, pp. 80–89, 1995.
- N. Olgac and M. Hosek, “Active vibration absorption using delayed resonator with relative position measurement,” Journal of Vibration and Acoustics, vol. 119, no. 1, pp. 131–136, 1997.
- Y.-Y. Zhao and J. Xu, “Effects of delayed feedback control on nonlinear vibration absorber system,” Journal of Sound and Vibration, vol. 308, no. 1-2, pp. 212–230, 2007.
- L. Chen, C. H. Hansen, F. He, and K. Sammut, “Active nonlinear vibration absorber design for flexible structures,” International Journal of Acoustics and Vibrations, vol. 12, no. 2, pp. 51–59, 2007.
- L. Chen, F. He, and K. Sammut, “Vibration suppression of a principal parametric resonance,” Journal of Vibration and Control, vol. 15, no. 3, pp. 439–463, 2009.
- K. A. Alhazza and M. A. Majeed, “Free vibrations control of a cantilever beam using combined time delay feedback,” Journal of Vibration and Control, vol. 18, no. 5, pp. 609–621, 2012.
- C. T. Chatzigogos, A. Pecker, and J. Salençon, “Macroelement modeling of shallow foundations,” Soil Dynamics and Earthquake Engineering, vol. 29, no. 5, pp. 765–781, 2009.
- A. Abe, Y. Kobayashi, and G. Yamada, “Nonlinear dynamic behaviors of clamped laminated shallow shells with one-to-one internal resonance,” Journal of Sound and Vibration, vol. 304, no. 3-5, pp. 957–968, 2007.
- J. V. Ferreira, Dynamics response analysis of structures with nonlinear components [Ph.D. thesis], Department of Mechanical Engineering, Imperial College, London, UK, 1998.
- A. Raghothama and S. Narayanan, “Bifurcation and chaos in geared rotor bearing system by incremental harmonic balance method,” Journal of Sound and Vibration, vol. 226, no. 3, pp. 469–492, 1999.
- L. Xu, M. W. Lu, and Q. J. Cao, “Nonlinear vibrations of dynamical systems with a general form of piecewise-linear viscous damping by incremental harmonic balance method,” Physics Letters. A, vol. 301, no. 1-2, pp. 65–73, 2002.