oofem::MPSDamMaterial Class Reference

#include <mpsdammat.h>

Inheritance diagram for oofem::MPSDamMaterial:

Collaboration diagram for oofem::MPSDamMaterial:

Public Member Functions
	MPSDamMaterial (int n, Domain *d)
bool	hasMaterialModeCapability (MaterialMode mode) const override
const char *	giveInputRecordName () const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
int	giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
void	giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) const override
double	givee0 (GaussPoint *gp) const
double	givegf (GaussPoint *gp) const
void	initDamaged (double kappa, FloatArray &totalStrainVector, GaussPoint *gp, TimeStep *tStep) const
void	initDamagedFib (GaussPoint *gp, TimeStep *tStep) const
double	computeDamage (double kappa, GaussPoint *gp) const
double	computeDamageForCohesiveCrack (double kappa, GaussPoint *gp) const
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
virtual double	computeTensileStrength (double equivalentTime) const
virtual double	computeFractureEnergy (double equivalentTime) const
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 3, 3 >	givePlaneStressStiffMtrx (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 4, 4 >	givePlaneStrainStiffMtrx (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 1, 1 >	give1dStressStiffMtrx (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const override
Public Member Functions inherited from oofem::MPSMaterial
	MPSMaterial (int n, Domain *d)
const char *	giveInputRecordName () const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
int	giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
void	giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) const override
void	giveShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override
double	computeCreepFunction (double t, double t_prime, GaussPoint *gp, TimeStep *tStep) const override
	Evaluation of the basic creep compliance function - can be used to compute elastic modulus in derived damage material.
double	giveEquivalentTime (GaussPoint *gp, TimeStep *tStep) const override
	By default returns equivalent time in the middle of the time step.
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
Public Member Functions inherited from oofem::KelvinChainSolidMaterial
	KelvinChainSolidMaterial (int n, Domain *d)
void	computeHiddenVars (GaussPoint *gp, TimeStep *tStep) const
Public Member Functions inherited from oofem::RheoChainMaterial
	RheoChainMaterial (int n, Domain *d)
virtual	~RheoChainMaterial ()
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.
FloatArrayF< 3 >	giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &strain, GaussPoint *gp, TimeStep *tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 1 >	giveRealStressVector_1d (const FloatArrayF< 1 > &strain, GaussPoint *gp, TimeStep *tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 2 >	giveRealStressVector_2dBeamLayer (const FloatArrayF< 2 > &strain, GaussPoint *gp, TimeStep *tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 5 >	giveRealStressVector_PlateLayer (const FloatArrayF< 5 > &strain, GaussPoint *gp, TimeStep *tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 6 >	giveThermalDilatationVector (GaussPoint *gp, TimeStep *tStep) const override
bool	hasMaterialModeCapability (MaterialMode mode) const override
bool	hasCastingTimeSupport () const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
int	giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
void	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp) override
void	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp) override
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 3, 3 >	givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 4, 4 >	givePlaneStrainStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 1, 1 >	give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
FloatArray	computeStressIndependentStrainVector (GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
double	giveAlphaOne () const
double	giveAlphaTwo () const
double	givePoissonsRatio () const
	Returns Poisson's ratio.
bool	isActivated (TimeStep *tStep) const override
	Extended meaning: returns true if the material is cast (target time > casting time) or the precasing time mat is defined.
Public Member Functions inherited from oofem::StructuralMaterial
	StructuralMaterial (int n, Domain *d)
bool	hasMaterialModeCapability (MaterialMode mode) const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
void	giveCharacteristicMatrix (FloatMatrix &answer, MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const override
	Returns characteristic matrix of the receiver.
void	giveCharacteristicVector (FloatArray &answer, FloatArray &flux, MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const override
	Returns characteristic vector of the receiver.
virtual void	giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArray	giveRealStressVector_StressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint *gp, TimeStep *tStep) const
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArray	giveRealStressVector_ShellStressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArrayF< 2 >	giveRealStressVector_Warping (const FloatArrayF< 2 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_Fiber (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_2dPlateSubSoil (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
	Default implementation is not provided.
virtual FloatArrayF< 6 >	giveRealStressVector_3dBeamSubSoil (const FloatArrayF< 6 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArrayF< 9 >	giveFirstPKStressVector_3d (const FloatArrayF< 9 > &vF, GaussPoint *gp, TimeStep *tStep) const
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
virtual FloatArrayF< 5 >	giveFirstPKStressVector_PlaneStrain (const FloatArrayF< 5 > &vF, GaussPoint *gp, TimeStep *tStep) const
	Default implementation relies on giveFirstPKStressVector_3d.
virtual FloatArray	giveFirstPKStressVector_StressControl (const FloatArray &reducedvF, const IntArray &FControl, GaussPoint *gp, TimeStep *tStep) const
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArrayF< 4 >	giveFirstPKStressVector_PlaneStress (const FloatArrayF< 4 > &vF, GaussPoint *gp, TimeStep *tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
virtual FloatArrayF< 1 >	giveFirstPKStressVector_1d (const FloatArrayF< 1 > &vF, GaussPoint *gp, TimeStep *tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
virtual void	giveCauchyStressVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void	giveCauchyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void	giveCauchyStressVector_PlaneStress (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void	giveCauchyStressVector_1d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void	giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
double	giveReferenceTemperature ()
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 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
void	printYourself () override
	Prints receiver state on stdout. Useful for debugging.
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	SofteningType { ST_Exponential_Cohesive_Crack , ST_Linear_Cohesive_Crack , ST_Disable_Damage }
Protected Types inherited from oofem::MPSMaterial
enum	coupledAnalysisType { Basic , MPS_full , MPS_humidity , MPS_temperature }

Protected Attributes
bool	timeDepFracturing = false
double	fib_s = 0.
double	fib_fcm28 = 0.
bool	isotropic = false
double	E = -1.
	dummy Young's modulus
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.
double	ft = 0.
	Equivalent strain at stress peak (or a similar parameter).
double	const_gf = 0.
int	checkSnapBack = 1
	Check possible snap back flag.
SofteningType	softType = ST_Exponential_Cohesive_Crack
	Parameter specifying the type of softening (damage law).
ElementCharSizeMethod	ecsMethod = ECSM_Projection
	Method used for evaluation of characteristic element size.
double	ft28 = 0.
	28-day value of tensile strength. Used only with "timedepfracturing"
double	gf28 = 0.
	28-day value of fracture energy. Used only with "timedepfracturing"
Protected Attributes inherited from oofem::MPSMaterial
double	t0 = 0.
	age when temperature or humidity starts to change
double	q1 = 0.
	compliances of the B3 model
double	q2 = 0.
double	q3 = 0.
double	q4 = 0.
double	lambda0 = 0.
	constant equal to one day in time units of analysis (eg. 86400 if the analysis runs in seconds)
enum oofem::MPSMaterial::coupledAnalysisType	CoupledAnalysis
double	EspringVal = 0.
double	kSh = 0.
	additional parameters for sorption isotherm (used to compute relative humidity from water content)
double	muS = 0.
	fluidity parameter used in viscosity evolution equation
double	k3 = 0.
double	kTm = 0.
	kTm replaces ln(h) on RHS of the differential equation describing evolution of MPS
double	kTc = 0.
	parameter reducing creep effects of thermal cycling (replaces kTm in such case)
double	khc = 0.
	parameter reducing creep effects of hygral cycling (0 = complete reduction, 1 = no reduction, default value)
double	roomTemperature = 0.
	parameter reducing creep effects of thermal cycling
double	QEtoR = 0.
	activation energies
double	QRtoR = 0.
double	QStoR = 0.
double	alphaE = 0.
	parameters that control the effect of humidity on rates of hydration, creep and microprestress relaxation
double	alphaR = 0.
double	alphaS = 0.
double	p = 0.
	exponent in the microprestress/viscosity governing equation
double	sh_a = 0.
	parameters for nonlinear shrinkage function
double	sh_hC = 0.
double	sh_n = 0.
bool	timeDependent_ksh = false
	flag activating equivalent-time-dependence of shrinkage coefficient
FloatArray	ksh_h
	nonlinear shrinkage function - user-defined piecewise-linear function
FloatArray	ksh_fh
double	eps_cas0 = 0.
	parameter for autogenous shrinkage according to fib MC 2010
double	b4_eps_au_infty = 0.
	parameters for autogenous shrinkage according to B4 model
double	b4_tau_au = 0.
double	b4_alpha = 0.
double	b4_r_t = 0.
double	stiffnessFactor = 0.
	scaling factor 1. for Pa, 1.e6 for MPa - only for empirical formulas - q1-q4 and ft and gf
double	temperScaleDifference = 0.
	0 for Kelvin, 273.15 for Celsius
int	hydrationTimescaleTF
int	autoShrinkageTF
double	tau_nano = 0.
Protected Attributes inherited from oofem::RheoChainMaterial
double	talpha = 0.
	thermal dilatation coeff.
int	nUnits = 0
	Number of (Maxwell or Kelvin) units in the rheologic chain.
double	relMatAge = 0.
	Physical age of the material at castingTime.
bool	lattice = false
double	nu = 0.
	Poisson's ratio (assumed to be constant, unaffected by creep).
double	alphaOne = 0.
	Parameters for the lattice model.
double	alphaTwo = 0.
double	EparValTime = -1.
	Time for which the partial moduli of individual units have been evaluated.
double	begOfTimeOfInterest = 0.
	Time from which the model should give a good approximation. Optional field. Default value is 0.1 [day].
double	endOfTimeOfInterest = 0.
	Time (age???) up to which the model should give a good approximation.
StructuralMaterial *	linearElasticMaterial = nullptr
	Associated linearElasticMaterial, with E = 1.
FloatArray	EparVal
	Partial moduli of individual units.
FloatArray	charTimes
	Characteristic times of individual units (relaxation or retardation times).
FloatArray	discreteTimeScale
	Times at which the errors are evaluated if the least-square method is used.
double	timeFactor = 0.
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.
Protected Member Functions inherited from oofem::MPSMaterial
void	predictParametersFrom (double, double, double, double)
void	computeCharTimes () override
	Evaluation of characteristic times.
FloatArray	computeCharCoefficients (double tPrime, GaussPoint *gp, TimeStep *tStep) const override
	Evaluation of characteristic moduli of the non-aging Kelvin chain.
double	giveEModulus (GaussPoint *gp, TimeStep *tStep) const override
	Evaluation of the incremental modulus.
double	computeSolidifiedVolume (GaussPoint *gp, TimeStep *tStep) const override
	Evaluation of the relative volume of the solidified material.
double	computeBetaMu (GaussPoint *gp, TimeStep *tStep, int Mu) const override
	factors for exponential algorithm
double	computeLambdaMu (GaussPoint *gp, TimeStep *tStep, int Mu) const override
double	computeFlowTermViscosity (GaussPoint *gp, TimeStep *tStep) const
	Evaluation of the flow term viscosity.
double	giveInitViscosity (TimeStep *tStep) const
	Returns initial value of the flow term viscosity.
void	giveEigenStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override
bool	hasIncrementalShrinkageFormulation () const override
void	computePointShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) const
	Evaluation of the shrinkageStrainVector - shrinkage is fully dependent on humidity rate in given GP.
void	computeFibAutogenousShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) const
	Evaluation of the autogenousShrinkageStrainVector according to fib MC 2010 - autogenous shrinkage is fully dependent on the equivalent age at given GP.
void	computeB4AutogenousShrinkageStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) const
	Evaluation of the autogenousShrinkageStrainVector according to Bazant's B4 model. In the model the evolution depends on temperature adjusted age, here on equivalent age (additional humidity influence).
void	computeAutogenousShrinkageDefinedByTF (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) const
	Evaluation of the autogenousShrinkageStrainVector given by an auxiliary time function (autoShrinkageTF parameter). The time scale in that time function corresponds to the equivalent material age.
double	giveHumidity (GaussPoint *gp, TimeStep *tStep, int option) const
double	giveHumidityNano (GaussPoint *gp, TimeStep *tStep, int option) const
double	giveTemperature (GaussPoint *gp, TimeStep *tStep, int option) const
double	computePsiR (GaussPoint *gp, TimeStep *tStep, int option) const
double	computePsiS (GaussPoint *gp, TimeStep *tStep) const
	Evaluation of the factor transforming real time to reduced time (effect on the evolution of microprestress).
double	computePsiE (GaussPoint *gp, TimeStep *tStep) const
	Evaluation of the factor transforming real time to equivalent time (effect on the solidified volume).
double	computeEquivalentTime (GaussPoint *gp, TimeStep *tStep, int option) const
Protected Member Functions inherited from oofem::RheoChainMaterial
const FloatArray &	giveDiscreteTimes () const
void	computeDiscreteRelaxationFunction (FloatArray &answer, const FloatArray &tSteps, double t0, double tr, GaussPoint *gp, TimeStep *tStep) const
void	giveUnitComplianceMatrix (FloatMatrix &answer, GaussPoint *gp, TimeStep *tStep) const
	Evaluation of elastic compliance matrix for unit Young's modulus.
void	giveUnitStiffnessMatrix (FloatMatrix &answer, GaussPoint *gp, TimeStep *tStep) const
	Evaluation of elastic stiffness matrix for unit Young's modulus.
virtual void	updateEparModuli (double tPrime, GaussPoint *gp, TimeStep *tStep) const
	Update of partial moduli of individual chain units.
double	giveEparModulus (int iChain) const
	Access to partial modulus of a given unit.
double	giveCharTime (int) const
	Access to the characteristic time of a given unit.
virtual double	giveCharTimeExponent (int i) const
	Exponent to be used with the char time of a given unit, usually = 1.0.
StructuralMaterial *	giveLinearElasticMaterial ()
	Access to the underlying linear elastic material with unit Young's modulus.
double	giveEndOfTimeOfInterest ()
	Access to the time up to which the response should be accurate.
void	computeTrueStressIndependentStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const
Static Protected Member Functions inherited from oofem::RheoChainMaterial
static FloatArray	generateLogTimeScale (double from, double to, int nsteps)

Detailed Description

This class extends the material model based on MPS theory (microprestress-solidification) for concrete creep and shrinkage by a simple isotropic damage model to take into account cracking in tension. Most of the methods are simply copied from "idm1" and "isodamagemodel" files. Since the damage is not the main feature of this model the methods were simplified to make the code brief. The following simplifications are assumed:

• only two softening laws are implemented: linear and exponential
• the equivalent strain measure uses Rankine (standard) criterion (max. principal strain)
• the damage parameters (tensile strength, fracture energy) can be either specified or they can be computed according to recommendations from fib Model Code 2010
• on contrary to idm1, the crack vector specifies the direction of crack opening and not the crack orientation
• snapback is checked by default
• element length is computed using ecsMethod = ECSM_Projection

Definition at line 180 of file mpsdammat.h.

Member Enumeration Documentation

SofteningType

enum oofem::MPSDamMaterial::SofteningType

protected

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

Enumerator
ST_Exponential_Cohesive_Crack
ST_Linear_Cohesive_Crack
ST_Disable_Damage

Definition at line 215 of file mpsdammat.h.

Constructor & Destructor Documentation

MPSDamMaterial()

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

Definition at line 193 of file mpsdammat.C.

References oofem::MPSMaterial::MPSMaterial().

Member Function Documentation

computeDamage()

double oofem::MPSDamMaterial::computeDamage	(	double	kappa,
		GaussPoint *	gp ) const

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

Parameters

kappa	Equivalent strain measure.
strain	Total strain in full form.
gp	Integration point.

Returns

Damage.

Definition at line 586 of file mpsdammat.C.

References computeDamageForCohesiveCrack(), softType, and ST_Disable_Damage.

Referenced by giveRealStressVector().

computeDamageForCohesiveCrack()

double oofem::MPSDamMaterial::computeDamageForCohesiveCrack	(	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

kappa	Equivalent strain measure.
gp	Integration point.

Returns

Damage.

Definition at line 596 of file mpsdammat.C.

References checkSnapBack, E, givee0(), oofem::GaussPoint::giveElement(), givegf(), oofem::Element::giveGlobalNumber(), oofem::FEMComponent::giveNumber(), oofem::Material::giveStatus(), MPSDAMMAT_ITERATION_LIMIT, OOFEM_ERROR, OOFEM_WARNING, softType, ST_Exponential_Cohesive_Crack, and ST_Linear_Cohesive_Crack.

Referenced by computeDamage().

computeFractureEnergy()

double oofem::MPSDamMaterial::computeFractureEnergy ( double equivalentTime ) const

virtual

Definition at line 505 of file mpsdammat.C.

References computeTensileStrength(), fib_fcm28, and gf28.

Referenced by initDamagedFib().

computeTensileStrength()

double oofem::MPSDamMaterial::computeTensileStrength ( double equivalentTime ) const

virtual

Definition at line 547 of file mpsdammat.C.

References fib_fcm28, fib_s, and ft28.

Referenced by computeFractureEnergy(), giveIPValue(), and initDamagedFib().

CreateStatus()

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

overridevirtual

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::KelvinChainSolidMaterial.

Definition at line 738 of file mpsdammat.C.

References oofem::RheoChainMaterial::nUnits.

give1dStressStiffMtrx()

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

overridevirtual

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 799 of file mpsdammat.C.

References oofem::Material::giveStatus(), isotropic, and oofem::min().

give3dMaterialStiffnessMatrix()

FloatMatrixF< 6, 6 > oofem::MPSDamMaterial::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 748 of file mpsdammat.C.

References oofem::Material::giveStatus(), isotropic, and oofem::min().

giveClassName()

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

inlineoverridevirtual

Returns

Class name of the receiver.

Reimplemented from oofem::KelvinChainSolidMaterial.

Definition at line 233 of file mpsdammat.h.

givee0()

double oofem::MPSDamMaterial::givee0

(

GaussPoint *

gp

)

const

Definition at line 481 of file mpsdammat.C.

References E, ft, oofem::MPSDamMaterialStatus::givee0(), oofem::Material::giveStatus(), and timeDepFracturing.

Referenced by computeDamageForCohesiveCrack(), giveRealStressVector(), and initDamaged().

givegf()

double oofem::MPSDamMaterial::givegf

(

GaussPoint *

gp

)

const

Definition at line 493 of file mpsdammat.C.

References const_gf, oofem::MPSDamMaterialStatus::givegf(), oofem::Material::giveStatus(), and timeDepFracturing.

Referenced by computeDamageForCohesiveCrack(), giveRealStressVector(), and initDamaged().

giveInputRecordName()

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

inlineoverridevirtual

Returns

Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 232 of file mpsdammat.h.

References _IFT_MPSDamMaterial_Name.

giveIPValue()

int oofem::MPSDamMaterial::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 816 of file mpsdammat.C.

References oofem::FloatArray::at(), computeTensileStrength(), ft, oofem::Material::giveStatus(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().

givePlaneStrainStiffMtrx()

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

overridevirtual

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

Parameters

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

Reimplemented from oofem::StructuralMaterial.

Definition at line 782 of file mpsdammat.C.

References oofem::Material::giveStatus(), isotropic, and oofem::min().

givePlaneStressStiffMtrx()

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

overridevirtual

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

Parameters

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

Reimplemented from oofem::StructuralMaterial.

Definition at line 765 of file mpsdammat.C.

References oofem::Material::giveStatus(), isotropic, and oofem::min().

giveRealStressVector()

void oofem::MPSDamMaterial::giveRealStressVector ( FloatArray & answer, GaussPoint * gp, const FloatArray & reducedStrain, TimeStep * tStep ) const

overridevirtual

Computes the real stress vector for given total strain and integration point. The total strain is defined as strain computed directly from displacement field at given time. The stress independent parts (temperature, eigenstrains) are subtracted in constitutive driver. The service should use previously reached equilibrium history variables. Also it should update temporary history variables in status according to newly reached state. The temporary history variables are moved into equilibrium ones after global structure equilibrium has been reached by iteration process.

Parameters

answer	Stress vector in reduced form. For large deformations it is treated as the second Piola-Kirchoff stress.
gp	Integration point.
reducedStrain	Strain vector in reduced form. For large deformations it is treated as the Green-Lagrange strain.
tStep	Current time step (most models are able to respond only when tStep is current time step).

Todo

Move this to StructuralCrossSection ?

Reimplemented from oofem::KelvinChainSolidMaterial.

Definition at line 277 of file mpsdammat.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeDamage(), oofem::StressVector::computePrincipalValDir(), oofem::StructuralMaterial::computePrincipalValues(), oofem::StressStrainBaseVector::convertFromFullForm(), E, oofem::MPSDamMaterialStatus::giveCharLength(), oofem::MPSDamMaterialStatus::giveDamage(), givee0(), givegf(), oofem::MPSDamMaterialStatus::giveKappa(), oofem::GaussPoint::giveMaterialMode(), oofem::StructuralMaterial::givePlaneStressVectorTranformationMtrx(), oofem::FloatArray::giveSize(), oofem::StructuralMaterial::giveSizeOfVoigtSymVector(), oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStressVectorTranformationMtrx(), oofem::MPSMaterialStatus::giveTempAutogenousShrinkageStrain(), oofem::MPSMaterialStatus::giveTempDryingShrinkageStrain(), oofem::RheoChainMaterialStatus::giveTempThermalStrain(), initDamaged(), initDamagedFib(), oofem::RheoChainMaterial::isActivated(), isotropic, oofem::MPSMaterial::lambda0, oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), oofem::MPSDamMaterialStatus::letTempViscoelasticStressVectorBe(), OOFEM_ERROR, oofem::principal_strain, oofem::FloatArray::resize(), oofem::FloatArray::resizeWithValues(), oofem::FloatArray::rotatedWith(), oofem::MPSDamMaterialStatus::setCrackWidth(), oofem::MPSDamMaterialStatus::setResidualTensileStrength(), oofem::MPSDamMaterialStatus::setTempDamage(), oofem::MPSDamMaterialStatus::setTempKappa(), softType, ST_Exponential_Cohesive_Crack, ST_Linear_Cohesive_Crack, timeDepFracturing, oofem::FloatArray::times(), and oofem::FloatArray::zero().

hasMaterialModeCapability()

bool oofem::MPSDamMaterial::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 198 of file mpsdammat.C.

initDamaged()

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

Abstract service allowing to perform some initialization, when damage first appear.

Parameters

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

Definition at line 698 of file mpsdammat.C.

References checkSnapBack, E, ecsMethod, oofem::Element::giveCharacteristicSize(), givee0(), oofem::GaussPoint::giveElement(), givegf(), oofem::Element::giveLabel(), oofem::Material::giveStatus(), initDamagedFib(), OOFEM_ERROR, OOFEM_WARNING, softType, ST_Disable_Damage, ST_Exponential_Cohesive_Crack, ST_Linear_Cohesive_Crack, and timeDepFracturing.

Referenced by giveRealStressVector().

initDamagedFib()

void oofem::MPSDamMaterial::initDamagedFib	(	GaussPoint *	gp,
		TimeStep *	tStep ) const

Definition at line 466 of file mpsdammat.C.

References oofem::MPSMaterial::computeEquivalentTime(), computeFractureEnergy(), computeTensileStrength(), E, and oofem::Material::giveStatus().

Referenced by giveRealStressVector(), and initDamaged().

initializeFrom()

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

Definition at line 205 of file mpsdammat.C.

References _IFT_MPSDamMaterial_checkSnapBack, _IFT_MPSDamMaterial_damageLaw, _IFT_MPSDamMaterial_fib_s, _IFT_MPSDamMaterial_ft, _IFT_MPSDamMaterial_ft28, _IFT_MPSDamMaterial_gf, _IFT_MPSDamMaterial_gf28, _IFT_MPSDamMaterial_isotropic, _IFT_MPSDamMaterial_timedepfracturing, _IFT_MPSMaterial_fc, checkSnapBack, const_gf, fib_fcm28, fib_s, ft, ft28, gf28, oofem::InputRecord::hasField(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, isotropic, OOFEM_ERROR, softType, ST_Disable_Damage, ST_Exponential_Cohesive_Crack, ST_Linear_Cohesive_Crack, and timeDepFracturing.

Member Data Documentation

checkSnapBack

int oofem::MPSDamMaterial::checkSnapBack = 1

protected

Check possible snap back flag.

Definition at line 210 of file mpsdammat.h.

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

const_gf

double oofem::MPSDamMaterial::const_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 207 of file mpsdammat.h.

Referenced by givegf(), and initializeFrom().

E

double oofem::MPSDamMaterial::E = -1.

mutableprotected

dummy Young's modulus

Definition at line 190 of file mpsdammat.h.

Referenced by computeDamageForCohesiveCrack(), givee0(), giveRealStressVector(), initDamaged(), and initDamagedFib().

ecsMethod

ElementCharSizeMethod oofem::MPSDamMaterial::ecsMethod = ECSM_Projection

protected

Method used for evaluation of characteristic element size.

Definition at line 220 of file mpsdammat.h.

Referenced by initDamaged().

fib_fcm28

double oofem::MPSDamMaterial::fib_fcm28 = 0.

protected

Definition at line 186 of file mpsdammat.h.

Referenced by computeFractureEnergy(), computeTensileStrength(), and initializeFrom().

fib_s

double oofem::MPSDamMaterial::fib_s = 0.

protected

Definition at line 185 of file mpsdammat.h.

Referenced by computeTensileStrength(), and initializeFrom().

ft

double oofem::MPSDamMaterial::ft = 0.

protected

Equivalent strain at stress peak (or a similar parameter).

constant tensile strength

Definition at line 199 of file mpsdammat.h.

Referenced by givee0(), giveIPValue(), and initializeFrom().

ft28

double oofem::MPSDamMaterial::ft28 = 0.

protected

28-day value of tensile strength. Used only with "timedepfracturing"

Definition at line 223 of file mpsdammat.h.

Referenced by computeTensileStrength(), and initializeFrom().

gf28

double oofem::MPSDamMaterial::gf28 = 0.

protected

28-day value of fracture energy. Used only with "timedepfracturing"

Definition at line 225 of file mpsdammat.h.

Referenced by computeFractureEnergy(), and initializeFrom().

isotropic

bool oofem::MPSDamMaterial::isotropic = false

protected

Definition at line 187 of file mpsdammat.h.

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

maxOmega

double oofem::MPSDamMaterial::maxOmega = 0.999999

protected

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.

Definition at line 193 of file mpsdammat.h.

softType

SofteningType oofem::MPSDamMaterial::softType = ST_Exponential_Cohesive_Crack

protected

Parameter specifying the type of softening (damage law).

Definition at line 217 of file mpsdammat.h.

Referenced by computeDamage(), computeDamageForCohesiveCrack(), giveRealStressVector(), initDamaged(), and initializeFrom().

timeDepFracturing

bool oofem::MPSDamMaterial::timeDepFracturing = false

protected

Definition at line 184 of file mpsdammat.h.

Referenced by givee0(), givegf(), giveRealStressVector(), initDamaged(), and initializeFrom().

---

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

• mpsdammat.h
• mpsdammat.C