Cortical thickness mapping to identify focal osteoporosis in patients with hip fracture.

BACKGROUND: Individuals with osteoporosis are predisposed to hip fracture during trips, stumbles or falls, but half of all hip fractures occur in those without generalised osteoporosis. By analysing ordinary clinical CT scans using a novel cortical thickness mapping technique, we discovered patches of markedly thinner bone at fracture-prone regions in the femurs of women with acute hip fracture compared with controls. METHODS: We analysed CT scans from 75 female volunteers with acute fracture and 75 age- and sex-matched controls. We classified the fracture location as femoral neck or trochanteric before creating bone thickness maps of the outer 'cortical' shell of the intact contra-lateral hip. After registration of each bone to an average femur shape and statistical parametric mapping, we were able to visualise and quantify statistically significant foci of thinner cortical bone associated with each fracture type, assuming good symmetry of bone structure between the intact and fractured hip. The technique allowed us to pinpoint systematic differences and display the results on a 3D average femur shape model. FINDINGS: The cortex was generally thinner in femoral neck fracture cases than controls. More striking were several discrete patches of statistically significant thinner bone of up to 30%, which coincided with common sites of fracture initiation (femoral neck or trochanteric). INTERPRETATION: Femoral neck fracture patients had a thumbnail-sized patch of focal osteoporosis at the upper head-neck junction. This region coincided with a weak part of the femur, prone to both spontaneous 'tensile' fractures of the femoral neck, and as a site of crack initiation when falling sideways. Current hip fracture prevention strategies are based on case finding: they involve clinical risk factor estimation to determine the need for single-plane bone density measurement within a standard region of interest (ROI) of the femoral neck. The precise sites of focal osteoporosis that we have identified are overlooked by current 2D bone densitometry methods.

Fracture toughness of fibrous membranes

Random fibrous networks exist in both natural biological and engineering materials. While the nonlinear deformation of fibrous networks has been extensively studied, the understanding of their fracture behaviour is still incomplete. To study the fracture toughness of fibrous materials, the near-tip region is crucial because failure mechanisms such as fibril rupture occur in this region. The consideration of this region in fracture studies is, however, a difficult task because it involves microscopic mechanical responses at a small length scale. This paper extends our previous finite element analysis by incorporating the microscopic responses into a macroscopic domain by using a submodeling technique. The detailed study of microstructures at crack tips show a stochastic toughness of membranes due to the random nature of fibrous networks. Further, the sizes of crack tip region, which are sufficient to provide a reasonable prediction of fracture behaviour in a specific type of fibrous network, were presented. Future work includes improving the current linear assumption in the macroscopic models to become nonlinear.

Deformation and fracture of impulsively loaded sandwich panels

Light metal sandwich panel structures with cellular cores have attracted interest for multifunctional applications which exploit their high bend strength and impact energy absorption. This concept has been explored here using a model 6061-T6 aluminum alloy system fabricated by friction stir weld joining extruded sandwich panels with a triangular corrugated core. Micro-hardness and miniature tensile coupon testing revealed that friction stir welding reduced the strength and ductility in the welds and a narrow heat affected zone on either side of the weld by approximately 30%. Square, edge clamped sandwich panels and solid plates of equal mass per unit area were subjected to localized impulsive loading by the impact of explosively accelerated, water saturated, sand shells. The hydrodynamic load and impulse applied by the sand were gradually increased by reducing the stand-off distance between the test charge and panel surfaces. The sandwich panels suffered global bending and stretching, and localized core crushing. As the pressure applied by the sand increased, face sheet fracture by a combination of tensile stretching and shear-off occurred first at the two clamped edges of the panels that were parallel with the corrugation and weld direction. The plane of these fractures always lay within the heat affected zone of the longitudinal welds. For the most intensively loaded panels additional cracks occurred at the other clamped boundaries and in the center of the panel. To investigate the dynamic deformation and fracture processes, a particle-based method has been used to simulate the impulsive loading of the panels. This has been combined with a finite element analysis utilizing a modified Johnson-Cook constitutive relation and a Cockcroft-Latham fracture criterion that accounted for local variation in material properties. The fully coupled simulation approach enabled the relationships between the soil-explosive test charge design, panel geometry, spatially varying material properties and the panel's deformation and dynamic failure responses to be explored. This comprehensive study reveals the existence of a strong instability in the loading that results from changes in sand particle reflection during dynamic evolution of the panel's surface topology. Significant fluid-structure interaction effects are also discovered at the sample sides and corners due to changes of the sand reflection angle by the edge clamping system. © 2012 Elsevier Ltd. All rights reserved.

Experimental and numerical investigation of fracture mechanics of bitumen beams

This work presents an experimental and numerical investigation to characterise the fracture properties of pure bitumen (the binder in asphalt paving materials). The paper is divided into two parts. The first part describes an experimental study of fracture characterisation parameters of pure bitumen as determined by three-point bend tests. The second part deals with modelling of fracture and failure of bitumen by Finite Element analysis. Fracture mechanics parameters, stress intensity factor, KIC, fracture energy, GIC, and J-integral, JIC, are used for evaluation of bitumen's fracture properties. The material constitutive model developed by Ossa et al. [40,41] which was implemented into a FE code by Costanzi [18] is combined with cohesive zone models (CZM) to simulate the fracture behaviour of pure bitumen. Experimental and numerical results are presented in the form of failure mechanism maps where ductile, brittle and brittle-ductile transition regimes of fracture behaviour are classified. The FE predictions of fracture behaviour match well with experimental results. © 2012 Elsevier Ltd.

Fracture energy of the concrete-FRP interface in strengthened beams

To determine the load at which FRPs debond from concrete beams using global-energy-balance-based fracture mechanics concepts, the single most important parameter is the fracture energy of the concrete-FRP interface, which is easy to define but difficult to determine. Debonding propagates in the narrow zone of concrete, between the FRP and the (tension) steel reinforcement bars in the beam, and the presence of nearby steel bars prevents the fracture process zone, which in concrete is normally extensive, from developing fully. The paper presents a detailed discussion of the mechanism of the FRP debonding, and shows that the initiation of debonding can be regarded as a Mode I (tensile) fracture in concrete, despite being loaded primarily in shear. It is shown that the incorporation of this fracture energy in the debonding model developed by the authors, details of which are presented elsewhere, gives predictions that match the test results reported in the literature. © 2013 Elsevier Ltd.

Modeling mode i fracture of bitumen films

The fracture behavior of thin films of bitumen in double cantilever beam (DCB) specimens was investigated over a wide range of temperature and loading rate conditions using finite-element analysis. The model includes a phenomenological model for the mechanical behavior of bitumen, implemented into a special-purpose finite-element user material subroutine, combined with a cohesive zone model (CZM) for simulating the fracture process. The finite-element model is validated against experimental results from laboratory tests of DCB specimens by comparing measured and predicted load-line deflection histories and fracture energy release rates. Computer simulation results agreed well with experimental data of DCB joints containing bitumen films in terms of peak stress, fracture toughness, and stress-strain history response. The predicted "normalized toughness," G=2h, was found to increase in a power-law manner with effective temperaturecompensated strain rate in the ductile region as previously observed experimentally. In the brittle regime, G=2h is virtually constant. The model successfully captured the ductile and brittle failure behavior of bitumen films in opening mode (tension) for stable crack growth conditions. © 2013 American Society of Civil Engineers.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.