#include <b3mat.h>

Inheritance diagram for oofem::B3Material:

Collaboration diagram for oofem::B3Material:

Protected Types
enum	b3ShModeType { B3_NoShrinkage , B3_AverageShrinkage , B3_PointShrinkage }

Protected Attributes
double	t0 = 0.
	Age when drying begins (in days).
double	w = 0.
double	E28 = 0.
double	q1 = 0.
double	q2 = 0.
double	q3 = 0.
double	q4 = 0.
double	q5 = 0.
enum oofem::B3Material::b3ShModeType	shMode = B3_NoShrinkage
Additional parameters for average cross section shrinkage.
double	EpsSinf = 0.
double	kt = 0.
double	ks = 0.
double	vs = 0.
double	hum = 0.
Additional parameters for free shrinkage at material point.
double	es0 = 0.
double	r = 0.
double	rprime = 0.
double	at = 0.
Protected Attributes inherited from oofem::RheoChainMaterial
double	talpha = 0.
	thermal dilatation coeff.
int	nUnits = 0
	Number of (Maxwell or Kelvin) units in the rheologic chain.
double	relMatAge = 0.
	Physical age of the material at castingTime.
bool	lattice = false
double	nu = 0.
	Poisson's ratio (assumed to be constant, unaffected by creep).
double	alphaOne = 0.
	Parameters for the lattice model.
double	alphaTwo = 0.
double	EparValTime = -1.
	Time for which the partial moduli of individual units have been evaluated.
double	begOfTimeOfInterest = 0.
	Time from which the model should give a good approximation. Optional field. Default value is 0.1 [day].
double	endOfTimeOfInterest = 0.
	Time (age???) up to which the model should give a good approximation.
StructuralMaterial *	linearElasticMaterial = nullptr
	Associated linearElasticMaterial, with E = 1.
FloatArray	EparVal
	Partial moduli of individual units.
FloatArray	charTimes
	Characteristic times of individual units (relaxation or retardation times).
FloatArray	discreteTimeScale
	Times at which the errors are evaluated if the least-square method is used.
double	timeFactor = 0.
Protected Attributes inherited from oofem::StructuralMaterial
double	referenceTemperature = 0.
	Reference temperature (temperature, when material has been built into structure).
MatResponseMode	SCStiffMode = TangentStiffness
	stifness mode used in stress control
double	SCRelTol = 1.e-3
	relative tolerance for stress control
double	SCAbsTol = 1.e-12
	absolute stress tolerance for stress control
int	SCMaxiter = 100000
	maximum iterations for stress-control
Protected Attributes inherited from oofem::Material
Dictionary	propertyDictionary
double	castingTime
int	preCastingTimeMat
	Material existing before casting time - optional parameter, zero by default.
Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number.
Domain *	domain
	Link to domain object, useful for communicating with other FEM components.

Additional parameters for absorption isotherm (used to compute relative humidity from water content)
double	w_h = 0.
	Constant water content (obtained from experiments) w_h [Pedersen, 1990].
double	n = 0.
	Constant-exponent (obtained from experiments) n [Pedersen, 1990].
double	a = 0.
	Constant (obtained from experiments) A [Pedersen, 1990].
	B3Material (int n, Domain *d)
void	giveShrinkageStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode) const override
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
void	initializeFrom (InputRecord &ir) override
double	computeCreepFunction (double t, double t_prime, GaussPoint gp, TimeStep tStep) const override
	Evaluation of the creep compliance function at time t when loading is acting from time t_prime.
bool	hasIncrementalShrinkageFormulation () const override
virtual void	computeTotalAverageShrinkageStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep) const
virtual void	computeShrinkageStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode) const
	Free shrinkage at material point, requires staggered analysis.
void	predictParametersFrom (double, double, double, double, double, double, double)
double	inverse_sorption_isotherm (double w) const

Additional Inherited Members
Public Member Functions inherited from oofem::MaxwellChainMaterial
	MaxwellChainMaterial (int n, Domain *d)
void	giveRealStressVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep) const override
void	computeHiddenVars (GaussPoint gp, TimeStep tStep) const
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
void	giveShrinkageStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode) const override
void	giveEigenStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode) const override
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
Public Member Functions inherited from oofem::RheoChainMaterial
	RheoChainMaterial (int n, Domain *d)
virtual	~RheoChainMaterial ()
void	giveRealStressVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep) const override
FloatArrayF< 6 >	giveRealStressVector_3d (const FloatArrayF< 6 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
FloatArrayF< 4 >	giveRealStressVector_PlaneStrain (const FloatArrayF< 4 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector_3d.
FloatArrayF< 3 >	giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 1 >	giveRealStressVector_1d (const FloatArrayF< 1 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 2 >	giveRealStressVector_2dBeamLayer (const FloatArrayF< 2 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 5 >	giveRealStressVector_PlateLayer (const FloatArrayF< 5 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 6 >	giveThermalDilatationVector (GaussPoint gp, TimeStep tStep) const override
bool	hasMaterialModeCapability (MaterialMode mode) const override
bool	hasCastingTimeSupport () const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
void	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp) override
void	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp) override
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 3, 3 >	givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 4, 4 >	givePlaneStrainStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 1, 1 >	give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
FloatArray	computeStressIndependentStrainVector (GaussPoint gp, TimeStep tStep, ValueModeType mode) const override
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
double	giveAlphaOne () const
double	giveAlphaTwo () const
double	givePoissonsRatio () const
	Returns Poisson's ratio.
bool	isActivated (TimeStep *tStep) const override
	Extended meaning: returns true if the material is cast (target time > casting time) or the precasing time mat is defined.
virtual double	giveEquivalentTime (GaussPoint gp, TimeStep tStep) const
	By default returns equivalent time in the middle of the time step.
Public Member Functions inherited from oofem::StructuralMaterial
	StructuralMaterial (int n, Domain *d)
bool	hasMaterialModeCapability (MaterialMode mode) const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
void	giveCharacteristicMatrix (FloatMatrix &answer, MatResponseMode type, GaussPoint gp, TimeStep tStep) const override
	Returns characteristic matrix of the receiver.
void	giveCharacteristicVector (FloatArray &answer, FloatArray &flux, MatResponseMode type, GaussPoint gp, TimeStep tStep) const override
	Returns characteristic vector of the receiver.
virtual void	giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep) const
virtual FloatArray	giveRealStressVector_StressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint gp, TimeStep tStep) const
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArray	giveRealStressVector_ShellStressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint gp, TimeStep tStep) const
virtual FloatArrayF< 2 >	giveRealStressVector_Warping (const FloatArrayF< 2 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_Fiber (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_2dPlateSubSoil (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation is not provided.
virtual FloatArrayF< 6 >	giveRealStressVector_3dBeamSubSoil (const FloatArrayF< 6 > &reducedE, GaussPoint gp, TimeStep tStep) const
virtual FloatArrayF< 9 >	giveFirstPKStressVector_3d (const FloatArrayF< 9 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
virtual FloatArrayF< 5 >	giveFirstPKStressVector_PlaneStrain (const FloatArrayF< 5 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_3d.
virtual FloatArray	giveFirstPKStressVector_StressControl (const FloatArray &reducedvF, const IntArray &FControl, GaussPoint gp, TimeStep tStep) const
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArrayF< 4 >	giveFirstPKStressVector_PlaneStress (const FloatArrayF< 4 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
virtual FloatArrayF< 1 >	giveFirstPKStressVector_1d (const FloatArrayF< 1 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
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	giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
double	giveReferenceTemperature ()
FloatArrayF< 6 >	computeStressIndependentStrainVector_3d (GaussPoint gp, TimeStep tStep, ValueModeType mode) const
virtual void	giveStiffnessMatrix_dPdF (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 9, 9 >	give3dMaterialStiffnessMatrix_dPdF (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const
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) const
int	setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type) override
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
virtual FloatMatrixF< 4, 4 >	givePlaneStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	givePlaneStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 5, 5 >	givePlaneStrainStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	givePlaneStrainStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 1, 1 >	give1dStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	give1dStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 2, 2 >	give2dBeamLayerStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 5, 5 >	givePlateLayerStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 3, 3 >	giveFiberStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 3, 3 >	give2dPlateSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 6, 6 >	give3dBeamSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
Public Member Functions inherited from oofem::Material
	Material (int n, Domain *d)
virtual	~Material ()=default
	Destructor.
virtual bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode) const
virtual double	giveCharacteristicValue (MatResponseMode type, GaussPoint gp, TimeStep tStep) const
	Returns characteristic value of the receiver.
virtual double	give (int aProperty, GaussPoint *gp) const
virtual bool	hasProperty (int aProperty, GaussPoint *gp) const
virtual void	modifyProperty (int aProperty, double value, GaussPoint *gp)
double	giveCastingTime () const
void	printYourself () override
	Prints receiver state on stdout. Useful for debugging.
int	checkConsistency () override
virtual void	restoreConsistency (GaussPoint *gp)
virtual int	initMaterial (Element *element)
virtual MaterialStatus *	giveStatus (GaussPoint *gp) const
virtual int	packUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
virtual int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
virtual int	estimatePackSize (DataStream &buff, GaussPoint *ip)
virtual double	predictRelativeComputationalCost (GaussPoint *gp)
virtual double	predictRelativeRedistributionCost (GaussPoint *gp)
virtual void	initTempStatus (GaussPoint *gp) const
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
virtual	~FEMComponent ()=default
	Virtual destructor.
Domain *	giveDomain () const
virtual void	setDomain (Domain *d)
int	giveNumber () const
void	setNumber (int num)
virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
virtual void	initializeFrom (InputRecord &ir, int priority)
virtual void	initializeFinish ()
virtual void	postInitialize ()
	Performs post initialization steps. Called after all components are created and initialized.
virtual void	printOutputAt (FILE file, TimeStep tStep)
virtual Interface *	giveInterface (InterfaceType t)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).
Static Public Member Functions inherited from oofem::StructuralMaterial
static int	giveSymVI (int ind1, int ind2)
static int	giveVI (int ind1, int ind2)
static FloatMatrixF< 9, 9 >	convert_dSdE_2_dPdF_3D (const FloatMatrixF< 6, 6 > &dSdE, const FloatArrayF< 6 > &S, const FloatArrayF< 9 > &F)
static FloatMatrixF< 5, 5 >	convert_dSdE_2_dPdF_PlaneStrain (const FloatMatrixF< 4, 4 > &dSdE, const FloatArrayF< 4 > &S, const FloatArrayF< 5 > &F)
static FloatMatrixF< 4, 4 >	convert_dSdE_2_dPdF_PlaneStress (const FloatMatrixF< 3, 3 > &dSdE, const FloatArrayF< 3 > &S, const FloatArrayF< 4 > &F)
static FloatMatrixF< 1, 1 >	convert_dSdE_2_dPdF_1D (const FloatMatrixF< 1, 1 > &dSdE, const FloatArrayF< 1 > &S, const FloatArrayF< 1 > &F)
static void	computePrincipalValues (FloatArray &answer, const FloatArray &s, stressStrainPrincMode mode)
	Common functions for convenience.
static FloatArrayF< 3 >	computePrincipalValues (const FloatMatrixF< 3, 3 > &s)
static FloatArrayF< 3 >	computePrincipalValues (double I1, double I2, double I3)
static void	computePrincipalValDir (FloatArray &answer, FloatMatrix &dir, const FloatArray &s, stressStrainPrincMode mode)
static std::pair< FloatArrayF< 3 >, FloatMatrixF< 3, 3 > >	computePrincipalValDir (const FloatMatrixF< 3, 3 > &s)
static FloatArrayF< 6 >	computeDeviator (const FloatArrayF< 6 > &s)
static std::pair< FloatArrayF< 6 >, double >	computeDeviatoricVolumetricSplit (const FloatArrayF< 6 > &s)
static FloatArrayF< 6 >	computeDeviatoricVolumetricSum (const FloatArrayF< 6 > &dev, double mean)
static FloatArrayF< 6 >	applyDeviatoricElasticCompliance (const FloatArrayF< 6 > &stress, double EModulus, double nu)
static FloatArrayF< 6 >	applyDeviatoricElasticCompliance (const FloatArrayF< 6 > &stress, double GModulus)
static FloatArrayF< 6 >	applyDeviatoricElasticStiffness (const FloatArrayF< 6 > &strain, double EModulus, double nu)
static FloatArrayF< 6 >	applyDeviatoricElasticStiffness (const FloatArrayF< 6 > &strain, double GModulus)
static FloatArrayF< 6 >	applyElasticStiffness (const FloatArrayF< 6 > &strain, double EModulus, double nu)
static FloatArrayF< 6 >	applyElasticCompliance (const FloatArrayF< 6 > &stress, double EModulus, double nu)
static double	computeStressNorm (const FloatArrayF< 6 > &stress)
static double	computeFirstInvariant (const FloatArrayF< 6 > &s)
static double	computeSecondStressInvariant (const FloatArrayF< 6 > &s)
static double	computeThirdStressInvariant (const FloatArrayF< 6 > &s)
static double	computeFirstCoordinate (const FloatArrayF< 6 > &s)
static double	computeSecondCoordinate (const FloatArrayF< 6 > &s)
static double	computeThirdCoordinate (const FloatArrayF< 6 > &s)
static int	giveVoigtVectorMask (IntArray &answer, MaterialMode mmode)
static int	giveVoigtSymVectorMask (IntArray &answer, MaterialMode mmode)
static void	giveInvertedVoigtVectorMask (IntArray &answer, MaterialMode mmode)
static int	giveSizeOfVoigtVector (MaterialMode mmode)
static int	giveSizeOfVoigtSymVector (MaterialMode mmode)
static void	giveFullVectorForm (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
	Converts the reduced symmetric Voigt vector (2nd order tensor) to full form.
static void	giveFullVectorFormF (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
	Converts the reduced deformation gradient Voigt vector (2nd order tensor).
static void	giveFullSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the reduced unsymmetric Voigt vector (2nd order tensor) to full form.
static void	giveReducedVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the full symmetric Voigt vector (2nd order tensor) to reduced form.
static void	giveReducedSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the full unsymmetric Voigt vector (2nd order tensor) to reduced form.
static void	giveFullSymMatrixForm (FloatMatrix &answer, const FloatMatrix &red, MaterialMode matMode)
	Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form.
static void	giveReducedMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
	Converts the full symmetric Voigt matrix (4th order tensor) to reduced form.
static void	giveReducedSymMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
	Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form.
static FloatArrayF< 6 >	transformStrainVectorTo (const FloatMatrixF< 3, 3 > &base, const FloatArrayF< 6 > &strain, bool transpose=false)
static FloatArrayF< 6 >	transformStressVectorTo (const FloatMatrixF< 3, 3 > &base, const FloatArrayF< 6 > &stress, bool transpose=false)
static double	computeVonMisesStress (const FloatArray &currentStress)
static double	computeVonMisesStress_3D (const FloatArrayF< 6 > &stress)
static double	computeVonMisesStress_PlaneStress (const FloatArrayF< 3 > &stress)
static FloatMatrixF< 6, 6 >	giveStrainVectorTranformationMtrx (const FloatMatrixF< 3, 3 > &base, bool transpose=false)
static FloatMatrixF< 3, 3 >	give2DStrainVectorTranformationMtrx (const FloatMatrixF< 2, 2 > &base, bool transpose=false)
static FloatMatrixF< 6, 6 >	giveStressVectorTranformationMtrx (const FloatMatrixF< 3, 3 > &base, bool transpose=false)
static FloatMatrixF< 3, 3 >	givePlaneStressVectorTranformationMtrx (const FloatMatrixF< 2, 2 > &base, bool transpose=false)
static void	sortPrincDirAndValCloseTo (FloatArray &pVal, FloatMatrix &pDir, const FloatMatrix &toPDir)
Static Public Attributes inherited from oofem::StructuralMaterial
static std::array< std::array< int, 3 >, 3 >	vIindex
	Voigt index map.
static std::array< std::array< int, 3 >, 3 >	svIndex
	Symmetric Voigt index map.
Protected Member Functions inherited from oofem::MaxwellChainMaterial
bool	hasIncrementalShrinkageFormulation () const override
FloatArray	computeCharCoefficients (double tPrime, GaussPoint gp, TimeStep tStep) const override
double	giveEModulus (GaussPoint gp, TimeStep tStep) const override
	Evaluation of the incremental modulus.
Protected Member Functions inherited from oofem::RheoChainMaterial
const FloatArray &	giveDiscreteTimes () const
void	computeDiscreteRelaxationFunction (FloatArray &answer, const FloatArray &tSteps, double t0, double tr, GaussPoint gp, TimeStep tStep) const
void	giveUnitComplianceMatrix (FloatMatrix &answer, GaussPoint gp, TimeStep tStep) const
	Evaluation of elastic compliance matrix for unit Young's modulus.
void	giveUnitStiffnessMatrix (FloatMatrix &answer, GaussPoint gp, TimeStep tStep) const
	Evaluation of elastic stiffness matrix for unit Young's modulus.
virtual void	updateEparModuli (double tPrime, GaussPoint gp, TimeStep tStep) const
	Update of partial moduli of individual chain units.
double	giveEparModulus (int iChain) const
	Access to partial modulus of a given unit.
virtual void	computeCharTimes ()
	Evaluation of characteristic times.
double	giveCharTime (int) const
	Access to the characteristic time of a given unit.
virtual double	giveCharTimeExponent (int i) const
	Exponent to be used with the char time of a given unit, usually = 1.0.
StructuralMaterial *	giveLinearElasticMaterial ()
	Access to the underlying linear elastic material with unit Young's modulus.
double	giveEndOfTimeOfInterest ()
	Access to the time up to which the response should be accurate.
void	computeTrueStressIndependentStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode) const
Static Protected Member Functions inherited from oofem::RheoChainMaterial
static FloatArray	generateLogTimeScale (double from, double to, int nsteps)

Detailed Description

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

Definition at line 75 of file b3mat.h.

Member Enumeration Documentation

◆ b3ShModeType

enum oofem::B3Material::b3ShModeType

protected

Enumerator
B3_NoShrinkage
B3_AverageShrinkage
B3_PointShrinkage

Definition at line 82 of file b3mat.h.

Constructor & Destructor Documentation

◆ B3Material()

oofem::B3Material::B3Material	(	int	n,
		Domain *	d )

inline

Definition at line 98 of file b3mat.h.

References oofem::MaxwellChainMaterial::MaxwellChainMaterial(), and n.

Member Function Documentation

◆ computeCreepFunction()

double oofem::B3Material::computeCreepFunction	(	double	t,
		double	t_prime,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Evaluation of the creep compliance function at time t when loading is acting from time t_prime.

Implements oofem::RheoChainMaterial.

Definition at line 185 of file b3mat.C.

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

Referenced by computeShrinkageStrainVector().

◆ computeShrinkageStrainVector()

void oofem::B3Material::computeShrinkageStrainVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep,
		ValueModeType	mode ) const

protectedvirtual

Free shrinkage at material point, requires staggered analysis.

Definition at line 316 of file b3mat.C.

References oofem::FloatArray::add(), at, oofem::FloatArray::at(), oofem::Material::castingTime, computeCreepFunction(), oofem::Element::computeGlobalCoordinates(), oofem::FEMComponent::domain, es0, oofem::GaussPoint::giveElement(), oofem::FieldManager::giveField(), oofem::GaussPoint::giveMaterialMode(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::RheoChainMaterialStatus::giveShrinkageStrainVector(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), oofem::TimeStep::giveTargetTime(), oofem::RheoChainMaterialStatus::giveViscoelasticStressVector(), inverse_sorption_isotherm(), OOFEM_ERROR, r, oofem::RheoChainMaterial::relMatAge, oofem::FloatArray::resize(), rprime, oofem::sgn(), t0, oofem::RheoChainMaterial::timeFactor, and oofem::FloatArray::zero().

Referenced by giveShrinkageStrainVector().

◆ computeTotalAverageShrinkageStrainVector()

void oofem::B3Material::computeTotalAverageShrinkageStrainVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

protectedvirtual

Definition at line 264 of file b3mat.C.

References oofem::FloatArray::at(), oofem::Material::castingTime, E28, EpsSinf, oofem::GaussPoint::giveMaterialMode(), oofem::TimeStep::giveTargetTime(), hum, ks, kt, oofem::RheoChainMaterial::relMatAge, t0, oofem::RheoChainMaterial::timeFactor, and vs.

Referenced by giveShrinkageStrainVector().

◆ giveClassName()

const char * oofem::B3Material::giveClassName ( ) const

inlineoverridevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 102 of file b3mat.h.

◆ giveInputRecordName()

const char * oofem::B3Material::giveInputRecordName ( ) const

inlineoverridevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 103 of file b3mat.h.

References _IFT_B3Material_Name.

◆ giveShrinkageStrainVector()

void oofem::B3Material::giveShrinkageStrainVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep,
		ValueModeType	mode ) const

overridevirtual

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

Parameters

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::RheoChainMaterial.

Definition at line 235 of file b3mat.C.

References B3_AverageShrinkage, B3_NoShrinkage, computeShrinkageStrainVector(), computeTotalAverageShrinkageStrainVector(), oofem::GaussPoint::giveMaterialMode(), oofem::TimeStep::givePreviousStep(), oofem::TimeStep::isTheFirstStep(), OOFEM_ERROR, oofem::FloatArray::resize(), shMode, oofem::FloatArray::subtract(), and oofem::FloatArray::zero().

◆ hasIncrementalShrinkageFormulation()

bool oofem::B3Material::hasIncrementalShrinkageFormulation ( ) const

inlineoverrideprotectedvirtual

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::RheoChainMaterial.

Definition at line 109 of file b3mat.h.

◆ initializeFrom()

void oofem::B3Material::initializeFrom ( InputRecord & ir )

overridevirtual

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. Note that initializeFrom may be called mutiple times.

See also: IR_GIVE_FIELD; IR_GIVE_OPTIONAL_FIELD

Parameters

ir	Input record to initialize from.
priority	Priority of the input record. This is used to determine the order of initialization

Reimplemented from oofem::FEMComponent.

Definition at line 46 of file b3mat.C.

References _IFT_B3Material_a, _IFT_B3Material_ac, _IFT_B3Material_alpha1, _IFT_B3Material_alpha2, _IFT_B3Material_at, _IFT_B3Material_cc, _IFT_B3Material_EpsSinf, _IFT_B3Material_es0, _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_r, _IFT_B3Material_rprime, _IFT_B3Material_shmode, _IFT_B3Material_t0, _IFT_B3Material_vs, _IFT_B3Material_wc, _IFT_B3Material_wh, a, at, B3_AverageShrinkage, B3_PointShrinkage, E28, EpsSinf, es0, hum, IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, ks, kt, n, predictParametersFrom(), q1, q2, q3, q4, q5, r, rprime, shMode, t0, vs, w, and w_h.

◆ inverse_sorption_isotherm()

double oofem::B3Material::inverse_sorption_isotherm ( double w ) const

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