pzthermicelast3d.h

Go to the documentation of this file.

#ifndef PZTHERMICELAST3D
#define PZTHERMICELAST3D

#include "pzelast3d.h"

class TPZThermicElast3D : public TPZElasticity3D {
                
private:
                
                bool fReferred;
                
                STATE fFinalTemperature;
                
public:
                
                STATE fThermalCoeff;
                
                STATE fRefTemperature;
                
                TPZThermicElast3D(int nummat, STATE ThermalCoeff, STATE RefTemp, STATE E, STATE poisson, TPZVec<STATE> &force);
                
                virtual ~TPZThermicElast3D();
                
                void SetConstantTemperatureField(STATE FinalTemp){
                                this->fFinalTemperature = FinalTemp;
                                this->fReferred = false;
                }
                
                void SetReferredTemperatureField(){
                                this->fFinalTemperature = 0.; //any value, it will not be used
                                this->fReferred = true;
                }
                
                bool IsReferredTemperatureField(){
                                return this->fReferred;
                }
                
                void ContributeThermalStress(TPZVec<STATE> &sol, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi, REAL weight, TPZFMatrix<STATE> &ef);
                
                virtual void Contribute(TPZMaterialData &data,
                                                                                                                REAL weight,
                                                                                                                TPZFMatrix<STATE> &ek,
                                                                                                                TPZFMatrix<STATE> &ef) override;
                
                virtual void Contribute(TPZMaterialData &data,
                                                                                                                REAL weight,
                                                                                                                TPZFMatrix<STATE> &ef) override
                {
                                TPZElasticity3D::Contribute(data,weight,ef);
                }
                
protected:
                virtual void Solution(TPZVec<STATE> &Sol, TPZFMatrix<STATE> &DSol,
                                                                                                  TPZFMatrix<REAL> &axes, int var, TPZVec<STATE> &Solout) override;
public:
                virtual void Solution(TPZMaterialData &data, int var, TPZVec<STATE> &Solout) override
                {
        int numbersol = data.sol.size();
        if (numbersol != 1) {
            DebugStop();
        }
                                Solution(data.sol[0],data.dsol[0],data.axes,var,Solout);
                }
                
                
    public:
int ClassId() const override;
 
};

#endif