# A New 3D Method for Reactor Core Vibration Based on Silicon Steel Lamination Rules and Application in UHV Shunt Reactors

Authors: Fangcheng Lü; Jiayi Guo; Leilei Niu; Jianghai Geng; Yirui Pan
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

A new three-dimensional (3D) analysis method is proposed as the existing two-dimensional (2D) method has low accuracy in analysing the vibration characteristics of oil-immersed shunt reactors, such as ultrahigh-voltage (UHV) shunt reactors. First of all, a set of 3D laminated coordinate systems was defined based on silicon steel lamination rules, in which the anisotropy of magnetic properties for laminated silicon steel in the rolling direction (RD), the transverse direction (TD), and the lamination direction (LD) were considered. Then, the mapping between laminated coordinate systems and space rectangular coordinate system was established to unify the parameters in different laminated coordinate systems. With the mapping, the anisotropy of the magnetic properties in the laminated coordinate systems was transformed into a rectangular coordinate system. Next, two sets of comparative studies between the new 3D method and the traditional 2D method were carried out, which show that the 3D method has high precision and a wide application range. Finally, the relationship between air gap number and core vibration of UHV shunt reactors was studied by the new 3D method. The results show that, as the number of air gaps increases, the magnetic flux density and the total force area of Maxwell force are increased, resulting in the intensification of core vibration. The conclusions of this paper are helpful for the design of large oil-immersed reactors.


---

## Body

## 1. Introduction

A shunt reactor, which compensates for long-distance line capacitance and suppressing the increase of power frequency, is one of the important items of equipment in AC transmission projects. Excessive vibration is one of the most common problems in the operation of shunt reactors due to the air gaps in the core limb that provide large inductance [1].The vibration of the core (iron yoke and core limb) is one of the most important causes of shunt reactor vibration [1–9]. The elastic deformation of the core in the magnetic field due to Maxwell force and magnetostriction is considered to be the main reason for shunt reactor core vibration [2]. 2D methods have been widely used to study this [3, 8]. In this type of method, the core limb is considered to have the same thickness in the z-direction. The difference in magnetic properties (magnetisation curves, magnetostrictive characteristic parameters, and so on) of the core material is called anisotropy. The 3D anisotropy (in the direction of RD, TD, and LD) has been simplified into 2D (RD and TD) in the 2D method. The 2D method has high accuracy in the study of small dry reactors with cuboid core limbs [3, 4]. However, the core limb of many large oil-immersed shunt reactors such as UHV shunt reactors is a cylinder with laminated silicon steel, which is very different from the core limb of a small dry reactor. Thus, a new and wider range of method is needed.The air gap in the core limb is one of the main reasons why the shunt reactor vibration is greater than in other similar items of equipment such as transformers. Adjusting the air gaps can be a useful method to reduce the vibration of shunt reactors. The method of adding air gaps in the yokes and adjusting the hardness of the air gap material was proposed by Gao et al. [4]. The method of filling the air gap by a novel method inserting negative giant magnetostrictive material (NGMM) was proposed by Yan et al. [5]. The aforementioned research has paid more attention to the relationship between air gap material and vibration. The influence of air gap structure on core vibration remains to be studied.Here, a set of 3D laminated coordinate systems was constructed based on the lamination structure of the shunt reactor core, and a new 3D analysis method for the core vibration of shunt reactor column structure was proposed. Then, the 3D method was compared with the traditional 2D method by two sets of comparative studies; this method has a wider range of applications than the 2D method. Finally, the 3D method was applied to study the relationship between air gap number and core vibration of a UHV shunt reactor, which is typically a large oil-immersed reactor.


## 2. Methods

### 2.1. Laminated Coordinate Systems

The vibration of the core (iron yoke and core limb) is one of the most important causes of the shunt reactor vibration. The core of the shunt reactor is composed of laminated silicon steel sheets, whose magnetic properties are related to the direction of lamination [6]: there are different lamination directions in different parts of the core. According to the rules governing the lamination of silicon steel sheets, three mutually perpendicular axes including RD, TD, and LD can be set; on this basis, a set of 3D coordinate systems is created. The anisotropy of magnetic properties for laminated silicon steel can be expressed more clearly thereby.There are two types of shunt reactor, dry type with small capacity in Figure1(a) and oil-immersed type with large capacity in Figure 1(b) [3, 7]. Core limb and iron yoke are the two main parts of the reactor core. The core limb includes core discs and air gaps (marble or ceramics). The core discs for small dry reactors are usually cuboid, and core discs for large oil-immersed shunt reactors are generally cylindrical. The mapping between the laminated coordinate systems and the spatial rectangular coordinate system is as shown in Table 1.Figure 1
Simplified diagrams of two types of reactors: (a) small dry reactor and (b) large oil-immersed shunt reactor.Table 1
The mapping between the two coordinate systems.
Part
RD
TD
LD
Large oil-immersed shunt reactors
Core disc
z-axis
Vertical to radial andz-axis
Radial
Top and bottom yokes
y-axis
z-axis
x-axis
Side yokes
z-axis
y-axis
x-axis
Small dry reactors
Core discs
y-axis
x-axis
z-axis
Yokes
x-axis
y-axis
z-axisThe lamination rules of the cores at different positions are different, so the number of lamination coordinate systems is large, which entails onerous computations. It is necessary to define a global coordinate system (a space rectangular coordinate system here), which is used to unify parameters of different coordinate systems. Any vector has two different representations belonging to the two coordinate systems. The mapping between these two representations is beneficial to the calculation of the vibration parameters of the reactor.


#### 2.1.1. Large Oil-Immersed Shunt Reactors

As shown in Table1, for any vector M (in 3D coordinates) in the iron yokes of large oil-immersed shunt reactors, the mapping between the two coordinate systems satisfying equation (1) (top and bottom yokes) and equation (2) (side yokes) is as follows:(1)Mld=Mx,Mtd=Mz,Mrd=My,(2)Mld=Mx,Mtd=My,Mrd=Mz.The core disc (Figure2) of a large oil-immersed shunt reactor is generally formed by splicing a plurality of fan-shaped laminated silicon steels. The fan-shaped laminated silicon steel is formed by stacking a plurality of rectangular silicon steel sheets with different lengths. The laminated coordinate system is related to the location of each sheet, and different laminated silicon steels belong to different 3D laminated coordinate systems. Taking the fanned laminated silicon steel marked with the border in Figure 2 as an example, the dotted lines indicate the LD (coinciding with the centreline of the sector), and the TD of the laminated coordinate system in which the laminated silicon steel is located. The solid lines indicate the x and y-axes in a 3D space rectangular coordinate system. The RD is recombined with the z-axis perpendicular to the core disc plane. The angle between the x-axis and the LD is α. The mapping of any vector M (in 3D coordinates) in the core discs between the two coordinate systems is as follows:(3)Mld=Mxcosα+Mysinα,Mtd=−Mxsinα+Mycosα,Mrd=Mz.Figure 2
A plan view of a core disc: the green arrow is the coordinate axis of the space rectangular coordinate system, and the red arrow represents the coordinate axis of the laminated coordinate system.


#### 2.1.2. Small Dry Reactors

As shown in Table1, the mapping for dry reactors can be performed by equations (4) (yokes) and (5) (core discs).(4)Mld=Mz,Mtd=My,Mrd=Mx,(5)Mld=Mz,Mtd=Mx,Mrd=My,where Mld, Mtd, and Mrd are the 3D components of the vector M in the lamination coordinate system, and Mx, My, and Mz are the 3D components of the vector M in the 3D space rectangular coordinate system in equations (1) to (5).


### 2.2. Vibration Analysis Method

The vibration of the core is caused by the Maxwell force and magnetostriction in the magnetic field, which are affected by the spatial distribution of the magnetic field.In the spatial rectangular coordinate system, the 3D electromagnetic equation of the reactor core under normal operating conditions is given by(6)∂∂yυz∂Ay∂x−∂Ax∂y−∂∂zυy∂Ax∂z−∂Az∂x=Jx,∂∂zυx∂Az∂y−∂Ay∂z−∂∂xυz∂Ay∂x−∂Ax∂y=Jy,∂∂xυy∂Ax∂z−∂Az∂x−∂∂yυx∂Az∂y−∂Ay∂z=Jz,where υx, υy, and υz are the components of the magnetoresistance ratio υ; Ax, Ay, and Az are the components of the vector magnetic potential A; and Jx, Jy, and Jz are the components of the winding current density J.The magnetic flux densityB in the spatial rectangular coordinate system can be calculated by(7)B=∇×A.The magnetic flux densityB in the lamination coordinate system was obtained through coordinate transformation (equations (1) to (5)) of B in the spatial rectangular coordinate system. The magnetic field strength H was calculated by the 3D electromagnetic field constitutive relationship as follows:(8)H=HrdHtdHld=υrd000υtd000υldBrdBtdBld.The 3D magnetisation curves (H-B curves, where H is the magnetic field strength in the lamination coordinate system and H is the magnetic flux density in the laminated coordinate system) of a type of laminated silicon steel were obtained by tests [10] (Figure 3).Figure 3
The 3D magnetisation curves (H-B curves) of the laminated silicon steel.The Maxwell forceFmax acting on the interfaces between core discs and air gaps in the magnetic field can be obtained from equation (9) [11–13]:(9)Fmax=∫sTmaxdS,where Tmax is the Maxwell stress tensor (equation (11)).The magnetostriction phenomenon is caused by the displacement of the magnetic domain inside the core (silicon steel) in the magnetic field, and its magnitude is closely related to the magnetic field. Magnetostrictive displacement and strain are obtained from the magnetic field strength by the magnetostrictive curves of the laminated silicon steel. The generalised form of Hooke’s law (equation (10)) can be used to describe the relationship between elastic stress and strain.(10)σe=Dεe,(11)Tmax=BrdHrd−12BHBrdHtdBrdHldBtdHrdBtdHtd−12BHBtdHldBldHrdBldHtdBldHld−12BHnrdntdnld,(12)D=E1−υ1+υ1−2υ⋅1υ1−υυ1−υ000υ1−υ1υ1−υ000υ1−υυ1−υ10000001−2υ21−υ0000001−2υ21−υ0000001−2υ21−υ,(13)εe=εrdεtdεldγrtγtlγlr,where D is the elastic tensor of core material (equation (12)) and εe is the strain in the core material due to magnetostriction (equation (13)). E and υ are Young’s modulus and Poisson’s ratio of the core material, respectively; εrd, εtd, and εld are the principal strains in the material; and γrt, γrt, and γtl are tangential strains. Regardless of the shear component, as in the following equation:(14)ε=εrdεtdεld=λ=λrdλtdλld,where(15)λ=λrdλtdλld=krd000ktd000kldBrdBtdBld.k
rd and ktd are experimentally measured magnetostrictive coefficients of the silicon steel laminate in the RD and TD. The magnetostrictive curves of a type of nonoriented silicon steel were obtained by tests [8] (Figure 4). The magnetostriction coefficient of the LD is currently not available. The magnetostriction of silicon steel is an elastic deformation in which the total volume of the material remains constant [14, 15]. Hence,(16)1+λrd1+λtd1+λld=1.Figure 4
The magnetostrictive curves of the laminated silicon steel.The vibration differential equation of the reactor core is given by(17)M∂2u∂t2+C∂u∂t+Ku=Ftotal,where M is the mass matrix of the core, C is the damping matrix, K is the stiffness matrix, u is the displacement vector, and Ftotal is the external force acted on the core, which is the resultant force of the Maxwell force and the magnetostrictive stress.(18)Ftotal=Fmax+Fmag,Fmag=−∇⋅σ.Regardless of the damping of the vibration process, equation (17) can be simplified to(19)M∂2u∂t2+Ku=Fmax+Fmag.


## 2.1. Laminated Coordinate Systems

The vibration of the core (iron yoke and core limb) is one of the most important causes of the shunt reactor vibration. The core of the shunt reactor is composed of laminated silicon steel sheets, whose magnetic properties are related to the direction of lamination [6]: there are different lamination directions in different parts of the core. According to the rules governing the lamination of silicon steel sheets, three mutually perpendicular axes including RD, TD, and LD can be set; on this basis, a set of 3D coordinate systems is created. The anisotropy of magnetic properties for laminated silicon steel can be expressed more clearly thereby.There are two types of shunt reactor, dry type with small capacity in Figure1(a) and oil-immersed type with large capacity in Figure 1(b) [3, 7]. Core limb and iron yoke are the two main parts of the reactor core. The core limb includes core discs and air gaps (marble or ceramics). The core discs for small dry reactors are usually cuboid, and core discs for large oil-immersed shunt reactors are generally cylindrical. The mapping between the laminated coordinate systems and the spatial rectangular coordinate system is as shown in Table 1.Figure 1
Simplified diagrams of two types of reactors: (a) small dry reactor and (b) large oil-immersed shunt reactor.Table 1
The mapping between the two coordinate systems.
Part
RD
TD
LD
Large oil-immersed shunt reactors
Core disc
z-axis
Vertical to radial andz-axis
Radial
Top and bottom yokes
y-axis
z-axis
x-axis
Side yokes
z-axis
y-axis
x-axis
Small dry reactors
Core discs
y-axis
x-axis
z-axis
Yokes
x-axis
y-axis
z-axisThe lamination rules of the cores at different positions are different, so the number of lamination coordinate systems is large, which entails onerous computations. It is necessary to define a global coordinate system (a space rectangular coordinate system here), which is used to unify parameters of different coordinate systems. Any vector has two different representations belonging to the two coordinate systems. The mapping between these two representations is beneficial to the calculation of the vibration parameters of the reactor.


### 2.1.1. Large Oil-Immersed Shunt Reactors

As shown in Table1, for any vector M (in 3D coordinates) in the iron yokes of large oil-immersed shunt reactors, the mapping between the two coordinate systems satisfying equation (1) (top and bottom yokes) and equation (2) (side yokes) is as follows:(1)Mld=Mx,Mtd=Mz,Mrd=My,(2)Mld=Mx,Mtd=My,Mrd=Mz.The core disc (Figure2) of a large oil-immersed shunt reactor is generally formed by splicing a plurality of fan-shaped laminated silicon steels. The fan-shaped laminated silicon steel is formed by stacking a plurality of rectangular silicon steel sheets with different lengths. The laminated coordinate system is related to the location of each sheet, and different laminated silicon steels belong to different 3D laminated coordinate systems. Taking the fanned laminated silicon steel marked with the border in Figure 2 as an example, the dotted lines indicate the LD (coinciding with the centreline of the sector), and the TD of the laminated coordinate system in which the laminated silicon steel is located. The solid lines indicate the x and y-axes in a 3D space rectangular coordinate system. The RD is recombined with the z-axis perpendicular to the core disc plane. The angle between the x-axis and the LD is α. The mapping of any vector M (in 3D coordinates) in the core discs between the two coordinate systems is as follows:(3)Mld=Mxcosα+Mysinα,Mtd=−Mxsinα+Mycosα,Mrd=Mz.Figure 2
A plan view of a core disc: the green arrow is the coordinate axis of the space rectangular coordinate system, and the red arrow represents the coordinate axis of the laminated coordinate system.


### 2.1.2. Small Dry Reactors

As shown in Table1, the mapping for dry reactors can be performed by equations (4) (yokes) and (5) (core discs).(4)Mld=Mz,Mtd=My,Mrd=Mx,(5)Mld=Mz,Mtd=Mx,Mrd=My,where Mld, Mtd, and Mrd are the 3D components of the vector M in the lamination coordinate system, and Mx, My, and Mz are the 3D components of the vector M in the 3D space rectangular coordinate system in equations (1) to (5).


## 2.1.1. Large Oil-Immersed Shunt Reactors

As shown in Table1, for any vector M (in 3D coordinates) in the iron yokes of large oil-immersed shunt reactors, the mapping between the two coordinate systems satisfying equation (1) (top and bottom yokes) and equation (2) (side yokes) is as follows:(1)Mld=Mx,Mtd=Mz,Mrd=My,(2)Mld=Mx,Mtd=My,Mrd=Mz.The core disc (Figure2) of a large oil-immersed shunt reactor is generally formed by splicing a plurality of fan-shaped laminated silicon steels. The fan-shaped laminated silicon steel is formed by stacking a plurality of rectangular silicon steel sheets with different lengths. The laminated coordinate system is related to the location of each sheet, and different laminated silicon steels belong to different 3D laminated coordinate systems. Taking the fanned laminated silicon steel marked with the border in Figure 2 as an example, the dotted lines indicate the LD (coinciding with the centreline of the sector), and the TD of the laminated coordinate system in which the laminated silicon steel is located. The solid lines indicate the x and y-axes in a 3D space rectangular coordinate system. The RD is recombined with the z-axis perpendicular to the core disc plane. The angle between the x-axis and the LD is α. The mapping of any vector M (in 3D coordinates) in the core discs between the two coordinate systems is as follows:(3)Mld=Mxcosα+Mysinα,Mtd=−Mxsinα+Mycosα,Mrd=Mz.Figure 2
A plan view of a core disc: the green arrow is the coordinate axis of the space rectangular coordinate system, and the red arrow represents the coordinate axis of the laminated coordinate system.


## 2.1.2. Small Dry Reactors

As shown in Table1, the mapping for dry reactors can be performed by equations (4) (yokes) and (5) (core discs).(4)Mld=Mz,Mtd=My,Mrd=Mx,(5)Mld=Mz,Mtd=Mx,Mrd=My,where Mld, Mtd, and Mrd are the 3D components of the vector M in the lamination coordinate system, and Mx, My, and Mz are the 3D components of the vector M in the 3D space rectangular coordinate system in equations (1) to (5).


## 2.2. Vibration Analysis Method

The vibration of the core is caused by the Maxwell force and magnetostriction in the magnetic field, which are affected by the spatial distribution of the magnetic field.In the spatial rectangular coordinate system, the 3D electromagnetic equation of the reactor core under normal operating conditions is given by(6)∂∂yυz∂Ay∂x−∂Ax∂y−∂∂zυy∂Ax∂z−∂Az∂x=Jx,∂∂zυx∂Az∂y−∂Ay∂z−∂∂xυz∂Ay∂x−∂Ax∂y=Jy,∂∂xυy∂Ax∂z−∂Az∂x−∂∂yυx∂Az∂y−∂Ay∂z=Jz,where υx, υy, and υz are the components of the magnetoresistance ratio υ; Ax, Ay, and Az are the components of the vector magnetic potential A; and Jx, Jy, and Jz are the components of the winding current density J.The magnetic flux densityB in the spatial rectangular coordinate system can be calculated by(7)B=∇×A.The magnetic flux densityB in the lamination coordinate system was obtained through coordinate transformation (equations (1) to (5)) of B in the spatial rectangular coordinate system. The magnetic field strength H was calculated by the 3D electromagnetic field constitutive relationship as follows:(8)H=HrdHtdHld=υrd000υtd000υldBrdBtdBld.The 3D magnetisation curves (H-B curves, where H is the magnetic field strength in the lamination coordinate system and H is the magnetic flux density in the laminated coordinate system) of a type of laminated silicon steel were obtained by tests [10] (Figure 3).Figure 3
The 3D magnetisation curves (H-B curves) of the laminated silicon steel.The Maxwell forceFmax acting on the interfaces between core discs and air gaps in the magnetic field can be obtained from equation (9) [11–13]:(9)Fmax=∫sTmaxdS,where Tmax is the Maxwell stress tensor (equation (11)).The magnetostriction phenomenon is caused by the displacement of the magnetic domain inside the core (silicon steel) in the magnetic field, and its magnitude is closely related to the magnetic field. Magnetostrictive displacement and strain are obtained from the magnetic field strength by the magnetostrictive curves of the laminated silicon steel. The generalised form of Hooke’s law (equation (10)) can be used to describe the relationship between elastic stress and strain.(10)σe=Dεe,(11)Tmax=BrdHrd−12BHBrdHtdBrdHldBtdHrdBtdHtd−12BHBtdHldBldHrdBldHtdBldHld−12BHnrdntdnld,(12)D=E1−υ1+υ1−2υ⋅1υ1−υυ1−υ000υ1−υ1υ1−υ000υ1−υυ1−υ10000001−2υ21−υ0000001−2υ21−υ0000001−2υ21−υ,(13)εe=εrdεtdεldγrtγtlγlr,where D is the elastic tensor of core material (equation (12)) and εe is the strain in the core material due to magnetostriction (equation (13)). E and υ are Young’s modulus and Poisson’s ratio of the core material, respectively; εrd, εtd, and εld are the principal strains in the material; and γrt, γrt, and γtl are tangential strains. Regardless of the shear component, as in the following equation:(14)ε=εrdεtdεld=λ=λrdλtdλld,where(15)λ=λrdλtdλld=krd000ktd000kldBrdBtdBld.k
rd and ktd are experimentally measured magnetostrictive coefficients of the silicon steel laminate in the RD and TD. The magnetostrictive curves of a type of nonoriented silicon steel were obtained by tests [8] (Figure 4). The magnetostriction coefficient of the LD is currently not available. The magnetostriction of silicon steel is an elastic deformation in which the total volume of the material remains constant [14, 15]. Hence,(16)1+λrd1+λtd1+λld=1.Figure 4
The magnetostrictive curves of the laminated silicon steel.The vibration differential equation of the reactor core is given by(17)M∂2u∂t2+C∂u∂t+Ku=Ftotal,where M is the mass matrix of the core, C is the damping matrix, K is the stiffness matrix, u is the displacement vector, and Ftotal is the external force acted on the core, which is the resultant force of the Maxwell force and the magnetostrictive stress.(18)Ftotal=Fmax+Fmag,Fmag=−∇⋅σ.Regardless of the damping of the vibration process, equation (17) can be simplified to(19)M∂2u∂t2+Ku=Fmax+Fmag.


## 3. Comparison of 2D and 3D Methods

2D methods have been widely used to study the vibration characteristics of small dry reactor cores, but the new 3D method also has a wide range of applications. Two sets of comparative studies were conducted to compare the differences between the two methods in analysing the vibration of small dry reactors and large oil-immersed shunt reactors. The magnetic flux density (B) and acceleration (ac) of the core vibration were selected as comparison indicators. The difference rate (η) was defined and used to compare the size of the difference between the two methods, as given by(20)ηX=X2D−X3DX3D×100%,where X refers to two contrast indicators (B and ac).In particular, the 2D method used here stems from other work [3, 8]. The parameters of reactors are provided by China TBEA Hengyang Transformer Co., Ltd and China Xi’an Xidian Transformer Co., Ltd (Table 2). The bottom of the reactor is generally connected to the ground and is set as a fixed boundary.Table 2
Parameters of two types of reactors.
Parameter
Large oil-immersed shunt reactor
Small dry reactor
Length of iron yokemm
2980
270
Width of iron yokemm
3045
60
Thickness of iron yokemm
1460
78
Height of windingmm
2650
142
Width of windingmm
372
14
Radius of windingmm
955
97.58
Height of core limbmm
2547
150
Winding turns
2147
10
Density of the core (kg/m3)
7650
7650
Young’s modulus of the core (GPa)
195
195
Poisson’s ratio of the core
0.27
0.27
Young’s modulus of gaps (GPa)
55
55
Poisson’s ratio of gaps
0.45
0.45
Height of an air gapmm
38
28
Number of air gaps
24
4
Height of a cakemm
80
78
Radius of core limbmm
405
\
Width of core limbmm
\
60
Simulation condition
Frequency domain
Frequency domain
Rated current (A)
377.9
300
Current frequency (Hz)
50
50The 2D and 3D methods were used separately, and the differences therein (Table3) were obtained. In terms of the calculation results of the small dry reactor, the distributions (Figure 5) of the two sets of comparative indicators obtained by the two methods are similar. The difference (4.5% and 10.8%) of the average of the two sets of comparative indicators is also practicable (Table 3). Therefore, the two methods are mutually replaceable for the study of small dry reactors. An experiment on the core vibration of a small dry reactor was carried out in reference [3], which verified the 2D method. The difference between the 3D method given here and the 2D method is small for the vibration parameters of the small dry reactor. Therefore, the calculation result of the 3D method proposed here is also correct; however, the 2D method and the 3D method are not equivalent for the large oil-immersed shunt reactors (Figure 6) because the differences (33.1% and 154.14%) in the two methods are too large to be ignored (Table 3). Compared with the 2D method, the 3D method yields the same analytical results for small dry reactors and is more accurate for large oil-immersed reactors. Therefore, the new 3D method should be adopted instead of the traditional 2D method for the large oil-immersed shunt reactors.Table 3
The comparison of 2D and 3D methods.
B
T
ac
m
/
s
2
Small dry reactor
2D
0.21
2.96
×
10
−
3
3D
0.22
3.32
×
10
−
3
η
%
4.5
10.8
Large oil-immersed shunt reactor
2D
1.57
16.57
3D
1.18
6.52
η
%
33.1
154.14Figure 5
The distributions of comparative indicators for small dry reactor. (a) The magnetic flux density. (b) The acceleration. The distributions of magnetic flux density and vibration acceleration of the small dry reactor core obtained by the two calculation methods are similar.
(a)
(b)Figure 6
The distributions of comparison indicators for large oil-immersed shunt reactor. (a) The magnetic flux density. (b) The acceleration. The distributions of magnetic flux density and vibration acceleration of the large oil-immersed reactor core obtained by the two calculation methods are very different.
(a)
(b)


## 4. Influence of Air Gap Number on UHV Shunt Reactor Core Vibration

Since 2006, China has established several UHV (1100 kV) transmission lines for long-distance, large-capacity transmission [16]. A UHV shunt reactor is an important item of equipment in UHV AC transmission projects. Excessive vibration is one of the most common problems in the operation thereof due to the air gaps in iron beam that provide large inductance and the large rated capacity, which is seven to eight times that of a 500 kV reactor [17]. Adjusting the air gaps can be a useful method to reduce vibrations therein. The influence of air gap number on core vibration of the UHV shunt reactor is of interest here. A UHV shunt reactor is a typical oil-immersed shunt reactor, so the new 3D method was adopted.The number of air gap is not allowed to vary widely because of other performance requirements of large oil-immersed shunt reactors. So, the number of air gaps (g) in the core limb is set to 20, 22, 24, 26, and 28. The number of air gaps affects the inductance of the reactor, which is one of the most important parameters of the reactor. The total length of the air gap is adjusted so that the inductance values remain the constant. The calculation method of the reactor inductance is from the reference [18], and the corresponding relationship between the number of air gaps and the total length of the air gap is shown in Table 4. The other reactor parameters are as listed in Table 2.Table 4
Total length of air gap with differentg.
g
20
22
24
26
28
d (mm)
952.52
816.18
707
628.75
566.82


### 4.1. Results of UHV Shunt Reactor Core Vibration

The displacement (d) and acceleration (ac) at each point of the core during its vibration were selected to characterise the vibration strength of the reactor core according to equation (19).The method of 3D finite element simulation was used, and the distribution map ofd was obtained, which shows a gradual increase with increasing number of air gaps (Figure 7). The average displacement of the core, the iron yoke, and the core disc with different air gap numbers was also calculated (Figure 8). They all show increase with increasing air gap numbers. Taking the average value of the core as an example, d=4.58×10−2mm when the air gap number is 20 and when g=24, d=6.06×10−2mm, which is approximately 1.5 times the value at g=20, the maximum value of d is 8.77×10−2mm, at g=28. The variation of the mean average value of d is similar in the iron yoke and the core disc.Figure 7
The distribution of vibration displacement (d) with different air gap numbers. The vibration displacement of the reactor core increases as the number of air gaps increases.Figure 8
The average of displacement (d) with different air gap numbers. The vibration displacement of each component (iron and cakes) of the reactor core increases as the number of air gaps increases.The distribution of vibration acceleration and the average value were also calculated with different air gap numbers (Figures9 and 10). The trend in acceleration is similar to that in the displacement (i.e., both increase).Figure 9
The distribution of vibration acceleration (ac) with different air gap numbers. The vibration acceleration of the reactor core increases as the number of air gaps increases.Figure 10
The average of acceleration (ac) with different air gap numbers. The vibration acceleration of each component (iron and cakes) of the reactor core increases as the number of air gaps increases.Both the displacement and the acceleration of the core vibration increase with the increase in the air gap number. It can be concluded that the increase in the air gap number will aggravate the vibration of the reactor core while maintaining the overall size (length, width, and height) and inductance of the core. In the future reactor design, the minimum number of core air gaps should be chosen to reduce core vibration.


### 4.2. Discussion of UHV Shunt Reactor Core Vibration

The stress on the core is the main cause of the vibration of the reactor core. The acceleration and displacement of the core vibration are positively related to the stress on it. Defining the average value of stress (σaverage) on the core as(21)σaverage=FtotalStotal=∫σtotaldsStotal,where S is total surface area of the core.σ
average on the core was calculated as shown in Figure 11. It increases with the increase in the number of air gaps, which causes the increase in core vibration indicators.Figure 11
σ
average on the core with different air gap numbers. The average value of stress (σaverage) increases as the number of air gaps increases.Magnetic flux density is one of the important factors affecting the stress on the core in the magnetic field. The distribution and average value of magnetic flux density in the reactor with different air gap numbers are as shown in Figures12 and 13. Magnetic flux density increases slightly with increasing air gap numbers. The inductance of the reactor decreases with the increase of the number of air gaps, on the premise that the magnetic flux density and the total length of air gap remain unchanged [18]. Here, the total length of air gap is reduced to maintain the inductance of the reactor at a constant level, when the number of air gaps increases.Figure 12
The distribution of flux density (B) with different air gap numbers. The magnetic flux density of the core increases as the number of air gaps increases.Figure 13
The average value of flux density (B) with different air gap numbers. The average value of the magnetic flux density increases as the number of air gaps increases.The magnetic flux can be divided into two parts: the main magnetic flux and the leakage magnetic flux. The material in the air gap of a reactor is generally marble or a ceramic. The permeability of these two materials is much lower than that of silicon steel, which is an important reason for the observed magnetic flux leakage. The leakage of reactor core is reduced when the length of air gap decreases, so the average magnetic flux density increases with the increase of air gap number.The force on the core is the resultant force of the Maxwell force on the core discs and the magnetostrictive force on the core. The Maxwell force is greater than the magnetostrictive force [9]. The Maxwell force is only applied to the core discs, and it is a surface stress that is related to the unit area stress and the total area acted upon by the applied force. The Maxwell force per unit area is positively correlated with the magnetic flux density. The total stress area of the core and the magnetic flux density both increase with the increase in the air gaps. So, the total Maxwell force increases as the number of air gaps increases. The magnetostrictive force can be calculated by the magnetic flux density and the magnetostriction coefficient, and the relationship between the magnetostriction coefficient and the magnetic flux density is shown in Figure 4. The magnetostrictive force also increases due to the increase in magnetic flux density. Therefore, the increase in stress on the core is attributed to the increase in magnetic flux density and the total force area upon which the Maxwell force acts.


## 4.1. Results of UHV Shunt Reactor Core Vibration

The displacement (d) and acceleration (ac) at each point of the core during its vibration were selected to characterise the vibration strength of the reactor core according to equation (19).The method of 3D finite element simulation was used, and the distribution map ofd was obtained, which shows a gradual increase with increasing number of air gaps (Figure 7). The average displacement of the core, the iron yoke, and the core disc with different air gap numbers was also calculated (Figure 8). They all show increase with increasing air gap numbers. Taking the average value of the core as an example, d=4.58×10−2mm when the air gap number is 20 and when g=24, d=6.06×10−2mm, which is approximately 1.5 times the value at g=20, the maximum value of d is 8.77×10−2mm, at g=28. The variation of the mean average value of d is similar in the iron yoke and the core disc.Figure 7
The distribution of vibration displacement (d) with different air gap numbers. The vibration displacement of the reactor core increases as the number of air gaps increases.Figure 8
The average of displacement (d) with different air gap numbers. The vibration displacement of each component (iron and cakes) of the reactor core increases as the number of air gaps increases.The distribution of vibration acceleration and the average value were also calculated with different air gap numbers (Figures9 and 10). The trend in acceleration is similar to that in the displacement (i.e., both increase).Figure 9
The distribution of vibration acceleration (ac) with different air gap numbers. The vibration acceleration of the reactor core increases as the number of air gaps increases.Figure 10
The average of acceleration (ac) with different air gap numbers. The vibration acceleration of each component (iron and cakes) of the reactor core increases as the number of air gaps increases.Both the displacement and the acceleration of the core vibration increase with the increase in the air gap number. It can be concluded that the increase in the air gap number will aggravate the vibration of the reactor core while maintaining the overall size (length, width, and height) and inductance of the core. In the future reactor design, the minimum number of core air gaps should be chosen to reduce core vibration.


## 4.2. Discussion of UHV Shunt Reactor Core Vibration

The stress on the core is the main cause of the vibration of the reactor core. The acceleration and displacement of the core vibration are positively related to the stress on it. Defining the average value of stress (σaverage) on the core as(21)σaverage=FtotalStotal=∫σtotaldsStotal,where S is total surface area of the core.σ
average on the core was calculated as shown in Figure 11. It increases with the increase in the number of air gaps, which causes the increase in core vibration indicators.Figure 11
σ
average on the core with different air gap numbers. The average value of stress (σaverage) increases as the number of air gaps increases.Magnetic flux density is one of the important factors affecting the stress on the core in the magnetic field. The distribution and average value of magnetic flux density in the reactor with different air gap numbers are as shown in Figures12 and 13. Magnetic flux density increases slightly with increasing air gap numbers. The inductance of the reactor decreases with the increase of the number of air gaps, on the premise that the magnetic flux density and the total length of air gap remain unchanged [18]. Here, the total length of air gap is reduced to maintain the inductance of the reactor at a constant level, when the number of air gaps increases.Figure 12
The distribution of flux density (B) with different air gap numbers. The magnetic flux density of the core increases as the number of air gaps increases.Figure 13
The average value of flux density (B) with different air gap numbers. The average value of the magnetic flux density increases as the number of air gaps increases.The magnetic flux can be divided into two parts: the main magnetic flux and the leakage magnetic flux. The material in the air gap of a reactor is generally marble or a ceramic. The permeability of these two materials is much lower than that of silicon steel, which is an important reason for the observed magnetic flux leakage. The leakage of reactor core is reduced when the length of air gap decreases, so the average magnetic flux density increases with the increase of air gap number.The force on the core is the resultant force of the Maxwell force on the core discs and the magnetostrictive force on the core. The Maxwell force is greater than the magnetostrictive force [9]. The Maxwell force is only applied to the core discs, and it is a surface stress that is related to the unit area stress and the total area acted upon by the applied force. The Maxwell force per unit area is positively correlated with the magnetic flux density. The total stress area of the core and the magnetic flux density both increase with the increase in the air gaps. So, the total Maxwell force increases as the number of air gaps increases. The magnetostrictive force can be calculated by the magnetic flux density and the magnetostriction coefficient, and the relationship between the magnetostriction coefficient and the magnetic flux density is shown in Figure 4. The magnetostrictive force also increases due to the increase in magnetic flux density. Therefore, the increase in stress on the core is attributed to the increase in magnetic flux density and the total force area upon which the Maxwell force acts.


## 5. Conclusion

A new three-dimensional (3D) analysis method is proposed. First, a set of 3D laminated coordinate systems was defined based on silicon steel lamination rules, in which the anisotropy of magnetic properties for laminated silicon steel in the rolling direction (RD), the transverse direction (TD), and the lamination direction (LD) were considered. Then, the mapping between laminated coordinate systems and space rectangular coordinate system was established to unify the parameters in different laminated coordinate systems. After mapping, the anisotropy of the magnetic properties in the laminated coordinate systems was transformed into a spatial rectangular coordinate system.The 3D method was compared with the traditional 2D method by two sets of comparative studies. Applied to small dry reactors, the two methods differ by less than 10.8%. The two methods are mutually replaceable for the study of small dry reactors. The difference is greater than 33.1% when they are used for large oil-immersed shunt reactors, so the new 3D method should be adopted instead of the traditional 2D method for large oil-immersed shunt reactors.The 3D method was used for a UHV shunt reactor, which is a type of a large oil-immersed shunt reactor. On the premise that the inductance of the reactor is constant, the increase in core air gap number will increase the displacement and acceleration of the core vibration. With the increase of the air gap number in the reactor, the flux density in the core increases to keep the inductance of the reactor constant. The total force area upon which the Maxwell force on the core acts also increases with increasing air gap number. They are the two reasons for the intensification of core vibration. To reduce the vibration of the reactor core, the minimum number of air gaps should be selected in the design of the reactor core.



---

*Source: 7290536-2019-10-23.xml*