elastic_isotropic_3d.h

Go to the documentation of this file.

// KRATOS  ___|  |                   |                   |
//       \___ \  __|  __| |   |  __| __| |   |  __| _` | |
//             | |   |    |   | (    |   |   | |   (   | |
//       _____/ \__|_|   \__,_|\___|\__|\__,_|_|  \__,_|_| MECHANICS
//
//  License:         BSD License
//                   license: StructuralMechanicsApplication/license.txt
//
//  Main authors:    Riccardo Rossi
//                   Vicente Mataix Ferrandiz
//
 
#pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "includes/constitutive_law.h"
 
namespace Kratos
{
 
 
 
 
 
class KRATOS_API(STRUCTURAL_MECHANICS_APPLICATION) ElasticIsotropic3D
    : public ConstitutiveLaw
{
public:
 
 
    typedef ProcessInfo      ProcessInfoType;
 
    typedef ConstitutiveLaw         BaseType;
 
    typedef std::size_t             SizeType;
 
    static constexpr SizeType Dimension = 3;
 
    static constexpr SizeType VoigtSize = 6;
 
    // Adding the respective using to avoid overload conflicts
    using BaseType::Has;
    using BaseType::GetValue;
 
    KRATOS_CLASS_POINTER_DEFINITION( ElasticIsotropic3D );
 
 
    ElasticIsotropic3D();
 
    ConstitutiveLaw::Pointer Clone() const override;
 
    ElasticIsotropic3D (const ElasticIsotropic3D& rOther);
 
    ~ElasticIsotropic3D() 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;
 
    bool RequiresInitializeMaterialResponse() override
    {
        return false;
    }
 
    bool RequiresFinalizeMaterialResponse() override
    {
        return false;
    }
 
    double& CalculateValue(
        ConstitutiveLaw::Parameters& rParameterValues,
        const Variable<double>& rThisVariable,
        double& rValue
        ) override;
 
    Vector& CalculateValue(
        ConstitutiveLaw::Parameters& rParameterValues,
        const Variable<Vector>& rThisVariable,
        Vector& rValue
        ) override;
 
    // /**
    //  * @brief It calculates the value of a specified variable (StrainVectorType case)
    //  * @param rParameterValues the needed parameters for the CL calculation
    //  * @param rThisVariable the variable to be returned
    //  * @param rValue a reference to the returned value
    //  * @return rValue output: the value of the specified variable
    //  */
    // ConstitutiveLaw::StrainVectorType& CalculateValue(
    //     ConstitutiveLaw::Parameters& rParameterValues,
    //     const Variable<StrainVectorType>& rThisVariable,
    //     ConstitutiveLaw::StrainVectorType& rValue
    //     ) override;
 
    // /**
    //  * @brief It calculates the value of a specified variable (StressVectorType case)
    //  * @param rParameterValues the needed parameters for the CL calculation
    //  * @param rThisVariable the variable to be returned
    //  * @param rValue a reference to the returned value
    //  * @return rValue output: the value of the specified variable
    //  */
    // ConstitutiveLaw::StressVectorType& CalculateValue(
    //     ConstitutiveLaw::Parameters& rParameterValues,
    //     const Variable<StressVectorType>& rThisVariable,
    //     ConstitutiveLaw::StressVectorType& rValue
    //     ) override;
 
    Matrix& CalculateValue(
        ConstitutiveLaw::Parameters& rParameterValues,
        const Variable<Matrix>& rThisVariable,
        Matrix& rValue
        ) override;
 
    // /**
    //  * @brief It calculates the value of a specified variable (VoigtSizeMatrixType case)
    //  * @param rParameterValues the needed parameters for the CL calculation
    //  * @param rThisVariable the variable to be returned
    //  * @param rValue a reference to the returned value
    //  * @return rValue output: the value of the specified variable
    //  */
    // ConstitutiveLaw::VoigtSizeMatrixType& CalculateValue(
    //     ConstitutiveLaw::Parameters& rParameterValues,
    //     const Variable<VoigtSizeMatrixType>& rThisVariable,
    //     ConstitutiveLaw::VoigtSizeMatrixType& rValue
    //     ) override;
 
    // /**
    //  * @brief It calculates the value of a specified variable (DeformationGradientMatrixType case)
    //  * @param rParameterValues the needed parameters for the CL calculation
    //  * @param rThisVariable the variable to be returned
    //  * @param rValue a reference to the returned value
    //  * @return rValue output: the value of the specified variable
    //  */
    // ConstitutiveLaw::DeformationGradientMatrixType& CalculateValue(
    //     ConstitutiveLaw::Parameters& rParameterValues,
    //     const Variable<DeformationGradientMatrixType>& rThisVariable,
    //     ConstitutiveLaw::DeformationGradientMatrixType& rValue
    //     ) override;
 
    int Check(
        const Properties& rMaterialProperties,
        const GeometryType& rElementGeometry,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
 
    std::string Info() const override;
 
    void PrintInfo(std::ostream& rOStream) const override;
 
protected:
 
 
 
 
 
    void CheckClearElasticMatrix(VoigtSizeMatrixType& rConstitutiveMatrix);
 
    virtual void CalculateElasticMatrix(
        ConstitutiveLaw::VoigtSizeMatrixType& rConstitutiveMatrix,
        ConstitutiveLaw::Parameters& rValues
        );
 
    virtual void CalculatePK2Stress(
        const ConstitutiveLaw::StrainVectorType& rStrainVector,
        ConstitutiveLaw::StressVectorType& rStressVector,
        ConstitutiveLaw::Parameters& rValues
        );
 
    virtual void CalculateCauchyGreenStrain(
        ConstitutiveLaw::Parameters& rValues,
        ConstitutiveLaw::StrainVectorType& 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 ElasticIsotropic3D
}  // namespace Kratos.