Skip to main content

Identification of diabetic neuropathic patients at risk of foot ulceration through finite element models and cluster analysis


Diabetic foot is an invalidating complication of diabetes mellitus that can lead to foot ulceration and amputations. While experimental analyses are limited solely to measurements of interfacial variables, three-dimensional (3D) patient specific finite element models (FEMs) of the foot can provide both the interfacial pressures and insight into internal stresses and strains tolerated by the plantar tissue [1]. FEMs allows quantifying the loads developed in the different anatomical structures of the foot and to understand how these affect foot tissue [2]. The aim of this study was to identify the neuropathic subjects at risk of ulceration with a cluster analysis classification of simulated plantar pressures and internal stresses. Simulations were ran with gait analysis data acquired 5 years prior to ulcerations.


Foot biomechanical analysis was carried out as in [3] on 16 diabetic neuropathic subjects by measn of a 6 cameras motion capture system (BTS, Padova), integrated and synchronized with 2 force plates (Bertec, USA), 2 plantar pressures systems (Imagortesi, Piacenza). For each patient the 3D kinematics, ground reaction forces and plantar pressures were calculated. Six of these subjects developed ulcers under metatarsals heads within 5 years after the acquisitions (ulcerated subjects (US)- age 62.3±4.1 years, BMI 26.3±2.0 kg/m2) while the other ten did not (non US - age 63.2±6.4 years, BMI 24.3±2.9 kg/m2).

In order to obtain the internal stresses (Von Mises and principal stresses) and the simulated plantar pressures (Figure 1), a recently developed 3D FEM [4] was adopted and the simulations were run adopting the experimental kinematic and kinetics as boundary conditions as in [4]. The midstance and the push-off phases of gait were considered as they are the instants when critical loads occur in the forefoot of the diabetic subjects.

Figure 1

Example of the results of a FEM simulation with ulcerated subject’s boundary conditions: Plantar pressures (left) and Von Mises internal stresses (right).

K-means and hierarchical cluster analysis were performed as in [5] with simulated plantar pressures and/or internal stresses as input.


The hierarchical method (Ward’s linkage) which led to the definition of three clusters (Table 1) gave the best result: 5 US were included in one cluster with only 3 non US.

Table 1 Results of the hierarchical cluster analysis: 3 clusters. Values are normalized over the subject’s weight. PP=plantar pressure.


A longer follow up is needed in order to verify whether the neuropathic subjects in cluster 2 and 3 will develop ulcers. A larger dataset is needed to further validate this methodology. Besides these limitations, results showed that combined FEMS and cluster analysis allowed to infer useful informations on the risk of ulceration even five years prior to the wound evolution.


  1. 1.

    Yarnitzky G, Yizhar Z, Gefen A: Real-time subject-specific monitoring of internal deformations and stresses in the soft tissues of the foot: a new approach in gait analysis. Journal of Biomechanics. 2006, 39: 2673-2689. 10.1016/j.jbiomech.2005.08.021.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Cavanagh P, Erdemir A, Petre M, Owings T, Botek G, Chokhandre S, Bafna R: Biomechanical factors in diabetic foot disease. Journal of Foot and Ankle Research. 2008, 1: K4-10.1186/1757-1146-1-S1-K4.

    PubMed Central  Article  Google Scholar 

  3. 3.

    Sawacha Z, Guarneri G, Cristoferi G, Guiotto A, Avogaro A, Cobelli C: Integrated kinematics–kinetics–plantar pressure data analysis: A useful tool for characterizing diabetic foot biomechanics. Gait & Posture. 2012, 36: 20-26. 10.1016/j.gaitpost.2011.12.007.

    Article  Google Scholar 

  4. 4.

    Guiotto A, Sawacha Z, Scarton A, Guarneri G, Avogaro A, Cobelli C: 3D finite element model of the diabetic neuropathic foot: a gait analysis driven approach. Proceedings ISB. 2013, Natal, Brasil

    Google Scholar 

  5. 5.

    Sawacha Z, Laudani L, Macaluso A, Vannozzi G: Identifying The Association Between Physical Activity Levels And Physiological Factors Underlying mobility: a descriptive data mining approach. Proceedings of IDAMAP. 2012, Pavia, Italy

    Google Scholar 

Download references

Author information



Corresponding author

Correspondence to Annamaria Guiotto.

Rights and permissions

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Cite this article

Guiotto, A., Sawacha, Z., Scarton, A. et al. Identification of diabetic neuropathic patients at risk of foot ulceration through finite element models and cluster analysis. J Foot Ankle Res 7, A27 (2014).

Download citation


  • Internal Stress
  • Ground Reaction Force
  • Plantar Pressure
  • Metatarsal Head
  • Neuropathic Patient