Damping in Structural Dynamics: Theory and Sources
March 14, 2019by Henrik Sönnerlind
If you strike a bowl made of glass or metal, you hear a tone with an intensity that decays with time. In a world without damping, the tone would linger forever. In reality, there are several physical processes through which the kinetic and elastic energy in the bowl dissipate into other energy forms. In this blog post, we will discuss how damping can be represented, and the physical phenomena that cause damping in vibrating structures.
How Is Damping Quantified?
There are several ways by which damping can be described from a mathematical point of view. Some of the more popular descriptions are summarized below.
One of the most obvious manifestations of damping is the amplitude decay during free vibrations, as in the case of a singing bowl. The rate of the decay depends on how large the damping is. It is most common that the vibration amplitude decreases exponentially with time. This is the case when the energy lost during a cycle is proportional to the amplitude of the cycle itself.
Let’s start out with the equation of motion for a system with a single degree of freedom (DOF) with viscous damping and no external loads,
After division with the mass, m, we get a normalized form, usually written as
Here, \omega_0 is the undamped natural frequency and \zeta is called the damping ratio.
In order for the motion to be periodic, the damping ratio must be limited to the range 0 \le \zeta < 1. The amplitude of the free vibration in this system will decay with the factor
where T0 is the period of the undamped vibration.
Decay of a free vibration for three different values of the damping ratio.
Another measure in use is the logarithmic decrement, δ. This is the logarithm of the ratio between the amplitudes of two subsequent peaks,
where T is the period.
The relation between the logarithmic decrement and the damping ratio is
Another case in which the effect of damping has a prominent role is when a structure is subjected to a harmonic excitation at a frequency that is close to a natural frequency. Exactly at resonance, the vibration amplitude tends to infinity, unless there is some damping in the system. The actual amplitude at resonance is controlled solely by the amount of damping.
Amplification for a single-DOF system for different frequencies and damping ratios.
In some systems, like resonators, the aim is to get as much amplification as possible. This leads to another popular damping measure: the quality factor or Q factor. It is defined as the amplification at resonance. The Q factor is related to the damping ratio by
Another starting point for the damping description is to assume that there is a certain phase shift between the applied force and resulting displacement, or between stress and strain. Talking about phase shifts is only meaningful for a steady-state harmonic vibration. If you plot the stress vs. strain for a complete period, you will see an ellipse describing a hysteresis loop.
You can think of the material properties as being complex-valued. Thus, for uniaxial linear elasticity, the complex-valued stress-strain relation can be written as
Here, the real part of Young’s modulus is called the storage modulus, and the imaginary part is called the loss modulus. Often, the loss modulus is described by a loss factor, η, so that
Here, E can be identified as the storage modulus E’. You may also encounter another definition, in which E is the ratio between the stress amplitude and strain amplitude, thus
in which case
The distinction is important only for high values of the loss factor.
An equivalent measure for loss factor damping is the loss tangent, defined as
The loss angle δ is the phase shift between stress and strain.
Damping defined by a loss factor behaves somewhat differently from viscous damping. Loss factor damping is proportional to the displacement amplitude, whereas viscous damping is proportional to the velocity. Thus, it is not possible to directly convert one number into the other.
In the figure below, the response of a single-DOF system is compared for the two damping models. It can be seen that viscous damping predicts higher damping than loss factor damping above the resonance and lower damping below it.
Comparison of dynamic response for viscous damping (solid lines) and loss factor damping (dashed lines).
Usually, the conversion between the damping ratio and loss factor damping is considered at a resonant frequency, and then \eta \approx 2 \zeta. However, this is only true at a single frequency. In the figure below, a two-DOF system is considered. The damping values have been matched at the first resonance, and it is clear that the predictions at the second resonance differ significantly.
Comparison of dynamic response for viscous damping and loss factor damping for a two-DOF system.
The loss factor concept can be generalized by defining the loss factor in terms of energy. It can be shown that for the material model described above, the energy dissipated during a load cycle is
where \varepsilon_a is the strain amplitude.
Similarly, the maximum elastic energy during the cycle is
The loss factor can thus be written in terms of energy as
This definition in terms of dissipated energy can be used irrespective of whether the hysteresis loop actually is a perfect ellipse or not — as long as the two energy quantities can be determined.
Sources of Damping
From the physical point of view, there are many possible sources of damping. Nature has a tendency to always find a way to dissipate energy.
Internal Losses in the Material
All real materials will dissipate some energy when strained. You can think of it as a kind of internal friction. If you look at a stress-strain curve for a complete load cycle, it will not trace a perfect straight line. Rather, you will see something that is more like a thin ellipse.
Often, loss factor damping is considered a suitable representation for material damping, since experience shows that the energy loss per cycle tends to have rather weak dependencies on frequency and amplitude. However, since the mathematical foundation for loss factor damping is based on complex-valued quantities, the underlying assumption is harmonic vibration. Thus, this damping model can only be used for frequency-domain analyses.
The loss factor for a material can have quite a large variation, depending on its detailed composition and which sources you consult. In the table below, some rough estimates are provided.
|Material||Loss Factor, η|
Loss factors and similar damping descriptions are mainly used when the exact physics of the damping in the material is not known or not important. In several material models, such as viscoelasticity, the dissipation is an inherent property of the model.
Friction in Joints
It is common that structures are joined by, for example, bolts or rivets. If the joined surfaces are sliding relative to each other during the vibration, the energy is dissipated through friction. As long as the value of the friction force itself does not change during the cycle, the energy lost per cycle is more or less frequency independent. In this sense, the friction is similar to internal losses in the material.
Bolted joints are common in mechanical engineering. The amount of dissipation that will be experienced in bolted joints can vary quite a lot, depending on the design. If low damping is important, then the bolts should be closely spaced and well-tightened so that macroscopic slip between the joined surfaces is avoided.
A vibrating surface will displace the surrounding air (or other surrounding medium) so that sound waves are emitted. These sound waves carry away some energy, which results in the energy loss from the point of view of the structure.
A plot of the sound emission in a Tonpilz transducer.
Often, a small component is attached to a larger structure that is not part of the simulation. When the component vibrates, some waves will be induced in the supporting structure and carried away. This phenomenon is often called anchor losses, particularly in the context of MEMS.
Even with pure elastic deformation without dissipation, straining a material will change its temperature slightly. Local stretching leads to a temperature decrease, while compression implies a local heating.
Fundamentally, this is a reversible process, so the temperature will return to the original value if the stress is released. Usually, however, there are gradients in the stress field with associated gradients in the temperature distribution. This will cause a heat flux from warmer to cooler regions. When the stress is removed during a later part of the load cycle, the temperature distribution is no longer the same as the one caused by the onloading. Thus, it is not possible to locally return to the original state. This becomes a source of dissipation.
The thermoelastic damping effect is mostly important when working with small length scales and high-frequency vibrations. For MEMS resonators, thermoelastic damping may give a significant decrease of the Q factor.
Sometimes, a structure contains intentional discrete dampers, like the shock absorbers in a wheel suspension.
Such components obviously have a large influence on the total damping in a structure, at least with respect to some vibration modes.
A particular case where much effort is spent on damping is in civil engineering structures in seismically active areas. It is of the utmost importance to reduce the vibration levels in buildings if hit by an earthquake. The purpose of such dampers can be both to isolate a structure from its foundation and to provide dissipation.
Read the follow-up to this blog post here: How to Model Different Types of Damping in COMSOL Multiphysics®
- COMSOL Now
- Today in Science
- COMSOL Now
- Today in Science