plastic.cpp File Reference

#include <MoFEM.hpp>
#include <MatrixFunction.hpp>
#include <IntegrationRules.hpp>
#include <cblas.h>
#include <lapack_wrap.h>
#include <quad.h>
#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

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
long double	hardening (long double tau)
long double	hardening_dtau (long double tau)
long double	hardening_dtau2 (long 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 = AssemblyType::PETSC
constexpr IntegrationType	G
PetscBool	is_large_strains = PETSC_TRUE
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	cn = 1
double	zeta = 5e-2
double	eta = 1e2
double	Qinf = 265
double	b_iso = 16.93
int	order = 2
	Order if fixed. More...
int	geom_order = 2
	Order if fixed. More...
static char	help [] = "...\n\n"
	[Solve] More...

Macro Definition Documentation

EXECUTABLE_DIMENSION

#define EXECUTABLE_DIMENSION   3

Examples

contact.cpp, eigen_elastic.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 57 of file plastic.cpp.

BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 53 of file plastic.cpp.

BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 54 of file plastic.cpp.

BoundaryNaturalBC

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

Definition at line 94 of file plastic.cpp.

DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 51 of file plastic.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 52 of file plastic.cpp.

DomainNaturalBC

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

Definition at line 90 of file plastic.cpp.

EntData

using EntData = EntitiesFieldData::EntData

Definition at line 50 of file plastic.cpp.

OpBodyForce

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

Definition at line 91 of file plastic.cpp.

OpBoundaryInternal

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

Examples

contact.cpp.

Definition at line 85 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 81 of file plastic.cpp.

OpBoundaryVec

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

Examples

contact.cpp.

Definition at line 83 of file plastic.cpp.

OpEssentialLhs

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

Definition at line 98 of file plastic.cpp.

OpEssentialRhs

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

Definition at line 100 of file plastic.cpp.

OpForce

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

Definition at line 95 of file plastic.cpp.

OpInertiaForce

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

Definition at line 69 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 62 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 76 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 60 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 74 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 67 of file plastic.cpp.

OpScaleL2

using OpScaleL2 = MoFEM::OpScaleBaseBySpaceInverseOfMeasure<DomainEleOp>

[Essential boundary conditions]

Definition at line 88 of file plastic.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 55 of file plastic.cpp.

Function Documentation

hardening()

long double hardening ( long double  tau )

inline

Examples

PlasticOpsGeneric.hpp, and plastic.cpp.

Definition at line 122 of file plastic.cpp.

                                              {
  return H * tau + Qinf * (1. - std::exp(-b_iso * tau)) + sigmaY;
}

hardening_dtau()

long double hardening_dtau ( long double  tau )

inline

Examples

PlasticOpsGeneric.hpp, and plastic.cpp.

Definition at line 126 of file plastic.cpp.

                                                   {
  return H + Qinf * b_iso * std::exp(-b_iso * tau);
}

hardening_dtau2()

long double hardening_dtau2 ( long double  tau )

inline

Examples

PlasticOpsGeneric.hpp, and plastic.cpp.

Definition at line 130 of file plastic.cpp.

                                                    {
  return -(Qinf * (b_iso * b_iso)) * std::exp(-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 insterface

[Create MoFEM]

[Load mesh]

[Load mesh]

[Example]

[Example]

Definition at line 983 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"));
  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 insterface
    //! [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 = AssemblyType::PETSC

constexpr

Examples

Remodeling.cpp, analytical_poisson_field_split.cpp, approx_sphere.cpp, build_large_problem.cpp, cell_forces.cpp, contact.cpp, dm_create_subdm.cpp, eigen_elastic.cpp, free_surface.cpp, minimal_surface_area.cpp, nonlinear_elastic.cpp, photon_diffusion.cpp, plastic.cpp, plate.cpp, prism_polynomial_approximation.cpp, quad_polynomial_approximation.cpp, simple_elasticity.cpp, simple_interface.cpp, test_jacobian_of_simple_contact_element.cpp, testing_jacobian_of_hook_element.cpp, testing_jacobian_of_hook_scaled_with_density_element.cpp, and thermo_elastic.cpp.

Definition at line 46 of file plastic.cpp.

b_iso

double b_iso = 16.93

Examples

plastic.cpp.

Definition at line 117 of file plastic.cpp.

cn

double cn = 1

Examples

ContactOps.hpp, PlasticOps.hpp, contact.cpp, and plastic.cpp.

Definition at line 113 of file plastic.cpp.

eta

double eta = 1e2

Examples

PlasticOpsGeneric.hpp, free_surface.cpp, plastic.cpp, and plot_base.cpp.

Definition at line 115 of file plastic.cpp.

G

constexpr IntegrationType G

constexpr

Initial value:

=
    IntegrationType::GAUSS

Examples

plastic.cpp.

Definition at line 47 of file plastic.cpp.

geom_order

int geom_order = 2

Order if fixed.

Examples

plastic.cpp.

Definition at line 120 of file plastic.cpp.

H

double H = 129

Examples

plastic.cpp.

Definition at line 111 of file plastic.cpp.

help

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

static

[Solve]

Definition at line 981 of file plastic.cpp.

is_large_strains

PetscBool is_large_strains = PETSC_TRUE

Examples

contact.cpp, and plastic.cpp.

Definition at line 103 of file plastic.cpp.

order

int order = 2

Order if fixed.

Definition at line 119 of file plastic.cpp.

poisson_ratio

double poisson_ratio = 0.29

Definition at line 108 of file plastic.cpp.

Qinf

double Qinf = 265

Examples

plastic.cpp.

Definition at line 116 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 109 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 105 of file plastic.cpp.

sigmaY

double sigmaY = 450

Examples

PlasticOps.hpp, and plastic.cpp.

Definition at line 110 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 44 of file plastic.cpp.

SPACE_DIM

constexpr int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

Definition at line 42 of file plastic.cpp.

visH

double visH = 0

Examples

PlasticOps.hpp, and plastic.cpp.

Definition at line 112 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 107 of file plastic.cpp.

zeta

double zeta = 5e-2

Examples

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

Definition at line 114 of file plastic.cpp.