#include <users_modules/linear_isotropic_hardening/src/PlDataAtGaussPts.hpp>

Collaboration diagram for CommonData:

Public Member Functions
	CommonData (MoFEM::Interface &m_field)

MoFEMErrorCode	getTags (const EntityHandle fe_ent, const int nb_gauss_pts, boost::shared_ptr< MatrixDouble > tag_ptr, Tag tag_ref)

MoFEMErrorCode	getTags (const EntityHandle fe_ent, const int nb_gauss_pts, boost::shared_ptr< VectorDouble > tag_ptr, Tag tag_ref)

template<class T >
MoFEMErrorCode	setTags (const EntityHandle fe_ent, const int nb_gauss_pts, T tag_ptr, Tag tag_ref)

MoFEMErrorCode	getParameters ()
	Get the Parameters object. More...

MoFEMErrorCode	calcReactionForces (VectorDouble &reactions, DM &dm, boost::shared_ptr< VolumeElementForcesAndSourcesCore > feRhs, int &reactionsSidesetId)

MoFEMErrorCode	atomTests (const int step_number, const VectorDouble &reactions)

Public Attributes
int	cOunter = 0

int	reactionsSidesetId

int	maxDiv

double	adaptPowCoef

double	stepRed

double	hardMod
	Isotropic Hardening modulus \( H^{i} \). More...

double	kinMod
	Kinematic Hardening modulus \( H^{k} \). More...

double	kAppa

double	yOung
	Young's modulus \( E \). More...

double	pOisson
	Poisson's ratio \( \nu \). More...

double	Mu
	shear modulus \( \mu = \frac{E}{2(1+\nu)} \) More...

double	cOefficient

double	yieldStress
	Prescribed value of the yield stress \( \sigma_{Y} \). More...

double	yieldCondition
	Yield criterion for Von Mises plasticity. More...

double	plasticMultiplier
	Plastic multiplier \( \dot{\gamma} \). More...

double	tOlerance

double	lAmbda

FTensor::Tensor4< double, 3, 3, 3, 3 >	tD

boost::shared_ptr< MatrixDouble >	gradDispPtr

MatrixDouble	consistentMatrix

boost::shared_ptr< MatrixDouble >	stressPtr

boost::shared_ptr< MatrixDouble >	plasticStrainPtr

boost::shared_ptr< VectorDouble >	hardParamPtr

boost::shared_ptr< VectorDouble >	kinParamPtr

boost::shared_ptr< VectorDouble >	kinParamPrevPtr

boost::shared_ptr< MatrixDouble >	backStressPtr
	Back stress tensor \( \boldsymbol{\beta} \). More...

PetscBool	calcReactions

PetscInt	testNb

PetscBool	isAtomTest

Tag	thPlasticStrain

Tag	thHardParam

Tag	thBackStress

Tag	thKinParam

Private Member Functions
MoFEMErrorCode	createTags (Tag &tag_ref, const string tag_string, const int &dim)

Private Attributes
MoFEM::Interface &	mField

Detailed Description

Examples: reaction_diffusion_equation.cpp, Remodeling.hpp, and SmallStrainPlasticity.hpp.

Definition at line 6 of file PlDataAtGaussPts.hpp.

Constructor & Destructor Documentation

◆ CommonData()

CommonData::CommonData ( MoFEM::Interface & m_field )

Definition at line 86 of file PlDataAtGaussPts.hpp.

                                        : mField(m_field) {
      // default base parameters
 
      kinMod            = 0.1;
      hardMod           = 0.1;
      plasticMultiplier = 0.0;
      yieldStress       = 1000;
      pOisson           = 0.1;
      yOung             = 50000;
      tOlerance         = 1e-9;
      maxDiv            = 100;
      stepRed           = 0.2;
      adaptPowCoef      = 1. / 4.;
 
      gradDispPtr      = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
      stressPtr        = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
      plasticStrainPtr = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
      hardParamPtr     = boost::shared_ptr<VectorDouble>(new VectorDouble());
      kinParamPtr      = boost::shared_ptr<VectorDouble>(new VectorDouble());
      kinParamPrevPtr  = boost::shared_ptr<VectorDouble>(new VectorDouble());
      backStressPtr    = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
 
      ierr = createTags(thPlasticStrain, "PLASTIC_STRAIN", 9);
      ierr = createTags(thHardParam, "HARDENING_PARAMETER", 1);
      ierr = createTags(thBackStress, "BACK_STRESS", 9);
      ierr = createTags(thKinParam, "KIN_HARDENING_PARAMETER", 1);
 
      CHKERRABORT(PETSC_COMM_WORLD, ierr);
   }

Member Function Documentation

◆ atomTests()

MoFEMErrorCode CommonData::atomTests	(	const int	step_number,
		const VectorDouble &	reactions
	)

Definition at line 632 of file linear_isotropic_hardening.cpp.

                                                                     {
 
   MoFEMFunctionBegin;
 
   if (testNb == 1) {
     if (step_number == 29) {
        int norm_reac;
        norm_reac =
            sqrt(reactions[0] * reactions[0] + reactions[1] * reactions[1] +
                 reactions[2] * reactions[2]);
        if (fabs(norm_reac - 3.8638e+04) > 1e-3) {
          CHKERR PetscPrintf(PETSC_COMM_WORLD, "norm error: %g \n",
                             norm_reac - 3.8638e+04);
          SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                  "atom test diverged!");
       }
     }
   }
 
   MoFEMFunctionReturn(0);
 }

◆ calcReactionForces()

MoFEMErrorCode CommonData::calcReactionForces	(	VectorDouble &	reactions,
		DM &	dm,
		boost::shared_ptr< VolumeElementForcesAndSourcesCore >	feRhs,
		int &	reactionsSidesetId
	)

Definition at line 587 of file linear_isotropic_hardening.cpp.

 {
   MoFEMFunctionBegin;
 
   Vec F_int;
   CHKERR DMCreateGlobalVector_MoFEM(dm, &F_int);
 
   feRhs->snes_ctx = FEMethod::CTX_SNESSETFUNCTION;
   feRhs->snes_f = F_int;
   CHKERR DMoFEMLoopFiniteElements(dm, "ELASTIC", feRhs);
 
   CHKERR VecAssemblyBegin(F_int);
   CHKERR VecAssemblyEnd(F_int);
   CHKERR VecGhostUpdateBegin(F_int, INSERT_VALUES, SCATTER_FORWARD);
   CHKERR VecGhostUpdateEnd(F_int, INSERT_VALUES, SCATTER_FORWARD);
 
   Reactions my_react(mField, "DM_ELASTIC", "U");
   CHKERR my_react.calculateReactions(F_int);
 
   const auto &reactions_map = my_react.getReactionsMap();
   const int id = reactionsSidesetId;
   reactions.clear();
   const bool check_key = reactions_map.find(id) == reactions_map.end();
   if (!check_key) {
     reactions = reactions_map.at(id);
   }
   Vec react_vec;
   CHKERR VecCreateMPI(mField.get_comm(), 1, PETSC_DECIDE, &react_vec);
   for (int dir : {0, 1, 2}) {
     CHKERR VecZeroEntries(react_vec);
     const int idx = mField.get_comm_rank();
     CHKERR VecSetValues(react_vec, 1, &idx, &reactions[dir], INSERT_VALUES);
     CHKERR VecAssemblyBegin(react_vec);
     CHKERR VecAssemblyEnd(react_vec);
     double sum = 0;
     CHKERR VecSum(react_vec, &sum);
     reactions[dir] = sum;
   }
   CHKERR VecDestroy(&react_vec);
   CHKERR VecDestroy(&F_int);
   MoFEMFunctionReturn(0);
 }

◆ createTags()

MoFEMErrorCode CommonData::createTags	(	Tag &	tag_ref,
		const string	tag_string,
		const int &	dim
	)

private

Definition at line 293 of file PlDataAtGaussPts.hpp.

                                                                     {
   
     MoFEMFunctionBegin;
     std::vector<double> def_val(dim, 0);
     CHKERR mField.get_moab().tag_get_handle(
         tag_string.c_str(), 1, MB_TYPE_DOUBLE, tag_ref,
         MB_TAG_CREAT | MB_TAG_VARLEN | MB_TAG_SPARSE, &def_val[0]);
     MoFEMFunctionReturn(0);
   }

◆ getParameters()

MoFEMErrorCode CommonData::getParameters ( )

Get the Parameters object.

Returns: MoFEMErrorCode

create the function to pass some command lines to the main file.

Definition at line 194 of file PlDataAtGaussPts.hpp.

                                  {
     MoFEMFunctionBegin;
 
     reactionsSidesetId = 0;
     calcReactions = PETSC_TRUE;
     testNb = 1;
     isAtomTest = PETSC_FALSE;
 
     CHKERR PetscOptionsBegin(PETSC_COMM_WORLD, "", "SimpleElasticProblem",
                              "none");
 
     CHKERR PetscOptionsScalar("-young", "Young modulus", "", yOung, &yOung,
                               PETSC_NULL);
 
     CHKERR PetscOptionsScalar("-poisson", "Poisson ratio", "", pOisson,
                               &pOisson, PETSC_NULL);
 
     CHKERR PetscOptionsScalar("-tolerance", "tolerance", "", tOlerance,
                               &tOlerance, PETSC_NULL);
 
     CHKERR PetscOptionsScalar("-yield", "yield stress", "", yieldStress,
                               &yieldStress, PETSC_NULL);
 
     CHKERR PetscOptionsScalar("-HardMod", "Hardening Modulus", "", hardMod,
                               &hardMod, PETSC_NULL);
 
     CHKERR PetscOptionsScalar("-KinMod", "Kinematic Hardening Modulus", "",
                               kinMod, &kinMod, PETSC_NULL);
 
     CHKERR PetscOptionsScalar("-step_red", "step reductcion while diverge", "",
                               stepRed, &stepRed, PETSC_NULL);
 
     CHKERR PetscOptionsScalar("-adapt_coef",
                               "code will be smarter depending on the value", "",
                               adaptPowCoef, &adaptPowCoef, PETSC_NULL);
 
     CHKERR PetscOptionsInt(
         "-calc_reactions", "calculate reactions for blocksets", "",
         reactionsSidesetId, &reactionsSidesetId, &calcReactions);
 
     CHKERR PetscOptionsInt("-is_atom_test", "Give number of the atom test", "",
                            testNb, &testNb, &isAtomTest);
 
     CHKERR PetscOptionsInt("-steps_max_div", "number of steps", "", maxDiv,
                            &maxDiv, PETSC_NULL);
 
     cOefficient = yOung / ((1 + pOisson) * (1 - 2 * pOisson));
     Mu = yOung / (2 * (1 + pOisson));
     kAppa = cOefficient * pOisson + (2. / 3.) * Mu;
      lAmbda      = cOefficient * pOisson;
      FTensor::Index<'i', 3> i;
      FTensor::Index<'j', 3> j;
      FTensor::Index<'k', 3> k;
      FTensor::Index<'l', 3> l;
      tD(i, j, k, l) = 0.;
   
     tD(0, 0, 0, 0) = lAmbda + 2 * Mu;
     tD(1, 1, 1, 1) = lAmbda + 2 * Mu;
     tD(2, 2, 2, 2) = lAmbda + 2 * Mu;
     tD(0, 0, 1, 1) = lAmbda;
     tD(0, 0, 2, 2) = lAmbda;
     tD(1, 1, 2, 2) = lAmbda;
     tD(2, 2, 1, 1) = lAmbda;
     tD(1, 1, 2, 2) = lAmbda;
     tD(2, 2, 0, 0) = lAmbda;
     tD(1, 1, 0, 0) = lAmbda;
     tD(0, 1, 0, 1) = Mu;
     tD(0, 2, 0, 2) = Mu;
     tD(1, 2, 1, 2) = Mu;
     tD(0, 1, 1, 0) = Mu;
     tD(0, 2, 2, 0) = Mu;
     tD(1, 2, 2, 1) = Mu;
     tD(2, 0, 0, 2) = Mu;
     tD(1, 0, 0, 1) = Mu;
     tD(2, 1, 1, 2) = Mu;
     tD(2, 0, 2, 0) = Mu;
     tD(1, 0, 1, 0) = Mu;
     tD(2, 1, 2, 1) = Mu;
 
 
     ierr           = PetscOptionsEnd();
     CHKERRQ(ierr);
 
     MoFEMFunctionReturn(0);
   }

◆ getTags() [1/2]

MoFEMErrorCode CommonData::getTags	(	const EntityHandle	fe_ent,
		const int	nb_gauss_pts,
		boost::shared_ptr< MatrixDouble >	tag_ptr,
		Tag	tag_ref
	)

Definition at line 116 of file PlDataAtGaussPts.hpp.

                                             {
     MoFEMFunctionBegin;
     double *tag_data;
     int tag_size;
     CHKERR mField.get_moab().tag_get_by_ptr(tag_ref, &fe_ent, 1,
                                             (const void **)&tag_data,
                                             &tag_size);
 
     if (tag_size == 1) {
       tag_ptr->resize(9, nb_gauss_pts , false);
       tag_ptr->clear();
       void const *tag_data[] = {&*tag_ptr->data().begin()};
        const int tag_size = tag_ptr->data().size();
       CHKERR mField.get_moab().tag_set_by_ptr(tag_ref, &fe_ent, 1,
                                               tag_data, &tag_size);
     } else if (tag_size != nb_gauss_pts * 9) {
       SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "Wrong size of the tag, data inconsistency");
     } else {
       MatrixAdaptor tag_vec = MatrixAdaptor(
           9, nb_gauss_pts,
           ublas::shallow_array_adaptor<double>(tag_size, tag_data));
 
       *tag_ptr = tag_vec;
     }
 
     
       MoFEMFunctionReturn(0);
   }

◆ getTags() [2/2]

MoFEMErrorCode CommonData::getTags	(	const EntityHandle	fe_ent,
		const int	nb_gauss_pts,
		boost::shared_ptr< VectorDouble >	tag_ptr,
		Tag	tag_ref
	)

Definition at line 149 of file PlDataAtGaussPts.hpp.

                                             {
     MoFEMFunctionBegin;
     double *tag_data;
     int tag_size;
     CHKERR mField.get_moab().tag_get_by_ptr(
         tag_ref, &fe_ent, 1, (const void **)&tag_data, &tag_size);
 
     if (tag_size == 1) {
       tag_ptr->resize(nb_gauss_pts);
       tag_ptr->clear();
       void const *tag_data[] = {&*tag_ptr->begin()};
       const int tag_size = tag_ptr->size();
       CHKERR mField.get_moab().tag_set_by_ptr(tag_ref, &fe_ent, 1,
                                               tag_data, &tag_size);
     } else if (tag_size != nb_gauss_pts) {
       SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "Wrong size of the tag, data inconsistency");
     } else {
       VectorAdaptor tag_vec = VectorAdaptor(
           tag_size, ublas::shallow_array_adaptor<double>(tag_size, tag_data));
       *tag_ptr = tag_vec;
     }
     MoFEMFunctionReturn(0);
   }

◆ setTags()

template<class T >

MoFEMErrorCode CommonData::setTags	(	const EntityHandle	fe_ent,
		const int	nb_gauss_pts,
		T	tag_ptr,
		Tag	tag_ref
	)

Definition at line 177 of file PlDataAtGaussPts.hpp.

                                                 {
     MoFEMFunctionBegin;
     void const *tag_data[] = {&*tag_ptr->data().begin()};
     const int tag_size = tag_ptr->data().size();
     CHKERR mField.get_moab().tag_set_by_ptr(tag_ref, &fe_ent, 1,
                                             tag_data, &tag_size);
     MoFEMFunctionReturn(0);
   }

Member Data Documentation

◆ adaptPowCoef

double CommonData::adaptPowCoef

Definition at line 12 of file PlDataAtGaussPts.hpp.

◆ backStressPtr

boost::shared_ptr<MatrixDouble> CommonData::backStressPtr

Back stress tensor \( \boldsymbol{\beta} \).

Definition at line 78 of file PlDataAtGaussPts.hpp.

◆ calcReactions

PetscBool CommonData::calcReactions

Definition at line 80 of file PlDataAtGaussPts.hpp.

◆ cOefficient

double CommonData::cOefficient

Definition at line 21 of file PlDataAtGaussPts.hpp.

◆ consistentMatrix

MatrixDouble CommonData::consistentMatrix

\( \mathbb{C} = \mathbb{D} < - \frac{\mathbb{D} : < \mathbf{n} \otimes < \mathbb{D} : < \mathbf{n}}{\mathbf{n} : < \mathbb{D} : \mathbf{n} + < \mathbf{K}^{2}} \)

Definition at line 44 of file PlDataAtGaussPts.hpp.

◆ cOunter

int CommonData::cOunter = 0

Definition at line 9 of file PlDataAtGaussPts.hpp.

◆ gradDispPtr

boost::shared_ptr<MatrixDouble> CommonData::gradDispPtr

Total strain tensor \( < \boldsymbol{\varepsilon} < \)

Definition at line 39 of file PlDataAtGaussPts.hpp.

◆ hardMod

double CommonData::hardMod

Isotropic Hardening modulus \( H^{i} \).

Definition at line 14 of file PlDataAtGaussPts.hpp.

◆ hardParamPtr

boost::shared_ptr<VectorDouble> CommonData::hardParamPtr

Vector containing the isotrtopic hardening parameter at each Gauss point

Definition at line 61 of file PlDataAtGaussPts.hpp.

◆ isAtomTest

PetscBool CommonData::isAtomTest

Definition at line 84 of file PlDataAtGaussPts.hpp.

◆ kAppa

double CommonData::kAppa

isostatic elasticity modulus \( \kappa = < \frac{E\nu}{(1-2\nu)(1 + \nu)} + \frac{2}{3}\mu \)

Definition at line 16 of file PlDataAtGaussPts.hpp.

◆ kinMod

double CommonData::kinMod

Kinematic Hardening modulus \( H^{k} \).

Definition at line 15 of file PlDataAtGaussPts.hpp.

◆ kinParamPrevPtr

boost::shared_ptr<VectorDouble> CommonData::kinParamPrevPtr

container that keeps the previous value of the kinetic hardening parameter

Definition at line 72 of file PlDataAtGaussPts.hpp.

◆ kinParamPtr

boost::shared_ptr<VectorDouble> CommonData::kinParamPtr

Vector containing the kinematic hardening parameter at each Gauss point

Definition at line 66 of file PlDataAtGaussPts.hpp.

◆ lAmbda

double CommonData::lAmbda

lame parameter \( \lambda = \frac{E \nu}{(1+\nu)(1-2\nu)} < \)

Definition at line 28 of file PlDataAtGaussPts.hpp.

◆ maxDiv

int CommonData::maxDiv

Definition at line 11 of file PlDataAtGaussPts.hpp.

◆ mField

MoFEM::Interface& CommonData::mField

private

Definition at line 291 of file PlDataAtGaussPts.hpp.

◆ Mu

double CommonData::Mu

shear modulus \( \mu = \frac{E}{2(1+\nu)} \)

Definition at line 20 of file PlDataAtGaussPts.hpp.

◆ plasticMultiplier

double CommonData::plasticMultiplier

Plastic multiplier \( \dot{\gamma} \).

Definition at line 25 of file PlDataAtGaussPts.hpp.

◆ plasticStrainPtr

boost::shared_ptr<MatrixDouble> CommonData::plasticStrainPtr

Plastic Strain \( < \boldsymbol{\varepsilon}^{p} = < \boldsymbol{\varepsilon} - < \boldsymbol{\varepsilon}^{e} \)

Definition at line 56 of file PlDataAtGaussPts.hpp.

◆ pOisson

double CommonData::pOisson

Poisson's ratio \( \nu \).

Definition at line 19 of file PlDataAtGaussPts.hpp.

◆ reactionsSidesetId

int CommonData::reactionsSidesetId

Definition at line 10 of file PlDataAtGaussPts.hpp.

◆ stepRed

double CommonData::stepRed

Definition at line 13 of file PlDataAtGaussPts.hpp.

◆ stressPtr

boost::shared_ptr<MatrixDouble> CommonData::stressPtr

Stress tensor \( < \boldsymbol{\sigma} \)

Definition at line 52 of file PlDataAtGaussPts.hpp.

◆ tD

FTensor::Tensor4<double, 3, 3, 3, 3> CommonData::tD

Elastic stiffness tensor (4th rank tensor with minor and major symmetry) \( \mathbb{D} = < \mu(\mathbb{I} + \bar{\mathbb{I}}) - < \lambda(\mathbf{I} \otimes < \mathbf{I}) \)

Definition at line 32 of file PlDataAtGaussPts.hpp.

◆ testNb

PetscInt CommonData::testNb

Definition at line 82 of file PlDataAtGaussPts.hpp.

◆ thBackStress

Tag CommonData::thBackStress

Definition at line 282 of file PlDataAtGaussPts.hpp.

◆ thHardParam

Tag CommonData::thHardParam

Definition at line 281 of file PlDataAtGaussPts.hpp.

◆ thKinParam

Tag CommonData::thKinParam

Definition at line 283 of file PlDataAtGaussPts.hpp.

◆ thPlasticStrain

Tag CommonData::thPlasticStrain

Definition at line 280 of file PlDataAtGaussPts.hpp.

◆ tOlerance

double CommonData::tOlerance

Definition at line 26 of file PlDataAtGaussPts.hpp.

◆ yieldCondition

double CommonData::yieldCondition

Yield criterion for Von Mises plasticity.

Definition at line 24 of file PlDataAtGaussPts.hpp.

◆ yieldStress

double CommonData::yieldStress

Prescribed value of the yield stress \( \sigma_{Y} \).

Definition at line 23 of file PlDataAtGaussPts.hpp.

◆ yOung

double CommonData::yOung

Young's modulus \( E \).

Definition at line 18 of file PlDataAtGaussPts.hpp.

The documentation for this struct was generated from the following files:

users_modules/linear_isotropic_hardening/src/PlDataAtGaussPts.hpp
users_modules/linear_isotropic_hardening/linear_isotropic_hardening.cpp

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes