[编程申请]：已知两点和其间的弧长画圆弧

caibaobao

*-*3 *-*3 *-*3
请大侠帮忙编程：已知两点座标以及其间的弧长，画圆弧——在画地形图中红线范围时常用到！！！

yshf

;请试用以下程序
;并注意已知两点输入的先后顺序不同会不同的结果

2. 
   
3. (defun c:lc()
   
4.     (setq cm (getvar "cmdecho") os (getvar "osmode"))
   
5.     (setvar "cmdecho" 0)(setvar "osmode" 33)
   
6.     (setq p1 (getpoint "\n请输入第一点:")
   
7.           p2 (getpoint p1 "\n请输入第二点:")
   
8.           l  (getdist "\n请输入圆弧长度L=")
   
9.           c  (distance p1 p2)
   
10.           r0 1e45      dr 1        b3 (/ 1.0 3.0)
    
11.           
    
12.    )
    
13.    (while (> dr 1e-6)
    
14.       (setq b  (/  l 2.0 r0)    bb (* b b)
    
15.             r1 (/ c 2.0 (* b (- 1 (* (/ bb 2.0) (- b3 (* (/ bb 20.0)(- b3 (/ bb 131.0))))))))
    
16.       )
    
17.       (setq dr (abs (- r1 r0))  r0 r1)
    
18.    )
    
19.    (setvar "osmode" 0)
    
20.    (command "_arc" p1 "e" p2 "r" r1 "chprop" (entlast) "" "c" 1 "")
    
21.    (setvar "cmdecho" cm)(setvar "osmode" os)(princ)
    
22. )

<br /> (defun c:lc()<br />     (setq cm (getvar "cmdecho") os (getvar "osmode"))<br />     (setvar "cmdecho" 0)(setvar "osmode" 33)<br />     (setq p1 (getpoint "\n请输入第一点:")<br />           p2 (getpoint p1 "\n请输入第二点:")<br />           l  (getdist "\n请输入圆弧长度L=")<br />           c  (distance p1 p2)<br />           r0 1e45      dr 1        b3 (/ 1.0 3.0)<br />           <br />    )<br />    (while (> dr 1e-6)<br />       (setq b  (/  l 2.0 r0)    bb (* b b)<br />             r1 (/ c 2.0 (* b (- 1 (* (/ bb 2.0) (- b3 (* (/ bb 20.0)(- b3 (/ bb 131.0))))))))<br />       )<br />       (setq dr (abs (- r1 r0))  r0 r1)<br />    )<br />    (setvar "osmode" 0)<br />    (command "_arc" p1 "e" p2 "r" r1 "chprop" (entlast) "" "c" 1 "")<br />    (setvar "cmdecho" cm)(setvar "osmode" os)(princ)<br /> )

xyp1964

楼上的结果不准确！

yshf

请问xyp1964，是如何的不准确?

xyp1964

测试结果：
请输入第一点:
请输入第二点:
请输入圆弧长度L=1200
Command: li nil

Select objects: 1 found

Select objects:
                  ARC       Layer: "0"
                            Space: Model space
                   Color: 1 (red)    Linetype: "BYLAYER"
                   Handle = 16A
            center point, X=2823.1930  Y=-1787.4167  Z=   0.0000
            radius  420.9321
             start angle    188
               end angle    351
            length 1200.1443

caibaobao

大侠们加油啊，小弟在这里给大家加油啦！！！

eachy

最初由 caibaobao 发布
[B]大侠们加油啊，小弟在这里给大家加油啦！！！ [/B]

你经常怎么算出来或者画出来？

民工

这里有个程序可以画

xyp1964

楼上的程序似乎是用逼近法得出的结果。

2. [FONT=courier new](while (> (abs (- x xx)) 0.0000000001)
   
3.         (setq x          xx
   
4.               fx  (- (/ (sin (/ x 2)) x) (/ (* 0.5 l) s))
   
5.               flx (/ (- (* 0.5 x (cos (* 0.5 x))) (sin (* 0.5 x))) (* x x))
   
6.               xx  (- x (/ fx flx))
   
7.         )
   
8.       )[/FONT]

[FONT=courier new](while (> (abs (- x xx)) 0.0000000001)<br />         (setq x          xx<br />               fx  (- (/ (sin (/ x 2)) x) (/ (* 0.5 l) s))<br />               flx (/ (- (* 0.5 x (cos (* 0.5 x))) (sin (* 0.5 x))) (* x x))<br />               xx  (- x (/ fx flx))<br />         )<br />       )[/FONT]

同心

不能下载

xyp1964

2. [FONT=courier new](load "xyp_lib")
   
3. ;|加载通用函数(可在签名栏直接下载)
   
4. 如果已经下载xyp_lib并放到搜索路径下可以不再下载!
   
5. 利用以下任何一种方式即可加载和运行通用函数内的所有子程序：
   
6. 1.在acad.lsp中增加(load"xyp_lib")
   
7. 2.在每个程序内增加(load"xyp_lib")
   
8. 3.在command下，输入(load"xyp_lib")
   
9. 4.在菜单.mnl中增加(load"xyp_lib")
   
10. 5.将xyp_lib.vlx文件直接拽到cad屏幕
    
11. [COLOR=red] ★通用函数下载地址：[/COLOR]
    
12. [url]http://www.xdcad.net/forum/attachment.php?s=&postid=1606661[/url]
    
13. [url]http://www.mjtd.com/bbs/dispbbs.asp?boardID=3&ID=37554&page=1[/url]|;
    
14. 
    
15. ;;;起点、终点和弧长画圆弧
    
16. (defun c:test212 ()
    
17.   (cmdla0)
    
18.   (setq        pt1 (getpoint "\圆弧端点 : ")
    
19.         pt2 (getpoint pt1 "\圆弧终点 : ")
    
20.         no1 (UREAL 7 "" "\n弧长" no1)
    
21.         a   (/ (distance pt1 pt2) 2.0)
    
22.         l   (/ no1 2.0)
    
23.         b   0.001
    
24.   )
    
25.   (if (> l b)
    
26.     (progn
    
27.       (while (< (/ b (sin b)) (/ l a 1.0))
    
28.         (setq b (+ b 0.000001))
    
29.       )
    
30.       (setq ang (rad2ang (* b 2)))
    
31.     )
    
32.   )
    
33.   (mkla"圆弧"1)
    
34.   (setvar "osmode"0)
    
35.   (command "arc" pt1 "e" pt2 "a" ang)
    
36.   (princ "弧长误差 ＝ ")
    
37.   (princ (rtos (- no1 (CURVE-LENG (entlast))) 2 8));测试线长
    
38.   (cmdla1)
    
39. )[/FONT]

[FONT=courier new](load "xyp_lib")<br /> ;|加载通用函数(可在签名栏直接下载)<br /> 如果已经下载xyp_lib并放到搜索路径下可以不再下载!<br /> 利用以下任何一种方式即可加载和运行通用函数内的所有子程序：<br /> 1.在acad.lsp中增加(load"xyp_lib")<br /> 2.在每个程序内增加(load"xyp_lib")<br /> 3.在command下，输入(load"xyp_lib")<br /> 4.在菜单.mnl中增加(load"xyp_lib")<br /> 5.将xyp_lib.vlx文件直接拽到cad屏幕<br /> [COLOR=red] ★通用函数下载地址：[/COLOR]<br /> [url]http://www.xdcad.net/forum/attachment.php?s=&postid=1606661[/url]<br /> [url]http://www.mjtd.com/bbs/dispbbs.asp?boardID=3&ID=37554&page=1[/url]|;<br /> <br /> ;;;起点、终点和弧长画圆弧<br /> (defun c:test212 ()<br />   (cmdla0)<br />   (setq        pt1 (getpoint "\圆弧端点 : ")<br />         pt2 (getpoint pt1 "\圆弧终点 : ")<br />         no1 (UREAL 7 "" "\n弧长" no1)<br />         a   (/ (distance pt1 pt2) 2.0)<br />         l   (/ no1 2.0)<br />         b   0.001<br />   )<br />   (if (> l b)<br />     (progn<br />       (while (< (/ b (sin b)) (/ l a 1.0))<br />         (setq b (+ b 0.000001))<br />       )<br />       (setq ang (rad2ang (* b 2)))<br />     )<br />   )<br />   (mkla"圆弧"1)<br />   (setvar "osmode"0)<br />   (command "arc" pt1 "e" pt2 "a" ang)<br />   (princ "弧长误差 ＝ ")<br />   (princ (rtos (- no1 (CURVE-LENG (entlast))) 2 8));测试线长<br />   (cmdla1)<br /> )[/FONT]

caibaobao

请教楼上斑竹，大概要用到您的函数库吧，那个XYP_LIB？

梦断江南

8楼的不错,精度和速度都比较满意.
是哈哈写的还是长不大版主写的?