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

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

Inheritance diagram for BCs_RVELagrange_Periodic::OpDmatRhs:

Collaboration diagram for BCs_RVELagrange_Periodic::OpDmatRhs:

Public Member Functions
	OpDmatRhs (const string field_name, const string lagrang_field_name, RVEBC_Data_Periodic &data, CommonDataPeriodic &common_data_periodic, CommonFunctionsPeriodic &common_functions_periodic, bool ho_geometry=false)

PetscErrorCode	calculateDmat (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
RVEBC_Data_Periodic &	dAta

bool	hoGeometry

CommonDataPeriodic &	commonDataPeriodic

CommonFunctionsPeriodic &	commonFunctionsPeriodic

ublas::vector< MatrixDouble >	X_mat

MatrixDouble	N_mat

Detailed Description

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

Definition at line 493 of file BCs_RVELagrange_Periodic.hpp.

Constructor & Destructor Documentation

◆ OpDmatRhs()

BCs_RVELagrange_Periodic::OpDmatRhs::OpDmatRhs	(	const string	field_name,
		const string	lagrang_field_name,
		RVEBC_Data_Periodic &	data,
		CommonDataPeriodic &	common_data_periodic,
		CommonFunctionsPeriodic &	common_functions_periodic,
		bool	ho_geometry = `false`
	)

Definition at line 500 of file BCs_RVELagrange_Periodic.hpp.

     :
    FlatPrismElementForcesAndSurcesCore::UserDataOperator(
      lagrang_field_name,UserDataOperator::OPROW
    ),
    dAta(data),
    hoGeometry(ho_geometry),
    commonDataPeriodic(common_data_periodic),
    commonFunctionsPeriodic(common_functions_periodic) {
    }

Member Function Documentation

◆ calculateDmat()

PetscErrorCode BCs_RVELagrange_Periodic::OpDmatRhs::calculateDmat	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

Definition at line 520 of file BCs_RVELagrange_Periodic.hpp.

                                                                                                 {
      PetscFunctionBegin;
//      if(data.getIndices().size()==0) PetscFunctionReturn(0);
//      if(type == MBEDGE && side >= 6) PetscFunctionReturn(0); //ignore second triangel edges (cananical number 6 to 8)
//      if(type == MBTRI && side == 4) PetscFunctionReturn(0);  //ignore second triangel face (cananical number 4)
      
      int rank = data.getFieldDofs()[0]->getNbOfCoeffs();
 
      double x,y,z,x1,y1,z1;
      for(unsigned int gg = 0;gg<data.getN().size1();gg++) {
        double area;
        VectorDouble normal_f3 = getNormalsAtGaussPtsF3(gg);
        area = cblas_dnrm2(3,&normal_f3[0],1)*0.5;
        double val = getGaussPts()(2,gg)*area;
 
        // Gauss point coordinates of the prism for both faces with HO geometry
        x = getHOCoordsAtGaussPtsF3()(gg,0);
        y = getHOCoordsAtGaussPtsF3()(gg,1);
        z = getHOCoordsAtGaussPtsF3()(gg,2);
 
        x1 = getHOCoordsAtGaussPtsF4()(gg,0);
        y1 = getHOCoordsAtGaussPtsF4()(gg,1);
        z1 = getHOCoordsAtGaussPtsF4()(gg,2);
 
//        cout<<"x "<<x << " y "<< y << " z   "<<  z <<endl;
//        cout<<"x1 "<<x1 << " y1 "<< y1 << " z1   "<<  z1 <<endl;
 
        commonDataPeriodic.D_mat.resize(2);
        X_mat.resize(2);
 
        switch(rank) {
          case 3: //mech problem
          X_mat[0].resize(rank,6,false);  X_mat[0].clear();
          X_mat[0](0,0)=2.0*x;  X_mat[0](0,3)=y;  X_mat[0](0,5)=z;
          X_mat[0](1,1)=2.0*y;  X_mat[0](1,3)=x;  X_mat[0](1,4)=z;
          X_mat[0](2,2)=2.0*z;  X_mat[0](2,4)=y;  X_mat[0](2,5)=x;
          X_mat[0]=0.5*X_mat[0];
 
          X_mat[1].resize(rank,6,false); X_mat[1].clear();
          X_mat[1](0,0)=2.0*x1;  X_mat[1](0,3)=y1;  X_mat[1](0,5)=z1;
          X_mat[1](1,1)=2.0*y1;  X_mat[1](1,3)=x1;  X_mat[1](1,4)=z1;
          X_mat[1](2,2)=2.0*z1;  X_mat[1](2,4)=y1;  X_mat[1](2,5)=x1;
          X_mat[1]=0.5*X_mat[1];
          break;
          case 1:  //moisture transport or thermal problem
          X_mat[0].resize(rank,3,false); X_mat[0].clear();
          X_mat[0](0,0)=x; X_mat[0](0,1)=y; X_mat[0](0,2)=z;
 
          X_mat[1].resize(rank,3,false);  X_mat[1].clear();
          X_mat[1](0,0)=x1;  X_mat[1](0,1)=y1;  X_mat[1](0,2)=z1;
          break;
          default:
          SETERRQ(PETSC_COMM_SELF,1,"not implemented");
        }
 
        //
        if(type == MBVERTEX) {
          ierr = commonFunctionsPeriodic.shapeMat(rank, gg, data, N_mat, 2); CHKERRQ(ierr);
        } else {
          ierr = commonFunctionsPeriodic.shapeMat(rank, gg, data, N_mat); CHKERRQ(ierr);
        }
 
        if(gg==0) {
          commonDataPeriodic.D_mat[0]=-1*val*prod(trans(N_mat),X_mat[0]);
          commonDataPeriodic.D_mat[1]=   val*prod(trans(N_mat),X_mat[1]);
        } else {
          commonDataPeriodic.D_mat[0]+=-1*val*prod(trans(N_mat),X_mat[0]);
          commonDataPeriodic.D_mat[1]+=   val*prod(trans(N_mat),X_mat[1]);
        }
      }
 
//      cout<<endl<<endl;
 
      PetscFunctionReturn(0);
    }