Next Article in Journal
Sensor Distribution Optimization for Composite Impact Monitoring Based on AR Model and LPP
Next Article in Special Issue
Replacing Induction Motors without Defined Efficiency Class by IE Class: Example of Energy, Economic, and Environmental Evaluation in 1.5 kW—IE3 Motors
Previous Article in Journal
Fault-Tolerant Predictive Current Control of Six-Phase PMSMs with a Single Isolated Neutral Configuration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stress-Dependent Magnetic Equivalent Circuit for Modeling Welding Effects in Electrical Steel Laminations

1
Department of Electromechanical, Systems and Metal Engineering, Ghent University, 9000 Ghent, Belgium
2
FlandersMake@UGent—Corelab EEDT-MP, 3001 Leuven, Belgium
3
Electrical Engineering Department, Kafrelshiekh University, Kafr El-Sheikh 33511, Egypt
*
Author to whom correspondence should be addressed.
Machines 2022, 10(12), 1153; https://doi.org/10.3390/machines10121153
Submission received: 7 November 2022 / Revised: 28 November 2022 / Accepted: 29 November 2022 / Published: 2 December 2022
(This article belongs to the Special Issue Modern Electrical Machines and Their Applications)

Abstract

:
Welding has a severe impact on the efficiency of electrical machines. The heat added during the welding process affects the microstructure of the material and causes residual stress. This results in local degradation of the magnetic permeability and facilitates additional iron losses in the machine core. With the purpose of modeling and simulating welding effects in electric machines, this paper proposes a stress-dependent magnetic equivalent circuit (MEC) model for welded non-grain-oriented electrical steel laminations. A modified iron loss model is proposed to accommodate these welding effects. Furthermore, the proposed MEC model is applied to a M270-35A stator core as a case study. It was demonstrated that the core losses increase by 25% when four welding joints are applied. With a limited number of magnetic measurements on a welded and unwelded core, the model can be fully parametrized. Finally, the model was successfully validated on a core with eight welding seams at 100 Hz. The proposed model can be integrated into the design of electric machines to consider the welding effects.

1. Introduction

Throughout the production process of electric machines, several techniques are applied to shape and stack the electrical steel laminations to make the stator and rotor core. After shaping and stacking the individual laminations, welding connects the laminations into a rigid stack and provides the torsion strength necessary for machine operation. An example of a welded stack of stator laminations is shown in Figure 1. Many researchers have studied the impact of welding on the core losses in electrical steel quantitatively, as reported in [1,2]. Others have studied the effects of welding on the microstructure of the material in and near the weld seam [3]. In [4], the residual stress near welding seams in electrical steel was measured and it was found that tensile stress appears within the weld seam while compressive stress appears in the heat-affected region near the weld seam. In [5], the relation between stress and the magnetic properties of electrical steel was reported. The variation of magnetic properties under stress can be linked to the impact on core losses of residual stress near welding seams. Additionally, in [6], it was noted that the welding process destroys the thin insulation layer between the lamination sheets, introducing short circuit paths, and resulting in increased classical losses. This effect is shown in the cross-section view in Figure 2.
It is clear from previous research that welding introduces residual mechanical stress in the area around the welding seam, which is known to change the local magnetic properties. Welding locally impairs the insulation layer between adjacent laminations, resulting in increased classical losses.
Recently, researchers have attempted to model the effects of welding in electrical steel in several ways. In [7], an analytical model is presented that shows the increase in classical losses when the weld bead radius increases. In [6], a 3D-FEM motor model was proposed where the yoke had uniform magnetic properties, which were identified using measurements on welded ring cores. In the [1] Epstein frame, measurements were used to parametrize a degradation model for magnetic material properties assuming homogenously damaged areas. It is considered that the model incorporates cutting effects adequately and allows incorporation of these effects in a finite-element method-based motor design software.
To the author’s knowledge, there are no reports of an analytical MEC model that accommodates local welding effects. The existing models apply either time-consuming finite element (FE) models or measurements with a standardized Epstein setup, which have different geometric properties than ring cores. Moreover, the existing models neglect the effects of residual stress or locally changing hysteretic material properties around the weld seam. A fast and accurate parametrized model for simulating welding effects on iron losses in electrical steel would be a valuable tool for studying the impact of each parameter. In this paper, a stress-dependent MEC model is presented that calculates the iron losses in welded stacks of electrical steel. The model is based on the known relations between mechanical stress and magnetic properties of non-grain oriented (NGO) electrical steel, with few geometrical and magnetic simplifications. Simultaneously, the model incorporates the effects of local eddy currents in the weld seam. Using Bertotti’s loss separation principle [8], a modified coarse-scale stress-dependent iron loss model is proposed. Finally, the model is applied to magnetic measurements on a welded stator core in the form of a case study.

2. Proposed Model

In this section, the specific aspects of the MEC-model will be presented. Further, the stress-dependency of the material model and loss model will be described in more detail. Afterward, in the Methodology section, it is explained how the model can be applied to calculate core losses in a welded stack of laminations. Table 1 contains the most important recurrent symbols in the model and Table 2 includes all the assumptions that were made.

2.1. MEC Model

In the model, the stator yoke of the machine is simplified as a ring core with a rectangular welded region and magnetically degraded region (Figure 2). Each welded section of the stator can be represented by a MEC, as shown in Figure 3. The weld seam geometry in this model was simplified using rectangular areas. The weld seam is located on the outer edge of the lamination. In the welding seam, the laminations are all electrically connected. In the region near the welding seam (region 2 in Figure 2), the heat from the welding process causes residual stress in the material, but the electrically insulating coating remains intact. For regions 2 and 3, a uniform stress state is assumed in the model. This stress state influences the magnetic permeability necessary for calculating the magnetic flux in each reluctance. In [9,10], measurement data shows the relation between the magnetic permeability and uni-axial (when magnetization and stress direction are parallel) stress state in NGO steels. It is assumed in this paper that these data are representative for each grade of NGO electrical steel. Based on these data stress-dependent correction factors for the relative permeability Kµ can be constructed as shown in Figure 4. However, these expressions are only valid when the direction of magnetization is parallel to the direction of stress, while it is known that the stress dependency of electrical steel is anisotropic. The influence of this stress anisotropy can be modeled using the equivalent stress approach. This approach was used by several researchers, such as in [11,12], and enables the calculation of magnetic properties under multi-axial stress when only uni-axial measurements are available. The equivalent stress model was adopted from [11] and is defined by:
σ e q = 1 K l n ( 2 exp ( K h T s h ) exp ( K t 1 T s t 1 ) + exp ( K t 2 T s t 2 ) )
where K is a material parameter (=4 × 10−9 for silicon-iron), s is the deviatoric part of the applied stress tensor (s = σ − (1/3)tr(σ)I), see [11] for more details, and I is the identity tensor. h, t1 and t2 are the direction vectors parallel to the applied field, orthogonal to the applied field and orthogonal to the sheet plane, respectively. In regions 2 and 3 (Figure 3), the stress is assumed to be directed radially away from the weld seam. This assumption corresponds to the results found in [2,4]. The stress in these areas is directed either 0°, 45° or 90° with respect to the magnetic flux direction (double-ended arrows in Figure 3). For the reluctances where the residual stress is not parallel to the magnetization direction, the corresponding σeq is calculated and used to determine the correction factors for each reluctance. In the proposed model, each reluctance in the network is defined by:
R ( σ e q , B ) = l r µ 0 K µ ( σ e q ) µ r ( B ) A r
where lr and Ar are the length and magnetic cross-section of each reluctance. The relative magnetic permeability µr is the stress-free permeability. In Figure 3, region 1 is the stress-free zone where Kµ = 1. Regions 2 and 3 are not stress-free and Kµ is determined by first calculating σeq and afterward applying Kµ based on Figure 4.

2.2. Coarse-Scale Stress Dependent Loss Coefficients

In the previous subsections, it is described how the MEC accounts for the residual stresses in and near the welding seam by assigning a stress-dependent magnetic permeability to each reluctance. Now, the MEC model can be solved for the flux density in each reluctance. However, the iron losses occurring in electrical steels also depend on stress [13]. In developing a stress-dependent iron loss model, a frequently used approach is to separate the core losses into hysteresis losses Physt, classical (eddy current) losses Pcl and excess losses Pexc and to investigate the relation to stress for each component:
P t o t = P h y s t + P c l + P e x c
When the flux density waveform is sinusoidal, the total core loss depends on the peak magnetic flux density Bp and frequency f following:
P t o t = C h B p 2 f + C c l B p 2 f 2 + C e x B p 1.5 f 1.5
where Ch, Ccl and Cex are grade-dependent coefficients. In [10,13,14,15], the stress-dependency of each loss coefficient was analyzed in more detail. These authors indicated the stress dependency of the hysteresis and excess loss coefficients clearly. Moreover, there appears to be a consensus concerning the independency of the classical loss to stress. In [14], the individual loss components are described as:
P h y s t = W h ( B p ) f
P c l = B p 2 f 2 d 2 π 2 ρ 6
P e x c = 8 ( 1 / ρ ) G A r V 0 ( B p ) ( B p f ) 1.5
where Wh is the hysteresis energy loss, ρ is the electrical resistivity, G = 0.1357, and d and S are the lamination thickness and cross-section, respectively. It can be assumed that the parameter V0 together with hysteresis energy loss Wh could describe the effect of the applied external mechanical stress on the microstructure and on the magnetic properties of NGO electrical steels, where V0 describes the influence on losses of various microstructural features such as (among others) grain size, crystallographic texture and residual stresses [14,16]. The models and data published by these authors enable the construction of a coarse-scale stress-dependent iron loss model. The stress-dependent hysteresis loss can be modeled using a correction factor Khyst as illustrated in Figure 5. Similarly, the stress dependency of the excess loss can be found in the V0 parameter, which can be modeled using a correction factor Kexc shown in Figure 6. However, B-dependency of V0 in Equation (7) is neglected in this model, which results in inaccuracies when the material goes into saturation. Therefore, in the case study, the model is validated on magnetic flux values up to 1.3 T to avoid modeling errors. Applying these correction factors to the iron loss model results in:
P h y s t ( σ e q ) = K h y s t ( σ e q ) C h B p 2 f
P e x c ( σ e q ) = K e x c ( σ e q ) C e x ( B p f ) 1.5

3. Local Classical Loss in Weld Seam

In the previous section, it was mentioned that the classical losses remain unaffected by the stress state. However, as the welding seam forms a connecting bar of electrical steel at the outer edge of the stack, the local eddy currents in this bar can be calculated using the classical loss equation but with thickness d equal to the weld seam radius Rweld (as illustrated in Figure 2). This way, the influence of the weld seam radius on the classical losses in the weld seam (region 3 in Figure 2) is included in the model:
P c l , w e l d = B p 2 f 2 R w e l d 2 π 2 ρ 6

3.1. Methodology

In this section, it is described in three steps how the model can be applied in any specific case of welded electrical steel. In Figure 7, a flowchart shows how the methodology can be applied step by step.

3.1.1. Model Construction

A parametrized MEC model should be constructed, as specified in Figure 3. In this figure, a MEC model is shown for a core with only one welding seam. Rweld, total length (l), total width (w) and thickness (d) should be set to the radius of the weld seam, the average yoke circumference, the average yoke width and lamination thickness, respectively. In the model, Rdeg is a parameter that will be optimized in the next steps. Consequently, the dimensions of each reluctance in the model should be a function of Rdeg. The reluctances are either located in the stress-free area, the degraded area or the welded area. The location of each reluctance will determine the permeability and core loss components of each reluctance. In the stress-free area, the reluctances have magnetic properties as defined by Equations (2), (5), (6) and (7). In the degraded zone, the reluctances have alternative magnetic properties which are defined by Equations (2), (6), (8) and (9). The reluctance which represents the welded zone is defined by the same equations as the degraded zone, except for the classical loss component, which is defined by Equation (10).

3.1.2. Magnetic Measurements

First, the single-valued BH-curve should be measured on an unwelded core. This curve defines the permeability of the reluctances in the stress-free zone. Secondly, the grade-dependent loss coefficients Ch, Ccl and Cex need to be identified using core loss measurements at different frequencies and peak flux densities on an unwelded core. Any fitting algorithm can be used, as long as the accuracy of the fitted model is sufficient. Once these coefficients are identified, they will represent the loss coefficients of the reluctances in the stress-free zone. Finally, core loss measurements should be done on the welded stack at several frequencies and peak flux densities. These measurements will be used for fitting the MEC model parameters.

3.1.3. Model Parameter Fitting

With the parametrized MEC model and welded core measurements available, the model parameters σ2, σ3 and Rdeg can be fitted using an optimization algorithm. The initial parameter values for σ2 and σ3 are set to 0 MPa. The initial value for Rdeg is set to Rweld + 1 mm, because Rdeg can never be smaller than or equal to Rweld. Then, the optimization is started using a constrained least squares algorithm. Once these parameters are optimized, the calculated core losses should correspond to the measured core losses on the welded core. Additionally, the model parameters give a coarse estimation of the stress near the weld seam and the size of the area which is damaged by the welding heat. The reluctances with their respective magnetic properties can now be copied into a full MEC machine model.

4. Case Study: Welded Stator Core Measurements

As a proof of concept, core losses were measured at different peak flux densities and excitation frequencies on a stack of M270-35A stator laminations before and after welding. The elastic limit of this material is 450 MPa and the inner and outer yoke diameter of the stator laminations was 147.6 mm and 193.9 mm, respectively. Consequently, the yoke section width is 23.15 mm. The laminations were laser cut and afterward stress relief annealed. This way, the cutting effect is reduced and a more accurate modeling of the welding effects is made possible. A stack of 15 laminations was welded on 4 equidistant locations around the edge of the yoke using tungsten inert gas (TIG) welding (compatible with assumption 5) with welding current 90 A at a welding speed of 4 mm/s. The weld seam radius Rweld was measured to be 3 mm. Global measurements of magnetic field strength Hs and magnetic flux density Ba were performed for a range of frequencies and peak flux densities similar as described in [17] using the setup displayed in Figure 8. The core losses are:
P c o r e = 1 T 0 T H s d B a d t d t
Core loss measurements were done on the unwelded and welded (with four seams) core over a range of frequencies (10, 50 and 200 Hz) and peak flux densities (0.1 T to 1.3 T in steps of 0.1 T). The measured losses in the unwelded stack were used to fit the parameters Ch, Ccl and Cexc of the undamaged material using Equation (4). The identification of these loss coefficients was made using the MATLAB surface fitting toolbox by plotting the measured and modeled core loss against the frequency and peak flux density.
Afterward, the unknown model parameters (Rdeg, σ2 and σ3) are fitted using a combination of the MATLAB genetic algorithm and fmincon toolbox. In this section, the measurements will be compared to the losses calculated by the fitted model and a validation of the model will be presented.
The identified core loss coefficients of undamaged material are Ch = 0.0126, Ccl = 1.195 × 10−5 and Cexc = 0.0011. After fitting the unknown model parameters, it was found that for a compressive stress state of σ2 = 20 MPa in region 2 and Rdeg =14 mm, the simulations correspond accurately with the measurements. The fitted stress state in the weld seam σ3 was 330 MPa tensile stress; however, it must be noted that this factor appeared to play a minor role in the loss calculation due to the comparatively small size of Rweld in the total width of 23 mm. Figure 9 displays the measured and simulated specific iron losses for different frequencies and peak flux densities for the welded and unwelded case. This graph shows a good correspondence between the fitted model and the measurements. The deviation between the measured and simulated values for each case can partially be explained by the noise-sensitive nature of the measurement setup. Moreover, the relation between the V0 and Bp was neglected in the model in order to simplify the stress dependency of V0. This could lead to inaccuracies when the material is saturated. Therefore, measurements and fitting were conducted for Bp-values up to 1.3 T. The eddy current losses at 1 T and 50 Hz were 15.8 W/kg in the welded zone and approximately 0.21 W/kg in the degraded and undamaged zone. The total core losses at 1 T and 50 Hz were 16.3 W/kg in the welded zone, 1.84 W/kg (on average) in the stress-free zone and 2.14 W/kg (on average) in the degraded zone.
In Figure 10, the model (which was fitted with loss measurements at 10, 50 and 200 Hz on a 4-weld seam stack) was validated in case the number of weld seams doubled to 8 at 100 Hz. Here, measurements on a welded stack with eight weld seams are compared to model calculations when 8-weld sections are simulated. This shows that the model accurately predicts the effect of core losses for varying frequencies and varying amounts of weld seams.

5. Discussion

In Figure 9, the core losses show an exponential increase with increasing peak flux density and frequency, which is in correspondence with Equation (4). When the laminations are welded, this exponential trend remains but appears to have shifted upwards. The total core losses increase by about 25% at 50 Hz and 1 T when four welding seams are added. In Table 3, this result is compared with other results found in the literature on similar welded cores. From this comparison, it is clear that the added core losses strongly depend on each specific case. Several specific parameters such as the steel grade, core geometry, welding speed and amount of welding seams all have an impact on the added losses.
The results obtained from the simulations demonstrate that the local eddy current losses in the welded zone are significantly higher than the other loss components. In this section, the effect of the weld radius Rweld on the total core losses is simulated and compared to the results presented in [7]. In this work, Wang and Zhang performed FE-simulations and analytical calculations on welded laminations where the total width was 5 mm and the sheet thickness was 0.27 mm. These results demonstrate that the eddy current losses increase exponentially when the weld radius increases from 0.8 mm to 1.3 mm at 50 Hz excitation frequency. With the optimized fitting parameters from the previous section kept constant, the effect of increasing the weld radius (from 3.7 mm to 6 mm, in order to keep the ratio of the weld radius over the total section width equal for both cases) on the total core losses is simulated and shown in Figure 11. This figure shows that the proposed stress-dependent MEC model adequately incorporates the effect of increasing core losses due to increasing weld radius, similar to the results in [7].

6. Conclusions

This paper presents a parametrized welding model that may be included in a complete MEC model of an electrical machine. The model accounts for local residual stress near the welded region and for local eddy currents in the weld seam. Although the model simplifies the welding geometry and magnetic sensitivity to stress and assumes uniform residual stress states in different sub-regions, it provides an interesting tool for studying the prominent effects of welding and allows for a qualitative study of the effect of different process parameters on the additional iron losses due to welding. When the model parameters are fitted correctly based on experimental observations on welded ring cores, it is possible to calculate additional core losses due to welding at varying amounts of welding seams, excitation frequencies and weld bead radii. The fitted parameters also give an estimation of the size of the degraded zone and the residual stress present in the area around the welded zone. The calculation time is low because the model does not require FE-simulations. In the case study, it was demonstrated that the core losses increase by approximately 25% when four welding joints are applied. The application of the model on cores of M270-35A grade steel has demonstrated that with a limited number of magnetic measurements on a welded and unwelded core, the model can be fully parametrized. The model was successfully validated on a core with eight welding seams at 100 Hz.

Author Contributions

Funding acquisition, L.D. and P.S.; Resources, M.N.I., L.D. and P.S.; Supervision, L.D. and P.S.; Writing—original draft, A.D.; Writing—review and editing, M.N.I. and P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Vlaio-NWO-project IWT.150457.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bourchas, K.; Stening, A.; Soulard, J.; Broddefalk, A.; Lindenmo, M.; Dahlén, M.; Gyllensten, F. Quantifying Effects of Cutting and Welding on Magnetic Properties of Electrical Steels. IEEE Trans. Ind. Appl. 2017, 53, 4269–4278. [Google Scholar] [CrossRef] [Green Version]
  2. Zhang, Y.; Wang, H.; Chen, K.; Li, S. Comparison of laser and TIG welding of laminated electrical steels. J. Mater. Process. Technol. 2017, 247, 55–63. [Google Scholar] [CrossRef]
  3. Schade, T.; Ramsayer, R.M.; Bergmann, J.P. Laser welding of electrical steel stacks investigation of the weldability. In Proceedings of the 2014 4th International, Proceedings of the Electric Drives Production Conference (EDPC), Nuremberg, Germany, 30 September–1 October 2014. [Google Scholar]
  4. Vourna, P.; Ktena, A.; Tsakiridis, P.E.; Hristoforou, E. A novel approach of accurately evaluating residual stress and microstructure of welded electrical steels. NDT E Int. 2015, 71, 33–42. [Google Scholar] [CrossRef]
  5. Leuning, N.; Steentjes, S.; Schulte, M.; Bleck, W.; Hameyer, K. Effect of elastic and plastic tensile mechanical loading on the magnetic properties of NGO electrical steel. J. Magn. Magn. Mater. 2016, 417, 42–48. [Google Scholar] [CrossRef]
  6. Krings, A.; Nategh, S.; Wallmark, O.; Soulard, J. Influence of the Welding Process on the Performance of Slotless PM Motors with SiFe and NiFe Stator Laminations. IEEE Trans. Ind. Appl. 2013, 50, 296–306. [Google Scholar] [CrossRef]
  7. Wang, H.; Zhang, Y. Modeling of Eddy-Current Losses of Welded Laminated Electrical Steels. IEEE Trans. Ind. Electron. 2016, 64, 2992–3000. [Google Scholar] [CrossRef]
  8. Bertotti, G. General properties of power losses in soft ferromagnetic materials. IEEE Trans. Magn. 1988, 24, 621–630. [Google Scholar] [CrossRef]
  9. Permiakov, V. 1D and 2D Magnetization in Electrical Steels under Uniaxial Stress. Doctoral dissertation, Ghent University, 2005. Available online: https://users.ugent.be/~ldupre/PhD_viatcheslav_permiakov.pdf (accessed on 6 November 2022).
  10. Pulnikov, A. Modification of magnetic properties of non oriented electrical steels by the production of electromagnetic devices. Doctoral dissertation, Ghent University, 2004. Available online: https://biblio.ugent.be/publication/8597986 (accessed on 6 November 2022).
  11. Aydin, U.; Rasilo, P.; Martin, F.; Belahcen, A.; Daniel, L.; Haavisto, A.; Arkkio, A. Effect of multi-axial stress on iron losses of electrical steel sheets. J. Magn. Magn. Mater. 2018, 469, 19–27. [Google Scholar] [CrossRef]
  12. Daniel, L.; Hubert, O. Equivalent Stress Criteria for the Effect of Stress on Magnetic Behavior. IEEE Trans. Magn. 2010, 46, 3089–3092. [Google Scholar] [CrossRef]
  13. LoBue, M.; Sasso, C.; Basso, V.; Fiorillo, F.; Bertotti, G. Power losses and magnetization process in Fe–Si non-oriented steels under tensile and compressive stress. J. Magn. Magn. Mater. 2000, 215-216, 124–126. [Google Scholar] [CrossRef]
  14. Permiakov, V.; Dupré, L.; Pulnikov, A.; Melkebeek, J. Loss separation and parameters for hysteresis modelling under compressive and tensile stresses. J. Magn. Magn. Mater. 2004, 272, E553–E554. [Google Scholar] [CrossRef]
  15. Karthaus, J.; Steentjes, S.; Leuning, N.; Hameyer, K. Effect of mechanical stress on different iron loss components up to high frequencies and magnetic flux densities. COMPEL-Int. J. Comput. Math. Electr. Electron. Eng. 2017, 36, 580–592. [Google Scholar] [CrossRef]
  16. Pulnikov, A.; Decocker, R.; Permiakov, V.; Dupré, L.; Vandevelde, L.; Petrov, R.; Melkebeek, J.; Houbaert, Y.; Gyselinck, J.; Wisselink, H. The relation between the magnetostriction and the hysteresis losses in the non-oriented electrical steels. J. Magn. Magn. Mater. 2005, 290, 1454–1456. [Google Scholar] [CrossRef]
  17. Abdallh, A.A.-E.; Sergeant, P.; Crevecoeur, G.; Vandenbossche, L.; Dupré, L.; Sablik, M. Magnetic material identification in geometries with non-uniform electromagnetic fields using global and local magnetic measurements. IEEE Trans. Magn. 2009, 45, 4157–4160. [Google Scholar] [CrossRef]
  18. Schoppa, A.; Schneider, J.; Wuppermann, C.-D.; Bakon, T. Influence of welding and sticking of laminations on the magnetic properties of non-oriented electrical steels. J. Magn. Magn. Mater. 2003, 254–255, 367–369. [Google Scholar] [CrossRef]
Figure 1. Welding seam on outer edge of stator laminations.
Figure 1. Welding seam on outer edge of stator laminations.
Machines 10 01153 g001
Figure 2. Top view and cross section view of a welded stack of laminations. Welded and degraded region are simplified to a rectangular area.
Figure 2. Top view and cross section view of a welded stack of laminations. Welded and degraded region are simplified to a rectangular area.
Machines 10 01153 g002
Figure 3. MEC model for a ring core with one weld seam where the flux perpendicular to the cross-section travels through parallel branches of healthy, degraded and welded material. The indices of the reluctances indicate different stress-dependent material properties corresponding to the regions described in Figure 2.
Figure 3. MEC model for a ring core with one weld seam where the flux perpendicular to the cross-section travels through parallel branches of healthy, degraded and welded material. The indices of the reluctances indicate different stress-dependent material properties corresponding to the regions described in Figure 2.
Machines 10 01153 g003
Figure 4. Approximated uni-axial relation between the stress correction factor of relative permeability Kµ and the present stress state according to [9,10]. A factor higher than 1 means increased permeability due to stress. The relative permeability peaks at 140% of the stress-free permeability at a stress state of 10% of the material elastic limit, denoted by e/10 in this figure. This figure is valid for a material with elastic limit e = 350 MPa. For materials with a different elastic limit, the graph should be scaled accordingly.
Figure 4. Approximated uni-axial relation between the stress correction factor of relative permeability Kµ and the present stress state according to [9,10]. A factor higher than 1 means increased permeability due to stress. The relative permeability peaks at 140% of the stress-free permeability at a stress state of 10% of the material elastic limit, denoted by e/10 in this figure. This figure is valid for a material with elastic limit e = 350 MPa. For materials with a different elastic limit, the graph should be scaled accordingly.
Machines 10 01153 g004
Figure 5. Approximated uni-axial relation between the stress correction factor of hysteresis loss coefficient Khyst and the stress state [14]. Khyst reaches a minimum of 0.3 at a stress of 10% of the material elastic limit e. This figure is valid for a material with elastic limit = 350 MPa.
Figure 5. Approximated uni-axial relation between the stress correction factor of hysteresis loss coefficient Khyst and the stress state [14]. Khyst reaches a minimum of 0.3 at a stress of 10% of the material elastic limit e. This figure is valid for a material with elastic limit = 350 MPa.
Machines 10 01153 g005
Figure 6. Approximated uni-axial relation between the stress correction factor of V0-parameter Kexc and the stress state [14]. Kexc reaches a minimum of 10% of the stress-free permeability at a stress state of 10% of the material elastic limit e. This figure is valid for a material with elastic limit = 350 MPa.
Figure 6. Approximated uni-axial relation between the stress correction factor of V0-parameter Kexc and the stress state [14]. Kexc reaches a minimum of 10% of the stress-free permeability at a stress state of 10% of the material elastic limit e. This figure is valid for a material with elastic limit = 350 MPa.
Machines 10 01153 g006
Figure 7. Flowchart describing the methodology for applying the MEC model.
Figure 7. Flowchart describing the methodology for applying the MEC model.
Machines 10 01153 g007
Figure 8. Setup used for obtaining loss measurements: amplifier (left), welded stator core (center) and data acquisition system (right).
Figure 8. Setup used for obtaining loss measurements: amplifier (left), welded stator core (center) and data acquisition system (right).
Machines 10 01153 g008
Figure 9. Comparison of the simulated specific losses and measured specific losses as a function of peak flux density for 50 Hz and 200 Hz excitation frequencies. The peak flux density values are averaged over the magnetic cross section, as they are known to vary locally near the welded region. The model parameters were fitted using these data.
Figure 9. Comparison of the simulated specific losses and measured specific losses as a function of peak flux density for 50 Hz and 200 Hz excitation frequencies. The peak flux density values are averaged over the magnetic cross section, as they are known to vary locally near the welded region. The model parameters were fitted using these data.
Machines 10 01153 g009
Figure 10. Comparison of the simulated specific losses and measured specific losses as a function of peak flux density for 100 Hz and 8 welding seams. The simulation curve is obtained by fitting the measured data at different frequencies of a stack with 4 weld seams.
Figure 10. Comparison of the simulated specific losses and measured specific losses as a function of peak flux density for 100 Hz and 8 welding seams. The simulation curve is obtained by fitting the measured data at different frequencies of a stack with 4 weld seams.
Machines 10 01153 g010
Figure 11. Effect of increasing weld radius on the total core losses when all other model parameters remain equal.
Figure 11. Effect of increasing weld radius on the total core losses when all other model parameters remain equal.
Machines 10 01153 g011
Table 1. Symbols commonly used in the model.
Table 1. Symbols commonly used in the model.
SymbolDescriptionUnit
RweldRadius of the welding seammm
RdegRadius of the degraded zonemm
σ1Stress in stress-free zone (always assumed 0 MPa)MPa
σ2Residual stress in degraded zoneMPa
σ3Residual stress in welded zoneMPa
σeqEquivalent stress calculated in each reluctance by Equation (1)MPa
ChCore loss coefficient for hysteresis loss
CclCore loss coefficient for classical loss
CexCore loss coefficient for excess loss
KµStress-dependent correction factor for magnetic permeability
KhStress-dependent correction factor for hysteresis loss component
KexStress-dependent correction factor for excess loss component
wtotAverage yoke widthmm
ltotAverage yoke circumferencemm
dLamination thicknessmm
BpPeak flux densityT
µrRelative permeability
lrLength of reluctancemm
ArCross section of reluctancemm2
fExcitation frequencyHz
ρElectrical resistivityΩm
Table 2. Model assumptions.
Table 2. Model assumptions.
1The weld seam area and its surrounding heat affected area are simplified as a rectangular area with dimensions specified in Figure 2.
2The curvature on the outer edge of the ring core is simplified to a straight line.
3A uniform isotropic stress state is assumed in each reluctance of the considered MEC, but stress anisotropy with respect to magnetization direction is considered.
4The skin effect is neglected which is a valid assumption when considering sufficiently low frequencies.
5The applied welding technique does not change the metallurgic composition of the weld seam. The electrical steel melts and bonds with adjacent laminations without any additional binding material during the welding process.
Table 3. Comparison of welded stator measurements with other examples found in literature.
Table 3. Comparison of welded stator measurements with other examples found in literature.
ReferenceSample GeometryAdded Core Loss Due to
Welding at Sinusoidal Flux Density of 50 Hz and Bp = 1 T
Case StudyRing Core, 4 weld seams25%
[1]Epstein Strips, 3 weld seams10%
[2]Ring Core, 6 weld seams42.5%
[6]Ring Core, 8 weld seams22%
[18]Ring Core, 2 weld seams10%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Daem, A.; Ibrahim, M.N.; Sergeant, P.; Dupré, L. Stress-Dependent Magnetic Equivalent Circuit for Modeling Welding Effects in Electrical Steel Laminations. Machines 2022, 10, 1153. https://doi.org/10.3390/machines10121153

AMA Style

Daem A, Ibrahim MN, Sergeant P, Dupré L. Stress-Dependent Magnetic Equivalent Circuit for Modeling Welding Effects in Electrical Steel Laminations. Machines. 2022; 10(12):1153. https://doi.org/10.3390/machines10121153

Chicago/Turabian Style

Daem, Andries, Mohamed N. Ibrahim, Peter Sergeant, and Luc Dupré. 2022. "Stress-Dependent Magnetic Equivalent Circuit for Modeling Welding Effects in Electrical Steel Laminations" Machines 10, no. 12: 1153. https://doi.org/10.3390/machines10121153

APA Style

Daem, A., Ibrahim, M. N., Sergeant, P., & Dupré, L. (2022). Stress-Dependent Magnetic Equivalent Circuit for Modeling Welding Effects in Electrical Steel Laminations. Machines, 10(12), 1153. https://doi.org/10.3390/machines10121153

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