cornea_elastic.cpp File Reference

#include <MoFEM.hpp>
#include <MatrixFunction.hpp>
#include <PostProcOnRefMesh.hpp>
#include <HenckyOps.hpp>
#include <OpPostProcElastic.hpp>

Go to the source code of this file.

Classes
struct	ElementsAndOps< DIM >
struct	ElementsAndOps< 3 >
struct	Example
	[Example] More...
struct	Monitor
	[Push operators to pipeline] More...

Typedefs
using	EntData = DataForcesAndSourcesCore::EntData
using	DomainEle = ElementsAndOps< SPACE_DIM >::DomainEle
using	DomainEleOp = ElementsAndOps< SPACE_DIM >::DomainEleOp
using	PostProcEle = ElementsAndOps< SPACE_DIM >::PostProcEle
using	OpK = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpGradTensorGrad< 1, SPACE_DIM, SPACE_DIM, 1 >
using	OpMass = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpMass< 1, SPACE_DIM >
using	OpInternalForce = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpGradTimesTensor< 1, SPACE_DIM, SPACE_DIM >
using	OpBodyForce = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpSource< 1, SPACE_DIM >
using	OpInertiaForce = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpBaseTimesVector< 1, SPACE_DIM, 1 >

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

Variables
constexpr int	SPACE_DIM = 3
constexpr bool	is_quasi_static = false
constexpr double	rho = 1
constexpr double	omega = 2.4
constexpr double	young_modulus = 1
constexpr double	poisson_ratio = 0.25
constexpr double	bulk_modulus_K = young_modulus / (3 * (1 - 2 * poisson_ratio))
constexpr double	shear_modulus_G = young_modulus / (2 * (1 + poisson_ratio))
static double *	ts_time_ptr
static double *	ts_aa_ptr
static char	help [] = "...\n\n"
	[Check] More...

Typedef Documentation

DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 40 of file cornea_elastic.cpp.

DomainEleOp

using DomainEleOp = ElementsAndOps<SPACE_DIM>::DomainEleOp

Definition at line 41 of file cornea_elastic.cpp.

EntData

using EntData = DataForcesAndSourcesCore::EntData

Definition at line 39 of file cornea_elastic.cpp.

OpBodyForce

using OpBodyForce = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::LinearForm< GAUSS>::OpSource<1, SPACE_DIM>

Definition at line 50 of file cornea_elastic.cpp.

OpInertiaForce

using OpInertiaForce = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<GAUSS>::OpBaseTimesVector<1, SPACE_DIM, 1>

Definition at line 52 of file cornea_elastic.cpp.

OpInternalForce

using OpInternalForce = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<GAUSS>::OpGradTimesTensor<1, SPACE_DIM, SPACE_DIM>

Definition at line 48 of file cornea_elastic.cpp.

OpK

using OpK = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm< GAUSS>::OpGradTensorGrad<1, SPACE_DIM, SPACE_DIM, 1>

Examples

PoissonOperators.hpp, and simple_elasticity.cpp.

Definition at line 44 of file cornea_elastic.cpp.

OpMass

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

Examples

scalar_check_approximation.cpp.

Definition at line 46 of file cornea_elastic.cpp.

PostProcEle

using PostProcEle = ElementsAndOps<SPACE_DIM>::PostProcEle

Definition at line 42 of file cornea_elastic.cpp.

Function Documentation

main()

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

Definition at line 449 of file cornea_elastic.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 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;
 
  CHKERR MoFEM::Core::Finalize();
}

Variable Documentation

bulk_modulus_K

constexpr double bulk_modulus_K = young_modulus / (3 * (1 - 2 * poisson_ratio))

constexpr

Examples

contact.cpp, plastic.cpp, and thermo_plastic.cpp.

Definition at line 60 of file cornea_elastic.cpp.

help

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

static

[Check]

Definition at line 447 of file cornea_elastic.cpp.

is_quasi_static

constexpr bool is_quasi_static = false

constexpr

Examples

NonlinearElasticElementInterface.hpp.

Definition at line 55 of file cornea_elastic.cpp.

omega

constexpr double omega = 2.4

constexpr

Definition at line 57 of file cornea_elastic.cpp.

poisson_ratio

constexpr double poisson_ratio = 0.25

constexpr

Examples

OpPostProcElastic.hpp, and Remodeling.cpp.

Definition at line 59 of file cornea_elastic.cpp.

rho

constexpr double rho = 1

constexpr

Examples

HookeElement.cpp, Remodeling.cpp, Remodeling.hpp, and heat_method.cpp.

Definition at line 56 of file cornea_elastic.cpp.

shear_modulus_G

constexpr double shear_modulus_G = young_modulus / (2 * (1 + poisson_ratio))

constexpr

Examples

contact.cpp, plastic.cpp, and thermo_plastic.cpp.

Definition at line 61 of file cornea_elastic.cpp.

SPACE_DIM

constexpr int SPACE_DIM = 3

constexpr

Definition at line 37 of file cornea_elastic.cpp.

ts_aa_ptr

double* ts_aa_ptr

static

Definition at line 69 of file cornea_elastic.cpp.

ts_time_ptr

double* ts_time_ptr

static

Definition at line 68 of file cornea_elastic.cpp.

young_modulus

constexpr double young_modulus = 1

constexpr

Examples

Remodeling.cpp.

Definition at line 58 of file cornea_elastic.cpp.