# Vibration Analysis and Experimental Study of GIS Busbar Enclosure under Electric Force

Authors: Dingge Yang; Yuzhao Li; Yanhua Han; Bin Ding; Liangshu Li; Jingfeng Wu
Journal: Shock and Vibration
 (2022)

---

## Abstract

To explore the vibration response of the GIS busbar enclosure in a strong electric field, the electric force on the busbar enclosure was solved by the voltage in the circuit and the principle of virtual work. The vibration mode was obtained by finite element technology. Accordingly, the vibration response of the busbar enclosure was obtained by vibration analysis in the frequency domain, and the vibration acceleration of the busbar was monitored in the substation. The simulation results showed that the busbar enclosure was subjected to the electric force because of the fundamental voltage and harmonic voltages in the conductor. By coupling the radial electric force and the vibration mode, it was found that the vibration responses of the busbar enclosure at 100 Hz and 2900 Hz were greater than those at other frequencies. The experimental result showed that the frequency with the highest vibration amplitude was the same as that of the simulation result, and the vibration acceleration amplitude in the experiment and simulation at 2900 Hz was basically consistent, which verified the accuracy of the simulation result. The study shows that the GIS enclosure shell will produce vibration at the two times power frequency and high frequency, which provides an explanation for the high-frequency vibration of busbar shell structure in the ultra-high-voltage substation and provides an objective basis for the design of the busbar enclosure in the ultra-high-voltage substation.


---

## Body

## 1. Introduction

Gas-insulated switchgear (GIS) is widely used in high-voltage power system substations because of its compact size and high reliability [1]. A busbar is a device in the GIS substation that can transmit electricity. The busbar enclosure generates mechanical vibration under different excitation during the operation of GIS [2, 3]. The high-frequency vibration of the equipment produces audible noise [4, 5], which can affect the physical and mental health of the staff. Low-frequency vibrations can loosen bolts, causing gas leaks and accidents [6]. It is of great engineering guiding significance to analyze the vibration of the busbar enclosure.There are many factors that cause the GIS busbar enclosure to vibrate. Many scholars studied the influence of partial discharges [7], poor contact [8–10], and temperature distribution change [11, 12] on the busbar enclosure vibration. These fault conditions can cause the failure of the substation operation, so the research studies are helpful to study the reliability and life of the equipment. However, the main excitation in the substation is current and voltage, which leads to the vibration of the electric equipment in normal operating condition. Many scholars studied the vibration responses of the busbar under the action of the electromagnetic force generated by the power frequency current [13–16], and they considered that the electromagnetic force is the main cause of the busbar vibration. Bian et al. [17] explored the vibration response of the GIS busbar under the condition of eccentricity using power frequency voltage. Zhang et al. [18] considered that the vibration is caused by current and voltage excitation at the same time and proposed that it is easier for the busbar enclosure to vibrate in no load because of the high voltage in the ultra-high-voltage substation. It is known that the electric force increases as the voltage rises, and the electric force is also related to the size of the busbar [18]. These research studies mainly focus on the vibration at 100 Hz. However, there are few research studies on the generation principle of high-frequency vibration of power equipment under normal operating conditions. Through the on-site test of the fault-free substation, the high-frequency vibration of the busbar shell appears at 2900 Hz, the cause of which is unknown. This phenomenon has not been encountered before. Therefore, it is needed to come up with an explanation for this phenomenon to provide an objective basis for the design of the busbar enclosure in the ultra-high-voltage substation.This paper proposes an interpretation method to explain the cause of high-frequency vibration of the busbar shell. This paper is organized as follows. Section1 first describes the busbar model and the radial electric force of the enclosure. Second, the vibration responses in the simulation and experiment are obtained in Section 2. The conclusions are obtained in Section 3. The results show that the vibration responses in the simulation and experiment are relatively consistent, indicating that the vibration of GIS busbar enclosure was caused by harmonic voltages in the circuit under the condition of no fault, which provides an explanation for the high-frequency vibration of busbar shell structure in the ultra-high-voltage substation.


### 1.1. Model of the GIS Busbar

#### 1.1.1. Busbar Model and Modal Analysis

The busbar model adopted in this paper is a three-phase divisional box used in 750 kV GIS substation in northwestern China. The busbar model is shown in Figure1, and the explosion model of the busbar is shown in Figure 2. The enclosure is outside the busbar, and the inside is a conductor. The two parts are connected by two basin-type insulators. The busbar enclosure is filled with SF6 gas inside and air outside. The material of the busbar enclosure and conductor is aluminum alloy 6063-T6. Young’s modulus of elasticity E = 69 GPa, density ρ = 2700 kg/m3, and Poisson’s ratio μ = 0.33.Figure 1
Busbar model.Figure 2
Explosion diagram of the busbar model.The size of the busbar enclosure is shown in Table1.Table 1
Dimension parameters of the busbar enclosure.
Conductor’s outer diameterd1 (mm)Conductor’s inner diameterd2 (mm)Busbar enclosure’s outer diameterD1 (mm)Busbar enclosure’s inner diameterD2 (mm)Busbar lengthL (mm)1801506005847500The busbar enclosure was simplified to shell structure to facilitate the calculation. The parts simplified had less influence on the dynamic characteristics of the busbar enclosure. The busbar shell is shown in Figure3.Figure 3
Structure model of the busbar shell.In Figure3, it can be known that the size transition characteristic of the busbar enclosure is preserved, which is essential to the dynamic behaviour of the busbar shell. The model can be used to analyze the dynamic feature of the busbar enclosure.


#### 1.1.2. Radial Electric Force

The voltage excitation state of GIS is the state of loading voltage but no current. The conductor is excited by fundamental voltage and harmonic voltages, busbar shell grounding. There is a potential difference between the conductor and busbar shell because of the no-load state of GIS, which causes electric field inside the busbar. Then, there are induced charges on the inner face of the shell, which are different from the charges on the conductor. The charges on the conductor and busbar shell cause the busbar to be subjected to electric force.The electric force was solved based on five basic assumptions: (1) ignore the charge relaxation—the charge relaxation time of good conductor is much less than 1; (2) no charge inside the conductor—the charge is only distributed on the surface of the conductor; (3) ignore the magnetic field generated by the displacement current; (4) ignore the electric field generated by the polarization charge of SF6; and (5) ignore the influence of shell and conductor deformation on electric field distribution. As a long cylindrical shell structure, it is simplified to a two-dimensional cylindrical shell structure, as shown in Figure 4, which can be used to calculate the radial electric force on the busbar shell.Figure 4
Two-dimensional busbar structure.The voltage between the two cylindrical conductors per unit length [19] is(1)U=τ2πεlnR2r1,where τ is the charge of the conductor per unit length; ε is dielectric constant of SF6, the relative dielectric constant is 1.002049, and the vacuum dielectric constant is 8.854 × 10−12 F/m; R2 is the inner radius of the busbar shell; and r1 is the outer radius of the conductor. According to the principle of virtual work, the electric force of the busbar shell per unit length can be calculated:(2)F=−πεU2R2lnR2/r12,where negative sign indicated that the direction of the radial electric force F was from the busbar shell to the conductor. We can calculate that the electric force increases as the voltage rises, and the electric force is also related to the size of the busbar. According to the assumption number two, we can know that the electric force on the inner face of the busbar shell is uniformed. The electric force per unit area, namely electric force density, can be obtained:(3)pr=εU22R22lnR2/r12,where U is the potential difference between the busbar shell and the conductor. Equation (3) indicated that the electric force density can be obtained once the internal gas and the size of busbar are determined.In this paper, the electric force density could be obtained using the measured harmonic voltage ratio in the substation. The ratio of first 60 harmonic voltages is shown in Figure5.Figure 5
The first 60 harmonic voltages.Figure5 shows that the power frequency (50 Hz) voltage ratio was 100%, and the harmonic voltage ratio was low (<1%). The total harmonic distortion rate of the first 60 harmonic voltages was 0.518%, the odd distortion rate was 0.507%, and the even distortion rate was 0.107%. There is no clear requirement for the harmonic voltage ratio rate of the 750 kV substation in the standard. According to the standard of the 110 kV substation, the limit value of the total harmonic voltage distortion rate is 2%, the limit value of the odd harmonic distortion rate is 1.6%, and the limit value of the even harmonic distortion rate is 0.8%. It was found that all test results met the requirements. The data can be used to calculate the vibration response of the busbar shell.The voltage in the time domain in conductor can be calculated using the ratio of first 60 harmonic voltages. Harmonic voltage ratio can be expressed as follows [20]:(4)HRUh=UhU1×100%,where U1=750×1/3=433kV and the root mean square value of h-th harmonic voltage can be obtained using equation (4). The voltage superposition in the time domain was carried out using the following equation:(5)Ut=∑h=1602Uhsinh×2πft.The voltage in the time domain can be obtained as shown in Figure6.Figure 6
The voltage in the time domain.Figure6 shows that the frequency of the voltage in the circuit was 50 Hz. The electric force density on the busbar shell can be obtained using the voltage in the time domain, as shown in Figure 7.Figure 7
The electric force density in the time domain.Figure7 shows that the frequency of the electric force density was 100 Hz, which was two times the power frequency. The maximum amplitude of the electric force density was 14.04 Pa and the minimum was 0 Pa. Fourier series expansion was used to transform the electric force density in the time domain to the electric force density in the frequency domain [21]:(6)prt=a02+∑n=160ancosnωt+bnsinωt,where(7)an=20.02∫00.02prtcosnωtdt,n=0,1,2,3,…,60,bn=20.02∫00.02prtsinnωtdt=0,n=1,2,3,…,60,where maximum value of n is 60, that is, the maximum frequency of the harmonic voltages is 3000 Hz. ω = 100π rad/s. The electric force density in the frequency domain can be obtained as follows.Figure8 shows that the maximum of the electric force density was at 0 Hz and 100 Hz, and the amplitude was 7.02 Pa. The amplitudes of the electric force density at other frequencies were mostly below 0.01 Pa. It was known that the busbar shell was subjected to electric force generated by the harmonic voltages of the conductor. We can solve the vibration response using the electric force and mode of the busbar shell.Figure 8
The electric force density in the frequency domain.


## 1.1. Model of the GIS Busbar

### 1.1.1. Busbar Model and Modal Analysis

The busbar model adopted in this paper is a three-phase divisional box used in 750 kV GIS substation in northwestern China. The busbar model is shown in Figure1, and the explosion model of the busbar is shown in Figure 2. The enclosure is outside the busbar, and the inside is a conductor. The two parts are connected by two basin-type insulators. The busbar enclosure is filled with SF6 gas inside and air outside. The material of the busbar enclosure and conductor is aluminum alloy 6063-T6. Young’s modulus of elasticity E = 69 GPa, density ρ = 2700 kg/m3, and Poisson’s ratio μ = 0.33.Figure 1
Busbar model.Figure 2
Explosion diagram of the busbar model.The size of the busbar enclosure is shown in Table1.Table 1
Dimension parameters of the busbar enclosure.
Conductor’s outer diameterd1 (mm)Conductor’s inner diameterd2 (mm)Busbar enclosure’s outer diameterD1 (mm)Busbar enclosure’s inner diameterD2 (mm)Busbar lengthL (mm)1801506005847500The busbar enclosure was simplified to shell structure to facilitate the calculation. The parts simplified had less influence on the dynamic characteristics of the busbar enclosure. The busbar shell is shown in Figure3.Figure 3
Structure model of the busbar shell.In Figure3, it can be known that the size transition characteristic of the busbar enclosure is preserved, which is essential to the dynamic behaviour of the busbar shell. The model can be used to analyze the dynamic feature of the busbar enclosure.


### 1.1.2. Radial Electric Force

The voltage excitation state of GIS is the state of loading voltage but no current. The conductor is excited by fundamental voltage and harmonic voltages, busbar shell grounding. There is a potential difference between the conductor and busbar shell because of the no-load state of GIS, which causes electric field inside the busbar. Then, there are induced charges on the inner face of the shell, which are different from the charges on the conductor. The charges on the conductor and busbar shell cause the busbar to be subjected to electric force.The electric force was solved based on five basic assumptions: (1) ignore the charge relaxation—the charge relaxation time of good conductor is much less than 1; (2) no charge inside the conductor—the charge is only distributed on the surface of the conductor; (3) ignore the magnetic field generated by the displacement current; (4) ignore the electric field generated by the polarization charge of SF6; and (5) ignore the influence of shell and conductor deformation on electric field distribution. As a long cylindrical shell structure, it is simplified to a two-dimensional cylindrical shell structure, as shown in Figure 4, which can be used to calculate the radial electric force on the busbar shell.Figure 4
Two-dimensional busbar structure.The voltage between the two cylindrical conductors per unit length [19] is(1)U=τ2πεlnR2r1,where τ is the charge of the conductor per unit length; ε is dielectric constant of SF6, the relative dielectric constant is 1.002049, and the vacuum dielectric constant is 8.854 × 10−12 F/m; R2 is the inner radius of the busbar shell; and r1 is the outer radius of the conductor. According to the principle of virtual work, the electric force of the busbar shell per unit length can be calculated:(2)F=−πεU2R2lnR2/r12,where negative sign indicated that the direction of the radial electric force F was from the busbar shell to the conductor. We can calculate that the electric force increases as the voltage rises, and the electric force is also related to the size of the busbar. According to the assumption number two, we can know that the electric force on the inner face of the busbar shell is uniformed. The electric force per unit area, namely electric force density, can be obtained:(3)pr=εU22R22lnR2/r12,where U is the potential difference between the busbar shell and the conductor. Equation (3) indicated that the electric force density can be obtained once the internal gas and the size of busbar are determined.In this paper, the electric force density could be obtained using the measured harmonic voltage ratio in the substation. The ratio of first 60 harmonic voltages is shown in Figure5.Figure 5
The first 60 harmonic voltages.Figure5 shows that the power frequency (50 Hz) voltage ratio was 100%, and the harmonic voltage ratio was low (<1%). The total harmonic distortion rate of the first 60 harmonic voltages was 0.518%, the odd distortion rate was 0.507%, and the even distortion rate was 0.107%. There is no clear requirement for the harmonic voltage ratio rate of the 750 kV substation in the standard. According to the standard of the 110 kV substation, the limit value of the total harmonic voltage distortion rate is 2%, the limit value of the odd harmonic distortion rate is 1.6%, and the limit value of the even harmonic distortion rate is 0.8%. It was found that all test results met the requirements. The data can be used to calculate the vibration response of the busbar shell.The voltage in the time domain in conductor can be calculated using the ratio of first 60 harmonic voltages. Harmonic voltage ratio can be expressed as follows [20]:(4)HRUh=UhU1×100%,where U1=750×1/3=433kV and the root mean square value of h-th harmonic voltage can be obtained using equation (4). The voltage superposition in the time domain was carried out using the following equation:(5)Ut=∑h=1602Uhsinh×2πft.The voltage in the time domain can be obtained as shown in Figure6.Figure 6
The voltage in the time domain.Figure6 shows that the frequency of the voltage in the circuit was 50 Hz. The electric force density on the busbar shell can be obtained using the voltage in the time domain, as shown in Figure 7.Figure 7
The electric force density in the time domain.Figure7 shows that the frequency of the electric force density was 100 Hz, which was two times the power frequency. The maximum amplitude of the electric force density was 14.04 Pa and the minimum was 0 Pa. Fourier series expansion was used to transform the electric force density in the time domain to the electric force density in the frequency domain [21]:(6)prt=a02+∑n=160ancosnωt+bnsinωt,where(7)an=20.02∫00.02prtcosnωtdt,n=0,1,2,3,…,60,bn=20.02∫00.02prtsinnωtdt=0,n=1,2,3,…,60,where maximum value of n is 60, that is, the maximum frequency of the harmonic voltages is 3000 Hz. ω = 100π rad/s. The electric force density in the frequency domain can be obtained as follows.Figure8 shows that the maximum of the electric force density was at 0 Hz and 100 Hz, and the amplitude was 7.02 Pa. The amplitudes of the electric force density at other frequencies were mostly below 0.01 Pa. It was known that the busbar shell was subjected to electric force generated by the harmonic voltages of the conductor. We can solve the vibration response using the electric force and mode of the busbar shell.Figure 8
The electric force density in the frequency domain.


## 1.1.1. Busbar Model and Modal Analysis

The busbar model adopted in this paper is a three-phase divisional box used in 750 kV GIS substation in northwestern China. The busbar model is shown in Figure1, and the explosion model of the busbar is shown in Figure 2. The enclosure is outside the busbar, and the inside is a conductor. The two parts are connected by two basin-type insulators. The busbar enclosure is filled with SF6 gas inside and air outside. The material of the busbar enclosure and conductor is aluminum alloy 6063-T6. Young’s modulus of elasticity E = 69 GPa, density ρ = 2700 kg/m3, and Poisson’s ratio μ = 0.33.Figure 1
Busbar model.Figure 2
Explosion diagram of the busbar model.The size of the busbar enclosure is shown in Table1.Table 1
Dimension parameters of the busbar enclosure.
Conductor’s outer diameterd1 (mm)Conductor’s inner diameterd2 (mm)Busbar enclosure’s outer diameterD1 (mm)Busbar enclosure’s inner diameterD2 (mm)Busbar lengthL (mm)1801506005847500The busbar enclosure was simplified to shell structure to facilitate the calculation. The parts simplified had less influence on the dynamic characteristics of the busbar enclosure. The busbar shell is shown in Figure3.Figure 3
Structure model of the busbar shell.In Figure3, it can be known that the size transition characteristic of the busbar enclosure is preserved, which is essential to the dynamic behaviour of the busbar shell. The model can be used to analyze the dynamic feature of the busbar enclosure.


## 1.1.2. Radial Electric Force

The voltage excitation state of GIS is the state of loading voltage but no current. The conductor is excited by fundamental voltage and harmonic voltages, busbar shell grounding. There is a potential difference between the conductor and busbar shell because of the no-load state of GIS, which causes electric field inside the busbar. Then, there are induced charges on the inner face of the shell, which are different from the charges on the conductor. The charges on the conductor and busbar shell cause the busbar to be subjected to electric force.The electric force was solved based on five basic assumptions: (1) ignore the charge relaxation—the charge relaxation time of good conductor is much less than 1; (2) no charge inside the conductor—the charge is only distributed on the surface of the conductor; (3) ignore the magnetic field generated by the displacement current; (4) ignore the electric field generated by the polarization charge of SF6; and (5) ignore the influence of shell and conductor deformation on electric field distribution. As a long cylindrical shell structure, it is simplified to a two-dimensional cylindrical shell structure, as shown in Figure 4, which can be used to calculate the radial electric force on the busbar shell.Figure 4
Two-dimensional busbar structure.The voltage between the two cylindrical conductors per unit length [19] is(1)U=τ2πεlnR2r1,where τ is the charge of the conductor per unit length; ε is dielectric constant of SF6, the relative dielectric constant is 1.002049, and the vacuum dielectric constant is 8.854 × 10−12 F/m; R2 is the inner radius of the busbar shell; and r1 is the outer radius of the conductor. According to the principle of virtual work, the electric force of the busbar shell per unit length can be calculated:(2)F=−πεU2R2lnR2/r12,where negative sign indicated that the direction of the radial electric force F was from the busbar shell to the conductor. We can calculate that the electric force increases as the voltage rises, and the electric force is also related to the size of the busbar. According to the assumption number two, we can know that the electric force on the inner face of the busbar shell is uniformed. The electric force per unit area, namely electric force density, can be obtained:(3)pr=εU22R22lnR2/r12,where U is the potential difference between the busbar shell and the conductor. Equation (3) indicated that the electric force density can be obtained once the internal gas and the size of busbar are determined.In this paper, the electric force density could be obtained using the measured harmonic voltage ratio in the substation. The ratio of first 60 harmonic voltages is shown in Figure5.Figure 5
The first 60 harmonic voltages.Figure5 shows that the power frequency (50 Hz) voltage ratio was 100%, and the harmonic voltage ratio was low (<1%). The total harmonic distortion rate of the first 60 harmonic voltages was 0.518%, the odd distortion rate was 0.507%, and the even distortion rate was 0.107%. There is no clear requirement for the harmonic voltage ratio rate of the 750 kV substation in the standard. According to the standard of the 110 kV substation, the limit value of the total harmonic voltage distortion rate is 2%, the limit value of the odd harmonic distortion rate is 1.6%, and the limit value of the even harmonic distortion rate is 0.8%. It was found that all test results met the requirements. The data can be used to calculate the vibration response of the busbar shell.The voltage in the time domain in conductor can be calculated using the ratio of first 60 harmonic voltages. Harmonic voltage ratio can be expressed as follows [20]:(4)HRUh=UhU1×100%,where U1=750×1/3=433kV and the root mean square value of h-th harmonic voltage can be obtained using equation (4). The voltage superposition in the time domain was carried out using the following equation:(5)Ut=∑h=1602Uhsinh×2πft.The voltage in the time domain can be obtained as shown in Figure6.Figure 6
The voltage in the time domain.Figure6 shows that the frequency of the voltage in the circuit was 50 Hz. The electric force density on the busbar shell can be obtained using the voltage in the time domain, as shown in Figure 7.Figure 7
The electric force density in the time domain.Figure7 shows that the frequency of the electric force density was 100 Hz, which was two times the power frequency. The maximum amplitude of the electric force density was 14.04 Pa and the minimum was 0 Pa. Fourier series expansion was used to transform the electric force density in the time domain to the electric force density in the frequency domain [21]:(6)prt=a02+∑n=160ancosnωt+bnsinωt,where(7)an=20.02∫00.02prtcosnωtdt,n=0,1,2,3,…,60,bn=20.02∫00.02prtsinnωtdt=0,n=1,2,3,…,60,where maximum value of n is 60, that is, the maximum frequency of the harmonic voltages is 3000 Hz. ω = 100π rad/s. The electric force density in the frequency domain can be obtained as follows.Figure8 shows that the maximum of the electric force density was at 0 Hz and 100 Hz, and the amplitude was 7.02 Pa. The amplitudes of the electric force density at other frequencies were mostly below 0.01 Pa. It was known that the busbar shell was subjected to electric force generated by the harmonic voltages of the conductor. We can solve the vibration response using the electric force and mode of the busbar shell.Figure 8
The electric force density in the frequency domain.


## 2. Results and Discussion

### 2.1. Vibration Response in the Simulation

In COMSOL Multiphysics, we can obtain the results by numerical simulation. Applying fixed boundary at the both ends of the busbar shell for modal analysis, the vibration mode of the busbar shell within the range of 0∼3000 Hz can be obtained. Among them, the mode shape of the busbar shell at 2883 Hz was special, which was the expansion-contraction deformation in the middle part, as shown in Figure9.Figure 9
Expansion-contraction mode shape of the busbar shell at 2883 Hz in COMSOL.The mode shape of the busbar shell at 2883 Hz was expansion-contraction deformation. When the frequencies of excitation and mode of the busbar shell were consistent, the busbar shell resonated. The mode shapes at other frequencies were complex and difficult to be excited, so they are not listed here.The vibration response of the busbar shell can be obtained using the electric force density in the frequency domain and the mode of the busbar shell. The voltage is accurate at an integer multiple of 50 Hz, so the frequency interval is 50 Hz, and the frequency range is 50∼3000 Hz in the process of solving. The modal damping loss factor considered in the calculations is 0.005. The measurement point in the simulation is shown in Figure3. The radial vibration response of the busbar shell can be obtained, as shown in Figure 10.Figure 10
Radial vibration response of the busbar shell in COMSOL.Figure10 shows that the vibration amplitudes at 100 Hz and 2900 Hz were greater than those at other frequencies. At 2900 Hz, the amplitude of radial vibration acceleration of the busbar shell was 0.02288 m/s2. The amplitudes at other frequencies were close to zero. The deformation of the busbar shell at 2900 Hz is shown in Figure 11.Figure 11
The acceleration of the busbar shell at 2900 Hz in COMSOL.The reason why the busbar shell vibrated at 100 Hz was that the electric force density at 100 Hz was larger than that at other frequencies. The electric force density at 2900 Hz was very small, but the frequency of expansion-contraction mode was 2883 Hz, which was near 2900 Hz, so the vibration amplitude at 2900 Hz was very large. It can be judged that the vibration at 100 Hz was forced vibration, and the vibration at 2900 Hz was resonance.Consistent results can be also obtained through numerical simulation in ANSYS Workbench. We can get the expansion-contraction mode shape of the busbar shell as shown in Figure12, radial vibration response of the busbar shell as shown in Figure 13, and the acceleration of the busbar shell at 2900 Hz as shown in Figure 14. The results in ANSYS Workbench were the same as the results in COMSOL Multiphysics.Figure 12
Expansion-contraction mode shape of the busbar shell at 2900 Hz in ANSYS Workbench.Figure 13
Radial vibration response of the busbar shell in ANSYS Workbench.Figure 14
The acceleration of the busbar shell at 2900 Hz in ANSYS Workbench.It can be known that the simulation results obtained by COMSOL Multiphysics and ANSYS Workbench are the same, so it can be known that the simulation results are correct. The vibration amplitude at 2900 Hz in COMSOL Multiphysics was 0.02288 m/s2, and the result in ANSYS Workbench was 0.02249 m/s2. The absolute deviation of vibration acceleration amplitude is 0.00039 m/s2, and the relative error is 1.7%, which is very small, so we can conclude that the results are correct.


### 2.2. Vibration Response in the Experiment

The radial vibration acceleration of the busbar enclosure in the actual operating condition was tested, and there was no fault inside the substation. The position of the measuring point is shown in Figure1.The equipment we used is the m + p software and hardware. All of these are provided by m + p international. The sensitivity of the sensor is 103.4 mV/g. The sample rate is 12800 Hz, and the useful bandwidth is 5000 Hz. The block size is 1024. The sample time is three seconds. The equipment in the scene is shown in Figure15.Figure 15
The equipment in the vibration test.The spectrum of radial vibration acceleration in the range of 0∼5000 Hz can be obtained, as shown in Figure16.Figure 16
Radial vibration response of the busbar shell in the test.Figure16 shows that the frequency with the highest vibration amplitude in the experiment was 2900 Hz, which was the same as the simulation result. The accuracy of the simulation result was verified. The simulation result showed that the vibration amplitude at 2900 Hz was 0.02288 m/s2 in COMSOL Multiphysics and 0.02249 m/s2 in ANSYS Workbench. In the experiment, the vibration amplitude at 2900 Hz was 0.02103 m/s2. The relative error between the simulation and the experiment was 8.80% and 6.94%, which was within the acceptable limit, indicating that the simulation was correct to a certain extent.There were several reasons for the error between the simulation and experimental results. (1) The model simplification error: the busbar enclosure was simplified to shell structure, and there may be some differences in dynamic characteristics between the busbar shell and the actual busbar enclosure. (2) The vibration test error: the vibration results may be affected by the external environment in different ways. (3) The harmonic voltage test error: the harmonic voltage data were more accurate near the low harmonic, but the high harmonic data were rough, which led to the error of electric force density. (4) The effect of current in the circuit: in this paper, only voltage was considered, and the magnetic force generated by the current in the operating condition may also affect the vibration response.Generally, the industrial equipment can generate the vibration at low frequencies because of the simple mode. At high frequencies, the mode shape is very complex, and the equipment is difficult to vibrate. However, the structure of the busbar enclosure is thin-shell structure, which has a special mode shape in the vicinity of 2900 Hz. The mode shape is the expansion-contraction deformation as shown in Figures9 and 12, which is similar to the electric force distribution. When the natural frequency of the expansion-contraction mode is close to the frequency of the electric force, the busbar enclosure can generate resonance, and the vibration amplitude is very large.In this case, the natural frequency of the expansion-contraction mode is in the vicinity of 2900 Hz, and the radial electric force at 2900 Hz is large, so the vibration acceleration is very large at 2900 Hz. By comparing the simulation with the experimental results, we can conclude that the electric force generated by the harmonic voltages can excite the mode of the busbar enclosure, causing the busbar enclosure to vibrate. When the amplitude of the electric force is large (100 Hz) or the electric force excitation frequency (2900 Hz) is close to the natural frequency with expansion–contraction deformation, the vibration response is large. We can know that the voltage is not the main factor, but the expansion and contraction mode of busbar shell structure is the main factor. The harmonic voltage content in the substation meets the standard requirements. The vibration occurred mainly because of the expansion and contraction mode at 2900 Hz.Using the same simulation steps, we can verify the vibration response of the product in the electric force to guide the product design. After the 3D model of the busbar shell is established, the mode is solved and the response can be obtained by applying electric field force to sweep frequency, and the existence of high-frequency vibration can be finally checked. If there is high-frequency vibration, the structure needs to be redesigned. If the high-frequency vibration is not obvious, the structure design is qualified.


## 2.1. Vibration Response in the Simulation

In COMSOL Multiphysics, we can obtain the results by numerical simulation. Applying fixed boundary at the both ends of the busbar shell for modal analysis, the vibration mode of the busbar shell within the range of 0∼3000 Hz can be obtained. Among them, the mode shape of the busbar shell at 2883 Hz was special, which was the expansion-contraction deformation in the middle part, as shown in Figure9.Figure 9
Expansion-contraction mode shape of the busbar shell at 2883 Hz in COMSOL.The mode shape of the busbar shell at 2883 Hz was expansion-contraction deformation. When the frequencies of excitation and mode of the busbar shell were consistent, the busbar shell resonated. The mode shapes at other frequencies were complex and difficult to be excited, so they are not listed here.The vibration response of the busbar shell can be obtained using the electric force density in the frequency domain and the mode of the busbar shell. The voltage is accurate at an integer multiple of 50 Hz, so the frequency interval is 50 Hz, and the frequency range is 50∼3000 Hz in the process of solving. The modal damping loss factor considered in the calculations is 0.005. The measurement point in the simulation is shown in Figure3. The radial vibration response of the busbar shell can be obtained, as shown in Figure 10.Figure 10
Radial vibration response of the busbar shell in COMSOL.Figure10 shows that the vibration amplitudes at 100 Hz and 2900 Hz were greater than those at other frequencies. At 2900 Hz, the amplitude of radial vibration acceleration of the busbar shell was 0.02288 m/s2. The amplitudes at other frequencies were close to zero. The deformation of the busbar shell at 2900 Hz is shown in Figure 11.Figure 11
The acceleration of the busbar shell at 2900 Hz in COMSOL.The reason why the busbar shell vibrated at 100 Hz was that the electric force density at 100 Hz was larger than that at other frequencies. The electric force density at 2900 Hz was very small, but the frequency of expansion-contraction mode was 2883 Hz, which was near 2900 Hz, so the vibration amplitude at 2900 Hz was very large. It can be judged that the vibration at 100 Hz was forced vibration, and the vibration at 2900 Hz was resonance.Consistent results can be also obtained through numerical simulation in ANSYS Workbench. We can get the expansion-contraction mode shape of the busbar shell as shown in Figure12, radial vibration response of the busbar shell as shown in Figure 13, and the acceleration of the busbar shell at 2900 Hz as shown in Figure 14. The results in ANSYS Workbench were the same as the results in COMSOL Multiphysics.Figure 12
Expansion-contraction mode shape of the busbar shell at 2900 Hz in ANSYS Workbench.Figure 13
Radial vibration response of the busbar shell in ANSYS Workbench.Figure 14
The acceleration of the busbar shell at 2900 Hz in ANSYS Workbench.It can be known that the simulation results obtained by COMSOL Multiphysics and ANSYS Workbench are the same, so it can be known that the simulation results are correct. The vibration amplitude at 2900 Hz in COMSOL Multiphysics was 0.02288 m/s2, and the result in ANSYS Workbench was 0.02249 m/s2. The absolute deviation of vibration acceleration amplitude is 0.00039 m/s2, and the relative error is 1.7%, which is very small, so we can conclude that the results are correct.


## 2.2. Vibration Response in the Experiment

The radial vibration acceleration of the busbar enclosure in the actual operating condition was tested, and there was no fault inside the substation. The position of the measuring point is shown in Figure1.The equipment we used is the m + p software and hardware. All of these are provided by m + p international. The sensitivity of the sensor is 103.4 mV/g. The sample rate is 12800 Hz, and the useful bandwidth is 5000 Hz. The block size is 1024. The sample time is three seconds. The equipment in the scene is shown in Figure15.Figure 15
The equipment in the vibration test.The spectrum of radial vibration acceleration in the range of 0∼5000 Hz can be obtained, as shown in Figure16.Figure 16
Radial vibration response of the busbar shell in the test.Figure16 shows that the frequency with the highest vibration amplitude in the experiment was 2900 Hz, which was the same as the simulation result. The accuracy of the simulation result was verified. The simulation result showed that the vibration amplitude at 2900 Hz was 0.02288 m/s2 in COMSOL Multiphysics and 0.02249 m/s2 in ANSYS Workbench. In the experiment, the vibration amplitude at 2900 Hz was 0.02103 m/s2. The relative error between the simulation and the experiment was 8.80% and 6.94%, which was within the acceptable limit, indicating that the simulation was correct to a certain extent.There were several reasons for the error between the simulation and experimental results. (1) The model simplification error: the busbar enclosure was simplified to shell structure, and there may be some differences in dynamic characteristics between the busbar shell and the actual busbar enclosure. (2) The vibration test error: the vibration results may be affected by the external environment in different ways. (3) The harmonic voltage test error: the harmonic voltage data were more accurate near the low harmonic, but the high harmonic data were rough, which led to the error of electric force density. (4) The effect of current in the circuit: in this paper, only voltage was considered, and the magnetic force generated by the current in the operating condition may also affect the vibration response.Generally, the industrial equipment can generate the vibration at low frequencies because of the simple mode. At high frequencies, the mode shape is very complex, and the equipment is difficult to vibrate. However, the structure of the busbar enclosure is thin-shell structure, which has a special mode shape in the vicinity of 2900 Hz. The mode shape is the expansion-contraction deformation as shown in Figures9 and 12, which is similar to the electric force distribution. When the natural frequency of the expansion-contraction mode is close to the frequency of the electric force, the busbar enclosure can generate resonance, and the vibration amplitude is very large.In this case, the natural frequency of the expansion-contraction mode is in the vicinity of 2900 Hz, and the radial electric force at 2900 Hz is large, so the vibration acceleration is very large at 2900 Hz. By comparing the simulation with the experimental results, we can conclude that the electric force generated by the harmonic voltages can excite the mode of the busbar enclosure, causing the busbar enclosure to vibrate. When the amplitude of the electric force is large (100 Hz) or the electric force excitation frequency (2900 Hz) is close to the natural frequency with expansion–contraction deformation, the vibration response is large. We can know that the voltage is not the main factor, but the expansion and contraction mode of busbar shell structure is the main factor. The harmonic voltage content in the substation meets the standard requirements. The vibration occurred mainly because of the expansion and contraction mode at 2900 Hz.Using the same simulation steps, we can verify the vibration response of the product in the electric force to guide the product design. After the 3D model of the busbar shell is established, the mode is solved and the response can be obtained by applying electric field force to sweep frequency, and the existence of high-frequency vibration can be finally checked. If there is high-frequency vibration, the structure needs to be redesigned. If the high-frequency vibration is not obvious, the structure design is qualified.


## 3. Conclusions

In this study, based on the actual harmonic voltages in the circuit, the vibration response of the busbar enclosure was explored considering the electric force and mode of the busbar enclosure. By comparing the numerical results with the experimental results, it can be shown that the GIS enclosure shell will produce vibration at the two times power frequency and high frequency. From the results presented in this study, the following conclusions can be drawn:(1)
The busbar enclosure was subjected to electric force due to the presence of voltage in the conductor. Based on the data of the fundamental voltage and harmonic voltage ratio, we can calculate the voltage in the time domain, the electric force density in time domain, and the electric force density in the frequency domain in turn. The excitation on the busbar enclosure can be obtained.(2)
The busbar enclosure vibrated under the action of electric force. When the amplitude of the electric force was large (100 Hz) or the frequency of the electric force (2900 Hz) was close to the frequency with expansion-contraction mode shape, the vibration response was large. The vibration at 100 Hz was forced vibration, and the vibration at 2900 Hz was resonance.(3)
The test results showed that the frequency with the highest vibration amplitude in the experiment was the same as that of the simulation results, and the amplitude between the simulation and test results at 2900 Hz was basically consistent. The comparison results verified the accuracy of the simulation in this paper. It provides an explanation for the high-frequency vibration of busbar shell structure and an objective basis for the design of the busbar enclosure in the ultra-high-voltage substation.



---

*Source: 3055139-2022-01-30.xml*