SimpleContactProblem::OpGetNormalSlaveALE Struct Reference

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

Inheritance diagram for SimpleContactProblem::OpGetNormalSlaveALE:

Collaboration diagram for SimpleContactProblem::OpGetNormalSlaveALE:

Public Member Functions
	OpGetNormalSlaveALE (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 slave surface vector based on material base coordinates. More...

Public Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact

Detailed Description

Definition at line 2645 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpGetNormalSlaveALE()

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

inline

Definition at line 2648 of file SimpleContact.hpp.

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

Member Function Documentation

doWork()

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

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

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

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

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

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

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

Definition at line 3337 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->normalVectorSlavePtr->resize(3, false);
 
  auto get_tensor_vec = [](VectorDouble &n) {
    return FTensor::Tensor1<double *, 3>(&n(0), &n(1), &n(2));
  };
 
  auto normal_original_slave =
      get_tensor_vec(*commonDataSimpleContact->normalVectorSlavePtr);
 
  commonDataSimpleContact->tangentOneVectorSlavePtr->resize(3, false);
  commonDataSimpleContact->tangentOneVectorSlavePtr->clear();
 
  commonDataSimpleContact->tangentTwoVectorSlavePtr->resize(3, false);
  commonDataSimpleContact->tangentTwoVectorSlavePtr->clear();
 
  auto tangent_0_slave =
      get_tensor_vec(*commonDataSimpleContact->tangentOneVectorSlavePtr);
  auto tangent_1_slave =
      get_tensor_vec(*commonDataSimpleContact->tangentTwoVectorSlavePtr);
 
  auto t_N = data.getFTensor1DiffN<2>(0, 0);
  auto t_dof = data.getFTensor1FieldData<3>();
 
  for (unsigned int dd = 0; dd != 3; ++dd) {
    tangent_0_slave(i) += t_dof(i) * t_N(0);
    tangent_1_slave(i) += t_dof(i) * t_N(1);
    ++t_dof;
    ++t_N;
  }
 
  normal_original_slave(i) =
      FTensor::levi_civita(i, j, k) * tangent_0_slave(j) * tangent_1_slave(k);
 
  const double normal_length =
      sqrt(normal_original_slave(i) * normal_original_slave(i));
  normal_original_slave(i) = normal_original_slave(i) / normal_length;
 
  commonDataSimpleContact->areaSlave = 0.5 * normal_length;
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataSimpleContact

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

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