OOFEM  2.3
Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes
oofem::B3SolidMaterial Class Reference

This class implements the B3 model for concrete creep and shrinkage based on the solidification theory. More...

#include <b3solidmat.h>

+ Inheritance diagram for oofem::B3SolidMaterial:
+ Collaboration diagram for oofem::B3SolidMaterial:

List of all members.

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.
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.
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.
virtual void giveThermalDilatationVector (FloatArray &answer, GaussPoint *, TimeStep *)
 Returns a vector of coefficients of thermal dilatation in direction of each material principal (local) axis.
virtual MaterialStatusCreateStatus (GaussPoint *gp) const
 Creates new copy of associated status and inserts it into given integration point.

Protected Types

enum  b3ShModeType { B3_NoShrinkage, B3_AverageShrinkage, B3_PointShrinkage, B3_PointShrinkageMPS }

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.
void computeTotalAverageShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)
void computePointShrinkageStrainVectorMPS (FloatArray &answer, GaussPoint *gp, TimeStep *tStep)
 Evaluation of the shrinkageStrainVector. Shrinkage is fully dependent on humidity rate in given GP.
void computeShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode)
void predictParametersFrom (double, double, double, double, double, double, double)
double computeNonAgingCreepFunction (double loadDuration)
 Evaluation of the compliance function of the non-aging solidifying constituent.
double computeSolidifiedVolume (TimeStep *tStep)
 Evaluation of the relative volume of the solidified material.
double computeFlowTermViscosity (GaussPoint *gp, TimeStep *tStep)
 Evaluation of the flow term viscosity.
virtual double computeCreepFunction (double tStep, double ofAge)
 Evaluation of the creep compliance function.
double inverse_sorption_isotherm (double w)
virtual void computeCharCoefficients (FloatArray &answer, double tStep)
 Evaluation of characteristic moduli of the non-aging Kelvin chain.
virtual void computeCharTimes ()
 Evaluation of characteristic times.
virtual double giveEModulus (GaussPoint *gp, TimeStep *tStep)
 Evaluation of the incremental modulus.
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).
double computeMicroPrestress (GaussPoint *gp, TimeStep *tStep, int option)
 Computes microprestress at given time step and GP.
double giveInitMicroPrestress (void)
 Computes initial value of the MicroPrestress.
double giveHumidity (GaussPoint *gp, TimeStep *tStep)
 Computes relative humidity at given time step and GP.
double giveHumidityIncrement (GaussPoint *gp, TimeStep *tStep)
 Computes relative humidity increment at given time step and GP.

Protected Attributes

double t0
double w
double E28
double q1
double q2
double q3
double q4
double q5
enum
oofem::B3SolidMaterial::b3ShModeType 
shMode
double EpsSinf
 Additional parameters for average cross section shrinkage.
double kt
double ks
double vs
double hum
double es0
 Additional parameters for free shrinkage at material point.
double r
double rprime
double at
double w_h
 Constant water content (obtained from experiments) w_h [Pedersen, 1990].
double n
 Constant-exponent (obtained from experiments) n [Pedersen, 1990].
double a
 Constant (obtained from experiments) A [Pedersen, 1990].
double talpha
 Thermal dilatation coeff.
double EspringVal
 elastic modulus of the aging spring (first member of Kelvin chain if retardation spectrum is used)
int EmoduliMode
 If 0, analysis of retardation spectrum is used for evaluation of Kelvin units moduli (default).
int MicroPrestress
 If 1, computation exploiting Microprestress solidification theory is done.
double c0
 MPS constant c0 [MPa^-1 * day^-1].
double c1
 MPS constant c1 (=C1*R*T/M)
double tS0
 MPS tS0 - necessary for the initial value of microprestress (age when the load is applied)
double kSh
 MPS shrinkage parameter. Either this or inithum and finalhum must be given in input record.

Detailed Description

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 75 of file b3solidmat.h.


Member Enumeration Documentation

Enumerator:
B3_NoShrinkage 
B3_AverageShrinkage 
B3_PointShrinkage 
B3_PointShrinkageMPS 

Definition at line 80 of file b3solidmat.h.


Constructor & Destructor Documentation

oofem::B3SolidMaterial::B3SolidMaterial ( int  n,
Domain d 
) [inline]

Definition at line 108 of file b3solidmat.h.

References B3_NoShrinkage, and shMode.

virtual oofem::B3SolidMaterial::~B3SolidMaterial ( ) [inline, virtual]

Definition at line 111 of file b3solidmat.h.


Member Function Documentation

void oofem::B3SolidMaterial::computeCharCoefficients ( FloatArray answer,
double  tStep 
) [protected, virtual]
void oofem::B3SolidMaterial::computeCharTimes ( ) [protected, virtual]
double oofem::B3SolidMaterial::computeCreepFunction ( double  ofAge,
double  tStep 
) [protected, virtual]

Evaluation of the creep compliance function.

Implements oofem::RheoChainMaterial.

Definition at line 440 of file b3solidmat.C.

References B3_AverageShrinkage, hum, ks, kt, n, q1, q2, q3, q4, q5, r, shMode, t0, and vs.

Referenced by computeShrinkageStrainVector().

double oofem::B3SolidMaterial::computeFlowTermViscosity ( GaussPoint gp,
TimeStep tStep 
) [protected]
double oofem::B3SolidMaterial::computeMicroPrestress ( GaussPoint gp,
TimeStep tStep,
int  option 
) [protected]

Computes microprestress at given time step and GP.

Parameters:
gpGauss point to compute at.
tStepTime step to compute for.
optionIf 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 813 of file b3solidmat.C.

References _error, c0, c1, giveHumidity(), giveHumidityIncrement(), giveInitMicroPrestress(), oofem::B3SolidMaterialStatus::giveMPS(), oofem::Material::giveStatus(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::isTheFirstStep(), S, and oofem::RheoChainMaterial::timeFactor.

Referenced by computeFlowTermViscosity(), and giveRealStressVector().

double oofem::B3SolidMaterial::computeNonAgingCreepFunction ( double  loadDuration) [protected]

Evaluation of the compliance function of the non-aging solidifying constituent.

Definition at line 424 of file b3solidmat.C.

References n, and q2.

Referenced by computeCharCoefficients().

Evaluation of the shrinkageStrainVector. Shrinkage is fully dependent on humidity rate in given GP.

Definition at line 572 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().

void oofem::B3SolidMaterial::computeShrinkageStrainVector ( FloatArray answer,
GaussPoint gp,
TimeStep tStep,
ValueModeType  mode 
) [protected]

Evaluation of the relative volume of the solidified material.

Definition at line 600 of file b3solidmat.C.

References oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), q2, q3, oofem::RheoChainMaterial::relMatAge, and oofem::RheoChainMaterial::timeFactor.

Referenced by giveEigenStrainVector(), and giveEModulus().

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

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

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 982 of file b3solidmat.C.

virtual const char* oofem::B3SolidMaterial::giveClassName ( ) const [inline, virtual]
Returns:
Class name of the receiver.

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 119 of file b3solidmat.h.

void oofem::B3SolidMaterial::giveEigenStrainVector ( FloatArray answer,
GaussPoint gp,
TimeStep tStep,
ValueModeType  mode 
) [protected, virtual]

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

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

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 642 of file b3solidmat.C.

References _error, oofem::FloatArray::add(), oofem::FloatArray::beProductOf(), computeFlowTermViscosity(), computeSolidifiedVolume(), oofem::FloatMatrix::giveNumberOfRows(), oofem::Material::giveStatus(), oofem::StructuralMaterialStatus::giveStressVector(), oofem::TimeStep::giveTimeIncrement(), oofem::RheoChainMaterial::giveUnitComplianceMatrix(), oofem::FloatArray::resize(), oofem::RheoChainMaterial::timeFactor, and oofem::FloatArray::times().

double oofem::B3SolidMaterial::giveEModulus ( GaussPoint gp,
TimeStep tStep 
) [protected, virtual]
double oofem::B3SolidMaterial::giveHumidity ( GaussPoint gp,
TimeStep tStep 
) [protected]
double oofem::B3SolidMaterial::giveHumidityIncrement ( GaussPoint gp,
TimeStep tStep 
) [protected]
double oofem::B3SolidMaterial::giveInitMicroPrestress ( void  ) [protected]

Computes initial value of the MicroPrestress.

Definition at line 907 of file b3solidmat.C.

References c0, and tS0.

Referenced by computeMicroPrestress().

virtual const char* oofem::B3SolidMaterial::giveInputRecordName ( ) const [inline, virtual]
Returns:
Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 120 of file b3solidmat.h.

References _IFT_B3SolidMaterial_Name.

void oofem::B3SolidMaterial::giveRealStressVector ( FloatArray answer,
GaussPoint gp,
const FloatArray reducedStrain,
TimeStep tStep 
) [virtual]

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

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

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

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 992 of file b3solidmat.C.

References computeMicroPrestress(), oofem::Material::giveStatus(), MicroPrestress, and oofem::B3SolidMaterialStatus::setMPS().

void oofem::B3SolidMaterial::giveShrinkageStrainVector ( FloatArray answer,
GaussPoint gp,
TimeStep tStep,
ValueModeType  mode 
) [virtual]

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

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

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 489 of file b3solidmat.C.

References _error, B3_AverageShrinkage, B3_NoShrinkage, B3_PointShrinkageMPS, computePointShrinkageStrainVectorMPS(), computeShrinkageStrainVector(), computeTotalAverageShrinkageStrainVector(), oofem::GaussPoint::giveMaterialMode(), oofem::TimeStep::givePreviousStep(), oofem::StructuralMaterial::giveSizeOfVoigtSymVector(), oofem::TimeStep::isTheFirstStep(), oofem::FloatArray::resize(), shMode, oofem::FloatArray::subtract(), and oofem::FloatArray::zero().

void oofem::B3SolidMaterial::giveThermalDilatationVector ( FloatArray answer,
GaussPoint gp,
TimeStep tStep 
) [virtual]

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

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

Reimplemented from oofem::RheoChainMaterial.

Definition at line 247 of file b3solidmat.C.

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

virtual int oofem::B3SolidMaterial::hasIncrementalShrinkageFormulation ( ) [inline, protected, virtual]

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 128 of file b3solidmat.h.

Initializes receiver according to object description stored in input record.

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

See also:
IR_GIVE_FIELD
IR_GIVE_OPTIONAL_FIELD
Parameters:
irInput record to initialize from.
Returns:
IRResultType

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 48 of file b3solidmat.C.

References _error, _IFT_B3Material_a, _IFT_B3Material_ac, _IFT_B3Material_alpha1, _IFT_B3Material_alpha2, _IFT_B3Material_at, _IFT_B3Material_c0, _IFT_B3Material_c1, _IFT_B3Material_cc, _IFT_B3Material_emodulimode, _IFT_B3Material_EpsSinf, _IFT_B3Material_es0, _IFT_B3Material_fc, _IFT_B3Material_finalhumidity, _IFT_B3Material_hum, _IFT_B3Material_initialhumidity, _IFT_B3Material_ks, _IFT_B3Material_ksh, _IFT_B3Material_kt, _IFT_B3Material_microprestress, _IFT_B3Material_mode, _IFT_B3Material_ncoeff, _IFT_B3Material_q1, _IFT_B3Material_q2, _IFT_B3Material_q3, _IFT_B3Material_q4, _IFT_B3Material_q5, _IFT_B3Material_r, _IFT_B3Material_rprime, _IFT_B3Material_shmode, _IFT_B3Material_t0, _IFT_B3Material_talpha, _IFT_B3Material_ts0, _IFT_B3Material_vs, _IFT_B3Material_wc, _IFT_B3Material_wh, a, at, B3_AverageShrinkage, B3_PointShrinkage, B3_PointShrinkageMPS, c0, c1, E28, EmoduliMode, EpsSinf, es0, hum, IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_OK, ks, kSh, kt, MicroPrestress, n, predictParametersFrom(), q1, q2, q3, q4, q5, r, rprime, shMode, t0, talpha, tS0, vs, w, and w_h.

double oofem::B3SolidMaterial::inverse_sorption_isotherm ( double  w) [protected]

Definition at line 788 of file b3solidmat.C.

References a, n, and w_h.

Referenced by computeShrinkageStrainVector(), giveHumidity(), and giveHumidityIncrement().

void oofem::B3SolidMaterial::predictParametersFrom ( double  fc,
double  c,
double  wc,
double  ac,
double  t0,
double  alpha1,
double  alpha2 
) [protected]

Definition at line 189 of file b3solidmat.C.

References B3_AverageShrinkage, EpsSinf, kt, oofem::FEMComponent::number, oofem::OOFEM_LOG_DEBUG(), q1, q2, q3, q4, q5, shMode, and w.

Referenced by initializeFrom().


Member Data Documentation

double oofem::B3SolidMaterial::a [protected]

Constant (obtained from experiments) A [Pedersen, 1990].

Definition at line 88 of file b3solidmat.h.

Referenced by initializeFrom(), and inverse_sorption_isotherm().

double oofem::B3SolidMaterial::at [protected]

Definition at line 84 of file b3solidmat.h.

Referenced by computeShrinkageStrainVector(), and initializeFrom().

double oofem::B3SolidMaterial::c0 [protected]

MPS constant c0 [MPa^-1 * day^-1].

Definition at line 101 of file b3solidmat.h.

Referenced by computeFlowTermViscosity(), computeMicroPrestress(), giveInitMicroPrestress(), and initializeFrom().

double oofem::B3SolidMaterial::c1 [protected]

MPS constant c1 (=C1*R*T/M)

Definition at line 102 of file b3solidmat.h.

Referenced by computeMicroPrestress(), and initializeFrom().

double oofem::B3SolidMaterial::E28 [protected]

Definition at line 79 of file b3solidmat.h.

Referenced by computeTotalAverageShrinkageStrainVector(), and initializeFrom().

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 95 of file b3solidmat.h.

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

double oofem::B3SolidMaterial::EpsSinf [protected]

Additional parameters for average cross section shrinkage.

Definition at line 82 of file b3solidmat.h.

Referenced by computeTotalAverageShrinkageStrainVector(), initializeFrom(), and predictParametersFrom().

double oofem::B3SolidMaterial::es0 [protected]

Additional parameters for free shrinkage at material point.

Definition at line 84 of file b3solidmat.h.

Referenced by computeShrinkageStrainVector(), and initializeFrom().

elastic modulus of the aging spring (first member of Kelvin chain if retardation spectrum is used)

Definition at line 90 of file b3solidmat.h.

Referenced by computeCharCoefficients(), and giveEModulus().

double oofem::B3SolidMaterial::hum [protected]
double oofem::B3SolidMaterial::ks [protected]
double oofem::B3SolidMaterial::kSh [protected]

MPS shrinkage parameter. Either this or inithum and finalhum must be given in input record.

Definition at line 104 of file b3solidmat.h.

Referenced by computePointShrinkageStrainVectorMPS(), and initializeFrom().

double oofem::B3SolidMaterial::kt [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 100 of file b3solidmat.h.

Referenced by computeFlowTermViscosity(), giveRealStressVector(), and initializeFrom().

double oofem::B3SolidMaterial::n [protected]

Constant-exponent (obtained from experiments) n [Pedersen, 1990].

Definition at line 87 of file b3solidmat.h.

Referenced by computeCreepFunction(), computeNonAgingCreepFunction(), initializeFrom(), and inverse_sorption_isotherm().

double oofem::B3SolidMaterial::q1 [protected]
double oofem::B3SolidMaterial::q2 [protected]
double oofem::B3SolidMaterial::q3 [protected]
double oofem::B3SolidMaterial::q4 [protected]
double oofem::B3SolidMaterial::q5 [protected]

Definition at line 79 of file b3solidmat.h.

Referenced by computeCreepFunction(), initializeFrom(), and predictParametersFrom().

double oofem::B3SolidMaterial::r [protected]
double oofem::B3SolidMaterial::rprime [protected]

Definition at line 84 of file b3solidmat.h.

Referenced by computeShrinkageStrainVector(), and initializeFrom().

double oofem::B3SolidMaterial::t0 [protected]
double oofem::B3SolidMaterial::talpha [protected]

Thermal dilatation coeff.

Definition at line 89 of file b3solidmat.h.

Referenced by giveThermalDilatationVector(), and initializeFrom().

double oofem::B3SolidMaterial::tS0 [protected]

MPS tS0 - necessary for the initial value of microprestress (age when the load is applied)

Definition at line 103 of file b3solidmat.h.

Referenced by giveInitMicroPrestress(), and initializeFrom().

double oofem::B3SolidMaterial::vs [protected]
double oofem::B3SolidMaterial::w [protected]

Definition at line 79 of file b3solidmat.h.

Referenced by initializeFrom(), and predictParametersFrom().

double oofem::B3SolidMaterial::w_h [protected]

Constant water content (obtained from experiments) w_h [Pedersen, 1990].

Definition at line 86 of file b3solidmat.h.

Referenced by initializeFrom(), and inverse_sorption_isotherm().


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

This page is part of the OOFEM documentation. Copyright (c) 2011 Borek Patzak
Project e-mail: info@oofem.org
Generated at Tue Feb 18 2014 10:53:50 for OOFEM by doxygen 1.7.6.1 written by Dimitri van Heesch, © 1997-2011