point_rigid_contact_penalty_3D_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_POINT_RIGID_CONTACT_PENALTY_3D_CONDITION_H_INCLUDED )
#define  KRATOS_POINT_RIGID_CONTACT_PENALTY_3D_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) PointRigidContactPenalty3DCondition
    : public PointRigidContactCondition
{
public:
 
 
    typedef BoundedVector<double, 3>     PointType;
 
    // Counted pointer of PointRigidContactCondition
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( PointRigidContactPenalty3DCondition );
 
 
    PointRigidContactPenalty3DCondition(){};
 
    PointRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry);
 
    PointRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
    PointRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties, SpatialBoundingBox::Pointer pRigidWall);
 
 
    PointRigidContactPenalty3DCondition( PointRigidContactPenalty3DCondition const& rOther);
 
 
    virtual ~PointRigidContactPenalty3DCondition();
 
 
 
 
 
    Condition::Pointer Create(IndexType NewId, NodesArrayType const&
                              ThisNodes,  PropertiesType::Pointer pProperties) const override;
 
    Condition::Pointer Clone(IndexType NewId,
                 NodesArrayType const& ThisNodes) const override;
 
 
protected:
 
 
 
 
    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& CalculateEffectiveNormalForceModulus( 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);
 
 
private:
 
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, PointRigidContactCondition )
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, PointRigidContactCondition )
    }
 
 
}; // Class PointRigidContactPenalty3DCondition
 
 
 
 
 
 
/*  inline std::istream& operator >> (std::istream& rIStream,
                    PointRigidContactPenalty3DCondition& rThis);
*/
/*  inline std::ostream& operator << (std::ostream& rOStream,
                    const PointRigidContactPenalty3DCondition& rThis)
    {
      rThis.PrintInfo(rOStream);
      rOStream << std::endl;
      rThis.PrintData(rOStream);
 
      return rOStream;
    }*/
 
}  // namespace Kratos.
 
#endif // KRATOS_POINT_RIGID_CONTACT_PENALTY_3D_CONDITION_H_INCLUDED  defined