Next Article in Journal
Predictive Energy Management of a Building-Integrated Microgrid: A Case Study
Previous Article in Journal
Distributed Energy Resources Management System (DERMS) and Its Coordination with Transmission System: A Review and Co-Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Dual-Phase Composite Magnetic Material and Its Application in Transformers

1
State Key Laboratory of Advanced Power Transmission Technology, Beijing 102209, China
2
State Grid Smart Grid Research Institute Co., Ltd., Beijing 102209, China
3
School of Electrical Engineering, Shenyang University of Technology, Shenyang 110870, China
*
Author to whom correspondence should be addressed.
Energies 2024, 17(6), 1354; https://doi.org/10.3390/en17061354
Submission received: 3 February 2024 / Revised: 1 March 2024 / Accepted: 6 March 2024 / Published: 12 March 2024
(This article belongs to the Section F: Electrical Engineering)

Abstract

:
Dual-phase composite magnetic materials have magnetic and permanent magnetic properties. They can realize the dual-phase conversion of soft magnetic and permanent magnetic composites with a small amount of excitation energy. They have the advantages of good control and conversion characteristics and save energy, and they have a wide range of application scenarios in regard to power equipment. In this paper, the magnetization modeling of dual-phase composite magnetic materials is carried out based on micromagnetic theory, and a specific mathematical expression is given. Secondly, the preparation process of the dual-phase composite magnetic material is studied, the dual-phase composite magnetic material is prepared, and the demagnetization curve of the dual-phase composite magnetic material is measured. Finally, the application of dual-phase composite magnetic materials in power equipment is carried out. Using the soft magnetic and permanent magnetic characteristics of dual-phase composite magnetic materials, their impact on DC bias suppression in transformers is assessed. Magnetic circuit reluctance theory is used to develop the structure and electromagnetic design of a transformer. A transformer prototype with DC bias suppression ability based on dual-phase composite magnetic materials is manufactured, and simulation and experimental research are carried out. The simulation and experimental results verify the correctness of the proposed scheme. Although this scheme requires a more complex core structure, the energy-saving effect is remarkable without changing the transformer’s neutral grounding. The indicators meet the actual requirements of the project.

1. Introduction

Dual-phase composite magnetic materials have magnetic and permanent magnetic properties. They can realize the dual-phase conversion of soft magnetic and permanent magnetic composites with a small amount of excitation energy. They have the advantages of good control and conversion characteristics and save energy, and they have a wide range of application scenarios in regard to power equipment [1,2,3,4].
The DC bias phenomenon in transformers, caused by HVDC and geomagnetic storms, can seriously affect the security and stable operation of the power system [5]. Consequently, researchers have paid attention to DC bias compensation technology [6,7,8]. In [9], a DC current-blocking device for transformer neutrals is developed. In [10], an investigation was conducted to reduce the flow of geomagnetically induced currents (GICs) via power transformers. This was achieved by implementing an AC grounding channel using a capacitor between the earth or ground and the neutral point of a power transformer, but only when the GIC was identified. In [11], when designing a DC bias device for transformers, it was found that thoroughly evaluating both the resistance and the withstand voltage was important. To ensure the safe operation of the transformer when the fault zero-sequence current was high, choosing a lower resistance was recommended whenever possible. In [12], an optimal study was conducted on the DC current-blocking technique using the existing dual protection topology. The selection and design parameters for the device were determined based on applicable concepts. In [9,10,11,12], the DC bias was suppressed by calculating and installing series capacitors or resistors, but this method will raise the neutral point potential of the transformer and may cause DC bias in other transformers. Neutral point compensation technology can suppress the grounding current [13]. However, according to a large number of previous studies, this compensation method has some limitations; for example, the real-time estimation of the injected current is imprecise, leading to inadequate compensation and overcompensation [14,15]. Utilizing the auxiliary windings of hybrid transformers can achieve full compensation of the core DC bias through the establishment of a modulation compensation current [16]. However, the increased size of the compensation coil may affect the normal operation of the transformer.
The existence of DC bias in power transformers seriously affects the power quality and normal operation of the equipment [16,17,18]. Therefore, DC bias has increasingly become the focus of research by scholars [19,20,21,22,23,24,25]. The main compensation technologies are used to connect the resistance with the neutral point, connect the capacitance, reverse the electric potential compensation, install the wound degaussing coil in the transformer core, and so on. However, these methods have drawbacks that pose a great threat to the security of power systems. Currently, devising effective measures to suppress the influence of DC bias in transformer design remains a challenge.
This research analyzed and completed the calculations for the preparation of dual-phase magnetic material, the structural and electromagnetic design of a new power transformer, and the control system, comprehensively starting from the design stage and the dual-phase magnetic conversion mechanism. It provides an effective method, reference, and technical means for manufacturing the new power transformer with a DC bias compensation function. The simulation and experimental results show that the compensation effect of the transformer on DC bias is evident. The results verify the correctness and feasibility of the scheme.

2. Modeling of Dual-Phase Composite Magnetic Materials

Figure 1 displays a structural schematic diagram of the new power transformer. The control system transfers the magnetic flux to the dual-phase magnetic material when DC bias is detected. The supplemented magnetic flux is equal to the DC bias magnetic flux but in the opposite direction. When the DC bias disappears, the control system demagnetizes the dual-phase magnetic material. The transformer returns to its normal state.

3. The Establishment of the Dual-Phase Magnetic Hysteresis Mathematical Model and the Research of the Preparation of Material

3.1. Micromagnetic Modeling of Dual-Phase Magnetic Materials

Due to the lack of consideration of magnetic material properties under mechanical stress conditions in this article, magnetoelasticity is not considered. In the theory of micromagnetics, when the magnetoelastic property is neglected, the total Gibbs free energy of any magnet can be expressed as:
E t o t = E e x + E a n + E d + E ap
where E e x is the exchange interaction energy; E a n is magnetocrystalline anisotropy; E d is demagnetization energy; and E ap is saimanneng.
(1)
Exchange energy
The exchange interaction energy refers to the electrostatic interaction energy between two adjacent atoms. For substances with strong magnetism, the magnetic moments of adjacent atoms are orderly arranged due to the existence of the exchange interaction energy. The pure quantum effect caused by the Coulomb force between atoms is the source of the exchange coupling of adjacent lattice points. Exchange interaction is a short-range force, so only the nearest-neighbor interaction is considered. According to Heisenberg’s exchange theory, the exchange energy between two nearest spin atoms is:
ε e x = 2 J i j S i · S j
where the i atomic spin angular momentum is S i , the j atomic spin angular momentum is S j , and J i j is the exchange integral constant.
In addition, it is generally believed that the angle φ i j between the spin magnetic moments S i and S j of the adjacent ferromagnetic atoms is very small, so:
S i · S j = S 2 cos φ i j S 2 ( 1 1 2 φ i j 2 )
ε e x = J S 2 cos φ j j
If the exchange interaction is zero when all spins are parallel, then the exchange interaction energy when spins are not parallel is expressed as:
ε e x = J S 2 ( 1 cos φ i j ) = J S φ i j 2
φ i j = | m i m j |
where m i and m j is the S j unit magnetization vector associated with S i , m = M / | M s | . r i j is used to express the position vector from atom i to atom j . | m i m j | is expanded to the first term of Taylor series, then:
| m i m j | = | ( r i j · ) m |
Finally, the exchange density is:
ε e x = J S 2 i i j [ r i j · m ] 2
If the sum is replaced by an integral, the exchange energy density of cubic crystal is:
ε e x = 1 2 C [ ( m x ) 2 + ( m y ) 2 + ( m z ) 2 ] d τ
where C = 2 J S 2 a C , C α is the lattice constant, and the values are 1, 2, and 4 for simple cubic, body-centered cubic, and face-centered cubic, respectively. For the dense hexagonal crystal system, the expression of the above enlightenment energy is also valid. The exchange energy constant C is expressed as C = 4 2 J S 2 a C α , where α refers to the distance between the nearest atoms.
When A = C / 2 , the total exchange energy can be expressed as:
E e x = A [ ( m x ) 2 + ( m x ) 2 + ( m x ) 2 ] d V
(2)
Zeeman energy
Zeeman energy refers to the energy generated by the magnet under the action of an external field H a p , expressed as:
E α p = Ω μ 0 M · H a p d V
(3)
Magnetocrystalline anisotropy
Anisotropy energy generally includes magnetocrystalline anisotropy, surface anisotropy, shape anisotropy, and anisotropy energy caused by magnetoelasticity. In general, in micromagnetic calculations, only magnetocrystalline anisotropy is considered, and it is introduced in this paper. Magnetic-crystal anisotropy is the energy that varies with the direction of the magnetization vector. Magnetic crystal anisotropy is usually very small compared with exchange energy, but the direction of the magnetic moment is determined by anisotropy. The magnetic anisotropy properties exhibit a symmetry relationship with the crystal structure of the material. Therefore, based on the crystal symmetry, magnetic crystal anisotropy can be expressed in a phenomenological form.
For a cubic crystal system, the density of magnetocrystalline anisotropy can be expressed as follows [26]:
ε a n = K 1 ( a 1 2 a 2 2 + a 3 2 a 2 2 + a 1 2 a 3 2 ) + K 2 a 1 2 a 2 2 a 3 2 + K 3 ( a 1 2 a 2 2 + a 3 2 a 2 2 + a 1 2 a 3 2 ) 2 +
In the above formula, a 1 , a 2 , and a 3 are the M directional cosine of the saturation magnetization vector and the rectangular coordinate system; K 1 , K 2 , and K 3 represent the anisotropy constants of the magnetic crystal of the cubic crystal system, which represent the degree of the magnetic crystal anisotropy.
Order: a 2 = m y , a 1 = m x , m = M / M s , a 3 = m z ; therefore:
ε a n = K 1 ( m x 2 m y 2 + m x 2 m z 2 + m z 2 m y 2 ) + K 2 m x 2 m y 2 m z +
Generally, K 2 is smaller than K 1 , so it can be ignored. For example, for Fe , K 1 > 0 .
The easily magnetized axis is (100); when K 1 < 0 , the direction of the easily magnetized axis is (111), such as N i (nickel). The density of the magnetocrystalline anisotropy of uniaxial crystals can be expressed as:
ε a n = K 1 sin 2 θ + K 2 sin 4 θ +
where θ is the included angle between the magnetization vector M and the axis of the magnetizable axis C . K 1 and K 2 , respectively, represent the corresponding anisotropy constants of the magnetic crystal, and the above formula can also be expressed as:
ε a n = K 1 [ 1 ( e · m ) 2 ] + K 2 [ 1 ( e · m ) 2 ] 2
where e is the unit vector in the direction of the magnetizable axis. When K 1 < 0 , K 1 > 0 , and the easy-magnetized axis is on the base plane. The concept of the anisotropic field is often used to describe the anisotropic properties of magnetic crystals, which is defined as:
ε a n = μ 0 H a n · M
For cubic crystals, the anisotropic field is expressed as:
H α n ( 100 ) = 2 K 1 μ 0 M s
For cubic H α n ( 111 ) = 4 K 1 3 μ 0 M s crystal, K 1 when it is M along the (111) direction, for uniaxial crystal H α n = 2 ( K 1 + K 2 ) μ 0 M s . In most cases, K 2 is much smaller than K 1 , and H α n = 2 K 1 μ 0 M s .
The total magnetocrystalline anisotropy is:
E a n = Ω ε a n d V
(4)
Demagnetization energy
Demagnetization energy originates from the interaction between classical magnetic dipoles and contributes to the stable existence of magnetic domains. Demagnetization energy is also known as static magnetic energy. In comparison, the calculation of demagnetization energy is relatively complicated and time-consuming. Solving the demagnetization problem is crucial in micromagnetism. The demagnetizing field can be calculated by magnetic scalar potential or magnetic vector potential. This paper simply introduces the method of magnetic scalar potential. The demagnetization field can be expressed as:
H d ( r ) = 1 4 π Ω [ M ( r ) [ 3 ( r r ) ( r r ) | r r | 5 1 | r r | 3 ] d V
Demagnetization energy is expressed as:
E d = 1 2 Ω μ 0 M · H d d V
The demagnetization energy is determined by both × H d = 0 and × B = 0 , and:
B = μ 0 ( M + H d )
Because the magnetic field of the vortex is 0, which means the H d irrotational property, the magnetic scale potential is introduced:
H d = U
Then, the magnetic scale potential U satisfies the following equation:
2 U i n = M
2 U o u t = 0
According to the boundary conditions of B and H d :
n ( B 2 B 1 ) = 0
n × ( H d 2 H d 1 ) = 0
where n is the normal unit vector from medium 1 to medium 2. The boundary conditions of U are:
{ U i n = U o u t U i n n U o u t n = M n
where n is the normal unit vector from the inside of the magnet to the outside of the magnet. In addition, when r , | r U | and | r 2 U | are limited. According to the uniqueness theorem, the U solution satisfying the above equations and boundary conditions is unique, and the general solution of H d and U can be obtained as follows:
H d ( r ) = 1 4 π Ω ( r r ) · M ( r ) | r r | 3 d V + 1 4 π Ω ( r r ) n · M ( r ) | r r | 3 d S
U ( r ) = 1 4 π Ω · M ( r ) | r r | 3 d V + 1 4 π Ω n · M ( r ) | r r | 3 d S
Due to the surface magnetic charge σ m = n M ( r ) and magnetic charge volume ρ m = n M ( r ) , the demagnetization field and magnetic scalar potential are expressed as follows:
H d = 1 4 π Ω ( r r ) ρ m ( r ) | r r | 3 d V + 1 4 π Ω ( r r ) σ m ( r ) | r r | 3 d S
U ( r ) = 1 4 π Ω ρ m ( r ) | r r | 3 d V + 1 4 π Ω σ m ( r ) | r r | 3 d S
From (30)–(31), the demagnetization field is the combination of the demagnetization field resulting from the surface magnetic charge and the body magnetic charge. Demagnetization energy can also be expressed as:
E d = 1 2 μ 0 Ω ρ m ( r ) U ( r ) d V + 1 2 μ 0 Ω σ m ( r ) U ( r ) d S
Equation (32) is the expression of demagnetization field and demagnetization energy, but it is difficult to complete such a large amount of calculation on an ordinarily configured computer. In view of this, researchers have adopted other acceleration methods to solve this problem, mainly the Hierarchical Calculation method and the Fast Fourier Transform method, which greatly improve the calculational speed.

3.2. Static and Dynamic Methods

(1)
Static calculation method
When the Gibbs free energy of the magnet is taken as a minimum, the equilibrium distribution of the magnetization vector can be obtained. The static micromagnetics calculation is based on the equilibrium condition of taking the minimum value from Etot or the condition that Etot’s variation is 0, and the Brown equation can be obtained:
{ m × [ 2 A m n + K s ( n · m ) n ] = 0 m × [ 2 A 2 m + μ 0 M s ( H o p + H d ) ω a n m ] = 0
When surface anisotropy is ignored, the boundary condition can be expressed as:
m n = n ( m ) = 0
where the concept of effective field H e f f is cited. Field effectiveness H e f f is described as the variation in free energy, which provides the actual moment acting on the magnetization vector.
H e f f = 1 μ 0 M s E t o t m = 1 μ 0 M s E e x m 1 μ 0 M s E a n m 1 μ 0 M s E a p m 1 μ 0 M s E d m = H e x + H a p + H d + H a n
Brown’s equation is written as:
μ 0 M s m × H e f f = 0
When the energy is very small, the directions of H e f f and m in the corresponding magnet are the same everywhere. The static micromagnetics calculation directly minimizes the energy according to the above principle, thus obtaining the distribution of the magnetization vector in the equilibrium state.
(2)
Dynamic calculation method
The expression of the total free energy of the magnet is written as (37):
E t o t = E e x + E d + E a n + E H
When the system is under a changing external field, the dynamic equation of the magnetic moment must be considered, which follows the Landau–Lifshitz–Gilbert (LLG) dynamic equation:
d M d t = ω 1 + a 2 M × H e f f a ω ( 1 + a 2 ) M s M × ( M × H e f f )
where M is the magnetization vector, α is the damping coefficient, ω is the rotational magnetic ratio, and the effective field H e f f is defined as the variation in free energy, H e f f = 1 μ 0 E t o t M , which provides the actual torque acting on the magnetization vector.
A small change in each magnetic moment in the system will cause a change in the energy state of the whole system. In other words, the dynamic equation of each magnetic moment is closely related, so the system will get the minimum energy, which is also the criterion for whether the calculation can be terminated. The boundary condition of the dynamic equation is the same as the static boundary condition, that is, when surface anisotropy is ignored, and the boundary condition is:
m n = 0
The dynamic micromagnetic calculation solves the above Gilbert or Landau–Lifshitz dynamic equation and then obtains the image of the magnetization vector moving with time.

3.3. Hysteresis Characteristics of Nanocomposite Magnetic Materials

This paper is based on the theory of micromagnetism, employing a three-dimensional dynamic model, and the dynamic calculation method of micromagnetism is used to calculate Nd2Fe14B/ α -Fe/Nd2Fe14B on the OOMMF 2.0 (beta) software platform. The magnetization process of Nd2Fe14B exchange-coupled magnetic three-layer film is simulated, and the hysteresis loops of nano-phase composite Nd2Fe14B three-layer film under different soft magnetic layer thickness and the relationship curve between soft magnetic layer thickness and coercivity are calculated. The three-layer film model is shown in Figure 2. n1 represents the thickness of the rigid magnetic stratum, n2 represents the thickness of the soft magnetic layer, and l 1 and l 2 , respectively, represent the width and length of the three-layer film model.
In the model calculated in this paper l 1 = l 2 = 300   nm , along the Z-axis, the external magnetic field is directed. The length and width of the split unit are 10 nm, and the height is 1 nm. The calculation parameters of the nano-phase composite Nd2Fe14B three-layer film are shown in Table 1.
In this paper, when the hard magnetic layer is 5 nm thick, the calculation of the soft magnetic layer is in the range of 1 nm to 20 nm, and the hysteresis loops under different soft magnetic layer thicknesses are obtained, as shown in Figure 3. The relationship curve between the coercivity and the soft magnetic layer thickness is obtained, as shown in Figure 4.
The relationship between the microscopic parameters and the macroscopic magnetic properties of the nanocomposite magnetic materials was established by using the micromagnetic theory and the OOMMF software platform to calculate the hysteresis properties of the nanocomposite magnetic materials. In this project, the microscopic parameters of the nanocomposite magnetic materials closest to the magnetic properties can be obtained according to the requirements of the material properties, which can be used to determine the preparation of the materials.

3.4. Research of the Influence of Element Composition and Heat Treatment Craft on the Dual-Phase Magnetic Material

Based on the above analysis, the remanence, coercivity, and hysteresis loop are defined. In order to manufacture the dual-phase magnetic material, the influence of elements and heat treatment craft to the dual-phase magnetic material are discussed below.
(1)
Influence of Element Composition on the Dual-phase Magnetic Material
The chemical formula for the dual-phase magnetic material is Fe-12Co-24Cr. The influence of the contents of the element on the dual-phase magnetic material is significant. A higher content of Co alloy leads to better magnetic performance. However, its workability is poor. If the content of the Co alloy is low, its workability is good, but the craft is complex and difficult to control. The adjustment range of the magnetic parameter will be narrow. The content fluctuation of the Cr will affect the coercivity of the material. Studies have shown that if the fluctuation of Cr is 3%, the fluctuation of coercivity will be 40%.
(2)
Influence of Heat Treatment Craft on the Dual-phase Magnetic Material
The heat treatment craft is closely related to the magnetic performance of the material. Temperature control has a significant impact on both the strong magnetic phase and the weak magnetic phase. Based on a large number of experiments, the five-rank tempering craft was adopted. The implementation process is listed as follows:
620 °C × 1 h + 610 °C × 1 h + 590 °C × 2 h + 570 °C × 3 h + 560 °C ×3 h.
(3)
Preparation of the Dual-phase Magnetic Material
Based on the above analysis, the dual-phase magnetic material can be prepared. Figure 5 shows the sample of the dual-phase magnetic material.

4. Structure Design of the New Power Transformer

According to the latest research results on dual-phase magnetic material, a more reasonable solution is to utilize the materials on both sides of the new power transformer adjacent to the column in order to achieve DC bias compensation. The equivalent magnetic circuit of the new power transformer was analyzed based on the reluctance theory under the condition of the DC magnetic potential, and the AC magnetic potential was considered, respectively. Figure 6 shows the structure of the new transformer.
The equivalent magnetic circuit of the new power transformer is shown in Figure 7. (IN)AC represents the AC potential, (IN)DC represents the DC potential, and (IN)Magnet represents the compensation magnetic potential.
Figure 8 shows the equivalent magnetic circuit of the new power transformer under the condition of considering the DC flux only. Equation (40) can be obtained based on Kirchhoff’s voltage and current laws.
{ ( I N ) D C = ϕ 1 [ R c + ( R c + R g 1 ) / / R g 2 ] ( I N ) M a g n e t = ϕ 2 [ R g 2 + ( R c + R g 1 ) / / R c ] ϕ 1 + ϕ 2 + ϕ 3 = 0
where R c is the reluctance of the silicon steel, and ϕ 1 and ϕ 3 are the DC flux. The DC flux is not expected to flow to the adjacent column, so ϕ 2 is set to zero. R g 1 and R g 2 can be obtained. That is to say, the gap can be obtained.
Figure 9 shows the equivalent magnetic circuit of the new power transformer under the condition of considering the AC flux only. Equation (41) can be obtained based on Kirchhoff’s voltage and current laws.
{ ( I N ) A C = ϕ 4 [ R c + ( R c + R g 1 ) / / R g 2 ] ϕ 5 R g 2 = ϕ 6 ( R c + R g 1 ) ϕ 4 + ϕ 5 + ϕ 6 = 0
where R c is the reluctance of the silicon steel, and ϕ 5 and ϕ 6 are the AC flux. The AC flux is not expected to flow the dual-phase magnetic material, so ϕ 4 is set to zero. R g 1 and R g 2 can be obtained. That is to say, the gap can be obtained. According to (40) and (41), the gap length of the 220 V/110 V/2 kVA new transformer is shown in Table 2. Figure 10 shows the prepared model.

5. Simulation and Experimental Study on the New Transformer with the DC Bias Compensation Function

In Section 4, the design of the transformer was carried out. This section will carry out simulation research on the transformer.

5.1. Simulation of the New Transformer with the DC Bias Compensation Function

Figure 11 shows the magnetic field distribution of the new transformer without DC bias.
Figure 12 shows the magnetic field distribution of the new transformer when the DC bias is 0.3 A. When the DC bias is 0.3 A, the magnetic density of the positive half cycle and the negative half cycle is not equal. The magnetic density of the positive half-cycle is about 1.3 T, and the magnetic density of the negative half-cycle is about 0.85 T.
Figure 13 shows the magnetic field distribution of the new transformer after compensation when the DC bias is 0.3 A. After compensation, the magnetic densities of the positive half-cycle and the negative half-cycle are close to each other. The magnetic density of the positive half-cycle is about 1.21 T, and the magnetic density of the negative half-cycle is about 0.86 T.
Figure 14 shows the magnetic field distribution of the new transformer when the DC bias is 0.6 A. When the DC bias is 0.6 A, the magnetic densities of the positive half-cycle and the negative half-cycle are not equal. The magnetic density of the positive half-cycle is about 1.38 T, and the magnetic density of the negative half-cycle is about 0.79 T.
Figure 15 shows the magnetic field distribution of the new transformer after compensation when the DC bias is 0.6 A. After compensation, the magnetic density of the positive half-cycle is about 1.32 T, and about 0.83 T in the negative half-cycle.

5.2. Experiment Simulation of the New Transformer with the DC Bias Compensation Function

The experimental site and experimental circuit diagram are shown in Figure 16.
Figure 17 is the excitation current waveform before and after the new transformer compensation when the bias current is 0.3 A. Figure 18 shows the excitation current waveform before and after the compensation of the new transformer when the bias current is 0.6 A.
The analysis indicates that the new transformer demonstrates effective compensation for bias current. When the DC bias phenomenon occurs, the positive peak of the excitation current increases, the negative peak decreases, and the working point of the transformer exciting current increases. When the bias current is compensated, the positive peak of the excitation current decreases and the negative peak increases. This compensation mechanism significantly contributes to addressing transformer faults caused by DC bias, aligning with the system design requirements.

6. Conclusions

Based on the dual-phase hysteresis mathematical model, preparation, structure, and electromagnetic design of dual-phase magnetic materials, and simulation and experiments, a new type of transformer was designed and analyzed in detail. The conclusions are as follows:
(1)
Based on the theory of micromagnetism, a hysteresis model of nano dual-phase composite magnetic materials was established using the OOMMF software platform and the micromagnetic simulation dynamic calculation method. The hysteresis loops of nano dual-phase composite magnetic materials under different soft magnetic layer thicknesses were calculated, and the relationship between soft magnetic layer thickness and coercivity was obtained.
(2)
The structure of the transformer was designed and used for simulation verification. The simulation results show that when DC bias occurs, almost all of the compensating magnetic flux of the compensating iron core passes through the main magnetic circuit of the transformer, canceling out the DC magnetic flux in the main magnetic circuit and achieving the goal of eliminating DC bias, and the working magnetic flux of the transformer does not demagnetize the compensating iron core. When there is no DC bias, the compensating iron core has no residual magnetism, and the magnetic circuit of the compensating iron core does not affect the normal operation of the transformer. This proves the feasibility of the voltage regulator structure.
(3)
A 2 kVA power transformer prototype was developed, and experiments were conducted on the prototype. The experimental results show that under different DC bias conditions, after compensating for DC bias, the excitation current waveform of the transformer tends to be symmetrical. The positive peak value decreases and the negative peak value increases. The positive half-cycle magnetic flux density decreases, the iron core returns to an unsaturated state, and the negative half-cycle magnetic flux density increases, proving the suppression effect of the power transformer on DC bias.
The simulation and experimental results prove that the structure of the power transformer is reasonable, the DC bias suppression method is feasible, and it has a good suppression effect on DC bias. It realizes the automatic suppression of DC bias by the transformer itself, contributing to further research and solution of DC bias problems.

Author Contributions

Conceptualization, Y.L. and F.Y.; methodology, Y.H.; validation, J.G.; formal analysis, D.C.; investigation, Y.L.; resources, F.Y. and Y.H.; data curation, D.C. and H.B.; writing—original draft preparation, F.Y.; writing—review and editing, Y.L.; visualization, H.B.; funding acquisition, J.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Science and Technology Project of State Grid Corporation of China (5500-202258312A-2-0-QZ).

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

Authors Yang Liu, Fuyao Yang, Yu Han and Jie Gao were employed by the State Grid Smart Grid Research Institute Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Cheng, S.; Cai, Z. Design and optimization of energy-saving wind power grid-connected contactor based on nano two-phase composite magnetic materials. In Proceedings of the 2017 4th International Conference on Electric Power Equipment-Switching Technology (ICEPE-ST), Xi’an, China, 22–25 October 2017; pp. 209–212. [Google Scholar]
  2. Kharitonskii, P.V.; Kosterov, A.A.; Gurylev, A.K.; Gareev, K.G.; Kirillova, S.A.; Zolotov, N.A.; Anikieva, Y.A. Magnetic States of Two-Phase Synthesized FemOn–Fe3–xTixO4 Particles: Experimental and Theoretical Analysis. Phys. Solid State 2020, 62, 1691–1694. [Google Scholar] [CrossRef]
  3. Kharitonskii, P.V.; Gareev, K.G.; Ralin, A.Y.; Sergienko, E.S. Superparamagnetism of Fe3O4–Fe3–xTixO4 Composites: Micromagnetic Modeling. Phys. Metals Metallogr. 2023, 124, 46–52. [Google Scholar] [CrossRef]
  4. Kramer, N.; Provino, A.; Hanley, C.; Kubiak, C.; Pathak, A.K. Magnetic properties of LaFe2−xCoxSi: A new two-phase naturally formed composite alloy. AIP Adv. 2023, 13, 025320. [Google Scholar] [CrossRef]
  5. Feng, H.; Wang, K.; Li, J.; Zhang, S.; Xie, Q.; Hu, X.; Li, B.; Wang, S.; Wang, S. Analysis of Exciting Current for Single-phase Four-Limb Ultra-High Voltage Transformer under DC Bias. In Proceedings of the 2018 Condition Monitoring and Diagnosis (CMD), Perth, WA, Australia, 23–26 September 2018; pp. 1–6. [Google Scholar]
  6. Cai, Y.; Yu, Q.; Guo, M.; Sun, G.; Sun, Q. Effect of heat treatments on the microstructure and magnetic properties of semi-hard magnetic 2J85 alloy. J. Magn. Mater. Devices 2013, 1, 47–51. [Google Scholar]
  7. Chapman, S.J. Electric Machinery Fundamentals; Publishing House of Electronics Industry: Beijing, China, 2012; pp. 1–30. [Google Scholar]
  8. Zhu, L.; Wei, C.; Yu, Y. Analysis of DC bias exciting current of the single-phase transformer and its effect on protection. Power Syst. Prot. Control 2010, 38, 158–162. [Google Scholar]
  9. Bolduc, L.; Granger, M.; Pare, G.; Saintonge, J.; Brophy, L. Development of a DC current-blocking device for transformer neutrals. IEEE Trans. Power Deliv. 2005, 20, 163–168. [Google Scholar] [CrossRef]
  10. Rajapakse, A.; Perera, N.; Faxvog, F.R.; Jensen, W.; Nordling, G.; Fuchs, G.; Jackson, D.B.; Volkmann, T.L.; Ruehl, N.; Groh, B. Power grid stability protection against GIC using a capacitive grounding circuit. In PES T&D; IEEE: Piscataway, NJ, USA, 2012; pp. 1–6. [Google Scholar]
  11. Hongtao, R. Analysis on Withstand Voltage of DC bias Device at 500kV Transformer Neutral Point. In Proceedings of the 2020 3rd International Conference on Power and Energy Applications (ICPEA), Busan, Republic of Korea, 9–11 October 2020; pp. 58–62. [Google Scholar]
  12. Li, C.; Zhang, Y.; Lv, Y.; Sun, S. Research on the design and experiment of power transformer DC bias blocking device. In Proceedings of the 2012 2nd International Conference on Consumer Electronics, Communications and Networks (CECNet), Yichang, China, 21–23 April 2012; pp. 2811–2814. [Google Scholar]
  13. Lasabi, O.; Swanson, A.; Jarvis, L. Heuristic Control of Neutral DC Compensation Method to Moderate DC bias in Power Transformer. In Proceedings of the 2022 IEEE PES/IAS PowerAfrica, Kigali, Rwanda, 22–26 August 2022; pp. 1–5. [Google Scholar]
  14. Pan, Z.; Wang, X.; Tan, B.; Zhu, L.; Liu, Y.; Liu, Y.; Wen, X. Potential compensation method for restraining the DC bias of transformers during HVDC monopolar operation. IEEE Trans. Power Deliv. 2015, 31, 103–111. [Google Scholar] [CrossRef]
  15. Bachinger, F.; Hackl, A.; Hamberger, P.; Leikermoser, A.; Leber, G.; Passath, H.; Stoessl, M. Direct current in transformers: Effects and compensation. E&I Elektrotechnik Informationstechnik 2013, 1–5. [Google Scholar] [CrossRef]
  16. Fei, X.; Zhiwei, C.; Linzi, D.; Linjia, H.; Chen, J.; Gang, W. DC Bias Suppression Scheme Based on Hybrid Transformer. In Proceedings of the 2022 25th International Conference on Electrical Machines and Systems (ICEMS), Chiang Mai, Thailand, 29 November–2 December 2022; pp. 1–5. [Google Scholar]
  17. Zhang, X.; Wu, W. Simulation analysis of converter transformer excitation current and reactive power loss in DC bias based on field circuit coupling. In Proceedings of the 2018 IEEE 4th Information Technology and Mechatronics Engineering Conference (ITOEC), Chongqing, China, 14–16 December 2018; pp. 861–864. [Google Scholar]
  18. Ni, Y.; Zeng, X.; Yu, K.; Sun, S.; Zhuo, C.; Peng, P.; Su, Z. Modeling calculation and influence analysis of transformer DC bias current caused by metro stray current. In Proceedings of the 2021 IEEE 4th International Electrical and Energy Conference (CIEEC), Wuhan, China, 28–30 May 2021; pp. 1–6. [Google Scholar]
  19. Zhao, J.; Fu, H.; Li, J.; Liu, Y.; Li, Z.; Pan, Z. The DC Bias Risk Assessment of Transformer Based on Nonlinear Magnetic Circuit Equations. In Proceedings of the 2019 IEEE 3rd Conference on Energy Internet and Energy System Integration (EI2), Changsha, China, 8–10 November 2019; pp. 2177–2181. [Google Scholar]
  20. Yan, X.; Dong, X.; Han, G.; Yin, X.; Yu, X.; Jiao, Z.; Wang, X. Study of Transformer Loss and Temperature Rise Under DC Bias Magnetism Based on Finite Element Method. In Proceedings of the 2023 IEEE/IAS Industrial and Commercial Power System Asia (I&CPS Asia), Chongqing, China, 14–16 July 2023; pp. 429–434. [Google Scholar]
  21. Subramanya, K.; Chelliah, T.R. Thermo Electrical Analysis and its Correlation in a DC Biased Power Transformer Subjected to Over and Under Excitations. In Proceedings of the 2023 13th International Conference on Power, Energy and Electrical Engineering (CPEEE), Tokyo, Japan, 25–27 February 2023; pp. 106–111. [Google Scholar]
  22. Wang, G.; Zhang, G.; Liu, H.; Yu, D.; Wang, Q.; Chen, L. Analysis of the Influence of Transformer Harmonics Caused by DC Bias on Reactive Power Compensation Capacitor Banks. In Proceedings of the 2022 IEEE 5th International Electrical and Energy Conference (CIEEC), Nangjing, China, 27–29 May 2022; pp. 4576–4581. [Google Scholar]
  23. Zheng, Y.; Pan, S.; Xia, Y.; Ouyang, X.; Shang, H.; Zhou, Q. Analysis on the Influence of the DC bias on spatial magnetic flux leakage in the transformer. In Proceedings of the 2022 IEEE International Conference on High Voltage Engineering and Applications (ICHVE), Chongqing, China, 25–29 September 2022; pp. 1–4. [Google Scholar]
  24. Zhou, Z.; Hui, X.; Zhu, X.; Zhu, Y.; Xia, R. Research Summary on Transformer DC Magnetic Bias and Suppression Technology. In Proceedings of the 2023 8th Asia Conference on Power and Electrical Engineering (ACPEE), Tianjin, China, 14–16 April 2023; pp. 2715–2719. [Google Scholar]
  25. Zheng, L.; Su, K.; Liu, X.; Zhao, J.; Guo, H.; Xun, H. Analysis of DC Magnetic Bias Influences on Transformers. In Proceedings of the 2020 IEEE 4th Conference on Energy Internet and Energy System Integration (EI2), Wuhan, China, 30 October–1 November 2020; pp. 3329–3332. [Google Scholar]
  26. Chikazumi, S. Physics of Ferromagnetism; Oxford University Press: New York, NY, USA, 1997. [Google Scholar]
Figure 1. Structure schematic diagram of the new power transformer.
Figure 1. Structure schematic diagram of the new power transformer.
Energies 17 01354 g001
Figure 2. Calculation model of nanometer dual-phase three-layer film.
Figure 2. Calculation model of nanometer dual-phase three-layer film.
Energies 17 01354 g002
Figure 3. Nano-phase composite Nd2Fe14B hysteresis loop family.
Figure 3. Nano-phase composite Nd2Fe14B hysteresis loop family.
Energies 17 01354 g003
Figure 4. Relationship curve between soft magnetic layer thickness and coercivity.
Figure 4. Relationship curve between soft magnetic layer thickness and coercivity.
Energies 17 01354 g004
Figure 5. Demonstration of the dual-phase magnetic material.
Figure 5. Demonstration of the dual-phase magnetic material.
Energies 17 01354 g005
Figure 6. Structure of the new transformer.
Figure 6. Structure of the new transformer.
Energies 17 01354 g006
Figure 7. Equivalent magnetic circuit model.
Figure 7. Equivalent magnetic circuit model.
Energies 17 01354 g007
Figure 8. Equivalent magnetic circuit model considering DC flux only.
Figure 8. Equivalent magnetic circuit model considering DC flux only.
Energies 17 01354 g008
Figure 9. Equivalent magnetic circuit model considering AC flux only.
Figure 9. Equivalent magnetic circuit model considering AC flux only.
Energies 17 01354 g009
Figure 10. Prototype of the new transformer.
Figure 10. Prototype of the new transformer.
Energies 17 01354 g010
Figure 11. Magnetic field distribution of the new transformer without DC bias: (a) positive half-cycle and (b) negative half-cycle.
Figure 11. Magnetic field distribution of the new transformer without DC bias: (a) positive half-cycle and (b) negative half-cycle.
Energies 17 01354 g011
Figure 12. Magnetic field distribution of the new transformer when the DC bias is 0.3 A: (a) positive half-cycle and (b) negative half-cycle.
Figure 12. Magnetic field distribution of the new transformer when the DC bias is 0.3 A: (a) positive half-cycle and (b) negative half-cycle.
Energies 17 01354 g012
Figure 13. Magnetic field distribution of the new transformer after compensation when the DC bias is 0.3 A: (a) positive half-cycle and (b) negative half-cycle.
Figure 13. Magnetic field distribution of the new transformer after compensation when the DC bias is 0.3 A: (a) positive half-cycle and (b) negative half-cycle.
Energies 17 01354 g013aEnergies 17 01354 g013b
Figure 14. Magnetic field distribution of the new transformer when the DC bias is 0.6 A: (a) positive half-cycle and (b) negative half-cycle.
Figure 14. Magnetic field distribution of the new transformer when the DC bias is 0.6 A: (a) positive half-cycle and (b) negative half-cycle.
Energies 17 01354 g014
Figure 15. Magnetic field distribution of the new transformer after compensation when the DC bias is 0.6 A: (a) positive half-cycle and (b) negative half-cycle.
Figure 15. Magnetic field distribution of the new transformer after compensation when the DC bias is 0.6 A: (a) positive half-cycle and (b) negative half-cycle.
Energies 17 01354 g015
Figure 16. The experimental site and experimental circuit diagram: (a) experimental site and (b) experimental circuit diagram.
Figure 16. The experimental site and experimental circuit diagram: (a) experimental site and (b) experimental circuit diagram.
Energies 17 01354 g016
Figure 17. The excitation current waveform before and after the new transformer compensation when the bias current is 0.3 A.
Figure 17. The excitation current waveform before and after the new transformer compensation when the bias current is 0.3 A.
Energies 17 01354 g017
Figure 18. The excitation current waveform before and after the new transformer compensation when the bias current is 0.6 A.
Figure 18. The excitation current waveform before and after the new transformer compensation when the bias current is 0.6 A.
Energies 17 01354 g018
Table 1. Parameters of nano-phase composite Nd2Fe14B three-layer film.
Table 1. Parameters of nano-phase composite Nd2Fe14B three-layer film.
ItemSaturated
Magnetization
(kA/m)
Anisotropic Constant
(MJ/m3)
Exchange Coupling
Constant (J/m)
Hard magnetic layer12884.37.7 × 10−12
Soft magnetosphere17200.0462.5 × 10−11
Table 2. Types sizes of the new transformer.
Table 2. Types sizes of the new transformer.
Transformer Typegap1 (mm)gap2 (mm)Length of the Material Magnetic Circuit (mm)
220 V/110 V/2 kVA0.3–0.70.1–0.4128
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, Y.; Yang, F.; Han, Y.; Gao, J.; Chen, D.; Bai, H. Modeling of Dual-Phase Composite Magnetic Material and Its Application in Transformers. Energies 2024, 17, 1354. https://doi.org/10.3390/en17061354

AMA Style

Liu Y, Yang F, Han Y, Gao J, Chen D, Bai H. Modeling of Dual-Phase Composite Magnetic Material and Its Application in Transformers. Energies. 2024; 17(6):1354. https://doi.org/10.3390/en17061354

Chicago/Turabian Style

Liu, Yang, Fuyao Yang, Yu Han, Jie Gao, Dezhi Chen, and Haonan Bai. 2024. "Modeling of Dual-Phase Composite Magnetic Material and Its Application in Transformers" Energies 17, no. 6: 1354. https://doi.org/10.3390/en17061354

APA Style

Liu, Y., Yang, F., Han, Y., Gao, J., Chen, D., & Bai, H. (2024). Modeling of Dual-Phase Composite Magnetic Material and Its Application in Transformers. Energies, 17(6), 1354. https://doi.org/10.3390/en17061354

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