Direct mechanics assessment of elastic symmetries and properties of trabecular bone architecture

A method is presented to find orthotropic elastic symmetries and constants directly from the elastic coefficients in the overall stiffness matrix of trabecular bone test specimens. Contrary to earlier developed techniques, this method does not require pure orthotropic behavior or additional fabric measurements. The method uses high-resolution computer reconstructions of trabecular bone specimens as input for large-scale FE-analyses to determine all the 21 elastic coefficients in the overall stiffness matrix of the specimen, using a direct mechanics approach. An optimization procedure is then used to find the coordinate transformation that yields the best orthotropic representation of this matrix. The method is illustrated here relative to two trabecular bone specimens. The techniques developed here can be used to obtain a complete characterization of the mechanical properties of trabecular architecture. With the development of in vivo reconstruction techniques, even in vivo measurements will be possible.     

Long-term study of bone remodelling after femoral stem: a comparison between dexa and finite element simulation

A hip replacement with a cemented or cementless femoral stem produces an effect on the bone called adaptive remodelling, attributable to mechanical and biological factors. The objective of all of cementless prostheses designs has been to achieve a perfect transfer of loads in order to avoid stress-shielding, which produces an osteopenia. In order to quantify this, the long term and mass-produced study with dual energy X-ray absorptiometry (DEXA) is necessary. Finite element (FE) simulation makes possible the explanation of the biomechanical changes which are produced in the femur after stem implantation. The good correlation obtained between the results of the FE simulation and the densitometric study allow, on one hand, to explain from the point of view of biomechanical performance the changes observed in bone density in the long-term, where it is clear that these are due to a different transfer of load in the implanted model compared to the healthy femur; on the other hand, it validates the simulation model, in a way that it can be used in different conditions and at different time periods, to carry out a sufficiently precise prediction of the evolution of the bone density from the biomechanical behaviour in the interaction between the prosthesis and femur.     

Relationship between the shape and density distribution of the femur and its natural frequencies of vibration

It has been recently suggested that mechanical loads applied at frequencies close to the natural frequencies of bone could enhance bone apposition due to the resonance phenomenon. Other applications of bone modal analysis are also suggested. For the above-mentioned applications, it is important to understand how patient-specific bone shape and density distribution influence the natural frequencies of bones. We used finite element models to study the effects of bone shape and density distribution on the natural frequencies of the femur in free boundary conditions. A statistical shape and appearance model that describes shape and density distribution independently was created, based on a training set of 27 femora. The natural frequencies were then calculated for different shape modes varied around the mean shape while keeping the mean density distribution, for different appearance modes around the mean density distribution while keeping the mean bone shape, and for the 27 training femora. Single shape or appearance modes could cause up to 15% variations in the natural frequencies with certain modes having the greatest impact. For the actual femora, shape and density distribution changed the natural frequencies by up to 38%. First appearance mode that describes the general cortical bone thickness and trabecular bone density had one of the strongest impacts. The first appearance mode could therefore provide a sensitive measure of general bone health and disease progression. Since shape and density could cause large variations in the calculated natural frequencies, patient-specific FE models are needed for accurate estimation of bone natural frequencies.     

Femoral neck bone adaptation to weight-bearing physical activity by computational analysis

Individual differences in bone mass distribution at the proximal femur may be determined by daily weight-bearing physical activity (PA) since bone self-adapts according to the mechanical loads that is submitted. The aim of this study was to analyse computationally the effect of different weight-bearing PA types in the adaptation of the femoral neck (FN) by analysing regional differences in bone mineral density (BMD) at the integral FN and its superior, inferior, anterior and posterior subregions. To achieve this, it was adopted a 3-D femoral finite element (FE) model coupled with a suitable bone remodeling model. Different PA types were determined based both on ordinary lifestyle and mechanically more demanding PA as low magnitude impacts (L–I), moderate-magnitude impacts from odd directions (O–I) and high-magnitude vertical impacts (H–I). It was observed that as time spent in weight-bearing PA increases, BMD augment around the integral FN, but with different bone mass gain rates between subregions depending on the magnitude and directions of the hip contact forces; H–I was the type of weight-bearing PA which structurally most favor the gain of bone mass superiorly at the FN while both the H–I and the O–I types of PA promoted the largest bone mass gain rates at the anterior and posterior subregions of the FN. Because these types of weight-bearing PA were associated with a more uniform bone mass spatial distribution at the FN, they should provide a potential basis for targeted PA-based intervention programs for improving hip strength.     

Thermal analysis of bone cement polymerisation at the cement-bone interface

The two major problems that have been reported with the use of polymethylmethacrylate (PMMA) cement are thermal necrosis of surrounding bone due to the high heat generation during polymerisation and chemical necrosis due to unreacted monomer release. Computer models have been used to study the temperature and monomer distribution after cementation. In most of these models, however, polymerisation is modelled as temperature independent and cancellous bone is modelled as a continuum. Such models thus cannot account for the expected important role of the trabecular bone micro-structure. The aim of this study is to investigate the distribution of temperature and monomer leftover at the cancellous bone–cement interface during polymerisation for a realistic trabecular bone—cement micro-structure and realistic temperature-dependent polymerisation kinetics behaviour.

A 3-D computer model of a piece of bovine cancellous bone that underwent pressurization with bone–cement was generated using a micro-computed tomography scanner. This geometry was used as the basis for a finite element model and a temperature-dependent problem for bone cement polymerisation kinetics was solved to simulate the bone cement polymerisation process in the vicinity of the interface. The transient temperature field throughout the interface was calculated, along with the polymerisation fraction distribution in the cement domain.

The calculations revealed that the tips of the bone trabeculae that are embedded in the cement attain temperatures much higher than the average temperature of the bone volume. A small fraction of the bone (10%) is exposed to temperatures exceeding 70°C, but the exposure time to these high temperatures is limited to 50 s. In the region near the bone, the cement polymerisation fraction (about 84%) is less than that in the centre (where it is reaching values of over 96%). An important finding of this study thus is the fact that the bone tissue that is subjected to the highest temperatures is also subjected to high leftover monomer concentration. Furthermore the maximum bone temperature is reached relatively early, when monomer content in the neighbouring cement is still quite high.     

Modelling external bone adaptation using evolutionary structural optimisation

External remodelling is significant in the bone healing process, and it is essential to predict the bone external shape in the design of artificial bone grafts. This paper demonstrates the effectiveness of the evolutionary structural optimisation (ESO) method for the simulation of bone morphology. A two-dimensional ESO strategy is developed which is capable of finding the modified bone topology beginning with any geometry under any loading conditions. The morphology of bone structure is described by the quantitative bone adaptation theory, which is integrated with the finite element method. The evolutionary topology optimisation process is introduced to find the bone shape. A rectangle, which occupies a larger space than the external shape of the bone structure, is specified as a design domain; the evolutionary process iteratively eliminates and redistributes material throughout the domain to obtain an optimum arrangement of bone materials. The technique has been tested on a wide range of examples. In this paper, the formation of trabecular bone architecture around an implant is studied; as another example, the growth of the coronal section of a vertebral body is predicted. The examples support the assertion that the external shape of bone structure can be successfully predicted by the proposed ESO procedure.     

Determining the most important cellular characteristics for fracture healing using design of experiments methods

Computational models are employed as tools to investigate possible mechanoregulation pathways for tissue differentiation and bone healing. However, current models do not account for the uncertainty in input parameters, and often include assumptions about parameter values that are not yet established. The objective of this study was to determine the most important cellular characteristics of a mechanoregulatory model describing both cell phenotype-specific and mechanobiological processes that are active during bone healing using a statistical approach. The computational model included an adaptive two-dimensional finite element model of a fractured long bone. Three different outcome criteria were quantified: (1) ability to predict sequential healing events, (2) amount of bone formation at early, mid and late stages of healing and (3) the total time until complete healing. For the statistical analysis, first a resolution IV fractional factorial design (L₆₄) was used to identify the most significant factors. Thereafter, a three-level Taguchi orthogonal array (L₂₇) was employed to study the curvature (non-linearity) of the 10 identified most important parameters. The results show that the ability of the model to predict the sequences of normal fracture healing was predominantly influenced by the rate of matrix production of bone, followed by cartilage degradation (replacement). The amount of bone formation at early stages was solely dependent on matrix production of bone and the proliferation rate of osteoblasts. However, the amount of bone formation at mid and late phases had the rate of matrix production of cartilage as the most influential parameter. The time to complete healing was primarily dependent on the rate of cartilage degradation during endochondral ossification, followed by the rate of cartilage formation. The analyses of the curvature revealed a linear response for parameters related to bone, where higher rates of formation were more beneficial to healing. In contrast, parameters related to fibrous tissue and cartilage showed optimum levels. Some fibrous connective tissue- and cartilage formation was beneficial to bone healing, but too much of either tissue delayed bone formation. The identified significant parameters and processes are further confirmed by in vivo animal experiments in the literature. This study illustrates the potential of design of experiments methods for evaluating computational mechanobiological model parameters and suggests that further experiments should preferably focus at establishing values of parameters related to cartilage formation and degradation.     

How plate positioning impacts the biomechanics of the open wedge tibial osteotomy; A finite element analysis

A numerical model of the medial open wedge tibial osteotomy based on the finite element method was developed. Two plate positions were tested numerically. In a configuration, (a), the plate was fixed in a medial position and (b) in an anteromedial position. The simulation took into account soft tissues preload, muscular tonus and maximal gait load.

The maximal stresses observed in the four structural elements (bone, plate, wedge, screws) of an osteotomy with plate in medial position were substantially higher (1.13–2.8 times more) than those observed in osteotomy with an anteromedial plate configuration. An important increase (1.71 times more) of the relative micromotions between the wedge and the bone was also observed. In order to avoid formation of fibrous tissue at the bone wedge interface, the osteotomy should be loaded under 18.8% (～50 kg) of the normal gait load until the osteotomy interfaces union is achieved.     

Reliable simulations of the human proximal femur by high-order finite element analysis validated by experimental observations

Background: The mechanical response of patient-specific bone to various load conditions is of major clinical importance in orthopedics. Herein we enhance the methods presented in Yosibash et al. [2007. A CT-based high-order finite element analysis of the human proximal femur compared to in-vitro experiments. ASME Journal of Biomechanical Engineering 129(3), 297–309.] for the reliable simulations of the human proximal femur by high-order finite elements (FEs) and validate the simulations by experimental observations.

Method of approach: A fresh-frozen human femur was scanned by quantitative computed tomography (QCT) and thereafter loaded (in vitro experiments) by a quasi-static force of up to 1250 N. QCT scans were manipulated to generate a high-order FE bone model with distinct cortical and trabecular regions having inhomogeneous isotropic elastic properties with Young's modulus represented by continuous spatial functions. Sensitivity analyses were performed to quantify parameters that mostly influence the mechanical response. FE results were compared to displacements and strains measured in the experiments.

Results: Young moduli correlated to QCT Hounsfield Units by relations in Keyak and Falkinstein [2003. Comparison of in situ and in vitro CT scan-based finite element model predictions of proximal femoral fracture load. Medical Engineering and Physics 25, 781–787.] were found to provide predictions that match the experimental results closely. Excellent agreement was found for both the displacements and strains. The presented study demonstrates that reliable and validated high-order patient-specific FE simulations of human femurs based on QCT data are achievable for clinical computer-aided decision making.     

股骨骨折术后1 年随访骨愈合模型的有限元分析

目的：对股骨骨折髓内钉术后1 年骨愈合模型快速建模,通过有限元分析研究对比术前术后模型,通过术前判定内固定取出后骨折断端是否断裂。方法：运用Mimics、Geomagic Studio、Abaqus等软件采用快速个体化建模方法对股骨骨折髓内钉术后1年内固定取出术前后的多层螺旋CT 数据进行快速建立模型,术前模型模拟剥除钢板后进行有限元分析,施加人体单腿站立时的静力载荷和约束,并将分析结果与术后模型进行对比,观察米塞斯应力分布情况、最大值及其所处部位。结果：按照材料属性进行区别显示米赛斯应力的最大值及最小值,在不同应力载荷下,手术前后各类型材料的米赛斯应力最大值及最小值部位相同,各类型材料中,最大值均没有位于骨折断端,不同方法的最大应力值部位相近,均在股骨中远端1/4 交界处,手术前后应力分布基本相同。结论：采用个体化建模方法可以对骨折内固定取出前的骨愈合模型进行运算分析,快速预判术后是否导致骨折断端断裂。     

股骨骨折术后1年随访骨愈合模型的有限元分析

目的：对股骨骨折髓内钉术后1年骨愈合模型快速建模,通过有限元分析研究对比术前术后模型,通过术前判定内固定取出后骨折断端是否断裂。方法：运用Mimics、Geomagic Studio、Abaqus等软件采用快速个体化建模方法对股骨骨折髓内钉术后1年内固定取出术前后的多层螺旋CT数据进行快速建立模型,术前模型模拟剥除钢板后进行有限元分析,施加人体单腿站立时的静力载荷和约束,并将分析结果与术后模型进行对比,观察米塞斯应力分布情况、最大值及其所处部位。结果：按照材料属性进行区别显示米赛斯应力的最大值及最小值,在不同应力载荷下,手术前后各类型材料的米赛斯应力最大值及最小值部位相同,各类型材料中,最大值均没有位于骨折断端,不同方法的最大应力值部位相近,均在股骨中远端1/4交界处,手术前后应力分布基本相同。结论：采用个体化建模方法可以对骨折内固定取出前的骨愈合模型进行运算分析,快速预判术后是否导致骨折断端断裂。     

Continuum remodeling revisited

Recent research effort in bone remodeling has been directed toward describing interstitial fluid flow in the lacuno-canalicular system and its potential as a cellular stimulus. Regardless of the precise contents of the mechanotransduction “black box”, it seems clear that the fluid flow on which the remodeling is predicated cannot occur under static loading conditions. In an attempt to help continuum remodeling simulations catch up with cellular and subcellular research, this paper presents a simple, strain rate driven remodeling algorithm for density allocation and principal material direction rotations. An explicit finite element code was written and deployed on a supercomputer which discretizes the remodeling process and uses an objective hypoelastic constitutive law to simulate trabecular realignment. Results indicate that a target strain rate for this dynamic approach is |D _I ^*| = 1.7% per second which seems reasonable when compared to observed strain rates. Simulations indicate that a morpho-mechanically realistic three-dimensional bone can be synthesized by applying a few dynamic loads at the envelope of common daily physiological rates, even with no static loading component.     

A novel approach to estimate trabecular bone anisotropy using a database approach

Continuum finite element (FE) models of bones have become a standard pre-clinical tool to estimate bone strength. These models are usually based on clinical CT scans and material properties assigned are chosen as isotropic based only on the density distribution. It has been shown, however, that trabecular bone elastic behavior is best described as orthotropic. Unfortunately, the use of orthotropic models in FE analysis derived from CT scans is hampered by the fact that the measurement of a trabecular orientation (fabric) is not possible from clinical CT images due to the low resolution of such images. In this study, we explore the concept of using a database (DB) of high-resolution bone models to derive the fabric information that is missing in clinical images. The goal of this study was to investigate if models with fabric derived from a relatively small database can already produce more accurate results than isotropic models.