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

Public Attributes
DEPRECATED typedef OpRhsEshelbyStress	OpRhsEshelbyStrees

MyVolumeFE	feRhs
	calculate right hand side for tetrahedral elements More...

MyVolumeFE	feLhs

MyVolumeFE	feEnergy
	calculate elastic energy More...

MoFEM::Interface &	mField

short int	tAg

std::map< int, BlockData >	setOfBlocks
	maps block set id with appropriate BlockData More...

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: cell_forces.cpp, HookeElement.hpp, nonlinear_dynamics.cpp, Remodeling.cpp, and simple_contact.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 81 of file NonLinearElasticElement.cpp.

83 : feRhs(m_field), feLhs(m_field), feEnergy(m_field), mField(m_field),

84 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`
	)

Examples: simple_contact.cpp.

Definition at line 1120 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);
 }

◆ setBlocks()

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

Examples: simple_contact.cpp.

Definition at line 1086 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

Examples: simple_contact.cpp.

Definition at line 1153 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);
 }