Added later: you also ask that the boundary of the triangulation be compatible with applying the same triangulation process to the boundaries of the original simplices. I think this construction above indeed has this property, if I understand correctly what you are asking for. To see this, it is convenient first (in terms of simplifying language) to replace the balls $\sum a_i^2 \le 1$ and $\sum b_i^2 \le 1$ above by the cubes $\max |a_i| \le 1$ and $\max |b_i| \le 1$, which we may do since the type B Coxeter complex will triangulate any full-dimensional convex body having the origin in its interior.
We are using a shuffled word such as $a_1a_2b_1a_3b_2b_3b_4$ to describe the simplex $0\le a_1 \le a_2 \le b_1 \le a_3\le b_2 \le b_3 \le b_4\le 1$. Now applying the boundary to this simplex means taking an alternating sum over all ways of turning one of these weak inequalities into an equality. If this equality is between some $a_i$ and some $b_j$, then the resulting face is an interior face of the triangulation -- because we can swap $a_i$ with $b_j$ to obtain a different shuffled word indexing a different simplex also having this face in its boundary. On the other hand, when a face results from an equality $a_i = a_{i+1}$ or $b_i = b_{i+1}$ or one from setting the first letter equal to 0 or the last letter equal to 1, then this is necessarily giving a boundary face in our triangulation of the product of simplices.
For comparison, when we take the boundary of our original simplex $\Delta^r $ (resp. $\Delta^s$), we take an alternating sum over ways of setting inequalities to equalities $a_i = a_{i+1}$ or $a_1=0$ or $a_r=1$ (resp. $b_j = b_{j+1}$ or $b_1=1$ or $b_s=1$). Applying our triangulation procedure after doing this, thus gives shuffled words possessing exactly one such equality of the type that puts you at the boundary of the original triangulation. So if I understand correctly what you are asking for, I think you do get your desired compatibility.

