v0.15.0
Loading...
Searching...
No Matches
PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp > Struct Template Reference

#include "tutorials/adv-0/src/PlasticOpsLargeStrains.hpp"

Inheritance diagram for PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >:
[legend]
Collaboration diagram for PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >:
[legend]

Public Member Functions

 OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl (const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< HenckyOps::CommonData > common_hencky_data_ptr, boost::shared_ptr< MatrixDouble > m_D_ptr)
 
MoFEMErrorCode iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
 

Private Attributes

boost::shared_ptr< CommonDatacommonDataPtr
 
boost::shared_ptr< HenckyOps::CommonDatacommonHenckyDataPtr
 
boost::shared_ptr< MatrixDouble > mDPtr
 
MatrixDouble locMat
 
MatrixDouble resDiff
 

Detailed Description

template<int DIM, typename AssemblyDomainEleOp>
struct PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >

Definition at line 8 of file PlasticOpsLargeStrains.hpp.

Constructor & Destructor Documentation

◆ OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl()

template<int DIM, typename AssemblyDomainEleOp >
PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl ( const std::string row_field_name,
const std::string col_field_name,
boost::shared_ptr< CommonData > common_data_ptr,
boost::shared_ptr< HenckyOps::CommonData > common_hencky_data_ptr,
boost::shared_ptr< MatrixDouble > m_D_ptr )

Definition at line 29 of file PlasticOpsLargeStrains.hpp.

35 : AssemblyDomainEleOp(row_field_name, col_field_name,
36 AssemblyDomainEleOp::OPROWCOL),
37 commonDataPtr(common_data_ptr),
38 commonHenckyDataPtr(common_hencky_data_ptr), mDPtr(m_D_ptr) {
39 AssemblyDomainEleOp::sYmm = false;
40}
FormsIntegrators< DomainEleOp >::Assembly< A >::OpBase AssemblyDomainEleOp

Member Function Documentation

◆ iNtegrate()

template<int DIM, typename AssemblyDomainEleOp >
MoFEMErrorCode PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::iNtegrate ( EntitiesFieldData::EntData & row_data,
EntitiesFieldData::EntData & col_data )

Definition at line 45 of file PlasticOpsLargeStrains.hpp.

46 {
48
49 FTensor::Index<'i', DIM> i;
50 FTensor::Index<'j', DIM> j;
51 FTensor::Index<'k', DIM> k;
52 FTensor::Index<'l', DIM> l;
53 FTensor::Index<'m', DIM> m;
54 FTensor::Index<'n', DIM> n;
55
56 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
58
59 auto &locMat = AssemblyDomainEleOp::locMat;
60
61 const size_t nb_integration_pts = row_data.getN().size1();
62 const size_t nb_row_base_functions = row_data.getN().size2();
63
64 if (AssemblyDomainEleOp::rowType == MBVERTEX &&
65 AssemblyDomainEleOp::rowSide == 0) {
66 resDiff.resize(DIM * DIM * size_symm, nb_integration_pts, false);
67 auto t_res_diff = getFTensor3FromMat<DIM, DIM, size_symm>(resDiff);
68 auto t_D = getFTensor4DdgFromMat<DIM, DIM, 0>(*mDPtr);
69 auto t_logC_dC =
70 getFTensor4DdgFromMat<DIM, DIM>(commonHenckyDataPtr->matLogCdC);
71 auto t_grad =
72 getFTensor2FromMat<DIM, DIM>(*(commonHenckyDataPtr->matGradPtr));
74 constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
75 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
77 t_F(i, j) = t_grad(i, j) + t_kd(i, j);
79 t_DLogC_dC(i, j, k, l) = t_D(m, n, k, l) * t_logC_dC(m, n, i, j);
81 t_FDLogC_dC(i, j, k, l) = t_F(i, m) * t_DLogC_dC(m, j, k, l);
83 t_res_diff(i, j, L) = t_FDLogC_dC(i, j, k, l) * t_L(k, l, L);
84 ++t_logC_dC;
85 ++t_grad;
86 ++t_res_diff;
87 }
88 }
89
90 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
91 auto t_row_diff_base = row_data.getFTensor1DiffN<DIM>();
92 auto t_res_diff = getFTensor3FromMat<DIM, DIM, size_symm>(resDiff);
93
94 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
95 double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
96
97 size_t rr = 0;
98 for (; rr != AssemblyDomainEleOp::nbRows / DIM; ++rr) {
99
100 auto t_mat =
102
104 t_tmp(i, L) = (t_res_diff(i, j, L) * (alpha * t_row_diff_base(j)));
105
106 auto t_col_base = col_data.getFTensor0N(gg, 0);
107 for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / size_symm; ++cc) {
108
109 t_mat(i, L) -= (t_col_base * t_tmp(i, L));
110
111 ++t_mat;
112 ++t_col_base;
113 }
114
115 ++t_row_diff_base;
116 }
117
118 for (; rr < nb_row_base_functions; ++rr)
119 ++t_row_diff_base;
120
121 ++t_w;
122 ++t_res_diff;
123 }
124
126}
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
const double n
refractive index of diffusive medium
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
auto symm_L_tensor(FTensor::Number< DIM >)
static FTensor::Tensor2< FTensor::PackPtr< double *, 3 >, 2, 3 > get_mat_vector_dtensor_sym(size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
constexpr auto size_symm
Definition plastic.cpp:42
FTensor::Index< 'm', 3 > m

Member Data Documentation

◆ commonDataPtr

template<int DIM, typename AssemblyDomainEleOp >
boost::shared_ptr<CommonData> PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::commonDataPtr
private

Definition at line 20 of file PlasticOpsLargeStrains.hpp.

◆ commonHenckyDataPtr

template<int DIM, typename AssemblyDomainEleOp >
boost::shared_ptr<HenckyOps::CommonData> PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::commonHenckyDataPtr
private

Definition at line 21 of file PlasticOpsLargeStrains.hpp.

◆ locMat

template<int DIM, typename AssemblyDomainEleOp >
MatrixDouble PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::locMat
private

Definition at line 23 of file PlasticOpsLargeStrains.hpp.

◆ mDPtr

template<int DIM, typename AssemblyDomainEleOp >
boost::shared_ptr<MatrixDouble> PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::mDPtr
private

Definition at line 22 of file PlasticOpsLargeStrains.hpp.

◆ resDiff

template<int DIM, typename AssemblyDomainEleOp >
MatrixDouble PlasticOps::OpCalculatePlasticInternalForceLhs_LogStrain_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::resDiff
private

Definition at line 24 of file PlasticOpsLargeStrains.hpp.


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