TPZVonMises.h

Go to the documentation of this file.

/*
 *  TPZVonMises.h
 *  ElastoPlasticModels
 *
 *  Created by Diogo Cecilio on 12/17/10.
 *  Copyright 2010 __MyCompanyName__. All rights reserved.
 *
 */

#ifndef TPZVONMISES_H
#define TPZVONMISES_H

#include "pzlog.h"
#include "TPZPlasticStep.h"
#include "TPZYCVonMises.h"
#include "TPZThermoForceA.h"
#include "TPZElasticResponse.h"
#include "pzvec_extras.h"
#include "TPZPlasticStepID.h"

#define VONMISESPARENT TPZPlasticStep<TPZYCVonMises, TPZThermoForceA, TPZElasticResponse>


class TPZVonMises : public VONMISESPARENT  {
                
public:
                
                enum {NYield =1};
                
public:
                
    TPZVonMises():VONMISESPARENT()
    {
                                fMaterialTensionSign  = 1; // internally in this material tension is negative
                                fInterfaceTensionSign =  1; // by default
    }
                
    TPZVonMises(const TPZVonMises & source):VONMISESPARENT(source)
    {
    }
                
    TPZVonMises & operator=(const TPZVonMises & source)
    {
                                VONMISESPARENT::operator=(source);                              
                                return *this;
    }
                
                virtual const char * Name() const override
                {
                                return "TPZVonMises";           
                }
                
    virtual void SetUp(const TPZTensor<REAL> & epsTotal) override {
        VONMISESPARENT::SetUp(epsTotal);
    }

                void SetUp()
                {
                                
                                //VONMISESPARENT::fYC.SetUp();
                                //VONMISESPARENT::fTFA.SetUp(210./* sigmaY = MPa */,    1164./*endurecimento para cada 0.1% de deformacao em MPa(290/0.249)*/);
                                //VONMISESPARENT::fER.SetUp(/*young*/ 210000., /*poisson*/ 0.3);
                                
                }
                
                REAL YieldRadius(TPZPlasticState<REAL> state)
                {
                                                                REAL radius = VONMISESPARENT::fTFA.fSigmaYield0;
                                                                return radius;
                }
                
                
                virtual void Print(std::ostream & out) const override
                {
                                out << "\n" << this->Name();
                                out << "\n Base Class Data:\n";
                                VONMISESPARENT::Print(out);                     
                }
                
                public:
int ClassId() const override;

    void Write(TPZStream &buf, int withclassid) const override{
                }
                
    void Read(TPZStream& buf, void* context) override {
                }               
    
    static void Steel(TPZVonMises & material)
    {
        material.fTFA.SetUp(210./* sigmaY = MPa */,   32222.2222/*endurecimento para cada 0.1% de deformacao em MPa(290/0.249)*/);
                                material.fER.SetEngineeringData(/*young*/ 210000., /*poisson*/ 0.);
    }
                
                
public:
                
    virtual int GetNYield() const {
        return as_integer(NYield);
    }
};


#endif //TPZVonMises_H