v0.16.0
Loading...
Searching...
No Matches
Public Member Functions | Private Attributes | List of all members
HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S > Struct Template Reference

#include "tutorials/vec-2_nonlinear_elasticity/src/HenckyOps.hpp"

Inheritance diagram for HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >:
[legend]
Collaboration diagram for HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >:
[legend]

Public Member Functions

 OpCalculateHenckyStressImpl (const std::string field_name, boost::shared_ptr< CommonData > common_data)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 
 OpCalculateHenckyStressImpl (const std::string field_name, boost::shared_ptr< CommonData > common_data)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 

Private Attributes

boost::shared_ptr< CommonDatacommonDataPtr
 

Detailed Description

template<int DIM, typename DomainEleOp, int S>
struct HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >

Definition at line 298 of file HenckyOps.hpp.

Constructor & Destructor Documentation

◆ OpCalculateHenckyStressImpl() [1/2]

template<int DIM, typename DomainEleOp , int S>
HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >::OpCalculateHenckyStressImpl ( const std::string  field_name,
boost::shared_ptr< CommonData common_data 
)
inline

Definition at line 301 of file HenckyOps.hpp.

303 : DomainEleOp(field_name, DomainEleOp::OPROW),
304 commonDataPtr(common_data) {
305 std::fill(&DomainEleOp::doEntities[MBEDGE],
306 &DomainEleOp::doEntities[MBMAXTYPE], false);
307 }
DomainEle::UserDataOperator DomainEleOp
Finire element operator type.
constexpr auto field_name

◆ OpCalculateHenckyStressImpl() [2/2]

template<int DIM, typename DomainEleOp , int S>
HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >::OpCalculateHenckyStressImpl ( const std::string  field_name,
boost::shared_ptr< CommonData common_data 
)
inline

Definition at line 309 of file HenckyOps.hpp.

311 : DomainEleOp(field_name, DomainEleOp::OPROW),
312 commonDataPtr(common_data) {
313 std::fill(&DomainEleOp::doEntities[MBEDGE],
314 &DomainEleOp::doEntities[MBMAXTYPE], false);
315 }

Member Function Documentation

◆ doWork() [1/2]

template<int DIM, typename DomainEleOp , int S>
MoFEMErrorCode HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >::doWork ( int  side,
EntityType  type,
EntData data 
)
inline

Definition at line 309 of file HenckyOps.hpp.

309 {
311
312 FTensor::Index<'i', DIM> i;
313 FTensor::Index<'j', DIM> j;
314 FTensor::Index<'k', DIM> k;
315 FTensor::Index<'l', DIM> l;
316
317 // const size_t nb_gauss_pts = matGradPtr->size2();
318 const size_t nb_gauss_pts = DomainEleOp::getGaussPts().size2();
319 auto t_D = getFTensor4DdgFromMat<DIM, DIM, S>(*commonDataPtr->matDPtr);
320 auto t_logC = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matLogC);
321 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
322 commonDataPtr->matHenckyStress.resize(nb_gauss_pts, size_symm, false);
323 auto t_T = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matHenckyStress);
324
325 for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
326 t_T(i, j) = t_D(i, j, k, l) * t_logC(k, l);
327 ++t_logC;
328 ++t_T;
329 ++t_D;
330 }
331
333 }
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
constexpr auto size_symm
Definition plastic.cpp:42

◆ doWork() [2/2]

template<int DIM, typename DomainEleOp , int S>
MoFEMErrorCode HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >::doWork ( int  side,
EntityType  type,
EntData data 
)
inline

Definition at line 317 of file HenckyOps.hpp.

317 {
319
320 FTensor::Index<'i', DIM> i;
321 FTensor::Index<'j', DIM> j;
322 FTensor::Index<'k', DIM> k;
323 FTensor::Index<'l', DIM> l;
324
325 // const size_t nb_gauss_pts = matGradPtr->size2();
326 const size_t nb_gauss_pts = DomainEleOp::getGaussPts().size2();
327 auto t_D = getFTensor4DdgFromMat<DIM, DIM, S, HenckyTensorLayout<S>>(
328 *commonDataPtr->matDPtr);
329 auto t_logC = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matLogC);
330 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
331 commonDataPtr->matHenckyStress.resize(size_symm, nb_gauss_pts, false);
332 auto t_T = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matHenckyStress);
333
334 for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
335 t_T(i, j) = t_D(i, j, k, l) * t_logC(k, l);
336 ++t_logC;
337 ++t_T;
338 ++t_D;
339 }
340
342 }

Member Data Documentation

◆ commonDataPtr

template<int DIM, typename DomainEleOp , int S>
boost::shared_ptr< CommonData > HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >::commonDataPtr
private

Definition at line 336 of file HenckyOps.hpp.


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