<figure><img src="image_1.png" /></figure>

HAL
open science

# Ground Reaction Forces and Moments Prediction of Challenging Motions: Fencing Lunges

Charles Pontonnier, Claire Livet, Antoine Muller, Anthony Sorel, Georges Dumont, Nicolas Bideau

► To cite this version:

Charles Pontonnier, Claire Livet, Antoine Muller, Anthony Sorel, Georges Dumont, et al.. Ground Reaction Forces and Moments Prediction of Challenging Motions: Fencing Lunges. Computer Methods in Biomechanics and Biomedical Engineering, 2019. hal-02142288

HAL Id: hal-02142288

https://hal.science/hal-02142288

Submitted on 28 May 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
# Ground Reaction Forces and Moments Prediction of Challenging Motions: Fencing Lunges

C. Pontonnier<sup>a1</sup>, C. Livet<sup>a</sup>, A. Muller<sup>b</sup>, A. Sorel<sup>c</sup>, G. Dumont<sup>a</sup> and N. Bideau<sup>c</sup>

<sup>a</sup>Univ Rennes, Inria, CNRS, IRISA, France; <sup>b</sup>Institut de recherche Robert-Sauvé en santé et en sécurité du travail (IRSST), Montréal, QC, Canada; <sup>c</sup>Univ Rennes, Inria, CNRS, M2S, France

**Keywords:** dynamics; optimization; GRF; motion; contact

## 1. Introduction

Classically, inverse dynamics studies of motion are based on a bottom-up approach measuring the Ground Reaction Forces and Moments (GRF&M) with force platforms. The field of study is constrained to the contact zones of these platforms, diminishing the ecological properties of the motion to be studied and constraining drastically the analysis. Such a limitation is critical for sports motion analysis, since motions to be studied can be quite large in amplitude and varied in their contacts with the ground.

GRF&M prediction offers the opportunity to enlarge the field of study without additional measure to be made. Optimization-based GRF&M prediction methods are using a set of contact points defined under the feet to generate forces acting on the body. Motion analysis is used to define at any moment the points in contact with the ground and an optimization routine finds the minimal set of forces to be applied on these points satisfying the dynamical equilibrium of the whole body. If such methods have been validated on daily living motions (Fluit et al. 2014) and some classical sports motions (Skals et al. 2017) such as running, vertical jump or side-cut motions, their capacity to predict GRF&M on complex sport motions is still to be proved.

Lunge motion is a fundamental attack of modern fencing, asking for a high level of coordination, speed and accuracy. It consists in an explosive extension of the front leg accompanying an extension of the sword arm. In such motions, the direction of action and the way feet are oriented – guard position - are particularly challenging for a GRF&M prediction method.

The current study aims at evaluating an optimization-based GRF&M prediction method on a large set of fencing lunge data.

## 2. Methods

### 2.1 Experimental data

The experimental data used in the current study is extracted from another one dealing with uncertainty during lunge motions (Sorel et al. 2019).

10 of 11 male regional to national level fencers (5 right handed, 5 left handed, aged 22±3 y.o., height

1.78±0.06m, weight 75±8kg) participated to this study, giving their informed consent after a presentation of the protocol (validated by the INRIA national ethics comitee). They performed a total of 120 fencing lunges each, under various conditions. 1 fencer of the original cohort was excluded from the present study due to unexploitable force platform data.

Whole body motion of the fencers was recorded thanks to a VICON optoelectronic system (24 cameras, 250 Hz). A set of 49 markers derived from the ISB recommendations was used. Each foot was placed under an AMTI Force platform (1000 Hz), enabling the measure of GRF&M under each foot independently during the whole motion. Lunges start-end were extracted with the method proposed in (Sorel et al. 2019). A whole body biomechanical model developed in the CusToM library (Muller et al. 2019) was used to perform the inverse kinematics analysis and the GRF&M prediction. Geometry of the model was scaled to the subjects (Muller et al. 2016) and masses and inertias were obtained from (Dumas et al. 2006).

### 2.2 GRF&M prediction and analysis

For the total of 1200 lunges extracted from the experimentation, a classical GRF&M prediction method was applied. A set of $N_f = 14$ contact points per foot was defined. At each frame, a Sequential Quadratic Programming (SQP) method found the solution of the following minimization problem:

$$ \min_F \sum_{i=1}^{2N_f} \|F_i\|^2 \\ \text{s.t.} \left\{ \begin{array}{l} M(q)\ddot{q} + C(q, \dot{q}) + G(q) + \lambda + E = 0 \\ \forall i \in [1, 2N_f], F_i < F_{i \max} \end{array} \right. $$

Where, by isolating the whole-body biomechanical model, $M(q)$ is the inertia matrix, $C(q, \dot{q})$ is the centrifugal and Coriolis force vector, $G(q)$ is the gravity force vector, $\lambda$ is the generalized internal force vector, $E$ is the generalized external force vector (computed from the contact force at each contact point $F_i$) and $F_{i \max}$ the vector containing the maximal force available for the contact point $i$. At each frame, a contact point was considered as active if it was close to

<footnote n="1">¹ Corresponding author: charles.pontonnier@ens-rennes.fr</footnote>
the floor and almost without motion. The distance and velocity thresholds were respectively 0.02 m and 0.8 m/s. When a contact point was respecting the thresholds, the associated force was limited to 0.4 Body Weight and had to respect the Coulomb's law of friction. A friction coefficient of 0.5 was used here (Skals et al., 2017). Predicted forces and moments of each foot were compared to the measured ones in terms of RMSE, RMSE normalized to body mass, rRMSE and Pearson's correlation coefficient.

### 3. Results and discussion

<table><caption>Table 1 – mean±std RMSE, rRMSE, RMSE/BW and r for all the trials and all the subjects.</caption><thead><tr><th rowspan="2"></th><th colspan="3">Left foot / Right foot</th><th rowspan="2">r</th></tr><tr><th>RMSE (N or N.m)</th><th>rRMSE (%)</th><th>RMSE/BM (N/kg or N.m/kg)</th></tr></thead><tbody><tr><td>F<sub>x</sub></td><td>12.1±3.8<br>11.2±3.6</td><td>21.8±10.8<br>20.0±7.7</td><td>0.16±0.05<br>0.15±0.04</td><td>0.64±0.24<br>0.69±0.23</td></tr><tr><td>F<sub>y</sub></td><td>71.9±11.9<br>71.0±11.2</td><td>17.2±3.4<br>15.5±2.6</td><td>0.97±0.20<br>0.96±0.18</td><td>0.89±0.04<br>0.90±0.04</td></tr><tr><td>F<sub>z</sub></td><td>66.5±68.6<br>65.8±68.5</td><td>7.0±6.0<br>6.8±4.8</td><td>0.85±0.77<br>0.84±0.77</td><td>0.96±0.07<br>0.96±0.08</td></tr><tr><td>M<sub>x</sub></td><td>89.1±78.5<br>81.1±76.5</td><td>5.7±2.8<br>7.4±5.9</td><td>1.14±0.88<br>1.05±0.87</td><td>0.98±0.03<br>0.94±0.14</td></tr><tr><td>M<sub>y</sub></td><td>15.2±7.7<br>17.8±9.5</td><td>7.3±2.7<br>6.4±2.1</td><td>0.20±0.08<br>0.23±0.12</td><td>0.97±0.03<br>0.98±0.02</td></tr><tr><td>M<sub>z</sub></td><td>31.3±12.0<br>29.7±11.1</td><td>32.2±13.5<br>26.3±20.7</td><td>0.41±0.13<br>0.40±0.17</td><td>0.50±0.26<br>0.68±0.30</td></tr></tbody></table>

Results presented Table 1 are very similar to some of the results presented in (Skals et al. 2017). The vertical forces $F_z$ are accurately estimated with a particularly low rRMSE. Most of correlation coefficients are high, excepted for the medio-lateral force $F_x$ and the vertical moment $M_z$. $F_x$ is particularly low during the whole motion - see RMSE/BM, explaining this issue, whereas $M_z$ is poorly predicted due to the bad estimation of the antero-posterior component $F_y$. This issue is particularly visible in Figure 1: at the beginning of the lunge, a static phase in guard position presents an important $F_y$ value of opposite value at each foot. Such a stabilization scheme, ensuring a good balance of the fencer, is not well predicted by the algorithm. Indeed, the minimization tends to find the minimal set of forces ensuring equilibrium, and therefore does not predict forces arising from specific motor control strategies, generally compensated in the whole-body dynamic equilibrium. This issue is a future way of improvement, by adding motor control features in the cost function/the constraints of the optimization scheme, such as balance or stiffness control.

In conclusion, such a method is already quite relevant to study specific and challenging motions in a less constrained field, even if improvements can be made in the definition of the optimization problem to be solved.

<figure><img src="image_2.png" /><figcaption>Figure 1 – Predicted vs. measured GRF : sample trial</figcaption></figure>

### References

* Fluit R, Andersen MS, Kolk S, Verdonschot N, Koopman HF. 2014. Prediction of ground reaction forces and moments during various activities of daily living. J Biomech

* Skals S, Jung MK, Damsgaard M., & Andersen MS 2017. Prediction of ground reaction forces and moments during sports-related movements. Multibody Syst Dyn

* Sorel A, Plantard P, Bideau N, Pontonnier C. 2019.Studying fencing lunge accuracy and response time in uncertain conditions with an innovative simulator, Submitted

* Muller A, Pontonnier C, Puchaud P, Dumont G. 2019. CusToM: a Matlab toolbox for musculoskeletal simulation, JOSS

* Muller A, Germain C, Pontonnier C, Dumont G. 2016. A simple method to calibrate kinematical invariants: application to overhead throwing. ISBS conference

* Dumas R, Cheze L, Verriest JP. 2007. Adjustments to McConville et al. and Young et al. body segment inertial parameters. J Biomech