Seismic performance of partially submerged R.C. caissons used in port structures

Proceedings of the 9th International Conference on Structural Dynamics, EURODYN 214 Porto, Portugal, 3 June - 2 July 214 A. Cunha, E. Caetano, P. Ribeiro, G. Müller (eds.) ISSN: 2311-92; ISBN: 978-972-752-165-4 Seismic performance of partially submerged R.C. caissons used in port structures Andrea Belleri 1 and Paolo Riva 1 1 Department of Engineering, University of Bergamo, Viale Marconi 5, 2444 Dalmine, Italy email: andrea.belleri@unibg.it, paolo.riva@unibg.it ABSTRACT: In this paper the transversal stability under seismic excitation of a possible structural solution for a pier container terminal is investigated. The considered solution adopts R.C. caisson structures resting on shallow foundations as vertical elements for the pier. The main aspect investigated is the caissons attitude to overturn due to excessive rocking at the base during earthquakes being the structure characterized by a shallow foundation. Being the caissons partially submerged bodies, the hydrodynamic effects associated to seismic excitation are included. The structure is verified by means of nonlinear time history analyses considering a set of spectrum compatible earthquake records. In the finite element model the caisson is schematized with mono-dimensional beam type elements. Soil-structure interaction is accounted for by placing compression-only axial springs with stiffness compatible to gravel soils. Both linearelastic and elasto-plastic springs are considered. The resulting non-linear dynamic behaviour is characterized by limited uplift of the foundation relative to the soil, this uplift is responsible of the non-linearity, mainly geometric, related to soil-structure interaction. The analyses showed limited displacements of the caisson top, with higher values in the case of elasto-plastic soil with low elastic stiffness. KEY WORDS: Rocking; Seismic performance; Soil structure interaction; Port structures. 1 INTRODUCTION In this paper the transversal stability of a R.C. precast caisson structures for a pier container terminal is assessed under seismic excitation. The main aspect investigated is the rocking stability of the selected caisson, characterized by a shallow foundation, and the possible overturning during earthquakes. Concrete caissons are cellular gravity structures built in a construction plant, towed to the site taking advantage of their buoyancy capacity and then sunk to the seabed. While sinking, the caisson is balanced by pumping or sucking water. Once the target position is achieved, the cells are filled up with ballast to stabilize the caisson as needed. During the construction phase, the scaffolding system consists of traditional formworks for the base slab and sliding formworks for the vertical walls. In the transportation phase the caisson leaves the construction plant towed by a tugboat. A typical caisson arrangement is shown in Figure 1. The caisson has vertical and horizontal cross sections according to Figure 2. Conservatively, the beneficial stability effect due to the connection of one caisson alignment to the adjacent ones through the deck is not considered. The permanent loads considered, in addition to the self-weight of the R.C. caisson and the weight of the filling material, are constituted by the deck structure (prestressed precast beams plus R.C. slab) and the finishing, while for live loads a portion of the load due to containers is considered (participation coefficient 2 =.8). The paper considers solely along the quay excitation, therefore rocking movements around y axis (Figure 2). 19,91 19,91 19,91 19,91 19,91 14, Selected Caisson Figure 1. Typical pier arrangement and selected caisson. BB y +. x B B AA y A A -21.1 x Figure 2. Cross sections of the selected caisson. 485

Proceedings of the 9th International Conference on Structural Dynamics, EURODYN 214 2 SEISMIC ACTION 2.1 Structure classification and seismic zone The structure is classified with an importance class II (EN 1998-1 par. 4.2.5) and the related importance factor I is equal to 1. Only the No-collapse Requirement (NCR) limit state is considered herein. The related probability of exceedance P NCR in the reference return period is 1% in 5 years, corresponding to a return period T NCR of 475 years. The site of interest is characterized by a maximum ground acceleration, for soil class A, of.2 g and the expected ground type is C. 2.2 Selected spectrum-compatible earthquake records A set of 7 scaled spectrum-compatible earthquake records has been selected from the European Ground Motion Database (Ambraseys et al. 24). The chosen minimum and maximum allowable scale factors are respectively.5 and 1.5. The selected records are reported in Table 1. The comparison between the pseudo acceleration spectra of the scaled records and the Eurocode 8 spectra is reported in Figure 3, 4 and 5. It is worth noting that the choice of spectrum compatible records in the horizontal direction does not guarantee a perfect match of the vertical acceleration spectrum. The time history analyses were carried out considering both the horizontal and vertical acceleration components acting simultaneously, therefore it is not convenient to choose uncorrelated vertical and horizontal acceleration components in order to obtain spectrum compatibility in both directions. Table 1. Selected earthquake records. Note: M = Magnitude; R = Epicentral distance. No. Earthquake Name Date M R (km) Scale factor 1 Biga 5/7/1983 6 95 1.398 2 Vrancea 3/5/199 7 162.67 3 Umbria - aftershock 3/1/1997 5 2.76 4 Ierissos - foreshock 14/6/1983 4 13 1.87 5 Kallithea 18/3/1993 6 41 1.488 6 Duzce 1 12/11/1999 7 174 1.81 7 Izmit - aftershock 31/8/1999 5 39.662 Sa (g).9.8.7.6.5.4.3.2.1 Horizontal Response Spectrum EC8 Mean GM Std GM Scaled GM.5 1 1.5 2 2.5 3 3.5 4 T Figure 3. Horizontal pseudo acceleration elastic spectrum. =.5. Sd (m).35.3.25.2.15.1.5 Displacement Spectrum EC8 Mean GM Std GM Scaled GM.5 1 1.5 2 2.5 3 3.5 4 T Figure 4. Horizontal displacement elastic spectrum. =.5. 1.4 1.2 1 Vertical Response Spectrum EC8 Mean GM Std GM Scaled GM Sv (g).8.6.4.2.5 1 1.5 2 2.5 3 3.5 4 T Figure 5. Vertical pseudo acceleration elastic spectrum. =.5. 486

Proceedings of the 9th International Conference on Structural Dynamics, EURODYN 214 3 HYDRODYNAMIC EFFECTS Being the caissons partially submerged bodies, they will be subjected to hydrodynamic effects (Newman 1977). The present report considers the hydrostatic pressure and the hydrodynamics effects associated to seismic excitation. The influence of wave loads, currents and tides are not considered. In the case of ground movement, earthquake type loading, it is possible to assume that the viscous effect and the velocity gradient of the fluid are limited to the base of the structure, therefore in this report the water is considered at rest above the layer in contact to the ground: displacement, velocity and acceleration of the caissons with respect to the surrounding water are obtained as the sum of the ground component and the caisson to ground relative values. The hydrodynamic effects considered herein are the hydrostatic buoyancy and the hydrodynamic added mass. The hydrostatic buoyancy needs to be considered because it will limit shear friction capacity at the foundation level. This load is not converted into seismic mass. The hydrodynamic added mass is related to inertial effects due to pressure gradient generation by the acceleration of the fluid as a consequence of structural movements. This effect is considered adding a hydrodynamic mass related to the crosssection shape and to the distance from the free surface (Figure 6). For structures with rectangular cross-section, the added mass m could be evaluated (EN 1998-2:25 Annex F) as: where: is the fluid density 2 m C M ay (1) a y is the semi-dimension of the cross-section in the direction perpendicular to the flow Table 2. Selected earthquake records. m / a m Section z z/h (m) (kg) 21 1 AA 2.95.37 AA 2773 18.86.64 AA 54977 16.76.77 AA 66359 14.67.84 AA 72158 12.57.88 AA 75594 1.48.91 AA 7817 8.38.93 AA 79888 6.29.77 BB 43717 4.19.78 BB 517435 2.1.79 BB 52449.8 BB 26454 4 FINITE ELEMENT MODEL In the finite element model, the caisson is schematized with mono-dimensional beam type elements. In order to properly evaluate the inertial properties of the caissons the rotational inertia of the filling material is considered, hence the model is constituted by two sets of overlapping elements, the first representing the structural properties of the R.C. caisson, the second representing the ballast material, whose stiffness is neglected by reducing the elastic modulus (E). The hydrodynamic added masses considered in the model act only in the horizontal plane and therefore are assigned as translational lumped masses on the caisson centre-line at the appropriate elevation (Figure 7). m a C M is a factor depending on the ratio a y /a x Figure 6. Geometry for added mass calculation For the considered caisson, the relevant values are 2 m 321m and 2 m 13.9m for cross-section type AA and BB. A correction factor f(z) is necessary in order to consider the reduction of hydrodynamic added mass related to the position of the considered cross-section compared to the water free surface (Goyal and Chopra 1989).The added mass distribution of the selected caissons as a function of elevation is reported in Table 2. The hydrodynamic inertial contribution to vertical translation is negligible. Figure 7 Finite element model scheme with lumped masses: hydrodynamic masses in the elevation and masses related to permanent and accidental loads at the top. 487

Proceedings of the 9th International Conference on Structural Dynamics, EURODYN 214 The caisson foundation is modelled with beam elements with cross-section 19.91x.8 m. The slab portion inside the caisson cross-section is considered rigid in order to take into account the caisson stiffness. Hence the elastic modulus of the corresponding beam elements is increased by 4 orders of magnitude. As mentioned before, the interaction between adjacent caissons and the flexural stiffness and capacity of the cast in place R.C. top slab are not considered. The soil is modelled with a series of linear elastic compression-only axial springs spaced 25cm centre to centre (Figure 8). The spring stiffness is assessed by considering typical gravel soils values for the Poisson coefficient (.15.35) and for the shear modulus G (25 75 MPa). The soil unit stiffness is taken as (FEMA 356 par. 4.4.2.1.4): 1.3 G k B(1 ) where B is the foundation cross section minimum dimension (14 m). Considering 5% reduction of the shear modulus due to the seismic actions, the stiffness of each axial spring ranges between 4' kn/m ( Soft ) and 5 kn/m ( Rigid ). These two extreme values are adopted in the analyses. For each value two soil conditions are considered: linear elastic soil and elasto-plastic soil with a single spring yielding of 5 kn. (2) Regarding non-linear time history analyses, the results are reported in terms of the mean of the maximum values obtained in each record set: the horizontal and vertical components of the selected records are applied simultaneously at the caisson foundation. The relevant data are the maximum base shear and overturning moment, the maximum displacement and rotation at the top of the caisson, the maximum rotation required at the edges of the cast in place R.C. slab, the foundation maximum shear and bending moment at the intersection with the caisson base. In order to highlight the dynamic non-linear behaviour of the caisson structure, Figures 1, 11, 12, 13 and 14 show in graphical terms the results of the caisson subjected to the earthquake record 1, considering live loads and soil axial springs with stiffness 4 kn/m in the case of linear elastic and elasto-plastic (F y = 5 kn) conditions. Shear (kn) 5 4 3 2 1-1 -2-3 -4-5 -.15 -.1 -.5..5.1.15 Tip displacement (m) Figure 1 Base shear tip displacement. Figure 8 Close up of the base of the finite element model: compression-only springs at foundation level. 8.E-3 5 ANALYSES RESULTS Figure 9 shows the non-linear behaviour of the caisson, with and without live loads (q 1 ), when a horizontal load at its tip is applied (pushover). 25 2 Tip rotation (rad) 6.E-3 4.E-3 2.E-3.E+ -2.E-3-4.E-3 Shear (kn) 15 1 5 Rigid - NoYield Rigid - Yield Soft - NoYield Soft - Yield -6.E-3 2 4 6 8 1 Figure 11 Rotation at the caisson tip..5.1.15.2.25.3.35 Tip displacement (m) Figure 9 Pushover considering live loads. 488

Proceedings of the 9th International Conference on Structural Dynamics, EURODYN 214 R.C. floor rotation (rad) 1.5E-2 1.E-2 5.E-3.E+ -5.E-3-1.E-2-1.5E-2 2 4 6 8 1 Figure 12 Imposed rotation of the cast in place R.C. slab. Table 3 Caisson with live loads. Foundation stiffness K f = 4 kn/m. soil soil x tip (mm) 8.5 ±38.9 81. ±15.1 tip (rad).33 ±.16.33 ±.6 slab (rad).62 ±.3.62 ±.11 base (rad).3 ±.16.31.6 V base (kn) 35 118 ±5 385 25 479 ±1 28 M base (knm) 645 961 ±113 718 451 29 ±15 44 V fond (kn) 74 318 ±18 917 3 246 ±1 99 M fond (knm) 66 317 ±18 418 24 841 ±896 1.5 x tip Residual (mm) // // 7.8 ±8.29 Horizontal tip acceleration (g) 1..5. -.5 Table 4 Caisson with live loads. Foundation stiffness K f = 5 kn/m. soil soil x tip (mm) 6.3 ±37.5 44.2 ±11.9 tip (rad).25 ±.15.18 ±.5-1. slab (rad).47 ±.29.34 ±.9-1.5 2 4 6 8 1 Figure 13 Horizontal acceleration at the caisson tip. base (rad).22 ±.15.16 ±.5 V base (kn) 41 395 ±1 973 26 89 ±1 575 M base (knm) 751 519 ±214 465 498 757 ±3 842 V fond (kn) 92 39 ±2 494 44 46 ±6 3 1.5 M fond (knm) 86 58 ±31 377 36 74 ±5 43 Vertical tip acceleration (g) 1..5. -.5-1. -1.5 2 4 6 8 1 Figure 14 Vertical acceleration at the caisson tip. Table 3 and 4 show the non-linear time history analyses results subdivided by foundation stiffness. The results correspond to the mean value of the maximum values obtained in each record (BS EN 1998-1:24); the error corresponds to the standard deviation of the maximum values. x tip Residual (mm) // // 1.91 ±1.43 6 CONCLUSION This paper investigated the dynamic behaviour under seismic action, No-Collapse-Requirement limit state, for a caisson as a possible solution of pier container terminals. The caisson was modelled with mono-dimensional beam elements and the interaction with adjacent caissons was not considered. Soilstructure interaction was accounted by placing compressiononly axial springs with stiffness compatible to gravel soils. Both linear-elastic and elasto-plastic springs were considered. Due to the low rocking motion energy radiation in stabilized soils, no added damping associated to the soil structure interaction was included. Hydrodynamic effects were considered in terms of hydrostatic loads and hydrodynamic added masses. The inertia associated to live loads was lumped at the top of the caisson in terms of horizontal and vertical translational masses. 489

Proceedings of the 9th International Conference on Structural Dynamics, EURODYN 214 The non-linear dynamic analyses showed limited displacements of the caisson top, with higher values in the case of elasto-plastic soil with low elastic stiffness. The higher demands were recorded for finite element models including live loads and with linear-elastic soil with high stiffness. The transverse rocking motion was characterized by selfcentring with negligible permanent top displacements in the case of elasto-plastic soils. The non-linear dynamic behaviour is characterized by limited uplift of the foundation relative to the soil, this uplift is responsible of the non-linearity, mainly geometric nonlinearity, related to soil-structure interaction. The analyses showed high frequency acceleration spikes associated to the closing of the base gap in subsequent cycles. It is worth noting that the amount and magnitude of these acceleration spikes is influenced by the finite element modelling and is emphasized by adopting elastic structural elements, linear-elastic soil base-springs with high stiffness and by lumping the seismic mass associated to live loads at the caisson s top. The amount and magnitude of high frequency acceleration spikes will decrease by modelling the non-linearity of structural elements and of the soil springs (especially radiation damping) and by distributing live loads seismic masses. On the basis of the analyses results it is possible to confirm the stability of the caissons in the transverse direction and the suitability of this solution for seismic applications. REFERENCES [1] Ambraseys N, Smit P, Douglas J, Margaris B, Sigbjornsson R, Olafsson S, Suhadolc P and Costa G, Internet-Site for European Strong-Motion Data, 24, Bollettino di Geofisica Teorica ed Applicata 45(3): 113-129. [2] BS EN 1998-1:24 - Eurocode 8 - Design of Structures for Earthquake Resistance. Part 1: General rules, seismic actions and rules for buildings. [3] BS EN 1998-2:25 - Eurocode 8 - Design of Structures for Earthquake Resistance. Part 2: Bridges. [4] FEMA 356 Prestandard and Commentary for the Seismic Rehabilitation of Buildings, 2, Federal Emergency Management Agency Washington, D.C. [5] Goyal A., Chopra A. K., Simplified Evaluation of Added Hydrodynamic Mass for Intake Towers, J. Eng. Mech. 1989.115:1393-1412. [6] Newman J.N., Marine Hydrodynamics, 1977, The MIT Press, Cambridge, MA - USA. 49