hencky_plastic_UP_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
//
 
 
#if !defined (KRATOS_HENCKY_PLASTIC_UP_3D_LAW_H_INCLUDED)
#define       KRATOS_HENCKY_PLASTIC_UP_3D_LAW_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_constitutive/hencky_plastic_3D_law.hpp"
#include "includes/ublas_interface.h"
 
namespace Kratos
{
 
class KRATOS_API(PARTICLE_MECHANICS_APPLICATION) HenckyElasticPlasticUP3DLaw : public HenckyElasticPlastic3DLaw
{
//protected:
 
    //struct MatrixSplit
    //{
    //Matrix  EigenValues;
    //Matrix  EigenVectors;
    //};
 
 
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(HenckyElasticPlasticUP3DLaw);
 
    HenckyElasticPlasticUP3DLaw();
 
 
    HenckyElasticPlasticUP3DLaw(MPMFlowRulePointer pMPMFlowRule, YieldCriterionPointer pYieldCriterion, HardeningLawPointer pHardeningLaw);
 
    HenckyElasticPlasticUP3DLaw (const HenckyElasticPlasticUP3DLaw& rOther)  ;
 
 
    //HenckyElasticPlastic3DLaw& operator=(const HenckyElasticPlastic3DLaw& rOther);
 
    ConstitutiveLaw::Pointer Clone() const override;
 
    ~HenckyElasticPlasticUP3DLaw() override;
 
    SizeType WorkingSpaceDimension() override
    {
        return 3;
    };
 
    SizeType GetStrainSize() const override
    {
        return 6;
    };
 
    void GetLawFeatures(Features& rFeatures) override;
 
    //virtual String Info() const;
    //virtual void PrintInfo(std::ostream& rOStream) const;
    //virtual void PrintData(std::ostream& rOStream) const;
 
protected:
 
 
 
 
 
    void CorrectDomainPressure( Matrix& rStressMatrix, const MaterialResponseVariables& rElasticVariables) override;
 
    void GetDomainPressure( double& rPressure, const MaterialResponseVariables& rElasticVariables);
 
    void CalculateElastoPlasticTangentMatrix( const ParticleFlowRule::RadialReturnVariables & rReturnMappingVariables, const Matrix& rNewElasticLeftCauchyGreen,const double& rAlpha, Matrix& rElastoPlasticTangentMatrix, const MaterialResponseVariables& rElasticVariables, const Properties& rProperties) override;
 
    void CalculateGreenLagrangeStrain( const Matrix & rRightCauchyGreen,
                                       Vector& rStrainVector ) override;
 
 
    void CalculateAlmansiStrain( const Matrix & rLeftCauchyGreen,
                                 Vector& rStrainVector ) override;
 
    void CalculatePrincipalStressTrial(const MaterialResponseVariables & rElasticVariables,Parameters & rValues,
                                       const ParticleFlowRule::RadialReturnVariables& rReturnMappingVariables,
                                       Matrix& rNewElasticLeftCauchyGreen, Matrix& rStressMatrix) override;
 
private:
 
 
 
 
 
 
 
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, HenckyElasticPlastic3DLaw )
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, HenckyElasticPlastic3DLaw )
    }
 
 
 
 
}; // Class HenckyElasticPlasticUP3DLaw
 
} //namespace Kratos
 
#endif  //KRATOS_HENCKY_PLASTIC_3D_LAW_H_INCLUDED