SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialMasterMaster Struct Reference

LHS-operator for the simple contact element with Augmented Lagrangian Method. More...

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

Inheritance diagram for SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialMasterMaster:

Collaboration diagram for SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialMasterMaster:

Public Member Functions
	OpCalContactAugmentedTractionOverSpatialMasterMaster (const string field_name, const string field_name_2, const double cn_value, boost::shared_ptr< CommonDataSimpleContact > common_data_contact)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
	Integrates virtual work on master side , \( \delta W_{\text c}\), derivative with respect to spatial positions of the master side and assembles its components to LHS global matrix. More...

Private Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact
const double	cN
MatrixDouble	NN

Detailed Description

LHS-operator for the simple contact element with Augmented Lagrangian Method.

Integrates Lagrange multipliers virtual work, \( \delta W_{\text c}\) derivative with respect to spatial positions on master side and assembles components of the LHS matrix.

Definition at line 1706 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalContactAugmentedTractionOverSpatialMasterMaster()

SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialMasterMaster::OpCalContactAugmentedTractionOverSpatialMasterMaster ( const string  field_name, const string  field_name_2, const double  cn_value, boost::shared_ptr< CommonDataSimpleContact >  common_data_contact  )

inline

Definition at line 1709 of file SimpleContact.hpp.

        : ContactPrismElementForcesAndSourcesCore::UserDataOperator(
              field_name, field_name_2, UserDataOperator::OPROWCOL,
              ContactPrismElementForcesAndSourcesCore::UserDataOperator::
                  FACEMASTERMASTER),
          cN(cn_value), commonDataSimpleContact(common_data_contact) {
      sYmm = false; // This will make sure to loop over all intities (e.g.
                    // for order=2 it will make doWork to loop 16 time)
    }

Member Function Documentation

doWork()

MoFEMErrorCode SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialMasterMaster::doWork	(	int	row_side,
		int	col_side,
		EntityType	row_type,
		EntityType	col_type,
		EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

Integrates virtual work on master side , \( \delta W_{\text c}\), derivative with respect to spatial positions of the master side and assembles its components to LHS global matrix.

Computes linearisation of integrated on slave side complementarity function and assembles Lagrange multipliers virtual work \( \delta W_{\text c}\) on master side with respect to spatial positions of the master side and assembles components to LHS global matrix

\[ {\text D} {\delta W^{(2)}_{\text c}(\lambda, \delta \mathbf{x}^{(2)}})[\Delta {\mathbf{x}^{(2)}}] \,\, = \left\{ \begin{array}{ll} \int_{{\gamma}^{(1)}_{\text c}} c_{\textrm n}\Delta {\mathbf{x}^{(2)}} \cdot [{\mathbf{n}}_{\rm c} \otimes {\mathbf{n}}_{\rm c}] \cdot \delta{\mathbf{x}^{(2)}}\,\,{ {\text d} {\gamma}} & \lambda + c_{\text n} g_{\textrm{n}}\leq 0 \\ 0 & \lambda + c_{\text n} g_{\textrm{n}}> 0 \\ \end{array} \right. \]

where \({\gamma}^{(1)}_{\text c}\) is the surface integration domain of the slave surface, \( \lambda\) is the Lagrange multiplier, \(\mathbf{x}^{(2)}\) are the coordinates of the overlapping gauss points at master triangles, \( c_{\textrm n}\) is the regularisation/augmentation parameter of stress dimensions and \( g_{\textrm{n}}\) is the gap evaluated on the corresponding slave side.

Definition at line 1168 of file SimpleContact.cpp.

                                               {
  MoFEMFunctionBegin;
 
  // Both sides are needed since both sides contribute their shape
  // function to the stiffness matrix
  const int nb_row = row_data.getIndices().size();
  if (!nb_row)
    MoFEMFunctionReturnHot(0);
  const int nb_col = col_data.getIndices().size();
  if (!nb_col)
    MoFEMFunctionReturnHot(0);
  const int nb_gauss_pts = row_data.getN().size1();
 
  int nb_base_fun_row = row_data.getFieldData().size() / 3;
  int nb_base_fun_col = col_data.getFieldData().size() / 3;
 
  auto get_tensor_from_mat = [](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));
  };
 
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
 
  NN.resize(3 * nb_base_fun_row, 3 * nb_base_fun_col, false);
  NN.clear();
 
  auto t_w = getFTensor0IntegrationWeightSlave();
 
  const double area_master = commonDataSimpleContact->areaSlave;
 
  auto normal =
      get_tensor_vec(commonDataSimpleContact->normalVectorSlavePtr.get()[0]);
 
  auto t_aug_lambda_ptr =
      getFTensor0FromVec(*commonDataSimpleContact->augmentedLambdasPtr);
 
  for (int gg = 0; gg != nb_gauss_pts; gg++) {
 
    if (t_aug_lambda_ptr > 0 && std::abs(t_aug_lambda_ptr) > ALM_TOL) {
      ++t_w;
      ++t_aug_lambda_ptr;
      continue;
    }
 
    const double val_m = t_w * area_master * cN;
 
    FTensor::Tensor0<double *> t_base_master_col(&col_data.getN()(gg, 0));
 
    for (int bbc = 0; bbc != nb_base_fun_col; bbc++) {
 
      FTensor::Tensor0<double *> t_base_master_row(&row_data.getN()(gg, 0));
      const double m = val_m * t_base_master_col;
 
      for (int bbr = 0; bbr != nb_base_fun_row; bbr++) {
 
        auto t_assemble_s = get_tensor_from_mat(NN, 3 * bbr, 3 * bbc);
 
        t_assemble_s(i, j) += m * normal(i) * normal(j) * t_base_master_row;
 
        ++t_base_master_row; // update rows
      }
      ++t_base_master_col; // update cols slave
    }
    ++t_w;
    ++t_aug_lambda_ptr;
  }
 
  CHKERR MatSetValues(getSNESB(), row_data, col_data, &*NN.data().begin(),
                      ADD_VALUES);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

cN

const double SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialMasterMaster::cN

private

Definition at line 1763 of file SimpleContact.hpp.

commonDataSimpleContact

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

private

Definition at line 1762 of file SimpleContact.hpp.

NN

MatrixDouble SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialMasterMaster::NN

private

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