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

Collaboration diagram for SimpleContactProblem::OpCalContactAugmentedTractionOverLambdaMasterSlave:

Public Member Functions
	OpCalContactAugmentedTractionOverLambdaMasterSlave (const string field_name, const string lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)
	Integrates virtual work on master side , \( \delta W_{\text c}\), derivative with respect to Lagrange multipliers on slave side and assembles its components to LHS global matrix. More...

Private Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact
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 Lagrange multipliers on master side and assembles components of the LHS global matrix.

Definition at line 1577 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalContactAugmentedTractionOverLambdaMasterSlave()

SimpleContactProblem::OpCalContactAugmentedTractionOverLambdaMasterSlave::OpCalContactAugmentedTractionOverLambdaMasterSlave ( const string  field_name, const string  lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact >  common_data_contact  )

inline

Definition at line 1580 of file SimpleContact.hpp.

        : ContactOp(field_name, lagrange_field_name, UserDataOperator::OPROWCOL,
                    ContactOp::FACEMASTERSLAVE),
          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::OpCalContactAugmentedTractionOverLambdaMasterSlave::doWork	(	int	row_side,
		int	col_side,
		EntityType	row_type,
		EntityType	col_type,
		EntData &	row_data,
		EntData &	col_data
	)

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

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

\[ {\text D} {\delta W^{(2)}_{\text c}(\lambda, \delta \mathbf{x}^{(2)}})[\Delta \lambda] \,\, = \left\{ \begin{array}{ll} \int_{{\gamma}^{(1)}_{\text c}} \Delta \lambda \mathbf{n}(\mathbf{x}^{(1)}) \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 803 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();
  const int nb_col = col_data.getIndices().size();
 
  if (nb_row && nb_col) {
    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();
 
    const double area_master =
        commonDataSimpleContact->areaSlave; // same area in master and slave
 
    auto get_tensor_from_mat = [](MatrixDouble &m, const int r, const int c) {
      return FTensor::Tensor1<double *, 3>(&m(r + 0, c), &m(r + 1, c),
                                           &m(r + 2, c));
    };
 
    auto get_tensor_vec = [](VectorDouble &n) {
      return FTensor::Tensor1<double *, 3>(&n(0), &n(1), &n(2));
    };
 
    FTensor::Index<'i', 3> i;
 
    NN.resize(3 * nb_base_fun_row, nb_base_fun_col, false);
    NN.clear();
 
    auto const_unit_n =
        get_tensor_vec(*(commonDataSimpleContact->normalVectorSlavePtr.get()));
 
    auto t_w = getFTensor0IntegrationWeightSlave();
 
    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;
      }
      double val_m = t_w * area_master;
      auto t_base_master = row_data.getFTensor0N(gg, 0);
 
      for (int bbr = 0; bbr != nb_base_fun_row; ++bbr) {
        auto t_assemble_m = get_tensor_from_mat(NN, 3 * bbr, 0);
        auto t_base_lambda = col_data.getFTensor0N(gg, 0);
        const double m = val_m * t_base_master;
        for (int bbc = 0; bbc != nb_base_fun_col; ++bbc) {
          const double n = m * t_base_lambda;
          t_assemble_m(i) += n * const_unit_n(i);
          ++t_assemble_m;
          ++t_base_lambda; // update cols slave
        }
        ++t_base_master; // update rows master
      }
      ++t_w;
      ++t_aug_lambda_ptr;
    }
 
    CHKERR MatSetValues(getSNESB(), row_data, col_data, &*NN.data().begin(),
                        ADD_VALUES);
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataSimpleContact

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

private

Definition at line 1629 of file SimpleContact.hpp.

NN

MatrixDouble SimpleContactProblem::OpCalContactAugmentedTractionOverLambdaMasterSlave::NN

private

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