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

Collaboration diagram for SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialSlaveSlave:

Public Member Functions
	OpCalContactAugmentedTractionOverSpatialSlaveSlave (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 slave spatial positions 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 Spatial position on slave side 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 1846 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalContactAugmentedTractionOverSpatialSlaveSlave()

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

inline

Definition at line 1849 of file SimpleContact.hpp.

        : ContactPrismElementForcesAndSourcesCore::UserDataOperator(
              field_name, field_name_2, UserDataOperator::OPROWCOL,
              ContactPrismElementForcesAndSourcesCore::UserDataOperator::
                  FACESLAVESLAVE),
          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::OpCalContactAugmentedTractionOverSpatialSlaveSlave::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 slave spatial positions and assembles its components to LHS global matrix.

Computes linearisation of virtual work on slave side integrated on the slave side and assembles the components of its derivative over Lagrange multipliers.

\[ {\text D} {\delta W^{(1)}_{\text c}(\lambda, \delta \mathbf{x}^{(1)}})[\Delta {\mathbf{x}^{(1)}}] \,\, = \left\{ \begin{array}{ll} \int_{{\gamma}^{(1)}_{\text c}} c_{\textrm n}\Delta {\mathbf{x}^{(1)}} \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}^{(1)}\) are the coordinates of the overlapping gauss points at 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 1337 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_slave_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_slave_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_slave_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::OpCalContactAugmentedTractionOverSpatialSlaveSlave::cN

private

Definition at line 1903 of file SimpleContact.hpp.

commonDataSimpleContact

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

private

Definition at line 1902 of file SimpleContact.hpp.

NN

MatrixDouble SimpleContactProblem::OpCalContactAugmentedTractionOverSpatialSlaveSlave::NN

private

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