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

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

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

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

Public Member Functions
	OpCalculateAdiabaticHeatingLhsdUImpl (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::OpCalculateAdiabaticHeatingLhsdUImpl< DIM, GAUSS, AssemblyDomainEleOp >

Definition at line 497 of file ThermoPlasticOps.hpp.

Constructor & Destructor Documentation

OpCalculateAdiabaticHeatingLhsdUImpl()

template<int DIM, typename AssemblyDomainEleOp >

ThermoPlasticOps::OpCalculateAdiabaticHeatingLhsdUImpl< DIM, GAUSS, AssemblyDomainEleOp >::OpCalculateAdiabaticHeatingLhsdUImpl ( 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 499 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::OpCalculateAdiabaticHeatingLhsdUImpl< DIM, GAUSS, AssemblyDomainEleOp >::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

protected

Definition at line 525 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;
  FTensor::Index<'o', DIM> o;
  FTensor::Index<'p', DIM> p;
 
  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_row_base = row_data.getFTensor0N();
 
  auto t_logC_dC =
      getFTensor4DdgFromMat<DIM, DIM>(elasticCommonDataPtr->matLogCdC);
  auto t_C_dF =
      getFTensor4FromMat<DIM, DIM, DIM, DIM, 1>(elasticCommonDataPtr->matCdF);
 
  // 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();
 
  for (auto gg = 0; gg != nb_integration_pts; ++gg) {
    auto t_vec = getFTensor1FromPtr<DIM>(&this->locMat(0, 0));
    double alpha = this->getMeasure() * t_w;
 
    auto rr = 0;
    for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
      auto t_col_grad_base = col_data.getFTensor1DiffN<DIM>(gg, 0);
      for (auto cc = 0; cc != AssemblyDomainEleOp::nbCols / DIM; cc++) {
#ifdef HENCKY_SMALL_STRAIN
        // TODO: implement correct small strain problem
        t_vec(m) +=
            0.5 * 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_logC_dC(k, l, o, p) * t_C_dF(o, p, m, n)) *
            t_col_grad_base(n);
#else
        // Multiply by a half for dEdC from dlogCdC
        t_vec(m) +=
            0.5 * 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_logC_dC(k, l, o, p) * t_C_dF(o, p, m, n)) *
            t_col_grad_base(n);
#endif
        ++t_col_grad_base;
        ++t_vec;
      }
 
      ++t_row_base;
    }
    for (; rr != nb_row_base_functions; ++rr)
      ++t_row_base;
 
    ++t_w;
    ++t_plastic_strain_rate;
    ++t_D;
    ++t_coords;
    ++t_logC_dC;
    ++t_C_dF;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

elasticCommonDataPtr

template<int DIM, typename AssemblyDomainEleOp >

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

protected

Definition at line 515 of file ThermoPlasticOps.hpp.

entsPtr

template<int DIM, typename AssemblyDomainEleOp >

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

protected

Definition at line 518 of file ThermoPlasticOps.hpp.

inelasticHeatingFunction

template<int DIM, typename AssemblyDomainEleOp >

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

protected

Definition at line 517 of file ThermoPlasticOps.hpp.

plasticCommonDataPtr

template<int DIM, typename AssemblyDomainEleOp >

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

protected

Definition at line 516 of file ThermoPlasticOps.hpp.

---

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

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