Technical Report: Coefficients Representing Damping Characteristics (2)
3. The effect of damping
Up to this point, we have explained how to determine the damping ratio and loss factor, assuming that the amplitude decreases exponentially. Now, we will briefly explain what the damping ratio means physically in actual vibrations. The loss factor and Q value can be easily converted from the damping ratio, so here we will focus on the damping ratio.
3-1 The effect of damping in free vibration
Free vibration is vibration that occurs "without any external force applied." Without any external force, it would remain stationary indefinitely, but if initial conditions such as an initial displacement and initial velocity are applied, it will begin to vibrate. For example, consider the spring-mass model shown in Figure 4. If you initially pull on a mass m to impart a displacement (initial displacement) to a spring k, and then suddenly release it, it will begin to vibrate; this is free vibration.
-
Figure 4: Free vibration model
If there is no damping force c, free vibration will continue indefinitely, and the vibration frequency ω0 in this case is expressed by the following equation.
(11)
ω0 is called the natural frequency. As is known empirically, this free vibration does not actually last forever, and a damping force c acts on it, causing the amplitude to gradually decrease as illustrated in Figure 1, until it eventually comes to a resting state. At this time, the behavior differs depending on whether the value of c is greater than or less than the following equation c.
(12)
c c is called the critical damping rate, and the ratio of c c to c is ζ (damping ratio), which is the subject of this paper.
(13)
As mentioned above, the amplitude of free vibration changes significantly depending on the value of ζ. An example is shown in Figure 5.
-
Figure 5. Response to free vibrations
The amplitude of damped free vibration when ζ < 1 is given by the following equation:
(14)
Here,
x 0: Initial displacement
v 0: initial speed
Furthermore, ωd is called the natural frequency of the damped system and is expressed by the following equation.
(15)
ωd is somewhat smaller than ω0, but in real-world vibration systems, the value of ζ is small, so ωd is close to ω0. As can be seen from equation (14), the behavior of a damped vibration system is determined by the initial conditions and the damping ratio ζ. Figure 5 shows the response for different damping ratios ζ when the initial velocity is 0 and the initial displacement is 1, and it can be seen that the behavior differs greatly depending on the damping ratio ζ.
Note that when ζ ≥ 1, the calculation cannot be performed using equation (14), and a different equation is required. While that calculation formula is omitted here, Figure 5 shows the response for comparison. Incidentally, the state where ζ = 1 is called critical damping, ζ > 1 is called overdamping, and 1 > ζ > 0 is called underdamping. In overdamping and critical damping, the motion is damped without oscillation. In Figure 5, ζ = 1 (critical damping) is emphasized for clarity, but this merely indicates the boundary between oscillation and non-oscillation, and does not imply that critical damping is particularly important.
The response for ζ = 0.707 (= 1/√2) is also shown, which is an important value in the following steady-state oscillation. Furthermore, although there is some overshoot (undershoot), this value is often used in control systems and other applications where response time is important, as it results in the shortest settling time (the time it takes for the response to converge to within 5% of the target value).
3-2 Effects of Damping in Steady-State Vibrations
Next, let's consider the case where a continuous external force is applied to a free-vibrating system.
Vibration caused by an external force is called forced vibration. However, vibration in a state where the external force is a sine wave and a sufficient amount of time has elapsed since the force was applied (steady state) is called steady vibration. In contrast, the process from the application of the external force to reaching the steady state is called the transient state, which will be explained in the next section.
As shown in Figure 6, consider a model in which an excitation force F cos ωt is applied to a one-degree-of-freedom vibration system.
First, let's review the symbols and formulas.
-
-
Figure 6: Forced Excitation Vibration Model
(11)
Since the excitation force is a repeating force of frequency ω, the steady-state vibration driven by it will also be a vibration of the same frequency. However, the amplitude and phase will be different, and here we will find the amplitude and phase.
The equation of motion for the system in Figure 6 is given by the following equation, and by solving this equation, the amplitude and phase of the steady-state oscillation can be determined.
(16)
The equation for steady-state oscillation is:
(17)
Therefore, the amplitude x a and phase φ are expressed by the following equations.
(18)
(19)
In equation (18), F/k is the static displacement when a static force F is applied, so if we let this be x s, then equation (18) becomes:
(20)
This can be expressed as follows: In other words, the amplitude of steady-state vibration is expressed as a frequency response function with static displacement x s, natural frequency ω 0, and damping ratio ζ.
Furthermore, equation (20) was transformed;
(21)
This is called the amplitude factor. Figure 7 shows the amplitude factor plotted on a logarithmic scale with ω/ ω0 on the horizontal axis and the amplitude factor on the vertical axis. This shows that steady-state vibrations resonate around ω0, and that the amplitude factor changes significantly depending on the damping ratio ζ.
-
Figure 7 Amplitude Magnification
To summarize the characteristics of amplitude magnification:
ζ < 1/√2 At that time, ω r = √1-2 ζ 2 ω o It has one resonance point with a resonant frequency of [frequency]. ω r teeth ζ The larger the value, the lower the value, but a low damping system, i.e. ζ When it is small (approximately) ζ < 0.05) ω r ≒ ω 0 In other words, it resonates at its natural frequency. ζ ≥ 1/√2 Then it won't resonate.
The amplitude magnification is maximized at the resonant frequency.
(22)
When ζ is small:
(23)
This becomes equal to the Q-factor. When ζ is small, i.e., when the resonance is sharp, the Q-factor is often used.
At frequencies below the resonance point, the amplitude magnification asymptotically approaches 1.
At frequencies higher than the resonance point, the amplitude multiplier asymptotically approaches 1/(ω/ω 0) 2, i.e., a slope of -40 dB/decade.
Equation (19) represents the phase difference between the excitation force and the steady-state vibration. Graphing this gives us Figure 8.
-
Figure 8 Phase lag of steady-state vibration
What we will learn from this is:
In the frequency range where the excitation force frequency is lower than ω0, the phase lag of the steady-state vibration asymptotically approaches 0 degrees, meaning it vibrates with a phase slightly delayed from the excitation force.
In the frequency range above ω0, it asymptotically approaches 180 degrees, meaning it oscillates with a phase close to the opposite phase to the excitation force.
At ω = ω0, the oscillation is delayed by 90 degrees, or 1/4 of a period.
When ζ is small, the phase changes abruptly around ω0, and as ζ increases, the change becomes more gradual.
3-3 Damping effect in transient vibrations
In the previous section, we considered a state where a sufficient amount of time had elapsed since an external force was applied to a steady state. Next, let's consider the transient state, that is, the state from the time an external force is applied until a steady state is reached.
Considering the vibration system in Figure 6, its equation of motion is given by equation (24). However, for simplicity, we will assume the external force is F sin ω t, and the initial condition is complete rest, meaning the initial displacement and initial velocity are zero.
(24)
Since this system is linear, the superposition principle holds, and the solution is obtained in the form of the sum of the vibration component due to the external force and the free vibration component, as we have seen so far.
(25)
Here;
(26)
(27)
The first term of equation (25) is the free vibration component, which decays over time until only the steady-state vibration component (the second term) remains. Figures 9, 1 to 4, show this process graphically. Here, with ζ = 0.01, the changes in the transient response are shown when ω/ ω0 is varied to 0.1, 0.9, 1.0, and 2.0.
When ω/ ω0 is small, free vibrations are simply superimposed on steady vibrations, and these free vibrations decay over time, transitioning to steady vibrations.
As ω/ ω0 approaches 1, that is, as the excitation frequency approaches the natural vibration frequency, the amplitude increases and a humming sound is produced.
ω/ ω₀ = 1, that is, when the excitation frequency matches the natural vibration frequency, the amplitude increases almost proportionally to time, reaching a very large amplitude, i.e., a resonance state.
When ω/ω 0 > 1, the amplitude decreases, but the waveform becomes more complex.
Here, we used a small value of ζ = 0.01 to clearly illustrate the transient state. However, a larger ζ value results in faster deceleration of the free vibration and a smaller amplitude of the steady-state vibration. This amplitude is shown in Figure 7. Conversely, a smaller ζ value results in a prolonged transient state and a longer period of instability. The amplitude of the steady-state vibration also increases, and particularly at frequencies around ω/ ω₀ = 1, even small vibrations will gradually increase in amplitude over time, growing into very large vibrations. (Note that Figures 9-1 to 9-4 have different vertical axis scales.)
ζ = 0.01
ω/ω0 = 0.1
ζ = 0.01
ω/ω0 = 0.9
ζ = 0.01
ω/ω0 = 1
ζ = 0.01
ω/ω0 = 2
Technical Report
FFT Analysis Solution
-
Engineering ServicesAcoustic vibration characteristics evaluation (JIS/ISO)
-
Portable 2/4 channel FFT analyzerCF-9200A/9400A
-
FFT ComparatorCF-4700A
-
Frequency response measurement softwareServo analyzer OS-4100
-
Sound and Vibration Analysis SystemO-Solution DS-5000
-
Experimental modal analysis system softwareMEscope
-
Engineering ServicesModel-based development support
-
Measurement and analysis softwareO-Solution
-
FFT analysis functionO-Solution OS-0522
-
Tracking analysis functionO-Solution OS-0523
-
Octave analysis functionO-Solution OS-0524
-
Engineering ServicesContract measurement and consulting
-
Engineering ServicesAutomotive powertrain performance evaluation