SimpleContactProblem::OpGetNormalMasterALE Struct Reference

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

Inheritance diagram for SimpleContactProblem::OpGetNormalMasterALE:

Collaboration diagram for SimpleContactProblem::OpGetNormalMasterALE:

Public Member Functions
	OpGetNormalMasterALE (const string field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
	Evaluates unit normal vector to the master surface vector based on material base coordinates. More...

Public Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact

Detailed Description

Definition at line 2689 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpGetNormalMasterALE()

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

inline

Definition at line 2692 of file SimpleContact.hpp.

        : ContactOp(field_name, UserDataOperator::OPCOL, ContactOp::FACEMASTER),
          commonDataSimpleContact(common_data_contact) {}

Member Function Documentation

doWork()

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

Evaluates unit normal vector to the master surface vector based on material base coordinates.

Computes normal vector based on material base coordinates based on mesh (moab vertices) coordinates:

\[ {\mathbf N}^{(2)}({\mathbf X}(\xi, \eta)) = \frac{\partial\mathbf{X}(\xi, \eta)}{\partial\xi}\times\frac{\partial \mathbf{X}(\xi, \eta)} {\partial\eta} \]

where \({\mathbf X}(\xi, \eta)\) is the vector of material coordinates at the gauss point on master surface with parent coordinates \(\xi\) and \(\eta\) evaluated according to

\[ {\mathbf X}(\xi, \eta) = \sum\limits^{3}_{i = 1} N_i(\xi, \eta){\bar{\mathbf X}}_i \]

where \( N_i \) is the shape function corresponding to the \( i-{\rm{th}}\) degree of freedom in the material configuration \({{\mathbf {\tilde X} }}_i\) corresponding to the 3 nodes of the triangular master face.

Definition at line 3394 of file SimpleContact.cpp.

                                                                  {
  MoFEMFunctionBegin;
 
  if (data.getFieldData().size() == 0)
    MoFEMFunctionReturnHot(0);
 
  if (type != MBVERTEX)
    MoFEMFunctionReturnHot(0);
 
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
 
  commonDataSimpleContact->normalVectorMasterPtr->resize(3, false);
 
  auto get_tensor_vec = [](VectorDouble &n) {
    return FTensor::Tensor1<double *, 3>(&n(0), &n(1), &n(2));
  };
 
  auto normal_original_master =
      get_tensor_vec(*commonDataSimpleContact->normalVectorMasterPtr);
 
  commonDataSimpleContact->tangentOneVectorMasterPtr->resize(3, false);
  commonDataSimpleContact->tangentOneVectorMasterPtr->clear();
 
  commonDataSimpleContact->tangentTwoVectorMasterPtr->resize(3, false);
  commonDataSimpleContact->tangentTwoVectorMasterPtr->clear();
 
  auto tangent_0_master =
      get_tensor_vec(*commonDataSimpleContact->tangentOneVectorMasterPtr);
  auto tangent_1_master =
      get_tensor_vec(*commonDataSimpleContact->tangentTwoVectorMasterPtr);
 
  auto t_N = data.getFTensor1DiffN<2>(0, 0);
  auto t_dof = data.getFTensor1FieldData<3>();
 
  for (unsigned int dd = 0; dd != 3; ++dd) {
    tangent_0_master(i) += t_dof(i) * t_N(0);
    tangent_1_master(i) += t_dof(i) * t_N(1);
    ++t_dof;
    ++t_N;
  }
 
  normal_original_master(i) =
      FTensor::levi_civita(i, j, k) * tangent_0_master(j) * tangent_1_master(k);
 
  const double normal_length =
      sqrt(normal_original_master(i) * normal_original_master(i));
  normal_original_master(i) = normal_original_master(i) / normal_length;
 
  commonDataSimpleContact->areaMaster = 0.5 * normal_length;
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataSimpleContact

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

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