OOFEM
2.4
OOFEM.org - Object Oriented Finite Element Solver
|
This class implements the B3 model for concrete creep and shrinkage based on the solidification theory. More...
#include <b3solidmat.h>
Public Member Functions | |
B3SolidMaterial (int n, Domain *d) | |
virtual | ~B3SolidMaterial () |
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 const char * | giveClassName () const |
virtual const char * | giveInputRecordName () const |
virtual IRResultType | initializeFrom (InputRecord *ir) |
Initializes receiver according to object description stored in input record. More... | |
virtual MaterialStatus * | CreateStatus (GaussPoint *gp) const |
Creates new copy of associated status and inserts it into given integration point. More... | |
virtual double | computeCreepFunction (double t, double t_prime, GaussPoint *gp, TimeStep *tStep) |
Evaluation of the compliance function of the non-aging solidifying constituent. More... | |
Public Member Functions inherited from oofem::KelvinChainMaterial | |
KelvinChainMaterial (int n, Domain *d) | |
virtual | ~KelvinChainMaterial () |
virtual int | hasNonLinearBehaviour () |
Returns nonzero if receiver is non linear. More... | |
void | computeHiddenVars (GaussPoint *gp, TimeStep *tStep) |
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 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... | |
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 | b3ShModeType { B3_NoShrinkage, B3_AverageShrinkage, B3_PointShrinkage } |
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... | |
void | computeTotalAverageShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) |
void | computePointShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) |
Evaluation of the shrinkageStrainVector. Shrinkage is fully dependent on humidity rate in given GP. More... | |
void | predictParametersFrom (double, double, double, double, double, double, double) |
double | computeSolidifiedVolume (TimeStep *tStep) |
Evaluation of the relative volume of the solidified material. More... | |
double | computeFlowTermViscosity (GaussPoint *gp, TimeStep *tStep) |
Evaluation of the flow term viscosity. More... | |
double | inverse_sorption_isotherm (double w) |
virtual void | computeCharCoefficients (FloatArray &answer, double tPrime, GaussPoint *gp, TimeStep *tStep) |
Evaluation of characteristic moduli of the non-aging Kelvin chain. More... | |
virtual void | updateEparModuli (double tPrime, GaussPoint *gp, TimeStep *tStep) |
Update of partial moduli of individual chain units. More... | |
virtual void | computeCharTimes () |
Evaluation of characteristic times. More... | |
virtual double | giveEModulus (GaussPoint *gp, TimeStep *tStep) |
Evaluation of the incremental modulus. 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... | |
double | computeMicroPrestress (GaussPoint *gp, TimeStep *tStep, int option) |
Computes microprestress at given time step and GP. More... | |
double | giveInitMicroPrestress (void) |
Computes initial value of the MicroPrestress. More... | |
double | giveHumidity (GaussPoint *gp, TimeStep *tStep) |
Computes relative humidity at given time step and GP. More... | |
double | giveHumidityIncrement (GaussPoint *gp, TimeStep *tStep) |
Computes relative humidity increment at given time step and GP. More... | |
Protected Member Functions inherited from oofem::KelvinChainMaterial | |
LinearElasticMaterial * | giveLinearElasticMaterial () |
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... | |
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 |
double | w |
double | E28 |
double | q1 |
double | q2 |
double | q3 |
double | q4 |
double | q5 |
double | lambda0 |
constant equal to one day in time units of analysis (eg. 86400 if the analysis runs in seconds) More... | |
enum oofem::B3SolidMaterial::b3ShModeType | shMode |
double | EpsSinf |
Additional parameters for average cross section shrinkage. More... | |
double | kt |
double | ks |
double | vs |
double | hum |
double | es0 |
Additional parameters for free shrinkage at material point. More... | |
double | r |
double | rprime |
double | at |
double | w_h |
Constant water content (obtained from experiments) w_h [Pedersen, 1990]. More... | |
double | n |
Constant-exponent (obtained from experiments) n [Pedersen, 1990]. More... | |
double | a |
Constant (obtained from experiments) A [Pedersen, 1990]. More... | |
double | EspringVal |
elastic modulus of the aging spring (first member of Kelvin chain if retardation spectrum is used) More... | |
int | EmoduliMode |
If 0, analysis of retardation spectrum is used for evaluation of Kelvin units moduli (default). More... | |
int | MicroPrestress |
If 1, computation exploiting Microprestress solidification theory is done. More... | |
double | c0 |
MPS constant c0 [MPa^-1 * day^-1]. More... | |
double | c1 |
MPS constant c1 (=C1*R*T/M) More... | |
double | tS0 |
MPS tS0 - necessary for the initial value of microprestress (age when the load is applied) More... | |
double | kSh |
MPS shrinkage parameter. Either this or inithum and finalhum must be given in input record. 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... | |
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 B3 model for concrete creep and shrinkage based on the solidification theory.
The implementation exploits a solidifying Kelvin chain.
Definition at line 86 of file b3solidmat.h.
|
protected |
Enumerator | |
---|---|
B3_NoShrinkage | |
B3_AverageShrinkage | |
B3_PointShrinkage |
Definition at line 95 of file b3solidmat.h.
|
inline |
Definition at line 122 of file b3solidmat.h.
|
inlinevirtual |
Definition at line 125 of file b3solidmat.h.
References oofem::RheoChainMaterialStatus::giveShrinkageStrainVector(), and oofem::IntegrationPointStatus::gp.
|
protectedvirtual |
Evaluation of characteristic moduli of the non-aging Kelvin chain.
Reimplemented from oofem::KelvinChainMaterial.
Definition at line 345 of file b3solidmat.C.
References oofem::FloatArray::at(), oofem::RheoChainMaterial::charTimes, oofem::KelvinChainMaterial::computeCharCoefficients(), EmoduliMode, EspringVal, oofem::RheoChainMaterial::giveCharTime(), lambda0, oofem::RheoChainMaterial::nUnits, q2, oofem::FloatArray::resize(), and oofem::FloatArray::zero().
|
protectedvirtual |
Evaluation of characteristic times.
Reimplemented from oofem::RheoChainMaterial.
Definition at line 298 of file b3solidmat.C.
References oofem::FloatArray::at(), oofem::RheoChainMaterial::begOfTimeOfInterest, oofem::RheoChainMaterial::charTimes, oofem::RheoChainMaterial::endOfTimeOfInterest, oofem::RheoChainMaterial::giveEndOfTimeOfInterest(), lambda0, oofem::RheoChainMaterial::nUnits, oofem::FloatArray::resize(), and oofem::FloatArray::zero().
|
virtual |
Evaluation of the compliance function of the non-aging solidifying constituent.
Implements oofem::RheoChainMaterial.
Definition at line 385 of file b3solidmat.C.
|
protected |
Evaluation of the flow term viscosity.
Definition at line 531 of file b3solidmat.C.
References c0, computeMicroPrestress(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), MicroPrestress, OOFEM_ERROR, q4, oofem::RheoChainMaterial::relMatAge, S, and oofem::RheoChainMaterial::timeFactor.
Referenced by giveEigenStrainVector(), and giveEModulus().
|
protected |
Computes microprestress at given time step and GP.
gp | Gauss point to compute at. |
tStep | Time step to compute for. |
option | If 0, microprestress is evaluated in the middle of the time step (used for stiffnesses). If 1, MPS is evaluated at the end of the time step. (Used for updating). |
Definition at line 618 of file b3solidmat.C.
References c0, c1, giveHumidity(), giveHumidityIncrement(), giveInitMicroPrestress(), oofem::B3SolidMaterialStatus::giveMPS(), oofem::Material::giveStatus(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::isTheFirstStep(), OOFEM_ERROR, S, and oofem::RheoChainMaterial::timeFactor.
Referenced by computeFlowTermViscosity(), and giveRealStressVector().
|
protected |
Evaluation of the shrinkageStrainVector. Shrinkage is fully dependent on humidity rate in given GP.
Definition at line 483 of file b3solidmat.C.
References oofem::FloatArray::at(), giveHumidityIncrement(), oofem::GaussPoint::giveMaterialMode(), oofem::StructuralMaterial::giveReducedSymVectorForm(), kSh, oofem::FloatArray::resize(), and oofem::FloatArray::zero().
Referenced by giveShrinkageStrainVector().
|
protected |
Evaluation of the relative volume of the solidified material.
Definition at line 511 of file b3solidmat.C.
References oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), lambda0, q2, q3, oofem::RheoChainMaterial::relMatAge, and oofem::RheoChainMaterial::timeFactor.
Referenced by giveEigenStrainVector(), and giveEModulus().
|
protected |
Definition at line 432 of file b3solidmat.C.
References oofem::FloatArray::at(), E28, EpsSinf, oofem::GaussPoint::giveMaterialMode(), oofem::StructuralMaterial::giveReducedSymVectorForm(), oofem::TimeStep::giveTargetTime(), hum, ks, kt, oofem::RheoChainMaterial::relMatAge, oofem::FloatArray::resize(), t0, oofem::RheoChainMaterial::timeFactor, vs, and oofem::FloatArray::zero().
Referenced by giveShrinkageStrainVector().
|
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::KelvinChainMaterial.
Definition at line 787 of file b3solidmat.C.
References oofem::FEMComponent::giveDomain(), and oofem::RheoChainMaterial::nUnits.
|
inlinevirtual |
Reimplemented from oofem::KelvinChainMaterial.
Definition at line 133 of file b3solidmat.h.
|
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::KelvinChainMaterial.
Definition at line 553 of file b3solidmat.C.
References oofem::FloatArray::add(), oofem::FloatArray::beProductOf(), computeFlowTermViscosity(), computeSolidifiedVolume(), oofem::KelvinChainMaterial::giveEigenStrainVector(), oofem::FloatMatrix::giveNumberOfRows(), oofem::Material::giveStatus(), oofem::TimeStep::giveTimeIncrement(), oofem::RheoChainMaterial::giveUnitComplianceMatrix(), oofem::RheoChainMaterialStatus::giveViscoelasticStressVector(), OOFEM_ERROR, oofem::FloatArray::resize(), oofem::RheoChainMaterial::timeFactor, oofem::FloatArray::times(), and oofem::FloatArray::zero().
|
protectedvirtual |
Evaluation of the incremental modulus.
Reimplemented from oofem::KelvinChainMaterial.
Definition at line 244 of file b3solidmat.C.
References computeFlowTermViscosity(), computeSolidifiedVolume(), EmoduliMode, EspringVal, oofem::KelvinChainMaterial::giveEModulus(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), q1, oofem::RheoChainMaterial::relMatAge, oofem::RheoChainMaterial::timeFactor, and updateEparModuli().
|
protected |
Computes relative humidity at given time step and GP.
Definition at line 721 of file b3solidmat.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(), inverse_sorption_isotherm(), and OOFEM_ERROR.
Referenced by computeMicroPrestress().
|
protected |
Computes relative humidity increment at given time step and GP.
Definition at line 752 of file b3solidmat.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(), inverse_sorption_isotherm(), and OOFEM_ERROR.
Referenced by computeMicroPrestress(), and computePointShrinkageStrainVector().
|
protected |
Computes initial value of the MicroPrestress.
Definition at line 712 of file b3solidmat.C.
Referenced by computeMicroPrestress().
|
inlinevirtual |
Implements oofem::FEMComponent.
Definition at line 134 of file b3solidmat.h.
References _IFT_B3SolidMaterial_Name, and oofem::MaterialStatus::initializeFrom().
|
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::KelvinChainMaterial.
Definition at line 797 of file b3solidmat.C.
References computeMicroPrestress(), oofem::KelvinChainMaterial::giveRealStressVector(), oofem::Material::giveStatus(), MicroPrestress, and oofem::B3SolidMaterialStatus::setMPS().
|
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::KelvinChainMaterial.
Definition at line 402 of file b3solidmat.C.
References B3_AverageShrinkage, B3_NoShrinkage, computePointShrinkageStrainVector(), computeTotalAverageShrinkageStrainVector(), oofem::GaussPoint::giveMaterialMode(), oofem::TimeStep::givePreviousStep(), oofem::StructuralMaterial::giveSizeOfVoigtSymVector(), oofem::TimeStep::isTheFirstStep(), OOFEM_ERROR, oofem::FloatArray::resize(), shMode, oofem::FloatArray::subtract(), and oofem::FloatArray::zero().
|
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::KelvinChainMaterial.
Definition at line 143 of file b3solidmat.h.
|
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. |
Reimplemented from oofem::KelvinChainMaterial.
Definition at line 49 of file b3solidmat.C.
References _IFT_B3Material_a, _IFT_B3Material_ac, _IFT_B3Material_alpha1, _IFT_B3Material_alpha2, _IFT_B3Material_cc, _IFT_B3Material_EpsSinf, _IFT_B3Material_fc, _IFT_B3Material_hum, _IFT_B3Material_ks, _IFT_B3Material_kt, _IFT_B3Material_mode, _IFT_B3Material_ncoeff, _IFT_B3Material_q1, _IFT_B3Material_q2, _IFT_B3Material_q3, _IFT_B3Material_q4, _IFT_B3Material_q5, _IFT_B3Material_shmode, _IFT_B3Material_t0, _IFT_B3Material_vs, _IFT_B3Material_wc, _IFT_B3Material_wh, _IFT_B3SolidMaterial_c0, _IFT_B3SolidMaterial_c1, _IFT_B3SolidMaterial_emodulimode, _IFT_B3SolidMaterial_finalhumidity, _IFT_B3SolidMaterial_initialhumidity, _IFT_B3SolidMaterial_ksh, _IFT_B3SolidMaterial_lambda0, _IFT_B3SolidMaterial_microprestress, a, B3_AverageShrinkage, B3_PointShrinkage, c0, c1, E28, EmoduliMode, EpsSinf, hum, oofem::KelvinChainMaterial::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_BAD_FORMAT, oofem::IRRT_OK, ks, kSh, kt, lambda0, MicroPrestress, n, OOFEM_ERROR, OOFEM_WARNING, predictParametersFrom(), q1, q2, q3, q4, q5, shMode, t0, vs, w, and w_h.
|
protected |
Definition at line 593 of file b3solidmat.C.
Referenced by giveHumidity(), and giveHumidityIncrement().
|
protected |
Definition at line 185 of file b3solidmat.C.
References B3_AverageShrinkage, EpsSinf, kt, oofem::FEMComponent::number, OOFEM_LOG_DEBUG, q1, q2, q3, q4, q5, shMode, and w.
Referenced by initializeFrom().
|
protectedvirtual |
Update of partial moduli of individual chain units.
Reimplemented from oofem::RheoChainMaterial.
Definition at line 286 of file b3solidmat.C.
References oofem::RheoChainMaterial::EparVal, oofem::FloatArray::isEmpty(), and oofem::RheoChainMaterial::updateEparModuli().
Referenced by giveEModulus().
|
protected |
Constant (obtained from experiments) A [Pedersen, 1990].
Definition at line 103 of file b3solidmat.h.
Referenced by initializeFrom(), and inverse_sorption_isotherm().
|
protected |
Definition at line 99 of file b3solidmat.h.
|
protected |
MPS constant c0 [MPa^-1 * day^-1].
Definition at line 115 of file b3solidmat.h.
Referenced by computeFlowTermViscosity(), computeMicroPrestress(), giveInitMicroPrestress(), and initializeFrom().
|
protected |
MPS constant c1 (=C1*R*T/M)
Definition at line 116 of file b3solidmat.h.
Referenced by computeMicroPrestress(), and initializeFrom().
|
protected |
Definition at line 90 of file b3solidmat.h.
Referenced by computeTotalAverageShrinkageStrainVector(), and initializeFrom().
|
protected |
If 0, analysis of retardation spectrum is used for evaluation of Kelvin units moduli (default).
If 1, least-squares method is used for evaluation of Kelvin units moduli.
Definition at line 109 of file b3solidmat.h.
Referenced by computeCharCoefficients(), giveEModulus(), and initializeFrom().
|
protected |
Additional parameters for average cross section shrinkage.
Definition at line 97 of file b3solidmat.h.
Referenced by computeTotalAverageShrinkageStrainVector(), initializeFrom(), and predictParametersFrom().
|
protected |
Additional parameters for free shrinkage at material point.
Definition at line 99 of file b3solidmat.h.
|
protected |
elastic modulus of the aging spring (first member of Kelvin chain if retardation spectrum is used)
Definition at line 104 of file b3solidmat.h.
Referenced by computeCharCoefficients(), and giveEModulus().
|
protected |
Definition at line 97 of file b3solidmat.h.
Referenced by computeTotalAverageShrinkageStrainVector(), and initializeFrom().
|
protected |
Definition at line 97 of file b3solidmat.h.
Referenced by computeTotalAverageShrinkageStrainVector(), and initializeFrom().
|
protected |
MPS shrinkage parameter. Either this or inithum and finalhum must be given in input record.
Definition at line 118 of file b3solidmat.h.
Referenced by computePointShrinkageStrainVector(), and initializeFrom().
|
protected |
Definition at line 97 of file b3solidmat.h.
Referenced by computeTotalAverageShrinkageStrainVector(), initializeFrom(), and predictParametersFrom().
|
protected |
constant equal to one day in time units of analysis (eg. 86400 if the analysis runs in seconds)
Definition at line 93 of file b3solidmat.h.
Referenced by computeCharCoefficients(), computeCharTimes(), computeCreepFunction(), computeSolidifiedVolume(), and initializeFrom().
|
protected |
If 1, computation exploiting Microprestress solidification theory is done.
Default value is 0 = without external fields it can be used for basic creep.
Definition at line 114 of file b3solidmat.h.
Referenced by computeFlowTermViscosity(), giveRealStressVector(), and initializeFrom().
|
protected |
Constant-exponent (obtained from experiments) n [Pedersen, 1990].
Definition at line 102 of file b3solidmat.h.
Referenced by computeCreepFunction(), initializeFrom(), and inverse_sorption_isotherm().
|
protected |
Definition at line 90 of file b3solidmat.h.
Referenced by giveEModulus(), initializeFrom(), and predictParametersFrom().
|
protected |
Definition at line 90 of file b3solidmat.h.
Referenced by computeCharCoefficients(), computeCreepFunction(), computeSolidifiedVolume(), initializeFrom(), and predictParametersFrom().
|
protected |
Definition at line 90 of file b3solidmat.h.
Referenced by computeSolidifiedVolume(), initializeFrom(), and predictParametersFrom().
|
protected |
Definition at line 90 of file b3solidmat.h.
Referenced by computeFlowTermViscosity(), initializeFrom(), and predictParametersFrom().
|
protected |
Definition at line 90 of file b3solidmat.h.
Referenced by initializeFrom(), and predictParametersFrom().
|
protected |
Definition at line 99 of file b3solidmat.h.
|
protected |
Definition at line 99 of file b3solidmat.h.
|
protected |
Referenced by giveShrinkageStrainVector(), initializeFrom(), and predictParametersFrom().
|
protected |
Definition at line 89 of file b3solidmat.h.
Referenced by computeTotalAverageShrinkageStrainVector(), and initializeFrom().
|
protected |
MPS tS0 - necessary for the initial value of microprestress (age when the load is applied)
Definition at line 117 of file b3solidmat.h.
Referenced by giveInitMicroPrestress().
|
protected |
Definition at line 97 of file b3solidmat.h.
Referenced by computeTotalAverageShrinkageStrainVector(), and initializeFrom().
|
protected |
Definition at line 90 of file b3solidmat.h.
Referenced by initializeFrom(), and predictParametersFrom().
|
protected |
Constant water content (obtained from experiments) w_h [Pedersen, 1990].
Definition at line 101 of file b3solidmat.h.
Referenced by initializeFrom(), and inverse_sorption_isotherm().