The effect of temperature variation on vibration monitoring of prestressed concrete (PSC) girders is experimentally analyzed. Firstly, vibration features such as autoregressive (AR) coefficient, correlation coefficient of power spectral density (CC of PSD), natural frequency, and mode shape are selected to estimate the effect of temperature variation on vibration characteristics of PSC girders. Secondly, vibration experiments on a lab-scale PSC girder are performed under the condition of temperature variation. Finally, the vibration features with respect to the temperature variation are analyzed to estimate the effect of temperature in vibration characteristics of the PSC girder.
1. Introduction
For the past decades, the interest in structural health monitoring of prestressed concrete (PSC) structures has been increased. In a variety of civil engineering structures, PSC girders are main components to resist against external loadings [1]. For a PSC girder bridge, flexural stiffness in girder and prestress force in tendon are two essential parameters that should be secured for its serviceability and safety against external loadings. Unless the PSC girder bridges are instrumented at the time of construction, the occurrence of damage cannot be directly monitored and other alternative methods should be sought.
Since as early as 1970s, many researchers have focused on the possibility of using vibration characteristics of a structure as an indication of its structural damage [2–6]. The most appealing feature associated with using vibration properties is that they are relatively simple to measure and to utilize for a prompt diagnosis. Recently, research efforts have been made to monitor the change in modal properties of the PSC structures in relation to the change in prestress forces [7], to investigate the dynamic behaviors of prestressed composite girder bridges [8], and to identify the change in prestress forces by measuring dynamic responses of prestressed beams [9–11].
However, the practicality of the so-called “bridge diagnosis via vibration monitoring” is limited for real PSC girder bridges since even significant damage sources may not be revealed as remarkable changes in vibration features due to temperature-induced uncertainty [12–17]. In a real-life situation of a PSC girder bridge, it is reasonable to say that temperature differences are about 10~20°C during day and night and about 20~50°C during a year. To account for the significant condition, therefore, temperature-induced changes in vibration properties should be accounted for the PSC girder bridge.
This study presents the effect of temperature on vibration responses of the prestressed concrete girder and also estimates the relative sensitivities of the selected vibration features. The following approaches are implemented to achieve the objective. Firstly, vibration features such as autoregressive coefficient, power spectral density, natural frequency, and mode shape are selected for estimating the effect of temperature variation on PSC girders. Secondly, experiments on a lab-scale PSC girder are described. The experiments are performed under the condition of temperature variation. Finally, vibration characteristics of the PSC girder are analyzed with respect to the temperature variation. Temperature-induced effects on vibration monitoring of the PSC girder are estimated for the selected vibration features.
2. Vibration Features for Estimation of Temperature Effect
Vibration features selected to estimate the effect of temperature variation on the PSC girder include autoregressive model (AR model), correlation coefficient of power spectral density (CC of PSD), natural frequency, and mode shape. Once acceleration responses are measured at distributed location, the vibration features are extracted for the structure. For each type of vibration features, the effect of temperature is estimated by regression analysis.
2.1. Autoregressive Model
AR model is one of the time series analysis methods. It forecasts future responses from past time history responses. In this study, the AR model is defined in terms of acceleration responses as follows [18]:(1)v¨t=∑j=1aφjv¨(t-j)+et,where v¨(t) is the time history of acceleration at time step t, φj is the jth AR coefficient, a is the order of AR model, and e(t) is the residual error. The change in AR coefficients represents the change in structural parameters.
Mahalanobis squared distance (MSD) is utilized to calculate the change of AR coefficients due to structural change. Consider(2)Dζ=φζ-φ-TS-1φζ-φ-,where Dζ is the MSD value, φζ is the AR coefficient vector (outlier), φ- is the mean vector of n sets of AR coefficient, and S is the covariance matrix of n sets of AR coefficient. The variation of MSD indicates the variation of AR coefficients due to structural change. The number of set, n, is determined by partial autocorrelation function of the first acceleration sample in intact state. It is noted that the number of set should be larger than the order of AR model.
2.2. Correlation Coefficient of Power Spectral Density
Assume there are two acceleration signals, xt and yt, measured before and after structural change, respectively. The corresponding power spectral densities, Sxx and Syy, are calculated from Welch’s procedure as [18](3)Sxxf=1ndT∑i=1ndXif,T2,(4)Syyf=1ndT∑i=1ndYif,T2,where X and Y are the frequency response transformed from correspondent acceleration signal; nd is the number of divided segments; and T is the data length of divided segment.
The correlation coefficient (CC) of PSDs represents the linear identity between the two PSDs obtained before and after structural change, as follows: (5)ρXY=ESxxSyy-μSxxμSyyσSxxσSyy,where E[·] is the expectation operator and σSxx and σSyy are the standard deviations of PSDs of acceleration signals measured before and after structural change, respectively. If any structural change occurs in target structure, its acceleration responses would be affected and, consequently, the decrement of CC can be a warning sign for the presence of the special cause in the structure.
For estimating CC of PSDs due to temperature variation, the excitation force should be consistently maintained during experiments. In the lab-scale experiment, the impact testing can be easily maintained to be consistent. However, the excitation magnitude can be different for each testing in the full-scale civil engineering structures. For such cases, the vibration responses should be normalized according to the excitation magnitude. In addition, the ambient vibration data seems to be difficult to apply as the input excitation varies over time, which results in different CC of PSDs although the temperature remains the same.
2.3. Modal Parameters
Frequency domain decomposition (FDD) method [19, 20] is used to extract modal parameters such as natural frequency and mode shape. The singular values of the PSD function matrix S(ω) are used to estimate the natural frequencies instead of the PSD functions themselves as follows [20]:(6)Sω=UωTΣωVω,where Σ is the diagonal matrix consisting of the singular values, σi’s (i=1,…,N), and U and V are unitary matrices. Since S(ω) is symmetric, U becomes equal to V. In this FDD method, the natural frequencies can be determined from the peak frequencies of the singular value and the mode shape from anyone of the column vectors of U(ω) at the corresponding peak frequencies. In this study, the first singular value σ1(ω) is used to estimate the modal parameters.
Mode shapes measured before and after structural change are examined by modal assurance criterion (MAC). The MAC provides a measure of the least squares deviation or scatter of the points from the linear correlation (7)MACϕi,ϕi∗=[ϕiTϕi∗]2ϕiTϕi[ϕi∗Tϕi∗],where ϕi and ϕi∗ are the modal vectors which are extracted from acceleration signals measured before and after structural change, respectively. The MAC value provides quantifying of the comparison between the two sets of mode shape data. In practice, a value of MAC close to the unity is expected if the two sets of mode shapes have no differences caused by structural change; otherwise, the MAC value would decrease. This method can be utilized to obtain the information on what vibration modes are more sensitive to the structural change.
2.4. Regression Analysis
Regression analysis gives the information on the relationship between a response variable and one or more independent variables. From regression analysis, the response variable is expressed as a function of the predictor variables. The relationship obtained from regression analysis can be used to predict values of the response variable and identify variables that most affect the response. The value of each predictor variable can be accessed through statistical tests on the estimated coefficients of the predictor variables.
An example of a regression model is the simple linear regression model which is a linear relationship between response variable (RV) and the predictor variable (T) of the form (8)RV=a+bTin which a and b are regression coefficients which are unknown modal parameters. In this study, the response variables are MSD of AR coefficient (2), CC of PSDs (3), and natural frequency and mode shape calculated by (6). Also, the predictor variable T is temperature measured from the experiment. The error term has to be equal to zero on average. In statistics, simple linear regression is the least squares estimator of a linear regression model with a single predictor variable. Simple linear regression fits a straight line through the set of n points in such a way that makes the sum of squared residuals of the model (i.e., vertical distances between the points of the data set and the fitted line) as small as possible.
3. Experiments on PSC Girder
As illustrated in Figure 1, dynamic tests were performed on a lab-scaled PSC girder. While room temperatures were handled to vary in the range of about 5~23°C, a set of prestress cases were simulated to the PSC girder from which vibration responses were measured to determine vibration features and modal parameters. In this study, vibration features of interest include autoregressive coefficient and power spectral density. Also, modal parameters of interest are natural frequencies and mode shape. Note that damping properties of the PSC girder are not dealt with due to the difficulties in measurement and interpretation.
Vibration tests on PSC girder.
Test setup
Test girder
Prestress force control
The PSC girder was simply supported and installed on a rigid testing frame. Two simple supports were modeled by steel rods between the girder and the rigid frame. As detailed in existing publications [21], the PSC girder model has the T-section reinforced in both longitudinal and transverse direction with 10 mm diameter reinforcing bars. As the prestressing tendon, a seven-wire monostrand with 15.2 mm diameter was embedded in a 25 mm diameter duct. During the tests, the prestress force was fixed as 98.0 kN.
Seven accelerometers noted as Sensors 1–7 in Figure 1(a) were placed in the girder with a constant 1 m interval. The impact excitation was applied in vertical direction by an electromagnetic shaker VTS100 at a location 1.7 m distanced from the right edge. The impact magnitude of the impact testing was consistently maintained during the experimental tests. Each accelerometer (PCB 393B04) has nominal sensitivity of 1 V/g and specified frequency range (±5%) of 0.06–450 Hz. The sampling frequency of 500 Hz was used to measure dynamic responses. The data acquisition system consists of a 16-channel PXI-4472 DAQ, a NI controller with LabVIEW, and MATLAB software. Temperature data were measured by using K-type thermocouple wires and KYOWA (EDX-100A) Dynamic Logger.
A series of tests were performed for 9 consecutive days. Figure 2 shows time history of room temperature during the test period. During the tests on the PSC girder, temperatures varied between 5°C and 23°C. Humidity in the laboratory was kept close to 40~45% in order to minimize the effect of humidity variation on the vibration characteristics. Temperature was controlled as designed for the tests by air conditioners and heaters (e.g., heater on and off). Vibration tests started at 18:00 hour of January 8 as the laboratory temperature reached up to 23°C. Then the room temperature was controlled to be decreased gradually for the remaining days. It is noted that the room temperature changed day and night.
Time history of temperature during experiments.
4. Temperature Effect on Vibration Features of PSC Girder4.1. AR Model
As the first vibration feature, AR model is used to estimate temperature effects on vibration characteristics of the PSC girder. The overall steps are as follows: firstly, a set of vibration responses are measured for various temperatures from the reference PSC girder; secondly, vibration features are extracted for the temperature variation by computing AR coefficient (1) and MSD (2); thirdly, time histories of temperature and MSD values are analyzed; finally, regression analysis is performed by fitting (8) to estimate the linear relationship between temperature variation and the extracted AR features.
As the first step, acceleration responses were measured from the PSC girder. Figure 3(a) shows the acceleration signals measured from Sensor 5. For temperature 22.5°C, 2048 acceleration data were sampled for 4 seconds by a single impulse excitation. As the second step, AR features were extracted from the measured acceleration signals. Figure 3(b) shows partial autocorrelation data for the acceleration signals of Sensor 5, from which the order of AR model was decided as 50. Using the decided order 50, AR coefficients for all acceleration signals measured from Sensor 5 were estimated for temperatures 5.4°C~22.5°C, as shown in Figure 3(c).
Acceleration signals and AR coefficients for estimation of AR model.
Acceleration signal of Sensor 5: temperature 22.5°C
Partial autocorrelation of Sensor 5: temperature 22.5°C
AR coefficient of Sensor 5: temperatures 5.4°C to 22.5°C
As the final step, the linear regression between temperature variation and AR coefficients was analyzed by the use of MSDs of AR coefficients. Figure 4(a) shows the time history for MSDs of AR coefficients calculated for temperatures 5.4°C~22.5°C. Figure 4(b) shows the linear relationship between the temperatures and the MSDs of AR coefficients. The empirical equation of the MSD as a function of temperature (T) is as follows: (9)MSD=-3283T+74569in which the MSD value changes as 3283 when temperature (T) changes by 1°C. It is noted that the correlation level was low. It is also noted that the AR model provided less consistent but rather poor estimation of the temperature effect on the vibration characteristics of the PSC girder. It is also noted that the quality of the AR model was dependent upon the amount of sampled data and the type of excitation sources.
Linear regression analysis result for acceleration-based AR model.
Time history of MSD of AR coefficient
Temperature versus MSD of AR model
4.2. CC of PSD
As the second feature, CC of PSDs was selected to estimate temperature effect on vibration characteristics of the PSC girder. The overall steps are as follows: firstly, acceleration data and temperature data are measured from selected distributed locations; secondly, CC of PSDs (3) and LCL (5) are computed for feature extraction; thirdly, time histories of temperature and CC of PSDs values are analyzed; and finally, a linear regression is estimated by fitting (8) for the relationship between temperature variation and the CC of PSDs.
For temperatures 5.4°C~22.5°C, as shown in Figure 5(a), PSDs were extracted from the acceleration signals of Sensor 5, at which 2048 acceleration data were sampled for 4 seconds by a single impulse excitation. As shown in Figure 5(b), CC of PSDs were computed for temperatures 5.4°C~22.5°C by using PSDs at the maximum 22.5°C as the reference. It is observed that CC of PSDs varies along with the variation of room temperature. Figure 5(c) shows the linear regression between the temperatures and the CC of PSDs. The empirical equation of the CC of PSDs as a function of temperature (T) is as follows: (10)CC=0.019T+0.676in which the CC of PSDs changes as 0.019 when temperature (T) changes by 1°C. It is noted that the correlation level was relatively moderate according to Dancey and Reidy’s categorization [22]. However, the CC of PSDs gave relatively consistent estimation of the temperature effect on the vibration characteristics of the PSC girder.
Linear regression analysis for correlation coefficient of power spectral density.
PSDs at 5.4°C~22.5°C
Time history of CC of PSD
Temperature versus CC of PSD
4.3. Natural Frequency
Next, temperature effect on vibration feature is estimated by natural frequency of the PSC girder. The overall steps are as follows: first, acceleration data and temperature data are measured from selected distributed locations; secondly, from the FDD method defined as (6), natural frequencies were estimated by the singular values of the PSD function matrix S(ω); thirdly, time histories of temperature and natural frequencies are analyzed; and finally, a linear regression is estimated for the relationship between temperature variation and the extracted natural frequencies.
Acceleration signals of Sensors 1–7 were utilized for the modal extraction. From each sensor, 2048 acceleration data were sampled for 4 seconds by a single impulse excitation. Figure 6 shows time histories of temperatures and natural frequencies of the first two modes of the PSC girder. Natural frequencies of the two modes were extracted as temperatures varied between 5.4°C and 20.4°C. In mode 1, changes in natural frequencies were very small as temperatures varied up to 15°C, as shown in Figure 6(a). In mode 2, however, changes in natural frequencies were relatively high due to the temperature variation, as shown in Figure 6(b).
Time history of temperature for natural frequencies.
Mode 1
Mode 2
From regression analysis, linear relationships between the two modes’ natural frequencies and temperatures were estimated as shown in Figure 7. For the two modes, the empirical equations of the natural frequencies as a function of temperature (T) are, respectively, as follows: (11a)Freq1=-0.027T+20.93,(11b)Freq2=-0.113T+105.08in which the first mode’s natural frequencies (Freq1) decrease as 0.027 as temperature (T) increases by 1°C; and the second mode’s natural frequencies (Freq2) decrease as 0.113 as temperature (T) increases by 1°C. It is noted that the correlation levels were relatively moderate for mode 1 but relatively strong for mode 2 [22]. It is also noted that natural frequencies gave relatively consistent estimation of the temperature effect on the vibration characteristics of the PSC girder.
Linear relationship between temperature and natural frequencies.
Mode 1
Mode 2
4.4. Mode Shape
Next, temperature effect on vibration feature is estimated by mode shapes of the PSC girder. The overall steps are as follows: first, acceleration data and temperature data are measured from selected distributed locations; secondly, mode shapes were estimated by the singular values of the PSD function matrix S(ω) from the FDD method (6); thirdly, MAC values are calculated by (7); fourthly, time histories of temperature and MAC values are analyzed; and finally, a linear regression is estimated for the relationship between temperature variation and the MAC values.
The same sets of acceleration signals and temperature signals as used for the analysis of natural frequencies were utilized for the analysis of mode shapes. Figure 8 shows mode shapes of the PSC girder measured by Sensors 1–7 as temperatures varies between 5.4°C and 22.5°C. It is observed that overall changes in mode shapes due to temperature variation look small in both mode 1 and mode 2. Figure 9 shows time histories of temperatures and mode shapes of the first two modes. In mode 1, the change in mode shape was very small as temperatures varied up to 15°. In mode 2, however, the change in mode shape was relatively high due to the temperature variation.
Natural frequency and mode shape on temperature conditions for PSC girder: prestress force 98 kN, temperatures 5.4°C–22.5°C.
Mode 1: 20.20 Hz–21.04 Hz
Mode 2: 103.33 Hz–104.65 Hz
Time history of temperature for mode shapes of PSC girder.
Mode 1
Mode 2
From regression analysis, linear relationships between the two mode shapes and temperatures were estimated as shown in Figure 10. For the two modes, the empirical equations of the mode shapes as a function of temperature (T) are, respectively, as follows: (12a)MAC1=-0.00002T+0.9999,(12b)MAC2=-0.0002T+0.9976in which the first mode’s modal assurance criterion (MAC1) almost does not change as temperature (T) changes; but the second mode’s modal assurance criterion (MAC2) changes as 0.0002 as temperature (T) increases by 1°C. It is noted that the correlation levels were relatively weak in mode 1 but relatively moderate in mode 2 [22]. It is also noted that modal assurance criteria of the two modes (MAC1 and MAC2) gave relatively consistent estimation of the temperature effect on the vibration characteristics of the PSC girder.
Linear relationship between temperature and mode shapes.
Mode 1
Mode 2
5. Summary and Conclusion
In this paper, the effect of temperature variation on vibration monitoring of prestressed concrete (PSC) girders was experimentally analyzed. Firstly, vibration features such as autoregressive (AR) coefficient, correlation coefficient of power spectral density (CC of PSD), natural frequency, and mode shape were selected to estimate the effect of temperature variation on PSC girders. Secondly, vibration experiments on a lab-scale PSC girder were performed under the condition of temperature variation. Finally, the selected vibration features were analyzed to estimate temperature-induced effects on vibration monitoring of the PSC girder.
Major results from the analysis on the four vibration features with respect to temperature variation are summarized as follows. Firstly, the AR model produced rather poor estimation of the temperature effect on the vibration characteristics of the PSC girder. Secondly, the CC of PSD gave relatively consistent estimation of the temperature effect on the vibration characteristics of the PSC girder. However, its correlation level of the linear regression analysis was relatively moderate. Thirdly, the natural frequencies produced relatively consistent estimation of the temperature effect on the vibration characteristics of the PSC girder. For the natural frequencies, the correlation levels of the linear regression analysis were relatively moderate in mode 1 but relatively strong in mode 2. Finally, modal assurance criteria of the first two modes gave relatively consistent estimation of the temperature effect on the vibration characteristics of the PSC girder. For the mode shapes, the correlation levels of the linear regression analysis were relatively weak in mode 1 but relatively moderate in mode 2. From the comparative analyses, it can be said that the natural frequency is highly correlated with temperature.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was supported by a grant from a strategic research project (Development of Smart Prestressing and Monitoring Technologies for Prestressed Concrete Bridges) funded by the Korea Institute of Construction Technology. Graduate students involved in this research were also partially supported by the Brain Korea 21 Plus program of Korean Government.
NawyE. G.1996New York, NY, USAPrentice HallAdamsR. D.CawleyP.PyeC. J.StoneB. J.A vibration technique for non-destructively assessing the integrity of structures19782029310010.1243/JMES_JOUR_1978_020_016_022-s2.0-0017951397StubbsN.OseguedaR.Global non-destructive damage evaluation in solids1990526779KimJ. T.StubbsN.Model-uncertainty impact and damage-detection accuracy in plate girder199512110140914172-s2.0-0029394559DoeblingS. W.FarrarC. R.PrimeM. B.A summary review of vibration-based damage identification methods19983029110510.1177/0583102498030002012-s2.0-0032023438KimJ.-T.RyuY.-S.ChoH.-M.StubbsN.Damage identification in beam-type structures: frequency-based method vs mode-shape-based method2003251576710.1016/S0141-0296(02)00118-92-s2.0-0037210189SaiidiM.DouglasB.FengS.Prestress force effect on vibration frequency of concrete bridges199412072233224110.1061/(ASCE)0733-9445(1994)120:7(2233)2-s2.0-0028467061MiyamotoA.TeiK.NakamuraH.BullJ. W.Behavior of prestressed beam strengthened with external tendons200012691033104410.1061/(ASCE)0733-9445(2000)126:9(1033)2-s2.0-0034276001KimJ.-T.YunC.-B.RyuY.-S.ChoH.-M.Identification of prestress-loss in PSC beams using modal information2004173-44674822-s2.0-154234891910.12989/sem.2004.17.3_4.467LawS. S.LuZ. R.Time domain responses of a prestressed beam and prestress identification20052884-51011102510.1016/j.jsv.2005.01.0452-s2.0-26444532832HuthO.FeltrinG.MaeckJ.KilicN.MotavalliM.Damage identification using modal data: experiences on a prestressed concrete bridge2005131121898191010.1061/(ASCE)0733-9445(2005)131:12(1898)2-s2.0-29044435147WahabM. A.de RoeckG.Effect of temperature on dynamic system parameters of a highway bridge1997742662702-s2.0-1542397331SohnH.DzwonczykM.StraserE. G.KiremidjianA. S.LawK.MengT.An experimental study of temperature effect on modal parameters of the Alamosa Canyon Bridge1999287-887989710.1002/(SICI)1096-9845(199908)28:8<879::AID-EQE845>3.0.CO;2-V2-s2.0-0032833429FarrarC. R.CornwellP. J.DoeblingS. W.Structural health monitoring studies of the Alamosa Canyon and I-40 bridges2000LA-13635-MSKimJ. T.YunC. B.YiJ. H.Temperature effects on frequency-based damage detection in plate-girder bridges200376725733NiY. Q.HuaX. G.FanK. Q.KoJ. M.Correlating modal properties with temperature using long-term monitoring data and support vector machine technique20052712176217732-s2.0-2384450452210.1016/j.engstruct.2005.02.020HongD.-S.NguyenK.-D.LeeI.-C.KimJ.-T.Temperature-compensated damage monitoring by using wireless acceleration-impedance sensor nodes in steel girder connection201220121216712010.1155/2012/1671202-s2.0-84867960039BendatJ. S.PiersolA. G.2003New York, NY, USAWiley-InterscienceOtteD.de PonseeleP. V.LeuridanJ.Operational shapes estimation as a function of dynamic loadsProceedings of the 8th International Modal Analysis ConferenceJanuary 1990Florida, Fla, USAYiJ.-H.YunC.-B.Comparative study on modal identification methods using output-only information2004173-444546610.12989/sem.2004.17.3_4.4452-s2.0-1542274235HoD.-D.KimJ.-T.StubbsN.ParkW.-S.Prestress-force estimation in PSC girder using modal parameters and system identification2012156997101210.1260/1369-4332.15.6.9972-s2.0-84864136373DanceyC. P.ReidyJ.2004Prentice Hall