#include <mazarsmodel.h>

Inheritance diagram for oofem::MazarsMaterial:

Collaboration diagram for oofem::MazarsMaterial:

Public Member Functions
	MazarsMaterial (int n, Domain *d)
	Constructor.
const char *	giveInputRecordName () const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
double	computeEquivalentStrain (const FloatArray &strain, GaussPoint gp, TimeStep tStep) const override
double	computeDamageParam (double kappa, const FloatArray &strain, GaussPoint *gp) const override
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
Public Member Functions inherited from oofem::IsotropicDamageMaterial1
	IsotropicDamageMaterial1 (int n, Domain *d)
	Constructor.
virtual	~IsotropicDamageMaterial1 ()
	Destructor.
void	giveInputRecord (DynamicInputRecord &input) override
bool	isCrackBandApproachUsed () const
void	computeEta (FloatArray &answer, const FloatArray &strain, GaussPoint gp, TimeStep tStep) const override
double	computeDamageParamForCohesiveCrack (double kappa, GaussPoint *gp) const
double	damageFunction (double kappa, GaussPoint *gp) const
double	damageFunctionPrime (double kappa, GaussPoint *gp) const override
double	complianceFunction (double kappa, GaussPoint *gp) const
double	evaluatePermanentStrain (double kappa, double omega) const override
Interface *	giveInterface (InterfaceType it) override
int	MMI_map (GaussPoint gp, Domain oldd, TimeStep *tStep) override
int	MMI_update (GaussPoint gp, TimeStep tStep, FloatArray *estrain=nullptr) override
int	MMI_finish (TimeStep *tStep) override
MaterialStatus *	giveStatus (GaussPoint *gp) const override
double	give (int aProperty, GaussPoint *gp) const override
bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode) const override
void	restoreContext (DataStream &stream, ContextMode mode) override
void	saveContext (DataStream &stream, ContextMode mode) override
Public Member Functions inherited from oofem::IsotropicDamageMaterial
	IsotropicDamageMaterial (int n, Domain *d)
	Constructor.
virtual	~IsotropicDamageMaterial ()
	Destructor.
bool	hasMaterialModeCapability (MaterialMode mode) const override
LinearElasticMaterial *	giveLinearElasticMaterial () const
	Returns reference to undamaged (bulk) material.
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const override
void	giveRealStressVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep) const override
FloatArrayF< 6 >	giveRealStressVector_3d (const FloatArrayF< 6 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
FloatArrayF< 4 >	giveRealStressVector_PlaneStrain (const FloatArrayF< 4 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector_3d.
FloatArray	giveRealStressVector_StressControl (const FloatArray &strain, const IntArray &strainControl, GaussPoint gp, TimeStep tStep) const override
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
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.
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
FloatArrayF< 6 >	giveThermalDilatationVector (GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 1, 1 >	give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) 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_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< 2 >	giveRealStressVector_2dBeamLayer (const FloatArrayF< 2 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 5 >	giveRealStressVector_PlateLayer (const FloatArrayF< 5 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_Fiber (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_2dPlateSubSoil (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation is not provided.
virtual FloatArrayF< 6 >	giveRealStressVector_3dBeamSubSoil (const FloatArrayF< 6 > &reducedE, GaussPoint gp, TimeStep tStep) const
virtual FloatArrayF< 9 >	giveFirstPKStressVector_3d (const FloatArrayF< 9 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
virtual FloatArrayF< 5 >	giveFirstPKStressVector_PlaneStrain (const FloatArrayF< 5 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_3d.
virtual FloatArray	giveFirstPKStressVector_StressControl (const FloatArray &reducedvF, const IntArray &FControl, GaussPoint gp, TimeStep tStep) const
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArrayF< 4 >	giveFirstPKStressVector_PlaneStress (const FloatArrayF< 4 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
virtual FloatArrayF< 1 >	giveFirstPKStressVector_1d (const FloatArrayF< 1 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
virtual void	giveCauchyStressVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveCauchyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveCauchyStressVector_PlaneStress (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveCauchyStressVector_1d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
double	giveReferenceTemperature ()
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 >	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 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)
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 ()
Public Member Functions inherited from oofem::MaterialModelMapperInterface
	MaterialModelMapperInterface ()
	Constructor.
virtual	~MaterialModelMapperInterface ()
	Destructor.

Protected Types
enum	mazarsModelVariant { maz_original , maz_modTension }
	Model variants. More...
Protected Types inherited from oofem::IsotropicDamageMaterial1
enum	EquivStrainType { EST_Mazars =0 , EST_Rankine_Smooth =1 , EST_ElasticEnergy =2 , EST_Mises =3 , EST_Rankine_Standard =4 , EST_ElasticEnergyPositiveStress =5 , EST_ElasticEnergyPositiveStrain =6 , EST_Griffith =7 , EST_Unknown = 100 }
enum	SofteningType { ST_Unknown , ST_Exponential , ST_Linear , ST_Mazars , ST_Smooth , ST_SmoothExtended , ST_Exponential_Cohesive_Crack , ST_Linear_Cohesive_Crack , ST_BiLinear_Cohesive_Crack , ST_Disable_Damage , ST_PowerExponential , ST_DoubleExponential , ST_Hordijk_Cohesive_Crack , ST_ModPowerExponential , ST_Trilinear_Cohesive_Crack }
Protected Types inherited from oofem::IsotropicDamageMaterial
enum	loaUnloCriterium { idm_strainLevelCR , idm_damageLevelCR }

Protected Member Functions
void	initDamaged (double kappa, FloatArray &totalStrainVector, GaussPoint *gp) const override
int	giveNumberOfSpatialDimensions (GaussPoint *gp) const
void	giveNormalBlockOfElasticCompliance (FloatMatrix &answer, GaussPoint *gp) const
double	computeGt (double kappa, GaussPoint *gp) const
double	computeGc (double kappa, GaussPoint *gp) const
Protected Member Functions inherited from oofem::IsotropicDamageMaterial
FloatMatrixF< 3, 3 >	givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 4, 4 >	givePlaneStrainStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
Protected Member Functions inherited from oofem::RandomMaterialExtensionInterface
void	_generateStatusVariables (GaussPoint *) const

Protected Attributes
double	E = 0.
	Elastic parameters.
double	nu = 0.
double	At = 0.
	Model parameters related to the shape of uniaxial stress-strain diagrams.
double	Bt = 0.
double	Ac = 0.
double	Bc = 0.
double	hReft = 0.
	Reference elem-length for objectivity.
double	hRefc = 0.
double	beta = 0.
	Beta coefficient reducing the effect of shear; default val = 1.06.
enum oofem::MazarsMaterial::mazarsModelVariant	modelVersion
Protected Attributes inherited from oofem::IsotropicDamageMaterial1
double	e0 = 0.
	Equivalent strain at stress peak (or a similar parameter).
double	ef = 0.
	Determines ductility -> corresponds to fracturing strain.
double	wf = 0.
	Determines ductility -> corresponds to crack opening in the cohesive crack model.
double	gf = 0.
double	gft = 0.
	Determines the softening for the bilinear law -> corresponds to the total fracture energy.
double	ek = 0.
	Determines the softening for the bilinear law -> corresponds to the strain at the knee point.
double	wk = 0.
	Determines the softening for the bilinear law -> corresponds to the crack opening at the knee point.
double	sk = 0.
	Determines the softening for the bilinear law -> corresponds to the stress at the knee point.
double	c1 = 3.
	Parameters used in Hordijk's softening law.
double	c2 = 6.93
double	w_k = 0.
	Parameters used in Trilinear_Cohesive_Crack softening law.
double	w_r = 0.
double	w_f = 0.
double	f_k = 0.
double	f_r = 0.
EquivStrainType	equivStrainType = EST_Unknown
	Parameter specifying the definition of equivalent strain.
double	k = 0.
	Parameter used in Mises definition of equivalent strain.
double	griff_n = 8.
	Parameter used in Griffith's criterion.
int	damageLaw = 0
	Temporary parameter reading type of softening law, used in other isotropic damage material models.
SofteningType	softType = ST_Unknown
	Parameter specifying the type of softening (damage law).
double	At = 0.
	Parameters used in Mazars damage law.
double	Bt = 0.
double	md = 1.
	Parameter used in "smooth damage law".
double	e1 = 0.
	Parameters used if softType = 7 (extended smooth damage law).
double	e2 = 0.
double	s1 = 0.
double	nd = 0.
int	checkSnapBack = 0
	Check possible snap back flag.
double	ep = 0.
	auxiliary input variablesfor softType == ST_SmoothExtended
double	ft = 0.
double	ps_alpha = 0.
	Parameters used by the model with permanent strain.
double	ps_H = 0.
ElementCharSizeMethod	ecsMethod = ECSM_Unknown
	Method used for evaluation of characteristic element size.
Set *	sourceElemSet = nullptr
	Cached source element set used to map internal variables (adaptivity), created on demand.
Protected Attributes inherited from oofem::IsotropicDamageMaterial
double	tempDillatCoeff = 0.
	Coefficient of thermal dilatation.
double	maxOmega = 0.999999
	Maximum limit on omega. The purpose is elimination of a too compliant material which may cause convergence problems. Set to something like 0.99 if needed.
int	permStrain = 0
	Indicator of the type of permanent strain formulation (0 = standard damage with no permanent strain).
LinearElasticMaterial *	linearElasticMaterial = nullptr
	Reference to bulk (undamaged) material.
enum oofem::IsotropicDamageMaterial::loaUnloCriterium	llcriteria = idm_strainLevelCR
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::IsotropicDamageMaterial1
static void	computeStrainInvariants (const FloatArray &strainVector, double &I1e, double &J2e)
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.
Static Protected Attributes inherited from oofem::IsotropicDamageMaterial1
static MMAContainingElementProjection	mapper
	Mapper used to map internal variables in adaptivity.

Detailed Description

This class implements a Mazars damage model form concrete. It introduces two damage parameters omega_t and omega_c that are computed from the same equivalent strain using two different damage functions g_t and g_c. The g_t is identified from the uniaxial tension tests, while g_c from compressive test. The damage parameter for general stress states omega is obtained as a linear combination of omega_t and omega_c.

Definition at line 92 of file mazarsmodel.h.

Member Enumeration Documentation

◆ mazarsModelVariant

enum oofem::MazarsMaterial::mazarsModelVariant

protected

Model variants.

Enumerator
maz_original
maz_modTension

Definition at line 105 of file mazarsmodel.h.

Constructor & Destructor Documentation

◆ MazarsMaterial()

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

Constructor.

Definition at line 51 of file mazarsmodel.C.

References oofem::IsotropicDamageMaterial::idm_damageLevelCR, oofem::IsotropicDamageMaterial1::IsotropicDamageMaterial1(), and oofem::IsotropicDamageMaterial::llcriteria.

Referenced by oofem::MazarsNLMaterial::MazarsNLMaterial().

Member Function Documentation

◆ computeDamageParam()

double oofem::MazarsMaterial::computeDamageParam	(	double	kappa,
		const FloatArray &	strain,
		GaussPoint *	gp ) const

overridevirtual

Computes the value of damage parameter omega, based on given value of equivalent strain.

Parameters

[out]	omega	Contains result.
	kappa	Equivalent strain measure.
	strain	Total strain in full form.
	gp	Integration point.

Reimplemented from oofem::IsotropicDamageMaterial1.

Definition at line 215 of file mazarsmodel.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::beInverseOf(), oofem::FloatArray::beProductOf(), beta, computeGc(), computeGt(), oofem::StructuralMaterial::computePrincipalValues(), oofem::IsotropicDamageMaterial1::e0, gc, giveNormalBlockOfElasticCompliance(), giveNumberOfSpatialDimensions(), oofem::FloatArray::isEmpty(), nu, oofem::principal_strain, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

◆ computeEquivalentStrain()

double oofem::MazarsMaterial::computeEquivalentStrain	(	const FloatArray &	strain,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Computes the equivalent strain measure from given strain vector (full form).

Parameters

[out]	kappa	Return parameter, containing the corresponding equivalent strain.
	strain	Total strain vector in full form.
	gp	Integration point.
	tStep	Time step.

Reimplemented from oofem::IsotropicDamageMaterial1.

Reimplemented in oofem::MazarsNLMaterial.

Definition at line 115 of file mazarsmodel.C.

References oofem::FloatArray::at(), oofem::StructuralMaterial::computePrincipalValues(), oofem::GaussPoint::giveMaterialMode(), giveNumberOfSpatialDimensions(), oofem::FloatArray::isEmpty(), nu, oofem::principal_strain, and oofem::FloatArray::resize().

◆ computeGc()

double oofem::MazarsMaterial::computeGc	(	double	kappa,
		GaussPoint *	gp ) const

protected

Definition at line 395 of file mazarsmodel.C.

References _MAZARS_MODEL_ITER_TOL, _MAZARS_MODEL_MAX_ITER, Ac, Bc, oofem::IsotropicDamageMaterial1::e0, gc, oofem::MazarsMaterialStatus::giveLec(), oofem::IsotropicDamageMaterial1::giveStatus(), hRefc, and OOFEM_ERROR.

Referenced by computeDamageParam().

◆ computeGt()

double oofem::MazarsMaterial::computeGt	(	double	kappa,
		GaussPoint *	gp ) const

protected

Definition at line 348 of file mazarsmodel.C.

References _MAZARS_MODEL_ITER_TOL, _MAZARS_MODEL_MAX_ITER, At, Bt, oofem::IsotropicDamageMaterial1::e0, oofem::IsotropicDamageMaterial1::ef, oofem::IsotropicDamageMaterialStatus::giveLe(), oofem::IsotropicDamageMaterial1::giveStatus(), hReft, maz_modTension, modelVersion, and OOFEM_ERROR.

Referenced by computeDamageParam().

◆ CreateStatus()

std::unique_ptr< MaterialStatus > oofem::MazarsMaterial::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::IsotropicDamageMaterial1.

Reimplemented in oofem::MazarsNLMaterial.

Definition at line 119 of file mazarsmodel.h.

◆ giveClassName()

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

inlineoverridevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::IsotropicDamageMaterial1.

Reimplemented in oofem::MazarsNLMaterial.

Definition at line 112 of file mazarsmodel.h.

◆ giveInputRecordName()

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

inlineoverridevirtual

Returns: Input record name of the receiver.

Reimplemented from oofem::IsotropicDamageMaterial1.

Definition at line 111 of file mazarsmodel.h.

References _IFT_MazarsMaterial_Name.

◆ giveNormalBlockOfElasticCompliance()

void oofem::MazarsMaterial::giveNormalBlockOfElasticCompliance	(	FloatMatrix &	answer,
		GaussPoint *	gp ) const

protected

Definition at line 199 of file mazarsmodel.C.

References oofem::FloatMatrix::at(), E, giveNumberOfSpatialDimensions(), nu, and oofem::FloatMatrix::resize().

Referenced by computeDamageParam().

◆ giveNumberOfSpatialDimensions()

int oofem::MazarsMaterial::giveNumberOfSpatialDimensions ( GaussPoint * gp ) const

protected

Definition at line 187 of file mazarsmodel.C.

References oofem::GaussPoint::giveMaterialMode().

Referenced by computeDamageParam(), computeEquivalentStrain(), and giveNormalBlockOfElasticCompliance().

◆ initDamaged()

void oofem::MazarsMaterial::initDamaged	(	double	kappa,
		FloatArray &	totalStrainVector,
		GaussPoint *	gp ) const

overrideprotectedvirtual

Perfoms initialization, when damage first appear. The Le characteristic length is computed from the direction of largest positive principal strain and stored in corresponding status.

Parameters

kappa	Scalar measure of strain level,
totalStrainVector	Current total strain vector,
gp	Integration point,

Reimplemented from oofem::IsotropicDamageMaterial1.

Reimplemented in oofem::MazarsNLMaterial.

Definition at line 426 of file mazarsmodel.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::StructuralMaterial::computePrincipalValDir(), oofem::IsotropicDamageMaterial1::e0, oofem::Element::giveCharacteristicLength(), oofem::IsotropicDamageMaterialStatus::giveDamage(), oofem::GaussPoint::giveElement(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatArray::giveSize(), oofem::IsotropicDamageMaterial1::giveStatus(), oofem::principal_strain, oofem::FloatMatrix::resize(), oofem::IsotropicDamageMaterialStatus::setLe(), oofem::MazarsMaterialStatus::setLec(), and oofem::FloatArray::zero().

◆ initializeFrom()

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