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

[Cleanup] Ignore positive normal stress

parent 9eac7877
No related branches found
No related tags found
No related merge requests found
......@@ -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;
}
}
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment