user_provided_linear_elastic_law.h

Go to the documentation of this file.

// KRATOS  ___|  |                   |                   |
//       \___ \  __|  __| |   |  __| __| |   |  __| _` | |
//             | |   |    |   | (    |   |   | |   (   | |
//       _____/ \__|_|   \__,_|\___|\__|\__,_|_|  \__,_|_| MECHANICS
//
//  License:         BSD License
//                   license: StructuralMechanicsApplication/license.txt
//
//  Main authors:    Ruben Zorrilla
//                   Riccardo Rossi
//
 
#pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "includes/constitutive_law.h"
 
namespace Kratos
{
 
 
 
 
 
 
template <unsigned int TDim>
class KRATOS_API(STRUCTURAL_MECHANICS_APPLICATION) UserProvidedLinearElasticLaw
    : public ConstitutiveLaw
{
public:
 
 
    typedef ProcessInfo      ProcessInfoType;
 
    typedef ConstitutiveLaw         BaseType;
 
    typedef std::size_t             SizeType;
 
    static constexpr SizeType Dimension = TDim;
 
    static constexpr SizeType StrainSize = (TDim * 3) - 3;
 
    KRATOS_CLASS_POINTER_DEFINITION( UserProvidedLinearElasticLaw );
 
 
    UserProvidedLinearElasticLaw();
 
    ConstitutiveLaw::Pointer Clone() const override;
 
    UserProvidedLinearElasticLaw (const UserProvidedLinearElasticLaw& rOther);
 
    ~UserProvidedLinearElasticLaw();
 
 
 
    void GetLawFeatures(Features& rFeatures) override;
 
    SizeType WorkingSpaceDimension() override
    {
        return Dimension;
    };
 
    SizeType GetStrainSize() const override
    {
        return StrainSize;
    };
 
    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 (Parameters& rValues) override {};
 
    void InitializeMaterialResponsePK2 (Parameters& rValues) override {};
 
    void InitializeMaterialResponseKirchhoff (Parameters& rValues) override {};
 
    void InitializeMaterialResponseCauchy (Parameters& rValues) override {};
 
    void FinalizeMaterialResponsePK1 (Parameters& rValues) override {};
 
    void FinalizeMaterialResponsePK2 (Parameters& rValues) override {};
 
    void FinalizeMaterialResponseKirchhoff (Parameters& rValues) override {};
 
    void FinalizeMaterialResponseCauchy (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;
 
    // /**
    //  * @brief It calculates the value of a specified variable (Vector 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 (Vector 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 (Matrix 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 (Matrix 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;
 
protected:
 
 
 
 
 
    virtual void CalculateElasticMatrix(
        ConstitutiveLaw::VoigtSizeMatrixType& rConstitutiveMatrix,
        ConstitutiveLaw::Parameters& rValues
        );
 
    virtual void CalculatePK2Stress(
        const ConstitutiveLaw::StrainVectorType& rStrainVector,
        ConstitutiveLaw::StressVectorType& rStressVector,
        ConstitutiveLaw::Parameters& rValues
        );
 
    virtual void CalculateGreenLagrangeStrainVector(
        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 UserProvidedLinearElasticLaw
}  // namespace Kratos.