Chapitre II Modélisation des bras de robots

II.1 Introduction
La conception et la commande des robots nécessitent le calcul de certains modèles
mathématiques, tels que :
– les modèles de transformation entre l'espace opérationnel (dans lequel est définie la
situation de l'organe terminal) et l'espace articulaire (dans lequel est définie la
configuration du robot). On distingue :
- les modèles géométriques direct et inverse qui expriment la situation de l'organe terminal
en fonction des variables articulaires du mécanisme et inversement.
- les modèles cinématiques direct et inverse qui expriment la vitesse de l'organe terminal
en fonction des vitesses articulaires et inversement.
– les modèles dynamiques définissant les équations du mouvement du robot, qui
permettent d'établir les relations entre les couples ou forces exercés par les actionneurs et
les positions, vitesses et accélérations des articulations [12].

II.2. Modélisation géométrique

II.2.1 La transformation homogène 
II.2.1.1 Système de coordonnées décrivant la position d’un solide dans l’espace
On choisit souvent le système de coordonnées qui convient le mieux à l’architecture du
robot et la nature de la tâche.
a) Système de coordonnées cartésiennes :
Il permet d’obtenir d’une manière directe les composantes du vecteur de position P.

Px
[]
P= Py
Pz

Figure II.1: Représentation cartésienne

Chapitre II Modélisation des bras de robots

b) Système de coordonnées cylindriques :

Cette représentation donne les composantes du vecteur de position tel que :
P = [r cos θ r sin θz]T
La relation qui lie les coordonnes cylindrique et cartésiennes et donnée par :

r=√ Px2 + Py 2

{ θ=arctg
z=Pz
Py
Px

c) Système de coordonnées sphériques :

Cette représentation donne le vecteur de position tel que :
P=[r cos θ sin ∅ r sinθ sin ∅ r cos ∅ ¿ ¿T
La relation qui lie les coordonnées sphériques et cartésiennes est tel que :

r =√ Px 2 + Py2 + Pz2

{ θ=arctg

∅=arctg √
Py
Px
Px2 + Py2
Pz

II.2.1.2 Matrice de transformation homogène 

Elle est utilisée pour décrire le lien géométrique entre les objets dans l’espace. Elle
permet d’associer à un point, une position et une orientation dans l’espace.
Une matrice de transformation nT m(4X4) permet le passage d’un repère Rm à un
repère Rn par rotation et translation :
S x N x A x Px Px
n

[S Ny Ay Py
T m= y
Sz N z A z P z
0 0 0 1
] []
avec P =
Py
Pz
1
Sx N x A x

[ ]
et Rot(r,θ)= S y N y A y
Sz N z Az

Rot(r,θ), est dite matrice d’orientation du référentiel Rm relative au référentiel Rn.

Le vecteur P est dit vecteur de position, les trois premières composantes
représentent les coordonnées cartésiennes de l’origine du référentiel Rm par rapport
référentiel Rn et la quatrième composante mise à 1 et dite facteur d’échelle.
Dans ce qui suit, on considère le repère fixe associé à la base du robot.
Chapitre II Modélisation des bras de robots

 Matrice de translation 
On définit une matrice de translation qui permet de représenter le repère R1 par rapport
au repère R0 comme suit :
1 0 0 Px

[ ]
0 10 P
Trans( P x , P y , P z) = 00 1 Py
0001
z

Avec P x , P y , P z les coordonnées de l’origine du repère R1 dans le repère R0 .

 Matrice de rotation 
Matrice de rotation autour de l’axe OX : on définit une matrice de rotation, d’un angle a
autour de l’axeO X 0, permettant de représenter le repère R1 par rapport au repère R0
comme suit :
1 0 0 0

[
0 cos α −sin α 0
Rot(x,α ) = 0 sin α cos α 0
¿000 1
]

Figure II.2: Représentation d’une rotation autour l’axeOx.

Chapitre II Modélisation des bras de robots

Matrice de rotation autour de l’axe OY : on définit la matrice de transformation de

rotation, d’un angle φ autour de l’axe OY0, permettant l’expression du repère R1 dans le
rpère R0 par :

cos φ 0 sin φ 0
Rot(y,φ ) =
[ 0 1 0 0
−sin φ 0 cos φ 0
000 1
]

Figure II.3: Représentation d’une rotation autour de l’axe Oy.

Matrice de rotation autour de l’axe OZ: on définit la matrice de transformation de

rotation, d’un angle ψautour de l’axe OY 0, permettant l’expression du repère R1 dans le
rpère R0 par :
cos ψ 0 −sin ψ 0
Rot(y,ψ ) =
[sin ψ
0
¿0
1 cos ψ 0
0 0 0
0 01
]
Chapitre II Modélisation des bras de robots

Figure II.4: Représentation d’une rotation.

II.2.2 Description géométrique 

La modélisation des robots de façon systématique et automatique exige une
méthode adéquate pour la description de leur morphologie. Plusieurs méthodes et notations
ont été proposées. La plus répandue est celle de Denavit-Hardenberg [12]. Dans les années
1950, les messieurs Jacques Denavit et Richard Hartenberg ont eu l'excellente idée de
proposer une méthode simple et systématique pour placer des référentiels sur chaque lien
d'un mécanisme sériel qui facilite énormément le calcul des matrices de transformation
[13].
II.2.3 Modèle géométrique direct (MGD)
Le modèle géométrique direct (MGD) est l'ensemble des relations qui permettent
d'exprimer la situation de l'organe terminal, c'est-à-dire les coordonnées opérationnelles du
robot, en fonction de ses coordonnées articulaires. Dans le cas d'une chaîne ouverte simple,
il peut être représenté par la matrice de transformation 0T n :
0 0 1 n−1
T n= T 1 ( q 1 ) T 2 ( q2 ) ⋯ T n ( qn ) (II.1)
Le modèle géométrique direct du robot peut aussi être représenté par la relation :
X =F( q) (II.2)
q étant le vecteur des variables articulaires tel que :
T
q=[ q 1 q 2 ⋯ qn ] (II.3)
Les coordonnées opérationnelles sont définies par :
X =[ X 1 X 2 ⋯ X m ]T (II.4)
Plusieurs possibilités existent pour la définition du vecteur X. Par exemple, avec les
éléments de la matrice0T n  :
Chapitre II Modélisation des bras de robots

X =[ n x n y n z o x o y o z a x a y a z p x p y pz ] (II.5)
Enfin, si nous avons un référentiel outil, et un référentiel atelier, la pose du premier par
rapport au dernier est définie par l'équation suivante:
atelier
T outil =atelierT outil 0T nl nT outil (II.6)
II.2.4 Modèle géométrique inverse(MGI)
On a vu que le modèle géométrique direct d'un robot permettait de calculer les
coordonnées opérationnelles donnant la situation de l'organe terminal en fonction des
coordonnées articulaires. Le problème inverse consiste à calculer les coordonnées
articulaires correspondant à une situation donnée de l'organe terminal. Lorsqu'elle existe,

la forme explicite qui donne toutes les solutions possibles (il y a rarement unicité de
solution) constitue ce que l'on appelle le modèle géométrique inverse. On peut distinguer
trois méthodes de calcul du MGI [12] :
 La méthode de Paul qui traite séparément chaque cas particulier et convient pour la
plupart des robots industriels.
 La méthode de Pieper qui permet de résoudre le problème pour les robots à six
degrés de liberté possédant trois articulations rotoïdes d'axes concourants ou trois
articulations prismatiques.
 La méthode générale de Raghavan et Roth, donnant la solution générale des robots
à six articulations à partir d'un polynôme de degré au plus égal à 16.
Quand le modèle géométrique inverse ne peut pas être obtenu, des techniques
numériques peuvent être employées. Ces techniques emploient la méthode de Newton-
Raphson ou des méthodes fondées sur la transpose de la matrice Jacobéenne.
Soit i-1TdEla matrice de transformation homogène représente la position désirée du
repère outil RE par rapport au repère atelier Rf. En général, on peut exprimer i-1TdE sous la
forme :
f
T dE=Z0T n (q) E (II.7)
Tels que :
 Z est la matrice de transformation définissant le repère de base du robot dans le
repère atelier R f .
0
 T n est la matrice de transformation du repère terminal dans le repère R0.
 E est la matrice de transformation du repère outil RE dans le repère terminal Rn.
Chapitre II Modélisation des bras de robots

Mettant tous les termes connus dans le membre gauche, on obtient :

U0 =0Tn(q) (II.8)
Avec : U0 = Z-1 f T dE E-1
Le MGD donnantX =F( q),q=¿ [q1 q2 ⋯ q n ]T X =[ X 1 X 2 ⋯ X m ]T , ou n est le nombre
de coordonnées opérationnelles et m le nombre de coordonnées articulaires, le problème
au-dessus (II.10) s'agit de résoudre un système de m équations a n inconnues, ce système
étant non linéaire.

Le nombre de solutions dépend de l'architecture du robot manipulateur et de

l'amplitude des articulations.
Trois cas se présentent pour calculer le MGI :
1. Solutions en nombre fini.
2. Absence de solution, lorsque la position de l'organe terminal désirée est en
dehors de la zone accessible du robot.
3. Infinité de solutions lorsque :
-Le robot est redondant vis-à-vis la tâche.
-Le robot se trouve dans certaines configurations singulières.
Lorsqu'il est possible de calculer toutes les configurations permettant d'atteindre
une situation donnée d'un robot manipulateur ce dernier est dit résoluble.
II.2.4.1 La méthode de Paul
Soit la matrice de transformation homogène d'un robot manipulateur (II.1) [12] :
0
T n=0T 1 ( q 1 ) 1T 2 ( q2 ) ⋯ n−1T n ( qn )
Soit U 0 la situation désirée telle que :
S x N x A x Px

[
S Ny Ay Py
U0 = y
Sz N z A z P z
0 0 0 1
On cherche à résoudre le système d'équations suivant :
] (II.9)

U 0 =0T 1 ( q1 ) 1T 2 ( q 2 ) ⋯ n−1T n ( q n ) (II.10)

Pour résoudre le système (II.10), Paul a proposé une méthode (Méthode de Paul)
qui consiste à pré multiplier successivement les deux membres de l'équation par les
Chapitre II Modélisation des bras de robots
i
matrices Ti-1 pour i allant de 1 à n-1.Ces opérations permettent d'isoler les variables
d'articulations l'une après l'autre.
U0=0T1 1T2 2T3…n-1Tn
1
T0 U0=0T1 1T2 2T3…n-1Tn
2
T1 U0=0T1 1T2 2T3…n-1Tn
3
T2 U0=0T1 1T2 2T3…n-1Tn
⋮
n-1
Tn-2 Un-2=n-1TnUj+1 = j+1Tn=j+1TnUj, pour j = 0,…,n.

Pour un robot à six degrés de liberté, on procède comme suit :

U0=0T1 1T2 2T3 3T44T5 5T6
1
T0U0=1T2 2T33T44T5 5T6
2
T11T0 U0=2T3 3T44T5 5T6
3
T2 2T11T0 U0 =3T4 4T55T6
4
T3 3T22T11T0 U0=4T55T6
5
T44T33T22T11T0 U0=5T6
II.2.5 Convention de Denvait-Hartenberg modifiée pour le calcul du MGD
Les paramètres de la convention Denavit-Hartenberg modifiée permettent de
disposer d'un paramétrage des liaisons tel que les matrices de passage aient toutes la même
forme littérale, ce qui facilite les calculs.
La méthode qui suit s'applique lorsque le robot correspond à une chaîne simple
ouverte et que ses articulations sont rotoïdes, ou prismatiques (ce qui est le cas en général).
Les corps constituant le robot sont supposés parfaitement rigides et connectés par des
articulations idéales (pas de jeu mécanique, pas d'élasticité) [6].

Figure II.5: Paramètres géométrique

Le passage du repère R j−1au repère R j se fait à travers 4 étapes :

Chapitre II Modélisation des bras de robots

- Une rotation d'angle α j autour de l'axe⃗

x j−1.
- Une translation de distance d jmesurée le long de l’axe ⃗
x j−1.
- Une rotation d'angleθ j autour del ' axe ⃗
Z j.
- Une translation de distancer jmesurée le long de l’axe⃗
Z j.

Avec :
 α j l'angle de rotation entre les axes ⃗
Z j−1 et ⃗
Z j autour de l’axe⃗
x j−1.

 d jla distance entre ⃗

Z j−1 et ⃗
Z j mesurée le longe de l’axe ⃗
x j−1.

 θ jl’angle de rotation entre les axes ⃗

x j−1 et ⃗
x j autour de l’axe ⃗
Z j.

 r jla distance entre ⃗

x j−1 et ⃗
x j mesurée le longe de l’axe ⃗
Zj .
La variable articulaireq j, associée à l'articulation j, est soit θ j, soit, selon que cette
articulation est respectivement de type rotoïde ou prismatique, ce qui se traduit par la
relation :
q j= (1-σ j)θ j+ σ j r j
La matrice de transformation homogène définissant le R j dans le repère R j−1 est
donnée par :
T j−1 , j= Rot( x j ,α j ).Trans ( x j−1,d j). Rot( z j ,θ j). Trans ( z j ,r j)

Cθ j −S θ j 0 dj
T j−1 , j=
[
Cα j S θ j C α j C θ j −Sα j −r j S α j
Sα j S θ j S α j C θ j Cα j −r j C α j
0 0 0 1
]
Chapitre II Modélisation des bras de robots

Figure II.6: Robot à structure ouverte simple.

II.3 Modélisation cinématique

II.3.1 Modèle cinématique direct
Le modèle cinématique direct d'un robot manipulateur décrit les vitesses des
coordonnées opérationnelles en fonction des vitesses articulaires. Il est noté [12] :
Ẋ =J (q) q̇ (II.11)
Où J(q) désigne la matrice jacobéenne de dimension (m× n) du mécanisme, égale à

∂X
et fonction de la configuration articulaireq. La même matrice jacobéenne intervient
∂q
dans le calcul du modèle différentiel direct qui donne les variations élémentaires dX des
coordonnées opérationnelles en fonction des variations élémentaires des coordonnées
articulairesdq, soit :
dX =J (q)dq
L'intérêt de la matrice jacobéenne est multiple :
– elle est à la base du modèle différentiel inverse, permettant de calculer une solution
locale des variables articulaires qconnaissant les coordonnées opérationnellesX .
– en statique, on utilise le jacobéen pour établir la relation liant les efforts exercés par
l'organe terminal sur l'environnement aux forces et couples des actionneurs.
– elle facilite le calcul des singularités et de la dimension de l'espace opérationnel
accessible du robot.
Chapitre II Modélisation des bras de robots

II.3.2 Modèle cinématique inverse

L'objectif du modèle cinématique inverse est de calculer, à partir d'une
configurationq donnée, les vitesses articulaires q̇ qui assurent au repère terminal une
vitesse opérationnelle Ẋ   imposée. Cette définition est analogue à celle du modèle
différentiel inverse. Ce dernier permet de déterminer la différentielle articulaire dq
correspondant à une différentielle des coordonnées opérationnellesdX spécifiée. Pour
obtenir le modèle cinématique inverse, on inverse le modèle cinématique direct en
résolvant un système d'équations linéaires. La mise en œuvre peut être faite de façon
analytique ou numérique :
1. la solution analytique a pour avantage de diminuer considérablement le nombre
d'opérations, mais on doit traiter séparément tous les cas singuliers.
2. les méthodes numériques sont plus générales, la plus répandue étant fondée sur
la notion de pseudo-inverse.
3. les algorithmes traitent de façon unifiée les cas réguliers, singuliers et
redondants. Ils nécessitent un temps de calcul relativement important.
Nous présentons dans ce paragraphe les techniques à mettre en œuvre pour établir
un modèle cinématique inverse dans les cas réguliers, singuliers et redondants.

II.3.3 Forme générale du modèle cinématique

Soit X =¿ [ X Tp X Tr ]T une représentation quelconque dans le repère R0 de la situation

du repère Rn fixé à un solide, les éléments X p et X r désignant respectivement la position et
l'orientation opérationnelles du solide. Lesrelations entre les vitesses Ẋ p et Ẋ r et les

vecteurs vitesses V 0n et ω 0n du repère Rn sont telles que :

Ẋ P Ω 0 3 V 0n V 0n
[ ] [ ][ ] [ ]
Ẋ r
=
03 Ω ω 0n
= Ω
ω0n
(II.12)

Les matrices Ω p et Ω r dépendant de la représentation choisie respectivement pour la

position et pour l'orientation.Le modèle cinématique direct a pour forme générale :
Ẋ =Jq̇
Dans le cas régulier la matrice jacobéenne J est carrée d'ordre n et son déterminant
est non nul. La méthode la plus générale consiste à calculer J-1, la matrice inverse de J,
qui permet de déterminer les vitesses articulaires q̇ grâce à la relation :
Chapitre II Modélisation des bras de robots

q̇ = J −1 Ẋ (II.13)

II.4 Modélisation dynamique

II.4.1 Modèle dynamique inverse :
Les équations dynamiques du mouvement d'un robot manipulateur sont l'ensemble
d'équations mathématiques décrivant le comportement dynamique des robots. De telles
équations sont utilisables pour la simulation du mouvement du robot manipulateur.
Le modèle dynamique est la relation entre les couples (et/ou forces) appliquée aux
actionneurs et les positions, vitesse et accélération articulaires. On représente donc le
modèle dynamique par une relation de la forme:
τ =f (q ,(q)̇ , q ̈ , f e )(II.14)
 τ : Vecteur des couples/ forces des actionneurs, selon que l'articulation est rotoïde
ou prismatique.
 q : Vecteur de position articulaire.
 q̇ : Vecteur de vitesse articulaire.
 q̈ : Vecteur d'accélération articulaire.
 fe: Vecteur représentant l’effort extérieur (force et moment) qu’exerce le robot sur
l’environnement.
II.4.2 Modèle dynamique direct
Le modèle dynamique direct exprime la position, la vitesse et l'accélération des
articulations en fonction du couple appliqué. Il est obtenu par inversion du modèle
précédent. Il s'écrit  [14] :
q̈ = (f (q, q̇,τ , f e) (II.15)
Le modèle dynamique des robots peut être obtenu par des lois mécaniques de
Newton et de Lagrange. Le formalisme d'E-L nous permet d'aboutir aux équations du
mouvement du robot.
II.4.3 L'approche de E-L
Le formalisme de E-L est utilisé pour modéliser le comportement dynamique d'un
robot. Cette approche particulière est assez simple à mettre en œuvre et elle est bien
adaptée de calcul manuel ainsi qu'aux méthodes de calcul assistées par ordinateur.
La dérivation des équations dynamiques d'un robot manipulateur à n degrés de
liberté est basée sur la compréhension de:
Chapitre II Modélisation des bras de robots
i-1
1. La matrice de transformation de coordonnées homogènes Ti qui décrit la
relation spatiale entre le i eme et le (i−1)eme système de coordonnées.Elle relie un
point fixe dans le corps i exprimé en coordonnées homogène par rapport au i eme
système de coordonnées au (i−1)eme système de coordonnées.
2. Formulation LE
Elle est basée sur l’équation du mouvement de Lagrange-Euler :
d ∂L ∂L
( ) − ¿ τ i=1,2 , ⋯ , n
dt ∂ q̇i ∂ qi i
(II.16)

Avec
n: Nombre de d.d.l du bras.
L = E-U : fonction de Lagrange.
E : énergie cinétique totale du bras.
U : énergie potentielle totale du bras.
τ i : force généralisée appliquée au système à l’articulation i.
q i: Variable articulaire généralisée.

II.4.4 Formalisme de Newton-Euler

Cette formulation est basée sur des équations récursives et décrit le mouvement du
bras en appliquant le principe de d’Alembert pour de chaque articulation. Ce dernier
considère les conditions d’équilibre statique dans des situations dynamiques en tenant
compte des forces externes appliquées et celles induites par réaction. Le lien entre les
forces, les inerties et les accélérations est exprimé par l’équation de Newton et celle
d’Euler [15].
 L’équation de Newton exprime la force appliquée au centre de masse d’un
corps et l’accélération induite :
F i=m. a⃗ i (II.17)
 L’équation d’Euler exprime le moment qui doit être appliqué au corps pour le
mettre en mouvement de rotation de vitesse ω iet d’accélération ω̇ i :
N i=I i ω̇i +ω i × I i ωi (II.18)
Les calculs s’effectuent en deux étapes :
 Propagation directe des vitesses et des accélérations de la base vers le poignet.
 Rétro propagation des forces et des couples exercés, du poignet vers la base.
Chapitre II Modélisation des bras de robots

III.4.5 Vitesses des articulations d’un bras manipulateur :

La formulation de Lagrange-Euler requière la détermination de l’énergie cinétique
du système physique, ce qui revient à déterminer les vitesses des différentes articulations.
Nous définissons les points suivants :
r ii : Point fixe qui repose sur le segment i de coordonnées homogénes définies par rapport au
repère Ri.
r 0i : définit le même point r ii de coordonnés homogènes définies par rapport au repère R0 .
Tel que :
xi

0 0 i
[]
r ii= y i =( x i , y i , z i , 1 )
zi
1
T

Ce qui nous permet d’écrire : r i =T i r i

La vitesse du point r ii est donnée par :

d
i
∂ T 0i i
v =v i = r 0i =
0
i
dt ( j=0 j
i
)
∑ ∂ q q̇ j r i=∑ U ij q̇ j r ii
j=1

Par définition, on a :

∂ T 0i T 0j−1 Q j T ij−1 pour j≤ i.
U ij =
∂qj
=
{ 0 pour j>1.

1 −1 0 0
Avec :
[ ]
Q j= 0 0 0 0
0 0 00
0 0 00
i=1 , 2, … , n .

II.5 Génération de trajectoires

La génération des trajectoires dans l'espace cartésien, permet au robot d'aller d'une
position actuelle Pa vers une position finale Pf appartenant à l'espace de travail. C’est
connu que l'espace de travail est bien limité. Cette limitation est due principalement aux
contraintes mécaniques du robot, telles que les marges de variation de chaque articulation
et les langueurs des segments [16].
II.5.1 Méthode des polynômes cubique :
Chapitre II Modélisation des bras de robots

Considérez le problème de déplacer l'outil de sa position initiale à une position

finale et ce durant un certain temps tf. La cinématique inverse permet de calculer
l’ensemble des angles articulaires qui correspondent à une posture désirée du poignet. La
posture initiale du manipulateur est également connue sous la forme d’un ensemble
d’angles des articulations. Ce qui est requis est une fonction pour chaque articulation dont
la valeur à t0est la position initiale de l’articulation et dont la valeur à tf est la position finale
souhaitée, il existe de nombreuses fonctions lisses,θ(t )qui peuvent être utilisées pour
interpoler la valeur de l’articulation.
En faisant un seul mouvement lisse, au moins quatre contraintes surθ(t ) sont
évidentes. Deux contraintes sur la valeur de la fonction proviennent de la valeur initiale et
la valeur finale [5].
θ ( 0 ) =θ0
θ(t f ) ¿ θf (II.19)

Figure II.7: Plusieurs formes de chemin possible pour une seule articulation.

Deux autres contraintes sont que la fonction soit continue en vitesse, ce qui dans ce
cas signifie que la vitesse initiale et finale sont nulles :
θ ( 0 ) =0 , θ ( t f ) =0 (II.20)
Ces quatre contraintes peuvent être satisfaites par un polynôme d’au moins un
troisième degré (polynôme cubique à quatre paramètres) de la forme :
θ ( t )=a0 + a1 t 1+ a2 t 2 +a 3 t 3 (II.21)
De sorte que la vitesse et l’accélération de l’articulation sur ce chemin sont définies
θ̇ ( t ) =a 1+2 a2 t 1+3 a 3 t 2
θ̈ ( t ) =2 a2 +6 a3 t (II.22)
Chapitre II Modélisation des bras de robots

La combinaison (II.22) et (II.21) avec les contraintes de quatre désirs donne quatre
équations en quatre inconnues
θ ( 0 ) =0, 0=a 1,
θ f =a0 +a1 t 1f +a 2 t 2f + a3 t 3f
0=a1 +2 a2 t 1f +3 a3 t 2f (II.23)
Résoudre les équations de a i, pour l’obtention
a 0=θ ( 0 ), a 1=0
3 −3
a 2= 2 (θf −θ0 ), a 3= 3 (θ f −θ0 ) (II.24)
tf tf
En utilisant (II.24), on peut calculer le polynôme cubique qui relie toute position
d’angle d’articulation initiale avec une position finale souhaitée. Cette solution est pour les
cas où l’articulation commence et se termine à une vitesse nulle.

II.6 Conclusion :
Dans ce chapitre nous avons étudié les concepts de base de la robotique et la
modélisation du robot à "n" degrés de liberté, en basant sur les paramètres de D-H. Cette
formulation nous a permis d'exploiter le formalisme d'Euler- Lagrange, afin d'établir le
modèle géométrique, cinématique et dynamique que nous utilisons pour le calcul des
commandes. Nous avons développé le calcul du modèle dynamique pour un robot à deux
degrés liberté.