CG  Version 25
Macros
viscoPlasticMacrosCpp.h File Reference
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Macros

#define defineViscoPlasticCoefficients(esr)
 
#define defineViscoPlasticCoefficientsAndTwoDerivatives(esr)
 
#define declareViscoPlasticParameters
 
#define strainRate2d(u0x, u0y, v0x, v0y)   (sqrt( (2./3.)*( SQR(u0x) + SQR(v0y) + .5*SQR( u0y + v0x ) ) )+epsViscoPlastic)
 
#define strainRate3d(u0x, u0y, u0z, v0x, v0y, v0z, w0x, w0y, w0z)   (sqrt( (2./3.)*( SQR(u0x) + SQR(v0y) + SQR(w0z) + .5*( SQR( u0y + v0x ) + SQR( u0z + w0x )+ SQR( v0z + w0y )) ) )+epsViscoPlastic)
 

Macro Definition Documentation

#define declareViscoPlasticParameters
Value:
/* visco-plastic parameters */ \
real nuViscoPlastic=1., etaViscoPlastic=1., yieldStressViscoPlastic=10., exponentViscoPlastic=10.; \
real epsViscoPlastic=1.e-10; /* small parameter used to offset the effective strain rate */ \
parameters.dbase.get<ListOfShowFileParameters>("pdeParameters").getParameter("nuViscoPlastic",nuViscoPlastic); \
parameters.dbase.get<ListOfShowFileParameters>("pdeParameters").getParameter("etaViscoPlastic",etaViscoPlastic); \
parameters.dbase.get<ListOfShowFileParameters>("pdeParameters").getParameter("yieldStressViscoPlastic",yieldStressViscoPlastic); \
parameters.dbase.get<ListOfShowFileParameters>("pdeParameters").getParameter("exponentViscoPlastic",exponentViscoPlastic); \
parameters.dbase.get<ListOfShowFileParameters>("pdeParameters").getParameter("epsViscoPlastic",epsViscoPlastic);

Referenced by InsParameters::getViscoPlasticVariables(), Cgins::saveShowFileComments(), and Cgins::setPlotTitle().

#define defineViscoPlasticCoefficients (   esr)
Value:
exp0 = exp(-exponentViscoPlastic*esr); \
nuT = (etaViscoPlastic + (yieldStressViscoPlastic/esr)*(1.-exp0)); \
nuTd = .5*( (-1./esr)*(1.-exp0) + exponentViscoPlastic*exp0 )*(yieldStressViscoPlastic/(esr*esr)); \
/* -- fake -- */ \
/* nuT = nu+ nuViscoPlastic*esr*esr; */ \
/* nuTd=nuViscoPlastic; */
#define defineViscoPlasticCoefficientsAndTwoDerivatives (   esr)
Value:
exp0 = exp(-exponentViscoPlastic*esr); \
nuT = (etaViscoPlastic + (yieldStressViscoPlastic/esr)*(1.-exp0)); \
nuTd = .5*( (-1./esr)*(1.-exp0) + exponentViscoPlastic*exp0 )*(yieldStressViscoPlastic/(esr*esr)); \
nuTdd= .25*( 3./(esr*esr)*(1.-exp0) -3./(esr)*exponentViscoPlastic*exp0 - exponentViscoPlastic*exponentViscoPlastic*exp0 )*(yieldStressViscoPlastic/(esr*esr*esr))
#define strainRate2d (   u0x,
  u0y,
  v0x,
  v0y 
)    (sqrt( (2./3.)*( SQR(u0x) + SQR(v0y) + .5*SQR( u0y + v0x ) ) )+epsViscoPlastic)
#define strainRate3d (   u0x,
  u0y,
  u0z,
  v0x,
  v0y,
  v0z,
  w0x,
  w0y,
  w0z 
)    (sqrt( (2./3.)*( SQR(u0x) + SQR(v0y) + SQR(w0z) + .5*( SQR( u0y + v0x ) + SQR( u0z + w0x )+ SQR( v0z + w0y )) ) )+epsViscoPlastic)