- Research article
- Open Access
In vitro study of foot bone kinematics via a custom-made cadaveric gait simulator
Journal of Orthopaedic Surgery and Research volume 15, Article number: 346 (2020)
Quantifying detailed kinematics of the intrinsic foot bone during gait is crucial for understanding biomechanical functions of the foot complex musculoskeletal structure and making appropriate surgery decisions.
The purpose of this experiment is to measure bone kinematic of the normal foot in a gait cycle via a custom-made cadaveric gait simulator.
In this experiment, we used a custom-made 6 degrees of freedom (DOF) of robotic gait simulator simulating normal human gait to measure the 3-dimensional (3D) kinematics of tibia, calcaneus, cuboid, navicular, medial cuneiform, first metatarsal, and fifth metatarsal through six cadaveric feet.
The results showed that the kinematic of the intrinsic foot bones in the stance phase of the gait was successfully quantified using a custom-made robotic gait simulator. During walking stance, the joints in the medial column of foot had less movement than those in the lateral column. And during the later portion of stance, no rotational cease was observed in the movement between navicular and cuboid, calcaneocuboid joint, or cuneonavicular joint.
This study described foot bone motion using a biomechanically near-physiological gait simulator with 6 DOF of the tibia. The kinematic data helps to clarify previous descriptions of several joint kinematics that are difficult to study in vivo. The methodology also provides a platform for researchers to explore more invasive foot biomechanics under dynamic and near-physiologic conditions.
Foot and ankle play an important role in bearing weight and managing normal daily activities such as walking. Given the fact that biomechanics is essential for foot and ankle surgery, it is a necessity for us, orthopedists, to improve our understanding of foot biomechanics, especially the 3D kinematics of the foot bones during human gait.
Previous studies conducted for the measurement of 3D kinematics of the foot bones during gait can be classified into four categories [1,2,3,4,5,6,7,8,9,10], including skin marker-based motion capture system, fluoroscopy study, in vivo measurement with bone pins, and measurement through cadaver gait simulators. However, some limitations are found in previous research methods. First, the skin marker-based motion capture system is less accurate due to skin marker artifacts [1, 11]. Next, fluoroscopy study was found less accurate in bony motion description due to the low frequency of observation and human-based inaccuracy in registration [8, 9, 12]. Some others placed bone pins to provide valuable foot kinematic data in vivo models, but the highly invasive procedure made it hard to apply in living subjects [4, 13].
Previous scholars [4,5,6, 14,15,16,17,18,19,20] have made many types of cadaver gait simulators for foot biomechanical research, and some have used these machines to measure foot skeletal motion during gait directly. However, previous machines have their limitations, like unphysiological ground reaction force (GRF) [6, 14, 17, 19], low simulation gait speed [6, 20], and simplified tibia kinematics [6, 14, 18]. Whittaker made a special gait simulator which achieved 6 DOF of tibia control by rotating the force plate while keeping the foot static. Although this machine achieved six-degree freedom tibia control, this simulation way differed from traditional walking pattern. Furthermore, the peak vertical GRF was about 0.8 bodyweight (BW) during the gait, which was less than the normal value .
We have developed a cadaveric gait simulator with 6 DOF of tibia in a physiological way that moved the cadaver foot while keeping the force plate static. Furthermore, the peak of vertical GRF achieved physiological value ranging from 1.1 BW and 1.3 BW [22, 23]. The aim of this work was to study the foot bony motion and joints move in a gait cycle with this foot and ankle gait simulator.
This study was approved by the Ethics Committee of our hospital. Six right-sided specimens were from three male and three female donators, whose death ages ranged from 45 to 69 years old with mixed death histories as shown in Table 1. None of the donators had evidence of prior foot surgery.
These specimens were fresh frozen cadavers, defrosted overnight for the experiment. And they were dissected to provide access to leg tendons, but all structures remained intact below the malleoli. The nine tendons in the foot were divided into four bundles, including anterior group comprised of tibialis anterior, extensor hallucis longus, and extensor digitorum longus. The posterior bundle was Achilles tendon. Medial group comprised of tibialis posterior, flexor hallucis longus, and flexor digitorum longus. The lateral group comprised of peroneus brevis and peroneus longus. The muscle tissue of above nine tendons was dissected, but the tendon tissue remained intact.
Artificial muscle forces were generated by four motors, which connected to each of four groups of tendons through custom-made clamps, as shown in Fig. 1. The tibia was sawed off 250 mm upper from pelma when ankle was in a neutral position, and a 1.27-cm diameter drill was used to hollow out the tibia intramedullary canal. We inserted the metal bar into the tibia intramedullary canal and drilled two screws to fix the bar and tibia together. Another end of the bar was fixed in the simulator machine’s tibia control place with screws. A band hoop was used to tie tibia and fibula together.
Each of tibia, calcaneus, cuboid, navicular, medial cuneiform, and 1-3-5 metatarsals were drilled with 1.6 mm Kirschner wires for supporting markers, as shown in Fig. 1. A combination of the real markers and virtual markers (generated from static standing trial) was used to define the x-axis (medial/lateral), y-axis (anterior/posterior), and z-axis (vertical axis) of the local coordinate frames for each bone. When the foot was weight-bearing and static, the local frames for each bone were fit with the axes of the global frame, and the virtual markers were generated. The midline of the foot was aligned parallel to the y (anterior/posterior) axis of the global coordinate system and the tibia vertical.
The simulation process started at heel contact and stopped when only toes were on the ground, and the metatarsal phalangeal joint was close to maximal dorsiflexion, as shown in Fig. 2. We collected movements of four joints (the cuneonavicular joint, the first tarsometatarsal joint, the calcaneocuboid joint, the fifth tarsometatarsal joint) relative motions of adjacent bones (including cuboid to navicular and calcaneus to tibia) during six gait simulation progress of each foot cadaver. The range of motions (ROMs) between adjacent bones in 3D was recorded, and the kinematic patterns were determined during the walking process.
The cadaveric gait simulator had another six actuating motors to control tibia to achieve a six-degree freedom parallel mechanism, including the direction of up-down, anterior-posterior, medial-lateral, and rotation in sagittal, in coronal, and in the transverse plane. For the control system, we used closed-loop control and iterative learning control to have a more real gait cycle. Lower limb cadaveric gait simulators aim to reproduce a normative vertical ground reaction force (vGRF) in vitro. During the cadaver gait simulations, the operator uses their expert knowledge of lower limb muscle and joint function to make educated trial-and-error guesses as to which muscle or kinematic input should be adjusted to achieve the desired vGRF. Closed-loop feedback control of the vGRF would likely improve the in vitro vGRF tracking accuracy and reduce the number of preliminary tuning simulations necessary to achieve vGRF tracking . An iterative learning mechanism was to optimize the desired position trajectory of the hydraulic cylinder which controls the movement of the tibia. After a gait simulation trial, the recorded cylinder loading force would be analyzed, and the desired position trajectory of the cylinder would be adjusted by an iterative learning mechanism. Then, the optimized trajectory would be saved and utilized in the next iteration. This process would be repeated until the tibia loading force could converge to the target curve through several iterations .
Ground reaction forces were synchronously collected using a force plate system of 1000 Hz (FP4060-15-2000Bertec, USA), and kinematic data was recorded (100 Hz) using seven cameras (Qualisys (Mocap Camera–Miqus M3) Motion Capture Systems, Sweden) which were positioned close to the specimens. For each tested bone, the local coordinate system was parallel to the global coordinate system (x-axis: medial/lateral, y-axis: anterior/posterior, z-axis: vertical axis) when the foot was load-bearing. Furthermore, the foot midline was parallel to the y-axis, and the tibia was vertical. Simultaneously, the initial virtual marker cluster was also defined. In the subsequent trials, the position and orientation of each local coordinate system can be calculated by combining the real marker cluster and the initial virtual marker cluster.
There was a total of 6 cadaver feet, and every cadaver had 6 times of gait simulation in this experiment. During 36 (6 × 6) times of gait simulations, the mean value and the standard deviation (SD) of the total ROMs of the relative motions of two pairs of bone blocks (cuboid to navicular and calcaneus to the tibia) and four joints (the medial cuneonavicular joint, the first tarsometatarsal joint, the calcaneocuboid joint, and the fifth tarsometatarsal joint) were calculated on the sagittal, coronal, and transverse planes respectively, as shown in Table 2.
In the sagittal plane, the medial cuneonavicular joint, the first tarsometatarsal joint, the calcaneocuboid joint, and the fifth tarsometatarsal joint all dorsiflexed from heel-strike followed by plantar flexion. The cuboid respect to navicular kept constant during heel-strike, then began to evert at gait phase 25%, and kept constant afterward. In the coronal plane, cuboid to calcaneus everted at heel-strike followed by inversion, remained relatively constant during foot flat, and inverted at push-off stage. The first tarsometatarsal joint, the medial cuneonavicular joint, and the fifth tarsometatarsal joint were relatively constant during heel-strike and midstance but everted at the push-off. The cuboid respecting to navicular everted at heel-strike, and followed by inversion during midstance, then everted at the push-off stage. In the transverse plane, the first tarsometatarsal joint began to adduct at heel-strike, kept constant in the midstance, but adducted at the push-off stage. The medial cuneonavicular joint, the fifth tarsometatarsal joint, and the calcaneocuboid joint began to adduct at heel-strike, kept constant during foot flat, but abducted at push-off. The cuboid respect to navicular kept constant during heel-strike and foot flat but abducted during push-off, as shown in Fig. 3.
It is of great importance to quantify the bony motion of the foot and ankle during gait. This knowledge is the foundation to understand its normal and pathologic function, to determine finite element boundary of the foot, and to give guidance of prosthetic joint design and replacement. Cadaver gait simulator was developed for foot bony kinematics measurement, including Whittaker et al.’s  bony motion measurement via a gait simulator, which was considered as the most advanced gait simulator. Although Baxter et al.  used their gait simulator constructed in the same way as Eric C. Whittaker et al.’s to assessed ankle, subtalar, and talonavicular kinematics, they did not study facet joints in the foot. Therefore, we compared our gait simulated results with Lundgren’s results and Whittaker’s results, as shown in Table 3.
In general, we found good agreement between our kinematic data and the data from in vivo studies by Lundgren, which is deemed as the gold standard. The total ROM reported here was within ± 1 SD (standard deviation) of the data reported by Lundgren et al.  for 13 out of 18 angles, while 15 of the 18 reported angles were within ± 2 SD of the data reported by Lundgren et al. .
Compared with the results of Whittaker et al. , our description of bone motion was similar to Whittaker’s work for the most part. However, we employed different simulation way from that of Whittaker et al.  and Baxter et al. . Our machine achieved the goal of six freedom motion control of tibia more physiologically, while Eric C. Whittaker and Baxter J.R employed an unphysiological way by motivating ground and keeping cadavers still to reproduce gait. Furthermore, the longitudinal axis of cadavers in their simulating process was parallel to the floor instead of vertical as a human being walking way. Furthermore, the peak of GRF of our machine was equal to 1.1 BW and 1.3 BW, while Whittaker et al.  measured the kinematics of foot bony only in 75% bodyweight and Baxter J.R in 25% bodyweight. From our simulated results, we found that the movement during gait between cuboid and navicular (sag, 4.7°; corn, 6.1°; trans, 6.9°) was close to the Lundgren’s results (sag, 7.2°; corn, 8.8°; trans, 8.9°) but not as large as Whittaker’s  results(sag, 18.7°; corn, 4.9°; trans, 20.1°). And we assumed the gravity of cadaver itself and the insufficient GRF in simulation might influence the correction of the results in Eric C. Whittaker’s cadaver gait simulator.
There are several meaningful results from our work. First, we found that several joints in foot cannot be regard as a rigid body during gait process, especially the movements in intertarsal joints and tarsometatarsal joints, including the medial cuneonavicular joint (sag, 7.4°; corn, 6.6°; trans, 5.4°), the first tarsometatarsal joint (sag, 4.0°; corn, 6.6°; trans, 4.8°), the fifth tarsometatarsal joint (sag, 9.1°; corn, 7.2°; trans, 6.8°), the calcaneocuboid joint (sag, 6.7°; corn, 9.9°; trans, 7.5°), and navicular to cuboid (sag, 4.7; corn, 6.1; trans, 6.9). These joints provided movement of 11.4° in the sagittal plane, 13.6° in coronel plane, and 10.2° in the transverse plane during gait. These movements confirmed that they could not be regarded as a rigid body during gait , instead was an important complementary portion of foot motion during gait.
Second, we found that the medial column had less ROM than that in the lateral during gait. During the simulated walking stance, the joints in the medial column of foot remained constant, like the medial cuneonavicular joint and the first tarsometatarsal joint. However, the lateral column joints had more movement, like calcaneocuboid joint and the fifth metatarsal joint. Compared the motion of the first tarsometatarsal joint (sag, 4.0°; corn, 6.6°; trans, 4.8°), the fifth tarsometatarsal joint had a greater total ROM (sag, 9.1°; corn, 7.2°; trans, 6.8°), which suggested that the medial column had less ROM than the lateral during gait to provide a firm support for the weight. These findings gave instructions to our clinical surgery that it is unwise to do joint fusion or joint movement restriction in the fifth tarsometatarsal joint. It was reported by Davitt and Morgan  that two flat-foot patients suffered the fifth metatarsal fatigue fractures after lateral column lengthening surgery. We thought that too much lengthening in the lateral necessarily might restrict joint motion, which caused stress concentration. Nevertheless, it was admitted to do fusion in medial joints, like treating severe Lisfranc injury.
Last, our results did not support the midtarsal locking mechanism proposed by Mann  that the relative midtarsal bones motion would cease to produce a rigid foot, which could effectively propel bodyweight during the later portion of walking stance. However, our in vitro kinematic results did not support the existence of the midtarsal locking mechanism during the stance phase. From Fig. 3, no rotational cease was observed in the calcaneocuboid joint and the medial cuneonavicular joint, as well as in the movement between navicular and cuboid during the later portion of stance as proposed by the precious pocking mechanism. Furthermore, we found greater rotation in the latter portion of the stance than early portion. Challenge to the traditional locking mechanism had also been reported by Okita et al.  and Chen et al. . They analyzed the midtarsal joint motion through a custom-made cadaveric gait simulator and fluoroscopic 3D–2D registration technique, respectively, and found the same phenomenon.
There are some limitations in the simulator, such as the reduced simulation velocity, faster increasing speed of the first peak of vertical GRF. And as other simulators machines, our cadaveric model did not simulate the intrinsic musculature force of the foot, which may be the reason why the results of bony motion were not correct enough compared to in vivo studies. Besides, some joints like subtalar joint and the first metatarsophalangeal joint were not included in the current study.
This study described foot bone motion using a biomechanically near-physical gait simulator of six-degree freedom of tibia. Our kinematic data generally agreed with invasive in vivo research and provided a realistic description of bony motion for an in vitro model. The kinematic data helped to clarify previous descriptions of the function of several joints that were difficult to study in vivo.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Degrees of freedom
- 3D kinematics of tibia:
Three dimensional kinematics of tibia
Ground reaction force
Range of motions
The first metatarsal
The fifth metatarsal
Leardini A, Benedetti MG, Berti L, Bettinelli D, Nativo R, Giannini S. Rear-foot, mid-foot and fore-foot motion during the stance phase of gait. Gait Posture. 2007;25(3):453–62.
de Asla RJ, Wan L, Rubash HE, Li G. Six DOF in vivo kinematics of the ankle joint complex: Application of a combined dual-orthogonal fluoroscopic and magnetic resonance imaging technique. J Orthop Res. 2006;24(5):1019–27.
Kaneda K, Harato K, Oki S, Ota T, Yamada Y, Yamada M, Matsumoto M, Nakamura M, Nagura T, Jinzaki M. Three-dimensional kinematic change of hindfoot during full weightbearing in standing: an analysis using upright computed tomography and 3D-3D surface registration. J Orthop Surg Res. 2019;14(1):355.
Lundgren P, Nester C, Liu A, Arndt A, Jones R, Stacoff A, Wolf P, Lundberg A. Invasive in vivo measurement of rear-, mid- and forefoot motion during walking. Gait Posture. 2008;28(1):93–100.
Peeters K, Natsakis T, Burg J, Spaepen P, Jonkers I, Dereymaeker G, Vander SJ. An in vitro approach to the evaluation of foot-ankle kinematics: performance evaluation of a custom-built gait simulator. Proc Inst Mech Eng H. 2013;227(9):955–67.
Nester CJ, Liu AM, Ward E, Howard D, Cocheba J, Derrick T, Patterson P. In vitro study of foot kinematics using a dynamic walking cadaver model. J Biomech. 2007;40(9):1927–37.
Baxter JR, Sturnick DR, Demetracopoulos CA, Ellis SJ, Deland JT. Cadaveric gait simulation reproduces foot and ankle kinematics from population-specific inputs. J Orthop Res. 2016;34(9):1663–8.
Kessler SE, Rainbow MJ, Lichtwark GA, Cresswell AG, D'Andrea SE, Konow N, Kelly LA. A Direct Comparison of Biplanar Videoradiography and Optical Motion Capture for Foot and Ankle Kinematics. Front Bioeng Biotechnol. 2019;7:199.
Phan CB, Shin G, Lee KM, Koo S. Skeletal kinematics of the midtarsal joint during walking: Midtarsal joint locking revisited. J Biomech. 2019;95:109287.
Zhang XB, Wu H, Zhang LG, Zhao JT, Zhang YZ. Calcaneal varus angle change in normal calcaneus: a three-dimensional finite element analysis. Med Biol Eng Comput. 2017;55(3):429–37.
Westblad P, Hashimoto T, Winson I, Lundberg A, Arndt A. Differences in ankle-joint complex motion during the stance phase of walking as measured by superficial and bone-anchored markers. Foot Ankle Int. 2002;23(9):856–63.
Chen WM, Geng X, Wang S, Xin MM, Xu WM, Jiazhang HM, Chao ZM, Li CM, Yang J, Wang K. In vivo kinematic study of the tarsal joints complex based on fluoroscopic 3D-2D registration technique. Gait Posture. 2016;49:54–60.
Arndt A, Westblad P, Winson I, Hashimoto T, Lundberg A. Ankle and subtalar kinematics measured with intracortical pins during the stance phase of walking. Foot Ankle Int. 2004;25(5):357–64.
Hurschler C, Emmerich J, Wulker N. In vitro simulation of stance phase gait part I: Model verification. Foot Ankle Int. 2003;24(8):614–22.
Sharkey NA, Hamel AJ. A dynamic cadaver model of the stance phase of gait: performance characteristics and kinetic validation. Clin Biomech (Bristol, Avon). 1998;13(6):420–33.
Sturnick DR, Demetracopoulos CA, Ellis SJ, Queen RM, Kolstov J, Deland JT, Baxter JR. Adjacent Joint Kinematics After Ankle Arthrodesis During Cadaveric Gait Simulation. Foot Ankle Int. 2017;38(11):1249–59.
Aubin PM WELW. A robotic cadaveric gait simulator with fuzzy logic vertical ground reaction force control. IEEE Trans Robotics 2012, 28(1): 246-255.
Kim KJ KHLZ. In vitro simulation of the stance phase in human gait. J Musculoskelet Res 2001, 5(2): 113-122.
Wulker N, Hurschler C, Emmerich J. In vitro simulation of stance phase gait part II: Simulated anterior tibial tendon dysfunction and potential compensation. Foot Ankle Int. 2003;24(8):623–9.
Aubin PM, Cowley MS, Ledoux WR. Gait simulation via a 6-DOF parallel robot with iterative learning control. IEEE Trans Biomed Eng. 2008;55(3):1237–40.
Whittaker EC, Aubin PM, Ledoux WR. Foot bone kinematics as measured in a cadaveric robotic gait simulator. Gait Posture. 2011;33(4):645–50.
Genrui Zhu ZWCY. Cadaver gait simulator imitative ability test and its biomechanics research application. Chin J Orthoped, 39(17): 1068-1074.
Qinyang Guo GSDW. Iterative learning based output feedback control for electro-hydraulic loading system of a gait simulator. Mechatronics 2018(54): 110-120.
Davitt JS, Morgan JM. Stress fracture of the fifth metatarsal after Evans’ calcaneal osteotomy: a report of two cases. Foot Ankle Int. 1998;19(10):710–2.
Mann RA. Biomechanics of the foot, AAOS Atlas of Orthotics: Biomechanical Principles and Application. St. Louis: C.V. Mosby Co; 1975.
Okita N, Meyers SA, Challis JH, Sharkey NA. Midtarsal joint locking: new perspectives on an old paradigm. J Orthop Res. 2014;32(1):110–5.
The authors would like to thank Anmin Liu (Centre for Rehabilitation and Human Performance Research, University of Salford, England, UK) and Dongmei Wang (School of Mechanical Engineering, Shanghai Jiaotong University, Shanghai 200240) for their supports. This work was funded by the National Natural Science Foundation of China (No. 81572176; No.8177090247).
Dr. Xu Wang is currently receiving grants from National Natural Science Foundation of China (Grant No. 81572176 and Grant No. 8177090247).
Ethics approval and consent to participate
This study was approved by the Institutional Review Board of Huashan Hospital, Fudan University.
Consent for publication
This manuscript is approved by all authors for publication in Journal of Orthopaedic Surgery and Research.
The authors declared that they have no conflicts of interest to this work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhu, G., Wang, Z., Yuan, C. et al. In vitro study of foot bone kinematics via a custom-made cadaveric gait simulator. J Orthop Surg Res 15, 346 (2020). https://doi.org/10.1186/s13018-020-01830-3
- Cadaveric gait simulator
- In vitro kinematics
- Foot bones