EP_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_EP_POINT_RIGID_CONTACT_PENALTY_3D_CONDITION_H_INCLUDED )
#define  KRATOS_EP_POINT_RIGID_CONTACT_PENALTY_3D_CONDITION_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_conditions/rigid_contact/point_rigid_contact_penalty_3D_condition.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
class KRATOS_API(CONTACT_MECHANICS_APPLICATION) EPPointRigidContactPenalty3DCondition
    : public PointRigidContactPenalty3DCondition
{
 
   protected:
      typedef struct
      {
         Vector PreviousStepForceVector;
         Vector t1;
         Vector t2;
         Vector n;
      } GeometricalInformation;
 
      typedef struct
      {
         // ConstitutiveInformation
         double TangentForceRatio;
         double  NormalTangentMatrix;
         double TangentTangentMatrix;
         Vector ForceDirection;
      }  ConstitutiveVariables;
 
   public:
 
 
    typedef BoundedVector<double, 3>     PointType;
 
    // Counted pointer of PointRigidContactCondition
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( EPPointRigidContactPenalty3DCondition );
 
 
    EPPointRigidContactPenalty3DCondition(){};
 
    EPPointRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry);
 
    EPPointRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
    EPPointRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties, SpatialBoundingBox::Pointer pRigidWall);
 
 
    EPPointRigidContactPenalty3DCondition( EPPointRigidContactPenalty3DCondition const& rOther);
 
 
    virtual ~EPPointRigidContactPenalty3DCondition();
 
 
 
 
 
    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 InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
    virtual void FinalizeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
 
    virtual void InitializeNonLinearIteration(const ProcessInfo& rCurrentProcessInfo) override;
 
 
protected:
 
 
    GeometricalInformation mCurrentInfo;
    GeometricalInformation mSavedInfo;
 
    double mElasticYoungModulus; // using MCC + IMPLEX, at the finalizeSolutionStep the Contact Forces are correctly ingrated with the finalized of the Young Modulus at the continuum elements, that is quite different from the one used during the implex step.
 
    bool mImplex;
 
 
 
    void CalculateAndAddKuugTangent(MatrixType& rLeftHandSideMatrix,
                                    ConditionVariables& rVariables,
                                    double& rIntegrationWeight) override;
 
 
 
 
    virtual void CalculateAndAddTangentContactForce(Vector& rRightHandSideVector, ConditionVariables& rVariables, double& rIntegrationWeight) override;
 
 
    bool CalculateFrictionLaw( ConditionVariables & rVariables, ConstitutiveVariables & rConstitutiveVariables, Vector & rTangentForce);
 
    virtual double CalculateSomeSortOfArea();
 
    double CalculateEffectiveNormalForceModulus( const double& rNormalForceModulus, const double & rArea);
 
    Matrix ConvertToTheAppropriateSize( const Matrix & rForceMatrix);
 
    virtual void CalculateContactFactors(ConditionVariables &rContact) override;
 
private:
 
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
    virtual void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, PointRigidContactPenalty3DCondition )
    }
 
    virtual void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, PointRigidContactPenalty3DCondition )
    }
 
 
}; // Class EPPointRigidContactPenalty3DCondition
 
 
 
 
 
 
/*  inline std::istream& operator >> (std::istream& rIStream,
                    EPPointRigidContactPenalty3DCondition& rThis);
*/
/*  inline std::ostream& operator << (std::ostream& rOStream,
                    const EPPointRigidContactPenalty3DCondition& 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