johnson_cook_thermal_plastic_3D_law.hpp

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:        BSD License
//                  Kratos default license: kratos/license.txt
//
//  Main authors:    Peter Wilson
//
 
#if !defined (KRATOS_JOHNSON_COOK_THERMAL_PLASTIC_3D_LAW_H_INCLUDED)
#define  KRATOS_JOHNSON_COOK_THERMAL_PLASTIC_3D_LAW_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
 
#include "custom_constitutive/hyperelastic_3D_law.hpp"
#include "includes/checks.h"
 
 
namespace Kratos
{
class KRATOS_API(PARTICLE_MECHANICS_APPLICATION) JohnsonCookThermalPlastic3DLaw : public HyperElastic3DLaw
{
public:
 
    typedef ProcessInfo      ProcessInfoType;
    typedef HyperElastic3DLaw         BaseType;
    typedef std::size_t             SizeType;
    typedef Properties::Pointer            PropertiesPointer;
 
    KRATOS_CLASS_POINTER_DEFINITION(JohnsonCookThermalPlastic3DLaw);
 
    JohnsonCookThermalPlastic3DLaw();
 
    JohnsonCookThermalPlastic3DLaw(const JohnsonCookThermalPlastic3DLaw& rOther);
 
    JohnsonCookThermalPlastic3DLaw& operator=(const JohnsonCookThermalPlastic3DLaw& rOther);
 
    ConstitutiveLaw::Pointer Clone() const override;
 
    ~JohnsonCookThermalPlastic3DLaw() override;
 
    void GetLawFeatures(Features& rFeatures) override;
 
    double& GetValue(const Variable<double>& rThisVariable, double& rValue) override;
 
    void SetValue(const Variable<double>& rVariable,
        const double& rValue,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    bool Has(const Variable<double>& rThisVariable) override;
 
    void InitializeMaterial(const Properties& rMaterialProperties,
        const GeometryType& rElementGeometry,
        const Vector& rShapeFunctionsValues) override;
 
    void CalculateMaterialResponseKirchhoff(Parameters& rValues) override;
 
    int Check(const Properties& rMaterialProperties, const GeometryType& rElementGeometry, const ProcessInfo& rCurrentProcessInfo) const override;//E:\Kratos\applications\SolidMechanicsApplication\custom_constitutive\hyperelastic_plastic_3D_law.hpp
 
protected:
 
 
    double mEquivalentStress;
    Vector mStrainOld;
    double mEquivalentPlasticStrainOld;
    double mPlasticStrainRateOld;
    double mTemperatureOld;
    double mGammaOld;
    double mEnergyInternal;
    double mEnergyDissipated;
    double mYieldStressOld;
    double mYieldStressVirgin;
    double mHardeningRatio;
 
 
    bool CheckParameters(Parameters& rValues) override;
 
    virtual void MakeStrainStressMatrixFromVector(const Vector& rInput, Matrix& rOutput);
 
    virtual void MakeStrainStressVectorFromMatrix(const Matrix& rInput, Vector& rOutput);
 
    void CheckIsExplicitTimeIntegration(const ProcessInfo& rCurrentProcessInfo);
 
private:
 
    double CalculateHardenedYieldStress(const Properties& MaterialProperties, const double EquivalentPlasticStrain,
        const double PlasticStrainRate, const double Temperature);
 
    double CalculateThermalHardeningFactor(const Properties& MaterialProperties, const double Temperature);
 
    double CalculateStrainRateHardeningFactor(const Properties& MaterialProperties, const double PlasticStrainRate);
 
    double CalculateThermalDerivative(const Properties& MaterialProperties, const double EquivalentPlasticStrain,
        const double PlasticStrainRate, const double Temperature);
 
    double CalculatePlasticStrainRateDerivative(const Properties& MaterialProperties, const double EquivalentPlasticStrain,
        const double PlasticStrainRate, const double Temperature);
 
    double CalculatePlasticStrainDerivative(const Properties& MaterialProperties, const double EquivalentPlasticStrain,
        const double PlasticStrainRate, const double Temperature);
 
    inline double GetSqrt32() { return 1.2247448713915900000; } //sqrt(3.0/2.0)
 
    inline double GetSqrt23() { return 0.8164965809277260000; } //sqrt(2.0/3.0)
 
    inline double GetSqrt6() { return 2.4494897427831800000; } //sqrt(6.0)
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, HyperElastic3DLaw);
 
        rSerializer.save("mEquivalentStress", mEquivalentStress);
        rSerializer.save("mStrainOld", mStrainOld);
        rSerializer.save("mEquivalentPlasticStrainOld", mEquivalentPlasticStrainOld);
        rSerializer.save("mPlasticStrainRateOld", mPlasticStrainRateOld);
        rSerializer.save("mTemperatureOld", mTemperatureOld);
        rSerializer.save("mGammaOld", mGammaOld);
        rSerializer.save("mEnergyInternal", mEnergyInternal);
        rSerializer.save("mEnergyDissipated", mEnergyDissipated);
        rSerializer.save("mYieldStressOld", mYieldStressOld);
        rSerializer.save("mYieldStressVirgin", mYieldStressVirgin);
        rSerializer.save("mHardeningRatio", mHardeningRatio);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, HyperElastic3DLaw);
 
        rSerializer.load("mEquivalentStress", mEquivalentStress);
        rSerializer.load("mStrainOld", mStrainOld);
        rSerializer.load("mEquivalentPlasticStrainOld", mEquivalentPlasticStrainOld);
        rSerializer.load("mPlasticStrainRateOld", mPlasticStrainRateOld);
        rSerializer.load("mTemperatureOld", mTemperatureOld);
        rSerializer.load("mGammaOld", mGammaOld);
        rSerializer.load("mEnergyInternal", mEnergyInternal);
        rSerializer.load("mEnergyDissipated", mEnergyDissipated);
        rSerializer.load("mYieldStressOld", mYieldStressOld);
        rSerializer.load("mYieldStressVirgin", mYieldStressVirgin);
        rSerializer.load("mHardeningRatio", mHardeningRatio);
    }
 
 
 
}; // Class JohnsonCookThermalPlastic3DLaw
}  // namespace Kratos.
#endif // KRATOS_JOHNSON_COOK_THERMAL_PLASTIC_3D_LAW_H_INCLUDED defined