elastic_isotropic_3d.h

Go to the documentation of this file.

// KRATOS  ___|  |                   |                   |
//       \___ \  __|  __| |   |  __| __| |   |  __| _` | |
//             | |   |    |   | (    |   |   | |   (   | |
//       _____/ \__|_|   \__,_|\___|\__|\__,_|_|  \__,_|_| MECHANICS
//
//  License:         BSD License
//                   license: structural_mechanics_application/license.txt
//
//  Main authors:    Riccardo Rossi
//                   Vicente Mataix Ferrandiz
//
 
#if !defined (KRATOS_ELASTIC_ISOTROPIC_3D_LAW_H_INCLUDED)
#define  KRATOS_ELASTIC_ISOTROPIC_3D_LAW_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "includes/constitutive_law.h"
 
namespace Kratos
{
 
 
 
 
 
class KRATOS_API(FEM_TO_DEM_APPLICATION) ElasticIsotropic3DFEMDEM
    : public ConstitutiveLaw
{
public:
 
 
    typedef ProcessInfo      ProcessInfoType;
 
    typedef ConstitutiveLaw         BaseType;
 
    typedef std::size_t             SizeType;
 
    static constexpr SizeType Dimension = 3;
 
    static constexpr SizeType VoigtSize = 6;
 
    KRATOS_CLASS_POINTER_DEFINITION( ElasticIsotropic3DFEMDEM );
 
 
    ElasticIsotropic3DFEMDEM();
 
    ConstitutiveLaw::Pointer Clone() const override;
 
    ElasticIsotropic3DFEMDEM (const ElasticIsotropic3DFEMDEM& rOther);
 
    ~ElasticIsotropic3DFEMDEM() override;
 
 
 
    void GetLawFeatures(Features& rFeatures) override;
 
    SizeType WorkingSpaceDimension() override
    {
        return Dimension;
    };
 
    SizeType GetStrainSize() const override
    {
        return VoigtSize;
    };
 
    StrainMeasure GetStrainMeasure() override
    {
        return StrainMeasure_Infinitesimal;
    }
 
    StressMeasure GetStressMeasure() override
    {
        return StressMeasure_Cauchy;
    }
 
    void CalculateMaterialResponsePK1 (ConstitutiveLaw::Parameters & rValues) override;
 
    void CalculateMaterialResponsePK2 (ConstitutiveLaw::Parameters & rValues) override;
 
    void CalculateMaterialResponseKirchhoff (ConstitutiveLaw::Parameters & rValues) override;
 
    void CalculateMaterialResponseCauchy (ConstitutiveLaw::Parameters & rValues) override;
 
    void InitializeMaterialResponsePK1 (ConstitutiveLaw::Parameters& rValues) override;
 
    void InitializeMaterialResponsePK2 (ConstitutiveLaw::Parameters& rValues) override;
 
    void InitializeMaterialResponseKirchhoff (ConstitutiveLaw::Parameters& rValues) override;
 
    void InitializeMaterialResponseCauchy (ConstitutiveLaw::Parameters& rValues) override;
 
    void FinalizeMaterialResponsePK1 (ConstitutiveLaw::Parameters & rValues) override;
 
    void FinalizeMaterialResponsePK2 (ConstitutiveLaw::Parameters & rValues) override;
 
    void FinalizeMaterialResponseKirchhoff (ConstitutiveLaw::Parameters & rValues)  override;
 
    void FinalizeMaterialResponseCauchy (ConstitutiveLaw::Parameters & rValues) override;
 
    double& CalculateValue(
        ConstitutiveLaw::Parameters& rParameterValues,
        const Variable<double>& rThisVariable,
        double& rValue
        ) override;
 
    Vector& CalculateValue(
        ConstitutiveLaw::Parameters& rParameterValues,
        const Variable<Vector>& rThisVariable,
        Vector& rValue
        ) override;
 
    Matrix& CalculateValue(
        ConstitutiveLaw::Parameters& rParameterValues,
        const Variable<Matrix>& rThisVariable,
        Matrix& rValue
        ) override;
 
    int Check(
        const Properties& rMaterialProperties,
        const GeometryType& rElementGeometry,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
protected:
 
 
 
 
 
    void CheckClearElasticMatrix(Matrix& rConstitutiveMatrix);
 
    virtual void CalculateElasticMatrix(
        Matrix& rConstitutiveMatrix,
        ConstitutiveLaw::Parameters& rValues
        );
 
    virtual void CalculatePK2Stress(
        const Vector& rStrainVector,
        Vector& rStressVector,
        ConstitutiveLaw::Parameters& rValues
        );
 
    virtual void CalculateCauchyGreenStrain(
        ConstitutiveLaw::Parameters& rValues,
        Vector& rStrainVector
        );
 
 
private:
 
 
 
 
 
 
    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 ElasticIsotropic3DFEMDEM
}  // namespace Kratos.
#endif // KRATOS_ELASTIC_ISOTROPIC_3D_LAW_H_INCLUDED  defined