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

This class implements nonassociated multisurface plasticity model. More...

#include <dustmat.h>

+ Inheritance diagram for oofem::DustMaterial:
+ Collaboration diagram for oofem::DustMaterial:

Public Member Functions

 DustMaterial (int n, Domain *d)
 Constructor. More...
 
virtual ~DustMaterial ()
 Destructor. More...
 
virtual IRResultType initializeFrom (InputRecord *ir)
 Initializes receiver according to object description stored in input record. More...
 
virtual const char * giveClassName () const
 
virtual const char * giveInputRecordName () const
 
virtual void giveRealStressVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &strainVector, TimeStep *tStep)
 Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. More...
 
virtual void give3dMaterialStiffnessMatrix (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point. More...
 
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 int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep)
 Returns the integration point corresponding value in Reduced form. More...
 
virtual bool isCharacteristicMtrxSymmetric (MatResponseMode rMode)
 Returns true if stiffness matrix of receiver is symmetric Default implementation returns true. 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 MaterialStatusCreateStatus (GaussPoint *gp) const
 Creates new copy of associated status and inserts it into given integration point. More...
 
double giveQ0 () const
 
- Public Member Functions inherited from oofem::StructuralMaterial
 StructuralMaterial (int n, Domain *d)
 Constructor. More...
 
virtual ~StructuralMaterial ()
 Destructor. More...
 
virtual int hasMaterialModeCapability (MaterialMode mode)
 Tests if material supports material mode. More...
 
virtual void giveInputRecord (DynamicInputRecord &input)
 Setups the input record string of receiver. 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 (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_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector_3d. 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_PlaneStress (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep)
 Default implementation relies on giveRealStressVector_StressControl. More...
 
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 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 (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing plane stress stiffness matrix of receiver. More...
 
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 (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
 Method for computing plane strain stiffness matrix of receiver. More...
 
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 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 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 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...
 
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 void initTempStatus (GaussPoint *gp)
 Initializes temporary variables stored in integration point status at the beginning of new time step. 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...
 
virtual InterfacegiveInterface (InterfaceType t)
 Interface requesting service. More...
 
std::string errorInfo (const char *func) const
 Returns string for prepending output (used by error reporting macros). More...
 

Protected Member Functions

double functionFe (double i1)
 Auxiliary equation Fe (7.8) More...
 
double functionFeDI1 (double i1)
 Derivative by i1 of auxiliary equation (7.8) More...
 
double functionFeDI1DI1 (double i1)
 Second derivative by i1 of auxiliary equation (7.8) More...
 
double functionFc (double rho, double i1, double q)
 Auxiliary equation Fc (7.9) More...
 
double functionX (double q)
 Auxiliary equation X (7.11) More...
 
double functionXDQ (double q)
 Derivative by q of auxiliary equation X (7.11) More...
 
double yieldFunction1 (double rho, double i1)
 Yield function 1 (shear dominant), equation 7.5. More...
 
double yieldFunction2 (double rho, double i1, double q)
 Yield function 2 (compression dominant), equation 7.6. More...
 
double yieldFunction3 (double i1)
 Yield function 3 (tension dominant), equation 7.7. More...
 
void solveQ0 (double &answer)
 Solves q0 according to given parameters, equation 7.12. More...
 
void computeAndSetBulkAndShearModuli (double &bulkModulus, double &shearModulus, GaussPoint *gp)
 Computes and sets all elastic moduli, with possible stiffening. More...
 
void performStressReturn (GaussPoint *gp, const FloatArray &strain)
 Perform stress return and update all internal variables. More...
 
void computePlastStrainDirM1 (FloatArray &answer, const FloatArray &stressDeviator, double rho, double i1, double q)
 Computes direction of plastic yielding m1, equation 7.17. More...
 
void computePlastStrainDirM2 (FloatArray &answer, const FloatArray &stressDeviator, double rho, double i1, double q)
 Computes direction of plastic yielding m2, equation 7.19. More...
 
void computePlastStrainDirM3 (FloatArray &answer, const FloatArray &stressDeviator, double rho, double i1, double q)
 Computes direction of plastic yielding m2, equation 7.18. More...
 
double functionH (double q, double tempQ)
 Auxiliary equation H (7.33 or 7.34) More...
 
double functionHDQ (double tempQ)
 Derivative by tempQ of auxiliary equation H (7.33 or 7.34) More...
 
double functionI1 (double q, double tempQ, double i1, double bulkModulus)
 Auxiliary equation I1 (7.32) More...
 
double functionI1DQ (double tempQ, double bulkModulus)
 Derivative by tempQ of auxiliary equation I1 (7.32) More...
 
void performF1return (double i1, double rho, GaussPoint *gp)
 Performs stress return of case of yield function F1, computes new value of tempQ and sets it to status. More...
 
void performF2return (double i1, double rho, GaussPoint *gp)
 Performs stress return of case of yield function F2, computes new value of tempQ and sets it to status. More...
 
void computeQFromPlastVolEps (double &answer, double q, double deltaVolumetricPlasticStrain)
 Computes tempQ from volumetric plastic strain increment, equation 7.44. More...
 
double computeDeltaGamma2 (double tempQ, double q, double i1, double bulkModulus)
 Computed value of plastic multiplier for F2 yield function, equation 7.39. More...
 
double computeDeltaGamma2DQ (double tempQ, double q, double i1, double bulkModulus)
 Computed derivative by tempQ of equation 7.39. More...
 
double fTempR2 (double tempQ, double q, double i1, double rho, double bulkModulus, double shearModulus)
 Equation 7.38. More...
 

Protected Attributes

IsotropicLinearElasticMaterialLEMaterial
 Pointer for linear elastic material. More...
 
double alpha
 Parameter determining shape of yield surface. More...
 
double beta
 Parameter determining shape of yield surface. More...
 
double lambda
 Parameter determining shape of yield surface. More...
 
double theta
 Parameter determining shape of yield surface. More...
 
double ft
 Parameter determining shape of yield surface (param T in original publication) More...
 
double rEll
 Parameter determining shape of yield surface (param R in original publication) More...
 
int hardeningType
 Parameter determining hardening type. More...
 
double mStiff
 Parameter increasing stiffness (parameter M in original publication) More...
 
double x0
 Parameter determining shape of yield surface (param X0 in original publication) More...
 
double q0
 Parameter determining shape of yield surface. More...
 
double wHard
 Parameter determining hardening law (parameter W in original publication) More...
 
double dHard
 Parameter determining hardening law (parameter D in original publication) More...
 
double newtonTol
 Tollerance for iterative methods. More...
 
int newtonIter
 Maximum number of iterations for iterative methods. 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...
 

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

This class implements nonassociated multisurface plasticity model.

Yield function depends on I1 and J2 invariants of stress tensor and internal variable q. The plasticity surface consists of three surfaces. Shear dominant surface (f1) is non linear Drucker Prager cone. Compressive dominant part (f2) is elliptic cap over the cone. Tension dominant part (f3) is plane, cutting the cone perpendicularly to hydrostatic axis. f1 and f3 are constant. The position of f2 surface is determined by internal variable q. During plastic yielding f1 or f3 (shear or tension dominant state), q is increasing and the admissible elastic domain is shrinking. During f2 plastic yielding (compressive dominant), the parameter q is decreasing and the admissible elastic domain is becoming larger.

This model is described in PhD thesis of Tobias Erhard Strategien zur numerischen Modellierung transierten Impaktvorgange bei nichtlinearem Materialverhalten Universitat Stuttgart 2004

Author
Jan Stransky

Definition at line 233 of file dustmat.h.

Constructor & Destructor Documentation

oofem::DustMaterial::DustMaterial ( int  n,
Domain d 
)

Constructor.

Definition at line 158 of file dustmat.C.

References LEMaterial.

oofem::DustMaterial::~DustMaterial ( )
virtual

Destructor.

Definition at line 163 of file dustmat.C.

References LEMaterial.

Member Function Documentation

void oofem::DustMaterial::computeAndSetBulkAndShearModuli ( double &  bulkModulus,
double &  shearModulus,
GaussPoint gp 
)
protected
double oofem::DustMaterial::computeDeltaGamma2 ( double  tempQ,
double  q,
double  i1,
double  bulkModulus 
)
protected

Computed value of plastic multiplier for F2 yield function, equation 7.39.

Parameters
tempQParameter tempQ
qParameter q from previous step
i1Trace of stress tensor
bulkModulusBulk modulus

Definition at line 739 of file dustmat.C.

References functionFe(), functionH(), and rEll.

Referenced by fTempR2(), and performStressReturn().

double oofem::DustMaterial::computeDeltaGamma2DQ ( double  tempQ,
double  q,
double  i1,
double  bulkModulus 
)
protected

Computed derivative by tempQ of equation 7.39.

Parameters
tempQParameter tempQ
qParameter q from previous step
i1Trace of stress tensor
bulkModulusBulk modulus

Definition at line 746 of file dustmat.C.

References functionFe(), functionFeDI1(), functionH(), functionHDQ(), and rEll.

void oofem::DustMaterial::computePlastStrainDirM1 ( FloatArray answer,
const FloatArray stressDeviator,
double  rho,
double  i1,
double  q 
)
protected

Computes direction of plastic yielding m1, equation 7.17.

Parameters
answerResult
stressDeviatorDeviator of stress tensor
rhoSecond Haigh-Westergaard coordinate
i1Trace of stress tensor
qParameter q

Definition at line 663 of file dustmat.C.

References oofem::FloatArray::at(), oofem::FloatArray::beScaled(), beta, ft, lambda, and theta.

Referenced by performStressReturn().

void oofem::DustMaterial::computePlastStrainDirM2 ( FloatArray answer,
const FloatArray stressDeviator,
double  rho,
double  i1,
double  q 
)
protected

Computes direction of plastic yielding m2, equation 7.19.

Parameters
answerResult
stressDeviatorDeviator of stress tensor
rhoSecond Haigh-Westergaard coordinate
i1Trace of stress tensor
qParameter q

Definition at line 674 of file dustmat.C.

References oofem::FloatArray::at(), oofem::FloatArray::beScaled(), functionFc(), and rEll.

Referenced by performStressReturn().

void oofem::DustMaterial::computePlastStrainDirM3 ( FloatArray answer,
const FloatArray stressDeviator,
double  rho,
double  i1,
double  q 
)
protected

Computes direction of plastic yielding m2, equation 7.18.

Parameters
answerResult
stressDeviatorDeviator of stress tensor
rhoSecond Haigh-Westergaard coordinate
i1Trace of stress tensor
qParameter q

Definition at line 686 of file dustmat.C.

References oofem::FloatArray::at(), oofem::FloatArray::beScaled(), ft, functionFe(), and functionFeDI1().

Referenced by performStressReturn().

void oofem::DustMaterial::computeQFromPlastVolEps ( double &  answer,
double  q,
double  deltaVolumetricPlasticStrain 
)
protected

Computes tempQ from volumetric plastic strain increment, equation 7.44.

Parameters
answerResult tempQ
qParameter q from previous step
deltaVolumetricPlasticStrainVolumetric plastic strain increment

Definition at line 476 of file dustmat.C.

References functionH(), functionHDQ(), newtonIter, newtonTol, OOFEM_ERROR, and OOFEM_LOG_DEBUG.

Referenced by performStressReturn().

MaterialStatus * oofem::DustMaterial::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 565 of file dustmat.C.

References oofem::FEMComponent::giveDomain(), and giveQ0().

double oofem::DustMaterial::fTempR2 ( double  tempQ,
double  q,
double  i1,
double  rho,
double  bulkModulus,
double  shearModulus 
)
protected

Equation 7.38.

Parameters
tempQParameter tempQ
qParameter q from previous step
i1Trace of stress tensor
rhoSecond Haigh-Westergaard coordinate
bulkModulusBulk modulus
shearModulusShear modulus

Definition at line 760 of file dustmat.C.

References computeDeltaGamma2(), functionFe(), and rEll.

Referenced by performF2return().

double oofem::DustMaterial::functionFc ( double  rho,
double  i1,
double  q 
)
protected

Auxiliary equation Fc (7.9)

Parameters
i1Trace of stress tensor
rhoSecond Haigh-Westergaard coordinate
qParameter q
Returns
Fc

Definition at line 589 of file dustmat.C.

References rEll.

Referenced by computePlastStrainDirM2(), and yieldFunction2().

double oofem::DustMaterial::functionFe ( double  i1)
protected

Auxiliary equation Fe (7.8)

Parameters
i1Trace of stress tensor
Returns
Fe

Definition at line 571 of file dustmat.C.

References alpha, beta, lambda, and theta.

Referenced by computeDeltaGamma2(), computeDeltaGamma2DQ(), computePlastStrainDirM3(), fTempR2(), performF1return(), performStressReturn(), yieldFunction1(), and yieldFunction2().

double oofem::DustMaterial::functionFeDI1 ( double  i1)
protected

Derivative by i1 of auxiliary equation (7.8)

Parameters
i1Trace of stress tensor
Returns
$ \frac{\partial Fe}{\partial i_1 } $

Definition at line 577 of file dustmat.C.

References beta, lambda, and theta.

Referenced by computeDeltaGamma2DQ(), computePlastStrainDirM3(), performF1return(), and performStressReturn().

double oofem::DustMaterial::functionFeDI1DI1 ( double  i1)
protected

Second derivative by i1 of auxiliary equation (7.8)

Parameters
i1Trace of stress tensor
Returns
$ \frac{\partial^2 Fe}{\partial i1^2 } $

Definition at line 583 of file dustmat.C.

References beta, and lambda.

Referenced by performF1return().

double oofem::DustMaterial::functionH ( double  q,
double  tempQ 
)
protected

Auxiliary equation H (7.33 or 7.34)

Parameters
qParameter q from previous step
tempQParameter tempQ
Returns
H

Definition at line 699 of file dustmat.C.

References dHard, functionX(), hardeningType, and wHard.

Referenced by computeDeltaGamma2(), computeDeltaGamma2DQ(), computeQFromPlastVolEps(), functionI1(), performF1return(), and performF2return().

double oofem::DustMaterial::functionHDQ ( double  tempQ)
protected

Derivative by tempQ of auxiliary equation H (7.33 or 7.34)

Parameters
tempQParameter tempQ
Returns
$ \frac{\partial X}{\partial tempQ} $

Definition at line 713 of file dustmat.C.

References dHard, functionX(), functionXDQ(), hardeningType, and wHard.

Referenced by computeDeltaGamma2DQ(), computeQFromPlastVolEps(), functionI1DQ(), performF1return(), and performF2return().

double oofem::DustMaterial::functionI1 ( double  q,
double  tempQ,
double  i1,
double  bulkModulus 
)
protected

Auxiliary equation I1 (7.32)

Parameters
qParameter q from previous step
tempQParameter tempQ
i1Trace of stress tensor
bulkModulusBulk modulus
Returns
I1

Definition at line 727 of file dustmat.C.

References functionH().

Referenced by performF1return(), and performStressReturn().

double oofem::DustMaterial::functionI1DQ ( double  tempQ,
double  bulkModulus 
)
protected

Derivative by tempQ of auxiliary equation I1 (7.32)

Parameters
tempQParameter tempQ
bulkModulusBulk modulus
Returns
$ \frac{\partial I1}{\partial tempQ} $

Definition at line 733 of file dustmat.C.

References functionHDQ().

Referenced by performF1return().

double oofem::DustMaterial::functionX ( double  q)
protected

Auxiliary equation X (7.11)

Parameters
qParameter q
Returns
X

Definition at line 613 of file dustmat.C.

References alpha, beta, lambda, rEll, and theta.

Referenced by functionH(), and functionHDQ().

double oofem::DustMaterial::functionXDQ ( double  q)
protected

Derivative by q of auxiliary equation X (7.11)

Parameters
qParameter q
Returns
$ \frac{\partial X }{\partial q} $

Definition at line 619 of file dustmat.C.

References beta, lambda, rEll, and theta.

Referenced by functionHDQ().

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

Definition at line 502 of file dustmat.C.

References oofem::IsotropicLinearElasticMaterial::give3dMaterialStiffnessMatrix(), oofem::Material::giveStatus(), oofem::DustMaterialStatus::giveVolumetricPlasticStrain(), oofem::IsotropicLinearElasticMaterial::giveYoungsModulus(), oofem::DustMaterialStatus::giveYoungsModulus(), LEMaterial, OOFEM_ERROR, and oofem::FloatMatrix::times().

virtual const char* oofem::DustMaterial::giveClassName ( ) const
inlinevirtual
Returns
Class name of the receiver.

Reimplemented from oofem::StructuralMaterial.

Definition at line 457 of file dustmat.h.

virtual const char* oofem::DustMaterial::giveInputRecordName ( ) const
inlinevirtual
Returns
Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 458 of file dustmat.h.

References _IFT_DustMaterial_Name, oofem::IntegrationPointStatus::gp, and oofem::StructuralMaterialStatus::strainVector.

int oofem::DustMaterial::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.
Todo:
This is actually the mean, not the volumetric part. / Mikael

Reimplemented from oofem::StructuralMaterial.

Definition at line 537 of file dustmat.C.

References oofem::FloatArray::at(), oofem::StructuralMaterial::computePrincipalValues(), oofem::StructuralMaterial::giveIPValue(), oofem::DustMaterialStatus::givePlasticStrain(), oofem::DustMaterialStatus::giveQ(), oofem::Material::giveStatus(), oofem::DustMaterialStatus::giveVolumetricPlasticStrain(), oofem::principal_strain, and oofem::FloatArray::resize().

double oofem::DustMaterial::giveQ0 ( ) const
inline

Definition at line 483 of file dustmat.h.

Referenced by CreateStatus().

void oofem::DustMaterial::giveRealStressVector_3d ( FloatArray answer,
GaussPoint gp,
const FloatArray reducedE,
TimeStep tStep 
)
virtual
virtual void oofem::DustMaterial::giveThermalDilatationVector ( FloatArray answer,
GaussPoint gp,
TimeStep tStep 
)
inlinevirtual

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 476 of file dustmat.h.

References oofem::IsotropicLinearElasticMaterial::giveThermalDilatationVector().

IRResultType oofem::DustMaterial::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::StructuralMaterial.

Definition at line 169 of file dustmat.C.

References _IFT_DustMaterial_alpha, _IFT_DustMaterial_beta, _IFT_DustMaterial_dHard, _IFT_DustMaterial_ft, _IFT_DustMaterial_hardeningType, _IFT_DustMaterial_lambda, _IFT_DustMaterial_mStiff, _IFT_DustMaterial_newtonIter, _IFT_DustMaterial_newtonTol, _IFT_DustMaterial_rEll, _IFT_DustMaterial_theta, _IFT_DustMaterial_wHard, _IFT_DustMaterial_x0, alpha, beta, dHard, ft, hardeningType, oofem::IsotropicLinearElasticMaterial::initializeFrom(), oofem::StructuralMaterial::initializeFrom(), IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_BAD_FORMAT, oofem::IRRT_OK, lambda, LEMaterial, mStiff, newtonIter, newtonTol, OOFEM_WARNING, q0, rEll, solveQ0(), theta, wHard, and x0.

virtual bool oofem::DustMaterial::isCharacteristicMtrxSymmetric ( MatResponseMode  rMode)
inlinevirtual

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

Reimplemented from oofem::Material.

Definition at line 474 of file dustmat.h.

void oofem::DustMaterial::performF1return ( double  i1,
double  rho,
GaussPoint gp 
)
protected

Performs stress return of case of yield function F1, computes new value of tempQ and sets it to status.

Equation 7.31

Parameters
i1Trace of stress tensor
rhoSecond Haigh-Westergaard coordinate
gpGauss point

Definition at line 390 of file dustmat.C.

References ft, functionFe(), functionFeDI1(), functionFeDI1DI1(), functionH(), functionHDQ(), functionI1(), functionI1DQ(), oofem::DustMaterialStatus::giveBulkModulus(), oofem::DustMaterialStatus::giveQ(), oofem::DustMaterialStatus::giveShearModulus(), oofem::Material::giveStatus(), oofem::DustMaterialStatus::giveTempQ(), oofem::DustMaterialStatus::letTempQBe(), newtonIter, newtonTol, OOFEM_ERROR, and OOFEM_LOG_DEBUG.

Referenced by performStressReturn().

void oofem::DustMaterial::performF2return ( double  i1,
double  rho,
GaussPoint gp 
)
protected

Performs stress return of case of yield function F2, computes new value of tempQ and sets it to status.

Equation 7.38

Parameters
i1Trace of stress tensor
rhoSecond Haigh-Westergaard coordinate
gpGauss point

Definition at line 436 of file dustmat.C.

References fTempR2(), functionH(), functionHDQ(), oofem::DustMaterialStatus::giveBulkModulus(), oofem::DustMaterialStatus::giveQ(), oofem::DustMaterialStatus::giveShearModulus(), oofem::Material::giveStatus(), oofem::DustMaterialStatus::letTempQBe(), newtonIter, newtonTol, OOFEM_ERROR, and q0.

Referenced by performStressReturn().

int oofem::DustMaterial::setIPValue ( const FloatArray value,
GaussPoint gp,
InternalStateType  type 
)
virtual

Sets the value of a certain variable at a given integration point to the given value.

Parameters
valueContains the value(s) to be set (in reduced form).
gpIntegration point.
typeDetermines the type of internal variable.
typeDetermines the type of internal variable.
Returns
Nonzero if ok, zero if var not supported.

Reimplemented from oofem::StructuralMaterial.

Definition at line 522 of file dustmat.C.

References oofem::FloatArray::at(), oofem::Material::giveStatus(), oofem::DustMaterialStatus::letPlasticStrainBe(), oofem::DustMaterialStatus::letQBe(), and oofem::StructuralMaterial::setIPValue().

void oofem::DustMaterial::solveQ0 ( double &  answer)
protected

Solves q0 according to given parameters, equation 7.12.

Parameters
answerResult

Definition at line 625 of file dustmat.C.

References alpha, beta, lambda, newtonIter, newtonTol, OOFEM_ERROR, rEll, theta, and x0.

Referenced by initializeFrom().

double oofem::DustMaterial::yieldFunction1 ( double  rho,
double  i1 
)
protected

Yield function 1 (shear dominant), equation 7.5.

Parameters
rhoSecond Haigh-Westergaard coordinate
i1Trace of stress tensor
Returns
Value

Definition at line 595 of file dustmat.C.

References functionFe().

Referenced by performStressReturn().

double oofem::DustMaterial::yieldFunction2 ( double  rho,
double  i1,
double  q 
)
protected

Yield function 2 (compression dominant), equation 7.6.

Parameters
rhoSecond Haigh-Westergaard coordinate
i1Trace of stress tensor
qParameter q
Returns
value

Definition at line 601 of file dustmat.C.

References functionFc(), and functionFe().

Referenced by performStressReturn().

double oofem::DustMaterial::yieldFunction3 ( double  i1)
protected

Yield function 3 (tension dominant), equation 7.7.

Parameters
i1Trace of stress tensor
Returns
value

Definition at line 607 of file dustmat.C.

References ft.

Referenced by performStressReturn().

Member Data Documentation

double oofem::DustMaterial::alpha
protected

Parameter determining shape of yield surface.

Definition at line 240 of file dustmat.h.

Referenced by functionFe(), functionX(), initializeFrom(), and solveQ0().

double oofem::DustMaterial::beta
protected

Parameter determining shape of yield surface.

Definition at line 242 of file dustmat.h.

Referenced by computePlastStrainDirM1(), functionFe(), functionFeDI1(), functionFeDI1DI1(), functionX(), functionXDQ(), initializeFrom(), and solveQ0().

double oofem::DustMaterial::dHard
protected

Parameter determining hardening law (parameter D in original publication)

Definition at line 262 of file dustmat.h.

Referenced by functionH(), functionHDQ(), and initializeFrom().

double oofem::DustMaterial::ft
protected

Parameter determining shape of yield surface (param T in original publication)

Definition at line 248 of file dustmat.h.

Referenced by computePlastStrainDirM1(), computePlastStrainDirM3(), initializeFrom(), performF1return(), performStressReturn(), and yieldFunction3().

int oofem::DustMaterial::hardeningType
protected

Parameter determining hardening type.

Definition at line 252 of file dustmat.h.

Referenced by functionH(), functionHDQ(), and initializeFrom().

double oofem::DustMaterial::lambda
protected

Parameter determining shape of yield surface.

Definition at line 244 of file dustmat.h.

Referenced by computePlastStrainDirM1(), functionFe(), functionFeDI1(), functionFeDI1DI1(), functionX(), functionXDQ(), initializeFrom(), performStressReturn(), and solveQ0().

IsotropicLinearElasticMaterial* oofem::DustMaterial::LEMaterial
protected

Pointer for linear elastic material.

Definition at line 237 of file dustmat.h.

Referenced by computeAndSetBulkAndShearModuli(), DustMaterial(), give3dMaterialStiffnessMatrix(), initializeFrom(), and ~DustMaterial().

double oofem::DustMaterial::mStiff
protected

Parameter increasing stiffness (parameter M in original publication)

Definition at line 254 of file dustmat.h.

Referenced by computeAndSetBulkAndShearModuli(), and initializeFrom().

int oofem::DustMaterial::newtonIter
protected

Maximum number of iterations for iterative methods.

Definition at line 266 of file dustmat.h.

Referenced by computeQFromPlastVolEps(), initializeFrom(), performF1return(), performF2return(), and solveQ0().

double oofem::DustMaterial::newtonTol
protected

Tollerance for iterative methods.

Definition at line 264 of file dustmat.h.

Referenced by computeQFromPlastVolEps(), initializeFrom(), performF1return(), performF2return(), and solveQ0().

double oofem::DustMaterial::q0
protected

Parameter determining shape of yield surface.

Definition at line 258 of file dustmat.h.

Referenced by initializeFrom(), and performF2return().

double oofem::DustMaterial::rEll
protected

Parameter determining shape of yield surface (param R in original publication)

Definition at line 250 of file dustmat.h.

Referenced by computeDeltaGamma2(), computeDeltaGamma2DQ(), computePlastStrainDirM2(), fTempR2(), functionFc(), functionX(), functionXDQ(), initializeFrom(), and solveQ0().

double oofem::DustMaterial::theta
protected

Parameter determining shape of yield surface.

Definition at line 246 of file dustmat.h.

Referenced by computePlastStrainDirM1(), functionFe(), functionFeDI1(), functionX(), functionXDQ(), initializeFrom(), and solveQ0().

double oofem::DustMaterial::wHard
protected

Parameter determining hardening law (parameter W in original publication)

Definition at line 260 of file dustmat.h.

Referenced by functionH(), functionHDQ(), initializeFrom(), and performStressReturn().

double oofem::DustMaterial::x0
protected

Parameter determining shape of yield surface (param X0 in original publication)

Definition at line 256 of file dustmat.h.

Referenced by initializeFrom(), and solveQ0().


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:34 for OOFEM by doxygen 1.8.11 written by Dimitri van Heesch, © 1997-2011