advanced_contact_search_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 "custom_processes/base_contact_search_process.h"
 
namespace Kratos
{
 
 
    using SizeType = std::size_t;
 
 
 
 
template<SizeType TDim, SizeType TNumNodes, SizeType TNumNodesMaster = TNumNodes>
class KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) AdvancedContactSearchProcess
    : public BaseContactSearchProcess<TDim, TNumNodes, TNumNodesMaster>
{
public:
 
    using BaseType = BaseContactSearchProcess<TDim, TNumNodes, TNumNodesMaster>;
 
    using NodesArrayType = typename BaseType::NodesArrayType;
    using ConditionsArrayType = typename BaseType::ConditionsArrayType;
    using GeometryType = typename BaseType::GeometryType;
 
    using IndexType= std::size_t;
 
    static constexpr double ZeroTolerance = std::numeric_limits<double>::epsilon();
 
    static constexpr double GapThreshold = 2.0e-4;
 
    KRATOS_CLASS_POINTER_DEFINITION( AdvancedContactSearchProcess );
 
 
 
    AdvancedContactSearchProcess(
        ModelPart& rMainModelPart,
        Parameters ThisParameters =  Parameters(R"({})"),
        Properties::Pointer pPairedProperties = nullptr
        );
 
    virtual ~AdvancedContactSearchProcess()= default;;
 
 
 
 
 
 
    /************************************ GET INFO *************************************/
    /***********************************************************************************/
 
    std::string Info() const override
    {
        return "AdvancedContactSearchProcess";
    }
 
    /************************************ PRINT INFO ***********************************/
    /***********************************************************************************/
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << Info();
    }
 
 
 
protected:
 
 
 
 
 
    void ComputeActiveInactiveNodes() override;
 
    void CheckPairing(
        ModelPart& rComputingModelPart,
        IndexType& rConditionId
        ) override;
 
 
 
 
 
private:
 
 
 
 
    void ComputeLinearRegressionGapPressure(
        double& a,
        double& b
        );
 
    void SetActiveNodeWithRegression(
        Node& rNode,
        const double a,
        const double b
        );
 
    void CorrectScalarMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
    void CorrectComponentsMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
    void CorrectALMFrictionlessMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
    void CorrectALMFrictionlessComponentsMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
    void CorrectALMFrictionalMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
    void PredictScalarMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
    void PredictComponentsMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
    void PredictALMFrictionlessMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
    void PredictALMFrictionlessComponentsMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
    void PredictALMFrictionalMortarLM(
        Node& rNode,
        const double a,
        const double b
        );
 
 
 
 
 
}; // Class AdvancedContactSearchProcess
 
 
 
 
 
/****************************** INPUT STREAM FUNCTION ******************************/
/***********************************************************************************/
 
template<SizeType TDim, SizeType TNumNodes, SizeType TNumNodesMaster = TNumNodes>
inline std::istream& operator >> (std::istream& rIStream,
                                  AdvancedContactSearchProcess<TDim, TNumNodes, TNumNodesMaster>& rThis);
 
/***************************** OUTPUT STREAM FUNCTION ******************************/
/***********************************************************************************/
 
template<SizeType TDim, SizeType TNumNodes, SizeType TNumNodesMaster = TNumNodes>
inline std::ostream& operator << (std::ostream& rOStream,
                                  const AdvancedContactSearchProcess<TDim, TNumNodes, TNumNodesMaster>& rThis)
{
    return rOStream;
}
 
 
}  // namespace Kratos.