Finite element prediction of proximal femur fracture pattern based on orthotropic behaviour law coupled to quasi-brittle damage

Simulation of cancer progression in bone by a virtual thermal flux with a case study on lumbar vertebrae with multiple myeloma

BACKGROUND

Two main problems examining the mechanism of cancer progression in the tissues using the computational models are lack of enough knowledge on the effective factors for such events in vivo environments and lack of specific parameters in the available computational models to simulate such complicated reactions.

METHODS

In this study, it was tried to simulate the progression of cancerous lesions in the bone tissues by an independent parameter from the anatomical and physiological characteristics of the tissues, so to degrade the orthotropic mechanical properties of the bone tissues, a virtual temperature was determined to be used by a well-known framework for simulation of damages in the composite materials. First, the reliability of the FE model to simulate hyperelastic response in the intervertebral discs (IVDs) and progressive failure in the bony components were verified by simulation of some In-Vitro tests, available in the literature. Then, the progression of the osteolytic damage was simulated in a clinical case with multiple myeloma in the lumbar vertebrae.

RESULTS

The FE model could simulate stress-shielding and diffusion of lesion in the posterior elements of the damaged vertebra which led to spinal stenosis. The load carrying shares associated with the anterior half and the posterior half of the examined vertebral body and the posterior elements were estimated equal to 41 %, 47 % and 12 %, respectively for the intact condition, that changed to 14 %, 16 % and 70 %, when lesion occupied one third of the vertebral body.

CONCLUSION

Correlation of the FE results with the deformation shapes, observed in the MRIs for the clinical case study, indicated appropriateness of the procedure, proposed for simulation of the progressive osteolytic damage in the vertebral segments. The future studies may follow simulation of tumor growth for various metastatic tissues using the method, established here.

Biomechanical evaluation of different strain judging criteria on the prediction precision of cortical bone fracture simulation under compression

Introduction: The principal strain or equivalent strain is mainly used in current numerical studies to determine the mechanical state of the element in the cortical bone finite element model and then perform fracture simulation. However, it is unclear which strain is more suitable for judging the element mechanical state under different loading conditions due to the lack of a general strain judging criterion for simulating the cortical bone fracture. Methods: This study aims to explore a suitable strain judging criterion to perform compressive fracture simulation on the rat femoral cortical bone based on continuum damage mechanics. The mechanical state of the element in the cortical bone finite element model was primarily assessed using the principal strain and equivalent strain separately to carry out fracture simulation. The prediction accuracy was then evaluated by comparing the simulated findings with different strain judging criteria to the corresponding experimental data. Results: The results showed that the fracture parameters predicted using the principal strain were closer to the experimental values than those predicted using the equivalent strain. Discussion: Therefore, the fracture simulation under compression was more accurate when the principal strain was applied to control the damage and failure state in the element. This finding has the potential to improve prediction accuracy in the cortical bone fracture simulation.

Investigation on the Differences in the Failure Processes of the Cortical Bone under Different Loading Conditions

Cortical bone is a transversely isotropic material, and the mechanical properties may be related to the loading direction on the osteon. Therefore, analyzing the differences in the failure processes of cortical bone under different loading conditions is necessary to explore the measures for reducing the incidence of fracture. In this study, to investigate the effects of different loading directions on the fracture performance in the cortical bone, a numerical method that could simultaneously simulate the failure processes in the cortical bone structure under compression and bending loads was established based on continuum damage mechanics theory. The prediction accuracy and feasibility of the numerical method were first verified by comparing with the corresponding experimental results. Then, the differences in the failure process and fracture performance of the same cortical bone structure under compression and bending loads were investigated. The simulation results indicated that for the same structure, the slip-open failure mode appeared under compression load, and the crack propagated along a certain angle to the loading direction; the tension-open failure mode appeared under bending load, and the crack propagated along the direction perpendicular to the loading direction. Meanwhile, the fracture load was greater and the fracture time was later in the compression than in the bending condition. These phenomena stated that discrepant failure processes and fracture patterns occurred in the same cortical bone structure under different loading conditions. The main reason may be related to the tension–compression asymmetry and transversely isotropic characteristics in the cortical bone material. The fracture simulations in the cortical bone under different loading conditions could improve the prediction accuracy in bone biomechanics and provide the prevention method for cortical bone damage and fracture.

Comparison of Cortical Bone Fracture Patterns Under Compression Loading Using Finite Element–Discrete Element Numerical Modeling Approach and Destructive Testing

Finite element analysis may not be the only method by which bone fracture initiation and propagation may be analyzed. This study compares fracture patterns generated from compression testing of bone to fracture patterns generated using a combination of both the finite element method (FEM) and discrete element method (DEM) as defined by the finite discrete element method (FDEM). Before testing, a three-dimensional bone model was developed using CT. Force and displacement data were collected during testing. The tested specimen was reimaged using CT. The solid model was discretized and material properties adjusted such that finite element-discrete element macro behavior matched the force-displacement data. A qualitative comparison of the fracture patterns demonstrates that FDEM can successfully be used to simulate and predict fracturing in bone, with this study representing the first time this has been done and reported.

Stochastic failure analysis of proximal femur using an isogeometric analysis based nonlocal gradient-enhanced damage model

BACKGROUND AND OBJECTIVE

Medical imaging-based finite element methods are more accurate tools for fracture risk prediction than the traditional aBMD based methods. However, these methods have drawbacks like geometric errors, high computational cost, mesh-dependent results, etc. In this article, the authors have proposed an isogeometric analysis-based nonlocal gradient-enhanced damage model to overcome some of these issues. Moreover, there are uncertainties in the values of input parameters for such analysis due to various measurement errors. Hence, stochastic analysis is performed to quantify the effect of these parametric uncertainties on the fracture behavior of the proximal femur.

METHODS

Computed Tomography images of a patient are used to create a 2D proximal femur model with a heterogeneous description of material properties. A numerical model based on gradient-enhanced nonlocal continuum damage mechanics is used for fracture analysis of proximal femur to overcome the issues related to mesh dependency in traditional continuum damage mechanics models. Further, a multipatch isogeometric solver is developed to solve the governing equations. Monte Carlo simulations are used to understand the effect of parametric uncertainties on the fracture behavior of the proximal femur.

RESULTS

The developed numerical framework is used to solve the fracture problem of proximal femur under single leg stance loading conditions. The obtained results are validated by comparing the load-displacement response and the crack path with that given in the literature. Stochastic analysis is performed by considering a ±5% variation in the elastic modulus, damage initiation strain, and the neck-shaft angle values.

CONCLUSION

The proposed numerical framework can correctly predict the damage initiation and propagation in a proximal femur. The results reveal that the heterogeneous nature of material properties of bone plays a significant role in determining the fracture characteristics of the proximal femur. Further, the results of the stochastic analysis reveal that the parametric uncertainties in the neck-shaft angle have a much more significant influence on the results of the analysis than the parametric uncertainties in the elastic modulus and damage initiation strain.

Subject specific finite element modelling of periprosthetic femoral fractures in different load cases

Periprosthetic femoral fractures (PFF) around total hip replacements are one of the biggest challenges for orthopaedic surgeons. To understand the risk factors and formation of these fractures, the development of a reliable finite element (FE) model incorporating bone failure is essential. Due to the anisotropic and complex hierarchical structure of bone, the mechanical behaviour under large strains is difficult to predict. In this study, a state-of-the-art subject specific FE modelling technique for bone is utilised to generate and investigate PFF. A bilinear constitutive law is applied to bone tissue in subject specific FE models of five human femurs which are virtually implanted with a straight hip stem to numerically analyse PFF. The material parameters of the models are expressed as a function of bone ash density and mapped node wise to the FE mesh. In this way the subject specific, heterogeneous structure of bone is mimicked. For material mapping of the parameters, computed tomography (CT) images of the original fresh-frozen femurs are used. Periprosthetic fractures are generated by deleting elements on the basis of a critical plastic strain failure criterion. The models are analysed under physiological and clinically relevant conditions in two different load cases re-enacting stumbling and a sideways fall on the hip. The results of the analyses are quantified with experimental data from previous work. With regard to fracture pattern, stiffness and failure load the simulations of the load case stumbling delivered the most stable and accurate results. In general, mapping of material properties was found to be an appropriate way to reproduce PFF with finite element models.

Damage Function of a Quasi-Brittle Material, Damage Rate, Acceleration and Jerk during Uniaxial Compression: Model and Application to Analysis of Trabecular Bone Tissue Destruction

A diversity of quasi-brittle materials can be observed in various engineering structures and natural objects (rocks, frozen soil, concrete, ceramics, bones, etc.). In order to predict the condition and safety of these objects, a large number of studies aimed at analyzing the strength of quasi-brittle materials has been conducted and presented in publications. However, at the modeling level, the problem of estimating the rate and acceleration of destruction of a quasi-brittle material under loading remains relevant. The purpose of the study was to substantiate the function of damage to a quasi-brittle material under uniaxial compression, determine the rate, acceleration and jerk of the damage process, and also to apply the results obtained to predicting the destruction of trabecular bone tissue. In accordance with the purpose of the study, the basic concepts of fracture mechanics and standard methods of mathematical modeling were used. The proposed model is based on the application of the previously obtained differentiable damage function without parameters. The results of the study are presented in the form of plots and analytical relations for computing the rate, acceleration and jerk of the damage process. Examples are given. The predicted peak of the combined effect of rate, acceleration and jerk of the damage process are found to be of practical interest as an additional criterion for destruction. The simulation results agree with the experimental data known from the available literature.

2D and 3D numerical models to evaluate trabecular bone damage

The comprehension of trabecular bone damage processes could be a crucial hint for understanding how bone damage starts and propagates. Currently, different approaches to bone damage identification could be followed. Clinical approaches start from dual X-ray absorptiometry (DXA) technique that can evaluate bone mineral density (BMD), an indirect indicator of fracture risk. DXA is, in fact, a two-dimensional technology, and BMD alone is not able to predict the effective risk of fractures. First attempts in overcoming this issue have been performed with finite element (FE) methods, combined with the use of three-dimensional high-resolution micro-computed tomographic images. The purpose of this work is to evaluate damage initiation and propagation in trabecular vertebral porcine samples using 2D linear-elastic FE models from DXA images and 3D linear FE models from micro-CT images. Results show that computed values of strains with 2D and 3D approaches (e.g., the minimum principal strain) are of the same order of magnitude. 2D DXA-based models still remain a powerful tool for a preliminary screening of trabecular regions that are prone to fracture, while from 3D micro-CT-based models, it is possible to reach details that permit the localization of the most strained trabecula.

Loss of longitudinal superiority marks the microarchitecture deterioration of osteoporotic cancellous bones

Osteoporosis (OP), a skeletal disease making bone mechanically deteriorate and easily fracture, is a global public health issue due to its high prevalence. It has been well recognized that besides bone loss, microarchitecture degradation plays a crucial role in the mechanical deterioration of OP bones, but the specific role of microarchitecture in OP has not been well clarified and quantified from mechanics perspective. Here, we successfully decoupled and identified the specific roles of microarchitecture, bone mass and tissue property in the failure properties of cancellous bones, through μCT-based digital modeling and finite element method simulations on bone samples from healthy and ovariectomy-induced osteoporotic mice. The results show that the microarchitecture of healthy bones exhibits longitudinal superiority in mechanical properties such as the effective stiffness, strength and toughness, which fits them well to bearing loads along their longitudinal direction. OP does not only reduce bone mass but also impair the microarchitecture topology. The former is mainly responsible for the mechanical degradation of bones in magnitude, wherever the latter accounts for the breakdown of their function-favorable anisotropy, the longitudinal superiority. Hence, we identified the microarchitecture-deterioration-induced directional mismatch between material and loading as a hazardous feature of OP and defined a longitudinal superiority index as measurement of the health status of bone microarchitecture. These findings provide useful insights and guidelines for OP diagnosis and treat assessment.

Mesh-independent damage model for trabecular bone fracture simulation and experimental validation

We propose in this study a two-dimensional constitutive model for trabecular bone combining continuum damage with embedded strong discontinuity. The model is capable of describing the three failure phases of trabecular bone tissue which is considered herein as a quasi-brittle material with strains and rotations assumed to be small and without viscous, thermal or other non-mechanical effects. The finite element implementation of the present model uses constant strain triangle (CST) elements. The displacement jump vector is implicitly solved through a return mapping algorithm at the local (finite element) level, while the global equilibrium equations are dealt with by Newton-Raphson method. The performance, accuracy and applicability of the proposed model for trabecular bone fracture are evaluated and validated against experimental measurements. These comparisons include both global and local aspects through numerical simulations of three-point bending tests performed on 10 single bovine trabeculae in the quasi-static regime.

Prediction of mechanical properties of trabecular bone in patients with type 2 diabetes using damage based finite element method

Type-2 diabetic (T2D) and osteoporosis (OP) suffered patients are more prone to fragile fracture though the nature of alteration in areal bone mineral density (aBMD) in these two cases are completely different. Therefore, it becomes crucial to compare the effect of T2D and OP on alteration in mechanical and structural properties of femoral trabecular bone. This study investigated the effect of T2D, OP, and osteopenia on bone structural and mechanical properties using micro-CT, nanoindentation and compression test. Further, a nanoscale finite element model (FEM) was developed to predict the cause of alteration in mechanical properties. Finally, a damage-based FEM was proposed to predict the pathological related alteration of bone's mechanical response. The obtained results demonstrated that the T2D group had lower volume fraction (-18.25%, p = 0.023), young's modulus (-23.47%, p = 0.124), apparent modulus (-37.15%, p = 0.02), and toughness (-40%, p = 0.001) than the osteoporosis group. The damage-based FE results were found in good agreement with the compression experiment results for all three pathological conditions. Also, nanoscale FEM results demonstrated that the elastic and failure properties of mineralised collagen fibril decreases with increase in crystal size. This study reveals that T2D patients are more prone to fragile fracture in comparison to OP and osteopenia patients. Also, the proposed damage-based FEM can help to predict the risk of fragility fracture for different pathological conditions.

Numerical analysis of hip fracture due to a sideways fall

The primary purpose of this paper is to outline a methodology for evaluating the likelihood of cortical bone fracture in the proximal femur in the event of a sideways fall. The approach includes conducting finite element (FE) analysis in which the cortical bone is treated as an anisotropic material, and the admissibility of the stress field is validated both in tension and compression regime. In assessing the onset of fracture, two methodologies are used, namely the Critical Plane approach and the Microstructure Tensor approach. The former is employed in the tension regime, while the latter governs the conditions at failure in compression. The propagation of localized damage is modeled using a constitutive law with embedded discontinuity (CLED). In this approach, the localized deformation is described by a homogenization procedure in which the average properties of cortical tissue intercepted by a macrocrack are established. The key material properties governing the conditions at failure are specified from a series of independent material tests conducted on cortical bone samples tested at different orientations relative to the loading direction. The numerical analysis deals with simulations of experiments involving the sideways fall, and the results are compared with the experimental data. This includes both the evolution of fracture pattern and the local load-displacement characteristics. The proposed approach is numerically efficient, and the results do not display a pathological mesh-dependency. Also, in contrast to the XFEM approach, the analysis does not require any extra degrees of freedom.

Image-based finite-element modeling of the human femur

Fracture is considered a critical clinical endpoint in skeletal pathologies including osteoporosis and bone metastases. However, current clinical guidelines are limited with respect to identifying cases at high risk of fracture, as they do not account for many mechanical determinants that contribute to bone fracture. Improving fracture risk assessment is an important area of research with clear clinical relevance. Patient-specific numerical musculoskeletal models generated from diagnostic images are widely used in biomechanics research and may provide the foundation for clinical tools used to quantify fracture risk. However, prior to clinical translation, in vitro validation of predictions generated from such numerical models is necessary. Despite adopting radically different models, in vitro validation of image-based finite element (FE) models of the proximal femur (predicting strains and failure loads) have shown very similar, encouraging levels of accuracy. The accuracy of such in vitro models has motivated their application to clinical studies of osteoporotic and metastatic fractures. Such models have demonstrated promising but heterogeneous results, which may be explained by the lack of a uniform strategy with respect to FE modeling of the human femur. This review aims to critically discuss the state of the art of image-based femoral FE modeling strategies, highlighting principal features and differences among current approaches. Quantitative results are also reported with respect to the level of accuracy achieved from in vitro evaluations and clinical applications and are used to motivate the adoption of a standardized approach/workflow for image-based FE modeling of the femur.

Prediction of failure in cancellous bone using extended finite element method

The objective of our study is to develop extended finite element method models of cancellous bone specimens that are capable of accurately predicting the onset and propagation of cracks under mechanical loading. In order to do so, previously published three-point bending test results of a single trabecula were replicated using two different extended finite element method approaches, namely, elastic-plastic-fracture and elastic-fracture that considered different configurations of the elasto-plastic properties of bone from which the best approach to fit the experimental data was identified. The behavior of a single trabecula was then used in 2D extended finite element method models to quantify the strength of the trabecular tissue of the forearm along three perpendicular anatomical axes. The results revealed that the elastic-plastic-fracture model better represented the experimental data in the model of a single trabecula. Considering the 2D trabecular specimens, the elastic fracture model predicted higher strength than the elastic-plastic-fracture model and there was no difference in stiffness between the two models. In general, the specimens exhibited higher failure strain and more ductile behavior in compression than in tension. In addition, strength and stiffness were found to be higher in tension than compression on average. It can be concluded that with proper parameters, extended finite element method is capable of simulating the ductile behavior of cancellous bone. The models are able to quantify the tensile strength of trabecular tissue in the various anatomical directions reporting an increased strength in the longitudinal direction of forearm cancellous bone tissue. Extended finite element method of cancellous bone proves to be a valuable tool to predict the mechanical characteristics of cancellous bones as a function of the microstructure.

Effect of orthopedic implants on canine long bone compression stiffness: a combined experimental and computational approach

Osteosynthesis for canine long bones is a complex process requiring knowledge of biology, surgical techniques and (bio)mechanical principles. Subject-specific finite element analysis constitutes a promising tool to evaluate the effect of surgical intervention on the global properties of a bone-implant construct, but suffers from a lack of validation. In this study, the biomechanical behavior of 10 canine humeri was compared before and after creation of a 10 mm bone defect stabilized with an eight-hole locking compression plate (Synthes^®) and two locking screws on each fragment. The response under compression of both intact and plated samples was measured experimentally and reproduced with a finite element model. The experimental stiffness ratio between plated and intact bone was equal to 0.39 ± 0.06. A subject-specific finite element analysis including density-dependent elasto-plastic material properties for canine bone and automatic generation of orthopedic implants was then conducted to recover these experimental results. The stiffness of intact and plated samples could be predicted, with no significant differences with experimental data. The simulated stiffness ratio between plated and intact canine bone was equal to 0.43 ± 0.03. This study constitutes a first step toward the building of a virtual database of pre-computed cases, aiming at helping the veterinary surgeons to make decisions regarding the most suited orthopedic solution for a given dog and a given fracture.

A novel phase field method for modeling the fracture of long bones

A proximal humerus fracture is an injury to the shoulder joint that necessitates medical attention. While it is one of the most common fracture injuries impacting the elder community and those who suffer from traumatic falls or forceful collisions, there are almost no validated computational methods that can accurately model these fractures. This could be due to the complex, inhomogeneous bone microstructure, complex geometries, and the limitations of current fracture mechanics methods. In this paper, we develop a novel phase field method to investigate the proximal humerus fracture. To model the fracture in the inhomogeneous domain, we propose a power-law relationship between bone mineral density and critical energy release rate. The method is validated by an in vitro experiment, in which a human humerus is constrained on both ends while subjected to compressive loads on its head, in the longitudinal direction, that lead to fracture at the anatomical neck. CT scans are employed to acquire the bone geometry and material parameters, from which detailed finite element meshes with inhomogeneous Young modulus distributions are generated. The numerical method, implemented in a high performance computing environment, is used to quantitatively predict the complex 3D brittle fracture of the bone and is shown to be in good agreement with experimental observations. Furthermore, our findings show that the damage is initiated in the trabecular bone-head and propagates outward towards the bone cortex. We conclude that the proposed phase field method is a promising approach to model bone fracture.

A round-robin finite element analysis of human femur mechanics between seven participating laboratories with experimental validation

Finite element analysis is a common tool that has been used for the past few decades to predict the mechanical behavior of bone. However, to our knowledge, there are no round-robin finite element analyses of long human bones with more than two participating biomechanics laboratories published yet, where the results of the experimental tests were not known in advance. We prepared a fresh-frozen human femur for a compression test in a universal testing machine measuring the strains at 10 bone locations as well as the deformation of the bone in terms of the displacement of the loading point at a load of 2 kN. The computed tomography data of the bone with a calibration phantom as well as the orientation of the bone in the testing machine with the according boundary conditions were delivered to seven participating laboratories. These were asked to perform a finite element analysis simulating the experimental setup and deliver their results to the coordinator without knowing the experimental results. Resultantly, four laboratories had deviations from the experimentally measured strains of less than 40%, and three laboratories had deviations of their numerically determined values compared to the experimental data of more than 120%. These deviations are thought to be based on different material laws and material data, as well as different material mapping methods. Investigations will be conducted to clarify and assess the reasons for the large deviations in the numerical data. It was shown that the precision of finite element models of the human femur is not yet as developed as desired by the biomechanics community.

Comparative Biomechanical Analysis of Stress-Strain State of the Elbow Joint After Displaced Radial Head Fractures

Radial head fractures are becoming a major public health problem and are an increasingly important target for both clinical and mechanical researchers. In this work, comparative biomechanical analyses of the stress–strain state of a healthy elbow joint and elbow joints with radial head compression from 2 to 5 mm due to injury are performed. Three-dimensional models of the elbow joint with cartilage surfaces and ligaments were constructed based on the results of computed tomography. This study is focused on an elbow joint range of motion ranging from 0° to 120° flexion. Analysis of the stress–strain state of cartilage and ligaments under the influence of functional loads is conducted using a finite element method (FEM) and the ABAQUS software package. The results show that with increasing compression of the radial head, contact stress increases at the olecranon, which can lead to cartilage damage. Analysis of displacement shows that compression of the radial head during full extension of the elbow joint leads to an increased humeral shift from 1.14° ± 0.22 in the healthy joint to 10.3° ± 2.13 during 5-mm compression of the radial head. Mathematical modeling performed in this study proved that reducing the height of the radial head and the contact area between the radial head and the humeral head led to increased medial collateral ligament stresses of up to 36 ± 3.8 MPa. This work confirmed that the head of the radius is the main stabilizing structure of the elbow joint and that the medial collateral ligament is the second structure responsible for valgus stability of the elbow joint.

Validation of a parametric finite element human femur model

OBJECTIVE

Finite element (FE) models with geometry and material properties that are parametric with subject descriptors, such as age and body shape/size, are being developed to incorporate population variability into crash simulations. However, the validation methods currently being used with these parametric models do not assess whether model predictions are reasonable in the space over which the model is intended to be used. This study presents a parametric model of the femur and applies a unique validation paradigm to this parametric femur model that characterizes whether model predictions reproduce experimentally observed trends.

METHODS

FE models of male and female femurs with geometries that are parametric with age, femur length, and body mass index (BMI) were developed based on existing statistical models that predict femur geometry. These parametric FE femur models were validated by comparing responses from combined loading tests of femoral shafts to simulation results from FE models of the corresponding femoral shafts whose geometry was predicted using the associated age, femur length, and BMI. The effects of subject variables on model responses were also compared with trends in the experimental data set by fitting similarly parameterized statistical models to both the results of the experimental data and the corresponding FE model results and then comparing fitted model coefficients for the experimental and predicted data sets.

RESULTS

The average error in impact force at experimental failure for the parametric models was 5%. The coefficients of a statistical model fit to simulation data were within one standard error of the coefficients of a similarly parameterized model of the experimental data except for the age parameter, likely because material properties used in simulations were not varied with specimen age. In simulations to explore the effects of femur length, BMI, and age on impact response, only BMI significantly affected response for both men and women, with increasing BMI producing higher impact forces.

CONCLUSIONS

Impactor forces from simulations, on average, matched experimental values at the time of failure. In addition, the simulations were able to match the trends in the experimental data set.

Hip joint geometry effects on cartilage contact stresses during a gait cycle

The cartilage surface geometry of natural human hip joint is commonly regarded as sphere. It has been widely applied in computational simulation and hip joint prosthesis design. Some new geometry models have been developed and the sphere assumption has been questioned recently. The objective of this study was to analyze joint geometry effects on cartilage contact stress distribution and investigate contact patterns during a whole gait cycle. Hip surface was reconstructed from CT data of a healthy volunteer. Three finite element (FE) models of hip joint were developed from different cartilage geometries: natural geometry, sphere and rotational ellipsoid. Loads at ten instants of gait cycle were applied to these models based on published in-vivo data. FE predictions of peak contact pressure during gait of natural hip were compared with sphere and rotational ellipsoid replaced hip joint. Contact occurs mainly in upper anterior region of both acetabulum and femur distributing along sagittal plane of human body. It moves towards inferolateral aspect as the resultant joint reaction force changes during walking for natural hip. Peak pressures at the instant with maximum contact force were 7.48 MPa, 14.97 MPa and 13.12 MPa for models with natural hip surface, sphere replaced and rotational ellipsoid replaced surface respectively. During the whole gait cycle, contact pressure of natural hip ranked lowest in most of the instants, followed by rotational ellipsoid replaced and sphere replaced hip. The results indicate that rotational ellipsoid is more consistent with natural hip cartilage geometry than sphere during normal walking. This means rotational ellipsoid prosthesis could give a better description of physiological structure compared with standard sphere prosthesis. Therefore, rotational ellipsoid would be a better choice for prosthesis design.

Assessment of compressive failure process of cortical bone materials using damage-based model

The main failure factors of cortical bone are aging or osteoporosis, accident and high energy trauma or physiological activities. However, the mechanism of damage evolution coupled with yield criterion is considered as one of the unclear subjects in failure analysis of cortical bone materials. Therefore, this study attempts to assess the structural response and progressive failure process of cortical bone using a brittle damaged plasticity model. For this reason, several compressive tests are performed on cortical bone specimens made of bovine femur, in order to obtain the structural response and mechanical properties of the material. Complementary finite element (FE) model of the sample and test is prepared to simulate the elastic-to-damage behavior of the cortical bone using the brittle damaged plasticity model. The FE model is validated in a comparative method using the predicted and measured structural response as load-compressive displacement through simulation and experiment. FE results indicated that the compressive damage initiated and propagated at central region where maximum equivalent plastic strain is computed, which coincided with the degradation of structural compressive stiffness followed by a vast amount of strain energy dissipation. The parameter of compressive damage rate, which is a function dependent on damage parameter and the plastic strain is examined for different rates. Results show that considering a similar rate to the initial slope of the damage parameter in the experiment would give a better sense for prediction of compressive failure.

Prediction of the mechanical response of canine humerus to three-point bending using subject-specific finite element modelling

Subject-specific finite element models could improve decision making in canine long-bone fracture repair. However, it preliminary requires that finite element models predicting the mechanical response of canine long bone are proposed and validated. We present here a combined experimental–numerical approach to test the ability of subject-specific finite element models to predict the bending response of seven pairs of canine humeri directly from medical images. Our results show that bending stiffness and yield load are predicted with a mean absolute error of 10.1% (±5.2%) for the 14 samples. This study constitutes a basis for the forthcoming optimization of canine long-bone fracture repair.

Modelling of bone fracture and strength at different length scales: a review

In this paper, we review analytical and computational models of bone fracture and strength. Bone fracture is a complex phenomenon due to the composite, inhomogeneous and hierarchical structure of bone. First, we briefly summarize the hierarchical structure of bone, spanning from the nanoscale, sub-microscale, microscale, mesoscale to the macroscale, and discuss experimental observations on failure mechanisms in bone at these scales. Then, we highlight representative analytical and computational models of bone fracture and strength at different length scales and discuss the main findings in the context of experiments. We conclude by summarizing the challenges in modelling of bone fracture and strength and list open topics for scientific exploration. Modelling of bone, accounting for different scales, provides new and needed insights into the fracture and strength of bone, which, in turn, can lead to improved diagnostic tools and treatments of bone diseases such as osteoporosis.

Linear Elastic Analysis of Bovine Cortical Bone under Compression Loading

Linear elastic response of the bovine cortical bone has been examined under compression load. Experimental and computational methods were used to observe and predict the response of cortical bone. In computational method, two mechanical behaviors of isotropic and orthotropic were considered to simulate the cortical bone deformation. In experimental process, the specimens were designed to show maximum stiffness and strength by specifying osteon direction along loading axis during tests. The tests were controlled by displacement rate of 0.5 mm/minute and the overall stiffness responses of the structures were recorded to extract mechanical properties and also for validation aims. Finite Element Method (FEM) was used to model the linear response of the structure by using ABAQUS6.9EF. The FE results using orthotropic definition shows a good correlation with experimental data. A discussion was given based on overall stiffness and effective stress variation for both mechanical behaviors. In order to design the optimal implant structure, the presented study was proposed for prediction of bone structure deformation that attached to the orthopedic implants.

Finite element prediction with experimental validation of damage distribution in single trabeculae during three-point bending tests

There is growing evidence that information on trabecular microarchitecture can improve the assessment of fracture risk. One current strategy is to exploit finite element (FE) analysis applied to experimental data of mechanically loaded single trabecular bone tissue obtained from non-invasive imaging techniques for the investigation of the damage initiation and growth of bone tissue. FE analysis of this type of bone has mainly focused on linear and non-linear analysis to evaluate the bone's failure properties. However, there is a lack of experimentally validated FE damage models at trabecular bone tissue level allowing for the simulation of the progressive damage process (initiation and growth) till complete fracture. Such models are needed to perform enhanced prediction of the apparent failure mechanical properties needed to assess the fracture risk of bone organs. In the current study, we develop a FE model based on a continuum damage mechanics (CDM) approach to simulate the damage initiation and propagation of a single trabecula till complete facture in quasi-static regime. Three-point bending experiments were performed on single bovine trabeculae and compared to FE results. In order to validate the proposed FE mode, (i) the force displacement curve was compared to the experimental one and (ii) the damage distribution was correlated to the measured one obtained by digital image correlation based on stress whitening in bone, reported to be correlated to microdamage. A very good agreement was obtained between the FE and experimental results, indicating that the proposed damage investigation protocol based on FE analysis and testing is reliable to assess the damage behavior of bone tissue and that the current damage model is able to accurately simulate the damaging and fracturing process of single trabeculae under quasi static load.

Mechanical competence of bone: a new parameter to grade trabecular bone fragility from tortuosity and elasticity

With the elderly population increase, osteoporosis and its consequences have become not just a health issue but also a serious economic burden. The trabecular bone structure plays a very important role for the bone quality and mechanical competence of the scaffold. Currently, it is claimed that the trabecular microarchitecture understanding can improve the fracture risk prediction above 65%. Several parameters seem to be correlated providing structural details of the trabecular bone network. However, the tortuosity of the trabeculae has not yet been systematically taken into account and its contribution has not been fully investigated and understood. In this paper, we discuss the relationship between the trabecular tortuosity, connectivity, volume fraction, and elasticity, and provide a unified parameter to estimate the mechanical competence of the structure. It is shown that the trabecular network tortuosity presents high linear correlation with the other parameters and that the trabeculae tend to get aligned in the direction where the structure is mostly submitted to stress, corresponding to higher stiffness orientation. This new parameter will help to integrate the relevant information of bone microarchitecture quality and assess more directly the real trabecular fragility in osteoporotic patients.