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

#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)
	Constructor. More...

virtual	~MDM ()
	Destructor. More...

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

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

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

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

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

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

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

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

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

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

virtual const char *	giveInputRecordName () const

virtual const char *	giveClassName () const

virtual void	giveRealMicroplaneStressVector (FloatArray &answer, Microplane mplane, const FloatArray &strain, TimeStep tStep)
	Computes real stress vector on given microplane (the meaning of values depends on particular implementation, e.g, can contain volumetric, deviatoric normal stresses and shear stresses on microplane) for given increment of microplane strains. More...

virtual void	initializeData (int numberOfMicroplanes)
	Initializes internal data (integration weights, microplane normals and computes projection tensors). More...

virtual void	updateBeforeNonlocAverage (const FloatArray &strainVector, GaussPoint gp, TimeStep tStep)
	Declares the service updating local variables in given integration points, which take part in nonlocal average process. More...

virtual double	computeWeightFunction (const FloatArray &src, const FloatArray &coord)
	Evaluates the basic nonlocal weight function for two points with given coordinates. More...

virtual int	hasBoundedSupport ()
	Determines, whether receiver has bounded weighting function (limited support). More...

virtual void	giveSupportRadius (double &radius)
	Determines the width (radius) of limited support of weighting function. More...

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

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

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

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

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

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

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

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

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

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

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

virtual	~MicroplaneMaterial ()
	Destructor. More...

double	computeNormalStrainComponent (Microplane *mplane, const FloatArray &macroStrain)
	Computes the length of normal strain vector on given microplane. More...

double	computeNormalVolumetricStrainComponent (Microplane *mplane, const FloatArray &macroStrain)
	Computes the normal volumetric component of macro strain on given microplane. More...

double	computeNormalDeviatoricStrainComponent (Microplane *mplane, const FloatArray &macroStrain)
	Computes the normal deviatoric component of macro strain on given microplane. More...

double	computeShearMStrainComponent (Microplane *mplane, const FloatArray &macroStrain)
	Computes the shear component (in m direction) of macro strain on given microplane. More...

double	computeShearLStrainComponent (Microplane *mplane, const FloatArray &macroStrain)
	Computes the shear component (in l direction) of macro strain on given microplane. More...

void	computeStrainVectorComponents (FloatArray &answer, Microplane *mplane, const FloatArray &macroStrain)
	Computes the vector of all micro stress components (Ev, En, Em, El) of macro strain vector on given microplane. More...

virtual void	giveMicroplaneNormal (FloatArray &answer, Microplane *mplane)
	Computes normal of given microplane. More...

virtual double	giveMicroplaneIntegrationWeight (Microplane *mplane)
	Returns microplane integration weight. More...

Microplane *	giveMicroplane (int i, GaussPoint *masterGp)
	Returns i-th microplane belonging to master-macro-integration point. )-based indexing. More...

virtual MaterialMode	giveCorrespondingSlaveMaterialMode (MaterialMode masterMode)
	Returns corresponding material mode for microplane according to macro integration mode. More...

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

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

virtual IntegrationPointStatus *	giveMicroplaneStatus (GaussPoint *gp)

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

virtual	~StructuralMaterial ()
	Destructor. More...

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

virtual	~Material ()
	Destructor. More...

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

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

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

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

double	giveCastingTime ()

virtual bool	isActivated (TimeStep *tStep)

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

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

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

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

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

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

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

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

Domain *	giveDomain () const

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

int	giveNumber () const

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

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

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

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

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

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

Public Member Functions inherited from oofem::StructuralNonlocalMaterialExtensionInterface
	StructuralNonlocalMaterialExtensionInterface (Domain *d)
	Constructor. More...

	~StructuralNonlocalMaterialExtensionInterface ()
	Destructor. More...

Public Member Functions inherited from oofem::NonlocalMaterialExtensionInterface
	NonlocalMaterialExtensionInterface (Domain *d)
	Constructor. More...

virtual	~NonlocalMaterialExtensionInterface ()
	Destructor. More...

void	updateDomainBeforeNonlocAverage (TimeStep *tStep)
	Updates data in all integration points before nonlocal average takes place. More...

void	buildNonlocalPointTable (GaussPoint *gp)
	Builds list of integration points which take part in nonlocal average in given integration point. More...

void	rebuildNonlocalPointTable (GaussPoint gp, IntArray contributingElems)
	Rebuild list of integration points which take part in nonlocal average in given integration point. More...

void	modifyNonlocalWeightFunctionAround (GaussPoint *gp)
	Recompute the nonlocal interaction weights based on the current solution (e.g., on the damage field). More...

void	modifyNonlocalWeightFunction_1D_Around (GaussPoint *gp)

virtual double	giveNonlocalMetricModifierAt (GaussPoint *gp)
	Provide the current value of the variable that affects nonlocal interaction (e.g., of damage) This method is used e.g. More...

double	computeDistanceModifier (double damage)
	Compute the factor that specifies how the interaction length should be modified, based on the current solution (e.g., on the damage field). More...

double	computeModifiedLength (double length, double dam1, double dam2)
	Compute the modified interaction length based on the current solution (e.g., on the damage field). More...

std::vector< localIntegrationRecord > *	giveIPIntegrationList (GaussPoint *gp)
	Returns integration list corresponding to given integration point. More...

virtual double	computeWeightFunction (double distance)
	Evaluates the basic nonlocal weight function for a given distance between interacting points. More...

double	giveIntegralOfWeightFunction (const int spatial_dimension)
	Provides the integral of the weight function over the contributing volume in 1, 2 or 3D. More...

virtual double	maxValueOfWeightFunction ()
	Determines the maximum value of the nonlocal weight function. More...

int	giveNumberOfRegions ()
	Determines the number of material regions of domain. More...

virtual double	evaluateSupportRadius ()
	Determines the width (radius) of limited support of weighting function. More...

Domain *	giveDomain ()
	Returns reference to domain. More...

IRResultType	initializeFrom (InputRecord *ir)

void	giveInputRecord (DynamicInputRecord &input)
	Stores receiver in an input record. More...

void	endIPNonlocalAverage (GaussPoint *gp)
	Notifies the receiver, that the nonlocal averaging has been finished for given ip. More...

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

virtual	~Interface ()

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

virtual	~MaterialModelMapperInterface ()
	Destructor. More...

Protected Types
enum	MDMMapperType { mdm_cpt =0, mdm_sft =1, mdm_lst =2 }

Protected Types inherited from oofem::NonlocalMaterialExtensionInterface
enum	WeightFunctionType { WFT_Unknown, WFT_Bell, WFT_Gauss, WFT_Green, WFT_Uniform, WFT_UniformOverElement, WFT_Green_21 }
	Type characterizing the nonlocal weight function. More...

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
StructuralMaterial *	giveLinearElasticMaterial ()
	Returns reference to undamaged (bulk) material. More...

virtual MaterialStatus *	CreateMicroplaneStatus (GaussPoint *gp)

void	computeDamageTensor (FloatMatrix &tempDamageTensor, const FloatArray &totalStrain, GaussPoint gp, TimeStep tStep)

void	computeLocalDamageTensor (FloatMatrix &tempDamageTensor, const FloatArray &totalStrain, GaussPoint gp, TimeStep tStep)

double	computeDamageOnPlane (GaussPoint gp, Microplane mplane, const FloatArray &strain)

void	computePDC (FloatMatrix &tempDamageTensor, FloatArray &tempDamageTensorEigenVals, FloatMatrix &tempDamageTensorEigenVec)

void	transformStrainToPDC (FloatArray &answer, FloatArray &strain, FloatMatrix &t, GaussPoint *gp)

void	applyDamageTranformation (FloatArray &strainPDC, const FloatArray &tempDamageTensorEigenVals)

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

void	computeEffectiveStress (FloatArray &stressPDC, const FloatArray &strainPDC, GaussPoint gp, TimeStep tStep)

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

void	applyDamageToStiffness (FloatMatrix &d, GaussPoint *gp)

void	transformStiffnessfromPDC (FloatMatrix &de, const FloatMatrix &t)

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

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

void	rotateTensor4 (FloatMatrix &Dlocal, const FloatMatrix &t)

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

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

Protected Member Functions inherited from oofem::MicroplaneMaterial
virtual void	initTempStatus (GaussPoint *gp)
	Initializes temporary variables stored in integration point status at the beginning of new time step. More...

Protected Member Functions inherited from oofem::NonlocalMaterialExtensionInterface
void	applyBarrierConstraints (const FloatArray &gpCoords, const FloatArray &jGpCoords, double &weight)

void	manipulateWeight (double &weight, GaussPoint gp, GaussPoint jGp)
	Manipulates weight on integration point in the element. More...

double	giveDistanceBasedInteractionRadius (const FloatArray &gpCoords)
	Provides the distance based interaction radius This function is called when nlvariation is set to 1. More...

int	mapToGridPoint (double x, double x0)

double	mapToGridCoord (double x, double x0)

double	dist2FromGridNode (double x, double y, int j, int i)

Protected Attributes
int	ndc
	Number of damage components. More...

int	nsd
	Number of spatial dimensions. More...

int	type_dam

int	type_soft

double	mdm_Ep
	Parameter controlling the elastic limit. More...

double	mdm_Efp
	Prescribed value of ef-ep. More...

double	ParMd
	(m/E*Ep) More...

double	tempDillatCoeff
	Temperature dilatation coeff. More...

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

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

MDMFormulatrionType	formulation

MDMModeType	mdmMode

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

int	nonlocal
	Flag indicating local or nonlocal mode. More...

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

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

MDMMapperType	mapperType

Protected Attributes inherited from oofem::MicroplaneMaterial
int	numberOfMicroplanes
	Number of microplanes. More...

double	microplaneWeights [MAX_NUMBER_OF_MICROPLANES]
	Integration weights of microplanes. More...

double	microplaneNormals [MAX_NUMBER_OF_MICROPLANES][3]
	Normals of microplanes. More...

double	Kronecker [6]
	Kronecker's delta. More...

double	N [MAX_NUMBER_OF_MICROPLANES][6]
	Normal projection tensors for all microplanes. More...

double	M [MAX_NUMBER_OF_MICROPLANES][6]
	Shear projection tensors (m direction) for all microplanes. More...

double	L [MAX_NUMBER_OF_MICROPLANES][6]
	Shear projection tensors (l direction) for all microplanes. More...

double	E
	Young's modulus. More...

double	nu
	Poisson's ratio. More...

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

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

double	castingTime
	Casting time. More...

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

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

Protected Attributes inherited from oofem::NonlocalMaterialExtensionInterface
Domain *	domain

IntArray	regionMap
	Map indicating regions to skip (region - cross section model). More...

bool	permanentNonlocTableFlag
	Flag indicating whether to keep nonlocal interaction tables of integration points cached. More...

WeightFunctionType	weightFun
	Parameter specifying the type of nonlocal weight function. More...

int	gridSize
	Grid on which the eikonal equation will be solved (used by eikonal nonlocal models) More...

Grid *	grid

FloatMatrix *	minDist2
	Auxiliary matrix to store minimum distances of grid points from Gauss points. More...

double	initDiag
	Optional parameters setting details of the fast marching method. More...

int	order

int	centDiff

double	cl
	Characteristic length of the nonlocal model (its interpretation depends on the type of weight function). More...

double	suprad
	Support radius. More...

double	mm
	For "undernonlocal" or "overnonlocal" formulation. More...

ScalingType	scaling
	Parameter specifying the type of scaling of nonlocal weight function. More...

AveragedVarType	averagedVar
	Parameter specifying the type of averaged (nonlocal) variable. More...

double	cl0
	Initial(user defined) characteristic length of the nonlocal model (its interpretation depends on the weight function) Is different to cl when a Stress-based or a Distance-based nonlocal variation is applied. More...

NlVariationType	nlvar
	Parameter specifying the type of nonlocal variation. More...

double	beta
	Parameter which multiplied with the interaction radius cl0 gives its minimum allowed value. More...

double	zeta
	Parameter used when Distance-based nonlocal variation is applied When it is multiplied with the interaction radius cl gives the maxinmum distance of the Gauss Point from the boundary. More...

double	px
	Parameter specifying the periodic shift in x-direction. More...

double	Rf
	Final value of interaction radius, for a model with evolving characteristic length. More...

double	exponent
	Parameter used as an exponent by models with evolving characteristic length. More...

int	averType
	Parameter specifying how the weight function should be adjusted due to damage. More...

Static Protected Attributes
static MMAShapeFunctProjection	mapperSFT
	Mapper used to map internal variables in adaptivity. More...

static MMALeastSquareProjection	mapperLST

static MMAClosestIPTransfer	mapper2
	Mapper used to map stresses in adaptivity. More...

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

static double	computeStressNorm (const FloatArray &stress)

static double	computeFirstInvariant (const FloatArray &s)

static double	computeSecondStressInvariant (const FloatArray &s)

static double	computeThirdStressInvariant (const FloatArray &s)

static double	computeFirstCoordinate (const FloatArray &s)

static double	computeSecondCoordinate (const FloatArray &s)

static double	computeThirdCoordinate (const FloatArray &s)

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

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

Detailed Description

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

Definition at line 161 of file mdm.h.

Member Enumeration Documentation

enum oofem::MDM::MDMFormulatrionType

Enumerator
COMPLIANCE_DAMAGE
STIFFNESS_DAMAGE

Definition at line 164 of file mdm.h.

enum oofem::MDM::MDMMapperType

protected

Enumerator
mdm_cpt
mdm_sft
mdm_lst

Definition at line 205 of file mdm.h.

enum oofem::MDM::MDMModeType

Enumerator
mdm_3d
mdm_2d

Definition at line 165 of file mdm.h.

Constructor & Destructor Documentation

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 232 of file mdm.h.

virtual oofem::MDM::~MDM ( )

inlinevirtual

Destructor.

Definition at line 242 of file mdm.h.

References oofem::IntegrationPointStatus::gp.

Member Function Documentation

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

protected

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

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

protected

Definition at line 416 of file mdm.C.

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

Referenced by giveRealStressVector().

double oofem::MDM::computeDamageOnPlane	(	GaussPoint *	gp,
		Microplane *	mplane,
		const FloatArray &	strain
	)

protected

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

protected

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

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

protected

Definition at line 465 of file mdm.C.

References oofem::FloatArray::beProductOf(), oofem::StructuralMaterial::give3dMaterialStiffnessMatrix(), giveLinearElasticMaterial(), oofem::StructuralMaterial::giveStiffnessMatrix(), mdm_3d, and mdmMode.

Referenced by giveRealStressVector().

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

protected

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

protected

Definition at line 349 of file mdm.C.

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

Referenced by giveRealStressVector().

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

virtual

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

References oofem::FloatArray::distance(), and R.

virtual MaterialStatus* oofem::MDM::CreateMicroplaneStatus ( GaussPoint * gp )

inlineprotectedvirtual

Implements oofem::MicroplaneMaterial.

Definition at line 313 of file mdm.h.

MaterialStatus * oofem::MDM::CreateStatus ( GaussPoint * gp ) const

virtual

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

References oofem::FEMComponent::giveDomain(), nsd, and oofem::MicroplaneMaterial::numberOfMicroplanes.

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

virtual

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

Reimplemented from oofem::Material.

Definition at line 1342 of file mdm.C.

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

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

protected

Definition at line 943 of file mdm.C.

References oofem::FloatMatrix::at(), FORMT33, FORMT44, FORMT66, OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by rotateTensor4().

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

virtual

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

Definition at line 641 of file mdm.C.

References giveMaterialStiffnessMatrix().

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

inlinevirtual

Implements oofem::Interface.

Definition at line 278 of file mdm.h.

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

virtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::MicroplaneMaterial.

Definition at line 877 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::MMALeastSquareProjection::giveInputRecord(), oofem::StructuralMaterial::giveInputRecord(), oofem::MaterialMappingAlgorithm::giveInputRecord(), oofem::MicroplaneMaterial::giveInputRecord(), oofem::NonlocalMaterialExtensionInterface::giveInputRecord(), linearElasticMaterial, mapper2, mapperLST, mapperSFT, mapperType, mdm_Efp, mdm_Ep, mdmMode, nonlocal, ParMd, R, oofem::DynamicInputRecord::setField(), and tempDillatCoeff.

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

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 277 of file mdm.h.

References _IFT_MDM_Name.

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

virtual

Interface requesting service.

See also: InterfaceType

Returns: Requested interface if implemented, otherwise NULL.

Reimplemented from oofem::FEMComponent.

Definition at line 1307 of file mdm.C.

References oofem::MaterialModelMapperInterfaceType, and oofem::NonlocalMaterialExtensionInterfaceType.

int oofem::MDM::giveIPValue	(	FloatArray &	answer,
		GaussPoint *	gp,
		InternalStateType	type,
		TimeStep *	tStep
	)

virtual

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

Definition at line 667 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::StructuralMaterial::giveIPValue(), oofem::MDMStatus::giveMicroplaneDamageValues(), oofem::Material::giveStatus(), oofem::MDMStatus::giveTempDamageTensor(), oofem::MDMStatus::giveTempDamageTensorEigenVals(), ndc, nsd, oofem::FloatArray::resize(), and oofem::FloatArray::zero().

StructuralMaterial* oofem::MDM::giveLinearElasticMaterial ( )

inlineprotected

Returns reference to undamaged (bulk) material.

Definition at line 311 of file mdm.h.

Referenced by computeEffectiveStress(), giveMaterialStiffnessMatrix(), and giveRawMDMParameters().

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

protected

Definition at line 525 of file mdm.C.

References applyDamageToStiffness(), giveLinearElasticMaterial(), oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStiffnessMatrix(), oofem::MDMStatus::giveTempDamageTensorEigenVec(), and transformStiffnessfromPDC().

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

void oofem::MDM::givePlaneStrainStiffMtrx	(	FloatMatrix &	answer,
		MatResponseMode	mmode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

protectedvirtual

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

References giveMaterialStiffnessMatrix().

void oofem::MDM::givePlaneStressStiffMtrx	(	FloatMatrix &	answer,
		MatResponseMode	mmode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

protectedvirtual

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

References giveMaterialStiffnessMatrix().

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

protected

virtual void oofem::MDM::giveRealMicroplaneStressVector	(	FloatArray &	answer,
		Microplane *	mplane,
		const FloatArray &	strain,
		TimeStep *	tStep
	)

inlinevirtual

Computes real stress vector on given microplane (the meaning of values depends on particular implementation, e.g, can contain volumetric, deviatoric normal stresses and shear stresses on microplane) for given increment of microplane strains.

Parameters

answer	Computed result.
mplane	Pointer to microplane object, for which response is computed.
strain	Strain vector.
tStep	Time step.

Implements oofem::MicroplaneMaterial.

Definition at line 280 of file mdm.h.

References oofem::StructuralMaterialStatus::strainVector.

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

virtual

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 95 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 MMI_update().

virtual void oofem::MDM::giveRealStressVector_3d	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedE,
		TimeStep *	tStep
	)

inlinevirtual

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

Reimplemented from oofem::StructuralMaterial.

Definition at line 256 of file mdm.h.

virtual void oofem::MDM::giveRealStressVector_PlaneStrain	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedE,
		TimeStep *	tStep
	)

inlinevirtual

Default implementation relies on giveRealStressVector_3d.

Reimplemented from oofem::StructuralMaterial.

Definition at line 258 of file mdm.h.

virtual void oofem::MDM::giveRealStressVector_PlaneStress	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	reducedE,
		TimeStep *	tStep
	)

inlinevirtual

Default implementation relies on giveRealStressVector_StressControl.

Reimplemented from oofem::StructuralMaterial.

Definition at line 260 of file mdm.h.

References oofem::FEMComponent::giveInputRecord(), and oofem::MaterialStatus::initializeFrom().

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

inlinevirtual

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

Definition at line 293 of file mdm.h.

References oofem::MDMStatus::giveInterface().

void oofem::MDM::giveThermalDilatationVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

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

References oofem::FloatArray::at(), oofem::FloatArray::resize(), tempDillatCoeff, and oofem::FloatArray::zero().

virtual int oofem::MDM::hasBoundedSupport ( )

inlinevirtual

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

Returns: True if weighting function bounded, zero otherwise.

Reimplemented from oofem::NonlocalMaterialExtensionInterface.

Definition at line 289 of file mdm.h.

int oofem::MDM::hasMaterialModeCapability ( MaterialMode mode )

virtual

Tests if material supports material mode.

Parameters

mode	Required material mode.

Returns: Nonzero if supported, zero otherwise.

Reimplemented from oofem::StructuralMaterial.

Definition at line 86 of file mdm.C.

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

virtual

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

Parameters

numberOfMicroplanes Number of required microplanes.

Reimplemented from oofem::MicroplaneMaterial.

Definition at line 1105 of file mdm.C.

References oofem::FloatArray::at(), oofem::MicroplaneMaterial::initializeData(), 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.

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

virtual

Initializes receiver according to object description stored in input record.

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

See also: IR_GIVE_FIELD; IR_GIVE_OPTIONAL_FIELD

Parameters

ir	Input record to initialize from.

Returns: IRResultType

Reimplemented from oofem::MicroplaneMaterial.

Definition at line 783 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::FEMComponent::giveDomain(), oofem::InputRecord::hasField(), oofem::MMALeastSquareProjection::initializeFrom(), oofem::StructuralMaterial::initializeFrom(), oofem::MaterialMappingAlgorithm::initializeFrom(), oofem::MicroplaneMaterial::initializeFrom(), oofem::NonlocalMaterialExtensionInterface::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_BAD_FORMAT, oofem::IRRT_OK, linearElasticMaterial, mapper2, mapperLST, mapperSFT, mapperType, mdm_2d, mdm_3d, mdm_Efp, mdm_Ep, mdmMode, ndc, nonlocal, nsd, OOFEM_LOG_INFO, OOFEM_WARNING, ParMd, R, and tempDillatCoeff.

int oofem::MDM::MMI_finish ( TimeStep * tStep )

virtual

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

References oofem::MMAClosestIPTransfer::finish(), oofem::MMAShapeFunctProjection::finish(), oofem::MMALeastSquareProjection::finish(), mapper2, mapperLST, and mapperSFT.

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

virtual

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

References oofem::IntArray::at(), oofem::FloatArray::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::MaterialMappingAlgorithm::init(), oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), mapper2, mapperLST, mapperSFT, mapperType, oofem::MMAShapeFunctProjection::mapVariable(), oofem::MaterialMappingAlgorithm::mapVariable(), mdm_cpt, mdm_lst, mdm_sft, OOFEM_ERROR, oofem::IntArray::resize(), oofem::Set::setElementList(), oofem::MDMStatus::setMicroplaneTempDamageValues(), sourceElemSet, and oofem::MDMStatus::updateYourself().

int oofem::MDM::MMI_update	(	GaussPoint *	gp,
		TimeStep *	tStep,
		FloatArray *	elemGPVec = `NULL`
	)

virtual

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

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

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

virtual

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

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

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

virtual

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

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

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

inlinevirtual

Returns the relative redistribution cost of the receiver.

Reimplemented from oofem::Material.

Definition at line 305 of file mdm.h.

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

protected

Definition at line 923 of file mdm.C.

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

Referenced by transformStiffnessfromPDC().

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

protected

Definition at line 634 of file mdm.C.

References rotateTensor4().

Referenced by giveMaterialStiffnessMatrix().

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

protected

Definition at line 371 of file mdm.C.

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

Referenced by giveRealStressVector().

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

protected

Definition at line 484 of file mdm.C.

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

Referenced by giveRealStressVector().

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

virtual

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

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

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

virtual

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