SimpleContactProblem::OpGapConstraintAugmentedOverSpatialMaster Struct Reference

LHS-operator for the simple contact element. More...

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

Inheritance diagram for SimpleContactProblem::OpGapConstraintAugmentedOverSpatialMaster:

Collaboration diagram for SimpleContactProblem::OpGapConstraintAugmentedOverSpatialMaster:

Public Member Functions
	OpGapConstraintAugmentedOverSpatialMaster (const string field_name, const string lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, const double cn)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)
	Integrates the conditions that fulfil KKT conditions at master face gauss points and assembles 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.

Integrates variation on the slave sid the conditions that fulfil KKT conditions with respect to Spatial positions on the master side and assembles components to LHS global matrix.

Definition at line 2051 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpGapConstraintAugmentedOverSpatialMaster()

SimpleContactProblem::OpGapConstraintAugmentedOverSpatialMaster::OpGapConstraintAugmentedOverSpatialMaster ( const string  field_name, const string  lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact >  common_data_contact, const double  cn  )

inline

Definition at line 2053 of file SimpleContact.hpp.

        : ContactOp(lagrange_field_name, field_name, UserDataOperator::OPROWCOL,
                    ContactOp::FACESLAVEMASTER),
          commonDataSimpleContact(common_data_contact), cN(cn) {
      sYmm = false; // This will make sure to loop over all entities (e.g.
                    // for order=2 it will make doWork to loop 16 time)
    }

Member Function Documentation

doWork()

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

Integrates the conditions that fulfil KKT conditions at master face gauss points and assembles components to LHS global matrix.

Integrates variation of the expresion that fulfils KKT conditions with respect to spatial positions in the integral sense and assembles components to LHS global matrix.

\[ {\text D}{\overline C(\lambda, \mathbf{x}^{(1)}, \delta \lambda)}[\Delta \mathbf{x}^{(2)}] = \left\{ \begin{array}{ll} \int_{{\gamma}^{(1)}_{\text c}} c_{\text n} \Delta \mathbf{x}^{(2)} \cdot {\mathbf{n}}_{\rm c} \delta{\lambda}\,\,{ {\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 and master triangles for \(i = 1\) and \(i = 2\), respectively. Furthermore, \( c_{\text n}\) works as an augmentation parameter and affects convergence, and \( g_{\textrm{n}}\) is the gap function evaluated at the slave triangle gauss points as:

\[ g_{\textrm{n}} = - \mathbf{n}(\mathbf{x}^{(1)}) \cdot \left( \mathbf{x}^{(1)} - \mathbf{x}^{(2)} \right) \]

Definition at line 1016 of file SimpleContact.cpp.

                                          {
  MoFEMFunctionBegin;
 
  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();
    int nb_base_fun_col = col_data.getFieldData().size() / 3;
 
    const double area_master =
        commonDataSimpleContact->areaSlave; // same area in master and slave
 
    const double area_slave =
        commonDataSimpleContact->areaSlave; // same area in master and slave
 
    NN.resize(nb_base_fun_row, 3 * nb_base_fun_col, false);
    NN.clear();
 
    auto get_tensor_from_vec = [](VectorDouble &n) {
      return FTensor::Tensor1<double *, 3>(&n(0), &n(1), &n(2));
    };
 
    FTensor::Index<'i', 3> i;
 
    auto t_const_unit_n =
        get_tensor_from_vec(*(commonDataSimpleContact->normalVectorSlavePtr));
 
    auto t_const_unit_master =
        get_tensor_from_vec(*(commonDataSimpleContact->normalVectorMasterPtr));
 
    auto t_aug_lambda_ptr =
        getFTensor0FromVec(*commonDataSimpleContact->augmentedLambdasPtr);
 
    auto t_w = getFTensor0IntegrationWeightSlave();
    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_slave;
 
      auto t_base_lambda = row_data.getFTensor0N(gg, 0);
      for (int bbr = 0; bbr != nb_base_fun_row; ++bbr) {
 
        auto t_base_master = col_data.getFTensor0N(gg, 0);
        FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3> t_mat{
            &NN(bbr, 0), &NN(bbr, 1), &NN(bbr, 2)};
        const double m = val_m * t_base_lambda;
 
        for (int bbc = 0; bbc != nb_base_fun_col; ++bbc) {
 
          t_mat(i) += t_const_unit_n(i) * m * t_base_master;
 
          ++t_base_master; // update rows
          ++t_mat;
        }
        ++t_base_lambda; // update cols master
      }
 
      ++t_aug_lambda_ptr;
      ++t_w;
    }
 
    CHKERR MatSetValues(getSNESB(), row_data, col_data, &*NN.data().begin(),
                        ADD_VALUES);
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

cN

const double SimpleContactProblem::OpGapConstraintAugmentedOverSpatialMaster::cN

private

Definition at line 2104 of file SimpleContact.hpp.

commonDataSimpleContact

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

private

Definition at line 2103 of file SimpleContact.hpp.

NN

MatrixDouble SimpleContactProblem::OpGapConstraintAugmentedOverSpatialMaster::NN

private

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