#include <users_modules/adolc-plasticity/src/ADOLCPlasticityMaterialModels.hpp>

Inheritance diagram for ADOLCPlasticity::ParaboloidalPlasticity:

Collaboration diagram for ADOLCPlasticity::ParaboloidalPlasticity:

Public Types
enum	ModelParams { LAMBDA, MU, NUP, SIGMA_YT, SIGMA_YC, HT, HC, NT, NC, LAST_PARAM }

using	ParamsArray = std::array< double, LAST_PARAM >

Public Types inherited from ADOLCPlasticity::ClosestPointProjection
enum	TypesTags { W = 0, Y, H, LAST_TAPE }

Public Member Functions
	ParaboloidalPlasticity ()

MoFEMErrorCode	getDefaultMaterialParameters ()

MoFEMErrorCode	addMatBlockOps (MoFEM::Interface &m_field, boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string block_name, Sev sev)
	Get model parameters from blocks. More...

MoFEMErrorCode	setParams (short tag, bool &recalculate_elastic_tangent)
	Set parameters for ADOL-C of constitutive relations. More...

MoFEMErrorCode	freeHemholtzEnergy ()
	[Paraboloidal free energy] More...

MoFEMErrorCode	evalF ()
	[Paraboloidal yield function] More...

MoFEMErrorCode	yieldFunction ()
	Set yield function. More...

MoFEMErrorCode	flowPotential ()
	[Paraboloidal yield function] More...

Public Member Functions inherited from ADOLCPlasticity::ClosestPointProjection
VectorAdaptor	getPlasticStrain ()

VectorAdaptor	getTotalStrain ()

VectorAdaptor	getInternalVariables ()

VectorAdaptor	getActiveVariablesYH ()

VectorAdaptor	getStress ()

VectorAdaptor	getInternalFluxes ()

	ClosestPointProjection ()

MoFEMErrorCode	recordW ()
	Record strain energy. More...

MoFEMErrorCode	recordY ()
	Record yield function. More...

MoFEMErrorCode	recordH ()
	Record flow potential. More...

MoFEMErrorCode	playW ()

MoFEMErrorCode	playW_NoHessian ()

MoFEMErrorCode	playY ()

MoFEMErrorCode	playY_NoGradient ()

MoFEMErrorCode	playH ()

MoFEMErrorCode	playH_NoHessian ()

MoFEMErrorCode	createMatAVecR ()

MoFEMErrorCode	evaluatePotentials ()

MoFEMErrorCode	playPotentials ()

MoFEMErrorCode	playPotentials_NoHessian ()

MoFEMErrorCode	calculateR (Vec R)

MoFEMErrorCode	calculateA ()

MoFEMErrorCode	snesCreate ()
	Create nested snes. More...

MoFEMErrorCode	solveClosetProjection ()
	Solve nonlinear system of equations to find stress, internal fluxes, and Lagrange plastic multiplier. More...

MoFEMErrorCode	consistentTangent ()
	Calculate consistent tangent matrix. More...

MoFEMErrorCode	recordTapes ()
	Record tapes. More...

Public Attributes
FTensor::Index< 'i', 3 >	i

FTensor::Index< 'j', 3 >	j

FTensor::Index< 'Z', 6 >	Z

FTensor::Tensor2_symmetric< adouble, 3 >	tTotalStrain

FTensor::Tensor2_symmetric< adouble, 3 >	tPlasticStrain

FTensor::Tensor2_symmetric< adouble, 3 >	tElasticStrain

FTensor::Tensor2_symmetric< adouble, 3 >	tStress

FTensor::Tensor2_symmetric< adouble, 3 >	tDeviator

adouble	tR

double	lambda

double	mu

double	nup

double	Ht

double	Hc

double	sigmaYt

double	sigmaYc

double	nt

double	nc

boost::shared_ptr< ParamsArray >	defaultParamsArrayPtr = nullptr

adouble	I1
	[Paraboloidal free energy] More...

adouble	J2

Public Attributes inherited from ADOLCPlasticity::ClosestPointProjection
int	nbInternalVariables

boost::function< int(int, int, int)>	integrationRule

VectorDouble	internalVariables0

VectorDouble	plasticStrain0

std::array< int, LAST_TAPE >	tapesTags

VectorAdaptor	activeVariablesW

VectorAdaptor	gradientW

double	w

double	y

double	h

double	deltaGamma

MatrixDouble	Ep

MatrixDouble	Cp

MatrixDouble	Cep

ublas::symmetric_matrix< double, ublas::lower >	C

ublas::symmetric_matrix< double, ublas::lower >	D

MatrixDouble	hessianW

VectorDouble	gradientY

VectorDouble	gradientH

MatrixDouble	hessianH

MatrixDouble	partialWStrainPlasticStrain

VectorAdaptor	partialYSigma

VectorAdaptor	partialYFlux

VectorAdaptor	partialHSigma

VectorAdaptor	partialHFlux

ublas::symmetric_matrix< double, ublas::lower >	partial2HSigma

ublas::symmetric_matrix< double, ublas::lower >	partial2HFlux

MatrixDouble	partial2HSigmaFlux

SmartPetscObj< Mat >	A

SmartPetscObj< Vec >	R

SmartPetscObj< Vec >	Chi

SmartPetscObj< Vec >	dChi

ublas::matrix< double, ublas::column_major >	dataA

SmartPetscObj< SNES >	sNes

ublas::vector< adouble >	a_sTrain

ublas::vector< adouble >	a_plasticStrain

ublas::vector< adouble >	a_internalVariables

ublas::vector< adouble >	a_sTress

ublas::vector< adouble >	a_internalFluxes

adouble	a_w

adouble	a_y

adouble	a_h

Detailed Description

[J2 2D]

Paraboloidal yield criterion

See paper for reference [71] and [68]

Linear hardening

\[ \psi = \frac{1}{2} \lambda \textrm{tr}[\varepsilon]^2 + \mu \varepsilon : \varepsilon + \sigma_{yt}\alpha_0 + \frac{1}{2} H_t \alpha_0^2 + \sigma_{yc}\alpha_1 + \frac{1}{2} H_c \alpha_1^2 \]

Exponential hardening

\[ \psi = \frac{1}{2} \lambda \textrm{tr}[\varepsilon]^2 + \mu \varepsilon : \varepsilon\\ + (\sigma_{yt}+H_t)\alpha_0 + \frac{H_t}{n_t} \exp{(-n_t \alpha_0)}\\ + (\sigma_{yc}+H_c)\alpha_0 + \frac{H_c}{n_c} \exp{(-n_c \alpha_1)} \]

\[ I_1 = \textrm{tr} (\boldsymbol{\sigma}) \]

\[ \eta=\textrm{dev}[\boldsymbol{\sigma}] \]

\[ J_2 = \frac{1}{2} \eta:\eta \]

\[ y = 6J_2 + 2I_1\left(\overline{\alpha_1}-\overline{\alpha_0}\right) - 2\overline{\alpha_0} \,\overline{\alpha_1} \]

where

\[ \overline{\alpha_0}=\frac{\partial \psi}{\partial \alpha_0}=\sigma_{yt} + H_t \alpha_0 \]

\[ \overline{\alpha_1}=\frac{\partial \psi}{\partial \alpha_1}=\sigma_{yc} + H_c \alpha_1 \]

\[ \Psi = 6J_2 + 2\alpha I_1 \left(\overline{\alpha_1}-\overline{\alpha_0}\right) - 2\overline{\alpha_0} \,\overline{\alpha_1} \]

\[ \alpha= \frac{1-2\nu_p}{1+\nu_p} \]

Definition at line 531 of file ADOLCPlasticityMaterialModels.hpp.

Member Typedef Documentation

◆ ParamsArray

using ADOLCPlasticity::ParaboloidalPlasticity::ParamsArray = std::array<double, LAST_PARAM>

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 570 of file ADOLCPlasticityMaterialModels.hpp.

Member Enumeration Documentation

◆ ModelParams

enum ADOLCPlasticity::ParaboloidalPlasticity::ModelParams

Enumerator
LAMBDA
MU
NUP
SIGMA_YT
SIGMA_YC
HT
HC
NT
NC
LAST_PARAM

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 557 of file ADOLCPlasticityMaterialModels.hpp.

                    {
     LAMBDA,
     MU,
     NUP,
     SIGMA_YT,
     SIGMA_YC,
     HT,
     HC,
     NT,
     NC,
     LAST_PARAM
   };

Constructor & Destructor Documentation

◆ ParaboloidalPlasticity()

ADOLCPlasticity::ParaboloidalPlasticity::ParaboloidalPlasticity ( )

inline

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 533 of file ADOLCPlasticityMaterialModels.hpp.

                            : ClosestPointProjection() {
     nbInternalVariables = 2;
     activeVariablesW.resize(12 + nbInternalVariables, false);
   }

Member Function Documentation

◆ addMatBlockOps()

MoFEMErrorCode ADOLCPlasticity::ParaboloidalPlasticity::addMatBlockOps	(	MoFEM::Interface &	m_field,
		boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &	pip,
		std::string	block_name,
		Sev	sev
	)

inlinevirtual

Get model parameters from blocks.

Reimplemented from ADOLCPlasticity::ClosestPointProjection.

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 641 of file ADOLCPlasticityMaterialModels.hpp.

                                                 {
     return 0;
   }

◆ evalF()

MoFEMErrorCode ADOLCPlasticity::ParaboloidalPlasticity::evalF ( )

inline

[Paraboloidal yield function]

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 681 of file ADOLCPlasticityMaterialModels.hpp.

                          {
     MoFEMFunctionBegin;
     auto t_voight_stress_op = voight_to_stress_op();
     auto t_vioght_stress = getFTensor1FromPtr<6>(&a_sTress[0]);
     tStress(i, j) = t_voight_stress_op(i, j, Z) * t_vioght_stress(Z);
     tR = tStress(i, i);
     I1 = tR;
     constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<double>();
     tR /= 3.;
     tDeviator(i, j) = tStress(i, j) - tR * t_kd(i, j);
     J2 = 0.5 * (tDeviator(i, j) * tDeviator(i, j));
     MoFEMFunctionReturn(0);
   }

◆ flowPotential()

MoFEMErrorCode ADOLCPlasticity::ParaboloidalPlasticity::flowPotential ( )

inlinevirtual

[Paraboloidal yield function]

[Paraboloidal flow potential]

Implements ADOLCPlasticity::ClosestPointProjection.

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 705 of file ADOLCPlasticityMaterialModels.hpp.

                                  {
     MoFEMFunctionBegin;
     CHKERR evalF();
     double alpha =
         (1 - 2 * nup) /
         (1 + nup); // relation between alpha and plastic poisson's ratio
     a_h = 6 * J2 +
           2 * alpha * I1 * (a_internalFluxes[1] - a_internalFluxes[0]) -
           2 * a_internalFluxes[0] * a_internalFluxes[1];
     MoFEMFunctionReturn(0);
   }

◆ freeHemholtzEnergy()

MoFEMErrorCode ADOLCPlasticity::ParaboloidalPlasticity::freeHemholtzEnergy ( )

inlinevirtual

[Paraboloidal free energy]

Implements ADOLCPlasticity::ClosestPointProjection.

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 653 of file ADOLCPlasticityMaterialModels.hpp.

                                       {
     MoFEMFunctionBegin;
  
     auto t_voight_strain_op = voight_to_strain_op();
     auto t_voight_total_strain = getFTensor1FromPtr<6>(&a_sTrain[0]);
     auto t_voight_plastic_strain = getFTensor1FromPtr<6>(&a_plasticStrain[0]);
  
     tTotalStrain(i, j) = t_voight_strain_op(i, j, Z) * t_voight_total_strain(Z);
     tPlasticStrain(i, j) =
         t_voight_strain_op(i, j, Z) * t_voight_plastic_strain(Z);
     tElasticStrain(i, j) = tTotalStrain(i, j) - tPlasticStrain(i, j);
     tR = tElasticStrain(i, i);
  
     a_w = 0.5 * lambda * tR * tR +
           mu * (tElasticStrain(i, j) * tElasticStrain(i, j));
  
     a_w += sigmaYt * a_internalVariables[0] +
            0.5 * Ht * a_internalVariables[0] * a_internalVariables[0];
     a_w += sigmaYc * a_internalVariables[1] +
            0.5 * Hc * a_internalVariables[1] * a_internalVariables[1];
  
     MoFEMFunctionReturn(0);
   }

◆ getDefaultMaterialParameters()

MoFEMErrorCode ADOLCPlasticity::ParaboloidalPlasticity::getDefaultMaterialParameters ( )

inline

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 573 of file ADOLCPlasticityMaterialModels.hpp.

                                                 {
     MoFEMFunctionBegin;
  
     double young_modulus = 1;
     double poisson_ratio = 0.25;
  
     sigmaYt = 1;
     sigmaYc = 1;
     Ht = 0.1;
     Hc = 0.1;
     nc = 1.;
     nt = 1.;
     nup = 0.;
  
     CHKERR PetscOptionsGetScalar(PETSC_NULL, "-young_modulus", &young_modulus,
                                  PETSC_NULL);
     CHKERR PetscOptionsGetScalar(PETSC_NULL, "-poisson_ratio", &poisson_ratio,
                                  0);
     CHKERR PetscOptionsGetScalar(PETSC_NULL, "-sigma_yt", &sigmaYt, PETSC_NULL);
     CHKERR PetscOptionsGetScalar(PETSC_NULL, "-sigma_yc", &sigmaYc, PETSC_NULL);
     CHKERR PetscOptionsGetScalar(PETSC_NULL, "-Ht", &Ht, PETSC_NULL);
     CHKERR PetscOptionsGetScalar(PETSC_NULL, "-Hc", &Hc, PETSC_NULL);
     CHKERR PetscOptionsGetScalar(PETSC_NULL, "-nt", &nt, PETSC_NULL);
     CHKERR PetscOptionsGetScalar(PETSC_NULL, "-nc", &nc, PETSC_NULL);
     CHKERR PetscOptionsGetScalar(PETSC_NULL, "-nup", &nup, PETSC_NULL);
  
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "Young modulus: " << young_modulus;
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "Poisson ratio: " << poisson_ratio;
  
     mu = MU(young_modulus, poisson_ratio);
     lambda = LAMBDA(young_modulus, poisson_ratio);
  
     defaultParamsArrayPtr = boost::make_shared<ParamsArray>();
     (*defaultParamsArrayPtr)[LAMBDA] = lambda;
     (*defaultParamsArrayPtr)[MU] = mu;
     (*defaultParamsArrayPtr)[NUP] = nup;
     (*defaultParamsArrayPtr)[SIGMA_YT] = sigmaYt;
     (*defaultParamsArrayPtr)[SIGMA_YC] = sigmaYc;
     (*defaultParamsArrayPtr)[HT] = Ht;
     (*defaultParamsArrayPtr)[HC] = Hc;
     (*defaultParamsArrayPtr)[NT] = nt;
     (*defaultParamsArrayPtr)[NC] = nc;
  
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "LAMBDA: " << (*defaultParamsArrayPtr)[LAMBDA];
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "MU: " << (*defaultParamsArrayPtr)[MU];
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "NUP: " << (*defaultParamsArrayPtr)[NUP];
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "SIGMA_YT: " << (*defaultParamsArrayPtr)[SIGMA_YT];
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "SIGMA_YC: " << (*defaultParamsArrayPtr)[SIGMA_YC];
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "HT: " << (*defaultParamsArrayPtr)[HT];
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "HC: " << (*defaultParamsArrayPtr)[HC];
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "NT: " << (*defaultParamsArrayPtr)[NT];
     MOFEM_LOG("ADOLCPlasticityWold", Sev::inform)
         << "NC: " << (*defaultParamsArrayPtr)[NC];
  
     MoFEMFunctionReturn(0);
   };

◆ setParams()

MoFEMErrorCode ADOLCPlasticity::ParaboloidalPlasticity::setParams	(	short	tag,
		bool &	recalculate_elastic_tangent
	)

inlinevirtual

Set parameters for ADOL-C of constitutive relations.

Parameters

tag	[in] - tag of the tape
recalculate_elastic_tangent	[out] - if setParam set it to true, tangent matrix for elastic domain should be recalculated

Reimplemented from ADOLCPlasticity::ClosestPointProjection.

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 647 of file ADOLCPlasticityMaterialModels.hpp.

                                                                          {
     recalculate_elastic_tangent = true;
     return 0;
   }

◆ yieldFunction()

MoFEMErrorCode ADOLCPlasticity::ParaboloidalPlasticity::yieldFunction ( )

inlinevirtual

Set yield function.

Implements ADOLCPlasticity::ClosestPointProjection.

Examples: ADOLCPlasticityMaterialModels.hpp.

Definition at line 695 of file ADOLCPlasticityMaterialModels.hpp.

                                  {
     MoFEMFunctionBegin;
     CHKERR evalF();
     a_y = 6 * J2 + 2 * I1 * (a_internalFluxes[1] - a_internalFluxes[0]) -
           2 * a_internalFluxes[0] * a_internalFluxes[1];
     MoFEMFunctionReturn(0);
   }