BCs_RVELagrange_Trac::CommonFunctions Struct Reference

#include "users_modules/mofem_um_homogenisation/src/BCs_RVELagrange_Trac.hpp"

Collaboration diagram for BCs_RVELagrange_Trac::CommonFunctions:

Public Member Functions
PetscErrorCode	shapeMat (int rank, unsigned int gg, DataForcesAndSurcesCore::EntData &col_data, MatrixDouble &n)
PetscErrorCode	hMat (VectorDouble face_normal, int rank, MatrixDouble &n_mat, MatrixDouble &h_mat)

Public Attributes
MatrixDouble	H_mat_1Node

Detailed Description

Definition at line 31 of file BCs_RVELagrange_Trac.hpp.

Member Function Documentation

hMat()

PetscErrorCode BCs_RVELagrange_Trac::CommonFunctions::hMat ( VectorDouble  face_normal, int  rank, MatrixDouble &  n_mat, MatrixDouble &  h_mat  )

inline

Definition at line 58 of file BCs_RVELagrange_Trac.hpp.

      {
      PetscFunctionBegin;
 
      switch(rank) {
        //Mechanical problem
        case 3: {
          H_mat_1Node.resize(6,3);
          H_mat_1Node.clear(); //for one node
          if(face_normal[0]>0) {     //+X face of the RVE
            H_mat_1Node(0,0) = 1.0;
            H_mat_1Node(3,1) = 1.0;
            H_mat_1Node(5,2) = 1.0;
          } else if(face_normal[0]<0) {    //-X face of the RVE
            H_mat_1Node(0,0) = -1.0;
            H_mat_1Node(3,1) = -1.0;
            H_mat_1Node(5,2) = -1.0;
          } else if(face_normal[1]>0) {     //+Y face of the RVE
            H_mat_1Node(1,1) = 1.0;
            H_mat_1Node(3,0) = 1.0;
            H_mat_1Node(4,2) = 1.0;
          } else if(face_normal[1]<0) {    //-Y face of the RVE
            H_mat_1Node(1,1) = -1.0;
            H_mat_1Node(3,0) = -1.0;
            H_mat_1Node(4,2) = -1.0;
          } else if(face_normal[2]>0) {    //+Z face of the RVE
            H_mat_1Node(2,2) = 1.0;
            H_mat_1Node(4,1) = 1.0;
            H_mat_1Node(5,0) = 1.0;
          } else if(face_normal[2]<0) {    //-Z face of the RVE
            H_mat_1Node(2,2) = -1.0;
            H_mat_1Node(4,1) = -1.0;
            H_mat_1Node(5,0) = -1.0;
          } else {
            SETERRQ(PETSC_COMM_SELF,MOFEM_IMPOSIBLE_CASE,"Can not be ?!");
          }
 
          int num_col = n_mat.size2();
          h_mat.resize(6,num_col);
 
          int cc1 = 0;
          for(int bb = 0; bb<num_col/3; bb++) {
            //blocks of 6x3
            for(int rr = 0; rr<6; rr++) {
              for(int cc = 0; cc<3; cc++) {
                h_mat(rr,(cc+cc1)) = H_mat_1Node(rr,cc);
              }
            }
            cc1+=3;
          }
 
        }
        break;
        case 1: //Moisture transport or thermal problem
        {
          H_mat_1Node.resize(3,1);
          H_mat_1Node.clear();
          if(face_normal[0]>0) {     //+X face of the RVE
            H_mat_1Node(0,0) = 1.0;
          }
          if(face_normal[0]<0) {    //-X face of the RVE
            H_mat_1Node(0,0) = -1.0;
          }
          if(face_normal[1]>0) {     //+Y face of the RVE
            H_mat_1Node(1,0) = 1.0;
          }
          if(face_normal[1]<0) {    //-Y face of the RVE
            H_mat_1Node(1,0) = -1.0;
          }
          if(face_normal[2]>0) {    //+Z face of the RVE
            H_mat_1Node(2,0) = 1.0;
          }
          if(face_normal[2]<0) {    //-Z face of the RVE
            H_mat_1Node(2,0) = -1.0;
          }
 
          int num_col = n_mat.size2();
          h_mat.resize(3,num_col);
 
          int cc1 = 0;
          for(int bb = 0; bb<num_col; bb++) {
            //blocks of 3x1
            for(int rr = 0; rr<3; rr++) {
              for(int cc = 0; cc<1; cc++) {
                h_mat(rr,(cc+cc1))=H_mat_1Node(rr,cc);
              }
            }
            cc1+=1;
          }
 
        }
        break;
        default:
        SETERRQ(PETSC_COMM_SELF,MOFEM_NOT_IMPLEMENTED,"not implemented");
      }
      PetscFunctionReturn(0);
    }

shapeMat()

PetscErrorCode BCs_RVELagrange_Trac::CommonFunctions::shapeMat ( int  rank, unsigned int  gg, DataForcesAndSurcesCore::EntData &  col_data, MatrixDouble &  n  )

inline

Definition at line 33 of file BCs_RVELagrange_Trac.hpp.

      {
      PetscFunctionBegin;
 
      int shape_size = col_data.getN().size2();
 
      n.resize(rank,shape_size*rank);
      n.clear();
 
      int kk = 0;
      for(int ii=0; ii<shape_size; ii++){
        double val = col_data.getN()(gg,ii);
        for(int jj=0; jj<rank; jj++){
          n(jj,kk) = val;
          kk++;
        }
      }
 
      PetscFunctionReturn(0);
    }

Member Data Documentation

H_mat_1Node

MatrixDouble BCs_RVELagrange_Trac::CommonFunctions::H_mat_1Node

Definition at line 56 of file BCs_RVELagrange_Trac.hpp.

---

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

• users_modules/mofem_um_homogenisation/src/BCs_RVELagrange_Trac.hpp