Numerical simulation of strain-adaptive bone remodelling in the ankle joint
© Bouguecha et al; licensee BioMed Central Ltd. 2011
Received: 11 April 2011
Accepted: 5 July 2011
Published: 5 July 2011
The use of artificial endoprostheses has become a routine procedure for knee and hip joints while ankle arthritis has traditionally been treated by means of arthrodesis. Due to its advantages, the implantation of endoprostheses is constantly increasing. While finite element analyses (FEA) of strain-adaptive bone remodelling have been carried out for the hip joint in previous studies, to our knowledge there are no investigations that have considered remodelling processes of the ankle joint. In order to evaluate and optimise new generation implants of the ankle joint, as well as to gain additional knowledge regarding the biomechanics, strain-adaptive bone remodelling has been calculated separately for the tibia and the talus after providing them with an implant.
FE models of the bone-implant assembly for both the tibia and the talus have been developed. Bone characteristics such as the density distribution have been applied corresponding to CT scans. A force of 5,200 N, which corresponds to the compression force during normal walking of a person with a weight of 100 kg according to Stauffer et al., has been used in the simulation. The bone adaptation law, previously developed by our research team, has been used for the calculation of the remodelling processes.
A total bone mass loss of 2% in the tibia and 13% in the talus was calculated. The greater decline of density in the talus is due to its smaller size compared to the relatively large implant dimensions causing remodelling processes in the whole bone tissue. In the tibia, bone remodelling processes are only calculated in areas adjacent to the implant. Thus, a smaller bone mass loss than in the talus can be expected. There is a high agreement between the simulation results in the distal tibia and the literature regarding.
In this study, strain-adaptive bone remodelling processes are simulated using the FE method. The results contribute to a better understanding of the biomechanical behaviour of the ankle joint and hence are useful for the optimisation of the implant geometry in the future.
Arthrodesis is a preferred and most used operative therapy in advanced symptomatic arthrosis of the ankle . The objective is pain relief with a stable osseous fusion. Disadvantages are the loss of movement in the joint accompanied by an increased mobility in the transversal joint (Chopart joint) as a compensation reaction. This can lead to secondary overloading and arthrosis. The activity of the patient is significantly limited by arthrodesis. Further disadvantages are the risk of non fusion as well as the long rehabilitation time [2–6].
The advantage of TAR is the maintenance of the ankle motion. Therefore, TAR should reduce the load on the adjacent joints and prevent the development of secondary arthrosis . However, the development of TAR is lagged behind that of the hip and knee . Difficulties are the smaller joint size  and the comparatively higher stresses applied to the ankle joint resulting from higher compression forces [10–14] and torques [15, 16]. Furthermore, patients requiring TAR are generally younger and therefore more active [1, 2].
Because of their full congruency without restriction of rotational motion, the modern three-component designs have mechanical and kinematical advantages compared to two-component designs. The rotational stresses at the bone-prosthesis interface are reduced by this development. Independent from the exact design, the 5-year survival rate of the third generation endoprostheses is up to 90% .
However, the results still indicate inferior patient satisfaction, compared to knee and hip arthroplasty, as mentioned previously. Aseptic loosening is the main reason for early failure of TAR . It results from stress shielding caused by the implant. The difference in mechanical properties of the bone tissue and the prosthesis and the changed bone loading condition due to TAR promote stress shielding. Bone tissue is in a permanent state of resorption and formation. Bone acts like a technical controller on load changes . Set point is the elastic deformation, which is regulated by a variation of the amount of bone tissue. Accordingly, a load increase leads to an increase in the elastic deformation and to the formation of new bone tissue. Contrarily, a load reduction leads to bone resorption . The implantation of an artificial joint leads to a change in biomechanics. This can lead to intense bone resorption around the implant and can cause the migration or loosening of the implant.
In order to increase long-term stability and to avoid implant loosening, it is important to estimate the adaptive bone remodelling prior to implantation in preclinical studies.
In our previous studies, bone remodelling processes were successfully simulated for canine  and human [25, 26] femora after total hip replacement (THR). The influence of the boundary conditions were also investigated . Furthermore, the significance of implant materials [28, 29] and implantation techniques [30, 31] of the remodelling processes in the femur were analysed. Recently, bone remodelling in the acetabulum was simulated as well .
FE simulations concerning TARs are quite rare. The following section provides an overview of simulations at the ankle joint found in the literature.
In the middle of the 1980s Falsig et al.  employed a 3D finite element stress analysis to calculate the stresses in the distal tibia. The goal was to reduce the cement-bone-interface and the stresses in the bone by varying the shape and the material of the tibial component. Prostheses with metal backing are advantageous compared to PE components without metal backing. Furthermore, long stems for the implant fixation proved unfavourable as well.
The stresses and wear of the PE inlay of TARs was examined by FEA from the research group of McIff [40, 41]. In the first study a two component design (AGILITY™) and a three component design (S.T.A.R.®) were compared. For the AGILITY™ system the loss of congruency between the PE inlay and the talar component resulted in high contact and internal stresses due to point and line contact. In a second study the design of the inlay and of the talar component was varied. Apparent insignificant geometrical specifications had an important influence on the contact mechanism.
A recent study of Espinosa et al.  confirmed the result of the first examination of McIff et al. . The research group around Galik  examined the PE wear and varied the thickness of the PE inlay.
More complex models were established by Reggiani et al  and Anderson et al. . Reggiani calculated the contact pressure on the components. They computed average values of 6.4 MPa for the tibial component and 10.3 MPa for the talar component .
Anderson et al.  developed a complex FE model for the simulation of the contact pressure distribution in the healthy ankle joint. Recently, Anderson et al.  performed a FEA with a new implant for defect resurfacing of the talus according clinical feasibility.
The aim of this study is to calculate the strain-adaptive bone remodelling after TAR by means of the finite element method (FEM).
First, FE models of the tibia and the talus should be built. Therefore, models of the intact bone and the bone provided with the implant are required. In order to obtain geometrical data, a left ankle joint from Sawbones (Malmö/Sweden) was optically measured with the 3D coordinate measurement system ATOS II (GOM mbH, Braunschweig/Germany). The bones from Sawbones were used to obtain a standard geometry model of the ankle joint and to eliminate patient-individual anatomical characteristics.
The digitalised STL data of the tibia and the talus were reconstructed by means of the 3D software Mimics (Materialise, Leuven/Belgium). CT scans (Philips Brilliance CT 64) of 20 cadaver ankle joints were performed to define the density distribution throughout the bone on the basis of the gray values. Accordingly, the tibial bone model was divided into three different density areas: cortical bone in the outer layer (1.7 g/cm3), cancellous bone in the inner layer (0.5 g/cm3) and an interface layer between cortical and cancellous bone (1.0 g/cm3). The talus was divided into only two different density areas namely cortical bone (1.7 g/cm3) and cancellous bone (0.5 g/cm3). An interface layer has not been observed on the CT scans.
In order to build the FE models of the bone-implant assembly, the prosthesis was integrated into the previously described bone models via the pre-processor software HyperMesh (Altair Engineering GmbH, Böblingen/Germany). The cementless three component S.T.A.R.® system was used, which is the most-implanted ankle joint system in Europe  (Figure 2).
The alignment of the implant was defined based on surgical experience.
For an optimal force transmission from the implant to the bone and for the modelling of the cementless fixation of the implant, a consistent meshing method has been used to realise these requirements. In addition, the models of the tibia and talus in the physiological state without implant already provide the implant geometry. Errors resulting from density transmission inaccuracies are excluded a priori.
The prosthesis is made of a cobalt-chromium-molybdenium alloy with titanium plasma spray coating. In the FE modelling a homogeneous, isotropic material law was used for the prosthesis (E = 210,000 N/mm2, ν = 0.3). The coating was considered using proper friction coefficients between the bone and the implant surface.
Loads and boundary conditions
The load was taken from the study of Stauffer et al. , in which a load of 5.2 times the body weight had been calculated for the normal ankle joint during the stance phase of gait. A person of 100 kg weight was assumed. Consequently, a static force F of 5,200 N was applied to the tibia and the talus by spreading it equally on ten nodes.
FE simulation of the bone remodelling
Herein, ε represents the strain vector and σ T the transposed stress vector.
The used bone adaptation law is a modification of the law of Huiskes . A limitation of the bone formation rate has been introduced. It is assumed, that severe overloading causes lyses of bone tissue and no unlimited bone formation.
On the basis of a constant element volume, the density values of each element were recalculated in every increment. The apparent bone density evolution rate was calculated depending on the stimulus for the bone remodelling. After the convergence of the calculation at a constant density value, a balance is adjusted between resorption and forming of bone tissue. The remodelling of bone tissue can be indicated according to the density distribution in the tibia and the talus at the end of the simulation.
Bone remodelling in the tibia
To the left image of Figure 7, the density distribution in the physiological state without implant is shown, whereas on the two right images, the density distribution of the prosthetic ankle is represented.
An increase in density is to be expected above both anchoring bolts, while a decrease in density is calculated centrally above the tibial plate and in the medial malleolus.
Bone remodelling in the talus
The bone mass loss is observed in the regions beneath the implant. A considerable increase in density is expected according to the FE analysis beneath the central fin of the component.
The number of TARs is constantly increasing . Thomas and Daniels  pointed out, that beside good clinical results with modern three-component systems, no statements can be made about the long-term stability of the implants.
This makes it even more important to simulate the bone remodelling in order to estimate the clinical performance of the endoprostheses prior to implantation.
In this work, a 3D FE simulation of the bone remodelling processes in the tibia and the talus was performed. The bone adaptation law used in the study has been calibrated for the hip joint via DEXA investigations in recent studies . The calculated bone mass loss was in good agreement to the conclusions of the DEXA analysis. The application of the adaptation law to the ankle joint is considered to be adequate.
Nevertheless, every FE simulation has to be validated in order to proof the accuracy of the calculation result. Therefore, further DEXA investigations of the ankle joint after TAR are still required.
The simulation was done using homogeneous elastic properties for the bones, although bone tissue consists of cortical and cancellous bone with complex trabecular architecture. This simplification was made due to the observation on the CT scans, where the patterns of the cancellous bone were not clearly visible. A nearly homogeneous density distribution for the cancellous and cortical bone in the talus has been observed. For the tibia, a third density area was added between cortical and cancellous bone to reflect the observations on the CT scans.
Finally, first results of the 3D FE simulation concerning the strain-adaptive bone remodelling processes in the whole ankle joint were presented in this study.
Good agreement was achieved between the simulation result in the tibia and the changes in density distribution after TAR described in the literature. An increase in density has been observed above the anchoring bolts of the implant and a density decrease was observed centrally above the tibial plate . This confirms the previously described areas of strain-adaptive bone remodelling. The fixation of the tibial component is achieved with two anchoring bolts. This results in force transmission from the two bars into the bone. This may lead to stress shielding between the anchoring bolts and above the tibial plate which has also been observed from Hintermann .
In the simulated model the fibula was disregarded because there is no direct fixation of an implant at the fibula. Furthermore, it does not contribute to the load distribution from the tibia to the talus. The fibula is irrelevant for bone remodelling processes in the ankle joint and hence neglecting its role is justified. Other researchers who examined the ankle joint by FEM [42–48] disregarded the fibula as well.
Another simplification was the use of only one static load case. In reality more complex loading conditions can be expected. Moreover, the muscle forces acting on the ankle joint were disregarded due to the lack of experimental study. For future investigations, force patterns, derived from multibody simulations, should be incorporated into the FE calculation to consider the whole gait cycle and in order to examine the influence of the muscle forces. Furthermore migration processes of the implants can be calculated by coupling the multibody simulation and the FE calculation.
Further investigations are planned regarding the previously described limitations.
To our knowledge the first finite element simulation of strain-adaptive bone remodelling of the ankle joint is presented in this study. High agreement between the simulation and the clinical results were shown in a qualitative comparison between the calculated bone remodelling in the distal tibia and the reported radiographic changes in literature . On the basis of the calculated density distributions, the design of total ankle prostheses can be evaluated and optimised.
This study was supported by the Leibniz Universität Hannover in the project "Experimental and numerical investigations for the optimisation of total ankle endoprostheses". We wish to thank the Leibniz Universität Hannover for the financial support. Furthermore, we thank the Small Animal Clinic of the University of Veterinary Medicine Hannover for carrying out the CT scans. We would also like to thank Small Bone Innovations (Donaueschingen/Germany) for providing the implants and technical drawings.
- Rao S, Ellis SJ, Deland JT, Hillstrom H: Nonmedicinal therapy in the management of ankle arthritis. Current Opinion in Rheumatology 2010, 22: 223–228. 10.1097/BOR.0b013e328335fcebView ArticleGoogle Scholar
- Hintermann B: Endoprothetik des Sprunggelenks. Historischer Überblick, aktuelle Therapiekonzepte und Entwicklungen. Wien New York, Springer Verlag; 2005.Google Scholar
- Bauer G, Eberhardt O, Rosenbaum D, Claes L: Total ankle replacement. Review and critical analysis of the current status. Foot and Ankle Surgery 1996, 2: 119–126. 10.1046/j.1460-9584.1996.00024.xView ArticleGoogle Scholar
- Doets HC, van Middelkoop M, Houdijk H, Nelissen RGHH, Veeger HEJD: Gait analysis after successful mobile bearing total ankle replacement. Foot and Ankle International 2007, 28: 313–322. 10.3113/FAI.2007.0313View ArticleGoogle Scholar
- Kofoed H, Stürup J: Comparison of ankle arthroplasty and arthrodesis. A prospective series with long-term follow-up. The Foot 1994, 4: 6–9. 10.1016/0958-2592(94)90003-5View ArticleGoogle Scholar
- McGuire MR, Kyle RF, Gustilo RB, Premer RF: Comparative analysis of ankle arthroplasty versus ankle arthrodesis. Clinical Orthopaedics and Related Research 1988, 226: 174–181.Google Scholar
- Gougoulias NE, Khanna A, Maffulli N: History and evolution in total ankle arthroplasty. British Medical Bulletin 2009, 89: 111–151.View ArticleGoogle Scholar
- Vickerstaff JA, Miles AW, Cunningham JL: A brief history of total ankle replacement and a review of the current status. Medical Engineering and Physics 2007, 29: 1056–1064. 10.1016/j.medengphy.2006.11.009View ArticleGoogle Scholar
- Leardini A, O'Connor JJ, Catani F, Giannini S: A geometric model of the human ankle joint. Journal of Biomechanics 1999, 32: 585–591. 10.1016/S0021-9290(99)00022-6View ArticleGoogle Scholar
- Stauffer RN, Chao EYS, Brewster RC: Force and motion analysis of the normal, diseased, and prosthetic ankle joint. Orthopaedics and Related Research 1977, 127: 189–196.Google Scholar
- Burdett RG: Forces predicted at the ankle during running. Medicine and Science in Sports and Exercise 1982, 14: 308–316. 10.1249/00005768-198204000-00010View ArticleGoogle Scholar
- Procter P, Paul JP: Ankle joint biomechanics. Journal of Biomechanics 1982, 15: 627–634. 10.1016/0021-9290(82)90017-3View ArticleGoogle Scholar
- Röhrle H, Scholten R, Sigolotto C, Sollbach W, Kellner H: Joint forces in the human pelvis-leg skeleton during walking. Journal of Biomechanics 1984, 17: 409–424. 10.1016/0021-9290(84)90033-2View ArticleGoogle Scholar
- Seireg A, Arvikar RJ: The prediction of muscular load sharing and joint forces in the lower extremities during walking. Journal of Biomechanics 1975, 8: 89–102. 10.1016/0021-9290(75)90089-5View ArticleGoogle Scholar
- Benedetti MG, Catani F, Leardini A, Pignotti E, Giannini S: Data management in gait analysis for clinical applications. Clinical Biomechanics 1998, 13: 204–215. 10.1016/S0268-0033(97)00041-7View ArticleGoogle Scholar
- Leardini A: Geometry and mechanics of the human ankle complex and ankle prosthesis design. Clinical Biomechanics 2001, 16: 706–709. 10.1016/S0268-0033(01)00022-5View ArticleGoogle Scholar
- Gougoulias N, Khanna A, Maffulli N: How successful are current ankle replacements? A systematic review. Clinical Orthopaedics and Related Research 2010, 468: 199–208. 10.1007/s11999-009-0987-3View ArticleGoogle Scholar
- Haddad SL, Coetzee JC, Estok R, Fahrbach K, Banel D, Nalysnyk L: Intermediate and long-term outcomes of total ankle arthroplasty and ankle arthrodesis. A systematic review of the literature. Journal of Bone and Joint Surgery Am 2007, 89: 1899–1905. 10.2106/JBJS.F.01149View ArticleGoogle Scholar
- S.T.A.R. ® Ankle Implants and Instruments Guide, Small Bone Innovations 2011. [http://www.totalsmallbone.com/us/pdfs/MKT16011RevE.pdf]
- Fevang BTS, Lie SA, Havelin LI, Skredderstuen A, Brun JG, Furnes O: 257 ankle arthroplasties performed in Norway between 1994 to 2005. Acta Orthopaedica 2007, 78: 575–83. 10.1080/17453670710014257View ArticleGoogle Scholar
- Henricson A, Skoog A, Carlsson A: The Swedish ankle arthroplasty register: An analysis of 531 arthroplasties between 1993 and 2005. Acta Orthopaedica 2007, 78: 569–574. 10.1080/17453670710014248View ArticleGoogle Scholar
- Kummer B: Biomechanik. Form und Funktion des Bewegungsapparates. Köln, Deutscher Ärzte-Verlag; 2005.Google Scholar
- Wolff J: Das Gesetz der Transformation der Knochen. Berlin, Hirschwald; 1892.Google Scholar
- Bouguecha A: Numerische und experimentelle Untersuchungen zum beanspruchungsadaptiven Knochenumbau im periprothetischen caninen Femur. Leibniz Universität Hannover, Institut für Umformtechnik und Umformmaschinen; 2007.Google Scholar
- Behrens B-A, Bouguecha A: Femorale Belastungsabschirmung nach einer Hüftarthroplastik. 57. Jahrestagung der Norddeutschen Orthopädenvereinigung, 12.-14.06. Hamburg 2008.Google Scholar
- Behrens B-A, Wirth CJ, Windhagen H, Nolte I, Meyer-Lindenberg A, Bouguecha A: Numerical investigations of stress shielding in total hip prostheses. Proceedings of the Institution of Mechanical Engineers, Part H, Journal of Engineering in Medicine 2008, 222: 593–600.View ArticleGoogle Scholar
- Behrens B-A, Nolte I, Wefstaedt P, Stukenborg-Colsman C, Bouguecha A: Numerical investigations on the strain-adaptive bone remodelling in the periprosthetic femur: Influence of the boundary conditions. BioMedical Engineering OnLine 2009, 8: 1–9. 10.1186/1475-925X-8-1View ArticleGoogle Scholar
- Behrens B-A, Bouguecha A, Stukenborg-Colsman C, Pressel T, Wefstaedt P, Nolte I: Numerische Untersuchungen zum Einfluss des Implantatwerkstoffs auf den postoperativen Knochenumbau im Femur. Jahrestagung der Deutschen Gesellschaft für Biomaterialien (DGBM), 22.-24.11. Hannover 2007.Google Scholar
- Bouguecha A, Behrens B-A, Meyer-Lindenberg A, Wefstaedt P, Stukenborg-Colsman C, Nolte I: Strain-adaptive bone modelling: Influence of the implant material. 16th Congress of the European Society of Biomechanics, 06.-09.07. Luzern, Schweiz 2008.Google Scholar
- Behrens B-A, Bouguecha A, Nolte I, Stukenborg-Colsman C, Pressel T: Strain adaptive bone remodelling: Influence of the implantation technique. BIOMECH Medicine meets Engineering, 13.-15.06. Regensburg 2007.Google Scholar
- Pressel T, Bouguecha A, Rittmann P, Wirth CJ, Nolte I, Meyer-Lindenberg A, Behrens B-A: Einfluss der Implantationstechnik auf die femorale Belastungsverteilung. Gemeinsamer Kongress der Orthopädie und der Unfallchirurgie, Berlin 2006.Google Scholar
- Bouguecha A, Elgaly I, Stukenborg-Colsman C, Lerch M, Nolte I, Wefstaedt P, Matthias T, Behrens B-A: Numerical investigations of the strain-adaptive bone remodelling in the prosthetic pelvis. IFMBE Proceedings 2010, 29: 562–565. 10.1007/978-3-642-13039-7_141View ArticleGoogle Scholar
- Ebbecke B: Theoretische und algorithmische Konzepte zur Beschreibung des beanspruchungsadaptiven Knochenwachstums. PhD thesis. Leibniz Universität Hannover, Institut für Baumechanik und Numerische Mechanik; 2006.Google Scholar
- Huiskes R, van Rietbergen B: Preclinical testing of total hip stems. The effects of coating placement. Clinical Orthopaedics and Related Research 1995, 319: 64–76.Google Scholar
- Nackenhorst U: Numerical simulation of stress stimulated bone remodeling. Technische Mechanik 1997, 17: 31–40.Google Scholar
- Prendergast PJ: Review paper. Finite element models in tissue mechanics and orthopaedic implant design. Clinical Biomechanics 1997, 12: 343–366. 10.1016/S0268-0033(97)00018-1View ArticleGoogle Scholar
- Van Rietbergen B, Huiskes R, Weinans H, Sumner DR, Turner TM, Galante JO: ESB Research Award 1992. The mechanism of bone remodeling and resorption around press-fitted THA stems. Journal of Biomechanics 1993, 26: 369–382. 10.1016/0021-9290(93)90001-UView ArticleGoogle Scholar
- Weinans H, Huiskes R, van Rietbergen B, Sumner DR, Turner TM, Galante JO: Adaptive bone remodelling around bonded noncemented total hip arthroplasty: A comparison between animal experiments and computer simulation. Journal of Orthopaedic Research 1993, 11: 500–513. 10.1002/jor.1100110405View ArticleGoogle Scholar
- Falsig J, Hvid I, Jensen NC: Finite element stress analysis of some ankle joint prostheses. Clinical Biomechanics 1986, 1: 71–76. 10.1016/0268-0033(86)90078-1View ArticleGoogle Scholar
- Mclff TE, Horton GA, Saltzman CL, Brown TD: Finite element modeling of total ankle replacements for constraint and stress analysis. Proceedings of the Fourth World Congress in Biomechanics, Calgary, Canada 2002.Google Scholar
- Mclff TE: Design factors affecting the contact stress patterns in a contemporary mobile bearing total ankle replacement. Proceedings of the Fourth World Congress in Biomechanics, Calgary, Canada 2002.Google Scholar
- Espinosa N, Walti M, Favre P, Snedeker JG: Misalignment of total ankle components can induce high joint contact pressures. Journal of Bone and Joint Surgery 2010, 92A: 1179–1187.View ArticleGoogle Scholar
- Galik K, Miller MC, Smolinski PJ: Effect of polyethylene thickness on stresses of an ankle joint implant. Proceedings of the Fourth World Congress in Biomechanics, Calgary, Canada 2002.Google Scholar
- Reggiani B, Leardini A, Corazza F, Taylor M: Finite element analysis of a total ankle replacement during the stance phase of gait. Journal of Biomechanics 2006, 39: 1435–1443. 10.1016/j.jbiomech.2005.04.010View ArticleGoogle Scholar
- Anderson DD, Goldsworthy JK, Shivanna K, Grosland NM, Pedersen DR, Thomas TP, Tochigi Y, Marsh JL, Brown TD: Intra-articular contact stress distributions at the ankle throughout stance phase-patient-specific finite element analysis as a metric of degeneration propensity. Biomechanics and Modeling in Mechanobiology 2006, 5: 82–89. 10.1007/s10237-006-0025-2View ArticleGoogle Scholar
- Anderson DD, Tochigi Y, Rudert MJ, Vaseenon T, Brown TD, Amendola A: Effect of implantation accuracy on ankle contact mechanics with a metallic focal resurfacing implant. Journal of Bone and Joint Surgery 2010, 92A: 1490–1500.View ArticleGoogle Scholar
- Carter DR, Hayas WC: The compressive behaviour of bone as a twophase porous structure. Journal of Bone and Joint Surgery 1977, 59A: 954–962.Google Scholar
- Thomas RH, Daniels TR: Ankle arthritis. Journal of Bone and Joint Surgery 2003, 85A: 923–936.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.