This class implements the B3 model for concrete creep and shrinkage. More...

#include <b3mat.h>

Inheritance diagram for oofem::B3Material:

Collaboration diagram for oofem::B3Material:

Public Member Functions
	B3Material (int n, Domain *d)
virtual	~B3Material ()
virtual void	giveShrinkageStrainVector (FloatArray &answer, MatResponseForm form, GaussPoint gp, TimeStep atTime, 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 gp, TimeStep tStep)
	Returns a vector of coefficients of thermal dilatation in direction of each material principal (local) axis.
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.
virtual void	computeTotalAverageShrinkageStrainVector (FloatArray &answer, MatResponseForm form, GaussPoint gp, TimeStep atTime)
virtual void	computeShrinkageStrainVector (FloatArray &answer, MatResponseForm form, GaussPoint gp, TimeStep atTime, ValueModeType mode)
	Free shrinkage at material point, requires staggered analysis.
void	predictParametersFrom (double, double, double, double, double, double, double)
virtual double	computeCreepFunction (GaussPoint *gp, double atTime, double ofAge)
	Evaluation of the creep compliance function.
double	inverse_sorption_isotherm (double w)
	Function calculates relative humidity from water content (inverse relation form sorption isotherm).
Protected Attributes
double	t0
	Age when drying begins (in days)
double	w
double	E28
double	q1
double	q2
double	q3
double	q4
double	q5
enum oofem::B3Material::b3ShModeType	shMode
Additional parameters for average cross section shrinkage.
double	EpsSinf
double	kt
double	ks
double	vs
double	hum
Additional parameters for free shrinkage at material point.
double	es0
double	r
double	rprime
double	at
Additional parameters for absorption isotherm (used to compute relative humidity from water content)
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.

Detailed Description

This class implements the B3 model for concrete creep and shrinkage.

Definition at line 44 of file b3mat.h.

Member Enumeration Documentation

enum oofem::B3Material::b3ShModeType [protected]

Enumerator:

B3_NoShrinkage
B3_AverageShrinkage
B3_PointShrinkage

Definition at line 51 of file b3mat.h.

Constructor & Destructor Documentation

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

Definition at line 68 of file b3mat.h.

References B3_NoShrinkage, and shMode.

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

Definition at line 69 of file b3mat.h.

Member Function Documentation

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

Evaluation of the creep compliance function.

Implements oofem::MaxwellChainMaterial.

Definition at line 205 of file b3mat.C.

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

Referenced by computeShrinkageStrainVector().

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

Free shrinkage at material point, requires staggered analysis.

Definition at line 350 of file b3mat.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().

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

virtual classType oofem::B3Material::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::MaxwellChainMaterial.

Definition at line 75 of file b3mat.h.

References oofem::B3MaterialClass.

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

Returns:: Class name of the receiver.

Reimplemented from oofem::MaxwellChainMaterial.

Definition at line 74 of file b3mat.h.

void oofem::B3Material::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::MaxwellChainMaterial.

Definition at line 258 of file b3mat.C.

References _error, B3_AverageShrinkage, B3_NoShrinkage, 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::B3Material::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 125 of file b3mat.C.

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

virtual int oofem::B3Material::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::MaxwellChainMaterial.

Definition at line 81 of file b3mat.h.

IRResultType oofem::B3Material::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.

See also:: IR_GIVE_FIELD; IR_GIVE_OPTIONAL_FIELD

Parameters:

ir	Input record to initialize from.

Returns:: IRResultType

Reimplemented from oofem::MaxwellChainMaterial.

Definition at line 43 of file b3mat.C.

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

Function calculates relative humidity from water content (inverse relation form sorption isotherm).

Relative humidity (phi) is from range 0.2 - 0.98 !!! Sorption isotherm by C. R. Pedersen (1990), Combined heat and moisture transfer in building constructions, PhD-thesis, Technical University of Denmark, Lingby.

Parameters: