In contrast with conventional PCA, a direct superposition and joint interpretation of loading plots is not possible in three-way PCA, since there may be data variance which is described by unequal components of different modes. The contributions to variance of all possible combinations of components are described in the core matrix. Body diagonalization, which is achieved by appropriate rotation of component matrices, is an essential tool for simplifying the core matrix structure. The maximum degree of body diagonality which may be obtained from such transformations is analysed from both the mathematical and simulation viewpoints. It is shown that, at least in the average case, high degrees can be expected, which makes the procedure reasonable for many practical applications. Furthermore, simulation as well as theoretical derivation show that the success of body diagonality depends on the socalled polarity of the core array. The methodology is illustrated by a three-way data example from environmental chemistry.