[编程申请]：条件有变化：如果X、Y轴比例不一致？由一点引出两直线a、b,如何求两线a、

Andyamy

由一点引出两直线a、b,如何求两线a、b的夹角A及R、L、E、T这五要素？R由用户输入。
在X、Y轴比例一致的情况下即都是1：1的条件下求得的A角度即为实际角度，如果X轴比例为1：2000，Y轴比例为1：100，那么该如何计算A角的实际角度？
谁能提供算法？

cy956

看题目还以为是搞公路的兄弟。
你这个问题约束不够，数学上说是无穷解。

Andyamy

即使是无穷解，也应该有一种算法。您能提供吗？

leochen

是否这样

Andyamy

OK，谢谢。基本上应该是这样。我试试。

leochen

不好意思，写错了，求角B应该
(setq B (abs (- (angle p1 p3) (angle p1 p2))))

coolzhb

最初由 Andyamy 发布
[B]由一点引出两直线a、b,如何求两线a、b的夹角A及R、L、E、T这五要素？R由用户输入。

谁能提供算法？ [/B]

我来试试！

2. 
   
3.   [FONT=courier new]
   
4. (defun dist (p1 p2 / dis)
   
5.   (setq dis (distance p1 p2))
   
6. )
   
7. 
   
8. 
   
9. (defun zhberr (msg)
   
10.   (setvar "osmode" os)
    
11.   (setvar "cmdecho" 1)
    
12.   (setvar "clayer" old_lay)
    
13.   (setvar "dimzin" dimzin_old)
    
14.   (setq *error* olderr)
    
15.   (princ)
    
16. )
    
17. 
    
18. (defun c:yh (/ pt1 pt2 p11 p12 p21 p22 scale e1 e2 ne1 ne2 text text1)
    
19.   (command "color" "")
    
20.   (setq os (getvar "osmode")
    
21.         olderr *error*
    
22.         *error* zhberr
    
23.   )
    
24.   (princ "\n*** 路线圆弧连接V1.0，上海宝钢工程技术有限公司 周海兵编 ***")
    
25.   (setvar "cmdecho" 0)
    
26.   (setq old_lay (getvar "clayer"))
    
27.   (setq dimzin_old (getvar "Dimzin"))
    
28.   (setvar "dimzin" 8)
    
29.   (if (= (tblobjname "LAYER" "T1") nil)
    
30.     (entmake (list '(0 . "LAYER") '(100 . "AcDbSymbolTableRecord") '
    
31.                    (100 . "AcDbLayerTableRecord") '(6 . "CONTINUOUS") '
    
32.                    (62 . 3) '(70 . 0) (cons 2 "T1")
    
33.              )
    
34.     )
    
35.   )
    
36.   (if (= (tblobjname "LAYER" "T2") nil)
    
37.     (entmake (list '(0 . "LAYER") '(100 . "AcDbSymbolTableRecord") '
    
38.                    (100 . "AcDbLayerTableRecord") '(6 . "CONTINUOUS") '
    
39.                    (62 . 5) '(70 . 0) (cons 2 "T2")
    
40.              )
    
41.     )
    
42.   )
    
43.   (if (= nil scaleqs) (setq scaleqs 1000))
    
44.         (setq scale (getreal (strcat "\n *** 初始化比例 <" (rtos scaleqs 2 0) "> *** :")))
    
45.         (if (= scale nil)
    
46.                 (setq scale scaleqs)
    
47.                 (setq scaleqs scale)
    
48.         )
    
49.   (setvar "osmode" 512)
    
50.   (if (= nil r0)
    
51.     (setq r0 (getvar "filletrad"))
    
52.   )
    
53.   (setq pt1 (getpoint "\n选择第一条线: "))
    
54.   (setq pt2 (getpoint "\n选择第二条线: "))
    
55.   (setvar "osmode" 0)
    
56.   (prompt (strcat "\n输入圆角半径<" (rtos r0 2 0) "> m :"))
    
57.   (initget 4)
    
58.   (setq r (getreal))
    
59.   (if (= nil r)
    
60.     (setq r r0)
    
61.     (setq r0 r)
    
62.   )
    
63.   (setvar "filletrad" r)
    
64.   (setq e1 (ssget pt1))
    
65.   (setq e2 (ssget pt2))
    
66.   (setq ne1 (ssname e1 0))
    
67.   (setq ne2 (ssname e2 0))
    
68.   (command "fillet" pt1 pt2)
    
69.   (setq h10 (cdr (assoc 10 (entget (entlast)))))
    
70.   (setq h40 (cdr (assoc 40 (entget (entlast)))))
    
71.   (setq h50 (cdr (assoc 50 (entget (entlast)))))
    
72.   (setq h51 (cdr (assoc 51 (entget (entlast)))))
    
73.   (setq hd  (list (+ (car h10) (* h40 (cos h50))) (+ (cadr h10) (* h40 (sin h50)))))
    
74.   (setvar "filletrad" 0)
    
75.   (command "fillet" pt1 pt2)
    
76.   (setq p11 (cdr (assoc 10 (entget ne1))))
    
77.   (setq p12 (cdr (assoc 11 (entget ne1))))
    
78.   (setq p21 (cdr (assoc 10 (entget ne2))))
    
79.   (setq p22 (cdr (assoc 11 (entget ne2))))
    
80.   (if (< (dist pt1 p11) (dist pt1 p12))
    
81.     (progn
    
82.       (setq p11a p11)
    
83.       (setq p12a p12)
    
84.     )
    
85.     (progn
    
86.       (setq p11a p12)
    
87.       (setq p12a p11)
    
88.     )
    
89.   )
    
90.   (if (< (dist pt2 p21) (dist pt2 p22))
    
91.     (progn
    
92.       (setq p21a p21)
    
93.       (setq p22a p22)
    
94.     )
    
95.     (progn
    
96.       (setq p21a p22)
    
97.       (setq p22a p21)
    
98.     )
    
99.   )
    
100.   (if (= qstext nil) (setq qstext 1))
     
101.   (setq text (getint (strcat "\n***交点序号TD <" (rtos qstext 2 0) "> : ")))
     
102.   (if (= text nil)
     
103.     (setq text qstext)
     
104.     (setq qstext text)
     
105.   )
     
106.   (setq text (strcat "TD" (rtos text 2 0)))
     
107.   (setvar "clayer" "T1")
     
108.   (command "text" p11a (* scale 3) "0" text)
     
109. 
     
110.   (setq ang1 (angle p11a p12a))
     
111.   (setq ang2 (angle p21a p22a))
     
112.   (setq ang (abs (/ (* 180.0 (- ang2 ang1)) pi)))
     
113.   (setq angstr (rtos ang 2 2))
     
114.   (princ (strcat "角度=" angstr ))
     
115.   (setq rstr (rtos r 2 3))
     
116.   (princ (strcat "r=" rstr ))
     
117.   (setq L (* (/ (* 2 pi r) 360) (- 180 ang)))
     
118.   (setq lstr (rtos l 2 3))
     
119.   (princ (strcat "弧长=" lstr ))
     
120.   
     
121.   (setq tt (distance hd p11a))
     
122.   (setq ttstr (rtos tt 2 3))
     
123.   (princ (strcat "t=" ttstr ))
     
124.   
     
125.   (setq delta (- (* 2 tt) l))
     
126.   (setq deltastr (rtos delta 2 3))
     
127.   (princ (strcat "delta=" deltastr ))
     
128.   (setq text1 (strcat angstr "  " rstr "  " lstr "  " ttstr "  " deltastr))
     
129.   (setvar "clayer" "T1")
     
130.   (command "text" h10 (* scale 3) "0" text1)
     
131.   
     
132.   (setvar "osmode" os)
     
133.   (setvar "cmdecho" 1)
     
134.   (setvar "clayer" old_lay)
     
135.   (setvar "dimzin" dimzin_old)
     
136.   (princ)
     
137. )
     
138. 
     
139.   [/FONT]
     

<br />   [FONT=courier new]<br /> (defun dist (p1 p2 / dis)<br />   (setq dis (distance p1 p2))<br /> )<br /> <br /> <br /> (defun zhberr (msg)<br />   (setvar "osmode" os)<br />   (setvar "cmdecho" 1)<br />   (setvar "clayer" old_lay)<br />   (setvar "dimzin" dimzin_old)<br />   (setq *error* olderr)<br />   (princ)<br /> )<br /> <br /> (defun c:yh (/ pt1 pt2 p11 p12 p21 p22 scale e1 e2 ne1 ne2 text text1)<br />   (command "color" "")<br />   (setq os (getvar "osmode")<br />         olderr *error*<br />         *error* zhberr<br />   )<br />   (princ "\n*** 路线圆弧连接V1.0，上海宝钢工程技术有限公司 周海兵编 ***")<br />   (setvar "cmdecho" 0)<br />   (setq old_lay (getvar "clayer"))<br />   (setq dimzin_old (getvar "Dimzin"))<br />   (setvar "dimzin" 8)<br />   (if (= (tblobjname "LAYER" "T1") nil)<br />     (entmake (list '(0 . "LAYER") '(100 . "AcDbSymbolTableRecord") '<br />                    (100 . "AcDbLayerTableRecord") '(6 . "CONTINUOUS") '<br />                    (62 . 3) '(70 . 0) (cons 2 "T1")<br />              )<br />     )<br />   )<br />   (if (= (tblobjname "LAYER" "T2") nil)<br />     (entmake (list '(0 . "LAYER") '(100 . "AcDbSymbolTableRecord") '<br />                    (100 . "AcDbLayerTableRecord") '(6 . "CONTINUOUS") '<br />                    (62 . 5) '(70 . 0) (cons 2 "T2")<br />              )<br />     )<br />   )<br />   (if (= nil scaleqs) (setq scaleqs 1000))<br />         (setq scale (getreal (strcat "\n *** 初始化比例 <" (rtos scaleqs 2 0) "> *** :")))<br />         (if (= scale nil)<br />                 (setq scale scaleqs)<br />                 (setq scaleqs scale)<br />         )<br />   (setvar "osmode" 512)<br />   (if (= nil r0)<br />     (setq r0 (getvar "filletrad"))<br />   )<br />   (setq pt1 (getpoint "\n选择第一条线: "))<br />   (setq pt2 (getpoint "\n选择第二条线: "))<br />   (setvar "osmode" 0)<br />   (prompt (strcat "\n输入圆角半径<" (rtos r0 2 0) "> m :"))<br />   (initget 4)<br />   (setq r (getreal))<br />   (if (= nil r)<br />     (setq r r0)<br />     (setq r0 r)<br />   )<br />   (setvar "filletrad" r)<br />   (setq e1 (ssget pt1))<br />   (setq e2 (ssget pt2))<br />   (setq ne1 (ssname e1 0))<br />   (setq ne2 (ssname e2 0))<br />   (command "fillet" pt1 pt2)<br />   (setq h10 (cdr (assoc 10 (entget (entlast)))))<br />   (setq h40 (cdr (assoc 40 (entget (entlast)))))<br />   (setq h50 (cdr (assoc 50 (entget (entlast)))))<br />   (setq h51 (cdr (assoc 51 (entget (entlast)))))<br />   (setq hd  (list (+ (car h10) (* h40 (cos h50))) (+ (cadr h10) (* h40 (sin h50)))))<br />   (setvar "filletrad" 0)<br />   (command "fillet" pt1 pt2)<br />   (setq p11 (cdr (assoc 10 (entget ne1))))<br />   (setq p12 (cdr (assoc 11 (entget ne1))))<br />   (setq p21 (cdr (assoc 10 (entget ne2))))<br />   (setq p22 (cdr (assoc 11 (entget ne2))))<br />   (if (< (dist pt1 p11) (dist pt1 p12))<br />     (progn<br />       (setq p11a p11)<br />       (setq p12a p12)<br />     )<br />     (progn<br />       (setq p11a p12)<br />       (setq p12a p11)<br />     )<br />   )<br />   (if (< (dist pt2 p21) (dist pt2 p22))<br />     (progn<br />       (setq p21a p21)<br />       (setq p22a p22)<br />     )<br />     (progn<br />       (setq p21a p22)<br />       (setq p22a p21)<br />     )<br />   )<br />   (if (= qstext nil) (setq qstext 1))<br />   (setq text (getint (strcat "\n***交点序号TD <" (rtos qstext 2 0) "> : ")))<br />   (if (= text nil)<br />     (setq text qstext)<br />     (setq qstext text)<br />   )<br />   (setq text (strcat "TD" (rtos text 2 0)))<br />   (setvar "clayer" "T1")<br />   (command "text" p11a (* scale 3) "0" text)<br /> <br />   (setq ang1 (angle p11a p12a))<br />   (setq ang2 (angle p21a p22a))<br />   (setq ang (abs (/ (* 180.0 (- ang2 ang1)) pi)))<br />   (setq angstr (rtos ang 2 2))<br />   (princ (strcat "角度=" angstr ))<br />   (setq rstr (rtos r 2 3))<br />   (princ (strcat "r=" rstr ))<br />   (setq L (* (/ (* 2 pi r) 360) (- 180 ang)))<br />   (setq lstr (rtos l 2 3))<br />   (princ (strcat "弧长=" lstr ))<br />   <br />   (setq tt (distance hd p11a))<br />   (setq ttstr (rtos tt 2 3))<br />   (princ (strcat "t=" ttstr ))<br />   <br />   (setq delta (- (* 2 tt) l))<br />   (setq deltastr (rtos delta 2 3))<br />   (princ (strcat "delta=" deltastr ))<br />   (setq text1 (strcat angstr "  " rstr "  " lstr "  " ttstr "  " deltastr))<br />   (setvar "clayer" "T1")<br />   (command "text" h10 (* scale 3) "0" text1)<br />   <br />   (setvar "osmode" os)<br />   (setvar "cmdecho" 1)<br />   (setvar "clayer" old_lay)<br />   (setvar "dimzin" dimzin_old)<br />   (princ)<br /> )<br /> <br />   [/FONT]<br />

yshf

你的目的是求圆曲线的要素，将“rtl.fas”下载后存入G盘中，在程序中或命令行作如下调用即可：
(if (null rtl) (load "g:rtl.fas"));存入其它地方请注意路径
(rtl)

Andyamy

OK，谢谢。基本上应该是这样。我试试。
试验结果不错！再次谢谢。

这是在X、Y轴比例一致的情况下即都是1：1的条件下求得的A角度，如果X轴比例为1：2000，Y轴比例为1：100，那么该如何计算A角的实际角度

请求帮助

yshf

;请将下载的“rtl.fas”存入G盘中，并试用以下程序
(defun c:cc()
   (if (null rtl)(load "g:rtl.fas"));如存入其它地方，请更改相应的路径名
   (setq xsca (getreal "\n请输入X轴比例尺分母<1>:")
         ysca (getreal "\n请输入Y轴比例尺分母<1>:")
   )
   (if (= xsca nil) (setq xsca 1))
   (if (= ysca nil) (setq ysca 1))
   (setq xybl (list (/ xsca 1000.0) (/ ysca 1000.0)))
   (while (setq r (getreal "\n请输入实际的转角半径R:"))
      (initget "L A")
      (setq fhb nil
            nbz (getkword "\n指定两直线(L)/指定圆弧(A)[L/A]:")
      )
      (cond ((= nbz "L")
             (if (and (setq en1  (car (entsel "\n请指定直线a:")))
                      (setq en2  (car (entsel "\n请指定直线b:")))
                 )
                 (setq fhb (rtl1 en1 en2 r xybl));fhb=(转角 半径 切线长 曲线长 外矢距 要素字符串)
                 (alert "只指定了一条直线或者两条直线均未指定！")
              )
            )
            ((= nbz "A")
             (if (setq en1  (car (entsel "\n请指定圆弧:")))
                 (setq fhb (rtl2 en1 r xybl));fhb=(转角 半径 切线长 曲线长 外矢距 要素字符串)
                 (alert "未指定圆弧！")
              )
               )
      )
      (if fhb (alert (last fhb)))
   )            
)