SimpleContactProblem::OpLhsConvectIntegrationPtsContactTraction Struct Reference

Calculate tangent operator for contact force for change of integration point positions, as result of change os spatial positions. More...

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

Inheritance diagram for SimpleContactProblem::OpLhsConvectIntegrationPtsContactTraction:

Collaboration diagram for SimpleContactProblem::OpLhsConvectIntegrationPtsContactTraction:

Public Member Functions
	OpLhsConvectIntegrationPtsContactTraction (const string row_field_name, const string col_field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, 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)

Public Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact

Private Attributes
MatrixDouble	matLhs
boost::shared_ptr< MatrixDouble >	diffConvect

Detailed Description

Calculate tangent operator for contact force for change of integration point positions, as result of change os spatial positions.

Definition at line 2361 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpLhsConvectIntegrationPtsContactTraction()

SimpleContactProblem::OpLhsConvectIntegrationPtsContactTraction::OpLhsConvectIntegrationPtsContactTraction ( const string  row_field_name, const string  col_field_name, boost::shared_ptr< CommonDataSimpleContact >  common_data_contact, const ContactOp::FaceType  face_type, boost::shared_ptr< MatrixDouble >  diff_convect  )

inline

Definition at line 2365 of file SimpleContact.hpp.

        : ContactOp(row_field_name, col_field_name, UserDataOperator::OPROWCOL,
                    face_type),
          commonDataSimpleContact(common_data_contact),
          diffConvect(diff_convect) {
      sYmm = false;
    }

Member Function Documentation

doWork()

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

Definition at line 2873 of file SimpleContact.cpp.

                                          {
  MoFEMFunctionBegin;
 
  const int nb_row_dofs = row_data.getIndices().size();
  const int nb_col_dofs = col_data.getIndices().size();
  if (nb_row_dofs && (nb_col_dofs && col_type == MBVERTEX)) {
 
    const int nb_gauss_pts = getGaussPtsSlave().size2();
    int nb_base_fun_row = nb_row_dofs / 3;
    int nb_base_fun_col = nb_col_dofs / 3;
 
    matLhs.resize(nb_row_dofs, nb_col_dofs, false);
    matLhs.clear();
 
    const double area_s =
        commonDataSimpleContact->areaSlave; // same area in master and slave
 
    FTensor::Index<'i', 3> i;
    FTensor::Index<'j', 3> j;
    FTensor::Index<'I', 2> I;
 
    auto get_tensor_vec = [](VectorDouble &n) {
      return FTensor::Tensor1<double *, 3>(&n[0], &n[1], &n[2]);
    };
 
    auto t_const_unit_n =
        get_tensor_vec(*(commonDataSimpleContact->normalVectorSlavePtr));
    auto t_diff_lambda_xi = getFTensor1FromMat<2>(
        *(commonDataSimpleContact->gradKsiLambdaAtGaussPtsPtr));
    auto t_w = getFTensor0IntegrationWeightSlave();
 
    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 t_base_row = row_data.getFTensor0N();
 
    for (int gg = 0; gg != nb_gauss_pts; ++gg) {
 
      double val_s = t_w * area_s;
      auto t_base_row = row_data.getFTensor0N(gg, 0);
 
      for (int rr = 0; rr != nb_base_fun_row; ++rr) {
 
        FTensor::Tensor2<FTensor::PackPtr<double *, 3>, 3, 3> t_mat{
            &matLhs(3 * rr + 0, 0), &matLhs(3 * rr + 0, 1),
            &matLhs(3 * rr + 0, 2),
 
            &matLhs(3 * rr + 1, 0), &matLhs(3 * rr + 1, 1),
            &matLhs(3 * rr + 1, 2),
 
            &matLhs(3 * rr + 2, 0), &matLhs(3 * rr + 2, 1),
            &matLhs(3 * rr + 2, 2)};
 
        auto t_diff_convect = get_diff_ksi(*diffConvect, 3 * gg);
 
        for (int cc = 0; cc != nb_base_fun_col; ++cc) {
          t_mat(i, j) -= val_s * t_base_row * t_const_unit_n(i) *
                         (t_diff_lambda_xi(I) * t_diff_convect(I, j));
 
          ++t_diff_convect;
          ++t_mat;
        }
 
        ++t_base_row;
      }
 
      ++t_diff_lambda_xi;
      ++t_w;
    } // for gauss points
 
    CHKERR MatSetValues(getSNESB(), row_data, col_data, &*matLhs.data().begin(),
                        ADD_VALUES);
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataSimpleContact

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

Definition at line 2363 of file SimpleContact.hpp.

diffConvect

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

private

Definition at line 2383 of file SimpleContact.hpp.

matLhs

MatrixDouble SimpleContactProblem::OpLhsConvectIntegrationPtsContactTraction::matLhs

private

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