SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialSlaveMaster 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::OpCalContactAugmentedTractionOverSpatialSlaveMaster:

Collaboration diagram for SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialSlaveMaster:

Public Member Functions
	OpCalContactAugmentedTractionOverSpatialSlaveMaster (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 slave side, \( \delta W_{\text c}\), derivative with respect to spatial positions on master side and assembles its components to the global LHS 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 virtual work of spatial position on slave side, \( \delta W_{\text c}\), derivative with respect to spatial positions on master side and assembles its components to the global LHS matrix.

Definition at line 1917 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalContactAugmentedTractionOverSpatialSlaveMaster()

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

inline

Definition at line 1920 of file SimpleContact.hpp.

        : ContactPrismElementForcesAndSourcesCore::UserDataOperator(
              field_name, field_name_2, UserDataOperator::OPROWCOL,
              ContactPrismElementForcesAndSourcesCore::UserDataOperator::
                  FACESLAVEMASTER),
          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::OpCalContactAugmentedTractionOverSpatialSlaveMaster::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 slave side, \( \delta W_{\text c}\), derivative with respect to spatial positions on master side and assembles its components to the global LHS matrix.

Computes linearisation of virtual work of spatial position on slave side , \( \delta W_{\text c}\), over master side spatial positions and assembles its components to LHS global matrix

\[ {\text D} {\delta W^{(1)}_{\text c}(\lambda, \Delta \mathbf{x}^{(2)}})[\delta {\mathbf{x}^{(1)}}] \,\, = \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}^{(1)}}\,\,{ {\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)}\) and \(\mathbf{x}^{(1)}\) are the coordinates of the overlapping gauss points at master and slave triangles, respectively. Also, \( 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 1420 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_slave = 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_s = t_w * area_slave * 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_slave_row(&row_data.getN()(gg, 0));
      const double s = val_s * 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) -= s * normal(i) * normal(j) * t_base_slave_row;
 
        ++t_base_slave_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::OpCalContactAugmentedTractionOverSpatialSlaveMaster::cN

private

Definition at line 1974 of file SimpleContact.hpp.

commonDataSimpleContact

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

private

Definition at line 1973 of file SimpleContact.hpp.

NN

MatrixDouble SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialSlaveMaster::NN

private

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