Open access peer-reviewed chapter

Simple Oscillating Systems

By Ivo Čáp, Klára Čápová, Milan Smetana and Štefan Borik

Reviewed: November 16th 2021Published: December 24th 2021

DOI: 10.5772/intechopen.101649

Downloaded: 37

1. Introduction

A wave is a disturbance (deviation from equilibrium) that propagates through space. This disturbance can be caused by an impulse excitation (shock wave) or by a time-varying excitation (e.g., the sound generated by vocal cords). The most frequent cause of wave excitation is a source with periodical time dependence—oscillations. The simplest one is harmonic oscillation. For example, we can express any periodic function of time as a sum of the harmonic functions (Fourier series). Thus, the excitation by the harmonic oscillations is a matter of specific interest. In this chapter, attention is, therefore, paid to the description of the physical nature of oscillations and their properties, see also Halliday [1].

Oscillations represent a very wide group of processes, which are generally characterised by their regular state repeating caused by the internal dynamics of a system. Such systems, whose internal couplings allow oscillations, are called oscillating systems. From the energy point of view, the oscillations are conditioned by the existence of two conservative forms of energy, which can reversiblyexchange due to the internal dynamics of the system. There is, for example, potential energy—kinetic energy (oscillations of mass on a spring) or electric field energy of capacitor—magnetic field energy of inductor (an oscillating LC circuit). A special case represents the ‘oscillations’ in a rotating system, such as the movement of a conical pendulum where energy exchanges between two perpendicular kinetic components of 12mvx2and 12mvy2, or the precession of a rotating body where energy exchanges between two perpendicular rotational components of kinetic energy 12Jxωx2and 12Jyωy2. Some of these cases will be described below as examples.

If the oscillating system is isolated from external influences, it oscillates spontaneously after the initial energy supply (excitation). Thus, we are talking about self-sustained oscillations. The oscillation amplitude remains almost constant if the energy losses of the oscillations in the system are negligibly small. The oscillations of the ideal lossless system are called undamped self-oscillationsand represent only theoretical idealisation. There exist loss mechanisms in each real system. They cause the irreversibletransformation of the conservative form of the system energy into another non-conservativeone, for example, friction, heat losses due to internal friction, energy losses of the electrical system by radiation, etc. As the total energy of the conservative components decreases, the amplitude of oscillations gradually decreases over time too. We call them damped self-oscillations. Due to the damping, these self-oscillations disappear after some time. Thus, they are a transient phenomenon in the system, for example, the vibrations of the string of the musical instrument fade; a swinging of pendulum stops after a certain time; oscillations of an LC circuit gradually disappear, etc.

The system may oscillate permanently without damping if there is a mechanism capable to cover the energy losses from an energy storage device. Such systems are different types of oscillators. Examples are pendulum clocks with weights or watches with a spring. In modern watches, a precisely sharpened crystal represents the oscillating system. In this case, a small lithium cell covers the energy losses caused by damping.

Specific phenomena arise when the system is exposed to periodic force. When such a force acts, the system, after attenuating the transient event, enters a steady state, characterised by oscillations with a constant amplitude and a period equal to the excitation period. These oscillations are called forced oscillationsof the system. The magnitude of the response to periodic excitation depends on the period or frequency of excitation. Significant is the resonancephenomenon that occurs when the excitation frequency is equal to the frequency of the system’s undamped oscillations.

The following sections focus on the different types of oscillations in simple systems, that is, in systems in which two conservative forms of energy occur in oscillations. This textbook presents a summary of the knowledge with an emphasis on application. A more detailed analysis of the mentioned phenomena can be found in physics textbooks.

1.1 Undamped self-oscillations

As the basic model of the oscillating system, we use a particle bound to the equilibrium position by the reversing conservative force of the springs (Figure 1).

Figure 1.

Particle bound to the equilibrium by the restoring force.

At the top of the figure, the particle is in equilibrium, and the resulting force acting on it is zero. If the particle shifts from the equilibrium position by the displacement of x, there arises a force of F(x) which depends on the xdisplacement, is reversible, and has the opposite direction as the displacement.

If the particle is displaced from the equilibrium position and released, it starts to move back to the equilibrium position. Its velocity is a derivative of the displacement

v=dxdt=ẋadotoverxdescribes the time derivative.E1

Displacing particle from equilibrium by x, we perform a work of W, which represents the potential energy of the particle

Epx=W=0xFξ.E2

A moving particle has a kinetic energy

Ekv=12mv2.E3

Reversible energy exchange occurs between the energy components of Epand Ek. It is a system with two degrees of freedom. If we consider no loss mechanisms, the sum of the Ep + Ekenergy components remains constant. The situation is in Figure 2. It indicates the potential energy as a function of the displacement x(solid line). We can see that in the equilibrium position (x = 0) the potential energy is minimal, and thus, at a constant sum, Etot = Ep + Ek, the kinetic energy is maximal. It means that the particle moves due to inertia until its kinetic energy drops to zero. The particle thus moves periodically between the extreme positions of Aand B, which are given by the total energy Etot.

Figure 2.

Total potential and kinetic energy of a particle moving along thex-axis under conservative forces (law of conservation of mechanical energy).

We can express a function of the potential energy Ep(x) near the minimum, that is, equilibrium position, by the Taylor power series

Epx=12!kx2+13!lx3+14!nx4+,E4

where k=d2Epdx2|x=0represents the stiffnessof the system, l=d3Epdx3|x=0expresses asymmetryof the potential energy function regarding the equilibrium position. Other coefficients such as n=d4Epdx4|x=0and higher (odd and even derivatives) have similar characteristics but they change function course in the larger distance of xfrom the equilibrium position. The first power term of xis zero because it is the local minimum of the Ep(x).

A negative potential energy gradient defines the force acting on a particle as follows

Fx=dEpxdx=kx12!lx213!nx3.E5

We distinguish the linearand non-linearoscillating systems depending on the number of force or energy terms that are considered in the motion. If the displacement xapproaches zero, the higher powers of xndecrease faster than the first one. Then, the higher terms of the function F(x) are negligibly small, and the system appears to be linear. However, the terms with higher powers apply, and the system behaves as non-linear if displacement x significantly increases.

1.1.1 Undamped self-oscillations of linear system

A system is linearif the restoring force is a linear function of the displacement of xfrom the equilibrium position. According to Eqs. (4) or (5), it follows

Fxkx,orEpx12kx2.E6

Potential energy is a quadratic function of the displacement xand is called a quadratic potential well. Its graph is a quadratic parabola. Figure 3 shows a replacement of the real function Ep(x) by a quadratic function. This replacement fits well only in the near vicinity of the minimum, that is, only for small variations of xaround the equilibrium. Oscillations within the range of the fitted region are sometimes called small oscillations.

Figure 3.

Quadratic fitting (dashed) of the potential energy function.

Equation of motion of the particle ma = F, where a=x¨is an acceleration, has form for the linear system as follows

x¨+ω02x=0,whereω0=km.E7

The solution of this equation is the function

xt=xmsinω0t+α,E8

where xmand αare integration constants and their values are determined from initial conditions x(0) = x0 and v0=ẋ0=v0at t = 0.

It is a harmonic motion where the xmis an amplitude of oscillations and αis a phase constant (initial phase).

Oscillations with harmonic time dependence are called harmonic oscillations. It follows from the previous description that harmonic oscillations occur when a particle (body) moves in a quadratic potential well.

The quantity

ω0=2πf0=2πT0E9

represents the angular frequency, where f0 is frequency and T0 is the period of the undamped self-oscillationsof the system. According to Eq. (7), these quantities depend on the stiffness kof the system, and the inertia given by the mass mof the particle. As the stiffness increases, the frequency f0 increases as well, and the period T0 decreases. With the mass increase, the frequency f0 decreases, and the period T0 increases. For example, as a body hung on the spring oscillates with the period of order seconds, an atom in the crystal lattice with the period of the order of 10−14 s.

Example 1. Oscillations in an electrical LC circuit.

Let us assume one loop electrical circuit consisting of an inductor Land acapacitor C. Electrical current i(t) flows through this circuit, and we can express the energy of an electrical field of the capacitor and a magnetic field of the inductor as follows:

E=121CQ2+12Li2,

where the electrical charge Qof the capacitor relates to the current iof the inductor i=Q̇. We can see the analogy between electrical and mechanical systems, in case x → Q, k → 1/Cand m → L. If we do not consider the power losses, the energy Eis constant, and by differentiating it, we get the equation

1CQdQdt+Lididt=0.

Dividing the equation by i = dQ/dt, we get

d2Qdt2+1LCQ=0,

which has the same form as (7). Then, the solution is

Qt=Qmsinω0t+α,whereω0=1LC.

The capacitor voltage is

uCt=QtC=QmCsinω0t+α=Umsinω0t+α

and the inductor current equals to

it=dQtdt=ω0Qmcosω0t+α=Imsinω0t+α+π2.

Thus, there are the harmonic undamped oscillations of the circuit quantities with the angular frequency of ω0.

Example 2. Pendulum.

Consider a small body suspended on a long fibre (Figure 4). After the initial excitation, the body oscillates around the equilibrium position, and thus performs a circular motion with a radius equal to the fibre length of the l. If we displace the fibre by an angle φfrom the equilibrium position, then the potential energy of the body changes as

Figure 4.

Pendulum.

Epφ=mgh=mgl1cosφ.

The Ep(φ) function is not quadratic, and therefore, we can use a decomposition using the power series

Epφ=mglφ22!φ44!+.

We can neglect the series terms of the higher order for φ < < 1, and then the potential energy is

Epφ=12mglφ2.

Additionally, the kinetic energy is

Ekv=12mv2=12ml2φ̇2.

We can see the analogy to (7) again, if x → φ, k → mgland m → ml2.

Pendulum displacement is described by the function

ϕt=ϕmsinω0t+α,whereω0=gl.

Hence, the body oscillates with the period

T0=2πlg.

By measuring the oscillation period, it is possible to determine the length of the pendulum if we do not have a measuring tool. Alternatively, with a pendulum of a certain length, we can realise a periodic movement with the required period, such as in the case of the pendulum clock.

Example 3. Cone pendulum.

Consider the same case as in the previous example, but let the body move along a circle in the horizontal plane (x, y). The pendulum copies a conical surface as it moves. Thus, the deviation angle from the vertical axis is φ, as shown in Figure 5. The radius of motion of the body is R = lsin φ. The kinetic energy of the body is

Figure 5.

Cone pendulum.

Ek=12mv2=12mvx2+vy2,

where vis the velocity of the circular motion.

The force acts on the body and equals to Fd = −2r, where the angular velocity is ω = v/r. The centrifugal force composes of gravitational force Fg = mg, while the resultant force has the direction of the pendulum fibre that means tgφ = Fd/Fg. For small displacement, when φ < < 1 rad, it follows that tgφ≈ sinφ = r/l. From the Fd/Fg = r/l, we obtain ω= g/l, which is similar to the previous example.

Potential energy connected with the centrifugal force is given as

Ep=12mω2r2=12mω2x2+y2.

And finally, the total energy can be expressed

Ek+Ep=12mvx2+12mω2x2+12mvy2+12mω2y2.

In this case, the motion can be considered as a superposition of two mutually perpendicular oscillations in the x- and y-direction, which are phase-shifted by π/2 rad. The total energy (Ek + Ep) is constant and is the sum of the total energy of oscillations in the x- and y-direction.

Example 4. Precession of magnetic dipole in the magnetic field.

Another specific case of periodic movement is the precession (rotating axis of a rotating body). We can observe it looking at a children’s toy, such as a spinning top. By spinning and laying it on the pad, the toy axis rotates, see the illustration. The precession occurs due to the gravitational force.

Similarly, the magnetic dipole, here the proton, is affected by an external magnetic field. One of the proton parameters is the angular momentum L, which describes its mechanical rotation. The rotation of the charged particle is associated with the accompanying magnetic field. Thus, the proton behaves like an elemental magnet (the magnetic dipole) with a magnetic moment m. The ratio of magnetic moment to mechanical angular momentum is called the gyromagnetic ratio γ = m/L(see Table 1). If the dipole is in an external magnetic field, then the moment of the force acting on it is.

Cores(spin)γ[×108 s−1⋅T−1]
1H (proton)1/22.68
13C1/20.67
19F1/22.52
31P1/21.08
Free electron1/2−1758

Table 1.

Properties of selected nuclei of atoms and electron.

M=m×B, where Bis the magnetic induction.

The moment of force determines the dynamics of the dipole movement. The basic equation of rotational motion (impulse theorem II) has the form

dLdt=M.

Combining both equations, we get dL=m×Bdt=γL×Bdt.

The dLvector is perpendicular to the vector L, and therefore its magnitude does not change but the direction only.

As shown in Figure 6, the end of the Lvector moves along a circle with a radius equal to Lsinα. The angle dφover the time dtdetermines the magnitude of change of dL = Lsinαdφ. The magnitude of the dLchange according to the equation of motion is dL = mBsinαdt. By comparing these two expressions, we get the angular velocity of the endpoint of the Lvector

Figure 6.

Precession movement of a magnetic dipole in a magnetic field.

ωL=dt=γB.

Thus, the dipole axis performs a circular (funnel) motion in the magnetic field, called the Larmor’s precession. The frequency of fL = ωL/2π of this motion depends on the type of particle represented by its gyromagnetic ratio γand the induction Bof the magnetic field but does not depend on the angle α. As we show later, this phenomenon is used in magnetic resonance imaging and magnetic resonance spectroscopy. The nature of the phenomenon is like that of a conical pendulum.

Example 5. Ion oscillations in the crystal.

Crystals represent a simple or more complex regular arrangement of atoms of solids. For example, aluminium consists of an arranged lattice of positive ions. There are Al+ ions and electron gas. The ions are subjected to electric forces by the surrounding particles. The equilibrium ion position is given by the zero resultant force or by the minimum value of the potential energy. If the ion deviates from the equilibrium position, it begins to oscillate around it.

As a simple model, consider three monovalent ions, of which two are fixed, and the third can move between the other two. In equilibrium, the distance of the central ion from the extreme ones is a(see Figure 7).

Figure 7.

Three monovalent ions.

Let us move the central ion displacing it from the equilibrium. Then, the force acting on the ion is

F=F1F2=ke2a+xke2ax=k2ae2a2x2x,

where k ≈ 9.0 × 109 m·F−1 is Coulomb’s law constant and e ≈ 1.6 × 10−19 C is the elementary charge.

If the displacement is x < < a, then we can express the resultant force by the linear approximation as

F2ke2ax=Kx.

As shown, if a particle with a mass mexerted by a reversing force proportional to the displacement x, the particle oscillates around an equilibrium position with a frequency

f=12πKm=12π2ke2ma.

The oscillating of a charged particle is the source of the electromagnetic wave at this frequency and the wavelength of this wave is

λ=cf=2πcma2ke2.

For example, if we use the typical values for aluminium: m ≈ 4.5 × 10−26 kg, a ≈ 2.8 × 10−10 m, while c ≈ 3.0 × 108 m·s−1, we get λ ≈ 31 μm.

The result corresponds to the wavelength of infrared (thermal) radiation. Oscillations of crystal lattice ions are the cause of the thermal radiation of the bodies.

There are many similar examples of oscillating systems, all of which have a similar physical nature. It is always a periodic exchange of energy between the various conservative forms of energy caused by the internal dynamics of the system.

1.1.2 Undamped self-oscillations of non-linear system

We find the system as non-linear if we cannot neglect its non-linearity. This means that we consider other higher terms in the expression of force by the power series [Eq. (5)]. Since the terms of the series generally gradually decrease with an increasing exponent of power, we can now consider the first higher non-zero member only. If the Ep(x) function is odd (it means asymmetric potential well), we consider the term with the lcoefficient, that is., quadratic term in the force expression. If the potential well is symmetrical, it is l = 0, the first non-linear term of the series is a cubic one. Accordingly, we are solving single cases by using this simplification.

In the following section, we analyse the case of oscillations in an asymmetric potential well, for which we express the force acting on a particle in the form

Fxkx12!lx2.E10

Then the equation of motion is

ma=kx12lx2.E11

We can rewrite the equation to the form

x¨+ω021+λxx=0,E12

where ω02=km.

The coefficient λ=l2kis the degree of asymmetry of the potential well.

Figure 8 shows an example of the asymmetric potential well and it illustrates the fitting of the well by a quadratic function (dashed line). This function fits the well only in the near vicinity of the equilibrium position. Cubic function correction is positive on the left side and negative on the right side, which means that the asymmetry coefficient is l < 0.

Figure 8.

Asymmetric potential well with asymmetry of the typel < 0 orλ < 0.

Equation (12) represents a non-linear differential equation. When solving it, we use the physical nature of the phenomenon, which means the particle motion is periodic with an unknown angular frequency ω. We know, the periodic function can be expressed in the form of a Fourier series. If we choose for the start time t = 0 the moment when the particle displacement crosses the extreme value, then we can describe the course of the time dependence as an even function (symmetrical around the beginning t = 0).

For an even function, the Fourier series contains only even (cosine) terms.

xt=a0+n=1ancosnωt.E13

The a0 value represents the mean value of the particle displacement, the anare amplitudes of the individual harmonics with frequencies of nω.

The solution procedure is such that we substitute the function (13) into the differential equation and arrange the terms according to the angular frequency. If the expression on the left side is to be equal to the right side of the equation (i.e., zero), all terms must be zero at corresponding frequencies—harmonics with angular frequencies ω, 2ω, etc. So, we get a set of equations for unknown parameters ω, a0, and anfor n = 1, 2, …

Assuming the weak non-linearity of the system, which is given by λ xm ≪ 1, where xmis the maximal displacement from the equilibrium position, we get

ω=ω0156λa12,a012λa1a1,a216λa1a1etc.,E14

where a1 is the amplitude of the first harmonic with the frequency of ω. Thus, non-linearity influences the frequency of the self-oscillations. It causes the shift of the mean value of the position a0, and it causes the higher harmonics involved in oscillations.

Example of the derivation:

After substituting into the differential equation, we get the equation

n=1ann2ω2cosnωt+ω02a0+ω02n=1ancosnωt+ω02λa0+n=1ancosnωta0+k=1akcoskωt=0

and then

ω02a0+ω02λa02+n=1ω02n2ω2ancosnωt+2ω02λa0n=1ancosnωt++ω02λn=1k=1akancoskωtcosnωt=0

and then

ω02a01+λa0+n=1ω02n2ω2+2ω02λa0ancosnωt++ω02λn=1k=112akancosk+nωt+cosknωt=0

Equality must be met separately for each harmonic component and for the constant component.

For constant terms of the equation, we have

ω02a01+λa0+12ω02λn=1an2=0.

For terms with a fundamental angular frequency ω, we get the equation

ω02ω2+2ω02λa0a1+ω02λn=1anan+1=0.

Then, the terms with the frequency of 2ω(second harmonic)

ω024ω2+2ω02λa0a2+12ω02λa12+ω02λn=1anan+2=0,etc.

From the second equation, we get after neglecting higher terms

ω02ω2+2ω02λa0+ω02λa2a1=0,

from where

ω2=ω021+λ2a0+a2ω02.

It yields from the first equation

a01+λa0+12λa12=0,and approximatelya012λa12.

From the third equation, we obtain

ω024ω2a2+12ω02λa12=0anda216λa12.

By substituting a0 and a2 to relation for ω2, we get a more precise result in the form

ω2=ω021+λ2a0+a2=ω02156λa12.

Example 6. Thermal expansion of substances.

Atoms or molecules of solids or liquids are arranged in ordered structures. Attractive electric forces ensure the consistency of the substance. Approaching or moving the molecules or atoms together causes repulsive forces, which, along with attractive forces, provide equilibrium distances. The potential energy of the particle relative to the adjacent particle is shown in the figure. We can see that the potential well is asymmetrical. The minimum potential energy corresponds to the equilibrium distance of the particles of the substance. If we supply the particles with energy (e.g., in the form of heat), the amplitude of the oscillations of the particles increases.

Moreover, due to the non-linearity of the binding potential, the mean interatomic distance also increases. It means the macroscopic elongation of the material. According to the Eq. (14), the displacement of the mean distance of a0 is proportional to the square a12 of the amplitude of the fundamental harmonic. This amplitude square is proportional to the energy of the oscillations and the temperature. Hence, the thermal expansion of the substances is

Δll0=αTT0,resp.l=l01+αTT0,

where αis the coefficient of the length thermal expansion.

1.2 Oscillations in the linear system with viscous damping

In real systems, oscillation damping occurs because of irreversible energy loss of the system during the oscillation process. The loss mechanism describes the force that depends on the movement state of the system. In mechanical systems, it is mainly friction or resistance of the environment. In electrical circuits, there are Joule losses when current is passing through a resistor or emitting EM waves to the surrounding space. Figure 9 shows an example of a damped oscillation model.

Figure 9.

Comparison of damped (a) and undamped oscillations (b).

Let us consider the loss mechanism that often occurs in oscillating systems, which is a viscous resistance. A resistive force proportional to the velocity of movement characterises it, or in other words, the viscous resistance depends on power dissipation proportional to the square of the velocity (in electric circuits, it is the square of the current). Thus,

Fo=rvE15

where ris the coefficient of resistance. In the case of mechanical resistance, the viscous resistance depends on the dimension and shape of the body. It depends on the surrounding medium viscosity in which the body moves.

Power of the resistive force (power dissipation)

Ps=Fov=rv2E16

is a quadratic function of the velocity. In the case of the electrical circuit, the power dissipation is expressed as P = Ri2. If the electrical current is analogous to the speed of motion, see Example 1, then this equation is analogous to Eq. (16) for viscous losses.

Motion equation ma = Ffor the linear system with viscous damping has a form

ma=kxrv,E17

which can be rearranged to

x¨+2bẋ+ω02x=0,E18

where v=ẋand b = r/(2 m) is damping coefficient.

This equation is a linear differential equation with constant coefficients, and we find the solution in a form of the exponential function eλt. We obtain the values of the λfrom the characteristic equation

λ2+2+ω02=0,E19

which solution is

λ1,2=b±b2ω02.E20

The type of motion of this oscillatory system depends on the ratio of the b, and ω0 values, which defines the quality factor

Q=ω02b.E21

1.2.1 Underdamped oscillation system

Underdamping occurs in systems with a quality factor of Q > 1/2. Characteristic equation solution corresponds to a complex number λ = −b ± , where ω=ω02b2. Then, the solution can be expressed as

xt=Aebtcosωt+α,E22

where Aand αare integration constants and they depend on the initial conditions of the movement, which are the initial particle displacement of x0 and initial velocity of v0 in time t = 0

x0=Acosαandv0=bAcosαωAsinα,E23

from where tanα=1ωv0x0+band A=x02+v0ω21+bx0v02.

See Figure 10 as example, where are underdamped oscillations for different values of the attenuation coefficients b = 0.1 s−1 and b = 0.5 s−1, or for quality factors Q = 5 and Q = 1, respectively, at ω0 = 1 rad∙s−1, and initial conditions x0 > 0 and v0 = 0 m∙s−1.

Figure 10.

Time course of subcritical damped oscillations for two values of the attenuation on the left isω0 = 1.0 rad∙s−1, b = 0.1 s−1, v0 = 0 m∙s−1, on the rightω0 = 1.0 rad∙s−1, b = 0.5 s−1, v0 = 0 m∙s−1.

In the case of the underdamped system, the particle displacement overshoots the zero value (see the negative values in the graphs).

The value of

τ=1b=2Qω0E24

is damping time constantand it indicates the time when the ebtfunction decreases to the value of 1/e ≈ 0.37(=37%). This constant provides information about the time when the oscillations disappear. Usually, we consider the disappearance time of 3τ, when the maximal particle displacement reaches e−3 ≈ 5% of its initial value, or the time of 5τ, at which the displacement drops to e−5 < 1%.

The ratio

τT=ω2πb=1πQ114Q2E25

represents the oscillation count during the time of τ. We can see that there are no oscillations in the system if Q ≤ 1/2.

1.2.2 Critical damping

Critical damping occurs if b = ω0, and Eq. (19) has only one double solution. In this case, the solution of the equation is

xt=A1+A2tebt,E26

where initial conditions are x = x0 and v = v0 at t = 0 s to determine A1 and A2.

Figure 11 shows typical time courses for different initial conditions.

Figure 11.

Time course of critical damped oscillations for ω0 = b = 1 s−1 (on the left x0 = 1 mm, v0 = 0, on the right x0 = 0, v0 = 1 mm∙s−1).

The importance of critical damping is that the system returns from the non-equilibrium state to the equilibrium fast and without overshooting. Various systems utilise critical damping, for example, shock absorbers for vehicles such as cars, motorcycles, etc. Critical damping is also used in the impulse electrical circuits to minimise distortion of the rising and falling edges of the impulse signal.

1.2.3 Overdamped oscillation system

Overdamping is given by b > ω0. If we denote a=b2ω02, then the solution of the Eq. (18) has a form as

xt=ebtA1eat+A2eat,E27

where A1 and A2 result from the initial conditions. The particle displacement over the time consists of two exponential functions while one function has a short relaxation time τ1 = 1/(b + a) and the second function has a time of τ2 = 1/(b − a). Figure 12 shows examples of critical damped systems for different initial conditions.

Figure 12.

Overdamped oscillation system for ω0 = 1.00 rad∙s−1, b = 1.10 s−1 (on the left x0 = 1 mm, v0 = 0 mm∙s−1, on the right x0 = 0 mm, v0 = 1 mm∙s−1).

Dashed lines in the graphs indicate both exponential components with different time constants. We can see that this is an aperiodic event with no overshoot through the equilibrium position.

Viscous damping occurs especially in the case of small oscillations of a mass in the liquid, when there is laminar flow, or in the case of capillary damping devices. Linear damping is also typical for oscillations of atoms due to heat exchange, or for damping of oscillations in electrical circuits. In the case of the mass movement in a gaseous medium, for example, the pendulum in the air, the aerodynamic drag force F ∼ v2 usually applies, which is characterised by a quadratic dependence on speed. This means that it is no longer a linear system, and the solution leads to a non-linear differential equation even at small oscillations.

Example 7. Pendulum in a liquid.

Consider a pendulum (Example 2), whereby the suspended ball moves in water in a dish. For low velocities, the viscous resistance force for the ball-shaped body is given by the Stokes relation

F=6πηrv,

where ηis the dynamic viscosity of the liquid, ris the ball radius and vis the velocity of the motion.

The attenuation coefficient follows from (18):

b=3πηrm=9η4ρr2,

where ρis the ball density.

For example, the water has η = 1.0 × 10−3 Pa·s (at 20°C), the density of steel is 7.8 × 103 kg·m−3 and the ball radius r = 5.0 mm. Then, we get b ≈ 1.2 × 10−2 s−1.

If the ball hangs on the thread of the length l = 1.0 m, then ω0 ≈ 3.1 s−1.

The Q-factor is Q ≈ 130. It is, therefore, subcritical damping and according to (25)

τT=ω2πb=1πQ114Q241.

It means that the oscillations are damped to the ratio of 1/e ≈ 37% after 41 periods.

Example 8. Oscillation damping in electrical RLCcircuit.

Consider a single loop of series-connected elements of an inductor L, a capacitor C, and a resistor R. Assume that initially, the capacitor was charged to a U0 voltage, and the current in the circuit was zero (RLconnection to the charged Ccapacitor).

The energy of the conservative energy components is then

E=12Q2C+12Li2.

The time change of this energy is equal to the power of Joule’s losses

dEdt=QCQ̇+Lii̇=Ri2,

where Q̇=i. If we divide the equation by the current i, and knowing the i̇=Q¨, we get

QC+LQ¨=RQ̇,resp.Q¨+2R2LQ̇+1LCQ=0,

where ω0=1LCand b=R2L.

For example, L = 50 mH, C = 20 μF, and R = 10 Ω, we obtain ω0 ≈ 1.0 × 103 s−1 and b ≈ 1.0 × 102 s−1.

Thus, there are subcritically damped oscillations.

The charge time response is then

Qt=Q0ebtcosωt,

where ω=ω02b2≈ 0.99 × 103 s−1.

τT=ω2πb=1πQ114Q21.6.

We can see that the angular frequency ωdiffers only slightly from the angular frequency ω0 of the non-attenuated oscillations. However, the motion is significantly attenuated. The relative decrease to 1/e ≈ 37% of the initial value occurs after 1.6 periods of oscillations.

1.3 Oscillation of damped system with harmonic excitation

If an external periodic excitation force acts on the oscillation system, the system responds, after the transient process has disappeared, with a periodic answer. If the excitation is harmonic and the system is linear, then the steady answer is also harmonic with the same frequency. Any periodic stimulus of the linear system represents a superposition of harmonic components in terms of the Fourier series. Therefore, we will pay special attention to the response of the linear oscillation system to the external harmonic excitation.

1.3.1 Spectral characteristics of linear system with harmonic excitation

The harmonic force acting on linear oscillation system with viscous damping is given by the equation of motion

Fv=FmsinΩt,E28

where Fmis force amplitude and Ωis its angular frequency.

Then the equation of motion has a form

ma=kxrv+Fv,E29

which we can rewrite to

x¨+2bẋ+ω02x=fmsinΩt,E30

where fm = Fm/mis external force amplitude related to the mass of the system.

The solution of the homogeneous equation corresponds to some of the results of the section 1.2 depending on the type of the system damping. A particular solution respects the right side. The homogeneous solution is a transient that fades out over time. The particular solution represents a process that lasts as long as the exciting force acts. There are steady harmonic oscillations in the system. In the case of harmonic excitation, the particular solution has a form

xpt=xmsinΩt+β,E31

where

xm=fmω02Ω22+22,andβ=arctg2ω02Ω2.E32

The xmis the amplitude of oscillations and βis the phase shift of the response compared to the phase of the excitation force (28). These results can be convinced by directly substituting the solution (31) into the Eq. (30). The linear oscillation system must respond to a harmonic response with the same angular frequency. As can be seen from the previous relationships, the amplitude and phase shift of the response depends on the Ωangular frequency of the excitation.

A special case is the excitation response with an angular frequency which is equal to the angular frequency Ωr = ω0 of the undamped system. This case is called resonance. For the resonance state, we get values from relations (32)

xmr=fm2bω0=x0Q,andβr=π2rad,E33

where x0 = Fm/kis the displacement from the equilibrium while the constant force Fmacts on the system (zero angular frequency Ω = 0). In the case of a system with a high Q-factor of Q≫ 1, the amplitude of the response in the resonance state is significantly greater than the displacement of x0 caused by the constant force. The response is phase-delayed by π/2 rad compared to excitation. Figure 13 shows the frequency response characteristics for different Q-factor values.

Figure 13.

Amplitude and phase shift response of the oscillating system versus a relative angular frequency for different values of the quality factor.

We can see from these characteristics that if the resonant amplification of the system oscillations is undesirable, it is necessary to choose critical or overcritical damping. In this case, however, considerable energy losses occur in the system because of the resistance force. On the other hand, there are systems with low internal losses and characterised by a very high Q-factor (in hundreds to thousands). In the case of high values of the quality factor (Q≫ 1), the frequency bandwidth of the resonant maximum can be determined at a level of 3db decrease relative to the maximum value (i.e., decrease to approximately 0.707 xmr)

ΔΩ=ω0Q.E34

The resonant maximum increases proportionally with the Q-factor and narrows inversely with it. Therefore, the systems with a very high Q-factor have high selectivity, and we can use them, for example, for the spectral analysis of an unknown signal or for controlling pendulum clocks, resonant crystal clocks, atomic clocks, etc. An undesirable consequence of resonance in mechanical devices can be the occurrence of vibrations, for example, when the engine rpm corresponds to the resonant frequency of the mechanical system. These phenomena are not limited to mechanical systems only. Similarly, resonance phenomena occur in electrical circuits or electromagnetic systems. An example, to be mentioned later, is magnetic resonance used in medical diagnostics. In a very simplified view, the human auditory organ is a complicated resonant system too that allows different sound frequencies (pitch of tones) to be distinguished.

1.3.2 Non-linear oscillating system with harmonic excitation

The situation is more complex in the case of a non-linear oscillating system exposed to external harmonic excitation. As an example, consider the non-linear system with the asymmetric potential well described in Section 1.1.2, with harmonic excitation and viscous damping. The equation of motion expressed in a standard form, see Eqs. (12) and (30), has a form

x¨+2bẋ+ω021+λxx=fmsinΩt.E35

We are interested again in the steady-state response of the system described by the particular solution of the differential equation. The response of a non-linear system to harmonic excitation is no longer harmonic but remains periodic with the same angular frequency. The periodic response function is expressed as a superposition of harmonic components using the Fourier series

xt=x0+n=1xmnsinnΩt+βn.E36

After substituting this assumed solution into the differential equation, we obtain the values of the individual quantities. In the case of weak non-linearity (λxm1 << 1), the results have the form of

xm1=fmω02Ω2+2λω02B02+2bΩ2fmω02Ω22+2bΩ2,E37
x0=12λxm12E38
xm2λω022xm12ω024Ω22+4bΩ2E39
xm3λω022Bm1ω029Ω22+6bΩ2xm2,etc.E40

The fundamental harmonic having amplitude xm1 dominates. Its properties are similar to those of the linear system. Also, the mid-position x0 is shifted due to the system’s non-linearity. Furthermore, the resonance occurs at subharmonic frequencies, an integer fraction of the fundamental harmonic frequency (Ωn = ω0/n). They are expressed by response amplitudes xmn. The subharmonics components have an origin caused by excitation having a specific subharmonic frequency Ωn. But the response has the fundamental resonance frequency ω0 since there is the response of specific harmonic defined as n = ω0. Figure 14 shows the frequency amplitude characteristics of the first and second harmonics.

Figure 14.

Frequency amplitude characteristics of the first and second harmonic components for values ofλ = 0.06 m−1, Q = 50.

The subharmonic resonance is important to explain the perception of musical chords by the non-linear system of the auditory organ. For example, if we hear two tones with frequencies in the ratio 2:1 (octave), the tone with the angular frequency ω01 produces a signal with the second harmonic of 2ω01. If this frequency is not equal to the frequency of ω02, the auditory organ sensitively detects the difference between ω02 and 2ω01 and evokes a feeling of non-tuned music interval.

1.3.3 Harmonic interaction in a non-linear oscillating system

In practice, we encounter cases, in which the oscillating system is simultaneously excited by several harmonic signals with different frequencies.

As a simple example, we will excite the system with two harmonic signals and determine its response to this excitation.

The equation of the response has the form

mx¨+rẋ+kx+l/2x2=Fm1sinΩ1t+Fm2sinΩ2t,E41

where the Ω1 and Ω2 are angular frequencies of the harmonic components of the excitation.

We can rewrite the equation to

x¨+2bẋ+ω0ω0+λxx=fm1sinΩ1t+fm2sinΩ2t,E42

where fm1 = Fm1/mand fm2 = Fm2/m.

Since the excitation signal is periodic, the response must also be periodic.

Considering the weak non-linearity when λω0xm, harmonic components with excitation angular frequencies dominate in response. Therefore, the steady response has the dominant components

x1=xm1sinΩ1t+β1+xm2sinΩ2t+β2.E43

If we substitute this function into a quadratic term in the Eq. (42), there are elements with combinational frequencies Ω1 ± Ω2 on the left side of the equation

x12=12xm121cos2Ω1t+2β1+12xm221cos2Ω2t+2β2++xm1xm2cosΩ1Ω2t+β1β2cosΩ1+Ω2t+β1+β2.

Due to non-linearity, components with frequencies 2 Ω1, 2 Ω2, Ω1 + Ω2, and Ω1 − Ω2 appear in the equation. Including these components together with the original components in the overall system response, the non-linearity (quadratic term) results in the second generation of components with twice the frequencies and with all combinations, for example, 3Ω1, 3Ω2, 2 Ω1 ± Ω2, Ω1 ± 2 Ω2. The solution is very complex, and therefore, we will focus on the approximate determination of combination components of the first generation.

In the time response of the system, we consider only the most significant components

xt=x0+xm1,0sinΩ1t+β1,0+xm0,1sinΩ2t+β0,1++xm2,0sin2Ω1t+β2,0+xm0,2sin2Ω2t+β0,2++xm1,1sinΩ1+Ω2t+β1,1+xm1,1sinΩ1Ω2t+β1,1,

where numbered indices correspond to frequency combinations, for example, Qmk,lrelates to the frequency of Ωk,l = 1 + 2.

If we substitute these components into the equation of motion and separate the corresponding harmonic elements on the left and right sides, we get a response for the amplitudes of the harmonic components

xm1,0=fm11ω02Ω12+2λω0Q02+2bΩ12xm0,1=fm21ω02Ω22+2λω0Q02+2bΩ22,

which corresponds to the frequency response of the linear system (resonant characteristics with a maximum at the resonant frequency of ω0).

Considering the dominant components with angular frequencies Ω1 and Ω2, the constant component is

x012λxm1,02+xm0,12ω0+λQ0λ2ω0xm1,02+xm0,12.

Components with double frequencies of 2 Ω1 and 2 Ω2 shall be determined as in the case of simple harmonic excitation and with the same results.

For the lowest combination frequencies, we get a relationship

xm1,±1λω0xm1,0xm0,1ω02Ω1±Ω22+2λω0Q02+2bΩ1±Ω22.

From this relationship, we can see that the cause of combination frequencies is non-linearity, which occurs in the result as λω0. Similarly, we can determine the amplitudes of the response components with higher combinational frequencies. As with resonance at subharmonicfrequencies, resonance occurs when combinational frequencies are

Ω1±Ω2=ω02+2λω0Q0ω0.

In this case, the response amplitude with an angular frequency of ω0 is given by equation

xm1,±1ω02bλω0xm1,0xm0,1.

System resonances also occur at higher combinational frequencies. Non-linearity and the resulting response components with combinational frequencies increase at higher excitation. There are systems where the combinational frequencies are undesirable. For example, the acoustic loudspeakers are load overrated, which means that the effects of system non-linearity under operating loads are negligible.

As indicated in Section 1.3.2, resonances with combinational frequencies or resonances at subharmonic frequencies are important, for example, in explaining the perception of musical chords by the non-linear system of the human auditory organ. For example, if we hear two tones with frequencies in the ratio of 2 (octave 2:1), the tone with the angular frequency of ω01 produces a signal with the second harmonic of 2ω01. If this frequency is not equal to the frequency of ω02, the auditory organ sensitively detects the difference of ω02 − 2ω01 and creates a sense of non-tuned music interval. It is like other music intervals such as small third 6:5, big third 5:4, fourth 4:3, fifth 3:2, small sixth 8:5, big sixth 5:3, small seventh 16:9, and big seventh 15:8. These music intervals have the ratio of frequencies equal to the integer ratio. Due to the non-linearity of the auditory organ, the music listener can distinguish a pure (harmonic) or impure (disharmonic) chord and thus perceive the beauty of musical compositions.

1.3.4 Power losses and the nature of spectroscopy

The oscillations relate to the exchange of energy between conservative elements of the system. The total energy of the system is equal to the sum of kinetic and potential energy or their equivalents. In the case of forced harmonic oscillations of a linear system with the fundamental frequency of ω0 and excitation frequency of Ω, the total energy of the system at a steady state is

E=12mv2+12kx2=12mvm2cos2Ωt+β+12kxm2sin2Ωt+β==14mΩ2+ω02xm2+14mΩ2ω02cos2Ωt+2β.

The mean value of this energy is

E=14mΩ2+ω02xm2=141+Ω2ω02kxm2.E44

The second (alternating) component is directly proportional to the difference Ω2 − ω02 and corresponds to the periodic energy exchange between the source and the system with an angular frequency of 2 Ω.

If we consider the viscous losses in the system, the energy losses in one period of Tare as follows

WT=0TPdt=0TFvdt=0TFmvmsinΩtcosΩt+βdt==T2Fmvmsinβ=πFmxmsinβ.

In the state of resonance at frequency Ω = ω0, β = −π/2 rad and the alternating component of energy Eis zero.

Then, active power supplied to the system in the case of the steady state of forced oscillations is

P=WTT=12FmxmΩsinβ=Fm22m2bΩ2ω02Ω22+2bΩ2E45

Figure 15 shows the graph of the active power spectral function [see Eq. (45)] for two Q-values. These figures show that an oscillating system with a high Q-factor absorbs the energy of the source only in a narrow interval around the resonant frequency and changes it most often into heat.

Figure 15.

Relative spectrum of the active power of oscillation system (on the right isQ = 10, on the leftQ = 100).

There are various bonds of atoms and molecules in biological tissues. These bonds represent microscopic oscillating systems with characteristic resonant frequencies. If these tissues are irradiated with monochromatic electromagnetic waves with a frequency corresponding to a resonant frequency of coupling, then energy is supplied to these coupled systems. This energy can stimulate tissues at low-power applications, for example, phototherapy. At higher power, the absorbed energy increases the temperature of the tissue structures, and thus, it can lead to their destruction used, for example, in the treatment of cancer.

Example 9. Spectroscopy.

The above phenomenon explains the physical nature of spectroscopy. In systems with a higher Q-factor, the resonance state relates either to dynamically increased oscillations or to power absorption of the source.

The conservative forces bond the atoms of the matter and determine their equilibrium position. The oscillations around the equilibrium position are at the natural frequency and depend on the properties of the particle (mass) and the features of the bond (stiffness). Thus, differently bound particles have different oscillation frequencies. Each material has characteristic frequencies according to its composition.

When an electromagnetic wave interacts with a material, it acts on its atoms. If the EM wave frequency equals one of the resonant frequencies of the substance, then it significantly absorbs and attenuates this wave. For example, if we observe the white light of the Sun on the surface of the Earth using a spectrometer, we find in the continuous visible light spectrum black lines that correspond to the absorption of light with the appropriate frequencies of gas molecules in the atmosphere. In this way, we can measure the concentration of greenhouse gases in the atmosphere.

Another example is optical spectroscopy used in biochemistry, pathology, or the investigation of blood plasma. As an example, let us pass the adjustable wavelength light through the liquid cuvette to search for wavelengths at which the liquid has a resonant absorption. Then, the found wavelengths or frequencies determine the presence of the individual substances of the material and their concentration in the solution.

Another example is magnetic resonance imaging, as discussed in the following paragraph.

1.3.5 Magnetic resonance

1.3.5.1 Nature of magnetic resonance

We talk about resonance if the frequency of external excitation on the oscillating system is the same as the frequency of its self-oscillations, and the mechanism of action can supply the oscillating system with energy. In the linear system, it is the frequency of its undamped oscillations. In the case of a magnetic dipole in a constant magnetic field, it is the Larmor frequencyof fL, see Example 4. Read also Vlaardingerbroek [2], Webb [3], or Hashemi [4].

If we create a rotating magnetic field in space with a rotation frequency fclose to the frequency fL, we can expect a resonance phenomenon. The external excitation magnetic field must be perpendicular to the precession axis (i.e., to the constant magnetic field B0) to interact with a magnetic dipole that performs a precession movement. We create a rotating magnetic field using two mutually perpendicular pairs of coils, which are fed by currents with the same frequency and with a mutual phase shift of π/2 rad.

Figure 16 illustrates the situation where perpendicular pairs of coils are on the left. If zis the direction of the constant magnetic field B0 and hence the axis of the dipole precession, the xand ydirections are perpendicular to the z-axis. One pair of coils creates the magnetic field of Bxin the x-axis direction, the other pair the field of Byin the y-axis direction. The coil currents and thus the magnetic field components are phase-shifted by π/2 rad, and thus

Figure 16.

Magnetic dipole in rotating magnetic field.

Bx=B1sinωt+ψ,andBy=B1cosωt+ψ.

Adding both components, we get the resulting B1 vector, which has a constant value of B1 and rotates in the x-yplane with the ωangular frequency of the coil current. The right part of the figure shows the direction of the dipole precession with the dipole moment of mand the direction of rotation of the rotating magnetic field.

Vector components of the mdipole moment are

mz=mcosα,mx=msinαsinωLt+ϕ,andmy=msinαcosωLt+ϕ.

Magnetic dipole in the magnetic field B0 has potential energy

Ep=mB0=mB0cosα.E46

If an external alternating magnetic field acts on the dipole, only the angle αcan change at constant values of mand B0. The following equation expresses the change of the potential energy dEpof the dipole

dEp=mB0sinαdα.E47

Power of external magnetic field torque is

P=MωL=m×B1ωL==mB1sinαsinωLt+ϕcosωt+ψcosωLt+ϕsinωt+ψωL==mB1ωLsinαsinωLωt+ϕψ.E48

We can see that the power is time-varying for ωL ≠ ωand the mean value of the power is zero. If ωL = ω, the power has a time-invariant component, which reaches the maximum at φ − ψ = π/2 rad. We describe this phenomenon as magnetic resonance.

Then the work over the time dtis

δW=Pdt=mB1ωLsinαdt.E49

Comparing (47) and (49), we obtain for dEp = δW

dαdt=B1B0ωL=γB1.E50

As a result, the αangle of the ‘precession funnel’ varies uniformly in the magnetic resonance state with an angular velocity of dα/dt, which depends on the amplitude of the induction of B1 of the alternating magnetic field. These angle αchanges are the periodic event, and therefore the magnetization of a substance changes periodically too. The magnetization inverts its value in the time

τ180=πγB1,E51

or it is perpendicular to the initial direction in the time

τ90=π2γB1,and similar.E52

It is typical for a forced oscillation of particles, and a forced precession of magnetic dipoles, that all particles oscillate synchronously with the same phase compared to the excitation signal.

1.3.5.2 FID signal origin

The paramagnetic material contains many magnetic dipoles randomly arranged due to particle thermal motion. Therefore, the resulting magnetic field of these dipoles is zero. If we insert the paramagnetic material into the B0 constant magnetic field, then the material magnetic dipoles partially arrange in the direction of the B0 vector. This behaviour better describes the magnetization vector (M0 = κμ0B0), where the κis the magnetic susceptibility of the substance. After switching on the B1 transverse alternating rotating magnetic field with an angular frequency ω = ωL, a resonance occurs, which causes a coherent precession of the oriented magnetic dipoles. Let us apply the field B1 during the τ90 time. Then the M0 constant magnetization vector, parallel to the B0 vector, changes to the Mvector, which has the same magnitude but rotates perpendicularly to the B0 with the angular frequency ω. The sample of a substance looks like a rotating magnet with a magnetic moment (m* = M0V), where Vis the sample volume. If we place a detection coil perpendicularly to the axis of rotation, then the voltage induces in it is

uFIDdmdtωLM0sinωLt.E53

Voltage induces in the coil only in the case of synchronous dipole precession, which results in rotating magnetization, and this can only happen if magnetic resonance conditions are met. For a given B0 and ω, the resonance occurs only for certain dipoles in the substance, which satisfy the condition ω = ωL = γB0. Thus, by measuring the induced signal, the presence of magnetic dipoles with a corresponding gyromagnetic moment γcan be detected, and their concentration determined.

If we switch off the B1 excitation field, the periodic event begins to damp due to the interaction of the dipoles with the surrounding particles of the substance. The detected signal is, therefore, attenuated (Figure 17). This damped signal calls the FID signal (free induction decay). In biomedicine, the magnetic resonance uses protons (with the γ = 2.68 × 108 s−1·T−1), which exist mainly as nuclei of hydrogen atoms and thus in water molecules.

Figure 17.

FID signal after magnetic dipoles excitation.

The organic compounds, such as biological tissues, contain hydrogen atoms too. In specific cases, the magnetic resonance uses nuclei of other biogenic elements such as isotopes of carbon 13C, fluorine 19F, phosphorus 31P, and so on (see Table 1, p. 8).

1.3.5.3 Relaxation

Perpendicular magnetization is an imbalance caused by the external source of the alternating magnetic field B1. If the excitation force stops to act on the system, the aligned movement of the dipole array decays. From the viewpoint of the FID signal, the decay of the in-phase periodic precession movement at first occurs due to the inhomogeneity of the magnetic field B0, and due to the influence of surrounding dipoles, so-called spin-spin interaction. After switching-off the exciting magnetic field, the precision movement of the dipoles remains for a short time in a plane perpendicular to the B0 direction, but due to a small change in the local magnetic field, the precession of the single dipoles is out-phased, which results in an exponential decrease of the transverse magnetization, and thus an FID signal. This decrease characterises the time constant T2. Its value is in the order of tenths of a second. The second slower mechanism of decay associates with the thermal relaxation of the imbalanced orientation of the magnetic dipoles and directs to the thermodynamic equilibrium of the dipoles, that is, to the equilibrium orientation of magnetization in the direction of the B0 vector. This process is approximately 10 times slower, and its time constant is denoted T1.

Different substances, and thus tissues, have different values of relaxation times of T1 and T2. In medical applications, protons (nuclei of hydrogen) are mostly used as magnetic dipoles since the body contains many of the hydrogen atoms (especially as part of water molecules).

For illustration, see Table 2, which contains values of relaxation times T1 and T2 for water and some tissues, as well as the relative concentration of hydrogen atoms in tissues compared to the concentration in pure water.

TissueT1 [ms]T2 [ms]Relative concentration (1H)
Water400020001.00
Cerebrospinal fluid25002800.98
Edema9001300.86
Grey matter760770.74
White matter510670.62
Muscle900500.50
Fat250601.00

Table 2.

Relaxation times and relative concentration of protons in water and selected tissues.

Chemical analyses also use nuclei of other paramagnetic atoms as magnetic dipoles. Thus, we can investigate the content of specific atoms or substances in the samples by measuring the FID signal and the relaxation times.

1.3.5.4 Magnetic resonance imaging (MRI)

One of the applications of the magnetic resonance phenomenon is the tomographic imagingof the morphological structure of the organism. The method lies in the use of the detection of hydrogen atom nuclei, which are mainly contained in water and thus in soft tissues. Consequently, we can obtain a two-dimensional image of tissue structures by identifying different types of tissue (see Figure 18).

Figure 18.

MRI of the cervical spine, part of the vascular system, thorax.

Using a relatively complicated device we call a tomography; it is possible to assign a specific T1 and T2 value or relative proton density PD to each point of the thin transverse layer of the examined object (body) and thus to distinguish individual tissues. Different values of these quantities are assigned a certain level of grey colour when displayed on the device monitor (see Figure 18). In this way, we can obtain different images such as T1-image, T2-image, and PD image. Each of them has a different contrast concerning tissue differentiation and, thus, different advantages in medical diagnostics.

1.3.5.5 Magnetic resonance spectroscopy

The second application of magnetic resonance is magnetic resonance spectroscopy(MRS). By variation frequency ω, it is possible to select the type of atomic nucleus with the Larmor frequency of ωL. Then, we identify the nucleus by the magnetic resonance FID signal at the frequency of ω = ωL. The magnetic fieldB0 at the location of a given nucleus, and thus the Larmor frequency of ωL, is slightly influenced by the magnetic field of the surrounding particles, such as electrons and other nuclei. The resonance frequency ωLof the atom nucleus is thus slightly influenced by the chemical bonds where the magnetic dipoles (e.g., nuclei of hydrogen atoms) occur. By examining the spectrum of resonances, it is possible to identify individual hydrogen bonds in the sample under investigation, for example, O-H, C-H, C-H2, C-H3, N-H2. Furthermore, we can identify the relevant organic substances (protein, enzyme, and metabolite) according to the measured resonance spectrum.

Figure 19 illustrates the organosilane spectrogram used in the manufacturing process of synthetic rubber. The horizontal axis is the offset of the resonant frequency in parts per million (ppm = 10−6) relative to the reference frequency. The reference substance could be tetramethylsilane (TMS) or another proper substance. For example, hydrogen in the = CH2 divalent group has a resonance frequency shifted by 1.3 ppm (A), in the -CH3 monovalent group, up to 4.0 ppm (B). Each substance has a characteristic spectrogram according to which we can identify it, even at a very low concentration.

Figure 19.

Magnetic resonance spectrogram of organosilane.

Magnetic resonance spectroscopy thus enables very sensitive biochemical diagnostics of different tissues or fluids and uses various biochemical markers to early diagnose a variety of diseases, such as epilepsy, Alzheimer’s disease, Parkinson’s disease, various cancers. Thus, magnetic resonance spectroscopy is a powerful diagnostic tool in medicine.

In specific cases, instead of hydrogen, the magnetic resonance spectroscopy uses the nuclei of other biogenic elements with an uncompensated magnetic moment such as 13C, 19F, 31P. The MRS apparatus is quite demanding, and therefore, a special investigation of the content of other nuclei is used only rarely. Thus, the MRS uses preferably only 1H (hydrogen-protons) for the determination of metabolite content, which in addition to MRI does not require additional devices and MRI and MRS images can be combined (see Figure 20). On the right side, it is an MRI image with defined the specific location of analysis, on the left side is an MRS spectrogram of the substance at that location. From the spectral peaks typical for certain substances (here Cr-chromatin, Cho-choline, NAA-N-acetyl aspartate) and their size, it is possible to diagnose possible health disorders.

Figure 20.

Combination of MRS and MRI.

1.3.5.6 Magnetic resonance therapy

Magnetic resonance therapy (MRT) is a treatment method that uses targeted stimulation of specific structures by providing them with energy through magnetic resonance. During a resonant RF excitation pulse of the τ180 length, the alternating magnetic field supplies the dipole with energy [see Eq. (49)]. This energy is transferred only to the nucleus of the atoms that are in resonance with an alternating magnetic field. We can supply the energy of the electromagnetic field to specific parts of the structure that contain the resonant nucleus of the atoms. Thus, we can stimulate intracellular processes such as cell nucleus growth.

The method of magnetic resonance therapy is successfully used in the treatment of osteoarthritis and osteoporosis as we supply the energy to help cartilage and bone regeneration, as well as recovery for spinal pain following surgery (see Figure 21).

Figure 21.

MRT—Hip joints on the left, post-hip treatment on the right.

1.4 Oscillators

We are using various sources of periodic signals or motions, which are commonly called oscillators. Oscillators, mechanical or electrical, are systems with high Q-factor value and low losses having a frequency f0 determined by the system parameters. However, each system always has, albeit small, losses that cause the oscillation to disappear at a time proportional to the quality factor [see Eq. (24)]. If the system is to oscillate continuously, we must balance its losses. This compensation consists of supplying energy equal to the losses in each period of oscillation, that is, the compensation process must be synchronous with the system’s oscillations. We can achieve this by periodic power supply directly controlled by system oscillations, which means a positive feedback method. The classic example shows a child on a swing. If a child sits on a swing and the parent pushes it, it will swing for a while, but it will soon hang in a steady position. Children almost intuitively understand to keep the swing in motion. They must compensate for the loss of energy by properly digging their legs in one extreme position and kicking in the other, utilising the energy of their muscle activity to increase the potential energy twice within one oscillation slightly. The child performs this activity intuitively. Thus, the child’s biological energy compensates for the energy losses of the swing.

Oscillators have a precisely defined frequency by their parameters. Therefore, we can use them as a reference time signal source. Thus, they represent the essential part of the clock (mechanical with pendulum, mechanical with the rotating flywheel on spiral spring, electrical with LC circuit, electrically controlled with crystal, atomically controlled with quantum transitions in caesium atoms). The electronic clock is a part of every computer and controls the operation of such components as the processor, data storage, and data exchange with peripherals.

1.4.1 Mechanical oscillator

A commonly known mechanical oscillator is a pendulum clock. Figure 22 shows a pendulum (dashed line) and a zoomed positive feedback step mechanism. The step wheel with inclined teeth is driven through the gearing by a force Fgenerated by a weight or a spring. In the picture, the pendulum moves to the right and the right inclined tooth ‘b’ pushes into the stop of the escapement and supports the right-hand rotation. After reaching the extreme position, the tooth is released to the right, and the wheel rotates so that the left step-stop rests on the left oblique tooth, which pushes into the stop of the escapement and supports the pendulum moving to the left. Thus, the inclined teeth of the wheel supply energy to the pendulum via a step mechanism. The wheel drive depends on the potential energy source of the weight or spring. The system is set up to maintain a stable pendulum operation.

Figure 22.

Pendulum.

There are many mechanical oscillators of analogous construction, for example, a flywheel on a spring in a mechanical wristwatch, a torsional pendulum of a decorative stand clock. The pendulum clock accuracy depends on the temperature regarding the thermal expansion of the mechanical parts. A special temperature-stabilised pendulum clock can achieve running stability of up to δT/T0 ≈ 10−6 (1-second deviation in 12 days).

1.4.2 LC oscillator

The electric oscillators commonly use the LC circuit with the frequency of natural oscillations of f0=12πLC. Due to the electrical resistance of the circuit, energy losses occur, which leads to oscillation damping. To cover energy losses and maintain the oscillations of the system, we must supply the LC circuit using a positive feedback method in connection with an amplifier. There are many LC circuit oscillators; Figure 23 shows some examples.

Figure 23.

Different types of electrical oscillators.

These oscillators use a transistor amplifier connected with a common emitter that changes the signal phase by 180°. We must connect the output voltage of the oscillator to the input with the same phase, respectively, with offset by 2 × 180° = 360°. The input part is an LC oscillating circuit with a split capacitor: (a) Colpitts circuit, or a split inductor and (b) Hartley circuit. As shown in the figures, there is an opposite phase on split elements regarding the amplifier input and output. The (c) case shows the Meissner circuit, where phase reversal is achieved by inductive coupling with oppositely oriented windings. Figure (d) shows an example of an RC oscillator that does not use an LC circuit.

We achieve positive feedback by a three-stage RC phase shifter. The elements have a total phase shift of 3 × 60° = 180° at the desired oscillation frequency. Since the phase shifter is frequency-dependent, positive feedback occurs at only one frequency.

1.4.3 Crystal controlled oscillators

The applications demanding higher frequency stability consider the circuits mentioned above as unsatisfactory due to the used circuit elements. For example, the inductors are highly temperature-dependent, or parasitic elements influence both the oscillator and amplifier circuits, and then voltage fluctuations, are applied. To suppress these parasitic effects, we are using piezoelectric crystals in the oscillating circuits instead of the inductors.

A piezoelectric crystal is an electromechanical oscillating system with a high Q-factor. This crystal is described by using the equivalent circuit diagram, as shown in Figure 24(b). The inductance of the crystal depends on the mass of the crystal. The capacity corresponds to its rigidity and the resistance to the internal power losses. The capacity C0 is the electrode capacity of C0C. For the crystal as a reactance electrical circuit, the imaginary part of the complex impedance is important, and we can express it as follows:

Figure 24.

Crystal controlled oscillator.

X=ωL11ω2LC1+CC0ω2LC0R2C01+C0Cω2LC02+ωRC02.E54

Figure 25.

Graph of the reactance X versus frequency f ot the crystal.

If we set the reactance equal to zero (X = 0), we can estimate the resonant frequency of the crystal. The reactance graph below (Figure 25) shows two resonant frequencies for the given values of the crystal (L = 100 μH, C = 100 pF, R = 1.0 Ω, C0 = 10 nF). There is the fs parallel and the fp series resonances. The equivalent circuit with a high Q-factor has the resonant frequencies as follows:

fs=12πLC,andfp=12πLCC0/C+C0.

In our case, fs ≈ 1.5916 MHz and fp ≈ 1.5994 MHz. As shown in the reactance graph, we can see a narrow interval between the fsand fpwhen X > 0. This means that the crystal has an inductive character. The crystal is connected in the oscillating circuit as an inductor with a parallel split capacitor with C1 and C2 capacitances. Thus,

Figure 24(a) shows Pierce’s circuit. The split capacitor is parallel connected to the C0, and therefore, the interval between resonant frequencies gets narrower. Thus, the oscillator can oscillate in the very narrow frequency range, which ensures high stability of the oscillator frequency. We realise positive feedback by connecting non-inverting output through the R2 resistor to the C2 capacitor. On the other hand, we can tune the oscillator in the range of several Hz. If we need to tune the frequency in the broader range, we must change the crystal. The crystal-controlled oscillator has high stability in order of 10−9, which means the time deviation of 1 s for 30 years. Achieving this stability, we use the thermostat to stabilise the temperature of the crystal. Due to the high-stability requirement, the computer clock uses only crystal-controlled oscillators as the clock pulse generator.

1.4.4 Multivibrators

In some cases, we require a harmonic signal for biomedical applications. There are diathermy, electrotherapy, sonography, or magnetic resonance. In other applications, we need to generate periodic, but non-harmonic voltages or currents. There are pacemakers or artificial lung ventilation. In these cases, we are using rectangular or sawtooth waveforms or short repetitive pacing pulses. The primary element of non-harmonic signal generators are multivibrators serving as sources of periodic rectangular pulses. Many mechanical and thermal devices switch between two states at regular intervals. There are electrical systems, which are the most important for biomedical applications. These systems serve as periodic and non-harmonic voltage sources. As an example, Figure 26 shows the circuit of an astable flip-flop multivibrator.

Figure 26.

Astable flip-flop circuit.

In the principle, the T1 transistor alternately switches between its ON/OFF states. If the T1 is open, the T2 is closed and vice versa. This process repeats periodically. The toggling period is given by time constants defined as R2C1 and R3C2. The output voltage is then rectangular. Connecting output to the differentiator circuit, we obtain short pulses, which can be used for the pacemakers. Using the integrator, we get a sawtooth waveform, which can be used for the generation of the linearly rising gradient field at magnetic resonance imaging.

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Ivo Čáp, Klára Čápová, Milan Smetana and Štefan Borik (December 24th 2021). Simple Oscillating Systems, Electromagnetic and Acoustic Waves in Bioengineering Applications, Ivo Čáp, Klára Čápová, Milan Smetana and Štefan Borik, IntechOpen, DOI: 10.5772/intechopen.101649. Available from:

chapter statistics

37total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Electromagnetic and Acoustic Waves in Bioengineering Applications

Authored by Ivo Čáp

Next chapter

Coupled Oscillators

By Ivo Čáp, Klára Čápová, Milan Smetana and Štefan Borik

Related Book

First chapter

Applications of Heat Transfer Enhancement Techniques: A State-of-the-Art Review

By Suvanjan Bhattacharyya, Devendra K. Vishwakarma, Sanghati Roy, Ranjib Biswas and Mohammad Moghimi Ardekani

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us