嗨，这个问题我之前帮同行处理过类似的，咱们一步步拆解来解决：

首先，你说的pyGCodeDecode偏向3D打印机确实是个小局限，但其实咱们完全可以自己针对性做适配，尤其是G2/G3的圆弧插补，核心是把圆弧的运动拆解成时间维度的坐标序列，再推导加速度。

一、G2/G3圆弧插补的核心思路 #

圆弧运动的本质是恒定线速度沿圆弧路径的运动，所以咱们先搞定几个关键参数：

1. 解析圆弧参数：先从G2/G3指令里提取出所有必要信息——
   • 起点坐标（当前机床位置，或者指令前的位置）
   • 终点坐标（指令里的X/Y值）
   • 圆心位置：要么是I/J/K（相对于起点的偏移量，比如I是X方向偏移，J是Y方向），要么是R（圆弧半径，注意有些系统用负R表示大于180度的大圆弧）
2. 计算圆弧的几何属性：
   • 算出圆弧半径R、起始角度θ₀、终止角度θ₁（用反正切函数atan2计算，注意坐标系方向，CNC通常是右手坐标系）
   • 圆弧总长度L = R × Δθ（Δθ是起始到终止的角度差，转成弧度制，注意G2顺时针、G3逆时针的角度变化方向）
3. 时间与坐标的映射：
   • 进给率F是沿路径的线速度，先把单位转成mm/s（F通常是mm/min，除以60即可）
   • 圆弧运动总时间T = L / v（v=F/60）
   • 按你需要的时间步长（比如每0.01秒采样一次），生成时间序列t₀,t₁,...,tn，每个时间点对应的角度θ(t) = θ₀ + (Δθ/T)×t
   • 最终坐标：

     X(t) = Xc + R×cos(θ(t))
     Y(t) = Yc + R×sin(θ(t))
     

     （Xc/Yc是圆心坐标，G2时Δθ是负的，G3是正的，对应顺时针/逆时针方向）

二、加速度的计算 #

加速度分两种情况，要结合CNC的实际运动逻辑：

• 匀速圆弧运动：如果机床在圆弧段全程保持恒定进给率F，那切向加速度为0，但存在向心加速度（法向）：
  aₙ = v²/R（v=F/60，单位mm/s²）
• 带加减速的圆弧运动：真实CNC都会有加减速（比如S曲线加减速），这时候要把圆弧段分成三个阶段：加速段、匀速段、减速段。每个阶段的切向加速度是变化的，向心加速度也会随着实时速度v(t)变化：aₙ(t) = v(t)²/R。你需要结合机床的最大加速度、加减速时间参数来模拟v(t)。

三、实操建议 #

• 自己封装轻量解析器：不用找现成工具，用Python的基础库（math、numpy）就能搞定，重点是处理不同格式的G2/G3指令（比如I/J/K vs R，大圆弧/小圆弧）
• 批量计算与绘图：用numpy生成时间序列和坐标序列，再用matplotlib画x/t、y/t、z/t的曲线，直观又方便
• 螺旋插补处理：如果G2/G3带Z轴指令（比如G2 X10 Y10 Z5 I5 J0 F1000），Z轴是线性变化的，直接按Z的总位移除以总时间，得到Z(t) = Z₀ + (ΔZ/T)×t即可

简单伪代码示例 #

import math
import numpy as np

# 示例G2指令：起点(X0,Y0)=(0,0)，终点(X1,Y1)=(10,0)，圆心偏移I5 J0，F=1000mm/min
start_pos = (0, 0)
end_pos = (10, 0)
center_offset = (5, 0)
F = 1000

# 计算圆心坐标
center_pos = (start_pos[0] + center_offset[0], start_pos[1] + center_offset[1])
R = math.hypot(center_offset[0], center_offset[1])

# 计算起始/终止角度（弧度）
theta_start = math.atan2(start_pos[1] - center_pos[1], start_pos[0] - center_pos[0])
theta_end = math.atan2(end_pos[1] - center_pos[1], end_pos[0] - center_pos[0])

# 处理G2顺时针的角度差（确保是负的）
delta_theta = theta_end - theta_start
if delta_theta > 0:
    delta_theta -= 2 * math.pi

# 计算总长度和总时间
L = abs(delta_theta) * R
v = F / 60  # 转成mm/s
total_time = L / v

# 生成时间点（1000个采样点）
time_points = np.linspace(0, total_time, num=1000)
# 计算每个时间点的角度
theta_points = theta_start + (delta_theta / total_time) * time_points
# 计算X/Y坐标
x_points = center_pos[0] + R * np.cos(theta_points)
y_points = center_pos[1] + R * np.sin(theta_points)

# 计算匀速时的向心加速度
a_n = (v ** 2) / R
print(f"匀速圆弧向心加速度：{a_n} mm/s²")

备注：内容来源于stack exchange，提问作者efish1824