# Research on NOFRF Entropy-Based Detection Method for Early Damage of Pillar Porcelain Insulator

Authors: Hanling Mao; Weili Tang; Yingxiang Huang; Hanying Mao; Zhenfeng Huang; Xinxin Li
Journal: Shock and Vibration
 (2020)

---

## Abstract

The pillar porcelain insulator is an important protection device that is related to the safety of the entire power grid. Small damage of it may even cause a disaster. Nonlinear output frequency response functions (NOFRFs) can well reflect the nonlinear characteristics of early damage in the system. This paper uses the NOFRF entropy-based harmonic excitation detection method to detect the damage of the structure. Its effect has been verified by using metal specimens. Then, the pulse hammer detection and harmonic detection are carried out to detect the early-stage damage of the pillar porcelain insulator, and they achieve a good result, which verifies that the harmonic detection method can detect the early-stage damage in the pillar porcelain insulator as well. In addition, the orthogonal test of the simulated breathing crack model is used to search the greatest influence of the parameters of crack on the detection index. Through the orthogonal analysis, the results show that among the lengths, positions, and angles of the crack, the length of crack is the main factor that affects the detection index.


---

## Body

## 1. Introduction

The pillar porcelain insulator is an insulating component that plays an important role in the power grid. During the service of the pillar porcelain insulator, cracks may occur due to the effects of wind load, sudden changes in temperature, and so on. The failure of the insulator would reduce the operational life of the entire line and even cause a disaster.Therefore, it is of great practical significance to study the online detection for early damage of insulators. The acoustic detection, ultraviolet detection, and laser Doppler vibration method introduced in [1–3] have the advantages of good real-time performance and high efficiency, but there are still problems such as being susceptible to environmental influences and expensive equipment. Because of the global nature of vibration detection, the vibration-based method is convenient and effective to detect structural fatigue damage. Vibration detection is widely used in structural nondestructive testing in the fields of construction and machinery [4–6]. In [7, 8], the vibration detection method of the pillar insulator based on two kinds of detection standards of natural frequency variation and characteristic frequency energy amplitude transfer was introduced, and the feasibility of vibration detection was verified by experiments. Furthermore, to find the damage of the structures in the early stage, some scholars proposed the vibration detection method based on nonlinear features to detect damage in the system [7–10]. Now, nonlinear methods have become a research hotspot because nonlinear algorithms are more sensitive to changes in damage.At present, there are many nonlinear analysis methods. The nonlinear output frequency response function (NOFRF) is an effective nonlinear algorithm, which was derived from the Volterra series. The Volterra series is widely used for damage detection [11–13]. The concept of NOFRFs was proposed by Lang and Billings in recent years [14, 15]. Peng and Lang [16–18] had done much work which has shown that NOFRFs can effectively describe the nonlinearity of the system and characterize the nonlinear characteristics of the damage, especially the early stage of the damage. After that, Mao et al. [19, 20] proposed a NOFRF estimation method based on an improved algorithm for the NARMAX model with rectangular pulse hammer detection based on NOFRF theory and built several detection indexes based on NOFRFs, which used the concepts of information entropy, complexity, and divergence, respectively. Through the experiment, the analysis results show that the pulse hammer detection and proposed damage detection indexes have a good performance for indicating the degree of fatigue damage, especially the index which used the concept of information entropy. The pulse hammer can excite a wide frequency range, which has a good effect on NOFRF estimation. Liu et al. [21, 22] established a weighted contribution rate as an index for NOFRFs. The index values of the rotor with different crack relative depths were estimated. The extracted fault features are enlarged in the method, which can be used in the diagnosis of the early-stage crack of the rotor. Liu et al. [23] also proposed a variable weighted contribution rate of NOFRFs, which is sensitive to the slight misalignment of the rotor system.However, when it comes to online monitoring, the shortcomings of the pulse hammer detection method are revealed. Pulse hammer excitation signal is not a continuous signal, so it is difficult to achieve real-time monitoring. Furthermore, in the actual detection process, the hammer excitation position and angle are difficult to control. So, a large number of repeated detections should be done to reduce the error.In this work, a NOFRF entropy-based harmonic excitation detection method was used to detect the damage of the structure. Its detection efficiency is higher than that of pulse hammer detection, and the detection cost is cheap as well. Unlike the pulse signal, the harmonic signal is a continuous excitation signal, which is more beneficial to online monitoring. Then, the harmonic excitation detection and hammering excitation detection were performed on the same specimens. The results of them were compared to verify that harmonic excitation detection is effective in estimating the NOFRF entropy index. When its effectiveness had been proved, it was applied to the pillar porcelain insulator. Finally, the greatest influence of crack on values of NOFRF detection indexes was explored by simulation.


## 2. Theoretical Background

### 2.1. Nonlinear Output Frequency Response Function

NOFRFs were derived from the Volterra series. Each order of NOFRFs can be estimated by using the frequency spectra of excitation signal and response signal, which are used as input and output, respectively. The definition [14] is given as(1)Gnjω=∫−∞∞⋯∫−∞∞Hnjω1,…,jωn∏i=1nUjω1dω1…dωn∫−∞∞⋯∫−∞∞∏i=1nUjω1dω1…dωn,1≤n≤N,where Hnjω1,⋯,jωn is the n-th order generalized frequency response function (GFRF); Ujω is the input; and Yjω is the output frequency response, which can be written as(2)Yjω=∑n=1NYnjω=∑n=1NGnjωUnjω.The above description can be expressed as Figure1.Figure 1
The expression of NOFRFs.In Figure1, uN(t) is the N power of u(t).G1(jω) represents the linear part of the system. And G2(jω), …, GN(jω) represent the nonlinear parts of the system.The general estimation method of NOFRFs was given by Peng and Lang [14] as follows.Ifn sets of excitation signals and response signals are obtained, equation (3) can be obtained.(3)Y1jω=G1jωα1U1∗jω+G2jωα12U2∗jω+⋯+Gnjωα1nUn∗jω,Y2jω=G1jωα2U1∗jω+G2jωα22U2∗jω+⋯+Gnjωα2nUn∗jω,⋮Ynjω=G1jωαnU1∗jω+G2jωαn2U2∗jω+⋯+Gnjωα2nUn∗jω,where Unjω represents the n-order nonlinear input of the system and U∗jω is the spectrum of specific input signal u∗t, in which ut=αu∗t.Then, the NOFRFs can be solved by the least squares method as shown in the following equation4:(4)G1jω,G2jω,…,GnjωT=αUjωTαUjω−1αUjωTYjω,where(5)αUjω=α1U1∗jω,α12U2∗jω,…,α1nUn∗jωα2U1∗jω,α22U2∗jω,⋯,α2nUn∗jω⋮αnU1∗jω,αn2U2∗jω,…,α2nUn∗jω,(6)Yjω=Y1jω,Y2jω,…,YnjωT.


### 2.2. Detection Index

Peng et al. [18] integrated the NOFRF values of each order to construct an index Fe with a comprehensive description of the characteristic changes of the nonlinear system, which is defined as(7)Fen=∫−∞∞Gnjωdω∑i=1N∫−∞∞Gijωdω,1≤n≤N.As can be seen from equation (7), the value of Fe(n) is in the range from zero to one. According to the meaning of G1(jω), Fe(1) indicates the linear degree of the system. Meanwhile, Fe(2), …, Fe(n) could indicate the nonlinear degree of the system, which represents the degree of damage. When the system was subjected to the damage, the ratio of higher-order Fe would increase.To build a more sensitive index, the concept of information entropy was used. Information entropy can describe the uncertainty of the numerical sequences. The definition of information entropy is(8)H=−∑i=1Npxilogpxi,where p(xi) represents the probability of a random variable xi.Based on the information entropy theory, the NOFRF entropy indexNE is proposed:(9)NE=−∑n=1NFennlogFenn.When the ratio of higher-order Fe increased, the value ofNE would increase as well. For damage detection analysis, the value of NE represents the severity of the damage.


## 2.1. Nonlinear Output Frequency Response Function

NOFRFs were derived from the Volterra series. Each order of NOFRFs can be estimated by using the frequency spectra of excitation signal and response signal, which are used as input and output, respectively. The definition [14] is given as(1)Gnjω=∫−∞∞⋯∫−∞∞Hnjω1,…,jωn∏i=1nUjω1dω1…dωn∫−∞∞⋯∫−∞∞∏i=1nUjω1dω1…dωn,1≤n≤N,where Hnjω1,⋯,jωn is the n-th order generalized frequency response function (GFRF); Ujω is the input; and Yjω is the output frequency response, which can be written as(2)Yjω=∑n=1NYnjω=∑n=1NGnjωUnjω.The above description can be expressed as Figure1.Figure 1
The expression of NOFRFs.In Figure1, uN(t) is the N power of u(t).G1(jω) represents the linear part of the system. And G2(jω), …, GN(jω) represent the nonlinear parts of the system.The general estimation method of NOFRFs was given by Peng and Lang [14] as follows.Ifn sets of excitation signals and response signals are obtained, equation (3) can be obtained.(3)Y1jω=G1jωα1U1∗jω+G2jωα12U2∗jω+⋯+Gnjωα1nUn∗jω,Y2jω=G1jωα2U1∗jω+G2jωα22U2∗jω+⋯+Gnjωα2nUn∗jω,⋮Ynjω=G1jωαnU1∗jω+G2jωαn2U2∗jω+⋯+Gnjωα2nUn∗jω,where Unjω represents the n-order nonlinear input of the system and U∗jω is the spectrum of specific input signal u∗t, in which ut=αu∗t.Then, the NOFRFs can be solved by the least squares method as shown in the following equation4:(4)G1jω,G2jω,…,GnjωT=αUjωTαUjω−1αUjωTYjω,where(5)αUjω=α1U1∗jω,α12U2∗jω,…,α1nUn∗jωα2U1∗jω,α22U2∗jω,⋯,α2nUn∗jω⋮αnU1∗jω,αn2U2∗jω,…,α2nUn∗jω,(6)Yjω=Y1jω,Y2jω,…,YnjωT.


## 2.2. Detection Index

Peng et al. [18] integrated the NOFRF values of each order to construct an index Fe with a comprehensive description of the characteristic changes of the nonlinear system, which is defined as(7)Fen=∫−∞∞Gnjωdω∑i=1N∫−∞∞Gijωdω,1≤n≤N.As can be seen from equation (7), the value of Fe(n) is in the range from zero to one. According to the meaning of G1(jω), Fe(1) indicates the linear degree of the system. Meanwhile, Fe(2), …, Fe(n) could indicate the nonlinear degree of the system, which represents the degree of damage. When the system was subjected to the damage, the ratio of higher-order Fe would increase.To build a more sensitive index, the concept of information entropy was used. Information entropy can describe the uncertainty of the numerical sequences. The definition of information entropy is(8)H=−∑i=1Npxilogpxi,where p(xi) represents the probability of a random variable xi.Based on the information entropy theory, the NOFRF entropy indexNE is proposed:(9)NE=−∑n=1NFennlogFenn.When the ratio of higher-order Fe increased, the value ofNE would increase as well. For damage detection analysis, the value of NE represents the severity of the damage.


## 3. Detection of Specimens

### 3.1. Specimens

In order to verify the effectiveness of the NOFRF entropy-based detection method for structure, a verification experiment was carried out. Four specimens were prepared. The specifications of the specimens are shown in Figure2 and Table 1. In the preparation of the specimens, the wire cutting was first performed. A wire cut was made in the middle of the specimen, as shown in Figure 2. Then, the surface of it was polished. Finally, three of them were subjected to a three-point bending fatigue loading test to make the prefabricated fatigue cracks.Figure 2
The schematic diagram of specimen.Table 1
The parameters of specimens.
Parameters of specimenValueSize130 mm × 30 mm × 15 mmMaterial45# steelWire cutting depth5 mmThe diameter of wire-cut wire0.18 mmLength of fatigue cracks0 mm, 3 mm, 10 mm, 15 mmIn order to choose the location of the suspension, the modal shapes of the specimen were analyzed by finite element software ABAQUS. The analysis results are shown in Figure3. There were different minimum displacement positions in different modes, which are shown as the blue area in Figure 3. The minimum displacement position of each mode was considered comprehensively, and two suspension positions were selected, as shown in Figure 4.Figure 3
Modal analysis of wire-cut notched metal specimens. (a) The first mode. (b) The second mode. (c) The third mode. (d) The fourth mode.
(a)(b)(c)(d)Figure 4
Schematic diagram of metal specimen suspension detection.


### 3.2. Harmonic Excitation Detection of Specimens

According to the result of modal analysis of the specimen, the harmonic excitation detection platform is built as shown in Figures4 and 5. In the platform, the sound vibration sensors could be used to realize both the excitation and acquisition of the given signal. Two sensors were pasted on both ends of the upper surface of the specimen. One of them was used to generate a vibration excitation signal, and the other one was used to acquire the response signal. The specimen was suspended by two rubber bands. The acquisition frequency range of sound vibration sensor is 0–800 kHz. The Tektronix function/arbitrary waveform generator (AFG3022C) is selected to generate a harmonic signal. Finally, both excitation and response signals were collected using the Tektronix mixed domain oscilloscope (MDO3022) at a sampling frequency of 50 kHz.Figure 5
Harmonic excitation detection platform.In order to calculate the first four orders of NOFRFs, four excitation signals with different amplitudes were required. In the experiment, harmonic signals with peak-to-peak values of 7 V, 8 V, 9 V, and 10 V were selected as input, respectively. And the excitation frequency was set to be 8963 Hz, the same as the natural frequency of the second-order modal frequency of the specimen. When the input and output data were collected, the values of detection indexes could be calculated through equations (3) and (4). After calculation, the experimental results are shown in Figure 6.Figure 6
Harmonic excitation damage detection index change curve.


### 3.3. Hammering Excitation Detection of Specimens

As a comparison with the results of harmonic detection, pulse hammer detection was performed. It had been proved that NOFRFs and the value ofNE index can be well estimated through the pulse hammer detection [19, 20].In this hammering excitation detection, the detection platform was built as shown in Figure7. The PULSE Lab shop general test system, 8206-002 impact hammer, 4508 one-way sensor, 3050-a-60 data acquisition module and 3099-a-n1 six-channel analyzer module from B&K Denmark were used in the detection. The sampling frequency was 5 kHz.Figure 7
Hammering excitation detection platform.The edge excitation signal was collected for 4 times using the force hammer axial hammer test specimen. Four pulse excitation signals with different amplitudes and their response signals were obtained. The least squares method was used to estimate NOFRF of each order as equations (3)–(6). The damage detection index value was calculated through equations (7) and (9). The calculation results are shown in Figure 8.Figure 8
Change curve of damage detection index value under hammer excitation.As can be seen from Figures6 and 8, under harmonic excitation and hammer excitation, the changes of specimen damage can be directly reflected by the detection indexes NE and Fe. When the damage became heavier, the value of Fe(1) decreased with the depth of crack, while the values of high-order Fe and the value of NE increased. It met the theoretical prediction.Through the above analysis, the validity of the NOFRF entropy-based harmonic detection method was verified.


## 3.1. Specimens

In order to verify the effectiveness of the NOFRF entropy-based detection method for structure, a verification experiment was carried out. Four specimens were prepared. The specifications of the specimens are shown in Figure2 and Table 1. In the preparation of the specimens, the wire cutting was first performed. A wire cut was made in the middle of the specimen, as shown in Figure 2. Then, the surface of it was polished. Finally, three of them were subjected to a three-point bending fatigue loading test to make the prefabricated fatigue cracks.Figure 2
The schematic diagram of specimen.Table 1
The parameters of specimens.
Parameters of specimenValueSize130 mm × 30 mm × 15 mmMaterial45# steelWire cutting depth5 mmThe diameter of wire-cut wire0.18 mmLength of fatigue cracks0 mm, 3 mm, 10 mm, 15 mmIn order to choose the location of the suspension, the modal shapes of the specimen were analyzed by finite element software ABAQUS. The analysis results are shown in Figure3. There were different minimum displacement positions in different modes, which are shown as the blue area in Figure 3. The minimum displacement position of each mode was considered comprehensively, and two suspension positions were selected, as shown in Figure 4.Figure 3
Modal analysis of wire-cut notched metal specimens. (a) The first mode. (b) The second mode. (c) The third mode. (d) The fourth mode.
(a)(b)(c)(d)Figure 4
Schematic diagram of metal specimen suspension detection.


## 3.2. Harmonic Excitation Detection of Specimens

According to the result of modal analysis of the specimen, the harmonic excitation detection platform is built as shown in Figures4 and 5. In the platform, the sound vibration sensors could be used to realize both the excitation and acquisition of the given signal. Two sensors were pasted on both ends of the upper surface of the specimen. One of them was used to generate a vibration excitation signal, and the other one was used to acquire the response signal. The specimen was suspended by two rubber bands. The acquisition frequency range of sound vibration sensor is 0–800 kHz. The Tektronix function/arbitrary waveform generator (AFG3022C) is selected to generate a harmonic signal. Finally, both excitation and response signals were collected using the Tektronix mixed domain oscilloscope (MDO3022) at a sampling frequency of 50 kHz.Figure 5
Harmonic excitation detection platform.In order to calculate the first four orders of NOFRFs, four excitation signals with different amplitudes were required. In the experiment, harmonic signals with peak-to-peak values of 7 V, 8 V, 9 V, and 10 V were selected as input, respectively. And the excitation frequency was set to be 8963 Hz, the same as the natural frequency of the second-order modal frequency of the specimen. When the input and output data were collected, the values of detection indexes could be calculated through equations (3) and (4). After calculation, the experimental results are shown in Figure 6.Figure 6
Harmonic excitation damage detection index change curve.


## 3.3. Hammering Excitation Detection of Specimens

As a comparison with the results of harmonic detection, pulse hammer detection was performed. It had been proved that NOFRFs and the value ofNE index can be well estimated through the pulse hammer detection [19, 20].In this hammering excitation detection, the detection platform was built as shown in Figure7. The PULSE Lab shop general test system, 8206-002 impact hammer, 4508 one-way sensor, 3050-a-60 data acquisition module and 3099-a-n1 six-channel analyzer module from B&K Denmark were used in the detection. The sampling frequency was 5 kHz.Figure 7
Hammering excitation detection platform.The edge excitation signal was collected for 4 times using the force hammer axial hammer test specimen. Four pulse excitation signals with different amplitudes and their response signals were obtained. The least squares method was used to estimate NOFRF of each order as equations (3)–(6). The damage detection index value was calculated through equations (7) and (9). The calculation results are shown in Figure 8.Figure 8
Change curve of damage detection index value under hammer excitation.As can be seen from Figures6 and 8, under harmonic excitation and hammer excitation, the changes of specimen damage can be directly reflected by the detection indexes NE and Fe. When the damage became heavier, the value of Fe(1) decreased with the depth of crack, while the values of high-order Fe and the value of NE increased. It met the theoretical prediction.Through the above analysis, the validity of the NOFRF entropy-based harmonic detection method was verified.


## 4. Detection of Pillar Porcelain Insulators

### 4.1. Experimental Setup and Process

In order to verify the validity of the detection method for damage detection of the pillar porcelain insulator, the harmonic and pulse excitation detection experiments were carried out. Two pillar porcelain insulators were prepared. One of them is without any defect, and the other one has a prefabricated crack inside the flange. The insulator sample is presented in Figure9.Figure 9
The sample of pillar porcelain insulator.During the experiments, the insulators were placed vertically on the ground, and the self-weight of the insulators was used to simulate the field installation state of the pillar porcelain insulator. The experimental platform is shown in Figures10 and 11, respectively. The devices used were the same as mentioned in Section 3.Figure 10
Harmonic excitation experiment.Figure 11
Hammering excitation experiment.The single-frequency harmonic excitation signal and the hammer excitation signal were, respectively, used as inputs. And the response signals were collected as output. Then, NOFRFs could be identified, and the nonlinear detection indexesNE and Fe could be calculated.In the hammering excitation experiment, the lower flange of the insulator was hammered four times by force hammer from the axial direction. Pulse hammering signals were chosen as input. The response signals were used as output signals. And the input and output signals were collected with a sampling frequency of 5 kHz.In the harmonic excitation experiment, harmonic signals with peak-to-peak values of 7 V, 8 V, 9 V, and 10 V were selected as input. The frequency of harmonic excitation signal was 1440 Hz which was close to the third modal frequency, and the sampling frequency was 16 kHz. Input and output signals were collected by the Tektronix mixed domain oscilloscope.The values of detection indexes were calculated as the algorithm described in Section2.


### 4.2. Result and Discussion of Experiment

The signal and its frequency spectrum collected from the accelerator in hammer detection are shown in Figure12.Figure 12
The response signal and its frequency spectrum in hammer detection. (a) The response signal. (b) The frequency spectrum of the response signal.
(a)(b)Then, the NOFRFs were calculated according to equations (3) and (4). The results are shown in Figures 13 and 14.Figure 13
The NOFRFs of healthy pillar porcelain insulator in hammer detection.Figure 14
The NOFRFs of cracked pillar porcelain insulator in hammer detection.The frequency spectra of collected signals in harmonic excitation detection are shown in Figure15.Figure 15
The spectrum of response signal in the harmonic detection. (a) The spectrum of response signal of healthy pillar porcelain insulator. (b) The spectrum of response signal of cracked pillar porcelain insulator.
(a)(b)The results of hammer and harmonic excitation experiments were obtained as shown in Tables2 and 3, respectively.Table 2
Damage detection index value of pillar porcelain insulator under harmonic excitation.
Fe (1)Fe (2)Fe (3)Fe (4)NEHealthy0.855030.134110.009560.001300.52484Prefabricated crack0.815840.172740.0096910.001740.61717Table 3
Damage detection index value of pillar porcelain insulator under hammer excitation.
Fe (1)Fe (2)Fe (3)Fe (4)NEHealthy0.952510.046950.000552.34E − 60.21664Prefabricated crack0.909880.088900.001215.10E − 060.35388From the tables, we can see that in harmonic excitation detection, the relative growth rate ofNE value after damage was 17.59%, while in hammer excitation detection, the relative growth rate of NE value after damage was 63.35%. It means that the sensitivity of index NE in hammer excitation detection was higher than that in harmonic excitation detection.Although their sensitivity is different, they can reflect the change of damage of the pillar porcelain insulator.


## 4.1. Experimental Setup and Process

In order to verify the validity of the detection method for damage detection of the pillar porcelain insulator, the harmonic and pulse excitation detection experiments were carried out. Two pillar porcelain insulators were prepared. One of them is without any defect, and the other one has a prefabricated crack inside the flange. The insulator sample is presented in Figure9.Figure 9
The sample of pillar porcelain insulator.During the experiments, the insulators were placed vertically on the ground, and the self-weight of the insulators was used to simulate the field installation state of the pillar porcelain insulator. The experimental platform is shown in Figures10 and 11, respectively. The devices used were the same as mentioned in Section 3.Figure 10
Harmonic excitation experiment.Figure 11
Hammering excitation experiment.The single-frequency harmonic excitation signal and the hammer excitation signal were, respectively, used as inputs. And the response signals were collected as output. Then, NOFRFs could be identified, and the nonlinear detection indexesNE and Fe could be calculated.In the hammering excitation experiment, the lower flange of the insulator was hammered four times by force hammer from the axial direction. Pulse hammering signals were chosen as input. The response signals were used as output signals. And the input and output signals were collected with a sampling frequency of 5 kHz.In the harmonic excitation experiment, harmonic signals with peak-to-peak values of 7 V, 8 V, 9 V, and 10 V were selected as input. The frequency of harmonic excitation signal was 1440 Hz which was close to the third modal frequency, and the sampling frequency was 16 kHz. Input and output signals were collected by the Tektronix mixed domain oscilloscope.The values of detection indexes were calculated as the algorithm described in Section2.


## 4.2. Result and Discussion of Experiment

The signal and its frequency spectrum collected from the accelerator in hammer detection are shown in Figure12.Figure 12
The response signal and its frequency spectrum in hammer detection. (a) The response signal. (b) The frequency spectrum of the response signal.
(a)(b)Then, the NOFRFs were calculated according to equations (3) and (4). The results are shown in Figures 13 and 14.Figure 13
The NOFRFs of healthy pillar porcelain insulator in hammer detection.Figure 14
The NOFRFs of cracked pillar porcelain insulator in hammer detection.The frequency spectra of collected signals in harmonic excitation detection are shown in Figure15.Figure 15
The spectrum of response signal in the harmonic detection. (a) The spectrum of response signal of healthy pillar porcelain insulator. (b) The spectrum of response signal of cracked pillar porcelain insulator.
(a)(b)The results of hammer and harmonic excitation experiments were obtained as shown in Tables2 and 3, respectively.Table 2
Damage detection index value of pillar porcelain insulator under harmonic excitation.
Fe (1)Fe (2)Fe (3)Fe (4)NEHealthy0.855030.134110.009560.001300.52484Prefabricated crack0.815840.172740.0096910.001740.61717Table 3
Damage detection index value of pillar porcelain insulator under hammer excitation.
Fe (1)Fe (2)Fe (3)Fe (4)NEHealthy0.952510.046950.000552.34E − 60.21664Prefabricated crack0.909880.088900.001215.10E − 060.35388From the tables, we can see that in harmonic excitation detection, the relative growth rate ofNE value after damage was 17.59%, while in hammer excitation detection, the relative growth rate of NE value after damage was 63.35%. It means that the sensitivity of index NE in hammer excitation detection was higher than that in harmonic excitation detection.Although their sensitivity is different, they can reflect the change of damage of the pillar porcelain insulator.


## 5. The Influence of Crack Parameters on the NOFRF Entropy Index

The vibration detection is a global detection method, which can detect cracks in almost all locations of the structure. However, different parameters of the crack in the structural system will affect the estimated results of the detection index. In order to investigate the influence degree of different crack parameters on the value of indexes in the detection, a simulation was carried out.


### 5.1. Simulation Settings

At present, research studies on damage detection of pillar porcelain insulator are mostly focused on the open defect [24]. In this case, the nonlinear contact effect caused by crack opening and closing behavior would be ignored. Therefore, a representative finite element model of breathing crack should be established. The contact area between two surfaces of the crack changes under the external excitation, and the one-piece combination was similar to the dynamic process of respiration [25].At first, the 3D model of pillar porcelain insulator was built by SolidWorks as shown in Figure16. Then, the finite element model of pillar porcelain insulator was imported into ABAQUS. The type of grid element is C3D8R. The material parameters and main dimensions of the pillar porcelain insulator are shown in Tables 4 and 5, respectively.Figure 16
3D model of pillar porcelain insulator.Table 4
Material parameters of pillar porcelain insulator.
MaterialE (MPa)Poisson’s ratioDensity,t (mm³)Porcelain1180000.162.6e − 9Cast iron1750000.287.8e − 9Cement310000.2152.3e − 9Table 5
Main dimensions of pillar porcelain insulator.
Stem diameter (mm)High (mm)Big umbrella diameter (mm)Small umbrella diameter (mm)Number of umbrellas59117512510827The breathing crack model was used to simulate the damage of the pillar porcelain insulator. The breathing crack would be set at one of three positions which are I, II, or III in Figure16. The big umbrellas and small umbrellas are also shown in Figure 16.The shape of breathing cracks is shown in Figure17. Δr is the crack depth. L is the crack width.Figure 17
Breathing crack model of pillar porcelain insulator.In ABAQUS, the crack was set in the porcelain body of the insulator model. The noncracked area was bounded by using the TIE tool. There was no relative motion and deformation between the bounded regions. The two opposing crack areas were set to be the normal hard contact by using the CONTACT tool. Crack faces were set to be independent. The tangential coefficient of friction was set to be 0.5. The Newmark implicit algorithm was used for time integration. When the external load was applied, the dynamic behavior of two crack faces was similar to the breathing action.


### 5.2. Result and Discussion of Simulation

In order to explore the influences on theNE value, different conditions were set. Different positions, lengths, and angles of breathing cracks were detected and analyzed through the orthogonal test. In the test, the triangular pulse signal was used as input to excite the pillar porcelain insulator with different conditions which are shown in Table 6. The calculation results are shown in Table 7.Table 6
Factor and level table of the test.
LevelCrack location (A)Crack length (with fixed depth 5 mm) (mm) (B)Angle of axial crack (°) (C)1Top flange3302Insulator middle5603Lower flange1590Table 7
Orthogonal test sets and visual analysis.
NumberA (location)B (mm)C (°)D (error term)NE111110.1379212220.1629313330.3089421230.0948522310.1275623120.3864731320.1141832130.2227933210.3350k10.20320.11560.24900.2001k20.20290.17100.19760.2211k30.22390.34340.18350.2088Range0.02070.22790.06550.0210In Table7, k1 is the average value of the factors under level 1, k2 is the average value of the factors under level 2, and k3 is the average value of the factors under level 3. The range is the difference between the maximum value and the minimum value of the average value.The result shows that the range ofB was larger than that of A, C, and D, which means that the length of breathing crack had the greatest influence on the NE value among these three factors. The axial angle and location of crack had a little influence on NE value.The probability values were obtained by the statistical method of a significance test, which are shown in Table8. When the probability value of a factor is less than 0.05, the influence of this factor is very significant.Table 8
Analysis of variance.
Sources of varianceDegrees of freedomDeviation sum of squaresMean sum of squareF-critProbability valueA20.0008690.0004351.300.434B20.0847070.042353126.850.008C20.0071290.00356510.680.086Error20.0006680.000334Summation80.093373As we can see from Table8, only the probability value of breathing crack length is less than 0.05, which means that the influence of crack length on the value of index NE was very significant.


## 5.1. Simulation Settings

At present, research studies on damage detection of pillar porcelain insulator are mostly focused on the open defect [24]. In this case, the nonlinear contact effect caused by crack opening and closing behavior would be ignored. Therefore, a representative finite element model of breathing crack should be established. The contact area between two surfaces of the crack changes under the external excitation, and the one-piece combination was similar to the dynamic process of respiration [25].At first, the 3D model of pillar porcelain insulator was built by SolidWorks as shown in Figure16. Then, the finite element model of pillar porcelain insulator was imported into ABAQUS. The type of grid element is C3D8R. The material parameters and main dimensions of the pillar porcelain insulator are shown in Tables 4 and 5, respectively.Figure 16
3D model of pillar porcelain insulator.Table 4
Material parameters of pillar porcelain insulator.
MaterialE (MPa)Poisson’s ratioDensity,t (mm³)Porcelain1180000.162.6e − 9Cast iron1750000.287.8e − 9Cement310000.2152.3e − 9Table 5
Main dimensions of pillar porcelain insulator.
Stem diameter (mm)High (mm)Big umbrella diameter (mm)Small umbrella diameter (mm)Number of umbrellas59117512510827The breathing crack model was used to simulate the damage of the pillar porcelain insulator. The breathing crack would be set at one of three positions which are I, II, or III in Figure16. The big umbrellas and small umbrellas are also shown in Figure 16.The shape of breathing cracks is shown in Figure17. Δr is the crack depth. L is the crack width.Figure 17
Breathing crack model of pillar porcelain insulator.In ABAQUS, the crack was set in the porcelain body of the insulator model. The noncracked area was bounded by using the TIE tool. There was no relative motion and deformation between the bounded regions. The two opposing crack areas were set to be the normal hard contact by using the CONTACT tool. Crack faces were set to be independent. The tangential coefficient of friction was set to be 0.5. The Newmark implicit algorithm was used for time integration. When the external load was applied, the dynamic behavior of two crack faces was similar to the breathing action.


## 5.2. Result and Discussion of Simulation

In order to explore the influences on theNE value, different conditions were set. Different positions, lengths, and angles of breathing cracks were detected and analyzed through the orthogonal test. In the test, the triangular pulse signal was used as input to excite the pillar porcelain insulator with different conditions which are shown in Table 6. The calculation results are shown in Table 7.Table 6
Factor and level table of the test.
LevelCrack location (A)Crack length (with fixed depth 5 mm) (mm) (B)Angle of axial crack (°) (C)1Top flange3302Insulator middle5603Lower flange1590Table 7
Orthogonal test sets and visual analysis.
NumberA (location)B (mm)C (°)D (error term)NE111110.1379212220.1629313330.3089421230.0948522310.1275623120.3864731320.1141832130.2227933210.3350k10.20320.11560.24900.2001k20.20290.17100.19760.2211k30.22390.34340.18350.2088Range0.02070.22790.06550.0210In Table7, k1 is the average value of the factors under level 1, k2 is the average value of the factors under level 2, and k3 is the average value of the factors under level 3. The range is the difference between the maximum value and the minimum value of the average value.The result shows that the range ofB was larger than that of A, C, and D, which means that the length of breathing crack had the greatest influence on the NE value among these three factors. The axial angle and location of crack had a little influence on NE value.The probability values were obtained by the statistical method of a significance test, which are shown in Table8. When the probability value of a factor is less than 0.05, the influence of this factor is very significant.Table 8
Analysis of variance.
Sources of varianceDegrees of freedomDeviation sum of squaresMean sum of squareF-critProbability valueA20.0008690.0004351.300.434B20.0847070.042353126.850.008C20.0071290.00356510.680.086Error20.0006680.000334Summation80.093373As we can see from Table8, only the probability value of breathing crack length is less than 0.05, which means that the influence of crack length on the value of index NE was very significant.


## 6. Conclusion

This paper proposed a NOFRF and information entropy-based harmonic detection method, and several experiments and simulation were carried out. Based on the result of the validation experiment, detection of pillar porcelain insulators, and orthogonal analysis in this paper, the following conclusions can be drawn:(1)
The result of NOFRF entropy index in harmonic detection was similar to the result of hammer detection, which proved that the harmonic detection method can effectively detect and evaluate the damage of the pillar porcelain insulator. Although their sensitivity was different, they can reflect the change of damage of the structural system.(2)
In the orthogonal analysis, the results of NOFRF index under different lengths, positions, and angles of the crack in the structure were obtained. It was found that the crack length was the main factor that affects the detection index.



---

*Source: 2841254-2020-07-06.xml*