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-f6ffd045-2ea1-4185-9259-3ffe17434efa");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-a79cc70a-7797-4760-90dd-4deee1b18527");katex.render("\\tau_{ij}", element, {displayMode:false,throwOnError:false});$ is the deviatoric stress tensor and $var element = document.getElementById("moose-equation-c5d352bf-28a7-4698-b555-c4e55ff2841e");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-b8191c1d-0645-41fa-bfb3-2da7d988ab5a");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-de773673-dce4-4a2a-96dd-4b1e44dbf613");katex.render("\\sigma_{xx} = - \\tau, \\quad \\sigma_{yy} = \\tau.", element, {displayMode:true,throwOnError:false});

var element = document.getElementById("moose-equation-7ac59807-eb56-4c00-821f-ccaa2e9bcc48");katex.render("\\dot{\\gamma} = \\frac{v}{h}", element, {displayMode:true,throwOnError:false});

With the following dimensionless quantities: $var element = document.getElementById("moose-equation-364038fe-bb08-4fb7-8d86-269f3cfda566");katex.render("\\overline{t} = \\frac{G}{\\eta} t", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-d0e55187-8672-4dfa-944a-aae12625ffdb");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-69db5999-90db-4819-b0be-ac7cd5c1629d");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