ElasticTieMeshExample Struct Reference

ElasticExample extension providing reusable TIE mesh behaviour. More...

#include "tutorials/vec-11_elastic_tie_mesh/src/ElasticTie.hpp"

Inheritance diagram for ElasticTieMeshExample:

Collaboration diagram for ElasticTieMeshExample:

Public Member Functions
	ElasticExample (MoFEM::Interface &m_field)
Public Member Functions inherited from ElasticExample
	ElasticExample (MoFEM::Interface &m_field)
virtual	~ElasticExample ()=default
MoFEMErrorCode	runProblem ()
	[Run problem]

Protected Member Functions
MoFEMErrorCode	readMesh () override
	[Run problem]
MoFEMErrorCode	setupProblem () override
	[Read mesh]
MoFEMErrorCode	setUpSolver (SmartPetscObj< KSP > solver) override
	[Push operators to pipeline]
MoFEMErrorCode	kspSetUpAndSolve (SmartPetscObj< KSP > solver) override
MoFEMErrorCode	checkResults () override
	[Postprocess results]
Protected Member Functions inherited from ElasticExample
virtual MoFEMErrorCode	boundaryCondition ()
	[Set up problem]
virtual MoFEMErrorCode	assembleSystem ()
virtual MoFEMErrorCode	solveSystem ()
	[Solve]
virtual MoFEMErrorCode	outputResults ()
	[Solve]
virtual MoFEMErrorCode	addDomainMaterialOps (Pipeline &lhs_pipeline, Pipeline &rhs_pipeline)
	Add constitutive operators to the domain assembly pipelines.

Protected Attributes
boost::shared_ptr< ElasticTie::CommonData >	tieData = nullptr
Range	physicalBoundaryFaces
Protected Attributes inherited from ElasticExample
MoFEM::Interface &	mField
boost::shared_ptr< MatrixDouble >	vectorFieldPtr = nullptr

Additional Inherited Members
Protected Types inherited from ElasticExample
using	Pipeline = boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator >

Detailed Description

ElasticExample extension providing reusable TIE mesh behaviour.

Examples

mofem/tutorials/vec-11_elastic_tie_mesh/elastic_tie_mesh.cpp.

Definition at line 1254 of file ElasticTie.hpp.

Member Function Documentation

checkResults()

MoFEMErrorCode ElasticTieMeshExample::checkResults ( )

inlineoverrideprotectedvirtual

[Postprocess results]

[Check]

Reimplemented from ElasticExample.

Reimplemented in ElasticTieMeshTutorial.

Definition at line 1355 of file ElasticTie.hpp.

                                                          {
  MoFEMFunctionBegin;
  MoFEMFunctionReturn(0);
}

ElasticExample()

ElasticExample::ElasticExample ( MoFEM::Interface &  m_field )

inline

Definition at line 20 of file ElasticExample.hpp.

: mField(m_field) {}

kspSetUpAndSolve()

MoFEMErrorCode ElasticTieMeshExample::kspSetUpAndSolve ( SmartPetscObj< KSP >  solver )

inlineoverrideprotectedvirtual

Reimplemented from ElasticExample.

Definition at line 1348 of file ElasticTie.hpp.

                                                                 {
  MoFEMFunctionBegin;
  CHKERR ElasticExample::kspSetUpAndSolve(solver);
  CHKERR ElasticTie::logTieDisplacementNorms(mField, tieData, "U");
  MoFEMFunctionReturn(0);
}

readMesh()

MoFEMErrorCode ElasticTieMeshExample::readMesh ( )

inlineoverrideprotectedvirtual

[Run problem]

[Read mesh]

Reimplemented from ElasticExample.

Definition at line 1269 of file ElasticTie.hpp.

                                                      {
  MoFEMFunctionBegin;
  CHKERR ElasticTie::loadTieMesh<SPACE_DIM>(mField, physicalBoundaryFaces);
  MoFEMFunctionReturn(0);
}

setupProblem()

MoFEMErrorCode ElasticTieMeshExample::setupProblem ( )

inlineoverrideprotectedvirtual

[Read mesh]

[Set up problem]

Reimplemented from ElasticExample.

Definition at line 1275 of file ElasticTie.hpp.

                                                          {
  MoFEMFunctionBegin;
  auto simple = mField.getInterface<Simple>();
 
  enum bases { AINSWORTH, DEMKOWICZ, LASBASETOPT };
  const char *list_bases[LASBASETOPT] = {"ainsworth", "demkowicz"};
  PetscInt choice_base_value = AINSWORTH;
  CHKERR PetscOptionsGetEList(PETSC_NULLPTR, PETSC_NULLPTR, "-base", list_bases,
                              LASBASETOPT, &choice_base_value, PETSC_NULLPTR);
 
  FieldApproximationBase base;
  switch (choice_base_value) {
  case AINSWORTH:
    base = AINSWORTH_LEGENDRE_BASE;
    break;
  case DEMKOWICZ:
    base = DEMKOWICZ_JACOBI_BASE;
    break;
  default:
    base = LASTBASE;
    break;
  }
 
  int order = 3;
  CHKERR PetscOptionsGetInt(PETSC_NULLPTR, "", "-order", &order,
                            PETSC_NULLPTR);
 
  if (SPACE_DIM == 3) {
    Range domain_ents;
    CHKERR mField.get_moab().get_entities_by_handle(simple->getMeshset(),
                                                    domain_ents, true);
    auto tet_hex_ents = ElasticTie::getPhysicalVolumeEntities(domain_ents);
 
    if (!tet_hex_ents.empty()) {
      EntityHandle domain_meshset;
      CHKERR ElasticTie::makeClosureMeshset(
          mField.get_moab(), tet_hex_ents, SPACE_DIM, domain_meshset);
      simple->getMeshset() = domain_meshset;
 
      Range boundary_faces;
      if (mField.get_comm_size() > 1) {
        Skinner skin(&mField.get_moab());
        Range local_skin;
        CHKERR skin.find_skin(0, tet_hex_ents, false, local_skin);
        boundary_faces = intersect(physicalBoundaryFaces, local_skin);
      } else if (!physicalBoundaryFaces.empty()) {
        boundary_faces = physicalBoundaryFaces;
      } else {
        Skinner skin(&mField.get_moab());
        CHKERR skin.find_skin(0, tet_hex_ents, false, boundary_faces);
      }
      EntityHandle boundary_meshset;
      CHKERR ElasticTie::makeClosureMeshset(mField.get_moab(), boundary_faces,
                                            SPACE_DIM - 1, boundary_meshset);
      simple->getBoundaryMeshSet() = boundary_meshset;
    }
  }
 
  CHKERR ElasticExample::setupProblem();
  CHKERR ElasticTie::initializeTieConstraints(mField, tieData, "U", base,
                                              order);
 
  MoFEMFunctionReturn(0);
}

setUpSolver()

MoFEMErrorCode ElasticTieMeshExample::setUpSolver ( SmartPetscObj< KSP >  solver )

inlineoverrideprotectedvirtual

[Push operators to pipeline]

Reimplemented from ElasticExample.

Definition at line 1341 of file ElasticTie.hpp.

                                                            {
  MoFEMFunctionBegin;
  CHKERR ElasticTie::setupTieSolver(mField, tieData, solver);
  MoFEMFunctionReturn(0);
}

Member Data Documentation

physicalBoundaryFaces

Range ElasticTieMeshExample::physicalBoundaryFaces

protected

Definition at line 1260 of file ElasticTie.hpp.

tieData

boost::shared_ptr<ElasticTie::CommonData> ElasticTieMeshExample::tieData = nullptr

protected

Definition at line 1259 of file ElasticTie.hpp.

---

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

• tutorials/vec-11_elastic_tie_mesh/src/ElasticTie.hpp