#include <trabbonegrad3d.h>

Inheritance diagram for oofem::TrabBoneGrad3D:

Collaboration diagram for oofem::TrabBoneGrad3D:

Public Member Functions
	TrabBoneGrad3D (int n, Domain *d)
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
void	initializeFrom (InputRecord &ir) override
bool	hasMaterialModeCapability (MaterialMode mode) const override
Interface *	giveInterface (InterfaceType t) override
void	giveGradientDamageStiffnessMatrix_uu (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep) override
	Left upper block.
void	giveGradientDamageStiffnessMatrix_du (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep) override
	Left lower block.
void	giveGradientDamageStiffnessMatrix_ud (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep) override
	Right upper block.
void	giveGradientDamageStiffnessMatrix_dd_BB (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep) override
void	giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode, GaussPoint gp, TimeStep tStep) const override
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)
void	computeLocalDamageDrivingVariable (double &answer, GaussPoint gp, TimeStep tStep) override
void	giveNonlocalInternalForces_N_factor (double &answer, double nlddv, GaussPoint gp, TimeStep tStep) override
void	giveNonlocalInternalForces_B_factor (FloatArray &answer, const FloatArray &nlddv, GaussPoint gp, TimeStep tStep) override
void	giveRealStressVectorGradientDamage (FloatArray &answer1, double &answer2, GaussPoint gp, const FloatArray &totalStrain, double nonlocalCumulatedStrain, TimeStep tStep) override
	gradient - based giveRealStressVector
double	computeCumPlastStrain (GaussPoint gp, TimeStep tStep) const override
void	performPlasticityReturn (GaussPoint *gp, const FloatArray &totalStrain)
Public Member Functions inherited from oofem::TrabBone3D
	TrabBone3D (int n, Domain *d)
bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode) const override
double	evaluateCurrentYieldStress (double kappa) const
double	evaluateCurrentPlasticModulus (double kappa) const
double	evaluateCurrentViscousStress (double deltaKappa, TimeStep *tStep) const
double	evaluateCurrentViscousModulus (double deltaKappa, TimeStep *tStep) const
bool	projectOnYieldSurface (double &tempKappa, FloatArrayF< 6 > &tempEffectiveStress, FloatArrayF< 6 > &tempPlasDef, const FloatArrayF< 6 > &trialEffectiveStress, const FloatMatrixF< 6, 6 > &elasticity, const FloatMatrixF< 6, 6 > &compliance, TrabBone3DStatus status, TimeStep tStep, GaussPoint *gp, int lineSearchFlag) const
void	performPlasticityReturn (GaussPoint gp, const FloatArrayF< 6 > &strain, TimeStep tStep) const
std::pair< FloatArrayF< 6 >, double >	constructPlasFlowDirec (FloatMatrixF< 6, 6 > &fabric, const FloatArrayF< 6 > &F, const FloatArrayF< 6 > &S) const
FloatMatrixF< 6, 6 >	constructDerivativeOfPlasFlowDirec (const FloatMatrixF< 6, 6 > &fabric, const FloatArrayF< 6 > &F, const FloatArrayF< 6 > &S) const
double	evaluatePlasCriterion (const FloatMatrixF< 6, 6 > &fabric, const FloatArrayF< 6 > &F, const FloatArrayF< 6 > &stress, double kappa, double deltaKappa, TimeStep *tStep) const
double	computeDamageParam (double kappa) const
double	computeDamageParamPrime (double kappa) const
double	computeDamage (GaussPoint gp, TimeStep tStep) const
void	computePlasStrainEnerDensity (GaussPoint *gp, const FloatArrayF< 6 > &strain, const FloatArrayF< 6 > &stress) const
FloatArrayF< 6 >	computeDensificationStress (GaussPoint gp, const FloatArrayF< 6 > &totalStrain, TimeStep tStep) const
FloatMatrixF< 6, 6 >	constructAnisoComplTensor () const
	Construct anisotropic compliance tensor.
FloatMatrixF< 6, 6 >	constructAnisoStiffnessTensor () const
	Construct anisotropic stiffness tensor.
FloatMatrixF< 6, 6 >	constructAnisoFabricTensor () const
	Construct anisotropic fabric tensor.
FloatArrayF< 6 >	constructAnisoFtensor () const
FloatMatrixF< 6, 6 >	constructStiffnessTransformationMatrix () const
FloatMatrixF< 6, 6 >	constructFabricTransformationMatrix () const
FloatMatrixF< 6, 6 >	constructNormAdjustTensor () const
	Construct Tensor to adjust Norm.
FloatArrayF< 6 >	giveRealStressVector_3d (const FloatArrayF< 6 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
double	predictRelativeComputationalCost (GaussPoint *gp) override
double	predictRelativeRedistributionCost (GaussPoint *gp) override
Public Member Functions inherited from oofem::StructuralMaterial
	StructuralMaterial (int n, Domain *d)
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	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)
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
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 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).

Protected Member Functions
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override

Protected Attributes
double	L = 0.
double	mParam = 0.
Protected Attributes inherited from oofem::TrabBone3D
double	m1 = 0.
double	m2 = 0.
double	rho = 0.
double	eps0 = 0.
double	nu0 = 0.
double	mu0 = 0.
double	expk = 0.
double	expl = 0.
double	sig0Pos = 0.
double	sig0Neg = 0.
double	chi0Pos = 0.
double	chi0 = 0.
double	chi0Neg = 0.
double	tau0 = 0.
double	expq = 0.
double	expp = 0.
double	plasHardFactor = 0.
double	expPlasHard = 0.
double	expDam = 0.
double	critDam = 0.
int	printflag = 0
int	max_num_iter = 0
double	rel_yield_tol = 0.
double	strain_tol = 0.
double	x1 = 1.
	Local coordinate system.
double	x2 = 0.
double	x3 = 0.
double	y1 = 0.
double	y2 = 1.
double	y3 = 0.
double	z1 = 0.
double	z2 = 0.
double	z3 = 1.
double	gammaL0 = 0.
	Densificator properties.
double	gammaP0 = 0.
double	tDens = 0.
double	densCrit = 0.
double	rL = 0.
double	rP = 0.
double	gammaL = 0.
double	gammaP = 0.
double	viscosity = 0.
	Viscosity parameter.
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.
Private Member Functions inherited from oofem::GradientDamageMaterialExtensionInterface
	GradientDamageMaterialExtensionInterface (Domain *d)
virtual	~GradientDamageMaterialExtensionInterface ()
	Destructor.
virtual void	giveGradientDamageStiffnessMatrix_dd_NN (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Right lower block.
virtual void	giveGradientDamageStiffnessMatrix_dd_BN (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
virtual void	giveFirstPKStressVectorGradientDamage (FloatArray &answer1, double &answer2, GaussPoint gp, const FloatArray &totalStrain, double nonlocalDamageDrivningVariable, TimeStep tStep)
virtual void	giveCauchyStressVectorGradientDamage (FloatArray &answer1, double &answer2, GaussPoint gp, const FloatArray &totalStrain, double nonlocalDamageDrivningVariable, TimeStep tStep)
virtual void	computeInternalForcesRegularizationTerm (double &answer, GaussPoint gp, TimeStep tStep)
virtual void	computeStiffnessRegularizationTerm (double &answer, GaussPoint gp, TimeStep tStep)
Private Member Functions inherited from oofem::Interface
	Interface ()
	Constructor.
virtual	~Interface ()
Private Attributes inherited from oofem::GradientDamageMaterialExtensionInterface
Domain *	dom = nullptr
double	internalLength = 0.

Detailed Description

Gradient bone damage-plastic material model.

Definition at line 79 of file trabbonegrad3d.h.

Constructor & Destructor Documentation

◆ TrabBoneGrad3D()

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

Definition at line 48 of file trabbonegrad3d.C.

References oofem::GradientDamageMaterialExtensionInterface::GradientDamageMaterialExtensionInterface(), and oofem::TrabBone3D::TrabBone3D().

Member Function Documentation

◆ computeCumPlastStrain()

double oofem::TrabBoneGrad3D::computeCumPlastStrain	(	GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Reimplemented from oofem::TrabBone3D.

Definition at line 297 of file trabbonegrad3d.C.

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

◆ computeLocalDamageDrivingVariable()

void oofem::TrabBoneGrad3D::computeLocalDamageDrivingVariable	(	double &	answer,
		GaussPoint *	gp,
		TimeStep *	tStep )

overridevirtual

Implements oofem::GradientDamageMaterialExtensionInterface.

Definition at line 135 of file trabbonegrad3d.C.

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

◆ CreateStatus()

std::unique_ptr< MaterialStatus > oofem::TrabBoneGrad3D::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::TrabBone3D.

Definition at line 123 of file trabbonegrad3d.h.

◆ give3dGprime()

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

Definition at line 237 of file trabbonegrad3d.C.

References oofem::FloatMatrix::at(), oofem::Material::giveStatus(), oofem::TrabBone3DStatus::giveTempKappa(), mParam, oofem::FloatMatrix::resize(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

Referenced by giveGradientDamageStiffnessMatrix_ud().

◆ give3dKappaMatrix()

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

Definition at line 209 of file trabbonegrad3d.C.

References oofem::FloatMatrix::at(), oofem::TrabBone3DStatus::giveKappa(), oofem::Material::giveStatus(), oofem::FloatMatrix::resize(), oofem::Tdot(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().

Referenced by giveGradientDamageStiffnessMatrix_du().

◆ give3dMaterialStiffnessMatrix()

FloatMatrixF< 6, 6 > oofem::TrabBoneGrad3D::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::TrabBone3D.

Definition at line 143 of file trabbonegrad3d.C.

References oofem::TrabBone3D::constructAnisoComplTensor(), oofem::dot(), oofem::dyad(), oofem::TrabBone3D::gammaL0, oofem::TrabBone3D::gammaP0, oofem::Material::giveStatus(), oofem::I6_I6, oofem::inv(), mParam, oofem::TrabBone3D::printflag, oofem::TrabBone3D::rho, oofem::TrabBone3D::rL, oofem::TrabBone3D::rP, oofem::TrabBone3D::tDens, and oofem::Tdot().

Referenced by giveGradientDamageStiffnessMatrix_uu().

◆ giveClassName()

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

inlineoverridevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::TrabBone3D.

Definition at line 88 of file trabbonegrad3d.h.

◆ giveGradientDamageStiffnessMatrix_dd_BB()

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

overridevirtual

Implements oofem::GradientDamageMaterialExtensionInterface.

Definition at line 108 of file trabbonegrad3d.C.

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

◆ giveGradientDamageStiffnessMatrix_du()

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

overridevirtual

Left lower block.

Implements oofem::GradientDamageMaterialExtensionInterface.

Definition at line 81 of file trabbonegrad3d.C.

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

◆ giveGradientDamageStiffnessMatrix_ud()

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

overridevirtual

Right upper block.

Implements oofem::GradientDamageMaterialExtensionInterface.

Definition at line 94 of file trabbonegrad3d.C.

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

◆ giveGradientDamageStiffnessMatrix_uu()

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

overridevirtual

Left upper block.

Implements oofem::GradientDamageMaterialExtensionInterface.

Definition at line 68 of file trabbonegrad3d.C.

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

◆ giveInputRecordName()

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

inlineoverridevirtual

Returns: Input record name of the receiver.

Reimplemented from oofem::TrabBone3D.

Definition at line 89 of file trabbonegrad3d.h.

References _IFT_TrabBoneGrad3D_Name.

◆ giveInterface()

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

inlineoverridevirtual

Interface requesting service.

See also: InterfaceType

Returns: Requested interface if implemented, otherwise NULL.

Reimplemented from oofem::FEMComponent.

Definition at line 93 of file trabbonegrad3d.h.

References oofem::GradientDamageMaterialExtensionInterface::GradientDamageMaterialExtensionInterface(), and oofem::GradientDamageMaterialExtensionInterfaceType.

◆ giveInternalLength()

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

Definition at line 263 of file trabbonegrad3d.C.

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

Referenced by giveGradientDamageStiffnessMatrix_dd_BB().

◆ giveNonlocalInternalForces_B_factor()

void oofem::TrabBoneGrad3D::giveNonlocalInternalForces_B_factor	(	FloatArray &	answer,
		const FloatArray &	nlddv,
		GaussPoint *	gp,
		TimeStep *	tStep )

overridevirtual

Implements oofem::GradientDamageMaterialExtensionInterface.

Definition at line 127 of file trabbonegrad3d.C.

References oofem::GradientDamageMaterialExtensionInterface::internalLength, and oofem::FloatArray::times().

◆ giveNonlocalInternalForces_N_factor()

void oofem::TrabBoneGrad3D::giveNonlocalInternalForces_N_factor	(	double &	answer,
		double	nlddv,
		GaussPoint *	gp,
		TimeStep *	tStep )

overridevirtual

Implements oofem::GradientDamageMaterialExtensionInterface.

Definition at line 121 of file trabbonegrad3d.C.

◆ giveRealStressVectorGradientDamage()

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

overridevirtual

gradient - based giveRealStressVector

Reimplemented from oofem::GradientDamageMaterialExtensionInterface.

Definition at line 270 of file trabbonegrad3d.C.

References oofem::FloatArray::add(), oofem::TrabBone3D::computeDamage(), oofem::TrabBone3D::computeDensificationStress(), oofem::TrabBone3D::densCrit, oofem::Material::giveStatus(), and oofem::Material::initTempStatus().

◆ giveStiffnessMatrix()

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

overridevirtual

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 58 of file trabbonegrad3d.C.

References OOFEM_ERROR.

◆ hasMaterialModeCapability()

bool oofem::TrabBoneGrad3D::hasMaterialModeCapability ( MaterialMode mode ) const

overridevirtual

Tests if material supports material mode.

Parameters

mode	Required material mode.

Returns: Nonzero if supported, zero otherwise.

Reimplemented from oofem::StructuralMaterial.

Definition at line 53 of file trabbonegrad3d.C.

◆ initializeFrom()

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

Todo: Why not use vectors?

Reimplemented from oofem::TrabBone3D.

Definition at line 307 of file trabbonegrad3d.C.

References _IFT_TrabBoneGrad3D_L, _IFT_TrabBoneGrad3D_m, IR_GIVE_OPTIONAL_FIELD, L, and mParam.

◆ performPlasticityReturn()

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

Member Data Documentation

◆ L

double oofem::TrabBoneGrad3D::L = 0.

protected

Definition at line 82 of file trabbonegrad3d.h.

Referenced by giveInternalLength(), and initializeFrom().

◆ mParam

double oofem::TrabBoneGrad3D::mParam = 0.

protected

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

◆ TrabBoneGrad3D()

Member Function Documentation

◆ computeCumPlastStrain()

◆ computeLocalDamageDrivingVariable()

◆ CreateStatus()

◆ give3dGprime()

◆ give3dKappaMatrix()

◆ give3dMaterialStiffnessMatrix()

◆ giveClassName()

◆ giveGradientDamageStiffnessMatrix_dd_BB()

◆ giveGradientDamageStiffnessMatrix_du()

◆ giveGradientDamageStiffnessMatrix_ud()

◆ giveGradientDamageStiffnessMatrix_uu()

◆ giveInputRecordName()

◆ giveInterface()

◆ giveInternalLength()

◆ giveNonlocalInternalForces_B_factor()

◆ giveNonlocalInternalForces_N_factor()

◆ giveRealStressVectorGradientDamage()

◆ giveStiffnessMatrix()

◆ hasMaterialModeCapability()

◆ initializeFrom()

◆ performPlasticityReturn()

Member Data Documentation

◆ L

◆ mParam