Cable Modal Parameter Identification. I: Theory

Cable modal parameters (natural frequencies and damping ratios) that represent the cable inherent dynamic characteristics play an important role in the construction, vibration control, condition assessment, and long-term health monitoring of cable-supported structures. The existing options to identify cable modal parameters through vibration measurements are somewhat limited. For this purpose, a cable dynamic stiffness based method is presented to effectively identify the cable modal parameters. In the first part of this two-part paper, the cable dynamic stiffness is analytically discussed for a viscously damped, uniform, inclined sagging cable supported at the lower end and subjected to a harmonically varying arbitrary angle displacement excitation in an arbitrary angle at the upper end when the cable is assumed to have a parabolic profile at its position of static equilibrium. Special attention is paid to the physical meaning and significance of every part of the frequency-dependent closed-form cable dynamic stiffness. Comprehensive numerical analyses have been carried out and a simplified cable dynamic stiffness is proposed for the purpose of identifying the cable modal parameters with a good accuracy over a wide range of frequencies.     

Bayesian Modal Updating using Complete Input and Incomplete Response Noisy Measurements

The problem of identification of the modal parameters of a structural model using complete input and incomplete response time histories is addressed. It is assumed that there exist both input error (due to input measurement noise) and output error (due to output measurement noise and modeling error). These errors are modeled by independent white noise processes, and contribute towards uncertainty in the identification of the modal parameters of the model. To explicitly treat these uncertainties, a Bayesian framework is adopted and a Bayesian time-domain methodology for modal updating based on an approximate conditional probability expansion is presented. The methodology allows one to obtain not only the optimal (most probable) values of the updated modal parameters but also their uncertainties, calculated from their joint probability distribution. Calculation of the uncertainties of the identified modal parameters is very important if one plans to proceed with the updating of a theoretical finite-element model based on these modal estimates. The proposed approach requires only one set of excitation and corresponding response data. It is found that the updated probability density function (PDF) can be well approximated by a Gaussian distribution centered at the optimal parameters at which the posterior PDF is maximized. Numerical examples using noisy simulated data are presented to illustrate the proposed method.     

Combined Experimental-Operational Modal Testing of Footbridges

In combined vibration testing, an artificial, measured force is used in operational conditions. This requires the identification of a system model that takes both the measured and the operational excitation into account. Advantages with respect to the classical operational modal analysis approach are the possibility of obtaining mass-normalized mode shapes and the increase of the excitation level and its frequency content. An advantage with respect to the classical experimental modal analysis approach, where the ambient excitation is not modeled, but considered as disturbing noise, is the possibility of using excitation levels that are of the same amplitude, or even smaller, than the ambient excitation levels. In this paper, combined modal testing of footbridges is explored using two case studies: a steel arch footbridge with spans of 75.2 m and 30.3 m and a concrete stress-ribbon footbridge with spans of 30 m and 28 m. The comparison of the modal parameters (eigenfrequencies, damping ratios, mode shapes, and modal scaling factors) obtained from a combined vibration test with the ones obtained from other modal tests and from a finite-element model, demonstrates the feasibility of using small and practical excitation devices for the modal testing of footbridges.     

Stochastic Response of an Inclined Shallow Cable with Linear Viscous Dampers under Stochastic Excitation

Considering the coupling between the in-plane and out-of-plane vibration, the stochastic response of an inclined shallow cable with linear viscous dampers subjected to Gaussian white noise excitation is investigated in this paper. Selecting the static deflection shape due to a concentrated force at the dampers location and the first sine term as shape functions, a reduced four-degree-of-freedom system of nonlinear stochastic ordinary differential equations are derived to describe dynamic response of the cable. Since only polynomial-type terms are contained, the fourth-order cumulant-neglect closure together with the C-type Gram-Charlier expansion with a fourth-order closure are applied to obtain statistical moments, power spectral density and probabilistic density function of the cable response, whose availability is verified by Monte Carlo method. Taking a typical cable as an example, the influence of several factors, which include excitation level and direction as well as damper size, on the dynamic response of the cable is extensively investigated. It is found that the sum of mean square in-plane and out-of-plane displacement is primarily independent of the load direction when the excitation level and viscous coefficient of the damper are fixed. Moreover, the peak frequency and half-band width of the spectra of both the in-plane and the out-of-plane displacements are increasing with excitation level when the damper size is constant. It is also observed that, even though the actual optimal damper size is slightly greater than the one obtained by the complex modal theory, the difference of statistical moment of the cable caused by these two damper size is negligible, so the vibration reduction effect provided by the theoretical optimal viscous coefficient is satisfactory.     

Vibration Control of Stay Cables of the Shandong Binzhou Yellow River Highway Bridge Using Magnetorheological Fluid Dampers

The Shandong Binzhou Yellow River Highway Bridge is a three-tower, cable-stayed bridge in Shandong Province, China. Because the stay cables are prone to vibration, 40 magnetorheological (MR) fluid dampers were attached to the 20 longest cables of this bridge to suppress possible vibration. An innovative control algorithm for active and semiactive control of mass-distributed dynamic systems, e.g., stay cables, was proposed. The frequencies and modal damping ratios of the unimpeded tested cable were identified through an ambient vibration test and free vibration tests, respectively. Subsequently, a series of field tests were carried out to investigate the control efficacy of the free cable vibrations achieved by semiactive MR dampers, “Passive-off” MR dampers and “Passive-on” MR dampers. The first three modal damping ratios of the cable incorporated with the MR dampers were also identified from the in situ experiments. The field experiment results indicated that the semiactive MR dampers can provide significantly greater supplemental damping for the cable than either the Passive-off or the Passive-on MR dampers because of the pseudonegative stiffness generated by the semiactive MR dampers.     

Dynamic Field Test, System Identification, and Modal Validation of an RC Minaret: Preprocessing and Postprocessing the Wind-Induced Ambient Vibration Data

The investigation of dynamic response for civil engineering structures largely depends on a detailed understanding of their dynamic characteristics, such as the natural frequencies, mode shapes, and modal damping ratios. Dynamic characteristics of structures may be obtained numerically and experimentally. The finite-element method is widely used to model structural systems numerically. However, there are some uncertainties in numerical models. Material properties and boundary conditions may not be modeled correctly. There may be some microcracks in the structures, and these cracks may directly affect the modeling parameters. Modal testing gives correct uncertain modeling parameters that lead to better predictions of the dynamic behavior of a target structure. Therefore, dynamic behavior of special structures, such as minarets, should be determined with ambient vibration tests. The vibration test results may be used to update numerical models and to detect microcracks distributed along the structure. The operational modal analysis procedure consists of several phases. First, vibration tests are carried out, spectral functions are produced from raw measured acceleration records, dynamic characteristics are determined by analyzing processed spectral functions, and finally analytical models are calibrated or updated depending on experimental analysis results. In this study, an ambient vibration test is conducted on the minaret under natural excitations, such as wind effects and human movement. The dynamic response of the minaret is measured through an array of four trixial force-balanced accelerometers deployed along the whole length of the minaret. The raw measured data obtained from ambient vibration testing are analyzed with the SignalCAD program, which was developed in MATLAB. The employed system identification procedures are based on output-only measurements because the forcing functions are not available during ambient vibration tests. The ModalCAD program developed in MATLAB is used for dynamic characteristic identification. A three-dimensional model of the minaret is constructed, and its modal analysis is performed to obtain analytical frequencies and mode shapes by using the ANSYS finite-element program. The obtained system identification results have very good agreement, thus providing a reliable set of identified modal properties (natural frequencies, damping ratios, and mode shapes) of the structure, which can be used to calibrate finite-element models and as a baseline in health monitoring studies.     

Stiffness Evaluation and Damage Detection of Ceramic Candle Filters

This paper presents a nondestructive evaluation method to identify the structural stiffness of ceramic candle filters. A ceramic candle filter is a hollow cylindrical structure made of a porous ceramic material used in advanced, coal-fired power generation systems. The candle filters need to sustain an extreme thermal and chemical environment over a great period of time to protect the gas turbine components from exposure to particulate matter. A total of 92 new candle filters and 29 used candle filters have been tested nondestructively using a dynamic characterization technique. All filters were subjected to an excitation force, and the response was picked up by an accelerometer in a free-free boundary condition. The frequency response function and vibration mode shapes of each filter were evaluated. Beam vibration equations and finite-element models were built to calculate the filter's dynamic response. Results indicate that the vibration signatures can be used as an index to quantify the structural properties of ceramic candle filters. The results also show estimations of the overall bending stiffness values for four different types of candle filters. The used filters show a trend of stiffness degradation, which was related to the filter's exposure time. Damage detection procedures using modal strain energy and finite-element simulation were studied for detection of a localized damage in the candle filter. The location and the size of the damaged section can be identified using the measured model strain energy.     

Fast Bayesian FFT Method for Ambient Modal Identification with Separated Modes

Previously a Bayesian theory for modal identification using the fast Fourier transform (FFT) of ambient data was formulated. That method provides a rigorous way for obtaining modal properties as well as their uncertainties by operating in the frequency domain. This allows a natural partition of information according to frequencies so that well-separated modes can be identified independently. Determining the posterior most probable modal parameters and their covariance matrix, however, requires solving a numerical optimization problem. The dimension of this problem grows with the number of measured channels; and its objective function involves the inverse of an ill-conditioned matrix, which makes the approach impractical for realistic applications. This paper analyzes the mathematical structure of the problem and develops efficient methods for computations, focusing on well-separated modes. A method is developed that allows fast computation of the posterior most probable values and covariance matrix. The analysis reveals a scientific definition of signal-to-noise ratio that governs the behavior of the solution in a characteristic manner. Asymptotic behavior of the modal identification problem is investigated for high signal-to-noise ratios. The proposed method is applied to modal identification of two field buildings. Using the proposed algorithm, Bayesian modal identification can now be performed in a few seconds even for a moderate to large number of measurement channels.     

General Realization Algorithm for Modal Identification of Linear Dynamic Systems

The general realization algorithm (GRA) is developed to identify modal parameters of linear multi-degree-of-freedom dynamic systems subjected to measured (known) arbitrary dynamic loading from known initial conditions. The GRA extends the well known eigensystem realization algorithm (ERA) based on Hankel matrix decomposition by allowing an arbitrary input signal in the realization algorithm. This generalization is obtained by performing a weighted Hankel matrix decomposition, where the weighting is determined by the loading. The state-space matrices are identified in a two-step procedure that includes a state reconstruction followed by a least-squares optimization to get the minimum prediction error for the response. The statistical properties (i.e., bias, variance, and robustness to added output noise introduced to model measurement noise and modeling errors) of the modal parameter estimators provided by the GRA are investigated through numerical simulation based on a benchmark problem with nonclassical damping.     

Modal Testing, Finite-Element Model Updating, and Dynamic Analysis of an Arch Type Steel Footbridge

This paper describes an arch type steel footbridge, its analytical modeling, modal testing, finite-element model updating, and dynamic analysis. A modern steel footbridge which has an arch type structural system and is located on the Karadeniz coast road in Trabzon, Turkey is selected as an application. An analytical modal analysis is performed on the developed three-dimensional finite-element model of footbridge to provide analytical frequencies and mode shapes. Field ambient vibration tests on the footbridge deck under natural excitation such as human walking and traffic loads are conducted. The output-only modal parameter identification is carried out by using peak picking of the average normalized power spectral densities in the frequency domain and stochastic subspace identification in the time domain, and dynamic characteristics such as natural frequencies, mode shapes, and damping ratios are determined. The finite-element model of the footbridge is updated to minimize the differences between analytically and experimentally estimated modal properties by changing some uncertain modeling parameters such as material properties. Dynamic analyses of the footbridge before and after finite-element model updating are performed using the 1992 Erzincan earthquake record. At the end of the study, maximum differences in the natural frequencies are reduced from 22 to only 5% and good agreement is found between analytical and experimental dynamic characteristics such as natural frequencies and mode shapes by model updating. Also, maximum displacements and principal stresses before and after model updating are compared with each other.     

Cross Correlations of Modal Responses of Tall Buildings in Wind-Induced Lateral-Torsional Motion

Recent trends towards developing increasingly taller and irregularly-shaped buildings imply that these complex structures are potentially more responsive to wind excitation. Making accurate predictions of wind loads and their effects on such structures is therefore a necessary step in the design synthesis process. This paper presents a framework for dynamic analysis of the wind-induced lateral-torsional response of tall buildings with three-dimensional (3D) mode shapes. The cross correlation reflecting the statistical coupling among modal responses under spatiotemporally varying dynamic wind excitations has been investigated in detail. The effects of intermodal correlations on the lateral-torsional response of tall buildings with 3D mode shapes and closely spaced natural frequencies are elucidated and a more accurate method for quantifying intermodal cross correlations is analytically developed. Utilizing the wind tunnel derived synchronous multipressure measurements, a full-scale 60-story asymmetric building of mixed steel and concrete construction is used to illustrate the proposed framework for the coupled dynamic analysis and highlight the intermodal correlation of modal responses on the accurate prediction of coupled building acceleration.     

Damping of Cables by a Transverse Force

A general asymptotic format is presented for the effect on the modal vibrations of a transverse damper close to the end of a cable. Complete locking of the damper leads to an increase of the natural frequencies, and it is demonstrated that the maximum attainable damping is a certain fraction of the relative frequency increase, depending on the type of damping device. The asymptotic format only includes a real and a complex nondimensional parameter, and it is demonstrated how these parameters can be determined from the frequency increase by locking and from an energy balance on the undamped natural vibration modes. It is shown how the asymptotic format can incorporate sag of the cable, and specific results are presented for viscous damping, the effect of stiffness and mass, fractional viscous damping, and a nonlinear viscous damper. The relation of the stiffness component to active and semiactive damping is discussed.     

Dynamic Behavior of Stay Cables with Rotational Dampers

Vibration reduction in stay cables by means of viscous dampers is of great interest in cable damage prevention and serviceability of structural system supported by such cables. The paper presents a study on the effectiveness, as well as the limits, of rotational viscous dampers and springs inserted at the two ends of a bending-stiff taut cable; influence of rotational stiffness of the springs is also investigated. After a nondimensional expression of the equation of motion has been obtained, as in other cases of nonproportionally damped continuous structures, complex modal analysis is pursued, obtaining complex eigenvalues and eigenfunctions. Comparison with intermediate dampers, widely used in bridge engineering, is performed showing the range of nondimensional parameters for which the proposed approach is of interest. Finally, a numerical technique based on complex mode superposition is presented in order to evaluate time domain responses for transversal distributed excitation. As an example, the procedure is applied to a wind-exposed cable.     

Two-Stage Structural Health Monitoring Approach for Phase I Benchmark Studies

This paper presents a two-stage structural health monitoring methodology and applies it to the Phase I benchmark study sponsored by the IASC-ASCE Task Group on Structural Health Monitoring. In the first stage, modal parameters are identified using measured structural response from the undamaged system and then from the (possibly) damaged system. In the second stage, these data are used to update a parametrized structural model of the system using Bayesian system identification. The approach allows one to obtain not only estimates of the stiffness parameters but also the probability that damage in any substructure exceeds any specified threshold expressed in terms of a fractional stiffness loss. It successfully identifies the location and severity of damage in all cases of the benchmark problem.     

Modal Contribution Coefficients in Bridge Condition Evaluation

Impact modal testing combined with finite element (FE) analyses is currently being used to evaluate the condition of steel bridges in the state of Ohio. Using modal testing techniques, it is relatively easy to measure the dynamic response of bridges, including mode shapes, frequencies, and modal scaling factors. These responses are compared to the results of the FE analyses and the model is iteratively updated until a good agreement is obtained. After a good agreement between experimental and analytical results has been achieved, the FE model is used to obtain stresses that are used to load rate the bridge. During the iterative calibration process, several quantities, including the fundamental mode shapes and frequencies, are used to evaluate the accuracy of the FE model. Since each mode shape plays a different role in the dynamic behavior of the structure, a more efficient calibration routine can be achieved if more emphasis is placed on obtaining good matches for the modes that are most influential. The aim of this paper is to develop a quantitative measure of the contribution of different modes to the overall dynamic response of a structure. The proposed measure, a series of contribution coefficients, is used to identify which modes are most critical in the process of modal testing and FE model calibration. Several applications of the contribution coefficients are identified.     

Dynamic Testing of Prestressed Concrete Bridges and Numerical Verification

This paper describes different dynamic tests carried out on three bridges, B14, B15, and B13, over highway E19, which connects Brussels and Antwerpen in Belgium. Different excitation types are considered: a drop weight, a heavy truck on the bridge, and ambient vibrations mainly due to the traffic under the bridge. Finite-element models are constructed to support and verify the dynamic measurements. The modal parameters are extracted from the response time series using the data dependent system approach. A vector autoregressive model is developed and successfully applied to the measured responses in time domain. Good correlation between the finite-element simulations and the experiments is obtained.     

Design of Mechanical Viscous Dampers for Stay Cables

External dampers have been utilized in a number of cable-stayed bridges to suppress transverse cable vibrations. However, simple and accurate damper design recommendations that concurrently consider all important cable parameters are lacking. Previous efforts have been based on the idealization of cables as taut strings. In this paper, the governing differential equation for vibration of cables containing a viscous damper was first converted to a complex eigenvalue problem containing nondimensional cable parameters. Then, a parametric study was conducted involving repeated solutions of the eigenvalue problem for a wide range of nondimensional parameters. Based on the results of the parametric study, the effects of dampers on first mode vibration frequencies and first mode cable damping ratios were presented in nondimensional format. It is shown that for the range of parameters involved in most stay cables, the influence of cable sag is insignificant, whereas the cable bending stiffness can have a significant influence on the resulting cable damping ratios. Simplified nondimensional relationships are proposed for calculating damper-induced changes in the first mode cable damping ratios. Results of laboratory tests on a scaled model cable are compared with the estimated values using the formulation presented. Finally, example problems are presented for comparison with other relationships, and for the design of mechanical viscous dampers for suppression of cable vibrations including rain-wind induced vibrations.     

Parameter Estimation for Multiple-Input Multiple-Output Modal Analysis of Large Structures

Experimental modal analysis (EMA) has been explored as a technology for condition assessment and damage identification of constructed structures. However, successful EMA applications such as damage detection to constructed systems pose certain difficulties. The properties of constructed systems are influenced by temperature changes as well as other natural influences such as movements in addition to any deterioration and damage. Writers were challenged in their attempts to measure the dynamic properties of an aged bridge by EMA due to inconsistencies within the data set due to short-term variations in ambient conditions. A complex interaction was observed between the dynamic properties of the bridge, hour-to-hour changes in temperature, and controlled damages applied to the bridge. Inconsistencies in the data set made curve fitting difficult for some common parameter estimation algorithms that have been designed to handle consistent data sets. Although the quality of measurements within the entire data set was affected by time variance and nonlinearity, increasing the number of reference measurements significantly improved the reliability of the information which could be extracted. In conjunction with the multiple-input multiple-output technique, a parameter estimation method using complex mode indicator function (CMIF) was developed and implemented in this study to determine the modal properties with proper scaling to obtain modal flexibility. This method proved to be very successful among many others with the data acquired from the aged and deteriorated highway bridge. In this paper, challenges in reliable identification of modal parameters from large structures are reviewed and the new CMIF based algorithm is documented. The method is evaluated on actual bridge data sets from a damage detection research study.     

Space Station Structure Modal Selection Criteria

In many dynamic analysis procedures the size of the problem is curtailed by truncation of the modes selected for study. A preferable approach in the case of a structure such as a large space frame is to identify and include those natural modes in the order of their increasing importance. Two of the many criteria which can be used to determine importance are examined in this paper. The first involves establishment of a “completeness index” reflecting the model identities satisfied by the natural frequencies and integrals of the mode shapes. The second involves displacements at a point of excitation and a response point. Both selection methods are implemented on a simplified model of the Initial Operating Configuration Space Station and comparative transient response analyses are made. It is shown that mode selection is a function of more than one variable. The first method performs well in reducing model error yet the second is necessary to predict those modes excited by various forcing functions when control of displacements is the chief concern. Different mode selection criteria will be necessary when other variables are involved.     

Finite-Element Analysis of Time Domain Reflectometry Cable-Grout-Soil Interaction

Three-dimensional finite element analysis is combined with field and laboratory measurement of time domain reflectometry (TDR) cable-grout response to analyze the interaction between the cable, grout, and surrounding soil mass during localized shearing. Finite element (FE) model parameters for the cable and cable-grout interface elements are back-calculated by matching results from laboratory shearing tests to FE calculated response. These parameters are employed in subsequent FE model geometries to model the behavior of TDR cable-grout composites in soft soils. Optimal grout and cable design is determined by analyzing the relationship between grout strength and stiffness and calculated cable shear stress.