Area of ellipse not in xy-plane

753 Views Asked by At

I've got a problem in which I'm trying to find the area of an ellipse which is given by the intersection of an elliptic cylinder with a plane. Nothing here is parallel to the coordinate axes, which is making it kind of annoying to work with.

The plane is given by the equation $x+ay+a^2z=0$, and the cylinder is given by $(x-a^2z)^2-a(x-a^2z)(y-az)+a^2(y-az)^2=L^2$.

I can think of some complicated ways to do it with integrals, and I'm wondering if there's something simpler that I'm missing. If I could transform into some coordinates in the plane I'm trying to work with, and I knew that the coordinate transformation would preserve areas, then I would be good, because I know how to get the area of an ellipse in the form $Ax^2+Bxy+Cy^2=1$. I'm really not sure how to do this transformation, though, and whether this is even the best way to proceed.

Maybe I should be using Langrange multipliers with two constraints to just obtain the lengths of the semi-major and semi-minor axes? That sounds like a pain, but doable. At least we're centered at the origin.

Thanks in advance for any assistance.


Edit: If I take $u=(x-a^2z)$ and $v=(y-az)$, then my cylinder becomes $u^2-auv+a^2v^2=L^2$, which is pretty great, but then I've basically cut the cylinder with a plane parallel to $z=0$. In that case, my given plane becomes $u+av+3a^2z=0$, and I'm not sure how that's helpful.

2

There are 2 best solutions below

4
On BEST ANSWER

Expanding the equation of the cylinder, we get $$ x^2-axy+a^2y^2-a^2xz-a^3yz+a^4z^2 = L^2 $$ Writing this as a quadratic form: $$ \begin{bmatrix}x\\y\\z\end{bmatrix}^T \begin{bmatrix}1 & -\frac{a}{2} & -\frac{a^2}{2} \\ -\frac{a}{2} & a^2 & -\frac{a^3}{2} \\ -\frac{a^2}{2} & -\frac{a^3}{2} & a^4 \end{bmatrix} \begin{bmatrix}x\\y\\z\end{bmatrix} = L^2 $$ Diagonalizing this matrix, we find that the eigenvalue $0$ is associated with the eigenvector $\vec{v} = (a^2,a,1)^T$, giving the axis of the cylinder. The other two eigenvalues $\lambda_1$ and $\lambda_2$ are related to the semi-major axes of the ellipse by $s_{1,2} = L/\sqrt{ \lambda_{1,2}}$. The other eigenvalues are: $$ \lambda_{1,2} = \frac{1}{2}\left(1+a^2+a^4\pm\sqrt{(a^2-1)^2(a^4+a^2+1)} \right)$$

The area of the cross section of this cylinder is $$A_0=\pi s_1s_2 = \frac{\pi L^2}{\sqrt{\lambda_1 \lambda_2}} = \frac{2\pi L^2}{a\sqrt{3(1+a^2+a^4)}}$$

The plane cuts the cylinder at an angle, and the area of the slanted cut is $A = A_0/\cos \theta$ where $\theta$ is the angle between the plane's normal axis and the cylinder axis. The plane's normal vector is $\vec{n} = (1,a,a^2)^T$. Using a dot product, $$\cos\theta = \frac{\vec{n}\cdot\vec{v}}{\left\lVert \vec{n} \right\rVert \left\lVert \vec{v} \right\rVert} = \frac{3a^2}{1+a^2+a^4}$$ So then, $$ A = \frac{2\pi L^2}{a\sqrt{3(1+a^2+a^4)}} \frac{1+a^2+a^4}{3a^2} = \frac{2\pi L^2}{3a^3\sqrt{3}} \sqrt{1+a^2+a^4}$$


Ignore this stuff: As OP pointed out, projecting out $x$ doesn't make any sense.

Substituting away $x$ using the plane equation $x=-ay-a^2 z$ into the cylinder equation gives $$ a^2(y+2az)^2+a^2(y+2az)(y-az)+a^2(y-az)^2=L^2 $$ $$ (y+2az)^2+(y+2az)(y-az)+(y-az)^2=\left(\frac{L}{a}\right)^2 $$ Let $u=y+2az$ and $v=y-az$, so $$ u^2+uv+v^2 = \left(\frac{L}{a}\right)^2 $$ The Jacobian of this transformation is $-3a$. The area of the ellipse in $u,v$ coordinates can be found using the second formula here. To obtain the area in the original $x,y,z$ coordinates, multiply by the Jacobian of the transformation.

0
On

I think this is less difficult than it appears. If you take any figure with area $S$ in a plane and projected it orthogonally onto another plane, the area of the image of your figure will be $S \cos \theta$, where $\theta$ is the angle between the planes. (Equivalently, $\theta$ is the angle your direction of projection makes with a vector normal to the plane you projected from.)

More generally, if you project the figure obliquely onto another plane, then the image has area $S \cos\theta\sec\phi$ where $\phi$ is the angle between the direction of projection and a vector normal to the plane you projected to.

Now consider the intersection of the cylinder with the $x,y$-plane. This is an ellipse whose formula you can derive simply by setting $z = 0$ in the formula of the cylinder. You know how to find the area $S$ of that ellipse. A projection parallel to the axis of the cylinder takes that ellipse onto your desired ellipse. I believe that projection is an orthogonal projection onto the plane $x+ay+a^2z=0$, in which case all you need to do is find the angle $\theta$ between the axis of the cylinder and the $z$ axis, and then the area of the ellipse you want is $S \cos\theta$.

If I'm wrong about the axis of the cylinder being orthogonal to the plane $x+ay+a^2z=0$, then you also have to find the angle $\phi$ between the axis and that plane, and then the answer is $S \cos\theta\sec\phi$.