contact.cpp File Reference

#include <MoFEM.hpp>
#include <MatrixFunction.hpp>
#include <GenericElementInterface.hpp>
#include <HenckyOps.hpp>
#include <ContactOps.hpp>
#include <PostProcContact.hpp>
#include <ContactNaturalBC.hpp>

Go to the source code of this file.

Classes
struct	ElementsAndOps< 2 >
struct	ElementsAndOps< 3 >
struct	Contact
struct	Contact::ScaledTimeScale
struct	SetUpSchur
	[Push operators to pipeline] More...
struct	SetUpSchurImpl

Namespaces
namespace	ContactOps

Macros
#define	EXECUTABLE_DIMENSION   3
#define	SCHUR_ASSEMBLE   0

Typedefs
using	DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle
using	BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle
using	SideEle = ElementsAndOps<SPACE_DIM>::FaceSideEle
using	OpSpringLhs
	[Operators used for contact]
using	OpSpringRhs
using	DomainRhsBCs = NaturalBC<DomainEleOp>::Assembly<AT>::LinearForm<IT>
using	OpDomainRhsBCs
using	BoundaryRhsBCs = NaturalBC<BoundaryEleOp>::Assembly<AT>::LinearForm<IT>
using	OpBoundaryRhsBCs
using	BoundaryLhsBCs = NaturalBC<BoundaryEleOp>::Assembly<AT>::BiLinearForm<IT>
using	OpBoundaryLhsBCs

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

Variables
constexpr AssemblyType	AT
constexpr IntegrationType	IT
constexpr int	SPACE_DIM
constexpr FieldSpace	CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE
	[Specialisation for assembly]
PetscBool	is_quasi_static = PETSC_TRUE
	[Operators used for contact]
int	order = 2
int	contact_order = 2
int	sigma_order = 1
int	geom_order = 1
double	young_modulus = 100
double	poisson_ratio = 0.25
double	rho = 0.0
double	spring_stiffness = 0.0
double	vis_spring_stiffness = 0.0
double	alpha_damping = 0
double	scale = 1.
PetscBool	is_axisymmetric = PETSC_FALSE
int	atom_test = 0
double	ContactOps::cn_contact = 0.1
static char	help [] = "...\n\n"
	[Check]

Detailed Description

\CONTACT contact.cpp

CONTACT of contact problem

Copyright

Copyright (c) 2023

Definition in file contact.cpp.

Macro Definition Documentation

EXECUTABLE_DIMENSION

#define EXECUTABLE_DIMENSION   3

Definition at line 14 of file contact.cpp.

SCHUR_ASSEMBLE

#define SCHUR_ASSEMBLE   0

Examples

plastic.cpp, and test_broken_space.cpp.

Definition at line 18 of file contact.cpp.

Typedef Documentation

BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 64 of file contact.cpp.

BoundaryLhsBCs

using BoundaryLhsBCs = NaturalBC<BoundaryEleOp>::Assembly<AT>::BiLinearForm<IT>

Definition at line 119 of file contact.cpp.

BoundaryRhsBCs

using BoundaryRhsBCs = NaturalBC<BoundaryEleOp>::Assembly<AT>::LinearForm<IT>

Definition at line 116 of file contact.cpp.

DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 62 of file contact.cpp.

DomainRhsBCs

using DomainRhsBCs = NaturalBC<DomainEleOp>::Assembly<AT>::LinearForm<IT>

Definition at line 113 of file contact.cpp.

OpBoundaryLhsBCs

using OpBoundaryLhsBCs

Initial value:

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

Definition at line 120 of file contact.cpp.

OpBoundaryRhsBCs

using OpBoundaryRhsBCs

Initial value:

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

Definition at line 117 of file contact.cpp.

OpDomainRhsBCs

using OpDomainRhsBCs

Initial value:

 
    DomainRhsBCs::OpFlux<ContactOps::DomainBCs, 1, SPACE_DIM>

Definition at line 114 of file contact.cpp.

OpSpringLhs

using OpSpringLhs

Initial value:

 FormsIntegrators<BoundaryEleOp>::Assembly<AT>::BiLinearForm<
    IT>::OpMass<1, SPACE_DIM>

[Operators used for contact]

Definition at line 72 of file contact.cpp.

OpSpringRhs

using OpSpringRhs

Initial value:

 FormsIntegrators<BoundaryEleOp>::Assembly<AT>::LinearForm<
    IT>::OpBaseTimesVector<1, SPACE_DIM, 1>

Definition at line 74 of file contact.cpp.

SideEle

using SideEle = ElementsAndOps<SPACE_DIM>::FaceSideEle

Definition at line 66 of file contact.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 926 of file contact.cpp.

                                 {
 
#ifdef PYTHON_SDF
  Py_Initialize();
  np::initialize();
#endif
 
  // 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(), "CONTACT"));
  LogManager::setLog("CONTACT");
  MOFEM_LOG_TAG("CONTACT", "Indent");
 
  try {
 
    //! [Register MoFEM discrete manager in PETSc]
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
    DMType dm_name_mg = "DMMOFEM_MG";
    CHKERR DMRegister_MGViaApproxOrders(dm_name_mg);
    //! [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]
    Contact ex(m_field);
    CHKERR ex.runProblem();
    //! [CONTACT]
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
 
#ifdef PYTHON_SDF
  if (Py_FinalizeEx() < 0) {
    exit(120);
  }
#endif
 
  return 0;
}

Variable Documentation

alpha_damping

double alpha_damping = 0

Definition at line 89 of file contact.cpp.

AT

AssemblyType AT

constexpr

Initial value:

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

Definition at line 36 of file contact.cpp.

atom_test

int atom_test = 0

Examples

thermo_elastic.cpp.

Definition at line 97 of file contact.cpp.

contact_order

int contact_order = 2

Definition at line 81 of file contact.cpp.

CONTACT_SPACE

FieldSpace CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE

constexpr

[Specialisation for assembly]

Definition at line 69 of file contact.cpp.

geom_order

int geom_order = 1

Definition at line 83 of file contact.cpp.

help

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

static

[Check]

Definition at line 924 of file contact.cpp.

is_axisymmetric

PetscBool is_axisymmetric = PETSC_FALSE

Examples

ContactOps.hpp.

Definition at line 93 of file contact.cpp.

is_quasi_static

PetscBool is_quasi_static = PETSC_TRUE

[Operators used for contact]

Definition at line 78 of file contact.cpp.

IT

IntegrationType IT

constexpr

Initial value:

=
    IntegrationType::GAUSS

Definition at line 39 of file contact.cpp.

order

int order = 2

Definition at line 80 of file contact.cpp.

poisson_ratio

double poisson_ratio = 0.25

Definition at line 85 of file contact.cpp.

rho

double rho = 0.0

Definition at line 86 of file contact.cpp.

scale

double scale = 1.

Definition at line 91 of file contact.cpp.

sigma_order

int sigma_order = 1

Definition at line 82 of file contact.cpp.

SPACE_DIM

int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

Definition at line 55 of file contact.cpp.

spring_stiffness

double spring_stiffness = 0.0

Definition at line 87 of file contact.cpp.

vis_spring_stiffness

double vis_spring_stiffness = 0.0

Definition at line 88 of file contact.cpp.

young_modulus

double young_modulus = 100

Definition at line 84 of file contact.cpp.