SimpleContactProblem::OpCalContactTractionOverLambdaMasterSlave Struct Reference

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

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

Inheritance diagram for SimpleContactProblem::OpCalContactTractionOverLambdaMasterSlave:

Collaboration diagram for SimpleContactProblem::OpCalContactTractionOverLambdaMasterSlave:

Public Member Functions
	OpCalContactTractionOverLambdaMasterSlave (const string field_name, const string lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)
	Integrates Lagrange multipliers virtual work, \( \delta W_{\text c}\) derivative with respect to Lagrange multipliers with respect to Lagrange multipliers on master side and assembles components to LHS global matrix. More...

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

Detailed Description

LHS-operator for the simple contact element.

Integrates Lagrange multipliers virtual work, \( \delta W_{\text c}\) derivative with respect to Lagrange multipliers on master side and assembles components of the RHS vector.

Definition at line 1277 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalContactTractionOverLambdaMasterSlave()

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

inline

Definition at line 1279 of file SimpleContact.hpp.

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

Member Function Documentation

doWork()

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

Integrates Lagrange multipliers virtual work, \( \delta W_{\text c}\) derivative with respect to Lagrange multipliers with respect to Lagrange multipliers on master side and assembles components to LHS global matrix.

Computes linearisation of integrated on slave side complementarity function and assembles derivative of Lagrange multipliers virtual work \( \delta W_{\text c}\) with respect to Lagrange multipliers and assembles components to LHS global matrix

\[ {\text D} {\delta W^{(2)}_{\text c}(\lambda, \delta \mathbf{x}^{(2)}})[\Delta \lambda] \,\, = \int_{{\gamma}^{(1)}_{\text c}} \Delta \lambda \delta{\mathbf{x}^{(2)}} \,\,{ {\text d} {\gamma}} \]

where \({\gamma}^{(1)}_{\text c}\) is the surface integration domain of the slave surface, \( \lambda\) is the Lagrange multiplier, \(\mathbf{x}^{(2)}\) are the coordinates of the overlapping gauss points at master triangles.

Definition at line 1738 of file SimpleContact.cpp.

                                          {
  MoFEMFunctionBegin;
 
  // Both sides are needed since both sides contribute their shape
  // function to the stiffness matrix
  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() / 3;
    int nb_base_fun_col = col_data.getFieldData().size();
 
    const double area_master =
        commonDataSimpleContact->areaMaster; // same area in master and slave
 
    auto get_tensor_from_mat = [](MatrixDouble &m, const int r, const int c) {
      return FTensor::Tensor1<double *, 3>(&m(r + 0, c), &m(r + 1, c),
                                           &m(r + 2, c));
    };
 
    auto get_tensor_vec = [](VectorDouble &n) {
      return FTensor::Tensor1<double *, 3>(&n(0), &n(1), &n(2));
    };
 
    FTensor::Index<'i', 3> i;
 
    NN.resize(3 * nb_base_fun_row, nb_base_fun_col, false);
    NN.clear();
 
    auto const_unit_n =
        get_tensor_vec(*(commonDataSimpleContact->normalVectorSlavePtr.get()));
 
    auto t_w = getFTensor0IntegrationWeightMaster();
 
    for (int gg = 0; gg != nb_gauss_pts; ++gg) {
 
      double val_m = t_w * area_master;
      auto t_base_master = row_data.getFTensor0N(gg, 0);
 
      for (int bbr = 0; bbr != nb_base_fun_row; ++bbr) {
        auto t_assemble_m = get_tensor_from_mat(NN, 3 * bbr, 0);
        auto t_base_lambda = col_data.getFTensor0N(gg, 0);
        const double m = val_m * t_base_master;
        for (int bbc = 0; bbc != nb_base_fun_col; ++bbc) {
          const double n = m * t_base_lambda;
          t_assemble_m(i) -= n * const_unit_n(i);
          ++t_assemble_m;
          ++t_base_lambda; // update cols slave
        }
        ++t_base_master; // update rows master
      }
      ++t_w;
    }
 
    CHKERR MatSetValues(getSNESB(), row_data, col_data, &*NN.data().begin(),
                        ADD_VALUES);
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataSimpleContact

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

private

Definition at line 1319 of file SimpleContact.hpp.

NN

MatrixDouble SimpleContactProblem::OpCalContactTractionOverLambdaMasterSlave::NN

private

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