#include <MoFEM.hpp>
#include <MatrixFunction.hpp>
#include <IntegrationRules.hpp>
#include <cstdlib>
#include <phg-quadrule/quad.h>
#include <HenckyOps.hpp>
#include <PlasticOps.hpp>
#include <PlasticNaturalBCs.hpp>
#include <ThermoElasticOps.hpp>
#include <FiniteThermalOps.hpp>
#include <ThermoPlasticOps.hpp>
#include <ThermalOps.hpp>
#include <EdgeFlippingOps.hpp>
#include <SolutionMapping.hpp>
#include <ThermalConvection.hpp>
#include <ThermalRadiation.hpp>

Classes
struct	ElementsAndOps< 2 >

struct	ElementsAndOps< 3 >

struct	TSPrePostProc
	Set of functions called by PETSc solver used to refine and update mesh. More...

struct	ResizeCtx

struct	Example
	[Example] More...

struct	Example::ScaledTimeScale

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

struct	MyTsCtx

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	SetPtsData = FieldEvaluatorInterface::SetPtsData

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

using	ScalerFunTwoArgs = boost::function< double(const double, const double)>

using	ScalerFunThreeArgs = boost::function< double(const double, const double, const double)>

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 >

using	DomainNaturalBCRhs = NaturalBC< DomainEleOp >::Assembly< AT >::LinearForm< IT >
	[Body and heat source]

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

using	OpHeatSource = DomainNaturalBCRhs::OpFlux< NaturalMeshsetType< BLOCKSET >, 1, 1 >

using	DomainNaturalBCLhs = NaturalBC< DomainEleOp >::Assembly< AT >::BiLinearForm< IT >

using	BoundaryNaturalBC = NaturalBC< BoundaryEleOp >::Assembly< AT >::LinearForm< IT >
	[Body and heat source]

using	OpForce = BoundaryNaturalBC::OpFlux< NaturalForceMeshsets, 1, SPACE_DIM >

using	OpTemperatureBC = BoundaryNaturalBC::OpFlux< NaturalTemperatureMeshsets, 3, SPACE_DIM >

using	OpEssentialFluxRhs = EssentialBC< BoundaryEleOp >::Assembly< AT >::LinearForm< IT >::OpEssentialRhs< HeatFluxCubitBcData, 3, SPACE_DIM >
	[Natural boundary conditions]

using	OpEssentialFluxLhs = EssentialBC< BoundaryEleOp >::Assembly< AT >::BiLinearForm< IT >::OpEssentialLhs< HeatFluxCubitBcData, 3, SPACE_DIM >

using	OpSetTemperatureRhs = DomainNaturalBCRhs::OpFlux< SetTargetTemperature, 1, 1 >
	[Essential boundary conditions (Least square approach)]

using	OpSetTemperatureLhs = DomainNaturalBCLhs::OpFlux< SetTargetTemperature, 1, 1 >

Enumerations
enum	ICType { IC_UNIFORM , IC_GAUSSIAN , IC_LINEAR }

enum	ThermoPlasticityBits { COMPUTATION_BIT = 0 , FLIPPED_BIT = 1 , FIRST_REF_BIT = 2 , REFINED_EDGES_BIT = BITREFLEVEL_SIZE - 4 , PREVIOUS_BIT = BITREFLEVEL_SIZE - 3 , VIRGIN_BIT = BITREFLEVEL_SIZE - 2 , STORAGE_BIT = BITREFLEVEL_SIZE - 1 }

Functions
double	H_thermal (double H, double omega_h, double temp_0, double temp)

double	dH_thermal_dtemp (double H, double omega_h)

double	y_0_thermal (double sigmaY, double omega_0, double temp_0, double temp)

double	dy_0_thermal_dtemp (double sigmaY, double omega_0)

double	Qinf_thermal (double Qinf, double omega_h, double temp_0, double temp)

double	dQinf_thermal_dtemp (double Qinf, double omega_h)

double	iso_hardening_exp (double tau, double b_iso)

double	temp_hat (double temp)

double	exp_softening (double temp_hat)

double	iso_hardening (double tau, double H, double omega_0, double Qinf, double omega_h, double b_iso, double sigmaY, double temp_0, double temp)

double	iso_hardening_dtau (double tau, double H, double omega_0, double Qinf, double omega_h, double b_iso, double sigmaY, double temp_0, double temp)

double	iso_hardening_dtemp (double tau, double H, double omega_0, double Qinf, double omega_h, double b_iso, double sigmaY, double temp_0, double temp)

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)

auto	postProcessHere (MoFEM::Interface &m_field, SmartPetscObj< DM > &dm, std::string output_name, int &counter= *(new int(0)))

auto	postProcessPETScHere (MoFEM::Interface &m_field, SmartPetscObj< DM > &dm, Vec sol, std::string output_name)

auto	printOnAllCores (MoFEM::Interface &m_field, const std::string &out_put_string, const auto &out_put_quantity)

PetscErrorCode	MyTSResizeTransfer (TS, PetscInt, Vec[], Vec[], void *)

MoFEMErrorCode	TSIterationPreStage (TS ts, PetscReal stagetime)

MoFEMErrorCode	SNESIterationMonitor (SNES snes, PetscInt its, PetscReal norm, void *ctx)

PetscErrorCode	MyTSResizeSetup (TS ts, PetscInt nstep, PetscReal time, Vec sol, PetscBool resize, void ctx)

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

Variables
constexpr int	SPACE_DIM

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

constexpr bool	IS_LARGE_STRAINS

constexpr AssemblyType	AT

constexpr IntegrationType	IT

constexpr FieldSpace	CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE

double	exp_C = -6.284

double	exp_D = -1.024

const bool	is_large_strains

PetscBool	set_timer = PETSC_FALSE
	Set timer.

PetscBool	do_eval_field = PETSC_FALSE
	Evaluate field.

char	restart_file [255]

PetscBool	restart_flg = PETSC_TRUE

double	restart_time = 0

int	restart_step = 0

PetscBool	is_distributed_mesh = PETSC_TRUE

double	scale = 1.

double	default_thermal_conductivity_scale = 1.

double	default_heat_capacity_scale = 1.

ScalerFunTwoArgs	thermal_conductivity_scaling

ScalerFunTwoArgs	heat_capacity_scaling

ScalerFunThreeArgs	inelastic_heat_fraction_scaling

double	young_modulus = 115000
	Young modulus.

double	poisson_ratio = 0.31
	Poisson ratio.

double	sigmaY = 936.2
	Yield stress.

double	H = 584.3
	Hardening.

double	visH = 0
	Viscous hardening.

double	zeta = 5e-2
	regularisation parameter

double	Qinf = 174.2
	Saturation yield stress.

double	b_iso = 63.69
	Saturation exponent.

double	C1_k = 0
	Kinematic hardening.

double	cn0 = 1

double	cn1 = 1

double	default_ref_temp = 20.0

const char *const	ICTypes []

ICType	ic_type = IC_UNIFORM

double	init_temp = 20.0

double	peak_temp = 1000.0

double	width = 10.0
	Width of Gaussian distribution.

double	default_coeff_expansion = 1e-5

double	default_heat_conductivity

double	default_heat_capacity = 10*2.332

double	omega_0 = 2e-3
	flow stress softening

double	omega_h = 2e-3
	hardening softening

double	inelastic_heat_fraction
	fraction of plastic dissipation converted to heat

double	temp_0 = 20
	reference temperature for thermal softening

int	order = 2
	Order displacement.

int	tau_order = order - 2
	Order of tau field.

int	ep_order = order - 1
	Order of ep field.

int	flux_order = order
	Order of tau field.

int	T_order = order - 1
	Order of ep field.

int	geom_order = 2
	Order if fixed.

int	atom_test = 0

PetscBool	order_edge = PETSC_FALSE

PetscBool	order_face = PETSC_FALSE

PetscBool	order_volume = PETSC_FALSE

PetscBool	is_quasi_static = PETSC_TRUE

double	rho = 0.0

double	alpha_damping = 0

double	init_dt = 0.05

double	min_dt = 1e-12

double	qual_tol = 0

double	qual_thresh = 0.1

double	edge_growth_thresh = 0.2

int	num_refinement_levels = 0

BitRefLevel	comp_mesh_bit = BitRefLevel().set(COMPUTATION_BIT)

BitRefLevel	flipped_bit = BitRefLevel().set(FLIPPED_BIT)

BitRefLevel	refined_bit = BitRefLevel().set(FIRST_REF_BIT)

BitRefLevel	prev_mesh_bit = BitRefLevel().set(PREVIOUS_BIT)

BitRefLevel	virgin_mesh_bit = BitRefLevel().set(VIRGIN_BIT)

BitRefLevel	storage_bit = BitRefLevel().set(STORAGE_BIT)

auto	Gaussian_distribution

auto	linear_distribution

auto	uniform_distribution

auto	init_T
	Initialisation function for temperature field.

int	post_processing_counter = 0

auto	save_range

auto	get_string_from_vector

static boost::weak_ptr< TSPrePostProc >	tsPrePostProc

static std::unordered_map< TS, ResizeCtx * >	ts_to_rctx

static std::unordered_map< TS, MyTsCtx * >	ts_to_ctx

Example *	thermoplastic_raw_ptr = nullptr

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

Macro Definition Documentation

◆ EXECUTABLE_DIMENSION

#define EXECUTABLE_DIMENSION 3

Definition at line 13 of file thermoplastic.cpp.

Typedef Documentation

◆ AssemblyDomainEleOp

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

Definition at line 71 of file thermoplastic.cpp.

◆ BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 64 of file thermoplastic.cpp.

◆ BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 65 of file thermoplastic.cpp.

◆ BoundaryLhsBCs

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

Definition at line 590 of file thermoplastic.cpp.

◆ BoundaryNaturalBC

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

[Body and heat source]

[Natural boundary conditions]

Definition at line 605 of file thermoplastic.cpp.

◆ BoundaryRhsBCs

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

Definition at line 587 of file thermoplastic.cpp.

◆ DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 62 of file thermoplastic.cpp.

◆ DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 63 of file thermoplastic.cpp.

◆ DomainNaturalBCLhs

using DomainNaturalBCLhs = NaturalBC<DomainEleOp>::Assembly<AT>::BiLinearForm<IT>

Definition at line 600 of file thermoplastic.cpp.

◆ DomainNaturalBCRhs

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

[Body and heat source]

Definition at line 595 of file thermoplastic.cpp.

◆ DomainRhsBCs

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

Definition at line 584 of file thermoplastic.cpp.

◆ EntData

using EntData = EntitiesFieldData::EntData

Definition at line 61 of file thermoplastic.cpp.

◆ OpBodyForce

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

Definition at line 596 of file thermoplastic.cpp.

◆ OpBoundaryLhsBCs

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

Definition at line 591 of file thermoplastic.cpp.

◆ OpBoundaryRhsBCs

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

Definition at line 588 of file thermoplastic.cpp.

◆ OpDomainRhsBCs

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

Definition at line 585 of file thermoplastic.cpp.

◆ OpEssentialFluxLhs

using OpEssentialFluxLhs = EssentialBC<BoundaryEleOp>::Assembly< AT>::BiLinearForm<IT>::OpEssentialLhs<HeatFluxCubitBcData, 3, SPACE_DIM>

Definition at line 615 of file thermoplastic.cpp.

◆ OpEssentialFluxRhs

using OpEssentialFluxRhs = EssentialBC<BoundaryEleOp>::Assembly<AT>::LinearForm< IT>::OpEssentialRhs<HeatFluxCubitBcData, 3, SPACE_DIM>

[Natural boundary conditions]

[Essential boundary conditions (Least square approach)]

Definition at line 613 of file thermoplastic.cpp.

◆ OpForce

using OpForce = BoundaryNaturalBC::OpFlux<NaturalForceMeshsets, 1, SPACE_DIM>

Definition at line 607 of file thermoplastic.cpp.

◆ OpHeatSource

using OpHeatSource = DomainNaturalBCRhs::OpFlux<NaturalMeshsetType<BLOCKSET>, 1, 1>

Definition at line 598 of file thermoplastic.cpp.

◆ OpSetTemperatureLhs

using OpSetTemperatureLhs = DomainNaturalBCLhs::OpFlux<SetTargetTemperature, 1, 1>

Definition at line 621 of file thermoplastic.cpp.

◆ OpSetTemperatureRhs

using OpSetTemperatureRhs = DomainNaturalBCRhs::OpFlux<SetTargetTemperature, 1, 1>

[Essential boundary conditions (Least square approach)]

Definition at line 619 of file thermoplastic.cpp.

◆ OpTemperatureBC

using OpTemperatureBC = BoundaryNaturalBC::OpFlux<NaturalTemperatureMeshsets, 3, SPACE_DIM>

Definition at line 608 of file thermoplastic.cpp.

◆ PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 66 of file thermoplastic.cpp.

◆ ScalerFunThreeArgs

using ScalerFunThreeArgs = boost::function<double(const double, const double, const double)>

Definition at line 74 of file thermoplastic.cpp.

◆ ScalerFunTwoArgs

using ScalerFunTwoArgs = boost::function<double(const double, const double)>

Definition at line 73 of file thermoplastic.cpp.

◆ SetPtsData

using SetPtsData = FieldEvaluatorInterface::SetPtsData

Definition at line 69 of file thermoplastic.cpp.

◆ SideEle

using SideEle = ElementsAndOps<SPACE_DIM>::SideEle

Definition at line 68 of file thermoplastic.cpp.

◆ SkinPostProcEle

using SkinPostProcEle = PostProcBrokenMeshInMoab<BoundaryEle>

Definition at line 67 of file thermoplastic.cpp.

Enumeration Type Documentation

◆ ICType

enum ICType

Enumerator
IC_UNIFORM
IC_GAUSSIAN
IC_LINEAR

Definition at line 279 of file thermoplastic.cpp.

279{ IC_UNIFORM, IC_GAUSSIAN, IC_LINEAR };

IC_UNIFORM

@ IC_UNIFORM

Definition thermoplastic.cpp:279

IC_LINEAR

@ IC_LINEAR

Definition thermoplastic.cpp:279

IC_GAUSSIAN

@ IC_GAUSSIAN

Definition thermoplastic.cpp:279

◆ ThermoPlasticityBits

enum ThermoPlasticityBits

Enumerator
COMPUTATION_BIT
FLIPPED_BIT
FIRST_REF_BIT
REFINED_EDGES_BIT
PREVIOUS_BIT
VIRGIN_BIT
STORAGE_BIT

Definition at line 323 of file thermoplastic.cpp.

                          {
  COMPUTATION_BIT = 0,
  FLIPPED_BIT = 1,
  FIRST_REF_BIT = 2,
  REFINED_EDGES_BIT = BITREFLEVEL_SIZE - 4, 
  PREVIOUS_BIT = BITREFLEVEL_SIZE - 3,
  VIRGIN_BIT = BITREFLEVEL_SIZE - 2,
  STORAGE_BIT = BITREFLEVEL_SIZE - 1
};

Function Documentation

◆ dH_thermal_dtemp()

double dH_thermal_dtemp	(	double	H,
		double	omega_h
	)

inline

Examples: thermoplastic.cpp.

Definition at line 148 of file thermoplastic.cpp.

                                                         {
  return -H * omega_h;
}

◆ dQinf_thermal_dtemp()

double dQinf_thermal_dtemp	(	double	Qinf,
		double	omega_h
	)

inline

Examples: thermoplastic.cpp.

Definition at line 166 of file thermoplastic.cpp.

                                                               {
  return -Qinf * omega_h;
}

◆ dy_0_thermal_dtemp()

double dy_0_thermal_dtemp	(	double	sigmaY,
		double	omega_0
	)

inline

Examples: thermoplastic.cpp.

Definition at line 157 of file thermoplastic.cpp.

                                                                {
  return -sigmaY * omega_0;
}

◆ exp_softening()

double exp_softening ( double temp_hat )

inline

Examples: thermoplastic.cpp.

Definition at line 184 of file thermoplastic.cpp.

                                             {
  return std::exp(exp_C * pow(temp_hat, 2) + exp_D * temp_hat);
}

◆ H_thermal()

double H_thermal	(	double	H,
		double	omega_h,
		double	temp_0,
		double	temp
	)

inline

Examples: thermoplastic.cpp.

Definition at line 144 of file thermoplastic.cpp.

                                                                              {
  return H * (1. - omega_h * (temp - temp_0));
}

◆ iso_hardening()

double iso_hardening	(	double	tau,
		double	H,
		double	omega_0,
		double	Qinf,
		double	omega_h,
		double	b_iso,
		double	sigmaY,
		double	temp_0,
		double	temp
	)

inline

Isotropic hardening

Definition at line 191 of file thermoplastic.cpp.

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

◆ iso_hardening_dtau()

double iso_hardening_dtau	(	double	tau,
		double	H,
		double	omega_0,
		double	Qinf,
		double	omega_h,
		double	b_iso,
		double	sigmaY,
		double	temp_0,
		double	temp
	)

inline

Definition at line 198 of file thermoplastic.cpp.

                                                                            {
  return (H + Qinf * b_iso * iso_hardening_exp(tau, b_iso)) *
         exp_softening(temp_hat(temp));
}

◆ iso_hardening_dtemp()

double iso_hardening_dtemp	(	double	tau,
		double	H,
		double	omega_0,
		double	Qinf,
		double	omega_h,
		double	b_iso,
		double	sigmaY,
		double	temp_0,
		double	temp
	)

inline

Examples: PlasticOpsGeneric.hpp, and thermoplastic.cpp.

Definition at line 205 of file thermoplastic.cpp.

                                                                             {
  double JC_ref_temp = 0.;
  double JC_melt_temp = 1650.;
  double T_hat = temp_hat(temp);
 
  return (sigmaY + H * tau + Qinf * (1. - iso_hardening_exp(tau, b_iso))) *
         (exp_D + 2 * T_hat * exp_C) *
         std::exp(exp_D * T_hat + pow(T_hat, 2) * exp_C) /
         (JC_melt_temp - JC_ref_temp);
}

◆ iso_hardening_exp()

double iso_hardening_exp	(	double	tau,
		double	b_iso
	)

inline

Definition at line 170 of file thermoplastic.cpp.

                                                          {
  return std::exp(-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

Definition at line 223 of file thermoplastic.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

Definition at line 237 of file thermoplastic.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]

Definition at line 5743 of file thermoplastic.cpp.

                                 {
 
#ifdef ADD_CONTACT
  #ifdef ENABLE_PYTHON_BINDING
  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(), "THERMAL"));
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "THERMOPLASTICITY"));
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "REMESHING"));
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "TIMER"));
  LogManager::setLog("PLASTICITY");
  LogManager::setLog("THERMAL");
  LogManager::setLog("THERMOPLASTICITY");
  LogManager::setLog("REMESHING");
  MOFEM_LOG_TAG("PLASTICITY", "Plasticity");
  MOFEM_LOG_TAG("THERMAL", "Thermal");
  MOFEM_LOG_TAG("THERMOPLASTICITY", "Thermoplasticity");
  MOFEM_LOG_TAG("REMESHING", "Remeshing");
 
#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();
 
    simple->getBitRefLevel() = comp_mesh_bit;
    simple->getBitRefLevelMask() = BitRefLevel().set();
 
    CHKERR PetscOptionsGetBool(PETSC_NULLPTR, "", "-is_distributed_mesh",
                               &is_distributed_mesh, PETSC_NULLPTR);
    if (is_distributed_mesh == PETSC_TRUE) {
      CHKERR simple->loadFile();
    } else {
      char meshFileName[255];
      CHKERR PetscOptionsGetString(PETSC_NULLPTR, PETSC_NULLPTR, "-file_name",
                                   meshFileName, 255, PETSC_NULLPTR);
      CHKERR simple->loadFile("", meshFileName);
    }
 
    Range ents;
    CHKERR m_field.get_moab().get_entities_by_dimension(0, SPACE_DIM, ents,
                                                        true);
    CHKERR m_field.getInterface<BitRefManager>()->setBitRefLevel(
        ents, comp_mesh_bit | prev_mesh_bit | virgin_mesh_bit | storage_bit,
        false);
    //! [Load mesh]
 
    //! [Example]
    Example ex(m_field);
    CHKERR ex.runProblem();
    //! [Example]
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
 
#ifdef ADD_CONTACT
  #ifdef ENABLE_PYTHON_BINDING
  if (Py_FinalizeEx() < 0) {
    exit(120);
  }
  #endif
#endif // ADD_CONTACT
 
  return 0;
}

◆ MyTSResizeSetup()

PetscErrorCode MyTSResizeSetup	(	TS	ts,
		PetscInt	nstep,
		PetscReal	time,
		Vec	sol,
		PetscBool *	resize,
		void *	ctx
	)

Examples: thermoplastic.cpp.

Definition at line 4216 of file thermoplastic.cpp.

                                                             {
  PetscFunctionBegin;
  ResizeCtx *rctx = static_cast<ResizeCtx *>(ctx);
  *resize = rctx->mesh_changed;
  PetscFunctionReturn(0);
}

◆ MyTSResizeTransfer()

PetscErrorCode MyTSResizeTransfer	(	TS	ts,
		PetscInt	nv,
		Vec	ts_vecsin[],
		Vec	ts_vecsout[],
		void *	ctx
	)

extern

Examples: thermoplastic.cpp.

Definition at line 4224 of file thermoplastic.cpp.

                                                               {
  PetscFunctionBegin;
  ResizeCtx *rctx = static_cast<ResizeCtx *>(ctx);
 
  if (auto ptr = rctx->ptr.lock()) {
 
    MoFEM::Interface &m_field = *(rctx->m_field);
 
    auto bit_mng = m_field.getInterface<BitRefManager>();
 
    auto set_prev_mesh_bit = [](EntityHandle ent, BitRefLevel &bit) {
      if (bit[COMPUTATION_BIT]) {
        bit.set(COMPUTATION_BIT, false);
        bit.set(PREVIOUS_BIT, true);
      }
    };
 
    CHKERR bit_mng->lambdaBitRefLevel(set_prev_mesh_bit);
 
    std::string improved_reasons = "";
 
    Range edges_to_not_refine;
 
    CHKERR ptr->ExRawPtr->doEdgeFlips(rctx->el_q_map, rctx->flipped_els,
                                      rctx->th_spatial_coords,
                                      rctx->new_connectivity);
    if (rctx->el_q_map.size() > 0) {
      improved_reasons += ", Edge flipping";
    }
 
    bool refined = false;
    int step;
    CHKERR TSGetStepNumber(ts, &step);
    CHKERR ptr->ExRawPtr->doEdgeSplits(refined, (step) ? true : false);
    if (refined) {
      improved_reasons += ", Edge splitting";
    }
 
    // Erase leading ", " from improved_reasons
    improved_reasons.erase(0, 2);
 
    if (improved_reasons != "")
      MOFEM_LOG("REMESHING", Sev::inform)
          << "Improved mesh quality by: " << improved_reasons;
 
#ifndef NDEBUG
 
    // comp_mesh_bit = BitRefLevel().set(COMPUTATION_BIT);
    // flipped_bit = BitRefLevel().set(FLIPPED_BIT);
    // refined_bit = BitRefLevel().set(FIRST_REF_BIT);
    // prev_mesh_bit = BitRefLevel().set(PREVIOUS_BIT);
    // virgin_mesh_bit = BitRefLevel().set(VIRGIN_BIT);
    // storage_bit = BitRefLevel().set(STORAGE_BIT);
 
    if (m_field.get_comm_rank() == 0) {
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          comp_mesh_bit, BitRefLevel().set(), 2,
          "before_mapping_comp_mesh_bit.vtk", "VTK", "");
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          flipped_bit, BitRefLevel().set(), 2, "before_mapping_flipped_bit.vtk",
          "VTK", "");
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          refined_bit, BitRefLevel().set(), 2, "before_mapping_refined_bit.vtk",
          "VTK", "");
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          prev_mesh_bit, BitRefLevel().set(), 2,
          "before_mapping_prev_mesh_bit.vtk", "VTK", "");
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          virgin_mesh_bit, BitRefLevel().set(), 2,
          "before_mapping_virgin_mesh_bit.vtk", "VTK", "");
    }
#endif // NDEBUG
 
    MOFEM_LOG("REMESHING", Sev::verbose) << "number of vectors to map = " << nv;
 
    for (PetscInt i = 0; i < nv; ++i) {
      double nrm;
      CHKERR VecNorm(ts_vecsin[i], NORM_2, &nrm);
      MOFEM_LOG("REMESHING", Sev::warning)
          << "Before remeshing: ts_vecsin[" << i << "] norm = " << nrm;
    }
 
    auto scatter_mng = m_field.getInterface<VecManager>();
    auto simple = m_field.getInterface<Simple>();
 
    // CHKERR m_field.getInterface<BitRefManager>()->writeBitLevel(
    //     BitRefLevel().set(1), BitRefLevel().set(), "before_resettingup.vtk",
    //     "VTK", "");
 
    // Rebuilding DM with old and new mesh entities
    auto bit = BitRefLevel().set();
    CHKERR SolutionMapping(m_field).reSetUp(
        {"TAU", "EP", "T"}, {tau_order, ep_order, T_order}, {"U"},
        {order}, {"FLUX"}, {flux_order}, bit);
 
    // CHKERR m_field.getInterface<BitRefManager>()->writeBitLevel(
    //     BitRefLevel().set(1), BitRefLevel().set(), "after_resettingup.vtk",
    //     "VTK", "");
 
    auto intermediate_dm = createDM(m_field.get_comm(), "DMMOFEM");
    CHKERR DMMoFEMCreateMoFEM(
        intermediate_dm, &m_field, "INTERMEDIATE_DM",
        prev_mesh_bit | BitRefLevel().set(FLIPPED_BIT + num_refinement_levels),
        BitRefLevel().set());
    CHKERR DMMoFEMSetDestroyProblem(intermediate_dm, PETSC_TRUE);
    CHKERR DMSetFromOptions(intermediate_dm);
    CHKERR DMMoFEMAddElement(intermediate_dm, simple->getDomainFEName());
    CHKERR DMMoFEMSetSquareProblem(intermediate_dm, PETSC_TRUE);
    CHKERR DMMoFEMSetIsPartitioned(intermediate_dm, is_distributed_mesh);
    CHKERR DMSetUp(intermediate_dm);
 
    Vec vec_in[nv], vec_out[nv];
    for (PetscInt i = 0; i < nv; ++i) {
      CHKERR DMCreateGlobalVector(intermediate_dm, &vec_in[i]);
      CHKERR VecDuplicate(vec_in[i], &vec_out[i]);
    }
 
    VecScatter scatter_to_intermediate;
 
    for (PetscInt i = 0; i < nv; ++i) {
      CHKERR scatter_mng->vecScatterCreate(
          ts_vecsin[i], simple->getProblemName(), RowColData::ROW, vec_in[i],
          "INTERMEDIATE_DM", RowColData::ROW, &scatter_to_intermediate);
      CHKERR VecScatterBegin(scatter_to_intermediate, ts_vecsin[i], vec_in[i],
                             INSERT_VALUES, SCATTER_FORWARD);
      CHKERR VecScatterEnd(scatter_to_intermediate, ts_vecsin[i], vec_in[i],
                           INSERT_VALUES, SCATTER_FORWARD);
    }
 
    CHKERR VecScatterDestroy(&scatter_to_intermediate);
 
    CHKERR SolutionMapping(m_field).mapFields(
      intermediate_dm, prev_mesh_bit,
      BitRefLevel().set(FLIPPED_BIT + num_refinement_levels), nv, vec_in,
      vec_out);
 
    // Build problem on new bit ref level
    auto ts_problem_name = simple->getProblemName();
    CHKERR m_field.modify_problem_ref_level_set_bit(
        ts_problem_name,
        BitRefLevel().set(FLIPPED_BIT + num_refinement_levels));
    CHKERR m_field.modify_problem_mask_ref_level_set_bit(ts_problem_name,
                                                         BitRefLevel().set());
    CHKERR DMMoFEMSetIsPartitioned(simple->getDM(), is_distributed_mesh);
    PetscBarrier(nullptr);
    CHKERR DMSetUp_MoFEM(simple->getDM());
 
    // Update meshsets for correct BCs
    MeshsetsManager *meshsets_manager_ptr;
    CHKERR m_field.getInterface(meshsets_manager_ptr);
    CHKERR meshsets_manager_ptr->updateAllMeshsetsByEntitiesChildren(
        BitRefLevel().set());
 
    CHKERR ptr->ExRawPtr->mechanicalBC(
        BitRefLevel().set(FLIPPED_BIT + num_refinement_levels),
        BitRefLevel().set());
    CHKERR ptr->ExRawPtr->thermalBC(
        BitRefLevel().set(FLIPPED_BIT + num_refinement_levels),
        BitRefLevel().set());
 
    VecScatter scatter_to_final;
 
    for (PetscInt i = 0; i < nv; ++i) {
      CHKERR DMCreateGlobalVector(simple->getDM(), &ts_vecsout[i]);
      CHKERR scatter_mng->vecScatterCreate(
          ts_vecsout[i], simple->getProblemName(), RowColData::ROW, vec_out[i],
          "INTERMEDIATE_DM", RowColData::ROW, &scatter_to_final);
      CHKERR VecScatterBegin(scatter_to_final, vec_out[i], ts_vecsout[i],
                             INSERT_VALUES, SCATTER_REVERSE);
      CHKERR VecScatterEnd(scatter_to_final, vec_out[i], ts_vecsout[i],
                           INSERT_VALUES, SCATTER_REVERSE);
      CHKERR VecScatterDestroy(&scatter_to_final);
    }
 
    for (PetscInt i = 0; i < nv; ++i) {
      CHKERR VecDestroy(&vec_in[i]);
      CHKERR VecDestroy(&vec_out[i]);
    }
 
    auto set_all_bit = [&](EntityHandle ent, BitRefLevel &bit) {
      for (auto l = 0; l <= FLIPPED_BIT + num_refinement_levels; ++l) {
        bit.set(STORAGE_BIT, true);
      }
      if (bit[VIRGIN_BIT]) {
        bit.set(VIRGIN_BIT, false);
      }
      if (bit[FLIPPED_BIT]) {
        bit.set(VIRGIN_BIT, true);
      }
      if (bit[PREVIOUS_BIT]) {
        bit.set(PREVIOUS_BIT, false);
      }
      if (bit[FLIPPED_BIT + num_refinement_levels]) {
        bit.set(COMPUTATION_BIT, true);
      }
    };
    CHKERR bit_mng->lambdaBitRefLevel(set_all_bit);
 
#ifndef NDEBUG
 
    // comp_mesh_bit = BitRefLevel().set(COMPUTATION_BIT);
    // flipped_bit = BitRefLevel().set(FLIPPED_BIT);
    // refined_bit = BitRefLevel().set(FIRST_REF_BIT);
    // prev_mesh_bit = BitRefLevel().set(PREVIOUS_BIT);
    // virgin_mesh_bit = BitRefLevel().set(VIRGIN_BIT);
    // storage_bit = BitRefLevel().set(STORAGE_BIT);
 
    if (m_field.get_comm_rank() == 0) {
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          comp_mesh_bit, BitRefLevel().set(), 2,
          "after_mapping_comp_mesh_bit.vtk", "VTK", "");
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          flipped_bit, BitRefLevel().set(), 2, "after_mapping_flipped_bit.vtk",
          "VTK", "");
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          refined_bit, BitRefLevel().set(), 2, "after_mapping_refined_bit.vtk",
          "VTK", "");
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          prev_mesh_bit, BitRefLevel().set(), 2,
          "after_mapping_prev_mesh_bit.vtk", "VTK", "");
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          virgin_mesh_bit, BitRefLevel().set(), 2,
          "after_mapping_virgin_mesh_bit.vtk", "VTK", "");
    }
#endif // NDEBUG
 
    CHKERR m_field.modify_problem_ref_level_set_bit(ts_problem_name,
                                                    comp_mesh_bit);
    CHKERR m_field.modify_problem_mask_ref_level_set_bit(ts_problem_name,
                                                         BitRefLevel().set());
 
    CHKERR TSSetDM(ts, simple->getDM());
 
    auto B = createDMMatrix(simple->getDM());
    CHKERR TSSetIJacobian(ts, B, B, nullptr, nullptr);
  }
 
  for (PetscInt i = 0; i < nv; ++i) {
    double nrm;
    CHKERR VecNorm(ts_vecsout[i], NORM_2, &nrm);
    MOFEM_LOG("REMESHING", Sev::warning)
        << "After remeshing: ts_vecsout[" << i << "] norm = " << nrm;
  }
 
  PetscFunctionReturn(0);
}

◆ postProcessHere()

auto postProcessHere	(	MoFEM::Interface &	m_field,
		SmartPetscObj< DM > &	dm,
		std::string	output_name,
		int &	counter = `*(new int(0))`
	)

Examples: thermoplastic.cpp.

Definition at line 413 of file thermoplastic.cpp.

                                                                          {
  MoFEMFunctionBegin;
 
  auto create_post_process_elements = [&]() {
    MoFEMFunctionBegin;
    auto pp_fe = boost::make_shared<PostProcEle>(m_field);
 
    auto push_vol_post_proc_ops = [&](auto &pp_fe) {
      MoFEMFunctionBegin;
 
      auto &pip = pp_fe->getOpPtrVector();
 
      auto TAU_ptr = boost::make_shared<VectorDouble>();
      pip.push_back(new OpCalculateScalarFieldValues("TAU", TAU_ptr));
      auto T_ptr = boost::make_shared<VectorDouble>();
      pip.push_back(new OpCalculateScalarFieldValues("T", T_ptr));
 
      auto T_grad_ptr = boost::make_shared<MatrixDouble>();
      pip.push_back(
          new OpCalculateScalarFieldGradient<SPACE_DIM>("T", T_grad_ptr));
      auto U_ptr = boost::make_shared<MatrixDouble>();
      pip.push_back(new OpCalculateVectorFieldValues<SPACE_DIM>("U", U_ptr));
      auto FLUX_ptr = boost::make_shared<MatrixDouble>();
      pip.push_back(
          new OpCalculateHVecVectorField<3, SPACE_DIM>("FLUX", FLUX_ptr));
 
      auto EP_ptr = boost::make_shared<MatrixDouble>();
      pip.push_back(
          new OpCalculateTensor2SymmetricFieldValues<SPACE_DIM>("EP", EP_ptr));
 
      using OpPPMap = OpPostProcMapInMoab<SPACE_DIM, SPACE_DIM>;
 
      pip.push_back(
 
          new OpPPMap(
 
              pp_fe->getPostProcMesh(), pp_fe->getMapGaussPts(),
 
              {{"TAU", TAU_ptr}, {"T", T_ptr}},
 
              {{"U", U_ptr}, {"FLUX", FLUX_ptr}, {"T_GRAD", T_grad_ptr}},
 
              {},
 
              {{"EP", EP_ptr}}
 
              )
 
      );
 
      MoFEMFunctionReturn(0);
    };
 
    CHKERR push_vol_post_proc_ops(pp_fe);
 
    if (m_field.get_comm_size() == 1) {
 
      CHKERR DMoFEMLoopFiniteElements(dm, "dFE", pp_fe);
      CHKERR pp_fe->writeFile("out_" + std::to_string(counter) + "_" +
                              output_name + ".h5m");
    } else {
 
      CHKERR DMoFEMLoopFiniteElementsUpAndLowRank(dm, "dFE", pp_fe, 0,
                                                  m_field.get_comm_size());
      auto &post_proc_moab = pp_fe->getPostProcMesh();
      auto file_name = "out_" + std::to_string(counter) + "_" + output_name +
                       "_" + std::to_string(m_field.get_comm_rank()) + ".vtk";
      MOFEM_LOG("WORLD", Sev::inform)
          << "Writing file " << file_name << std::endl;
      CHKERR post_proc_moab.write_file(file_name.c_str(), "VTK");
    }
    counter++;
 
    MoFEMFunctionReturn(0);
  };
 
  CHKERR create_post_process_elements();
 
  MoFEMFunctionReturn(0);
}

◆ postProcessPETScHere()

auto postProcessPETScHere	(	MoFEM::Interface &	m_field,
		SmartPetscObj< DM > &	dm,
		Vec	sol,
		std::string	output_name
	)

Examples: thermoplastic.cpp.

Definition at line 495 of file thermoplastic.cpp.

                                                          {
  MoFEMFunctionBegin;
 
  auto smart_sol = SmartPetscObj<Vec>(sol, true);
 
  auto create_post_process_elements = [&]() {
    MoFEMFunctionBegin;
    auto pp_fe = boost::make_shared<PostProcEle>(m_field);
 
    auto push_vol_post_proc_ops = [&](auto &pp_fe) {
      MoFEMFunctionBegin;
 
      auto &pip = pp_fe->getOpPtrVector();
 
      auto TAU_ptr = boost::make_shared<VectorDouble>();
      pip.push_back(
          new OpCalculateScalarFieldValues("TAU", TAU_ptr, smart_sol));
      auto T_ptr = boost::make_shared<VectorDouble>();
      pip.push_back(new OpCalculateScalarFieldValues("T", T_ptr, smart_sol));
 
      auto U_ptr = boost::make_shared<MatrixDouble>();
      pip.push_back(
          new OpCalculateVectorFieldValues<SPACE_DIM>("U", U_ptr, smart_sol));
      auto FLUX_ptr = boost::make_shared<MatrixDouble>();
      pip.push_back(new OpCalculateHVecVectorField<3, SPACE_DIM>(
          "FLUX", FLUX_ptr, smart_sol));
 
      auto EP_ptr = boost::make_shared<MatrixDouble>();
      pip.push_back(new OpCalculateTensor2SymmetricFieldValues<SPACE_DIM>(
          "EP", EP_ptr, smart_sol));
 
      using OpPPMap = OpPostProcMapInMoab<SPACE_DIM, SPACE_DIM>;
 
      pip.push_back(
 
          new OpPPMap(
 
              pp_fe->getPostProcMesh(), pp_fe->getMapGaussPts(),
 
              {{"TAU", TAU_ptr}, {"T", T_ptr}},
 
              {{"U", U_ptr}, {"FLUX", FLUX_ptr}},
 
              {},
 
              {{"EP", EP_ptr}}
 
              )
 
      );
 
      MoFEMFunctionReturn(0);
    };
 
    CHKERR push_vol_post_proc_ops(pp_fe);
 
    CHKERR DMoFEMLoopFiniteElements(dm, "dFE", pp_fe);
    CHKERR pp_fe->writeFile("out_" + output_name + ".h5m");
 
    MoFEMFunctionReturn(0);
  };
 
  CHKERR create_post_process_elements();
 
  MoFEMFunctionReturn(0);
}

◆ printOnAllCores()

auto printOnAllCores	(	MoFEM::Interface &	m_field,
		const std::string &	out_put_string,
		const auto &	out_put_quantity
	)

Examples: thermoplastic.cpp.

Definition at line 563 of file thermoplastic.cpp.

                                                   {
  MoFEMFunctionBegin;
  auto rank = m_field.get_comm_rank();
  auto size = m_field.get_comm_size();
 
  // Loop through all processes and print one by one
  for (int i = 0; i < size; ++i) {
    // If it is this process's turn, print the data
    if (i == rank) {
      std::cout << "Proc " << rank << " " + out_put_string + " "
                << out_put_quantity << std::endl;
    }
    // All processes wait here until the current one has finished
    // printing
    CHKERR PetscBarrier(NULL);
  }
  MoFEMFunctionReturn(0);
};

◆ Qinf_thermal()

double Qinf_thermal	(	double	Qinf,
		double	omega_h,
		double	temp_0,
		double	temp
	)

inline

Examples: thermoplastic.cpp.

Definition at line 161 of file thermoplastic.cpp.

                                        {
  return Qinf * (1. - omega_h * (temp - temp_0));
}

◆ SNESIterationMonitor()

MoFEMErrorCode SNESIterationMonitor	(	SNES	snes,
		PetscInt	its,
		PetscReal	norm,
		void *	ctx
	)

Examples: thermoplastic.cpp.

Definition at line 4208 of file thermoplastic.cpp.

                                               {
  MoFEMFunctionBegin;
  MyTsCtx *my_ctx = static_cast<MyTsCtx *>(ctx);
  my_ctx->snes_iter_counter++;
  MoFEMFunctionReturn(0);
}

◆ temp_hat()

double temp_hat ( double temp )

inline

Examples: thermoplastic.cpp.

Definition at line 177 of file thermoplastic.cpp.

                                    {
  double JC_ref_temp = 0.;
  double JC_melt_temp = 1650.;
 
  return (temp - JC_ref_temp) / (JC_melt_temp - JC_ref_temp);
}

◆ TSIterationPreStage()

MoFEMErrorCode TSIterationPreStage	(	TS	ts,
		PetscReal	stagetime
	)

Examples: thermoplastic.cpp.

Definition at line 4200 of file thermoplastic.cpp.

                                                               {
  MoFEMFunctionBegin;
  MyTsCtx *my_ctx = ts_to_ctx[ts];
  my_ctx->snes_iter_counter = 0;
  MoFEMFunctionReturn(0);
}

◆ y_0_thermal()

double y_0_thermal	(	double	sigmaY,
		double	omega_0,
		double	temp_0,
		double	temp
	)

inline

Examples: thermoplastic.cpp.

Definition at line 152 of file thermoplastic.cpp.

                                       {
  return sigmaY * (1. - omega_0 * (temp - temp_0));
}

Variable Documentation

◆ alpha_damping

double alpha_damping = 0

Definition at line 315 of file thermoplastic.cpp.

◆ AT

constexpr AssemblyType AT

constexpr

Initial value:

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

Examples: thermoplastic.cpp.

Definition at line 51 of file thermoplastic.cpp.

◆ atom_test

int atom_test = 0

Definition at line 307 of file thermoplastic.cpp.

◆ b_iso

double b_iso = 63.69

Saturation exponent.

Definition at line 272 of file thermoplastic.cpp.

◆ C1_k

double C1_k = 0

Kinematic hardening.

Definition at line 273 of file thermoplastic.cpp.

◆ cn0

double cn0 = 1

Definition at line 275 of file thermoplastic.cpp.

◆ cn1

double cn1 = 1

Definition at line 276 of file thermoplastic.cpp.

◆ comp_mesh_bit

BitRefLevel comp_mesh_bit = BitRefLevel().set(COMPUTATION_BIT)

Examples: thermoplastic.cpp.

Definition at line 333 of file thermoplastic.cpp.

◆ CONTACT_SPACE

constexpr FieldSpace CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE

constexpr

Definition at line 59 of file thermoplastic.cpp.

◆ default_coeff_expansion

double default_coeff_expansion = 1e-5

Definition at line 287 of file thermoplastic.cpp.

◆ default_heat_capacity

double default_heat_capacity = 10*2.332

Definition at line 291 of file thermoplastic.cpp.

◆ default_heat_capacity_scale

double default_heat_capacity_scale = 1.

Examples: thermoplastic.cpp.

Definition at line 260 of file thermoplastic.cpp.

◆ default_heat_conductivity

double default_heat_conductivity

Initial value:

=

11.4

Definition at line 288 of file thermoplastic.cpp.

◆ default_ref_temp

double default_ref_temp = 20.0

Definition at line 278 of file thermoplastic.cpp.

◆ default_thermal_conductivity_scale

double default_thermal_conductivity_scale = 1.

Examples: ThermoElasticOps.hpp, and thermoplastic.cpp.

Definition at line 259 of file thermoplastic.cpp.

◆ do_eval_field

PetscBool do_eval_field = PETSC_FALSE

Evaluate field.

Definition at line 251 of file thermoplastic.cpp.

◆ edge_growth_thresh

double edge_growth_thresh = 0.2

Examples: thermoplastic.cpp.

Definition at line 320 of file thermoplastic.cpp.

◆ ep_order

int ep_order = order - 1

Order of ep field.

Definition at line 302 of file thermoplastic.cpp.

◆ exp_C

double exp_C = -6.284

Examples: thermoplastic.cpp.

Definition at line 174 of file thermoplastic.cpp.

◆ exp_D

double exp_D = -1.024

Examples: thermoplastic.cpp.

Definition at line 175 of file thermoplastic.cpp.

◆ flipped_bit

BitRefLevel flipped_bit = BitRefLevel().set(FLIPPED_BIT)

Examples: thermoplastic.cpp.

Definition at line 334 of file thermoplastic.cpp.

◆ flux_order

int flux_order = order

Order of tau field.

Examples: thermoplastic.cpp.

Definition at line 303 of file thermoplastic.cpp.

◆ Gaussian_distribution

auto Gaussian_distribution

Initial value:

= [](double init_temp, double peak_temp, double x,
                                double y, double z) {
  double s =
      width /
      std::sqrt(-2. * std::log(0.01 * peak_temp / (peak_temp - init_temp)));
  return (peak_temp - init_temp) * exp(-pow(y, 2) / (2. * pow(s, 2))) +
         init_temp;
}

Examples: thermoplastic.cpp.

Definition at line 340 of file thermoplastic.cpp.

                                                    {
  double s =
      width /
      std::sqrt(-2. * std::log(0.01 * peak_temp / (peak_temp - init_temp)));
  return (peak_temp - init_temp) * exp(-pow(y, 2) / (2. * pow(s, 2))) +
         init_temp;
};

◆ geom_order

int geom_order = 2

Order if fixed.

Examples: thermoplastic.cpp.

Definition at line 305 of file thermoplastic.cpp.

◆ get_string_from_vector

auto get_string_from_vector

Initial value:

= [](auto vec) {
  std::string output_string = "";
  for (auto el : vec) {
    output_string += std::to_string(el) + " ";
  }
  return output_string;
}

Examples: thermoplastic.cpp.

Definition at line 633 of file thermoplastic.cpp.

                                           {
  std::string output_string = "";
  for (auto el : vec) {
    output_string += std::to_string(el) + " ";
  }
  return output_string;
};

◆ H

double H = 584.3

Hardening.

Examples: thermoplastic.cpp.

Definition at line 268 of file thermoplastic.cpp.

◆ heat_capacity_scaling

ScalerFunTwoArgs heat_capacity_scaling

Examples: ThermoElasticOps.hpp, ThermoPlasticOps.hpp, and thermoplastic.cpp.

Definition at line 262 of file thermoplastic.cpp.

◆ help

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

static

[Solve]

Definition at line 5741 of file thermoplastic.cpp.

◆ ic_type

ICType ic_type = IC_UNIFORM

Examples: thermoplastic.cpp.

Definition at line 283 of file thermoplastic.cpp.

◆ ICTypes

const char* const ICTypes[]

Initial value:

= {"uniform", "gaussian", "linear",

"ICType", "IC_", nullptr}

Examples: thermoplastic.cpp.

Definition at line 280 of file thermoplastic.cpp.

280 {"uniform", "gaussian", "linear",

281 "ICType", "IC_", nullptr};

◆ inelastic_heat_fraction

double inelastic_heat_fraction

Initial value:

=

0.9

fraction of plastic dissipation converted to heat

Examples: thermoplastic.cpp.

Definition at line 296 of file thermoplastic.cpp.

◆ inelastic_heat_fraction_scaling

ScalerFunThreeArgs inelastic_heat_fraction_scaling

Examples: ThermoPlasticOps.hpp, and thermoplastic.cpp.

Definition at line 263 of file thermoplastic.cpp.

◆ init_dt

double init_dt = 0.05

Examples: thermoplastic.cpp.

Definition at line 316 of file thermoplastic.cpp.

◆ init_T

auto init_T

Initial value:

= [](double init_temp, double peak_temp, double x, double y,
                 double z) {
  switch (ic_type) {
  case IC_GAUSSIAN:
    return -Gaussian_distribution(init_temp, peak_temp, x, y, z);
  case IC_LINEAR:
    return -linear_distribution(init_temp, peak_temp, x, y, z);
  case IC_UNIFORM:
  default:
    return -uniform_distribution(init_temp, peak_temp, x, y, z);
  }
}

Initialisation function for temperature field.

Examples: EdgeFlippingOps.hpp, ThermoElasticOps.hpp, and thermoplastic.cpp.

Definition at line 360 of file thermoplastic.cpp.

                           {
  switch (ic_type) {
  case IC_GAUSSIAN:
    return -Gaussian_distribution(init_temp, peak_temp, x, y, z);
  case IC_LINEAR:
    return -linear_distribution(init_temp, peak_temp, x, y, z);
  case IC_UNIFORM:
  default:
    return -uniform_distribution(init_temp, peak_temp, x, y, z);
  }
};

◆ init_temp

double init_temp = 20.0

Examples: thermoplastic.cpp.

Definition at line 284 of file thermoplastic.cpp.

◆ is_distributed_mesh

PetscBool is_distributed_mesh = PETSC_TRUE

Examples: SolutionMapping.hpp, and thermoplastic.cpp.

Definition at line 256 of file thermoplastic.cpp.

◆ IS_LARGE_STRAINS

constexpr bool IS_LARGE_STRAINS

constexpr

Initial value:

=

FINITE_DEFORMATION_FLAG

#define FINITE_DEFORMATION_FLAG

Definition thermo_elastic.cpp:14

Definition at line 48 of file thermoplastic.cpp.

◆ is_large_strains

const bool is_large_strains

Initial value:

=

IS_LARGE_STRAINS

constexpr bool IS_LARGE_STRAINS

Definition thermoplastic.cpp:48

Definition at line 248 of file thermoplastic.cpp.

◆ is_quasi_static

PetscBool is_quasi_static = PETSC_TRUE

Definition at line 313 of file thermoplastic.cpp.

◆ IT

constexpr IntegrationType IT

constexpr

Initial value:

=

IntegrationType::GAUSS

Examples: thermoplastic.cpp.

Definition at line 54 of file thermoplastic.cpp.

◆ linear_distribution

auto linear_distribution

Initial value:

= [](double init_temp, double peak_temp, double x,
                              double y, double z) {
  return ((peak_temp - init_temp) / width) * y + (peak_temp + init_temp) / 2.0;
}

Examples: thermoplastic.cpp.

Definition at line 349 of file thermoplastic.cpp.

                                                  {
  return ((peak_temp - init_temp) / width) * y + (peak_temp + init_temp) / 2.0;
};

◆ min_dt

double min_dt = 1e-12

Examples: thermoplastic.cpp.

Definition at line 317 of file thermoplastic.cpp.

◆ num_refinement_levels

int num_refinement_levels = 0

Examples: SolutionMapping.hpp, and thermoplastic.cpp.

Definition at line 321 of file thermoplastic.cpp.

◆ omega_0

double omega_0 = 2e-3

flow stress softening

Examples: thermoplastic.cpp.

Definition at line 294 of file thermoplastic.cpp.

◆ omega_h

double omega_h = 2e-3

hardening softening

Examples: thermoplastic.cpp.

Definition at line 295 of file thermoplastic.cpp.

◆ order

int order = 2

Order displacement.

Definition at line 300 of file thermoplastic.cpp.

◆ order_edge

PetscBool order_edge = PETSC_FALSE

Examples: mofem/tutorials/adv-0_plasticity/plastic.cpp, plastic.cpp, and thermoplastic.cpp.

Definition at line 309 of file thermoplastic.cpp.

◆ order_face

PetscBool order_face = PETSC_FALSE

Examples: mofem/tutorials/adv-0_plasticity/plastic.cpp, plastic.cpp, and thermoplastic.cpp.

Definition at line 310 of file thermoplastic.cpp.

◆ order_volume

PetscBool order_volume = PETSC_FALSE

Examples: mofem/tutorials/adv-0_plasticity/plastic.cpp, plastic.cpp, and thermoplastic.cpp.

Definition at line 311 of file thermoplastic.cpp.

◆ peak_temp

double peak_temp = 1000.0

Examples: thermoplastic.cpp.

Definition at line 285 of file thermoplastic.cpp.

◆ poisson_ratio

double poisson_ratio = 0.31

Poisson ratio.

Definition at line 266 of file thermoplastic.cpp.

◆ post_processing_counter

int post_processing_counter = 0

Examples: thermoplastic.cpp.

Definition at line 411 of file thermoplastic.cpp.

◆ prev_mesh_bit

BitRefLevel prev_mesh_bit = BitRefLevel().set(PREVIOUS_BIT)

Examples: thermoplastic.cpp.

Definition at line 336 of file thermoplastic.cpp.

◆ Qinf

double Qinf = 174.2

Saturation yield stress.

Definition at line 271 of file thermoplastic.cpp.

◆ qual_thresh

double qual_thresh = 0.1

Examples: thermoplastic.cpp.

Definition at line 319 of file thermoplastic.cpp.

◆ qual_tol

double qual_tol = 0

Examples: thermoplastic.cpp.

Definition at line 318 of file thermoplastic.cpp.

◆ refined_bit

BitRefLevel refined_bit = BitRefLevel().set(FIRST_REF_BIT)

Examples: thermoplastic.cpp.

Definition at line 335 of file thermoplastic.cpp.

◆ restart_file

char restart_file[255]

Examples: ep.cpp, mofem/users_modules/eshelbian_plasticity/ep.cpp, and thermoplastic.cpp.

Definition at line 252 of file thermoplastic.cpp.

◆ restart_flg

PetscBool restart_flg = PETSC_TRUE

Examples: ep.cpp, mofem/users_modules/eshelbian_plasticity/ep.cpp, and thermoplastic.cpp.

Definition at line 253 of file thermoplastic.cpp.

◆ restart_step

int restart_step = 0

Examples: PlasticOpsMonitor.hpp, and thermoplastic.cpp.

Definition at line 255 of file thermoplastic.cpp.

◆ restart_time

double restart_time = 0

Examples: thermoplastic.cpp.

Definition at line 254 of file thermoplastic.cpp.

◆ rho

double rho = 0.0

Definition at line 314 of file thermoplastic.cpp.

◆ save_range

auto save_range

Initial value:

= [](moab::Interface &moab, const std::string name,
                     const Range r) {
  MoFEMFunctionBegin;
  auto out_meshset = get_temp_meshset_ptr(moab);
  CHKERR moab.add_entities(*out_meshset, r);
  CHKERR moab.write_file(name.c_str(), "VTK", "", out_meshset->get_ptr(), 1);
  MoFEMFunctionReturn(0);
}

Examples: thermoplastic.cpp.

Definition at line 624 of file thermoplastic.cpp.

                                    {
  MoFEMFunctionBegin;
  auto out_meshset = get_temp_meshset_ptr(moab);
  CHKERR moab.add_entities(*out_meshset, r);
  CHKERR moab.write_file(name.c_str(), "VTK", "", out_meshset->get_ptr(), 1);
  MoFEMFunctionReturn(0);
};

◆ scale

double scale = 1.

Examples: thermoplastic.cpp.

Definition at line 258 of file thermoplastic.cpp.

◆ set_timer

PetscBool set_timer = PETSC_FALSE

Set timer.

Definition at line 250 of file thermoplastic.cpp.

◆ sigmaY

double sigmaY = 936.2

Yield stress.

Definition at line 267 of file thermoplastic.cpp.

◆ size_symm

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

constexpr

Definition at line 46 of file thermoplastic.cpp.

◆ SPACE_DIM

constexpr int SPACE_DIM

constexpr

Initial value:

=

3

Examples: thermoplastic.cpp.

Definition at line 44 of file thermoplastic.cpp.

◆ storage_bit

BitRefLevel storage_bit = BitRefLevel().set(STORAGE_BIT)

Examples: thermoplastic.cpp.

Definition at line 338 of file thermoplastic.cpp.

◆ T_order

int T_order = order - 1

Order of ep field.

Examples: thermoplastic.cpp.

Definition at line 304 of file thermoplastic.cpp.

◆ tau_order

int tau_order = order - 2

Order of tau field.

Definition at line 301 of file thermoplastic.cpp.

◆ temp_0

double temp_0 = 20

reference temperature for thermal softening

Examples: thermoplastic.cpp.

Definition at line 298 of file thermoplastic.cpp.

◆ thermal_conductivity_scaling

ScalerFunTwoArgs thermal_conductivity_scaling

Examples: ThermoElasticOps.hpp, ThermoPlasticOps.hpp, and thermoplastic.cpp.

Definition at line 261 of file thermoplastic.cpp.

◆ thermoplastic_raw_ptr

Example* thermoplastic_raw_ptr = nullptr

Examples: thermoplastic.cpp.

Definition at line 4472 of file thermoplastic.cpp.

◆ ts_to_ctx

std::unordered_map<TS, MyTsCtx *> ts_to_ctx

static

Examples: thermoplastic.cpp.

Definition at line 4198 of file thermoplastic.cpp.

◆ ts_to_rctx

std::unordered_map<TS, ResizeCtx *> ts_to_rctx

static

Examples: thermoplastic.cpp.

Definition at line 703 of file thermoplastic.cpp.

◆ tsPrePostProc

boost::weak_ptr<TSPrePostProc> tsPrePostProc

static

Examples: mofem/tutorials/adv-4_dynamic_first_order_con_law/dynamic_first_order_con_law.cpp, mofem/tutorials/vec-5_free_surface/free_surface.cpp, and thermoplastic.cpp.

Definition at line 689 of file thermoplastic.cpp.

◆ uniform_distribution

auto uniform_distribution

Initial value:

= [](double init_temp, double peak_temp, double x,

double y, double z) { return init_temp; }

Examples: thermoplastic.cpp.

Definition at line 354 of file thermoplastic.cpp.

355 { return init_temp; };

◆ virgin_mesh_bit

BitRefLevel virgin_mesh_bit = BitRefLevel().set(VIRGIN_BIT)

Examples: thermoplastic.cpp.

Definition at line 337 of file thermoplastic.cpp.

◆ visH

double visH = 0

Viscous hardening.

Definition at line 269 of file thermoplastic.cpp.

◆ width

double width = 10.0

Width of Gaussian distribution.

Examples: thermoplastic.cpp.

Definition at line 286 of file thermoplastic.cpp.

◆ young_modulus

double young_modulus = 115000

Young modulus.

Definition at line 265 of file thermoplastic.cpp.

◆ zeta

double zeta = 5e-2

regularisation parameter

Definition at line 270 of file thermoplastic.cpp.

Classes

Macros

Typedefs

Enumerations

Functions

Variables

Macro Definition Documentation

◆ EXECUTABLE_DIMENSION

Typedef Documentation

◆ AssemblyDomainEleOp

◆ BoundaryEle

◆ BoundaryEleOp

◆ BoundaryLhsBCs

◆ BoundaryNaturalBC

◆ BoundaryRhsBCs

◆ DomainEle

◆ DomainEleOp

◆ DomainNaturalBCLhs

◆ DomainNaturalBCRhs

◆ DomainRhsBCs

◆ EntData

◆ OpBodyForce

◆ OpBoundaryLhsBCs

◆ OpBoundaryRhsBCs

◆ OpDomainRhsBCs

◆ OpEssentialFluxLhs

◆ OpEssentialFluxRhs

◆ OpForce

◆ OpHeatSource

◆ OpSetTemperatureLhs

◆ OpSetTemperatureRhs

◆ OpTemperatureBC

◆ PostProcEle

◆ ScalerFunThreeArgs

◆ ScalerFunTwoArgs

◆ SetPtsData

◆ SideEle

◆ SkinPostProcEle

Enumeration Type Documentation

◆ ICType

◆ ThermoPlasticityBits

Function Documentation

◆ dH_thermal_dtemp()

◆ dQinf_thermal_dtemp()

◆ dy_0_thermal_dtemp()

◆ exp_softening()

◆ H_thermal()

◆ iso_hardening()

◆ iso_hardening_dtau()

◆ iso_hardening_dtemp()

◆ iso_hardening_exp()

◆ kinematic_hardening()

◆ kinematic_hardening_dplastic_strain()

◆ main()

◆ MyTSResizeSetup()

◆ MyTSResizeTransfer()

◆ postProcessHere()

◆ postProcessPETScHere()

◆ printOnAllCores()

◆ Qinf_thermal()

◆ SNESIterationMonitor()

◆ temp_hat()

◆ TSIterationPreStage()

◆ y_0_thermal()

Variable Documentation

◆ alpha_damping

◆ AT

◆ atom_test

◆ b_iso

◆ C1_k

◆ cn0

◆ cn1

◆ comp_mesh_bit

◆ CONTACT_SPACE

◆ default_coeff_expansion

◆ default_heat_capacity

◆ default_heat_capacity_scale

◆ default_heat_conductivity

◆ default_ref_temp

◆ default_thermal_conductivity_scale