Skip to content

Commit

Permalink
fix for materials namespace
Browse files Browse the repository at this point in the history
  • Loading branch information
henrij22 committed Nov 12, 2024
1 parent eaea696 commit 312ce21
Show file tree
Hide file tree
Showing 2 changed files with 6 additions and 7 deletions.
6 changes: 2 additions & 4 deletions src/iks006_nonlinear2DSolid.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -15,8 +15,6 @@
#include <dune/iga/nurbsbasis.hh>
#include <dune/iga/nurbsgrid.hh>

#include "spdlog/spdlog.h"

#include <Eigen/Core>

#include <ikarus/assembler/simpleassemblers.hh>
Expand Down Expand Up @@ -151,8 +149,8 @@ auto run() {

auto matParameter = Ikarus::toLamesFirstParameterAndShearModulus({.emodul = 1000, .nu = 0.3});

Ikarus::StVenantKirchhoff matSVK(matParameter);
auto reducedMat = planeStress(matSVK, 1e-8);
Materials::StVenantKirchhoff matSVK(matParameter);
auto reducedMat = Materials::planeStress(matSVK, 1e-8);

auto vL = [](auto& globalCoord, auto& lamb) {
Eigen::Vector2d fext;
Expand Down
7 changes: 4 additions & 3 deletions src/iks008_cooksMembrane.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -25,6 +25,7 @@
#include <ikarus/finiteelements/mechanics/linearelastic.hh>
#include <ikarus/finiteelements/mechanics/loads.hh>
#include <ikarus/finiteelements/mechanics/materials/linearelasticity.hh>
#include <ikarus/finiteelements/mechanics/materials/vanishingstress.hh>
#include <ikarus/solver/linearsolver/linearsolver.hh>
#include <ikarus/utils/basis.hh>
#include <ikarus/utils/dirichletvalues.hh>
Expand Down Expand Up @@ -136,9 +137,9 @@ int main(int argc, char** argv) {

BoundaryPatch<decltype(gridView)> neumannBoundary(gridView, neumannVertices);

auto linMat = Ikarus::LinearElasticity(Ikarus::toLamesFirstParameterAndShearModulus({.emodul = E, .nu = nu}));
auto sk = skills(linearElastic(planeStress(linMat)), eas(numberOfEASParameters), volumeLoad<2>(vL),
neumannBoundaryLoad(&neumannBoundary, neumannBl));
auto linMat = Materials::LinearElasticity(Ikarus::toLamesFirstParameterAndShearModulus({.emodul = E, .nu = nu}));
auto sk = skills(linearElastic(Materials::planeStress(linMat)), eas(numberOfEASParameters), Ikarus::volumeLoad<2>(vL),
neumannBoundaryLoad(&neumannBoundary, neumannBl));
using FEType = decltype(makeFE(basis, sk));
std::vector<FEType> fes;
for (auto&& ge : elements(gridView)) {
Expand Down

0 comments on commit 312ce21

Please sign in to comment.