oofem::AnisotropicDamageMaterial Class Reference

#include <anisodamagemodel.h>

Inheritance diagram for oofem::AnisotropicDamageMaterial:

Collaboration diagram for oofem::AnisotropicDamageMaterial:

Public Member Functions
	AnisotropicDamageMaterial (int n, Domain *d)
	Constructor.
bool	hasMaterialModeCapability (MaterialMode mode) const override
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
FloatArrayF< 3 >	giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const override
	Plane-stress version of the stress evaluation algorithm.
void	computePrincValDir2D (double &D1, double &D2, double &c, double &s, double Dx, double Dy, double Dxy) const
bool	checkPrincVal2D (double Dx, double Dy, double Dxy) const
void	computeDamage (FloatMatrix &tempDamage, const FloatMatrix &damage, double kappa, double eps1, double eps2, double ceps, double seps, double epsZ) const
double	computeTraceD (double equivStrain) const
double	computeOutOfPlaneStrain (const FloatArray &inplaneStrain, const FloatMatrix &dam, bool tens_flag) const
double	computeDimensionlessOutOfPlaneStress (const FloatArray &inplaneStrain, double epsZ, const FloatMatrix &dam) const
void	computeInplaneStress (FloatArray &inplaneStress, const FloatArray &inplaneStrain, double epsZ, const FloatMatrix &dam) const
double	obtainAlpha1 (FloatMatrix tempDamageTensor, double deltaLambda, FloatMatrix positiveStrainTensor, double damageThreshold) const
	Obtains the proportion of the damage tensor that is needed to get the first eigenvalue equal to the damage threshold.
double	obtainAlpha2 (FloatMatrix tempDamageTensor, double deltaLambda, FloatMatrix positiveStrainTensor, FloatMatrix ProjMatrix, double damageThreshold) const
	Obtains the proportion of the damage tensor that is needed to get the second eigenvalue equal to the damage threshold.
double	obtainAlpha3 (FloatMatrix tempDamageTensor, double deltaLambda, FloatMatrix positiveStrainTensor, FloatArray vec3, double damageThreshold) const
	Obtains the proportion of the damage tensor that is needed to get the third eigenvalue equal to the damage threshold.
double	checkSymmetry (FloatMatrix matrix) const
void	correctBigValues (FloatMatrix &matrix) const
double	computeTraceD (FloatMatrix tempDamageTensor, FloatMatrix strainTensor, GaussPoint *gp) const
double	computeCorrectionFactor (FloatMatrix tempDamageTensor, FloatMatrix strainTensor, GaussPoint *gp) const
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< 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 *atTime) override
void	computeEquivalentStrain (double &kappa, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep) const
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
void	computeDamageTensor (FloatMatrix &answer, GaussPoint *gp, const FloatArray &totalStrain, double equivStrain, TimeStep *atTime) const
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_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)
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< 4, 4 >	givePlaneStrainStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
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 bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode) const
virtual double	giveCharacteristicValue (MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const
	Returns characteristic value of the receiver.
virtual double	give (int aProperty, GaussPoint *gp) const
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.
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 Types
enum	EquivStrainType {   EST_Unknown , EST_Mazars , EST_Rankine_Smooth , EST_Rankine_Standard ,   EST_ElasticEnergy , EST_ElasticEnergyPositiveStress , EST_ElasticEnergyPositiveStrain , EST_Mises ,   EST_Griffith }
	Type characterizing the algorithm used to compute equivalent strain measure. More...
enum	DamLawType {   DLT_Unknown , DLT_Desmorat1 , DLT_Desmorat2 , DLT_Linear ,   DLT_Exponential }
	Type characterizing the damage law. More...

Protected Member Functions
FloatMatrixF< 3, 3 >	givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 1, 1 >	give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
void	computePlaneStressStrain (FloatMatrix &answer, FloatMatrix damageTensor, FloatArray totalStrain, GaussPoint *gp, TimeStep *atTime) const
void	computePlaneStressSigmaZ (double &answer, FloatMatrix damageTensor, FloatArray reducedTotalStrainVector, double epsZ, GaussPoint *gp, TimeStep *atTime) const
virtual void	computeSecantOperator (FloatMatrix &answer, FloatMatrix strainTensor, FloatMatrix damageTensor, GaussPoint *gp) const
double	computeK (GaussPoint *gp)
double	computeKappa (FloatMatrix damageTensor) const

Protected Attributes
IsotropicLinearElasticMaterial	linearElasticMaterial
	Reference to bulk (undamaged) material.
double	E = 0.
	Young's modulus.
double	nu = 0.
	Poisson's ratio.
double	kappa0 = 0.
	Damage threshold kappa0, as defined in the paper mentioned above.
double	kappaf = 0.
	Damage parameter kappa_f (in the paper denoted as "a").
double	aA = 0.
	Damage parameter a*A, needed to obtain Kappa(trD), according to eq. 33 in the paper mentioned above.
EquivStrainType	equivStrainType = EST_Unknown
	Parameter specifying the definition of equivalent strain.
DamLawType	damageLawType = DLT_Unknown
	Parameter specifying the damage law.
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

Class representing anisotropic damage model. Based on the paper : "Nonlocal anisotropic damage model and related computational aspects for quasi-brittle materials" by R. Desmorat, F. Gatuingt, and F. Ragueneau Plane stress case implemented by the algorithm described in a report by M. Jirasek & F. Suarez, 25 April 2014

Author

: Martin Horak : nitra.nosp@m.mkar.nosp@m.oh@se.nosp@m.znam.nosp@m..cz

: Fernando Suarez : ferna.nosp@m.ndo..nosp@m.suare.nosp@m.z@fs.nosp@m.v.cvu.nosp@m.t.cz

: Milan Jirasek : milan.nosp@m..jir.nosp@m.asek@.nosp@m.fsv..nosp@m.cvut..nosp@m.cz

Definition at line 177 of file anisodamagemodel.h.

Member Enumeration Documentation

DamLawType

enum oofem::AnisotropicDamageMaterial::DamLawType

protected

Type characterizing the damage law.

Enumerator
DLT_Unknown
DLT_Desmorat1
DLT_Desmorat2
DLT_Linear
DLT_Exponential

Definition at line 197 of file anisodamagemodel.h.

EquivStrainType

enum oofem::AnisotropicDamageMaterial::EquivStrainType

protected

Type characterizing the algorithm used to compute equivalent strain measure.

Enumerator
EST_Unknown
EST_Mazars
EST_Rankine_Smooth
EST_Rankine_Standard
EST_ElasticEnergy
EST_ElasticEnergyPositiveStress
EST_ElasticEnergyPositiveStrain
EST_Mises
EST_Griffith

Definition at line 193 of file anisodamagemodel.h.

Constructor & Destructor Documentation

AnisotropicDamageMaterial()

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

Constructor.

Definition at line 53 of file anisodamagemodel.C.

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

Referenced by giveRealStressVector_1d(), giveRealStressVector_3d(), giveRealStressVector_PlaneStrain(), and giveRealStressVector_StressControl().

Member Function Documentation

checkPrincVal2D()

bool oofem::AnisotropicDamageMaterial::checkPrincVal2D	(	double	Dx,
		double	Dy,
		double	Dxy ) const

Definition at line 214 of file anisodamagemodel.C.

Referenced by computeDamage().

checkSymmetry()

double oofem::AnisotropicDamageMaterial::checkSymmetry

(

FloatMatrix

matrix

)

const

Definition at line 1028 of file anisodamagemodel.C.

References oofem::FloatMatrix::at(), and oofem::FloatMatrix::giveNumberOfRows().

computeCorrectionFactor()

double oofem::AnisotropicDamageMaterial::computeCorrectionFactor	(	FloatMatrix	tempDamageTensor,
		FloatMatrix	strainTensor,
		GaussPoint *	gp ) const

Definition at line 1096 of file anisodamagemodel.C.

References oofem::FloatMatrix::at(), oofem::AnisotropicDamageMaterialStatus::giveFlag(), and oofem::Material::giveStatus().

Referenced by giveRealStressVector().

computeDamage()

void oofem::AnisotropicDamageMaterial::computeDamage	(	FloatMatrix &	tempDamage,
		const FloatMatrix &	damage,
		double	kappa,
		double	eps1,
		double	eps2,
		double	ceps,
		double	seps,
		double	epsZ ) const

Definition at line 229 of file anisodamagemodel.C.

References checkPrincVal2D(), computePrincValDir2D(), computeTraceD(), and oofem::macbra().

Referenced by giveRealStressVector_PlaneStress().

computeDamageTensor()

void oofem::AnisotropicDamageMaterial::computeDamageTensor	(	FloatMatrix &	answer,
		GaussPoint *	gp,
		const FloatArray &	totalStrain,
		double	equivStrain,
		TimeStep *	atTime ) const

Definition at line 1419 of file anisodamagemodel.C.

References oofem::FloatMatrix::add(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beProductOf(), computeKappa(), computePlaneStressStrain(), computeTraceD(), correctBigValues(), oofem::AnisotropicDamageMaterialStatus::giveDamage(), oofem::Material::giveStatus(), oofem::AnisotropicDamageMaterialStatus::giveTempStrainZ(), oofem::FloatMatrix::jaco_(), obtainAlpha1(), obtainAlpha2(), obtainAlpha3(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::subtract(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

Referenced by giveRealStressVector().

computeDimensionlessOutOfPlaneStress()

double oofem::AnisotropicDamageMaterial::computeDimensionlessOutOfPlaneStress	(	const FloatArray &	inplaneStrain,
		double	epsZ,
		const FloatMatrix &	dam ) const

Definition at line 373 of file anisodamagemodel.C.

References computePrincValDir2D(), oofem::macbra(), and nu.

Referenced by giveRealStressVector_PlaneStress().

computeEquivalentStrain()

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

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.

Definition at line 650 of file anisodamagemodel.C.

References oofem::FloatArray::at(), oofem::StructuralMaterial::computePrincipalValues(), equivStrainType, EST_Mazars, oofem::GaussPoint::giveMaterialMode(), oofem::Material::giveStatus(), oofem::FloatArray::isEmpty(), nu, OOFEM_ERROR, and oofem::principal_strain.

Referenced by giveRealStressVector().

computeInplaneStress()

void oofem::AnisotropicDamageMaterial::computeInplaneStress	(	FloatArray &	inplaneStress,
		const FloatArray &	inplaneStrain,
		double	epsZ,
		const FloatMatrix &	dam ) const

Definition at line 410 of file anisodamagemodel.C.

References computePrincValDir2D(), E, oofem::macbra(), and nu.

Referenced by giveRealStressVector_PlaneStress().

computeK()

double oofem::AnisotropicDamageMaterial::computeK ( GaussPoint * gp )

protected

computeKappa()

double oofem::AnisotropicDamageMaterial::computeKappa ( FloatMatrix damageTensor ) const

protected

Definition at line 794 of file anisodamagemodel.C.

References oofem::FloatMatrix::giveTrace(), kappa0, kappaf, and oofem::trace().

Referenced by computeDamageTensor(), and giveRealStressVector().

computeOutOfPlaneStrain()

double oofem::AnisotropicDamageMaterial::computeOutOfPlaneStrain	(	const FloatArray &	inplaneStrain,
		const FloatMatrix &	dam,
		bool	tens_flag ) const

Definition at line 337 of file anisodamagemodel.C.

References computePrincValDir2D(), oofem::macbra(), and nu.

Referenced by giveRealStressVector_PlaneStress().

computePlaneStressSigmaZ()

void oofem::AnisotropicDamageMaterial::computePlaneStressSigmaZ ( double & answer, FloatMatrix damageTensor, FloatArray reducedTotalStrainVector, double epsZ, GaussPoint * gp, TimeStep * atTime ) const

protected

Definition at line 1341 of file anisodamagemodel.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatArray::beProductOf(), E, oofem::FloatMatrix::jaco_(), nu, oofem::FloatArray::resize(), and oofem::FloatMatrix::resize().

computePlaneStressStrain()

void oofem::AnisotropicDamageMaterial::computePlaneStressStrain ( FloatMatrix & answer, FloatMatrix damageTensor, FloatArray totalStrain, GaussPoint * gp, TimeStep * atTime ) const

protected

Definition at line 1258 of file anisodamagemodel.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatArray::beProductOf(), oofem::FloatMatrix::jaco_(), nu, oofem::FloatArray::resize(), and oofem::FloatMatrix::resize().

Referenced by computeDamageTensor(), and givePlaneStressStiffMtrx().

computePrincValDir2D()

void oofem::AnisotropicDamageMaterial::computePrincValDir2D	(	double &	D1,
		double &	D2,
		double &	c,
		double &	s,
		double	Dx,
		double	Dy,
		double	Dxy ) const

Definition at line 185 of file anisodamagemodel.C.

Referenced by computeDamage(), computeDimensionlessOutOfPlaneStress(), computeInplaneStress(), computeOutOfPlaneStrain(), and giveRealStressVector_PlaneStress().

computeSecantOperator()

void oofem::AnisotropicDamageMaterial::computeSecantOperator ( FloatMatrix & answer, FloatMatrix strainTensor, FloatMatrix damageTensor, GaussPoint * gp ) const

protectedvirtual

Definition at line 1730 of file anisodamagemodel.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeTraceD(), E, oofem::FloatMatrix::jaco_(), nu, oofem::FloatMatrix::resize(), oofem::FloatMatrix::subtract(), and oofem::FloatMatrix::zero().

computeTraceD() [1/2]

double oofem::AnisotropicDamageMaterial::computeTraceD

(

double

equivStrain

)

const

Definition at line 308 of file anisodamagemodel.C.

References aA, damageLawType, DLT_Desmorat1, DLT_Desmorat2, DLT_Exponential, DLT_Linear, kappa0, kappaf, nu, and OOFEM_ERROR.

Referenced by computeDamage(), computeDamageTensor(), computeSecantOperator(), and giveRealStressVector().

computeTraceD() [2/2]

double oofem::AnisotropicDamageMaterial::computeTraceD	(	FloatMatrix	tempDamageTensor,
		FloatMatrix	strainTensor,
		GaussPoint *	gp ) const

Definition at line 1063 of file anisodamagemodel.C.

References oofem::FloatMatrix::at(), oofem::AnisotropicDamageMaterialStatus::giveFlag(), and oofem::Material::giveStatus().

correctBigValues()

void oofem::AnisotropicDamageMaterial::correctBigValues

(

FloatMatrix &

matrix

)

const

Definition at line 1048 of file anisodamagemodel.C.

References oofem::FloatMatrix::at(), and oofem::FloatMatrix::giveNumberOfRows().

Referenced by computeDamageTensor(), giveRealStressVector(), obtainAlpha1(), obtainAlpha2(), and obtainAlpha3().

CreateStatus()

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

Definition at line 284 of file anisodamagemodel.h.

give1dStressStiffMtrx()

FloatMatrixF< 1, 1 > oofem::AnisotropicDamageMaterial::give1dStressStiffMtrx ( MatResponseMode mmode, GaussPoint * gp, TimeStep * tStep ) const

overrideprotectedvirtual

Method for computing 1d stiffness matrix of receiver. Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to 1d 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

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

Returns

Stiffness matrix.

Reimplemented from oofem::StructuralMaterial.

Definition at line 1240 of file anisodamagemodel.C.

References oofem::FloatArray::at(), oofem::Material::giveStatus(), linearElasticMaterial, and OOFEM_ERROR.

give3dMaterialStiffnessMatrix()

FloatMatrixF< 6, 6 > oofem::AnisotropicDamageMaterial::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 1136 of file anisodamagemodel.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStressDependentPartOfStrainVector(), and linearElasticMaterial.

giveClassName()

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

inlineoverridevirtual

Returns

Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 207 of file anisodamagemodel.h.

giveInputRecord()

void oofem::AnisotropicDamageMaterial::giveInputRecord ( DynamicInputRecord & input )

overridevirtual

Setups the input record string of receiver.

Parameters

input

Dynamic input record to be filled by receiver.

Reimplemented from oofem::FEMComponent.

Definition at line 2056 of file anisodamagemodel.C.

References _IFT_AnisotropicDamageMaterial_kappa0, kappa0, and oofem::DynamicInputRecord::setField().

giveInputRecordName()

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

inlineoverridevirtual

Returns

Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 208 of file anisodamagemodel.h.

References _IFT_AnisotropicDamageMaterial_Name.

giveIPValue()

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

Definition at line 1935 of file anisodamagemodel.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::AnisotropicDamageMaterialStatus::giveDamage(), oofem::AnisotropicDamageMaterialStatus::giveDissWork(), oofem::AnisotropicDamageMaterialStatus::giveKappa(), oofem::Material::giveStatus(), oofem::AnisotropicDamageMaterialStatus::giveStressWork(), oofem::AnisotropicDamageMaterialStatus::giveTempDamage(), oofem::FloatMatrix::jaco_(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

givePlaneStressStiffMtrx()

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

overrideprotectedvirtual

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 1177 of file anisodamagemodel.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrixF< N, M >::at(), computePlaneStressStrain(), E, oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStressDependentPartOfStrainVector(), linearElasticMaterial, nu, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

giveRealStressVector()

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

check if computeTraceD is necessary and check its implementation

Todo

Move this to StructuralCrossSection ?

Reimplemented from oofem::StructuralMaterial.

Definition at line 488 of file anisodamagemodel.C.

References oofem::FloatMatrix::add(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beMatrixForm(), oofem::FloatArray::beProductOf(), oofem::FloatMatrix::beProductOf(), oofem::FloatArray::beSymVectorForm(), computeCorrectionFactor(), computeDamageTensor(), computeEquivalentStrain(), computeKappa(), computeTraceD(), oofem::AnisotropicDamageMaterialStatus::computeWork(), correctBigValues(), oofem::AnisotropicDamageMaterialStatus::giveDamage(), oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStiffnessMatrix(), oofem::StructuralMaterial::giveStressDependentPartOfStrainVector(), oofem::Material::initTempStatus(), oofem::FloatMatrix::jaco_(), kappa0, oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), linearElasticMaterial, oofem::FloatMatrix::resize(), oofem::AnisotropicDamageMaterialStatus::setTempDamage(), oofem::AnisotropicDamageMaterialStatus::setTempKappa(), oofem::FloatMatrix::subtract(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

Referenced by giveRealStressVector_1d(), giveRealStressVector_3d(), giveRealStressVector_PlaneStrain(), and giveRealStressVector_StressControl().

giveRealStressVector_1d()

FloatArrayF< 1 > oofem::AnisotropicDamageMaterial::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 259 of file anisodamagemodel.h.

References AnisotropicDamageMaterial(), and giveRealStressVector().

giveRealStressVector_3d()

FloatArrayF< 6 > oofem::AnisotropicDamageMaterial::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 241 of file anisodamagemodel.h.

References AnisotropicDamageMaterial(), and giveRealStressVector().

giveRealStressVector_PlaneStrain()

FloatArrayF< 4 > oofem::AnisotropicDamageMaterial::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 247 of file anisodamagemodel.h.

References AnisotropicDamageMaterial(), and giveRealStressVector().

giveRealStressVector_PlaneStress()

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

overridevirtual

Plane-stress version of the stress evaluation algorithm.

Reimplemented from oofem::StructuralMaterial.

Definition at line 71 of file anisodamagemodel.C.

References AD_MAXITER, AD_TOLERANCE, oofem::FloatArray::at(), computeDamage(), computeDimensionlessOutOfPlaneStress(), computeInplaneStress(), computeOutOfPlaneStrain(), computePrincValDir2D(), oofem::AnisotropicDamageMaterialStatus::computeWork(), oofem::AnisotropicDamageMaterialStatus::giveDamage(), oofem::AnisotropicDamageMaterialStatus::giveKappa(), oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStressDependentPartOfStrainVector(), oofem::Material::initTempStatus(), kappa0, oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), oofem::macbra(), OOFEM_ERROR, oofem::AnisotropicDamageMaterialStatus::setTempDamage(), oofem::AnisotropicDamageMaterialStatus::setTempKappa(), and oofem::AnisotropicDamageMaterialStatus::setTempStrainZ().

giveRealStressVector_StressControl()

FloatArray oofem::AnisotropicDamageMaterial::giveRealStressVector_StressControl ( const FloatArray & reducedE, const IntArray & strainControl, GaussPoint * gp, TimeStep * tStep ) const

inlineoverridevirtual

Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·

Todo

We need a relative tolerance here! A relative tolerance like this could work, but if a really small increment is performed it won't work (it will be limited by machine precision)

Reimplemented from oofem::StructuralMaterial.

Definition at line 253 of file anisodamagemodel.h.

References AnisotropicDamageMaterial(), and giveRealStressVector().

hasMaterialModeCapability()

bool oofem::AnisotropicDamageMaterial::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 59 of file anisodamagemodel.C.

initializeFrom()

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

Definition at line 2002 of file anisodamagemodel.C.

References _IFT_AnisotropicDamageMaterial_aA, _IFT_AnisotropicDamageMaterial_damageLawType, _IFT_AnisotropicDamageMaterial_kappa0, _IFT_AnisotropicDamageMaterial_kappaf, aA, damageLawType, DLT_Desmorat1, DLT_Desmorat2, DLT_Exponential, DLT_Linear, E, equivStrainType, EST_ElasticEnergy, EST_ElasticEnergyPositiveStrain, EST_ElasticEnergyPositiveStress, EST_Griffith, EST_Mazars, EST_Mises, EST_Rankine_Smooth, EST_Rankine_Standard, IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, kappa0, kappaf, linearElasticMaterial, and nu.

obtainAlpha1()

double oofem::AnisotropicDamageMaterial::obtainAlpha1	(	FloatMatrix	tempDamageTensor,
		double	deltaLambda,
		FloatMatrix	positiveStrainTensor,
		double	damageThreshold ) const

Obtains the proportion of the damage tensor that is needed to get the first eigenvalue equal to the damage threshold.

Definition at line 806 of file anisodamagemodel.C.

References oofem::FloatMatrix::add(), oofem::FloatArray::at(), oofem::FloatMatrix::beProductOf(), correctBigValues(), oofem::FloatArray::giveSize(), oofem::FloatMatrix::jaco_(), and oofem::FloatMatrix::times().

Referenced by computeDamageTensor().

obtainAlpha2()

double oofem::AnisotropicDamageMaterial::obtainAlpha2	(	FloatMatrix	tempDamageTensor,
		double	deltaLambda,
		FloatMatrix	positiveStrainTensor,
		FloatMatrix	ProjMatrix,
		double	damageThreshold ) const

Obtains the proportion of the damage tensor that is needed to get the second eigenvalue equal to the damage threshold.

Definition at line 886 of file anisodamagemodel.C.

References oofem::FloatMatrix::add(), oofem::FloatArray::at(), correctBigValues(), oofem::FloatArray::giveSize(), oofem::FloatMatrix::jaco_(), and oofem::FloatMatrix::times().

Referenced by computeDamageTensor().

obtainAlpha3()

double oofem::AnisotropicDamageMaterial::obtainAlpha3	(	FloatMatrix	tempDamageTensor,
		double	deltaLambda,
		FloatMatrix	positiveStrainTensor,
		FloatArray	vec3,
		double	damageThreshold ) const

Obtains the proportion of the damage tensor that is needed to get the third eigenvalue equal to the damage threshold.

Definition at line 969 of file anisodamagemodel.C.

References oofem::FloatMatrix::add(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beDyadicProductOf(), oofem::FloatMatrix::beProductOf(), correctBigValues(), oofem::FloatMatrix::jaco_(), and oofem::FloatMatrix::times().

Referenced by computeDamageTensor().

Member Data Documentation

aA

double oofem::AnisotropicDamageMaterial::aA = 0.

protected

Damage parameter a*A, needed to obtain Kappa(trD), according to eq. 33 in the paper mentioned above.

Definition at line 191 of file anisodamagemodel.h.

Referenced by computeTraceD(), and initializeFrom().

damageLawType

DamLawType oofem::AnisotropicDamageMaterial::damageLawType = DLT_Unknown

protected

Parameter specifying the damage law.

Definition at line 199 of file anisodamagemodel.h.

Referenced by computeTraceD(), and initializeFrom().

E

double oofem::AnisotropicDamageMaterial::E = 0.

protected

Young's modulus.

Definition at line 183 of file anisodamagemodel.h.

Referenced by computeInplaneStress(), computePlaneStressSigmaZ(), computeSecantOperator(), givePlaneStressStiffMtrx(), and initializeFrom().

equivStrainType

EquivStrainType oofem::AnisotropicDamageMaterial::equivStrainType = EST_Unknown

protected

Parameter specifying the definition of equivalent strain.

Definition at line 195 of file anisodamagemodel.h.

Referenced by computeEquivalentStrain(), and initializeFrom().

kappa0

double oofem::AnisotropicDamageMaterial::kappa0 = 0.

protected

Damage threshold kappa0, as defined in the paper mentioned above.

Definition at line 187 of file anisodamagemodel.h.

Referenced by computeKappa(), computeTraceD(), giveInputRecord(), giveRealStressVector(), giveRealStressVector_PlaneStress(), and initializeFrom().

kappaf

double oofem::AnisotropicDamageMaterial::kappaf = 0.

protected

Damage parameter kappa_f (in the paper denoted as "a").

Definition at line 189 of file anisodamagemodel.h.

Referenced by computeKappa(), computeTraceD(), and initializeFrom().

linearElasticMaterial

IsotropicLinearElasticMaterial oofem::AnisotropicDamageMaterial::linearElasticMaterial

protected

Reference to bulk (undamaged) material.

Definition at line 181 of file anisodamagemodel.h.

Referenced by AnisotropicDamageMaterial(), give1dStressStiffMtrx(), give3dMaterialStiffnessMatrix(), givePlaneStressStiffMtrx(), giveRealStressVector(), and initializeFrom().

nu

double oofem::AnisotropicDamageMaterial::nu = 0.

protected

Poisson's ratio.

Definition at line 185 of file anisodamagemodel.h.

Referenced by computeDimensionlessOutOfPlaneStress(), computeEquivalentStrain(), computeInplaneStress(), computeOutOfPlaneStrain(), computePlaneStressSigmaZ(), computePlaneStressStrain(), computeSecantOperator(), computeTraceD(), givePlaneStressStiffMtrx(), and initializeFrom().

---

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

• anisodamagemodel.h
• anisodamagemodel.C