ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp > Struct Template Reference

#include "users_modules/multifield-thermoplasticity-private/src/ThermoPlasticOps.hpp"

Inheritance diagram for ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp >:

Collaboration diagram for ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp >:

Public Member Functions
	OpCalculateAdiabaticHeatingLhsdEPImpl (const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< HenckyOps::CommonData > elastic_common_data_ptr, boost::shared_ptr< PlasticOps::CommonData > plastic_common_data_ptr, ScalarFun inelastic_heating_function, boost::shared_ptr< Range > ents_ptr=nullptr)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Protected Attributes
boost::shared_ptr< HenckyOps::CommonData >	elasticCommonDataPtr
boost::shared_ptr< PlasticOps::CommonData >	plasticCommonDataPtr
ScalarFun	inelasticHeatingFunction
boost::shared_ptr< Range >	entsPtr

Detailed Description

template<int DIM, typename AssemblyDomainEleOp>
struct ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp >

Definition at line 348 of file ThermoPlasticOps.hpp.

Constructor & Destructor Documentation

OpCalculateAdiabaticHeatingLhsdEPImpl()

template<int DIM, typename AssemblyDomainEleOp >

ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::OpCalculateAdiabaticHeatingLhsdEPImpl ( const std::string  row_field_name, const std::string  col_field_name, boost::shared_ptr< HenckyOps::CommonData >  elastic_common_data_ptr, boost::shared_ptr< PlasticOps::CommonData >  plastic_common_data_ptr, ScalarFun  inelastic_heating_function, boost::shared_ptr< Range >  ents_ptr = nullptr  )

inline

Definition at line 350 of file ThermoPlasticOps.hpp.

      : AssemblyDomainEleOp(row_field_name, col_field_name,
                            AssemblyDomainEleOp::OPROWCOL),
        elasticCommonDataPtr(elastic_common_data_ptr),
        plasticCommonDataPtr(plastic_common_data_ptr),
        inelasticHeatingFunction(inelastic_heating_function),
        entsPtr(ents_ptr) {
    this->sYmm = false;
  }

Member Function Documentation

iNtegrate()

template<int DIM, typename AssemblyDomainEleOp >

MoFEMErrorCode ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

protected

Definition at line 376 of file ThermoPlasticOps.hpp.

                                                  {
  MoFEMFunctionBegin;
 
  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;
  constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
  FTensor::Index<'L', size_symm> L;
 
  const auto nb_integration_pts = row_data.getN().size1();
  const auto nb_row_base_functions = row_data.getN().size2();
  auto t_w = this->getFTensor0IntegrationWeight();
 
  auto t_L = PlasticOps::symm_L_tensor(FTensor::Number<DIM>());
 
  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
  auto t_row_base = row_data.getFTensor0N();
 
  auto t_stress = getFTensor2SymmetricFromMat<DIM>(
      *(elasticCommonDataPtr->getMatHenckyStress()));
 
  // get plastic strain rate values
  auto t_plastic_strain_rate = getFTensor2SymmetricFromMat<DIM>(
      *(plasticCommonDataPtr->getPlasticStrainDotPtr()));
 
  // get elasticity tensor
  auto t_D =
      getFTensor4DdgFromMat<DIM, DIM, 0>(*(elasticCommonDataPtr->matDPtr));
 
  // get coordinate at integration points
  auto t_coords = AssemblyDomainEleOp::getFTensor1CoordsAtGaussPts();
 
  auto t_diff = PlasticOps::diff_tensor(FTensor::Number<DIM>());
 
  for (auto gg = 0; gg != nb_integration_pts; ++gg) {
    // get the local vector to store the symmetric tensor of size_symm
    auto t_vec = getFTensor1FromPtr<size_symm>(&this->locMat(0, 0));
    double alpha = this->getMeasure() * t_w;
    auto rr = 0;
    for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
      auto t_col_base = col_data.getFTensor0N(gg, 0);
      for (auto cc = 0; cc != AssemblyDomainEleOp::nbCols / size_symm; cc++) {
#ifdef HENCKY_SMALL_STRAIN
        t_vec(L) -=
            alpha * t_row_base *
            inelasticHeatingFunction(t_coords(0), t_coords(1), t_coords(2)) *
            (t_plastic_strain_rate(i, j) * t_D(i, j, k, l) *
                 t_diff(k, l, m, n) -
             this->getTSa() * t_stress(i, j) * t_diff(i, j, m, n)) *
            t_col_base * t_L(m, n, L);
#else
        t_vec(L) -=
            alpha * t_row_base *
            inelasticHeatingFunction(t_coords(0), t_coords(1), t_coords(2)) *
            (t_plastic_strain_rate(i, j) * t_D(i, j, k, l) *
                 t_diff(k, l, m, n) -
             this->getTSa() * t_stress(i, j) * t_diff(i, j, m, n)) *
            t_col_base * t_L(m, n, L);
#endif
        ++t_vec;
        ++t_col_base;
      }
 
      ++t_row_base;
    }
    for (; rr != nb_row_base_functions; ++rr)
      ++t_row_base;
 
    ++t_w;
    ++t_plastic_strain_rate;
    ++t_D;
    ++t_coords;
    ++t_stress;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

elasticCommonDataPtr

template<int DIM, typename AssemblyDomainEleOp >

boost::shared_ptr<HenckyOps::CommonData> ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::elasticCommonDataPtr

protected

Definition at line 366 of file ThermoPlasticOps.hpp.

entsPtr

template<int DIM, typename AssemblyDomainEleOp >

boost::shared_ptr<Range> ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::entsPtr

protected

Definition at line 369 of file ThermoPlasticOps.hpp.

inelasticHeatingFunction

template<int DIM, typename AssemblyDomainEleOp >

ScalarFun ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::inelasticHeatingFunction

protected

Definition at line 368 of file ThermoPlasticOps.hpp.

plasticCommonDataPtr

template<int DIM, typename AssemblyDomainEleOp >

boost::shared_ptr<PlasticOps::CommonData> ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::plasticCommonDataPtr

protected

Definition at line 367 of file ThermoPlasticOps.hpp.

---

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

• users_modules/multifield-thermoplasticity-private/src/ThermoPlasticOps.hpp