pzviscoelastic.h

Go to the documentation of this file.

#ifndef TPZVISCOELASTIC_H
#define TPZVISCOELASTIC_H

const int _XX_ = 0;
const int _XY_ = 1;
const int _XZ_ = 2;
const int _YY_ = 3;
const int _YZ_ = 4;
const int _ZZ_ = 5;


#include <iostream>
#include "pzfmatrix.h"
#include "pzvec.h"
#include <vector>

#include "pzelast3d.h"
#include "TPZMatWithMem.h"

class TPZViscoelastic : public TPZMatWithMem<TPZFMatrix<STATE>, TPZElasticity3D>
{
                                
public:
                enum SOLUTIONVARS{ENone = -1, EViscoStressX = 30, EViscoStressY = 31, EViscoStressZ = 32};
  
                TPZViscoelastic();
                
                TPZViscoelastic(int id);

                TPZViscoelastic(int id,STATE ElaE, STATE poissonE, STATE lambdaV, STATE muV, STATE alpha, STATE deltaT, TPZVec <STATE> &force);

                void SetMaterialDataHooke(STATE ElaE, STATE poissonE, STATE ElaV, STATE poissonV, STATE alpha, STATE deltaT, TPZVec <STATE> &force);
                
                
                virtual void Contribute(TPZMaterialData &data,
                                                                                                                REAL weight,
                                                                                                                TPZFMatrix<STATE> &ek,
                                                                                                                TPZFMatrix<STATE> &ef) override;
                
                void UpdateQsi(TPZMaterialData &data);
                
                virtual int VariableIndex(const std::string &name) override;
                
                virtual int NSolutionVariables(int var) override;
                
                /*
                 Computes the stress. 
                 Remember you cant update qsi if you want to calculate stress 
                 */
                virtual void ComputeStressTensor(TPZFMatrix<STATE> &Stress, TPZMaterialData &data) const override;
                
                virtual void Solution(TPZMaterialData &data, int var, TPZVec<STATE> &Solout) override;

                virtual void FillDataRequirements(TPZMaterialData &data) override;
                
                void Write(TPZStream &buf, int withclassid) const override;
                
                void Read(TPZStream &buf, void *context) override;

                public:
int ClassId() const override;

                
protected:
                
                STATE fAlpha;
                
                STATE fDeltaT;

                STATE fLambdaV;
                STATE fmuV;
};

#endif