Linear Maxwell viscoelasticity

This problem considers a squared medium subject to external velocities leading to a stress build-up. The material deforms following a linear Maxwell viscoelastic constitutive model. This problem was first introduced by Gerya and Yuen (2007).

Setup

The squared medium is subject to a compressive velocity horizontally and tension velocity vertically resulting in a deviatoric stress build-up. The setup is sketched in Figure 1.

Figure 1: Setup for the linear Maxwell viscoelastic medium.

Solutions

The stress evolution is governed by the following constitutive model:

var element = document.getElementById("moose-equation-b77155bb-7cc2-449a-b270-c33e7563516a");katex.render("\\dot{\\tau}_{ij} = 2 G \\left(\\dot{e}_{ij} - \\frac{\\tau_{ij}}{2 \\eta}\\right)", element, {displayMode:true,throwOnError:false});

where $var element = document.getElementById("moose-equation-f56af10a-c981-4142-90a1-1cc9c2a12253");katex.render("\\tau_{ij}", element, {displayMode:false,throwOnError:false});$ is the deviatoric stress tensor and $var element = document.getElementById("moose-equation-c703e6a5-7da0-41c6-8bc0-b881c9c42894");katex.render("e_{ij}", element, {displayMode:false,throwOnError:false});$ the deviatoric strain tensor.

The solution for this problem is given by Gerya and Yuen (2007). The stress solution is given as:

var element = document.getElementById("moose-equation-2c77834e-d8e1-4504-8f8e-4a3f7c9c2f93");katex.render("\\tau = \\eta \\dot{\\gamma} \\left(1 - \\exp\\left(- \\frac{G t}{\\eta}\\right)\\right)", element, {displayMode:true,throwOnError:false});

where:

var element = document.getElementById("moose-equation-5dac889e-2f18-4653-b843-cee4139b63ea");katex.render("\\sigma_{xx} = - \\tau, \\quad \\sigma_{yy} = \\tau.", element, {displayMode:true,throwOnError:false});

var element = document.getElementById("moose-equation-56ae84f4-b608-45e6-977c-d597a462e0fa");katex.render("\\dot{\\gamma} = \\frac{v}{h}", element, {displayMode:true,throwOnError:false});

With the following dimensionless quantities: $var element = document.getElementById("moose-equation-ce45ad60-098a-4a22-af40-6c3c2a450582");katex.render("\\overline{t} = \\frac{G}{\\eta} t", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-a1751584-8b54-42a3-b668-28b4b1a2c75c");katex.render("\\overline{\\tau} = \\frac{\\tau}{\\eta \\dot{\\gamma}}", element, {displayMode:false,throwOnError:false});$ , the solution can be expressed as:

(1)var element = document.getElementById("moose-equation-ff55ace8-8822-4cf0-bc24-ff856c253557");katex.render("\\overline{\\tau} = 1 - \\exp\\left(- t\\right)", element, {displayMode:true,throwOnError:false});

Figure 2: Stress evolution in a linear viscoelastic Maxwell medium.

Complete Source Files

linear_maxwell.i

References

Taras V. Gerya and David A. Yuen. Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems. Physics of the Earth and Planetary Interiors, 163(1-4):83–105, 2007. doi:10.1016/j.pepi.2007.04.015.

@article{Gerya2007,
    author = "Gerya, Taras V. and Yuen, David A.",
    year = "2007",
    title = "{Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems}",
    journal = "Physics of the Earth and Planetary Interiors",
    issn = "0031-9201",
    doi = "10.1016/j.pepi.2007.04.015",
    pages = "83--105",
    number = "1-4",
    volume = "163"
}

Linear Kelvin visc. . .

Setup
Solutions
Complete Source Files
References