HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S > Struct Template Reference

[Op Hencky Tangent impl] More...

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

Inheritance diagram for HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >:

Collaboration diagram for HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >:

Public Member Functions
	OpCalculateHenckyThermalStressdTImpl (const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< CommonData > elastic_common_data_ptr, boost::shared_ptr< VectorDouble > coeff_expansion_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Private Attributes
boost::shared_ptr< CommonData >	elasticCommonDataPtr
boost::shared_ptr< VectorDouble >	coeffExpansionPtr

Detailed Description

template<int DIM, typename AssemblyDomainEleOp, int S>
struct HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >

[Op Hencky Tangent impl]

Definition at line 595 of file HenckyOps.hpp.

Constructor & Destructor Documentation

OpCalculateHenckyThermalStressdTImpl()

template<int DIM, typename AssemblyDomainEleOp , int S>

HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::OpCalculateHenckyThermalStressdTImpl	(	const std::string	row_field_name,
		const std::string	col_field_name,
		boost::shared_ptr< CommonData >	elastic_common_data_ptr,
		boost::shared_ptr< VectorDouble >	coeff_expansion_ptr
	)

Definition at line 611 of file HenckyOps.hpp.

    : AssemblyDomainEleOp(row_field_name, col_field_name,
                          AssemblyDomainEleOp::OPROWCOL),
      elasticCommonDataPtr(elastic_common_data_ptr),
      coeffExpansionPtr(coeff_expansion_ptr) {
  this->sYmm = false;
}

Member Function Documentation

iNtegrate()

template<int DIM, typename AssemblyDomainEleOp , int S>

MoFEMErrorCode HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::iNtegrate	(	EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

Definition at line 625 of file HenckyOps.hpp.

                                                  {
  MoFEMFunctionBegin;
 
  auto &locMat = AssemblyDomainEleOp::locMat;
 
  const auto nb_integration_pts = row_data.getN().size1();
  const auto nb_row_base_functions = row_data.getN().size2();
  auto t_w = this->getFTensor0IntegrationWeight();
 
  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
  auto t_row_diff_base = row_data.getFTensor1DiffN<DIM>();
  auto t_D = getFTensor4DdgFromMat<DIM, DIM, 0>(*elasticCommonDataPtr->matDPtr);
  auto t_grad =
      getFTensor2FromMat<DIM, DIM>(*(elasticCommonDataPtr->matGradPtr));
  auto t_logC_dC =
      getFTensor4DdgFromMat<DIM, DIM>(elasticCommonDataPtr->matLogCdC);
  FTensor::Tensor2_symmetric<double, DIM> t_eigen_strain;
  FTensor::Tensor2<double, DIM, DIM> t_F;
 
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
  FTensor::Index<'m', DIM> m;
  FTensor::Index<'n', DIM> n;
  FTensor::Index<'o', DIM> o;
 
  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];
  }
 
  t_eigen_strain(i, j) = (t_D(i, j, k, l) * t_coeff_exp(k, l));
 
  for (auto gg = 0; gg != nb_integration_pts; ++gg) {
 
    t_F(i, j) = t_grad(i, j) + t_kd(i, j);
 
    double alpha = this->getMeasure() * t_w;
    auto rr = 0;
    for (; rr != AssemblyDomainEleOp::nbRows / DIM; ++rr) {
      auto t_mat =
          getFTensor1FromMat<DIM, 1,
                             DataLayoutTraits<DataLayout::CoeffsByGauss>>(
              locMat, rr * DIM);
      auto t_col_base = col_data.getFTensor0N(gg, 0);
      for (auto cc = 0; cc != AssemblyDomainEleOp::nbCols; cc++) {
#ifdef HENCKY_SMALL_STRAIN
        t_mat(i) -=
            (t_row_diff_base(j) * t_eigen_strain(i, j)) * (t_col_base * alpha);
#else
        t_mat(i) -= (t_row_diff_base(j) *
                     (t_F(i, o) * ((t_D(m, n, k, l) * t_coeff_exp(k, l)) *
                                   t_logC_dC(m, n, o, j)))) *
                    (t_col_base * alpha);
#endif
 
        ++t_mat;
        ++t_col_base;
      }
 
      ++t_row_diff_base;
    }
    for (; rr != nb_row_base_functions; ++rr)
      ++t_row_diff_base;
 
    ++t_w;
    ++t_grad;
    ++t_logC_dC;
    ++t_D;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

coeffExpansionPtr

template<int DIM, typename AssemblyDomainEleOp , int S>

boost::shared_ptr<VectorDouble> HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::coeffExpansionPtr

private

Definition at line 607 of file HenckyOps.hpp.

elasticCommonDataPtr

template<int DIM, typename AssemblyDomainEleOp , int S>

boost::shared_ptr<CommonData> HenckyOps::OpCalculateHenckyThermalStressdTImpl< DIM, GAUSS, AssemblyDomainEleOp, S >::elasticCommonDataPtr

private

Definition at line 606 of file HenckyOps.hpp.

---

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

• tutorials/vec-2_nonlinear_elasticity/src/HenckyOps.hpp