contact_domain_LM_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_CONTACT_DOMAIN_LM_3D_CONDITION_H_INCLUDED )
#define  KRATOS_CONTACT_DOMAIN_LM_3D_CONDITION_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_conditions/deformable_contact/contact_domain_condition.hpp"
 
namespace Kratos
{
 
 
 
 
class KRATOS_API(CONTACT_MECHANICS_APPLICATION) ContactDomainLM3DCondition
    : public ContactDomainCondition
{
public:
 
 
    typedef ConstitutiveLaw ConstitutiveLawType;
    typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
    typedef GeometryData::IntegrationMethod IntegrationMethod;
 
    typedef Node NodeType;
    typedef Geometry<NodeType> GeometryType;
    typedef Element::ElementType ElementType;
 
    typedef ContactDomainUtilities::PointType             PointType;
    typedef ContactDomainUtilities::SurfaceVector      SurfaceVector;
    typedef ContactDomainUtilities::SurfaceScalar      SurfaceScalar;
    typedef ContactDomainUtilities::BaseLengths          BaseLengths;
 
    typedef ContactDomainUtilities::SurfaceBase           SurfaceBase;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( ContactDomainLM3DCondition );
 
 
 
    ContactDomainLM3DCondition() : ContactDomainCondition() {};
 
    ContactDomainLM3DCondition(IndexType NewId, GeometryType::Pointer pGeometry);
 
    ContactDomainLM3DCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
    ContactDomainLM3DCondition(ContactDomainLM3DCondition const& rOther);
 
 
    virtual ~ContactDomainLM3DCondition();
 
 
    ContactDomainLM3DCondition& operator=(ContactDomainLM3DCondition const& rOther);
 
 
 
    Condition::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override;
 
    Condition::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override;
 
 
    //************************************************************************************
    //************************************************************************************
    //int Check(const ProcessInfo& rCurrentProcessInfo);
 
    //std::string Info() const;
 
 
 
    //      String Info() const override;
 
    //      void PrintInfo(std::ostream& rOStream) const override;
 
    //      void PrintData(std::ostream& rOStream) const override;
 
protected:
 
 
    void ResolveElementType();
 
 
    void SetMasterGeometry() override;
 
 
    void CalculateContactFactor(const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void CalculatePreviousGap() override;
 
    void CalculatePreviousGapEdgeType();
 
    void CalculatePreviousGapFaceType();
 
    void CalculateExplicitFactors(ConditionVariables& rVariables,
                                  const ProcessInfo& rCurrentProcessInfo) override;
 
    virtual void CalculateExplicitFactorsEdgeType(ConditionVariables& rVariables,
                          const ProcessInfo& rCurrentProcessInfo);
 
    virtual void CalculateExplicitFactorsFaceType(ConditionVariables& rVariables,
                          const ProcessInfo& rCurrentProcessInfo);
 
 
    void CalculateDomainShapeN(ConditionVariables& rVariables) override;
 
 
    double& CalculateIntegrationWeight(double& rIntegrationWeight) override;
 
    void CalculateContactStiffness (double &Kcont,ConditionVariables& rVariables,
                                    unsigned int& ndi,unsigned int& ndj,
                                    unsigned int& idir,unsigned int& jdir) override;
 
 
    void CalculateNormalForce       (double &F,ConditionVariables& rVariables,
                         unsigned int& ndi,unsigned int& idir) override;
 
    void CalculateTangentStickForce (double &F,ConditionVariables& rVariables,
                         unsigned int& ndi,unsigned int& idir) override;
    void CalculateTangentSlipForce  (double &F,ConditionVariables& rVariables,
                         unsigned int& ndi,unsigned int& idir) override;
 
 
    inline bool CheckFictiousContacts(ConditionVariables& rVariables);
 
    PointType& CalculateCurrentTangent(PointType &rTangent) override;
 
    void FSigmaP(ConditionVariables& rVariables, std::vector<Vector >& rSigmaP, PointType& rDirVector,unsigned int &ndi,unsigned int &ndj,unsigned int &ndk,unsigned int &ndl,unsigned int &ndm,unsigned int &ndn);
 
    void FSigmaPnd(ConditionVariables& rVariables, std::vector<Vector >& rSigmaP, PointType& rDirVector,unsigned int &ndi,unsigned int &ndj);
 
 
 
    void TransformCovariantToContravariantBase(SurfaceBase& Covariant,SurfaceBase& Contravariant);
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
private:
 
 
}; // Class ContactDomainLM3DCondition
 
/*  inline std::istream& operator >> (std::istream& rIStream,
    ContactDomainLM3DCondition& rThis);
*/
/*  inline std::ostream& operator << (std::ostream& rOStream,
    const ContactDomainLM3DCondition& rThis)
    {
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
    }*/
 
} // namespace Kratos.
#endif // KRATOS_CONTACT_DOMAIN_LM_3D_CONDITION_H_INCLUDED  defined