SimpleContactProblem::OpCalIntCompFunSlave Struct Reference

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

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

Inheritance diagram for SimpleContactProblem::OpCalIntCompFunSlave:

Collaboration diagram for SimpleContactProblem::OpCalIntCompFunSlave:

Public Member Functions
	OpCalIntCompFunSlave (const string lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, boost::shared_ptr< double > cn)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
	Integrates the complementarity function at slave face gauss points and assembles components to the RHS global vector. More...

Private Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact
boost::shared_ptr< double >	cNPtr
VectorDouble	vecR

Detailed Description

RHS-operator for the simple contact element.

Integrates complementarity function that fulfills KKT conditions over slave contact area and assembles components of the RHS vector.

Definition at line 1160 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalIntCompFunSlave()

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

inline

Definition at line 1162 of file SimpleContact.hpp.

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

Member Function Documentation

doWork()

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

Integrates the complementarity function at slave face gauss points and assembles components to the RHS global vector.

Integrates the complementarity function to fulfil KKT conditions in the integral sense and assembles components to the RHS global vector

\[ {\overline C(\lambda, \mathbf{x}^{(i)}, \delta \lambda)} = \int_{{\gamma}^{(1)}_{\text c}} \left( \lambda + c_{\text n} g_{\textrm{n}} - \dfrac{1}{r}{\left| \lambda - c_{\text n} g_{\textrm{n}}\right|}^{r}\right) \delta{{\lambda}} \,\,{ {\text d} {\gamma}} \]

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 1643 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();
  const double cn_value = *cNPtr.get();
  for (int gg = 0; gg != nb_gauss_pts; ++gg) {
    const double val_s = t_w * area_s *
                         SimpleContactProblem::ConstrainFunction(
                             cn_value, t_gap_gp, t_lagrange_slave);
    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_w;
  } // for gauss points
 
  CHKERR VecSetValues(getSNESf(), data, &vecR[0], ADD_VALUES);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

cNPtr

boost::shared_ptr<double> SimpleContactProblem::OpCalIntCompFunSlave::cNPtr

private

Definition at line 1201 of file SimpleContact.hpp.

commonDataSimpleContact

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

private

Definition at line 1200 of file SimpleContact.hpp.

vecR

VectorDouble SimpleContactProblem::OpCalIntCompFunSlave::vecR

private

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