hencky_plastic_3D_law.hpp

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:        BSD License
//                  Kratos default license: kratos/license.txt
//
//  Main authors:    Ilaria Iaconeta, Bodhinanda Chandra
//
 
 
#if !defined (KRATOS_HENCKY_PLASTIC_3D_LAW_H_INCLUDED)
#define       KRATOS_HENCKY_PLASTIC_3D_LAW_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_constitutive/hyperelastic_3D_law.hpp"
#include "custom_constitutive/flow_rules/particle_flow_rule.hpp"
#include "includes/ublas_interface.h"
 
namespace Kratos
{
 
class KRATOS_API(PARTICLE_MECHANICS_APPLICATION) HenckyElasticPlastic3DLaw : public HyperElastic3DLaw
{
protected:
 
    struct VectorSplit
    {
        Vector  Isochoric;
        Vector  Volumetric;
    };
 
    struct MatrixSplit
    {
        Matrix  EigenValues;
        Matrix  EigenVectors;
        Matrix  Isochoric;
        Matrix  Volumetric;
        Matrix  Plastic;
    };
 
    struct PlasticMaterialResponseVariables
    {
        Matrix  TrialLeftStretchTensor;
        Matrix  InverseTrialLeftStretchTensor;
    };
 
 
public:
    typedef ProcessInfo      ProcessInfoType;
    typedef ConstitutiveLaw         BaseType;
    typedef std::size_t             SizeType;
 
    typedef ParticleFlowRule::Pointer                MPMFlowRulePointer;
    typedef ParticleYieldCriterion::Pointer    YieldCriterionPointer;
    typedef ParticleHardeningLaw::Pointer        HardeningLawPointer;
    typedef Properties::Pointer            PropertiesPointer;
 
    KRATOS_CLASS_POINTER_DEFINITION(HenckyElasticPlastic3DLaw);
 
    HenckyElasticPlastic3DLaw();
 
 
    HenckyElasticPlastic3DLaw(MPMFlowRulePointer pMPMFlowRule, YieldCriterionPointer pYieldCriterion, HardeningLawPointer pHardeningLaw);
 
    HenckyElasticPlastic3DLaw (const HenckyElasticPlastic3DLaw& rOther)  ;
 
 
    //HyperElasticPlastic3DLaw& operator=(const HyperElasticPlastic3DLaw& rOther);
 
    ConstitutiveLaw::Pointer Clone() const override;
 
    ~HenckyElasticPlastic3DLaw() override;
 
    SizeType WorkingSpaceDimension() override
    {
        return 3;
    };
 
    SizeType GetStrainSize() const override
    {
        return 6;
    };
 
 
    void GetLawFeatures(Features& rFeatures) override;
 
 
    bool Has( const Variable<double>& rThisVariable ) override;
    bool Has( const Variable<Vector>& rThisVariable ) override;
    bool Has( const Variable<Matrix>& rThisVariable ) override;
 
    double& GetValue( const Variable<double>& rThisVariable, double& rValue ) override;
    Vector& GetValue( const Variable<Vector>& rThisVariable, Vector& rValue ) override;
    Matrix& GetValue( const Variable<Matrix>& rThisVariable, Matrix& rValue ) override;
 
 
    void SetValue( const Variable<double>& rVariable,
                   const double& rValue,
                   const ProcessInfo& rCurrentProcessInfo ) override;
    void SetValue( const Variable<Vector>& rThisVariable,
                   const Vector& rValue,
                   const ProcessInfo& rCurrentProcessInfo ) override;
    void SetValue( const Variable<Matrix>& rThisVariable,
                   const Matrix& rValue,
                   const ProcessInfo& rCurrentProcessInfo ) override;
 
 
    void InitializeMaterial( const Properties& rProps,
                             const GeometryType& rGeometry,
                             const Vector& rShapeFunctionsValues ) override;
 
    void CalculateMaterialResponseKirchhoff (Parameters & rValues) override;
 
 
    int Check(const Properties& rMaterialProperties, const GeometryType& rElementGeometry, const ProcessInfo& rCurrentProcessInfo) const override;
 
 
 
    //virtual String Info() const;
    //virtual void PrintInfo(std::ostream& rOStream) const;
    //virtual void PrintData(std::ostream& rOStream) const;
 
protected:
 
    unsigned int mPlasticRegion;
    Matrix mPlasticDeformationGradient;
    Matrix mElasticLeftCauchyGreen;
 
    MPMFlowRulePointer mpMPMFlowRule;
 
    YieldCriterionPointer mpYieldCriterion;
 
    HardeningLawPointer   mpHardeningLaw;
 
 
 
    void GetDomainPressure( double& rPressure, const MaterialResponseVariables& rElasticVariables);
 
    virtual Matrix SetConstitutiveMatrixToAppropiateDimension(Matrix& rConstitutiveMatrix, const Matrix& rElastoPlasticTangentMatrix);
 
    virtual Vector SetStressMatrixToAppropiateVectorDimension(Vector& rStressVector, const Matrix& rStressMatrix );
 
    virtual void CorrectDomainPressure( Matrix& rStressMatrix, const MaterialResponseVariables& rElasticVariables);
 
    virtual void CalculateElastoPlasticTangentMatrix( const ParticleFlowRule::RadialReturnVariables & rReturnMappingVariables, const Matrix& rNewElasticLeftCauchyGreen,const double& rAlpha, Matrix& rElastoPlasticMatrix, const MaterialResponseVariables& rElasticVariables, const Properties& rProperties);
 
    double& TensorComponent(double & rCabcd,
                            const Matrix& rMA, const Matrix& rMB,
                            const unsigned int& a, const unsigned int& b,
                            const unsigned int& c, const unsigned int& d);
 
    virtual void MyTensorProduct(const Matrix& rMA, const Matrix& rMB,
                                 Matrix& rEigenbasesProductMatrix);
 
    double& TensorComponent2(double & rCabcd,
                             const Matrix& rMA, const Matrix& rMB,
                             const unsigned int& a, const unsigned int& b,
                             const unsigned int& c, const unsigned int& d);
 
    virtual void MyTensorProduct2(const Matrix& rMA, const Matrix& rMB,
                                  Matrix& rEigenbasesProductMatrix);
 
    double& TensorComponent3(double & rCabcd,
                             const Matrix& rMA,
                             const unsigned int& a, const unsigned int& b,
                             const unsigned int& c, const unsigned int& d);
 
    virtual void MyTensorProduct3(const Matrix& rMA,
                                  Matrix& rEigenbasesProductMatrix);
 
    double& TensorComponent4(double & rCabcd,
                             const Matrix& rMA,
                             const unsigned int& a, const unsigned int& b,
                             const unsigned int& c, const unsigned int& d);
 
    virtual void MyTensorProduct4(const Matrix& rMA,
                                  Matrix& rEigenbasesProductMatrix);
 
 
    virtual Matrix CalculateEigenbases(const ParticleFlowRule::RadialReturnVariables& rReturnMappingVariables, Matrix& rEigenbasesMatrix);
 
 
 
    Vector& GetStressVectorFromMatrix(const Matrix& rStressMatrix,
                                      Vector& rMainStress,
                                      const Matrix& rEigenVectors);
 
    virtual void CalculateHenckyMainStrain(const Matrix& rCauchyGreeMatrix,
                                           ParticleFlowRule::RadialReturnVariables& rReturnMappingVariables,
                                           Vector& rMainStrain);
 
    virtual void CalculatePrincipalStressTrial(const MaterialResponseVariables & rElasticVariables,Parameters & rValues, const ParticleFlowRule::RadialReturnVariables& rReturnMappingVariables,
            Matrix& rNewElasticLeftCauchyGreen, Matrix& rStressMatrix);
 
 
private:
 
 
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, HyperElastic3DLaw )
 
        rSerializer.save("mElasticLeftCauchyGreen",mElasticLeftCauchyGreen);
        rSerializer.save("mpMPMFlowRule",mpMPMFlowRule);
        rSerializer.save("mpYieldCriterion",mpYieldCriterion);
        rSerializer.save("mpHardeningLaw",mpHardeningLaw);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, HyperElastic3DLaw )
 
        rSerializer.load("mElasticLeftCauchyGreen",mElasticLeftCauchyGreen);
        rSerializer.load("mpMPMFlowRule",mpMPMFlowRule);
        rSerializer.load("mpYieldCriterion",mpYieldCriterion);
        rSerializer.load("mpHardeningLaw",mpHardeningLaw);
    }
 
 
 
 
}; // Class HenckyElasticPlastic3DLaw
 
} //namespace Kratos
 
#endif  //KRATOS_HENCKY_PLASTIC_3D_LAW_H_INCLUDED