diff --git a/src/assemblers.cc b/src/assemblers.cc index c2a9fae984b1aa9154ab6f17c001b36ed4346a45..b2fea9b0d992a8239483ee6dc924c99a0251f72c 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; + } } }