APMR
Acoustical Porous Material Recipes

Mechanical damping

The mechanical damping tends to reduce the amplitude of a vibrating structure by molecular interactions inside the solid phase of a porous material. Mainly two models accounting for this damping are used, the the specific damping, usually characterized with letter ξ and the hysteretic one, usually characterized by letter η.

Hysteretic damping

This model of damping is only defined in Fourier's domain.

(mω2+k(1+iη))u~(ω)=f~(ω) (1)

dividing eq. (1) by the stiffness of the spring, k, leads to the following frequency response (ration of applied force over resulting displacement):

H~=1/kmkω2+(1+iη) (2)

Introducing the angular frequency of the undamped spring-mass, ω0=k/m, in the previous expression leads to:

H~=1/k1(ωω0)2+iη (3)

From this last equation, one can deduce the magnitude of the frequency response:

|H~|=1/k(1(ωω0)2)2+η2 (4)

The resonance angular frequency of the damped system (3), ωr, is obtained solving:

|H~|2ω=0 (5)

The left hand side of the previous equation is calculating from eq. (4):

|H~|2ω=1/k2(2(2ωω02)(1(ωω0)2))[(1(ωω0)2)2+η2]2 (6)

The right hand side of eq. (6) vanishes when ω, equals ω0 (for ω0).
Thus, the resonance angular frequency of the damped system (3), ωr, is equal to resonance angular frequency of the undamped system ω0 as the driving angular frequency, ω, does not appear in the imaginary part of the frequency response (3).

Specific damping

The so-called fundamental equation of the damped harmonic oscillator in the time domain can be written as:

md2u(t)dt2+cdu(t)dt+ku(t)=f(t) (7)

In Fourier's domain, this expression becomes:

(mω2+icω+k)u~(ω)=f~(ω) (8)

From equation (8) the frequency response (ratio of force applied over resulting displacement) is:

H~=1/k1(ωω0)2+ickω (9)

Introducing the specific damping, ξ, as c/k=2ξ/ω0 one rewrites eq. (9) as:

H~=1/k1(ωω0)2+2iξωω0 (10)

From eq. (10), the expression of magnitude of the frequency response is deduced:

|H~|=1/k(1(ωω0)2)2+4ξ2(ωω0)2 (11)

The resonance angular frequency ωr is calculated using the same reasoning as for the hysteric case above, i.e. finding ω so that:

|H~|2ω=0 (12)

The left hand side of the previous equation is calculating from eq. (11):

|H~|2ω=1/k2(2(2ωω02)(1(ωω0)2)+8ξ2ωω02)[(1(ωω0)2)2+4ξ2(ωω0)2]2 (13)

Except for frequency 0 Hz, ωr=ω0(12ξ2)1/2 because the driving frequency appears in the imaginary part. With this model, damping appears in the resonance angular frequency.
Note that ξ is well defined only for the harmonic oscillator because ξ depends on ω0 and for a system with multiple freedom degrees ξ must be defined for each frequency mode.

Comparison between specific and hysteretic damping

ξ and η have both non-dimensional quantities but the two FRF have different behaviours regarding first their resonance frequency.

Damping model Magnitude of frequency response
η |H~|=1/k(1(ωω0)2)2+η2
ξ |H~|=1/k(1(ωω0)2)2+4ξ2(ωω0)2



At ω=ω0 it is possible to write:

Damping model Magnitude of frequency response at ω=ω0
η 1/kη2
ξ 1/k4ξ2



If the two frequency reponses must be the same at ω=ω0, then η=2ξ. This is true at only one frequency. In practice, these two models are similar for very small damping.



The content of this page is copyleft under :
the creative commons license Attribution 3.0 Unported (CC BY 3.0).
Christophe Heinkelé, Luc Jaouen (@ljaouen), ISSN 2606-4138.
Site hosted by Matelys.