thermoplastic.cpp File Reference

#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>

Go to the source code of this file.

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	DomainEleOnRange = ElementsAndOps< SPACE_DIM >::DomainEleOnRange
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 }

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)
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)
PetscErrorCode	MyTSResizeTransfer (TS ts, PetscInt nv, Vec ts_vecsin[], Vec ts_vecsout[], 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 = 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
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
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 76 of file thermoplastic.cpp.

BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 69 of file thermoplastic.cpp.

BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 70 of file thermoplastic.cpp.

BoundaryLhsBCs

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

Definition at line 510 of file thermoplastic.cpp.

BoundaryNaturalBC

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

[Body and heat source]

[Natural boundary conditions]

Definition at line 525 of file thermoplastic.cpp.

BoundaryRhsBCs

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

Definition at line 507 of file thermoplastic.cpp.

DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 66 of file thermoplastic.cpp.

DomainEleOnRange

using DomainEleOnRange = ElementsAndOps<SPACE_DIM>::DomainEleOnRange

Definition at line 67 of file thermoplastic.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 68 of file thermoplastic.cpp.

DomainNaturalBCLhs

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

Definition at line 520 of file thermoplastic.cpp.

DomainNaturalBCRhs

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

[Body and heat source]

Definition at line 515 of file thermoplastic.cpp.

DomainRhsBCs

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

Definition at line 504 of file thermoplastic.cpp.

EntData

using EntData = EntitiesFieldData::EntData

Definition at line 65 of file thermoplastic.cpp.

OpBodyForce

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

Definition at line 516 of file thermoplastic.cpp.

OpBoundaryLhsBCs

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

Definition at line 511 of file thermoplastic.cpp.

OpBoundaryRhsBCs

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

Definition at line 508 of file thermoplastic.cpp.

OpDomainRhsBCs

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

Definition at line 505 of file thermoplastic.cpp.

OpEssentialFluxLhs

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

Definition at line 535 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 533 of file thermoplastic.cpp.

OpForce

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

Definition at line 527 of file thermoplastic.cpp.

OpHeatSource

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

Definition at line 518 of file thermoplastic.cpp.

OpSetTemperatureLhs

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

Definition at line 541 of file thermoplastic.cpp.

OpSetTemperatureRhs

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

[Essential boundary conditions (Least square approach)]

Definition at line 539 of file thermoplastic.cpp.

OpTemperatureBC

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

Definition at line 528 of file thermoplastic.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 71 of file thermoplastic.cpp.

ScalerFunThreeArgs

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

Definition at line 79 of file thermoplastic.cpp.

ScalerFunTwoArgs

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

Definition at line 78 of file thermoplastic.cpp.

SetPtsData

using SetPtsData = FieldEvaluatorInterface::SetPtsData

Definition at line 74 of file thermoplastic.cpp.

SideEle

using SideEle = ElementsAndOps<SPACE_DIM>::SideEle

Definition at line 73 of file thermoplastic.cpp.

SkinPostProcEle

using SkinPostProcEle = PostProcBrokenMeshInMoab<BoundaryEle>

Definition at line 72 of file thermoplastic.cpp.

Enumeration Type Documentation

ICType

enum ICType

Enumerator
IC_UNIFORM
IC_GAUSSIAN
IC_LINEAR

Definition at line 217 of file thermoplastic.cpp.

{ IC_UNIFORM, IC_GAUSSIAN, IC_LINEAR };

Function Documentation

dH_thermal_dtemp()

double dH_thermal_dtemp ( double  H, double  omega_h  )

inline

Examples

thermoplastic.cpp.

Definition at line 86 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 104 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 95 of file thermoplastic.cpp.

                                                                {
  return -sigmaY * omega_0;
}

exp_softening()

double exp_softening ( double  temp_hat )

inline

Examples

thermoplastic.cpp.

Definition at line 122 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 82 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 129 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 136 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 143 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 108 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 161 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 175 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]

[Load mesh]

[Example]

[Example]

Definition at line 5339 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();
    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);
    }
    //! [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 3858 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
	)

Examples

thermoplastic.cpp.

Definition at line 3866 of file thermoplastic.cpp.

                                                               {
  PetscFunctionBegin;
  ResizeCtx *rctx = static_cast<ResizeCtx *>(ctx);
 
  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;
  }
 
  if (auto ptr = tsPrePostProc.lock()) {
 
    MoFEM::Interface &m_field = *(rctx->m_field);
 
    MOFEM_LOG("REMESHING", Sev::verbose)
        << "rctx->all_old_els = " << rctx->all_old_els;
    MOFEM_LOG("REMESHING", Sev::verbose)
        << "rctx->new_elements = " << rctx->new_elements;
    MOFEM_LOG("REMESHING", Sev::verbose)
        << "rctx->flipped_els = " << rctx->flipped_els;
 
    auto scatter_mng = m_field.getInterface<VecManager>();
    auto simple = m_field.getInterface<Simple>();
 
    // Rebuilding DM with old and new mesh entities
    auto bit = BitRefLevel().set();
    CHKERR SolutionMapping(m_field).reSetUp(
        {"T", "TAU", "EP"}, {T_order, tau_order, ep_order}, bit);
 
    auto intermediate_dm = createDM(m_field.get_comm(), "DMMOFEM");
    CHKERR DMMoFEMCreateMoFEM(intermediate_dm, &m_field, "INTERMEDIATE_DM",
                              BitRefLevel().set(),
                              BitRefLevel().set(BITREFLEVEL_SIZE - 1).flip());
    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 DMMoFEMSetIsPartitioned(intermediate_dm, is_distributed_mesh);
    CHKERR DMSetUp_MoFEM(intermediate_dm);
 
    constexpr bool do_fake_mapping = false;
    if (do_fake_mapping) {
      auto tmp_D = createDMVector(intermediate_dm);
      CHKERR DMoFEMMeshToLocalVector(intermediate_dm, tmp_D, INSERT_VALUES,
                                     SCATTER_FORWARD);
      CHKERR VecCopy(vec_in[0], vec_out[0]);
 
      auto sub_dm = createDM(m_field.get_comm(), "DMMOFEM");
      CHKERR DMSetType(sub_dm, "DMMOFEM");
      CHKERR DMMoFEMCreateSubDM(sub_dm, intermediate_dm, "SUB_MAPPING");
      CHKERR DMMoFEMAddElement(sub_dm, simple->getDomainFEName());
      CHKERR DMMoFEMSetSquareProblem(sub_dm, PETSC_TRUE);
      CHKERR DMMoFEMSetDestroyProblem(sub_dm, PETSC_TRUE);
 
      MOFEM_LOG("REMESHING", Sev::inform) << "Created sub DpoM";
 
      for (auto f : {"T", "TAU", "EP"}) {
        CHKERR DMMoFEMAddSubFieldRow(
            sub_dm, f, boost::make_shared<Range>(rctx->new_elements));
        CHKERR DMMoFEMAddSubFieldCol(
            sub_dm, f, boost::make_shared<Range>(rctx->new_elements));
      }
 
      auto bit = BitRefLevel().set();
      CHKERR m_field.modify_problem_ref_level_set_bit("SUB_MAPPING", bit);
      CHKERR m_field.modify_problem_mask_ref_level_set_bit("SUB_MAPPING", bit);
 
      CHKERR DMMoFEMSetIsPartitioned(sub_dm, is_distributed_mesh);
      CHKERR DMSubDMSetUp_MoFEM(sub_dm);
 
      auto fe_method = boost::make_shared<DomainEle>(m_field);
      CHKERR DMoFEMLoopFiniteElements(sub_dm, simple->getDomainFEName(),
                                      fe_method);
 
      CHKERR DMoFEMMeshToLocalVector(intermediate_dm, tmp_D, INSERT_VALUES,
                                     SCATTER_REVERSE);
    } else {
      CHKERR SolutionMapping(m_field).mapFields(
          intermediate_dm, rctx->all_old_els, rctx->new_elements,
          rctx->flipped_els, 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(1));
    CHKERR m_field.modify_problem_mask_ref_level_set_bit(
        ts_problem_name, BitRefLevel().set(BITREFLEVEL_SIZE - 1).flip());
    CHKERR DMMoFEMSetIsPartitioned(simple->getDM(), is_distributed_mesh);
    CHKERR DMSetUp_MoFEM(simple->getDM());
 
#ifndef NDEBUG
 
    // e10 -> e11 (old to new) no flipped
    // e10 -> e10 old flipped
    // -> e01 new flliped
    // e11, e01 -> bit01 mask 11 (all)
 
    if (m_field.get_comm_rank() == 0) {
      MOFEM_LOG("REMESHING", Sev::verbose)
          << "Writing debug VTK files for remeshing verification";
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          BitRefLevel().set(0), BitRefLevel().set(0), 2, "bit_0_mask_0.vtk",
          "VTK", ""); // that are elements on old which are flipped
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          BitRefLevel().set(0), BitRefLevel().set(), 2, "bit_0_mask_all.vtk",
          "VTK", ""); // all old elements
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          BitRefLevel().set(1), BitRefLevel().set(1), 2, "bit_1_mask_1.vtk",
          "VTK", ""); // that are new elements
      CHKERR m_field.getInterface<BitRefManager>()->writeBitLevelByDim(
          BitRefLevel().set(1), BitRefLevel().set(), 2, "bit_1_mask_all.vtk",
          "VTK", ""); // all new elements
    }
#endif // NDEBUG
 
    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]);
    }
 
    // Squash and change bit ref level back
    auto bit_mng = m_field.getInterface<BitRefManager>();
 
    Range flipped_ents;
    CHKERR bit_mng->getEntitiesByRefLevel(BitRefLevel().set(0),
                                          BitRefLevel().set(0), flipped_ents);
    CHKERR bit_mng->addBitRefLevel(flipped_ents,
                                   BitRefLevel().set(BITREFLEVEL_SIZE - 1));
    BitRefLevel shift_mask = BitRefLevel().set(BITREFLEVEL_SIZE - 1);
    shift_mask.flip();
    CHKERR bit_mng->shiftRightBitRef(1, shift_mask);
    CHKERR m_field.modify_problem_ref_level_set_bit(ts_problem_name,
                                                    BitRefLevel().set(0));
    // CHKERR DMMoFEMSetIsPartitioned(simple->getDM(), is_distributed_mesh);
    // CHKERR DMSetUp_MoFEM(simple->getDM());
 
    CHKERR TSSetDM(ts, simple->getDM());
    auto B = createDMMatrix(simple->getDM());
    CHKERR TSSetIJacobian(ts, B, B, nullptr, nullptr);
    rctx->mesh_changed = PETSC_FALSE;
  }
 
  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 331 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 413 of file thermoplastic.cpp.

                                                          {
  MoFEMFunctionBegin;
 
  auto smart_sol = SmartPetscObj<Vec>(sol, true);
 
  auto create_post_process_elements = [&]() {
    MoFEMFunctionBegin;
    Simple *simple = m_field.getInterface<Simple>();
    auto pp_fe = boost::make_shared<PostProcEle>(m_field);
    auto &pip = pp_fe->getOpPtrVector();
 
    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 483 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 99 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 3850 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 115 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 3842 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 90 of file thermoplastic.cpp.

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

Variable Documentation

alpha_damping

double alpha_damping = 0

Definition at line 253 of file thermoplastic.cpp.

AT

constexpr AssemblyType AT

constexpr

Initial value:

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

Examples

thermoplastic.cpp.

Definition at line 55 of file thermoplastic.cpp.

atom_test

int atom_test = 0

Definition at line 245 of file thermoplastic.cpp.

b_iso

double b_iso = 63.69

Saturation exponent.

Definition at line 210 of file thermoplastic.cpp.

C1_k

double C1_k = 0

Kinematic hardening.

Definition at line 211 of file thermoplastic.cpp.

cn0

double cn0 = 1

Definition at line 213 of file thermoplastic.cpp.

cn1

double cn1 = 1

Definition at line 214 of file thermoplastic.cpp.

CONTACT_SPACE

constexpr FieldSpace CONTACT_SPACE = ElementsAndOps<SPACE_DIM>::CONTACT_SPACE

constexpr

Definition at line 63 of file thermoplastic.cpp.

default_coeff_expansion

double default_coeff_expansion = 1e-5

Definition at line 225 of file thermoplastic.cpp.

default_heat_capacity

double default_heat_capacity = 2.332

Definition at line 229 of file thermoplastic.cpp.

default_heat_capacity_scale

double default_heat_capacity_scale = 1.

Examples

thermoplastic.cpp.

Definition at line 198 of file thermoplastic.cpp.

default_heat_conductivity

double default_heat_conductivity

Initial value:

=
    11.4

Definition at line 226 of file thermoplastic.cpp.

default_ref_temp

double default_ref_temp = 20.0

Definition at line 216 of file thermoplastic.cpp.

default_thermal_conductivity_scale

double default_thermal_conductivity_scale = 1.

Examples

ThermoElasticOps.hpp, and thermoplastic.cpp.

Definition at line 197 of file thermoplastic.cpp.

do_eval_field

PetscBool do_eval_field = PETSC_FALSE

Evaluate field.

Definition at line 189 of file thermoplastic.cpp.

ep_order

int ep_order = order - 1

Order of ep field.

Definition at line 240 of file thermoplastic.cpp.

exp_C

double exp_C = -6.284

Examples

thermoplastic.cpp.

Definition at line 112 of file thermoplastic.cpp.

exp_D

double exp_D = -1.024

Examples

thermoplastic.cpp.

Definition at line 113 of file thermoplastic.cpp.

flux_order

int flux_order = order

Order of tau field.

Examples

thermoplastic.cpp.

Definition at line 241 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 258 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 243 of file thermoplastic.cpp.

H

double H = 584.3

Hardening.

Examples

thermoplastic.cpp.

Definition at line 206 of file thermoplastic.cpp.

heat_capacity_scaling

ScalerFunTwoArgs heat_capacity_scaling

Examples

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

Definition at line 200 of file thermoplastic.cpp.

help

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

static

[Solve]

Definition at line 5337 of file thermoplastic.cpp.

ic_type

ICType ic_type = IC_UNIFORM

Examples

thermoplastic.cpp.

Definition at line 221 of file thermoplastic.cpp.

ICTypes

const char* const ICTypes[]

Initial value:

= {"uniform", "gaussian", "linear",
                               "ICType",  "IC_",      nullptr}

Examples

thermoplastic.cpp.

Definition at line 218 of file thermoplastic.cpp.

                              {"uniform", "gaussian", "linear",
                               "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 234 of file thermoplastic.cpp.

inelastic_heat_fraction_scaling

ScalerFunThreeArgs inelastic_heat_fraction_scaling

Examples

ThermoPlasticOps.hpp, and thermoplastic.cpp.

Definition at line 201 of file thermoplastic.cpp.

init_dt

double init_dt = 0.05

Examples

thermoplastic.cpp.

Definition at line 254 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 278 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 222 of file thermoplastic.cpp.

is_distributed_mesh

PetscBool is_distributed_mesh = PETSC_TRUE

Examples

SolutionMapping.hpp, and thermoplastic.cpp.

Definition at line 194 of file thermoplastic.cpp.

IS_LARGE_STRAINS

constexpr bool IS_LARGE_STRAINS

constexpr

Initial value:

=
    FINITE_DEFORMATION_FLAG

Definition at line 52 of file thermoplastic.cpp.

is_large_strains

const bool is_large_strains

Initial value:

=
    IS_LARGE_STRAINS

Definition at line 186 of file thermoplastic.cpp.

is_quasi_static

PetscBool is_quasi_static = PETSC_TRUE

Definition at line 251 of file thermoplastic.cpp.

IT

constexpr IntegrationType IT

constexpr

Initial value:

=
    IntegrationType::GAUSS

Examples

thermoplastic.cpp.

Definition at line 58 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 267 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 255 of file thermoplastic.cpp.

omega_0

double omega_0 = 2e-3

flow stress softening

Examples

thermoplastic.cpp.

Definition at line 232 of file thermoplastic.cpp.

omega_h

double omega_h = 2e-3

hardening softening

Examples

thermoplastic.cpp.

Definition at line 233 of file thermoplastic.cpp.

order

int order = 2

Order displacement.

Definition at line 238 of file thermoplastic.cpp.

order_edge

PetscBool order_edge = PETSC_FALSE

Examples

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

Definition at line 247 of file thermoplastic.cpp.

order_face

PetscBool order_face = PETSC_FALSE

Examples

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

Definition at line 248 of file thermoplastic.cpp.

order_volume

PetscBool order_volume = PETSC_FALSE

Examples

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

Definition at line 249 of file thermoplastic.cpp.

peak_temp

double peak_temp = 1000.0

Examples

thermoplastic.cpp.

Definition at line 223 of file thermoplastic.cpp.

poisson_ratio

double poisson_ratio = 0.31

Poisson ratio.

Definition at line 204 of file thermoplastic.cpp.

post_processing_counter

int post_processing_counter = 0

Examples

thermoplastic.cpp.

Definition at line 329 of file thermoplastic.cpp.

Qinf

double Qinf = 174.2

Saturation yield stress.

Definition at line 209 of file thermoplastic.cpp.

qual_tol

double qual_tol = 0

Examples

thermoplastic.cpp.

Definition at line 256 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 190 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 191 of file thermoplastic.cpp.

restart_step

int restart_step = 0

Examples

PlasticOpsMonitor.hpp, and thermoplastic.cpp.

Definition at line 193 of file thermoplastic.cpp.

restart_time

double restart_time = 0

Examples

thermoplastic.cpp.

Definition at line 192 of file thermoplastic.cpp.

rho

double rho = 0.0

Definition at line 252 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 544 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 196 of file thermoplastic.cpp.

set_timer

PetscBool set_timer = PETSC_FALSE

Set timer.

Definition at line 188 of file thermoplastic.cpp.

sigmaY

double sigmaY = 936.2

Yield stress.

Definition at line 205 of file thermoplastic.cpp.

size_symm

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

constexpr

Definition at line 50 of file thermoplastic.cpp.

SPACE_DIM

constexpr int SPACE_DIM

constexpr

Initial value:

=
     3

Examples

thermoplastic.cpp.

Definition at line 48 of file thermoplastic.cpp.

T_order

int T_order = order - 1

Order of ep field.

Examples

thermoplastic.cpp.

Definition at line 242 of file thermoplastic.cpp.

tau_order

int tau_order = order - 2

Order of tau field.

Definition at line 239 of file thermoplastic.cpp.

temp_0

double temp_0 = 20

reference temperature for thermal softening

Examples

thermoplastic.cpp.

Definition at line 236 of file thermoplastic.cpp.

thermal_conductivity_scaling

ScalerFunTwoArgs thermal_conductivity_scaling

Examples

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

Definition at line 199 of file thermoplastic.cpp.

thermoplastic_raw_ptr

Example* thermoplastic_raw_ptr = nullptr

Examples

thermoplastic.cpp.

Definition at line 4060 of file thermoplastic.cpp.

ts_to_ctx

std::unordered_map<TS, MyTsCtx *> ts_to_ctx

static

Examples

thermoplastic.cpp.

Definition at line 3840 of file thermoplastic.cpp.

ts_to_rctx

std::unordered_map<TS, ResizeCtx *> ts_to_rctx

static

Examples

thermoplastic.cpp.

Definition at line 613 of file thermoplastic.cpp.

tsPrePostProc

boost::weak_ptr<TSPrePostProc> tsPrePostProc

static

Examples

mofem/tutorials/adv-4/dynamic_first_order_con_law.cpp, mofem/tutorials/vec-5/free_surface.cpp, and thermoplastic.cpp.

Definition at line 601 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 272 of file thermoplastic.cpp.

                                                   { return init_temp; };

visH

double visH = 0

Viscous hardening.

Definition at line 207 of file thermoplastic.cpp.

width

double width = 10.0

Width of Gaussian distribution.

Examples

thermoplastic.cpp.

Definition at line 224 of file thermoplastic.cpp.

young_modulus

double young_modulus = 115000

Young modulus.

Definition at line 203 of file thermoplastic.cpp.

zeta

double zeta = 5e-2

regularisation parameter

Definition at line 208 of file thermoplastic.cpp.