This paper presents a comprehensive approach to design of series of permanent magnet assisted synchronous reluctance motors using combined analytical and finite element calculations. A global optimization metaheuristic algorithm (Differential Evolution) is utilized in order to achieve optimal design in terms of maximum torque per volume with numerous specific boundaries imposed on motor geometry and performance. A novel approach to calculation of the specified constant power speed range and demagnetization effect in sudden symmetrical short circuit using iterative finite element magnetostatic simulations is presented. Based on the results of optimized design, a 100 kW prototype was built and tested. 7. β/β 0-angle of the slanted magnet relative to the maximum allowed angle of the slanted magnet (β 0 = 0.5π(1 − 1/p)), 8. λp-angular span of the inner cavity relative to the pole pitch. Lower and upper boundary constraints of the variables are: 0.45 ≤Ds/Dout ≤ 0.75, 0.2 ≤hys/[(Dout − Ds)/2] ≤ 0.6, 0.3 ≤bt/τu ≤ 0.7, 0.05 ≤λm ≤ 0.5, 0.2 ≤λ md1 ≤ 0.6, 0.05 ≤λ md2 ≤ 0.4, 0.5 ≤β/β 0 ≤ 1.0, 0.75 ≤λp ≤ 0.95. Variables with constant value, i.e. design parameters that do not change during the optimization procedure (preset parameters) according to Fig. 1 are: 1. δ = 0.8 mm-air-gap length,