In situ micropillar compression reveals superior strength and ductility but an absence of damage in lamellar bone

Ageing 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 hierarchical organization of bone. Here, we investigate the mechanical response under monotonic and cyclic compression of both single osteonal lamellae and macroscopic samples containing numerous osteons. Micropillar compression tests in a scanning electron microscope, microindentation and macroscopic compression tests were performed on dry ovine bone to identify the elastic modulus, yield stress, plastic deformation, damage accumulation and failure mechanisms. We found that isolated lamellae exhibit a plastic behaviour, with higher yield stress and ductility but no damage. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behaviour of bone at the microscale to a quasi-brittle response driven by the growth of cracks along interfaces or in the vicinity of pores at the macroscale.

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.

Determining gypsum growth temperatures using monophase fluid inclusions-Application to the giant gypsum crystals of Naica, Mexico

Determining the formation temperature of minerals using fluid inclusions is a crucial step in understanding rock-forming scenarios. Unfortunately, fluid inclusions in minerals formed at low temperature, such as gypsum, are commonly in a metastable monophase liquid state. To overcome this problem, ultra-short laser pulses can be used to induce vapor bubble nucleation, thus creating a stable two-phase fluid inclusion appropriate for subsequent measurements of the liquid-vapor homogenization temperature, T-h. In this study we evaluate the applicability of T-h data to accurately determine gypsum formation temperatures. We used fluid inclusions in synthetic gypsum crystals grown in the laboratory at different temperatures between 40 degrees C and 80 degrees C under atmospheric pressure conditions. We found an asymmetric distribution of the T-h values, which are systematically lower than the actual crystal growth temperatures, T-g; this is due to (1) the effect of surface tension on liquid-vapor homogenization, and (2) plastic deformation of the inclusion walls due to internal tensile stress occurring in the metastable state of the inclusions. Based on this understanding, we have determined growth temperatures of natural giant gypsum crystals from Naica (Mexico), yielding 47 +/- 1.5 degrees C for crystals grown in the Cave of Swords (120 m below surface) and 54.5 +/- 2 degrees C for giant crystals grown in the Cave of Crystals (290 m below surface). These results support the earlier hypothesis that the population and the size of the Naica crystals were controlled by temperature. In addition, this experimental method opens a door to determining the growth temperature of minerals forming in low-temperature environments.

In vivo compatibility of Dynesys(®) spinal implants: a case series of five retrieved periprosthetic tissue samples and corresponding implants

PURPOSE To determine whether particulate debris is present in periprosthetic tissue from revised Dynesys(®) devices, and if present, elicits a biological tissue reaction. METHODS Five Dynesys(®) dynamic stabilization systems consisting of pedicle screws (Ti alloy), polycarbonate-urethane (PCU) spacers and a polyethylene-terephthalate (PET) cord were explanted for pain and screw loosening after a mean of 2.86 years (1.9-5.3 years). Optical microscopy and scanning electron microscopy were used to evaluate wear, deformation and surface damage, and attenuated total reflectance Fourier transform infrared spectroscopy to assess surface chemical composition of the spacers. Periprosthetic tissue morphology and wear debris were determined using light microscopy, and PCU and PET wear debris by polarized light microscopy. RESULTS All implants had surface damage on the PCU spacers consistent with scratches and plastic deformation; 3 of 5 exhibited abrasive wear zones. In addition to fraying of the outer fibers of the PET cords in five implants, one case also evidenced cord fracture. The pedicle screws were unremarkable. Patient periprosthetic tissues around the three implants with visible PCU damage contained wear debris and a corresponding macrophage infiltration. For the patient revised for cord fracture, the tissues also contained large wear particles (>10 μm) and giant cells. Tissues from the other two patients showed comparable morphologies consisting of dense fibrous tissue with no inflammation or wear debris. CONCLUSIONS This is the first study to evaluate wear accumulation and local tissue responses for explanted Dynesys(®) devices. Polymer wear debris and an associated foreign-body macrophage response were observed in three of five cases.

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.

Temperature- and strain-rate-dependent microfabric evolution in monomineralic mylonite: evidence from in situ deformation of norcamphor

Norcamphor (C7H10O) was subjected to plane strain simple shear in a see-through deformation rig at four different strain rate and temperature conditions. Two transient stages in the microfabric evolution to steady state are distinguished. The grain scale mechanisms associated with the microstructural and textural evolution vary with the applied temperature, strain rate and strain. In high-temperature-low-strain-rate experiments, computer integrated polarization microscopy reveals that the texture evolution is closely related to the crystallographic rotation paths and rotation rates of individual grains. High c-axis rotation rates at low to intermediate shear strains are related to the development of a symmetrical c-axis cross girdle by the end of the first transient stage (γ = 1.5 to 2). During the second transient stage (γ = 1.5 to 6), the cross girdle yields to an oblique c-axis single girdle as c-axis rotation rates decrease and the relative activity of grain boundary migration recrystallization increases. Steady state (γ > 8) is characterized by a stable end orientation of the sample texture and the cyclic growth, rotation and consumption of individual grains within the aggregate.

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.