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. More...

virtual MoFEMErrorCode	calculateElasticEnergy (const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	calculate density of strain energy More...

Public Member Functions inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
	FunctionsToCalculatePiolaKirchhoffI ()

virtual	~FunctionsToCalculatePiolaKirchhoffI ()=default

MoFEMErrorCode	calculateC_CauchyDeformationTensor ()

MoFEMErrorCode	calculateE_GreenStrain ()

MoFEMErrorCode	calculateS_PiolaKirchhoffII ()

virtual MoFEMErrorCode	calculateP_PiolaKirchhoffI (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Function overload to implement user material. More...

virtual MoFEMErrorCode	calculateCauchyStress (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Function overload to implement user material. More...

virtual MoFEMErrorCode	setUserActiveVariables (int &nb_active_variables)
	add additional active variables More...

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. More...

virtual MoFEMErrorCode	calculateElasticEnergy (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Calculate elastic energy density. More...

virtual MoFEMErrorCode	calculatesIGma_EshelbyStress (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Calculate Eshelby stress. More...

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. More...

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. More...

CommonData *	commonDataPtr

MoFEM::VolumeElementForcesAndSourcesCore::UserDataOperator *	opPtr
	pointer to finite element tetrahedral operator More...

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.

Examples: cell_forces.cpp, and simple_contact.cpp.

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.

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

NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI

Implementation of elastic (non-linear) St. Kirchhoff equation.

Definition: NonLinearElasticElement.hpp:130

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

Public Member Functions

Public Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Hooke()

Member Function Documentation

◆ calculateElasticEnergy()

◆ calculateP_PiolaKirchhoffI()

Member Data Documentation

◆ ePs

◆ tR