OOFEM 3.0
Loading...
Searching...
No Matches
oofem::Eurocode2CreepMaterial Class Reference

#include <eurocode2creep.h>

Inheritance diagram for oofem::Eurocode2CreepMaterial:
Collaboration diagram for oofem::Eurocode2CreepMaterial:

Public Member Functions

 Eurocode2CreepMaterial (int n, Domain *d)
void giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) const override
void giveShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override
const char * giveClassName () const override
const char * giveInputRecordName () const override
void initializeFrom (InputRecord &ir) override
std::unique_ptr< MaterialStatusCreateStatus (GaussPoint *gp) const override
virtual double computeConcreteStrengthAtAge (double age) const
 evaluates concrete strength at given age
virtual double computeMeanElasticModulusAtAge (double age) const
 evaluates concrete mean elastic modulus at given age
double computeCreepFunction (double t, double t_prime, GaussPoint *gp, TimeStep *tStep) const override
 Evaluation of the compliance function.
virtual double computeCreepCoefficient (double t, double t_prime, GaussPoint *gp, TimeStep *tStep) const
 Evaluation of the compliance function (according to appendix B from the EC).
Public Member Functions inherited from oofem::KelvinChainMaterial
 KelvinChainMaterial (int n, Domain *d)
const char * giveClassName () const override
void initializeFrom (InputRecord &ir) override
void giveShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override
void giveEigenStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override
std::unique_ptr< MaterialStatusCreateStatus (GaussPoint *gp) const override
void giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) const override
void computeHiddenVars (GaussPoint *gp, TimeStep *tStep) const
Public Member Functions inherited from oofem::RheoChainMaterial
 RheoChainMaterial (int n, Domain *d)
virtual ~RheoChainMaterial ()
void giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) const override
FloatArrayF< 6 > giveRealStressVector_3d (const FloatArrayF< 6 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
FloatArrayF< 4 > giveRealStressVector_PlaneStrain (const FloatArrayF< 4 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_3d.
FloatArrayF< 3 > giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 1 > giveRealStressVector_1d (const FloatArrayF< 1 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 2 > giveRealStressVector_2dBeamLayer (const FloatArrayF< 2 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 5 > giveRealStressVector_PlateLayer (const FloatArrayF< 5 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 6 > giveThermalDilatationVector (GaussPoint *gp, TimeStep *tStep) const override
bool hasMaterialModeCapability (MaterialMode mode) const override
bool hasCastingTimeSupport () const override
const char * giveClassName () const override
void initializeFrom (InputRecord &ir) override
int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
void saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp) override
void restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp) override
FloatMatrixF< 6, 6 > give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 3, 3 > givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 4, 4 > givePlaneStrainStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 1, 1 > give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
FloatArray computeStressIndependentStrainVector (GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override
std::unique_ptr< MaterialStatusCreateStatus (GaussPoint *gp) const override
double giveAlphaOne () const
double giveAlphaTwo () const
double givePoissonsRatio () const
 Returns Poisson's ratio.
bool isActivated (TimeStep *tStep) const override
 Extended meaning: returns true if the material is cast (target time > casting time) or the precasing time mat is defined.
virtual double giveEquivalentTime (GaussPoint *gp, TimeStep *tStep) const
 By default returns equivalent time in the middle of the time step.
Public Member Functions inherited from oofem::StructuralMaterial
 StructuralMaterial (int n, Domain *d)
bool hasMaterialModeCapability (MaterialMode mode) const override
const char * giveClassName () const override
void initializeFrom (InputRecord &ir) override
void giveInputRecord (DynamicInputRecord &input) override
void giveCharacteristicMatrix (FloatMatrix &answer, MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const override
 Returns characteristic matrix of the receiver.
void giveCharacteristicVector (FloatArray &answer, FloatArray &flux, MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const override
 Returns characteristic vector of the receiver.
virtual void giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArray giveRealStressVector_StressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint *gp, TimeStep *tStep) const
 Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArray giveRealStressVector_ShellStressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArrayF< 2 > giveRealStressVector_Warping (const FloatArrayF< 2 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 > giveRealStressVector_Fiber (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 > giveRealStressVector_2dPlateSubSoil (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation is not provided.
virtual FloatArrayF< 6 > giveRealStressVector_3dBeamSubSoil (const FloatArrayF< 6 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArrayF< 9 > giveFirstPKStressVector_3d (const FloatArrayF< 9 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
virtual FloatArrayF< 5 > giveFirstPKStressVector_PlaneStrain (const FloatArrayF< 5 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveFirstPKStressVector_3d.
virtual FloatArray giveFirstPKStressVector_StressControl (const FloatArray &reducedvF, const IntArray &FControl, GaussPoint *gp, TimeStep *tStep) const
 Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArrayF< 4 > giveFirstPKStressVector_PlaneStress (const FloatArrayF< 4 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveFirstPKStressVector_StressControl.
virtual FloatArrayF< 1 > giveFirstPKStressVector_1d (const FloatArrayF< 1 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveFirstPKStressVector_StressControl.
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 giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
double giveReferenceTemperature ()
FloatArrayF< 6 > computeStressIndependentStrainVector_3d (GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const
virtual void giveStiffnessMatrix_dPdF (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep)
virtual FloatMatrixF< 9, 9 > give3dMaterialStiffnessMatrix_dPdF (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const
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) const
int setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type) override
int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
virtual FloatMatrixF< 4, 4 > givePlaneStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
virtual void givePlaneStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
virtual FloatMatrixF< 5, 5 > givePlaneStrainStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
virtual void givePlaneStrainStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
virtual FloatMatrixF< 1, 1 > give1dStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
virtual void give1dStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
virtual FloatMatrixF< 2, 2 > give2dBeamLayerStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
virtual FloatMatrixF< 5, 5 > givePlateLayerStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
virtual FloatMatrixF< 3, 3 > giveFiberStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
virtual FloatMatrixF< 3, 3 > give2dPlateSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
virtual FloatMatrixF< 6, 6 > give3dBeamSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
Public Member Functions inherited from oofem::Material
 Material (int n, Domain *d)
virtual ~Material ()=default
 Destructor.
virtual bool isCharacteristicMtrxSymmetric (MatResponseMode rMode) const
virtual double giveCharacteristicValue (MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const
 Returns characteristic value of the receiver.
virtual double give (int aProperty, GaussPoint *gp) const
virtual bool hasProperty (int aProperty, GaussPoint *gp) const
virtual void modifyProperty (int aProperty, double value, GaussPoint *gp)
double giveCastingTime () const
void printYourself () override
 Prints receiver state on stdout. Useful for debugging.
int checkConsistency () override
virtual void restoreConsistency (GaussPoint *gp)
virtual int initMaterial (Element *element)
virtual MaterialStatusgiveStatus (GaussPoint *gp) const
virtual int packUnknowns (DataStream &buff, TimeStep *tStep, GaussPoint *ip)
virtual int unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep, GaussPoint *ip)
virtual int estimatePackSize (DataStream &buff, GaussPoint *ip)
virtual double predictRelativeComputationalCost (GaussPoint *gp)
virtual double predictRelativeRedistributionCost (GaussPoint *gp)
virtual void initTempStatus (GaussPoint *gp) const
void saveContext (DataStream &stream, ContextMode mode) override
void restoreContext (DataStream &stream, ContextMode mode) override
Public Member Functions inherited from oofem::FEMComponent
 FEMComponent (int n, Domain *d)
virtual ~FEMComponent ()=default
 Virtual destructor.
DomaingiveDomain () const
virtual void setDomain (Domain *d)
int giveNumber () const
void setNumber (int num)
virtual void updateLocalNumbering (EntityRenumberingFunctor &f)
virtual void initializeFrom (InputRecord &ir, int priority)
virtual void initializeFinish ()
virtual void postInitialize ()
 Performs post initialization steps. Called after all components are created and initialized.
virtual void printOutputAt (FILE *file, TimeStep *tStep)
virtual InterfacegiveInterface (InterfaceType t)
std::string errorInfo (const char *func) const
 Returns string for prepending output (used by error reporting macros).

Protected Types

enum  ec2ShrinkageType { EC2_NoShrinkage , EC2_TotalShrinkage , EC2_DryingShrinkage , EC2_AutogenousShrinkage }
 shrinkage option More...

Protected Member Functions

bool hasIncrementalShrinkageFormulation () const override
double giveEModulus (GaussPoint *gp, TimeStep *tStep) const override
 Evaluation of the incremental modulus.
virtual double computeEquivalentAge (GaussPoint *gp, TimeStep *tStep) const
 implements B.9
virtual double computeEquivalentMaturity (GaussPoint *gp, TimeStep *tStep) const
 implements B.10
void computeElasticityStrengthParams (int)
 sets parameters for elasticity and strength according to formulas from EC
void computeShrinkageParams (int, double)
 sets parameters for shrinkage according to formulas from EC
void computeCreepParams (int, double)
 sets parameters for creep according to formulas from EC
void computeCharTimes () override
 computes retardation times of the aging Kelvin chain
double computeRetardationTimeCorrection (int mu) const
 computes correction factor which multiplies the retardation times
double evaluateSpectrumAt (double tau) const
 evaluates retardation spectrum at given time (t-t')
FloatArray computeCharCoefficients (double tPrime, GaussPoint *gp, TimeStep *tStep) const override
 Evaluation of characteristic moduli of the Kelvin chain.
void computeIncrementOfDryingShrinkageVector (FloatArray &answer, GaussPoint *gp, double tNow, double tThen) const
 computes increment of drying shrinkage - the shrinkage strain is isotropic
void computeIncrementOfAutogenousShrinkageVector (FloatArray &answer, GaussPoint *gp, double tNow, double tThen) const
 computes increment of autogenous shrinkage - the shrinkage strain is isotropic
Protected Member Functions inherited from oofem::KelvinChainMaterial
bool hasIncrementalShrinkageFormulation () const override
FloatArray computeCharCoefficients (double tPrime, GaussPoint *gp, TimeStep *tStep) const override
 Evaluation of the moduli of individual units.
double giveEModulus (GaussPoint *gp, TimeStep *tStep) const override
 Evaluation of the incremental modulus.
Protected Member Functions inherited from oofem::RheoChainMaterial
const FloatArraygiveDiscreteTimes () const
void computeDiscreteRelaxationFunction (FloatArray &answer, const FloatArray &tSteps, double t0, double tr, GaussPoint *gp, TimeStep *tStep) const
void giveUnitComplianceMatrix (FloatMatrix &answer, GaussPoint *gp, TimeStep *tStep) const
 Evaluation of elastic compliance matrix for unit Young's modulus.
void giveUnitStiffnessMatrix (FloatMatrix &answer, GaussPoint *gp, TimeStep *tStep) const
 Evaluation of elastic stiffness matrix for unit Young's modulus.
virtual void updateEparModuli (double tPrime, GaussPoint *gp, TimeStep *tStep) const
 Update of partial moduli of individual chain units.
double giveEparModulus (int iChain) const
 Access to partial modulus of a given unit.
double giveCharTime (int) const
 Access to the characteristic time of a given unit.
virtual double giveCharTimeExponent (int i) const
 Exponent to be used with the char time of a given unit, usually = 1.0.
StructuralMaterialgiveLinearElasticMaterial ()
 Access to the underlying linear elastic material with unit Young's modulus.
double giveEndOfTimeOfInterest ()
 Access to the time up to which the response should be accurate.
void computeTrueStressIndependentStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const

Protected Attributes

double tau1 = 0.
 fixed retardation time of the first unit
double EspringVal = 0.
 stiffness of the zeroth Kelvin unit
double fcm28 = 0.
 mean compressive strength at 28 days default - to be specified in units of the analysis (e.g. 30.e6 + stiffnessFacotr 1. or 30. + stiffnessFactor 1.e6)
double Ecm28 = 0.
 Young's modulus at 28 days default [MPa].
double stiffnessFactor = 0.
 factor unifying stiffnesses (Ecm is predicted from fcm...)
double s = 0.
 parameter determined by cement type
double phi_RH = 0.
 drying creep coefficient
double beta_fcm = 0.
 drying creep coefficient
double beta_H = 0.
 drying creep coefficient
double h0 = 0.
 effective thickness [mm]
double alpha_T_cement = 0.
 influence of cement type on concrete equivalent age, B.9 in EC2
double t0 = 0.
 duration of curing [day, sec, ...]
double begOfDrying = 0.
 age (absolute) when concrete started drying, must be in days!
double kh = 0.
 drying shrinkage coefficient
double eps_cd_0 = 0.
 asymptotic value of drying shrinkage at zero relative humidity, B.11 in EC2
double eps_ca_infty = 0.
 asymptotic value of autogenous shrinakge, 3.12 in EC2
enum oofem::Eurocode2CreepMaterial::ec2ShrinkageType shType = EC2_NoShrinkage
bool retardationSpectrumApproximation = 0.
bool temperatureDependent = 0.
 switch for temperature dependence of concrete maturity (default option is off)
Protected Attributes inherited from oofem::RheoChainMaterial
double talpha = 0.
 thermal dilatation coeff.
int nUnits = 0
 Number of (Maxwell or Kelvin) units in the rheologic chain.
double relMatAge = 0.
 Physical age of the material at castingTime.
bool lattice = false
double nu = 0.
 Poisson's ratio (assumed to be constant, unaffected by creep).
double alphaOne = 0.
 Parameters for the lattice model.
double alphaTwo = 0.
double EparValTime = -1.
 Time for which the partial moduli of individual units have been evaluated.
double begOfTimeOfInterest = 0.
 Time from which the model should give a good approximation. Optional field. Default value is 0.1 [day].
double endOfTimeOfInterest = 0.
 Time (age???) up to which the model should give a good approximation.
StructuralMateriallinearElasticMaterial = nullptr
 Associated linearElasticMaterial, with E = 1.
FloatArray EparVal
 Partial moduli of individual units.
FloatArray charTimes
 Characteristic times of individual units (relaxation or retardation times).
FloatArray discreteTimeScale
 Times at which the errors are evaluated if the least-square method is used.
double timeFactor = 0.
Protected Attributes inherited from oofem::StructuralMaterial
double referenceTemperature = 0.
 Reference temperature (temperature, when material has been built into structure).
MatResponseMode SCStiffMode = TangentStiffness
 stifness mode used in stress control
double SCRelTol = 1.e-3
 relative tolerance for stress control
double SCAbsTol = 1.e-12
 absolute stress tolerance for stress control
int SCMaxiter = 100000
 maximum iterations for stress-control
Protected Attributes inherited from oofem::Material
Dictionary propertyDictionary
double castingTime
int preCastingTimeMat
 Material existing before casting time - optional parameter, zero by default.
Protected Attributes inherited from oofem::FEMComponent
int number
 Component number.
Domaindomain
 Link to domain object, useful for communicating with other FEM components.

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 FloatMatrixF< 9, 9 > convert_dSdE_2_dPdF_3D (const FloatMatrixF< 6, 6 > &dSdE, const FloatArrayF< 6 > &S, const FloatArrayF< 9 > &F)
static FloatMatrixF< 5, 5 > convert_dSdE_2_dPdF_PlaneStrain (const FloatMatrixF< 4, 4 > &dSdE, const FloatArrayF< 4 > &S, const FloatArrayF< 5 > &F)
static FloatMatrixF< 4, 4 > convert_dSdE_2_dPdF_PlaneStress (const FloatMatrixF< 3, 3 > &dSdE, const FloatArrayF< 3 > &S, const FloatArrayF< 4 > &F)
static FloatMatrixF< 1, 1 > convert_dSdE_2_dPdF_1D (const FloatMatrixF< 1, 1 > &dSdE, const FloatArrayF< 1 > &S, const FloatArrayF< 1 > &F)
static void computePrincipalValues (FloatArray &answer, const FloatArray &s, stressStrainPrincMode mode)
 Common functions for convenience.
static FloatArrayF< 3 > computePrincipalValues (const FloatMatrixF< 3, 3 > &s)
static FloatArrayF< 3 > computePrincipalValues (double I1, double I2, double I3)
static void computePrincipalValDir (FloatArray &answer, FloatMatrix &dir, const FloatArray &s, stressStrainPrincMode mode)
static std::pair< FloatArrayF< 3 >, FloatMatrixF< 3, 3 > > computePrincipalValDir (const FloatMatrixF< 3, 3 > &s)
static FloatArrayF< 6 > computeDeviator (const FloatArrayF< 6 > &s)
static std::pair< FloatArrayF< 6 >, double > computeDeviatoricVolumetricSplit (const FloatArrayF< 6 > &s)
static FloatArrayF< 6 > computeDeviatoricVolumetricSum (const FloatArrayF< 6 > &dev, double mean)
static FloatArrayF< 6 > applyDeviatoricElasticCompliance (const FloatArrayF< 6 > &stress, double EModulus, double nu)
static FloatArrayF< 6 > applyDeviatoricElasticCompliance (const FloatArrayF< 6 > &stress, double GModulus)
static FloatArrayF< 6 > applyDeviatoricElasticStiffness (const FloatArrayF< 6 > &strain, double EModulus, double nu)
static FloatArrayF< 6 > applyDeviatoricElasticStiffness (const FloatArrayF< 6 > &strain, double GModulus)
static FloatArrayF< 6 > applyElasticStiffness (const FloatArrayF< 6 > &strain, double EModulus, double nu)
static FloatArrayF< 6 > applyElasticCompliance (const FloatArrayF< 6 > &stress, double EModulus, double nu)
static double computeStressNorm (const FloatArrayF< 6 > &stress)
static double computeFirstInvariant (const FloatArrayF< 6 > &s)
static double computeSecondStressInvariant (const FloatArrayF< 6 > &s)
static double computeThirdStressInvariant (const FloatArrayF< 6 > &s)
static double computeFirstCoordinate (const FloatArrayF< 6 > &s)
static double computeSecondCoordinate (const FloatArrayF< 6 > &s)
static double computeThirdCoordinate (const FloatArrayF< 6 > &s)
static int giveVoigtVectorMask (IntArray &answer, MaterialMode mmode)
static int giveVoigtSymVectorMask (IntArray &answer, MaterialMode mmode)
static void giveInvertedVoigtVectorMask (IntArray &answer, MaterialMode mmode)
static int giveSizeOfVoigtVector (MaterialMode mmode)
static int giveSizeOfVoigtSymVector (MaterialMode mmode)
static void giveFullVectorForm (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
 Converts the reduced symmetric Voigt vector (2nd order tensor) to full form.
static void giveFullVectorFormF (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
 Converts the reduced deformation gradient Voigt vector (2nd order tensor).
static void giveFullSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
 Converts the reduced unsymmetric Voigt vector (2nd order tensor) to full form.
static void giveReducedVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
 Converts the full symmetric Voigt vector (2nd order tensor) to reduced form.
static void giveReducedSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
 Converts the full unsymmetric Voigt vector (2nd order tensor) to reduced form.
static void giveFullSymMatrixForm (FloatMatrix &answer, const FloatMatrix &red, MaterialMode matMode)
 Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form.
static void giveReducedMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
 Converts the full symmetric Voigt matrix (4th order tensor) to reduced form.
static void giveReducedSymMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
 Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form.
static FloatArrayF< 6 > transformStrainVectorTo (const FloatMatrixF< 3, 3 > &base, const FloatArrayF< 6 > &strain, bool transpose=false)
static FloatArrayF< 6 > transformStressVectorTo (const FloatMatrixF< 3, 3 > &base, const FloatArrayF< 6 > &stress, bool transpose=false)
static double computeVonMisesStress (const FloatArray &currentStress)
static double computeVonMisesStress_3D (const FloatArrayF< 6 > &stress)
static double computeVonMisesStress_PlaneStress (const FloatArrayF< 3 > &stress)
static FloatMatrixF< 6, 6 > giveStrainVectorTranformationMtrx (const FloatMatrixF< 3, 3 > &base, bool transpose=false)
static FloatMatrixF< 3, 3 > give2DStrainVectorTranformationMtrx (const FloatMatrixF< 2, 2 > &base, bool transpose=false)
static FloatMatrixF< 6, 6 > giveStressVectorTranformationMtrx (const FloatMatrixF< 3, 3 > &base, bool transpose=false)
static FloatMatrixF< 3, 3 > givePlaneStressVectorTranformationMtrx (const FloatMatrixF< 2, 2 > &base, bool transpose=false)
static void sortPrincDirAndValCloseTo (FloatArray &pVal, FloatMatrix &pDir, const FloatMatrix &toPDir)
Static Public Attributes inherited from oofem::StructuralMaterial
static std::array< std::array< int, 3 >, 3 > vIindex
 Voigt index map.
static std::array< std::array< int, 3 >, 3 > svIndex
 Symmetric Voigt index map.
Static Protected Member Functions inherited from oofem::RheoChainMaterial
static FloatArray generateLogTimeScale (double from, double to, int nsteps)

Detailed Description

This class implements a model for concrete creep and shrinkage according to EuroCode 2 The creep is assumed to be linear (formula from section 3.7 is not considered here) with aging. The linearity assumption is considered OK for stresses below 0.45 of the compressive strength. The shrinkage deformation is additively split into two parts: drying shrinkage $\varepsilon_{sh,d}$ and autogenous shrinkage $\varepsilon_{sh,a}$. A detailed description of the material model is given in the material manual. See the tests "/sm/EC2creep.in" and "/sm/EC2shrinkage.in" for examples.

Definition at line 99 of file eurocode2creep.h.

Member Enumeration Documentation

◆ ec2ShrinkageType

shrinkage option

Enumerator
EC2_NoShrinkage 
EC2_TotalShrinkage 
EC2_DryingShrinkage 
EC2_AutogenousShrinkage 

Definition at line 160 of file eurocode2creep.h.

Constructor & Destructor Documentation

◆ Eurocode2CreepMaterial()

oofem::Eurocode2CreepMaterial::Eurocode2CreepMaterial ( int n,
Domain * d )
inline

Member Function Documentation

◆ computeCharCoefficients()

FloatArray oofem::Eurocode2CreepMaterial::computeCharCoefficients ( double tPrime,
GaussPoint * gp,
TimeStep * tStep ) const
overrideprotectedvirtual

◆ computeCharTimes()

void oofem::Eurocode2CreepMaterial::computeCharTimes ( )
overrideprotectedvirtual

◆ computeConcreteStrengthAtAge()

double oofem::Eurocode2CreepMaterial::computeConcreteStrengthAtAge ( double age) const
virtual

evaluates concrete strength at given age

Definition at line 230 of file eurocode2creep.C.

References fcm28, s, and oofem::RheoChainMaterial::timeFactor.

Referenced by computeMeanElasticModulusAtAge().

◆ computeCreepCoefficient()

double oofem::Eurocode2CreepMaterial::computeCreepCoefficient ( double t,
double t_prime,
GaussPoint * gp,
TimeStep * tStep ) const
virtual

Evaluation of the compliance function (according to appendix B from the EC).

Definition at line 264 of file eurocode2creep.C.

References beta_fcm, beta_H, computeEquivalentAge(), phi_RH, temperatureDependent, and oofem::RheoChainMaterial::timeFactor.

Referenced by computeCreepFunction().

◆ computeCreepFunction()

double oofem::Eurocode2CreepMaterial::computeCreepFunction ( double t,
double t_prime,
GaussPoint * gp,
TimeStep * tStep ) const
overridevirtual

Evaluation of the compliance function.

Implements oofem::RheoChainMaterial.

Definition at line 250 of file eurocode2creep.C.

References computeCreepCoefficient(), computeMeanElasticModulusAtAge(), and Ecm28.

◆ computeCreepParams()

void oofem::Eurocode2CreepMaterial::computeCreepParams ( int cemType,
double henv )
protected

sets parameters for creep according to formulas from EC

Definition at line 172 of file eurocode2creep.C.

References alpha_T_cement, beta_fcm, beta_H, fcm28, h0, oofem::min(), OOFEM_ERROR, phi_RH, and stiffnessFactor.

Referenced by initializeFrom().

◆ computeElasticityStrengthParams()

void oofem::Eurocode2CreepMaterial::computeElasticityStrengthParams ( int cemType)
protected

sets parameters for elasticity and strength according to formulas from EC

Definition at line 101 of file eurocode2creep.C.

References Ecm28, fcm28, OOFEM_ERROR, s, and stiffnessFactor.

Referenced by initializeFrom().

◆ computeEquivalentAge()

double oofem::Eurocode2CreepMaterial::computeEquivalentAge ( GaussPoint * gp,
TimeStep * tStep ) const
protectedvirtual

implements B.9

Definition at line 325 of file eurocode2creep.C.

References alpha_T_cement, computeEquivalentMaturity(), and oofem::max().

Referenced by computeCharCoefficients(), and computeCreepCoefficient().

◆ computeEquivalentMaturity()

◆ computeIncrementOfAutogenousShrinkageVector()

void oofem::Eurocode2CreepMaterial::computeIncrementOfAutogenousShrinkageVector ( FloatArray & answer,
GaussPoint * gp,
double tNow,
double tThen ) const
protected

computes increment of autogenous shrinkage - the shrinkage strain is isotropic

Definition at line 619 of file eurocode2creep.C.

References oofem::FloatArray::at(), eps_ca_infty, and oofem::GaussPoint::giveMaterialMode().

Referenced by giveShrinkageStrainVector().

◆ computeIncrementOfDryingShrinkageVector()

void oofem::Eurocode2CreepMaterial::computeIncrementOfDryingShrinkageVector ( FloatArray & answer,
GaussPoint * gp,
double tNow,
double tThen ) const
protected

computes increment of drying shrinkage - the shrinkage strain is isotropic

Definition at line 591 of file eurocode2creep.C.

References oofem::FloatArray::at(), eps_cd_0, oofem::GaussPoint::giveMaterialMode(), h0, and kh.

Referenced by giveShrinkageStrainVector().

◆ computeMeanElasticModulusAtAge()

double oofem::Eurocode2CreepMaterial::computeMeanElasticModulusAtAge ( double age) const
virtual

evaluates concrete mean elastic modulus at given age

Definition at line 238 of file eurocode2creep.C.

References computeConcreteStrengthAtAge(), Ecm28, and fcm28.

Referenced by computeCreepFunction(), and giveEModulus().

◆ computeRetardationTimeCorrection()

double oofem::Eurocode2CreepMaterial::computeRetardationTimeCorrection ( int mu) const
protected

computes correction factor which multiplies the retardation times

Definition at line 451 of file eurocode2creep.C.

References beta_H, tau1, and oofem::RheoChainMaterial::timeFactor.

Referenced by computeCharTimes().

◆ computeShrinkageParams()

void oofem::Eurocode2CreepMaterial::computeShrinkageParams ( int cemType,
double henv )
protected

sets parameters for shrinkage according to formulas from EC

Definition at line 121 of file eurocode2creep.C.

References EC2_AutogenousShrinkage, EC2_DryingShrinkage, EC2_TotalShrinkage, eps_ca_infty, eps_cd_0, fcm28, h0, kh, oofem::max(), OOFEM_ERROR, shType, and stiffnessFactor.

Referenced by initializeFrom().

◆ CreateStatus()

std::unique_ptr< MaterialStatus > oofem::Eurocode2CreepMaterial::CreateStatus ( GaussPoint * gp) const
overridevirtual

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.

Definition at line 647 of file eurocode2creep.C.

References oofem::RheoChainMaterial::nUnits.

◆ evaluateSpectrumAt()

double oofem::Eurocode2CreepMaterial::evaluateSpectrumAt ( double tau) const
protected

evaluates retardation spectrum at given time (t-t')

Definition at line 461 of file eurocode2creep.C.

References beta_H, and oofem::RheoChainMaterial::timeFactor.

Referenced by computeCharCoefficients().

◆ giveClassName()

const char * oofem::Eurocode2CreepMaterial::giveClassName ( ) const
inlineoverridevirtual
Returns
Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 180 of file eurocode2creep.h.

◆ giveEModulus()

double oofem::Eurocode2CreepMaterial::giveEModulus ( GaussPoint * gp,
TimeStep * tStep ) const
overrideprotectedvirtual

◆ giveInputRecordName()

const char * oofem::Eurocode2CreepMaterial::giveInputRecordName ( ) const
inlineoverridevirtual
Returns
Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 181 of file eurocode2creep.h.

References _IFT_Eurocode2CreepMaterial_Name.

◆ giveRealStressVector()

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

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.

Definition at line 654 of file eurocode2creep.C.

References oofem::RheoChainMaterial::isActivated().

◆ giveShrinkageStrainVector()

void oofem::Eurocode2CreepMaterial::giveShrinkageStrainVector ( FloatArray & answer,
GaussPoint * gp,
TimeStep * tStep,
ValueModeType mode ) const
overridevirtual

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 from oofem::RheoChainMaterial.

Definition at line 531 of file eurocode2creep.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::Material::castingTime, computeIncrementOfAutogenousShrinkageVector(), computeIncrementOfDryingShrinkageVector(), EC2_AutogenousShrinkage, EC2_DryingShrinkage, EC2_NoShrinkage, EC2_TotalShrinkage, oofem::GaussPoint::giveMaterialMode(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), oofem::RheoChainMaterial::isActivated(), oofem::max(), OOFEM_ERROR, oofem::RheoChainMaterial::relMatAge, oofem::FloatArray::resize(), shType, t0, oofem::RheoChainMaterial::timeFactor, and oofem::FloatArray::zero().

◆ hasIncrementalShrinkageFormulation()

bool oofem::Eurocode2CreepMaterial::hasIncrementalShrinkageFormulation ( ) const
inlineoverrideprotectedvirtual

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 from oofem::RheoChainMaterial.

Definition at line 201 of file eurocode2creep.h.

◆ initializeFrom()

void oofem::Eurocode2CreepMaterial::initializeFrom ( InputRecord & ir)
overridevirtual

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. Note that initializeFrom may be called mutiple times.

See also
IR_GIVE_FIELD
IR_GIVE_OPTIONAL_FIELD
Parameters
irInput record to initialize from.
priorityPriority of the input record. This is used to determine the order of initialization

Reimplemented from oofem::FEMComponent.

Definition at line 48 of file eurocode2creep.C.

References _IFT_Eurocode2CreepMaterial_cemType, _IFT_Eurocode2CreepMaterial_fcm28, _IFT_Eurocode2CreepMaterial_h0, _IFT_Eurocode2CreepMaterial_henv, _IFT_Eurocode2CreepMaterial_shType, _IFT_Eurocode2CreepMaterial_spectrum, _IFT_Eurocode2CreepMaterial_stiffnessFactor, _IFT_Eurocode2CreepMaterial_t0, _IFT_Eurocode2CreepMaterial_temperatureDependent, computeCreepParams(), computeElasticityStrengthParams(), computeShrinkageParams(), fcm28, h0, oofem::InputRecord::hasField(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, OOFEM_ERROR, retardationSpectrumApproximation, shType, stiffnessFactor, t0, and temperatureDependent.

Member Data Documentation

◆ alpha_T_cement

double oofem::Eurocode2CreepMaterial::alpha_T_cement = 0.
protected

influence of cement type on concrete equivalent age, B.9 in EC2

Definition at line 138 of file eurocode2creep.h.

Referenced by computeCreepParams(), and computeEquivalentAge().

◆ begOfDrying

double oofem::Eurocode2CreepMaterial::begOfDrying = 0.
protected

age (absolute) when concrete started drying, must be in days!

Definition at line 146 of file eurocode2creep.h.

◆ beta_fcm

double oofem::Eurocode2CreepMaterial::beta_fcm = 0.
protected

drying creep coefficient

Definition at line 129 of file eurocode2creep.h.

Referenced by computeCharCoefficients(), computeCreepCoefficient(), and computeCreepParams().

◆ beta_H

double oofem::Eurocode2CreepMaterial::beta_H = 0.
protected

drying creep coefficient

Definition at line 132 of file eurocode2creep.h.

Referenced by computeCreepCoefficient(), computeCreepParams(), computeRetardationTimeCorrection(), and evaluateSpectrumAt().

◆ Ecm28

double oofem::Eurocode2CreepMaterial::Ecm28 = 0.
protected

Young's modulus at 28 days default [MPa].

Definition at line 116 of file eurocode2creep.h.

Referenced by computeCharCoefficients(), computeCreepFunction(), computeElasticityStrengthParams(), and computeMeanElasticModulusAtAge().

◆ eps_ca_infty

double oofem::Eurocode2CreepMaterial::eps_ca_infty = 0.
protected

asymptotic value of autogenous shrinakge, 3.12 in EC2

Definition at line 157 of file eurocode2creep.h.

Referenced by computeIncrementOfAutogenousShrinkageVector(), and computeShrinkageParams().

◆ eps_cd_0

double oofem::Eurocode2CreepMaterial::eps_cd_0 = 0.
protected

asymptotic value of drying shrinkage at zero relative humidity, B.11 in EC2

Definition at line 152 of file eurocode2creep.h.

Referenced by computeIncrementOfDryingShrinkageVector(), and computeShrinkageParams().

◆ EspringVal

double oofem::Eurocode2CreepMaterial::EspringVal = 0.
mutableprotected

stiffness of the zeroth Kelvin unit

Definition at line 109 of file eurocode2creep.h.

Referenced by computeCharCoefficients(), and giveEModulus().

◆ fcm28

double oofem::Eurocode2CreepMaterial::fcm28 = 0.
protected

mean compressive strength at 28 days default - to be specified in units of the analysis (e.g. 30.e6 + stiffnessFacotr 1. or 30. + stiffnessFactor 1.e6)

Definition at line 113 of file eurocode2creep.h.

Referenced by computeConcreteStrengthAtAge(), computeCreepParams(), computeElasticityStrengthParams(), computeMeanElasticModulusAtAge(), computeShrinkageParams(), and initializeFrom().

◆ h0

double oofem::Eurocode2CreepMaterial::h0 = 0.
protected

effective thickness [mm]

Definition at line 135 of file eurocode2creep.h.

Referenced by computeCreepParams(), computeIncrementOfDryingShrinkageVector(), computeShrinkageParams(), and initializeFrom().

◆ kh

double oofem::Eurocode2CreepMaterial::kh = 0.
protected

drying shrinkage coefficient

Definition at line 149 of file eurocode2creep.h.

Referenced by computeIncrementOfDryingShrinkageVector(), and computeShrinkageParams().

◆ phi_RH

double oofem::Eurocode2CreepMaterial::phi_RH = 0.
protected

drying creep coefficient

Definition at line 126 of file eurocode2creep.h.

Referenced by computeCharCoefficients(), computeCreepCoefficient(), and computeCreepParams().

◆ retardationSpectrumApproximation

bool oofem::Eurocode2CreepMaterial::retardationSpectrumApproximation = 0.
protected

If true, analysis of retardation spectrum is used for evaluation of Kelvin units moduli If false, least-squares method is used for evaluation of Kelvin units moduli (default)

Definition at line 166 of file eurocode2creep.h.

Referenced by computeCharCoefficients(), computeCharTimes(), giveEModulus(), and initializeFrom().

◆ s

double oofem::Eurocode2CreepMaterial::s = 0.
protected

parameter determined by cement type

Definition at line 122 of file eurocode2creep.h.

Referenced by computeConcreteStrengthAtAge(), and computeElasticityStrengthParams().

◆ shType

◆ stiffnessFactor

double oofem::Eurocode2CreepMaterial::stiffnessFactor = 0.
protected

factor unifying stiffnesses (Ecm is predicted from fcm...)

Definition at line 119 of file eurocode2creep.h.

Referenced by computeCreepParams(), computeElasticityStrengthParams(), computeShrinkageParams(), and initializeFrom().

◆ t0

double oofem::Eurocode2CreepMaterial::t0 = 0.
protected

duration of curing [day, sec, ...]

Definition at line 143 of file eurocode2creep.h.

Referenced by giveShrinkageStrainVector(), and initializeFrom().

◆ tau1

double oofem::Eurocode2CreepMaterial::tau1 = 0.
protected

fixed retardation time of the first unit

Definition at line 106 of file eurocode2creep.h.

Referenced by computeCharCoefficients(), computeCharTimes(), and computeRetardationTimeCorrection().

◆ temperatureDependent

bool oofem::Eurocode2CreepMaterial::temperatureDependent = 0.
protected

switch for temperature dependence of concrete maturity (default option is off)

Definition at line 169 of file eurocode2creep.h.

Referenced by computeCharCoefficients(), computeCreepCoefficient(), and initializeFrom().


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

This page is part of the OOFEM-3.0 documentation. Copyright Copyright (C) 1994-2025 Borek Patzak Bořek Patzák
Project e-mail: oofem@fsv.cvut.cz
Generated at for OOFEM by doxygen 1.15.0 written by Dimitri van Heesch, © 1997-2011