oofem::ConcreteFCMViscoElastic Class Reference

#include <concretefcmviscoelastic.h>

Inheritance diagram for oofem::ConcreteFCMViscoElastic:

Collaboration diagram for oofem::ConcreteFCMViscoElastic:

Public Member Functions
	ConcreteFCMViscoElastic (int n, Domain *d)
void	initializeFrom (InputRecord &ir) override
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
double	give (int aProperty, GaussPoint *gp) const override
void	giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) const override
FloatArray	computeStressIndependentStrainVector (GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override
int	giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
MaterialStatus *	giveStatus (GaussPoint *gp) const override
Public Member Functions inherited from oofem::ConcreteFCM
	ConcreteFCM (int n, Domain *d)
Public Member Functions inherited from oofem::FCMMaterial
	FCMMaterial (int n, Domain *d)
bool	hasMaterialModeCapability (MaterialMode mode) const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
double	give (int aProperty, GaussPoint *gp) const 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
void	giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) const override
virtual void	initializeCrack (GaussPoint *gp, TimeStep *tStep, FloatMatrix &base, int nCrack) const
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.
int	giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
virtual double	computeNormalCrackOpening (GaussPoint *gp, int i) const
	uses temporary cracking strain and characteristic length to obtain the crack opening
virtual double	computeMaxNormalCrackOpening (GaussPoint *gp, TimeStep *tStep, int i) const
	uses maximum equilibrated cracking strain and characteristic length to obtain the maximum reached crack opening
virtual double	computeShearSlipOnCrack (GaussPoint *gp, TimeStep *tStep, int i) const
	computes total shear slip on a given crack plane (i = 1, 2, 3); the slip is computed from the temporary cracking strain
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
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)
virtual FloatArrayF< 6 >	giveThermalDilatationVector (GaussPoint *gp, TimeStep *tStep) const
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 >	give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
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 bool	hasProperty (int aProperty, GaussPoint *gp) const
virtual void	modifyProperty (int aProperty, double value, GaussPoint *gp)
double	giveCastingTime () const
virtual bool	isActivated (TimeStep *tStep) const
virtual bool	hasCastingTimeSupport () const
void	printYourself () override
	Prints receiver state on stdout. Useful for debugging.
virtual void	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
virtual void	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
int	checkConsistency () override
virtual void	restoreConsistency (GaussPoint *gp)
virtual int	initMaterial (Element *element)
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).
Public Member Functions inherited from oofem::RandomMaterialExtensionInterface
	RandomMaterialExtensionInterface ()
	Constructor.
virtual	~RandomMaterialExtensionInterface ()
	Destructor.
void	initializeFrom (InputRecord &ir)
void	giveInputRecord (DynamicInputRecord &ir)
bool	give (int key, GaussPoint *gp, double &value) const
Public Member Functions inherited from oofem::Interface
	Interface ()
	Constructor.
virtual	~Interface ()

Protected Member Functions
double	giveTensileStrength (GaussPoint *gp, TimeStep *tStep) const override
	comutes tensile strength
double	giveFractureEnergy (GaussPoint *gp, TimeStep *tStep) const override
double	computeOverallElasticStiffness (GaussPoint *gp, TimeStep *tStep) const override
	returns overall Young's modulus
double	computeOverallElasticShearModulus (GaussPoint *gp, TimeStep *tStep) const override
	returns overall shear modulus
int	checkConsistency (void) override
virtual double	giveEquivalentTime (GaussPoint *gp, TimeStep *tStep) const
	returns equivalent time (used to compute time-dependent ft and gf)
Protected Member Functions inherited from oofem::ConcreteFCM
double	giveCrackingModulus (MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep, int i) const override
	returns stiffness in the normal direction of the i-th crack
double	giveCrackingModulusInTension (MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep, int i) const override
	returns stiffness in the normal direction of the i-th crack
double	computeEffectiveShearModulus (GaussPoint *gp, TimeStep *tStep, int i) const override
	returns Geff which is necessary in the global stiffness matrix
double	computeD2ModulusForCrack (GaussPoint *gp, TimeStep *tStep, int icrack) const override
	shear modulus for a given crack plane (1, 2, 3)
double	computeNumerD2ModulusForCrack (GaussPoint *gp, TimeStep *tStep, int icrack) const override
	shear modulus for numerical purpose (stiffness matrix) for a given crack plane (1, 2, 3)
double	giveNormalCrackingStress (GaussPoint *gp, TimeStep *tStep, double eps_cr, int i) const override
	computes normal stress associated with i-th crack direction
double	maxShearStress (GaussPoint *gp, TimeStep *tStep, int i) const override
	computes the maximum value of the shear stress; if the shear stress exceeds this value, it is cropped
virtual double	computeResidualTensileStrength (GaussPoint *gp, TimeStep *tStep) const
	based on the maximum crack opening evaluates the residual strength
void	checkSnapBack (GaussPoint *gp, TimeStep *tStep, int crack) const override
	checks possible snap-back
Protected Member Functions inherited from oofem::FCMMaterial
virtual void	updateCrackStatus (GaussPoint *gp) const
	updates crack statuses
virtual double	giveCharacteristicElementLength (GaussPoint *gp, const FloatArray &crackPlaneNormal) const
	returns characteristic element length in given direction
virtual double	computeTotalD2Modulus (GaussPoint *gp, TimeStep *tStep, int i) const
	shear modulus for a given shear direction (4, 5, 6)
virtual double	computeNumerD2Modulus (GaussPoint *gp, TimeStep *tStep, int i) const
	shear modulus in a STIFFNESS MATRIX for a given shear direction (4, 5, 6)
virtual bool	isIntactForShear (GaussPoint *gp, int i) const
	returns true for closed or no cracks associated to given shear direction (i = 4, 5, 6)
virtual bool	isIntact (GaussPoint *gp, int icrack) const
	returns true for closed or no crack (i = 1, 2, 3)
virtual bool	isThisShearComponent (GaussPoint *gp, int component) const
	returns true if current component is associated with shear
virtual bool	checkStrengthCriterion (FloatMatrix &newBase, const FloatArray &globalStress, GaussPoint *gp, TimeStep *tStep, int nCrack) const
	checks if the globalStress does not exceed strength in the direction of newBase for n-th crack
virtual bool	isStrengthExceeded (const FloatMatrix &base, GaussPoint *gp, TimeStep *tStep, int iCrack, double trialStress) const
	compares trial stress with strength. Returns true if the strength is exceeded. Function oveloaded in the nonlocal approach for the fiber reinforced composites
virtual double	computeShearStiffnessRedistributionFactor (GaussPoint *gp, TimeStep *tStep, int ithCrackPlane, int jthCrackDirection) const
	function calculating ratio used to split shear slips on two crack planes
virtual double	giveCrackSpacing (void) const
	returns either user-provided value of crack spacing or a value computed from composition
virtual double	giveNumberOfCracksInDirection (GaussPoint *gp, int iCrack) const
	returns number of fictiotious parallel cracks in the direction of i-th crack
virtual double	giveNumberOfCracksForShearDirection (GaussPoint *gp, int i) const
	returns number of cracks for given shear direction (i = 4, 5, 6) which is treated as the maximum of the two associated normal directions
virtual void	giveMaterialStiffnessMatrix (FloatMatrix &answer, MatResponseMode, GaussPoint *gp, TimeStep *tStep) const
virtual void	giveTotalLocalCrackedStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep) const
	returns local stiffness matrix in the cracks' direction (total according to the material mode)
virtual void	giveNormalLocalCrackedStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep) const
	returns local stiffness matrix in the cracks' direction (only normal components)
virtual double	givePoissonsRatio () const
	returns Poisson's ratio
Protected Member Functions inherited from oofem::RandomMaterialExtensionInterface
void	_generateStatusVariables (GaussPoint *) const

Protected Attributes
int	viscoMat = 0
	number of the viscoelastic material
bool	fib = false
double	fib_s = 0.
double	fib_fcm28 = 0.
double	timeFactor = 0.
double	stiffnessFactor = 0.
Protected Attributes inherited from oofem::ConcreteFCM
double	Gf = 0.
	Fracture energy.
double	Ft = 0.
	Tensile strength.
double	beta = 0.
	shear retention factor
double	sf = 0.
	shear factor
double	sf_numer = 0.
	shear factor for numerical purpose
double	fc = 0.
	Collins' aggregate interlock: compressive strength in MPa.
double	ag = 0.
	Collins' aggregate interlock: aggregate diameter in appropriate units (same as FE mesh).
double	lengthScale = 0.
	Collins' aggregate interlock: 1 for meter, 1000 for analysis in mm.
FloatArray	soft_w
	user-defined softening (traction-COD)
FloatArray	soft_function_w
FloatArray	soft_eps
	user-defined softening (traction-strain)
FloatArray	soft_function_eps
FloatArray	beta_w
	user-defined shear retention factor (with respect to crack opening)
FloatArray	beta_function
double	H = 0.
	hardening modulus
double	eps_f = 0.
	strain at failure
SofteningType	softType = ST_NONE
ShearRetentionType	shearType = SHR_Unknown
ShearStrengthType	shearStrengthType = SHS_Unknown
Protected Attributes inherited from oofem::FCMMaterial
IsotropicLinearElasticMaterial	linearElasticMaterial
int	nAllowedCracks
	allowed number of cracks (user-defined)
ElementCharSizeMethod	ecsMethod
	Method used for evaluation of characteristic element size.
bool	multipleCrackShear
	if true = takes shear compliance of all cracks, false = only dominant crack contribution, default value is false
int	iterLimitGlobal = 20
double	crackSpacing
	value of crack spacing (allows to "have" more parallel cracks in one direction if the element size exceeds user-defined or computed crack spacing).
double	shearCoeffNumer
	minimum ratio of effective shear modulus vs. shear modulus - just for numerical purpose
double	normalCoeffNumer
	minimum ratio of effective normal stiffness vs. overall modulus - just for numerical purpose
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.
Protected Attributes inherited from oofem::RandomMaterialExtensionInterface
IntArray	randVariables
	Array of randomized variables (identified by a key).
IntArray	randomVariableGenerators
	Array of generators id's for corresponding randomized variables.

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 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 Types inherited from oofem::ConcreteFCM
enum	SofteningType {   ST_NONE , ST_Exponential , ST_Linear , ST_Hordijk ,   ST_UserDefinedCrack , ST_LinearHardeningStrain , ST_UserDefinedStrain , ST_Unknown }
	type of post-peak behavior in the normal direction to the crack plane More...
enum	ShearRetentionType {   SHR_NONE , SHR_Const_ShearRetFactor , SHR_Const_ShearFactorCoeff , SHR_UserDefined_ShearRetFactor ,   SHR_Unknown }
	type of reduction of the shear stiffness caused by cracking More...
enum	ShearStrengthType {   SHS_NONE , SHS_Const_Ft , SHS_Collins_Interlock , SHS_Residual_Ft ,   SHS_Unknown }
	defines the maximum value of shear stress More...

Detailed Description

add description

Definition at line 96 of file concretefcmviscoelastic.h.

Constructor & Destructor Documentation

ConcreteFCMViscoElastic()

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

Definition at line 45 of file concretefcmviscoelastic.C.

References oofem::ConcreteFCM::ConcreteFCM().

Member Function Documentation

checkConsistency()

int oofem::ConcreteFCMViscoElastic::checkConsistency ( void )

overrideprotectedvirtual

Allows programmer to test some internal data, before computation begins. For example, one may use this function, to ensure that element has material with required capabilities is assigned to element. This must be done after all mesh components are instanciated.

Returns

Nonzero if receiver is consistent.

Reimplemented from oofem::FEMComponent.

Definition at line 358 of file concretefcmviscoelastic.C.

References oofem::FEMComponent::checkConsistency(), oofem::FEMComponent::domain, oofem::RheoChainMaterial::givePoissonsRatio(), oofem::FCMMaterial::linearElasticMaterial, OOFEM_ERROR, tAlpha, and viscoMat.

computeOverallElasticShearModulus()

double oofem::ConcreteFCMViscoElastic::computeOverallElasticShearModulus ( GaussPoint * gp, TimeStep * tStep ) const

overrideprotectedvirtual

returns overall shear modulus

Reimplemented from oofem::FCMMaterial.

Definition at line 341 of file concretefcmviscoelastic.C.

References oofem::FEMComponent::domain, give(), oofem::RheoChainMaterial::giveEModulus(), oofem::GaussPoint::giveMaterialStatus(), oofem::ConcreteFCMViscoElasticStatus::giveSlaveGaussPointVisco(), oofem::FCMMaterial::linearElasticMaterial, and viscoMat.

computeOverallElasticStiffness()

double oofem::ConcreteFCMViscoElastic::computeOverallElasticStiffness ( GaussPoint * gp, TimeStep * tStep ) const

overrideprotectedvirtual

returns overall Young's modulus

Reimplemented from oofem::FCMMaterial.

Definition at line 328 of file concretefcmviscoelastic.C.

References oofem::FEMComponent::domain, oofem::RheoChainMaterial::giveEModulus(), oofem::GaussPoint::giveMaterialStatus(), oofem::ConcreteFCMViscoElasticStatus::giveSlaveGaussPointVisco(), and viscoMat.

computeStressIndependentStrainVector()

FloatArray oofem::ConcreteFCMViscoElastic::computeStressIndependentStrainVector ( GaussPoint * gp, TimeStep * tStep, ValueModeType mode ) const

overridevirtual

Computes reduced strain vector in given integration point, generated by internal processes in material, which are independent on loading in particular integration point. Default implementation takes into account temperature induced strains and eigenstrains.

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

Todo

Hack for loose gausspoints. We shouldn't ask for "gp->giveElement()". FIXME

Reimplemented from oofem::StructuralMaterial.

Definition at line 122 of file concretefcmviscoelastic.C.

References oofem::RheoChainMaterial::computeStressIndependentStrainVector(), oofem::FEMComponent::domain, oofem::GaussPoint::giveMaterialStatus(), oofem::ConcreteFCMViscoElasticStatus::giveSlaveGaussPointVisco(), oofem::Material::isActivated(), and viscoMat.

CreateStatus()

std::unique_ptr< MaterialStatus > oofem::ConcreteFCMViscoElastic::CreateStatus ( GaussPoint * gp ) const

inlineoverridevirtual

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

Definition at line 105 of file concretefcmviscoelastic.h.

Referenced by giveStatus().

give()

double oofem::ConcreteFCMViscoElastic::give ( int aProperty, GaussPoint * gp ) const

overridevirtual

Returns the value of material property 'aProperty'. Property must be identified by unique int id. Integration point also passed to allow for materials with spatially varying properties

Parameters

aProperty	ID of property requested.
gp	Integration point,

Returns

Property value.

Reimplemented from oofem::ConcreteFCM.

Definition at line 84 of file concretefcmviscoelastic.C.

References oofem::ConcreteFCM::give().

Referenced by computeOverallElasticShearModulus().

giveClassName()

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

inlineoverridevirtual

Returns

Class name of the receiver.

Reimplemented from oofem::ConcreteFCM.

Definition at line 102 of file concretefcmviscoelastic.h.

giveEquivalentTime()

double oofem::ConcreteFCMViscoElastic::giveEquivalentTime ( GaussPoint * gp, TimeStep * tStep ) const

protectedvirtual

returns equivalent time (used to compute time-dependent ft and gf)

Definition at line 379 of file concretefcmviscoelastic.C.

References oofem::FEMComponent::domain, oofem::RheoChainMaterial::giveEquivalentTime(), oofem::GaussPoint::giveMaterialStatus(), oofem::ConcreteFCMViscoElasticStatus::giveSlaveGaussPointVisco(), and viscoMat.

Referenced by giveTensileStrength().

giveFractureEnergy()

double oofem::ConcreteFCMViscoElastic::giveFractureEnergy ( GaussPoint * gp, TimeStep * tStep ) const

overrideprotectedvirtual

Reimplemented from oofem::ConcreteFCM.

Definition at line 277 of file concretefcmviscoelastic.C.

References fib, fib_fcm28, oofem::ConcreteFCM::Ft, oofem::ConcreteFCM::Gf, oofem::ConcreteFCM::giveFractureEnergy(), oofem::ConcreteFCMViscoElasticStatus::giveFractureEnergy(), oofem::GaussPoint::giveMaterialStatus(), oofem::FCMMaterialStatus::giveNumberOfTempCracks(), giveTensileStrength(), oofem::ConcreteFCMViscoElasticStatus::setFractureEnergy(), and stiffnessFactor.

giveInputRecordName()

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

inlineoverridevirtual

Returns

Input record name of the receiver.

Reimplemented from oofem::ConcreteFCM.

Definition at line 103 of file concretefcmviscoelastic.h.

References _IFT_ConcreteFCMViscoElastic_Name.

giveIPValue()

int oofem::ConcreteFCMViscoElastic::giveIPValue ( FloatArray & answer, GaussPoint * gp, InternalStateType type, TimeStep * tStep )

overridevirtual

Returns the integration point corresponding value in Reduced form.

Parameters

answer	Contain corresponding ip value, zero sized if not available.
gp	Integration point to which the value refers.
type	Determines the type of internal variable.
tStep	Determines the time step.

Returns

Nonzero if the assignment can be done, zero if this type of variable is not supported.

Reimplemented from oofem::ConcreteFCM.

Definition at line 139 of file concretefcmviscoelastic.C.

References oofem::FloatArray::at(), oofem::StructuralMaterial::computePrincipalValues(), oofem::FEMComponent::domain, oofem::StructuralMaterial::giveFullSymVectorForm(), oofem::ConcreteFCM::giveIPValue(), oofem::RheoChainMaterial::giveIPValue(), oofem::GaussPoint::giveMaterialMode(), oofem::GaussPoint::giveMaterialStatus(), oofem::FCMMaterialStatus::giveMaxCrackStrain(), oofem::ConcreteFCM::giveNormalCrackingStress(), oofem::FCMMaterialStatus::giveNumberOfCracks(), oofem::FCMMaterialStatus::giveNumberOfTempCracks(), oofem::ConcreteFCMViscoElasticStatus::giveSlaveGaussPointVisco(), oofem::StructuralMaterialStatus::giveStressVector(), giveTensileStrength(), oofem::ConcreteFCMViscoElasticStatus::giveTensileStrength(), oofem::max(), oofem::min(), oofem::principal_stress, oofem::FloatArray::resize(), viscoMat, and oofem::FloatArray::zero().

giveRealStressVector()

void oofem::ConcreteFCMViscoElastic::giveRealStressVector ( FloatArray & answer, GaussPoint * gp, const FloatArray & reducedStrain, TimeStep * tStep ) const

overridevirtual

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

Parameters

answer	Stress vector in reduced form. For large deformations it is treated as the second Piola-Kirchoff stress.
gp	Integration point.
reducedStrain	Strain vector in reduced form. For large deformations it is treated as the Green-Lagrange strain.
tStep	Current time step (most models are able to respond only when tStep is current time step).

Todo

Move this to StructuralCrossSection ?

Reimplemented from oofem::StructuralMaterial.

Definition at line 91 of file concretefcmviscoelastic.C.

References oofem::FloatArray::beProductOf(), oofem::FEMComponent::domain, oofem::FCMMaterialStatus::giveL2GStrainVectorTransformationMtrx(), oofem::GaussPoint::giveMaterialStatus(), oofem::FCMMaterialStatus::giveNumberOfTempCracks(), oofem::FCMMaterial::giveRealStressVector(), oofem::RheoChainMaterial::giveRealStressVector(), oofem::ConcreteFCMViscoElasticStatus::giveSlaveGaussPointVisco(), oofem::FCMMaterialStatus::giveTempCrackStrainVector(), oofem::FloatArray::subtract(), and viscoMat.

giveStatus()

MaterialStatus * oofem::ConcreteFCMViscoElastic::giveStatus ( GaussPoint * gp ) const

overridevirtual

Returns material status of receiver in given integration point. If status does not exist yet, it is created using CreateStatus member function.

Parameters

gp	Returns reference to material status belonging to integration point gp.

Returns

Material status associated with given integration point.

Reimplemented from oofem::ConcreteFCM.

Definition at line 207 of file concretefcmviscoelastic.C.

References oofem::RandomMaterialExtensionInterface::_generateStatusVariables(), CreateStatus(), oofem::GaussPoint::giveMaterialStatus(), oofem::GaussPoint::hasMaterialStatus(), and oofem::GaussPoint::setMaterialStatus().

giveTensileStrength()

double oofem::ConcreteFCMViscoElastic::giveTensileStrength ( GaussPoint * gp, TimeStep * tStep ) const

overrideprotectedvirtual

comutes tensile strength

Reimplemented from oofem::ConcreteFCM.

Definition at line 220 of file concretefcmviscoelastic.C.

References fib, fib_fcm28, fib_s, oofem::ConcreteFCM::Ft, giveEquivalentTime(), oofem::GaussPoint::giveMaterialStatus(), oofem::FCMMaterialStatus::giveNumberOfTempCracks(), oofem::ConcreteFCM::giveTensileStrength(), oofem::ConcreteFCMViscoElasticStatus::giveTensileStrength(), oofem::ConcreteFCMViscoElasticStatus::setTensileStrength(), stiffnessFactor, and timeFactor.

Referenced by giveFractureEnergy(), and giveIPValue().

initializeFrom()

void oofem::ConcreteFCMViscoElastic::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::ConcreteFCM.

Definition at line 50 of file concretefcmviscoelastic.C.

References _IFT_ConcreteFCMViscoElastic_fcm28, _IFT_ConcreteFCMViscoElastic_fib_s, _IFT_ConcreteFCMViscoElastic_ft28, _IFT_ConcreteFCMViscoElastic_gf28, _IFT_ConcreteFCMViscoElastic_stiffnessFactor, _IFT_ConcreteFCMViscoElastic_timedepfracturing, _IFT_ConcreteFCMViscoElastic_timeFactor, _IFT_ConcreteFCMViscoElastic_viscoMat, fib, fib_fcm28, fib_s, oofem::ConcreteFCM::Ft, oofem::ConcreteFCM::Gf, oofem::InputRecord::hasField(), oofem::ConcreteFCM::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, OOFEM_ERROR, oofem::Material::propertyDictionary, stiffnessFactor, tAlpha, timeFactor, and viscoMat.

Member Data Documentation

fib

bool oofem::ConcreteFCMViscoElastic::fib = false

protected

Definition at line 123 of file concretefcmviscoelastic.h.

Referenced by giveFractureEnergy(), giveTensileStrength(), and initializeFrom().

fib_fcm28

double oofem::ConcreteFCMViscoElastic::fib_fcm28 = 0.

protected

Definition at line 125 of file concretefcmviscoelastic.h.

Referenced by giveFractureEnergy(), giveTensileStrength(), and initializeFrom().

fib_s

double oofem::ConcreteFCMViscoElastic::fib_s = 0.

protected

Definition at line 124 of file concretefcmviscoelastic.h.

Referenced by giveTensileStrength(), and initializeFrom().

stiffnessFactor

double oofem::ConcreteFCMViscoElastic::stiffnessFactor = 0.

protected

Definition at line 131 of file concretefcmviscoelastic.h.

Referenced by giveFractureEnergy(), giveTensileStrength(), and initializeFrom().

timeFactor

double oofem::ConcreteFCMViscoElastic::timeFactor = 0.

protected

Definition at line 128 of file concretefcmviscoelastic.h.

Referenced by giveTensileStrength(), and initializeFrom().

viscoMat

int oofem::ConcreteFCMViscoElastic::viscoMat = 0

protected

number of the viscoelastic material

Definition at line 121 of file concretefcmviscoelastic.h.

Referenced by checkConsistency(), computeOverallElasticShearModulus(), computeOverallElasticStiffness(), computeStressIndependentStrainVector(), giveEquivalentTime(), giveIPValue(), giveRealStressVector(), and initializeFrom().

---

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

• concretefcmviscoelastic.h
• concretefcmviscoelastic.C