Purpose – Multiphysical models are often useful for the design of electrical devices such as electrical machines. In this way, the modelling of thermal, magnetic and electrical phenomena by using an equivalent circuit approach is often used in sizing problems. The coupling of such models with other models is difficult to take into account, partly because it adds complexity to the process. The purpose of this paper is to propose an automatic modelling of thermal and magnetic aspects from an equivalent circuit approach, with its computation of gradients, using selectivity on the variables. Then, it discusses the coupling of various physical models, for the sizing by optimization algorithms. Sensibility analyses are discussed and the multiphysical approach is applied on a permanent magnet synchronous machine. Design/methodology/approach – The paper allows one to describe thermal and magnetic models by equivalent circuits. Magnetic aspects are represented by reluctance networks and thermal aspects by thermal equivalent circuits. From circuit modelling and analytical equations, models are generated, coupled and translated into computational codes (Java, C), including the computation of their Jacobians. To do so, model generators are used: CADES, Reluctool, Thermotool. The paper illustrates the modelling and automatic programming aspects with Thermotool. The generated codes are directly available for optimization algorithms. Then, the formulation of the coupling with other models is studied in the case of a multiphysical sizing by optimization of the Toyota PRIUS electrical motor. Findings – A main specificity of the approach is the ability to easily deal with the selectivity of the inputs and outputs of the generated model according to the problem specifications, thus reducing drastically the size of the Jacobian matrix and the computational complexity. Another specificity is the coupling of the models using analytical equations, possibly implicit equations. Research limitations/implications – At the present time, the multiphysical modelling is considered only for static phenomena. However, this limit is not important for numerous sizing applications. Originality/value – The analytical approach with the selectivity gives fast models, well adapted for optimization. The use of model generators allows robust programming of the models and their Jacobians. The automatic calculation of the gradients allows the use of determinist algorithms, such as sequential quadratic programming, well adapted to deal with numerous constraints. [ABSTRACT FROM AUTHOR]