Inheritance diagram for OpCalculateLameStress< DIM >:

Collaboration diagram for OpCalculateLameStress< DIM >:

Public Member Functions
	OpCalculateLameStress (double m_u, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< VectorDouble > pressure_ptr)

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

Private Attributes
double	mU

boost::shared_ptr< MatrixDouble >	stressPtr

boost::shared_ptr< MatrixDouble >	strainPtr

boost::shared_ptr< VectorDouble >	pressurePtr

Detailed Description

template<int DIM>
struct OpCalculateLameStress< DIM >

Definition at line 208 of file incompressible_elasticity.cpp.

Constructor & Destructor Documentation

◆ OpCalculateLameStress()

template<int DIM>

OpCalculateLameStress< DIM >::OpCalculateLameStress	(	double	m_u,
		boost::shared_ptr< MatrixDouble >	stress_ptr,
		boost::shared_ptr< MatrixDouble >	strain_ptr,
		boost::shared_ptr< VectorDouble >	pressure_ptr )

inline

Definition at line 209 of file incompressible_elasticity.cpp.

      : ForcesAndSourcesCore::UserDataOperator(NOSPACE, OPLAST), mU(m_u),
        stressPtr(stress_ptr), strainPtr(strain_ptr),
        pressurePtr(pressure_ptr) {}

Member Function Documentation

◆ doWork()

template<int DIM>

MoFEMErrorCode OpCalculateLameStress< DIM >::doWork	(	int	side,
		EntityType	type,
		EntitiesFieldData::EntData &	data )

inline

Definition at line 216 of file incompressible_elasticity.cpp.

                                                        {
    MoFEMFunctionBegin;
    // Define Indicies
    FTensor::Index<'i', SPACE_DIM> i;
    FTensor::Index<'j', SPACE_DIM> j;
 
    // Define Kronecker Delta
    constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<double>();
 
    // Number of Gauss points
    const size_t nb_gauss_pts = getGaussPts().size2();
 
    stressPtr->resize((DIM * (DIM + 1)) / 2, nb_gauss_pts);
    auto t_stress = getFTensor2SymmetricFromMat<DIM>(*(stressPtr));
    auto t_strain = getFTensor2SymmetricFromMat<DIM>(*(strainPtr));
    auto t_pressure = getFTensor0FromVec(*(pressurePtr));
 
    const double l_mu = mU;
    // Extract matrix from data matrix
    for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
 
      t_stress(i, j) = t_pressure * t_kd(i, j) + 2. * l_mu * t_strain(i, j);
 
      ++t_strain;
      ++t_stress;
      ++t_pressure;
    }
 
    MoFEMFunctionReturn(0);
  }