beam_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 2013 $
//   Revision:            $Revision:                    0.0 $
//
//
 
#if !defined(KRATOS_BEAM_POINT_RIGID_CONTACT_PENALTY_3D_CONDITION_H_INCLUDED )
#define  KRATOS_BEAM_POINT_RIGID_CONTACT_PENALTY_3D_CONDITION_H_INCLUDED
 
 
// System includes
 
// External includes
 
// Project includes
#include "custom_conditions/thermal_contact/beam_point_rigid_contact_condition.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
class KRATOS_API(CONTACT_MECHANICS_APPLICATION) BeamPointRigidContactPenalty3DCondition
    : public BeamPointRigidContactCondition
{
protected:
 
    typedef struct
     {
        bool   Slip;
        double Sign;
 
        double DeltaTime;
        double PreviousTangentForceModulus;
 
        double FrictionCoefficient;
        double DynamicFrictionCoefficient;
        double StaticFrictionCoefficient;
 
     } TangentialContactVariables;
 
 
 
public:
 
 
    //typedef BoundedVector<double, 3>     PointType;
    typedef array_1d<double, 3>             PointType;
 
    // Counted pointer of BeamPointRigidContactCondition
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( BeamPointRigidContactPenalty3DCondition );
 
 
    BeamPointRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry);
 
    BeamPointRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
    BeamPointRigidContactPenalty3DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties, SpatialBoundingBox::Pointer pRigidWall);
 
 
    BeamPointRigidContactPenalty3DCondition( BeamPointRigidContactPenalty3DCondition const& rOther);
 
 
    virtual ~BeamPointRigidContactPenalty3DCondition();
 
 
 
 
 
    Condition::Pointer Create(IndexType NewId, NodesArrayType const&
                              ThisNodes,  PropertiesType::Pointer pProperties) const;
 
    virtual void InitializeSolutionStep(ProcessInfo& rCurrentProcessInfo);
 
 
    virtual void InitializeNonLinearIteration(ProcessInfo& rCurrentProcessInfo);
 
 
protected:
 
    // A protected default constructor necessary for serialization
    BeamPointRigidContactPenalty3DCondition() {};
 
 
    TangentialContactVariables mTangentialVariables;
 
 
 
    void InitializeConditionVariables(ConditionVariables& rVariables,
                    const ProcessInfo& rCurrentProcessInfo);
 
    virtual void CalculateKinematics(ConditionVariables& rVariables,
                     const ProcessInfo& rCurrentProcessInfo,
                     const double& rPointNumber);
 
 
    virtual void CalculateAndAddKuug(MatrixType& rLeftHandSideMatrix,
                     ConditionVariables& rVariables,
                     double& rIntegrationWeight);
 
    virtual void CalculateAndAddKuugTangent(MatrixType& rLeftHandSideMatrix,
                     ConditionVariables& rVariables,
                     double& rIntegrationWeight);
 
    virtual void CalculateAndAddContactForces(Vector& rRightHandSideVector,
                          ConditionVariables& rVariables,
                          double& rIntegrationWeight );
 
 
    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& CalculateTangentRelativeMovement( double& rTangentRelativeMovement, ConditionVariables& rVariables );
 
    double CalculateCoulombsFrictionLaw( double& rTangentForceModulus, double& rNormalForceModulus, ConditionVariables& rVariables );
 
    double CalculateFrictionCoefficient(double & rTangentRelativeMovement);
 
 
    virtual void CalculateContactFactors(ConditionVariables &rContact);
 
 
private:
 
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
    virtual void save(Serializer& rSerializer) const
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, BeamPointRigidContactCondition )
    }
 
    virtual void load(Serializer& rSerializer)
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, BeamPointRigidContactCondition )
    }
 
 
}; // Class BeamPointRigidContactPenalty3DCondition
 
 
 
 
 
 
/*  inline std::istream& operator >> (std::istream& rIStream,
                    BeamPointRigidContactPenalty3DCondition& rThis);
*/
/*  inline std::ostream& operator << (std::ostream& rOStream,
                    const BeamPointRigidContactPenalty3DCondition& 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