ADOLCPlasticity::J2Plasticity< 2 > Struct Reference

[J2 2D] More...

#include "users_modules/adolc-plasticity/src/ADOLCPlasticityMaterialModels.hpp"

Inheritance diagram for ADOLCPlasticity::J2Plasticity< 2 >:

Collaboration diagram for ADOLCPlasticity::J2Plasticity< 2 >:

Public Member Functions
MoFEMErrorCode	freeHemholtzEnergy ()
	Set Hemholtz energy.
MoFEMErrorCode	codedHessianW (vector< double * > hessian)
Public Member Functions inherited from ADOLCPlasticity::J2Plasticity< 3 >
	J2Plasticity ()
MoFEMErrorCode	getDefaultMaterialParameters ()
MoFEMErrorCode	addMatBlockOps (MoFEM::Interface &m_field, boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pip, std::string block_name, Sev sev)
	Get material parameters from block.
MoFEMErrorCode	setParams (short tag, bool &recalculate_elastic_tangent)
	Set material parameters at integration point.
MoFEMErrorCode	freeHemholtzEnergy ()
	[Free energy J2]
MoFEMErrorCode	evalF ()
	[Yield function J2]
MoFEMErrorCode	yieldFunction ()
	Set yield function.
MoFEMErrorCode	flowPotential ()
	[Yield function J2]
MoFEMErrorCode	codedHessianW (vector< double * > hessian)
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.
MoFEMErrorCode	recordY ()
	Record yield function.
MoFEMErrorCode	recordH ()
	Record flow potential.
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.
MoFEMErrorCode	solveClosestProjection ()
	Solve nonlinear system of equations to find stress, internal fluxes, and Lagrange plastic multiplier.
MoFEMErrorCode	consistentTangent ()
	Calculate consistent tangent matrix.
MoFEMErrorCode	recordTapes ()
	Record tapes.

Additional Inherited Members
Public Types inherited from ADOLCPlasticity::J2Plasticity< 3 >
enum	ModelParams {   LAMBDA , MU , ISOH , KINK ,   PHI , SIGMA_Y , LAST_PARAM }
using	ParamsArray = std::array< double, LAST_PARAM >
Public Types inherited from ADOLCPlasticity::ClosestPointProjection
enum	TypesTags { W = 0 , Y , H , LAST_TAPE }
Public Attributes inherited from ADOLCPlasticity::J2Plasticity< 3 >
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 >	tInternalTensor
FTensor::Tensor2_symmetric< adouble, 3 >	tStress
FTensor::Tensor2_symmetric< adouble, 3 >	tBackStress
FTensor::Tensor2_symmetric< adouble, 3 >	tDeviator
adouble	tR
double	mu
double	lambda
double	H
double	K
double	phi
double	sigmaY
boost::shared_ptr< ParamsArray >	defaultParamsArrayPtr = nullptr
boost::shared_ptr< ParamsArray >	paramsArrayPtr = nullptr
boost::shared_ptr< ParamsArray >	oldParamsArrayPtr = nullptr
adouble	f
	[Free energy 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
PetscBool	implHessianW
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]

Definition at line 467 of file ADOLCPlasticityMaterialModels.hpp.

Member Function Documentation

codedHessianW()

MoFEMErrorCode ADOLCPlasticity::J2Plasticity< 2 >::codedHessianW ( vector< double * >  hessian )

inlinevirtual

Reimplemented from ADOLCPlasticity::ClosestPointProjection.

Definition at line 518 of file ADOLCPlasticityMaterialModels.hpp.

                                                         {
    MoFEMFunctionBegin;
    // Not implemented yet
    SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
            "2D Energy Hessian not yet implemented");
    MoFEMFunctionReturn(0);
  }

freeHemholtzEnergy()

MoFEMErrorCode ADOLCPlasticity::J2Plasticity< 2 >::freeHemholtzEnergy ( )

inlinevirtual

Set Hemholtz energy.

[Plane stress]

[Plane stress]

Implements ADOLCPlasticity::ClosestPointProjection.

Definition at line 470 of file ADOLCPlasticityMaterialModels.hpp.

                                      {
    MoFEMFunctionBegin;
 
    auto p_lambda = mkparam(lambda);  // 0
    auto p_mu = mkparam(mu);          // 1
    auto p_H = mkparam(H);            // 2
    auto p_K = mkparam(K);            // 3
    auto p_phi = mkparam(phi);        // 4
    auto p_sigma_y = mkparam(sigmaY); // 5
 
    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]);
    auto t_voight_internal_tensor =
        getFTensor1FromPtr<6>(&a_internalVariables[1]);
 
    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);
    tInternalTensor(i, j) =
        t_voight_strain_op(i, j, Z) * t_voight_internal_tensor(Z);
 
    tElasticStrain(i, j) = tTotalStrain(i, j) - tPlasticStrain(i, j);
 
    //! [Plane stress]
    // Calculate st  rain at ezz enforcing plane strain case
    tElasticStrain(2, 2) =
        (-tElasticStrain(0, 0) * p_lambda - tElasticStrain(1, 1) * p_lambda) /
        (p_lambda + 2 * p_mu);
    //! [Plane stress]
 
    tR = tElasticStrain(i, i);
 
    a_w = 0.5 * p_lambda * tR * tR +
          p_mu * (tElasticStrain(i, j) * tElasticStrain(i, j));
 
    // plastic part
    // isotropic
    a_w +=
        a_internalVariables[0] * p_sigma_y +
        0.5 * p_phi * p_H * (a_internalVariables[0] * a_internalVariables[0]);
    // kinematic
    a_w += (0.5 * (1 - p_phi) * p_K) *
           (tInternalTensor(i, j) * tInternalTensor(i, j));
 
    MoFEMFunctionReturn(0);
  }

---

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

• users_modules/adolc-plasticity/src/ADOLCPlasticityMaterialModels.hpp