Volumen de la intersección de dos cilindros

Question

Tengo dos cilindros infinitos de radio unitario en $\mathbb{R}^3$ cuyos ejes son líneas oblicuas. Digamos que el eje de una está centrado en la $x$ -eje, y el eje del otro está determinado por los dos puntos $a$ y $b$ . ¿Existe una fórmula para el volumen de su intersección, en función de $a$ y $b$ ?

---

         

---

Answer 1

Dejemos que

• $r = 1$ sean los radios comunes de los dos cilindros.
• $2d$ sea la distancia más cercana entre los dos ejes de los cilindros.
• $\alpha = 2\beta$ sea el ángulo entre los dos ejes de los cilindros.

Elige el sistema de coordenadas para que el eje de los cilindros pase por $(0, \pm d, 0)$ con vectores tangentes $( \pm\sin\beta, 0, \cos\beta)$ respectivamente. No es difícil ver que los dos cilindros se dan por.

$$\mathcal{C}_{\pm} = \left\{ (x,y,z) \in \mathbb{R}^3 : ( x\cos\beta \mp z\sin\beta )^2 + (y \mp d)^2 \le r^2 \right\}$$

Si se cruza $\mathcal{C}_{\pm}$ con un plano de constante $z$ las intersecciones serán elipses de dos ejes alineados con semieje mayor $\frac{r}{\cos\beta}$ en el $x$ -dirección y semieje menor $r$ en el $y$ -dirección. Introducir nuevas coordenadas $$(u,v,w) = (x\cos\beta, y, z\sin\beta) \quad\iff\quad (x,y,z) = \left(\frac{u}{\cos\beta}, v, \frac{w}{\sin\beta}\right) $$ Las intersecciones de $\mathcal{C}_{\pm}$ con un plano de constante $w$ se convierten en dos círculos de radios $r$ . $$(u \mp w)^2 + (v \mp d)^2 \le r^2$$ y estos dos círculos se cruzan cuando y sólo cuando $d^2 + w^2 \le r^2$ .

Cuando se cruzan, que $\theta = \cos^{-1}\left(\frac{\sqrt{d^2 + w^2}}{r}\right)$ . En el diagrama siguiente, la intersección es la delimitada por los dos arcos circulares en rojo. Su área es igual a

$$\begin{align} & \verb/Area/(\text{sector } ABD) + \verb/Area/(\text{sector } CDB) - \verb/Area/(\text{rhombus } ABCD)\\ = & 2\verb/Area/(\text{sector } ABD) - 4\verb/Area/(\text{triangle } OAB)\\ = & 2\left(\frac12 r^2(2\theta)\right) - 4\left(\frac12 (r\sin\theta)(r\cos\theta)\right)\\ = & r^2(2\theta - \sin(2\theta)) \end{align} $$

$\hspace1.25in$

Dejemos que $\displaystyle\;\theta_0 = \cos^{-1}\left(\frac{d}{r}\right)\;$ , el volumen que queremos se convierte en

$$\begin{align}\verb/Volume/ &= \overbrace{\frac{1}{\cos\beta\sin\beta}}^{\color{blue}{\text{Jacobian}}} \int_{-\sqrt{r^2 - d^2}}^{\sqrt{r^2 - d^2}} r^2 (2\theta - \sin(2\theta)) d w\\ &= \frac{-4r^3}{\sin\alpha} \int_{0}^{\theta_0} (2\theta - \sin(2\theta)) d\sqrt{\cos^2\theta - \cos^2\theta_0}\\ \color{blue}{\text{ integrate by part } \rightarrow} &= \frac{16 r^3}{\sin\alpha}\int_0^{\theta_0} \sin^2\theta \sqrt{\cos^2\theta - \cos^2\theta_0} d\theta\\ \color{blue}{ m = \sin^2\theta_0 \rightarrow} &= \frac{16 r^3}{\sin\alpha}\int_0^\sqrt{m} \frac{\sin^2\theta}{\cos\theta}\sqrt{m - \sin\theta^2} d\sin\theta\\ \color{blue}{ \sin\theta = \sqrt{m} t \rightarrow} &= \frac{16 r^3 m^2}{\sin\alpha}\int_0^1 t^2\sqrt{\frac{1-t^2}{1-mt^2}} d t\\ &= \frac{16 r^3 m^2}{\sin\alpha}\int_0^1 \frac{t^2 - t^4}{\sqrt{(1-t^2)(1-mt^2)}} dt \end{align} $$

Dejemos que $Q(t) = (1-t^2)(1-mt^2)$ es fácil de comprobar

$$\frac{d}{dt}\left(t \sqrt{Q(t)}\right) = \frac{1}{\sqrt{Q(t)}}\left(3m t^4 - 2(m+1) t^2 + 1\right)$$ Podemos usar esto para transformar la última integral como

$$\begin{align}\verb/Volume/ &= \frac{16 r^3 m}{3\sin\alpha}\int_0^1 \frac{3m t^2 - 2(m+1) t^2 + 1}{\sqrt{Q(t)}} dt\\ &= \frac{16 r^3}{3\sin\alpha}\int_0^1 \frac{(2-m)(1-mt^2) - 2(1-m)}{\sqrt{Q(t)}} dt \end{align} $$ Compárelo con la definición del completo integrales elípticas de la $1^{st}$ y $2^{nd}$ tipo.

$$\begin{align} K(m) &= \int_0^{\pi/2} \frac{d\theta}{\sqrt{1-m\sin^2\theta}} = \int_0^1 \frac{dt}{\sqrt{(1-mt^2)(1-t^2)}}\\ E(m) &= \int_0^{\pi/2} \sqrt{1-m\sin^2\theta} d\theta = \int_0^1 \sqrt{\frac{1-mt^2}{1-t^2}} dt \end{align} $$ Obtenemos

$$\verb/Volume/ = \frac{16 r^3}{3\sin\alpha}\left((2-m)E(m) - 2(1-m)K(m)\right)$$

Como comprobación de cordura, cuando $d \to 0$ , $m \to 1^{-}$ tenemos

$$(1-m)K(m) \to 0\quad\text{ and }\quad (2-m)E(m) \to E(1) = 1$$ Esto conduce a un resultado conocido para el volumen de intersección cuando los ejes se cruzan.

$$\verb/Volume/_{d = 0} = \frac{16r^3}{3\sin\alpha}$$