rigid_body_point_rigid_contact_condition.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosContactMechanicsApplication $
//   Created by:          $Author:              JMCarbonell $
//   Last modified by:    $Co-Author:                       $
//   Date:                $Date:                  July 2016 $
//   Revision:            $Revision:                    0.0 $
//
//
 
#if !defined(KRATOS_RIGID_BODY_POINT_RIGID_CONTACT_CONDITION_H_INCLUDED )
#define  KRATOS_RIGID_BODY_POINT_RIGID_CONTACT_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) RigidBodyPointRigidContactCondition
    : public PointRigidContactCondition
{
public:
 
 
    //typedef BoundedVector<double, 3>     PointType;
    typedef array_1d<double, 3>             PointType;
 
    // Counted pointer of RigidBodyPointRigidContactCondition
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( RigidBodyPointRigidContactCondition );
 
 
    RigidBodyPointRigidContactCondition( IndexType NewId, GeometryType::Pointer pGeometry );
 
    RigidBodyPointRigidContactCondition( IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties );
 
    RigidBodyPointRigidContactCondition( IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties, SpatialBoundingBox::Pointer pRigidWall );
 
    RigidBodyPointRigidContactCondition( RigidBodyPointRigidContactCondition const& rOther);
 
    virtual ~RigidBodyPointRigidContactCondition();
 
 
 
 
    Condition::Pointer Create(IndexType NewId,
                  NodesArrayType const& ThisNodes,
                  PropertiesType::Pointer pProperties ) const override;
 
 
    Condition::Pointer Clone(IndexType NewId,
                 NodesArrayType const& ThisNodes) const 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;
 
 
    void AddExplicitContribution(const VectorType& rRHSVector,
                     const Variable<VectorType>& rRHSVariable,
                     const Variable<array_1d<double,3> >& rDestinationVariable,
                     const ProcessInfo& rCurrentProcessInfo) override;
 
    //************************************************************************************
    //************************************************************************************
    int Check( const ProcessInfo& rCurrentProcessInfo ) const override;
 
 
protected:
    RigidBodyPointRigidContactCondition() {};
 
 
    SpatialBoundingBox::Pointer mpRigidWall;
 
 
 
 
    void InitializeSystemMatrices(MatrixType& rLeftHandSideMatrix,
                      VectorType& rRightHandSideVector,
                      Flags& rCalculationFlags) override;
 
    void CalculateKinematics(ConditionVariables& rVariables,
                     const ProcessInfo& rCurrentProcessInfo,
                     const double& rPointNumber) override;
 
 
    void CalculateAndAddKuug(MatrixType& rLeftHandSideMatrix,
                     ConditionVariables& rVariables,
                     double& rIntegrationWeight) override;
 
    virtual void CalculateAndAddKuugTangent(MatrixType& rLeftHandSideMatrix,
                     ConditionVariables& rVariables,
                     double& rIntegrationWeight);
 
    void CalculateAndAddContactForces(Vector& rRightHandSideVector,
                          ConditionVariables& rVariables,
                          double& rIntegrationWeight ) override;
 
 
 
    virtual void CalculateAndAddNormalContactForce(Vector& rRightHandSideVector, ConditionVariables& rVariables, double& rIntegrationWeight);
 
 
    virtual void CalculateAndAddTangentContactForce(Vector& rRightHandSideVector, ConditionVariables& rVariables, double& rIntegrationWeight);
 
 
    double& CalculateNormalForceModulus( double& rNormalForceModulus, ConditionVariables& rVariables );
 
    double CalculateCoulombsFrictionLaw( double& rTangentForceModulus, double& rNormalForceModulus, ConditionVariables& rVariables );
 
    double CalculateFrictionCoefficient( const double& rTangentRelativeMovement, const double& rDeltaTime );
 
 
    virtual void CalculateContactFactors(ConditionVariables &rContact);
 
 
     void VectorToSkewSymmetricTensor( const Vector& rVector,
                       Matrix& rSkewSymmetricTensor );
 
 
 
 
 
 
 
 
 
private:
 
 
 
    ElementWeakPtrVectorType mMasterElements;
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
    void save( Serializer& rSerializer ) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, Condition )
    }
 
    void load( Serializer& rSerializer ) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, Condition )
    }
 
 
}; // class RigidBodyPointRigidContactCondition.
 
} // namespace Kratos.
 
#endif // KRATOS_RIGID_BODY_POINT_RIGID_CONTACT_CONDITION_H_INCLUDED defined