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

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

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

Public Member Functions

 OpHenckyTangentImpl (const std::string field_name, boost::shared_ptr< CommonData > common_data, boost::shared_ptr< MatrixDouble > mat_D_ptr=nullptr)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 
 OpHenckyTangentImpl (const std::string field_name, boost::shared_ptr< CommonData > common_data, boost::shared_ptr< MatrixDouble > mat_D_ptr=nullptr)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 

Private Attributes

boost::shared_ptr< CommonData > commonDataPtr
 
boost::shared_ptr< MatrixDouble > matDPtr
 

Detailed Description

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

Definition at line 495 of file HenckyOps.hpp.

Constructor & Destructor Documentation

◆ OpHenckyTangentImpl() [1/2]

template<int DIM, typename DomainEleOp , int S>
HenckyOps::OpHenckyTangentImpl< DIM, GAUSS, DomainEleOp, S >::OpHenckyTangentImpl ( const std::string  field_name,
boost::shared_ptr< CommonData >  common_data,
boost::shared_ptr< MatrixDouble >  mat_D_ptr = nullptr 
)
inline

Definition at line 496 of file HenckyOps.hpp.

499 : DomainEleOp(field_name, DomainEleOp::OPROW),
500 commonDataPtr(common_data) {
501 std::fill(&DomainEleOp::doEntities[MBEDGE],
502 &DomainEleOp::doEntities[MBMAXTYPE], false);
503 if (mat_D_ptr)
504 matDPtr = mat_D_ptr;
505 else
506 matDPtr = commonDataPtr->matDPtr;
507 }
DomainEleOp
DomainEle::UserDataOperator DomainEleOp
Finire element operator type.
Definition child_and_parent.cpp:37
field_name
constexpr auto field_name
Definition poisson_2d_homogeneous.cpp:13

◆ OpHenckyTangentImpl() [2/2]

template<int DIM, typename DomainEleOp , int S>
HenckyOps::OpHenckyTangentImpl< DIM, GAUSS, DomainEleOp, S >::OpHenckyTangentImpl ( const std::string  field_name,
boost::shared_ptr< CommonData >  common_data,
boost::shared_ptr< MatrixDouble >  mat_D_ptr = nullptr 
)
inline

Definition at line 651 of file HenckyOps.hpp.

654 : DomainEleOp(field_name, DomainEleOp::OPROW),
655 commonDataPtr(common_data) {
656 std::fill(&DomainEleOp::doEntities[MBEDGE],
657 &DomainEleOp::doEntities[MBMAXTYPE], false);
658 if (mat_D_ptr)
659 matDPtr = mat_D_ptr;
660 else
661 matDPtr = commonDataPtr->matDPtr;
662 }

Member Function Documentation

◆ doWork() [1/2]

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

Definition at line 509 of file HenckyOps.hpp.

509 {
510 MoFEMFunctionBegin;
511
512 FTensor::Index<'i', DIM> i;
513 FTensor::Index<'j', DIM> j;
514 FTensor::Index<'k', DIM> k;
515 FTensor::Index<'l', DIM> l;
516 FTensor::Index<'m', DIM> m;
517 FTensor::Index<'n', DIM> n;
518 FTensor::Index<'o', DIM> o;
519 FTensor::Index<'p', DIM> p;
520
521 [[maybe_unused]] constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
522 // const size_t nb_gauss_pts = matGradPtr->size2();
523 const size_t nb_gauss_pts = DomainEleOp::getGaussPts().size2();
524 commonDataPtr->matTangent.resize(nb_gauss_pts, DIM * DIM * DIM * DIM);
525 auto dP_dF =
526 getFTensor4FromMat<DIM, DIM, DIM, DIM>(commonDataPtr->matTangent);
527
528 auto t_D = getFTensor4DdgFromMat<DIM, DIM, S>(*matDPtr);
529 auto t_eig_val = getFTensor1FromMat<DIM>(commonDataPtr->matEigVal);
530 auto t_eig_vec = getFTensor2FromMat<DIM, DIM>(commonDataPtr->matEigVec);
531 auto t_T = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matHenckyStress);
532 auto t_S =
533 getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matSecondPiolaStress);
534 auto t_grad = getFTensor2FromMat<DIM, DIM>(*(commonDataPtr->matGradPtr));
535 auto t_logC_dC = getFTensor4DdgFromMat<DIM, DIM>(commonDataPtr->matLogCdC);
536
537 for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
538
539#ifdef HENCKY_SMALL_STRAIN
540 dP_dF(i, j, k, l) = t_D(i, j, k, l);
541#else
542
543 FTensor::Tensor2<double, DIM, DIM> t_F;
544 t_F(i, j) = t_grad(i, j) + t_kd(i, j);
545
546 // rare case when two eigen values are equal
547 auto nb_uniq = getUniqNb<DIM>(getVectorAdaptor(&t_eig_val(0), DIM),
548 FTensor::Number<DIM>{});
549 FTensor::Tensor4<double, DIM, DIM, DIM, DIM> dC_dF;
550 dC_dF(i, j, k, l) = (t_kd(i, l) * t_F(k, j)) + (t_kd(j, l) * t_F(k, i));
551
552 auto TL = EigenMatrix::getDiffDiffMat(t_eig_val, t_eig_vec, f, d_f, dd_f,
553 t_T, nb_uniq);
554 TL(i, j, k, l) *= 4;
555 FTensor::Ddg<double, DIM, DIM> P_D_P_plus_TL;
556 P_D_P_plus_TL(i, j, k, l) =
557 TL(i, j, k, l) +
558 (t_logC_dC(i, j, o, p) * t_D(o, p, m, n)) * t_logC_dC(m, n, k, l);
559 P_D_P_plus_TL(i, j, k, l) *= 0.5;
560 dP_dF(i, j, m, n) = t_kd(i, m) * (t_kd(k, n) * t_S(k, j));
561 dP_dF(i, j, m, n) +=
562 t_F(i, k) * (P_D_P_plus_TL(k, j, o, p) * dC_dF(o, p, m, n));
563
564#endif
565
566 ++dP_dF;
567
568 ++t_grad;
569 ++t_eig_val;
570 ++t_eig_vec;
571 ++t_logC_dC;
572 ++t_S;
573 ++t_T;
574 ++t_D;
575 }
576
577 MoFEMFunctionReturn(0);
578 }
FTensor::Kronecker_Delta
Kronecker Delta class.
Definition Kronecker_Delta.hpp:15
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition definitions.h:359
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition definitions.h:432
t_kd
constexpr auto t_kd
Definition free_surface.cpp:161
i
FTensor::Index< 'i', SPACE_DIM > i
Definition hcurl_divergence_operator_2d.cpp:27
n
const double n
refractive index of diffusive medium
Definition initial_diffusion.cpp:39
l
FTensor::Index< 'l', 3 > l
Definition matrix_function.cpp:23
j
FTensor::Index< 'j', 3 > j
Definition matrix_function.cpp:21
k
FTensor::Index< 'k', 3 > k
Definition matrix_function.cpp:22
EigenMatrix::getDiffDiffMat
auto getDiffDiffMat(A &&t_val, B &&t_vec, Fun< double > f, Fun< double > d_f, Fun< double > dd_f, C &&t_S, const int nb)
Get the Diff Diff Mat object.
Definition MatrixFunction.hpp:414
MoFEM::getVectorAdaptor
auto getVectorAdaptor(T1 ptr, const size_t n)
Get Vector adaptor.
Definition Templates.hpp:49
m
FTensor::Index< 'm', 3 > m
Definition shallow_wave.cpp:80

◆ doWork() [2/2]

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

Definition at line 664 of file HenckyOps.hpp.

664 {
665 MoFEMFunctionBegin;
666
667 FTensor::Index<'i', DIM> i;
668 FTensor::Index<'j', DIM> j;
669 FTensor::Index<'k', DIM> k;
670 FTensor::Index<'l', DIM> l;
671 FTensor::Index<'m', DIM> m;
672 FTensor::Index<'n', DIM> n;
673 FTensor::Index<'o', DIM> o;
674 FTensor::Index<'p', DIM> p;
675
676 constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
677 // const size_t nb_gauss_pts = matGradPtr->size2();
678 const size_t nb_gauss_pts = DomainEleOp::getGaussPts().size2();
679 commonDataPtr->matTangent.resize(nb_gauss_pts, DIM * DIM * DIM * DIM);
680 auto dP_dF =
681 getFTensor4FromMat<DIM, DIM, DIM, DIM>(commonDataPtr->matTangent);
682
683 auto t_D = getFTensor4DdgFromMat<
684 DIM, DIM, S, DataLayoutTraits<DataLayout::GaussByCoeffs>>(*matDPtr);
685 auto t_eig_val = getFTensor1FromMat<DIM>(commonDataPtr->matEigVal);
686 auto t_eig_vec = getFTensor2FromMat<DIM, DIM>(commonDataPtr->matEigVec);
687 auto t_T = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matHenckyStress);
688 auto t_S =
689 getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matSecondPiolaStress);
690 auto t_grad = getFTensor2FromMat<DIM, DIM>(*(commonDataPtr->matGradPtr));
691 auto t_logC_dC = getFTensor4DdgFromMat<DIM, DIM>(commonDataPtr->matLogCdC);
692 commonDataPtr->matCdF.resize(nb_gauss_pts, DIM * DIM * DIM * DIM);
693 auto dC_dF =
694 getFTensor4FromMat<DIM, DIM, DIM, DIM>(commonDataPtr->matCdF);
695
696 for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
697
698#ifdef HENCKY_SMALL_STRAIN
699 dP_dF(i, j, k, l) = t_D(i, j, k, l);
700#else
701
702 FTensor::Tensor2<double, DIM, DIM> t_F;
703 t_F(i, j) = t_grad(i, j) + t_kd(i, j);
704
705 // rare case when two eigen values are equal
706 auto nb_uniq = get_uniq_nb<DIM>(&t_eig_val(0));
707
708 dC_dF(i, j, k, l) = (t_kd(i, l) * t_F(k, j)) + (t_kd(j, l) * t_F(k, i));
709
710 auto TL = EigenMatrix::getDiffDiffMat(t_eig_val, t_eig_vec, f, d_f, dd_f,
711 t_T, nb_uniq);
712 TL(i, j, k, l) *= 4;
713 FTensor::Ddg<double, DIM, DIM> P_D_P_plus_TL;
714 P_D_P_plus_TL(i, j, k, l) =
715 TL(i, j, k, l) +
716 (t_logC_dC(i, j, o, p) * t_D(o, p, m, n)) * t_logC_dC(m, n, k, l);
717 P_D_P_plus_TL(i, j, k, l) *= 0.5;
718 dP_dF(i, j, m, n) = t_kd(i, m) * (t_kd(k, n) * t_S(k, j));
719 dP_dF(i, j, m, n) +=
720 t_F(i, k) * (P_D_P_plus_TL(k, j, o, p) * dC_dF(o, p, m, n));
721
722#endif
723
724 ++dP_dF;
725
726 ++t_grad;
727 ++t_eig_val;
728 ++t_eig_vec;
729 ++t_logC_dC;
730 ++t_S;
731 ++t_T;
732 ++t_D;
733 ++dC_dF;
734 }
735
736 MoFEMFunctionReturn(0);
737 }

Member Data Documentation

◆ commonDataPtr

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

Definition at line 581 of file HenckyOps.hpp.

◆ matDPtr

template<int DIM, typename DomainEleOp , int S>
boost::shared_ptr< MatrixDouble > HenckyOps::OpHenckyTangentImpl< DIM, GAUSS, DomainEleOp, S >::matDPtr
private

Definition at line 582 of file HenckyOps.hpp.

The documentation for this struct was generated from the following files:
Generated by Doxygen 1.9.8 and hosted at University of Glasgow