Closest point projection algorithm. More...

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

Inheritance diagram for ADOLCPlasticity::ClosestPointProjection:

Collaboration diagram for ADOLCPlasticity::ClosestPointProjection:

Public Types
enum	TypesTags { W = 0, Y, H, LAST_TAPE }

Public Member Functions
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...

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

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

virtual MoFEMErrorCode	freeHemholtzEnergy ()=0
	Set Hemholtz energy. More...

virtual MoFEMErrorCode	yieldFunction ()=0
	Set yield function. More...

virtual MoFEMErrorCode	flowPotential ()=0
	Set flow potential. More...

Public Attributes
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

Friends
MoFEMErrorCode	ADOLCPlasticityRes (SNES, Vec, Vec, void *ctx)
	Function executed by nested SENES at evaluationg residual. More...

MoFEMErrorCode	ADOLCPlasticityJac (SNES, Vec, Mat, Mat, void *ctx)
	Function executed by nested SENES at evaluationg tangent matrix. More...

Detailed Description

Closest point projection algorithm.

Definition at line 139 of file ADOLCPlasticity.hpp.

Member Enumeration Documentation

◆ TypesTags

enum ADOLCPlasticity::ClosestPointProjection::TypesTags

Enumerator
W
Y
H
LAST_TAPE

Examples: ADOLCPlasticity.hpp.

Definition at line 168 of file ADOLCPlasticity.hpp.

168 { W = 0, Y, H, LAST_TAPE }; //< W - energy, Y - yield, H - flow

Constructor & Destructor Documentation

◆ ClosestPointProjection()

ADOLCPlasticity::ClosestPointProjection::ClosestPointProjection ( )

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 15 of file ClosetPointProjection.cpp.

                                                {
  
   if (!LogManager::checkIfChannelExist("ADOLCPlasticityWold")) {
     auto core_log = logging::core::get();
  
     core_log->add_sink(LogManager::createSink(LogManager::getStrmWorld(),
                                               "ADOLCPlasticityWold"));
     core_log->add_sink(LogManager::createSink(LogManager::getStrmSync(),
                                               "ADOLCPlasticitySync"));
     core_log->add_sink(LogManager::createSink(LogManager::getStrmSelf(),
                                               "ADOLCPlasticitySelf"));
  
     LogManager::setLog("ADOLCPlasticityWold");
     LogManager::setLog("ADOLCPlasticitySync");
     LogManager::setLog("ADOLCPlasticitySelf");
  
     MOFEM_LOG_TAG("ADOLCPlasticityWold", "ADOL-C Plasticity");
     MOFEM_LOG_TAG("ADOLCPlasticitySync", "ADOL-C Plasticity");
     MOFEM_LOG_TAG("ADOLCPlasticitySelf", "ADOL-C Plasticity");
   }
  
   tapesTags = {0, 1, 2}; //< set tapes default tags
  
 }

Member Function Documentation

◆ addMatBlockOps()

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

inlinevirtual

Get model parameters from blocks.

Reimplemented in ADOLCPlasticity::ParaboloidalPlasticity, and ADOLCPlasticity::J2Plasticity< 3 >.

Examples: ADOLCPlasticity.hpp.

Definition at line 261 of file ADOLCPlasticity.hpp.

                                                 {
     return 0;
   }

◆ calculateA()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::calculateA ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 382 of file ClosetPointProjection.cpp.

                                                   {
   MoFEMFunctionBegin;
   // matrix A
   CHKERR MatZeroEntries(A);
   // row 0 (Strain)
   MatrixDouble a00 = prod(partial2HSigma, partialWStrainPlasticStrain);
   MatrixDouble a01 = prod(partial2HSigmaFlux, D);
   for (int ii = 0; ii < 6; ii++) {
     for (int jj = 0; jj < 6; jj++) {
       dataA(ii, jj) = deltaGamma * a00(ii, jj);
     }
     for (unsigned int jj = 0; jj < nbInternalVariables; jj++) {
       dataA(ii, 6 + jj) = deltaGamma * a01(ii, jj);
     }
     dataA(ii, ii) -= 1;
   }
   for (int ii = 0; ii < 6; ii++) {
     dataA(ii, 6 + nbInternalVariables) = partialHSigma[ii];
   }
   // row 1 (Fluxes)
   MatrixDouble a11 = prod(partial2HFlux, D);
   MatrixDouble a10 =
       prod(trans(partial2HSigmaFlux), partialWStrainPlasticStrain);
   for (unsigned int ii = 0; ii < nbInternalVariables; ii++) {
     for (int jj = 0; jj < 6; jj++) {
       dataA(6 + ii, jj) += deltaGamma * a10(ii, jj);
     }
     for (unsigned int jj = 0; jj < nbInternalVariables; jj++) {
       dataA(6 + ii, 6 + jj) += deltaGamma * a11(ii, jj);
     }
     dataA(6 + ii, 6 + ii) -= 1;
   }
   for (unsigned int ii = 0; ii < nbInternalVariables; ii++) {
     dataA(6 + ii, 6 + nbInternalVariables) = partialHFlux[ii];
   }
   // row 3 (Plastic multiplier)
   VectorDouble partial_y_sigma_c =
       prod(trans(partialWStrainPlasticStrain), partialYSigma);
   VectorDouble partial_y_flux_d = prod(trans(D), partialYFlux);
   for (unsigned int dd = 0; dd < partial_y_sigma_c.size(); dd++) {
     dataA(6 + nbInternalVariables, dd) = partial_y_sigma_c[dd];
   }
   for (unsigned int dd = 0; dd < partial_y_flux_d.size(); dd++) {
     dataA(6 + nbInternalVariables, 6 + dd) = partial_y_flux_d[dd];
   }
   dataA(6 + nbInternalVariables, 6 + nbInternalVariables) = 0;
   CHKERR MatAssemblyBegin(A, MAT_FINAL_ASSEMBLY);
   CHKERR MatAssemblyEnd(A, MAT_FINAL_ASSEMBLY);
   MoFEMFunctionReturn(0);
 }

◆ calculateR()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::calculateR ( Vec R )

Examples: ADOLCPlasticity.hpp.

Definition at line 360 of file ClosetPointProjection.cpp.

                                         {
   MoFEMFunctionBegin;
  
   auto plastic_strain = getPlasticStrain();
   auto internal_variables = getInternalVariables();
  
   // Residual
   double *array;
   CHKERR VecGetArray(R, &array);
   for (int ii = 0; ii < 6; ii++) {
     array[ii] = -plastic_strain[ii] + plasticStrain0[ii] +
                 deltaGamma * partialHSigma[ii];
   }
   for (unsigned int ii = 0; ii < nbInternalVariables; ii++) {
     array[6 + ii] = -internal_variables[ii] + internalVariables0[ii] +
                     deltaGamma * partialHFlux[ii];
   }
   array[6 + nbInternalVariables] = y;
   CHKERR VecRestoreArray(R, &array);
   MoFEMFunctionReturn(0);
 }

◆ consistentTangent()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::consistentTangent ( )

Calculate consistent tangent matrix.

Examples: ADOLCPlasticity.hpp.

Definition at line 509 of file ClosetPointProjection.cpp.

                                           {
   MoFEMFunctionBegin;
   CHKERR playPotentials();
   CHKERR calculateA();
   Ep.resize(6, 6, false);
   MatrixDouble ep_row; // inrement internal varaibles as result of incremen of
                        // total strain
   ep_row = deltaGamma * prod(partial2HSigma, C);
   MatrixDouble alpha_row;
   alpha_row = deltaGamma * prod(trans(partial2HSigmaFlux), C);
   VectorDouble y_row;
   y_row = prod(C, partialYSigma);
   const int n = 6 + nbInternalVariables;
   // Iterate strains
   for (int dd = 0; dd < 6; dd++) {
     CHKERR VecZeroEntries(R);
     double *array;
     CHKERR VecGetArray(R, &array);
     {
       for (auto ii = 0; ii < 6; ii++) {
         array[ii] = ep_row(ii, dd);
       }
       for (auto ii = 0; ii < nbInternalVariables; ii++) {
         array[6 + ii] = alpha_row(ii, dd);
       }
       array[n] = y_row[dd];
     }
     CHKERR VecRestoreArray(R, &array);
     CHKERR VecAssemblyBegin(R);
     CHKERR VecAssemblyEnd(R);
  
     KSP ksp;
     CHKERR SNESGetKSP(sNes, &ksp);
     CHKERR KSPSolve(ksp, R, dChi);
     CHKERR VecGetArray(dChi, &array);
     {
       for (auto ii = 0; ii < 6; ii++) {
         Ep(ii, dd) = array[ii];
       }
     }
     CHKERR VecRestoreArray(dChi, &array);
   }
   Cp.resize(6, 6, false);
   noalias(Cp) = prod(C, Ep);
   Cep.resize(6, 6, false);
   noalias(Cep) = C + Cp;
   MoFEMFunctionReturn(0);
 }

◆ createMatAVecR()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::createMatAVecR ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 314 of file ClosetPointProjection.cpp.

                                                       {
   MoFEMFunctionBegin;
   PetscInt n;
   n = 1 + 6 + nbInternalVariables;
   dataA.resize(n, n, false);
   Mat A_tmp;
   Vec R_tmp, Chi_tmp, dChi_tmp;
   CHKERR MatCreateSeqDense(PETSC_COMM_SELF, n, n, &dataA(0, 0), &A_tmp);
   CHKERR VecCreateSeq(PETSC_COMM_SELF, n, &R_tmp);
   CHKERR VecCreateSeq(PETSC_COMM_SELF, n, &Chi_tmp);
   CHKERR VecCreateSeq(PETSC_COMM_SELF, n, &dChi_tmp);
  
   A = SmartPetscObj<Mat>(A_tmp);
   R = SmartPetscObj<Vec>(R_tmp);
   Chi = SmartPetscObj<Vec>(Chi_tmp);
   dChi = SmartPetscObj<Vec>(dChi_tmp);
  
   MoFEMFunctionReturn(0);
 }

◆ evaluatePotentials()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::evaluatePotentials ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 335 of file ClosetPointProjection.cpp.

                                            {
   MoFEMFunctionBegin;
   CHKERR recordW();
   CHKERR recordY();
   CHKERR recordH();
   MoFEMFunctionReturn(0);
 }

◆ flowPotential()

virtual MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::flowPotential ( )

pure virtual

Set flow potential.

Implemented in ADOLCPlasticity::ParaboloidalPlasticity, and ADOLCPlasticity::J2Plasticity< 3 >.

Examples: ADOLCPlasticity.hpp.

◆ freeHemholtzEnergy()

virtual MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::freeHemholtzEnergy ( )

pure virtual

Set Hemholtz energy.

Implemented in ADOLCPlasticity::ParaboloidalPlasticity, ADOLCPlasticity::J2Plasticity< 2 >, and ADOLCPlasticity::J2Plasticity< 3 >.

Examples: ADOLCPlasticity.hpp.

◆ getActiveVariablesYH()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getActiveVariablesYH ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 158 of file ADOLCPlasticity.hpp.

                                               {
     return getVectorAdaptor(&(gradientW[6]), 6 + nbInternalVariables);
   }

◆ getInternalFluxes()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getInternalFluxes ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 164 of file ADOLCPlasticity.hpp.

                                            {
     return getVectorAdaptor(&(gradientW[12]), nbInternalVariables);
   }

◆ getInternalVariables()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getInternalVariables ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 155 of file ADOLCPlasticity.hpp.

                                               {
     return getVectorAdaptor(&(activeVariablesW[12]), nbInternalVariables);
   }

◆ getPlasticStrain()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getPlasticStrain ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 149 of file ADOLCPlasticity.hpp.

                                           {
     return getVectorAdaptor(&(activeVariablesW[0]), 6);
   }

◆ getStress()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getStress ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 161 of file ADOLCPlasticity.hpp.

                                    {
     return getVectorAdaptor(&(gradientW[6]), 6);
   }

◆ getTotalStrain()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getTotalStrain ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 152 of file ADOLCPlasticity.hpp.

                                         {
     return getVectorAdaptor(&(activeVariablesW[6]), 6);
   }

◆ playH()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playH ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 230 of file ClosetPointProjection.cpp.

                                              {
   MoFEMFunctionBegin;
   auto active_variables_yh = getActiveVariablesYH();
   int r;
   int adloc_return_value = 0;
   r = ::function(tapesTags[ClosestPointProjection::H], 1,
                  active_variables_yh.size(), &active_variables_yh[0], &h);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   gradientH.resize(active_variables_yh.size());
   r = gradient(tapesTags[ClosestPointProjection::H], active_variables_yh.size(),
                &active_variables_yh[0], &gradientH[0]);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   hessianH.resize(6 + nbInternalVariables, 6 + nbInternalVariables, false);
   hessianH.clear();
   vector<double *> hessian_h(active_variables_yh.size());
   for (int dd = 0; dd < active_variables_yh.size(); dd++) {
     hessian_h[dd] = &hessianH(dd, 0);
   }
   r = hessian(tapesTags[ClosestPointProjection::H], active_variables_yh.size(),
               &active_variables_yh[0], &hessian_h[0]);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   partialHSigma =
       VectorAdaptor(6, ublas::shallow_array_adaptor<double>(6, &gradientH[0]));
   partialHFlux = VectorAdaptor(
       nbInternalVariables,
       ublas::shallow_array_adaptor<double>(nbInternalVariables, &gradientH[6]));
   partialHFlux *= -1;
   partial2HSigma.resize(6, 6, false);
   for (int ii = 0; ii < 6; ii++) {
     for (int jj = 0; jj <= ii; jj++) {
       partial2HSigma(ii, jj) = hessianH(ii, jj);
     }
   }
   partial2HFlux.resize(nbInternalVariables, nbInternalVariables, false);
   for (unsigned int ii = 0; ii < nbInternalVariables; ii++) {
     for (unsigned int jj = 0; jj <= ii; jj++) {
       partial2HFlux(ii, jj) = -hessianH(6 + ii, 6 + jj);
     }
   }
   partial2HSigmaFlux.resize(6, nbInternalVariables, false);
   partial2HSigmaFlux.clear();
   for (int ii = 0; ii < 6; ii++) {
     for (unsigned int jj = 0; jj < nbInternalVariables; jj++) {
       partial2HSigmaFlux(ii, jj) = -hessianH(6 + jj, ii);
     }
   }
   MoFEMFunctionReturn(0);
 }

◆ playH_NoHessian()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playH_NoHessian ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 287 of file ClosetPointProjection.cpp.

                                                        {
   MoFEMFunctionBegin;
   auto active_variables_yh = getActiveVariablesYH();
   int r;
   int adloc_return_value = 0;
   r = ::function(tapesTags[ClosestPointProjection::H], 1,
                  active_variables_yh.size(), &active_variables_yh[0], &h);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   gradientH.resize(active_variables_yh.size());
   r = gradient(tapesTags[ClosestPointProjection::H], active_variables_yh.size(),
                &active_variables_yh[0], &gradientH[0]);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   partialHSigma =
       VectorAdaptor(6, ublas::shallow_array_adaptor<double>(6, &gradientH[0]));
   partialHFlux = VectorAdaptor(
       nbInternalVariables,
       ublas::shallow_array_adaptor<double>(nbInternalVariables, &gradientH[6]));
   partialHFlux *= -1;
   MoFEMFunctionReturn(0);
 }

◆ playPotentials()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playPotentials ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 343 of file ClosetPointProjection.cpp.

                                                       {
   MoFEMFunctionBegin;
   CHKERR playW();
   CHKERR playY();
   CHKERR playH();
   MoFEMFunctionReturn(0);
 }

◆ playPotentials_NoHessian()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playPotentials_NoHessian ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 351 of file ClosetPointProjection.cpp.

                                                                 {
   MoFEMFunctionBegin;
   CHKERR playW_NoHessian();
   CHKERR playY_NoGradient();
   CHKERR playH_NoHessian();
   MoFEMFunctionReturn(0);
 }

◆ playW()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playW ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 122 of file ClosetPointProjection.cpp.

                                              {
   MoFEMFunctionBegin;
   int r;
   int adloc_return_value = 0;
   r = ::function(tapesTags[ClosestPointProjection::W], 1,
                  activeVariablesW.size(), &activeVariablesW[0], &w);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   r = gradient(tapesTags[ClosestPointProjection::W], activeVariablesW.size(),
                &activeVariablesW[0], &gradientW[0]);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   hessianW.resize(activeVariablesW.size(), activeVariablesW.size(), false);
   hessianW.clear();
   vector<double *> hessian_w(activeVariablesW.size());
   for (unsigned int dd = 0; dd < activeVariablesW.size(); dd++) {
     hessian_w[dd] = &hessianW(dd, 0);
   }
   r = hessian(tapesTags[ClosestPointProjection::W], activeVariablesW.size(),
               &activeVariablesW[0], &hessian_w[0]);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   partialWStrainPlasticStrain.resize(6, 6, false);
   for (int ii = 0; ii < 6; ii++) {
     for (int jj = 0; jj < 6; jj++) {
       partialWStrainPlasticStrain(ii, jj) = hessianW(6 + ii, jj);
     }
   }
   C.resize(6, 6, false);
   for (int ii = 0; ii < 6; ii++) {
     for (int jj = 0; jj <= ii; jj++) {
       C(ii, jj) = hessianW(6 + ii, 6 + jj);
     }
   }
   D.resize(nbInternalVariables, nbInternalVariables, false);
   for (unsigned int ii = 0; ii < nbInternalVariables; ii++) {
     for (unsigned int jj = 0; jj <= ii; jj++) {
       D(ii, jj) = hessianW(12 + ii, 12 + jj);
     }
   }
   MoFEMFunctionReturn(0);
 }

◆ playW_NoHessian()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playW_NoHessian ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 172 of file ClosetPointProjection.cpp.

                                         {
   MoFEMFunctionBegin;
   int r;
   int adloc_return_value = 0;
   r = ::function(tapesTags[ClosestPointProjection::W], 1,
                  activeVariablesW.size(), &activeVariablesW[0], &w);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   r = gradient(tapesTags[ClosestPointProjection::W], activeVariablesW.size(),
                &activeVariablesW[0], &gradientW[0]);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   MoFEMFunctionReturn(0);
 }

◆ playY()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playY ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 191 of file ClosetPointProjection.cpp.

                                              {
   MoFEMFunctionBegin;
   auto active_variables_yh = getActiveVariablesYH();
   int r;
   int adloc_return_value = 0;
   r = ::function(tapesTags[ClosestPointProjection::Y], 1,
                  active_variables_yh.size(), &active_variables_yh[0], &y);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   gradientY.resize(active_variables_yh.size());
   r = gradient(tapesTags[ClosestPointProjection::Y], active_variables_yh.size(),
                &active_variables_yh[0], &gradientY[0]);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   partialYSigma =
       VectorAdaptor(6, ublas::shallow_array_adaptor<double>(6, &gradientY[0]));
   partialYFlux = VectorAdaptor(
       nbInternalVariables,
       ublas::shallow_array_adaptor<double>(nbInternalVariables, &gradientY[6]));
   MoFEMFunctionReturn(0);
 }

◆ playY_NoGradient()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playY_NoGradient ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 217 of file ClosetPointProjection.cpp.

                                          {
   MoFEMFunctionBegin;
   auto active_variables_yh = getActiveVariablesYH();
   int r;
   int adloc_return_value = 0;
   r = ::function(tapesTags[ClosestPointProjection::Y], 1,
                  active_variables_yh.size(), &active_variables_yh[0], &y);
   if (r < adloc_return_value) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
             "ADOL-C function evaluation with error");
   }
   MoFEMFunctionReturn(0);
 }

◆ recordH()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::recordH ( )

Record flow potential.

Examples: ADOLCPlasticity.hpp.

Definition at line 96 of file ClosetPointProjection.cpp.

                                                {
   MoFEMFunctionBegin;
  
   auto stress = getStress();
   auto internal_fluxes = getInternalFluxes();
  
   trace_on(tapesTags[ClosestPointProjection::H]);
   {
     // active variables
     a_sTress.resize(6, false);
     for (int ii = 0; ii < 6; ii++) {
       a_sTress[ii] <<= stress[ii];
     }
     a_internalFluxes.resize(nbInternalVariables, false);
     for (unsigned int ii = 0; ii < nbInternalVariables; ii++) {
       a_internalFluxes[ii] <<= internal_fluxes[ii];
     }
     // evaluete functions
     CHKERR flowPotential();
     // return variables
     a_h >>= h;
   }
   trace_off();
   MoFEMFunctionReturn(0);
 }

◆ recordTapes()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::recordTapes ( )

Record tapes.

Examples: ADOLCPlasticity.hpp.

Definition at line 610 of file ClosetPointProjection.cpp.

                                                    {
   MoFEMFunctionBegin;
   VectorDouble active(12 + nbInternalVariables);
   VectorDouble gradient(12 + nbInternalVariables);
   auto tmp_active = getVectorAdaptor(&(active[0]), 12 + nbInternalVariables);
   auto tmp_gradient =
       getVectorAdaptor(&(gradient[0]), 12 + nbInternalVariables);
   activeVariablesW.swap(tmp_active);
   gradientW.swap(tmp_gradient);
   CHKERR evaluatePotentials();
   MoFEMFunctionReturn(0);
 }

◆ recordW()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::recordW ( )

Record strain energy.

Examples: ADOLCPlasticity.hpp.

Definition at line 40 of file ClosetPointProjection.cpp.

                                                {
   MoFEMFunctionBegin;
  
   auto plastic_strain = getPlasticStrain();
   auto total_strain = getTotalStrain();
   auto internal_variables = getInternalVariables();
  
   trace_on(tapesTags[ClosestPointProjection::W]);
   {
     // active variables
     a_plasticStrain.resize(6, false);
     for (int ii = 0; ii < 6; ii++) {
       a_plasticStrain[ii] <<= plastic_strain[ii];
     }
     a_sTrain.resize(6, false);
     for (int ii = 0; ii < 6; ii++) {
       a_sTrain[ii] <<= total_strain[ii];
     }
     a_internalVariables.resize(nbInternalVariables, false);
     for (unsigned int ii = 0; ii < nbInternalVariables; ii++) {
       a_internalVariables[ii] <<= internal_variables[ii];
     }
     // evaluete functions
     CHKERR freeHemholtzEnergy();
     a_w >>= w;
   }
   trace_off();
   MoFEMFunctionReturn(0);
 }

◆ recordY()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::recordY ( )

Record yield function.

Examples: ADOLCPlasticity.hpp.

Definition at line 70 of file ClosetPointProjection.cpp.

                                                {
   MoFEMFunctionBegin;
  
   auto stress = getStress();
   auto internal_fluxes = getInternalFluxes();
  
   trace_on(tapesTags[ClosestPointProjection::Y]);
   {
     // active variables
     a_sTress.resize(6, false);
     for (int ii = 0; ii < 6; ii++) {
       a_sTress[ii] <<= stress[ii];
     }
     a_internalFluxes.resize(nbInternalVariables, false);
     for (unsigned int ii = 0; ii < nbInternalVariables; ii++) {
       a_internalFluxes[ii] <<= internal_fluxes[ii];
     }
     // evaluete functions
     CHKERR yieldFunction();
     // return variables
     a_y >>= y;
   }
   trace_off();
   MoFEMFunctionReturn(0);
 }

◆ setParams()

virtual MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::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 in ADOLCPlasticity::ParaboloidalPlasticity, and ADOLCPlasticity::J2Plasticity< 3 >.

Examples: ADOLCPlasticity.hpp.

Definition at line 274 of file ADOLCPlasticity.hpp.

                                                                       {
     return 0;
   }

◆ snesCreate()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::snesCreate ( )

Create nested snes.

Examples: ADOLCPlasticity.hpp.

Definition at line 433 of file ClosetPointProjection.cpp.

                                                   {
   MoFEMFunctionBegin;
   sNes = createSNES(PETSC_COMM_SELF);
   CHKERR SNESSetFunction(sNes, R, ADOLCPlasticityRes, (void *)this);
   CHKERR SNESSetJacobian(sNes, A, A, ADOLCPlasticityJac, (void *)this);
  
   KSP ksp;
   PC pc;
   SNESLineSearch line_search;
   CHKERR SNESGetKSP(sNes, &ksp);
   CHKERR KSPGetPC(ksp, &pc);
   CHKERR KSPSetType(ksp, KSPPREONLY);
   CHKERR PCSetType(pc, PCLU);
   CHKERR SNESSetTolerances(sNes, 1e-8, 1e-6, 1e-12, 20, 1000);
   CHKERR SNESGetLineSearch(sNes, &line_search);
   // CHKERR SNESLineSearchSetType(line_search, SNESLINESEARCHBASIC);
   CHKERR SNESLineSearchSetType(line_search, SNESLINESEARCHBT);
   // CHKERR SNESLineSearchSetType(line_search, SNESLINESEARCHL2);
   
   CHKERR SNESAppendOptionsPrefix(sNes, "cp_");
   CHKERR SNESSetFromOptions(sNes);
  
   MoFEMFunctionReturn(0);
 }

◆ solveClosetProjection()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::solveClosetProjection ( )

Solve nonlinear system of equations to find stress, internal fluxes, and Lagrange plastic multiplier.

Examples: ADOLCPlasticity.hpp.

Definition at line 459 of file ClosetPointProjection.cpp.

                                               {
   MoFEMFunctionBegin;
  
   auto plastic_strain = getPlasticStrain();
   auto internal_variables = getInternalVariables();
  
   {
     double *array;
     CHKERR VecGetArray(Chi, &array);
     for (int ii = 0; ii < 6; ii++) {
       array[ii] = plastic_strain[ii];
     }
     int nb_internal_variables = nbInternalVariables;
     for (unsigned int ii = 0; ii < nb_internal_variables; ii++) {
       array[6 + ii] = internal_variables[ii];
     }
     deltaGamma = 0;
     CHKERR VecRestoreArray(Chi, &array);
   }
  
   CHKERR SNESSolve(sNes, PETSC_NULL, Chi);
   SNESConvergedReason reason;
   CHKERR SNESGetConvergedReason(sNes, &reason);
   if (reason < 0) {
     int its;
     CHKERR SNESGetIterationNumber(sNes, &its);
     MOFEM_LOG_C(
         "ADOLCPlasticitySelf", Sev::warning,
         "Plasticity Closet Point Projection - number of Newton iterations = %d",
         its);
     plastic_strain = plasticStrain0;
     internal_variables = internalVariables0;
     deltaGamma = 0;
   } else {
     double *array;
     CHKERR VecGetArray(Chi, &array);
     for (int ii = 0; ii < 6; ii++) {
       plastic_strain[ii] = array[ii];
     }
     int nb_internal_variables = nbInternalVariables;
     for (unsigned int ii = 0; ii < nb_internal_variables; ii++) {
       internal_variables[ii] = array[6 + ii];
     }
     deltaGamma = array[6 + nb_internal_variables];
     CHKERR VecRestoreArray(Chi, &array);
   }
   MoFEMFunctionReturn(0);
 }

◆ yieldFunction()

virtual MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::yieldFunction ( )

pure virtual

Set yield function.

Implemented in ADOLCPlasticity::ParaboloidalPlasticity, and ADOLCPlasticity::J2Plasticity< 3 >.

Examples: ADOLCPlasticity.hpp.

Friends And Related Function Documentation

◆ ADOLCPlasticityJac

MoFEMErrorCode ADOLCPlasticityJac	(	SNES	snes,
		Vec	chi,
		Mat	A,
		Mat	,
		void *	ctx
	)

friend

Function executed by nested SENES at evaluationg tangent matrix.

Examples: ADOLCPlasticity.hpp.

Definition at line 584 of file ClosetPointProjection.cpp.

                                                     {
   MoFEMFunctionBegin;
   ClosestPointProjection *cp;
   cp = (ClosestPointProjection *)ctx;
  
   auto plastic_strain = cp->getPlasticStrain();
   auto internal_variables = cp->getInternalVariables();
  
   {
     const double *array;
     CHKERR VecGetArrayRead(chi, &array);
     for (auto ii = 0; ii < 6; ii++) {
       plastic_strain[ii] = array[ii];
     }
     for (auto ii = 0; ii < cp->nbInternalVariables; ii++) {
       internal_variables[ii] = array[6 + ii];
     }
     cp->deltaGamma = array[6 + cp->nbInternalVariables];
     CHKERR VecRestoreArrayRead(chi, &array);
   }
   CHKERR cp->playPotentials();
   CHKERR cp->calculateA();
   MoFEMFunctionReturn(0);
 }

◆ ADOLCPlasticityRes

MoFEMErrorCode ADOLCPlasticityRes	(	SNES	snes,
		Vec	chi,
		Vec	r,
		void *	ctx
	)

friend

Function executed by nested SENES at evaluationg residual.

Examples: ADOLCPlasticity.hpp.

Definition at line 558 of file ClosetPointProjection.cpp.

                                                                         {
   MoFEMFunctionBegin;
   ClosestPointProjection *cp;
   cp = (ClosestPointProjection *)ctx;
  
   auto plastic_strain = cp->getPlasticStrain();
   auto internal_variables = cp->getInternalVariables();
  
   {
     const double *array;
     CHKERR VecGetArrayRead(chi, &array);
     for (auto ii = 0; ii < 6; ii++) {
       plastic_strain[ii] = array[ii];
     }
     for (auto ii = 0; ii < cp->nbInternalVariables; ii++) {
       internal_variables[ii] = array[6 + ii];
     }
     cp->deltaGamma = array[6 + cp->nbInternalVariables];
     CHKERR VecRestoreArrayRead(chi, &array);
   }
   CHKERR cp->playPotentials_NoHessian();
   CHKERR cp->calculateR(r);
  
   MoFEMFunctionReturn(0);
 }

Member Data Documentation

◆ A

SmartPetscObj<Mat> ADOLCPlasticity::ClosestPointProjection::A

Examples: ADOLCPlasticity.hpp.

Definition at line 317 of file ADOLCPlasticity.hpp.

◆ a_h

adouble ADOLCPlasticity::ClosestPointProjection::a_h

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 330 of file ADOLCPlasticity.hpp.

◆ a_internalFluxes

ublas::vector<adouble> ADOLCPlasticity::ClosestPointProjection::a_internalFluxes

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 327 of file ADOLCPlasticity.hpp.

◆ a_internalVariables

ublas::vector<adouble> ADOLCPlasticity::ClosestPointProjection::a_internalVariables

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 326 of file ADOLCPlasticity.hpp.

◆ a_plasticStrain

ublas::vector<adouble> ADOLCPlasticity::ClosestPointProjection::a_plasticStrain

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 325 of file ADOLCPlasticity.hpp.

◆ a_sTrain

ublas::vector<adouble> ADOLCPlasticity::ClosestPointProjection::a_sTrain

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 324 of file ADOLCPlasticity.hpp.

◆ a_sTress

ublas::vector<adouble> ADOLCPlasticity::ClosestPointProjection::a_sTress

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 327 of file ADOLCPlasticity.hpp.

◆ a_w

adouble ADOLCPlasticity::ClosestPointProjection::a_w

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 328 of file ADOLCPlasticity.hpp.

◆ a_y

adouble ADOLCPlasticity::ClosestPointProjection::a_y

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 329 of file ADOLCPlasticity.hpp.

◆ activeVariablesW

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::activeVariablesW

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 173 of file ADOLCPlasticity.hpp.

◆ C

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

Examples: ADOLCPlasticity.hpp.

Definition at line 182 of file ADOLCPlasticity.hpp.

◆ Cep

MatrixDouble ADOLCPlasticity::ClosestPointProjection::Cep

Examples: ADOLCPlasticity.hpp.

Definition at line 181 of file ADOLCPlasticity.hpp.

◆ Chi

SmartPetscObj<Vec> ADOLCPlasticity::ClosestPointProjection::Chi

Examples: ADOLCPlasticity.hpp.

Definition at line 319 of file ADOLCPlasticity.hpp.

◆ Cp

MatrixDouble ADOLCPlasticity::ClosestPointProjection::Cp

Examples: ADOLCPlasticity.hpp.

Definition at line 181 of file ADOLCPlasticity.hpp.

◆ D

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

Examples: ADOLCPlasticity.hpp.

Definition at line 182 of file ADOLCPlasticity.hpp.

◆ dataA

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

Examples: ADOLCPlasticity.hpp.

Definition at line 321 of file ADOLCPlasticity.hpp.

◆ dChi

SmartPetscObj<Vec> ADOLCPlasticity::ClosestPointProjection::dChi

Examples: ADOLCPlasticity.hpp.

Definition at line 320 of file ADOLCPlasticity.hpp.

◆ deltaGamma

double ADOLCPlasticity::ClosestPointProjection::deltaGamma

Examples: ADOLCPlasticity.hpp.

Definition at line 180 of file ADOLCPlasticity.hpp.

◆ Ep

MatrixDouble ADOLCPlasticity::ClosestPointProjection::Ep

Examples: ADOLCPlasticity.hpp.

Definition at line 181 of file ADOLCPlasticity.hpp.

◆ gradientH

VectorDouble ADOLCPlasticity::ClosestPointProjection::gradientH

Examples: ADOLCPlasticity.hpp.

Definition at line 207 of file ADOLCPlasticity.hpp.

◆ gradientW

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::gradientW

Examples: ADOLCPlasticity.hpp.

Definition at line 174 of file ADOLCPlasticity.hpp.

◆ gradientY

VectorDouble ADOLCPlasticity::ClosestPointProjection::gradientY

Examples: ADOLCPlasticity.hpp.

Definition at line 203 of file ADOLCPlasticity.hpp.

◆ h

double ADOLCPlasticity::ClosestPointProjection::h

Examples: ADOLCPlasticity.hpp.

Definition at line 178 of file ADOLCPlasticity.hpp.

◆ hessianH

MatrixDouble ADOLCPlasticity::ClosestPointProjection::hessianH

Examples: ADOLCPlasticity.hpp.

Definition at line 208 of file ADOLCPlasticity.hpp.

◆ hessianW

MatrixDouble ADOLCPlasticity::ClosestPointProjection::hessianW

Examples: ADOLCPlasticity.hpp.

Definition at line 199 of file ADOLCPlasticity.hpp.

◆ integrationRule

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

Initial value:

= [](int, int, int p) {
    return 2 * (p - 1);
  }

Examples: ADOLCPlasticity.hpp.

Definition at line 142 of file ADOLCPlasticity.hpp.

◆ internalVariables0

VectorDouble ADOLCPlasticity::ClosestPointProjection::internalVariables0

Examples: ADOLCPlasticity.hpp.

Definition at line 146 of file ADOLCPlasticity.hpp.

◆ nbInternalVariables

int ADOLCPlasticity::ClosestPointProjection::nbInternalVariables

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 141 of file ADOLCPlasticity.hpp.

◆ partial2HFlux

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

Second partial derivative of flow potential with respect to internal fluxes

Examples: ADOLCPlasticity.hpp.

Definition at line 312 of file ADOLCPlasticity.hpp.

◆ partial2HSigma

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

Second partial derivative of flow potential with respect to stresses and internal

Examples: ADOLCPlasticity.hpp.

Definition at line 309 of file ADOLCPlasticity.hpp.

◆ partial2HSigmaFlux

MatrixDouble ADOLCPlasticity::ClosestPointProjection::partial2HSigmaFlux

Mixed decond partial derivative of flow potential with respect to stresses and internal fluxes

Examples: ADOLCPlasticity.hpp.

Definition at line 315 of file ADOLCPlasticity.hpp.

◆ partialHFlux

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::partialHFlux

Examples: ADOLCPlasticity.hpp.

Definition at line 304 of file ADOLCPlasticity.hpp.

◆ partialHSigma

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::partialHSigma

Examples: ADOLCPlasticity.hpp.

Definition at line 302 of file ADOLCPlasticity.hpp.

◆ partialWStrainPlasticStrain

MatrixDouble ADOLCPlasticity::ClosestPointProjection::partialWStrainPlasticStrain

Examples: ADOLCPlasticity.hpp.

Definition at line 296 of file ADOLCPlasticity.hpp.

◆ partialYFlux

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::partialYFlux

Examples: ADOLCPlasticity.hpp.

Definition at line 300 of file ADOLCPlasticity.hpp.

◆ partialYSigma

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::partialYSigma

Examples: ADOLCPlasticity.hpp.

Definition at line 298 of file ADOLCPlasticity.hpp.

◆ plasticStrain0

VectorDouble ADOLCPlasticity::ClosestPointProjection::plasticStrain0

Examples: ADOLCPlasticity.hpp.

Definition at line 147 of file ADOLCPlasticity.hpp.

◆ R

SmartPetscObj<Vec> ADOLCPlasticity::ClosestPointProjection::R

Examples: ADOLCPlasticity.hpp.

Definition at line 318 of file ADOLCPlasticity.hpp.

◆ sNes

SmartPetscObj<SNES> ADOLCPlasticity::ClosestPointProjection::sNes

Examples: ADOLCPlasticity.hpp.

Definition at line 322 of file ADOLCPlasticity.hpp.

◆ tapesTags

std::array<int, LAST_TAPE> ADOLCPlasticity::ClosestPointProjection::tapesTags

Examples: ADOLCPlasticity.hpp.

Definition at line 169 of file ADOLCPlasticity.hpp.

◆ w

double ADOLCPlasticity::ClosestPointProjection::w

Examples: ADOLCPlasticity.hpp.

Definition at line 176 of file ADOLCPlasticity.hpp.

◆ y

double ADOLCPlasticity::ClosestPointProjection::y

Examples: ADOLCPlasticity.hpp.

Definition at line 177 of file ADOLCPlasticity.hpp.

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

users_modules/adolc-plasticity/src/ADOLCPlasticity.hpp
users_modules/adolc-plasticity/src/impl/ClosetPointProjection.cpp

Public Types

Public Member Functions

Public Attributes

Friends

Detailed Description

Member Enumeration Documentation

◆ TypesTags

Constructor & Destructor Documentation

◆ ClosestPointProjection()

Member Function Documentation

◆ addMatBlockOps()

◆ calculateA()

◆ calculateR()

◆ consistentTangent()

◆ createMatAVecR()

◆ evaluatePotentials()

◆ flowPotential()

◆ freeHemholtzEnergy()

◆ getActiveVariablesYH()

◆ getInternalFluxes()

◆ getInternalVariables()

◆ getPlasticStrain()

◆ getStress()

◆ getTotalStrain()

◆ playH()

◆ playH_NoHessian()

◆ playPotentials()

◆ playPotentials_NoHessian()

◆ playW()

◆ playW_NoHessian()

◆ playY()

◆ playY_NoGradient()

◆ recordH()

◆ recordTapes()

◆ recordW()

◆ recordY()

◆ setParams()

◆ snesCreate()

◆ solveClosetProjection()

◆ yieldFunction()

Friends And Related Function Documentation

◆ ADOLCPlasticityJac

◆ ADOLCPlasticityRes

Member Data Documentation

◆ A

◆ a_h

◆ a_internalFluxes

◆ a_internalVariables

◆ a_plasticStrain

◆ a_sTrain

◆ a_sTress

◆ a_w

◆ a_y

◆ activeVariablesW

◆ C

◆ Cep

◆ Chi

◆ Cp

◆ D

◆ dataA

◆ dChi

◆ deltaGamma

◆ Ep

◆ gradientH

◆ gradientW

◆ gradientY

◆ h

◆ hessianH

◆ hessianW

◆ integrationRule

◆ internalVariables0

◆ nbInternalVariables

◆ partial2HFlux

◆ partial2HSigma

◆ partial2HSigmaFlux

◆ partialHFlux

◆ partialHSigma

◆ partialWStrainPlasticStrain

◆ partialYFlux

◆ partialYSigma