v0.15.0
Loading...
Searching...
No Matches
HenckyOps::OpCalculatePiolaStressImpl< DIM, GAUSS, DomainEleOp, S > Struct Template Reference

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

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

Public Member Functions

 OpCalculatePiolaStressImpl (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::OpCalculatePiolaStressImpl< DIM, GAUSS, DomainEleOp, S >

Definition at line 474 of file HenckyOps.hpp.

Constructor & Destructor Documentation

◆ OpCalculatePiolaStressImpl()

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

Definition at line 477 of file HenckyOps.hpp.

479 : DomainEleOp(field_name, DomainEleOp::OPROW),
480 commonDataPtr(common_data) {
481 std::fill(&DomainEleOp::doEntities[MBEDGE],
482 &DomainEleOp::doEntities[MBMAXTYPE], false);
483 }
DomainEle::UserDataOperator DomainEleOp
Finire element operator type.
constexpr auto field_name

Member Function Documentation

◆ doWork()

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

Definition at line 485 of file HenckyOps.hpp.

485 {
487
488 FTensor::Index<'i', DIM> i;
489 FTensor::Index<'j', DIM> j;
490 FTensor::Index<'k', DIM> k;
491 FTensor::Index<'l', DIM> l;
492
493 constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
494
495 // const size_t nb_gauss_pts = matGradPtr->size2();
496 const size_t nb_gauss_pts = DomainEleOp::getGaussPts().size2();
497#ifdef HENCKY_SMALL_STRAIN
498 auto t_D = getFTensor4DdgFromMat<DIM, DIM, S>(*commonDataPtr->matDPtr);
499#endif
500 auto t_logC = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matLogC);
501 auto t_logC_dC = getFTensor4DdgFromMat<DIM, DIM>(commonDataPtr->matLogCdC);
502 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
503 commonDataPtr->matFirstPiolaStress.resize(DIM * DIM, nb_gauss_pts, false);
504 commonDataPtr->matSecondPiolaStress.resize(size_symm, nb_gauss_pts, false);
505 auto t_P = getFTensor2FromMat<DIM, DIM>(commonDataPtr->matFirstPiolaStress);
506 auto t_T = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matHenckyStress);
507 auto t_S =
508 getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matSecondPiolaStress);
509 auto t_grad = getFTensor2FromMat<DIM, DIM>(*(commonDataPtr->matGradPtr));
510
511 for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
512
513#ifdef HENCKY_SMALL_STRAIN
514 t_P(i, j) = t_D(i, j, k, l) * t_grad(k, l);
515#else
517 t_F(i, j) = t_grad(i, j) + t_kd(i, j);
518 t_S(k, l) = t_T(i, j) * t_logC_dC(i, j, k, l);
519 t_P(i, l) = t_F(i, k) * t_S(k, l);
520#endif
521
522 ++t_grad;
523 ++t_logC;
524 ++t_logC_dC;
525 ++t_P;
526 ++t_T;
527 ++t_S;
528#ifdef HENCKY_SMALL_STRAIN
529 ++t_D;
530#endif
531 }
532
534 }
Kronecker Delta class.
#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()
constexpr auto t_kd
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

Member Data Documentation

◆ commonDataPtr

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

Definition at line 537 of file HenckyOps.hpp.


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