\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`
	)

inline

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
	)

inline

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);
     }