RHS - operator for the contact element (material configuration) Integrates virtual work of contact traction in the material configuration on master surface. More...

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

Inheritance diagram for SimpleContactProblem::OpCalMatForcesALEMaster:

[legend]

Collaboration diagram for SimpleContactProblem::OpCalMatForcesALEMaster:

[legend]

Public Member Functions
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
	Integrates virtual work of contact traction in the material configuration. More...

MoFEMErrorCode	iNtegrate (EntData &row_data)

MoFEMErrorCode	aSsemble (EntData &row_data)

	OpCalMatForcesALEMaster (const string field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact)

Public Attributes
VectorInt	rowIndices

int	nbRows
	number of dofs on rows More...

int	nbIntegrationPts
	number of integration points More...

boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact

VectorDouble	vecF

Detailed Description

RHS - operator for the contact element (material configuration) Integrates virtual work of contact traction in the material configuration on master surface.

Definition at line 2514 of file SimpleContact.hpp.

Constructor & Destructor Documentation

◆ OpCalMatForcesALEMaster()

SimpleContactProblem::OpCalMatForcesALEMaster::OpCalMatForcesALEMaster	(	const string	field_name,
		boost::shared_ptr< CommonDataSimpleContact >	common_data_contact
	)

inline

Definition at line 2569 of file SimpleContact.hpp.

2572 : ContactOp(field_name, UserDataOperator::OPROW, ContactOp::FACEMASTER),

2573 commonDataSimpleContact(common_data_contact) {}

Member Function Documentation

◆ aSsemble()

MoFEMErrorCode SimpleContactProblem::OpCalMatForcesALEMaster::aSsemble ( EntData & row_data )

Definition at line 3218 of file SimpleContact.cpp.

                                                                        {
   MoFEMFunctionBegin;
  
   // get pointer to first global index on row
   const int *row_indices = &*row_data.getIndices().data().begin();
   auto &data = *commonDataSimpleContact;
   if (data.forcesOnlyOnEntitiesRow.empty())
     MoFEMFunctionReturnHot(0);
  
   rowIndices.resize(nbRows, 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;
     }
   }
  
   CHKERR VecSetValues(getSNESf(), nbRows, row_indices, &*vecF.data().begin(),
                       ADD_VALUES);
   MoFEMFunctionReturn(0);
 }

◆ doWork()

MoFEMErrorCode SimpleContactProblem::OpCalMatForcesALEMaster::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Integrates virtual work of contact traction in the material configuration.

Virtual work of the contact traction corresponding to a test function of the material configuration \((\delta{\mathbf{X}}^{(2)})\):

\[ \delta W^\text{material}_p({\mathbf{x}}^{(2)}, {\mathbf{X}}^{(2)}, \delta{\mathbf{X}}^{(2)}, \lambda) = -\int\limits_\mathcal{T^{\rm{(2)}}} \lambda \left\{{\left( \mathbf{F}^{(2)} \right)}^{\intercal}\cdot \mathbf{N}({\mathbf{X}}^{(2)}) \right\} \cdot \delta{\mathbf{X}}^{(2)}\, \textrm{d}\mathcal{T^{\rm {(2)}}} = -\int\limits_{\mathcal{T^{\rm{(2)}}}_{\xi}} \lambda \left\{{\left( \mathbf{F}^{(2)} \right)}^{\intercal}\cdot \left(\frac{\partial\mathbf{X}^{(2)}} {\partial\xi}\times\frac{\partial {\mathbf{X}}^{(2)}} {\partial\eta}\right) \right\} \cdot \delta{\mathbf{X}}^{(2)}\, \textrm{d}\xi\textrm{d}\eta \]

where \((2)\) denotes that variables are evaluated on master side, \( \lambda \) is contact traction on master surface, \({\mathbf{N}}({\mathbf{X}}^{(2)})\) is a normal to the face in the material configuration, \(\xi, \eta\) are coordinates in the parent space \(({\mathcal{T}}^{(2)}_\xi)\) and \(\mathbf{F}^{(2)}\) is the deformation gradient:

\[ \mathbf{F}^{(2)} = \mathbf{h}({\mathbf{x}}^{(2)})\,\mathbf{H}({\mathbf{X}}^{(2)})^{-1} = \frac{\partial{\mathbf{x}}^{(2)}}{\partial{\boldsymbol{\chi}}^{(2)}} \frac{\partial{\boldsymbol{\chi}}^{(2)}}{\partial{\mathbf{X}}^{(2)}} \]

where \(\mathbf{h}\) and \(\mathbf{H}\) are the gradients of the spatial and material maps, respectively, and \(\boldsymbol{\chi}\) are the reference coordinates.

Definition at line 3142 of file SimpleContact.cpp.

                                                                          {
  
   MoFEMFunctionBegin;
  
   // get number of dofs on row
   nbRows = row_data.getIndices().size();
   // if no dofs on row, exit that work, nothing to do here
   if (!nbRows)
     MoFEMFunctionReturnHot(0);
  
   vecF.resize(nbRows, false);
   vecF.clear();
  
   // get number of integration points
   nbIntegrationPts = getGaussPtsMaster().size2();
  
   // integrate local matrix for entity block
   CHKERR iNtegrate(row_data);
  
   // assemble local matrix
   CHKERR aSsemble(row_data);
  
   MoFEMFunctionReturn(0);
 }

◆ iNtegrate()

MoFEMErrorCode SimpleContactProblem::OpCalMatForcesALEMaster::iNtegrate ( EntData & row_data )

Definition at line 3169 of file SimpleContact.cpp.

                                                                     {
   MoFEMFunctionBegin;
  
   if (data.getIndices().size() == 0)
     MoFEMFunctionReturnHot(0);
   // in case the tet is not in database
   if (commonDataSimpleContact->FMat->size1() != 9)
     MoFEMFunctionReturnHot(0);
  
   const int nb_base_fun_col = nbRows / 3;
  
   auto get_tensor_vec = [](VectorDouble &n, const int r) {
     return FTensor::Tensor1<double *, 3>(&n(r + 0), &n(r + 1), &n(r + 2));
   };
  
   FTensor::Index<'i', 3> i;
   FTensor::Index<'j', 3> j;
   auto lagrange_slave =
       getFTensor0FromVec(*commonDataSimpleContact->lagMultAtGaussPtsPtr);
  
   auto t_F = getFTensor2FromMat<3, 3>(*commonDataSimpleContact->FMat);
   auto normal_at_gp =
       get_tensor_vec(*(commonDataSimpleContact->normalVectorMasterPtr), 0);
  
   auto t_w = getFTensor0IntegrationWeightMaster();
   double &area_master = commonDataSimpleContact->areaMaster;
   for (int gg = 0; gg != nbIntegrationPts; ++gg) {
     const double val_m = area_master * t_w;
     FTensor::Tensor0<double *> t_base_master(&data.getN()(gg, 0));
  
     for (int bbc = 0; bbc != nb_base_fun_col; ++bbc) {
  
       const double s = val_m * t_base_master * lagrange_slave;
  
       auto t_assemble_s = get_tensor_vec(vecF, 3 * bbc);
  
       t_assemble_s(i) -= s * t_F(j, i) * normal_at_gp(j);
  
       ++t_base_master;
     }
     ++t_F;
     ++lagrange_slave;
     ++t_w;
   } // for gauss points
  
   MoFEMFunctionReturn(0);
 }

Member Data Documentation

◆ commonDataSimpleContact

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

Definition at line 2521 of file SimpleContact.hpp.

◆ nbIntegrationPts

int SimpleContactProblem::OpCalMatForcesALEMaster::nbIntegrationPts

number of integration points

Definition at line 2519 of file SimpleContact.hpp.

◆ nbRows

int SimpleContactProblem::OpCalMatForcesALEMaster::nbRows

number of dofs on rows

Definition at line 2518 of file SimpleContact.hpp.

◆ rowIndices

VectorInt SimpleContactProblem::OpCalMatForcesALEMaster::rowIndices

Definition at line 2516 of file SimpleContact.hpp.

◆ vecF

VectorDouble SimpleContactProblem::OpCalMatForcesALEMaster::vecF

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