This material model is an extension of the ConcreteFCM described in Section 1.6.11.
It is possible to choose from three different “classes” of fibers. This choice is controlled by keyword fiberType: 0 = continuous aligned fibers (CAF), 1 = short aligned fibers (SAF) and 2 = short random fibers (SRF). Currently, it is not possible to combine more classes of fibers in one material model.
All of the above-mentioned fiber types are further defined by the material properties and geometry. Fiber quantity is captured by the dimensionless volume fraction Vf (as decimal). All fibers are assumed to have a circular cross-section (shape factor for shear kfib with default value 0.9) and to possess the same geometry characterized by the diameter Df and length Lf.
The overall elastic stiffness of the fiber-reinforced composite is calculated as a weighted average of the moduli of matrix and fibers and the Poisson's ratio is considered to be equal to the Poisson's ratio of the matrix. Similarly to ConcreteFCM, cracking is initiated once the tensile stress in matrix reaches the tensile strength .
The nominal bridging stress for continuous aligned fibers which are not perpendicular to the crack plane can be very easily obtained by multiplying the nominal bridging stress for perpendicular fibers by two terms: the first one reflecting lower volume of inclined fibers passing through the crack plane and the second one capturing the snubbing effect:
For the CAF perpendicular to the crack, the nominal bridging stress can be derived as
For SRF the nominal bridging stress is
Larger pull-out displacements can lead to significant physical changes in the fiber surface which can result into changes in the bond shear stress. This phenomenon is captured by function relating the frictional bond to the crack opening and is implemented in three alternative formulations. (In order to keep use fssType = 0.) In conventional FRC with ordinary concrete matrix, the frictional bond usually decreases with increasing slip. To capture this type of behavior we adopt the function proposed by Sajdlová (activated with fssType = 1) reads
The bridging stress for SRF is defined as
It is possible to delay the activation of the stress in fibers using prameteter fibreActivationOpening which can be imagined as a “lag” of the fiber-related crack opening behind the matrix-related crack opening.
During unloading the stress in fibers does not decrease linearly to origin. Current implementation uses a power function
The influence of crack opening and sliding on the bridging shear stress only due to fibers is expressed as
It has been found that in some high performance fiber reinforced cement composites, fibers rupture when cracks are exposed to shearing. This phenomenon is modeled by damage parameter , which accounts for the ratio of ruptured fibers and varies between the values of 0 and 1. It is assumed that depends on the maximum shear strain sustained by the protruding portions of bridging fibers throughout the loading history. This crack shear strain can be expressed as:
where is the crack sliding displacement (CSD) and and is the maximum value of the crack opening displacement of the i-th crack. This means that the damage does not grow if the crack closes (crack opening decreases). If more cracks exist, the maximum contribution is considered.
Two different one-parameter damage evolution laws are currently implemented. For fDamType = 0 the damage is deactivated, with fDamType = 1 damage is described by
Since damage reduces the number of crack-bridging fibers, which is proportional to the fiber volume fraction, its effect can be suitably implemented by introducing the effective volume fraction
The material parameters are summarized in Tables 42 (matrix) and 43 (fiber extension).
Sample syntax for a fixed crack model reinforced with fibers with volume density 24 kN/m, thermal dilation coefficient
K, Young's modulus of the matrix 20 GPa, Poisson's ratio of matrix 0.2, fracture energy of matrix 100 N/m, tensile strength of matrix 2 MPa, linear tension softening, constant shear retention factor
, unlimited shear strength (shearStrengthType = 0), continuous aligned fibers, fiber volume 2%, fiber diameter 0.04 mm, Young's modulus of fibers 20 GPa, shear modulus of fibers 1 GPa, fiber-matrix bond strength 1 MPa, snubbing coefficient 0.7, shear correction coefficient 0.9, deactivated fiber damage, fiber act if COD exceeds 10 m (with smoothing from
m), fiber orientation at 45 degrees in x-y plane, automatic evaluation of crack spacing from composition; the analysis uses [m], [MPa] and [MN]:
FRCFCM 1 d 24.e-3 talpha 12.e-6 E 20000. n 0.2 Gf 100e-6 ft 2.0
softType 2 shearType 1 beta 0.05 FiberType 0 Vf 0.02 Df 0.04e-3
Ef 20000. Gfib 1000. tau_0 1. FSStype 0 f 0.7 kfib 0.9 fDamType 0
fibreactivationopening 10.e-6 dw0 2.e-6 dw1 1.e-6 orientationVector 3 1. 1. 0. computeCrackSpacing
|Description||Fixed crack model for FRC|
|Record Format||FRCFCM input record of ConcreteFCM Vf(rn) # Lf(rn) # Df(rn) # Ef(rn) # [ nuf(rn) #] [ Gfib(rn) #] [ kfib(rn) #] tau_0(rn) # b0(rn) # b1(rn) # b2(rn) # b3(rn) # f(rn) # [ M(in) #] [ fibreOrientationVector(ra) #] [ fssType(in) #] [ fDamType(in) #] [ fiberType(in) #] [ gammaCrack(rn) #] [ computeCrackSpacing ] [ fibreActivationOpening(rn) #] [ dw0(rn) #] [ dw1(rn) #]|
|Parameters||- Vf fiber content expressed as decimal|
|- Lf fiber length|
|- Df fiber diameter|
|- Ef fiber Young's modulus|
|- nuf fiber Poisson's ratio|
|- Gfib fiber shear modulus (read when nuf is not provided)|
|- kfib fiber cross-sectional shape correction factor|
|- tau_0 bond shear strength at zero slip|
|- b0 micromechanical parameter for fiber shear according to Sajdlová|
|- b1, b2, b3 micromechanical parameter for fiber shear according to Kabele|
|- f snubbing friction coefficient|
|- M exponent related to fiber unloading|
|- fibreOrientationVector vector specifying orientation for CAF and SAF fibers|
|- fssType type of Fiber bond Shear Strength (bond shear strength vs. crack opening)
|- fDamType type of damage law for fibers
|- fiberType type of reinforcing fibers
|- gammaCrack crack shear strain parameter applicable with fDamType = 1 or 2 (here the crack shear strain is understood as the crack slip divided by the crack opening )|
|- computeCrackSpacing crack spacing is evaluated automatically based on provided composition|
|- fibreActivationOpening crack opening at which the fibers begin transferring bridging stress|
|- dw0, dw1 applicable only if fibreActivationOpening , then it allows to smooth the traction-separation law for fibers; lower bound is fibreActivationOpening - dw0 and the upper bound is fibreActivationOpening - dw1|
|Supported modes||3dMat, PlaneStress, PlaneStrain|