Bài s 2
Cho m t c c u robot 2DOF nh hình v ơ ư
S d ng matlab- Simulink
1. Xây d ng mô hình toán h c robot
2. Thi t k b đi u khi n PID bù tr ng tr ng. Ch y mô ph ngế ế ườ
S li u c a robot:
Chi u dài thanh n i 1: l 1 = 0,5 m
Chi u dài thanh n i 2: l 2 = 0,4 m
Kh i l ng thanh n i 1: m ượ 1 = 1 Kg
Kh i l ng thanh n i 2: m ượ 2 = 1 Kg
Momen quán tính kh p 1 quay quanh tâm kh i: J 1 = 0,5 Kgm2
Momen quán tính kh p 2 quay quanh tâm kh i: J 2 = 0,5 Kgm2
Kho ng cách t kh p 1 đ n tâm thanh n i 1: l ế g1 = 0,25 m
Kho ng cách t kh p 2 đ n tâm thanh n i 2: l ế g2 = 0,2 m
Góc ban đ u c a các thanh n i: θ 1 = 0 rad ; θ2 = Л/6 rad
Góc đ t c a các thanh n i: θ 1 = Л/3 rad ; θ2 = Л/2 rad
1. Xây d ng mô hình robot
Gi thi t kh p 1 sinh ra momen M ế 1 tác d ng gi a b thanh n i 1;
kh p 2 sinh ra momen M2 tác d ng gi a thanh n i 12; tr ng l c có h ng ướ
theo tr c y
Ký hi u đ ng năng và th năng c a thanh n i i là K ế i và Pi (i=1,2)
Ta có:
Đ ng năng thanh n i 1:
Th năng thanh n i 1:ế
T a đ tâm kh i thanh n i 2:
T c đ tâm kh i thanh n i 2:
→Bình ph ng t c đ dài tâm thanh n i 2:ươ
Đ ng năng thanh n i 2:
Th năng thanh n i 2:ế
Hàm lagrange c a robot:
Mômen c a kh p 1:
Mômen c a kh p 2:
Đ t ,
Ta có mô hình toán h c robot:
Trong đó:
: ma tr n quán tính
V i:
: thành ph n momen nh t và h ng tâm ướ
: thành ph n momen tr ng l c
V i:
Th c hi n mô ph ng trên matlab:
Mô hình trên matlab:
Trong đó:
- Hàm tính G:
function y = G(q)
% thong so ro bot
m1 = 1; m2 = 1;
l1 = 0.5; l2 = 0.4;
g = 9.8;
J1 = 0.5; J2 = 0.5;
lg1= l1/2; lg2 = l2/2;
% ma tran G
g1 = m1*g*lg1*cos(q(1))+m2*g*[l1*cos(q(1))+lg2*cos(q(1)+q(2))];
g2 = m2*g*lg2*cos(q(1)+q(2));
y = [g1;g2];
- Hàm tính V:
function y = V(q,dq)
% thong so robot
m1 = 1; m2 = 1;
l1 = 0.5; l2 = 0.4;
g = 9.8;
J1 = 0.5; J2 = 0.5;
lg1= l1/2; lg2 = l2/2;
% ma tran V
V1 = -m2*l1*lg2*sin(q(2))*[dq(2)^2+2*dq(1)*dq(2)];
V2 = m2*l1*lg2*sin(q(2))*dq(1)^2;
y = [V1;V2];
- Kh i
Hàm tính
function y = H(q,M)
% thong so robot
m1 = 1; m2 = 1;
l1 = 0.5; l2 = 0.4;
g = 9.8;
J1 = 0.5; J2 = 0.5;
lg1 = l1/2; lg2 = l2/2;
% bien khop
theta1=q(1);
theta2=q(2);
% ma tran H(q)
H11=m1*lg1^2+J1+m2*(l1^2+lg2^2+2*l1*lg2*cos(theta2))+J2;
H12=m2*(lg2^2+l1*lg2*cos(theta2))+J2;
H21=H12;
H22=m2*lg2^2+J2;
H = [H11 H12;H21 H22];
% ma tran nghich dao H^(-1)(q)
Hinv = inv(H);
% Gia toc khop
y = Hinv * M;
T o kh i robot (Create Subsystem)