SimpleContactProblem::OpCalDerIntCompFunOverLambdaSlaveSlave Struct Reference

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

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

Inheritance diagram for SimpleContactProblem::OpCalDerIntCompFunOverLambdaSlaveSlave:

Collaboration diagram for SimpleContactProblem::OpCalDerIntCompFunOverLambdaSlaveSlave:

Public Member Functions
	OpCalDerIntCompFunOverLambdaSlaveSlave (const string lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, boost::shared_ptr< double > cn)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)
	Integrates the complementarity function at slave face gauss points and assembles components to LHS global matrix. More...

Private Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact
boost::shared_ptr< double >	cNPtr
MatrixDouble	NN

Detailed Description

LHS-operator for the simple contact element.

Integrates variation of the complementarity function with respect to Lagrange multipliers to fulfils KKT conditions in the integral sense on slave side and assembles components to LHS global matrix.

Definition at line 1386 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalDerIntCompFunOverLambdaSlaveSlave()

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

inline

Definition at line 1388 of file SimpleContact.hpp.

        : ContactOp(lagrange_field_name, UserDataOperator::OPROWCOL,
                    ContactOp::FACESLAVESLAVE),
          commonDataSimpleContact(common_data_contact), cNPtr(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::OpCalDerIntCompFunOverLambdaSlaveSlave::doWork	(	int	row_side,
		int	col_side,
		EntityType	row_type,
		EntityType	col_type,
		EntData &	row_data,
		EntData &	col_data
	)

Integrates the complementarity function at slave face gauss points and assembles components to LHS global matrix.

Integrates variation of the complementarity function with respect to Lagrange multipliers to fulfils KKT conditions in the integral sense nd assembles components to LHS global matrix.

\[ {\text D}{\overline C(\lambda, \mathbf{x}^{(i)}, \delta \lambda)}[\Delta \lambda] = \int_{{\gamma}^{(1)}_{\text c}} \Delta \lambda \left( 1 - {\text {sign}}\left( \lambda - c_{\text n} g_{\textrm{n}} \right) {\left| \lambda - c_{\text n} g_{\textrm{n}}\right|}^{r-1}\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 1992 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();
    const double area_slave =
        commonDataSimpleContact->areaSlave; // same area in master and slave
 
    NN.resize(nb_row, nb_col, false);
    NN.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 = SimpleContactProblem::ConstrainFunction_dl(
                               cn_value, t_gap_gp, t_lagrange_slave) *
                           t_w * area_slave;
 
      FTensor::Tensor0<FTensor::PackPtr<double *, 1>> t_mat(
          &*NN.data().begin());
 
      auto t_base_lambda_row = row_data.getFTensor0N(gg, 0);
      for (int bbr = 0; bbr != nb_row; ++bbr) {
        auto t_base_lambda_col = col_data.getFTensor0N(gg, 0);
        const double s = val_s * t_base_lambda_row;
        for (int bbc = 0; bbc != nb_col; ++bbc) {
 
          t_mat += s * t_base_lambda_col;
 
          ++t_mat;
          ++t_base_lambda_col; // update cols
        }
        ++t_base_lambda_row; // update rows
      }
 
      ++t_lagrange_slave;
      ++t_gap_gp;
      ++t_w;
    }
 
    CHKERR MatSetValues(getSNESB(), row_data, col_data, &*NN.data().begin(),
                        ADD_VALUES);
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

cNPtr

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

private

Definition at line 1435 of file SimpleContact.hpp.

commonDataSimpleContact

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

private

Definition at line 1434 of file SimpleContact.hpp.

NN

MatrixDouble SimpleContactProblem::OpCalDerIntCompFunOverLambdaSlaveSlave::NN

private

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