TPZElasticResponse.h

Go to the documentation of this file.

//
//  TPZElasticResponse.h
//  pz
//
//  Created by Erick Slis Raggio Santos on 04/07/2009.
//


#ifndef TPZElasticResponse_H
#define TPZElasticResponse_H

#include "TPZTensor.h"
#include "pzreal.h"


class TPZElasticResponse : public TPZSavable {
    
    REAL m_lambda;
    
    REAL m_mu;
    
    TPZTensor<REAL> m_eps_star;
    
public:
    
    int ClassId() const override;
    
    TPZElasticResponse();
    
    TPZElasticResponse(const TPZElasticResponse & other);
    
    TPZElasticResponse & operator=(const TPZElasticResponse & other);
    
    void Write(TPZStream &buf, int withclassid) const override;
    
    void Read(TPZStream &buf, void *context) override;
    
    
    const char * Name() const;
    
    void Print(std::ostream & out) const;
    
    template<class T>
    void ComputeStress(const TPZTensor<T> & epsilon, TPZTensor<T> & sigma) const {

        TPZTensor<T> delta_epsilon;
        delta_epsilon.XX() = -m_eps_star.XX();
        delta_epsilon.XY() = -m_eps_star.XY();
        delta_epsilon.XZ() = -m_eps_star.XZ();
        delta_epsilon.YY() = -m_eps_star.YY();
        delta_epsilon.YZ() = -m_eps_star.YZ();
        delta_epsilon.ZZ() = -m_eps_star.ZZ();
        delta_epsilon += epsilon; // Substract residual strain
        T trace = delta_epsilon.I1();
        sigma.Identity();
        sigma.Multiply(trace, m_lambda);
        sigma.Add(delta_epsilon, 2. * m_mu);
    }
    
    template<class T>
    void ComputeStrain(const TPZTensor<T> & sigma, TPZTensor<T> & epsilon) const {
        const T fac = T((1 / 3.)*(1. / (3. * m_lambda + 2. * m_mu) - 1. / (2. * m_mu)));
        REAL trace = sigma.I1();
        epsilon.Identity();
        epsilon.Multiply(trace, fac);
        epsilon.Add(sigma, 1. / (2. * m_mu));
        epsilon += m_eps_star;// Adding residual strain
    }
    
    
    void De(TPZFMatrix<REAL> & De);
    
    
    void SetEngineeringData(REAL Eyoung, REAL Poisson);
    
    
    void SetLameData(REAL lambda, REAL mu);
    
    REAL Lambda() const;
    
    REAL K() const;
    
    REAL Mu() const;
    
    REAL G() const;
    
    REAL E() const;
    
    REAL Poisson() const;
    
    void SetResidualStrainData(TPZTensor<REAL> & eps_res);
    
    TPZTensor<REAL> & ResidualStrainData();
};

#endif /* TPZElasticResponse_h */