elastic_cohesive_3D_law.hpp

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Ignasi de Pouplana and Danilo Cavalcanti
//
 
#if !defined (KRATOS_ELASTIC_COHESIVE_3D_LAW_H_INCLUDED)
#define  KRATOS_ELASTIC_COHESIVE_3D_LAW_H_INCLUDED
 
// System includes
#include <cmath>
 
// Project includes
#include "includes/serializer.h"
#include "includes/checks.h"
#include "includes/constitutive_law.h"
 
// Application includes
#include "poromechanics_application_variables.h"
 
namespace Kratos
{
 
class KRATOS_API(POROMECHANICS_APPLICATION) ElasticCohesive3DLaw : public ConstitutiveLaw
{
 
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(ElasticCohesive3DLaw);
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    // Default Constructor
    ElasticCohesive3DLaw()
    {
    }
 
    ConstitutiveLaw::Pointer Clone() const override
    {
        return Kratos::make_shared<ElasticCohesive3DLaw>(ElasticCohesive3DLaw(*this));
    }
 
    // Copy Constructor
    ElasticCohesive3DLaw (const ElasticCohesive3DLaw& rOther) : ConstitutiveLaw(rOther)
    {
    }
 
    // Destructor
    ~ElasticCohesive3DLaw() override
    {
    }
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    void GetLawFeatures(Features& rFeatures) override;
 
    int Check(const Properties& rMaterialProperties, const GeometryType& rElementGeometry, const ProcessInfo& rCurrentProcessInfo) const override;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    void CalculateMaterialResponseCauchy(Parameters & rValues) override;
 
    void FinalizeMaterialResponseCauchy(Parameters & rValues) override;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
protected:
 
    struct ConstitutiveLawVariables
    {
        double NormalStiffness;
        double ShearStiffness;
        double PenaltyStiffness; 
    };
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    virtual void InitializeConstitutiveLawVariables(ConstitutiveLawVariables& rVariables, Parameters& rValues);
 
    virtual void ComputeConstitutiveMatrix(Matrix& rConstitutiveMatrix,
                                            ConstitutiveLawVariables& rVariables,
                                            Parameters& rValues);
 
    virtual void ComputeStressVector(Vector& rStressVector,
                                        ConstitutiveLawVariables& rVariables,
                                        Parameters& rValues);
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
private:
 
    // Serialization
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, ConstitutiveLaw )
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, ConstitutiveLaw )
    }
 
}; // Class ElasticCohesive3DLaw
}  // namespace Kratos.
#endif // KRATOS_ELASTIC_COHESIVE_3D_LAW_H_INCLUDED  defined