pmsm-control/test_sw/main_pmsm.c

   1 /**
   2  * \file main_pmsm.c
   3  * \author Martin Prudek
   4  * \brief Mainfile pro pmsm control.
   5  */
   6
   7 #ifndef NULL
   8 #define NULL (void*) 0
   9 #endif /*NULL*/
  10
  11
  12 #include <stdlib.h>     /*exit*/
  13 #include <signal.h>     /*signal handler Ctrl+C*/
  14 #include <stdio.h>      /*printf*/
  15 #include <sched.h>      /*sheduler*/
  16 #include <unistd.h>     /*usleep*/
  17 #include <pthread.h>    /*threads*/
  18
  19 #include "rpin.h"       /*gpclk*/
  20 #include "rp_spi.h"     /*spi*/
  21 #include "misc.h"       /*structure for priorities*/
  22 #include "pxmc_sin_fixed.h"     /*to test sin commutation */
  23
  24
  25
  26 #define PRUM_PROUD      2061
  27 #define PRUM_SOUC       6183
  28 #define MAX_DUTY        128
  29 #define PID_P           0.1
  30
  31 #define PRIOR_KERN      50
  32 #define PRIOR_HIGH      49
  33 #define PRIOR_LOW       20
  34
  35 #define THREAD_SHARED   0
  36 #define INIT_VALUE      0       /*init value for semaphor*/
  37
  38
  39 #define PXMC_SIN_FIX_TAB_BITS 9
  40 #define PXMC_SIN_FIX_IDX_SLR  23
  41 #define PXMC_SIN_FIX_XD_MASK  0x007fffff
  42 #define PXMC_SIN_FIX_XD_SLR   8
  43 #define PXMC_SIN_FIX_A_MASK   0xffffc000
  44 #define PXMC_SIN_FIX_B_SLL    19
  45 #define PXMC_SIN_FIX_B_SAR    16
  46 #define PXMC_SIN_FIX_B_XD_SAR 6
  47 #define PXMC_SIN_FIX_ZIC_MASK 0x00002000
  48 #define PXMC_SIN_FIX_ZIC_BIT  13
  49
  50 #define PXMC_SIN_FIX_PI2      0x40000000
  51 #define PXMC_SIN_FIX_2PI3     0x55555555
  52
  53 struct sigaction sighnd; /*struktura pro signal handler*/
  54 struct rpi_in data;
  55 struct rpi_state{
  56         uint8_t test;
  57         uint16_t pwm1, pwm2, pwm3;
  58         uint16_t t_pwm1, t_pwm2, t_pwm3;
  59         char commutate;
  60         int duty;                       /* duty cycle of pwm */
  61         uint16_t index_dist;            /* distance to index position */
  62         unsigned char index_ok;
  63         uint32_t tf_count;              /*number of transfer*/
  64         int desired_pos;                /* desired position */
  65 }rps;
  66
  67
  68 /**
  69  * \brief Initilizes GPCLK.
  70  */
  71 int clk_init()
  72 {
  73         initialise(); /*namapovani gpio*/
  74         initClock(PLLD_500_MHZ, 10, 0);
  75         gpioSetMode(4, FSEL_ALT0);
  76         return 0;
  77 }
  78 /*
  79  * \brief Terminates GPCLK.
  80  */
  81
  82 inline void clk_disable(){
  83         termClock(0);
  84 }
  85
  86 /**
  87  * \brief Signal handler pro Ctrl+C
  88  */
  89 void sighnd_fnc(){
  90         spi_disable();
  91         clk_disable();
  92         /*muzeme zavrit semafor*/
  93         sem_destroy(&thd_par_sem);
  94         printf("\nprogram bezpecne ukoncen\n");
  95         exit(0);
  96 }
  97
  98 void substractOffset(struct rpi_in* data, struct rpi_in* offset){
  99         data->pozice_raw=data->pozice;
 100         data->pozice-=offset->pozice;
 101         return;
 102 }
 103 /*
 104  * pocita procentualni odchylku od prumerneho proudu
 105  */
 106 float diff_p(float value){
 107         return ((float)value-PRUM_PROUD)*100/PRUM_PROUD;
 108 }
 109 /*
 110  * pocita procentualni odchylku od prumerneho souctu proudu
 111  */
 112 float diff_s(float value){
 113         return ((float)value-PRUM_SOUC)*100/PRUM_SOUC;
 114 }
 115 /*
 116  * tiskne potrebna data
 117  */
 118 void printData(){
 119         struct rpi_in data_p;
 120         struct rpi_state s;     /*state*/
 121         float cur0, cur1, cur2;
 122         int i;
 123         /* copy the data */
 124         sem_wait(&thd_par_sem);
 125         data_p = data;
 126         s=rps;
 127         sem_post(&thd_par_sem);
 128
 129         if (data_p.adc_m_count){
 130                 cur0=data_p.ch0/data_p.adc_m_count;
 131                 cur1=data_p.ch1/data_p.adc_m_count;
 132                 cur2=data_p.ch2/data_p.adc_m_count;
 133         }
 134         for (i = 0; i < 16; i++) {
 135                         if (!(i % 6))
 136                                 puts("");
 137                         printf("%.2X ", data_p.debug_rx[i]);
 138         }
 139         puts("");
 140         printf("\npozice=%d\n",(int32_t)data_p.pozice);
 141         printf("chtena pozice=%d\n",s.desired_pos);
 142         printf("transfer count=%u\n",s.tf_count);
 143         printf("raw_pozice=%d\n",(int32_t)data_p.pozice_raw);
 144         printf("raw_pozice last12=%u\n",(data_p.pozice_raw&0x0FFF));
 145         printf("index position=%u\n",data_p.index_position);
 146         printf("hal1=%d, hal2=%d, hal3=%d\n",data_p.hal1,data_p.hal2,data_p.hal3);
 147         printf("en1=%d, en2=%d, en3=%d (Last sent)\n",!!(0x40&s.test),!!(0x20&s.test),!!(0x10&s.test));
 148         printf("shdn1=%d, shdn2=%d, shdn3=%d (L.s.)\n",!!(0x08&s.test),!!(0x04&s.test),!!(0x02&s.test));
 149         printf("PWM1=%u(L.s.)\n",s.pwm1);
 150         printf("PWM2=%u(L.s.)\n",s.pwm2);
 151         printf("PWM3=%u(L.s.)\n",s.pwm3);
 152         printf("distance to index=%u\n",s.index_dist);
 153         printf("T_PWM1=%u T_PWM2=%u T_PWM3=%u\n",s.t_pwm1,s.t_pwm2, s.t_pwm3);
 154         printf("Pocet namerenych proudu=%u\n",data_p.adc_m_count);
 155         printf("(pwm1) (ch1)=%d (avg=%4.0f) (%2.2f%%)\n",data_p.ch1,cur1,diff_p(cur1));
 156         printf("(pwm2) (ch2)=%d (avg=%4.0f)(%2.2f%%)\n",data_p.ch2,cur2,diff_p(cur2));
 157         printf("(pwm3) (ch0)=%d (avg=%4.0f)(%2.2f%%)\n",data_p.ch0,cur0,diff_p(cur0));
 158         printf("soucet prumeru=%5.0f (%2.2f%%)\n",cur0+cur1+cur2,diff_s(cur0+cur1+cur2));
 159         printf("duty=%d\n",s.duty);
 160         if (s.index_ok) printf("index ok\n");
 161         if (s.commutate) printf("commutation in progress\n");
 162 }
 163 void prepare_tx(uint8_t * tx){
 164
 165         /*Data format:
 166          * tx[4] - bity 95 downto 88 - bits that are sent first
 167          * tx[5] - bity 87 downto 80
 168          * tx[6] - bity 79 downto 72
 169          * tx[7] - bity 71 downto 64
 170          * tx[8] - bity 63 downto 56
 171          * tx[9] - bity 55 downto 48
 172          * tx[10] - bity 47 downto 40
 173          * tx[11] - bity 39 downto 32
 174          * tx[12] - bity 31 downto 24
 175          * tx[13] - bity 23 downto 16
 176          * tx[14] - bity 15 downto 8
 177          * tx[15] - bity 7 downto 0
 178          *
 179          * bit 95 - ADC reset
 180          * bit 94 - enable PWM1
 181          * bit 93 - enable PWM2
 182          * bit 92 - enable PWM3
 183          * bit 91 - shutdown1
 184          * bit 90 - shutdown2
 185          * bit 89 - shutdown3
 186          *      .
 187          *      .
 188          *      .
 189          * bits 66 .. 56 - match PWM1
 190          * bits 55 .. 45 - match PWM2
 191          * bit 11,12 - Unused
 192          * bits 42 .. 32  - match PWM3
 193          */
 194
 195
 196         uint16_t tmp;
 197
 198         /* keep the cap*/
 199
 200         if (rps.pwm1>2047) rps.pwm1=2047;
 201         if (rps.pwm2>2047) rps.pwm2=2047;
 202         if (rps.pwm3>2047) rps.pwm3=2047;
 203
 204         tx[0]=rps.test; /*bit 94 - enable PWM1*/
 205
 206         /*pwm1*/
 207         tx[7]=(tx[7] & 0xF8) | (0x07 & ((uint8_t*)&rps.pwm1)[1]); /*MSB*/
 208         tx[8]=((uint8_t*)&rps.pwm1)[0]; /*LSB*/
 209
 210         /*pwm2*/
 211         tmp=rps.pwm2;
 212         tmp<<=5;
 213         tx[9]=((uint8_t*)&tmp)[1]; /*MSB*/
 214         tx[10]=(tx[10] & 0x1F) | (0xE0 & ((uint8_t*)&tmp)[0]); /*LSB*/
 215
 216         /*pwm3*/
 217         tx[10]=(tx[10] & 0xF8) | (0x07 & ((uint8_t*)&rps.pwm3)[1]); /*MSB*/
 218         tx[11]=((uint8_t*)&rps.pwm3)[0]; /*LSB*/
 219
 220
 221 }
 222 /**
 223  * Funkce pravidelne vypisuje posledni zjistenou pozici lokalniho motoru
 224  */
 225 void * pos_monitor(void* param){
 226         set_priority(param);            /*set priority*/
 227         while(1){
 228                 printData();
 229                 usleep(1000000);        /*1 Hz*/
 230         }
 231         return (void*)0;
 232 }
 233
 234 inline
 235 int sin_commutator(int duty){
 236         #define DEGREE_60        715827883
 237         #define DEGREE_120      1431655765
 238         #define DEGREE_180      2147483648
 239         #define DEGREE_240      2863311531
 240         #define DEGREE_300      3579139413
 241         uint32_t j,pos;
 242         int32_t sin;
 243         uint16_t pwm;
 244         pos=rps.index_dist*4294967;
 245         if (duty>=0){   /*clockwise rotation*/
 246                 /* 1st phase */
 247                 sin = pxmc_sin_fixed_inline(pos+DEGREE_240,10); /*10+1 bity*/ /*-120*/
 248                 pwm=0;
 249                 for(j=0;j!=11;j++){
 250                         /* multiplicate as if maximum sinus value was unity */
 251                         pwm+=(!(sin & 0x10000000))*(((1 << j) & sin)>>j)*(duty>>(10-j));
 252                 }
 253                 rps.pwm1=pwm;
 254
 255                 /* 2nd phase */
 256                 sin = pxmc_sin_fixed_inline(pos+DEGREE_120,10); /*10+1 bity*/ /*-240*/
 257                 pwm=0;
 258                 for(j=0;j!=11;j++){
 259                         /* multiplicate as if maximum sinus value was unity */
 260                         pwm+=(!(sin & 0x10000000))*(((1 << j) & sin)>>j)*(duty>>(10-j));
 261                 }
 262                 rps.pwm2=pwm;
 263
 264                 /* 3rd phase */
 265                 sin = pxmc_sin_fixed_inline(pos,10); /*10+1 bity*/
 266                 pwm=0;
 267                 for(j=0;j!=11;j++){
 268                         /* multiplicate as if maximum sinus value was unity */
 269                         pwm+=(!(sin & 0x10000000))*(((1 << j) & sin)>>j)*(duty>>(10-j));
 270                 }
 271                 rps.pwm3=pwm;
 272         }else{
 273                 duty=-duty;
 274
 275                 /* 1st phase */
 276                 sin = pxmc_sin_fixed_inline(pos+DEGREE_60,10); /*10+1 bity*/ /*-300*/
 277                 pwm=0;
 278                 for(j=0;j!=11;j++){
 279                         /* multiplicate as if maximum sinus value was unity */
 280                         pwm+=(!(sin & 0x10000000))*(((1 << j) & sin)>>j)*(duty>>(10-j));
 281                 }
 282                 rps.pwm1=pwm;
 283
 284                 /* 2nd phase */
 285                 sin = pxmc_sin_fixed_inline(pos+DEGREE_300,10); /*10+1 bity*/ /*-60-*/
 286                 pwm=0;
 287                 for(j=0;j!=11;j++){
 288                         /* multiplicate as if maximum sinus value was unity */
 289                         pwm+=(!(sin & 0x10000000))*(((1 << j) & sin)>>j)*(duty>>(10-j));
 290                 }
 291                 rps.pwm2=pwm;
 292
 293                 /* 3rd phase */
 294                 sin = pxmc_sin_fixed_inline(pos+DEGREE_180,10); /*10+1 bity*/ /*-180*/
 295                 pwm=0;
 296                 for(j=0;j!=11;j++){
 297                         /* multiplicate as if maximum sinus value was unity */
 298                         pwm+=(!(sin & 0x10000000))*(((1 << j) & sin)>>j)*(duty>>(10-j));
 299                 }
 300                 rps.pwm3=pwm;
 301         }
 302         return 0;
 303 }
 304 /*
 305  * \brief
 306  * Test function to be placed in controll loop.
 307  * Switches PWM's at point where they produce same force.
 308  * This points are found thanks to IRC position,
 309  */
 310 inline
 311 void simple_ind_dist_commutator(int duty){
 312         if (duty>=0){ /* clockwise - so that position increase */
 313                 /* pwm3 */
 314                 if ((rps.index_dist>=45 && rps.index_dist<=373) ||
 315                 (rps.index_dist>=1048 && rps.index_dist<=1377)){
 316                         rps.pwm1=0;
 317                         rps.pwm2=0;
 318                         rps.pwm3=duty;
 319                         /* pwm1 */
 320                 }else if ((rps.index_dist>=373 && rps.index_dist<=711) ||
 321                 (rps.index_dist>=1377 && rps.index_dist<=1711)){
 322                         rps.pwm1=duty;
 323                         rps.pwm2=0;
 324                         rps.pwm3=0;
 325                         /* pwm2 */
 326                 }else if ((rps.index_dist>=0 && rps.index_dist<=45) ||
 327                 (rps.index_dist>=711 && rps.index_dist<=1048) ||
 328                 (rps.index_dist>=1711 && rps.index_dist<=1999)){
 329                         rps.pwm1=0;
 330                         rps.pwm2=duty;
 331                         rps.pwm3=0;
 332                 }
 333         }else{  /*counter-clockwise - position decrease */
 334                 /* pwm3 */
 335                 if ((rps.index_dist>=544 && rps.index_dist<=881) ||
 336                 (rps.index_dist>=1544 && rps.index_dist<=1878)){
 337                         rps.pwm1=0;
 338                         rps.pwm2=0;
 339                         rps.pwm3=-duty;
 340                         /* pwm1 */
 341                 }else if ((rps.index_dist>=0 && rps.index_dist<=211) ||
 342                 (rps.index_dist>=881 && rps.index_dist<=1210) ||
 343                 (rps.index_dist>=1878 && rps.index_dist<=1999)){
 344                         rps.pwm1=-duty;
 345                         rps.pwm2=0;
 346                         rps.pwm3=0;
 347                         /* pwm2 */
 348                 }else if ((rps.index_dist>=211 && rps.index_dist<=544) ||
 349                 (rps.index_dist>=1210 && rps.index_dist<=1544)){
 350                         rps.pwm1=0;
 351                         rps.pwm2=-duty;
 352                         rps.pwm3=0;
 353                 }
 354         }
 355 }
 356 /*
 357  * \brief
 358  * Test function to be placed in controll loop.
 359  * Switches PWM's at point where they produce same force
 360  */
 361 inline void simple_hall_commutator(int duty){
 362         if (duty>=0){ /* clockwise - so that position increase */
 363                 /* pwm3 */
 364                 if (data.hal2 && !data.hal3){
 365                         rps.pwm1=0;
 366                         rps.pwm2=0;
 367                         rps.pwm3=duty;
 368                         /* pwm1 */
 369                 }else if (data.hal1 && !data.hal2){
 370                         rps.pwm1=duty;
 371                         rps.pwm2=0;
 372                         rps.pwm3=0;
 373                         /* pwm2 */
 374                 }else if (!data.hal1 && data.hal3){
 375                         rps.pwm1=0;
 376                         rps.pwm2=duty;
 377                         rps.pwm3=0;
 378                 }
 379         }else{  /*counter-clockwise - position decrease */
 380                 /* pwm3 */
 381                 if (!data.hal2 && data.hal3){
 382                         rps.pwm1=0;
 383                         rps.pwm2=0;
 384                         rps.pwm3=-duty;
 385                         /* pwm1 */
 386                 }else if (!data.hal1 && data.hal2){
 387                         rps.pwm1=-duty;
 388                         rps.pwm2=0;
 389                         rps.pwm3=0;
 390                         /* pwm2 */
 391                 }else if (data.hal1 && !data.hal3){
 392                         rps.pwm1=0;
 393                         rps.pwm2=-duty;
 394                         rps.pwm3=0;
 395                 }
 396         }
 397 }
 398 /**
 399  * Funkce pravidelne vycita data z motoru
 400  */
 401 inline void comIndDist(){
 402         rps.index_dist=0x0FFF & (data.pozice_raw - data.index_position);
 403         /*
 404          * if distance is bigger than 2047, the distance underflown
 405          * -> if 12th bit is set, substract 2096
 406          */
 407         rps.index_dist-=((rps.index_dist & 0x0800)>>11)*2096;
 408 }
 409 /*
 410  * \brief
 411  * Very simple PID regulator.
 412  * Now only with P-part so that the error doesnt go to zero.
 413  * TODO: add anti-wind up and I and D parts
 414  */
 415 inline void pid(){
 416         int duty_tmp;
 417         duty_tmp = PID_P*(rps.desired_pos - (int32_t)data.pozice);
 418         if (duty_tmp>MAX_DUTY){
 419                 rps.duty=MAX_DUTY;
 420         }else if (duty_tmp<-MAX_DUTY){
 421                 rps.duty=-MAX_DUTY;
 422         }else{
 423                 rps.duty = duty_tmp;
 424         }
 425 }
 426 /*
 427  * \brief
 428  * Feedback loop.
 429  * TODO: replace usleep with real-time wait
 430  *      measure times
 431  */
 432 void * read_data(void* param){
 433         int i;
 434         struct rpi_in pocatek;
 435         uint8_t tx[]={0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ;
 436         char first=1;
 437         uint16_t last_index;                            /*we have index up-to date*/
 438         set_priority(param);                            /*set priority*/
 439         pocatek = spi_read(tx);
 440                 while(1){
 441                         prepare_tx(tx);                 /*save the data to send*/
 442                         sem_wait(&thd_par_sem);         /*---take semaphore---*/
 443                         data = spi_read(tx);            /*exchange data*/
 444                         /*subtract initiate postion */
 445                         rps.tf_count++;
 446                         substractOffset(&data,&pocatek);
 447                         comIndDist();
 448                         if (!rps.index_ok){
 449                                 if (first){
 450                                         last_index=data.index_position;
 451                                         first=0;
 452                                 }else if (last_index!=data.index_position){
 453                                         rps.index_ok=1;
 454                                 }
 455                         }
 456                         pid();
 457                         if (rps.index_ok && rps.commutate){
 458                                 /*simple_ind_dist_commutator(rps.duty);*/
 459                                 sin_commutator(rps.duty);
 460                         }else if(!rps.index_ok && rps.commutate){
 461                                 simple_hall_commutator(rps.duty);
 462                         }
 463                         sem_post(&thd_par_sem);         /*--post semaphore---*/
 464                         usleep(1000);                           /*1kHz*/
 465                 }
 466 }
 467
 468 /**
 469  * \brief Main function.
 470  */
 471
 472 int main(){
 473         unsigned int tmp;
 474
 475         /*nastaveni priorit vlaken*/
 476         struct thread_param tsp;
 477         tsp.sch_policy = SCHED_FIFO;
 478
 479         /*nastaveni signalu pro vypnuti pomoci Ctrl+C*/
 480         sighnd.sa_handler=&sighnd_fnc;
 481         sigaction(SIGINT, &sighnd, NULL );
 482
 483         clk_init();             /* inicializace gpio hodin */
 484         spi_init();             /* iniicializace spi*/
 485
 486         /*semafor pro detekci zpracovani parametru vlaken*/
 487         sem_init(&thd_par_sem,THREAD_SHARED,INIT_VALUE);
 488
 489         /*vlakna*/
 490         pthread_t tid;                  /*identifikator vlakna*/
 491         pthread_attr_t attr;            /*atributy vlakna*/
 492         pthread_attr_init(&attr);       /*inicializuj implicitni atributy*/
 493
 494
 495
 496         /*ziskavani dat z motoru*//*vysoka priorita*/
 497         tsp.sch_prior = PRIOR_HIGH;
 498         pthread_create(&tid, &attr, read_data, (void*)&tsp);
 499
 500         /*vypisovani lokalni pozice*//*nizka priorita*/
 501         tsp.sch_prior = PRIOR_LOW;
 502         sem_wait(&thd_par_sem);
 503         pthread_create(&tid, &attr, pos_monitor, (void*)&tsp);
 504
 505
 506
 507         /*
 508          * Note:
 509          * pri pouziti scanf("%u",&simple_hall_duty); dochazelo
 510          * k preukladani hodnot na promenne test. Dost divne.
 511          */
 512         while (1){
 513                 scanf("%u",&tmp);
 514                 printf("volba=%u\n",tmp);
 515                 switch (tmp){
 516                 case 1:
 517                         scanf("%u",&tmp);
 518                         sem_wait(&thd_par_sem);
 519                         rps.pwm1=tmp&0xFFF;
 520                         sem_post(&thd_par_sem);
 521                         break;
 522                 case 2:
 523                         scanf("%u",&tmp);
 524                         sem_wait(&thd_par_sem);
 525                         rps.pwm2=tmp&0xFFF;
 526                         sem_post(&thd_par_sem);
 527                         break;
 528                 case 3:
 529                         scanf("%u",&tmp);
 530                         sem_wait(&thd_par_sem);
 531                         rps.pwm3=tmp&0xFFF;
 532                         sem_post(&thd_par_sem);
 533                         break;
 534                 case 4:
 535                         scanf("%u",&tmp);
 536                         sem_wait(&thd_par_sem);
 537                         rps.test=tmp&0xFF;
 538                         sem_post(&thd_par_sem);
 539                         break;
 540                 case 5:
 541                         sem_wait(&thd_par_sem);
 542                         rps.commutate=!rps.commutate;
 543                         /* switch off pwms at the end of commutation */
 544                         rps.pwm1&=rps.commutate*0xFFFF;
 545                         rps.pwm2&=rps.commutate*0xFFFF;
 546                         rps.pwm3&=rps.commutate*0xFFFF;
 547                         sem_post(&thd_par_sem);
 548                         break;
 549                 case 6:
 550                         scanf("%d",&tmp);
 551                         sem_wait(&thd_par_sem);
 552                         rps.duty=tmp;
 553                         sem_post(&thd_par_sem);
 554                         break;
 555                 case 7:
 556                         scanf("%d",&tmp);
 557                         sem_wait(&thd_par_sem);
 558                         rps.desired_pos=tmp;
 559                         sem_post(&thd_par_sem);
 560                         break;
 561
 562                 default:
 563                         break;
 564                 }
 565
 566         }
 567         return 0;
 568 }
 569