interface_preprocess.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
#include <iostream>
 
// External includes
 
// Project includes
#include "includes/model_part.h"
#include "includes/define.h"
#include "includes/kratos_parameters.h"
 
namespace Kratos
{
 
 
 
 
 
class KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) InterfacePreprocessCondition
{
public:
 
    using PointType = Point;
    using GeometryType = Geometry<Node>;
    using GeometryPointType = Geometry<PointType>;
 
    using IndexType = std::size_t;
 
    using SizeType = std::size_t;
 
    using NodesArrayType = ModelPart::NodesContainerType;
    using ElementsArrayType = ModelPart::ElementsContainerType;
    using ConditionsArrayType = ModelPart::ConditionsContainerType;
    
    KRATOS_CLASS_POINTER_DEFINITION(InterfacePreprocessCondition);
 
 
 
    InterfacePreprocessCondition(ModelPart& rMainModelPrt)
    :mrMainModelPart(rMainModelPrt)
    {
    }
 
    virtual ~InterfacePreprocessCondition() = default;
 
 
 
    void GenerateInterfacePart(
        ModelPart& rInterfacePart,
        Parameters ThisParameters =  Parameters(R"({})")
        );
 
protected:
 
 
 
 
 
 
private:
 
 
    ModelPart& mrMainModelPart; 
 
 
 
    void CheckAndCreateProperties(ModelPart& rInterfacePart);
 
    bool CheckOnTheFace(
        const std::vector<std::size_t>& rIndexVector,
        GeometryType& rElementGeometry
        );
 
    std::unordered_map<IndexType, Properties::Pointer> CreateNewProperties();
 
    template<class TClass>
    void CopyProperties(
        Properties::Pointer pOriginalProperty,
        Properties::Pointer pNewProperty,
        const Variable<TClass>& rVariable,
        const bool AssignZero = true
        )
    {
        if(pOriginalProperty->Has(rVariable)) {
            const TClass& value = pOriginalProperty->GetValue(rVariable);
            pNewProperty->SetValue(rVariable, value);
        } else if (AssignZero) {
            KRATOS_INFO("InterfacePreprocessCondition") << "Property " << rVariable.Name() << " not available. Assigning zero value" << std::endl;
            pNewProperty->SetValue(rVariable, rVariable.Zero());
        }
    }
 
    void CreateNewCondition(
        Properties::Pointer prThisProperties,
        const GeometryType& rGeometry,
        const IndexType ConditionId,
        Condition const& rCondition
        );
 
    void AssignMasterSlaveCondition(Condition::Pointer pCond);
 
    void PrintNodesAndConditions(
        const IndexType NodesCounter,
        const IndexType rCondCounter
        );
 
    IndexType ReorderConditions();
 
    inline void GenerateEdgeCondition(
        ModelPart& rInterfacePart,
        Properties::Pointer prThisProperties,
        const GeometryType& rEdgeGeometry,
        const bool SimplestGeometry,
        IndexType& rCondCounter,
        IndexType& rConditionId
        );
 
    inline void GenerateFaceCondition(
        ModelPart& rInterfacePart,
        Properties::Pointer prThisProperties,
        const GeometryType& rFaceGeometry,
        const bool SimplestGeometry,
        IndexType& rCondCounter,
        IndexType& rConditionId
        );
 
 
 
 
}; // Class InterfacePreprocessCondition
}