This class implements a Mazars damage model form concrete. More...

#include <mazarsmodel.h>

Inheritance diagram for oofem::MazarsMaterial:

Collaboration diagram for oofem::MazarsMaterial:

Public Member Functions
	MazarsMaterial (int n, Domain *d)
	Constructor. More...

virtual	~MazarsMaterial ()
	Destructor. More...

virtual const char *	giveInputRecordName () const

virtual const char *	giveClassName () const

virtual IRResultType	initializeFrom (InputRecord *ir)
	Initializes receiver according to object description stored in input record. More...

virtual void	computeEquivalentStrain (double &kappa, const FloatArray &strain, GaussPoint gp, TimeStep tStep)
	Computes the equivalent strain measure from given strain vector (full form). More...

virtual void	computeDamageParam (double &omega, double kappa, const FloatArray &strain, GaussPoint *gp)
	Computes the value of damage parameter omega, based on given value of equivalent strain. More...

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

Public Member Functions inherited from oofem::IsotropicDamageMaterial1
	IsotropicDamageMaterial1 (int n, Domain *d)
	Constructor. More...

virtual	~IsotropicDamageMaterial1 ()
	Destructor. More...

virtual void	giveInputRecord (DynamicInputRecord &input)
	Setups the input record string of receiver. More...

bool	isCrackBandApproachUsed ()

virtual void	computeEta (FloatArray &answer, const FloatArray &strain, GaussPoint gp, TimeStep tStep)
	Computes derivative of the equivalent strain with regards to strain. More...

void	computeDamageParamForCohesiveCrack (double &omega, double kappa, GaussPoint *gp)
	computes the value of damage parameter omega, based on a given value of equivalent strain, using iterations to achieve objectivity, based on the crack band concept (effective element size used) More...

double	damageFunction (double kappa, GaussPoint *gp)
	Returns the value of damage parameter corresponding to a given value of the damage-driving variable kappa, depending on the type of selected damage law, using a simple dependence (no adjustment for element size). More...

double	damageFunctionPrime (double kappa, GaussPoint *gp)
	Returns the value of compliance parameter corresponding to a given value of the damage-driving variable kappa, depending on the type of selected damage law, using a simple dependence (no adjustment for element size). More...

double	complianceFunction (double kappa, GaussPoint *gp)
	Returns the value of compliance parameter corresponding to a given value of the damage-driving variable kappa, depending on the type of selected damage law, using a simple dependence (no adjustment for element size). More...

double	evaluatePermanentStrain (double kappa, double omega)

virtual Interface *	giveInterface (InterfaceType it)
	Interface requesting service. More...

virtual int	MMI_map (GaussPoint gp, Domain oldd, TimeStep *tStep)
	Maps the required internal state variables from old mesh oldd to given ip. More...

virtual int	MMI_update (GaussPoint gp, TimeStep tStep, FloatArray *estrain=NULL)
	Updates the required internal state variables from previously mapped values. More...

virtual int	MMI_finish (TimeStep *tStep)
	Finishes the mapping for given time step. More...

virtual MaterialStatus *	giveStatus (GaussPoint *gp) const
	Returns material status of receiver in given integration point. More...

virtual double	give (int aProperty, GaussPoint *gp)
	Returns the value of material property 'aProperty'. More...

virtual bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode)
	Returns true if stiffness matrix of receiver is symmetric Default implementation returns true. More...

Public Member Functions inherited from oofem::IsotropicDamageMaterial
	IsotropicDamageMaterial (int n, Domain *d)
	Constructor. More...

virtual	~IsotropicDamageMaterial ()
	Destructor. More...

virtual int	hasNonLinearBehaviour ()
	Returns nonzero if receiver is non linear. More...

virtual int	hasMaterialModeCapability (MaterialMode mode)
	Tests if material supports material mode. More...

LinearElasticMaterial *	giveLinearElasticMaterial ()
	Returns reference to undamaged (bulk) material. More...

virtual void	give3dMaterialStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point. More...

virtual void	giveRealStressVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep)
	Computes the real stress vector for given total strain and integration point. More...

virtual void	giveRealStressVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. More...

virtual void	giveRealStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)
	Default implementation relies on giveRealStressVector_3d. More...

virtual void	giveRealStressVector_StressControl (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, const IntArray &strainControl, TimeStep tStep)
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero· More...

virtual void	giveRealStressVector_PlaneStress (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)
	Default implementation relies on giveRealStressVector_StressControl. More...

virtual void	giveRealStressVector_1d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)
	Default implementation relies on giveRealStressVector_StressControl. More...

virtual int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)
	Returns the integration point corresponding value in Reduced form. More...

virtual void	giveThermalDilatationVector (FloatArray &answer, GaussPoint , TimeStep )
	Returns a vector of coefficients of thermal dilatation in direction of each material principal (local) axis. More...

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

Public Member Functions inherited from oofem::StructuralMaterial
	StructuralMaterial (int n, Domain *d)
	Constructor. More...

virtual	~StructuralMaterial ()
	Destructor. More...

virtual void	giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Computes the stiffness matrix for giveRealStressVector of receiver in given integration point, respecting its history. More...

virtual void	giveRealStressVector_ShellStressControl (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, const IntArray &strainControl, TimeStep tStep)

virtual void	giveRealStressVector_Warping (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)
	Default implementation relies on giveRealStressVector_StressControl. More...

virtual void	giveRealStressVector_2dBeamLayer (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)
	Default implementation relies on giveRealStressVector_StressControl. More...

virtual void	giveRealStressVector_PlateLayer (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)
	Default implementation relies on giveRealStressVector_StressControl. More...

virtual void	giveRealStressVector_Fiber (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)
	Default implementation relies on giveRealStressVector_StressControl. More...

virtual void	giveRealStressVector_Lattice2d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)

virtual void	giveRealStressVector_Lattice3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)

virtual void	giveRealStressVector_2dPlateSubSoil (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)
	Default implementation is not provided. More...

virtual void	giveRealStressVector_3dBeamSubSoil (FloatArray &answer, GaussPoint gp, const FloatArray &reducedE, TimeStep tStep)

virtual void	giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
	Prototype for computation of Eshelby stress. More...

void	give_dPdF_from (const FloatMatrix &dSdE, FloatMatrix &answer, GaussPoint *gp)

void	convert_dSdE_2_dPdF (FloatMatrix &answer, const FloatMatrix &dSdE, const FloatArray &S, const FloatArray &F, MaterialMode matMode)

double	giveReferenceTemperature ()
	Returns the reference temperature of receiver. More...

virtual void	computeStressIndependentStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode)
	Computes reduced strain vector in given integration point, generated by internal processes in material, which are independent on loading in particular integration point. More...

virtual void	computeStressIndependentStrainVector_3d (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode)

virtual void	give3dMaterialStiffnessMatrix_dPdF (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

virtual void	give3dMaterialStiffnessMatrix_dCde (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

void	giveStressDependentPartOfStrainVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrainVector, TimeStep tStep, ValueModeType mode)
	Method for subtracting from reduced space strain vector its stress-independent parts (caused by temperature, shrinkage, creep and possibly by other phenomena). More...

void	giveStressDependentPartOfStrainVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrainVector, TimeStep tStep, ValueModeType mode)

virtual int	setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type)
	Sets the value of a certain variable at a given integration point to the given value. More...

virtual void	give2dBeamLayerStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 2d beam layer stiffness matrix of receiver. More...

virtual void	givePlateLayerStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 2d plate layer stiffness matrix of receiver. More...

virtual void	giveFiberStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 1d fiber stiffness matrix of receiver. More...

virtual void	give2dLatticeStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 2d lattice stiffness matrix of receiver. More...

virtual void	give3dLatticeStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 3d lattice stiffness matrix of receiver. More...

virtual void	give2dPlateSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing stiffness matrix of plate subsoil model. More...

virtual void	give3dBeamSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing stiffness matrix of beam3d subsoil model. More...

virtual void	giveFirstPKStressVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. More...

virtual void	giveFirstPKStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
	Default implementation relies on giveFirstPKStressVector_3d. More...

virtual void	giveFirstPKStressVector_PlaneStress (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
	Default implementation relies on giveFirstPKStressVector_3d. More...

virtual void	giveFirstPKStressVector_1d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
	Default implementation relies on giveFirstPKStressVector_3d. More...

virtual void	giveCauchyStressVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)

virtual void	giveCauchyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)

virtual void	giveCauchyStressVector_PlaneStress (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)

virtual void	giveCauchyStressVector_1d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)

virtual void	givePlaneStressStiffMtrx_dPdF (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)

virtual void	givePlaneStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)

virtual void	givePlaneStrainStiffMtrx_dPdF (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)

virtual void	givePlaneStrainStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)

virtual void	give1dStressStiffMtrx_dPdF (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)

virtual void	give1dStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)

Public Member Functions inherited from oofem::Material
	Material (int n, Domain *d)
	Constructor. More...

virtual	~Material ()
	Destructor. More...

virtual bool	hasProperty (int aProperty, GaussPoint *gp)
	Returns true if 'aProperty' exists on material. More...

virtual void	modifyProperty (int aProperty, double value, GaussPoint *gp)
	Modify 'aProperty', which already exists on material. More...

double	giveCastingTime ()

virtual bool	isActivated (TimeStep *tStep)

virtual int	hasCastingTimeSupport ()
	Tests if material supports casting time. More...

virtual void	printYourself ()
	Prints receiver state on stdout. Useful for debugging. More...

virtual contextIOResultType	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
	Stores integration point state to output stream. More...

virtual contextIOResultType	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
	Reads integration point state to output stream. More...

virtual int	checkConsistency ()
	Allows programmer to test some internal data, before computation begins. More...

virtual int	initMaterial (Element *element)
	Optional function to call specific procedures when initializing a material. More...

virtual int	packUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
	Pack all necessary data of integration point (according to element parallel_mode) into given communication buffer. More...

virtual int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
	Unpack and updates all necessary data of given integration point (according to element parallel_mode) into given communication buffer. More...

virtual int	estimatePackSize (DataStream &buff, GaussPoint *ip)
	Estimates the necessary pack size to hold all packed data of receiver. More...

virtual double	predictRelativeComputationalCost (GaussPoint *gp)
	Returns the weight representing relative computational cost of receiver The reference material model is linear isotropic material - its weight is set to 1.0 The other material models should compare to this reference model. More...

virtual double	predictRelativeRedistributionCost (GaussPoint *gp)
	Returns the relative redistribution cost of the receiver. More...

virtual void	initTempStatus (GaussPoint *gp)
	Initializes temporary variables stored in integration point status at the beginning of new time step. More...

Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
	Regular constructor, creates component with given number and belonging to given domain. More...

virtual	~FEMComponent ()
	Virtual destructor. More...

Domain *	giveDomain () const

virtual void	setDomain (Domain *d)
	Sets associated Domain. More...

int	giveNumber () const

void	setNumber (int num)
	Sets number of receiver. More...

virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
	Local renumbering support. More...

virtual contextIOResultType	saveContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Stores receiver state to output stream. More...

virtual contextIOResultType	restoreContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Restores the receiver state previously written in stream. More...

virtual void	printOutputAt (FILE file, TimeStep tStep)
	Prints output of receiver to stream, for given time step. More...

std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros). More...

Public Member Functions inherited from oofem::RandomMaterialExtensionInterface
	RandomMaterialExtensionInterface ()
	Constructor. More...

virtual	~RandomMaterialExtensionInterface ()
	Destructor. More...

IRResultType	initializeFrom (InputRecord *ir)
	Initializes receiver according to object description stored in input record. More...

void	giveInputRecord (DynamicInputRecord &ir)

bool	give (int key, GaussPoint *gp, double &value)
	Returns the property in associated status of given integration point if defined. More...

Public Member Functions inherited from oofem::Interface
	Interface ()
	Constructor. More...

virtual	~Interface ()

Public Member Functions inherited from oofem::MaterialModelMapperInterface
	MaterialModelMapperInterface ()
	Constructor. More...

virtual	~MaterialModelMapperInterface ()
	Destructor. More...

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 }
	Type characterizing the algorithm used to compute equivalent strain measure. More...

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 }
	Type characterizing the formula for the damage law. More...

Protected Types inherited from oofem::IsotropicDamageMaterial
enum	loaUnloCriterium { idm_strainLevelCR, idm_damageLevelCR }
	Variable controlling type of loading/unloading law, default set to idm_strainLevel defines the two two possibilities: More...

Protected Member Functions
virtual void	initDamaged (double kappa, FloatArray &totalStrainVector, GaussPoint *gp)
	Perfoms initialization, when damage first appear. More...

int	giveNumberOfSpatialDimensions (GaussPoint *gp)

void	giveNormalBlockOfElasticCompliance (FloatMatrix &answer, GaussPoint *gp)

double	computeGt (double kappa, GaussPoint *gp)

double	computeGc (double kappa, GaussPoint *gp)

Protected Member Functions inherited from oofem::IsotropicDamageMaterial
virtual void	givePlaneStressStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing plane stress stiffness matrix of receiver. More...

virtual void	givePlaneStrainStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing plane strain stiffness matrix of receiver. More...

virtual void	give1dStressStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 1d stiffness matrix of receiver. More...

Protected Member Functions inherited from oofem::RandomMaterialExtensionInterface
void	_generateStatusVariables (GaussPoint *) const
	Sets up (generates) the variables identified in randVariables array using generators given in randomVariableGenerators and stores them in given status. More...

Protected Attributes
double	E
	Elastic parameters. More...

double	nu

double	At
	Model parameters related to the shape of uniaxial stress-strain diagrams. More...

double	Bt

double	Ac

double	Bc

double	hReft
	Reference elem-length for objectivity. More...

double	hRefc

double	beta
	Beta coefficient reducing the effect of shear; default val = 1.06. More...

enum oofem::MazarsMaterial::mazarsModelVariant	modelVersion

Protected Attributes inherited from oofem::IsotropicDamageMaterial1
double	e0
	Equivalent strain at stress peak (or a similar parameter). More...

double	ef
	Determines ductility -> corresponds to fracturing strain. More...

double	wf
	Determines ductility -> corresponds to crack opening in the cohesive crack model. More...

double	gf
	Determines the softening -> corresponds to the initial fracture energy. More...

double	gft
	Determines the softening for the bilinear law -> corresponds to the total fracture energy. More...

double	ek
	Determines the softening for the bilinear law -> corresponds to the strain at the knee point. More...

double	wk
	Determines the softening for the bilinear law -> corresponds to the crack opening at the knee point. More...

double	sk
	Determines the softening for the bilinear law -> corresponds to the stress at the knee point. More...

double	c1
	Parameters used in Hordijk's softening law. More...

double	c2

EquivStrainType	equivStrainType
	Parameter specifying the definition of equivalent strain. More...

double	k
	Parameter used in Mises definition of equivalent strain. More...

double	griff_n
	Parameter used in Griffith's criterion. More...

int	damageLaw
	Temporary parameter reading type of softening law, used in other isotropic damage material models. More...

SofteningType	softType
	Parameter specifying the type of softening (damage law). More...

double	At
	Parameters used in Mazars damage law. More...

double	Bt

double	md
	Parameter used in "smooth damage law". More...

double	e1
	Parameters used if softType = 7 (extended smooth damage law) More...

double	e2

double	s1

double	nd

int	checkSnapBack
	Check possible snap back flag. More...

double	ep
	auxiliary input variablesfor softType == ST_SmoothExtended More...

double	ft

double	ps_alpha
	Parameters used by the model with permanent strain. More...

double	ps_H

ElementCharSizeMethod	ecsMethod
	Method used for evaluation of characteristic element size. More...

Set *	sourceElemSet
	Cached source element set used to map internal variables (adaptivity), created on demand. More...

Protected Attributes inherited from oofem::IsotropicDamageMaterial
double	tempDillatCoeff
	Coefficient of thermal dilatation. More...

double	maxOmega
	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. More...

int	permStrain
	Indicator of the type of permanent strain formulation (0 = standard damage with no permanent strain) More...

LinearElasticMaterial *	linearElasticMaterial
	Reference to bulk (undamaged) material. More...

enum oofem::IsotropicDamageMaterial::loaUnloCriterium	llcriteria

Protected Attributes inherited from oofem::StructuralMaterial
double	referenceTemperature
	Reference temperature (temperature, when material has been built into structure). More...

Protected Attributes inherited from oofem::Material
Dictionary	propertyDictionary
	Property dictionary. More...

double	castingTime
	Casting time. More...

Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number. More...

Domain *	domain
	Link to domain object, useful for communicating with other FEM components. More...

Protected Attributes inherited from oofem::RandomMaterialExtensionInterface
IntArray	randVariables
	Array of randomized variables (identified by a key). More...

IntArray	randomVariableGenerators
	Array of generators id's for corresponding randomized variables. More...

Additional Inherited Members
Static Public Member Functions inherited from oofem::IsotropicDamageMaterial1
static void	computeStrainInvariants (const FloatArray &strainVector, double &I1e, double &J2e)
	Computes invariants I1 and J2 of the strain tensor from the strain components stored in a vector. More...

Static Public Member Functions inherited from oofem::StructuralMaterial
static int	giveSymVI (int ind1, int ind2)

static int	giveVI (int ind1, int ind2)

static int	giveVoigtVectorMask (IntArray &answer, MaterialMode mmode)
	Returns a mask of the vector indicies corresponding to components in a general (non-symmetric) second order tensor of some stress/strain/deformation measure that performes work. More...

static int	giveVoigtSymVectorMask (IntArray &answer, MaterialMode mmode)
	The same as giveVoigtVectorMask but returns a mask corresponding to a symmetric second order tensor. More...

static void	giveInvertedVoigtVectorMask (IntArray &answer, MaterialMode mmode)
	Gives the inverted version of giveVoigtVectorMask. More...

static int	giveSizeOfVoigtVector (MaterialMode mmode)
	Returns the size of reduced stress/strain vector according to given mode. More...

static int	giveSizeOfVoigtSymVector (MaterialMode mmode)
	Returns the size of symmetric part of a reduced stress/strain vector according to given mode. More...

static void	giveFullVectorForm (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
	Converts the reduced symmetric Voigt vector (2nd order tensor) to full form. More...

static void	giveFullVectorFormF (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
	Converts the reduced deformation gradient Voigt vector (2nd order tensor). More...

static void	giveFullSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the reduced unsymmetric Voigt vector (2nd order tensor) to full form. More...

static void	giveReducedVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the full symmetric Voigt vector (2nd order tensor) to reduced form. More...

static void	giveReducedSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the full unsymmetric Voigt vector (2nd order tensor) to reduced form. More...

static void	giveFullSymMatrixForm (FloatMatrix &answer, const FloatMatrix &red, MaterialMode matMode)
	Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form. More...

static void	giveReducedMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
	Converts the full symmetric Voigt matrix (4th order tensor) to reduced form. More...

static void	giveReducedSymMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
	Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form. More...

static void	transformStrainVectorTo (FloatArray &answer, const FloatMatrix &base, const FloatArray &strainVector, bool transpose=false)
	Transforms 3d strain vector into another coordinate system. More...

static void	transformStressVectorTo (FloatArray &answer, const FloatMatrix &base, const FloatArray &stressVector, bool transpose=false)
	Transforms 3d stress vector into another coordinate system. More...

static double	computeVonMisesStress (const FloatArray *currentStress)
	Computes equivalent of von Mises stress. More...

static void	giveStrainVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false)
	Computes 3d strain vector transformation matrix from standard vector transformation matrix. More...

static void	give2DStrainVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false)
	Computes 2d strain vector transformation matrix from standard vector transformation matrix. More...

static void	giveStressVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false)
	Computes 3d stress vector transformation matrix from standard vector transformation matrix. More...

static void	givePlaneStressVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false)
	Computes 2d stress vector transformation matrix from standard vector transformation matrix. More...

static void	sortPrincDirAndValCloseTo (FloatArray pVal, FloatMatrix pDir, FloatMatrix *toPDir)
	Method for sorting newly computed principal values (pVal) and corresponding principal directions (pDir) to be closed to some (often previous) principal directions (toPDir). More...

static void	computePrincipalValues (FloatArray &answer, const FloatArray &s, stressStrainPrincMode mode)
	Common functions for convenience. More...

static void	computePrincipalValDir (FloatArray &answer, FloatMatrix &dir, const FloatArray &s, stressStrainPrincMode mode)
	Computes principal values and directions of stress or strain vector. More...

static double	computeDeviatoricVolumetricSplit (FloatArray &dev, const FloatArray &s)
	Computes split of receiver into deviatoric and volumetric part. More...

static void	computeDeviatoricVolumetricSum (FloatArray &s, const FloatArray &dev, double mean)

static void	applyDeviatoricElasticCompliance (FloatArray &strain, const FloatArray &stress, double EModulus, double nu)

static void	applyDeviatoricElasticCompliance (FloatArray &strain, const FloatArray &stress, double GModulus)

static void	applyDeviatoricElasticStiffness (FloatArray &stress, const FloatArray &strain, double EModulus, double nu)

static void	applyDeviatoricElasticStiffness (FloatArray &stress, const FloatArray &strain, double GModulus)

static void	applyElasticStiffness (FloatArray &stress, const FloatArray &strain, double EModulus, double nu)

static void	applyElasticCompliance (FloatArray &strain, const FloatArray &stress, double EModulus, double nu)

static double	computeStressNorm (const FloatArray &stress)

static double	computeFirstInvariant (const FloatArray &s)

static double	computeSecondStressInvariant (const FloatArray &s)

static double	computeThirdStressInvariant (const FloatArray &s)

static double	computeFirstCoordinate (const FloatArray &s)

static double	computeSecondCoordinate (const FloatArray &s)

static double	computeThirdCoordinate (const FloatArray &s)

Static Public Attributes inherited from oofem::StructuralMaterial
static std::vector< std::vector< int > >	vIindex
	Voigt index map. More...

static std::vector< std::vector< int > >	svIndex
	Symmetric Voigt index map. More...

Static Protected Attributes inherited from oofem::IsotropicDamageMaterial1
static MMAContainingElementProjection	mapper
	Mapper used to map internal variables in adaptivity. More...

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 95 of file mazarsmodel.h.

Member Enumeration Documentation

enum oofem::MazarsMaterial::mazarsModelVariant

protected

Model variants.

Enumerator
maz_original
maz_modTension

Definition at line 108 of file mazarsmodel.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 51 of file mazarsmodel.C.

References oofem::IsotropicDamageMaterial::idm_damageLevelCR, and oofem::IsotropicDamageMaterial::llcriteria.

oofem::MazarsMaterial::~MazarsMaterial ( )

virtual

Destructor.

Definition at line 61 of file mazarsmodel.C.

Member Function Documentation

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

virtual

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

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

virtual

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 130 of file mazarsmodel.C.

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

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

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

protected

Definition at line 412 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().

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

protected

Definition at line 364 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().

virtual MaterialStatus* oofem::MazarsMaterial::CreateStatus ( GaussPoint * gp ) const

inlinevirtual

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 125 of file mazarsmodel.h.

References oofem::FEMComponent::domain, and oofem::MazarsMaterialStatus::MazarsMaterialStatus().

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

inlinevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::IsotropicDamageMaterial1.

Reimplemented in oofem::MazarsNLMaterial.

Definition at line 118 of file mazarsmodel.h.

References oofem::IntegrationPointStatus::gp, oofem::MaterialStatus::initializeFrom(), and oofem::IsotropicDamageMaterialStatus::kappa.

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

inlinevirtual

Returns: Input record name of the receiver.

Reimplemented from oofem::IsotropicDamageMaterial1.

Definition at line 117 of file mazarsmodel.h.

References _IFT_MazarsMaterial_Name.

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

protected

Definition at line 215 of file mazarsmodel.C.

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

Referenced by computeDamageParam().

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

protected

Definition at line 203 of file mazarsmodel.C.

References oofem::GaussPoint::giveMaterialMode().

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

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

protectedvirtual

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

IRResultType oofem::MazarsMaterial::initializeFrom ( InputRecord * ir )

virtual

Initializes receiver according to object description stored in input record.

This function is called immediately after creating object using constructor. Input record can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record.