Skip to content
Snippets Groups Projects
Commit 6e057788 authored by Dmitri Naumov's avatar Dmitri Naumov
Browse files

[PL] Add ViscoElasticBurgers model to SD Process.

parent 4265da68
No related branches found
No related tags found
No related merge requests found
Visco elastic material model based on Burgers' rheological model. Viscosities
Visco elastic material model based on Burgers rheological model. Viscosities
and Kelvin shear modulus are stress dependent according to LUBBY2
model. Heusermann 1983.
model. See \cite Heusermann1983 and \cite Heusermann2003 for details.
......@@ -54,3 +54,27 @@
Month = may,
Url = {http://arxiv.org/abs/1605.09606}
}
@Book{Heusermann1983,
Title = {{Entwicklung mathematisch-mechanischer Modelle zur Beschreibung des Stoffverhaltens von Salzgestein in Abh{\"a}ngigkeit von der Zeit und der Temperatur auf der Grundlage von Laborversuchen mit begleitenden kontinuumsmechanischen Berechnungen nach der Methode der finiten Elemente}},
Author = {Heusermann, Stefan and Lux, Karl-Heinz and Rokahr, Reinhard},
Publisher = {Fachinformationszentrum Energie, Physik, Mathematik Karlsruhe},
Year = {1983}
}
@Article{Heusermann2003,
Title = {Nonlinear finite-element analysis of solution mined storage caverns in rock salt using the {LUBBY2} constitutive model },
Author = {Stefan Heusermann and Olaf Rolfs and Uwe Schmidt},
Journal = {Computers \& Structures },
Year = {2003},
Note = {K.J Bathe 60th Anniversary Issue },
Number = {8–11},
Pages = {629 - 638},
Volume = {81},
Doi = {http://dx.doi.org/10.1016/S0045-7949(02)00415-7},
ISSN = {0045-7949},
Keywords = {Rock salt},
Owner = {nagelt},
Timestamp = {2015.06.01},
Url = {http://www.sciencedirect.com/science/article/pii/S0045794902004157}
}
......@@ -12,6 +12,7 @@
#include <cassert>
#include "MaterialLib/SolidModels/CreateLinearElasticIsotropic.h"
#include "MaterialLib/SolidModels/CreateLubby2.h"
#include "ProcessLib/Utils/ParseSecondaryVariables.h"
#include "SmallDeformationProcess.h"
......@@ -77,6 +78,11 @@ createSmallDeformationProcess(
MaterialLib::Solids::createLinearElasticIsotropic<DisplacementDim>(
parameters, constitutive_relation_config);
}
else if (type == "Lubby2")
{
material = MaterialLib::Solids::createLubby2<DisplacementDim>(
parameters, constitutive_relation_config);
}
else
{
OGS_FATAL(
......
......@@ -10,6 +10,7 @@
#ifndef PROCESS_LIB_CREATESMALLDEFORMATIONPROCESS_H_
#define PROCESS_LIB_CREATESMALLDEFORMATIONPROCESS_H_
#include "MaterialLib/SolidModels/CreateLubby2.h"
#include "ProcessLib/Process.h"
namespace ProcessLib
......
......@@ -14,6 +14,7 @@
#include <vector>
#include "MaterialLib/SolidModels/LinearElasticIsotropic.h"
#include "MaterialLib/SolidModels/Lubby2.h"
#include "MathLib/LinAlg/Eigen/EigenMapTools.h"
#include "NumLib/Extrapolation/ExtrapolatableElement.h"
#include "NumLib/Fem/FiniteElement/TemplateIsoparametric.h"
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment