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

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

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

Public Member Functions

 OpCalculateConstraintsLhs_LogStrain_dUImpl (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 resDiff
 

Detailed Description

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

Definition at line 256 of file PlasticOpsLargeStrains.hpp.

Constructor & Destructor Documentation

◆ OpCalculateConstraintsLhs_LogStrain_dUImpl()

template<int DIM, typename AssemblyDomainEleOp >
PlasticOps::OpCalculateConstraintsLhs_LogStrain_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >::OpCalculateConstraintsLhs_LogStrain_dUImpl ( 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 274 of file PlasticOpsLargeStrains.hpp.

280 : AssemblyDomainEleOp(row_field_name, col_field_name,
281 DomainEleOp::OPROWCOL),
282 commonDataPtr(common_data_ptr),
283 commonHenckyDataPtr(common_hencky_data_ptr), mDPtr(m_D_ptr) {
284 AssemblyDomainEleOp::sYmm = false;
285}
FormsIntegrators< DomainEleOp >::Assembly< A >::OpBase AssemblyDomainEleOp

Member Function Documentation

◆ iNtegrate()

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

Definition at line 289 of file PlasticOpsLargeStrains.hpp.

291 {
293
294 FTensor::Index<'i', DIM> i;
295 FTensor::Index<'j', DIM> j;
296 FTensor::Index<'k', DIM> k;
297 FTensor::Index<'l', DIM> l;
298
299 constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
300
301 AssemblyDomainEleOp::locMat.resize(AssemblyDomainEleOp::nbRows, AssemblyDomainEleOp::nbCols,
302 false);
303 AssemblyDomainEleOp::locMat.clear();
304
305 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
306 const auto nb_row_base_functions = row_data.getN().size2();
307
308 if (AssemblyDomainEleOp::colType == MBVERTEX &&
309 AssemblyDomainEleOp::colSide == 0) {
310
311 resDiff.resize(DIM * DIM, nb_integration_pts, false);
312 auto t_res_diff = getFTensor2FromMat<DIM, DIM>(resDiff);
313 auto t_grad = getFTensor2FromMat<DIM, DIM>(*(commonDataPtr->mGradPtr));
314 auto t_logC_dC =
315 getFTensor4DdgFromMat<DIM, DIM>(commonHenckyDataPtr->matLogCdC);
316 auto t_c_dstrain =
317 getFTensor2SymmetricFromMat<DIM>(commonDataPtr->resCdStrain);
318
319 auto next = [&]() {
320 ++t_grad;
321 ++t_logC_dC;
322 ++t_c_dstrain;
323 ++t_res_diff;
324 };
325
326 auto t_diff_grad_symmetrise = diff_symmetrize(FTensor::Number<DIM>());
327
328 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
332
333 t_diff_ls_dlog_c(k, l) =
334 (t_c_dstrain(i, j)) * (t_logC_dC(i, j, k, l) / 2);
335 t_F(i, j) = t_grad(i, j) + t_kd(i, j);
336 t_dC_dF(i, j, k, l) = (t_kd(i, l) * t_F(k, j)) + (t_kd(j, l) * t_F(k, i));
337
339 t_res_diff(i, j) = (t_diff_ls_dlog_c(k, l) * t_dC_dF(k, l, i, j));
340 next();
341 }
342 }
343
344 auto t_res_diff = getFTensor2FromMat<DIM, DIM>(resDiff);
345
346 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
347 auto t_row_base = row_data.getFTensor0N();
348 for (auto gg = 0; gg != nb_integration_pts; ++gg) {
349 double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
350 ++t_w;
351
352 auto t_mat =
353 getFTensor1FromPtr<DIM>(AssemblyDomainEleOp::locMat.data().data());
354 size_t rr = 0;
355 for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
356 const auto row_base = alpha * t_row_base;
357 auto t_col_diff_base = col_data.getFTensor1DiffN<DIM>(gg, 0);
358 for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / DIM; cc++) {
359 t_mat(i) += row_base * (t_res_diff(i, j) * t_col_diff_base(j));
360 ++t_mat;
361 ++t_col_diff_base;
362 }
363 ++t_row_base;
364 }
365 for (; rr != nb_row_base_functions; ++rr)
366 ++t_row_base;
367
368 ++t_res_diff;
369 }
370
372}
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
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_LogStrain_dUImpl< DIM, GAUSS, AssemblyDomainEleOp >::commonDataPtr
private

Definition at line 267 of file PlasticOpsLargeStrains.hpp.

◆ commonHenckyDataPtr

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

Definition at line 268 of file PlasticOpsLargeStrains.hpp.

◆ mDPtr

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

Definition at line 269 of file PlasticOpsLargeStrains.hpp.

◆ resDiff

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

Definition at line 270 of file PlasticOpsLargeStrains.hpp.


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