From dc3abdeb60ccbd83db1c2c8882fae80d34e444b8 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 22 Jul 2014 16:20:58 +0200
Subject: [PATCH] [Cleanup] Ignore positive normal stress
---
src/assemblers.cc | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/src/assemblers.cc b/src/assemblers.cc
index c2a9fae9..b2fea9b0 100644
--- a/src/assemblers.cc
+++ b/src/assemblers.cc
@@ -116,7 +116,11 @@ void MyAssembler<GridView, dimension>::assembleNormalStress(
frictionalBoundary.getNormals(normals);
for (size_t i = 0; i < traction.size(); ++i) {
normalStress[i] = normals[i] * traction[i];
- assert(normalStress[i] <= 0.0);
+ if (normalStress[i] > 0.0) {
+ normalStress[i] = 0.0;
+ std::cout << "Warning: Manually reducing positive normal stress to zero."
+ << std::endl;
+ }
}
}
--
GitLab