Gradient bone damage-plastic material model. More...

#include <trabbonegrad3d.h>

Inheritance diagram for oofem::TrabBoneGrad3D:

Collaboration diagram for oofem::TrabBoneGrad3D:

Public Member Functions
	TrabBoneGrad3D (int n, Domain *d)

virtual	~TrabBoneGrad3D ()

virtual const char *	giveClassName () const

virtual const char *	giveInputRecordName () const

virtual IRResultType	initializeFrom (InputRecord *ir)
	Initializes receiver according to object description stored in input record. More...

virtual int	hasMaterialModeCapability (MaterialMode mode)
	Tests if material supports material mode. More...

virtual Interface *	giveInterface (InterfaceType t)
	Interface requesting service. More...

virtual void	givePDGradMatrix_uu (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Left upper block. More...

virtual void	givePDGradMatrix_ku (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Left lower block. More...

virtual void	givePDGradMatrix_uk (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Right upper block. More...

virtual void	givePDGradMatrix_kk (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Right lower block. More...

virtual void	givePDGradMatrix_LD (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Stress-based averaging. More...

virtual void	giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep)
	Computes the stiffness matrix for giveRealStressVector of receiver in given integration point, respecting its history. More...

virtual void	give3dMaterialStiffnessMatrix (FloatMatrix &answer, MatResponseMode, GaussPoint gp, TimeStep tStep)
	Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point. More...

void	give3dKappaMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

void	give3dGprime (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

void	giveInternalLength (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

virtual void	giveRealStressVectorGrad (FloatArray &answer1, double &answer2, GaussPoint gp, const FloatArray &totalStrain, double nonlocalCumulatedStrain, TimeStep tStep)
	gradient - based giveRealStressVector More...

virtual void	computeCumPlastStrain (double &kappa, GaussPoint gp, TimeStep tStep)

void	performPlasticityReturn (GaussPoint *gp, const FloatArray &totalStrain)

Public Member Functions inherited from oofem::TrabBone3D
	TrabBone3D (int n, Domain *d)

bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode)
	Returns true if stiffness matrix of receiver is symmetric Default implementation returns true. More...

double	evaluateCurrentYieldStress (const double kappa)

double	evaluateCurrentPlasticModulus (const double kappa)

double	evaluateCurrentViscousStress (const double deltaKappa, TimeStep *tStep)

double	evaluateCurrentViscousModulus (const double deltaKappa, TimeStep *tStep)

bool	projectOnYieldSurface (double &tempKappa, FloatArray &tempEffectiveStress, FloatArray &tempPlasDef, const FloatArray &trialEffectiveStress, const FloatMatrix &elasticity, const FloatMatrix &compliance, TrabBone3DStatus status, TimeStep tStep, GaussPoint *gp, int lineSearchFlag)

void	performPlasticityReturn (GaussPoint gp, const FloatArray &totalStrain, TimeStep tStep)

void	constructPlasFlowDirec (FloatArray &answer, double &norm, FloatMatrix &fabric, FloatArray &F, FloatArray &S)

void	constructDerivativeOfPlasFlowDirec (FloatMatrix &answer, FloatMatrix &fabric, FloatArray &F, FloatArray &S)

double	evaluatePlasCriterion (FloatMatrix &fabric, FloatArray &F, FloatArray &stress, double kappa, double deltaKappa, TimeStep *tStep)

double	computeDamageParam (double kappa)

double	computeDamageParamPrime (double kappa)

double	computeDamage (GaussPoint gp, TimeStep tStep)

void	computePlasStrainEnerDensity (GaussPoint *gp, const FloatArray &totalStrain, const FloatArray &totalStress)

void	computeDensificationStress (FloatArray &answer, GaussPoint gp, const FloatArray &totalStrain, TimeStep tStep)

void	constructAnisoComplTensor (FloatMatrix &answer)
	Construct anisotropic compliance tensor. More...

void	constructAnisoStiffnessTensor (FloatMatrix &answer)
	Construct anisotropic stiffness tensor. More...

void	constructAnisoFabricTensor (FloatMatrix &answer)
	Construct anisotropic fabric tensor. More...

void	constructAnisoFtensor (FloatArray &answer)

void	constructStiffnessTransformationMatrix (FloatMatrix &answer)

void	constructFabricTransformationMatrix (FloatMatrix &answer)

void	constructNormAdjustTensor (FloatMatrix &answer)
	Construct Tensor to adjust Norm. More...

virtual void	giveRealStressVector_3d (FloatArray &answer, GaussPoint , const FloatArray &, TimeStep )
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. More...

virtual int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)
	Returns the integration point corresponding value in Reduced form. 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...

Public Member Functions inherited from oofem::StructuralMaterial
	StructuralMaterial (int n, Domain *d)
	Constructor. More...

virtual	~StructuralMaterial ()
	Destructor. More...

virtual void	giveInputRecord (DynamicInputRecord &input)
	Setups the input record string of receiver. 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)

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

double	giveReferenceTemperature ()
	Returns the reference temperature of receiver. More...

virtual void	computeStressIndependentStrainVector (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode)
	Computes reduced strain vector in given integration point, generated by internal processes in material, which are independent on loading in particular integration point. More...

virtual void	computeStressIndependentStrainVector_3d (FloatArray &answer, GaussPoint gp, TimeStep tStep, ValueModeType mode)

virtual void	give3dMaterialStiffnessMatrix_dPdF (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

virtual void	give3dMaterialStiffnessMatrix_dCde (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

void	giveStressDependentPartOfStrainVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrainVector, TimeStep tStep, ValueModeType mode)
	Method for subtracting from reduced space strain vector its stress-independent parts (caused by temperature, shrinkage, creep and possibly by other phenomena). More...

void	giveStressDependentPartOfStrainVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrainVector, TimeStep tStep, ValueModeType mode)

virtual int	setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type)
	Sets the value of a certain variable at a given integration point to the given value. More...

virtual void	give2dBeamLayerStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 2d beam layer stiffness matrix of receiver. More...

virtual void	givePlateLayerStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 2d plate layer stiffness matrix of receiver. More...

virtual void	giveFiberStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 1d fiber stiffness matrix of receiver. More...

virtual void	give2dLatticeStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 2d lattice stiffness matrix of receiver. More...

virtual void	give3dLatticeStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing 3d lattice stiffness matrix of receiver. More...

virtual void	give2dPlateSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing stiffness matrix of plate subsoil model. More...

virtual void	give3dBeamSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
	Method for computing stiffness matrix of beam3d subsoil model. More...

virtual void	giveFirstPKStressVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. More...

virtual void	giveFirstPKStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
	Default implementation relies on giveFirstPKStressVector_3d. More...

virtual void	giveFirstPKStressVector_PlaneStress (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
	Default implementation relies on giveFirstPKStressVector_3d. More...

virtual void	giveFirstPKStressVector_1d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
	Default implementation relies on giveFirstPKStressVector_3d. More...

virtual void	giveCauchyStressVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)

virtual void	giveCauchyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)

virtual void	giveCauchyStressVector_PlaneStress (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)

virtual void	giveCauchyStressVector_1d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)

virtual void	givePlaneStressStiffMtrx (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 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...

std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros). More...

Protected Member Functions
MaterialStatus *	CreateStatus (GaussPoint *gp) const
	Creates new copy of associated status and inserts it into given integration point. More...

Protected Attributes
double	L

double	mParam

Protected Attributes inherited from oofem::TrabBone3D
double	m1

double	m2

double	rho

double	eps0

double	nu0

double	mu0

double	expk

double	expl

double	sig0Pos

double	sig0Neg

double	chi0Pos

double	chi0

double	chi0Neg

double	tau0

double	expq

double	expp

double	plasHardFactor

double	expPlasHard

double	expDam

double	critDam

int	printflag

int	max_num_iter

double	rel_yield_tol

double	strain_tol

double	x1
	Local coordinate system. More...

double	x2

double	x3

double	y1

double	y2

double	y3

double	z1

double	z2

double	z3

double	gammaL0
	Densificator properties. More...

double	gammaP0

double	tDens

double	densCrit

double	rL

double	rP

double	gammaL

double	gammaP

double	viscosity
	Viscosity parameter. 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...

Private Member Functions inherited from oofem::GradDpMaterialExtensionInterface
	GradDpMaterialExtensionInterface (Domain *d)
	Constructor. More...

virtual	~GradDpMaterialExtensionInterface ()
	Destructor. More...

virtual void	giveFirstPKStressVectorGrad (FloatArray &answer1, double &answer2, GaussPoint gp, const FloatArray &totalStrain, double nonlocalDamageDrivningVariable, TimeStep tStep)

virtual void	giveCauchyStressVectorGrad (FloatArray &answer1, double &answer2, GaussPoint gp, const FloatArray &totalStrain, double nonlocalDamageDrivningVariable, TimeStep tStep)

int	giveAveragingType ()

void	giveDistanceBasedCharacteristicLength (const FloatArray &gpCoords)
	Provides the distance based interaction radius This function is called when averType is set to 1. More...

Private Member Functions inherited from oofem::Interface
	Interface ()
	Constructor. More...

virtual	~Interface ()

Private Attributes inherited from oofem::GradDpMaterialExtensionInterface
Domain *	dom

double	cl0
	Initial(user defined) characteristic length of the nonlocal model (its interpretation depends on the weight function) Is different to cl when a Stress-based or a Distance-based nonlocal variation is applied. More...

int	averType
	Parameter which defines the averaging type When averType is equal to zereo classical approach is used for averType equal to one, distance-based apprach is used and avetType equal to two corresponds to stress-based averaging. More...

double	beta
	Parameter which multiplied with the interaction radius cl0 gives its minimum allowed value. More...

double	zeta
	Parameter used when Distance-based nonlocal variation is applied When it is multiplied with the interaction radius cl gives the maxinmum distance of the Gauss Point from the boundary. More...

double	cl
	Characteristic length of the nonlocal model (its interpretation depends on the type of weight function). More...

Detailed Description

Gradient bone damage-plastic material model.

Definition at line 81 of file trabbonegrad3d.h.

Constructor & Destructor Documentation

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

Definition at line 48 of file trabbonegrad3d.C.

References L.

oofem::TrabBoneGrad3D::~TrabBoneGrad3D ( )

virtual

Definition at line 53 of file trabbonegrad3d.C.

Member Function Documentation

void oofem::TrabBoneGrad3D::computeCumPlastStrain	(	double &	kappa,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Reimplemented from oofem::TrabBone3D.

Definition at line 315 of file trabbonegrad3d.C.

References oofem::TrabBoneGrad3DStatus::giveNonlocalCumulatedStrain(), oofem::Material::giveStatus(), oofem::TrabBone3DStatus::giveTempKappa(), and mParam.

MaterialStatus* oofem::TrabBoneGrad3D::CreateStatus ( GaussPoint * gp ) const

inlineprotectedvirtual

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

Definition at line 122 of file trabbonegrad3d.h.

References oofem::FEMComponent::domain, and oofem::TrabBoneGrad3DStatus::TrabBoneGrad3DStatus().

void oofem::TrabBoneGrad3D::give3dGprime	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

Definition at line 246 of file trabbonegrad3d.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::TrabBone3D::computeDamageParam(), oofem::TrabBone3D::computeDamageParamPrime(), oofem::TrabBone3DStatus::giveDam(), oofem::TrabBoneGrad3DStatus::giveNonlocalCumulatedStrain(), oofem::Material::giveStatus(), oofem::TrabBone3DStatus::giveTempEffectiveStress(), oofem::TrabBone3DStatus::giveTempKappa(), mParam, oofem::FloatMatrix::resize(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

Referenced by givePDGradMatrix_uk().

void oofem::TrabBoneGrad3D::give3dKappaMatrix	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

Definition at line 215 of file trabbonegrad3d.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::GradDpMaterialExtensionInterface::beta, oofem::FloatArray::beTProductOf(), oofem::TrabBone3DStatus::giveBeta(), oofem::TrabBone3DStatus::giveKappa(), oofem::TrabBone3DStatus::givePlasFlowDirec(), oofem::TrabBone3DStatus::giveSSaTensor(), oofem::Material::giveStatus(), oofem::TrabBone3DStatus::giveTempKappa(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

Referenced by givePDGradMatrix_ku().

void oofem::TrabBoneGrad3D::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::TrabBone3D.

Definition at line 132 of file trabbonegrad3d.C.

Referenced by givePDGradMatrix_uu().

virtual const char* oofem::TrabBoneGrad3D::giveClassName ( ) const

inlinevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::TrabBone3D.

Definition at line 91 of file trabbonegrad3d.h.

virtual const char* oofem::TrabBoneGrad3D::giveInputRecordName ( ) const

inlinevirtual

Returns: Input record name of the receiver.

Reimplemented from oofem::TrabBone3D.

Definition at line 92 of file trabbonegrad3d.h.

References _IFT_TrabBoneGrad3D_Name, and oofem::MaterialStatus::initializeFrom().

virtual Interface* oofem::TrabBoneGrad3D::giveInterface ( InterfaceType t )

inlinevirtual

Interface requesting service.

See also: InterfaceType

Returns: Requested interface if implemented, otherwise NULL.

Reimplemented from oofem::FEMComponent.

Definition at line 96 of file trabbonegrad3d.h.

References oofem::IntegrationPointStatus::gp, oofem::GradDpMaterialExtensionInterfaceType, oofem::TrabBone3DStatus::kappa, and oofem::GradDpMaterialStatusExtensionInterface::nonlocalCumulatedStrain.

void oofem::TrabBoneGrad3D::giveInternalLength	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

Definition at line 277 of file trabbonegrad3d.C.

References oofem::FloatMatrix::at(), L, and oofem::FloatMatrix::resize().

Referenced by givePDGradMatrix_kk().

void oofem::TrabBoneGrad3D::givePDGradMatrix_kk	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Right lower block.

Implements oofem::GradDpMaterialExtensionInterface.

Definition at line 111 of file trabbonegrad3d.C.

References oofem::__MaterialModeToString(), giveInternalLength(), oofem::GaussPoint::giveMaterialMode(), and OOFEM_ERROR.

void oofem::TrabBoneGrad3D::givePDGradMatrix_ku	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Left lower block.

Implements oofem::GradDpMaterialExtensionInterface.

Definition at line 85 of file trabbonegrad3d.C.

References oofem::__MaterialModeToString(), give3dKappaMatrix(), oofem::GaussPoint::giveMaterialMode(), and OOFEM_ERROR.

void oofem::TrabBoneGrad3D::givePDGradMatrix_LD	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Stress-based averaging.

Implements oofem::GradDpMaterialExtensionInterface.

Definition at line 124 of file trabbonegrad3d.C.

References oofem::__MaterialModeToString(), oofem::GaussPoint::giveMaterialMode(), and OOFEM_ERROR.

void oofem::TrabBoneGrad3D::givePDGradMatrix_uk	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Right upper block.

Implements oofem::GradDpMaterialExtensionInterface.

Definition at line 98 of file trabbonegrad3d.C.

References oofem::__MaterialModeToString(), give3dGprime(), oofem::GaussPoint::giveMaterialMode(), and OOFEM_ERROR.

void oofem::TrabBoneGrad3D::givePDGradMatrix_uu	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Left upper block.

Implements oofem::GradDpMaterialExtensionInterface.

Definition at line 72 of file trabbonegrad3d.C.

References oofem::__MaterialModeToString(), give3dMaterialStiffnessMatrix(), oofem::GaussPoint::giveMaterialMode(), and OOFEM_ERROR.

void oofem::TrabBoneGrad3D::giveRealStressVectorGrad	(	FloatArray &	answer1,
		double &	answer2,
		GaussPoint *	gp,
		const FloatArray &	totalStrain,
		double	nonlocalDamageDrivningVariable,
		TimeStep *	tStep
	)

virtual

gradient - based giveRealStressVector

Reimplemented from oofem::GradDpMaterialExtensionInterface.

Definition at line 284 of file trabbonegrad3d.C.

References oofem::FloatArray::add(), oofem::FloatArray::beScaled(), oofem::TrabBone3D::computeDamage(), oofem::TrabBone3D::computeDensificationStress(), oofem::TrabBone3D::densCrit, oofem::Material::giveStatus(), oofem::TrabBone3DStatus::giveTempEffectiveStress(), oofem::TrabBone3DStatus::giveTempKappa(), oofem::Material::initTempStatus(), oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), oofem::TrabBone3D::performPlasticityReturn(), oofem::TrabBoneGrad3DStatus::setNonlocalCumulatedStrain(), oofem::TrabBone3DStatus::setTempDam(), and oofem::TrabBone3DStatus::setTempEffectiveStress().

void oofem::TrabBoneGrad3D::giveStiffnessMatrix	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Computes the stiffness matrix for giveRealStressVector of receiver in given integration point, respecting its history.

The algorithm should use temporary or equilibrium history variables stored in integration point status to compute and return required result.

Parameters

answer	Contains result.
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 62 of file trabbonegrad3d.C.

References OOFEM_ERROR.

int oofem::TrabBoneGrad3D::hasMaterialModeCapability ( MaterialMode mode )

virtual

Tests if material supports material mode.

Parameters

mode	Required material mode.

Returns: Nonzero if supported, zero otherwise.

Reimplemented from oofem::StructuralMaterial.

Definition at line 57 of file trabbonegrad3d.C.

IRResultType oofem::TrabBoneGrad3D::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

Todo:: Why not use vectors?

Reimplemented from oofem::TrabBone3D.

Definition at line 325 of file trabbonegrad3d.C.

References _IFT_TrabBoneGrad3D_L, _IFT_TrabBoneGrad3D_m, oofem::TrabBone3D::initializeFrom(), IR_GIVE_OPTIONAL_FIELD, L, and mParam.

void oofem::TrabBoneGrad3D::performPlasticityReturn	(	GaussPoint *	gp,
		const FloatArray &	totalStrain
	)

Member Data Documentation

double oofem::TrabBoneGrad3D::L

protected

Definition at line 84 of file trabbonegrad3d.h.

Referenced by giveInternalLength(), initializeFrom(), and TrabBoneGrad3D().

double oofem::TrabBoneGrad3D::mParam

protected

Definition at line 85 of file trabbonegrad3d.h.

Referenced by computeCumPlastStrain(), give3dGprime(), give3dMaterialStiffnessMatrix(), 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

Member Function Documentation

Member Data Documentation