plastic.cpp File Reference

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

Go to the source code of this file.

Classes
struct	ElementsAndOps< 2 >
struct	ElementsAndOps< 3 >
struct	PlasticOps::AddHOOps< 2, 3, 3 >
struct	PlasticOps::AddHOOps< 1, 2, 2 >
struct	PlasticOps::AddHOOps< 3, 3, 3 >
struct	PlasticOps::AddHOOps< 2, 2, 2 >
struct	Example
	[Example] More...
struct	Example::ScaledTimeScale
struct	SetUpSchur
	[Push operators to pipeline] More...
struct	SetUpSchurImpl

Namespaces
namespace	PlasticOps

Macros
#define	EXECUTABLE_DIMENSION   3

Typedefs
using	DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle
using	BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle
using	PostProcEle = PostProcBrokenMeshInMoab<DomainEle>
using	SkinPostProcEle = PostProcBrokenMeshInMoab<BoundaryEle>
using	SideEle = ElementsAndOps<SPACE_DIM>::SideEle
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
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[])
template<int FE_DIM, int PROBLEM_DIM, int SPACE_DIM>
MoFEMErrorCode	PlasticOps::scaleL2 (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pipeline, std::string geom_field_name)

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.
PetscBool	set_timer = PETSC_FALSE
	Set timer.
double	scale = 1.
double	young_modulus = 206913
	Young modulus.
double	poisson_ratio = 0.29
	Poisson ratio.
double	sigmaY = 450
	Yield stress.
double	H = 129
	Hardening.
double	visH = 0
	Viscous hardening.
double	zeta = 5e-2
	Viscous hardening.
double	Qinf = 265
	Saturation yield stress.
double	b_iso = 16.93
	Saturation exponent.
double	C1_k = 0
	Kinematic hardening.
double	cn0 = 1
double	cn1 = 1
int	order = 2
	Order displacement.
int	tau_order = order - 2
	Order of tau files.
int	ep_order = order - 1
	Order of ep files.
int	geom_order = 2
	Order if fixed.
PetscBool	is_quasi_static = PETSC_TRUE
double	rho = 0.0
double	alpha_damping = 0
static char	help [] = "...\n\n"
	[TestOperators]

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, seepage.cpp, tensor_divergence_operator.cpp, test_broken_space.cpp, and thermo_elastic.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.

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.

DomainRhsBCs

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

Definition at line 165 of file plastic.cpp.

OpBoundaryLhsBCs

using OpBoundaryLhsBCs

Initial value:

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

Definition at line 172 of file plastic.cpp.

OpBoundaryRhsBCs

using OpBoundaryRhsBCs

Initial value:

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

Definition at line 169 of file plastic.cpp.

OpDomainRhsBCs

using OpDomainRhsBCs

Initial value:

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

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

PlasticOpsGeneric.hpp, and plastic.cpp.

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

PlasticOpsGeneric.hpp, and plastic.cpp.

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

PlasticOpsGeneric.hpp, and plastic.cpp.

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

PlasticOpsGeneric.hpp, and 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]

[Load mesh]

[Example]

[Example]

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

AssemblyType AT

constexpr

Initial value:

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

Definition at line 44 of file plastic.cpp.

b_iso

double b_iso = 16.93

Saturation exponent.

Examples

PlasticOps.hpp, and plastic.cpp.

Definition at line 128 of file plastic.cpp.

C1_k

double C1_k = 0

Kinematic hardening.

Examples

PlasticOps.hpp, and plastic.cpp.

Definition at line 129 of file plastic.cpp.

cn0

double cn0 = 1

Examples

PlasticOpsGeneric.hpp, and plastic.cpp.

Definition at line 131 of file plastic.cpp.

cn1

double cn1 = 1

Examples

PlasticOpsGeneric.hpp, and plastic.cpp.

Definition at line 132 of file plastic.cpp.

CONTACT_SPACE

FieldSpace CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE

constexpr

Examples

adolc_plasticity.cpp, and 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

plastic.cpp.

Definition at line 137 of file plastic.cpp.

H

double H = 129

Hardening.

Examples

ADOLCPlasticityLargeStrain.hpp, ADOLCPlasticityMaterialModels.hpp, HenckyOps.hpp, PlasticOps.hpp, and plastic.cpp.

Definition at line 124 of file plastic.cpp.

help

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

static

[TestOperators]

Definition at line 1388 of file plastic.cpp.

is_large_strains

PetscBool is_large_strains = PETSC_TRUE

Large strains.

Examples

PlasticOps.hpp, and plastic.cpp.

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

IntegrationType IT

constexpr

Initial value:

=
    IntegrationType::GAUSS

Definition at line 47 of file plastic.cpp.

order

int order = 2

Order displacement.

Definition at line 134 of file plastic.cpp.

poisson_ratio

double poisson_ratio = 0.29

Poisson ratio.

Examples

ADOLCPlasticityMaterialModels.hpp, HenckyOps.hpp, PlasticOps.hpp, Remodeling.cpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, plastic.cpp, seepage.cpp, and thermo_elastic.cpp.

Definition at line 122 of file plastic.cpp.

Qinf

double Qinf = 265

Saturation yield stress.

Examples

PlasticOps.hpp, and plastic.cpp.

Definition at line 127 of file plastic.cpp.

rho

double rho = 0.0

Examples

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

Definition at line 140 of file plastic.cpp.

scale

double scale = 1.

Examples

EshelbianOperators.cpp, HenckyOps.hpp, PlasticOps.hpp, UnsaturatedFlow.hpp, adolc_plasticity.cpp, cell_forces.cpp, phase.cpp, plastic.cpp, seepage.cpp, and test_jacobian_of_simple_contact_element.cpp.

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, PlasticOps.hpp, and plastic.cpp.

Definition at line 123 of file plastic.cpp.

size_symm

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

constexpr

Examples

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

Definition at line 42 of file plastic.cpp.

SPACE_DIM

int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

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

PlasticOps.hpp, and plastic.cpp.

Definition at line 125 of file plastic.cpp.

young_modulus

double young_modulus = 206913

Young modulus.

Examples

ADOLCPlasticityMaterialModels.hpp, HenckyOps.hpp, PlasticOps.hpp, Remodeling.cpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, plastic.cpp, seepage.cpp, and thermo_elastic.cpp.

Definition at line 121 of file plastic.cpp.

zeta

double zeta = 5e-2

Viscous hardening.

Examples

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

Definition at line 126 of file plastic.cpp.