CreatePermeabilityMohrCoulombFailureIndexModel.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "CreatePermeabilityMohrCoulombFailureIndexModel.h"
 
#include <string>
 
#include "BaseLib/ConfigTree.h"
#include "Parameter.h"
#include "ParameterLib/CoordinateSystem.h"
#include "ParameterLib/Parameter.h"
#include "ParameterLib/Utils.h"
#include "PermeabilityMohrCoulombFailureIndexModel.h"
 
namespace MaterialPropertyLib
{
std::unique_ptr<Property> createPermeabilityMohrCoulombFailureIndexModel(
    int const geometry_dimension,
    BaseLib::ConfigTree const& config,
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    ParameterLib::CoordinateSystem const* const local_coordinate_system)
{
    if ((geometry_dimension != 2) && (geometry_dimension != 3))
    {
        OGS_FATAL(
            "The PermeabilityMohrCoulombFailureIndexModel is implemented only "
            "for 2D or 3D problems");
    }

    config.checkConfigParameter("type",
                                "PermeabilityMohrCoulombFailureIndexModel");
 
    // Second access for storage.
    auto property_name = config.peekConfigParameter<std::string>("name");
 
    DBUG("Create PermeabilityMohrCoulombFailureIndexModel property {:s}.",
         property_name);
 
    std::string const& parameter_name =
        config.getConfigParameter<std::string>("initial_permeability");
    auto const& parameter_k0 = ParameterLib::findParameter<double>(
        parameter_name, parameters, 0, nullptr);
 
    auto const kr =
        config.getConfigParameter<double>("reference_permeability");
    auto const b =
        config.getConfigParameter<double>("fitting_factor");
    auto const c =
        config.getConfigParameter<double>("cohesion");
    auto const phi =
        config.getConfigParameter<double>("friction_angle");
    auto const max_k =
        config.getConfigParameter<double>("maximum_permeability");
    auto const t_sigma_max =
        config.getConfigParameter<double>("tensile_strength_parameter");
 
    if (geometry_dimension == 2)
    {
        return std::make_unique<PermeabilityMohrCoulombFailureIndexModel<2>>(
            std::move(property_name), parameter_k0, kr, b, c, phi, max_k,
            t_sigma_max, local_coordinate_system);
    }
 
    return std::make_unique<PermeabilityMohrCoulombFailureIndexModel<3>>(
        std::move(property_name), parameter_k0, kr, b, c, phi, max_k,
        t_sigma_max, local_coordinate_system);
}
}  // namespace MaterialPropertyLib