Specimen specific parameter identification of ovine lumbar intervertebral discs: On the influence of fibre-matrix and fibre-fibre shear interactions.

Numerical models of the intervertebral disc, which address mechanical questions commonly make use of the difference in water content between annulus and nucleus, and thus fluid and solid parts are separated. Despite this simplification, models remain complex due to the anisotropy and nonlinearity of the annulus and regional variations of the collagen fibre density. Additionally, it has been shown that cross-links make a large contribution to the stiffness of the annulus. Because of this complex composite structure, it is difficult to reproduce several sets of experimental data with one single set of material parameters. This study addresses the question to which extent the ultrastructure of the intervertebral disc should be modelled so that its moment-angle behaviour can be adequately described. Therefore, a hyperelastic constitutive law, based on continuum mechanical principles was derived, which does not only consider the anisotropy from the collagen fibres, but also interactions among the fibres and between the fibres and the ground substance. Eight ovine lumbar intervertebral discs were tested on a custom made spinal loading simulator in flexion/extension, lateral bending and axial rotation. Specimen-specific geometrical models were generated using CT images and T2 maps to distinguish between annulus fibrosus and nucleus pulposus. For the identification of the material parameters the annulus fibrosus was described with two scenarios: with and without fibre-matrix and fibre-fibre interactions. Both scenarios showed a similar behaviour on a load displacement level. Comparing model predictions to the experimental data, the mean RMS of all specimens and all load cases was 0.54±0.15° without the interaction and 0.54±0.19° when the fibre-matrix and fibre-fibre interactions were included. However, due to the increased stiffness when cross-links effects were included, this scenario showed more physiological stress-strain relations in uniaxial and biaxial stress states. Thus, the present study suggests that fibre-matrix and fibre-fibre interactions should be considered in the constitutive law when the model addresses questions concerning the stress field of the annulus fibrosus.

Experimental validation of a nonlinear μ FE model based on cohesive-frictional plasticity for trabecular bone

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of μFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

Functionalization of deproteinized bovine bone with a coating-incorporated depot of BMP-2 renders the material efficiently osteoinductive and suppresses foreign-body reactivity

The repair of critical-sized bony defects remains a challenge in the fields of implantology, maxillofacial surgery and orthopaedics. As an alternative bone-defect filler to autologous bone grafts, deproteinized bovine bone (DBB) is highly osteoconductive and clinically now widely used. However, this product suffers from the disadvantage of not being intrinsically osteoinductive. In the present study, this property was conferred by coating DBB with a layer of calcium phosphate into which bone morphogenetic protein 2 (BMP-2) was incorporated. Granules of DBB bearing a coating-incorporated depot of BMP-2--together with the appropriate controls (DBB bearing a coating but no BMP-2; uncoated DBB bearing adsorbed BMP-2; uncoated DBB bearing no BMP-2)--were implanted subcutaneously in rats. Five weeks later, the implants were withdrawn for a histomorphometric analysis of the volume densities of (i) bone, (ii) bone marrow, (iii) foreign-body giant cells and (iv) fibrous capsular tissue. Parameters (i) and (ii) were highest, whilst parameters (iii) and (iv) were lowest in association with DBB bearing a coating-incorporated depot of BMP-2. Hence, this mode of functionalization not only confers DBB with the property of osteoinductivity but also improves its biocompatibility--thus dually enhancing its clinical potential in the repair of bony defects.

Conversion of CT numbers into tissue parameters for Monte Carlo dose calculations: a multi-centre study

The conversion of computed tomography (CT) numbers into material composition and mass density data influences the accuracy of patient dose calculations in Monte Carlo treatment planning (MCTP). The aim of our work was to develop a CT conversion scheme by performing a stoichiometric CT calibration. Fourteen dosimetrically equivalent tissue subsets (bins), of which ten bone bins, were created. After validating the proposed CT conversion scheme on phantoms, it was compared to a conventional five bin scheme with only one bone bin. This resulted in dose distributions D(14) and D(5) for nine clinical patient cases in a European multi-centre study. The observed local relative differences in dose to medium were mostly smaller than 5%. The dose-volume histograms of both targets and organs at risk were comparable, although within bony structures D(14) was found to be slightly but systematically higher than D(5). Converting dose to medium to dose to water (D(14) to D(14wat) and D(5) to D(5wat)) resulted in larger local differences as D(5wat) became up to 10% higher than D(14wat). In conclusion, multiple bone bins need to be introduced when Monte Carlo (MC) calculations of patient dose distributions are converted to dose to water.

Maxillary sinus floor elevation using the (transalveolar) osteotome technique with or without grafting material. Part I: Implant survival and patients' perception

OBJECTIVES: To analyze the survival and success rates of implants installed utilizing the (transalveolar) osteotome technique, to compare peri-implant soft tissue parameters and marginal bone levels of osteotome-installed implants with implants placed using standard surgical procedures, and to evaluate patient-centered outcomes. MATERIAL AND METHODS: During 2000 to 2005, 252 Straumann dental implants were inserted in 181 patients. The surgical technique was a modification of the original osteotome technique presented by Summers. In addition to the clinical examination, the patients were asked to give their perception of the surgical procedure, utilizing a visual analogue scale. RESULTS: The cumulative survival rate of the osteotome-installed implants after a mean follow-up time of 3.2 years, was 97.4% (95% confidence intervals: 94.4-98.8%). From the 252 implants inserted, three were lost before loading and another three were lost in the first and second year. According to residual bone height the survival was 91.3% for implant sites with < or =4 mm residual bone height, and 90% for sites with 4 mm and 5 mm, when compared with that of 100% in sites with bone height of above 5 mm. According to implant length the survival rates were 100% for 12 mm, 98.7% for 10 mm, 98.7% for 8 mm and only 47.6% for 6 mm implants. Soft tissue parameters (pocket probing depth, probing attachment level, bleeding on probing and marginal bone levels) did not yield any differences between the osteotome-installed and the conventionally placed implants. More than 90% of the patients were satisfied with the implant therapy and would undergo similar therapy again if necessary. The cost associated with implant therapy was considered to be justified. CONCLUSION: In conclusion, the osteotome technique was a reliable method for implant insertion in the posterior maxilla, especially at sites with 5 mm or more of preoperative residual bone height and a relatively flat sinus floor.

High-Pressure Behavior and Phase Stability of Al5BO9, a Mullite-Type Ceramic Material

Phase stability, elastic behavior, and pressure-induced structural evolution of synthetic boron-mullite Al5BO9 (a = 5.6780(7), b = 15.035(6), and c =7.698(3) Å, space group Cmc21, Z = 4) were investigated up to 25.6(1) GPa by in situ single-crystal synchrotron X-ray diffraction with a diamond anvil cell (DAC) under hydrostatic conditions. No evidence of phase transition was observed up to 21.7(1) GPa. At 25.6(1) GPa, the refined unit-cell parameters deviated significantly from the compressional trend, and the diffraction peaks appeared broader than at lower pressure. At 26.7(1) GPa, the diffraction pattern was not indexable, suggesting amorphization of the material or a phase transition to a high-pressure polymorph. Fitting the P–V data up to 21.7(1) GPa with a second-order Birch–Murnaghan Equation-of-State, we obtained a bulk modulus KT0 = 164(1) GPa. The axial compressibilities, here described as linearized bulk moduli, are as follows: KT0(a) = 244(9), KT0(b) = 120(4), and KT0(c) = 166(11) GPa (KT0(a):KT0(b):KT0(c) = 2.03:1:1.38). The structure refinements allowed a description of the main deformation mechanisms in response to the applied pressure. The stiffer crystallographic direction appears to be controlled by the infinite chains of edge-sharing octahedra running along [100], making the structure less compressible along the a-axis than along the b- and c-axis.

Microstructural Parameters of Bone Evaluated Using HR-pQCT Correlate with the DXA-Derived Cortical Index and the Trabecular Bone Score in a Cohort of Randomly Selected Premenopausal Women.

BACKGROUND Areal bone mineral density is predictive for fracture risk. Microstructural bone parameters evaluated at the appendicular skeleton by high-resolution peripheral quantitative computed tomography (HR-pQCT) display differences between healthy patients and fracture patients. With the simple geometry of the cortex at the distal tibial diaphysis, a cortical index of the tibia combining material and mechanical properties correlated highly with bone strength ex vivo. The trabecular bone score derived from the scan of the lumbar spine by dual-energy X-ray absorptiometry (DXA) correlated ex vivo with the micro architectural parameters. It is unknown if these microstructural correlations could be made in healthy premenopausal women. METHODS Randomly selected women between 20-40 years of age were examined by DXA and HR-pQCT at the standard regions of interest and at customized sub regions to focus on cortical and trabecular parameters of strength separately. For cortical strength, at the distal tibia the volumetric cortical index was calculated directly from HR-pQCT and the areal cortical index was derived from the DXA scan using a Canny threshold-based tool. For trabecular strength, the trabecular bone score was calculated based on the DXA scan of the lumbar spine and was compared with the corresponding parameters derived from the HR-pQCT measurements at radius and tibia. RESULTS Seventy-two healthy women were included (average age 33.8 years, average BMI 23.2 kg/m(2)). The areal cortical index correlated highly with the volumetric cortical index at the distal tibia (R  =  0.798). The trabecular bone score correlated moderately with the microstructural parameters of the trabecular bone. CONCLUSION This study in randomly selected premenopausal women demonstrated that microstructural parameters of the bone evaluated by HR-pQCT correlated with the DXA derived parameters of skeletal regions containing predominantly cortical or cancellous bone. Whether these indexes are suitable for better predictions of the fracture risk deserves further investigation.

Provenance analysis using Raman spectroscopy of carbonaceous material: A case study in the Southern Alps of New Zealand

Detrital provenance analyses in orogenic settings, in which sediments are collected at the outlet of a catchment, have become an important tool to estimate how erosion varies in space and time. Here we present how Raman Spectroscopy on Carbonaceous Material (RSCM) can be used for provenance analysis. RSCM provides an estimate of the peak temperature (RSCM-T) experienced during metamorphism. We show that we can infer modern erosion patterns in a catchment by combining new measurements on detrital sands with previously acquired bedrock data. We focus on the Whataroa catchment in the Southern Alps of New Zealand and exploit the metamorphic gradient that runs parallel to the main drainage direction. To account for potential sampling biases, we also quantify abrasion properties using flume experiments and measure the total organic carbon content in the bedrock that produced the collected sands. Finally, we integrate these parameters into a mass-conservative model. Our results first demonstrate that RSCM-T can be used for detrital studies. Second, we find that spatial variations in tracer concentration and erosion have a first-order control on the RSCM-T distributions, even though our flume experiments reveal that weak lithologies produce substantially more fine particles than do more durable lithologies. This result implies that sand specimens are good proxies for mapping spatial variations in erosion when the bedrock concentration of the target mineral is quantified. The modeling suggests that highest present-day erosion rates (in Whataroa catchment) are not situated at the range front but around 10 km into the mountain belt.