#include <MoFEM.hpp>
#include <MatrixFunction.hpp>
#include <IntegrationRules.hpp>
#include <HenckyOps.hpp>
#include <PlasticOps.hpp>
#include <PlasticNaturalBCs.hpp>

Classes
struct	ElementsAndOps< DIM >

struct	ElementsAndOps< 2 >

struct	ElementsAndOps< 3 >

struct	Example
	[Example] More...

struct	Example::ScaledTimeScale

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

struct	SetUpSchurImpl

Macros
#define	EXECUTABLE_DIMENSION 3

Typedefs
using	EntData = EntitiesFieldData::EntData

using	DomainEle = ElementsAndOps< SPACE_DIM >::DomainEle

using	DomainEleOp = DomainEle::UserDataOperator

using	BoundaryEle = ElementsAndOps< SPACE_DIM >::BoundaryEle

using	BoundaryEleOp = BoundaryEle::UserDataOperator

using	PostProcEle = PostProcBrokenMeshInMoab< DomainEle >

using	SkinPostProcEle = PostProcBrokenMeshInMoab< BoundaryEle >

using	SideEle = ElementsAndOps< SPACE_DIM >::SideEle

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

using	OpDomainRhsBCs = DomainRhsBCs::OpFlux< PlasticOps::DomainBCs, 1, SPACE_DIM >

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

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

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

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

Functions
double	iso_hardening_exp (double tau, double b_iso)

double	iso_hardening (double tau, double H, double Qinf, double b_iso, double sigmaY)

double	iso_hardening_dtau (double tau, double H, double Qinf, double b_iso)

template<typename T , int DIM>
auto	kinematic_hardening (FTensor::Tensor2_symmetric< T, DIM > &t_plastic_strain, double C1_k)

template<int DIM>
auto	kinematic_hardening_dplastic_strain (double C1_k)

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

Variables
constexpr int	SPACE_DIM

constexpr auto	size_symm = (SPACE_DIM * (SPACE_DIM + 1)) / 2

constexpr AssemblyType	AT

constexpr IntegrationType	IT

constexpr FieldSpace	CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE

PetscBool	is_large_strains = PETSC_TRUE
	Large strains. More...

PetscBool	set_timer = PETSC_FALSE
	Set timer. More...

double	scale = 1.

double	young_modulus = 206913
	Young modulus. More...

double	poisson_ratio = 0.29
	Poisson ratio. More...

double	sigmaY = 450
	Yield stress. More...

double	H = 129
	Hardening. More...

double	visH = 0
	Viscous hardening. More...

double	zeta = 5e-2
	Viscous hardening. More...

double	Qinf = 265
	Saturation yield stress. More...

double	b_iso = 16.93
	Saturation exponent. More...

double	C1_k = 0
	Kinematic hardening. More...

double	cn0 = 1

double	cn1 = 1

int	order = 2
	Order displacement. More...

int	tau_order = order - 2
	Order of tau files. More...

int	ep_order = order - 1
	Order of ep files. More...

int	geom_order = 2
	Order if fixed. More...

PetscBool	is_quasi_static = PETSC_TRUE

double	rho = 0.0

double	alpha_damping = 0

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

Macro Definition Documentation

◆ EXECUTABLE_DIMENSION

#define EXECUTABLE_DIMENSION 3

Examples: adolc_plasticity.cpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, level_set.cpp, plastic.cpp, plot_base.cpp, poisson_2d_homogeneous.cpp, scalar_check_approximation.cpp, tensor_divergence_operator.cpp, and test_broken_space.cpp.

Definition at line 13 of file plastic.cpp.

Typedef Documentation

◆ BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 57 of file plastic.cpp.

◆ BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 58 of file plastic.cpp.

◆ BoundaryLhsBCs

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

Definition at line 171 of file plastic.cpp.

◆ BoundaryRhsBCs

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

Definition at line 168 of file plastic.cpp.

◆ DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 55 of file plastic.cpp.

◆ DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 56 of file plastic.cpp.

◆ DomainRhsBCs

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

Definition at line 165 of file plastic.cpp.

◆ EntData

using EntData = EntitiesFieldData::EntData

Examples: plastic.cpp.

Definition at line 54 of file plastic.cpp.

◆ OpBoundaryLhsBCs

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

Definition at line 173 of file plastic.cpp.

◆ OpBoundaryRhsBCs

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

Definition at line 170 of file plastic.cpp.

◆ OpDomainRhsBCs

using OpDomainRhsBCs = DomainRhsBCs::OpFlux<PlasticOps::DomainBCs, 1, SPACE_DIM>

Definition at line 167 of file plastic.cpp.

◆ PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 59 of file plastic.cpp.

◆ SideEle

using SideEle = ElementsAndOps<SPACE_DIM>::SideEle

Definition at line 61 of file plastic.cpp.

◆ SkinPostProcEle

using SkinPostProcEle = PostProcBrokenMeshInMoab<BoundaryEle>

Definition at line 60 of file plastic.cpp.

Function Documentation

◆ iso_hardening()

double iso_hardening	(	double	tau,
		double	H,
		double	Qinf,
		double	b_iso,
		double	sigmaY
	)

inline

Isotropic hardening

Examples: plastic.cpp, and PlasticOpsGeneric.hpp.

Definition at line 72 of file plastic.cpp.

                                            {
   return H * tau + Qinf * (1. - iso_hardening_exp(tau, b_iso)) + sigmaY;
 }

◆ iso_hardening_dtau()

double iso_hardening_dtau	(	double	tau,
		double	H,
		double	Qinf,
		double	b_iso
	)

inline

Examples: plastic.cpp, and PlasticOpsGeneric.hpp.

Definition at line 77 of file plastic.cpp.

                                                {
   auto r = [&](auto tau) {
     return H + Qinf * b_iso * iso_hardening_exp(tau, b_iso);
   };
   constexpr double eps = 1e-12;
   return std::max(r(tau), eps * r(0));
 }

◆ iso_hardening_exp()

double iso_hardening_exp	(	double	tau,
		double	b_iso
	)

inline

Examples: plastic.cpp.

Definition at line 63 of file plastic.cpp.

                                                           {
   return std::exp(
       std::max(static_cast<double>(std::numeric_limits<float>::min_exponent10),
                -b_iso * tau));
 }

◆ kinematic_hardening()

template<typename T , int DIM>

auto kinematic_hardening	(	FTensor::Tensor2_symmetric< T, DIM > &	t_plastic_strain,
		double	C1_k
	)

inline

Kinematic hardening

Examples: plastic.cpp, and PlasticOpsGeneric.hpp.

Definition at line 91 of file plastic.cpp.

                                  {
   FTensor::Index<'i', DIM> i;
   FTensor::Index<'j', DIM> j;
   FTensor::Tensor2_symmetric<double, DIM> t_alpha;
   if (C1_k < std::numeric_limits<double>::epsilon()) {
     t_alpha(i, j) = 0;
     return t_alpha;
   }
   t_alpha(i, j) = C1_k * t_plastic_strain(i, j);
   return t_alpha;
 }

◆ kinematic_hardening_dplastic_strain()

template<int DIM>

auto kinematic_hardening_dplastic_strain ( double C1_k )

inline

Examples: plastic.cpp.

Definition at line 105 of file plastic.cpp.

                                                              {
   FTensor::Index<'i', DIM> i;
   FTensor::Index<'j', DIM> j;
   FTensor::Index<'k', DIM> k;
   FTensor::Index<'l', DIM> l;
   FTensor::Ddg<double, DIM, DIM> t_diff;
   constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
   t_diff(i, j, k, l) = C1_k * (t_kd(i, k) ^ t_kd(j, l)) / 4.;
   return t_diff;
 }

◆ 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]

[Example]

Examples: plastic.cpp.

Definition at line 1319 of file plastic.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(), "PLASTICITY"));
   core_log->add_sink(
       LogManager::createSink(LogManager::getStrmWorld(), "TIMER"));
   LogManager::setLog("PLASTICITY");
   MOFEM_LOG_TAG("PLASTICITY", "Plasticity");
  
 #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]
  
     //! [Load mesh]
     Simple *simple = m_field.getInterface<Simple>();
     CHKERR simple->getOptions();
     CHKERR simple->loadFile();
     //! [Load mesh]
  
     //! [Example]
     Example ex(m_field);
     CHKERR ex.runProblem();
     //! [Example]
   }
   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

◆ alpha_damping

double alpha_damping = 0

Examples: plastic.cpp.

Definition at line 141 of file plastic.cpp.

◆ AT

constexpr AssemblyType AT

constexpr

Initial value:

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

Examples: plastic.cpp.

Definition at line 44 of file plastic.cpp.

◆ b_iso

double b_iso = 16.93

Saturation exponent.

Examples: plastic.cpp, and PlasticOps.hpp.

Definition at line 128 of file plastic.cpp.

◆ C1_k

double C1_k = 0

Kinematic hardening.

Examples: plastic.cpp, and PlasticOps.hpp.

Definition at line 129 of file plastic.cpp.

◆ cn0

double cn0 = 1

Examples: plastic.cpp, and PlasticOpsGeneric.hpp.

Definition at line 131 of file plastic.cpp.

◆ cn1

double cn1 = 1

Examples: plastic.cpp, and PlasticOpsGeneric.hpp.

Definition at line 132 of file plastic.cpp.

◆ CONTACT_SPACE

constexpr FieldSpace CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE

constexpr

Examples: plastic.cpp.

Definition at line 52 of file plastic.cpp.

◆ ep_order

int ep_order = order - 1

Order of ep files.

Examples: plastic.cpp.

Definition at line 136 of file plastic.cpp.

◆ geom_order

int geom_order = 2

Order if fixed.

Examples: adolc_plasticity.cpp, and plastic.cpp.

Definition at line 137 of file plastic.cpp.

◆ H

double H = 129

Hardening.

Examples: ADOLCPlasticityMaterialModels.hpp, HenckyOps.hpp, plastic.cpp, and PlasticOps.hpp.

Definition at line 124 of file plastic.cpp.

◆ help

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

static

[TestOperators]

Examples: plastic.cpp.

Definition at line 1317 of file plastic.cpp.

◆ is_large_strains

PetscBool is_large_strains = PETSC_TRUE

Large strains.

Examples: plastic.cpp, and PlasticOps.hpp.

Definition at line 116 of file plastic.cpp.

◆ is_quasi_static

PetscBool is_quasi_static = PETSC_TRUE

Examples: NonlinearElasticElementInterface.hpp, and plastic.cpp.

Definition at line 139 of file plastic.cpp.

◆ IT

constexpr IntegrationType IT

constexpr

Initial value:

=

IntegrationType::GAUSS

Examples: plastic.cpp.

Definition at line 47 of file plastic.cpp.

◆ order

int order = 2

Order displacement.

Examples: plastic.cpp.

Definition at line 134 of file plastic.cpp.

◆ poisson_ratio

double poisson_ratio = 0.29

Poisson ratio.

Examples: ADOLCPlasticityMaterialModels.hpp, HenckyOps.hpp, plastic.cpp, PlasticOps.hpp, Remodeling.cpp, seepage.cpp, and thermo_elastic.cpp.

Definition at line 122 of file plastic.cpp.

◆ Qinf

double Qinf = 265

Saturation yield stress.

Examples: plastic.cpp, and PlasticOps.hpp.

Definition at line 127 of file plastic.cpp.

◆ rho

double rho = 0.0

Examples: dynamic_first_order_con_law.cpp, heat_method.cpp, HookeElement.cpp, plastic.cpp, Remodeling.cpp, and Remodeling.hpp.

Definition at line 140 of file plastic.cpp.

◆ scale

double scale = 1.

Examples: cell_forces.cpp, EshelbianOperators.cpp, HenckyOps.hpp, phase.cpp, plastic.cpp, PlasticOps.hpp, test_jacobian_of_simple_contact_element.cpp, and UnsaturatedFlow.hpp.

Definition at line 119 of file plastic.cpp.

◆ set_timer

PetscBool set_timer = PETSC_FALSE

Set timer.

Examples: plastic.cpp.

Definition at line 117 of file plastic.cpp.

◆ sigmaY

double sigmaY = 450

Yield stress.

Examples: ADOLCPlasticityMaterialModels.hpp, plastic.cpp, and PlasticOps.hpp.

Definition at line 123 of file plastic.cpp.

◆ size_symm

constexpr auto size_symm = (SPACE_DIM * (SPACE_DIM + 1)) / 2

constexpr

Examples: eigen_elastic.cpp, EshelbianOperators.cpp, HenckyOps.hpp, plastic.cpp, PlasticOps.hpp, PlasticOpsGeneric.hpp, PlasticOpsLargeStrains.hpp, PlasticOpsSmallStrains.hpp, seepage.cpp, and thermo_elastic.cpp.

Definition at line 42 of file plastic.cpp.

◆ SPACE_DIM

constexpr int SPACE_DIM

constexpr

Initial value:

=

EXECUTABLE_DIMENSION

Examples: plastic.cpp.

Definition at line 40 of file plastic.cpp.

◆ tau_order

int tau_order = order - 2

Order of tau files.

Examples: plastic.cpp.

Definition at line 135 of file plastic.cpp.

◆ visH

double visH = 0

Viscous hardening.

Examples: plastic.cpp, and PlasticOps.hpp.

Definition at line 125 of file plastic.cpp.

◆ young_modulus

double young_modulus = 206913

Young modulus.

Examples: ADOLCPlasticityMaterialModels.hpp, HenckyOps.hpp, plastic.cpp, PlasticOps.hpp, Remodeling.cpp, seepage.cpp, and thermo_elastic.cpp.

Definition at line 121 of file plastic.cpp.

◆ zeta

double zeta = 5e-2

Viscous hardening.

Examples: plastic.cpp, PlasticOpsGeneric.hpp, and plot_base.cpp.

Definition at line 126 of file plastic.cpp.

Classes

Macros

Typedefs

Functions

Variables

Macro Definition Documentation

◆ EXECUTABLE_DIMENSION

Typedef Documentation

◆ BoundaryEle

◆ BoundaryEleOp

◆ BoundaryLhsBCs

◆ BoundaryRhsBCs

◆ DomainEle

◆ DomainEleOp

◆ DomainRhsBCs

◆ EntData

◆ OpBoundaryLhsBCs

◆ OpBoundaryRhsBCs

◆ OpDomainRhsBCs

◆ PostProcEle

◆ SideEle

◆ SkinPostProcEle

Function Documentation

◆ iso_hardening()

◆ iso_hardening_dtau()

◆ iso_hardening_exp()

◆ kinematic_hardening()

◆ kinematic_hardening_dplastic_strain()

◆ main()

Variable Documentation

◆ alpha_damping

◆ AT

◆ b_iso

◆ C1_k

◆ cn0

◆ cn1

◆ CONTACT_SPACE

◆ ep_order

◆ geom_order

◆ H

◆ help

◆ is_large_strains

◆ is_quasi_static

◆ IT

◆ order

◆ poisson_ratio

◆ Qinf

◆ rho

◆ scale

◆ set_timer

◆ sigmaY

◆ size_symm

◆ SPACE_DIM

◆ tau_order

◆ visH

◆ young_modulus

◆ zeta