thermal_contact_domain_condition.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosMachiningApplication $
//   Created by:          $Author:       JMCarbonell $
//   Last modified by:    $Co-Author:                $
//   Date:                $Date:        October 2013 $
//   Revision:            $Revision:             0.0 $
//
//
 
#if !defined(KRATOS_THERMAL_CONTACT_DOMAIN_CONDITION_H_INCLUDED )
#define  KRATOS_THERMAL_CONTACT_DOMAIN_CONDITION_H_INCLUDED
 
// System includes
#include <iostream>
#include <iomanip>
 
// External includes
#include "boost/smart_ptr.hpp"
 
 
// Project includes
#include "includes/define.h"
#include "includes/element.h"
#include "includes/serializer.h"
#include "includes/condition.h"
#include "includes/ublas_interface.h"
#include "includes/variables.h"
#include "includes/constitutive_law.h"
#include "includes/condition.h"
 
#include "custom_utilities/contact_domain_utilities.hpp"
 
namespace Kratos
{
 
class KRATOS_API(CONTACT_MECHANICS_APPLICATION) ThermalContactDomainCondition
    : public Condition
{
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 GlobalPointersVector<Node > NodeWeakPtrVectorType;
    typedef GlobalPointersVector<Element> ElementWeakPtrVectorType;
    typedef GlobalPointersVector<Condition> ConditionWeakPtrVectorType;
 
protected:
 
    typedef struct
    {
        Flags           Options;      //calculation options
 
        double          ThermalGap;   //thermal gap
 
 
        //Geometric variables
    Vector           ProjectionsVector;
 
    SurfaceVector        CurrentSurface;      //Current Normal and Tangent
        SurfaceVector        ReferenceSurface;    //Reference Normal and Tangent
 
        std::vector<BaseLengths>   CurrentBase;    //Current Base Lengths variables
        std::vector<BaseLengths>   ReferenceBase;  //Reference Base Lengths variables
 
 
        //Axisymmetric
        double  CurrentRadius;
        double  ReferenceRadius;
 
 
        //friction:
        double          RelativeVelocityNorm;
    double          FrictionForceNorm;
 
 
    } GeneralVariables;
 
 
    struct ContactVariables
    {
 
       //The stabilization parameter and penalty parameter
        double          StabilizationFactor;
 
        //Contact condition conectivities
        std::vector<unsigned int> nodes;
        std::vector<unsigned int> order;
        std::vector<unsigned int> slaves;
 
    //Pointer Variables
        GeometryType*         mpMasterGeometry;
        ElementType*          mpMasterElement;
        Condition*            mpMasterCondition;
        NodeType*             mpMasterNode;
 
    void SetMasterGeometry  (GeometryType& rGeometry){ mpMasterGeometry = &rGeometry; }
        void SetMasterElement   (ElementType& rElement){ mpMasterElement = &rElement; }
        void SetMasterCondition (ConditionType& rCondition){ mpMasterCondition = &rCondition; }
        void SetMasterNode      (NodeType& rNode){ mpMasterNode = &rNode; }
 
        GeometryType& GetMasterGeometry()   { return (*mpMasterGeometry); }
    ElementType& GetMasterElement()     { return (*mpMasterElement); }
    ConditionType& GetMasterCondition() { return (*mpMasterCondition); }
    NodeType& GetMasterNode()           { return (*mpMasterNode); }
 
 
    };
 
 
 
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(ThermalContactDomainCondition);
 
 
 
    ThermalContactDomainCondition() : Condition() {};
 
    ThermalContactDomainCondition(IndexType NewId, GeometryType::Pointer pGeometry);
 
    ThermalContactDomainCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
    ThermalContactDomainCondition(ThermalContactDomainCondition const& rOther);
 
 
    virtual ~ThermalContactDomainCondition();
 
 
    ThermalContactDomainCondition& operator=(ThermalContactDomainCondition 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;
 
    //************* GETTING METHODS
 
    IntegrationMethod GetIntegrationMethod() const override;
 
    void GetDofList(DofsVectorType& rConditionalDofList, 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;
 
 
 
    //on integration points:
    //SET
    void SetValuesOnIntegrationPoints(const Variable<double>& rVariable, const std::vector<double>& rValues, const ProcessInfo& rCurrentProcessInfo) override;
    void SetValuesOnIntegrationPoints(const Variable<Vector>& rVariable, const std::vector<Vector>& rValues, const ProcessInfo& rCurrentProcessInfo) override;
 
    void SetValuesOnIntegrationPoints(const Variable<Matrix>& rVariable, const std::vector<Matrix>& rValues, const ProcessInfo& rCurrentProcessInfo) override;
 
 
    //************* STARTING - ENDING  METHODS
 
    void Initialize(const ProcessInfo& CurrentProcessInfo) override;
 
 
    void InitializeSolutionStep(const ProcessInfo& CurrentProcessInfo) override;
 
    void InitializeNonLinearIteration(const ProcessInfo& CurrentProcessInfo) override;
 
    void FinalizeSolutionStep(const ProcessInfo& CurrentProcessInfo) override;
 
 
 
    //************* COMPUTING  METHODS
 
    void CalculateLocalSystem(MatrixType& rLeftHandSideMatrix, VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateLeftHandSide (MatrixType& rLeftHandSideMatrix, const ProcessInfo& rCurrentProcessInfo) override;
 
 
    //on integration points:
    void CalculateOnIntegrationPoints(const Variable<double>& rVariable, std::vector<double>& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(const Variable<Vector>& rVariable, std::vector<Vector>& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(const Variable<Matrix >& rVariable, std::vector< Matrix >& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
 
 
 
    //************************************************************************************
    //************************************************************************************
    int Check(const ProcessInfo& rCurrentProcessInfo) const override;
 
    //std::string Info() const;
 
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "Thermal Contact Domain Condition #" << Id();
        return buffer.str();
 
    }
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "Thermal Contact Domain Condition #" << Id();
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
      GetGeometry().PrintData(rOStream);
    }
 
 
protected:
 
    IntegrationMethod mThisIntegrationMethod;
 
 
    ContactVariables mContactVariables;
 
 
    ContactDomainUtilities  mContactUtilities;
 
 
    virtual void SetMasterGeometry()
    {
      KRATOS_THROW_ERROR( std::invalid_argument, "Calling base class in contact domain", "" );
    };
 
    void CalculateHeatConductivity();
 
 
    virtual void CalculateConditionalSystem(MatrixType& rLeftHandSideMatrix,
                                            VectorType& rRightHandSideVector,
                                            const ProcessInfo& rCurrentProcessInfo,
                        Flags& rCalculationFlags);
 
 
 
    virtual void CalculateKinematics(GeneralVariables& rVariables,
                     const ProcessInfo& rCurrentProcessInfo,
                     const unsigned int& rPointNumber)
    {
      KRATOS_THROW_ERROR( std::invalid_argument, "Calling base class in contact domain", "" );
    };
 
 
 
    void InitializeVariables ();
 
 
 
    virtual void InitializeSystemMatrices(MatrixType& rLeftHandSideMatrix,
                      VectorType& rRightHandSideVector,
                      Flags& rCalculationFlags);
 
 
    virtual void CalculateAndAddLHS(MatrixType& rLeftHandSideMatrix,
                    GeneralVariables& rVariables,
                    double& rIntegrationWeight);
 
    virtual void CalculateAndAddRHS(VectorType& rRightHandSideVector,
                                    GeneralVariables& rVariables,
                                    double& rIntegrationWeight);
 
    virtual void CalculateAndAddThermalKm(MatrixType& rLeftHandSideMatrix,
                      GeneralVariables& rVariables,
                      double& rIntegrationWeight);
 
    virtual void CalculateAndAddThermalContactForces(VectorType& rRightHandSideVector,
                             GeneralVariables& rVariables,
                             double& rIntegrationWeight);
 
 
    virtual double& CalculateIntegrationWeight(double& rIntegrationWeight);
 
 
    virtual void CalculateThermalFrictionForce(double &F, GeneralVariables& rVariables, unsigned int& ndi)
    {
        KRATOS_THROW_ERROR( std::invalid_argument, "Calling base class in contact domain", "" );
 
    };
 
 
    virtual void CalculateThermalConductionForce(double &F, GeneralVariables& rVariables, unsigned int& ndi)
    {
        KRATOS_THROW_ERROR( std::invalid_argument, "Calling base class in contact domain", "" );
 
    };
 
 
    void CalculateRelativeVelocity(GeneralVariables& rVariables, PointType& TangentVelocity, const ProcessInfo& rCurrentProcessInfo);
 
 
    void CalculateRelativeDisplacement(GeneralVariables& rVariables, PointType & TangentDisplacement, const ProcessInfo& rCurrentProcessInfo);
 
 
    virtual PointType & CalculateCurrentTangent(PointType &rTangent)
    {
        KRATOS_THROW_ERROR( std::invalid_argument, "Calling base class in contact domain", "" );
 
    };
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
private:
 
 
}; // Class ThermalContactDomainCondition
 
 
} // namespace Kratos.
#endif // KRATOS_THERMAL_CONTACT_DOMAIN_CONDITION_H_INCLUDED  defined