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:

Public Member Functions
	B3SolidMaterial (int n, Domain *d)
virtual	~B3SolidMaterial ()
virtual void	updateYourself (GaussPoint gp, TimeStep tStep)
	Update of MatStatus to the newly reached (equilibrium) state.
virtual void	giveShrinkageStrainVector (FloatArray &answer, MatResponseForm form, 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 classType	giveClassID () const
	Returns classType id of receiver.
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 MaterialStatus *	CreateStatus (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, MatResponseForm form, GaussPoint gp, TimeStep tStep)
void	computePointShrinkageStrainVectorMPS (FloatArray &answer, MatResponseForm form, GaussPoint gp, TimeStep tStep)
	Evaluation of the shrinkageStrainVector. Shrinkage is fully dependent on humidity rate in given GP.
void	computeShrinkageStrainVector (FloatArray &answer, MatResponseForm form, GaussPoint gp, TimeStep tStep, ValueModeType mode)
void	predictParametersFrom (double, double, double, double, double, double, double)
double	computeNonAgingCreepFunction (GaussPoint *gp, double loadDuration)
	Evaluation of the compliance function of the non-aging solidifying constituent.
double	computeSolidifiedVolume (GaussPoint gp, 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 (GaussPoint *gp, double tStep, double ofAge)
	Evaluation of the creep compliance function.
double	inverse_sorption_isotherm (double w)
virtual void	computeCharCoefficients (FloatArray &answer, GaussPoint *gp, double)
	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, MatResponseForm form, 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 (=C1RT/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 74 of file b3solidmat.h.

Member Enumeration Documentation

enum oofem::B3SolidMaterial::b3ShModeType [protected]

Enumerator:

B3_NoShrinkage
B3_AverageShrinkage
B3_PointShrinkage
B3_PointShrinkageMPS

Definition at line 79 of file b3solidmat.h.

Constructor & Destructor Documentation

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

Definition at line 107 of file b3solidmat.h.

References B3_NoShrinkage, and shMode.

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

Definition at line 108 of file b3solidmat.h.

Member Function Documentation

void oofem::B3SolidMaterial::computeCharCoefficients	(	FloatArray &	answer,
		GaussPoint *	gp,
		double	atTime
	)		`[protected, virtual]`

Evaluation of characteristic moduli of the non-aging Kelvin chain.

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 342 of file b3solidmat.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::RheoChainMaterial::charTimes, computeNonAgingCreepFunction(), EmoduliMode, EspringVal, oofem::RheoChainMaterial::giveCharTime(), oofem::RheoChainMaterial::giveDiscreteTimes(), oofem::FloatArray::giveSize(), oofem::RheoChainMaterial::nUnits, q2, r, oofem::FloatArray::resize(), and oofem::FloatMatrix::solveForRhs().

void oofem::B3SolidMaterial::computeCharTimes ( ) [protected, virtual]

Evaluation of characteristic times.

Reimplemented from oofem::RheoChainMaterial.

Definition at line 299 of file b3solidmat.C.

References oofem::FloatArray::at(), oofem::RheoChainMaterial::begOfTimeOfInterest, oofem::RheoChainMaterial::charTimes, oofem::RheoChainMaterial::endOfTimeOfInterest, oofem::RheoChainMaterial::giveEndOfTimeOfInterest(), oofem::RheoChainMaterial::nUnits, and oofem::FloatArray::resize().

double oofem::B3SolidMaterial::computeCreepFunction	(	GaussPoint *	gp,
		double	ofAge,
		double	atTime
	)		`[protected, virtual]`

Evaluation of the creep compliance function.

Implements oofem::KelvinChainMaterial.

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]`

Evaluation of the flow term viscosity.

Definition at line 641 of file b3solidmat.C.

References _error, c0, computeMicroPrestress(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), MicroPrestress, q4, oofem::RheoChainMaterial::relMatAge, S, and oofem::RheoChainMaterial::timeFactor.

Referenced by giveEigenStrainVector(), and giveEModulus().

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

Computes microprestress at given time step and GP.

Parameters:

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 854 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 updateYourself().

double oofem::B3SolidMaterial::computeNonAgingCreepFunction	(	GaussPoint *	gp,
		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().

void oofem::B3SolidMaterial::computePointShrinkageStrainVectorMPS	(	FloatArray &	answer,
		MatResponseForm	form,
		GaussPoint *	gp,
		TimeStep *	tStep
	)		`[protected]`

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

Definition at line 586 of file b3solidmat.C.

References oofem::FloatArray::at(), oofem::FullForm, oofem::GaussPoint::giveCrossSection(), giveHumidityIncrement(), oofem::GaussPoint::giveMaterialMode(), oofem::StructuralMaterial::giveReducedCharacteristicVector(), kSh, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

Referenced by giveShrinkageStrainVector().

void oofem::B3SolidMaterial::computeShrinkageStrainVector	(	FloatArray &	answer,
		MatResponseForm	form,
		GaussPoint *	gp,
		TimeStep *	tStep,
		ValueModeType	mode
	)		`[protected]`

Definition at line 708 of file b3solidmat.C.

References _error, _error2, oofem::FloatArray::add(), at, oofem::FloatArray::at(), computeCreepFunction(), oofem::Element::computeGlobalCoordinates(), oofem::FEMComponent::domain, es0, oofem::Field::evaluateAt(), oofem::FullForm, oofem::EngngModel::giveContext(), oofem::GaussPoint::giveCoordinates(), oofem::GaussPoint::giveCrossSection(), oofem::GaussPoint::giveElement(), oofem::Domain::giveEngngModel(), oofem::FieldManager::giveField(), oofem::EngngModelContext::giveFieldManager(), oofem::StructuralMaterial::giveFullCharacteristicVector(), oofem::GaussPoint::giveMaterialMode(), oofem::StructuralMaterial::giveReducedCharacteristicVector(), oofem::RheoChainMaterialStatus::giveShrinkageStrainVector(), oofem::FloatArray::giveSize(), oofem::StructuralMaterial::giveSizeOfReducedStressStrainVector(), oofem::Material::giveStatus(), oofem::StructuralMaterialStatus::giveStressVector(), oofem::TimeStep::giveTargetTime(), inverse_sorption_isotherm(), r, oofem::RheoChainMaterial::relMatAge, oofem::FloatArray::resize(), rprime, oofem::sgn(), t0, oofem::RheoChainMaterial::timeFactor, and oofem::FloatArray::zero().

Referenced by giveShrinkageStrainVector().

double oofem::B3SolidMaterial::computeSolidifiedVolume	(	GaussPoint *	gp,
		TimeStep *	tStep
	)		`[protected]`

Evaluation of the relative volume of the solidified material.

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

void oofem::B3SolidMaterial::computeTotalAverageShrinkageStrainVector	(	FloatArray &	answer,
		MatResponseForm	form,
		GaussPoint *	gp,
		TimeStep *	tStep
	)		`[protected]`

MaterialStatus * oofem::B3SolidMaterial::CreateStatus ( GaussPoint * gp ) const [virtual]

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

Parameters:

gp	Integration point where newly created status will be stored.

Returns:: Reference to new status.

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 1023 of file b3solidmat.C.

virtual classType oofem::B3SolidMaterial::giveClassID ( ) const [inline, virtual]

Returns classType id of receiver.

Intended for run time type checking. Every derived class have to overload this method.

See also:: classType.

Returns:: Class type of receiver.

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 116 of file b3solidmat.h.

References oofem::B3SolidMaterialClass.

virtual const char* oofem::B3SolidMaterial::giveClassName ( ) const [inline, virtual]

Returns:: Class name of the receiver.

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 115 of file b3solidmat.h.

void oofem::B3SolidMaterial::giveEigenStrainVector	(	FloatArray &	answer,
		MatResponseForm	form,
		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:

answer	Computed strains.
form	Material response form.
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 662 of file b3solidmat.C.

References _error, oofem::FloatArray::add(), oofem::FloatArray::beProductOf(), computeFlowTermViscosity(), computeSolidifiedVolume(), oofem::GaussPoint::giveCrossSection(), oofem::StructuralMaterial::giveFullCharacteristicVector(), oofem::Matrix::giveNumberOfRows(), oofem::Material::giveStatus(), oofem::StructuralMaterialStatus::giveStressVector(), oofem::TimeStep::giveTimeIncrement(), oofem::RheoChainMaterial::giveUnitComplianceMatrix(), oofem::ReducedForm, oofem::FloatArray::resize(), oofem::RheoChainMaterial::timeFactor, and oofem::FloatArray::times().

double oofem::B3SolidMaterial::giveEModulus	(	GaussPoint *	gp,
		TimeStep *	atTime
	)		`[protected, virtual]`

Evaluation of the incremental modulus.

Reimplemented from oofem::KelvinChainMaterial.

Definition at line 261 of file b3solidmat.C.

References computeFlowTermViscosity(), computeSolidifiedVolume(), EmoduliMode, oofem::RheoChainMaterial::EparVal, EspringVal, oofem::FloatArray::giveSize(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), q1, oofem::RheoChainMaterial::relMatAge, oofem::RheoChainMaterial::timeFactor, and oofem::RheoChainMaterial::updateEparModuli().

double oofem::B3SolidMaterial::giveHumidity	(	GaussPoint *	gp,
		TimeStep *	tStep
	)		`[protected]`

Computes relative humidity at given time step and GP.

Definition at line 957 of file b3solidmat.C.

References _error, _error2, oofem::FloatArray::at(), oofem::Element::computeGlobalCoordinates(), oofem::FEMComponent::domain, oofem::Field::evaluateAt(), oofem::EngngModel::giveContext(), oofem::GaussPoint::giveCoordinates(), oofem::GaussPoint::giveElement(), oofem::Domain::giveEngngModel(), oofem::FieldManager::giveField(), oofem::EngngModelContext::giveFieldManager(), and inverse_sorption_isotherm().

Referenced by computeMicroPrestress().

double oofem::B3SolidMaterial::giveHumidityIncrement	(	GaussPoint *	gp,
		TimeStep *	tStep
	)		`[protected]`

Computes relative humidity increment at given time step and GP.

Definition at line 988 of file b3solidmat.C.

References _error, _error2, oofem::FloatArray::at(), oofem::Element::computeGlobalCoordinates(), oofem::FEMComponent::domain, oofem::Field::evaluateAt(), oofem::EngngModel::giveContext(), oofem::GaussPoint::giveCoordinates(), oofem::GaussPoint::giveElement(), oofem::Domain::giveEngngModel(), oofem::FieldManager::giveField(), oofem::EngngModelContext::giveFieldManager(), and inverse_sorption_isotherm().

Referenced by computeMicroPrestress(), and computePointShrinkageStrainVectorMPS().

double oofem::B3SolidMaterial::giveInitMicroPrestress ( void ) [protected]

Computes initial value of the MicroPrestress.

Definition at line 948 of file b3solidmat.C.

References c0, and tS0.

Referenced by computeMicroPrestress().

void oofem::B3SolidMaterial::giveShrinkageStrainVector	(	FloatArray &	answer,
		MatResponseForm	form,
		GaussPoint *	gp,
		TimeStep *	tStep,
		ValueModeType	mode
	)		`[virtual]`

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

Parameters:

answer	Returned strain vector.
form	Material response form.
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 489 of file b3solidmat.C.

References _error, B3_AverageShrinkage, B3_NoShrinkage, B3_PointShrinkageMPS, computePointShrinkageStrainVectorMPS(), computeShrinkageStrainVector(), computeTotalAverageShrinkageStrainVector(), oofem::FullForm, oofem::GaussPoint::giveMaterialMode(), oofem::TimeStep::givePreviousStep(), oofem::StructuralMaterial::giveSizeOfReducedStressStrainVector(), 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:

answer	Vector of thermal dilatation coefficients.
gp	Integration point.
tStep	Time step (most models are able to respond only when atTime is current time step).

Reimplemented from oofem::RheoChainMaterial.

Definition at line 244 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 124 of file b3solidmat.h.

IRResultType oofem::B3SolidMaterial::initializeFrom ( InputRecord * ir ) [virtual]

Initializes receiver according to object description stored in input record.

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