#include <fcm.h>

Inheritance diagram for oofem::FCMMaterial:

Collaboration diagram for oofem::FCMMaterial:

Public Member Functions
	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 ()
virtual FloatArray	computeStressIndependentStrainVector (GaussPoint gp, TimeStep tStep, ValueModeType mode) const
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 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.
virtual const char *	giveInputRecordName () const =0
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).

Protected Member Functions
virtual double	giveTensileStrength (GaussPoint gp, TimeStep tStep) const =0
	comutes tensile strength
virtual void	checkSnapBack (GaussPoint gp, TimeStep tStep, int crack) const =0
	checks possible snap-back
virtual void	updateCrackStatus (GaussPoint *gp) const
	updates crack statuses
virtual double	giveNormalCrackingStress (GaussPoint gp, TimeStep tStep, double eps_cr, int i) const =0
	computes normal stress associated with i-th crack direction
virtual double	giveCharacteristicElementLength (GaussPoint *gp, const FloatArray &crackPlaneNormal) const
	returns characteristic element length in given direction
virtual double	giveCrackingModulus (MatResponseMode rMode, GaussPoint gp, TimeStep tStep, int i) const =0
	returns stiffness in the normal direction of the i-th crack
virtual double	giveCrackingModulusInTension (MatResponseMode rMode, GaussPoint gp, TimeStep tStep, int i) const =0
	returns stiffness in the normal direction of the i-th crack
virtual double	computeEffectiveShearModulus (GaussPoint gp, TimeStep tStep, int i) const =0
	returns Geff which is necessary in the global stiffness matrix
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 double	computeD2ModulusForCrack (GaussPoint gp, TimeStep tStep, int icrack) const =0
	shear modulus for a given crack plane (1, 2, 3)
virtual double	computeNumerD2ModulusForCrack (GaussPoint gp, TimeStep tStep, int icrack) const =0
	shear modulus for numerical purpose (stiffness matrix) for a given crack plane (1, 2, 3)
virtual double	maxShearStress (GaussPoint gp, TimeStep tStep, int i) const =0
	computes the maximum value of the shear stress; if the shear stress exceeds this value, it is cropped
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	computeOverallElasticStiffness (GaussPoint gp, TimeStep tStep) const
	returns overall Young's modulus
virtual double	computeOverallElasticShearModulus (GaussPoint gp, TimeStep tStep) const
	returns overall shear modulus
virtual double	givePoissonsRatio () const
	returns Poisson's ratio

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

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.

Detailed Description

This class implements a Fixed Crack Model for fracture (after initiation the crack directions cannot rotate). After the onset of cracking, the cracks can still transfer tractions (i.e. cohesive crack approach) in both normal and shear directions. Cracks can develop only in mutually perpendicular directions. In elastic state this model is isotropic. This class is purely abstract, it can be used only in the derived classes (e.g. ConcreteFCM)

Definition at line 196 of file fcm.h.

Constructor & Destructor Documentation

◆ FCMMaterial()

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

Definition at line 46 of file fcm.C.

References oofem::ECSM_Unknown, ecsMethod, linearElasticMaterial, and oofem::StructuralMaterial::StructuralMaterial().

Referenced by oofem::ConcreteFCM::ConcreteFCM(), give3dMaterialStiffnessMatrix(), givePlaneStrainStiffMtrx(), givePlaneStressStiffMtrx(), giveRealStressVector_1d(), giveRealStressVector_2dBeamLayer(), giveRealStressVector_3d(), giveRealStressVector_PlaneStrain(), giveRealStressVector_PlaneStress(), and giveRealStressVector_PlateLayer().

Member Function Documentation

◆ checkSnapBack()

virtual void oofem::FCMMaterial::checkSnapBack	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	crack ) const

protectedpure virtual

checks possible snap-back

Implemented in oofem::ConcreteFCM, and oofem::FRCFCM.

Referenced by initializeCrack().

◆ checkStrengthCriterion()

bool oofem::FCMMaterial::checkStrengthCriterion	(	FloatMatrix &	newBase,
		const FloatArray &	globalStress,
		GaussPoint *	gp,
		TimeStep *	tStep,
		int	nCrack ) const

protectedvirtual

checks if the globalStress does not exceed strength in the direction of newBase for n-th crack

Definition at line 1301 of file fcm.C.

Referenced by giveRealStressVector().

◆ computeD2ModulusForCrack()

virtual double oofem::FCMMaterial::computeD2ModulusForCrack	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	icrack ) const

protectedpure virtual

shear modulus for a given crack plane (1, 2, 3)

Implemented in oofem::ConcreteFCM, and oofem::FRCFCM.

Referenced by computeShearStiffnessRedistributionFactor(), and computeTotalD2Modulus().

◆ computeEffectiveShearModulus()

virtual double oofem::FCMMaterial::computeEffectiveShearModulus	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	i ) const

protectedpure virtual

returns Geff which is necessary in the global stiffness matrix

Implemented in oofem::ConcreteFCM, and oofem::FRCFCM.

Referenced by giveMaterialStiffnessMatrix().

◆ computeMaxNormalCrackOpening()

double oofem::FCMMaterial::computeMaxNormalCrackOpening	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	i ) const

virtual

uses maximum equilibrated cracking strain and characteristic length to obtain the maximum reached crack opening

Definition at line 1149 of file fcm.C.

References oofem::FCMMaterialStatus::giveCharLength(), oofem::FCMMaterialStatus::giveMaxCrackStrain(), giveNumberOfCracksInDirection(), oofem::FCMMaterialStatus::giveNumberOfTempCracks(), oofem::Material::giveStatus(), oofem::max(), and N.

Referenced by oofem::ConcreteFCM::computeD2ModulusForCrack(), oofem::FRCFCM::computeTempDamage(), and oofem::ConcreteFCM::maxShearStress().

◆ computeNormalCrackOpening()

double oofem::FCMMaterial::computeNormalCrackOpening	(	GaussPoint *	gp,
		int	i ) const

virtual

uses temporary cracking strain and characteristic length to obtain the crack opening

Definition at line 1130 of file fcm.C.

References oofem::FCMMaterialStatus::giveCharLength(), giveNumberOfCracksInDirection(), oofem::FCMMaterialStatus::giveNumberOfTempCracks(), oofem::Material::giveStatus(), oofem::FCMMaterialStatus::giveTempCrackStrain(), oofem::max(), and N.

Referenced by oofem::ConcreteFCM::computeD2ModulusForCrack(), oofem::FRCFCMNL::computeNonlocalStressInFibers(), oofem::FRCFCMNL::computeNonlocalStressInFibersInUncracked(), and oofem::ConcreteFCM::maxShearStress().

◆ computeNumerD2Modulus()

double oofem::FCMMaterial::computeNumerD2Modulus	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	i ) const

protectedvirtual

shear modulus in a STIFFNESS MATRIX for a given shear direction (4, 5, 6)

Definition at line 1927 of file fcm.C.

References computeNumerD2ModulusForCrack(), computeOverallElasticStiffness(), fcm_BIGNUMBER, isIntact(), isIntactForShear(), oofem::min(), multipleCrackShear, and OOFEM_ERROR.

Referenced by oofem::ConcreteFCM::computeEffectiveShearModulus().

◆ computeNumerD2ModulusForCrack()

virtual double oofem::FCMMaterial::computeNumerD2ModulusForCrack	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	icrack ) const

protectedpure virtual

shear modulus for numerical purpose (stiffness matrix) for a given crack plane (1, 2, 3)

Implemented in oofem::ConcreteFCM.

Referenced by computeNumerD2Modulus().

◆ computeOverallElasticShearModulus()

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

inlineprotectedvirtual

returns overall shear modulus

Reimplemented in oofem::ConcreteFCMViscoElastic, and oofem::FRCFCM.

Definition at line 380 of file fcm.h.

References linearElasticMaterial.

Referenced by oofem::ConcreteFCM::computeD2ModulusForCrack(), oofem::ConcreteFCM::computeEffectiveShearModulus(), giveMaterialStiffnessMatrix(), and giveRealStressVector().

◆ computeOverallElasticStiffness()

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

inlineprotectedvirtual

returns overall Young's modulus

Reimplemented in oofem::ConcreteFCMViscoElastic, and oofem::FRCFCM.

Definition at line 377 of file fcm.h.

References linearElasticMaterial.

Referenced by oofem::ConcreteFCM::checkSnapBack(), oofem::ConcreteFCM::computeD2ModulusForCrack(), computeNumerD2Modulus(), oofem::ConcreteFCM::computeNumerD2ModulusForCrack(), computeTotalD2Modulus(), oofem::ConcreteFCM::giveCrackingModulus(), oofem::ConcreteFCM::giveCrackingModulusInTension(), giveMaterialStiffnessMatrix(), and oofem::ConcreteFCM::giveNormalCrackingStress().

◆ computeShearSlipOnCrack()

double oofem::FCMMaterial::computeShearSlipOnCrack	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	i ) const

virtual

computes total shear slip on a given crack plane (i = 1, 2, 3); the slip is computed from the temporary cracking strain

Definition at line 1168 of file fcm.C.

References computeShearStiffnessRedistributionFactor(), oofem::FCMMaterialStatus::giveCharLength(), oofem::GaussPoint::giveMaterialMode(), giveNumberOfCracksInDirection(), oofem::Material::giveStatus(), oofem::FCMMaterialStatus::giveTempCrackStatus(), oofem::FCMMaterialStatus::giveTempCrackStrain(), OOFEM_ERROR, and pscm_NONE.

Referenced by oofem::FRCFCM::computeTempDamage().

◆ computeShearStiffnessRedistributionFactor()

double oofem::FCMMaterial::computeShearStiffnessRedistributionFactor	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	ithCrackPlane,
		int	jthCrackDirection ) const

protectedvirtual

function calculating ratio used to split shear slips on two crack planes

Reimplemented in oofem::FRCFCM.

Definition at line 1287 of file fcm.C.

References computeD2ModulusForCrack().

Referenced by computeShearSlipOnCrack().

◆ computeTotalD2Modulus()

double oofem::FCMMaterial::computeTotalD2Modulus	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	i ) const

protectedvirtual

shear modulus for a given shear direction (4, 5, 6)

Definition at line 1856 of file fcm.C.

References computeD2ModulusForCrack(), computeOverallElasticStiffness(), fcm_BIGNUMBER, isIntact(), isIntactForShear(), oofem::min(), multipleCrackShear, and OOFEM_ERROR.

Referenced by giveRealStressVector(), and giveTotalLocalCrackedStiffnessMatrix().

◆ CreateStatus()

std::unique_ptr< MaterialStatus > oofem::FCMMaterial::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::Material.

Reimplemented in oofem::FRCFCM, and oofem::FRCFCMNL.

Definition at line 270 of file fcm.h.

◆ give()

double oofem::FCMMaterial::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::Material.

Definition at line 2017 of file fcm.C.

References linearElasticMaterial.

Referenced by giveMaterialStiffnessMatrix().

◆ give3dMaterialStiffnessMatrix()

FloatMatrixF< 6, 6 > oofem::FCMMaterial::give3dMaterialStiffnessMatrix	(	MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point.

Parameters

answer	Computed results.
mode	Material response mode.
gp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Definition at line 2267 of file fcm.C.

References FCMMaterial(), and giveMaterialStiffnessMatrix().

◆ giveCharacteristicElementLength()

double oofem::FCMMaterial::giveCharacteristicElementLength	(	GaussPoint *	gp,
		const FloatArray &	crackPlaneNormal ) const

protectedvirtual

returns characteristic element length in given direction

Definition at line 1501 of file fcm.C.

References oofem::Element::giveCharacteristicLength(), and oofem::GaussPoint::giveElement().

Referenced by initializeCrack().

◆ giveClassName()

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

inlineoverridevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Reimplemented in oofem::FRCFCMNL.

Definition at line 206 of file fcm.h.

◆ giveCrackingModulus()

virtual double oofem::FCMMaterial::giveCrackingModulus	(	MatResponseMode	rMode,
		GaussPoint *	gp,
		TimeStep *	tStep,
		int	i ) const

protectedpure virtual

returns stiffness in the normal direction of the i-th crack

Implemented in oofem::ConcreteFCM, and oofem::FRCFCM.

Referenced by giveMaterialStiffnessMatrix(), giveNormalLocalCrackedStiffnessMatrix(), giveRealStressVector(), and giveTotalLocalCrackedStiffnessMatrix().

◆ giveCrackingModulusInTension()

virtual double oofem::FCMMaterial::giveCrackingModulusInTension	(	MatResponseMode	rMode,
		GaussPoint *	gp,
		TimeStep *	tStep,
		int	i ) const

protectedpure virtual

returns stiffness in the normal direction of the i-th crack

Implemented in oofem::ConcreteFCM.

◆ giveCrackSpacing()

double oofem::FCMMaterial::giveCrackSpacing ( void ) const

protectedvirtual

returns either user-provided value of crack spacing or a value computed from composition

Definition at line 2025 of file fcm.C.

References crackSpacing.

Referenced by giveNumberOfCracksInDirection(), and oofem::ConcreteFCM::initializeFrom().

◆ giveIPValue()

int oofem::FCMMaterial::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::Material.

Reimplemented in oofem::FRCFCM, and oofem::FRCFCMNL.

Definition at line 2076 of file fcm.C.

◆ giveMaterialStiffnessMatrix()

void oofem::FCMMaterial::giveMaterialStiffnessMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

protectedvirtual

Reimplemented in oofem::FRCFCMNL.

Definition at line 1643 of file fcm.C.

References oofem::FloatMatrix::add(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beInverseOf(), oofem::FloatMatrix::beProductOf(), computeEffectiveShearModulus(), computeOverallElasticShearModulus(), computeOverallElasticStiffness(), fcm_BIGNUMBER, fcm_THRESHOLD_CRACK_STRAIN, give(), giveCrackingModulus(), oofem::FCMMaterialStatus::giveCrackStrain(), oofem::FCMMaterialStatus::giveG2LStrainVectorTransformationMtrx(), oofem::FCMMaterialStatus::giveMaxNumberOfCracks(), oofem::FCMMaterialStatus::giveNumberOfTempCracks(), givePoissonsRatio(), oofem::Material::giveStatus(), oofem::FCMMaterialStatus::giveTempCrackStrain(), linearElasticMaterial, oofem::max(), normalCoeffNumer, OOFEM_ERROR, oofem::FloatMatrix::resize(), oofem::FloatMatrix::resizeWithData(), oofem::FloatMatrix::rotatedWith(), oofem::FloatMatrix::subtract(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

Referenced by give3dMaterialStiffnessMatrix(), givePlaneStrainStiffMtrx(), and givePlaneStressStiffMtrx().

◆ giveNormalCrackingStress()

virtual double oofem::FCMMaterial::giveNormalCrackingStress	(	GaussPoint *	gp,
		TimeStep *	tStep,
		double	eps_cr,
		int	i ) const

protectedpure virtual

computes normal stress associated with i-th crack direction

Implemented in oofem::ConcreteFCM, and oofem::FRCFCM.

◆ giveNormalLocalCrackedStiffnessMatrix()

void oofem::FCMMaterial::giveNormalLocalCrackedStiffnessMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

protectedvirtual

returns local stiffness matrix in the cracks' direction (only normal components)

Definition at line 1622 of file fcm.C.

References oofem::FloatMatrix::at(), giveCrackingModulus(), oofem::FCMMaterialStatus::giveMaxNumberOfCracks(), oofem::Material::giveStatus(), and oofem::FloatMatrix::resize().

Referenced by giveRealStressVector().

◆ giveNumberOfCracksForShearDirection()

double oofem::FCMMaterial::giveNumberOfCracksForShearDirection	(	GaussPoint *	gp,
		int	i ) const

protectedvirtual

returns number of cracks for given shear direction (i = 4, 5, 6) which is treated as the maximum of the two associated normal directions

Definition at line 2051 of file fcm.C.

References giveNumberOfCracksInDirection(), oofem::max(), N, and OOFEM_ERROR.

Referenced by oofem::ConcreteFCM::computeEffectiveShearModulus().

◆ giveNumberOfCracksInDirection()

double oofem::FCMMaterial::giveNumberOfCracksInDirection	(	GaussPoint *	gp,
		int	iCrack ) const

protectedvirtual

returns number of fictiotious parallel cracks in the direction of i-th crack

Definition at line 2032 of file fcm.C.

References oofem::FCMMaterialStatus::giveCharLength(), giveCrackSpacing(), oofem::Material::giveStatus(), and N.

Referenced by oofem::ConcreteFCM::computeD2ModulusForCrack(), oofem::ConcreteFCM::computeEffectiveShearModulus(), computeMaxNormalCrackOpening(), computeNormalCrackOpening(), oofem::ConcreteFCM::computeNumerD2ModulusForCrack(), oofem::ConcreteFCM::computeResidualTensileStrength(), computeShearSlipOnCrack(), oofem::FRCFCM::computeStressInFibersInCracked(), oofem::FRCFCM::giveCrackingModulus(), oofem::ConcreteFCM::giveCrackingModulusInTension(), giveIPValue(), oofem::ConcreteFCM::giveNormalCrackingStress(), and giveNumberOfCracksForShearDirection().

◆ givePlaneStrainStiffMtrx()

FloatMatrixF< 4, 4 > oofem::FCMMaterial::givePlaneStrainStiffMtrx	(	MatResponseMode	mmode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Method for computing plane strain stiffness matrix of receiver. Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to plane strain stiffness using reduce method described above. However, this reduction is quite time consuming and if it is possible, it is recommended to overload this method and provide direct method for computing particular stiffness matrix. Note: as already described, if zero strain component is imposed (Plane strain, ..) this condition must be taken into account in geometrical relations, and corresponding component has to be included in reduced vector. (So plane strain conditions are \( \epsilon_z = \gamma_{xz} = \gamma_{yz} = 0 \), but relations for \( \epsilon_z\) and \(\sigma_z\) are included).

Parameters

answer	Stiffness matrix.
mmode	Material response mode.
gp	Integration point, which load history is used.
tStep	Time step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Definition at line 2289 of file fcm.C.

References FCMMaterial(), and giveMaterialStiffnessMatrix().

◆ givePlaneStressStiffMtrx()

FloatMatrixF< 3, 3 > oofem::FCMMaterial::givePlaneStressStiffMtrx	(	MatResponseMode	mmode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Method for computing plane stress stiffness matrix of receiver. Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to plane stress stiffness using reduce method described above. However, this reduction is quite time consuming and if it is possible, it is recommended to overload this method and provide direct method for computing particular stiffness matrix.

Parameters

answer	Stiffness matrix.
mmode	Material response mode.
gp	Integration point, which load history is used.
tStep	Time step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Definition at line 2278 of file fcm.C.

References FCMMaterial(), and giveMaterialStiffnessMatrix().

◆ givePoissonsRatio()

virtual double oofem::FCMMaterial::givePoissonsRatio ( ) const

inlineprotectedvirtual

returns Poisson's ratio

Definition at line 383 of file fcm.h.

References linearElasticMaterial.

Referenced by giveMaterialStiffnessMatrix().

◆ giveRealStressVector()

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

Reimplemented in oofem::FRCFCMNL.

Definition at line 60 of file fcm.C.

Referenced by oofem::ConcreteFCMViscoElastic::giveRealStressVector(), giveRealStressVector_1d(), giveRealStressVector_2dBeamLayer(), giveRealStressVector_3d(), giveRealStressVector_PlaneStrain(), giveRealStressVector_PlaneStress(), and giveRealStressVector_PlateLayer().

◆ giveRealStressVector_1d()

FloatArrayF< 1 > oofem::FCMMaterial::giveRealStressVector_1d	(	const FloatArrayF< 1 > &	reducedE,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

inlineoverridevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 241 of file fcm.h.

References FCMMaterial(), and giveRealStressVector().

◆ giveRealStressVector_2dBeamLayer()

FloatArrayF< 2 > oofem::FCMMaterial::giveRealStressVector_2dBeamLayer	(	const FloatArrayF< 2 > &	reducedE,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

inlineoverridevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 247 of file fcm.h.

References FCMMaterial(), and giveRealStressVector().

◆ giveRealStressVector_3d()

FloatArrayF< 6 > oofem::FCMMaterial::giveRealStressVector_3d	(	const FloatArrayF< 6 > &	strain,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

inlineoverridevirtual

Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.

Reimplemented from oofem::StructuralMaterial.

Definition at line 223 of file fcm.h.

References FCMMaterial(), and giveRealStressVector().

◆ giveRealStressVector_PlaneStrain()

FloatArrayF< 4 > oofem::FCMMaterial::giveRealStressVector_PlaneStrain	(	const FloatArrayF< 4 > &	strain,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

inlineoverridevirtual

Default implementation relies on giveRealStressVector_3d.

Reimplemented from oofem::StructuralMaterial.

Definition at line 229 of file fcm.h.

References FCMMaterial(), and giveRealStressVector().

◆ giveRealStressVector_PlaneStress()

FloatArrayF< 3 > oofem::FCMMaterial::giveRealStressVector_PlaneStress	(	const FloatArrayF< 3 > &	reducedE,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

inlineoverridevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 235 of file fcm.h.

References FCMMaterial(), and giveRealStressVector().

◆ giveRealStressVector_PlateLayer()

FloatArrayF< 5 > oofem::FCMMaterial::giveRealStressVector_PlateLayer	(	const FloatArrayF< 5 > &	reducedE,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

inlineoverridevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 253 of file fcm.h.

References FCMMaterial(), and giveRealStressVector().

◆ giveTensileStrength()

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

protectedpure virtual

comutes tensile strength

Implemented in oofem::ConcreteFCM, and oofem::ConcreteFCMViscoElastic.

Referenced by giveRealStressVector(), and isStrengthExceeded().

◆ giveTotalLocalCrackedStiffnessMatrix()

void oofem::FCMMaterial::giveTotalLocalCrackedStiffnessMatrix	(	FloatMatrix &	answer,
		MatResponseMode	rMode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

protectedvirtual

returns local stiffness matrix in the cracks' direction (total according to the material mode)

Definition at line 1593 of file fcm.C.

References oofem::FloatMatrix::at(), oofem::IntArray::at(), computeTotalD2Modulus(), giveCrackingModulus(), oofem::GaussPoint::giveMaterialMode(), oofem::IntArray::giveSize(), and oofem::FloatMatrix::resize().

Referenced by giveRealStressVector().

◆ hasMaterialModeCapability()

bool oofem::FCMMaterial::hasMaterialModeCapability ( MaterialMode mode ) const

overridevirtual

Tests if material supports material mode.

Parameters

mode	Required material mode.

Returns: Nonzero if supported, zero otherwise.

Reimplemented from oofem::Material.

Definition at line 53 of file fcm.C.

◆ initializeCrack()

void oofem::FCMMaterial::initializeCrack	(	GaussPoint *	gp,
		TimeStep *	tStep,
		FloatMatrix &	base,
		int	nCrack ) const

virtual

Definition at line 972 of file fcm.C.

Referenced by giveRealStressVector().

◆ initializeFrom()

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

Reimplemented in oofem::FRCFCM, and oofem::FRCFCMNL.

Definition at line 1974 of file fcm.C.

References _IFT_FCM_crackSpacing, _IFT_FCM_ecsm, _IFT_FCM_multipleCrackShear, _IFT_FCM_nAllowedCracks, _IFT_FCM_normalCoeffNumer, _IFT_FCM_shearCoeffNumer, crackSpacing, oofem::ECSM_Oliver1, oofem::ECSM_Oliver1modified, oofem::ECSM_Projection, oofem::ECSM_ProjectionCentered, oofem::ECSM_SquareRootOfArea, ecsMethod, fcm_BIGNUMBER, oofem::InputRecord::hasField(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, iterLimitGlobal, linearElasticMaterial, multipleCrackShear, nAllowedCracks, normalCoeffNumer, and shearCoeffNumer.

◆ isIntact()

bool oofem::FCMMaterial::isIntact	(	GaussPoint *	gp,
		int	icrack ) const

protectedvirtual

returns true for closed or no crack (i = 1, 2, 3)

Definition at line 1042 of file fcm.C.

References oofem::FCMMaterialStatus::giveCrackStatus(), oofem::Material::giveStatus(), oofem::FCMMaterialStatus::giveTempCrackStatus(), OOFEM_ERROR, pscm_CLOSED, and pscm_NONE.

Referenced by oofem::ConcreteFCM::computeD2ModulusForCrack(), oofem::FRCFCM::computeD2ModulusForCrack(), oofem::FRCFCM::computeEffectiveShearModulus(), computeNumerD2Modulus(), oofem::ConcreteFCM::computeNumerD2ModulusForCrack(), oofem::FRCFCM::computeTempDamage(), computeTotalD2Modulus(), oofem::FRCFCM::estimateD2ModulusForCrack(), and oofem::FRCFCM::maxShearStress().

◆ isIntactForShear()

bool oofem::FCMMaterial::isIntactForShear	(	GaussPoint *	gp,
		int	i ) const

protectedvirtual

returns true for closed or no cracks associated to given shear direction (i = 4, 5, 6)

Definition at line 1060 of file fcm.C.

References oofem::FCMMaterialStatus::giveCrackStatus(), oofem::Material::giveStatus(), OOFEM_ERROR, and pscm_NONE.

Referenced by oofem::ConcreteFCM::computeEffectiveShearModulus(), oofem::FRCFCM::computeEffectiveShearModulus(), computeNumerD2Modulus(), computeTotalD2Modulus(), and oofem::ConcreteFCM::maxShearStress().

◆ isStrengthExceeded()

bool oofem::FCMMaterial::isStrengthExceeded	(	const FloatMatrix &	base,
		GaussPoint *	gp,
		TimeStep *	tStep,
		int	iCrack,
		double	trialStress ) const

protectedvirtual

compares trial stress with strength. Returns true if the strength is exceeded. Function oveloaded in the nonlocal approach for the fiber reinforced composites

Reimplemented in oofem::FRCFCM, and oofem::FRCFCMNL.

Definition at line 1275 of file fcm.C.

References giveTensileStrength().

Referenced by checkStrengthCriterion(), and giveRealStressVector().

◆ isThisShearComponent()

bool oofem::FCMMaterial::isThisShearComponent	(	GaussPoint *	gp,
		int	component ) const

protectedvirtual

returns true if current component is associated with shear

Definition at line 1094 of file fcm.C.

References oofem::GaussPoint::giveMaterialMode(), and OOFEM_ERROR.

Referenced by giveRealStressVector().

◆ maxShearStress()

virtual double oofem::FCMMaterial::maxShearStress	(	GaussPoint *	gp,
		TimeStep *	tStep,
		int	i ) const

protectedpure virtual

computes the maximum value of the shear stress; if the shear stress exceeds this value, it is cropped

Implemented in oofem::ConcreteFCM, and oofem::FRCFCM.

Referenced by giveRealStressVector().

◆ updateCrackStatus()

void oofem::FCMMaterial::updateCrackStatus ( GaussPoint * gp ) const

protectedvirtual

updates crack statuses

Definition at line 1508 of file fcm.C.

References oofem::FloatArray::at(), oofem::FCMMaterialStatus::giveCrackStatus(), oofem::FCMMaterialStatus::giveMaxCrackStrainVector(), oofem::FCMMaterialStatus::giveMaxNumberOfCracks(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), oofem::FCMMaterialStatus::giveTempCrackStatus(), oofem::FCMMaterialStatus::giveTempCrackStrainVector(), OOFEM_ERROR, pscm_CLOSED, pscm_JUST_INIT, pscm_NONE, pscm_SOFTENING, pscm_UNLO_RELO, oofem::FCMMaterialStatus::setTempCrackStatus(), and oofem::FCMMaterialStatus::setTempMaxCrackStrain().

Referenced by giveRealStressVector().

Member Data Documentation

◆ crackSpacing

double oofem::FCMMaterial::crackSpacing

protected

value of crack spacing (allows to "have" more parallel cracks in one direction if the element size exceeds user-defined or computed crack spacing).

Definition at line 343 of file fcm.h.

Referenced by giveCrackSpacing(), initializeFrom(), and oofem::FRCFCM::initializeFrom().

◆ ecsMethod

ElementCharSizeMethod oofem::FCMMaterial::ecsMethod

protected

Method used for evaluation of characteristic element size.

Definition at line 278 of file fcm.h.

Referenced by FCMMaterial(), and initializeFrom().

◆ iterLimitGlobal

int oofem::FCMMaterial::iterLimitGlobal = 20

protected

Definition at line 283 of file fcm.h.

Referenced by giveRealStressVector(), and initializeFrom().

◆ linearElasticMaterial

IsotropicLinearElasticMaterial oofem::FCMMaterial::linearElasticMaterial

protected

Definition at line 199 of file fcm.h.

Referenced by oofem::ConcreteFCMViscoElastic::checkConsistency(), oofem::FRCFCM::computeD2ModulusForCrack(), oofem::ConcreteFCMViscoElastic::computeOverallElasticShearModulus(), computeOverallElasticShearModulus(), oofem::FRCFCM::computeOverallElasticShearModulus(), computeOverallElasticStiffness(), oofem::FRCFCM::computeOverallElasticStiffness(), oofem::FRCFCM::estimateD2ModulusForCrack(), FCMMaterial(), give(), giveMaterialStiffnessMatrix(), givePoissonsRatio(), initializeFrom(), oofem::FRCFCM::isStrengthExceeded(), and oofem::FRCFCM::maxShearStress().

◆ multipleCrackShear

bool oofem::FCMMaterial::multipleCrackShear

protected

if true = takes shear compliance of all cracks, false = only dominant crack contribution, default value is false

Definition at line 281 of file fcm.h.

Referenced by oofem::ConcreteFCM::computeEffectiveShearModulus(), oofem::FRCFCM::computeEffectiveShearModulus(), computeNumerD2Modulus(), computeTotalD2Modulus(), and initializeFrom().

◆ nAllowedCracks

int oofem::FCMMaterial::nAllowedCracks

protected

allowed number of cracks (user-defined)

Definition at line 275 of file fcm.h.

Referenced by giveRealStressVector(), and initializeFrom().

◆ normalCoeffNumer

double oofem::FCMMaterial::normalCoeffNumer

protected

minimum ratio of effective normal stiffness vs. overall modulus - just for numerical purpose

Definition at line 349 of file fcm.h.

Referenced by giveMaterialStiffnessMatrix(), and initializeFrom().

◆ shearCoeffNumer

double oofem::FCMMaterial::shearCoeffNumer

protected

minimum ratio of effective shear modulus vs. shear modulus - just for numerical purpose

Definition at line 346 of file fcm.h.

Referenced by initializeFrom().

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

fcm.h
fcm.C

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FCMMaterial()

Member Function Documentation

◆ checkSnapBack()

◆ checkStrengthCriterion()

◆ computeD2ModulusForCrack()

◆ computeEffectiveShearModulus()

◆ computeMaxNormalCrackOpening()

◆ computeNormalCrackOpening()

◆ computeNumerD2Modulus()

◆ computeNumerD2ModulusForCrack()

◆ computeOverallElasticShearModulus()

◆ computeOverallElasticStiffness()

◆ computeShearSlipOnCrack()

◆ computeShearStiffnessRedistributionFactor()

◆ computeTotalD2Modulus()

◆ CreateStatus()

◆ give()

◆ give3dMaterialStiffnessMatrix()

◆ giveCharacteristicElementLength()

◆ giveClassName()

◆ giveCrackingModulus()

◆ giveCrackingModulusInTension()

◆ giveCrackSpacing()

◆ giveIPValue()

◆ giveMaterialStiffnessMatrix()

◆ giveNormalCrackingStress()

◆ giveNormalLocalCrackedStiffnessMatrix()

◆ giveNumberOfCracksForShearDirection()

◆ giveNumberOfCracksInDirection()

◆ givePlaneStrainStiffMtrx()

◆ givePlaneStressStiffMtrx()

◆ givePoissonsRatio()

◆ giveRealStressVector()

◆ giveRealStressVector_1d()

◆ giveRealStressVector_2dBeamLayer()

◆ giveRealStressVector_3d()

◆ giveRealStressVector_PlaneStrain()

◆ giveRealStressVector_PlaneStress()

◆ giveRealStressVector_PlateLayer()

◆ giveTensileStrength()

◆ giveTotalLocalCrackedStiffnessMatrix()

◆ hasMaterialModeCapability()

◆ initializeCrack()

◆ initializeFrom()

◆ isIntact()

◆ isIntactForShear()

◆ isStrengthExceeded()

◆ isThisShearComponent()

◆ maxShearStress()

◆ updateCrackStatus()

Member Data Documentation

◆ crackSpacing

◆ ecsMethod

◆ iterLimitGlobal

◆ linearElasticMaterial

◆ multipleCrackShear

◆ nAllowedCracks

◆ normalCoeffNumer

◆ shearCoeffNumer