[求助]：如何解决这个程序的运行时间问题

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))
前几次圆弧插补的时候，速度还可以，但到后面的时候，每走一步就要十秒左右，整个程序运行完，起码要四五天。
大家有没有办法解决这个运行时间问题

wujimmy

我建议你在程序运行的过程中设置记录时间的函数
然后分析一下是哪些函数用的时间最长,找出原因后问题就好解决了
你的函数说老实话看不懂,三维的没怎么用:)
记录时间的函数可参考我的一篇文章:ssget函数....

xuweiwei417

(command "subtract" e1 "" e2 "")
这一步花费的时间最长了

wujimmy

既然用人家现成的,那没有办法,
速度不好岣?
只有想想用别的函数或者别的算法
干脆用二维看会不会快一点

xuweiwei417

我模拟的就是要和加工出来的一模一样，我的最终目的是要模拟出来的工件，在yz平面上的截面，分析这一截面的形状。
因为砂轮和工件接触不是固定点，所以只能用三维进行模拟

zxq0220

提高速度的另一个办法是：
1、将砂轮画好后放在固定的位置，用拷贝后的砂轮进行差运算。
2、对坐标系进行变化，取消ROTATE3D命令。
3、关闭UNDO。

dandee

用布尔运算当然是花时间的

aeo

(command "rotate3d" e2 "" "x" cen 90)
应该为 (command "rotate3d" e2 "" "x" cen "90")

(command "subtract" e1 "" e2 "")后 用(gc)一下不知会不会改善
我的电脑上速度还可以(1G内存)

因为每次都裁掉一点,越到后面,实体越复杂,所以会很慢

应该考虑,作一份,再拷贝(一份全部制作完整,包括插补),再一起合并就比较快了.

xuweiwei417

楼上的意思是不是把整个加工过程分成多次，再用并集进行呢
但如果这样的话，包络迹高度就可能发生变化了呀