Digital Volume Correlation Analysis Research Articles

Sensitivity Study Using Synthetic 3D Image Datasets to Investigate the Effect of Noise Artefacts on Digital Volume Correlation

BackgroundThe potential effect of image noise artefacts on Digital Volume Correlation (DVC) analysis has not been thoroughly studied and, more particularly quantified, even though DVC is an emerging technique widely used in life and material science over the last decade.ObjectiveThis paper presents the results of a sensitivity study to shed light on the effect of various noise artefacts on the full-field kinematic fields generated by DVC, both in zero and rigid body motion.MethodsVarious noise artefacts were studied, including the Gaussian, Salt & Pepper, Speckle noise and embedded Ring Artefacts. A noise-free synthetic microstructure was generated using Discrete Element Modelling (DEM), representing an idealistic case, and acting as the reference dataset for the DVC analysis. Noise artefacts of various intensities (including selected extreme cases) were added to the reference image datasets using MATLAB (R2022) to form the outline of the parametric study. DVC analyses were subsequently conducted employing AVIZO (Thermo Fisher). A subset-based local approach was adopted. A three-dimensional version of the Structural Similarity Index Measure (SSIM) was used to define the similarity between the compared image datasets on each occasion. Sub-pixel rigid body motion was applied on the DEM-generated microstructure and subsequently “poisoned” with noise artefacts to evaluate mean bias and random error of the DVC analysis.ResultsWhen the local approach is implemented, the sensitivity study on zero motion data revealed the insignificant effect of the Gaussian, Salt & Pepper, and Speckle noise on the DVC-computed kinematic field. Therefore, the presence of such noise artefacts can be neglected when DVC is executed. On the contrary, Ring Artefacts can pose a considerable challenge and therefore, DVC results need to be evaluated cautiously. A linear relationship between SSIM and the correlation index is observed for the same noise artefacts. Gaussian noise has a pronounced effect on the mean bias error for sub-pixel rigid body motion.ConclusionsGenerating synthetic image datasets using DEM enabled the investigation of a variety of noise artefacts that potentially affect a DVC analysis. Given that, any microstructure – resembling the material studied – can be simulated and used for a DVC sensitivity analysis, supporting the user in appropriately evaluating the computed kinematic field. Even though the study is conducted for a two-phase material, the method elaborated in this paper also applies to heterogeneous multi-phase materials also. The conclusions drawn are valid within the environment of the AVIZO DVC extension module. Alternative DVC algorithms, utilising different approaches for the cross-correlation and the sub-pixel interpolation methods, need to be investigated.

Full-field strain distribution in hierarchical electrospun nanofibrous poly-L(lactic) acid/collagen scaffolds for tendon and ligament regeneration: A multiscale study

Regeneration of injured tendons and ligaments (T/L) is a worldwide need. In this study electrospun hierarchical scaffolds made of a poly-L (lactic) acid/collagen blend were developed reproducing all the multiscale levels of aggregation of these tissues. Scanning electron microscopy, microCT and tensile mechanical tests were carried out, including a multiscale digital volume correlation analysis to measure the full-field strain distribution of electrospun structures. The principal strains (εp1 and εp3) described the pattern of strains caused by the nanofibers rearrangement, while the deviatoric strains (εD) revealed the related internal sliding of nanofibers and bundles. The results of this study confirmed the biomimicry of such electrospun hierarchical scaffolds, paving the way to further tissue engineering and clinical applications.

Quantifying the immediate post-implantation strain field of cadaveric tibiae implanted with cementless tibial trays: A time-elapsed micro-CT and digital volume correlation analysis during stair descent

Primary stability, the mechanical fixation between implant and bone prior to osseointegration, is crucial for the long-term success of cementless tibial trays. However, little is known about the mechanical interplay between the implant and bone internally, as experimental studies quantifying internal strain are limited. This study employed digital volume correlation (DVC) to quantify the immediate post-implantation strain field of five cadaveric tibiae implanted with a commercially available cementless titanium tibial tray (Attune, DePuy Synthes). The tibiae were subjected to a five-step loading sequence (0–2.5 bodyweight, BW) replicating stair descent, with concomitant time-elapsed micro-CT imaging. With progressive loads, increased compression of trabecular bone was quantified, with the highest strains directly under the posterior region of the tibial tray implant, dissipating with increasing distance from the bone-implant interface. After load removal of the last load step (2.5BW), residual strains were observed in all of the five tibiae, with residual strains confined within 3.14 mm from the bone-implant interface. The residual strain is reflective of the observed initial migration of cementless tibial trays reported in clinical studies. The presence of strains above the yield strain of bone accepted in literature suggests that inelastic properties should be included within finite element models of the initial mechanical environment. This study provides a means to experimentally quantify the internal strain distribution of human tibia with cementless trays, increasing the understanding of the mechanical interaction between bone and implant.

Adduction induces large optic nerve head deformations in subjects with normal-tension glaucoma

PurposeTo assess intraocular pressure (IOP)-induced and gaze-induced optic nerve head (ONH) strains in subjects with high-tension glaucoma (HTG) and normal-tension glaucoma (NTG).DesignClinic-based cross-sectional study.MethodsThe ONH from one eye of 228...

Damage characteristics and constitutive model of phosphogypsum/fly ash/slag recycled aggregate concrete under uniaxial compression

The recycling of construction and industrial waste is essential for environmental protection and carbon emissions reduction. A phosphogypsum/fly ash/slag (PFS) composite was therefore employed in this study to fully replace ordinary Portland cement, producing PFS recycled aggregate concrete (PFSRAC). The mechanical properties of various PFSRAC mixes were first investigated using a series of uniaxial compression tests, showing the combination of PFS and recycled aggregate has an obvious effect on the compressive strength of PFSRAC. Next, the change in the mean and variance of grayscale values in cross-sectional images obtained by in-situ X-ray computed tomography were found to effectively reflect the damage characteristics of PFSRAC under uniaxial compression. A digital volume correlation analysis of these images indicated no deformation localization area in the interfacial transition zone (ITZ) between the PFS matrix and recycled aggregate, demonstrating effective densification of the ITZ that increased the compressive strength of the PFSRAC. Finally, the constitutive model for PFSRAC was established using a strain statistics analysis that was shown to effectively reflect the observed stress–strain relationship.

Effect of radiation-induced damage of trabecular bone tissue evaluated using indentation and digital volume correlation

Exposure to X-ray radiation for an extended amount of time can cause damage to the bone tissue and therefore affect its mechanical properties. Specifically, high-resolution X-ray Computed Tomography (XCT), in both synchrotron and lab-based systems, has been employed extensively for evaluating bone micro-to-nano architecture. However, to date, it is still unclear how long exposures to X-ray radiation affect the mechanical properties of trabecular bone, particularly in relation to lab-XCT systems. Indentation has been widely used to identify local mechanical properties such as hardness and elastic modulus of bone and other biological tissues. The purpose of this study is therefore, to use indentation and XCT-based investigative tools such as digital volume correlation (DVC) to assess the microdamage induced by long exposure of trabecular bone tissue to X-ray radiation and how this affects its local mechanical properties. Trabecular bone specimens were indented before and after X-ray exposures of 33 and 66 h, where variation of elastic modulus was evaluated at every stage. The resulting elastic modulus was decreased, and micro-cracks appeared in the specimens after the first long X-ray exposure and crack formation increased after the second exposure. High strain concentration around the damaged tissue exceeding 1% was also observed from DVC analysis. The outcomes of this study show the importance of designing appropriate XCT-based experiments in lab systems to avoid degradation of the bone tissue mechanical properties due to radiation and these results will help to inform future studies that require long X-ray exposure for in situ experiments or generation of reliable subject-specific computational models.

Strain localization in 4D imaged by X-ray computed tomography and digital volume correlation: Discrete compaction bands in Leitha limestone

Understanding compaction localization in porous limestone in the laboratory is significantly more challenging than in sandstone because of the lack of consistent acoustic emission activity in carbonate samples. Previous studies have therefore relied on X-ray Computed Tomography imaging (CT). The first unambiguous evidence of compaction band development in limestone was provided by Huang et al. (2019), who performed synchrotron in situ CT imaging during shear-enhanced compaction in a sample of Leitha limestone. In this study, we analyzed this data set using Digital Volume Correlation (DVC). Our new results confirm that inelastic compaction occurred in two stages in Leitha limestone: macropore collapse first and then sequential growth of compaction bands. The displacement field inferred from DVC revealed that the bands in Leitha limestone showed mostly normal displacement discontinuities, as expected for compaction bands. DVC analysis also gave more constraints on band geometric attributes, as well as on the spatiotemporal development of strain. Analysis of the autocorrelation function for the volumetric strain suggests spatiotemporal correlation of strain precursory to the band's formation, their width, spacing and their complex sequential growth.

In situ micromechanical analysis of discontinuous fiber-reinforced composite material based on DVC strain and fiber orientation fields

The influence of fiber orientation and its distribution on the global and local mechanical behavior of a short fiber-reinforced thermoplastic (SFT) composite is studied based on an in situ tensile test using synchrotron X-ray tomography. The in situ test data are utilized to compute three-dimensional (3D) strain fields inside the material at various loading steps through an in-house digital volume correlation (DVC) code. The DVC analysis finds that, although the global strain fields seem relatively uniform, the internal strain fields are locally non-uniform because of the complex reinforcement architecture. The local non-uniformity is further investigated at a single fiber level to find a correlation between the orientation of fibers, fiber volume fraction, and strain field. The investigation is performed on micro-scale subvolumes, in which all the fibers are individually segmented from the 3D tomographic data. It is found that the fibers aligned in the direction transverse to loading impede transversal shrinkage, resulting in a locally small compressive strain in that direction. The in situ test has revealed that the average fiber orientation does not change noticeably even after the global strain reaches 18%. This small change implies that initially-determined fiber orientation is crucial to the locally anisotropic mechanical behavior of a SFT-like material even undergoing large deformation. On the other hand, the loading-directional strain is more affected by the local volume fraction than the fiber orientation. The local volume fraction also mainly determines the stiffness of the subvolume that can be estimated from an analytical homogenization model. The effect of the local volume fraction on the subvolume stiffness is identified from multivariate analysis.

Application of phase-field fracture theories and digital volume correlation to synchrotron X-ray monitored fractures in human trabecular bone: A case study

Fracture processes of trabecular bone have been studied using various approaches over the years. However, reliable methods to analyse fracture at the single trabecula level are limited. In this study, a digital volume correlation (DVC) and a phase-field fracture model are applied and contrasted for human trabecular bone to analyse its failure under global compression at high resolution.A human trabecular bone sample was fractured in situ under synchrotron-based X-ray micro computed tomography (CT). Reconstructed CT data was then used in DVC algorithms to obtain high-resolution displacement fields in the bone at different load steps. A high-resolution specimen-specific structural mesh was discretized from the CT data and used for the phase-field simulation of the fracturing bone.The DVC analysis showed opening mode cracks as well as shear mode cracks. Strains in cracked regions were analysed. The load distribution in the trabecular structure resulted in two completely separated fracture regions in the sample body. A phenomenon that was also captured in the phase-field model. The results encourage us to believe improvements in boundary conditions and material models are worthwhile pursuing. Findings in this study support further development of a phase-field method to analyse fracture in samples with complex morphology, such as trabecular bone, and the capacity of DVC to quantify strains and slowly growing stable fractures during step-wise loading of trabecular bone.

Debonding analysis via digital volume correlation during in-situ pull-out tests on fractal fibers

The quantification of debonding was performed for additively manufactured “fractal” fibers embedded within two brittle matrices. Three pull-out tests were carried out inside of an X-ray tomograph allowing for Digital Volume Correlation analyses. Relative motions at the interfaces were measured thanks to adapted meshes with split nodes. Profiles of normal, tangential and vertical displacement jumps as well as vertical strains in the fibers were used to study interfacial debonding. An articulated load transfer mechanism between the fiber and the matrix was observed in the examined tests, as demonstrated by zigzagged distributions of vertical displacement jumps and vertical strain profiles in the fibers at the initial stages of pull-out. Vertical strain concentrations were observed in correspondence to lateral protrusions (or ribs) of the reinforcing fibers. These results suggest that fiber–matrix interlocking may be affected by geometry-driven tensile stiffening effects between the ribs. For larger values of pull-out displacements, more diffuse damage of the fiber–matrix interface was observed between the ribs, especially in plaster matrices.

Micro-CT scan optimisation for mechanical loading of tibia with titanium tibial tray: A digital volume correlation zero strain error analysis

Primary stability of press-fit tibial trays is achieved by introducing an interference fit between bone and implant. The internal cancellous bone strains induced during this process and during loading have yet to be quantified experimentally. Advancements in large-gantry micro-CT imaging and digital volume correlation (DVC) allow quantification of such strains. However, before undertaking such a test, experimental requirements and DVC performance need to be examined, particularly considering the presence of a large orthopaedic implant (tibial tray). The aim of this study was to assess the DVC zero-strain accuracy (mean absolute error: MAER) and precision (standard deviation of error: SDER) on a cadaveric human tibia implanted with a titanium press-fit tray across four plausible scanning configurations, using a cabinet micro-CT system (Nikon XT H 225 ST). These varied in rotation step and resulting scanning time (106 min vs. 66 min), presence or absence of a 2 mm-thick aluminium cylinder for mechanical testing, and X-ray tube voltage (150 kVp vs. 215 kVp). One proximal tibia was implanted and micro-CT scanned (42 μm/pixel), with repeated scanning and specimen repositioning in between. DVC (DaVis, LaVision, direct correlation) was performed on nine cubic volumes of interest (VOIs: 13.4 mm-side) and across the entire proximal tibia. Strain errors were comparable across the four scanning configurations and sufficiently low for assessing bone within its elastic region in VOIs (MAER=223-540 με; SDER=88-261 με) and at organ level (MAER=536 με; SDER=473 με). Whilst the investigated experimental conditions, including a large titanium implant, present added complexity for DVC analysis, scans of sufficient quality can be achieved, reaching a compromise between the DVC requirements and the wanted application. The approach used for choosing the X-ray source settings considering the transmitted X-ray signal intensity and source power, is also discussed.

Pore and permeability properties of reservoir sandstone under a uniaxial compression CT test

For underground gas storage projects that use depleted gas reservoirs, cyclic gas injection will change the stress fields and cause rock deformation. In-depth knowledge of the stress-induced deformation of reservoir rocks is important for secure gas storage. In this study, we performed a uniaxial compression computed tomography (CT) test with a reservoir sandstone sample to investigate the influence of rock deformation on the pore and permeability properties. Based on scanned 3D images, we analyzed the axial strain evolution by using the digital volume correlation (DVC) technique and the variations in pore properties with the microstructures and permeability properties that were extracted by the pore network modeling method. The DVC analysis shows that strain localization develops gradually with increasing axial stress. The initial porosity distribution with sample height is inhomogeneous. The local porosity variations are not completely consistent with the global porosity variations under increasing axial stress. The level of inconsistency will decrease with increasing axial stress in both the volumetric compression and dilation stages. Increases in axial stress mainly induce the variations in pore numbers and pore volumes. The variations in pore connectivity are not significant. For the permeability evolution, X-ray CT scanning combined with pore network modeling can simulate the permeability decreases during volumetric compression. For high porosity sandstones (>15%), the permeability still decreases in the volumetric dilation stage due to grain size reductions and tortuosity increases. However, this technique has limitations in analyzing the permeability reductions in the volumetric dilation stage. The microcracks that are induced inside the grain particles are too small to be recognized. The variations in the relative gas and water permeabilities at different axial stresses are not significant.

Digital Volume Correlation analyses to study deformation and damage mechanisms of teak in torsion

Digital Volume Correlation analyses to study deformation and damage mechanisms of teak in torsion

Deep learning-based digital volume correlation

With the rapid development of 3D volume imaging technology, digital volume correlation (DVC) has already made great progresses and become an indispensable tool in various engineering kingdoms. Inspired by artificial intelligence-related technologies, we develop a novel deep learning-based DVC methodology for 3D deformation analyses with the aid of a set of well-trained neural networks, hereafter referred to as DVC-Net, to implement accurate, precise and robust 3D deformation characterization in an efficient manner in this work. This proposed DVC-Net comprises three sub-networks, i.e., DVC-NetI, DVC-NetS and DVC-NetD, which are used in turn for integer-voxel searching, sub-voxel registration and displacement field denoising. Unlike traditional DVC methods, e.g., subvolume-based local DVC or finite-element method-based global DVC, the proposed DVC-Net neither explicitly relies on correlation criteria nor requires numerical iterative calculations, which greatly reduces computational complexity of the DVC analyses and therefore improves computational efficiency by at least 2∼3 orders of magnitude. Besides, it is suitable for dealing with various practical problems arising in the DVC-related studies such as volumetric images that experience changes in lighting conditions or noisy images. It is anticipated that the proposed DVC-Net can provide an alternative method for high-performance volumetric deformation characterization.

Synchrotron 4D X-Ray Imaging Reveals Strain Localization at the Onset of System-Size Failure in Porous Reservoir Rocks

Understanding the mechanisms of strain localization leading to brittle failure in reservoir rocks can shed light on geomechanical processes such as porosity and permeability evolution during rock deformation, induced seismicity, fracturing, and subsidence in geological reservoirs. We perform triaxial compression tests on three types of porous reservoir rocks to reveal the local deformation mechanisms that control system-size failure. We deformed cylindrical samples of Adamswiller sandstone (23% porosity), Bentheim sandstone (23% porosity), and Anstrude limestone (20% porosity), using an X-ray transparent triaxial deformation apparatus. This apparatus enables the acquisition of three-dimensional synchrotron X-ray images, under in situ stress conditions. Analysis of the tomograms provide 3D distributions of the microfractures and dilatant pores from which we calculated the evolving macroporosity. Digital volume correlation analysis reveals the dominant strain localization mechanisms by providing the incremental strain components of pairs of tomograms. In the three rock types, damage localized as a single shear band or by the formation of conjugate bands at failure. The porosity evolution closely matches the evolution of the incremental strain components of dilation, contraction, and shear. With increasing confinement, the dominant strain in the sandstones shifts from dilative strain (Bentheim sandstone) to contractive strain (Adamswiller sandstone). Our study also links the formation of compactive shear bands with porosity variations in Anstrude limestone, which is characterized by a complex pore geometry. Scanning electron microscopy images indicate that the microscale mechanisms guiding strain localization are pore collapse and grain crushing in sandstones, and pore collapse, pore-emanated fractures and cataclasis in limestones. Our dynamic X-ray microtomography data brings unique insights on the correlation between the evolutions of rock microstructure, porosity evolution, and macroscopic strain during the approach to brittle failure in reservoir rocks.

High resolution three-dimensional strain measurements in human articular cartilage

An unresolved challenge in osteoarthritis research is characterising the localised intra-tissue mechanical response of articular cartilage. The aim of this study was to explore whether laboratory micro-computed tomography (micro-CT) and digital volume correlation (DVC) permit non-destructive quantification of three-dimensional (3D) strain fields in human articular cartilage. Human articular cartilage specimens were harvested from the knee, mounted into a loading device and imaged in the unloaded and loaded states using a micro-CT scanner. Strain was measured throughout the cartilage volume using the micro-CT image data and DVC analysis. The volumetric DVC-measured strain was within 5% of the known applied strain. Variation in strain distribution between the superficial, middle and deep zones was observed, consistent with the different architecture of the material in these locations. These results indicate DVC method may be suitable for calculating strain in human articular cartilage.

Using In Situ Synchrotron‐Radiation‐Based Microtomography to Investigate 3D Structure‐Dependent Material Properties of Tension Wood

Synchrotron‐radiation‐based microtomography enables the 3D analysis of biological samples in situ beyond simple visualization, providing accurate measurements and recording of temporal data. The microtomography end station at P05 (PETRA III, Hamburg, Germany), operated by the Helmholtz Zentrum Hereon, can accommodate complex sample environments such as a load frame for mechanical testing. Herein, the strain analysis of a volumetric time series dataset of small hardwood samples is presented. For this in situ mechanical testing, the load frame is operated in tensile mode, and the biogenic material samples are subdivided into those specialized for tensile forces and those of regular anatomical structure. Digital volume correlation analysis allows for the prediction of strain development at any position in the sample. The tissue specialized for tensile strength can align dislocations formed under tension such that the deformation is stronger and more ordered. The results show how high‐resolution imaging of sequential loading and the subsequent localization of strain can reveal functional morphological relationships. These methods can be extended to other biological materials and may prove to be extremely relevant for the analyses of fibrous and/or layered composite materials.

INTERNAL SPECKLE PATTERN QUALITY ASSESSMENT AND OPTIMAL SELECTION OF VOXEL POINTS FOR DIGITAL VOLUME CORRELATION

Digital volume correlation (DVC) is a powerful experimental tool for quantitative 3D internal deformation measurement throughout the interior of materials or biological tissues. By comparing the volumetric images acquired at the reference and deformed states by a volumetric imaging facility (e.g., X-ray CT), DVC provides full-field displacements with subvoxel accuracy and full-field strain maps. When using DVC method for internal deformation measurement, the quality of internal speckle pattern of a test sample has an important impact on the accuracy and precision of the measured displacements. In this work, theoretical analysis and numerical simulation experiments are firstly performed to assess the quality of internal speckle patterns. The results show that the displacement measurement error of DVC is negative correlation with the sum of square of subvolume intensity gradient (SSSIG) of a subvolume. Thus, the SSSIG value of a subvolume can be used as a simple and effective metric for quality assessment of its internal speckle pattern. In practical DVC analyses, increasing the SSSIG parameter of a subvolume helps to improve the measurement accuracy and precision of DVC. Despite the SSSIG in a subvolume can be increased by enlarging its subvolume size, this simple approach, however, would lead to increased calculation burden. With a purpose to increase SSSIG without significantly increase the amount of calculation, this paper further proposes a method for optimal selection of calculation voxel points in subvolumes. Specifically, the voxel points with small gray gradients in a subvolume are excluded from DVC calculation to decrease the amount of calculation when increasing the size of the subvolume. The efficacy of the presented voxel point optimal selection method is validated by analyzing computer simulated and experimentally obtained volume images. Experimental results reveal that the proposed voxel point selection method is capable of reducing the calculation burden when increasing the subvolume size to increase the SSSIG parameters.

Digital volume correlation analysis of polylactic acid based fused filament fabrication printed composites.

Fused filament fabrication (FFF) has rapidly begun to see implementation in industrial fields as a method of rapid manufacturing. Traditional FFF parts are made from a single thermoplastic polymer. The polymer is heated to its melting point and deposited on a work bed where a model is gradually built from the base up. While traditional FFF parts have low mechanical properties, a reinforcing phase allows for improved mechanical properties. The addition of a reinforcing material to the base polymer and complex internal microstructure of the 3 D printed party leads to anisotropic mechanical properties. Thus, these materials’ mechanical properties become challenging to characterize using traditional measurement techniques due to the previously mentioned factors. Therefore, it is essential to develop a method in which mechanical properties can be measured and analyzed. This study aims to characterize the mechanical behaviour under a uniaxial tensile load of an FFF produced polylactic acid (PLA)-copper particulate composite. The internal response of the FFF sample was imaged using micro-computed tomography at predetermined loads. The μ-CT images were input into an open-source digital volume correlation (DVC) software to measure the internal displacements and strain tensor fields. The study results show the development of different strain fields and interior features of the FFF parts.

The Brittle‐Ductile Transition in Porous Limestone: Failure Mode, Constitutive Modeling of Inelastic Deformation and Strain Localization

AbstractUnderstanding of the mechanics of the brittle‐ductile transition (BDT) in porous limestone is significantly more challenging than for sandstone because of the lack of consistent acoustic emission activity in limestone, meaning that one must rely on alternative techniques. In this paper, we investigate systematically the failure modes in Indiana limestone using X‐ray microComputed Tomography imaging (μCT) and Digital Volume Correlation (DVC). Our new mechanical data show that the envelope for the onset of shear‐enhanced compaction can be well approximated by an elliptical cap. The DVC analysis revealed the development of shear bands through the BDT, but no evidence of compaction bands. The shear band angles were between 29° and 46° with respect to the maximum principal stress. Compiling these new results with published data on Purbeck and Leitha limestones, we showed that inelastic compaction in each of these dual porosity allochemical limestones was in a good agreement with the normality condition, as defined in plasticity theory. Comparison of the observed failure modes with predictions based on bifurcation analysis showed that the shear band angles are consistently smaller than the theoretical predictions.