Open access peer-reviewed chapter

Introduction to Parametric and Autoparametric Resonance

Written By

Lukas C. Kurmann

Submitted: 12 April 2017 Reviewed: 04 August 2017 Published: 29 November 2017

DOI: 10.5772/intechopen.70521

From the Edited Volume

Resonance

Edited by Jan Awrejcewicz

Chapter metrics overview

2,554 Chapter Downloads

View Full Metrics

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].

TranslationalRotational
SymbolDescriptionSI UnitSymbolDescriptionSI Unit
sDisplacementmϕAnglerad
v=dsdtVelocitymsω=ϕddtAngular velocityrads
a=dvdtms2α=ωddtAngular accelerationrads2
mMasskgJInertiakg m2
F = m aForceNT = J αTorqueNm
I = m vMomentumNsL = J ωAngular momentumNms
T=12mv2Kinetic energyNmT=12Jω2Kinetic EnergyNm
U=12ky2Potential energyNmU=12cϕ2Potential energyNm
W =  ∫ FdsWorkJW =  ∫ TdϕWorkJ
P = F vPowerWP = J ωAngular powerW

Table 1.

Comparison of translational and rotational motion parameter characteristics.

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 systems1. 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 T and the potential energy U.

L=TUE1

In (2) the Lagrangian equation is given with generalized coordinates qi of a dynamic system and dissipative generalized forces Qi.

ddtLqi̇Lqi=QiE2

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 U 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.

Looking at translational (classic) mechanical springs, the displacement dependent force F(qt) can be written as shown in (3) using spring stiffness k1Nm, the generalized translational coordinates qt and having also introduced a nonlinear spring term kn and an exponent 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.

Fqt=k1qt+knqtnE3

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

Eelastics=0sFqtdqt=0sk1qt+knqtn=12k1+kn1+nqt1+ns0
Eelastics=12k1s2+kn1+ns1+nwithn>1E4

The elastic energy Eelastic (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 D1Nmrad, the generalized translational coordinates qr and having also introduced a nonlinear torsion spring term Dn.

Tqr=D1qr+DnqrnE5

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

Eelasticϕ=0ϕTqrdqr=0ϕD1qr+Dnqrndqr=12D1qr2+kn1+nDnqr1+nϕ0
Eelasticϕ=12D1ϕ2+11+nDnϕ1+nwithn>1E6

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 kn > 0 – forming so called hardening springs if kn < 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.

Figure 1.

Sketches of translational and rotational spring systems.

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.

Figure 2.

Exemplarily displacement-force signals (l) and angular displacement-torque signals (r) of the shown spring systems of Figure 1.

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 q and current i on the right-hand side. Kinetic energies are denoted with TT (translational kinetic energy) and TL (inductive kinetic energy), potentials are written as US (spring potential) and UC (capacitive potential) and the gravitational potential denoted as UG and the DC battery voltage is UB. Non-conservative components are in the mechanical system the viscous damping force and in the electrical system the electrical resistor.

Table 2.

Equivalence of electrical and mechanical systems.

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

Lmech=TTUS+UG=12my2+12ky2mɡyE7
Lel=TLUC+UB=12Lq2+121Cq2U0qE8

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 q 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.

my+dy+ky=ɡmE9
Lq+Rq+1Cq=U0E10

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 ωmech2=km and for the electrical system (10) ωel2=1LC.

Advertisement

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 d, a linear spring with a spring rate k1 and an external basepoint excited harmonic force with amplitude y0(t). System coordinate origin is placed at the basepoint excitation.

Figure 3.

Linear single degree of freedom (SDoF) spring mass damper model of a resonant harmonic basepoint excited oscillator.

The kinetic energy (11) of this system and the potential energy (12) form the non-dissipative energy of this system. The dissipative force Fd of this system – we consider only viscous friction – is shown in (13) and the driving force F0 (14) assuming a harmonic basepoint excitation with amplitude A and driving frequency ω.

T=12my2E11
U=mɡy+12k1y2E12
Fd=dyE13
F0=md2dt2Acosωt=mAω2cosωtE14

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 F0, as it is a harmonic signal, can be introduced with a positive or a negative sign, resulting in a phase shift of 180°.

myk1y=FdF0E15
my+dy+k1y+=mAω2cosωtE16

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 u (18).

τ=tω1;Ω=ωω1;ω12=k1m;ξ1=d2mω1;ϱ=ɡAω12E17
path uτ=ytAE18

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 u (20).

u+2ξ1u+u+ϱ=Ω2cosΩτE19
u+2ξ1u+u=Ω2cosΩτE20

Frequency domain behavior is obtained by applying the Laplace Transformation (2124) and by replacing s = jΩ we obtain the frequency response (25).

Us=Luτ=0uτedτE21
Lu+2ξ1u+u=d2dt2cosΩτE22
Uss2+2ξ1Uss+Us=s2Y0sE23
Gs=UsY0s=s21+2ξ1s+s2E24
Gξ1Ω=UjΩY0jΩ=Ω21+Ω2+j2ξ1ΩE25

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 Ω.

Figure 4.

Representation of frequency response of a linear SDoF system using (24) Bode diagram (left) and absolute value representation of (25) (right).

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 m1 and m2, damping factors d1 and d2, linear springs with spring rates k1 and k2 and an external basepoint excited harmonic force with amplitude y0(t). System coordinate origin is placed at the basepoint excitation.

Figure 5.

Linear 2DoF spring mass damper model of a resonant harmonic basepoint excited oscillator.

T=12m1y12+12m2y22E26
U=m1ɡy1+m2ɡy2+12k1y12+12k2y2y12E27
Fd=d1y1+d2y2E28
F0=m1+m2d2dt2Acosωt=m1+m2Aω2cosωtE29

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

m1y1+d1y1+k1y1k2y2y1+m1g=m1+m2Aω2cosωtE30
m2y2+d2y2+k2y2y1+m2ɡ=0E31

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

uτ+2ξ1uτ+uτ+λmΩ02uτλmΩ02vτ+ϱ=Ω2cosτΩE32
vτ+2ξ2Ω0vτ+Ω02vτΩ02uτ+ϱ=0E33

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 u and v.

τ=tω1;Ω=ωω1;ω12=k1m1;ω22=k2m2;ξ1=d12m1ω1;ξ2=d22m2ω2;Ω0=ω2ω1;λm=m2m1;ϱ=ɡAω12E34
path uτ=y1tA and path vτ=y2tAE35

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).

Uss2+2ξ1Uss+Us1+λmΩ02λmΩ02Vs+ϱ=s2Y0sE36
Vss2+2ξ2Ω0Vss+VsΩ02UsΩ02+ϱ=0E37

As (36) and (37) represent two algebraic equations, U(s) and V(s) can be separated, resulting in (38) and (39). Note that the gravity term ϱ in the numerator will introduce an additional damping of the transfer function.

Us=ρλmΩ02+s2+2sξ2Ω0+Ω02ρ+s2Y0sλmΩ04s2+2sξ2Ω0+Ω021+s2+2sξ1+λmΩ02E38
Vs=ρ1+s2+2sξ1+1+λmΩ02+s2Ω02Y0ss4+Ω02+2s3ξ1+ξ2Ω0+2sΩ0ξ2+ξ1Ω0+λmξ2Ω02+s21+Ω04ξ1ξ2+Ω0+λmΩ0E39

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, λm = 0.42 and Ω0 = 1. As we have two resonators with same system frequencies ω1=ω2=169rads, two resonances will occur.This system reaches resonances at 0.71 Ω (19Hz) and at 1.41 Ω (38Hz) for Ω0 = 1 and 0.82 Ω and at 6.1 Ω for Ω0 = 5 (dashed lines).

Figure 6.

Bode diagram of a linear 2DoF system represented by (38) and (39) (left); first oscillator with mass m1 in blue and second with m2 in red with Ω0 = 1 and Ω0 = 5 (dashed lines) and constructed frequency response using time domain signals (right).

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 ŷ1start<1μmandŷ2start<1μm). The main simulation parameters are shown in the heading, a variant of this setup using nonlinear springs is given in the appendix A.3.

Figure 7.

Numerical simulation results of the linear basepoint excited 2DoF system shown in (32) and (33) using a constant acceleration of 0.5g and a basepoint excitation of ω = 25 Hz.

Advertisement

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.

Figure 8.

Examples of physical systems exhibiting potential parametric resonance effects, adapted after [2].

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).

q+qa+bcosΩt=0E40

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

q+qa+fpt=0E41

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 y0 the primary system y, compare also the lumped parameter model in Figure 12 (in an experiment we would simply make the stiffness k of the system very large, for example, replacing the spring with a fixed stiff rod). Now the new induced basepoint excitation y 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=Acosωty=Aω2cosωtE42

Writing (63) as an acceleration DE (dividing by m2 l) and inserting it the acceleration of (42), it is read (43).

lϕ+Dm2lϕ+ɡsinϕ+Aω2cosωtsinϕ=0E43

Rearranging the terms and setting sinϕ ≅ ϕ, we get (44), which is a Mathieu type DE with parameters a=gl and b=Aω2l.

ϕ+Dm2ϕ+sinϕɡl+Aω2lcosωt=0E44

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 ω22=gl. For generating parametric resonances for which the angle ϕ(t) is growing exponentially, a frequency ratio ω : ω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.

Figure 9.

Nonlinear single degree of freedom (SDoF) spring mass damper model of a resonant harmonic basepoint excited oscillator.

Figure 10.

Response signals of a parametrically excited pendulum examining DE (47) with keeping damping term D = 0, l = 108.1 mm, A = 100 mm and ω1 = 10 rad.

Note that the response signal in orange on the left-hand side of Figure 10 is using an approximated linear displacement function sin(ϕ(t)) → ϕ(t) and is scaled down by factor 10−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 k3 > 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 k1 and nonlinear spring rate k3 of an oscillator mass m, a viscous damping factor d and an external basepoint excited harmonic force with amplitude y0(t). System coordinate origin is placed at the basepoint excitation.

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

Eelasticy=12k1y2+14k3y4E45

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).

T=12my2E46
U=mgy+12k1y2+14k3y4E47
Fd=dyE48
F0=md2dt2Acosωt=mAω2cosωtE49

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.

my+dy+k1y+k3y3+=mAω2cosωtE50

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 u (52).

τ=tω1;Ω=ωω1;ω12=k1m;β=k3k1A2;ξ1=d2mω1;ϱ=ɡAω12E51
path uτ=ytAE52

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

u+2ξ1u+u+βu3+ϱ=Ω2cosΩτE53

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 βun 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 βu3, 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 ŷ<1μm). The top row depicts one simulation point in the bottom row, where a frequency sweep has been done, sweeping the basepoint excitation from ω = 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.

Figure 11.

Numerical simulation results of the nonlinear basepoint excited SDoF system shown in (53) using a constant acceleration of 0.5 g and a basepoint excitation of ω = 25 Hz (top row) and its sweep behavior ω = 5…40 Hz (bottom row).

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 m2 suspended on a spring damper system with mass m1 and stiffness k and damping factor d. Mass m1 can move only in the depicted y direction and pendulum only in the X-Y plane.

Figure 12.

Nonlinear two degree of freedom (2DoF) spring mass damper model.

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 ry and rϕ (54).

ry=0yt and rϕ=lsinϕtytlcosϕtE54

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

Ty=12m1ddtry.x2+ddtry.y2=12m1y2E55
Tϕ=12m2ddtrϕ.x2+ddtrϕ.y2=12m2l2cosϕ2ϕ2+12m2y+lsinϕϕ2E56

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

Uy=m1ɡy+12ky2E57
Uϕ=m2ɡylcosϕE58

The Lagrange energy function L becomes:

L=Ty+Tϕ(Uy+Uϕ)E59

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

Fdy=dy and T=DϕE60
F0=m1+m2d2dt2Acosωt=m1+m2Aω2cosωtE61

Applying the Lagrange formalism (2) for both degrees q1 = y and q2 = ϕ 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.

m1+m2y+ɡm1+m2+ky+lm2ϕ2cosϕ+ϕsinϕ=dy+F0E62
m2l2ϕ+ɡlm2sinϕ+lm2ysinϕ=DϕE63

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 ω.

τ=tω1;λm=m2m1+m2=m2m;Ω=ωω1;ω12=km;ω22=ɡl;Ω02=ω22ω12;λl=lA;
ξ1=d2m1ω1;ξ2=D2m2l2ω2;ϱ=ɡAω12E64
path uτ=ytA and angle θτ=ϕtϕ0E65
u+2ξ1u+u+λmλlθ2cosθ+θsinθ+ϱ=Ω2cosΩτE66
θ+2ξ2Ω0θ+sinθΩ02+λl1u=0E67

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.

Figure 13.

Numerical simulation results of the nonlinear basepoint excited 2DoF system shown in (66, 67) using a constant acceleration of 0.2 g and a constant basepoint excitation of fexc = 22.22 Hz.

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.

Advertisement

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.

Advertisement

A. Appendix

Advertisement

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.

Figure A1.

Spring force behavior of ring magnets using different distances of non-movable magnets.

Figure A2.

Response signals of a 2DoF system using lumped parameter model in Figure 5: Instead of having linear springs, also nonlinear springs are present.

Advertisement

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 u3 and β2 v3). The relative response signals of such a system is depicted in Figure A2.

Figure A3.

Response signals of a one-sided bound magnetic spring.

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 App = 7.9 mm at 25.5 Hz.

Advertisement

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. 1. Feynman RP, Leighton RB, Sands ML. The Feynman Lectures on Physics. New York: Basic Books, 2011; 2010
  2. 2. M. Cartmell, Introduction to linear, parametric and nonlinear vibrations. London, UK: Chapman and Hall, 1990
  3. 3. Calkin MG. Lagrangian and Hamiltonian Mechanics: Hamiltonian and Lagrangian Formalism. Singapore, River Edge NJ: World Scientific; 1996
  4. 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. 5. Spreemann D, Manoli Y. Electromagnetic Vibration Energy Harvesting Devices: Architectures, Design, Modeling and Optimization. Dordrecht: Springer Netherlands; 2012
  6. 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. 7. Tondl A. Autoparametric Resonance in Mechanical Systems. Cambridge, UK: Cambridge University Press, 2000
  8. 8. Nayfeh AH, Mook DT. Nonlinear Oscillations. Weinheim: Wiley; 2007
  9. 9. Kurmann L et al. Autoparametric resonance Systems for vibration-based energy harvesters. Journal of Physics: Conference Series. 2015;660:12070
  10. 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. 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. 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.

Written By

Lukas C. Kurmann

Submitted: 12 April 2017 Reviewed: 04 August 2017 Published: 29 November 2017