Next Article in Journal
Laser Powder-Bed Fusion as an Alloy Development Tool: Parameter Selection for In-Situ Alloying Using Elemental Powders
Next Article in Special Issue
Microstructural Analysis and Mechanical Properties of TiMo20Zr7Ta15Six Alloys as Biomaterials
Previous Article in Journal
The Effect of the Type and Amount of Synthetic Fibers on the Effectiveness of Dispersed Reinforcement in Soil-Cements
Previous Article in Special Issue
Contactless Measurements of Carrier Concentrations in InGaAs Layers for Utilizing in InP-Based Quantum Cascade Lasers by Employing Optical Spectroscopy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detection and Quantification of Cracking in Concrete Aggregate through Virtual Data Fusion of X-Ray Computed Tomography Images

1
Bundesanstalt für Materialforschung und–prüfung, BAM (Federal Institute for Materials Research and Testing), 12205 Berlin, Germany
2
Institute of Physics and Astronomy, University of Potsdam, Karl-Liebknecht-Str.24-25, 14476 Potsdam, Germany
*
Author to whom correspondence should be addressed.
Materials 2020, 13(18), 3921; https://doi.org/10.3390/ma13183921
Submission received: 24 July 2020 / Revised: 31 August 2020 / Accepted: 3 September 2020 / Published: 4 September 2020
(This article belongs to the Special Issue Micro Non-destructive Testing and Evaluation)

Abstract

:
In this work, which is part of a larger research program, a framework called “virtual data fusion” was developed to provide an automated and consistent crack detection method that allows for the cross-comparison of results from large quantities of X-ray computed tomography (CT) data. A partial implementation of this method in a custom program was developed for use in research focused on crack quantification in alkali-silica reaction (ASR)-sensitive concrete aggregates. During the CT image processing, a series of image analyses tailored for detecting specific, individual crack-like characteristics were completed. The results of these analyses were then “fused” in order to identify crack-like objects within the images with much higher accuracy than that yielded by any individual image analysis procedure. The results of this strategy demonstrated the success of the program in effectively identifying crack-like structures and quantifying characteristics, such as surface area and volume. The results demonstrated that the source of aggregate has a very significant impact on the amount of internal cracking, even when the mineralogical characteristics remain very similar. River gravels, for instance, were found to contain significantly higher levels of internal cracking than quarried stone aggregates of the same mineralogical type.

1. Introduction

1.1. Alkali-Silica Reaction (ASR)

Despite decades of research, the problem of harmful alkali-silica reaction (ASR) in the field of concrete construction has not yet been satisfactorily solved. For the first time in 1940, Stanton [1] reported damaging strains within concrete due to chemical reactions of cement and aggregate. In the 1950s, Powers and Steinour [2,3] developed initial models of ASR’s damage mechanism. In the 1970s, Locher and Sprung [4] identified opal and porous flint as alkali-sensitive aggregates and developed theories on their reaction mechanisms. In the 1980s, various researchers conducted in-depth studies on the influence of alkali metal salts on the swelling pressures of the ASR gel [5,6]. The current state of knowledge in the field of ASR has also been extensively described in number of recent publications [7,8].
During the ASR process, the reactive SiO2 within aggregates reacts with alkalis (supplied from the cement paste or from an external source) in the presence of water to form expansive alkali silicate hydrates (ASR gels). Because the tensile strength of road surface concrete is often significantly lower than the swelling pressures caused by ASR gels, cracking can be induced [9]. The progression and the extent of the resulting cracking processes are determined to a large extent by the type of aggregate. For example, fast-reacting aggregates (among others, flint, opaline sandstone and mudstone) are characterized by gel and crack formation emanating from the transition zone between the grain and the mortar matrix [7]. On the other hand, in the case of the slow-reacting aggregates (for instance greywacke or quartz porphyry), which were of primary of interest in this project, the gel formation takes place above all inside the aggregate itself, which results in the formation of internal aggregate cracks [7].
The severity of the ASR degradation process is thought to be partially dependent on the amount of porosity within a given aggregate that is accessible to liquid penetrating from the sample surface. Against this background and to evaluate the alkali sensitivity of the aggregate, the influence of the specific surface area on the solubility behaviour of four different aggregates in 0.1 M potassium hydroxide solution without and with defined addition of NaCl at a temperature of 80 °C was thoroughly investigated in a joint project [10]. To quantify this relationship, a non-destructive method is needed for measuring both the external surface area of aggregates and the internal surface area of aggregate cracks and pores, including a differentiation of internal voids connected to the sample surface from those isolated from the surface. The primary focus of this publication is devoted to the crack detection method. A detailed analysis of the implications for ASR damage, including a comparison of the CT results with those from other porosity measurement methods, such as mercury porosimetry and the Brunauer–Emmett–Teller (BET) method, can be found in Oesch et al. 2020 [11] and in Weise et al. 2019 [10].

1.2. X-Ray Computed Tomography (CT)

The development of X-ray computed tomography (CT) began in the 1960s and clinical X-ray CT investigations have been widely conducted since the 1970s [12,13]. Since that time, many different reconstruction algorithms have been developed for clinical use, including algorithms based on the algebraic reconstruction technique (ART), filtered back projection (FBP), and iterative reconstruction (IR) [14]. Medical X-ray CT scanning systems are, however, unsuitable for many applications in materials science given their lower X-ray energy characteristics and coarser resolution compared to specialized laboratory-based X-ray CT systems used for materials research [15]. These differences occur both because the size of the intended scanning objects tends to significantly differ between clinical and materials science applications and because the X-ray absorption characteristics of live tissues are much lower than those of materials such as concrete and steel.
X-ray CT has been used in non-destructive concrete research applications for more than 30 years [16,17]. In this scanning method, a sample is placed on a rotating table between an X-ray source and an X-ray detector [18]. By adjusting the distances between the X-ray source, the sample and the X-ray detector, it is possible to vary the voxel (i.e., 3D pixel) resolution in the resulting images. The penetration of the sample by the X-ray beam causes an X-ray attenuation image of the sample to be projected upon the detector. By recording these projected images during the 360° rotation of the sample, the projections can be inverted using volume reconstruction algorithms, which produce a 3D representation of X-ray attenuation within the sample [19]. X-ray attenuation is approximately proportional to local material density and can be used to identify single objects within a material (or structure) and to individually separate and analyze those objects.
Previous research has shown that X-ray CT scans can be taken during incremental testing. This includes mechanical testing (such as unconfined compression, split cylinder, triaxial, and reinforcing bar pull-out testing [20,21,22,23,24]), chemical testing (such as the measurement of progressive corrosion in reinforced concrete during repeated exposure to chloride [25] and the transport of water [26,27,28]) and thermal testing (such as water migration in heated concrete [29]).

1.3. Crack Detection and Quantification

Crack detection and quantification is important for understanding and modelling a series of material behaviours. Precise measurement of crack surface area is needed, for instance, in order to calculate the fracture energy expended during damage processes using basic fracture mechanics relationships [30]. The crack orientations have also been observed to exhibit a behaviour that is highly dependent on the anisotropy of the material structure, such as fibre orientation within fibre-reinforced concretes [31].
Crack detection and quantification within X-ray CT images have been the subject of extensive past research. Most of these crack detection methods have leveraged one or more unique characteristics of cracks, which differentiate them from the surrounding material. Possibly the most popular research approaches have focused on the use of template-matching methods in order to separate cracks from the surrounding materials [32,33]. This method of crack detection relies on the similarity of cracking structures to certain template shapes, such as small planes or discs. Although impressive results have been demonstrated using the template-matching method, the template parameters are not universal and must generally be tailored for each material and each imaging scenario.
Research was recently carried out by Paetsch (2019) [34] with the goal of partially overcoming these challenges related to using template matching approaches. This research indicated that the results of a series of analyses carried out using different template shapes can be combined in order to obtain a greater accuracy of the detected cracks. However, Paetsch (2019) [34] has underlined that further problems remain to be solved that are common to most template-matching methods, such as difficulties detecting cracks in areas where significant crack branching or widening occurs.
Another method that takes advantage of the narrow shape characteristics of cracks is a Hessian-based approach [35,36], which identifies regions that exhibit sharp changes in image intensity. Percolation methods have also shown significant promise in detecting cracks across a range of materials [37]. These methods leverage the fact that most cracks are continuous, narrow objects with relatively consistent (low) density. Percolation methods have proven insufficient, however, to accurately detect complex cracks of varying size in most materials. Impressive results have also been obtained through a combination of the Hessian and percolation-based methods into the Hessian-driven percolation approach, although the processing time required for such an analysis remains prohibitive for most high-resolution CT images [32].
Many methods have also been employed that leverage the unique characteristics of a specific material or damage scenario. One excellent example of this approach is the use of digital volume correlation (DVC) to detect cracks in samples subjected to in-situ loading [38,39,40]. DVC is used to measure strains within samples by calculating voxel movements between subsequent CT images (such as images calculated before and after a loading increment). Cracks can, thus, be identified as areas of either high strain or poor correlation within DVC images or through more complex analysis methods, such as phase-congruency analysis. Although these methods appear to exhibit a relatively high accuracy (even sub-pixel, see [41]), they are only useful for detecting cracks caused by progressive in-situ testing with simultaneous CT. They do not provide any benefit for detecting cracks that are already present in specimens prior to testing. Phase congruency can, however, also be used to detect edge features (including cracks) within the greyscale images without direct reference to DVC [42].
Another example of a material-specific crack detection approach is the leveraging of typical wood structure within logs to identify cracks. These cracks typically run perpendicular to the growth rings of the logs and can, thus, be easily identified by their orientation characteristics [43].
Despite the many promising crack detection methods outlined here, a series of obstacles remain that have prevented the implementation of consistent, accurate, and quantitative crack analyses as part of CT scanning. First, most of these methods require some amount of tailoring for specific material properties (as in the case of template matching and percolation) or specific crack conditions (such as measuring only those arising from in-situ testing using DVC methods). Second, none of these methods has been successfully validated for the quantitative determination of crack properties (such as surface area) based on other standard measurement techniques.
In this study, we will show, on the example of different aggregate types, how the obstacles mentioned above can be circumvented by a novel data fusion strategy.

2. Materials and Methods

2.1. Sample Selection

For this testing series, a group of aggregates from four different categories were selected and analysed. These categories were selected in order to include stones with a variety of different mineralogical compositions, deterioration conditions, and alkali-sensitivity characteristics (Table 1 and Figure 1) [44]. All stones included in this analysis were in the 8 mm to 16 mm size range: During the sieving process, the individual stones passed successfully through a sieve with a 16 mm mesh but were unable to pass through a sieve with an 8 mm mesh [45].
For the categories GK1 and GK4, the selected stones were sieved from a natural river gravel. Such river gravels are typically characterized by significant mineralogical deterioration due to naturally occurring weathering processes. In order to characterize these river gravels, which have a heterogeneous mineralogical composition, the primary types of rock (7 types in total) occurring in GK1 and GK4 were determined and 10 individual grains of each type were selected for CT analysis.
For the categories GK2 and GK3, the selected stones were sieved from crushed stone chips that were quarried from solid rock deposits. Such quarried stones are typically characterized by greater mineralogical integrity than river gravels because they have not been exposed to significant weathering. Due to the homogeneity of the quarried stone, the selection of individual grains was limited to 10 each for GK2 and GK3. Thus, combining all the individual grains from each of the four stone categories included in this research study, a total of 90 different individual grains were investigated.

2.2. CT Scanning

During this research program, an acceleration voltage of 130 kV and current of 180 µA were used for the X-ray source. The X-ray beam was also filtered using a 0.5 mm thick Copper plate immediately upon leaving the source in order to remove (unwanted) photons of small energies from the X-ray beam, thereby increasing the contrast of the resulting images. The flat panel detector used for this scanning contained a 2048 × 2048 pixel field.
Individual aggregates were sorted based on mineralogical characteristics and placed within corresponding plastic tubes with small pieces of foam separating the aggregates from one another. As a result of heating and deterioration of the target material within the X-ray tube as well as changes in detector sensitivity over time, significant variations in the measured X-ray beam intensity and distribution can occur. To compensate for these variations, “dark-field” and “bright-field” images, which correspond to blank images (i.e., containing no sample) acquired with no illumination and full illumination, respectively, were acquired prior to the scanning of each plastic tube of samples. These images were then used to calibrate the X-ray images of the samples.
The CT machine was pre-programmed to collect a complete scan of each aggregate before repositioning the plastic tube using a manipulator and beginning the scan of the next aggregate. Thus, the scanning conditions for all stones within any given plastic tube were identical. All scanning conditions other than resolution were also held constant for all plastic tubes.
It was important to maintain scanning conditions that were as consistent as possible to ensure that the results of the crack analysis would be comparable. In spite of this, the resolution was maximised for each stone type, if they significantly varied in size and shape. This was done by adjusting the distance between the plastic tube and the X-ray source. The corresponding voxel sizes were then calculated directly from the measured distances between the X-ray source, the sample holder and the X-ray detector for each individual set of scans.
Although these variations in image resolution are known to directly affect measurements of crack properties, such as surface area (increased surface area is generally detected with decreasing voxel sizes), some estimation of the magnitude of this effect can already be accounted for based on a recent study [46]. The results of this study indicate that even a doubling of voxel size does not appear to generally change the measured crack surface area by more than a factor around two. We will see below that the small voxel size variations had little influence on the much larger variations in measured crack surface area among the different aggregate types, so that comparisons could be made, and conclusions drawn.

2.3. Image Analysis

2.3.1. Data Fusion Approach

In such an environment, where many analysis methods are readily available, but no individual method is sufficiently accurate to provide the needed measurements, the use of a data-fusion inspired approach becomes very promising and attractive. Data fusion enables researchers to combine data from multiple sources in order to produce more consistent and accurate results than those provided by any single source [47]. In non-destructive testing, the different sources used in data fusion typically result from differing non-destructive measurement approaches. Such data-fusion based approaches have even been successfully applied to the application of crack detection [48].
It is, however, often the case that researchers only possess meaningful data from a single measurement technique, such as CT (often indeed used as a reference), or do not have access to further non-destructive testing equipment. Even in this case, it should be possible to use the theory behind the data fusion approach to improve the overall quality of the quantitative determinations resulting from image analysis. Rather than relying on a wide variety of physical measurement techniques, data fusion in this approach would be carried out by using the output from an array of different image analysis techniques. Since such an approach relies on the fusion of different computationally generated data sources resulting from the same original CT image rather than on a variety of different physical measurement techniques, it can be more appropriately described as virtual data fusion. This is similar to the process that is thought to occur when an expert identifies cracks using the human eye. The expert has only one measurement technique (a visual image), but by considering its many different aspects (such as coloration, shape, and relationship to surrounding objects or planes of stress), the expert is able to quite easily and accurately detect a crack on the surface of a specimen.
Figure 2 shows a diagram of this virtual data fusion approach. In the diagram, the characteristics of a crack are displayed in dark grey boxes and the image processing steps for identifying objects with those characteristics are depicted in light grey boxes. The list provided here is only meant to serve as an example and is by no means exhaustive. It is clear, however, that results obtained through a fusion of the results from twelve such independent analyses will be much more accurate and resilient to varying material conditions than the results from any single analysis technique. An example of how such a step-by-step approach can identify individual cracks within a generic concrete sample is shown in Figure 3.

2.3.2. Implementation

Although, in theory, twelve or more different analysis methods could be used in the virtual data fusion implementation for analysing the ASR-sensitive aggregate dataset, in practice such a full implementation was impractical. We selected a partial implementation of the scheme including only three analysis approaches. The choice included the approaches that would produce cracking data with sufficient accuracy for this specific application while simultaneously testing the effectiveness of the virtual data fusion concept for obtaining accurate and consistent results. A successful demonstration of the virtual data fusion method with only three analysis components would then give justification for the further development of the algorithm to gradually include additional analysis modules, simultaneously growing in accuracy and resiliency. Moreover, the modular architecture of the strategy could more easily be adapted to different problems than a full but rigid analysis.
In fact, for the specific ASR-sensitive aggregate analysis described in this paper, it was not important to separate internal pores from internal cracks since the surface area of both pores and cracks was vulnerable to ASR degradation. Thus, modules related to this differentiation could be left out of the analysis. Furthermore, since scans were not available at varying levels of degradation, all methods relying on changes in sample state relative to time (such as DVC-based methods) were also left out of this analysis. All image analysis described in this paper was completed using custom algorithms developed and implemented in MATLAB [49].

Module 1: Identification of Objects with Low Density

When cracks or pores have widths exceeding two voxels, the voxels in their centres are completely filled with air. Thus, these central voxels are characterized by a particularly low greyscale value. In order to separate these voxels, a threshold has to be selected and subsequently used for image binarization. All voxels darker than the threshold would then be transformed to white and all voxels lighter than the threshold transformed to black. Given that this threshold must be consistently applied for a large range of datasets, an automated selection method was implemented.
Using the triangle selection algorithm [50,51], a virtual line is drawn from the origin of the image histogram to the top of the largest histogram peak (excluding the initial peak at zero, which represents voxels from the air around the specimen) (Figure 4). A calculation is then conducted to determine which point on the histogram is furthest from the virtual line along an intersecting, perpendicular line. The location of that point is identified as a potential threshold. It was found, however, that such a threshold is rather over-encompassing, leading to the introduction of considerable noise into the resulting binarized images. Thus, a slightly more conservative threshold value was also implemented for this module that is 25% lower than the original threshold. An example of the results from this analysis can be observed in Figure 5.

Module 2: Identification of Objects with High Gradient

When an object in a CT-image has either a much higher or a much lower greyscale value compared to the surrounding material, its edges can be identified as regions of high gradient. For this purpose, a gradient-magnitude image is calculated from the original CT-image through the use of the Sobel operator (Figure 6) [52]. Since the resulting gradient image consists of a wide range of grey values, it must also be binarized. For this purpose, a histogram of the gradient image grey values is computed and another automated triangular threshold selection is completed, this time from the right edge of the histogram (Figure 7). The resulting binary image contains the edges of both bright and dark objects in the original CT image (Figure 8). In order to remove the edges associated with the bright (high-density) objects in the CT-image as well as some of the noise, the binarized gradient image can be multiplied by a binarized CT-image (this time using the “over-encompassing void-solid threshold” identified in Figure 4) (Figure 9).

Module 3: Identification of Objects above a Given Size

Prior to noise identification and removal, the results from Modules 1 and 2 were combined into a single image (Figure 10). This ensured that the centres and edges of the cracks/pores would both be present and accounted for during object size analysis. The object size analysis was completed using a connected components algorithm. During this analysis, individual voxels are assessed to determine whether they are part of a larger object in the binary image by analysing whether they adjoin other voxels of the same color. For this analysis, white voxels with touching faces, edges or corners (also known as “26-connected”) were identified as connected. Objects with voxel volumes smaller than 125 voxels were subsequently identified and eliminated (Figure 11).
It was found that the use of smaller cubes than 125 voxels (with five voxel long sides) tended to leave considerable noise within the image while the use of larger cubes resulted in the loss of a considerable number of voxels along the path of the cracks. Given that the full CT images each contained over seven billion voxels, such an object with a volume of 125 voxels represented less than one ten-millionth of the total image volume. A flowchart of the entire crack-detection process is provided in Figure 12.

3. Results

After the cracks and internal voids were identified within each of the aggregates, their characteristics could be quantified and compared. Of primary interest for this analysis was the determination of volume and surface area characteristics. In particular, in order to compare the results of the CT-analysis with those of other non-destructive measurement techniques, it was important to separate surface-connected cracks and voids (referred here to as “open voids”) from internally isolated cracks and voids (referred to here as “closed voids”). This is because measurements of internal surface area using the Brunauer–Emmett–Teller (BET) method [53,54] only account for internal voids accessible from outside of the sample. Similarly, the measurement of void volume through mercury porosimetry is thought to depend primarily on the saturation of internal voids that are connected to the stone surface. Such a separation of surface-connected cracks could be completed using a connected components analysis in which only cracks and voids containing voxels that touched the stone surface were retained.
The crack/pore surface area measurements obtained using this approach are provided for all of the river-gravel type aggregates in Figure 13 and for both river-gravel and quarried-stone aggregates of the minerals rhyolite and greywacke in Figure 14. Note that the measurements for the greywacke (GK4) and rhyolite (GK1) river gravels appear in both figures for comparison purposes. Figure 15 also provides 3D images of cracking distributions for two selected individual grains of the same mineral (greywacke), where one grain has been extracted from a quarry and the other has been taken from river gravel. For images of all analysed aggregate types, see Appendix A. Tabulated values of the surface area measurements are also provided for each individual sample in Appendix B.
Error bars in Figure 13 and Figure 14 could not be calculated. This is because there is still no universally accepted method for estimating the combined error introduced by CT measurement systems and image processing algorithms. Although numerous approaches for estimating such error bars have been proposed [55,56], these tend to be rather computationally intensive and time consuming and remain an active area of research. For pure dimensional measurements, it is common to use either the voxel size or the focal spot size of the X-ray tube as an estimation of error. Given that the focal spot size of this scanning system was significantly smaller than the voxel sizes obtained during these investigations, the voxel sizes for each scan (also listed in Figure 13 and Figure 14) can be taken as an estimation of possible error in the dimensional measurements.
From Figure 14 it is clear that the amount of internal cracking (including both surface-connected cracking and non-surface-connected cracking) in the river gravel aggregates was much higher than that in the quarried aggregates, even when their mineralogical characteristics were similar. The magnitude of this effect is also much too large to be attributed to variations in CT resolution. The underlying basis for these differences in quantitative crack measurements can also be estimated through visual observation of CT images, such as those displayed in Figure 15. Clear, layered cracking is visible within river gravel greywacke stones; this is not present within the greywacke stones extracted from quarries. This is thought to result from the aggressive weathering process that river gravel is subjected to during its lifecycle prior to construction use. This indicates that the selection of high-quality aggregate based on mineralogical characteristics alone may be insufficient.
It is also clear from Figure 13 that even for stones from a single source and with a single mineral composition, a significant amount of variability in internal porosity and cracking is present. The variation between individual stones of a single type (such as rhyolite (GK1)) is often much larger than the average difference between two entirely different stone types (such as between rhyolite (GK1) and granite (GK4)). Thus, we recommend the use of large statistical samples to properly characterize each stone type for ASR sensitivity.

4. Discussion and Conclusions

This research clearly demonstrates the need for universal, automated, and consistent crack detection methods that allow the cross comparison of results from large quantities of CT-scan data from different sample types. A framework, called “virtual data fusion“, was developed that has the potential to successfully provide such a method. A partial implementation of this method in a custom program was developed for use in research focused on crack measurement in ASR-sensitive aggregates. Our results demonstrated the success of the program in effectively identifying crack-like structures and measuring their characteristics such as crack extension (relative surface area) and surface connectivity.
These results demonstrate the significant impact that the source of extraction can have on the characteristics of aggregates. Even for aggregates of the same mineral type, river gravels contain significantly higher levels of internal porosity and cracking than quarried stone. This is thought to result from the aggressive weathering process that river gravel is subjected to prior to its selection and use for construction. This indicates that the selection of high-quality aggregate based on mineralogical characteristics alone may be insufficient. It is also clear from these results that there is a significant amount of variability in internal porosity and cracking even for stones with the same mineralogical characteristics and extraction source. Thus, large statistical samples will be necessary to properly characterize each stone type for ASR sensitivity.

Author Contributions

Conceptualization, T.O.; methodology, T.O. and F.W.; software, T.O.; investigation, T.O. and F.W.; resources, T.O., F.W., and G.B.; writing—original draft preparation, T.O. and F.W.; writing—review and editing, T.O., F.W., and G.B.; visualization, T.O.; supervision, F.W. and G.B.; project administration, F.W.; funding acquisition, F.W. All authors have read and agreed to the published version of the manuscript.

Funding

This presentation is based on parts of the research project carried out at the request of the Federal Ministry of Transport and Digital Infrastructure, represented by the Federal Highway Research Institute, under research project No. 06.0108/2014/BRB. The author is solely responsible for the content. (Dieser Präsentation liegen Teile der im Auftrag des Bundesministeriums für Verkehr und digitale Infrastruktur, vertreten durch die Bundesanstalt für Straßenwesen, unter FE-Nr. 06.0108/2014/BRB durchgeführten Forschungsarbeit zugrunde. Die Verantwortung für den Inhalt liegt allein beim Autor).

Acknowledgments

The authors would like to thank Dietmar Meinel from the Bundesanstalt für Materialforschung und—prüfung (BAM, Federal Institute for Materials Research and Testing) for his guidance and support during the CT scanning and data analysis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A. Visualization of Individual Grains

A visual impression of the internal microstructure of the individual grains is provided by selected CT-based visualizations (Figure A1, Figure A2, Figure A3 and Figure A4). In these images, the solid material of the individual grains is shown semi-transparently. This allows a better spatial visualization of the cracks and pores emanating from the outer surface (i.e., open voids—colored red) and the cracks and pores not accessible from the outside (i.e., closed voids—colored yellow).
Figure A1. CT-visualizations for individual river gravel aggregate of type GK1 (the left and right images for each aggregate show views from the 0 and 90 degrees, respectively) (based on figures in [10,11]).
Figure A1. CT-visualizations for individual river gravel aggregate of type GK1 (the left and right images for each aggregate show views from the 0 and 90 degrees, respectively) (based on figures in [10,11]).
Materials 13 03921 g0a1aMaterials 13 03921 g0a1b
Figure A2. CT-visualizations for individual river gravel aggregate of type GK4 (the left and right images for each aggregate show views from the 0 and 90 degrees, respectively) (based on figures in [10,11]).
Figure A2. CT-visualizations for individual river gravel aggregate of type GK4 (the left and right images for each aggregate show views from the 0 and 90 degrees, respectively) (based on figures in [10,11]).
Materials 13 03921 g0a2aMaterials 13 03921 g0a2b
Figure A3. CT-visualizations for individual greywacke aggregate from quarried stone (GK2) and river gravel (GK4) (based on figures in [10,11]).
Figure A3. CT-visualizations for individual greywacke aggregate from quarried stone (GK2) and river gravel (GK4) (based on figures in [10,11]).
Materials 13 03921 g0a3
Figure A4. CT-visualizations for individual rhyolite aggregate from quarried stone (GK3) and river gravel (GK1) (based on figures in [10,11]).
Figure A4. CT-visualizations for individual rhyolite aggregate from quarried stone (GK3) and river gravel (GK1) (based on figures in [10,11]).
Materials 13 03921 g0a4
The strong fluctuation in the amount of open and closed voids in the individual quartz/quartzite grains of aggregate type GK1 is clearly visible in Figure A1. For the rhyolite of aggregate type GK1, Figure A1 includes a CT-based visualization of the single grain with the highest surface areas of open and closed voids. The broad spectrum of open and closed void surface characteristics for plutonite is demonstrated by the CT-based visualizations of the two individual grains shown in Figure A1.
In contrast to the GK1 aggregate, the surface characteristics of the individual grains of aggregate type GK4 tend to vary much less. In Figure A2, a CT-based visualization of the individual sandstone grain with the highest surface area of open voids of all examined sandstone grains is provided. As expected, this image shows a relatively high content of open voids and a low content of closed voids. The spatial arrangement of the open voids suggests a layered structure. The examined single grains of mudstone and greywacke also show a stratification due to the similar formation history. In the latter case, however, the closed voids predominate over open voids. The final CT-based visualization in Figure A2 is a single grain of granite and demonstrates that this type of rock can also have high contents of open voids.
As mentioned in the paper, the individual grains of the quarried stone have very small open and closed void surface areas compared to the river gravel grains. This is impressively documented by the CT-based visualizations for greywacke and rhyolite, which are shown in Figure A3 and Figure A4 and which compare individual grains originating from quarried stone and river gravel. Aside from those shown in Figure A3 and Figure A4, further visualizations of the quarried stone aggregate (GK2 and GK3) have not been included in the appendix. This is because these stone types have little to no visible porosity in the CT-based visualizations.

Appendix B. Tabulated Results of Cracking Analysis

Table A1. Individual grain results from the cracking analysis of stone category GK1 (river gravel; alkali sensitivity EIII-S).
Table A1. Individual grain results from the cracking analysis of stone category GK1 (river gravel; alkali sensitivity EIII-S).
Stone TypeMass of the Individual Grain [g]CT Scan NumberVoxel Size [µm]Surface Areas Measured Using CT [m2]
Closed PorosityOuter SurfaceOpen PorosityOuter Surface + Open Porosity
Absolute [m2]Relative [m2/g]
Quartz/Quartzite3.37703016.55.69 × 10−42.17 × 10−32.20 × 10−52.18 × 10−36.47 × 10−4
4.4970311.42 × 10−41.10 × 10−36.20 × 10−51.12 × 10−32.49 × 10−4
2.8970323.85 × 10−51.44 × 10−39.65 × 10−51.47 × 10−35.09 × 10−4
5.6370333.94 × 10−69.50 × 10−41.43 × 10−69.50 × 10−41.69 × 10−4
6.2670341.10 × 10−31.55 × 10−31.46 × 10−41.57 × 10−32.51 × 10−4
8.0470359.13 × 10−42.39 × 10−33.24 × 10−32.79 × 10−33.47 × 10−4
7.3370362.62 × 10−41.95 × 10−31.88 × 10−41.97 × 10−32.69 × 10−4
8.3270372.33 × 10−32.38 × 10−38.53 × 10−42.48 × 10−32.98 × 10−4
6.5070383.64 × 10−51.84 × 10−32.89 × 10−61.84 × 10−32.83 × 10−4
11.1770393.10 × 10−42.09 × 10−33.71 × 10−42.13 × 10−31.91 × 10−4
Mean5.71 × 10−41.79 × 10−34.98 × 10−41.85 × 10−33.21 × 10−4
Standard Deviation6.88 × 10−44.82 × 10−49.46 × 10−45.52 × 10−41.40 × 10−4
Rhyolite0.77688011.20.005.18 × 10−40.005.18 × 10−46.73 × 10−4
0.7468818.60 × 10−64.87 × 10−41.71 × 10−64.88 × 10−46.60 × 10−4
1.2768821.59 × 10−47.50 × 10−47.07 × 10−58.21 × 10−46.47 × 10−4
1.2468836.55 × 10−49.78 × 10−44.87 × 10−41.46 × 10−31.18 × 10−3
1.9668849.57 × 10−68.32 × 10−41.49 × 10−58.47 × 10−44.32 × 10−4
1.1668851.30 × 10−65.65 × 10−41.33 × 10−65.67 × 10−44.89 × 10−4
2.9068868.63 × 10−51.21 × 10−33.28 × 10−51.25 × 10−34.30 × 10−4
2.7468871.79 × 10−41.20 × 10−33.35 × 10−51.24 × 10−34.52 × 10−4
2.8568883.67 × 10−51.03 × 10−31.80 × 10−51.05 × 10−33.68 × 10−4
5.7068896.38 × 10−51.62 × 10−33.81 × 10−61.62 × 10−32.85 × 10−4
Mean1.20 × 10−49.20 × 10−46.64 × 10−59.86 × 10−45.62 × 10−4
Standard Deviation1.89 × 10−43.44 × 10−41.42 × 10−43.82 × 10−42.41 × 10−4
Plutonite2.43697011.21.21 × 10−49.24 × 10−41.17 × 10−59.35 × 10−43.85 × 10−4
1.5069711.57 × 10−49.03 × 10−43.64 × 10−41.27 × 10−38.44 × 10−4
1.6269724.53 × 10−56.74 × 10−49.49 × 10−66.84 × 10−44.22 × 10−4
2.0969732.12 × 10−44.72 × 10−39.71 × 10−45.69 × 10−32.72 × 10−3
2.6769741.23 × 10−31.07 × 10−32.26 × 10−41.30 × 10−34.87 × 10−4
1.9669753.46 × 10−48.53 × 10−41.02 × 10−49.54 × 10−44.87 × 10−4
5.5869761.66 × 10−33.43 × 10−32.35 × 10−35.79 × 10−31.04 × 10−3
2.4769773.58 × 10−41.00 × 10−37.83 × 10−51.08 × 10−34.38 × 10−4
4.0869787.40 × 10−41.65 × 10−32.67 × 10−41.91 × 10−34.69 × 10−4
4.4469790.001.24 × 10−39.91 × 10−81.24 × 10−32.79 × 10−4
Mean4.86 × 10−41.65 × 10−34.38 × 10−42.08 × 10−37.57 × 10−4
Standard Deviation5.27 × 10−41.27 × 10−36.95 × 10−41.85 × 10−36.89 × 10−4
Table A2. Individual grain results from the cracking analysis of stone category GK2 (quarried stone; alkali sensitivity EIII-S).
Table A2. Individual grain results from the cracking analysis of stone category GK2 (quarried stone; alkali sensitivity EIII-S).
Stone TypeMass of the Individual Grain [g]CT Scan NumberVoxel Size [µm]Surface Areas Measured Using CT [m2]
Closed PorosityOuter SurfaceOpen PorosityOuter Surface + Open Porosity
Absolute [m2]Relative [m2/g]
Greywacke1.03686011.22.94 × 10−75.31 × 10−42.67 × 10−65.34 × 10−45.18 × 10−4
1.1268615.87 × 10−86.23 × 10−40.006.23 × 10−45.56 × 10−4
0.6668627.67 × 10−74.74 × 10−48.48 × 10−74.74 × 10−47.19 × 10−4
2.2368632.59 × 10−79.14 × 10−45.49 × 10−79.15 × 10−44.10 × 10−4
2.0468640.009.55 × 10−43.03 × 10−79.55 × 10−44.68 × 10−4
1.968650.009.78 × 10−41.72 × 10−69.80 × 10−45.16 × 10−4
1.3968661.14 × 10−77.21 × 10−40.007.21 × 10−45.19 × 10−4
1.0868670.007.26 × 10−41.48 × 10−77.26 × 10−46.72 × 10−4
2.1168689.99 × 10−81.03 × 10−36.28 × 10−71.03 × 10−34.87 × 10−4
3.9868694.92 × 10−81.61 × 10−31.21 × 10−61.61 × 10−34.04 × 10−4
Mean1.64 × 10−78.56 × 10−48.07 × 10−78.57 × 10−45.27 × 10−4
Standard Deviation2.24 × 10−73.11 × 10−48.09 × 10−73.11 × 10−49.63 × 10−5
Table A3. Individual grain results from the cracking analysis of stone category GK3 (quarried stone; alkali sensitivity EI-S).
Table A3. Individual grain results from the cracking analysis of stone category GK3 (quarried stone; alkali sensitivity EI-S).
Stone TypeMass of the Individual Grain [g]CT Scan NumberVoxel Size [µm]Surface Areas Measured Using CT [m2]
Closed PorosityOuter SurfaceOpen PorosityOuter Surface + Open Porosity
Absolute [m2]Relative [m2/g]
Rhyolite1.53689016.55.14 × 10−69.00 × 10−41.86 × 10−59.18 × 10−46.00 × 10−4
1.6668914.52 × 10−79.12 × 10−48.54 × 10−69.21 × 10−45.55 × 10−4
2.6868926.18 × 10−61.32 × 10−34.16 × 10−51.36 × 10−35.07 × 10−4
2.5668936.92 × 10−61.13 × 10−32.93 × 10−51.16 × 10−34.52 × 10−4
2.8968941.36 × 10−71.30 × 10−31.24 × 10−51.32 × 10−34.56 × 10−4
3.0868954.42 × 10−61.24 × 10−36.88 × 10−61.25 × 10−34.05 × 10−4
5.1868967.97 × 10−71.97 × 10−36.97 × 10−61.98 × 10−33.82 × 10−4
3.4468979.38 × 10−61.42 × 10−32.21 × 10−51.44 × 10−34.19 × 10−4
5.0768981.44 × 10−61.96 × 10−36.57 × 10−61.96 × 10−33.87 × 10−4
3.4768998.69 × 10−71.50 × 10−33.59 × 10−61.50 × 10−34.32 × 10−4
Mean3.57 × 10−61.36 × 10−31.57 × 10−51.38 × 10−34.59 × 10−4
Standard Deviation3.10 × 10−63.52 × 10−41.16 × 10−53.48 × 10−46.92 × 10−5
Table A4. Individual grain results from the cracking analysis of stone category GK4 (river gravel; alkali sensitivity EI-S).
Table A4. Individual grain results from the cracking analysis of stone category GK4 (river gravel; alkali sensitivity EI-S).
Stone TypeMass of the Individual Grain [g]CT Scan NumberVoxel Size [µm]Surface Areas Measured Using CT [m2]
Closed PorosityOuter SurfaceOpen PorosityOuter Surface + Open Porosity
Absolute [m2]Relative [m2/g]
Sandstone0.84698021.10.005.67 × 10−40.005.67 × 10−46.75 × 10−4
1.9569810.007.24 × 10−40.007.24 × 10−43.71 × 10−4
1.3569821.06 × 10−56.95 × 10−48.69 × 10−67.03 × 10−45.21 × 10−4
3.6369835.04 × 10−61.07 × 10−32.26 × 10−71.07 × 10−32.96 × 10−4
3.3569847.44 × 10−51.28 × 10−35.34 × 10−41.82 × 10−35.42 × 10−4
4.1969857.09 × 10−51.56 × 10−31.09 × 10−41.67 × 10−34.00 × 10−4
4.7869861.01 × 10−31.42 × 10−34.23 × 10−51.47 × 10−33.07 × 10−4
5.2869873.25 × 10−41.66 × 10−32.01 × 10−51.68 × 10−33.18 × 10−4
6.2669881.80 × 10−41.91 × 10−34.91 × 10−51.95 × 10−33.12 × 10−4
8.369893.80 × 10−42.20 × 10−31.59 × 10−52.21 × 10−32.67 × 10−4
Mean2.05 × 10−41.31 × 10−37.79 × 10−51.39 × 10−34.01 × 10−4
Standard Deviation2.97 × 10−45.17 × 10−41.55 × 10−45.52 × 10−41.28 × 10−4
Mudstone3.6696021.11.72 × 10−41.30 × 10−32.56 × 10−51.33 × 10−33.68 × 10−4
1.5969614.28 × 10−77.14 × 10−49.34 × 10−77.15 × 10−44.49 × 10−4
1.6169625.36 × 10−59.01 × 10−45.17 × 10−41.42 × 10−38.81 × 10−4
3.2769630.001.15 × 10−31.26 × 10−61.15 × 10−33.51 × 10−4
1.1369640.001.45 × 10−30.001.45 × 10−31.28 × 10−3
2.7869657.04 × 10−41.02 × 10−34.66 × 10−41.49 × 10−35.36 × 10−4
2.3969664.35 × 10−68.54 × 10−40.008.54 × 10−43.57 × 10−4
2.6369672.16 × 10−41.12 × 10−34.55 × 10−51.16 × 10−34.42 × 10−4
4.3469680.001.70 × 10−38.01 × 10−71.70 × 10−33.93 × 10−4
2.5169691.13 × 10−51.06 × 10−31.56 × 10−61.06 × 10−34.23 × 10−4
Mean1.16 × 10−41.13 × 10−31.06 × 10−41.23 × 10−35.48 × 10−4
Standard Deviation2.10 × 10−42.79 × 10−41.94 × 10−42.88 × 10−42.86 × 10−4
Greywacke8.19687016.55.61 × 10−42.33 × 10−31.25 × 10−42.46 × 10−33.00 × 10−4
5.2868716.17 × 10−51.99 × 10−33.93 × 10−52.03 × 10−33.84 × 10−4
8.6768725.56 × 10−42.70 × 10−33.04 × 10−43.01 × 10−33.47 × 10−4
5.9968731.71 × 10−63.13 × 10−33.88 × 10−53.17 × 10−35.29 × 10−4
4.8968741.68 × 10−31.69 × 10−34.92 × 10−42.18 × 10−34.46 × 10−4
2.6168753.26 × 10−61.43 × 10−34.02 × 10−61.44 × 10−35.50 × 10−4
6.368767.10 × 10−42.16 × 10−31.11 × 10−42.27 × 10−33.61 × 10−4
3.8668778.41 × 10−51.51 × 10−37.54 × 10−61.51 × 10−33.92 × 10−4
3.0968784.12 × 10−42.31 × 10−32.89 × 10−42.60 × 10−38.41 × 10−4
2.6368799.88 × 10−61.24 × 10−32.88 × 10−61.24 × 10−34.72 × 10−4
Mean4.08 × 10−42.05 × 10−31.41 × 10−42.19 × 10−34.62 × 10−4
Standard Deviation4.96 × 10−45.68 × 10−41.58 × 10−46.19 × 10−41.47 × 10−4
Granite1.53695011.22.14 × 10−58.15 × 10−46.23 × 10−58.77 × 10−45.74 × 10−4
0.7269511.58 × 10−65.08 × 10−41.41 × 10−55.22 × 10−47.26 × 10−4
1.0969523.63 × 10−56.80 × 10−41.84 × 10−48.64 × 10−47.93 × 10−4
0.9869531.17 × 10−55.76 × 10−45.53 × 10−56.31 × 10−46.44 × 10−4
1.9269541.01 × 10−68.71 × 10−41.59 × 10−58.87 × 10−44.62 × 10−4
1.7669557.23 × 10−59.39 × 10−46.90 × 10−51.01 × 10−35.73 × 10−4
3.1769561.36 × 10−41.43 × 10−33.65 × 10−41.80 × 10−35.68 × 10−4
2.7769571.37 × 10−41.28 × 10−35.26 × 10−41.81 × 10−36.53 × 10−4
2.6969584.95 × 10−51.62 × 10−35.91 × 10−51.67 × 10−36.23 × 10−4
6.1569591.09 × 10−41.87 × 10−39.91 × 10−81.87 × 10−33.04 × 10−4
Mean5.76 × 10−51.06 × 10−31.35 × 10−41.19 × 10−35.92 × 10−4
Standard Deviation5.06 × 10−54.42 × 10−41.66 × 10−45.04 × 10−41.29 × 10−4

References

  1. Stanton, T. Expansion of Concrete Through Reaction Between Cement and Aggregate. Proc. Am. Soc. Civ. Eng. 1940, 66, 1781–1812. [Google Scholar]
  2. Powers, T.; Steinour, H. An Interpretation of Some Published Researches on the Alkali-Aggregate Reaction Part 1—The Chemical Reactions and Mechanism of Expansion. J. Am. Concr. I 1955, 26, 497–516. [Google Scholar]
  3. Powers, T.; Steinour, H. An Interpretation of Some Published Researches on the Alkali-Aggregate Reaction Part 2—A Hypothesis Concerning Safe and Unsafe Reactions with Reactive Silica in Concrete. J. Am. Concr. I 1955, 51, 785–812. [Google Scholar]
  4. Locher, F.; Sprung, S. Ursache und Wirkungsweise der Alkalireaktion; Betontechnische Berichte; Forschungsinstitut der Zementindustrie: Düsseldorf, Germany, 1973; pp. 101–123. [Google Scholar]
  5. Chatterji, S.; Thaulow, N.; Jensen, A. Studies of alkali-silica-reaction. Part 5. Verification of a newly proposed reaction mechanism. Cem. Concr. Res. 1989, 19, 177–183. [Google Scholar] [CrossRef]
  6. Chatterji, S.; Thaulow, N.; Jensen, A. Studies of alkali-silica-reaction. Part 6. Practical implications of a proposed reaction mechanism. Cem. Concr. Res. 1988, 18, 363–366. [Google Scholar] [CrossRef]
  7. Stark, J.; Erfurt, D.; Freyburg, E.; Giebson, C.; Seyfarth, K.; Wicht, B. Alkali-Kieselsäure-Reaktion; F. A. Finger-Instituts für Baustoffkunde: Weimar, Germany, 2008. [Google Scholar]
  8. Thomas, M.D.A.; Fournier, B.; Folliard, K.J. Alkali-Aggregate Reactivity (AAR) Facts Book; U.S. Department of Transportation, Federal Highway Administration (FHWA): Washington, DC, USA, 2013.
  9. Qi, Y.; Ziyun, W. Study of expansion mechanism of A.S.R. using sol-gel expansion method. In Proceedings of the 12th ICAAR, Beijing, China, 15–19 October 2004; pp. 226–229. [Google Scholar]
  10. Weise, F.; Kositz, M.; Oesch, T.; Huenger, K.-J.; Wilsch, G.; Sigmund, S. Analyse des Gefügeabhängigen Löslichkeitsverhaltens Potenziell AKR-Empfindlicher Gesteinskörnungen; Bundesanstalt für Straßenwesen (BASt), Schlussbericht FE 06.0108/2014/BRB: Bergisch Gladbach, Germany, 2019. [Google Scholar]
  11. Oesch, T.; Weise, F.; Marx, H.; Kositz, M.; Huenger, K.-J. Analysis of the porosity of alkali-sensitive aggregates for the assessment of microstructure-dependent solubility in the context of ASR (ACCEPTED). In Proceedings of the 4th International RILEM Conference on Microstructure Related Durability of Cementitious Composites, Den Haag, The Netherlands, 28–30 April 2021. [Google Scholar]
  12. Buzug, T.M. Computed Tomography. In Springer Handbook of Medical Technology; Kramme, R., Hoffmann, K.-P., Pozos, R.S., Eds.; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar] [CrossRef]
  13. Kalender, W.A. Computed Tomography: Fundamentals, System Technology, Image Quality, Applications, 3rd ed.; Publicis Publishing: Erlangen, Germany, 2011. [Google Scholar]
  14. Willemink, M.J.; Noël, P.B. The evolution of image reconstruction for CT-from filtered back projection to artificial intelligence. Eur. Radiol. 2019, 29, 2185–2195. [Google Scholar] [CrossRef] [Green Version]
  15. du Plessis, A.; le Roux, S.G.; Guelpa, A. Comparison of medical and industrial X-ray computed tomography for non-destructive testing. Case Stud. Nondestruct. Test. Eval. 2016, 6, 17–25. [Google Scholar] [CrossRef] [Green Version]
  16. Martz, H.E.; Schneberk, D.J.; Roberson, G.P.; Monteiro, P.J.M. Computerized-Tomography Analysis of Reinforced-Concrete. ACI Mater. J. 1993, 90, 259–264. [Google Scholar]
  17. Morgan, I.L.; Ellinger, H.; Klinksiek, R.; Thompson, J.N. Examination of Concrete by Computerized-Tomography. J. Am. Concr. I 1980, 77, 23–27. [Google Scholar]
  18. Flannery, B.P.; Deckman, H.W.; Roberge, W.G.; D’Amico, K.L. Three-Dimensional X-ray Microtomography. Science 1987, 237, 1439–1444. [Google Scholar] [CrossRef]
  19. Feldkamp, L.A.; Davis, L.C.; Kress, J.W. Practical cone-beam algorithm. J. Opt. Soc. Am. A 1984, 1, 612–619. [Google Scholar] [CrossRef] [Green Version]
  20. Landis, E.N.; Nagy, E.N.; Keane, D.T. Microtomographic measurements of internal damage in portland-cement-based composites. J. Aerosp. Eng. 1997, 10, 2–6. [Google Scholar] [CrossRef]
  21. Asahina, D.; Landis, E.N.; Bolander, J.E. Modeling of phase interfaces during pre-critical crack growth in concrete. Cem. Concr. Compos. 2011, 33, 966–977. [Google Scholar] [CrossRef]
  22. Poinard, C.; Piotrowska, E.; Malecot, Y.; Daudeville, L.; Landis, E.N. Compression triaxial behavior of concrete: The role of the mesostructure by analysis of X-ray tomographic images. Eur. J. Environ. Civ. Eng. 2012, 16, S115–S136. [Google Scholar] [CrossRef]
  23. Oesch, T.; Landis, E.; Kuchma, D. Conventional Concrete and UHPC Performance-Damage Relationships Identified Using Computed Tomography. J. Eng. Mechan. 2016, 142. [Google Scholar] [CrossRef]
  24. Oesch, T. In-Situ CT Investigation of Pull-Out Failure for Reinforcing Bars Embedded in Conventional and High-Performance Concretes. In Proceedings of the 6th Conference on Industrial Computed Tomography (ICT), Wels, Austria, 9–12 February 2016. [Google Scholar]
  25. Paetsch, O.; Baum, D.; Prohaska, S.; Ehrig, K.; Meinel, D.; Ebell, G. 3D Corrosion Detection in Time-dependent CT Images of Concrete. In Proceedings of the Digital Industrial Radiology and Computed Tomography Conference, Ghent, Belgium, 22–25 June 2015; p. 10. [Google Scholar]
  26. Yang, L.; Zhang, Y.S.; Liu, Z.Y.; Zhao, P.; Liu, C. In-situ tracking of water transport in cement paste using X-ray computed tomography combined with CsCl enhancing. Mater. Lett. 2015, 160, 381–383. [Google Scholar] [CrossRef]
  27. Boone, M.A.; De Kock, T.; Bultreys, T.; De Schutter, G.; Vontobel, P.; Van Hoorebeke, L.; Cnudde, V. 3D mapping of water in oolithic limestone at atmospheric and vacuum saturation using X-ray micro-CT differential imaging. Mater. Charact. 2014, 97, 150–160. [Google Scholar] [CrossRef]
  28. Oesch, T.; Weise, F.; Meinel, D.; Gollwitzer, G. Quantitative In-situ Analysis of Water Transport in Concrete Completed Using X-ray Computed Tomography. Trans. Porous Med. 2019, 127, 371–389. [Google Scholar] [CrossRef]
  29. Powierza, B.; Stelzner, L.; Oesch, T.; Gollwitzer, C.; Weise, F.; Bruno, G. Water Migration in One-Side Heated Concrete: 4D In-Situ CT Monitoring of the Moisture-Clog-Effect. J. Nondestruct. Eval. 2019, 38, 15. [Google Scholar] [CrossRef]
  30. Bažant, Z.P.; Planas, J. Fracture and Size Effect in Concrete and Other Quasibrittle Materials; CRC Press: Boca Raton, FL, USA, 1998; p. 640. [Google Scholar]
  31. Oesch, T. Investigation of Fiber and Cracking Behavior for Conventional and Ultra-High Performance Concretes Using X-ray Computed Tomography. Ph.D. Thesis, University of Illinois, Urbana, IL, USA, 2015. [Google Scholar]
  32. Ehrig, K.; Goebbels, J.; Meinel, D.; Paetsch, O.; Prohaska, S.; Zobel, V. Comparison of Crack Detection Methods for Analyzing Damage Processes in Concrete with Computed Tomography. In Proceedings of the International Symposium on Digital Industrial Radiology and Computed Tomography, Berlin, Germany, 20–22 June 2011. [Google Scholar]
  33. Lebbink, M.N.; Geerts, W.J.C.; van der Krift, T.P.; Bouwhuis, M.; Hertzberger, L.O.; Verkleij, A.J.; Koster, A.J. Template matching as a tool for annotation of tomograms of stained biological structures. J. Struct. Biol. 2007, 158, 327–335. [Google Scholar] [CrossRef]
  34. Paetsch, O. Possibilities and Limitations of Automatic Feature Extraction shown by the Example of Crack Detection in 3D-CT Images of Concrete Specimen. In Proceedings of the 9th Conference on Industrial Computed Tomography, Padova, Italy, 13–15 February 2019. [Google Scholar]
  35. Sato, Y.; Westin, C.F.; Bhalerao, A.; Nakajima, S.; Shiraga, N.; Tamura, S.; Kikinis, R. Tissue classification based on 3D local intensity structures for volume rendering. IEEE Trans. Vis. Comput. Gr. 2000, 6, 160–180. [Google Scholar] [CrossRef] [Green Version]
  36. Watanabe, S.; Ohtake, Y.; Nagai, Y.; Suzuki, H. Detection of Narrow Gaps Using Hessian Eigenvalues for Shape Segmentation of a CT Volume of Assembled Parts. In Proceedings of the 8th Conference on Industrial Computed Tomography, Wels, Austria, 6–9 February 2018. [Google Scholar]
  37. Picard, D.; Lauzon-Gauthier, J.; Duchesne, C.; Alamdari, H.; Fafard, M.; Ziegler, D.P. Crack Detection Method Applied to 3D Computed Tomography Images of Baked Carbon Anodes. Metals 2016, 6, 272. [Google Scholar] [CrossRef]
  38. Cinar, A.F.; Hollis, D.; Tomlinson, R.A.; Marrow, T.J.; Mostafavi, M. Application of 3D phase congruency in crack identification within materials. In Proceedings of the 11th International Conference on Advances in Experimental Mechanics, Exeter, UK, 5–7 September 2016. [Google Scholar]
  39. Cinar, A.F.; Barhli, S.M.; Hollis, D.; Flansbjer, M.; Tomlinson, R.A.; Marrow, T.J.; Mostafavi, M. An autonomous surface discontinuity detection and quantification method by digital image correlation and phase congruency. Opt. Laser Eng. 2017, 96, 94–106. [Google Scholar] [CrossRef]
  40. Barhli, S.M.; Saucedo-Mora, L.; Jordan, M.S.L.; Cinar, A.F.; Reinhard, C.; Mostafavi, M.; Marrow, T.J. Synchrotron X-ray characterization of crack strain fields in polygranular graphite. Carbon 2017, 124, 357–371. [Google Scholar] [CrossRef] [Green Version]
  41. Powierza, B.; Gollwitzer, C.; Wolgast, D.; Staude, A.; Bruno, G. Fully experiment-based evaluation of few digital volume correlation techniques. Rev. Sci. Instrum. 2019, 90, 115105. [Google Scholar] [CrossRef] [PubMed]
  42. Kovesi, P. Phase congruency: A low-level image invariant. Psychol. Res. 2000, 64, 136–148. [Google Scholar] [CrossRef] [PubMed]
  43. Bhandarkar, S.M.; Luo, X.Z.; Daniels, R.; Tollner, E.W. Detection of cracks in computer tomography images of logs. Pattern Recogn. Lett. 2005, 26, 2282–2294. [Google Scholar] [CrossRef]
  44. DAfStb. Vorbeugende Maßnahmen Gegen Schädigende Alkalireaktion im Beton (Alkali-Richtlinie); Deutscher Ausschuss für Stahlbeton; Beuth Verlag GmbH: Berlin, Germany, 2013. [Google Scholar]
  45. DIN. Tests for Geometrical Properties of Aggregates—Part 1: Determination of Particle Size Distribution—Sieving Method; German Version EN 933-1:2012; Deutsches Institut für Normung; Beuth Verlag GmbH: Berlin, Germany, 2012. [Google Scholar]
  46. Latief, F.D.; Fauzi, U.; Irayani, Z.; Dougherty, G. The effect of X-ray micro computed tomography image resolution on flow properties of porous rocks. J. Microsc. 2017, 266, 69–88. [Google Scholar] [CrossRef]
  47. Krumm, A.; Means, B.; Bienkowski, M. Learning Analytics Goes to School: A Collaborative Approach to Improving Education; Routledge, Taylor & Francis Group: New York, NY, USA, 2018. [Google Scholar]
  48. Heideklang, R.; Shokouhi, P. Multi-sensor image fusion at signal level for improved near-surface crack detection. Ndt E Int. 2015, 71, 16–22. [Google Scholar] [CrossRef]
  49. The Mathworks. MATLAB, R2017b; The Mathworks: Natick, MA, USA, 2018. [Google Scholar]
  50. Young, I.T.; Gerbrands, J.J.; van Vliet, L.J. Fundamentals of Image Processing, 2.2 ed.; Delft University of Technology: Delft, The Netherlands, 1998. [Google Scholar]
  51. Zack, G.W.; Rogers, W.E.; Latt, S.A. Automatic measurement of sister chromatid exchange frequency. J. Histochem. Cytochem. 1977, 25, 741–753. [Google Scholar] [CrossRef]
  52. Pringle, K.K. Visual Perception by a computer. In Automatic Interpretation and Classification of Images; Grasselli, A., Ed.; Academic Press: New York, NY, USA, 1969; pp. 277–284. [Google Scholar]
  53. Brunauer, S.; Emmett, P.H.; Teller, E. Adsorption of gases in multimolecular layers. J. Am. Chem. Soc. 1938, 60, 309–319. [Google Scholar] [CrossRef]
  54. Langmuir, I. The Adsorption of Gases on Plane Surfaces of Glass, Mica and Platinum. J. Am. Chem. Soc. 1918, 40, 1361–1403. [Google Scholar] [CrossRef] [Green Version]
  55. Cooper, R.C.; Bruno, G.; Onel, Y.; Lange, A.; Watkins, T.R.; Shyam, A. Young’s modulus and Poisson’s ratio changes due to machining in porous microcracked cordierite. J. Mater. Sci. 2016, 51, 9749–9760. [Google Scholar] [CrossRef]
  56. Mishurova, T.; Artzt, K.; Haubrich, J.; Requena, G.; Bruno, G. New aspects about the search for the most relevant parameters optimizing SLM materials. Addit. Manuf. 2019, 25, 325–334. [Google Scholar] [CrossRef]
Figure 1. Photographs of typical individual grains from each of the stone categories.
Figure 1. Photographs of typical individual grains from each of the stone categories.
Materials 13 03921 g001
Figure 2. Virtual data fusion web linking cracks (black) with their characteristics (dark grey) and image analysis techniques for identifying objects with those characteristics (light grey).
Figure 2. Virtual data fusion web linking cracks (black) with their characteristics (dark grey) and image analysis techniques for identifying objects with those characteristics (light grey).
Materials 13 03921 g002
Figure 3. Step-by-step isolation of cracking by means of virtual data fusion. CT image of a concrete sample (a), voids identified using a greyscale threshold (b), voids identified in (b) added to interfacial zones identified using a gradient-based analysis (c), removal of small connected components (noise and isolated pores) from image (d), and elimination of components with small specific surface areas from image (e). Raw data courtesy of U.S. Army Engineer Research and Development Center (ERDC).
Figure 3. Step-by-step isolation of cracking by means of virtual data fusion. CT image of a concrete sample (a), voids identified using a greyscale threshold (b), voids identified in (b) added to interfacial zones identified using a gradient-based analysis (c), removal of small connected components (noise and isolated pores) from image (d), and elimination of components with small specific surface areas from image (e). Raw data courtesy of U.S. Army Engineer Research and Development Center (ERDC).
Materials 13 03921 g003
Figure 4. Automated threshold selection using the triangle approach (based on figure in [10]).
Figure 4. Automated threshold selection using the triangle approach (based on figure in [10]).
Materials 13 03921 g004
Figure 5. Binarization of greyscale image using the conservative and over-encompassing thresholds identified using the automated selection method (based on figure in [10]).
Figure 5. Binarization of greyscale image using the conservative and over-encompassing thresholds identified using the automated selection method (based on figure in [10]).
Materials 13 03921 g005
Figure 6. Generation of the gradient-magnitude image (based on figure in [10]).
Figure 6. Generation of the gradient-magnitude image (based on figure in [10]).
Materials 13 03921 g006
Figure 7. Automated gradient threshold selection using the triangle approach (based on figure in [10]).
Figure 7. Automated gradient threshold selection using the triangle approach (based on figure in [10]).
Materials 13 03921 g007
Figure 8. Binarization of the gradient image (based on figure in [10]).
Figure 8. Binarization of the gradient image (based on figure in [10]).
Materials 13 03921 g008
Figure 9. Generation of an image containing the outer edges of only low-density objects (based on figure in [10]).
Figure 9. Generation of an image containing the outer edges of only low-density objects (based on figure in [10]).
Materials 13 03921 g009
Figure 10. Fusion of images resulting from Modules 1 and 2 (based on figure in [10]).
Figure 10. Fusion of images resulting from Modules 1 and 2 (based on figure in [10]).
Materials 13 03921 g010
Figure 11. Removal of noise from binary crack and void image (based on figure in [10]).
Figure 11. Removal of noise from binary crack and void image (based on figure in [10]).
Materials 13 03921 g011
Figure 12. Flowchart of crack-detection process (based on figure in [10]).
Figure 12. Flowchart of crack-detection process (based on figure in [10]).
Materials 13 03921 g012
Figure 13. Measured surface area for river-gravel type aggregates (based on figure in [10]).
Figure 13. Measured surface area for river-gravel type aggregates (based on figure in [10]).
Materials 13 03921 g013
Figure 14. Surface area measurement comparison between quarried-stone (GK2 and GK3) and river-gravel (GK1 and GK4) aggregates of the same mineral types (based on figures in [10,11]).
Figure 14. Surface area measurement comparison between quarried-stone (GK2 and GK3) and river-gravel (GK1 and GK4) aggregates of the same mineral types (based on figures in [10,11]).
Materials 13 03921 g014
Figure 15. Example images of cracking within a quarried greywacke aggregate (a) and a river-gravel type greywacke aggregate (b). Red pixels denote externally accessible cracks/pores and yellow pixels denote externally inaccessible (closed) cracks/pores (reproduction of figures from [10,11]).
Figure 15. Example images of cracking within a quarried greywacke aggregate (a) and a river-gravel type greywacke aggregate (b). Red pixels denote externally accessible cracks/pores and yellow pixels denote externally inaccessible (closed) cracks/pores (reproduction of figures from [10,11]).
Materials 13 03921 g015
Table 1. Categories of stone selected for cracking analysis.
Table 1. Categories of stone selected for cracking analysis.
CategoryStone TypeAlkali Sensitivity [40 °C-BV]
GK1River GravelEIII-S
GK2Quarried Stone (Greywacke)EIII-S
GK3Quarried Stone (Rhyolite)EI-S
GK4River GravelEI-S

Share and Cite

MDPI and ACS Style

Oesch, T.; Weise, F.; Bruno, G. Detection and Quantification of Cracking in Concrete Aggregate through Virtual Data Fusion of X-Ray Computed Tomography Images. Materials 2020, 13, 3921. https://doi.org/10.3390/ma13183921

AMA Style

Oesch T, Weise F, Bruno G. Detection and Quantification of Cracking in Concrete Aggregate through Virtual Data Fusion of X-Ray Computed Tomography Images. Materials. 2020; 13(18):3921. https://doi.org/10.3390/ma13183921

Chicago/Turabian Style

Oesch, Tyler, Frank Weise, and Giovanni Bruno. 2020. "Detection and Quantification of Cracking in Concrete Aggregate through Virtual Data Fusion of X-Ray Computed Tomography Images" Materials 13, no. 18: 3921. https://doi.org/10.3390/ma13183921

APA Style

Oesch, T., Weise, F., & Bruno, G. (2020). Detection and Quantification of Cracking in Concrete Aggregate through Virtual Data Fusion of X-Ray Computed Tomography Images. Materials, 13(18), 3921. https://doi.org/10.3390/ma13183921

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop