Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

4.3: The Ellipsoid

( \newcommand{\kernel}{\mathrm{null}\,}\)

Consider the Equation

x2a2+z2c2=1,

with a>c, in the xz-plane. The length of the semi major axis is a and the length of the semi minor axis is c. If this figure is rotated through 360 about its minor (z-) axis, the three- dimensional figure so obtained is called an oblate spheroid. The figure of the Earth is not exactly spherical; it approximates to a very slightly oblate spheroid, the ellipticity (ca)/a being only 0.00335. (The actual figure of the Earth, mean sea level, is often referred to as the geoid.)

The Equation to the oblate spheroid referred to above is

x2a2+y2a2+z2c2=1.

If the ellipse 4.3.1 is rotated through 360 about its major (x-) axis, the figure so obtained is called a prolate spheroid. A rugby football (or, to a lesser extent, a North American football, which is a bit too pointed) is a good approximation to a prolate spheroid.

The Equation to the prolate spheroid just described is

x2a2+y2c2+z2c2=1.

Either type of spheroid can be referred to as an "ellipsoid of revolution".

The figure described by the Equation

x2a2+y2b2+z2c2=1

is a tri-axial ellipsoid. Unless stated otherwise, I shall adopt the convention a>b>c, and choose the coordinate axes such that the major, intermediate and minor axes are along the x-, y- and z-axes respectively. A tri-axial ellipsoid is not an ellipsoid of revolution; it cannot be obtained by rotating an ellipse about an axis.

The special case a=b=c:

x2+y2+z2=a2

is, of course, a sphere.

Figure IV.4 shows the cross-section of a tri-axial ellipse in the yz- plane (a), the xz-plane (b) and (twice - (c), (d)) the xy-plane. If you imagine your eye wandering in the xz-plane from the x-axis (a) to the z-axis (c), you will be convinced that there is a direction in the xz-plane from which the

Figure 4.4.png
FIGURE IV.4

cross-section of the ellipse is a circle. There are actually two such directions, symmetrically situated on either side of the z-axis, but there are no such directions in either the xy- or the yz-planes from which the cross-section of the ellipsoid appears as a circle. Expressed otherwise, there are two planes that intersect the ellipsoid in a circle. This fact is of some importance in the description of the propagation of light in a bi-axial crystal, in which one of the wavefronts is a tri-axial ellipsoid.

Let us refer the ellipsoid 4.3.4 to a set of axes Oxyz such that the angles zOz and xOx are each θ, and the y- and y-axes are identical. The Equation of the ellipsoid referred to the new axes is (by making use of the usual formulas for the rotation of axes)

(zsinθ+xcosθ)2a2+y2b2+(zcosθxsinθ)2c2=1.

The cross-section of the ellipsoid in the xy-plane (i.e. normal to the z-axis) is found by putting z=0:

(xcosθ)2a2+y2b2+(xsinθ)2c2=1.

This is a circle if the coefficients of x and y are equal. Thus it is a circle if

cos2θ=a2(b2c2)b2(a2c2).

Thus, a plane whose normal is in the xz-plane (i.e. between the major and minor axis) and inclined at an angle θ to the minor (z-) axis, cuts the tri-axial ellipsoid in a circle. As viewed from either of these directions, the cross-section of the ellipsoid is a circle of radius b.

As an asteroid tumbles over and over, its brightness varies, for several reasons, such as its changing phase angle, the directional reflective properties of its regolith, and, of course, the cross-sectional area presented to the observer. The number of factors that affect the light-curve of a rotating asteroid is, in fact, so large that it is doubtful if it is possible, from the light-curve alone, to deduce with much credibility or accuracy the true shape of the asteroid. However, it is obviously of some interest for a start in any such investigation to be able to calculate the cross-sectional area of the ellipsoid 4.3.3 as seen from some direction (θ,ϕ).

Let us erect a set of coordinate axes Oxyz such that Oz is in the direction (θ,ϕ), first by a rotation through ϕ about Oz to form intermediate axes Ox1y1z1, followed by a rotation through θ about Oy1. The (x,y,z) coordinates are related to the (x,y,z) coordinates by

(xyz)=(cosϕsinϕ0sinϕcosϕ0001)(cosθ0sinθ010sinθ0cosθ)(xyz)

If we substitute for x, y, z in Equation 4.3.4 from Equation 4.3.9, we obtain the Equation to the ellipsoid referred to the Oxyz coordinate systems. And if we put z=0, we see the elliptical crosssection of the ellipsoid in the plane normal to Oz. This will be of the form

Ax2+2Hxy+By2=1,

where A=cos2θ(cos2ϕa2+sin2ϕb2)+sin2θc2

2H=2cosθsinϕcosϕ(1b21a2),

B=sin2ϕa2+cos2ϕb2.

This is an ellipse whose axes are inclined at an angle ψ from Ox given by

tan2ψ=2HAB.

By replacing x and y by x and y, where

(xy)=(cosψsinψsinψcosψ)(xy)

we shall be able to describe the ellipse in a coordinate system Oxy whose axes are along the axes of the ellipse, and the Equation will be of the form

x2a2+y2b2=1

and the area of the cross-section is πab.

For example, suppose the semi axes of the ellipsoid are a=3, b=2, y=1, and we look at it from the direction θ=60, ϕ=45. Following Equations 4.4.9,10,11,12, we obtain for the Equation of the elliptical cross-section referred to the system Oxyz

0.79513˙8x2+0.069˙4xy+0.180˙5y2=1.

From Equation 4.4.13 we find ψ=3 .22338. Equation 4.4.14 then transforms Equation 4.4.16 to

0.797094x2+0.178600y2=1

or x2(1.1201)2+y2(2.3662)2=1.

The area is π×1.1201×2.3662=8.362.

It is suggested here that the reader could write a computer program in the language of his or her choice for calculating the cross-sectional area of an ellipsoid as seen from any direction. As an example, I reproduce below a Fortran program for an ellipse with (a,b,c)=(3,2,1). It is by no means the fastest and most efficient Fortran program that could be written, but is sufficiently straightforward that anyone familiar with Fortran and probably many who are not should be able to follow the steps.

A=3.
B=2.
C=1.
A2=A*A
B2=B*B
C2=C*C
READ(5,*)TH,PH
TH=TH/57.29578
PH=PH/57.29578
STH=SIN(TH)
CTH=COS(TH)
SPH=SIN(PH)
CPH=COS(PH)
STH2=STH*STH
CTH2=CTH*CTH
SPH2=SPH*SPH
CPH2=CPH*CPH
AA=CTH2*(CPH2/A2+SPH2/B2)+STH2/C2
TWOHH=2.*CTH*SPH*CPH*(1./B2−1./A2)
BB=SPH2/A2+CPH2/B2
PS=.5*ATAN2(TWOHH,AA−BB)
SPS=SIN(PS)
CPS=COS(PS)
AAA=CPS*(AA*CPS+TWOHH*SPS)+BB*SPS*SPS
BBB=SPS*(AA*SPS−TWOHH*CPS)+BB*CPS*CPS
SEMAX1=1./SQRT(AAA)
SEMAX2=1./SQRT(BBB)
AREA=3.1415927*SEMAX1*SEMAX2
WRITE(6,1)AREA
1 FORMAT(' Area = ',F7.3)
STOP
END


This page titled 4.3: The Ellipsoid is shared under a CC BY-NC 4.0 license and was authored, remixed, and/or curated by Jeremy Tatum via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?