Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access March 2, 2023

Development of an adaptive coaxial concrete rheometer and rheological characterisation of fresh concrete

  • Sebastian Josch , Steffen Jesinghausen EMAIL logo and Hans-Joachim Schmid
From the journal Applied Rheology

Abstract

The accessibility to rheological parameters for concrete is becoming more and more relevant. This is mainly related to the constantly emerging challenges, such as not only the development of high-strength concretes is progressing very fast but also the simulation of the flow behaviour is of high importance. The main problem, however, is that the rheological characterisation of fresh concrete is not possible via commercial rheometers. The so-called concrete rheometers provide valuable relative values for comparing different concretes, but they cannot measure absolute values. Therefore, we developed an adaptive coaxial concrete rheometer (ACCR) that allows the measurement of fresh concrete with particles up to d max = 5.5 mm . The comparison of the ACCR with a commercial rheometer showed very good agreement for selected test materials (Newtonian fluid, shear thinning fluid, suspension, and yield stress fluid), so that self-compacting concrete was subsequently measured. Since these measurements showed a very high reproducibility, the rheological properties of the fresh concrete could be determined with high accuracy. The common flow models (Bingham (B), Herschel–Bulkley, modified Bingham (MB) models) were also tested for their applicability, with the Bingham and the modified Bingham model proving to be the best suitable ones.

1 Introduction

In the last decades, the development of new concretes gained more and more importance, which is mainly related to emerging challenges such as high-performance concretes, 3D printable concretes, or predictive flow simulations. To achieve this and superior concrete performance, the profound understanding of the rheology of fresh concrete is an essential requirement and therefore of high significance [1,2,3,4].

However, determining the exact rheological properties of fresh concrete is still very difficult today. This is mainly due to the fact that the measuring gap of commercial rheometers is too small for big aggregates. Therefore, the so-called “concrete rheometers” were developed, which are mostly similar to a coaxial system [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20], but there are also some that work on the principle of a parallel-plate system [21,22,23,24] or characterise the concrete via spreading tests [25]. In many cases, these rheometers are research projects and thus not commercially available. However, there are some commercial “concrete rheometers” offered by Schleibinger with the Viskomat XL, Viskomat NT, and eBT-V models. For determining the pumpability of concrete, Schleibinger also offers a very practice-oriented measuring device called SLIPER, which, however, cannot be compared with a rheometer purely in terms of principle.

All these devices are usually referred to as “concrete rheometers,” which, with a few exceptions, is due to the fact that the characteristic measured values are rotational speed and torque (analogous to common rheometers). In most cases, the Reiner–Riwlin equation [26], which is based on the Bingham model, is used to determine the flow properties, regardless of the actual flow behaviour. Therefore, even if the designation “rheometer” implies that absolute rheological parameters can be determined, this is only possible with difficulty or not at all possible with almost all “concrete rheometers” developed in the past [18]. This is due to the measuring geometries, which do not provide a mandatory analytical describable shear field (e.g. vane-like geometries, paddles, or spheres). This inevitably also means that different concrete rheometers are not comparable.

To counteract this problem at least for selected concrete rheometers, round-robin tests were carried out in 2,000 at LCPC in France and in 2003 at Master Builders Solutions in the USA [14,27]. The aim was to develop correlations between the rheometers by measuring identical concrete mixes, so that in future it would be possible to compare different concrete mixes measured in different concrete rheometers. The rheometers used were the BTRheom [21,22], the IBB [11], the Two-point rheometer [6], the ConTec BML [8,9], and the CEMAGREF-IMG [12]. Although different rheological properties, such as yield stress and plastic viscosity, were analysed, it was not possible to find a general solution for correlating these devices. Therefore, different suggestions for future investigations were raised, e.g. the establishment of standardised measurement methods for concrete and the use of suitable reference materials [14,27]. One of the goals pursued with these suggestions was the development of better rheometers, with which absolute values of concrete can be determined and used in simulations [27]. Since such rheometers have not been developed in the recent years, we developed an innovative concrete rheometer, the adaptive coaxial concrete rheometer (ACCR), to tackle the aforementioned challenges, which were met with great success.

2 Theory

2.1 Rheometry of concentric coaxial cylinder systems

The measuring gap of concentric coaxial cylinder systems is formed by two cylinders with the same symmetry axis, the inner cylinder (bob, radius R i , length L ) and the outer cylinder (cup, radius R o ). Although standardised systems have a rather small ratio δ = R o / R i , typically with values of δ 1.2 , the increasing cylindrical plane in the measuring gap cannot be neglected, leading to a dependency of the shear stress τ on the radius r , as shown in equation (1).

(1) τ ( r ) = M 2 π r 2 L .

For Newtonian fluids, the shear rate can directly be expressed by equation (2), where ω i and ω o are the angular velocity at the inner and outer cylinder, respectively, and β = ( R i / R o ) ² . Since most commercial rheometers and also the ACCR are searle type cylinder systems, meaning only the inner cylinder is rotating, and thus ω o = 0 , the viscosity can be calculated by equation (3) [28,29].

(2) γ ̇ ( r ) = 2 R i 2 ( ω i ω o ) r 2 ( 1 β ) ,

(3) η = γ ̇ ( r ) τ ( r ) = M ( 1 β ) 4 π ω i R i 2 L .

2.2 Correction of unknown flow behaviour

Since equations (1)–(3) are only valid for Newtonian fluids, different correction methods have been developed for unknown flow behaviour.

A rather simple method is the Schümmer correction, which uses the representative shear rate γ ̇ rep and the representative shear stress τ rep . This method is based on the work of Schümmer [30], who was able to define a representative radius r rep in the measuring gap at which the shear rate and the shear stress are equal to those of a Newtonian fluid. Giesekus and Langer [31] extended Schümmer’s approach and developed equations for further calculation of γ ̇ rep and τ rep , equations (4) and (5).

(4) γ ̇ rep = 1 + β 2 γ ̇ i = ε γ ̇ i ,

(5) τ rep = 1 + β 2 τ i = ε τ i .

However, this correction method is mainly accurate for fluids with certain flow behaviour, so that a more general and accurate correction method developed by Krieger and Elrod should be applied, if possible [28,31].

The Krieger–Elrod correction is based on the general formulation of the shear rate, which was developed in a power series, equation (6). The parameter s represents the slope of the curve obtained by pointwise differentiation in the ln ( ω ) ln ( τ ) -diagram.

(6) γ ̇ ( τ ) = ω i ln δ 1 + s ln δ + 1 3 ( s ln δ ) 2 + 1 3 ln δ 2 d s dln τ .

The accuracy of this correction depends on the factor s ln δ . For s ln δ < 0.1 , it is sufficient to only consider the first two terms, since the resulting error is below 1 % . If s ln δ > 0.1 , the remaining terms of the series must also be included in the correction. As the derivation for the development of the series assumes that γ ̇ ( τ = 0 ) = 0 , this correction is not suitable for materials with yield stress [28,32].

2.3 Flow behaviour of fresh concrete

The most established model for describing the flow behaviour of cement-based suspensions is the Bingham model [14,27,33,34,35]. Bingham fluids are characterised by a yield stress with subsequent Newtonian flow behaviour and can be described mathematically by equation (7) [34].

(7) τ = τ y + η pl γ . ̇

The slope of the straight line in the τ / γ ̇ -diagram is referred to as the plastic viscosity η pl in order to distinguish it from the common Newtonian viscosity.

However, nonlinear flow behaviour of especially self-compacting concretes has been observed in numerous investigations [36,37,38,39,40] and the application of the Bingham model led to physically impossible negative yield stress values [36,37,41]. Therefore, the so-called Herschel–Bulkley model was increasingly used to describe the rheological properties, equation (8).

(8) τ = τ y + k γ ̇ n .

This model replaces the original linear Newtonian relationship between τ and γ ̇ with the approach for power law fluids developed by Ostwald and de Waele. The Parameter k is the consistency of a fluid and n describes the flow exponent, which is n > 1 for shear thickening and n < 1 for shear thinning fluids.

Even though the Herschel–Bulkley model initially appears to be suitable for characterizing nonlinear behaviour, some further investigations have concluded that the yield stress values determined are strongly dependent on the flow exponent. The comparison of this model with other models has consistently shown a yield stress that is too low for shear-thinning materials [42] and too high for shear-thickening materials [43,44]. To overcome this problem, the Bingham model was extended by the quadratic term c γ ̇ ² , from which the so-called modified Bingham model was derived, equation (9) [42]. The parameter c is an empirical coefficient for the quadratic term and a more precise definition of this value is missing in the literature.

(9) τ = τ y + η pl γ ̇ + c γ ̇ 2 .

When applying the modified Bingham model, it has been observed that the yield stress determined via this model always lies between the yield stress determined via Herschel–Bulkley and Bingham and therefore, seems to be a more accurate estimation [42]. In addition, according to Feys et al., the yield stress is determined independent of the degree of shear thickening, which can further support the accuracy of this model [40].

3 Development and construction of the ACCR

In principle, rheometry offers a wide range of rheometer types, but the choice of a suitable measuring system for determining absolute flow properties as accurately as possible is limited. The most common rheometers are parallel-plate systems, coaxial systems, and capillary rheometers. However, for a highly concentrated, coarse-grained suspension such as concrete, which, depending on its composition, tends to sediment, parallel-plate rheometers and capillary rheometers are not suitable. For a parallel-plate system, even a slight sedimentation would cause the upper plate to apparently slip and the sheer dimension of the measuring gap, due to the particle size, would most probably lead to a flow-out of the concrete during or even before the measurement itself. In a capillary rheometer, there would be no possibility to prevent wall slip or to perform time dependent or yield stress measurements. Thus, only a coaxial rheometer is suitable. Figure 1 shows the schematic structure of the ACCR.

Figure 1 
               Schematic diagram of the ACCR with corresponding dimensions (in mm); DU = drive unit and CU = control unit.
Figure 1

Schematic diagram of the ACCR with corresponding dimensions (in mm); DU = drive unit and CU = control unit.

3.1 Basic dimensions

When dimensioning a coaxial system for suspensions, the maximum particle size is the first thing to be determined. The maximum measurable particle size d max is limited by the size of the measuring gap a of the rheometer.

(10) d max a / 5 .

If this relationship is neglected, measurement errors due to friction and the influence of the particles on the flow field cannot be excluded [29,45].

Second, according to ISO 3219, the radius ratio δ should be chosen as follows [46]:

(11) δ = R o / R i 1 . 2 .

To maintain this ratio is mandatory to achieve an uniform shear rate and stress distribution among the gap [29,46]. For manufacturing and handling reasons, a maximum measurable particle size of d max = 5.5 mm was chosen and accordingly, the radii of the ACCR are R o = 165 mm and R i = 137.5 mm , resulting in a gap width of a = 27.5 mm .

With regard to the rheometer dimensions, the ISO 3219 further specifies that the following relationship must apply to the ratio of the inner cylinder length L and the inner radius R i to neglect edge effects [46].

(12) L / R i 3 .

However, since there is no rheological justification to exact this ratio in ISO 3219 and a corresponding length of L = 412.5 mm has not been possible for manufacturing reasons, the cylinder length of the ACCR is L = 300 mm , which corresponds to a ratio of L / R i = 2.18 . Considering all the components of the rheometer, this results in a measuring volume of V = 13.87 l .

3.2 Rotating bottom plate

In all coaxial systems, front face influences occur on the upper and lower front faces of the inner cylinder, which, without correction or other measures, lead to erroneous results. The upper front face influence can be eliminated by filling the measuring gap to the upper edge of the inner cylinder only. Since the upper face and other components of the rheometer above it are only in contact with air, the influence is in general negligibly small.

The situation is different on the lower front face of the inner cylinder. An additional parallel-plate system with a theoretical measuring gap L' would form between the lower front face of the inner cylinder and the bottom of the outer cylinder, which can have a considerable influence on the measuring results. For this reason, the ISO 3219 also specifies the value for the gap with L / R i 1 , with the aim of keeping the influence of the lower front face negligibly small, since the shear stress decreases linearly with the distance for a parallel-plate system.

In the case of the ACCR, a design with this specification is not expedient because of handling issues and a significantly higher, impractical sample volume. Instead, the effect was inhibited with a rotating bottom plate below the inner cylinder moving at the same speed, Figure 2. Thus, the material between the rotating bottom plate and the inner cylinder is not sheared at all and therefore the torque is not increased. In order to ensure that the rotational speed of the inner cylinder and the bottom plate match as closely as possible, both components are driven by identical motors and gears.

Figure 2 
                  Sectional view to illustrate the emptying of rheometer by moving the bottom unit, closed on the left, and opened on the right.
Figure 2

Sectional view to illustrate the emptying of rheometer by moving the bottom unit, closed on the left, and opened on the right.

Furthermore, the whole bottom part of the rheometer is, unlike commercial rheometers, detachable to empty and clean the rheometer after measurement, Figure 2.

3.3 Adaptive measuring profiles with integrated cooling

When measuring suspensions, apparent wall slip is a common problem, so that wall adhesion can no longer be guaranteed and the use of roughened or profiled measuring surfaces is usually practiced [47,48,49,50,51]. To allow a slip inhibition for the ACCR, an adaptive approach with several, exchangeable measuring profiles was implemented. The measuring surface of the coaxial system therefore consists of six individual measuring profiles at the inner and outer cylinders, which, when lined up together, form a circular ring. The standard measuring profiles for materials prone to wall slip are made of stainless steel with a length of 300 mm and have square like longitudinal grooves with a dimension of 5.5 mm , which corresponds to the maximum measurable particle size, Figure 3. For non-reactive, wall-adhering fluids, sandblasted profiles made of aluminium are also available, which were mainly used for reference measurements.

Figure 3 
                  CAD-construction of the cooling profiles (outer cylinder); profiles of inner cylinder are similar, but without cooling-channel.
Figure 3

CAD-construction of the cooling profiles (outer cylinder); profiles of inner cylinder are similar, but without cooling-channel.

In addition, the standard profiles (stainless steel) for the outer cylinder cool the measured fluid in the gap. A copper coil is integrated in each profile to be fed with a temperature-controlled cooling fluid, Figure 3. The cooling fluid is fed in parallel through the profiles, ensuring a homogeneous temperature. Alternatively, a cooling coil wrapped around the outer cylinder can be used for cooling, which is sufficient for non-exothermally reactive materials and was used in the context of comparative measurements together with the smooth aluminium profiles, see Section 4, and partly during the first series of tests, see Section 5.1. It is important to mention that a cooling device is always required because of the self-heating caused by the shearing. In commercial rheometers, cooling can be realized with little effort, since the sample volume is very small. However, in the ACCR with a much larger sample volume, heat is generated not only by the shear process, but also by the reaction of the fresh concrete itself. Therefore, the cooling of the concrete is generally a very difficult challenge, but one that has been mastered very well with the cooling profiles when the temperature values achieved are considered, see Section 5.

3.4 Overall design and key parameters

Figure 4 shows the setup of the whole rheometer. The inner cylinder and the bottom plate are each driven by an AC asynchronous motor from Lenze (type MXXMA), which has a nominal rotational speed of N = 3,480 1 / min (controlled via an incremental encoder) and a nominal torque of M = 10.9 Nm . In order to be able to guarantee low rotational speeds with a high degree of running smoothness, a Bonfiglioli gearbox (type 32 2 P 52.4 HS V5) with a ratio of i = 52.4 is installed between the motor and the corresponding drive shaft. This combination allows shear rates in the range of γ ̇ = 1 45 s 1 . The torque is measured between the gear unit and the drive shaft of the inner cylinder via a sensor from burster (type 8661). The sensor has a measuring range of M = 0 50 Nm with a relative linearity deviation of 0.05 % in relation to the end value of the measuring range. An encoder disk with 2,000 increments for speed and angle of rotation measurement is also integrated in the sensor. Alternatively, a smaller torque sensor from Lorenz can be installed for fluids with significantly lower viscosity. It has a measuring range of M d = 0 5 Nm with a relative linearity deviation of 0.25 % in relation to the measuring range end value. However, the large sensor is used as a standard, the small sensor was used only for the initial comparative measurements of the rheometer, see Section 4. The rheometer is controlled via a self-written LabVIEW program.

Figure 4 
                  ACCR setup.
Figure 4

ACCR setup.

3.5 Error analysis

To assess the validity of the results presented in this study, it is necessary to estimate the measurement accuracy of the ACCR. Therefore, all possible sources of error were identified and evaluated. The relevant errors for the rheological parameters shear rate and shear stress were characterised in more detail in a Gaussian error propagation, which lead to relative inaccuracies of u ̅ rel , γ ̇ = 2.18 2.73 % and u ̅ rel , τ = 3.54 % . Further information on the calculation is given in the supplementary material S1.

4 Comparative measurements

For evaluation purposes, the ACCR was tested first with simple, controllable fluids and compared to a commercial Anton Paar MCR 501 using a CC27 standard system ( R i = 13.333 mm , R o = 14.462 mm , and L = 40.026 mm ). A highly concentrated sugar solution (Newtonian flow behaviour), a cream bath (shear thinning properties), a suspension (shear thinning properties and low yield stress), and mayonnaise (shear thinning properties and high yield stress) were chosen as test fluids to comprise a great range of rheological behaviour.

For the sugar solution and the cream bath, a positive shear rate ramp (PSR) starting with γ ̇ i , s = 1 s 1 and raising it to γ ̇ i , s = 45 s 1 was applied. For the suspension and the mayonnaise, a negative shear rate ramp (NSR) lowering γ ̇ i , s from 45 to 1 s 1 . This is recommended for yield stress fluids with time-dependent structural build-up like the later investigated fresh concrete suspensions [52]. Thus, the (reversible) structure formation is inhibited till the end of the measurement, enhancing the reproducibility. This would not have been necessary for the fluids investigated at this point but was done to facilitate a comparison to the concrete measurements. In addition, a pre-shearing step ( γ ̇ i , s = 45 s 1 for 30 s ) was performed on the suspension and mayonnaise after the materials were filled into the gap to homogenise the fluids, and a direct pre-shearing step ( γ ̇ i , s = 45 s 1 for 30 s ) was performed before each individual measurement to inhibit structural build-up over time. Each measurement point was held for 10 s and averaged over the last 5  s, since the torque curves appeared to be constant in this time range. There were different measuring series which always comprised three individual measurements without changing the specimen (Table 1). Before every series, the friction of the system was measured and considered during evaluation. Following the above protocols, the mean standard deviation between measurements was below 1 % (sugar solution, cream bath, and suspension) and 4 % (mayonnaise) and therefore, too small to be shown in Figure 5. To calculate the deviation between the commercial system and the ACCR, the results were fitted using an adequate fit and a mean deviation over the measurement range was calculated. For all measurements the ACCR was equipped with the smooth aluminium surface and the appropriate torque sensor, see Section 3.4.

Table 1

Analysed measurement samples with associated measurement parameters and relative deviations

Fluid MS Temperature (°C) Correction Fit R 2 Deviation (%)
Sugar solution 3 20.65 ± 0.15 Schümmer 1 1 < 4.8
Cream bath 2 × 3 19.75 ± 0.05 Krieger–Elrod 5 1 < 4.0
20.30 ± 0.10
Suspension 3 23.14 ± 0.26 Krieger–Elrod 6 1 < 4.7
Mayonnaise 3 22.83 ± 0.14 Schümmer 1 < 12.7

The fit-numbers represent the used polynomial degree, except for the mayonnaise, where the Herschel-Bulkley-model was used.

Figure 5 
               Flow curves measured with the ACCR and the MCR 501; (a) sugar solution with 
                     
                        
                        
                           
                              
                                 ω
                              
                              
                                 m
                              
                           
                           =
                           0.74
                        
                        {\omega }_{{\rm{m}}}=0.74
                     
                  ; (b) cream bath; (c) suspension consisting of cream bath and polyamide particles with 
                     
                        
                        
                           
                              
                                 ϕ
                              
                              
                                 PA
                              
                           
                           =
                           0.3
                        
                        {\phi }_{{\rm{PA}}}=0.3
                     
                  ; and (d) mayonnaise (mean values of the associated individual measurements are marked as X.1, X.2, and X.3).
Figure 5

Flow curves measured with the ACCR and the MCR 501; (a) sugar solution with ω m = 0.74 ; (b) cream bath; (c) suspension consisting of cream bath and polyamide particles with ϕ PA = 0.3 ; and (d) mayonnaise (mean values of the associated individual measurements are marked as X.1, X.2, and X.3).

4.1 Newtonian fluid–sugar solution

The aqueous sugar solution contained a sugar mass fraction of ω m = 0.736 and a citric acid mass fraction of ω m = 0.0025 to inhibit crystallisation. It can be clearly seen from Figure 5(a) that the results of both the rheometers agree well and clearly show a Newtonian behaviour at a viscosity of η MCR = 1.992 ± 0.002 Pa s and η ACCR = 2.088 ± 0.084 Pa s . The values of the ACCR tend to be slightly higher than those of the MCR 501.

4.2 Shear thinning fluid–cream bath

For the cream bath measurements, an additional pre-shearing step for 30 s at γ ̇ i , s = 5 s 1 with a subsequent pause of 10 s was added to the schedule to homogenise the fluid in the measuring gap. Since the flow clearly shows a shear thinning behaviour, Figure 5(b), a correction according to Krieger–Elrod was applied. Again, the results agree quite well and the values of the ACCR tend to be slightly higher than those of the MCR 501 in the medium shear rate range ( γ ̇ 10 40 s 1 ), although this behaviour seems to reverse in the higher shear rate range ( γ ̇ > 40 s 1 ).

4.3 Shear thinning suspension–cream bath with polyamide particles

To evaluate the ACCR regarding suspensions, it had to be first decided whether the exact same suspension or a similar suspension with the same particle diameter to gap ratio should be used, while both the variants promise different advantages. Since preliminary measurements showed that there is nearly no tendency for wall slip (plate-plate and coaxial systems), it was decided to use the exact same suspension (particle volume fraction ϕ PA = 0.3 , x 50,3 = 60.87 ± 0.02 μ m , cream bath as matrix fluid) as it was easier to prepare, knowing that it might behave different in both the rheometers because of the particle to gap ratio, even though the recommended particle to gap ratio smaller than 1 / 5 was met, equation (10). The mean particle diameter was chosen according to the commercial system. With a value of ρ PA = 1,010 kg / m 3 , the polyamide particles have a density close to that of the cream bath with ρ CB = 1,040 kg / m 3 and thus flotation effects can be neglected, at least during the time of measurement. Since a suspension also behaves non-Newtonian, the Krieger–Elrod correction was applied, Figure 5(c). Although materials with yield stress should not be corrected via the Krieger–Elrod method, this influence was considered negligible, since the uncorrected values indicated that a yield stress, if present at all, is close to zero. The measured values of both the rheometers also agree very well and the shear thinning behaviour of the suspension becomes evident. An extrapolation of the polynomials leads to yield stresses of τ y = 1.4 Pa (ACCR) and τ y = 1.6 Pa (MCR), although these values depend on the fit and have to be considered with care. The slight tendency of the ACCR values lying above the values of the MCR is also obvious here. However, since the preliminary investigation has shown a minute tendency to wall slip in the MCR, this deviation might be even smaller due to the less pronounced wall slip in the ACCR (smaller particle to gap ratio and sandblasted aluminium surface).

4.4 Yield stress fluid–mayonnaise

For mayonnaise, a decrease in the measured values was observed with the increase in the individual measurements for both the rheometers, so that the individual measurements (X.1–X.3) are presented instead of a single mean curve. The reason for this behaviour might be rising air bubbles during the measurements that had been entrapped in the filling process, structural changes in the mayonnaise due to drying processes, or even an influence of the different loading processes. The effect is stronger for the ACCR, which corresponds with the mentioned assumptions. Wall slip can be excluded as a reason since preliminary investigations showed no tendency. Due to the high yield stress, a correction via Krieger–Elrod is not possible, so that the Schümmer-correction was applied. However, the Herschel–Bulkley model could be fitted to the values for calculating the deviations and analysing the yield stresses, Figure 5(d). The results show that the values of the ACCR are below the values of the MCR. There are deviations in the range between 9.83 % (X.1) and 12.69 % (X.3). But again, a good agreement was achieved. For the yield stresses derived from the Herschel–Bulkley model, deviations between 6.56 % (X.1) and 3.21 % were measured.

4.5 Conclusion

The first three comparisons clearly show that the ACCR delivers consistent absolute rheological values within a close error margin to commercial measuring systems. On a closer look, a constant offset of about 5 % is present. The fourth measurement of the more complex mayonnaise shows slightly stronger deviations in the other directions. This can be accredited to the more complex loading process and the bigger volume of the ACCR compared to the MCR, which seems to enforce certain errors. Because of the general good agreement, no further calibration was carried out and the ACCR was used for the investigation of real concrete.

5 Rheological characterisation of fresh concrete

One of the most important reasons why a generally valid rheological classification of fresh concrete is much more difficult than for other materials is the fact that the term “concrete” does not refer to a clearly defined substance, but rather to an entire class of materials. Concrete consists of four components (cement, water, aggregates, and various additives) in different ratios. Even the mixing time and the time to the measurement are relevant, since it is a reactive material. Therefore, it is necessary to define a mixture as well as an accurate mixing and measuring procedure to allow for comparable measurements. The exact composition and handling are described in the supplementary material S2. All measurement results are approximated by representative values according to Schümmer. From this point, the following wordings are used to distinguish between different sets of measurements:

  • Individual measurements

  • Measurement series (consisting of multiple individual measurements)

  • Investigation series (consisting of multiple measurement series)

  • Series of tests (consisting of multiple investigation series)

5.1 First series of tests

During the first series of tests, the self compacting concrete (SCC) was investigated in two investigation series, which are presented and discussed below. The first investigation series consisted of measuring the SCC using originally smooth aluminium profiles, which were fitted with brass flat bars ( 4 mm thick and 6 mm wide) to prevent wall slip. The reaction occurring between aluminium and concrete with the formation of hydrogen [53] was considered negligible, since no visible wear could be detected when the profiles were cleaned. The second investigation series comprised hysteresis measurements using stainless steel profiles. During these first series of tests, the concrete was still measured with the alternative tempering system (cooling hose on the outside) and was filled into the rheometer via suitable funnels.

5.1.1 Standard measuring procedure – aluminium profiles

In the first investigation series, the SCC was measured in three measurement series, with one measurement series consisting of five individual measurements without changing the specimen. The results of all 15 individual measurements are shown in Figure 6. The mean temperature of the first measurement series was T m = 22.17 ± 0.27 ° C , for the second measurement series T m = 22.20 ± 0.11 ° C , and for the third measurement series T m = 19.77 ± 0.12 ° C , whereby the temperatures as expected increased slightly with the increase in the individual measurements for each measurement series. Regarding the reproducibility, the measured values, with the exception of the individual measurements 3 and 5 of the third measurement series, agree quite well. A deviation in the third measurement series can be explained by the lower temperature as an indication of a comparatively reduced reaction of the concrete and hence a less pronounced structural build-up, causing a lower shear stress compared to measurement series 1 and 2 .

Figure 6 
                     Representative values of SCC, measured in the standard process with modified aluminium profiles.
Figure 6

Representative values of SCC, measured in the standard process with modified aluminium profiles.

It can be seen that the shear stress increases significantly with the increase in the individual measurement number. This can be attributed to the reaction of the concrete and the structural build-up over time. This behaviour is particularly important in the context of the reproducibility of future measurements, as it highlights the relevance of adhering to a strict schedule when measuring concrete.

Furthermore, the flow behaviour of the concrete can be characterised in more detail. For low shear rates of approx. γ ̇ rep = 5 s 1 , the curves indicate a shear thinning behaviour, which changes to Newtonian flow behaviour in the medium shear rate range of approx. γ ̇ rep = 5 35 s 1 . If the shear rate increases above γ ̇ rep = 35 s 1 , dilatant behaviour seems to arise. These characteristics could be explained by the properties of highly concentrated suspensions. As the shear rate increases, the particles form sliding planes, which lead to a reduced resistance, i.e. shear thinning behaviour, and then converges to a steady state with Newtonian behaviour. If the relative velocities of the particles become too high at high shear rates, the impact probability of the particles and thus the energy dissipation increases, which could explain the dilatant flow behaviour. But the alleged dilatant behaviour at high shear rates may also be the result of too less pre-shearing and an already existing structural build-up at the beginning.

Furthermore, the shear thinning behaviour at low shear rates must be considered in more detail. During some measurements, a completely unsheared area was formed, probably due to a drop below the yield stress, resulting in a reduction in the measuring gap starting from the outer cylinder, Figure 7. In such a state, the shear rate is effectively higher than assumed and the influence of coarse-grained particles in the reduced gap is bigger, so that the measured shear stresses lose validity. However, this phenomenon occurred primarily in the later individual measurements, so that at least the first individual measurements are valid for classical rheological investigation. A correction of this block building phenomenon is almost impossible due to the large number of influencing parameters. Not only the yield stress of the fresh concrete, which depends on the time and the load history, is relevant, but also the parameters already mentioned, such as the particle-to-gap ratio, the particle shape, the orientation of the particles during shear, and the particle migration. It is also unclear to what extent the local concrete composition changes with this reduction in the measuring gap, since, at least apparently, an increased water content could be observed in the reduced gap. Furthermore, the occurrence of the gap reduction can only be detected by observing the measuring gap and not from the measurement values themselves.

Figure 7 
                     Schematic representation of the angular velocity and yield stress dependent gap reduction from 
                           
                              
                              
                                 
                                    
                                       R
                                    
                                    
                                       o
                                       ,
                                       0
                                    
                                 
                              
                              {R}_{{\rm{o}},0}
                           
                         to 
                           
                              
                              
                                 
                                    
                                       R
                                    
                                    
                                       o
                                    
                                 
                                 (
                                 ω
                                 ,
                                 
                                    
                                       τ
                                    
                                    
                                       y
                                    
                                 
                                 )
                              
                              {R}_{{\rm{o}}}(\omega \left,{\tau }_{y})
                           
                         in the presence of 
                           
                              
                              
                                 
                                    
                                       τ
                                    
                                    
                                       y
                                    
                                 
                                 >
                                 
                                    
                                       τ
                                    
                                    
                                       o
                                    
                                 
                              
                              {\tau }_{y}\gt {\tau }_{{\rm{o}}}
                           
                         from the sectional view with shear stress profile (left) and the plan view (right).
Figure 7

Schematic representation of the angular velocity and yield stress dependent gap reduction from R o , 0 to R o ( ω , τ y ) in the presence of τ y > τ o from the sectional view with shear stress profile (left) and the plan view (right).

A possible solution to this problem would be the increase in the geometry parameter β , since this minimises the difference between τ o and τ i , and the yield stress on the outer cylinder would only be undercut at very low shear rates, if at all. Although β can be increased in the ACCR by using thicker-walled measurement profiles, the particle-to-gap ratio would decrease and make further measurements doubtful in this regard. A scale up of the whole system would come with other problems like general manufacturing, measurement operability, or positioning accuracy.

In the following, the applicability of the common flow models (Bingham (B), Herschel-Bulkley (HB), modified Bingham model (MB)) for concrete, see Section 2.3, will be used for further characterisation. Since the values in the high shear rate range ( γ ̇ rep > 35 s 1 ) are considered problematic due to a lack of knowledge regarding sufficient pre-shearing, only measured values up to γ ̇ rep 35 s 1 are examined. The measured values in the low shear rate range ( γ ̇ rep < 5 s 1 ) are also discussed carefully due to the abovementioned gap reduction, which occurred only in the later individual measurements 4 and 5 . Therefore, only the first three individual measurements are considered in the following evaluation. In addition, since the third measurement series only matched the other measurement series to a limited extent due to the deviating temperature and strong deviation of the individual measurements 3 and 5 , only the mean values of the first two measurement series are included in the evaluation. Figure 8 shows the mean values of the considered individual measurements and the corresponding fitting parameters can be found in Table 2. The classical concrete models describe the flow behaviour of the SCC well, but there are some deviations in the low shear rate region ( γ ̇ rep < 5 s 1 ). No significant difference can be seen between the Bingham and the modified Bingham model, only the Herschel–Bulkley model seems to be able to somehow reproduce the shear thinning region for γ ̇ rep < 5 s 1 .

Figure 8 
                     Mean values of the individual measurements of SCC from the first two measurement series fitted with the Bingham model, the Herschel–Bulkley model, and the modified Bingham model.
Figure 8

Mean values of the individual measurements of SCC from the first two measurement series fitted with the Bingham model, the Herschel–Bulkley model, and the modified Bingham model.

Table 2

Model parameters of the Bingham, Herschel–Bulkley, and modified Bingham model for the fits of the mean values of the first three individual measurements of the first two measurement series

1st individual measurement 2nd individual measurement 3rd individual measurement
Bingham R 2 = 0.9990 R 2 = 0.9977 R 2 = 0.9961
τ y / Pa 72.24 114.57 158.97
η pl / Pa s 11.67 11.13 10.44
Herschel–Bulkley R 2 = 0.9991 R 2 = 0.9987 R 2 = 0.9986
τ y / Pa 69.78 98.85 135.20
k / Pa s n 13.05 16.36 18.96
n / 0.97 0.90 0.84
Mod. Bingham R 2 = 0.9990 R 2 = 0.9982 R 2 = 0.9975
τ y / Pa 73.69 108.40 149.24
η pl / Pa s 11.75 12.12 12.05
c / Pa s 2 0.002 0.027 0.045

Taking a closer look at Table 2, it can be seen that the yield stress increases significantly for all models with the increase in the individual measurement and is about twice as high for the third individual measurement. A reason for this might be the structural build-up. However, the yield stresses are of the same order of magnitude for all models with very low deviations.

The parameters describing the viscosity ( η pl and k ) are also of the same order of magnitude for all models in the first individual measurement, but deviate more strongly from one another for the further individual measurements. It is striking that for the Bingham and the modified Bingham model, the plastic viscosity η pl changes only slightly with the increase in the individual measurement, but the consistency k of the Herschel–Bulkley model increases strongly. While both the Bingham models thus indicate a flow behaviour in which only the yield stress increases with time, but the plastic viscosity remains almost constant, the Herschel–Bulkley model indicates a basic change in the flow behaviour. This is shown not only by the increase in consistency, but also by the flow index n , which by a steady decrease indicates an increasingly pronounced shear thinning behaviour with the increase in the individual measurement.

To this extent, it is not possible to clearly identify which of the models is best suited to describe the flow behaviour of the SCC based on the model parameters listed in Table 2 and Figure 8. If, due to the gap reduction, only the first individual measurement is considered, the Herschel–Bulkley model seems to be the best in terms of both error sum of squares and visual matching with the measured values. However, since the differences are minute and this evaluation is based on only two measurements, further investigation is necessary, see Section 5.2.

The reaction of the aluminium was assumed to be negligible in this investigation because the contact area of the aluminium to the concrete had been very small compared to the total sample material. Nevertheless, during this investigation series, it was found that the gas bubbles had risen. For this reason, all subsequent investigations have been carried out with stainless steel profiles.

5.1.2 Hysteresis measurements

In the context of this work, the flow behaviour of SCC is generally characterised by conducting an NSR starting at high shear rates and then decreasing the shear rate step by step. Nevertheless, investigation via hysteresis measurements can provide a deeper insight into the rheology of SCC. Therefore, a PSR starting at low shear rates and then increasing the shear rate was added to the standard measurement procedure, which is performed directly after the negative ramp. Figure 9 presents the results of the hysteresis measurements where two measurement series were performed with two individual measurements each. Due to the longer investigation period, only two measurement cycles were possible before stronger hardening became relevant. The mean temperature for the first measurement series was T m = 23.80 ± 0.15 ° C and for the second measurement series T m = 22.20 ± 0.3 ° C . As before, the temperatures increased slightly for each individual measurement. Again, maintaining the same temperature for all measurements is very difficult and the deviations in the results can be ascribed to the insufficient temperature control. Nevertheless, a qualitative analyzation is possible. First of all, it can be noted that the curve of the PSR shows a clearly different characteristic than the NSR. While the NSR at the beginning, i.e. at shear rates in the range of γ ̇ rep = 35 40 s 1 , indicates a dilatant behaviour, this cannot be recognised for the PSR at any time. This supports the assumption of insufficient pre-shearing already made in Section 5.1.1. Furthermore, in the range of low shear rates ( γ ̇ rep < 10 s 1 ), the PSR shows a slightly more pronounced shear thinning behaviour than the NSR. However, it is unclear whether this behaviour could also be due to the phenomenon of gap reduction and/or structure build-up, since the values of the PSR for the low shear rate range show much larger fluctuations than those of the NSR. An interesting effect becomes obvious analysing the NSR curves recorded immediately after the PSR curve, i.e. measurements 1.1 PSR → 1.2 NSR and 2.1 PSR → 2.2 NSR, as they lie below the values of the PSR curve and intersect them in both cases at a shear rate of about γ ̇ rep = 10 s 1 . This behaviour suggests that more structures were present in the earlier timed PSR measurements. This leads to two conclusions. There is a structure formation even if low shear is present and the structure destruction may be strongly dependent on the shear rate and the time.

Figure 9 
                     Representative values of the hysteresis measurements of SCC for performing two measurement series with two individual measurements each.
Figure 9

Representative values of the hysteresis measurements of SCC for performing two measurement series with two individual measurements each.

Therefore, the dilatant behaviour at the beginning of the NSR might occur due to insufficient pre-shearing since this behaviour could not be observed in the PSR. These results clearly emphasize how important a sufficient pre-shearing and the use of NSR is.

5.2 Second series of tests

The first series of tests have provided valuable results, especially with regard to the rheological properties of fresh concrete and the test procedure. However, many sources of error were also identified, which must be considered for a more reliable rheological characterisation of fresh concrete.

As the temperature was a major problem in both the measurement series of the previous tests with a deviation of 1 or 2 degrees, the directly cooled profiles, see Section 3.3, were used in the following second series of tests. Additionally, the pre-shearing step was intensified due to the observations made in Sections 5.1.1 and 5.1.2, Table S2.3 (supplementary material S2). Since the reaction with aluminium could not be excluded and even little gas bubbles emerged during earlier measurements, stainless steel profiles were used for further investigation, as already done for the hysteresis measurements. To reduce the influence of the filling process, a filling aid was introduced. The filling aid essentially consists of a large hopper in which the concrete is force-fed through an auger and filled into the rheometer gap via a hose. This makes it possible to fill the rheometer more reproducible than using the funnels, and additionally, shear can be applied in the measuring gap during the filling process, which not only serves to homogenise and destroy structures, but also contributes to deaeration. A preliminary investigation with the new setup showed that, despite being stored in vacuum, the cement mix appeared to have changed a little and no coherent measurements were possible. Therefore, it was necessary to adjust the concrete composition. To change the mixture as less as possible, a rheological thickener was added (Sika® Stabiliser-165, 0.165 % by weight of cement), which successfully reduced the spread flow to a value originally achieved i.e. approx. 65 cm , resulting in a similar behaviour as before. The mixing process was adapted to add the thickener. During the initial mixing of all components, the thickener was distributed over a large area while the mixing was paused for 30 s . This process was added to the mixing procedure, see supplementary material S2, Table S2.2, in the time 04 : 30 05 : 00 . However, the corresponding stage was not extended by 30 s , so that all other time steps remained identical. Further measurements were performed using this adjustment.

The results presented below are divided into three sections: reproducibility, rheological characterisation, and additional measurements. The additional measurements focus on sedimentation. Since the concrete already showed a high degree of consolidation in the individual measurements 4 and 5 in the previous series of tests, only the first two individual measurements were investigated from here on.

5.2.1 Reproducibility

All 15 measurement series conducted can be seen in Figure 10 (first individual measurements) and Figure 11 (second individual measurements). When comparing the two diagrams, it is obvious that the fluctuations during one measurement are more pronounced in the first individual measurements. Even if the pre-shearing has already been greatly increased, this, again, could be an indication that various accumulations of structures were still present. However, a clear cause for these curves could not be determined. Especially for the second individual measurements, it is noticeable that some individual measurements strongly deviate from the majority and were excluded from further investigation. These are 1.2, 9.2, 12.2, 14.2, and 15.2. These deviations coincide with the macroscopic concrete behaviour observed during the measurements, since either an increased temperature was detected (1) or the formation of shear bands (unsheared, circular ring-shaped areas moving with the direction of rotation) occurred (9, 12, 14, 15). Even though in a good agreement with the other individual measurements, measurements 2 and 6 were also excluded from consideration because of the same reasons. Figure 12(a) and (b) shows the remaining measurements and the good reproducibility. As an objective indication for the reproducibility, the coefficient of variation v , the quotient of standard deviation and the arithmetic mean, was calculated and plotted against the shear rate in Figure 13. As expected, the values of v are smaller for the second individual measurements. It is also noticeable that the coefficients in both cases decrease with the increase in the shear rate, which can be explained by the concrete behaviour at low shear rates already described.

Figure 10 
                     Representative values of the adjusted concrete mix for the first individual measurements of the 
                           
                              
                              
                                 15
                              
                              15
                           
                         measurement series.
Figure 10

Representative values of the adjusted concrete mix for the first individual measurements of the 15 measurement series.

Figure 11 
                     Representative values of the adjusted concrete mix for the second individual measurements of the 
                           
                              
                              
                                 15
                              
                              15
                           
                         measurement series.
Figure 11

Representative values of the adjusted concrete mix for the second individual measurements of the 15 measurement series.

Figure 12 
                     Selected representative values of the adjusted concrete mix of the first individual measurements (a) and the second individual measurements (b).
Figure 12

Selected representative values of the adjusted concrete mix of the first individual measurements (a) and the second individual measurements (b).

Figure 13 
                     Values of the coefficient of variation 
                           
                              
                              
                                 v
                              
                              v
                           
                         vs the shear rate for the first individual measurements (
                           
                              
                              
                                 
                                    
                                       v
                                    
                                    
                                       1
                                    
                                 
                              
                              {v}_{1}
                           
                        ) and the second individual measurements (
                           
                              
                              
                                 
                                    
                                       v
                                    
                                    
                                       2
                                    
                                 
                              
                              {v}_{2}
                           
                        ) with mean values of 
                           
                              
                              
                                 
                                    
                                       v
                                    
                                    
                                       1
                                    
                                 
                                 =
                                 0.033
                                 ±
                                 0.011
                              
                              {v}_{1}=0.033\pm 0.011
                           
                         and 
                           
                              
                              
                                 
                                    
                                       v
                                    
                                    
                                       2
                                    
                                 
                                 =
                                 0.026
                                 ±
                                 0.007
                              
                              {v}_{2}=0.026\pm 0.007
                           
                        .
Figure 13

Values of the coefficient of variation v vs the shear rate for the first individual measurements ( v 1 ) and the second individual measurements ( v 2 ) with mean values of v 1 = 0.033 ± 0.011 and v 2 = 0.026 ± 0.007 .

In principle, the evaluation of the coefficient of variation depends on the considered area, but in general, coefficients below 0.1 can be described as good [54], which shows the high reproducibility of the ACCR with v max < 0.05 .

5.2.2 Rheological characterisation of fresh concrete

Since the previous obtained values showed a very good reproducibility, they were averaged for the individual measurements 1 and 2 to obtain only one curve for each individual measurement, Figure 14(a). The mean temperatures are T m , 1 = 21.99 ± 0.27 ° C and T m , 2 = 21.95 ± 0.21 ° C . The standard deviation is smaller than the line thickness and the thickness of the measuring points so that it is not presented in the diagrams. Additionally, some values are excluded (see marked values) and the curves are fitted with the common concrete models in Figure 14(a), which will be discussed later.

Figure 14 
                     (a) Representative mean values of the first and second individual measurements of the selected measurement series with excluded shear rates and fitted models. (b) Representative measured values of an additional measurement with increased shear rate.
Figure 14

(a) Representative mean values of the first and second individual measurements of the selected measurement series with excluded shear rates and fitted models. (b) Representative measured values of an additional measurement with increased shear rate.

While the course of the first individual measurements is almost linear for shear rates γ ̇ rep > 10 s 1 and for the previously critically considered high shear rates ( γ ̇ rep > 35 s 1 ), the course of the second individual measurements shows stronger curvatures both for γ ̇ rep > 10 s 1 and for γ ̇ rep > 35 s 1 . Thus, the first individual measurements indicate a Newtonian flow behaviour from γ ̇ rep > 10 s 1 , while the second individual measurements indicate a slightly shear thinning behaviour for 10 s 1 < γ ̇ rep < 30 s 1 , which appears to reverse the dilatant behaviour from γ ̇ rep > 30 s 1 . For shear rates, γ ̇ rep < 10 s 1 , the curves are similar, at least with respect to curvature, and suggest shear thinning behaviour. The shear stresses are in the range of τ rep 150 650 Pa throughout the measurement range for the first individual measurements, while for the second individual measurements, it is τ rep 200 700 Pa , which again can be explained by time hardening.

Since the course of the first individual measurements is now also almost linear for γ ̇ rep > 35 s 1 , the previously interpreted dilatant behaviour (c.f. Section 5.1) can be attributed to insufficient pre-shear (time and/or intensity). However, this is different for the second individual measurements with a more aged concrete mix. To get a better insight in the high shear region, further investigations were performed with an increased maximum shear rate from γ ̇ i > 45 to 47 s 1 , Figure 14(b).

Comparing Figure 14(a) and (b), it can be clearly seen that the onset of the apparent shear thickening shifts to higher shear rates, which again is an indication for insufficient pre-shearing, especially for the second individual measurements. The hysteresis measurements, see Section 5.1.2, also showed no increase in this shear rate range for the PSR, while showed an increase for the NSR. Therefore, the described effect can clearly be attributed to insufficient pre-shearing and not to dilatant flow behaviour. This also becomes obvious when looking at the course of the torque, see supplementary material S3, since the torque for the first measurement value has not yet reached a pure steady state value in the second individual measurement, but has already done so in the first individual measurement. Therefore, the measured value at the highest shear rate of the second individual measurements is not included in the further evaluation. Again, the effect of gap reduction occurred at low shear rates, as discussed in detail in Section 5.1.1. The start of the gap reduction can only be determined to a limited extent via the measured values, which makes it a critical problem for the measurement. It is only possible to identify a kind of “step course” in both the individual measurements from a shear rate of approx. γ ̇ rep < 5 s 1 , which >somehow coincides with the optical observation of the gap reduction starting from a shear rate of γ ̇ rep 5 s 1 for the first individual measurements and from γ ̇ rep 8 s 1 for the second individual measurements. Therefore, this low shear region was also excluded in the following characterisation, c.f. Figure 14(a).

Looking at the first individual measurements, the values show a linear and thus Newtonian behaviour. Although the first point ( γ ̇ rep = 38.13 s 1 ) suggests a slightly increasing tendency, this cannot be assessed with certainty at this point and can be interpreted as a fluctuation in the measured values. The Newtonian behaviour is also obvious in the very good model fits (Bingham, Herschel–Bulkley, modified Bingham model), see Table 3, although the Herschel–Bulkley model ( n = 1.08 ) and the modified Bingham model ( c = 0.029 ) indicate a very slight dilatancy. However, these parameters are close to one and zero, respectively, and thus have little to no influence, see equations (8) and (9). Furthermore, when comparing the rheological parameters, it is noticeable that both the yield stresses and the characteristic values describing the viscosities are in a similar range ( τ y 200 Pa , k and η pl 10 Pa s ( n ) ). Thus, to this extent, all models resemble the concrete behaviour in a very good manner in a medium shear rate range.

Table 3

Model parameters of the Bingham, Herschel–Bulkley, and modified Bingham models for the fits of the mean values of the first two individual measurements

First individual measurement Second individual measurement
Bingham R 2 = 0.9983 R 2 = 0.9956
τ y / Pa 198.24 300.85
η pl / Pa s 11.78 10.36
Herschel–Bulkley R 2 = 0.9987 R 2 = 0.9983
τ y / Pa 213.44 216.08
k / Pa s n 8.61 35.12
n / 1.08 0.71
Mod. Bingham R 2 = 0.9988 R 2 = 0.9977
τ y / Pa 208.89 272.95
η pl / Pa s 10.54 13.20
c / Pa s 2 0.029 0.062

The situation is different for the second individual measurements. As already described in the consideration of the entire shear rate range, a slight tendency to shear thinning up to γ ̇ rep 25 s 1 is evident here, which then changes to Newtonian behaviour. Thus, the yield stress for the Bingham model is clearly above the value of the first individual measurement and also marks the maximum value when compared with the values of the other models. The plastic viscosity, on the other hand, deviates only slightly from the value of the first individual measurement, so that the Bingham model, roughly speaking, would only mark an increase in the yield stress from the first to the second individual measurement. The situation is similar for the modified Bingham model. There is an increase in yield stress, but less pronounced and also a slight rise in plastic viscosity can be seen. The negative coefficient c = 0.062 Pa s 2 indicates an overall slightly shear thinning behaviour, but, as in the first individual measurements, it is very close to zero. The parameters of the Herschel–Bulkley model have developed in a purely different way. There is a strong increase in consistency, which is about four times that of the first individual measurements. The flow exponent has also changed from an approximately Newtonian behaviour ( n = 1.08 ) to a clearly pronounced shear thinning behaviour ( n = 0.71 ). However, the yield stress remains almost identical. This is a very interesting aspect, since it basically argues against the application of the Herschel–Bulkley model. In all previously performed series of measurements, there has been an increased shear stress with the increase in the individual measurement, which has been justified by the increasing structure formation over time.

Usually, the yield stress is related to exactly such phenomena, which means that an increased structure formation must also result in an increase in the yield stress. Due to the NSR, it could be argued that an increasing destruction of the structures takes place over the measurement period and that both the individual measurements have the same initial state and thus the same yield stress at the end of the measurement ( γ ̇ rep = 0 ). However, this would require that even low shear rates can prevent structure formation, which can be excluded due to the general trend of higher values of the second individual measurement. Therefore, it is to conclude that the Herschel–Bulkley model tends to underestimate the yield stress and the other models are better suited to describe the rheological behaviour of fresh concrete in the investigated shear rate range.

In future investigations, the characterisation of fresh concrete should be extended to shear rates γ ̇ rep < 1 , which would allow the applicability of the flow models and the flow behaviour of the concrete to be analysed even more. However, due to the onset of gap reduction, this was not possible in this work and would only be feasible by constructing a much larger rheometer with a radius ratio of β 1 .

5.2.3 Investigation of sedimentation tendency

In order to assess the sedimentation and segregation tendency, additional measurements were carried out during the adjustment of the concrete mix and during the second test series. These consisted of a so-called cylinder sedimentation test [55] and representative sample taking after individual measurements.

The execution of the cylinder sedimentation test is shown schematically in Figure 15. For this purpose, three stacked cylinders ( d i = 150 mm , h cylinder = 150 mm , h total = 450 mm ), which have a total volume of approx. V total = 8 l , are filled with concrete directly after the mixing process. Following a rest period of t = 30 min , the cylinders are individually analysed and the grains are separated via washing and sieving. The mesh sizes used for sieving were 2.0 and 2.3 mm . After the grains have been dried, the mass is determined and the maximum degree of sedimentation m is defined according to equation (13) [55]. In the case of standard concretes with largest grains 8 mm , sufficient sedimentation resistance is said to exist if m does not exceed 20 % [55].

(13) m = max m i m ̅ 1 i = 1 3 = max 3 m i i = 1 3 m i 1 i = 1 3 .

Figure 15 
                     Schematic representation of the cylinder sedimentation test according to ref. [55].
Figure 15

Schematic representation of the cylinder sedimentation test according to ref. [55].

For the concrete investigated in this work with largest aggregates 4 mm , the results for both the mesh sizes are given in Table 4. It can be seen that the degrees of sedimentation of m = 3.25 % (mesh size 2.3 mm ) and m = 4.16 % (mesh size 2.0 mm ) are achieved, marking the system as sedimentation stable. The weight distribution is random and no trend for higher grain content at the bottom could be seen. However, the cylinder sedimentation test is only suitable for characterizing concretes that have not been stressed yet and therefore has a limited significance for sedimentation or segregation in a rheometer.

Table 4

Results of the cylinder sedimentation test

Cylinder Mass of aggregates ( g ) Deviation to mean value ( % )
Mesh size 2.3 mm Upper 785.40 3.25
Middle 749.50 1.47
Lower 747.10 1.78
m ̅ = 760.67 m = 3.25
Mesh size 2.0 mm Upper 896.89 4.16
Middle 829.90 3.61
Lower 856.20 0.55
m ̅ = 860.97 m = 4.16

Therefore, a lance was constructed to take samples directly from the rheometer gap at different heights. The lance has four sealable cups (A1 (top) to A4 (bottom)) with a volume of 11.65 c m 3 each, which are evenly distributed over the measuring gap length of L = 300 mm , see Figure 16. For sampling, the closed lance was inserted into the gap after a complete series of measurements. As soon as it was in position, the cups were opened and the lance was moved slightly back and forth so that the concrete could flow into the cups. The cups were then closed and the lance pulled out of the gap.

Figure 16 
                     CAD construction of sample lance; left = bottom and right = top.
Figure 16

CAD construction of sample lance; left = bottom and right = top.

After the individual samples were washed over a sieve with a mesh size of 700 µ m , the particle size distributions have been determined optically via microscope (Carl Zeiss International AG (Axiocam 506 colour)) using the software “ZEN core V3.1,” Figure 17. The particle diameter d E corresponds to the diameter of a sphere equal to the projection area and based on this diameter, the number distributions q 0 ( d E ) were calculated in size classes of 100 µ m . The number distributions clearly show that there are no segregation and sedimentation effects, since no trend can be detected as to whether a certain particle size is predominantly located at a certain height of the measuring gap. Thus, the concrete is not only sedimentation-stable before it is sheared, but also shows no segregation during shear.

Figure 17 
                     Particle size distributions of the sample cups (A1–A4) for a representative measurement.
Figure 17

Particle size distributions of the sample cups (A1–A4) for a representative measurement.

6 Conclusion

The ACCR presented in this study was developed for the purpose of determining the rheological properties of suspensions, primarily fresh concrete, with a maximum particle size of d max = 5.5 mm . Comparative measurements with a commercial rheometer (Anton Paar MCR 501) have shown that the ACCR is able to generate precise measurement results within a common accuracy. It was possible to create a SCC concrete mixture with nearly no tendency of segregation and sedimentation, that was used for investigation.

After the difficulties of concrete measurements encountered in the preliminary tests had been solved, the fresh concrete could be rheologically characterised with extremely high reproducibility. The applicability of the common flow models for concrete (Bingham, Herschel–Bulkley, and modified Bingham model) could be satisfactorily analysed. It was shown that especially the Bingham and modified Bingham model are suitable for the description of the concrete used in this work, while the Herschel–Bulkley model seems to underestimate the yield stress and represents shear-thinning behaviour. Furthermore, only the values of the first individual measurements are valid for further evaluation of rheological properties, since the influence of structural build-up in the second individual measurements is too high to be neglected and can lead to a false interpretation of the rheological behaviour.

However, a limited shear rate range has resulted due to the yield stress and the structural build-up of the concrete, which in particular have led to a reduced measurement gap. This mainly affected the low shear rate range. Within the scope of this work, this phenomenon could not be avoided and can also only be reliably fixed in the future by redesigning the rheometer with a β closer to 1 . The investigation of the low shear rate range would be of great interest in the future, as it would not only allow the applicability of the current rheological models to be better investigated, but also the yield stress and therefore the real flow behaviour to be determined much more accurately.



  1. Funding information: We would like to thank the German Research Foundation (DFG) for supporting the research within the priority program SPP 2005 “Opus Fluidum Futurum.”

  2. Author contributions: Conceptualisation: S.J. and S.J.; data curation: S.J.; formal analysis: S.J. and S.J.; funding acquisition: H.-J.S. and S.J.; investigation: S.J.; methodology: S.J. and S.J.; project administration: S.J. and H.-J.S.; supervision: H.-J.S.; validation: S.J. and S.J.; visualisation: S.J.; writing – original draft: S.J.; writing – review and editing: S.J. and H.-J.S.

  3. Conflict of interest: The authors declare no conflict of interest.

  4. Ethical approval: The conducted research is not related to either human or animal use.

  5. Data availability statement: The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

References

[1] Roussel N, Geiker MR, Dufour F, Thrane LN, Szabo P. Computational modeling of concrete flow: General overview. Cem Concr Res. 2007;37(9):1298–307. 10.1016/j.cemconres.2007.06.007.Search in Google Scholar

[2] Ferraris CF, Billberg P, Ferron R, Feys D, Hu J, Kawashima S, et al. Role of rheology in achieving successful concrete performance. Concr Int. 2017;39(6):48.Search in Google Scholar

[3] Roussel N. Rheological requirements for printable concretes. Cem Concr Res. 2018;112:76–85. 10.1016/j.cemconres.2018.04.005.Search in Google Scholar

[4] Vasilic K, Gram A, Wallevik JE. Numerical simulation of fresh concrete flow: insight and challenges. RILEM Tech Lett. 2019;4:57–66. 10.21809/rilemtechlett.2019.92.Search in Google Scholar

[5] Uzomaka OJ. A concrete rheometer and its application to a rheological study of concrete mixes. Rheol Acta. 1974;13:12–21.10.1007/978-3-662-41458-3_78Search in Google Scholar

[6] Tattersall GH, Bloomer SJ. Further development of the two-point test for workability and extension of its range. Mag Concr Res. 1979;31(109):205.Search in Google Scholar

[7] Weber M, Dedegil MY. Rotationsviskosimeter. Utility patent 0166332; 1985.Search in Google Scholar

[8] Wallevik OH. The Rheology of Fresh Concrete and its Application on Concrete with and without Silica Fume: Dissertation, NTH Trondheim; 1990.Search in Google Scholar

[9] Wallevik OH, Gjorv OE. Development of a coaxial cylinder viscometer for fresh concrete. RILEM Proceedings; 1990. p. 213–24.Search in Google Scholar

[10] Coussot P. Rheologie des boues et laves torrentielles – Etudes de dispersions et suspensions concentrées. PhD thesis. Versailles: l’institut National Polytechnique de Grenoble; 1993.Search in Google Scholar

[11] Beaupre D. Rheology of High Performance Shotcrete. PhD thesis. Vancouver: University of British Colombia; 1994.Search in Google Scholar

[12] Hu C. Rhéologie des bétons fluides. PhD thesis. Versailles: Ecole Nationale des Ponts et Chaussées; 1995.Search in Google Scholar

[13] Bartos P, Marrs DL, Cleland DJ, editors. Production methods and workability of concrete. Proceedings of the International RILEM Conference, Paisley, Scotland, June 3–5, 1996. London: E & FN Spon; 1996.10.1201/9781482271782Search in Google Scholar

[14] Ferraris CF, Banfill PFG, Chapdelaine F, Larrard FD, Domone PLJ, Nachbaur L, et al. Comparison of concrete rheometers: International tests at LCPC (Nantes, France) in October, 2000. NISTIR 6819 2001.10.6028/NIST.IR.6819Search in Google Scholar

[15] Koehler EP, Fowler DW. Development of a portable rheometer for fresh portland cement concrete. Texas: International Center for Aggregates and Research; 2004.Search in Google Scholar

[16] Jau W-C, Yang C-T. Development of a modified concrete rheometer to measure the rheological behavior of conventional and self-consolidating concretes. Cem Concr Compos. 2010;32(6):450–60. 10.1016/j.cemconcomp.2010.01.001.Search in Google Scholar

[17] Estellé P, Lanos C. High torque vane rheometer for concrete: principle and validation from rheological measurements. Appl Rheol. 2012;22(1):3. 10.3933/ApplRheol-22-12881.Search in Google Scholar

[18] Roussel N. Understanding the rheology of concrete. Cambridge, Philadelphia, PA: Woodhead Publishing; 2012.10.1533/9780857095282Search in Google Scholar

[19] Soualhi H, Kadri EH, Ngo T-T, Bouvet A, Cussigh F, Kenai S. A new vane rheometer for fresh mortar: development and validation. Appl Rheol. 2014;24(2):1–7. 10.3933/APPLRHEOL-24-22594.Search in Google Scholar

[20] Coussot P, Piau J-M. A large‐scale field coaxial cylinder rheometer for the study of the rheology of natural coarse suspensions. J Rheol. 1995;39(1):105–24. 10.1122/1.550693.Search in Google Scholar

[21] Larrard FD, Szitkar J-C, Hu C, Joly M, Derkx F. Design of a rheometer for fluid concretes. RILEM Workshop Spec Concr – Workability Mixing. 1993;186:201–8.Search in Google Scholar

[22] Larrard FD, Sedran T, Hu C, Szitkar J-C, Joly M. Evolution of the workability of superplasticized concretes: assessment with BTRHEOM rheometer. RILEM International Conference on Production Methods and Workability of Concrete. RILEM Proceedings 32; 1996. p. 377–88.Search in Google Scholar

[23] Szecsy RE. Concrete Rheology. PhD thesis. Illinois: University of Illinois at Urbana Champaign; 1997.Search in Google Scholar

[24] Struble L, Puri U, Ji X. Concrete rheometer. Adv Cem Res. 2001;13(2):53–63.10.1680/adcr.2001.13.2.53Search in Google Scholar

[25] Thrane LN, Pade C, Nielsen CV, Jeknavorian AA, Schemmel JJ, Dean SW. Determination of rheology of self-consolidating concrete using the 4C-Rheometer and how to make use of the results. J ASTM Int. 2010;7(1):102003. 10.1520/JAI102003.Search in Google Scholar

[26] Reiner M. Deformation and flow. An elementary introduction to theoretical rheology. Great Britain: H.K. Lewis & Co, Limited; 1949.Search in Google Scholar

[27] Ferraris CF, Beaupre D, Chapdelaine F, Domone P, Koehler E, Shen L, et al. Comparison of concrete rheometers: International tests at MB (Cleveland OH, USA) in May, 2003. NISTIR 7154 2004.10.6028/NIST.IR.7154Search in Google Scholar

[28] Pahl M, Gleißle W, Laun H-M. Praktische Rheologie der Kunststoffe und Elastomere. Düsseldorf: VDI-Verlag; 1991.Search in Google Scholar

[29] Mezger T. Das Rheologie-Handbuch: Für Anwender von Rotations- und Oszillations-Rheometern. 4. Aufl. Hannover: Vincentz Network; 2012.Search in Google Scholar

[30] Schümmer P. Zur Darstellung der Durchflußcharakteristik, viskoelastischer Flüssigkeiten in Rohrleitungen. Chem Ing Technik. 1970;42(19):1239. 10.1002/cite.330421916.Search in Google Scholar

[31] Giesekus H, Langer G. Die Bestimmung der wahren Fließkurven nicht-newtonscher Flüssigkeiten und plastischer Stoffe mit der Methode der repräsentativen Viskosität. Rheol Acta. 1977;16(1):1–22. 10.1007/BF01516925.Search in Google Scholar

[32] Krieger IM, Elrod H. Direct determination of the flow curves of non‐Newtonian fluids. II. shearing rate in the concentric cylinder viscometer. J Appl Phys. 1953;24(2):134–6. 10.1063/1.1721226.Search in Google Scholar

[33] Banfill PFG. A viscometric study of cement pastes containing superplasticizers with a note on experimental techniques. Mag Concr Res. 1981;33(114):37–47. 10.1680/macr.1981.33.114.37.Search in Google Scholar

[34] Flatt RJ. Towards a prediction of superplasticized concrete rheology. Mater Struct. 2004;37(5):289–300. 10.1007/BF02481674.Search in Google Scholar

[35] Yahia A, Mantellato S, Flatt RJ. Concrete rheology. In: Science and Technology of Concrete Admixtures. Amsterdam, Netherlands: Elsevier; 2016. p. 97–127.10.1016/B978-0-08-100693-1.00007-2Search in Google Scholar

[36] Heirman G, Vandewalle L, van Gemert D, Wallevik Ó. Integration approach of the Couette inverse problem of powder type self-compacting concrete in a wide-gap concentric cylinder rheometer. J Non-Newtonian Fluid Mech. 2008; 150(2–3):93–103. 10.1016/j.jnnfm.2007.10.003.Search in Google Scholar

[37] Feys D, Verhoeven R, De Schutter G. Fresh self compacting concrete, a shear thickening material. Cem Concr Res. 2008;38(7):920–9. 10.1016/j.cemconres.2008.02.008.Search in Google Scholar

[38] Feys D, Verhoeven R, De Schutter G. Why is fresh self-compacting concrete shear thickening? Cem Concr Res. 2009;39(6):510–23. 10.1016/j.cemconres.2009.03.004.Search in Google Scholar

[39] Yahia A. Shear-thickening behavior of high-performance cement grouts – Influencing mix-design parameters. Cem Concr Res. 2011;41(3):230–5. 10.1016/j.cemconres.2010.11.004.Search in Google Scholar

[40] Feys D, Wallevik JE, Yahia A, Khayat KH, Wallevik OH. Extension of the Reiner–Riwlin equation to determine modified Bingham parameters measured in coaxial cylinders rheometers. Mater Struct. 2013;46(1–2):289–311. 10.1617/s11527-012-9902-6.Search in Google Scholar

[41] Larrard FD, Ferraris CF, Sedran T. Fresh concrete: A Herschel-Bulkley material. Mater Struct. 1998;31(7):494–8. 10.1007/BF02480474.Search in Google Scholar

[42] Yahia A, Khayat K. Analytical models for estimating yield stress of high-performance pseudoplastic grout. Cem Concr Res. 2001;31(5):731–8. 10.1016/S0008-8846(01)00476-8.Search in Google Scholar

[43] Feys D, Heirman G, Schutter G, de Verhoeven R, Vandewalle L, van Gemert D. Comparison of two concrete rheometers for shear thickening behaviour of SCC. Proc 5th Int RILEM Symp SCC. 2007;2007:365–70.Search in Google Scholar

[44] Feys D, Verhoeven R, de Schutter G. Evaluation of time independent rheological models applicable to fresh self-compacting concrete. Appl Rheol. 2007;17(5):56244-1–10. 10.1515/arh-2007-0018.Search in Google Scholar

[45] Ferraris CF, Martys N. De la pâte de ciment au béton: modélisation et mesures expérimentales des propriétés rhéologiques. Proc. Rhéologie Génie Civil et Environment; 2001. p. 226–30.Search in Google Scholar

[46] Deutsches Institut für Normung e.V. DIN EN ISO 3219:1994-10. Kunststoffe - Polymere/Harze in flüssigem, emulgiertem oder dispergiertem Zustand - Bestimmung der Viskosität mit einem Rotationsviskosimeter bei definiertem Geschwindigkeitsgefälle (ISO 3219:1993); Deutsche Fassung EN ISO 3219:1994. Berlin: Beuth Verlag GmbH.Search in Google Scholar

[47] Mooney M. Explicit formulas for slip and fluidity. J Rheol. 1931;2(2):210–22. 10.1122/1.2116364.Search in Google Scholar

[48] Yoshimura A, Prud'homme RK. Wall slip corrections for couette and parallel disk viscometers. J Rheol. 1988;32(1):53–67. 10.1122/1.549963.Search in Google Scholar

[49] Kiljański T. A method for correction of the wall-slip effect in a Couette rheometer. Rheol Acta. 1989;28(1):61–4. 10.1007/bf01354770.Search in Google Scholar

[50] Jesinghausen S, Weiffen R, Schmid H-J. Direct measurement of wall slip and slip layer thickness of non-Brownian hard-sphere suspensions in rectangular channel flows. Exp Fluids. 2016;57(9):159. 10.1007/s00348-016-2241-6.Search in Google Scholar

[51] Pawelczyk S, Kniepkamp M, Jesinghausen S, Schmid H-J. Absolute rheological measurements of model suspensions: Influence and correction of wall slip prevention measures. Mater (Basel). 2020;13(2):31963746. 10.3390/ma13020467 [PMID: 31963746].Search in Google Scholar PubMed PubMed Central

[52] Weston JS, Harwell JH, Grady BP. Rheological characterization of yield stress gels formed via electrostatic heteroaggregation of metal oxide nanoparticles. Soft Matter. 2017;13(38):6743–55. 10.1039/c7sm01035d [PMID: 28819663].Search in Google Scholar PubMed

[53] Narayanan N, Ramamurthy K. Structure and properties of aerated concrete: a review. Cem Concr Compos. 2000;22(5):321–9. 10.1016/S0958-9465(00)00016-0.Search in Google Scholar

[54] Gomez KA, Gomez AA. Statistical procedures for agricultural research. 2 edn New York: Wiley; 1984.Search in Google Scholar

[55] Lowke D, Wiegrink K-H, Schiessl P A simple and significant segregation test for SCC. Proceedings of 3rd International RILEM Symposium on Self Compacting Concrete; 2003.Search in Google Scholar

Received: 2022-11-28
Revised: 2023-01-12
Accepted: 2023-01-19
Published Online: 2023-03-02

© 2023 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 11.12.2023 from https://www.degruyter.com/document/doi/10.1515/arh-2022-0140/html
Scroll to top button