SimpleContactProblem::OpCalAugmentedTractionRhsSlave Struct Reference

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

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

Inheritance diagram for SimpleContactProblem::OpCalAugmentedTractionRhsSlave:

Collaboration diagram for SimpleContactProblem::OpCalAugmentedTractionRhsSlave:

Public Member Functions
	OpCalAugmentedTractionRhsSlave (const string field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact)
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
	Integrates Lagrange multipliers virtual work on slave surface and assembles components to the RHS global vector. More...

Private Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact
VectorDouble	vecF

Detailed Description

RHS-operator for the simple contact element.

Integrates Augmented Lagrange multipliers virtual work on slave surface and assembles components to the RHS global vector.

Definition at line 1104 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalAugmentedTractionRhsSlave()

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

inline

Definition at line 1107 of file SimpleContact.hpp.

        : ContactPrismElementForcesAndSourcesCore::UserDataOperator(
              field_name, UserDataOperator::OPCOL,
              ContactPrismElementForcesAndSourcesCore::UserDataOperator::
                  FACESLAVE),
          commonDataSimpleContact(common_data_contact) {}

Member Function Documentation

doWork()

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

Integrates Lagrange multipliers virtual work on slave surface and assembles components to the RHS global vector.

Integrates Lagrange multipliers virtual work \( \delta W_{\text c}\) on slave surface and assembles components to the RHS global vector

\[ {\delta W^{(1)}_{\text c}(\lambda, \delta \mathbf{x}^{(1)}}) \,\, = - \left\{ \begin{array}{ll} \int_{{\gamma}^{(1)}_{\text c}} (\lambda + c_{\textrm n} g_{\textrm{n}}) \mathbf{n}(\mathbf{x}^{(1)}) \cdot \delta{\mathbf{x}^{(1)}} \,\,{ {\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, \( c_{\textrm n}\) is the regularisation/augmentation parameter of stress dimensions, \( g_{\textrm{n}}\) is the value of gap at the associated gauss point, \(\mathbf{x}^{(1)}\) are the coordinates of the overlapping gauss points at slave triangles.

Definition at line 1803 of file SimpleContact.cpp.

                                                               {
  MoFEMFunctionBegin;
 
  if (data.getIndices().size() == 0)
    MoFEMFunctionReturnHot(0);
 
  const int nb_dofs = data.getIndices().size();
  if (nb_dofs) {
    const int nb_gauss_pts = data.getN().size1();
    int nb_base_fun_col = data.getFieldData().size() / 3;
 
    vecF.resize(nb_dofs, false);
    vecF.clear();
 
    const double area_s =
        commonDataSimpleContact->areaSlave; // same area in master and slave
 
    auto get_tensor_vec = [](VectorDouble &n, const int r) {
      return FTensor::Tensor1<double *, 3>(&n(r + 0), &n(r + 1), &n(r + 2));
    };
 
    FTensor::Index<'i', 3> i;
 
    auto const_unit_n = get_tensor_vec(
        commonDataSimpleContact->normalVectorSlavePtr.get()[0], 0);
 
    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_aug_lambda_ptr;
        ++t_w;
        continue;
      }
      const double val_s = t_w * area_s;
 
      FTensor::Tensor0<double *> t_base_slave(&data.getN()(gg, 0));
 
      for (int bbc = 0; bbc != nb_base_fun_col; ++bbc) {
 
        auto t_assemble_s = get_tensor_vec(vecF, 3 * bbc);
 
        t_assemble_s(i) -=
            val_s * const_unit_n(i) * t_aug_lambda_ptr * t_base_slave;
 
        ++t_base_slave;
      }
      ++t_aug_lambda_ptr;
      ++t_w;
    } // for gauss points
 
    CHKERR VecSetValues(getSNESf(), data, &*vecF.begin(), ADD_VALUES);
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataSimpleContact

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

private

Definition at line 1148 of file SimpleContact.hpp.

vecF

VectorDouble SimpleContactProblem::OpCalAugmentedTractionRhsSlave::vecF

private

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