Next Article in Journal
Cr2P2O7 as a Novel Anode Material for Sodium and Lithium Storage
Previous Article in Journal
Experimental Investigation and ANFIS-Based Modelling During Machining of EN31 Alloy Steel
Previous Article in Special Issue
Investigating the Effect of Zn Ferrite Nanoparticles on the Thermomechanical, Dielectric and Magnetic Properties of Polymer Nanocomposites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Anisotropic Vector Hysteresis Simulation of Soft Magnetic Composite Materials Based on a Hybrid Algorithm of PSO–Powell

1
Department of Electrical and Electronic Engineering, North China Electric Power University, Baoding 071003, China
2
Hebei Provincial Key Laboratory of Electromagnetic & Structural Performance of Power Transmission and Transformation Equipment, Baoding 071056, China
3
State Key Lab of Reliability and Intelligence of Electrical Equipment, Hebei University of Technology, Tianjin 300130, China
*
Author to whom correspondence should be addressed.
Materials 2020, 13(14), 3138; https://doi.org/10.3390/ma13143138
Submission received: 29 May 2020 / Revised: 9 July 2020 / Accepted: 12 July 2020 / Published: 14 July 2020

Abstract

:
To simulate the anisotropic hysteresis characteristics of soft magnetic composite (SMC) materials accurately, an improved vector hysteresis model was proposed and utilized to adjust the shape of hysteresis curves by introducing two parameters. These two parameters are correlated with the amplitude of the vector Everett function and the projection of magnetic flux density along different directions. An experimental platform was built to measure the two-dimensional (2-D) magnetic properties of the SMC material under rotational magnetizations. The scalar and vector Everett functions were constructed by the measured limiting hysteresis loops. A hybrid optimization strategy based on the particle swarm optimization (PSO) and Powell technique was proposed to identify the parameters of the improved model efficiently and precisely, which significantly improved the local optimization ability of the PSO algorithm. The simulated results strongly agree with the measured ones, and thus the effectiveness of the improved vector model and the parameter identification method proposed in this paper was verified.

1. Introduction

Soft magnetic composite (SMC) materials have been extensively used in motors and power electronic devices due to their unique electromagnetic properties such as magnetic and thermal isotropy, diversified processing shapes, and low eddy current loss at medium and higher frequencies [1]. As a kind of complex characteristic inherent in SMC materials, hysteresis has a great impact on the optimal design and analysis of electrical equipment, so accurate measurement and simulation of the hysteresis properties of these materials is critical to their research and application [2].
Traditional scalar simulation of the magnetic performance enables precise results under alternating magnetic excitations. However, the electrical devices usually operate under alternating and rotational excitations in practice, and the magnetic flux density B and magnetic field strength H do not always align in the same direction. The iron loss caused by the rotational excitations is greater than that induced by the unidirectional alternating fields [3,4,5]. Therefore, it is of great importance to measure and model the vector hysteresis properties of SMC materials under rotational excitations [6].
Different vector hysteresis models are presented in previous studies to obtain the rotational magnetic characteristics by using some newly developed measuring methods [7]. Stoner and Wohlfarth proposed a vector hysteresis model which was designed as an ensemble of single-domain and uniaxial particles. Nevertheless, it was impossible to fit the asymmetric hysteresis loops and accomplish the parameter identification completely [8]. In view of the defects of the model mentioned above, Mayergoyz presented a classical vector Preisach model to simulate the isotropic vector hysteresis characteristics of magnetic materials [9]. However, the anisotropy magnetic properties under rotational excitations are found even at a relatively low frequency [10,11], so the classical model is no longer applicable in most cases, leading to imperative demands for further improvement. Miklós Kuczmann proposed an improved vector Preisach model to simulate the slight anisotropic behavior of isotropic materials, whereas it indicated great errors at low amplitudes of the magnetic flux density [12]. In addition, M. Enokizono and N. Soda presented the E&S vector hysteresis model considering the anisotropic properties under rotational magnetization. It requires large amounts of experimental data to obtain numerous parameters, which severely limits its application [13].
In this paper, an experimental platform for the measurement of vector magnetic properties was employed to obtain the 2-D magnetic properties of the SMC material under rotational excitations [14,15,16]. By introducing the correlation parameters, the classic vector hysteresis model was modified to simulate the vector hysteresis curves of the SMC material, considering the anisotropy property. Based on the measured limiting hysteresis loops in two orthogonal directions, the first-order reversal curves (FORCs) were numerically generated to construct the scalar and vector Everett functions. A hybrid optimization strategy combining the particle swarm optimization (PSO) algorithm and the Powell technique was proposed for parameter identification in the improved model with good efficiency and precision, which significantly improved the local optimization ability of the PSO algorithm. The validity of the proposed method was proved by comparing the measured results with the simulated ones.

2. Material and Methods

2.1. Material

The SMC material, SOMALOYTM 500 (developed by Höganäs in Höganäs, Sweden), was employed, which has a maximum magnetic flux density of 2.1 T at a magnetic field strength of 100 kA/m, residual magnetization of 0.25 T, a coercive magnetic field strength of 250 A/m, conductivity of 30 μS/m, a coefficient of thermal conductivity of 17 W/(m∙K) and an initial relative magnetic permeability of 130.

2.2. Vector Hysteresis Measurement of SMC Material

The vector magnetic property testing system consisted of a main measuring device, a power amplifier, a magnetic sensing system and a signal-processing unit. The main measuring devices are shown in Figure 1, including the “C-type” core-poles and core-yokes, excitation windings, and three B-H composite sensing coils.
The measuring system was used to get the 2-D magnetic properties of the material. During the experiment, a cubic specimen made of SMC with a size of 22 × 22 × 22 mm3 was fixed in a sensing box mounted with B-H sensing coils on six sides to measure the corresponding B and H.
Figure 2 is the structure diagram of the measuring system. The signal generating system was used as the computing core to provide the analog excitation signals for three high power amplifiers. Then the excitation voltage signals were selected by the impedance matching loop and the resonant circuit to obtain a large excitation current, which acted as the input of the main excitation magnetic circuit so as to generate a rotational magnetic field inside the measured specimen. The sensing coils on the specimen can pick up weak electrical signals, and then these signals are input into the signal processing circuit and collected by a LabView embedded control system. The closed-loop feedback of the system is realized in such way.
In order to ensure the circular trajectory of the 2-D magnetic flux density B, a circular rotational excitation at a low frequency (f = 5 Hz) was set on the x-o-y plane. The locus of the vector B and H with the increase of excitation current are demonstrated in Figure 3, and the locus of H exhibits a special feature in an asymmetric and irregular shape due to the slight anisotropic property.
The measured concentric hysteresis loops along two orthogonal directions, i.e., x direction and y direction, are shown in Figure 4. The hysteresis curves at the saturated magnetic flux density are the limiting hysteresis loops required for identification of the Preisach model.

2.3. Hysteresis Modeling Based on Improved Preisach Model

2.3.1. Improved Vector Preisach Model Based on Classical Model

The inverse form of the classical vector Preisach model, which predicts the magnetic field strength H from the magnetic flux density B, is expressed as follows [9],
H ( t ) = π / 2 π / 2 e φ ( H φ ( B φ ) ) d φ
where Hφ(Bφ) is the scalar magnetic field strength of H(t) in the direction eφ.
In the numerical calculation, the angle φ ϵ [−π/2, π/2] is evenly divided into n cells:
φ i = π 2 + i 1 n π
where i = 1,…,n and n is the number of directions.
The vector magnetic field strength H(t) in Equation (1) can be expressed as the sum of the scalar magnetic field strength Hφi(Bφi) along all the n directions, described as follows,
H ( t ) = i = 1 n e φ i ( H φ i ( B φ i ) ) Δ φ
H φ i ( B φ i ) = α β v ( α , β ) γ α β B φ i d α d β
where ∆φ = π/n. γαβ is the simplest hysteresis operator with α and β corresponding to “up” and “down” switching values of the input, respectively. v is the vector Preisach distribution function.
The magnetic flux density in the direction φi can be calculated from the two components of the magnetic flux density:
B φ i = B x cos φ i + B y sin φ i
As shown in Equation (4), the key to identifying the vector model is to obtain the vector distribution function v. There are two main difficulties in identifying v directly. On one hand, it is time-consuming to evaluate the double integral in Equation (4) numerically. On the other hand, the determination of the distribution function v requires differentiations of experimentally obtained data, which may amplify the inherent experimental errors significantly [9].
In order to eliminate the double integral of α and β, the discrete Everett function E [9] is commonly used to identify the vector hysteresis model numerically. The specific expression of E is as follows,
E ( α , β ) = T ( α , β ) v ( α , β ) d α d β
where T(α,β) is the limit triangle surrounded by α and β in the Preisach plane.
According to Equation (6), Hφi(Bφi) can be expressed by E in the following equation:
H φ i ( B φ i ) = E ( B 0 , b 0 ) + 2 k = 1 n φ i ( t ) [ E ( B φ i , k , b φ i , k 1 ) E ( B φ i , k , b φ i , k ) ]
where B0 and b0 are the positive and negative saturation values of B. Bφi,k and bφi,k−1 represent a series of maximum and minimum values of Bφi. nφi(t) is the number of reversal points of the first-order reversal curves (FORCs) [17,18,19] in the direction φi.
The function E is related to the scalar Everett function F, as described in the following equation:
F ( α , β ) = π / 2 π / 2 cos φ E ( α cos φ , β cos φ ) d φ
where function F can be generated from the measured data by numerical methods.
The classical vector Preisach hysteresis model given above implies that for a uniformly and rotationally applied B, the locus of H is circular [20], as shown in Figure 3a. However, the locus of H exhibits a bit petal-like shape because of the slight magnetic anisotropy property of the SMC material itself, as shown in Figure 3b. It is obvious that the classical vector model is no longer applicable. In order to circumvent the limitations of the classical model above, an improved Preisach model was proposed by introducing two parameters w and z. These two parameters are connected to the amplitude of Everett function E and the projection of magnetic flux density along different directions.
H ( t ) = π / 2 π / 2 e φ ( H φ ( B φ ) ) d φ i = 1 n e φ i ( H φ i ( B φ i ) ) Δ φ
B φ i = B x sign ( cos φ i ) | cos φ i | w + B y sign ( sin φ i ) | sin φ i | w
F ( α , β ) = π / 2 π / 2 cos φ E ( α cos w φ , β cos w φ ) z ( B m ) d φ
where Bm is the magnitude of B.
From Equations (10) and (11), it can be seen that the parameter w generalizes the projection of B and controls the projection of H along each discrete direction. It enables the slight anisotropic property of the model. The parameter z was used to adjust the function E and the amplitude of H under different Bm to get more accurate fitted results. When 0 < w < 1, the improved model can simulate slight anisotropic hysteresis characteristics of the material, thus making the locus of H an elliptical or a petal-like shape. In particular, it is in accordance with the classic vector model as w = 1.
The two orthogonal components of the magnetic field strength are given by
H x = i = 1 n H φ i cos φ i
H y = i = 1 n H φ i sin φ i

2.3.2. Identification Procedure of Improved Vector Preisach Model

According to Equation (11), the numerical identification of scalar Everett function F [21] can be represented by function E as follows,
F ( α , β ) = i = 1 n cos φ i E ( α cos w φ i , β cos w φ i ) z ( B m ) Δ φ
The function F in Equation (14) can be constructed numerically based on the measured limiting hysteresis loops in Figure 4. Their descending branches are used for the generation of the FORCs by the numerical method proposed by Dlala [22], as shown in Figure 5.
According to the identification procedure presented in [9], the scalar Everett functions can be obtained from the interpolation based on FORCs,
F ( b ¯ u , b ¯ v ) = 1 2 ( H forc ( b ¯ u , b ¯ v ) H forc ( b ¯ u ) )
where b ¯ u and b ¯ v respectively correspond to the discrete increasing and decreasing values of B on the FORCs. Hforc( b ¯ u , b ¯ v ) refers to the value of H on the FORCs, while Hforc( b ¯ u ) represents the value of H at the reversal points.
According to Equation (15), the corresponding scalar Everett functions Fx and Fy are depicted in Figure 6, respectively.
The scalar Everett function in the direction φ is expressed as the elliptic interpolation of Fx and Fy to approximate the smooth angular behavior,
F ( α , β , φ ) = F x 2 ( α , β ) cos 2 φ + F y 2 ( α , β ) sin 2 φ
Therefore, the scalar Everett function F under the circular rotational magnetization is represented as follows [21],
F ( α , β ) = Δ φ π ( F ( α , β , 0 ) + F ( α , β , π 2 ) + 2 j = 1 n 1 F ( α , β , φ j ) )
The generated scalar Everett function F is depicted in Figure 7a.
As shown in Equation (14), it should be noticed that the function E can be obtained from the scaler Everett function F at specific values of Bm, w, and z. Assuming that Bm = 1.398 T, w = 0.7 and z = 2, the corresponding vector Everett function is shown in Figure 7b.

2.3.3. Parameter Extraction of Improved Model Based on Hybrid Optimization Algorithm

The improved model can simulate the anisotropy property of SMC materials by introducing parameters w and z. The extraction of these two parameters is essential to the simulation accuracy. For full utilization of the strong global search ability of the random optimization algorithm as well as the fast local convergence of the deterministic optimization algorithm, a hybrid optimization strategy combining PSO with the Powell algorithm was presented to implement the parameter extraction.
The mean absolute percent error (MAPE) was used to evaluate the simulation accuracy of the improved vector hysteresis model. Thereafter, the parameter extraction of the improved model can be performed by searching for the minimum value of the objective function f given in Equation (18).
min f = 1 N i = 1 N | H cal ( i ) H mea ( i ) H mea ( i ) | × 100 %
where Hcal and Hmea refer to the calculated and measured magnetic field strength, respectively. N is the total amount of data.
During the initial iteration of the hybrid algorithm, the PSO algorithm was first used to perform a wide-ranging optimization to quickly lock the solution region. Specific steps are as follows.
(1)
Set the parameters of the PSO algorithm. The particle number N, the acceleration factors c1 and c2, the inertia factor ω, the maximum number of iterations T and the initial iteration number k are set to 5, 0.3, 0.3, 1, 100 and 1 respectively.
(2)
Initialize the population. The position X(wi0, zi0) of each initial particle i is generated randomly within wi0 ϵ [0.5,1], zi0 ϵ [1,1.5]; the range of the initial particles’ velocity Vi0 ϵ [Vmin, Vmax] is set to [0,0.3].
(3)
Calculate the objective function fi0 of each initial particle to obtain the historical optimal position Pi0 = Xi0 and the global optimal position Pg0 = P(wg0, zg0).
(4)
Update the particle velocity Vik and position Xik by Pik−1 and Pgk−1
{ V i k = ω V i k 1 + c 1 r 1 ( P i k 1 X i k 1 ) + c 2 r 2 ( P g k 1 X i k 1 ) X i k = X i k 1 + V i k
where r1 and r2 are generated randomly within the interval of [0,1].
(5)
Evaluate the fik of each particle and update the historical optimal position Pik of each particle as well as global optimal objective position Pgk as follows
{ f i k = f ( P i k ) f g k = f ( P g k ) = min i = 1 , , N { f i k } = min i = 1 , , N { f ( P i k ) }
where fgk is the global optimal objective function at the kth iteration.
(6)
Determine whether the switching criteria as Inequation (21) is satisfied. If satisfied, the current optimal solution Pgk= P(wgk, zgk), and the corresponding objective function fgk, are transferred to the Powell algorithm and the calculation process is ended. Otherwise, set k = k + 1 and repeat from step (3).
n = 1 n 0 | f g k + 1 n f g k n | f g k < ε
where ε = 0.01. n0 is the number of consecutive iterations, and was set to 10 in this paper.
After receiving the optimal solution provided by the PSO, the Powell algorithm which can converge to the optimal solution efficiently was utilized to optimize the parameters of the vector Preisach model. The specific process is as follows:
(1)
Initial the basic point of the Powell algorithm: x0(1) = x(0) = x(w(0), z(0)) = P(wgk, zgk).
(2)
Set the parameters of the Powell algorithm: the iteration accuracy e, the initial direction S1(1), S2(1), and the initial iteration number t, are set to 0.001, (1,0), (0,1) and 1, respectively.
(3)
Basic search: start from x0(t) and do a 1-D search along S1(t) and S2(t) to obtain the extreme points x1(t) and x2(t) for f.
(4)
Accelerated search: start from x0(t), perform a 1-D search along the conjugate direction S(t) = x2(t) - x0(t) to get the extreme point x3(t).
(5)
Determine whether the termination condition as Inequality (22) is met. If satisfied, the current optimal solution x* = x3(t) = x(w3(t), z3(t)), and the corresponding optimal value f(x*), are obtained. Otherwise, go to step (6).
x 3 ( t ) x 0 ( t ) e
(6)
Calculate the maximum drop Δ m ( t ) of f and the corresponding direction Sm(t) as follows.
{ Δ m ( t ) = max i = 1 , 2 { Δ i ( t ) } = max i = 1 , 2 { f ( x i - 1 ( t ) ) f ( x i ( t ) ) } S m ( t ) = x m ( t ) x m - 1 ( t )
Calculate the mapping point xmap(t) = 2x2(t)x0(t) along direction S(t) and set f1 = f(x0(t)), f2 = f(x2(t)), f3 = f(xmap(t)). Update the initial point x0(t + 1) = x3(t) and the search direction Sm(t) = S(t) if the Powell condition is satisfied as Inequality (24), and then repeat from step (3). Otherwise, go to step (7).
{ f 3 < f 1 ( f 1 2 f 2 + f 3 ) ( f 1 f 2 Δ m ( t ) ) 2 < 0.5 Δ m ( t ) ( f 1 f 3 ) 2
(7)
Update the initial point: x0(t + 1) = x2(t) if f 2 < f 3 . Otherwise, update the initial point: x0(t + 1) = xmap(t). Set t = t + 1 and repeat from step (3).
The specific calculation flowchart of the parameter extraction based on hybrid algorithm of PSO–Powell is demonstrated in Figure 8.
By using the PSO–Powell hybrid algorithm, the numerical simulation of the anisotropic vector hysteresis properties under circular rotational magnetization was completed by the improved vector hysteresis model, and the detailed procedures are given in Figure 9.

3. Results and Discussion

The PSO algorithm and the hybrid optimization strategy of PSO–Powell, respectively, give the same initial values of w and z to optimize the two parameters, so we can compare the convergence performance of the two algorithms intuitively. The values of the maximum iteration number are set appropriately for the two algorithms, and the termination criteria of the hybrid algorithm is determined by the convergence condition of Powell method.
Based on the PSO algorithm, the variation of MAPE with the number of iterations at Bm = 1.398 T is demonstrated in Figure 10, which shows strong global optimization ability but worse local search ability. Although it takes only four iterations to reduce MAPE from 21.3698% to 7.3943% in the beginning, the solutions are trapped in local optimum when the values of MAPE vary little. The aforementioned defect indicates that the convergence performance of the PSO algorithm tends to deteriorate when approaching the solution region, which makes it difficult to extract the optimal parameters accurately.
Table 1 lists the parameters w and z, corresponding values of MAPE, as well as the bounds of relative error extracted by PSO algorithm at Bm = 0.178 T, 0.423 T, 0.704 T, 1.053 T and 1.398 T. As shown in Figure 11, the comparison between the measured hysteresis curves of vector H and the simulated results obtained from the extracted parameters were made, which showed considerable discrepancy. This brings a significant obstacle to the magnetic hysteresis property simulation of the SMC material.
In contrast to the convergence performance of the PSO algorithm in Figure 10, Figure 12 as below shows the corresponding variation of MAPE with the iteration steps based on the PSO-Powell hybrid algorithm.
The trend of MAPE variation shows that the PSO algorithm meets the switching criteria automatically when it reaches the 12th iteration, thereafter the Powell algorithm is started with the previously optimized results as the initial solution. Then it takes only four iterations to achieve the convergence. These results show that the hybrid algorithm of PSO–Powell exhibits a faster convergence speed than the PSO algorithm.
In contrast to the values in Table 1, the values of the parameters and corresponding error extracted by the hybrid algorithm are listed in Table 2. The values of MAPE and bounds of relative error are, respectively, adjusted below 5.9571% and 10.7819%, which shows a significant reduction in the error of the improved model and a higher accuracy than that of the PSO algorithm.
Based on the extracted parameters in Table 2, the comparison between the measured hysteresis curves of vector H and simulated results optimized by the hybrid algorithm was made and is depicted in Figure 13. The simulated results strongly agree with the measured ones when Bm varies within a wide range, which can reflect the slight anisotropy hysteresis properties of the specimen. The effectiveness and accuracy of the improved vector model connected with the hybrid algorithm of PSO–Powell were verified.
By applying the hybrid algorithm, the slight error of the hysteresis loops simulated by the extracted parameters tended to increase at high Bm (more than 1 T), which was mainly attributed to the limitations of the proposed model. The model makes some improvements on the isotropic vector Preisach hysteresis model. It is noticed that the motion of the domain wall under the circular rotational excitation was more complicated with the increase of magnetic flux density (not only translation, but also rotation). The anisotropy properties of the material were more obvious, which makes accurate prediction by the proposed isotropic vector model difficult. Further improvements based on the classical vector hysteresis model, or application of the anisotropic vector model, are required for better accuracy.

4. Conclusions

The rotational magnetic characteristics of the SMC material are measured based on the experimental platform. Considering the anisotropic property of the material, an improved vector hysteresis model was put forward by introducing two parameters correlated with the amplitude of vector Everett function and the projection of magnetic flux density along different directions, to improve the shape of hysteresis curves.
Benefiting from the strong global search ability of the random optimization algorithm, in collaboration with fast local convergence of the deterministic optimization algorithm, a hybrid optimization strategy of PSO–Powell was proposed to extract the model parameters precisely and efficiently.
The simulated magnetic hysteresis characteristics of the SMC material under circular rotational magnetic excitation were basically consistent with the measured ones, which verified the effectiveness of the parameter extraction process by the hybrid optimization strategy and the improved model proposed in this paper.

Author Contributions

Methodology, X.Z., H.X. and Z.D.; validation, H.X. and Z.D.; formal analysis, X.Z., H.X. and Z.D.; resources, Z.D., Y.L. and L.L.; data curation, X.Z. and Z.Z.; Writing—Original draft preparation, X.Z. and H.X.; Writing—Review and editing, Z.D.; visualization, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key Research and Development Program of China, grant number 2017YFB0902703, National Natural Science Foundation of China, grant number 51777073, Hebei Province Natural Science Foundation, grant number E2017502061, State Key Laboratory of Reliability and Intelligence of Electrical Equipment, Hebei University of Technology, grant number EERIKF2018011 and Fundamental Research Funds for the Central Universities, grant number 2019MS078.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhu, J.; Guo, Y.; Lin, Z.; Zhong, J. Measurement and modeling of SMC materials under vector magnetizations. In Proceedings of the 2005 International Conference on Electrical Machines and Systems, Nanjing, China, 27–29 September 2005; Volume 3, pp. 2354–2359. [Google Scholar]
  2. Guo, Y.; Zhu, J.; Watterson, P.A.; Wu, W. Development of a PM transverse flux motor with soft magnetic composite core. IEEE Trans. Energy Convers. 2006, 21, 426–434. [Google Scholar] [CrossRef]
  3. Coulson, M.A.; Slater, R.D.; Simpson, R.R.S. Representation of magnetic characteristic, including hysteresis, using Preisach’s theory. Proc. Inst. Electr. Eng. 1977, 124, 895–898. [Google Scholar] [CrossRef]
  4. Yang, Q.; Li, Y.; Zhao, Z.; Zhu, L.; Luo, Y.; Zhu, J. Design of a 3-D rotational magnetic properties measurement structure for soft magnetic materials. IEEE Trans. Appl. Supercond. 2014, 24, 1–4. [Google Scholar] [CrossRef]
  5. Moses, A.J. Importance of rotational losses in rotating machines and transformers. J. Mater. Eng. Perform. 1992, 1, 235–244. [Google Scholar] [CrossRef]
  6. Hernandez-Aramburo, C.A.; Green, T.C.; Smith, A.C. Estimating rotational iron losses in an induction machine. IEEE Trans. Magn. 2003, 39, 3527–3533. [Google Scholar] [CrossRef] [Green Version]
  7. Li, Y.; Yang, Q.; Zhu, J.; Lin, Z.; Guo, Y.; Sun, J. Research of Three-Dimensional Magnetic Reluctivity Tensor Based on Measurement of Magnetic Properties. IEEE Trans. Appl. Supercond. 2010, 20, 1932–1935. [Google Scholar]
  8. Stoner, E.; Wohlfarth, E. A mechanism of magnetic hysteresis in heterogeneous alloys. IEEE Trans. Magn. 1991, 27, 3475–3518. [Google Scholar] [CrossRef]
  9. Mayergoyz, I.D. Mathematical Model of Hysteresis and Their Applications; Academic Press: New York, NY, USA, 2003. [Google Scholar]
  10. Zhu, L.; Wu, W.; Xu, X.; Guo, Y.; Li, W.; Lu, K.; Koh, C. An Improved Anisotropic Vector Preisach Hysteresis Model Taking Account of Rotating Magnetic Fields. IEEE Trans. Magn. 2019, 55, 1–4. [Google Scholar] [CrossRef]
  11. Li, Y.; Zhu, J.; Yang, Q.; Lin, Z.W.; Wang, Y.; Xu, W. Magnetic properties measurement of soft magnetic composite materials over wide range of excitation frequency. IEEE Trans. Ind. Appl. 2012, 48, 88–97. [Google Scholar] [CrossRef]
  12. Kuczmann, M. Measurement and Simulation of Vector Hysteresis Characteristics. IEEE Trans. Magn. 2009, 45, 5188–5191. [Google Scholar] [CrossRef]
  13. Soda, N.; Enokizono, M. Utilization method of electrical steel sheets on stator of self-propelled rotary actuator. In Proceedings of the XXII International Conference on Electrical Machines (ICEM), Lausanne, Switzerland, 4–7 September 2016; pp. 918–923. [Google Scholar]
  14. Zhu, J.; Ramsden, V.S. Two dimensional measurement of magnetic field and core loss using a square specimen tester. IEEE Trans. Magn. 1993, 29, 2995–2997. [Google Scholar] [CrossRef]
  15. Guo, Y.; Zhu, J.; Zhong, J.; Lu, H.; Jin, J. Measurement and modeling of rotational core losses of soft magnetic materials used in electrical machines: A review. IEEE Trans. Magn. 2008, 44, 279–291. [Google Scholar]
  16. Sievert, J. The measurement of magnetic properties of electrical sheet steel—Survey on methods and situation of standards. J. Magn. Magn. Mater. 2000, 215, 647–651. [Google Scholar] [CrossRef]
  17. Dlala, E.; Saitz, J.; Arkkio, A. Inverted and forward Preisach models for numerical analysis of electromagnetic field problems. IEEE Trans. Magn. 2006, 42, 1963–1973. [Google Scholar] [CrossRef]
  18. Lin, J.; Chiang, M. Tracking Control of a magnetic shape memory actuator using an inverse Preisach model with modified fuzzy sliding mode control. Sensors 2016, 16, 1368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Naidu, S.R. Simulation of the hysteresis phenomenon using Preisach’s theory. IEE Proc. Part A 1990, 137, 73–79. [Google Scholar] [CrossRef]
  20. Dlala, E.; Belahcen, A.; Fonteyn, K.A.; Belkasim, M. Improving loss properties of the Mayergoyz vector hysteresis model. IEEE Trans. Magn. 2010, 46, 918–924. [Google Scholar] [CrossRef]
  21. Iványi, A.; Kuczmann, M. Vector hysteresis model based on neural network. COMPEL Int. J. Comput. Math. Electr. Electron. Eng. 2003, 22, 730–743. [Google Scholar]
  22. Dlala, E. Efficient Algorithms for inclusion of the Preisach hysteresis model in nonlinear finite-element methods. IEEE Trans. Magn. 2011, 47, 395–408. [Google Scholar] [CrossRef]
Figure 1. Experimental platform: (a) Main measuring system; (b) Diagrammatic sketch of magnetic circuit for measurement.
Figure 1. Experimental platform: (a) Main measuring system; (b) Diagrammatic sketch of magnetic circuit for measurement.
Materials 13 03138 g001
Figure 2. Structure block diagram of the measuring system of vector magnetic property.
Figure 2. Structure block diagram of the measuring system of vector magnetic property.
Materials 13 03138 g002
Figure 3. Measured hysteresis curves: (a) Locus of vector B; (b) Locus of vector H.
Figure 3. Measured hysteresis curves: (a) Locus of vector B; (b) Locus of vector H.
Materials 13 03138 g003
Figure 4. Measured concentric hysteresis loops in two orthogonal directions: (a) In x direction; (b) In y direction.
Figure 4. Measured concentric hysteresis loops in two orthogonal directions: (a) In x direction; (b) In y direction.
Materials 13 03138 g004
Figure 5. First-order reversal curves (FORCs) in two orthogonal directions: (a) In x direction; (b) In y direction.
Figure 5. First-order reversal curves (FORCs) in two orthogonal directions: (a) In x direction; (b) In y direction.
Materials 13 03138 g005
Figure 6. Normalized inverted scalar Everett function constructed from FORCs: (a) Fx; (b) Fy.
Figure 6. Normalized inverted scalar Everett function constructed from FORCs: (a) Fx; (b) Fy.
Materials 13 03138 g006
Figure 7. Normalized inverted scalar Everett function F and vector Everett function E: (a) Function F; (b) Function E.
Figure 7. Normalized inverted scalar Everett function F and vector Everett function E: (a) Function F; (b) Function E.
Materials 13 03138 g007
Figure 8. Flowchart of parameter extraction based on a hybrid algorithm of particle swarm optimization (PSO)–Powell.
Figure 8. Flowchart of parameter extraction based on a hybrid algorithm of particle swarm optimization (PSO)–Powell.
Materials 13 03138 g008
Figure 9. Flow diagram of simulation procedures.
Figure 9. Flow diagram of simulation procedures.
Materials 13 03138 g009
Figure 10. Variation of mean absolute percent error (MAPE) with number of iterations based on the PSO algorithm at 1.398T.
Figure 10. Variation of mean absolute percent error (MAPE) with number of iterations based on the PSO algorithm at 1.398T.
Materials 13 03138 g010
Figure 11. Comparison between measured H and simulated H based on the PSO algorithm.
Figure 11. Comparison between measured H and simulated H based on the PSO algorithm.
Materials 13 03138 g011
Figure 12. Variation of MAPE with iteration steps based on the PSO–Powell hybrid algorithm at 1.398 T: (a) Global optimization by PSO algorithm; (b) Local optimization by Powell algorithm.
Figure 12. Variation of MAPE with iteration steps based on the PSO–Powell hybrid algorithm at 1.398 T: (a) Global optimization by PSO algorithm; (b) Local optimization by Powell algorithm.
Materials 13 03138 g012
Figure 13. Comparison between experimental H and simulated H based on hybrid algorithm of PSO–Powell.
Figure 13. Comparison between experimental H and simulated H based on hybrid algorithm of PSO–Powell.
Materials 13 03138 g013
Table 1. Values of parameters and corresponding error extracted by the PSO algorithm.
Table 1. Values of parameters and corresponding error extracted by the PSO algorithm.
Bm/TwzMAPE/%Bounds of Relative Error/%
1.3980.78931.95677.394321.0209
1.0530.70042.12194.97228.9200
0.7040.69972.22314.07739.6547
0.4230.57832.21173.38018.8839
0.1780.68912.15832.89248.0391
Table 2. Values of parameters and error extracted by hybrid algorithm of PSO–Powell.
Table 2. Values of parameters and error extracted by hybrid algorithm of PSO–Powell.
Bm/TwzMAPE/%Bounds of Relative Error/%
1.3980.69512.16555.947110.7819
1.0530.69602.16664.16728.2838
0.7040.57552.17153.02897.0621
0.4230.59392.16862.92726.3889
0.1780.71822.19952.68486.5886

Share and Cite

MDPI and ACS Style

Zhao, X.; Xu, H.; Du, Z.; Li, Y.; Liu, L.; Zhao, Z. Anisotropic Vector Hysteresis Simulation of Soft Magnetic Composite Materials Based on a Hybrid Algorithm of PSO–Powell. Materials 2020, 13, 3138. https://doi.org/10.3390/ma13143138

AMA Style

Zhao X, Xu H, Du Z, Li Y, Liu L, Zhao Z. Anisotropic Vector Hysteresis Simulation of Soft Magnetic Composite Materials Based on a Hybrid Algorithm of PSO–Powell. Materials. 2020; 13(14):3138. https://doi.org/10.3390/ma13143138

Chicago/Turabian Style

Zhao, Xiaojun, Huawei Xu, Zhenbin Du, Yongjian Li, Lanrong Liu, and Zhigang Zhao. 2020. "Anisotropic Vector Hysteresis Simulation of Soft Magnetic Composite Materials Based on a Hybrid Algorithm of PSO–Powell" Materials 13, no. 14: 3138. https://doi.org/10.3390/ma13143138

APA Style

Zhao, X., Xu, H., Du, Z., Li, Y., Liu, L., & Zhao, Z. (2020). Anisotropic Vector Hysteresis Simulation of Soft Magnetic Composite Materials Based on a Hybrid Algorithm of PSO–Powell. Materials, 13(14), 3138. https://doi.org/10.3390/ma13143138

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