penalty_frictional_mortar_contact_condition.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_utilities/contact_utilities.h"
#include "custom_conditions/mortar_contact_condition.h"
 
namespace Kratos
{
 
 
 
    using PointType = Point;
    using GeometryType = Geometry<Node>;
    using GeometryPointType = Geometry<PointType>;
    using IntegrationMethod = GeometryData::IntegrationMethod;
 
 
 
 
template< std::size_t TDim, std::size_t TNumNodes, bool TNormalVariation, std::size_t TNumNodesMaster = TNumNodes>
class KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) PenaltyMethodFrictionalMortarContactCondition
    : public MortarContactCondition<TDim, TNumNodes, FrictionalCase::FRICTIONAL_PENALTY, TNormalVariation, TNumNodesMaster>
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( PenaltyMethodFrictionalMortarContactCondition );
 
    using BaseType = MortarContactCondition<TDim, TNumNodes, FrictionalCase::FRICTIONAL_PENALTY, TNormalVariation, TNumNodesMaster>;
 
    using ConditionBaseType = Condition;
 
    using PairedConditionBaseType = PairedCondition;
 
    using MortarConditionMatrices = typename BaseType::MortarConditionMatrices;
 
    using GeneralVariables = typename BaseType::GeneralVariables;
 
    using IntegrationUtility = typename BaseType::IntegrationUtility;
 
    using DerivativesUtilitiesType = typename BaseType::DerivativesUtilitiesType;
 
    using BelongType = typename BaseType::BelongType;
 
    using ConditionArrayListType = typename BaseType::ConditionArrayListType;
 
    using MortarBaseConditionMatrices = MortarOperator<TNumNodes, TNumNodesMaster>;
 
    using VectorType = typename ConditionBaseType::VectorType;
 
    using MatrixType = typename ConditionBaseType::MatrixType;
 
    using IndexType = typename ConditionBaseType::IndexType;
 
    using GeometryPointerType = typename ConditionBaseType::GeometryType::Pointer;
 
    using NodesArrayType = typename ConditionBaseType::NodesArrayType;
 
    using PropertiesType = typename ConditionBaseType::PropertiesType;
 
    using PropertiesPointerType = typename ConditionBaseType::PropertiesType::Pointer;
 
    using EquationIdVectorType = typename ConditionBaseType::EquationIdVectorType;
 
    using DofsVectorType = typename ConditionBaseType::DofsVectorType;
 
    using LineType = Line2D2<Point>;
 
    using TriangleType = Triangle3D3<Point>;
 
    using DecompositionType = typename std::conditional<TDim == 2, LineType, TriangleType>::type;
 
    using DerivativeDataType = DerivativeDataFrictional<TDim, TNumNodes, TNumNodesMaster>;
 
    static constexpr IndexType MatrixSize = TDim * (TNumNodes + TNumNodesMaster);
 
    static constexpr IndexType StepSlip = TNormalVariation ? 0 : 1;
 
 
    PenaltyMethodFrictionalMortarContactCondition()
        : BaseType()
    {
    }
 
    // Constructor 1
    PenaltyMethodFrictionalMortarContactCondition(
        IndexType NewId,
        GeometryPointerType pGeometry
        ):BaseType(NewId, pGeometry)
    {
    }
 
    // Constructor 2
    PenaltyMethodFrictionalMortarContactCondition(
        IndexType NewId,
        GeometryPointerType pGeometry,
        PropertiesPointerType pProperties
        ):BaseType( NewId, pGeometry, pProperties )
    {
    }
 
    // Constructor 3
    PenaltyMethodFrictionalMortarContactCondition(
        IndexType NewId,
        GeometryPointerType pGeometry,
        PropertiesPointerType pProperties,
        GeometryType::Pointer pMasterGeometry
        ):BaseType( NewId, pGeometry, pProperties, pMasterGeometry )
    {
    }
 
    PenaltyMethodFrictionalMortarContactCondition( PenaltyMethodFrictionalMortarContactCondition const& rOther)
    {
    }
 
    ~PenaltyMethodFrictionalMortarContactCondition() override;
 
 
 
 
    void Initialize(const ProcessInfo& rCurrentProcessInfo) override;
 
    void InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
    void FinalizeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
    Condition::Pointer Create(
        IndexType NewId,
        NodesArrayType const& rThisNodes,
        PropertiesPointerType pProperties
        ) const override;
 
    Condition::Pointer Create(
        IndexType NewId,
        GeometryPointerType pGeom,
        PropertiesPointerType pProperties
        ) const override;
 
    Condition::Pointer Create(
        IndexType NewId,
        GeometryPointerType pGeom,
        PropertiesPointerType pProperties,
        GeometryPointerType pMasterGeom
        ) const override;
 
    void AddExplicitContribution(const ProcessInfo& rCurrentProcessInfo) override;
 
    void AddExplicitContribution(
        const VectorType& rRHSVector,
        const Variable<VectorType>& rRHSVariable,
        const Variable<double >& rDestinationVariable,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
    void AddExplicitContribution(const VectorType& rRHSVector,
        const Variable<VectorType>& rRHSVariable,
        const Variable<array_1d<double, 3> >& rDestinationVariable,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
    /********************************************************************************/
    /**************** METHODS TO CALCULATE MORTAR CONDITION MATRICES ****************/
    /********************************************************************************/
 
    static void StaticCalculateLocalRHS(
        PairedCondition* pCondition,
        const MortarBaseConditionMatrices& rPreviousMortarOperators,
        const array_1d<double, TNumNodes>& mu,
        Vector& rLocalRHS,
        const MortarConditionMatrices& rMortarConditionMatrices,
        const DerivativeDataType& rDerivativeData,
        const IndexType rActiveInactive,
        const ProcessInfo& rCurrentProcessInfo
        );
 
    /******************************************************************/
    /********** AUXILLIARY METHODS FOR GENERAL CALCULATIONS ***********/
    /******************************************************************/
 
    void EquationIdVector(
        EquationIdVectorType& rResult,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
    void GetDofList(
        DofsVectorType& rConditionalDofList,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
    int Check(const ProcessInfo& rCurrentProcessInfo) const override;
 
 
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "PenaltyMethodFrictionalMortarContactCondition #" << this->Id();
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "PenaltyMethodFrictionalMortarContactCondition #" << this->Id();
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
        PrintInfo(rOStream);
        this->GetParentGeometry().PrintData(rOStream);
        this->GetPairedGeometry().PrintData(rOStream);
    }
 
 
 
protected:
 
 
    bool mPreviousMortarOperatorsInitialized = false;     
 
    MortarBaseConditionMatrices mPreviousMortarOperators; 
 
    // TODO: Define the "CL" or friction law to compute this. Or do it nodally
 
 
 
    /********************************************************************************/
    /**************** METHODS TO CALCULATE MORTAR CONDITION MATRICES ****************/
    /********************************************************************************/
 
    void CalculateLocalLHS(
        Matrix& rLocalLHS,
        const MortarConditionMatrices& rMortarConditionMatrices,
        const DerivativeDataType& rDerivativeData,
        const IndexType rActiveInactive,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
    void CalculateLocalRHS(
        Vector& rLocalRHS,
        const MortarConditionMatrices& rMortarConditionMatrices,
        const DerivativeDataType& rDerivativeData,
        const IndexType rActiveInactive,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
    /******************************************************************/
    /********** AUXILLIARY METHODS FOR GENERAL CALCULATIONS ***********/
    /******************************************************************/
 
    IndexType GetActiveInactiveValue(const GeometryType& rCurrentGeometry) const override
    {
        IndexType value = 0;
        for (IndexType i_node = 0; i_node < TNumNodes; ++i_node) {
            if (rCurrentGeometry[i_node].Is(ACTIVE) == true) {
                if (rCurrentGeometry[i_node].Is(SLIP) == true)
                    value += std::pow(3, i_node);
                else
                    value += 2 * std::pow(3, i_node);
            }
        }
 
        return value;
    }
 
    array_1d<double, TNumNodes> GetFrictionCoefficient()
    {
        // The friction coefficient
        array_1d<double, TNumNodes> friction_coeffient_vector;
        auto& r_geometry = this->GetParentGeometry();
 
        for (std::size_t i_node = 0; i_node < TNumNodes; ++i_node) {
            friction_coeffient_vector[i_node] = r_geometry[i_node].GetValue(FRICTION_COEFFICIENT);
        }
 
        // TODO: Define the "CL" or friction law to compute this
 
        return friction_coeffient_vector;
    }
 
 
 
 
private:
 
 
 
 
    void ComputePreviousMortarOperators(const ProcessInfo& rCurrentProcessInfo);
 
    static inline BoundedMatrix<double, TNumNodes, TDim> ComputeTangentMatrixSlip(const GeometryType& rGeometry)
    {
        return MortarUtilities::ComputeTangentMatrix<TDim, TNumNodes>(rGeometry);
//         return ContactUtilities::ComputeTangentMatrixSlip<TDim, TNumNodes>(rGeometry, StepSlip);
    }
 
 
 
 
    // Serialization
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, BaseType );
        rSerializer.save("PreviousMortarOperators", mPreviousMortarOperators);
        rSerializer.save("PreviousMortarOperatorsInitialized", mPreviousMortarOperatorsInitialized);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, BaseType );
        rSerializer.load("PreviousMortarOperators", mPreviousMortarOperators);
        rSerializer.load("PreviousMortarOperatorsInitialized", mPreviousMortarOperatorsInitialized);
    }
 
 
}; // Class PenaltyMethodFrictionalMortarContactCondition
 
 
 
 
 
}// namespace Kratos.