Delay Harmonic Oscillator

The motion of two point objects at the end of a spring is analyzed. The objects interact by an elastic wave propagating through the spring. A new comprehensive method, Reaction Mechanics, for the analysis of this motion is used. This analysis is valid when the propagation of the interaction through the spring wire takes less time than the period of the oscillating frequency. The propagation delay couples the oscillating and center of mass motions. If the masses are equal, the center of mass velocity is a constant, and the objects oscillate with a frequency which is a modification of the oscillation frequency with no delay. If the masses are not equal, the center of mass also oscillates. In the case of zero delay, the motion of the objects reverts to the motion of a Simple Harmonic Oscillator.

• Lagrangian
• Equation of Motion
• Hamiltonian
• Delay
• Oscillation

1. Introduction

The motion of two objects connected by a spring is analyzed. The interaction propagates between the objects as a mechanical wave through the spring as shown in Figure 1. The mechanical wave propagates through the spring with a velocity [1] that depends on the stiffness of the spring and the mass of the spring per unit length. Therefore, there is a time delay τ required for the deformation applied to one end of the spring to reach the other end. The interaction is not instantaneous. This is seldom taken into consideration when systems of springs and masses are calculated.

In this paper, a new comprehensive method called Reaction Mechanics [2] is used for the analysis of the motion of a system that takes the propagation delay of the interaction into account.

The interaction travels between the two objects, through the spring, with a velocity v_S equal to the square root of the ratio of the spring stiffness constant in Newtons, divided by the mass density of the spring in kg per m. The spring stiffness is equal to an elastic constant times a cross sectional area. The analysis is valid when the propagation of the interaction through the spring wire takes less time than the period of the oscillating frequency. The propagation delay through the spring couples the Oscillatory and Center of Mass motions of the objects.

It is demonstrated that the delayed interaction [2] is completely causal. A force F_R that was radiated by object TWO in the past is received currently by object ONE. Object ONE also radiates a force F_A currently that object TWO might or might not receive in the future.

By Newton’s law of action and reaction object ONE experiences a recoil force F_Recoil equal in magnitude and opposed in direction to the force F_A it radiated. The recoil force is independent whether or not object TWO receives the radiated force in the future. The forces are continually radiated and received.

Two forces act on object ONE. The force F_R radiated by object TWO in the past and the recoil force F_Recoil. These two forces are equal to the mass of object ONE times its acceleration. Thus, the formulation is causal [2] .

The calculations show that the two masses and the spring oscillate. But, if the two masses are not equal, the Center of Mass of the objects will also oscillate.

The results of this paper will be useful for the analysis of machines with coil springs and masses. The most obvious application is the vehicle suspension [3] . The vehicle suspension has two masses connected by coil springs and vibration dampers. However, the vibration dampers are not part of this analysis. The upper mass is the body, chassis and engine, and the lower mass is the axels and wheels. The upper mass is much heavier than the lower mass. Therefore, the center of mass will vibrate. If the vibration dampers were not there, it would cause a large vibration of the upper mass.

If one uses the standard calculation, neglecting the propagation delay through the spring, one would then obtain that the two masses vibrate about the stationary center of mass.

This type of motion also occurs when an electric motor is mounted on a platform that is suspended on elastic supports. Even if the rotation of the motor is well balanced, the platform might still vibrate because there are different masses on both sides of the elastic support.

The intake and exhaust valves with their springs of an internal combustion engine can cause oscillations in the engine block.

2. Equation of Motion

In this analysis, the masses are treated as point objects and the model of the spring includes the propagation delay. The propagation delay is modeled by considering that the forces, implicit in Equation (1), with components K x μ k and K y μ k − 2 at the two ends of the spring, do not occur simultaneously. It is assumed that the time difference t k − t k − 2 between the forces occurring is equal to the propagation delay through the spring. This eliminates the necessity of describing the spring by a Continuum Mechanics model. It is assumed that the mechanical signal propagates through the spring, from one object to the other, in a shorter time than the period of oscillation of the masses and the spring. In this case, the causal Lagrangian [4] for the Delay Harmonic Oscillator (DHO) can be approximated as follows:

L k = 1 2 m x ˙ μ k x ˙ μ k + 1 2 M y ˙ μ k y ˙ μ k − 1 4 K ( x μ k − y μ k − 2 ) ( x μ k − y μ k − 2 )     − 1 4 K ( y μ k − x μ k − 1 ) ( y μ k − x μ k − 1 ) (1)

where summation over repeated Greek indices is implied. The subscripts µ = 1, 2, 3, range over the tree dimensions of the coordinates. Single superior dots x ˙ μ k denote time derivatives and double superior dots x ¨ μ k denote second time derivatives. Here x μ k at time step t_k is a coordinate component, and m is the mass of Object ONE, and y μ k at time step t_k is a coordinate component, and M is the mass of Object TWO. K is a spring constant in Newtons per m.

Forming a sum S of Lagrangians.

S = ∑ k = − ∞ k = ∞ L k (2)

Only the Lagrangians L_k and L_k+1 contain the coordinate components x μ k and y μ k occurring at time step t_k. The Lagrangian L_k is advanced to the next time step L_k+1 by advancing the coordinates to their next nodal points, x μ k advances to x μ k + 1 , x μ k − 1 advances to x μ k , y μ k advances to y μ k + 2 and y μ k − 2 advances to y μ k .

As shown in Equation (1) and Equation (3), t k − 1 is a time step of object ONE in the past and t k − 2 is a time step of object TWO in the past. The time t k − 2 in the past, at which object TWO radiated an interaction that object ONE is receiving currently at time t_k, is not the same as the time t k − 1 in the past at which Object ONE was radiating an interaction that object TWO is receiving currently. A similar condition exists for the future coordinate times. The future forces radiated by the objects do not act on the objects. Recoil forces due to the objects radiating interactions act on the objects. The Lagrangian at time step t k + 1 is:

L k + 1 = 1 2 m x ˙ μ k + 1 x ˙ μ k + 1 + 1 2 M y ˙ μ k + 2 y ˙ μ k + 2 − 1 4 K ( x μ k + 1 − y μ k ) ( x μ k + 1 − y μ k )     − 1 4 K ( y μ k + 2 − x μ k ) ( y μ k + 2 − x μ k ) (3)

Making a transformation of variables [4] .

a ) r μ k = x μ k − y μ k b ) R μ k = m m + M x μ k + M m + M y μ k (4)

where r μ k and R μ k are also in rectangular coordinates. Here r μ k is a oscillation coordinate component and R μ k is a center of mass coordinate component. Inverting Equation (4).

a ) x μ k = M m + M r μ k + R μ k b ) y μ k = − m m + M r μ k + R μ k (5)

Substituting the transformation of variables, Equation (5) into the Lagrangian of Equation (1).

L k = 1 2 m R r ˙ μ k r ˙ μ k + 1 2 ( m + M ) R ˙ μ k R ˙ μ k     − 1 4 K ( M m + M r μ k + R μ k + m m + M r μ k − 1 − R μ k − 1 )     × ( M m + M r μ k + R μ k + m m + M r μ k − 1 − R μ k − 1 )     − 1 4 K ( − m m + M r μ k + R μ k − M m + M r μ k − 1 − R μ k − 1 )     × ( − m m + M r μ k + R μ k − M m + M r μ k − 1 − R μ k − 1 ) (6)

where m R is the reduced mass.

m R ≡ m M m + M (7)

As described previously, the coordinate components x μ k − 1 and y μ k − 2 do not occur simultaneously. However, a simplification is made in Equation (6) and Equation (8) that assumes that those coordinate components do occur simultaneously. Substituting the transformation of variables, Equation (5) into the Lagrangian of Equation (3)

L k + 1 = 1 2 m R r ˙ μ k + 1 r ˙ μ k + 1 + 1 2 ( m + M ) R ˙ μ k + 1 R ˙ μ k + 1     − 1 4 K ( M m + M r μ k + 1 + R μ k + 1 + m m + M r μ k − R μ k )     × ( M m + M r μ k + 1 + R μ k + 1 + m m + M r μ k − R μ k )     − 1 4 K ( − m m + M r μ k + 1 + R μ k + 1 − M m + M r μ k − R μ k )     × ( − m m + M r μ k + 1 + R μ k + 1 − M m + M r μ k − R μ k ) (8)

The Euler Lagrange equations of motion of Reaction Mechanics [2] for the oscillating and center of mass motions are:

a ) d d t k ∂ S ∂ r ˙ μ k − ∂ S ∂ r μ k = 0 b ) d d t k ∂ S ∂ R ˙ μ k − ∂ S ∂ R μ k = 0 (9)

Substituting the Lagangians of Equation (6) and Equation (8) into Equation (2) and the resulting sum of Lagrangians into the Euler Lagrange equation of motion, Equation (9a) to obtain the oscillation equation of motion.

m R r ¨ μ k + M m + M K 2 ( M m + M r μ k + R μ k + m m + M r μ k − 1 − R μ k − 1 ) + m m + M K 2 ( m m + M r μ k − R μ k + M m + M r μ k − 1 + R μ k − 1 ) + m m + M K 2 ( M m + M r μ k + 1 + R μ k + 1 + m m + M r μ k − R μ k ) + M m + M K 2 ( m m + M r μ k + 1 − R μ k + 1 + M m + M r μ k + R μ k ) = 0 (10)

Performing a first step of collecting terms:

m R r ¨ μ k + K r μ k − M m K ( m + M ) 2 ( 2 r μ k + r μ k − 1 + r μ k + 1 ) + M − m M + m K 2 ( 2 R μ k − R μ k − 1 − R μ k + 1 ) = 0 (11)

Performing a second step of collecting terms:

m R r ¨ μ k + K r μ k − m R K m + M ( r μ k r μ k r ˙ μ k r ˙ μ k ) 2 r μ k − r μ k − 1 − r μ k + 1 τ 2 + M − m 2 ( m + M ) K ( r μ k r μ k r ˙ μ k r ˙ μ k ) 2 R μ k − R μ k − 1 − R μ k + 1 τ 2 = 0 (12)

where the delay time τ is approximately:

τ 2 = r μ k r μ k r ˙ μ k r ˙ μ k (13)

Substituting the Lagrangians of Equation (6) and Equation (8) into Equation (2) and the resulting sum of Lagrangians into the Euler Lagrange equation of motion, Equation (9b) to obtain the center of mass equation of motion.

( m + M ) R ¨ μ k + 1 2 K ( M m + M r μ k + R μ k + m m + M r μ k − 1 − R μ k − 1 ) + 1 2 K ( − m m + M r μ k + R μ k − M m + M r μ k − 1 − R μ k − 1 ) + 1 2 K ( − M m + M r μ k + 1 − R μ k + 1 − m m + M r μ k + R μ k ) + 1 2 K ( m m + M r μ k + 1 − R μ k + 1 + M m + M r μ k + R μ k ) = 0 (14)

Collecting terms:

( m + M ) R ¨ μ k + K τ 2 2 R μ k − R μ k − 1 − R μ k + 1 τ 2 + K τ 2 2 M − m m + M 2 r μ k − r μ k − 1 − r μ k + 1 τ 2 = 0 (15)

where the second and third terms have been multiplied and divided by the square of the delay time τ². Approximating the second differences divided by τ² such as

2 r μ k − r μ k − 1 − r μ k + 1 τ 2

by the second derivatives times the square τ² of the propagation delay time τ.

a ) r μ k + 1 − 2 r μ k + r μ k − 1 τ 2 ≈ r ¨ μ k b ) R μ k + 1 − 2 R μ k + R μ k - 1 τ 2 ≈ R ¨ μ k (16)

Substituting Equation (16) into Equation (12) and dividing by the reduced mass m_R.

( 1 − K τ 2 m + M ) r ¨ μ + K m R r μ + M − m 2 ( m + M ) K τ 2 m R R ¨ μ = 0 (17)

Substituting Equation (16) into Equation (15).

( 1 − K τ 2 m + M ) R ¨ μ − K τ 2 2 M − m ( m + M ) 2 r ¨ μ = 0 (18)

Solving Equation (18) for R ¨ μ and substituting the resulting equation into Equation (17) to obtain the oscillation equation of motion of the Delay Harmonic Oscillator.

[ ( 1 − K τ 2 m + M ) 2 + K 2 τ 4 4 m M ( M − m ) 2 ( m + M ) 2 ] r ¨ μ + K m R ( 1 − K τ 2 m + M ) r μ = 0 (19)

Solving Equation (17) for r ¨ μ and substituting the resulting equation into Equation (18) to obtain the center of mass equation of motion of the Delay Harmonic Oscillator.

[ ( 1 − K τ 2 m + M ) 2 + K 2 τ 4 4 m M ( M − m ) 2 ( m + M ) 2 ] R ¨ μ + K 2 τ 4 2 m M M − m m + M r μ = 0 (20)

The final results are Equation (19) and Equation (20), the Equations of Motion of the Delay Harmonic Oscillator.

In the limit when both masses are equal, m = M, the center of mass acceleration components R ¨ μ are equal to zero, and the center of mass, at most, moves at a constant velocity and the objects oscillate with an approximate modified angular frequency.

ω m = M 2 ≈ K m R ( 1 − K τ 2 m + M ) (21)

Note from Equation (13) that the delay time τ is not constant.

In the limit when the delay time τ goes to zero, the equations of motions 19 and 20 revert to the equations of motion of a Simple Harmonic Oscillator [4] .

a ) m R r ¨ μ + K r μ = 0 b ) R ¨ μ = 0 (22)

3. Conservation of Energy

In this conservative system, it is assumed that the energy is equal to the Hamiltonian H of the system. In order to calculate the Hamiltonian, we must first calculate the momentum component p μ of the system. The total momentum component p μ is:

p μ = ∂ S ∂ r ˙ μ + ∂ S ∂ R ˙ μ (23)

Substituting Equation (6) and Equation (8) into Equation (2) and the resulting equation into Equation (23),

p μ k = m R r ˙ μ k + ( m + M ) R ˙ μ k (24)

The Hamiltonian sum is equal to the Legender Transform of Sum of Lagrangians with respect to the velocities r ˙ μ k and R ˙ μ k [2] :

H = m R r ˙ μ k r ˙ μ k + ( m + M ) R ˙ μ k R ˙ μ k − L k − L k + 1 (25)

where the Lagrangian L_k at time step t_k is given by equation 6 and the Lagrangian L_k+1 at time step t_k+1 is given by Equation (8). Taking the derivative with respect to the time t_k of the Hamiltonian of Equation (25) using Equation (6) and Equation (8) one obtains:

d H d t k = m R r ¨ μ k r ˙ μ k + ( m + M ) R ¨ μ k R ˙ μ k     + 1 2 K ( M m + M r μ k + R μ k + m m + M r μ k − 1 − R μ k − 1 ) ( M m + M r ˙ μ k + R ˙ μ k )     + 1 2 K ( − m m + M r μ k + R μ k − M m + M r μ k − 1 − R μ k − 1 ) ( − m m + M r ˙ μ k + R ˙ μ k )     + 1 2 K ( M m + M r μ k + 1 + R μ k + 1 + m m + M r μ k − R μ k ) ( m m + M r ˙ μ k − R ˙ μ k )     + 1 2 K ( − m m + M r μ k + 1 + R μ k + 1 − M m + M r μ k − R μ k ) ( M m + M r ˙ μ k − R ˙ μ k ) (26)

By multiplying Equation (10) by r ˙ μ k and Equation (14) by and R ˙ μ k and adding the result one obtains an equation identical to Equation (26). Since Equation (10) and Equation (14) are equal to zero, the change of the Hamiltonian H with time is equal to zero, and therefore the Hamiltonian and thus the energy are constant.

a ) d H d t k = 0 b ) H = constant c ) Energy = constant (27)

4. Conclusions

In this paper, the propagation delay of a deformation traveling through the spring is considered, and therefore, is called the Delay Harmonic Oscillator. The Equations of Motion of two masses at the end of a spring through which a deformation propagates with a finite velocity are derived. The effect of the delay is modeled by assuming that the forces at both ends of the spring are not simultaneous. The difference in time of the occurrence of these forces is equal to the propagation delay of the deformation signal traveling through the spring. This effect couples the oscillatory and center of mass motions. But this derivation is valid when the propagation of the interaction through the spring wire takes less time than the period of the oscillating frequency. It demonstrates that if the masses are unequal, the center of mass of this system oscillates. Also, the effect of the finite propagation delay causes the oscillation to be non sinusoidal. In the case of zero delay, the motion of the masses reverts to the motion of a Simple Harmonic Oscillator. The effect of the propagation delay through the spring is rarely discussed in literature.

The Equations of Motion derived here are useful for the design of mechanical systems with coil springs and masses. The vehicle suspention [3] is the most obvious application. It has two masses connected by coil springs and vibration dampers. However, the vibration dampers are not part of this analysis. The upper mass consists of the body, chassis and engine and the lower mass consists of the axels, differential and wheels. The upper mass is much heavier than the lower mass. Therefore, the center of mass will vibrate. If the vibration dampers were not there, this would cause a large vibration of the upper mass.

If one uses the standard calculation neglecting the propagation delay through the spring, one would obtain that the two masses vibrate about the stationary center of mass.

Another example of this type of motion occurs when an electric motor is mounted on a platform that is suspended on elastic supports. Even if the rotation of the motor is well balanced, the platform may still vibrate because there are different masses on both sides of the elastic support.

The Equations of Motion of the Delay Harmonic Oscillator are derived in this paper.

Acknowledgements

I thank my wife Marlene Danzig Kornreich for her suggestions to the text and for making the text more understandable to a reader. I also thank her for proofreading and working together on the manuscript.

References

1. Achenbach, J.D. (1984) Wave Propagation in Elastic Solids. North Holland. Elsevier Science, Amsterdam.

2. Kornreich, P. (2017) Reaction Mechanics of Point Objects. Journal of Applied Mathematics and Physics, 5, 137-152. https://doi.org/10.4236/jamp.2017.51015

3. Lavanayal, N., Rao Sampath, P. and Redd Promod, M. (2014) Design and Analysis of a Suspension Coil Spring for Automotive Vehicle. Journal of Engineering Research and Applications, 4, 151-157.

4. Goldstein, H. (1981) Classical Mechanics. 2nd Edition, Addison-Wesley Publishing Company, Boston.