This paper presents the improved methodology for the direct calculation of steady-state periodic solutions for electromagnetic devices, as described by nonlinear differential equations, in the time domain. A novel differential operator is developed for periodic functions and the iterative algorithm determining periodic steady-state solutions in a selected set of time instants is identified. Its application to steady-state analysis is verified by an elementary example. The modified algorithm reduces the complexity of steady-state analysis, particularly for electromagnetic devices described by high-dimensional nonlinear differential equations.

This paper describes an algorithm for finding steady states in AC machines for the cases of their two-periodic nature. The algorithm enables to specify the steady-state solution identified directly in time domain despite of the fact that two-periodic waveforms are not repeated in any finite time interval. The basis for such an algorithm is a discrete differential operator that specifies the temporary values of the derivative of the twoperiodic function in the selected set of points on the basis of the values of that function in the same set of points. It allows to develop algebraic equations defining the steady state solution reached in a chosen point set for the nonlinear differential equations describing the AC machines when electrical and mechanical equations should be solved together. That set of those values allows determining the steady state solution at any time instant up to infinity. The algorithm described in this paper is competitive with respect to the one known in literature an approach based on the harmonic balance method operated in frequency domain.

This paper presents the concept of using algorithms for reducing the dimensions of finite-difference equations of two-dimensional (2D) problems, for second-order partial differential equations. Solutions are predicted as two-variable functions over the rectangular domain, which are periodic with respect to each variable and which repeat outside the domain. Novel finite-difference operators, of both the first and second orders, are developed for such functions. These operators relate the value of derivatives at each point to the values of the function at all points distributed uniformly over the function domain. A specific feature of the novel operators follows from the arrangement of the function values as well as the values of derivatives, which are rectangular matrices instead of vectors. This significantly reduces the dimensions of the finite-difference operators to the numbers of points in each direction of the 2D area. The finite-difference equations are created exemplary elliptic equations. An original iterative algorithm is proposed for reducing the process of solving finite-difference equations to the multiplication of matrices.