BCs_RVELagrange_Disp::OpRVEHomoStress Struct Reference

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

#include <users_modules/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, DataForcesAndSourcesCore::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, DataForcesAndSourcesCore::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 610 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 614 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, DataForcesAndSourcesCore::EntData &  data  )

inline

Definition at line 627 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;
      
 
      auto weak_ptr_dof =
          getNumeredEntFiniteElementPtr()->getRowDofsByPetscGlobalDofIdx(
              data.getIndices()[0]);
      const FENumeredDofEntity *dof_ptr;
      if (auto ptr = weak_ptr_dof.lock())
        dof_ptr = ptr.get();
      else
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Dof not found");
 
      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 612 of file BCs_RVELagrange_Disp.hpp.

Stress_Homo_elem

VectorDouble BCs_RVELagrange_Disp::OpRVEHomoStress::Stress_Homo_elem

Definition at line 625 of file BCs_RVELagrange_Disp.hpp.

---

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

• users_modules/homogenisation/src/BCs_RVELagrange_Disp.hpp