# Standard linear solid model

The standard linear solid (SLS) model, also known as the Zener model, is a method of modeling the behavior of a viscoelastic material using a linear combination of springs and dashpots to represent elastic and viscous components, respectively. Often, the simpler Maxwell model and the Kelvin–Voigt model are used. These models often prove insufficient, however; the Maxwell model does not describe creep or recovery, and the Kelvin–Voigt model does not describe stress relaxation. SLS is the simplest model that predicts both phenomena.

## Definition

Materials undergoing strain are often modeled with mechanical components, such as springs (restorative force component) and dashpots (damping component).

Connecting a spring and damper in series yields a model of a Maxwell material while connecting a spring and damper in parallel yields a model of a Kelvin–Voigt material.[1] In contrast to the Maxwell and Kelvin–Voigt models, the SLS is slightly more complex, involving elements both in series and in parallel. Springs, which represent the elastic component of a viscoelastic material, obey Hooke's Law:

${\displaystyle \sigma _{s}=E\varepsilon }$

where σ is the applied stress, E is the Young's Modulus of the material, and ε is the strain. The spring represents the elastic component of the model's response.[1]

Dashpots represent the viscous component of a viscoelastic material. In these elements, the applied stress varies with the time rate of change of the strain:

${\displaystyle \sigma _{D}=\eta {\frac {d\varepsilon }{dt}}}$

where η is viscosity of the dashpot component.

These elements are connected as shown on the picture:

Standard Linear Solid model

This model consists of two systems in parallel. The first, referred to as the Maxwell arm, contains a spring (${\displaystyle E=E_{2}}$) and dashpot (viscosity ${\displaystyle \eta }$) in series.[1] The other system contains only a spring (${\displaystyle E=E_{1}}$).

## Solving the model

In order to model this system, the following physical relations must be realized:

For parallel components: ${\displaystyle \sigma _{tot}=\sigma _{1}+\sigma _{2}}$, and ${\displaystyle \varepsilon _{tot}=\varepsilon _{1}=\varepsilon _{2}}$.[1]

For series components: ${\displaystyle \sigma _{tot}=\sigma _{1}=\sigma _{2}}$, and ${\displaystyle \varepsilon _{tot}=\varepsilon _{1}+\varepsilon _{2}}$.[1]

These relationships help relate the various stresses and strains in the overall system and the Maxwell arm:

${\displaystyle \sigma _{tot}=\sigma _{m}+\sigma _{S_{1}}}$

${\displaystyle \varepsilon _{tot}=\varepsilon _{m}=\varepsilon _{S_{1}}}$

${\displaystyle \sigma _{m}=\sigma _{D}=\sigma _{S_{2}}}$

${\displaystyle \varepsilon _{m}=\varepsilon _{D}+\varepsilon _{S_{2}}}$

where the subscripts ${\displaystyle m}$, ${\displaystyle D}$, ${\displaystyle S_{1}}$,and ${\displaystyle S_{2}}$ refer to Maxwell, dashpot, spring one, and spring two, respectively.

Using these relationships, their time derivatives, and the above stress-strain relationships for the spring and dashpot elements, the system can be modeled as follows:

${\displaystyle {\frac {d\varepsilon (t)}{dt}}={\frac {{\frac {E_{2}}{\eta }}\left({\frac {\eta }{E_{2}}}{\frac {d\sigma (t)}{dt}}+\sigma (t)-E_{1}\varepsilon (t)\right)}{E_{1}+E_{2}}}}$ [2]

The equation can also be expressed as:

${\displaystyle {\frac {d\varepsilon (t)}{dt}}=\left(E_{1}+E_{2}\right)^{-1}\cdot \left[{\frac {d\sigma (t)}{dt}}+{\frac {E_{2}}{\eta }}\sigma (t)-{\frac {E_{1}E_{2}}{\eta }}\varepsilon (t)\right]}$

The relaxation time, ${\displaystyle \tau }$, is different for each material and is equal to

${\displaystyle {\frac {\eta }{E_{2}}}=\tau }$

## Model characteristics

The standard linear solid model combines aspects of the Maxwell and Kelvin–Voigt models to accurately describe the overall behavior of a system under a given set of loading conditions. The behavior of a material applied to an instantaneous stress is shown as having an instantaneous component of the response. Instantaneous release of a stress also results in a discontinuous decrease in strain, as is expected. The shape of the time-dependent strain curve is true to the type of equation that characterizes the behavior of the model over time, depending upon how the model is loaded.

Although this model can be used to accurately predict the general shape of the strain curve, as well as behavior for long time and instantaneous loads, the model lacks the ability to accurately model material systems numerically.

## References

1. David Roylance, "Engineering Viscoelasticity" (October 24, 2001) http://ocw.mit.edu/courses/materials-science-and-engineering/3-11-mechanics-of-materials-fall-1999/modules/visco.pdf
2. ^ Krystyn J. Van Vliet, MIT course 3.032 Lecture, October 23, 2006 http://stellar.mit.edu/S/course/3/fa06/3.032/index.html