Register to reply 
Rotated Spheroid 
Share this thread: 
#1
Jul1505, 11:52 AM

P: 28

If I take a spheroid (slightly flattened sphere) with the following equation:
x2/a2 + y2/b2 + z2/a2  1 = 0 and rotate it by arbitrary angles about the x, y, and zaxes, the new equation for the rotated spheroid is: Ax2 + By2 + Cz2 + Dxy + Exz + Fyz + Gx + Hy + Iz + J = 0 How do I find the outline of its projection on the xyplane? 


#3
Jul1505, 02:33 PM

P: 28

I derive all equations by hand. I'm looking for a stepbystep solution to this problem.



#4
Jul1505, 06:17 PM

Sci Advisor
HW Helper
P: 9,495

Rotated Spheroid
ok. but sometimes equations for the image of a projection are a little hard to come by. thats why the mechanism of grobner bases and the program macaulay were written.



#5
Jul1605, 06:56 AM

P: 255

2Ax dx + 2By dy + 2Cz dz + Dx dy + Dy dx + Ex dz + Ez dx +Fy dz + Fz dy + G dx + H dy + I dz So, we need to find all of the points where, if we plug e_3 into the above 1form, we get zero. This boils down to: 2Cz + Ex + Fy + I = 0. So, you need to find the intersection of the surface with the above plane in 3space. This will be the curve (in 3space) that projects to the boundary of the projection on the xyplane  clearly some sort of ellipse, since it will be a conic section (its equation is quadratic in x and y) and it is bounded. 


#6
Jul1805, 08:55 AM

P: 28

Ok. But in your solution, the tangent plane has just as many unknowns as the original equation. Not a complete solution I'm afraid.



#7
Jul1805, 01:33 PM

P: 255

[tex]2Ax v_1 + 2By v_2 + 2Cz v_3 + Dx v_2 + Dy v_1 + Ex v_3 + Ez v_1 +Fy v_3 + Fz v_2 + G v_1 + H v_2 + I v_3 = 0 [/tex] which is in fact a 2diml. subspace of [tex]R^3[/tex] given the definitions of the constants A, B, ... , I (keep in mind that you have to regard (x,y,z) as fixed). The point we're looking for are those points on the surface which will have a tangent vector of (0,0,1). Plugging that into the above eqn. gives the plane equation that I got before. So, in other words, to find all of the points on your surface that project to the boundary of the domain in the xyplane, you need to solve this set of equations: Ax2 + By2 + Cz2 + Dxy + Exz + Fyz + Gx + Hy + Iz + J = 0 2Cz + Ex + Fy + I = 0 This is the intersection of two transverse surfaces, so the result will be a curve in 3space. Once you find those, the curve in the xyplane is just the same points, ignoring the zcoordinates of each. Not a very easy task. But I would like to point out that your original equation defines a much broader surface than just the spheroid. Any quadric surface (e.g. hyperboloids of one or two sheets, general ellipsoids, cones, etc.) satisfies this equation: Ax2 + By2 + Cz2 + Dxy + Exz + Fyz + Gx + Hy + Iz + J = 0. So, you might benefit by finding out some constraints on those constants. 


#8
Jul1905, 09:46 PM

P: 28

Ok. The surface normal is always perpendicular to the zaxis, and rotates about it, such that the partial derivative of the rotated spheroid equation F with respect to z is zero.
Fz(x,y,z) = Ex + Fy + 2Cz + I = 0 What you said. Solving for z: z = (I + Ex + Fy) / (2C) and substituting into F yields an ellipse of the form: Kx2 + Ly2 + Mxy + Nx + Oy + P = 0 So I obviously found the outline of the projection of the rotated spheroid. Converting to parametric form should allow me to find the extrema and plot the curve. Not so hard after all! 


#9
Jul2005, 05:00 AM

P: 255

You can easily argue that the result will be an ellipse without using any calculus whatsoever. However, you have not answered what type of ellipse you end up with. Clearly, the major and minor axes are limited by the values of a and b. But, is there a restriction on the eccentricity of the resulting ellipse? 


#10
Jul2005, 01:58 PM

P: 28

The constants would be known exactly by the original spheroid equation prior to rotation, i.e. semimajor and semiminor axes. These would be changed in a known and exact way by the angles of rotation, resulting in exact values for the constants of the final equations. Unless there is a problem with floating point errors (which is a distinct possibility) then constraints on the range of allowed values for the constants can be ignored. When converting to parametric form, the only constraint would be 0 <= t <= 2*pi. But that has nothing to do with eccentricity, more with the range of the independent variable t. Otherwise it wouldn't make sense to plot a curve outside of the range in which it exists. For example, you can't plot a unit circle for a range of xvalues greater than 1.
What were you suggesting? 


#11
Jul2005, 04:43 PM

P: 255

However, I haven't done it myself. So, it might indeed be easy. 


#12
Jul2005, 11:14 PM

P: 28

Looks like I made a few notsoobvious typos, and that was the cause of my problem. My methods were correct. Now that the typos are fixed, everything checks out, the results plot nicely. Cool. This problem had me stumped for quite a while.



Register to reply 
Related Discussions  
Angular momentum of an oblate spheroid  Astronomy & Astrophysics  4  
Volume of a rotated graph?  Calculus & Beyond Homework  8  
Equation of spheroid  Calculus & Beyond Homework  4  
Volume of rotated graph...  Calculus & Beyond Homework  4  
Deriving the Implicit Equation of a Rotated Spheroid  Differential Geometry  4 