mueslimaterials.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2021-2025 The Ikarus Developers mueller@ibb.uni-stuttgart.de
// SPDX-License-Identifier: LGPL-3.0-or-later
 
#pragma once
 
#if ENABLE_MUESLI
 
  #include <ikarus/finiteelements/mechanics/materials/muesli/mueslifinite.hh>
  #include <ikarus/finiteelements/mechanics/materials/muesli/mueslihelpers.hh>
  #include <ikarus/finiteelements/mechanics/materials/muesli/mueslismall.hh>
 
namespace Ikarus::Materials {
 
template <Concepts::MPTuple MPT>
inline auto makeMuesliLinearElasticity(const MPT& mpt) {
  auto muesliParameters = propertiesFromIkarusMaterialParameters(mpt);
  return SmallStrain<muesli::elasticIsotropicMaterial>(muesliParameters);
}
 
template <Concepts::MPTuple MPT>
inline auto makeMuesliNeoHooke(const MPT& mpt, bool useDeviatoricStretches = false) {
  auto muesliParameters = propertiesFromIkarusMaterialParameters(mpt);
  if (useDeviatoricStretches)
    addTag(muesliParameters, "subtype regularized");
  return FiniteStrain<muesli::neohookeanMaterial>(muesliParameters);
}
 
template <Concepts::MPTuple MPT>
inline auto makeMuesliSVK(const MPT& mpt) {
  auto muesliParameters = propertiesFromIkarusMaterialParameters(mpt);
  return FiniteStrain<muesli::svkMaterial>(muesliParameters);
}
 
inline auto makeMuesliArrudaBoyce(double C1, double lambda_m, double K, bool compressible = true) {
  auto muesliParameters = muesli::materialProperties{};
  muesliParameters.insert({"c1", C1});
  muesliParameters.insert({"lambdam", lambda_m});
  muesliParameters.insert({"bulk", K});
  if (compressible)
    addTag(muesliParameters, "compressible");
  return FiniteStrain<muesli::arrudaboyceMaterial>(muesliParameters);
}
 
inline auto makeMuesliYeoh(std::array<double, 3> C, double K, bool compressible = true) {
  auto muesliParameters = muesli::materialProperties{};
  muesliParameters.insert({"c1", C[0]});
  muesliParameters.insert({"c2", C[1]});
  muesliParameters.insert({"c3", C[2]});
  muesliParameters.insert({"bulk", K});
  if (compressible)
    addTag(muesliParameters, "compressible");
  return FiniteStrain<muesli::yeohMaterial>(muesliParameters);
}
 
inline auto makeMooneyRivlin(std::array<double, 3> alpha, bool incompressible = false) {
  auto muesliParameters = muesli::materialProperties{};
  muesliParameters.insert({"alpha0", alpha[0]});
  muesliParameters.insert({"alpha1", alpha[1]});
  muesliParameters.insert({"alpha2", alpha[2]});
  if (incompressible)
    addTag(muesliParameters, "incompressible");
  return FiniteStrain<muesli::mooneyMaterial>(muesliParameters);
}
 
} // namespace Ikarus::Materials
 
#else
  #error Muesli materials depends on the Muesli library, which is not included
#endif