coupling_penalty_condition.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Tobias Tescheamacher
//                   Michael Breitenberger
//                   Riccardo Rossi
//
 
#if !defined(KRATOS_COUPLING_PENALTY_CONDITION_H_INCLUDED )
#define  KRATOS_COUPLING_PENALTY_CONDITION_H_INCLUDED
 
// System includes
#include "includes/define.h"
#include "includes/condition.h"
 
// External includes
 
// Project includes
#include "iga_application_variables.h"
#include "custom_utilities/iga_flags.h"
 
#include "geometries/coupling_geometry.h"
 
namespace Kratos
{
 
 
class CouplingPenaltyCondition
    : public Condition
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(CouplingPenaltyCondition);
 
    typedef std::size_t SizeType;
    typedef std::size_t IndexType;
 
 
    CouplingPenaltyCondition(
        IndexType NewId,
        GeometryType::Pointer pGeometry)
        : Condition(NewId, pGeometry)
    {};
 
    CouplingPenaltyCondition(
        IndexType NewId,
        GeometryType::Pointer pGeometry,
        PropertiesType::Pointer pProperties)
        : Condition(NewId, pGeometry, pProperties)
    {};
 
    CouplingPenaltyCondition()
        : Condition()
    {};
 
    virtual ~CouplingPenaltyCondition() = default;
 
 
    Condition::Pointer Create(
        IndexType NewId,
        GeometryType::Pointer pGeom,
        PropertiesType::Pointer pProperties
    ) const override
    {
        return Kratos::make_intrusive<CouplingPenaltyCondition>(
            NewId, pGeom, pProperties);
    };
 
    Condition::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        PropertiesType::Pointer pProperties
    ) const override
    {
        return Kratos::make_intrusive< CouplingPenaltyCondition >(
            NewId, GetGeometry().Create(ThisNodes), pProperties);
    };
 
 
    void CalculateRightHandSide(
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override
    {
        MatrixType left_hand_side_matrix = Matrix(0, 0);
 
        CalculateAll(left_hand_side_matrix, rRightHandSideVector,
            rCurrentProcessInfo, false, true);
    }
 
    void CalculateLeftHandSide(
        MatrixType& rLeftHandSideMatrix,
        const ProcessInfo& rCurrentProcessInfo) override
    {
        VectorType right_hand_side_vector = Vector(0);
 
        CalculateAll(rLeftHandSideMatrix, right_hand_side_vector,
            rCurrentProcessInfo, true, false);
    }
 
    void CalculateLocalSystem(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override
    {
        CalculateAll(rLeftHandSideMatrix, rRightHandSideVector,
            rCurrentProcessInfo, true, true);
    }
 
    void EquationIdVector(
        EquationIdVectorType& rResult,
        const ProcessInfo& rCurrentProcessInfo
    ) const override;
 
    void GetDofList(
        DofsVectorType& rElementalDofList,
        const ProcessInfo& rCurrentProcessInfo
    ) const override;
 
    void CalculateAll(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo,
        const bool CalculateStiffnessMatrixFlag,
        const bool CalculateResidualVectorFlag
    );
 
    void DeterminantOfJacobianInitial(
        const GeometryType& rGeometry,
        Vector& rDeterminantOfJacobian);
 
 
    int Check(const ProcessInfo& rCurrentProcessInfo) const override;
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "\"CouplingPenaltyCondition\" #" << Id();
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "\"CouplingPenaltyCondition\" #" << Id();
    }
 
    void PrintData(std::ostream& rOStream) const override {
        pGetGeometry()->PrintData(rOStream);
    }
 
 
private:
    // Compute rotational shape functions
    void CalculateRotationalShapeFunctions(
        IndexType IntegrationPointIndex,
        Vector& phi_r,
        Matrix& phi_rs,
        array_1d<double, 2>& diff_phi);
 
    // Compute rotation
    void CalculateRotation(
        IndexType IntegrationPointIndex,
        const Matrix &rShapeFunctionGradientValues,
        Vector &phi_r,
        Matrix &phi_rs,
        array_1d<double, 2> &phi,
        array_1d<double, 3> &trim_tangent,
        const Vector &local_tangent,
        const bool master);
 
 
    friend class Serializer;
 
    virtual void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, Condition);
    }
 
    virtual void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, Condition);
    }
 
 
}; // Class CouplingPenaltyCondition
 
}  // namespace Kratos.
 
#endif // KRATOS_COUPLING_PENALTY_CONDITION_H_INCLUDED  defined