pzelasAXImat.h

Go to the documentation of this file.

#ifndef ELASAXIMATHPP
#define ELASAXIMATHPP

#include <iostream>

#include "TPZMaterial.h"
#include "pzdiscgal.h"

#include <vector>
#include <math.h>

class TPZElasticityAxiMaterial : public TPZDiscontinuousGalerkin {
                
                public :

                TPZElasticityAxiMaterial();
                TPZElasticityAxiMaterial(int num, REAL E, REAL nu, REAL fx, REAL fy);
                
                TPZElasticityAxiMaterial(int num, REAL E, REAL nu, REAL fx, REAL fy, REAL coefTheta,  REAL coefAlpha);
                
    TPZElasticityAxiMaterial(const TPZElasticityAxiMaterial &copy);
    
    virtual TPZMaterial * NewMaterial()  override { return new TPZElasticityAxiMaterial(*this);}
    
    virtual ~TPZElasticityAxiMaterial();
                
                void SetOrigin(TPZManVector<REAL> &Orig, TPZManVector<REAL> &AxisZ, TPZManVector<REAL> &AxisR);
                
                REAL ComputeR(TPZVec<REAL> &x);
                
                void SetMohrCoulomb(double c, double phi)
                {
                                f_c = c;
                                f_phi = phi;
                }
    
    void SetThermalExpansionCoefficient(REAL alpha)
    {
        fAlpha = alpha;
    }
    
    void SetTemperature(REAL delt)
    {
        fDelTemperature = delt;
    }
                
                virtual int Dimension() const  override { return 2;}
                
                virtual  int NStateVariables() const override;
                                                
                virtual void Print(std::ostream & out = std::cout) override;
                
                std::string Name() override { return "TPZElasticityAxiMaterial"; }
                
                virtual short NumberOfFluxes(){return 3;}
                
                virtual int NFluxes() override { return 3;}
                
                virtual void Contribute(TPZMaterialData &data, REAL weight,TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef) override;
                
                virtual void ContributeBC(TPZMaterialData &data,REAL weight,
                                                                                                                  TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef,TPZBndCond &bc) override;
                
                virtual void ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, 
                                     REAL weight, TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef) override;
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft, 
                                       REAL weight, TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef,TPZBndCond &bc) override {
                                
                                PZError << "\nFATAL ERROR - Method not implemented: " << __PRETTY_FUNCTION__ << "\n";
        DebugStop();
                }
                
                virtual int VariableIndex(const std::string &name) override;
                
                virtual int NSolutionVariables(int var) override;
                
                virtual void Solution(TPZMaterialData &data, int var, TPZVec<STATE> &Solout) override;
                
                virtual void Flux(TPZVec<REAL> &x, TPZVec<STATE> &Sol, TPZFMatrix<STATE> &DSol, TPZFMatrix<REAL> &axes, TPZVec<STATE> &flux) override;
                
                void Errors(TPZVec<REAL> &x,TPZVec<STATE> &u,
                                                                TPZFMatrix<STATE> &dudx, TPZFMatrix<REAL> &axes, TPZVec<STATE> &flux,
                                                                TPZVec<STATE> &u_exact,TPZFMatrix<STATE> &du_exact,TPZVec<REAL> &values) override;//Cedric
                
                REAL E() {return fE;}
                
                REAL Nu() {return fnu;}
                
                void SetPreStress(REAL Sigxx, REAL Sigyy, REAL Sigxy);
    
    void SetTemperatureFunction(void (*func)(const TPZVec<REAL> &rz, REAL &temperature))
    {
        fTemperatureFunction = func;
    }
                
                public:
virtual int ClassId() const override;

                
                virtual void Read(TPZStream &buf, void *context) override;
                
                virtual void Write(TPZStream &buf, int withclassid) const override;
                
                TPZManVector<REAL> GetAxisR();
                TPZManVector<REAL> GetAxisZ();
                TPZManVector<REAL> GetOrigin();
                
                REAL fIntegral;
                
private:
                
                REAL f_phi;
                REAL f_c;
                
                REAL fE;
                
                REAL fnu;
                
    REAL fAlpha;
    
                REAL ff[3];
                
    REAL fDelTemperature;
    
                REAL fEover21PlusNu;
                
                REAL fEover1MinNu2;
                
                TPZManVector<REAL> f_AxisR;
                
                TPZManVector<REAL> f_AxisZ;
                
                TPZManVector<REAL> f_Origin;
                
                REAL fSymmetric;
                
                REAL fPenaltyConstant;
    
    void (*fTemperatureFunction)(const TPZVec<REAL> &rz, REAL &temperature);
                
};

#endif