kom op

Dependencies:   QEI mbed

Committer:
JornJan
Date:
Thu Nov 02 09:45:25 2017 +0000
Revision:
2:0686b2132556
Parent:
1:7969189824ff
Child:
3:3b5b85a32c9a
eerste kinematica implementatie;

Who changed what in which revision?

UserRevisionLine numberNew contents of line
JornJan 0:921402a95675 1 #include "mbed.h"
JornJan 0:921402a95675 2 #include "Serial.h"
JornJan 0:921402a95675 3 #include "math.h"
JornJan 0:921402a95675 4 #include "QEI.h"
JornJan 0:921402a95675 5
JornJan 0:921402a95675 6 Serial pc(USBTX, USBRX); //Serial PC connectie
JornJan 0:921402a95675 7 DigitalOut motor1DirectionPin(D4); //Motorrichting op D4 (connected op het bord)
JornJan 0:921402a95675 8 PwmOut motor1MagnitudePin(D5); //Motorkracht op D5 (connected op het bord)
JornJan 0:921402a95675 9 DigitalOut motor2DirectionPin(D7); //Motorrichting op D4 (connected op het bord)
JornJan 0:921402a95675 10 PwmOut motor2MagnitudePin(D6); //Motorkracht op D5 (connected op het bord)
JornJan 0:921402a95675 11 QEI q1_enc(D13, D12, NC, 32); //encoder motor 1 instellen
JornJan 0:921402a95675 12 QEI q2_enc(D11, D10, NC, 32); // encoder motor 2 instellen
JornJan 0:921402a95675 13 const double pi = 3.1415926535897; // waarde voor pi aanmaken
JornJan 0:921402a95675 14
JornJan 0:921402a95675 15 Timer t;
JornJan 0:921402a95675 16 Ticker aansturing;
JornJan 0:921402a95675 17 Ticker stepres;
JornJan 0:921402a95675 18
JornJan 0:921402a95675 19 double ref1 = 0;
JornJan 0:921402a95675 20 double refrad1;
JornJan 0:921402a95675 21 double refrad2;
JornJan 0:921402a95675 22 double ref2 = 0;
JornJan 0:921402a95675 23 double Kp1 = 0.002; //kp motor 2 = 0.0001
JornJan 0:921402a95675 24 double Ki1 = 0;
JornJan 0:921402a95675 25 double Kp2 = 0.0001;
JornJan 0:921402a95675 26 double Ki2 = 0;
JornJan 0:921402a95675 27 int q1_puls;
JornJan 0:921402a95675 28 int q2_puls;
JornJan 0:921402a95675 29 double rad2pulses=4200/(2*pi);
JornJan 2:0686b2132556 30 double pulses2rad = (2*pi)/4200;
JornJan 0:921402a95675 31 double ts = 0.001;
JornJan 0:921402a95675 32 double PI1;
JornJan 0:921402a95675 33 int n;
JornJan 0:921402a95675 34
JornJan 0:921402a95675 35 double error1_1;
JornJan 0:921402a95675 36 double error2_1 = 0;
JornJan 0:921402a95675 37 double error_I_1;
JornJan 0:921402a95675 38 double error_I2_1 = 0;
JornJan 0:921402a95675 39 double error1_2;
JornJan 0:921402a95675 40 double error2_2 = 0;
JornJan 0:921402a95675 41 double error_I_2;
JornJan 0:921402a95675 42 double error_I2_2 = 0;
JornJan 0:921402a95675 43 double PI2;
JornJan 0:921402a95675 44
JornJan 2:0686b2132556 45 // kinematica gegevens
JornJan 2:0686b2132556 46 // lengte armen
JornJan 2:0686b2132556 47 double L1 = 0.250;
JornJan 2:0686b2132556 48 double L2 = 0.355;
JornJan 2:0686b2132556 49 double L3 = 0.150;
JornJan 2:0686b2132556 50
JornJan 2:0686b2132556 51 // reference position
JornJan 2:0686b2132556 52 double q1=0; // positie q1 in radialen
JornJan 2:0686b2132556 53 double q2=0; // positie q2 in radialen
JornJan 2:0686b2132556 54 double q1_pos;
JornJan 2:0686b2132556 55 double q2_pos;
JornJan 2:0686b2132556 56
JornJan 2:0686b2132556 57 // EMG Input_k
JornJan 2:0686b2132556 58 double vx;
JornJan 2:0686b2132556 59 double vy;
JornJan 2:0686b2132556 60
JornJan 2:0686b2132556 61
JornJan 2:0686b2132556 62 void kinematica()
JornJan 2:0686b2132556 63 {
JornJan 2:0686b2132556 64 // encoders uitlezen
JornJan 2:0686b2132556 65 q1_puls = q1_enc.getPulses();
JornJan 2:0686b2132556 66 q1_pos = q1_puls*pulses2rad; // berekent positie q1 in radialen
JornJan 2:0686b2132556 67 q2_puls = q2_enc.getPulses();
JornJan 2:0686b2132556 68 q2_pos = q2_puls*pulses2rad; // berekent positie q2 in radialen
JornJan 2:0686b2132556 69
JornJan 2:0686b2132556 70 vx = 4*sin( 2*t.read() );
JornJan 2:0686b2132556 71 vy = 0;
JornJan 2:0686b2132556 72
JornJan 2:0686b2132556 73 q1 = ((-(L3 + L2*cos(q1_pos - q2_pos))/(L1*L2*cos(2*q1_pos - q2_pos) + L1*L3*cos(q1_pos)))*vx + ((L2*sin(q1_pos - q2_pos))/(L1*L2*cos(2*q1_pos - q2_pos) + L1*L3*cos(q1_pos)))*vy) * ts + q1_pos;
JornJan 2:0686b2132556 74 q2 = (((L3 - L1*sin(q1_pos) + L2*cos(q1_pos - q2_pos))/(L1*L2*cos(2*q1_pos - q2_pos) + L1*L3*cos(q1_pos))*vx) + ((L1*cos(q1_pos) - L2*sin(q1_pos - q2_pos))/(L1*L2*cos(2*q1_pos - q2_pos) + L1*L3*cos(q1_pos)))*vy) * ts + q2_pos;
JornJan 2:0686b2132556 75
JornJan 2:0686b2132556 76
JornJan 2:0686b2132556 77 ref1 = q1*rad2pulses;
JornJan 2:0686b2132556 78 ref2 = q2*rad2pulses;
JornJan 2:0686b2132556 79 }
JornJan 2:0686b2132556 80
JornJan 2:0686b2132556 81
JornJan 0:921402a95675 82 void controller()
JornJan 0:921402a95675 83 {
JornJan 2:0686b2132556 84 kinematica();
JornJan 0:921402a95675 85
JornJan 0:921402a95675 86 //PID 1
JornJan 0:921402a95675 87 error1_1 = ref1 - q1_puls;
JornJan 0:921402a95675 88 error_I_1 = error_I2_1 + ts*((error1_1 + error2_1)/2);
JornJan 0:921402a95675 89
JornJan 0:921402a95675 90 PI1 = Kp1*error1_1 + Ki1*error_I_1;
JornJan 0:921402a95675 91
JornJan 0:921402a95675 92 error2_1 = error1_1; // opslaan variabelen voor integraal onderdeel
JornJan 0:921402a95675 93 error_I2_1 = error_I_1;
JornJan 0:921402a95675 94
JornJan 0:921402a95675 95 //PID 2
JornJan 0:921402a95675 96 error1_2 = ref2 - q2_puls;
JornJan 0:921402a95675 97 error_I_2 = error_I2_2 + ts*((error1_2 + error2_2)/2);
JornJan 0:921402a95675 98
JornJan 0:921402a95675 99 PI2 = Kp2*error1_2 + Ki2*(error_I_2);
JornJan 0:921402a95675 100
JornJan 0:921402a95675 101
JornJan 0:921402a95675 102 if (PI1<=0)
JornJan 0:921402a95675 103 {
JornJan 0:921402a95675 104 motor1DirectionPin = 1;
JornJan 0:921402a95675 105 motor1MagnitudePin = fabs(PI1);
JornJan 0:921402a95675 106 }
JornJan 0:921402a95675 107 else if (PI1>0)
JornJan 0:921402a95675 108 {
JornJan 0:921402a95675 109 motor1DirectionPin = 0;
JornJan 0:921402a95675 110 motor1MagnitudePin = fabs(PI1);
JornJan 0:921402a95675 111 }
JornJan 0:921402a95675 112
JornJan 0:921402a95675 113 if (PI2>=0)
JornJan 0:921402a95675 114 {
JornJan 0:921402a95675 115 motor2DirectionPin = 1;
JornJan 1:7969189824ff 116 //motor2MagnitudePin = fabs(PI2);
JornJan 0:921402a95675 117 }
JornJan 0:921402a95675 118 else if(PI2<0)
JornJan 0:921402a95675 119 {
JornJan 0:921402a95675 120 motor2DirectionPin = 0;
JornJan 1:7969189824ff 121 //motor2MagnitudePin = fabs(PI2);
JornJan 0:921402a95675 122 }
JornJan 0:921402a95675 123
JornJan 0:921402a95675 124 if(n==500){
JornJan 2:0686b2132556 125 printf("PI1 = %f PI2 = %f\n\r", ref1, ref2);
JornJan 0:921402a95675 126 n=0;
JornJan 0:921402a95675 127 }
JornJan 0:921402a95675 128 else{
JornJan 0:921402a95675 129 n=n+1;
JornJan 0:921402a95675 130 }
JornJan 0:921402a95675 131
JornJan 0:921402a95675 132 }
JornJan 0:921402a95675 133
JornJan 0:921402a95675 134 void stapfunc()
JornJan 0:921402a95675 135 {
JornJan 0:921402a95675 136 if (ref1==0){
JornJan 0:921402a95675 137 ref1 = 1000;
JornJan 0:921402a95675 138 }
JornJan 0:921402a95675 139 else if (ref1==1000)
JornJan 0:921402a95675 140 {
JornJan 0:921402a95675 141 ref1=0;
JornJan 0:921402a95675 142 }
JornJan 0:921402a95675 143 }
JornJan 0:921402a95675 144
JornJan 0:921402a95675 145 int main()
JornJan 0:921402a95675 146 {
JornJan 0:921402a95675 147 pc.baud(115200);
JornJan 0:921402a95675 148 t.start();
JornJan 0:921402a95675 149 aansturing.attach_us(&controller, 1000);
JornJan 0:921402a95675 150 //stepres.attach(&stapfunc, 3);
JornJan 0:921402a95675 151
JornJan 0:921402a95675 152 while(true){
JornJan 0:921402a95675 153
JornJan 0:921402a95675 154 }
JornJan 0:921402a95675 155
JornJan 0:921402a95675 156
JornJan 0:921402a95675 157 }