OOFEM 3.0
Loading...
Searching...
No Matches
oofem::IsotropicDamageMaterial1 Class Reference

#include <idm1.h>

Inheritance diagram for oofem::IsotropicDamageMaterial1:
Collaboration diagram for oofem::IsotropicDamageMaterial1:

Public Member Functions

 IsotropicDamageMaterial1 (int n, Domain *d)
 Constructor.
virtual ~IsotropicDamageMaterial1 ()
 Destructor.
const char * giveClassName () const override
const char * giveInputRecordName () const override
void initializeFrom (InputRecord &ir) override
void giveInputRecord (DynamicInputRecord &input) override
bool isCrackBandApproachUsed () const
double computeEquivalentStrain (const FloatArray &strain, GaussPoint *gp, TimeStep *tStep) const override
void computeEta (FloatArray &answer, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep) const override
double computeDamageParam (double kappa, const FloatArray &strain, GaussPoint *gp) 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
InterfacegiveInterface (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
std::unique_ptr< MaterialStatusCreateStatus (GaussPoint *gp) const override
MaterialStatusgiveStatus (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
LinearElasticMaterialgiveLinearElasticMaterial () 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.
DomaingiveDomain () 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.

Static Public Member Functions

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)

Protected Types

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
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 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.
SetsourceElemSet = 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).
LinearElasticMateriallinearElasticMaterial = 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.
Domaindomain
 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.

Static Protected Attributes

static MMAContainingElementProjection mapper
 Mapper used to map internal variables in adaptivity.

Additional Inherited Members

Static Public Attributes inherited from oofem::StructuralMaterial
static std::array< std::array< int, 3 >, 3 > vIindex
 Voigt index map.
static std::array< std::array< int, 3 >, 3 > svIndex
 Symmetric Voigt index map.

Detailed Description

This class implements a simple local isotropic damage model for concrete in tension. A model is based on isotropic damage concept, assuming that damage evolution law is postulated in explicit form, relation damage parameter (omega) to scalar measure of the largest strain level ever reached in material (kappa).

Definition at line 139 of file idm1.h.

Member Enumeration Documentation

◆ EquivStrainType

Type characterizing the algorithm used to compute equivalent strain measure. Note that the assigned numbers to enum values have to correspond to values used in initializeFrom to resolve EquivStrainType. If not, the consistency between initializeFrom and giveInputRecord methods is lost.

Enumerator
EST_Mazars 
EST_Rankine_Smooth 
EST_ElasticEnergy 
EST_Mises 
EST_Rankine_Standard 
EST_ElasticEnergyPositiveStress 
EST_ElasticEnergyPositiveStrain 
EST_Griffith 
EST_Unknown 

Definition at line 182 of file idm1.h.

◆ SofteningType

Type characterizing the formula for the damage law. For example, linear softening can be specified with fracturing strain or crack opening.

Enumerator
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 

Definition at line 208 of file idm1.h.

Constructor & Destructor Documentation

◆ IsotropicDamageMaterial1()

◆ ~IsotropicDamageMaterial1()

oofem::IsotropicDamageMaterial1::~IsotropicDamageMaterial1 ( )
virtual

Destructor.

Definition at line 78 of file idm1.C.

References sourceElemSet.

Member Function Documentation

◆ complianceFunction()

double oofem::IsotropicDamageMaterial1::complianceFunction ( double kappa,
GaussPoint * gp ) const

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). The compliance parameter gamma is defined as gamma = omega/(1-omega) where omega is the damage.

Parameters
kappaEquivalent strain measure.
gpIntegration point.

Definition at line 1134 of file idm1.C.

References damageFunction().

Referenced by oofem::IDNLMaterial::updateBeforeNonlocAverage().

◆ computeDamageParam()

double oofem::IsotropicDamageMaterial1::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]omegaContains result.
kappaEquivalent strain measure.
strainTotal strain in full form.
gpIntegration point.

Implements oofem::IsotropicDamageMaterial.

Reimplemented in oofem::MazarsMaterial.

Definition at line 796 of file idm1.C.

References computeDamageParamForCohesiveCrack(), damageFunction(), isCrackBandApproachUsed(), softType, and ST_Disable_Damage.

Referenced by oofem::IsotropicGradientDamageMaterial::giveRealStressVectorGradientDamage().

◆ computeDamageParamForCohesiveCrack()

double oofem::IsotropicDamageMaterial1::computeDamageParamForCohesiveCrack ( double kappa,
GaussPoint * gp ) const

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)

Parameters
[out]omegaContains the resulting damage.
kappaEquivalent strain measure.
gpIntegration point.

Definition at line 808 of file idm1.C.

References checkSnapBack, E, e0, e0_ID, ef, ek, ek_ID, f_k, f_r, gf, gf_ID, gft, gft_ID, give(), oofem::GaussPoint::giveElement(), oofem::Element::giveGlobalNumber(), oofem::IsotropicDamageMaterialStatus::giveLe(), oofem::FEMComponent::giveNumber(), giveStatus(), IDM1_ITERATION_LIMIT, oofem::IsotropicDamageMaterial::linearElasticMaterial, oofem::IsotropicDamageMaterial::maxOmega, OOFEM_ERROR, OOFEM_WARNING, sk, softType, ST_BiLinear_Cohesive_Crack, ST_Exponential_Cohesive_Crack, ST_Linear_Cohesive_Crack, ST_Trilinear_Cohesive_Crack, w_f, w_k, w_r, wf, wf_ID, and wk.

Referenced by computeDamageParam().

◆ computeEquivalentStrain()

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

◆ computeEta()

void oofem::IsotropicDamageMaterial1::computeEta ( FloatArray & answer,
const FloatArray & strain,
GaussPoint * gp,
TimeStep * tStep ) const
overridevirtual

◆ computeStrainInvariants()

void oofem::IsotropicDamageMaterial1::computeStrainInvariants ( const FloatArray & strainVector,
double & I1e,
double & J2e )
static

Computes invariants I1 and J2 of the strain tensor from the strain components stored in a vector.

Parameters
strainVectorInput strain components.
[out]I1eOutput value of strain invariant I1.
[out]J2eOutput value of strain invariant J2.

Definition at line 786 of file idm1.C.

References oofem::FloatArray::at(), and s1.

Referenced by computeEquivalentStrain().

◆ CreateStatus()

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

Creates new copy of associated status and inserts it into given integration point.

Parameters
gpIntegration point where newly created status will be stored.
Returns
Reference to new status.

Reimplemented from oofem::IsotropicDamageMaterial.

Reimplemented in oofem::IsotropicGradientDamageMaterial, oofem::MazarsMaterial, and oofem::MazarsNLMaterial.

Definition at line 1386 of file idm1.C.

Referenced by giveStatus().

◆ damageFunction()

double oofem::IsotropicDamageMaterial1::damageFunction ( double kappa,
GaussPoint * gp ) const

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

Parameters
kappaEquivalent strain measure.
gpIntegration point.

Definition at line 954 of file idm1.C.

References At, Bt, c2, e0, e0_ID, e1, e2, ef, ef_ID, give(), md, nd, OOFEM_WARNING, oofem::IsotropicDamageMaterial::permStrain, ps_alpha, s1, softType, ST_DoubleExponential, ST_Exponential, ST_Linear, ST_Mazars, ST_ModPowerExponential, ST_PowerExponential, ST_Smooth, and ST_SmoothExtended.

Referenced by complianceFunction(), oofem::IDNLMaterial::computeDamageParam(), computeDamageParam(), and oofem::IDNLMaterial::updateBeforeNonlocAverage().

◆ damageFunctionPrime()

double oofem::IsotropicDamageMaterial1::damageFunctionPrime ( double kappa,
GaussPoint * gp ) const
overridevirtual

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). The compliance parameter gamma is defined as gamma = omega/(1-omega) where omega is the damage.

Parameters
kappaEquivalent strain measure.
gpIntegration point. Returns the value of derivative of damage function wrt damage-driving variable kappa corresponding to a given value of the kappa, depending on the type of selected damage law.
kappaEquivalent strain measure.
gpIntegration point.

Reimplemented from oofem::IsotropicDamageMaterial.

Definition at line 1024 of file idm1.C.

References At, Bt, E, e0, e0_ID, e1, e2, ef, ef_ID, f_k, f_r, gf, give(), oofem::IsotropicDamageMaterialStatus::giveLe(), giveStatus(), oofem::IsotropicDamageMaterialStatus::giveTempDamage(), oofem::IsotropicDamageMaterial::linearElasticMaterial, md, nd, OOFEM_ERROR, s1, softType, ST_Exponential, ST_Exponential_Cohesive_Crack, ST_Linear, ST_Linear_Cohesive_Crack, ST_Mazars, ST_Smooth, ST_SmoothExtended, ST_Trilinear_Cohesive_Crack, w_f, w_k, w_r, and wf.

Referenced by oofem::IsotropicGradientDamageMaterial::giveGradientDamageStiffnessMatrix_dd_BN(), oofem::IsotropicGradientDamageMaterial::giveGradientDamageStiffnessMatrix_dd_NN(), and oofem::IsotropicGradientDamageMaterial::giveGradientDamageStiffnessMatrix_ud().

◆ evaluatePermanentStrain()

double oofem::IsotropicDamageMaterial1::evaluatePermanentStrain ( double kappa,
double omega ) const
overridevirtual

Reimplemented from oofem::IsotropicDamageMaterial.

Definition at line 1141 of file idm1.C.

References oofem::IsotropicDamageMaterial::permStrain, and ps_alpha.

◆ give()

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

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

Parameters
aPropertyID of property requested.
gpIntegration point,
Returns
Property value.

Reimplemented from oofem::IsotropicDamageMaterial.

Definition at line 1324 of file idm1.C.

References e0, e0_ID, ef, ef_ID, ek, ek_ID, gf, gf_ID, gft, gft_ID, giveStatus(), wf, and wf_ID.

Referenced by computeDamageParamForCohesiveCrack(), damageFunction(), damageFunctionPrime(), oofem::IDNLMaterial::giveLocalNonlocalStiffnessContribution(), initDamaged(), and oofem::IDNLMaterial::NonlocalMaterialStiffnessInterface_showSparseMtrxStructure().

◆ giveClassName()

const char * oofem::IsotropicDamageMaterial1::giveClassName ( ) const
inlineoverridevirtual
Returns
Class name of the receiver.

Reimplemented from oofem::IsotropicDamageMaterial.

Reimplemented in oofem::IsotropicGradientDamageMaterial, oofem::MazarsMaterial, and oofem::MazarsNLMaterial.

Definition at line 258 of file idm1.h.

◆ giveInputRecord()

void oofem::IsotropicDamageMaterial1::giveInputRecord ( DynamicInputRecord & input)
overridevirtual

◆ giveInputRecordName()

const char * oofem::IsotropicDamageMaterial1::giveInputRecordName ( ) const
inlineoverridevirtual
Returns
Input record name of the receiver.

Implements oofem::FEMComponent.

Reimplemented in oofem::IsotropicGradientDamageMaterial, and oofem::MazarsMaterial.

Definition at line 259 of file idm1.h.

References _IFT_IsotropicDamageMaterial1_Name.

◆ giveInterface()

Interface * oofem::IsotropicDamageMaterial1::giveInterface ( InterfaceType t)
overridevirtual

Interface requesting service.

See also
InterfaceType
Returns
Requested interface if implemented, otherwise NULL.

Reimplemented from oofem::FEMComponent.

Reimplemented in oofem::IsotropicGradientDamageMaterial, and oofem::MazarsNLMaterial.

Definition at line 1348 of file idm1.C.

References oofem::MaterialModelMapperInterface::MaterialModelMapperInterface(), and oofem::MaterialModelMapperInterfaceType.

◆ giveStatus()

MaterialStatus * oofem::IsotropicDamageMaterial1::giveStatus ( GaussPoint * gp) const
overridevirtual

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

Parameters
gpReturns reference to material status belonging to integration point gp.
Returns
Material status associated with given integration point.

Reimplemented from oofem::Material.

Definition at line 1392 of file idm1.C.

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

Referenced by oofem::IDNLMaterial::computeAngleAndSigmaRatio(), computeDamageParamForCohesiveCrack(), oofem::IsotropicGradientDamageMaterial::computeEikonalInternalLength_a(), oofem::IsotropicGradientDamageMaterial::computeEikonalInternalLength_aPrime(), oofem::IsotropicGradientDamageMaterial::computeEikonalInternalLength_b(), oofem::IsotropicGradientDamageMaterial::computeEikonalInternalLength_bPrime(), oofem::IDNLMaterial::computeEquivalentStrain(), oofem::MazarsNLMaterial::computeEquivalentStrain(), oofem::MazarsMaterial::computeGc(), oofem::MazarsMaterial::computeGt(), oofem::IsotropicGradientDamageMaterial::computeInternalLength(), damageFunctionPrime(), give(), oofem::IsotropicGradientDamageMaterial::giveGradientDamageStiffnessMatrix_dd_BN(), oofem::IsotropicGradientDamageMaterial::giveGradientDamageStiffnessMatrix_dd_NN(), oofem::IsotropicGradientDamageMaterial::giveGradientDamageStiffnessMatrix_du(), oofem::IsotropicGradientDamageMaterial::giveGradientDamageStiffnessMatrix_ud(), oofem::IsotropicGradientDamageMaterial::giveGradientDamageStiffnessMatrix_uu(), oofem::IDNLMaterial::giveIPValue(), oofem::IDNLMaterial::giveLocalNonlocalStiffnessContribution(), oofem::IsotropicGradientDamageMaterial::giveNonlocalInternalForces_N_factor(), oofem::IDNLMaterial::giveNonlocalMetricModifierAt(), oofem::IsotropicGradientDamageMaterial::giveRealStressVectorGradientDamage(), oofem::IDNLMaterial::giveRemoteNonlocalStiffnessContribution(), initDamaged(), oofem::MazarsMaterial::initDamaged(), oofem::MazarsNLMaterial::initDamaged(), MMI_map(), MMI_update(), oofem::IDNLMaterial::NonlocalMaterialStiffnessInterface_addIPContribution(), oofem::IDNLMaterial::NonlocalMaterialStiffnessInterface_giveIntegrationDomainList(), oofem::IDNLMaterial::NonlocalMaterialStiffnessInterface_showSparseMtrxStructure(), oofem::IDNLMaterial::packUnknowns(), oofem::MazarsNLMaterial::packUnknowns(), oofem::IDNLMaterial::predictRelativeComputationalCost(), oofem::IDNLMaterial::unpackAndUpdateUnknowns(), oofem::MazarsNLMaterial::unpackAndUpdateUnknowns(), oofem::IDNLMaterial::updateBeforeNonlocAverage(), and oofem::MazarsNLMaterial::updateBeforeNonlocAverage().

◆ initDamaged()

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

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

Parameters
kappaScalar measure of strain level.
totalStrainVectorCurrent total strain vector.
gpIntegration point.

Reimplemented from oofem::IsotropicDamageMaterial.

Reimplemented in oofem::MazarsMaterial, and oofem::MazarsNLMaterial.

Definition at line 1168 of file idm1.C.

References oofem::FloatMatrix::add(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatArray::beColumnOf(), oofem::FloatArray::beProductOf(), checkSnapBack, oofem::StructuralMaterial::computePrincipalValDir(), oofem::FEMComponent::domain, E, e0, e0_ID, ecsMethod, ef, ef_ID, equivStrainType, EST_Griffith, gf, gf_ID, give(), oofem::Element::giveCharacteristicSize(), oofem::GaussPoint::giveElement(), oofem::FloatArray::giveIndexMaxElem(), oofem::FloatArray::giveIndexMinElem(), oofem::Element::giveLabel(), oofem::GaussPoint::giveMaterialMode(), giveStatus(), oofem::StructuralMaterial::giveStiffnessMatrix(), isCrackBandApproachUsed(), oofem::IsotropicDamageMaterial::linearElasticMaterial, M_PI, OOFEM_ERROR, OOFEM_WARNING, oofem::principal_strain, oofem::principal_stress, oofem::FloatArray::rotatedWith(), sk, softType, ST_BiLinear_Cohesive_Crack, ST_Disable_Damage, ST_Exponential_Cohesive_Crack, ST_Linear_Cohesive_Crack, ST_Trilinear_Cohesive_Crack, w_f, wf, wf_ID, wk, and oofem::FloatMatrix::zero().

Referenced by oofem::IsotropicGradientDamageMaterial::giveRealStressVectorGradientDamage().

◆ initializeFrom()

void oofem::IsotropicDamageMaterial1::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
irInput record to initialize from.
priorityPriority of the input record. This is used to determine the order of initialization

Reimplemented from oofem::IsotropicDamageMaterial.

Reimplemented in oofem::IsotropicGradientDamageMaterial, oofem::MazarsMaterial, and oofem::MazarsNLMaterial.

Definition at line 86 of file idm1.C.

References _IFT_IsotropicDamageMaterial1_alphaps, _IFT_IsotropicDamageMaterial1_At, _IFT_IsotropicDamageMaterial1_Bt, _IFT_IsotropicDamageMaterial1_c1, _IFT_IsotropicDamageMaterial1_c2, _IFT_IsotropicDamageMaterial1_checkSnapBack, _IFT_IsotropicDamageMaterial1_damageLaw, _IFT_IsotropicDamageMaterial1_e0, _IFT_IsotropicDamageMaterial1_e1, _IFT_IsotropicDamageMaterial1_e2, _IFT_IsotropicDamageMaterial1_ecsm, _IFT_IsotropicDamageMaterial1_ef, _IFT_IsotropicDamageMaterial1_ek, _IFT_IsotropicDamageMaterial1_ep, _IFT_IsotropicDamageMaterial1_equivstraintype, _IFT_IsotropicDamageMaterial1_f_k, _IFT_IsotropicDamageMaterial1_f_r, _IFT_IsotropicDamageMaterial1_ft, _IFT_IsotropicDamageMaterial1_gf, _IFT_IsotropicDamageMaterial1_gft, _IFT_IsotropicDamageMaterial1_k, _IFT_IsotropicDamageMaterial1_md, _IFT_IsotropicDamageMaterial1_n, _IFT_IsotropicDamageMaterial1_nd, _IFT_IsotropicDamageMaterial1_sk, _IFT_IsotropicDamageMaterial1_skft, _IFT_IsotropicDamageMaterial1_w_f, _IFT_IsotropicDamageMaterial1_w_k, _IFT_IsotropicDamageMaterial1_w_r, _IFT_IsotropicDamageMaterial1_wf, _IFT_IsotropicDamageMaterial1_wk, _IFT_IsotropicDamageMaterial1_wkwf, _IFT_IsotropicLinearElasticMaterial_e, At, Bt, c1, c2, checkSnapBack, damageLaw, E, e0, e1, e2, oofem::ECSM_Oliver1, oofem::ECSM_Oliver1modified, oofem::ECSM_Projection, oofem::ECSM_ProjectionCentered, oofem::ECSM_SquareRootOfArea, ecsMethod, ef, ek, ep, equivStrainType, EST_ElasticEnergy, EST_ElasticEnergyPositiveStrain, EST_ElasticEnergyPositiveStress, EST_Griffith, EST_Mazars, EST_Mises, EST_Rankine_Smooth, EST_Rankine_Standard, f_k, f_r, ft, gf, gft, griff_n, oofem::InputRecord::hasField(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, k, oofem::IsotropicDamageMaterial::linearElasticMaterial, mapper, md, nd, OOFEM_WARNING, oofem::IsotropicDamageMaterial::permStrain, ps_alpha, s1, sk, softType, ST_BiLinear_Cohesive_Crack, ST_Disable_Damage, ST_DoubleExponential, ST_Exponential, ST_Exponential_Cohesive_Crack, ST_Hordijk_Cohesive_Crack, ST_Linear, ST_Linear_Cohesive_Crack, ST_Mazars, ST_ModPowerExponential, ST_PowerExponential, ST_Smooth, ST_SmoothExtended, ST_Trilinear_Cohesive_Crack, w_f, w_k, w_r, wf, and wk.

Referenced by restoreContext().

◆ isCharacteristicMtrxSymmetric()

bool oofem::IsotropicDamageMaterial1::isCharacteristicMtrxSymmetric ( MatResponseMode rMode) const
inlineoverridevirtual

Returns true if stiffness matrix of receiver is symmetric Default implementation returns true.

Reimplemented from oofem::Material.

Definition at line 344 of file idm1.h.

◆ isCrackBandApproachUsed()

bool oofem::IsotropicDamageMaterial1::isCrackBandApproachUsed ( ) const
inline

◆ MMI_finish()

int oofem::IsotropicDamageMaterial1::MMI_finish ( TimeStep * tStep)
overridevirtual

Finishes the mapping for given time step. Used to perform cleanup. Typically some mappers require to compute some global mesh data related to current step, which are valid for example to all IPs - so they are computed only once for all IPs, stored and they need to be deallocated. These mappers are typically class variables, but their finish is invoked by all members.

Implements oofem::MaterialModelMapperInterface.

Definition at line 1485 of file idm1.C.

References mapper.

◆ MMI_map()

int oofem::IsotropicDamageMaterial1::MMI_map ( GaussPoint * gp,
Domain * oldd,
TimeStep * tStep )
overridevirtual

◆ MMI_update()

int oofem::IsotropicDamageMaterial1::MMI_update ( GaussPoint * gp,
TimeStep * tStep,
FloatArray * elemGPVec = nullptr )
overridevirtual

Updates the required internal state variables from previously mapped values. The result is stored in gp status. This map and update splitting is necessary, for example for nonlocal models that local quantity to be averaged must be mapped in all integration points and then update can happen, because it may depend on nonlocal variable, which is computed from local values.

Parameters
gpIntegration point belonging to new domain which values will be mapped.
tStepTime step.
elemGPVecVector passed to MMI_update at material level, probably computed from primary unknowns (for structural elements this represent strain vector).
Returns
Nonzero if o.k.

Implements oofem::MaterialModelMapperInterface.

Definition at line 1466 of file idm1.C.

References oofem::IsotropicDamageMaterial::giveRealStressVector(), giveStatus(), oofem::StructuralMaterialStatus::giveStrainVector(), and oofem::GaussPoint::updateYourself().

◆ restoreContext()

void oofem::IsotropicDamageMaterial1::restoreContext ( DataStream & stream,
ContextMode mode )
overridevirtual

Restores the receiver state previously written in stream.

See also
saveContext
Parameters
streamInput stream.
modeDetermines amount of info available in stream (state, definition, ...).
Exceptions
throwsan ContextIOERR exception if error encountered.

Reimplemented from oofem::IsotropicDamageMaterial.

Definition at line 1372 of file idm1.C.

References oofem::CIO_IOERR, CM_Definition, initializeFrom(), oofem::DataStream::read(), and THROW_CIOERR.

◆ saveContext()

void oofem::IsotropicDamageMaterial1::saveContext ( DataStream & stream,
ContextMode mode )
overridevirtual

Stores receiver state to output stream.

Parameters
streamOutput stream.
modeDetermines amount of info required in stream (state, definition, ...).
Exceptions
throwsan ContextIOERR exception if error encountered.

Reimplemented from oofem::IsotropicDamageMaterial.

Definition at line 1358 of file idm1.C.

References oofem::CIO_IOERR, CM_Definition, giveInputRecord(), oofem::DynamicInputRecord::giveRecordInTXTFormat(), THROW_CIOERR, and oofem::DataStream::write().

Member Data Documentation

◆ At

double oofem::IsotropicDamageMaterial1::At = 0.
protected

Parameters used in Mazars damage law.

Definition at line 214 of file idm1.h.

Referenced by damageFunction(), damageFunctionPrime(), giveInputRecord(), and initializeFrom().

◆ Bt

double oofem::IsotropicDamageMaterial1::Bt = 0.
protected

Definition at line 214 of file idm1.h.

Referenced by damageFunction(), damageFunctionPrime(), giveInputRecord(), and initializeFrom().

◆ c1

double oofem::IsotropicDamageMaterial1::c1 = 3.
protected

Parameters used in Hordijk's softening law.

Definition at line 172 of file idm1.h.

Referenced by initializeFrom().

◆ c2

double oofem::IsotropicDamageMaterial1::c2 = 6.93
protected

Definition at line 172 of file idm1.h.

Referenced by damageFunction(), and initializeFrom().

◆ checkSnapBack

int oofem::IsotropicDamageMaterial1::checkSnapBack = 0
protected

Check possible snap back flag.

Definition at line 221 of file idm1.h.

Referenced by computeDamageParamForCohesiveCrack(), initDamaged(), and initializeFrom().

◆ damageLaw

int oofem::IsotropicDamageMaterial1::damageLaw = 0
protected

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

Definition at line 203 of file idm1.h.

Referenced by giveInputRecord(), and initializeFrom().

◆ e0

◆ e1

double oofem::IsotropicDamageMaterial1::e1 = 0.
protected

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

Definition at line 219 of file idm1.h.

Referenced by oofem::IDNLMaterial::computeAngleAndSigmaRatio(), damageFunction(), damageFunctionPrime(), giveInputRecord(), and initializeFrom().

◆ e2

double oofem::IsotropicDamageMaterial1::e2 = 0.
protected

◆ ecsMethod

ElementCharSizeMethod oofem::IsotropicDamageMaterial1::ecsMethod = ECSM_Unknown
protected

Method used for evaluation of characteristic element size.

Definition at line 230 of file idm1.h.

Referenced by giveInputRecord(), initDamaged(), and initializeFrom().

◆ ef

double oofem::IsotropicDamageMaterial1::ef = 0.
protected

◆ ek

double oofem::IsotropicDamageMaterial1::ek = 0.
protected

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

Definition at line 163 of file idm1.h.

Referenced by computeDamageParamForCohesiveCrack(), give(), giveInputRecord(), and initializeFrom().

◆ ep

double oofem::IsotropicDamageMaterial1::ep = 0.
protected

auxiliary input variablesfor softType == ST_SmoothExtended

Definition at line 224 of file idm1.h.

Referenced by giveInputRecord(), and initializeFrom().

◆ equivStrainType

EquivStrainType oofem::IsotropicDamageMaterial1::equivStrainType = EST_Unknown
protected

◆ f_k

double oofem::IsotropicDamageMaterial1::f_k = 0.
protected

◆ f_r

double oofem::IsotropicDamageMaterial1::f_r = 0.
protected

◆ ft

double oofem::IsotropicDamageMaterial1::ft = 0.
protected

Definition at line 224 of file idm1.h.

Referenced by giveInputRecord(), and initializeFrom().

◆ gf

double oofem::IsotropicDamageMaterial1::gf = 0.
protected

Determines the softening -> corresponds to the initial fracture energy. For a linear law, it is the area under the stress/strain curve. For an exponential law, it is the area bounded by the elastic range and a tangent to the softening part of the curve at the peak stress. For a bilinear law, gf corresponds to area bounded by elasticity and the first linear softening line projected to zero stress.

Definition at line 157 of file idm1.h.

Referenced by computeDamageParamForCohesiveCrack(), damageFunctionPrime(), give(), giveInputRecord(), initDamaged(), and initializeFrom().

◆ gft

double oofem::IsotropicDamageMaterial1::gft = 0.
protected

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

Definition at line 160 of file idm1.h.

Referenced by computeDamageParamForCohesiveCrack(), give(), giveInputRecord(), and initializeFrom().

◆ griff_n

double oofem::IsotropicDamageMaterial1::griff_n = 8.
protected

Parameter used in Griffith's criterion.

Definition at line 200 of file idm1.h.

Referenced by computeEquivalentStrain(), and initializeFrom().

◆ k

double oofem::IsotropicDamageMaterial1::k = 0.
protected

Parameter used in Mises definition of equivalent strain.

Definition at line 197 of file idm1.h.

Referenced by computeEquivalentStrain(), and initializeFrom().

◆ mapper

MMAContainingElementProjection oofem::IsotropicDamageMaterial1::mapper
staticprotected

Mapper used to map internal variables in adaptivity.

Definition at line 241 of file idm1.h.

Referenced by giveInputRecord(), initializeFrom(), oofem::IsotropicGradientDamageMaterial::initializeFrom(), oofem::MazarsMaterial::initializeFrom(), MMI_finish(), and MMI_map().

◆ md

double oofem::IsotropicDamageMaterial1::md = 1.
protected

Parameter used in "smooth damage law".

Definition at line 216 of file idm1.h.

Referenced by damageFunction(), damageFunctionPrime(), giveInputRecord(), and initializeFrom().

◆ nd

double oofem::IsotropicDamageMaterial1::nd = 0.
protected

Definition at line 219 of file idm1.h.

Referenced by damageFunction(), damageFunctionPrime(), giveInputRecord(), and initializeFrom().

◆ ps_alpha

double oofem::IsotropicDamageMaterial1::ps_alpha = 0.
protected

Parameters used by the model with permanent strain.

Definition at line 227 of file idm1.h.

Referenced by damageFunction(), evaluatePermanentStrain(), and initializeFrom().

◆ ps_H

double oofem::IsotropicDamageMaterial1::ps_H = 0.
protected

Definition at line 227 of file idm1.h.

◆ s1

double oofem::IsotropicDamageMaterial1::s1 = 0.
protected

◆ sk

double oofem::IsotropicDamageMaterial1::sk = 0.
protected

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

Definition at line 169 of file idm1.h.

Referenced by computeDamageParamForCohesiveCrack(), giveInputRecord(), initDamaged(), and initializeFrom().

◆ softType

SofteningType oofem::IsotropicDamageMaterial1::softType = ST_Unknown
protected

◆ sourceElemSet

Set* oofem::IsotropicDamageMaterial1::sourceElemSet = nullptr
protected

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

Definition at line 233 of file idm1.h.

Referenced by MMI_map(), and ~IsotropicDamageMaterial1().

◆ w_f

double oofem::IsotropicDamageMaterial1::w_f = 0.
protected

◆ w_k

double oofem::IsotropicDamageMaterial1::w_k = 0.
protected

Parameters used in Trilinear_Cohesive_Crack softening law.

Definition at line 175 of file idm1.h.

Referenced by computeDamageParamForCohesiveCrack(), damageFunctionPrime(), giveInputRecord(), and initializeFrom().

◆ w_r

double oofem::IsotropicDamageMaterial1::w_r = 0.
protected

◆ wf

double oofem::IsotropicDamageMaterial1::wf = 0.
protected

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

Definition at line 149 of file idm1.h.

Referenced by computeDamageParamForCohesiveCrack(), damageFunctionPrime(), give(), giveInputRecord(), initDamaged(), and initializeFrom().

◆ wk

double oofem::IsotropicDamageMaterial1::wk = 0.
protected

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

Definition at line 166 of file idm1.h.

Referenced by computeDamageParamForCohesiveCrack(), giveInputRecord(), initDamaged(), and initializeFrom().


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

This page is part of the OOFEM-3.0 documentation. Copyright Copyright (C) 1994-2025 Borek Patzak Bořek Patzák
Project e-mail: oofem@fsv.cvut.cz
Generated at for OOFEM by doxygen 1.15.0 written by Dimitri van Heesch, © 1997-2011