Abstract

In order to consider the influence of nonlinear characteristics of porcelain pillar electrical equipment on the dynamic response under seismic excitation, a theoretical analysis method of nonlinear dynamics was raised to define the nonlinear parameter of the flange connection and establish a dynamic model of porcelain pillar electrical equipment. The theoretical analysis and the test results have a good degree of fitting, which verifies the correctness of the dynamic model and reveals the nonlinear seismic response law of the porcelain pillar equipment. According to the results, both the nonlinear calculation results of the displacement at the top end of the porcelain pillar and the stress at the bottom end of the porcelain pillar are smaller than the results in the linear calculation. The difference between them increases gradually with the increase of seismic excitation. The differences in displacement and stress at a seismic excitation of 0.5 g are 14.58% and 23.25%, respectively. When the nonlinear parameter increases to a certain value, the impact of maximum stress on the bottom of the porcelain pillar is very small and the change is not obvious. The research provides a theoretical reference for the seismic design of pillar electrical equipment.

1. Introduction

The electric power system is an important part of the lifeline project, which is related to the normal life of the people and the development of industry and manufacturing. The earthquakes have caused serious damage to the electrical equipment. An electrical substation, as the node of an electrical grid, is extremely important for the reliable service of an electrical grid system [1]. A wide variety of porcelain pillar electrical equipment, such as post insulators and surge arresters, is a major component of an electrical substation [2]. Due to the brittle nature of porcelain, equipment made of porcelain is vulnerable to earthquake shaking. For porcelain pillar electrical equipment when it is subjected to earthquake, due to the small rigidity of the cemented part of the flange connection, it is prone to damage at the flange connection. The post-earthquake field investigation and shaking table tests (Figures 1 and 2) showed failure characteristics of fracture during the earthquake damage [35].

Many scholars have conducted extensive researches on the seismic response mechanism of porcelain pillars. Chinese codes are given in [69] combined with experimental research. In the recommended practice in the relevant code, the calculation equation for the bending stiffness of the joint between the flange and the porcelain bushing isin which Kc is the bending stiffness of the connection between the flange and the porcelain bushing (in N·m/rad); dc is the porcelain bushing outer diameter (in m); hc is the height of porcelain bushing at cemented part (in m); te is the thickness of the cemented part between the flange and the porcelain bushing (in m) (refer to Figure 3(a)); and β is the calculation parameter of the bending stiffness of the joint between the porcelain bushing and the flange, which is recommended as β = 6.54 × 107 in the Chinese code. The physical concept of the above equation is relatively clear, which basically reflects the mechanical characteristics of the cement layer between the flange and the porcelain bushing. For the bolted connection, the bolts are usually pretensioned and there is no opening between the two faces when the equipment is subjected to bending moment. The stiffness of bolted connection can be assumed infinitely large. Therefore, the stiffness between metal cap and porcelain pillar is weaker and the bending capacity is determined by stiffness between metal cap and porcelain pillar. However, its accuracy depends on the empirical parameter of the bending stiffness of the porcelain bushing and the flange.

Zhang et al. [10] obtained the bending rigidity calculation parameter for the connections between UHV porcelain bushing and flange by conducting the bending rigidity test and, accordingly, applied this calculation parameter to the finite element calculation of the simulated UHV porcelain bushings. Singh et al. [11] introduced the response spectrum analysis to evaluate linear control systems for seismic inputs defined by code-prescribed or site-specific ground response spectra. Gilani et al. [12] carried out a shaking table test for two 230 kV, 3000 A, Type-U transformer bushings and found that both bushings survived earthquake shaking compatible with the IEEE spectrum for high-level qualification when mounted on a rigid frame. When mounted on a flexible frame, one 230 kV bushing was unable to sustain the high-level qualification shaking without oil leakage and a slip of the porcelain units. Filiatrault and Matt [13] investigated numerically the dynamic response of porcelain bushings mounted on transformer tanks. The results of the numerical study show that large amplification occurs when the fundamental frequency of the porcelain bushing is tuned with the fundamental frequency of the transformer tank. Bender and Farid [14] added additional evidence that bushing amplifications greater than 2.0 are possible, and they demonstrated that the complexity of these systems makes generalizations of design. Villaverde et al. [15] undertook the study to assess the adequacy of the amplification factor of 2.0 specified by the Institute of Electrical and Electronics Engineers in Standard IEEE 693-1997 for the seismic qualification of transformer bushings. Schiff [16] proposed a new method to seismically qualify power transformers and their bushings to meet the requirements of IEEE Standard 693. He et al. [17] analyzed the influence of different types of connections on the seismic responses of the oil-side bushings and a simplified method to rapidly estimate the seismic responses of transformer bushings was implemented.

Cao [18] conducted the characteristic research on the recovery capacity of pillar insulators and discussed the impact of various factors such as the insulator cementing parts between insulator flange and porcelain body on the seismic performance of insulator pillars. Schiff [19] presented a method for evaluating the reliability of ceramic structural members subjected to earthquake-induced vibrations which used the Weibull distribution to represent the probability of failure of the brittle ceramic material. The authors of [2022] introduce the experimental study of seismic performance of pillar electrical equipment of different voltage classes and different types. In order to improve the seismic capacity of pillar electrical equipment, composite materials are often used in the porcelain pillar equipment. IEC 61462 [23] and IEEE Std 693-2005 [24] introduce the research on the improvement of bearing capacity of porcelain pillar equipment by using new composite materials. Mohammadi et al. [25] proposed a Four Degrees-of-Freedom (4-DOF) system in order to assess the dynamic properties of substation support structures and verified through the finite element method (FEM).

In the above researches, during the seismic analysis of porcelain pillar electrical equipment, the beams with equivalent rigidity are often used directly to represent the flange connections. This paper focuses on the nonlinear factors at the flange connections by introducing a nonlinear rigidity parameter. The nonlinear dynamic model is then established. After that, the seismic performance of multiple porcelain pillars is studied and compared with the shaking table test results, which proved the correctness of the nonlinear dynamic model. Finally, the maximum stresses at the bottom of porcelain pillars are analyzed when the nonlinear rigidity parameter changes. The maximum stress at the bottom of the porcelain pillar with a different seismic excitation can be predicted through this way.

2. Nonlinear Dynamic Model of Porcelain Pillar Equipment

In an electrical substation, some porcelain pillar equipment is directly installed on the foundations. Take a piece of typical flanged 4-section porcelain post insulator as an example, as illustrated in Figure 4. Figure 4 also shows the key structural and typical sections of a solid core post insulator with sheds. The equipment has structural characteristics (1): the main body of the equipment is made of porcelain core as annotated in Figure 4; (2) the porcelain component is connected through metal caps at both ends with bolted flange connections. The connections between porcelain pillars contain two parts, cemented connections between metal caps and flange, and connection between two flanges with bolts.

Figure 5 is an analysis diagram of the porcelain pillar equipment structure, where xb is the seismic excitation, and the sequential numbers of porcelain pillars are 1, 2, 3, and 4, respectively, from bottom to top. The equipment is mainly in a bending state when it is subjected to seismic and the flange connection between two pieces of porcelain pillars rotates by a certain angle θ. The bottom end of the upper porcelain pillar rotates by θi, and the upper end of the lower porcelain rotates by θi − 1. So, the relative angle θ can be expressed as θ = θi − θi − 1. The angles of each flange connection are θ1, θ2, θ3, and θ4, respectively, as shown in Figure 5.

The analysis diagram of flange connections is illustrated in Figure 3. Figure 3(a) shows the deformation of the cemented part. Figure 3(b) shows the mechanical simplified diagram of the flange connection. According to the deformation characteristics, the connection between porcelain and metal cap can be regarded as a spring with rigidity kid. Correspondingly, the connection between lower porcelain and flange can be regarded as a spring with rigidity k(i − 1)u. The cemented part between metal caps and flange can be equivalent to a series spring with the equivalent rigidity ki, and ki, can be expressed as follows:where kid and k(i − 1)u are calculated as follows [2]:in which du and dd are the porcelain diameters at cemented parts of (i − 1)th and ith porcelain pillar, respectively; hu and hd are the cemented heights of (i − 1)th and ith porcelain pillar, respectively; tu and td are the thicknesses of cement layer of (i − 1)th and ith porcelain pillar, respectively; βu and βd are bending rigidity parameters at the connections of (i − 1)th and ith porcelain pillar, respectively, the value of β is 6.54 × 107 when the porcelain diameter at the cemented part is smaller than 275 mm, the value of β is 5.0 × 107 when the porcelain diameter at the cemented part is larger than 375 mm, and the value of β can be obtained by linear interpolation when the porcelain diameter at the cemented part is between 275 mm to 375 mm.

As can be seen from the above analysis, the rotational bending moment caused by the relative rotation of the up and down flanges can be expressed using

Equation (4) expresses rotational bending moment under ideal conditions. The exact mechanics of the rotation stiffness ki is complex, which involves nonlinear mechanical properties of the cemented material and the mechanics at the contact interfaces. These factors should be included in Mi to make it more realistic. In most previous studies, the flange joint was equivalent to a beam with equivalent rotational stiffness [2, 26, 27]. Figure 6 shows a diagram of the force of the microelement. According to the elastic theory, when the beam subjected to large deformation, the strain can be expressed asin which u represents the horizontal displacement, θ represents the rotational angle, and represents the curvature.

Then, the stress can be calculated usingin which E represents Young’s modulus.

Further, the bending moment on the cross section of the porcelain bushing isin which A represents the area of cross section.

In addition, In equation (8), I represents the moment of inertia.

Consequently, equation (7) is reduced to

θ can be solved by methods of separating variable; the expression of θ is shown in where represents the vibration shape function of porcelain pillar, and is the rotational angle of flange connection.

When the height of the cemented part is h, the bending moment provided by the flange is

Let,; then,

Referring to equation (12), equation (4) is rewritten as follows after considering the nonlinear factors:

In equation (13), the parameter η is the nonlinear rigidity parameter combined with discontinuous deformation at the flange when subjected to large deformation. It is closely related to the vibration shape function of porcelain pillar assembly. It is difficult to get through a theoretical solution, so the value of η is based on an experimental solution. When η = 0, equation (13) is used to calculate the rotational bending moment under the ideal conditions represented by equation (4). The different values of η correspond to different degrees of the nonlinearity of the flange cemented part.

If there are Nc porcelain pillars, there will be Nc flange connections. The elastic potential energy at a flange connection is calculated as follows:

Since the rotation of porcelain pillars causes the change of gravitational potential energy of porcelain pillars, the total gravitational potential energy is calculated as follows:in which Lj represents the length of jth porcelain pillar.

The kinetic energy of porcelain pillars consists of average kinetic energy and rotational kinetic energy:

in which is the centroid velocity of the ith porcelain pillar, and mi is the mass of the ith porcelain pillar. For the four-section porcelain pillar system, the speed relationship of the ith porcelain pillar is illustrated in Figure 7:

As illustrated in Figure 7, the velocity of 1st porcelain pillar bottom is

The centroid velocity of the 1st porcelain pillar is calculated as follows:

The centroid velocity of the 2nd porcelain pillar is calculated as follows:

The centroid velocity of the 3rd porcelain pillar is calculated as follows:

The centroid velocity of the 4th porcelain pillar is calculated as follows:

The Hamilton variational equation [28] is as follows:

In equation (21), U = Uk + . Then,

After equations (23) and (24) are substituted into equation (22) and further simplification, equation (22) is rewritten in matrix form as follows:where

In addition,

3. Dynamic Response of Porcelain Pillar Electrical Equipment under Seismic Excitation

3.1. Shaking Table Test

In this section, the accuracy and effectiveness of the proposed modeling method will be validated and demonstrated by comparing the analysis results with the experimental results.

The shaking table test set-up of a 1000 kV porcelain surge arrester is shown in Figure 8(a). This equipment contains 4 insulators with a similar configuration. The mass, dimension, and material parameters are listed in Table 1.

In the shaking table tests, the strain gauges were placed in the X direction of the root of each section of the casing (two on the opposite side); the accelerometer is placed in the X direction of the top of the vibrating table and the top of each section of the casing, which is shown in Figure 8(b).

The amplitudes of the input artificial waves are 0.1 g, 0.15 g, 0.2 g, and 0.3 g, and the direction of the input artificial wave is a single X direction (refer to Figure 8(b)). Figure 9 shows the input seismic wave of the artificial ground motions, a broad resonant period range from 0.1 s to 0.9 s, and the stress responses at the bottom with peak ground acceleration of 0.15 g. The vibration table test obtains the strain at the bottom end of the porcelain pillar. The stress can be obtained by multiplying the elastic modulus of the porcelain bushing (110 GPa) by the strain. Figure 9(b) shows the maximum stress at the bottom of porcelain.

As seen in Figure 9, the maximum acceleration of the ground is about 2.3 m/s2 and the maximum stress at the bottom end of the porcelain pillar is 11.14 MPa. Figure 10 shows the results when the peak value of the input wave is 0.3 g. As seen in Figure 10, the maximum acceleration of the ground is about 4 m/s2 and the maximum stress at the bottom end of the porcelain pillar is 22.02 MPa.

3.2. Seismic Wave Input in Theoretical Analysis

The seismic wave data are those for artificial waves jointly proposed by China Electric Power Research Institute and China Earthquake Disaster Prevention Center. This artificial wave is fitted by the seismic acceleration response spectrum with a characteristic cycle of 0.9 s that can envelop the characteristic cycles of Category I∼Category III soil sites, and it is suitable for seismic design, test, and seismic performance evaluation of HV and UHV electrical equipment. The artificial wave time history and the fitting data for artificial wave response spectrum and target spectrum at a peak seismic acceleration of 1g are illustrated in Figure 11.

The data Xb can be fitted with a periodic function :

in which

When the value of NF is 500, the fitting result is illustrated in Figure 11. As seen from Figure 11, the blue fitting curve has a high matching degree with the black seismic wave curve.

3.3. Method Validation

In order to verify the correctness of the theoretical model, the numerical method was used to solve equation (25) with the fitting seismic wave proposed in Section 3.2, and the results were compared with the curve obtained by the shaking table test in Section 3.1.

Based on the calculation by the kinetic model established according to equation (25), the first-order frequency of the four-section porcelain pillar equipment is 2.28 Hz, and the result obtained by shaking table test is 2.18 Hz, so there is a difference of 4.6% between the test value and the theoretical value. It indicates that the theoretical analysis is correct.

Figure 12 is a comparison curve of stress time histories at the bottom end of the porcelain pillars obtained by theoretical analysis and vibration table test at a = 0.15 g, a = 0.3 g respectively, when the value of nonlinear rigidity parameter η is taken as 20,000. The maximum stresses at the bottom of the 1st section of the equipment obtained by test and theoretical analysis are 11.14 MPa and 11.00 MPa, respectively, at a = 0.15 g, indicating that the difference between maximum stress values is 1.2%. The theoretical and test curves also show a similar change trend.

The maximum stresses at the bottom of the 1st section of the equipment obtained by test and theoretical analysis are 22.02 MPa and 20.30 MPa, respectively, at a = 0.3 g, indicating that the difference between maximum stress values is 9.04%. Seismic excitation usually has nonlinear characteristics. With the increase of seismic excitation, the nonlinear characteristics of the porcelain pillar become more significant, which leads to the fact that the difference between test and analytical value becomes larger.

Figure 13 shows a comparative diagram of maximum stresses at the bottom of the 2nd section pillar. When the peak acceleration of the input seismic wave is 0.15 g, the maximum value of the nonlinear theoretical result is 9.41 MPa and the maximum value of the test is 9.6 MPa. The difference between them is 2.0%. When the peak acceleration of the input seismic wave is 0.3 g, the maximum value of the nonlinear theoretical result is 19.25 MPa and the maximum value of the test is 21.00 MPa. The difference between them is 8.3%.

Through the above analysis, the theoretical calculation results have a rather high goodness of fit with the test results, thus verifying the correctness of the analytical model. When the value of η is changed, new nonlinear theoretical analysis results are obtained. The differences between nonlinear analysis results and test results increase. According to the comparison with the test results, it can be assumed that for the four-section porcelain pillar equipment, η = 20,000 can reflect the nonlinear factors such as materials of cemented parts and large deformation in the actual applications.

4. Seismic Response of 1,000 kV Porcelain Pillar Arrester Equipment

Equation (25) is solved to obtain the linear and nonlinear dynamic response of the 1,000 kV arrester at η = 0 and η = 20,000, respectively. Figures 1419 illustrate the curves of the displacement at the top end of the porcelain pillar and the maximum stress at the bottom end of the porcelain pillar under different seismic excitations with the linear values used. Figures 2025 illustrate the curves of the displacement at the top end of the porcelain pillar and the maximum stress at the bottom end of the porcelain pillar under different seismic excitations with the nonlinear values used. From Figures 1425, the maximum stresses and displacements of the porcelain pillar with linear and nonlinear conditions can be obtained, as shown in Table 2. Figure 26 illustrates the dynamic response of the equipment under different seismic excitations with the linear and nonlinear values used.

As seen from the analysis results in Table 2 and Figure 26, for the 1000 kV arrester equipment, when the flange connection is linear, the maximum displacement and stress grow linearly with the increase of peak ground acceleration; when the nonlinear factor is considered, the maximum displacement and stress increase nonlinearly, and the magnitude of the increase gradually decreases. It can be seen from Figure 25 that the nonlinear seismic response of the equipment is smaller than the linear seismic response under seismic excitation. The greater the ground motion excitation, the larger the difference between nonlinear response results and linear response results of the equipment. It indicates that equipment exhibits more obvious nonlinear characteristics as seismic excitation increases.

When the peak ground acceleration increases from 0.1 g to 0.5 g, the maximum displacement of the top of the porcelain pillar increases from 22.26 mm to 111.28 mm without consideration of nonlinear factor; when the nonlinear factor is considered, the range of this value is from 22 mm to 95.06 mm. Correspondingly, the variation range of the maximum stress at the bottom end of the porcelain pillar is from 7.68 MPa to 38.41 MPa and 7.48 MPa to 29.48 MPa, respectively. When the peak ground acceleration is 0.1 g, 0.2 g, 0.3 g, 0.4 g, and 0.5 g, the percentages of the difference between nonlinear analysis results and linear analysis results of displacement are 1.17%, 2.88%, 5.72%, 9.66%, and 14.58% respectively, and the percentages of the difference between nonlinear analysis results and linear analysis results of stress are 2.60%, 6.77%, 11.93%, 17.51%, and 23.25%, respectively.

5. Characteristics of Nonlinear Seismic Response within Structural Parameter Domains

Figure 27 illustrates the curves for the influence of changes of nonlinear rigidity parameter within a certain range on the maximum stress at the bottom end of porcelain pillar when the seismic excitations are 0.1 g, 0.2 g, 0.3 g, 0.4 g, and 0.5 g. Table 3 shows the corresponding maximum stress value at the bottom end of the porcelain pillar with the change of nonlinear rigidity η under different seismic excitations.

As seen from Figure 27, for different seismic excitations, as the nonlinear rigidity increases, the maximum stress at the bottom end of the porcelain pillar exhibits change characteristic of decreasing first and then becoming stable. When the seismic excitation is 0.1 g and the nonlinear rigidity is 12 × 104, the maximum stress at the bottom end of porcelain is 6.67 MPa, which has a decrease of about 9.5% compared with that when the nonlinear rigidity is 0. The nonlinear rigidity change has a rather little impact on the maximum stress at the bottom end of the porcelain pillar. As seen from Table 3, the impact of nonlinear rigidity on the maximum stress at the bottom end of the porcelain pillar increases gradually when the seismic excitation increases. When the seismic excitation is 0.2 g, 0.3 g, 0.4 g, and 0.5 g, respectively, and the nonlinear rigidity is 12 × 104, the maximum stress at the bottom end of the porcelain pillar decreases by 23.8%, 26.6%, 35.5%, and 42.8% compared with that when the nonlinear rigidity is 0.

According to Section 3.3, η = 20,000 can reflect the nonlinear factors, and the maximum stress at the bottom end of the porcelain pillar is analyzed with η changes in the range of [0, 4 × 104] in order to consider the influence of fluctuation of η on the maximum stress at the bottom end of porcelain pillar. When the seismic excitation is 0.1 g, 0.2 g, 0.3 g, 0.4 g, and 0.5 g, respectively, and η increases from 0 to 2 × 104, the maximum stress at the bottom end of the porcelain pillar decreases by 1.9%, 6.7%, 13.43%, 21.2%, and 23.7%; when η increases from 2 × 104 to 4 × 104, the maximum stress at the bottom end of porcelain pillar decreases by 1.8%, 5.9%, 11.33%, 3.7%, and 4.01%, respectively. As seen from the above analysis, the decrease rate of the maximum stress at the bottom end of the porcelain pillar when η increases from 2 × 104 to 4 × 104 is smaller than that when η increases from 0 to 2 × 104.

As seen from Figure 27, as the seismic excitation increases, the nonlinear rigidity change ranges corresponding to rapid decrease stages of maximum stress at the bottom end of porcelain pillar are [0, 12 × 104], [0, 10 × 104], [0,8 × 104], [0, 4 × 104], and [0, 2 × 104], and the change ranges have a trend of gradual decrease. It indicates that the influence of nonlinear rigidity on the maximum stress at the bottom end of the porcelain pillar under high seismic excitation is larger and the maximum stress changes rapidly. When the nonlinear rigidity increases to a certain value, the impact of nonlinear rigidity on the maximum stress at the bottom end of the porcelain pillar is rather small, and the stress change of the porcelain pillar has entered stationary change.

6. Conclusions

(1)According to elastic mechanics, in the range of elasticity, it is assumed that the glued part of flange connection assembly before and after deformation follows the plain-section assumption; that is, the cross section after deformation at the glued part assembly is perpendicular to the neutral plane. A dynamic model of porcelain pillar electrical equipment was established for considering the nonlinear characteristics of flange connections. The theoretical analysis results and the test results have a rather good degree of the fitting, which verifies the correctness of the dynamic model. When the nonlinear rigidity parameter is 20,000, the influence of nonlinear factors at the cementing parts can be reflected(2)Both the nonlinear calculation results of the displacement at the top end of the porcelain pillar and the stress at the bottom end of the porcelain pillar are smaller than the linear calculation results. The difference between such two kinds of results increases gradually with the increase of seismic excitation. The differences in displacement and stress at a seismic excitation of 0.5 g are 14.58% and 23.25%, respectively(3)When the nonlinear rigidity increases to a certain value, the influence of nonlinear rigidity on the maximum stress at the bottom end of the porcelain pillar is rather small, the stress change of the porcelain pillar has entered stationary change. As the seismic excitation increases, the nonlinear rigidity change ranges corresponding to rapid decrease stages of maximum stress at the bottom end of the porcelain pillar decrease gradually

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was funded by the National Key R&D Program of China (2018YFC0809400).