mortar_explicit_contribution_utilities.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 "utilities/math_utils.h"
#include "custom_conditions/paired_condition.h"
#include "includes/mortar_classes.h"
 
/* Utilities */
#include "utilities/exact_mortar_segmentation_utility.h"
#include "custom_utilities/derivatives_utilities.h"
 
/* Geometries */
#include "geometries/line_2d_2.h"
#include "geometries/triangle_3d_3.h"
 
namespace Kratos
{
 
 
 
 
 
template< const SizeType TDim, const SizeType TNumNodes, const FrictionalCase TFrictional, const bool TNormalVariation, const SizeType TNumNodesMaster = TNumNodes>
class KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) MortarExplicitContributionUtilities
{
public:
 
    using SizeType = std::size_t;
 
    using IndexType = std::size_t;
 
    using PointType = Point;
 
    using GeometryType = Geometry<Node>;
 
    // Type definition for integration methods
    using IntegrationPointsType = GeometryType::IntegrationPointsArrayType;
 
    using BelongType = typename std::conditional<TNumNodes == 2, PointBelongsLine2D2N, typename std::conditional<TNumNodes == 3, typename std::conditional<TNumNodesMaster == 3, PointBelongsTriangle3D3N, PointBelongsTriangle3D3NQuadrilateral3D4N>::type, typename std::conditional<TNumNodesMaster == 3, PointBelongsQuadrilateral3D4NTriangle3D3N, PointBelongsQuadrilateral3D4N>::type>::type>::type;
 
    using PointBelongType = PointBelong<TNumNodes, TNumNodesMaster>;
 
    using GeometryPointBelongType = Geometry<PointBelongType>;
 
    using ConditionArrayType = array_1d<PointBelongType, TDim>;
 
    using ConditionArrayListType = typename std::vector<ConditionArrayType>;
 
    using LineType = Line2D2<PointType>;
 
    using TriangleType = Triangle3D3<PointType>;
 
    using DecompositionType = typename std::conditional<TDim == 2, LineType, TriangleType>::type;
 
    using DerivativeDataType = typename std::conditional<TFrictional == FrictionalCase::FRICTIONAL || TFrictional == FrictionalCase::FRICTIONAL_PENALTY, DerivativeDataFrictional<TDim, TNumNodes, TNumNodesMaster>, DerivativeData<TDim, TNumNodes, TNumNodesMaster>>::type;
 
    static constexpr IndexType MatrixSize = (TFrictional == FrictionalCase::FRICTIONLESS) ? TDim * (TNumNodesMaster + TNumNodes) + TNumNodes : (TFrictional == FrictionalCase::FRICTIONLESS_COMPONENTS || TFrictional == FrictionalCase::FRICTIONAL) ? TDim * (TNumNodesMaster + TNumNodes + TNumNodes) : TDim * (TNumNodesMaster + TNumNodes);
 
    static constexpr bool IsFrictional = (TFrictional == FrictionalCase::FRICTIONAL || TFrictional == FrictionalCase::FRICTIONAL_PENALTY) ? true : false;
 
    using GeneralVariables = MortarKinematicVariablesWithDerivatives<TDim, TNumNodes, TNumNodesMaster>;
 
    using AeData = DualLagrangeMultiplierOperatorsWithDerivatives<TDim, TNumNodes, IsFrictional, TNumNodesMaster>;
 
    using MortarConditionMatrices = MortarOperatorWithDerivatives<TDim, TNumNodes, IsFrictional, TNumNodesMaster>;
 
    using IntegrationUtility = ExactMortarIntegrationUtility<TDim, TNumNodes, true, TNumNodesMaster>;
 
    using DerivativesUtilitiesType = DerivativesUtilities<TDim, TNumNodes, IsFrictional, TNormalVariation, TNumNodesMaster>;
 
    // The threshold coefficient considered for checking
    static constexpr double CheckThresholdCoefficient = 1.0e-12;
 
 
 
 
    static MortarConditionMatrices AddExplicitContributionOfMortarCondition(
        PairedCondition* pCondition,
        const ProcessInfo& rCurrentProcessInfo,
        const IndexType IntegrationOrder = 2,
        const bool AxisymmetricCase = false,
        const bool ComputeNodalArea = false,
        const bool ComputeDualLM = true,
        const Variable<double>& rAreaVariable = NODAL_AREA
        );
    static MortarConditionMatrices AddExplicitContributionOfMortarFrictionalCondition(
        PairedCondition* pCondition,
        const ProcessInfo& rCurrentProcessInfo,
        const MortarOperator<TNumNodes, TNumNodesMaster>& rPreviousMortarOperators,
        const IndexType IntegrationOrder = 2,
        const bool AxisymmetricCase = false,
        const bool ComputeNodalArea = false,
        const bool ComputeDualLM = true,
        const Variable<double>& rAreaVariable = NODAL_AREA,
        const bool ConsiderObjetiveFormulation = false
        );
 
    static bool ExplicitCalculateAe(
        const GeometryType& rSlaveGeometry,
        GeneralVariables& rVariables,
        const ConditionArrayListType& rConditionsPointsSlave,
        BoundedMatrix<double, TNumNodes, TNumNodes>& rAe,
        const IntegrationMethod& rIntegrationMethod,
        const double AxiSymCoeff = 1.0
        );
 
    static void ExplicitCalculateKinematics(
        const PairedCondition* pCondition,
        GeneralVariables& rVariables,
        const BoundedMatrix<double, TNumNodes, TNumNodes>& rAe,
        const array_1d<double, 3>& rNormalMaster,
        const PointType& rLocalPointDecomp,
        const PointType& rLocalPointParent,
        const GeometryPointType& rGeometryDecomp,
        const bool DualLM = true
        );
 
    static void ComputeNodalArea(
        PairedCondition* pCondition,
        const ProcessInfo& rCurrentProcessInfo,
        const Variable<double>& rAreaVariable = NODAL_AREA,
        const IndexType IntegrationOrder = 2,
        const bool AxisymmetricCase = false
        );
 
    static bool ComputePreviousMortarOperators(
        PairedCondition* pCondition,
        const ProcessInfo& rCurrentProcessInfo,
        MortarOperator<TNumNodes, TNumNodesMaster>& rPreviousMortarOperators,
        const IndexType IntegrationOrder = 2,
        const bool AxisymmetricCase = false,
        const bool ComputeNodalArea = false,
        const bool ComputeDualLM = true,
        const Variable<double>& rAreaVariable = NODAL_AREA
        );
 
    static void CalculateKinematics(
        const PairedCondition* pCondition,
        GeneralVariables& rVariables,
        const DerivativeDataType& rDerivativeData,
        const array_1d<double, 3>& rNormalMaster,
        const PointType& rLocalPointDecomp,
        const PointType& rLocalPointParent,
        const GeometryPointType& rGeometryDecomp,
        const bool DualLM = true
        );
 
    static void MasterShapeFunctionValue(
        const PairedCondition* pCondition,
        GeneralVariables& rVariables,
        const array_1d<double, 3>& rNormalMaster,
        const PointType& rLocalPoint
        );
 
}; // class MortarExplicitContributionUtilities
 
namespace AuxiliaryOperationsUtilities
{
    double KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) GetAxisymmetricCoefficient(
        const PairedCondition* pCondition,
        const Vector& rNSlave
        );
 
    double KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) CalculateRadius(
        const PairedCondition* pCondition,
        const Vector& rNSlave
        );
}
 
}  // namespace Kratos