decomposehyperelastic.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.hh>
#include <ikarus/utils/concepts.hh>
 
namespace Ikarus {
namespace Materials {
 
  template <typename MAT>
  struct DecomposedMaterialTypes
  {
    using TupleType = decltype(decomposeHyperelasticAndGetMaterialParameters(std::declval<MAT>()));
    using DevType   = std::tuple_element_t<0, TupleType>; 
    using VolType   = std::tuple_element_t<1, TupleType>; 
    using DevParams = std::tuple_element_t<2, TupleType>; 
    using VolParams = std::tuple_element_t<3, TupleType>; 
  };
 
  template <Concepts::Material MAT>
  requires(MAT::isHyperelastic)
  auto decomposeHyperelasticAndGetMaterialParameters(const MAT& mat) {
    static_assert(
        Materials::UnderlyingMaterial_t<MAT>::hasVolumetricPart,
        "decomposeHyperelastic is only implemented for the case where both deviatoric and volumetric part exists.");
    const auto& umat                  = Materials::underlyingMaterial(mat);
    const auto dev                    = umat.deviatoricFunction();
    const auto vol                    = umat.volumetricFunction().volumetricFunction();
    const auto [devParams, volParams] = umat.materialParametersImpl();
    const auto hdev                   = Materials::Hyperelastic(dev);
    const auto hvol                   = Materials::makePureVolumetric(vol, 1.0); // material parameter set to unity.
    if constexpr (not MAT::isReduced) {
      return std::make_tuple(hdev, hvol, devParams, volParams);
    } else {
      constexpr auto fixedPairs = MAT::fixedPairs;
      if constexpr (MAT::isStrainVanished) {
        const auto reducedDev = Materials::VanishingStrain<fixedPairs, std::remove_cvref_t<decltype(hdev)>>(hdev);
        const auto reducedVol = Materials::VanishingStrain<fixedPairs, std::remove_cvref_t<decltype(hvol)>>(hvol);
        return std::make_tuple(reducedDev, reducedVol, devParams, volParams);
      } else {
        const auto reducedDev = Materials::VanishingStress<fixedPairs, std::remove_cvref_t<decltype(hdev)>>(hdev);
        const auto reducedVol = Materials::VanishingStress<fixedPairs, std::remove_cvref_t<decltype(hvol)>>(hvol);
        return std::make_tuple(reducedDev, reducedVol, devParams, volParams);
      }
    }
  }
} // namespace Materials
} // namespace Ikarus