Effect of beam hardening on arterial enhancement in thoracoabdominal CT angiography with increasing patient size: an in vitro and in vivo study

To assess the effect of beam hardening on arterial enhancement in thoracoabdominal computed tomographic (CT) angiography in various body sizes in a phantom and in a clinical study.

Monoenergetic computed tomography reconstructions reduce beam hardening artifacts from dental restorations

The aim of this study was to assess the potential of monoenergetic computed tomography (CT) images to reduce beam hardening artifacts in comparison to standard CT images of dental restoration on dental post-mortem CT (PMCT). Thirty human decedents (15 male, 58 ± 22 years) with dental restorations were examined using standard single-energy CT (SECT) and dual-energy CT (DECT). DECT data were used to generate monoenergetic CT images, reflecting the X-ray attenuation at energy levels of 64, 69, 88 keV, and at an individually adjusted optimal energy level called OPTkeV. Artifact reduction and image quality of SECT and monoenergetic CT were assessed objectively and subjectively by two blinded readers. Subjectively, beam artifacts decreased visibly in 28/30 cases after monoenergetic CT reconstruction. Inter- and intra-reader agreement was good (k = 0.72, and k = 0.73 respectively). Beam hardening artifacts decreased significantly with increasing monoenergies (repeated-measures ANOVA p < 0.001). Artifact reduction was greatest on monoenergetic CT images at OPTkeV. Mean OPTkeV was 108 ± 17 keV. OPTkeV yielded the lowest difference between CT numbers of streak artifacts and reference tissues (-163 HU). Monoenergetic CT reconstructions significantly reduce beam hardening artifacts from dental restorations and improve image quality of post-mortem dental CT.

Variation of the mechanical properties of PMMA to suit osteoporotic cancellous bone

Poly(methyl methacrylate) (PMMA) is by far the most frequently used bone substitute material for vertebroplasty. However, there are serious complications, such as cement leakage and an increased fracture rate of the adjacent vertebral bodies. The latter may be related to the mechanical properties of the augmented segment within the osteoporotic spine. A possible counter-measure is prophylactic augmentation at additional levels, but this aggravates the risk for the patient. Introduction of pores is a possible method to reduce the inherent high stiffness of PMMA. This study investigates the effect of porosity on the mechanical properties of PMMA bone cement. Different fractions of a highly viscous liquid were mixed into the PMMA during preparation. An open-porous material with adjustable mechanical properties resulted after removal of the aqueous phase. Different radiopacifiers were admixed to investigate their suitability for vertebroplasty. The final material was characterized mechanically by compressive testing, microscopically and radiologically. In addition, the monomer release subsequent to hardening was measured by means of gas chromatography. The Young's modulus in compression could be varied between 2800 +/- 70 MPa and 120 +/- 150 MPa, and the compression ultimate strength between 170 +/- 5 MPa and 8 +/- 9 MPa for aqueous fractions ranging between 0 and 50% of volume. Only a slight decrease of the Young's modulus and small changes of ultimate strength were found when the mixing time was increased. An organic hydrophilic and lipophilic radiopacifier led to a higher Young's modulus of the porous material; however, the ultimate strength was not significantly affected by adding different radiopacifiers to the porous cement. The radiopacity was lost after washing the aqueous phase out of the pores. No separation occurred between the aqueous and the PMMA phase during injection into an open porous ceramic material. The monomer released was found to increase for increasing aqueous fractions, but remained comparable in magnitude to standard PMMA. This study demonstrates that a conventional PMMA can be modified to obtain a range of mechanical properties, including those of osteoporotic bone.

Clinical investigations of polymethylmethacrylate cement viscosity during vertebroplasty and related in vitro measurements

Percutaneous vertebroplasty, comprising an injection of polymethylmethacrylate (PMMA) into vertebral bodies, is a practical procedure for the stabilization of osteoporotic compression fractures as well as other weakening lesions. Cement leakage is considered to be one of the major and most severe complications during percutaneous vertebroplasty. The viscosity of the material plays a key role in this context. In order to enhance the safety for the patient, a rheometer system was developed to measure the cement viscosity intraoperatively. For this development, it is of great importance to know the proper viscosity to start the procedure determined by experienced surgeons and the relation between the time period when different injection devices are used and the cement viscosity. The purpose of the study was to investigate the viscosity ranges for different injection systems during conventional vertebroplasty. Clinically observed viscosity values and related time periods showed high scattering. In order to get a better understanding of the clinical observations, cement viscosity during hardening at different ambient temperatures and by simulation of the body temperature was investigated in vitro. It could be concluded, that the direct viscosity assessment with a rheometer during vertebroplasty can help clinicians to define a lower threshold viscosity and thereby decrease the risk of leakage and make adjustments to their injection technique in real time. Secondly, the acceleration in hardening of PMMA-based cements at body temperature can be useful in minimizing leakages by addressing them with a short injection break.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

From transient to steady state deformation and grain size: A thermodynamic approach using elasto-visco-plastic numerical modeling

Numerical simulation experiments give insight into the evolving energy partitioning during high-strain torsion experiments of calcite. Our numerical experiments are designed to derive a generic macroscopic grain size sensitive flow law capable of describing the full evolution from the transient regime to steady state. The transient regime is crucial for understanding the importance of micro structural processes that may lead to strain localization phenomena in deforming materials. This is particularly important in geological and geodynamic applications where the phenomenon of strain localization happens outside the time frame that can be observed under controlled laboratory conditions. Ourmethod is based on an extension of the paleowattmeter approach to the transient regime. We add an empirical hardening law using the Ramberg-Osgood approximation and assess the experiments by an evolution test function of stored over dissipated energy (lambda factor). Parameter studies of, strain hardening, dislocation creep parameter, strain rates, temperature, and lambda factor as well asmesh sensitivity are presented to explore the sensitivity of the newly derived transient/steady state flow law. Our analysis can be seen as one of the first steps in a hybrid computational-laboratory-field modeling workflow. The analysis could be improved through independent verifications by thermographic analysis in physical laboratory experiments to independently assess lambda factor evolution under laboratory conditions.

Experimental, theoretical and numerical investigation of the nonlinear micromechanical properties of bone

Aging societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the bone hierarchical organization. A good understanding has been reached for elastic properties on several length scales, but up to now there is a lack of reliable postyield data on the lower length scales. In order to be able to describe the behavior of bone at the microscale, an anisotropic elastic-viscoplastic damage model was developed using an eccentric generalized Hill criterion and nonlinear isotropic hardening. The model was implemented as a user subroutine in Abaqus and verified using single element tests. A FE simulation of microindentation in lamellar bone was finally performed show-ing that the new constitutive model can capture the main characteristics of the indentation response of bone. As the generalized Hill criterion is limited to elliptical and cylindrical yield surfaces and the correct shape for bone is not known, a new yield surface was developed that takes any convex quadratic shape. The main advantage is that in the case of material identification the shape of the yield surface does not have to be anticipated but a minimization results in the optimal shape among all convex quadrics. The generality of the formulation was demonstrated by showing its degeneration to classical yield surfaces. Also, existing yield criteria for bone at multiple length scales were converted to the quadric formulation. Then, a computational study to determine the influence of yield surface shape and damage on the in-dentation response of bone using spherical and conical tips was performed. The constitutive model was adapted to the quadric criterion and yield surface shape and critical damage were varied. They were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic to total work ratio were found to be very well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not a significant fac-tor, while for spherical tips damage was insignificant. All inverse methods based on microindentation suffer from a lack of uniqueness of the found material properties in the case of nonlinear material behavior. Therefore, monotonic and cyclic micropillar com-pression tests in a scanning electron microscope allowing a straightforward interpretation comple-mented by microindentation and macroscopic uniaxial compression tests were performed on dry ovine bone to identify modulus, yield stress, plastic deformation, damage accumulation and failure mecha-nisms. While the elastic properties were highly consistent, the postyield deformation and failure mech-anisms differed between the two length scales. A majority of the micropillars showed a ductile behavior with strain hardening until failure by localization in a slip plane, while the macroscopic samples failed in a quasi-brittle fashion with microcracks coalescing into macroscopic failure surfaces. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behavior of bone at the microscale to a quasi-brittle response driven by the growth of preexisting cracks along interfaces or in the vicinity of pores at the macroscale. Subsequently, a study was undertaken to quantify the topological variability of indentations in bone and examine its relationship with mechanical properties. Indentations were performed in dry human and ovine bone in axial and transverse directions and their topography measured by AFM. Statistical shape modeling of the residual imprint allowed to define a mean shape and describe the variability with 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was highly consistent and free of any pile up. A few of the topological parameters, in particular depth, showed significant correlations to variations in mechanical properties, but the cor-relations were not very strong or consistent. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small de-viations from the ideal case. As the uniaxial properties measured by micropillar compression are in conflict with the current literature on bone indentation, another dissipative mechanism has to be present. The elastic-viscoplastic damage model was therefore extended to viscoelasticity. The viscoelastic properties were identified from macroscopic experiments, while the quasistatic postelastic properties were extracted from micropillar data. It was found that viscoelasticity governed by macroscale properties has very little influence on the indentation curve and results in a clear underestimation of the creep deformation. Adding viscoplasticity leads to increased creep, but hardness is still highly overestimated. It was possible to obtain a reasonable fit with experimental indentation curves for both Berkovich and spherical indenta-tion when abandoning the assumption of shear strength being governed by an isotropy condition. These results remain to be verified by independent tests probing the micromechanical strength prop-erties in tension and shear. In conclusion, in this thesis several tools were developed to describe the complex behavior of bone on the microscale and experiments were performed to identify its material properties. Micropillar com-pression highlighted a size effect in bone due to the presence of preexisting cracks and pores or inter-faces like cement lines. It was possible to get a reasonable fit between experimental indentation curves using different tips and simulations using the constitutive model and uniaxial properties measured by micropillar compression. Additional experimental work is necessary to identify the exact nature of the size effect and the mechanical role of interfaces in bone. Deciphering the micromechanical behavior of lamellar bone and its evolution with age, disease and treatment and its failure mechanisms on several length scales will help preventing fractures in the elderly in the future.

Interpretation of fluid inclusions in quartz deformed by weak ductile shearing: reconstruction of differential stress magnitudes and pre-deformation fluid properties

A well developed theoretical framework is available in which paleofluid properties, such as chemical composition and density, can be reconstructed from fluid inclusions in minerals that have undergone no ductile deformation. The present study extends this framework to encompass fluid inclusions hosted by quartz that has undergone weak ductile deformation following fluid entrapment. Recent experiments have shown that such deformation causes inclusions to become dismembered into clusters of irregularly shaped relict inclusions surrounded by planar arrays of tiny, new-formed (neonate) inclusions. Comparison of the experimental samples with a naturally sheared quartz vein from Grimsel Pass, Aar Massif, Central Alps, Switzerland, reveals striking similarities. This strong concordance justifies applying the experimentally derived rules of fluid inclusion behaviour to nature. Thus, planar arrays of dismembered inclusions defining cleavage planes in quartz may be taken as diagnostic of small amounts of intracrystalline strain. Deformed inclusions preserve their pre-deformation concentration ratios of gases to electrolytes, but their H2O contents typically have changed. Morphologically intact inclusions, in contrast, preserve the pre-deformation composition and density of their originally trapped fluid. The orientation of the maximum principal compressive stress (σ1σ1) at the time of shear deformation can be derived from the pole to the cleavage plane within which the dismembered inclusions are aligned. Finally, the density of neonate inclusions is commensurate with the pressure value of σ1σ1 at the temperature and time of deformation. This last rule offers a means to estimate magnitudes of shear stresses from fluid inclusion studies. Application of this new paleopiezometer approach to the Grimsel vein yields a differential stress (σ1–σ3σ1–σ3) of ∼300 MPa∼300 MPa at View the MathML source390±30°C during late Miocene NNW–SSE orogenic shortening and regional uplift of the Aar Massif. This differential stress resulted in strain-hardening of the quartz at very low total strain (<5%<5%) while nearby shear zones were accommodating significant displacements. Further implementation of these experimentally derived rules should provide new insight into processes of fluid–rock interaction in the ductile regime within the Earth's crust.

Continuum damage interactions between tension and compression in osteonal bone

Skeletal diseases such as osteoporosis impose a severe socio-economic burden to ageing societies. Decreasing mechanical competence causes a rise in bone fracture incidence and mortality especially after the age of 65 y. The mechanisms of how bone damage is accumulated under different loading modes and its impact on bone strength are unclear. We hypothesise that damage accumulated in one loading mode increases the fracture risk in another. This study aimed at identifying continuum damage interactions between tensile and compressive loading modes. We propose and identify the material constants of a novel piecewise 1D constitutive model capable of describing the mechanical response of bone in combined tensile and compressive loading histories. We performed several sets of loading–reloading experiments to compute stiffness, plastic strains, and stress-strain curves. For tensile overloading, a stiffness reduction (damage) of 60% at 0.65% accumulated plastic strain was detectable as stiffness reduction of 20% under compression. For compressive overloading, 60% damage at 0.75% plastic strain was detectable as a stiffness reduction of 50% in tension. Plastic strain at ultimate stress was the same in tension and compression. Compression showed softening and tension exponential hardening in the post-yield regime. The hardening behaviour in compression is unaffected by a previous overload in tension but the hardening behaviour in tension is affected by a previous overload in compression as tensile reloading strength is significantly reduced. This paper demonstrates how damage accumulated under one loading mode affects the mechanical behaviour in another loading mode. To explain this and to illustrate a possible implementation we proposed a theoretical model. Including such loading mode dependent damage and plasticity behaviour in finite element models will help to improve fracture risk analysis of whole bones and bone implant structures.

Properties of granular analogue model materials: A community wide survey

We report the material properties of 26 granular analogue materials used in 14 analogue modelling laboratories. We determined physical characteristics such as bulk density, grain size distribution, and grain shape, and performed ring shear tests to determine friction angles and cohesion, and uniaxial compression tests to evaluate the compaction behaviour. Mean grain size of the materials varied between c. 100 and 400 μm. Analysis of grain shape factors shows that the four different classes of granular materials (14 quartz sands, 5 dyed quartz sands, 4 heavy mineral sands and 3 size fractions of glass beads) can be broadly divided into two groups consisting of 12 angular and 14 rounded materials. Grain shape has an influence on friction angles, with most angular materials having higher internal friction angles (between c. 35° and 40°) than rounded materials, whereas well-rounded glass beads have the lowest internal friction angles (between c. 25° and 30°). We interpret this as an effect of intergranular sliding versus rolling. Most angular materials have also higher basal friction angles (tested for a specific foil) than more rounded materials, suggesting that angular grains scratch and wear the foil. Most materials have an internal cohesion in the order of 20–100 Pa except for well-rounded glass beads, which show a trend towards a quasi-cohesionless (C < 20 Pa) Coulomb-type material. The uniaxial confined compression tests reveal that rounded grains generally show less compaction than angular grains. We interpret this to be related to the initial packing density after sifting, which is higher for rounded grains than for angular grains. Ring-shear test data show that angular grains undergo a longer strain-hardening phase than more rounded materials. This might explain why analogue models consisting of angular grains accommodate deformation in a more distributed manner prior to strain localisation than models consisting of rounded grains.