From bd8a598e99df305ecff74b1a8e4580311efe7ab1 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 28 Sep 2007 11:56:38 +0000
Subject: [PATCH] add stress output
[[Imported from SVN: r1647]]
---
dirneucoupling.cc | 21 ++++++++++++++++-----
1 file changed, 16 insertions(+), 5 deletions(-)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 021deb8b..2176c806 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -25,6 +25,7 @@
#include "../common/prolongboundarypatch.hh"
#include "../common/sampleonbitfield.hh"
#include "../common/neumannassembler.hh"
+#include "../common/computestress.hh"
#include "src/quaternion.hh"
#include "src/rodassembler.hh"
@@ -37,6 +38,9 @@
// Space dimension
const int dim = 3;
+const double E = 2.5e5;
+const double nu = 0.3;
+
using namespace Dune;
using std::string;
@@ -82,10 +86,10 @@ int main (int argc, char *argv[]) try
const double initialTrustRegionRadius = parameterSet.get<double>("initialTrustRegionRadius");
const double dirichletDamping = parameterSet.get<double>("dirichletDamping");
FieldVector<double,3> neumannDamping;
- neumannDamping = parameterSet.get<double>("neumannDamping");
- neumannDamping[0] = parameterSet.get<double>("neumannDampingT");
- neumannDamping[1] = parameterSet.get<double>("neumannDampingT");
- neumannDamping[2] = parameterSet.get<double>("neumannDampingL");
+ neumannDamping = parameterSet.get("neumannDamping", double(1));
+ neumannDamping[0] = parameterSet.get("neumannDampingT", neumannDamping[0]);
+ neumannDamping[1] = parameterSet.get("neumannDampingT", neumannDamping[1]);
+ neumannDamping[2] = parameterSet.get("neumannDampingL", neumannDamping[2]);
std::string resultPath = parameterSet.get("resultPath", "");
// Problem settings
@@ -193,7 +197,7 @@ int main (int argc, char *argv[]) try
// Assemble 3d linear elasticity problem
// //////////////////////////////////////////
LeafP1Function<GridType,double,dim> u(grid),f(grid);
- LinearElasticityLocalStiffness<GridType,double> lstiff(2.5e5, 0.3);
+ LinearElasticityLocalStiffness<GridType,double> lstiff(E, nu);
LeafP1OperatorAssembler<GridType,double,dim> hessian3d(grid);
hessian3d.assemble(lstiff,u,f);
@@ -626,8 +630,15 @@ int main (int argc, char *argv[]) try
// //////////////////////////////
LeafAmiraMeshWriter<GridType> amiraMeshWriter(grid);
amiraMeshWriter.addVertexData(x3d, grid.leafIndexSet());
+
+ BlockVector<FieldVector<double,1> > stress;
+ Stress<GridType,dim>::getStress(grid, x3d, stress, E, nu);
+ amiraMeshWriter.addVertexData(stress, grid.leafIndexSet());
+
amiraMeshWriter.write(resultPath + "grid.result");
+
+
writeRod(rodX, resultPath + "rod3d.result");
} catch (Exception e) {
--
GitLab