NonlinearElasticElementInterface Struct Reference

Set of functions declaring elements and setting operators for generic element interface. More...

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

Inheritance diagram for NonlinearElasticElementInterface:

Collaboration diagram for NonlinearElasticElementInterface:

Public Member Functions
	NonlinearElasticElementInterface (MoFEM::Interface &m_field, string postion_field, string mesh_posi_field_name="MESH_NODE_POSITIONS", bool is_displacement_field=true, PetscBool is_quasi_static=PETSC_TRUE)
	~NonlinearElasticElementInterface ()
MoFEMErrorCode	getCommandLineParameters ()
MoFEMErrorCode	addElementFields ()
MoFEMErrorCode	createElements ()
MoFEMErrorCode	setOperators ()
BitRefLevel	getBitRefLevel ()
MoFEMErrorCode	addElementsToDM (SmartPetscObj< DM > dm)
MoFEMErrorCode	setupSolverJacobianSNES ()
MoFEMErrorCode	setupSolverFunctionSNES ()
MoFEMErrorCode	setupSolverJacobianTS (const TSType type)
MoFEMErrorCode	setupSolverFunctionTS (const TSType type)
MoFEMErrorCode	updateElementVariables ()
MoFEMErrorCode	postProcessElement (int step)
Public Member Functions inherited from GenericElementInterface
	GenericElementInterface ()
virtual	~GenericElementInterface ()
virtual MoFEMErrorCode	setGlobalBoundaryMarker (BcMarkerPtr mark)
virtual BcMarkerPtr	getGlobalBoundaryMarker ()
virtual MoFEMErrorCode	setMonitorPtr (boost::shared_ptr< MoFEM::FEMethod > monitor_ptr)
virtual BitRefLevel	getBitRefLevelMask ()
virtual MoFEMErrorCode	opFactoryDomainRhs (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip)
virtual MoFEMErrorCode	opFactoryDomainLhs (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip)
virtual MoFEMErrorCode	setUpdateElementVariablesOperators ()
virtual MoFEMErrorCode	updateElementVariables (SmartPetscObj< DM > dm, string fe_name)
virtual MoFEMErrorCode	postProcessElement (int step, SmartPetscObj< DM > dm, string fe_name)
	GenericElementInterface ()
virtual	~GenericElementInterface ()
virtual MoFEMErrorCode	setGlobalBoundaryMarker (BcMarkerPtr mark)
virtual BcMarkerPtr	getGlobalBoundaryMarker ()
virtual MoFEMErrorCode	setMonitorPtr (boost::shared_ptr< FEMethod > monitor_ptr)
virtual BitRefLevel	getBitRefLevelMask ()
virtual MoFEMErrorCode	opFactoryDomainRhs (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip)
virtual MoFEMErrorCode	opFactoryDomainLhs (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip)
virtual MoFEMErrorCode	setUpdateElementVariablesOperators ()
virtual MoFEMErrorCode	updateElementVariables (SmartPetscObj< DM > dm, string fe_name)
virtual MoFEMErrorCode	postProcessElement (int step, SmartPetscObj< DM > dm, string fe_name)

Public Attributes
MoFEM::Interface &	mField
SmartPetscObj< DM >	dM
PetscBool	isQuasiStatic
PetscInt	oRder
bool	isDisplacementField
BitRefLevel	bIt
boost::shared_ptr< NonlinearElasticElement >	elasticElementPtr
boost::shared_ptr< ElasticMaterials >	elasticMaterialsPtr
boost::shared_ptr< PostProcVolumeOnRefinedMesh >	postProcMeshPtr
string	positionField
string	meshNodeField
Public Attributes inherited from GenericElementInterface
BcMarkerPtr	mBoundaryMarker
int	atomTest
int	restartRunStep
boost::shared_ptr< MoFEM::FEMethod >	monitorPtr
boost::shared_ptr< FEMethod >	monitorPtr

Additional Inherited Members
Public Types inherited from GenericElementInterface
enum	TSType {   EX , IM , IM2 , IMEX ,   DEFAULT , EX , IM , IM2 ,   IMEX , DEFAULT }
enum	TSType {   EX , IM , IM2 , IMEX ,   DEFAULT , EX , IM , IM2 ,   IMEX , DEFAULT }
using	BcMarkerPtr = boost::shared_ptr<std::vector<char unsigned>>

Detailed Description

Set of functions declaring elements and setting operators for generic element interface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 15 of file NonlinearElasticElementInterface.hpp.

Constructor & Destructor Documentation

NonlinearElasticElementInterface()

NonlinearElasticElementInterface::NonlinearElasticElementInterface ( MoFEM::Interface & m_field, string postion_field, string mesh_posi_field_name = "MESH_NODE_POSITIONS", bool is_displacement_field = true, PetscBool is_quasi_static = PETSC_TRUE )

inline

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 31 of file NonlinearElasticElementInterface.hpp.

      : mField(m_field), positionField(postion_field),
        meshNodeField(mesh_posi_field_name),
        isDisplacementField(is_displacement_field),
        isQuasiStatic(is_quasi_static) {
    oRder = 1;
  }

~NonlinearElasticElementInterface()

NonlinearElasticElementInterface::~NonlinearElasticElementInterface ( )

inline

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 43 of file NonlinearElasticElementInterface.hpp.

{}

Member Function Documentation

addElementFields()

MoFEMErrorCode NonlinearElasticElementInterface::addElementFields ( )

inlinevirtual

Implements GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 56 of file NonlinearElasticElementInterface.hpp.

                                    {
    MoFEMFunctionBeginHot;
    auto simple = mField.getInterface<Simple>();
    if (!mField.check_field(positionField)) {
      CHKERR simple->addDomainField(positionField, H1, AINSWORTH_LEGENDRE_BASE,
                                    3);
      CHKERR simple->addBoundaryField(positionField, H1,
                                      AINSWORTH_LEGENDRE_BASE, 3);
      CHKERR simple->setFieldOrder(positionField, oRder);
    }
    if (!mField.check_field(meshNodeField)) {
      CHKERR simple->addDataField(meshNodeField, H1, AINSWORTH_LEGENDRE_BASE,
                                  3);
      CHKERR simple->setFieldOrder(meshNodeField, 2);
    }
 
    MoFEMFunctionReturnHot(0);
  };

addElementsToDM()

MoFEMErrorCode NonlinearElasticElementInterface::addElementsToDM ( SmartPetscObj< DM > dm )

inlinevirtual

Implements GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 112 of file NonlinearElasticElementInterface.hpp.

                                                       {
    MoFEMFunctionBeginHot;
    this->dM = dm;
    CHKERR DMMoFEMAddElement(dM, "ELASTIC");
    mField.getInterface<Simple>()->getOtherFiniteElements().push_back(
        "ELASTIC");
        
    MoFEMFunctionReturnHot(0);
  };

createElements()

MoFEMErrorCode NonlinearElasticElementInterface::createElements ( )

inlinevirtual

Implements GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 75 of file NonlinearElasticElementInterface.hpp.

                                  {
    MoFEMFunctionBeginHot;
 
    elasticElementPtr = boost::make_shared<NonlinearElasticElement>(mField, 2);
    elasticMaterialsPtr = boost::make_shared<ElasticMaterials>(mField, "HOOKE");
    CHKERR elasticMaterialsPtr->setBlocks(elasticElementPtr->setOfBlocks);
 
    CHKERR addHOOpsVol(meshNodeField, elasticElementPtr->getLoopFeRhs(), true,
                       false, false, false);
    CHKERR addHOOpsVol(meshNodeField, elasticElementPtr->getLoopFeLhs(), true,
                       false, false, false);
    CHKERR addHOOpsVol(meshNodeField, elasticElementPtr->getLoopFeEnergy(),
                       true, false, false, false);
    CHKERR elasticElementPtr->addElement("ELASTIC", positionField,
                                         meshNodeField, false);
 
 
    MoFEMFunctionReturnHot(0);
  };

getBitRefLevel()

BitRefLevel NonlinearElasticElementInterface::getBitRefLevel ( )

inlinevirtual

Reimplemented from GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 111 of file NonlinearElasticElementInterface.hpp.

{ return bIt; };

getCommandLineParameters()

MoFEMErrorCode NonlinearElasticElementInterface::getCommandLineParameters ( )

inlinevirtual

Reimplemented from GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 45 of file NonlinearElasticElementInterface.hpp.

                                            {
    MoFEMFunctionBegin;
    isQuasiStatic = PETSC_FALSE;
    oRder = 2;
    CHKERR PetscOptionsGetBool(PETSC_NULLPTR, "-is_quasi_static", &isQuasiStatic,
                               PETSC_NULLPTR);
    CHKERR PetscOptionsGetInt(PETSC_NULLPTR, "-order", &oRder, PETSC_NULLPTR);
 
    MoFEMFunctionReturn(0);
  };

postProcessElement()

MoFEMErrorCode NonlinearElasticElementInterface::postProcessElement ( int step )

inlinevirtual

Implements GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 181 of file NonlinearElasticElementInterface.hpp.

                                              {
    MoFEMFunctionBegin;
    
    if (elasticElementPtr->setOfBlocks.empty())
      MoFEMFunctionReturnHot(0);
 
    if (!postProcMeshPtr) {
      postProcMeshPtr = boost::make_shared<PostProcVolumeOnRefinedMesh>(mField);
 
      if (mField.check_field("MESH_NODE_POSITIONS"))
        CHKERR addHOOpsVol("MESH_NODE_POSITIONS", *postProcMeshPtr, true, false,
                           false, false);
      CHKERR postProcMeshPtr->generateReferenceElementMesh();
      CHKERR postProcMeshPtr->addFieldValuesPostProc(positionField);
      CHKERR postProcMeshPtr->addFieldValuesPostProc(meshNodeField);
      CHKERR postProcMeshPtr->addFieldValuesGradientPostProc(positionField);
 
      for (auto &sit : elasticElementPtr->setOfBlocks) {
        postProcMeshPtr->getOpPtrVector().push_back(new PostProcStress(
            postProcMeshPtr->postProcMesh, postProcMeshPtr->mapGaussPts,
            positionField, sit.second, postProcMeshPtr->commonData,
            isDisplacementField));
      }
    }
 
    elasticElementPtr->getLoopFeEnergy().snes_ctx = SnesMethod::CTX_SNESNONE;
    elasticElementPtr->getLoopFeEnergy().eNergy = 0;
    // MOFEM_LOG("WORLD", Sev::inform) << "Loop energy\n";
    CHKERR DMoFEMLoopFiniteElements(dM, "ELASTIC",
                                    &elasticElementPtr->getLoopFeEnergy());
 
    auto E = elasticElementPtr->getLoopFeEnergy().eNergy;
    // Print elastic energy
    MOFEM_LOG_C("WORLD", Sev::inform, "%d Time %3.2e Elastic energy %3.2e",
                step, elasticElementPtr->getLoopFeRhs().ts_t, E);
 
    CHKERR DMoFEMLoopFiniteElements(dM, "ELASTIC", postProcMeshPtr);
    auto out_name = "out_vol_" + to_string(step) + ".h5m";
 
    CHKERR postProcMeshPtr->writeFile(out_name);
 
    MoFEMFunctionReturn(0);
  };

setOperators()

MoFEMErrorCode NonlinearElasticElementInterface::setOperators ( )

inlinevirtual

Implements GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 95 of file NonlinearElasticElementInterface.hpp.

                                {
    MoFEMFunctionBegin;
    auto &pipeline_rhs = elasticElementPtr->feRhs.getOpPtrVector();
    auto &pipeline_lhs = elasticElementPtr->feLhs.getOpPtrVector();
 
    pipeline_rhs.push_back(new OpSetBc(positionField, true, mBoundaryMarker));
    pipeline_lhs.push_back(new OpSetBc(positionField, true, mBoundaryMarker));
 
    CHKERR elasticElementPtr->setOperators(positionField, meshNodeField, false,
                                           isDisplacementField);
                                           
    pipeline_rhs.push_back(new OpUnSetBc(positionField));
    pipeline_lhs.push_back(new OpUnSetBc(positionField));
    MoFEMFunctionReturn(0);
  }

setupSolverFunctionSNES()

MoFEMErrorCode NonlinearElasticElementInterface::setupSolverFunctionSNES ( )

inlinevirtual

Reimplemented from GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 130 of file NonlinearElasticElementInterface.hpp.

                                           {
    MoFEMFunctionBegin;
    CHKERR DMMoFEMSNESSetFunction(dM, "ELASTIC",
                                  &elasticElementPtr->getLoopFeRhs(),
                                  PETSC_NULLPTR, PETSC_NULLPTR);
    MoFEMFunctionReturn(0);
  };

setupSolverFunctionTS()

MoFEMErrorCode NonlinearElasticElementInterface::setupSolverFunctionTS ( const TSType type )

inlinevirtual

Implements GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 159 of file NonlinearElasticElementInterface.hpp.

                                                          {
    MoFEMFunctionBegin;
    auto &method = elasticElementPtr->getLoopFeRhs();
    switch (type) {
    case IM:
      CHKERR DMMoFEMTSSetIFunction(dM, "ELASTIC", &method, &method, &method);
      break;
    case IM2:
      CHKERR DMMoFEMTSSetI2Function(dM, "ELASTIC", &method, &method, &method);
      break;
    case EX:
      CHKERR DMMoFEMTSSetRHSFunction(dM, "ELASTIC", &method, &method, &method);
      break;
    default:
      SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "Not implemented");
      break;
    }
 
    MoFEMFunctionReturn(0);
  };

setupSolverJacobianSNES()

MoFEMErrorCode NonlinearElasticElementInterface::setupSolverJacobianSNES ( )

inlinevirtual

Reimplemented from GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 122 of file NonlinearElasticElementInterface.hpp.

                                           {
    MoFEMFunctionBegin;
 
    CHKERR DMMoFEMSNESSetJacobian(
        dM, "ELASTIC", &elasticElementPtr->getLoopFeLhs(), NULL, NULL);
 
    MoFEMFunctionReturn(0);
  };

setupSolverJacobianTS()

MoFEMErrorCode NonlinearElasticElementInterface::setupSolverJacobianTS ( const TSType type )

inlinevirtual

Implements GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 138 of file NonlinearElasticElementInterface.hpp.

                                                          {
    MoFEMFunctionBegin;
    auto &method = elasticElementPtr->getLoopFeLhs();
    switch (type) {
    case IM:
      CHKERR DMMoFEMTSSetIJacobian(dM, "ELASTIC", &method, &method, &method);
      break;
    case IM2:
      CHKERR DMMoFEMTSSetI2Jacobian(dM, "ELASTIC", &method, &method, &method);
      break;
    case EX:
      CHKERR DMMoFEMTSSetRHSJacobian(dM, "ELASTIC", &method, &method, &method);
      break;
    default:
      SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
              "This TS is not yet implemented");
      break;
    }
    MoFEMFunctionReturn(0);
  };

updateElementVariables()

MoFEMErrorCode NonlinearElasticElementInterface::updateElementVariables ( )

inlinevirtual

Reimplemented from GenericElementInterface.

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 180 of file NonlinearElasticElementInterface.hpp.

{ return 0; };

Member Data Documentation

bIt

BitRefLevel NonlinearElasticElementInterface::bIt

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 23 of file NonlinearElasticElementInterface.hpp.

dM

SmartPetscObj<DM> NonlinearElasticElementInterface::dM

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 18 of file NonlinearElasticElementInterface.hpp.

elasticElementPtr

boost::shared_ptr<NonlinearElasticElement> NonlinearElasticElementInterface::elasticElementPtr

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 24 of file NonlinearElasticElementInterface.hpp.

elasticMaterialsPtr

boost::shared_ptr<ElasticMaterials> NonlinearElasticElementInterface::elasticMaterialsPtr

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 25 of file NonlinearElasticElementInterface.hpp.

isDisplacementField

bool NonlinearElasticElementInterface::isDisplacementField

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 22 of file NonlinearElasticElementInterface.hpp.

isQuasiStatic

PetscBool NonlinearElasticElementInterface::isQuasiStatic

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 19 of file NonlinearElasticElementInterface.hpp.

meshNodeField

string NonlinearElasticElementInterface::meshNodeField

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 29 of file NonlinearElasticElementInterface.hpp.

mField

MoFEM::Interface& NonlinearElasticElementInterface::mField

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 17 of file NonlinearElasticElementInterface.hpp.

oRder

PetscInt NonlinearElasticElementInterface::oRder

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 21 of file NonlinearElasticElementInterface.hpp.

positionField

string NonlinearElasticElementInterface::positionField

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 28 of file NonlinearElasticElementInterface.hpp.

postProcMeshPtr

boost::shared_ptr<PostProcVolumeOnRefinedMesh> NonlinearElasticElementInterface::postProcMeshPtr

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 26 of file NonlinearElasticElementInterface.hpp.

---

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

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