#include <mdm.h>

Inheritance diagram for oofem::MDM:

Collaboration diagram for oofem::MDM:

Public Types
enum	MDMFormulatrionType { COMPLIANCE_DAMAGE , STIFFNESS_DAMAGE }
enum	MDMModeType { mdm_3d , mdm_2d }

Public Member Functions
	MDM (int n, Domain *d)
bool	hasMaterialModeCapability (MaterialMode mode) 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.
FloatArrayF< 3 >	giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const override
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
FloatArrayF< 6 >	giveThermalDilatationVector (GaussPoint gp, TimeStep tStep) const override
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
const char *	giveInputRecordName () const override
const char *	giveClassName () const override
void	initializeData (int numberOfMicroplanes) override
void	updateBeforeNonlocAverage (const FloatArray &strainVector, GaussPoint gp, TimeStep tStep) const override
double	computeWeightFunction (const double cl, const FloatArray &src, const FloatArray &coord) const override
int	hasBoundedSupport () const override
virtual void	giveSupportRadius (double &radius)
Interface *	giveInterface (InterfaceType it) override
int	MMI_map (GaussPoint gp, Domain oldd, TimeStep *tStep) override
int	MMI_update (GaussPoint gp, TimeStep tStep, FloatArray *estrain=nullptr) override
int	MMI_finish (TimeStep *tStep) override
int	packUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip) override
int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip) override
int	estimatePackSize (DataStream &buff, GaussPoint *ip) override
double	predictRelativeComputationalCost (GaussPoint *gp) override
double	predictRelativeRedistributionCost (GaussPoint *gp) override
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
Public Member Functions inherited from oofem::MicroplaneMaterial
	MicroplaneMaterial (int n, Domain *d)
double	computeNormalStrainComponent (int mnumber, const FloatArray &macroStrain) const
double	computeNormalVolumetricStrainComponent (const FloatArray &macroStrain) const
double	computeNormalDeviatoricStrainComponent (int mnumber, const FloatArray &macroStrain) const
double	computeShearMStrainComponent (int mnumber, const FloatArray &macroStrain) const
double	computeShearLStrainComponent (int mnumber, const FloatArray &macroStrain) const
MicroplaneState	computeStrainVectorComponents (int mnumber, const FloatArray &macroStrain) const
double	giveMicroplaneIntegrationWeight (int mnumber) const
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const override
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) 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_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< 1 >	giveRealStressVector_1d (const FloatArrayF< 1 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
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 >	give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 1, 1 >	give1dStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	give1dStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 2, 2 >	give2dBeamLayerStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 5, 5 >	givePlateLayerStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 3, 3 >	giveFiberStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 3, 3 >	give2dPlateSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 6, 6 >	give3dBeamSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
Public Member Functions inherited from oofem::Material
	Material (int n, Domain *d)
virtual	~Material ()=default
	Destructor.
virtual bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode) const
virtual double	giveCharacteristicValue (MatResponseMode type, GaussPoint gp, TimeStep tStep) const
	Returns characteristic value of the receiver.
virtual double	give (int aProperty, GaussPoint *gp) const
virtual bool	hasProperty (int aProperty, GaussPoint *gp) const
virtual void	modifyProperty (int aProperty, double value, GaussPoint *gp)
double	giveCastingTime () const
virtual bool	isActivated (TimeStep *tStep) const
virtual bool	hasCastingTimeSupport () const
void	printYourself () override
	Prints receiver state on stdout. Useful for debugging.
virtual void	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
virtual void	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
int	checkConsistency () override
virtual void	restoreConsistency (GaussPoint *gp)
virtual int	initMaterial (Element *element)
virtual MaterialStatus *	giveStatus (GaussPoint *gp) const
virtual void	initTempStatus (GaussPoint *gp) const
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
virtual	~FEMComponent ()=default
	Virtual destructor.
Domain *	giveDomain () const
virtual void	setDomain (Domain *d)
int	giveNumber () const
void	setNumber (int num)
virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
virtual void	initializeFrom (InputRecord &ir, int priority)
virtual void	initializeFinish ()
virtual void	postInitialize ()
	Performs post initialization steps. Called after all components are created and initialized.
virtual void	printOutputAt (FILE file, TimeStep tStep)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).
Public Member Functions inherited from oofem::StructuralNonlocalMaterialExtensionInterface
	StructuralNonlocalMaterialExtensionInterface (Domain *d)
	~StructuralNonlocalMaterialExtensionInterface ()
	Destructor.
Public Member Functions inherited from oofem::NonlocalMaterialExtensionInterface
	NonlocalMaterialExtensionInterface (Domain *d)
virtual	~NonlocalMaterialExtensionInterface ()
	Destructor.
void	updateDomainBeforeNonlocAverage (TimeStep *tStep) const
void	buildNonlocalPointTable (GaussPoint *gp) const
void	rebuildNonlocalPointTable (GaussPoint gp, IntArray contributingElems) const
void	modifyNonlocalWeightFunctionAround (GaussPoint *gp) const
void	modifyNonlocalWeightFunction_1D_Around (GaussPoint *gp) const
virtual double	giveNonlocalMetricModifierAt (GaussPoint *gp) const
double	computeDistanceModifier (double cl, double damage) const
double	computeModifiedLength (double length, double dam1, double dam2) const
std ::vector< localIntegrationRecord > *	giveIPIntegrationList (GaussPoint *gp) const
virtual double	computeWeightFunction (const double cl, const double distance) const
double	giveIntegralOfWeightFunction (double cl, const int spatial_dimension) const
virtual double	maxValueOfWeightFunction ()
	Determines the maximum value of the nonlocal weight function.
int	giveNumberOfRegions () const
virtual double	evaluateSupportRadius (double cl) const
Domain *	giveDomain ()
	Returns reference to domain.
void	initializeFrom (InputRecord &ir)
void	giveInputRecord (DynamicInputRecord &input)
void	endIPNonlocalAverage (GaussPoint *gp) const
Public Member Functions inherited from oofem::Interface
	Interface ()
	Constructor.
virtual	~Interface ()
Public Member Functions inherited from oofem::MaterialModelMapperInterface
	MaterialModelMapperInterface ()
	Constructor.
virtual	~MaterialModelMapperInterface ()
	Destructor.

Protected Types
enum	MDMMapperType { mdm_cpt =0 , mdm_sft =1 , mdm_lst =2 }
Protected Types inherited from oofem::NonlocalMaterialExtensionInterface
enum	ScalingType { ST_Unknown , ST_Standard , ST_Noscaling , ST_Borino }
	Type characterizing the scaling approach. More...
enum	AveragedVarType { AVT_Unknown , AVT_EqStrain , AVT_Compliance , AVT_Damage }
	Type characterizing the averaged (nonlocal) variable. More...
enum	NlVariationType { NLVT_Standard , NLVT_DistanceBasedLinear , NLVT_StressBased , NLVT_DistanceBasedExponential }
	Type characterizing the Nonlocal variation. More...

Protected Member Functions
void	computeDamageTensor (FloatMatrix &tempDamageTensor, const FloatArray &totalStrain, GaussPoint gp, TimeStep tStep) const
void	computeLocalDamageTensor (FloatMatrix &tempDamageTensor, const FloatArray &totalStrain, GaussPoint gp, TimeStep tStep) const
double	computeDamageOnPlane (GaussPoint *gp, int mnumber, const FloatArray &strain) const
void	computePDC (FloatMatrix &tempDamageTensor, FloatArray &tempDamageTensorEigenVals, FloatMatrix &tempDamageTensorEigenVec) const
void	transformStrainToPDC (FloatArray &answer, FloatArray &strain, FloatMatrix &t, GaussPoint *gp) const
void	applyDamageTranformation (FloatArray &strainPDC, const FloatArray &tempDamageTensorEigenVals) const
void	transformStressFromPDC (FloatArray &answer, const FloatArray &stressPDC, const FloatMatrix &t, GaussPoint *gp) const
void	computeEffectiveStress (FloatArray &stressPDC, const FloatArray &strainPDC, GaussPoint gp, TimeStep tStep) const
void	giveMaterialStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep) const
void	applyDamageToStiffness (FloatMatrix &d, GaussPoint *gp) const
void	transformStiffnessfromPDC (FloatMatrix &de, const FloatMatrix &t) const
FloatMatrixF< 3, 3 >	givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 4, 4 >	givePlaneStrainStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
void	rotateTensor4 (FloatMatrix &Dlocal, const FloatMatrix &t) const
void	formTransformationMatrix (FloatMatrix &answer, const FloatMatrix &t, int n) const
void	giveRawMDMParameters (double &Efp, double &Ep, const FloatArray &reducedStrain, GaussPoint *gp) const
Protected Member Functions inherited from oofem::NonlocalMaterialExtensionInterface
void	applyBarrierConstraints (const FloatArray &gpCoords, const FloatArray &jGpCoords, double &weight) const
void	manipulateWeight (double &weight, GaussPoint gp, GaussPoint jGp) const
double	giveDistanceBasedInteractionRadius (const FloatArray &gpCoords) const
int	mapToGridPoint (double x, double x0) const
double	mapToGridCoord (double x, double x0) const
double	dist2FromGridNode (double x, double y, int j, int i) const

Protected Attributes
int	ndc = 0
	Number of damage components.
int	nsd = 0
	Number of spatial dimensions.
int	type_dam = 0
int	type_soft = 0
double	mdm_Ep = 0.
	Parameter controlling the elastic limit.
double	mdm_Efp = 0.
	Prescribed value of ef-ep.
double	ParMd = 0.
	THREAD UNSAFE. These are modified when evaluating the material. Why? This must be avoided at all cost.
double	tempDillatCoeff = 0.
	Temperature dilatation coeff.
double	Gf = 0.
	Fracture energy (necessary to determine Ep and Efp if not given).
double	Ft = 0.
	Macroscopic tensile strength (necessary to determine Ep and Efp if not given).
MDMFormulatrionType	formulation = COMPLIANCE_DAMAGE
MDMModeType	mdmMode = mdm_3d
IsotropicLinearElasticMaterial	linearElasticMaterial
	Reference to bulk (undamaged) material.
int	nonlocal = 0
	Flag indicating local or nonlocal mode.
double	R = 0.
	Interaction radius, related to the nonlocal characteristic length of material.
std::unique_ptr< Set >	sourceElemSet
	cached source element set used to map internal variables (adaptivity), created on demand
MDMMapperType	mapperType = mdm_cpt
Protected Attributes inherited from oofem::MicroplaneMaterial
int	numberOfMicroplanes = 0
	Number of microplanes.
FloatArray	microplaneWeights
	Integration weights of microplanes.
std::vector< FloatArrayF< 3 > >	microplaneNormals
	Normals of microplanes.
FloatArrayF< 6 >	Kronecker
	Kronecker's delta.
std::vector< FloatArrayF< 6 > >	N
std::vector< FloatArrayF< 6 > >	M
std::vector< FloatArrayF< 6 > >	L
double	E = 0.
	Young's modulus.
double	nu = 0.
	Poisson's ratio.
Protected Attributes inherited from oofem::StructuralMaterial
double	referenceTemperature = 0.
	Reference temperature (temperature, when material has been built into structure).
MatResponseMode	SCStiffMode = TangentStiffness
	stifness mode used in stress control
double	SCRelTol = 1.e-3
	relative tolerance for stress control
double	SCAbsTol = 1.e-12
	absolute stress tolerance for stress control
int	SCMaxiter = 100000
	maximum iterations for stress-control
Protected Attributes inherited from oofem::Material
Dictionary	propertyDictionary
double	castingTime
int	preCastingTimeMat
	Material existing before casting time - optional parameter, zero by default.
Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number.
Domain *	domain
	Link to domain object, useful for communicating with other FEM components.
Protected Attributes inherited from oofem::NonlocalMaterialExtensionInterface
Domain *	domain = nullptr
IntArray	regionMap
	Map indicating regions to skip (region - cross section model).
bool	permanentNonlocTableFlag = false
	Flag indicating whether to keep nonlocal interaction tables of integration points cached.
WeightFunctionType	weightFun
	Parameter specifying the type of nonlocal weight function.
int	gridSize = 0
	Grid on which the eikonal equation will be solved (used by eikonal nonlocal models).
double	initDiag = 0.
	Optional parameters setting details of the fast marching method.
int	order = 0
int	centDiff = 0
double	cl = 0.
double	suprad = 0.
	Support radius.
double	mm = 0.
	For "undernonlocal" or "overnonlocal" formulation.
ScalingType	scaling
	Parameter specifying the type of scaling of nonlocal weight function.
AveragedVarType	averagedVar
	Parameter specifying the type of averaged (nonlocal) variable.
double	cl0 = 0.
NlVariationType	nlvar
	Parameter specifying the type of nonlocal variation.
double	beta = 0.
double	zeta = 0.
double	px = 0.
double	Rf = 0.
	Final value of interaction radius, for a model with evolving characteristic length.
double	exponent = 0.
	Parameter used as an exponent by models with evolving characteristic length.
int	averType = 0
	Parameter specifying how the weight function should be adjusted due to damage.

Static Protected Attributes
static MMAShapeFunctProjection	mapperSFT
	Mapper used to map internal variables in adaptivity.
static MMALeastSquareProjection	mapperLST
static MMAClosestIPTransfer	mapper2
	Mapper used to map stresses in adaptivity.

Additional Inherited Members
Static Public Member Functions inherited from oofem::StructuralMaterial
static int	giveSymVI (int ind1, int ind2)
static int	giveVI (int ind1, int ind2)
static FloatMatrixF< 9, 9 >	convert_dSdE_2_dPdF_3D (const FloatMatrixF< 6, 6 > &dSdE, const FloatArrayF< 6 > &S, const FloatArrayF< 9 > &F)
static FloatMatrixF< 5, 5 >	convert_dSdE_2_dPdF_PlaneStrain (const FloatMatrixF< 4, 4 > &dSdE, const FloatArrayF< 4 > &S, const FloatArrayF< 5 > &F)
static FloatMatrixF< 4, 4 >	convert_dSdE_2_dPdF_PlaneStress (const FloatMatrixF< 3, 3 > &dSdE, const FloatArrayF< 3 > &S, const FloatArrayF< 4 > &F)
static FloatMatrixF< 1, 1 >	convert_dSdE_2_dPdF_1D (const FloatMatrixF< 1, 1 > &dSdE, const FloatArrayF< 1 > &S, const FloatArrayF< 1 > &F)
static void	computePrincipalValues (FloatArray &answer, const FloatArray &s, stressStrainPrincMode mode)
	Common functions for convenience.
static FloatArrayF< 3 >	computePrincipalValues (const FloatMatrixF< 3, 3 > &s)
static FloatArrayF< 3 >	computePrincipalValues (double I1, double I2, double I3)
static void	computePrincipalValDir (FloatArray &answer, FloatMatrix &dir, const FloatArray &s, stressStrainPrincMode mode)
static std::pair< FloatArrayF< 3 >, FloatMatrixF< 3, 3 > >	computePrincipalValDir (const FloatMatrixF< 3, 3 > &s)
static FloatArrayF< 6 >	computeDeviator (const FloatArrayF< 6 > &s)
static std::pair< FloatArrayF< 6 >, double >	computeDeviatoricVolumetricSplit (const FloatArrayF< 6 > &s)
static FloatArrayF< 6 >	computeDeviatoricVolumetricSum (const FloatArrayF< 6 > &dev, double mean)
static FloatArrayF< 6 >	applyDeviatoricElasticCompliance (const FloatArrayF< 6 > &stress, double EModulus, double nu)
static FloatArrayF< 6 >	applyDeviatoricElasticCompliance (const FloatArrayF< 6 > &stress, double GModulus)
static FloatArrayF< 6 >	applyDeviatoricElasticStiffness (const FloatArrayF< 6 > &strain, double EModulus, double nu)
static FloatArrayF< 6 >	applyDeviatoricElasticStiffness (const FloatArrayF< 6 > &strain, double GModulus)
static FloatArrayF< 6 >	applyElasticStiffness (const FloatArrayF< 6 > &strain, double EModulus, double nu)
static FloatArrayF< 6 >	applyElasticCompliance (const FloatArrayF< 6 > &stress, double EModulus, double nu)
static double	computeStressNorm (const FloatArrayF< 6 > &stress)
static double	computeFirstInvariant (const FloatArrayF< 6 > &s)
static double	computeSecondStressInvariant (const FloatArrayF< 6 > &s)
static double	computeThirdStressInvariant (const FloatArrayF< 6 > &s)
static double	computeFirstCoordinate (const FloatArrayF< 6 > &s)
static double	computeSecondCoordinate (const FloatArrayF< 6 > &s)
static double	computeThirdCoordinate (const FloatArrayF< 6 > &s)
static int	giveVoigtVectorMask (IntArray &answer, MaterialMode mmode)
static int	giveVoigtSymVectorMask (IntArray &answer, MaterialMode mmode)
static void	giveInvertedVoigtVectorMask (IntArray &answer, MaterialMode mmode)
static int	giveSizeOfVoigtVector (MaterialMode mmode)
static int	giveSizeOfVoigtSymVector (MaterialMode mmode)
static void	giveFullVectorForm (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
	Converts the reduced symmetric Voigt vector (2nd order tensor) to full form.
static void	giveFullVectorFormF (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
	Converts the reduced deformation gradient Voigt vector (2nd order tensor).
static void	giveFullSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the reduced unsymmetric Voigt vector (2nd order tensor) to full form.
static void	giveReducedVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the full symmetric Voigt vector (2nd order tensor) to reduced form.
static void	giveReducedSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the full unsymmetric Voigt vector (2nd order tensor) to reduced form.
static void	giveFullSymMatrixForm (FloatMatrix &answer, const FloatMatrix &red, MaterialMode matMode)
	Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form.
static void	giveReducedMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
	Converts the full symmetric Voigt matrix (4th order tensor) to reduced form.
static void	giveReducedSymMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
	Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form.
static FloatArrayF< 6 >	transformStrainVectorTo (const FloatMatrixF< 3, 3 > &base, const FloatArrayF< 6 > &strain, bool transpose=false)
static FloatArrayF< 6 >	transformStressVectorTo (const FloatMatrixF< 3, 3 > &base, const FloatArrayF< 6 > &stress, bool transpose=false)
static double	computeVonMisesStress (const FloatArray &currentStress)
static double	computeVonMisesStress_3D (const FloatArrayF< 6 > &stress)
static double	computeVonMisesStress_PlaneStress (const FloatArrayF< 3 > &stress)
static FloatMatrixF< 6, 6 >	giveStrainVectorTranformationMtrx (const FloatMatrixF< 3, 3 > &base, bool transpose=false)
static FloatMatrixF< 3, 3 >	give2DStrainVectorTranformationMtrx (const FloatMatrixF< 2, 2 > &base, bool transpose=false)
static FloatMatrixF< 6, 6 >	giveStressVectorTranformationMtrx (const FloatMatrixF< 3, 3 > &base, bool transpose=false)
static FloatMatrixF< 3, 3 >	givePlaneStressVectorTranformationMtrx (const FloatMatrixF< 2, 2 > &base, bool transpose=false)
static void	sortPrincDirAndValCloseTo (FloatArray &pVal, FloatMatrix &pDir, const FloatMatrix &toPDir)
Static Public Attributes inherited from oofem::StructuralMaterial
static std::array< std::array< int, 3 >, 3 >	vIindex
	Voigt index map.
static std::array< std::array< int, 3 >, 3 >	svIndex
	Symmetric Voigt index map.

Detailed Description

Implementation of microplane damage material (According to M.Jirasek).

Definition at line 158 of file mdm.h.

Member Enumeration Documentation

◆ MDMFormulatrionType

enum oofem::MDM::MDMFormulatrionType

Enumerator
COMPLIANCE_DAMAGE
STIFFNESS_DAMAGE

Definition at line 161 of file mdm.h.

◆ MDMMapperType

enum oofem::MDM::MDMMapperType

protected

Enumerator
mdm_cpt
mdm_sft
mdm_lst

Definition at line 202 of file mdm.h.

◆ MDMModeType

enum oofem::MDM::MDMModeType

Enumerator
mdm_3d
mdm_2d

Definition at line 162 of file mdm.h.

Constructor & Destructor Documentation

◆ MDM()

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

inline

Constructor. Creates Microplane Material belonging to domain d, with number n.

Parameters

n	Material number.
d	Domain to which newly created material belongs.

Definition at line 229 of file mdm.h.

References linearElasticMaterial, oofem::MaterialModelMapperInterface::MaterialModelMapperInterface(), mdm_Efp, mdm_Ep, oofem::MicroplaneMaterial::MicroplaneMaterial(), and oofem::StructuralNonlocalMaterialExtensionInterface::StructuralNonlocalMaterialExtensionInterface().

Referenced by give3dMaterialStiffnessMatrix(), giveMaterialStiffnessMatrix(), givePlaneStrainStiffMtrx(), givePlaneStressStiffMtrx(), giveRealStressVector_3d(), giveRealStressVector_PlaneStrain(), and giveRealStressVector_PlaneStress().

Member Function Documentation

◆ applyDamageToStiffness()

void oofem::MDM::applyDamageToStiffness	(	FloatMatrix &	d,
		GaussPoint *	gp ) const

protected

Definition at line 538 of file mdm.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), COMPLIANCE_DAMAGE, formulation, oofem::FloatMatrix::giveNumberOfRows(), oofem::Material::giveStatus(), oofem::MDMStatus::giveTempDamageTensorEigenVals(), MAX_REL_COMPL_TRESHOLD, mdm_2d, mdm_3d, mdmMode, oofem::min(), OOFEM_ERROR, and STIFFNESS_DAMAGE.

Referenced by giveMaterialStiffnessMatrix().

◆ applyDamageTranformation()

void oofem::MDM::applyDamageTranformation	(	FloatArray &	strainPDC,
		const FloatArray &	tempDamageTensorEigenVals ) const

protected

Definition at line 403 of file mdm.C.

References oofem::FloatArray::at(), COMPLIANCE_DAMAGE, formulation, mdm_2d, mdm_3d, mdmMode, OOFEM_ERROR, and STIFFNESS_DAMAGE.

Referenced by giveRealStressVector().

◆ computeDamageOnPlane()

double oofem::MDM::computeDamageOnPlane	(	GaussPoint *	gp,
		int	mnumber,
		const FloatArray &	strain ) const

protected

◆ computeDamageTensor()

void oofem::MDM::computeDamageTensor	(	FloatMatrix &	tempDamageTensor,
		const FloatArray &	totalStrain,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

protected

Definition at line 153 of file mdm.C.

References oofem::FloatMatrix::at(), oofem::NonlocalMaterialExtensionInterface::buildNonlocalPointTable(), computeLocalDamageTensor(), oofem::NonlocalMaterialExtensionInterface::endIPNonlocalAverage(), oofem::NonlocalMaterialStatusExtensionInterface::giveIntegrationScale(), oofem::NonlocalMaterialExtensionInterface::giveIPIntegrationList(), oofem::MDMStatus::giveLocalDamageTensorForAverage(), oofem::Material::giveStatus(), ndc, nonlocal, nsd, oofem::FloatMatrix::symmetrized(), oofem::FloatMatrix::times(), and oofem::NonlocalMaterialExtensionInterface::updateDomainBeforeNonlocAverage().

Referenced by giveRealStressVector().

◆ computeEffectiveStress()

void oofem::MDM::computeEffectiveStress	(	FloatArray &	stressPDC,
		const FloatArray &	strainPDC,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

protected

Definition at line 452 of file mdm.C.

References oofem::FloatArray::beProductOf(), linearElasticMaterial, mdm_3d, and mdmMode.

Referenced by giveRealStressVector().

◆ computeLocalDamageTensor()

void oofem::MDM::computeLocalDamageTensor	(	FloatMatrix &	tempDamageTensor,
		const FloatArray &	totalStrain,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

protected

◆ computePDC()

void oofem::MDM::computePDC	(	FloatMatrix &	tempDamageTensor,
		FloatArray &	tempDamageTensorEigenVals,
		FloatMatrix &	tempDamageTensorEigenVec ) const

protected

Definition at line 336 of file mdm.C.

References oofem::FloatMatrix::jaco_(), nsd, oofem::FloatArray::resize(), and oofem::FloatMatrix::resize().

Referenced by giveRealStressVector().

◆ computeWeightFunction()

double oofem::MDM::computeWeightFunction	(	const double	cl,
		const FloatArray &	src,
		const FloatArray &	coord ) const

overridevirtual

Evaluates the basic nonlocal weight function for two points with given coordinates. This function is NOT normalized by the condition of unit integral.

Parameters

src	Coordinates of source point.
coord	Coordinates of receiver point.

Returns: Value of weight function.

Reimplemented from oofem::NonlocalMaterialExtensionInterface.

Definition at line 1154 of file mdm.C.

References oofem::distance(), and R.

◆ CreateStatus()

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

Definition at line 74 of file mdm.C.

References nsd, and oofem::MicroplaneMaterial::numberOfMicroplanes.

◆ estimatePackSize()

int oofem::MDM::estimatePackSize	(	DataStream &	buff,
		GaussPoint *	ip )

overridevirtual

Estimates the necessary pack size to hold all packed data of receiver.

Reimplemented from oofem::Material.

Definition at line 1323 of file mdm.C.

References oofem::FloatMatrix::givePackSize(), nonlocal, and nsd.

◆ formTransformationMatrix()

void oofem::MDM::formTransformationMatrix	(	FloatMatrix &	answer,
		const FloatMatrix &	t,
		int	n ) const

protected

Definition at line 921 of file mdm.C.

References FORMT33, FORMT44, FORMT66, and OOFEM_ERROR.

Referenced by rotateTensor4().

◆ give3dMaterialStiffnessMatrix()

FloatMatrixF< 6, 6 > oofem::MDM::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 630 of file mdm.C.

References giveMaterialStiffnessMatrix(), and MDM().

◆ giveClassName()

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

inlineoverridevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 271 of file mdm.h.

◆ giveInputRecord()

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

overridevirtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::Material.

Definition at line 855 of file mdm.C.

References _IFT_MDM_efp, _IFT_MDM_ep, _IFT_MDM_formulation, _IFT_MDM_ft, _IFT_MDM_gf, _IFT_MDM_mapper, _IFT_MDM_mode, _IFT_MDM_nonloc, _IFT_MDM_parmd, _IFT_MDM_r, _IFT_MDM_talpha, formulation, Ft, Gf, linearElasticMaterial, mapper2, mapperLST, mapperSFT, mapperType, mdm_Efp, mdm_Ep, mdmMode, nonlocal, ParMd, R, oofem::DynamicInputRecord::setField(), and tempDillatCoeff.

◆ giveInputRecordName()

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

inlineoverridevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 270 of file mdm.h.

References _IFT_MDM_Name.

◆ giveInterface()

Interface * oofem::MDM::giveInterface ( InterfaceType t )

overridevirtual

Interface requesting service.

See also: InterfaceType

Returns: Requested interface if implemented, otherwise NULL.

Reimplemented from oofem::FEMComponent.

Definition at line 1288 of file mdm.C.

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

◆ giveIPValue()

int oofem::MDM::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 659 of file mdm.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beInverseOf(), oofem::FloatMatrix::beProductOf(), COMPLIANCE_DAMAGE, formulation, oofem::MDMStatus::giveDamageTensor(), oofem::MDMStatus::giveDamageTensorEigenVals(), oofem::MDMStatus::giveMicroplaneDamageValues(), oofem::Material::giveStatus(), oofem::MDMStatus::giveTempDamageTensor(), oofem::MDMStatus::giveTempDamageTensorEigenVals(), ndc, nsd, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

◆ giveMaterialStiffnessMatrix()

void oofem::MDM::giveMaterialStiffnessMatrix	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

protected

FIXME: ask for 3d stiffness?

Definition at line 513 of file mdm.C.

References applyDamageToStiffness(), oofem::Material::giveStatus(), linearElasticMaterial, MDM(), and transformStiffnessfromPDC().

Referenced by give3dMaterialStiffnessMatrix(), givePlaneStrainStiffMtrx(), and givePlaneStressStiffMtrx().

◆ givePlaneStrainStiffMtrx()

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

overrideprotectedvirtual

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 650 of file mdm.C.

References giveMaterialStiffnessMatrix(), and MDM().

◆ givePlaneStressStiffMtrx()

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

overrideprotectedvirtual

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

Parameters

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

Reimplemented from oofem::StructuralMaterial.

Definition at line 641 of file mdm.C.

References giveMaterialStiffnessMatrix(), and MDM().

◆ giveRawMDMParameters()

void oofem::MDM::giveRawMDMParameters	(	double &	Efp,
		double &	Ep,
		const FloatArray &	reducedStrain,
		GaussPoint *	gp ) const

protected

Definition at line 972 of file mdm.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::StructuralMaterial::computePrincipalValDir(), oofem::MicroplaneMaterial::E, Ex, Ft, Gf, oofem::Element::giveCharacteristicLength(), oofem::GaussPoint::giveElement(), oofem::GaussPoint::giveMaterialMode(), linearElasticMaterial, oofem::max(), mdm_Efp, mdm_Ep, nonlocal, nsd, OOFEM_ERROR, ParMd, oofem::principal_strain, and R.

Referenced by computeDamageOnPlane().

◆ giveRealStressVector()

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

Definition at line 87 of file mdm.C.

References applyDamageTranformation(), oofem::FloatArray::at(), computeDamageTensor(), computeEffectiveStress(), computePDC(), oofem::GaussPoint::giveElement(), oofem::FEMComponent::giveNumber(), oofem::GaussPoint::giveNumber(), oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStressDependentPartOfStrainVector(), oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), nsd, OOFEM_ERROR, oofem::MDMStatus::setTempDamageTensor(), oofem::MDMStatus::setTempDamageTensorEigenVals(), oofem::MDMStatus::setTempDamageTensorEigenVec(), transformStrainToPDC(), and transformStressFromPDC().

Referenced by giveRealStressVector_3d(), giveRealStressVector_PlaneStrain(), giveRealStressVector_PlaneStress(), and MMI_update().

◆ giveRealStressVector_3d()

FloatArrayF< 6 > oofem::MDM::giveRealStressVector_3d	(	const FloatArrayF< 6 > &	strain,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

inlineoverridevirtual

Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.

Reimplemented from oofem::StructuralMaterial.

Definition at line 242 of file mdm.h.

References giveRealStressVector(), and MDM().

◆ giveRealStressVector_PlaneStrain()

FloatArrayF< 4 > oofem::MDM::giveRealStressVector_PlaneStrain	(	const FloatArrayF< 4 > &	strain,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

inlineoverridevirtual

Default implementation relies on giveRealStressVector_3d.

Reimplemented from oofem::StructuralMaterial.

Definition at line 248 of file mdm.h.

References giveRealStressVector(), and MDM().

◆ giveRealStressVector_PlaneStress()

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

inlineoverridevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 254 of file mdm.h.

References giveRealStressVector(), and MDM().

◆ giveSupportRadius()

virtual void oofem::MDM::giveSupportRadius ( double & radius )

inlinevirtual

Determines the width (radius) of limited support of weighting function.

Definition at line 283 of file mdm.h.

◆ giveThermalDilatationVector()

FloatArrayF< 6 > oofem::MDM::giveThermalDilatationVector	(	GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Returns a vector of coefficients of thermal dilatation in direction of each material principal (local) axis.

Parameters

answer	Vector of thermal dilatation coefficients.
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 763 of file mdm.C.

References tempDillatCoeff.

◆ hasBoundedSupport()

int oofem::MDM::hasBoundedSupport ( ) const

inlineoverridevirtual

Determines, whether receiver has bounded weighting function (limited support).

Returns: True if weighting function bounded, zero otherwise.

Reimplemented from oofem::NonlocalMaterialExtensionInterface.

Definition at line 279 of file mdm.h.

◆ hasMaterialModeCapability()

bool oofem::MDM::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 81 of file mdm.C.

◆ initializeData()

void oofem::MDM::initializeData ( int numberOfMicroplanes )

overridevirtual

Initializes internal data (integration weights, microplane normals and computes projection tensors).

Parameters

numberOfMicroplanes Number of required microplanes.

Reimplemented from oofem::MicroplaneMaterial.

Definition at line 1083 of file mdm.C.

References oofem::FloatArray::at(), oofem::MicroplaneMaterial::L, oofem::MicroplaneMaterial::M, M_PI, MAX_NUMBER_OF_MICROPLANES, mdm_2d, mdm_3d, mdmMode, oofem::MicroplaneMaterial::microplaneNormals, oofem::MicroplaneMaterial::microplaneWeights, oofem::MicroplaneMaterial::N, oofem::MicroplaneMaterial::numberOfMicroplanes, and OOFEM_ERROR.

◆ initializeFrom()

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

Definition at line 775 of file mdm.C.

References _IFT_MDM_efp, _IFT_MDM_ep, _IFT_MDM_formulation, _IFT_MDM_ft, _IFT_MDM_gf, _IFT_MDM_mapper, _IFT_MDM_mode, _IFT_MDM_nonloc, _IFT_MDM_parmd, _IFT_MDM_r, _IFT_MDM_talpha, formulation, Ft, Gf, oofem::InputRecord::hasField(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, linearElasticMaterial, mapper2, mapperLST, mapperSFT, mapperType, mdm_2d, mdm_3d, mdm_Efp, mdm_Ep, mdmMode, ndc, nonlocal, nsd, OOFEM_LOG_INFO, ParMd, R, and tempDillatCoeff.

◆ MMI_finish()

int oofem::MDM::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 1274 of file mdm.C.

References mapper2, mapperLST, and mapperSFT.

◆ MMI_map()

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

overridevirtual

Maps the required internal state variables from old mesh oldd to given ip. The result is stored in gp status.

Parameters

gp	Integration point belonging to new domain which values will be mapped
oldd	Old mesh reference.
tStep	Time step.

Returns: Nonzero if o.k.

Implements oofem::MaterialModelMapperInterface.

Definition at line 1170 of file mdm.C.

References oofem::FloatArray::at(), oofem::IntArray::at(), COMPLIANCE_DAMAGE, oofem::IntArray::followedBy(), formulation, oofem::Domain::giveElement(), oofem::Element::giveMaterial(), oofem::FEMComponent::giveNumber(), oofem::Domain::giveNumberOfElements(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), mapper2, mapperLST, mapperSFT, mapperType, mdm_cpt, mdm_lst, mdm_sft, OOFEM_ERROR, oofem::IntArray::resize(), oofem::MDMStatus::setMicroplaneTempDamageValues(), sourceElemSet, and oofem::MDMStatus::updateYourself().

◆ MMI_update()

int oofem::MDM::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

gp	Integration point belonging to new domain which values will be mapped.
tStep	Time step.
elemGPVec	Vector 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 1259 of file mdm.C.

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

◆ packUnknowns()

int oofem::MDM::packUnknowns	(	DataStream &	buff,
		TimeStep *	tStep,
		GaussPoint *	ip )

overridevirtual

Pack all necessary data of integration point (according to element parallel_mode) into given communication buffer. The nature of packed data is material model dependent. Typically, for material of "local" response (response depends only on integration point local state) no data are exchanged. For "nonlocal" constitutive models the send/receive of local values which undergo averaging is performed between local and corresponding remote elements.

Parameters

buff	Communication buffer.
tStep	Solution step.
ip	Integration point.

Reimplemented from oofem::Material.

Definition at line 1300 of file mdm.C.

References oofem::NonlocalMaterialExtensionInterface::buildNonlocalPointTable(), oofem::CIO_OK, oofem::MDMStatus::giveLocalDamageTensorForAveragePtr(), oofem::Material::giveStatus(), oofem::FloatMatrix::storeYourself(), and oofem::NonlocalMaterialExtensionInterface::updateDomainBeforeNonlocAverage().

◆ predictRelativeComputationalCost()

double oofem::MDM::predictRelativeComputationalCost ( GaussPoint * gp )

overridevirtual

Returns the weight representing relative computational cost of receiver The reference material model is linear isotropic material - its weight is set to 1.0 The other material models should compare to this reference model.

Reimplemented from oofem::Material.

Definition at line 1338 of file mdm.C.

References oofem::NonlocalMaterialStatusExtensionInterface::giveIntegrationDomainList(), oofem::Material::giveStatus(), nonlocal, and nsd.

◆ predictRelativeRedistributionCost()

double oofem::MDM::predictRelativeRedistributionCost ( GaussPoint * gp )

inlineoverridevirtual

Returns the relative redistribution cost of the receiver

Reimplemented from oofem::Material.

Definition at line 295 of file mdm.h.

◆ rotateTensor4()

void oofem::MDM::rotateTensor4	(	FloatMatrix &	Dlocal,
		const FloatMatrix &	t ) const

protected

Definition at line 901 of file mdm.C.

References oofem::FloatMatrix::beProductOf(), formTransformationMatrix(), and oofem::FloatMatrix::giveNumberOfRows().

Referenced by transformStiffnessfromPDC().

◆ transformStiffnessfromPDC()

void oofem::MDM::transformStiffnessfromPDC	(	FloatMatrix &	de,
		const FloatMatrix &	t ) const

protected

Definition at line 623 of file mdm.C.

References rotateTensor4().

Referenced by giveMaterialStiffnessMatrix().

◆ transformStrainToPDC()

void oofem::MDM::transformStrainToPDC	(	FloatArray &	answer,
		FloatArray &	strain,
		FloatMatrix &	t,
		GaussPoint *	gp ) const

protected

Definition at line 358 of file mdm.C.

References oofem::FloatArray::at(), E, oofem::GaussPoint::giveMaterialMode(), mdm_2d, mdm_3d, mdmMode, N, and oofem::FloatArray::resize().

Referenced by giveRealStressVector().

◆ transformStressFromPDC()

void oofem::MDM::transformStressFromPDC	(	FloatArray &	answer,
		const FloatArray &	stressPDC,
		const FloatMatrix &	t,
		GaussPoint *	gp ) const

protected

Definition at line 472 of file mdm.C.

References oofem::FloatArray::at(), oofem::GaussPoint::giveMaterialMode(), mdm_2d, mdm_3d, mdmMode, Nt, oofem::FloatArray::resize(), and S.

Referenced by giveRealStressVector().

◆ unpackAndUpdateUnknowns()

int oofem::MDM::unpackAndUpdateUnknowns	(	DataStream &	buff,
		TimeStep *	tStep,
		GaussPoint *	ip )

overridevirtual

Unpack and updates all necessary data of given integration point (according to element parallel_mode) into given communication buffer.

See also: packUnknowns service.

Parameters

buff	Communication buffer.
tStep	Solution step.
ip	Integration point.

Reimplemented from oofem::Material.

Definition at line 1311 of file mdm.C.

References oofem::CIO_OK, oofem::Material::giveStatus(), oofem::FloatMatrix::restoreYourself(), and oofem::MDMStatus::setLocalDamageTensorForAverage().

◆ updateBeforeNonlocAverage()

void oofem::MDM::updateBeforeNonlocAverage	(	const FloatArray &	strainVector,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Declares the service updating local variables in given integration points, which take part in nonlocal average process. Because value of single integration point influences nonlocal variables in several near integration points, it is suitable to compute these variables only once. These should be stored in integration point associated statuses. The implementation is left on derived classes. Provide material local strain increment - as is provided to computeRealStresVector. This allows to update internal vars to be averaged to new state

Parameters

strainVector	Total strain vector in given integration point.
gp	Integration point to update.
tStep	Solution step indicating time of update.

Implements oofem::StructuralNonlocalMaterialExtensionInterface.

Definition at line 1138 of file mdm.C.

References computeLocalDamageTensor(), and oofem::Material::giveStatus().

Member Data Documentation

◆ formulation

MDMFormulatrionType oofem::MDM::formulation = COMPLIANCE_DAMAGE

protected

Definition at line 182 of file mdm.h.

Referenced by applyDamageToStiffness(), applyDamageTranformation(), computeLocalDamageTensor(), giveInputRecord(), giveIPValue(), initializeFrom(), and MMI_map().

◆ Ft

double oofem::MDM::Ft = 0.

protected

Macroscopic tensile strength (necessary to determine Ep and Efp if not given).

Definition at line 180 of file mdm.h.

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

◆ Gf

double oofem::MDM::Gf = 0.

protected

Fracture energy (necessary to determine Ep and Efp if not given).

Definition at line 178 of file mdm.h.

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

◆ linearElasticMaterial

IsotropicLinearElasticMaterial oofem::MDM::linearElasticMaterial

protected

Reference to bulk (undamaged) material.

Definition at line 186 of file mdm.h.

Referenced by computeDamageOnPlane(), computeEffectiveStress(), giveInputRecord(), giveMaterialStiffnessMatrix(), giveRawMDMParameters(), initializeFrom(), and MDM().

◆ mapper2

MMAClosestIPTransfer oofem::MDM::mapper2

staticprotected

Mapper used to map stresses in adaptivity.

Definition at line 221 of file mdm.h.

Referenced by giveInputRecord(), initializeFrom(), MMI_finish(), and MMI_map().

◆ mapperLST

MMALeastSquareProjection oofem::MDM::mapperLST

staticprotected

Definition at line 199 of file mdm.h.

Referenced by giveInputRecord(), initializeFrom(), MMI_finish(), and MMI_map().

◆ mapperSFT

MMAShapeFunctProjection oofem::MDM::mapperSFT

staticprotected

Mapper used to map internal variables in adaptivity.

Definition at line 198 of file mdm.h.

Referenced by giveInputRecord(), initializeFrom(), MMI_finish(), and MMI_map().

◆ mapperType

MDMMapperType oofem::MDM::mapperType = mdm_cpt

protected

Definition at line 203 of file mdm.h.

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

◆ mdm_Efp

double oofem::MDM::mdm_Efp = 0.

mutableprotected

Prescribed value of ef-ep.

Definition at line 173 of file mdm.h.

Referenced by giveInputRecord(), giveRawMDMParameters(), initializeFrom(), and MDM().

◆ mdm_Ep

double oofem::MDM::mdm_Ep = 0.

mutableprotected

Parameter controlling the elastic limit.

Definition at line 171 of file mdm.h.

Referenced by giveInputRecord(), giveRawMDMParameters(), initializeFrom(), and MDM().

◆ mdmMode

MDMModeType oofem::MDM::mdmMode = mdm_3d

protected

Definition at line 183 of file mdm.h.

Referenced by applyDamageToStiffness(), applyDamageTranformation(), computeEffectiveStress(), giveInputRecord(), initializeData(), initializeFrom(), transformStrainToPDC(), and transformStressFromPDC().

◆ ndc

int oofem::MDM::ndc = 0

protected

Number of damage components.

Definition at line 165 of file mdm.h.

Referenced by computeDamageTensor(), computeLocalDamageTensor(), giveIPValue(), and initializeFrom().

◆ nonlocal

int oofem::MDM::nonlocal = 0

protected

Flag indicating local or nonlocal mode.

Definition at line 189 of file mdm.h.

Referenced by computeDamageTensor(), estimatePackSize(), giveInputRecord(), giveRawMDMParameters(), initializeFrom(), and predictRelativeComputationalCost().

◆ nsd

int oofem::MDM::nsd = 0

protected

Number of spatial dimensions.

Definition at line 166 of file mdm.h.

Referenced by computeDamageOnPlane(), computeDamageTensor(), computePDC(), CreateStatus(), estimatePackSize(), giveIPValue(), giveRawMDMParameters(), giveRealStressVector(), initializeFrom(), and predictRelativeComputationalCost().

◆ ParMd

double oofem::MDM::ParMd = 0.

protected

THREAD UNSAFE. These are modified when evaluating the material. Why? This must be avoided at all cost.

(m/E*Ep)

Definition at line 174 of file mdm.h.

Referenced by computeDamageOnPlane(), giveInputRecord(), giveRawMDMParameters(), and initializeFrom().

◆ R

double oofem::MDM::R = 0.

protected

Interaction radius, related to the nonlocal characteristic length of material.

Definition at line 191 of file mdm.h.

Referenced by computeWeightFunction(), giveInputRecord(), giveRawMDMParameters(), and initializeFrom().

◆ sourceElemSet

std::unique_ptr<Set> oofem::MDM::sourceElemSet

protected

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

Definition at line 194 of file mdm.h.

Referenced by MMI_map().

◆ tempDillatCoeff

double oofem::MDM::tempDillatCoeff = 0.

protected

Temperature dilatation coeff.

Definition at line 175 of file mdm.h.

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

◆ type_dam

int oofem::MDM::type_dam = 0

protected

Definition at line 167 of file mdm.h.

Referenced by computeDamageOnPlane().

◆ type_soft

int oofem::MDM::type_soft = 0

protected

Definition at line 168 of file mdm.h.

Referenced by computeDamageOnPlane().

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

mdm.h
mdm.C

Public Types

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Static Protected Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ MDMFormulatrionType

◆ MDMMapperType

◆ MDMModeType

Constructor & Destructor Documentation

◆ MDM()

Member Function Documentation

◆ applyDamageToStiffness()

◆ applyDamageTranformation()

◆ computeDamageOnPlane()

◆ computeDamageTensor()

◆ computeEffectiveStress()

◆ computeLocalDamageTensor()

◆ computePDC()

◆ computeWeightFunction()

◆ CreateStatus()

◆ estimatePackSize()

◆ formTransformationMatrix()

◆ give3dMaterialStiffnessMatrix()

◆ giveClassName()

◆ giveInputRecord()

◆ giveInputRecordName()

◆ giveInterface()

◆ giveIPValue()

◆ giveMaterialStiffnessMatrix()

◆ givePlaneStrainStiffMtrx()

◆ givePlaneStressStiffMtrx()

◆ giveRawMDMParameters()

◆ giveRealStressVector()

◆ giveRealStressVector_3d()

◆ giveRealStressVector_PlaneStrain()

◆ giveRealStressVector_PlaneStress()

◆ giveSupportRadius()

◆ giveThermalDilatationVector()

◆ hasBoundedSupport()

◆ hasMaterialModeCapability()

◆ initializeData()

◆ initializeFrom()

◆ MMI_finish()

◆ MMI_map()

◆ MMI_update()

◆ packUnknowns()

◆ predictRelativeComputationalCost()

◆ predictRelativeRedistributionCost()

◆ rotateTensor4()

◆ transformStiffnessfromPDC()

◆ transformStrainToPDC()

◆ transformStressFromPDC()

◆ unpackAndUpdateUnknowns()

◆ updateBeforeNonlocAverage()

Member Data Documentation

◆ formulation

◆ Ft

◆ Gf

◆ linearElasticMaterial

◆ mapper2

◆ mapperLST

◆ mapperSFT

◆ mapperType

◆ mdm_Efp

◆ mdm_Ep

◆ mdmMode

◆ ndc

◆ nonlocal

◆ nsd

◆ ParMd

◆ R

◆ sourceElemSet

◆ tempDillatCoeff

◆ type_dam

◆ type_soft