normal_gap_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 "processes/process.h"
#include "includes/model_part.h"
#include "processes/simple_mortar_mapper_process.h"
 
namespace Kratos
{
 
 
 
 
 
template<SizeType TDim, SizeType TNumNodes, SizeType TNumNodesMaster = TNumNodes>
class KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) NormalGapProcess
    : public Process
{
public:
 
    using MapperType = SimpleMortarMapperProcess<TDim, TNumNodes, Variable<array_1d<double, 3>>, TNumNodesMaster>;
 
    using NodesArrayType = ModelPart::NodesContainerType;
 
    static constexpr double ZeroTolerance = std::numeric_limits<double>::epsilon();
 
    KRATOS_CLASS_POINTER_DEFINITION( NormalGapProcess );
 
 
 
    NormalGapProcess(
        ModelPart& rMasterModelPart,
        ModelPart& rSlaveModelPart,
        const bool SearchOrientation = true
        ) : mrMasterModelPart(rMasterModelPart),
            mrSlaveModelPart(rSlaveModelPart),
            mSearchOrientation(SearchOrientation)
    {
    }
 
    virtual ~NormalGapProcess()= default;;
 
 
 
    void operator()()
    {
        Execute();
    }
 
 
    void Execute() override;
 
 
 
 
    /************************************ GET INFO *************************************/
    /***********************************************************************************/
 
    std::string Info() const override
    {
        return "NormalGapProcess";
    }
 
    /************************************ PRINT INFO ***********************************/
    /***********************************************************************************/
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << Info();
    }
 
 
 
protected:
 
 
 
    ModelPart& mrMasterModelPart;  
    ModelPart& mrSlaveModelPart;   
    const bool mSearchOrientation; 
 
 
 
    static inline void SwitchFlagNodes(NodesArrayType& rNodes)
    {
        block_for_each(rNodes, [&](Node& rNode) {
            rNode.Flip(SLAVE);
            rNode.Flip(MASTER);
        });
    }
 
    void ComputeNormalGap(NodesArrayType& rNodes);
 
 
 
 
 
private:
 
 
 
 
 
 
 
 
}; // Class NormalGapProcess
 
 
 
 
 
/****************************** INPUT STREAM FUNCTION ******************************/
/***********************************************************************************/
 
template<SizeType TDim, SizeType TNumNodes, SizeType TNumNodesMaster>
inline std::istream& operator >> (std::istream& rIStream,
                                  NormalGapProcess<TDim, TNumNodes, TNumNodesMaster>& rThis);
 
/***************************** OUTPUT STREAM FUNCTION ******************************/
/***********************************************************************************/
 
template<SizeType TDim, SizeType TNumNodes, SizeType TNumNodesMaster>
inline std::ostream& operator << (std::ostream& rOStream,
                                  const NormalGapProcess<TDim, TNumNodes, TNumNodesMaster>& rThis)
{
    return rOStream;
}
 
 
}  // namespace Kratos.