Skip to content
Snippets Groups Projects
Commit d36286e5 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Re-arrange boundary parameters

parent 2d9eaad2
No related branches found
No related tags found
No related merge requests found
...@@ -220,12 +220,12 @@ int main(int argc, char *argv[]) { ...@@ -220,12 +220,12 @@ int main(int argc, char *argv[]) {
std::vector<double> coefficientOfFriction; std::vector<double> coefficientOfFriction;
coefficientOfFriction.resize(grid.size(grid.maxLevel(), dim)); coefficientOfFriction.resize(grid.size(grid.maxLevel(), dim));
std::fill(coefficientOfFriction.begin(), coefficientOfFriction.end(), std::fill(coefficientOfFriction.begin(), coefficientOfFriction.end(),
parset.get<double>("boundary.mu")); parset.get<double>("boundary.friction.mu"));
std::vector<double> normalStress; std::vector<double> normalStress;
normalStress.resize(grid.size(grid.maxLevel(), dim)); normalStress.resize(grid.size(grid.maxLevel(), dim));
std::fill(normalStress.begin(), normalStress.end(), std::fill(normalStress.begin(), normalStress.end(),
parset.get<double>("boundary.normalstress")); parset.get<double>("boundary.friction.normalstress"));
Dune::GlobalNonlinearity<dim> *myGlobalNonlinearity; Dune::GlobalNonlinearity<dim> *myGlobalNonlinearity;
...@@ -234,11 +234,13 @@ int main(int argc, char *argv[]) { ...@@ -234,11 +234,13 @@ int main(int argc, char *argv[]) {
if (friction_model == std::string("Ruina")) { if (friction_model == std::string("Ruina")) {
std::vector<double> a; std::vector<double> a;
a.resize(grid.size(grid.maxLevel(), dim)); a.resize(grid.size(grid.maxLevel(), dim));
std::fill(a.begin(), a.end(), parset.get<double>("boundary.a")); std::fill(a.begin(), a.end(),
parset.get<double>("boundary.friction.ruina.a"));
std::vector<double> eta; std::vector<double> eta;
eta.resize(grid.size(grid.maxLevel(), dim)); eta.resize(grid.size(grid.maxLevel(), dim));
std::fill(eta.begin(), eta.end(), parset.get<double>("boundary.eta")); std::fill(eta.begin(), eta.end(),
parset.get<double>("boundary.friction.eta"));
myGlobalNonlinearity = new Dune::GlobalRuinaNonlinearity<dim>( myGlobalNonlinearity = new Dune::GlobalRuinaNonlinearity<dim>(
nodalIntegrals, a, coefficientOfFriction, eta, normalStress); nodalIntegrals, a, coefficientOfFriction, eta, normalStress);
......
...@@ -11,11 +11,11 @@ nu = 0.3 ...@@ -11,11 +11,11 @@ nu = 0.3
maxiterations = 100000 maxiterations = 100000
tolerance = 1e-6 tolerance = 1e-6
[boundary] [boundary.friction]
a = 0.0015
normalstress = 0.1 normalstress = 0.1
mu = 0.75 mu = 0.75
eta = 1 eta = 1
[boundary.friction]
model = Laursen model = Laursen
[boundary.friction.ruina]
a = 0.0015
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment