diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index fe3c4590fa63f39e9871147503420c1976c1a760..36827736004729f89784a3129f477e6e7fa84c2f 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -51,8 +51,6 @@ #include <dune/fufem/assemblers/localassemblers/stvenantkirchhoffassembler.hh> #include <dune/fufem/assemblers/localassemblers/vonmisesstressassembler.hh> #include <dune/fufem/assemblers/operatorassembler.hh> -#include <dune/fufem/boundarypatch.hh> -#include <dune/fufem/computestress.hh> #include <dune/fufem/dunepython.hh> #include <dune/fufem/functions/basisgridfunction.hh> #include <dune/fufem/functions/constantfunction.hh> @@ -272,7 +270,6 @@ int main(int argc, char *argv[]) { auto const nodalIntegrals = assemble_frictional<GridType, GridView, SmallVector, P1Basis>( leafView, p1Basis, frictionalNodes); - BoundaryPatch<GridView> const frictionalBoundary(leafView, frictionalNodes); // {{{ Initialise vectors VectorType u(finestSize); @@ -399,21 +396,6 @@ int main(int argc, char *argv[]) { timeSteppingScheme->extractDisplacement(u); timeSteppingScheme->extractVelocity(ud); - // Update the normal stress - if (parset.get<bool>("updateNormalStress")) { - VectorType tmp(finestSize); - Stress<GridView>::getElasticSurfaceNormalStress(frictionalBoundary, - u, tmp, E, nu); - for (size_t i = 0; i < frictionalNodes.size(); ++i) - if (frictionalNodes[i][0]) { - if (tmp[i][1] < 0) { - DUNE_THROW(Dune::Exception, - "Encountered negative normal stress"); - } else - surfaceNormalStress[i] = tmp[i][1]; - } - } - // Update the state for (size_t i = 0; i < frictionalNodes.size(); ++i) { if (frictionalNodes[i][0]) { diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset index 3c90028836e6cff9a2401148d8e5961f10edfec2..bf591145a08d15dc7b89ad58eba32f120022cdfc 100644 --- a/src/one-body-sample.parset +++ b/src/one-body-sample.parset @@ -11,8 +11,6 @@ enableTimer = false timeSteppingScheme = implicitEuler -updateNormalStress = false - [grid] refinements = 4