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.

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.

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

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

virtual MoFEMErrorCode	freeHemholtzEnergy ()=0
	Set Hemholtz energy.

virtual MoFEMErrorCode	yieldFunction ()=0
	Set yield function.

virtual MoFEMErrorCode	flowPotential ()=0
	Set flow potential.

virtual MoFEMErrorCode	codedHessianW (vector< double * >)

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

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

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

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

Detailed Description

Closest point projection algorithm.

Definition at line 145 of file ADOLCPlasticity.hpp.

Member Enumeration Documentation

◆ TypesTags

enum ADOLCPlasticity::ClosestPointProjection::TypesTags

Enumerator
W
Y
H
LAST_TAPE

Examples: ADOLCPlasticity.hpp.

Definition at line 174 of file ADOLCPlasticity.hpp.

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

ADOLCPlasticity::ClosestPointProjection::H

@ H

Definition ADOLCPlasticity.hpp:174

ADOLCPlasticity::ClosestPointProjection::W

@ W

Definition ADOLCPlasticity.hpp:174

ADOLCPlasticity::ClosestPointProjection::Y

@ Y

Definition ADOLCPlasticity.hpp:174

ADOLCPlasticity::ClosestPointProjection::LAST_TAPE

@ LAST_TAPE

Definition ADOLCPlasticity.hpp:174

Constructor & Destructor Documentation

◆ ClosestPointProjection()

ADOLCPlasticity::ClosestPointProjection::ClosestPointProjection ( )

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 15 of file ClosestPointProjection.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
 
  implHessianW = PETSC_FALSE;
  CHKERR PetscOptionsGetBool(PETSC_NULLPTR, "-impl_hessianW", &implHessianW, 0);
 
}

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::HeterogeneousParaboloidalPlasticity, ADOLCPlasticity::J2Plasticity< 3 >, and ADOLCPlasticity::ParaboloidalPlasticity.

Examples: ADOLCPlasticity.hpp.

Definition at line 269 of file ADOLCPlasticity.hpp.

                                                {
    return 0;
  }

◆ calculateA()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::calculateA ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 393 of file ClosestPointProjection.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 371 of file ClosestPointProjection.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);
}

◆ codedHessianW()

virtual MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::codedHessianW ( vector< double * > )

inlinevirtual

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

Examples: ADOLCPlasticity.hpp.

Definition at line 302 of file ADOLCPlasticity.hpp.

                                                         {
    return 0;
  }

◆ consistentTangent()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::consistentTangent ( )

Calculate consistent tangent matrix.

Examples: ADOLCPlasticity.hpp.

Definition at line 520 of file ClosestPointProjection.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 325 of file ClosestPointProjection.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 346 of file ClosestPointProjection.cpp.

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

◆ flowPotential()

virtual MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::flowPotential ( )

pure virtual

Set flow potential.

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

Examples: ADOLCPlasticity.hpp.

◆ freeHemholtzEnergy()

virtual MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::freeHemholtzEnergy ( )

pure virtual

Set Hemholtz energy.

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

Examples: ADOLCPlasticity.hpp.

◆ getActiveVariablesYH()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getActiveVariablesYH ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 164 of file ADOLCPlasticity.hpp.

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

◆ getInternalFluxes()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getInternalFluxes ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 170 of file ADOLCPlasticity.hpp.

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

◆ getInternalVariables()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getInternalVariables ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 161 of file ADOLCPlasticity.hpp.

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

◆ getPlasticStrain()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getPlasticStrain ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 155 of file ADOLCPlasticity.hpp.

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

◆ getStress()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getStress ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 167 of file ADOLCPlasticity.hpp.

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

◆ getTotalStrain()

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::getTotalStrain ( )

inline

Examples: ADOLCPlasticity.hpp.

Definition at line 158 of file ADOLCPlasticity.hpp.

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

◆ playH()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playH ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 241 of file ClosestPointProjection.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 298 of file ClosestPointProjection.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 354 of file ClosestPointProjection.cpp.

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

◆ playPotentials_NoHessian()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playPotentials_NoHessian ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 362 of file ClosestPointProjection.cpp.

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

◆ playW()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::playW ( )

Examples: ADOLCPlasticity.hpp.

Definition at line 125 of file ClosestPointProjection.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);
  }
  if (!implHessianW) {
    // call adolc hessian (energy, active variables)
    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");
    }
  } else {
    // call hard coded hessian depending on material model
    CHKERR codedHessianW(hessian_w);
  }
  
  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 183 of file ClosestPointProjection.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 202 of file ClosestPointProjection.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 228 of file ClosestPointProjection.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 99 of file ClosestPointProjection.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 621 of file ClosestPointProjection.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 43 of file ClosestPointProjection.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 73 of file ClosestPointProjection.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::J2Plasticity< 3 >, and ADOLCPlasticity::ParaboloidalPlasticity.

Examples: ADOLCPlasticity.hpp.

Definition at line 282 of file ADOLCPlasticity.hpp.

                                                                      {
    return 0;
  }

◆ snesCreate()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::snesCreate ( )

Create nested snes.

Examples: ADOLCPlasticity.hpp.

Definition at line 444 of file ClosestPointProjection.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);
}

◆ solveClosestProjection()

MoFEMErrorCode ADOLCPlasticity::ClosestPointProjection::solveClosestProjection ( )

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

Examples: ADOLCPlasticity.hpp.

Definition at line 470 of file ClosestPointProjection.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_NULLPTR, Chi);
  SNESConvergedReason reason;
  CHKERR SNESGetConvergedReason(sNes, &reason);
  if (reason < 0) {
    int its;
    CHKERR SNESGetIterationNumber(sNes, &its);
    MOFEM_LOG_C(
        "ADOLCPlasticitySelf", Sev::warning,
        "Plasticity Closest 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::J2Plasticity< 3 >, and ADOLCPlasticity::ParaboloidalPlasticity.

Examples: ADOLCPlasticity.hpp.

Friends And Related Symbol 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 595 of file ClosestPointProjection.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 569 of file ClosestPointProjection.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 329 of file ADOLCPlasticity.hpp.

◆ a_h

adouble ADOLCPlasticity::ClosestPointProjection::a_h

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 342 of file ADOLCPlasticity.hpp.

◆ a_internalFluxes

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

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 339 of file ADOLCPlasticity.hpp.

◆ a_internalVariables

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

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 338 of file ADOLCPlasticity.hpp.

◆ a_plasticStrain

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

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 337 of file ADOLCPlasticity.hpp.

◆ a_sTrain

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

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 336 of file ADOLCPlasticity.hpp.

◆ a_sTress

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

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 339 of file ADOLCPlasticity.hpp.

◆ a_w

adouble ADOLCPlasticity::ClosestPointProjection::a_w

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 340 of file ADOLCPlasticity.hpp.

◆ a_y

adouble ADOLCPlasticity::ClosestPointProjection::a_y

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 341 of file ADOLCPlasticity.hpp.

◆ activeVariablesW

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::activeVariablesW

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 179 of file ADOLCPlasticity.hpp.

◆ C

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

Examples: ADOLCPlasticity.hpp.

Definition at line 188 of file ADOLCPlasticity.hpp.

◆ Cep

MatrixDouble ADOLCPlasticity::ClosestPointProjection::Cep

Examples: ADOLCPlasticity.hpp.

Definition at line 187 of file ADOLCPlasticity.hpp.

◆ Chi

SmartPetscObj<Vec> ADOLCPlasticity::ClosestPointProjection::Chi

Examples: ADOLCPlasticity.hpp.

Definition at line 331 of file ADOLCPlasticity.hpp.

◆ Cp

MatrixDouble ADOLCPlasticity::ClosestPointProjection::Cp

Examples: ADOLCPlasticity.hpp.

Definition at line 187 of file ADOLCPlasticity.hpp.

◆ D

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

Examples: ADOLCPlasticity.hpp.

Definition at line 188 of file ADOLCPlasticity.hpp.

◆ dataA

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

Examples: ADOLCPlasticity.hpp.

Definition at line 333 of file ADOLCPlasticity.hpp.

◆ dChi

SmartPetscObj<Vec> ADOLCPlasticity::ClosestPointProjection::dChi

Examples: ADOLCPlasticity.hpp.

Definition at line 332 of file ADOLCPlasticity.hpp.

◆ deltaGamma

double ADOLCPlasticity::ClosestPointProjection::deltaGamma

Examples: ADOLCPlasticity.hpp.

Definition at line 186 of file ADOLCPlasticity.hpp.

◆ Ep

MatrixDouble ADOLCPlasticity::ClosestPointProjection::Ep

Examples: ADOLCPlasticity.hpp.

Definition at line 187 of file ADOLCPlasticity.hpp.

◆ gradientH

VectorDouble ADOLCPlasticity::ClosestPointProjection::gradientH

Examples: ADOLCPlasticity.hpp.

Definition at line 215 of file ADOLCPlasticity.hpp.

◆ gradientW

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::gradientW

Examples: ADOLCPlasticity.hpp.

Definition at line 180 of file ADOLCPlasticity.hpp.

◆ gradientY

VectorDouble ADOLCPlasticity::ClosestPointProjection::gradientY

Examples: ADOLCPlasticity.hpp.

Definition at line 211 of file ADOLCPlasticity.hpp.

◆ h

double ADOLCPlasticity::ClosestPointProjection::h

Examples: ADOLCPlasticity.hpp.

Definition at line 184 of file ADOLCPlasticity.hpp.

◆ hessianH

MatrixDouble ADOLCPlasticity::ClosestPointProjection::hessianH

Examples: ADOLCPlasticity.hpp.

Definition at line 216 of file ADOLCPlasticity.hpp.

◆ hessianW

MatrixDouble ADOLCPlasticity::ClosestPointProjection::hessianW

Examples: ADOLCPlasticity.hpp.

Definition at line 207 of file ADOLCPlasticity.hpp.

◆ implHessianW

PetscBool ADOLCPlasticity::ClosestPointProjection::implHessianW

Examples: ADOLCPlasticity.hpp.

Definition at line 190 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 148 of file ADOLCPlasticity.hpp.

                                                                          {
    return 2 * (p - 1);
  };

◆ internalVariables0

VectorDouble ADOLCPlasticity::ClosestPointProjection::internalVariables0

Examples: ADOLCPlasticity.hpp.

Definition at line 152 of file ADOLCPlasticity.hpp.

◆ nbInternalVariables

int ADOLCPlasticity::ClosestPointProjection::nbInternalVariables

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 147 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 324 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 321 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 327 of file ADOLCPlasticity.hpp.

◆ partialHFlux

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::partialHFlux

Examples: ADOLCPlasticity.hpp.

Definition at line 316 of file ADOLCPlasticity.hpp.

◆ partialHSigma

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::partialHSigma

Examples: ADOLCPlasticity.hpp.

Definition at line 314 of file ADOLCPlasticity.hpp.

◆ partialWStrainPlasticStrain

MatrixDouble ADOLCPlasticity::ClosestPointProjection::partialWStrainPlasticStrain

Examples: ADOLCPlasticity.hpp.

Definition at line 308 of file ADOLCPlasticity.hpp.

◆ partialYFlux

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::partialYFlux

Examples: ADOLCPlasticity.hpp.

Definition at line 312 of file ADOLCPlasticity.hpp.

◆ partialYSigma

VectorAdaptor ADOLCPlasticity::ClosestPointProjection::partialYSigma

Examples: ADOLCPlasticity.hpp.

Definition at line 310 of file ADOLCPlasticity.hpp.

◆ plasticStrain0

VectorDouble ADOLCPlasticity::ClosestPointProjection::plasticStrain0

Examples: ADOLCPlasticity.hpp.

Definition at line 153 of file ADOLCPlasticity.hpp.

◆ R

SmartPetscObj<Vec> ADOLCPlasticity::ClosestPointProjection::R

Examples: ADOLCPlasticity.hpp.

Definition at line 330 of file ADOLCPlasticity.hpp.

◆ sNes

SmartPetscObj<SNES> ADOLCPlasticity::ClosestPointProjection::sNes

Examples: ADOLCPlasticity.hpp.

Definition at line 334 of file ADOLCPlasticity.hpp.

◆ tapesTags

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

Examples: ADOLCPlasticity.hpp, and ADOLCPlasticityMaterialModels.hpp.

Definition at line 175 of file ADOLCPlasticity.hpp.

◆ w

double ADOLCPlasticity::ClosestPointProjection::w

Examples: ADOLCPlasticity.hpp.

Definition at line 182 of file ADOLCPlasticity.hpp.

◆ y

double ADOLCPlasticity::ClosestPointProjection::y

Examples: ADOLCPlasticity.hpp.

Definition at line 183 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/ClosestPointProjection.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()

◆ codedHessianW()

◆ 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()

◆ solveClosestProjection()

◆ yieldFunction()

Friends And Related Symbol 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

◆ implHessianW

◆ integrationRule

◆ internalVariables0

◆ nbInternalVariables

◆ partial2HFlux

◆ partial2HSigma

◆ partial2HSigmaFlux

◆ partialHFlux

◆ partialHSigma

◆ partialWStrainPlasticStrain