No Matches
Public Member Functions | Private Attributes | List of all members
EshelbianPlasticity::OpCalculateStrainEnergy Struct Reference

#include <users_modules/eshelbian_plasticty/src/EshelbianPlasticity.hpp>

Inheritance diagram for EshelbianPlasticity::OpCalculateStrainEnergy:
Collaboration diagram for EshelbianPlasticity::OpCalculateStrainEnergy:

Public Member Functions

 OpCalculateStrainEnergy (const std::string field_name, boost::shared_ptr< DataAtIntegrationPts > data_at_pts, boost::shared_ptr< double > &e)
MoFEMErrorCode doWork (int side, EntityType type, EntitiesFieldData::EntData &data)

Private Attributes

boost::shared_ptr< DataAtIntegrationPtsdataAtPts
boost::shared_ptr< doubleenergy

Detailed Description

Definition at line 807 of file EshelbianPlasticity.hpp.

Constructor & Destructor Documentation

◆ OpCalculateStrainEnergy()

EshelbianPlasticity::OpCalculateStrainEnergy::OpCalculateStrainEnergy ( const std::string  field_name,
boost::shared_ptr< DataAtIntegrationPts data_at_pts,
boost::shared_ptr< double > &  e 

Definition at line 809 of file EshelbianPlasticity.hpp.

813 dataAtPts(data_at_pts), energy(e) {}
VolumeElementForcesAndSourcesCore::UserDataOperator VolUserDataOperator
constexpr auto field_name
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
operator doWork function is executed on FE rows

Member Function Documentation

◆ doWork()

MoFEMErrorCode EshelbianPlasticity::OpCalculateStrainEnergy::doWork ( int  side,
EntityType  type,
EntitiesFieldData::EntData data 

Definition at line 2074 of file EshelbianOperators.cpp.

2075 {
2077 if (type == MBTET) {
2078 int nb_integration_pts = data.getN().size1();
2079 auto v = getVolume();
2080 auto t_w = getFTensor0IntegrationWeight();
2081 auto t_P = getFTensor2FromMat<3, 3>(dataAtPts->approxPAtPts);
2082 auto t_h = getFTensor2FromMat<3, 3>(dataAtPts->hAtPts);
2087 for (int gg = 0; gg != nb_integration_pts; ++gg) {
2088 const double a = t_w * v;
2089 // FIXME: this is wrong, energy should be calculated in material
2090 (*energy) += a * t_P(i, J) * t_h(i, J);
2091 ++t_w;
2092 ++t_P;
2093 ++t_h;
2094 }
2095 }
static Index< 'J', 3 > J
constexpr double a
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:346
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:416
FTensor::Index< 'i', SPACE_DIM > i
const double v
phase velocity of light in medium (cm/ns)
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....

Member Data Documentation

◆ dataAtPts

boost::shared_ptr<DataAtIntegrationPts> EshelbianPlasticity::OpCalculateStrainEnergy::dataAtPts

Definition at line 819 of file EshelbianPlasticity.hpp.

◆ energy

boost::shared_ptr<double> EshelbianPlasticity::OpCalculateStrainEnergy::energy

Definition at line 820 of file EshelbianPlasticity.hpp.

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