CreateLubby2.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "Lubby2.h"
#include "NumLib/CreateNewtonRaphsonSolverParameters.h"
#include "ParameterLib/Utils.h"
 
namespace MaterialLib
{
namespace Solids
{
namespace Lubby2
{
template <int DisplacementDim>
std::unique_ptr<Lubby2<DisplacementDim>> createLubby2(
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    BaseLib::ConfigTree const& config)
{
    config.checkConfigParameter("type", "Lubby2");
    DBUG("Create Lubby2 material");
 
    // Kelvin shear modulus.
    auto& kelvin_shear_modulus = ParameterLib::findParameter<double>(
        config, "kelvin_shear_modulus", parameters, 1);
 
    DBUG("Use '{:s}' as kelvin shear modulus parameter.",
         kelvin_shear_modulus.name);
 
    // Kelvin viscosity.
    auto& kelvin_viscosity = ParameterLib::findParameter<double>(
        config, "kelvin_viscosity", parameters, 1);
 
    DBUG("Use '{:s}' as kelvin viscosity parameter.", kelvin_viscosity.name);
 
    // Maxwell shear modulus.
    auto& maxwell_shear_modulus = ParameterLib::findParameter<double>(
        config, "maxwell_shear_modulus", parameters, 1);
 
    DBUG("Use '{:s}' as maxwell shear modulus parameter.",
         maxwell_shear_modulus.name);
 
    // Maxwell bulk modulus.
    auto& maxwell_bulk_modulus = ParameterLib::findParameter<double>(
        config, "maxwell_bulk_modulus", parameters, 1);
 
    DBUG("Use '{:s}' as maxwell bulk modulus parameter.",
         maxwell_bulk_modulus.name);
 
    // Maxwell viscosity.
    auto& maxwell_viscosity = ParameterLib::findParameter<double>(
        config, "maxwell_viscosity", parameters, 1);
 
    DBUG("Use '{:s}' as maxwell viscosity parameter.", maxwell_viscosity.name);
 
    // Dependency parameter for mK.
    auto& dependency_parameter_mK = ParameterLib::findParameter<double>(
        config, "dependency_parameter_mk", parameters, 1);
 
    DBUG("Use '{:s}' as dependency parameter mK.",
         dependency_parameter_mK.name);
 
    // Dependency parameter for mvK.
    auto& dependency_parameter_mvK = ParameterLib::findParameter<double>(
        config, "dependency_parameter_mvk", parameters, 1);
 
    DBUG("Use '{:s}' as dependency parameter mvK.",
         dependency_parameter_mvK.name);
 
    // Dependency parameter for mvM.
    auto& dependency_parameter_mvM = ParameterLib::findParameter<double>(
        config, "dependency_parameter_mvm", parameters, 1);
 
    DBUG("Use '{:s}' as dependency parameter mvM.",
         dependency_parameter_mvM.name);
 
    Lubby2MaterialProperties mp{
        kelvin_shear_modulus,     maxwell_shear_modulus,
        maxwell_bulk_modulus,     kelvin_viscosity,
        maxwell_viscosity,        dependency_parameter_mK,
        dependency_parameter_mvK, dependency_parameter_mvM};
 
    auto const& nonlinear_solver_config =
        config.getConfigSubtree("nonlinear_solver");
    auto const nonlinear_solver_parameters =
        NumLib::createNewtonRaphsonSolverParameters(nonlinear_solver_config);
 
    return std::unique_ptr<Lubby2<DisplacementDim>>{
        new Lubby2<DisplacementDim>{nonlinear_solver_parameters, mp}};
}
 
}  // namespace Lubby2
}  // namespace Solids
}  // namespace MaterialLib