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

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

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

Public Member Functions

 OpCalculatePlasticFlowLhs_dUImpl (const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< CommonData > common_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< MatrixDouble > mDPtr
 

Detailed Description

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

Definition at line 117 of file PlasticOpsSmallStrains.hpp.

Constructor & Destructor Documentation

◆ OpCalculatePlasticFlowLhs_dUImpl()

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

Definition at line 132 of file PlasticOpsSmallStrains.hpp.

137 : AssemblyDomainEleOp(row_field_name, col_field_name,
138 AssemblyDomainEleOp::OPROWCOL),
139 commonDataPtr(common_data_ptr), mDPtr(m_D_ptr) {
140 AssemblyDomainEleOp::sYmm = false;
141}
FormsIntegrators< DomainEleOp >::Assembly< A >::OpBase AssemblyDomainEleOp

Member Function Documentation

◆ iNtegrate()

template<int DIM, typename AssemblyDomainEleOp >
MoFEMErrorCode PlasticOps::OpCalculatePlasticFlowLhs_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >::iNtegrate ( EntitiesFieldData::EntData & row_data,
EntitiesFieldData::EntData & col_data )
Examples
PlasticOpsSmallStrains.hpp.

Definition at line 145 of file PlasticOpsSmallStrains.hpp.

147 {
149
150 FTensor::Index<'i', DIM> i;
151 FTensor::Index<'j', DIM> j;
152 FTensor::Index<'k', DIM> k;
153 FTensor::Index<'l', DIM> l;
154 FTensor::Index<'m', DIM> m;
155 FTensor::Index<'n', DIM> n;
156
157 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
159
160 auto &locMat = AssemblyDomainEleOp::locMat;
161
162 const size_t nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
163 const size_t nb_row_base_functions = row_data.getN().size2();
164
165 auto t_res_flow_dstrain =
166 getFTensor4DdgFromMat<DIM, DIM>(commonDataPtr->resFlowDstrain);
167
168 auto next = [&]() { ++t_res_flow_dstrain; };
169
171 constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
173 t_diff_grad(i, j, k, l) = t_kd(i, k) * t_kd(j, l);
174
175 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
176 auto t_row_base = row_data.getFTensor0N();
177 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
178
179 double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
180 ++t_w;
182 t_res_tens(L, i, j) =
183 alpha * ((t_L(m, n, L) * (t_res_flow_dstrain(m, n, k, l))) *
184 t_diff_grad(k, l, i, j));
185 next();
186
187 size_t rr = 0;
188 for (; rr != AssemblyDomainEleOp::nbRows / size_symm; ++rr) {
189 auto t_mat =
191 auto t_col_diff_base = col_data.getFTensor1DiffN<DIM>(gg, 0);
192 for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / DIM; ++cc) {
193 t_mat(L, l) += t_row_base * (t_res_tens(L, l, k) * t_col_diff_base(k));
194 ++t_mat;
195 ++t_col_diff_base;
196 }
197 ++t_row_base;
198 }
199
200 for (; rr < nb_row_base_functions; ++rr)
201 ++t_row_base;
202 }
203
205}
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 auto get_mat_tensor_sym_dvector(size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
[Lambda functions]
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::OpCalculatePlasticFlowLhs_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >::commonDataPtr
private

Definition at line 127 of file PlasticOpsSmallStrains.hpp.

◆ mDPtr

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

Definition at line 128 of file PlasticOpsSmallStrains.hpp.


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