pzbctension.h

Go to the documentation of this file.

#ifndef BCTENSIONHPP
#define BCTENSIONHPP

#include "pzbndcond.h"
#include "TPZPlacaOrthotropic.h"
#include "TPZMulticamadaOrtho.h"


template <class T, int N>
class TPZManVector;
class TPZInterpolatedElement;
class TPZMulticamadaOrthotropic;


class TPZBCTension : public TPZBndCond {
                
                TPZMulticamadaOrthotropic *fMultCam;
                int fCamada;
                REAL fSign;
                
private:
                
                public :
    
    ~TPZBCTension(){}
                
                TPZBCTension(TPZMaterial * &material,int id,int type,TPZFMatrix<STATE> &val1,TPZFMatrix<STATE> &val2, REAL sign, TPZMulticamadaOrthotropic *mult, int camada);
                
                virtual int NFluxes() override { return Material()->NFluxes(); }
                
                virtual int NStateVariables() const override { return Material()->NStateVariables(); }
                
                void Contribute(TPZMaterialData &data, REAL weight,TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef)  override {

                                TPZFMatrix<REAL> &dphi = data.dphix;
                                TPZFMatrix<REAL> &phi = data.phi;
                                TPZManVector<REAL,3> &x = data.x;
                                TPZVec<STATE> &sol=data.sol[0];
                                TPZFMatrix<STATE> &dsol=data.dsol[0];
                                TPZFMatrix<REAL> &jacinv = data.jacinv;
                                TPZFMatrix<REAL> &axes = data.axes;

                                int typekeep = fType;
                                if(fType == 4) {
                                                TPZManVector<REAL,3> normal(3);
                                                normal[0] = fSign*(axes(0,1)*axes(1,2)-axes(0,2)*axes(1,1));
                                                normal[1] = fSign*(axes(0,2)*axes(1,0)-axes(0,0)*axes(1,2));
                                                normal[2] = fSign*(axes(0,0)*axes(1,1)-axes(0,1)*axes(1,0));
                                                TPZFNMatrix<9> tensor(3,3);
                                                fMultCam->Tensor(x,fCamada,tensor);
                                                int i,j;
                                                Val2().Zero();
                                                for(i=0; i<3; i++) {
                                                                for(j=0; j<3; j++) {
                                                                                Val2()(i,0) += tensor(i,j)*normal[j];
                                                                }
                                                }
                                                fType = 1;
                                                Material()->ContributeBC(data,weight,ek,ef,*this);
                                                fType = typekeep;
                                } else {
                                                TPZMaterialData data;
                                                data.x = x;
                                                data.jacinv = jacinv;
                                                data.sol[0] = sol;
                                                data.dsol[0] = dsol;
                                                data.axes = axes;
                                                data.phi = phi;
                                                data.dphix = dphi;
                                                TPZBndCond::Contribute(data,weight,ek,ef);
                                }
                                
                }
                void Contribute(TPZMaterialData &data, REAL weight,TPZFMatrix<STATE> &ef) override
                {
                                TPZBndCond::Contribute(data,weight,ef);
                }
public:
virtual int ClassId() const override;

};

#endif