A bone remodelling model coupling micro-damage growth and repair by 3D BMU-activity.

Bone as most of living tissues is able, during its entire lifetime, to adapt its internal microstructure and subsequently its associated mechanical properties to its specific mechanical and physiological environment in a process commonly known as bone remodelling. Bone is therefore continuously renewed and micro-damage, accumulated by fatigue or creep, is removed minimizing the risk of fracture. Nevertheless, bone is not always able to repair itself completely. Actually, if bone repairing function is slower than micro-damage accumulation, a type of bone fracture, usually known as "stress fracture", can finally evolve. In this paper, we propose a bone remodelling continuous model able to simulate micro-damage growth and repair in a coupled way and able therefore to predict the occurrence of "stress fractures". The biological bone remodelling process is modelled in terms of equations that describe the activity of basic multicellular units. The predicted results show a good correspondence with experimental and clinical data. For example, in disuse, bone porosity increases until an equilibrium situation is achieved. In overloading, bone porosity decreases unless the damage rate is so high that causes resorption or "stress fracture".

Transplantation of human fetal blood stem cells in the osteogenesis imperfecta mouse leads to improvement in multiscale tissue properties.

Osteogenesis imperfecta (OI or brittle bone disease) is a disorder of connective tissues caused by mutations in the collagen genes. We previously showed that intrauterine transplantation of human blood fetal stem/stromal cells in OI mice (oim) resulted in a significant reduction of bone fracture. This work examines the cellular mechanisms and mechanical bone modifications underlying these therapeutic effects, particularly examining the direct effects of donor collagen expression on bone material properties. In this study, we found an 84% reduction in femoral fractures in transplanted oim mice. Fetal blood stem/stromal cells engrafted in bones, differentiated into mature osteoblasts, expressed osteocalcin, and produced COL1a2 protein, which is absent in oim mice. The presence of normal collagen decreased hydroxyproline content in bones, altered the apatite crystal structure, increased the bone matrix stiffness, and reduced bone brittleness. In conclusion, expression of normal collagen from mature osteoblast of donor origin significantly decreased bone brittleness by improving the mechanical integrity of the bone at the molecular, tissue, and whole bone levels.

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.

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.

Using mesoscopic models to design strong and tough biomimetic polymer networks

Using computational modeling, we investigate the mechanical properties of polymeric materials composed of coiled chains, or "globules", which encompass a folded secondary structure and are cross-linked by labile bonds to form a macroscopic network. In the presence of an applied force, the globules can unfold into linear chains and thereby dissipate energy as the network is deformed; the latter attribute can contribute to the toughness of the material. Our goal is to determine how to tailor the labile intra- and intermolecular bonds within the network to produce material exhibiting both toughness and strength. Herein, we use the lattice spring model (LSM) to simulate the globules and the cross-linked network. We also utilize our modified Hierarchical Bell model (MHBM) to simulate the rupture and reforming of N parallel bonds. By applying a tensile deformation, we demonstrate that the mechanical properties of the system are sensitive to the values of N in and N out, the respective values of N for the intra- and intermolecular bonds. We find that the strength of the material is mainly controlled by the value of N out, with the higher value of N out providing a stronger material. We also find that, if N in is smaller than N out, the globules can unfold under the tensile load before the sample fractures and, in this manner, can increase the ductility of the sample. Our results provide effective strategies for exploiting relatively weak, labile interactions (e.g., hydrogen bonding or the thiol/disulfide exchange reaction) in both the intra- and intermolecular bonds to tailor the macroscopic performance of the materials. © 2011 American Chemical Society.

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.

Fracturing pressure in clay

Hydraulic fracturing in clayey soils can be triggered by either tensile or shear failure. In this paper, the physical meanings of various equations to predict fracture initiation pressure proposed in the past are discussed using the cavity expansion theory. In particular, when fracturing pressure is plotted against initial confining pressure, published laboratory test results as well as analytical models show a linear relationship. When the slope is close to 2, fracture is initiated by tensile failure of the clay, whereas when the slope is close to 1, it is initiated by shear failure of the clay. In this study, the analytical models, validated only on laboratory test data to date, were applied to unique data from field grouting work in which extensive soil fracturing was carried out to improve the mechanical characteristics of the soft silty clay underlying a bell tower in Venice, Italy. By a careful assessment of initial confining pressure in the field, the variation in recorded injection pressures with confining pressure was examined. Results suggest that the fractures at this site were likely to be initiated by shear failure of the clay, and the values were similar to what was predicted by the model with the shear failure criterion. © 2013 American Society of Civil Engineers.