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< PostProcBrokenMeshInMoab< VolumeElementForcesAndSourcesCore > >	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 > >
using	BcMarkerPtr = boost::shared_ptr< std::vector< char unsigned > >

Detailed Description

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

Examples

mofem/users_modules/basic_finite_elements/nonlinear_elastic_materials/src/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

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

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

Definition at line 45 of file NonlinearElasticElementInterface.hpp.

{}

Member Function Documentation

addElementFields()

MoFEMErrorCode NonlinearElasticElementInterface::addElementFields ( )

inlinevirtual

Implements GenericElementInterface.

Examples

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

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

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

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

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

Definition at line 77 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 AddHOOps<3, 3, 3>::add(
        elasticElementPtr->getLoopFeRhs().getOpPtrVector(), {H1},
        meshNodeField);
    CHKERR AddHOOps<3, 3, 3>::add(
        elasticElementPtr->getLoopFeLhs().getOpPtrVector(), {H1},
        meshNodeField);
    CHKERR AddHOOps<3, 3, 3>::add(
        elasticElementPtr->getLoopFeEnergy().getOpPtrVector(), {H1},
        meshNodeField);
    CHKERR elasticElementPtr->addElement("ELASTIC", positionField,
                                         meshNodeField, false);
 
 
    MoFEMFunctionReturnHot(0);
  };

getBitRefLevel()

BitRefLevel NonlinearElasticElementInterface::getBitRefLevel ( )

inlinevirtual

Reimplemented from GenericElementInterface.

Examples

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

Definition at line 116 of file NonlinearElasticElementInterface.hpp.

{ return bIt; };

getCommandLineParameters()

MoFEMErrorCode NonlinearElasticElementInterface::getCommandLineParameters ( )

inlinevirtual

Reimplemented from GenericElementInterface.

Examples

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

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

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

Definition at line 186 of file NonlinearElasticElementInterface.hpp.

                                              {
    MoFEMFunctionBegin;
    
    if (elasticElementPtr->setOfBlocks.empty())
      MoFEMFunctionReturnHot(0);
 
    if (!postProcMeshPtr) {
      postProcMeshPtr = boost::make_shared<
          PostProcBrokenMeshInMoab<VolumeElementForcesAndSourcesCore>>(mField);
 
      if (mField.check_field("MESH_NODE_POSITIONS"))
        CHKERR AddHOOps<3, 3, 3>::add(postProcMeshPtr->getOpPtrVector(), {H1},
                                      "MESH_NODE_POSITIONS");
      auto position_ptr = boost::make_shared<MatrixDouble>();
      auto mesh_pos_ptr = boost::make_shared<MatrixDouble>();
      auto position_grad_ptr = boost::make_shared<MatrixDouble>();
      postProcMeshPtr->getOpPtrVector().push_back(
          new OpCalculateVectorFieldValues<3>(positionField, position_ptr));
      postProcMeshPtr->getOpPtrVector().push_back(
          new OpCalculateVectorFieldValues<3>(meshNodeField, mesh_pos_ptr));
      postProcMeshPtr->getOpPtrVector().push_back(
          new OpCalculateVectorFieldGradient<3, 3>(positionField,
                                                   position_grad_ptr));
 
      for (auto &sit : elasticElementPtr->setOfBlocks) {
        postProcMeshPtr->getOpPtrVector().push_back(new PostProcStress(
            postProcMeshPtr->getPostProcMesh(),
            postProcMeshPtr->getMapGaussPts(),
            postProcMeshPtr->getPostProcElements(), positionField, sit.second,
            position_ptr, mesh_pos_ptr, position_grad_ptr, meshNodeField,
            isDisplacementField));
      }
      using OpPPMap = OpPostProcMapInMoab<3, 3>;
      postProcMeshPtr->getOpPtrVector().push_back(new OpPPMap(
          postProcMeshPtr->getPostProcMesh(),
          postProcMeshPtr->getMapGaussPts(), {},
          {{positionField, position_ptr}, {meshNodeField, mesh_pos_ptr}},
          {{positionField + "_GRAD", position_grad_ptr}}, {}));
    }
 
    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

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

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

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

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

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

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

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

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

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

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

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

Definition at line 185 of file NonlinearElasticElementInterface.hpp.

{ return 0; };

Member Data Documentation

bIt

BitRefLevel NonlinearElasticElementInterface::bIt

Examples

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

Definition at line 23 of file NonlinearElasticElementInterface.hpp.

dM

SmartPetscObj<DM> NonlinearElasticElementInterface::dM

Examples

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

Definition at line 18 of file NonlinearElasticElementInterface.hpp.

elasticElementPtr

boost::shared_ptr<NonlinearElasticElement> NonlinearElasticElementInterface::elasticElementPtr

Examples

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

Definition at line 24 of file NonlinearElasticElementInterface.hpp.

elasticMaterialsPtr

boost::shared_ptr<ElasticMaterials> NonlinearElasticElementInterface::elasticMaterialsPtr

Examples

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

Definition at line 25 of file NonlinearElasticElementInterface.hpp.

isDisplacementField

bool NonlinearElasticElementInterface::isDisplacementField

Examples

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

Definition at line 22 of file NonlinearElasticElementInterface.hpp.

isQuasiStatic

PetscBool NonlinearElasticElementInterface::isQuasiStatic

Examples

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

Definition at line 19 of file NonlinearElasticElementInterface.hpp.

meshNodeField

string NonlinearElasticElementInterface::meshNodeField

Examples

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

Definition at line 31 of file NonlinearElasticElementInterface.hpp.

mField

MoFEM::Interface& NonlinearElasticElementInterface::mField

Examples

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

Definition at line 17 of file NonlinearElasticElementInterface.hpp.

oRder

PetscInt NonlinearElasticElementInterface::oRder

Examples

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

Definition at line 21 of file NonlinearElasticElementInterface.hpp.

positionField

string NonlinearElasticElementInterface::positionField

Examples

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

Definition at line 30 of file NonlinearElasticElementInterface.hpp.

postProcMeshPtr

boost::shared_ptr< PostProcBrokenMeshInMoab<VolumeElementForcesAndSourcesCore> > NonlinearElasticElementInterface::postProcMeshPtr

Examples

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

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