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

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

Inheritance diagram for BCs_RVELagrange_Disp::OpRVEHomoStress:

Collaboration diagram for BCs_RVELagrange_Disp::OpRVEHomoStress:

Public Member Functions
	OpRVEHomoStress (const string field_name, const string lagrang_field_name, Vec stress_homo, RVEBC_Data &data, bool ho_geometry=false)

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

Public Member Functions inherited from BCs_RVELagrange_Disp::OpDmatRhs
	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
Vec	Stress_Homo

VectorDouble	Stress_Homo_elem

Public Attributes inherited from BCs_RVELagrange_Disp::OpDmatRhs
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

\biref operator to calculate the RVE homogenised stress

Definition at line 586 of file BCs_RVELagrange_Disp.hpp.

Constructor & Destructor Documentation

◆ OpRVEHomoStress()

BCs_RVELagrange_Disp::OpRVEHomoStress::OpRVEHomoStress	(	const string	field_name,
		const string	lagrang_field_name,
		Vec	stress_homo,
		RVEBC_Data &	data,
		bool	ho_geometry = `false`
	)

inline

Definition at line 590 of file BCs_RVELagrange_Disp.hpp.

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

Member Function Documentation

◆ doWork()

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

inline

Definition at line 603 of file BCs_RVELagrange_Disp.hpp.

                                                                                         {
      PetscFunctionBegin;
      if(data.getIndices().size()==0) PetscFunctionReturn(0);
      if(dAta.tRis.find(getNumeredEntFiniteElementPtr()->getEnt())==dAta.tRis.end()) PetscFunctionReturn(0);
      //        cout<<"OpRVEBCsRhs "<<endl;
      
      
      const FENumeredDofEntity *dof_ptr;
      ierr = getNumeredEntFiniteElementPtr()->getRowDofsByPetscGlobalDofIdx(data.getIndices()[0],&dof_ptr); CHKERRQ(ierr);
      rank = dof_ptr->getNbOfCoeffs();
      nb_row_dofs = data.getIndices().size()/rank;
      
      ierr = calculateDmat(side,type,data); CHKERRQ(ierr);
      
      Stress_Homo_elem.resize(D_mat.size2());
      
      // data.getFieldData() is value of Lagrange multply
      // data.getFieldData() is reaction force (so multiply for -1 to get the force)
      noalias(Stress_Homo_elem) = prod(trans(D_mat), -data.getFieldData());
      
      const int Indices6[6] = {0, 1, 2, 3, 4, 5};
      const int Indices3[3] = {0, 1, 2};
      switch(rank) {
        case 3:
          ierr = VecSetValues(Stress_Homo,6,Indices6,&(Stress_Homo_elem.data())[0],ADD_VALUES); CHKERRQ(ierr);
          break;
        case 1:
          ierr = VecSetValues(Stress_Homo,3,Indices3,&(Stress_Homo_elem.data())[0],ADD_VALUES); CHKERRQ(ierr);
          break;
        default:
          SETERRQ(PETSC_COMM_SELF,1,"not implemented");
      }
      
      PetscFunctionReturn(0);
    }

Member Data Documentation

◆ Stress_Homo

Vec BCs_RVELagrange_Disp::OpRVEHomoStress::Stress_Homo

Definition at line 588 of file BCs_RVELagrange_Disp.hpp.

◆ Stress_Homo_elem

VectorDouble BCs_RVELagrange_Disp::OpRVEHomoStress::Stress_Homo_elem

Definition at line 601 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

Public Member Functions

Public Attributes