truss_plasticity_constitutive_law.h

Go to the documentation of this file.

// KRATOS ___                _   _ _         _   _             __                       _
//       / __\___  _ __  ___| |_(_) |_ _   _| |_(_)_   _____  / /  __ ___      _____   /_\  _ __  _ __
//      / /  / _ \| '_ \/ __| __| | __| | | | __| \ \ / / _ \/ /  / _` \ \ /\ / / __| //_\\| '_ \| '_  |
//     / /__| (_) | | | \__ \ |_| | |_| |_| | |_| |\ V /  __/ /__| (_| |\ V  V /\__ \/  _  \ |_) | |_) |
//     \____/\___/|_| |_|___/\__|_|\__|\__,_|\__|_| \_/ \___\____/\__,_| \_/\_/ |___/\_/ \_/ .__/| .__/
//                                                                                         |_|   |_|
//
//  License:         BSD License
//                   license: structural_mechanics_application/license.txt
//
//  Main authors:    Klaus B. Sautter
//
 
#pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "includes/constitutive_law.h"
#include "includes/checks.h"
 
namespace Kratos
{
 
class KRATOS_API(CONSTITUTIVE_LAWS_APPLICATION) TrussPlasticityConstitutiveLaw : public ConstitutiveLaw
{
public:
    typedef ProcessInfo      ProcessInfoType;
    typedef ConstitutiveLaw         BaseType;
    typedef std::size_t             SizeType;
    KRATOS_CLASS_POINTER_DEFINITION( TrussPlasticityConstitutiveLaw );
 
    TrussPlasticityConstitutiveLaw();
 
    ConstitutiveLaw::Pointer Clone() const override;
 
    TrussPlasticityConstitutiveLaw (const TrussPlasticityConstitutiveLaw& rOther);
 
 
    ~TrussPlasticityConstitutiveLaw() override;
 
    void GetLawFeatures(Features& rFeatures) override;
 
    void SetValue(const Variable<double>& rThisVariable,
                          const double& rValue,
                          const ProcessInfo& rCurrentProcessInfo) override;
 
    double& GetValue(const Variable<double>& rThisVariable,
                    double& rValue) override;
 
    array_1d<double, 3 > & GetValue(const Variable<array_1d<double, 3 > >& rThisVariable,
        array_1d<double, 3 > & rValue) override;
 
    double& CalculateValue(ConstitutiveLaw::Parameters& rParameterValues,
                    const Variable<double>& rThisVariable,
                    double& rValue) override;
 
    double TrialYieldFunction(const Properties& rMaterialProperties,
     const double& rCurrentStress);
 
    bool CheckIfIsPlasticRegime(Parameters& rValues,const double& rCurrentStress);
 
    void FinalizeMaterialResponsePK2(Parameters& rValues) override;
 
    void CalculateMaterialResponsePK2(Parameters& rValues) override;
 
    void CalculateMaterialResponsePK2Custom(Parameters& rValues, double& rCurrentAccumulatedPlasticStrain, double& rCurrentPlasticAlpha);
 
    bool RequiresInitializeMaterialResponse() override
    {
        return false;
    }
 
    SizeType GetStrainSize() const override
    {
        return 1;
    }
 
    int Check(
        const Properties& rMaterialProperties,
        const GeometryType& rElementGeometry,
        const ProcessInfo& rCurrentProcessInfo
    ) const override;
 
protected:
 
 
 
 
 
private:
 
 
    //the following members are
    bool mCurrentInElasticFlag = false;
    double mPlasticAlpha = 0.0;
    double mAccumulatedPlasticStrain = 0.0;
 
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, ConstitutiveLaw);
        rSerializer.save("PlasticAlpha", mPlasticAlpha);
        rSerializer.save("AccumulatedPlasticStrain", mAccumulatedPlasticStrain);
        rSerializer.save("CurrentInElasticFlag", mCurrentInElasticFlag);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, ConstitutiveLaw);
        rSerializer.load("PlasticAlpha", mPlasticAlpha);
        rSerializer.load("AccumulatedPlasticStrain", mAccumulatedPlasticStrain);
        rSerializer.load("CurrentInElasticFlag", mCurrentInElasticFlag);
    }
 
 
}; // Class TrussPlasticityConstitutiveLaw
}  // namespace Kratos.