title: Standard linear solid model
text: # Standard linear solid model
  
  The standard linear solid (SLS), 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. 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:
  
      $\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.
  
  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:
  
      $\sigma _{D}=\eta {\frac {d\varepsilon }{dt}}$
  
  where η is viscosity of the dashpot component.
  
  ## Solving the model
  
  In order to model this system, the following physical relations must be realized:
  
  For parallel components: $\sigma _{tot}=\sigma _{1}+\sigma _{2}$, and $\varepsilon _{tot}=\varepsilon _{1}=\varepsilon _{2}$.
  
  For series components: $\sigma _{tot}=\sigma _{1}=\sigma _{2}$, and $\varepsilon _{tot}=\varepsilon _{1}+\varepsilon _{2}$.
  
  ### Maxwell representation
  
  This model consists of two systems in parallel. The first, referred to as the Maxwell arm, contains a spring ($E=E_{2}$) and dashpot (viscosity $\eta$) in series. The other system contains only a spring ($E=E_{1}$).
  
  These relationships help relate the various stresses and strains in the overall system and the Maxwell arm:
  
  $\sigma _{tot}=\sigma _{m}+\sigma _{S_{1}}$
  
  $\varepsilon _{tot}=\varepsilon _{m}=\varepsilon _{S_{1}}$
  
  $\sigma _{m}=\sigma _{D}=\sigma _{S_{2}}$
  
  $\varepsilon _{m}=\varepsilon _{D}+\varepsilon _{S_{2}}$
  
  where the subscripts $m$, $D$, $S_{1}$ and $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:
  
      ${\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}}}$
  
  The equation can also be expressed as:
  
      $\sigma (t)+{\frac {\eta }{E_{2}}}{\frac {d\sigma (t)}{dt}}=E_{1}\varepsilon (t)+{\frac {\eta (E_{1}+E_{2})}{E_{2}}}{\frac {d\varepsilon (t)}{dt}}$
  
  or, in dot notation:
  
      $\sigma +{\frac {\eta }{E_{2}}}{\dot {\sigma }}=E_{1}\varepsilon +{\frac {\eta (E_{1}+E_{2})}{E_{2}}}{\dot {\varepsilon }}$
  
  The relaxation time, $\tau$, is different for each material and is equal to
  
      $\tau ={\frac {\eta }{E_{2}}}$
  
  ### Kelvin-Voigt representation
  
  This model consists of two systems in series. The first, referred to as the Kelvin arm, contains a spring ($E=E_{2}$) and dashpot (viscosity $\eta$) in parallel. The other system contains only a spring ($E=E_{1}$).
  
  These relationships help relate the various stresses and strains in the overall system and the Kelvin arm:
  
  $\sigma _{tot}=\sigma _{k}=\sigma _{S_{1}}$
  
  $\varepsilon _{tot}=\varepsilon _{k}+\varepsilon _{S_{1}}$
  
  $\sigma _{k}=\sigma _{D}+\sigma _{S_{2}}$
  
  $\varepsilon _{k}=\varepsilon _{D}=\varepsilon _{S_{2}}$
  
  where the subscripts $k$, $D$, $S_{1}$,and $S_{2}$ refer to Kelvin, 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:
  
      $\sigma (t)+{\frac {\eta }{E_{1}+E_{2}}}{\frac {d\sigma (t)}{dt}}={\frac {E_{1}E_{2}}{E_{1}+E_{2}}}\varepsilon (t)+{\frac {E_{1}\eta }{E_{1}+E_{2}}}{\frac {d\varepsilon (t)}{dt}}$
  
  or, in dot notation:
  
      $\sigma +{\frac {\eta }{E_{1}+E_{2}}}{\dot {\sigma }}={\frac {E_{1}E_{2}}{E_{1}+E_{2}}}\varepsilon +{\frac {E_{1}\eta }{E_{1}+E_{2}}}{\dot {\varepsilon }}$
  
  The retardation time, ${\bar {\tau }}$, is different for each material and is equal to
  
      ${\bar {\tau }}={\frac {\eta }{E_{2}}}$
  
  ## 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.
  
  The fluid model equivalent to the standard linear solid model includes a dashpot in series with the Kelvin–Voigt model and is called the Jeffreys model.
  
  