element_structure.hpp

Go to the documentation of this file.

#pragma once

#include "mpi_structure.hpp"

#include <cmath>
#include <iostream>
#include <cstdlib>

#include "config_structure.hpp"
#include "geometry_structure.hpp"
#include "gauss_structure.hpp"

using namespace std;

class CElement {
protected:
    unsigned short nGaussPoints;        
    unsigned short nGaussPointsP;       
    unsigned short nNodes;              
    static unsigned short nDim;     
    su2double **CurrentCoord,               
    **RefCoord;                         
    su2double **GaussCoord,             
    *GaussWeight;                       
    su2double   **GaussCoordP,              
    *GaussWeightP;                      
    su2double **NodalExtrap;                
    su2double **NodalStress;                
    CGaussVariable **GaussPoint;        
    CGaussVariable **GaussPointP;       
    su2double **Mab;                        
    su2double ***Kab;                       
    su2double **Ks_ab;                      
    su2double ***Kk_ab;                 
    su2double **Kt_a;                       
    su2double **FDL_a;                      
    su2double el_Pressure;                  
public:
    CElement(void);

    CElement(unsigned short val_nDim, CConfig *config);

    virtual ~CElement(void);

    unsigned short GetnNodes(void);

    unsigned short GetnGaussPoints(void);

    unsigned short GetnGaussPointsP(void);

    void SetRef_Coord(su2double val_CoordRef, unsigned short iNode, unsigned short iDim);

    void SetCurr_Coord(su2double val_CoordCurr, unsigned short iNode, unsigned short iDim);

    void SetElement_Pressure(su2double val_ElPressure);

    su2double GetRef_Coord(unsigned short iNode, unsigned short iDim);

    su2double GetCurr_Coord(unsigned short iNode, unsigned short iDim);

    su2double GetWeight(unsigned short iGauss);

    su2double GetWeight_P(unsigned short iGaussP);

    su2double GetJ_X(unsigned short iGauss);

    su2double GetJ_x(unsigned short iGauss);

    su2double GetJ_X_P(unsigned short iGauss);

    su2double GetJ_x_P(unsigned short iGauss);

    su2double GetElement_Pressure(void);

    void Add_Mab(su2double val_Mab, unsigned short nodeA, unsigned short nodeB);

    void Add_Kab(su2double **val_Kab, unsigned short nodeA, unsigned short nodeB);

    void Add_Kab_T(su2double **val_Kab, unsigned short nodeA, unsigned short nodeB);


    void Add_Ks_ab(su2double val_Ks_ab, unsigned short nodeA, unsigned short nodeB);

    void Add_Kt_a(su2double *val_Kt_a, unsigned short nodeA);

    void Add_FDL_a(su2double *val_FDL_a, unsigned short nodeA);

    void Set_Kk_ab(su2double **val_Kk_ab, unsigned short nodeA, unsigned short nodeB);

    void clearElement(void);

    void clearStress(void);

    su2double Get_Mab(unsigned short nodeA, unsigned short nodeB);

    su2double *Get_Kab(unsigned short nodeA, unsigned short nodeB);

    su2double Get_Ks_ab(unsigned short nodeA, unsigned short nodeB);

    su2double *Get_Kk_ab(unsigned short nodeA, unsigned short nodeB);

    su2double *Get_Kt_a(unsigned short nodeA);

    su2double *Get_FDL_a(unsigned short nodeA);

    su2double GetNi(unsigned short iNode, unsigned short iGauss);

    su2double GetGradNi_X(unsigned short iNode, unsigned short iGauss, unsigned short iDim);

    su2double GetGradNi_x(unsigned short iNode, unsigned short iGauss, unsigned short iDim);

    su2double GetGradNi_x_P(unsigned short iNode, unsigned short iGaussP, unsigned short iDim);

    su2double GetNi_Extrap(unsigned short iNode, unsigned short iGauss);

    void Add_NodalStress(su2double val_Stress, unsigned short iNode, unsigned short iVar);

    su2double Get_NodalStress(unsigned short iNode, unsigned short iVar);

    virtual void ComputeGrad_Linear(void);

    virtual void ComputeGrad_NonLinear(void);

    virtual void ComputeGrad_Pressure(void);


};

class CTRIA1 : public CElement {

protected:

public:

    CTRIA1(void);

    CTRIA1(unsigned short val_nDim, CConfig *config);

    ~CTRIA1(void);

    void ComputeGrad_Linear(void);

    void ComputeGrad_NonLinear(void);

};


class CQUAD4 : public CElement {

protected:

public:

    CQUAD4(void);

    CQUAD4(unsigned short val_nDim, CConfig *config);

    virtual ~CQUAD4(void);

    void ComputeGrad_Linear(void);

    void ComputeGrad_NonLinear(void);

    virtual void ComputeGrad_Pressure(void);


};

class CQUAD4P1 : public CQUAD4 {

protected:

public:

    CQUAD4P1(void);

    CQUAD4P1(unsigned short val_nDim, CConfig *config);

    ~CQUAD4P1(void);

    void ComputeGrad_Pressure(void);


};

class CTETRA1 : public CElement {

protected:

public:

    CTETRA1(void);

    CTETRA1(unsigned short val_nDim, CConfig *config);

    ~CTETRA1(void);

    void ComputeGrad_Linear(void);

    void ComputeGrad_NonLinear(void);

};

class CHEXA8 : public CElement {

protected:

public:

    CHEXA8(void);

    CHEXA8(unsigned short val_nDim, CConfig *config);

    virtual ~CHEXA8(void);

    void ComputeGrad_Linear(void);

    void ComputeGrad_NonLinear(void);

    virtual void ComputeGrad_Pressure(void);


};

class CHEXA8P1 : public CHEXA8 {

protected:

public:

    CHEXA8P1(void);

    CHEXA8P1(unsigned short val_nDim, CConfig *config);

    ~CHEXA8P1(void);

    void ComputeGrad_Pressure(void);


};

class CBOUND2D : public CElement {

protected:

public:

  CBOUND2D(void);

  CBOUND2D(unsigned short val_nDim, CConfig *config);

  ~CBOUND2D(void);

  void ComputeGrad_Linear(void);

  void ComputeGrad_NonLinear(void);

};


#include "element_structure.inl"