BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal Struct Reference

\biref operator to calculate the RHS of the constrain for the RVE boundary conditions More...

#include <users_modules/homogenisation/src/BCs_RVELagrange_Trac.hpp>

Inheritance diagram for BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal:

Collaboration diagram for BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal:

Public Member Functions
	OpRVEBCsRhs_Cal (const string field_name, RVEBC_Data &data, CommonData &common_data, CommonFunctions &common_functions, bool ho_geometry=false)
PetscErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
RVEBC_Data &	dAta
bool	hoGeometry
CommonData &	commonData
CommonFunctions &	commonFunctions
MatrixDouble	X_mat
MatrixDouble	N_mat
MatrixDouble	NTX
MatrixDouble	H_mat

Detailed Description

\biref operator to calculate the RHS of the constrain for the RVE boundary conditions

Definition at line 263 of file BCs_RVELagrange_Trac.hpp.

Constructor & Destructor Documentation

OpRVEBCsRhs_Cal()

BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::OpRVEBCsRhs_Cal ( const string  field_name, RVEBC_Data &  data, CommonData &  common_data, CommonFunctions &  common_functions, bool  ho_geometry = false  )

inline

Definition at line 270 of file BCs_RVELagrange_Trac.hpp.

     :
    FaceElementForcesAndSourcesCore::UserDataOperator(field_name, UserDataOperator::OPCOL),
    dAta(data),
    hoGeometry(ho_geometry),
    commonData(common_data),
    commonFunctions(common_functions) {
    }

Member Function Documentation

doWork()

PetscErrorCode BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::doWork ( int  side, EntityType  type, DataForcesAndSourcesCore::EntData &  data  )

inline

Definition at line 282 of file BCs_RVELagrange_Trac.hpp.

                                                                                          {
      PetscFunctionBegin;
 
      if(type!=MBVERTEX) PetscFunctionReturn(0);
      if(data.getIndices().size()==0) PetscFunctionReturn(0);
 
      auto weak_ptr_dof =
          getNumeredEntFiniteElementPtr()->getRowDofsByPetscGlobalDofIdx(
              data.getIndices()[0]);
      const FENumeredDofEntity *dof_ptr;
      if (auto ptr = weak_ptr_dof.lock())
        dof_ptr = ptr.get();
      else
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Dof not found");
 
      int rank = dof_ptr->getNbOfCoeffs();
 
      double x,y,z;
      for(unsigned int gg = 0;gg<data.getN().size1();gg++) {
        double area;
        if(hoGeometry) {
          area = norm_2(getNormalsAtGaussPts(gg))*0.5;
        } else {
          area = getArea();
        }
        double val = getGaussPts()(2,gg)*area;
 
        x = getCoordsAtGaussPts()(gg,0);
        y = getCoordsAtGaussPts()(gg,1);
        z = getCoordsAtGaussPts()(gg,2);
 
        switch(rank) {
          case 3: //mech problem
          X_mat.resize(rank,6,false);
          X_mat.clear();
          X_mat(0,0) = 2.0*x;  X_mat(0,3) = y;  X_mat(0,5) = z;
          X_mat(1,1) = 2.0*y;  X_mat(1,3) = x;  X_mat(1,4) = z;
          X_mat(2,2) = 2.0*z;  X_mat(2,4) = y;  X_mat(2,5) = x;
          X_mat=0.5*X_mat;
          break;
          case 1:  //moisture transport or thermal problem
          X_mat.resize(rank,3,false);
          X_mat.clear();
          X_mat(0,0) = x;
          X_mat(0,1) = y;
          X_mat(0,2) = z;
          break;
          default:
          SETERRQ(PETSC_COMM_SELF,MOFEM_NOT_IMPLEMENTED,"not implemented");
        }
        ierr = commonFunctions.shapeMat(rank, gg,  data, N_mat); CHKERRQ(ierr);
        ierr = commonFunctions.hMat(getNormalsAtGaussPts(gg), rank, N_mat, H_mat); CHKERRQ(ierr);
 
        if(gg==0) {
          NTX = val*prod(trans(N_mat),X_mat);
          commonData.D_mat = prod(H_mat, NTX);
        } else {
          noalias(NTX) = val*prod(trans(N_mat),X_mat);
          commonData.D_mat += prod(H_mat, NTX);
        }
      }
 
      PetscFunctionReturn(0);
    }

Member Data Documentation

commonData

CommonData& BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::commonData

Definition at line 267 of file BCs_RVELagrange_Trac.hpp.

commonFunctions

CommonFunctions& BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::commonFunctions

Definition at line 268 of file BCs_RVELagrange_Trac.hpp.

dAta

RVEBC_Data& BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::dAta

Definition at line 265 of file BCs_RVELagrange_Trac.hpp.

H_mat

MatrixDouble BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::H_mat

Definition at line 280 of file BCs_RVELagrange_Trac.hpp.

hoGeometry

bool BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::hoGeometry

Definition at line 266 of file BCs_RVELagrange_Trac.hpp.

N_mat

MatrixDouble BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::N_mat

Definition at line 280 of file BCs_RVELagrange_Trac.hpp.

NTX

MatrixDouble BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::NTX

Definition at line 280 of file BCs_RVELagrange_Trac.hpp.

X_mat

MatrixDouble BCs_RVELagrange_Trac::OpRVEBCsRhs_Cal::X_mat

Definition at line 280 of file BCs_RVELagrange_Trac.hpp.

---

The documentation for this struct was generated from the following file:

• users_modules/homogenisation/src/BCs_RVELagrange_Trac.hpp