Development and validation of a subject-specific finite element model of a human clavicle

This study aimed to develop and validate a finite element (FE) model of a human clavicle which can predict the structural response and bone fractures under both axial compression and anterior–posterior three-point bending loads. Quasi-static non-injurious axial compression and three-point bending tests were first conducted on a male clavicle followed by a dynamic three-point bending test to fracture. Then, two types of FE models of the clavicle were developed using bone material properties which were set to vary with the computed tomography image density of the bone. A volumetric solid FE model comprised solely of hexahedral elements was first developed. A solid-shell FE model was then created which modelled the trabecular bone as hexahedral elements and the cortical bone as quadrilateral shell elements. Finally, simulations were carried out using these models to evaluate the influence of variations in cortical thickness, mesh density, bone material properties and modelling approach on the biomechanical responses of the clavicle, compared with experimental data. The FE results indicate that the inclusion of density-based bone material properties can provide a more accurate reproduction of the force–displacement response and bone fracture timing than a model with uniform bone material properties. Inclusion of a variable cortical thickness distribution also slightly improves the ability of the model to predict the experimental response. The methods developed in this study will be useful for creating subject-specific FE models to better understand the biomechanics and injury mechanism of the clavicle.     

A meshless microscale bone tissue trabecular remodelling analysis considering a new anisotropic bone tissue material law

In this work, a novel anisotropic material law for the mechanical behaviour of the bone tissue is proposed. This new law, based on experimental data, permits to correlate the bone apparent density with the obtained level of stress. Combined with the proposed material law, a biomechanical model for predicting bone density distribution was developed, based on the assumption that the bone structure is a gradually self-optimising anisotropic biological material that maximises its own structural stiffness. The strain and the stress field required in the iterative remodelling process are obtained by means of an accurate meshless method, the Natural Neighbour Radial Point Interpolation Method (NNRPIM). Comparing with other numerical approaches, the inclusion of the NNRPIM presents numerous advantages such as the high accuracy and the smoother stress and strain field distribution. The natural neighbour concept permits to impose organically the nodal connectivity and facilitates the analysis of convex boundaries and extremely irregular meshes. The viability and efficiency of the model were tested on several trabecular benchmark patch examples. The results show that the pattern of the local bone apparent density distribution and the anisotropic bone behaviour predicted by the model for the microscale analysis are in good agreement with the expected structural architecture and bone apparent density distribution.     

Anisotropic post-yield response of cancellous bone simulated by stress–strain curves of bulk equivalent structures

During the last decade, finite element (FE) modelling has become ubiquitous in understanding complex mechanobiological phenomena, e.g. bone–implant interactions. The extensive computational effort required to achieve biorealistic results when modelling the post-yield behaviour of microstructures like cancellous bone is a major limitation of these techniques. This study describes the anisotropic biomechanical response of cancellous bone through stress–strain curves of equivalent bulk geometries. A cancellous bone segment, reverse engineered by micro computed tomography, was subjected to uniaxial compression. The material's constitutive law, obtained by nano-indentations, was considered during the simulation of the experimental process. A homodimensionally bulk geometry was employed to determine equivalent properties, resulting in a similar anisotropic response to the trabecular structure. The experimental verification of our model sustained that the obtained stress–strain curves can adequately reflect the post-yield behaviour of the sample. The introduced approach facilitates the consideration of nonlinearity and anisotropy of the tissue, while reducing the geometrical complexity of the model to a minimum.     

Determination of orthotropic bone elastic constants using FEA and modal analysis

Finite element models have been widely employed in an effort to quantify the stress and strain distribution around implanted prostheses and to explore the influence of these distributions on their long-term stability. In order to provide meaningful predictions, such models must contain an appropriate reflection of mechanical properties. Detailed geometrical and density information is now readily available from CT scanning. However, despite the use of phantoms, a method of determining mechanical properties (or elastic constants) from bone density has yet to be made available in a usable form.In this study, a cadaveric bone was CT scanned and its natural frequencies were measured using modal analysis. Using the geometry obtained from the CT scan data, a finite element mesh was created with the distribution of density established by matching the mass of the FE bone model with the mass of the cadaveric bone. The maximum values of the orthotropic elastic constants were then established by matching the predictions from FE modal analyses to the experimental natural frequencies, giving a maximum error of 7.8% over 4 modes of vibration. Finally, the elastic constants of the bone derived from the analyses were compared with those measured using ultrasound techniques. This produced a difference of <1% for both the maximum density and axial Young's Modulus. This study has thereby produced an orthotropic finite element model of a human femur. More importantly, however, is the implication that it is possible to create a valid FE model by simply comparing the FE results with the measured resonant frequency of the CT scanned bone.     

Influence of material properties on the mechanical behaviour of the L5-S1 intervertebral disc in compression: a nonlinear finite element study

A simple finite element model of the L₅-S₁ intervertebral disc body has been constructed; it is circular and symmetrical about the sagittal plane. The annulus fibrosus of the model was idealized as an inhomogeneous composite of an isotropic ground substance, reinforced by helically oriented collagen fibres so that the model has six different structural components namely: cortical bone, cancellous bone, cartilaginous endplates, nucleus pulposus, ground substance and collagen fibres. A sensitivity analysis of the material properties of each structural component was carried out by varying those properties for one structural component at a time and evaluating the changes in the biomechanical response to compressive displacements. Experimentally available relations between the applied compressive force and the vertical displacements, the nucleus pulposus pressure increase and the disc lateral bulge were used to evaluate the biomechanical responses for each set of material properties. Results showed that both the Poisson's ratio and the Young's modulus of the ground substance play an important role in the prediction of the biomechanical response.     

Homogenization theory and digital imaging: A basis for studying the mechanics and design principles of bone tissue

Bone tissue is a complex multilevel composite which has the ability to sense ad respond to its mechanical environment. It is believed that bone cells called osteocytes within the bone matrix sense the mechanical environment and determine whether structural alterations are needed. At present it is not known, however, how loads are transferred from the whole bone level to cells. A computational procedure combining representative volume element (RVE) based homogenization theory with digital imaging is proposed to estimate strains at various levels of bone structure. Bone tissue structural organization and RVE based analysis are briefly reviewed. The digital image based computational procedure was applied to estimate strains in individual trabeculae (first-level microstructure). Homogenization analysis of an idealized model was used to estimate strains at one level of bone structure around osteocyte lacunae (second-level trabecular microstructure). The results showed that strain at one level of bone structure is amplified to a broad range at the next microstructural level. In one case, a zeor-level tensile principal strain of 495 muE engendered strains ranging between -1000 and 7000 muE in individual trabeculae (first-level microstructure). Subsequently, a first-level tensile principal strains of 1325 muE within an inidividual trabecula engendered strains ranging between 782 and 2530 muE around osteocyte lacunae. Lacunar orientation was found to influence strains around osteocyte lacunae much more than lacunar ellipticity. In conclusion, the computational procedure combining homogenization theory with digital imaging can proveide estimates of cell level strains within whole bones. Such results may be used to bridge experimental studies of bone adaptation at the whole bone and cell culture level. (c) 1994 John Wiley & Sons, Inc.     

A Three-Dimensional Finite Element Scheme to Investigate the Apparent Mechanical Properties of Trabecular Bone

A computational technique is described for investigating the apparent mechanical properties of trabecular bone based on tissue geometry obtained from the marching cubes volume rendering scheme. Using this scheme, a 3D representation of the trabecular bone was extracted from two-dimensional cross-sections of the tissue originating from a quantitative serial sectioning procedure. Surface information consists of node coordinates and polygon connectivity in a 3D space. A custom, adaptive mesh generation technique using a normal offset was used to prepare 3D finite element volume meshes (4-node tetrahedral elements) of variable mesh density from the extracted surface geometry. Nine target mesh resolutions (32 μm to 107 μm) were examined for a (1.5 mmx 1.5 mmx 2 mm) volume of trabecular bone. A mesh density of 50,000 elements/mm(3) of bone tissue was found to be adequate for convergence of apparent (bulk) modulus for 1% uniaxial compression. For this convergent case, the maximum local normal compressive tissue stress was 400 MPa which was six hundred-fold greater than the computed apparent stress. Variation in the apparent modulus was less than 5% when Poisson's ratio values were varied between 0.1 and 0.4. Poisson's ratio values greater than 0.4 had a more marked effect on the apparent modulus. Based upon these results, approximately 1 million, 4-node tetrahedral elements are required to analyze a continuum scale model of trabecular bone (5 mm cube).     

Development of a three-dimensional finite element model of a human tibia using experimental modal analysis.

The modal analysis of a human tibia consisted of characterizing its dynamic behavior by determining natural frequency, damping ratio and mode shapes. Two methods were used to perform the modal analysis: (1) a finite element method (structural model); (2) an experimental modal analysis (modal model). The experimental modal model was used to optimize the structural model. After optimization, differences in results between the two models were found to be due only to mechanical properties and mass distribution. The influences of boundary conditions and geometric properties (such as inertia and length) were eliminated by the finite element model itself. The percent relative error between the two methods was approximately 3%, corresponding to the standard deviation of the measured frequencies. For the frequency range considered, the mode shapes were bending modes in two different vibration planes (latero-medial and sagittal), with a slight torsion effect due to the twisted geometry of the tibia.     

Comparative analysis of bone remodelling models with respect to computerised tomography-based finite element models of bone

Subject-specific finite element models are an extensively used tool for the numerical analysis of the biomechanical behaviour of human bones. However, bone modelling is not an easy task due to the complex behaviour of bone tissue, involving non-homogeneous and anisotropic mechanical properties. Moreover, bone is a living tissue and therefore its microstructure and mechanical properties evolve with time in a known process called bone remodelling. This phenomenon has been widely studied, many being the numerical models that have been formulated to predict density distribution and its evolution in several bones. The aim of the present study is to assess the capability of a bone remodelling model to predict the bone density distribution of different types of human bone (femur, tibia and mandible) comparing the obtained results with the bone density estimated by means of computerised tomography. Good accuracy was observed for the bone remodelling predictions including the thickness of the cortical layer.     

Prediction of shape and internal structure of the proximal femur using a modified level set method for structural topology optimisation

A computational framework was developed to simulate the bone remodelling process as a structural topology optimisation problem. The mathematical formulation of the Level Set technique was extended and then implemented into a coronal plane model of the proximal femur to simulate the remodelling of internal structure and external geometry of bone into the optimal state. Results indicated that the proposed approach could reasonably mimic the major geometrical and material features of the natural bone. Simulation of the internal bone remodelling on the typical gross shape of the proximal femur, resulted in a density distribution pattern with good consistency with that of the natural bone. When both external and internal bone remodelling were simulated simultaneously, the initial rectangular design domain with a regularly distributed mass reduced gradually to an optimal state with external shape and internal structure similar to those of the natural proximal femur.     

A new cortical thickness mapping method with application to an in vivo finite element model

Finite element modelling of musculoskeletal systems, with geometrical structures constructed from computed tomography (CT) scans, is a useful and powerful tool for biomechanical studies. The use of CT scans from living human subjects, however, is still limited. Accurate reconstruction of thin cortical bone structures from CT scans of living human subjects is especially problematic, due to low CT resolution that results from mandatory low radiation doses and/or involuntary movements of the subject. In this study, a new method for mapping cortical thickness is described. Using the method, cortical thickness measurements of a coxal (pelvis) bone obtained from CT scans of a cadaver were mapped to the coxal geometry as obtained through CT scans of a live human subject, resulting in accurate cortical thickness while maintaining geometric fidelity of the live subject. The mapping procedure includes shape-preserving parameterisation, mesh movement and interpolation of thickness using a search algorithm. The methodology is applicable to modelling of other bones where accurate cortical thickness is needed and for which such data exist.     

Biomaterial-loaded allograft threaded cage for the treatment of femoral head osteonecrosis in a goat model

The purpose of this study was to evaluate the efficacy of core decompression with a biomaterial-loaded allograft threaded cage (ATC) for the treatment of femoral head osteonecrosis in an established goat model. First, bilateral early-stage osteonecrosis was induced. After core decompression, the remaining goats were randomly divided into three groups: Group A, the goats were left without any treatment; Group B, the goats were treated with implanting a composite of autologous bone and decalcified bone matrix (DBM); Group C, the goats were treated using insertion of ATC loaded with DBM and autogenous bone graft. Then radiographic, histological and biomechanical analysis were taken in each group at 5, 10, and 20 weeks postoperation. In Group A, the classical signs of osteonecrosis of the femoral head were identified 10 weeks after the induction. Twenty weeks later, the density, surface and biomechanical stability of the femoral head were normal in Group C, while an irregular surface and an inhomogeneous microstructure or variation of density/hardness were identified in Group B. The specimens revealed a continuous trabaecular bone structure throughout the cage and extensive bone ingrowth and remodeling in Group C, while fibrous tissue was evident in Group B. Core decompression with a biomaterial loaded ATC almost uniformly delays or arrests the progression of the disease before articular collapse, and it could help to get the balance between bone resorption and new bone formation, strengthen structural mechanics of the femoral head, provide structure support of articular cartilage.     

Factors affecting subject-specific finite element models of implant-fitted rat bone specimens: critical analysis of a technical protocol

The authors propose a protocol to derive finite element (FE) models from micro computer tomography scans of implanted rat bone. A semi-automatic procedure allows segmenting the images using specimen-specific bone mineral density (BMD) thresholds. An open-source FE model generator processes the segmented images to a quality tetrahedral mesh. The material properties assigned to each element are integrated from the BMD field. Piecewise, threshold-dependent density–elasticity relationships are implemented to limit the effects of metal artefacts. A detailed sensitivity study highlights the coherence of the generated models and quantifies the influence of the modelling parameters on the results. Two applications of the protocol are proposed. The stiffness of bare and implanted rat tibiae specimens is predicted by simulating three-point bending and inter-implant displacement, respectively. Results are compared with experimental tests. The mean value and the variability between the specimens are well captured in both tests.     

A contact model with ingrowth control for bone remodelling around cementless stems.

This work presents a computational model for bone remodelling around cementless stems. The problem is formulated as a material optimisation problem considering the bone and stem surfaces to be in contact. To emphasise the behaviour of the bone/stem interface, the computer model detects the existence of bone ingrowth during the remodelling; consequently, the contact conditions are changed for a better interface simulation. The trabecular bone is modelled as a strictly orthotropic material with equivalent properties computed by homogenisation. The distribution of bone relative density is obtained by the minimisation of a function that considers both the bone structural stiffness and the biological cost associated with metabolic maintenance of bone tissue. The situation of multiple load conditions is considered. The remodelling law, obtained from the necessary conditions for an optimum, is derived analytically from the optimisation problem and solved numerically using a suitable finite element mesh. The formulation is applied to an implanted femur. Results of bone density and ingrowth distribution are obtained for different coating conditions. Bone ingrowth does not occur over the entire coated surfaces. Indeed, we observed regions where separation or high relative displacement occurs that preclude bone ingrowth attachment. This prediction of the model is consistent with clinical observations of bone ingrowth. Thus, this model, which detect bone ingrowth and allow modification of the interface conditions, are useful for analysis of existing stems as well as design optimisation of coating extent and location on such stems.     

Computer aided stress analysis of long bones utilizing computed tomography

A computer aided analysis method has been developed which utilizes computed tomography (CT) and a finite element (FE) computer program to determine the stress-displacement pattern in a long bone section. The CT data file provides the geometry, the apparent density and the elastic properties for the three-dimensional FE model. A developed pre-processor generates the FE model of a human diaphyseal tibia section which is then analyzed by the SAP IV finite element program. The results obtained are sorted and displayed by a developed post-processor and compared with stresses and deformations from the literature. The model generation method was verified by applying it to a model of simple geometry and boundary conditions, then comparing the results with the analytical solution of the same problem. The convergence behavior of nodal displacements was tested as a function of mesh refinement. This method provides an automatic, versatile, non-invasive and accurate tool of long bone modeling for finite element stress analysis.     

Finite element model development of a child pelvis with optimization-based material identification

A finite element (FE) model of a 10-years-old child pelvis was developed and validated against experimental data from lateral impacts of pediatric pelves. The pelvic bone geometry was reconstructed from a set of computed tomography images, and a hexahedral mesh was generated using a new octree-based hexahedral meshing technique. Lateral impacts to the greater trochanter and iliac wing of the seated pelvis were simulated. Sensitivity analysis was conducted to identify material parameters that substantially affected the model response. An optimization-based material identification method was developed to obtain the most favorable material property set by minimizing differences in biomechanical responses between experimental and simulation results. This study represents a pilot effort in the development and validation of age-dependent musculoskeletal FE models for children, which may ultimately serve to evaluate injury mechanisms and means of protection for the pediatric population.     

Experimental validation of a finite element model of the proximal femur using digital image correlation and a composite bone model

Computational biomechanical models are useful tools for supporting orthopedic implant design and surgical decision making, but because they are a simplification of the clinical scenario they must be carefully validated to ensure that they are still representative. The goal of this study was to assess the validity of the generation process of a structural finite element model of the proximal femur employing the digital image correlation (DIC) strain measurement technique. A finite element analysis model of the proximal femur subjected to gait loading was generated from a CT scan of an analog composite femur, and its predicted mechanical behavior was compared with an experimental model. Whereas previous studies have employed strain gauging to obtain discreet point data for validation, in this study DIC was used for full field quantified comparison of the predicted and experimentally measured strains. The strain predicted by the computational model was in good agreement with experimental measurements, with R(2) correlation values from 0.83 to 0.92 between the simulation and the tests. The sensitivity and repeatability of the strain measurements were comparable to or better than values reported in the literature for other DIC tests on tissue specimens. The experimental-model correlation was in the same range as values obtained from strain gauging, but the DIC technique produced more detailed, full field data and is potentially easier to use. As such, the findings supported the validity of the model generation process, giving greater confidence in the model's predictions, and digital image correlation was demonstrated as a useful tool for the validation of biomechanical models.     

Influences of geometrical and mechanical properties of bone tissues in mandible behaviour – experimental and numerical predictions

The properties and geometry of bone in the mandible play a key role in mandible behaviour during a person’s lifetime, and attention needs to be paid to the influence of bone properties. We analysed the effect of bone geometry, size and bone properties in mandible behaviour, experimenting on cadaveric mandibles and FE models. The study was developed using the geometry of a cadaveric mandible without teeth. Three models of cadaveric condyles were experimentally tested with instrumented with four rosettes, and a condyle reaction of 300 N. Four finite element models were considered to validate the experiments and analyse mandible behaviour. One numeric model was simulated with 10 muscles in a quasi-static condition. The experimental results present different condyle stiffness’s, of 448, 215 and 254 N/mm. The values presented in the rosettes are influenced by bone geometry and bone thickness; maximum value was ?600 με in rosette #4, and the maximum strain difference between mandibles was 111%. The numerical results show that bone density decreases and strain distribution increases in the thinner mandible regions. Nevertheless, the global behaviour of the structure remains similar, but presents different strain magnitudes. The study shows the need to take into account bone characteristics and their evolutions in order to improve implant design and fixation throughout the patient life. The change in bone stiffness promotes a change in maximum strain distribution with same global behaviour.     

An improved method to assess torsional properties of rodent long bones

Torsion is an important testing modality commonly used to calculate structural properties of long bones. However, the effects of size and geometry must be excluded from the overall structural response in order to compare material properties of bones of different size, age and species. We have developed a new method to analyze torsional properties of bones using actual cross-sectional information and length-wise geometrical variations obtained by micro-computed topographic (μCT) imaging. The proposed method was first validated by manufacturing three rat femurs through rapid prototyping using a plastic with known material properties. The observed variations in calculated torsional shear modulus of the hollow elliptical model of mid-shaft cross-section (Ekeland et al.), multi-prismatic model of five true cross-sections (Levenston et al.) and multi-slice model presented in this study were 96%, ?7% and 6% from the actual properties of the plastic, respectively. Subsequently, we used this method to derive relationships expressing torsional properties of rat cortical bone as a function of μCT-based bone volume fraction or apparent density over a range of normal and pathologic bone densities. Results indicate that a regression model of shear modulus or shear strength and bone volume fraction or apparent density described at least 81% of the variation in torsional properties of normal and pathologic bones. Coupled with the structural rigidity analysis technique introduced by the authors, the relationships reported here can provide a non-invasive tool to assess fracture risk in bones affected by pathologies and/or treatment options.