OOFEM
2.4
OOFEM.org - Object Oriented Finite Element Solver
|
This class implements the extended B3 model for concrete creep and shrinkage based on the microprestress-solidification theory. More...
#include <mps.h>
Public Member Functions | |
MPSMaterial (int n, Domain *d) | |
virtual | ~MPSMaterial () |
virtual const char * | giveInputRecordName () const |
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 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 | 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 double | computeCreepFunction (double t, double t_prime, GaussPoint *gp, TimeStep *tStep) |
Evaluation of the basic creep compliance function - can be used to compute elastic modulus in derived damage material. More... | |
virtual MaterialStatus * | CreateStatus (GaussPoint *gp) const |
Creates new copy of associated status and inserts it into given integration point. More... | |
Public Member Functions inherited from oofem::KelvinChainSolidMaterial | |
KelvinChainSolidMaterial (int n, Domain *d) | |
virtual | ~KelvinChainSolidMaterial () |
void | computeHiddenVars (GaussPoint *gp, TimeStep *tStep) |
virtual int | hasNonLinearBehaviour () |
Returns nonzero if receiver is non linear. More... | |
Public Member Functions inherited from oofem::RheoChainMaterial | |
RheoChainMaterial (int n, Domain *d) | |
virtual | ~RheoChainMaterial () |
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 int | hasMaterialModeCapability (MaterialMode mode) |
Tests if material supports material mode. More... | |
virtual int | hasCastingTimeSupport () |
Tests if material supports casting time. 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... | |
double | giveAlphaOne () const |
double | giveAlphaTwo () const |
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... | |
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 Types | |
enum | coupledAnalysisType { Basic, MPS_full, MPS_humidity, MPS_temperature } |
Protected Member Functions | |
void | predictParametersFrom (double, double, double, double) |
virtual void | computeCharTimes () |
Evaluation of characteristic times. More... | |
virtual void | computeCharCoefficients (FloatArray &answer, double, GaussPoint *gp, TimeStep *tStep) |
Evaluation of characteristic moduli of the non-aging Kelvin chain. More... | |
virtual double | giveEModulus (GaussPoint *gp, TimeStep *tStep) |
Evaluation of the incremental modulus. More... | |
virtual double | computeSolidifiedVolume (GaussPoint *gp, TimeStep *tStep) |
Evaluation of the relative volume of the solidified material. More... | |
virtual double | computeBetaMu (GaussPoint *gp, TimeStep *tStep, int Mu) |
factors for exponential algorithm More... | |
virtual double | computeLambdaMu (GaussPoint *gp, TimeStep *tStep, int Mu) |
double | computeFlowTermViscosity (GaussPoint *gp, TimeStep *tStep) |
Evaluation of the flow term viscosity. More... | |
double | giveInitViscosity (TimeStep *tStep) |
Returns initial value of the flow term viscosity. 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 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... | |
void | computePointShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) |
Evaluation of the shrinkageStrainVector - shrinkage is fully dependent on humidity rate in given GP. More... | |
void | computeFibAutogenousShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) |
Evaluation of the autogenousShrinkageStrainVector according to fib MC 2010 - autogenous shrinkage is fully dependent on the equivalent age at given GP. More... | |
void | computeB4AutogenousShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) |
Evaluation of the autogenousShrinkageStrainVector according to Bazant's B4 model. In the model the evolution depends on temperature adjusted age, here on equivalent age (additional humidity influence) More... | |
double | giveHumidity (GaussPoint *gp, TimeStep *tStep, int option) |
Gives value of humidity at given GP and timestep option = 0 ... More... | |
double | giveTemperature (GaussPoint *gp, TimeStep *tStep, int option) |
Gives value of temperature at given GP and timestep option = 0 ... More... | |
double | computePsiR (GaussPoint *gp, TimeStep *tStep, int option) |
Evaluation of the factor transforming real time to reduced time (effect on the flow term) option = 0 ... More... | |
double | computePsiS (GaussPoint *gp, TimeStep *tStep) |
Evaluation of the factor transforming real time to reduced time (effect on the evolution of microprestress) More... | |
double | computePsiE (GaussPoint *gp, TimeStep *tStep) |
Evaluation of the factor transforming real time to equivalent time (effect on the solidified volume) More... | |
double | computeEquivalentTime (GaussPoint *gp, TimeStep *tStep, int option) |
Computes equivalent time at given time step and GP. More... | |
Protected Member Functions inherited from oofem::RheoChainMaterial | |
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... | |
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... | |
Protected Attributes | |
double | t0 |
age when temperature or humidity starts to change More... | |
double | q1 |
compliances of the B3 model More... | |
double | q2 |
double | q3 |
double | q4 |
double | lambda0 |
constant equal to one day in time units of analysis (eg. 86400 if the analysis runs in seconds) More... | |
enum oofem::MPSMaterial::coupledAnalysisType | CoupledAnalysis |
double | EspringVal |
double | kSh |
additional parameters for sorption isotherm (used to compute relative humidity from water content) More... | |
double | muS |
fluidity parameter used in viscosity evolution equation More... | |
double | k3 |
double | kTm |
kTm replaces ln(h) on RHS of the differential equation describing evolution of MPS More... | |
double | kTc |
parameter reducing creep effects of thermal cycling (replaces kTm in such case) More... | |
double | ct |
parameter reducing creep effects of thermal cycling More... | |
double | roomTemperature |
reference room temperature for MPS algorithm [K] More... | |
double | QEtoR |
activation energies More... | |
double | QRtoR |
double | QStoR |
double | alphaE |
parameters that control the effect of humidity on rates of hydration, creep and microprestress relaxation More... | |
double | alphaR |
double | alphaS |
double | p |
exponent in the microprestress/viscosity governing equation More... | |
double | sh_a |
parameters for nonlinear shrinkage function More... | |
double | sh_hC |
double | sh_n |
double | eps_cas0 |
parameter for autogenous shrinkage according to fib MC 2010 More... | |
double | b4_eps_au_infty |
parameters for autogenous shrinkage according to B4 model More... | |
double | b4_tau_au |
double | b4_alpha |
double | b4_r_t |
double | stiffnessFactor |
scaling factor 1. for Pa, 1.e6 for MPa - only for empirical formulas - q1-q4 and ft and gf More... | |
double | temperScaleDifference |
0 for Kelvin, 273.15 for Celsius More... | |
Protected Attributes inherited from oofem::RheoChainMaterial | |
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... | |
Friends | |
class | RankineMPSmat |
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... | |
Static Protected Member Functions inherited from oofem::RheoChainMaterial | |
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... | |
This class implements the extended B3 model for concrete creep and shrinkage based on the microprestress-solidification theory.
The implementation exploits a solidifying Kelvin chain. Creep is affected by variable temperature and humidity. Comparing to other material models for concrete creep MPS material model is unit independent (except parameters corresponding to concrete strength and composition: fc, c, wc, ac ). In the input record is is necessary to specify time factor equal to 1!!!
|
protected |
|
protected |
Evaluation of the autogenousShrinkageStrainVector according to Bazant's B4 model. In the model the evolution depends on temperature adjusted age, here on equivalent age (additional humidity influence)
Definition at line 1270 of file mps.C.
References oofem::FloatArray::at(), oofem::MPSMaterialStatus::giveEquivalentTime(), oofem::TimeStep::giveIntrinsicTime(), oofem::GaussPoint::giveMaterialMode(), oofem::StructuralMaterial::giveReducedSymVectorForm(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::isTheFirstStep(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().
|
protectedvirtual |
factors for exponential algorithm
Reimplemented from oofem::KelvinChainSolidMaterial.
Definition at line 794 of file mps.C.
References oofem::TimeStep::giveTimeIncrement().
|
protectedvirtual |
Evaluation of characteristic moduli of the non-aging Kelvin chain.
Implements oofem::RheoChainMaterial.
Definition at line 653 of file mps.C.
References oofem::FloatArray::at(), oofem::RheoChainMaterialStatus::nUnits, oofem::FloatArray::resize(), and oofem::FloatArray::zero().
|
protectedvirtual |
Evaluation of characteristic times.
Reimplemented from oofem::RheoChainMaterial.
Definition at line 609 of file mps.C.
References oofem::RheoChainMaterial::giveEndOfTimeOfInterest(), and oofem::RheoChainMaterialStatus::nUnits.
|
virtual |
Evaluation of the basic creep compliance function - can be used to compute elastic modulus in derived damage material.
Reimplemented from oofem::KelvinChainSolidMaterial.
Definition at line 582 of file mps.C.
Referenced by oofem::MPSDamMaterial::giveRealStressVector().
|
protected |
Computes equivalent time at given time step and GP.
If option == 0, equivalentTime is evaluated in the middle of the time step (to determine solidified ratio). If option == 1, equivalentTime is evaluated at the end of the time step. (for updating).
Definition at line 1542 of file mps.C.
References oofem::MPSMaterialStatus::giveEquivalentTime(), oofem::TimeStep::giveIntrinsicTime(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::isTheFirstStep(), and OOFEM_ERROR.
Referenced by oofem::MPSDamMaterial::initDamagedFib().
|
protected |
Evaluation of the autogenousShrinkageStrainVector according to fib MC 2010 - autogenous shrinkage is fully dependent on the equivalent age at given GP.
Definition at line 1213 of file mps.C.
References oofem::FloatArray::at(), oofem::MPSMaterialStatus::giveEquivalentTime(), oofem::TimeStep::giveIntrinsicTime(), oofem::GaussPoint::giveMaterialMode(), oofem::StructuralMaterial::giveReducedSymVectorForm(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::isTheFirstStep(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().
|
protected |
Evaluation of the flow term viscosity.
Definition at line 844 of file mps.C.
References oofem::MPSMaterialStatus::giveFlowTermViscosity(), oofem::MPSMaterialStatus::giveFlowTermViscosityTemp(), oofem::TimeStep::giveIntrinsicTime(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), oofem::MPSMaterialStatus::giveTmax(), oofem::TimeStep::isTheFirstStep(), OOFEM_ERROR, oofem::MPSMaterialStatus::setFlowTermViscosityTemp(), and oofem::MPSMaterialStatus::setTmax().
|
protectedvirtual |
Reimplemented from oofem::KelvinChainSolidMaterial.
Definition at line 818 of file mps.C.
References oofem::TimeStep::giveTimeIncrement().
|
protected |
Evaluation of the shrinkageStrainVector - shrinkage is fully dependent on humidity rate in given GP.
Definition at line 1172 of file mps.C.
References oofem::FloatArray::at(), oofem::GaussPoint::giveMaterialMode(), oofem::StructuralMaterial::giveReducedSymVectorForm(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().
|
protected |
Evaluation of the factor transforming real time to equivalent time (effect on the solidified volume)
Definition at line 1515 of file mps.C.
References OOFEM_ERROR.
|
protected |
Evaluation of the factor transforming real time to reduced time (effect on the flow term) option = 0 ...
beginning of the time step option = 1 ... end of the time step option = 2 ... average value
Definition at line 1463 of file mps.C.
References OOFEM_ERROR, and oofem::MPSMaterialStatus::T.
|
protected |
Evaluation of the factor transforming real time to reduced time (effect on the evolution of microprestress)
Definition at line 1489 of file mps.C.
References OOFEM_ERROR.
|
protectedvirtual |
Evaluation of the relative volume of the solidified material.
Implements oofem::KelvinChainSolidMaterial.
Definition at line 772 of file mps.C.
References oofem::TimeStep::giveTargetTime(), and oofem::TimeStep::giveTimeIncrement().
|
virtual |
Creates new copy of associated status and inserts it into given integration point.
gp | Integration point where newly created status will be stored. |
Reimplemented from oofem::KelvinChainSolidMaterial.
Reimplemented in oofem::MPSDamMaterial.
Definition at line 536 of file mps.C.
References oofem::FEMComponent::giveDomain(), oofem::IntegrationPointStatus::gp, oofem::MPSMaterialStatus::MPSMaterialStatus(), and oofem::RheoChainMaterialStatus::nUnits.
|
inlinevirtual |
Reimplemented from oofem::KelvinChainSolidMaterial.
Reimplemented in oofem::MPSDamMaterial.
Definition at line 262 of file mps.h.
References oofem::RheoChainMaterialStatus::giveShrinkageStrainVector(), oofem::IntegrationPointStatus::gp, and oofem::MaterialStatus::initializeFrom().
|
protectedvirtual |
Computes, for the given integration point, the strain vector induced by the stress history (typically creep strain).
answer | Computed strains. |
gp | Integration point. |
tStep | Time step (most models are able to respond only when tStep is the current time step). |
mode | Determines response mode. |
Reimplemented from oofem::KelvinChainSolidMaterial.
Definition at line 1095 of file mps.C.
References oofem::FloatArray::add(), oofem::FloatArray::beProductOf(), oofem::KelvinChainSolidMaterial::giveEigenStrainVector(), oofem::MPSMaterialStatus::giveFlowTermViscosity(), oofem::TimeStep::giveIntrinsicTime(), oofem::TimeStep::giveTimeIncrement(), oofem::RheoChainMaterialStatus::giveViscoelasticStressVector(), oofem::TimeStep::isTheFirstStep(), OOFEM_ERROR, oofem::MPSMaterialStatus::setCreepStrainIncrement(), and oofem::FloatArray::times().
|
protectedvirtual |
Evaluation of the incremental modulus.
Reimplemented from oofem::KelvinChainSolidMaterial.
Definition at line 681 of file mps.C.
References oofem::KelvinChainSolidMaterial::giveEModulus(), oofem::MPSMaterialStatus::giveFlowTermViscosity(), oofem::TimeStep::giveIntrinsicTime(), oofem::MPSMaterialStatus::giveStoredEmodulus(), oofem::MPSMaterialStatus::giveStoredEmodulusFlag(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::isTheFirstStep(), OOFEM_ERROR, oofem::MPSMaterialStatus::setEmodulusFlag(), and oofem::MPSMaterialStatus::storeEmodulus().
|
protected |
Gives value of humidity at given GP and timestep option = 0 ...
beginning of the time step option = 1 ... end of the time step option = 2 ... average values option = 3 ... incremental values
Definition at line 1344 of file mps.C.
References oofem::FloatArray::at(), oofem::Element::computeGlobalCoordinates(), oofem::FEMComponent::domain, oofem::EngngModel::giveContext(), oofem::GaussPoint::giveElement(), oofem::Domain::giveEngngModel(), oofem::FieldManager::giveField(), oofem::EngngModelContext::giveFieldManager(), oofem::MPSMaterialStatus::giveHum(), oofem::MPSMaterialStatus::giveHumIncrement(), oofem::GaussPoint::giveNaturalCoordinates(), OOFEM_ERROR, oofem::MPSMaterialStatus::setHum(), and oofem::MPSMaterialStatus::setHumIncrement().
|
protected |
Returns initial value of the flow term viscosity.
Definition at line 1084 of file mps.C.
References oofem::TimeStep::giveTimeIncrement(), and OOFEM_ERROR.
|
inlinevirtual |
Implements oofem::FEMComponent.
Reimplemented in oofem::MPSDamMaterial.
Definition at line 261 of file mps.h.
References _IFT_MPSMaterial_Name.
|
virtual |
Returns the integration point corresponding value in Reduced form.
answer | Contain corresponding ip value, zero sized if not available. |
gp | Integration point to which the value refers. |
type | Determines the type of internal variable. |
tStep | Determines the time step. |
Reimplemented from oofem::RheoChainMaterial.
Reimplemented in oofem::MPSDamMaterial.
Definition at line 1577 of file mps.C.
References oofem::FloatArray::at(), oofem::MPSMaterialStatus::giveAutogenousShrinkageStrain(), oofem::MPSMaterialStatus::giveCreepStrain(), oofem::MPSMaterialStatus::giveDryingShrinkageStrain(), oofem::StructuralMaterial::giveFullSymVectorForm(), oofem::RheoChainMaterial::giveIPValue(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().
Referenced by oofem::MPSDamMaterial::giveIPValue().
|
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.
answer | Stress vector in reduced form. For large deformations it is treated as the second Piola-Kirchoff stress. |
gp | Integration point. |
reducedStrain | Strain vector in reduced form. For large deformations it is treated as the Green-Lagrange strain. |
tStep | Current time step (most models are able to respond only when tStep is current time step). |
Reimplemented from oofem::KelvinChainSolidMaterial.
Reimplemented in oofem::MPSDamMaterial.
Definition at line 430 of file mps.C.
References oofem::KelvinChainSolidMaterial::giveRealStressVector(), and oofem::MPSMaterialStatus::setEquivalentTimeTemp().
Referenced by oofem::MPSDamMaterial::giveRealStressVector().
|
virtual |
Computes, for the given integration point, the strain vector induced by stress-independent shrinkage.
answer | Returned strain vector. |
gp | Integration point. |
tStep | Time step (most models are able to respond only when tStep is current time step). |
mode | Determines response mode (Total or incremental). |
Reimplemented from oofem::KelvinChainSolidMaterial.
Definition at line 447 of file mps.C.
References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::MPSMaterialStatus::giveAutogenousShrinkageStrain(), oofem::MPSMaterialStatus::giveDryingShrinkageStrain(), oofem::GaussPoint::giveMaterialMode(), oofem::StructuralMaterial::giveReducedSymVectorForm(), oofem::FloatArray::giveSize(), oofem::StructuralMaterial::giveSizeOfVoigtSymVector(), OOFEM_ERROR, oofem::FloatArray::resize(), oofem::MPSMaterialStatus::setTempAutogenousShrinkageStrain(), oofem::MPSMaterialStatus::setTempDryingShrinkageStrain(), and oofem::FloatArray::zero().
|
protected |
Gives value of temperature at given GP and timestep option = 0 ...
beginning of the time step option = 1 ... end of the time step option = 2 ... average values option = 3 ... incremental values
Definition at line 1405 of file mps.C.
References oofem::FloatArray::at(), oofem::Element::computeGlobalCoordinates(), oofem::FEMComponent::domain, oofem::EngngModel::giveContext(), oofem::GaussPoint::giveElement(), oofem::Domain::giveEngngModel(), oofem::FieldManager::giveField(), oofem::EngngModelContext::giveFieldManager(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::MPSMaterialStatus::giveT(), oofem::MPSMaterialStatus::giveTIncrement(), OOFEM_ERROR, oofem::MPSMaterialStatus::setT(), and oofem::MPSMaterialStatus::setTIncrement().
|
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 from oofem::KelvinChainSolidMaterial.
|
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.
ir | Input record to initialize from. |
if sortion parameters are provided then it is assumed that the transport analysis uses moisture content/ratio; if not then it exports relative humidity
flag - if true, external fields and reference temperature are in Celsius
auxiliary parameters for autogenous shrinkage according to B4 model
Reimplemented from oofem::KelvinChainSolidMaterial.
Reimplemented in oofem::MPSDamMaterial.
Definition at line 192 of file mps.C.
References _IFT_MPSMaterial_ac, _IFT_MPSMaterial_alpha_as, _IFT_MPSMaterial_alphae, _IFT_MPSMaterial_alphar, _IFT_MPSMaterial_alphas, _IFT_MPSMaterial_B4_alpha, _IFT_MPSMaterial_B4_cem_type, _IFT_MPSMaterial_B4_eps_au_infty, _IFT_MPSMaterial_B4_r_t, _IFT_MPSMaterial_B4_tau_au, _IFT_MPSMaterial_cc, _IFT_MPSMaterial_coupledanalysistype, _IFT_MPSMaterial_eps_cas0, _IFT_MPSMaterial_fc, _IFT_MPSMaterial_k3, _IFT_MPSMaterial_ksh, _IFT_MPSMaterial_ktc, _IFT_MPSMaterial_ktm, _IFT_MPSMaterial_lambda0, _IFT_MPSMaterial_mode, _IFT_MPSMaterial_mus, _IFT_MPSMaterial_p, _IFT_MPSMaterial_p_tilde, _IFT_MPSMaterial_q1, _IFT_MPSMaterial_q2, _IFT_MPSMaterial_q3, _IFT_MPSMaterial_q4, _IFT_MPSMaterial_qetor, _IFT_MPSMaterial_qrtor, _IFT_MPSMaterial_qstor, _IFT_MPSMaterial_sh_a, _IFT_MPSMaterial_sh_hC, _IFT_MPSMaterial_sh_n, _IFT_MPSMaterial_stiffnessfactor, _IFT_MPSMaterial_t0, _IFT_MPSMaterial_temperInCelsius, _IFT_MPSMaterial_wc, _IFT_RheoChainMaterial_timefactor, oofem::InputRecord::hasField(), oofem::KelvinChainSolidMaterial::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_BAD_FORMAT, oofem::IRRT_OK, OOFEM_ERROR, and OOFEM_WARNING.
Referenced by oofem::MPSDamMaterial::initializeFrom().
|
protected |
Definition at line 546 of file mps.C.
References oofem::FEMComponent::number, and OOFEM_LOG_DEBUG.
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
constant equal to one day in time units of analysis (eg. 86400 if the analysis runs in seconds)
Definition at line 214 of file mps.h.
Referenced by oofem::MPSDamMaterial::computeFractureEnergy(), and oofem::MPSDamMaterial::computeTensileStrength().
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
scaling factor 1. for Pa, 1.e6 for MPa - only for empirical formulas - q1-q4 and ft and gf
Definition at line 251 of file mps.h.
Referenced by oofem::MPSDamMaterial::computeFractureEnergy(), and oofem::MPSDamMaterial::computeTensileStrength().
|
protected |
|
protected |