#include <tutorials/adv-2/src/ThermoElasticOps.hpp>

Inheritance diagram for ThermoElasticOps::OpStressThermal:

Collaboration diagram for ThermoElasticOps::OpStressThermal:

Public Member Functions
DEPRECATED	OpStressThermal (const std::string field_name, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< VectorDouble > temp_ptr, boost::shared_ptr< MatrixDouble > m_D_ptr, boost::shared_ptr< VectorDouble > coeff_expansion_ptr, boost::shared_ptr< MatrixDouble > stress_ptr)

	OpStressThermal (boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< VectorDouble > temp_ptr, boost::shared_ptr< MatrixDouble > m_D_ptr, boost::shared_ptr< VectorDouble > coeff_expansion_ptr, boost::shared_ptr< MatrixDouble > stress_ptr)

MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
	[Calculate stress] More...

Private Attributes
boost::shared_ptr< MatrixDouble >	strainPtr

boost::shared_ptr< VectorDouble >	tempPtr

boost::shared_ptr< MatrixDouble >	mDPtr

boost::shared_ptr< VectorDouble >	coeffExpansionPtr

boost::shared_ptr< MatrixDouble >	stressPtr

Detailed Description

Examples: thermo_elastic.cpp.

Definition at line 23 of file ThermoElasticOps.hpp.

Constructor & Destructor Documentation

◆ OpStressThermal() [1/2]

ThermoElasticOps::OpStressThermal::OpStressThermal	(	const std::string	field_name,
		boost::shared_ptr< MatrixDouble >	strain_ptr,
		boost::shared_ptr< VectorDouble >	temp_ptr,
		boost::shared_ptr< MatrixDouble >	m_D_ptr,
		boost::shared_ptr< VectorDouble >	coeff_expansion_ptr,
		boost::shared_ptr< MatrixDouble >	stress_ptr
	)

Deprecated:: do not use this constructor

Examples: ThermoElasticOps.hpp.

Definition at line 115 of file ThermoElasticOps.hpp.

     : DomainEleOp(field_name, DomainEleOp::OPROW), strainPtr(strain_ptr),
       tempPtr(temp_ptr), mDPtr(m_D_ptr), coeffExpansionPtr(coeff_expansion_ptr),
       stressPtr(stress_ptr) {
   // Operator is only executed for vertices
   std::fill(&doEntities[MBEDGE], &doEntities[MBMAXTYPE], false);
 }

◆ OpStressThermal() [2/2]

ThermoElasticOps::OpStressThermal::OpStressThermal	(	boost::shared_ptr< MatrixDouble >	strain_ptr,
		boost::shared_ptr< VectorDouble >	temp_ptr,
		boost::shared_ptr< MatrixDouble >	m_D_ptr,
		boost::shared_ptr< VectorDouble >	coeff_expansion_ptr,
		boost::shared_ptr< MatrixDouble >	stress_ptr
	)

Definition at line 128 of file ThermoElasticOps.hpp.

     : DomainEleOp(NOSPACE, DomainEleOp::OPSPACE), strainPtr(strain_ptr),
       tempPtr(temp_ptr), mDPtr(m_D_ptr), coeffExpansionPtr(coeff_expansion_ptr),
       stressPtr(stress_ptr) {}

Member Function Documentation

◆ doWork()

MoFEMErrorCode ThermoElasticOps::OpStressThermal::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

[Calculate stress]

Examples: ThermoElasticOps.hpp.

Definition at line 139 of file ThermoElasticOps.hpp.

                                                       {
   MoFEMFunctionBegin;
   const auto nb_gauss_pts = getGaussPts().size2();
   stressPtr->resize((SPACE_DIM * (SPACE_DIM + 1)) / 2, nb_gauss_pts);
  
   auto t_D = getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*mDPtr);
   auto t_strain = getFTensor2SymmetricFromMat<SPACE_DIM>(*strainPtr);
   auto t_stress = getFTensor2SymmetricFromMat<SPACE_DIM>(*stressPtr);
   auto t_temp = getFTensor0FromVec(*tempPtr);
  
   FTensor::Index<'i', SPACE_DIM> i;
   FTensor::Index<'j', SPACE_DIM> j;
   FTensor::Index<'k', SPACE_DIM> k;
   FTensor::Index<'l', SPACE_DIM> l;
  
   auto t_coeff_exp = FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM>();
   t_coeff_exp(i, j) = 0;
   for (auto d = 0; d != SPACE_DIM; ++d) {
     t_coeff_exp(d, d) = (*coeffExpansionPtr)[d];
   }
  
   for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
     t_stress(i, j) = t_D(i, j, k, l) *
                      (t_strain(k, l) - t_coeff_exp(k, l) * (t_temp - ref_temp));
     ++t_strain;
     ++t_stress;
     ++t_temp;
   }
   MoFEMFunctionReturn(0);
 }