Inheritance diagram for OpCalculatePiolaIncompressibleNH< DIM_0, DIM_1 >:

Collaboration diagram for OpCalculatePiolaIncompressibleNH< DIM_0, DIM_1 >:

Public Member Functions
	OpCalculatePiolaIncompressibleNH (double shear_modulus, double bulk_modulus, double m_u, double lambda_lamme, boost::shared_ptr< MatrixDouble > first_piola_ptr, boost::shared_ptr< MatrixDouble > def_grad_ptr, boost::shared_ptr< MatrixDouble > inv_def_grad_ptr, boost::shared_ptr< VectorDouble > det)

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

Private Attributes
double	shearModulus

double	bulkModulus

double	mU

double	lammeLambda

boost::shared_ptr< MatrixDouble >	firstPiolaPtr

boost::shared_ptr< MatrixDouble >	defGradPtr

boost::shared_ptr< MatrixDouble >	invDefGradPtr

boost::shared_ptr< VectorDouble >	dEt

Detailed Description

template<int DIM_0, int DIM_1>
struct OpCalculatePiolaIncompressibleNH< DIM_0, DIM_1 >

Examples: mofem/tutorials/adv-4_dynamic_first_order_con_law/dynamic_first_order_con_law.cpp.

Definition at line 267 of file dynamic_first_order_con_law.cpp.

Constructor & Destructor Documentation

◆ OpCalculatePiolaIncompressibleNH()

template<int DIM_0, int DIM_1>

OpCalculatePiolaIncompressibleNH< DIM_0, DIM_1 >::OpCalculatePiolaIncompressibleNH	(	double	shear_modulus,
		double	bulk_modulus,
		double	m_u,
		double	lambda_lamme,
		boost::shared_ptr< MatrixDouble >	first_piola_ptr,
		boost::shared_ptr< MatrixDouble >	def_grad_ptr,
		boost::shared_ptr< MatrixDouble >	inv_def_grad_ptr,
		boost::shared_ptr< VectorDouble >	det
	)

inline

Definition at line 269 of file dynamic_first_order_con_law.cpp.

      : ForcesAndSourcesCore::UserDataOperator(NOSPACE, OPLAST),
        shearModulus(shear_modulus), bulkModulus(bulk_modulus), mU(m_u),
        lammeLambda(lambda_lamme), firstPiolaPtr(first_piola_ptr),
        defGradPtr(def_grad_ptr), invDefGradPtr(inv_def_grad_ptr), dEt(det) {}

Member Function Documentation

◆ doWork()

template<int DIM_0, int DIM_1>

MoFEMErrorCode OpCalculatePiolaIncompressibleNH< DIM_0, DIM_1 >::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

inline

Examples: mofem/tutorials/adv-4_dynamic_first_order_con_law/dynamic_first_order_con_law.cpp.

Definition at line 280 of file dynamic_first_order_con_law.cpp.

                                                               {
    MoFEMFunctionBegin;
    // Define Indicies
    FTensor::Index<'i', SPACE_DIM> i;
    FTensor::Index<'j', SPACE_DIM> j;
    FTensor::Index<'k', SPACE_DIM> k;
    FTensor::Index<'l', SPACE_DIM> l;
 
    // Define Kronecker Delta
    constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
 
    // Number of Gauss points
    const size_t nb_gauss_pts = getGaussPts().size2();
 
    // Resize Piola
    firstPiolaPtr->resize(nb_gauss_pts, DIM_0 * DIM_1, false); // ignatios check
    firstPiolaPtr->clear();
 
    // Extract matrix from data matrix
    auto t_P = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*firstPiolaPtr);
    auto t_F = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*defGradPtr);
    auto t_inv_F = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*invDefGradPtr);
    auto t_det = getFTensor0FromVec<1>(*dEt);
    const double two_o_three = 2. / 3.;
    const double one_o_three = 1. / 3.;
    const double bulk_mod = bulkModulus;
    const double shear_mod = shearModulus;
    for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
 
      // Nearly incompressible NH
      // volumetric part
      t_P(i, j) = bulk_mod * (t_det - 1.) * t_det * t_inv_F(j, i);
      // deviatoric part
      t_P(i, j) +=
          shear_mod * pow(t_det, two_o_three) *
          (t_F(i, j) - one_o_three * (t_F(l, k) * t_F(l, k)) * t_inv_F(j, i));
 
      ++t_F;
      ++t_P;
      ++t_inv_F;
      ++t_det;
    }
 
    MoFEMFunctionReturn(0);
  }