\biref operator to calculate the RVE homogenised stress More...

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

Inheritance diagram for BCs_RVELagrange_Periodic::OpRVEHomoStress:

Collaboration diagram for BCs_RVELagrange_Periodic::OpRVEHomoStress:

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

PetscErrorCode	doWork (int side, EntityType type, DataForcesAndSurcesCore::EntData &data)

Public Member Functions inherited from 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)

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

Public Attributes
Vec	Stress_Homo

ublas::vector< VectorDouble >	Stress_Homo_elem

VectorDouble	field_data

Public Attributes inherited from BCs_RVELagrange_Periodic::OpDmatRhs
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 RVE homogenised stress

Definition at line 1227 of file BCs_RVELagrange_Periodic.hpp.

Constructor & Destructor Documentation

◆ OpRVEHomoStress()

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

inline

Definition at line 1231 of file BCs_RVELagrange_Periodic.hpp.

     :
    OpDmatRhs(
      field_name,lagrang_field_name,data,common_data_periodic,common_functions_periodic,ho_geometry
    ),
    Stress_Homo(stress_homo) {
    }

Member Function Documentation

◆ doWork()

PetscErrorCode BCs_RVELagrange_Periodic::OpRVEHomoStress::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSurcesCore::EntData &	data
	)

inline

Definition at line 1249 of file BCs_RVELagrange_Periodic.hpp.

                                                                                         {
      PetscFunctionBegin;
 
      if(data.getIndices().size()==0) PetscFunctionReturn(0);
      if(type == MBEDGE && side >= 3) 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)
 
      
      ierr = calculateDmat(side,type,data); CHKERRQ(ierr);
 
//      cerr << commonDataPeriodic.D_mat[0] << endl;
//      cerr << data.getFieldData() << endl;
//      cerr << rowFieldName << endl;
//      cerr << data.getFieldDofs()[0]->getName() << endl;
//      cerr << data.getFieldDofs()[0]->getNbOfCoeffs() << endl;
//      cerr << data.getFieldDofs().size() << endl;
//      cerr << data.getIndices() << endl;
//      */
 
 
      Stress_Homo_elem.resize(2);
      if(type == MBVERTEX) {
        int nb=data.getFieldData().size()/2;  //for one face for nodes only  (total indieces are for both faces)
        field_data.resize(nb);
 
        for(int ii=0; ii<nb; ii++){
          field_data[ii]=data.getFieldData()[ii];
        }
        Stress_Homo_elem[0] = prod(trans(commonDataPeriodic.D_mat[0]), -1*field_data);   //Lamda=data.getFieldData() is reaction force (so multiply for -1 to get the force)
        Stress_Homo_elem[1] = prod(trans(commonDataPeriodic.D_mat[1]), -1*field_data);   //Lamda=data.getFieldData() is reaction force (so multiply for -1 to get the force)
      }
      else{
        Stress_Homo_elem[0] = prod(trans(commonDataPeriodic.D_mat[0]), -1*data.getFieldData());   //Lamda=data.getFieldData() is reaction force (so multiply for -1 to get the force)
        Stress_Homo_elem[1] = prod(trans(commonDataPeriodic.D_mat[1]), -1*data.getFieldData());   //Lamda=data.getFieldData() is reaction force (so multiply for -1 to get the force)
      }
 
 
//      cerr << "field_data "<<field_data << endl;
//      cerr << "Stress_Homo_elem[0] "<<Stress_Homo_elem[0] << endl;
//      cerr << "Stress_Homo_elem[1] "<<Stress_Homo_elem[1] << endl;
//      cerr << "commonDataPeriodic.D_mat[0] "<<commonDataPeriodic.D_mat[0] << endl;
 
      int Indices6[6]={0, 1, 2, 3, 4, 5};
      int Indices3[3]={0, 1, 2};
 
      int rank = data.getFieldDofs()[0]->getNbOfCoeffs();
      switch(rank) {
        case 3:
        ierr = VecSetValues(Stress_Homo,6,Indices6,&Stress_Homo_elem[0][0],ADD_VALUES); CHKERRQ(ierr);
        ierr = VecSetValues(Stress_Homo,6,Indices6,&Stress_Homo_elem[1][0],ADD_VALUES); CHKERRQ(ierr);
        break;
        case 1:
        ierr = VecSetValues(Stress_Homo,3,Indices3,&Stress_Homo_elem[0][0],ADD_VALUES); CHKERRQ(ierr);
        ierr = VecSetValues(Stress_Homo,3,Indices3,&Stress_Homo_elem[1][0],ADD_VALUES); CHKERRQ(ierr);
        break;
        default:
        SETERRQ(PETSC_COMM_SELF,1,"not implemented");
      }
 
      PetscFunctionReturn(0);
    }