The aim of the article is to demonstrate an algorithm for constructing measure-preserving three-dimensional chaotic maps defined in domains formed by rotation bodies. On the one hand, the class of multidimensional chaotic mappings is expended, and on the other hand, we obtain formulas for simulating pseudo-random quantities that are in demand in problems solving by the Monte Carlo method.