Predicting the Fracture Toughness of Human Cancellous Bone in Fractured Neck of Femur Patients Using Bone Volume and Micro-Architecture.

The current protocol used to determine if an individual is osteoporotic relies on assessment of the individual's bone mineral density (BMD), which allows clinicians to judge the condition of a patient with respect to their peers. This, in essence, evaluates a person's fracture risk, because BMD is a good surrogate measure for strength and stiffness. In recent studies, the authors were the first to produce fracture toughness (FT) data from osteoporotic (OP) and osteoarthritic (OA) patients, by using a testing technique which basically analyzes the prerequisite stress conditions for the onset of growth of a major crack through cancellous bone tissue. FT depends mainly on bone quantity (BV/TV, bone volume/tissue volume), but also on bone micro-architecture (mArch), the inner trabecular design of the bone. The working research hypothesis of the present study is that mArch offers added prediction power to BV/TV in determining FT parameters. Consequently, our aim was to investigate the use of predictive models for fracture toughness and also to investigate if there are any significant differences between the models produced from samples loaded across (AC, transverse to) the main trabecular orientation and along (AL, in parallel) the trabeculae. In multilinear regression analysis, we found that the strength of the relationships varied for a crack growing in these two orthogonal directions. Adding mArch variables in the Ac direction helped to increase the R2 to 0.798. However, in the AL direction, adding the mArch parameters did not add any predictive power to using BV/TV alone; BV/TV on its own could produce R2 = 0.730. The present results also imply that the anisotropic layout of the trabeculae makes it more difficult for a major crack to grow transversely across them. Cancellous bone models and remodels itself in a certain way to resist fracture in a specific direction, and thus, we should be mindful that architectural quality as well as bone quantity are needed to understand the resistance to fracture.

Parametrization of the Calcaneus and Medial Cuneiform to Aid Potential Advancements in Flatfoot Surgery.

INTRODUCTION: Flatfoot is a condition commonly seen in children; however, there is general disagreement over its incidence, characterization and correction. Painful flatfoot accompanied with musculoskeletal and soft tissue problems requires surgery to avoid arthritis in adulthood, the most common surgical approach being two osteotomies to the calcaneus and medial cuneiform bones of the foot. OBJECTIVES: This study focuses on the parametrization of these two bones to understand their bone morphology differences in a population sample among 23 normal subjects. Population differences could help in understanding whether bone shape may be an important factor in aiding surgical planning and outcomes. METHODS: A total of 45 sets of CT scans of these subjects were used to generate surface meshes of the two bones and converted to be iso-topological meshes, simplifying the application of Generalized Procrustes Analysis and Principal Component Analysis, allowing the main sources of variation between the subjects to be quantified. RESULTS: For the calcaneus, 16 Principal Components (PCs) and, for the medial cuneiform, 12 PCs were sufficient to describe 90% of the dataset variability. The quantitative and qualitative analyses confirm that for the calcaneus PC1 describes the Achilles attachment location and PC2 largely describes the anterior part of the bone. For the medial cuneiform, PC1 describes the medial part of the bone, while PC2 mainly describes the superior part. CONCLUSION: Most importantly, the PCs did not seem to describe the osteotomy sites for both bones, suggesting low population variability at the bone cutting points. Further studies are needed to evaluate how shape variability impacts surgical outcomes. Future implications could include better surgical planning and may pave the way for complex robotic surgeries to become a reality.

Quantifying microcracks on fractured bone surfaces - Potential use in forensic anthropology.

Bone fracture surface morphology (FSM) can provide valuable information on the cause of failure in forensic and archaeological applications and it depends primarily on three factors, the loading conditions (like strain rate), the ambient conditions (wet or dry bone material) and the quality of bone material itself. The quality of bone material evidently changes in taphonomy as a result of the decomposition process and that in turn is expected to affect FSM. Porcine bones were fractured by a standardised impact during the course of soft tissue decomposition, at 28-day intervals, over 140 days (equivalent to 638 cooling degree days). Measurements of the associated microcracks on the fractured cortical bone surfaces indicated a progressive increase in mean length during decomposition from around 180 µm-375 µm. The morphology of these microcracks also altered, from multiple intersecting microcracks emanating from a central point at 0-28 cumulative cooling degree days, to longer linear cracks appearing to track lamellae as soft tissue decomposition progressed. The implications of these findings are that taphonomic changes of bone may offer the real possibility of distinguishing perimortem and taphonomic damage and also provide a new surrogate parameter for estimation of post-mortem interval (PMI) in forensics.

Association between nanoscale strains and tissue level nanoindentation properties in age-related hip-fractures.

Measurement of the properties of bone as a material can happen in various length scales in its hierarchical and composite structure. The aim of this study was to test the tissue level properties of clinically-relevant human bone samples which were collected from donors belonging to three groups: ageing donors who suffered no fractures (Control); untreated fracture patients (Fx-Untreated) and patient who experienced hip fracture despite being treated with bisphosphonates (Fx-BisTreated). Tissue level properties were assessed by (a) nanoindentation and (b) synchrotron tensile tests (STT) where strains were measured at the 'tissue', 'fibril' and 'mineral' levels by using simultaneous Wide-angle - (WAXD) and Small angle- X-ray diffraction (SAXD). The composition was analysed by thermogravimetric analysis and material level endo- and exo-thermic reactions by differential scanning calorimetry (TGA/DSC3+). Irrespective of treatment fracture donors exhibited significantly lower tissue, fibril and mineral strain at the micro and nanoscale respectively and had a higher mineral content than controls. In nanoindentation only nanohardness was significantly greater for Controls and Fx-BisTreated versus Fx-Untreated. The other nanoindentation parameters did not vary significantly across the three groups. There was a highly significant positive correlation (p < 0.001) between organic content and tissue level strain behaviour. Overall hip-fractures were associated with lower STT nanostrains and it was behaviour measured by STT which proved to be a more effective approach for predicting fracture risk because evidently it was able to demonstrate the mechanical deficit for the bone tissue of the donors who had experienced fractures.

The interplay between BMU activity linked to mechanical stress, specific surface and inhibitory theory dictate bone mass distribution: Predictions from a 3D computational model.

Bone mechanical and biological properties are closely linked to its internal tissue composition and mass distribution, which are in turn governed by the purposeful action of the basic multicellular units (BMUs). The orchestrated action of osteoclasts and osteoblasts, the resorbing and forming tissue cells respectively, in BMUs is responsible for tissue maintenance, repair and adaptation to changing load demands through the phenomenon known as remodelling. In this work, a computational mechano-biological model of bone remodelling based on the inhibitory theory and a new scheme of bone resorption introduced previously in a 2D model, is extended to a 3D model of the real external geometry of a femur under normal walking loads. Starting from a uniform apparent density (ratio of tissue local mass to total local volume) distribution, the BMU action can be shown to lead naturally to an internal density distribution similar to that of a real bone, provided that the initial density value is high enough to avoid unrealistic final mass deposition in zones of high energy density and excessive damage. Physiological internal density values are reached throughout the whole 3D geometry, and at the same time a 'boomerang'-like relationship between apparent and material density (ratio of tissue mass to tissue volume) emerges naturally under the proposed remodelling scheme. It is also shown here that bone-specific surface is a key parameter that determines the intensity of BMU action linked to the mechanical and biological requirements. Finally, by engaging in simulations of bone in disuse, we were able to confirm the appropriate selection of the model parameters. As an example, our results show good agreement with experimental measurements of bone mass on astronauts a fact that strengthens our belief in the insightful nature of our novel 3D computational model.

Assessing bone maturity: Compositional and mechanical properties of rib cortical bone at different ages.

Understanding what maturity entails for bone, when it arrives, and its pre- and post-maturity traits and properties are very important for understanding its evolution and physiology. There is a clear but fine distinction between the chronological age of bone (the age of its donor) and the tissue age of the bone packets it comprises at the microscopic level. Whole bone fragility changes with age due to mass and architecture effects, but so do the properties of bone at the tissue level. Tissue age and tissue-level properties are therefore increasingly attracting a great deal of attention recently. The present study investigated compositional and material changes in the hydroxyapatite crystals, the collagenous phase, changes in bone matrix composition and its nanoindentation properties and their decline with chronological age in later life. The aim was to track the age threshold at which cortical bone arrives at maturity and what happens following that threshold. To do so FTIR, DSC/TGA, XRD, nanoindentation and microindentation were used to investigate rib cortical bone material across a cohort of 86 individuals from one ethnic group with age spanning between 17 and 82 years. Results of this cross-sectional study showed a clear increase in mineral content relative to the organic and water contents across all ages. Furthermore, an increase in crystal size and consequent decrease in strain (coherence length) was detected associated with secondary mineralisation and an increase in carbonate substitution. Overall, we observe a number of modifications which contribute to a typical functional behaviour of bone showing an increase in both indentation modulus and hardness until the age of about 35 after which both of these properties decline gradually and concomitantly to other physicochemical changes and seemingly until the end of one's life.

Age related changes of rib cortical bone matrix and the application to forensic age-at-death estimation.

Forensic anthropology includes, amongst other applications, the positive identification of unknown human skeletal remains. The first step in this process is an assessment of the biological profile, that is: sex, age, stature and ancestry. In forensic contexts, age estimation is one of the main challenges in the process of identification. Recently established admissibility criteria are driving researchers towards standardisation of methodological procedures. Despite these changes, experience still plays a central role in anthropological examinations. In order to avoid this issue, age estimation procedures (i) must be presented to the scientific community and published in peer reviewed journals, (ii) accurately explained in terms of procedure and (iii) present clear information about the accuracy of the estimation and possible error rates. In order to fulfil all these requirements, a number of methods based on physiological processes which result in biochemical changes in various tissue structures at the molecular level, such as modifications in DNA-methylation and telomere shortening, racemization of proteins and stable isotopes analysis, have been developed. The current work proposes a new systematic approach in age estimation based on tracing physicochemical and mechanical degeneration of the rib cortical bone matrix. This study used autopsy material from 113 rib specimens. A set of 33 parameters were measured by standard bio-mechanical (nanoindentation and microindentation), physical (TGA/DSC, XRD and FTIR) and histomorphometry (porosity-ImageJ) methods. Stepwise regressions were used to create equations that would produce the best 'estimates of age at death' vs real age of the cadavers. Five equations were produced; in the best of cases an equation counting 7 parameters had an R2 = 0.863 and mean absolute error of 4.64 years. The present method meets all the admissibility criteria previously described. Furthermore, the method is experience-independent and as such can be performed without previous expert knowledge of forensic anthropology and human anatomy.

Evaluation of bone excision effects on a human skull model - I: Mechanical testing and digital image correlation.

The mechanisms of skull impact loading may change following surgical interventions such as the removal of bone lesions, but little is known about the consequences in the event of subsequent head trauma. We, therefore, prepared acrylonitrile butadiene styrene human skull models based on clinical computed tomography skull data using a three-dimensional printer. Six replicate physical skull models were tested, three with bone excisions and three without. A drop tower was used to simulate the impact sustained by falling backwards onto the occipital lobe region. The impacts were recorded with a high-speed camera, and the occipital strain response was determined by digital image correlation. Although the hole affected neither the magnitude nor the sequence of the fracture pattern, the digital image correlation analysis highlighted an increase in strain around the excised area (0.45%-16.4% of the principal strain). Our approach provides a novel method that could improve the quality of life for patients on many fronts, including protection against trauma, surgical advice, post-operative care, advice in litigation cases, as well as facilitating general biomechanical research in the area of trauma injuries.

Evaluation of bone excision effects on a human skull model-II: Finite element analysis.

Patient-specific computational models are powerful tools which may assist in predicting the outcome of invasive surgery on the musculoskeletal system, and consequently help to improve therapeutic decision-making and post-operative care. Unfortunately, at present the use of personalized models that predict the effect of biopsies and full excisions is so specialized that tends to be restricted to prominent individuals, such as high-profile athletes. We have developed a finite element analysis model to determine the influence of the location of an ellipsoidal excision (14.2 mm × 11.8 mm) on the structural integrity of a human skull when exposed to impact loading, representing a free fall of an adult male from standing height. The finite element analysis model was compared to empirical data based on the drop-tower testing of three-dimensional-printed physical skull models where deformations were recorded by digital image correlation. In this bespoke example, we found that the excision site did not have a major effect on the calculated stress and strain magnitudes unless the excision was in the temporal region, where the reduction in stiffness around the excision caused failure within the neighboring area. The finite element analysis model allowed meaningful conclusions to be drawn for the implications of using such a technique based on what we know about such conditions indicating that the approach could be both clinically beneficial and also cost-effective for wider use.

Ageing bone fractures: The case of a ductile to brittle transition that shifts with age.

Human bone becomes increasingly brittle with ageing. Bones also fracture differently under slow and fast loadings, being ductile and brittle, respectively. The effects of a combination of these two factors have never been examined before. Here we show that cortical bone is most fracture-resistant at the physiologically prevalent intermediate strain rates of 10-3 s-1 to 10-2 s-1 such as they occur in walking or running, slightly weaker at slower quasistatic and much weaker at fast impact loading rates. In young cortical bone (15 years of age) the ductile-to-brittle transition (DBT) occurs at strain rates of 10-2 s-1, in old cortical bone (85 yrs) at speeds lower by a factor of 10 to 40. Other research has shown that the energy required to break bone (per unit of fracture surface) drops as much as 60% between these two ages. Therefore, DBT seems to compound the well-known phenomenon of 'brittle old bones'. Old bones can only cope with slow movement, young ones with both slow and fast movement. These observed material characteristics of (i) a shift of the DBT and (ii) a reduced energy absorption capacity appear to contribute at least as much to the loss of bone quality as the various quantity based (lowered bone density and mineral content) explanations of the past. They also provide a new powerful paradigm, which allows us to demonstrate mechanically, and uniquely, how human bone becomes increasingly brittle with age.

Age-Related Trends in the Trabecular Micro-Architecture of the Medial Clavicle: Is It of Use in Forensic Science?

The mechanical and structural properties of bone are known to change significantly with age. Within forensic and archaeological investigations, the medial end of the clavicle is typically used for estimating the age-at-death of an unknown individual. Although, this region of the skeleton is of interest to forensic and clinical domains, alterations beyond the macro-scale have not been fully explored. For this study, non-destructive micro-computed tomography (µ-CT) was employed to characterize structural alterations to the cancellous bone of the medial clavicle. Fresh human cadaveric specimens (12-59 years) obtained at autopsy were utilized for this study, and were scanned with a voxel size of ~83 µm. Morphometric properties were quantified and indicated that the bone volume, connectivity density, mineral density, and number of trabeculae decreased with age, while the spacing between the trabeculae increased with age. In contrast to other sub-regions of the skeleton, trabecular thickness, and degree of anisotropy did not correlate with age. Collectively, this could suggest that the network is becoming increasingly perforated with age rather than exhibiting trabecular thinning. These results are used in the context of deriving a potential protocol for forensic investigations by using this particular and largely unexplored region of the skeleton, and provide inspiration for future experiments concerning micro-architectural and small scale changes in other regions of the human skeleton.

Age-Related Changes in Femoral Head Trabecular Microarchitecture.

Osteoporosis is a prevalent bone condition, characterised by low bone mineral density and increased fracture risk. Currently, the gold standard for identifying osteoporosis and increased fracture risk is through quantification of bone mineral density using dual energy X-ray absorption. However, many studies have shown that bone strength, and consequently the probability of fracture, is a combination of both bone mass and bone 'quality' (architecture and material chemistry). Although the microarchitecture of both non-fracture and osteoporotic bone has been previously investigated, many of the osteoporotic studies are constrained by factors such as limited sample number, use of ovariectomised animal models, and lack of male and female discrimination. This study reports significant differences in bone quality with respect to the microarchitecture between fractured and non-fractured human femur specimens. Micro-computed tomography was utilised to investigate the microarchitecture of femoral head trabecular bone from a relatively large cohort of non-fracture and fracture human donors. Various microarchitectural parameters have been determined for both groups, providing an understanding of the differences between fracture and non -fracture material. The microarchitecture of non-fracture and fracture bone tissue is shown to be significantly different for many parameters. Differences between sexes also exist, suggesting differences in remodelling between males and females in the fracture group. The results from this study will, in the future, be applied to develop a fracture model which encompasses bone density, architecture and material chemical properties for both female and male tissues.

Early arthritis induces disturbances at bone nanostructural level reflected in decreased tissue hardness in an animal model of arthritis.

INTRODUCTION: Arthritis induces joint erosions and skeletal bone fragility. OBJECTIVES: The main goal of this work was to analyze the early arthritis induced events at bone architecture and mechanical properties at tissue level. METHODS: Eighty-eight Wistar rats were randomly housed in experimental groups, as follows: adjuvant induced arthritis (AIA) (N = 47) and a control healthy group (N = 41). Rats were monitored during 22 days for the inflammatory score, ankle perimeter and body weight and sacrificed at different time points (11 and 22 days post disease induction). Bone samples were collected for histology, micro computed tomography (micro-CT), 3-point bending and nanoindentation. Blood samples were also collected for bone turnover markers and systemic cytokine quantification. RESULTS: At bone tissue level, measured by nanoindentation, there was a reduction of hardness in the arthritic group, associated with an increase of the ratio of bone concentric to parallel lamellae and of the area of the osteocyte lacuna. In addition, increased bone turnover and changes in the microstructure and mechanical properties were observed in arthritic animals, since the early phase of arthritis, when compared with healthy controls. CONCLUSION: We have shown in an AIA rat model that arthritis induces very early changes at bone turnover, structural degradation and mechanical weakness. Bone tissue level is also affected since the early phase of arthritis, characterized by decreased tissue hardness associated with changes in bone lamella organization and osteocyte lacuna surface. These observations highlight the pertinence of immediate control of inflammation in the initial stages of arthritis.

Effects of tofacitinib in early arthritis-induced bone loss in an adjuvant-induced arthritis rat model.

Objectives: The main goal of this work was to analyse how treatment intervention with tofacitinib prevents the early disturbances of bone structure and mechanics in the rat model of adjuvant-induced arthritis. This is the first study to access the impact of tofacitinib on the skeletal bone effects of inflammation. Methods: Fifty Wistar rats with adjuvant-induced arthritis were randomly housed in experimental groups, as follows: non-arthritic healthy group (n = 20); arthritic non-treated group (n = 20); and 10 animals undergoing tofacitinib treatment. Rats were monitored during 22 days after disease induction for the inflammatory score, ankle perimeter and body weight. Healthy non-arthritic rats were used as controls for comparison. After 22 days of disease progression, rats were killed and bone samples collected for histology, micro-CT, three-point bending and nanoindentation analysis. Blood samples were also collected for quantification of bone turnover markers and systemic cytokines. Results: At the tissue level, measured by nanoindentation, tofacitinib increased bone cortical and trabecular hardness. However, micro-CT and three-point bending tests revealed that tofacitinib did not reverse the effects of arthritis on the cortical and trabecular bone structure and on mechanical properties. Conclusion: Possible reasons for these observations might be related to the mechanism of action of tofacitinib, which leads to direct interactions with bone metabolism, and/or to the kinetics of its bone effects, which might need longer exposure.

'Phoenix polymers': fire induced nanohardness in fibril-forming aromatic cyanate esters.

For the first time we present nanoindentation analysis of charred, cured aromatic cyanate esters, which exhibit outstanding mechanical properties when analysed under applied loads of 0.1-300 mN. Following charring (900 °C for 10 minutes to achieve graphitised structures), the samples display a remarkable combination of a modulus of elasticity of around 25 GPa and nanohardness of 300 kgf mm-2, making them some 30-40% stiffer than bone and practically as hard as tooth enamel. At the same time we find that under the same conditions the chars are highly resilient, displaying complete elastic recovery with very little plastic deformation. When cured in the presence of copper(ii) acetylacetonate (200 ppm) in dodecylphenol (1% w/v active copper suspension) to form a polycyanurate, compound (2) forms a dense, consolidated structure compared with compound (1) under the same conditions. At high magnification, the presence of a nanoscale, fibrillar structure is observed, accounting for the high resilience.

Rib biomechanical properties exhibit diagnostic potential for accurate ageing in forensic investigations.

Age estimation remains one of the most challenging tasks in forensic practice when establishing a biological profile of unknown skeletonised remains. Morphological methods based on developmental markers of bones can provide accurate age estimates at a young age, but become highly unreliable for ages over 35 when all developmental markers disappear. This study explores the changes in the biomechanical properties of bone tissue and matrix, which continue to change with age even after skeletal maturity, and their potential value for age estimation. As a proof of concept we investigated the relationship of 28 variables at the macroscopic and microscopic level in rib autopsy samples from 24 individuals. Stepwise regression analysis produced a number of equations one of which with seven variables showed an R2 = 0.949; a mean residual error of 2.13 yrs ±0.4 (SD) and a maximum residual error value of 2.88 yrs. For forensic purposes, by using only bench top machines in tests which can be carried out within 36 hrs, a set of just 3 variables produced an equation with an R2 = 0.902 a mean residual error of 3.38 yrs ±2.6 (SD) and a maximum observed residual error 9.26yrs. This method outstrips all existing age-at-death methods based on ribs, thus providing a novel lab based accurate tool in the forensic investigation of human remains. The present application is optimised for fresh (uncompromised by taphonomic conditions) remains, but the potential of the principle and method is vast once the trends of the biomechanical variables are established for other environmental conditions and circumstances.

Sharp and blunt force trauma concealment by thermal alteration in homicides: An in-vitro experiment for methodology and protocol development in forensic anthropological analysis of burnt bones.

Burning of human remains is one method used by perpetrators to conceal fatal trauma and expert opinions regarding the degree of skeletal evidence concealment are often disparate. This experiment aimed to reduce this incongruence in forensic anthropological interpretation of burned human remains and implicitly contribute to the development of research methodologies sufficiently robust to withstand forensic scrutiny in the courtroom. We have tested the influence of thermal alteration on pre-existing sharp and blunt trauma on twenty juvenile sheep radii in the laboratory using an automated impact testing system and an electric furnace. The testing conditions simulated a worst-case scenario where remains with pre-existing sharp or blunt trauma were exposed to burning with an intentional vehicular fire scenario in mind. All impact parameters as well as the burning conditions were based on those most commonly encountered in forensic cases and maintained constant throughout the experiment. The results have shown that signatures associated with sharp and blunt force trauma were not masked by heat exposure and highlights the potential for future standardization of fracture analysis in burned bone. Our results further emphasize the recommendation given by other experts on handling, processing and recording burned remains at the crime scene and mortuary.

Localized tissue mineralization regulated by bone remodelling: A computational approach.

Bone is a living tissue whose main mechanical function is to provide stiffness, strength and protection to the body. Both stiffness and strength depend on the mineralization of the organic matrix, which is constantly being remodelled by the coordinated action of the bone multicellular units (BMUs). Due to the dynamics of both remodelling and mineralization, each sample of bone is composed of structural units (osteons in cortical and packets in cancellous bone) created at different times, therefore presenting different levels of mineral content. In this work, a computational model is used to understand the feedback between the remodelling and the mineralization processes under different load conditions and bone porosities. This model considers that osteoclasts primarily resorb those parts of bone closer to the surface, which are younger and less mineralized than older inner ones. Under equilibrium loads, results show that bone volumes with both the highest and the lowest levels of porosity (cancellous and cortical respectively) tend to develop higher levels of mineral content compared to volumes with intermediate porosity, thus presenting higher material densities. In good agreement with recent experimental measurements, a boomerang-like pattern emerges when plotting apparent density at the tissue level versus material density at the bone material level. Overload and disuse states are studied too, resulting in a translation of the apparent-material density curve. Numerical results are discussed pointing to potential clinical applications.

Development of a strain rate dependent material model of human cortical bone for computer-aided reconstruction of injury mechanisms.

Computer-aided methods such as finite-element simulation offer a great potential in the forensic reconstruction of injury mechanisms. Numerous studies have been performed on understanding and analysing the mechanical properties of bone and the mechanism of its fracture. Determination of the mechanical properties of bones is made on the same basis used for other structural materials. The mechanical behaviour of bones is affected by the mechanical properties of the bone material, the geometry, the loading direction and mode and of course the loading rate. Strain rate dependency of mechanical properties of cortical bone has been well demonstrated in literature studies, but as many of these were performed on animal bones and at non-physiological strain rates it is questionable how these will apply in the human situations. High strain-rates dominate in a lot of forensic applications in automotive crashes and assault scenarios. There is an overwhelming need to a model which can describe the complex behaviour of bone at lower strain rates as well as higher ones. Some attempts have been made to model the viscoelastic and viscoplastic properties of the bone at high strain rates using constitutive mathematical models with little demonstrated success. The main objective of the present study is to model the rate dependent behaviour of the bones based on experimental data. An isotropic material model of human cortical bone with strain rate dependency effects is implemented using the LS-DYNA material library. We employed a human finite element model called THUMS (Total Human Model for Safety), developed by Toyota R&D Labs and the Wayne State University, USA. The finite element model of the human femur is extracted from the THUMS model. Different methods have been employed to develop a strain rate dependent material model for the femur bone. Results of one the recent experimental studies on human femur have been employed to obtain the numerical model for cortical femur. A forensic application of the model is explained in which impacts to the arm have been reconstructed using the finite element model of THUMS. The advantage of the numerical method is that a wide range of impact conditions can be easily reconstructed. Impact velocity has been changed as a parameter to find the tolerance levels of injuries to the lower arm. The method can be further developed to study the assaults and the injury mechanism which can lead to severe traumatic injuries in forensic cases.