elastoplastic_mohr_coulomb_cohesive_3D_law.hpp

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Danilo Cavalcanti and Ignasi de Pouplana
//
 
#if !defined (KRATOS_ELASTOPLASTIC_MOHRCOULOMB_COHESIVE_3D_LAW_H_INCLUDED)
#define  KRATOS_ELASTOPLASTIC_MOHRCOULOMB_COHESIVE_3D_LAW_H_INCLUDED
 
// System includes
#include <cmath>
 
// Project includes
#include "utilities/math_utils.h"
#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) ElastoPlasticMohrCoulombCohesive3DLaw : public ConstitutiveLaw
{
 
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(ElastoPlasticMohrCoulombCohesive3DLaw);
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    // Default Constructor
    ElastoPlasticMohrCoulombCohesive3DLaw()
    {
    }
 
    ConstitutiveLaw::Pointer Clone() const override
    {
        return Kratos::make_shared<ElastoPlasticMohrCoulombCohesive3DLaw>(ElastoPlasticMohrCoulombCohesive3DLaw(*this));
    }
 
    // Copy Constructor
    ElastoPlasticMohrCoulombCohesive3DLaw (const ElastoPlasticMohrCoulombCohesive3DLaw& rOther) : ConstitutiveLaw(rOther)
    {
    }
 
    // Destructor
    ~ElastoPlasticMohrCoulombCohesive3DLaw() override
    {
    }
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    void GetLawFeatures(Features& rFeatures) override;
 
    int Check(const Properties& rMaterialProperties, const GeometryType& rElementGeometry, const ProcessInfo& rCurrentProcessInfo) const override;
 
    void InitializeMaterial( const Properties& rMaterialProperties,const GeometryType& rElementGeometry,const Vector& rShapeFunctionsValues ) override;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    void CalculateMaterialResponseCauchy(Parameters & rValues) override;
 
    void FinalizeMaterialResponseCauchy(Parameters & rValues) override;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
    Vector& GetValue( const Variable<Vector>& rThisVariable, Vector& rValue ) override;
 
    void SetValue( const Variable<Vector>& rThisVariable,
                   const Vector& rValue,
                   const ProcessInfo& rCurrentProcessInfo ) override;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
protected:
 
    struct ConstitutiveLawVariables
    {
        // Material parameters
        double NormalStiffness;
        double ShearStiffness;
        double PenaltyStiffness;
        double TensileStrength; 
        double FrictionAngle;
        double DilatancyAngle;
        double Cohesion;
    };
 
    struct ElastoPlasticConstitutiveLawVariables
    {
        // Yield surfaces
        double YieldFunction_MC;
        double YieldFunction_CutOff;
 
        // Vector normal to the yield surface
        Vector n_MC;
        Vector n_TC;
 
        // Vector normal to the plastic potential surface
        Vector np_MC;
        Vector np_TC;
    };
 
    // Member Variables
    Vector mPlasticStrainVector;    
    Vector mOldPlasticStrainVector;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    virtual void InitializeConstitutiveLawVariables(ConstitutiveLawVariables& rVariables, Parameters& rValues);
 
    virtual void InitializeElastoPlasticConstitutiveLawVariables(ElastoPlasticConstitutiveLawVariables& rEPlasticVariables, const double VoigtSize);
 
    virtual void ComputeYieldFunction(Vector& StressVector, ConstitutiveLawVariables& rVariables, ElastoPlasticConstitutiveLawVariables& rEPlasticVariables, Parameters& rValues);
 
    virtual double GetShearResultantStressVector(Vector& StressVector);
 
    virtual void ReturnMapping(Vector& rStressVector, Matrix& rConstitutiveMatrix, Vector& TrialStressVector, Matrix& ElasticConstitutiveMatrix, ConstitutiveLawVariables& rVariables, ElastoPlasticConstitutiveLawVariables& rEPlasticVariables, Parameters& rValues);
 
    virtual void DerivativesPlasticPotentialSurface(Vector& StressVector, ConstitutiveLawVariables& rVariables, ElastoPlasticConstitutiveLawVariables& rEPlasticVariables, Parameters& rValues);
 
    virtual void DerivativesYieldSurface(Vector& StressVector, ConstitutiveLawVariables& rVariables, ElastoPlasticConstitutiveLawVariables& rEPlasticVariables, Parameters& rValues);
 
    virtual void StressVectorInstersectionYieldSurfaces(Vector& rStressVector, const double ts, const double ts_intersection, const double ft);
 
    virtual void ConstitutiveMatrixInstersectionYieldSurfaces(Vector& StressVector, Matrix& rConstitutiveMatrix, ConstitutiveLawVariables& rVariables);
 
    virtual void GetElasticConstitutiveMatrix(Matrix& rElasticConstitutiveMatrix, 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 ElastoPlasticMohrCoulombCohesive3DLaw
}  // namespace Kratos.
#endif // KRATOS_ELASTOPLASTIC_MOHRCOULOMB_COHESIVE_3D_LAW_H_INCLUDED  defined