contact_structural_mechanics_application.h

Go to the documentation of this file.

// KRATOS    ______            __             __  _____ __                  __                   __
//          / ____/___  ____  / /_____ ______/ /_/ ___// /________  _______/ /___  ___________ _/ /
//         / /   / __ \/ __ \/ __/ __ `/ ___/ __/\__ \/ __/ ___/ / / / ___/ __/ / / / ___/ __ `/ / 
//        / /___/ /_/ / / / / /_/ /_/ / /__/ /_ ___/ / /_/ /  / /_/ / /__/ /_/ /_/ / /  / /_/ / /  
//        \____/\____/_/ /_/\__/\__,_/\___/\__//____/\__/_/   \__,_/\___/\__/\__,_/_/   \__,_/_/  MECHANICS
//
//  License:         BSD License
//                   license: ContactStructuralMechanicsApplication/license.txt
//
//  Main authors:    Vicente Mataix
//
 
#pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "includes/kratos_application.h"
#include "includes/variables.h"
 
/* GEOMETRIES */
#include "geometries/triangle_3d_3.h"
#include "geometries/quadrilateral_3d_4.h"
#include "geometries/line_2d_2.h"
 
/* CONDITIONS */
// Mortar conditions
#include "custom_conditions/mesh_tying_mortar_condition.h"
#include "custom_conditions/ALM_frictionless_mortar_contact_condition.h"
#include "custom_conditions/ALM_frictionless_components_mortar_contact_condition.h"
#include "custom_conditions/penalty_frictionless_mortar_contact_condition.h"
#include "custom_conditions/ALM_frictionless_mortar_contact_axisym_condition.h"
#include "custom_conditions/penalty_frictionless_mortar_contact_axisym_condition.h"
#include "custom_conditions/ALM_frictional_mortar_contact_condition.h"
#include "custom_conditions/penalty_frictional_mortar_contact_condition.h"
#include "custom_conditions/ALM_frictional_mortar_contact_axisym_condition.h"
#include "custom_conditions/penalty_frictional_mortar_contact_axisym_condition.h"
#include "custom_conditions/mpc_mortar_contact_condition.h"
 
/* CONSTRAINTS */
#include "custom_master_slave_constraints/contact_master_slave_constraint.h"
 
namespace Kratos
{
 
 
 
 
 
 
class KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) KratosContactStructuralMechanicsApplication
    : public KratosApplication
{
public:
 
    using GeometryType = Geometry<Node>;
 
    using GeometryPointerType = typename GeometryType::Pointer;
 
    using PointsArrayType = typename GeometryType::PointsArrayType;
 
    using LineType = Line2D2<Node>;
 
    using TriangleType = Triangle3D3<Node>;
 
    using QuadrilateralType = Quadrilateral3D4<Node>;
 
    KRATOS_CLASS_POINTER_DEFINITION(KratosContactStructuralMechanicsApplication);
 
 
    KratosContactStructuralMechanicsApplication();
 
    ~KratosContactStructuralMechanicsApplication() override = default;
 
 
 
 
    void Register() override;
 
 
 
 
 
 
    std::string Info() const override
    {
        return "KratosContactStructuralMechanicsApplication";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << Info();
        PrintData(rOStream);
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
        KRATOS_WATCH(KratosComponents<VariableData>::GetComponents().size() );
        rOStream << "Variables:" << std::endl;
        KratosComponents<VariableData>().PrintData(rOStream);
        rOStream << std::endl;
        rOStream << "Elements:" << std::endl;
        KratosComponents<Element>().PrintData(rOStream);
        rOStream << std::endl;
        rOStream << "Conditions:" << std::endl;
        KratosComponents<Condition>().PrintData(rOStream);
    }
 
 
 
protected:
 
 
 
 
 
 
 
private:
 
 
    /* CONDITIONS*/
    // Mesh tying mortar condition
    const MeshTyingMortarCondition<2, 2> mMeshTyingMortarCondition2D2N;      // 2D Line/Line
    const MeshTyingMortarCondition<3, 3> mMeshTyingMortarCondition3D3N;      // 3D Triangle/Triangle
    const MeshTyingMortarCondition<3, 4> mMeshTyingMortarCondition3D4N;      // 3D Quadrilateral/Quadrilateral
    const MeshTyingMortarCondition<3, 3, 4> mMeshTyingMortarCondition3D3N4N; // 3D Triangle/Quadrilateral
    const MeshTyingMortarCondition<3, 4, 3> mMeshTyingMortarCondition3D4N3N; // 3D Quadrilateral/Triangle
 
    // ALM Mortar contact conditions
    // Frictionless cases
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<2, 2, false> mALMFrictionlessMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<2, 2, true> mALMNVFrictionlessMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionlessMortarContactAxisymCondition<2, false> mALMFrictionlessAxisymMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionlessMortarContactAxisymCondition<2, true> mALMNVFrictionlessAxisymMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<3, 3, false, 3> mALMFrictionlessMortarContactCondition3D3N;
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<3, 3, true,  3> mALMNVFrictionlessMortarContactCondition3D3N;
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<3, 4, false, 4> mALMFrictionlessMortarContactCondition3D4N;
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<3, 4, true,  4> mALMNVFrictionlessMortarContactCondition3D4N;
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<3, 3, false, 4> mALMFrictionlessMortarContactCondition3D3N4N;
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<3, 3, true,  4> mALMNVFrictionlessMortarContactCondition3D3N4N;
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<3, 4, false, 3> mALMFrictionlessMortarContactCondition3D4N3N;
    const AugmentedLagrangianMethodFrictionlessMortarContactCondition<3, 4, true,  3> mALMNVFrictionlessMortarContactCondition3D4N3N;
    // Frictionless components cases
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<2, 2, false> mALMFrictionlessComponentsMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<2, 2, true> mALMNVFrictionlessComponentsMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<3, 3, false, 3> mALMFrictionlessComponentsMortarContactCondition3D3N;
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<3, 3, true,  3> mALMNVFrictionlessComponentsMortarContactCondition3D3N;
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<3, 3, false, 4> mALMFrictionlessComponentsMortarContactCondition3D3N4N;
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<3, 3, true,  4> mALMNVFrictionlessComponentsMortarContactCondition3D3N4N;
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<3, 4, false, 4> mALMFrictionlessComponentsMortarContactCondition3D4N;
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<3, 4, true,  4> mALMNVFrictionlessComponentsMortarContactCondition3D4N;
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<3, 4, false, 3> mALMFrictionlessComponentsMortarContactCondition3D4N3N;
    const AugmentedLagrangianMethodFrictionlessComponentsMortarContactCondition<3, 4, true,  3> mALMNVFrictionlessComponentsMortarContactCondition3D4N3N;
    // Frictional cases
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<2, 2, false> mALMFrictionalMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<2, 2, true> mALMNVFrictionalMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionalMortarContactAxisymCondition<2, false> mALMFrictionalAxisymMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionalMortarContactAxisymCondition<2, true> mALMNVFrictionalAxisymMortarContactCondition2D2N;
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<3, 3, false, 3> mALMFrictionalMortarContactCondition3D3N;
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<3, 3, true,  3> mALMNVFrictionalMortarContactCondition3D3N;
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<3, 4, false, 4> mALMFrictionalMortarContactCondition3D4N;
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<3, 4, true,  4> mALMNVFrictionalMortarContactCondition3D4N;
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<3, 3, false, 4> mALMFrictionalMortarContactCondition3D3N4N;
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<3, 3, true,  4> mALMNVFrictionalMortarContactCondition3D3N4N;
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<3, 4, false, 3> mALMFrictionalMortarContactCondition3D4N3N;
    const AugmentedLagrangianMethodFrictionalMortarContactCondition<3, 4, true,  3> mALMNVFrictionalMortarContactCondition3D4N3N;
    // Frictionless penalty cases
    const PenaltyMethodFrictionlessMortarContactCondition<2, 2, false> mPenaltyFrictionlessMortarContactCondition2D2N;
    const PenaltyMethodFrictionlessMortarContactCondition<2, 2, true> mPenaltyNVFrictionlessMortarContactCondition2D2N;
    const PenaltyMethodFrictionlessMortarContactAxisymCondition<2, false> mPenaltyFrictionlessAxisymMortarContactCondition2D2N;
    const PenaltyMethodFrictionlessMortarContactAxisymCondition<2, true> mPenaltyNVFrictionlessAxisymMortarContactCondition2D2N;
    const PenaltyMethodFrictionlessMortarContactCondition<3, 3, false, 3> mPenaltyFrictionlessMortarContactCondition3D3N;
    const PenaltyMethodFrictionlessMortarContactCondition<3, 3, true,  3> mPenaltyNVFrictionlessMortarContactCondition3D3N;
    const PenaltyMethodFrictionlessMortarContactCondition<3, 4, false, 4> mPenaltyFrictionlessMortarContactCondition3D4N;
    const PenaltyMethodFrictionlessMortarContactCondition<3, 4, true,  4> mPenaltyNVFrictionlessMortarContactCondition3D4N;
    const PenaltyMethodFrictionlessMortarContactCondition<3, 3, false, 4> mPenaltyFrictionlessMortarContactCondition3D3N4N;
    const PenaltyMethodFrictionlessMortarContactCondition<3, 3, true,  4> mPenaltyNVFrictionlessMortarContactCondition3D3N4N;
    const PenaltyMethodFrictionlessMortarContactCondition<3, 4, false, 3> mPenaltyFrictionlessMortarContactCondition3D4N3N;
    const PenaltyMethodFrictionlessMortarContactCondition<3, 4, true,  3> mPenaltyNVFrictionlessMortarContactCondition3D4N3N;
    // Frictional penalty cases
    const PenaltyMethodFrictionalMortarContactCondition<2, 2, false> mPenaltyFrictionalMortarContactCondition2D2N;
    const PenaltyMethodFrictionalMortarContactCondition<2, 2, true> mPenaltyNVFrictionalMortarContactCondition2D2N;
    const PenaltyMethodFrictionalMortarContactAxisymCondition<2, false> mPenaltyFrictionalAxisymMortarContactCondition2D2N;
    const PenaltyMethodFrictionalMortarContactAxisymCondition<2, true> mPenaltyNVFrictionalAxisymMortarContactCondition2D2N;
    const PenaltyMethodFrictionalMortarContactCondition<3, 3, false, 3> mPenaltyFrictionalMortarContactCondition3D3N;
    const PenaltyMethodFrictionalMortarContactCondition<3, 3, true,  3> mPenaltyNVFrictionalMortarContactCondition3D3N;
    const PenaltyMethodFrictionalMortarContactCondition<3, 4, false, 4> mPenaltyFrictionalMortarContactCondition3D4N;
    const PenaltyMethodFrictionalMortarContactCondition<3, 4, true,  4> mPenaltyNVFrictionalMortarContactCondition3D4N;
    const PenaltyMethodFrictionalMortarContactCondition<3, 3, false, 4> mPenaltyFrictionalMortarContactCondition3D3N4N;
    const PenaltyMethodFrictionalMortarContactCondition<3, 3, true,  4> mPenaltyNVFrictionalMortarContactCondition3D3N4N;
    const PenaltyMethodFrictionalMortarContactCondition<3, 4, false, 3> mPenaltyFrictionalMortarContactCondition3D4N3N;
    const PenaltyMethodFrictionalMortarContactCondition<3, 4, true,  3> mPenaltyNVFrictionalMortarContactCondition3D4N3N;
 
    // MPC Conditions
    const MPCMortarContactCondition<2, 2> mMPCMortarContactCondition2D2N;
    const MPCMortarContactCondition<3, 3, 3> mMPCMortarContactCondition3D3N;
    const MPCMortarContactCondition<3, 4, 4> mMPCMortarContactCondition3D4N;
    const MPCMortarContactCondition<3, 3, 4> mMPCMortarContactCondition3D3N4N;
    const MPCMortarContactCondition<3, 4, 3> mMPCMortarContactCondition3D4N3N;
 
    const ContactMasterSlaveConstraint mContactMasterSlaveConstraint;
 
 
 
 
 
 
 
 
 
 
    KratosContactStructuralMechanicsApplication& operator=(KratosContactStructuralMechanicsApplication const& rOther);
 
    KratosContactStructuralMechanicsApplication(KratosContactStructuralMechanicsApplication const& rOther);
 
 
 
}; // Class KratosContactStructuralMechanicsApplication
 
 
 
 
 
}  // namespace Kratos.