SimpleContactProblem::OpLhsConvectIntegrationPtsConstrainMasterGap Struct Reference

Tangent opeerator for contrains equation for change of spatial positions on master and slave. More...

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

Inheritance diagram for SimpleContactProblem::OpLhsConvectIntegrationPtsConstrainMasterGap:

Collaboration diagram for SimpleContactProblem::OpLhsConvectIntegrationPtsConstrainMasterGap:

Public Member Functions
	OpLhsConvectIntegrationPtsConstrainMasterGap (const string lagrange_field_name, const string field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, boost::shared_ptr< double > cn, const ContactOp::FaceType face_type, boost::shared_ptr< MatrixDouble > diff_convect)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)

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

Detailed Description

Tangent opeerator for contrains equation for change of spatial positions on master and slave.

Definition at line 2427 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpLhsConvectIntegrationPtsConstrainMasterGap()

SimpleContactProblem::OpLhsConvectIntegrationPtsConstrainMasterGap::OpLhsConvectIntegrationPtsConstrainMasterGap ( const string  lagrange_field_name, const string  field_name, boost::shared_ptr< CommonDataSimpleContact >  common_data_contact, boost::shared_ptr< double >  cn, const ContactOp::FaceType  face_type, boost::shared_ptr< MatrixDouble >  diff_convect  )

inline

Definition at line 2429 of file SimpleContact.hpp.

        : ContactOp(lagrange_field_name, field_name, UserDataOperator::OPROWCOL,
                    face_type),
          commonDataSimpleContact(common_data_contact), cNPtr(cn),
          diffConvect(diff_convect) {
      sYmm = false;
    }

Member Function Documentation

doWork()

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

Definition at line 2987 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 && col_type == MBVERTEX) {
 
    const int nb_gauss_pts = row_data.getN().size1();
    int nb_base_fun_row = nb_row;
    int nb_base_fun_col = nb_col / 3;
 
    const double area_slave =
        commonDataSimpleContact->areaSlave; // same area in master and slave
 
    matLhs.resize(nb_row, nb_col, false);
    matLhs.clear();
 
    auto get_diff_ksi = [](auto &m, auto gg) {
      return FTensor::Tensor2<FTensor::PackPtr<double *, 1>, 2, 3>(
          &m(0, gg), &m(1, gg), &m(2, gg), &m(3, gg), &m(4, gg), &m(5, gg));
    };
 
    auto get_tensor_from_vec = [](VectorDouble &n) {
      return FTensor::Tensor1<double *, 3>(&n(0), &n(1), &n(2));
    };
 
    FTensor::Index<'i', 3> i;
    FTensor::Index<'j', 3> j;
    FTensor::Index<'I', 2> I;
 
    auto t_const_unit_n =
        get_tensor_from_vec(*(commonDataSimpleContact->normalVectorSlavePtr));
 
    auto t_lagrange_slave =
        getFTensor0FromVec(*commonDataSimpleContact->lagMultAtGaussPtsPtr);
    auto t_gap_gp = getFTensor0FromVec(*commonDataSimpleContact->gapPtr);
    auto &xi_grad_mat =
        *(commonDataSimpleContact->gradKsiPositionAtGaussPtsPtr);
    auto t_grad = getFTensor2FromMat<3, 2>(xi_grad_mat);
 
    auto t_w = getFTensor0IntegrationWeightSlave();
    const double cn_value = *cNPtr.get();
    for (int gg = 0; gg != nb_gauss_pts; ++gg) {
      const double val_m = SimpleContactProblem::ConstrainFunction(
                               cn_value, t_gap_gp, t_lagrange_slave) *
                           t_w * area_slave;
      const double val_diff_m_l = SimpleContactProblem::ConstrainFunction_dl(
                                      cn_value, t_gap_gp, t_lagrange_slave) *
                                  t_w * area_slave;
      const double val_diff_m_g = SimpleContactProblem::ConstrainFunction_dg(
                                      cn_value, t_gap_gp, t_lagrange_slave) *
                                  t_w * area_slave;
 
      auto t_base_lambda = row_data.getFTensor0N(gg, 0);
      auto t_base_diff_lambda = row_data.getFTensor1DiffN<2>(gg, 0);
 
      FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3> t_mat{
          &matLhs(0, 0), &matLhs(0, 1), &matLhs(0, 2)};
 
      for (int bbr = 0; bbr != nb_base_fun_row; ++bbr) {
 
        auto t_base_diff_disp = col_data.getFTensor1DiffN<2>(gg, 0);
        auto t_diff_convect = get_diff_ksi(*diffConvect, 3 * gg);
 
        for (int bbc = 0; bbc != nb_base_fun_col; ++bbc) {
 
          t_mat(i) += t_base_lambda * val_diff_m_g * t_const_unit_n(j) *
                      t_grad(j, I) * t_diff_convect(I, i);
 
          ++t_base_diff_disp;
          ++t_diff_convect;
          ++t_mat;
        }
 
        ++t_base_lambda;
        ++t_base_diff_lambda;
      }
 
      ++t_gap_gp;
      ++t_lagrange_slave;
      ++t_grad;
      ++t_w;
    }
 
    CHKERR MatSetValues(getSNESB(), row_data, col_data, &*matLhs.data().begin(),
                        ADD_VALUES);
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

cNPtr

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

private

Definition at line 2448 of file SimpleContact.hpp.

commonDataSimpleContact

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

private

Definition at line 2447 of file SimpleContact.hpp.

diffConvect

boost::shared_ptr<MatrixDouble> SimpleContactProblem::OpLhsConvectIntegrationPtsConstrainMasterGap::diffConvect

private

Definition at line 2449 of file SimpleContact.hpp.

matLhs

MatrixDouble SimpleContactProblem::OpLhsConvectIntegrationPtsConstrainMasterGap::matLhs

private

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