Abstract: To assess whether the transient stresses of foot impact with the ground are similar to those found during midstance loading and if the location of high stress correlate with the sites most commonly associated with mechanically induced osteoarthritis (OA). We compared impact stresses in subchondral bone between two subject-specific, three-dimensional, finite-element models of the equine metacarpophalangeal (MCP) joint-one with advanced OA and one healthy, and with similar published data on the stresses that occur at midstance. Methods: Two right MCP joints (third metacarpal and proximal phalanx) were scanned using micro-computed tomography (μCT). Images were segmented, and meshed using modified 10-node quadratic tetrahedral elements. Bone material properties were assigned based on the bone density. An impact velocity of 3.55 m/s was applied to each model and contact pressures and stress distribution were calculated for each. In a separate iteration, the third metacarpal was loaded statically. A sampling grid of 160 equidistant points was superimposed over selected slices, and average peak stresses were calculated for 6 anatomical regions. Within-region maximal peak and average von Mises stresses were compared between healthy and OA bones in both midstance and impact loading. Results: Average impact stresses across all regions, in both locations (palmar and dorsal) were greater in the OA model. Highest impact stresses were located in the dorsal medial condyle in the healthy (12.8 MPa) and OA (14.1MPa) models, and were lowest in the palmar medial and lateral parasagittal grooves in the healthy (5.94 MPa) and OA (7.07 MPa) models. The healthy static model had higher peak (up to 49.7% greater) and average (up to 38.6% greater) stresses in both locations and across all regions compared to the OA static model. Conclusions: Under simulated footfall a trot, loading on the dorsal aspect of the third metacarpal at impact created stresses similar to those found during midstance. The high accelerations that occur under impact loading are likely responsible for creating the high stresses, as opposed to midstance loading where the high stresses are the result of high mass loading. Although the stress magnitudes were found to be similar among the two loading conditions, the location of the high stress loading occurred in sites that are not typically associated with osteoarthritic changes.
The Equine Research Bank provides access to a large database of publicly available scientific literature. Inclusion in the Research Bank does not imply endorsement of study methods or findings by Mad Barn.
This research summary has been generated with artificial intelligence and may contain errors and omissions. Refer to the original study to confirm details provided. Submit correction.
The research focuses on investigating whether the temporary stresses experienced when a horse’s foot hits the ground are similar to those during midstance loading and if these high-stress locations relate to the common sites for developing osteoarthritis (OA). The study used computer-based modeling to compare two conditions: one with advanced OA and one healthy.
Methods Utilized
The researchers used two right Metacarpophalangeal (MCP) joints from equine specimens for the study. MCP joints used were scanned using micro-computed tomography (μCT), which generated detailed images of the structure.
These images were then segmented and meshed using modified 10-node quadratic tetrahedral elements. Here “meshing” relates to the process of subdividing a larger, complex area (like MCP joint) into smaller, simpler parts for easier analysis.
Material properties of the bone were assigned based on the bone density, which allow accurate simulations in the model.
Each model was subjected to an impact velocity of 3.55 m/s to simulate a horse’s foot hitting the ground, with contact pressures and stress distribution calculated for each model.
In a separate iteration, the third metacarpal was loaded statically, that is, without movement, to simulate midstance loading.
The stress levels were sampled at 160 equidistant points selected over the area of the bone and average peak stresses were calculated for 6 anatomical regions.
The highest peak stresses and average von Mises stresses were then compared between healthy and OA bones under midstance and impact loading.
Results Gathered
Results revealed that average impact stresses across all regions, in both locations (palmar and dorsal) were higher in the OA model.
The highest stress areas during impact were located in the dorsal medial condyle in both the healthy and OA models, and were lowest in the palmar medial and lateral parasagittal grooves.
The healthy static model (representing a horse during midstance) showed higher peak and average stresses across all regions compared to the OA static model.
Conclusion
Under simulated footfall at a trot, loading on the dorsal aspect of the third metacarpal bone at impact created stresses similar to those found during midstance.
The high accelerations that occur under impact loading could be responsible for creating high stresses, as opposed to midstance loading where the high stresses result from high mass loading.
Even though the stress magnitudes were similar among the two loading conditions, the locations of high stress loading occurred in sites that are not typically associated with osteoarthritic changes. This offers new insights into understanding the development of OA conditions in equine models.
Cite This Article
APA
McCarty CA, Thomason JJ, Gordon KD, Burkhart TA, Milner JS, Holdsworth DW.
(2016).
Finite-Element Analysis of Bone Stresses on Primary Impact in a Large-Animal Model: The Distal End of the Equine Third Metacarpal.
PLoS One, 11(7), e0159541.
https://doi.org/10.1371/journal.pone.0159541
Radin EL, Martin RB, Burr DB, Caterson B, Boyd RD, Goodwin C. Effects of mechanical loading on the tissues of the rabbit knee.. J Orthop Res 1984;2(3):221-34.
Gustas P, Johnston C, Drevemo S. Ground reaction force and hoof deceleration patterns on two different surfaces at the trot. Equine and Comparative Exercise Physiology 3(4):209–16.
McClinchey HL, Thomason JJ, Jofriet JC. Isolating the effects of equine hoof shape measurements on capsule strain with finite element analysis. Veterinary and Comparative Orthopaedics And Traumatology 16: 67–75.
Oloyede A, Flachsmann R, Broom ND. The dramatic influence of loading velocity on the compressive response of articular cartilage.. Connect Tissue Res 1992;27(4):211-24.
Shepherd DE, Seedhom BB. The 'instantaneous' compressive modulus of human articular cartilage in joints of the lower limb.. Rheumatology (Oxford) 1999 Feb;38(2):124-32.
Li LP, Buschmann MD, Shirazi-Adl A. Strain-rate dependent stiffness of articular cartilage in unconfined compression.. J Biomech Eng 2003 Apr;125(2):161-8.
Langelier E, Buschmann MD. Increasing strain and strain rate strengthen transient stiffness but weaken the response to subsequent compression for articular cartilage in unconfined compression.. J Biomech 2003 Jun;36(6):853-9.
Austman RL, Milner JS, Holdsworth DW, Dunning CE. The effect of the density-modulus relationship selected to apply material properties in a finite element model of long bone.. J Biomech 2008 Nov 14;41(15):3171-6.
McErlain DD, Milner JS, Ivanov TG, Jencikova-Celerin L, Pollmann SI, Holdsworth DW. Subchondral cysts create increased intra-osseous stress in early knee OA: A finite element analysis using simulated lesions.. Bone 2011 Mar 1;48(3):639-46.
Burr DB, Radin EL. Microfractures and microcracks in subchondral bone: are they relevant to osteoarthrosis?. Rheum Dis Clin North Am 2003 Nov;29(4):675-85.
Colahan P, Turner TA, Poulos P, Piotrowski G. Mechanical functions and sources of injury in the fetlock and carpus. Proceedings of the annual convention of the American Association of Equine Practitioners 689–699.
Lopas LA, Shen H, Zhang N, Jang Y, Tawfik VL, Goodman SB, Natoli RM. Clinical Assessments of Fracture Healing and Basic Science Correlates: Is There Room for Convergence?. Curr Osteoporos Rep 2023 Apr;21(2):216-227.
Pramudita JA, Hiroki W, Yoda T, Tanabe Y. Variations in Strain Distribution at Distal Radius under Different Loading Conditions. Life (Basel) 2022 May 16;12(5).
Sakashita M, Yamasaki S, Yaji K, Kawamoto A, Kondo S. Three-dimensional topology optimization model to simulate the external shapes of bone. PLoS Comput Biol 2021 Jun;17(6):e1009043.
Paul D, Sta Maria D, Sadat SA, Rahman MA, Shahzad H, Apu EH, Shaikh MH. Exploring the Intricacies of Finite Element Modeling of 3D-Printed Scaffolds for Musculoskeletal Applications: An In-Depth Review. Cureus 2026 Jan;18(1):e101070.
Caffrey JM, Thomas PK, Appt SE, Burkart HB, Weaver CM, Kleinberger M, Gayzik FS. Contrast enhanced computed tomography of small ruminants: Caprine and ovine. PLoS One 2023;18(12):e0287529.