v0.15.0
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Private Attributes | List of all members
OpFaceMaterialForce Struct Reference

#include "users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp"

Inheritance diagram for OpFaceMaterialForce:
[legend]
Collaboration diagram for OpFaceMaterialForce:
[legend]

Public Types

using OP = FaceElementForcesAndSourcesCore::UserDataOperator
 

Public Member Functions

 OpFaceMaterialForce (boost::shared_ptr< DataAtIntegrationPts > data_ptr)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 

Private Attributes

boost::shared_ptr< DataAtIntegrationPtsdataAtPts
 data at integration pts
 

Detailed Description

Examples
mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianPlasticity.cpp.

Definition at line 937 of file EshelbianOperators.hpp.

Member Typedef Documentation

◆ OP

using OpFaceMaterialForce::OP = FaceElementForcesAndSourcesCore::UserDataOperator

Definition at line 941 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

◆ OpFaceMaterialForce()

OpFaceMaterialForce::OpFaceMaterialForce ( boost::shared_ptr< DataAtIntegrationPts data_ptr)
inline

Definition at line 943 of file EshelbianOperators.hpp.

944 : OP(NOSPACE, OPSPACE), dataAtPts(data_ptr) {}
@ NOSPACE
Definition definitions.h:83
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
data at integration pts
FaceElementForcesAndSourcesCore::UserDataOperator OP

Member Function Documentation

◆ doWork()

MoFEMErrorCode OpFaceMaterialForce::doWork ( int  side,
EntityType  type,
EntData data 
)
Examples
mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 3507 of file EshelbianOperators.cpp.

3508 {
3510
3511#ifndef NDEBUG
3512 auto fe_mi_ptr = getFEMethod()->numeredEntFiniteElementPtr;
3513 auto pstatus = fe_mi_ptr->getPStatus();
3514 if (pstatus) {
3515 auto owner = fe_mi_ptr->getOwnerProc();
3516 MOFEM_LOG("SELF", Sev::noisy)
3517 << "OpFaceMaterialForce: owner proc is not 0, owner proc: " << owner
3518 << " " << getPtrFE()->mField.get_comm_rank();
3519 }
3520#endif // NDEBUG
3521
3523
3525 t_face_T(I) = 0.;
3526 double face_pressure = 0.;
3527 auto t_T = getFTensor1FromMat<SPACE_DIM>(
3528 dataAtPts->faceMaterialForceAtPts); //< face material force
3529 auto t_p =
3530 getFTensor0FromVec(dataAtPts->normalPressureAtPts); //< normal pressure
3531 auto t_w = getFTensor0IntegrationWeight();
3532 for (auto gg = 0; gg != getGaussPts().size2(); ++gg) {
3533 t_face_T(I) += t_w * t_T(I);
3534 face_pressure += t_w * t_p;
3535 ++t_w;
3536 ++t_T;
3537 ++t_p;
3538 }
3539 t_face_T(I) *= getMeasure();
3540 face_pressure *= getMeasure();
3541
3542 auto get_tag = [&](auto name, auto dim) {
3543 auto &moab = getPtrFE()->mField.get_moab();
3544 Tag tag;
3545 double def_val[] = {0., 0., 0.};
3546 CHK_MOAB_THROW(moab.tag_get_handle(name, dim, MB_TYPE_DOUBLE, tag,
3547 MB_TAG_CREAT | MB_TAG_SPARSE, def_val),
3548 "create tag");
3549 return tag;
3550 };
3551
3552 auto set_tag = [&](auto &&tag, auto ptr) {
3553 auto &moab = getPtrFE()->mField.get_moab();
3554 auto face = getPtrFE()->getFEEntityHandle();
3555 CHK_MOAB_THROW(moab.tag_set_data(tag, &face, 1, ptr), "set tag");
3556 };
3557
3558 set_tag(get_tag("MaterialForce", 3), &t_face_T(0));
3559 set_tag(get_tag("FacePressure", 1), &face_pressure);
3560
3562}
#define FTENSOR_INDEX(DIM, I)
constexpr int SPACE_DIM
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define CHK_MOAB_THROW(err, msg)
Check error code of MoAB function and throw MoFEM exception.
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MOFEM_LOG(channel, severity)
Log.
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
constexpr IntegrationType I

Member Data Documentation

◆ dataAtPts

boost::shared_ptr<DataAtIntegrationPts> OpFaceMaterialForce::dataAtPts
private

data at integration pts

Definition at line 950 of file EshelbianOperators.hpp.


The documentation for this struct was generated from the following files: