plasticity_isotropic_kinematic_j2.h

Go to the documentation of this file.

// KRATOS ___                _   _ _         _   _             __                       _
//       / __\___  _ __  ___| |_(_) |_ _   _| |_(_)_   _____  / /  __ ___      _____   /_\  _ __  _ __
//      / /  / _ \| '_ \/ __| __| | __| | | | __| \ \ / / _ \/ /  / _` \ \ /\ / / __| //_\\| '_ \| '_  |
//     / /__| (_) | | | \__ \ |_| | |_| |_| | |_| |\ V /  __/ /__| (_| |\ V  V /\__ \/  _  \ |_) | |_) |
//     \____/\___/|_| |_|___/\__|_|\__|\__,_|\__|_| \_/ \___\____/\__,_| \_/\_/ |___/\_/ \_/ .__/| .__/
//                                                                                         |_|   |_|
//
//  License:         BSD License
//                   license: structural_mechanics_application/license.txt
//
//  Main authors:    Fernando Rastellini
//  Collaborators:   Riccardo Rossi
//                   Vicente Mataix Ferrandiz
 
#pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "includes/checks.h"
#include "includes/constitutive_law.h"
 
namespace Kratos
{
 
 
 
 
 
class KRATOS_API(CONSTITUTIVE_LAWS_APPLICATION) PlasticityIsotropicKinematicJ2
    : public ConstitutiveLaw
{
public:
 
 
    typedef ProcessInfo      ProcessInfoType;
    typedef ConstitutiveLaw         BaseType;
    typedef std::size_t             SizeType;
 
    static constexpr SizeType VoigtSize = 6;   //3D only (for the moment)
 
    static constexpr SizeType Dimension = VoigtSize == 6 ? 3 : 2;
 
    typedef array_1d<double, VoigtSize> BoundedArrayType;
 
    typedef Matrix MatrixType;
 
    // Counted pointer of PlasticityIsotropicKinematicJ2
    KRATOS_CLASS_POINTER_DEFINITION(PlasticityIsotropicKinematicJ2);
 
 
    PlasticityIsotropicKinematicJ2();
 
    PlasticityIsotropicKinematicJ2(const PlasticityIsotropicKinematicJ2& rOther);
 
    ~PlasticityIsotropicKinematicJ2() override;
 
    ConstitutiveLaw::Pointer Clone() const override;
 
 
 
    void GetLawFeatures(Features& rFeatures) override;
 
    SizeType WorkingSpaceDimension() override
    {
        return Dimension;
    };
 
    SizeType GetStrainSize() const override
    {
        return VoigtSize;
    };
 
    bool Has(const Variable<double>& rThisVariable) override;
 
    double& GetValue(
        const Variable<double>& rThisVariable,
        double& rValue
        ) override;
 
    void SetValue(
        const Variable<double>& rThisVariable,
        const double& rValue,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
    void InitializeMaterial(const Properties& rMaterialProperties,
                            const GeometryType& rElementGeometry,
                            const Vector& rShapeFunctionsValues) override;
 
    void FinalizeSolutionStep(const Properties& rMaterialProperties,
                            const GeometryType& rElementGeometry,
                            const Vector& rShapeFunctionsValues,
                            const ProcessInfo& rCurrentProcessInfo) 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 CalculateMaterialResponsePK1(ConstitutiveLaw::Parameters& rValues) override;
 
    void CalculateMaterialResponsePK2(ConstitutiveLaw::Parameters& rValues) override;
 
    void CalculateMaterialResponseKirchhoff(ConstitutiveLaw::Parameters& rValues) override;
 
    void CalculateMaterialResponseCauchy(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;
 
    int Check(
        const Properties& rMaterialProperties,
        const GeometryType& rElementGeometry,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
 
    void PrintData(std::ostream& rOStream) const override {
        rOStream << "Linear J2 Plasticity 3D constitutive law\n";
    };
 
protected:
 
 
 
    BoundedArrayType mPlasticStrain;           
    double           mEquivalentPlasticStrain; 
 
 
 
    double YieldFunction(
        const double      NormDeviatoricStress,
        const Properties& rMaterialProperties
        );
 
    void CalculateResponse6(
        ConstitutiveLaw::Parameters& rValues,
        BoundedArrayType&            rPlasticStrain,
        double&                      rAccumulatedPlasticStrain
        );
 
    void CalculateTangentTensor6(
        const double            DeltaGamma,
        const double            NormStressTrial,
        const BoundedArrayType& rYieldFunctionNormalVector,
        const Properties&       rMaterialProperties,
        MatrixType&             rTangentTensor
        );
 
    void CalculateElasticMatrix6(
        const Properties& rMaterialProperties,
        MatrixType&       rElasticityTensor
        );
 
 
 
 
 
 
private:
 
 
 
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
}; // Class PlasticityIsotropicKinematicJ2
} // namespace Kratos.