Identification des paramètres inertiels segmentaires humains

Identification of human body segment inertial parameters

• Auteur / Autrice : Marien Couvertier
• Direction : Patrick Lacouture, Tony Monnet
• Type : Thèse de doctorat
• Discipline(s) : Biomécanique et Bio-ingénierie
• Date : Soutenance le 18/12/2018
• Etablissement(s) : Poitiers
• Ecole(s) doctorale(s) : École doctorale Sciences et ingénierie des matériaux, mécanique, énergétique et aéronautique (Poitiers ; 2009-2018)
• Partenaire(s) de recherche : Laboratoire : Pôle poitevin de recherche pour l'ingénieur en mécanique, matériaux et énergétique - PPRIMME (Poitiers) - XLIM / XLIM
• Partenaire(s) de recherche : faculte : Université de Poitiers. UFR des sciences fondamentales et appliquées
• Jury : Président / Présidente : Hélène Pillet
• Jury : Examinateurs / Examinatrices : Patrick Lacouture, Tony Monnet, Thomas Robert
• Jury : Rapporteurs / Rapporteuses : Paul Fisette, Georges Dumont

Mots clés

Mots clés contrôlés

Filtres numériques (mathématiques)

Mécanique humaine

Identification des systèmes

Analyse numérique

Biomécanique

Mots clés libres

Biomécanique

Identification

Paramètres inertiels segmentaires

Méthodes numériques

Rigidification.

Biomechanics

Identification

Body segment inertial parameters

Numerical methods

Rigidification.

Résumé

L'objectif de ce travail est d'identifier les paramètres inertiels segmentaires humains, c’est-à-dire la masse, la position du centre de masse et la matrice d’inertie de ces segments. Ces paramètres, au nombre de dix par segment, constituent une donnée d'entrée indispensable aux calculs de dynamique inverse utilisés dans les études de biomécanique. Bien qu'il existe des méthodes pour avoir accès à ces paramètres par le biais de tables anthropométriques ou le calcul des volumes segmentaires, l'identification apparaît nécessaire dès lors que les sujets étudiés sont atypiques (handicapés physiques, femmes enceintes, sportifs présentant des hypertrophies musculaires). L'originalité de ce travail est de proposer une approche mixte dans l'écriture du problème d'identification combinant une formulation vectorielle et matricielle des équations du mouvement d’un système poly-articulé supposé rigide, qui ont déjà été établies au sein de l'équipe dans les travaux de thèse de Tony Monnet. La première permet d’identifier les masses et les centres de masse segmentaires. La deuxième permet d’identifier, elle, les matrices d’inertie segmentaires. Les paramètres d’entrée de cette méthode d’identification sont les matrices rotations segmentaires, leurs dérivées secondes, les accélérations segmentaires, ainsi que le torseur externe. Si ce dernier est directement mesuré par une plateforme de force, les autres grandeurs sont obtenues après des opérations sur la mesure de la cinématique segmentaire du sujet obtenue par un système opto-électronique. Ce système mesurant la cinématique du sujet grâce à des marqueurs cutanés, cette cinématique diffère de la cinématique théorique obtenue si les segments sont rigides, du fait des mouvements des masses molles. Ce travail a donc porté sur le calcul d’une matrice rotation optimale, basée sur une transformation matérielle décrivant le mouvement segmentaire.De plus les mouvements des masses molles ainsi que les instruments de mesure induisent un bruit dans les signaux cinématiques. Du fait de la double dérivation de ces signaux pour le calcul des accélérations segmentaires et des dérivées secondes des matrices rotations, ce bruit devient prépondérant sur le signal porteur. Ce travail a donc également porté sur le filtrage à effectuer pour atténuer ces bruits. Cinq filtres utilisés dans la littérature (filtre de Butterworth, lissage de Savitsky-Golay, moyenne glissante, lissage par spline et analyse spectrale) ont été implémentés et leurs effets sur les paramètres inertiels identifiés ont été comparés. Les résultats montrent que les paramètres identifiés avec la méthode vectorielle ne nécessitent pas de traitement. L’identification des matrices d’inertie nécessite, elle, un traitement et le lissage optimal est obtenu avec le moyenne glissante.Enfin, une modélisation du membre supérieur par une chaîne cinématique a également été implémentée afin de rigidifier la cinématique acquise. Les premiers résultats ne sont pas satisfaisants mais le modèle retenu peut être affiné avant de conclure sur l’intérêt de cette modélisation pour l’identification des paramètres inertiels. Finalement, l’approche mixte développée permet l'identification des dix paramètres inertiels des segments du corps humain. La méthode a été validée en identifiant les paramètres inertiels des segments constituant le membre supérieur de dix huit sujets. Les paramètres obtenus ont ensuite été comparés à ceux issus d’une table anthropométrique. Les résultats montrent que les paramètres identifiés sont très proches de ceux estimés. Cela montre donc que l’identification des paramètres inertiels est fiable et permet d’avoir accès aux paramètres inertiels de sujets atypiques, pour qui les tables anthropométriques ne sont pas disponibles.

The aim of this thesis is the identification of body segment inertial parameters (BSIP), i.e. the segment mass, center of mass location and inertia tensor. Those ten parameters per segment are a mandatory input for inverse kinetics methods which are widely used in biomechanics studies. Despite the fact that methods exist to estimate them from anthropometric tables or segment volumes measurements, identification is useful when subjects are atypical (such as disabled people, pregnant women or athletes with muscular hypertrophies). The originality of this work is to use a mixed approach to write the identification problem, combining a vectorial and a matrix formulations of rigid multi-body motion equations, based on previous work did in the RoBioSS axis by Tony Monnet during his PhD. The first one permit to identify segmental masses and center of mass locations. The second one identifies segmental inertia tensors.Inputs of identification algorithm are rotation matrices, their second derivatives, segmental accelerations, and external torsor. Even though this external torsor is directly measured with a force plate, the others inputs are derived from kinematics measurements performed by an optoelectronical device. This device measures kinematics with skin mounted markers tracked by cameras, and the obtained kinematics deviate from the theoretical kinematics of rigid bodies, because of the soft tissues artefacts. In order to deal with these artefacts an optimal rotation matrix computation, based on material transformation, has been performed.Also, noise appears during measurement because of the soft tissues artefacts and the measure device. When double numerical derivatives are applied, this noise becomes greater than the carrier signal. In order to deal with it, five filters, i.e. Butterworth filter, Savitsky-Golay smoothing, sliding average window, spline smoothing and singular spectrum analysis, taken from literature have been implemented and compared. Results show that BSIP identify from vectorial formulation didn’t need any filtering. On the other hand, inertia tensors identification needed smoothed inputs and the best way to smooth them was the sliding average window.Finally, a kinematic chain model of the upper limb has been implemented to rigidify the kinematics. Preliminary results aren’t satisfying but the chain model can be improved before assuming kinematic chain aren’t well suited to enhance BSIP identification. Ultimately, the developed mixed approach has been validated by upper limb inertial parameters identification of eighteen subjects. Identified inertial parameters have also been compared with ones estimated with an anthropometric table. The conclusion is that the identified parameters were very close to the estimated ones, which shows that identification will be reliable to estimate inertial parameters of atypical subjects for whom anthropometric tables aren’t available.