plastic.cpp File Reference

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

Go to the source code of this file.

Classes
struct	ElementsAndOps< DIM >
struct	ElementsAndOps< 2 >
struct	ElementsAndOps< 3 >
struct	Example
	[Example] More...
struct	Example::PlasticityTimeScale
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	AssemblyDomainEleOp = FormsIntegrators< DomainEleOp >::Assembly< A >::OpBase
using	OpKCauchy = FormsIntegrators< DomainEleOp >::Assembly< A >::BiLinearForm< GAUSS >::OpGradSymTensorGrad< 1, SPACE_DIM, SPACE_DIM, 0 >
	[Only used with Hooke equation (linear material model)] More...
using	OpInternalForceCauchy = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< G >::OpGradTimesSymTensor< 1, SPACE_DIM, SPACE_DIM >
using	OpMass = FormsIntegrators< DomainEleOp >::Assembly< A >::BiLinearForm< GAUSS >::OpMass< 1, SPACE_DIM >
	[Only used with Hooke equation (linear material model)] More...
using	OpInertiaForce = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< G >::OpBaseTimesVector< 1, SPACE_DIM, 1 >
using	OpKPiola = FormsIntegrators< DomainEleOp >::Assembly< A >::BiLinearForm< GAUSS >::OpGradTensorGrad< 1, SPACE_DIM, SPACE_DIM, 1 >
	[Only used for dynamics] More...
using	OpInternalForcePiola = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< G >::OpGradTimesTensor< 1, SPACE_DIM, SPACE_DIM >
using	OpBoundaryMass = FormsIntegrators< BoundaryEleOp >::Assembly< PETSC >::BiLinearForm< G >::OpMass< 1, SPACE_DIM >
	[Only used with Hencky/nonlinear material] More...
using	OpBoundaryVec = FormsIntegrators< BoundaryEleOp >::Assembly< PETSC >::LinearForm< G >::OpBaseTimesVector< 1, SPACE_DIM, 0 >
using	OpBoundaryInternal = FormsIntegrators< BoundaryEleOp >::Assembly< PETSC >::LinearForm< G >::OpBaseTimesVector< 1, SPACE_DIM, 1 >
using	OpScaleL2 = MoFEM::OpScaleBaseBySpaceInverseOfMeasure< DomainEleOp >
	[Essential boundary conditions] More...
using	DomainNaturalBC = NaturalBC< DomainEleOp >::Assembly< A >::LinearForm< G >
using	OpBodyForce = DomainNaturalBC::OpFlux< NaturalMeshsetType< BLOCKSET >, 1, SPACE_DIM >
using	BoundaryNaturalBC = NaturalBC< BoundaryEleOp >::Assembly< A >::LinearForm< G >
using	OpForce = BoundaryNaturalBC::OpFlux< NaturalMeshsetType< BLOCKSET >, 1, SPACE_DIM >
using	OpEssentialLhs = EssentialBC< BoundaryEleOp >::Assembly< A >::BiLinearForm< GAUSS >::OpEssentialLhs< DisplacementCubitBcData, 1, SPACE_DIM >
using	OpEssentialRhs = EssentialBC< BoundaryEleOp >::Assembly< A >::LinearForm< GAUSS >::OpEssentialRhs< DisplacementCubitBcData, 1, SPACE_DIM >

Functions
double	hardening_exp (double tau)
double	hardening (double tau)
double	hardening_dtau (double tau)
double	hardening_dtau2 (double tau)
int	main (int argc, char *argv[])

Variables
constexpr int	SPACE_DIM
constexpr auto	size_symm = (SPACE_DIM * (SPACE_DIM + 1)) / 2
constexpr AssemblyType	A
constexpr IntegrationType	G
PetscBool	is_large_strains = PETSC_TRUE
PetscBool	set_timer = PETSC_FALSE
double	scale = 1.
double	young_modulus = 206913
double	poisson_ratio = 0.29
double	rho = 0
double	sigmaY = 450
double	H = 129
double	visH = 0
double	cn0 = 1
double	cn1 = 1
double	zeta = 5e-2
double	Qinf = 265
double	b_iso = 16.93
int	order = 2
	Order if fixed. More...
int	geom_order = 2
	Order if fixed. More...
constexpr size_t	activ_history_sise = 1
static char	help [] = "...\n\n"
	[Solve] More...

Macro Definition Documentation

EXECUTABLE_DIMENSION

#define EXECUTABLE_DIMENSION   3

Examples

contact.cpp, eigen_elastic.cpp, level_set.cpp, plastic.cpp, plot_base.cpp, scalar_check_approximation.cpp, and thermo_elastic.cpp.

Definition at line 10 of file plastic.cpp.

Typedef Documentation

AssemblyDomainEleOp

using AssemblyDomainEleOp = FormsIntegrators<DomainEleOp>::Assembly<A>::OpBase

Examples

PlasticOpsGeneric.hpp, PlasticOpsLargeStrains.hpp, PlasticOpsSmallStrains.hpp, and ThermoElasticOps.hpp.

Definition at line 48 of file plastic.cpp.

BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 44 of file plastic.cpp.

BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 45 of file plastic.cpp.

BoundaryNaturalBC

using BoundaryNaturalBC = NaturalBC<BoundaryEleOp>::Assembly<A>::LinearForm<G>

Definition at line 85 of file plastic.cpp.

DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 42 of file plastic.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 43 of file plastic.cpp.

DomainNaturalBC

using DomainNaturalBC = NaturalBC<DomainEleOp>::Assembly<A>::LinearForm<G>

Definition at line 81 of file plastic.cpp.

EntData

using EntData = EntitiesFieldData::EntData

Definition at line 41 of file plastic.cpp.

OpBodyForce

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

Definition at line 82 of file plastic.cpp.

OpBoundaryInternal

using OpBoundaryInternal = FormsIntegrators<BoundaryEleOp>::Assembly< PETSC>::LinearForm<G>::OpBaseTimesVector<1, SPACE_DIM, 1>

Examples

contact.cpp.

Definition at line 76 of file plastic.cpp.

OpBoundaryMass

using OpBoundaryMass = FormsIntegrators<BoundaryEleOp>::Assembly< PETSC>::BiLinearForm<G>::OpMass<1, SPACE_DIM>

[Only used with Hencky/nonlinear material]

[Essential boundary conditions]

Examples

contact.cpp, and helmholtz.cpp.

Definition at line 72 of file plastic.cpp.

OpBoundaryVec

using OpBoundaryVec = FormsIntegrators<BoundaryEleOp>::Assembly< PETSC>::LinearForm<G>::OpBaseTimesVector<1, SPACE_DIM, 0>

Examples

contact.cpp.

Definition at line 74 of file plastic.cpp.

OpEssentialLhs

using OpEssentialLhs = EssentialBC<BoundaryEleOp>::Assembly<A>::BiLinearForm< GAUSS>::OpEssentialLhs<DisplacementCubitBcData, 1, SPACE_DIM>

Definition at line 89 of file plastic.cpp.

OpEssentialRhs

using OpEssentialRhs = EssentialBC<BoundaryEleOp>::Assembly<A>::LinearForm< GAUSS>::OpEssentialRhs<DisplacementCubitBcData, 1, SPACE_DIM>

Definition at line 91 of file plastic.cpp.

OpForce

using OpForce = BoundaryNaturalBC::OpFlux<NaturalMeshsetType<BLOCKSET>, 1, SPACE_DIM>

Definition at line 86 of file plastic.cpp.

OpInertiaForce

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

Definition at line 60 of file plastic.cpp.

OpInternalForceCauchy

using OpInternalForceCauchy = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<G>::OpGradTimesSymTensor<1, SPACE_DIM, SPACE_DIM>

Examples

contact.cpp, and plastic.cpp.

Definition at line 53 of file plastic.cpp.

OpInternalForcePiola

using OpInternalForcePiola = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<G>::OpGradTimesTensor<1, SPACE_DIM, SPACE_DIM>

Examples

contact.cpp, and plastic.cpp.

Definition at line 67 of file plastic.cpp.

OpKCauchy

using OpKCauchy = FormsIntegrators<DomainEleOp>::Assembly<A>::BiLinearForm< GAUSS>::OpGradSymTensorGrad<1, SPACE_DIM, SPACE_DIM, 0>

[Only used with Hooke equation (linear material model)]

Examples

contact.cpp, and plastic.cpp.

Definition at line 51 of file plastic.cpp.

OpKPiola

using OpKPiola = FormsIntegrators<DomainEleOp>::Assembly<A>::BiLinearForm< GAUSS>::OpGradTensorGrad<1, SPACE_DIM, SPACE_DIM, 1>

[Only used for dynamics]

[Only used with Hencky/nonlinear material]

Examples

contact.cpp, and plastic.cpp.

Definition at line 65 of file plastic.cpp.

OpMass

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

[Only used with Hooke equation (linear material model)]

[Only used for dynamics]

Definition at line 58 of file plastic.cpp.

OpScaleL2

using OpScaleL2 = MoFEM::OpScaleBaseBySpaceInverseOfMeasure<DomainEleOp>

[Essential boundary conditions]

Definition at line 79 of file plastic.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 46 of file plastic.cpp.

Function Documentation

hardening()

double hardening ( double  tau )

inline

Examples

PlasticOpsGeneric.hpp, and plastic.cpp.

Definition at line 122 of file plastic.cpp.

                                    {
  return H * tau + Qinf * (1. - hardening_exp(tau)) + sigmaY;
}

hardening_dtau()

double hardening_dtau ( double  tau )

inline

Examples

PlasticOpsGeneric.hpp, and plastic.cpp.

Definition at line 126 of file plastic.cpp.

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

hardening_dtau2()

double hardening_dtau2 ( double  tau )

inline

Examples

plastic.cpp.

Definition at line 132 of file plastic.cpp.

                                          {
  return -(Qinf * (b_iso * b_iso)) * hardening_exp(tau);
}

hardening_exp()

double hardening_exp ( double  tau )

inline

Examples

plastic.cpp.

Definition at line 116 of file plastic.cpp.

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

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]

[Example]

[Example]

Definition at line 1186 of file plastic.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"));
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "TIMER"));
  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 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();
}

Variable Documentation

A

constexpr AssemblyType A

constexpr

Initial value:

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

Definition at line 35 of file plastic.cpp.

activ_history_sise

constexpr size_t activ_history_sise = 1

constexpr

Examples

plastic.cpp.

Definition at line 114 of file plastic.cpp.

b_iso

double b_iso = 16.93

Examples

plastic.cpp.

Definition at line 109 of file plastic.cpp.

cn0

double cn0 = 1

Examples

PlasticOps.hpp, and plastic.cpp.

Definition at line 105 of file plastic.cpp.

cn1

double cn1 = 1

Examples

PlasticOps.hpp, and plastic.cpp.

Definition at line 106 of file plastic.cpp.

G

constexpr IntegrationType G

constexpr

Initial value:

=
    IntegrationType::GAUSS

Examples

plastic.cpp.

Definition at line 38 of file plastic.cpp.

geom_order

int geom_order = 2

Order if fixed.

Examples

plastic.cpp.

Definition at line 112 of file plastic.cpp.

H

double H = 129

Examples

plastic.cpp.

Definition at line 103 of file plastic.cpp.

help

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

static

[Solve]

Definition at line 1184 of file plastic.cpp.

is_large_strains

PetscBool is_large_strains = PETSC_TRUE

Examples

contact.cpp, and plastic.cpp.

Definition at line 94 of file plastic.cpp.

order

int order = 2

Order if fixed.

Definition at line 111 of file plastic.cpp.

poisson_ratio

double poisson_ratio = 0.29

Definition at line 100 of file plastic.cpp.

Qinf

double Qinf = 265

Examples

plastic.cpp.

Definition at line 108 of file plastic.cpp.

rho

double rho = 0

Examples

HookeElement.cpp, Remodeling.cpp, Remodeling.hpp, contact.cpp, eigen_elastic.cpp, heat_method.cpp, and plastic.cpp.

Definition at line 101 of file plastic.cpp.

scale

double scale = 1.

Examples

PlasticOps.hpp, PlasticOpsGeneric.hpp, UnsaturatedFlow.hpp, cell_forces.cpp, phase.cpp, plastic.cpp, and test_jacobian_of_simple_contact_element.cpp.

Definition at line 97 of file plastic.cpp.

set_timer

PetscBool set_timer = PETSC_FALSE

Examples

plastic.cpp.

Definition at line 95 of file plastic.cpp.

sigmaY

double sigmaY = 450

Examples

PlasticOps.hpp, and plastic.cpp.

Definition at line 102 of file plastic.cpp.

size_symm

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

constexpr

Examples

PlasticOpsGeneric.hpp, PlasticOpsLargeStrains.hpp, PlasticOpsSmallStrains.hpp, contact.cpp, eigen_elastic.cpp, nonlinear_elastic.cpp, plastic.cpp, and thermo_elastic.cpp.

Definition at line 33 of file plastic.cpp.

SPACE_DIM

constexpr int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

Definition at line 31 of file plastic.cpp.

visH

double visH = 0

Examples

PlasticOps.hpp, and plastic.cpp.

Definition at line 104 of file plastic.cpp.

young_modulus

double young_modulus = 206913

Examples

Remodeling.cpp, contact.cpp, eigen_elastic.cpp, nonlinear_elastic.cpp, plastic.cpp, and thermo_elastic.cpp.

Definition at line 99 of file plastic.cpp.

zeta

double zeta = 5e-2

Examples

PlasticOps.hpp, plastic.cpp, and plot_base.cpp.

Definition at line 107 of file plastic.cpp.