Insulation Performance and Simulation Analysis of SiO2-Aramid Paper under High-Voltage Bushing

The long-term safe and stable operation of oil-impregnated paper (OIP) bushings is of great significance to the operation of power systems. With the growth of OIP bushing, its internal insulation will gradually decay. Aramid insulation paper has excellent thermal aging characteristics and its insulation performance can be improved by using nano-modification technology. In this paper, the nano-SiO2 particles were used as the modified additives, and the modified aramid insulation paper was prepared through four steps: ultrasonic stirring, fiber dissociation, paper sample copying and superheated calendering. The microscopic physical morphology and chemical components of the insulation specimens before and after modification were analyzed by atomic force microscopy (AFM), scanning electron microscopy (SEM) and X-ray photoelectron spectroscopy (XPS), and an OIP bushing model based on the modified aramid insulation paper was constructed and its electric field distribution was analyzed. The simulation results show that the use of SiO2-modified aramid insulation paper can improve the electric field distribution of OIP bushings and increase the operating life of power transformers.


Introduction
As an important part of power transformers, oil-impregnated paper (OIP) bushings assume the important roles of leading out the high and low voltage leads of power transformers, fixing the internal lead structure, and achieving insulation between the leads and the transformer tank [1]. Currently, OIP bushings for power transformers usually adopt a structure of multi-layer metal laminations, while the oil-impregnated paper insulation is widely used inside power transformers due to its high thermal stability.
In long-term operation of transformer bushings, due to the combined effect of high temperature and strong electric fields, oil-impregnated paper insulation will be aged with the operation of electrical equipment, thus causing a decline in electrical and mechanical properties. In actual operation, aging transformer oil can be filtered or even changed to remove aging products to restore insulation performance, but insulation paper in the bushing cannot be replaced, so the aging state of the insulation paper determines the aging life of the transformer.
The main components of aramid paper are two polymers, aramid short-cut fiber and aramid pulp [2]. It has the excellent properties of aramid fiber, high strength and high modulus, low deformation, high temperature insulation, corrosion resistance, and is an excellent insulating material due to its light weight, flame retardancy, and abrasion resistance. Aramid insulating paper not only has good electrical properties, but also can operate at 180 • C for a long time, has excellent anti-heat aging properties, and can be used safely and stably for more than 10 years [3]. Therefore, replacing the original cable paper in oil-impregnated paper bushing by using aramid paper was considered to improve the thermal aging resistance characteristics of the oil-impregnated paper insulation. However, the very small radial insulation size of the OIP bushings leads to a large spatial electric field strength, and the dielectric constant of the aramid paper needs to be modified to homogenize electric field distribution inside the bushing. At present, domestic and foreign scholars have made relevant results on the application of nano-modification technology to insulating paper .
In the field of using nanoparticles to improve the electrical properties of traditional insulating paper, scholars have tried to modify insulating paper by adding nano-SiO 2 , nano-Al 2 O 3 , nano-TiO 2 , nano-MgO, and nano-ZnO particles, and the performance of these nanoparticles can be different. The nano-SiO 2 particles can reduce the dielectric constant of insulating paper, and at a content of 5%, the dielectric constant is 34% lower than that before modification [12]; the thermal conductivity of insulating paper increases with the addition of nano-Al 2 O 3 compared with the previous one [13,14]; the addition of nano-TiO 2 increases the frequency breakdown voltage of insulating paper [15]; nano-MgO modified insulating paper has better anti-aging performance [19]; nano-ZnO has better suppressive effect on space charge of insulating paper [19]. Some scholars also used nano-SiC to modify insulating paper for the transformer outlet, and the study showed that field strength in the modified insulating paper was reduced under the action of DC and polarity reversal voltage, and the electric field distribution in the oil-impregnated paper tended to be more uniform [16]. The literature [23,24] shows that nanoparticles with low dielectric constants are more beneficial to improve the breakdown strength of composite nanomaterials, and common nanoparticles for cellulose insulating paper modification are mainly nano-Al 2 O 3 , nano-TiO 2 and nano-SiO 2 , whose dielectric constants are 9.8, 110 and 4.5, respectively. Meanwhile, scholars proposed a method to prepare low dielectric constant insulating paper using nano-SiO 2 materials, which also provided ideas for this paper [10]. Therefore, in this paper, nano-SiO 2 particles were used as the modified additive of aramid insulating paper, and the nano-modified aramid insulating paper was prepared based on nano-modification technology, then the modified aramid insulating paper was tested by atomic force microscopy (AFM), scanning electron microscopy (SEM) and X-ray photoelectron spectroscopy (XPS) to further analyze the effect of nano-SiO 2 particles on the physical and chemical properties of the aramid insulating paper. Based on the results, an OIP bushing model was developed in Comsol to analyze the effect of modified aramid insulating paper on the actual electric field distribution of power transformers.

Material Preparation
In this paper, nano-modification technology was used to modify aramid insulating paper. The nano-modified aramid paper samples were prepared using nano-SiO 2 particles as the modified additives of aramid insulating paper through four steps: ultrasonic stirring, fiber dissociation, paper sample copying and superheated calendering.
Raw materials used in the whole modified aramid sample preparation process include aramid staple fiber, aramid pulp and nano-SiO 2 particles, among which aramid staple fiber and aramid pulp, as important components of aramid insulating paper, play a decisive role in the material properties of aromatic insulating paper. The selected aramid short-cut fiber is 12 mm in length, about 10 µm in diameter and light yellow in color; the selected aramid pulp is 2 mm in length and light yellow in color; the selected SiO 2 nanoparticles are about 12 nm in diameter, with a specific surface area of about 200 m 2 /g and are a white solid powder. The aramid short-cut fiber and aramid pulp selected for this paper were purchased from Yantai Teplon Advanced Manufacturing Technology Co. The nano-SiO 2 powder was purchased from Guangzhou Wanjing Co. (Guangzhou Wanjing Co., Guangzhou, China) The whole modified aramid insulating paper sample preparation process is shown in Figure 1. The specific process is as follows: 1. The SiO2 nanoparticles with a diameter of about 12 nm are added to deionized water together with a 12% dispersant solution, and the mixed solution is dispersed ultrasonically with the aid of an ultrasonic cleaner; 2. Take a fiber dissociator, add an appropriate amount of aramid short-cut fibers with a length of 12 mm to it, add the previously prepared nano-SiO2 particles dispersion solution after being dispersed evenly, and stir it fully for 8 min. Subsequently add an appropriate amount of aramid pulp with a length of 2 mm, and stir it fully for 8 min to obtain a mixed original pulp containing nano-SiO2 particles; 3. Use a paper sample copying machine to copy and form the mixed original pulp containing the nano-SiO2 particles to obtain dry paper containing the mixture of nano-SiO2 particles; 4. The obtained dry paper containing the mixture of nano-SiO2 particles is superheated and calendered to obtain an improved aramid insulating paper containing nano-SiO2 particles.

Test Methods
In this paper, a Dimension FastScan atomic force microscope (AFM) (Suzhou Feisman Precision Instruments Co., Suzhou, China) and Quanta FEG 250 scanning electron microscope (SEM) ( Carl Zeiss AG, Oberkochen, Germany)were used to test the external physical morphology of modified aramid insulating paper samples containing SiO2 nanoparticles to investigate the effect of SiO2 nanoparticles on the microscopic morphology of aramid insulating paper. Thermo SCIENTIFIC ESCALAB 250Xi X-ray photoelectron spectroscopy (XPS)(Thermo Fisher Scientific, Waltham, MA, USA) was used to analyze the surface chemistry of the modified aramid insulating paper samples to obtain the composition elements of the modified aramid insulating paper samples.
By controlling the amount of SiO2 nanoparticles added during the preparation of the modified aramid insulating paper samples, the optimal modification ratio was investigated, and three groups of modified aramid insulating paper samples with different contents of SiO2 nanoparticles were produced. The modified aramid insulating paper samples with 0.5 wt%, 1 wt% and 3 wt% of nano-SiO2 particles were labeled as F1, F2 and F3, respectively, and the original aramid insulating paper samples without nano-SiO2 particles were labeled as F0.

Physical Test Results
The AFM test results of the modified aramid insulating paper samples are shown in Figure 2.
From Figure 2, it can be seen that the surface of the aramid insulating paper samples without nano-modification technology treatment is relatively flat, without obvious concave and convex changes, and overall surface roughness is low. With a gradual increase of nano-SiO2 particle filler, grooves on the surface of the modified aramid insulating paper samples increased significantly and surface roughness also increased, but when nano-SiO2 particle filler increased to a certain extent, grooves on the surface of the modified aramid insulating paper samples decreased again and surface roughness also decreased. The specific process is as follows: 1.
The SiO 2 nanoparticles with a diameter of about 12 nm are added to deionized water together with a 12% dispersant solution, and the mixed solution is dispersed ultrasonically with the aid of an ultrasonic cleaner; 2.
Take a fiber dissociator, add an appropriate amount of aramid short-cut fibers with a length of 12 mm to it, add the previously prepared nano-SiO 2 particles dispersion solution after being dispersed evenly, and stir it fully for 8 min. Subsequently add an appropriate amount of aramid pulp with a length of 2 mm, and stir it fully for 8 min to obtain a mixed original pulp containing nano-SiO 2 particles; 3.
Use a paper sample copying machine to copy and form the mixed original pulp containing the nano-SiO 2 particles to obtain dry paper containing the mixture of nano-SiO 2 particles; 4.
The obtained dry paper containing the mixture of nano-SiO 2 particles is superheated and calendered to obtain an improved aramid insulating paper containing nano-SiO 2 particles.

Test Methods
In this paper, a Dimension FastScan atomic force microscope (AFM) (Suzhou Feisman Precision Instruments Co., Suzhou, China) and Quanta FEG 250 scanning electron microscope (SEM) ( Carl Zeiss AG, Oberkochen, Germany)were used to test the external physical morphology of modified aramid insulating paper samples containing SiO 2 nanoparticles to investigate the effect of SiO 2 nanoparticles on the microscopic morphology of aramid insulating paper. Thermo SCIENTIFIC ESCALAB 250Xi X-ray photoelectron spectroscopy (XPS)(Thermo Fisher Scientific, Waltham, MA, USA) was used to analyze the surface chemistry of the modified aramid insulating paper samples to obtain the composition elements of the modified aramid insulating paper samples.
By controlling the amount of SiO 2 nanoparticles added during the preparation of the modified aramid insulating paper samples, the optimal modification ratio was investigated, and three groups of modified aramid insulating paper samples with different contents of SiO 2 nanoparticles were produced. The modified aramid insulating paper samples with 0.5 wt%, 1 wt% and 3 wt% of nano-SiO 2 particles were labeled as F1, F2 and F3, respectively, and the original aramid insulating paper samples without nano-SiO 2 particles were labeled as F0.

Physical Test Results
The AFM test results of the modified aramid insulating paper samples are shown in Figure 2.
From Figure 2, it can be seen that the surface of the aramid insulating paper samples without nano-modification technology treatment is relatively flat, without obvious concave and convex changes, and overall surface roughness is low. With a gradual increase of nano-SiO 2 particle filler, grooves on the surface of the modified aramid insulating paper samples increased significantly and surface roughness also increased, but when nano-SiO 2 particle filler increased to a certain extent, grooves on the surface of the modified aramid insulating paper samples decreased again and surface roughness also decreased. The results of the surface roughness Rq of the modified aramid paper samples are shown in Table 1. SEM tests were conducted on aramid insulating paper samples treated with the nanomodification technology, and test results are shown in Figure 3. From test results, it can be seen that with a gradual increase of nano-SiO2 particle filler, the surface of the aramid insulating paper sample is gradually becoming etched by nano-SiO2 particles. The reason is that with an increase of nano-SiO2 particle modifier, the degree of etching of nano-SiO2 particle modifier on the surface of aramid insulating paper sample gradually increases, the surface of the aramid insulating paper sample is gradually changed, the surface of the substrate on the surface of the sample is gradually peeled off resulting in gradual exposure of the nano-SiO2 particle modifier, and the quantity shows an increasing trend. The results of the surface roughness Rq of the modified aramid paper samples are shown in Table 1. SEM tests were conducted on aramid insulating paper samples treated with the nanomodification technology, and test results are shown in Figure 3. From test results, it can be seen that with a gradual increase of nano-SiO 2 particle filler, the surface of the aramid insulating paper sample is gradually becoming etched by nano-SiO 2 particles. The reason is that with an increase of nano-SiO 2 particle modifier, the degree of etching of nano-SiO 2 particle modifier on the surface of aramid insulating paper sample gradually increases, the surface of the aramid insulating paper sample is gradually changed, the surface of the substrate on the surface of the sample is gradually peeled off resulting in gradual exposure of the nano-SiO 2 particle modifier, and the quantity shows an increasing trend.  At the same time, as the number of SiO2 particle modifiers increased, the surface of the aramid paper sample was not further peeled off by the intermolecular force of the aggregation of SiO2 particle modifiers, which led to the reduction of surface roughness. The test results are shown in Figure 4.

Chemical Test Results
XPS tests were performed on aramid insulation paper samples treated with nanomodification technology, and the results are shown in Figure 5. Prepared modified aramid insulating paper samples were treated under a high vacuum environment for 6 h, and then were analyzed by XPS scanning on an ESCALABMKII XPS instrument, setting the vacuum level of the analysis chamber used for the experiment to 5 × 10 −6 Pa, setting the At the same time, as the number of SiO 2 particle modifiers increased, the surface of the aramid paper sample was not further peeled off by the intermolecular force of the aggregation of SiO 2 particle modifiers, which led to the reduction of surface roughness. The test results are shown in Figure 4. At the same time, as the number of SiO2 particle modifiers increased, the surface of the aramid paper sample was not further peeled off by the intermolecular force of the aggregation of SiO2 particle modifiers, which led to the reduction of surface roughness. The test results are shown in Figure 4.

Chemical Test Results
XPS tests were performed on aramid insulation paper samples treated with nanomodification technology, and the results are shown in Figure 5. Prepared modified aramid insulating paper samples were treated under a high vacuum environment for 6 h, and then were analyzed by XPS scanning on an ESCALABMKII XPS instrument, setting the vacuum level of the analysis chamber used for the experiment to 5 × 10 −6 Pa, setting the

Chemical Test Results
XPS tests were performed on aramid insulation paper samples treated with nanomodification technology, and the results are shown in Figure 5. Prepared modified aramid insulating paper samples were treated under a high vacuum environment for 6 h, and then were analyzed by XPS scanning on an ESCALABMKII XPS instrument, setting the vacuum level of the analysis chamber used for the experiment to 5 × 10 −6 Pa, setting the power to 500 W, using C1s as the energy counterpart for the whole XPS test process, and using MgKa rays as the energy source for the XPS test. The results of the whole XPS analysis are shown in Figure 4, where C1s, N1s, and O1s indicate the inner energy levels of C, N, and O. power to 500 W, using C1s as the energy counterpart for the whole XPS test process, and using MgKa rays as the energy source for the XPS test. The results of the whole XPS analysis are shown in Figure 4, where C1s, N1s, and O1s indicate the inner energy levels of C, N, and O.
(a) (b) It can be seen from test results that with a gradual increase of SiO2 nanoparticle filler, the characteristic peaks of C1s, N1s and O1s spectra changed, the chemical state of the surface of the aramid insulating paper samples treated by the nano-modification technology changed accordingly, and the contents of C, N and O all increased respectively. The content of the groups on the surface of the aramid fibers changed significantly, with the content of group -C-OH increasing and the content of the group -C=O showing an increasing and then decreasing trend, both of which were beneficial to improving the electrical properties of the aramid insulation paper samples, verifying the effectiveness of this modification method. It can be seen from test results that with a gradual increase of SiO 2 nanoparticle filler, the characteristic peaks of C1s, N1s and O1s spectra changed, the chemical state of the surface of the aramid insulating paper samples treated by the nano-modification technology changed accordingly, and the contents of C, N and O all increased respectively. The content of the groups on the surface of the aramid fibers changed significantly, with the content of group -C-OH increasing and the content of the group -C=O showing an increasing and then decreasing trend, both of which were beneficial to improving the electrical properties of the aramid insulation paper samples, verifying the effectiveness of this modification method.

Dielectric Constant Test Results
In order to understand the specific aramid modification results, a dielectric constant test was performed on this modified sample, and the test platform is shown in Figure 6. Four modified aramid samples were placed on the platform for testing. All dielectric constant test results are shown in Table 2.
Nanomaterials 2022, 12, x FOR PEER REVIEW 7 of

Dielectric Constant Test Results
In order to understand the specific aramid modification results, a dielectric constan test was performed on this modified sample, and the test platform is shown in Figure  Four modified aramid samples were placed on the platform for testing. All dielectric con stant test results are shown in Table 2.

Structure of the Oil-Filled Paper Condenser Type Bushing
Oil-filled paper condenser type bushings are one of the common types of high-vol age bushings. Although the bushing structure under different voltage levels is slightl different, the structure of the internal capacitor packet is basically the same [7]. This paper takes a 72.5 kV OIP bushing as the research object and uses COMSO multiphysics 5.2a simulation software to simulate the electric field distribution before an after replacement of the capacitor packet. An oil-filled paper condenser type bushing mainly composed of a connecting terminal, conservator, fastening spring, central guid rod, upper porcelain bushing, grounding flange, flange cylinder, lower porcelain bushin and capacitor packet. The bushing is filled with insulating oil, which together with th capacitor packet composed of insulating paper and aluminum foil forms the oil-impreg nated paper insulation system. The length of the bushing is 2245 mm, the length of th upper porcelain bushing is 850 mm, the length of the lower porcelain bushing is 210 mm the lead-in length of the cable is 1865 mm, the inner diameter of the conductive rod is 3 mm, and the maximum radius of the capacitor packet is 43.6 mm. The structure of th 72.5 kV OIP bushing is shown in Figure 7.  Oil-filled paper condenser type bushings are one of the common types of high-voltage bushings. Although the bushing structure under different voltage levels is slightly different, the structure of the internal capacitor packet is basically the same [7]. This paper takes a 72.5 kV OIP bushing as the research object and uses COMSOL multiphysics 5.2a simulation software to simulate the electric field distribution before and after replacement of the capacitor packet. An oil-filled paper condenser type bushing is mainly composed of a connecting terminal, conservator, fastening spring, central guide rod, upper porcelain bushing, grounding flange, flange cylinder, lower porcelain bushing and capacitor packet. The bushing is filled with insulating oil, which together with the capacitor packet composed of insulating paper and aluminum foil forms the oil-impregnated paper insulation system. The length of the bushing is 2245 mm, the length of the upper porcelain bushing is 850 mm, the length of the lower porcelain bushing is 210 mm, the lead-in length of the cable is 1865 mm, the inner diameter of the conductive rod is 35 mm, and the maximum radius of the capacitor packet is 43.6 mm. The structure of this 72.5 kV OIP bushing is shown in Figure 7. In actual operation, the central guide rod of the bushing carries current and is fixed to the top of the transformer through a grounding flange. The capacitor packet uniformizes the electric field distribution between the central guide rod and the ground potential through an internal insulation system composed of insulating oil, insulating paper and aluminum foil. The thin capacitor packet bears the transformer phase voltage, the electric field strength is large, the operating environment is harsh, and manufacturing process requirements are high. Due to the existence of air bubbles in the capacitor packet, poor impregnation, and air gaps between oil and paper and other process problems, bushing failures and even transformer explosion accidents occur from time to time [25][26][27].

Insulation Principle of Capacitor Packet
The oil-impregnated paper insulation system composed of insulating oil, insulating paper and aluminum foil is the main insulation for the transformer bushing. By alternately winding 0.08 mm-0.12 mm-thick insulating paper and 0.01 mm-thick aluminum foil on the central guide rod to form a plurality of capacitors with the same capacitance value, uniform voltage and electric fields are achieved [28]. Generally, the number of layers of aluminum foil (also known as capacitor plate) in the 72.5 kV OIP bushing is 18, the outermost layer of which is the end screen. The end screen is grounded during operation, which is at the ground potential. The voltage between the center guide rod and the end screen is the sum of the voltages between the plates.
If the error caused by the edge effect of the capacitor plate is ignored, it can be considered the electric flux passing through each level of the plate does not change: where D is the electric flux, S is the area of the electrode plate, ε0 is the dielectric constant in vacuum, εr is the relative permittivity of the material, Er is the radial electric field strength here, r is the radius of the electrode, l is the length of the electrode, and c1 is a constant. Structure of the bushing capacitor packet is shown in Figure 8. Suppose the voltage between two adjacent layers of plates is dU, and the axial field strength is El, there are:  In actual operation, the central guide rod of the bushing carries current and is fixed to the top of the transformer through a grounding flange. The capacitor packet uniformizes the electric field distribution between the central guide rod and the ground potential through an internal insulation system composed of insulating oil, insulating paper and aluminum foil. The thin capacitor packet bears the transformer phase voltage, the electric field strength is large, the operating environment is harsh, and manufacturing process requirements are high. Due to the existence of air bubbles in the capacitor packet, poor impregnation, and air gaps between oil and paper and other process problems, bushing failures and even transformer explosion accidents occur from time to time [25][26][27].

Insulation Principle of Capacitor Packet
The oil-impregnated paper insulation system composed of insulating oil, insulating paper and aluminum foil is the main insulation for the transformer bushing. By alternately winding 0.08 mm-0.12 mm-thick insulating paper and 0.01 mm-thick aluminum foil on the central guide rod to form a plurality of capacitors with the same capacitance value, uniform voltage and electric fields are achieved [28]. Generally, the number of layers of aluminum foil (also known as capacitor plate) in the 72.5 kV OIP bushing is 18, the outermost layer of which is the end screen. The end screen is grounded during operation, which is at the ground potential. The voltage between the center guide rod and the end screen is the sum of the voltages between the plates.
If the error caused by the edge effect of the capacitor plate is ignored, it can be considered the electric flux passing through each level of the plate does not change: where D is the electric flux, S is the area of the electrode plate, ε 0 is the dielectric constant in vacuum, ε r is the relative permittivity of the material, E r is the radial electric field strength here, r is the radius of the electrode, l is the length of the electrode, and c 1 is a constant. Structure of the bushing capacitor packet is shown in Figure 8. Suppose the voltage between two adjacent layers of plates is dU, and the axial field strength is E l , there are:  In actual operation, the central guide rod of the bushing carries current and is fixed to the top of the transformer through a grounding flange. The capacitor packet uniformizes the electric field distribution between the central guide rod and the ground potential through an internal insulation system composed of insulating oil, insulating paper and aluminum foil. The thin capacitor packet bears the transformer phase voltage, the electric field strength is large, the operating environment is harsh, and manufacturing process requirements are high. Due to the existence of air bubbles in the capacitor packet, poor impregnation, and air gaps between oil and paper and other process problems, bushing failures and even transformer explosion accidents occur from time to time [25][26][27].

Insulation Principle of Capacitor Packet
The oil-impregnated paper insulation system composed of insulating oil, insulating paper and aluminum foil is the main insulation for the transformer bushing. By alternately winding 0.08 mm-0.12 mm-thick insulating paper and 0.01 mm-thick aluminum foil on the central guide rod to form a plurality of capacitors with the same capacitance value, uniform voltage and electric fields are achieved [28]. Generally, the number of layers of aluminum foil (also known as capacitor plate) in the 72.5 kV OIP bushing is 18, the outermost layer of which is the end screen. The end screen is grounded during operation, which is at the ground potential. The voltage between the center guide rod and the end screen is the sum of the voltages between the plates.
If the error caused by the edge effect of the capacitor plate is ignored, it can be considered the electric flux passing through each level of the plate does not change: where D is the electric flux, S is the area of the electrode plate, ε0 is the dielectric constant in vacuum, εr is the relative permittivity of the material, Er is the radial electric field strength here, r is the radius of the electrode, l is the length of the electrode, and c1 is a constant. Structure of the bushing capacitor packet is shown in Figure 8. Suppose the voltage between two adjacent layers of plates is dU, and the axial field strength is El, there are:  Assuming the axial field strength as a constant, it can be obtained from Equation (1) that: In the formula, c 2 and c 3 are constants. Integrate formula (4) to obtain: In the formula, c 4 is a constant. Substituting formula (3) and boundary conditions r = r 0 , l = l 0 ; r = r n , l = l n into formula (5), it can be solved as follows: From formula (7), it can be seen when E l is a constant, the product of E r and rl is inversely proportional. The radial field strength of the capacitor packet presents a U-shaped distribution, the electric field strength is higher near the end screen and the center guide rod, and the electric field strength in the middle part is lower. Through the equal-capacitance design method, the length of each plate in the packet is solved by using the fact that the capacitance of each plate is equal. Suppose the interlayer capacitance is C k , then: 2πε 0 ε r l 1 ln(r 1 /r 0 ) = 2πε 0 ε r l 2 ln(r 2 /r 1 ) = · · · = 2πε 0 ε r l n ln(r n /r n−1 ) (8) It can be seen from formulas (7) and (8) that replacing the insulating material of the capacitor packet causes a change in the dielectric constant, resulting in a change in the electric field strength.

Finite Element Simulation Theory
The finite element method is a numerical calculation method based on the difference method and the variational principle. It was used in the structural mechanics in early days and was then gradually applied in the fields of electromagnetic field and fluid mechanics [8].
The electric field of the high-voltage bushing is generated by the industrial-frequency alternating current passing through the center guide rod. The frequency is 50 Hz, which can be reduced to a two-dimensional axisymmetric electrostatic field problem. The Laplace equation of the potential is: The boundary conditions are the first and second types of boundary conditions, which are: Solving the electric field of the transformer bushing in the COMSOL finite element software can be divided into six steps, including selecting the electric field module according to the physical field of the problem to be solved; building the model according to the structure; setting the material parameters of each part; setting the load and boundary conditions; performing meshing according to the characteristics of the electric field and selecting a solver to solve the model [29].

Parameters and Settings of Finite Element Simulation
Due to highly axisymmetric characteristics of the transformer bushing, the threedimensional model can be simplified into a two-dimensional axisymmetric model for calculation, and this has little effect on calculation accuracy. According to the structure of the transformer bushing, a two-dimensional symmetrical model is established in COMSOL, as shown in Figure 9.
Solving the electric field of the transformer bushing in the COMSOL finite software can be divided into six steps, including selecting the electric field module ing to the physical field of the problem to be solved; building the model accordin structure; setting the material parameters of each part; setting the load and bound ditions; performing meshing according to the characteristics of the electric field lecting a solver to solve the model [29].

Parameters and Settings of Finite Element Simulation
Due to highly axisymmetric characteristics of the transformer bushing, the t mensional model can be simplified into a two-dimensional axisymmetric model culation, and this has little effect on calculation accuracy. According to the structu transformer bushing, a two-dimensional symmetrical model is established in CO as shown in Figure 9. We add material properties to the solution domain according to bushing st In the simulation calculation, the relative permittivity of each material is set as sh Table 3: Table 3. Relative permittivity of each material.

Material
Insulating Oil Oil-Impregnated Paper According to actual working conditions of the transformer bushing, a worki age of 41.9 kV is applied to one end of the center guide rod, and a working curren A is applied to the other end. The pole plate in the capacitor packet is set to a potential, and the end screen is connected to the flange by a grounding lead, whi to the ground potential in the model, and the flange is also set to the ground poten air area at the length of 1.5 m and a fuel tank area at the length of 0.5 m are set out bushing to simulate actual working conditions.
Since the thickness of the aluminum foil is only 0.01 mm, the mesh size wi tremely small, the number of meshes will be greatly increased, the quality of th will be low, and the calculation accuracy will be affected when the finite element is used for meshing. Setting it as a single-layer material in the simulation can avoi ing problems.

Analysis of the Simulation Results
The simulation yields a cloud diagram of electric field distribution of the tran bushing as shown in Figure 10. A cut-off line is taken at the middle of the bushin itor packet (near the flange) and the voltage distribution along the line is obta shown in Figure 11. We add material properties to the solution domain according to bushing structure. In the simulation calculation, the relative permittivity of each material is set as shown in Table 3: Table 3. Relative permittivity of each material.

Material
Insulating Oil

Oil-Impregnated Paper
Oil-Impregnated Aramid Paper According to actual working conditions of the transformer bushing, a working voltage of 41.9 kV is applied to one end of the center guide rod, and a working current of 630 A is applied to the other end. The pole plate in the capacitor packet is set to a floating potential, and the end screen is connected to the flange by a grounding lead, which is set to the ground potential in the model, and the flange is also set to the ground potential. An air area at the length of 1.5 m and a fuel tank area at the length of 0.5 m are set outside the bushing to simulate actual working conditions.
Since the thickness of the aluminum foil is only 0.01 mm, the mesh size will be extremely small, the number of meshes will be greatly increased, the quality of the mesh will be low, and the calculation accuracy will be affected when the finite element method is used for meshing. Setting it as a single-layer material in the simulation can avoid meshing problems.

Analysis of the Simulation Results
The simulation yields a cloud diagram of electric field distribution of the transformer bushing as shown in Figure 10. A cut-off line is taken at the middle of the bushing capacitor packet (near the flange) and the voltage distribution along the line is obtained as shown in Figure 11.
The radial electric field distribution of 72.5 kV bushing can be calculated using Equation (7). As shown in Figure 12, it can be seen that when the bushing is operated at the rated phase voltage, the electric field intensity in the bushing capacitor packet shows a "U"-shaped distribution, with the highest electric field intensity at the outermost end screen of the capacitor packet and the lowest field intensity at the middle part. At the same time, Figure 12 also shows the simulation of the electric field intensity; you can see the simulation value and the theoretical value are quite close. The maximum value of the electric field intensity through simulation is 2.68 kV/mm, and the electric field intensity distribution of the capacitor packet is quite uniform, which is also a curve showing a "U" type distribution, the same as the theoretical value. The difference between the simulated value and the theoretical value at the innermost side of the capacitor packet is large, probably due to the large difference between the parameters of the metal guide rod and the insulation paper, resulting in a large calculation error.  The radial electric field distribution of 72.5 kV bushing can be calculated using Equa tion (7). As shown in Figure 12, it can be seen that when the bushing is operated at th rated phase voltage, the electric field intensity in the bushing capacitor packet shows "U"-shaped distribution, with the highest electric field intensity at the outermost en screen of the capacitor packet and the lowest field intensity at the middle part. At the sam time, Figure 12 also shows the simulation of the electric field intensity; you can see th simulation value and the theoretical value are quite close. The maximum value of the ele tric field intensity through simulation is 2.68 kV/mm, and the electric field intensity di tribution of the capacitor packet is quite uniform, which is also a curve showing a "U  The radial electric field distribution of 72.5 kV bushing can be calculated using Equation (7). As shown in Figure 12, it can be seen that when the bushing is operated at the rated phase voltage, the electric field intensity in the bushing capacitor packet shows a "U"-shaped distribution, with the highest electric field intensity at the outermost end screen of the capacitor packet and the lowest field intensity at the middle part. At the same time, Figure 12 also shows the simulation of the electric field intensity; you can see the simulation value and the theoretical value are quite close. The maximum value of the electric field intensity through simulation is 2.68 kV/mm, and the electric field intensity distribution of the capacitor packet is quite uniform, which is also a curve showing a "U" type distribution, the same as the theoretical value. The difference between the simulated value and the theoretical value at the innermost side of the capacitor packet is large, probably due to the large difference between the parameters of the metal guide rod and the insulation paper, resulting in a large calculation error. To investigate the effect of the relative permittivity of the capacitor packet on the electric field intensity of the OIP bushing, a parametric scan is performed using COMSOL simulation software. A variety of electric field distributions were obtained by setting different dielectric constants. By comparing each electric field distribution, the dielectric constant for optimal electric field intensity distribution can be found. Set the dielectric constant variation interval as (2,10) and the variation step as 0.2. The simulation solves the electric field intensity at each value, and takes the electric field intensity of the same position of the capacitor packet for comparison, as shown in Figure 13. From Figure 13, it can be seen the electric field intensity exhibits a rotation around a point in the center as the dielectric constant changes. In the inner part of the capacitor packet, the electric field intensity shows a decrease with increasing dielectric constant; near the outer part of the capacitor packet, the electric field intensity increases with increasing dielectric constant.  To investigate the effect of the relative permittivity of the capacitor packet on the electric field intensity of the OIP bushing, a parametric scan is performed using COMSOL simulation software. A variety of electric field distributions were obtained by setting different dielectric constants. By comparing each electric field distribution, the dielectric constant for optimal electric field intensity distribution can be found. Set the dielectric constant variation interval as (2,10) and the variation step as 0.2. The simulation solves the electric field intensity at each value, and takes the electric field intensity of the same position of the capacitor packet for comparison, as shown in Figure 13. From Figure 13, it can be seen the electric field intensity exhibits a rotation around a point in the center as the dielectric constant changes. In the inner part of the capacitor packet, the electric field intensity shows a decrease with increasing dielectric constant; near the outer part of the capacitor packet, the electric field intensity increases with increasing dielectric constant. value and the theoretical value at the innermost side of the capacitor packet is large, probably due to the large difference between the parameters of the metal guide rod and the insulation paper, resulting in a large calculation error. To investigate the effect of the relative permittivity of the capacitor packet on the electric field intensity of the OIP bushing, a parametric scan is performed using COMSOL simulation software. A variety of electric field distributions were obtained by setting different dielectric constants. By comparing each electric field distribution, the dielectric constant for optimal electric field intensity distribution can be found. Set the dielectric constant variation interval as (2,10) and the variation step as 0.2. The simulation solves the electric field intensity at each value, and takes the electric field intensity of the same position of the capacitor packet for comparison, as shown in Figure 13. From Figure 13, it can be seen the electric field intensity exhibits a rotation around a point in the center as the dielectric constant changes. In the inner part of the capacitor packet, the electric field intensity shows a decrease with increasing dielectric constant; near the outer part of the capacitor packet, the electric field intensity increases with increasing dielectric constant.  The maximum and variance of the electric field intensity in the capacitor packet under each electric field distribution are calculated, as shown in Figure 14. High local electric field intensity will cause partial discharge, which will lead to insulation degradation and even insulation breakdown, so maximum electric field intensity in the capacitor packet needs to be controlled. The capacitor packet electric field intensity variance is used to characterize uniformity of the electric field intensity inside the capacitor packet, and a uniform electric field environment is beneficial to the long-term work of the insulation material.
needs to be controlled. The capacitor packet electric field intensity variance is used to characterize uniformity of the electric field intensity inside the capacitor packet, and a uniform electric field environment is beneficial to the long-term work of the insulation material.
As shown in Figure 14, when relative permittivity is small, the electric field is uniform and the maximum electric field intensity is small, and where the relative permittivity is 2.8, the electric field intensity variance is smallest and the electric field is most uniformly distributed. However, fabrication of the material does not guarantee exact control of the dielectric constant to 2.8. Controlling the relative dielectric constant between 2.4 and 3.4 can achieve a more uniform electric field and a smaller maximum electric field intensity.

Figure 14.
Maximum and variance of the electric field intensity.

Conclusions
In this paper, nano-modification technology was used to modify aramid insulation paper, and the physical and chemical properties of the aramid insulation paper before and after modification were tested and analyzed. Based on the technical parameters of the modified aramid insulation paper, a model of OIP bushing was constructed, and its effect on the electric field distribution was studied. The main conclusions of this paper are as follows: 1. The surface of aramid insulation paper sample can be effectively stripped by nanomodification technology, which increases the degree of etching of a nano-SiO2 particle modifier on its surface gradually and the surface roughness of the sample, making the surface gouge situation increase significantly. However, with the increase of nano-SiO2 particle modifier content, it is influenced by the aggregation of nano-particles into clusters, which leads to the reduction of its surface roughness instead. 2. After comparing the samples of aramid insulation paper modified with nano-SiO2 particles with those of aramid insulating paper without nano-SiO2 particles, it can be seen that with the gradual increase of nano-SiO2 particles filler, the chemical state of its surface changed accordingly, the content of groups on the surface of aramid fiber changed significantly, the content of group-C-OH increased, and the content of group-C=O showed a trend of first increasing and then decreasing. As shown in Figure 14, when relative permittivity is small, the electric field is uniform and the maximum electric field intensity is small, and where the relative permittivity is 2.8, the electric field intensity variance is smallest and the electric field is most uniformly distributed. However, fabrication of the material does not guarantee exact control of the dielectric constant to 2.8. Controlling the relative dielectric constant between 2.4 and 3.4 can achieve a more uniform electric field and a smaller maximum electric field intensity.

Conclusions
In this paper, nano-modification technology was used to modify aramid insulation paper, and the physical and chemical properties of the aramid insulation paper before and after modification were tested and analyzed. Based on the technical parameters of the modified aramid insulation paper, a model of OIP bushing was constructed, and its effect on the electric field distribution was studied. The main conclusions of this paper are as follows: 1.
The surface of aramid insulation paper sample can be effectively stripped by nanomodification technology, which increases the degree of etching of a nano-SiO 2 particle modifier on its surface gradually and the surface roughness of the sample, making the surface gouge situation increase significantly. However, with the increase of nano-SiO 2 particle modifier content, it is influenced by the aggregation of nano-particles into clusters, which leads to the reduction of its surface roughness instead.

2.
After comparing the samples of aramid insulation paper modified with nano-SiO 2 particles with those of aramid insulating paper without nano-SiO 2 particles, it can be seen that with the gradual increase of nano-SiO 2 particles filler, the chemical state of its surface changed accordingly, the content of groups on the surface of aramid fiber changed significantly, the content of group-C-OH increased, and the content of group-C=O showed a trend of first increasing and then decreasing.

3.
Electric intensity has the smallest variance and the electric field distribution is the most uniform. However, the manufacturing process of the material does not guarantee that the dielectric constant is accurately controlled at 2.8. Controlling the relative dielectric constant between 2.4 and 3.4 can achieve a more uniform electric field and a smaller maximum intensity. In this paper, the modified aramid paper insulation dielectric constant is basically close to the range of 2.4-3.4, which can achieve a more uniform electric field and a smaller maximum electric intensity. Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available in Insulation Performance and Simulation Analysis of Nano-modified Aramid Applied in High-Voltage Bushing.