SimpleContactProblem::OpCalMatForcesALEMaster Struct Reference

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:

Collaboration diagram for SimpleContactProblem::OpCalMatForcesALEMaster:

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.

        : ContactOp(field_name, UserDataOperator::OPROW, ContactOp::FACEMASTER),
          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