Synchrotron Micro-computed Tomography Research Articles

Influence of crystal structure on constitutive anisotropy of silica sand at particle-scale

In spite of recent advanced geotechnical research regarding the influence of particle-scale morphology and fabric on the constitutive behavior of silica sands, there is still no published literature that investigated the influence of crystal structure of sand particles on their constitutive behavior. The main constituent of silica sand particles is quartz mineral, which has a trigonal crystal system and hexagonal lattice. This paper investigates the influence of quartz crystal structure on the constitutive response of silica sand particles using 3D X-ray diffraction (3DXRD), synchrotron micro-computed tomography (SMT), and 3D finite element (FE) analysis. Unconfined uniaxial compression experiments on single particles of synthetic silica cubes exposed the notion of constitutive anisotropy caused by the crystal structure of quartz. Since the silica cubes have simple regular geometry 1mm3, they were loaded in different directions with respect to their crystal local planes. The crystal structure was also identified for a natural silica sand specimen composed of 2705 particles and loaded in confined uniaxial compression. Physical limitations of the experimental measurements made a strong case for the use of 3D FE analysis to study how the change in crystal local orientation of individual sand particles can fundamentally produce anisotropy in the constitutive response of natural sand specimens. The 3D FE analysis was validated to accurately model the crystal-based constitutive anisotropy in silica particles using the results of the silica cube experiments. The analysis was then executed on 3D FE meshes that closely matched the physical shape and fabric of the 2705 particles in the natural sand specimen. Both the experiments and 3D FE analysis showed that the crystal structure of quartz essentially causes a directional anisotropy in the constitutive behavior of silica particles.

Crack initiation and propagation in braided SiC/SiC composite tubes: Effect of braiding angle

Crack initiation and propagation in three braided SiC/SiC composite tubes with different braiding angles are investigated by in situ tensile tests with synchrotron micro-computed tomography. Crack networks are precisely detected after an image subtraction procedure based on Digital Volume Correlation. FFT based simulations are performed on the full-resolution 3D images to assess elastic stress/strain fields. Quantitative measurements of the crack geometries are performed using a novel method based on grey levels. The results show that braiding angle has no obvious effect on the location of crack onsets (initiation always occurs at tow interfaces), whereas it significantly affects the paths of crack propagation. This work provides an explicit demonstration of the crack propagation scenarios with respect to the mesoscopic fibre architectures.

The application of synchrotron micro-computed tomography to characterize the three-dimensional microstructure in irradiated nuclear fuel

This research demonstrated the first time that synchrotron X-ray micro-computed tomography (μ-CT) was utilized to analyze the three-dimensional (3-D) morphology of the microstructure of irradiated nuclear fuel. Both absorption contrast and propagation-based phase contrast enhanced μ-CT techniques were evaluated. Porosity and phase region volumes were determined in neutron irradiated uranium 10 wt% zirconium (U–10Zr) using propagation-based phase contrast enhanced μ-CT. Volumes were calculated using image processing software and algorithms, with the porosity volume subsequently compared to the ASTM E562 method. The 3-D porosity of a small region (∼8 × 105 μm3) located in the intermediate fuel redistribution zone that was removed via a focused ion beam was determined to be 7.2% via image processing. A detailed 3-D analysis of the porosity with respect to volume and morphology was provided. Five growth stages of pore evolution were observed including: nucleation, growth, coalescence, interconnected porosity, and extended/interconnected porosity. Four distinct phase regions were identified, including the pores, via both absorption contrast and propagation-based phase contrast enhanced μ-CT. A quantitative comparison of fuel swelling between this localized microscopic U–10Zr fuel region that was previously irradiated in the Fast Flux Test Facility, and historic U–10Zr fuels irradiated in the Experimental Breeder Reactor-II (EBR-II), was conducted. It was determined that for similar burnups, the localized swelling in this microscopic fuel region was 7.7%, which was approximately one-half of the value of the macroscopic swelling reported in EBR-II irradiated fuel. With the advent of new specimen preparation tools for irradiated materials, as well as advances in characterization techniques, researchers can now obtain previously inaccessible microstructural insight into irradiated fuels. This possibility naturally lends itself to improved fuel performance models for the nuclear community.

3D measurements of hydrate surface area during hydrate dissociation in porous media using dynamic 3D imaging

A better understanding of the kinetics of hydrate dissociation is essential to reliably predict gas production potential from natural hydrate reservoirs. Most hydrate dissociation models assume hydrates to be a constant number of equal-sized spheres dissociating at a constant rate. This paper uses dynamic 3D synchrotron micro-computed tomography (SMT) imaging to study hydrate surface area evolution during Xenon hydrate dissociation. Hydrates are formed inside a high-pressure low-temperature cell filled with partially saturated ASTM 20-30 Ottawa sand. Hydrate dissociation is initiated through depressurization in the first experiment and through thermal stimulation in the second experiment. During dissociation, continuous full 3D SMT images were acquired where each scan took 45 s to complete. A combination of cementing, pore-filling, and surface coating pore habits were observed for the depressurization experiment and pore-filling for the thermal stimulation experiment. Surface coating hydrates dissociate faster than hydrates with pore-filling pore habit due to the higher specific area which allows for more surface for hydrates to dissociate it. Direct measurements of hydrate volume and hydrate surface area suggest that even with a combination of hydrate pore habits formed within the 3D porous media, estimation of hydrate surface area as a linear relation with (hydrate volume)2/3 is best for hydrate saturation less than a threshold value depending on the dissociation method and driving force. (hydrate volume)2 and (hydrate volume)3 were found to better estimate hydrate interfacial area in comparison to (hydrate volume)2/3 for the depressurization experiment and thermal stimulation experiment, respectively.

Dynamic 3D imaging of gas hydrate kinetics using synchrotron computed tomography

The availability of natural gas hydrates and the continuing increase in energy demand, motivated researchers to consider gas hydrates as a future source of energy. Fundamental understanding of hydrate dissociation kinetics is essential to improve techniques of gas production from natural hydrates reservoirs. During hydrate dissociation, bonds between water (host molecules) and gas (guest molecules) break and free gas is released. This paper investigates the evolution of hydrate surface area, pore habit, and tortuosity using in-situ imaging of Xenon (Xe) hydrate formation and dissociation in porous media with dynamic three-dimensional synchrotron microcomputed tomography (SMT). Xe hydrate was formed inside a high- pressure, low-temperature cell and then dissociated by thermal stimulation. During formation and dissociation, full 3D SMT scans were acquired continuously and reconstructed into 3D volume images. Each scan took only 45 seconds to complete, and a total of 60 scans were acquired. Hydrate volume and surface area evolution were directly measured from the SMT scans. At low hydrate saturation, the predominant pore habit was surface coating, while the predominant pore habit at high hydrate saturation was pore filling. A second-degree polynomial can be used to predict variation of tortuosity with hydrate saturation with an R2 value of 0.997.

3D Finite Element Modelling of Fracture of Sand Particles Subjected to High Strain Loading Rate

Fracture is a common failure mode of sand particles when they are subjected to high strain rate (HSR) loading conditions such as blast, impact, or projectile penetration. The constitutive behavior and the failure mode of sand particles are influenced by the loading rate; therefore, a particle-scale constitutive model is necessary to address the effects of HSR loading and particle fracture within a sandy material. To that end, a Kolsky test (i.e., a HSR 1D compression test), 3D X-ray computed tomography (CT), and finite element methods were employed in this study to investigate the failure mode of individual natural sand particles when they are subjected to a HSR loading. Individual particles were first imaged using CT technique followed by testing them using Kolsky bar at an approximate strain rate of $$10^{4} S^{ - 1}$$ . The fragments of fractured particles were collected and imaged using synchrotron micro computed tomography (SMT) for further evaluation of the fracture mechanisms within individual particles. A brittle fracture model was adopted to perform 3D Finite element (FE) modelling to capture the fracture of individual sand particles. 3D CT images of the particles were used to generate 3D meshes with similar morphology as the actual sand particles and the particles were loaded similar to Kolsky experiments. The paper discusses the calibration and the validation of the model and compares the fracture mechanisms within the sand particles based on FE simulations with experimental measurements. The brittle fracture model successfully simulated the fracture mechanisms of the experiments. The effect of loading mechanisms on the mechanisms of the fracture of the particles is discussed in detail.

Silver-Albumin Tissue Staining Protocol to Visualize Odontogenesis in Whole Embryos.

Visualizing tooth organs from their earliest inception as they actually appear in three dimensions has, until recently, been difficult due to the technical obstacle of imaging these tiny, translucent, low-density embryonic craniodental tissues. Related to this obstacle, quantifying craniodental morphology has been confounded by the time consuming need to physically section and then digitally photograph and reconstruct these images of tissues into 3D volumes. Here we provide a simple solution in the form of an overnight silver albumin tissue stain for whole embryos. Because it is differentially absorbed by embryonic tissues, this stain generates the contrast needed to detect and visualize unmineralized dental tissues. Stained specimens can be scanned using either desktop or synchrotron micro-computed tomography systems, generating digital 3D datasets of whole embryos that can immediately be used to assess dental morphology and histology. Craniodental structures can then be measured with high precision and accuracy using 3D image analysis software.

Disentangling isolated dental remains of Asian Pleistocene hominins and pongines.

Scholars have debated the taxonomic identity of isolated primate teeth from the Asian Pleistocene for over a century, which is complicated by morphological and metric convergence between orangutan (Pongo) and hominin (Homo) molariform teeth. Like Homo erectus, Pongo once showed considerable dental variation and a wide distribution throughout mainland and insular Asia. In order to clarify the utility of isolated dental remains to document the presence of hominins during Asian prehistory, we examined enamel thickness, enamel-dentine junction shape, and crown development in 33 molars from G. H. R. von Koenigswald's Chinese Apothecary collection (11 Sinanthropus officinalis [= Homo erectus], 21 “Hemanthropus peii,” and 1 “Hemanthropus peii” or Pongo) and 7 molars from Sangiran dome (either Homo erectus or Pongo). All fossil teeth were imaged with non-destructive conventional and/or synchrotron micro-computed tomography. These were compared to H. erectus teeth from Zhoukoudian, Sangiran and Trinil, and a large comparative sample of fossil Pongo, recent Pongo, and recent human teeth. We find that Homo and Pongo molars overlap substantially in relative enamel thickness; molar enamel-dentine junction shape is more distinctive, with Pongo showing relatively shorter dentine horns and wider crowns than Homo. Long-period line periodicity values are significantly greater in Pongo than in H. erectus, leading to longer crown formation times in the former. Most of the sample originally assigned to S. officinalis and H. erectus shows greater affinity to Pongo than to the hominin comparative sample. Moreover, enamel thickness, enamel-dentine junction shape, and a long-period line periodicity value in the “Hemanthropus peii” sample are indistinguishable from fossil Pongo. These results underscore the need for additional recovery and study of associated dentitions prior to erecting new taxa from isolated teeth.

The effect of growth rate on the three-dimensional orientation of vascular canals in the cortical bone of broiler chickens.

Vascular canals in cortical bone during growth and development typically show an anisotropic pattern with canals falling into three main categories: circumferential, radial, and longitudinal. Two major hypotheses attempt to explain the preferred orientations in bone: that vascular canal orientation is optimized to resist a predominant strain direction from functional loading, or that it reflects growth requirements and velocity. We use a controlled growth experiment in broiler chickens to investigate the effect of growth rate on vascular canal orientation. Using feed restriction we set up a fast growing control group and a slow growing restricted group. We compared the microstructure in the humerus and the femur at 42 days of age using synchrotron micro‐computed tomography (micro‐CT), a three‐dimensional (3D) method that visualizes the full canal network. We measured the 3D orientation of each canal in the whole cross‐section of the bone cortex using a set of custom imagej scripts. Using these orientations we compute laminar, radial, and longitudinal indices that measure the proportion of circumferential, radial, and longitudinal canals, by unit of length, in the cortex. Following previous studies we hypothesized that vascular canal orientation is related to growth, with radial canals linked to a faster growth rate and related to functional loading through a high laminar index in flight bones which reflects torsional loading resulting from active flight. The control group had final body weights that were nearly twice the final weights of the restricted group and higher absolute growth rates. We found consistent patterns in the comparison between the humerus and the femur in both groups, with the humerus having higher laminar and longitudinal indices, and a lower radial index than the femur. The control group had higher radial indices and lower laminar and longitudinal indices in both the humerus and the femur than the restricted group. The higher radial indices in our control group point to a link between radial canals and faster growth, and between laminar canals and slower growth, while the higher laminar indices in the humerus point to a link between circumferential canals and torsional loading. Overall, our results indicate that the orientation of the cortical canal network in a bone is the consequence of a complex interaction between the growth rate of that bone and functional loading environment.

Construction of poly-ellipsoidal grain shapes from SMT imaging on sand, and the development of a new DEM contact detection algorithm

PurposeThis paper aims to construct smooth poly-ellipsoid shapes from synchrotron microcomputed tomography (SMT) images on sand and to develop a new discrete element method (DEM) contact detection algorithm.Design/methodology/approachVoxelated images generated by SMT on Colorado Mason sand are processed to construct smooth poly-ellipsoidal particle approximations. For DEM contact detection, cuboidal shape approximations to the poly-ellipsoids are used to speed up contact detection.FindingsThe poly-ellipsoid particle shape approximation to Colorado Mason sand grains is better than a simpler ellipsoidal approximation. The new DEM contact algorithm leads to significant speedup and accuracy is maintained.Research limitations/implicationsThe paper limits particle shape approximation to smooth poly-ellipsoids.Practical implicationsPoly-ellipsoids provide asymmetry of particle shapes as compared to ellipsoids, thus allowing closer representation of real sand grain shapes that may be angular and unsymmetric. When incorporated in a DEM for computation, the poly-ellipsoids allow better representation of particle rolling, sliding and interlocking phenomena.Originality/valueMethod to construct poly-ellipsoid particle shapes from SMT data on real sands and computationally efficient DEM contact detection algorithm for poly-ellipsoids.

3D experimental quantification of fabric and fabric evolution of sheared granular materials using synchrotron micro-computed tomography

Many experimental studies have demonstrated that mechanical response of granular materials is highly influenced by micro-structural fabric and its evolution. In the current literature, quantification of fabric and its evolution has been developed based on micro-structural observations using Discrete Element Method or 2D experiments with simple particle shapes. The emergence of X-ray computed tomography technique has made quantification of such experimental micro-structural properties possible using 3D high-resolution images. In this paper, synchrotron micro-computed tomography was used to acquire 3D images during in-situ conventional triaxial compression experiments on granular materials with different morphologies. 3D images were processed to quantify fabric and its evolution based on experimental measurements of contact normal vectors between particles. Overall, the directional distribution of contact normals exhibited the highest degree of isotropy at initial state (i.e., zero global axial strain). As compression progressed, contact normals evolved in the direction of loading until reaching a constant fabric when experiments approached the critical state condition. Further assessment of the influence of confining pressure, initial density state, and particle-level morphology on fabric and its evolution was formed. Results show that initial density state and applied confining pressure significantly influence the fabric-induced internal anisotropy of tested specimens at initial states. Relatively, a higher applied confining pressure and a looser initial density state resulted in a higher degree of fabric-induced internal anisotropy. Influence of particle-level morphology was also found to be significant particularly on fabric evolution.

Endothiodon cf. bathystoma (Synapsida: Dicynodontia) bony labyrinth anatomy, variation and body mass estimates.

The semicircular canal (SC) system of the inner ear detects head angular accelerations and is essential for navigation and spatial awareness in vertebrates. Because the bony labyrinth encloses the membranous labyrinth SCs, it can be used as a proxy for animal behavior. The bony labyrinth of dicynodonts, a clade of herbivorous non-mammalian synapsids, has only been described in a handful of individuals and remains particularly obscure. Here we describe the bony labyrinth anatomy of three Endothiodon cf. bathystoma specimens from Mozambique based on digital reconstructions from propagation phase-contrast synchrotron micro-computed tomography. We compare these findings with the bony labyrinth anatomy of their close relative Niassodon. The bony labyrinths of Endothiodon and Niassodon are relatively similar and show only differences in the shape of the horizontal SCs and the orientation of the vertical SCs. When compared to extant mammals, Endothiodon and Niassodon have highly eccentric SCs. In addition, the Endothiodon SCs are nearly orthogonal. An eccentric and orthogonal SC morphology is consistent with a specialization in rapid head movements, which are typical of foraging or feeding behaviors. Furthermore, we estimate the body mass of these Endothiodon specimens at ~116 to 182 kg, based on the average SC radii calculated using a linear regression model optimized by the Amemiya Prediction Criterion. Our findings provide novel insights into the paleobiology of Endothiodon which are consistent with the peculiar feeding mechanism among dicynodonts presumed from their multiple postcanine toothrows.

Interpreting the three-dimensional orientation of vascular canals and cross-sectional geometry of cortical bone in birds and bats.

Cortical bone porosity and specifically the orientation of vascular canals is an area of growing interest in biomedical research and comparative/paleontological anatomy. The potential to explain microstructural adaptation is of great interest. However, the determinants of the development of canal orientation remain unclear. Previous studies of birds have shown higher proportions of circumferential canals (called laminarity) in flight bones than in hindlimb bones, and interpreted this as a sign that circumferential canals are a feature for resistance to the torsional loading created by flight. We defined the laminarity index as the percentage of circumferential canal length out of the total canal length. In this study we examined the vascular canal network in the humerus and femur of a sample of 31 bird and 24 bat species using synchrotron micro‐computed tomography (micro‐CT) to look for a connection between canal orientation and functional loading. The use of micro‐CT provides a full three‐dimensional (3D) map of the vascular canal network and provides measurements of the 3D orientation of each canal in the whole cross‐section of the bone cortex. We measured several cross‐sectional geometric parameters and strength indices including principal and polar area moments of inertia, principal and polar section moduli, circularity, buckling ratio, and a weighted cortical thickness index. We found that bat cortices are relatively thicker and poorly vascularized, whereas those of birds are thinner and more highly vascularized, and that according to our cross‐sectional geometric parameters, bird bones have a greater resistance to torsional stress than the bats; in particular, the humerus in birds is more adapted to resist torsional stresses than the femur. Our results show that birds have a significantly (P = 0.031) higher laminarity index than bats, with birds having a mean laminarity index of 0.183 in the humerus and 0.232 in the femur, and bats having a mean laminarity index of 0.118 in the humerus and 0.119 in the femur. Counter to our expectation, the birds had a significantly higher laminarity index in the femur than in the humerus (P = 0.035). To evaluate whether this discrepancy was a consequence of methodology we conducted a comparison between our 3D method and an analogue to two‐dimensional (2D) histological measurements. This comparison revealed that 2D methods significantly underestimate (P < 0.001) the amount of longitudinal canals by an average of 20% and significantly overestimate (P < 0.001) the laminarity index by an average of 7.7%, systematically mis‐estimating indices of vascular canal orientations. In comparison with our 3D results, our approximated 2D measurement had the same results for comparisons between the birds and bats but found significant differences only in the longitudinal index between the humerus and the femur for both groups. The differences between our 3D and pseudo‐2D results indicate that differences between our findings and the literature may be partially based in methodology. Overall, our results do not support the hypothesis that the bones of flight are more laminar, suggesting a complex relation between functional loading and microstructural adaptation.

Experimental Fracture Analysis of Individual Sand Particles at High Loading Rates

AbstractHigh strain rate loading conditions such as blast, impact, or projectile penetration cause major damages to infrastructure and soils. Sand particles will likely fracture if they are exposed to such loading conditions. Modeling the constitutive behavior of sand deposits when they are subjected to dynamic loading requires a high-fidelity particle-scale analysis, calibration, and validation using experimental measurements. In this paper, natural sand particles with complex morphology, mineralogy, and grain size were randomly selected from the bulk material and were separated based on particle size and mineralogy. Intact particles were first imaged using a desktop X-ray computed tomography (CT) scanner to calculate morphology and provide further input to the mineralogy and internal structure. Each particle was placed between incident and transmitter bars of a Kolsky bar setup, and two optical images (side and top views) were collected to visualize the loading geometry and direction. Each particle was then dynamically loaded to failure, at which load and compressive displacement were measured. Some of recovered particle fragments were imaged using synchrotron microcomputed tomography (SMT) to determine the fracture mode and fracture surface. Weibull statistical analyses were performed and multivariable nonlinear regression was implemented, using particle characteristics as predictors. The paper discusses the influence of morphology, mineralogy, internal structure, and size of particles on dynamic particle failure strength.

3D finite element modelling of force transmission and particle fracture of sand

Global compressive loading of granular media causes rearrangements of particles into a denser configuration. Under 1D compression, researchers observed that particles initially translate and rotate which lead to more contacts between particles and the development of force chains to resist applied loads. Particles within force chains resist most of the applied loads while neighbor particles provide lateral support to prevent particles within force chains from buckling. Several experimental and numerical models have been proposed in the literature to characterize force chains within granular materials. This paper presents a 3D finite element (FE) model that simulates 1D compression experiment on F-75 Ottawa sand. The FE mesh of particles closely matched 3D physical shape of sand particles that were acquired using 3D synchrotron micro-computed tomography (SMT) technique. The paper presents a quantitative assessment of the model, in which evolution of force chains, fracture modes, and stress-strain relationships showed an excellent agreement with experimental measurements reported by Cil et al. Alshibli (2017).

3D Experimental Measurement of Lattice Strain and Fracture Behavior of Sand Particles Using Synchrotron X-Ray Diffraction and Tomography

AbstractThree-dimensional synchrotron X-ray diffraction (3DXRD) and synchrotron microcomputed tomography (SMT) techniques were used to measure and monitor the lattice strain evolution and fracture ...

Synchrotron Micro-computed Tomography Reveals Increased Osteocyte Volume in Subchondral Bone with Osteoarthritis Severity

Synchrotron Micro-computed Tomography Reveals Increased Osteocyte Volume in Subchondral Bone with Osteoarthritis Severity

Influence of Particle Morphology on 3D Kinematic Behavior and Strain Localization of Sheared Sand

AbstractThe constitutive behavior of sheared sand is highly influenced by particle morphology, gradation, mineralogy, specimen density, loading condition, stress path, and boundary conditions. The current literature lacks a three-dimensional (3D) systematic experimental study that investigates the influence of particle morphology, confining pressure, and specimen density on the failure mode of sheared sand. In this paper, surface texture, roundness, and sphericity of three uniform sands and glass beads with similar grain size were quantified by using 3D images of particles. In situ nondestructive 3D synchrotron microcomputed tomography (SMT) was used to monitor the deformation of medium-dense and very dense dry sand specimens that were tested under axisymmetric triaxial loading condition at 15 and 400 kPa confining pressures. The particles were identified and tracked in 3D as shearing progressed within the specimens, and maps of incremental particle translation and rotation were developed and used to unco...

Multiscale alterations in bone matrix quality increased fragility in steroid induced osteoporosis

A serious adverse clinical effect of glucocorticoid steroid treatment is secondary osteoporosis, enhancing fracture risk in bone. This rapid increase in bone fracture risk is largely independent of bone loss (quantity), and must therefore arise from degradation of the quality of the bone matrix at the micro- and nanoscale. However, we lack an understanding of both the specific alterations in bone quality n steroid-induced osteoporosis as well as the mechanistic effects of these changes. Here we demonstrate alterations in the nanostructural parameters of the mineralized fibrillar collagen matrix, which affect bone quality, and develop a model linking these to increased fracture risk in glucocorticoid induced osteoporosis. Using a mouse model with an N-ethyl-N-nitrosourea (ENU)-induced corticotrophin releasing hormone promoter mutation (Crh−120/+) that developed hypercorticosteronaemia and osteoporosis, we utilized in situ mechanical testing with small angle X-ray diffraction, synchrotron micro-computed tomography and quantitative backscattered electron imaging to link altered nano- and microscale deformation mechanisms in the bone matrix to abnormal macroscopic mechanics. We measure the deformation of the mineralized collagen fibrils, and the nano-mechanical parameters including effective fibril modulus and fibril to tissue strain ratio. A significant reduction (51%) of fibril modulus was found in Crh−120/+ mice. We also find a much larger fibril strain/tissue strain ratio in Crh−120/+ mice (~1.5) compared to the wild-type mice (~0.5), indicative of a lowered mechanical competence at the nanoscale. Synchrotron microCT show a disruption of intracortical architecture, possibly linked to osteocytic osteolysis. These findings provide a clear quantitative demonstration of how bone quality changes increase macroscopic fragility in secondary osteoporosis.

Failure and failure mechanisms of wood during longitudinal compression monitored by synchrotron micro-computed tomography

Abstract The lasting plastic deformation of the cellular elements of beech, fir, and spruce wood under uniaxial compression has been investigated by the combination of in situ loading and synchrotron micro-computed tomography. The deformation of singular elements embedded in the tissue and their influence on the deformation lines and surrounding tissue was examined by tomographic reconstructions. An automatic observation of the failure area was applied in the case of the softwoods, which permitted the determination of the densification degree. The development of failure lines differed from the expected pattern by often showing branching. However, the 3D observations confirmed the effects known from 2D examinations in many cases.