$\newcommand\R{\mathbb R}$Now concerning any dimension $m$:
Without loss of generality, we have two open balls: $B_1=B_0(r_1)$ and $B_2=B_D(r_2)$, for some $D\in\R^m$ with $\|D\|=d$, where $B_x(r)$ is the open ball in $\R^m$ of radius $r$ centered at $x$ and $\|\cdot\|$ is the Euclidean norm.
Take any unit vector $u$. Take any $x\in B_1$ such that $x\cdot u=0$, where $\cdot$ denotes the dot product. Then the condition that the line through $x$ in the direction of $u$ intersects $B_2$ is that $x\in B_D(R_u)$, where $R_u:=\sqrt{(D\cdot u)^2+r_2^2}$; this can be obtained by minimizing $\|x-tu-D\|^2$ in real $t$ and taking into account the condition $x\cdot u=0$.
So, $$\mathcal A(C,u)=L_{m-1}(B_0(r_1)\cap B_D(R_u)\cap H_u),$$ where $H_u:=\{x\in\R^m\colon\, x\cdot u=0\}$ and $L_{m-1}$ is the $(m-1)$-dimensional Lebesgue measure over the hyperplane $H_u$. So, $\mathcal A(C,u)$ is the $(m-1)$-volume of the intersection of two balls, $B_0(r_1)\cap H_u$ and $B_D(R_u)\cap H_u$, in $H_u$.
So, $\mathcal A(C,u)$ is the $(m-1)$-volume of the intersection of two balls in $\R^{m-1}$, of radii $r_1$ and $r_{2,u}:=\sqrt{\max[0,r_2^2-(D\cdot u)^2]}$, with the centers of the two balls at distance $d_u:=\sqrt{d^2-(D\cdot u)^2}$ from each other.
So, $\mathcal A(C,u)$ is an ordinary integral of the ($(m-2)$-dimensional) volumes of $(m-2)$-dimensional balls, which is actually known -- but is expressed via the incomplete beta function, rather than in elementary functions.
Now it remains to integrate $\mathcal A(C,u)$ over all $u$ on the $(m-1)$-dimensional unit sphere $S^{m-1}$. When doing that, by the spherical symmetry we can replace $D\cdot u$ in the above expressions for $r_{2,u}$ and $d_u$ by $u_1d$, where $u_1$ is the first coordinate of the unit vector $u$. So, $$\mathcal A(C,u)=g(u_1^2),$$ for a certain function $g=g_{r_1,r_2,d}$, expressed explicitly via an incomplete beta function. So, the étendue measure of $C$ is $$\frac12\,\int_{S^{m-1}}\mathcal A(C,u)\,du= \frac12\,\int_{S^{m-1}}g(u_1^2)\,du \\ =\frac12\,Eg(U_1^2)=\frac{m-1}2\,\int_0^1 g_{r_1,r_2,d}(t)(1-t)^{m-2}\,dt,$$ where $U_1^2$ is a (real-valued) random variable with the beta distribution with parameters $1,m-1$.