Next Article in Journal
Novel Antibacterial Materials and Coatings—A Perspective by the Editors
Next Article in Special Issue
Experimental Study of Drilling Damage Outcomes in Hybrid Composites with Waste Micro-Inclusions
Previous Article in Journal
Emergent Magnonic Materials: Challenges and Opportunities
Previous Article in Special Issue
Micro-Scale Numerical Simulation of Fatigue Failure for CFRP Subjected to Multiple-Amplitude Cyclic Loadings Based on Entropy Damage Criterion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Damage Detection in a Polymer Matrix Composite from 4D Displacement Field Measurements

1
Faculty of Mechanical Engineering and Naval Architecture, University of Zagreb, 10000 Zagreb, Croatia
2
Université Paris-Saclay, CentraleSupélec, ENS Paris-Saclay, CNRS, LMPS–Laboratoire de Mécanique Paris-Saclay, 91190 Gif-sur-Yvette, France
3
Institut für Angewandte Mathematik (IfAM), Leibniz Universität Hannover, 30167 Hannover, Germany
4
GE HealthCare, 78530 Buc, France
*
Authors to whom correspondence should be addressed.
Materials 2023, 16(18), 6300; https://doi.org/10.3390/ma16186300
Submission received: 7 August 2023 / Revised: 13 September 2023 / Accepted: 15 September 2023 / Published: 20 September 2023

Abstract

:
Standard Digital Volume Correlation (DVC) approaches enable quantitative analyses of specimen deformation to be performed by measuring displacement fields between discrete states. Such frameworks are thus limited by the number of scans (due to acquisition duration). Considering only one projection per loading step, Projection-based Digital Volume Correlation (P-DVC) allows 4D (i.e., space and time) full-field measurements to be carried out over entire loading histories. The sought displacement field is decomposed over a basis of separated variables, namely, temporal and spatial modes. In the present work, the spatial modes are constructed via scan-wise DVC, and only the temporal amplitudes are sought via P-DVC. The proposed method is applied to a glass fiber mat reinforced polymer specimen containing a machined notch, subjected to in situ cyclic tension and imaged via X-ray Computed Tomography. The P-DVC enhanced DVC method employed herein enables for the quantification of damage growth over the entire loading history up to failure.

1. Introduction

As it provides 3D images of scanned microstructures in a non-destructive way, X-ray Computed Tomography (XCT) has given a major impetus to the field of mechanics of materials [1,2]. With the recent developments of testing machines, 3D full-field displacement fields can be quantified in situ by coupling XCT and Digital Volume Correlation (DVC) into unique frameworks [3,4]. As DVC is based upon matching the gray levels between fully reconstructed volumes, the major limitations of such approach are its low temporal sampling and experiment duration. Each reconstructed volume requires from few minutes up to few hours of radiograph series to be acquired. This process governs the number of possible acquired scans. As the material response between two consecutive scans is not accessible, the quantification of time-dependent phenomena (e.g., stress-relaxation or crack propagation) is restricted. Furthermore, due to relaxation during in situ experiments, the reconstructed volumes may be impacted by motion artifacts. This issue can be avoided by starting the radiograph acquisition only after the measured force has stabilized. The aforementioned restrictions may be overcome by two different routes. The first route is to perform fast acquisitions, say in synchrotron facilities, with extremely bright X-ray beams using high-speed cameras. However, synchrotron imaging presents some drawbacks such as poor accessibility, high doses and the reconstruction quality that may degrade due to blur induced by vibrations [5,6]. Another route is to utilize the recently developed Projection-based Digital Volume Correlation (P-DVC), which aims to surpass the low temporal resolution attributed to classical (i.e., scan-wise) DVC approaches [7]. This method enables for the measurement of 4D (i.e., space and time) displacement fields from series of 2D radiographs acquired at different angles and load levels (instead of working with series of reconstructed 3D volumes). In such a way, huge gains in acquisition time and testing duration can be reached [4,8,9].
The ever-increasing application of fiber-reinforced polymers (FRPs) across various industrial branches is attributed to their advantageous properties such as high stiffness-to-weight ratios, which outperform conventional engineering materials [10]. Due to the occurrence of various damage mechanisms at different scales under mechanical loading, XCT is appealing in damage analyses of FRPs. The latter has been employed to quantify damage growth in FRPs under various loading conditions by, e.g., analyzing volumes [11,12,13,14]. However, it is worth noting that with this approach, only microstructural changes were evaluated, while it cannot provide access to the bulk kinematics, nor can it correlate it with microstructural changes. This limitation can be surpassed by DVC, which enables for 3D full-field displacement measurements [15,16]. Local approaches to DVC, which do not assume continuity of the displacement fields as they split the investigated Region of Interest (ROI) into smaller zones and register them independently, have been employed in the analysis of FRPs [17,18,19]. Conversely, global approaches based on, e.g., finite element (FE) discretizations were introduced more recently [20,21]. Such approaches are based on the assumption of continuity of displacement fields. A significant advantage of FE-DVC is the availability of correlation residuals. They correspond to the gray level differences between the reference volume and that of the deformed material corrected by the measured displacement field [4]. The conservation of gray levels is the underlying hypothesis of DVC [16,22,23]. Therefore, the correlation residuals can be used to check the quality of the registration. The initiation and growth of damage violates gray level conservation, thus the analysis of correlation residual fields (for converged displacement fields) reveal discontinuities corresponding to damaged zones [4,24,25]. Furthermore, FRPs are likely to undergo stress relaxation during tests [26,27]. As reconstructed volumes may last for up to hours, the specimen should be in a steady state to avoid poor image quality due to motion artifacts. Therefore, performing full-field displacement measurements in space and time is crucial especially for FRPs. The full insight into the material behavior is then provided in situ, and it is fully coupled with changes in the underlying microstructure.
In general, P-DVC relies on spacetime discretizations of the sought displacement fields [9]. A new spacetime framework was recently proposed to measure 4D displacements by combining DVC and P-DVC analyses [28]. By using spatiotemporal separations of variables, the spatial components were obtained via scan-wise DVC, and the instantaneous temporal amplitude was computed via P-DVC on each individual radiograph acquired on-the-fly during the whole loading history. In the following study, a new P-DVC enhanced DVC algorithm is introduced in which spatial modes are measured via DVC, and the temporal modes are sought with P-DVC. For the first time, the temporal modes are constructed such that they are compatible with the loading history of the experiment. In this regard, temporal shape functions are introduced. The displacement fields are sought in a 4D (i.e., 3D in space and 1D in time) vector space generated by a reduced spacetime kinematic basis. Only one projection per loading level is still needed. Additionally, different temporal interpolations are investigated, and the number of spatial modes selected to describe the specimen kinematics is analyzed and justified.
The aim of this work is the 4D characterization of a mat glass fiber reinforced polyester resin subjected to in situ cyclic tension, and imaged via XCT. A dogbone specimen contained a machined rectangular notch to induce high strain gradients. The properties and fabrication procedure of the investigated composite material are presented in Section 2, followed by a description of experimental setup. In Section 3, the principle of the proposed procedure is presented by introducing DVC and P-DVC formulations. The method is then applied to an in situ cyclic tensile test of a single-notched specimen (Section 4). The reconstruction error was first quantified by evaluating the difference between the corresponding projections and re-projections of the reconstructed reference volume. Last, the methodology to quantify damage growth from projection residual fields is presented. The implemented methodology is employed to directly assess damage from radiographs, thereby representing a completely novel approach.

2. Material and Methods

2.1. Material and Experimental Setup

In the present work, a polyester resin reinforced with a continuous glass fiber mat is investigated. The composite plate comprised 12 layers of R-glass fiber mat and was produced by manual lay-up, followed by compression molding. In terms of volume fraction, the composition had 40% of fibers, 55% of resin and 5% of voids (due to the fact that vacuum was not applied during molding). Figure 1 shows an optical micrograph of the investigated material. Due to the contrast between constituents, this material is well-suited for DVC, which relies on natural contrast [4].
The thickness of the investigated specimen was 5.2 mm. The specimen contained a machined rectangular notch (Figure 2a), which induced high strain gradients. The depth of the notch was 1.6 mm, while its width was 0.6 mm. The in situ cyclic tensile test was performed with the TTC Deben testing machine, i.e., the specimen was continuously loaded, rotated and imaged in the X50+ scanner (North Star Imaging) of LMPS. Before prescribing any loading, a reference scan (denoted with the black dots in Figure 2b), comprising 800 projections was acquired at equally spaced angles spanning over a complete 360° revolution with high-quality (HQ) scanning parameters (Table A1 in Appendix A). As this scan took approximately 2 h to be completed, the number of averaging frames was subsequently reduced from 20 to 1 to lower the experiment duration and mitigate, e.g., stress-relaxation. Such scanning parameters are called continuous scanning, and each scan took 4 min to be completed. After the acquisition of the reference scan, cyclic tension was applied at a constant stroke velocity of 4 µm/s. The radiographs were continuously acquired during the entire loading history at 4 frames per second with angular increments of 0.47º. These settings led to 768 radiographs per turn. It is worth highlighting that the P-DVC methodology presented herein is not dependent on the rate of the conducted experiments. However, there may be limitations related to the temporal resolution of each imaging device.
The deformed scans were reconstructed from the stages where a constant stroke was applied (marked with red dots in Figure 2b). They were used to measure 3D displacement fields via FE-DVC within the Correli 3.0 framework [29]. The FE mesh used in all DVC analyses (Figure 2a) was composed of first-order tetrahedral (T4) elements with piecewise linear (P1) shape functions. It was made of 88 nodes (with three DOFs per node), and the mean element length was equal to 12 vx.
A total of 17,442 projections was acquired during the prescribed loading history. The radiographs had an initial definition of 1944 × 1536 px. The focus was put on the notched-region (as the majority of damaged zones concentrated within this region [27]). Coarse graining of 4 × 4 elementary pixels into one superpixel was performed. This led to 486 × 384  px radiographs (called images at scale 4). The physical length of one voxel was equal to 58 µm. The volumes were reconstructed within the ASTRA toolbox with the simultaneous iterative reconstruction (SIRT) algorithm suited for cone beams [30], and employing geometrical parameters given by the tomograph calibration.

2.2. Projection-Based Digital Volume Correlation

This section outlines the basic principles of DVC and P-DVC. The notations used herein are introduced as well. The radiographs (i.e., the sum of absorption coefficient along each material point x of the beam ray hitting the detector at position r) are denoted as p ( r , t ) at different time t. The reconstructed volume is such that its projection should match the recorded projection p ( r , t ) , which results in a linear relationship between f and p
Π θ ( t ) [ f ( x ) ] = p ( r , t ) ,
where Π θ ( t ) is the projection operator at angle θ ( t ) . The reconstruction is thus the inversion of this linear system for a large sequence of angles [31]. By reconstructing a series of such tomographic images, the 3D kinematics of a medium can be quantified. The displacement field is estimated via DVC, which relies on the conservation of gray levels
f T 0 ( x ^ ) = f T ( x ^ + U ^ ( x ^ , T ) ) ,
where x ^ denotes the position of any material point in the reference configuration, and U ^ the measured displacement field that describes the change of the reference volume f T 0 (usually captured in the undeformed state) to the deformed configuration (here denoted as f T ). Two different time scales have to be highlighted here. First, the 3D images are reconstructed over the duration needed for one turn T. Hundreds up to thousands of instantaneous radiographs acquired at θ ( t ) are needed over a full revolution (or turn T) of the specimen. The turn parameter T thus denotes the different stages at which a 3D scan was captured, while T 0 denotes that of the reference scan. The displacement field to be measured is obtained from the minimization of Γ DVC , which is the quadratic differences between the corrected deformed and the reference volumes
Γ DVC = x ^ ( f T ( x ^ + U ^ DVC ( x ^ , T ) ) f T 0 ( x ^ ) ) 2
with respect to the parameterization of the trial displacement fields U ^ DVC . A global (i.e., FE-based) approach was utilized in this work. The displacement field is then expressed in a kinematic basis consisting of FE shape functions Ψ ^ j
U ^ DVC ( x ^ , T ) = j v j ( T ) Ψ ^ j ( x ^ ) ,
where v j are the nodal displacements. As the sought displacement fields are determined with respect to the reference configuration, the DVC framework presented herein is Lagrangian.
Instead of working with fully reconstructed volumes, P-DVC measures 3D displacement fields from a series of 2D radiographs captured at different angles θ ( t ) and loading steps F ( t ) [32]. The cost function Γ P - DVC to be minimized is defined as follows
Γ P - DVC = r , t ( Π θ ( t ) [ f T 0 ( x u ( x , t ) ) ] p ( r , t ) ) 2 .
According to Equation (5), the problem is no longer stated in the Lagrangian framework, but in the Eulerian system (contrary to Equation (4)). To complete the notations needed to describe the motion in Lagrangian and Eulerian settings, it is worth highlighting that the displacement at any spatial location x reads
u ( x , t ) = x x ^ = u ^ ( x ^ , t ) .
The present P-DVC framework requires the acquisition of one reference volume f T 0 from which the microstructure of the sample is known. The remaining unknowns are the displacements sought in the spatiotemporal framework. In the present study, the spatial components are constructed with all scan-wise DVC results, while only temporal modes are sought via P-DVC
u ( x , t ) = τ α τ U ^ DVC ( x ^ , τ ) σ τ ( t ) ,
where σ τ ( t ) denote the temporal modes, and α τ the temporal amplitudes. The presented methodology is 4D in the sense that it gives access to the entire displacement field u ( x , t ) in (3D) space and (1D) time. Furthermore, DVC fields are the result of the minimization (Equation (3)) when the displacement fields are parameterized with FE shape functions (Equation (4)).
The displacement fields are sought in a vector space generated by a reduced kinematic basis. The interpolation functions σ τ naturally introduce temporal regularization [9]. The nodal displacements v j ( τ ) are measured via FE-DVC, and the temporal amplitudes gathered in column vector { α } are sought via P-DVC and iteratively updated
{ α l + 1 } = { α l } { δ α l } ,
where the corrections are calculated by performing Gauss-Newton minimizations of the functional Γ P - DVC
{ δ α l } = [ M l ] 1 { m l } .
In Equation (9), the Hessian matrix M τ k is updated for each iteration l
M τ k l = r , t ( S τ l ( r , t ) σ τ ( t ) ) · ( S k l ( r , t ) σ k ( t ) ) ,
where S τ l is the projected sensitivity at iteration l
S τ l ( r , t ) = Π θ ( t ) [ U ^ DVC ( x ^ , τ ) · f T 0 ( x u l ( x , t ) ) ] .
Furthermore, the second member vector m τ l is based on the projection residual fields
m τ l = r , t S τ l ( r , t ) σ τ ( t ) φ C l ( r , t ) ,
where the projection residual field per angle at each iteration l reads
φ C l ( r , t ) = p ( r , t ) Π θ ( t ) [ f T 0 ( x u l ( x , t ) ) ] .
According to Equation (13), the reference volume, corrected by the displacement field u l ( x , t ) , and projected for each angle θ ( t ) should match the acquired projection p ( r , t ) . The projection residual field for all angles is to be minimized during the iterative procedure. As in DVC, these residuals reveal what was not captured by the 4D corrections, herein projected onto the 2D detector plane. A general overview of 4D analyses in which spatial modes are controlled by DVC fields is shown in Algorithm 1.
Algorithm 1: P-DVC with DVC spatial modes
 Select spatial modes U ^ DVC ( x ^ , τ )
 Choose an initial guess of { α 0 }
while { δ α l } > 10 3 do (
  Corrections f l ( x , t ) f T 0 ( x u l ( x , t ) )
  Update [ M l ] and { m l }
  Solve { δ α l } = [ M l ] 1 { m l }
  Update temporal amplitudes { α l + 1 } = { α l } { δ α l }
  Update displacement fields u l + 1 ( x , t ) = τ α τ l + 1 U ^ DVC ( x ^ , τ ) σ τ ( t )
  Update projection residuals φ C l + 1 ( r , t ) = p ( r , t ) Π θ ( t ) [ f T 0 ( x u l + 1 ( x , t ) ) ]
end (

3. Results

In this section, the previous P-DVC enhanced DVC framework is applied to an in situ cyclic tensile test on a single-notched specimen. First, the reconstruction error of tomographic volumes is evaluated. The results of 4D measurements are then presented. All presented analyses converged, i.e., the norm of the change of amplitude corrections { δ α } (Equation (9)) between two iterations became less than 10 3 . Convergence and the obtained temporal amplitudes are discussed.

3.1. Evaluation of Reconstruction Error

According to Equation (13), the P-DVC algorithm minimizes the projection residual fields (i.e., differences between the acquired radiograph p ( r , t ) and the reference volume deformed by the measured displacement field, and projected according to the angle θ ( t ) ). This procedure is herein extended to the initial step (i.e., the reference HQ scan 0), where u = 0 . The reconstruction error was quantified by evaluating the difference between the corresponding projections and re-projections of the reconstructed reference volume 0. Due to acquisition noise, uncertainties in the geometric parameters used for the reconstruction and assumptions in the projection operator (e.g., pixel/voxel integration, beam hardening), the initial projection residual field was not equal to 0.
The root mean square (rms) of projection residuals ( φ C ), expressed in arbitrary units, evaluated for the first 360° revolution (scan 0) is shown in Figure 3a. The projection residuals within the whole ROI are low with respect to the mean dynamic range of the original radiographs (i.e., 0.7 a.u.). The highest residual levels are reached for θ = 90 and θ = 270 . Figure 3b,c shows the projection residuals for θ = 90 and θ = 270 . These initial projection residual maps are shown with a divergent color map to highlight positive and negative values. The highest rms( φ C ) levels (i.e., θ = 90 and θ = 270 ) are due to the specimen orientation with respect to the X-ray source. The thickness of the specimen with respect to the X-ray source is the highest at those angles, which may introduce nonlinearity in the X-ray attenuation. The initial projection residual field for θ = 0 reveals elevated levels close to the edges of the notch. The present volumes were reconstructed with the SIRT algorithm, which iterates between forward and backprojections until it reaches convergence. The algorithm was initialized with the Feldkamp–Davis–Kress (FDK) procedure. SIRT is computationally more expensive compared to FDK. However, it provides images with reduced noise. In the present case, the volumes were reconstructed with 100 SIRT iterations. The reconstruction residuals may be reduced by further adjusting the geometric parameters of the tomographic setup, and increasing the number of SIRT iterations. Other sources for such errors could be small motions during scanning and beam hardening, which were not corrected in the present work.

3.2. Full-Field Measurements over the Entire Loading History

All the results presented in the sequel are obtained with a temporal basis consisting of 11 plateau-like functions interpolating the cyclic loading history (Figure 4). The selection of this temporal basis is detailed in Appendix B. The temporal shape functions are linked to 11 spatial modes. The latter ones are constructed via DVC, i.e., by performing volumetric correlations between scans 0 (reference scan) and all subsequent reconstructed volumes (Figure 2b). It is important to emphasize that all P-DVC measurements are performed until specimen failure (i.e., after acquiring scan 009). The last temporal mode is related to the DVC field of scan 007. The (small) force drop between scans 007 and 008 was due to the fact that the continuous acquisition between these two steps was interrupted by acquiring one high quality scan, which was excluded from the present analyses.
The projections to be registered with P-DVC are selected so that the angle between two consecutive acquisitions is approximately equal to 4.6°. This choice results in 1491 projections to be analyzed (out of a total of 17,442 acquired during the complete loading history). The chosen temporal sampling leads to 77 projections per turn (360° revolution). The sampling is governed by the computation time. In P-DVC procedures, the heaviest operations are the computation of the Hessian matrices and the deformation of the reference volume at each time step. To perform measurements in reasonable time, such trade-off had to be made between temporal resolution and computation time. Compared to standard DVC analyses that would yield 11 fields in the present analyses, the full P-DVC analysis probes 1491 radiographs corresponding to different temporal states (i.e., more than 2 orders of magnitude higher than regular DVC).
The initial values in the amplitude vector { α 0 } are set to 1 (Equation (7)). The P-DVC analysis takes 4 iterations to converge. Each iteration lasts approximately 5 h. The measured displacements yield a global decrease of the residuals within the investigated ROI as shown in Figure 5.
The converged temporal amplitudes are gathered in Table 1. They are all close to 1, except for the last one where the amplitude is greater than 1. This is due to the fact that the DVC displacement field of scan 007 (Figure 2b) is used as a spatial component of the last mode, while 4D kinematic measurements are performed up to specimen failure.
Furthermore, the rms differences between the 3D nodal displacements obtained with FE-based DVC and projection-based measurements are calculated. Since large RBRs in x and y directions are observed during the experiment, these differences are only calculated for z-displacements as the mechanical component is dominant in that direction. The rms differences are computed between the DVC displacement fields and converged P-DVC displacement fields for the time-steps corresponding to the beginning of each constant stroke stage. The corresponding values are gathered in Table 2. As the projections acquired at the maximum load levels are affected by time-dependent motions (i.e., stress-relaxation phenomena), which influence the reconstruction of the full volume, differences are expected. From the values reported in Table 2, it is concluded that, when employing 11 plateau-like temporal functions, the rms differences are very small. With such settings, the kinematics of the specimen was properly captured.

4. Discussion: Quantification of Damage Growth

In this section, the methodology to quantify damage growth from projection residual fields is outlined and discussed for one particular case. The validation of the proposed procedure is detailed in Appendix C. The starting point of the method is to exploit the fact that DVC residuals revealed damaged zones within the specimen. Figure 6a shows the DVC residual field of scan 007 (last loading plateau prior to specimen failure). This map reveals numerous damaged zones in the specimen bulk, the majority of which are concentrated within the notched region. Next, the absolute maximum value of correlation residuals within the whole ROI is determined. The following criterion is employed to find finite elements that have correlation residuals greater than 70% of the maximum value integrated over each element volume
rms ( φ C , DVC ) Mask > 0.7 · rms ( φ C , DVC ) ROI , max
Based on the previous criterion, 14 elements are found. These elements compose the mask shown in Figure 6b, which is used to inspect the projection residuals. The mask is projected onto the detector plane with the projection operator Π θ ( t ) . The quantification of damage growth is performed within each element of the mask.
The analysis of damage growth is performed within element #89, which is located below the notch root (marked with red contour in Figure 6a. This element contains the zone of elevated correlation residuals, which correspond to a crack propagating from the notch root (Figure 6c).
As 1491 radiographs are analyzed with P-DVC, not all angles are inspected in the damage procedure. To find for which angles damage (i.e., mesocrack) is most pronounced, the DVC residuals within element #89 are projected onto the detector plane. First, the rms residuals are computed for the full specimen revolution for both loaded and unloaded plateaus. The residuals of the loaded scans are shown in Figure 7a, whereas those of the unloaded acquisitions in Figure 7b. The rms residual levels of scan 000 (acquired in the undeformed state) is included in order to track the increases in DVC residuals with respect to the base level. For scan 001, there is an increase in correlation residuals with respect to levels of scan 000. When the specimen is unloaded (scan 002), the residuals remain close to the levels of scan 001. This is trend is also visible when plotting the rms DVC residuals of element #89 during the entire returns for the acquisitions of scans 001 and 002. These increases in correlation residuals compared to the base level are attributed to the specimen kinematics [27] and damage still did not occur at this stage.
As the specimen is once again loaded (scan 003), the correlation residuals keep increasing. Upon subsequent unloading, the residual levels decrease and are very close to those of scan 002. Damage growth did not yet occur, or if it had, it remained very low and at a scale lower than the scan resolution [33]. The significant increase in the correlation residuals compared to previous acquisitions is visible for scan 005. This is especially pronounced in Figure 7c. At this stage, the crack originating from the notch root and propagating below it is visible on DVC correlation residual field [34]. The subsequent unloading (scan 006) led to partial crack closure (as the correlation residuals were higher than those of the previous unloaded acquisitions, Figure 7c). At the peak of the last loading plateau prior to specimen failure (scans 007 and 008), the correlation residuals once again increase, and upon specimen unloading (scan 009), partial crack closure appears. From the reported data, the first interval of residual increase is within the range of 50° to 110°. In this range, damage is expected to be pronounced on residual fields. Thus, the projected DVC residuals within element #89 are shown in Figure 7d for an angle of 75 . The area containing elevated residuals corresponding to the crack is marked with the red arrow.
With the previous observations, the projection residuals are only inspected for specific angles, namely, 0 (to prove that this angle is damage-insensitive for this specific case), 56°, 75° and 105° acquired during the loading history. First, damage growth is analyzed by tracking the increase in rms ( φ C P DVC ) with respect to the base level. The rms projection residual levels are divided by that of scan 000 (prior to loading). This is the so-called amplification factor shown in Figure 8.
For θ = 0 , there is no increase that may reveal damage growth, except for the very last loading step prior to specimen failure, which may also be attributed to the extrapolation using the DVC displacement field of scan 007. The same conclusions are drawn for θ = 56 . Furthermore, the amplification levels when θ = 75 reveal the first significant increase in projection residuals compared to the base level for time step 941, which corresponds to the range of scan 007 (Figure 2b). This observation is consistent with DVC results at a finer scale [34].
The accumulation of damage and the activation and deactivation of cracks (due to the cyclic loading history) is most visible for θ = 105 . It is concluded that damage within element #89 first initiated at time step 331, which corresponds to the range of scan 003 (Figure 2b). During the subsequent unloading (time step 485), it appears that damage is deactivated. Damage within element #89 is once again active upon loading the specimen to the peak of the third loading cycle (time step 639). However, when the specimen is once again unloaded (time step 793), the crack within element #89 is only partially closed (the amplification of time step 793 is greater than that of time step 485). Within the range of scans 007 and 008 (time steps 947 to 1255) the crack is again active and more pronounced compared to the previous loading stages. The subsequent unloading (time step 1332), once again leads to partial crack closure.
A singular value decomposition (SVD) is applied to the projection residuals constructing the matrix [ U ] for each analyzed angle ( θ = 0 °, 56°, 75° and 105°) as explained in Appendix B. The singular values are shown in Figure 9. For all four analyzed angles, the first singular value is significantly higher than the second one, thereby indicating a good separation between them. Thus, only the first temporal mode of each angle is analyzed for the purpose of quantification of damage growth.
Figure 10 shows the first temporal modes calculated when θ = 0 °, 56°, 75° and 105°. For θ = 0 ° and 56°, the values for different loading steps (in this case marked with turns) are distributed within the same range. Thus, for these angles, damage is not detected. The increases in values for the last turn prior to specimen failure mainly stem from the extrapolation of the DVC field of scan 007. Furthermore, the values obtained for θ = 75 ° reveal the first increase for the loading stage corresponding to scan 005 (turns 9 and 10). However, for subsequent unloading the crack is only partially closed. This remark is consistent with DVC results on a finer scale [34]. Damage growth is pronounced for scans 007 and 008 (turns 13 and 14). The subsequent unloading resulted in partial crack closure. The values reported for θ = 105 ° clearly highlight the accumulation of damage within element #89. These values also reveal the activation and deactivation of damage. Thus, this angle is deemed the most sensitive to the presence of cracks.
The first spatial modes of the projection residuals computed for the selected angles are shown in Figure 11. For θ = 0 °, 56° and 105°, the spatial distribution is homogeneous. The area of elevated values is visible near the element angle when θ = 75 °. The reported results demonstrate that the analysis of projection residuals enables for the quantification of damage growth within the investigated specimen. However, due to the heterogeneous and random mesostructure, the damage analysis is highly dependent on the analyzed angles.

5. Conclusions

In the present work, a 4D characterization of a glass fiber mat reinforced polyester resin was performed by employing DVC enhanced P-DVC that relied on spacetime discretization of the measured displacement fields. An alternative approach in which spatial modes were fully constructed with DVC, and only temporal modes were sought via P-DVC was implemented. All the analyses were performed with a temporal basis consisting of 11 plateau-like functions interpolating the cyclic loading history. These temporal shape functions were linked to 11 spatial modes constructed via DVC (i.e., by performing volumetric correlations between scan 0 and all subsequent reconstructed volumes).
A methodology to quantify damage growth was proposed by analyzing projection residuals for specific angles. First, zones in which damage growth was to be quantified were selected based on the DVC gray level residuals. Finite elements exhibiting correlation residuals greater than 70% of the maximum value within the entire ROI composed the mask. This mask was then projected onto the detector plane. In the present study, the analysis of damage growth was performed for one element that contained elevated correlation residuals corresponding to crack propagating from the notch root. The quantification of damage growth was performed by multiplying the projection residuals of this element with the mask of this element. As the temporal sampling was such that 1491 radiographs were analyzed with P-DVC, not all angles were inspected. To find for which angles damage was most pronounced, the DVC gray level residuals of the same element were projected onto the detector plane. Based upon the root mean square residual, specific angles were analyzed. Damage growth was first analyzed by tracking the increase in terms of rms values with respect to the base level. An amplification factor that revealed the increase in residual levels during the cyclic loading history indicated the activation and deactivation of damage. Last, an SVD procedure was also applied to the projection residuals. It was shown that the analysis of projection residuals may enable for the quantification of damage growth within the investigated specimen.
One of the principal aspects of future investigations will be to modify the temporal sampling in a way to find angles most sensitive to cracks. In such a way, the angles that do not display damage would be excluded from the analysis. Further, the presented results were obtained using coarse grained images (i.e., the resulting definition was down-sampled over a 4 × 4 regular grid). This binning level may lead to a loss of information beneficial for mechanical characterization. Another aspect will thus be to perform 4D measurements at the original scale. In addition, the measurements may be improved by further adjusting the geometrical parameters of the tomography setup, as well as performing beam hardening corrections (which was not carried out in the current work).
The P-DVC procedure employed in the present work offers the opportunity to continuously measure 4D (i.e., 3D space and 1D time) displacement fields with simultaneous loading and rotation of specimens. In the proposed framework, the spatial modes were constructed with scan-wise DVC analyses, and only temporal modes were sought via P-DVC. The P-DVC enhanced DVC method enabled for the quantification of damage growth over the entire loading history up to failure. The continuous acquisition of radiographs over the entire loading history and the analysis of only one radiograph per loading step significantly enriched the temporal sampling compared to scan-wise DVC (i.e., by more than two orders of magnitude in the present case). Due to its rich temporal sampling, the proposed framework may also enable for the assessment of time-dependent phenomena such as, e.g., stress-relaxation and crack propagation. It is important to highlight that the employed P-DVC methodology is generic in a way that it can be applied to different materials with various damage mechanisms and behavior, as well as different in situ experimental protocols (when specimens are simultaneously loaded and imaged). In addition to enabling experimental tests to be performed very quickly, it also allows different in situ experimental protocols (designed for standard full-field measurements) to be optimized. This procedure may significantly reduce the quantity of data required for volumetric full-field measurements, while enriching the temporal sampling compared to standard DVC approaches.

Author Contributions

Conceptualization, A.M., F.H. and Z.T.; methodology, C.J., B.S. and F.H.; software, A.M., C.J. and V.K.; validation, V.K., C.J. and F.H.; formal analysis, A.M., V.K. and F.H.; data curation, B.S.; writing—original draft preparation, A.M.; writing—review and editing, V.K., C.J., Z.T. and F.H.; supervision, F.H. and Z.T.; project administration, Z.T.; funding acquisition, Z.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was performed within the FULLINSPECT project supported by the Croatian Science Foundation (UIP-2019-04-5460 Grant). This work was also partially supported by the French “Agence Nationale de la Recherche” through the “Investissements d’avenir” program (ANR10-EQPX-37 MATMECA Grant). The financial support of the French-German University through the French-German Doctoral college “Sophisticated Numerical and Testing Approaches” (CDFA-DFDK 19-04) is acknowledged.

Data Availability Statement

The dataset is currently unavailable.

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.

Abbreviations

The following abbreviations are used in this manuscript:
DVCDigital Volume Correlation
FEFinite Element
FDKFeldkamp–Davis–Kress
FRPFiber Reinforced Polymer
HQHigh Quality
P-DVCProjection-based Digital Volume Correlation
RBMRigid Body Motion
RBRRigid Body Rotation
rmsRoot mean square
ROIRegion Of Interest
SIRTSimultaneous Iterative Reconstruction
SVDSingular Value Decomposition
TTCTension, Torsion and Compression
XCTX-ray Computed Tomography
4DFour-dimensional (3D space and 1D time)

Appendix A. Hardware Parameters

Table A1. Scanning parameters.
Table A1. Scanning parameters.
TomographNorth Star Imaging X50+
X-ray sourceXRayWorX XWT-240-CT
Target/AnodeW (reflection mode)
FilterNone
Voltage150 kV
Current7.8 µA
Focal spot size5 µm
Tube to detector910 mm
Tube to object53 mm
DetectorDexela 2923
Definition1944 × 1536 pixels (2 × 2 binning)
Scanning settingsHigh qualityContinuous
Number of projections800768
Angular amplitude360°360°
Frame average20 per projectionContinuous (1 per step)
Frame rate3 fps3 fps
Acquisition duration1 h 46 min 26 s4 min 26 s
Reconstruction algorithmSIRTSIRT
Field of view7.3 × 10 × 20.5 mm37.3 × 10 × 20.5 mm3
Image scale14.6 µm/voxel14.6 µm/voxel

Appendix B. Selection of Temporal Bases

In this appendix, the selection of a suited temporal basis is discussed. This analysis is first conducted by performing a singular value decomposition (SVD) on the DVC displacement fields. An SVD is applied to rectangular matrices [ U ] m × n , with m > n in this work. To construct this matrix, the nodal displacements of each DVC result are gathered in a column vector, and then concatenated to construct a rectangular matrix [ U ] that corresponds to a space/turn separation of the measured data.
The mesh consists of 88 nodes (with 3 DOFs per node), and 10 DVC fields are used (namely, scans 000 to 009). Therefore, a 264 × 10 [ U ] matrix is created. Three matrices are generated after conducting SVD: the singular values matrix [ S ] m × n , and the modal matrices [ V ] m × m and [ T ] n × n . These matrices allow for the reconstruction of matrix [ U ]
[ U ] = [ V ] [ S ] [ T ] ,
where [ S ] is a rectangular matrix, and the non-zero terms of this matrix are the singular values s i of the matrix ordered by importance ( s 1 s 2 s 3 s n ). Due to their ranking, this matrix indicates that certain modes are more significant than others by several orders of magnitude. The spatial modes (SM) are gathered in matrix [ V ] , while the temporal modes (TM) are gathered in [ T ] . These modes can be truncated to obtain approximations [ U k ] with the first k instead of all n modes
[ U k ] = j = 1 k s j { V j } { T j } ,
where s j is the j-th singular value { V j } and { T j } is the corresponding modal vectors.
Furthermore, the displacement fields measured with DVC consisted of rigid body motions (RBMs) U rbm DVC ( x ^ , T ) and the part causing the deformation called mechanical displacement U mec DVC ( x ^ , T )
U ^ tot DVC ( x ^ , T ) = U ^ rbm DVC ( x ^ , T ) + U ^ mec DVC ( x ^ , T ) .
The SVD procedure is applied to both U ^ tot DVC (i.e., the displacement fields containing RBMs) and U ^ mec DVC (i.e., the displacement fields from which RBMs are subtracted). This analysis is performed by using the displacement fields of scans 000 to 009. Figure A1a shows the spectrum of singular values normalized by the highest one (i.e., s 1 ). From the reported data, it is shown that for both analyzed cases s 1 is greater than s 2 by about one order of magnitude, which indicates that Mode 1 captures most of the 4D kinematics. In the case of U ^ tot DVC , the difference between the first two singular values was smaller compared to U ^ mec DVC . This may be attributed to the influence of RBMs. It is worth highlighting that large RBMs are observed during the experiment, and hence they participate actively to the first mode of U ^ tot DVC .
The temporal modes obtained with U ^ tot DVC are shown in Figure A1b. From T tot 1 , a clear temporal signature resembling the cyclic loading history is observed (Figure 2b). Because of the orthogonality requirement of SVD, the fluctuations for higher order modes are more difficult to interpret.
Figure A1. (a) Spectrum of normalized singular values. (b) First three temporal modes.
Figure A1. (a) Spectrum of normalized singular values. (b) First three temporal modes.
Materials 16 06300 g0a1
The first three spatial modes of U ^ tot DVC are shown in Figure A2. In the first two spatial modes, large rigid body rotations (RBRs) are observed in the in-plane components, which occurred during the experiment. For the loading direction, the displacements of the first two modes correspond to tensile loading of the notched sample. The third mode is much more difficult to elucidate.
Figure A2. First (a), second (b), and third (c) spatial modes of U ^ tot DVC .
Figure A2. First (a), second (b), and third (c) spatial modes of U ^ tot DVC .
Materials 16 06300 g0a2
The first two modes are compared to the displacement field of scan 007 (Figure 2b), as shown in Figure A3. The first mode is easily interpreted as the dominant kinematic mode.
Figure A3. Displacement field of scan 007 containing RBMs.
Figure A3. Displacement field of scan 007 containing RBMs.
Materials 16 06300 g0a3
To confirm the previous statement, Table A2 gathers rms differences between the DVC displacement field of scan 007 and rescaled 1st, 2nd, and 3rd spatial modes. The difference for the first spatial modes is very small compared to those of the two higher modes.
Table A2. Rms difference expressed in vx between the DVC displacement field of scan 007 and rescaled 1st, 2nd, and 3rd spatial mode.
Table A2. Rms difference expressed in vx between the DVC displacement field of scan 007 and rescaled 1st, 2nd, and 3rd spatial mode.
xyz
1st mode0.340.080.09
2nd mode4.62.01.2
3rd mode4.42.77.8
Last, the truncation error evaluates the loss of information caused by using a subgroup of modes to reconstruct the kinematics. It is computed as
Truncation error = [ U ] [ U k ] 2 [ U ] 2 .
Figure A4 shows the truncation error as a function of mode number. If only one mode is kept, the truncation error is significantly larger compared to the error if more modes are kept. For the first mode, when the total displacement field is considered, the truncation error is significantly higher than that of the mechanical field. This difference is due to significant RBRs. Since the truncation error when only 1 mode is kept is not too high for the total displacement fields, it is decided to investigate two reduced bases consisting of one single mode. Both of them are compared to results obtained with the full set of DVC displacements.
Figure A4. Truncation error as a function of mode number.
Figure A4. Truncation error as a function of mode number.
Materials 16 06300 g0a4
Different temporal interpolations are exploited on the chosen temporal sampling. First, the 4D kinematics is measured using the normalized force signal as a unique temporal mode (Figure A5a). The measured uniaxial force history reveals rapid but limited relaxations during the acquisition of scans 001, 003, and 005 (the drop of the peak force is approximately 160 N, see Figure 2b). During the acquisition of scans 007 and 008, the force decrease is approximately 50 N and 20 N, respectively.
Next, the normalized stroke is used as a unique temporal mode (Figure A5a). In both cases, only one DVC field is used as spatial mode, namely, the total displacement field of scan 007 (corresponding to the last loading cycle prior to specimen failure). This field is obtained by running DVC between HQ scan 0 (reference) and scan 007. It is important to emphasize that 4D measurements are performed until specimen failure (i.e., after acquiring scan 007). This extrapolation is justified by the fact that scan 007 is acquired at the peak of the last loading cycle prior to specimen failure.
The third approach consists in using 11 plateau-like functions interpolating the cyclic loading history (Figure A5). These temporal shape functions are linked to 11 spatial modes. The latter ones are constructed via DVC, i.e., by performing volumetric correlations between scans 0 (reference scan) and all subsequent reconstructed volumes (Figure 2b). It is important to emphasize that all P-DVC measurements are performed until specimen failure (i.e., after acquiring scan 009). The last temporal mode is related to the DVC field of scan 007. The (small) force drop between scans 007 and 008 is due to the fact that the continuous acquisition between these two steps was interrupted by acquiring one high-quality scan, which is excluded from the present analyses.
Figure A5. (a) Temporal bases consisting of either normalized force history or stroke signal. (b) Temporal modes consisting of 11 plateau-like functions. The black dashed line represents the measured force signal normalized by its maximum value.
Figure A5. (a) Temporal bases consisting of either normalized force history or stroke signal. (b) Temporal modes consisting of 11 plateau-like functions. The black dashed line represents the measured force signal normalized by its maximum value.
Materials 16 06300 g0a5
The initial value of the amplitude vector { α 0 } is set to 1 (Equation (7)). In all three cases, the P-DVC analysis takes four iterations to converge. Each iteration lasts approximately 5 h. The measured displacements yield a global decrease of the residuals within the investigated ROI, as shown in Figure A6. The rms projection residuals within the ROI are the lowest when using 11 plateau-like functions. From the reported results, it is shown that utilizing only one DVC displacement field (as a spatial mode) is sufficient to properly capture the kinematics up to specimen failure.
Figure A6. Change in rms projection residuals over the entire Region of Interest during the iterative procedure for different temporal interpolations.
Figure A6. Change in rms projection residuals over the entire Region of Interest during the iterative procedure for different temporal interpolations.
Materials 16 06300 g0a6
Furthermore, the rms differences between the 3D nodal displacements obtained with FE-based DVC and projection-based measurements are calculated. Since large RBRs in x- and y-directions are observed during the experiment, these differences are only calculated for the z-direction displacements as the mechanical component is dominant in that direction. The rms differences are calculated between DVC displacement fields and converged P-DVC displacement fields for the time-steps corresponding to the beginning of each constant stroke stage. The obtained values are gathered in Table A3. As the projections acquired at the maximum load levels are affected by time-dependent motions (i.e., stress-relaxation phenomena), which influence the reconstruction of the full volume, differences are expected. From the values reported in Table A3, it is concluded that, when using measured force or prescribed stroke history as a temporal basis in combination with a single spatial mode, the rms differences are higher for the unloaded stages. This is especially pronounced when using the measured uniaxial force as the temporal interpolation basis, due to the fact that part of the deformation mechanisms is deactivated [34]. For the loaded stages, the differences are very small. When employing 11 plateau-like temporal functions, the rms differences are very small. The reported values confirm that, despite using only one spatial mode (i.e., the displacement field of scan 007) when performing 4D measurements, the kinematics is very well captured. In addition, employing the stroke history as the temporal interpolation better captured the temporal fluctuations (as the rms differences are smaller compared to those obtained when using the measured uniaxial force). However, the rms projection residuals within the ROI are the lowest when using 11 plateau-like functions. With such settings, the kinematics of the specimen is properly captured. This observation is also consistent with the truncation analysis where it was shown that the truncation error is significantly lower if more modes are kept (compared to the error if only one mode is used, see Figure A4).
Table A3. Comparison of rms differences in z direction between the 3D nodal displacements measured with FE-DVC and projection-based measurements. The comparison is performed for the converged P-DVC displacement fields for the angles at the beginning of each constant stroke stage.
Table A3. Comparison of rms differences in z direction between the 3D nodal displacements measured with FE-DVC and projection-based measurements. The comparison is performed for the converged P-DVC displacement fields for the angles at the beginning of each constant stroke stage.
Fu σ
ScanMechanical State z z z
000Unloaded0.050.050.005
001Loaded0.810.050.007
002Unloaded0.350.060.01
003Loaded0.880.060.007
004Unloaded1.560.220.012
005Loaded0.730.080.002
006Unloaded2.770.30.011
007Loaded0.780.090.012
008Loaded0.760.080.011
009Unloaded3.760.250.011
Table A4 gathers the converged amplitudes. The amplitudes obtained with normalized force or stroke signals are very close. In both cases, the temporal amplitudes are greater than 1, as the DVC displacement field of the last loading cycle prior to specimen failure is used as spatial mode, while the 4D characterization is performed beyond this stage (i.e., up to specimen failure). Despite this extrapolation, the displacement fields for all analyzed radiographs converge and are consistent. When using 11 plateau-like functions, the amplitudes are close to 1 for all loading steps, except for the last one where the amplitude is greater than 1. Once again, this result stems from the fact that the DVC displacement field of scan 007 (Figure 2) is used as spatial component of the last mode, whereas 4D kinematic measurements are performed up to specimen failure (beyond scan 007). The displacement fields measured with P-DVC beyond scan 009 converge and are consistent. From the reported data, it is concluded that the kinematics is very well captured with the chosen temporal sampling, even though the displacement field becomes more complex with the increase in load level. Based on the reported data, 11 plateau-like functions constrained to 11 spatial modes obtained with FE-DVC are chosen as a suited temporal basis (Section 3).
Table A4. Temporal amplitudes. The * highlights an extrapolation, namely, when employing 11-plateau like functions as temporal basis, the temporal interpolation function describing the last loading stage prior to failure was constrained to the spatial mode of scan 007. When using the measured uniaxial force or stroke signals, only the displacement field of scan 007 is considered as the spatial mode.
Table A4. Temporal amplitudes. The * highlights an extrapolation, namely, when employing 11-plateau like functions as temporal basis, the temporal interpolation function describing the last loading stage prior to failure was constrained to the spatial mode of scan 007. When using the measured uniaxial force or stroke signals, only the displacement field of scan 007 is considered as the spatial mode.
Fu σ
ScanMechanical State α α α
000Unloaded 1.005
001Loaded 1.007
002Unloaded 0.983
003Loaded 1.007
004Unloaded 1.012
005Loaded 1.002
006Unloaded 0.998
007Loaded 1.011
008Loaded 1.006
009Unloaded 0.996
007 *Loaded1.2061.1971.128

Appendix C. Validation of the Methodology for Quantification of Damage Growth

To validate the damage detection procedure (Section 4), an element within which damage does not occur is analyzed in the sequel. It corresponds to element #3, which is located at the specimen bottom surface (Figure A7).
Figure A7. DVC residuals of scan 007. The red contour denotes element #3.
Figure A7. DVC residuals of scan 007. The red contour denotes element #3.
Materials 16 06300 g0a7
The DVC residuals within element #3 are projected onto the detector plane (Figure A8a,b). The first interval of residual increase for the loaded scans is within the range θ = 50 ° to 110°. The reported values are lower than those obtained within element #89 (Figure 7). There is an increase in correlation residuals with respect to levels of scan 000. However, it is not as pronounced as for element #89. It is also important to note that the residual levels obtained for the scans acquired in the loaded steps lie within the same range (Figure A8a). The values obtained for unloaded scans are within the range of scan 000 (Figure A8b). When plotting the rms DVC residuals of element #3 during the entire turns of scan acquisitions, the correlation residuals remain within the same range (Figure A8c). The projected residual map within element #3 calculated for scan 007 when θ = 90 is shown in Figure A8d. Their distribution is homogeneous, and the levels are low. It is concluded that damage did not occur within this element. This hypothesis will also be inspected with the projection residuals.
Figure A8. (a) Rms DVC residuals for element #3 projected onto the detector plane computed for the scans acquired at the peak of the loading cycles (Figure 2b). (b) Corresponding residuals computed for the scans acquired at the unloaded steps (Figure 2b). (c) Rms DVC residuals of element #3 during the entire turns of scan acquisitions. (d) Projected DVC residuals of element #3 when θ = 90 .
Figure A8. (a) Rms DVC residuals for element #3 projected onto the detector plane computed for the scans acquired at the peak of the loading cycles (Figure 2b). (b) Corresponding residuals computed for the scans acquired at the unloaded steps (Figure 2b). (c) Rms DVC residuals of element #3 during the entire turns of scan acquisitions. (d) Projected DVC residuals of element #3 when θ = 90 .
Materials 16 06300 g0a8
As the first interval of DVC residual increase is within the range θ = 50 ° to 110° (Figure A8), the projection residuals are thus only inspected for specific angles, namely, θ = 0 , 42°, 60° and 90° acquired during the loading history. Figure A9 shows the amplification factor for these angles. For all observed angles, the values for the whole loading history are within the same range. There are no significant increases in amplification levels that may highlight damage growth, except for the very last part of the experiment.
Figure A9. Amplification factor of projection residuals within element #3 when (a) θ = 0 °, (b) θ = 42 °, (c) θ = 60 ° and (d) θ = 90 °.
Figure A9. Amplification factor of projection residuals within element #3 when (a) θ = 0 °, (b) θ = 42 °, (c) θ = 60 ° and (d) θ = 90 °.
Materials 16 06300 g0a9
The SVD procedure is also applied to the projection residuals for each analyzed angle ( θ = 0 °, 42°, 60°, and 90°) as explained in Section 4. The singular values are shown in Figure A10. For all four analyzed angles, the first singular value is higher than the second one, thereby indicating a good separation between them. Thus, only the first temporal mode of each angle is analyzed.
Figure A11 shows the first temporal modes computed for θ = 0 °, 42°, 60°, and 90°. When θ = 42 ° and 90°, the levels are within the range of the unloaded stage (here marked with turn 1, corresponding to the range of scan 000). However, increases in levels compared to the first turn are observed when θ = 0 ° and 60°. This trend is due to the fact that these angles may be impacted by damage occurring in other zones/elements. Still, these increases are deemed not significant. The reported results demonstrate that, in this specific case, the damage analysis is highly dependent on the analyzed angles.
Last, the first spatial modes of the projection residuals computed for the chosen angles are shown in Figure A12. For all analyzed angles, the spatial distribution within element #3 is homogeneous.
Figure A10. Normalized singular values for element #3 when (a) θ = 0 °, (b) θ = 42 °, (c) θ = 60 °, and (d) θ = 90 °.
Figure A10. Normalized singular values for element #3 when (a) θ = 0 °, (b) θ = 42 °, (c) θ = 60 °, and (d) θ = 90 °.
Materials 16 06300 g0a10
Figure A11. The first temporal mode of projection residuals within element #3 computed when (a) θ = 0 °, (b) θ = 42 °, (c) θ = 60 °, and (d) θ = 90 °.
Figure A11. The first temporal mode of projection residuals within element #3 computed when (a) θ = 0 °, (b) θ = 42 °, (c) θ = 60 °, and (d) θ = 90 °.
Materials 16 06300 g0a11
Figure A12. First spatial mode of the projection residuals within element #3 computed when (a) θ = 0 °, (b) θ = 42 °, (c) θ = 60 °, and (d) θ = 90 °.
Figure A12. First spatial mode of the projection residuals within element #3 computed when (a) θ = 0 °, (b) θ = 42 °, (c) θ = 60 °, and (d) θ = 90 °.
Materials 16 06300 g0a12

References

  1. Salvo, L.; Cloetens, P.; Maire, E.; Zabler, S.; Blandin, J.J.; Buffière, J.Y.; Ludwig, W.; Boller, E.; Bellet, D.; Josserond, C. X-ray micro-tomography an attractive characterisation technique in materials science. Nucl. Instrum. Methods Phys. Res. Sect.-Beam Interact. Mater. At. 2003, 200, 273–286. [Google Scholar] [CrossRef]
  2. Maire, E.; Buffière, J.Y.; Salvo, L.; Blandin, J.J.; Ludwig, W.; Létang, J.M. On the application of x-ray microtomography in the field of materials science. Adv. Eng. Mater. 2001, 3, 539–546. [Google Scholar] [CrossRef]
  3. Bay, B. Methods and applications of digital volume correlation. J. Strain Anal. 2008, 43, 745–760. [Google Scholar] [CrossRef]
  4. Buljac, A.; Jailin, C.; Mendoza, A.; Neggers, J.; Taillandier-Thomas, T.; Bouterf, A.; Smaniotto, B.; Hild, F.; Roux, S. Digital Volume Correlation: Review of Progress and Challenges. Exp. Mech. 2018, 58, 661–708. [Google Scholar] [CrossRef]
  5. Maire, E.; Bourlot, C.L.; Adrien, J.; Mortensen, A.; Mokso, R. 20 Hz X-ray tomography during an in situ tensile test. Int. J. Fract. 2016, 200, 3–12. [Google Scholar] [CrossRef]
  6. García-Moreno, F.; Kamm, P.H.; Neu, T.R.; Bülk, F.; Mokso, R.; Schlepütz, C.M.; Stampanoni, M.F.M.; Banhart, J. Using X-ray tomoscopy to explore the dynamics of foaming metal. Nat. Commun. 2019, 10, 3762. [Google Scholar] [CrossRef]
  7. Leclerc, H.; Roux, S.; Hild, F. Projection Savings in CT-based Digital Volume Correlation. Exp. Mech. 2015, 55, 275–287. [Google Scholar] [CrossRef]
  8. Jailin, C.; Bouterf, A.; Poncelet, M.; Roux, S. In situ μ CT-scan Mechanical Tests: Fast 4D Mechanical Identification. Exp. Mech. 2017, 57, 1327–1340. [Google Scholar] [CrossRef]
  9. Jailin, C.; Buljac, A.; Bouterf, A.; Hild, F.; Roux, S. Fast four-dimensional tensile test monitored via X-ray computed tomography: Single projection–based digital volume correlation dedicated to slender samples. J. Strain Anal. Eng. Des. 2018, 53, 473–484. [Google Scholar] [CrossRef]
  10. Prashanth, S.; Km, S.; Nithin, K.; Sachhidan, S. Fiber Reinforced Composites—A Review. J. Mater. Sci. Eng. 2017, 6, 1–6. [Google Scholar]
  11. Cosmi, F.; Ravalico, C. Threshold Identification for Micro-Tomographic Damage Characterisation in a Short-Fibre-Reinforced Polymer. Strain 2015, 51, 171–179. [Google Scholar] [CrossRef]
  12. Sket, F.; Enfedaque, A.; Lopez, C.D.; González, C.; Molina-Aldareguia, J.M.; Llorca, J. X-ray computed tomography analysis of damage evolution in open hole carbon fiber-reinforced laminates subjected to in-plane shear. Compos. Sci. Technol. 2016, 133, 40–50. [Google Scholar] [CrossRef]
  13. Rolland, H.; Saintier, N.; Wilson, P.; Merzeau, J.; Robert, G. In situ X-ray tomography investigation on damage mechanisms in short glass fibre reinforced thermoplastics: Effects of fibre orientation and relative humidity. Compos. Part B-Eng. 2017, 109, 170–186. [Google Scholar]
  14. Ortiz-Morales, A.M.; Hanhan, I.; Solano, J.J.; Sangid, M.D. Observing progressive damage in carbon fiber epoxy laminate composites via 3D in situ X-ray tomography. Eng. Fract. Mech. 2021, 246, 107626. [Google Scholar] [CrossRef]
  15. Buljac, A.; Navas, V.M.T.; Shakoor, M.; Bouterf, A.; Neggers, J.; Bernacki, M.; Bouchard, P.; Morgeneyer, T.; Hild, F. On the calibration of elastoplastic parameters at the microscale via X-ray microtomography and digital volume correlation for the simulation of ductile damage. Eur. J. Mech. A-Solids 2018, 72, 287–297. [Google Scholar] [CrossRef]
  16. Tomičević, Z.; Bouterf, A.; Surma, R.; Hild, F. Damage observation in glass fiber reinforced composites via μ-tomography. Mater. Today Proc. 2019, 12, 185–191. [Google Scholar] [CrossRef]
  17. Brault, R.; Germaneau, A.; Dupré, J.; Doumalin, P.; Mistou, S.; Fazzini, M. In-situ Analysis of Laminated Composite Materials by X-ray Micro-Computed Tomography and Digital Volume Correlation. Exp. Mech. 2013, 53, 1143–1151. [Google Scholar] [CrossRef]
  18. Croom, B.; Jin, H.; Mills, B.E.; Carroll, J.; Long, K.; Brown, J.A.; Li, X. Damage mechanisms in elastomeric foam composites: Multiscale X-ray computed tomography and finite element analyses. Compos. Sci. Technol. 2019, 169, 195–202. [Google Scholar] [CrossRef]
  19. Wang, K.; Pei, S.; Li, Y.; Li, J.; Zeng, D.; Su, X.; Xiao, X.; Chen, N. In-situ 3D fracture propagation of short carbon fiber reinforced polymer composites. Compos. Sci. Technol. 2019, 182, 107788. [Google Scholar]
  20. Roux, S.; Hild, F.; Viot, P.; Bernard, D. Three dimensional image correlation from X-ray computed tomography of solid foam. Compos. Part A-Appl. Sci. Manuf. 2008, 39, 1253–1265. [Google Scholar] [CrossRef]
  21. Réthoré, J.; Tinnes, J.; Roux, S.; Buffière, J.; Hild, F. Extended three-dimensional digital image correlation (X3D-DIC). Comptes Rendus Mécanique 2008, 336, 643–649. [Google Scholar] [CrossRef]
  22. Leclerc, H.; Périé, J.; Roux, S.; Hild, F. Voxel-scale digital volume correlation. Exp. Mech. 2011, 51, 479–490. [Google Scholar] [CrossRef]
  23. Réthoré, J.; Limodin, N.; Buffière, J.; Hild, F.; Ludwig, W.; Roux, S. Digital volume correlation analyses of synchrotron tomographic images. J. Strain Anal. Eng. Des. 2011, 46, 683–695. [Google Scholar] [CrossRef]
  24. Tomičević, Z.; Kodvanj, J.; Hild, F. Characterization of the nonlinear behavior of nodular graphite cast iron via inverse identification. Analysis of uniaxial tests. Eur. J. Mech.-A/Solids 2016, 59, 330–343. [Google Scholar] [CrossRef]
  25. Vrgoč, A.; Tomičević, Z.; Smaniotto, B.; Hild, F. Application of different imaging techniques for the characterization of damage in fiber reinforced polymer. Compos. Part A Appl. Sci. Manuf. 2021, 150, 106576. [Google Scholar] [CrossRef]
  26. Maurer, J.; Jerabek, M.; Salaberger, D.; Thor, M.; Kastner, J.; Major, Z. Stress relaxation behaviour of glass fibre reinforced thermoplastic composites and its application to the design of interrupted in situ tensile tests for investigations by X-ray computed tomography. Polym. Test. 2022, 109, 107551. [Google Scholar] [CrossRef]
  27. Vrgoč, A.; Tomičević, Z.; Smaniotto, B.; Hild, F. Characterization of glass fiber reinforced polymer via Digital Volume Correlation: Investigation of notch sensitivity. Mech. Mater. 2023, 177, 104552. [Google Scholar] [CrossRef]
  28. Kosin, V.; Fau, A.; Jailin, C.; Smaniotto, B.; Wick, T.; Hild, F. A projection-based approach to extend digital volume correlation for 4D spacetime measurements. Comptes Rendus Mécanique 2023, 351, 265–280. [Google Scholar] [CrossRef]
  29. Leclerc, H.; Neggers, J.; Mathieu, F.; Hild, F.; Roux, S. Correli 3.0; IDDN.FR.001.520008.000.S.P.2015.000.31500; Agence de Protection des Programmes: Paris, France, 2015. [Google Scholar]
  30. Van Aarle, W.; Palenstijn, W.J.; De Beenhouwer, J.; Altantzis, T.; Bals, S.; Batenburg, K.J.; Sijbers, J. The ASTRA Toolbox: A platform for advanced algorithm development in electron tomography. Ultramicroscopy 2015, 157, 35–47. [Google Scholar] [CrossRef]
  31. Kak, A.; Slaney, M. Principles of Computerized Tomographic Imaging; IEEE Press: New York, NY, USA, 1988. [Google Scholar]
  32. Jailin, C.; Buljac, A.; Bouterf, A.; Poncelet, M.; Hild, F.; Roux, S. Self-calibration for lab-μCT using space-time regularized projection-based DVC and model reduction. Meas. Sci. Technol. 2018, 29, 024003. [Google Scholar] [CrossRef]
  33. Hallo, G.; Lacombe, C.; Parreault, R.; Roquin, N.; Donval, T.; Lamaignère, L.; Néauport, J.; Hild, F. Sub-pixel detection of laser-induced damage and its growth on fused silica optics using registration residuals. Opt. Express 2021, 29 22, 35820–35836. [Google Scholar] [CrossRef]
  34. Vrgoč, A.; Tomičević, Z.; Smaniotto, B.; Hild, F. Characterization of glass fiber reinforced polymer via digital volume correlation: Quantification of strain activity and damage growth. Compos. Sci. Technol. 2023, 234, 109932. [Google Scholar] [CrossRef]
Figure 1. Optical micrograph of the glass fiber mat reinforced polyester resin composite studied herein.
Figure 1. Optical micrograph of the glass fiber mat reinforced polyester resin composite studied herein.
Materials 16 06300 g001
Figure 2. (a) Geometry of the investigated single notched dogbone specimen together with the Region of Interest (depicted with red contour) and finite element mesh employed herein (superimposed over the reconstructed reference volume). The dimensions are expressed in mm. The size of Region of Interest was 384 × 384 × 486 vx, while the physical length of one voxel was 58 µm. (b) Measured uniaxial force (blue) and stroke history (green) of the studied in situ tensile test. The black dots mark the acquisition of the high-quality reference scan (0). The red dots depict the load levels at which full volumes were reconstructed. These volumes were employed in the scan-wise DVC analyses.
Figure 2. (a) Geometry of the investigated single notched dogbone specimen together with the Region of Interest (depicted with red contour) and finite element mesh employed herein (superimposed over the reconstructed reference volume). The dimensions are expressed in mm. The size of Region of Interest was 384 × 384 × 486 vx, while the physical length of one voxel was 58 µm. (b) Measured uniaxial force (blue) and stroke history (green) of the studied in situ tensile test. The black dots mark the acquisition of the high-quality reference scan (0). The red dots depict the load levels at which full volumes were reconstructed. These volumes were employed in the scan-wise DVC analyses.
Materials 16 06300 g002
Figure 3. (a) Root mean square (rms) projection residuals evaluated for the first 360° revolution (scan 0). Initial projection residual field for θ = 0 (b) and θ = 90 (c).
Figure 3. (a) Root mean square (rms) projection residuals evaluated for the first 360° revolution (scan 0). Initial projection residual field for θ = 0 (b) and θ = 90 (c).
Materials 16 06300 g003
Figure 4. Temporal modes consisting of 11 plateau-like functions. The black dashed line represents the measured force signal normalized by its maximum value.
Figure 4. Temporal modes consisting of 11 plateau-like functions. The black dashed line represents the measured force signal normalized by its maximum value.
Materials 16 06300 g004
Figure 5. Change of rms projection residuals over the entire Region of Interest during the iterative procedure.
Figure 5. Change of rms projection residuals over the entire Region of Interest during the iterative procedure.
Materials 16 06300 g005
Figure 6. (a) DVC residuals of scan 007. The red contour denotes element #89. (b) Finite elements containing the highest DVC residual levels. (c) DVC gray level residuals laid over the mesostructure section of the front specimen surface. The red circle marks the area exhibiting elevated correlation residuals within element #89.
Figure 6. (a) DVC residuals of scan 007. The red contour denotes element #89. (b) Finite elements containing the highest DVC residual levels. (c) DVC gray level residuals laid over the mesostructure section of the front specimen surface. The red circle marks the area exhibiting elevated correlation residuals within element #89.
Materials 16 06300 g006
Figure 7. (a) Rms DVC residuals of element #89 projected onto the detector plane computed for the scans acquired at the peak of the loading cycles (Figure 2b). (b) Corresponding residuals computed for the scans acquired at the unloaded steps (Figure 2b). (c) Rms DVC residuals of element #89 during the entire turns for scan acquisitions. (d) Projected DVC residuals of element #89 when θ = 75 . The red arrow marks the crack visible for this specific angle.
Figure 7. (a) Rms DVC residuals of element #89 projected onto the detector plane computed for the scans acquired at the peak of the loading cycles (Figure 2b). (b) Corresponding residuals computed for the scans acquired at the unloaded steps (Figure 2b). (c) Rms DVC residuals of element #89 during the entire turns for scan acquisitions. (d) Projected DVC residuals of element #89 when θ = 75 . The red arrow marks the crack visible for this specific angle.
Materials 16 06300 g007
Figure 8. Amplification of projection residuals within element #89 when (a) θ = 0 °, (b) θ = 56 °, (c) θ = 75 ° and (d) θ = 105 °.
Figure 8. Amplification of projection residuals within element #89 when (a) θ = 0 °, (b) θ = 56 °, (c) θ = 75 ° and (d) θ = 105 °.
Materials 16 06300 g008
Figure 9. Normalized singular values for element #89 when (a) θ = 0 °, (b) θ = 56 °, (c) θ = 75 ° and (d) θ = 105 °.
Figure 9. Normalized singular values for element #89 when (a) θ = 0 °, (b) θ = 56 °, (c) θ = 75 ° and (d) θ = 105 °.
Materials 16 06300 g009
Figure 10. First temporal mode of projection residuals within element #89 computed when (a) θ = 0 °, (b) θ = 56 °, (c) θ = 75 ° and (d) θ = 105 °.
Figure 10. First temporal mode of projection residuals within element #89 computed when (a) θ = 0 °, (b) θ = 56 °, (c) θ = 75 ° and (d) θ = 105 °.
Materials 16 06300 g010
Figure 11. First spatial mode of projection residuals within element #89 calculated when (a) θ = 0 °, (b) θ = 56 °, (c) θ = 75 ° and (d) θ = 105 °.
Figure 11. First spatial mode of projection residuals within element #89 calculated when (a) θ = 0 °, (b) θ = 56 °, (c) θ = 75 ° and (d) θ = 105 °.
Materials 16 06300 g011
Table 1. Converged temporal amplitudes obtained when using 11 plateau-like functions as the temporal basis. The * highlights an extrapolation, namely, the temporal interpolation function describing the last loading stage prior to failure is related to the spatial mode of scan 007.
Table 1. Converged temporal amplitudes obtained when using 11 plateau-like functions as the temporal basis. The * highlights an extrapolation, namely, the temporal interpolation function describing the last loading stage prior to failure is related to the spatial mode of scan 007.
ScanMechanical State α
000Unloaded1.005
001Loaded1.007
002Unloaded0.983
003Loaded1.007
004Unloaded1.012
005Loaded1.002
006Unloaded0.998
007Loaded1.011
008Loaded1.006
009Unloaded0.996
007 *Loaded1.128
Table 2. Comparison of rms differences expressed in vx in z-direction between the nodal displacements measured with FE-DVC and projection-based measurements. The comparison is performed for the converged P-DVC displacement fields for the angles at the beginning of each constant stroke stage.
Table 2. Comparison of rms differences expressed in vx in z-direction between the nodal displacements measured with FE-DVC and projection-based measurements. The comparison is performed for the converged P-DVC displacement fields for the angles at the beginning of each constant stroke stage.
ScanMechanical Statez
000Unloaded0.005
001Loaded0.007
002Unloaded0.01
003Loaded0.007
004Unloaded0.012
005Loaded0.002
006Unloaded0.011
007Loaded0.012
008Loaded0.011
009Unloaded0.011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mandić, A.; Kosin, V.; Jailin, C.; Tomičević, Z.; Smaniotto, B.; Hild, F. Damage Detection in a Polymer Matrix Composite from 4D Displacement Field Measurements. Materials 2023, 16, 6300. https://doi.org/10.3390/ma16186300

AMA Style

Mandić A, Kosin V, Jailin C, Tomičević Z, Smaniotto B, Hild F. Damage Detection in a Polymer Matrix Composite from 4D Displacement Field Measurements. Materials. 2023; 16(18):6300. https://doi.org/10.3390/ma16186300

Chicago/Turabian Style

Mandić, Ana, Viktor Kosin, Clément Jailin, Zvonimir Tomičević, Benjamin Smaniotto, and François Hild. 2023. "Damage Detection in a Polymer Matrix Composite from 4D Displacement Field Measurements" Materials 16, no. 18: 6300. https://doi.org/10.3390/ma16186300

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