NonlinearElasticElement Struct Reference

structure grouping operators and data used for calculation of nonlinear elastic element More...

#include "users_modules/basic_finite_elements/src/HookeElement.hpp"

Collaboration diagram for NonlinearElasticElement:

Classes
struct	BlockData
	data for calculation heat conductivity and heat capacity elements More...
struct	CommonData
	common data used by volume elements More...
struct	FunctionsToCalculatePiolaKirchhoffI
	Implementation of elastic (non-linear) St. Kirchhoff equation. More...
struct	MyVolumeFE
	definition of volume element More...
struct	OpEnergy
struct	OpGetCommonDataAtGaussPts
struct	OpGetDataAtGaussPts
struct	OpJacobianEnergy
	Calculate explicit derivative of free energy. More...
struct	OpJacobianEshelbyStress
struct	OpJacobianPiolaKirchhoffStress
	Operator performs automatic differentiation. More...
struct	OpLhsEshelby_dx
struct	OpLhsEshelby_dX
struct	OpLhsPiolaKirchhoff_dx
struct	OpLhsPiolaKirchhoff_dX
struct	OpRhsEshelbyStress
struct	OpRhsPiolaKirchhoff

Public Member Functions
MyVolumeFE &	getLoopFeRhs ()
	get rhs volume element
MyVolumeFE &	getLoopFeLhs ()
	get lhs volume element
MyVolumeFE &	getLoopFeEnergy ()
	get energy fe
	NonlinearElasticElement (MoFEM::Interface &m_field, short int tag)
virtual	~NonlinearElasticElement ()=default
MoFEMErrorCode	setBlocks (boost::shared_ptr< FunctionsToCalculatePiolaKirchhoffI< double > > materialDoublePtr, boost::shared_ptr< FunctionsToCalculatePiolaKirchhoffI< adouble > > materialAdoublePtr)
MoFEMErrorCode	addElement (const std::string element_name, const std::string spatial_position_field_name, const std::string material_position_field_name="MESH_NODE_POSITIONS", const bool ale=false)
MoFEMErrorCode	setOperators (const std::string spatial_position_field_name, const std::string material_position_field_name="MESH_NODE_POSITIONS", const bool ale=false, const bool field_disp=false)
	Set operators to calculate left hand tangent matrix and right hand residual.

Public Attributes
DEPRECATED typedef OpRhsEshelbyStress	OpRhsEshelbyStrees
MyVolumeFE	feRhs
	calculate right hand side for tetrahedral elements
MyVolumeFE	feLhs
MyVolumeFE	feEnergy
	calculate elastic energy
MoFEM::Interface &	mField
short int	tAg
std::map< int, BlockData >	setOfBlocks
	maps block set id with appropriate BlockData
CommonData	commonData
FTensor::Index< 'i', 3 >	i
FTensor::Index< 'j', 3 >	j
FTensor::Index< 'k', 3 >	k

Detailed Description

structure grouping operators and data used for calculation of nonlinear elastic element

In order to assemble matrices and right hand vectors, the loops over elements, entities over that elements and finally loop over integration points are executed.

Following implementation separate those three categories of loops and to each loop attach operator.

Examples

mofem/users_modules/basic_finite_elements/nonlinear_elasticity/nonlinear_dynamics.cpp, mofem/users_modules/basic_finite_elements/src/HookeElement.hpp, and mofem/users_modules/bone_remodelling/src/impl/Remodeling.cpp.

Definition at line 27 of file HookeElement.hpp.

Constructor & Destructor Documentation

NonlinearElasticElement()

NonlinearElasticElement::NonlinearElasticElement	(	MoFEM::Interface &	m_field,
		short int	tag
	)

Definition at line 78 of file NonLinearElasticElement.cpp.

    : feRhs(m_field), feLhs(m_field), feEnergy(m_field), mField(m_field),
      tAg(tag) {}

~NonlinearElasticElement()

virtual NonlinearElasticElement::~NonlinearElasticElement ( )

virtualdefault

Member Function Documentation

addElement()

MoFEMErrorCode NonlinearElasticElement::addElement	(	const std::string	element_name,
		const std::string	spatial_position_field_name,
		const std::string	material_position_field_name = "MESH_NODE_POSITIONS",
		const bool	ale = false
	)

Definition at line 1118 of file NonLinearElasticElement.cpp.

                                                                  {
  MoFEMFunctionBegin;
 
  CHKERR mField.add_finite_element(element_name, MF_ZERO);
  CHKERR mField.modify_finite_element_add_field_row(
      element_name, spatial_position_field_name);
  CHKERR mField.modify_finite_element_add_field_col(
      element_name, spatial_position_field_name);
  CHKERR mField.modify_finite_element_add_field_data(
      element_name, spatial_position_field_name);
  if (mField.check_field(material_position_field_name)) {
    if (ale) {
      CHKERR mField.modify_finite_element_add_field_row(
          element_name, material_position_field_name);
      CHKERR mField.modify_finite_element_add_field_col(
          element_name, material_position_field_name);
    }
    CHKERR mField.modify_finite_element_add_field_data(
        element_name, material_position_field_name);
  }
 
  std::map<int, BlockData>::iterator sit = setOfBlocks.begin();
  for (; sit != setOfBlocks.end(); sit++) {
    CHKERR mField.add_ents_to_finite_element_by_type(sit->second.tEts, MBTET,
                                                     element_name);
  }
 
  MoFEMFunctionReturn(0);
}

getLoopFeEnergy()

MyVolumeFE & NonlinearElasticElement::getLoopFeEnergy ( )

inline

get energy fe

Definition at line 71 of file NonLinearElasticElement.hpp.

getLoopFeLhs()

MyVolumeFE & NonlinearElasticElement::getLoopFeLhs ( )

inline

get lhs volume element

Definition at line 68 of file NonLinearElasticElement.hpp.

getLoopFeRhs()

MyVolumeFE & NonlinearElasticElement::getLoopFeRhs ( )

inline

get rhs volume element

Definition at line 66 of file NonLinearElasticElement.hpp.

setBlocks()

MoFEMErrorCode NonlinearElasticElement::setBlocks	(	boost::shared_ptr< FunctionsToCalculatePiolaKirchhoffI< double > >	materialDoublePtr,
		boost::shared_ptr< FunctionsToCalculatePiolaKirchhoffI< adouble > >	materialAdoublePtr
	)

Definition at line 1084 of file NonLinearElasticElement.cpp.

                            {
  MoFEMFunctionBegin;
 
  if (!materialDoublePtr) {
    SETERRQ(mField.get_comm(), MOFEM_DATA_INCONSISTENCY,
            "Pointer for materialDoublePtr not allocated");
  }
  if (!materialAdoublePtr) {
    SETERRQ(mField.get_comm(), MOFEM_DATA_INCONSISTENCY,
            "Pointer for materialAdoublePtr not allocated");
  }
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
           mField, BLOCKSET | MAT_ELASTICSET, it)) {
    Mat_Elastic mydata;
    CHKERR it->getAttributeDataStructure(mydata);
    int id = it->getMeshsetId();
    EntityHandle meshset = it->getMeshset();
    CHKERR mField.get_moab().get_entities_by_type(meshset, MBTET,
                                                  setOfBlocks[id].tEts, true);
    setOfBlocks[id].iD = id;
    setOfBlocks[id].E = mydata.data.Young;
    setOfBlocks[id].PoissonRatio = mydata.data.Poisson;
    setOfBlocks[id].materialDoublePtr = materialDoublePtr;
    setOfBlocks[id].materialAdoublePtr = materialAdoublePtr;
  }
 
  MoFEMFunctionReturn(0);
}

setOperators()

MoFEMErrorCode NonlinearElasticElement::setOperators	(	const std::string	spatial_position_field_name,
		const std::string	material_position_field_name = "MESH_NODE_POSITIONS",
		const bool	ale = false,
		const bool	field_disp = false
	)

Set operators to calculate left hand tangent matrix and right hand residual.

Parameters

fun	class needed to calculate Piola Kirchhoff I Stress tensor
spatial_position_field_name	name of approximation field
material_position_field_name	name of field to define geometry
ale	true if arbitrary Lagrangian Eulerian formulation
field_disp	true if approximation field represents displacements otherwise it is field of spatial positions

Definition at line 1151 of file NonLinearElasticElement.cpp.

                           {
  MoFEMFunctionBegin;
 
  commonData.spatialPositions = spatial_position_field_name;
  commonData.meshPositions = material_position_field_name;
 
  // Rhs
  feRhs.getOpPtrVector().push_back(
      new OpGetCommonDataAtGaussPts(spatial_position_field_name, commonData));
  if (mField.check_field(material_position_field_name)) {
    feRhs.getOpPtrVector().push_back(new OpGetCommonDataAtGaussPts(
        material_position_field_name, commonData));
  }
  std::map<int, BlockData>::iterator sit = setOfBlocks.begin();
  for (; sit != setOfBlocks.end(); sit++) {
    feRhs.getOpPtrVector().push_back(new OpJacobianPiolaKirchhoffStress(
        spatial_position_field_name, sit->second, commonData, tAg, false, ale,
        field_disp));
    feRhs.getOpPtrVector().push_back(new OpRhsPiolaKirchhoff(
        spatial_position_field_name, sit->second, commonData));
  }
 
  // Energy
  feEnergy.getOpPtrVector().push_back(
      new OpGetCommonDataAtGaussPts(spatial_position_field_name, commonData));
  if (mField.check_field(material_position_field_name)) {
    feEnergy.getOpPtrVector().push_back(new OpGetCommonDataAtGaussPts(
        material_position_field_name, commonData));
  }
  sit = setOfBlocks.begin();
  for (; sit != setOfBlocks.end(); sit++) {
    feEnergy.getOpPtrVector().push_back(
        new OpEnergy(spatial_position_field_name, sit->second, commonData,
                     feEnergy.V, field_disp));
  }
 
  // Lhs
  feLhs.getOpPtrVector().push_back(
      new OpGetCommonDataAtGaussPts(spatial_position_field_name, commonData));
  if (mField.check_field(material_position_field_name)) {
    feLhs.getOpPtrVector().push_back(new OpGetCommonDataAtGaussPts(
        material_position_field_name, commonData));
  }
  sit = setOfBlocks.begin();
  for (; sit != setOfBlocks.end(); sit++) {
    feLhs.getOpPtrVector().push_back(new OpJacobianPiolaKirchhoffStress(
        spatial_position_field_name, sit->second, commonData, tAg, true, ale,
        field_disp));
    feLhs.getOpPtrVector().push_back(new OpLhsPiolaKirchhoff_dx(
        spatial_position_field_name, spatial_position_field_name, sit->second,
        commonData));
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonData

CommonData NonlinearElasticElement::commonData

Definition at line 121 of file NonLinearElasticElement.hpp.

feEnergy

MyVolumeFE NonlinearElasticElement::feEnergy

calculate elastic energy

Definition at line 70 of file NonLinearElasticElement.hpp.

feLhs

MyVolumeFE NonlinearElasticElement::feLhs

Definition at line 67 of file NonLinearElasticElement.hpp.

feRhs

MyVolumeFE NonlinearElasticElement::feRhs

calculate right hand side for tetrahedral elements

Definition at line 65 of file NonLinearElasticElement.hpp.

i

FTensor::Index<'i', 3> NonlinearElasticElement::i

Definition at line 123 of file NonLinearElasticElement.hpp.

j

FTensor::Index<'j', 3> NonlinearElasticElement::j

Definition at line 124 of file NonLinearElasticElement.hpp.

k

FTensor::Index<'k', 3> NonlinearElasticElement::k

Definition at line 125 of file NonLinearElasticElement.hpp.

mField

MoFEM::Interface& NonlinearElasticElement::mField

Definition at line 73 of file NonLinearElasticElement.hpp.

OpRhsEshelbyStrees

DEPRECATED typedef OpRhsEshelbyStress NonlinearElasticElement::OpRhsEshelbyStrees

Deprecated:

name with spelling mistake

Definition at line 631 of file NonLinearElasticElement.hpp.

setOfBlocks

std::map<int, BlockData> NonlinearElasticElement::setOfBlocks

maps block set id with appropriate BlockData

Definition at line 100 of file NonLinearElasticElement.hpp.

tAg

short int NonlinearElasticElement::tAg

Definition at line 74 of file NonLinearElasticElement.hpp.

---

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

• users_modules/basic_finite_elements/src/HookeElement.hpp
• users_modules/basic_finite_elements/src/NonLinearElasticElement.hpp
• users_modules/basic_finite_elements/src/impl/NonLinearElasticElement.cpp