SimpleContactProblem::OpCalContactTractionOverLambdaSlaveSlave Struct Reference

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

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

Inheritance diagram for SimpleContactProblem::OpCalContactTractionOverLambdaSlaveSlave:

Collaboration diagram for SimpleContactProblem::OpCalContactTractionOverLambdaSlaveSlave:

Public Member Functions
	OpCalContactTractionOverLambdaSlaveSlave (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 and assembles Lagrange multipliers virtual work, \( \delta W_{\text c}\) derivative with respect to Lagrange multipliers with respect to Lagrange multipliers on slave 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 with respect to Lagrange multipliers on slave side side and assembles components to LHS global matrix.

Definition at line 1332 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpCalContactTractionOverLambdaSlaveSlave()

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

inline

Definition at line 1334 of file SimpleContact.hpp.

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

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

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

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

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

Definition at line 1927 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_slave =
        commonDataSimpleContact->areaSlave; // 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 = getFTensor0IntegrationWeightSlave();
 
    for (int gg = 0; gg != nb_gauss_pts; ++gg) {
 
      double val_m = t_w * area_slave;
      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::OpCalContactTractionOverLambdaSlaveSlave::commonDataSimpleContact

private

Definition at line 1373 of file SimpleContact.hpp.

NN

MatrixDouble SimpleContactProblem::OpCalContactTractionOverLambdaSlaveSlave::NN

private

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