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

[Hencky Stress] More...

#include "tutorials/vec-2_nonlinear_elasticity/src/HenckyOps.hpp"

Inheritance diagram for HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >:
[legend]
Collaboration diagram for HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >:
[legend]

Public Member Functions

 OpCalculateHenckyThermalStressImpl (const std::string field_name, boost::shared_ptr< VectorDouble > temperature, boost::shared_ptr< CommonData > common_data, boost::shared_ptr< VectorDouble > coeff_expansion_ptr, boost::shared_ptr< double > ref_temp_ptr)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 
 OpCalculateHenckyThermalStressImpl (const std::string field_name, boost::shared_ptr< VectorDouble > temperature, boost::shared_ptr< CommonData > common_data, boost::shared_ptr< VectorDouble > coeff_expansion_ptr, boost::shared_ptr< double > ref_temp_ptr)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 

Private Attributes

boost::shared_ptr< CommonData > commonDataPtr
 
boost::shared_ptr< VectorDouble > tempPtr
 
boost::shared_ptr< VectorDouble > coeffExpansionPtr
 
boost::shared_ptr< double > refTempPtr
 

Detailed Description

template<int DIM, typename DomainEleOp, int S>
struct HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >

[Hencky Stress]

Definition at line 345 of file HenckyOps.hpp.

Constructor & Destructor Documentation

◆ OpCalculateHenckyThermalStressImpl() [1/2]

template<int DIM, typename DomainEleOp , int S>
HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >::OpCalculateHenckyThermalStressImpl ( const std::string  field_name,
boost::shared_ptr< VectorDouble >  temperature,
boost::shared_ptr< CommonData >  common_data,
boost::shared_ptr< VectorDouble >  coeff_expansion_ptr,
boost::shared_ptr< double >  ref_temp_ptr 
)
inline

Definition at line 348 of file HenckyOps.hpp.

353 : DomainEleOp(field_name, DomainEleOp::OPROW), tempPtr(temperature),
354 commonDataPtr(common_data), coeffExpansionPtr(coeff_expansion_ptr),
355 refTempPtr(ref_temp_ptr) {
356 std::fill(&DomainEleOp::doEntities[MBEDGE],
357 &DomainEleOp::doEntities[MBMAXTYPE], false);
358 }
DomainEleOp
DomainEle::UserDataOperator DomainEleOp
Finire element operator type.
Definition child_and_parent.cpp:37
field_name
constexpr auto field_name
Definition poisson_2d_homogeneous.cpp:13

◆ OpCalculateHenckyThermalStressImpl() [2/2]

template<int DIM, typename DomainEleOp , int S>
HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >::OpCalculateHenckyThermalStressImpl ( const std::string  field_name,
boost::shared_ptr< VectorDouble >  temperature,
boost::shared_ptr< CommonData >  common_data,
boost::shared_ptr< VectorDouble >  coeff_expansion_ptr,
boost::shared_ptr< double >  ref_temp_ptr 
)
inline

Definition at line 358 of file HenckyOps.hpp.

363 : DomainEleOp(field_name, DomainEleOp::OPROW), tempPtr(temperature),
364 commonDataPtr(common_data), coeffExpansionPtr(coeff_expansion_ptr),
365 refTempPtr(ref_temp_ptr) {
366 std::fill(&DomainEleOp::doEntities[MBEDGE],
367 &DomainEleOp::doEntities[MBMAXTYPE], false);
368 }

Member Function Documentation

◆ doWork() [1/2]

template<int DIM, typename DomainEleOp , int S>
MoFEMErrorCode HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >::doWork ( int  side,
EntityType  type,
EntData &  data 
)
inline

Definition at line 360 of file HenckyOps.hpp.

360 {
361 MoFEMFunctionBegin;
362
363 FTensor::Index<'i', DIM> i;
364 FTensor::Index<'j', DIM> j;
365 FTensor::Index<'k', DIM> k;
366 FTensor::Index<'l', DIM> l;
367
368 constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
369
370 // const size_t nb_gauss_pts = matGradPtr->size2();
371 const size_t nb_gauss_pts = DomainEleOp::getGaussPts().size2();
372 auto t_D = getFTensor4DdgFromMat<DIM, DIM, S>(*commonDataPtr->matDPtr);
373 auto t_logC = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matLogC);
374 auto t_logC_dC = getFTensor4DdgFromMat<DIM, DIM>(commonDataPtr->matLogCdC);
375 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
376 commonDataPtr->matHenckyStress.resize(nb_gauss_pts, size_symm, false);
377 commonDataPtr->matFirstPiolaStress.resize(nb_gauss_pts, DIM * DIM, false);
378 commonDataPtr->matSecondPiolaStress.resize(nb_gauss_pts, size_symm, false);
379 auto t_T = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matHenckyStress);
380 auto t_P = getFTensor2FromMat<DIM, DIM>(commonDataPtr->matFirstPiolaStress);
381 auto t_S =
382 getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matSecondPiolaStress);
383 auto t_grad = getFTensor2FromMat<DIM, DIM>(*(commonDataPtr->matGradPtr));
384 auto t_temp = getFTensor0FromVec(*tempPtr);
385
386 auto t_coeff_exp = FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM>();
387 t_coeff_exp(i, j) = 0;
388 for (auto d = 0; d != SPACE_DIM; ++d) {
389 t_coeff_exp(d, d) = (*coeffExpansionPtr)[d];
390 }
391
392 for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
393#ifdef HENCKY_SMALL_STRAIN
394 t_P(i, j) = t_D(i, j, k, l) *
395 (t_grad(k, l) - t_coeff_exp(k, l) * (t_temp - (*refTempPtr)));
396#else
397 t_T(i, j) = t_D(i, j, k, l) *
398 (t_logC(k, l) - t_coeff_exp(k, l) * (t_temp - (*refTempPtr)));
399 FTensor::Tensor2<double, DIM, DIM> t_F;
400 t_F(i, j) = t_grad(i, j) + t_kd(i, j);
401 t_S(k, l) = t_T(i, j) * t_logC_dC(i, j, k, l);
402 t_P(i, l) = t_F(i, k) * t_S(k, l);
403#endif
404 ++t_grad;
405 ++t_logC;
406 ++t_logC_dC;
407 ++t_P;
408 ++t_T;
409 ++t_S;
410 ++t_D;
411 ++t_temp;
412 }
413
414 MoFEMFunctionReturn(0);
415 }
SPACE_DIM
constexpr int SPACE_DIM
Definition child_and_parent.cpp:17
FTensor::Kronecker_Delta
Kronecker Delta class.
Definition Kronecker_Delta.hpp:15
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition definitions.h:359
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition definitions.h:432
t_kd
constexpr auto t_kd
Definition free_surface.cpp:161
i
FTensor::Index< 'i', SPACE_DIM > i
Definition hcurl_divergence_operator_2d.cpp:27
l
FTensor::Index< 'l', 3 > l
Definition matrix_function.cpp:23
j
FTensor::Index< 'j', 3 > j
Definition matrix_function.cpp:21
k
FTensor::Index< 'k', 3 > k
Definition matrix_function.cpp:22
sdf_hertz_2d_axisymm_plane.d
float d
Definition sdf_hertz_2d_axisymm_plane.py:4
size_symm
constexpr auto size_symm
Definition plastic.cpp:42

◆ doWork() [2/2]

template<int DIM, typename DomainEleOp , int S>
MoFEMErrorCode HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >::doWork ( int  side,
EntityType  type,
EntData &  data 
)
inline

Definition at line 370 of file HenckyOps.hpp.

370 {
371 MoFEMFunctionBegin;
372
373 FTensor::Index<'i', DIM> i;
374 FTensor::Index<'j', DIM> j;
375 FTensor::Index<'k', DIM> k;
376 FTensor::Index<'l', DIM> l;
377
378 constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
379
380 // const size_t nb_gauss_pts = matGradPtr->size2();
381 const size_t nb_gauss_pts = DomainEleOp::getGaussPts().size2();
382 auto t_D =
383 getFTensor4DdgFromMat<DIM, DIM, S,
384 DataLayoutTraits<DataLayout::GaussByCoeffs>>(
385 *commonDataPtr->matDPtr);
386 auto t_logC = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matLogC);
387 auto t_logC_dC = getFTensor4DdgFromMat<DIM, DIM>(commonDataPtr->matLogCdC);
388 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
389 commonDataPtr->matHenckyStress.resize(size_symm, nb_gauss_pts, false);
390 commonDataPtr->matFirstPiolaStress.resize(DIM * DIM, nb_gauss_pts, false);
391 commonDataPtr->matSecondPiolaStress.resize(size_symm, nb_gauss_pts, false);
392 auto t_T = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matHenckyStress);
393 auto t_P = getFTensor2FromMat<DIM, DIM>(commonDataPtr->matFirstPiolaStress);
394 auto t_S =
395 getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matSecondPiolaStress);
396 auto t_grad = getFTensor2FromMat<DIM, DIM>(*(commonDataPtr->matGradPtr));
397 auto t_temp = getFTensor0FromVec(*tempPtr);
398
399 auto t_coeff_exp = FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM>();
400 t_coeff_exp(i, j) = 0;
401 for (auto d = 0; d != SPACE_DIM; ++d) {
402 t_coeff_exp(d, d) = (*coeffExpansionPtr)[d];
403 }
404
405 for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
406#ifdef HENCKY_SMALL_STRAIN
407 t_P(i, j) = t_D(i, j, k, l) *
408 (t_grad(k, l) - t_coeff_exp(k, l) * (t_temp - (*refTempPtr)));
409#else
410 t_T(i, j) = t_D(i, j, k, l) *
411 (t_logC(k, l) - t_coeff_exp(k, l) * (t_temp - (*refTempPtr)));
412 FTensor::Tensor2<double, DIM, DIM> t_F;
413 t_F(i, j) = t_grad(i, j) + t_kd(i, j);
414 t_S(k, l) = t_T(i, j) * t_logC_dC(i, j, k, l);
415 t_P(i, l) = t_F(i, k) * t_S(k, l);
416#endif
417 ++t_grad;
418 ++t_logC;
419 ++t_logC_dC;
420 ++t_P;
421 ++t_T;
422 ++t_S;
423 ++t_D;
424 ++t_temp;
425 }
426
427 MoFEMFunctionReturn(0);
428 }

Member Data Documentation

◆ coeffExpansionPtr

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

Definition at line 420 of file HenckyOps.hpp.

◆ commonDataPtr

template<int DIM, typename DomainEleOp , int S>
boost::shared_ptr< CommonData > HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >::commonDataPtr
private

Definition at line 418 of file HenckyOps.hpp.

◆ refTempPtr

template<int DIM, typename DomainEleOp , int S>
boost::shared_ptr< double > HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >::refTempPtr
private

Definition at line 421 of file HenckyOps.hpp.

◆ tempPtr

template<int DIM, typename DomainEleOp , int S>
boost::shared_ptr< VectorDouble > HenckyOps::OpCalculateHenckyThermalStressImpl< DIM, GAUSS, DomainEleOp, S >::tempPtr
private

Definition at line 419 of file HenckyOps.hpp.

The documentation for this struct was generated from the following files:
Generated by Doxygen 1.9.8 and hosted at University of Glasgow