Operator used as base struct for OpContactMaterialVolOnSideLhs_dX_dx OpContactMaterialVolOnSideLhs_dX_dX operators that use side volume element adjacent to current contact prism needed to evaluate of deformation gradient tensor derivative. More...

#include <users_modules/mortar_contact/src/SimpleContact.hpp>

Inheritance diagram for SimpleContactProblem::OpContactMaterialVolOnSideLhs:

[legend]

Collaboration diagram for SimpleContactProblem::OpContactMaterialVolOnSideLhs:

[legend]

Public Member Functions
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

virtual MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

	OpContactMaterialVolOnSideLhs (const string field_name_1, const string field_name_2, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, const bool is_master)
	LHS-operator for the contact element (material configuration) More...

Public Attributes
MatrixDouble	matLhs

boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact

VectorInt	rowIndices

VectorInt	colIndices

int	row_nb_dofs

int	col_nb_dofs

int	nb_gauss_pts

int	nb_base_fun_row

int	nb_base_fun_col

boost::shared_ptr< VectorDouble >	tangentOne

boost::shared_ptr< VectorDouble >	tangentTwo

boost::shared_ptr< VectorDouble >	normalVector

double	aRea

Private Attributes
bool	isMaster

Detailed Description

Operator used as base struct for OpContactMaterialVolOnSideLhs_dX_dx OpContactMaterialVolOnSideLhs_dX_dX operators that use side volume element adjacent to current contact prism needed to evaluate of deformation gradient tensor derivative.

Definition at line 594 of file SimpleContact.hpp.

Constructor & Destructor Documentation

◆ OpContactMaterialVolOnSideLhs()

SimpleContactProblem::OpContactMaterialVolOnSideLhs::OpContactMaterialVolOnSideLhs	(	const string	field_name_1,
		const string	field_name_2,
		boost::shared_ptr< CommonDataSimpleContact >	common_data_contact,
		const bool	is_master
	)

inline

LHS-operator for the contact element (material configuration)

Parameters

field_name_1	String of field name for spatial positions for rows
field_name_2	String of field name for spatial positions for columns
common_data_contact	Pointer to the common data for simple contact element
is_master	Bool parameter to distinguish between master and slave sides

Definition at line 644 of file SimpleContact.hpp.

         : MoFEM::VolumeElementForcesAndSourcesCoreOnContactPrismSide::
               UserDataOperator(field_name_1, field_name_2,
                                UserDataOperator::OPROWCOL),
           commonDataSimpleContact(common_data_contact), isMaster(is_master) {
       sYmm = false; // This will make sure to loop over all entities
       normalVector = boost::make_shared<VectorDouble>();
       tangentOne = boost::make_shared<VectorDouble>();
       tangentTwo = boost::make_shared<VectorDouble>();
     }

Member Function Documentation

◆ aSsemble()

MoFEMErrorCode SimpleContactProblem::OpContactMaterialVolOnSideLhs::aSsemble	(	EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

Definition at line 4402 of file SimpleContact.cpp.

                                           {
  
   MoFEMFunctionBegin;
  
   // get pointer to first global index on row
   const int *row_indices = &*row_data.getIndices().data().begin();
   // get pointer to first global index on column
   const int *col_indices = &*col_data.getIndices().data().begin();
  
   auto &data = *commonDataSimpleContact;
   if (!data.forcesOnlyOnEntitiesRow.empty()) {
     rowIndices.resize(row_nb_dofs, false);
     noalias(rowIndices) = row_data.getIndices();
     row_indices = &rowIndices[0];
     VectorDofs &dofs = row_data.getFieldDofs();
     VectorDofs::iterator dit = dofs.begin();
     for (int ii = 0; dit != dofs.end(); ++dit, ++ii) {
       if (data.forcesOnlyOnEntitiesRow.find((*dit)->getEnt()) ==
           data.forcesOnlyOnEntitiesRow.end()) {
         rowIndices[ii] = -1;
       }
     }
   }
  
   // assemble local matrix
   CHKERR MatSetValues(getSNESB(), row_nb_dofs, row_indices, col_nb_dofs,
                       col_indices, &*matLhs.data().begin(), ADD_VALUES);
  
   MoFEMFunctionReturn(0);
 }

◆ doWork()

MoFEMErrorCode SimpleContactProblem::OpContactMaterialVolOnSideLhs::doWork	(	int	row_side,
		int	col_side,
		EntityType	row_type,
		EntityType	col_type,
		EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

Definition at line 4356 of file SimpleContact.cpp.

                                           {
  
   MoFEMFunctionBegin;
   if (row_type != MBVERTEX)
     MoFEMFunctionReturnHot(0);
   row_nb_dofs = row_data.getIndices().size();
   if (!row_nb_dofs)
     MoFEMFunctionReturnHot(0);
   col_nb_dofs = col_data.getIndices().size();
   if (!col_nb_dofs)
     MoFEMFunctionReturnHot(0);
   nb_gauss_pts = row_data.getN().size1();
  
   nb_base_fun_row = row_data.getFieldData().size() / 3;
   nb_base_fun_col = col_data.getFieldData().size() / 3;
  
   matLhs.resize(3 * nb_base_fun_row, 3 * nb_base_fun_col, false);
   matLhs.clear();
  
   normalVector->resize(3, false);
   tangentOne->resize(3, false);
   tangentTwo->resize(3, false);
  
   if (isMaster) {
     normalVector = commonDataSimpleContact->normalVectorMasterPtr;
     tangentOne = commonDataSimpleContact->tangentOneVectorMasterPtr;
     tangentTwo = commonDataSimpleContact->tangentOneVectorMasterPtr;
     aRea = commonDataSimpleContact->areaMaster;
   } else {
     normalVector = commonDataSimpleContact->normalVectorSlavePtr;
     tangentOne = commonDataSimpleContact->tangentOneVectorSlavePtr;
     tangentTwo = commonDataSimpleContact->tangentOneVectorSlavePtr;
     aRea = commonDataSimpleContact->areaSlave;
   }
  
   // integrate local matrix for entity block
   CHKERR iNtegrate(row_data, col_data);
  
   // assemble local matrix
   CHKERR aSsemble(row_data, col_data);
  
   MoFEMFunctionReturn(0);
 }

◆ iNtegrate()

virtual MoFEMErrorCode SimpleContactProblem::OpContactMaterialVolOnSideLhs::iNtegrate	(	EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

inlinevirtual

Reimplemented in SimpleContactProblem::OpContactMaterialVolOnSideLhs_dX_dX.

Definition at line 623 of file SimpleContact.hpp.

                                                   {
       MoFEMFunctionBegin;
       MoFEMFunctionReturn(0);
     }

Member Data Documentation

◆ aRea

double SimpleContactProblem::OpContactMaterialVolOnSideLhs::aRea

Definition at line 616 of file SimpleContact.hpp.

◆ col_nb_dofs

int SimpleContactProblem::OpContactMaterialVolOnSideLhs::col_nb_dofs

Definition at line 606 of file SimpleContact.hpp.

◆ colIndices

VectorInt SimpleContactProblem::OpContactMaterialVolOnSideLhs::colIndices

Definition at line 603 of file SimpleContact.hpp.

◆ commonDataSimpleContact

boost::shared_ptr<CommonDataSimpleContact> SimpleContactProblem::OpContactMaterialVolOnSideLhs::commonDataSimpleContact

Definition at line 600 of file SimpleContact.hpp.

◆ isMaster

bool SimpleContactProblem::OpContactMaterialVolOnSideLhs::isMaster

private

Definition at line 659 of file SimpleContact.hpp.

◆ matLhs

MatrixDouble SimpleContactProblem::OpContactMaterialVolOnSideLhs::matLhs

Definition at line 598 of file SimpleContact.hpp.

◆ nb_base_fun_col

int SimpleContactProblem::OpContactMaterialVolOnSideLhs::nb_base_fun_col

Definition at line 610 of file SimpleContact.hpp.

◆ nb_base_fun_row

int SimpleContactProblem::OpContactMaterialVolOnSideLhs::nb_base_fun_row

Definition at line 609 of file SimpleContact.hpp.

◆ nb_gauss_pts

int SimpleContactProblem::OpContactMaterialVolOnSideLhs::nb_gauss_pts

Definition at line 607 of file SimpleContact.hpp.

◆ normalVector

boost::shared_ptr<VectorDouble> SimpleContactProblem::OpContactMaterialVolOnSideLhs::normalVector

Definition at line 615 of file SimpleContact.hpp.

◆ row_nb_dofs

int SimpleContactProblem::OpContactMaterialVolOnSideLhs::row_nb_dofs

Definition at line 605 of file SimpleContact.hpp.

◆ rowIndices

VectorInt SimpleContactProblem::OpContactMaterialVolOnSideLhs::rowIndices

Definition at line 602 of file SimpleContact.hpp.

◆ tangentOne

boost::shared_ptr<VectorDouble> SimpleContactProblem::OpContactMaterialVolOnSideLhs::tangentOne

Definition at line 612 of file SimpleContact.hpp.

◆ tangentTwo

boost::shared_ptr<VectorDouble> SimpleContactProblem::OpContactMaterialVolOnSideLhs::tangentTwo

Definition at line 613 of file SimpleContact.hpp.

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

users_modules/mortar_contact/src/SimpleContact.hpp
users_modules/mortar_contact/src/impl/SimpleContact.cpp

Public Member Functions

Public Attributes

Private Attributes