BCs_RVELagrange_Disp::OpDmatRhs Struct Reference

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

Inheritance diagram for BCs_RVELagrange_Disp::OpDmatRhs:

Collaboration diagram for BCs_RVELagrange_Disp::OpDmatRhs:

Public Member Functions
	OpDmatRhs (const string field_name, const string lagrang_field_name, RVEBC_Data &data, bool ho_geometry=false)
PetscErrorCode	calculateDmat (int side, EntityType type, DataForcesAndSurcesCore::EntData &data)

Public Attributes
RVEBC_Data &	dAta
bool	hoGeometry
VectorDouble	f
VectorDouble	applied_strain
MatrixDouble	X_mat
MatrixDouble	N_mat
MatrixDouble	D_mat
int	rank
int	nb_row_dofs

Detailed Description

Definition at line 212 of file BCs_RVELagrange_Disp.hpp.

Constructor & Destructor Documentation

OpDmatRhs()

BCs_RVELagrange_Disp::OpDmatRhs::OpDmatRhs ( const string  field_name, const string  lagrang_field_name, RVEBC_Data &  data, bool  ho_geometry = false  )

inline

Definition at line 217 of file BCs_RVELagrange_Disp.hpp.

     :
    FaceElementForcesAndSourcesCore::UserDataOperator(
      lagrang_field_name,UserDataOperator::OPROW
    ),
    dAta(data),hoGeometry(ho_geometry) {
    }

Member Function Documentation

calculateDmat()

PetscErrorCode BCs_RVELagrange_Disp::OpDmatRhs::calculateDmat ( int  side, EntityType  type, DataForcesAndSurcesCore::EntData &  data  )

inline

Definition at line 232 of file BCs_RVELagrange_Disp.hpp.

                                                                                                {
      PetscFunctionBegin;
 
      double x,y,z;
      for(unsigned int gg = 0;gg<data.getN().size1();gg++) {
 
        double area;
        if(hoGeometry) {
          area = norm_2(getNormalsAtGaussPt(gg))*0.5;
        } else {
          area = getArea();
        }
        double val = getGaussPts()(2,gg)*area;
 
        x = getHoCoordsAtGaussPts()(gg,0);
        y = getHoCoordsAtGaussPts()(gg,1);
        z = getHoCoordsAtGaussPts()(gg,2);
 
        switch(rank) {
          case 3: //mech problem
          X_mat.resize(3,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(3,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,1,"not implemented");
        }
 
        int shape_size = data.getN().size2();
        N_mat.resize(rank,shape_size*rank);
        N_mat.clear();
 
        // FIXME: this is inefficient implementation
        {
          int kk=0;
          for(int ii=0; ii<shape_size; ii++){
            //number of shape functions
            double val = data.getN()(gg,ii);
            for(int jj=0; jj<rank; jj++) {
              N_mat(jj,kk) = val;
              kk++;
            }
          }
          if(gg==0) {
            D_mat = val*prod(trans(N_mat),X_mat);
          } else{
            D_mat += val*prod(trans(N_mat),X_mat);
          }
        }
 
      }
 
      PetscFunctionReturn(0);
    }

Member Data Documentation

applied_strain

VectorDouble BCs_RVELagrange_Disp::OpDmatRhs::applied_strain

Definition at line 227 of file BCs_RVELagrange_Disp.hpp.

D_mat

MatrixDouble BCs_RVELagrange_Disp::OpDmatRhs::D_mat

Definition at line 228 of file BCs_RVELagrange_Disp.hpp.

dAta

RVEBC_Data& BCs_RVELagrange_Disp::OpDmatRhs::dAta

Definition at line 214 of file BCs_RVELagrange_Disp.hpp.

f

VectorDouble BCs_RVELagrange_Disp::OpDmatRhs::f

Definition at line 227 of file BCs_RVELagrange_Disp.hpp.

hoGeometry

bool BCs_RVELagrange_Disp::OpDmatRhs::hoGeometry

Definition at line 215 of file BCs_RVELagrange_Disp.hpp.

N_mat

MatrixDouble BCs_RVELagrange_Disp::OpDmatRhs::N_mat

Definition at line 228 of file BCs_RVELagrange_Disp.hpp.

nb_row_dofs

int BCs_RVELagrange_Disp::OpDmatRhs::nb_row_dofs

Definition at line 230 of file BCs_RVELagrange_Disp.hpp.

rank

int BCs_RVELagrange_Disp::OpDmatRhs::rank

Definition at line 230 of file BCs_RVELagrange_Disp.hpp.

X_mat

MatrixDouble BCs_RVELagrange_Disp::OpDmatRhs::X_mat

Definition at line 228 of file BCs_RVELagrange_Disp.hpp.

---

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

• users_modules/mofem_um_homogenisation/src/BCs_RVELagrange_Disp.hpp