factory.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2021-2025 The Ikarus Developers ikarus@ibb.uni-stuttgart.de
// SPDX-License-Identifier: LGPL-3.0-or-later
 
#pragma once
 
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/deviatoric/arrudaboyce.hh>
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/deviatoric/blatzko.hh>
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/deviatoric/gent.hh>
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/deviatoric/interface.hh>
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/deviatoric/invariantbased.hh>
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/deviatoric/nodeviatoricfunction.hh>
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/deviatoric/ogden.hh>
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/interface.hh>
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/neohooke.hh>
#include <ikarus/finiteelements/mechanics/materials/hyperelastic/volumetric/volumetricfunctions.hh>
#include <ikarus/finiteelements/mechanics/materials/tags.hh>
 
namespace Ikarus::Materials {
 
inline auto makeBlatzKo(double mu) {
  auto bk  = BlatzKo(mu);
  auto dev = Deviatoric<decltype(bk)>(bk);
 
  return Hyperelastic(dev);
}
 
template <int n, PrincipalStretchTags tag, Concepts::VolumetricFunction VolumetricFunction = VF0>
inline auto makeOgden(const typename Ogden<n, tag>::MaterialParameters& mu,
                      const typename Ogden<n, tag>::MaterialExponents alpha, double K = 0.0,
                      const VolumetricFunction& vf = VolumetricFunction{}) {
  auto ogPre = Ogden<n, tag>(mu, alpha);
  auto dev   = Deviatoric(ogPre);
  auto vol   = Volumetric(K, vf);
 
  return Hyperelastic(dev, vol);
}
template <int n, Concepts::VolumetricFunction VolumetricFunction = VF0>
inline auto makeInvariantBased(const typename InvariantBased<n>::MaterialParameters& mu,
                               const typename InvariantBased<n>::Exponents pex,
                               const typename InvariantBased<n>::Exponents qex, double K = 0.0,
                               const VolumetricFunction& vf = VolumetricFunction{}) {
  auto invariantBasedPre = InvariantBased<n>(pex, qex, mu);
  auto dev               = Deviatoric(invariantBasedPre);
  auto vol               = Volumetric(K, vf);
 
  return Hyperelastic(dev, vol);
}
 
template <Concepts::VolumetricFunction VolumetricFunction = VF0>
inline auto makeMooneyRivlin(const typename InvariantBased<2>::MaterialParameters& mu, double K = 0.0,
                             const VolumetricFunction& vf = VolumetricFunction{}) {
  typename InvariantBased<2>::Exponents pex = {1, 0};
  typename InvariantBased<2>::Exponents qex = {0, 1};
 
  return makeInvariantBased<2, VolumetricFunction>(mu, pex, qex, K, vf);
}
 
template <Concepts::VolumetricFunction VolumetricFunction = VF0>
inline auto makeYeoh(const typename InvariantBased<3>::MaterialParameters& mu, double K = 0.0,
                     const VolumetricFunction& vf = VolumetricFunction{}) {
  typename InvariantBased<3>::Exponents pex = {1, 2, 3};
  typename InvariantBased<3>::Exponents qex = {0, 0, 0};
 
  return makeInvariantBased<3, VolumetricFunction>(mu, pex, qex, K, vf);
}
 
template <Concepts::VolumetricFunction VolumetricFunction = VF0>
inline auto makeArrudaBoyce(const ArrudaBoyceMatParameters& matPar, double K = 0.0,
                            const VolumetricFunction& vf = VolumetricFunction{}) {
  auto abPre = ArrudaBoyce(matPar);
  auto dev   = Deviatoric(abPre);
  auto vol   = Volumetric(K, vf);
 
  return Hyperelastic(dev, vol);
}
 
template <Concepts::VolumetricFunction VolumetricFunction = VF0>
inline auto makeGent(const GentMatParameters& matPar, double K = 0.0,
                     const VolumetricFunction& vf = VolumetricFunction{}) {
  auto gentPre = Gent(matPar);
  auto dev     = Deviatoric(gentPre);
  auto vol     = Volumetric(K, vf);
 
  return Hyperelastic(dev, vol);
}
 
auto makePureVolumetric(const Concepts::VolumetricFunction auto& vf, double K) {
  auto noDev = Deviatoric(NoDev{});
  auto vol   = Volumetric(K, vf);
  return Hyperelastic(noDev, vol);
}
 
} // namespace Ikarus::Materials