OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
Public Member Functions | Protected Member Functions | Static Protected Member Functions | Protected Attributes | List of all members
oofem::RheoChainMaterial Class Referenceabstract

This class implements a rheologic chain model describing a viscoelastic material. More...

#include <rheoChM.h>

+ Inheritance diagram for oofem::RheoChainMaterial:
+ Collaboration diagram for oofem::RheoChainMaterial:

Public Member Functions

 RheoChainMaterial (int n, Domain *d)
 
virtual ~RheoChainMaterial ()
 
virtual void giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep)
 Computes the real stress vector for given total strain and integration point. More...
 
virtual void giveRealStressVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. More...
 
virtual void giveRealStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector_3d. More...
 
virtual void giveRealStressVector_PlaneStress (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector_StressControl. More...
 
virtual void giveRealStressVector_1d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector_StressControl. More...
 
virtual void giveRealStressVector_2dBeamLayer (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector_StressControl. More...
 
virtual void giveRealStressVector_PlateLayer (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector_StressControl. More...
 
virtual void giveThermalDilatationVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)
 Returns a vector of coefficients of thermal dilatation in direction of each material principal (local) axis. More...
 
virtual double giveEModulus (GaussPoint *gp, TimeStep *tStep)=0
 Evaluation of the incremental modulus. More...
 
virtual void computeCharCoefficients (FloatArray &answer, double tPrime, GaussPoint *gp, TimeStep *tStep)=0
 Evaluation of the moduli of individual units. More...
 
virtual int hasNonLinearBehaviour ()
 Returns nonzero if receiver is non linear. More...
 
virtual int hasMaterialModeCapability (MaterialMode mode)
 Tests if material supports material mode. More...
 
virtual int hasCastingTimeSupport ()
 Tests if material supports casting time. More...
 
virtual const char * giveClassName () const
 
virtual IRResultType initializeFrom (InputRecord *ir)
 Initializes receiver according to object description stored in input record. More...
 
virtual int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep)
 Returns the integration point corresponding value in Reduced form. More...
 
virtual contextIOResultType saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
 Stores integration point state to output stream. More...
 
virtual contextIOResultType restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
 Reads integration point state to output stream. More...
 
virtual void give3dMaterialStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep)
 Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point. More...
 
virtual void givePlaneStressStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing plane stress stiffness matrix of receiver. More...
 
virtual void givePlaneStrainStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing plane strain stiffness matrix of receiver. More...
 
virtual void give1dStressStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing 1d stiffness matrix of receiver. More...
 
virtual void give2dLatticeStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing 2d lattice stiffness matrix of receiver. More...
 
virtual void give3dLatticeStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing 3d lattice stiffness matrix of receiver. More...
 
virtual void computeStressIndependentStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode)
 Computes reduced strain vector in given integration point, generated by internal processes in material, which are independent on loading in particular integration point. More...
 
virtual void giveShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode)
 Computes, for the given integration point, the strain vector induced by stress-independent shrinkage. More...
 
virtual void giveEigenStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode)
 Computes, for the given integration point, the strain vector induced by the stress history (typically creep strain). More...
 
virtual MaterialStatus * CreateStatus (GaussPoint *gp) const
 Creates new copy of associated status and inserts it into given integration point. More...
 
double giveAlphaOne () const
 
double giveAlphaTwo () const
 
virtual double computeCreepFunction (double t, double t_prime, GaussPoint *gp, TimeStep *tStep)=0
 Evaluation of the creep compliance function at time t when loading is acting from time t_prime. More...
 
virtual bool isActivated (TimeStep *tStep)
 
- Public Member Functions inherited from oofem::StructuralMaterial
 StructuralMaterial (int n, Domain *d)
 Constructor. More...
 
virtual ~StructuralMaterial ()
 Destructor. More...
 
virtual void giveInputRecord (DynamicInputRecord &input)
 Setups the input record string of receiver. More...
 
virtual void giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep)
 Computes the stiffness matrix for giveRealStressVector of receiver in given integration point, respecting its history. More...
 
virtual void giveRealStressVector_StressControl (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, const IntArray &strainControl, TimeStep *tStep)
 Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero· More...
 
virtual void giveRealStressVector_ShellStressControl (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, const IntArray &strainControl, TimeStep *tStep)
 
virtual void giveRealStressVector_Warping (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector_StressControl. More...
 
virtual void giveRealStressVector_Fiber (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector_StressControl. More...
 
virtual void giveRealStressVector_Lattice2d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 
virtual void giveRealStressVector_Lattice3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 
virtual void giveRealStressVector_2dPlateSubSoil (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation is not provided. More...
 
virtual void giveRealStressVector_3dBeamSubSoil (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 
virtual void giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
 Prototype for computation of Eshelby stress. More...
 
void give_dPdF_from (const FloatMatrix &dSdE, FloatMatrix &answer, GaussPoint *gp)
 
void convert_dSdE_2_dPdF (FloatMatrix &answer, const FloatMatrix &dSdE, const FloatArray &S, const FloatArray &F, MaterialMode matMode)
 
double giveReferenceTemperature ()
 Returns the reference temperature of receiver. More...
 
virtual void computeStressIndependentStrainVector_3d (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode)
 
virtual void give3dMaterialStiffnessMatrix_dPdF (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep)
 
virtual void give3dMaterialStiffnessMatrix_dCde (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep)
 
void giveStressDependentPartOfStrainVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrainVector, TimeStep *tStep, ValueModeType mode)
 Method for subtracting from reduced space strain vector its stress-independent parts (caused by temperature, shrinkage, creep and possibly by other phenomena). More...
 
void giveStressDependentPartOfStrainVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrainVector, TimeStep *tStep, ValueModeType mode)
 
virtual int setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type)
 Sets the value of a certain variable at a given integration point to the given value. More...
 
virtual void give2dBeamLayerStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing 2d beam layer stiffness matrix of receiver. More...
 
virtual void givePlateLayerStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing 2d plate layer stiffness matrix of receiver. More...
 
virtual void giveFiberStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing 1d fiber stiffness matrix of receiver. More...
 
virtual void give2dPlateSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing stiffness matrix of plate subsoil model. More...
 
virtual void give3dBeamSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing stiffness matrix of beam3d subsoil model. More...
 
virtual void giveFirstPKStressVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
 Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. More...
 
virtual void giveFirstPKStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
 Default implementation relies on giveFirstPKStressVector_3d. More...
 
virtual void giveFirstPKStressVector_PlaneStress (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
 Default implementation relies on giveFirstPKStressVector_3d. More...
 
virtual void giveFirstPKStressVector_1d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
 Default implementation relies on giveFirstPKStressVector_3d. More...
 
virtual void giveCauchyStressVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
 
virtual void giveCauchyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
 
virtual void giveCauchyStressVector_PlaneStress (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
 
virtual void giveCauchyStressVector_1d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
 
virtual void givePlaneStressStiffMtrx_dPdF (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 
virtual void givePlaneStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 
virtual void givePlaneStrainStiffMtrx_dPdF (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 
virtual void givePlaneStrainStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 
virtual void give1dStressStiffMtrx_dPdF (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 
virtual void give1dStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 
- Public Member Functions inherited from oofem::Material
 Material (int n, Domain *d)
 Constructor. More...
 
virtual ~Material ()
 Destructor. More...
 
virtual bool isCharacteristicMtrxSymmetric (MatResponseMode rMode)
 Returns true if stiffness matrix of receiver is symmetric Default implementation returns true. More...
 
virtual double give (int aProperty, GaussPoint *gp)
 Returns the value of material property 'aProperty'. More...
 
virtual bool hasProperty (int aProperty, GaussPoint *gp)
 Returns true if 'aProperty' exists on material. More...
 
virtual void modifyProperty (int aProperty, double value, GaussPoint *gp)
 Modify 'aProperty', which already exists on material. More...
 
double giveCastingTime ()
 
virtual void printYourself ()
 Prints receiver state on stdout. Useful for debugging. More...
 
virtual int checkConsistency ()
 Allows programmer to test some internal data, before computation begins. More...
 
virtual int initMaterial (Element *element)
 Optional function to call specific procedures when initializing a material. More...
 
virtual MaterialStatus * giveStatus (GaussPoint *gp) const
 Returns material status of receiver in given integration point. More...
 
virtual int packUnknowns (DataStream &buff, TimeStep *tStep, GaussPoint *ip)
 Pack all necessary data of integration point (according to element parallel_mode) into given communication buffer. More...
 
virtual int unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep, GaussPoint *ip)
 Unpack and updates all necessary data of given integration point (according to element parallel_mode) into given communication buffer. More...
 
virtual int estimatePackSize (DataStream &buff, GaussPoint *ip)
 Estimates the necessary pack size to hold all packed data of receiver. More...
 
virtual double predictRelativeComputationalCost (GaussPoint *gp)
 Returns the weight representing relative computational cost of receiver The reference material model is linear isotropic material - its weight is set to 1.0 The other material models should compare to this reference model. More...
 
virtual double predictRelativeRedistributionCost (GaussPoint *gp)
 Returns the relative redistribution cost of the receiver. More...
 
virtual void initTempStatus (GaussPoint *gp)
 Initializes temporary variables stored in integration point status at the beginning of new time step. More...
 
- Public Member Functions inherited from oofem::FEMComponent
 FEMComponent (int n, Domain *d)
 Regular constructor, creates component with given number and belonging to given domain. More...
 
virtual ~FEMComponent ()
 Virtual destructor. More...
 
virtual const char * giveInputRecordName () const =0
 
Domain * giveDomain () const
 
virtual void setDomain (Domain *d)
 Sets associated Domain. More...
 
int giveNumber () const
 
void setNumber (int num)
 Sets number of receiver. More...
 
virtual void updateLocalNumbering (EntityRenumberingFunctor &f)
 Local renumbering support. More...
 
virtual contextIOResultType saveContext (DataStream &stream, ContextMode mode, void *obj=NULL)
 Stores receiver state to output stream. More...
 
virtual contextIOResultType restoreContext (DataStream &stream, ContextMode mode, void *obj=NULL)
 Restores the receiver state previously written in stream. More...
 
virtual void printOutputAt (FILE *file, TimeStep *tStep)
 Prints output of receiver to stream, for given time step. More...
 
virtual Interface * giveInterface (InterfaceType t)
 Interface requesting service. More...
 
std::string errorInfo (const char *func) const
 Returns string for prepending output (used by error reporting macros). More...
 

Protected Member Functions

virtual int hasIncrementalShrinkageFormulation ()
 If only incremental shrinkage strain formulation is provided, then total shrinkage strain must be tracked in status in order to be able to compute total value. More...
 
const FloatArray & giveDiscreteTimes ()
 
void computeDiscreteRelaxationFunction (FloatArray &answer, const FloatArray &tSteps, double t0, double tr, GaussPoint *gp, TimeStep *tSte)
 Evaluation of the relaxation function at given times. More...
 
void giveUnitComplianceMatrix (FloatMatrix &answer, GaussPoint *gp, TimeStep *tStep)
 Evaluation of elastic compliance matrix for unit Young's modulus. More...
 
void giveUnitStiffnessMatrix (FloatMatrix &answer, GaussPoint *gp, TimeStep *tStep)
 Evaluation of elastic stiffness matrix for unit Young's modulus. More...
 
virtual void updateEparModuli (double tPrime, GaussPoint *gp, TimeStep *tStep)
 Update of partial moduli of individual chain units. More...
 
double giveEparModulus (int iChain)
 Access to partial modulus of a given unit. More...
 
virtual void computeCharTimes ()
 Evaluation of characteristic times. More...
 
double giveCharTime (int) const
 Access to the characteristic time of a given unit. More...
 
virtual double giveCharTimeExponent (int i)
 Exponent to be used with the char time of a given unit, usually = 1.0. More...
 
LinearElasticMaterial * giveLinearElasticMaterial ()
 Access to the underlying linear elastic material with unit Young's modulus. More...
 
double giveEndOfTimeOfInterest ()
 Access to the time up to which the response should be accurate. More...
 
void computeTrueStressIndependentStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode)
 Computes, for the given integration point, the strain vector induced by stress-independent internal processes in the material. More...
 

Static Protected Member Functions

static void generateLogTimeScale (FloatArray &answer, double from, double to, int nsteps)
 Generates discrete times starting from time "from" to time "to" uniformly distributed in log time scale. More...
 

Protected Attributes

double talpha
 thermal dilatation coeff. More...
 
int nUnits
 Number of (Maxwell or Kelvin) units in the rheologic chain. More...
 
double relMatAge
 Physical age of the material at simulation time = 0. More...
 
bool lattice
 
double nu
 Poisson's ratio (assumed to be constant, unaffected by creep). More...
 
double alphaOne
 Parameters for the lattice model. More...
 
double alphaTwo
 
double EparValTime
 Time for which the partial moduli of individual units have been evaluated. More...
 
double begOfTimeOfInterest
 Time from which the model should give a good approximation. Optional field. Default value is 0.1 [day]. More...
 
double endOfTimeOfInterest
 Time (age???) up to which the model should give a good approximation. More...
 
LinearElasticMaterial * linearElasticMaterial
 Associated linearElasticMaterial, with E = 1. More...
 
FloatArray EparVal
 Partial moduli of individual units. More...
 
FloatArray charTimes
 Characteristic times of individual units (relaxation or retardation times). More...
 
FloatArray discreteTimeScale
 Times at which the errors are evaluated if the least-square method is used. More...
 
double timeFactor
 Scaling factor transforming the simulation time units into days (gives the number of simulation time units in one day, e.g. More...
 
int preCastingTimeMat
 Stiffness at time less than casting time - optional parameter, negative by default. More...
 
- Protected Attributes inherited from oofem::StructuralMaterial
double referenceTemperature
 Reference temperature (temperature, when material has been built into structure). More...
 
- Protected Attributes inherited from oofem::Material
Dictionary propertyDictionary
 Property dictionary. More...
 
double castingTime
 Casting time. More...
 
- Protected Attributes inherited from oofem::FEMComponent
int number
 Component number. More...
 
Domain * domain
 Link to domain object, useful for communicating with other FEM components. More...
 

Additional Inherited Members

- Static Public Member Functions inherited from oofem::StructuralMaterial
static int giveSymVI (int ind1, int ind2)
 
static int giveVI (int ind1, int ind2)
 
static int giveVoigtVectorMask (IntArray &answer, MaterialMode mmode)
 Returns a mask of the vector indicies corresponding to components in a general (non-symmetric) second order tensor of some stress/strain/deformation measure that performes work. More...
 
static int giveVoigtSymVectorMask (IntArray &answer, MaterialMode mmode)
 The same as giveVoigtVectorMask but returns a mask corresponding to a symmetric second order tensor. More...
 
static void giveInvertedVoigtVectorMask (IntArray &answer, MaterialMode mmode)
 Gives the inverted version of giveVoigtVectorMask. More...
 
static int giveSizeOfVoigtVector (MaterialMode mmode)
 Returns the size of reduced stress/strain vector according to given mode. More...
 
static int giveSizeOfVoigtSymVector (MaterialMode mmode)
 Returns the size of symmetric part of a reduced stress/strain vector according to given mode. More...
 
static void giveFullVectorForm (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
 Converts the reduced symmetric Voigt vector (2nd order tensor) to full form. More...
 
static void giveFullVectorFormF (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
 Converts the reduced deformation gradient Voigt vector (2nd order tensor). More...
 
static void giveFullSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
 Converts the reduced unsymmetric Voigt vector (2nd order tensor) to full form. More...
 
static void giveReducedVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
 Converts the full symmetric Voigt vector (2nd order tensor) to reduced form. More...
 
static void giveReducedSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
 Converts the full unsymmetric Voigt vector (2nd order tensor) to reduced form. More...
 
static void giveFullSymMatrixForm (FloatMatrix &answer, const FloatMatrix &red, MaterialMode matMode)
 Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form. More...
 
static void giveReducedMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
 Converts the full symmetric Voigt matrix (4th order tensor) to reduced form. More...
 
static void giveReducedSymMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
 Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form. More...
 
static void transformStrainVectorTo (FloatArray &answer, const FloatMatrix &base, const FloatArray &strainVector, bool transpose=false)
 Transforms 3d strain vector into another coordinate system. More...
 
static void transformStressVectorTo (FloatArray &answer, const FloatMatrix &base, const FloatArray &stressVector, bool transpose=false)
 Transforms 3d stress vector into another coordinate system. More...
 
static double computeVonMisesStress (const FloatArray *currentStress)
 Computes equivalent of von Mises stress. More...
 
static void giveStrainVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false)
 Computes 3d strain vector transformation matrix from standard vector transformation matrix. More...
 
static void give2DStrainVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false)
 Computes 2d strain vector transformation matrix from standard vector transformation matrix. More...
 
static void giveStressVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false)
 Computes 3d stress vector transformation matrix from standard vector transformation matrix. More...
 
static void givePlaneStressVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false)
 Computes 2d stress vector transformation matrix from standard vector transformation matrix. More...
 
static void sortPrincDirAndValCloseTo (FloatArray *pVal, FloatMatrix *pDir, FloatMatrix *toPDir)
 Method for sorting newly computed principal values (pVal) and corresponding principal directions (pDir) to be closed to some (often previous) principal directions (toPDir). More...
 
static void computePrincipalValues (FloatArray &answer, const FloatArray &s, stressStrainPrincMode mode)
 Common functions for convenience. More...
 
static void computePrincipalValDir (FloatArray &answer, FloatMatrix &dir, const FloatArray &s, stressStrainPrincMode mode)
 Computes principal values and directions of stress or strain vector. More...
 
static double computeDeviatoricVolumetricSplit (FloatArray &dev, const FloatArray &s)
 Computes split of receiver into deviatoric and volumetric part. More...
 
static void computeDeviatoricVolumetricSum (FloatArray &s, const FloatArray &dev, double mean)
 
static void applyDeviatoricElasticCompliance (FloatArray &strain, const FloatArray &stress, double EModulus, double nu)
 
static void applyDeviatoricElasticCompliance (FloatArray &strain, const FloatArray &stress, double GModulus)
 
static void applyDeviatoricElasticStiffness (FloatArray &stress, const FloatArray &strain, double EModulus, double nu)
 
static void applyDeviatoricElasticStiffness (FloatArray &stress, const FloatArray &strain, double GModulus)
 
static void applyElasticStiffness (FloatArray &stress, const FloatArray &strain, double EModulus, double nu)
 
static void applyElasticCompliance (FloatArray &strain, const FloatArray &stress, double EModulus, double nu)
 
static double computeStressNorm (const FloatArray &stress)
 
static double computeFirstInvariant (const FloatArray &s)
 
static double computeSecondStressInvariant (const FloatArray &s)
 
static double computeThirdStressInvariant (const FloatArray &s)
 
static double computeFirstCoordinate (const FloatArray &s)
 
static double computeSecondCoordinate (const FloatArray &s)
 
static double computeThirdCoordinate (const FloatArray &s)
 
- Static Public Attributes inherited from oofem::StructuralMaterial
static std::vector< std::vector< int > > vIindex
 Voigt index map. More...
 
static std::vector< std::vector< int > > svIndex
 Symmetric Voigt index map. More...
 

Detailed Description

This class implements a rheologic chain model describing a viscoelastic material.

It serves as the parent class for Maxwell and Kelvin chains.

Definition at line 130 of file rheoChM.h.

Constructor & Destructor Documentation

oofem::RheoChainMaterial::RheoChainMaterial ( int  n,
Domain *  d 
)

Definition at line 48 of file rheoChM.C.

References EparValTime, linearElasticMaterial, nUnits, preCastingTimeMat, and relMatAge.

oofem::RheoChainMaterial::~RheoChainMaterial ( )
virtual

Definition at line 59 of file rheoChM.C.

References linearElasticMaterial.

Member Function Documentation

virtual void oofem::RheoChainMaterial::computeCharCoefficients ( FloatArray &  answer,
double  tPrime,
GaussPoint *  gp,
TimeStep *  tStep 
)
pure virtual

Evaluation of the moduli of individual units.

Implemented in oofem::MPSMaterial, oofem::Eurocode2CreepMaterial, oofem::B3SolidMaterial, oofem::MaxwellChainMaterial, and oofem::KelvinChainMaterial.

Referenced by updateEparModuli().

void oofem::RheoChainMaterial::computeCharTimes ( )
protectedvirtual

Evaluation of characteristic times.

Reimplemented in oofem::MPSMaterial, oofem::Eurocode2CreepMaterial, and oofem::B3SolidMaterial.

Definition at line 403 of file rheoChM.C.

References oofem::FloatArray::at(), begOfTimeOfInterest, charTimes, giveEndOfTimeOfInterest(), nUnits, OOFEM_ERROR, relMatAge, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by initializeFrom().

virtual double oofem::RheoChainMaterial::computeCreepFunction ( double  t,
double  t_prime,
GaussPoint *  gp,
TimeStep *  tStep 
)
pure virtual

Evaluation of the creep compliance function at time t when loading is acting from time t_prime.

Implemented in oofem::MPSMaterial, oofem::Eurocode2CreepMaterial, oofem::B3SolidMaterial, oofem::B3Material, oofem::KelvinChainSolidMaterial, oofem::CebFip78Material, and oofem::DoublePowerLawMaterial.

Referenced by oofem::KelvinChainMaterial::computeCharCoefficients(), and computeDiscreteRelaxationFunction().

void oofem::RheoChainMaterial::computeDiscreteRelaxationFunction ( FloatArray &  answer,
const FloatArray &  tSteps,
double  t0,
double  tr,
GaussPoint *  gp,
TimeStep *  tSte 
)
protected

Evaluation of the relaxation function at given times.

This functions solves numerically an integral equation of the form

\[ \varepsilon(t) = \int_{0}^{t} J(t, \tau) \mathrm{d}\sigma(\tau) + \varepsilon_n(t) \]

where $ \varepsilon_n(t) $ is stress-independent deformation, for the case where $ \varepsilon(t) = 1 $ is kept at constant value in time.

Parameters
[out]answerArray with evaluated relaxation function.
t0Age of material when load is applied.
trAge of material when relaxation has begun ???
tStepsAt which times the relaxation function will be evaluated.
Warning
tSteps should be uniformly distributed in log time scale and relatively dense (100 intervals) in order to achieve a reasonable accuracy.

Definition at line 193 of file rheoChM.C.

References oofem::FloatArray::at(), computeCreepFunction(), oofem::FloatArray::giveSize(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by oofem::MaxwellChainMaterial::computeCharCoefficients().

void oofem::RheoChainMaterial::computeStressIndependentStrainVector ( FloatArray &  answer,
GaussPoint *  gp,
TimeStep *  tStep,
ValueModeType  mode 
)
virtual

Computes reduced strain vector in given integration point, generated by internal processes in material, which are independent on loading in particular integration point.

Default implementation takes into account temperature induced strains and eigenstrains.

Parameters
answerReturned strain vector.
gpIntegration point.
tStepTime step (most models are able to respond only when tStep is current time step).
modeDetermines response mode (Total or incremental).
Todo:
Hack for loose gausspoints. We shouldn't ask for "gp->giveElement()". FIXME

Reimplemented from oofem::StructuralMaterial.

Definition at line 371 of file rheoChM.C.

References oofem::FloatArray::add(), oofem::Material::castingTime, computeTrueStressIndependentStrainVector(), giveEigenStrainVector(), oofem::TimeStep::giveIntrinsicTime(), and oofem::FloatArray::giveSize().

Referenced by oofem::KelvinChainSolidMaterial::computeHiddenVars(), and oofem::KelvinChainMaterial::computeHiddenVars().

void oofem::RheoChainMaterial::computeTrueStressIndependentStrainVector ( FloatArray &  answer,
GaussPoint *  gp,
TimeStep *  tStep,
ValueModeType  mode 
)
protected

Computes, for the given integration point, the strain vector induced by stress-independent internal processes in the material.

Takes into account only temperature and shrinkage-induced strains.

Parameters
answerReturned strain vector.
gpIntegration point.
tStepTime step (most models are able to respond only when tStep is current time step).
modeDetermines response mode (Total or incremental).

Definition at line 344 of file rheoChM.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::StructuralMaterial::computeStressIndependentStrainVector(), giveShrinkageStrainVector(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), oofem::RheoChainMaterialStatus::giveThermalStrain(), and oofem::RheoChainMaterialStatus::setTempThermalStrain().

Referenced by oofem::MaxwellChainMaterial::computeHiddenVars(), and computeStressIndependentStrainVector().

MaterialStatus * oofem::RheoChainMaterial::CreateStatus ( GaussPoint *  gp) const
virtual

Creates new copy of associated status and inserts it into given integration point.

Parameters
gpIntegration point where newly created status will be stored.
Returns
Reference to new status.

Reimplemented from oofem::Material.

Reimplemented in oofem::MPSDamMaterial, oofem::MPSMaterial, oofem::Eurocode2CreepMaterial, oofem::B3SolidMaterial, oofem::MaxwellChainMaterial, oofem::KelvinChainSolidMaterial, and oofem::KelvinChainMaterial.

Definition at line 576 of file rheoChM.C.

References oofem::FEMComponent::giveDomain(), and nUnits.

void oofem::RheoChainMaterial::generateLogTimeScale ( FloatArray &  answer,
double  from,
double  to,
int  nsteps 
)
staticprotected

Generates discrete times starting from time "from" to time "to" uniformly distributed in log time scale.

The time interval (to-from) is divided to nsteps intervals. We return times starting from ("from" + first increment)

Parameters
answerResulting array of discrete times.
fromStarting time
toEnd time
nstepsNumber of discrete steps.

Definition at line 247 of file rheoChM.C.

References oofem::FloatArray::at(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by initializeFrom().

void oofem::RheoChainMaterial::give1dStressStiffMtrx ( FloatMatrix &  answer,
MatResponseMode  mmode,
GaussPoint *  gp,
TimeStep *  tStep 
)
virtual

Method for computing 1d stiffness matrix of receiver.

Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to 1d stiffness using reduce method described above. However, this reduction is quite time consuming and if it is possible, it is recommended to overload this method and provide direct method for computing particular stiffness matrix.

Parameters
answerStiffness matrix.
mmodeMaterial response mode.
gpIntegration point, which load history is used.
tStepTime step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::MPSDamMaterial.

Definition at line 521 of file rheoChM.C.

References oofem::Material::castingTime, oofem::FEMComponent::domain, oofem::StructuralMaterial::give1dStressStiffMtrx(), giveEModulus(), oofem::TimeStep::giveIntrinsicTime(), giveLinearElasticMaterial(), oofem::Domain::giveMaterial(), oofem::Material::giveStatus(), preCastingTimeMat, and oofem::FloatMatrix::times().

Referenced by oofem::MPSDamMaterial::give1dStressStiffMtrx().

void oofem::RheoChainMaterial::give2dLatticeStiffMtrx ( FloatMatrix &  answer,
MatResponseMode  mmode,
GaussPoint *  gp,
TimeStep *  tStep 
)
virtual

Method for computing 2d lattice stiffness matrix of receiver.

Parameters
answerStiffness matrix.
mmodeMaterial response mode.
gpIntegration point, which load history is used.
tStepTime step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Definition at line 538 of file rheoChM.C.

References oofem::Material::castingTime, oofem::FEMComponent::domain, oofem::StructuralMaterial::give2dLatticeStiffMtrx(), giveEModulus(), oofem::TimeStep::giveIntrinsicTime(), giveLinearElasticMaterial(), oofem::Domain::giveMaterial(), oofem::Material::giveStatus(), preCastingTimeMat, and oofem::FloatMatrix::times().

void oofem::RheoChainMaterial::give3dLatticeStiffMtrx ( FloatMatrix &  answer,
MatResponseMode  mmode,
GaussPoint *  gp,
TimeStep *  tStep 
)
virtual

Method for computing 3d lattice stiffness matrix of receiver.

Parameters
answerStiffness matrix.
mmodeMaterial response mode.
gpIntegration point, which load history is used.
tStepTime step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Definition at line 557 of file rheoChM.C.

References oofem::Material::castingTime, oofem::FEMComponent::domain, oofem::StructuralMaterial::give3dLatticeStiffMtrx(), giveEModulus(), oofem::TimeStep::giveIntrinsicTime(), giveLinearElasticMaterial(), oofem::Domain::giveMaterial(), oofem::Material::giveStatus(), preCastingTimeMat, and oofem::FloatMatrix::times().

void oofem::RheoChainMaterial::give3dMaterialStiffnessMatrix ( FloatMatrix &  answer,
MatResponseMode  mode,
GaussPoint *  gp,
TimeStep *  tStep 
)
virtual

Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point.

Parameters
answerComputed results.
modeMaterial response mode.
gpIntegration point.
tStepTime step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::MPSDamMaterial.

Definition at line 459 of file rheoChM.C.

References oofem::Material::castingTime, oofem::FEMComponent::domain, oofem::StructuralMaterial::give3dMaterialStiffnessMatrix(), giveEModulus(), oofem::TimeStep::giveIntrinsicTime(), giveLinearElasticMaterial(), oofem::Domain::giveMaterial(), oofem::Material::giveStatus(), preCastingTimeMat, and oofem::FloatMatrix::times().

Referenced by oofem::MPSDamMaterial::give3dMaterialStiffnessMatrix().

double oofem::RheoChainMaterial::giveAlphaOne ( ) const
inline

Definition at line 286 of file rheoChM.h.

double oofem::RheoChainMaterial::giveAlphaTwo ( ) const
inline

Definition at line 287 of file rheoChM.h.

double oofem::RheoChainMaterial::giveCharTime ( int  i) const
protected

Access to the characteristic time of a given unit.

Definition at line 397 of file rheoChM.C.

References oofem::FloatArray::at(), and charTimes.

Referenced by oofem::KelvinChainSolidMaterial::computeBetaMu(), oofem::KelvinChainMaterial::computeCharCoefficients(), oofem::MaxwellChainMaterial::computeCharCoefficients(), oofem::B3SolidMaterial::computeCharCoefficients(), oofem::MaxwellChainMaterial::computeHiddenVars(), oofem::KelvinChainMaterial::computeHiddenVars(), oofem::KelvinChainSolidMaterial::computeLambdaMu(), oofem::KelvinChainMaterial::giveEigenStrainVector(), oofem::MaxwellChainMaterial::giveEigenStrainVector(), oofem::KelvinChainMaterial::giveEModulus(), and oofem::MaxwellChainMaterial::giveEModulus().

virtual double oofem::RheoChainMaterial::giveCharTimeExponent ( int  i)
inlineprotectedvirtual

Exponent to be used with the char time of a given unit, usually = 1.0.

Definition at line 356 of file rheoChM.h.

Referenced by oofem::MaxwellChainMaterial::computeCharCoefficients(), oofem::MaxwellChainMaterial::computeHiddenVars(), oofem::MaxwellChainMaterial::giveEigenStrainVector(), and oofem::MaxwellChainMaterial::giveEModulus().

virtual const char* oofem::RheoChainMaterial::giveClassName ( ) const
inlinevirtual
Returns
Class name of the receiver.

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::MPSMaterial, oofem::MPSDamMaterial, oofem::Eurocode2CreepMaterial, oofem::B3SolidMaterial, oofem::B3Material, oofem::MaxwellChainMaterial, oofem::KelvinChainSolidMaterial, oofem::KelvinChainMaterial, oofem::CebFip78Material, and oofem::DoublePowerLawMaterial.

Definition at line 219 of file rheoChM.h.

References oofem::MaterialStatus::initializeFrom().

const FloatArray & oofem::RheoChainMaterial::giveDiscreteTimes ( )
protected

Definition at line 260 of file rheoChM.C.

References discreteTimeScale.

Referenced by oofem::KelvinChainMaterial::computeCharCoefficients(), oofem::MaxwellChainMaterial::computeCharCoefficients(), and oofem::KelvinChainMaterial::initializeFrom().

virtual void oofem::RheoChainMaterial::giveEigenStrainVector ( FloatArray &  answer,
GaussPoint *  gp,
TimeStep *  tStep,
ValueModeType  mode 
)
inlinevirtual

Computes, for the given integration point, the strain vector induced by the stress history (typically creep strain).

Parameters
answerComputed strains.
gpIntegration point.
tStepTime step (most models are able to respond only when tStep is the current time step).
modeDetermines response mode.

Reimplemented in oofem::MPSMaterial, oofem::B3SolidMaterial, oofem::MaxwellChainMaterial, oofem::KelvinChainSolidMaterial, and oofem::KelvinChainMaterial.

Definition at line 282 of file rheoChM.h.

Referenced by computeStressIndependentStrainVector().

virtual double oofem::RheoChainMaterial::giveEModulus ( GaussPoint *  gp,
TimeStep *  tStep 
)
pure virtual

Evaluation of the incremental modulus.

Implemented in oofem::MPSMaterial, oofem::Eurocode2CreepMaterial, oofem::B3SolidMaterial, oofem::MaxwellChainMaterial, oofem::KelvinChainMaterial, and oofem::KelvinChainSolidMaterial.

Referenced by give1dStressStiffMtrx(), give2dLatticeStiffMtrx(), give3dLatticeStiffMtrx(), give3dMaterialStiffnessMatrix(), givePlaneStrainStiffMtrx(), givePlaneStressStiffMtrx(), and giveRealStressVector().

double oofem::RheoChainMaterial::giveEndOfTimeOfInterest ( )
protected

Access to the time up to which the response should be accurate.

Definition at line 647 of file rheoChM.C.

References endOfTimeOfInterest, oofem::FEMComponent::giveDomain(), oofem::EngngModel::giveEndOfTimeOfInterest(), oofem::Domain::giveEngngModel(), and timeFactor.

Referenced by oofem::B3SolidMaterial::computeCharTimes(), oofem::MPSMaterial::computeCharTimes(), computeCharTimes(), and initializeFrom().

double oofem::RheoChainMaterial::giveEparModulus ( int  iChain)
protected

Access to partial modulus of a given unit.

Definition at line 303 of file rheoChM.C.

References oofem::FloatArray::at(), and EparVal.

Referenced by oofem::MaxwellChainMaterial::computeHiddenVars(), oofem::KelvinChainMaterial::computeHiddenVars(), oofem::KelvinChainSolidMaterial::giveEigenStrainVector(), oofem::KelvinChainSolidMaterial::giveEModulus(), oofem::KelvinChainMaterial::giveEModulus(), and oofem::MaxwellChainMaterial::giveEModulus().

int oofem::RheoChainMaterial::giveIPValue ( FloatArray &  answer,
GaussPoint *  gp,
InternalStateType  type,
TimeStep *  tStep 
)
virtual

Returns the integration point corresponding value in Reduced form.

Parameters
answerContain corresponding ip value, zero sized if not available.
gpIntegration point to which the value refers.
typeDetermines the type of internal variable.
tStepDetermines the time step.
Returns
Nonzero if the assignment can be done, zero if this type of variable is not supported.

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::MPSMaterial, and oofem::MPSDamMaterial.

Definition at line 701 of file rheoChM.C.

References oofem::FloatArray::at(), oofem::StructuralMaterial::giveIPValue(), oofem::Material::giveStatus(), oofem::RheoChainMaterialStatus::giveThermalStrain(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by oofem::MPSMaterial::giveIPValue().

LinearElasticMaterial * oofem::RheoChainMaterial::giveLinearElasticMaterial ( )
protected

Access to the underlying linear elastic material with unit Young's modulus.

Definition at line 629 of file rheoChM.C.

References oofem::FEMComponent::giveDomain(), oofem::FEMComponent::giveNumber(), lattice, linearElasticMaterial, and nu.

Referenced by give1dStressStiffMtrx(), give2dLatticeStiffMtrx(), give3dLatticeStiffMtrx(), give3dMaterialStiffnessMatrix(), givePlaneStrainStiffMtrx(), givePlaneStressStiffMtrx(), giveUnitComplianceMatrix(), giveUnitStiffnessMatrix(), restoreIPContext(), and saveIPContext().

void oofem::RheoChainMaterial::givePlaneStrainStiffMtrx ( FloatMatrix &  answer,
MatResponseMode  mmode,
GaussPoint *  gp,
TimeStep *  tStep 
)
virtual

Method for computing plane strain stiffness matrix of receiver.

Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to plane strain stiffness using reduce method described above. However, this reduction is quite time consuming and if it is possible, it is recommended to overload this method and provide direct method for computing particular stiffness matrix. Note: as already described, if zero strain component is imposed (Plane strain, ..) this condition must be taken into account in geometrical relations, and corresponding component has to be included in reduced vector. (So plane strain conditions are $ \epsilon_z = \gamma_{xz} = \gamma_{yz} = 0 $, but relations for $ \epsilon_z$ and $\sigma_z$ are included).

Parameters
answerStiffness matrix.
mmodeMaterial response mode.
gpIntegration point, which load history is used.
tStepTime step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::MPSDamMaterial.

Definition at line 503 of file rheoChM.C.

References oofem::Material::castingTime, oofem::FEMComponent::domain, giveEModulus(), oofem::TimeStep::giveIntrinsicTime(), giveLinearElasticMaterial(), oofem::Domain::giveMaterial(), oofem::StructuralMaterial::givePlaneStrainStiffMtrx(), oofem::Material::giveStatus(), preCastingTimeMat, and oofem::FloatMatrix::times().

Referenced by oofem::MPSDamMaterial::givePlaneStrainStiffMtrx().

void oofem::RheoChainMaterial::givePlaneStressStiffMtrx ( FloatMatrix &  answer,
MatResponseMode  mmode,
GaussPoint *  gp,
TimeStep *  tStep 
)
virtual

Method for computing plane stress stiffness matrix of receiver.

Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to plane stress stiffness using reduce method described above. However, this reduction is quite time consuming and if it is possible, it is recommended to overload this method and provide direct method for computing particular stiffness matrix.

Parameters
answerStiffness matrix.
mmodeMaterial response mode.
gpIntegration point, which load history is used.
tStepTime step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::MPSDamMaterial.

Definition at line 486 of file rheoChM.C.

References oofem::Material::castingTime, oofem::FEMComponent::domain, giveEModulus(), oofem::TimeStep::giveIntrinsicTime(), giveLinearElasticMaterial(), oofem::Domain::giveMaterial(), oofem::StructuralMaterial::givePlaneStressStiffMtrx(), oofem::Material::giveStatus(), preCastingTimeMat, and oofem::FloatMatrix::times().

Referenced by oofem::MPSDamMaterial::givePlaneStressStiffMtrx().

void oofem::RheoChainMaterial::giveRealStressVector ( FloatArray &  answer,
GaussPoint *  gp,
const FloatArray &  reducedStrain,
TimeStep *  tStep 
)
virtual

Computes the real stress vector for given total strain and integration point.

The total strain is defined as strain computed directly from displacement field at given time. The stress independent parts (temperature, eigenstrains) are subtracted in constitutive driver. The service should use previously reached equilibrium history variables. Also it should update temporary history variables in status according to newly reached state. The temporary history variables are moved into equilibrium ones after global structure equilibrium has been reached by iteration process.

Parameters
answerStress vector in reduced form. For large deformations it is treated as the second Piola-Kirchoff stress.
gpIntegration point.
reducedStrainStrain vector in reduced form. For large deformations it is treated as the Green-Lagrange strain.
tStepCurrent time step (most models are able to respond only when tStep is current time step).
Todo:
Move this to StructuralCrossSection ?

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::MPSMaterial, oofem::MPSDamMaterial, oofem::Eurocode2CreepMaterial, oofem::B3SolidMaterial, oofem::KelvinChainMaterial, oofem::MaxwellChainMaterial, and oofem::KelvinChainSolidMaterial.

Definition at line 81 of file rheoChM.C.

References oofem::FloatArray::add(), oofem::FloatArray::beDifferenceOf(), oofem::FloatArray::beProductOf(), oofem::Material::castingTime, oofem::FEMComponent::domain, giveEModulus(), oofem::TimeStep::giveIntrinsicTime(), oofem::Domain::giveMaterial(), oofem::GaussPoint::giveMaterialMode(), giveShrinkageStrainVector(), oofem::FloatArray::giveSize(), oofem::StructuralMaterial::giveSizeOfVoigtSymVector(), oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStiffnessMatrix(), oofem::StructuralMaterialStatus::giveStrainVector(), oofem::StructuralMaterial::giveStressDependentPartOfStrainVector(), giveUnitStiffnessMatrix(), oofem::RheoChainMaterialStatus::giveViscoelasticStressVector(), hasIncrementalShrinkageFormulation(), oofem::Material::initTempStatus(), isActivated(), oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), preCastingTimeMat, oofem::FloatArray::resize(), oofem::RheoChainMaterialStatus::setShrinkageStrainVector(), oofem::FloatMatrix::times(), and oofem::FloatArray::zero().

Referenced by oofem::KelvinChainSolidMaterial::giveRealStressVector(), oofem::MaxwellChainMaterial::giveRealStressVector(), and oofem::KelvinChainMaterial::giveRealStressVector().

virtual void oofem::RheoChainMaterial::giveRealStressVector_1d ( FloatArray &  answer,
GaussPoint *  gp,
const FloatArray &  reducedE,
TimeStep *  tStep 
)
inlinevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 187 of file rheoChM.h.

virtual void oofem::RheoChainMaterial::giveRealStressVector_2dBeamLayer ( FloatArray &  answer,
GaussPoint *  gp,
const FloatArray &  reducedE,
TimeStep *  tStep 
)
inlinevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 189 of file rheoChM.h.

virtual void oofem::RheoChainMaterial::giveRealStressVector_3d ( FloatArray &  answer,
GaussPoint *  gp,
const FloatArray &  reducedE,
TimeStep *  tStep 
)
inlinevirtual

Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.

Reimplemented from oofem::StructuralMaterial.

Definition at line 181 of file rheoChM.h.

virtual void oofem::RheoChainMaterial::giveRealStressVector_PlaneStrain ( FloatArray &  answer,
GaussPoint *  gp,
const FloatArray &  reducedE,
TimeStep *  tStep 
)
inlinevirtual

Default implementation relies on giveRealStressVector_3d.

Reimplemented from oofem::StructuralMaterial.

Definition at line 183 of file rheoChM.h.

virtual void oofem::RheoChainMaterial::giveRealStressVector_PlaneStress ( FloatArray &  answer,
GaussPoint *  gp,
const FloatArray &  reducedE,
TimeStep *  tStep 
)
inlinevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 185 of file rheoChM.h.

virtual void oofem::RheoChainMaterial::giveRealStressVector_PlateLayer ( FloatArray &  answer,
GaussPoint *  gp,
const FloatArray &  reducedE,
TimeStep *  tStep 
)
inlinevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 191 of file rheoChM.h.

References oofem::IntegrationPointStatus::gp.

virtual void oofem::RheoChainMaterial::giveShrinkageStrainVector ( FloatArray &  answer,
GaussPoint *  gp,
TimeStep *  tStep,
ValueModeType  mode 
)
inlinevirtual

Computes, for the given integration point, the strain vector induced by stress-independent shrinkage.

Parameters
answerReturned strain vector.
gpIntegration point.
tStepTime step (most models are able to respond only when tStep is current time step).
modeDetermines response mode (Total or incremental).

Reimplemented in oofem::MPSMaterial, oofem::Eurocode2CreepMaterial, oofem::B3SolidMaterial, oofem::B3Material, oofem::MaxwellChainMaterial, oofem::KelvinChainSolidMaterial, and oofem::KelvinChainMaterial.

Definition at line 267 of file rheoChM.h.

References oofem::FloatArray::clear().

Referenced by computeTrueStressIndependentStrainVector(), and giveRealStressVector().

void oofem::RheoChainMaterial::giveThermalDilatationVector ( FloatArray &  answer,
GaussPoint *  gp,
TimeStep *  tStep 
)
virtual

Returns a vector of coefficients of thermal dilatation in direction of each material principal (local) axis.

Parameters
answerVector of thermal dilatation coefficients.
gpIntegration point.
tStepTime step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Definition at line 169 of file rheoChM.C.

References oofem::FloatArray::at(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatArray::resize(), talpha, and oofem::FloatArray::zero().

void oofem::RheoChainMaterial::giveUnitComplianceMatrix ( FloatMatrix &  answer,
GaussPoint *  gp,
TimeStep *  tStep 
)
protected

Evaluation of elastic compliance matrix for unit Young's modulus.

Definition at line 286 of file rheoChM.C.

References oofem::FloatMatrix::beInverseOf(), giveLinearElasticMaterial(), and oofem::StructuralMaterial::giveStiffnessMatrix().

Referenced by oofem::KelvinChainSolidMaterial::giveEigenStrainVector(), oofem::MaxwellChainMaterial::giveEigenStrainVector(), and oofem::B3SolidMaterial::giveEigenStrainVector().

void oofem::RheoChainMaterial::giveUnitStiffnessMatrix ( FloatMatrix &  answer,
GaussPoint *  gp,
TimeStep *  tStep 
)
protected

Evaluation of elastic stiffness matrix for unit Young's modulus.

Definition at line 267 of file rheoChM.C.

References giveLinearElasticMaterial(), and oofem::StructuralMaterial::giveStiffnessMatrix().

Referenced by oofem::KelvinChainSolidMaterial::computeHiddenVars(), oofem::MaxwellChainMaterial::computeHiddenVars(), and giveRealStressVector().

virtual int oofem::RheoChainMaterial::hasCastingTimeSupport ( )
inlinevirtual

Tests if material supports casting time.

Returns
Nonzero if supported, zero otherwise.

Reimplemented from oofem::Material.

Definition at line 217 of file rheoChM.h.

virtual int oofem::RheoChainMaterial::hasIncrementalShrinkageFormulation ( )
inlineprotectedvirtual

If only incremental shrinkage strain formulation is provided, then total shrinkage strain must be tracked in status in order to be able to compute total value.

Reimplemented in oofem::MPSMaterial, oofem::Eurocode2CreepMaterial, oofem::B3SolidMaterial, oofem::B3Material, oofem::KelvinChainSolidMaterial, oofem::MaxwellChainMaterial, and oofem::KelvinChainMaterial.

Definition at line 306 of file rheoChM.h.

Referenced by giveRealStressVector().

int oofem::RheoChainMaterial::hasMaterialModeCapability ( MaterialMode  mode)
virtual

Tests if material supports material mode.

Parameters
modeRequired material mode.
Returns
Nonzero if supported, zero otherwise.

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::MPSDamMaterial.

Definition at line 68 of file rheoChM.C.

virtual int oofem::RheoChainMaterial::hasNonLinearBehaviour ( )
inlinevirtual

Returns nonzero if receiver is non linear.

Reimplemented from oofem::Material.

Reimplemented in oofem::MPSDamMaterial, oofem::MaxwellChainMaterial, oofem::KelvinChainSolidMaterial, and oofem::KelvinChainMaterial.

Definition at line 214 of file rheoChM.h.

IRResultType oofem::RheoChainMaterial::initializeFrom ( InputRecord *  ir)
virtual

Initializes receiver according to object description stored in input record.

This function is called immediately after creating object using constructor. Input record can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record.

See also
IR_GIVE_FIELD
IR_GIVE_OPTIONAL_FIELD
Parameters
irInput record to initialize from.
Returns
IRResultType
Warning
Stiffness is time dependant, so the variable changes with time.

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::MPSMaterial, oofem::MPSDamMaterial, oofem::Eurocode2CreepMaterial, oofem::B3SolidMaterial, oofem::B3Material, oofem::MaxwellChainMaterial, oofem::KelvinChainSolidMaterial, oofem::KelvinChainMaterial, oofem::CebFip78Material, and oofem::DoublePowerLawMaterial.

Definition at line 583 of file rheoChM.C.

References _IFT_RheoChainMaterial_alphaOne, _IFT_RheoChainMaterial_alphaTwo, _IFT_RheoChainMaterial_begoftimeofinterest, _IFT_RheoChainMaterial_endoftimeofinterest, _IFT_RheoChainMaterial_lattice, _IFT_RheoChainMaterial_n, _IFT_RheoChainMaterial_preCastingTimeMat, _IFT_RheoChainMaterial_relmatage, _IFT_RheoChainMaterial_talpha, _IFT_RheoChainMaterial_timefactor, alphaOne, alphaTwo, begOfTimeOfInterest, computeCharTimes(), discreteTimeScale, endOfTimeOfInterest, generateLogTimeScale(), giveEndOfTimeOfInterest(), oofem::InputRecord::hasField(), oofem::StructuralMaterial::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_OK, lattice, MNC_NPOINTS, nu, preCastingTimeMat, relMatAge, talpha, and timeFactor.

Referenced by oofem::KelvinChainMaterial::initializeFrom(), oofem::KelvinChainSolidMaterial::initializeFrom(), and oofem::MaxwellChainMaterial::initializeFrom().

virtual bool oofem::RheoChainMaterial::isActivated ( TimeStep *  tStep)
inlinevirtual
Parameters
tStepTime step to check activity for.
Returns
True if material is activated for given solution step.

Reimplemented from oofem::Material.

Definition at line 292 of file rheoChM.h.

References oofem::Material::isActivated().

Referenced by oofem::KelvinChainSolidMaterial::computeHiddenVars(), oofem::Eurocode2CreepMaterial::giveEModulus(), giveRealStressVector(), oofem::Eurocode2CreepMaterial::giveRealStressVector(), oofem::MPSDamMaterial::giveRealStressVector(), and oofem::Eurocode2CreepMaterial::giveShrinkageStrainVector().

contextIOResultType oofem::RheoChainMaterial::restoreIPContext ( DataStream &  stream,
ContextMode  mode,
GaussPoint *  gp 
)
virtual

Reads integration point state to output stream.

Parameters
streamOutput stream.
modeDetermines amount of info required in stream (state, definition, ...).
gpintegration point.
Returns
contextIOResultType.
Exceptions
throwsan ContextIOERR exception if error encountered.

Reimplemented from oofem::Material.

Definition at line 679 of file rheoChM.C.

References oofem::CIO_OK, giveLinearElasticMaterial(), oofem::Material::restoreIPContext(), and THROW_CIOERR.

contextIOResultType oofem::RheoChainMaterial::saveIPContext ( DataStream &  stream,
ContextMode  mode,
GaussPoint *  gp 
)
virtual

Stores integration point state to output stream.

Parameters
streamOutput stream.
modeDetermines amount of info required in stream (state, definition, ...).
gpintegration point.
Returns
contextIOResultType.
Exceptions
throwsan ContextIOERR exception if error encountered.

Reimplemented from oofem::Material.

Definition at line 659 of file rheoChM.C.

References oofem::CIO_OK, giveLinearElasticMaterial(), oofem::Material::saveIPContext(), and THROW_CIOERR.

void oofem::RheoChainMaterial::updateEparModuli ( double  tPrime,
GaussPoint *  gp,
TimeStep *  tStep 
)
protectedvirtual

Update of partial moduli of individual chain units.

Reimplemented in oofem::B3SolidMaterial.

Definition at line 313 of file rheoChM.C.

References computeCharCoefficients(), EparVal, EparValTime, and TIME_DIFF.

Referenced by oofem::MaxwellChainMaterial::computeHiddenVars(), oofem::KelvinChainSolidMaterial::giveEigenStrainVector(), oofem::KelvinChainSolidMaterial::giveEModulus(), oofem::KelvinChainMaterial::giveEModulus(), oofem::MaxwellChainMaterial::giveEModulus(), and oofem::B3SolidMaterial::updateEparModuli().

Member Data Documentation

double oofem::RheoChainMaterial::alphaOne
protected

Parameters for the lattice model.

Definition at line 144 of file rheoChM.h.

Referenced by initializeFrom().

double oofem::RheoChainMaterial::alphaTwo
protected

Definition at line 144 of file rheoChM.h.

Referenced by initializeFrom().

double oofem::RheoChainMaterial::begOfTimeOfInterest
protected

Time from which the model should give a good approximation. Optional field. Default value is 0.1 [day].

Definition at line 149 of file rheoChM.h.

Referenced by oofem::B3SolidMaterial::computeCharTimes(), oofem::Eurocode2CreepMaterial::computeCharTimes(), computeCharTimes(), and initializeFrom().

FloatArray oofem::RheoChainMaterial::charTimes
protected

Characteristic times of individual units (relaxation or retardation times).

Definition at line 158 of file rheoChM.h.

Referenced by oofem::B3SolidMaterial::computeCharCoefficients(), oofem::B3SolidMaterial::computeCharTimes(), oofem::Eurocode2CreepMaterial::computeCharTimes(), computeCharTimes(), and giveCharTime().

FloatArray oofem::RheoChainMaterial::discreteTimeScale
protected

Times at which the errors are evaluated if the least-square method is used.

Definition at line 160 of file rheoChM.h.

Referenced by giveDiscreteTimes(), and initializeFrom().

double oofem::RheoChainMaterial::endOfTimeOfInterest
protected

Time (age???) up to which the model should give a good approximation.

Definition at line 151 of file rheoChM.h.

Referenced by oofem::B3SolidMaterial::computeCharTimes(), oofem::Eurocode2CreepMaterial::computeCharTimes(), giveEndOfTimeOfInterest(), and initializeFrom().

FloatArray oofem::RheoChainMaterial::EparVal
protected

Partial moduli of individual units.

Definition at line 155 of file rheoChM.h.

Referenced by oofem::KelvinChainSolidMaterial::giveEigenStrainVector(), oofem::KelvinChainSolidMaterial::giveEModulus(), giveEparModulus(), oofem::B3SolidMaterial::updateEparModuli(), and updateEparModuli().

double oofem::RheoChainMaterial::EparValTime
protected

Time for which the partial moduli of individual units have been evaluated.

Definition at line 146 of file rheoChM.h.

Referenced by RheoChainMaterial(), and updateEparModuli().

bool oofem::RheoChainMaterial::lattice
protected

Definition at line 140 of file rheoChM.h.

Referenced by giveLinearElasticMaterial(), and initializeFrom().

LinearElasticMaterial* oofem::RheoChainMaterial::linearElasticMaterial
protected

Associated linearElasticMaterial, with E = 1.

Definition at line 153 of file rheoChM.h.

Referenced by giveLinearElasticMaterial(), RheoChainMaterial(), and ~RheoChainMaterial().

double oofem::RheoChainMaterial::nu
protected

Poisson's ratio (assumed to be constant, unaffected by creep).

Definition at line 142 of file rheoChM.h.

Referenced by giveLinearElasticMaterial(), and initializeFrom().

int oofem::RheoChainMaterial::nUnits
protected

Number of (Maxwell or Kelvin) units in the rheologic chain.

Definition at line 136 of file rheoChM.h.

Referenced by oofem::KelvinChainMaterial::computeCharCoefficients(), oofem::MaxwellChainMaterial::computeCharCoefficients(), oofem::B3SolidMaterial::computeCharCoefficients(), oofem::Eurocode2CreepMaterial::computeCharCoefficients(), oofem::B3SolidMaterial::computeCharTimes(), oofem::Eurocode2CreepMaterial::computeCharTimes(), computeCharTimes(), oofem::KelvinChainSolidMaterial::computeHiddenVars(), oofem::MaxwellChainMaterial::computeHiddenVars(), oofem::KelvinChainMaterial::computeHiddenVars(), oofem::KelvinChainMaterial::CreateStatus(), oofem::KelvinChainSolidMaterial::CreateStatus(), oofem::MaxwellChainMaterial::CreateStatus(), oofem::B3SolidMaterial::CreateStatus(), oofem::Eurocode2CreepMaterial::CreateStatus(), oofem::MPSDamMaterial::CreateStatus(), CreateStatus(), oofem::KelvinChainMaterial::giveEigenStrainVector(), oofem::KelvinChainSolidMaterial::giveEigenStrainVector(), oofem::MaxwellChainMaterial::giveEigenStrainVector(), oofem::KelvinChainSolidMaterial::giveEModulus(), oofem::KelvinChainMaterial::giveEModulus(), oofem::MaxwellChainMaterial::giveEModulus(), and RheoChainMaterial().

int oofem::RheoChainMaterial::preCastingTimeMat
protected

Stiffness at time less than casting time - optional parameter, negative by default.

Definition at line 172 of file rheoChM.h.

Referenced by give1dStressStiffMtrx(), give2dLatticeStiffMtrx(), give3dLatticeStiffMtrx(), give3dMaterialStiffnessMatrix(), givePlaneStrainStiffMtrx(), givePlaneStressStiffMtrx(), giveRealStressVector(), initializeFrom(), and RheoChainMaterial().

double oofem::RheoChainMaterial::relMatAge
protected

Physical age of the material at simulation time = 0.

Definition at line 138 of file rheoChM.h.

Referenced by computeCharTimes(), oofem::Eurocode2CreepMaterial::computeEquivalentMaturity(), oofem::B3SolidMaterial::computeFlowTermViscosity(), oofem::MaxwellChainMaterial::computeHiddenVars(), oofem::B3Material::computeShrinkageStrainVector(), oofem::B3SolidMaterial::computeSolidifiedVolume(), oofem::B3Material::computeTotalAverageShrinkageStrainVector(), oofem::B3SolidMaterial::computeTotalAverageShrinkageStrainVector(), oofem::KelvinChainMaterial::giveEModulus(), oofem::MaxwellChainMaterial::giveEModulus(), oofem::B3SolidMaterial::giveEModulus(), oofem::Eurocode2CreepMaterial::giveEModulus(), oofem::Eurocode2CreepMaterial::giveShrinkageStrainVector(), initializeFrom(), and RheoChainMaterial().

double oofem::RheoChainMaterial::talpha
protected

thermal dilatation coeff.

Definition at line 134 of file rheoChM.h.

Referenced by giveThermalDilatationVector(), and initializeFrom().

double oofem::RheoChainMaterial::timeFactor
protected

Scaling factor transforming the simulation time units into days (gives the number of simulation time units in one day, e.g.

86400 if the simulation works with seconds as the time units)

Definition at line 167 of file rheoChM.h.

Referenced by oofem::Eurocode2CreepMaterial::computeCharCoefficients(), oofem::Eurocode2CreepMaterial::computeCharTimes(), oofem::Eurocode2CreepMaterial::computeConcreteStrengthAtAge(), oofem::Eurocode2CreepMaterial::computeCreepCoefficient(), oofem::B3SolidMaterial::computeFlowTermViscosity(), oofem::MaxwellChainMaterial::computeHiddenVars(), oofem::B3SolidMaterial::computeMicroPrestress(), oofem::Eurocode2CreepMaterial::computeRetardationTimeCorrection(), oofem::B3Material::computeShrinkageStrainVector(), oofem::B3SolidMaterial::computeSolidifiedVolume(), oofem::B3Material::computeTotalAverageShrinkageStrainVector(), oofem::B3SolidMaterial::computeTotalAverageShrinkageStrainVector(), oofem::Eurocode2CreepMaterial::evaluateSpectrumAt(), oofem::MaxwellChainMaterial::giveEigenStrainVector(), oofem::B3SolidMaterial::giveEigenStrainVector(), oofem::MaxwellChainMaterial::giveEModulus(), oofem::B3SolidMaterial::giveEModulus(), giveEndOfTimeOfInterest(), oofem::Eurocode2CreepMaterial::giveShrinkageStrainVector(), and initializeFrom().

The documentation for this class was generated from the following files:
This page is part of the OOFEM documentation. Copyright (c) 2011 Borek Patzak
Project e-mail: info@oofem.org
Generated at Tue Jan 2 2018 20:07:41 for OOFEM by doxygen 1.8.11 written by Dimitri van Heesch, © 1997-2011