OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
oofem::MDM Class Reference

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 InterfacegiveInterface (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 MaterialStatusCreateStatus (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...
 
MicroplanegiveMicroplane (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 IntegrationPointStatusgiveMicroplaneStatus (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 MaterialStatusgiveStatus (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...
 
DomaingiveDomain () 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...
 
DomaingiveDomain ()
 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

StructuralMaterialgiveLinearElasticMaterial ()
 Returns reference to undamaged (bulk) material. More...
 
virtual MaterialStatusCreateMicroplaneStatus (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
 
StructuralMateriallinearElasticMaterial
 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...
 
SetsourceElemSet
 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...
 
Domaindomain
 Link to domain object, useful for communicating with other FEM components. More...
 
- Protected Attributes inherited from oofem::NonlocalMaterialExtensionInterface
Domaindomain
 
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...
 
Gridgrid
 
FloatMatrixminDist2
 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

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.

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
nMaterial number.
dDomain 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::applyDamageTranformation ( FloatArray strainPDC,
const FloatArray tempDamageTensorEigenVals 
)
protected
void oofem::MDM::computeEffectiveStress ( FloatArray stressPDC,
const FloatArray strainPDC,
GaussPoint gp,
TimeStep tStep 
)
protected
void oofem::MDM::computePDC ( FloatMatrix tempDamageTensor,
FloatArray tempDamageTensorEigenVals,
FloatMatrix tempDamageTensorEigenVec 
)
protected
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
srcCoordinates of source point.
coordCoordinates 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
gpIntegration 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
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
answerComputed results.
modeMaterial response mode.
gpIntegration point.
tStepTime 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.

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
answerContain corresponding ip value, zero sized if not available.
gpIntegration point to which the value refers.
typeDetermines the type of internal variable.
tStepDetermines 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::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
answerStiffness matrix.
mmodeMaterial response mode.
gpIntegration point, which load history is used.
tStepTime 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
answerStiffness matrix.
mmodeMaterial response mode.
gpIntegration point, which load history is used.
tStepTime 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().

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
answerComputed result.
mplanePointer to microplane object, for which response is computed.
strainStrain vector.
tStepTime 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
answerStress vector in reduced form. For large deformations it is treated as the second Piola-Kirchoff stress.
gpIntegration point.
reducedStrainStrain vector in reduced form. For large deformations it is treated as the Green-Lagrange strain.
tStepCurrent 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
answerVector of thermal dilatation coefficients.
gpIntegration point.
tStepTime 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
modeRequired 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
numberOfMicroplanesNumber 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
irInput 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_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
gpIntegration point belonging to new domain which values will be mapped.
tStepTime step.
elemGPVecVector passed to MMI_update at material level, probably computed from primary unknowns (for structural elements this represent strain vector).
Returns
Nonzero if o.k.

Implements oofem::MaterialModelMapperInterface.

Definition at line 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
buffCommunication buffer.
tStepSolution step.
ipIntegration 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
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::transformStressFromPDC ( FloatArray answer,
const FloatArray stressPDC,
const FloatMatrix t,
GaussPoint gp 
)
protected
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
buffCommunication buffer.
tStepSolution step.
ipIntegration 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
strainVectorTotal strain vector in given integration point.
gpIntegration point to update.
tStepSolution step indicating time of update.

Implements oofem::StructuralNonlocalMaterialExtensionInterface.

Definition at line 1156 of file mdm.C.

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

Member Data Documentation

double oofem::MDM::Ft
protected

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

Definition at line 183 of file mdm.h.

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

double oofem::MDM::Gf
protected

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

Definition at line 181 of file mdm.h.

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

StructuralMaterial* oofem::MDM::linearElasticMaterial
protected

Reference to bulk (undamaged) material.

Definition at line 189 of file mdm.h.

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

MMAClosestIPTransfer oofem::MDM::mapper2
staticprotected

Mapper used to map stresses in adaptivity.

Definition at line 224 of file mdm.h.

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

MMALeastSquareProjection oofem::MDM::mapperLST
staticprotected

Definition at line 202 of file mdm.h.

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

MMAShapeFunctProjection oofem::MDM::mapperSFT
staticprotected

Mapper used to map internal variables in adaptivity.

Definition at line 201 of file mdm.h.

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

MDMMapperType oofem::MDM::mapperType
protected

Definition at line 206 of file mdm.h.

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

double oofem::MDM::mdm_Efp
protected

Prescribed value of ef-ep.

Definition at line 176 of file mdm.h.

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

double oofem::MDM::mdm_Ep
protected

Parameter controlling the elastic limit.

Definition at line 174 of file mdm.h.

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

int oofem::MDM::ndc
protected

Number of damage components.

Definition at line 168 of file mdm.h.

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

int oofem::MDM::nonlocal
protected

Flag indicating local or nonlocal mode.

Definition at line 192 of file mdm.h.

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

double oofem::MDM::ParMd
protected

(m/E*Ep)

Definition at line 177 of file mdm.h.

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

double oofem::MDM::R
protected

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

Definition at line 194 of file mdm.h.

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

Set* oofem::MDM::sourceElemSet
protected

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

Definition at line 197 of file mdm.h.

Referenced by MMI_map().

double oofem::MDM::tempDillatCoeff
protected

Temperature dilatation coeff.

Definition at line 178 of file mdm.h.

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

int oofem::MDM::type_dam
protected

Definition at line 170 of file mdm.h.

Referenced by computeDamageOnPlane().

int oofem::MDM::type_soft
protected

Definition at line 171 of file mdm.h.

Referenced by computeDamageOnPlane().


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

This page is part of the OOFEM documentation. Copyright (c) 2011 Borek Patzak
Project e-mail: info@oofem.org
Generated at Tue Jan 2 2018 20:07:38 for OOFEM by doxygen 1.8.11 written by Dimitri van Heesch, © 1997-2011