What is the distribution for various polyhedral dice all rolled at once? Take the 5 Platonic solids from a set of Dungeons&Dragons dice.  These consist of a 4-sided, 6-sided (conventional), 8-sided, 12-sided, and 20-sided dice.  All start at the number 1 and count upwards by 1 to their total.
Roll them all, and take their sum (minimum sum is 5, max is 50).  Do so multiple times.  What is the distribution?
Obviously they will tend towards the low end [Edit: this is wrong.], since there are more lower numbers than higher.  But will there be notable inflection points at each boundary of the individual die?
[Edit:  Apparently, what seemed obvious isn't.  According to one of the commentators, the average is (5+50)/2=27.5.  I wasn't expecting this.  I'd still like to see a graph.]
 A: So I made this code:
d4 <- 1:4  #the faces on a d4
d6 <- 1:6  #the faces on a d6
d8 <- 1:8  #the faces on a d8
d10 <- 1:10 #the faces on a d10 (not used)
d12 <- 1:12 #the faces on a d12
d20 <- 1:20 #the faces on a d20

N <- 2000000  #run it 2 million times
mysum <- numeric(length = N)

for (i in 1:N){
     mysum[i] <- sample(d4,1)+
                 sample(d6,1)+
                 sample(d8,1)+
                 sample(d12,1)+
                 sample(d20,1)
}

#make the plot
hist(mysum,breaks = 1000,freq = FALSE,ylim=c(0,1))
grid()

The result is this plot.

It is quite Gaussian looking.  I think we (again) may have demonstrated a variation on the central limit theorem.
A: A little help to your intuition:
First, consider what happens if you add one to all the faces of one die, e.g. the d4.  So, instead of 1,2,3,4, the faces now show 2,3,4,5.
Comparing this situation to the original, it is easy to see that the total sum is now one higher than it used to be. This means that the shape of the distribution is unchanged, it is just moved one step to the side.
Now subtract the average value of each die from every side of that die.
This gives dice marked


*

*$-{3\over 2}$,$-{1\over 2}$,${1\over 2}$,${3\over 2}$

*$-{5\over 2}$,$-{3\over 2}$,$-{1\over 2}$,${1\over 2}$,${3\over 2}$,${5\over 2}$

*$-{7\over 2}$,$-{5\over 2}$,$-{3\over 2}$,$-{1\over 2}$,${1\over 2}$,${3\over 2}$,${5\over 2}$,${7\over 2}$


etc.
Now, the sum of these dice should still have the same shape as the original, only shifted downwards.  It should be clear that this sum is symmetrical around zero. Therefore the original distribution is also symmetrical.
A: Using the R software I posted earlier at https://stats.stackexchange.com/a/116913/919 for solving problems like this, you can compute the solution in one line:
(all <- d(1,4) + d(1,6) + d(1,8) + d(1,12) + d(1,20))

The output gives all 46 probabilities (not shown), which can be plotted with another line:
plot(all, xlab="Value", yaxp=c(0,1,2), main=expression(d[4]+d[6]+d[8]+d[12]+d[20]))

To this plot I have added the graph of the Normal distribution with the same variance and mean (employing a continuity correction),
curve(pnorm(x, mean(all)-1/2, sqrt(var.die(all))), add=TRUE, col="Red")


If you prefer to see the probability function, here it is:
with(all, plot(value, prob, type="h", main="Probability Function", cex.main=1))


Clearly the Normal approximation is already good, so we may continue to use it to describe the sum of many rolls of this combination.  But if you want to see it precisely computed, you may do so.  For instance, here is the sum of four trials with its Normal approximation superimposed (no continuity correction needed),
with(all+all+all+all, plot(value, prob, type="h", main="Sum of Four Trials", cex.main=1))
curve(dnorm(x, mean(all)*4, sqrt(4*var.die(all))), add=TRUE, col="Red", lwd=2)


A: The Central Limit Theorem answers your question. Though its details and its proof (and that Wikipedia article) are somewhat brain-bending, the gist of it is simple. Per Wikipedia, it states that

the sum of a number of independent and identically distributed random variables with finite variances will tend to a normal distribution as the number of variables grows.

Sketch of a proof for your case:
When you say “roll all the dice at once,” each roll of all the dice is a random variable.
Your dice have finite numbers printed on them. The sum of their values therefore has finite variance.
Every time you roll all the dice, the probability distribution of the outcome is the same. (The dice don’t change between rolls.)
If you roll the dice fairly, then every time you roll them, the outcome is independent. (Previous rolls don’t affect future rolls.)
Independent? Check. Identically distributed? Check. Finite variance? Check. Therefore the sum tends toward a normal distribution.
It wouldn’t even matter if the distribution for one roll of all dice were lopsided toward the low end. I wouldn’t matter if there were cusps in that distribution. All the summing smooths it out and makes it a symmetrical gaussian. You don’t even need to do any algebra or simulation to show it! That’s the surprising insight of the CLT.
