Volume of the hemisphere cut at an angle

Volume of hemisphere asymmetrically cut by a flat plane under a given angle.

Cette page existe grâce aux efforts des personnes suivantes :

Anton

Créé: 2018-07-18 21:21:18, Dernière mise à jour: 2020-11-03 14:19:35

The calculator below gives volume of a body formed by a hemisphere cut by a flat surface under a given angle at a given distance from the hemisphere edge.

1)Hemisphere cut at a given length under a given angle
1)Hemisphere cut at a given length under a given angle



The cut angle is between (0..180) degrees, the height must not exceed the hemisphere diameter. You may find calculation formulas just below the calculator.

PLANETCALC, Volume of hemisphere cut by a plane under a given angle

Volume of hemisphere cut by a plane under a given angle

The distance between hemisphere base edge and the cutting plane.
Angle in degrees between hemisphere base and the cutting plane.
Digits after the decimal point: 5
Volume
 

Formulas to calculate hemisphere cut volume

Let's start with simple cases.

Hemisphere cut over center h=R

2)Hemisphere center cut.
2)Hemisphere center cut.

This is simplest case - we just need slightly modify hemisphere volume formula:
V=\frac{2}{3}R^3\pi
In this way:
V=\frac{2}{3}R^3\alpha, where a is the cut angle in radians.

Hemisphere cut under 90 degrees

3)90 degrees hemisphere cut.
3)90 degrees hemisphere cut.



If the hemisphere is cut under right angle, we simply take half spherical cap volume.

V=\frac{\pi H^2(3R-H)}{6}

Hemisphere cut below center h > R

4)Hemisphere cut below the hemisphere center.
4)Hemisphere cut below the hemisphere center.



We observe the hardest case when the hemisphere is cut below the center, so remained base edge exceeds the hemisphere radius. In this case we calculate discarded part volume and subtract it from the hemisphere volume to obtain remained part volume:
V=\frac{2}{3}R^3\pi-V_d
Vd is calculated by the formulas below, assuming \alpha=\pi-\alpha

Hemisphere cut above center h < R

5)Hemisphere cut above center
5)Hemisphere cut above center



This general case we'll solve by volume integral.
To simplify integration limits we place our hemisphere cut parallel to Z-O-Y plane, so the hemisphere center point is at 0 and the cutting plane is at sin(α)*r0 distance from the Z-O-Y plane (where r0=R-h is a distance from the cutting plane to the hemisphere center):

6)Hemisphere cut positioning in the coordinate system
6)Hemisphere cut positioning in the coordinate system



Cutting plane equation in linear coordinate system is x=r_0 \sin \alpha.
In spherical (ISO) coordinate system: r=\frac{r_0 \sin \alpha}{\sin \theta \cos \varphi}(1)
Now we solve volume integral for our hemisphere cut in spherical coordinate system:
V=2\int_{\frac{\pi}{2}-\alpha}^{\frac{\pi}{2}-{\varphi}_0}d{\varphi}\int_{\theta_{\varphi}}^{\frac{\pi}{2}}\sin{\theta}d{\theta} \int_{r_{\varphi,\theta}}^R{r^2}d{r} (2)
We integrate over one hemisphere half, and multiply volume integral by two, since we put the cut hemisphere symmetrically against to the X-O-Y plane.

Integration intervals

  1. φ angle interval is constant:[\frac{\pi}{2}-\alpha-\frac{\pi}{2}-{\varphi}_0]
    Where a is the cut angle in radians
    and {\varphi}_0=\arcsin(\frac{r_0 \sin \alpha}{R}) (3)
  2. Radius integration interval
    The radius upper bound is limited by sphere surface - R. Cutting plane limits the lower radius bound, according (1)
    r_{\varphi,\theta}=\frac{r_0 \sin \alpha}{\sin \theta \cos \varphi}(4)
  3. θ angle intervals
    We'll integrate only over upper half of θ interval, so upper limit of θ is \frac{\pi}{2} and lower limit equation can be evaluated from (1) with r=R in this way:
    \theta_{\varphi}=\arcsin\left(\frac{r_0 \sin{\alpha}}{R\cos\varphi}\right)(5)
    Or using (3) we can write it like this:
    \theta_{\varphi}=\arcsin\left(\frac{\sin {\varphi}_{0}}{\cos\varphi}\right)(6)

Volume formula

To obtain Volume formula for hemisphere cut we solve integral (2) using intervals (3),(4),(6):

  1. Integrate over r

    V=2\int_{\frac{\pi}{2}-\alpha}^{\frac{\pi}{2}-{\varphi}_0}d{\varphi}\int_{\theta_{\varphi}}^{\frac{\pi}{2}}\sin{\theta}d{\theta}\left.\frac{1}{3}r^3\right|_{\frac{r_0 \sin \alpha}{\sin \theta \cos \varphi}}^R

    Result:
    V=\frac{2}{3}\int_{\frac{\pi}{2}-\alpha}^{\frac{\pi}{2}-{\varphi}_0}d{\varphi}\int_{\theta_{\varphi}}^{\frac{\pi}{2}}R^3\sin{\theta}-\frac{r_0^3 \sin^3 \alpha}{\cos^3 \varphi}d{\theta}

  2. Integrate over θ
    V=\frac{2}{3}\int_{\frac{\pi}{2}-\alpha}^{\frac{\pi}{2}-{\varphi}_0}d{\varphi}\left. \cot \theta  \frac{r_0^3\sin^3 \alpha}{\cos^3 \varphi} - R^3 \cos\theta\right|_{\arcsin\left( \frac{\sin {\varphi}_{0}}{\cos \varphi}\right)}^{\frac{\pi}{2}}
    Result:
    V=\frac{2}{3}\int_{\frac{\pi}{2}-\alpha}^{\frac{\pi}{2}-{\varphi}_0} R^2 r_0\sin\alpha \frac{\sqrt{\frac{\cos^2\varphi}{\sin^2{\varphi}_0}-1}}{\cos\varphi} - r_0^3\sin^3\alpha\frac{\sqrt{\frac{\cos^2\varphi}{\sin^2{\varphi}_0}-1}}{\cos^3\varphi} d{\varphi}

  3. Integrate over φ
    V=\frac{2}{3}R^2 r_0 \sin \alpha \left( \csc{\varphi}_0 \arcsin\frac{\sin\varphi}{\cos{\varphi}_0} - \arctan\frac{\sin\varphi \sin {\varphi}_0}{\sqrt{\cos^2\varphi-\sin^2\varphi_0}} \right)
     - \frac{r_0^3\sin^3 \alpha}{3}\left( \cot^2{\varphi}_0  \arctan\frac{\sin\varphi \sin {\varphi}_0}{\sqrt{\cos^2\varphi-\sin^2\varphi_0}} +\csc {\varphi}_0 \tan\varphi \sec\varphi \sqrt{\cos^2\varphi-\sin^2\varphi_0} \right)\Bigg|_{\frac{\pi}{2}-\alpha}^{\frac{\pi}{2}-{\varphi}_0}
    So the final volume formula is:
    V= \frac{r_0^3\sin^3 \alpha}{3}\left( \cot^2{\varphi}_{0}(K_2- \frac{\pi}{2}) +K_1\csc {\varphi}_0 \cot\alpha \csc\alpha\right)
     -\frac{2}{3}R^2 r_0 \sin \alpha \left( \csc{\varphi}_0 \left(\arcsin\left[\frac{\cos\alpha}{\cos{\varphi}_0}\right]-\frac{\pi}{2}\right) - K_2+\frac{\pi}{2}\right)
    where:
    K_1=\sqrt{\sin^2\alpha-\sin^2\varphi_0},
    K_2=\arctan\frac{\cos\alpha \sin {\varphi}_0}{K_1}.
URL copiée dans le presse-papiers
PLANETCALC, Volume of the hemisphere cut at an angle

commentaires