SimpleContactProblem::OpContactMaterialVolOnSideLhs_dX_dx Struct Reference

LHS-operator (material configuration) on the side volume of either master or slave side. More...

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

Inheritance diagram for SimpleContactProblem::OpContactMaterialVolOnSideLhs_dX_dx:

Collaboration diagram for SimpleContactProblem::OpContactMaterialVolOnSideLhs_dX_dx:

Public Member Functions
MoFEMErrorCode	iNtegrate (EntData &row_data, EntData &col_data)
	Integrates over a face contribution from a side volume. More...
	OpContactMaterialVolOnSideLhs_dX_dx (const string field_name_1, const string field_name_2, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, const bool is_master)
Public Member Functions inherited from SimpleContactProblem::OpContactMaterialVolOnSideLhs
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...

Additional Inherited Members
Public Attributes inherited from SimpleContactProblem::OpContactMaterialVolOnSideLhs
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

Detailed Description

LHS-operator (material configuration) on the side volume of either master or slave side.

Computes the linearisation of the material component with respect to a variation of material coordinates on the side volume for either master or slave side.

Definition at line 3341 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpContactMaterialVolOnSideLhs_dX_dx()

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

inline

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 3379 of file SimpleContact.hpp.

        : OpContactMaterialVolOnSideLhs(field_name_1, field_name_2,
                                        common_data_contact, is_master) {
      sYmm = false; // This will make sure to loop over all entities
    };

Member Function Documentation

iNtegrate()

MoFEMErrorCode SimpleContactProblem::OpContactMaterialVolOnSideLhs_dX_dx::iNtegrate	(	EntData &	row_data,
		EntData &	col_data
	)

Integrates over a face contribution from a side volume.

Computes linearisation of the material component with respect to a variation of material coordinates:

\[ \textrm{D} \delta W^\text{(i)}_{\rm{material}}({\mathbf{x}}^{(i)}, {\mathbf{X}}^{(i)}, \delta{\mathbf{x}}^{(i)}) [\Delta{\mathbf{X}}^{(i)}] = \int\limits_{{\mathcal{T}}^{(i)}_{\xi}} \lambda \left\{\left[ {\mathbf{h}({\mathbf{x}}^{(i)})}\,{\mathbf{H}^{-1}({\mathbf{X}}^{(i)})}\,\frac{\partial\Delta\mathbf{X}} {\partial{\boldsymbol{\chi}}^{(i)}}\,\mathbf{H}^{-1}({\mathbf{X}}^{(i)}) \right]^{\intercal}\cdot\left(\frac{\partial{\mathbf{X}}^{(i)}}{\partial\xi} \times\frac{\partial{\mathbf{X}}^{(i)}}{\partial\eta}\right)\right\} \cdot \delta{\mathbf{X}}^{(i)}\, \textrm{d}\xi\textrm{d}\eta \]

where \(i\) denoted either master or slave side.

Definition at line 3974 of file SimpleContact.cpp.

                                          {
 
  MoFEMFunctionBegin;
 
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
 
  auto get_tensor2 = [](MatrixDouble &m, const int r, const int c) {
    return FTensor::Tensor2<double *, 3, 3>(
        &m(r + 0, c + 0), &m(r + 0, c + 1), &m(r + 0, c + 2), &m(r + 1, c + 0),
        &m(r + 1, c + 1), &m(r + 1, c + 2), &m(r + 2, c + 0), &m(r + 2, c + 1),
        &m(r + 2, c + 2));
  };
 
  auto get_tensor_vec = [](VectorDouble &n) {
    return FTensor::Tensor1<double *, 3>(&n(0), &n(1), &n(2));
  };
 
  auto t_w = getFTensor0IntegrationWeight();
 
  auto t_inv_H = getFTensor2FromMat<3, 3>(*commonDataSimpleContact->invHMat);
 
  auto lagrange_slave =
      getFTensor0FromVec(*commonDataSimpleContact->lagMultAtGaussPtsPtr);
 
  auto normal_at_gp = get_tensor_vec(*normalVector);
 
  const double area_m = aRea;
 
  for (int gg = 0; gg != nb_gauss_pts; ++gg) {
 
    double a = -t_w * lagrange_slave * area_m;
 
    auto t_col_diff_base = col_data.getFTensor1DiffN<3>(gg, 0);
 
    int bbc = 0;
    for (; bbc != nb_base_fun_col; ++bbc) {
 
      FTensor::Tensor0<double *> t_row_base(&row_data.getN()(gg, 0));
 
      int bbr = 0;
      for (; bbr != nb_base_fun_row; ++bbr) {
 
        auto t_assemble = get_tensor2(matLhs, 3 * bbr, 3 * bbc);
        // TODO: handle hoGeometry
 
        t_assemble(i, j) += a * t_row_base * t_inv_H(k, i) *
                            t_col_diff_base(k) * normal_at_gp(j);
 
        ++t_row_base;
      }
      ++t_col_diff_base;
    }
    ++t_w;
    ++t_inv_H;
    ++lagrange_slave;
  }
 
  MoFEMFunctionReturn(0);
}

---

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