From 3bb03bf614280af75c54291b5ac75772549166ed Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 28 Jun 2013 13:52:57 +0000
Subject: [PATCH] boundary conditions for the L-shape problem
[[Imported from SVN: r9277]]
---
cosserat-continuum.cc | 11 ++++++++++-
1 file changed, 10 insertions(+), 1 deletion(-)
diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc
index e03b1c59..cff9b6c8 100644
--- a/cosserat-continuum.cc
+++ b/cosserat-continuum.cc
@@ -163,7 +163,7 @@ int main (int argc, char *argv[]) try
for (; vIt!=vEndIt; ++vIt) {
-#if 1 // Boundary conditions for the cantilever example
+#if 0 // Boundary conditions for the cantilever example
if (vIt->geometry().corner(0)[0] < 1.0+1e-3 /* or vIt->geometry().corner(0)[0] > upper[0]-1e-3*/ ) {
// Only translation dofs are Dirichlet
for (int j=0; j<3; j++)
@@ -178,6 +178,15 @@ int main (int argc, char *argv[]) try
for (int j=0; j<3; j++)
dirichletNodes[grid->leafIndexSet().index(*vIt)][j] = true;
}
+#endif
+#if 1 // Boundary conditions for the cantilever example
+ if (vIt->geometry().corner(0)[0] < 1.0) {
+ // Only translation dofs are Dirichlet
+ for (int j=0; j<3; j++)
+ dirichletNodes[grid->leafIndexSet().index(*vIt)][j] = true;
+ }
+ if (vIt->geometry().corner(0)[0] < -239 )
+ neumannNodes[grid->leafIndexSet().index(*vIt)][0] = true;
#endif
}
--
GitLab