Vol50,2,2007 259 ANNALS OF GEOPHYSICS, VOL. 50, N. 2, April 2007 Key words damping – fundamental mode – single- degree-of-freedom – ambient vibration 1. Introduction The need for a simple and reliable method to estimate building damping arose during the proj- ect NATO Science for Peace 980857 (ASSAS- BV) whose aim is to mitigate seismic risk in the cities of Skopje, Zagreb and Ljubljana. One of the tasks of the project is to perform estimates of the fundamental frequency of buildings and of underlying soils, in order to identify possible res- onance phenomena. This resonance approach was recently introduced for microzonation stud- ies of large cities (e.g., Granada, Spain; Navarro et al., 2004). Obtaining fundamental frequency and damping is the first step to validate vulnera- bility models of buildings. Standard and detailed approaches can be used only for a limited number of selected buildings: the challenge faced in the ASSASBV project was to obtain simple but reli- able estimates of dynamic parameters for thou- sands of building in a two-year time. This has to be achieved with the same instrumentation used on soils and without induced forcing. The proposed methodology does not require filtering, pre-selection or assumptions on signal Non-parametric analysis of a single seismometric recording to obtain building dynamic parameters Marco Mucciarelli (1) and Maria Rosaria Gallipoli (1)(2) (1) Dipartimento di Strutture, Geotecnica, Geologia Applicata all’Ingegneria (DiSGG), Facoltà di Ingegneria, Università degli Studi della Basilicata, Potenza, Italy (2) Istituto di Metodologie per l’Analisi Ambientale (IMAA), CNR, Tito Scalo (PZ), Italy Abstract We propose a simple time-domain, non-parametric method to estimate the damping at the fundamental frequen- cy of a building. The method aims at obtaining quick-and-dirty data on large sets of buildings, at the expenses of the accuracy provided by other, more complex and resource-demanding techniques. The analysis of a 10 min recording of ambient vibration with a single high-resolution seismometer atop the building can provide a good estimate of the required parameters for the first flexural modes on orthogonal components. The proposed methodology does not require complex filtering and assumptions on signal structure, nor multiple measurement points or clear single transients induced by shakers, shocks or release tests. We checked the stability of the pro- posed method in terms of duration and characteristics of the signal, and compared the results obtained by others with standard techniques. Then, we tested the ability of the proposed technique to identify damping and frequen- cy variations due to large displacements, damage or changes in the structural characteristics. The proposed methodology provides a satisfactory agreement when compared with other techniques. Even if it is not always possible to obtain higher modes, the advantage is that it is possible to study with limited resources the funda- mental parameters for a large number of buildings. This is useful to include experimental data on building be- haviour in microzonation studies. Mailing address: Dr. Marco Mucciarelli, Dipartimento di Strutture, Geotecnica, Geologia Applicata all’Ingegneria (DiSGG), Facoltà di Ingegneria, Università degli Studi del- la Basilicata, Viale dell’Ateneo Lucano 10, 85100 Potenza Italy; e-mail: mucciarelli@unibas.it 260 Marco Mucciarelli and Maria Rosaria Gallipoli structure, nor multiple measurement points or clear single transients induced by shakers, shocks or release test. This overcomes some of the operational disadvantages of well established approaches like RandomDec (Vandiver et al., 1982), half-power method (Clough and Penzien, 1975), cross-correlation or wavelet-logarithmic decrement (Lamarque et al., 2000). For a de- tailed review of pros and cons of different ambi- ent and forced vibration tests see Huang et al. (1999, and references therein), that also present a method applicable to multiple measurement points. A generally accepted assumption is that the first approximation of the dynamic behaviour of a building is provided by the Single-Degree- Of-Freedom damped oscillator (SDOF) (1.1) where γ is the damping and ω0 is the fundamen- tal frequency, in its turn related to the mass m and stiffness k of the system as . (1.2) Even with this simple formulation, the system identification becomes a problem when the forc- ing term f(t) is unknown. With ambient vibra- tion excitation, a further warning comes from the difficulty to separate the noise filtered from the structure from the instrumental noise that is unavoidable during the measurements. The same problem is well known from microtremor measurements on soils, see e.g., Mucciarelli (1998) and Mucciarelli et al. (2005). The solution of eq. (1.1) is given by (1.3) where x0 is the initial amplitude, ϕ the phase and ε a perturbation term. The damping is commonly derived as the logarithmic decrement of the amplitude over a pseudo-period T according to the following equation: . (1.4) Now we face the problem of numerical estima- ( ) / ( )ln T x t t Tx 1 γ = + ( ) ( )cosx t x e tt0 ω ϕ ε= + +γ− m k 0ω = ( )x x x f t2 0 2γ ω+ + =p o tion of the above parameters from a discrete time series. If some form of excitation is ap- plied, the damping can be straightforwardly es- timated picking the amplitudes of subsequent maxima soon after the excitation is stopped and the system returns to its equilibrium position. For ambient vibration recordings, the Ran- domDec approach solves this problem using a stacking of pre-selected segments exceeding a fixed amplitude threshold xl with common ini- tial conditions 2. Proposed method and numerical simulations Assuming that the ground ambient vibration is close to white noise, a simple idealisation of the ambient vibration recording atop a building is composed by three factors (2.1) where xω0(ω) is the oscillation associated to the fundamental period (or x (t) in eq. (1.3)), that cannot be observed in its full development be- cause of the superposition of xεg(ω) that is the perturbation induced by the unknown, continu- ous forcing induced by ground ambient vibra- tion (or ε in eq. (1.3)). The term xεi(ω) repre- sents the effect of instrumental noise during the measurement, as well as other sources of noise unfiltered by the building (e.g., wind impacting on unprotected sensors on a terrace roof). If we suppose that xω0(ω) > xεg(ω) and xω0(ω) >> >>xεi (ω), then we can try to retrieve γ directly from the recordings following this procedure: i) Use a high resolution seismometer for taking measurements of ambient noise time his- tories. Accelerometers provide insufficient re- solving power, as shown for measurements on soil by Mucciarelli (1998) and confirmed recent- ly by the SESAME project (Duval et al., 2004). ii) After a standard correction procedure (baseline correction and high-pass filtering), in- tegrate the time-history to obtain displacements. iii) Select only positive values that are local maxima. Build a 2 × m matrix: amplitudes x(ti) of maximum in one column and time of occur- rence ti in the other. ( ) ( ) ( ) ( )X x x xg i0ω ω ω ω= + +ω ε ε ( ( ) ; ( ) ( ) ) .x t x x t t x t 0>∆= + −l 261 Non-parametric analysis of a single seismometric recording to obtain building dynamic parameters iv) For each subsequent row in the above matrix, verify if x(ti)>x(ti+1); if yes, calculate piecewise pseudo-frequency and damping as (2.2) . (2.3) v) Build a 2-histogram, binning piecewise pseudo-frequencies and damping. vi) Produce a graphic output to visualise the maximum of the empirical distribution, that is the main result of this non-parametric damp- ing analysis (hereafter NonPaDAn), providing the fundamental frequency and its associated damping. We implemented the above procedure using MatLab: the relevant m-file is freely available for download at the address . To verify the above procedure and hypothe- sis, we built a numerical model using SimuLink. Figure 1 shows the relevant block diagram. The input is the sum of repeated, large pulses and ( / ) ( ( ) / ( ))x t x t1 2i i i 1:γ π= + / ( )t t1 i i1 1υ = −+ Gaussian white noise with variable amplitude. The mass and stiffness of the building model are tuned to obtain a fundamental frequency equal to 3 Hz. The damping is set to 3%. Before the recording, more Gaussian white noise is added to simulate the presence of instrumental noise. Figure 2 reports the output of different sim- ulations with increasing noise, both at the base and at the top of the building. It is possible to see how the NonPaDAn method retrieves cor- rectly both damping and frequency if the noise is absent (maximum is the black dot inside the red area). The increase of input noise leaves the frequency unchanged but makes the damping peak broader. Nevertheless, the maximum in the 2D histograms shows the correct values for damping and fundamental frequency. On the contrary, sparse, spurious pseudo-frequencies arise when we add instrumental noise, and the damping is not correctly preserved. This simu- lation shows that the initial hypothesis is cor- rect: to obtain a reliable estimate it is necessary to keep at a minimum the instrumental noise. Fig. 1. SimuLink block model implemented to verify the assumptions underlying the NonPaDAn method. 262 Fig. 2. Effects of increasing input and instrumental noise on the simulated recording atop a SDOF building. Marco Mucciarelli and Maria Rosaria Gallipoli Fig. 3. NonPaDAn estimate of frequency and damping as a function of recording length (left: 30 min; right: 10 min). 263 Non-parametric analysis of a single seismometric recording to obtain building dynamic parameters 3. Empirical data analysis The first check performed on the reliability of NonPaDAn is relevant to the duration of the sig- nal. If using RandomDec, a minimum of 30 min of recording is required to obtain a stable esti- mate, according to Dunand et al. (2002). A 30 min recording at the 5th floor of the Engineering Faculty building at the University of Basilicata was performed using a digital tromometer (Mi- cromed Tromino). Then, the damping was esti- mated using both the whole record and a section of 10 min. Figure 3 shows the comparison: the fundamental frequency is the same and the damp- ing remains centred at around 3%. In this case it is also possible to identify a higher mode, but one can see in the following that it is not always so. The second check was about the compari- son with published results obtained using dif- ferent methods to estimate frequency and damping. Di Giulio et al. (2005) used the Horizontal to Vertical Spectral Ratio (HVSR) and the Hor- izontal to Horizontal Spectral Ratio (HHSR) technique to estimate the fundamental frequen- cy of some buildings in the town of Palermo, Italy. The recordings were collected with a Lennartz 5 s seismometer and a MarsLite digi- tiser. We analysed 30 min of ambient vibration they recorded at the 5th floor of the R.C. build- ing code-named PAL8. Figure 4 shows the re- sult, with a clear estimate of fundamental fre- quency (3.2 Hz) and damping (3.5%). The esti- mation of frequency is in a good agreement with the result given by Di Giulio et al. (2005). During 2005 a series of tests were conduct- ed on a R.C. scale model on the ENEA shaking table at the Casaccia (Rome) facility to study retrofit systems (Cardone et al., 2003). The in- put ranged from random noise to real strong motion. Cardone et al. (2005) estimated from clear transients a damping around 3% and a fundamental frequency of about 2 Hz. Figure 5 reports the results obtained with NonPaDAn, in good agreement with the estimate from Car- done et al. (2005). Fig. 4. NonPaDAn estimate of frequency and damping for a R.C. building (PAL8 in Di Giulio et al., 2004). Fig. 5. NonPaDAn estimate of frequency and damping for a scale model tested with random vibration at the ENEA Casaccia shaking table (Rome). 4 5 264 Marco Mucciarelli and Maria Rosaria Gallipoli The last check is relevant to the ability of NonPaDAn to estimate the variations of frequen- cy and damping due to large displacements, damage or changes in the structural characteris- tics. For the first two items, the recordings ob- tained during the 2002 Molise earthquake se- quence inside a building in Bonefro, Italy are re-analysed. The record comprises 1 min of am- bient vibration, the 1st November shock (Mw = = 5.6) that damaged the building and whose vi- bration lasted 2.5 min and finally 1 min of ambi- ent vibration. Mucciarelli et al. (2004) analysed this records using different techniques (Gabor transform, wavelets, moving window HVSR). The building has an initial frequency of 2.5 Hz. During the large oscillations and subsequent damage the frequency lowers to 1.5 Hz, to reach a final value for the damaged building at about 2 Hz. Figure 6 reports the result of NonPaDAn ob- tained for three 1 min windows: the first minute of ambient vibration, the coda of the shock and the last minute of ambient vibration. The fre- quency shift estimated with NonPaDAn agrees with the result of Mucciarelli et al. (2004). Moreover, it is worth noting that during the dam- aging large oscillations, a secondary, isolated damping emerges with a value of 10%. This is Fig. 7. NonPaDAn estimate of frequency and damping for a 4 stories R.C. Building in potenza, Italy with the bare frame only (left) and completed with infills (right). Fig. 6. NonPaDAn estimate of frequency and damping for a R.C. building in Bonefro, Italy, during the dam- age induced by the 1st November 2002 Molise earthquake. 265 Non-parametric analysis of a single seismometric recording to obtain building dynamic parameters 2.5 times the starting damping, in agreement with observation on damage in laboratory tests (see for example Cardone et al., 2004; and Dolce et al., 2004). The estimated damping factors are confirmed by a numerical modelling of the building, as described in Dolce et al. (2007). The reliability of NonPaDAn to estimate the damping variation due to different building char- acteristics was checked taking advantage of a 4 storey, R.C. building under construction in Po- tenza, Italy. The first measurement was taken when the bare frame was completed and the bricks for the infills were loaded at each floor. The second measurement took place after the in- fills were completed. In this way, a variation of the stiffness occurred while the total mass of the building remained constant. In both case 10 min of ambient vibration were sampled with a Mi- cromed digital tromometer. Figure 7 shows that NonPaDAn is able to capture the increase of fre- quency and the decrease of damping that was ex- pected as a consequence of the increase in stiff- ness. Using this last example, we take the proposed method a step further, suggesting a possible para- meterisation of the model and thus the possibility of performing quantitative statistical tests. The best way of observing the true damping of a structure is after forcing, when the free oscillations’ amplitude decreases. The ambient noise is a particular case of output-only, un- known-input system: the forcing is random in amplitude and phase, so the structure cannot ful- ly develop its damped free oscillation. Neverthe- less, an improved strategy is to look for damp- ing on more than two relative maxima separated by the same pseudo-period. These consecutive, decreasing local maxima are the best approxi- mation of damped free oscillation that ambient noise measurement can produce. Of course a compromise is needed, between the search for long sequences (that would provide the best, more stable result) and the fact that the number of available sequences decrease exponentially when increasing their length. If we select only the damping values calculated on three consec- utive maxima, the result is less dispersed and the size of the sample still permits to estimate their statistical distribution, which is log-normal. For the R.C. building with and without infills the plot of the two distributions (fig. 8) reveals that there is a constant shift of the damping by a fac- tor 1.3. Moreover, a Kolmogorov-Smirnov test allows rejecting the hypothesis that the two CDFs derive from the same distribution with a 97.5% confidence level. This approach could be useful for measurements perfomed before and after an earthquake to establish if there is a sig- nificant frequency and damping shift due to damage: the ratio between CDFs appears to be more stable than the comparison between two couples of values. 4. Conclusions We proposed a simple method for the non- parametric analysis of damping of buildings (NonPaDAn). It is able to capture the funda- mental mode and its associated damping with a short, single measurement of ambient vibration atop a building. The resolving power is less than that provided by more advanced tech- niques and the higher modes and associated damping are rarely visible. However, the advan- tage of the proposed method is related to the limited resources it requires in terms of time, instrumentation and manpower. This candidates NonPaDAn as a tool to be used for dense meas- urement campaigns, when it is more important Fig. 8. Comparison of CDFs of damping derived by more than two subsequent maxima for the same building of fig. 7. 266 Marco Mucciarelli and Maria Rosaria Gallipoli to know the fundamental dynamic parameters of large sets of buildings rather than focusing on the detailed behaviour of a single item. Acknowledgements This work was funded by NATO grant SfP 980857. Thanks are due to Antonio Rovelli and Giuseppe Di Giulio for providing recordings of PAL8, Felice Carlo Ponzo for providing record- ings of ENEA test and Mauro Dolce for com- ments and suggestions. REFERENCES CARDONE D., G. DE CANIO, M. DOLCE, R. MARNETTO, C. MORONI, M. NICOLETTI, D. NIGRO, A. PIZZARI, F.C. PONZO, E. RENZI, G. SANTARSIERO and D. SPINA (2003): Comparison of different passive control techniques through shaking table tests, in Proceedings of the 8th World Seminar on Seismic Isolation, Energy Dissipa- tion and Active Vibration Control of Structures, Octo- ber 6-10, 2003, Yerevan, Armenia. CARDONE, D., E. COELHO, M. DOLCE and F. PONZO (2004): Experimental behaviour of R/C frames retrofitted with dissipating and re-centring braces, J. Earthq. Eng., 8 (3), 361-396. CARDONE, D., A. DI CESARE, M. DOLCE, C. MORONI, D. NI- GRO, F.C. PONZO and M. NICOLETTI (2005): Compari- son of shaking table and pseudodynamic test on a base isolated R/C structure, in Proceedings of the 9th World Seminar on Seismic Isolation, Energy Dissipation and Active Vibration Control of Structures, June 13-16, 2005, Kobe, Japan. CLOUGH, R.W. and J. PENZIEN (1975): Dynamics of Struc- tures (Mac Graw-Hill, New York), 2nd edition, pp. 576. DI GIULIO, G., R M. AZZARA, G. CULTRERA, M.S. GI- AMMARINARO, P. VALLONE and A. ROVELLI (2005): Ef- fect of local geology on ground motion in the city of Palermo, Italy, as inferred from aftershocks of the 6 September 2002 Mw 5.9 earthquake, Bull. Seismol. Soc. Am., 95, 2328-2341. DOLCE, M., D. CARDONE, F.C. PONZO and C. VALENTE (2004): Shaking table tests on reinforced concrete frames without and with passive control systems, Earthq. Eng. Struct. Dyn., 34, 1687-1717. DOLCE, M., F.C. PONZO, A. MASI, M. MUCCIARELLI, M.R. GALLIPOLI, A. DI CESARE, M. TETTA and M. VONA (2007): Identification of the structural model and analysis of the global seismic behaviour of the IACP Bonefro building, Bull. Earthq. Eng. (submitted). DUNAND, F., P.-Y. BARD, J. L. CHATELAIN, PH. GUÉGUEN, T. VASSAIL and M.N. FARSI (2002): Damping and fre- quency from RandomDec method applied to in situ measurement of ambient vibrations: evidence for effec- tive soil structure interaction, in Proceedings of the 12th European Conference on Earth Engineering, Lon- don, U.K., CD-ROM edition, Pap. no. 869. DUVAL, A.-M., J.-L. CHATELAIN, B. GUILLIER and SESAME PROJECT WP02 TEAM (2004): Influence of experimen- tal conditions on H/V determination using ambient vi- brations (noise), in Proceedings of the 13th World Con- ference on Earthquake Engineering, August 1-6, 2004, Vancouver, B.C., Canada, Pap. no. 306. HUANG, C.S., Y.B. YANG, L.Y. LU and C.H. CHEN (1999): Dy- namic testing and system identification of a multi-span highway bridge, Earthq. Eng. Struct. Dyn., 28, 857-878. LAMARQUE, C.-H., S. PERNOT and A. CUER (2000): Damp- ing identification in multi-degree-of-freedom systems via a wavelet-logarithmic decrement, Part 1. Theory, J. Sound Vib., 235, 361-374. MUCCIARELLI, M. (1998): Reliability and applicability range of the Nakamura’s technique, J. Earthq. Eng., 2 (4), 625- 638. MUCCIARELLI, M., A. MASI, M.R. GALLIPOLI, P. HARABA- GLIA, M. VONA, F. PONZO and M. DOLCE (2004): Analy- sis of a R.C. building dynamic response and soil-build- ing resonance based on data recorded during a damag- ing earthquake (Molise, Italy 2002), Bull. Seismol. Soc. Am., 94 (5), 1943-1953. MUCCIARELLI, M., M. R. GALLIPOLI, D. DI GIACOMO, F. DI NOTA and E. NINO (2005): Wind influence on seismic noise measurements, Geophys. J. Int., 161, 303-308. NAVARRO, M., F. VIDAL, M. FERICHE, T. ENOMOTO, F.J. SÁNCHEZ and I. MATSUDA (2004): Expected ground-RC building structures resonance phenomena in Granada city (Southern Spain), in Proceedings of the 13th World Conference on Earthquake Engineering, August 1-6, 2004, Vancouver, B.C., Canada, Pap. no. 3308. VANDIVER, J.K., A.B. DUNWOODY, R.B. CAMPBELL and M F. COOK (1982): A mathematical basis for the random decrement vibration signature analysis technique, J. Mech. Des., 104, 307-313. (received October 27, 2006; accepted April 4, 2007)