#include <misesmatnl.h>

Inheritance diagram for oofem::MisesMatNl:

Collaboration diagram for oofem::MisesMatNl:

Public Member Functions
	MisesMatNl (int n, Domain *d)
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
Interface *	giveInterface (InterfaceType) override
virtual double	computeCumPlasticStrain (GaussPoint gp, TimeStep tStep) const
double	computeDamage (GaussPoint gp, TimeStep tStep) const
void	modifyNonlocalWeightFunctionAround (GaussPoint *gp) const
double	computeDistanceModifier (double damage) const
double	computeLocalCumPlasticStrain (GaussPoint gp, TimeStep tStep) const
FloatMatrixF< 1, 1 >	give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
void	NonlocalMaterialStiffnessInterface_addIPContribution (SparseMtrx &dest, const UnknownNumberingScheme &s, GaussPoint gp, TimeStep tStep) override
	Computes and adds IP contributions to destination matrix.
std ::vector< localIntegrationRecord > *	NonlocalMaterialStiffnessInterface_giveIntegrationDomainList (GaussPoint *gp) override
int	giveLocalNonlocalStiffnessContribution (GaussPoint gp, IntArray &loc, const UnknownNumberingScheme &s, FloatArray &lcontrib, TimeStep tStep)
void	giveRemoteNonlocalStiffnessContribution (GaussPoint gp, IntArray &rloc, const UnknownNumberingScheme &s, FloatArray &rcontrib, TimeStep tStep)
FloatArrayF< 6 >	giveRealStressVector_3d (const FloatArrayF< 6 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
FloatArrayF< 1 >	giveRealStressVector_1d (const FloatArrayF< 1 > &strainVector, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
void	updateBeforeNonlocAverage (const FloatArray &strainVector, GaussPoint gp, TimeStep tStep) const override
int	hasBoundedSupport () const override
int	packUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip) override
int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip) override
int	estimatePackSize (DataStream &buff, GaussPoint *ip) override
Public Member Functions inherited from oofem::MisesMat
	MisesMat (int n, Domain *d)
void	performPlasticityReturn (const FloatArray &totalStrain, GaussPoint gp, TimeStep tStep) const
void	performPlasticityReturn_PlaneStress (const FloatArrayF< 3 > &totalStrain, GaussPoint gp, TimeStep tStep) const
double	checkYieldStress (double &dKappa, double kappa, GaussPoint gp, TimeStep tStep) const
double	computeYieldStress (double kappa, GaussPoint gp, TimeStep tStep) const
double	computeYieldStressPrime (double kappa) const
double	computeDamage (GaussPoint gp, TimeStep tStep) const
double	computeDamageParam (double tempKappa) const
double	computeDamageParamPrime (double tempKappa) const
virtual double	computeCumPlastStrain (GaussPoint gp, TimeStep tStep) const
bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode) const override
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 3, 3 >	givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
FloatArrayF< 3 >	giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &totalStrain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector_StressControl.
double	giveS (GaussPoint gp, TimeStep tStep) const
double	giveTemperature (GaussPoint gp, TimeStep tStep) 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< 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)
virtual FloatArrayF< 6 >	giveThermalDilatationVector (GaussPoint gp, TimeStep tStep) const
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< 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 >	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 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)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).
Public Member Functions inherited from oofem::StructuralNonlocalMaterialExtensionInterface
	StructuralNonlocalMaterialExtensionInterface (Domain *d)
	~StructuralNonlocalMaterialExtensionInterface ()
	Destructor.
Public Member Functions inherited from oofem::NonlocalMaterialExtensionInterface
	NonlocalMaterialExtensionInterface (Domain *d)
virtual	~NonlocalMaterialExtensionInterface ()
	Destructor.
void	updateDomainBeforeNonlocAverage (TimeStep *tStep) const
void	buildNonlocalPointTable (GaussPoint *gp) const
void	rebuildNonlocalPointTable (GaussPoint gp, IntArray contributingElems) const
void	modifyNonlocalWeightFunctionAround (GaussPoint *gp) const
void	modifyNonlocalWeightFunction_1D_Around (GaussPoint *gp) const
virtual double	giveNonlocalMetricModifierAt (GaussPoint *gp) const
double	computeDistanceModifier (double cl, double damage) const
double	computeModifiedLength (double length, double dam1, double dam2) const
std ::vector< localIntegrationRecord > *	giveIPIntegrationList (GaussPoint *gp) const
virtual double	computeWeightFunction (const double cl, const double distance) const
virtual double	computeWeightFunction (const double cl, const FloatArray &src, const FloatArray &coord) const
double	giveIntegralOfWeightFunction (double cl, const int spatial_dimension) const
virtual double	maxValueOfWeightFunction ()
	Determines the maximum value of the nonlocal weight function.
int	giveNumberOfRegions () const
virtual double	evaluateSupportRadius (double cl) const
Domain *	giveDomain ()
	Returns reference to domain.
void	initializeFrom (InputRecord &ir)
void	giveInputRecord (DynamicInputRecord &input)
void	endIPNonlocalAverage (GaussPoint *gp) const
Public Member Functions inherited from oofem::Interface
	Interface ()
	Constructor.
virtual	~Interface ()
Public Member Functions inherited from oofem::NonlocalMaterialStiffnessInterface
	NonlocalMaterialStiffnessInterface ()
	Constructor.
virtual void	NonlocalMaterialStiffnessInterface_showSparseMtrxStructure (GaussPoint gp, oofegGraphicContext &gc, TimeStep )

Protected Member Functions
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
Protected Member Functions inherited from oofem::MisesMat
void	computeGLPlasticStrain (const FloatMatrix &F, FloatMatrix &Ep, FloatMatrix b, double J)
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
Protected Member Functions inherited from oofem::NonlocalMaterialExtensionInterface
void	applyBarrierConstraints (const FloatArray &gpCoords, const FloatArray &jGpCoords, double &weight) const
void	manipulateWeight (double &weight, GaussPoint gp, GaussPoint jGp) const
double	giveDistanceBasedInteractionRadius (const FloatArray &gpCoords) const
int	mapToGridPoint (double x, double x0) const
double	mapToGridCoord (double x, double x0) const
double	dist2FromGridNode (double x, double y, int j, int i) const

Protected Attributes
double	Rf = 0.
double	exponent = 1.
int	averType = 0
Protected Attributes inherited from oofem::MisesMat
IsotropicLinearElasticMaterial	linearElasticMaterial
	Reference to the basic elastic material.
double	G = 0.
	Elastic shear modulus.
double	K = 0.
	Elastic bulk modulus.
double	H = 0.
	Hardening modulus.
ScalarFunction	sig0
	Initial (uniaxial) yield stress.
int	hType
	type of hardening function
FloatArray	h_eps
	user-defined hardening (yield stress - kappa)
FloatArray	h_function_eps
double	omega_crit = 0.
	critical(maximal) damage.
double	a = 0.
	exponent in damage function.
double	yieldTol = 0.
	tolerance for the yield function in RRM algorithm.
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.
Protected Attributes inherited from oofem::NonlocalMaterialExtensionInterface
Domain *	domain = nullptr
IntArray	regionMap
	Map indicating regions to skip (region - cross section model).
bool	permanentNonlocTableFlag = false
	Flag indicating whether to keep nonlocal interaction tables of integration points cached.
WeightFunctionType	weightFun
	Parameter specifying the type of nonlocal weight function.
int	gridSize = 0
	Grid on which the eikonal equation will be solved (used by eikonal nonlocal models).
double	initDiag = 0.
	Optional parameters setting details of the fast marching method.
int	order = 0
int	centDiff = 0
double	cl = 0.
double	suprad = 0.
	Support radius.
double	mm = 0.
	For "undernonlocal" or "overnonlocal" formulation.
ScalingType	scaling
	Parameter specifying the type of scaling of nonlocal weight function.
AveragedVarType	averagedVar
	Parameter specifying the type of averaged (nonlocal) variable.
double	cl0 = 0.
NlVariationType	nlvar
	Parameter specifying the type of nonlocal variation.
double	beta = 0.
double	zeta = 0.
double	px = 0.
double	Rf = 0.
	Final value of interaction radius, for a model with evolving characteristic length.
double	exponent = 0.
	Parameter used as an exponent by models with evolving characteristic length.
int	averType = 0
	Parameter specifying how the weight function should be adjusted due to damage.

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.
Protected Types inherited from oofem::NonlocalMaterialExtensionInterface
enum	ScalingType { ST_Unknown , ST_Standard , ST_Noscaling , ST_Borino }
	Type characterizing the scaling approach. More...
enum	AveragedVarType { AVT_Unknown , AVT_EqStrain , AVT_Compliance , AVT_Damage }
	Type characterizing the averaged (nonlocal) variable. More...
enum	NlVariationType { NLVT_Standard , NLVT_DistanceBasedLinear , NLVT_StressBased , NLVT_DistanceBasedExponential }
	Type characterizing the Nonlocal variation. More...

Detailed Description

Mises nonlocal material.

Author: Milan

Definition at line 88 of file misesmatnl.h.

Constructor & Destructor Documentation

◆ MisesMatNl()

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

Definition at line 51 of file misesmatnl.C.

References oofem::MisesMat::MisesMat(), oofem::NonlocalMaterialStiffnessInterface::NonlocalMaterialStiffnessInterface(), and oofem::StructuralNonlocalMaterialExtensionInterface::StructuralNonlocalMaterialExtensionInterface().

Referenced by NonlocalMaterialStiffnessInterface_addIPContribution().

Member Function Documentation

◆ computeCumPlasticStrain()

double oofem::MisesMatNl::computeCumPlasticStrain	(	GaussPoint *	gp,
		TimeStep *	tStep ) const

virtual

Computes the nonlocal cumulated plastic strain from its local form.

Parameters

[out]	kappa	Return param, containing the corresponding cumulated plastic strain.
	gp	Integration point.
	tStep	Time step.

Definition at line 233 of file misesmatnl.C.

References oofem::NonlocalMaterialExtensionInterface::buildNonlocalPointTable(), oofem::NonlocalMaterialExtensionInterface::giveIPIntegrationList(), oofem::MisesMatNlStatus::giveLocalCumPlasticStrainForAverage(), oofem::Material::giveStatus(), oofem::NonlocalMaterialExtensionInterface::mm, oofem::NonlocalMaterialExtensionInterface::scaling, oofem::NonlocalMaterialExtensionInterface::ST_Borino, oofem::NonlocalMaterialExtensionInterface::ST_Standard, and oofem::NonlocalMaterialExtensionInterface::updateDomainBeforeNonlocAverage().

Referenced by computeDamage(), give1dStressStiffMtrx(), and giveLocalNonlocalStiffnessContribution().

◆ computeDamage()

double oofem::MisesMatNl::computeDamage	(	GaussPoint *	gp,
		TimeStep *	tStep ) const

Definition at line 326 of file misesmatnl.C.

References computeCumPlasticStrain(), oofem::MisesMat::computeDamageParam(), and oofem::Material::giveStatus().

Referenced by giveRealStressVector_1d().

◆ computeDistanceModifier()

double oofem::MisesMatNl::computeDistanceModifier ( double damage ) const

Definition at line 213 of file misesmatnl.C.

References averType, oofem::NonlocalMaterialExtensionInterface::cl, exponent, M_PI, and Rf.

Referenced by modifyNonlocalWeightFunctionAround().

◆ computeLocalCumPlasticStrain()

double oofem::MisesMatNl::computeLocalCumPlasticStrain	(	GaussPoint *	gp,
		TimeStep *	tStep ) const

inline

Definition at line 117 of file misesmatnl.h.

Referenced by updateBeforeNonlocAverage().

◆ CreateStatus()

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

inlineoverrideprotectedvirtual

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

Definition at line 166 of file misesmatnl.h.

◆ estimatePackSize()

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

overridevirtual

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

Reimplemented from oofem::Material.

Definition at line 526 of file misesmatnl.C.

References oofem::DataStream::givePackSizeOfDouble().

◆ give1dStressStiffMtrx()

FloatMatrixF< 1, 1 > oofem::MisesMatNl::give1dStressStiffMtrx	(	MatResponseMode	mmode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Method for computing 1d stiffness matrix of receiver. Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to 1d stiffness using reduce method described above. However, this reduction is quite time consuming and if it is possible, it is recommended to overload this method and provide direct method for computing particular stiffness matrix.

Parameters

mmode	Material response mode.
gp	Integration point, which load history is used.
tStep	Time step (most models are able to respond only when tStep is current time step).

Returns: Stiffness matrix.

Reimplemented from oofem::MisesMat.

Definition at line 79 of file misesmatnl.C.

References computeCumPlasticStrain(), oofem::MisesMat::computeDamageParamPrime(), E, oofem::Material::giveStatus(), oofem::MisesMat::H, oofem::MisesMat::linearElasticMaterial, oofem::NonlocalMaterialExtensionInterface::mm, and oofem::sgn().

◆ giveClassName()

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

inlineoverridevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::MisesMat.

Definition at line 99 of file misesmatnl.h.

◆ giveInputRecord()

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

overridevirtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::Material.

Definition at line 308 of file misesmatnl.C.

References _IFT_MisesMatNl_averagingtype, _IFT_MisesMatNl_exp, _IFT_MisesMatNl_rf, averType, exponent, Rf, and oofem::DynamicInputRecord::setField().

◆ giveInputRecordName()

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

inlineoverridevirtual

Returns: Input record name of the receiver.

Reimplemented from oofem::MisesMat.

Definition at line 100 of file misesmatnl.h.

References _IFT_MisesMatNl_Name.

◆ giveInterface()

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

overridevirtual

Interface requesting service.

See also: InterfaceType

Returns: Requested interface if implemented, otherwise NULL.

Reimplemented from oofem::FEMComponent.

Definition at line 269 of file misesmatnl.C.

References oofem::NonlocalMaterialExtensionInterfaceType, oofem::NonlocalMaterialStiffnessInterface::NonlocalMaterialStiffnessInterface(), oofem::NonlocalMaterialStiffnessInterfaceType, and oofem::StructuralNonlocalMaterialExtensionInterface::StructuralNonlocalMaterialExtensionInterface().

◆ giveLocalNonlocalStiffnessContribution()

int oofem::MisesMatNl::giveLocalNonlocalStiffnessContribution	(	GaussPoint *	gp,
		IntArray &	loc,
		const UnknownNumberingScheme &	s,
		FloatArray &	lcontrib,
		TimeStep *	tStep )

Computes the "local" part of nonlocal stiffness contribution assembled for given integration point.

Parameters

gp	Source integration point
loc	Local code numbers
lcontrib	"Local" contribution
s	Numbering scheme that determines assembly.
tStep	Time step.

Returns: Nonzero if local point contributes (loading) or zero if not (unloading in elastic range, elastic)

Definition at line 381 of file misesmatnl.C.

References oofem::FloatArray::clear(), computeCumPlasticStrain(), oofem::MisesMat::computeDamageParamPrime(), oofem::GaussPoint::giveElement(), oofem::Material::giveStatus(), oofem::NonlocalMaterialExtensionInterface::mm, and oofem::FloatArray::plusProduct().

Referenced by NonlocalMaterialStiffnessInterface_addIPContribution().

◆ giveRealStressVector_1d()

FloatArrayF< 1 > oofem::MisesMatNl::giveRealStressVector_1d	(	const FloatArrayF< 1 > &	reducedE,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Default implementation relies on giveRealStressVector_StressControl.

Note: : One should obtain the same answer using the iterations in the default implementation (this is verified for this model).

Reimplemented from oofem::MisesMat.

Definition at line 64 of file misesmatnl.C.

References computeDamage(), oofem::Material::giveStatus(), and oofem::MisesMat::performPlasticityReturn().

◆ giveRealStressVector_3d()

FloatArrayF< 6 > oofem::MisesMatNl::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::MisesMat.

Definition at line 57 of file misesmatnl.C.

References OOFEM_ERROR.

◆ giveRemoteNonlocalStiffnessContribution()

void oofem::MisesMatNl::giveRemoteNonlocalStiffnessContribution	(	GaussPoint *	gp,
		IntArray &	rloc,
		const UnknownNumberingScheme &	s,
		FloatArray &	rcontrib,
		TimeStep *	tStep )

Computes the "remote" part of nonlocal stiffness contribution assembled for given integration point.

Parameters

gp	Remote integration point.
rloc	Remote element code numbers.
rcontrib	"Remote" contribution.
s	Numbering scheme that determines assembly.
tStep	Time step.

Definition at line 405 of file misesmatnl.C.

References oofem::FloatArray::at(), oofem::FloatArray::clear(), oofem::StructuralElement::computeBmatrixAt(), E, oofem::MisesMatStatus::giveCumulativePlasticStrain(), oofem::GaussPoint::giveElement(), oofem::Element::giveLocationArray(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::Material::giveStatus(), oofem::MisesMatStatus::giveTempCumulativePlasticStrain(), oofem::MisesMatStatus::giveTempEffectiveStress(), oofem::MisesMat::H, oofem::MisesMat::linearElasticMaterial, oofem::FloatArray::plusProduct(), oofem::FloatArray::resize(), oofem::sgn(), and oofem::FloatArray::zero().

Referenced by NonlocalMaterialStiffnessInterface_addIPContribution().

◆ hasBoundedSupport()

int oofem::MisesMatNl::hasBoundedSupport ( ) const

inlineoverridevirtual

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

Returns: True if weighting function bounded, zero otherwise.

Reimplemented from oofem::NonlocalMaterialExtensionInterface.

Definition at line 159 of file misesmatnl.h.

◆ initializeFrom()

void oofem::MisesMatNl::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::MisesMat.

Definition at line 282 of file misesmatnl.C.

References _IFT_MisesMatNl_averagingtype, _IFT_MisesMatNl_exp, _IFT_MisesMatNl_rf, averType, exponent, IR_GIVE_OPTIONAL_FIELD, and Rf.

◆ modifyNonlocalWeightFunctionAround()

void oofem::MisesMatNl::modifyNonlocalWeightFunctionAround ( GaussPoint * gp ) const

Definition at line 135 of file misesmatnl.C.

References oofem::FloatArray::at(), oofem::NonlocalMaterialExtensionInterface::cl, computeDistanceModifier(), oofem::Element::computeGlobalCoordinates(), oofem::Element::computeVolumeAround(), oofem::NonlocalMaterialExtensionInterface::computeWeightFunction(), oofem::distance(), oofem::GaussPoint::giveElement(), oofem::NonlocalMaterialExtensionInterface::giveIPIntegrationList(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::Material::giveStatus(), oofem::MisesMatStatus::giveTempDamage(), and oofem::NonlocalMaterialStatusExtensionInterface::setIntegrationScale().

Referenced by updateBeforeNonlocAverage().

◆ NonlocalMaterialStiffnessInterface_addIPContribution()

void oofem::MisesMatNl::NonlocalMaterialStiffnessInterface_addIPContribution	(	SparseMtrx &	dest,
		const UnknownNumberingScheme &	s,
		GaussPoint *	gp,
		TimeStep *	tStep )

overridevirtual

Computes and adds IP contributions to destination matrix.

Implements oofem::NonlocalMaterialStiffnessInterface.

Definition at line 341 of file misesmatnl.C.

References oofem::SparseMtrx::assemble(), oofem::FloatMatrix::clear(), oofem::Element::computeVolumeAround(), oofem::GaussPoint::giveElement(), oofem::NonlocalMaterialStatusExtensionInterface::giveIntegrationDomainList(), oofem::NonlocalMaterialStatusExtensionInterface::giveIntegrationScale(), giveLocalNonlocalStiffnessContribution(), giveRemoteNonlocalStiffnessContribution(), oofem::Material::giveStatus(), MisesMatNl(), and oofem::FloatMatrix::plusDyadUnsym().

◆ NonlocalMaterialStiffnessInterface_giveIntegrationDomainList()

std::vector< localIntegrationRecord > * oofem::MisesMatNl::NonlocalMaterialStiffnessInterface_giveIntegrationDomainList ( GaussPoint * gp )

overridevirtual

Returns integration list of receiver. Contains localIntegrationRecord structures, containing references to integration points and their weights that influence to nonlocal average in receiver's associated integration point.

Implements oofem::NonlocalMaterialStiffnessInterface.

Definition at line 372 of file misesmatnl.C.

References oofem::NonlocalMaterialExtensionInterface::buildNonlocalPointTable(), oofem::NonlocalMaterialStatusExtensionInterface::giveIntegrationDomainList(), and oofem::Material::giveStatus().

◆ packUnknowns()

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

overridevirtual

Pack all necessary data of integration point (according to element parallel_mode) into given communication buffer. The nature of packed data is material model dependent. Typically, for material of "local" response (response depends only on integration point local state) no data are exchanged. For "nonlocal" constitutive models the send/receive of local values which undergo averaging is performed between local and corresponding remote elements.

Parameters

buff	Communication buffer.
tStep	Solution step.
ip	Integration point.

Reimplemented from oofem::Material.

Definition at line 501 of file misesmatnl.C.

References oofem::NonlocalMaterialExtensionInterface::buildNonlocalPointTable(), oofem::MisesMatNlStatus::giveLocalCumPlasticStrainForAverage(), oofem::Material::giveStatus(), oofem::NonlocalMaterialExtensionInterface::updateDomainBeforeNonlocAverage(), and oofem::DataStream::write().

◆ unpackAndUpdateUnknowns()

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

overridevirtual

Unpack and updates all necessary data of given integration point (according to element parallel_mode) into given communication buffer.

See also: packUnknowns service.

Parameters

buff	Communication buffer.
tStep	Solution step.
ip	Integration point.

Reimplemented from oofem::Material.

Definition at line 513 of file misesmatnl.C.

References oofem::Material::giveStatus(), oofem::DataStream::read(), and oofem::MisesMatNlStatus::setLocalCumPlasticStrainForAverage().

◆ updateBeforeNonlocAverage()

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

overridevirtual

Declares the service updating local variables in given integration points, which take part in nonlocal average process. Because value of single integration point influences nonlocal variables in several near integration points, it is suitable to compute these variables only once. These should be stored in integration point associated statuses. The implementation is left on derived classes. Provide material local strain increment - as is provided to computeRealStresVector. This allows to update internal vars to be averaged to new state

Parameters

strainVector	Total strain vector in given integration point.
gp	Integration point to update.
tStep	Solution step indicating time of update.

Implements oofem::StructuralNonlocalMaterialExtensionInterface.

Definition at line 110 of file misesmatnl.C.

References averType, computeLocalCumPlasticStrain(), oofem::Material::giveStatus(), oofem::Material::initTempStatus(), modifyNonlocalWeightFunctionAround(), and oofem::MisesMat::performPlasticityReturn().

Member Data Documentation

◆ averType

int oofem::MisesMatNl::averType = 0

protected

Definition at line 94 of file misesmatnl.h.

Referenced by computeDistanceModifier(), giveInputRecord(), initializeFrom(), and updateBeforeNonlocAverage().

◆ exponent

double oofem::MisesMatNl::exponent = 1.

protected

Definition at line 93 of file misesmatnl.h.

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

◆ Rf

double oofem::MisesMatNl::Rf = 0.

protected

Definition at line 92 of file misesmatnl.h.

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

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

◆ MisesMatNl()

Member Function Documentation

◆ computeCumPlasticStrain()

◆ computeDamage()

◆ computeDistanceModifier()

◆ computeLocalCumPlasticStrain()

◆ CreateStatus()

◆ estimatePackSize()

◆ give1dStressStiffMtrx()

◆ giveClassName()

◆ giveInputRecord()

◆ giveInputRecordName()

◆ giveInterface()

◆ giveLocalNonlocalStiffnessContribution()

◆ giveRealStressVector_1d()

◆ giveRealStressVector_3d()

◆ giveRemoteNonlocalStiffnessContribution()

◆ hasBoundedSupport()

◆ initializeFrom()

◆ modifyNonlocalWeightFunctionAround()

◆ NonlocalMaterialStiffnessInterface_addIPContribution()

◆ NonlocalMaterialStiffnessInterface_giveIntegrationDomainList()

◆ packUnknowns()

◆ unpackAndUpdateUnknowns()

◆ updateBeforeNonlocAverage()

Member Data Documentation

◆ averType

◆ exponent

◆ Rf