# Stress Measurement of Steel Strands Based on System Identification Method

Authors: Ji Qian; Jianbin Li; Changchun Li; Jianting Zhou; Guowen Yao
Journal: Shock and Vibration
 (2019)

---

## Abstract

Steel strands are crucial components of a prestressed structure, but their performance degrades inevitably with time. Among the existing nondestructive testing methods for measuring the stress state of steel strands, ultrasonic guided wave methods have received the most research attention and are most likely to be applied to actual engineering. This work views the propagation of ultrasonic guided wave as an individual system, and the variations in the stress state in steel strands can be related to the parameters gained by the system identification model. The finite element simulation and experimental results show the propagation characteristics of the guided wave are influenced distinctly by the stress state. The index constructed by the system parameters exhibits good monotonic linear with stress level variation, and the slope of the fitted line from experimental data is similar to the simulation result. The sensor placement and loading path have a limited effect on the stress state and the proposed index.


---

## Body

## 1. Introduction

Prestressed steel strands are used widely in civil engineering buildings and are the most important mechanical components of long-span bridges, acting as the main prestressed tendons of girder bridges, suspenders of arch bridges, and cables of cable-stayed bridges. Being in the high-stress state for a long time, the steel strands are very sensitive to environmental erosion and structural creep and the actual retained working stress always decreases, resulting in structural cracking, deflection, and other failures, which finally reduces the bearing capacity and durability of long-span structures. Steel strands are usually protected with concrete or sheaths of high-density polyethylene, which are used as anticorrosion requirements. However, while the anticorrosion performance is enhanced, inspecting and monitoring the steel strands during periodical inspections can be difficult. The detection and evaluation of the stress of in-service steel strands have long been a technical problem in the field of civil engineering. The relatively mature technology for measuring the stress of steel strands focuses on the detection of stress increment and should adopt embedments in its early stage, making it difficult to be applied to in-service prestressed and cable-supported bridges [1].The use of ultrasonic guided wave (UGW) methods for monitoring the stress status of steel strands has elicited considerable attention in recent years. The change in guided wave velocity has been applied to monitor stress in plates with biaxial loads [2], bolts (to verify tightening) [3], rails [4–6], steel strands [7–9], grouted tendons [10], and pipes [11]. Kwun et al. [12] performed early experimental studies on the effect of stress on guided waves in a strand. They excited guided waves as a whole throughout the cross section of a strand using magnetostrictive sensors. Chen et al. [13] conducted an early study on the fundamental longitudinal guided model (0, 1) excited within the individual wires of a strand.Most studies on stress estimation in the experimental and theoretical aspects are based on the acoustoelastic effect. Group velocity in low-order modes will decrease with the increase in the stress of steel strands [14]. However, conducting acoustoelastic measurements is difficult in this case because the acoustoelastic effect is inherently low. The velocity of the first-order longitudinal wave only exhibits 5% variation when axial stress reaches its maximum value [15]. The acoustoelastic effect also presents a nonlinear behavior, particularly under the conditions of assembled strands and low stress [16]. Di Scalea et al. [17] proposed a similar description that was formulated under the same test conditions used in Washer [16]. However, a sound explanation for these irregularities has yet to be provided.Another phenomenon in which a portion of the frequency band will disappear from the guided wave frequency spectrum with increasing tensile stress has been discussed. Kwun et al. [12] and Wu et al. [18, 19] found from their experiment results that the longitudinal modes disappear in certain bands and the center frequency (notch frequency) of the band is logarithmically linear to tension when a steel strand is loaded with an axial force. This phenomenon was also reported in Laguerre et al. [20] when guided wave propagation was compared between stressed and unstressed seven-wire strands.A semianalytical finite element (SAFE) method has been used recently to study guided waves propagating as a whole within seven-wire strands. Treyssède and Laguerre [21] presented the first SAFE study applied to seven-wire strands. They explained the notch frequency discovered by Kwun et al. [12]. Nucera and Di Scalea [22] studied the effect of stress on higher harmonics in which the generation of a second harmonic from a primary excitation was correlated with stress. Schaal et al. [23] and Treyssède [24] examined energy leakage between adjacent wires in a stressed strand. SAFE was adopted as an effective approximation method for dealing with a regular waveguide. Accordingly, the use of enormous computing resources in a minuscule meshing size was avoided.The characteristics of guided waves are affected undoubtedly by the stress status of steel strands. However, convincing explanations for these phenomena remain lacking. In this work, seven-core steel strands commonly used in structures and bridges are taken as subjects. Having considered the propagation of the guided wave in steel strands to be an individual system, the parameters of this system are recognized by the system identification method and a relationship between the proposed index composed by model parameters and stress level is constructed through finite element simulation and experiments. Finally, this work discusses the influence of the placement of sensors and loading paths.


## 2. System Identification Model and Stress Index

The propagation of elastic waves in solids is an elastodynamic problem, and the dispersion curve can be solved by incorporating the geometrical and mechanical boundary conditions into the wave equation. For single high-strength steel wire, the classical theory of the cylindrical waveguide can explain the longitudinal waves, bending waves, torsional waves, and other modal curves. However, establishing the mathematical model of guided wave propagation in strands is difficult because of the contact coupling problem of multiple steel wires and no analytical solution exists at present.When the guided wave propagating within the steel strand is regarded as an individual system, the change of stress status can lead to the variation of system performance (Figure1). The stress level can be evaluated by analyzing the parameters of the system model. The measured guided wave is a weakly stationary process, and it can be described by an auto-regressive model in time series.Figure 1
System of guided wave propagation.The auto-regressive model is given by(1)xt=q1xt−1+q2xt−2+⋯+qpxt−p+et,where xt is the sample observation, et is a random noise with zero mean, and qjj=1,2,…,p are regression coefficients.When a time seriesxt is plugged into equation (1), the results can be expressed as follows:(2)xn+1=q1x+nq2x+n−1⋯+qpx1−en+1,xn+2=q1x+n+1q2x+n⋯+qpx2−en+2,⋮xN=q1x+N−1q2x+N−2⋯+qpxN−n−eN.Then, equation (2) can be written into matrix form(3)Y=Xφ−ε,where(4)Y=xn+1xn+2⋯xNT,φ=q1q2⋯qpT,ε=en+1en+2⋯eNT,X=xnxn−1⋯x1xn+1xn⋯x2⋮⋮⋮⋮xN−1xN−2⋯xN−n.The sum of squares of residuals is taken as the criterion function:(5)Jφ=∑i=1Nε2=∑i=1NY−Xφ2=YTY+φTXTXφ−φTXTY−YTXφ.The necessary condition for finding the least-squares estimateφ is(6)∂Jφ/∂φφ=φls=0.By substituting equation (5) in (6), the following can be obtained:(7)XTXφls=XTY.IfXTX is full in rank (7), then XTX−1 exists, and it can be derived as follows:(8)φls=XTX−1XTY.Theφls is the least-squares estimate of the parameter qjj=1,2,…,p.According to the least-squares identification, the most probable value of the unknown model parameterφ is obtained when the sum of the cumulative errors between the actual and calculated values reaches the minimum.The process of measuring the stress status of steel strands based on the auto-regressive model is to extract the auto-regressive parameters of the system model as feature vectors and to calculate the distance between different stress levels. While the model in stress level 70%fpk (fpk is the ultimate tensile strength) is considered as the reference, the index Istf is defined as the distance between a certain stress level to reference:(9)Istfi=∑j=1pqjref−qji2,where i is a certain stress level, qjref is the parameter in the reference model, and qji is the parameter at a certain stress level.


## 3. Finite Element Simulation

### 3.1. Finite Element Model

ABAQUS/Explicit is used to simulate the propagation of ultrasonic guided waves in steel strands. The length of the steel strands is 0.52 m. The geometric and material parameters are shown in Table1.Table 1
Geometric and material parameters of steel strands.
Geometric parameter
Material parameter
Steel wire diameter,dc (mm)
5.08
Young’s modulus,E (GPa)
196
Strand diameter,d (mm)
15.20
Poisson’s ratio,v
0.29
Peripheral pitch,h (mm)
260
Density,ρ (kg/m3)
7850
Peripheral twist angle,β∘
7.90
Ultimate tensile strength,fpk (MPa)
1860The finite element model needs eight nodes in a wavelength [25] and a maximum frequency of 500 kHz is considered in this paper to capture accurately the fluctuation effect which indicates Δl≤λmin/8−1=cT/7×fmax=3222.6/7×500000=0.92mm, of which cT is the velocity of transverse waves. The unit size along the axis of the strands is 1 mm (slightly larger than 0.92 mm) and the minimum size of the contact area is 0.1 mm. The model is composed of 1,745,623 8-node hexahedral elements after meshing (Figure 2). Apart from the dimension of grids, time integral step is another important factor in controlling the accuracy of finite element calculation. The response of the structure can be regarded as a combination of modes of each order, and the minimum integral step should be able to solve the highest order modes in the combination of structural response modes. In the process of solving transient dynamics with the finite element method, the time integration step size is generally less than 1/20 of the maximum frequency and should meet the requirements of computing stability at the same time with all free time integral step.Figure 2
Finite element calculation model of steel strand.
(a)
(b)Normal contact between steel wires in strands is simulated with “hard” contact and the tangential contact is simulated with “penalty” friction formula with a friction coefficient of 0.6 [22]. One end of the strand is fixed completely and the other only releases axial displacement, thereby restricting the other direction of freedom.The entire process of simulation is divided into three stages: axial tension force infliction, guided wave excitation, and guided wave propagation.Stage 1: tension is imposed at one end of the axial displacement, which is a quasi-static loading process. The amplitude curve of tension should be as smooth as possible and the applied time should be as long as possible to prevent interference signal. In this paper, the smooth amplitude curve is used to load and the loading time is 300μs. The loading amplitude curve is shown in Figure 3.Stage 2: guided wave excitation is carried out at the center of the central steel wire. The excitation is a triangular pulse with a duration of 3μs, as shown in Figure 4. The energy of the guided wave generated by the excitation should be much larger than that of the prestressed one to prevent the guided wave signal caused by the excitation load from being submerged by the disturbance signal generated during the prestressed loading process.Stage 3: after exerting the excitation pulse, the propagation of elastic waves is simulated in the steel strands with 697μs.Figure 3
Tension load curve of tension.Figure 4
Excitation pulse.


### 3.2. Parameters Extraction of Auto-Regression Model

The received waves under different stress level are shown in Figure5.Figure 5
Finite element simulation of the guided wave signal.The received guided waves are used to establish the auto-regression model after the mean and trend terms are removed. The next important step is finding an appropriate order of the model to ensure accuracy. The Akaike information criterion (AIC) assumes that when the density function of a parameter model withkk≤p parameters is denoted as gyθk and the maximum likelihood function is denoted gθ^ky, where θk is an unknown parameter and θ^k is its maximum likelihood estimation, the AIC is defined as AIC=−2lngθ^ky+2k. The first term reflects the goodness of model fit, while the second term represents a “penalty” for the model with more variables. Finally, the final order is the result of weighing the effects of model fitting and the number of variables.The optimal order of the system is determined asp=10 when the AIC is applied to handle the data in this work. Taking the triangular pulse in Figure 4 as the excitation and the received signal under 20%–70% fpk as the output of the system, the parameters of the auto-regression model are obtained and shown in Table 2.The increase of tension in steel strands causes the parametersq1⋯q10 of the system model to change significantly and monotonously. In order to verify the accuracy of the identification model of finite element data, the fitting degree is calculated according to te following equation:(10)fit=1−Q∑i=1ny20.5,where Q=∑i=1ny−y∗2 is the residual sum of squares, in which y is the original data and y∗ is calculated by the system model.The waveform shown in Figure3 is taken as the input of the system, and q1⋯q10 corresponding to 70% fpk is taken as the model parameter. The original data and predicted data are shown in Figure 6(a), and the residual error is shown in Figure 6(b). The relative value of residual errors is less than 2% and the fitting degree reached 98.41% as calculated by (10), which shows that the auto-regression model can obtain a reasonable order and reliable recognition results.Table 2
Identification model parameters in different tension force.
Stress
Parameter
q
1
q
2
q
3
q
4
q
5
q
6
q
7
q
8
q
9
q
10
0.2fpk
−2.525
1.977
−0.042
−0.476
−0.171
0.521
−0.316
0.068
−0.057
0.060
0.3fpk
−2.491
1.955
−0.158
−0.364
−0.075
0.348
−0.205
0.030
−0.083
0.082
0.4fpk
−2.518
2.083
−0.358
−0.249
−0.088
0.278
−0.056
−0.046
−0.118
0.110
0.5fpk
−2.558
2.22
−0.493
−0.286
0.057
0.181
−0.043
0.005
−0.170
0.124
0.6fpk
−2.653
2.546
−0.852
−0.226
0.241
−0.123
0.346
−0.332
0.016
0.069
0.7fpk
−2.656
2.652
−1.111
−0.013
0.223
−0.161
0.261
−0.143
−0.124
0.106Figure 6
Fit degree of the recognition model.
(a)
(b)


### 3.3. Relationship between Stress Level and Index

The parameters of the auto-regression model were calculated by using the optimized order in stepwise loading procedure. The model parameters under the stress level of 70%fpk are taken as the reference and the index Istf is calculated according to (9). The results are shown in Figure 7.Figure 7
Relationship curves between the recognition index and the tensile force.The index of stress evaluationIstf shows an obvious monotonous variation with the increase in tensile stress. Considering the actual stress state of the structure during the actual operation period, the data do not include those of the stress state under 20% fpk. The linear fitting of the data is judged by the coefficient R2=1−∑i=1nyi−fxi2/∑i=1nyi−y¯2, where fxi is the model predicted stress value, yi is the actual stress value, and the average value is y¯. The calculated slope of the fitting line is −5.85 and the definite coefficient is R2=0.9738, which presents a clear linear rule.Taking the recognition indexIstf as an independent variable and the tension force of strands as a dependent variable, the functional relationship between the stress level σ and the index Istf is established as(11)σfpk=−0.171Istf+0.662.


## 3.1. Finite Element Model

ABAQUS/Explicit is used to simulate the propagation of ultrasonic guided waves in steel strands. The length of the steel strands is 0.52 m. The geometric and material parameters are shown in Table1.Table 1
Geometric and material parameters of steel strands.
Geometric parameter
Material parameter
Steel wire diameter,dc (mm)
5.08
Young’s modulus,E (GPa)
196
Strand diameter,d (mm)
15.20
Poisson’s ratio,v
0.29
Peripheral pitch,h (mm)
260
Density,ρ (kg/m3)
7850
Peripheral twist angle,β∘
7.90
Ultimate tensile strength,fpk (MPa)
1860The finite element model needs eight nodes in a wavelength [25] and a maximum frequency of 500 kHz is considered in this paper to capture accurately the fluctuation effect which indicates Δl≤λmin/8−1=cT/7×fmax=3222.6/7×500000=0.92mm, of which cT is the velocity of transverse waves. The unit size along the axis of the strands is 1 mm (slightly larger than 0.92 mm) and the minimum size of the contact area is 0.1 mm. The model is composed of 1,745,623 8-node hexahedral elements after meshing (Figure 2). Apart from the dimension of grids, time integral step is another important factor in controlling the accuracy of finite element calculation. The response of the structure can be regarded as a combination of modes of each order, and the minimum integral step should be able to solve the highest order modes in the combination of structural response modes. In the process of solving transient dynamics with the finite element method, the time integration step size is generally less than 1/20 of the maximum frequency and should meet the requirements of computing stability at the same time with all free time integral step.Figure 2
Finite element calculation model of steel strand.
(a)
(b)Normal contact between steel wires in strands is simulated with “hard” contact and the tangential contact is simulated with “penalty” friction formula with a friction coefficient of 0.6 [22]. One end of the strand is fixed completely and the other only releases axial displacement, thereby restricting the other direction of freedom.The entire process of simulation is divided into three stages: axial tension force infliction, guided wave excitation, and guided wave propagation.Stage 1: tension is imposed at one end of the axial displacement, which is a quasi-static loading process. The amplitude curve of tension should be as smooth as possible and the applied time should be as long as possible to prevent interference signal. In this paper, the smooth amplitude curve is used to load and the loading time is 300μs. The loading amplitude curve is shown in Figure 3.Stage 2: guided wave excitation is carried out at the center of the central steel wire. The excitation is a triangular pulse with a duration of 3μs, as shown in Figure 4. The energy of the guided wave generated by the excitation should be much larger than that of the prestressed one to prevent the guided wave signal caused by the excitation load from being submerged by the disturbance signal generated during the prestressed loading process.Stage 3: after exerting the excitation pulse, the propagation of elastic waves is simulated in the steel strands with 697μs.Figure 3
Tension load curve of tension.Figure 4
Excitation pulse.


## 3.2. Parameters Extraction of Auto-Regression Model

The received waves under different stress level are shown in Figure5.Figure 5
Finite element simulation of the guided wave signal.The received guided waves are used to establish the auto-regression model after the mean and trend terms are removed. The next important step is finding an appropriate order of the model to ensure accuracy. The Akaike information criterion (AIC) assumes that when the density function of a parameter model withkk≤p parameters is denoted as gyθk and the maximum likelihood function is denoted gθ^ky, where θk is an unknown parameter and θ^k is its maximum likelihood estimation, the AIC is defined as AIC=−2lngθ^ky+2k. The first term reflects the goodness of model fit, while the second term represents a “penalty” for the model with more variables. Finally, the final order is the result of weighing the effects of model fitting and the number of variables.The optimal order of the system is determined asp=10 when the AIC is applied to handle the data in this work. Taking the triangular pulse in Figure 4 as the excitation and the received signal under 20%–70% fpk as the output of the system, the parameters of the auto-regression model are obtained and shown in Table 2.The increase of tension in steel strands causes the parametersq1⋯q10 of the system model to change significantly and monotonously. In order to verify the accuracy of the identification model of finite element data, the fitting degree is calculated according to te following equation:(10)fit=1−Q∑i=1ny20.5,where Q=∑i=1ny−y∗2 is the residual sum of squares, in which y is the original data and y∗ is calculated by the system model.The waveform shown in Figure3 is taken as the input of the system, and q1⋯q10 corresponding to 70% fpk is taken as the model parameter. The original data and predicted data are shown in Figure 6(a), and the residual error is shown in Figure 6(b). The relative value of residual errors is less than 2% and the fitting degree reached 98.41% as calculated by (10), which shows that the auto-regression model can obtain a reasonable order and reliable recognition results.Table 2
Identification model parameters in different tension force.
Stress
Parameter
q
1
q
2
q
3
q
4
q
5
q
6
q
7
q
8
q
9
q
10
0.2fpk
−2.525
1.977
−0.042
−0.476
−0.171
0.521
−0.316
0.068
−0.057
0.060
0.3fpk
−2.491
1.955
−0.158
−0.364
−0.075
0.348
−0.205
0.030
−0.083
0.082
0.4fpk
−2.518
2.083
−0.358
−0.249
−0.088
0.278
−0.056
−0.046
−0.118
0.110
0.5fpk
−2.558
2.22
−0.493
−0.286
0.057
0.181
−0.043
0.005
−0.170
0.124
0.6fpk
−2.653
2.546
−0.852
−0.226
0.241
−0.123
0.346
−0.332
0.016
0.069
0.7fpk
−2.656
2.652
−1.111
−0.013
0.223
−0.161
0.261
−0.143
−0.124
0.106Figure 6
Fit degree of the recognition model.
(a)
(b)


## 3.3. Relationship between Stress Level and Index

The parameters of the auto-regression model were calculated by using the optimized order in stepwise loading procedure. The model parameters under the stress level of 70%fpk are taken as the reference and the index Istf is calculated according to (9). The results are shown in Figure 7.Figure 7
Relationship curves between the recognition index and the tensile force.The index of stress evaluationIstf shows an obvious monotonous variation with the increase in tensile stress. Considering the actual stress state of the structure during the actual operation period, the data do not include those of the stress state under 20% fpk. The linear fitting of the data is judged by the coefficient R2=1−∑i=1nyi−fxi2/∑i=1nyi−y¯2, where fxi is the model predicted stress value, yi is the actual stress value, and the average value is y¯. The calculated slope of the fitting line is −5.85 and the definite coefficient is R2=0.9738, which presents a clear linear rule.Taking the recognition indexIstf as an independent variable and the tension force of strands as a dependent variable, the functional relationship between the stress level σ and the index Istf is established as(11)σfpk=−0.171Istf+0.662.


## 4. Experiments

### 4.1. Experiments Equipment

The experiment system includes the reaction wall, hydraulic jack, anchor devices, guided wave apparatus, preamplifier, and sensors (Figure8). The loading step controlled by a background terminal was set as 10% of the ultimate tensile strength (fpk) and the entire loading process was from 20% to 70% fpk (52 kN–182 kN). Every step was sustained for 5 min to excite and collect guided waves. The experiments were performed on eight steel strands.Figure 8
Design of the experiment system (unit: mm). (a) Layout of the experiment system. (b) Equipment of the system.
(a)
(b)The material and geometric parameters of steel strands are shown in Table1. The end faces of the steel strands were polished to ensure good contact between the steel strands and the sensors. The sensors used in this experiment have a wideband frequency response in the range of 100–1000 kHz, and the sampling frequency was 2 MHz.The system excitation is a series of single-period sinusoidal pulses with step frequencyΔf = 50 kHz and is shown in the following equations: (12)ft=∑i=0f−f0/ΔfVitsin2πf0+i⋅Δf,(13)Vit=10,∑j=0iTj≤t<∑j=0iTj+1,0,else,where f0 = 100 kHz is the origin frequency, f = 1000 kHz is the cut-off frequency, Vit is a rectangular window function, Tj is the period, and T0 is zero.


### 4.2. Analysis of Experimental Results

The input excitation and measured guide wave under different stress levels are shown in Figures9 and 10. The excitation waveform remains constant, but the amplitude of the measured waveform decreased obviously with the increase in the stress level. This condition may be attributed to wires having more contact surfaces as stress is enhanced. Accordingly, more guided waves propagate back and forth among wires, causing the energy to be attenuated quickly.Figure 9
Input excitation.Figure 10
Measured guided waves.Using the mentioned modeling methods above, the optimal model order of the measured signal is determined to bep=10. Eight strands are loaded separately stepwise in experiments, and the same model order is used to calculate the parameters q1⋯q10. The average value q¯1⋯q¯10 of eight strands at 70% fpk is taken as the reference value, and the index Istf is calculated according to equation (9). The results are shown in Figure 11, and the parameters of the fitted line are shown in Table 3.Figure 11
Relationship curves between the recognition index and tensile force.Table 3
The parameters of the fitted line.
Number of strands
1
2
3
4
5
6
7
8
Slope
−5.765
−6.141
−5.749
−6.462
−5.738
−5.829
−5.699
−5.767
Intercept
3.920
4.093
4.087
4.389
3.846
3.886
3.921
3.920The indexIstf calculated using the eight independent loading procedures had evident monotonous variation with the increase of the stress level. Using the linear fit, the slope of the fitted line is −5.81 and the determination coefficient is 0.9585. The experimental results indicate good consistency, reproducibility, and stability.Taking the indexIstf as an independent variable and the stress level as a dependent variable, the functional relationship between them is established as follows:(14)σfpk=−0.172Istf+0.671.Although the slopes of the fitted line have little fluctuation, the line fitted by all experimental data collectively is similar to the finite element result. Thus, the system identification model reveals distinctly the propagation performance of the steel strand and that the parameters of the system can be affected by the corresponding stress status.


## 4.1. Experiments Equipment

The experiment system includes the reaction wall, hydraulic jack, anchor devices, guided wave apparatus, preamplifier, and sensors (Figure8). The loading step controlled by a background terminal was set as 10% of the ultimate tensile strength (fpk) and the entire loading process was from 20% to 70% fpk (52 kN–182 kN). Every step was sustained for 5 min to excite and collect guided waves. The experiments were performed on eight steel strands.Figure 8
Design of the experiment system (unit: mm). (a) Layout of the experiment system. (b) Equipment of the system.
(a)
(b)The material and geometric parameters of steel strands are shown in Table1. The end faces of the steel strands were polished to ensure good contact between the steel strands and the sensors. The sensors used in this experiment have a wideband frequency response in the range of 100–1000 kHz, and the sampling frequency was 2 MHz.The system excitation is a series of single-period sinusoidal pulses with step frequencyΔf = 50 kHz and is shown in the following equations: (12)ft=∑i=0f−f0/ΔfVitsin2πf0+i⋅Δf,(13)Vit=10,∑j=0iTj≤t<∑j=0iTj+1,0,else,where f0 = 100 kHz is the origin frequency, f = 1000 kHz is the cut-off frequency, Vit is a rectangular window function, Tj is the period, and T0 is zero.


## 4.2. Analysis of Experimental Results

The input excitation and measured guide wave under different stress levels are shown in Figures9 and 10. The excitation waveform remains constant, but the amplitude of the measured waveform decreased obviously with the increase in the stress level. This condition may be attributed to wires having more contact surfaces as stress is enhanced. Accordingly, more guided waves propagate back and forth among wires, causing the energy to be attenuated quickly.Figure 9
Input excitation.Figure 10
Measured guided waves.Using the mentioned modeling methods above, the optimal model order of the measured signal is determined to bep=10. Eight strands are loaded separately stepwise in experiments, and the same model order is used to calculate the parameters q1⋯q10. The average value q¯1⋯q¯10 of eight strands at 70% fpk is taken as the reference value, and the index Istf is calculated according to equation (9). The results are shown in Figure 11, and the parameters of the fitted line are shown in Table 3.Figure 11
Relationship curves between the recognition index and tensile force.Table 3
The parameters of the fitted line.
Number of strands
1
2
3
4
5
6
7
8
Slope
−5.765
−6.141
−5.749
−6.462
−5.738
−5.829
−5.699
−5.767
Intercept
3.920
4.093
4.087
4.389
3.846
3.886
3.921
3.920The indexIstf calculated using the eight independent loading procedures had evident monotonous variation with the increase of the stress level. Using the linear fit, the slope of the fitted line is −5.81 and the determination coefficient is 0.9585. The experimental results indicate good consistency, reproducibility, and stability.Taking the indexIstf as an independent variable and the stress level as a dependent variable, the functional relationship between them is established as follows:(14)σfpk=−0.172Istf+0.671.Although the slopes of the fitted line have little fluctuation, the line fitted by all experimental data collectively is similar to the finite element result. Thus, the system identification model reveals distinctly the propagation performance of the steel strand and that the parameters of the system can be affected by the corresponding stress status.


## 5. Discussion

### 5.1. Positions of Sensor Placement

Whether the sensor is arranged at the end section or the side surface, the measured guided wave has different modes, such as longitudinal wave, and bending wave, or the energy proportion among different modes changes. Under the same load condition and excitation input, the indexIstf calculated from the sensors at the end and side positions are shown in Figure 12.Figure 12
Recognition index curve at different sensors position.In Figure12, index Istf shows the same monotonic variation, regardless of sensors placement. Compared with the sensor arranged on the side surface, the absolute slope of the fitting line at the end section is very close but the determination coefficient is slightly larger. Hence, regardless of whether the sensors are placed at the end or on the side of steel strands, the index can reflect the stress state with a good linear law; however, the linear fit has a better effect in the former case.


### 5.2. Effects of Stress Loading Path

Evaluation of the bearing capacity of the prestressed structure only considers the stress state at some point, which has nothing to do with the loading path. Therefore, the stress index should avoid the effect of stress loading path as far as possible. The index is calculated on the condition of loading and unloading stage, which is shown in Figure13. The parameters of the fitting line are shown in Table 4.Figure 13
Recognition index curve for loading and unloading process.Table 4
Parameters of the fitting line.
k
R
2
First loading
−5.93
0.9842
First unloading
−5.34
0.9861
Second loading
−5.85
0.9713
Second unloading
−5.29
0.9798Figure13 shows that the index Istf in the process of loading or unloading is a monotonous linear change. Comparing the change of slope in two loading and unloading cycles, the difference between the two loading stages is 1.3% while that between two unloading stages is 0.9%. Relatively speaking, the slope between the loading and unloading stages change significantly, showing a decrease of 9.9% and 9.6%, respectively, mainly due to the reading errors caused by the hydraulic lag in load-unload switching.


## 5.1. Positions of Sensor Placement

Whether the sensor is arranged at the end section or the side surface, the measured guided wave has different modes, such as longitudinal wave, and bending wave, or the energy proportion among different modes changes. Under the same load condition and excitation input, the indexIstf calculated from the sensors at the end and side positions are shown in Figure 12.Figure 12
Recognition index curve at different sensors position.In Figure12, index Istf shows the same monotonic variation, regardless of sensors placement. Compared with the sensor arranged on the side surface, the absolute slope of the fitting line at the end section is very close but the determination coefficient is slightly larger. Hence, regardless of whether the sensors are placed at the end or on the side of steel strands, the index can reflect the stress state with a good linear law; however, the linear fit has a better effect in the former case.


## 5.2. Effects of Stress Loading Path

Evaluation of the bearing capacity of the prestressed structure only considers the stress state at some point, which has nothing to do with the loading path. Therefore, the stress index should avoid the effect of stress loading path as far as possible. The index is calculated on the condition of loading and unloading stage, which is shown in Figure13. The parameters of the fitting line are shown in Table 4.Figure 13
Recognition index curve for loading and unloading process.Table 4
Parameters of the fitting line.
k
R
2
First loading
−5.93
0.9842
First unloading
−5.34
0.9861
Second loading
−5.85
0.9713
Second unloading
−5.29
0.9798Figure13 shows that the index Istf in the process of loading or unloading is a monotonous linear change. Comparing the change of slope in two loading and unloading cycles, the difference between the two loading stages is 1.3% while that between two unloading stages is 0.9%. Relatively speaking, the slope between the loading and unloading stages change significantly, showing a decrease of 9.9% and 9.6%, respectively, mainly due to the reading errors caused by the hydraulic lag in load-unload switching.


## 6. Conclusions

Stress measuring of steel strands based on the ultrasonic guided wave is conducted mostly using the theory of cylindrical waveguide, but the difference between those two waveguides is obvious, such as the contact condition among wires and energy leaks at the contact area, which enhances the difficulty in establishing a theoretical resolution in steel strands. In this work, the propagation of the guided wave in steel strands is considered as an individual system and the parameters of this system are extracted using the system identification model. These parameters are used to propose a stress index.The propagation of the guided wave in steel strands is simulated by the finite element method, and the auto-regressive model in time series is adopted to process the collected signal. A functional relation between the stress state and identification index, which presents a good linear rule, is constructed by optimizing the model order and verifying the accuracy of the system parameters.The modeling method used in simulation analysis is adopted to process the data from stepwise loading experiments in eight steel strands. Although the slope of the fitted lines displays little fluctuation, the line fitted by all experimental data collectively is close to the simulation result.The index can reflect the stress state with a good linear law regardless of whether the sensor is placed at the end or on the side of the steel strands; however, the linear fit has a better effect in the former. The slope of the fitted line remains unchanged under the condition of loading and unloading two cycles, which helps in eliminating the effect of loading path.



---

*Source: 5653902-2019-12-20.xml*