hydraulic_rigid_contact_penalty_3D_condition.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosContactMechanicsApplication $
//   Created by:          $Author:                LMonforte $
//   Last modified by:    $Co-Author:                       $
//   Date:                $Date:               January 2018 $
//   Revision:            $Revision:                    0.0 $
//
//
 
#if !defined(KRATOS_HYDRAULIC_RIGID_CONTACT_PENALTY_CONDITION_H_INCLUDED )
#define  KRATOS_HYDRAULIC_RIGID_CONTACT_PENALTY_CONDITION_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_conditions/rigid_contact/point_rigid_contact_condition.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
class KRATOS_API(CONTACT_MECHANICS_APPLICATION) HydraulicRigidContactPenalty3DCondition
    : public PointRigidContactCondition
{
public:
 
 
    typedef BoundedVector<double, 3>     PointType;
 
    // Counted pointer of WaterPointRigidContactCondition
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( HydraulicRigidContactPenalty3DCondition );
 
 
    HydraulicRigidContactPenalty3DCondition(){};
 
    HydraulicRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry);
 
    HydraulicRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
    HydraulicRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties, SpatialBoundingBox::Pointer pRigidWall);
 
 
    HydraulicRigidContactPenalty3DCondition( HydraulicRigidContactPenalty3DCondition const& rOther);
 
 
    virtual ~HydraulicRigidContactPenalty3DCondition();
 
 
 
 
 
    Condition::Pointer Create(IndexType NewId, NodesArrayType const&
                              ThisNodes,  PropertiesType::Pointer pProperties) const override;
 
    Condition::Pointer Clone(IndexType NewId,
                 NodesArrayType const& ThisNodes) const override;
 
 
 
    virtual void InitializeNonLinearIteration(const ProcessInfo& rCurrentProcessInfo) override;
 
    //************* GETTING METHODS
 
    void GetDofList(DofsVectorType& rConditionDofList,
            const ProcessInfo& rCurrentProcessInfo ) const override;
 
    void EquationIdVector(EquationIdVectorType& rResult,
              const ProcessInfo& rCurrentProcessInfo ) const override;
 
    void GetValuesVector(Vector& rValues,
             int Step = 0 ) const override;
 
    void GetFirstDerivativesVector(Vector& rValues,
                   int Step = 0 ) const override;
 
    void GetSecondDerivativesVector(Vector& rValues,
                    int Step = 0 ) const override;
 
 
 
 
protected:
 
 
 
 
    virtual void CalculateKinematics(ConditionVariables& rVariables,
                     const ProcessInfo& rCurrentProcessInfo,
                     const double& rPointNumber) override;
 
 
 
    virtual void CalculateAndAddKuug(MatrixType& rLeftHandSideMatrix,
                     ConditionVariables& rVariables,
                     double& rIntegrationWeight) override;
 
    virtual void CalculateAndAddContactForces(Vector& rRightHandSideVector,
                          ConditionVariables& rVariables,
                          double& rIntegrationWeight ) override;
 
 
    void CalculateAndAddNormalContactForce(Vector& rRightHandSideVector, ConditionVariables& rVariables, double& rIntegrationWeight);
 
    double& CalculateNormalForceModulus( double& rNormalForceModulus, ConditionVariables& rVariables );
 
 
    virtual void CalculateContactFactors(ConditionVariables &rContact);
 
 
private:
 
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
    virtual void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, PointRigidContactCondition )
    }
 
    virtual void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, PointRigidContactCondition )
    }
 
 
}; // Class HydraulicRigidContactPenalty3DCondition
 
 
 
 
 
 
/*  inline std::istream& operator >> (std::istream& rIStream,
                    HydraulicRigidContactPenalty3DCondition& rThis);
*/
/*  inline std::ostream& operator << (std::ostream& rOStream,
                    const HydraulicRigidContactPenalty3DCondition& rThis)
    {
      rThis.PrintInfo(rOStream);
      rOStream << std::endl;
      rThis.PrintData(rOStream);
 
      return rOStream;
    }*/
 
}  // namespace Kratos.
 
#endif // KRATOS_HYDRAULIC_RIGID_CONTACT_PENALTY_CONDITION_H_INCLUDED  defined