#include <MoFEM.hpp>
#include <MatrixFunction.hpp>
#include <ContactNaturalBC.hpp>
#include <ADOLCPlasticity.hpp>
#include <ADOLCPlasticityMaterialModels.hpp>

Classes
struct	ElementsAndOps< DIM >

struct	ElementsAndOps< 2 >

struct	ElementsAndOps< 3 >

struct	PlasticProblem

struct	Monitor
	[Push operators to pipeline] More...

Typedefs
using	EntData = EntitiesFieldData::EntData

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

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

using	DomainEleOp = DomainEle::UserDataOperator

using	BoundaryEleOp = BoundaryEle::UserDataOperator

using	DomainNaturalBC = NaturalBC< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >

using	OpBodyForce = DomainNaturalBC::OpFlux< NaturalMeshsetType< BLOCKSET >, 1, SPACE_DIM >

using	BoundaryRhsBCs = NaturalBC< BoundaryEleOp >::Assembly< PETSC >::LinearForm< GAUSS >

using	OpBoundaryRhsBCs = BoundaryRhsBCs::OpFlux< ContactOps::BoundaryBCs, 1, SPACE_DIM >

using	BoundaryLhsBCs = NaturalBC< BoundaryEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >

using	OpBoundaryLhsBCs = BoundaryLhsBCs::OpFlux< ContactOps::BoundaryBCs, 1, SPACE_DIM >

using	PostProcEleDomain = ElementsAndOps< SPACE_DIM >::PostProcEleDomain

using	SideEle = ElementsAndOps< SPACE_DIM >::SideEle

using	PostProcEleBdy = ElementsAndOps< SPACE_DIM >::PostProcEleBdy

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

Variables
constexpr int	SPACE_DIM

constexpr FieldSpace	CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE

double	scale = 1.

static boost::shared_ptr< TSUpdate >	ts_update_ptr = nullptr

static char	help [] = "...\n\n"
	[Check] More...

Typedef Documentation

◆ BoundaryEle

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

Select finite element type for integrate on boundary based on problem dimension

Definition at line 28 of file adolc_plasticity.cpp.

◆ BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Operators used to assemble boundary integrals

Definition at line 35 of file adolc_plasticity.cpp.

◆ BoundaryLhsBCs

using BoundaryLhsBCs = NaturalBC<BoundaryEleOp>::Assembly<PETSC>::BiLinearForm<GAUSS>

Use Gauss quadrature rule and PETSc assembly to integrate neural boundary conditions. Select bi-linear forms operators.

Definition at line 70 of file adolc_plasticity.cpp.

◆ BoundaryRhsBCs

using BoundaryRhsBCs = NaturalBC<BoundaryEleOp>::Assembly<PETSC>::LinearForm<GAUSS>

Use Gauss quadrature rule and PETSc assembly to integrate neural boundary conditions. Select linear forms operators.

Definition at line 58 of file adolc_plasticity.cpp.

◆ DomainEle

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

Select finite element type for integration on domain based on problem dimension

Examples: adolc_plasticity.cpp.

Definition at line 22 of file adolc_plasticity.cpp.

◆ DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Operators used to assemble domain integrals

Definition at line 31 of file adolc_plasticity.cpp.

◆ DomainNaturalBC

using DomainNaturalBC = NaturalBC<DomainEleOp>::Assembly<PETSC>::LinearForm<GAUSS>

Linear forms used to integrate body forces

Definition at line 40 of file adolc_plasticity.cpp.

◆ EntData

using EntData = EntitiesFieldData::EntData

Examples: adolc_plasticity.cpp.

Definition at line 17 of file adolc_plasticity.cpp.

◆ OpBodyForce

using OpBodyForce = DomainNaturalBC::OpFlux<NaturalMeshsetType<BLOCKSET>, 1, SPACE_DIM>

Select linear froms reading data from blockest (e.g. "BODY_FORCE") and applying body force.

Definition at line 46 of file adolc_plasticity.cpp.

◆ OpBoundaryLhsBCs

using OpBoundaryLhsBCs = BoundaryLhsBCs::OpFlux<ContactOps::BoundaryBCs, 1, SPACE_DIM>

Use specialization from adv-1 integrating boundary conditions on forces and with springs

Definition at line 76 of file adolc_plasticity.cpp.

◆ OpBoundaryRhsBCs

using OpBoundaryRhsBCs = BoundaryRhsBCs::OpFlux<ContactOps::BoundaryBCs, 1, SPACE_DIM>

Definition at line 64 of file adolc_plasticity.cpp.

◆ PostProcEleBdy

using PostProcEleBdy = ElementsAndOps<SPACE_DIM>::PostProcEleBdy

Definition at line 104 of file adolc_plasticity.cpp.

◆ PostProcEleDomain

using PostProcEleDomain = ElementsAndOps<SPACE_DIM>::PostProcEleDomain

Definition at line 102 of file adolc_plasticity.cpp.

◆ SideEle

using SideEle = ElementsAndOps<SPACE_DIM>::SideEle

Definition at line 103 of file adolc_plasticity.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]

[PlasticProblem]

Examples: adolc_plasticity.cpp.

Definition at line 1189 of file adolc_plasticity.cpp.

                                  {
  
 #ifdef ADD_CONTACT
 #ifdef PYTHON_SDF
   Py_Initialize();
   np::initialize();
 #endif
 #endif // ADD_CONTACT
  
   // 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(), "PlasticPrb"));
   LogManager::setLog("PlasticPrb");
   MOFEM_LOG_TAG("PlasticPrb", "PlasticPrb");
   MOFEM_LOG("PlasticPrb", Sev::inform) << "SPACE_DIM " << SPACE_DIM;
 #ifdef ADD_CONTACT
   core_log->add_sink(
       LogManager::createSink(LogManager::getStrmWorld(), "CONTACT"));
   LogManager::setLog("CONTACT");
   MOFEM_LOG_TAG("CONTACT", "Contact");
 #endif // ADD_CONTACT
  
   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]
  
     //! [PlasticProblem]
     PlasticProblem ex(m_field);
     CHKERR ex.runProblem();
     //! [PlasticProblem]
   }
   CATCH_ERRORS;
  
   CHKERR MoFEM::Core::Finalize();
 #ifdef ADD_CONTACT
 #ifdef PYTHON_SDF
   if (Py_FinalizeEx() < 0) {
     exit(120);
   }
 #endif
 #endif // ADD_CONTACT
   return 0;
 }

Variable Documentation

◆ CONTACT_SPACE

constexpr FieldSpace CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE

constexpr

Examples: adolc_plasticity.cpp.

Definition at line 100 of file adolc_plasticity.cpp.

◆ help

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

static

[Check]

Examples: adolc_plasticity.cpp.

Definition at line 1187 of file adolc_plasticity.cpp.

◆ scale

double scale = 1.

Examples: adolc_plasticity.cpp.

Definition at line 109 of file adolc_plasticity.cpp.

◆ SPACE_DIM

constexpr int SPACE_DIM