Viscoelasticity & TMD Design

Physical Intuition of Damping and Viscoelasticity

In mechanical and structural systems, energy is continually converted between kinetic and potential forms. Without a mechanism to dissipate this energy, a system excited at its natural frequency would experience unbounded vibration amplitudes, leading to catastrophic fatigue failure, severe acoustic radiation, and general noise issues. Damping is the physical process of removing mechanical energy from a vibrating system, typically converting it into heat.

In Noise, Vibration, and Harshness (NVH) engineering, we rely heavily on two primary damping techniques:

Material Damping (Viscoelasticity): Utilizing polymer materials that exhibit both viscous and elastic behaviors to dissipate energy throughout a structure's surface.

Structural Damping Devices (Tuned Mass Dampers): Attaching secondary mass-spring-damper systems that absorb and dissipate energy from a specific structural resonance mode.

The Physics of Viscoelastic Materials (VEMs)

Viscoelastic materials (such as rubbers, polyurethanes, and synthetic polymers) exhibit a unique phase lag between stress and strain.

•

Under purely elastic deformation (like a metal spring), stress and strain are perfectly in phase, meaning energy is stored during loading and fully recovered during unloading.

•

Under purely viscous deformation (like a fluid dashpot), stress is proportional to the strain *rate*, creating a $90^\circ$ phase lag. All input energy is dissipated as heat.

•

Viscoelastic materials fall in between. The stress response lags the applied strain by a phase angle $\delta$ ($0 < \delta < 90^\circ$). During cyclic loading, this phase lag creates a hysteresis loop in the stress-strain plane. The area enclosed by the loop is the energy dissipated per cycle.

Unconstrained vs. Constrained Layer Damping (ULD vs. CLD)

To damp a vibrating sheet metal panel (such as a car door, floor pan, or roof), VEMs are bonded directly to the metal sheet. There are two primary application architectures:

•

Unconstrained Layer Damping (ULD / Free Layer Damping):

A viscoelastic layer is bonded directly to the base structure. As the structure bends, the VEM undergoes extension and compression along the structural surface. Damping is achieved through the tensile/compressive strain of the VEM. ULD treatments are simple and light, but their damping performance is moderate.

•

Constrained Layer Damping (CLD):

A viscoelastic layer is sandwiched between the base structure and a thin, stiff sheet of metal (typically aluminum or steel) called the constraining layer. When the assembly bends, the stiff outer layers force the soft viscoelastic core to deform in shear. Because VEMs are highly efficient at dissipating energy under shear deformation, CLD treatments provide significantly higher damping capacity than ULD treatments for a given mass.

Core Concepts

1. The Complex Modulus Formulation

For a viscoelastic material undergoing harmonic excitation, we represent the stress-strain relationship using complex variables. If strain is $\epsilon(t) = \epsilon_0 e^{i\omega t}$, the stress response is $\sigma(t) = \sigma_0 e^{i(\omega t + \delta)}$. The ratio of stress to strain defines the complex Young's modulus $E^*$:

$$E^* = \frac{\sigma(t)}{\epsilon(t)} = \frac{\sigma_0}{\epsilon_0} e^{i\delta} = E' + i E''$$

where:

•

$E'$ is the storage modulus (representing elastic stiffness and stored energy).

•

$E''$ is the loss modulus (representing viscous dissipation and heat generation).

•

$\delta$ is the loss angle.

We define the loss factor $\eta$ as the ratio of the loss modulus to the storage modulus:

$$\eta = \tan \delta = \frac{E''}{E'}$$

This allows us to write the complex modulus as:

$$E^* = E' (1 + i \eta)$$

Similarly, for shear deformation, we define the complex shear modulus:

$$G^* = G' (1 + i \eta)$$

2. Environmental Dependency: Temperature and Frequency

Viscoelastic properties are highly sensitive to temperature and excitation frequency:

•

Glassy Region (Low Temp / High Freq): The polymer chains are locked. The storage modulus $E'$ is very high (the material is stiff and brittle), but the loss factor $\eta$ is low because the chains cannot slide to dissipate energy.

•

Transition Region: The polymer undergoes a glass transition. The storage modulus drops rapidly, and the loss factor $\eta$ reaches its peak. This is the optimal operating range for damping.

•

Rubbery Region (High Temp / Low Freq): The chains are highly mobile. The storage modulus $E'$ is low (soft and rubbery), and the loss factor is low.

Engineers use the Williams-Landel-Ferry (WLF) equation to map viscoelastic properties across different temperatures and frequencies using a frequency-temperature shift factor $a_T$, facilitating the creation of Reduced Frequency Nomograms.

3. Oberst vs. Ross-Kerwin-Ungar (RKU) Models

To predict the damping of treated structures, we must calculate the composite system loss factor $\eta_s$.

Unconstrained Layer Damping (Oberst Model)

For a two-layer beam (base metal sheet of thickness $h_1$ and VEM of thickness $h_2$), Oberst derived the composite loss factor $\eta_s$ as:

$$\eta_s = \frac{\eta_2 e_2 h_2^3 + 12 \eta_2 e_2 h_2 (1 + h_2)^2}{1 + e_2 h_2^3 + 3 e_2 h_2 (2 + h_2)^2}$$

where:

•

$e_2 = E_2' / E_1$ is the ratio of VEM to base metal storage modulus.

•

$h_2 = h_2 / h_1$ is the thickness ratio.

•

$\eta_2$ is the loss factor of the VEM.

Constrained Layer Damping (RKU Model)

For a three-layer sandwich beam (base plate, VEM core, constraining plate), the Ross-Kerwin-Ungar (RKU) theory defines a dimensionless shear parameter $g$ that represents the shear coupling between the outer sheets:

$$g = \frac{G_2' L^2}{\pi^2 E_3 h_3 h_2} \cdot \Phi$$

where $G_2'$ is the VEM shear modulus, $L$ is the bending wavelength, $E_3$ and $h_3$ are the modulus and thickness of the constraining layer, and $\Phi$ is a geometric factor. The system loss factor is:

$$\eta_s = \frac{\eta_2 g Y}{1 + (2+Y)g + (1+\eta_2^2)(1+Y)g^2}$$

where $Y$ is the geometric stiffness parameter (or structural coupling parameter), representing the bending stiffness ratio of the uncoupled to coupled system.

4. Tuned Mass Damper (TMD) Design & Optimization

A Tuned Mass Damper (TMD) is a device consisting of a mass, spring, and damper attached to a primary structure to reduce its vibration amplitude at a specific resonance frequency.

                  Primary Force F(t)
                         │
                         ▼
                 ┌───────────────┐
                 │ Primary Mass  │ ◄─── Displacement x_s(t)
                 │     M_s       │
                 └───────────────┘
                   │           │
           Spring  └───[C_d]───┘ Damper
            K_d         │
                        ▼
                 ┌───────────────┐
                 │   TMD Mass    │ ◄─── Displacement x_d(t)
                 │     m_d       │
                 └───────────────┘

Den Hartog's Optimum Tuning under Harmonic Force

For a primary SDOF structure ($M_s, K_s$) excited by a harmonic force, Den Hartog derived the optimum parameters for a TMD ($m_d, k_d, c_d$) based on the mass ratio $\mu$:

$$\mu = \frac{m_d}{M_s}$$

Optimum Tuning Frequency Ratio ($g_{opt}$):

The ratio of the TMD natural frequency $\omega_d = \sqrt{k_d / m_d}$ to the primary structure natural frequency $\omega_s = \sqrt{K_s / M_s}$ must be:

$$g_{opt} = \frac{1}{1 + \mu}$$

Optimum TMD Damping Ratio ($\zeta_{opt}$):

The optimum damping ratio for the TMD spring-damper is:

$$\zeta_{opt} = \sqrt{\frac{3\mu}{8(1 + \mu)}}$$

These tuning parameters minimize the peak response of the primary structure, splitting the single resonance peak into two lower, equal-height peaks.

High-Quality Schematic Diagrams

Diagram 1: ULD vs. CLD Viscoelastic Deformation Mechanisms

This diagram compares how Unconstrained Layer Damping (ULD) and Constrained Layer Damping (CLD) deform under bending. In ULD, the viscoelastic layer undergoes extension/compression. In CLD, the VEM is sandwiched, forcing it to undergo severe shear strain.

Diagram 2: Tuned Mass Damper Effect on SDOF Frequency Response

This diagram illustrates the frequency response function (FRF) of the primary structure before and after adding a tuned mass damper. It shows how the single high-amplitude peak splits into two much lower peaks.