Abstract

In this work, the strain measurement accuracy of single-mode fibre (SMF) under thermal and vibration loads is investigated by strain-frequency shift coefficient analyses. This research allows for the application of SMF sensors for structural health monitoring in real operational conditions. The strain measurement accuracy under combined static and thermal load is investigated experimentally, which demonstrated that temperature fluctuations induce non-negligible errors in the strain measurement, even with temperature compensation applied. The temperature fluctuation range which can induce measurement errors is quantified as less than −20C or higher than 55C. In addition, a fatigue experiment is conducted to investigate the measurement accuracy under low-frequency vibration load. The results of the fatigue experiment demonstrate that the vibrations mainly increase the ratio of null values in strain measurements. Findings from experiments can be applied to enhance structural health monitoring accuracy and reduce false positives. This study has important implications for the service application of distributed optical fibre sensing for composite structure health monitoring.

1. Introduction

Carbon fibre reinforced polymers (CFRPs) have been widely used in various aircraft structures [1, 2], due to their outstanding advantages including light weightedness, superior strength, and low costs. However, durability, unforeseen, and complex damage scenarios are challenges that require attention. Damage types in CFRP include delamination, debonding, fibre breakage, and matrix cracking and are difficult to characterize and detect [3, 4], which makes the monitoring of the integrity of composite parts particularly important. Strain measurement is an important method for structural health monitoring (SHM) because the strain distribution usually changes when damage is initiated in a structure [5, 6]. An ideal sensor for strain measurement in CFRP is the single-mode fibre (SMF) which is interrogated by an optical frequency domain reflectometry (OFDR) system [7]. This is because SMFs are light, cheap, easily available, and enable distributed strain sensing with high spatial accuracy.

Various fibre optical (FO) sensors have been adopted in different SHM applications [79]. For example, buckling detection based on distributed FO sensors [10], a hybrid piezoelectric-fibre optic system based on fibre Bragg gratings (FBG) [11, 12], embedded FBG sensors in SHM of complex composite structures [13], and sandwich panel [14]. In addition, FO sensors are adopted in composite manufacturing monitoring [15, 16] and damage detection [17, 18]. For in-field applications, however, many more FBG sensors are required than distributed FO sensors to obtain sufficient measurements, which results in higher system weight and cost.

In OFDR systems, Rayleigh backscattering of SMFs is used to interfere with the incident light to produce a beat signal. The frequency and amplitude of this beat signal correspond to a specific location on the SMF and the Rayleigh backscattering intensity at that location, respectively, since the light frequency of Rayleigh backscattering does not change during the previous processes. Moreover, the intensity of this backscattering is related to defects introduced during the manufacturing process. The beat signal is adopted as a unique code to identify a position on the SMF [19]. When this position is subjected to an external load, the resulting strain causes a frequency shift of the beat signal which is then calculated by matching the beat signals before and after loading. The strain is finally measured using a determined strain-frequency shift coefficient. This coefficient is the ratio between strain shift and corresponding frequency shift.

SMFs interrogated by OFDR have been adopted as strain sensors in many studies for various analyses of composite structures [5]. For example, shape sensing of a beam-like composite structure [20], damage detection in composite wind turbines [21], multiparameter strain sensing of composites [22], hygrothermal degradation of single fibre composite [23], and health state assessment of composite laminates during manufacturing [24]. Furthermore, the small size of SMF makes it suitable for application to complex composite structures including cylinders [25], bonded lap joints [26], and CFRP plates with circular holes [27]. As the measurement results of SMFs are influenced by both mechanical and thermal loads, many studies have been conducted on temperature compensation to improve strain measurement accuracy [28]. However, the temperature compensation cannot eliminate the strain measurement error induced by strain transfer coefficient change due to temperature fluctuations [29]. In addition, the adhesion and embedding quality of SMF is also proved to be affected by mechanical loading [30, 31]. Factors including temperature fluctuations and mechanical vibrations are inevitable in field applications, which induce non-negligible errors in strain measurement even with temperature compensation applied. Therefore, the investigation of strain measurement accuracy under thermal and vibration loads is essential for the field applications of SMF-based distributed sensing in CFRP.

Several researchers have validated the FO sensors for strain measurement under quasistatic loading [26, 28, 32]. However, the operational conditions of aircraft results in low-frequency vibration in-service in which the SMF sensors will also be exposed to. Therefore, the accuracy of distributed strain sensing under vibration load needs to be demonstrated. When calculating the strain based on distributed SMFs, the strain-frequency shift coefficient is usually obtained in a straightforward way by reading the frequency shift under a known load [3336]. This method is unable to obtain high-precision strain information in a complex environment such as during flight. To address this problem, this paper has studied the strain measurement accuracy of distributed SMF under thermal and vibration loads for the first time. The results show that environmental condition is one of the critical factors that need to be considered for higher accuracy of strain sensing with distributed SMFs. It also demonstrates that at a certain temperature range, the strain-frequency shift coefficient can be considered constant, but outside its range, it needs to be compensated for. The results of this work can be used as important references to service applications using distributed SMFs in composite structure health monitoring.

In this paper, the measurement error of SMF under thermal and vibration loads is investigated by strain-frequency shift coefficient analyses. The SMF and strain gauges are surface mounted on a composite laminate coupon. The strain gauge measurements are used as the reference strain. The fluctuation of the strain-frequency shift coefficient under thermal load is investigated by a static loading experiment first, followed by a fatigue experiment to investigate the fluctuation of this coefficient under low-frequency vibration load. The influence of temperature and vibration loads on the frequency shift is studied in combination and separately to isolate their effect on the strain measurement accuracy.

2. Experimental Setup

2.1. Specimen and Sensors

The size of the specimen manufactured and tested is  mm. This specimen is made of CYTEC noncrimp fabric. The layout of this specimen is . The thickness is 3.2 mm. The elastic properties of the lamina used in this paper are shown in Table 1.

The single mode fibre (SM1500 (9/125)P manufactured by FIBERCORE) and strain gauges are both adopted in the experiment for strain measurement. The SM1500 (9/125)P with polyimide coating has a variety of features such as high-temperature survival, bend insensitivity, and enhanced photosensitivity, which makes it very attractive for distributed strain sensing. In this work, the single mode fibres were adopted for strain measurement of composite laminates under both mechanical load (static and low-frequency dynamic load) and thermal load. The frequency shift of SMF can be related to a temperature or strain change. This relationship is expressed in equation as follows:where is the frequency shift; and are temperature and strain changes, respectively; express the temperature-frequency shift coefficient (or known as conversion factors from frequency shift to temperature change); is the strain-frequency shift coefficient (or known as conversion factors from frequency shift to strain). The detailed expressions of these two coefficients are reported in the study [37]. Both and are usually regarded as constant in the reported applications. For example, this equation works in a temperature range of −40C to 220C or in a strain range of −10000 –10000  for the ODiSI-B system manufactured by LUNA Ltd. However, the experiment results in this paper demonstrate that is not a constant under different ambient temperatures, which may induce strain measurement errors in real operational applications.

The SMF and strain gauges are bonded to the surface of the specimen by a thermoplastic adhesive film based on modified polyolefins (Pontacol 22.100). This type of thermoplastic adhesive film is adopted because of advantages including easy maintenance, high peel strength, and controllable thickness [17, 3840]. As shown in Figure 1, two “S” shape SMF, called SMF1 and SMF2, are mounted on each side of the specimen independently, to provide strain readings in different directions. A strain gauge is mounted on each side of the specimen: on side 1 at with respect to the fibre direction, defined as Gauge1 (as shown in Figure 1(a)), and at on side 2 (as shown in Figure 1(b)) defined at Gauge2. The data acquired by these two strain gauges are used as a reference strain for the strain-frequency shift coefficient calculation.

2.2. Static Loading Experiment

In the static loading experiment, a simultaneous loading of temperature and tension was applied to the specimen. This load was realized by the combination of an Instron dual-column machine (Instron 50 kN 5960 system) and a fitted environmental chamber (Instron 3119–600) shown in Figure 2(a). Figures 2(b) and 2(c) show sides 1 and 2 of the specimen in the test setup. The SMF data were acquired by the ODiSI-B system manufactured by LUNA Ltd.

Two types of tests were carried out as follows: the first test scenario is a measurement of frequency shift with a tension load, under different temperatures. In this test, the temperature changed from −30C to 80C in steps of 5C. Each time when the temperature stabilised, the reference spectrum of the SMF sensor was reset. A 0 to 10 kN (625 MPa) tensile load was then applied to the specimen, and the frequency shift of SMF was acquired. Therefore, this frequency shift is only caused by tension load under different temperatures. Another test scenario was to continuously measure the frequency change with tension and temperature variation. In this test, the temperature changed from −30C to 80C in steps of 5C. The reference spectrum of SMF is only set once at −30C. A 0 to 10 kN (625 MPa) tension load is directly applied to the specimen each time after the temperature became stable. Therefore, the measured frequency shift in this test scenario contains the combined effects of both tension and temperature variation.

2.3. Fatigue Experiment

In the fatigue experiment, a 5 Hz sine wave signal was used as a vibration load. The range of this load was set between 0 and 3 kN (187.5 MPa). This vibration load was realized by an Instron 8032 servo-hydraulic universal testing machine. The fatigue experiment lasted for a whole week to reach more than 1.6 million fatigue cycles. This experiment was carried out at room temperature. The strain gauge data were measured automatically every 1000 fatigue cycles. As for the SMF data, only SMF2 was measured manually about every hour during the day. This is because the ODiSI-B system cannot support multiple channels during dynamic measurements. The main purpose of the fatigue experiment is the applicability of the SMF under continuous dynamic measurement and frequency-strain coefficient stability under low frequency vibration load.

3. Data Analyses of Static Loading Experiment

3.1. Data Cleaning and Reconstruction

In the static loading experiment, many outliers and null values are observed in the measurement data acquired by the ODiSI-B system because of the similarity degenerate between measurement and reference Rayleigh backscattering under high load, which can significantly affect the results [41, 42]. Therefore, an algorithm is developed for data cleaning and reconstruction. The data of SMF2 at −30C in an independent test are used to illustrate the data reconstruction algorithm. The first step is outlier elimination. This is realized by the local outlier factor (LOF) algorithm [43]. This method is adopted because of the continuity of strain measurements. A constant LOF ( in equation (2)) is sufficient to achieve the outlier elimination in this experiment since the outliers are orders of magnitude greater than other data points. Assuming that the measurements of SMF are , the data after the elimination with the LOF algorithm is shown in equation as follows:

After the elimination, the data of the surface-mounted area is extracted according to peak alignment. Since the spatial resolution of the ODiSI-B system is 2.6 mm, 73 sensing segments are included along the fibre length according to its layout. The peak value of measurement data indicates the sensing segment along the loading direction in the static loading experiment (around the 25th and 49th sensing segment as shown in Figure 3(a)). Therefore, the data along the fibre orientation are extracted according to the locations of the peak values. The data after filtering and extraction are shown in Figure 3(b). The x-axis shows the spatial distribution of the measurement segments of SMF, which range from 0 to 0.2 m; the y-axis shows the static load from 0 to 10 kN; the z-axis plots the frequency shift of ODiSI-B at each measurement segment under each static load.

Null values are observed after outlier elimination, which implies a significant information loss. Temporal-based and spatial-based reconstruction methods [44] are adopted to address this problem since the measurement data have a clear spatial and temporal pattern. In the time domain, the linear fit between load and frequency shift is adopted first for data reconstruction. This is because of the linear relationship between the tension load and the measurements for every measurement segment on SMF. All missing data were reconstructed after the linear fit in the time domain. The results are shown in Figure 4(a). However, because of the low signal-noise ratio and measurement errors, there still exist errors after the linear fit. These errors (discontinuous value) are clearly indicated in Figure 4(a). Therefore, the robust locally weighted scatter-plot smoothing (RLOWESS) algorithm is adopted to eliminate these errors and to smooth the data in the temporal and spatial domain. The discrete measurements of SMF measurements after outlier elimination and linear fit reconstruction are . The measurements within a local sliding window at are adopted for the RLOWESS algorithm. The weighted function adopted in this method is given by the symmetric tricube function shown as follows:where is the weight of the measurement in the linear regression; is the centre of the local sliding window, which has the highest weight in linear regression; is the measurements within the local sliding window; is the Euclidean length of the local sliding window. Therefore, the sum of the residual squares after weighted linear regression is expressed by equation as follows:where and are the estimated coefficients of weighted linear regression. Therefore, the optimal solution of and is obtained by minimising , as expressed in equation as follows:

Following by the weighted linear regression mentioned previously, a modified weighted function based on regression residual is adopted to enhance the robustness of this algorithm. This modified weighted function is shown in equation as follows:where expresses the regression residuals; is the median absolute deviation of the residuals. The result after data cleaning and reconstruction is shown in Figure 4(b).

It is worth noting that the smoothing process mentioned previously is intended to restore the outliers and null values that are a result of systematic errors of ODiSI-B. However, this process may inadvertently mask crucial strain features in proximity to the damaged region, thereby negatively impacting certain applications such as damage detection.

3.2. Temperature Compensation

The data cleaning and reconstruction algorithm successfully restored the results in both test scenarios. As mentioned previously, the recorded data in both tests are different: the continuous measurements from the second test scenario contain tension and temperature information, whereas measurements from the first test scenario are mainly influenced by tension load. To illustrate the measurement consistency of both scenarios, temperature compensation is applied to the continuous measurement test based on the data recorded at the part of the SMF, which is only exposed to thermal load. The continuous measurements from the second test scenario are compensated and then compared with the measurements from the first test scenario. The comparison result under 80C is shown in Figure 5, where Figure 5(a) shows the results of distributed measurements at a 6 kN (375 MPa) tension load. Figure 5(b) shows the frequency shift along the load direction, under different tension loads. Note that the data of the 49th segment on SMF2 (Figure 3(a)) are adopted in this figure. This comparison illustrates that the compensated results of both test scenarios are highly consistent. Therefore, the measurements from the first test scenario are regarded as continuous measurements with temperature compensation in the following analyses.

3.3. Strain Measurement Accuracy under Thermal Load

The measurements from the first test scenario after data cleaning and reconstruction are adopted in strain calculation based on the previous analyses. Since these measurements are in the form of frequency shift of the SMF, the selection of the strain-frequency shift coefficient directly influences the strain results. The constant strain-frequency shift coefficient is usually used under different measurement states in real applications. However, the question raises that whether variability in operational conditions will also influence this coefficient. In other words, if the strain-frequency shift is calculated at room temperature, e.g., 20C, can this coefficient be used to convert the SFM readings recorded at a different temperature into strain? Therefore, in this work, the strain-frequency shift coefficient at different ambient temperatures is investigated in detail.

The strain-frequency shift coefficient is calculated based on measurements of SMF2 and Gauge2, with Gauge2 measurements being adopted as the reference strain. Therefore, the reliability of these strain gauge measurements at different temperatures is first investigated. Under the experimental conditions in this paper, this stability can be directly reflected by the variation in the strain-load coefficient of Gauge2. This strain-load coefficient is the ratio between the strain shift and the corresponding load shift. This coefficient should be stable because it is not expected to see any reduction in the CFRP modulus at −30C to 80C [45]. However, as shown in Figure 6(a), a significant decreasing trend is observed at temperatures greater than 45C. This decrease is mainly caused by the deterioration of the thermoplastic adhesive film used for bonding the sensors.

To obtain a reliable reference strain, the average value of the strain-load coefficient at temperatures below 45C is adopted to calibrate the measurement strain of Gauge2 above 45C. With the calibrated reference strain and SMF2 measurements near Gauge2 (49th segment shown in Figure 3(a)), the strain-frequency shift coefficient is directly calculated and shown in Figure 6(b). This coefficient is stable when the temperature stays between −20C and 55C. However, when the temperature becomes higher than 55C or lower than −20C, a significant decrease is observed. This decrease induces non-negligible errors in SMF strain even with temperature compensation when a constant strain-frequency coefficient is adopted. This conclusion is demonstrated by Figure 7(a) which shows a comparison of Gauge2 reference strain and the calculated strain from the SMF2 at 50C and 80C when the strain-frequency shift coefficient at 20C is adopted. Note that the data of the 49th segment on SMF2 (Figure 3(a)) are adopted in this figure. Based on this conclusion, a further analysis of SMF strain error is conducted. In this analysis, the frequency shift under the current state and the strain-frequency shift coefficient under the baseline state are adopted to calculate the SMF strain. The error between the calculated SMF strain and the reference strain is defined as as follows:where is the temperature of the current state; is the temperature of the baseline state; expresses the frequency shift of the SMF under the current state; expresses the strain-frequency coefficient for the baseline state; denotes the reference strain under the current state. The errors for different current and baseline state combinations are shown in Figure 7(b). In this figure, the x-axis is the temperature of the current state (the frequency shift data of SMF2 are adopted under this temperature); the y-axis shows the temperature of the baseline state (the strain-frequency shift coefficient is adopted under this temperature); the z-axis express the percent error of the calculated strain compared with the reference strain. For comparison with the reference strain, only the data of the 49th segment on SMF2 (Figure 3(a)) are adopted in this figure. This figure illustrates that increases rapidly when the temperature of the current or the baseline state is below −20C or above 55C. This result suggests that a constant strain-frequency shift coefficient cannot be adopted for strain measurement at all temperature conditions. However, there exists a temperature insensitive area according to Figure7(b). When the temperature of both current and baseline states are between −20C and 55C, is less than 2% and a constant strain-frequency shift coefficient can be adopted. Therefore, the calibration of the strain-frequency shift coefficient is only necessary when the measurement temperature is out of this range.

3.4. Numerical Verification

The SMF strain is calculated and analyzed in the previous sections. For example, the calculated strain along SMF2 under different loads at 20C is shown in Figure 8(a). To validate the reliability of the calculated strain, a numerical model is built based on the elastic parameters of the carbon fibre which is shown in Table 1. Figure 8(b) shows the model of the composite laminate. A 10 kN (625 MPa) tension load is applied to the short edge of the composite laminate. The simulation result is shown in Figure 8(c). The maximum displacement is , which means that the strain of the specimen under 10 KN (625 MPa) tension is .

As shown in Figure 8(a), the calculated strain of SMF2 under 10 kN (625 MPa) tension at 0 direction is (C), which is less than the numerical result. This is mainly because of the strain transfer rate between the substrate and SMF. In the work of [46], a model is built to estimate this strain transfer rate by equation (8). The expression of in this equation is shown in equation (9) as follows:where is the strain transfer rate; represents the sensing gauge length; expresses the radius of SMF; denotes bonding thickness; is the bonding length; is the shear modulus of composite laminate; is the elastic modulus of SMF; represents the influence depth of the shear stress; denotes the elastic modulus of the adhesive film; is the shear modulus of the adhesive film.

The parameters of the thermoplastic adhesive film and SMF adopted in this experiment are summarized in Table 2. With these parameters, the strain transfer rate is therefore obtained using equations (8) and (9), which was calculated to be 0.842 (set ). In this experiment, the strain transfer rate calculated from the modelling result and SMF measurement is 0.824. This result agrees fairly well with the value of the strain transfer rate calculated from the literature [46]. This result demonstrates that the calculated strain of SMF2 based on the proposed data cleaning and reconstruction, temperature compensation, and strain-frequency shift coefficient calibration is highly accurate and reliable.

4. Data Analyses of Fatigue Experiment

4.1. Dynamic Measurements and Data Resampling

In the fatigue experiment, the data of SMF2 and Gauge2 are measured under cyclic load with a frequency of 5 Hz and an amplitude of 3 kN (625 MPa). The SMF2 measurements at around 54000 cycles are shown in Figure 9(a). This diagram illustrates high-precision continuous sensing by the SMF under low-frequency vibration load. Figure 9(b) shows the strain data of Gauge2 at around 54000 cycles. As with the static measurements, SMF-based distributed sensing obtains more strain information on specimens than strain gauges.

The sampling frequency of SMF2 and Gauge2 is 50 Hz and 100 Hz, respectively. To compare the measurement results of SMF2 and Gauge2, the SMF2 measurements near Gauge2 are adopted (49th segment shown in Figure 3(a)). An interpolation method is then adopted to unify the sampling frequency. The results of SMF2 and Gauge2 are shown in Figures 10(a) and 10(b), respectively. The sampling frequency after the interpolation is 5 kHz.

4.2. Strain Measurement Errors under Low Frequency Vibration Load

As with the static analysis, the Gauge2 measurements are adopted as the reference strain in the strain-frequency shift coefficient calculation of SMF2. Therefore, the reliability of Gauge2 measurements throughout the fatigue experiment is first investigated based on the analysis of the strain-load coefficient. The calculation result of this coefficient is shown in Figure 11(a). This figure illustrates that the ratio between the strain shift and the corresponding load shift is stable during the whole fatigue experiment. In addition, the value of this coefficient is consistent with the result in the static loading experiment as shown in Figure 6(a). However, the reading of Gauge2 creeps with the increase of the fatigue cycle, even though the strain-load coefficient stays stable. This creep behaviour is illustrated in Figure 11(b) which shows the reading at 0 kN and 3 kN (187.5 MPa). In this figure, the reading of Gauge2 at 0 kN gradually creeps from to throughout the whole fatigue experiment. Therefore, the Gauge2 reading at the beginning of the fatigue experiment was adopted as the reference strain during the whole fatigue experiment to avoid this creep effect.

The strain-frequency coefficient is then obtained based on the reference strain and reading of the 49th segment on SMF2 (Figure 3(a)) at 3 kN (187.5 MPa). The result of this coefficient is shown in Figure 12(a). The strain-frequency coefficient stays stable and is consistent with the result in the static loading experiment shown in Figure 6(b) during the whole fatigue experiment. Therefore, the strain under 0 to 3 kN (0 to 187.5 MPa) during the fatigue experiment measured by SMF2 is calculated based on the frequency shift of SMF2 and the strain-frequency coefficient. The peak values of the calculated strains are shown in Figure 12(b). A constant load is induced during the fatigue experiment because the reference spectrum of the SMF is set before the experiment, which is the reason for nonzero values in Figure 12(b) when the tension load is 0 kN. These two diagrams demonstrate that strain measurement errors of the SMF-based distributed sensing under low-frequency vibration load is negligible when a constant strain-frequency shift coefficient is adopted. Although the fatigue experiments do not show a significant influence on the strain-frequency shift coefficient, there is a clear tendency that the proportion of null values in the SMF2 measurements increases with the increase of fatigue cycles. These null values in strain measurements are shown as a white area in Figure 13(a). The proportion of null value with an increase in the fatigue cycle is shown in Figure 13(b).

5. Conclusions

This study investigates the robustness of single-mode fibres, as strain sensors, for in-field applications, where operational and environmental conditions such as thermal fluctuations and vibration loads can influence the reliability of the strain distribution. Because SMF does not measure strains directly, there needs to be a conversion coefficient used for transferring the recorded frequency shifts to strain output. This paper investigates the influence of temperature and vibration (fatigue load) on the calculation of this conversion coefficient. To this aim, two experiments were conducted separately on a CFRP specimen under thermal and low-frequency vibration load. The strain measurement was compared with strain gauge readings, and the errors of the SMF-based distributed sensing were investigated by strain-frequency shift coefficient analyses under temperature fluctuations and mechanical vibrations. The conclusions of this work can be summarized as follows:(i)Temperature fluctuations induce non-negligible errors in strain measurement even with temperature compensation. This error is mainly caused by the mechanical deterioration of adhesives. However, there exists a temperature insensitive range in which this error is negligible. For example, the temperature insensitive range in this study is −20C to 55C. Calibration of the strain-frequency shift coefficient is necessary for high accurate strain measurements when the measurement temperature is out of this range for in-field applications.(ii)SMF-based distributed sensing obtained more strain information compared to strain gauges in both the static loading experiment and fatigue experiment. A reading creep of Gauge2 is observed even though the strain-load coefficient is stable during the entirety of the fatigue experiment. In contrast, the calculated strain and strain-frequency coefficient of SMF2 remained stable during the whole fatigue experiment, which means that strain measurement errors of SMF-based distributed sensing under low-frequency vibration load is negligible. This conclusion illustrates the advantage of SMF for accurate distributed sensing in low-frequency vibration environments compared to strain gauges.(iii)Although the fatigue experiments do not show a significant influence on SMF2 measurement values, there is a clear tendency that the number of null values in SMF measurements will increase as the fatigue cycles increase. Therefore, data reconstruction algorithms are required for recovering the strain information when SMF is adopted for strain measurement of composite laminate under low-frequency vibration load. It is worth noting that excessive data reconstruction may lead to the suppression of important features in the strain field, which can have adverse effects on various applications, including damage detection. Specifically, while the primary goal of the smoothing process is to restore the null values arising under low-frequency vibration load, it can also unintentionally dampen critical strain features near the damaged region, thereby hampering potential uses such as damage detection.

Data Availability

The MAT-File 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.

Authors’ Contributions

Yingwu Li was involved in conceptualization, software, methodology, visualization, investigation, and writing the original draft. Zahra Sharif-Khodaei was involved in conceptualization, supervision, methodology, resources, writing, reviewing, and editing.