[ADN开发资料]：椭圆弧组码41的意义...

XDSoft

Explanation of elliptical arcs group code 41  
ID    467  
Applies to:    AutoCAD

Date    6/22/2000  

This document is part of    Visual Lisp   DXF Issues   Ellipse   Group Codes     


Question
In elliptical arcs, group code 41 is described as the start parameter. However,
it doesn't seem to be associated with any of the known values for the arc. How
is this value derived?
Answer
An elliptical arc is a special version of an arc that follows the eccentricity
of the ellipse.  One way to generate this type of arc is to find the parametric
normal of the starting point. To do so, you must specify a start angle that is
different from the actual start angle of the drawn arc. Group code 41 contains
this parametric angle expressed in radians.

WHAT IS THE PARAMETRIC ANGLE?

The parametric angle is generated from two concentric circles whose center is
the center of the ellipse and whose radii are the major and minor axes,
respectively. Every point of the ellipse lies either between or on these two
circles, and each elliptical point can be defined by a unique relation to them.
To discover this relationship, draw a line perpendicular to the major axis from
a point on the ellipse to the closest intersection with the circle described by
the major axis. Then do the same with the minor axis, starting from the
elliptical point and drawing perpendicular to the minor axis until the line
intersects the circle described by the minor axis. The two points of
intersection with the circles are colinear with the center of the ellipse. The
angle between the line containing these three points and the major axis is the
parametric angle specified by group code 41.   

HOW DO I CALCULATE IT FROM THE TRUE START ANGLE?

To calculate the parametric angle from the true start angle, you must first find
the start point on the ellipse. This requires a simultaneous solution to the
equations for the line and the ellipse. In this example we assume that the major
axis of the ellipse lies on the x-axis with the origin at the center of the
ellipse. When this point is found, you can use its y-value and the minor axis to
solve the equation for the circle whose radius is the minor axis value and whose
center is the center of the ellipse. This will provide the x,y point on the
circle that dictates the parametric angle from the center of the ellipse.

The following is an AutoLisp example that demonstrates how to use trigonometric
functions to determine the parametric angle:

2. 
   
3. [FONT=courier new]
   
4. (defun c:e_arc( / a b slope ang q1 q2 q3 q4 qmode x y a2)
   
5. ;assuming 0,0 is at the center of the ellipse, major axis in x direction
   
6.         (setq ang (getangle '(0.0 0.0) "Choose start angle: "))
   
7.         (setq a 1
   
8.                 b 0.5
   
9.                 slope (/ (sin ang) (cos ang))
   
10.                 q1 (/ pi 2.0)
    
11.                 q2 pi
    
12.                 q3 (/ (* 3 pi) 2.0)
    
13.                 q4 (* 2.0 pi)
    
14.                 qmode 'q1
    
15.         );setq
    
16.         (entmake (setq ent '((0 . "ELLIPSE")
    
17.                         (100 . "AcDbEntity")
    
18.                         (100 . "AcDbEllipse")
    
19.                         (10 0.0 0.0 0.0)
    
20.                         (11 1.0 0.0 0.0)
    
21.                         (40 . 0.5)
    
22.                         (62 . 1))
    
23.                 );setq
    
24.         );entmake
    
25. 
    
26. ;line equation is y = mx + 0, where m is the slope and 0 is the y-intercept
    
27. ;ellipse equation is x^2/a^2 + y^2/b^2 = 1
    
28. ;solve line and ellipse equations simultaneously to find x and y values
    
29. 
    
30.         (setq y (/ (* a b slope) (sqrt (+ (* (* slope slope) (* a a)) (* b b))))
    
31.         );setq
    
32. 
    
33. ;minor axis circle equation is x^2 + y^2 = b^2
    
34. ;solve circle equation where y = value calculated above
    
35.         (setq x (sqrt (- (* b b) (* y y))))
    
36. 
    
37. ;calculate start angle trigonometrically
    
38.         (setq cos_a2 (/ x b)
    
39.               sin_a2 (/ y b)
    
40.         );setq
    
41.         (if (/= cos_a2 0)
    
42.                 (setq a2 (atan (/ sin_a2 cos_a2)))
    
43.                 (setq a2 q1 qmode 'q1)                         
    
44.         );if
    
45. 
    
46. 
    
47. ;make a2 insensitive to quadrant
    
48.         (cond ((and (> ang q1) (< ang q2))
    
49.                                 (setq a2 (- pi (abs a2)) qmode 'q2
    
50.                                 );setq
    
51.                         );statement 1
    
52.                         ((and (> ang q2) (< ang q3))
    
53.                                 (setq a2 (+ (abs a2) pi) qmode 'q3
    
54.                                 );setq
    
55.                         );statement 2
    
56.                         ((and (> ang q3) (< ang q4))
    
57.                                 (setq a2 (abs (- (* 2 pi) (abs a2))) qmode 'q4
    
58.                                 );setq
    
59.                         );statement 3
    
60. 
    
61.         ;special cases: angle = 0, 90, 180, 270 or 360 deg
    
62.                         ((or (= ang 0) (= ang q1))
    
63.                                 (setq qmode 'q1)
    
64.                         );statement 4
    
65.                         ((= ang q2)
    
66.                                 (setq a2 pi qmode 'q1)
    
67.                         );statement 5
    
68.                         ((= ang q3)
    
69.                                 (setq a2 (- (/ pi 2.0) pi) qmode 'q1)
    
70.                         );statement 6
    
71.                         (t nil);default statement
    
72.         );cond
    
73. 
    
74.         (command "zoom" "c" "0,0" 3)               
    
75.         (setq ent (append  ent (list (cons 41 a2) (cons 42 (+ a2 (/ pi 2.0))))))
    
76.         (setq ent (subst '(62 . 5) (assoc 62 ent) ent))
    
77.         (entmake ent)
    
78.        
    
79.         (setq a2
    
80.                 (cond ((= qmode 'q1) a2)
    
81.                         ((= qmode 'q2) a2)
    
82.                         ((= qmode 'q3) (- a2 q4))
    
83.                         ((= qmode 'q4) (- a2 q4))
    
84.                         (t nil)
    
85.                 );cond
    
86.         );setq
    
87.        
    
88.         (princ "\nParametric angle in radians: ")
    
89.         (princ a2)
    
90.         (princ "\nParametric angle in degrees: ")
    
91.         (princ (/ (* 180 a2) pi))
    
92.         (princ)
    
93.        
    
94. );e_arc[/FONT]

<br /> [FONT=courier new]<br /> (defun c:e_arc( / a b slope ang q1 q2 q3 q4 qmode x y a2)<br /> ;assuming 0,0 is at the center of the ellipse, major axis in x direction<br />         (setq ang (getangle '(0.0 0.0) "Choose start angle: "))<br />         (setq a 1 <br />                 b 0.5 <br />                 slope (/ (sin ang) (cos ang))<br />                 q1 (/ pi 2.0)<br />                 q2 pi<br />                 q3 (/ (* 3 pi) 2.0)<br />                 q4 (* 2.0 pi)<br />                 qmode 'q1<br />         );setq<br />         (entmake (setq ent '((0 . "ELLIPSE") <br />                         (100 . "AcDbEntity") <br />                         (100 . "AcDbEllipse") <br />                         (10 0.0 0.0 0.0) <br />                         (11 1.0 0.0 0.0)<br />                         (40 . 0.5) <br />                         (62 . 1))<br />                 );setq<br />         );entmake<br /> <br /> ;line equation is y = mx + 0, where m is the slope and 0 is the y-intercept<br /> ;ellipse equation is x^2/a^2 + y^2/b^2 = 1<br /> ;solve line and ellipse equations simultaneously to find x and y values<br /> <br />         (setq y (/ (* a b slope) (sqrt (+ (* (* slope slope) (* a a)) (* b b))))<br />         );setq<br /> <br /> ;minor axis circle equation is x^2 + y^2 = b^2<br /> ;solve circle equation where y = value calculated above<br />         (setq x (sqrt (- (* b b) (* y y))))<br /> <br /> ;calculate start angle trigonometrically<br />         (setq cos_a2 (/ x b)<br />               sin_a2 (/ y b)<br />         );setq<br />         (if (/= cos_a2 0)<br />                 (setq a2 (atan (/ sin_a2 cos_a2)))<br />                 (setq a2 q1 qmode 'q1)                         <br />         );if<br /> <br /> <br /> ;make a2 insensitive to quadrant<br />         (cond ((and (> ang q1) (< ang q2))<br />                                 (setq a2 (- pi (abs a2)) qmode 'q2<br />                                 );setq<br />                         );statement 1<br />                         ((and (> ang q2) (< ang q3))<br />                                 (setq a2 (+ (abs a2) pi) qmode 'q3<br />                                 );setq<br />                         );statement 2<br />                         ((and (> ang q3) (< ang q4))<br />                                 (setq a2 (abs (- (* 2 pi) (abs a2))) qmode 'q4<br />                                 );setq<br />                         );statement 3<br /> <br />         ;special cases: angle = 0, 90, 180, 270 or 360 deg<br />                         ((or (= ang 0) (= ang q1))<br />                                 (setq qmode 'q1)<br />                         );statement 4<br />                         ((= ang q2) <br />                                 (setq a2 pi qmode 'q1)<br />                         );statement 5<br />                         ((= ang q3)<br />                                 (setq a2 (- (/ pi 2.0) pi) qmode 'q1)<br />                         );statement 6<br />                         (t nil);default statement<br />         );cond<br /> <br />         (command "zoom" "c" "0,0" 3)                <br />         (setq ent (append  ent (list (cons 41 a2) (cons 42 (+ a2 (/ pi 2.0))))))<br />         (setq ent (subst '(62 . 5) (assoc 62 ent) ent))<br />         (entmake ent)<br />         <br />         (setq a2 <br />                 (cond ((= qmode 'q1) a2)<br />                         ((= qmode 'q2) a2)<br />                         ((= qmode 'q3) (- a2 q4))<br />                         ((= qmode 'q4) (- a2 q4))<br />                         (t nil)<br />                 );cond<br />         );setq<br />         <br />         (princ "\nParametric angle in radians: ")<br />         (princ a2)<br />         (princ "\nParametric angle in degrees: ")<br />         (princ (/ (* 180 a2) pi))<br />         (princ)<br />         <br /> );e_arc[/FONT]