TPZMatElastoPlastic.h

Go to the documentation of this file.

#ifndef PZELASTOPLASTIC_H
#define PZELASTOPLASTIC_H


#include "TPZMaterial.h"
#include "TPZMatWithMem.h"
#include "TPZElastoPlasticMem.h"
#include "pzporoelastoplasticmem.h"
#include "TPZPorousElasticResponse.h"

template <class T, class TMEM = TPZElastoPlasticMem>
class  TPZMatElastoPlastic : public TPZMatWithMem<TMEM>
{

public:
    
    TPZMatElastoPlastic();

    TPZMatElastoPlastic(int id);

    TPZMatElastoPlastic(const TPZMatElastoPlastic &other);

    virtual ~TPZMatElastoPlastic();

    virtual void SetPlasticityModel(T & plasticity);

    virtual void UpdateMaterialCoeficients(TPZVec<REAL> &x,T & plasticity);

    virtual void SetBulkDensity(REAL & RhoB);

    virtual std::string Name() override;

    virtual int Dimension() const override { return 3; }

    virtual int NStateVariables() const override{ return 3; }

    virtual void Print(std::ostream &out, const int memory);

    virtual void Print(std::ostream &out) override;

    virtual int VariableIndex(const std::string &name) override;

    virtual int NSolutionVariables(int var) override;

    virtual void Solution(TPZMaterialData &data, int var, TPZVec<REAL> &Solout) override;

    virtual int NFluxes() override
    {
         PZError << "TPZMatElastoPlastic::NFluxes() - Method not implemented\n";
         return 0;
    }

    virtual void Flux(TPZVec<REAL> &x, TPZVec<REAL> &Sol, TPZFMatrix<REAL> &DSol, TPZFMatrix<REAL> &axes, TPZVec<REAL> &flux) override
    {
        PZError << "TPZMatElastoPlastic::Flux - Method not implemented\n";
    }

    virtual void Errors(TPZVec<REAL> &x,TPZVec<REAL> &u, TPZFMatrix<REAL> &dudx,
                      TPZFMatrix<REAL> &axes, TPZVec<REAL> &flux,
                      TPZVec<REAL> &u_exact,TPZFMatrix<REAL> &du_exact,TPZVec<REAL> &values) override;
    virtual int NEvalErrors()  override {return 3;}

    virtual void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<REAL> &ek, TPZFMatrix<REAL> &ef) override;

    virtual void ContributeBC(TPZMaterialData &data, REAL weight, TPZFMatrix<REAL> &ek, TPZFMatrix<REAL> &ef, TPZBndCond &bc) override;

    virtual void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<REAL> &ef) override;

    virtual void ContributeBC(TPZMaterialData &data, REAL weight, TPZFMatrix<REAL> &ef, TPZBndCond &bc) override;

    void ComputeStrainVector(TPZMaterialData & data, TPZFMatrix<REAL> &Strain);

    void ComputeDeltaStrainVector(TPZMaterialData & data, TPZFMatrix<REAL> &DeltaStrain);

    void ComputeStressVector(TPZMaterialData & data, TPZFMatrix<REAL> &Stress);

    void CheckConvergence(TPZMaterialData & data,TPZFMatrix<REAL> & DeltaStrain);

    void ApplyDeltaStrainComputeDep(TPZMaterialData & data, TPZFMatrix<REAL> & DeltaStrain,
                                            TPZFMatrix<REAL> & Stress, TPZFMatrix<REAL> & Dep);

    void ApplyDeltaStrain(TPZMaterialData & data, TPZFMatrix<REAL> & DeltaStrain,
                                            TPZFMatrix<REAL> & Stress);

    void ApplyDirection(TPZFMatrix<REAL> &vectorTensor, TPZVec<REAL> &Out);

    void vectorToTensor(const TPZFMatrix<REAL> & vectorTensor, TPZFMatrix<REAL> & Tensor);

    void EigenValues(TPZFMatrix<REAL> & vectorTensor, TPZVec<REAL> & ev);

    void EigenVectors(TPZFMatrix<REAL> &vectorTensor, TPZVec< REAL > &Solout, int direction);

    virtual TPZMaterial * NewMaterial() override;

    virtual int ClassId() const override;

    virtual void Write(TPZStream &buf, int withclassid) const override;

    virtual void Read(TPZStream &buf, void *context) override;

    virtual void SetTol(const REAL & tol);

    virtual void SetBulkDensity(const REAL & bulk);
    
    virtual void SetPorousElasticity(TPZPorousElasticResponse & PER);
    
    void SetPlasticModel(T & plasticity_model);
    
    virtual T & GetPlasticModel();
    
    virtual TPZPorousElasticResponse & GetPorousElasticity(TPZPorousElasticResponse & PER);

    virtual void FillDataRequirements(TPZMaterialData &data) override;

    virtual void FillBoundaryConditionDataRequirement(int type,TPZMaterialData &data) override;

    enum ESolutionVar {
        ENone = -1,
        EDisplacementDoF  = 0,
        EDisplacement     = 1,
        EStrain           = 2,
        EStress           = 3,
        EStrainElastic    = 4,
        EStrainPlastic    = 5,
        EYield            = 6,
        EVolHardening     = 7,
        EStrainPValues    = 8,
        EStressPValues    = 9,
        EStrainElasticPValues    = 10,
        EStrainPlasticPValues    = 11,
        EStrainI1           = 12,
        EStressI1           = 13,
        EStrainElasticI1    = 14,
        EStrainPlasticI1    = 15,
        EStrainJ2           = 16,
        EStressJ2           = 17,
        EStrainElasticJ2    = 18,
        EStrainPlasticJ2    = 19,
        EFailureType    = 20,
    };

protected:

    TPZManVector<REAL, 3> m_force;

    REAL m_rho_bulk;

    TPZManVector<REAL,3> m_PostProcessDirection;

    T m_plasticity_model;

    REAL m_tol;

    bool m_use_non_linear_elasticity_Q;
    
    TPZPorousElasticResponse m_PER;
                
};

template <class T, class TMEM>
int TPZMatElastoPlastic<T,TMEM>::ClassId() const{
    return Hash("TPZMatElastoPlastic") ^ TPZMatWithMem<TMEM>::ClassId() << 1 ^ T().ClassId() << 2;
}

#endif