\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

Public Member Functions

Public Attributes