eigen_elastic.cpp File Reference

#include <MoFEM.hpp>
#include <slepceps.h>

Go to the source code of this file.

Classes
struct	ElementsAndOps< DIM >
	[Operators_definition] More...
struct	ElementsAndOps< 2 >
struct	ElementsAndOps< 3 >
struct	Example
	[Example] More...

Typedefs
using	DomainEle = ElementsAndOps< SPACE_DIM >::DomainEle
using	DomainEleOp = DomainEle::UserDataOperator
using	PostProcEle = PostProcBrokenMeshInMoab< DomainEle >
using	OpDomainGradGrad = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpGradSymTensorGrad< 1, SPACE_DIM, SPACE_DIM, 0 >
using	OpDomainMass = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpMass< 1, SPACE_DIM >

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

Variables
constexpr int	SPACE_DIM
double	rho = 7829e-9
	[Operators_definition]
double	young_modulus = 2.07e8
double	poisson_ratio = 0.33
double	bulk_modulus_K
double	shear_modulus_G
int	order = 2
	[Physical_parameters]
static char	help [] = "...\n\n"
	[Check]

Typedef Documentation

DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 29 of file eigen_elastic.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 30 of file eigen_elastic.cpp.

OpDomainGradGrad

using OpDomainGradGrad = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm< GAUSS>::OpGradSymTensorGrad<1, SPACE_DIM, SPACE_DIM, 0>

Definition at line 33 of file eigen_elastic.cpp.

OpDomainMass

using OpDomainMass = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::BiLinearForm<GAUSS>::OpMass<1, SPACE_DIM>

Definition at line 36 of file eigen_elastic.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 31 of file eigen_elastic.cpp.

Function Documentation

main()

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

[Main]

[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 insterface

[Create MoFEM]

[Example]

[Example]

Definition at line 474 of file eigen_elastic.cpp.

                                 {
 
  // Initialisation of MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  SlepcInitialize(&argc, &argv, param_file, help);
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  // Add logging channel for example
  auto core_log = logging::core::get();
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "EXAMPLE"));
  LogManager::setLog("EXAMPLE");
  MOFEM_LOG_TAG("EXAMPLE", "example");
 
  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 insterface
    //! [Create MoFEM]
 
    //! [Example]
    Example ex(m_field);
    CHKERR ex.runProblem();
    //! [Example]
  }
  CATCH_ERRORS;
 
  SlepcFinalize();
  CHKERR MoFEM::Core::Finalize();
}

Variable Documentation

bulk_modulus_K

double bulk_modulus_K

Definition at line 45 of file eigen_elastic.cpp.

help

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

static

[Check]

Definition at line 471 of file eigen_elastic.cpp.

order

int order = 2

[Physical_parameters]

Definition at line 49 of file eigen_elastic.cpp.

poisson_ratio

double poisson_ratio = 0.33

Definition at line 43 of file eigen_elastic.cpp.

rho

double rho = 7829e-9

[Operators_definition]

[Physical_parameters]

Definition at line 41 of file eigen_elastic.cpp.

shear_modulus_G

double shear_modulus_G

Definition at line 46 of file eigen_elastic.cpp.

SPACE_DIM

constexpr int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

Definition at line 26 of file eigen_elastic.cpp.

young_modulus

double young_modulus = 2.07e8

Definition at line 42 of file eigen_elastic.cpp.