BCs_RVELagrange_Disp::OpRVEBCsRhs Struct Reference

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

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

Inheritance diagram for BCs_RVELagrange_Disp::OpRVEBCsRhs:

Collaboration diagram for BCs_RVELagrange_Disp::OpRVEBCsRhs:

Public Member Functions
	OpRVEBCsRhs (const string field_name, const string lagrang_field_name, Vec f, VectorDouble given_strain, boost::ptr_vector< MethodForForceScaling > &methods_op, 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	F
VectorDouble	givenStrain
boost::ptr_vector< MethodForForceScaling > &	methodsOp
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 RHS of the constrain for the RVE boundary conditions

Definition at line 305 of file BCs_RVELagrange_Disp.hpp.

Constructor & Destructor Documentation

OpRVEBCsRhs()

BCs_RVELagrange_Disp::OpRVEBCsRhs::OpRVEBCsRhs ( const string  field_name, const string  lagrang_field_name, Vec  f, VectorDouble  given_strain, boost::ptr_vector< MethodForForceScaling > &  methods_op, RVEBC_Data &  data, bool  ho_geometry = false  )

inline

Definition at line 310 of file BCs_RVELagrange_Disp.hpp.

     :
    OpDmatRhs(field_name,lagrang_field_name,data,ho_geometry),F(f),givenStrain(given_strain),methodsOp(methods_op) {
    }

Member Function Documentation

doWork()

PetscErrorCode BCs_RVELagrange_Disp::OpRVEBCsRhs::doWork ( int  side, EntityType  type, DataForcesAndSourcesCore::EntData &  data  )

inline

Definition at line 317 of file BCs_RVELagrange_Disp.hpp.

                                                                                          {
      PetscFunctionBegin;
 
      nb_row_dofs = data.getIndices().size();
 
      if(nb_row_dofs ==0) PetscFunctionReturn(0);
      if(dAta.tRis.find(getNumeredEntFiniteElementPtr()->getEnt())==dAta.tRis.end()) PetscFunctionReturn(0);
 
      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 /= rank;
 
      ierr = calculateDmat(side,type,data); CHKERRQ(ierr);
      
      //This is to scale givenStrain vector to apply the strian in steps
      VectorDouble scaled_given_strain;
      scaled_given_strain.resize(6);
      scaled_given_strain=givenStrain;
//      cout <<"givenStrain Before = "<<scaled_given_strain<<endl;
      ierr = MethodForForceScaling::applyScale(getFEMethod(),methodsOp,scaled_given_strain); CHKERRQ(ierr);
//      cout <<"givenStrain  After = "<<scaled_given_strain<<endl;
      f.resize(D_mat.size1(),false);
      noalias(f) = prod(D_mat, scaled_given_strain);
      f*=-1;
 
//      cout <<"f   = "<<f<<endl;
//      string wait;
//      cin>>wait;
 
      Vec myF = F;
      if(F == PETSC_NULL) {
        myF = getFEMethod()->snes_f;
      }
 
      ierr = VecSetValues(myF,data.getIndices().size(),&data.getIndices()[0],&f[0],ADD_VALUES); CHKERRQ(ierr);
      PetscFunctionReturn(0);
    }

Member Data Documentation

F

Vec BCs_RVELagrange_Disp::OpRVEBCsRhs::F

Definition at line 306 of file BCs_RVELagrange_Disp.hpp.

givenStrain

VectorDouble BCs_RVELagrange_Disp::OpRVEBCsRhs::givenStrain

Definition at line 307 of file BCs_RVELagrange_Disp.hpp.

methodsOp

boost::ptr_vector<MethodForForceScaling>& BCs_RVELagrange_Disp::OpRVEBCsRhs::methodsOp

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