CellEngineering::OpElasticEnergy Struct Reference

Inheritance diagram for CellEngineering::OpElasticEnergy:

Collaboration diagram for CellEngineering::OpElasticEnergy:

Public Member Functions
	OpElasticEnergy (boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< MatrixDouble > stress_ptr, double &energy)
MoFEMErrorCode	doWork (int, EntityType, EntitiesFieldData::EntData &) override

Private Attributes
boost::shared_ptr< MatrixDouble >	strainPtr
boost::shared_ptr< MatrixDouble >	stressPtr
double &	energy

Detailed Description

Examples

cell_forces.cpp.

Definition at line 429 of file cell_forces.cpp.

Constructor & Destructor Documentation

OpElasticEnergy()

CellEngineering::OpElasticEnergy::OpElasticEnergy ( boost::shared_ptr< MatrixDouble >  strain_ptr, boost::shared_ptr< MatrixDouble >  stress_ptr, double &  energy  )

inline

Definition at line 432 of file cell_forces.cpp.

      : VolumeElementForcesAndSourcesCore::UserDataOperator(
            NOSPACE, UserDataOperator::OPSPACE),
        strainPtr(strain_ptr), stressPtr(stress_ptr), energy(energy) {}

Member Function Documentation

doWork()

MoFEMErrorCode CellEngineering::OpElasticEnergy::doWork ( int  , EntityType  , EntitiesFieldData::EntData &    )

inlineoverride

Examples

cell_forces.cpp.

Definition at line 438 of file cell_forces.cpp.

                                                             {
    MoFEMFunctionBegin;
 
    const size_t nb_gauss_pts = getGaussPts().size2();
    using DL = DataLayoutTraits<DataLayout::GaussByCoeffs>;
    auto get_strain =
        MatrixSizeHelper<GetFTensor2SymmetricFromMatType<SPACE_DIM, -1, DL>,
                         DL>::get(*strainPtr, nb_gauss_pts);
    auto get_stress =
        MatrixSizeHelper<GetFTensor2SymmetricFromMatType<SPACE_DIM, -1, DL>,
                         DL>::get(*stressPtr, nb_gauss_pts);
    auto t_strain = get_strain();
    auto t_stress = get_stress();
    FTensor::Index<'i', SPACE_DIM> i;
    FTensor::Index<'j', SPACE_DIM> j;
 
    for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
      const double alpha = getMeasure() * getGaussPts()(SPACE_DIM, gg);
      energy += 0.5 * alpha * t_strain(i, j) * t_stress(i, j);
      ++t_strain;
      ++t_stress;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

energy

double& CellEngineering::OpElasticEnergy::energy

private

Examples

cell_forces.cpp.

Definition at line 468 of file cell_forces.cpp.

strainPtr

boost::shared_ptr<MatrixDouble> CellEngineering::OpElasticEnergy::strainPtr

private

Examples

cell_forces.cpp.

Definition at line 466 of file cell_forces.cpp.

stressPtr

boost::shared_ptr<MatrixDouble> CellEngineering::OpElasticEnergy::stressPtr

private

Examples

cell_forces.cpp.

Definition at line 467 of file cell_forces.cpp.

---

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

• users_modules/cell_engineering/cell_forces.cpp