2016-07-07 22:52:44 +03:00
|
|
|
# Code for handling the kinematics of cartesian robots
|
2016-05-25 18:37:40 +03:00
|
|
|
#
|
|
|
|
# Copyright (C) 2016 Kevin O'Connor <kevin@koconnor.net>
|
|
|
|
#
|
|
|
|
# This file may be distributed under the terms of the GNU GPLv3 license.
|
2016-07-07 22:52:44 +03:00
|
|
|
import logging
|
|
|
|
import stepper, homing
|
2016-05-25 18:37:40 +03:00
|
|
|
|
2016-07-10 19:23:35 +03:00
|
|
|
StepList = (0, 1, 2)
|
2016-05-25 18:37:40 +03:00
|
|
|
|
2016-07-07 22:20:18 +03:00
|
|
|
class CartKinematics:
|
|
|
|
def __init__(self, printer, config):
|
2016-07-10 19:23:35 +03:00
|
|
|
steppers = ['stepper_x', 'stepper_y', 'stepper_z']
|
2016-07-07 22:20:18 +03:00
|
|
|
self.steppers = [stepper.PrinterStepper(printer, config.getsection(n))
|
|
|
|
for n in steppers]
|
2016-09-30 23:00:32 +03:00
|
|
|
self.limits = [(1.0, -1.0)] * 3
|
2016-07-10 19:23:35 +03:00
|
|
|
self.stepper_pos = [0, 0, 0]
|
2016-07-07 22:20:18 +03:00
|
|
|
def build_config(self):
|
2016-07-22 22:38:13 +03:00
|
|
|
for stepper in self.steppers[:2]:
|
|
|
|
stepper.set_max_jerk(0.005 * stepper.max_accel) # XXX
|
2016-07-07 22:20:18 +03:00
|
|
|
for stepper in self.steppers:
|
|
|
|
stepper.build_config()
|
|
|
|
def set_position(self, newpos):
|
|
|
|
self.stepper_pos = [int(newpos[i]*self.steppers[i].inv_step_dist + 0.5)
|
|
|
|
for i in StepList]
|
2016-09-30 21:47:45 +03:00
|
|
|
def get_max_speed(self):
|
2016-07-07 23:57:25 +03:00
|
|
|
max_xy_speed = min(s.max_velocity for s in self.steppers[:2])
|
|
|
|
max_xy_accel = min(s.max_accel for s in self.steppers[:2])
|
2016-07-07 22:20:18 +03:00
|
|
|
return max_xy_speed, max_xy_accel
|
2016-07-27 05:06:14 +03:00
|
|
|
def get_homed_position(self):
|
2016-10-13 17:04:30 +03:00
|
|
|
return [s.position_endstop + s.get_homed_offset()*s.step_dist
|
|
|
|
for s in self.steppers]
|
2016-07-26 06:47:30 +03:00
|
|
|
def home(self, toolhead, axes):
|
2016-07-07 22:20:18 +03:00
|
|
|
# Each axis is homed independently and in order
|
2016-07-26 06:47:30 +03:00
|
|
|
homing_state = homing.Homing(toolhead, axes)
|
|
|
|
for axis in axes:
|
|
|
|
s = self.steppers[axis]
|
2016-09-30 23:00:32 +03:00
|
|
|
self.limits[axis] = (s.position_min, s.position_max)
|
2016-07-26 06:47:30 +03:00
|
|
|
# Determine moves
|
|
|
|
if s.homing_positive_dir:
|
2016-09-15 19:20:16 +03:00
|
|
|
pos = s.position_endstop - 1.5*(
|
|
|
|
s.position_endstop - s.position_min)
|
2016-07-26 06:47:30 +03:00
|
|
|
rpos = s.position_endstop - s.homing_retract_dist
|
|
|
|
r2pos = rpos - s.homing_retract_dist
|
|
|
|
else:
|
|
|
|
pos = s.position_endstop + 1.5*(
|
|
|
|
s.position_max - s.position_endstop)
|
|
|
|
rpos = s.position_endstop + s.homing_retract_dist
|
|
|
|
r2pos = rpos + s.homing_retract_dist
|
|
|
|
# Initial homing
|
|
|
|
homepos = [None, None, None, None]
|
|
|
|
homepos[axis] = s.position_endstop
|
|
|
|
coord = [None, None, None, None]
|
|
|
|
coord[axis] = pos
|
|
|
|
homing_state.plan_home(list(coord), homepos, [s], s.homing_speed)
|
|
|
|
# Retract
|
|
|
|
coord[axis] = rpos
|
|
|
|
homing_state.plan_move(list(coord), s.homing_speed)
|
|
|
|
# Home again
|
|
|
|
coord[axis] = r2pos
|
|
|
|
homing_state.plan_home(list(coord), homepos, [s], s.homing_speed/2.0)
|
2016-07-07 22:20:18 +03:00
|
|
|
return homing_state
|
|
|
|
def motor_off(self, move_time):
|
2016-09-30 23:00:32 +03:00
|
|
|
self.limits = [(1.0, -1.0)] * 3
|
2016-07-07 22:20:18 +03:00
|
|
|
for stepper in self.steppers:
|
|
|
|
stepper.motor_enable(move_time, 0)
|
2016-09-22 18:09:20 +03:00
|
|
|
def query_endstops(self, move_time):
|
|
|
|
return homing.QueryEndstops(["x", "y", "z"], self.steppers)
|
2016-09-30 23:00:32 +03:00
|
|
|
def check_endstops(self, move):
|
|
|
|
for i in StepList:
|
|
|
|
if (move.axes_d[i]
|
|
|
|
and (move.pos[i] < self.limits[i][0]
|
|
|
|
or move.pos[i] > self.limits[i][1])):
|
|
|
|
if self.limits[i][0] > self.limits[i][1]:
|
|
|
|
raise homing.EndstopError(move.pos, "Must home axis first")
|
|
|
|
raise homing.EndstopError(move.pos)
|
2016-09-30 21:47:45 +03:00
|
|
|
def check_move(self, move):
|
2016-09-30 23:00:32 +03:00
|
|
|
limits = self.limits
|
|
|
|
xpos, ypos = move.pos[:2]
|
|
|
|
if (xpos < limits[0][0] or xpos > limits[0][1]
|
|
|
|
or ypos < limits[1][0] or ypos > limits[1][1]):
|
|
|
|
self.check_endstops(move)
|
2016-09-30 21:47:45 +03:00
|
|
|
if not move.axes_d[2]:
|
|
|
|
# Normal XY move - use defaults
|
|
|
|
return
|
|
|
|
# Move with Z - update velocity and accel for slower Z axis
|
2016-09-30 23:00:32 +03:00
|
|
|
self.check_endstops(move)
|
2016-09-30 21:47:45 +03:00
|
|
|
axes_d = move.axes_d
|
|
|
|
move_d = move.move_d
|
|
|
|
velocity_factor = min([self.steppers[i].max_velocity / abs(axes_d[i])
|
|
|
|
for i in StepList if axes_d[i]])
|
|
|
|
accel_factor = min([self.steppers[i].max_accel / abs(axes_d[i])
|
|
|
|
for i in StepList if axes_d[i]])
|
|
|
|
move.limit_speed(velocity_factor * move_d, accel_factor * move_d)
|
2016-07-07 22:20:18 +03:00
|
|
|
def move(self, move_time, move):
|
|
|
|
inv_accel = 1. / move.accel
|
|
|
|
inv_cruise_v = 1. / move.cruise_v
|
|
|
|
for i in StepList:
|
|
|
|
new_step_pos = int(move.pos[i]*self.steppers[i].inv_step_dist + 0.5)
|
|
|
|
steps = new_step_pos - self.stepper_pos[i]
|
|
|
|
if not steps:
|
|
|
|
continue
|
|
|
|
self.stepper_pos[i] = new_step_pos
|
|
|
|
sdir = 0
|
|
|
|
if steps < 0:
|
|
|
|
sdir = 1
|
|
|
|
steps = -steps
|
2016-08-24 23:42:25 +03:00
|
|
|
mcu_time, so = self.steppers[i].prep_move(move_time, sdir)
|
2016-07-07 22:20:18 +03:00
|
|
|
|
|
|
|
step_dist = move.move_d / steps
|
|
|
|
step_offset = 0.5
|
|
|
|
|
|
|
|
# Acceleration steps
|
|
|
|
#t = sqrt(2*pos/accel + (start_v/accel)**2) - start_v/accel
|
2016-08-24 23:42:25 +03:00
|
|
|
accel_time_offset = move.start_v * inv_accel
|
|
|
|
accel_sqrt_offset = accel_time_offset**2
|
|
|
|
accel_multiplier = 2.0 * step_dist * inv_accel
|
2016-07-07 22:20:18 +03:00
|
|
|
accel_steps = move.accel_r * steps
|
|
|
|
step_offset = so.step_sqrt(
|
2016-08-24 23:42:25 +03:00
|
|
|
mcu_time - accel_time_offset, accel_steps, step_offset
|
2016-07-07 22:20:18 +03:00
|
|
|
, accel_sqrt_offset, accel_multiplier)
|
2016-08-24 23:42:25 +03:00
|
|
|
mcu_time += move.accel_t
|
2016-07-07 22:20:18 +03:00
|
|
|
# Cruising steps
|
|
|
|
#t = pos/cruise_v
|
2016-08-24 23:42:25 +03:00
|
|
|
cruise_multiplier = step_dist * inv_cruise_v
|
2016-07-07 22:20:18 +03:00
|
|
|
cruise_steps = move.cruise_r * steps
|
|
|
|
step_offset = so.step_factor(
|
2016-08-24 23:42:25 +03:00
|
|
|
mcu_time, cruise_steps, step_offset, cruise_multiplier)
|
|
|
|
mcu_time += move.cruise_t
|
2016-07-07 22:20:18 +03:00
|
|
|
# Deceleration steps
|
|
|
|
#t = cruise_v/accel - sqrt((cruise_v/accel)**2 - 2*pos/accel)
|
2016-08-24 23:42:25 +03:00
|
|
|
decel_time_offset = move.cruise_v * inv_accel
|
|
|
|
decel_sqrt_offset = decel_time_offset**2
|
2016-07-07 22:20:18 +03:00
|
|
|
decel_steps = move.decel_r * steps
|
|
|
|
so.step_sqrt(
|
2016-08-24 23:42:25 +03:00
|
|
|
mcu_time + decel_time_offset, decel_steps, step_offset
|
2016-07-07 22:20:18 +03:00
|
|
|
, decel_sqrt_offset, -accel_multiplier)
|