v0.14.0
Public Member Functions | Public Attributes | List of all members
SimpleContactProblem::OpCalculateDeformation Struct Reference

Operator for computing deformation gradients in side volumes. More...

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

Inheritance diagram for SimpleContactProblem::OpCalculateDeformation:
[legend]
Collaboration diagram for SimpleContactProblem::OpCalculateDeformation:
[legend]

Public Member Functions

MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 
 OpCalculateDeformation (const string field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, bool ho_geometry=false)
 

Public Attributes

boost::shared_ptr< CommonDataSimpleContactcommonDataSimpleContact
 
bool hoGeometry
 

Detailed Description

Operator for computing deformation gradients in side volumes.

Definition at line 2456 of file SimpleContact.hpp.

Constructor & Destructor Documentation

◆ OpCalculateDeformation()

SimpleContactProblem::OpCalculateDeformation::OpCalculateDeformation ( const string  field_name,
boost::shared_ptr< CommonDataSimpleContact common_data_contact,
bool  ho_geometry = false 
)
inline

Definition at line 2465 of file SimpleContact.hpp.

2471  commonDataSimpleContact(common_data_contact),
2472  hoGeometry(ho_geometry) {
2473  doEdges = false;
2474  doQuads = false;
2475  doTris = false;
2476  doTets = false;
2477  doPrisms = false;
2478  sYmm = false;
2479  };

Member Function Documentation

◆ doWork()

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

Definition at line 3081 of file SimpleContact.cpp.

3082  {
3083 
3088 
3089  // get number of integration points
3090  const int nb_integration_pts = getGaussPts().size2();
3091 
3092  auto t_h = getFTensor2FromMat<3, 3>(*commonDataSimpleContact->hMat);
3093  auto t_H = getFTensor2FromMat<3, 3>(*commonDataSimpleContact->HMat);
3094 
3095  commonDataSimpleContact->detHVec->resize(nb_integration_pts, false);
3096  commonDataSimpleContact->invHMat->resize(9, nb_integration_pts, false);
3097  commonDataSimpleContact->FMat->resize(9, nb_integration_pts, false);
3098 
3099  commonDataSimpleContact->detHVec->clear();
3100  commonDataSimpleContact->invHMat->clear();
3101  commonDataSimpleContact->FMat->clear();
3102 
3103  auto t_detH = getFTensor0FromVec(*commonDataSimpleContact->detHVec);
3104  auto t_invH = getFTensor2FromMat<3, 3>(*commonDataSimpleContact->invHMat);
3105  auto t_F = getFTensor2FromMat<3, 3>(*commonDataSimpleContact->FMat);
3106 
3107  for (int gg = 0; gg != nb_integration_pts; ++gg) {
3108  CHKERR determinantTensor3by3(t_H, t_detH);
3109  CHKERR invertTensor3by3(t_H, t_detH, t_invH);
3110  t_F(i, j) = t_h(i, k) * t_invH(k, j);
3111  ++t_h;
3112  ++t_H;
3113  ++t_detH;
3114  ++t_invH;
3115  ++t_F;
3116  }
3117 
3119 }

Member Data Documentation

◆ commonDataSimpleContact

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

Definition at line 2460 of file SimpleContact.hpp.

◆ hoGeometry

bool SimpleContactProblem::OpCalculateDeformation::hoGeometry

Definition at line 2461 of file SimpleContact.hpp.


The documentation for this struct was generated from the following files:
SimpleContactProblem::OpCalculateDeformation::hoGeometry
bool hoGeometry
Definition: SimpleContact.hpp:2461
MoFEM::invertTensor3by3
MoFEMErrorCode invertTensor3by3(ublas::matrix< T, L, A > &jac_data, ublas::vector< T, A > &det_data, ublas::matrix< T, L, A > &inv_jac_data)
Calculate inverse of tensor rank 2 at integration points.
Definition: Templates.hpp:1588
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:548
MoFEM::getFTensor0FromVec
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
Definition: Templates.hpp:135
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
field_name
constexpr auto field_name
Definition: poisson_2d_homogeneous.cpp:13
FTensor::Index< 'i', 3 >
MoFEM::determinantTensor3by3
static auto determinantTensor3by3(T &t)
Calculate the determinant of a 3x3 matrix or a tensor of rank 2.
Definition: Templates.hpp:1540
UserDataOperator
ForcesAndSourcesCore::UserDataOperator UserDataOperator
Definition: HookeElement.hpp:75
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
SimpleContactProblem::OpCalculateDeformation::commonDataSimpleContact
boost::shared_ptr< CommonDataSimpleContact > commonDataSimpleContact
Definition: SimpleContact.hpp:2460
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:429
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:359
MoFEM::ForcesAndSourcesCore::UserDataOperator::OPROW
@ OPROW
operator doWork function is executed on FE rows
Definition: ForcesAndSourcesCore.hpp:567