Hooke< TYPE > Struct Template Reference

Hook equation. More...

#include "users_modules/basic_finite_elements/nonlinear_elastic_materials/src/Hooke.hpp"

Inheritance diagram for Hooke< TYPE >:

Collaboration diagram for Hooke< TYPE >:

Public Member Functions
	Hooke ()
virtual MoFEMErrorCode	calculateP_PiolaKirchhoffI (const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Hooke equation.
virtual MoFEMErrorCode	calculateElasticEnergy (const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	calculate density of strain energy
Public Member Functions inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
	FunctionsToCalculatePiolaKirchhoffI ()
virtual	~FunctionsToCalculatePiolaKirchhoffI ()=default
MoFEMErrorCode	calculateC_CauchyDeformationTensor ()
MoFEMErrorCode	calculateE_GreenStrain ()
MoFEMErrorCode	calculateS_PiolaKirchhoffII ()
virtual MoFEMErrorCode	calculateCauchyStress (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Function overload to implement user material.
virtual MoFEMErrorCode	setUserActiveVariables (int &nb_active_variables)
	add additional active variables
virtual MoFEMErrorCode	setUserActiveVariables (VectorDouble &activeVariables)
	Add additional independent variables More complex physical models depend on gradient of defamation and some additional variables. For example can depend on temperature. This function adds additional independent variables to the model.
virtual MoFEMErrorCode	calculatesIGma_EshelbyStress (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Calculate Eshelby stress.
virtual MoFEMErrorCode	getDataOnPostProcessor (std::map< std::string, std::vector< VectorDouble > > &field_map, std::map< std::string, std::vector< MatrixDouble > > &grad_map)
	Do operations when pre-process.

Public Attributes
ublas::matrix< TYPE >	ePs
TYPE	tR
Public Attributes inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
FTensor::Index< 'i', 3 >	i
FTensor::Index< 'j', 3 >	j
FTensor::Index< 'k', 3 >	k
double	lambda
double	mu
MatrixBoundedArray< TYPE, 9 >	F
MatrixBoundedArray< TYPE, 9 >	C
MatrixBoundedArray< TYPE, 9 >	E
MatrixBoundedArray< TYPE, 9 >	S
MatrixBoundedArray< TYPE, 9 >	invF
MatrixBoundedArray< TYPE, 9 >	P
MatrixBoundedArray< TYPE, 9 >	sIGma
MatrixBoundedArray< TYPE, 9 >	h
MatrixBoundedArray< TYPE, 9 >	H
MatrixBoundedArray< TYPE, 9 >	invH
MatrixBoundedArray< TYPE, 9 >	sigmaCauchy
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_F
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_C
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_E
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_S
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_invF
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_P
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_sIGma
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_h
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_H
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_invH
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_sigmaCauchy
TYPE	J
TYPE	eNergy
TYPE	detH
TYPE	detF
int	gG
	Gauss point number.
CommonData *	commonDataPtr
MoFEM::VolumeElementForcesAndSourcesCore::UserDataOperator *	opPtr
	pointer to finite element tetrahedral operator

Additional Inherited Members
Static Protected Member Functions inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
static auto	resizeAndSet (MatrixBoundedArray< TYPE, 9 > &m)

Detailed Description

template<typename TYPE>
struct Hooke< TYPE >

Hook equation.

Definition at line 19 of file Hooke.hpp.

Constructor & Destructor Documentation

Hooke()

template<typename TYPE >

Hooke< TYPE >::Hooke ( )

inline

Definition at line 23 of file Hooke.hpp.

      : NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI<TYPE>() {}

Member Function Documentation

calculateElasticEnergy()

template<typename TYPE >

virtual MoFEMErrorCode Hooke< TYPE >::calculateElasticEnergy ( const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr )

inlinevirtual

calculate density of strain energy

\(\Psi = \frac{1}{2}\lambda(\textrm{tr}[\varepsilon])^2+\mu\varepsilon:\varepsilon\)

Reimplemented from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >.

Definition at line 65 of file Hooke.hpp.

                                                             {
    MoFEMFunctionBeginHot;
    this->lambda = LAMBDA(block_data.E, block_data.PoissonRatio);
    this->mu = MU(block_data.E, block_data.PoissonRatio);
    ePs.resize(3, 3);
    noalias(ePs) = this->F;
    for (int dd = 0; dd < 3; dd++) {
      ePs(dd, dd) -= 1;
    }
    ePs += trans(ePs);
    ePs *= 0.5;
    this->eNergy = 0;
    tR = 0;
    for (int dd = 0; dd < 3; dd++) {
      tR += ePs(dd, dd);
      for (int jj = 0; jj < 3; jj++) {
        this->eNergy += this->mu * ePs(dd, jj) * ePs(dd, jj);
      }
    }
    this->eNergy += 0.5 * this->lambda * tR * tR;
    MoFEMFunctionReturnHot(0);
  }

calculateP_PiolaKirchhoffI()

template<typename TYPE >

virtual MoFEMErrorCode Hooke< TYPE >::calculateP_PiolaKirchhoffI ( const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr )

inlinevirtual

Hooke equation.

\(\sigma = \lambda\textrm{tr}[\varepsilon]+2\mu\varepsilon\)

Reimplemented from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >.

Definition at line 34 of file Hooke.hpp.

                                                             {
    MoFEMFunctionBeginHot;
    this->lambda = LAMBDA(block_data.E, block_data.PoissonRatio);
    this->mu = MU(block_data.E, block_data.PoissonRatio);
    ePs.resize(3, 3, false);
    noalias(ePs) = this->F;
    for (int dd = 0; dd < 3; dd++) {
      ePs(dd, dd) -= 1;
    }
    ePs += trans(ePs);
    ePs *= 0.5;
    this->P.resize(3, 3, false);
    noalias(this->P) = 2 * this->mu * ePs;
    tR = 0;
    for (int dd = 0; dd < 3; dd++) {
      tR += ePs(dd, dd);
    }
    for (int dd = 0; dd < 3; dd++) {
      this->P(dd, dd) += this->lambda * tR;
    }
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

ePs

template<typename TYPE >

ublas::matrix<TYPE> Hooke< TYPE >::ePs

Definition at line 26 of file Hooke.hpp.

tR

template<typename TYPE >

TYPE Hooke< TYPE >::tR

Definition at line 27 of file Hooke.hpp.

---

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

• users_modules/basic_finite_elements/nonlinear_elastic_materials/src/Hooke.hpp