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

#include "users_modules/multifield-thermoplasticity-private/tutorials/vec-2/src/HenckyOps.hpp"

Inheritance diagram for HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >:
[legend]
Collaboration diagram for HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >:
[legend]

Public Member Functions

 OpCalculateHenckyThermalStressdTImpl (const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< CommonData > elastic_common_data_ptr, boost::shared_ptr< VectorDouble > coeff_expansion_ptr)
 
MoFEMErrorCode iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
 
 OpCalculateHenckyThermalStressdTImpl (const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< CommonData > elastic_common_data_ptr, boost::shared_ptr< VectorDouble > coeff_expansion_ptr)
 
MoFEMErrorCode iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
 

Private Attributes

boost::shared_ptr< CommonDataelasticCommonDataPtr
 
boost::shared_ptr< VectorDouble > coeffExpansionPtr
 

Detailed Description

template<int DIM, typename AssemblyDomainEleOp, int S>
struct HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >

Definition at line 631 of file HenckyOps.hpp.

Constructor & Destructor Documentation

◆ OpCalculateHenckyThermalStressdTImpl() [1/2]

template<int DIM, typename AssemblyDomainEleOp , int S>
HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::OpCalculateHenckyThermalStressdTImpl ( const std::string  row_field_name,
const std::string  col_field_name,
boost::shared_ptr< CommonData elastic_common_data_ptr,
boost::shared_ptr< VectorDouble >  coeff_expansion_ptr 
)

Definition at line 647 of file HenckyOps.hpp.

◆ OpCalculateHenckyThermalStressdTImpl() [2/2]

template<int DIM, typename AssemblyDomainEleOp , int S>
HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::OpCalculateHenckyThermalStressdTImpl ( const std::string  row_field_name,
const std::string  col_field_name,
boost::shared_ptr< CommonData elastic_common_data_ptr,
boost::shared_ptr< VectorDouble >  coeff_expansion_ptr 
)

Member Function Documentation

◆ iNtegrate() [1/2]

template<int DIM, typename AssemblyDomainEleOp , int S>
MoFEMErrorCode HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::iNtegrate ( EntitiesFieldData::EntData &  row_data,
EntitiesFieldData::EntData &  col_data 
)

Definition at line 661 of file HenckyOps.hpp.

673 {
674
675#ifdef HENCKY_SMALL_STRAIN
676 dP_dF(i, j, k, l) = t_D(i, j, k, l);
677#else
678
680 t_F(i, j) = t_grad(i, j) + t_kd(i, j);
681
682 // rare case when two eigen values are equal
683 auto nb_uniq = get_uniq_nb<DIM>(&t_eig_val(0));
684
685 dC_dF(i, j, k, l) = (t_kd(i, l) * t_F(k, j)) + (t_kd(j, l) * t_F(k, i));
686
687 auto TL = EigenMatrix::getDiffDiffMat(t_eig_val, t_eig_vec, f, d_f, dd_f,
688 t_T, nb_uniq);
689 TL(i, j, k, l) *= 4;
690 FTensor::Ddg<double, DIM, DIM> P_D_P_plus_TL;
691 P_D_P_plus_TL(i, j, k, l) =
692 TL(i, j, k, l) +
693 (t_logC_dC(i, j, o, p) * t_D(o, p, m, n)) * t_logC_dC(m, n, k, l);
694 P_D_P_plus_TL(i, j, k, l) *= 0.5;
695 dP_dF(i, j, m, n) = t_kd(i, m) * (t_kd(k, n) * t_S(k, j));
696 dP_dF(i, j, m, n) +=
697 t_F(i, k) * (P_D_P_plus_TL(k, j, o, p) * dC_dF(o, p, m, n));
698
699#endif
700
701 ++dP_dF;
702
703 ++t_grad;
704 ++t_eig_val;
705 ++t_eig_vec;
706 ++t_logC_dC;
707 ++t_S;
708 ++t_T;
709 ++t_D;
710 ++dC_dF;
711 }
712
714 }
715
716private:
717 boost::shared_ptr<CommonData> commonDataPtr;
718 boost::shared_ptr<MatrixDouble> matDPtr;
719};
720
721template <int DIM, typename AssemblyDomainEleOp, int S>
723 : public AssemblyDomainEleOp {
724 OpCalculateHenckyThermalStressdTImpl(
725 const std::string row_field_name, const std::string col_field_name,
726 boost::shared_ptr<CommonData> elastic_common_data_ptr,
727 boost::shared_ptr<VectorDouble> coeff_expansion_ptr);
728
729 MoFEMErrorCode iNtegrate(EntitiesFieldData::EntData &row_data,
730 EntitiesFieldData::EntData &col_data);
731
732private:
733 boost::shared_ptr<CommonData> elasticCommonDataPtr;
734 boost::shared_ptr<VectorDouble> coeffExpansionPtr;
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
constexpr auto t_kd
@ GAUSS
Gaussian quadrature integration.
FTensor::Index< 'i', SPACE_DIM > i
const double n
refractive index of diffusive medium
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
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.
PetscErrorCode MoFEMErrorCode
MoFEM/PETSc error code.
FTensor::Index< 'm', 3 > m
OpCalculateHenckyThermalStressdTImpl(const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< CommonData > elastic_common_data_ptr, boost::shared_ptr< VectorDouble > coeff_expansion_ptr)

◆ iNtegrate() [2/2]

template<int DIM, typename AssemblyDomainEleOp , int S>
MoFEMErrorCode HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::iNtegrate ( EntitiesFieldData::EntData &  row_data,
EntitiesFieldData::EntData &  col_data 
)

Member Data Documentation

◆ coeffExpansionPtr

template<int DIM, typename AssemblyDomainEleOp , int S>
boost::shared_ptr< VectorDouble > HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::coeffExpansionPtr
private

Definition at line 643 of file HenckyOps.hpp.

◆ elasticCommonDataPtr

template<int DIM, typename AssemblyDomainEleOp , int S>
boost::shared_ptr< CommonData > HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::elasticCommonDataPtr
private

Definition at line 642 of file HenckyOps.hpp.


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