[求助]：帮我看一下这个程序，为什么运行到最后速度会越来越慢

xuweiwei417

原程序如下：
;;;程序目的为对数控加工椭圆进行模拟
;;;变量说明：a,b,c为椭圆轴，l为工件加工的椭圆长度，h包络迹，r为砂轮半径，rr为砂轮小圆弧半径，y为砂轮起始点的y轴坐标
;;; d为相邻两次插补时砂轮轴心的距离，（x0,y0）为切点坐标，（x00,y00)为砂轮轴心坐标
(d defun c:jiagong()
  (setq orig_osnp(getvar "osmode"))
  (setq orig_cmd(getvar "cmdecho"))
  (setvar "osmode" 0)
  (setvar "cmdecho" 0)
  (setq a 170.0)
  (setq b 68.0)
  (setq c 68.0)
  (setq l 100.0)
  (setq h 0.001)
  (setq r 87.08)
  (setq rr 14.29)
  (setq y -30)
  (setq m 0)
  (setq d(* 2 (sqrt (- (* r r) (* (- r h) (- r h))))))
  (setq x0(- 0 (/ l 2)))                                      ;给定起始切点坐标
  (setq y0(sqrt ( * (- 1 (/ (* x0 x0) (* a a))) (* b b))))
  (setq g(- b (* c (sqrt(- 1 (/ (* (/ l 2) (/ l 2)) (* a a)))))))   ;椭圆高度计算
  (setq k1(- 0 (/ (* b b x0) (* a a y0))))                          ;切点的斜率
  (setq k2(- 0 (/ 1 k1)))                                            
  (setq hyy0(- y0 (/ (* r k2) (sqrt(+ 1 (* k2 k2))))))              ;起始点的砂轮轴心的y轴坐标
  (setq hyy0 (abs hyy0))
  (setq gy(- (+ hyy0 g) (+ r 15)))                                  ;计算工件的底面的z轴坐标
  (setq p1(list -50 -30 gy))                                        ;工件底面的两个对角点的坐标
  (setq p2(list 50  30 gy))
  (command "box" p1 p2 16 "")                                       ;画出工件
  (setq e1 (entlast))
  (while (<= x0 (/ l 2))
    (setq y0(sqrt ( * (- 1 (/ (* x0 x0) (* a a))) (* b b))))
    (setq k1(- 0 (/ (* b b x0) (* a a y0))))
  (setq k2(- 0 (/ 1 k1)))
  (if(<= x0 0)
      (progn
        (setq xx0(- x0 (/ r (sqrt(+ 1 (* k2 k2))))))
        (setq yy0(- y0 (/ (* r k2) (sqrt(+ 1 (* k2 k2))))))
        )
      )
  (if(> x0 0)
      (progn
        (setq xx0(+ x0 (/ r (sqrt(+ 1 (* k2 k2))))))
        (setq yy0(+ y0 (/ (* r k2) (sqrt(+ 1 (* k2 k2))))))
        (princ yy0)
        )
      )
(if(= x0 (- 0 (/ l 2)))
      (progn       
       (setq hyy0(- (abs yy0) r))                           ;计算开始时的插补半径，这个半径就为椭圆球的轴心线到工件平面的高度
       (setq ys(+ hyy0 rr))
       )
      )
(setq x1(+ x0 (/ d (sqrt(+ 1 (* k1 k1))))))                 ;x0增加
(setq r1(- (abs yy0) r))
(setq r2(+ r1 rr))                                          ;圆弧插补半径
(setq w(sqrt(- (* r2 r2) (* ys ys))))                       ;圆弧宽度
(setq z1(sqrt(- (* yy0 yy0) (* w w))))                      ;圆弧端点时的砂轮轴心z轴坐标
  (setq cen(list xx0 -30 z1))                                ;在工件端面时的砂轮轴心坐标
  (setq r5(- r rr))                                          ;圆弧半径
  (while (<= y 30)
    (if(or(<= y (- 0 w)) (>= y w))
      (progn
        (setq cen1 cen)
        (setq y(+ y 0.2))                                    ;圆弧插补步长
        (setq cen(list xx0 y z1))
        (command "torus" cen r5 rr)                          ;生成砂轮
        (setq e2 (entlast))
        (command "rotate3d" e2 "" "x" cen 90)                ;把砂轮旋转到xz平面
        (setq e2 (entlast))
        (command "subtract" e1 "" e2 "")                     ;切除工件
        (setq e1(entlast))
        
      )
      )
    (if(and(> y (- 0 w)) (< y w))
      (progn
        (setq cen1 cen)
        (setq y(+ y 0.2))
        (setq z1(sqrt(- (* yy0 yy0) (* y y))))             ;计算圆弧插补时，y变时，z的坐标
        (setq cen(list xx0 y z1))
        (command "torus" cen r5 rr)
        (setq e2 (entlast))
        (command "rotate3d" e2 "" "x" cen 90)
        (setq e2 (entlast))
        (command "subtract" e1 "" e2 "")
        (command e1(entlast))
      
        )
      ))
    (setq x0 x1)
    (setq y -30)                                        ;砂轮回到工件端面
    (setq m(+ m 1))
    )
  (setvar "osmode" orig_osnp)
  (setvar "cmdecho" orig_cmd))
前几次圆弧插补的时候，速度还可以，但到后面的时候，每走一步就要十秒左右，整个程序运行完，起码要四五天。
大家有没有办法解决这个运行时间问题

zxq0220

程序反复用实体的差运算，当然越来越慢。建议解决的办法：可以将砂轮形状改为圆柱形，不要沿Y轴加工，直接进行Z-X平面加工。

xuweiwei417

改成圆柱就不可能与实际加工相符了。