adolc_plasticity.cpp File Reference

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

Go to the source code of this file.

Classes
struct	ElementsAndOps< 2 >
struct	ElementsAndOps< 3 >
struct	PlasticProblem
struct	Monitor
	[Push operators to pipeline] More...

Typedefs
using	DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle
using	BoundaryEle
using	DomainNaturalBC
using	OpBodyForce
using	BoundaryRhsBCs
using	OpBoundaryRhsBCs
using	BoundaryLhsBCs
using	OpBoundaryLhsBCs
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]

Typedef Documentation

BoundaryEle

using BoundaryEle

Initial value:

 
    PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::BoundaryEle

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

Definition at line 27 of file adolc_plasticity.cpp.

BoundaryLhsBCs

using BoundaryLhsBCs

Initial value:

 
    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 69 of file adolc_plasticity.cpp.

BoundaryRhsBCs

using BoundaryRhsBCs

Initial value:

 
    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 57 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.

DomainNaturalBC

using DomainNaturalBC

Initial value:

 
    NaturalBC<DomainEleOp>::Assembly<PETSC>::LinearForm<GAUSS>

Linear forms used to integrate body forces

Definition at line 39 of file adolc_plasticity.cpp.

OpBodyForce

using OpBodyForce

Initial value:

 
    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 45 of file adolc_plasticity.cpp.

OpBoundaryLhsBCs

using OpBoundaryLhsBCs

Initial value:

 
    BoundaryLhsBCs::OpFlux<ContactOps::BoundaryBCs, 1, SPACE_DIM>

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

Definition at line 75 of file adolc_plasticity.cpp.

OpBoundaryRhsBCs

using OpBoundaryRhsBCs

Initial value:

 
    BoundaryRhsBCs::OpFlux<ContactOps::BoundaryBCs, 1, SPACE_DIM>

Definition at line 63 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]

[PlasticProblem]

Definition at line 1342 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

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]

Definition at line 1340 of file adolc_plasticity.cpp.

scale

double scale = 1.

Examples

adolc_plasticity.cpp.

Definition at line 106 of file adolc_plasticity.cpp.

SPACE_DIM

int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

Examples

adolc_plasticity.cpp, dynamic_first_order_con_law.cpp, and free_surface.cpp.

Definition at line 14 of file adolc_plasticity.cpp.

ts_update_ptr

boost::shared_ptr<TSUpdate> ts_update_ptr = nullptr

static

Examples

adolc_plasticity.cpp.

Definition at line 792 of file adolc_plasticity.cpp.