Re: Trouble calculating a volume
- To: mathgroup at smc.vnet.net
- Subject: [mg33840] Re: [mg33806] Trouble calculating a volume
- From: Rob Pratt <rpratt at email.unc.edu>
- Date: Fri, 19 Apr 2002 02:28:56 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
On Tue, 16 Apr 2002, Jay D. Martin wrote: > Respected Colleagues: > > I am attempting to calculate the volume of the frustum of a cone. > I have been able to calculate the volume of a cone using a triple > integral as follows: > > Integrate[1, {x, 0, h}, {y, (-r)*(x/h), r*(x/h)}, > {z, -Sqrt[(r*(x/h))^2 - y^2], Sqrt[(r*(x/h))^2 - y^2]}] > > This integral does evaluate to the correct answer > > (1/3)*h*Pi*r^2 > > where r is the base diameter and h is the height. > > In order to calculate the volume of the frustum of a cone, I have used the > following equation : > > Simplify[Integrate[1, {x, 0, h}, > {y, -(r1 + ((r2 - r1)*x)/h), r1 + ((r2 - r1)*x)/h}, > {z, -Sqrt[(r1 + ((r2 - r1)*x)/h)^2 - y^2], > Sqrt[(r1 + ((r2 - r1)*x)/h)^2 - y^2]}]] > > > where r1 and r2 are the top and bottom radii of the section and h > is the height. This integral evaluates to include Logarithms and > imaginary numbers > > (I*Sqrt[h^2*r2]*(r1^2 + r1*r2 + r2^2)* > (Log[-2*I*h*r2] - Log[2*I*h*r2]))/(3*Sqrt[r2]) > > when the solution should be > > (1/3)*h*Pi*(r1^2+r1*r2+r2^2) > > Which, if you substitute r1=0 and r2=r you get the volume of > a cone. > > I am using Mathematica Version 4.1.2.0 > running on Windows 2000 > > Thank you for any help > > Jay Martin > > PS: what I am ultimately after are the moments of inertia of the section > That is why I am using the triple integral in terms of x, y, and z Using Simplify with assumptions {h >= 0, r2 > 0} yields your desired answer. Rob Pratt Department of Operations Research The University of North Carolina at Chapel Hill rpratt at email.unc.edu http://www.unc.edu/~rpratt/