v0.16.0
Loading...
Searching...
No Matches
Public Types | Public Member Functions | List of all members
EshelbianPlasticity::OpCohesive_dJ_dkappa Struct Reference
Inheritance diagram for EshelbianPlasticity::OpCohesive_dJ_dkappa:
[legend]
Collaboration diagram for EshelbianPlasticity::OpCohesive_dJ_dkappa:
[legend]

Public Types

using OP = OpCohesiveRhs
 
- Public Types inherited from EshelbianPlasticity::OpCohesiveRhs
using OP = OpBrokenBaseCohesive
 
- Public Types inherited from EshelbianPlasticity::OpBrokenBaseCohesive
using OP = FormsIntegrators< FaceElementForcesAndSourcesCore::UserDataOperator >::Assembly< A >::OpBrokenBase
 

Public Member Functions

MoFEMErrorCode iNtegrate (EntData &data)
 
MoFEMErrorCode aSsemble (EntData &data)
 
 OpCohesiveRhs (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< double > gc_ptr, boost::shared_ptr< VectorDouble > kappa_ptr, boost::shared_ptr< VectorDouble > kappa_delta_ptr, boost::shared_ptr< std::array< MatrixDouble, 2 > > lambda_ptr=nullptr, Tag dissipation_tags=Tag(), Tag grad_dissipation_tags=Tag(), SmartPetscObj< Vec > vec_dJ_dx=SmartPetscObj< Vec >(), boost::shared_ptr< Range > ents_ptr=nullptr)
 
- Public Member Functions inherited from EshelbianPlasticity::OpCohesiveRhs
 OpCohesiveRhs (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< double > gc_ptr, boost::shared_ptr< VectorDouble > kappa_ptr, boost::shared_ptr< VectorDouble > kappa_delta_ptr, boost::shared_ptr< std::array< MatrixDouble, 2 > > lambda_ptr=nullptr, Tag dissipation_tags=Tag(), Tag grad_dissipation_tags=Tag(), SmartPetscObj< Vec > vec_dJ_dx=SmartPetscObj< Vec >(), boost::shared_ptr< Range > ents_ptr=nullptr)
 
MoFEMErrorCode iNtegrate (EntData &data)
 
- Public Member Functions inherited from EshelbianPlasticity::OpBrokenBaseCohesive
 OpBrokenBaseCohesive (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_flux_data_ptr, boost::shared_ptr< Range > ents_ptr=nullptr)
 
MoFEMErrorCode doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)
 

Additional Inherited Members

- Protected Attributes inherited from EshelbianPlasticity::OpCohesiveRhs
boost::shared_ptr< MatrixDouble > uGammaPtr
 
boost::shared_ptr< double > gcPtr
 
boost::shared_ptr< VectorDouble > kappaPtr
 
boost::shared_ptr< VectorDouble > kappaDeltaPtr
 
boost::shared_ptr< std::array< MatrixDouble, 2 > > lambdaPtr
 
boost::shared_ptr< double > totalDissipation
 
boost::shared_ptr< double > tatalDissipationGrad
 
SmartPetscObj< Vec > vec_dJdu
 
Tag dissipationTag
 
Tag gradDissipationTag
 
- Protected Attributes inherited from EshelbianPlasticity::OpBrokenBaseCohesive
boost::shared_ptr< MatrixDouble > fluxMatPtr
 
int faceSense = 0
 

Detailed Description

Definition at line 622 of file EshelbianCohesive.cpp.

Member Typedef Documentation

◆ OP

using EshelbianPlasticity::OpCohesive_dJ_dkappa::OP = OpCohesiveRhs

Definition at line 624 of file EshelbianCohesive.cpp.

Member Function Documentation

◆ aSsemble()

MoFEMErrorCode EshelbianPlasticity::OpCohesive_dJ_dkappa::aSsemble ( EntData &  data)
inline

Definition at line 723 of file EshelbianCohesive.cpp.

723 {
724 MoFEMFunctionBegin;
725 MoFEMFunctionReturn(0);
726 }
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

◆ iNtegrate()

MoFEMErrorCode EshelbianPlasticity::OpCohesive_dJ_dkappa::iNtegrate ( EntData &  data)
inline

Definition at line 627 of file EshelbianCohesive.cpp.

627 {
628 MoFEMFunctionBegin;
629
630 auto get_sense_index = [this]() { return (faceSense == 1) ? 0 : 1; };
631
632 FTENSOR_INDEX(3, i);
633 FTENSOR_INDEX(3, J);
634 FTENSOR_INDEX(3, K);
635
636 int nb_integration_pts = OP::getGaussPts().size2();
637
638 auto t_P = getFTensor2FromMat<3, 3>(*(OP::fluxMatPtr));
639 auto t_lambda = getFTensor2FromMat<3, 3>(lambdaPtr->at(get_sense_index()));
640 auto t_kappa = getFTensor0FromVec<0>(*kappaPtr);
641 auto t_delta_kappa = getFTensor0FromVec<0>(*kappaDeltaPtr);
642 auto t_face_normal = getFTensor1NormalsAtGaussPts();
643 auto t_w = OP::getFTensor0IntegrationWeight();
644
645 auto next = [&]() {
646 ++t_P;
647 ++t_face_normal;
648 ++t_kappa;
649 ++t_delta_kappa;
650 ++t_lambda;
651 ++t_w;
652 };
653
654 double face_dissipation = 0.0;
655 double face_grad_dissipation = 0.0;
656 auto face_handle = OP::getFEEntityHandle();
657 CHKERR getMoab().tag_get_data(dissipationTag, &face_handle, 1,
658 &face_dissipation);
659 CHKERR getMoab().tag_get_data(gradDissipationTag, &face_handle, 1,
660 &face_grad_dissipation);
661
662 for (auto gg = 0; gg != nb_integration_pts; ++gg) {
663 FTensor::Tensor1<double, 3> t_normal;
664 t_normal(J) = t_face_normal(J) * (faceSense / 2.);
665 FTensor::Tensor1<double, 3> t_normalized_normal;
666 t_normalized_normal(J) = t_normal(J);
667 t_normalized_normal.normalize();
668 FTensor::Tensor1<double, 3> t_traction;
669 t_traction(i) = t_P(i, J) * t_normalized_normal(J);
670
671 auto dJ_dkappa = [](auto &t_delta_kappa, auto &t_traction,
672 auto &t_normalized_normal, auto &t_kappa, auto gc) {
673 double kappa = static_cast<double>(t_kappa);
674 double delta_kappa = static_cast<double>(t_delta_kappa);
675 double teff = GriffithCohesiveLaw::calculateEffectiveTraction(
676 t_traction, t_normalized_normal, OpGetParameters::beta);
677 double m = GriffithCohesiveLaw::calculateDissipationSurplus(
678 delta_kappa, teff, kappa, gc, OpGetParameters::min_stiffness);
679 double m_grad =
680 GriffithCohesiveLaw::calculateDissipationSurplusDiffKappa(
681 delta_kappa, teff, kappa, gc, OpGetParameters::min_stiffness);
682 return boost::make_tuple(m, m_grad);
683 };
684
685 auto dr_kappa = [](auto &t_delta_kappa, auto &t_traction,
686 auto &t_normalized_normal, auto &t_kappa, auto gc) {
687 double kappa = static_cast<double>(t_kappa);
688 double delta_kappa = static_cast<double>(t_delta_kappa);
689 double kappa_plus_delta = kappa + delta_kappa;
690 double diff_alpha = GriffithCohesiveLaw::getDiffAlpha(
691 kappa_plus_delta, gc, OpGetParameters::min_stiffness);
692 auto [teff, t_gap] = GriffithCohesiveLaw::calculateGap(
693 t_traction, t_normalized_normal, diff_alpha, OpGetParameters::beta);
694 FTENSOR_INDEX(3, i);
695 FTensor::Tensor1<double, 3> t_gap_double;
696 t_gap_double(i) = -t_gap(i);
697 return t_gap_double;
698 };
699
700 auto [J, dJ] = dJ_dkappa(t_delta_kappa, t_traction, t_normalized_normal,
701 t_kappa, *gcPtr);
702 face_dissipation += t_w * J * t_normal.l2();
703 face_grad_dissipation += t_w * dJ * t_normal.l2();
704
705 auto t_dr = dr_kappa(t_delta_kappa, t_traction, t_normalized_normal,
706 t_kappa, *gcPtr);
707 FTensor::Tensor1<double, 3> t_l;
708 t_l(i) = t_lambda(i, K) * t_normal(K);
709 face_grad_dissipation -= t_w * t_l(i) * t_dr(i);
710
711 next();
712 }
713
714 CHKERR getMoab().tag_set_data(dissipationTag, &face_handle, 1,
715 &face_dissipation);
716 CHKERR getMoab().tag_set_data(gradDissipationTag, &face_handle, 1,
717 &face_grad_dissipation);
718
719
720 MoFEMFunctionReturn(0);
721 };
FTENSOR_INDEX
#define FTENSOR_INDEX(DIM, I)
Definition Templates.hpp:2673
FTensor::Tensor1::normalize
Tensor1< T, Tensor_Dim > normalize()
Definition Tensor1_value.hpp:77
CHKERR
#define CHKERR
Inline error check.
Definition definitions.h:551
kappa
double kappa
Definition free_surface.cpp:205
i
FTensor::Index< 'i', SPACE_DIM > i
Definition hcurl_divergence_operator_2d.cpp:27
J
FTensor::Index< 'J', DIM1 > J
Definition level_set.cpp:30
m
FTensor::Index< 'm', 3 > m
Definition shallow_wave.cpp:80
EshelbianPlasticity::GriffithCohesiveLaw::calculateDissipationSurplus
static auto calculateDissipationSurplus(const T &delta_kappa, const T t_eff, const T &kappa, double gc, double min_stiffness)
Definition EshelbianCohesive.cpp:109
EshelbianPlasticity::GriffithCohesiveLaw::calculateGap
static auto calculateGap(const FTensor::Tensor1< T, 3 > &t_traction, FTensor::Tensor1< double, 3 > &t_n_normalize, double alpha, double beta, bool sign_sensitive=true)
Definition EshelbianCohesive.cpp:155
EshelbianPlasticity::GriffithCohesiveLaw::calculateEffectiveTraction
static auto calculateEffectiveTraction(const FTensor::Tensor1< T, 3 > &t_traction, FTensor::Tensor1< double, 3 > &t_n_normalize, double beta)
Definition EshelbianCohesive.cpp:227
EshelbianPlasticity::GriffithCohesiveLaw::getDiffAlpha
static T getDiffAlpha(const T &kappa, double gc, double min_stiffness)
Definition EshelbianCohesive.cpp:51
EshelbianPlasticity::GriffithCohesiveLaw::calculateDissipationSurplusDiffKappa
static auto calculateDissipationSurplusDiffKappa(const T &delta_kappa, const T &t_eff, const T &kappa, double gc, double min_stiffness)
Definition EshelbianCohesive.cpp:118
EshelbianPlasticity::OpBrokenBaseCohesive::fluxMatPtr
boost::shared_ptr< MatrixDouble > fluxMatPtr
Definition EshelbianCohesive.cpp:414
EshelbianPlasticity::OpCohesiveRhs::gcPtr
boost::shared_ptr< double > gcPtr
Definition EshelbianCohesive.cpp:510
EshelbianPlasticity::OpCohesiveRhs::kappaDeltaPtr
boost::shared_ptr< VectorDouble > kappaDeltaPtr
Definition EshelbianCohesive.cpp:512
EshelbianPlasticity::OpCohesiveRhs::kappaPtr
boost::shared_ptr< VectorDouble > kappaPtr
Definition EshelbianCohesive.cpp:511
EshelbianPlasticity::OpCohesiveRhs::lambdaPtr
boost::shared_ptr< std::array< MatrixDouble, 2 > > lambdaPtr
Definition EshelbianCohesive.cpp:513

◆ OpCohesiveRhs()

EshelbianPlasticity::OpCohesiveRhs::OpCohesiveRhs ( boost::shared_ptr< std::vector< BrokenBaseSideData > >  broken_base_side_data,
boost::shared_ptr< double >  gc_ptr,
boost::shared_ptr< VectorDouble >  kappa_ptr,
boost::shared_ptr< VectorDouble >  kappa_delta_ptr,
boost::shared_ptr< std::array< MatrixDouble, 2 > >  lambda_ptr = nullptr,
Tag  dissipation_tags = Tag(),
Tag  grad_dissipation_tags = Tag(),
SmartPetscObj< Vec >  vec_dJ_dx = SmartPetscObj<Vec>(),
boost::shared_ptr< Range >  ents_ptr = nullptr 
)
inline

Definition at line 421 of file EshelbianCohesive.cpp.

430 : OpBrokenBaseCohesive(broken_base_side_data, ents_ptr), gcPtr(gc_ptr),
431 kappaPtr(kappa_ptr), kappaDeltaPtr(kappa_delta_ptr),
432 lambdaPtr(lambda_ptr), dissipationTag(dissipation_tags),
433 gradDissipationTag(grad_dissipation_tags), vec_dJdu(vec_dJ_dx) {}
EshelbianPlasticity::OpBrokenBaseCohesive::OpBrokenBaseCohesive
OpBrokenBaseCohesive(boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_flux_data_ptr, boost::shared_ptr< Range > ents_ptr=nullptr)
Definition EshelbianCohesive.cpp:338
The documentation for this struct was generated from the following file:
Generated by Doxygen 1.9.8 and hosted at University of Glasgow