Add vanishing strain wrapper (#44)

ikarus-project · Sep 27, 2024 · eaea696 · eaea696
1 parent e531151
commit eaea696
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 1 deletion.
diff --git a/.gitignore b/.gitignore
@@ -15,3 +15,6 @@ cmake-build-release-dockerfull/
 cmake-build-relwithdebinfo/
 docs/website/drawio-exporter
 *.pyc
+
+.devcontainer/* 
+docker-compose.yml
diff --git a/src/iks008_cooksMembrane.cpp b/src/iks008_cooksMembrane.cpp
@@ -24,6 +24,7 @@
 #include <ikarus/finiteelements/mechanics/enhancedassumedstrains.hh>
 #include <ikarus/finiteelements/mechanics/linearelastic.hh>
 #include <ikarus/finiteelements/mechanics/loads.hh>
+#include <ikarus/finiteelements/mechanics/materials/linearelasticity.hh>
 #include <ikarus/solver/linearsolver/linearsolver.hh>
 #include <ikarus/utils/basis.hh>
 #include <ikarus/utils/dirichletvalues.hh>
@@ -134,7 +135,9 @@ int main(int argc, char** argv) {
       }
 
       BoundaryPatch<decltype(gridView)> neumannBoundary(gridView, neumannVertices);
-      auto sk      = skills(linearElastic({E, nu}), eas(numberOfEASParameters), volumeLoad<2>(vL),
+
+      auto linMat  = Ikarus::LinearElasticity(Ikarus::toLamesFirstParameterAndShearModulus({.emodul = E, .nu = nu}));
+      auto sk      = skills(linearElastic(planeStress(linMat)), eas(numberOfEASParameters), volumeLoad<2>(vL),
                             neumannBoundaryLoad(&neumannBoundary, neumannBl));
       using FEType = decltype(makeFE(basis, sk));
       std::vector<FEType> fes;