按推荐建议,继续执行

结论:已将 LinuxCNC TRT xyzac/xyzbc 运动学源码纳入 vendored manifest、native probe 和文档验证链,完整 native 验证通过。
This commit is contained in:
2026-06-07 19:21:54 +08:00
parent 1b2f95b3d1
commit e634adfa56
13 changed files with 804 additions and 15 deletions

View File

@@ -0,0 +1,366 @@
/**************************************************************************
* Copyright 2016 Rudy du Preez <rudy@asmsa.co.za>
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
**************************************************************************/
/********************************************************************
* Kinematics functions (forward,inverse) for:
* 1) 5 axis mill (XYZAC)
* This mill has a tilting table (A axis) and horizontal rotary
* mounted to the table (C axis).
* 2) 5 axis mill (XYZBC)
* This mill has a tilting table (B axis) and horizontal rotary
* mounted to the table (C axis).
*
* Note: The directions of the rotational axes are the opposite of the
* conventional axis directions. See
* https://linuxcnc.org/docs/html/gcode/machining-center.html
********************************************************************/
#include <rtapi_math.h>
#include <rtapi_string.h>
#include <rtapi_ctype.h>
#include <hal.h>
#include <emcmotcfg.h>
#include <kinematics.h>
static int trtfuncs_max_joints;
// joint number assignments (-1 ==> not assigned)
static int JX = -1;
static int JY = -1;
static int JZ = -1;
static int JA = -1;
static int JB = -1;
static int JC = -1;
static int JU = -1;
static int JV = -1;
static int JW = -1;
struct haldata {
hal_float_t *x_rot_point;
hal_float_t *y_rot_point;
hal_float_t *z_rot_point;
hal_float_t *x_offset;
hal_float_t *y_offset;
hal_float_t *z_offset;
hal_float_t *tool_offset;
hal_bit_t *conventional_directions; // default: false
} *haldata;
int trtKinematicsSetup(const int comp_id,
const char* coordinates,
kparms* kp)
{
int i,jno,res=0;
int axis_idx_for_jno[EMCMOT_MAX_JOINTS];
int rqdjoints = strlen(kp->required_coordinates);
if (rqdjoints > kp->max_joints) {
rtapi_print_msg(RTAPI_MSG_ERR,
"ERROR %s: supports %d joints, <%s> requires %d\n",
kp->kinsname,
kp->max_joints,
coordinates,
rqdjoints);
goto error;
}
trtfuncs_max_joints = kp->max_joints;
if (map_coordinates_to_jnumbers(coordinates,
kp->max_joints,
kp->allow_duplicates,
axis_idx_for_jno)) {
goto error;
}
// require all chars in reqd_coords (order doesn't matter)
for (i=0; i < rqdjoints; i++) {
char reqd_char;
reqd_char = *(kp->required_coordinates + i);
if ( !strchr(coordinates,toupper(reqd_char))
&& !strchr(coordinates,tolower(reqd_char)) ) {
rtapi_print_msg(RTAPI_MSG_ERR,
"ERROR %s:\nrequired coordinates:%s\n"
"specified coordinates:%s\n",
kp->kinsname, kp->required_coordinates, coordinates);
goto error;
}
}
// assign principal joint numbers (first found in coordinates map)
// duplicates are handled by position_to_mapped_joints()
for (jno=0; jno < EMCMOT_MAX_JOINTS; jno++) {
if (axis_idx_for_jno[jno] == 0 && JX==-1) {JX = jno;}
if (axis_idx_for_jno[jno] == 1 && JY==-1) {JY = jno;}
if (axis_idx_for_jno[jno] == 2 && JZ==-1) {JZ = jno;}
if (axis_idx_for_jno[jno] == 3 && JA==-1) {JA = jno;}
if (axis_idx_for_jno[jno] == 4 && JB==-1) {JB = jno;}
if (axis_idx_for_jno[jno] == 5 && JC==-1) {JC = jno;}
if (axis_idx_for_jno[jno] == 6 && JU==-1) {JU = jno;}
if (axis_idx_for_jno[jno] == 7 && JV==-1) {JV = jno;}
if (axis_idx_for_jno[jno] == 8 && JW==-1) {JW = jno;}
}
rtapi_print("%s coordinates=%s assigns:\n", kp->kinsname,coordinates);
for (jno=0; jno<EMCMOT_MAX_JOINTS; jno++) {
if (axis_idx_for_jno[jno] == -1) break; //fini
rtapi_print(" Joint %d ==> Axis %c\n",
jno,"XYZABCUVW"[axis_idx_for_jno[jno]]);
}
haldata = hal_malloc(sizeof(struct haldata));
if (!haldata) {goto error;}
res += hal_pin_float_newf(HAL_IN, &(haldata->x_rot_point), comp_id,
"%s.x-rot-point",kp->halprefix);
res += hal_pin_float_newf(HAL_IN, &(haldata->y_rot_point), comp_id,
"%s.y-rot-point",kp->halprefix);
res += hal_pin_float_newf(HAL_IN, &(haldata->z_rot_point), comp_id,
"%s.z-rot-point",kp->halprefix);
res += hal_pin_float_newf(HAL_IN, &(haldata->x_offset), comp_id,
"%s.x-offset",kp->halprefix);
res += hal_pin_float_newf(HAL_IN, &(haldata->y_offset), comp_id,
"%s.y-offset",kp->halprefix);
res += hal_pin_float_newf(HAL_IN, &(haldata->z_offset), comp_id,
"%s.z-offset",kp->halprefix);
res += hal_pin_float_newf(HAL_IN, &(haldata->tool_offset), comp_id,
"%s.tool-offset",kp->halprefix);
res += hal_pin_bit_newf(HAL_IN, &(haldata->conventional_directions), comp_id,
"%s.conventional-directions", kp->halprefix);
if (res) {goto error;}
return 0;
error:
rtapi_print_msg(RTAPI_MSG_ERR,"trtKinematicsSetup() FAIL\n");
return -1;
} // trtKinematicsSetup()
int xyzacKinematicsForward(const double *joints,
EmcPose * pos,
const KINEMATICS_FORWARD_FLAGS * fflags,
KINEMATICS_INVERSE_FLAGS * iflags)
{
(void)fflags;
(void)iflags;
const double x_rot_point = *(haldata->x_rot_point);
const double y_rot_point = *(haldata->y_rot_point);
const double z_rot_point = *(haldata->z_rot_point);
const double dt = *(haldata->tool_offset);
const double dy = *(haldata->y_offset);
const double dz = *(haldata->z_offset) + dt;
const double a_rad = joints[JA]*TO_RAD;
const double c_rad = joints[JC]*TO_RAD;
const real_t con = *(haldata->conventional_directions) ? 1.0 : -1.0;
pos->tran.x = + cos(c_rad) * (joints[JX] - x_rot_point)
- con * sin(c_rad) * cos(a_rad) * (joints[JY] - dy - y_rot_point)
+ sin(c_rad) * sin(a_rad) * (joints[JZ] - dz - z_rot_point)
- con * sin(c_rad) * dy
+ x_rot_point;
pos->tran.y = + con * sin(c_rad) * (joints[JX] - x_rot_point)
+ cos(c_rad) * cos(a_rad) * (joints[JY] - dy - y_rot_point)
- con * cos(c_rad) * sin(a_rad) * (joints[JZ] - dz - z_rot_point)
+ cos(c_rad) * dy
+ y_rot_point;
pos->tran.z = + 0
+ con * sin(a_rad) * (joints[JY] - dy - y_rot_point)
+ cos(a_rad) * (joints[JZ] - dz - z_rot_point)
+ dz
+ z_rot_point;
pos->a = joints[JA];
pos->c = joints[JC];
// optional letters (specify with coordinates module parameter)
pos->b = (JB != -1)? joints[JB] : 0;
pos->u = (JU != -1)? joints[JU] : 0;
pos->v = (JV != -1)? joints[JV] : 0;
pos->w = (JW != -1)? joints[JW] : 0;
return 0;
} // xyzacKinematicsForward()
int xyzacKinematicsInverse(const EmcPose * pos,
double *joints,
const KINEMATICS_INVERSE_FLAGS * iflags,
KINEMATICS_FORWARD_FLAGS * fflags)
{
(void)iflags;
(void)fflags;
const double x_rot_point = *(haldata->x_rot_point);
const double y_rot_point = *(haldata->y_rot_point);
const double z_rot_point = *(haldata->z_rot_point);
const double dy = *(haldata->y_offset);
const double dt = *(haldata->tool_offset);
const double dz = *(haldata->z_offset) + dt;
const double a_rad = pos->a*TO_RAD;
const double c_rad = pos->c*TO_RAD;
const real_t con = *(haldata->conventional_directions) ? 1.0 : -1.0;
EmcPose P; // computed position
P.tran.x = + cos(c_rad) * (pos->tran.x - x_rot_point)
+ con * sin(c_rad) * (pos->tran.y - y_rot_point)
+ x_rot_point;
P.tran.y = - con * sin(c_rad) * cos(a_rad) * (pos->tran.x - x_rot_point)
+ cos(c_rad) * cos(a_rad) * (pos->tran.y - y_rot_point)
+ con * sin(a_rad) * (pos->tran.z - z_rot_point)
- cos(a_rad) * dy
- con * sin(a_rad) * dz
+ dy
+ y_rot_point;
P.tran.z = + sin(c_rad) * sin(a_rad) * (pos->tran.x - x_rot_point)
- con * cos(c_rad) * sin(a_rad) * (pos->tran.y - y_rot_point)
+ cos(a_rad) * (pos->tran.z - z_rot_point)
+ con * sin(a_rad) * dy
- cos(a_rad) * dz
+ dz
+ z_rot_point;
P.a = pos->a;
P.c = pos->c;
// optional letters (specify with coordinates module parameter)
P.b = (JB != -1)? pos->b : 0;
P.u = (JU != -1)? pos->u : 0;
P.v = (JV != -1)? pos->v : 0;
P.w = (JW != -1)? pos->w : 0;
// update joints with support for
// multiple-joints per-coordinate letter:
// based on computed position
position_to_mapped_joints(trtfuncs_max_joints,
&P,
joints);
return 0;
} // xyzacKinematicsInverse()
int xyzbcKinematicsForward(const double *joints,
EmcPose * pos,
const KINEMATICS_FORWARD_FLAGS * fflags,
KINEMATICS_INVERSE_FLAGS * iflags)
{
(void)fflags;
(void)iflags;
// Note: 'principal' joints are used
const double x_rot_point = *(haldata->x_rot_point);
const double y_rot_point = *(haldata->y_rot_point);
const double z_rot_point = *(haldata->z_rot_point);
const double dx = *(haldata->x_offset);
const double dt = *(haldata->tool_offset);
const double dz = *(haldata->z_offset) + dt;
const double b_rad = joints[JB]*TO_RAD;
const double c_rad = joints[JC]*TO_RAD;
const real_t con = *(haldata->conventional_directions) ? 1.0 : -1.0;
pos->tran.x = cos(c_rad) * cos(b_rad) * (joints[JX] - dx - x_rot_point)
- con * sin(c_rad) * (joints[JY] - y_rot_point)
+ con * cos(c_rad) * sin(b_rad) * (joints[JZ] - dz - z_rot_point)
+ cos(c_rad) * dx
+ x_rot_point;
pos->tran.y = + con * sin(c_rad) * cos(b_rad) * (joints[JX] - dx - x_rot_point)
+ cos(c_rad) * (joints[JY] - y_rot_point)
+ sin(c_rad) * sin(b_rad) * (joints[JZ] - dz - z_rot_point)
+ con * sin(c_rad) * dx
+ y_rot_point;
pos->tran.z = - con * sin(b_rad) * (joints[JX] - dx - x_rot_point)
+ cos(b_rad) * (joints[JZ] - dz - z_rot_point)
+ dz
+ z_rot_point;
pos->b = joints[JB];
pos->c = joints[JC];
// optional letters (specify with coordinates module parameter)
pos->a = (JA != -1)? joints[JA] : 0;
pos->u = (JU != -1)? joints[JU] : 0;
pos->v = (JV != -1)? joints[JV] : 0;
pos->w = (JW != -1)? joints[JW] : 0;
return 0;
} // xyzbcKinematicsForward()
int xyzbcKinematicsInverse(const EmcPose * pos,
double *joints,
const KINEMATICS_INVERSE_FLAGS * iflags,
KINEMATICS_FORWARD_FLAGS * fflags)
{
(void)iflags;
(void)fflags;
const double x_rot_point = *(haldata->x_rot_point);
const double y_rot_point = *(haldata->y_rot_point);
const double z_rot_point = *(haldata->z_rot_point);
const double dx = *(haldata->x_offset);
const double dt = *(haldata->tool_offset);
const double dz = *(haldata->z_offset) + dt;
const double b_rad = pos->b*TO_RAD;
const double c_rad = pos->c*TO_RAD;
const double dpx = -cos(b_rad)*dx + sin(b_rad)*dz + dx;
const double dpz = -sin(b_rad)*dx - cos(b_rad)*dz + dz;
const real_t con = *(haldata->conventional_directions) ? 1.0 : -1.0;
EmcPose P; // computed position
P.tran.x = + cos(c_rad) * cos(b_rad) * (pos->tran.x - x_rot_point)
+ con * sin(c_rad) * cos(b_rad) * (pos->tran.y - y_rot_point)
- con * sin(b_rad) * (pos->tran.z - z_rot_point)
+ dpx
+ x_rot_point;
P.tran.y = - con * sin(c_rad) * (pos->tran.x - x_rot_point)
+ cos(c_rad) * (pos->tran.y - y_rot_point)
+ y_rot_point;
P.tran.z = + con * cos(c_rad) * sin(b_rad) * (pos->tran.x - x_rot_point)
+ sin(c_rad) * sin(b_rad) * (pos->tran.y - y_rot_point)
+ cos(b_rad) * (pos->tran.z - z_rot_point)
+ dpz
+ z_rot_point;
P.b = pos->b;
P.c = pos->c;
// optional letters (specify with coordinates module parameter)
P.a = (JA != -1)? pos->a : 0;
P.u = (JU != -1)? pos->u : 0;
P.v = (JV != -1)? pos->v : 0;
P.w = (JW != -1)? pos->w : 0;
// update joints with support for
// multiple-joints per-coordinate letter:
// based on computed position
position_to_mapped_joints(trtfuncs_max_joints,
&P,
joints);
return 0;
} // xyzbcKinematicsInverse()

View File

@@ -0,0 +1,57 @@
/********************************************************************
* xyzac-trt-kins.c employing switchkins.[ch]
* License: GPL Version 2
*
* NOTEs:
* 1) specify all kparms items
* 2) specify 3 KS,KF,KI functions for switchkins_type=0,1,2
* 3) the 0th switchkins_type is the startup default
* 4) sparm is a module string parameter for configuration
* 5) The directions of the rotational axes are the opposite of the
* conventional axis directions.
*/
#include <rtapi.h>
#include <rtapi_string.h>
#include <emcmotcfg.h>
#include "switchkins.h"
int switchkinsSetup(kparms* kp,
KS* kset0, KS* kset1, KS* kset2,
KF* kfwd0, KF* kfwd1, KF* kfwd2,
KI* kinv0, KI* kinv1, KI* kinv2
)
{
kp->kinsname = "xyzac-trt-kins"; // !!! must agree with filename
kp->halprefix = "xyzac-trt-kins"; // hal pin names
kp->required_coordinates = "xyzac";
kp->allow_duplicates = 1;
kp->max_joints = EMCMOT_MAX_JOINTS;
if (kp->sparm && strstr(kp->sparm,"identityfirst")) {
rtapi_print("\n!!! switchkins-type 0 is IDENTITY\n");
*kset0 = identityKinematicsSetup;
*kfwd0 = identityKinematicsForward;
*kinv0 = identityKinematicsInverse;
*kset1 = trtKinematicsSetup; // trt: xyzac,xyzbc
*kfwd1 = xyzacKinematicsForward;
*kinv1 = xyzacKinematicsInverse;
} else {
rtapi_print("\n!!! switchkins-type 0 is %s\n",kp->kinsname);
*kset0 = trtKinematicsSetup; // trt: xyzac,xyzbc
*kfwd0 = xyzacKinematicsForward;
*kinv0 = xyzacKinematicsInverse;
*kset1 = identityKinematicsSetup;
*kfwd1 = identityKinematicsForward;
*kinv1 = identityKinematicsInverse;
}
*kset2 = userkKinematicsSetup;
*kfwd2 = userkKinematicsForward;
*kinv2 = userkKinematicsInverse;
return 0;
}

View File

@@ -0,0 +1,57 @@
/********************************************************************
* xyzbc-trt-kins.c employing switchkins.[ch]
* License: GPL Version 2
*
* NOTEs:
* 1) specify all kparms items
* 2) specify 3 KS,KF,KI functions for switchkins_type=0,1,2
* 3) the 0th switchkins_type is the startup default
* 4) sparm is a module string parameter for configuration
* 5) The directions of the rotational axes are the opposite of the
* conventional axis directions.
*/
#include <rtapi.h>
#include <rtapi_string.h>
#include <emcmotcfg.h>
#include "switchkins.h"
int switchkinsSetup(kparms* kp,
KS* kset0, KS* kset1, KS* kset2,
KF* kfwd0, KF* kfwd1, KF* kfwd2,
KI* kinv0, KI* kinv1, KI* kinv2
)
{
kp->kinsname = "xyzbc-trt-kins"; // !!! must agree with filename
kp->halprefix = "xyzbc-trt-kins"; // hal pin names
kp->required_coordinates = "xyzbc";
kp->allow_duplicates = 1;
kp->max_joints = EMCMOT_MAX_JOINTS;
if (kp->sparm && strstr(kp->sparm,"identityfirst")) {
rtapi_print("\n!!! switchkins-type 0 is IDENTITY\n");
*kset0 = identityKinematicsSetup;
*kfwd0 = identityKinematicsForward;
*kinv0 = identityKinematicsInverse;
*kset1 = trtKinematicsSetup; // trt: xyzac,xyzbc
*kfwd1 = xyzbcKinematicsForward;
*kinv1 = xyzbcKinematicsInverse;
} else {
rtapi_print("\n!!! switchkins-type 0 is %s\n",kp->kinsname);
*kset0 = trtKinematicsSetup; // trt: xyzac,xyzbc
*kfwd0 = xyzbcKinematicsForward;
*kinv0 = xyzbcKinematicsInverse;
*kset1 = identityKinematicsSetup;
*kfwd1 = identityKinematicsForward;
*kinv1 = identityKinematicsInverse;
}
*kset2 = userkKinematicsSetup;
*kfwd2 = userkKinematicsForward;
*kinv2 = userkKinematicsInverse;
return 0;
}