Benefit of Fibres Addition to Normal Strength Concrete

19668 words (79 pages) Dissertation

13th Dec 2019 Dissertation Reference this

Tags: Environmental StudiesBuilding Materials

Disclaimer: This work has been submitted by a student. This is not an example of the work produced by our Dissertation Writing Service. You can view samples of our professional work here.

Any opinions, findings, conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of

Benefit of Fibres Addition to Normal Strength Concrete: An Empirical Model for the determination of Tensile Capacity and Fracture Energy


Fibre-reinforced concrete is widely used in construction to enhance the durability of concrete by controlling the formation and development of the material cracking. Its ability can be attributed to the ductility induced by the addition of fibres to the material matrix. The quantification of this improvement is investigated in this paper through the analysis of laboratory, computer and theoretical results. Testing of specifically designed dog bone shape specimens have shown that a 1.0% volume fibre addition allows for a 26% increase in tensile capacity under uniaxial stress conditions. While a general relation describing the induced improvement was defined as a function of the material characteristic cylinder strength, a general formulation involving the volume fraction of fibre was not possible due to the other fibre admixture being defective. A graphical resolution was, however, proposed for the available data.

The material behaviour was successfully modelled up to the maximum tensile stress as an adjustment of Hognestad’s model (1952) for concrete in compression. This was checked against the finite element analysis results showing a good fit limited to low stress values ( < 30% maximum stress) as the model appeared non-convergent outside this range.

Eventually, from the graphical analysis of the experimental results of the fibre admixtures, the material behaviour following the first cracking development was described. A model, based on the bi-linear method proposed by Hillerborg (1976), was adopted to estimate the total fracture energy dissipated by the cracking development. This outlined a 10-fold increase in the energy necessary for the initial cracking formation and a total fracture energy two orders of magnitude greater than the corresponding plain mixture.




List of tables

List of figures


1 / Introduction

1.1 / Aim and objectives

1.2 / Methodology

2 / Literature review

2.1 / Specimen profile design

2.2 / Concrete mix design

3 / Background

3.1 / Stress-strain curve

3.1.1 / Hognestad model

3.1.2 / Eurocode model

3.2 / Fracture energy

4 / Finite element modelling process

4.1 / Introduction

4.2 / Finite element analysis by ANSYS

4.2.1 / Concrete model

4.2.2 / Fibre model

4.3 / Material model

4.3.1 / Fibre distribution approximation

4.3.2 / Reduced model

4.3.3 / Computer model generation

4.3.4 / Reduced model results

4.3.5 / Validity of finite element model

4.4 / Tensile specimen model

4.4.1 / Model 1 – SOLID65 with reinforcement values

4.4.2 / Model 2 – SOLID65 with modelled fibres

4.4.3 / Model limitations

5 / Experimental process

5.1 / Mould design

5.1.1 / Design process

5.2 / Materials and mix details

5.3 / Test specimen and procedure

5.3.1 / Control cube

5.3.2 / Control cylinder

5.3.3 / Tensile specimen

5.4 / Results. Interpretation and reduction

5.4.1 / Readings up to failure

5.4.2 / Reading after failure

6 / Discussion

6.1 / Tensile capacity. Direct and splitting tests

6.2 / Modified Hognestad model

6.3 / Finite element model results

6.4 / Linearity of strain variation

6.5 / Formulation of tensile strength

6.6 / Fracture energy for fibre reinforced concrete

7 / Conclusion

7.1 / Recommendations for further studies



List of tables

Table 3.1 – Applied model values

Table 4.1 – Concrete properties

Table 4.2 – Fibre properties – PROPEX Enduro HPP45

Table 4.3  – Derived fibre’s properties

Table 4.4 – Real constant parameters

Table 5.1 – Results of moisture content test

Table 5.2 – Concrete mixes details – water cement ratio 0.5

Table 5.3 – Cube testing results

Table 5.4 – Cylinder testing results

Table 6.1 – Direct tensile test comparison

Table 6.2 – Results variation from Eurocodes’ values

Table 6.3 – Fitting lines’ characteristics

Table 6.4 – Characteristic strength interpolation

Table 6.5 – Fitting lines’ constant terms

Table 6.6 – Generalised linear equations

Table 6.7 – Base fracture energies per aggregate size

Table 6.8 – Comparison between maximum and low-peak stress

Table 6.9 – Comparison between total and initial energy

Table 6.10 – Improvement in initial fracture energy

Table 6.11 – Improvement in total fracture energy

List of figures

Figure 3.1 – BS flexural test configuration

Figure 3.2 – Uniaxial tensile test concept

Figure 3.3 – Comparison of Hognestad and Eurocode model

Figure 4.1 – Solid65 element (ANSYS, release 5.5)

Figure 4.2 – Solid285 Element (Ansys, released 14.0)

Figure 4.3 – Fibre generation diagram

Figure 4.4 – Model generation sequence

Figure 4.5 – Stress and strain variation along YY axis Fibre Content 1.0% in weight

Figure 4.6 – Stress and strain variation along YY axis Fibre Content 2.0% in weight

Figure 4.7 – Stress and Strain Distribution along YY axis

Figure 4.8 – Model 1 and 2, low stress results and their projection

Figure 4.9 – Model 2 configuration

Figure 5.1 – Particle size distribution chart. Left) Fine aggregate, Right) Coarse aggregate

Figure 5.2 – Propex ENDURO HPP45 (

Figure 5.3 – Gauges’ set-up

Figure 5.4 – General test arrangement

Figure 5.5 – Specimen T1, raw data and fitting line

Figure 5.6 – Specimen T1, BoxPlot

Figure 5.7 – Specimen T1, refined data with improved fitting line

Figure 5.8 – Data refinement – Specimen T2

Figure 5.9 – Data refinement – Specimen T3

Figure 5.10 – Data refinement – Specimen T4

Figure 5.11 – Data refinement – Specimen T5

Figure 5.12 – Data refinement – Specimen T6

Figure 5.13 – Specimen T4, stress-separation curve with fitting curve

Figure 5.14 – Specimen T5, stress-separation curve with fitting curve

Figure 6.1 – Fitting of linear model to fit T1-3 data

Figure 6.2 – First modification of Hognestad model fitting T1-3 data

Figure 6.3 – First modification of Hognestad model fitting T1-3 data

Figure 6.4 – First modification of Hognestad model fitting T4-6 data

Figure 6.5 – Finite element against experimental results

Figure 6.6 – Linearity of characteristic results

Figure 6.7 – Tensile capacity graphical resolution

Figure 6.8 – T4 and T5 raw data and fitting lines

Figure 6.9 – Fracture energy model, specimen T4

Figure 6.10 – Fracture energy model, specimen T5


Tensile strength from flexural test (MPa)
Derived maximum tensile strength (MPa)
Derived ultimate tensile strength (MPa)
Derived maximum tensile strength of composite (MPa)
Mean experimental tensile strength (MPa)
Interpolated characteristic strength for observed
Mean target compressive strength (MPa)
Stress-separation curve
Gradient to initial tangent (-)
Standard mean tensile strength (MPa)
Applied Stress (MPa)
Ultimate compressive strength (MPa)
Induced strain (-)
Strain at ultimate compressive stress (-)
Strain at maximum tensile stress (-)
Mean experimental strain at ultimate stress (-)
Mean experimental strain at stress (-)
Mean experimental strain at maximum stress of composite (-)
Total fracture energy (N / m)
Initial fracture energy (N / m)
Secant elastic modulus (GPa)
Tangent or dynamic elastic modulus (GPa)
Crack opening, separation (mm)

1 /                      Introduction

Composite materials found their strength on the concept of the juxtaposition of two or more materials characterised by complementary properties. The potential of such a simple technique is fascinating and at the same time a challenging research topic.

Countless examples of the application of this technique are available in everyday life. The concept of composite material or object has been in use for centuries and it spreads across most engineering fields. If carried to extremes, everything can be considered to be a composite formed by basic chemical elements.

Moving back into the macroscopic realm, the project analysis is focused on a very common example of composite material, reinforced concrete (RC). Concrete reinforced with metal bars has been, over the last two centuries, a common method to improve the behaviour of concrete in tension and widen its applications.

While the principle of reinforced concrete was patented by Joseph Monier in 1867, the first example of a concrete structure reinforced with iron bars was built by Francoise Coignet in 1853. Over the last decades, the concept of “reinforcing” has seemed to shift and to focus on improving concrete at an elementary level. The most recent trend of reinforcing concrete mixtures with fibre-like addition can, in fact, be considered an attempt to enhance the fundamental material properties.

Despite the fact that information on this topic is still limited if compared to other materials, the use of fibre-reinforced concrete (FRC) has seen a quick application growth and it is now regularly employed in numerous structural applications. For instance, FRC is employed in the field of bridges, pavements, and architectural structures (Hassan et al. 2011).

Despite the potentials of FRC and its wide range of applications, a standardised test method to describe the material’s behaviour in tension and compression has not yet been found (Hassan et al. 2011).

1.1 /      Aim and objectives

The aim of the project is to determine and quantify the benefit due to fibrous materials to cement mixtures and to observe the variation of fracture energy between a plain concrete mix and the corresponding modification.

A quantification of the  properties’ enhancement will be obtained from the analysis and modelling of the experimental results. Two models will be proposed respectively to describe the tensile capacity and the fracture energy of a plain concrete and its fibre composite.

A generalised law will eventually be defined to estimate the improved tensile capacity as a function of the mix characteristic compressive strength (cylinder). The validity of the law will be limited to C30/37 concrete with a content of fibre of 1.0% and 2.0% in weight. Similarly, a relation for the concrete strain at maximum stress will be derived.

While the first model is considered to produce an estimate close to the improved tensile capacity, the second model proposes a graphical resolution to simulate the behaviour of the composite and to estimate the order of magnitude of the fibre addition benefit.

1.2 /      Methodology

An experimental approach will be adopted to achieve the project aim. A non-standardised direct test will be designed for the testing of the material along with the shape and size of the specimen to be tested.

The validity and the ground for comparison of the experimental data will be based on the results obtained from the control cubes and cylinders prepared from the same concrete batch and tested according to the British Standards “Testing of Hardened Concrete” (BS EN 12390:2009).

Prior to the study of the results a simple statistical analysis of the data will also be performed in order to reduce potential error due to equipment and testing procedures inaccuracies.

The purged data will then be used as the starting point for the definition of the models.

Eventually, a computer model, prepared with ANSYS finite element analysis software, will be introduced to determine whether simple modelling techniques based on the macroscopic structure of the material can be considered sufficient to explain the material response to stress.

2 /                      Literature review

In this chapter, a critical analysis of previous works on the topic of fibre-reinforced concrete will be presented. All aspects analysed derived from either the necessity to verify certain assumptions or from general research on the subject. This review will keep into account both supporting as well as opposing facts.

The existence of composite material is not new. Numerous examples of composites materials can be found throughout history. A famous composite example can be dated back to the Mongols, inventors of the so-called modern composite bow (approximately 1200 BC).

The genius of the Mongols was to understand that two different material characterised by opposing properties could have been juxtaposed to produce a third material with enhanced properties. The composite bow was, in fact, the result of a combination of horn and sinew to respectively work in compression and tension.

Despite the topic of composite materials having been the focus of many research papers and publications, material on its testing in direct tension test are limited.

Several publications can be found on the modified indirect tension test – or Barcelona test – such as “Splitting tensile test for fibre-reinforced concrete” (Denneman et al., 2011).

But, as outlined by van Vliet et al. (2000), when some difficulties are experienced in the uniaxial tensile test, this should not imply one reverting to an apparently easier experiment to be performed, but with more difficult results to be interpreted. Denneman’s experiment, besides its set-up complexity, produces results difficult to interpret.

For example, Denneman’s model aims to calculate the load-deformation curve in order to derive the peak stress on the main failure plane while also keeping into consideration the secondary cracking. The reliability of a test aiming to identify two or more parameters is almost impossible due to the material variability, scatter, error of testing method, and mathematical ill-conditioning (Bazant, 2002).

Besides the complexity of the results analysis, an important aspect that doesn’t appear to be clear from Denneman’s method is the observation of FRC behaviour after cracking and thus the benefit of the fibre reinforcement to the concrete composite.

The modified test proposed by Denneman does not consider the modified nature of the material, and he forgets to account for an important aspect such as the improved ductility of FRC. Testing a ductile material through a standardised splitting test configuration is not possible (Denneman et al., 2011) as it will only test the strength of the composite matrix. In addition, Rots et al. (1989), while analysing the concrete fracture behaviour in direct tension, has noticed that it is difficult to obtain realistic results due to the composite nature that can neither be classified as elastic-perfectly brittle nor as elastic-perfectly plastic.

The author agrees with the vision proposed by Marcel et al. (2000) and understands the need of deriving a standardised test aims to compare the results obtained from between the conventional Barcelona test and the direct tension test, limited to Mode I failure (Splitting). This approach should provide reasonable findings

2.1 /      Specimen profile design

The original mould design previously used in this typology of experiment within the London South Bank University was a prismatic 100 x 100 x 400 mm specimen with a one-directional necking at its middle.

Major findings on the size-effect of the specimen was found in the extensive work of Bazânt (1997) and Carpineti (1994). It was, in fact, seen that a size-effect does yet exist, and it is due to a number of factors amongst which the heterogeneity of concrete.

From visual observation of the existing mould, no particular relations were found between the necking’s – the narrower section of the specimen – geometry and the fibre’s length. The necking was reduced to a cross section of to 50 x 100 mm2 to maximise the likelihood of failure in such area.

Several theories have been proposed about the specimen size. According to Bazânt (1997), large specimens are required to guarantee a homogeneous fibre distribution in the area where the failure needs to be induced (“ligament area”). This specimen size, however, is also subjected to large vertical flexure due to the irregular crack development, which can be negligible on small specimens (Barragan et al., 2003). Small specimen results are in turn hindered by the boundary conditions. Limitations due to the testing equipment and available materials must also be accounted for in the design of the specimen.

As a compromise, the specimen size was defined both in relation to the fibre length and the aggregate’s maximum size. The smallest necking dimension was taken as the average fibre length in order to reduce fibre balling and maximise random fibre distribution (Hassan et al., 2012). At the same time, the average of the necking dimensions was seven times the aggregate maximum diameter as per guidelines proposed by van Vliet et al. (2002) who suggested a minimum specimen size eight times the maximum aggregate size.

This was further implemented by limiting the necking to one dimension only. The depth of the necking was fixed at 100 mm (twice the fibre length) allowing for free fibre rotation in that direction.

Reducing the area of the necking in order to increase the chance of failure within this region is not a new technique. The method does exploit the definition of stress of force per unit area and by decreasing the available area it increases the stress. The critical aspect lays in the procedure adopted to produce the necking. The two main methods are by scaling of a cured prism or by casting.

The scaling of a notch, according to Marcel et al. (2000) besides changing the effect of the notch can also cause damage to the specimen. The ideal solution would involve a modified geometry with cast production process. The bays forming the necking would have a circular profile that, as proven by van Vliet (2002), reduces the stress concentration at the ends increasing the probability of failure within this area.

As a circular bay profile can involve geometry discrepancies due to manual production processes, a modified trapezoidal bay will be used instead.

2.2 /      Concrete mix design

Concrete is a heterogeneous material (van Mier et al. 2002) and as such it is possible to distinguish its aggregates. By understanding that the heterogeneous concrete internal structure can influence the stiffness distribution within the material, Van Mier predicts a possible non-uniform stress distribution, thus a non-homogeneous stretch during testing.

However, it might be objected that this technical problem is valid for small ratio of aggregate size to specimen’s dimensions. And, more importantly, this might not be the case for fibre-reinforced concrete where a more complex matrix is employed.

In order to avoid problems related to the particle size, as previously observed, a coarse aggregate of nominal maximum size of 10mm is used for all the concrete mixtures, as it is also suggested in “Fibre Cements and Fibre Concretes” (Hannant, D.J., 1978).

No particular considerations were found about the water-cement ratio and the workability of the mixture.

3 /                      Background

Before the advent of fibres as a concrete enhancing material, a limited amount of research was carried out to understand and model the behaviour of concrete in tension.

Common testing to determine concrete’s tensile strength are:

  • Three point flexural test,
  • Splitting test,
  • Direct (or axial) test.

The flexural test according to BS EN 12390-5:2009 involves the testing up to failure of a prismatic specimen with dimensions meeting BS EN 12390-1:2012 requirements. The test is carried out through the application of a load at the two third-points of the specimen while being simply supported at both ends. The test configuration is shown in figure 3.1.


Figure 3.1 – BS flexural test configuration

From the load causing failure of the prism the flexural strength is obtained as


Where F is the failure load, I is the distance between the supporting rollers and d the cross sectional characteristic dimension for a squared section prism.

Following the Eurocodes guidelines, the flexural tensile strength can be used to calculate the axial tensile strength as


It can be easily calculated that for a standard specimen 100 x 100 x 350 mm, the tensile axial strength is 60-70% of the flexural tensile strength.

The second method for the determination of the tensile capacity of concrete is the cylinder splitting test. The test involves the crushing of a d x 2d cylinder, where d is the diameter in mm, as per BS EN 12390-1:2012 requirements.

The specimen is tested to failure through the application of a load perpendicularly to its longitudinal axis.

As this testing method is used in the experiment as a control mean, its set-up and procedure will be described in the experimental section of this paper.

The last method is the direct axial test. Despite the test being mentioned in the British Standards and Eurocodes, at this time there is no standard testing procedure defined.

However, as the test’s name suggests, the result obtained from the test are likely to be the actual direct stress value. The force is, in fact, applied along the longitudinal axis of the specimen inducing a stress on the sample cross-sectional area. This stress, by definition, is the direct (or normal) stress due to the tensile force applied.

On this very concept, the test set-up adopted for the main laboratory experience was based, the initial experiment configuration can be seen in figure 3.2.

Tensile Concept

Figure 3.2 – Uniaxial tensile test concept

3.1 /      Stress-strain curve

Plain concrete is a material with a complex behaviour (Mosley, B., Bungey, J. and Hulse, R., 2007). It does not behave linearly, and its properties vary depending on the cement, the aggregates and the amount of water used to bound everything together in a macroscopically heterogeneous material.

Several models describing the behaviour of concrete under stress have been described by different scientists and engineers.

The most noticeable ones are:

  • Hognestad (1951). His model described both the ascending and the descending part of the curve. Up to the maximum compressive strength, the curve is described by a parabolic function thereafter by a straight line.
  • Kent and Park (1971). In their model, they generalised Hognestad’s model to obtained a more practical model.
  • Popovics (1973). He aimed to describe the stress-strain curve as a single function depending on three parameters (f’c, εco, Ec)

The model that will be adopted for the modelling of concrete in compression are:

  • Hognestad (1951) with the parametric values to meet EC,
  • and, the proposed EC model for non-linear structural analysis.

3.1.1 /                        Hognestad model

Experimental results (Rots, J. G. et al., 1989) show that concrete is characterised by a non-linear behaviour in compression. Hogenstad’s method aims to describe the material behaviour through the use of simple functions. A parabolic shape defines the stress strain relationship of the concrete as the ratio between the applied stress and the maximum compressive stress, and it is defined as

     for any    3.3

This relation is valid only up to the failure strain (strain at first crack development) i.e. for the ascending part of the stress-strain curve.

For the descending part, or softening region, the above ratio is instead described as a straight line down to failure.

     for any    3.4

Where the above notations stand for:

σ applied stress,

σcu peak compressive strength,

ε induced strain due to stress σ,

ε’0 strain at peak compressive strength, failure strain,

εcu ultimate strain

The initial elastic modulus, or dynamic modulus, can also be determined as the tangent at the origin to the curve.


The EC values for the designed concrete mix – C30/37 – were adopted to calculate Hognestad stress-strain curve, summarised in table 3.1 (EC extract available in appendix A)

Table 3.1 – Applied model values
Description (Hognestad, EC2) C30/37
Peak compressive strength (σcu, fcm) – MPa

Strain at peak compressive strength (ε’0, εC1) – ‰

Ultimate Strain (εCU, εCU1) – ‰




3.1.2 /                        Eurocode model

Differently from Hognestad’s, the method proposed in the EC guidelines is described by a single equation for both the ascending and descending part of the curve. The EC curve equation recalls the European Concrete Committee (CEB) equation for short-term loading.




The Eurocodes instead define η and k as


Inserting the EC values (Table 3.1) for the designed concrete mix into equations 3.3 and 3.4 for the Hognestad’s model and 3.5 for the Eurocode model, the graph in figure 3.3 was obtained.


Figure 3.3 – Comparison of Hognestad and Eurocode model

3.2 /      Fracture energy

The fictitious crack model by Hillerborg (1976) is the most popular model adopted for description of the behaviour of a finite-size fracture process zone. This model is herewith proposed as it does not require knowledge of fracture mechanism (Bazânt, 2002).

The applicability of the Hillerborg model presumes certain assumptions: (Bazânt, 2002)

  • the failure is due to the tensile stress  (failure mode I),
  • stress is applied normally the crack surface and no other external stresses are involved,
  • the fracture process zone is of a finite width (localized cracking) and it can be described by a fictitious crack surface perpendicular to the stress.

An indicator of the total energy dissipated (GF) by the fracture mechanism is the area subtended by the entire stress-separation curve of the softening region. The value obtained is the energy dissipated per unit area of the crack until its full separation.

Hillerborg also defines the area under the initial tangent to the stress-separation curve and names it Gf. This value according to Bazânt, is the main parameter that determines the maximum allowable load of the structure. GF and Gf can mathematically be described as

     ,    3.11, 3.12

Where the above notations stand for:

GF total fracture energy

Gf area under initial tangent

f’t tensile peak stress

w crack separation

f(w) stress-separation curve

σ’0 gradient of initial tangent

According to Planas et al. (1992) and Guinea et al. (1994) an experimental value of the ratio between these two values can be seen to be


and it can be used to determine the shape of the bilinear softening curve describing the fictitious crack model.

4 /                      Finite element modelling process

4.1 /      Introduction

Within polymer fibre reinforced concrete, as well as for all the other typology of fibre-reinforced cementitious composites, fibres and matrix are mixed together to achieve an enhanced material.

Many factors are known to be determining the behaviour of the composite material. For instance, several studies can be found in the literature on the importance of fibre-matrix interface (Owen  et al., 1968, Paramisivan, 1984, and Pedziwiart, 2009), the typology and the geometry of the fibre (Uomoto, 2002) or the influence of the concrete mix used for the composite’s matrix.

This chapter will describe the process behind the design of a computer model aiming to study the behaviour of the matrix depending on the fibre content.

The purpose of the models is to estimate the material behaviour for all the three concrete mixes – Plain, 1% and 2% of content fibre in weight. Two estimates will be derived from a model based on a one-element structure that accounts for both the concrete and the fibre behaviour, the second will instead use a two-element structure simulating the fibre distribution in a matrix.

The models prepared for the analysis are as follows:

  • Model 1: SOLID65 with reinforcement values and plain concrete elastic modulus;
  • Model 2: SOLID65 and SOLID285 respectively modelling the matrix and the fibre components.

Due to the complexity of the geometry of a composite material model, a reduced simulation will be carried out to determine the validity of the model results. This simulation will be limited to 10% of the fibre content to reduce the computing time for each operation.

4.2 /      Finite element analysis by ANSYS

4.2.1 /                        Concrete model

The properties of concrete vary depending on the design mix characteristics and they are different whether the material is subjected to tension or compression.

Although the behaviour of concrete can be assumed elastic for low stress values, up to 30-40% of tensile and compressive capacity (Bhatt, MacGinley and Choo, 2005), the same is not true for higher stress values at which the concrete shows a behaviour completely non-linear.

Creating a verisimilar model for complex material such as concrete can be challenging for a professional modeller. For this reason, ANSYS finite element software was chosen for the modelling as it already incorporates an element able to simulate concrete behaviour.

Figure 4.1 – Solid65 element (ANSYS, release 5.5)

The ANSYS designated element for modelling is SOLID65 (Figure 4.1). This element is an 8-node solid in which all nodes are defined by three degrees of freedom, translation parallel to the X, Y, and Z axis. In addition to allowing for approximation of  the plain material behaviour, this element can also account for the reinforcement presence through the material “real constants”. This will be discussed further on in the design of the first model.

SOLID65, as a solid material, requires the definition of isotropic material both in its elastic (linear) and plastic region (multilinear). The material’s properties adopted throughout all computer models for SOLID65 are listed in table 4.1 and they are obtained from the stress-strain model previously presented.

Table 4.1 – Concrete properties
Isotropic linear   Isotropic multilinear  
Young’s Modulus – Ec

Poisson ratio – ν


3.3 x 1010
Set Strain – εc
Stress – σc












2.9 x 106


9.0 x 106

16.775 x 106

27.859 x 106

32.581 x 106

33.300 x 106

4.2.2 /                        Fibre model

The model for the fibrous component of the composite was modelled using SOLID285, a solid element with only four nodes (Figure4.2), in order to reduce the number of equations and boundary conditions necessary for the computer simulation. The element definition according to the National Agency for Finite Element Method and Standards (NAFEMS, London) states that the element has a linear displacement and hydrostatic pressure behaviour.


Figure 4.2 – Solid285 Element (Ansys, released 14.0)

This is ideal for the modelling of the fibres as it allows for linear deformation of the fibres up to matrix cracking.  The hydrostatic pressure can be replaced with the value of volume change rate for quasi-incompressible materials. This can be assumed to be accurate up to the critical stress due to the large difference between the concrete and polymer fibre tensile strength.

The material properties adopted for SOLID285 were obtained from the supplier technical specifications for the product here listed in table 4.2.

Table 4.2 – Fibre properties – PROPEX Enduro HPP45
Material properties   Shape properties
Young’s Modulus – Ef


Tensile Strength

Poisson ratio* – ν


Polymer type

Specific gravity


Not available
0.3 (Assumed)









* – Poisson ratio cannot be uniquely determined. Its values changes depending on the fibres’ orientation

4.3 /      Material model

4.3.1 /                        Fibre distribution approximation

The model consists of two separate volumes: the matrix and the group of fibres. The added fibres can be assumed to be homogeneously distributed within the matrix as the concrete paste is mixed for five minutes as per supplier specifications.

The distribution of fibre within the matrix was determined using the RANDfunction available in Office Excel by Microsoft. The function, which is designed to randomly generate a number within a given range, was used to determine five values, respectively:

  • the x, y, z coordinates of the fibre centroid, within a 0.1 x 0.1 x 0.1 m volume;
  • the angle θ for the rotation about the Z-Z axis, between 0 and π radiant;
  • and, the angle γ for the rotation about the Y-Y axis, between 0 and π radiant.

Using simple geometrical properties, and the values previously determined in excel the coordinates of the fibre’s extremes were determined as two points equidistant from the centroid along the direction defined by the angles of rotation θ and γ (Figure 4.3).


Figure 4.3 – Fibre generation diagram

Once the method to control the position of the fibres in the matrix was determined, the number of fibres had to be estimated from the properties of the fibres for each fibre content fraction.

The fibre’s shape can be described to have a two-dimensional wavelength shape with a circular cross section of diameter 0.95mm and 45mm length. As the effective length of the fibre can not be easily calculated, the fibre was idealised as a cylinder of diameter 1.5mm and length 45mm. This approximation is possible only because the bond element between fibres and matrix is not kept into account for the purpose of this project.

The number of fibres for the two fibre contents was determined as follow (Table 4.3).

Table 4.3  – Derived fibre’s properties
Fibre Dimensions   Fibre properties









7.952 x 10-8

7.236 x 10-5

Specific gravity


Density *

Kg / m3



Number of fibres per concrete mix per m3
Mix type Weight

kg / m3

Fibre content

(% of weight)

Total weight of fibres

Kg / m3

Estimated number of fibres

( – / m3)

Plain Concrete (PC)

Fibre 0.5%


Fibre 1.0%














Number of fibre per 0.1 x 0.1 x 0.1 m3 volume
Mix type Number of fibres
Plain Concrete (PC)

Fibre 0.5%


Fibre 1.0%





* – water density: 1,000 kg / m3

4.3.2 /                        Reduced model

Finite elements models with complex volumes derived by random factors, as in this experiment, can  produce a substantial amount of elements because of the lack of regular patterns in the geometry. A complex model also implies a considerable amount of computing time for each solution.

In order to prove the validity of the model strategy that is going to be adopted throughout the finite element analysis, a reduced model has been designed to determine the typology of results.

The following modifications were adopted:

  • the sample volume was taken as an eighth of the initial volume. The new test volume was 0.05 x 0.05 x 0.05 m;
  • The fibres cross-sectional area was assumed square with area equal to the circular one;
  • and, only 10% of the fibre content was modelled.

4.3.3 /                        Computer model generation

The generation of the material model followed three basic steps for all the material samples. Firstly, the fibre volumes where lofted along the lines drawn using the coordinates of the fibres’ extremes obtained from Excel; these were then copied and offset outside of the 0.05 x 0.05 x 0.05 m volume of the matrix (Figure 4.4.a). Secondly, the volume of the matrix was created (Figure 4.4.b) and the fibre’s volumes were subtracted to the matrix volume (Figure 4.4.c); the two were superimposed to match the original fibre position within the matrix (Figure 4.4.d).

REDUCED_FibreOffset REDUCED_FibreMatrix
a) b)
REDUCED_FibreMatrix_Subtracted REDUCED_FibreMatrix_Superimposed
c) d)

Figure 4.4 – Model generation sequence

a) Fibre modelling and offset, b) matrix modelling, c) fibre volume subtraction,

d) Superposition of Fibre and matrix volumes

Once the volumes were modelled, the material’s properties, previously listed, were applied accordingly.

The model was then further optimised through the application of the following load and boundary conditions:

  • symmetry conditions were applied on the volume faces parallel to the applied load, this allowed for a more refined result without increasing the number of model elements;
  • the load was applied as a surface pressure on the top surface. Particular attention was given to the sign convention. Positive defines a compressive and negative a tensile strength,
  • and, the translation in the direction of the applied load was assumed to be zero to simulate a vertical restraint.

4.3.4 /                        Reduced model results

Several combinations of loads and fibre content were performed for the reduced model. For a chosen load of -50,000 Pa the sets of results from the  different sets were graphically compared (Figure 4.5).

a) b)

Figure 4.5 – Stress and strain variation along YY axis
Fibre Content 1.0% in weight

From figure 4.5, it appears that the model is behaving as it was assumed it would have.

The distribution of stress does, in fact, show an increase of stress concentration around the fibre. In the graphs, this can be recognised in a minor peak before the major ones of opposite magnitude. The major stress peaks occur at the fibre location, only within the fibre element, this is due to a change of material and, thus, of the material elastic modulus.

A similar and opposite behaviour is seen in the strain variation along the YY axis (Figure 4.5.b). By recalling Hooke’s law of linear behaviour for small stresses (Equation 4.1)


It can be observed that the elastic modulus is close to 30 GPa as per plain concrete. Similarly, the stress and strain behaviour for the 2% fibre content model were obtained (Figure 4.6).

a) b)

Figure 4.6 – Stress and strain variation along YY axis
Fibre Content 2.0% in weight

4.3.5 /                        Validity of finite element model

From the analysis of the results just obtained, the model can be considered acceptable. Hooke’s law is respected throughout the tests for the relatively small stress applied and, furthermore, the expected behaviour has been satisfied as is shown in the previous graphs (Figure 4.5-6) and the stress distribution below, plotted via ANSYS on four cross sections parallel to the load direction. (Figure 4.7)

Set1_1000N_YStrain_Section001 Set1_1000N_YStress_Section001
Set1_1000N_YStrain_Section002 Set1_1000N_YStress_Section002
Set1_1000N_YStrain_Section003 Set1_1000N_YStress_Section003

Figure 4.7 – Stress and Strain Distribution along YY axis

4.4 /      Tensile specimen model

The model of the tensile specimen, different from the one employed in the determination of the validity of the method, is designed following the real dimensions of the specimen to be tested. As anticipated, two models will be here introduced.

4.4.1 /                        Model 1 – SOLID65 with reinforcement values

For the simulation of the fibre addition to the matrix, modified “real constants” valueswere adopted for the experiment. From the concrete design properties and the fibre characteristics, the following parameters were used respectively for the 1.0% and 2.0% fibre content in weight.

Table 4.4 – Real constant parameters
Material geometry
ANSYS notation 1.0% 2.0%
Volume ratio

Orientation angle

Orientation angle










Material Properties
Young’s Modulus – Ef


Tensile Strength

Poisson ratio* – ν


Polymer type

Specific gravity


Not available
0.3 (Assumed)



The model results for low value of stress and their projection is shown in figure 4.8.

Figure 4.8 – Model 1 and 2, low stress results and their projection

4.4.2 /                        Model 2 – SOLID65 with modelled fibres

The second model uses the material property values available in table 4.1 and 4.4. The symmetry of the model was, however, exploited to reduce the computational time. The configuration of the model was as per figure 4.9.

The simulation was carried out for the same loads applied to model 1.

4.4.3 /                        Model limitations

Despite the fact that the model appeared to fairly approximate the material response to applied stress in the validation of the reduced model, this has been observed not be adequate for high stress values close to the plastic region of the material. A possible solution to this restriction will be proposed in the conclusion chapter.

Figure 4.9 – Model 2 configuration

5 /                      Experimental process

5.1 /      Mould design

Substantial influence on the quality of the specimen to be tested comes from the quality of the casting mould.

Alternatively from the mould used in previous tests, modifications were brought to the design in order to minimise damage to the sample and to simplify the de-moulding process.

The main differences were as follow:

  • mould sides were designed as one piece with the base in order to reduce the assembling and breakdown time,
  • the inside of the mould was lined with High Impact Polystyrene (HIPS) to eliminate the gap between the sides and notch elements, and also to reduce concrete adhesion to the mould,
  • and, the elements marking the clamps’ groves, previously made of timber, were replaced with a silicone rubber equivalent that can easily be removed without damaging the clamping area.

5.1.1 /                        Design process

The design steps were as follows:

  • the base of the mould was divided along its longitudinal axis, the sides were fixed to the base through a 5 mm slot and 50 x 50 x 8 mm steel angled profile;
  • the two sides were then reconnected through a 300 mm wide board placed beneath the side elements and connected with M6 bolts;
  • the gap between the two sides was made watertight with a 4mm silicone strip as well as all other small gaps due to production defects;
  • end plates and a mid-divider were then inserted in order to restrain side’s movement due to concrete pressure;
  • and, the mould interior was lined with HIPS to eliminate any gaps between mould components and to facilitate the de-moulding process.

Detailed drawings of the mould design are available in appendix B, C, and D.

5.2 /      Materials and mix details

The materials used for the test were: Portland-Limestone cement (Mastercrete – LAFARGE) that meets BS EN 12390-2:2009, 10 mm maximum size coarse aggregates, natural sand, and polypropylene fibre with a density of 910 kg / m3.

Moisture content and sieving tests were performed for both aggregate types to confirm the suitability of the material and their initial condition. From the sieving test results (Figure 5.1) and observation, both the sand and the gravel appeared to be uncrushed, non-angular and well graded minimising the opportunity of void forming.


Figure 5.1 – Particle size distribution chart. Left) Fine aggregate, Right) Coarse aggregate

An increase in water content can be assumed from the findings of the moisture content test. Fine and coarse aggregates were, in fact, found to have a water content – in grams – respectively of 0.37% and 0.51% (Table 5.1). These characteristics can be said to increase the water/cement ratio of approximately two percentage points for the designed concrete mix.

This factor was, however, not included due to the difficulty on monitoring the water content of the aggregates during their handling.

Table 5.1Results of moisture content test
Fine aggregate sample   Coarse aggregate sample
Weight of beaker (BK) – grams








Wet sample – grams
(+ BK)



(Mean) 1517.57 (- BK) 1117.76 (+BK)







24h over-dried sample – grams


1513.39 1113.58 1443.72


1443.75 1049.97
Moisture content g 4.18


Moisture content g 5.43


The fibres adopted were 0.95 mm in diameter and 45 mm long (Full supplier specifications available in table 4.2), characterised by a two-dimensional sinusoidal shape (Figure 5.2)

EnduroHPP45 macro

Figure 5.2 – Propex ENDURO HPP45 (

The mix used was valid for normal weight concrete with normal compressive of 37 MPa at 28 days, the design followed the guidelines proposed by the Building Research Establishment (BRE) in Design of Normal Concrete Mixes.

Details of the mixes can be found in table 5.2.

Table 5.2 – Concrete mixes details – water cement ratio 0.5
Mix Cement

Kg / m3


Kg / m3

Fine Aggregate

Kg / m3

Coarse Aggregate

Kg / m3

Fibre content

Kg / m3

Plain Concrete


Fibre 0.5%


Fibre 1%

















Fibre contents of 1% and 2% of the concrete weight were used throughout the experiment as well as a plain one for the control of the results. The respective volume fraction for the two mixes was calculated to be approximately 0.5% and 1.0%.

5.3 /      Test specimen and procedure

The materials were mixed using a laboratory concrete mixer adding in order: coarse aggregates, fine aggregates, and cement. The dry components were mixed until they reached an even tone. Then, water was gradually added in, and mixing was carried on until the paste looked homogeneous, approximately 2 minutes.

Eventually, the fibres were added to the concrete paste. In this step particular care was given to the addition of fibres as adding the fibres in bunches could have caused the balling of fibre within the mix.

Along with the tensile specimens, for which the special mould was produced, two control cylinders (150 diameter x 300 mm), and three control cubes (100 x 100 x 100 mm) were cast for each batch. The number of control specimens was limited due to the limited number of moulds available in the testing facilities.

5.3.1 /                        Control cube

The control cubes were cast and tested following BS EN12390-3:2009. Each 100 x 100 x 100mm was tested in a compressive test plant with a maximum load capacity of 3,000 kN. A constant load rate of 3,500 N/sec was applied throughout the experiment.

From the testing, the following results were found.

Table 5.3Cube testing results
Mix Cube #


Mean Strength


Standard deviation from Coefficient of Variation – Cv
1 2 3 Mean 47
PC 49.92 52.96 44.53 49.14 4.269 5.007 8.7%
FRC_1.0 52.43 55.27 50.97 52.89 2.187 4.1%
FRC_0.5 44.02 46.11 43.47 44.53 1.393 3.1%

According to the mix design, the target strength for the designed mix was 43.1 MPa, however, for the determination of the concrete mix quality the cube target strength of 47 MPa, proposed in the Eurocodes, was considered.

As it immediately appears, the concrete strength of the mix with 0.5% fibre content is 9.4% lower than the same mix without fibres. As this result does not meet the expected strength value, it can be assumed defective.

Further discussion on the validity of the FRC_0.5 batch will be discussed further on along with the splitting test and tensile test results.

The plain and the 1.0% fibre content mixes do, instead, behave as per design estimate. The standard deviation of 5.0 N/mm2 from the target value is, in fact, lower than the 8.0 N/mm2 suggested by the BRE (1988). Furthermore, the coefficient of variation for both PC and FRC_0.5 shows that the data lay within a small range from each other. Hence, despite the low number of samples, they can be assumed to be reasonably characteristic of the design mix.

5.3.2 /                        Control cylinder

The splitting cylinder test was also carried out to determine the tensile strength of concrete. Testing procedure followed BS EN 12390-6:2009 and it was performed in the same testing equipment used for the control cubes with the same load rate of 3,500 N/sec.

The result obtained confirmed the possibility of a mix defect for the FRC_0.5.

The splitting cylinder test results are presented in table 5.4.

Table 5.4Cylinder testing results
Mix Cylinder #


Mean strength MPa
1 2
PC 2.521 3.108 2.81
FRC_1.0 3.766 4.032 3.90
FRC_0.5 2.917 2.927 2.92

Due to the limited amount of samples tested it would be erroneous to determine the variation of data from standard deviation and coefficient of variation.

On the other hand, it can still be observed that the values for PC are close to the target splitting tensile strength of 3.2 MPa proposed in the EC, and a 38.8% strength increase in the FRC_1.0 can also be noted.

5.3.3 /                        Tensile specimen

Although a standardised method outlining the testing procedure results is not yet available, direct testing of the tensile strength of concrete is considered a possible method to determine the tensile capacity of concrete by both Eurocode and British Standards.

In order to allow for a more reliable future comparison, a detailed experiment set-up will be herewith described.

The tensile specimens to be tested were produced using the concrete mix previously described and cast in the specially designed mould. They were allowed to rest for 48h prior de-moulding, and then left to cure in water for approximately 26 days.

On the testing day the samples were removed from the curing tank and once dry they were fitted with two displacement sensors, also know as Linear Variable Differential Transformers (LVDT). The sensors were placed in the middle on the narrow side of the necking covering a 75mm length. A diagrammatic explanation of the gauges set-up is shown in figure 5.3.

Figure 5.3 – Gauges’ set-up

Each test specimen was equipped with two clamping systems, one at either end of the specimen, after that the instant adhesive used to install the gauge mounting blocks and angled plate had set.

The clamping system was composed of two metal plates slotting into the 7 x 11 mm grooves bridged at their outer edge by a shaft. The horizontal shaft was in turn connected to the primary shaft via a hollow cylindrical element.

The two steel plates were eventually connected together to clamp the specimen with four M12 tightened to the tightening torque of 20 Nm. The machine set-up is shown in figure 5.4. It is important to notice that the horizontal shaft was free to rotate inside the hollow cylindrical element.

Figure 5.4 – General test arrangement

The testing equipment employed for the experiment was a high-capacity static materials testing machine by Zwick/Roell The machine was fitted with tension wedge grips both on the bottom and top cross heads.

The wedge grips were used to clamp the primary shafts from the clamping systems allowing for the test to be carried out.

Once the set-up of the specimen was completed the test was started.

Two computers were collecting the readings respectively from the machine and the two gauges. The first computer, which was the main control unit for the testing equipment, was set to record the cross head translation and the induced load; the second was to collect the gauges’ displacement readings for each load recorded by the main control unit.

Although the testing equipment was able to collect value of translation, these were not used due to the unavoidable slipping of the primary shafts out of wedge grips. Please note, that no alternative testing equipment was available.

Hence the LVDTs’ vertical translation was, along with the correspondent load value, the only data considered. The testing for the plain concrete specimens was carried out up to fracture while the fibre-reinforced concrete specimens up to a machine translation of 5 mm.

5.4 /      Results. Interpretation and reduction

The raw data collected were divided in three arrays:

  • Imposed load – resolution 0.01 kN
  • Extension 1 – resolution 0.001 mm
  • Extension 2 – resolution 0.001 mm

5.4.1 /                        Readings up to failure

A first conversion involved averaging the effective extension values (original value minus the observed reading at beginning of the testing) and dividing it by the gauge length of 75mm. Then, the induced stress in the necking was obtained as the applied load divided by the necking area (0.046 x 0.1 m2).

From the plot of these results, the following behaviour was observed.


Figure 5.5 – Specimen T1, raw data and fitting line

Despite the fact that several results seem to be scattered without any correlation, it also appears a strong concentration of results along, what seems to look like, a curve with a steep initial gradient as proposed by Hognestad for the compression stress-strain curve.

A defined vertical arrays distribution arise from the plot; this is due to the low resolution of the LVDT, in respect to the very small strain of concrete. There are, in fact, different load readings for the same strain value.

Exploiting the built-in curve fitting tool of Microsoft Excel it can be seen that a second order polynomial intersecting the origin is close to the data concentration earlier observed, the value of the coefficient of determination (R2) is low. The validity of the curve is probably weakened by the several values falling out of the concentration area, such values will be hereafter referred to as outliers.

In order to remove most of the outliers and obtain an improved fitting curve, it was decided to consider only the data around each array’s median, i.e., the stress value dividing the array in the higher and lower half.

To do this, the statistical tool of quartiles was introduced. Quartiles are, as the word suggests, those values that divide an array of data in four equal groups. They can be summarised as follow:

  • Zero quartile (Q0), the lowest value of the set, i.e the minimum,
  • First quartile (Q1) divides the lower 25% of data from the upper 75%,
  • Second quartile (Q2), also referred to as median,
  • Third quartile (Q3), divides the lower 75% of data from the upper 25%,
  • Fourth quartile (Q4), the highest value of the set, i.e., the maximum.

From the box plot of the quartiles, it appears that the 50% of each data array around the median was scattered around the observed pattern.


Figure 5.6 – Specimen T1, BoxPlot

The elimination of the lower and higher 25% of the data resulted in a 46 percentage point increase in the R2 value of the fitting curve (Figure 5.7)


Figure 5.7 – Specimen T1, refined data with improved fitting line

The method was, therefore, considered correct, and it was applied to the other 6 sets of values. The following graphical results were obtained.(Figure 5.8-12)


Figure 5.8 – Data refinement – Specimen T2
Top) Raw data, Middle) BoxPlot, Bottom) Refined data and best fitting line


Figure 5.9 – Data refinement – Specimen T3
Top) Raw Data, Middle) Box plot, Bottom) Final refined data and best fitting line


Figure 5.10 – Data refinement – Specimen T4
Top) Raw Data, Middle) Box plot, Bottom) Final refined data and best fitting line


Figure 5.11 – Data refinement – Specimen T5
Top) Raw Data, Middle) Box plot, Bottom) Final refined data and best fitting line


Figure 5.12 – Data refinement – Specimen T6
Top) Raw Data, Middle) Box plot, Bottom) Final refined data and best fitting line

5.4.2 /                        Reading after failure

Following the failure of the plain concrete specimen, no readings were collected due the instantaneous separation of the specimen at the crack and hence its inability of carrying any stresses. (Appendix F)

This is, however, not true for reinforced concrete where the reinforcement comes into play after the first cracking development. Readings describing this phenomenon, previously introduced as strain-softening, were once again collected from the gauges. However as the crack needs to develop in the range of the gauge for it to collect relevant data, only two sets of data out of three were admissible (Apendix G). T6 did, in fact, fail at the upper bound of the notch causing the gauge blocks to fall, the value from the 0.5% fibre content batch was also ignored as defected.

Differently from the data prior cracking, the actual extension (mm) was used to describe the behaviour of the sample against the applied stress. In fact, following Hillerborg previously introduced definition of fracture energy, it can be defined as the area subtended by the stress (N/mm2, or MPa) separation (mm) curve.

In this instance, no reduction of the data was necessary.

The results were plotted from the peak tensile stress onwards and a second order polynomial passing through the peak stress was fitted using MATLAB curve fitting tool (Figure 5.13-14).


Figure 5.13 – Specimen T4, stress-separation curve with fitting curve


Figure 5.14 – Specimen T5, stress-separation curve with fitting curve

6 /                      Discussion

6.1 /      Tensile capacity. Direct and splitting tests

Limited to the data collected during the experimental experience of this project, the tensile capacity of the designed concrete mix have been seen to benefit of a considerable 26.1% improvement.

Table 6.1 – Direct tensile test comparison
Plain Concrete (PC) 1% of Fibres (FRC_1.0)
ID Strength


ID Strength














Mean 2.117 Mean 2.670
Increase +26.1%

This material improvement can be assumed to be entirely due to the addition of fibres as the two batches were prepared at the same time with the same materials and were cured under the same conditions.

The collected data does, however, not match the expected value of tensile strength of 2.9 MPa proposed in the Eurocodes, the result for plain concrete is, in fact, 37.0% lower than expected.

On the other hand, the splitting tensile strength and the compressive strength were found to be respectively only 12.2% lower and 4.6% higher than the proposed EC values.

Table 6.2 – Results variation from Eurocodes’ values
Strength type Experimental







Splitting tensile

Axial tensile










The high discrepancy between the values of axial strength would normally be an indicator of mistakes in either the testing procedure or the sample production. However, being the EC value theoretically derived from the splitting test (BS EN 1992-1-1:2004, 3.1.2(8)) and keeping into account that there are currently no standards for concrete testing in direct tension, this value is considered a reliable indicator of the actual tensile strength of the designed concrete mix.

6.2 /      Modified Hognestad model

From the comparison of the linear model and the experimental results (Figure 6.1) it appears that the assumed linear response of concrete in tension is not an accurate describer of the material behaviour.

Figure 6.1 – Fitting of linear model to fit T1-3 data

From observation, the pattern described by the collected results recalls the shape of the ascending region of Hognestad’s stress-strain model for concrete in compression. The Hognestad model can, hence, be amended to model the tensile region (Figure 6.2). The modification involved changing the maximum compressive stress value to match the experimental mean tensile strength (2.117 Mpa).

Figure 6.2 – First modification of Hognestad model fitting T1-3 data

The adjusted model shows a stronger correlation with the experimental results than the linear model. It yet shows a moderate gradient at the beginning of the curve closer to the secant value of elastic modulus – EC – than the expected dynamic elastic modulus – E0.

A further amendment to Hognestad model involves matching the failure strain with the strain values corresponding to the maximum tensile stress. This modification considerably improved the quality of the fitting curve up to the maximum stress value, as it is shown in figure 6.3.

This model does, however, not explain the material behaviour after first cracking development.

Figure 6.3 – First modification of Hognestad model fitting T1-3 data

A similar model can be applied to the results of the fibre reinforce concrete samples.

Figure 6.4 – First modification of Hognestad model fitting T4-6 data

6.3 /      Finite element model results

The results obtained from the finite element models, despite the fact that are limited to stress values below 30% of the maximum tensile capacity, show a good approximation of the experimental results. (Figure 6.5)

Figure 6.5 – Finite element against experimental results

On the other hand, considering the extensions of the two fitting lines an approximation of the material behaviour at higher stresses, only model 2 display a good match.

This is very interesting and suggests that the basis of the model simulating the fibre addition to the concrete matrix might be correct, and that further development of the model could produce convergence of the model results.

6.4 /      Linearity of strain variation

Furthermore, a linear correlation appears between both the failure strain and the maximum strain for both the plain and reinforced concrete (Figure 6.6). A value close to half of the characteristic strength was found to be the linearity constant for all mixes and strain values.

Figure 6.6 – Linearity of characteristic results

From the best fitting lines the following characteristics were found

Table 6.3 – Fitting lines’ characteristics  
Line Gradient Constant terms









Average 15,263 (Cv=2.98%)

The average gradient is only 1.8% higher than half of the 30MPa concrete characteristic strength and 2.6% higher than the effective characteristic strength (31.34 MPa) interpolated from Eurocode’s values (Table 6.4).

Table 6.4 – Characteristic strength interpolation
Eurocode values














Please note that this value is arbitrary and it might not be correct for different concrete mixes and it has to be used in relation to the indicated units. Further investigation is necessary to confirm the hypothesis.

6.5 /      Formulation of tensile strength

An indication on the magnitude of the tensile capacity of concrete from direct tensile test can be useful to determine the actual strength of either plain or fibre-reinforced concrete.

The aim of the formulation is to define the tensile strength of a concrete mix on the only basis of its characteristic cylinder strength. For the project purposes, the interpolation of the cylinder strength from Table 6.4 was adopted as the effective characteristic strength of the mix (31.14 MPa).

The generalised formula was based on the following assumptions:

  • a second order polynomial model is considered to be characteristic of the stress-strain behaviour of concrete up to the maximum stress,
  • at zero strain corresponds zero stress,
  • the model reaches a maximum (first derivative equal zero) at the maximum stress.

On this assumption the model was generally defined as it follows.


Assuming the first derivative equal zero maximum stress the following relationship between coefficients a and b can be defined.


ε’0 can be calculated from the experimental linear equation observed earlier for the average maximum tensile strength


that can be approximated to


The corresponding strain for the maximum tensile strength of 2.117 MPa is

millistrain 6.5

For the obtained value of strain the coefficient b can be calculated from 6.2 as


and by forcing the polynomial through (ε’0 ; ft*) the coefficient a and b are univocally determined as

,   6.8

The model can, therefore, be summarised as


Re-calculating the first derivative at peak stress the generalised equation of strain can be derived


By substituting EC elastic modulus formula into 6.11 for the given characteristic strength the strain at maximum stress can be found as per equation 6.13.


Eventually, combining 6.13 with 6.4 the generalised equation for the tensile capacity of plain concrete can be formulated as per equation 6.14.


Equation 6.14 can be considered a good generalised formula for the determination of the tensile strength of the designed concrete mix along with equation 6.13 for the corresponding strain.

Recalling figure 6.6, another observation can be made, the linearity of the results for the 1% fibre-reinforced samples appears to lay in between the relations of the fracture strain to maximum stress and the ultimate stress to the ultimate strain. From the analysis of the three equations’ constant terms, it was in fact observed a proportionality.

Table 6.5 – Fitting lines’ constant terms
Line Constant term



0.7576 (3 γ)

0.2732 (1 γ)

0.4955 (2 γ)

The value of γ was found equal to 0.2544 [MPa] from which a more generalised layout of the previous equations can be derived (Table 6.6)

Table 6.6 – Generalised linear equations
Line Equation
PC_Peak 6.15
PC_Ultimate 6.16
FRC_1.0_Peak 6.17

As the 0.5% fibre content mix was not considered acceptable, it is not possible to describe a direct relation from the mean peak tensile strengths of plain to fibre-reinforced concrete.

A graphical solution to determine the benefit of the fibre addition was hence proposed. The graphical resolution is presented in the figure below.

Figure 6.7 – Tensile capacity graphical resolution

The method can be summarised as it follows:

  • project the peak strength onto the ultimate stress line,
  • find corresponding stress,
  • and, project on to FRC_1.0 line.

Mathematically, the relation can be expressed as the plain concrete tensile stress plus half the vertical distance between the two lines, i.e γ. From observation


Inserting this value into PC_Ultimate line (Equation 6.16)


Inserting the above value of strain into FRC_1.0_Peak (Equation 6.17) allows for the final formulation of tensile capacity of fibre-reinforced concrete to be defined.


The above solution delivers a value very close to the one found in practice with a 4% difference.

6.6 /      Fracture energy for fibre reinforced concrete

In accordance with Hillerborg definition, the area subtended by the stress-separation curve represents the total fracture energy dissipated by the system with the development of the crack.

As a record of the behaviour post fracture could not be collected for the plain concrete specimens as the gauges were not sensitive enough to record the material behaviour during the quasi-instantaneous fracture failure. The fracture energy for plain concrete was calculated according to the equation proposed by Hilsdorf and Brameshuber (1991) (FIB, 2009).

for 6.21
for 6.22


GF:  Fracture energy [N / mm]

GF0: Base value of fracture energy, table 6.6 [N / mm]

FCM: Mean compressive strength [MPa]

FCM0: 10 [Mpa]

Table 6.7 – Base fracture energies per aggregate size
dMax – mm 8 16 32
GF0 – N /mm 0.025 0.030 0.058



The mean compressive strength was taken as the proposed EC mean cylinder strength for C30/37 concrete.


On the other hand, the fracture energy for the fibre-reinforced specimens cannot be calculated using a simple formula, and a model based on Hillerborg’s Fictitious Crack model was adopted. As previously introduced, the fracture energy can be defined as the area subtended by the stress-separation curve.

Figure 6.8 – T4 and T5 raw data and fitting lines

Although the collected results are limited to a local extension of 2.5 mm, the curves appear to fairly approximate the behaviour of the material. It does yet show a gradual decrease in gradient following the initial drop, as observed by Grassl, P. et al. (2010), Barrang, B. E. et al. (2003), and Bazânt, Z. P. (2002).

The curve was then simplified to meet the bilinear model proposed by Hillerborg (1976). The final model was taken as two lines intersecting at the projection of the first low-peak onto the fitting curve and

  • Intersecting the maximum tensile stress value,
  • And, tangent to the fitting curve at the low-peak projection.

The model proposed is shown in figure 6.8 and 6.9 for the two sets of data in consideration.

Figure 6.9 – Fracture energy model, specimen T4

Figure 6.10 – Fracture energy model, specimen T5

It is interesting to notice that in average the low peak stress is 29.7% lower than the maximum tensile strength (Table 6.8)

Table 6.8 – Comparison between maximum and low-peak stress
Test Maximum Stress (MS)


First Low-Peak (FLP)












Mean 70.3%

These values are more than two folds out of the range observed by Bazânt (2002) of 15 – 33 % of the maximum stress. This variation is assumed to be due to the fibre reinforcement. Bazânt’s model does, in fact, not account for the presence of fibres that are theoretically assumed to improve the ductility of the material, hence increasing the amount of energy necessary to obtain ductile fracture.

Using the two bi-linear models proposed for T4 and T5 specimens an estimate of the fracture energy for fibre-reinforced concrete can be determined.

Table 6.9 – Comparison between total and initial energy

N / m


N / m

GF / Gf









Being the ratio GF / Gf constant allows to considered the proposed model a good indicator of the material behaviour after failure, further test should be performed to verify this statement.

Gf,H according to Hillerborg definition can also be calculated. Assuming the line through the maximum stress and the first low-peak a good approximation of the initial tangent to the curve the following values can be defined (Table 6.10).

Table 6.10 – Improvement in initial fracture energy

N / m


N / m

GF,H / Gf,PC









The ratio GF.H / Gf,PC can be considered to be representative of the benefit due to the fibre addition and its contribution to the energy necessary for the first cracking development. From literature (FIB, 2009) the total fracture energy for fibre-reinforced concrete has been seen to be three orders of magnitude greater than the corresponding plain mix.

Having observed that the tail of the stress-separation curve considerably extends with the addition of fibre (Bazânt, 2002), i.e more energy is now dissipated following the initial crack opening, it is reasonable to assume that roughly a 10 times greater energy is necessary only for the initial crack opening.

The total fracture energy for the system is calculated as it follows and its estimate appears to be two orders of magnitude greater than the corresponding plain mix.

Table 6.11 – Improvement in total fracture energy

N / m


N / m

GF,FRC_1.0 / GF,PC









Despite the estimate is lower than the values presented in literature, it is also necessary to mention that due to the material variability, possible testing errors and model approximations, it is difficult to verify the model proposed.

7 /                      Conclusions

The project results have shown that, limited to the experimental values collected, the fibre addition can improve the tensile capacity of concrete of approximately 20% in tension (1% fibre content) and 10% in compression, similarly to the ones observed by Kakooei et al. (2012) and Choi and Yuan (2005). These results are yet limited to a 1% polypropylene fibre content.

On the other hand, the results’ comparison of the standard splitting test and the designed direct tensile test have shown a large discrepancy from the tensile capacity proposed in the Eurocodes (2.9 MPa). Further research should, therefore, be undertaken to understand whether the 20% difference from the EC value is due to an inappropriate direct testing method, or to the inaccuracy of the splitting test.

It is also established that the splitting test in the determination of the FRC behaviour following the maximum stress is not adequate due to the ductility of the material (Denneman, 2011). A new standardised test should then be introduced on the line of the direct test proposed in the project.

Conclusively, despite the fact that the research is currently focused on the mechanics of the cracking development and the prediction of material behaviour, it is important to notice that the data previously analysed shows a significant decrease in the results range. The results gathered for the FRC specimens are, in fact, characterised by a coefficient of variation between 4.1 and 4.8% in contrast to the 8.7 and 18.4% of the corresponding plain mixture, as it was mentioned by Carmona et al (2013) and Hasewander (1995).

The power of this observation should not be underestimated because, whether the material behaviour is modelled from theory or from extensive experimental investigation, the factor that the stress response of concrete can be accurately predicted can represent a considerable economic and environmental effect for the construction sector. The current 1.5 partial safety factor for concrete could potentially be reduced, and the concrete tensile capacity could also be accounted for in the design process allowing for the reduction of materials and costs.

7.1 /      Recommendations for further studies

In its entirety, the project can be considered to have successfully presented a reliable set of data and theoretical basis. Nonetheless, improvements can still be made both for theoretical and experimental processes.

A possible theoretical improvement could be made within the analysis and determination of the fracture energy. Such revision would require advanced knowledge of fracture mechanism to understand and apply alternative theoretical models. An alternative experiment configuration, which focuses on the test of a notched prism as proposed by Zhang et al (2001) and Cunha et al. (2011), should also be considered for the gathering of data describing the material’s behaviour at fracture.

Concerning the experimental side of the project, modification can be applied to:

  • size of the specimen’s range,
  • number of mixture with different water/cement ratio as well as fibre content,
  • test set-up,
  • specimen geometry ,
  • typology of testing equipment,
  • resolution and range of the gauges,

An ideal evolution of this experiment capable of producing statistically significant results of the material strength should involve a minimum of 20 samples for each mixture and specimen typology (BRE, 1988). It should be carried out with a testing machine equipped with hydraulic wedge grips to eliminate the data corruption due to ends’ slipping. So lastly, the clamping mechanism should be designed to restrain any possible rotation due to non-uniform stress distribution following the first cracking development.

In the case these amendments were not feasible, the modification of the specimen geometry should be considered as the first possible experiment variation.

The current geometry does, in fact, define an area of stress concentration that causes an unwanted failure where the necking opens up to the end sections (Appendix E). As the gauges range over the middle of the necking, every failure outside this range makes the readings after peak stress unusable. A possible alternative to the current geometry is a simple 100 x 100 x 400mm prism with a 5mm notch 25mm deep, similar to the layout proposed by Zhang et al (2001) and Cunha et al. (2011), with 4 gauges (one per face) bridging the notch over a minimum range of three times the maximum coarse aggregate diameter.

By exploiting the localised stress concentration, this set-up should guarantee failure at the notch and thus suitable data for future studies.


  1. Aiello, M. A. and Ombres, L. (2002) Structural Performances of Concrete Beams with Hybrid (Fiber-Reinforced Polymer-Steel) Reinforcements, Journal of Materials in Civil Engineering (May 2002), pp. 133 – 140
  2. Alani, A. M. and Aboutalebi, M. (2013) Mechanical properties of fibre reinforced concrete – a comparative experimental study, International Journal of Civil, Architectural Science and Engineering, 7, pp. 111 – 116
  3. Avram, C., Bob, C., Friedrich, R. and Stoian V. (1993) Numerical Analysis of Reinforced Concrete Structures. Amsterdam: Elsevier
  4. Barragan, B. E., Gettu, R., Martin, M.A. and Zerbino, R.L. (2003) Uniaxial tension test for steel reinforced concrete – a parametric study, Cement & Concrete Composites, 23, pp. 767 – 777
  5. Bazant, Z. P. (2002) Concrete fracture models: testing and practice, Engineering Fracture Mechanics, 69, pp. 165 – 205
  6. Bhatt, P., MacGinley, T. and Choo, B. (2005) Reinforced Concrete: Design Theory and Examples. 3rd ed. Oxon: Taylor & Francis
  7. Bischoff, P. H. (2003) Tension Stiffening and Cracking of Steel Fiber-Reinforced Concrete, Journal of Materials in Civil Engineering (March/April 2003), pp. 174 – 182
  8. British Standards Institute (2009) BS EN 12390-1:2009: Testing hardened concrete – Part 1: Shape, dimensions and other requirements for specimens and moulds. London: British Standard Institute
  9. British Standards Institute (2009) BS EN 12390-2:2009: Testing hardened concrete – Part 2: Making and curing specimens for strength tests. London: British Standard Institute
  10. British Standards Institute (2009) BS EN 12390-3:2009: Testing hardened concrete – Part 3: Compressive strength of test specimens. London: British Standard Institute
  11. British Standards Institute (2009) BS EN 12390-5:2009: Testing hardened concrete – Part 5: Flexural strength of test specimens. London: British Standard Institute
  12. British Standards Institute (2009) BS EN 12390-6:2009: Testing hardened concrete – Part 6: Tensile splitting strength of test specimens. London: British Standard Institute
  13. British Standards Institute (2009) BS EN 1992-1-1:2004:Eurocode 2: Design of concrete structures – Part 1-1: General rules and rules for buildings. London: British Standard Institute
  14. Carmona, S., Aguado, A. and Molins, C. (2013) Characterization of the properties of steel fiber reinforced concrete by means of generalized Barcelona test, Construction and Building Materials, 48, pp. 592 – 600
  15. Choi, Y. and Yuan, R. L. (2005) Experimental relationship between splitting tensile strength and compressive strength of GFRC and PFRC, Cement and Concrete Research, 35, pp. 1587 – 1591
  16. Cunha, V. M. C. F., Barros, J. A. O. and Sena-Cruz, J. M. (2012) A finite element model with discrete embedded elements of fibre reinforced composites, Computers and Structures, 94-95, pp. 22 -33
  17. Dennema, E., Kearsley, E. P. and Visser, A.T. (2011) Splitting tensile test for fibre reinforced concrete, Materials and Structures, 44, pp. 1441 – 1449
  18. Grassl, P. and Jirasek, M. (2010) Meso-scale approach to modelling the fracture process zone of concrete subjected to uniaxial tension, International Journal of Solids and Structures, 47, pp. 957 – 968
  19. Hannant, D.J. (1978) Fibre Cements and Fibre Concretes. Chichester: John Wiley and Sons
  20. Hasalwander, B., Jonas, W. and Riech, H. (1995) Material equations for steel fibre reinforced concrete members, Nuclear Engineering and Design, 156, pp. 235 – 248
  21. Hassan, A. M. T., Jones, S. W. and Mahmud, G. H. (2012) Experimental test methods to determine the uniaxial tensile and compressive behaviour of ultra high performance fibre reinforced concrete (UHPFRC), Construction and Building Materials, 37 (December 2012), pp. 874 – 882
  22. Hull, D. (1981) An Introduction to Composite Material. Cambridge: Press Syndicate of the University of Cambridge
  23. Jenq, Y. S. and Shah, S. P. (1986) Crack Propagation in Fiber-Reinforced Concrete, Journal of Structural Engineer, 112, pp. 19 – 34
  24. Kakooei, S., Akil, H. M., Jamshidi, M. and Rouhi, J. (2012) The effect of polypropylene fibers on the properties of reinforced concrete structures, Construction and Building Materials, 27 (1), pp. 73 – 77
  25. Li, F. and Li, Z. (2001) Continuum damage mechanics based modelling of fiber reinforced concrete in tension, International Journal of Solids and Structures, 38 (5), pp. 777 – 793
  26. Owen, D. R. J., Holbeche, J. and Zienkiewicz, O. C. (1968) Elastic-plastic analysis of fibre-reinforced material, Fibre Science and Technology, 1, pp. 185 – 207
  27. Paramasivam, P., Curiskis, J. I. and Valliappan, S. (1984) Micromechanics analysis of fibre reinforced cement composites, Fibre Science and Technology, 20, pp. 99 – 120
  28. Pedziwiatr, J. (2009) The influence of the bond between concrete and reinforcement on tension stiffening effect, Magazine of Concrete Research, 61 (6), pp. 437 – 443
  29. Peng, X. and Meyer, C. (2000) A continuum damage mechanics model for reinforced concrete reinforced with randomly distributed short fibers, Computer & Structures, 78 (4), pp. 505 – 515
  30. Reddiar, M. K. M. (2009) Stress-strain model of unconfined and confined concrete and stress-block parameters, M.A. Thesis, Pondicherry Engineering College, Puducherry, India.
  31. Rots, J. G. and de Borst, R. (1989) Analysis of concrete fracture in “direct” tension, International Journal Solids Structures, 25 (12), pp. 1381 – 1394
  32. San-Jose, J. T., Vegas, I. and Meyer F. (2006) Structural analysis of FRC reinforced polymer concrete material, Construction and building materials, 20, pp. 971 – 981
  33. The Concrete Centre (2008) Properties of Concrete for use in Eurocode 2. Camberley: The Concrete Centre
  34. The International Federation of Structural Concrete (2009) Structural Concrete Textbook on Behaviour, Design and Performance , Volume 1. 2nd ed. Lausanne: fib
  35. Uomoto, T., Mutsuyoshi, H., Katsuki, F. and Misra, S. (2002) Use of fiber reinforced polymers composites as reinforcing material for concrete, Journal of Materials in Civil Engineering (May/June 2002), pp. 191- 209
  36. van Mier, J.G.M. and van Vliet, M.R.A. (2002) Uniaxial tension test for the determination of fracture parameters of concrete: state of the art, Engineering Fracture Mechanics, 69, pp. 235 – 247
  37. van Vliet, M.R.A., and van Mier, J.G.M. (2000) Experimental investigation of size effect in concrete and sandstone under uniaxial tension, Engineering Fracture Mechanics, 65, pp. 165 – 188
  38. Zhang, J., Stang, H. and Li, V. C. (2001) Crack bridging model for fiber reinforced concrete under fatigue tension, International Journal of Fatigue, 23 (8), pp. 655 – 670
  39. Bhatt, P., Macginley, J. and Choo, B. S. (2006) Reinforce Concrete. Design Theory and Examples. 3rd ed. Oxon: Taylor & Francis
  40. Mosley, W. H. Bungey, J. H. and Hulse, R. (2007) Reinforced Concrete Design to Eurocode 2. 6th ed. Basingstoke: Palgrave MacMillan


Appendix A – Table 3.1 / Strength and deformation characteristics for concrete (BS EN 1992-1-1:2004)

Appendix B – Mould construction drawings – Top view (not at scale)

Appendix C – Mould construction drawings – Side view and section BB (not at scale)

Appendix D – Mould construction drawings – Section AA (not at scale)

Appendix E – FRC_1.0 Fracture out of gauges’ range


Appendix F – Plain concrete brittle fracture

Appendix G – FRC_1.0 Fracture within gauges’ range

Cite This Work

To export a reference to this article please select a referencing stye below:

Reference Copied to Clipboard.
Reference Copied to Clipboard.
Reference Copied to Clipboard.
Reference Copied to Clipboard.
Reference Copied to Clipboard.
Reference Copied to Clipboard.
Reference Copied to Clipboard.

Related Services

View all

DMCA / Removal Request

If you are the original writer of this dissertation and no longer wish to have your work published on the website then please: