StructurePlasticDiscretization< T, Diag, OffDiag > Class Template Reference

#include <StructurePlasticDiscretization.h>

Inheritance diagram for StructurePlasticDiscretization< T, Diag, OffDiag >:

Collaboration diagram for StructurePlasticDiscretization< T, Diag, OffDiag >:

Public Types
typedef Array< T >	TArray
typedef Vector< T, 3 >	VectorT3
typedef Array< VectorT3 >	VectorT3Array
typedef Gradient< VectorT3 >	VGradType
typedef Array< Gradient < VectorT3 > >	VGradArray
typedef CRMatrix< Diag, OffDiag, VectorT3 >	CCMatrix
typedef CCMatrix::DiagArray	DiagArray
typedef CCMatrix::PairWiseAssembler	CCAssembler
typedef GradientModel< VectorT3 >	VGradModelType
typedef VGradModelType::GradMatrixType	VGradMatrix

Public Member Functions
	StructurePlasticDiscretization (const MeshList &meshes, const GeomFields &geomFields, Field &varField, const Field &muField, const Field &lambdaField, const Field &alphaField, const Field &varGradientField, const Field &temperatureField, const T &referenceTemperature, const T &residualXXStress, const T &residualYYStress, const T &residualZZStress, const bool &thermo, const bool &residualStress, const Field &devStressField, const Field &VMStressField, Field &plasticStrainField, Field &creepConstant, T &A, const T &B, const T &mm, const T &nn, const T &Sy0, const T &timeStep, const int &creepModel, bool fullLinearization=true)
void	discretize (const Mesh &mesh, MultiFieldMatrix &mfmatrix, MultiField &xField, MultiField &rField)
void	discretizeFaces (const Mesh &mesh, const StorageSite &faces, MultiFieldMatrix &mfmatrix, MultiField &xField, MultiField &rField, const bool isBoundary, const bool isSymmetry)
Public Member Functions inherited from Discretization
	Discretization (const MeshList &meshes)
virtual	~Discretization ()
	DEFINE_TYPENAME ("Discretization")

Private Attributes
const GeomFields &	_geomFields
Field &	_varField
const Field &	_muField
const Field &	_lambdaField
const Field &	_alphaField
const Field &	_varGradientField
const Field &	_temperatureField
const T	_referenceTemperature
const T	_residualXXStress
const T	_residualYYStress
const T	_residualZZStress
const bool	_thermo
const bool	_residualStress
const Field &	_devStressField
const Field &	_VMStressField
Field &	_plasticStrainField
Field &	_creepConstant
T	_A
const T	_B
const T	_m
const T	_n
const T	_Sy0
const T	_timeStep
const int	_creepModel
const bool	_fullLinearization

Additional Inherited Members
Protected Attributes inherited from Discretization
const MeshList &	_meshes

Detailed Description

template<class T, class Diag, class OffDiag>
class StructurePlasticDiscretization< T, Diag, OffDiag >

Definition at line 21 of file StructurePlasticDiscretization.h.

Member Typedef Documentation

template<class T , class Diag , class OffDiag >

typedef CCMatrix::PairWiseAssembler StructurePlasticDiscretization< T, Diag, OffDiag >::CCAssembler

Definition at line 33 of file StructurePlasticDiscretization.h.

template<class T , class Diag , class OffDiag >

typedef CRMatrix<Diag,OffDiag,VectorT3> StructurePlasticDiscretization< T, Diag, OffDiag >::CCMatrix

Definition at line 31 of file StructurePlasticDiscretization.h.

template<class T , class Diag , class OffDiag >

typedef CCMatrix::DiagArray StructurePlasticDiscretization< T, Diag, OffDiag >::DiagArray

Definition at line 32 of file StructurePlasticDiscretization.h.

template<class T , class Diag , class OffDiag >

typedef Array<T> StructurePlasticDiscretization< T, Diag, OffDiag >::TArray

Definition at line 25 of file StructurePlasticDiscretization.h.

template<class T , class Diag , class OffDiag >

typedef Vector<T,3> StructurePlasticDiscretization< T, Diag, OffDiag >::VectorT3

Definition at line 26 of file StructurePlasticDiscretization.h.

template<class T , class Diag , class OffDiag >

typedef Array<VectorT3> StructurePlasticDiscretization< T, Diag, OffDiag >::VectorT3Array

Definition at line 27 of file StructurePlasticDiscretization.h.

template<class T , class Diag , class OffDiag >

typedef Array<Gradient<VectorT3> > StructurePlasticDiscretization< T, Diag, OffDiag >::VGradArray

Definition at line 29 of file StructurePlasticDiscretization.h.

template<class T , class Diag , class OffDiag >

typedef VGradModelType::GradMatrixType StructurePlasticDiscretization< T, Diag, OffDiag >::VGradMatrix

Definition at line 36 of file StructurePlasticDiscretization.h.

template<class T , class Diag , class OffDiag >

typedef GradientModel<VectorT3> StructurePlasticDiscretization< T, Diag, OffDiag >::VGradModelType

Definition at line 35 of file StructurePlasticDiscretization.h.

template<class T , class Diag , class OffDiag >

typedef Gradient<VectorT3> StructurePlasticDiscretization< T, Diag, OffDiag >::VGradType

Definition at line 28 of file StructurePlasticDiscretization.h.

Constructor & Destructor Documentation

template<class T , class Diag , class OffDiag >

StructurePlasticDiscretization< T, Diag, OffDiag >::StructurePlasticDiscretization ( const MeshList &  meshes, const GeomFields &  geomFields, Field &  varField, const Field &  muField, const Field &  lambdaField, const Field &  alphaField, const Field &  varGradientField, const Field &  temperatureField, const T &  referenceTemperature, const T &  residualXXStress, const T &  residualYYStress, const T &  residualZZStress, const bool &  thermo, const bool &  residualStress, const Field &  devStressField, const Field &  VMStressField, Field &  plasticStrainField, Field &  creepConstant, T &  A, const T &  B, const T &  mm, const T &  nn, const T &  Sy0, const T &  timeStep, const int &  creepModel, bool  fullLinearization = true  )

inline

Definition at line 38 of file StructurePlasticDiscretization.h.

                                                              :
    Discretization(meshes),
    _geomFields(geomFields),
    _varField(varField),
    _muField(muField),
    _lambdaField(lambdaField),
    _alphaField(alphaField),
    _varGradientField(varGradientField),
    _temperatureField(temperatureField),
    _referenceTemperature(referenceTemperature),
    _residualXXStress(residualXXStress),
    _residualYYStress(residualYYStress),
    _residualZZStress(residualZZStress),
    _thermo(thermo),
    _residualStress(residualStress),
    _devStressField(devStressField),
    _VMStressField(VMStressField),
    _plasticStrainField(plasticStrainField),
    _creepConstant(creepConstant),
    _A(A),
    _B(B),
    _m(mm),
    _n(nn),
    _Sy0(Sy0),
    _timeStep(timeStep),
    _creepModel(creepModel),
    _fullLinearization(fullLinearization)
   {}

Member Function Documentation

template<class T , class Diag , class OffDiag >

void StructurePlasticDiscretization< T, Diag, OffDiag >::discretize ( const Mesh &  mesh, MultiFieldMatrix &  mfmatrix, MultiField &  xField, MultiField &  rField  )

inlinevirtual

Implements Discretization.

Definition at line 92 of file StructurePlasticDiscretization.h.

References StructurePlasticDiscretization< T, Diag, OffDiag >::_A, StructurePlasticDiscretization< T, Diag, OffDiag >::_B, StructurePlasticDiscretization< T, Diag, OffDiag >::_creepConstant, StructurePlasticDiscretization< T, Diag, OffDiag >::_creepModel, StructurePlasticDiscretization< T, Diag, OffDiag >::_devStressField, StructurePlasticDiscretization< T, Diag, OffDiag >::_m, StructurePlasticDiscretization< T, Diag, OffDiag >::_n, StructurePlasticDiscretization< T, Diag, OffDiag >::_plasticStrainField, StructurePlasticDiscretization< T, Diag, OffDiag >::_Sy0, StructurePlasticDiscretization< T, Diag, OffDiag >::_timeStep, StructurePlasticDiscretization< T, Diag, OffDiag >::_VMStressField, StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces(), Mesh::getAllFaceGroups(), Mesh::getCells(), StorageSite::getCount(), Mesh::getInteriorFaceGroup(), FaceGroup::groupType, FaceGroup::site, and sqrt().

  {
    const StorageSite& iFaces = mesh.getInteriorFaceGroup().site;
    
    const StorageSite& cells = mesh.getCells();
    const VGradArray& devStress =
      dynamic_cast<const VGradArray&>(_devStressField[cells]);
    const TArray& VMStress =
      dynamic_cast<const TArray&>(_VMStressField[cells]);
    VGradArray& plasticStrain =
      dynamic_cast<VGradArray&>(_plasticStrainField[cells]);
    const TArray& creepConstant = 
      dynamic_cast<const TArray&> (_creepConstant[cells]);
    const int nCells = cells.getCount();
    const T zero(0.0);
    const T onethird(1.0/3.0);
    const T half(0.5);
    const T twothirds(2.0/3.0);
    const T one(1.0);
    const T onepointfive(1.5);
    const T two(2.0);
    const T three(3.0);
    const T six(6.0);   

    if (_creepModel==1)
    {
        for(int n=0; n<nCells; n++)
        {
            _A = creepConstant[n];
            T VMPlasticStrain = sqrt(half*
                                     (pow(((plasticStrain[n])[0][0]-(plasticStrain[n])[1][1]),2.0) +
                                      pow(((plasticStrain[n])[1][1]-(plasticStrain[n])[2][2]),2.0) +
                                      pow(((plasticStrain[n])[2][2]-(plasticStrain[n])[0][0]),2.0) +
                                      six*(pow(((plasticStrain[n])[0][1]),2.0) +
                                           pow(((plasticStrain[n])[1][2]),2.0) +
                                           pow(((plasticStrain[n])[2][0]),2.0))));
            //VMPlasticStrain = sqrt(pow(plasticStrain[n][0][0],two));
            //VMPlasticStrain = zero;
            T Sy = _Sy0*(one + _B*pow(VMPlasticStrain,_n));
            T mult = _A*(pow((VMStress[n]/Sy),_m))/VMStress[n];
            if (isnan(mult)) mult = 0.0;
            for(int i=0;i<3;i++)
              for(int j=0;j<3;j++)
                (plasticStrain[n])[i][j] += mult*((devStress[n])[i][j])*_timeStep;
        }
    }
    else if(_creepModel==2)
    {
        for(int n=0; n<nCells; n++)
        {
            _A = creepConstant[n];
            if (n==0)
              cout<<"creep initiated"<<endl;
            T mult = (VMStress[n]/_B)-1.;
            if (mult<=zero)
              mult = zero;
            if (mult!=zero)
              cout<<"yielding has begun"<<endl;

            mult = _A*pow(mult,_n)*onepointfive/VMStress[n];

            for(int i=0;i<3;i++)
              for(int j=0;j<3;j++)
                (plasticStrain[n])[i][j] += mult*((devStress[n])[i][j])*_timeStep;
        }
    }


    discretizeFaces(mesh, iFaces, mfmatrix, xField, rField, false, false);

    /*
    foreach(const FaceGroupPtr fgPtr, mesh.getInterfaceGroups())
    {
        const FaceGroup& fg = *fgPtr;
        const StorageSite& faces = fg.site;
        discretizeFaces(mesh, faces, mfmatrix, xField, rField, false, false);
    }
    */
        
    // boundaries and interfaces
    foreach(const FaceGroupPtr fgPtr, mesh.getAllFaceGroups())
    {
        const FaceGroup& fg = *fgPtr;
        const StorageSite& faces = fg.site;
        if (fg.groupType!="interior")
        {
            discretizeFaces(mesh, faces, mfmatrix, xField, rField,
                            fg.groupType!="interface",
                            fg.groupType=="symmetry");
        }
    }
  }

template<class T , class Diag , class OffDiag >

void StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces ( const Mesh &  mesh, const StorageSite &  faces, MultiFieldMatrix &  mfmatrix, MultiField &  xField, MultiField &  rField, const bool  isBoundary, const bool  isSymmetry  )

inline

Definition at line 187 of file StructurePlasticDiscretization.h.

References StructurePlasticDiscretization< T, Diag, OffDiag >::_alphaField, StructurePlasticDiscretization< T, Diag, OffDiag >::_fullLinearization, StructurePlasticDiscretization< T, Diag, OffDiag >::_geomFields, StructurePlasticDiscretization< T, Diag, OffDiag >::_lambdaField, StructurePlasticDiscretization< T, Diag, OffDiag >::_muField, StructurePlasticDiscretization< T, Diag, OffDiag >::_plasticStrainField, StructurePlasticDiscretization< T, Diag, OffDiag >::_referenceTemperature, StructurePlasticDiscretization< T, Diag, OffDiag >::_residualStress, StructurePlasticDiscretization< T, Diag, OffDiag >::_residualXXStress, StructurePlasticDiscretization< T, Diag, OffDiag >::_residualYYStress, StructurePlasticDiscretization< T, Diag, OffDiag >::_residualZZStress, StructurePlasticDiscretization< T, Diag, OffDiag >::_temperatureField, StructurePlasticDiscretization< T, Diag, OffDiag >::_thermo, StructurePlasticDiscretization< T, Diag, OffDiag >::_varField, StructurePlasticDiscretization< T, Diag, OffDiag >::_varGradientField, GeomFields::area, GeomFields::areaMag, GeomFields::coordinate, dot(), Mesh::getCells(), GradientMatrix< T_Scalar >::getCoeff(), CRMatrix< T_Diag, T_OffDiag, X >::getCoeff(), CRMatrix< T_Diag, T_OffDiag, X >::PairWiseAssembler::getCoeff01(), CRMatrix< T_Diag, T_OffDiag, X >::PairWiseAssembler::getCoeff10(), CRConnectivity::getCol(), GradientMatrix< T_Scalar >::getConnectivity(), StorageSite::getCount(), CRMatrix< T_Diag, T_OffDiag, X >::getDiag(), Mesh::getDimension(), Mesh::getFaceCells(), GradientModel< X >::getGradientMatrix(), MultiFieldMatrix::getMatrix(), CRMatrix< T_Diag, T_OffDiag, X >::getPairWiseAssembler(), CRConnectivity::getRow(), StorageSite::getSelfCount(), harmonicAverage(), CRMatrix< T_Diag, T_OffDiag, X >::hasCoeff(), and GeomFields::volume.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

  {
    const StorageSite& cells = mesh.getCells();

    const MultiField::ArrayIndex cVarIndex(&_varField,&cells);

    const VectorT3Array& faceArea =
      dynamic_cast<const VectorT3Array&>(_geomFields.area[faces]);

    const TArray& faceAreaMag =
      dynamic_cast<const TArray&>(_geomFields.areaMag[faces]);

    const VectorT3Array& cellCentroid =
      dynamic_cast<const VectorT3Array&>(_geomFields.coordinate[cells]);

    const VectorT3Array& faceCentroid =
      dynamic_cast<const VectorT3Array&>(_geomFields.coordinate[faces]);


    const TArray& cellVolume =
      dynamic_cast<const TArray&>(_geomFields.volume[cells]);

    const CRConnectivity& faceCells = mesh.getFaceCells(faces);

    VectorT3Array& rCell = dynamic_cast<VectorT3Array&>(rField[cVarIndex]);
    const VectorT3Array& xCell = dynamic_cast<const VectorT3Array&>(xField[cVarIndex]);

    const VGradArray& vGradCell =
      dynamic_cast<const VGradArray&>(_varGradientField[cells]);

    VGradArray& plasticStrain =
      dynamic_cast<VGradArray&>(_plasticStrainField[cells]);

    const TArray& muCell =
      dynamic_cast<const TArray&>(_muField[cells]);

    const TArray& lambdaCell =
      dynamic_cast<const TArray&>(_lambdaField[cells]);

    const TArray& alphaCell =
      dynamic_cast<const TArray&>(_alphaField[cells]);

    const TArray& temperatureCell =
      dynamic_cast<const TArray&>(_temperatureField[cells]);

    CCMatrix& matrix = dynamic_cast<CCMatrix&>(mfmatrix.getMatrix(cVarIndex,
                                                             cVarIndex));
    CCAssembler& assembler = matrix.getPairWiseAssembler(faceCells);
    DiagArray& diag = matrix.getDiag();


    const int nFaces = faces.getCount();

    const VGradMatrix& vgMatrix = VGradModelType::getGradientMatrix(mesh,_geomFields);
    const CRConnectivity& cellCells = vgMatrix.getConnectivity();
    const Array<int>& ccRow = cellCells.getRow();
    const Array<int>& ccCol = cellCells.getCol();

    const int nInteriorCells = cells.getSelfCount();

    const T two(2.0);
    const T three(3.0);
    
    for(int f=0; f<nFaces; f++)
    {
        const int c0 = faceCells(f,0);
        const int c1 = faceCells(f,1);

        const VectorT3& Af = faceArea[f];
        const VectorT3 en = Af/faceAreaMag[f];
        
        VectorT3 ds=cellCentroid[c1]-cellCentroid[c0];

        T vol0 = cellVolume[c0];
        T vol1 = cellVolume[c1];

        T wt0 = vol0/(vol0+vol1);
        T wt1 = vol1/(vol0+vol1);

        if (isBoundary && !isSymmetry)
        {  
            wt0 = T(1.0);
            wt1 = T(0.);
        }
        
        T faceMu(1.0);
        T faceLambda(1.0);
        T faceAlpha(1.0);
        T faceTemperature(1.0);

        Diag& a00 = diag[c0];
        Diag& a11 = diag[c1];
        OffDiag& a01 = assembler.getCoeff01(f);
        OffDiag& a10 = assembler.getCoeff10(f);

        if (vol0 == 0.)
        {
            faceMu = muCell[c1];
            faceLambda = lambdaCell[c1];
            faceAlpha = alphaCell[c1];
            faceTemperature = temperatureCell[c1];
        }
        else if (vol1 == 0.)
        {
            faceMu = muCell[c0];
            faceLambda = lambdaCell[c0];
            faceAlpha = alphaCell[c0];
            faceTemperature = temperatureCell[c0];
        }
        else
        {
            faceMu = harmonicAverage(muCell[c0],muCell[c1]);
            faceLambda = harmonicAverage(lambdaCell[c0],lambdaCell[c1]);
            faceAlpha = harmonicAverage(alphaCell[c0],alphaCell[c1]);
            faceTemperature = harmonicAverage(temperatureCell[c0],temperatureCell[c1]);
        }

        faceMu = muCell[c0]*wt0 + muCell[c1]*wt1;
        faceLambda = lambdaCell[c0]*wt0 + lambdaCell[c1]*wt1;
        faceAlpha = alphaCell[c0]*wt0 + alphaCell[c1]*wt1;
        faceTemperature = temperatureCell[c0]*wt0 + temperatureCell[c1]*wt1;
        
        const VGradType gradF = (vGradCell[c0]*wt0 + vGradCell[c1]*wt1);
        const VGradType plasticStrainF = (plasticStrain[c0]*wt0 + plasticStrain[c1]*wt1);

        VectorT3 source(NumTypeTraits<VectorT3>::getZero());
        VectorT3 thermalSource(NumTypeTraits<VectorT3>::getZero());
        VectorT3 residualSource(NumTypeTraits<VectorT3>::getZero());
        const T divU = (gradF[0][0] + gradF[1][1] + gradF[2][2]);

        const T diffMetric = faceAreaMag[f]*faceAreaMag[f]/dot(faceArea[f],ds);
        const VectorT3 secondaryCoeff = faceMu*(faceArea[f]-ds*diffMetric);

        // mu*grad U ^ T + lambda * div U I
        source[0] = faceMu*(gradF[0][0]*Af[0] + gradF[0][1]*Af[1] + gradF[0][2]*Af[2])
          + faceLambda*divU*Af[0];
        
        source[1] = faceMu*(gradF[1][0]*Af[0] + gradF[1][1]*Af[1] + gradF[1][2]*Af[2])
          + faceLambda*divU*Af[1];
        
        source[2] = faceMu*(gradF[2][0]*Af[0] + gradF[2][1]*Af[1] + gradF[2][2]*Af[2])
          + faceLambda*divU*Af[2];

        if(_thermo)
        {
            if(mesh.getDimension()==2)
              thermalSource -= (two*faceLambda+two*faceMu)*faceAlpha*(faceTemperature-_referenceTemperature)*Af;
            else
              thermalSource -= (three*faceLambda+two*faceMu)*faceAlpha*(faceTemperature-_referenceTemperature)*Af;
        }

        if(_residualStress)
        {
            residualSource[0] = _residualXXStress*Af[0];
            residualSource[1] = _residualYYStress*Af[1];
            residualSource[2] = _residualZZStress*Af[2];
        }

        VectorT3 creepSource(NumTypeTraits<VectorT3>::getZero());
        T divEP = (plasticStrainF[0][0] + plasticStrainF[1][1] + plasticStrainF[2][2]);
        if (mesh.getDimension() == 2)
          {
            divEP = plasticStrainF[0][0] + plasticStrainF[1][1];
          }

        creepSource[0] -= (two*faceMu*(plasticStrainF[0][0]*Af[0] + plasticStrainF[1][0]*Af[1] + 
                                       plasticStrainF[2][0]*Af[2]) + faceLambda*divEP*Af[0]);
        creepSource[1] -= (two*faceMu*(plasticStrainF[0][1]*Af[0] + plasticStrainF[1][1]*Af[1] + 
                                       plasticStrainF[2][1]*Af[2]) + faceLambda*divEP*Af[1]);
        creepSource[2] -= (two*faceMu*(plasticStrainF[0][2]*Af[0] + plasticStrainF[1][2]*Af[1] + 
                                       plasticStrainF[2][2]*Af[2]) + faceLambda*divEP*Af[2]);
        VectorT3 s0(NumTypeTraits<VectorT3>::getZero());
        VectorT3 s1(NumTypeTraits<VectorT3>::getZero());
        if (_fullLinearization)
        {
            // loop over level 0 neighbors
            for(int nnb = ccRow[c0]; nnb<ccRow[c0+1]; nnb++)
            {
                const int nb = ccCol[nnb];

                // get the coefficient from the gradient matrix that uses modified cellCells
                // getting a copy rather than reference since we might modify it
                VectorT3 g_nb = vgMatrix.getCoeff(c0,nb);
#if 1
                if (isSymmetry)
                {
                    const T gnb_dot_nx2 = T(2.0)*dot(en,g_nb);
                    g_nb = gnb_dot_nx2*en;
                }
#endif
                Diag coeff;

                for(int i=0; i<3; i++)
                {
                    for(int j=0; j<3; j++)
                    {
                        coeff(i,j) = wt0*(faceMu*Af[j]*g_nb[i] 
                                          + faceLambda*Af[i]*g_nb[j]
                                          );
                    }

                    for(int k=0; k<3; k++)
                      coeff(i,i) += wt0*secondaryCoeff[k]*g_nb[k];
                }
#if 0
                if (isSymmetry)
                {
                    SquareTensor<T,3> R;
                    
                    for(int i=0;i<3;i++)
                      for(int j=0; j<3; j++)
                      {
                          if (i==j)
                            R(i,j) = 1.0 - 2*en[i]*en[j];
                          else
                            R(i,j) = - 2*en[i]*en[j];
                      }

                    Diag coeff1(R*coeff*R);
                    coeff += coeff1;
                      
                }
#endif
                OffDiag& a0_nb = matrix.getCoeff(c0,nb);

                a0_nb += coeff;
                a00 -= coeff;
                a10 += coeff;
                
                if (c1 != nb)
                {
                    // if the coefficient does not exist we assume c1
                    // is a ghost cell for which we don't want an
                    // equation in this matrix
                    if (matrix.hasCoeff(c1,nb))
                    {
                        OffDiag& a1_nb = matrix.getCoeff(c1,nb);
                        a1_nb -= coeff;
                    }
                }
                else
                  a11 -= coeff;
            }


            if (!isBoundary)
            {
                for(int nnb = ccRow[c1]; nnb<ccRow[c1+1]; nnb++)
                {
                    const int nb = ccCol[nnb];
                    const VectorT3& g_nb = vgMatrix.getCoeff(c1,nb);
                    
                    Diag coeff;
                    
                    for(int i=0; i<3; i++)
                    {
                        for(int j=0; j<3; j++)
                        {
                            coeff(i,j) = wt1*(faceMu*Af[j]*g_nb[i] 
                                              + faceLambda*Af[i]*g_nb[j]
                                              );
                        }
                        
                        for(int k=0; k<3; k++)
                          coeff(i,i) += wt1*secondaryCoeff[k]*g_nb[k];
                    }
                    
                    
                    if (matrix.hasCoeff(c1,nb))
                    {
                        OffDiag& a1_nb = matrix.getCoeff(c1,nb);
                        a1_nb -= coeff;
                        a11 += coeff;
                    }
                    a01 -= coeff;
                    
                    
                    if (c0 != nb)
                    {
                        OffDiag& a0_nb = matrix.getCoeff(c0,nb);
                        
                        a0_nb += coeff;
                    }
                    else
                      a00 += coeff;
                    
                }
            }
        }
        
        // mu*gradU, primary part
        

        source += faceMu*diffMetric*(xCell[c1]-xCell[c0]);

        // mu*gradU, secondart part


        source += gradF*secondaryCoeff;

        // add flux to the residual of c0 and c1
        rCell[c0] += source;
        rCell[c1] -= source;

        // add flux due to thermal Stress to the residual of c0 and c1
        rCell[c0] += thermalSource;
        rCell[c1] -= thermalSource;

        // add flux due to residual Stress to the residual of c0 and c1
        rCell[c0] += residualSource;
        rCell[c1] -= residualSource;

        // add flux due to creep to the residual of c0 and c1
        rCell[c0] += creepSource;
        rCell[c1] -= creepSource;

        // for Jacobian, use 2*mu + lambda as the diffusivity
        const T faceDiffusivity = faceMu;
        const T diffCoeff = faceDiffusivity*diffMetric;

        
        a01 +=diffCoeff;
        a10 +=diffCoeff;

        
        a00 -= diffCoeff;
        a11 -= diffCoeff;

    }
  }

Member Data Documentation

template<class T , class Diag , class OffDiag >

T StructurePlasticDiscretization< T, Diag, OffDiag >::_A

private

Definition at line 540 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

template<class T , class Diag , class OffDiag >

const Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_alphaField

private

Definition at line 527 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const T StructurePlasticDiscretization< T, Diag, OffDiag >::_B

private

Definition at line 541 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

template<class T , class Diag , class OffDiag >

Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_creepConstant

private

Definition at line 539 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

template<class T , class Diag , class OffDiag >

const int StructurePlasticDiscretization< T, Diag, OffDiag >::_creepModel

private

Definition at line 546 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

template<class T , class Diag , class OffDiag >

const Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_devStressField

private

Definition at line 536 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

template<class T , class Diag , class OffDiag >

const bool StructurePlasticDiscretization< T, Diag, OffDiag >::_fullLinearization

private

Definition at line 547 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const GeomFields& StructurePlasticDiscretization< T, Diag, OffDiag >::_geomFields

private

Definition at line 523 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_lambdaField

private

Definition at line 526 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const T StructurePlasticDiscretization< T, Diag, OffDiag >::_m

private

Definition at line 542 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

template<class T , class Diag , class OffDiag >

const Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_muField

private

Definition at line 525 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const T StructurePlasticDiscretization< T, Diag, OffDiag >::_n

private

Definition at line 543 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

template<class T , class Diag , class OffDiag >

Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_plasticStrainField

private

Definition at line 538 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize(), and StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const T StructurePlasticDiscretization< T, Diag, OffDiag >::_referenceTemperature

private

Definition at line 530 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const bool StructurePlasticDiscretization< T, Diag, OffDiag >::_residualStress

private

Definition at line 535 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const T StructurePlasticDiscretization< T, Diag, OffDiag >::_residualXXStress

private

Definition at line 531 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const T StructurePlasticDiscretization< T, Diag, OffDiag >::_residualYYStress

private

Definition at line 532 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const T StructurePlasticDiscretization< T, Diag, OffDiag >::_residualZZStress

private

Definition at line 533 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const T StructurePlasticDiscretization< T, Diag, OffDiag >::_Sy0

private

Definition at line 544 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

template<class T , class Diag , class OffDiag >

const Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_temperatureField

private

Definition at line 529 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const bool StructurePlasticDiscretization< T, Diag, OffDiag >::_thermo

private

Definition at line 534 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const T StructurePlasticDiscretization< T, Diag, OffDiag >::_timeStep

private

Definition at line 545 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

template<class T , class Diag , class OffDiag >

Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_varField

private

Definition at line 524 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_varGradientField

private

Definition at line 528 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretizeFaces().

template<class T , class Diag , class OffDiag >

const Field& StructurePlasticDiscretization< T, Diag, OffDiag >::_VMStressField

private

Definition at line 537 of file StructurePlasticDiscretization.h.

Referenced by StructurePlasticDiscretization< T, Diag, OffDiag >::discretize().

---

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

• StructurePlasticDiscretization.h