SmallStrainPlasticity.hpp File Reference

Operators and data structures for small strain plasticity. More...

Go to the source code of this file.

Classes
struct	SmallStrainPlasticity
	Small strain finite element implementation. More...
struct	SmallStrainPlasticity::MyVolumeFE
	definition of volume element More...
struct	SmallStrainPlasticity::CommonData
	common data used by volume elements More...
struct	SmallStrainPlasticity::OpGetDataAtGaussPts
struct	SmallStrainPlasticity::OpGetCommonDataAtGaussPts
struct	SmallStrainPlasticity::OpUpdate
struct	SmallStrainPlasticity::MakeB
struct	SmallStrainPlasticity::OpAssembleRhs
struct	SmallStrainPlasticity::OpAssembleLhs
struct	SmallStrainPlasticity::ClosestPointProjection
	Closest point projection algorithm. More...
struct	SmallStrainPlasticity::OpCalculateStress

Functions
PetscErrorCode	SmallStrainPlasticityfRes (SNES snes, Vec chi, Vec r, void *ctx)
PetscErrorCode	SmallStrainPlasticityfJac (SNES, Vec, Mat, Mat, void *ctx)

Detailed Description

Operators and data structures for small strain plasticity.

Definition in file SmallStrainPlasticity.hpp.

Function Documentation

SmallStrainPlasticityfJac()

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

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	snes,
		Vec	chi,
		Vec	r,
		void *	ctx
	)

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);
}