• The Stepper Motor Driver Board

    agp.cooper12/14/2018 at 22:52 0 comments

    The Stepper Motor Driver Board

    Here is the Stepper Motor Driver  Board:

    Almost no documentation in the Internet on this board! It is cheap and I have used it before.

    Mapping/tracing the Nano pins I get:

    • RX - Serial Receive (pin provided but not used)
    • TX - Serial Transmit (pin provided but not used)
    • D2 - X Direction
    • D3 - Y Direction
    • D4 - Z Direction
    • D5 - X Step
    • D6 - Y Step
    • D7 - Z Step
    • D8 - Enable
    • D9 - +/- X End Stops (pin provided but not used)
    • D10 - +/-Y End Stops (pin provided but not used)
    • D11 - +/-Z End Stops (pin provided but not used)
    • D12 - Free (pin provided)
    • D13 - Free (pin provided)
    • A0 - Abort (pin provided but not used)
    • A1 - Hold (pin provided but not used)
    • A2 - Resume (pin provided but not used)
    • A3 - Coolant Enable (pin provided but not used)
    • A4 - SDA (pin provided but not used)
    • A5 - SCL (pin provided but not used)
    • A6 - Free (pin provided)
    • A7 - Free (pin provided)

    There are also:

    • +5v0 and ground pins
    • +3v3 and ground pins (uses an on board regulator)
    • Reset button (for Nano)
    • E-Stop pins (same as reset)
    • The power plug (Motor Power 8-12v):
      • Powers the Nano and via a shunt powers the A4988 stepper motor boards.
      • Otherwise the A4988 boards and the Nano are powered separately.
    • Under each A4988 stepper motor board are pins to shunt (S) the micro-stepping:
      • MS1 MS2 MS3 (up view)
      •      -        -       -  Full Step
      •      S       -       -  Half Step
      •      -       S       -    1/4 Step
      •      S      S       -   1 /8 Step
      •      S      S      S  1/16 Step   

    WATCH OUT!

    One fault with these boards is that the shunts to control the micro stepping does not work. I added these 10k pull-ups to fix the problem:

    Grbl

    Pre-processing "cartesian" gCode with gCode2SCARA allows the use of Grbl for a SCARA arm (UVZ type). This what I would recommend to get started.

    Not Grbl

    Is some code that I wrote that will also do the job.

    Motion Controller

    Is some more code that I wrote (but needs some work of the interpreter) that could also do the job. I intend to finish this code off to work with gCodeFilter and gCode2SCARA at some point of time.

    gCode2SCARA Command Line Options

    If you have not realised, gCode2SCARA is a command line program program designed to be used via a batch file. Here are the current options:

    • -i       Input File

    • -o     Output File

    • -t      XY Space Tolerance (mm)

    • -s     UV Space Tolerance (full steps)

    • -xc   gCode X Offset (mm)

    • -yc   gCode Y Offset (mm)

    • -L1   Arm 1 Length (mm)

    • -L2  Arm 2 Length (mm)

    • -M1  Motor 1 Steps

    • -M2 Motor 2 Steps

    This project could benefit from a "Visual Windows" programming approach, so probably time to look at what is available for Linux that is simple to use.

    AlanX

  • gCode2SCARA

    agp.cooper12/13/2018 at 00:08 0 comments

    gCode to SCARA

    The code is basically done, some check work remains.

    Here is the XY space (filtered from 45984 points down to 2029 points, i.e. -t 0.01):

    And here is the UV Space (-s 0.01 and 2029 points, no additional inserted point):

    The UV plot (i.e. the output) however does not consider the step resolution (i.e. UV space is floating point). The plot seems to be in the wrong quadrant (it is possible as I have changed the reverse kinematics mathematics).

    I will will have to create another plot window to show the effect of step resolution.

    So just some tidy up remains.

    Code Checks Out

    Yesterday I was helping my partner clean her house after the floor boards were sanded and varnished. So no work!

    The code checks out, need to do something better for when the arms will not reach, currently just set the location to (0,0). Better to make it obvious in the imagery.

    Checked that the reverse kinematics works in all quadrants. Here is an example with the arm movements shown:

    Added a Bressenham trace (yellow/green), in this case for M1 and M2 equal to 400 steps each (note, the image has been shifted to exceed the limits of the reach of the arms):

    Here is the case for M1 and M2 equal to 1600 steps (1/4 micro-stepping):

    And finally M12800 (1/32 micro stepping):

    Note that G0 (rapid movement is an arc rather than a straight line).

    So other than modifications required (discovered) when I actually use the code on a machine, it is done.

    AlanX

  • Short Cut - Converting gCode Cartesian to Polar

    agp.cooper12/12/2018 at 01:44 0 comments

    Short Cut

    Its a rather "long road" to write a Grbl replacement (perhaps later), but converting gCode from XYZ to UVZ is not. I have seen that someone has already done this on the Internet (sorry I have lost the link). So basically after conversion you can use Grbl instead.

    Perhaps I am wrong but Grbl does not appear to handle UVZ coordinate systems.

    I already have a gCode Filter project that can be co-opted into service.

    After filtering the gCode to remove "redundant" segments (i.e. the main function of gCode Filter), I only need map the UVZ space and to insert the critical points.

    I have made a start, set up two graphical results windows, one for XY space and the other for UV space. Still about a days coding left to complete (wishful thinking perhaps?).

    Bearings

    One of the issues for construction of a SCARA is the bearing. Nearly all the projects I have looked at are DIY. It would be nice to get something reasonalbe (lazy susan bearings are pretty bad!) off the self.

    For a big project a car one piece stub axle looks like an option (~AUD50):

    After a bit of searching I found good candidate, a bicycle front wheel hub:

    Are you excited yet? AlanX

  • Divide and Conquer Wins

    agp.cooper12/11/2018 at 12:44 0 comments

    Divide and Conquer (Option 1)

    After playing with the Options presented in my first post, I have decided Option 1 is the best. Why? It  should have the least number of the expensive reverse kinematic calculations. Consider a very coarse tolerance (=128):

    Note that the UV space (blue lines) has only 6 points (the straight lines between them are Bressenham lines).

    Now the medium tolerance case (=32):

    Here there is only 13 reverse kinematic calculations.

    Now the very fine tolerance (=8):

    Here there is only 25 reverse kinematic calculations.

    Finally ultra fine tolerance (=1), has 67 calculations:

    Fast Trigonometry

    Although the number of reverse kinematic calculations are a minimum, there still may be advantage is fast (integer?) methods. The CORDIC algorithm comes to mind, especially for the ATAN2() and sqrt() functions. ArcSin/ArcCos should also be possible.

    Here is my version of Rect2Polar (i.e ATAN2() and sqrt()) using short (i.e int16) integers:

    // Include libraries
    #include <stdio.h>
    #include <stdlib.h>
    #include <stdint.h>
    #include <string.h>
    #include <math.h>
    #include <stdbool.h>
    
    void rect2polar(short X,short Y,short *R,short *A)
    {
      // ATAN_Table is the values of ATAN(1/(2^i))*400*32/2/Pi micro-steps
      short ATAN[12]={1600,945,499,253,127,64,32,16,8,4,2,1};
      short i,Xnew,Ynew;
    
      (*A)=0;
      if (X<0) {
        (*A)=200*32; // 180 degrees
        X=-X;
        Y=-Y;
      } else if (Y<0) {
        (*A)=400*32; // =360 degrees
      }
      for (i=0;i<12;i++) {
        if (Y<0) {
          // Rotate CCW
          Xnew=X-(Y>>i);
          Ynew=Y+(X>>i);
          (*A)=(*A)-ATAN[i];
        } else {
          // Rotate CW
          Xnew=X+(Y>>i);
          Ynew=Y-(X>>i);
          (*A)=(*A)+ATAN[i];
        }
        X=Xnew;
        Y=Ynew;
      }
      // Adjust for gain (=X*0.607252935)
      X=X>>1;
      (*R)=(((((((((((((((((X>>2)+X)>>2)+X)>>1)+X)>>1)+X)>>2)+X)>>1)+X)>>2)+X)>>1)+X)>>3)+X;
    }
    int main(void) {
      short A,R,X,Y;
    
      X=-439*32;
      Y=-439*32;
      // R=621*32 max
      rect2polar(X,Y,&R,&A);
      printf("Rect2Polar x %8.3f y %8.3f ang %8.3f hyp %8.3f\n",X/32.0,Y/32.0,A*360.0/400.0/32.0,R/32.0);
      printf("atan2()    x %8.3f y %8.3f ang %8.3f hyp %8.3f\n",X/32.0,Y/32.0,360+atan2(Y,X)*180/M_PI,sqrt(X*X+Y*Y)/32.0);
    
      return 0;
    }

    Here is the result:

    $ ./ShortCAtan
    Rect2Polar x -439.000 y -439.000 ang  225.028 hyp  620.812
    atan2()        x -439.000 y -439.000 ang  225.000 hyp  620.840
    $

    So about  4 digits of accuracy with 16 bit integers.

    Reverse Kinematica

    I reworked the mathematics to suit CORDIC:

    • R=sqrt(X*X+Y*Y);
    • A0=atan2(X,Y);
    • A1=A0-PI/2+asin((L1*L1+R*R-L2*L2)/L1/R/2);
    • A2=PI/2-asin((R*R-L1*L1-L2*L2)/L1/L2/2);
    • X=cos(A1)*L1+cos(A1+A2)*L2;
    • Y=sin(A1)*L1+sin(A1+A2)*L2;

    ArcSin CORDIC

    Here is my ArcSin CORDIC:

    // This version works properly but has a multiplication in the main loop 
    short arcSin(short R, short S) {
      // ATAN_Table is the values of ATAN(1/(2^i))*400*32/2/Pi micro-steps
      short ATAN[12]={1600,945,499,253,127,64,32,16,8,4,2,1};
      // SEC[i]=(int)((1/COS(ATAN(1/(2^i)))-1)*(2>>14)+0.5);
      short SEC[15]={6787,1935,505,128,32,8,2,1,1,1,1,1,1,1,1};
      short i,Rnew,Y,Ynew,A;
    
      Y=0;
      A=0;
      for (i=0;i<=11;i++) {
        if ((Y<S)&&(A<100*32)) { // 90 degrees
          // Rotate CCW
          Ynew=Y+(R>>i);
          Rnew=R-(Y>>i);
          A=A+ATAN[i];
        } else {
          // Rotate CW
          Ynew=Y-(R>>i);
          Rnew=R+(Y>>i);
          A=A-ATAN[i];
        }
        Y=Ynew;
        R=Rnew;
        S=S+((S*SEC[i])>>14);
      }
      return A;
    }

     The down side of my ArcSin CORDIC is the multiplication inside the main loop. I remember I had to do the multiplication to avoid instability.

    AlanX

  • SCARA Reverse Kinematics

    agp.cooper12/10/2018 at 02:37 0 comments

    SCARA Model

    Here is an example of the type of SCARA that I want to model:

    Two X and Y arms and linear Z (say):

    L1 = 200mm

    L2 = 150mm

    Z = 150mm

    Two stepper motors:

    M1 = 400 steps (assume direct drive for the moment)

    M2 = 400 steps(assume direct drive for the moment)

    Also I want to use an Arduino Nano for the motion controller.

    Forward Kinematics

    Ignoring the Z axis for the time being, the forward kinematics represented by (U,V) to (X,Y), where (U,V) are polar motor steps and (X,Y) are the "real world" cartesian coordinates:

    • U is the polar (angular) position (steps) for motor 1
    • V is the polar (angular) position (steps) for motor 2
    • X is the "x cartesian" position in mm
    • Y is the "y cartesian" position in mm

    The forward mathematics is:

    • X=sin(2*Pi*U/M1)+sin(2*Pi*(U/M1+V/M2))
    • Y=cos(2*Pi*U/M1)+cos(2*Pi*(U/M1+V/M2))

    ReverseKinematics

    Ignoring the Z axis for the time being, the reverse kinematics represented by (X,Y) to (U,V):

    • R=sqrt(X*X+Y*Y)
    • U=atan2(X,Y)-acos((L1*L1+R*R-L2*L2)/L1/R/2)
    • V=PI-acos((L1*L1+L2*L2-R*R)/L1/L2/2)

    or

    • V=acos((R*R-L1*L1-L2*L2)/L1/L2/2)

    I have reworked the equations to:

    • R=sqrt(X*X+Y*Y)

    • CosA=(L1*L1+R*R-L2*L2)

    • SinA=sqrt((2*L1*R-CosA)*(2*L1*R+CosA))

    • U=atan2(Y,X)-atan2(SinA,CosA)

    • CosA=(R*R-L1*L1-L2*L2)

    • SinA=sqrt((2*L1*L2-CosA)*(2*L1*L2+CosA))

    • V=atan2(SinA,CosA);

    The reason is that sqrt() and atan2() is a well known for embedded applications.

    Problems with Reverse Kinematics

    Two issues:

    1. The equations assume either left or right elbow operation. The mathematics is right elbow at the moment.
    2. The cartesian coordinates do not map neatly with polar coordinates. That is, x and y positions do not usually match (map) an exact stepper motor position.

    For example:

    Note how straight lines in XY space are curves in UV space, and motor step resolution (in UV Space) maps into jagged steps in XY space.

    (The plotted square is 125 mm x 125 mm, the arms are 200 mm and 150 mm, and the motors are both 400 steps.)

    Okay, the third problem is the need for gearing as the motor step resolution is too low.

    ---

    Mapping Options

    I have tried three methods for accurate XY space to UV space.

    Option 1

    Segment the UV space curves into straight lines (by inserting points) using a tolerance (i.e. if the tolerance is not exceeded then a point need not be inserted):

    This works well for feeding the mapped UV space directly into a Bressenham's Line Algorithm (i.e. Grbl).

    This option could be pre-processing or on-the-fly (i.e. on the motion controller).

    Option 2

    Nearest UV space point mapping (path search):

    The method seems to be efficient (with look up tables) and path search algorithm can certainly be improved. This algorithm would replace the Bressenham's Line Algorithm (i.e. exports to the motor controller hardware directly).

    This option would be on-the-fly (i.e. replaces Grbl).

    Option 3

    This options breaks the XY space lines into very small segments that are then mapped to UV Space. Duplicate UV space points can be removed. Simple but inefficient.

    Pre-processing only.

    AlanX