Simple finite element models for use in the design of therapeutic footwear

Therapeutic footwear is frequently prescribed in cases of rheumatoid arthritis and diabetes to relieve or redistribute high plantar pressures in the region of the metatarsal heads. Few guidelines exist as to how these interventions should be designed and what effect such interventions actually have on the plantar pressure distribution. Finite element analysis has the potential to assist in the design process by refining a given intervention or identifying an optimal intervention without having to actually build and test each condition. However, complete and detailed foot models based on medical image segmentation have proven time consuming to build and computationally expensive to solve, hindering their utility in practice. Therefore, the goal of the current work was to determine if a simplified patient-specific model could be used to assist in the design of foot orthoses to reduce the plantar pressure in the metatarsal head region. The approach is illustrated by a case study of a diabetic patient experiencing high pressures and pain over the fifth metatarsal head. The simple foot model was initially calibrated by adjusting the individual loads on the metatarsals to approximate measured peak plantar pressure distributions in the barefoot condition to within 3%. This loading was used in various shod conditions to identify an effective orthosis. Model results for metatarsal pads were considerably higher than measured values but predictions for uniform surfaces were generally within 16% of measured values. The approach enabled virtual prototyping of the orthoses, identifying the most favorable approach to redistribute the patient’s plantar pressures.     

Reduction of plantar heel pressures: Insole design using finite element analysis

Plantar heel pain is a common condition that is often exacerbated by the repetitive stresses of walking. Treatment usually includes an in-shoe intervention designed to reduce plantar pressure under the heel by using insoles and a variety of off-the-shelf products. The design process for these products is often intuitive in nature and does not always rely on scientifically derived guidelines. Finite element analysis provides an efficient computational framework to investigate the performance of a large number of designs for optimal plantar pressure reduction. In this study, we used two-dimensional plane strain finite element modeling to investigate 27 insole designs. Combinations of three insole conformity levels (flat, half conforming, full conforming), three insole thickness values (6.3, 9.5 and 12.7 mm) and three insole materials (Poron Cushioning, Microcel Puff Lite and Microcel Puff) were simulated during the early support phase of gait. Plantar pressures predicted by the model were validated by experimental trials conducted in the same subject whose heel was modeled by loading the bare foot on a rigid surface and on foam mats. Conformity of the insole was the most important design variable, whereas peak pressures were relatively insensitive to insole material selection. The model predicted a 24% relief in pressure compared to barefoot conditions when using flat insoles; the reduction increased up to 44% for full conforming insoles.     

Validation of plantar pressure simulations using finite and discrete element modelling in healthy and diabetic subjects

Plantar pressure simulation driven by integrated 3D motion capture data, using both a finite element and a discrete element model, is compared for ten healthy and ten diabetic neuropathic subjects. The simulated peak pressure deviated on average between 16.7 and 34.2% from the measured peak pressure. The error in the position of the peak pressure was on average smaller than 4.2 cm. No method was more accurate than the other although statistical differences were found between them. Both techniques are thus complementary and useful tools to better understand the alteration of diabetic foot biomechanics during gait.     

Comparison of hexahedral and tetrahedral elements in finite element analysis of the foot and footwear

Finite element analysis has been widely used in the field of foot and footwear biomechanics to determine plantar pressures as well as stresses and strains within soft tissue and footwear materials. When dealing with anatomical structures such as the foot, hexahedral mesh generation accounts for most of the model development time due to geometric complexities imposed by branching and embedded structures. Tetrahedral meshing, which can be more easily automated, has been the approach of choice to date in foot and footwear biomechanics. Here we use the nonlinear finite element program Abaqus (Simulia, Providence, RI) to examine the advantages and disadvantages of tetrahedral and hexahedral elements under compression and shear loading, material incompressibility, and frictional contact conditions, which are commonly seen in foot and footwear biomechanics. This study demonstrated that for a range of simulation conditions, hybrid hexahedral elements (Abaqus C3D8H) consistently performed well while hybrid linear tetrahedral elements (Abaqus C3D4H) performed poorly. On the other hand, enhanced quadratic tetrahedral elements with improved stress visualization (Abaqus C3D10I) performed as well as the hybrid hexahedral elements in terms of contact pressure and contact shear stress predictions. Although the enhanced quadratic tetrahedral element simulations were computationally expensive compared to hexahedral element simulations in both barefoot and footwear conditions, the enhanced quadratic tetrahedral element formulation seems to be very promising for foot and footwear applications as a result of decreased labor and expedited model development, all related to facilitated mesh generation.     

Cardiovascular stent design and vessel stresses: a finite element analysis

Intravascular stents of various designs are currently in use to restore patency in atherosclerotic coronary arteries and it has been found that different stents have different in-stent restenosis rates. It has been hypothesized that the level of vascular injury caused to a vessel by a stent determines the level of restenosis. Computational studies may be used to investigate the mechanical behaviour of stents and to determine the biomechanical interaction between the stent and the artery in a stenting procedure. In this paper, we test the hypothesis that two different stent designs will provoke different levels of stress within an atherosclerotic artery and hence cause different levels of vascular injury. The stents analysed using the finite-element method were the S7 (Medtronic AVE) and the NIR (Boston Scientific) stent designs. An analysis of the arterial wall stresses in the stented arteries indicates that the modular S7 stent design causes lower stress to an atherosclerotic vessel with a localized stenotic lesion compared to the slotted tube NIR design. These results correlate with observed clinical restenosis rates, which have found higher restenosis rates in the NIR compared with the S7 stent design. Therefore, the testing methodology outlined here is proposed as a pre-clinical testing tool, which could be used to compare and contrast existing stent designs and to develop novel stent designs.     

Density-based load estimation using two-dimensional finite element models: a parametric study

A parametric investigation was conducted to determine the effects on the load estimation method of varying: (1) the thickness of back-plates used in the two-dimensional finite element models of long bones, (2) the number of columns of nodes in the outer medial and lateral sections of the diaphysis to which the back-plate multipoint constraints are applied and (3) the region of bone used in the optimization procedure of the density-based load estimation technique. The study is performed using two-dimensional finite element models of the proximal femora of a chimpanzee, gorilla, lion and grizzly bear. It is shown that the density-based load estimation can be made more efficient and accurate by restricting the stimulus optimization region to the metaphysis/epiphysis. In addition, a simple method, based on the variation of diaphyseal cortical thickness, is developed for assigning the thickness to the back-plate. It is also shown that the number of columns of nodes used as multipoint constraints does not have a significant effect on the method.     

Density-based load estimation using two-dimensional finite element models: a parametric study

A parametric investigation was conducted to determine the effects on the load estimation method of varying: (1) the thickness of back-plates used in the two-dimensional finite element models of long bones, (2) the number of columns of nodes in the outer medial and lateral sections of the diaphysis to which the back-plate multipoint constraints are applied and (3) the region of bone used in the optimization procedure of the density-based load estimation technique. The study is performed using two-dimensional finite element models of the proximal femora of a chimpanzee, gorilla, lion and grizzly bear. It is shown that the density-based load estimation can be made more efficient and accurate by restricting the stimulus optimization region to the metaphysis/epiphysis. In addition, a simple method, based on the variation of diaphyseal cortical thickness, is developed for assigning the thickness to the back-plate. It is also shown that the number of columns of nodes used as multipoint constraints does not have a significant effect on the method.     

Elastic properties of cancellous bone derived from finite element models of parameterized microstructure cells

Evaluation of constitutive properties of cancellous bone and their relationships to microstructural parameters is a crucial issue in analysis of stresses and strains in bone tissues and simulation of their remodelling. Known limitations of experimental methods as well as of the micro-FE techniques make the analysis and homogenization of 'equivalent' trabecular microstructures an advantageous tool for this task. In this study, parameterized orthotropic constitutive models of cancellous bone are derived from finite element analysis of repeatable microstructure cells. Two cell types are analysed: cube- and prism-based. The models are fully three-dimensional, have realistic curvilinear shapes and are parameterized with three shape parameters. Variation of the parameters allows to imitate most of the typical microstructure patterns observed in real bones, along with variety of intermediate geometries. Finite element models of cells are generated by a special-purpose structured mesh generator for any arbitrary set of shape parameter values. Six static numerical tests are performed for an exhaustive number of parameter value sets (microstructure instances). Multi-point boundary conditions imposed on the models ensure mutual fitting of deformed neighbouring cells. Values of computed stresses allow to determine all coefficients of elastic orthotropic stiffness matrix. Results have a form of tabularized functions of elastic constants versus the shape parameters. Comparison of the results with micro-FE data obtained for a large set of cancellous bone specimens proves a good agreement, though evidently better in the case of the prism-based cell model.     

Non-invasive prediction of the mouse tibia mechanical properties from microCT images: comparison between different finite element models

New treatments for bone diseases require testing in animal models before clinical translation, and the mouse tibia is among the most common models. In vivo micro-Computed Tomography (microCT)-based micro-Finite Element (microFE) models can be used for predicting the bone strength non-invasively, after proper validation against experimental data. Different modelling techniques can be used to estimate the bone properties, and the accuracy associated with each is unclear. The aim of this study was to evaluate the ability of different microCT-based microFE models to predict the mechanical properties of the mouse tibia under compressive load. Twenty tibiae were microCT scanned at 10.4 µm voxel size and subsequently compressed at 0.03 mm/s until failure. Stiffness and failure load were measured from the load–displacement curves. Different microFE models were generated from each microCT image, with hexahedral or tetrahedral mesh, and homogeneous or heterogeneous material properties. Prediction accuracy was comparable among models. The best correlations between experimental and predicted mechanical properties, as well as lower errors, were obtained for hexahedral models with homogeneous material properties. Experimental stiffness and predicted stiffness were reasonably well correlated (R² = 0.53–0.65, average error of 13–17%). A lower correlation was found for failure load (R² = 0.21–0.48, average error of 9–15%). Experimental and predicted mechanical properties normalized by the total bone mass were strongly correlated (R² = 0.75–0.80 for stiffness, R² = 0.55–0.81 for failure load). In conclusion, hexahedral models with homogeneous material properties based on in vivo microCT images were shown to best predict the mechanical properties of the mouse tibia.

     

Prediction of femoral strength using 3D finite element models reconstructed from DXA images: validation against experiments

Computed tomography (CT)-based finite element (FE) models may improve the current osteoporosis diagnostics and prediction of fracture risk by providing an estimate for femoral strength. However, the need for a CT scan, as opposed to the conventional use of dual-energy X-ray absorptiometry (DXA) for osteoporosis diagnostics, is considered a major obstacle. The 3D shape and bone mineral density (BMD) distribution of a femur can be reconstructed using a statistical shape and appearance model (SSAM) and the DXA image of the femur. Then, the reconstructed shape and BMD could be used to build FE models to predict bone strength. Since high accuracy is needed in all steps of the analysis, this study aimed at evaluating the ability of a 3D FE model built from one 2D DXA image to predict the strains and fracture load of human femora. Three cadaver femora were retrieved, for which experimental measurements from ex vivo mechanical tests were available. FE models were built using the SSAM-based reconstructions: using only the SSAM-reconstructed shape, only the SSAM-reconstructed BMD distribution, and the full SSAM-based reconstruction (including both shape and BMD distribution). When compared with experimental data, the SSAM-based models predicted accurately principal strains (coefficient of determination >0.83, normalized root-mean-square error <16%) and femoral strength (standard error of the estimate 1215 N). These results were only slightly inferior to those obtained with CT-based FE models, but with the considerable advantage of the models being built from DXA images. In summary, the results support the feasibility of SSAM-based models as a practical tool to introduce FE-based bone strength estimation in the current fracture risk diagnostics.     

A new approach for assigning bone material properties from CT images into finite element models

Generation of subject-specific finite element (FE) models from computed tomography (CT) datasets is of significance for application of the FE analysis to bone structures. A great challenge that remains is the automatic assignment of bone material properties from CT Hounsfield Units into finite element models. This paper proposes a new assignment approach, in which material properties are directly assigned to each integration point. Instead of modifying the dataset of FE models, the proposed approach divides the assignment procedure into two steps: generating the data file of the image intensity of a bone in a MATLAB program and reading the file into ABAQUS via user subroutines. Its accuracy has been validated by assigning the density of a bone phantom into a FE model. The proposed approach has been applied to the FE model of a sheep tibia and its applicability tested on a variety of element types. The proposed assignment approach is simple and illustrative. It can be easily modified to fit users’ situations.     

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.     

Subject-specific finite element models of long bones: An in vitro evaluation of the overall accuracy

The determination of the mechanical stresses induced in human bones is of great importance in both research and clinical practice. Since the stresses in bones cannot be measured non-invasively in vivo, the only way to estimate them is through subject-specific finite element modelling. Several methods exist for the automatic generation of these models from CT data, but before bringing them in the clinical practice it is necessary to assess their accuracy in the predictions of the bone stresses. Particular attention should be paid to those regions, like the epiphyseal and metaphyseal parts of long bones, where the automatic methods are typically less accurate. Aim of the present study was to implement a general procedure to automatically generate subject-specific finite element models of bones from CT data and estimate the accuracy of this general procedure by applying it to one real femur. This femur was tested in vitro under five different loading scenarios and the results of these tests were used to verify how the adoption of a simplified two-material homogeneous model would change the accuracy with respect to the density-based inhomogeneous one, with special attention paid to the epiphyseal and metaphyseal proximal regions of the bone. The results showed that the density-based inhomogeneous model predicts with a very good accuracy the measured stresses (R(2)=0.91, RMSE=8.6%, peak error=27%), while the two-material model was less accurate (R(2)=0.89, RMSE=9.6%, peak error=35%). The results showed that it is possible to automatically generate accurate finite element models of bones from CT data and that the strategy of material properties mapping has a significant influence on its accuracy.     

Deformation and stress distribution of the human foot after plantar ligaments release: A cadaveric study and finite element analysis

The majority of foot deformities are related to arch collapse or instability, especially the longitudinal arch. Although the relationship between the plantar fascia and arch height has been previously investigated, the stress distribution remains unclear. The aim of this study was to explore the role of the plantar ligaments in foot arch biomechanics. We constructed a geometrical detailed three-dimensional (3-D) finite element (FE) model of the human foot and ankle from computer tomography images. The model comprised the majority of joints in the foot as well as bone segments, major ligaments, and plantar soft tissue. Release of the plantar fascia and other ligaments was simulated to evaluate the corresponding biomechanical effects on load distribution of the bony and ligamentous structures. These intrinsic ligaments of the foot arch were sectioned to simulate different pathologic situations of injury to the plantar ligaments, and to explore bone segment displacement and stress distribution. The validity of the 3-D FE model was verified by comparing results with experimentally measured data via the displacement and von Mise stress of each bone segment. Plantar fascia release decreased arch height, but did not cause total collapse of the foot arch. The longitudinal foot arch was lost when all the four major plantar ligaments were sectioned simultaneously. Plantar fascia release was compromised by increased strain applied to the plantar ligaments and intensified stress in the midfoot and metatarsal bones. Load redistribution among the centralized metatarsal bones and focal stress relief at the calcaneal insertion were predicted. The 3-D FE model indicated that plantar fascia release may provide relief of focal stress and associated heel pain. However, these operative procedures may pose a risk to arch stability and clinically may produce dorsolateral midfoot pain. The initial strategy for treating plantar fasciitis should be non-operative.     

Correction to: Prediction of femoral strength using 3D finite element models reconstructed from DXA images: validation against experiments

Biomechanics and Modeling in Mechanobiology - In the original publication of the article, Fig.&nbsp;3 and Tables&nbsp;2, 4 and 5 were published with errors. The issue was caused by an error...     

Fracture prediction for the proximal femur using finite element models: Part I--Linear analysis.

Over 90 percent of the more than 250,000 hip fractures that occur annually in the United States are the result of falls from standing height. Despite this, the stresses associated with femoral fracture from a fall have not been investigated previously. Our objectives were to use three-dimensional finite element models of the proximal femur (with geometries and material properties based directly on quantitative computed tomography) to compare predicted stress distributions for one-legged stance and for a fall to the lateral greater trochanter. We also wished to test the correspondence between model predictions and in vitro strain gage data and failure loads for cadaveric femora subjected to these loading conditions. An additional goal was to use the model predictions to compare the sensitivity of several imaging sites in the proximal femur which are used for the in vivo prediction of hip fracture risk. In this first of two parts, linear finite element models of two unpaired human cadaveric femora were generated. In Part II, the models were extended to include nonlinear material properties for the cortical and trabecular bone. While there was poor correspondence between strain gage data and model predictions, there was excellent agreement between the in vitro failure data and the linear model, especially using a von Mises effective strain failure criterion. Both the onset of structural yielding (within 22 and 4 percent) and the load at fracture (within 8 and 5 percent) were predicted accurately for the two femora tested. For the simulation of one-legged stance, the peak stresses occurred in the primary compressive trabeculae of the subcapital region.(ABSTRACT TRUNCATED AT 250 WORDS)     

From high-resolution CT data to finite element models: development of an integrated modular framework

This work introduces a novel method of automating the process of patient-specific finite element (FE) model development using a mapped mesh technique. The objective is to map a predefined mesh (template) of high quality directly onto a new bony surface (target) definition, thereby yielding a similar mesh with minimal user interaction. To bring the template mesh into correspondence with the target surface, a deformable registration technique based on the FE method has been adopted. The procedure has been made hierarchical allowing several levels of mesh refinement to be used, thus reducing the time required to achieve a solution. Our initial efforts have focused on the phalanx bones of the human hand. Mesh quality metrics, such as element volume and distortion were evaluated. Furthermore, the distance between the target surface and the final mapped mesh were measured. The results have satisfactorily proven the applicability of the proposed method.     

Discretization error when using finite element models: Analysis and evaluation of an underestimated problem

Mesh convergence tests are often insufficiently performed in finite element analyses. There are many parameters which may have an effect on the mesh convergence behavior. The aim of this study was to identify the influence of different parameters on the mesh convergence behavior.For this purpose we used a simplified axis-symmetrical model of a single pedicle screw flank with surrounding bone to simulate a pull-out test. In parameter studies, the flank radii and the contact conditions at the bone–screw interface were varied. These parameter studies were carried out using an implicit and explicit solver. Thereby, the convergence criteria and the number of the substeps for the implicit nonlinear iteration process as well as the velocity and the material density for the explicit approach were considered.The mesh convergence behavior was influenced by varying the flank radii and the contact conditions. The implicit calculations led to a reaction force, which converged rapidly to a certain value with increasing mesh density, whereas the maximum von-Mises stress showed substantial convergence problems. The number of substeps and the convergence criteria of the iteration process strongly influenced the implicit solutions. In contrast, the maximum von-Mises stresses resulting from explicit calculations converged to a certain value after only a few refinement steps. Different pull-out velocities substantially affected the mesh convergence behavior, while the material density showed only a negligible influence.The results indicated the need to perform an appropriate mesh convergence test when using finite element methods. We were able to show that different parameters strongly influence the mesh convergence behavior and we demonstrated that convergence tests do not always lead to a satisfactory or acceptable solution.     

Fracture prediction for the proximal femur using finite element models: Part II--Nonlinear analysis.

In Part I we reported the results of linear finite element models of the proximal femur generated using geometric and constitutive data collected with quantitative computed tomography. These models demonstrated excellent agreement with in vitro studies when used to predict ultimate failure loads. In Part II, we report our extension of those finite element models to include nonlinear behavior of the trabecular and cortical bone. A highly nonlinear material law, originally designed for representing concrete, was used for trabecular bone, while a bilinear material law was used for cortical bone. We found excellent agreement between the model predictions and in vitro fracture data for both the onset of bone yielding and bone fracture. For bone yielding, the model predictions were within 2 percent for a load which simulated one-legged stance and 1 percent for a load which simulated a fall. For bone fracture, the model predictions were within 1 percent and 17 percent, respectively. The models also demonstrated different fracture mechanisms for the two different loading configurations. For one-legged stance, failure within the primary compressive trabeculae at the subcapital region occurred first, leading to load transfer and, ultimately, failure of the surrounding cortical shell. However, for a fall, failure of the cortical and trabecular bone occurred simultaneously within the intertrochanteric region. These results support our previous findings that the strength of the subcapital region is primarily due to trabecular bone whereas the strength of the intertrochanteric region is primarily due to cortical bone.