AXY:Element Coordination/One Dimensional Tracking Algorithm

One Dimensional Tracking Algorithm

This algorithm is used always in simple axis based MF tracking:

The algorithm consist of a state machine with the following states:

• Tracking process

State: Full Synchronization follow the master

• track ← master

• Check if master velocity or acceleration exceeds max master frame's values if yes return an error

• flopstime = flops = 0;

State: Tracking process

Initial Setup and testing

Scaled value of AccMaxTran (or AccMaxRot) is used of checking of acceleration limit.

Scaled value of VelMaxTran (or VelMaxRot) is used of checking of velocity limit.

• Check if master velocity or acceleration exceeds max master frame's values if yes return an error
  

• Compute differences:
  

$pos_{diff}=master_{pos}-track_{pos}$ pure delta-position
$vel_{diff}=master_{vel}-track_{vel}$ Pure delta-velocity

Time-to-Reach initial estimation

${sync_{jrk}}$ defined as a scaled value SyncJerkRot or SyncJerkRot

${sync_{acc}}$ defined as a scaled value SyncAcclerationRot or SyncAcclerationRot

This is an estimation of time needed to reach the target. It is not an exact formula
if the velocity difference is positive, it means the robot velocity is lower then conveyors
therefore we need more time once to get to a higher velocity then the conveyers and once
to get from that (higher) velocity to conveyers - therefore the multiplication by 2

$time_{toreach}={\sqrt {\frac {K*|{vel_{diff}}|}{sync_{jrk}}}}+{\sqrt {\frac {4*|{pos_{diff}}|}{sync_{acc}}}}$

where K is 4 if ($vel_{diff}>0$) else it is 2.

DampingFactor - user defined by: MF.DAMPINGFACTOR

Additionaly we mutliply by a prediction factor - to decrease acc/jerk
(rounding to integer number of samples):

time_to_reach *= DampingFactor
time_to_reach = T*int(time_to_reach/T+1) 

Prediction of the Rendezvous-Point

The predicted position of the rendezvous (minus T, if it will happen in this sample!)
pred_pos: position prediction assuming constant acceleration
pred_diff:Predicted position difference
pred_vel: velocity prediction

t = (time_to_reach-T);
pred_pos = master_pos + t*(master_vel + t*master_acc/2);                    
pred_diff = pred_pos - track_pos;		
pred_vel = master_vel + t*master_acc;	

acceleration will be zero

pred_acc = 0;																					  

Polynomial coefficients of 5-degree polynom connecting the current p,v,a to the randevous p,v,a
$a0=track_{pos}$ - not used, directly written into code
$a1=track_{vel}$ - not used, directly written into code
$a2=0.5*track_{acc}$ - not used, directly written into code
$a3=(0.5*(pred_{acc}-3*track_{acc})~+~(10*pred_{diff}/time_{toreach}$$-~6*track_{vel}~-~4*pred_{vel})/time_{toreach})/time_{toreach}$
$a4=(((5*pred_{diff}/time_{toreach}~-~(pred_{vel}+4*track_{vel}))/time_{toreach}$$-1.5*track_{acc})/time_{toreach}~-~2*a3)/time_{toreach}$
$a5=((((pred_{diff}/time_{toreach}~-~track_{vel})/time_{toreach}$$-0.5*track_{acc})/time_{toreach}~-~a3)/time_{toreach}~-~a4)/time_{toreach}$

Next SamplePosition

The position at the next sample will be:
dp = ((((a5*T + a4)*T + a3)*T + 0.5*track_acc)*T + track_vel)*T - Initial increment

The accleration and jerk increments could be computed in this way also: ....

ddp    =  (((5*4*a5*T + 4*3*a4)*T + 3*2*a3)*T + track_acc)*T*T; 			
dddp   =  ((5*4*3*a5*T + 4*3*2*a4)*T + 3*2*a3)*T*T*T*T; 				

But they are computed as average values for the time of sample duration, this approach doeas fit more to
the tools we are using to test this feature, (BMDS...) as they all compute accleration & jerk by numeric differentiation
inside one sample.
The side effect of all this is that we can have a change in jerk/acc sign inside one sample and the average value we are using
ddp,dddp could have an opposite sign, so we are seeing effects like jerk or accleeration is limited from a "wrong side" i.e. instead
using jmax -jmax is used. There is not much we can do against it, using exact jerk/acc values (from polynom) returns than wrong avergae values
we observe as jerk/acc exceed over theri max values. So the avergae value approach is not perefect but gives results within the given limits.

Initial setting of predicted time and the filter limit flag ( if the filter did work - this flag will be on)

T1           = T 
limit        = false 
state.satVel = false 
state.satAcc = false  
state.satJrk = false 

Velocity, Acceleration and Jerk Filter

sync_vel defined as a scaled value VelSyncTran or VelSyncRot

sync_acc defined as a scaled value AccSyncTran or AccSyncRot

sync_jrk defined as a scaled value JrkSyncTran or JrkSyncRot

T is sampling time

if (dp > sync_velT)

dp =  sync_velT; 
limit = state.satVel = true; 

else if (dp < -sync_velT)

dp = -sync_velT; 
limit = state.satVel = true;

Setting Initial velocity delta

ddp = dp - dp0

if (ddp > sync_accT2)

	
dp  =  sync_accT2 + dp0;	
ddp =  sync_accT2;  
limit = state.satAcc = true;

else if (ddp < -sync_accT2)

	
dp  = -sync_accT2 + dp0;	
ddp = -sync_accT2;  
limit = state.satAcc = true;

Setting Initial acceleration delta

dddp = ddp - ddp0

if (dddp > sync_jrkT3 )

 
ddp =  sync_jrkT3 + ddp0; 
dp = ddp + dp0;  
limit = state.satJrk = true;

else if (dddp < -sync_jrkT3 )

 
ddp = -sync_jrkT3 + ddp0; 
dp = ddp + dp0;  
limit = state.satJrk = true;

New position is:

track_pos += dp;

If the filter was activated then:

if the velocity,accleration or jerk is limited then the velocity computed using polynomials does not fit any more, this is observable as PCMD/VCDM missmatch (Issue 2596) if(limit)

track_vel   = dp/T;
track_acc   = ddp/T/T;

else Velocity & Accleration will be computed exactly not using dp and ddp, it has been shown that this approach is much better - means more stable

track_vel   = (((5*a5*T1  + 4*a4)*T1  + 3*a3)*T1  + track_acc)*T1 + track_vel;
track_acc   = ((20*a5*T1 + 12*a4)*T1  + 6*a3)*T1   + track_acc;

keep the old values for next sample

ddp0 = ddp; dp0 = dp;

System is synchronized if NEW delta's are under certesian thresholds:

dv0 = dv; - keep the previous dv
dv = master_vel - track_vel

Successful tracking completion

if((fabs(dv) <= syncGain * sync_accT) && (fabs(master_pos - track_pos) <= syncGain * sync_accT2)) Synchronization
Take the middle for an additional smoothing:

track_pos	   =	0.5*(track_pos + master_pos);
track_vel	   =	0.5*(track_vel + master_vel);
ResetFilter();		// Reset the filter

state.tracking = MOT_MASTER_SYNC; // System synchronized else if( dv0*dv < 0) Check against over-oscillation

if (flops == -1) flops = 0; // ignore the first sample
else	         flops++;
flopstime = 0;

Over-Oscillation Check

maxFlops defined by the user with: MF.MaxFlops

if (flops > maxFlops && (flopstime > maxFlopsTime && fabs(dv) > syncGain * sync_accT2))

state.tracking 	= MOT_MASTER_OSC; // OScilationss

else

flopstime++;

State: Stopping (de-synchronization process)

flopstime = flops = 0;

• De-Syncing profile, follow the given profile, ignore the real source

desync.Profile();
track_pos = track_pos_0+master_direction*desync.path.curr_pos;   	
track_vel = master_direction*desync.path.vel;			
track_acc = master_direction*desync.path.acc;			
if(desync.path.Status.stop)			// if stopped,change the state to OFF
	 SetState(MOT_MASTER_OFF);