SimpleContactProblem::OpContactMaterialVolOnSideLhs_dX_dX Struct Reference

LHS-operator for the contact element (material configuration) 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 (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
	Compute part of the left-hand side. 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)
	LHS-operator for the contact element (material configuration) More...
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)
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 for the contact element (material configuration)

Computes linearisation of normal vector from the expression for material traction contribution with respect to material coordinates on slave or master side.

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

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 3441 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 ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

virtual

Compute part of the left-hand side.

Computes the linearisation of the material component with respect to a variation of material coordinates \((\Delta{\mathbf{X}}^{(i)})\):

\[ \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 \, \mathbf{F}^{\intercal}\cdot \left[ \frac{\partial{\mathbf{X}}^{(i)}} {\partial\xi} \cdot \left(\frac{\partial\Delta {\mathbf{X}}^{(i)}}{\partial\eta}\times\delta{\mathbf{x}}^{(i)}\right) -\frac{\partial{\mathbf{X}}^{(i)}} {\partial\eta} \cdot \left(\frac{\partial\Delta {\mathbf{X}}^{(i)}}{\partial\xi}\times \delta{\mathbf{x}}^{(i)}\right)\right] \textrm{d}\xi\textrm{d}\eta \]

Here superscript \((i)\) is either equal to 1 or 2 denoting slave or master side coordinates and surfaces, respectively. Moreover, \(\lambda\) is the lagrange multiplier.

Reimplemented from SimpleContactProblem::OpContactMaterialVolOnSideLhs.

Definition at line 4038 of file SimpleContact.cpp.

                                          {
 
  MoFEMFunctionBegin;
 
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
  FTensor::Index<'l', 3> l;
  FTensor::Index<'m', 3> m;
 
  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_h = getFTensor2FromMat<3, 3>(*commonDataSimpleContact->hMat);
  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) {
 
    const 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);
 
        t_assemble(i, j) -= a * t_row_base * t_inv_H(l, j) *
                            t_col_diff_base(m) * t_inv_H(m, i) * t_h(k, l) *
                            normal_at_gp(k);
 
        ++t_row_base;
      }
      ++t_col_diff_base;
    }
    ++t_w;
    ++t_h;
    ++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