Consistency Results Separating Three Cardinal Characteristics Simultaneously
Bumping an old question, because the situation has changed dramatically in the last year or so:
The most well-known cardinal characteristics of the continuum are those appearing in Cichoń's Diagram, which also presents all $ZFC$-provable relations between pairs of these characteristics: besides the dominating number $\mathfrak{d}$ and the bounding number $\mathfrak{b}$, these are all of the form $inv(\mathcal{I})$, for $\mathcal{I}$ the ideal of meager sets $(\mathcal{M})$ or null sets $(\mathcal{N})$, and $inv$ one of $add, cov, non, cof$:
$add$ of an ideal $\mathcal{I}$ is the least cardinality of a family of sets from $\mathcal{I}$ whose union is not in $\mathcal{I}$;
$cov$ of an ideal $\mathcal{I}$ is the least cardinality of a family of sets from $\mathcal{I}$ covering all of $\mathbb{R}$;
$non$ of an ideal $\mathcal{I}$ is the least cardinality of a set of reals not in $\mathcal{I}$; and
$cof$ of an ideal $\mathcal{I}$ is the cofinality of the partial order $\langle\mathcal{I}, \subseteq\rangle$.
As far as I know, the first results separating more than two cardinal characteristics from Cichon's diagram came out in the last two years:
In late 2012, Mejia constructed several examples of models separating multiple characteristics at once (see section 6).
In a paper arxived today(!), Fischer/Goldstern/Kellner/Shelah produced a model of $\aleph_1=\mathfrak{d}=cov(\mathcal{N})<non(\mathcal{M})<non(\mathcal{N})<cof(\mathcal{N})<2^{\aleph_0}$.
Now, I don't understand the proofs at all, but it seems the proofs by Mejia and by F/G/K/S are fundamentally different. The F/G/K/S paper has this to say about the two different approaches to separating multiple characteristics simultaneously:
"We cannot use the two best understood methods [to separate $\ge 3$ characteristics simultaneously], countable support iterations of proper forcings (as it forces $2^{\aleph_0}\le\aleph_2$) and, at least for the "right hand side" of the diagram, we cannot use finite support iterations of ccc forcings in the straightforward way (as it adds lots of Cohen reals, and thus increases $cov(\mathcal{M})$ to $2^{\aleph_0}$).
There are ways to overcome this obstacle. One way would be to first increase the continuum in a "long" finite support iteration, resulting in $cov(\mathcal{M})=2^{\aleph_0}$, and then "collapsing" $cov(\mathcal{M})$ in another, "short" finite support iteration. In a much more sophisticated version of this idea, Mejia [Mej13] recently constructed several models with many simultaneously different cardinal characteristics in Cichon's Diagram (building on work of Brendle [Bre91], Blass-Shelah [BS89] and Brendle-Fischer [BF11]).
We take a different approach, completely avoiding finite support, and use something in between a countable and finite support product (or: a form of iteration with very "restricted memory").
This construction avoids Cohen reals, in fact it is $\omega^\omega$-bounding, resulting in $\mathfrak{d}=\aleph_1$. This way we get an independence result "orthogonal" to the ccc/finite-support results of Mejia.
The fact that our construction is $\omega^\omega$-bounding is not incidental, but rather a necessary consequence of the two features which, in our construction, are needed to guarantee properness: a "compact" or "finite splitting" version of pure decision, and fusion . . ."
[Creature forcing and five cardinal characteristics in Cichon's diagram], pg. 2
I can't resist bumping this to add one fascinating result (from this morning!) of Goldstern, Kellner, and Shelah: that as many as possible of the cardinal characteristics in Cichon's diagram can be separated simultaneously! Now granted, this isn't a consistency result relative to ZFC - at present, it takes four strongly compact cardinals, which is no small assumption - but still, this is absolutely amazing.
EDIT: And as of a couple weeks ago, the other shoe has dropped (thanks to Jakob for bringing this to my attention below): Goldstern, Kellner, Mejia, and Shelah got it in ZFC.
There is a paper of Shelah and Goldstern devoted to the separation of many simple cardinal invariants (this is a technical term). There are more recent papers on this subject by Kellner and Shelah, if I remember correctly.
An easy case that I am very familiar with are the so called localization numbers.
A closed set $S\subseteq\omega^\omega$ is $n$-ary if in the tree $T(S)$ of finite initial segments of elements of $S$ every node has at most $n$ immediate successors.
For $n\geq 2$ let $\ell_n$ be the least size of a family of $(n-1)$-ary sets that covers all of $n^\omega$.
Any finitely many $\ell_n$ can be separated from each other simultaneously.
This is shown in [Geschke, Kojman, Convexity numbers of closed subsets in R^n,
Proc. Am. Math. Soc. 130, No. 10, 2871-2881 (2002)], which is here (Wayback Machine).
Proofs of such statements usually involve forcing with large countable support products over a model of GCH rather than iterated forcing. However, there are also some examples that use iterated forcing. See for example this paper by Shelah and Steprans.