CreateCoulomb.cpp

Go to the documentation of this file.

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