Cup and wedge product in singular and de Rham cohomology
Solution 1:
The deRham map gives an arrow $\theta:H_{dR}(X) \to H(X)^*$, and the algebra structure on the right is induced from the coalgebra structue on $H(X)$ given by the Alexander-Whitney coproduct: $\Delta(\sigma) = \sum \sigma_i \otimes \sigma^i$ where $\sigma_i$ is the front $i$-face and $\sigma^i$ is the back $i$-face. Denote by $\langle -,-\rangle : H_{dR}(X)\otimes H(X)^*\to \mathbb R$ pairing that gives rise to $\theta$.
In terms of the pairing, you are correct this is the same as stating that, for any two forms $\omega,\eta$ of degee $p,q$ and any singular chain $\sigma$ of degree $p+q$, we have that $$\tag{1}\langle \omega\wedge \eta,\sigma\rangle = \langle \omega,\sigma_p\rangle\langle \eta,\sigma^p\rangle$$ which is the same as $\langle \mu,-\rangle =\mu_\mathbb R\langle - \otimes -,\Delta\rangle$. You also note that there is a problem of writing simplices as products. There is a way of fixing this problem, which is (as Spivak does and in another context, Serre does), for example, done by working with singular cubes instead of singular chains.
Concretely, compute $H_*(X)$ using maps $\sigma : [0,1]^*\to X$. The differential is another, but the coproduct is simply given by restricting to coordinates, as in the case of simplices (though more terms appear). Concretely, if $\sigma : [0,1]^p\to M$ is a singular cube, the coproduct of $\sigma$ is given by a sum
$$\sum_{A\sqcup B = [p]}\sigma_A^0\otimes \sigma_B^1 $$ where $(-)_S^\varepsilon$ means you set the $S$-coordinates of the subset $S$ of $[p] = \{1,\ldots,p\}$ to $\varepsilon$ and keep the coordinates of the complement. See page 17 in Serre's thesis. In this case, if $\sigma$ is a singular cube and $\omega = fdx_I$ is a form of the same degree $p$, it is straightforward to check that
$$\langle \omega,\sigma\rangle = \int_{[0,1]^p} f(\sigma) J_\sigma dx_1\cdots dx_p$$
and then $(1)$ is an immediate consequence of Fubini's theorem, as you wanted. There is a natural isomorphism between $H(X)$ computed in two different ways, either from cubical or simplicial singular chains, coming from a map from the cubical complex of $X$ to the singular complex of $X$ that needs to respect the coproduct (this is the nontrivial part of you question, I think, but you can find this map say in Serre's thesis), and this completes the proof that the deRham isomoprhism is a map of algebras.
Solution 2:
I just wanted to point out that $$\int_{\sigma^{q+p}}\omega^q\wedge \eta^p \neq \int_{\sigma|_{[e_0,\dots,e_q]}}\omega^q \cdot \int_{\sigma|_{[e_{q},\dots,e_{q+p}]}}\eta^p$$ in general, since simplices are not closed. The formula $$\omega^q\cup \eta^p (\sigma^{p+q}) = \omega(\sigma|_{[e_0,\dots,e_q]})\cdot \eta(\sigma|_{[e_{q},\dots,e_{q+p}]})$$ is valid only for singular simplex $\sigma:\Delta^{p+q}\to M$.