linear_elastic_plane_strain_2D_law.h

Go to the documentation of this file.

// KRATOS___
//     //   ) )
//    //         ___      ___
//   //  ____  //___) ) //   ) )
//  //    / / //       //   / /
// ((____/ / ((____   ((___/ /  MECHANICS
//
//  License:         geo_mechanics_application/license.txt
//
//  Main authors:    Vahid Galavi
//
 
#pragma once
 
// Project includes
#include "custom_constitutive/linear_elastic_plane_strain_K0_law.h"
 
namespace Kratos
{
 
class KRATOS_API(GEO_MECHANICS_APPLICATION) GeoLinearElasticPlaneStrain2DLaw
    : public LinearPlaneStrainK0Law
{
public:
    using BaseType = LinearPlaneStrainK0Law;
 
    using SizeType = std::size_t;
 
    static constexpr SizeType Dimension = N_DIM_2D;
 
    // for the time being
    static constexpr SizeType VoigtSize = VOIGT_SIZE_2D_PLANE_STRAIN;
 
    KRATOS_CLASS_POINTER_DEFINITION( GeoLinearElasticPlaneStrain2DLaw );
 
    ConstitutiveLaw::Pointer Clone() const override;
 
    bool RequiresInitializeMaterialResponse() override;
    void InitializeMaterialResponseCauchy(ConstitutiveLaw::Parameters& rValues) override;
 
    bool RequiresFinalizeMaterialResponse() override;
    void FinalizeMaterialResponseCauchy(ConstitutiveLaw::Parameters & rValues) override;
    void FinalizeMaterialResponsePK2(ConstitutiveLaw::Parameters& rValues) override;
 
    void GetLawFeatures(Features& rFeatures) override;
 
    SizeType WorkingSpaceDimension() override;
 
    SizeType GetStrainSize() const override;
 
    bool IsIncremental() override;
 
    bool& GetValue(const Variable<bool>& rThisVariable, bool& rValue) override;
    using ConstitutiveLaw::GetValue;
 
protected:
    void CalculateElasticMatrix(Matrix& C, ConstitutiveLaw::Parameters& rValues) override;
 
    void CalculatePK2Stress(const Vector& rStrainVector,
                            Vector& rStressVector,
                            ConstitutiveLaw::Parameters& rValues) override;
 
 
 
 
private:
    Vector mStressVector          = ZeroVector(VoigtSize);
    Vector mStressVectorFinalized = ZeroVector(VoigtSize);
    Vector mDeltaStrainVector     = ZeroVector(VoigtSize);
    Vector mStrainVectorFinalized = ZeroVector(VoigtSize);
    bool   mIsModelInitialized    = false;
 
    friend class Serializer;
    void save(Serializer& rSerializer) const override;
    void load(Serializer& rSerializer) override;
}; // Class GeoLinearElasticPlaneStrain2DLaw
 
}