#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< double > ref_temp_ptr, boost::shared_ptr< MatrixDouble > stress_ptr)

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

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< double >	refTempPtr

boost::shared_ptr< MatrixDouble >	stressPtr

Detailed Description

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: mofem/tutorials/adv-2/src/ThermoElasticOps.hpp.

Definition at line 117 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< double >	ref_temp_ptr,
		boost::shared_ptr< MatrixDouble >	stress_ptr
	)

Definition at line 130 of file ThermoElasticOps.hpp.

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

Member Function Documentation

◆ doWork()

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

[Calculate stress]

Examples: mofem/tutorials/adv-2/src/ThermoElasticOps.hpp.

Definition at line 142 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 - (*refTempPtr)));
    ++t_strain;
    ++t_stress;
    ++t_temp;
  }
  MoFEMFunctionReturn(0);
}