J2 plasticity (Kinematic Isotropic (Linear) Hardening) More...

#include <users_modules/small_strain_plasticity/src/SmallStrainPlasticityMaterialModels.hpp>

Inheritance diagram for SmallStrainJ2Plasticity:

Collaboration diagram for SmallStrainJ2Plasticity:

Public Member Functions
	SmallStrainJ2Plasticity ()

virtual PetscErrorCode	freeHemholtzEnergy ()
	Available energy. More...

PetscErrorCode	evalF ()

virtual PetscErrorCode	yieldFunction ()
	Evaluate yield function. More...

virtual PetscErrorCode	flowPotential ()
	Associated flow potential. More...

Public Member Functions inherited from SmallStrainPlasticity::ClosestPointProjection
	ClosestPointProjection ()

virtual PetscErrorCode	setActiveVariablesW ()

virtual PetscErrorCode	setActiveVariablesYH ()

virtual PetscErrorCode	rEcordW ()

virtual PetscErrorCode	rEcordY ()

virtual PetscErrorCode	rEcordH ()

virtual PetscErrorCode	pLayW ()

virtual PetscErrorCode	pLayW_NoHessian ()

virtual PetscErrorCode	pLayY ()

virtual PetscErrorCode	pLayY_NoGradient ()

virtual PetscErrorCode	pLayH ()

virtual PetscErrorCode	createMatAVecR ()

virtual PetscErrorCode	destroyMatAVecR ()

virtual PetscErrorCode	evaluatePotentials ()

virtual PetscErrorCode	cAlculateR (Vec R)

virtual PetscErrorCode	cAlculateA ()

PetscErrorCode	snesCreate ()

PetscErrorCode	snesDestroy ()

PetscErrorCode	solveColasetProjection ()

PetscErrorCode	consistentTangent ()

Public Attributes
double	mu

double	lambda

double	H

double	K

double	phi

double	sIgma_y

adouble	t
	Auxiliary function. More...

adouble	f

Public Attributes inherited from SmallStrainPlasticity::ClosestPointProjection
VectorDouble	sTrain

VectorDouble	internalVariables

VectorDouble	internalVariables0

VectorDouble	plasticStrain0

VectorDouble	plasticStrain

double	deltaGamma

double	tOl

int	gG

VectorAdaptor	sTress

VectorAdaptor	internalFluxes

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

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

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

vector< int >	tAgs

VectorDouble	activeVariablesW

VectorDouble	activeVariablesYH

double	w

double	y

double	h

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

VectorDouble	gradientW

MatrixDouble	hessianW

VectorDouble	gradientY

VectorDouble	gradientH

MatrixDouble	hessianH

Mat	A

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

Vec	R

Vec	Chi

Vec	dChi

SNES	sNes

KSP	kSp

PC	pC

SNESLineSearch	lineSearch

MatrixDouble	Ep

MatrixDouble	Lp

MatrixDouble	Cp

MatrixDouble	Cep

VectorDouble	Yp

Detailed Description

J2 plasticity (Kinematic Isotropic (Linear) Hardening)

Examples: SmallStrainPlasticityMaterialModels.hpp.

Definition at line 32 of file SmallStrainPlasticityMaterialModels.hpp.

Constructor & Destructor Documentation

◆ SmallStrainJ2Plasticity()

SmallStrainJ2Plasticity::SmallStrainJ2Plasticity ( )

inline

Examples: SmallStrainPlasticityMaterialModels.hpp.

Definition at line 34 of file SmallStrainPlasticityMaterialModels.hpp.

                            :
   SmallStrainPlasticity::ClosestPointProjection() {
     internalVariables.resize(7,false);
     internalVariables.clear();
   }

Member Function Documentation

◆ evalF()

PetscErrorCode SmallStrainJ2Plasticity::evalF ( )

inline

Examples: SmallStrainPlasticityMaterialModels.hpp.

Definition at line 107 of file SmallStrainPlasticityMaterialModels.hpp.

                          {
     PetscFunctionBegin;
     adouble tR = 0;
     for(int dd = 0;dd<3;dd++) {
       tR += a_sTress[dd];
     }
     // Fix the constant to mach uniaxial test
     const double c = 3./2.;
     tR /= 3.;
     f = 0;
     for(int dd = 0;dd<3;dd++) {
       t = (a_sTress[dd]-tR)-a_internalFluxes[dd+1];
       f += c*t*t;
     }
     for(int dd = 3;dd<6;dd++) {
       t = a_sTress[dd]-a_internalFluxes[dd+1];
       f += c*2.*t*t; // s:s off diagonal terms need to go twice
     }
     PetscFunctionReturn(0);
   }

◆ flowPotential()

virtual PetscErrorCode SmallStrainJ2Plasticity::flowPotential ( )

inlinevirtual

Associated flow potential.

See yieldFunction.

Reimplemented from SmallStrainPlasticity::ClosestPointProjection.

Examples: SmallStrainPlasticityMaterialModels.hpp.

Definition at line 148 of file SmallStrainPlasticityMaterialModels.hpp.

                                          {
     PetscFunctionBegin;
     ierr = evalF(); CHKERRQ(ierr);
     a_h = sqrt(f) - a_internalFluxes[0];
     PetscFunctionReturn(0);
   }

◆ freeHemholtzEnergy()

virtual PetscErrorCode SmallStrainJ2Plasticity::freeHemholtzEnergy ( )

inlinevirtual

Available energy.

\[ \psi = \frac{1}{2} \lambda \textrm{tr}[\varepsilon]^2 + \mu \varepsilon : \varepsilon + \sigma_y\alpha + \frac{1}{2} \phi H \alpha^2 + \frac{1}{2} (1-\phi)K \beta^2 \]

Isotropic hardening variable \(\alpha\) is in first index of internal variables vector. Kinematic hardening variable is in the remaining indices of internal variables vector.

Reimplemented from SmallStrainPlasticity::ClosestPointProjection.

Examples: SmallStrainPlasticityMaterialModels.hpp.

Definition at line 65 of file SmallStrainPlasticityMaterialModels.hpp.

                                               {
     PetscFunctionBegin;
     //elastic part
     adouble tR = 0;
     for(int dd = 0;dd<3;dd++) {
       tR += a_sTrain[dd]-a_plasticStrain[dd];
     }
     a_w = 0.5*lambda*tR*tR;
     adouble t;
     int dd = 0;
     for(;dd<3;dd++) {
       t = a_sTrain[dd]-a_plasticStrain[dd];
       a_w += mu*t*t;
     }
     for(;dd<6;dd++) {
       t = a_sTrain[dd]-a_plasticStrain[dd];
       a_w += 0.5*mu*t*t;
     }
     //plastic part
     //isotropic
     a_w += a_internalVariables[0]*sIgma_y+0.5*phi*H*a_internalVariables[0]*a_internalVariables[0];
     //kinematic
     for(unsigned int dd = 1;dd<a_internalVariables.size();dd++) {
      a_w += 0.5*(1-phi)*K*a_internalVariables[dd]*a_internalVariables[dd];
     }
     PetscFunctionReturn(0);
   }

◆ yieldFunction()

virtual PetscErrorCode SmallStrainJ2Plasticity::yieldFunction ( )

inlinevirtual

Evaluate yield function.

\[ y = \sqrt{f} - \overline{\alpha} \]

where \(f\) is defined in evalF.

Reimplemented from SmallStrainPlasticity::ClosestPointProjection.

Examples: SmallStrainPlasticityMaterialModels.hpp.

Definition at line 136 of file SmallStrainPlasticityMaterialModels.hpp.

                                          {
     PetscFunctionBegin;
     ierr = evalF(); CHKERRQ(ierr);
     a_y = sqrt(f) - a_internalFluxes[0];
     PetscFunctionReturn(0);
   }