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, 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 double	computeCreepFunction (double t, double t_prime, GaussPoint gp, TimeStep tStep)
	Evaluation of the creep compliance function at time t when loading is acting from time t_prime. More...

Public Member Functions inherited from oofem::MaxwellChainMaterial
	MaxwellChainMaterial (int n, Domain *d)

virtual	~MaxwellChainMaterial ()

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

void	computeHiddenVars (GaussPoint gp, TimeStep tStep)

virtual int	hasNonLinearBehaviour ()
	Returns nonzero if receiver is non linear. 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...

virtual MaterialStatus *	CreateStatus (GaussPoint *gp) const
	Creates new copy of associated status and inserts it into given integration point. More...

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

virtual void	computeTotalAverageShrinkageStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep)

virtual void	computeShrinkageStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode)
	Free shrinkage at material point, requires staggered analysis. More...

void	predictParametersFrom (double, double, double, double, double, double, double)

double	inverse_sorption_isotherm (double w)
	Function calculates relative humidity from water content (inverse relation form sorption isotherm). More...

Protected Member Functions inherited from oofem::MaxwellChainMaterial
virtual void	computeCharCoefficients (FloatArray &answer, double tPrime, GaussPoint gp, TimeStep tStep)
	This function computes the moduli of individual Maxwell units such that the corresponding Dirichlet series gives the best approximation of the actual relaxation function. More...

virtual double	giveEModulus (GaussPoint gp, TimeStep tStep)
	Evaluation of the incremental modulus. More...

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

virtual void	updateEparModuli (double tPrime, GaussPoint gp, TimeStep tStep)
	Update of partial moduli of individual chain units. More...

double	giveEparModulus (int iChain)
	Access to partial modulus of a given unit. More...

virtual void	computeCharTimes ()
	Evaluation of characteristic times. 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
	Age when drying begins (in days) More...

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]. More...

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

double	a
	Constant (obtained from experiments) A [Pedersen, 1990]. 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...

Detailed Description

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

Definition at line 75 of file b3mat.h.

Member Enumeration Documentation

enum oofem::B3Material::b3ShModeType

protected

Enumerator
B3_NoShrinkage
B3_AverageShrinkage
B3_PointShrinkage

Definition at line 82 of file b3mat.h.

Constructor & Destructor Documentation

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

inline

Definition at line 98 of file b3mat.h.

References B3_NoShrinkage.

virtual oofem::B3Material::~B3Material ( )

inlinevirtual

Definition at line 101 of file b3mat.h.

References giveShrinkageStrainVector().

Member Function Documentation

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

virtual

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

Implements oofem::RheoChainMaterial.

Definition at line 190 of file b3mat.C.

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

Referenced by computeShrinkageStrainVector(), and giveInputRecordName().

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

protectedvirtual

Free shrinkage at material point, requires staggered analysis.

Definition at line 319 of file b3mat.C.

References oofem::FloatArray::add(), at, oofem::FloatArray::at(), computeCreepFunction(), oofem::Element::computeGlobalCoordinates(), oofem::FEMComponent::domain, es0, oofem::EngngModel::giveContext(), oofem::GaussPoint::giveElement(), oofem::Domain::giveEngngModel(), oofem::FieldManager::giveField(), oofem::EngngModelContext::giveFieldManager(), oofem::StructuralMaterial::giveFullSymVectorForm(), oofem::GaussPoint::giveMaterialMode(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::StructuralMaterial::giveReducedSymVectorForm(), oofem::RheoChainMaterialStatus::giveShrinkageStrainVector(), oofem::FloatArray::giveSize(), oofem::StructuralMaterial::giveSizeOfVoigtSymVector(), 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(), and hasIncrementalShrinkageFormulation().

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

protectedvirtual

Definition at line 271 of file b3mat.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(), and hasIncrementalShrinkageFormulation().

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

inlinevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::MaxwellChainMaterial.

Definition at line 105 of file b3mat.h.

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

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 106 of file b3mat.h.

References _IFT_B3Material_Name, computeCreepFunction(), and initializeFrom().

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

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

Definition at line 241 of file b3mat.C.

References B3_AverageShrinkage, B3_NoShrinkage, computeShrinkageStrainVector(), 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().

Referenced by ~B3Material().

virtual int oofem::B3Material::hasIncrementalShrinkageFormulation ( )

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

Definition at line 112 of file b3mat.h.

References computeShrinkageStrainVector(), computeTotalAverageShrinkageStrainVector(), inverse_sorption_isotherm(), and predictParametersFrom().

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 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, oofem::MaxwellChainMaterial::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_OK, ks, kt, n, predictParametersFrom(), q1, q2, q3, q4, q5, r, rprime, shMode, t0, vs, w, and w_h.

Referenced by giveInputRecordName().

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