TPZElasticCriterion.h

Go to the documentation of this file.

/*
 *  TPZElasticCriterion
 *
 *  Created by Nathan Shauer on 5/4/14.
 *  Copyright 2014 __LabMeC__. All rights reserved.
 *
 */

#ifndef TPZElasticResponse_h
#define TPZElasticResponse_h

#include "TPZPlasticState.h"
#include "TPZPlasticStep.h"
#include "TPZElasticResponse.h"

class TPZElasticCriterion : public TPZPlasticBase, public TPZPlasticCriterion
{
    
public:
    
    enum {NYield=3};
    
    typedef TPZElasticCriterion fNYields;
    
    TPZPlasticState<STATE> fN;
    
    TPZElasticResponse fER;
    
public:
    
    virtual int ClassId() const override;
    
    TPZElasticCriterion();
    
    TPZElasticCriterion(const TPZElasticCriterion &cp);
    
    TPZElasticCriterion & operator=(const TPZElasticCriterion &cp);
    
    
    void Read(TPZStream& buf, void* context) override;
    
    
    void Write(TPZStream& buf, int withclassid) const override;
    
    virtual int IntegrationSteps()const override;
    
    
    virtual const char * Name() const override
    {
        return "TPZElasticCriterion";
    }
    
    
    virtual void Print(std::ostream & out) const override
    {
        out << "Classe: " << this->Name();
        fN.Print(out);
    }
    
    virtual void ApplyStrain(const TPZTensor<REAL> &epsTotal) override;
    
    virtual void ApplyStrainComputeSigma(const TPZTensor<REAL> &epsTotal, TPZTensor<REAL> &sigma, TPZFMatrix<REAL> * De = NULL) override;
    
    
    virtual void ApplyStrainComputeDep(const TPZTensor<REAL> &epsTotal, TPZTensor<REAL> &sigma, TPZFMatrix<REAL> &De) override {
        std::cerr << "Deprecated gradient calculation is incorporated on ApplyStrainComputeSigma method." << std::endl;
        DebugStop();
    }
    
    virtual void ApplyLoad(const TPZTensor<REAL> & sigma, TPZTensor<REAL> &epsTotal) override;
    

    virtual void SetState(const TPZPlasticState<REAL> &state) override;
    
    virtual TPZPlasticState<REAL> GetState() const override;
    
    
    virtual void Phi(const TPZTensor<REAL> &epsTotal, TPZVec<REAL> &phi) const override;
    
    
    virtual void SetElasticResponse(TPZElasticResponse &ER) override
    {
        fER = ER;
    }
    
    
    virtual TPZElasticResponse GetElasticResponse() const override
    {
        return fER;
    }
    
    
    TPZPlasticCriterion & GetYC() override{
        return *this;
    }
    
    
    void YieldFunction(const TPZVec<STATE>& sigma, STATE hardening, TPZVec<STATE>& yield) const override{
        TPZTensor<STATE> sigmaTensor;
        sigmaTensor.XX() = sigma[0];
        sigmaTensor.YY() = sigma[1];
        sigmaTensor.ZZ() = sigma[2];
        Phi(sigmaTensor, yield);
    }
    
    virtual int GetNYield() const override {
        return as_integer(NYield);
    }
    
    
};


#endif /* TPZElasticCriterion_h */