【越飞越高讲堂14】椭圆论

crtrccrt

一个字 太牛了 太牛了

海盗曹

向前辈学习啦

highflybir

更新了一些问题的解。如：十一题 最似圆椭圆问题

MENGZE

版主，弱弱问下在CAD模型空间中椭圆形椎体怎么展开？求回复！！！

wangdaobin

严重支持高师傅！！

highflybir

MENGZE 发表于 2013-5-4 00:28
版主，弱弱问下在CAD模型空间中椭圆形椎体怎么展开？求回复！！！

1. ;;;=====================================================================
   
2. ;;; 功能: 展开椭圆锥体。                                   
   
3. ;;; 输入: 椭圆锥底椭圆的长轴a，短轴b，椭圆锥的高h和分弧精度n            
   
4. ;;; 输出: 展开图的曲线的坐标点(可据此形成展开面)                        
   
5. ;;;=====================================================================
   
6. (defun ELL:ExpandCone (a b h n / Pts Ang D Param DivAng P0 P1 R0 R1 X Y)
   
7.   (setq DivAng (/ pi n 0.5))      ;等分角度
   
8.   (setq Param 0)       ;开始参数为0
   
9.   (setq p0 (list a 0))       ;从椭圆最右边的点开始
   
10.   (setq r0 (distance '(0 0 0) (list a 0 h)))    ;开始的锥顶最右边点的距离
    
11.   (setq Pts (list (list r0 0) '(0 0)))     ;把最开始的两点加入到数据表
    
12.   (setq Ang 0)        ;开始角度也为0
    
13.   (repeat n
    
14.     (setq Param  (+ Param DivAng))
    
15.     (setq x   (* a (cos Param)))     ;椭圆上的x坐标
    
16.     (setq y   (* b (sin Param)))     ;椭圆上的y坐标
    
17.     (setq p1  (list x y))      ;点的坐标
    
18.     (setq r1  (distance '(0 0 0) (list x y h)))    ;展开面的锥线长度
    
19.     (setq d   (distance p0 p1))      ;椭圆上的上一点到这点的距离
    
20.     (setq Ang (+ Ang (car (TRI:CosinesLaw d r0 r1))))   ;坐标角度
    
21.     (setq Pts (cons (polar '(0 0) Ang r1) Pts))    ;得到了坐标，并把它加入到表中
    
22.     (setq p0 p1)       ;用新的坐标位置替换旧位置
    
23.     (setq r0 r1)       ;用新的锥线长度替换旧长度
    
24.   )
    
25.   (reverse Pts)        ;逆转表，得到正序的数据
    
26. )
    
27. ;;; 测试椭圆锥的展开
    
28. (defun c:ExpandCone(/ h N e d r c m a b l e o)
    
29.   (initget 1)
    
30.   (setq h (getdist "\n输入椭圆锥高: "))
    
31.   (initget 7)
    
32.   (setq N (getint "\n等分数:"))
    
33.   (prompt "\n选择椭圆:")
    
34.   (if (setq s (ssget ":S" '((0 . "ELLIPSE"))))
    
35.     (progn
    
36.       (setq e (ssname s 0))
    
37.       (setq d (entget e))
    
38.       (setq r (cdr (assoc 40 d)))
    
39.       (setq c (cdr (assoc 10 d)))
    
40.       (setq m (cdr (assoc 11 d)))
    
41.       (setq a (distance '(0 0 0) m))
    
42.       (setq b (* r a))
    
43.       (setq l (ELL:ExpandCone a b h n))
    
44.       (setq o (vlax-ename->vla-object (Ent:Make_LWPoly l T)))
    
45.       (vla-move o (vlax-3d-point '(0 0 0)) (vlax-3d-point c))
    
46.       (princ)
    
47.     )
    
48.   )
    
49. )
    
50. ;;;----------------------------------------------------;
    
51. ;;; 功能: 根据三边求三个角(利用余弦定理)               ;
    
52. ;;; 输入: 构成三角形的三边a,b,c                        ;
    
53. ;;; 输出: 返会三条边对应的三个角                       ;
    
54. ;;;----------------------------------------------------;
    
55. (defun TRI:CosinesLaw (a b c / cc sc a1 a2)
    
56.   (if (TRI:IsTriangle a b c)
    
57.     (progn
    
58.       (setq a  (float a))     ;为了防止整除
    
59.       (setq cc (/ (+ (* a a) (* (+ b c) (- b c))) (+ a a) b)) ;角C的余弦
    
60.       (setq sc (sqrt (* (- 1 cc) (1+ cc))))    ;角C的正弦
    
61.       (setq a1 (atan (* a sc) (- b (* a cc))))    ;角A
    
62.       (setq a2 (atan sc cc))     ;角C
    
63.       (list a1 (- pi a1 a2) a2)     ;返回三个角度的列表
    
64.     )
    
65.   )
    
66. )
    
67. ;;;----------------------------------------------------;
    
68. ;;;创建一个三维多段线                                  ;
    
69. ;;;输入: 三维的点集                                    ;
    
70. ;;;输出: 三维多段线实体                                ;
    
71. ;;;----------------------------------------------------;
    
72. (defun Ent:Make_Poly (pts / e)
    
73.   (setq e (Entmake (list '(0 . "POLYLINE") '(70 . 9))))
    
74.   (foreach p pts
    
75.     (entmake (list '(0 . "VERTEX") '(70 . 32) (cons 10 p)))
    
76.   )
    
77.   (entmake '((0 . "SEQEND")))
    
78.   (entlast)
    
79. )
    

;;;=====================================================================<br /> ;;; 功能: 展开椭圆锥体。                                   <br /> ;;; 输入: 椭圆锥底椭圆的长轴a，短轴b，椭圆锥的高h和分弧精度n            <br /> ;;; 输出: 展开图的曲线的坐标点(可据此形成展开面)                        <br /> ;;;=====================================================================<br /> (defun ELL:ExpandCone (a b h n / Pts Ang D Param DivAng P0 P1 R0 R1 X Y)<br />   (setq DivAng (/ pi n 0.5))      ;等分角度<br />   (setq Param 0)       ;开始参数为0<br />   (setq p0 (list a 0))       ;从椭圆最右边的点开始<br />   (setq r0 (distance '(0 0 0) (list a 0 h)))    ;开始的锥顶最右边点的距离<br />   (setq Pts (list (list r0 0) '(0 0)))     ;把最开始的两点加入到数据表<br />   (setq Ang 0)        ;开始角度也为0<br />   (repeat n<br />     (setq Param  (+ Param DivAng))<br />     (setq x   (* a (cos Param)))     ;椭圆上的x坐标<br />     (setq y   (* b (sin Param)))     ;椭圆上的y坐标<br />     (setq p1  (list x y))      ;点的坐标<br />     (setq r1  (distance '(0 0 0) (list x y h)))    ;展开面的锥线长度<br />     (setq d   (distance p0 p1))      ;椭圆上的上一点到这点的距离<br />     (setq Ang (+ Ang (car (TRI:CosinesLaw d r0 r1))))   ;坐标角度<br />     (setq Pts (cons (polar '(0 0) Ang r1) Pts))    ;得到了坐标，并把它加入到表中<br />     (setq p0 p1)       ;用新的坐标位置替换旧位置<br />     (setq r0 r1)       ;用新的锥线长度替换旧长度<br />   )<br />   (reverse Pts)        ;逆转表，得到正序的数据<br /> )<br /> ;;; 测试椭圆锥的展开<br /> (defun c:ExpandCone(/ h N e d r c m a b l e o)<br />   (initget 1)<br />   (setq h (getdist "\n输入椭圆锥高: "))<br />   (initget 7)<br />   (setq N (getint "\n等分数:"))<br />   (prompt "\n选择椭圆:")<br />   (if (setq s (ssget ":S" '((0 . "ELLIPSE"))))<br />     (progn<br />       (setq e (ssname s 0))<br />       (setq d (entget e))<br />       (setq r (cdr (assoc 40 d)))<br />       (setq c (cdr (assoc 10 d)))<br />       (setq m (cdr (assoc 11 d)))<br />       (setq a (distance '(0 0 0) m))<br />       (setq b (* r a))<br />       (setq l (ELL:ExpandCone a b h n))<br />       (setq o (vlax-ename->vla-object (Ent:Make_LWPoly l T)))<br />       (vla-move o (vlax-3d-point '(0 0 0)) (vlax-3d-point c))<br />       (princ)<br />     )<br />   )<br /> )<br /> ;;;----------------------------------------------------;<br /> ;;; 功能: 根据三边求三个角(利用余弦定理)               ;<br /> ;;; 输入: 构成三角形的三边a,b,c                        ;<br /> ;;; 输出: 返会三条边对应的三个角                       ;<br /> ;;;----------------------------------------------------;<br /> (defun TRI:CosinesLaw (a b c / cc sc a1 a2)<br />   (if (TRI:IsTriangle a b c)<br />     (progn<br />       (setq a  (float a))     ;为了防止整除<br />       (setq cc (/ (+ (* a a) (* (+ b c) (- b c))) (+ a a) b)) ;角C的余弦<br />       (setq sc (sqrt (* (- 1 cc) (1+ cc))))    ;角C的正弦<br />       (setq a1 (atan (* a sc) (- b (* a cc))))    ;角A<br />       (setq a2 (atan sc cc))     ;角C<br />       (list a1 (- pi a1 a2) a2)     ;返回三个角度的列表<br />     )<br />   )<br /> )<br /> ;;;----------------------------------------------------;<br /> ;;;创建一个三维多段线                                  ;<br /> ;;;输入: 三维的点集                                    ;<br /> ;;;输出: 三维多段线实体                                ;<br /> ;;;----------------------------------------------------;<br /> (defun Ent:Make_Poly (pts / e)<br />   (setq e (Entmake (list '(0 . "POLYLINE") '(70 . 9))))<br />   (foreach p pts<br />     (entmake (list '(0 . "VERTEX") '(70 . 32) (cons 10 p)))<br />   )<br />   (entmake '((0 . "SEQEND")))<br />   (entlast)<br /> )<br />

MENGZE

谢谢老师！！！

MENGZE

受教了！！！

自贡黄明儒

C:test4中(setq ptb (vlax-curve-getPointAtParam ent (* pi 0.5)))与成了(setq ptb (vlax-curve-getPointAtParam en (* pi 0.5)))

NJZX05

鸟哥的东西，不赞都不行啊！！！！！要是我早看到，要少走多少弯路啊！！！