Research Papers

A Three-Dimensional Ankle Kinetostatic Model to Simulate Loaded and Unloaded Joint Motion

[+] Author and Article Information
Margherita Forlani

DIN-Department of Industrial Engineering,
Health Sciences and Technologies,
Interdepartmental Centre for Industrial
Research (HST-ICIR),
University of Bologna,
Bologna 40136, Italy
e-mail: margherita.forlani2@unibo.it

Nicola Sancisi

DIN-Department of Industrial Engineering,
Health Sciences and Technologies,
Interdepartmental Centre for Industrial
Research (HST-ICIR),
University of Bologna,
Bologna 40136, Italy
e-mail: nicola.sancisi@unibo.it

Vincenzo Parenti-Castelli

DIN-Department of Industrial Engineering,
Health Sciences and Technologies,
Interdepartmental Centre for Industrial
Research (HST-ICIR),
University of Bologna,
Bologna 40136, Italy
e-mail: vincenzo.parenti@unibo.it

Manuscript received March 5, 2014; final manuscript received February 6, 2015; published online March 25, 2015. Assoc. Editor: Tammy L. Haut Donahue.

J Biomech Eng 137(6), 061005 (Jun 01, 2015) (12 pages) Paper No: BIO-14-1105; doi: 10.1115/1.4029978 History: Received March 05, 2014; Revised February 06, 2015; Online March 25, 2015

A kinetostatic model able to replicate both the natural unloaded motion of the tibiotalar (or ankle) joint and the joint behavior under external loads is presented. The model is developed as the second step of a sequential procedure, which allows the definition of a kinetostatic model as a generalization of a kinematic model of the joint defined at the first step. Specifically, this kinematic model taken as the starting point of the definition procedure is a parallel spatial mechanism which replicates the ankle unloaded motion. It features two rigid bodies (representing the tibia–fibula and the talus–calcaneus complexes) interconnected by five rigid binary links, that mimic three articular contacts and two nearly isometric fibers (IFs) of the tibiocalcaneal ligament (TiCaL) and calcaneofibular ligament (CaFiL). In the kinetostatic model, the five links are considered as compliant; moreover, further elastic structures are added to represent all the main ankle passive structures of the joint. Thanks to this definition procedure, the kinetostatic model still replicates the ankle unloaded motion with the same accuracy as the kinematic model. In addition, the model can replicate the behavior of the joint when external loads are applied. Finally, the structures that guide these motions are consistent with the anatomical evidence. The parameters of the model are identified for two specimens from both subject-specific and published data. Loads are then applied to the model in order to simulate two common clinical tests. The model-predicted ankle motion shows good agreement with results from the literature.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Corazza, F., O’Connor, J. J., Leardini, A., and Parenti-Castelli, V., 2003, “Ligament Fibre Recruitment and Forces for the Anterior Drawer Test at the Human Ankle Joint,” J. Biomech., 36(3), pp. 363–372. [CrossRef] [PubMed]
Delp, S., Anderson, F., Arnold, A., Loan, P., Habib, A., John, C., Guendelman, E., and Thelem, D., 2007, “Opensim: Open-Source Software to Create and Analyze Dynamic Simulations of Movement,” IEEE Trans. Biomed. Eng., 54(11), pp. 1940–1950. [CrossRef] [PubMed]
Siegler, S., Block, J., and Schneck, C. D., 1988, “The Mechanical Characteristics of the Collateral Ligaments of the Human Ankle Joint,” Foot Ankle, 8(5), pp. 234–242. [CrossRef] [PubMed]
Mommersteeg, T. J. A., Blankevoort, L., Huiskes, R., Kooloos, J. G. M., Kauer, J. M. G., and Hendriks, J., 1995, “The Effect of Variable Relative Insertion Orientation of Human Knee Bone-Ligament-Bone Complexes on the Tensile Stiffness,” J. Biomech., 28(6), pp. 745–752. [CrossRef] [PubMed]
Guess, T. M., Thiagarajan, G., Kia, M., and Mishra, M., 2010, “A Subject Specific Multibody Model of the Knee With Menisci,” Med. Eng. Phys., 32(5), pp. 505–515. [CrossRef] [PubMed]
Liacouras, P. C., and Wayne, J. S., 2007, “Computational Modeling to Predict Mechanical Function of Joints: Application to the Lower Leg With Simulation of Two Cadaver Studies,” ASME J. Biomech. Eng., 129(6), pp. 811–817. [CrossRef]
Imhauser, C., Siegler, S., Udupa, J., and Toy, J., 2008, “Subject-Specific Models of the Hindfoot Reveal a Relationship Between Morphology and Passive Mechanical Properties,” J. Biomech., 41(6), pp. 1341–1349. [CrossRef] [PubMed]
Wei, F., Hunley, S., Powell, J., and Haut, R., 2011, “Development and Validation of a Computational Model to Study the Effect of Foot Constraint on Ankle Injury Due to External Rotation,” Ann. Biomed. Eng., 39(2), pp. 756–765. [CrossRef] [PubMed]
Shelburne, K. B., Torry, M. R., and Pandy, M. G., 2005, “Muscle, Ligament, and Joint-Contact Forces at the Knee During Walking,” Med. Sci. Sports Exercise, 37(11), pp. 1948–1956. [CrossRef]
Franci, R., Parenti-Castelli, V., and Sancisi, N., 2009, “A Three-Step Procedure for the Modelling of Human Diarthrodial Joints,” Int. J. Mech. Control, 10(1), pp. 3–12.
Sancisi, N., and Parenti-Castelli, V., 2011, “A Sequentially-Defined Stiffness Model of the Knee,” Mech. Mach. Theory, 46(12), pp. 1920–1928. [CrossRef]
Forlani, M., Baldisserri, B., Sancisi, N., and Parenti-Castelli, V., 2011, “On the Modelling of the Ankle Motion Under Static Loads by a Sequential Procedure: Model Definition and Preliminary Results,” Atti del XX CONGRESSO dell’Associazione Italiana di Meccanica Teorica e Applicata, pp. 1–10.
Leardini, A., O’Connor, J. J., Catani, F., and Giannini, S., 1999, “A Geometric Model of the Human Ankle Joint,” J. Biomech., 32(6), pp. 585–591. [CrossRef] [PubMed]
Di Gregorio, R., Parenti-Castelli, V., O’Connor, J. J., and Leardini, A., 2007, “Mathematical Models of Passive Motion at the Human Ankle Joint by Equivalent Spatial Parallel Mechanisms,” Med. Biol. Eng. Comput., 45(3), pp. 305–313. [CrossRef] [PubMed]
Franci, R., and Parenti-Castelli, V., 2011, “A One-Degree-of-Freedom Spherical Wrist for the Modelling of Passive Motion of the Human Ankle Joint,” Interdisciplinary Applications of Kinematics, Springer, Dordrecht, The Netherlands, pp. 183–195.
Franci, R., Parenti-Castelli, V., Belevedere, C., and Leardini, A., 2009, “A New One-DOF Fully Parallel Mechanism for Modelling Passive Motion at the Human Tibiotalar Joint,” J. Biomech., 42(10), pp. 1403–1408. [CrossRef] [PubMed]
Sancisi, N., Baldisserri, B., Parenti-Castelli, V., Belvedere, C., and Leardini, A., 2014, “One-Degree-of-Freedom Spherical Model for the Passive Motion of the Human Ankle Joint,” Med. Biol. Eng. Comput., 52(4), pp. 363–373. [CrossRef] [PubMed]
Leardini, A., O’Connor, J. J., Catani, F., and Giannini, S., 2000, “The Role of the Passive Structures in the Mobility and Stability of the Human Ankle Joint: A Literature Review,” Foot Ankle Int., 21(7), pp. 602–615. [PubMed]
Sancisi, N., and Parenti-Castelli, V., 2010, “A 1-DOF Parallel Spherical Wrist for the Modelling of the Knee Passive Motion,” Mech. Mach. Theory, 45(4), pp. 658–665. [CrossRef]
Parenti-Castelli, V., and Sancisi, N., 2013, “Synthesis of Spatial Mechanisms to Model Human Joints,” 21st Century Kinematics, M.McCarthy, ed., Springer, London, pp. 49–84.
Wu, G., Siegler, S., Allard, P., Kirtley, C., Leardini, A., Rosenbaum, D., Whittle, M., DLima, D. D., Cristofolini, L., Witte, H., Schmid, O., and Stokes, I., 2002, “ISB Recommendation on Definitions of Joint Coordinate System of Various Joints for the Reporting of Human Joint Motion—Part I: Ankle, Hip, and Spine,” J. Biomech., 35(4), pp. 543–548. [CrossRef] [PubMed]
Pankovich, A. M., and Shivaram, M. S., 1989, “Anatomical Basis of Variability in Injuries of the Medial Malleolus and the Deltoid Ligament. Part 1: Anatomical Studies,” Acta Orthop., 50(2), pp. 217–223. [CrossRef]
Hefzy, M. S., and Grood, E. S., 1983, “An Analytical Technique for Modeling Knee Joint Stiffness—Part II: Ligamentous Geometric Nonlinearities,” ASME J. Biomech. Eng., 105(2), pp. 145–153. [CrossRef]
Zavatsky, A. B., and O’Connor, J. J., 1992, “A Model of Human Ligaments in the Sagittal Plane. Part 2: Fibre Recruitment Under Load,” Proc. Inst. Mech. Eng., Part H, 206(3), pp. 135–145. [CrossRef]
Sindel, M., Demir, S., Yildirim, A., and Ucar, Y., 1998, “Anatomy of the Lateral Ankle Ligaments,” Turk. J. Med. Sci., 28(1), pp. 53–56.
Taser, F., Shafiq, Q., and Ebraheim, N. A., 2006, “Anatomy of Lateral Ankle Ligaments and Their Relationship to Bony Landmarks,” Surg. Radiol. Anat., 28(4), pp. 391–397. [CrossRef] [PubMed]
Primal Pictures, 2003, Primal 3D Interactive Series, Primal Pictures Ltd., London.
Netter, F. H., 2010, Atlas of Human Anatomy, Saunders Elsevier, Philadelphia.
Siegler, S., Udupa, J. K., Ringleb, S. I., Imhauser, C. W., Hirsch, B. E., Odhner, D., Saha, P. K., Okereke, E., and Roach, N., 2005, “Mechanics of the Ankle and Subtalar Joints Revealed Through a 3D Quasi-Static Stress MRI Technique,” J. Biomech., 38(3), pp. 567–578. [CrossRef] [PubMed]
Leardini, A., O’Connor, J. J., Catani, F., and Giannini, S., 1999, “Kinematics of the Human Ankle Complex in Passive Flexion: A Single Degree of Freedom System,” J. Biomech., 32(2), pp. 111–118. [CrossRef] [PubMed]
Moissenet, F., Chèze, L., and Dumas, R., 2012, “Anatomical Kinematic Constraints: Consequences on Musculo-Tendon Forces and Joint Reactions,” Multibody Syst. Dyn., 28(1–2), pp. 125–141. [CrossRef]
Jend, H. H., Ney, R., and Heller, M., 1985, “Evaluation of Tibiofibular Motion Under Load Conditions by Computed Tomography,” J. Orthop. Res., 3(4), pp. 418–423. [CrossRef] [PubMed]
Beumer, A., Valstar, E., Garling, E., Niesing, R., Ranstam, J., Löfvenberg, R., and Swierstra, B., 2003, “Kinematics of the Distal Tibiofibular Syndesmosis: Radiostereometry in 11 Normal Ankles,” Acta Orthop., 74(3), pp. 337–343. [CrossRef]
Leardini, A., Stagni, R., and O’Connor, J. J., 2001, “Mobility of the Subtalar Joint in the Intact Ankle Complex,” J. Biomech., 34(6), pp. 805–809. [CrossRef] [PubMed]
Beimers, L., Maria Tuijthof, G. J., Blankevoort, L., Jonges, R., Maas, M., and van Dijk, C. N., 2008, “In-Vivo Range of Motion of the Subtalar Joint Using Computed Tomography,” J. Biomech., 41(7), pp. 1390–1397. [CrossRef] [PubMed]
Bahr, R., Pena, F., Shine, J., Lew, W. D., Lindquist, C., Tyrdal, S., and Engebretsen, L., 1997, “Mechanics of the Anterior Drawer and Talar Tilt Tests: A Cadaveric Study of Lateral Ligament Injuries of the Ankle,” Acta Orthop., 68(5), pp. 435–441. [CrossRef]
Lapointe, S. J., Siegler, S., Hillstrom, H., Nobilini, R. R., Mlodzienski, A., and Techner, L., 1997, “Changes in the Flexibility Characteristics of the Ankle Complex Due to Damage to the Lateral Collateral Ligaments: An in vitro and in vivo Study,” J. Orthop. Res., 15(3), pp. 331–341. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Definition of the kinetostatic model: (a) and (b) represent the definition of the kinematic model (first step of the sequential approach); (c) and (d) show the generalization of this model to obtain the kinetostatic one (second step of the procedure). In particular, the main structures guiding the passive motion are shown in (a): the three bone contacts and the TiCaL and CaFiL. A schematic representation of the kinematic model is reported in (b), where the three articular contacts are modeled as sphere-on-sphere contacts, and the two IFs of TiCaL and CaFiL are modeled by two rigid links. The IFs are considered as compliant (black lines) in (c) and elastic fibers are added (gray lines) to model other ligaments. The 26 fibers identified on the 3D bone surfaces obtained from CT scans of one representative ankle specimen are shown in (d): in particular, the IFs at the TiCaL (IF1) and CaFiL (IF2), and the CFs at the lateral malleolus (C1), at the internal region of the inferior surface of the distal tibia (C2), and at the medial malleolus (C3) are highlighted: these fibers correspond to the five constraints in (b).

Grahic Jump Location
Fig. 2

The joint coordinate system used for description of the ankle motion. The reference frames Sf of the tibia–fibula and Sc of the talus–calcaneus complexes are defined according to the ISB recommendations [21]. The six components of the relative motion between the two bone complexes are described as rotations about and displacements along the axes e1, e2, e3 of the joint coordinate system: axis e1 is coincident to zf, axis e3 is coincident to xc, and e2 is a floating axis perpendicular to e1 and e3.

Grahic Jump Location
Fig. 3

Pose components of the unloaded motion for the specimens A (left) and B (right). Rotations are: D/P = dorsiflexion(+)/plantarflexion(−), Int/Ext = internal(+)/external(−), and Inv/Ev = inversion(+)/eversion(−). Displacements are: L/M = lateral(+)/medial(−), P/D = proximal(+)/distal(−), and A/P = anterior(+)/posterior(−). Crosses are experimental data, while continuous lines represent the model motion.

Grahic Jump Location
Fig. 4

Results of the two kinetostatic models determined for the specimens A (circles) and B (squares) compared with the mean values (crosses) and SDs (lines) of the reference experimental data [29]. The flexion is fixed at the mean value reported in the reference paper, and the other five motion components are rotations and displacements from the neutral pose. In the tables, the absolute and relative (with respect to the corresponding SD) differences between the computed and the mean reference values of the pose components are reported.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In