OOFEM 3.0
Loading...
Searching...
No Matches
oofem::TrabBone3D Class Reference

#include <trabbone3d.h>

Inheritance diagram for oofem::TrabBone3D:
Collaboration diagram for oofem::TrabBone3D:

Public Member Functions

 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
virtual double computeCumPlastStrain (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.
FloatMatrixF< 6, 6 > give3dMaterialStiffnessMatrix (MatResponseMode, GaussPoint *gp, TimeStep *tStep) const override
FloatArrayF< 6 > giveRealStressVector_3d (const FloatArrayF< 6 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
const char * giveInputRecordName () const override
const char * giveClassName () const override
void initializeFrom (InputRecord &ir) override
std::unique_ptr< MaterialStatusCreateStatus (GaussPoint *gp) const override
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)
bool hasMaterialModeCapability (MaterialMode mode) const override
void giveInputRecord (DynamicInputRecord &input) override
void giveCharacteristicMatrix (FloatMatrix &answer, MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const override
 Returns characteristic matrix of the receiver.
void giveCharacteristicVector (FloatArray &answer, FloatArray &flux, MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const override
 Returns characteristic vector of the receiver.
virtual void giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const
virtual void giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) const
virtual FloatArrayF< 4 > giveRealStressVector_PlaneStrain (const FloatArrayF< 4 > &strain, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_3d.
virtual FloatArray giveRealStressVector_StressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint *gp, TimeStep *tStep) const
 Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArray giveRealStressVector_ShellStressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArrayF< 3 > giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 1 > giveRealStressVector_1d (const FloatArrayF< 1 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 2 > giveRealStressVector_Warping (const FloatArrayF< 2 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 2 > giveRealStressVector_2dBeamLayer (const FloatArrayF< 2 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 5 > giveRealStressVector_PlateLayer (const FloatArrayF< 5 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 > giveRealStressVector_Fiber (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 > giveRealStressVector_2dPlateSubSoil (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation is not provided.
virtual FloatArrayF< 6 > giveRealStressVector_3dBeamSubSoil (const FloatArrayF< 6 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArrayF< 9 > giveFirstPKStressVector_3d (const FloatArrayF< 9 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
virtual FloatArrayF< 5 > giveFirstPKStressVector_PlaneStrain (const FloatArrayF< 5 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveFirstPKStressVector_3d.
virtual FloatArray giveFirstPKStressVector_StressControl (const FloatArray &reducedvF, const IntArray &FControl, GaussPoint *gp, TimeStep *tStep) const
 Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArrayF< 4 > giveFirstPKStressVector_PlaneStress (const FloatArrayF< 4 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveFirstPKStressVector_StressControl.
virtual FloatArrayF< 1 > giveFirstPKStressVector_1d (const FloatArrayF< 1 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveFirstPKStressVector_StressControl.
virtual void giveCauchyStressVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void giveCauchyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void giveCauchyStressVector_PlaneStress (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void giveCauchyStressVector_1d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
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 MaterialStatusgiveStatus (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.
DomaingiveDomain () const
virtual void setDomain (Domain *d)
int giveNumber () const
void setNumber (int num)
virtual void updateLocalNumbering (EntityRenumberingFunctor &f)
virtual void initializeFrom (InputRecord &ir, int priority)
virtual void initializeFinish ()
virtual void postInitialize ()
 Performs post initialization steps. Called after all components are created and initialized.
virtual void printOutputAt (FILE *file, TimeStep *tStep)
virtual InterfacegiveInterface (InterfaceType t)
std::string errorInfo (const char *func) const
 Returns string for prepending output (used by error reporting macros).

Protected Attributes

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.
Domaindomain
 Link to domain object, useful for communicating with other FEM components.

Additional Inherited Members

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

Detailed Description

Definition at line 169 of file trabbone3d.h.

Constructor & Destructor Documentation

◆ TrabBone3D()

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

Member Function Documentation

◆ computeCumPlastStrain()

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

◆ computeDamage()

◆ computeDamageParam()

double oofem::TrabBone3D::computeDamageParam ( double kappa) const

Definition at line 482 of file trabbone3d.C.

References critDam, and expDam.

Referenced by computeDamage().

◆ computeDamageParamPrime()

double oofem::TrabBone3D::computeDamageParamPrime ( double kappa) const

Definition at line 496 of file trabbone3d.C.

References critDam, and expDam.

◆ computeDensificationStress()

FloatArrayF< 6 > oofem::TrabBone3D::computeDensificationStress ( GaussPoint * gp,
const FloatArrayF< 6 > & totalStrain,
TimeStep * tStep ) const

◆ computePlasStrainEnerDensity()

void oofem::TrabBone3D::computePlasStrainEnerDensity ( GaussPoint * gp,
const FloatArrayF< 6 > & strain,
const FloatArrayF< 6 > & stress ) const

◆ constructAnisoComplTensor()

◆ constructAnisoFabricTensor()

FloatMatrixF< 6, 6 > oofem::TrabBone3D::constructAnisoFabricTensor ( ) const

Construct anisotropic fabric tensor.

Definition at line 669 of file trabbone3d.C.

References oofem::FloatMatrixF< N, M >::at(), chi0, constructFabricTransformationMatrix(), expp, expq, m1, m2, rho, oofem::rotate(), sig0Neg, sig0Pos, and tau0.

Referenced by projectOnYieldSurface().

◆ constructAnisoFtensor()

FloatArrayF< 6 > oofem::TrabBone3D::constructAnisoFtensor ( ) const

◆ constructAnisoStiffnessTensor()

FloatMatrixF< 6, 6 > oofem::TrabBone3D::constructAnisoStiffnessTensor ( ) const

Construct anisotropic stiffness tensor.

Definition at line 617 of file trabbone3d.C.

References oofem::FloatMatrixF< N, M >::at(), constructStiffnessTransformationMatrix(), eps0, expk, expl, m1, m2, mu0, nu0, rho, and oofem::rotate().

◆ constructDerivativeOfPlasFlowDirec()

FloatMatrixF< 6, 6 > oofem::TrabBone3D::constructDerivativeOfPlasFlowDirec ( const FloatMatrixF< 6, 6 > & fabric,
const FloatArrayF< 6 > & F,
const FloatArrayF< 6 > & S ) const

Definition at line 435 of file trabbone3d.C.

References constructNormAdjustTensor(), oofem::dot(), oofem::dyad(), oofem::norm(), and S.

Referenced by projectOnYieldSurface().

◆ constructFabricTransformationMatrix()

FloatMatrixF< 6, 6 > oofem::TrabBone3D::constructFabricTransformationMatrix ( ) const

◆ constructNormAdjustTensor()

FloatMatrixF< 6, 6 > oofem::TrabBone3D::constructNormAdjustTensor ( ) const

Construct Tensor to adjust Norm.

Definition at line 768 of file trabbone3d.C.

References oofem::FloatMatrixF< N, M >::at().

Referenced by constructDerivativeOfPlasFlowDirec(), constructPlasFlowDirec(), and projectOnYieldSurface().

◆ constructPlasFlowDirec()

std::pair< FloatArrayF< 6 >, double > oofem::TrabBone3D::constructPlasFlowDirec ( FloatMatrixF< 6, 6 > & fabric,
const FloatArrayF< 6 > & F,
const FloatArrayF< 6 > & S ) const

Definition at line 415 of file trabbone3d.C.

References constructNormAdjustTensor(), oofem::dot(), oofem::norm(), and S.

Referenced by projectOnYieldSurface().

◆ constructStiffnessTransformationMatrix()

FloatMatrixF< 6, 6 > oofem::TrabBone3D::constructStiffnessTransformationMatrix ( ) const

Definition at line 717 of file trabbone3d.C.

References oofem::FloatMatrixF< N, M >::at(), x1, x2, x3, y1, y2, y3, z1, z2, and z3.

Referenced by constructAnisoComplTensor(), and constructAnisoStiffnessTensor().

◆ CreateStatus()

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

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

Parameters
gpIntegration point where newly created status will be stored.
Returns
Reference to new status.

Reimplemented from oofem::Material.

Reimplemented in oofem::TrabBoneGrad3D, and oofem::TrabBoneNL3D.

Definition at line 1121 of file trabbone3d.C.

◆ evaluateCurrentPlasticModulus()

double oofem::TrabBone3D::evaluateCurrentPlasticModulus ( double kappa) const

Definition at line 144 of file trabbone3d.C.

References expPlasHard, and plasHardFactor.

Referenced by projectOnYieldSurface().

◆ evaluateCurrentViscousModulus()

double oofem::TrabBone3D::evaluateCurrentViscousModulus ( double deltaKappa,
TimeStep * tStep ) const

Definition at line 165 of file trabbone3d.C.

References oofem::TimeStep::giveTimeIncrement(), and viscosity.

Referenced by projectOnYieldSurface().

◆ evaluateCurrentViscousStress()

double oofem::TrabBone3D::evaluateCurrentViscousStress ( double deltaKappa,
TimeStep * tStep ) const

Definition at line 151 of file trabbone3d.C.

References oofem::TimeStep::giveTimeIncrement(), and viscosity.

Referenced by evaluatePlasCriterion().

◆ evaluateCurrentYieldStress()

double oofem::TrabBone3D::evaluateCurrentYieldStress ( double kappa) const

Definition at line 138 of file trabbone3d.C.

References expPlasHard, and plasHardFactor.

Referenced by evaluatePlasCriterion(), and projectOnYieldSurface().

◆ evaluatePlasCriterion()

double oofem::TrabBone3D::evaluatePlasCriterion ( const FloatMatrixF< 6, 6 > & fabric,
const FloatArrayF< 6 > & F,
const FloatArrayF< 6 > & stress,
double kappa,
double deltaKappa,
TimeStep * tStep ) const

◆ give3dMaterialStiffnessMatrix()

FloatMatrixF< 6, 6 > oofem::TrabBone3D::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
answerComputed results.
modeMaterial response mode.
gpIntegration point.
tStepTime step (most models are able to respond only when tStep is current time step).

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::TrabBoneGrad3D, and oofem::TrabBoneNL3D.

Definition at line 78 of file trabbone3d.C.

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

◆ giveClassName()

const char * oofem::TrabBone3D::giveClassName ( ) const
inlineoverridevirtual
Returns
Class name of the receiver.

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::TrabBoneGrad3D, and oofem::TrabBoneNL3D.

Definition at line 242 of file trabbone3d.h.

◆ giveInputRecordName()

const char * oofem::TrabBone3D::giveInputRecordName ( ) const
inlineoverridevirtual
Returns
Input record name of the receiver.

Implements oofem::FEMComponent.

Reimplemented in oofem::TrabBoneGrad3D, and oofem::TrabBoneNL3D.

Definition at line 241 of file trabbone3d.h.

References _IFT_TrabBone3D_Name.

◆ giveIPValue()

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

Returns the integration point corresponding value in Reduced form.

Parameters
answerContain corresponding ip value, zero sized if not available.
gpIntegration point to which the value refers.
typeDetermines the type of internal variable.
tStepDetermines the time step.
Returns
Nonzero if the assignment can be done, zero if this type of variable is not supported.

Reimplemented from oofem::StructuralMaterial.

Definition at line 924 of file trabbone3d.C.

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

◆ giveRealStressVector_3d()

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

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

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::TrabBoneNL3D.

Definition at line 547 of file trabbone3d.C.

References computeDamage(), computeDensificationStress(), computePlasStrainEnerDensity(), densCrit, oofem::Material::giveStatus(), oofem::Material::initTempStatus(), and performPlasticityReturn().

◆ initializeFrom()

void oofem::TrabBone3D::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
irInput record to initialize from.
priorityPriority of the input record. This is used to determine the order of initialization
Todo
Why not use vectors?

Reimplemented from oofem::StructuralMaterial.

Reimplemented in oofem::TrabBoneGrad3D, and oofem::TrabBoneNL3D.

Definition at line 835 of file trabbone3d.C.

References _IFT_TrabBone3D_chi0Pos, _IFT_TrabBone3D_critDam, _IFT_TrabBone3D_eps0, _IFT_TrabBone3D_expDam, _IFT_TrabBone3D_expk, _IFT_TrabBone3D_expl, _IFT_TrabBone3D_expp, _IFT_TrabBone3D_expPlasHard, _IFT_TrabBone3D_expq, _IFT_TrabBone3D_m1, _IFT_TrabBone3D_m2, _IFT_TrabBone3D_max_num_iter, _IFT_TrabBone3D_mu0, _IFT_TrabBone3D_nu0, _IFT_TrabBone3D_plasHardFactor, _IFT_TrabBone3D_printflag, _IFT_TrabBone3D_rel_yield_tol, _IFT_TrabBone3D_rho, _IFT_TrabBone3D_sig0Neg, _IFT_TrabBone3D_sig0Pos, _IFT_TrabBone3D_strain_tol, _IFT_TrabBone3D_tau0, _IFT_TrabBone3D_viscosity, _IFT_TrabBone3D_x1, _IFT_TrabBone3D_x2, _IFT_TrabBone3D_x3, _IFT_TrabBone3D_y1, _IFT_TrabBone3D_y2, _IFT_TrabBone3D_y3, chi0, chi0Neg, chi0Pos, critDam, eps0, expDam, expk, expl, expp, expPlasHard, expq, IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, m1, m2, max_num_iter, mu0, nu0, plasHardFactor, printflag, rel_yield_tol, rho, sig0Neg, sig0Pos, strain_tol, tau0, viscosity, x1, x2, x3, y1, y2, y3, z1, z2, and z3.

◆ isCharacteristicMtrxSymmetric()

bool oofem::TrabBone3D::isCharacteristicMtrxSymmetric ( MatResponseMode rMode) const
inlineoverridevirtual

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

Reimplemented from oofem::Material.

Definition at line 195 of file trabbone3d.h.

◆ performPlasticityReturn()

◆ predictRelativeComputationalCost()

double oofem::TrabBone3D::predictRelativeComputationalCost ( GaussPoint * gp)
overridevirtual

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.

Reimplemented from oofem::Material.

Definition at line 960 of file trabbone3d.C.

References oofem::Material::giveStatus().

◆ predictRelativeRedistributionCost()

double oofem::TrabBone3D::predictRelativeRedistributionCost ( GaussPoint * gp)
overridevirtual

Returns the relative redistribution cost of the receiver

Reimplemented from oofem::Material.

Definition at line 972 of file trabbone3d.C.

◆ projectOnYieldSurface()

Member Data Documentation

◆ chi0

double oofem::TrabBone3D::chi0 = 0.
protected

Definition at line 177 of file trabbone3d.h.

Referenced by constructAnisoFabricTensor(), and initializeFrom().

◆ chi0Neg

double oofem::TrabBone3D::chi0Neg = 0.
protected

Definition at line 177 of file trabbone3d.h.

Referenced by initializeFrom().

◆ chi0Pos

double oofem::TrabBone3D::chi0Pos = 0.
protected

Definition at line 177 of file trabbone3d.h.

Referenced by initializeFrom().

◆ critDam

◆ densCrit

◆ eps0

double oofem::TrabBone3D::eps0 = 0.
protected

◆ expDam

◆ expk

double oofem::TrabBone3D::expk = 0.
protected

◆ expl

double oofem::TrabBone3D::expl = 0.
protected

◆ expp

double oofem::TrabBone3D::expp = 0.
protected

Definition at line 179 of file trabbone3d.h.

Referenced by constructAnisoFabricTensor(), constructAnisoFtensor(), and initializeFrom().

◆ expPlasHard

double oofem::TrabBone3D::expPlasHard = 0.
protected

◆ expq

double oofem::TrabBone3D::expq = 0.
protected

Definition at line 179 of file trabbone3d.h.

Referenced by constructAnisoFabricTensor(), constructAnisoFtensor(), and initializeFrom().

◆ gammaL

double oofem::TrabBone3D::gammaL = 0.
protected

Definition at line 188 of file trabbone3d.h.

◆ gammaL0

double oofem::TrabBone3D::gammaL0 = 0.
protected

◆ gammaP

double oofem::TrabBone3D::gammaP = 0.
protected

Definition at line 188 of file trabbone3d.h.

◆ gammaP0

double oofem::TrabBone3D::gammaP0 = 0.
protected

◆ m1

double oofem::TrabBone3D::m1 = 0.
protected

◆ m2

double oofem::TrabBone3D::m2 = 0.
protected

◆ max_num_iter

int oofem::TrabBone3D::max_num_iter = 0
protected

Definition at line 181 of file trabbone3d.h.

Referenced by initializeFrom(), and projectOnYieldSurface().

◆ mu0

double oofem::TrabBone3D::mu0 = 0.
protected

◆ nu0

double oofem::TrabBone3D::nu0 = 0.
protected

◆ plasHardFactor

double oofem::TrabBone3D::plasHardFactor = 0.
protected

◆ printflag

int oofem::TrabBone3D::printflag = 0
protected

◆ rel_yield_tol

double oofem::TrabBone3D::rel_yield_tol = 0.
protected

Definition at line 182 of file trabbone3d.h.

Referenced by initializeFrom(), and projectOnYieldSurface().

◆ rho

◆ rL

double oofem::TrabBone3D::rL = 0.
protected

◆ rP

double oofem::TrabBone3D::rP = 0.
protected

◆ sig0Neg

double oofem::TrabBone3D::sig0Neg = 0.
protected

Definition at line 176 of file trabbone3d.h.

Referenced by constructAnisoFabricTensor(), constructAnisoFtensor(), and initializeFrom().

◆ sig0Pos

double oofem::TrabBone3D::sig0Pos = 0.
protected

Definition at line 176 of file trabbone3d.h.

Referenced by constructAnisoFabricTensor(), constructAnisoFtensor(), and initializeFrom().

◆ strain_tol

double oofem::TrabBone3D::strain_tol = 0.
protected

Definition at line 182 of file trabbone3d.h.

Referenced by initializeFrom(), and projectOnYieldSurface().

◆ tau0

double oofem::TrabBone3D::tau0 = 0.
protected

Definition at line 178 of file trabbone3d.h.

Referenced by constructAnisoFabricTensor(), and initializeFrom().

◆ tDens

double oofem::TrabBone3D::tDens = 0.
protected

◆ viscosity

double oofem::TrabBone3D::viscosity = 0.
protected

Viscosity parameter.

Definition at line 190 of file trabbone3d.h.

Referenced by evaluateCurrentViscousModulus(), evaluateCurrentViscousStress(), and initializeFrom().

◆ x1

double oofem::TrabBone3D::x1 = 1.
protected

Local coordinate system.

Definition at line 184 of file trabbone3d.h.

Referenced by constructFabricTransformationMatrix(), constructStiffnessTransformationMatrix(), and initializeFrom().

◆ x2

double oofem::TrabBone3D::x2 = 0.
protected

◆ x3

double oofem::TrabBone3D::x3 = 0.
protected

◆ y1

double oofem::TrabBone3D::y1 = 0.
protected

◆ y2

double oofem::TrabBone3D::y2 = 1.
protected

◆ y3

double oofem::TrabBone3D::y3 = 0.
protected

◆ z1

double oofem::TrabBone3D::z1 = 0.
protected

◆ z2

double oofem::TrabBone3D::z2 = 0.
protected

◆ z3

double oofem::TrabBone3D::z3 = 1.
protected

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

This page is part of the OOFEM-3.0 documentation. Copyright Copyright (C) 1994-2025 Borek Patzak Bořek Patzák
Project e-mail: oofem@fsv.cvut.cz
Generated at for OOFEM by doxygen 1.15.0 written by Dimitri van Heesch, © 1997-2011