1. Introduction
To understand and describe the structural behavior of concrete, in-depth knowledge of the damage mechanisms is required. A crack formed in concrete is the result of a damage process, which starts with the formation and merging of micro-cracks. The evaluation of the resulting crack propagation velocity can provide valuable information on the process and rate of the damage. Such information is crucial to explaining the behavior of concrete at high loading rates. Several researchers have tried to explain the rate-dependent behavior of concrete based on the combined effects of high loading rate and low crack propagation velocity [
1,
2,
3]. Given that crack formation does not occur instantaneously, but rather requires a certain amount of time for the development of damage, it can be deduced that under high loading rates, a strength-enhancing effect can emerge. Thus, a partial cause of the so-called strain rate effect [
4,
5,
6,
7] is likely attributable to the limited crack propagation velocity. The theoretical maximum of crack propagation velocity is the Rayleigh or surface wave velocity, that is the velocity at which energy can be transported along a surface. In concrete, the Rayleigh wave propagation velocity is approximately 2200 m/s [
8]. However, experimental investigations by various researchers have been able to prove only significantly lower crack propagation velocities.
Bhargava/Rehnström [
9] estimated crack propagation velocities of about 180 m/s based on film recordings taken in blast tests on concrete prisms. Mindess and Bentur resp. Mindess [
10,
11] also used camera images to evaluate impact-loaded beams and determined an approximate crack propagation velocity of about 100 m/s by comparing the crack progress interactively, frame-to-frame. Curbach [
12] used notched concrete specimens with applied conductive lacquer barriers, which were destroyed one after the other depending on the crack propagation progress. Based on the measurement data, he postulated an upper limit of the crack velocity of about 500 m/s in concrete. Zhang et al. [
13] investigated the crack velocity in high-strength concretes under different loading rates using strain gauges. They observed an increase in crack propagation velocity with the increasing loading rate. At low loading rates (10
−4 to 10
1 mm/s), significant differences in the velocity were found during cracking. Late-stage velocity was observed to be one order of magnitude higher when compared to that observed in earlier stages. However, under high loading rates (10
2 to 10
3 mm/s), a steadier crack propagation velocity was observed, with velocities in the range of 200 to 400 m/s. Acoustic emission [
14] and digital image correlation (DIC) [
15] can be considered as the most recent methods for measuring crack propagation speed. Obtaining the temporal position of the crack tip using high-speed camera shots in combination with displacement field analysis using a commercial DIC software [
16] shows promise in the application of such novel methods. A crack propagation velocity of about 1300 m/s was obtained in spalling experiments on notched concrete specimens.
An obvious advantage of high-speed, camera-based techniques in detecting cracks and determining their crack propagation velocity is the fact that the entire surface of a sample can be analyzed simultaneously, thus enabling in turn the simultaneous observation of multiple cracks. However, interactive image sequence processing is rather tedious. Hence, automated techniques are required for the efficient processing of the experimental data. Reliable and precise image processing routines should also exclude observer errors and thus contribute to increasing the reproducibility of results. The article at hand presents an approach to crack velocity measurement based on automated high-speed camera image processing. A high-speed camera was used to monitor dynamic tests on compact-tension specimens at a frame rate of 160,000 images per second. The images were processed using an automated crack-detection and measurement procedure able to detect cracks with a width on the order of 0.1 pixel. In the following, the image sequence processing chain for crack detection is outlined. Then, the detection of crack tips (and their propagation velocities) from the image processing results are described. Finally, the results of applying the methods developed for high-speed camera image processing will be compared to those obtained using conventional conductive lacquer barrier measurements.
2. Photogrammetric Techniques for Crack Detection and Crack Width Measurement in High-Speed Camera Image Sequences
The image sequence processing chain, which is used for crack detection and crack measurement in the work presented here, bears the name “cascaded image analysis” in ref. [
17]. It has been described in detail by Liebold and Maas [
18,
19,
20] and is briefly summarized here. We limit ourselves to the analysis of specimens with planar surfaces as recorded by monocular image sequences.
As a first step, a grid of points is defined in the camera image of the first epoch under zero load. The grid points are tracked through the consecutive images of the sequence using a sub-pixel accuracy least-squares image matching technique [
21] in order to obtain displacement fields for each epoch with regard to the first epoch.
Figure 1 shows the displacement field from a bend test on a concrete beam with a thinned-out point set, wherein cracks cause discontinuities in the displacement vector field.
As an alternative to defining a regular grid of points to be tracked through the image sequence, an interest operator may be used to detect surface points with good image contrast, thus guaranteeing precise least-squares matching results. These points, which are not distributed in perfect regularity, are then triangulated into a mesh. The triangles of all images after the zero-load epoch are subsequently analyzed for deformations. Triangles with cracks penetrating them entirely show significant deformations, while the remaining triangles are unchanged. For that purpose, the relative translation vector
can be computed for each triangle based on the relative movement of its vertices [
19,
20].
Figure 2a schematically depicts a divided triangle used as the assumed model for the computation of the relative translation vector.
The norm of this vector
can be used as a deformation quantity for brittle material.
Figure 2b shows a color-coded map of this quantity. The norm is checked for exceeding a threshold value
, that is on an order of magnitude of the precision of the displacement field, i.e., approximately a tenth of a pixel. Due to the development of some distributional noise in the deformation field, filter methods such as the bilateral filter technique can be applied [
18]. In classifying crack regions, the connected components of triangles where
are determined by applying a region-growing technique including the hysteresis method with a second threshold, as performed by Canny [
22].
Figure 3 shows an example of the connected components of deformed triangles from a load test with a concrete beam.
Furthermore, it is also possible to derive crack widths and crack opening vectors for deformed crack triangles with the help of the relative translation vector
as shown by Liebold and Maas [
19,
20]. Crack widths can be computed as the projection of the relative translation vector
onto the crack normal
, which has to be determined first; see
Figure 4a. For each deformed triangle, the crack normal can be estimated by a line fitted through neighboring deformed triangles; see
Figure 4b.
The result of the image sequence processing chain described is a complete crack pattern for each image of the sequence plus subpixel accuracy information on local crack width. If a crack propagates, the crack pattern on the specimen’s surface and the measured crack widths will change. In the next step, these changes are analyzed to quantify the crack tip motion.
It is also possible to use and extend the algorithms presented here to measurements with stereo camera systems in order to analyze non-planar surfaces as shown by Liebold et al. [
23,
24].
3. Crack Tip Detection in Multi-Temporal Crack Pattern Images
To estimate the propagation velocity of cracks in an image sequence, it is necessary to detect the crack tips in each image. Direct detection of the crack tip in the images is restricted by the subpixel width of a crack near its tip. As a more reliable alternative, the crack tip position and motion can be determined indirectly using the algorithms presented in the previous section. Accordingly, the displacement field and the relative translation vectors are computed for each triangle of the mesh. Crack candidates are found by thresholding where
, and connected components are determined with a region-growing technique using the hysteresis method with a second threshold (δ
2 = 0.5 × δ
1); see
Figure 5.
In the following, the crack tip must be found in the connected component of the deformed triangles. Two possible methods are presented here, an approach based on triangle neighborhood analysis and a principal component analysis approach.
The former method considers the neighborhood of each crack triangle.
Figure 6 shows four possible cases in the neighborhood analysis of the crack candidates, shown as bold red triangles. In
Figure 6a–c, connected components of non-deformed triangles, i.e., darker blue triangles with edges marked in thicker, colored lines, in the set of the neighbor triangles are determined. The number of these components are counted.
Figure 6a shows the standard case with two components,
Figure 6b represents a crack junction with more than two connected components around the dark red crack candidate (sets of triangles with yellow or green edge marking). In
Figure 6c, only one component of non-deformed neighbor triangles is determined, implying a possible crack tip.
Figure 6d shows a crack candidate at the boundary of the mesh, which can be analyzed with the help of the halfedge data structure [
25]. A boundary triangle has at least one halfedge with no opposite halfedge partner.
Sometimes, the classification of crack tips can fail as shown in
Figure 7a. This might also be the case if the thickness of the area of deformed triangles is larger than one triangle, which can be caused by filtering methods, the smearing of
.
Figure 7b depicts the connected component of crack triangles (blue) in an experiment. Mesh boundary triangles are colored yellow and triangles that are identified as crack tips are depicted in green.
The latter method in crack tip identification is the principal component analysis (PCA), as shown in
Figure 8. The set of center points of the deformed triangles from the connected component are analyzed by computing the mean and the covariance matrix; see Equations (1) and (2) as well as
Figure 8a.
where
= center of triangle
i of the connected component,
= mean of the
values,
n = number of triangles of the connected component,
Σ = covariance matrix of the
values.
After this, the covariance matrix can be decomposed using an eigenvalue decomposition:
where
V = rotation matrix with the eigenvectors,
Λ = eigenvalue matrix (diagonal).
The center points of the triangles of the connected components can then be transformed:
The minimum and maximum values in the x and y directions belong to possible crack tip candidates. The dimension with the greater extension is chosen; see Equation (5) and
Figure 8b.
Figure 8c illustrates the extension of the transformed coordinates of the triangle centers as a rotated bounding box. If the crack tip candidate is a mesh boundary triangle, the candidate should be rejected.
Due to possible movements of the entire specimen between the images, the crack tip position of the previous image should be transformed into the current time step. The correction can be done by adding the mean vertex displacement changes to the previous crack tip triangle center.
The PCA method works well for cracks without branches, but may fail in cases of branching. During the experimentation, the authors could observe only one crack without branches. Here, the PCA method produced fewer outliers than the first presented method and consequently was used.
4. Crack Velocity Determination
Crack velocities can be obtained as derivatives of the crack lengths over time. This is often done by numerical differentiation with, for instance, finite differences, i.e., forward and central differences. As this may lead to high fluctuations due to uncertainties, the authors decided to compute regression lines. The least-squares method offers the advantage of providing standard deviations as internal error measures. The derivatives are computed at each data point, including the four nearest neighbors. The mathematic model for the regression is:
where
d = crack length,
t = time,
m = velocity,
c =
y–intercept.
This can also be written in matrix notation:
where
A = design matrix,
L = observation vector,
x = vector of unknowns,
v = residual vector.
For the i
th data point, the matrices are:
The measurements are weighted using binomial weights (1, 4, 6, 4, 1). The parameter vector
x can be obtained by minimizing the residuals and solving the normal equations:
where
P = weight matrix.
Then, the residual vector
v can be calculated:
The square of the
a posteriori standard deviation
of the unit weight can be computed as follows:
where
n = number of observations,
u = number of unknowns.
The covariance matrix
of the parameters can be computed as follows:
In addition to the actual velocity, the standard deviation of the velocity can be derived from Equation (15):
7. Experimental Results
For the detection of the deformed triangles, the algorithm described in
Section 3 and
Section 4 was applied.
Figure 11 shows color-coded maps of the relative translation vector lengths for each triangle for different time steps of one of the experiments. Then, the thresholding was done for the norm of the relative translation vector with different threshold values
(0.08, 0.10, 0.12 px). The connected components and the crack tips were determined by means of the PCA method.
Figure 12 shows an example with a threshold of 0.10 px.
The propagation length between the consecutive images, and accordingly its velocity, were determined automatically using the image processing chain described in the previous sections. The distance from the notch tip, that is the beginning of the crack, to the first crack tip was measured manually at the left border of the image since the notch tip was not included in the measuring field of the optical sensor. This measure is not necessary for the velocity determination but is useful for the comparison to the lacquer barriers in the crack length over the course of time.
Figure 13 shows the crack propagation lengths as a function of time for the different thresholds δ for
as well as for the first three conductive lacquer barriers for the four specimens tested. The last barrier was ignored due to crack branching which happened between the third and fourth barriers. In addition, the optical measurement range does not cover the last barrier; see also
Figure 10b. As expected, the photogrammetric method with the threshold of
= 0.08 px is the most sensitive, detecting crack tips earlier in the image sequence than if higher thresholds are used. The magenta-dotted line shows the crack length as a function of time for the conductive lacquer barriers. Comparing the methods, both techniques show monotonically increasing behavior of similar intensity. There is a slight mismatch in the time domain, which might be caused by synchronization errors but is not relevant to velocity determination.
Next, the velocities were computed as derivatives of the crack lengths with respect to time.
Figure 14 shows the results for the four specimens. The 5-point least-squares velocities obtained by the photogrammetric method for the different thresholds and the two velocities calculated between the conductive lacquer barriers are presented. In addition to the velocities, errors bars are plotted, showing the range of ± one standard deviation. The photogrammetric curves show some dependence on the threshold. This leads to an offset of the determined velocities on the time axis and to a certain fluctuation in the velocities as determined. In fact, as noted above, lower threshold values lead to an earlier detection of the crack tip. However, the resulting shift on the time axis is irrelevant to the velocity calculations. For the conductive lacquer barriers, only two velocity values can be derived. This is not sufficient for a comparison over the entire time span in which the photogrammetric method was used. However, the velocities are on the same order of magnitude. The standard deviations of the velocity values obtained with the photogrammetric method are in a range between 40 and 170 m/s, the standard deviations of the lacquer barrier based velocities reach several hundred m/s.
Due to the relatively high standard deviations, the image sequence based instantaneous velocities can also be translated into average velocities using the least-squares method for all data points.
Figure 15 shows the corresponding results: The average velocities obtained with the photogrammetric techniques vary between 650 and 900 m/s with standard deviations in the range of 30 to 90 m/s. The velocities between the conductive lacquer barriers are lower, in a range between 490 and 600 m/s, with higher standard deviations of up to 240 m/s. It is obvious that the velocities obtained from the high-speed camera data are consistently higher than the values obtained from the conductive lacquer barriers. However, it should be noted also that the standard deviations of the lacquer barrier measurements are much higher, and thus less reliable as only three data points are used in the computation. The number of conductive lacquer barriers in the measuring area of the camera could be increased in order to have more data points and to compare the methods, but this may also lead to increasing measurement errors due to the small distances between the barriers. This indicates that the quality of the conductive lacquer barriers measurements cannot be considered as a reference for the results obtained from the photogrammetric image processing chain due to their lower quality.
8. Conclusions
The article at hand presents a novel technique for determining crack propagation velocity on the basis of recorded high-speed camera image sequences. A precise, automatic, photogrammetric image sequence processing chain was developed, implemented, and validated. While using a high-speed camera means a relatively high instrumental effort, it provides the obvious advantage of offering spatio-temporally resolved, non-contact measurements rather than single-point measurements. The technique is able to detect cracks with a width on the order of 0.1 pixel and may thus be applied even in tracking fissures not visible to the human eye. The crack propagation velocities obtained from four laboratory experiments conducted in this study were on the order of 800 m/s.
For the experiments, a camera with a frame rate of 160,000 fps was used. The high frame rate comes with the drawback of an image format limited to 256 × 184 pixels. Future high-speed cameras will offer higher frame rates and larger image formats, which will be beneficial both for the number of cracks to be analyzed simultaneously and for the precision of the velocity values obtained by the method developed.