CreateLinearElasticIsotropic.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "CreateLinearElasticIsotropic.h"
 
#include "LinearElasticIsotropic.h"
#include "ParameterLib/Utils.h"
 
namespace MaterialLib
{
namespace Fracture
{
template <int DisplacementDim>
std::unique_ptr<FractureModelBase<DisplacementDim>>
createLinearElasticIsotropic(
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    BaseLib::ConfigTree const& config)
{
    config.checkConfigParameter("type", "LinearElasticIsotropic");
    DBUG("Create LinearElasticIsotropic material");
 
    auto const& Kn = ParameterLib::findParameter<double>(
        config, "normal_stiffness", parameters, 1);
 
    auto const& Ks = ParameterLib::findParameter<double>(
        config, "shear_stiffness", parameters, 1);
 
    auto const penalty_aperture_cutoff =
        config.getConfigParameter<double>("penalty_aperture_cutoff");
 
    auto const tension_cutoff =
        config.getConfigParameter<bool>("tension_cutoff");
 
    typename LinearElasticIsotropic<DisplacementDim>::MaterialProperties mp{Kn,
                                                                            Ks};
 
    return std::make_unique<LinearElasticIsotropic<DisplacementDim>>(
        penalty_aperture_cutoff, tension_cutoff, mp);
}
 
template std::unique_ptr<FractureModelBase<2>> createLinearElasticIsotropic(
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    BaseLib::ConfigTree const& config);
 
template std::unique_ptr<FractureModelBase<3>> createLinearElasticIsotropic(
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    BaseLib::ConfigTree const& config);
 
}  // namespace Fracture
}  // namespace MaterialLib