SimpleContactProblem::OpCalMatForcesALESlave Struct Reference

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

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

Inheritance diagram for SimpleContactProblem::OpCalMatForcesALESlave:

Collaboration diagram for SimpleContactProblem::OpCalMatForcesALESlave:

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)
	OpCalMatForcesALESlave (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 slave surface.

Definition at line 2581 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalMatForcesALESlave()

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

inline

Definition at line 2636 of file SimpleContact.hpp.

        : ContactOp(field_name, UserDataOperator::OPROW, ContactOp::FACESLAVE),
          commonDataSimpleContact(common_data_contact) {}

Member Function Documentation

aSsemble()

MoFEMErrorCode SimpleContactProblem::OpCalMatForcesALESlave::aSsemble

(

EntData &

row_data

)

Definition at line 3309 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::OpCalMatForcesALESlave::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}}^{(1)})\):

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

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

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

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 3245 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 = getGaussPtsSlave().size2();
 
  // integrate local matrix for entity block
  CHKERR iNtegrate(row_data);
 
  // assemble local matrix
  CHKERR aSsemble(row_data);
 
  MoFEMFunctionReturn(0);
}

iNtegrate()

MoFEMErrorCode SimpleContactProblem::OpCalMatForcesALESlave::iNtegrate

(

EntData &

row_data

)

Definition at line 3272 of file SimpleContact.cpp.

                                                                   {
  MoFEMFunctionBegin;
 
  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->normalVectorSlavePtr), 0);
  auto t_w = getFTensor0IntegrationWeightSlave();
  double &area_slave = commonDataSimpleContact->areaSlave;
  for (int gg = 0; gg != nbIntegrationPts; ++gg) {
    double val_s = t_w * area_slave;
    FTensor::Tensor0<double *> t_base_master(&data.getN()(gg, 0));
    for (int bbc = 0; bbc != nb_base_fun_col; ++bbc) {
      const double s = val_s * 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;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataSimpleContact

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

Definition at line 2588 of file SimpleContact.hpp.

nbIntegrationPts

int SimpleContactProblem::OpCalMatForcesALESlave::nbIntegrationPts

number of integration points

Definition at line 2586 of file SimpleContact.hpp.

nbRows

int SimpleContactProblem::OpCalMatForcesALESlave::nbRows

number of dofs on rows

Definition at line 2585 of file SimpleContact.hpp.

rowIndices

VectorInt SimpleContactProblem::OpCalMatForcesALESlave::rowIndices

Definition at line 2583 of file SimpleContact.hpp.

vecF

VectorDouble SimpleContactProblem::OpCalMatForcesALESlave::vecF

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