SimpleContactProblem::OpGapConstraintAugmentedRhs Struct Reference

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

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

Inheritance diagram for SimpleContactProblem::OpGapConstraintAugmentedRhs:

Collaboration diagram for SimpleContactProblem::OpGapConstraintAugmentedRhs:

Public Member Functions
	OpGapConstraintAugmentedRhs (const string lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, const double cn)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact
	Integrates KKT conditions for Augmented Lagrangian formulation at slave face gauss points and assembles components to the RHS global vector. More...
const double	cN
VectorDouble	vecR

Detailed Description

RHS-operator for the simple contact element for Augmented Lagrangian Method.

Integrates ALM constraints that fulfill KKT conditions over slave contact area and assembles components of the RHS vector.

Definition at line 1214 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpGapConstraintAugmentedRhs()

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

inline

Definition at line 1217 of file SimpleContact.hpp.

        : ContactPrismElementForcesAndSourcesCore::UserDataOperator(
              lagrange_field_name, UserDataOperator::OPCOL,
              ContactOp::FACESLAVE),
          commonDataSimpleContact(common_data_contact), cN(cn) {}

Member Function Documentation

doWork()

MoFEMErrorCode SimpleContactProblem::OpGapConstraintAugmentedRhs::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Definition at line 1685 of file SimpleContact.cpp.

                                              {
  MoFEMFunctionBegin;
 
  if (data.getIndices().size() == 0)
    MoFEMFunctionReturnHot(0);
 
  const int nb_gauss_pts = data.getN().size1();
 
  int nb_base_fun_col = data.getFieldData().size();
  const double area_s =
      commonDataSimpleContact->areaSlave; // same area in master and slave
 
  vecR.resize(nb_base_fun_col, false);
  vecR.clear();
 
  auto t_lagrange_slave =
      getFTensor0FromVec(*commonDataSimpleContact->lagMultAtGaussPtsPtr);
  auto t_gap_gp = getFTensor0FromVec(*commonDataSimpleContact->gapPtr);
  auto t_w = getFTensor0IntegrationWeightSlave();
 
  auto t_aug_lambda_ptr =
      getFTensor0FromVec(*commonDataSimpleContact->augmentedLambdasPtr);
 
  for (int gg = 0; gg != nb_gauss_pts; ++gg) {
 
    double branch_gg;
    if (t_aug_lambda_ptr > 0 && std::abs(t_aug_lambda_ptr) > ALM_TOL) {
      branch_gg = -t_lagrange_slave / cN;
    } else {
      branch_gg = t_gap_gp;
    }
 
    const double val_s = t_w * area_s * branch_gg;
    auto t_base_lambda = data.getFTensor0N(gg, 0);
    for (int bbr = 0; bbr != nb_base_fun_col; ++bbr) {
      vecR[bbr] += val_s * t_base_lambda;
 
      ++t_base_lambda; // update rows
    }
 
    ++t_lagrange_slave;
    ++t_gap_gp;
    ++t_aug_lambda_ptr;
    ++t_w;
  } // for gauss points
 
  CHKERR VecSetValues(getSNESf(), data, &vecR[0], ADD_VALUES);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

cN

const double SimpleContactProblem::OpGapConstraintAugmentedRhs::cN

private

Definition at line 1265 of file SimpleContact.hpp.

commonDataSimpleContact

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

private

Integrates KKT conditions for Augmented Lagrangian formulation at slave face gauss points and assembles components to the RHS global vector.

Integrates the Augmented Lagrangian multipliers formulation to fulfil KKT conditions in the integral sense and assembles components to the RHS global vector

\[ {\overline C(\lambda, \mathbf{x}^{(i)}, \delta \lambda)} = \left\{ \begin{array}{ll} \int_{{\gamma}^{(1)}_{\text c}} g_{\textrm{n}} \delta{{\lambda}} \,\,{ {\text d} {\gamma}} & \lambda + c_{\text n} g_{\textrm{n}}\leq 0 \\ \int_{{\gamma}^{(1)}_{\text c}} -\dfrac{1}{c_{\text n}} \lambda \delta{{\lambda}} \,\,{ {\text d} {\gamma}} & \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}^{(i)}\) 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, \(r\) is regularisation parameter that can be chosen in \([1, 1.1]\) ( \(r = 1\)) is the default value) 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 1264 of file SimpleContact.hpp.

vecR

VectorDouble SimpleContactProblem::OpGapConstraintAugmentedRhs::vecR

private

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