CreateLinearElasticTransverseIsotropic.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "CreateLinearElasticTransverseIsotropic.h"
 
#include "ParameterLib/Utils.h"
 
namespace MaterialLib
{
namespace Solids
{
template <int DisplacementDim>
std::unique_ptr<LinearElasticTransverseIsotropic<DisplacementDim>>
createLinearElasticTransverseIsotropic(
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    std::optional<ParameterLib::CoordinateSystem> const&
        local_coordinate_system,
    BaseLib::ConfigTree const& config, const bool skip_type_checking)
{
    if (!skip_type_checking)
    {
        config.checkConfigParameter("type", "LinearElasticTransverseIsotropic");
        DBUG("Create LinearElasticTransverseIsotropic material");
    }
 
    auto const& E_i = ParameterLib::findParameter<double>(
        config, "youngs_modulus_i", parameters, 1);
    DBUG("Use '{}' as the in-plane Young’s modulus, E_i.", E_i.name);
 
    auto const& E_a = ParameterLib::findParameter<double>(
        config, "youngs_modulus_a", parameters, 1);
    DBUG(
        "Use '{}' as the Young’s modulus w.r.t. the direction of anisotropy, "
        "E_a.",
        E_a.name);
 
    auto const& nu_ii = ParameterLib::findParameter<double>(
        config, "poissons_ratio_ii", parameters, 1);
    DBUG("Use '{}' as the in-plane Poisson’s ratio, nu_ii.", nu_ii.name);
 
    auto const& nu_ia = ParameterLib::findParameter<double>(
        config, "poissons_ratio_ia", parameters, 1);
    DBUG(
        "Use '{}' as the Poisson ratio perpendicular to the plane of isotropy, "
        "due to strain in the plane of isotropy, nu_ia.",
        nu_ia.name);
 
    auto const& G_ia = ParameterLib::findParameter<double>(
        config, "shear_modulus_ia", parameters, 1);
    DBUG(
        "Use '{}' as the shear modulus between directions of isotropy and "
        "anisotropy , G_ia.",
        G_ia.name);
 
    return std::make_unique<LinearElasticTransverseIsotropic<DisplacementDim>>(
        E_i, E_a, nu_ii, nu_ia, G_ia, local_coordinate_system);
}
 
template std::unique_ptr<LinearElasticTransverseIsotropic<2>>
createLinearElasticTransverseIsotropic(
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    std::optional<ParameterLib::CoordinateSystem> const&
        local_coordinate_system,
    BaseLib::ConfigTree const& config, const bool skip_type_checking);
 
template std::unique_ptr<LinearElasticTransverseIsotropic<3>>
createLinearElasticTransverseIsotropic(
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    std::optional<ParameterLib::CoordinateSystem> const&
        local_coordinate_system,
    BaseLib::ConfigTree const& config, const bool skip_type_checking);
}  // namespace Solids
}  // namespace MaterialLib