Finite Element Assessment of Bone Fragility from Clinical Images.

PURPOSE OF REVIEW: We re-evaluated clinical applications of image-to-FE models to understand if clinical advantages are already evident, which proposals are promising, and which questions are still open. RECENT FINDINGS: CT-to-FE is useful in longitudinal treatment evaluation and groups discrimination. In metastatic lesions, CT-to-FE strength alone accurately predicts impending femoral fractures. In osteoporosis, strength from CT-to-FE or DXA-to-FE predicts incident fractures similarly to DXA-aBMD. Coupling loads and strength (possibly in dynamic models) may improve prediction. One promising MRI-to-FE workflow may now be tested on clinical data. Evidence of artificial intelligence usefulness is appearing. CT-to-FE is already clinical in opportunistic CT screening for osteoporosis, and risk of metastasis-related impending fractures. Short-term keys to improve image-to-FE in osteoporosis may be coupling FE with fall risk estimates, pool FE results with other parameters through robust artificial intelligence approaches, and increase reproducibility and cross-validation of models. Modeling bone modifications over time and bone fracture mechanics are still open issues.

Bone adaptation of a biologically reconstructed femur after Ewing sarcoma: Long-term morphological and densitometric evolution.

Combining bone allografts and vascularized fibular autografts in intercalary reconstructions after resection of bone sarcomas is of particular interest in young patients as it facilitates bone healing and union and helps reduce fractures. However, adverse events related to bone adaptation still occur. Bone adaptation is driven by mechanical loading, but no quantitative biomechanical studies exist that would help surgical planning and rehabilitation. We analyzed the bone adaptation of a successful femoral reconstruction after Ewing sarcoma during 76-month follow-up using a novel methodology that allows CT-based quantification of morphology and density. The results indicated that the vital allograft promoted bone adaptation in the reconstruction. However, an overall negative balance of bone remodeling and a progressive mineral density decrease in the femoral neck might threaten long-term bone safety. These concerns seem related to both surgical technique and mechanical stimuli, where a stiff metal implant may determine load sharing, which negatively affects bone remodeling.

Navigation in Orthognathic Surgery: 3D Accuracy.

This article aims to determine the absolute accuracy of maxillary repositioning during orthognathic surgery according to simulation-guided navigation, that is, the combination of navigation and three-dimensional (3D) virtual surgery. We retrospectively studied 15 patients treated for asymmetric dentofacial deformities at the Oral and Maxillofacial Surgery Unit of the S.Orsola-Malpighi University Hospital in Bologna, Italy, from January 2010 to January 2012. Patients were scanned with a cone-beam computed tomography before and after surgery. The virtual surgical simulation was realized with a dedicated software and loaded on a navigation system to improve intraoperative reproducibility of the preoperative planning. We analyzed the outcome following two protocols: (1) planning versus postoperative 3D surface analysis; (2) planning versus postoperative point-based analysis. For 3D surface comparison, the mean Hausdorff distance was measured, and median among cases was 0.99 mm. Median reproducibility < 1 mm was 61.88% and median reproducibility < 2 mm was 85.46%. For the point-based analysis, with sign, the median distance was 0.75 mm in the frontal axis, -0.05 mm in the caudal-cranial axis, -0.35 mm in the lateral axis. In absolute value, the median distance was 1.19 mm in the frontal axis, 0.59 mm in the caudal-cranial axis, and 1.02 mm in the lateral axis. We suggest that simulation-guided navigation makes accurate postoperative outcomes possible for maxillary repositioning in orthognathic surgery, if compared with the surgical computer-designed project realized with a dedicated software, particularly for the vertical dimension, which is the most challenging to manage.

Tensile Yield Strain of Human Cortical Bone from the Femoral Diaphysis Is Constant among Healthy Adults and across the Anatomical Quadrants.

The literature suggests that the yield strain of cortical bone is invariant to its stiffness (elastic modulus) and strength (yield stress). However, data about intra-individual variations, e.g., the influence of different collagen/mineral organisations observed in bone aspects withstanding different habitual loads, are lacking. The hypothesis that the yield strain of human cortical bone tissue, retrieved from femoral diaphyseal quadrants subjected to different habitual loads, is invariant was tested. Four flat dumbbell-shaped specimens were machined from each quadrant of the proximal femoral diaphysis of five adult donors for a total of 80 specimens. Two extensometers attached to the narrow specimen region were used to measure deformation during monotonic tensile testing. The elastic modulus (linear part of the stress-strain curve) and yield strain/stress at a 0.2% offset were obtained. Elastic modulus and yield stress values were, respectively, in the range of 12.2-20.5 GPa and 75.9-136.6 MPa and exhibited a positive linear correlation. All yield strain values were in the narrow range of 0.77-0.87%, regardless of the stiffness and strength of the tissue and the anatomical quadrant. In summary, the results corroborate the hypothesis that tensile yield strain in cortical bone is invariant, irrespective also of the anatomical quadrant. The mean yield strain value found in this study is similar to what was reported by inter-species and evolution studies but slightly higher than previous reports in humans, possibly because of the younger age of our subjects. Further investigations are needed to elucidate a possible dependence of yield strain on age.

Does a novel bridging collar in endoprosthetic replacement optimise the mechanical environment for osseointegration? A finite element study.

Introduction: Limb-salvage surgery using endoprosthetic replacements (EPRs) is frequently used to reconstruct segmental bone defects, but the reconstruction longevity is still a major concern. In EPRs, the stem-collar junction is the most critical region for bone resorption. We hypothesised that an in-lay collar would be more likely to promote bone ongrowth in Proximal Femur Reconstruction (PFR), and we tested this hypothesis through validated Finite Element (FE) analyses simulating the maximum load during walking. Methods: We simulated three different femur reconstruction lengths (proximal, mid-diaphyseal, and distal). For each reconstruction length one in-lay and one traditional on-lay collar model was built and compared. All reconstructions were virtually implanted in a population-average femur. Personalised Finite Element models were built from Computed Tomography for the intact case and for all reconstruction cases, including contact interfaces where appropriate. We compared the mechanical environment in the in-lay and on-lay collar configurations, through metrics of reconstruction safety, osseointegration potential, and risk of long-term bone resorption due to stress-shielding. Results: In all models, differences with respect to intact conditions were localized at the inner bone-implant interface, being more marked in the collar-bone interface. In proximal and mid-diaphyseal reconstructions, the in-lay configuration doubled the area in contact at the bone-collar interface with respect to the on-lay configuration, showed less critical values and trends of contact micromotions, and consistently showed higher (roughly double) volume percentages of predicted bone apposition and reduced (up to one-third) percentages of predicted bone resorption. In the most distal reconstruction, results for the in-lay and on-lay configurations were generally similar and showed overall less favourable maps of the bone remodelling tendency. Discussion: In summary, the models corroborate the hypothesis that an in-lay collar, by realising a more uniform load transfer into the bone with a more physiological pattern, creates an advantageous mechanical environment at the bone-collar interface, compared to an on-lay design. Therefore, it could significantly increase the survivorship of endo-prosthetic replacements.

Applying a homogeneous pressure distribution to the upper vertebral endplate: Validation of a new loading system, pilot application to human vertebral bodies, and finite element predictions of DIC measured displacements and strains.

Image-based personalized Finite Element Models (pFEM) could detect alterations in physiological deformation of human vertebral bodies, but their accuracy has been seldom reported. Meaningful validation experiments should allow vertebral endplate deformability and ensure well-controlled boundary conditions. This study aimed to (i) validate a new loading system to apply a homogeneous pressure on the vertebral endplate during vertebral body compression regardless of endplate deformation; (ii) perform a pilot study on human vertebral bodies measuring surface displacements and strains with Digital Image Correlation (DIC); (iii) determine the accuracy of pFEM of the vertebral bodies. Homogeneous pressure application was achieved by pressurizing a fluid silicone encased in a rubber silicone film acting on the cranial endplate. The loading system was validated by comparing DIC-measured longitudinal strains and lower-end contact pressures, measured on three homogeneous pseudovertebrae of constant transversal section at 2.0 kN, against theoretically calculated values. Longitudinal strains and contact pressures were rather homogeneous, and their mean values close to theoretical calculations (5% underestimation). DIC measurements of surface longitudinal and circumferential displacements and strains were obtained on three human vertebral bodies at 2.0 kN. Complete displacement and strain maps were achieved for anterolateral aspects with random errors ≤0.2 µm and ≤30 µstrain, respectively. Venous plexus and double curvatures limited the completeness and accuracy of DIC data in posterior aspects. pFEM of vertebral bodies, including cortical bone mapping, were built from computed tomography images. In anterolateral aspects, pFEM accuracy of the three vertebrae was: (i) comparable to literature in terms of longitudinal displacements (R2>0.8); (ii) extended to circumferential displacements (pooled data: R2>0.9) and longitudinal strains (zero median error, 95% error: <27%). Circumferential strains were overestimated (median error: 39%). The new methods presented may permit to study how physiological and pathologic conditions influence the ability of vertebral endplates/bodies to sustain loads.

High-resolution peripheral quantitative computed tomography: research or clinical practice?

High-resolution peripheral quantitative CT (HR-pQCT) is a low-dose three-dimensional imaging technique, originally developed for in vivo assessment of bone microarchitecture at the distal radius and tibia in osteoporosis. HR-pQCT has the ability to discriminate trabecular and cortical bone compartments, providing densitometric and structural parameters. At present, HR-pQCT is mostly used in research settings, despite evidence showing that it may be a valuable tool in osteoporosis and other diseases. This review summarizes the main applications of HR-pQCT and addresses the limitations that currently prevent its integration into routine clinical practice. In particular, the focus is on the use of HR-pQCT in primary and secondary osteoporosis, chronic kidney disease (CKD), endocrine disorders affecting bone, and rare diseases. A section on novel potential applications of HR-pQCT is also present, including assessment of rheumatic diseases, knee osteoarthritis, distal radius/scaphoid fractures, vascular calcifications, effect of medications, and skeletal muscle. The reviewed literature seems to suggest that a more widespread implementation of HR-pQCT in clinical practice would offer notable opportunities. For instance, HR-pQCT can improve the prediction of incident fractures beyond areal bone mineral density provided by dual-energy X-ray absorptiometry. In addition, HR-pQCT may be used for the monitoring of anti-osteoporotic therapy or for the assessment of mineral and bone disorder associated with CKD. Nevertheless, several obstacles currently prevent a broader use of HR-pQCT and would need to be targeted, such as the small number of installed machines worldwide, the uncertain cost-effectiveness, the need for improved reproducibility, and the limited availability of reference normative data sets.

Individual Trajectories of Bone Mineral Density Reveal Persistent Bone Loss in Bone Sarcoma Patients: A Retrospective Study.

Multiagent chemotherapy offers an undoubted therapeutic benefit to cancer patients, but is also associated with chronic complications in survivors. Osteoporosis affects the quality of life of oncologic patients, especially at the paediatric age. However, very few studies have described the extent of loss of bone mineral density (BMD) in bone sarcoma patients. We analysed a retrospective series of children and adolescents with primary malignant bone tumours (52 osteosarcoma and 31 Ewing sarcoma) and retrieved their BMD at diagnosis and follow-up as Hounsfield units (HU). We studied their individual BMD trajectories before and after chemotherapy up to 5 years, using routine chest CT scan and attenuation thresholds on T12 vertebrae ROI. At one year, bone sarcoma patients showed significant bone loss compared to diagnosis: 17.6% and 17.1% less for OS and EW, respectively. Furthermore, a bone loss of more than 49.2 HU at one-year follow-up was predictive of the persistence of a reduced bone mass over the following 4 years, especially in patients with EW. At 4 years, only 26% and 12.5% of OS and EW, respectively, had recovered or improved their BMD with respect to the onset, suggesting a risk of developing morbidities related to a low BMD in those subjects.

A Novel Quantitative Computer-Assisted Score Can Improve Repeatability in the Estimate of Vascular Calcifications at the Abdominal Aorta.

In CKD and in the elderly, Vascular Calcifications (VC) are associated to cardiovascular events and bone fractures. VC scores at the abdominal aorta (AA) from lateral spine radiographs are widely applied (the 0-24 semiquantitative discrete visual score (SV) being the most used). We hypothesised that a novel continuum score based on quantitative computer-assisted tracking of calcifications (QC score) can improve the precision of the SV score. This study tested the repeatability and reproducibility of QC score and SV score. In forty-four patients with VC from an earlier study, five experts from four specialties evaluated the data twice using a dedicated software. Test-retest was performed on eight subjects. QC results were reported in a 0-24 scale to readily compare with SV. The QC score showed higher intra-operator repeatability: the 95% CI of Bland-Altman differences was almost halved in QC; intra-operator R2 improved from 0.67 for SV to 0.79 for QC. Inter-observer repeatability was higher for QC score in the first (Intraclass Correlation Coefficient 0.78 vs. 0.64), but not in the second evaluation (0.84 vs. 0.82), indicating a possible heavier learning artefact for SV. The Minimum Detectable Difference (MDD) was smaller for QC (2.98 vs. 4 for SV, in the 0-24 range). Both scores were insensitive to test-retest procedure. Notably, QC and SV scores were discordant: SV showed generally higher values, and an increasing trend of differences with VC severity. In summary, the new QC score improved the precision of lateral spine radiograph scores in estimating VC. We reported for the first time an estimate of MDD in VC assessment that was 25% lower for the new QC score with respect to the usual SV score. An ongoing study will determine whether this lower MDD may reduce follow-up times to check for VC progression.

A taper-fit junction to improve long bone reconstruction: A parametric In Silico model.

PURPOSE: Critical size long bone defects represent a clinical challenge in orthopaedic surgery. Various grafting techniques have been developed through the years, but they all present several downsides. A key requirement of all grafting techniques is the achievement of a continuous interface between host bone and graft to enhance both biological processes and mechanical stability. This study used a parametric in silico model to quantify the biomechanical effect of the inaccuracies inherent to current osteotomy techniques, and to test a new concept of accurate taper-fit junction that may improve the biomechanical parameters of the reconstruction under load. METHODS: A population-based in-silico 3D model of the reconstruction of a long bone defect was built to represent a defect of the femoral mid-diaphysis. To fix the reconstruction a titanium plate was placed on the lateral aspect of the reconstruction. The model was modified to (i) quantify the biomechanical consequences of actual inaccuracies in the realization of a flat host-graft interface, (ii) compare the contact behaviour and bone strains among different taper angles of the new design and the current host-graft flat interface, (iii) evaluate the robustness of the taper-fit design to inter-subject variability in bone geometry and defect length. RESULTS: The influence of 2° single-plane misalignments of the host-graft interface is highly dependent on the misalignment orientation with respect to the metal plate. For some misalignment orientations, tangential micromotions of contact interfaces exceeded alert thresholds. When the angle of the taper-fit host-graft junction is changed from 10° to 30° and the results obtained are compared with the planar case, the overall stiffness is almost preserved, the bone strains are almost unchanged with safety factors higher than five, and full contact closure around the host-graft junction is achieved at 20°. Similarly, contact pressures decrease almost linearly with a 20% decrease at 30°. The host-graft micro motions are almost unchanged in both value and distribution up to 20° and never exceed the warning threshold of 50 µm. CONCLUSIONS: The present in silico study developed quantitative biomechanical evidence that an osteotomy performed with attention to the perpendicularity of the cut planes is needed to reduce the risk of mismatch and possible complications of long bone reconstructions, and that a new concept of a taper-fit junction may improve the biomechanical environment of the interface between the graft and the host bone. The optimal taper-fit configuration is suggested to be around a 20° taper angle. These results will serve as an input to conduct exvivo experiments to further corroborate the proposed taper-fit junction concept and to refine its surgical implementation.

The Vessels-Bone Axis: Iliac Artery Calcifications, Vertebral Fractures and Vitamin K from VIKI Study.

Vascular calcification and fragility fractures are associated with high morbidity and mortality, especially in end-stage renal disease. We evaluated the relationship of iliac arteries calcifications (IACs) and abdominal aortic calcifications (AACs) with the risk for vertebral fractures (VFs) in hemodialysis patients. The VIKI study was a multicenter cross-sectional study involving 387 hemodialysis patients. The biochemical data included bone health markers, such as vitamin K levels, vitamin K-dependent proteins, vitamin 25(OH)D, alkaline phosphatase, parathormone, calcium, and phosphate. VF, IACs and AACs was determined through standardized spine radiograms. VF was defined as >20% reduction of vertebral body height, and VC were quantified by measuring the length of calcium deposits along the arteries. The prevalence of IACs and AACs were 56.1% and 80.6%, respectively. After adjusting for confounding variables, the presence of IACs was associated with 73% higher odds of VF (p = 0.028), whereas we found no association (p = 0.294) for AACs. IACs were associated with VF irrespective of calcification severity. Patients with IACs had lower levels of vitamin K2 and menaquinone 7 (0.99 vs. 1.15 ng/mL; p = 0.003), and this deficiency became greater with adjustment for triglycerides (0.57 vs. 0.87 ng/mL; p < 0.001). IACs, regardless of their extent, are a clinically relevant risk factor for VFs. The association is enhanced by adjusting for vitamin K, a main player in bone and vascular health. To our knowledge these results are the first in the literature. Prospective studies are needed to confirm these findings both in chronic kidney disease and in the general population.

Cortical bone mapping improves finite element strain prediction accuracy at the proximal femur.

Despite evidence of the biomechanical role of cortical bone, current state of the art finite element models of the proximal femur built from clinical CT data lack a subject-specific representation of the bone cortex. Our main research hypothesis is that the subject-specific modelling of cortical bone layer from CT images, through a deconvolution procedure known as Cortical Bone Mapping (CBM, validated for cortical thickness and density estimates) can improve the accuracy of CT-based FE models of the proximal femur, currently limited by partial volume artefacts. Our secondary hypothesis is that a careful choice of cortical-specific density-elasticity relationship may improve model accuracy. We therefore: (i) implemented a procedure to include subject-specific CBM estimates of both cortical thickness and density in CT-based FE models. (ii) defined alternative models that included CBM estimates and featured a cortical-specific or an independently optimised density-elasticity relationship. (iii) tested our hypotheses in terms of elastic strain estimates and failure load and location prediction, by comparing with a published cohort of 14 femurs, where strain and strength in stance and fall loading configuration were experimentally measured, and estimated through reference FE models that did not explicitly model the cortical compartment. Our findings support the main hypothesis: an explicit modelling of the proximal femur cortical bone layer including CBM estimates of cortical bone thickness and density increased the FE strains prediction, mostly by reducing peak errors (average error reduced by 30%, maximum error and 95th percentile of error distribution halved) and especially when focusing on the femoral neck locations (all error metrics at least halved). We instead rejected the secondary hypothesis: changes in cortical density-elasticity relationship could not improve validation performances. From these improved baseline strain estimates, further work is needed to achieve accurate strength predictions, as models incorporating cortical thickness and density produced worse estimates of failure load and equivalent estimates of failure location when compared to reference models. In summary, we recommend including local estimates of cortical thickness and density in FE models to estimate bone strains in physiological conditions, and especially when designing exercise studies to promote bone strength.

Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro.

No agreement on the choice of the failure criterion to adopt for the bone tissue can be found in the literature among the finite element studies aiming at predicting fracture risk of bones. The use of stress-based criteria seems to prevail on strain-based ones, while basic bone biomechanics suggest using strain parameters to describe failure. The aim of the present combined experimental-numerical study was to verify, using subject-specific finite element models able to accurately predict strains, if a strain-based failure criterion could identify the failure patterns of bones. Three cadaver femurs were CT-scanned and subsequently fractured in a clinically relevant single-stance loading scenario. Load-displacement curves and high-speed movies were acquired to define the failure load and the location of fracture onset, respectively. Subject-specific finite element models of the three femurs were built from CT data following a validated procedure. A maximum principal strain criterion was implemented in the finite element models, and two stress-based criteria selected for comparison. The failure loads measured were applied to the models, and the computed risks of fracture were compared to the results of the experimental tests. The proposed principal strain criterion managed to correctly identify the level of failure risk and the location of fracture onset in all the modelled specimens, while Von Mises or maximum principal stress criteria did not give significant information. A maximum principal strain criterion can thus be defined a suitable candidate for the in vivo risk factor assessment on long bones.

An accurate estimation of bone density improves the accuracy of subject-specific finite element models.

An experimental-numerical study was performed to investigate the relationships between computed tomography (CT)-density and ash density, and between ash density and apparent density for bone tissue, to evaluate their influence on the accuracy of subject-specific FE models of human bones. Sixty cylindrical bone specimens were examined. CT-densities were computed from CT images while apparent and ash densities were measured experimentally. The CT/ash-density and ash/apparent-density relationships were calculated. Finite element models of eight human femurs were generated considering these relationships to assess their effect on strain prediction accuracy. CT and ash density were linearly correlated (R(2)=0.997) over the whole density range but not equivalent (intercep t <0, slope >1). A constant ash/apparent-density ratio (0.598+/-0.004) was found for cortical bone. A lower ratio, with a larger dispersion, was found for trabecular bone (0.459+/-0.100), but it became less dispersed, and equal to that of cortical tissue, when testing smaller trabecular specimens (0.598+/-0.036). This suggests that an experimental error occurred in apparent-density measurements for large trabecular specimens and a constant ratio can be assumed valid for the whole density range. Introducing the obtained relationships in the FE modelling procedure improved strain prediction accuracy (R(2)=0.95, RMSE=7%). The results suggest that: (i) a correction of the densitometric calibration should be used when evaluating bone ash-density from clinical CT scans, to avoid ash-density underestimation and overestimation for low- and high-density bone tissue, respectively; (ii) the ash/apparent-density ratio can be assumed constant in human femurs and (iii) the correction improves significantly the model accuracy and should be considered in subject-specific bone modelling.

A modified method for assigning material properties to FE models of bones.

The aim of the present study is to compare the results from subject-specific finite element analysis (FEA) of a human femur to experimental measurements, using two different methods for assigning material properties to the FE models. A modified material mapping strategy allowing for spatial variation of material properties within the elements and Young's modulus surface corrections is presented and compared to a more conventional strategy, whereby constant material properties are assigned to each element. The accuracy of the superficial stress-strain predictions was evaluated against experimental results from 13 strain gauges and five different load cases. Both methods predicted stresses with acceptable accuracy (R(2) = 0.92, root mean square error, RMSE < 10%), with the conventional method performing slightly better. The modified method performed better in strain prediction (R(2) = 0.85, RMSE = 23% versus R(2) = 0.79, RMSE = 31%).

Mathematical relationships between bone density and mechanical properties: a literature review.

BACKGROUND: In many published studies, elastic properties of bone are correlated to the bone density, in order to derive an empirical elasticity-density relationship. The most common use of these relationships is the prediction of the bone local properties from medical imaging data in subject-specific numerical simulation studies. The proposed relationships are substantially different one from the other. It is unclear whether such differences in elasticity-density relationships can be entirely explained in terms of methodological discrepancies among studies. METHODS: All relevant literature was reviewed. Only elasticity-density relationships derived from similarly controlled experiments were included and properly normalized. The resulting relationships were grouped according to the most important methodological differences: type of end support during testing, specimen geometry, and anatomical sampling location. FINDINGS: Even after normalization with respect to strain rate and densitometric measurement unit, substantial inter-study differences do exist, and they can only be partially explained by the methodological differences between studies. INTERPRETATION: Some recommendations are made for the application of elasticity-density relationships to subject-specific finite element studies. The importance of defining a standardized mechanical testing methodology for bone specimens is stressed, and some guidelines that emerged from the literature are proposed. To identify density-elasticity relationships suitable for use in subject-specific FE studies, the development of a benchmark study is also proposed, where the elasticity-density relationship is taken as the variable under study, and a numerical model of known numerical accuracy predicts experimental strain measurements.

Subject-specific finite element models can accurately predict strain levels in long bones.

The prediction of the stress-state and fracture risk induced in bones by various loading conditions in individual patients using subject-specific finite element models still represents a challenge in orthopaedic biomechanics. The accuracy of the strain predictions reported in the literature is variable and generally not satisfactory. The aim of the present study was to evaluate if a proper choice of the density-elasticity relationship can lead to accurate strain predictions in the frame of an automatic subject-specific model generation strategy. To this aim, a combined numerical-experimental study was performed comparing finite element predicted strains with strain-gauges measurements obtained on eight cadaver proximal femurs, each instrumented with 15 rosettes mostly concentrated in the bone metaphyses, tested non-destructively in vitro under six different loading scenarios. Three different density-elasticity power relationships were selected from the literature and implemented in the finite element models derived from computed tomography data. The results of the present study confirm the great influence of the density-elasticity relationship used on the accuracy of numerical predictions. One of the tested constitutive laws provided a very good agreement (R(2)=0.91, RMSE lower than 10% of the maximum measured value) between numerical calculations and experimental measurements. The presented results show, in addition, that the adoption of a single density-elasticity relationship over the whole bone density range is adequate to obtain an accuracy that is already suitable for many applications.

The material mapping strategy influences the accuracy of CT-based finite element models of bones: an evaluation against experimental measurements.

Aim of the present study was to evaluate the influence on the global model's accuracy of the strategy adopted to define the average element Young's modulus in subject-specific finite element models of bones from computed tomography data. The classic strategy of calculating the Young's modulus from an average element density and the one that averages the Young's moduli directly derived from each voxel Hounsfield Unit were considered. These strategies were applied to the finite element model of a real human femur. The accuracy of the superficial stress and strain predictions was evaluated against experimentally measured values in 13 strain-gauge locations for five different loading conditions. The results obtained for the two material distributions were statistically different. Both models predicted very accurately the superficial stresses, with regression coefficients higher than 0.9 and slopes not significantly different from unity. The second strategy definitely improved the strains prediction accuracy: the regression coefficient raised from 0.69 to 0.79; the average and peak errors decreased from 45.1% to 31.3% and from 228% to 134% of the maximum measured strain, respectively. The stress fields predicted inside the bone were also significantly different. A new software implementing the second strategy was made available in the public domain.

nmsBuilder: Freeware to create subject-specific musculoskeletal models for OpenSim.

BACKGROUND AND OBJECTIVE: Musculoskeletal modeling and simulations of movement have been increasingly used in orthopedic and neurological scenarios, with increased attention to subject-specific applications. In general, musculoskeletal modeling applications have been facilitated by the development of dedicated software tools; however, subject-specific studies have been limited also by time-consuming modeling workflows and high skilled expertise required. In addition, no reference tools exist to standardize the process of musculoskeletal model creation and make it more efficient. Here we present a freely available software application, nmsBuilder 2.0, to create musculoskeletal models in the file format of OpenSim, a widely-used open-source platform for musculoskeletal modeling and simulation. nmsBuilder 2.0 is the result of a major refactoring of a previous implementation that moved a first step toward an efficient workflow for subject-specific model creation. METHODS: nmsBuilder includes a graphical user interface that provides access to all functionalities, based on a framework for computer-aided medicine written in C++. The operations implemented can be used in a workflow to create OpenSim musculoskeletal models from 3D surfaces. A first step includes data processing to create supporting objects necessary to create models, e.g. surfaces, anatomical landmarks, reference systems; and a second step includes the creation of OpenSim objects, e.g. bodies, joints, muscles, and the corresponding model. RESULTS: We present a case study using nmsBuilder 2.0: the creation of an MRI-based musculoskeletal model of the lower limb. The model included four rigid bodies, five degrees of freedom and 43 musculotendon actuators, and was created from 3D surfaces of the segmented images of a healthy subject through the modeling workflow implemented in the software application. CONCLUSIONS: We have presented nmsBuilder 2.0 for the creation of musculoskeletal OpenSim models from image-based data, and made it freely available via nmsbuilder.org. This application provides an efficient workflow for model creation and helps standardize the process. We hope this would help promote personalized applications in musculoskeletal biomechanics, including larger sample size studies, and might also represent a basis for future developments for specific applications.

Relationship between bone adaptation and in-vivo mechanical stimulus in biological reconstructions after bone tumor: A biomechanical modeling analysis.

BACKGROUND: Biomechanical interpretations of bone adaptation in biological reconstructions following bone tumors would be crucial for orthopedic oncologists, particularly if based on quantitative observations. This would help plan for surgical treatments, rehabilitative programs and communication with the patients. We aimed to analyze the biomechanical adaptation of a femoral reconstruction after Ewing sarcoma according to an increasingly-used surgical technique, and to relate in-progress bone resorption to the mechanical stimulus induced by different motor activities. METHODS: We created a multiscale musculoskeletal and finite element model from CT scans and motion analysis data at a 76-month follow-up of a patient, to analyze muscle and joint loads, and to compare the mechanical competence of the reconstructed bone with the contralateral limb, in the current real condition and in a possible revision surgery that removed proximal screws. FINDINGS: Our results showed strategies of muscle coordination that led to differences in joint loads between limbs more marked in more demanding motor activities, and generally larger in the contralateral limb. The operated femur presented a markedly low ratio of physiological strain due to load-sharing with the metal implant, particularly in the lateral aspect. The possible revision surgery would help restore a physiological strain configuration, while the safety of the reconstruction would not be threatened. INTERPRETATION: We suggest that bone resorption is related to load-sharing and to the internal forces exerted during movement, and the mechanical stimulus should be improved by adopting modifications in the surgical treatment and by promoting physical therapy aimed at specific muscle strengthening.