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 MaterialStatus *	CreateStatus (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 MaterialStatus *	giveStatus (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...

Domain *	giveDomain () const

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

int	giveNumber () const

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

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

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

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

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

virtual Interface *	giveInterface (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
IsotropicLinearElasticMaterial *	LEMaterial
	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...

Domain *	domain
	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

Computes and sets all elastic moduli, with possible stiffening.

Equation 7.4

Parameters

bulkModulus	Bulk modulus
shearModulus	Shear modulus
gp	Gauss point

Definition at line 645 of file dustmat.C.

References oofem::IsotropicLinearElasticMaterial::computeBulkModulusFromYoungAndPoisson(), oofem::IsotropicLinearElasticMaterial::computeShearModulusFromYoungAndPoisson(), oofem::IsotropicLinearElasticMaterial::givePoissonsRatio(), oofem::Material::giveStatus(), oofem::DustMaterialStatus::giveVolumetricPlasticStrain(), oofem::IsotropicLinearElasticMaterial::giveYoungsModulus(), LEMaterial, mStiff, oofem::DustMaterialStatus::setBulkModulus(), oofem::DustMaterialStatus::setShearModulus(), and oofem::DustMaterialStatus::setYoungsModulus().

Referenced by performStressReturn().

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

tempQ	Parameter tempQ
q	Parameter q from previous step
i1	Trace of stress tensor
bulkModulus	Bulk 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

tempQ	Parameter tempQ
q	Parameter q from previous step
i1	Trace of stress tensor
bulkModulus	Bulk 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

answer	Result
stressDeviator	Deviator of stress tensor
rho	Second Haigh-Westergaard coordinate
i1	Trace of stress tensor
q	Parameter 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

answer	Result
stressDeviator	Deviator of stress tensor
rho	Second Haigh-Westergaard coordinate
i1	Trace of stress tensor
q	Parameter 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

answer	Result
stressDeviator	Deviator of stress tensor
rho	Second Haigh-Westergaard coordinate
i1	Trace of stress tensor
q	Parameter 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

answer	Result tempQ
q	Parameter q from previous step
deltaVolumetricPlasticStrain	Volumetric 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

gp	Integration 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

tempQ	Parameter tempQ
q	Parameter q from previous step
i1	Trace of stress tensor
rho	Second Haigh-Westergaard coordinate
bulkModulus	Bulk modulus
shearModulus	Shear 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

i1	Trace of stress tensor
rho	Second Haigh-Westergaard coordinate
q	Parameter 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

i1	Trace 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

i1	Trace 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

i1	Trace 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

q	Parameter q from previous step
tempQ	Parameter 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

tempQ Parameter 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

q	Parameter q from previous step
tempQ	Parameter tempQ
i1	Trace of stress tensor
bulkModulus	Bulk 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

tempQ	Parameter tempQ
bulkModulus	Bulk 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

q	Parameter 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

q	Parameter 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

answer	Computed results.
mode	Material response mode.
gp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Definition at line 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

answer	Contain corresponding ip value, zero sized if not available.
gp	Integration point to which the value refers.
type	Determines the type of internal variable.
tStep	Determines the time step.

Returns: Nonzero if the assignment can be done, zero if this type of variable is not supported.

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

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

Reimplemented from oofem::StructuralMaterial.

Definition at line 256 of file dustmat.C.

References oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStressDependentPartOfStrainVector_3d(), oofem::StructuralMaterialStatus::giveTempStressVector(), oofem::Material::initTempStatus(), oofem::StructuralMaterialStatus::letTempStrainVectorBe(), and performStressReturn().

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

answer	Vector of thermal dilatation coefficients.
gp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Definition at line 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

ir	Input record to initialize from.

Returns: IRResultType

Reimplemented from oofem::StructuralMaterial.

Definition at line 169 of file dustmat.C.

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

i1	Trace of stress tensor
rho	Second Haigh-Westergaard coordinate
gp	Gauss 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

i1	Trace of stress tensor
rho	Second Haigh-Westergaard coordinate
gp	Gauss 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().

void oofem::DustMaterial::performStressReturn	(	GaussPoint *	gp,
		const FloatArray &	strain
	)

protected

Perform stress return and update all internal variables.

Parameters

gp	Gauss point
strain	Strain

Definition at line 282 of file dustmat.C.

Referenced by giveRealStressVector_3d().

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

value	Contains the value(s) to be set (in reduced form).
gp	Integration point.
type	Determines the type of internal variable.
type	Determines the type of internal variable.