In a ball with random thread/strings, how does the density of threads/strings change with radius? A large plastic ball full of holes is given. (So the holes are in a plastic shell.) Straight threads connect these holes randomly, by passing through the interior of the ball/shell.
For a big ball or shell, say a meter in size, with thousands of holes, this makes (1/2 times) thousands of straight threads inside it. (Each hole has the diameter of the thread, so that each hole can only have one string passing through it.)
Now the question: Inside the ball/shell (assumed to be large), is the density of the random threads homogeneous, or does it depend on the radius?
 A: Added. Now, I think that G. Smith gave the right answer to the initial question. And I was solving a different problem. My former solution implies that we chose any thread with equal probability and then uniformly chose the point of this thread. This procedure is not equivalent to the finding of mass distribution. To find the distribution of mass we should choose threads with probabilities proportional to their length. Just because of the lengthy thread contains more mass. Doing like this one obtains a mass distribution with constant density.
Former solution. I obtained the following expression for the density of "matter" inside the sphere of radius 1
$$
\rho(r) = \frac{A}{r}\log\left(\frac{\sqrt{2}(1+r)}{\sqrt{|2r^2+\cos(\varphi(r))-1|}+\sqrt{2}\cos(\varphi(r)/2)}\right).
$$
Here $A$ is constant and $\varphi(r) = 2\arcsin(r)$. The valiue of this density at $r=0$ is equal to $A$, and it diverges as $r$ tends to 1.
Upd. This expression is obtained in the following way.
For any pair of holes let's draw z-axis through one of them and the center of sphere. Then a position of the second is defined by a polar angle $\varphi\in[0,\pi]$. The angle is random 
and the corresponding pdf is $w_1(\varphi)=\sin(\varphi)/2$. Uniform distribution of "matter" along the line connecting two holes leads to the following distribution of radius:
$$
w_2(r|\varphi) = \frac{r}{\cos(\varphi/2)\sqrt{r^2-\sin^2(\varphi/2)}},
$$
where $r\in[\sin(\varphi/2),1]$. The minimal value of radius along the line is equal to $\sin(\varphi/2)$, hence the definition of $\varphi(r)$. Averaging with respect to angles gives the radius pdf:
$$
w_3(r) = \int_0^{\varphi(r)} w_1(\varphi)w_2(r|\varphi)d\varphi .
$$
And the density of "matter" is proportional to $w_3(r)/r^2$.
A: I believe the density is homogeneous throughout the ball.
I did a numerical simulation of this in Mathematica. I assumed the sphere had radius 1 and generated 100,000 pairs of random points on it, each pair to be connected with string. Then I analyzed this set of random strings to see how much total mass (i.e., length of string) lay between $r$ and $r+dr$ in various spherical shells with radii (0.1, 0.2, …, 0.9), using a bit of geometry. Then I divided by the square of the radius of the spherical shell to get the volume density and plotted it. The 9 points lay almost on a horizontal line:

The horizonal axis is the radial coordinate and the vertical axis is the mass density.
ADDENDUM
Here is an analytic proof that the density is homogeneous, based on @Gec's answer. I agree with his approach but not his former result.
Take the sphere to have unit radius and the strings to have unit linear mass density so that the mass of a small segment is just the length of that segment.
As Gec points out, a string can be characterized by the angle it subtends, which I'm going to call $\theta$. A string has a minimum radial distance of $\cos{(\theta/2)}\equiv a$ and a length of $2\sin{(\theta/2)}=2\sqrt{1-a^2}$.
Introduce a linear coordinate $s$ along the string, measured from its midpoint. Then one has $a^2+s^2=r^2$ so
$$s=\sqrt{r^2-\cos^2{(\theta/2)}}.$$
Differentiating with respect to $r$, we find
$$ds=\frac{r\,dr}{\sqrt{r^2-\cos^2{(\theta/2)}}}=\frac{r\,dr}{\sqrt{r^2-a^2}}.$$
This tells us that the mass of this string that lies within a spherical shell between $r$ and $r+dr$ is
$$dm=\frac{2r\,dr}{\sqrt{r^2-\cos^2{(\theta/2)}}}=\frac{2r\,dr}{\sqrt{r^2-a^2}}.$$
(The string passes through the shell on both sides of its center.)
We can check that this is correct by integrating it over $r$ from $a$ to $1$:
$$m=2\int_a^1\frac{r\,dr}{\sqrt{r^2-a^2}}=2\sqrt{1-a^2},$$
which agrees with the length of the string.
Now we need to integrate over strings between random points on a sphere.
As Gec pointed out, the spherical symmetry means that we can consider just strings with one endpoint at the north pole, and the other end at polar angle $\theta$ and azimuthal angle $\phi$. To randomly average a quantity $f$ over the randomly placed other end, we compute $\langle f \rangle=\frac{1}{4\pi}\iint f\,\sin{\theta}\,d\theta\,d\phi$. By azimuthal symmetry, this simplifies to $\frac{1}{2}\int f\,\sin{\theta}\,d\theta$.
To compute the averaged mass $dM$ in a spherical shell between $r$ and $r+dr$, we integrate $dm$ over $\theta$, but only between $2\cos^{-1}r$ and $\pi$. For smaller angles, the string would not pass through the shell and thus would not contribute any mass. So
$$\frac{dM}{dr}=\int_{2\cos^{-1}r}^\pi \frac{r\sin{\theta}\,d\theta}{\sqrt{r^2-\cos^2{(\theta/2)}}}$$
The substitution $u=\cos{(\theta/2)}$ simplifies this integral to
$$\frac{dM}{dr}=4r\int_0^r\frac{u\,du}{\sqrt{r^2-u^2}}=4r^2.$$.
To get the volume mass density $\rho=dM/dV$, we divide by the area of the spherical shell, $4\pi r^2$, to get a homogeneous density of
$$\rho=\frac{1}{\pi}.$$
My numerical simulation gave $2$ rather than $1/\pi$ because (1) I didn’t multiply by 2 to take into account that a string passes through a shell on both sides of its midpoint, and (2) at the end I divided by $r^2$ rather than $4\pi r^2$.
A: While the number of strings passing nearby the shell is higher than that of the strings passing through the center, also the constant radius surfaces near the shell are larger than those near the center of the sphere.
We can define the string density $\rho $ through
$$4\pi r^2\rho(r)= N (r) $$
where $N (r) $ is the number of times that the strings intersect the constant radius $r$ surface.
Note that, assuming the holes follow a homogeneous distribution, in the limit of large number of holes you are just connecting random points of the sphere with lines that cross the sphere.
Fixed a point from which the line is drawn, you have equal probability of connecting it to any other point of the sphere.
The line (string) will pass through the center only if the opposite point is chosen.
Conversely every line will pass through the sphere surface and almost every line will pass at a slightly smaller radius. 
You can compute the number of lines of a certain length $L $ that can be drawn from a chosen point; even better you can express this using the angle formed by the two points connected and the center of the sphere:
$$L=2Rsin\theta\;,\qquad
N_L= 2\pi R sin 2\theta$$
The segments corresponding to an angle $\theta$ will contribute to the density for radiai in the range $[Rcos\theta, R]$ with 2 points each except in the case of the minimum radius value (here the string passes only once).
Now $N (r) $ will be proportional to
$$ \int_{\theta*}^{\pi/2} 2\pi R sin (2\theta)d\theta $$
where $cos\theta*=\frac {r}{R}$.
The proportionality constant is basically the number of the endpoints since you integrate their distribution on the $R $ shell (you also have a factor 2 because each string is counted twice almost everywhere and a factor 1/2 to avoid the overcounting when integrating over endpoints).
The integral gives $$2\pi R \left(\frac {r}{R}\right)^2$$
so that when you compute $\rho (r)$ you indeed get a term which is independent from $r $.
If we were to stop here, the density would be uniform.
One could think that we still need to remove the overcounting of pieces of string in the minimum radius each segment reaches:
do we have to subtract from $N (r)$ one counting of the intersection at minimum radius, i.e. the quantity $$2\pi R sin (2\theta*)$$
This would give a part which is dependent on $r $ in the distribution: $$\rho (r)\sim const+\frac {\sqrt{1-(r/R)^2}}{r}$$
The truth is that the term must be subtracted into the integral of $N (r) $ and there gives zero contribute, since its a modification on a set of zero measure.
So in conclusion there is no term to be subtracted and the density is indeed constant.
It would be nice to see if there are other endpoints distributions that are mimicked by the string density..
