SmallStrainPlasticity::ClosestPointProjection Struct Reference

Closest point projection algorithm. More...

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

Inheritance diagram for SmallStrainPlasticity::ClosestPointProjection:

Collaboration diagram for SmallStrainPlasticity::ClosestPointProjection:

Public Member Functions
	ClosestPointProjection ()
virtual PetscErrorCode	setActiveVariablesW ()
virtual PetscErrorCode	setActiveVariablesYH ()
virtual PetscErrorCode	rEcordW ()
virtual PetscErrorCode	rEcordY ()
virtual PetscErrorCode	rEcordH ()
virtual PetscErrorCode	pLayW ()
virtual PetscErrorCode	pLayW_NoHessian ()
virtual PetscErrorCode	pLayY ()
virtual PetscErrorCode	pLayY_NoGradient ()
virtual PetscErrorCode	pLayH ()
virtual PetscErrorCode	createMatAVecR ()
virtual PetscErrorCode	destroyMatAVecR ()
virtual PetscErrorCode	evaluatePotentials ()
virtual PetscErrorCode	cAlculateR (Vec R)
virtual PetscErrorCode	cAlculateA ()
PetscErrorCode	snesCreate ()
PetscErrorCode	snesDestroy ()
PetscErrorCode	solveColasetProjection ()
PetscErrorCode	consistentTangent ()
virtual PetscErrorCode	freeHemholtzEnergy ()
virtual PetscErrorCode	yieldFunction ()
virtual PetscErrorCode	flowPotential ()

Public Attributes
VectorDouble	sTrain
VectorDouble	internalVariables
VectorDouble	internalVariables0
VectorDouble	plasticStrain0
VectorDouble	plasticStrain
double	deltaGamma
double	tOl
int	gG
VectorAdaptor	sTress
VectorAdaptor	internalFluxes
ublas::symmetric_matrix< double, ublas::lower >	C
ublas::symmetric_matrix< double, ublas::lower >	D
MatrixDouble	partialWStrainPlasticStrain
VectorAdaptor	partialYSigma
VectorAdaptor	partialYFlux
VectorAdaptor	partialHSigma
VectorAdaptor	partialHFlux
ublas::symmetric_matrix< double, ublas::lower >	partial2HSigma
ublas::symmetric_matrix< double, ublas::lower >	partial2HFlux
MatrixDouble	partial2HSigmaFlux
vector< int >	tAgs
VectorDouble	activeVariablesW
VectorDouble	activeVariablesYH
double	w
double	y
double	h
ublas::vector< adouble >	a_sTrain
ublas::vector< adouble >	a_plasticStrain
ublas::vector< adouble >	a_internalVariables
ublas::vector< adouble >	a_sTress
ublas::vector< adouble >	a_internalFluxes
adouble	a_w
adouble	a_y
adouble	a_h
VectorDouble	gradientW
MatrixDouble	hessianW
VectorDouble	gradientY
VectorDouble	gradientH
MatrixDouble	hessianH
Mat	A
ublas::matrix< double, ublas::column_major >	dataA
Vec	R
Vec	Chi
Vec	dChi
SNES	sNes
KSP	kSp
PC	pC
SNESLineSearch	lineSearch
MatrixDouble	Ep
MatrixDouble	Lp
MatrixDouble	Cp
MatrixDouble	Cep
VectorDouble	Yp

Friends
PetscErrorCode	SmallStrainPlasticityfRes (SNES, Vec, Vec, void *ctx)
PetscErrorCode	SmallStrainPlasticityfJac (SNES, Vec, Mat, Mat, void *ctx)

Detailed Description

Closest point projection algorithm.

Examples

SmallStrainPlasticityMaterialModels.hpp.

Definition at line 482 of file SmallStrainPlasticity.hpp.

Constructor & Destructor Documentation

ClosestPointProjection()

SmallStrainPlasticity::ClosestPointProjection::ClosestPointProjection ( )

inline

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 508 of file SmallStrainPlasticity.hpp.

{}

Member Function Documentation

cAlculateA()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::cAlculateA ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 978 of file SmallStrainPlasticity.cpp.

                                                                       {
  PetscFunctionBegin;
  // matrix A
  ierr = MatZeroEntries(A); CHKERRQ(ierr);
  // row 0
  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<internalVariables.size();jj++) {
      dataA(ii,6+jj) = deltaGamma*a01(ii,jj);
    }
    dataA(ii,ii) -= 1;
  }
  for(int ii = 0;ii<6;ii++) {
    dataA(ii,6+internalVariables.size()) = partialHSigma[ii];
  }
  // row 1
  MatrixDouble a11 = prod(partial2HFlux,D);
  MatrixDouble a10 = prod(trans(partial2HSigmaFlux),partialWStrainPlasticStrain);
  for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
    for(int jj = 0;jj<6;jj++) {
      dataA(6+ii,jj) += deltaGamma*a10(ii,jj);
    }
    for(unsigned int jj = 0;jj<internalVariables.size();jj++) {
      dataA(6+ii,6+jj) += deltaGamma*a11(ii,jj);
    }
    dataA(6+ii,6+ii) -= 1;
  }
  for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
    dataA(6+ii,6+internalVariables.size()) = partialHFlux[ii];
  }
  // row 3
  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+internalVariables.size(),dd) = partial_y_sigma_c[dd];
  }
  for(unsigned int dd = 0;dd<partial_y_flux_d.size();dd++) {
    dataA(6+internalVariables.size(),6+dd) = partial_y_flux_d[dd];
  }
  dataA(6+internalVariables.size(),6+internalVariables.size()) = 0;
  ierr = MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY); CHKERRQ(ierr);
  ierr = MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

cAlculateR()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::cAlculateR ( Vec  R )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 962 of file SmallStrainPlasticity.cpp.

                                                                            {
  PetscFunctionBegin;
  // Residual
  double *array;
  ierr = VecGetArray(R,&array); CHKERRQ(ierr);
  for(int ii = 0;ii<6;ii++) {
    array[ii] = -plasticStrain[ii] + plasticStrain0[ii] + deltaGamma*partialHSigma[ii];
  }
  for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
    array[6+ii] = -internalVariables[ii] + internalVariables0[ii] + deltaGamma*partialHFlux[ii];
  }
  array[6+internalVariables.size()] = y;
  ierr = VecRestoreArray(R,&array); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

consistentTangent()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::consistentTangent

(

)

Examples

SmallStrainPlasticity.hpp.

Definition at line 1095 of file SmallStrainPlasticity.cpp.

                                                                              {
  PetscFunctionBegin;
  ierr = evaluatePotentials(); CHKERRQ(ierr);
  ierr = cAlculateA(); CHKERRQ(ierr);
  Ep.resize(6,6,false);
  Lp.resize(internalVariables.size(),6,false);
  Yp.resize(6,false);
  MatrixDouble ep_row;
  ep_row = deltaGamma*prod(partial2HSigma,C);
  MatrixDouble alpha_row;
  alpha_row = deltaGamma*prod(trans(partial2HSigmaFlux),C);
  VectorDouble y_row;
  y_row = prod(trans(C),partialYSigma);
  const int n = 6+internalVariables.size();
  for(int dd = 0;dd<6;dd++) {
    ierr = VecZeroEntries(R); CHKERRQ(ierr);
    double *array;
    ierr = VecGetArray(R,&array); CHKERRQ(ierr);
    {
      for(int ii = 0;ii<6;ii++) {
        array[ii] = ep_row(ii,dd);
      }
      for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
        array[6+ii] = alpha_row(ii,dd);
      }
      array[n] = y_row[dd];
    }
    ierr = VecRestoreArray(R,&array); CHKERRQ(ierr);
    ierr = VecAssemblyBegin(R); CHKERRQ(ierr);
    ierr = VecAssemblyEnd(R); CHKERRQ(ierr);
    ierr = KSPSolve(kSp,R,dChi); CHKERRQ(ierr);
    ierr = VecGetArray(dChi,&array); CHKERRQ(ierr);
    {
      for(int ii = 0;ii<6;ii++) {
        Ep(ii,dd) = array[ii];
      }
      for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
        Lp(ii,dd) = array[6+ii];
      }
      Yp[dd] = array[n];
    }
    ierr = VecRestoreArray(dChi,&array); CHKERRQ(ierr);
  }
  Cp.resize(6,6,false);
  noalias(Cp) = prod(C,Ep);
  Cep.resize(6,6,false);
  noalias(Cep) = C+Cp;
  // cout << endl;
  // cout << "A " << dataA << endl;
  // cout << "C " << C << endl;
  // cout << "Ep " << Ep << endl;
  // // cout << "Lp " << Lp << endl;
  // // cout << "Yp " << Yp << endl;
  // cout << "Cp " << Cp << endl;
  // cout << "Cep " << Cep << endl;
  PetscFunctionReturn(0);
}

createMatAVecR()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::createMatAVecR ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 924 of file SmallStrainPlasticity.cpp.

                                                                           {
  PetscFunctionBegin;
  PetscInt n;
  n = 1+6+internalVariables.size();
  dataA.resize(n,n,false);
  ierr = MatCreateSeqDense(PETSC_COMM_SELF,n,n,&dataA(0,0),&A); CHKERRQ(ierr);
  ierr = VecCreateSeq(PETSC_COMM_SELF,n,&R); CHKERRQ(ierr);
  ierr = VecCreateSeq(PETSC_COMM_SELF,n,&Chi); CHKERRQ(ierr);
  ierr = VecCreateSeq(PETSC_COMM_SELF,n,&dChi); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

destroyMatAVecR()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::destroyMatAVecR ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 936 of file SmallStrainPlasticity.cpp.

                                                                            {
  PetscFunctionBegin;
  ierr = MatDestroy(&A); CHKERRQ(ierr);
  ierr = VecDestroy(&R); CHKERRQ(ierr);
  ierr = VecDestroy(&Chi); CHKERRQ(ierr);
  ierr = VecDestroy(&dChi); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

evaluatePotentials()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::evaluatePotentials ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 945 of file SmallStrainPlasticity.cpp.

                                                                               {
  PetscFunctionBegin;
  if(gG == 0) {
    ierr = rEcordW(); CHKERRQ(ierr);
  }
  ierr = setActiveVariablesW(); CHKERRQ(ierr);
  ierr = pLayW(); CHKERRQ(ierr);
  if(gG == 0) {
    ierr = rEcordY(); CHKERRQ(ierr);
    ierr = rEcordH(); CHKERRQ(ierr);
  }
  ierr = setActiveVariablesYH(); CHKERRQ(ierr);
  ierr = pLayY(); CHKERRQ(ierr);
  ierr = pLayH(); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

flowPotential()

virtual PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::flowPotential ( )

inlinevirtual

Reimplemented in SmallStrainParaboloidalPlasticity, and SmallStrainJ2Plasticity.

Examples

SmallStrainPlasticity.hpp.

Definition at line 591 of file SmallStrainPlasticity.hpp.

                                           {
      PetscFunctionBegin;
      SETERRQ(
        PETSC_COMM_SELF,
        MOFEM_NOT_IMPLEMENTED,
        "Not implemented (overload class and implement this function)"
      );
      PetscFunctionReturn(0);
    }

freeHemholtzEnergy()

virtual PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::freeHemholtzEnergy ( )

inlinevirtual

Reimplemented in SmallStrainParaboloidalPlasticity, and SmallStrainJ2Plasticity.

Examples

SmallStrainPlasticity.hpp.

Definition at line 571 of file SmallStrainPlasticity.hpp.

                                                {
      PetscFunctionBegin;
      SETERRQ(
        PETSC_COMM_SELF,
        MOFEM_NOT_IMPLEMENTED,
        "Not implemented (overload class and implement this function)"
      );
      PetscFunctionReturn(0);
    }

pLayH()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::pLayH ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 870 of file SmallStrainPlasticity.cpp.

                                                                  {
  PetscFunctionBegin;
  int r;
  int adloc_return_value = 0;
  r = ::function(tAgs[2],1,activeVariablesYH.size(),&activeVariablesYH[0],&h);
  if(r<adloc_return_value) {
    SETERRQ(PETSC_COMM_SELF,MOFEM_OPERATION_UNSUCCESSFUL,"ADOL-C function evaluation with error");
  }
  gradientH.resize(activeVariablesYH.size());
  r = gradient(tAgs[2],activeVariablesYH.size(),&activeVariablesYH[0],&gradientH[0]);
  if(r<adloc_return_value) {
    SETERRQ(PETSC_COMM_SELF,MOFEM_OPERATION_UNSUCCESSFUL,"ADOL-C function evaluation with error");
  }
  hessianH.resize(6+internalVariables.size(),6+internalVariables.size(),false);
  hessianH.clear();
  vector<double*> hessian_h(activeVariablesYH.size());
  for(int dd = 0;dd<activeVariablesYH.size();dd++) {
    hessian_h[dd] = &hessianH(dd,0);
  }
  r = hessian(tAgs[2],activeVariablesYH.size(),&activeVariablesYH[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(
    internalVariables.size(),
    ublas::shallow_array_adaptor<double>(
      internalVariables.size(),
      &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(internalVariables.size(),internalVariables.size(),false);
  for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
    for(unsigned int jj = 0;jj<=ii;jj++) {
      partial2HFlux(ii,jj) = -hessianH(6+ii,6+jj);
    }
  }
  partial2HSigmaFlux.resize(6,internalVariables.size(),false);
  partial2HSigmaFlux.clear();
  for(int ii = 0;ii<6;ii++) {
    for(unsigned int jj = 0;jj<internalVariables.size();jj++) {
      partial2HSigmaFlux(ii,jj) = -hessianH(6+jj,ii);
    }
  }
  PetscFunctionReturn(0);
}

pLayW()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::pLayW ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 762 of file SmallStrainPlasticity.cpp.

                                                                  {
  PetscFunctionBegin;
  int r;
  int adloc_return_value = 0;
  r = ::function(tAgs[0],1,activeVariablesW.size(),&activeVariablesW[0],&w);
  if(r<adloc_return_value) {
    SETERRQ(PETSC_COMM_SELF,MOFEM_OPERATION_UNSUCCESSFUL,"ADOL-C function evaluation with error");
  }
  gradientW.resize(activeVariablesW.size(),false);
  r = gradient(tAgs[0],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(tAgs[0],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");
  }
  sTress = VectorAdaptor(6,ublas::shallow_array_adaptor<double>(6,&gradientW[0]));
  internalFluxes = VectorAdaptor(
    internalVariables.size(),
    ublas::shallow_array_adaptor<double>(
      internalVariables.size(),
      &gradientW[12]
    )
  );
  C.resize(6,6,false);
  for(int ii = 0;ii<6;ii++) {
    for(int jj = 0;jj<=ii;jj++) {
      C(ii,jj) = hessianW(ii,jj);
    }
  }
  partialWStrainPlasticStrain.resize(6,6,false);
  for(int ii = 0;ii<6;ii++) {
    for(int jj = 0;jj<6;jj++) {
      partialWStrainPlasticStrain(ii,jj) = hessianW(6+jj,ii);
    }
  }
  D.resize(internalVariables.size(),internalVariables.size(),false);
  for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
    for(unsigned int jj = 0;jj<=ii;jj++) {
      D(ii,jj) = hessianW(12+ii,12+jj);
    }
  }
  PetscFunctionReturn(0);
}

pLayW_NoHessian()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::pLayW_NoHessian ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 813 of file SmallStrainPlasticity.cpp.

                                                                            {
  PetscFunctionBegin;
  int r;
  int adloc_return_value = 0;
  r = ::function(tAgs[0],1,activeVariablesW.size(),&activeVariablesW[0],&w);
  if(r<adloc_return_value) {
    SETERRQ(PETSC_COMM_SELF,MOFEM_OPERATION_UNSUCCESSFUL,"ADOL-C function evaluation with error");
  }
  gradientW.resize(activeVariablesW.size(),false);
  r = gradient(tAgs[0],activeVariablesW.size(),&activeVariablesW[0],&gradientW[0]);
  if(r<adloc_return_value) {
    SETERRQ(PETSC_COMM_SELF,MOFEM_OPERATION_UNSUCCESSFUL,"ADOL-C function evaluation with error");
  }
  sTress = VectorAdaptor(6,ublas::shallow_array_adaptor<double>(6,&gradientW[0]));
  internalFluxes = VectorAdaptor(
    internalVariables.size(),
    ublas::shallow_array_adaptor<double>(
      internalVariables.size(),
      &gradientW[12]
    )
  );
  PetscFunctionReturn(0);
}

pLayY()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::pLayY ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 837 of file SmallStrainPlasticity.cpp.

                                                                  {
  PetscFunctionBegin;
  int r;
  int adloc_return_value = 0;
  r = ::function(tAgs[1],1,activeVariablesYH.size(),&activeVariablesYH[0],&y);
  if(r<adloc_return_value) {
    SETERRQ(PETSC_COMM_SELF,MOFEM_OPERATION_UNSUCCESSFUL,"ADOL-C function evaluation with error");
  }
  gradientY.resize(activeVariablesYH.size());
  r = gradient(tAgs[1],activeVariablesYH.size(),&activeVariablesYH[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(
    internalVariables.size(),
    ublas::shallow_array_adaptor<double>(
      internalVariables.size(),
      &gradientY[6]
    )
  );
  PetscFunctionReturn(0);
}

pLayY_NoGradient()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::pLayY_NoGradient ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 860 of file SmallStrainPlasticity.cpp.

                                                                             {
  PetscFunctionBegin;
  int r;
  int adloc_return_value = 0;
  r = ::function(tAgs[1],1,activeVariablesYH.size(),&activeVariablesYH[0],&y);
  if(r<adloc_return_value) {
    SETERRQ(PETSC_COMM_SELF,MOFEM_OPERATION_UNSUCCESSFUL,"ADOL-C function evaluation with error");
  }
  PetscFunctionReturn(0);
}

rEcordH()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::rEcordH ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 740 of file SmallStrainPlasticity.cpp.

                                                                    {
  PetscFunctionBegin;
  trace_on(tAgs[2]);
  {
    //active variables
    a_sTress.resize(6,false);
    for(int ii = 0;ii<6;ii++) {
      a_sTress[ii] <<= sTress[ii];
    }
    a_internalFluxes.resize(internalVariables.size(),false);
    for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
      a_internalFluxes[ii] <<= internalFluxes[ii];
    }
    //evaluete functions
    ierr = flowPotential(); CHKERRQ(ierr);
    //return variables
    a_h >>= h;
  }
  trace_off();
  PetscFunctionReturn(0);
}

rEcordW()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::rEcordW ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 692 of file SmallStrainPlasticity.cpp.

                                                                    {
  PetscFunctionBegin;
  trace_on(tAgs[0]);
  {
    //active variables
    a_sTrain.resize(6,false);
    for(int ii = 0;ii<6;ii++) {
      a_sTrain[ii] <<= sTrain[ii];
    }
    a_plasticStrain.resize(6,false);
    for(int ii = 0;ii<6;ii++) {
      a_plasticStrain[ii] <<= plasticStrain[ii];
    }
    a_internalVariables.resize(internalVariables.size(),false);
    for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
      a_internalVariables[ii] <<= internalVariables[ii];
    }
    //evaluete functions
    ierr = freeHemholtzEnergy(); CHKERRQ(ierr);
    //return variables
    a_w >>= w;
  }
  trace_off();
  PetscFunctionReturn(0);
}

rEcordY()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::rEcordY ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 718 of file SmallStrainPlasticity.cpp.

                                                                    {
  PetscFunctionBegin;
  trace_on(tAgs[1]);
  {
    //active variables
    a_sTress.resize(6,false);
    for(int ii = 0;ii<6;ii++) {
      a_sTress[ii] <<= sTress[ii];
    }
    a_internalFluxes.resize(internalVariables.size(),false);
    for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
      a_internalFluxes[ii] <<= internalFluxes[ii];
    }
    //evaluete functions
    ierr = yieldFunction(); CHKERRQ(ierr);
    //return variables
    a_y >>= y;
  }
  trace_off();
  PetscFunctionReturn(0);
}

setActiveVariablesW()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::setActiveVariablesW ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 665 of file SmallStrainPlasticity.cpp.

                                                                                {
  PetscFunctionBegin;
  activeVariablesW.resize(12+internalVariables.size(),false);
  for(int ii = 0;ii<6;ii++) {
    activeVariablesW[ii] = sTrain[ii];
  }
  for(int ii = 0;ii<6;ii++) {
    activeVariablesW[6+ii] = plasticStrain[ii];
  }
  for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
    activeVariablesW[12+ii] = internalVariables[ii];
  }
  PetscFunctionReturn(0);
}

setActiveVariablesYH()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::setActiveVariablesYH ( )

virtual

Examples

SmallStrainPlasticity.hpp.

Definition at line 680 of file SmallStrainPlasticity.cpp.

                                                                                 {
  PetscFunctionBegin;
  activeVariablesYH.resize(6+internalVariables.size(),false);
  for(int ii = 0;ii<6;ii++) {
    activeVariablesYH[ii] = sTress[ii];
  }
  for(unsigned int ii = 0;ii<internalVariables.size();ii++) {
    activeVariablesYH[6+ii] = internalFluxes[ii];
  }
  PetscFunctionReturn(0);
}

snesCreate()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::snesCreate

(

)

Examples

SmallStrainPlasticity.hpp.

Definition at line 1027 of file SmallStrainPlasticity.cpp.

                                                                       {
  PetscFunctionBegin;
  ierr = SNESCreate(PETSC_COMM_SELF,&sNes); CHKERRQ(ierr);
  ierr = SNESSetFunction(sNes,R,SmallStrainPlasticityfRes,(void*)this); CHKERRQ(ierr);
  ierr = SNESSetJacobian(sNes,A,A,SmallStrainPlasticityfJac,(void*)this); CHKERRQ(ierr);
  ierr = SNESGetKSP(sNes,&kSp); CHKERRQ(ierr);
  ierr = KSPGetPC(kSp,&pC); CHKERRQ(ierr);
  ierr = KSPSetType(kSp,KSPPREONLY); CHKERRQ(ierr);
  ierr = PCSetType(pC,PCLU); CHKERRQ(ierr);
  ierr = SNESSetTolerances(sNes,tOl,1e-12,0,20,1000); CHKERRQ(ierr);
  ierr = SNESGetLineSearch(sNes,&lineSearch); CHKERRQ(ierr); CHKERRQ(ierr);
  //ierr = SNESLineSearchSetType(lineSearch,SNESLINESEARCHBASIC); CHKERRQ(ierr);
  ierr = SNESLineSearchSetType(lineSearch,SNESLINESEARCHBT); CHKERRQ(ierr);
 
  PetscFunctionReturn(0);
}

snesDestroy()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::snesDestroy

(

)

Examples

SmallStrainPlasticity.hpp.

Definition at line 1044 of file SmallStrainPlasticity.cpp.

                                                                        {
  PetscFunctionBegin;
  ierr = SNESDestroy(&sNes); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

solveColasetProjection()

PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::solveColasetProjection

(

)

Examples

SmallStrainPlasticity.hpp.

Definition at line 1050 of file SmallStrainPlasticity.cpp.

                                                                                   {
  PetscFunctionBegin;
  {
    double *array;
    ierr = VecGetArray(Chi,&array); CHKERRQ(ierr);
    for(int ii = 0;ii<6;ii++) {
      array[ii] = plasticStrain[ii];
    }
    int nb_internal_variables = internalVariables.size();
    for(unsigned int ii = 0;ii<nb_internal_variables;ii++) {
      array[6+ii] = internalVariables[ii];
    }
    deltaGamma = 0;
    ierr = VecRestoreArray(Chi,&array); CHKERRQ(ierr);
  }
  ierr = SNESSolve(sNes,PETSC_NULL,Chi); CHKERRQ(ierr);
  SNESConvergedReason reason;
  ierr = SNESGetConvergedReason(sNes,&reason); CHKERRQ(ierr);
  if(reason < 0) {
    int its;
    ierr = SNESGetIterationNumber(sNes,&its); CHKERRQ(ierr);
//    ierr = PetscPrintf(
//      PETSC_COMM_SELF,
//      "** WARNING Plasticity Closet Point Projection - number of Newton iterations = %D\n",
//      its
//    ); CHKERRQ(ierr);
    noalias(plasticStrain) = plasticStrain0;
    noalias(internalVariables) = internalVariables0;
    deltaGamma = 0;
  } else {
    double *array;
    ierr = VecGetArray(Chi,&array); CHKERRQ(ierr);
    for(int ii = 0;ii<6;ii++) {
      plasticStrain[ii] = array[ii];
    }
    int nb_internal_variables = internalVariables.size();
    for(unsigned int ii = 0;ii<nb_internal_variables;ii++) {
      internalVariables[ii] = array[6+ii];
    }
    deltaGamma = array[6+nb_internal_variables];
    ierr = VecRestoreArray(Chi,&array); CHKERRQ(ierr);
  }
  PetscFunctionReturn(0);
}

yieldFunction()

virtual PetscErrorCode SmallStrainPlasticity::ClosestPointProjection::yieldFunction ( )

inlinevirtual

Reimplemented in SmallStrainParaboloidalPlasticity, and SmallStrainJ2Plasticity.

Examples

SmallStrainPlasticity.hpp.

Definition at line 581 of file SmallStrainPlasticity.hpp.

                                           {
      PetscFunctionBegin;
      SETERRQ(
        PETSC_COMM_SELF,
        MOFEM_NOT_IMPLEMENTED,
        "Not implemented (overload class and implement this function)"
      );
      PetscFunctionReturn(0);
    }

Friends And Related Function Documentation

SmallStrainPlasticityfJac

PetscErrorCode SmallStrainPlasticityfJac ( SNES  , Vec  , Mat  , Mat  , void *  ctx  )

friend

Examples

SmallStrainPlasticity.hpp.

Definition at line 1207 of file SmallStrainPlasticity.cpp.

                                                                                {
  PetscFunctionBegin;
  SmallStrainPlasticity::ClosestPointProjection *cp;
  cp = (SmallStrainPlasticity::ClosestPointProjection*)ctx;
  
  {
    #if PETSC_VERSION_GE(3,6,0)
    const double *array;
    ierr = VecGetArrayRead(chi,&array); CHKERRQ(ierr);
    #else
    double *array;
    ierr = VecGetArray(chi,&array); CHKERRQ(ierr);
    #endif
    for(int ii = 0;ii<6;ii++) {
      cp->plasticStrain[ii] = array[ii];
    }
    for(unsigned int ii = 0;ii<cp->internalVariables.size();ii++) {
      cp->internalVariables[ii] = array[6+ii];
    }
    cp->deltaGamma = array[6+cp->internalVariables.size()];
    #if PETSC_VERSION_GE(3,6,0)
    ierr = VecRestoreArrayRead(chi,&array); CHKERRQ(ierr);
    #else
    ierr = VecRestoreArray(chi,&array); CHKERRQ(ierr);
    #endif
  }
  ierr = cp->evaluatePotentials(); CHKERRQ(ierr);
  ierr = cp->cAlculateA(); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

SmallStrainPlasticityfRes

PetscErrorCode SmallStrainPlasticityfRes ( SNES  , Vec  , Vec  , void *  ctx  )

friend

Examples

SmallStrainPlasticity.hpp.

Definition at line 1154 of file SmallStrainPlasticity.cpp.

                                                                            {
  PetscFunctionBegin;
  SmallStrainPlasticity::ClosestPointProjection *cp;
  cp = (SmallStrainPlasticity::ClosestPointProjection*)ctx;
  
  //ierr = VecView(chi,PETSC_VIEWER_STDOUT_WORLD); CHKERRQ(ierr);
  {
    #if PETSC_VERSION_GE(3,6,0)
    const double *array;
    ierr = VecGetArrayRead(chi,&array); CHKERRQ(ierr);
    #else
    const double *array;
    ierr = VecGetArray(chi,&array); CHKERRQ(ierr);
    #endif
    for(int ii = 0;ii<6;ii++) {
      cp->plasticStrain[ii] = array[ii];
    }
    for(unsigned int ii = 0;ii<cp->internalVariables.size();ii++) {
      cp->internalVariables[ii] = array[6+ii];
    }
    cp->deltaGamma = array[6+cp->internalVariables.size()];
    #if PETSC_VERSION_GE(3,6,0)
    ierr = VecRestoreArrayRead(chi,&array); CHKERRQ(ierr);
    #else
    ierr = VecRestoreArray(chi,&array); CHKERRQ(ierr);
    #endif
  }
  ierr = cp->evaluatePotentials(); CHKERRQ(ierr);
  ierr = cp->cAlculateR(r); CHKERRQ(ierr);
  // ierr = VecView(r,PETSC_VIEWER_STDOUT_WORLD); CHKERRQ(ierr);
  // cout << "sTress " << cp->sTress << endl;
  // cout << "internalVariables " << cp->internalVariables << endl;
  // cout << "C " << cp->C << endl;
  // cout << "D " << cp->D << endl;
  // cout << "y " << cp->y << endl;
  // cout << "h " << cp->h << endl;
  // cout << "partialYSigma " << cp->partialYSigma << endl;
  // cout << "partialYFlux " << cp->partialYFlux << endl;
  // cout << "partialHSigma " << cp->partialHSigma << endl;
  // cout << "partialHFlux " << cp->partialHFlux << endl;
  // cout << "partial2HSigma " << cp->partial2HSigma << endl;
  // cout << "partial2HFlux " << cp->partial2HFlux << endl;
  // cout << "partial2HSigmaFlux " << cp->partial2HSigmaFlux << endl;
  // cout << "C " << cp->C << endl;
  // cout << "D " << cp->D << endl;
  // cout << "partialWStrainPlasticStrain " << cp->partialWStrainPlasticStrain << endl;
  // cout << "plasticStrain " << cp->plasticStrain << endl;
  // cout << "internalVariables " << cp->internalVariables << endl;
  // cout << "delta plasticStrain " << cp->plasticStrain-cp->plasticStrain0 << endl;
  // cout << "delta internalVariables " << cp->internalVariables-cp->internalVariables << endl;
  PetscFunctionReturn(0);
}

Member Data Documentation

A

Mat SmallStrainPlasticity::ClosestPointProjection::A

Examples

SmallStrainPlasticity.hpp.

Definition at line 540 of file SmallStrainPlasticity.hpp.

a_h

adouble SmallStrainPlasticity::ClosestPointProjection::a_h

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 519 of file SmallStrainPlasticity.hpp.

a_internalFluxes

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

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 518 of file SmallStrainPlasticity.hpp.

a_internalVariables

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

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 517 of file SmallStrainPlasticity.hpp.

a_plasticStrain

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

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 517 of file SmallStrainPlasticity.hpp.

a_sTrain

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

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 517 of file SmallStrainPlasticity.hpp.

a_sTress

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

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 518 of file SmallStrainPlasticity.hpp.

a_w

adouble SmallStrainPlasticity::ClosestPointProjection::a_w

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 519 of file SmallStrainPlasticity.hpp.

a_y

adouble SmallStrainPlasticity::ClosestPointProjection::a_y

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 519 of file SmallStrainPlasticity.hpp.

activeVariablesW

VectorDouble SmallStrainPlasticity::ClosestPointProjection::activeVariablesW

Examples

SmallStrainPlasticity.hpp.

Definition at line 510 of file SmallStrainPlasticity.hpp.

activeVariablesYH

VectorDouble SmallStrainPlasticity::ClosestPointProjection::activeVariablesYH

Examples

SmallStrainPlasticity.hpp.

Definition at line 513 of file SmallStrainPlasticity.hpp.

C

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

Examples

SmallStrainPlasticity.hpp.

Definition at line 496 of file SmallStrainPlasticity.hpp.

Cep

MatrixDouble SmallStrainPlasticity::ClosestPointProjection::Cep

Examples

SmallStrainPlasticity.hpp.

Definition at line 566 of file SmallStrainPlasticity.hpp.

Chi

Vec SmallStrainPlasticity::ClosestPointProjection::Chi

Examples

SmallStrainPlasticity.hpp.

Definition at line 542 of file SmallStrainPlasticity.hpp.

Cp

MatrixDouble SmallStrainPlasticity::ClosestPointProjection::Cp

Examples

SmallStrainPlasticity.hpp.

Definition at line 566 of file SmallStrainPlasticity.hpp.

D

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

Examples

SmallStrainPlasticity.hpp.

Definition at line 496 of file SmallStrainPlasticity.hpp.

dataA

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

Examples

SmallStrainPlasticity.hpp.

Definition at line 541 of file SmallStrainPlasticity.hpp.

dChi

Vec SmallStrainPlasticity::ClosestPointProjection::dChi

Examples

SmallStrainPlasticity.hpp.

Definition at line 542 of file SmallStrainPlasticity.hpp.

deltaGamma

double SmallStrainPlasticity::ClosestPointProjection::deltaGamma

Examples

SmallStrainPlasticity.hpp.

Definition at line 490 of file SmallStrainPlasticity.hpp.

Ep

MatrixDouble SmallStrainPlasticity::ClosestPointProjection::Ep

Examples

SmallStrainPlasticity.hpp.

Definition at line 566 of file SmallStrainPlasticity.hpp.

gG

int SmallStrainPlasticity::ClosestPointProjection::gG

Examples

SmallStrainPlasticity.hpp.

Definition at line 492 of file SmallStrainPlasticity.hpp.

gradientH

VectorDouble SmallStrainPlasticity::ClosestPointProjection::gradientH

Examples

SmallStrainPlasticity.hpp.

Definition at line 536 of file SmallStrainPlasticity.hpp.

gradientW

VectorDouble SmallStrainPlasticity::ClosestPointProjection::gradientW

Examples

SmallStrainPlasticity.hpp.

Definition at line 527 of file SmallStrainPlasticity.hpp.

gradientY

VectorDouble SmallStrainPlasticity::ClosestPointProjection::gradientY

Examples

SmallStrainPlasticity.hpp.

Definition at line 532 of file SmallStrainPlasticity.hpp.

h

double SmallStrainPlasticity::ClosestPointProjection::h

Examples

SmallStrainPlasticity.hpp.

Definition at line 516 of file SmallStrainPlasticity.hpp.

hessianH

MatrixDouble SmallStrainPlasticity::ClosestPointProjection::hessianH

Examples

SmallStrainPlasticity.hpp.

Definition at line 537 of file SmallStrainPlasticity.hpp.

hessianW

MatrixDouble SmallStrainPlasticity::ClosestPointProjection::hessianW

Examples

SmallStrainPlasticity.hpp.

Definition at line 528 of file SmallStrainPlasticity.hpp.

internalFluxes

VectorAdaptor SmallStrainPlasticity::ClosestPointProjection::internalFluxes

Examples

SmallStrainPlasticity.hpp.

Definition at line 495 of file SmallStrainPlasticity.hpp.

internalVariables

VectorDouble SmallStrainPlasticity::ClosestPointProjection::internalVariables

Examples

SmallStrainPlasticity.hpp, and SmallStrainPlasticityMaterialModels.hpp.

Definition at line 486 of file SmallStrainPlasticity.hpp.

internalVariables0

VectorDouble SmallStrainPlasticity::ClosestPointProjection::internalVariables0

Examples

SmallStrainPlasticity.hpp.

Definition at line 487 of file SmallStrainPlasticity.hpp.

kSp

KSP SmallStrainPlasticity::ClosestPointProjection::kSp

Examples

SmallStrainPlasticity.hpp.

Definition at line 557 of file SmallStrainPlasticity.hpp.

lineSearch

SNESLineSearch SmallStrainPlasticity::ClosestPointProjection::lineSearch

Examples

SmallStrainPlasticity.hpp.

Definition at line 559 of file SmallStrainPlasticity.hpp.

Lp

MatrixDouble SmallStrainPlasticity::ClosestPointProjection::Lp

Examples

SmallStrainPlasticity.hpp.

Definition at line 566 of file SmallStrainPlasticity.hpp.

partial2HFlux

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

Examples

SmallStrainPlasticity.hpp.

Definition at line 503 of file SmallStrainPlasticity.hpp.

partial2HSigma

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

Examples

SmallStrainPlasticity.hpp.

Definition at line 502 of file SmallStrainPlasticity.hpp.

partial2HSigmaFlux

MatrixDouble SmallStrainPlasticity::ClosestPointProjection::partial2HSigmaFlux

Examples

SmallStrainPlasticity.hpp.

Definition at line 504 of file SmallStrainPlasticity.hpp.

partialHFlux

VectorAdaptor SmallStrainPlasticity::ClosestPointProjection::partialHFlux

Examples

SmallStrainPlasticity.hpp.

Definition at line 501 of file SmallStrainPlasticity.hpp.

partialHSigma

VectorAdaptor SmallStrainPlasticity::ClosestPointProjection::partialHSigma

Examples

SmallStrainPlasticity.hpp.

Definition at line 500 of file SmallStrainPlasticity.hpp.

partialWStrainPlasticStrain

MatrixDouble SmallStrainPlasticity::ClosestPointProjection::partialWStrainPlasticStrain

Examples

SmallStrainPlasticity.hpp.

Definition at line 497 of file SmallStrainPlasticity.hpp.

partialYFlux

VectorAdaptor SmallStrainPlasticity::ClosestPointProjection::partialYFlux

Examples

SmallStrainPlasticity.hpp.

Definition at line 499 of file SmallStrainPlasticity.hpp.

partialYSigma

VectorAdaptor SmallStrainPlasticity::ClosestPointProjection::partialYSigma

Examples

SmallStrainPlasticity.hpp.

Definition at line 498 of file SmallStrainPlasticity.hpp.

pC

PC SmallStrainPlasticity::ClosestPointProjection::pC

Examples

SmallStrainPlasticity.hpp.

Definition at line 558 of file SmallStrainPlasticity.hpp.

plasticStrain

VectorDouble SmallStrainPlasticity::ClosestPointProjection::plasticStrain

Examples

SmallStrainPlasticity.hpp.

Definition at line 489 of file SmallStrainPlasticity.hpp.

plasticStrain0

VectorDouble SmallStrainPlasticity::ClosestPointProjection::plasticStrain0

Examples

SmallStrainPlasticity.hpp.

Definition at line 488 of file SmallStrainPlasticity.hpp.

R

Vec SmallStrainPlasticity::ClosestPointProjection::R

Examples

SmallStrainPlasticity.hpp.

Definition at line 542 of file SmallStrainPlasticity.hpp.

sNes

SNES SmallStrainPlasticity::ClosestPointProjection::sNes

Examples

SmallStrainPlasticity.hpp.

Definition at line 556 of file SmallStrainPlasticity.hpp.

sTrain

VectorDouble SmallStrainPlasticity::ClosestPointProjection::sTrain

Examples

SmallStrainPlasticity.hpp.

Definition at line 485 of file SmallStrainPlasticity.hpp.

sTress

VectorAdaptor SmallStrainPlasticity::ClosestPointProjection::sTress

Examples

SmallStrainPlasticity.hpp.

Definition at line 494 of file SmallStrainPlasticity.hpp.

tAgs

vector<int> SmallStrainPlasticity::ClosestPointProjection::tAgs

Examples

SmallStrainPlasticity.hpp.

Definition at line 506 of file SmallStrainPlasticity.hpp.

tOl

double SmallStrainPlasticity::ClosestPointProjection::tOl

Examples

SmallStrainPlasticity.hpp.

Definition at line 491 of file SmallStrainPlasticity.hpp.

w

double SmallStrainPlasticity::ClosestPointProjection::w

Examples

SmallStrainPlasticity.hpp.

Definition at line 516 of file SmallStrainPlasticity.hpp.

y

double SmallStrainPlasticity::ClosestPointProjection::y

Examples

SmallStrainPlasticity.hpp.

Definition at line 516 of file SmallStrainPlasticity.hpp.

Yp

VectorDouble SmallStrainPlasticity::ClosestPointProjection::Yp

Examples

SmallStrainPlasticity.hpp.

Definition at line 567 of file SmallStrainPlasticity.hpp.

---

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

• users_modules/small_strain_plasticity/src/SmallStrainPlasticity.hpp
• users_modules/small_strain_plasticity/src/impl/SmallStrainPlasticity.cpp