#include <dustmat.h>

Inheritance diagram for oofem::DustMaterial:

Collaboration diagram for oofem::DustMaterial:

Public Member Functions
	DustMaterial (int n, Domain *d)
	Constructor.
void	initializeFrom (InputRecord &ir) override
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
FloatArrayF< 6 >	giveRealStressVector_3d (const FloatArrayF< 6 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
int	setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type) override
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode) const override
FloatArrayF< 6 >	giveThermalDilatationVector (GaussPoint gp, TimeStep tStep) const override
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
double	giveQ0 () const
Public Member Functions inherited from oofem::StructuralMaterial
	StructuralMaterial (int n, Domain *d)
bool	hasMaterialModeCapability (MaterialMode mode) const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
void	giveCharacteristicMatrix (FloatMatrix &answer, MatResponseMode type, GaussPoint gp, TimeStep tStep) const override
	Returns characteristic matrix of the receiver.
void	giveCharacteristicVector (FloatArray &answer, FloatArray &flux, MatResponseMode type, GaussPoint gp, TimeStep tStep) const override
	Returns characteristic vector of the receiver.
virtual void	giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep) const
virtual void	giveRealStressVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep) const
virtual FloatArrayF< 4 >	giveRealStressVector_PlaneStrain (const FloatArrayF< 4 > &strain, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_3d.
virtual FloatArray	giveRealStressVector_StressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint gp, TimeStep tStep) const
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArray	giveRealStressVector_ShellStressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint gp, TimeStep tStep) const
virtual FloatArrayF< 3 >	giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 1 >	giveRealStressVector_1d (const FloatArrayF< 1 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 2 >	giveRealStressVector_Warping (const FloatArrayF< 2 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 2 >	giveRealStressVector_2dBeamLayer (const FloatArrayF< 2 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 5 >	giveRealStressVector_PlateLayer (const FloatArrayF< 5 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_Fiber (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_2dPlateSubSoil (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation is not provided.
virtual FloatArrayF< 6 >	giveRealStressVector_3dBeamSubSoil (const FloatArrayF< 6 > &reducedE, GaussPoint gp, TimeStep tStep) const
virtual FloatArrayF< 9 >	giveFirstPKStressVector_3d (const FloatArrayF< 9 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
virtual FloatArrayF< 5 >	giveFirstPKStressVector_PlaneStrain (const FloatArrayF< 5 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_3d.
virtual FloatArray	giveFirstPKStressVector_StressControl (const FloatArray &reducedvF, const IntArray &FControl, GaussPoint gp, TimeStep tStep) const
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArrayF< 4 >	giveFirstPKStressVector_PlaneStress (const FloatArrayF< 4 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
virtual FloatArrayF< 1 >	giveFirstPKStressVector_1d (const FloatArrayF< 1 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
virtual void	giveCauchyStressVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveCauchyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveCauchyStressVector_PlaneStress (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveCauchyStressVector_1d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
double	giveReferenceTemperature ()
virtual FloatArray	computeStressIndependentStrainVector (GaussPoint gp, TimeStep tStep, ValueModeType mode) const
FloatArrayF< 6 >	computeStressIndependentStrainVector_3d (GaussPoint gp, TimeStep tStep, ValueModeType mode) const
virtual void	giveStiffnessMatrix_dPdF (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 9, 9 >	give3dMaterialStiffnessMatrix_dPdF (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const
virtual void	give3dMaterialStiffnessMatrix_dCde (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
void	giveStressDependentPartOfStrainVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrainVector, TimeStep tStep, ValueModeType mode) const
int	setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type) override
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
virtual FloatMatrixF< 3, 3 >	givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 4, 4 >	givePlaneStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	givePlaneStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 4, 4 >	givePlaneStrainStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 5, 5 >	givePlaneStrainStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	givePlaneStrainStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 1, 1 >	give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 1, 1 >	give1dStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	give1dStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 2, 2 >	give2dBeamLayerStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 5, 5 >	givePlateLayerStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 3, 3 >	giveFiberStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 3, 3 >	give2dPlateSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 6, 6 >	give3dBeamSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
Public Member Functions inherited from oofem::Material
	Material (int n, Domain *d)
virtual	~Material ()=default
	Destructor.
virtual double	giveCharacteristicValue (MatResponseMode type, GaussPoint gp, TimeStep tStep) const
	Returns characteristic value of the receiver.
virtual double	give (int aProperty, GaussPoint *gp) const
virtual bool	hasProperty (int aProperty, GaussPoint *gp) const
virtual void	modifyProperty (int aProperty, double value, GaussPoint *gp)
double	giveCastingTime () const
virtual bool	isActivated (TimeStep *tStep) const
virtual bool	hasCastingTimeSupport () const
void	printYourself () override
	Prints receiver state on stdout. Useful for debugging.
virtual void	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
virtual void	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
int	checkConsistency () override
virtual void	restoreConsistency (GaussPoint *gp)
virtual int	initMaterial (Element *element)
virtual MaterialStatus *	giveStatus (GaussPoint *gp) const
virtual int	packUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
virtual int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
virtual int	estimatePackSize (DataStream &buff, GaussPoint *ip)
virtual double	predictRelativeComputationalCost (GaussPoint *gp)
virtual double	predictRelativeRedistributionCost (GaussPoint *gp)
virtual void	initTempStatus (GaussPoint *gp) const
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
virtual	~FEMComponent ()=default
	Virtual destructor.
Domain *	giveDomain () const
virtual void	setDomain (Domain *d)
int	giveNumber () const
void	setNumber (int num)
virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
virtual void	initializeFrom (InputRecord &ir, int priority)
virtual void	initializeFinish ()
virtual void	postInitialize ()
	Performs post initialization steps. Called after all components are created and initialized.
virtual void	printOutputAt (FILE file, TimeStep tStep)
virtual Interface *	giveInterface (InterfaceType t)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).

Protected Member Functions
double	functionFe (double i1) const
double	functionFeDI1 (double i1) const
double	functionFeDI1DI1 (double i1) const
double	functionFc (double rho, double i1, double q) const
double	functionX (double q) const
double	functionXDQ (double q) const
double	yieldFunction1 (double rho, double i1) const
double	yieldFunction2 (double rho, double i1, double q) const
double	yieldFunction3 (double i1) const
void	solveQ0 (double &answer) const
void	computeAndSetBulkAndShearModuli (double &bulkModulus, double &shearModulus, GaussPoint *gp) const
void	performStressReturn (GaussPoint *gp, const FloatArrayF< 6 > &strain) const
FloatArrayF< 6 >	computePlastStrainDirM1 (const FloatArrayF< 6 > &stressDeviator, double rho, double i1, double q) const
FloatArrayF< 6 >	computePlastStrainDirM2 (const FloatArrayF< 6 > &stressDeviator, double rho, double i1, double q) const
FloatArrayF< 6 >	computePlastStrainDirM3 (const FloatArrayF< 6 > &stressDeviator, double rho, double i1, double q) const
double	functionH (double q, double tempQ) const
double	functionHDQ (double tempQ) const
double	functionI1 (double q, double tempQ, double i1, double bulkModulus) const
double	functionI1DQ (double tempQ, double bulkModulus) const
void	performF1return (double i1, double rho, GaussPoint *gp) const
void	performF2return (double i1, double rho, GaussPoint *gp) const
void	computeQFromPlastVolEps (double &answer, double q, double deltaVolumetricPlasticStrain) const
double	computeDeltaGamma2 (double tempQ, double q, double i1, double bulkModulus) const
double	computeDeltaGamma2DQ (double tempQ, double q, double i1, double bulkModulus) const
double	fTempR2 (double tempQ, double q, double i1, double rho, double bulkModulus, double shearModulus) const

Protected Attributes
IsotropicLinearElasticMaterial	LEMaterial
	Pointer for linear elastic material.
double	alpha = 0.
	Parameter determining shape of yield surface.
double	beta = 0.
	Parameter determining shape of yield surface.
double	lambda = 0.
	Parameter determining shape of yield surface.
double	theta = 0.
	Parameter determining shape of yield surface.
double	ft = 0.
	Parameter determining shape of yield surface (param T in original publication).
double	rEll = 0.
	Parameter determining shape of yield surface (param R in original publication).
int	hardeningType = 0
	Parameter determining hardening type.
double	mStiff = 0.
	Parameter increasing stiffness (parameter M in original publication).
double	x0 = 0.
	Parameter determining shape of yield surface (param X0 in original publication).
double	q0 = 0.
	Parameter determining shape of yield surface.
double	wHard = 0.
	Parameter determining hardening law (parameter W in original publication).
double	dHard = 0.
	Parameter determining hardening law (parameter D in original publication).
double	newtonTol = 0.
	Tollerance for iterative methods.
int	newtonIter = 0
	Maximum number of iterations for iterative methods.
Protected Attributes inherited from oofem::StructuralMaterial
double	referenceTemperature = 0.
	Reference temperature (temperature, when material has been built into structure).
MatResponseMode	SCStiffMode = TangentStiffness
	stifness mode used in stress control
double	SCRelTol = 1.e-3
	relative tolerance for stress control
double	SCAbsTol = 1.e-12
	absolute stress tolerance for stress control
int	SCMaxiter = 100000
	maximum iterations for stress-control
Protected Attributes inherited from oofem::Material
Dictionary	propertyDictionary
double	castingTime
int	preCastingTimeMat
	Material existing before casting time - optional parameter, zero by default.
Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number.
Domain *	domain
	Link to domain object, useful for communicating with other FEM components.

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

Detailed Description

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

Constructor & Destructor Documentation

◆ DustMaterial()

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

Constructor.

Definition at line 140 of file dustmat.C.

References LEMaterial, and oofem::StructuralMaterial::StructuralMaterial().

Member Function Documentation

◆ computeAndSetBulkAndShearModuli()

void oofem::DustMaterial::computeAndSetBulkAndShearModuli	(	double &	bulkModulus,
		double &	shearModulus,
		GaussPoint *	gp ) const

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 601 of file dustmat.C.

References oofem::Material::giveStatus(), LEMaterial, and mStiff.

Referenced by performStressReturn().

◆ computeDeltaGamma2()

double oofem::DustMaterial::computeDeltaGamma2	(	double	tempQ,
		double	q,
		double	i1,
		double	bulkModulus ) const

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 686 of file dustmat.C.

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

Referenced by fTempR2(), and performStressReturn().

◆ computeDeltaGamma2DQ()

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

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 693 of file dustmat.C.

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

◆ computePlastStrainDirM1()

FloatArrayF< 6 > oofem::DustMaterial::computePlastStrainDirM1	(	const FloatArrayF< 6 > &	stressDeviator,
		double	rho,
		double	i1,
		double	q ) const

protected

Computes direction of plastic yielding m1, equation 7.17

Parameters

stressDeviator	Deviator of stress tensor
rho	Second Haigh-Westergaard coordinate
i1	Trace of stress tensor
q	Parameter q
return	Result

Definition at line 619 of file dustmat.C.

References beta, oofem::StructuralMaterial::computeDeviatoricVolumetricSum(), ft, lambda, and theta.

Referenced by performStressReturn().

◆ computePlastStrainDirM2()

FloatArrayF< 6 > oofem::DustMaterial::computePlastStrainDirM2	(	const FloatArrayF< 6 > &	stressDeviator,
		double	rho,
		double	i1,
		double	q ) const

protected

Computes direction of plastic yielding m2, equation 7.19

Parameters

stressDeviator	Deviator of stress tensor
rho	Second Haigh-Westergaard coordinate
i1	Trace of stress tensor
q	Parameter q
return	Result

Definition at line 627 of file dustmat.C.

References oofem::StructuralMaterial::computeDeviatoricVolumetricSum(), functionFc(), and rEll.

Referenced by performStressReturn().

◆ computePlastStrainDirM3()

FloatArrayF< 6 > oofem::DustMaterial::computePlastStrainDirM3	(	const FloatArrayF< 6 > &	stressDeviator,
		double	rho,
		double	i1,
		double	q ) const

protected

Computes direction of plastic yielding m2, equation 7.18

Parameters

stressDeviator	Deviator of stress tensor
rho	Second Haigh-Westergaard coordinate
i1	Trace of stress tensor
q	Parameter q
return	Result

Definition at line 636 of file dustmat.C.

References oofem::StructuralMaterial::computeDeviatoricVolumetricSum(), ft, functionFe(), and functionFeDI1().

Referenced by performStressReturn().

◆ computeQFromPlastVolEps()

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

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 436 of file dustmat.C.

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

Referenced by performStressReturn().

◆ CreateStatus()

std::unique_ptr< MaterialStatus > oofem::DustMaterial::CreateStatus ( GaussPoint * gp ) const

overridevirtual

Creates new copy of associated status and inserts it into given integration point.

Parameters

gp	Integration point where newly created status will be stored.

Returns: Reference to new status.

Reimplemented from oofem::Material.

Definition at line 522 of file dustmat.C.

References giveQ0().

◆ fTempR2()

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

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 707 of file dustmat.C.

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

Referenced by performF2return().

◆ functionFc()

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

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 546 of file dustmat.C.

References rEll.

Referenced by computePlastStrainDirM2(), and yieldFunction2().

◆ functionFe()

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

protected

Auxiliary equation Fe (7.8)

Parameters

i1	Trace of stress tensor

Returns: Fe

Definition at line 528 of file dustmat.C.

References alpha, beta, lambda, and theta.

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

◆ functionFeDI1()

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

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 534 of file dustmat.C.

References beta, lambda, and theta.

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

◆ functionFeDI1DI1()

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

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 540 of file dustmat.C.

References beta, and lambda.

Referenced by performF1return().

◆ functionH()

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

protected

Auxiliary equation H (7.33 or 7.34)

Parameters

q	Parameter q from previous step
tempQ	Parameter tempQ

Returns: H

Definition at line 646 of file dustmat.C.

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

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

◆ functionHDQ()

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

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 660 of file dustmat.C.

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

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

◆ functionI1()

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

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 674 of file dustmat.C.

References functionH().

Referenced by performF1return(), and performStressReturn().

◆ functionI1DQ()

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

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 680 of file dustmat.C.

References functionHDQ().

Referenced by performF1return().

◆ functionX()

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

protected

Auxiliary equation X (7.11)

Parameters

q	Parameter q

Returns: X

Definition at line 570 of file dustmat.C.

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

Referenced by functionH(), and functionHDQ().

◆ functionXDQ()

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

protected

Derivative by q of auxiliary equation X (7.11)

Parameters

q	Parameter q

Returns: \( \frac{\partial X }{\partial q} \)

Definition at line 576 of file dustmat.C.

References beta, lambda, rEll, and theta.

Referenced by functionHDQ().

◆ give3dMaterialStiffnessMatrix()

FloatMatrixF< 6, 6 > oofem::DustMaterial::give3dMaterialStiffnessMatrix	(	MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point.

Parameters

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

Reimplemented from oofem::StructuralMaterial.

Definition at line 461 of file dustmat.C.

References oofem::Material::giveStatus(), LEMaterial, and OOFEM_ERROR.

◆ giveClassName()

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

inlineoverridevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 452 of file dustmat.h.

◆ giveInputRecordName()

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

inlineoverridevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 453 of file dustmat.h.

References _IFT_DustMaterial_Name.

◆ giveIPValue()

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

overridevirtual

Returns the integration point corresponding value in Reduced form.

Parameters

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

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

Todo: This is actually the mean, not the volumetric part. / Mikael

Reimplemented from oofem::Material.

Definition at line 494 of file dustmat.C.

References oofem::FloatArray::at(), oofem::StructuralMaterial::computePrincipalValues(), oofem::Material::giveStatus(), oofem::principal_strain, and oofem::FloatArray::resize().

◆ giveQ0()

double oofem::DustMaterial::giveQ0 ( ) const

inline

Definition at line 472 of file dustmat.h.

References q0.

Referenced by CreateStatus().

◆ giveRealStressVector_3d()

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

overridevirtual

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

Reimplemented from oofem::StructuralMaterial.

Definition at line 220 of file dustmat.C.

References oofem::StructuralMaterial::computeStressIndependentStrainVector_3d(), oofem::Material::giveStatus(), oofem::Material::initTempStatus(), and performStressReturn().

◆ giveThermalDilatationVector()

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

inlineoverridevirtual

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

References LEMaterial.

◆ initializeFrom()

void oofem::DustMaterial::initializeFrom ( InputRecord & ir )

overridevirtual

Initializes receiver according to object description stored in input record. This function is called immediately after creating object using constructor. Input record can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record. Note that initializeFrom may be called mutiple times.

See also: IR_GIVE_FIELD; IR_GIVE_OPTIONAL_FIELD

Parameters

ir	Input record to initialize from.
priority	Priority of the input record. This is used to determine the order of initialization

Reimplemented from oofem::FEMComponent.

Definition at line 147 of file dustmat.C.

◆ isCharacteristicMtrxSymmetric()

bool oofem::DustMaterial::isCharacteristicMtrxSymmetric ( MatResponseMode rMode ) const

inlineoverridevirtual

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

Reimplemented from oofem::Material.

Definition at line 463 of file dustmat.h.

◆ performF1return()

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

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 354 of file dustmat.C.

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

Referenced by performStressReturn().

◆ performF2return()

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

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 398 of file dustmat.C.

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

Referenced by performStressReturn().

◆ performStressReturn()

void oofem::DustMaterial::performStressReturn	(	GaussPoint *	gp,
		const FloatArrayF< 6 > &	strain ) const

protected

Perform stress return and update all internal variables

Parameters

gp	Gauss point
strain	Strain

Definition at line 243 of file dustmat.C.

References computeAndSetBulkAndShearModuli(), computeDeltaGamma2(), oofem::StructuralMaterial::computeDeviatoricVolumetricSplit(), oofem::StructuralMaterial::computeDeviatoricVolumetricSum(), computePlastStrainDirM1(), computePlastStrainDirM2(), computePlastStrainDirM3(), computeQFromPlastVolEps(), oofem::StructuralMaterial::computeSecondCoordinate(), ft, functionFe(), functionFeDI1(), functionI1(), oofem::Material::giveStatus(), lambda, OOFEM_ERROR, performF1return(), performF2return(), wHard, yieldFunction1(), yieldFunction2(), and yieldFunction3().

Referenced by giveRealStressVector_3d().

◆ setIPValue()

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

overridevirtual

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.

Returns: Nonzero if ok, zero if var not supported.

Reimplemented from oofem::Material.

Definition at line 479 of file dustmat.C.

References oofem::FloatArray::at(), and oofem::Material::giveStatus().

◆ solveQ0()

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

protected

Solves q0 according to given parameters, equation 7.12

Parameters

answer Result

Definition at line 582 of file dustmat.C.

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

Referenced by initializeFrom().

◆ yieldFunction1()

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

protected

Yield function 1 (shear dominant), equation 7.5

Parameters

rho	Second Haigh-Westergaard coordinate
i1	Trace of stress tensor

Returns: Value

Definition at line 552 of file dustmat.C.

References functionFe().

Referenced by performStressReturn().

◆ yieldFunction2()

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

protected

Yield function 2 (compression dominant), equation 7.6

Parameters

rho	Second Haigh-Westergaard coordinate
i1	Trace of stress tensor
q	Parameter q

Returns: value

Definition at line 558 of file dustmat.C.

References functionFc(), and functionFe().

Referenced by performStressReturn().

◆ yieldFunction3()

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

protected

Yield function 3 (tension dominant), equation 7.7

Parameters

i1	Trace of stress tensor

Returns: value

Definition at line 564 of file dustmat.C.

References ft.

Referenced by performStressReturn().

Member Data Documentation

◆ alpha

double oofem::DustMaterial::alpha = 0.

protected

Parameter determining shape of yield surface.

Definition at line 237 of file dustmat.h.

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

◆ beta

double oofem::DustMaterial::beta = 0.

protected

Parameter determining shape of yield surface.

Definition at line 239 of file dustmat.h.

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

◆ dHard

double oofem::DustMaterial::dHard = 0.

protected

Parameter determining hardening law (parameter D in original publication).

Definition at line 259 of file dustmat.h.

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

◆ ft

double oofem::DustMaterial::ft = 0.

protected

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

Definition at line 245 of file dustmat.h.

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

◆ hardeningType

int oofem::DustMaterial::hardeningType = 0

protected

Parameter determining hardening type.

Definition at line 249 of file dustmat.h.

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

◆ lambda

double oofem::DustMaterial::lambda = 0.

protected

Parameter determining shape of yield surface.

Definition at line 241 of file dustmat.h.

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

◆ LEMaterial

IsotropicLinearElasticMaterial oofem::DustMaterial::LEMaterial

protected

Pointer for linear elastic material.

Definition at line 234 of file dustmat.h.

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

◆ mStiff

double oofem::DustMaterial::mStiff = 0.

protected

Parameter increasing stiffness (parameter M in original publication).

Definition at line 251 of file dustmat.h.

Referenced by computeAndSetBulkAndShearModuli(), and initializeFrom().

◆ newtonIter

int oofem::DustMaterial::newtonIter = 0

protected

Maximum number of iterations for iterative methods.

Definition at line 263 of file dustmat.h.

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

◆ newtonTol

double oofem::DustMaterial::newtonTol = 0.

protected

Tollerance for iterative methods.

Definition at line 261 of file dustmat.h.

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

◆ q0

double oofem::DustMaterial::q0 = 0.

protected

Parameter determining shape of yield surface.

Definition at line 255 of file dustmat.h.

Referenced by giveQ0(), initializeFrom(), and performF2return().

◆ rEll

double oofem::DustMaterial::rEll = 0.

protected

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

Definition at line 247 of file dustmat.h.

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

◆ theta

double oofem::DustMaterial::theta = 0.

protected

Parameter determining shape of yield surface.

Definition at line 243 of file dustmat.h.

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

◆ wHard

double oofem::DustMaterial::wHard = 0.

protected

Parameter determining hardening law (parameter W in original publication).

Definition at line 257 of file dustmat.h.

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

◆ x0

double oofem::DustMaterial::x0 = 0.

protected

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

Definition at line 253 of file dustmat.h.

Referenced by initializeFrom(), and solveQ0().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DustMaterial()

Member Function Documentation

◆ computeAndSetBulkAndShearModuli()

◆ computeDeltaGamma2()

◆ computeDeltaGamma2DQ()

◆ computePlastStrainDirM1()

◆ computePlastStrainDirM2()

◆ computePlastStrainDirM3()

◆ computeQFromPlastVolEps()

◆ CreateStatus()

◆ fTempR2()

◆ functionFc()

◆ functionFe()

◆ functionFeDI1()

◆ functionFeDI1DI1()

◆ functionH()

◆ functionHDQ()

◆ functionI1()

◆ functionI1DQ()

◆ functionX()

◆ functionXDQ()

◆ give3dMaterialStiffnessMatrix()

◆ giveClassName()

◆ giveInputRecordName()

◆ giveIPValue()

◆ giveQ0()

◆ giveRealStressVector_3d()

◆ giveThermalDilatationVector()

◆ initializeFrom()

◆ isCharacteristicMtrxSymmetric()

◆ performF1return()

◆ performF2return()

◆ performStressReturn()

◆ setIPValue()

◆ solveQ0()

◆ yieldFunction1()

◆ yieldFunction2()

◆ yieldFunction3()

Member Data Documentation

◆ alpha

◆ beta

◆ dHard

◆ ft

◆ hardeningType

◆ lambda

◆ LEMaterial

◆ mStiff

◆ newtonIter

◆ newtonTol

◆ q0

◆ rEll

◆ theta

◆ wHard

◆ x0