TPZWillamWarnke.h

Go to the documentation of this file.

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


/* Generated by Together */// $Id: TPZDruckerPrager.h,v 1.24 2009-09-24 04:19:52 Diogo Exp $

#ifndef TZPWILLAMWARNKE_H
#define TZPWILLAMWARNKE_H

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


#define WILLAMWARNKEPARENT TPZPlasticStep<TPZYCWillamWarnke, TPZThermoForceA, TPZElasticResponse>


class TPZWillamWarnke : public WILLAMWARNKEPARENT  {
                
public:
                
                enum {NYield = TPZYCWillamWarnke::NYield};
                
public:
                
    TPZWillamWarnke():WILLAMWARNKEPARENT(), faPa(0.), fInitialEps()
    {
                                fMaterialTensionSign  = 1; // internally in this material tension is negative
                                fInterfaceTensionSign =  1; // by default
    }
                
    TPZWillamWarnke(const TPZWillamWarnke & source):WILLAMWARNKEPARENT(source)
    {
                                faPa                            = source.faPa;
                                fInitialEps = source.fInitialEps;
    }
                
    TPZWillamWarnke & operator=(const TPZWillamWarnke & source)
    {
                                WILLAMWARNKEPARENT::operator=(source);
                                faPa                            = source.faPa;
                                fInitialEps = source.fInitialEps;
                                
                                return *this;
    }
                
                virtual const char * Name() const override
                {
                                return "TPZWillamWarnke";       
                }
                
                //              void SetUp()
                //              {
                //                              
                //                              REAL pi = 3.141592653589793;
                //                              DRUCKERPARENT::fYC.SetUp(/*phi=20*/ 20./180. * pi ,/*innerMCFit*/0);
                //                              DRUCKERPARENT::fTFA.SetUp(/*yield- coesao inicial correspondeno a fck igual 32 Mpa */ 9.2376, /*k Modulo de hardening da coesao equivante 1 Mpa a cada 0.1% de deformacao */1000.);
                //                              DRUCKERPARENT::fER.SetUp(/*young*/ 20000., /*poisson*/ 0.2);
                //                              TPZTensor<REAL> nullSigma,epsA;
                //                              fInitialEps = DRUCKERPARENT::GetState();
                //                              
                //              }
                
                virtual void Print(std::ostream & out) const override
                {
                                out << "\n" << this->Name();
                                out << "\n Base Class Data:\n";
                                WILLAMWARNKEPARENT::Print(out);
                                out << "\nTPZDruckerPrager internal members:";
                                out << "\n a*Pa = " << faPa;
                                out << "\n InitialEps = " << fInitialEps;
                                
                }
                
    public:
    int ClassId() const override;

    void Write(TPZStream &buf, int withclassid) const override{
        buf. Write(&faPa, 1);
        fInitialEps.Write(buf, withclassid);
    }

    void Read(TPZStream& buf, void* context) override {
        buf. Read(&faPa, 1);
        fInitialEps.Read(buf, context);
    }
    
    virtual int GetNYield() const {
        return as_integer(NYield);
    }
                
private:
                
                

    REAL faPa;
                

    TPZPlasticState<STATE> fInitialEps;
                
                
};


#endif //TPZDruckerPrager_H