Advanced Dynamics and Vibrations: Lagrange’s equations applied to dynamic systems
Analytical Mechanics – Lagrange’s Equations
Up to the present we have formulated problems using newton’s laws in which the main disadvantage of this approach is that we must consider individual rigid body components and as a result, we must deal with interaction forces that we really have no interest in. These forces come from constraints which one part of the system puts on another.
An alternate approach referred to as analytical mechanics, consider the system as a whole and formulate the equations of motion starting from scalar quantities – the kinetic and potential energies and an expression of the virtual work associated with non-conservative forces. If relies heavily on the concept of virtual displacements.
Any set of numbers that serve to specify the configuration of the system are examples of generalized coordinates. To these we can add constraints. The number of DOF is then the number of generalized coordinates minus the number of independent equations of constraint.
e.g a point moving on the surface of a sphere of radius R with center must obey:
There are 3 coordinates and 1 equation of constraint. Therefore it has 2 DOF. If we used two angles, say & then we would not need the constraint.
These constraints are classified as:
If the constraint can be written in finite form (as above) they are called holonomic if not (then in differential form) then they are non-holonomic.
If they are equalities then they are bilateral and unilateral if inequalities
If there is no explicit dependence on time then they are scleronomic and if there is explicit time dependence they are rheonomic.
Therefore, 2 DOF holonomic, scleronomic bilateral system.
however if we allow the point of attachment to move (in the x-y plane)
Vertical displacement & Virtual Work
The work done during the motion is
where F & dr are vectors and is a scalar. For a particle where may be the sum of several forces. Therefore,
where this is a true differential of the kinetic energy.
If the particle moves from positive to then
if depends only on ,
where is a scalar function that depends only on the position and not on or time. is term the potential energy.
The total energy is constant.
When we have a conservative system we can use Rayleigh’s to estimate the natural frequencies.
Consider a simple problem and use visual work to solve.
Two bars of length . What is the equilibrium position
Consider a rigid beam (uniform) of mass located by a uniform time dependent load and supported by a spring & dashpot as shown
Using determine the equation of motion.
Virtual work of force applied Virtual work of metric forces
We can use virtual work to calculate the energy in a multi-degree of freedom system.
For a displacement of only
and similarly for , alone
So for the combination (assuming superposition)
The work done by force in is energy stored due to that force:
If there are n forces and the displacements are (i = 1,2, … n) then the total energy is:
We can also show that the kinetic energy is
& are examples of so-called quadratic forms. Quadratic forms are said to be positive definite if they are never negative and only zero when the variables are zero ()
They are said to be positively semi-definite if they are never negative but can be zero for non-negative variables.
This is useful to know because if we have a stiffness matrix that is positive semi-definite it means there is a rigid body mode with zero natural frequency [m] is positive definite
For example a submarine:
has 2 rigid body modes
Principle of Virtual Work
We will consider a system of particles moving in a 3D space and define the virtual displacements as for particle 1 for particle 2 etc.
These are not true displacements but small variations in the so-called generalized coordinates. These are imagined displacements (slight variations) which are at the same time and which satisfy the constraints given by
Then for virtual changes in these coordinates:
Expanding this in a Taylor series gives:
So that there are arbitrary displacements
On each of these particles the forces acting are
Where the resultant force is composed of the external forces and the constraint forces.
When all the particles are in static equilibrium
where this represents the virtual work performed by the resultant force on the particle.
If we rule out work by the constraint forces ( as the constraint force is perpendicular to the motion), therefore cannot handle friction.
Consider two particles held by a rigid massless rod.
But because the rod is rigid
Therefore for the constraint forces
and the Principal of Virtual Work say
that the work performed by the applied forces through infinitesimal virtual displacements compatible with the constraints is zero.
To make the virtual work equations more useful we write it in terms of the generalized coordinates
The generalized forces. These are not necessarily forces as they can be moments. It is only necessary that the quantity has the units of work.
For the situation in which we have a conservative system the work can be expressed differently.
as we can select the virtual displacement arbitrarily (i.e. only one non-zero at a time)
This states that has a stationary value at the equilibrium position (can use this for stability consideration). To extend this to dynamical situations we use D’Alembert’s principle.
Where the are the inertia forces. The virtual work for a system of particles is:
and since we showed the virtual work of the interned forces is zero.
This is called the generalized D ‘Alembert’s principle and is really not useful as it is still in vector form. To get it in a scalar form we rewrite the second term in terms of any generalized coordinates and this leads to Lagrange’s equations of motion.
Consider a coordinate transformation where radius vector is written in terms of instead .
Note: that since do not depend explicitly on time that.
now consider a virtual change in
and apply this to the second term of D ‘Alembert’s principle
JASON SECTION STARTS HERE
Consider a typical term in the bracket:
Now interchange the order of differentiation in the last term.
As a result, the total second term becomes:
As is the kinetic energy of the entire system. To complete the conversion, we must reunite the forces.
And we call:
Some of the forces will be conservative (derivable from a potential and some non-conservative).
So that the rewritten D’Alembert’s principle becomes:
And since the displacements are arbitrary and independent.
However, since the potential energy does not depend on the velocities, this can be written as:
- Does not have to be only for “small” oscillations.
- Does not matter whether the conservative forces are included in the .
- Sometimes the are able to be put into a “potential like” function.
- Can include reaction if wanted using Lagrange multipliers.
- Determine the using virtual work.
Find the equations of motion:
NOTE: Sometimes write the damping forces as a dissipation function.
Find the equations of motion:
The potential energy is:
And the damping force is the same. Note that:
The equation of motion comes from:
For small deflection and :
The uniform rod of length and mass carries a slider of mass which is attached by a spring of stiffness . The spring is unstretched at .
Also, note that:
The equilibrium positions are:
We could now study the stability of these configurations.
Example with External Forces
Using , we can say:
Where is the datum. Now we must find the generalized forces. For a virtual displacement only:
For a virtual displacement only:
Then, note that:
Therefore, (1) becomes:
And (2) becomes:
For small angles these become:
If we have support motion instead of forcing,
then we cannot neccesairly use a potential, so now give the system a and . For a :
as no work is done for alone. Thus, T is the same, and:
For small motions:
Consider the double pendulum:
Consider only small angles. Thus:
Now try a different set of generalized coordinates for the same problem.
If we consider small angles:
For mode shapes:
These are different values of ratios than in the first case. However, when plotted:
Now use the newtonain approach to the same problem:
Note both mass and stiffness matrices are non-symmetric!
Calculating the eigenvalues gives the same result and the mode shapes are identical with the other 2 formulations.
NOTE: These sets(first and third formulation) are linearly dependent.
Lagrange’s Equation for Small Oscillations
Our main interest is in the motion of MDOF systems in the neighborhood of equilibrium positions. WLOG we assume the equilibrium position is and that the displacements are sufficiently small that the linear force- displacement and force-velocity relations hold. This means that the generalized coordinates are their time derivatives appear in the differential equations to only the first power.
And are constants.
Set as the generalized masses.
Similarly we can write the potential energy as and consider a Taylor expansion of V about the equilibrium configuration.
Where the partial derivatives are all evaluated at .
As the equilibrium point is a condition in which:
Where are the stiffness coefficients. If we now use Lagrange’s equation the equations of motion become:
While this is very straightforward, we can use this result to obtain approximate formulations of the equations of motion as continuous systems it replaces them with a degree of freedom system. It can also extend the technique to find approximate response to a specific forcing function. This is called the assumed modes method (AMM).
For discrete system it replaces them with a smaller number of DOF, while for continuous system it replaces them with an degree of freedom system. We can also easily extend it to a forced system.
Then the method assumes:
Where are trial functions and are the generalized coordinates.
needs only satisfy the geometric boundary conditions.
For continuous systems:
While for discrete systems
The potential energy can be determined is a similar manner.
Consider a generalized beam:
This will allow calculation of a consistent mass and stiffness matrices.
Example: Discrete System
We want to estimate the lowest natural frequency using assumed modes method.
Now calculate the consistent mass and stiffness using this assumed shape:
What if we choose .
Then and are the same.
So we have replace the system by:
Now calculate the ‘s.
The answer is dependent on the assumption of mode shape
EXAMPLE: Use static deflection as mode.
Add masses to the diagram
It is sometimes useful to consider these approximate solutions to finding the eigenvalues and eigenvectors. The UAMM gives upper bounds as do most techniques. We can get lower bounds as well.
Consider again the equations and motion.
To get to an eigenvalue problem we can do it in two ways.
- Premultiply by
- Premultiply by
either gives the same eigenvalues and equivalent eigenvectors.
One can use these to find the eigenvalues/vectors through matrix iteration.
One clever idea that gives a lower bound is Dunkerleys formula. Consider the Dynamic Matrix formulation
and expand the determinant of the coefficients
Now consider that we have the roots of this polynomial . This equation is equivalent to (factored)
and expand this to get
Now equate the coefficients
However the largest term is the first one and often dominates
A special case occurs when the mass matrix is diagonal. Then tr D is
and we need calculate only
In this case there is a physical interpretation of the expressions
namely that is the natural frequency of the system with only as the mass (all other masses set to zero).
This is also sometimes shown as Dunkerley’s formula.
Consider previous example
- Apply unit load to
Apply unit load to
Apply unit load to
Consider an example of a continuous system
Therefore the eigenvalue problem is
The mode shape is
What if we apply some dynamic loading to the structure? Then we must calculate the generalized form for our model. Consider a virtual displacement
Consider the previous example
Whirling of Rotating Shafts
Rotating shafts tend to bow out at certain speeds and whirl in a complicated manner. Whirling is defined as the rotation of the bent shaft and the line of centers of the bearing. This phenomenon results from various causes as mass unbalance, hysteresis dangers, gyroscopic effects, fluid friction, etc. This whirling can take place in the same or opposite direction as that of the rotation of the shaft. The whirling speed may or may not be equal to the rotation speed.
Let us assume a simple system
is the stiffing of the system to lateral of flection.
We will assume the shaft to be rotating at an instant angular velocity and OS to be whirling at angular speed
These are nonlinear equations and there are more than one solution.
Consider the simplest case of synchronous whirl i.e. try
This looks exactly the same as rotating inbalance excitation
So that synchronous whirl occurs when
Self Excited Vibrations
Almost all of the vibrations considered so far were either free or forced vibrations. A fundamentally different class are the so-called self-excited vibrations.
“In a self-excited vibration the alternating force that sustains the motion is created or controlled by the motion itself ; when the motion stops the alternating force disappears”
In a forced vibration the sustaining alternating force exists independently of the motion and persists even when the vibratory motion is stopped.
- Use vibration phenomena for the display of classic ideas that ask the question “Why does that occur”
- Search kinetic art
- Naum Gabo Stationary Wave (1919)