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

#include <frcfcm.h>

Inheritance diagram for oofem::FRCFCM:
Collaboration diagram for oofem::FRCFCM:

Public Member Functions

 FRCFCM (int n, Domain *d)
void initializeFrom (InputRecord &ir) override
std::unique_ptr< MaterialStatusCreateStatus (GaussPoint *gp) const override
int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
Public Member Functions inherited from oofem::ConcreteFCM
 ConcreteFCM (int n, Domain *d)
const char * giveClassName () const override
const char * giveInputRecordName () const override
double give (int aProperty, GaussPoint *gp) const override
MaterialStatusgiveStatus (GaussPoint *gp) const override
Public Member Functions inherited from oofem::FCMMaterial
 FCMMaterial (int n, Domain *d)
bool hasMaterialModeCapability (MaterialMode mode) const override
const char * giveClassName () const override
double give (int aProperty, GaussPoint *gp) const override
FloatMatrixF< 6, 6 > give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 3, 3 > givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
FloatMatrixF< 4, 4 > givePlaneStrainStiffMtrx (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const override
void giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) const override
virtual void initializeCrack (GaussPoint *gp, TimeStep *tStep, FloatMatrix &base, int nCrack) const
FloatArrayF< 6 > giveRealStressVector_3d (const FloatArrayF< 6 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
FloatArrayF< 4 > giveRealStressVector_PlaneStrain (const FloatArrayF< 4 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_3d.
FloatArrayF< 3 > giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 1 > giveRealStressVector_1d (const FloatArrayF< 1 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 2 > giveRealStressVector_2dBeamLayer (const FloatArrayF< 2 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_StressControl.
FloatArrayF< 5 > giveRealStressVector_PlateLayer (const FloatArrayF< 5 > &strain, GaussPoint *gp, TimeStep *tStep) const override
 Default implementation relies on giveRealStressVector_StressControl.
virtual double computeNormalCrackOpening (GaussPoint *gp, int i) const
 uses temporary cracking strain and characteristic length to obtain the crack opening
virtual double computeMaxNormalCrackOpening (GaussPoint *gp, TimeStep *tStep, int i) const
 uses maximum equilibrated cracking strain and characteristic length to obtain the maximum reached crack opening
virtual double computeShearSlipOnCrack (GaussPoint *gp, TimeStep *tStep, int i) const
 computes total shear slip on a given crack plane (i = 1, 2, 3); the slip is computed from the temporary cracking strain
Public Member Functions inherited from oofem::StructuralMaterial
 StructuralMaterial (int n, Domain *d)
bool hasMaterialModeCapability (MaterialMode mode) const override
const char * giveClassName () const override
void initializeFrom (InputRecord &ir) override
void giveInputRecord (DynamicInputRecord &input) override
void giveCharacteristicMatrix (FloatMatrix &answer, MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const override
 Returns characteristic matrix of the receiver.
void giveCharacteristicVector (FloatArray &answer, FloatArray &flux, MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const override
 Returns characteristic vector of the receiver.
virtual void giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArray giveRealStressVector_StressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint *gp, TimeStep *tStep) const
 Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArray giveRealStressVector_ShellStressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArrayF< 2 > giveRealStressVector_Warping (const FloatArrayF< 2 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 > giveRealStressVector_Fiber (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 > giveRealStressVector_2dPlateSubSoil (const FloatArrayF< 3 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
 Default implementation is not provided.
virtual FloatArrayF< 6 > giveRealStressVector_3dBeamSubSoil (const FloatArrayF< 6 > &reducedE, GaussPoint *gp, TimeStep *tStep) const
virtual FloatArrayF< 9 > giveFirstPKStressVector_3d (const FloatArrayF< 9 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
virtual FloatArrayF< 5 > giveFirstPKStressVector_PlaneStrain (const FloatArrayF< 5 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveFirstPKStressVector_3d.
virtual FloatArray giveFirstPKStressVector_StressControl (const FloatArray &reducedvF, const IntArray &FControl, GaussPoint *gp, TimeStep *tStep) const
 Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArrayF< 4 > giveFirstPKStressVector_PlaneStress (const FloatArrayF< 4 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveFirstPKStressVector_StressControl.
virtual FloatArrayF< 1 > giveFirstPKStressVector_1d (const FloatArrayF< 1 > &vF, GaussPoint *gp, TimeStep *tStep) const
 Default implementation relies on giveFirstPKStressVector_StressControl.
virtual void giveCauchyStressVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void giveCauchyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void giveCauchyStressVector_PlaneStress (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void giveCauchyStressVector_1d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual void giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep)
virtual FloatArrayF< 6 > giveThermalDilatationVector (GaussPoint *gp, TimeStep *tStep) const
double giveReferenceTemperature ()
virtual FloatArray computeStressIndependentStrainVector (GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const
FloatArrayF< 6 > computeStressIndependentStrainVector_3d (GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const
virtual void giveStiffnessMatrix_dPdF (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep)
virtual FloatMatrixF< 9, 9 > give3dMaterialStiffnessMatrix_dPdF (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const
virtual void give3dMaterialStiffnessMatrix_dCde (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep)
void giveStressDependentPartOfStrainVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrainVector, TimeStep *tStep, ValueModeType mode) const
int setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type) override
int giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override
virtual FloatMatrixF< 4, 4 > givePlaneStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const
virtual void givePlaneStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep)
virtual FloatMatrixF< 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 bool isCharacteristicMtrxSymmetric (MatResponseMode rMode) const
virtual double giveCharacteristicValue (MatResponseMode type, GaussPoint *gp, TimeStep *tStep) const
 Returns characteristic value of the receiver.
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 int packUnknowns (DataStream &buff, TimeStep *tStep, GaussPoint *ip)
virtual int unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep, GaussPoint *ip)
virtual int estimatePackSize (DataStream &buff, GaussPoint *ip)
virtual double predictRelativeComputationalCost (GaussPoint *gp)
virtual double predictRelativeRedistributionCost (GaussPoint *gp)
virtual void initTempStatus (GaussPoint *gp) const
void saveContext (DataStream &stream, ContextMode mode) override
void restoreContext (DataStream &stream, ContextMode mode) override
Public Member Functions inherited from oofem::FEMComponent
 FEMComponent (int n, Domain *d)
virtual ~FEMComponent ()=default
 Virtual destructor.
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).
Public Member Functions inherited from oofem::RandomMaterialExtensionInterface
 RandomMaterialExtensionInterface ()
 Constructor.
virtual ~RandomMaterialExtensionInterface ()
 Destructor.
void initializeFrom (InputRecord &ir)
void giveInputRecord (DynamicInputRecord &ir)
bool give (int key, GaussPoint *gp, double &value) const
Public Member Functions inherited from oofem::Interface
 Interface ()
 Constructor.
virtual ~Interface ()

Protected Types

enum  FiberShearStrengthType {
  FSS_NONE , FSS_Sajdlova , FSS_Kabele , FSS_Havlasek ,
  FSS_Unknown
}
enum  FiberDamageType { FDAM_NONE , FDAM_GammaCrackLin , FDAM_GammaCrackExp , FDAM_Unknown }
 Type of fibre damage which is triggered by crack shearing strain = w / u. More...
enum  FiberType {
  FT_CAF , FT_SAF , FT_SRF , FT_SRF2D ,
  FT_Unknown
}
Protected Types inherited from oofem::ConcreteFCM
enum  SofteningType {
  ST_NONE , ST_Exponential , ST_Linear , ST_Hordijk ,
  ST_UserDefinedCrack , ST_LinearHardeningStrain , ST_UserDefinedStrain , ST_Unknown
}
 type of post-peak behavior in the normal direction to the crack plane More...
enum  ShearRetentionType {
  SHR_NONE , SHR_Const_ShearRetFactor , SHR_Const_ShearFactorCoeff , SHR_UserDefined_ShearRetFactor ,
  SHR_Unknown
}
 type of reduction of the shear stiffness caused by cracking More...
enum  ShearStrengthType {
  SHS_NONE , SHS_Const_Ft , SHS_Collins_Interlock , SHS_Residual_Ft ,
  SHS_Unknown
}
 defines the maximum value of shear stress More...

Protected Member Functions

double giveCrackingModulus (MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep, int i) const override
 returns stiffness in the normal direction of the i-th crack
virtual double computeFiberBond (double w) const
 evaluates the fiber bond if w > w*
double giveNormalCrackingStress (GaussPoint *gp, TimeStep *tStep, double eps_cr, int i) const override
 computes normal stress associated with i-th crack direction
virtual double computeStressInFibersInCracked (GaussPoint *gp, TimeStep *tStep, double eps_cr, int i) const
 compute the nominal stress in fibers in the i-th crack
double computeEffectiveShearModulus (GaussPoint *gp, TimeStep *tStep, int i) const override
 returns Geff which is necessary in the global stiffness matrix
double computeD2ModulusForCrack (GaussPoint *gp, TimeStep *tStep, int icrack) const override
 shear modulus for a given crack plane (1, 2, 3)
virtual double estimateD2ModulusForCrack (GaussPoint *gp, TimeStep *tStep, int icrack) const
 estimate shear modulus for a given crack plane (1, 2, 3). Uses equilibrated value of damage.
double maxShearStress (GaussPoint *gp, TimeStep *tStep, int i) const override
 computes the maximum value of the shear stress; if the shear stress exceeds this value, it is cropped
virtual double computeTempDamage (GaussPoint *gp, TimeStep *tStep) const
 evaluates temporary value of damage caused by fibre shearing
virtual double computeCrackSpacing ()
 computes crack spacing based on composition of the fibre composite
virtual double computeCrackFibreAngle (GaussPoint *gp, int i) const
 compute the angle between the fibre and i-th crack normal
void checkSnapBack (GaussPoint *gp, TimeStep *tStep, int crack) const override
 checks possible snap-back
bool isStrengthExceeded (const FloatMatrix &base, GaussPoint *gp, TimeStep *tStep, int iCrack, double trialStress) const override
 compares trial stress with strength. Returns true if the strength is exceeded. Function oveloaded in the nonlocal approach for the fiber reinforced composites
double computeShearStiffnessRedistributionFactor (GaussPoint *gp, TimeStep *tStep, int ithCrackPlane, int jthCrackDirection) const override
 function calculating ratio used to split shear slips on two crack planes
double computeOverallElasticStiffness (GaussPoint *gp, TimeStep *tStep) const override
 returns overall Young's modulus
double computeOverallElasticShearModulus (GaussPoint *gp, TimeStep *tStep) const override
 returns overall shear modulus
Protected Member Functions inherited from oofem::ConcreteFCM
double giveTensileStrength (GaussPoint *gp, TimeStep *tStep) const override
 comutes tensile strength
virtual double giveFractureEnergy (GaussPoint *gp, TimeStep *tStep) const
double giveCrackingModulusInTension (MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep, int i) const override
 returns stiffness in the normal direction of the i-th crack
double computeNumerD2ModulusForCrack (GaussPoint *gp, TimeStep *tStep, int icrack) const override
 shear modulus for numerical purpose (stiffness matrix) for a given crack plane (1, 2, 3)
virtual double computeResidualTensileStrength (GaussPoint *gp, TimeStep *tStep) const
 based on the maximum crack opening evaluates the residual strength
Protected Member Functions inherited from oofem::FCMMaterial
virtual void updateCrackStatus (GaussPoint *gp) const
 updates crack statuses
virtual double giveCharacteristicElementLength (GaussPoint *gp, const FloatArray &crackPlaneNormal) const
 returns characteristic element length in given direction
virtual double computeTotalD2Modulus (GaussPoint *gp, TimeStep *tStep, int i) const
 shear modulus for a given shear direction (4, 5, 6)
virtual double computeNumerD2Modulus (GaussPoint *gp, TimeStep *tStep, int i) const
 shear modulus in a STIFFNESS MATRIX for a given shear direction (4, 5, 6)
virtual bool isIntactForShear (GaussPoint *gp, int i) const
 returns true for closed or no cracks associated to given shear direction (i = 4, 5, 6)
virtual bool isIntact (GaussPoint *gp, int icrack) const
 returns true for closed or no crack (i = 1, 2, 3)
virtual bool isThisShearComponent (GaussPoint *gp, int component) const
 returns true if current component is associated with shear
virtual bool checkStrengthCriterion (FloatMatrix &newBase, const FloatArray &globalStress, GaussPoint *gp, TimeStep *tStep, int nCrack) const
 checks if the globalStress does not exceed strength in the direction of newBase for n-th crack
virtual double giveCrackSpacing (void) const
 returns either user-provided value of crack spacing or a value computed from composition
virtual double giveNumberOfCracksInDirection (GaussPoint *gp, int iCrack) const
 returns number of fictiotious parallel cracks in the direction of i-th crack
virtual double giveNumberOfCracksForShearDirection (GaussPoint *gp, int i) const
 returns number of cracks for given shear direction (i = 4, 5, 6) which is treated as the maximum of the two associated normal directions
virtual void giveMaterialStiffnessMatrix (FloatMatrix &answer, MatResponseMode, GaussPoint *gp, TimeStep *tStep) const
virtual void giveTotalLocalCrackedStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep) const
 returns local stiffness matrix in the cracks' direction (total according to the material mode)
virtual void giveNormalLocalCrackedStiffnessMatrix (FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep) const
 returns local stiffness matrix in the cracks' direction (only normal components)
virtual double givePoissonsRatio () const
 returns Poisson's ratio
Protected Member Functions inherited from oofem::RandomMaterialExtensionInterface
void _generateStatusVariables (GaussPoint *) const

Protected Attributes

double tau_0 = 0.
 fiber shear strength at zero slip
double b0 = 0.
 micromechanical parameter for fiber shear according to Sajdlova
double b1 = 0.
 micromechanical parameter for fiber shear according to Kabele
double b2 = 0.
double b3 = 0.
double f = 0.
 snubbing factor "f"
double g = 0.
 auxiliary parameter computed from snubbing factor "f"
double Vf = 0.
 volume fraction of fibers
double Lf = 0.
 fiber length
double Df = 0.
 fiber diameter
double Ef = 0.
 fiber Young's modulus
double Gfib = 0.
 fiber shear modulus
double kfib = 0.
 fiber cross-sectional shear factor
double w_star = 0.
 transitional opening
double eta = 0.
 aux. factor
double gammaCrackFail = 0.
 shear strain at full fibers rupture
double minDamageOpening = 0.
 minimum opening at which damage can start
int M = 0
FloatArray orientationVector
 orientation of fibres
double fibreActivationOpening = 0.
 crack opening at which the crossing fibers begin to be activated
double dw0 = 0.
double dw1 = 0.
bool smoothen = false
FiberShearStrengthType fiberShearStrengthType = FSS_Unknown
FiberDamageType fiberDamageType = FDAM_Unknown
FiberType fiberType = FT_Unknown
Protected Attributes inherited from oofem::ConcreteFCM
double Gf = 0.
 Fracture energy.
double Ft = 0.
 Tensile strength.
double beta = 0.
 shear retention factor
double sf = 0.
 shear factor
double sf_numer = 0.
 shear factor for numerical purpose
double fc = 0.
 Collins' aggregate interlock: compressive strength in MPa.
double ag = 0.
 Collins' aggregate interlock: aggregate diameter in appropriate units (same as FE mesh).
double lengthScale = 0.
 Collins' aggregate interlock: 1 for meter, 1000 for analysis in mm.
FloatArray soft_w
 user-defined softening (traction-COD)
FloatArray soft_function_w
FloatArray soft_eps
 user-defined softening (traction-strain)
FloatArray soft_function_eps
FloatArray beta_w
 user-defined shear retention factor (with respect to crack opening)
FloatArray beta_function
double H = 0.
 hardening modulus
double eps_f = 0.
 strain at failure
SofteningType softType = ST_NONE
ShearRetentionType shearType = SHR_Unknown
ShearStrengthType shearStrengthType = SHS_Unknown
Protected Attributes inherited from oofem::FCMMaterial
IsotropicLinearElasticMaterial linearElasticMaterial
int nAllowedCracks
 allowed number of cracks (user-defined)
ElementCharSizeMethod ecsMethod
 Method used for evaluation of characteristic element size.
bool multipleCrackShear
 if true = takes shear compliance of all cracks, false = only dominant crack contribution, default value is false
int iterLimitGlobal = 20
double crackSpacing
 value of crack spacing (allows to "have" more parallel cracks in one direction if the element size exceeds user-defined or computed crack spacing).
double shearCoeffNumer
 minimum ratio of effective shear modulus vs. shear modulus - just for numerical purpose
double normalCoeffNumer
 minimum ratio of effective normal stiffness vs. overall modulus - just for numerical purpose
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.
Protected Attributes inherited from oofem::RandomMaterialExtensionInterface
IntArray randVariables
 Array of randomized variables (identified by a key).
IntArray randomVariableGenerators
 Array of generators id's for corresponding randomized variables.

Additional Inherited Members

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

Detailed Description

This class implements a FRCFCM material (Fiber Reinforced Concrete base on Fixed Crack Model) in a finite element problem. This class provides an extension to the ConcreteFCM which serves as a material model for matrix while the present class FRCFCM adds the contribution of fibers. The contribution of the two constituents are defined by Vf parameter = volume of fibers.

Definition at line 110 of file frcfcm.h.

Member Enumeration Documentation

◆ FiberDamageType

Type of fibre damage which is triggered by crack shearing strain = w / u.

Enumerator
FDAM_NONE 
FDAM_GammaCrackLin 
FDAM_GammaCrackExp 
FDAM_Unknown 

Definition at line 197 of file frcfcm.h.

◆ FiberShearStrengthType

Type strength of the shear bond. This is activated only for short fibers (SAF or SRF) once the maximum crack opening exceeds w*

Enumerator
FSS_NONE 
FSS_Sajdlova 
FSS_Kabele 
FSS_Havlasek 
FSS_Unknown 

Definition at line 193 of file frcfcm.h.

◆ FiberType

enum oofem::FRCFCM::FiberType
protected

Type fo fibers in the composite. CAF = continuous aligned fibers SAF = short aligned fibers SRF = short random fibers SRF2D = short fibers random in 2D

Enumerator
FT_CAF 
FT_SAF 
FT_SRF 
FT_SRF2D 
FT_Unknown 

Definition at line 207 of file frcfcm.h.

Constructor & Destructor Documentation

◆ FRCFCM()

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

Definition at line 45 of file frcfcm.C.

References oofem::ConcreteFCM::ConcreteFCM().

Referenced by oofem::FRCFCMNL::FRCFCMNL().

Member Function Documentation

◆ checkSnapBack()

void oofem::FRCFCM::checkSnapBack ( GaussPoint * gp,
TimeStep * tStep,
int crack ) const
inlineoverrideprotectedvirtual

checks possible snap-back

Reimplemented from oofem::ConcreteFCM.

Definition at line 227 of file frcfcm.h.

◆ computeCrackFibreAngle()

◆ computeCrackSpacing()

double oofem::FRCFCM::computeCrackSpacing ( )
protectedvirtual

computes crack spacing based on composition of the fibre composite

Definition at line 1093 of file frcfcm.C.

References Df, f, fiberType, oofem::ConcreteFCM::Ft, FT_CAF, FT_SAF, FT_SRF, FT_SRF2D, g, Lf, M_PI, OOFEM_ERROR, tau_0, and Vf.

Referenced by initializeFrom().

◆ computeD2ModulusForCrack()

double oofem::FRCFCM::computeD2ModulusForCrack ( GaussPoint * gp,
TimeStep * tStep,
int icrack ) const
overrideprotectedvirtual

◆ computeEffectiveShearModulus()

double oofem::FRCFCM::computeEffectiveShearModulus ( GaussPoint * gp,
TimeStep * tStep,
int i ) const
overrideprotectedvirtual

◆ computeFiberBond()

double oofem::FRCFCM::computeFiberBond ( double w) const
protectedvirtual

◆ computeOverallElasticShearModulus()

double oofem::FRCFCM::computeOverallElasticShearModulus ( GaussPoint * gp,
TimeStep * tStep ) const
inlineoverrideprotectedvirtual

returns overall shear modulus

Reimplemented from oofem::FCMMaterial.

Definition at line 231 of file frcfcm.h.

References computeOverallElasticStiffness(), and oofem::FCMMaterial::linearElasticMaterial.

Referenced by computeEffectiveShearModulus().

◆ computeOverallElasticStiffness()

double oofem::FRCFCM::computeOverallElasticStiffness ( GaussPoint * gp,
TimeStep * tStep ) const
overrideprotectedvirtual

returns overall Young's modulus

Reimplemented from oofem::FCMMaterial.

Definition at line 1183 of file frcfcm.C.

References Ef, fiberType, FT_CAF, FT_SAF, FT_SRF, FT_SRF2D, oofem::FCMMaterial::linearElasticMaterial, OOFEM_ERROR, and Vf.

Referenced by computeOverallElasticShearModulus().

◆ computeShearStiffnessRedistributionFactor()

double oofem::FRCFCM::computeShearStiffnessRedistributionFactor ( GaussPoint * gp,
TimeStep * tStep,
int ithCrackPlane,
int jthCrackDirection ) const
overrideprotectedvirtual

function calculating ratio used to split shear slips on two crack planes

Reimplemented from oofem::FCMMaterial.

Definition at line 1144 of file frcfcm.C.

References estimateD2ModulusForCrack().

◆ computeStressInFibersInCracked()

◆ computeTempDamage()

◆ CreateStatus()

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

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

Reimplemented in oofem::FRCFCMNL.

Definition at line 117 of file frcfcm.h.

◆ estimateD2ModulusForCrack()

double oofem::FRCFCM::estimateD2ModulusForCrack ( GaussPoint * gp,
TimeStep * tStep,
int icrack ) const
protectedvirtual

◆ giveCrackingModulus()

◆ giveIPValue()

int oofem::FRCFCM::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::ConcreteFCM.

Reimplemented in oofem::FRCFCMNL.

Definition at line 1159 of file frcfcm.C.

References oofem::FloatArray::at(), computeStressInFibersInCracked(), oofem::FCMMaterialStatus::giveCrackStrain(), oofem::FRCFCMStatus::giveDamage(), oofem::ConcreteFCM::giveStatus(), and oofem::FloatArray::resize().

◆ giveNormalCrackingStress()

double oofem::FRCFCM::giveNormalCrackingStress ( GaussPoint * gp,
TimeStep * tStep,
double eps_cr,
int i ) const
overrideprotectedvirtual

computes normal stress associated with i-th crack direction

Reimplemented from oofem::ConcreteFCM.

Definition at line 602 of file frcfcm.C.

References computeStressInFibersInCracked(), and Vf.

◆ initializeFrom()

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

Reimplemented from oofem::ConcreteFCM.

Reimplemented in oofem::FRCFCMNL.

Definition at line 49 of file frcfcm.C.

References _IFT_FRCFCM_b0, _IFT_FRCFCM_b1, _IFT_FRCFCM_b2, _IFT_FRCFCM_b3, _IFT_FRCFCM_computeCrackSpacing, _IFT_FRCFCM_Df, _IFT_FRCFCM_dw0, _IFT_FRCFCM_dw1, _IFT_FRCFCM_Ef, _IFT_FRCFCM_f, _IFT_FRCFCM_fDamType, _IFT_FRCFCM_fiberType, _IFT_FRCFCM_fibreActivationOpening, _IFT_FRCFCM_fssType, _IFT_FRCFCM_gammaCrack, _IFT_FRCFCM_Gfib, _IFT_FRCFCM_kfib, _IFT_FRCFCM_Lf, _IFT_FRCFCM_M, _IFT_FRCFCM_nuf, _IFT_FRCFCM_orientationVector, _IFT_FRCFCM_tau_0, _IFT_FRCFCM_Vf, _IFT_IsotropicLinearElasticMaterial_e, b0, b1, b2, b3, computeCrackSpacing(), oofem::FCMMaterial::crackSpacing, Df, dw0, dw1, Ef, eta, f, FDAM_GammaCrackExp, FDAM_GammaCrackLin, FDAM_NONE, FDAM_Unknown, fiberDamageType, fiberShearStrengthType, fiberType, fibreActivationOpening, FSS_Havlasek, FSS_Kabele, FSS_NONE, FSS_Sajdlova, FSS_Unknown, FT_CAF, FT_SAF, FT_SRF, FT_SRF2D, FT_Unknown, g, gammaCrackFail, Gfib, oofem::InputRecord::hasField(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, kfib, length(), Lf, M, M_PI, OOFEM_ERROR, orientationVector, smoothen, tau_0, Vf, and w_star.

◆ isStrengthExceeded()

bool oofem::FRCFCM::isStrengthExceeded ( const FloatMatrix & base,
GaussPoint * gp,
TimeStep * tStep,
int iCrack,
double trialStress ) const
overrideprotectedvirtual

compares trial stress with strength. Returns true if the strength is exceeded. Function oveloaded in the nonlocal approach for the fiber reinforced composites

Reimplemented from oofem::FCMMaterial.

Reimplemented in oofem::FRCFCMNL.

Definition at line 1117 of file frcfcm.C.

References Ef, oofem::ConcreteFCM::giveTensileStrength(), oofem::FCMMaterial::linearElasticMaterial, and Vf.

◆ maxShearStress()

double oofem::FRCFCM::maxShearStress ( GaussPoint * gp,
TimeStep * tStep,
int i ) const
overrideprotectedvirtual

Member Data Documentation

◆ b0

double oofem::FRCFCM::b0 = 0.
protected

micromechanical parameter for fiber shear according to Sajdlova

Definition at line 126 of file frcfcm.h.

Referenced by computeFiberBond(), and initializeFrom().

◆ b1

double oofem::FRCFCM::b1 = 0.
protected

micromechanical parameter for fiber shear according to Kabele

Definition at line 128 of file frcfcm.h.

Referenced by computeFiberBond(), and initializeFrom().

◆ b2

double oofem::FRCFCM::b2 = 0.
protected

Definition at line 128 of file frcfcm.h.

Referenced by computeFiberBond(), and initializeFrom().

◆ b3

double oofem::FRCFCM::b3 = 0.
protected

Definition at line 128 of file frcfcm.h.

Referenced by computeFiberBond(), and initializeFrom().

◆ Df

◆ dw0

double oofem::FRCFCM::dw0 = 0.
protected

smooth transition of the bridging stress if fibreActivationOpening is applied dw0 = distance from the fibreActivationOpening where the smooth transition starts dw1 = distance from the fibreActivationOpening where the smooth transition ends smoothen = flag

Definition at line 186 of file frcfcm.h.

Referenced by computeStressInFibersInCracked(), giveCrackingModulus(), and initializeFrom().

◆ dw1

double oofem::FRCFCM::dw1 = 0.
protected

Definition at line 186 of file frcfcm.h.

Referenced by computeStressInFibersInCracked(), giveCrackingModulus(), and initializeFrom().

◆ Ef

◆ eta

double oofem::FRCFCM::eta = 0.
protected

◆ f

◆ fiberDamageType

FiberDamageType oofem::FRCFCM::fiberDamageType = FDAM_Unknown
protected

Definition at line 198 of file frcfcm.h.

Referenced by computeTempDamage(), and initializeFrom().

◆ fiberShearStrengthType

FiberShearStrengthType oofem::FRCFCM::fiberShearStrengthType = FSS_Unknown
protected

Definition at line 194 of file frcfcm.h.

Referenced by computeFiberBond(), and initializeFrom().

◆ fiberType

◆ fibreActivationOpening

double oofem::FRCFCM::fibreActivationOpening = 0.
protected

◆ g

double oofem::FRCFCM::g = 0.
protected

◆ gammaCrackFail

double oofem::FRCFCM::gammaCrackFail = 0.
protected

shear strain at full fibers rupture

Definition at line 161 of file frcfcm.h.

Referenced by computeTempDamage(), and initializeFrom().

◆ Gfib

double oofem::FRCFCM::Gfib = 0.
protected

fiber shear modulus

Definition at line 149 of file frcfcm.h.

Referenced by computeD2ModulusForCrack(), estimateD2ModulusForCrack(), initializeFrom(), and maxShearStress().

◆ kfib

double oofem::FRCFCM::kfib = 0.
protected

fiber cross-sectional shear factor

Definition at line 152 of file frcfcm.h.

Referenced by computeD2ModulusForCrack(), estimateD2ModulusForCrack(), initializeFrom(), and maxShearStress().

◆ Lf

◆ M

int oofem::FRCFCM::M = 0
protected

Exponent in the unloading-reloading constitutive law. the function is defined by two points - the origin and the stress at maximum cracking strain sigma(w) = sig_max * w^M / w_max^M

Definition at line 172 of file frcfcm.h.

Referenced by computeStressInFibersInCracked(), giveCrackingModulus(), and initializeFrom().

◆ minDamageOpening

double oofem::FRCFCM::minDamageOpening = 0.
protected

minimum opening at which damage can start

Definition at line 164 of file frcfcm.h.

◆ orientationVector

FloatArray oofem::FRCFCM::orientationVector
protected

◆ smoothen

bool oofem::FRCFCM::smoothen = false
protected

Definition at line 187 of file frcfcm.h.

Referenced by computeStressInFibersInCracked(), giveCrackingModulus(), and initializeFrom().

◆ tau_0

double oofem::FRCFCM::tau_0 = 0.
protected

◆ Vf

◆ w_star

double oofem::FRCFCM::w_star = 0.
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