mueslihelpers.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 <muesli/muesli.h>
 
  #include <ikarus/finiteelements/physicshelper.hh>
  #include <ikarus/utils/tensorutils.hh>
 
namespace Ikarus::Concepts {
template <typename MAT>
concept MuesliMaterialImplementation = std::is_base_of_v<muesli::material, MAT>;
} // namespace Ikarus::Concepts
 
namespace Ikarus::Materials {
 
template <Concepts::MPTuple MPT>
inline muesli::materialProperties propertiesFromIkarusMaterialParameters(const MPT& mpt) {
  auto converter = convertLameConstants(mpt);
 
  auto mpm = muesli::materialProperties{};
  mpm.insert({"lambda", converter.toLamesFirstParameter()});
  mpm.insert({"mu", converter.toShearModulus()});
 
  return mpm;
}
 
inline void addTag(muesli::materialProperties& mpm, const std::string& tagName, double tagValue = 0) {
  mpm.insert({tagName, tagValue});
}
 
inline istensor toistensor(const Eigen::Matrix<double, 3, 3>& C) {
  return istensor(C(0, 0), C(1, 1), C(2, 2), C(1, 2), C(2, 0), C(0, 1));
}
 
inline itensor toitensor(const Eigen::Matrix<double, 3, 3>& C) {
  return itensor(C(0, 0), C(0, 1), C(0, 2), C(1, 0), C(1, 1), C(1, 2), C(2, 0), C(2, 1), C(2, 2));
}
 
inline auto toEigenMatrix(const istensor& it) {
  auto S = Eigen::Matrix<double, 3, 3>{};
  for (auto i : Dune::range(3))
    for (auto j : Dune::range(3))
      S(i, j) = it(i, j);
  return S;
}
 
inline auto toEigenTensor(const itensor4& it) {
  Eigen::TensorFixedSize<double, Eigen::Sizes<3, 3, 3, 3>> moduli{};
  for (auto i : Dune::range(3))
    for (auto j : Dune::range(3))
      for (auto k : Dune::range(3))
        for (auto l : Dune::range(3))
          moduli(i, j, k, l) = it(i, j, k, l);
 
  return moduli;
}
 
template <Concepts::MuesliMaterialImplementation MAT>
constexpr std::string materialName() {
  std::string matName = Dune::className<MAT>();
 
  // Find and remove the "muesli::" namespace
  std::string prefix = "muesli::";
  matName            = matName.substr(prefix.size());
 
  // Capitalize the first letter
  if (!matName.empty() && std::islower(matName[0])) {
    matName[0] = std::toupper(matName[0]);
  }
  return matName;
}
 
} // namespace Ikarus::Materials
 
#else
  #error Muesli materials depends on the Muesli library, which is not included
#endif