Next Article in Journal
Features and Practicability of the Next-Generation Sensors and Monitors for Exposure Assessment to Airborne Pollutants: A Systematic Review
Next Article in Special Issue
Fostering Cyber-Physical Social Systems through an Ontological Approach to Personality Classification Based on Social Media Posts
Previous Article in Journal
IoT-Enabled Smart Cities: Evolution and Outlook
Previous Article in Special Issue
A Coalitional Distributed Model Predictive Control Perspective for a Cyber-Physical Multi-Agent Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Combined Model and Data-Driven Approach for the Determination of Rotor Temperature in an Induction Machine

by
Razvan Mocanu
*,
Alexandru Onea
and
Constantin Catalin Dosoftei
Department of Automatic Control and Applied Informatics, Gheorghe Asachi Technical University of Iasi, Bulevardul Profesor Dimitrie Mangeron 67, 700050 Iași, Romania
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(13), 4512; https://doi.org/10.3390/s21134512
Submission received: 10 May 2021 / Revised: 21 June 2021 / Accepted: 25 June 2021 / Published: 30 June 2021
(This article belongs to the Special Issue Cyber-Physical Systems - from Perception to Action)

Abstract

:
The need for protection of electrical machines comes as a demand of safety regulations in the automotive industry as well as a result of the general desire to obtain a robust and reliable electric powertrain. This paper introduces a hybrid method for estimating the temperature of the rotor of an Induction Machine (IM) based on a Nonlinear Autoregressive Network with Exogenous inputs (NARX) used as a prediction function within a particle filter. The temperature of the stator case is measured, and the information is used as an input to a NARX network and as a variable to a thermal process with first-order dynamics which serves as an observation function. Uncertainties of the NARX and thermal model are determined and used to correct the posterior estimate. Experimental data are used from a real IM test-bench and the results prove the applicability and good performance.

1. Introduction

The topic discussed in this paper is contextualized in the field of electric machine protection. The goal is to protect the rotor of the Induction Machine (IM) from overheating to increase the lifetime of the electrical machine. Improper thermal protection can lead to serious defects and even life-threatening situations. The method presented here requires a collection of data representing a sparse range of operating points of the electrical machine (i.e., phase currents, voltages, angular speed, etc.). The RMS-current and the temperature of the stator are measured with onboard sensors along with the angular speed while the rotor current is derived from the current controller. With the a priori data we parametrized the thermal model which serves as an observation function and trained a Nonlinear Autoregressive Network with Exogenous inputs (NARX) neural network to predict state transitions of the temperature. We tested the algorithm with input data acquired from a Belt-Driven Booster (BDB) test-bench for 48 V Mild-hybrid Electrical Vehicles (Mild-HEVs).
To justify the importance of this work we must first contextualize it in the automotive industry. It can be difficult to include in the end product expensive components such as infrared sensors, mainly due to two reasons. First is the extra cost, which is hardly acceptable in a tough competitive industry. Secondly, the reliability and accessibility of the sensors play a key role. For example, often the sensor must be positioned in hard-to-reach locations to provide useful data. This is the case of infrared sensors used in performance electric motors with fast thermal dynamics. Impurities, long-term exposure to mechanical vibrations can affect the precision of the measurements. We meet these rough conditions often in vehicles.
Acquiring information about the temperature of the IM is important because the thermal behavior of the IM influences the performance of the electrical drive train. Electrical machines should be monitored and diagnosed as per safety criteria [1,2,3].
The IM can be exposed to temperatures ranging from 50 C to more than 150 C in Mild-HEVs. Asymmetric power supply caused by inverter faults can lead to thermal overload of the IM. Thermal overload can also be caused by inadequate cooling, air-gap eccentricity, or overstress of the electrical machine (i.e., requesting high power from the IM for an extended time) [4]. Aluminum, copper, or their alloys are primary materials for the rotor bars of IM, while laminated steel is used for the inner portion of the stator. To ensure the protection of the electrical machine, we need to know the melting points and understand the deformation characteristics of the components. If an early failure is not detected, deformation of the rotor bars will worsen, leading to the raising of the bars to where the stator core and stator windings are damaged. Implicitly, the motor performance is altered, and in the worst-case scenario, torque at the shaft becomes erratic. The algorithm points to solve the problem of determining the temperature of the rotor of an IM used in Mild-HEVs without using a rotor temperature sensor. For validation purposes, we acquire the temperature of the rotor through an infrared sensor. Towards this, a hole was drilled into the case of the IM and the infrared sensor was placed with the direct sight of the rotor. On the final product, namely the BDB, this temperature sensor will not be present due to reliability issues mentioned above and to keep the production costs low. The target was to place the infrared sensor as close as possible to the assumed hot-spot considering that in an IM most of the high temperatures are located at the joints of the aluminum/copper bars-end with the end rings. Previously, we implemented a particle filter for estimation of the hot-spot temperature on the stator windings of a Permanent Magnet Synchronous Machine (PMSM) [5]. Different approaches to solve this problem have been investigated. The Nelder-mead method is employed towards nonlinear parameter estimation of a thermal model and the impedance response over a temperature range is determined. The impedance is measured with a specialized instrument (N4L precision impedance/LCR analyzer) and relates to the temperature of the IM [6]. Thermal Impulse Response modeling is employed to forecast the temperature in different areas of an electric machine to localize the hot spots. A heat narrow impulse is applied convolving with the loss of the machine and the temperature is calculated. The method is suitable in the design phase of the electrical machines but is also useful because it provides prior information about the most critical areas [7]. The paper describes an attempt to estimate the spatial distribution of temperature of a permanent magnet in an electrical machine based on the estimation of BEMF harmonics. The authors propose an empirical quadratic law to capture the relation between the maximum and the minimum temperatures across an interior permanent magnet synchronous machine and they emphasize the need for a ground-based approach for the mathematical model [8]. A multiple-scale method is applied to obtain the homogenized heat equation [9]. The case of an automotive switched reluctance machine is considered, and the parameters of a lumped thermal network are identified [10]. In another study, it is pointed out the effect of disturbance of the resolver signal on the temperature of an IM. The authors propose a method of estimation of stator temperature starting from the fixed frame where a DC signal is injected, and the stator resistance is estimated in a field-oriented control scheme [11]. Pulse injection with zero voltage average is used to estimate rotor parameters [12]. Contrary to invasive methods, optimal or sub-optimal filters can be used to estimate the temperature of electrical machines considering the dependence between the magnetic flux and rotor temperature. These approaches require prior knowledge of the electrical subsystem and often assume Gaussian uncertainties [13,14]. A method is described for estimating the temperature of the magnet of a PMSM using the d q frame model within a Kalman filter. The VSI distortion term is included in the estimator [15]. The equivalent circuit parameters of an IM are estimated by a particle swarm optimization algorithm in steady-state conditions. The authors emphasis the applicability of the algorithm for continuous working conditions but mention the lack of precision in transient operations [16]. Estimation of internal resistance of the rotor through model reference adaptive is a technique used widely; however, it requires the knowledge of the machine nominal parameters to provide good results [17]. Additional search coils are mounted inside a reluctance machine on the stator windings to provide a direct measurement of stator flux and further provide an estimation of temperature and shaft position [18]. In open-end wingdings PMSM applications where small torque ripples are acceptable, zero current sequences are used to determine the stator resistance and further the temperature [19]. The time constant of a first-order thermal model is computed based on a polynomial function of the speed and stator temperature of a PMSM [20]. A design procedure for guaranteed state estimation of unmeasurable states with practical application to induction machines is discussed [21]. Related to the approach involving particle filter we mention a few remarks. The particle filter is based on an effective sequential Monte Carlo method to solve the recursive Bayesian filtering problem [22,23]. The computational power allows now to use the particle filter in a wider range of applications. In the automotive world is used to estimate the traffic conditions [24]. Temperature estimation of a satellite is implemented considering a double lumped thermal model [25]. Particle filters are used for fault prediction and diagnosis. IGBTs faults are predicted with application in power electronics [26]. A simulation framework with practical results for temperature estimation of IGBT is detailed in the work of authors [27]. Throughout the mentioned publications in this field we could practically divide the described methods into three categories: (a) methods that rely on additional hardware (e.g., search coils, real-time impedance measurement devices, magnetic flux sensors, etc.). (b) methods that are invasive and can alter the control command of the electrical machine (e.g., DC injection, impulse injection) and methods that rely on a precise thermal model and precise knowledge of nominal parameters.
The method proposed in this paper has the advantage of being simple despite the complicated-look of the particle filter, it is straightforward and provides good results. It is decoupled and can be modularized, allowing an update or change to the prediction and observation model in a ‘plug and play’ manner. It can be seen as three individual parts: two information channels and one merging strategy. Except for the fairly common particle filter, the other two information channels can be adjusted or replaced without any impact on the particle filter implementation, thus the overall implementation is somewhat simple. Furthermore, we mention the most important features:
  • The algorithm consumes data that is already available on most standard architectures of power electronics used in electric machine control (i.e., phase currents sensors, stator temperature sensor, position, and speed encoder).
  • The method does not require additional components (e.g., real-time precision impedance measurement, search coils, etc.)
  • The method is not invasive and does not alter the control commands to the electrical machine.
  • In the common NARX network, the state is updated with the actual output of the network. In the proposed algorithm, we update the state of the NARX network with the posterior estimate. The posterior estimate is obtained after merging with the particle filter the two information channels: the NARX and the thermal model. This corrects the prediction of the neural network in the prediction stage. Thus, is obtained a conditioned output of the NARX model with the thermal model which in the end improves the precision of the estimated value. This approach is distinctive, and it was not investigated before to the problem of temperature estimation.
The paper continues with a brief description of the self-contained BDB with IM in Section 2, while Section 3 introduces the mathematical models considered in the particle filter. In Section 4, the equations of the recursive Bayesian filter and particle filter are provided with the application for rotor temperature estimation. In Section 5 are presented the experimental results from the motor-load test-bench. Finally, Section 6 concludes this work and points to some possible improvements.

2. System Description—Belt-Driven Booster

The BDB consists of a three-phase IM with an integrated inverter and control unit which is capable of offering a continuous power of 6KW and a torque of 60NM on a 48 V power net system. The system is especially efficient in diesel hybrid cars. The nitrogen oxide emissions are reduced during accelerations from low speed when some of the torque is produced by the BDB. This system allows a fast engine start and energy recuperation. Additionally, engine-off coasting can be permitted (i.e., the internal combustion engine is shut down after the vehicle reaches a higher speed and the propulsion is ensured by the BDB). Additionally, the IM can be maintenance-free mainly, due to the lack of slip rings, thus reducing the costs. The commands to the BDB are sent through a Controller Area Network (CAN)/Flexray channel by the electrical vehicle main controller, thus the system presents itself as a self-contained, compact solution Overall, the 48 V systems (DC/DC, inverter) are, at the moment, one of the most affordable and quick solutions for electrification in the automotive industry. Figure 1 shows a conceptual representation of a 48 V traction system.
The dynamics of the electrical quantities of the IM are modeled by a 4th order system. The variables and parameters are expressed in the two-axis reference frame. The voltage equilibrium equations of the stator ( V d s , V q s ) and of the rotor ( V d r , V q r ) are defined by (1):
V q s = R s i q s + d ψ q s d t + ω ψ d s V d s = R s i d s + d ψ d s d t ω ψ q s V q r = R r i q r + d ψ q r d t + ( ω ω r ) ψ d r V d r = R r i d r + d ψ d r d t ( ω ω r ) ψ q r
where ω is the mechanical angular velocity, ω r is the electrical angular velocity, i d s , i q s are the stator currents and i d r , i q r are the rotor currents. The rotor is squirrel-cage type and V q r = 0 , V d r = 0 . R s and R r are the stator and rotor resistances, respectively. The stator fluxes ψ q s , ψ d s and the rotor fluxes ψ q r , ψ d r are defined by:
ψ q s = L s i q s + L m i q r ψ d s = L s i d s + L m i d r ψ q r = L r i q r + L m i q s ψ d r = L r i d r + L m i d s
where L s = L l s + L m , L r = L l r + L m define the total inductance of the stator and of the rotor, respectively. L m is the mutual inductance and L l s and L l r are the leakage inductance of the stator and rotor, respectively. i r is an estimated value of the magnetization current. In our case i r was approximated by low-pass filtering i d s with the time constant of the rotor electrical dynamics τ r = L r / R r . The speed-torque characteristics is plotted concerning the DC voltage bus in Figure 2. The parameters of the IM are listed in Table 1.
f s is the inverter switching frequency, η denotes the angular speed, V is the nominal DC voltage bus, P n is the nominal mechanical power. Index n is used for nominal values.

3. The Mathematical Models of the State Transition and Observation Functions

We propose to use the model defined by (3) considering that the major heat source in an IM is the rotor and the stator will follow to a certain extent the rotor temperature curve. The model is a simple heat transfer function of first-order which includes a thermal resistance through which heat is dissipated and a thermal capacitance that absorbs heat. τ is a time constant introduced for realizability. The model is not meant to be a precise description of the thermal process but rather to act as a filter/preconditioner of the NARX network.
τ d T α ( t ) d t + T α ( t ) = α 1 d T s ( t ) d t + α 2 T s ( t )
T α is the rotor temperature computed from the model (3) and T s is the measured stator temperature. The coefficients α R 2 x 1 are calculated offline with the least-squares equation:
α = ( C α * T C α * ) 1 C α * T T r *
where C α * R N x 2 is defined as:
C α * = [ Δ T s * T s * ]
and Δ T s * is the gradient of the stator temperature sensor:
Δ T s * = ( T s * [ 1 ] T s * [ 0 ] ) / h ( T s * [ N ] T s * [ N 1 ] ) / h
τ is the time constant of the low-pass filter introduced for realizability and h is the sampling time of the algorithm.
The NARX is a recurrent network with feedback connections encircling several layers of the network. The NARX structure is built on the linear ARX model, which is generally used in time-series modeling [28]. We use the structure of the recursive neural network with delay action because it accounts for the intrinsic memory of thermal processes.
The defining equation for the NARX model is:
T n e t [ k ] = n a r x ( i r [ k ] , i s [ k ] , η [ k ] , T s [ k ] , T n e t [ k 1 ] )
where i s the stator current and η is the angular speed. We added one hidden layer, fully connected, to the network structure which is enough to capture the dynamics of the predicted temperature T n e t . Increasing the size of the hidden layer or the number of layers did not contribute to significant improvements in the precision but rather to the increase of the computational run-time. Thus, f 1 is the hidden layer’s sigmoid activation function, while f 2 is the output layer’s linear activation function. The hidden layer’s weights and biases are W l i , j and b L i , whereas the output layer’s weights and biases are W L i and b 0 . The indexes are in the range of i = 1…10 and j = 1…4.

Data Acquisition and Network Training

Figure 3 illustrates the data acquisition chain. The data available at the level of the Electronic Control Unit (ECU) is the RMS stator current ( i s ) acquired through the phase current sensors, the stator temperature T s acquired through the stator thermocouple, the mechanical position of the shaft acquired through the position sensor. The ECU estimates the rotor temperature and calculates the angular speed from the position sensor. The role of the power analyzer is to provide analysis of the active and reactive power and can provide a redundant measurement of the phase currents, voltages, and torque. In the training process, we used data received from the ECU through the CAN bus transducer ( i r , i s , T s , η ) and the real rotor temperature from the infrared-sensor acquisition board. The currents are acquired with a sampling time of 100 μs, the angular speed is acquired with 10 ms, and the temperature with a sampling time of 100 ms. Because of the different sampling times, we performed linear interpolation for the low-sampling variables to have the same size of data before the training process. Figure 4 illustrates the structure of the neural network.
For training the network we recorded different operating points with ranging speeds from 0 to 15,000 RPM and torque from 0 to 50 Nm summing a total of approximately 5 h of recordings.
The model was validated on the test set by comparing the estimated temperature T ^ n e t with the real temperature acquired with the infrared rotor temperature sensor T r  Figure 5 shows the error distribution of the NARX network.
The training was carried out using M a t l a b ’s neural network utility (trainlm) and included two steps of preparation: first, converting data from D A T representation (test-bench data acquired via XCP protocol) to M A T format ( M a t l a b ’s particular format), and second, normalizing all variables prior to the training process. The weights W l i , j , W L i and biases b L i , b 0 are updated during the backpropagation training with the Levenberg-Marquardt (LM) optimization algorithm. For training, we used the open-loop architecture in which the real temperature of the rotor is used as feedback in the backpropagation algorithm instead of the estimated temperature. This offers two benefits. The first is that the feedforward network’s input is more accurate. The second benefit is that the resulting network has a feedforward design, allowing for the usage of a more efficient training procedure such as the LM method. Comprehensive literature and one of the earliest applications of the LM method to neural network training is documented by the authors of [29,30]. This approach is still one of the fastest for training intermediate-size feedforward neural networks of up to several hundred weights [30]. The implementation of the LM algorithm is fairly common, and a Matlab-specific example may be consulted here [31].
The error Probability Density Function (PDF) of the neural network and observation model as derived from the test-bench data are depicted in Figure 5 and Figure 6.

4. Particle Filter for Rotor Temperature Estimation

Figure 7 is a block representation of the estimation strategy. First, the inputs stator temperature T s , rotor current i r , stator current i s , angular speed η and the previous posterior temperature set T p o s t feed the recursive neural network block and several N p particles are obtained. The PF block pulls a measurement out of the observation model and merges the uncertainties PDF (i.e., Gaussian distributed R m and R n e t ) around the newly estimated temperature.
The motivation of using the particle filter for the task of merging data from models (3) and (7) is that the predicted temperature with the NARX network is nonlinear and both models can reveal non-Gaussian noise. In our thermal process the uncertainties can be roughly approximated with the Gaussian distribution (Figure 5 and Figure 6). However, this can vary with the positioning and precision of the sensors. Except for the few mentioned papers, the particle filter was not used in the temperature estimation but it can undertake this task well since the heat exchange is a slow process and computational power is enough for this task. In the prediction step of the recursive Bayesian filter, the temperature T n e t is calculated as in (7). The corresponding PDF of the prior is given by the Chapman-Kolmogorov equation [22,23] adapted to our set of data d:
p T n e t k / d 1 : k 1 prior PDF = p ( T n e t k / T n e t k 1 , T s k 1 , η k 1 , . . . . i r k 1 , i s k 1 ) transition model ( NARX net function ) · p T n e t , k 1 / d 1 : k 1 posterior PDF d T n e t k 1
where p denotes the probability, d defines the set of data (will be define below) and k is the iteration index. Thus, to calculate the prior rotor temperature the NARX net is sampled (i.e., Equation (7) is evaluated) and the weighting is done according to the prediction of the observation model (3). Consequently, the prior precision depends on the precision of the NARX prediction and posterior precision. Following, the posterior PDF is derived from the prior density function of the neural network (8) and from the observation function (3) [32,33] with a normalization factor:
p T n e t k / d 1 : k posterior PDF = p T α k / T n e t k observation model p T n e t k / d 1 : k 1 prior PDF / p T α k / d 1 : k 1 normalization PDF
The normalization PDF is defined as [22]:
p T α k / d 1 : k 1 = p T α k / T n e t k p T n e t k / d 1 : k 1
In the previous equations d 1 : k and d 1 : k 1 define the set of data (i.e., observation model predicted temperature T α , stator temperature T s , angular speed η , rotor current i r and stator current i s ) available at recurrence k and k 1 , respectively:
d 1 : k = T α 1 : k , T s 1 : k , η 1 : k , i r 1 : k , i s 1 : k d 1 : k 1 = T α 1 : k 1 , T s 1 : k 1 , η 1 : k 1 , i r 1 : k 1 , i s 1 : k 1
The PF approximates the PDFs with a finite number of samples, rather than computing (10) at each sample time which can be time consuming in the case of complex analytical forms. As a result, a series of samples (particles) characterizes the PDF of the posterior temperature p ( T n e t k / d 1 : k ) . An approximation of the functional form of the posterior PDF is obtained for an infinitely large number of particles. The posterior PDF in its discrete form is defined by [32]:
p T n e t k / d 1 : k i = 1 N p w k ( i ) δ ( T n e t k T n e t k ( i ) )
The PF generates a set of rotor temperatures and their corresponding weights { ( w k ( i ) , T n e t k ( i ) ) : i = 1 N p } . N p is the number of temperature points in the set. In this work N p = 60 was determined experimentally and provides satisfactory results. The weights are chosen at each sample time according to the principle of ‘importance sampling’. Hence, the prior PDF is represented as [32]:
p T n e t k / d 1 : k 1 prior PDF i = 1 N p w k 1 ( i ) p T n e t k / T n e t k 1 ( i ) , T s k , η k , i r k , i s k transition model ( NARX net function )
The PDF of the posterior temperature is given by [22,34]:
p T n e t k / d 1 : k posterior PDF ψ p T α k / T n e t k observation model p T n e t k / d 1 : k 1 prior PDF
In (14) ψ is the normalization constant analogous to (9).
A three-step process is followed to implement the PF: predict, update, re-sample. The NARX function (7) with inputs T s , i r , i s and η is evaluated in the prediction step for each temperature point from the set of posterior estimates. There is a total of N p temperature points in the posterior set. The obtained vales at the output of the network represent the projected rotor temperatures which are distributed conforming to the prior PDF (13) [22,34].
T n e t k ( i ) p T n e t k / d 1 : k 1
Following, the observation model (3) is evaluated in the update step and a set of N p rotor temperature points are weighted based on the PDF of the observation function. For each projected temperature, the weight is determined as follows [22]:
w k ( i ) = p T α k / T n e t k
The weights w k ( i ) approximate the densities of the rotor temperature points and are normalized so that:
i = 1 N p w k ( i ) = 1
Consequently, the  N p particles will be dispersed conforming to:
T n e t k ( i ) p T n e t k / d 1 : k
The effect of degeneracy is avoided in the re-sample step where a new group of N p particles are generated from the posterior PDF (14). In this step we integrate (i.e., sum) the values in the set w for each particle and we obtain a cumulative real sum denoted in Algorithm 1 as c s w [ 0 , 1 ] . In the continuous-time domain, this is analog to the cumulative distribution function. We generate a set X of N p real values distributed uniformly X [ 0 , 1 ] . We used a built-in function that will not be covered in this paper. However, a good reference for uniformly distributed numbers generator can be found here [35]. In the next step, we searched the indices of X for each particle that correspond to the c s w . These represent the re-sampled indices of the prior estimate T n e t . Following, the posterior estimation of the rotor temperature is determined as the mean value of new posterior particles:
T ^ r k = 1 N p i = 1 N p T n e t k ( i )
The re-sampling part of the algorithm is the most computational expensive from the particle filter since it requires the generation of random variables from given distributions. It is more expensive than inferring the neural network. The computational complexity of the re-sampling stage is linear increase with the number of particles. This complexity could be decreased if, for example, we would use a hardware generator of uniform distributions.
The implementation of the estimation algorithm can be synthesized as shown in Algorithm 1. The algorithm is illustrated for 1 sample time.
The parameters of the observation model and particle filter are listed in Table 2.
Figure 8 details the test-bench setup.
Algorithm 1 Particle filter for rotor temperature estimation
Initialization
Initialize posterior distribution: p T n e t , k 1 / d 1 : k 1 p 0
Generate N p particles with initial distribution p 0 and place them in a set called T p o s t .
Run
Compute  T α [ k ]
T α [ k ] = h ( h + τ ) h τ T α [ k 1 ] α 1 h T s [ k 1 ] + ( α 1 + h α 2 ) h T s [ k ]
while i N p do
   Predict
   T n e t [ i ] = n a r x ( i r [ k ] , i s [ k ] , η [ k ] , T s [ k ] , T p o s t [ i ] )
   Update
  Calculate the weight associated with each projected state:
   w ( i ) = 1 / 2 π R m e x p [ ( T α [ k ] T n e t ( i ) ) 2 / ( 2 R m ) ] ;
  i = i + 1
end while
Normalize weights : w ( i ) = w ( i ) i = 1 N p w ( i )
Re - sample
Sample the discrete distribution w to generate N p random samples which will represent the posterior T p o s t :
X = r a n d ( N p , 1 ) ; fetch N p values from a uniform distribution X [ 0 , 1 ] ;
for j = 1: N p  do
   c s w ( j ) = L = 1 j w ( L ) ; Cumulative Sum of W; c s w [ 0 , 1 ]
end for
for j = 1: N p  do
   i p ( j ) = a r g m i n i | X ( i ) c s w ( j ) | ; i p will be the index of the posterior after re-sampling
end for
for j = 1: N p  do
   T p o s t ( j ) = T n e t [ i p ( j ) ]
end for
T ^ r ( k ) = 1 N p i = 1 N p T p o s t ( i )

5. Results and Discussions

To simulate and validate the temperature estimation method, recorded data of rotor current, stator current, and stator temperature along with angular speed information is provided to the estimation algorithm. The IM is set in torque control mode while the load machine is controlled in speed mode. The angular velocity of the load machine is controlled by a dedicated voltage source inverter with own control unit. Torque steps with various amplitudes are requested from the BDB. Different ranges of torque and speed setpoints are measured. Validation is done on separate data sets.
In the first scenario represented in Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 the IM is accelerated up to 9000 RPM and decelerated to 0 RPM, 0 A to capture part of the cooling behavior. The phase current peaks are up to 350 A and are caused by peak torque requests as can be seen in Figure 12. The magnetization current has a peak of 250 A. The raw (unfiltered) output of the NARX is seen in Figure 9. With the first-order thermal model process the effect is of filtering high frequencies at the output of the NARX network. This is accomplished by the particle filter in the update step where the prior estimated is smoothed considering the uncertainty model of the recursive network. As expected according to the Chapman-Kolmogorov the effect is of filtering ripples in the prior estimate (according to the transition PDF). In all experiments, the initial temperature guess is the stator temperature.
Higher speed interval is tested in the case represented in Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18. The speed is varied in the range of 1000 RPM to 15,000 RPM with different steps. We allow partial cooling of the machine, then we repeat the speed and torque requests. The RMS stator current has a peak of 360 A, while the magnetization current reaches approximately 250 A. Cooling of the machine is captured well by both the thermal model and NARX network and overall we obtain a precision with less than 4 error.
In the test case represented in Figure 19, Figure 20, Figure 21, Figure 22 and Figure 23 the IM is rotated nearly to the maximum speed. Over the entire test set the estimated temperature shows an average deviation of 3.3 C (Figure 24).
The performance metrics for the considered test cases are summarized in Table 3.
Figure 24 depicts the error histogram for the entire test set. The parameters of the neural network are provided in Table 4.

6. Conclusions

I. There are different development perspectives of this estimation strategy. A recurrent neural network can offer a solution when the thermal modeling capabilities are missing but a high quantity of data are available. Training a neural network involves time for acquiring collections of data. At no supplementary development time, the same data are used to identify the parameters of the observation model and to obtain the uncertainty of the model. Then, we merge both channels through the particle filter. The observation or the transition models can be improved but this means investing more time in the development of a mathematical description with a deeper insight into the thermal physics. The transition and observation functions can be a well-defined mathematical characterization of the thermal process of multiple variables instead of an empirically structured neural network. Thus, one who must decide between these two paths must evaluate the capabilities of the development team and the time constraints. Obtaining a detailed thermal model could involve more time and effort than a simple job of acquiring data on a test-bench, train a network and identify coefficients of a simple first-order process. We saw that a simple first-order process captures well to some extent the cause-effect of the relation rotor temperature → stator temperature. Not enough, but this is compensated by the particle filtering with the neural network’s prediction.
II. It is important to highlight that to obtain a precise online estimation of the rotor temperature, the appropriate experiment should be performed built around the targeted application. Therefore, it is advised that the data measurement for identification be done on the actual environment of the IM (e.g., with the BDB mounted inside the vehicle). This is due to possibly different temperature behavior and other heating sources. The external temperature sensor can be placed not only on the stator windings but also on the stator case or at proximity to the electrical machine (such that it exists somewhat thermal conductivity). The sensor could also play the role of an environment temperature sensor and we expect that the NARX network can accommodate during training and can capture the contribution of the external sensor to the estimated temperature of the rotor. We believe that this method is general enough to apply to a wide range of electrical machines and other types of thermal processes. Further development of this work will be in the direction of (a) increasing the precision of the thermal model for this particular application (b) develop a theoretical generic framework for the combined model and data-driven estimation approach to be used in other applications.

Author Contributions

Conceptualization, R.M. and A.O.; methodology, R.M., A.O. and C.C.D.; software, R.M.; validation, R.M., A.O. and C.C.D.; formal analysis, R.M., A.O. and C.C.D.; investigation, R.M., A.O. and C.C.D.; resources, R.M., A.O. and C.C.D.; data curation, R.M.; writing—original draft preparation, R.M.; writing—review and editing, R.M., C.C.D.; visualization, R.M.; project administration, C.C.D.; funding acquisition, C.C.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant of the Romanian National Authority for Scientific Research and Innovation, CNCS/CCCDI – UEFISCDI, project no. PN-III-P2-2.1-PTE-2019-0731, contract no. 19/2020, within PNCDI III.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BDBBelt-Driven Booster
CANController Area Network
ECUElectronic Control Unit
HEVHybrid Electrical Vehicle
IMInduction Machine
LMLevenberg-Marquardt
Mild-HEVsMild-hybrid Electrical Vehicles
NARXNonlinear AutoRegressive network with eXogenous inputs
PDFProbability Density Function
PFParticle Filter
PMSMPermanent Magnet Synchronous Machine

List of Symbols

The following nomenclature is used throughout the manuscript:
V d s , V q s [V] voltage equilibrium of the stator in rotating frame;
V d r , V q r [V] voltage equilibrium of the rotor in rotating frame;
ω [rad/s] mechanical angular velocity;
ω r [rad/s] electrical angular velocity;
i d s , i q s stator currents in rotating frame;
i d r , i q r [A] rotor currents in rotating frame;
ψ q s , ψ d s [Wb] stator fluxes in rotating frame;
ψ q r , ψ d r [Wb] rotor fluxes in rotating frame;
L s [H] inductance of the stator;
L r [H] inductance of the rotor;
L m [H] mutual inductance;
L l s [H] leakage inductance of the stator;
L l r [H] leakage inductance of the rotor;
i r [A] estimated value of the magnetization current;
i s [A] stator current;
τ r [s] time constant of the rotor electrical dynamics;
f s [Hz] inverter switching frequency;
R r [ Ω ] nominal rotor resistance;
R s [ Ω ] nominal stator resistance;
R m , R n e t normal distribution variances of observation and transition model;
τ [s] time constant of the low-pass filter introduced for realizability;
α 1 , α 2 coefficients of the thermal model;
T s [ C] stator temperature;
T r [ C] rotor temperature;
T α [ C] predicted temperature with the thermal model;
T n e t [ C] predicted temperature with the NARX model;
T ^ r [ C] posterior-estimated value of the rotor temperature;
h[s] sample time;
η [RPM] mechanical speed of the rotor measured in revolutions per minute;
Δ gradient;
reference value;
w[-] particle filter sample weights;
f 1 sigmoid activation function of the hidden layer;
f 2 linear activation function of the output layer;

References

  1. Sharma, V.K.; Usman, A.; Rajpurohit, B.S. Fault Diagnosis of BLDC Motor Drive using Vibration Spectrum Analysis: Part II. In Proceedings of the 2020 IEEE International Power and Renewable Energy Conference, Karunagappally, India, 30 Octomber–1 November 2020; pp. 1–6. [Google Scholar] [CrossRef]
  2. Usman, A.; Rajpurohit, B.S. Time-Efficient Fault Diagnosis of a BLDC Motor Drive Deployed in Electric Vehicle Applications. In Proceedings of the 2020 IEEE Global Humanitarian Technology Conference (GHTC), Seattle, WA, USA, 29 Octomber–1 November 2020; pp. 1–5. [Google Scholar] [CrossRef]
  3. Park, H.; Suh, Y. Fault-Tolerant Control Strategy for Reduced Torque Ripple of Independent Twelve-phase BLDC Motor Drive System under Open-Circuit Faults. In Proceedings of the 2020 IEEE Energy Conversion Congress and Exposition (ECCE), Detroit, MI, USA, 11–15 October 2020; pp. 3370–3375. [Google Scholar] [CrossRef]
  4. Emadi, A.; Lee, Y.J.; Rajashekara, K. Power electronics and motor drives in electric, hybrid electric, and plug-in hybrid electric vehicles. IEEE Trans. Ind. Electron. 2008, 55, 2237–2245. [Google Scholar] [CrossRef]
  5. Mocanu, R.; Onea, A. Determination of stator temperature for thermal protection in a Permanent Magnet Synchronous Machine. In Proceedings of the 2017 25th Mediterranean Conference on Control and Automation (MED), Valletta, Malta, 3–6 July 2017; pp. 1321–1325. [Google Scholar]
  6. Agbaje, O.; Kavanagh, D.; Sumisławska, M.; Howey, D.; McCulloch, M.; Burnham, K. TEstimation of Temperature Dependent Equivalent Circuit Parameters for Traction-Based Electric Machines. In Proceedings of the IET Hybrid and Electric Vehicles Conference 2013 (HEVC 2013), London, UK, 6–7 November 2013. [Google Scholar]
  7. Rahman, N.A.; Gu, L.; Bostanci, E.; Fahimi, B. Temperature estimation of switched reluctance machines using thermal impulse response technique. In Proceedings of the 2016 IEEE Conference on Electromagnetic Field Computation (CEFC), Miami, FL, USA, 13–16 November 2016; p. 1. [Google Scholar]
  8. Reigosa, D.D.; Fernandez, D.; Tanimoto, T.; Kato, T.; Briz, F. Permanent-magnet temperature distribution estimation in permanent-magnet synchronous machines using back electromotive force harmonics. IEEE Trans. Ind. Appl. 2016, 52, 3093–3103. [Google Scholar] [CrossRef] [Green Version]
  9. Romanazzi, P.; Bruna, M.; Howey, D.A. Thermal homogenization of electrical machine windings applying the multiple-scales method. J. Heat Transf. 2017, 139, 012101. [Google Scholar] [CrossRef] [Green Version]
  10. Qi, F.; Ralev, I.; Klein-Hessling, A.; De Doncker, R.W. Online temperature estimation of an automotive switched reluctance motor using space-resolved lumped parameter network. In Proceedings of the 2016 19th International Conference on Electrical Machines and Systems (ICEMS), Chiba, Japan, 13–16 November 2016; pp. 1–6. [Google Scholar]
  11. Baneira, F.; Asiminoaei, L.; Doval-Gandoy, J.; Delpino, H.A.M.; Yepes, A.G.; Godbersen, J. Estimation method of stator winding resistance for induction motor drives based on dc-signal injection suitable for low inertia. IEEE Trans. Power Electron. 2018, 34, 5646–5654. [Google Scholar] [CrossRef]
  12. Fahrner, W.; Vogelsberger, M.A.; Wolbank, T. A New Technique to Identify Induction Machine Rotor Parameters During Dynamic Operation and Low Speed. In Proceedings of the 2018 IEEE 18th International Power Electronics and Motion Control Conference (PEMC), Budapest, Hungary, 26–30 August 2018; pp. 471–476. [Google Scholar]
  13. Min, Y.; Huang, W.; Yang, J.; Zhao, Y. On-line Estimation of Permanent-magnet Flux and Temperature Rise in Stator Winding for PMSM. In Proceedings of the 2019 22nd International Conference on Electrical Machines and Systems (ICEMS), Harbin, China, 11–14 August 2019; pp. 1–5. [Google Scholar]
  14. Tran, T.V.; Nègre, E. Efficient Estimator of Rotor Temperature Designing for Electric and Hybrid Powertrain Platform. Electronics 2020, 9, 1096. [Google Scholar] [CrossRef]
  15. Feng, G.; Lai, C.; Tjong, J.; Kar, N.C. Noninvasive Kalman filter based permanent magnet temperature estimation for permanent magnet synchronous machines. IEEE Trans. Power Electron. 2018, 33, 10673–10682. [Google Scholar] [CrossRef]
  16. Bijan, M.G.; Pillay, P. Efficiency estimation of the induction machine by particle swarm optimization using rapid test data with range constraints. IEEE Trans. Ind. Electron. 2018, 66, 5883–5894. [Google Scholar] [CrossRef]
  17. Zhao, H.; Eldeeb, H.H.; Wang, J.; Zhan, Y.; Xu, G.; Mohammed, O.A. Online Estimation of Rotor Temperature in Induction Motors Based on Parameter Identification. In Proceedings of the 2019 IEEE Energy Conversion Congress and Exposition (ECCE), Baltimore, MD, USA, 29 September–3 October 2019; pp. 1629–1634. [Google Scholar]
  18. Tornello, L.D.; Scelba, G.; Scarcella, G.; Cacciato, M.; Testa, A.; Foti, S.; de Caro, S.; Pulvirenti, M. Combined Rotor-Position Estimation and Temperature Monitoring in Sensorless, Synchronous Reluctance Motor Drives. IEEE Trans. Ind. Appl. 2019, 55, 3851–3862. [Google Scholar] [CrossRef]
  19. Hunter, N.; Cox, T.; Zanchetta, P.; Odhano, S.A.; Rovere, L. Non-intrusive online stator temperature estimation for open-end winding PMSM. In Proceedings of the 2018 IEEE Energy Conversion Congress and Exposition (ECCE), Portland, OR, USA, 23–27 September 2018; pp. 6845–6851. [Google Scholar]
  20. Carpiuc, S.C. Rotor temperature detection in permanent magnet synchronous machine-based automotive electric traction drives. IEEE Trans. Power Electron. 2016, 32, 2090–2097. [Google Scholar] [CrossRef]
  21. Schwartz, M.; Krebs, S.; Hohmann, S. Guaranteed State Estimation Using a Bundle of Interval Observers with Adaptive Gains Applied to the Induction Machine. Sensors 2021, 21, 2584. [Google Scholar] [CrossRef] [PubMed]
  22. Thrun, S. Particle filters in robotics. In Proceedings of the Eighteenth Conference on Uncertainty in Artificial Intelligence, Edmonton, AB, Canada, 1–4 August 2002; Morgan Kaufmann Publishers Inc.: San Francisco, CA, USA, 2002; pp. 511–518. [Google Scholar]
  23. Marconato, A.; Schoukens, M.; Schoukens, J. Filter-based regularisation for impulse response modelling. IET Control Theory Appl. 2016, 11, 194–204. [Google Scholar] [CrossRef] [Green Version]
  24. Sutarto, H.Y.; Boel, R.K.; Joelianto, E. Parameter estimation for stochastic hybrid model applied to urban traffic flow estimation. IET Control Theory Appl. 2015, 9, 1683–1691. [Google Scholar] [CrossRef] [Green Version]
  25. Li-ping, P.; Hongquan, Q. Particle Filtering approach to parameter estimate and temperature prediction of satellite. In Proceedings of the 2008 7th World Congress on Intelligent Control and Automation, Chongqing, China, 25–27 June 2008; pp. 3001–3004. [Google Scholar]
  26. Haque, M.S.; Choi, S.; Baek, J. Auxiliary particle filtering-based estimation of remaining useful life of IGBT. IEEE Trans. Ind. Electron. 2017, 65, 2693–2703. [Google Scholar] [CrossRef]
  27. Lim, H.; Hwang, J.; Kwon, S.; Baek, H.; Uhm, J.; Lee, G. A Study on Real Time IGBT Junction Temperature Estimation Using the NTC and Calculation of Power Losses in the Automotive Inverter System. Sensors 2021, 21, 2454. [Google Scholar] [CrossRef] [PubMed]
  28. Diaconescu, E. The use of NARX neural networks to predict chaotic time series. Wseas Trans. Comput. Res. 2008, 3, 182–191. [Google Scholar]
  29. Hagan, M.T.; Menhaj, M.B. Training feedforward networks with the Marquardt algorithm. IEEE Trans. Neural Netw. 1994, 5, 989–993. [Google Scholar] [CrossRef] [PubMed]
  30. Hagan, M.; Demuth, H.; Beale, M.; De Jesus, O. Neural Network Design; PWS Pub. Co.: Boston, MA, USA, 1996. [Google Scholar]
  31. Mathworks. Implementation of Levenberg-Marquardt backpropagation—MATLAB. Levenberg-Marquardt Backpropagation—MATLAB. Available online: https://www.mathworks.com/help/deeplearning/ref/trainlm.html (accessed on 25 June 2021).
  32. Fox, V.; Hightower, J.; Liao, L.; Schulz, D.; Borriello, G. Bayesian filtering for location estimation. IEEE Pervasive Comput. 2003, 2, 24–33. [Google Scholar] [CrossRef] [Green Version]
  33. Ryan, A.; Hedrick, J.K. Particle filter based information-theoretic active sensing. Robot. Auton. Syst. 2010, 58, 574–584. [Google Scholar] [CrossRef]
  34. Thallas, A.; Tsardoulias, E.; Petrou, L. Particle filter—Scan matching hybrid SLAM employing topological information. In Proceedings of the 2016 24th Mediterranean Conference on Control and Automation (MED), Athens, Greece, 21–24 June 2016; pp. 226–231. [Google Scholar]
  35. L’Ecuyer, P. Uniform random number generation. Ann. Oper. Res. 1994, 53, 77–120. [Google Scholar] [CrossRef]
Figure 1. A conceptual illustration of a 48 V electrification system.
Figure 1. A conceptual illustration of a 48 V electrification system.
Sensors 21 04512 g001
Figure 2. The speed-toque characteristic w.r.t to the DC voltage bus.
Figure 2. The speed-toque characteristic w.r.t to the DC voltage bus.
Sensors 21 04512 g002
Figure 3. Illustration of the data acquisition chain.
Figure 3. Illustration of the data acquisition chain.
Sensors 21 04512 g003
Figure 4. NARX Neural Network.
Figure 4. NARX Neural Network.
Sensors 21 04512 g004
Figure 5. The error distribution of the neural network.
Figure 5. The error distribution of the neural network.
Sensors 21 04512 g005
Figure 6. The PDF of the observation model.
Figure 6. The PDF of the observation model.
Sensors 21 04512 g006
Figure 7. Block representation of the estimator.
Figure 7. Block representation of the estimator.
Sensors 21 04512 g007
Figure 8. View of the test-bench.
Figure 8. View of the test-bench.
Sensors 21 04512 g008
Figure 9. Estimated and actual rotor temperature (test scenario I).
Figure 9. Estimated and actual rotor temperature (test scenario I).
Sensors 21 04512 g009
Figure 10. Post estimate temperature deviation (test scenario I).
Figure 10. Post estimate temperature deviation (test scenario I).
Sensors 21 04512 g010
Figure 11. Shaft angular speed (test scenario I).
Figure 11. Shaft angular speed (test scenario I).
Sensors 21 04512 g011
Figure 12. Stator current (RMS) (test scenario I).
Figure 12. Stator current (RMS) (test scenario I).
Sensors 21 04512 g012
Figure 13. Rotor current (RMS) (test scenario I).
Figure 13. Rotor current (RMS) (test scenario I).
Sensors 21 04512 g013
Figure 14. The estimated and the actual rotor temperature (test case II).
Figure 14. The estimated and the actual rotor temperature (test case II).
Sensors 21 04512 g014
Figure 15. Post estimate temperature deviation (test scenario II).
Figure 15. Post estimate temperature deviation (test scenario II).
Sensors 21 04512 g015
Figure 16. Shaft angular speed (test scenario II).
Figure 16. Shaft angular speed (test scenario II).
Sensors 21 04512 g016
Figure 17. Stator current (RMS) (test scenario II).
Figure 17. Stator current (RMS) (test scenario II).
Sensors 21 04512 g017
Figure 18. Rotor current (RMS) (test scenario II).
Figure 18. Rotor current (RMS) (test scenario II).
Sensors 21 04512 g018
Figure 19. Estimated and actual rotor temperature (test scenario III).
Figure 19. Estimated and actual rotor temperature (test scenario III).
Sensors 21 04512 g019
Figure 20. Post estimate temperature deviation (test scenario III).
Figure 20. Post estimate temperature deviation (test scenario III).
Sensors 21 04512 g020
Figure 21. Shaft angular speed (test scenario III).
Figure 21. Shaft angular speed (test scenario III).
Sensors 21 04512 g021
Figure 22. Stator current (RMS) (test scenario III).
Figure 22. Stator current (RMS) (test scenario III).
Sensors 21 04512 g022
Figure 23. Rotor current (RMS) (test scenario III).
Figure 23. Rotor current (RMS) (test scenario III).
Sensors 21 04512 g023
Figure 24. The error histogram of the estimated temperatures over the entire test set.
Figure 24. The error histogram of the estimated temperatures over the entire test set.
Sensors 21 04512 g024
Table 1. The parameters of the induction machine.
Table 1. The parameters of the induction machine.
R r [ Ω ] @ 20 °C R s [ Ω ] @ 20 °C L ls [ H ] L lr [ H ] L m [ H ] Pole Pairs
0.00230.0024 5.23 × 10 6 5.23 × 10 6 1.586 × 10 4 4
f s [ K h z ] P n [ K W ] V n [ V ] T m a x [ N m ] η m a x [ R P M ] η n [ R P M ]
2006486018,0004500
Table 2. Table of parameters for the observation model and particle filter.
Table 2. Table of parameters for the observation model and particle filter.
Model CoefficientModel CoefficientFilter Time Constant
α 1 α 2 τ [ s ]
14.80521.33320.01
Sample TimeMeasurement NoiseNumber of Particles
h [ s ] R m [°C] N p
0.00154.364160
Table 3. Performance metrics for test scenarios 1–3.
Table 3. Performance metrics for test scenarios 1–3.
Test ScenarioMSEMAEMAXVAF
I3.421.30 °C6.90 °C98.15%
II2.861.16 °C6.31 °C98.52%
III3.801.79 °C7.21 °C97.56%
Table 4. Table of weights and biases of the NARX network for 100 epochs.
Table 4. Table of weights and biases of the NARX network for 100 epochs.
W l 1 , 1 W l 1 , 2 W l 1 , 3 W l 1 , 4 W l 1 , 5 W l 1 , 6 W l 1 , 7 W l 1 , 8 W l 1 , 9 W l 1 , 10
0.14031.21660.1194 1.4536 1.1580 1.4596 0.0907 0.97940.08611.1758
W l 2 , 1 W l 2 , 2 W l 2 , 3 W l 2 , 4 W l 2 , 5 W l 2 , 6 W l 2 , 7 W l 2 , 8 W l 2 , 9 W l 2 , 10
2.0016 0.5212 1.8543 1.4572 0.1994 1.4414 0.0046 1.7634 0.6388 0.4281
W l 3 , 1 W l 3 , 2 W l 3 , 3 W l 3 , 4 W l 3 , 5 W l 3 , 6 W l 3 , 7 W l 3 , 8 W l 3 , 9 W l 3 , 10
0.5346 1.2230 0.8087 0.1372 1.4326 1.0812 1.6871 0.6863 0.70850.3366
W l 4 , 1 W l 4 , 2 W l 4 , 3 W l 4 , 4 W l 4 , 5 W l 4 , 6 W l 4 , 7 W l 4 , 8 W l 4 , 9 W l 4 , 10
1.9551 0.7925 1.0801 0.6204 0.3679 0.2235 0.5951 0.85612.05621.2913
W D 1 W D 2 W D 3 W D 4 W D 5 W D 6 W D 7 W D 8 W D 9 W D 10
0.88330.02231.0860 2.0130 0.1146 0.9997 1.4808 0.3406 1.8001 1.6212
W L 1 W L 2 W L 3 W L 4 W L 5 W L 6 W L 7 W L 8 W L 9 W L 10
0.6887 0.4060 0.3728 0.4766 0.7304 0.6666 0.5783 0.5917 0.7903 0.7827
b L 1 b L 2 b L 3 b L 4 b L 5 b L 6 b L 7 b L 8 b L 9 b L 10
2.0631 2.0804 0.7474 0.1600 0.5823 0.0132 0.2237 0.4601 1.74742.1912
b o ---------
0.0661 ---------
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mocanu, R.; Onea, A.; Dosoftei, C.C. A Combined Model and Data-Driven Approach for the Determination of Rotor Temperature in an Induction Machine. Sensors 2021, 21, 4512. https://doi.org/10.3390/s21134512

AMA Style

Mocanu R, Onea A, Dosoftei CC. A Combined Model and Data-Driven Approach for the Determination of Rotor Temperature in an Induction Machine. Sensors. 2021; 21(13):4512. https://doi.org/10.3390/s21134512

Chicago/Turabian Style

Mocanu, Razvan, Alexandru Onea, and Constantin Catalin Dosoftei. 2021. "A Combined Model and Data-Driven Approach for the Determination of Rotor Temperature in an Induction Machine" Sensors 21, no. 13: 4512. https://doi.org/10.3390/s21134512

APA Style

Mocanu, R., Onea, A., & Dosoftei, C. C. (2021). A Combined Model and Data-Driven Approach for the Determination of Rotor Temperature in an Induction Machine. Sensors, 21(13), 4512. https://doi.org/10.3390/s21134512

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop