|
|
|||
|
|
|
|
|
| Unfactored Dead Load (tonnes) |
|
|
|
| No. of Bearings |
|
|
|
| Size (mm) |
|
|
|
| Expected Mean Deflection (mm) |
|
|
|
A popular method of dynamic testing involves hitting the structure with a hammer. Clearly, the magnitude of impact is a function of the mass of the hammer, and its velocity when it hits the structure. Variability in the magnitude of impact is not considered a problem, as the results are normalised. There arises a system frequency characterised by the stiffness of the contacting surface, and the mass of the hammerhead, above which it is difficult to deliver energy into the structure (Ewins, 1984). The impact should ideally produce a half sine pulse, the duration of which controls this upper frequency. The stiffer the contact element, the shorter the pulse duration, which leads to a higher frequency range that is excited. However, lower frequencies were of greater interest for the testing of this structure.
A large impulse was desired, to provide an adequate signal at remote measurement points. For a given impacter mass, it is possible to increase the impulse by increasing the drop height. But this invariably causes a rebound, which should be kept well away from the first impulse, to avoid signal processing difficulties. A suitable method has been devised by Hunt (1988), where a weight with a rubber tip is dropped on an arm from a given height, where the arm is pivoted at one point. This method is however limited in its ability to excite the structure adequately at low frequencies (Cryer, 1994).
The Author wanted to investigate the possibility of using a larger mass, falling a smaller height, but on a layer of sand, as the contacting resilient element. Such a method has been used by the Author to excite bridge decks (Sharif, 1996b). This method of using a sand layer in place of a rubber pad would eliminate a rebound, and should be adequate to excite the low frequency modes.
The impact method used here meant that the mass would fall on a yielding layer of sand. The energy is imparted both into the structure, and work done in compacting the layer of sand. It was anticipated that the layer of sand would vary between impacts, and effect the delivery of high frequency energy. But it was interesting to examine if this variation between impacts would alter the low frequency response that was of interest.
It was desirable to arrange the impacts on the pile cap to assess transmissibility across the isolators. But this was not possible, as the void between the pile cap and the isolated building was too small to deploy any impact testing equipment. Therefore impacts could only be applied above the isolators, from within the basement car park.
An impacter mass of 150 kg was used and made up of ten steel plates, 300mm x 300mm x 21mm thick which were bolted together on site. The impacter mass was suspended in position via a rope off a mobile crane. A bag of sand was laid flat in the desired position of impact to give a nominal thickness after initial impacts of 50mm. Plate 8.5 shows a set up of the steel impacter mass with the sandbag in position at the base of isolated column B15. The suspension rope was severed using a hand cutter.
The drop height was varied on site until adequate signals were recorded at all measurement locations. This drop height (nominally 200mm) was then repeated throughout the tests.
8.13.1 Impact Tests - Signal processing
It was necessary to obtain the force of impact, to calculate the 'point' inertance, defined as the ratio of acceleration response at a given point with the applied force at the same point. The 'transfer' inertance is obtained when the response at a given point is divided by a force applied elsewhere, for example, the response of the top of the column to a force applied at the base.
The force of impact was to be deduced using accelerometer recordings taken from the centre of the impacter mass. The deceleration recorded, together with the impacter mass, would give force from Newton's Second Law, and Newtons Third Law would give the force imparted to the structure. This would provide a very simple means of deducing force (also used by Cryer, 1994), but problems were anticipated.
Diagram 8.1 shows a schematic of a Brüel & Kjær DeltaShear® accelerometer. It was recognised that under static conditions the seismic mass would be in equilibrium by shear forces generated by the piezoelectric element.
Schematic of B&K DeltaShear® accelerometer (B&K Vibration handbook) Diagram 8.1
However, during free fall the body of the accelerometer case would fall away from the seismic mass due to its inertia, leading to a shear force on the piezoelectric material. A charge is thus produced, and reflects the force required to accelerate the seismic mass by 'g' (9.81 m/s2). This would result in a step in the acceleration output.
An ideal impulse would be a half sine wave. A dc offset and low frequency trend would cause errors in the auto spectra used to calculate inertance. Unfortunately due to cost restrictions this was the only method of force detection available to the Author. It will be apparent from the subsequent discussion that this financial restriction proved to compromise the experiment in terms of the accuracy of the force spectrum. The impulse varied in magnitude of force, and the period varied from 10msec to 20msec. This variation in period is important as it affects the upper frequency at which energy can be imparted to the structure. However, as only axial modes up to 50Hz were of primary interest, this method of impact testing was considered acceptable.
Figure 8.34(b) shows a time history of force, which shows a dc offset at the start and immediately after the impulse. The auto spectrum of force could be obtained using a hanning window of 2 seconds, centred on the impulse, and zero padded to 12 seconds, to reduce the effect of noise in the part of the time history that was of no interest. It was felt that this window would not distort the impulse or structures' response, when it was centred on the event and of 2 seconds length (used for Figures 8.45 and 8.46 which compare vibration measurements). However this method distorted the force auto spectrum (see Sharif, 1999). An exponential window was therefore used, as it had the advantage of reducing side lobes in Nyquist plots (see Chapter 5; Mercer, Turner and Kandianis, 1969; Bhaskar, 1992), and in this case would also eliminate the dc offset in force measurement at start of impulse, and weight the signal in a way that reduced the contribution of noise downstream of the impulse. The original and windowed force and vibration measurements are shown in Figure 8.34.
The shape of the impulse obtained (Figure 8.34(b)) was not in keeping with the ideal half sine pulse shown in Figure 8.35(a). There is a sloping edge and a sharper peak to the measured impulse, due to the non-linear character of loading the sand. Cryer (1994) has also observed strong ringing after such impulses, known as Gibb's phenomena.
Figure 8.35(c) shows a theoretical auto spectrum of force for the idealised impulse (Figure 8.35(a)), along with the result for exponential window applied to measured data in Figure 8.34(b). There is clearly still some distortion at low frequencies, which could be eliminated by using a boxcar window that centres on the impulse, and that matches the impulse period, as shown in Figure 8.35(b). However it would not be appropriate to apply this same short boxcar window on the vibration measurements of the impulse.
It was therefore decided for ease and consistency to apply the exponential window for impact test results involving both force and vibration signals
We can note that the accelerometer could not be placed directly on the centroid of the impacter, as it would clash with the hook used for suspending it. This would therefore cause a small error, as the deceleration would include rocking components on impact.
However, despite these points, it is interesting to examine the results, although the force results and those derived from it are to be interpreted with some caution.
8.13.2 Impact Tests - Isolated column B15 / Un-isolated column G15
Impact tests were undertaken on isolated column B15 and the identical but un-isolated column G15.
The impacter steel block was positioned directly adjacent to column B15, for impacts on the basement slab, representing the base of the isolated column. Figure 8.36 shows the impulse, along with the structure's response. Trace (f) shows the force time history, with traces (b), (c) and (d) giving the structure's response on the same vertical scale. Magnified vertical scales for the top of the isolated column and pile cap responses are shown in traces (a) and (e). We can see that the impacter mass of 150kg with a 200mm drop height produces in this example a peak force of 52kN with a pulse duration of some 10msec. This impulse produced an adequate response signal above background levels at all locations. Clearly building response measurements do not show any step seen in the force measurements, and therefore any spectral analysis of these is reliable.
There exists a delay in the response of the top of the column compared to the base at which the impulse is applied. From this delay and the distance between the two sensors (25.5m), we can simply determine the compression wave propagation velocity as 3336 m/sec. On other columns this varied down to 2000 m/sec, representing a large range.
Impacts were again applied on the basement slab at the base of the un-isolated column G15, using the same drop height of 200mm.
The ratio of the auto spectra estimates (accelerometer response / input force) were used to derive inertance. The 'point' inertance of the base and 'transfer' inertance of the top of the isolated column and un-isolated column are shown in Figure 8.37. We can see that it is easier to excite the base of the column, when applying the force at the base. This is not surprising as the impulse at the base is attenuated three dimensionally by the structure as it travels up the building.
These results also show that an isolated column is easier to excite than its un-isolated counterpart. This has implications both to external sources of vibration, but also to internally generated forces. Therefore, the evidence that the isolated structure shakes significantly more than the un-isolated structure, due to building services plant under background conditions, shown earlier (Figures 8.29), is to be expected.
This difference in behaviour of isolated and un-isolated columns is likely to be due to a combination of factors. The impulse more easily excites the base of the isolated column which is relatively 'free', whereas it is harder to excite the base of the un-isolated column which is more effectively restrained by the stiffness of the pile. The impulse energy should also radiate into the soil more readily for the un-isolated column, whereas this path should be more inhibited by the isolators.
The Nyquist plots of the cross spectra between force and structure response were processed to allow the identification of natural frequencies and circle fits used to establish modal damping ratios, as described in Chapter 5. Figure 8.38 shows an example of the point mobility of the isolated column together with estimates of modal damping. Table 8.2 summarises a selection of modal damping ratios obtained for the isolated and un-isolated columns from an impact test.
Table 8.2 Comparison of modal damping ratios from circle fits of impact tests
|
|
|
|||
|
(Hz) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|
|
|
|
|
The damping values are on the high side at low frequencies, although not unrealistic. We can see that the damping ratios reduce for higher natural frequencies. The values in the isolated and un-isolated columns show no dramatic differences, although we cannot make direct comparisons for each value as they can pertain to different modes. The trend of a reduction in critical damping ratio with increasing frequency would imply that a hysteretic type damping model would be appropriate (see chapter 5).
8.14 Shaker Tests
Another technique to identify the natural modes of a structure is to observe the steady state response under harmonic excitation, across the frequency range of interest. This technique was used to excite the ground to understand foundation dynamics, as early as 1936 by the research organisation Degebo based in Germany. The results of similar early work are described in full by Lorenz (1960).
There are a number of ways to apply this harmonic force in such applications. The earlier systems utilised rotating eccentric masses, where force is simply proportional to the frequency squared. A disadvantage with this method is that non-linear effects can arise under different levels of strain caused by the large variation in force. Later systems utilise electrodynamic (Grootenhuis & Awojobi, 1966) or hydraulic shakers, where force can be controlled independent of frequency. Another advantage with these two methods is that they can provide any shaped forcing function. Although, for the large forces required to drive a structure, such systems are often more expensive than the simpler rotating eccentric mass exciters.
Rather than building such a unit from scratch, the Author decided to modify an existing device that is readily available to the building industry. This consists of a 'wacker' plate, used for compacting hardcore and road sub-bases. The standard unit consists of a steel bedplate on which is rigidly mounted an eccentric mass on a shaft, which is driven at a constant speed by a petrol engine, via a drive belt.
The Author replaced the engine with a dc motor, where the armature voltage is proportional to frequency and can therefore be adjusted by a dc motor controller. The dc controller was also modified to provide accurate control of frequency using a ten-turn micrometer reading potentiometer. This allowed frequency to be adjusted in 1Hz steps from 4Hz to 50Hz.
The frequency range of the shaker was limited to 4Hz at the low end, below which the forces were too small, and limited to 50Hz at the high end, above which the forces would have been too large.
Appendix 8.1 describes the shaker in more detail along with the simple test to establish the delivered force as a function of frequency. This force was derived by measuring the acceleration of the shaker, when suspended in air. From the mass of the shaker, forces giving rise to the measured accelerations were simple to deduce. The test showed the unit would deliver a quadratic force reaching 7200N at 50Hz.
The unit as conceived, produced a rotating vector force, with vertical and horizontal components. It is possible to eliminate the horizontal component by arranging for counter rotating masses, which balance the horizontal components, but amplify vertical ones. This is the principle of operation that was used in the exciters described by Lorenz (1960), and used as the basis of design in vibrators used for driving piles. However this would have required detailed modifications not considered essential.
Attaching the unit on site did present a problem, as it needed to be achieved without introducing any permanent fixings, yet given the large forces that arise both vertically and horizontally, a strong method of attachment was required.
A simple method of locating the shaker without interfering with the structure in any permanent way was devised. This entailed using three acro props, which were expanded in position, with one end pushing the baseplate of the shaker down onto the basement floor, whilst the other end derived a reaction off a concrete beam supporting the floor above.
The shaker was deployed at the base of three columns. On the isolated columns A15 and B15, to compare similar isolated columns, and the un-isolated column G15, which is identical in every other respect to B15.
Plate 8.6 shows the shaker in position at the base of isolated column B15.
8.14.1 Shaker Tests - Signal Processing
Figure 8.39 shows time histories for measurements on the pile cap, base and top of isolated column B15. The shaker is operating at 50Hz. The driving force applied at the base is attenuated significantly across the isolators to the pile cap (note different scales), but more so up to the top of the column.
There are higher frequencies evident at the base of the isolated column and pile cap, which is due to bearing noise of the shaft. There is also an indirect component coming from the flexing of the acro props, which are put into transverse motion in varying degrees, as the tests sweep across their natural frequencies.
The complete run of shaker tests was very time consuming. The charge amplifiers frequently required their gain settings to be adjusted to optimise the signal, as did the input range on the PC based monitoring system, to optimise quantisation of the signal. At each frequency, the test was sustained for at least 45 seconds, from which the end section was recorded to obtain steady state conditions. The tests continued in approximately 1Hz steps, all the way down to 4Hz.
It is important that the vibration signals from the shaker tests are significantly above background levels. Yet at low speeds, the forces from the shaker are much smaller (46N), and therefore the signals were less distinct against background levels. As a result, the test results should only be scrutinised down to 10Hz, even though the tests were undertaken down to 4Hz.
The post-processing was automated and undertaken in two ways, using programs written by the Author in MATLAB®.
One method entailed calculating inertance of the pile cap, base and top of the column from a ratio of the peak response at the driving frequency using auto spectrum estimates, divided by the respective force deduced for the shaker at each test frequency.
A simple search algorithm was devised to obtain magnitude of the peak at the relevant frequency for each test, as otherwise a peak from a neighbouring frequency could be mistaken for the forced signal. However, despite the use of this method of analysis, the inertance so derived was unable to resolve the detailed character of response at lower frequencies, owing to dominance of noise in neighbouring frequencies.
The alternative method used and presented here, entailed applying a filter which would selectively bandpass the driving frequency. The peak acceleration at various measurement points was then simply divided by the peak harmonic force at each test frequency.
The digital filtering could not be achieved using a single bandpass filter, as the high sampling rate used prevented a narrow and accurate bandpass to be designed. This could be overcome by resampling the data, at a lower rate, which would be sufficient to characterise the low frequencies. However, it was decided to achieve the bandpass filter on the original data by in effect applying two filters sequentially, a low pass and a high pass filter. This would provide the desired accuracy of 0.1dB ripple in the passband and (- 50dB) in the stopband.
An Elliptic filter was used as it provides a minimum transition width compared with the conventional Butterworth filter. Filtering causes phase distortion, and these were minimised by filtering the original sequence and its time reversed version, which averages out these phase errors, as described by Krauss et all (1994).
The harmonic response testing also allows us to estimate the degree of damping using the half power bandwidth method described in Chapter 5. This method can be used to estimate modal damping for any well separated mode.
8.14.2 Shaker Tests - Isolated Columns A15 & B15 / Un-isolated Column G15
Figure 8.40 compares the 'point' inertance of the base of isolated columns A15 and B15, which shows that it is easier to excite the outer isolated column A15, given an excitation at its base. The 'transfer' inertance is also shown to the top of isolated columns showing the same trend. Shaker tests on isolated column A15 were not extended to the full 50Hz range owing to some initial concern over the high levels of vibration that arose.
The fact that it is easier to excite the outer column A15 compared to the inner column B15, should not be surprising when we consider that column A15 supports about half the mass of column B15, by virtue of it carrying a smaller area of floor loading.
We can now uniquely compare the inertance for isolated (B15) and un-isolated (G15) columns that are otherwise identical. Figure 8.41 shows the 'point' inertance at the base of the columns, where it is seen to be clearly easier to excite an isolated column. This figure also shows a similar comparison for the 'transfer' inertance to the top, which shows that it is easier to excite the top of an isolated column at some frequencies, although at 15Hz, 26Hz and approaching 50Hz, the inertance of both isolated and un-isolated columns are almost identical. This may be due to axial modes that are common to both isolated and un-isolated columns.
We can use the half power bandwidth method to estimate modal damping ratios for the peaks obtained in the response of the top of the column, as described in Chapter 5. Where the peak response is not fully defined, we can by eye estimate straight line extensions. Some of the peaks in the response of the top of the columns and their respective damping values are summarised in Table 8.3.
Table 8.3 Peak response and damping for isolated and un-isolated columns
|
Isolated
Column B15 |
Isolated
Column A15 |
Un-isolated
Column G15 |
||||
|
Frequency of selected peaks at top of column
|
23 Hz
|
35 Hz
|
16Hz
|
22 Hz
|
27 Hz
|
43 Hz
|
|
Damping for the mode
|
4%
|
9%
|
8%
|
6%
|
5%
|
12 %
|
The sample of selected peaks shown in Table 8.3 are only an example where the modes were separated enough to make an estimate using the half power bandwidth method, although many other peaks in the response did arise.
Even so, the damping values were obtained for modes that are not all ideally independent, but are to a degree affected by nearby modes. The frequency testing in nominal 1Hz steps also did not resolve the peaks clearly. This makes use of the half power bandwidth method for estimating damping in these cases unreliable, even though the values obtained are not unrealistic.
A more reliable method of estimating modal damping, where modes are close by, is to plot the real and imaginary components of the complex response at equal frequency intervals. This method was described in Chapter 5 and used to obtain damping estimates from the impact tests. It could not be applied in this test sequence as the frequency intervals were nominally selected around 1Hz, controlled manually, and therefore not equally spaced.
8.15 'Point' Inertance of Pile Cap supporting Isolated and Un-isolated Columns
It was desirable to compare the 'point' inertance of the pile caps beneath an isolated column with that of an un-isolated column. This would provide a unique insight into how the interaction of the supported mass on a pile cap is affected by base isolation.
Comparisons of inertance using a ratio of pile cap response with the shaker applied force would not be meaningful. As in the case of the pile cap beneath the isolated column it would be the 'transfer' inertance, whereas in the case of the un-isolated column, it would be the 'point' inertance. This is because in the case of the isolated column, the force from the shaker is not the actual force applied on the pile cap below the isolator.
To obtain the pile cap force in an indirect way, the relative displacement (deduced from accelerometer measurements) between the pile cap and the base of the isolated column, was multiplied by the group isolator dynamic stiffness.
The group dynamic stiffness of the rubber bearings was simply taken from production load-deflection test curves. The tangent stiffness was calculated at the working load and multiplied by an empirical factor of 1.8, provided by the manufacturer (BTR) to deduce the dynamic stiffness. This is intended to allow for the higher stiffness that would result from the low strains expected in service (i.e. small strain cycling about a point on the overall load deflection curve leads to a steeper gradient of the stress strain loop c.f. tangent stiffness). Dynamic stiffness of natural rubber bearings loaded with carbon black are reported to be more sensitive to strain amplitude than frequency of loading (see Payne, 1962; Ahmadi & Muhr, 1999). Whilst dynamic stiffness increase slightly with frequency, it reduces with increasing strain amplitude. Sensitivity is dependent upon amount of carbon black used, being larger at higher proportions of carbon black.
Whilst dynamic tests on rubber bearings are simple to perform, the manufacturers did not have a large enough test rig to establish the actual dynamic stiffness of the bearings used. In the experiment with the shaker, strain increases as frequency is increased, and therefore would tend to reduce dynamic stiffness, which would therefore offset any increase in stiffness that might arise with increasing frequency of the loading. Given that we cannot quantify these small effects, it becomes an experimental assumption that the dynamic stiffness deduced from an empirical correction of the tangent to the load deflection curve, at the working load, is constant over the small frequency range (5Hz to 50Hz) of the test, and appropriate to the small level of strain imposed (<0.002%).
The 'point' inertance for the pile cap of the isolated column is therefore only considered an estimate, as it assumes that the dynamic stiffness of the rubber bearings is both correctly estimated, and constant in the frequency range of the test.
Figure 8.42 compares the 'point' inertance of the pile cap below isolated columns A15 and B15. The two curves are similar, although it is easier to excite the pile cap below column A15. This is not surprising as the pile cap for column A15 is formed using only one pile and carries a smaller load, whereas column B15 is heavier, and supported on a pile cap with two piles.
Figure 8.43 shows the unique comparison for 'point' inertance of the pile cap below the otherwise identical, but isolated and un-isolated columns. There is a distinct difference, increasing with frequency, whereby it is easier to excite the pile cap below the isolated column. At some low frequencies, the 'point' inertance is similar, where the presence of the isolators has little or no effect. But particularly at higher frequencies, the decoupling effect of the isolators increases. This indicates that the isolators decouple the inertia/stiffness of the supported column from the pile cap, and this decoupling is frequency dependent.
The transmissibilities for train events, reported in preceding sections, between pile cap and isolated column tells us just that, and not the insertion loss for base isolation. The insertion loss for base isolation is that benefit or disbenefit that arises between the isolated and un-isolated situation. If the pile cap beneath an isolated column were to respond more than the equivalent pile cap beneath the un-isolated column, then clearly an adjustment is required on the transmissibility results, if they are to be used to properly indicate insertion loss.
Figure 8.44(a) shows that for an applied dynamic force, the pile cap below the isolated column, would respond by in general 5dB to 10dB more than the pile cap below the un-isolated column. This result is significant, because whilst we know that the transmissibility deduced from pile cap to isolated column does not yield actual insertion loss, we can for the first time realise the size of the correction needed to more meaningfully assess the actual insertion loss from such transmissibility measurements.
This indicative adjustment has been estimated for a harmonic loading case, over a small frequency range and is shown along with the mean total transmissibility curve for isolated column B15 in Figure 8.44(b). The actual correction necessary will be more complicated, as the rail source is clearly not simply harmonic. However, an upward adjustment of 5dB to 10dB would clearly reduce the little evidence of isolation seen from such transmissibility curves, and would then show similarity to the comparative response of neighbouring isolated and un-isolated structures seen in Figure 8.28.
8.16 Comparison of Shaker and Impact Testing
We can compare the total transmissibilities from the base to the top of the isolated and un-isolated columns obtained under both shaker and impact testing.
Figure 8.45 shows a comparison for isolated column B15. Clearly we can only compare results up to the shaker test range of 50Hz, which shows a great deal of similarity with the impact testing. The size and position of the peaks are all similar, except that there is a greater divergence at 50Hz. This may in part be due to the high level of background vibration at 50Hz, which may not have been adequately exceeded by the impact test, whereas was in the case of the shaker results. Also, 50Hz was potentially a resonant mode of the isolated column (free-fixed), which under steady state conditions imposed by the shaker produces more pronounced response than did the impact test. Similar correlation between the two methods of testing arises in the un-isolated column results shown in Figure 8.46.
This confirms that the impact testing method used is acceptable for mode characterisation, is quicker to implement and therefore a very desirable technique. However, there was some variability in the results (see Sharif, 1999) and so it is important that a number of impact tests are sampled.
Even though this impact testing method has been validated, it cannot be used to assess the performance of a large base isolated building due to a train source. This is because the impacts cause vibration sources that are local and this decays three dimensionally into the structure. The rail source would in fact excite all the columns in varying degrees and with varying correlation of inputs (see section 8.10). The three dimensional decay of input from one pile under train vibration excitation would be offset by vibration coming from three dimensional decay from neighbouring columns as they shed energy through the structure. Such a return of energy from neighbouring columns would not arise as well in the impact testing case, and so would not represent a realistic simulation of a rail source.
8.17 Theoretical Analysis and Comparison with Measurements
This section describes theoretical analysis of the axial modes of a column using Finite Elements, which were modelled using the Ansys program version 5.2, running on a Pentium Pro computer. The modal characteristics obtained from simple bar models, lumped mass column models are compared with measurements.
The finite stiffness of the piled foundation was also considered, and calculated using the Dyna4 program, which is based upon the methods due to Novak et al (1974,1978a,1978b,1980,1993).
This foundation analysis program was also used to calculate the theoretical point inertance of the piled foundations, which were treated as unloaded and loaded. These were then compared with measured values.
8.17.1 Simple Bar Models of Column
The concrete column modelled as a simple bar, with fixed support produces a fixed-free mode with a frequency of 34Hz. With resilient supports for the isolated case (target rigid body natural frequency of 7.5Hz) it produces a free-free mode at 69Hz.
It will be recalled that site column measurements showed many resonance peaks in the 100Hz range. Whilst we would not expect to model the 'real' three dimensional structure using such a simple column model, we can significantly improve this bar representation by modelling the column with a series of lumped masses at each storey height to represent the inertia of transverse beams, floors and finishes.
It is recognised that lumping all such dead loads onto the column, storey by storey, may be unrealistic. During column oscillation the entire mass of floors and beams may not be mobilised and is likely to be frequency (wavelength) dependent. Nonetheless, it is considered more realistic than just modelling the column as a bar.
8.17.2 Concrete Column - Lumped Masses - (30m length)
The concrete column was therefore modelled allowing an appropriate lumped mass per storey height (25.5 tonnes for columns on Gridline A, 51.18 tonnes for columns on Gridline B and 39 tonnes for columns on Gridline C), representing the unfactored dead loads. The concrete column was also modelled to represent its varying cross section, which becomes more slender with height. Cross sectional area is in this case relevant, as it affects the columns' axial stiffness, storey by storey under varying load conditions.
A fixed-free mode occurs at 10.9Hz for columns on Gridline A, 7.0Hz for columns on Gridline B and at 8.3Hz for columns on Gridline C, which compares with 34.0Hz for the column modelled as a bar without any lumped masses. It may be surprising that the first mode axial frequency of the fixed base column with lumped masses is so low, and in the region often targeted for rigid body natural frequencies in base isolation. For the same columns, but with an isolated base set to a rigid body vertical natural frequency of 7.5Hz, the first mode of the column is lowered to 6.0Hz for columns on gridline A, 5.3Hz for columns on Gridline B and 5.6Hz for columns on gridline C.
For the isolated case there now also exists a "free-free" mode at 19.9Hz for columns on gridline A, 15.2Hz for columns on gridline B and 17.6Hz for columns on gridline C. This compares with a free-free mode at 69Hz when the column is modelled as a simple bar. Free-free modes should be important, as the anti-node occurs at the point of excitation, and therefore are easy to excite.
An indication of the variation in axial modes that arise for different columns in the same building, due to different section sizes and storey loading is summarised in Table 8.4.
Table 8.4 Axial Modes for Columns on Gridlines "A", "B" & "C"
|
Un-isolated column
|
|
||
|
Grid "A"
|
|
|
|
|
Grid "B"
|
|
|
|
|
Grid "C"
|
|
|
|
|
Isolated column
|
|
||
|
Grid "A"
|
|
|
|
|
Grid "B"
|
|
|
|
|
Grid "C"
|
|
|
|
Differences vary and depend upon which mode is being compared. It will be shown that these differences are significant in limiting isolation obtained for a given column.
Whilst single column models, with lumped mass loading per storey height are instructive, spatially extended models are still necessary. This is because they can allow the correct inclusion of mass/stiffness per storey height from floors and transverse beams, and would model transfer of vibration from one column and another.
8.17.3 Sensitivity of 'free-free' Axial Column Modes to Isolator Support Stiffness
The natural frequency of a 'free-free' mode of a column will be affected by the magnitude of the support stiffness. We can examine this by considering three values of isolator support stiffness. The dynamic stiffness of the group of isolators under the relevant column was used in the isolated cases so far described, but two softer situations were considered, at 1/10th and 1/100th the actual stiffness values. These values do not represent any practical variation that would ever occur, but are used to assess how the isolators influence the 'free-free' modal characteristics of the column.
The 'free-free' modes vary as a function of isolator support stiffness, which are summarised in Table 8.5.
Table 8.5 Sensitivity 'free-free' column natural frequency to support stiffness (kv)
|
|
|
|
|
|
|
|
|
|
|
|
From this we can see that the 'free-free' mode established using the actual isolator stiffness is not the actual 'free-free' mode of the column were it to be suspended freely in space. This would be closer to a value of 18Hz as seen by the softer mounts case. Clearly the isolators permit a 'free-free' mode to exist, but their stiffness effects the natural frequency obtained.
It is therefore interesting to examine the 'fixed-free' case again, but this time using the finite stiffness provided by piled foundations, rather than assuming a fully fixed base. To enable this analysis, a dynamic analysis of the pile foundation was required.
8.17.4 Sub-Structure model: Pile / Soil System
The Dyna4 computer program was used to model the first axial mode of a lumped mass on a single pile, in a layered soil. A pile of 1m diameter, 18m length was modelled, with its end bearing conditions resembling the chalk stratum into which the piles were toed. Soil dynamic properties should ideally be based upon site tests, which usually involve cross-hole seismic techniques. However, in this case values were assumed based upon the Author's experience of other sites. An arbitrary vertical force was applied across a broad frequency range to identify the first axial mode. All the properties assumed for this model are shown in Figure 8.47.
The natural frequency of the pile loaded with a lump mass (appropriate to the total foundation load at pile cap level) is 19Hz, which is within the range of dominant frequency spectra caused by trains, and also close to the axial modes of the isolated column.
The finite element model of the column with appropriate lump masses per storey height was then modified by replacing the fixed base assumption of the column with the finite stiffness for the piled foundation, taken from this Dyna4 model, and implemented using a simple spring element.
Whilst stiffness and damping are frequency dependent, a single value for stiffness at 20Hz was used and taken to be 3 x 109 N/m. The piled foundation is therefore at this frequency seven times the dynamic stiffness of the group of rubber bearings used beneath this column.
A column on gridline A was modelled with a piled base, and had the first mode at 9.6Hz which is broadly similar to the first mode for the fixed base at 10.9Hz. However the second mode, with piled base is 24.9Hz compared with 29.3Hz for the fixed base.
Therefore a piled based alters the axial modes from the fully fixed base assumption. A 'free-free' mode of the column can now begin to emerge, as the pile introduces resilience and therefore a degree of freedom at the base of the column.
8.17.5 Theoretical 'Point' Inertance of Pile Cap for Isolated and Un-Isolated Columns
We can also use the Dyna4 program to calculate the 'point' inertance of the pile cap, which we can model as unloaded and loaded.
Piled foundations were modelled in such a way as to match those which were tested, namely; isolated column pile caps A15 & B15, and the un-isolated column pile cap G15. These are schematically shown in Figure 8.47, along with the concrete and soil properties that were assumed.
The pile cap below isolated column A15 is provided with a single pile, which is loaded concentrically. The pile caps beneath isolated column B15 and un-isolated column G15 are in fact identical, being rectangular and served by two piles. The column loads are applied near the centroid of the rectangular pile cap, but these are necessarily eccentric to the two piles below.
Figure 8.48 compares two versions of the theoretical point inertance of the pile cap beneath isolated column A15, along with the measured result. One of the theoretical predictions assumes that the pile cap is unloaded (solid line), whilst the other assumes that the total mass of the column is lumped as a point load on the pile cap (dashed line). We can see that the theoretical prediction, which assumes that the pile cap is unloaded, matches the measured results.
This is an important result, as it indicates that the pile cap beneath an isolated column does not mobilise the inertia of the supported building. This arises because the isolators dynamically detach the inertia of the building from the motion of the pile cap.
We can examine if this correspondence also arises for the pile cap beneath isolated column B15, shown in Figure 8.49. When we model this rectangular pile cap served by two piles as being rigid, as we did in the case of the pile cap served by a single pile for column A15, we find there is poor correlation with the measured result. However, when we model the rectangular pile cap as flexible, we find that the correspondence is greatly improved. There are two curves which straddle the measurements, and refer to the inertance at different points of the flexible pile cap, one of which corresponds better with the measurement position used on site (see footnote to Figure 8.49).
This is a significant result, because firstly it tells us that even pile caps, which are traditionally thought of as rigid entities, must clearly be treated as elastic when multiple piles are involved, or where a pile or piles are not loaded concentrically. Secondly, the result again confirms that the measured inertance follows the theoretical inertance, when the pile cap is assumed to be unloaded. Therefore, the interaction of the mass of the supported building with the pile cap is affected by the isolators.
We can recall from Figure 8.42 that the measured inertance of the pile cap below column A15 was slightly higher than the inertance of the pile cap below column B15. This small difference was also borne out by the theoretical inertance of the two unloaded pile caps (see Sharif, 1999).
We can compare the theoretical inertance of an unloaded and loaded pile cap for columns B15 and G15 along with the measured results in Figure 8.50. We see as before, that the pile cap below the isolated column closely follows the theoretical inertance where the pile cap is assumed to be unloaded. Significantly, we see the measured inertance of the pile cap for the un-isolated column becoming increasingly divergent from this prediction as frequency increases.
The theoretical prediction of inertance for the loaded pile cap is also shown. Whilst the measured inertance for the pile cap below the un-isolated column does not follow this prediction either, it does lie between the two theoretical curves. This implies that the inertance of a pile cap below an un-isolated column is affected by the inertia of the building. But as the un-isolated column is itself elastic, its interaction with the pile cap is frequency dependent, and so cannot be simply modelled as a rigid mass on the pile cap.
8.17.6 Theoretical Transmissibilities
Column models were subjected to harmonic support motion to obtain their frequency response. Chapter 5 showed how such results are very sensitive to the nature and magnitude of damping assumed. In this case a constant damping ratio of 5% critical was assumed. That is, the solution at each frequency is repeated using a different set of constant coefficients to obtain the same 5% damping ratio through the frequency range.
Theoretical transmissibilities from a rigid pile cap to the top of an un-isolated column, and to the base and top of an isolated column are shown in Figure 8.51. With trace (a) comparing otherwise equivalent columns on gridline A and H and trace (b) comparing otherwise equivalent columns on gridline B and G (see Figure 8.1). This shows that the top of an isolated column does theoretically respond less than the equivalent but un-isolated column, except for the situation being either comparable or worse at certain frequencies. We also see that there are differences in the relative behaviour of isolated and un-isolated columns according to section size and column loading, compare trace(a) with trace(b). Note the small difference in behaviour between un-isolated and isolated columns at low frequencies (trace (b)), which was borne out by measurements (c.f. Figure 8.31).
We should note that these theoretical transmissibilities do not provide the actual insertion loss for base isolation, but aid in understanding the relative response of the structures compared to an input on the pile cap.
The measured mean (from 5 train events) of direct transmissibilities is compared to theoretical prediction for both un-isolated and isolated columns in Figures 8.52 and 8.53. Comparisons of direct as opposed to total transmissibilities are made with the theoretical model, because theory only considered the behaviour of one column. The measured direct transmissibility describes the response that is related to the input that arises on the pile cap, and should in theory exclude any contribution that arise from other paths (if un-correlated), or due to noise.
The measured mean direct transmissibility from pile cap to the top of the un-isolated column G15 shows correspondence with theory at the first mode resonance is excellent (Figure 8.52), which clearly confirms that a lumped mass column model is an adequate representation of the fixed-free mode, and the assumption of 5% critical damping ratio is appropriate. But there is otherwise very poor correspondence between the two. The theoretical prediction falls away at higher frequencies because the top of the column becomes a node to other axial and flexural modes of the column.
The measured mean direct transmissibility from pile cap to base and top of isolated column B15 are also compared with a theoretical prediction in Figure 8.53 which also assumes a constant damping ratio of 5% critical. Were we to have used a SDOF model, we would expect to see reduction begin to arise above Ö 2 times the rigid body natural frequency (see Figure 5.3 of chapter 5) equivalent to 10.6Hz for this comparison, but this is not the case here. However the MDOF model used shows correspondence up to 20Hz is fair, although the similarity is lost at higher frequencies.
The loss in similarity at higher frequencies in both cases can be partly explained as due to a lumped mass column model not being appropriate, since the interaction of floor mass is likely to be frequency (wavelength) dependent.
It is significant to note that we theoretically predict deep anti-resonances and troughs, yet these features are not borne out by the transmissibility measurements. The reason is that our theoretical model considers a lumped mass model of a single column. Whereas the actual structure is 3-dimensional, where many more intermediate modes exist, which arise from both variations column to column (e.g. Table 8.4), the presence of flexural modes of the frame (due to network of floors and beams that link the columns), and the influence of the piled foundations. Another reason for not seeing strong evidence of isolation is if the strength of train signals are low, leaving background response to dominate, which is particularly the case at high frequencies.
It is suggested that varying correlation of inputs (section 8.10), the transfer of vibration from neighbouring columns and the range of frame/foundation modes have a net effect that deep troughs and anti-resonances do not materialise. Improved correspondence can therefore only be realised from extended models of a type proposed by Cryer (1994).
8.18 Conclusions
This case study examined the relative behaviour of isolated (BTR @ 7.5Hz & 9.5Hz) and un-isolated buildings on a site, which are otherwise similar or identical.
Short circuits were found in some columns as grout leaks that bridged part of the pile cap to isolated column. Comparative measurements with neighbouring columns 'free' of short circuits did not show a dramatic difference, although there was evidence that free columns performed slightly better. It was suggested that the grout might have cracked by differential movements between the isolated structure and pile cap.
Transmissibilities from pile cap to base and top of isolated structure only began to show reduction with respect to pile cap levels above 20Hz, and therefore does not address a significant fraction of the spectrum associated with railway sources at grade. Furthermore a SDOF approach to modelling, where reduction is expected to arise above Ö 2 times rigid body natural frequency (º 10.6Hz in this case) is misleading.
Transmissibilities from pile cap to top of isolated and un-isolated columns that are otherwise identical were similar. Column models with lump masses to represent floor loading reproduced some of the measured characteristics at low frequencies, and are an improvement over SDOF or uniform bar models. Such models predict deep anti resonances and troughs not borne out by measurements. This could in part be due to weak train signal, where background response can dominate. However spatially extended models would be required to consider the range of modes that actually exist.
Comparisons relative to pile cap do not provide insertion loss for base isolation, obtained as difference in response between otherwise identical isolated and un-isolated structures. This showed in spectral terms that on the whole an isolated column responded more than an adjoining un-isolated core area. But we cannot generalise, as neighbouring structures were not absolutely identical and more significantly whilst at the same distance from the railway, the isolated column would receive vibration from neighbouring isolated columns nearer to the tracks, whereas other parts of the un-isolated structure were altogether further from the tracks.
Impact and shaker testing showed inertance of an isolated column to be greater than an equivalent but un-isolated column. This implies an isolated structure would be easier to excite, which was borne out by comparisons of building response to services noise.
Impact testing was used to estimate modal damping, which was found to be broadly similar in the isolated and un-isolated columns, and was of a hysteretic form.
Shaker testing was used to determine point inertance of pile caps which was larger for a pile cap below an isolated column compared to an un-isolated column. Measured inertance of pile caps below an isolated column correlated best with theory when the pile was assumed to be unloaded rather than loaded with a lump mass representing the dead load of the column above. In contrast, the measured inertance of the pile cap below the un-isolated column lay between the theoretical prediction under the two loading cases. Isolators can therefore decouple building inertia from interacting with the pile, and makes this pile easier to excite compared with a pile under an un-isolated building.
For low frequency vibration from railways at grade, base isolation of a large framed building using BTR was in this case inappropriate. However there have been no adverse comments, but clearly this in itself should not infer base isolation performance.