alm_variables_calculation_process.h

Go to the documentation of this file.

// KRATOS    ______            __             __  _____ __                  __                   __
//          / ____/___  ____  / /_____ ______/ /_/ ___// /________  _______/ /___  ___________ _/ /
//         / /   / __ \/ __ \/ __/ __ `/ ___/ __/\__ \/ __/ ___/ / / / ___/ __/ / / / ___/ __ `/ / 
//        / /___/ /_/ / / / / /_/ /_/ / /__/ /_ ___/ / /_/ /  / /_/ / /__/ /_/ /_/ / /  / /_/ / /  
//        \____/\____/_/ /_/\__/\__,_/\___/\__//____/\__/_/   \__,_/\___/\__/\__,_/_/   \__,_/_/  MECHANICS
//
//  License:         BSD License
//                   license: ContactStructuralMechanicsApplication/license.txt
//
//  Main authors:    Vicente Mataix Ferrandiz
//
 
#pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "includes/kratos_parameters.h"
#include "processes/process.h"
#include "includes/model_part.h"
#include "geometries/point.h"
 
namespace Kratos
{
 
 
 
 
class KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) ALMVariablesCalculationProcess
    : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(ALMVariablesCalculationProcess);
 
    using SizeType = std::size_t;
 
    using IndexType= std::size_t;
 
    using GeometryType = Geometry<Node>;
 
    using NodesArrayType = ModelPart::NodesContainerType;
 
    using ConditionsArrayType = ModelPart::ConditionsContainerType;
 
 
    ALMVariablesCalculationProcess(
        ModelPart& rThisModelPart,
        Variable<double>& rNodalLengthVariable = NODAL_H,
        Parameters ThisParameters =  Parameters(R"({})")
        ):mrThisModelPart(rThisModelPart),
          mrNodalLengthVariable(rNodalLengthVariable)
    {
        KRATOS_TRY
 
        ThisParameters.ValidateAndAssignDefaults(GetDefaultParameters());
 
        // Assign members variables from parameters
        mFactorStiffness = ThisParameters["stiffness_factor"].GetDouble();
        mPenaltyScale = ThisParameters["penalty_scale_factor"].GetDouble();
        mComputeScaleFactor = ThisParameters["compute_scale_factor"].GetBool();
        mComputePenalty = ThisParameters["compute_penalty"].GetBool();
 
        KRATOS_ERROR_IF_NOT(rThisModelPart.HasNodalSolutionStepVariable( rNodalLengthVariable )) << "Missing variable " << rNodalLengthVariable;
 
        KRATOS_CATCH("")
    }
 
    ~ALMVariablesCalculationProcess() override = default;
 
 
 
 
    void operator()()
    {
        Execute();
    }
 
 
    void Execute() override;
 
    const Parameters GetDefaultParameters() const override
    {
        const Parameters default_parameters = Parameters(R"(
        {
            "stiffness_factor"     : 10.0,
            "penalty_scale_factor" : 1.0,
            "compute_scale_factor" : true,
            "compute_penalty"      : true
        })" );
 
        return default_parameters;
    }
 
 
 
 
    std::string Info() const override
    {
        return "ALMVariablesCalculationProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "ALMVariablesCalculationProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
private:
 
 
    ModelPart& mrThisModelPart;              
    Variable<double>& mrNodalLengthVariable; 
    double mFactorStiffness;                 
    double mPenaltyScale;                    
    bool mComputeScaleFactor;                
    bool mComputePenalty;                    
 
 
 
 
 
 
    ALMVariablesCalculationProcess& operator=(ALMVariablesCalculationProcess const& rOther) = delete;
 
    //ALMVariablesCalculationProcess(ALMVariablesCalculationProcess const& rOther);
 
 
}; // Class ALMVariablesCalculationProcess
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  ALMVariablesCalculationProcess& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const ALMVariablesCalculationProcess& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
}