Comparison of translational and rotational motion parameter characteristics.

## Abstract

This chapter will give an introduction to linear and nonlinear oscillators and will propose literature to this topic. Most importantly, hands on examples with numerical simulations are illustrating oscillations and resonance phenomena and where useful, also analytical methods to treat nonlinear behavior are given.

### Keywords

- parametric resonance
- autoparametric resonance
- nonlinear vibration
- Mathieu equation
- Hopf bifurcation
- Strutt diagram
- nonlinear natural frequency
- instability domain
- basepoint excited primary and secondary system

## 1. Introduction

When a mechanical system has at least two vibrating components, the vibration of one of the components may influence the other component. This influence effect which might stabilize or destabilize the system is called autoparametric resonance. This chapter will introduce autoparametric resonance by examining hands on examples for such systems. In particular, basepoint excited systems are analyzed. Beside purely mechanical systems, also examples of an electrical system with two coupled resonators are investigated.

There are three main types of oscillation: (1) free oscillation, (2) forced excited oscillation and (3) self-excited oscillation.

Free oscillation is defined as temporal fluctuations of the state variables of a system. Such temporal fluctuations can be defined as deviations from a mean value. Vibrations are present in many mechanical systems and occur always in feedback systems. The concept of free oscillation is misleading since nearly all physical systems are subject to attenuation. However, it depends on the size (and thus the time). Exceptions could be, for example, orbit oscillations of planets (macroscopic) or oscillations of electrons (microscopic). The two systems mentioned are also subjected to a type of damping, since both systems cannot remain stable indefinitely, but for an extremely long time.

A forced excited spring mass system might be a mechanically forced oscillator. Such systems of translational motions are discussed in Sections 2 and 3. Beside translatory oscillations, rotatory oscillations and resonance is of vital interest to design engineers of aircraft turbines, etc. Unbalanced rotating machine parts are sources of unwanted vibrations and might resonate when excited accordingly.

Self-excited oscillation, also called as self-oscillation, self-induced, maintained or autonomous oscillation is known in electronics as parasitic oscillation and in mechanical engineering literature as hunting. Such systems are discussed in Section 3.

Table 1 depicts relevant parameters for characterization motion in translational and rotational structures. The parameters for displacement, velocity and acceleration have been written as absolute values – knowing that depending on the application, they might be vectors, depending on the chosen frame of reference. In the most general case, they form a four vector. The force is written as mass times acceleration (Newton’s second law) and therefore force is also a vector. That brings us to Newton’s first law, which states that an object that is at rest will stay at rest unless a force acts upon it or inversely an object will not change its velocity unless a force acts upon it. For completeness, also Newton’s third law shall be given: Actio et Reactio – all forces between two objects exist in equal magnitude and opposite direction. A treaty to Newton’s laws of dynamics can be found, for example in chapter 9 of volume I [1].

Translational | Rotational | ||||
---|---|---|---|---|---|

Symbol | Description | SI Unit | Symbol | Description | SI Unit |

Displacement | m | Angle | rad | ||

Velocity | Angular velocity | ||||

Angular acceleration | |||||

Mass | kg | Inertia | kg m^{2} | ||

= | Force | N | = | Torque | Nm |

= | Momentum | Ns | = | Angular momentum | Nms |

Kinetic energy | Nm | Kinetic Energy | Nm | ||

Potential energy | Nm | Potential energy | Nm | ||

= ∫ | Work | J | = ∫ | Work | J |

= | Power | W | = | Angular power | W |

D'Alembert’s principle is a statement of the fundamental classical laws of motion. It is the dynamic analogue to the principle of virtual work for applied forces in a static system and in fact is more general than Hamilton’s principle, avoiding restriction to holonomic systems^{1}. If the negative terms in accelerations are recognized as inertial forces, the statement of d'Alembert's principle becomes “the total virtual work of the impressed forces plus the inertial forces vanishes for reversible displacements”. The principle does not apply for irreversible displacements, such as sliding friction, and more general specification of the irreversibility is required. A derivation of the Lagrangian equation of motion is well explained in Chapter 1 of [2] or [3]. In (1), the non-conservative energy term defined as the Lagrangian (* L*) is, composed of the kinetic energy

*and the potential energy*T

*.*U

In (2) the Lagrangian equation is given with generalized coordinates * q* of a dynamic system and dissipative generalized forces

_{i}

*.*Q

_{i}

The sum of all kinetic energies * T* in the system, whether translational or rotational character (see also Table 1) needs to be included. The sum of all potential energies

*in the system, whether it stems from the gravitation or energy from linear or nonlinear springs or whatever scalar field. Elastic potential energy from any linear or nonlinear spring can be obtained calculating its potential energy. The Lagrangian formalism can also be used for mechanical systems with mass explicitly dependent on position, see for example [4]. In the following chapters, all analyzed dynamical systems are derived using this elegant and powerful method.*U

Looking at translational (classic) mechanical springs, the displacement dependent force * F*(

*) can be written as shown in (3) using spring stiffness*q

_{t}

*and having also introduced a nonlinear spring term*q

_{t}

*and an exponent*k

_{n}

*for setting nonlinearity of spring. For linear springs, where the force is proportional to the displacement (Hooke’s law), this nonlinear spring term will be zero.*n

The elastic energy * E* of the spring is obtained by integrating the exerted force over its covered path s.

_{elastic}

The elastic energy * E* (a potential U) for a linear translational spring is given also in Table 1. For (rotational) torsion springs, the procedure is the same and a given spring torque can be expressed as shown in (5) with torsion coefficient

_{elastic}

*and having also introduced a nonlinear torsion spring term*q

_{r}

*.*D

_{n}

The elastic rotational energy can be expressed as follows (6):

In Figure 1, sketches on the top show translational springs (from left to right linear, nonlinear and nonlinear unsymmetrical) and sketches on the bottom depict rotational springs (from left to right linear, nonlinear and nonlinear unsymmetrical). The origin is depicted with an O and the spring displacement is depicted as * y* and

*, respectively. For the translational magnetic spring systems, the lower and upper magnets are fixed to the reference frame with origin O. The nonlinear symmetric magnetic spring uses three identical block or disk magnets. Such magnetic springs have a nonlinear term*ϕ

*> 0 – forming so called hardening springs if*k

_{n}

*< 0 in literature referred to softening springs, see for example [5]). The nonlinear unsymmetrical translational magnetic spring is also shown in the neutral position and the displacement around the origin is unsymmetrical.*k

_{n}

The rotational magnetic spring systems have a ferromagnetic stator fixed to the reference frame and a rotating hollow shaft carrying two permanent magnets, in this scenario also made of ferromagnetic material. The spring is drawn in the unstable equilibrium position. Exemplarily spring characteristics for such translational and rotational spring systems are depicted in Figure 2.

Note that integrating the force or torque response will add up to zero for rotational and translational spring systems. The given exemplarily spring characteristics of drawn spring systems in Figure 1 are shown in Figure 2. On the left-hand side, exemplarily translational spring characteristics are shown and on the right-hand side rotational spring characteristics are depicted. The linear case where the force and torque are proportional to the displacement is shown in blue. In red, symmetric nonlinear – here for translational and rotational systems a permanent magnet system is shown, but also mechanical spring systems could be envisaged. The curves of the nonlinear asymmetric cases are depicted in orange. In the appendix A.1, more simulations have been depicted for translational symmetric spring systems using ring magnets.

Equivalence of electrical and mechanical systems are shown in Table 2. On the left-hand side, a mechanical system with only one degree of freedom in * y* direction is shown and its equivalent electrical structure with charge

*and current*q

*on the right-hand side. Kinetic energies are denoted with*i

*(translational kinetic energy) and*T

_{T}

*(inductive kinetic energy), potentials are written as*T

_{L}

*(spring potential) and*U

_{S}

*(capacitive potential) and the gravitational potential denoted as*U

_{C}

*and the DC battery voltage is*U

_{G}

*. Non-conservative components are in the mechanical system the viscous damping force and in the electrical system the electrical resistor.*U

_{B}

The Lagrange energy function is shown for the mechanical system in (7) and its equivalent electrical system in (8).

Applying the Lagrangian formalism (1) and (2) to these SDoF systems, will lead to the resulting DE’s as shown in (9) and (10). Note that the sign of the (viscous) damping must be introduced always with a negative sign using the generalized coordinates – as its velocity is always opposing the system velocity. In the electrical circuit, having the flowing charge velocity ^{’} for example, current * i* defined in clockwise direction, the battery voltage, as it is a source, must act in the opposite direction and therefore this potential energy must be introduced with a negative sign.

Both systems (9) a force DE, (10) a voltage DE, belong to the same class of ordinary linear second-order DE. Resonance frequency for the mechanical system (9) is

## 2. Linear resonance systems

### 2.1. Linear single degree of freedom systems

In this section, a linear basepoint excited single degree of freedom systems is discussed. The lumped parameter model for the examined system (Figure 3) consists of a linear oscillator with mass * m* damping factor

*, a linear spring with a spring rate*d

k

_{1}and an external basepoint excited harmonic force with amplitude

y

_{0}(

*). System coordinate origin is placed at the basepoint excitation.*t

The kinetic energy (11) of this system and the potential energy (12) form the non-dissipative energy of this system. The dissipative force * F* of this system – we consider only viscous friction – is shown in (13) and the driving force

_{d}

F

_{0}(14) assuming a harmonic basepoint excitation with amplitude

*and driving frequency*A

*.*ω

Applying the Lagrangian formalism (1) and (2), we deal with SDoF system, will lead to the resulting DE shown in (15). Note that the sign of the viscous damping must be introduced with a negative sign using the generalized coordinates. The driving force _{0}, as it is a harmonic signal, can be introduced with a positive or a negative sign, resulting in a phase shift of 180°.

Introducing dimensionless notation, by using a dimensionless time * τ*, a dimensionless system resonance frequency Ω, the damping factor

ξ

_{1}and the gravity offset term

*(17) and setting the dimensionless displacement*ϱ

*(18).*u

By replacing parameters of (16) with (17, 18), we obtain (19). We can drop the gravity offset term * ϱ*, as it will only add a non-time dependent offset to the solution

*(20).*u

Frequency domain behavior is obtained by applying the Laplace Transformation (21–24) and by replacing * s* =

*Ω we obtain the frequency response (25).*j

* G* represents the relative motion of the oscillation. As long as the excitation frequency can be represented by a Fourier series of harmonic functions, this obtained solution is valid and a very powerful result. (24) and (25) are represented in Figure 4. The advantage of the Bode Plot representation is to have also the phase shown. As smaller the dimensionless damping

ξ

_{1}become, as larger becomes the scaled resonance at the dimensionless frequency ration Ω.

### 2.2. Linear two degree of freedom systems (2DoF systems)

In this section, a linear basepoint excited two degree of freedom systems is discussed. The lumped parameter model for the examined system (Figure 5) consists of two linear oscillators with mass _{1} and _{2}, damping factors _{1} and _{2}, linear springs with spring rates _{1} and _{2} and an external basepoint excited harmonic force with amplitude _{0}(* t*). System coordinate origin is placed at the basepoint excitation.

Applying the Lagrangian formalism (1) and (2) to this 2DoF problem, we obtain the coupled DE system shown in (30) and (31).

DE system shown in (30), (31) is represented dimensionless in equation DE system (32), (33) using the dimensionless parameters of (34), (35).

Similar to dimensionless parameters of (17), (18), the dimensionless time * τ*, a dimensionless system resonance frequency Ω, damping factors

ξ

_{1}and

ξ

_{2}, a gravity offset term ϱ, system resonance frequencies of each oscillator

ω

_{1}and

ω

_{2}plus a mass ratio λ_m and an oscillator frequency ratio Ω

_{0}plus dimensionless displacements

*and*u

*.*v

The frequency response of this coupled oscillator system can again be obtained using the Laplace transformation introduced in (21). The system in the frequency domain is shown in (36) and (37) using the same steps as shown in the SDoF system (22)–(25).

As (36) and (37) represent two algebraic equations, * U*(

*) and*s

*(*V

*) can be separated, resulting in (38) and (39). Note that the gravity term*s

*in the numerator will introduce an additional damping of the transfer function.*ϱ

Figure 6 depicts the relative oscillation response in the frequency (left) and time domain (right) of the derived 2DoF system. The frequency response is given as a dimensionless ratio Ω, see also (34). The dimensionless simulation parameters have been set exemplarily to _{1} = _{2} = 0.021, * λ* = 0.42 and Ω

_{m}

_{0}= 1. As we have two resonators with same system frequencies

_{0}= 1 and 0.82 Ω and at 6.1 Ω for Ω

_{0}= 5 (dashed lines).

The time-domain response from this coupled DE system with lumped parameter model Figure 5 and (32) and (33) is shown in Figure 7. On the left-hand side, the ^dimensionless basepoint acceleration signal is given and its dimensionless response signals of first (blue) and second (red) DoF, simulating 50 periods and starting with settled initial conditions (amplitude of

## 3. Nonlinear resonance systems

In the introduction, Section 1, we distinguished three cases of vibration. The class forced excitation will be further investigated in this section. In Figure 8 five systems are depicted that can potentially exhibit parametric resonance effects. The term parametric means that of cases where the external excitation appears as a time varying modification of a system parameter. A “normal” forced excitation system whether linear or nonlinear, will respond to the excitation with or without resonance using the energy fed into it and no time varying modification of a system parameter might excite additionally the system.

The five depicted systems in Figure 8 might show an exponential amplitude growth when excited externally in presence of a system damping factor. In the two electrical systems on the right-hand side, any of the three components R (here not drawn), L or C that is parametrically excited will respond with an exponential amplitude growth, if the mathematical physical system model has at least one degree of freedom of the Mathieu DE (40) or the Hill DE (41).

The Hill differential equation is a generalized form of (40), in which the harmonic function is replaced with any periodic function, shown in (41).

It is most interestingly that any system parameter including also damping factors with time varying influence of a system parameter will result in an exponential growth of the response amplitude. To give a concrete example of this behavior, we consider here the example from Section 3.2 and inspect the resulting (dimensioned) DE system with (62) as primary system and (63) as secondary system of such a behavior.

The primary system has no such configuration, but the secondary system (63) is of Mathieu type. To simplify the treated system, we use instead of the basepoint excitation _{0} the primary system * y*, compare also the lumped parameter model in Figure 12 (in an experiment we would simply make the stiffness

*of the system very large, for example, replacing the spring with a fixed stiff rod). Now the new induced basepoint excitation*k

*will excite the secondary system directly. As*y

*is appearing in (63) as acceleration, we adjust this basepoint excitation simply in form of an acceleration (42).*y

Writing (63) as an acceleration DE (dividing by _{2} * l*) and inserting it the acceleration of (42), it is read (43).

Rearranging the terms and setting sin* ϕ* ≅

*, we get (44), which is a Mathieu type DE with parameters*ϕ

For generating parametric resonances, the (natural) system frequency needs to be coupled with the excitation frequency ω. Using the same nomenclature as in (64), we define the pendulum system frequency * ϕ*(

*) is growing exponentially, a frequency ratio*t

*:*ω

ω

_{2}= 1 : 1 is sufficient (as well as the ratio

*:*ω

ω

_{2}= 2 : 1), see also left-hand side of Figure 9. In case of letting the displacement term be harmonic (left-and right-hand side of Figures 9), the frequency ratio must be very close to a 2:1 ratio to have a large amplitude response. The ratio tolerance for having a large growth has a band width of ca. 1 rad to keep a large amplitude growth going.

Note that the response signal in orange on the left-hand side of Figure 10 is using an approximated linear displacement function sin(* ϕ*(

*)) →*t

*(*ϕ

*) and is scaled down by factor 10*t

^{−4}.The generated beat frequency signal is obtained using the exact harmonic displacement function.

### 3.1. Nonlinear single degree of freedom systems

Similar to the case in Section 2.1, also a SDoF system will be discussed, but this time a linear and a nonlinear spring will be present. The nonlinearity of this spring shall have the form shown in (3) having a nonlinear exponent * n* = 2 and

k

_{3}> 0 , a parameterization like that is generally used for a magnetic spring (see also top middle sketch in Figure 1 and appendix A.1). The lumped parameter model for the examined system (Figure 10) consists beside this spring system with linear spring rate

k

_{1}and nonlinear spring rate

k

_{3}of an oscillator mass

*, a viscous damping factor d and an external basepoint excited harmonic force with amplitude*m

y

_{0}(

*). System coordinate origin is placed at the basepoint excitation.*t

The elastic energy of this nonlinear spring system with * n* = 2 will lead to the following spring energy, see also derivation in (4).

Adding up all kinetic energies and all potential energies, disturbances in form of a viscous damping and a basepoint excited force is given in (46)–(49).

Applying the Lagrangian formalism (1) and (2), we deal again with a SDoF system, will lead to the resulting DE shown in (50), similar to the result derived in Section 2.1 – but here we have now introduced a nonlinear spring system.

Introducing dimensionless notation, by using a dimensionless time * τ* , a dimensionless system resonance frequency Ω, the damping factor

ξ

_{1}and the gravity offset term

*(51) and setting the dimensionless displacement*ϱ

*(52).*u

By replacing parameters of (50) with (51), (52) we obtain (53), including the dimensionless gravity term * ϱ* as well.

DE (53) is nonlinear, as we have also a path dependent function to the power of 3 and a dimensionless factor * β* which is generally small. If this factor

*is positive, we deal with a nonlinear spring hardening system, if*β

*is negative, it is a spring softening system. Unfortunately, such a system cannot be examined using the Laplace or Laplace-like transformation, such as [6], as this transformation can deal only with linear functions, respectively, nonlinear quadratic functions. For solving this nonlinear so-called Duffing DE, there are several methods available, such as averaging method or the harmonic balancing method. The averaging method assumes that a solution of the DE can be obtained using harmonic functions In Ref. [7], chapter 9.3, a general solution for nonlinearity terms with a positive integer exponent*β

*is obtained using the averaging method. Another method to get analytic solutions is, as said, the harmonic balance method, which is well explained in the textbook [8], chapter 2.3.4; the DE case of (53) is discussed in the same book, chapter 4.1 and there are many research papers to discuss this nonlinear DE, see for example [9, 10]. Note that this case is of nonlinear nature, as it includes the nonlinear term*βu

^{n}

βu

^{3}, but cannot exert parametric resonance. However, there are also many research papers where such nonlinearities coupled with a Mathieu DE are discussed, see for example [11].

The depicted Figure 11 shows simulation results of DE (53). The time domain behavior (top left) and its dimensionless phase space behavior (top right) is shown with simulating 50 periods with settled initial conditions (amplitude of * ω* = 5…40Hz and keeping the acceleration signal constant at 0.5 g.To make sure that only non-transient amplitudes are selected to create the frequency response, only in the last 5 periods (out of 50) the maximal and minimal value is selected). The top row is using a constant angular excitation of 25 Hz and depicts only one simulation point of the generated frequency response. The main simulation parameters are shown in the heading, a variant of this simulation is given in the appendix A.3. Note the shown simulated data are taken from a validated electromagnetic SDoF vibration energy harvester system by the author.

### 3.2. Nonlinear two degree of freedom systems

Let us consider the lumped parameter model in Figure 12. A pendulum with a stiff rod of length * l* and mass

m

_{2}suspended on a spring damper system with mass

m

_{1}and stiffness

*and damping factor*k

*. Mass*d

m

_{1}can move only in the depicted

*direction and pendulum only in the X-Y plane.*y

Governing equations are derived using again the Lagrange formalism. Considering the frame of reference at the origin shown in Figure 12 and defining in Cartesian coordinates first the two degrees of freedom vector _{y} and _{ϕ} (54).

The kinetic energy for both degrees of freedom are shown in (55, 56).

The potential energies derived from the same vectors lead to (57) and (58).

The Lagrange energy function * L* becomes:

The viscous friction for both degree of freedoms is given in (60) and the basepoint excited driving force is given in (61).

Applying the Lagrange formalism (2) for both degrees _{1} = * y* and

q

_{2}=

*lead to the coupled DE system of (62), (63), representing a force DE respectively a torque DE. On the right-hand side of those DE’s, the defined viscous frictions of (60) and the basepoint excitation (61) is present.*ϕ

The parameters for non dimensionalization are given in (64), (65). Note that the reference system frequency _{1} is set to the 1. DoF (also called primary, the mass spring system) and the second system frequency ω_{2} to the 2. DoF (the pendulum – also called secondary system). The excitation frequency is associated to * ω*.

Figure 13 depicts left the excitation (magenta) and the time response signals of * y* (blue) and

*(red) and its phase space behavior on the right-hand side.The parameters are chosen in such a way,that the pendulum starts rotate.The excitation primary system has a resonance ratio of*ϕ

*:*ω

ω

_{1}= 1 : 1 and secondary primary frequency ratio is

ω

_{2}:

ω

_{1}= 2 : 1. The stability of such a pendulum is described in chapter 4.4 of [7], where so called semi trivial and nontrivial solutions for this system are discussed.

A treaty of such a system, a kinetic energy harvesting device, with additionally a nonlinear spring system on the primary and an electromagnetic harvester on the secondary system is given in Ref. [12]. Note that the derivation of the system equations there have been made without the Lagrangian formalism and the found system equations are equivalent.

## 4. Conclusions

In the introduction, we showed the equivalence of rotary and translatory mechanical systems as well as the equivalence of mechanical and electrical resonance systems. Also, a brief introduction to the Lagrangian formalism is given. In preparation to nonlinear resonance systems, also rotational and translational springs are discussed. Three classes of spring systems have been identified: linear springs nonlinear symmetric springs and nonlinear asymmetric springs. Throughout the chapter further readings are proposed.

In Section 2, linear resonance systems with one and two degrees of freedom have been investigated using basepoint excited systems. Using the Laplace Transformation is most useful to analyze any linear resonance system with a periodic excitation.

Section 3 deals with nonlinear resonance systems. When in such a dynamical system one of the resulting DE’s is of Mathieu or Hill type, the response amplitude of such a system might grow exponentially. This is exemplary demonstrated in Section 3.2 identifying the system differential equations of a basepoint excited two degree of freedom system. Some dynamic properties of such a system is demonstrated.

## A. Appendix

## A.1. Nonlinear symmetric spring systems

Using instead of disk magnets ring magnets, strong nonlinearities can be generated. The following series in Figure A1 depicts a few simulation cases. Some of shown simulation cases have been validated and proven experimentally.

## A.2. Variant of linear 2DoF system

Instead of using linear springs, magnetic nonlinear springs can be used (see also a selection of such spring characteristics in A.1). Using nonlinear springs and making the system nonlinear (instead of having only a linear spring term, we have for each spring also a term of the form _{1} ^{3} and _{2} ^{3}). The relative response signals of such a system is depicted in Figure A2.

It is interesting, that the relative motion of the 2. DoF is responding with resonance between 19…25.5 Hz.The first degree of freedom has a nonlinear spring hardening behavior, reaching * A* = 7.9 mm at 25.5 Hz.

_{pp}

## A.3. Variant of nonlinear SDoF system

A created tool by the author in Matlab/Simulink has been used to simulate many basepoint excited SDoF or nDoF systems with rotational or translational or mixed structures.

It allows to simulate such systems with constant amplitude or constant acceleration, can handle hard or soft-impact of the oscillating proof mass(es). In addition, one sided spring characteristics can be simulated, see also Figure A3 – a feature that is especially interesting in relation with magnetic springs. Main disadvantage of such one-sided bound springs is the fact, that they need to be installed upright. The behavior of such a one-sided magnetic spring is depicted in Figure A3. It has a frequency response similar to a softening spring. The maximal amplitude of 3.65 mm occurs at 17.5 Hz (the two-sided classical hardening magnetic spring reaches an amplitude maximum of 4.3 mm at 27 Hz). Such a spring system could also be analytically described, by introducing for example continuous piecewise functions.

## References

- 1.
Feynman RP, Leighton RB, Sands ML. The Feynman Lectures on Physics. New York: Basic Books, 2011; 2010 - 2.
M. Cartmell, Introduction to linear, parametric and nonlinear vibrations. London, UK: Chapman and Hall, 1990 - 3.
Calkin MG. Lagrangian and Hamiltonian Mechanics: Hamiltonian and Lagrangian Formalism. Singapore, River Edge NJ: World Scientific; 1996 - 4.
Pesce CP. The application of Lagrange equations to mechanical systems with mass explicitly dependent on position. Journal of Applied Mechanics. 2003; 70 (5):751 - 5.
Spreemann D, Manoli Y. Electromagnetic Vibration Energy Harvesting Devices: Architectures, Design, Modeling and Optimization. Dordrecht: Springer Netherlands; 2012 - 6.
Elzaki TM. On the New Integral Transform “ELzaki Transform” Fundamental Properties Investigations and Applications. Global Journal of Mathematical Sciences. 2012; 4 (1):1-13 - 7.
Tondl A. Autoparametric Resonance in Mechanical Systems. Cambridge, UK: Cambridge University Press, 2000 - 8.
Nayfeh AH, Mook DT. Nonlinear Oscillations. Weinheim: Wiley; 2007 - 9.
Kurmann L et al. Autoparametric resonance Systems for vibration-based energy harvesters. Journal of Physics: Conference Series. 2015; 660 :12070 - 10.
Mann BP, Sims ND. Energy harvesting from the nonlinear oscillations of magnetic levitation. Journal of Sound and Vibration. 2009; 319 (1–2):515-530 - 11.
A. C. Luo and D. O’Connor, Stable and unstable periodic solutions to the mathieu-duffing oscillator: 6 - 11 Aug. 2012, Budapest, Hungary ; [proceedings]. Piscataway, NJ: IEEE, 2012 - 12.
Kurmann L, Jia Y, Manoli Y, Woias P. Magnetically levitated autoparametric broadband vibration energy harvesting. Journal of Physics: Conference Series. 2016; 773 :12006

## Notes

- A holonomic constraint depends only on the coordinates and time and does not depend on velocities.