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

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

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

Public Member Functions

 OpCalculateConstraintsLhs_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::OpCalculateConstraintsLhs_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >

Definition at line 1064 of file PlasticOpsGeneric.hpp.

Constructor & Destructor Documentation

◆ OpCalculateConstraintsLhs_dUImpl()

template<int DIM, typename AssemblyDomainEleOp >
PlasticOps::OpCalculateConstraintsLhs_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >::OpCalculateConstraintsLhs_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 208 of file PlasticOpsSmallStrains.hpp.

213 : AssemblyDomainEleOp(row_field_name, col_field_name,
214 DomainEleOp::OPROWCOL),
215 commonDataPtr(common_data_ptr), mDPtr(m_D_ptr) {
216 AssemblyDomainEleOp::sYmm = false;
217}
FormsIntegrators< DomainEleOp >::Assembly< A >::OpBase AssemblyDomainEleOp

Member Function Documentation

◆ iNtegrate()

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

Definition at line 221 of file PlasticOpsSmallStrains.hpp.

223 {
225
226 FTensor::Index<'i', DIM> i;
227 FTensor::Index<'j', DIM> j;
228 FTensor::Index<'k', DIM> k;
229 FTensor::Index<'l', DIM> l;
230
231 auto &locMat = AssemblyDomainEleOp::locMat;
232
233 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
234 const auto nb_row_base_functions = row_data.getN().size2();
235
236 auto t_c_dstrain =
237 getFTensor2SymmetricFromMat<DIM>(commonDataPtr->resCdStrain);
238 auto t_diff_grad_symmetrise = diff_symmetrize(FTensor::Number<DIM>());
239
240 auto next = [&]() { ++t_c_dstrain; };
241
242 auto get_mat_scalar_dvector = [&]() {
243 if constexpr (DIM == 2)
244 return FTensor::Tensor1<FTensor::PackPtr<double *, 2>, 2>{&locMat(0, 0),
245 &locMat(0, 1)};
246 else
248 &locMat(0, 0), &locMat(0, 1), &locMat(0, 2)};
249 };
250
251 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
252 auto t_row_base = row_data.getFTensor0N();
253 for (auto gg = 0; gg != nb_integration_pts; ++gg) {
254 double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
255 ++t_w;
256
258 t_res_mat(i, j) =
259 ((t_c_dstrain(k, l)) * t_diff_grad_symmetrise(k, l, i, j));
260 next();
261
262 auto t_mat = get_mat_scalar_dvector();
263 size_t rr = 0;
264 for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
265 const double row_base = alpha * t_row_base;
266 auto t_col_diff_base = col_data.getFTensor1DiffN<DIM>(gg, 0);
267 for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / DIM; cc++) {
268 t_mat(i) += row_base * (t_res_mat(i, j) * t_col_diff_base(j));
269 ++t_mat;
270 ++t_col_diff_base;
271 }
272 ++t_row_base;
273 }
274 for (; rr != nb_row_base_functions; ++rr)
275 ++t_row_base;
276 }
277
279}
#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()
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
auto diff_symmetrize(FTensor::Number< DIM >)

Member Data Documentation

◆ commonDataPtr

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

Definition at line 1074 of file PlasticOpsGeneric.hpp.

◆ mDPtr

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

Definition at line 1075 of file PlasticOpsGeneric.hpp.


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