incompressible_elasticity.cpp File Reference

Incompressible elasticity problem. More...

#include <MoFEM.hpp>

Go to the source code of this file.

Classes
struct	MonitorIncompressible
struct	DomainEleOpStab
	Element used to specialise assembly. More...
struct	Incompressible
struct	OpCalculateLameStress< DIM >
struct	SetUpSchur
	[Push operators to pipeline] More...
struct	SetUpSchurImpl

Typedefs
using	EntData = EntitiesFieldData::EntData
using	DomainEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >::DomainEle
using	DomainEleOp = DomainEle::UserDataOperator
using	BoundaryEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >::BoundaryEle
using	BoundaryEleOp = BoundaryEle::UserDataOperator
using	PostProcEle = PostProcBrokenMeshInMoab< DomainEle >
using	SkinPostProcEle = PostProcBrokenMeshInMoab< BoundaryEle >
using	SetPtsData = FieldEvaluatorInterface::SetPtsData

Functions
int	main (int argc, char *argv[])

Variables
constexpr int	SPACE_DIM
constexpr AssemblyType	AT
constexpr IntegrationType	IT
constexpr CoordinateTypes	coord_type = CARTESIAN
PetscBool	doEvalField
int	order = 2
	[Specialisation for assembly]
int	geom_order = 1
static double	young_modulus = 100
static double	poisson_ratio = 0.25
static double	mu
static double	lambda
PetscBool	isDiscontinuousPressure = PETSC_FALSE
static char	help [] = "...\n\n"
	[Check]

Detailed Description

Incompressible elasticity problem.

Definition in file incompressible_elasticity.cpp.

Typedef Documentation

BoundaryEle

using BoundaryEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::BoundaryEle

Definition at line 24 of file incompressible_elasticity.cpp.

BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 26 of file incompressible_elasticity.cpp.

DomainEle

using DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle

Definition at line 22 of file incompressible_elasticity.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 23 of file incompressible_elasticity.cpp.

EntData

using EntData = EntitiesFieldData::EntData

Definition at line 21 of file incompressible_elasticity.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 27 of file incompressible_elasticity.cpp.

SetPtsData

using SetPtsData = FieldEvaluatorInterface::SetPtsData

Definition at line 29 of file incompressible_elasticity.cpp.

SkinPostProcEle

using SkinPostProcEle = PostProcBrokenMeshInMoab<BoundaryEle>

Definition at line 28 of file incompressible_elasticity.cpp.

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

[Register MoFEM discrete manager in PETSc]

[Register MoFEM discrete manager in PETSc

[Create MoAB]

< mesh database

< mesh database interface

[Create MoAB]

[Create MoFEM]

< finite element database

< finite element database interface

[Create MoFEM]

[Load mesh]

[Load mesh]

[CONTACT]

[CONTACT]

Definition at line 940 of file incompressible_elasticity.cpp.

                                 {
 
  // Initialisation of MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  // Add logging channel for CONTACT
  auto core_log = logging::core::get();
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "INCOMP_ELASTICITY"));
  LogManager::setLog("INCOMP_ELASTICITY");
  MOFEM_LOG_TAG("INCOMP_ELASTICITY", "Indent");
 
  try {
 
    //! [Register MoFEM discrete manager in PETSc]
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
    //! [Register MoFEM discrete manager in PETSc
 
    //! [Create MoAB]
    moab::Core mb_instance;              ///< mesh database
    moab::Interface &moab = mb_instance; ///< mesh database interface
    //! [Create MoAB]
 
    //! [Create MoFEM]
    MoFEM::Core core(moab);           ///< finite element database
    MoFEM::Interface &m_field = core; ///< finite element database interface
    //! [Create MoFEM]
 
    //! [Load mesh]
    Simple *simple = m_field.getInterface<Simple>();
    CHKERR simple->getOptions();
    CHKERR simple->loadFile("");
    //! [Load mesh]
 
    //! [CONTACT]
    Incompressible ex(m_field);
    CHKERR ex.runProblem();
    //! [CONTACT]
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
 
  return 0;
}

Variable Documentation

AT

constexpr AssemblyType AT

constexpr

Initial value:

=
    (SCHUR_ASSEMBLE) ? AssemblyType::SCHUR
                     : AssemblyType::PETSC

Definition at line 13 of file incompressible_elasticity.cpp.

coord_type

constexpr CoordinateTypes coord_type = CARTESIAN

constexpr

Examples

free_surface.cpp.

Definition at line 19 of file incompressible_elasticity.cpp.

doEvalField

PetscBool doEvalField

Examples

seepage.cpp, and thermo_elastic.cpp.

Definition at line 32 of file incompressible_elasticity.cpp.

geom_order

int geom_order = 1

Definition at line 177 of file incompressible_elasticity.cpp.

help

char help[] = "...\n\n"

static

[Check]

Definition at line 938 of file incompressible_elasticity.cpp.

isDiscontinuousPressure

PetscBool isDiscontinuousPressure = PETSC_FALSE

Definition at line 183 of file incompressible_elasticity.cpp.

IT

constexpr IntegrationType IT

constexpr

Initial value:

=
    IntegrationType::GAUSS

Definition at line 17 of file incompressible_elasticity.cpp.

lambda

double lambda

inlinestatic

Examples

ElasticityMixedFormulation.hpp, bernstein_bezier_generate_base.cpp, free_surface.cpp, plate.cpp, and shallow_wave.cpp.

Definition at line 181 of file incompressible_elasticity.cpp.

mu

double mu

inlinestatic

Definition at line 180 of file incompressible_elasticity.cpp.

order

int order = 2

[Specialisation for assembly]

Definition at line 176 of file incompressible_elasticity.cpp.

poisson_ratio

double poisson_ratio = 0.25

inlinestatic

Definition at line 179 of file incompressible_elasticity.cpp.

SPACE_DIM

constexpr int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

Definition at line 10 of file incompressible_elasticity.cpp.

young_modulus

double young_modulus = 100

inlinestatic

Definition at line 178 of file incompressible_elasticity.cpp.