Ice-Structure Interaction with Artificially Generated Force Records

A. S. J. Swamidas; D. V. Reddy; G. Purcell

doi:10.3189/S0022143000029348

Ice-Structure Interaction with Artificially Generated Force Records

Published online by Cambridge University Press: 30 January 2017

A. S. J. Swamidas,

D. V. Reddy and

G. Purcell

Show author details

A. S. J. Swamidas: Affiliation:
Faculty of Engineering and Applied Science, Memorial University of Newfoundland, St. John’s, Newfoundland, Canada
D. V. Reddy: Affiliation:
Faculty of Engineering and Applied Science, Memorial University of Newfoundland, St. John’s, Newfoundland, Canada
G. Purcell: Affiliation:
Faculty of Engineering and Applied Science, Memorial University of Newfoundland, St. John’s, Newfoundland, Canada

Article contents

Abstract
Introduction and Review of Literature
The Problem
Analysis
Results and Discussion
References

Rights & Permissions

Abstract

The paper describes the dynamic ice-structure interaction study of monopod platforms (fixed type) which are of considerable interest to the off-shore oil industry in cold regions. In view of the scantiness of continuous ice-force records, artificial-force records are generated from available small-scale field ice-force records using stochastic concepts. The method, based on the similarity between the fluctuating parts of randomly varying ice-force records and the seismic records uses a non-stationary random process obtained from filtering a shot noise through a second-order filter. The dynamic response analysis of a monopod structure located at Cook Inlet, Alaska, is carried out by hybridizing the relevant sub-routines of the computer programmes, SAP IV (Structural Analysis Programme) and EATSW (Earthquake response of Axisymmetric Tower Structures surrounded by Water). The work includes a study of the influence of soil properties on the frequencies and responses, and a comparison between the responses of fixed and elastically supported (by the soil foundation) monopod platforms.

Résumé

Enregistrements de l’interaction entre la structure de la glace el des farces produites artificiellement. Cet article décrit I’étude dynamique de l’action sur la structure de la glace de plates-formes monopodes (d’un type fixé) qui sont d’un intérêt considérable dans l’industric pétrolièrc des régions froides. En ratson du petit nombre d’enregistrements continus des efforts dans la glace, on engendre artificiellement des efforts distribués de manière aléatoire à partir des enregistrements disponibles des efforts constatés sur le terrain à petite échelle. La méthode basée sur la similitude entre les parties fluctuantes des enregistrements d’efforts dans la glace variant au hasard, et les enregistrements sismiques utilise une succession aléatoire non stationnaire obtenue en filtrant un bruit de détonation à travers un filtre du second ordre. L’analyse de la reponse dynamique d’une structure monopode, située à Cook Inlet, Alaska, est conduite en combinant les boucles utilisables des programmes d’ordinateurs SAP IV (Programme d’Analyse Structurale) et EATSW (Réponse aux Tremblements de Terre de Tours en Béton Asymétriques Entourées d’Eau), Le travail inclut une étude de l’influence des propriétés du sol sur les fréquences et les réponses et une comparaison entre les réponses des plates-formes monopodes fixes et à support élastique (par la fondation sur le sol).

Zusammenfassung

Wechselwirkungen zwischen Eis und Bauwerken mit simulierten Kräften. Die Arbeit handelt von der Untersuchung der dynamischen Wechselwirkung zwischen Eis und Bauwerk für einstützige Plattformen (ortsfester Typ), die von beträchtlichem Interesse für die küstennahe Ölsuche in kalten Gebieten ist. Wegen des Mangels an kontinuierlichen Aufzeichnungen der Eiskräfte werden künstliche Kraftaufzeich-nungen aus verfügbaren Daten von Eiskraftfeldern kleinen Massstabs unter Verwendung stochastischer Konzepte erzeugt. Die Methode beruht auf der Ähnlichkeit zwischen den fluktuierenden Anteilen von Daten zufällig veränderlicher Eiskräfte mit seismischen Aufzeichnungen und benutzt einen nichtstationären Zufallsprozess, der aus der Filterung eines Schussgeräusches durch ein Filter zweiter Ordnung hergeleitet wird. Die Analyse der dynamischen Reaktion eines einstützigen Bauwerks am Cook Inlet, Alaska, ergibt sich durch Koppelung der einschlägigen Prozeduren in den Rechnerprogrammen SAP IV (Struktur-Analyse-Programm) und EATSW (Erdbebenwirkung auf achssymmetrische, von Wasser umgebene Turmbauwerke). Die Untersuchung erstreckt sich auch auf den Einfluss von Bodeneigenschaften auf die Frequenzen und Reaktionen, Weiter werden die Reaktionen von ortsfesten und elastisch am Fundament befestigten einstützigen Plattformen miteinander verglichen.

Information

Type: Research Article
Information: Journal of Glaciology , Volume 19 , Issue 81 , 1977, pp. 265 - 283
DOI: https://doi.org/10.3189/S0022143000029348 [Opens in a new window]
Copyright: Copyright © International Glaciological Society 1977

Introduction and Review of Literature

The development of off-shore exploration for oil and gas in the Arctic and Subarctic regions has motivated the interaction studies of stationary and mobile ice-floes and ice fields with off-shore structures. Various structural forms are used for off-shore towers depending on the location—framed and trussed towers, monopods, and tripods (fixed and gravity-type structures). In this study, the transient interaction of a moving ice field with a stationary monopod structure located at Cook Inlet, Alaska, is considered. Since ice-force records at the site were not available, artificial ice-force records have been generated based on the actual field measurements of Reference BlenkarnBlenkarn (1970). The method used for generating these artificial ice-force records is also described.

Reference Blumberg and StraderBlumberg and Stradcr (1969) presented the first analysis of off-shore structures subjected to ice forces for the monopod structure at Cook Inlet, Alaska, by replacing the effect of the moving ice-floes by a static load at the highest tide level. Reference MatlockMatlock and others (1971) represented the interaction of a moving ice sheet with the structure by a physical model of vertical, elastic-plastic cantilevers, mounted on a smoothly moving rigid base, impacting a spring-mass (single-degree-of-freedom) system. The dynamic response, determined for a “saw-tooth” type of deterministic loading simulated by the model, showed good agreement with earlier field measurements of Reference PeytonPeyton (1968[a]). Sundararajan and Reference Reddy and CheemaReddy ([1974]) used spectral techniques to study the stochastic response of a single-degree-of-freedom (SDOF) model subjected to ice-loading obtained by field measurements (Blenkarn, 1970) at Cook Inlet, Alaska. They noted that the system responded to excitation for ω ≤ ω₀ (the fundamental frequency of vibration of the structure), but filtered out the excitations for ω » ω₀. This agrees with the independent conclusions arrived by Reference KivisildKivisild (1969), and Reference PeytonPeyton (1966, Reference Peyton1968[a], Reference Peyton[b]) that the force oscillations observed were due to the velocities and vibrations of the ice-floes (whose natural frequencies were usually less than one hertz) rather than due to the structure. Reference KivisildKivisild (1969) estimated that the fundamental frequencies of fluid-supported ice sheets, deforming laterally in dish form, were less than one hertz. Reference RechReeh (1972) derived relationships to determine the natural frequencies and mode shapes of ice sheets vibrating in water and observed that the frequencies of rigid-body motion of floating ice sheets were of the order of one hertz. The power-spectral-density technique for ice-force loading was extended to multi-degree-of-freedom (MDOF) systems by Reddy and Cheema (1974). Reference Reddy and Frankenstein(Reddy and others 1975, Reference Reddy[^e1976]) used response-spectrum techniques, widely used in earthquake analysis, to calculate the probable maximum dynamic response of a three-dimensional off-shore structure (idealized as a two-dimensional frame) subjected to impact loads from moving ice sheets. This was the first time response spectra had been developed for ice forces. The added water mass was assumed equal to the mass of water displaced (Reference KingKing, 1974). The ice-force records were digitized using the Nyquist criterion and the stationarity of the chosen ice-force records was verified by the Kolmogorov-Smyrnov test. The work included an extensive literature review on ice-structure interaction. Reference Reddy and FrankensteinReddy and others (1975) extended the power-spectral-density method to three-dimensional, framed, off-shore structures subjected to ice forces and determined the mean-square responses.

Monopod off-shore platforms, either fixed to the bottom or resting on it, have been suggested for exploratory and production drilling purposes in the technical brochures of Selmer (1974, 1975) and Reference Hochtief.Hochtief (1975). The Condeep (Norwegian Contractors, 1975) and Brent Field C (Offshore, 1975) concrete gravity platforms arc similar to monopod structures in their response to ice forces. Imperial Oil Limited examined different types of monopod platforms for off-shore drilling in the Canadian Beaufort Sea; Reference Jazrawi and DaviesJazrawi and Davies (unpublished) presented a detailed procedure for the analysis and design of such structures. Reference Durning and EngleDurning and Engle (1976) described the continuous monitoring techniques used for estimating the response of the Cook Inlet monopod structure to ice forces during three successive winter seasons. Reference UttUtt and others (1976) analysed the recorded dynamic response of the Cook Inlet monopod platform, and simulated the response with computer models. The purpose of the study was to determine the reason for observed changes in the fundamental frequency from 0.91 Hz in 1967 to 0.77 Hz in 1972. It was found that a change in soil properties was the main reason for alteration in the frequencies. This indicates the need for modifying Reference VandiverVandiver’s (1975) recommendation for evaluating the superstructure integrity by frequency measurements.

The dynamic response of axisymmetric towers partly submerged in water and subjected to earthquake ground motion, has been analysed by Reference Liaw and ChopraLiaw and Chopra (1973, Reference Liaw and Chopra1974). The analysis was based on the “substructure concept” for structure-water interaction described by Chakrabarti and Reference Chakrabarti and ChopraChopra (1972, Reference Chakrabarti and Chopra1973) in which the terms associated with hydrodynamic interaction in the modal equations of motion of the structure were determined as solutions of the boundary-value problems for the fluid domain. The findings were: (a) interaction with surrounding water increases the fundamental period of vibration and decreases the modal damping ratio, (b) compressibility of water is essentially of no consequence in the hydrodynamic effects on slender towers, and (c) the traditional definition of added mass, though conceptually deficient, is simple to employ and gives negligible errors for the fundamental mode of vibration. From hydrodynamic interaction studies of the cylindrical legs of deep-water platforms, Reference Liaw and ReimerLiaw and Reimer (1975) concluded that the commonly used “added mass” approach is very reasonable for slender-legged structures subjected to low-frequency excitation, but for other types of structure more sophisticated analyses are needed. A finite-element dynamic analysis of the Cook Inlet monopod tower subjected to impact by moving ice-floes was presented by Reference Swamidas and ReddySwamidas and Reddy (1976). The water-structure interaction was taken into account by the substructure approach. The ice forces were generated artificially by using a linear stochastic model with filtered Gaussian non-stationary shot-noise excitation similar to the one proposed by Reference Ruiz. and PenzienRuiz and Penzien (1969). The “semi-analytical finite-element process” of (Reference ZienkiewiczZienkiewicz 1971, ch. 13) was used in the solution.

The Problem

Procedure

The structure analysed is the monopod tower shown in Figures 1 and 2, constructed by Brown and Root Co. for Union Oil Company of California and Marathon Oil Company, and located at Cook Inlet, Alaska. The structure is composed of large-diameter (2.14 m to 8.69 m) tubular members built from steel plates. The base is rectangular in plan and is supported by piles at its four corners, as shown in Figure 2; the base rests on silty clays and in turn supports a pyramidal tubular system. A single 8.69 m diameter vertical column, framed to the base, is braced by the four radial braces that constitute the pyramidal tubular system. This vertical tubular member, which extends through the water surface and supports the platform made up of box-section plate girders, is subjected to impact forces of tidal-driven rafted pack ice, sometimes up to 1.83 m in thickness (Blumberg and Strader, 1969) as shown in Figure 2. These ice-floes impact the monopod leg twice daily with a maximum velocity of 3.05 m S^-1 and give rise to large dynamic displacements. Owing to the non-availability of detailed dimensions, the structure is modelled as shown in Figure 3, taking into consideration equivalent static and dynamic properties (viz. static deflection and frequencies). The material is assumed to be homogeneously orthotropic to take into consideration the axial and circumferential stiffeners of the cylindrical shell structure (vertical column). The radial tubular members bracing into the vertical cylindrical column at the bottom are represented by a truncated cone-type shell gripping the main column. The inclined, platform-carrying, box-section plate girders at the top are modelled by a composite truncated cone-cylinder type shell. The shell properties are “massaged” based on the frequencies and static deflections of the actual structure given by Blumberg and Strader (1969). The deck structure at the top is replaced by concentrated loads and masses at the topmost nodal points of the concentric cylinders. The mass of the water inside the tower is lumped at the inside nodal points of the finite-element grid. To study the effect of the foundation on the dynamic response of the tower, the soil below and around the tower is also taken into consideration in the analysis. The specification of the finite depth and width of the soil is based on computing costs. The soil-structure system is idealized as an axisymmetric system and broken into finite elements as shown in Figure 3. Since continuous records of ice force versus time were not available, artificial records were generated, as explained later; these records were based on the actual field measurements of ice forces reported by Reference BlenkarnBlenkarn (1970), Reference Afanas’yevAfanas’yev’s formula (1972) and the field measurements of ice thickness in Alaska reported by Reference Bilello and BatesBilello and Bates (1971, Reference Bilello and Bates1975). The frequencies, the eigenvectors, and the added masses were obtained by using the EATSW programme (a programme for Earthquake response of Axisymmetric Tower Structures surrounded by Water) developed by Reference Liaw and ChopraLiaw and Chopra (1973). The dynamic responses were determined by incorporating the relevant subroutines of SAP IV (Reference BatheBathe and others, 1974) (a Structural Analysis Program for static and dynamic response of linear systems) to form a hybrid programme.

Fig. 1.

Union Oil monopod platform at Cook Inlet, Alaska (photograph by a member of Brown & Root Company, reproduced by permission from Marine Technology Society Journal, Vol. 8, No. 7, 1974).

Fig. 2.

Cook Inlet monopod platform with ice loads (from Blumberg and Strader, 1969).

Fig. 3.

Details of sectional elevation and finite-element idealization.

Analysis

Equations of motion

The equations of motion for the tower, including hydrodynamic effects are written in terms of nodal coordinates as

(1)

where {u} is the nodal displacement matrix, [M], [C], [K] the mass (lumped structural mass), damping, and stiffness matrices respectively, {P (t)} the ice forces at various nodes, and {R (t)} the nodal loads associated with hydrodynamic pressures.

The axisymmetric tower, modelled as shown in Figure 3, is represented as an assemblage of quadrilateral toroid elements interconnected along nodal circles. The soil is broken into finite elements as shown in Figure 3(a); triangular elements have been used sparingly in order to have compatible displacement modes and to reduce the number of nodes (i.e. number of dynamic equations). The principal features of the finite-element analysis are given in the work reported by Reference Liaw and ChopraLiaw and Chopra (1973).

The hydrodynamic nodal loads are related to the accelerations of the outer surface of the tower, and hence modify the inertial properties of the structure. The normal modes of the tower without the surrounding water will not uncouple Equation (1), but can still be used to reduce the number of equations. The displacements of the tower including the effects of surrounding water are expressed in terms of the first F mode shapes as

(2)

in which Y_j(t) is the generalized displacement, and {φ_j} the shape of the jth vibration mode of the tower (with no surrounding water), φ_js are obtained as solutions of the eigenvalue problem,

(3)

where ω_j is the jth vibration frequency (with no surrounding water).

Since these nodal hydrodynamic pressures act only at the outer surface of the tower, the elements in {R (t)} corresponding to the inner nodal circles are zero. The nodal hydro-dynamic pressures associated with (i) the circumferential (θ) and vertical (z) degrees of freedom of a vertical structure-water interface and (ii) all the three degrees of freedom of the nodal points above the water are zero. The jth decoupled equation of motion of Equation (1) obtained by the normal mode method (Reference Clough and PenzienClough and Penzien, 1975, p. 191-207) is

(4)

in which

where ζj is the modal damping ratio and {φ_j^f} are the modal vectors associated with the hydrodynamic interaction.

The hydrodynamic loads are related to the accelerations of the outer surface of the tower. Neglecting the compressibility of water and the effect of surface waves, the hydrodynamic pressures on the tower are expressed as

(5)

Hence, the hydrodynamic force matrix (R(t)} is written (when θ=0) as

(6)

Introducing Equation (6) into Equation (4) and rearranging, the modal equations of motion in matrix form are written as

(7)

where

(8)

The matrix [B], known as the added water-mass matrix, is due to the interaction of the surrounding water with the tower. These quantities are obtained by solving the wave equation governing the dynamics of compressible fluids, subject to appropriate boundary conditions. For cylindrical towers, they can be obtained as closed-form solutions of the boundary value problem. For more complex shapes with non-vertical outside surfaces, the finite-element method can be used for the fluid medium to solve the equations. The small amplitude, irrotational motion of a compressible fluid is governed by the wave equation

(9)

in which p(r, z, θ, t) is the hydrodynamic pressure (in excess of hydrostatic pressure), and c is the velocity of sound in water. For water the compressibility could be neglected. Since the portion of the tower in contact with water is cylindrical and circular in plan, Equation (9) could be solved explicitly. The hydrodynamic pressure terms are (Reference Liaw and ChopraLiaw and Chopra, 1973).

(10)

where

K ₀(λ_m′r₀), K ₂, (λ_m′r₀) are modified Bessel functions of the second kind, w is the specific weight of water, g the acceleration due to gravity, r ₀ the outer radius of the partly submerged cylinder, H the height of the submerged portion of the cylindrical tower, and e_j(z) the continuous analogue of {φ_jf} given in Equation (4).

Only the radial accelerations of the outside surface of the tower need to be considered in the above solutions, since the circumferential and vertical motions of the tower do not produce any hydrodynamic pressures if water is treated as an ideal fluid. Substitution of Equation (10) into Equation (4) gives the added water mass in the modal of motion. The frequencies of the tower-water system are approximately obtained as

(11)

Generation of artificial ice-force records

Since records of ice force versus time are not available for such large structures, artificial ice-force records are generated as explained below. The method makes use of the similarity between the randomly oscillating ice-force records and seismic records to generate the ice-force records.

Available field measurements of continuous ice forces by Reference PeytonPeyton (1968[a]), Reference BlenkarnBlenkarn (1970), (Reference NeillNeill[1971], [1972], unpublished), and Reference SchwarzSchwarz (1970, Reference Schwarz[1971]) indicate resemblance to random time functions. The force variations are aperiodic and highly oscillatory about a varying mean force. Once the trend is removed, the ice-force records are similar to seismic records. The randomly oscillating part of an ice force record builds up from small amplitudes to a certain intensity that satisfies the stationarity condition for some time and then decays steadily.

The simulation of ice-force records can be satisfactorily achieved with a stationary process for the maximum response of a linear, well-damped system with a single degree of freedom, However, for systems with low damping as well as those which are non-linear or have many degrees of freedom the effect of non-stationarity becomes important when the main factor to be considered is the energy dissipated in plastic deformation of ice. The method detailed herein uses a non-stationary random process obtained from filtering a shot noise through a second-order linear filter. Both the filter and the shot noise are selected so that the artificial ice-force records generated will, on average, resemble the most relevant features of the records obtained when relatively strong ice sheets impact off-shore structures. In passing it must be noted that the ice-sheet strength is a function of temperature, ice structure, brine volume porosity, and rate of loading.

The second-order linear filter is simulated by the physical model of ice shown in Figure 4. The governing differential equation is

(12)

where a_i(t) is the acceleration of the ice before impacting the structure, and the interactive acceleration at the ice-structure interface ( times unit mass is the ice force on the structure).

Fig. 4.

Transfer-Junction model for randomly-fluctuating tee loading.

This equation is referred to as a second-order linear filter due to its analogue in electrical circuits. The details of the theory associated with the artificial generation of random ice-force records are given in the report of Reference Ruiz. and PenzienRuiz and Penzien (1969). The filter parameters ω₀ and ζ, and the variance intensity function p(t) that governs the acceleration function a_i(t), are estimated using the available ice-force records of Reference BlenkarnBlenkarn (1970) shown in Figure 5. The records are divided by unit mass, the constant means removed, and the fluctuating parts normalized to unit spectral intensities (Reference HousnerHousner, 1959). The estimation of w₀, ζ and ρ(t) is done in the frequency domain and transformed back into the time domain. Also, since it is not possible to determine the distribution of the estimators, the estimated quantities are considered to be true values of the parameters. The estimates are obtained using a Hewlett-Packard Fourier Analyser.

Fig. 5.

Blenkarn’s ice-force records (approximate velocity of ice flow = 1 m/s).

The values of ω₀ and ζ that give the best fit with the available field records are ω₀ = 1.23 Hz and ζ = 0.22, as shown in Figure 6. The variance intensity function is obtained as indicated in Figure 7; the sharp irregularities in the curve due to the reduced number of samples are smoothed out as shown. The constants t ₀, p ₀ and c are assumed to be t ₀ = T—3 s, p _a = 5.1 X 10^-5 and c = 50, where T is the total duration of the ice-force record. The durations of the ice-force records are governed by the lengths and velocities of the ice-floes traversing the structure. Moreover, the expected maximum deviation of the ice force from the mean is found to be a function of temperature, previous thermal and stress history, porosity and brine volume. In this study the deviations arc assumed to be 20% to 40% of the mean.

Fig. 6.

Filter characteristics.

Fig. 7.

Variance intensity function.

The numerical procedure used for generating records is as follows: First, samples of a white noise with flat spectral densities of power are generated with random numbers (of uniform distribution) in the range (0 → 1) with time intervals Δτ. The white noise is multiplied by a shaping function ρ(t) to get the desired intensity. The time spacing Δτ gives an indication of the frequency variation in the generated ice-force records. In this study the chosen time spacings are 0.025, 0.037 5 and 0.050 s. These waveforms are passed through a second-order linear filter to obtain the interactive ice accelerations (and hence the ice-force records). Base-line corrections are made to let the interactive ice accelerations tend to zero at the end of the process.

In estimating the varying mean of the ice-force records, the actual field ice-cover measurements of Reference Bilello and BatesBilello and Bates (1971, Reference Bilello and Bates1975) are incorporated in the formula suggested by Reference Chakrabarti and ChopraAfanas’yev (1972). Afanas’yev’s formula for the maximum ice-force developed when an ice cover is cut through by the vertical column of the tower is

(13)

where m is the shape coefficient, equal to 0.90 for a semi-circle, k = [(5h/B) +1]^1/2 for 1 ≤ B/h ≤ 6, [for B|h 1,the values are to be obtained from his graph], B is the width of the resisting structure, h the thickness of the ice sheet, and σ₀ the uniaxial compressive strength of ice.

The uniaxial compressive strengths of ice are assumed to be 1 380, 1 724 and 2 069 kN/m². The ice thickness profile used in this analysis is given in Figure 8 and the artificially generated ice-force records in Figures 9 and 10.

Fig. 8.

Thickness profile of an ice-floe.

Fig. 9.

Typical artificially generated ice-force records (amplitude variation)

Fig. 10.

Typical artificially generated ice-force records (amplitude and frequency variation).

Response

The ice-force record used in the determination of the dynamic responses is given in Figure 11; a shorter length of the ice-force record is considered to reduce the computation time. As mentioned before, the relevant portions of the SAP IV response history analysis arc isolated from the main programme to constitute a separate programme. The dynamic analysis is carried out using this segmented programme. The decoupled equations of motion are integrated by using Reference WilsonWilson’s (1972)θ-method which is an unconditionally stable integration scheme.

Fig. 11.

Ice forces acting on the structure.

Results and Discussion

The tower is modelled as an axisymmetric structure and the restraint provided by the piles is simulated by restraining the conical base of the bottom of the tower against lift while allowing horizontal movement. The material of the tower is steel with a modulus of elasticity, E = 2 000 kN/m². In order to account for the circumferential and axial stiffeners present in the vertical cylindrical shell, the structure is assumed to be homogenously orthotropic. The soil properties for silty clay are assumed to be as given by (Reference RichartRichart and others 1970, p. 354). In order to determine the influence of the soil properties on the structural frequencies and responses, a different set of soil properties was assumed for a second computation. The assumed soil and structural properties are given in Table I. The deck structure is replaced by concentrated masses (loads) at the nodal points as indicated in Table II.

TABLE I.

PHYSICAL PROPERTIES OF THE TOWER-SOIL SYSTEM

TABLE II.

ADDED MODAL MASSES

The maximum static stresses due to self weight, deck loads and water pressure are presented in Table III. The influence of the elasticity of the foundation on the structural/soil stresses appears significant. Stress redistributions occur as a result of the yielding foundation; the maximum stress in the tower with a clay foundation is found to be more than twice that for a fixed-base tower.

TABLE III.

STATIC STRESSES DUE TO SELF WEIGHT, DECK LOADS AND WATER PRESSURE

[Refer Figure 3]

Figures 12 and 13, and Table IV summarize the influence of the soil properties on the frequencies of the structure. The frequencies of the tower—soil system are given in Table IV and the modal shapes in Figures 12 and 13. Soil-structure interaction decreases the frequencies of the structure; for the soil properties considered in this study, even in the fundamental frequency the change is 31.21% from that for a fixed-base tower. This finding seems to be contrary to an earlier observation made by Reference KhannaKhanna (1969) that the fundamental frequencies are not very sensitive to changes in soil properties. For higher frequencies the differences between corresponding values are much higher. This calls for rigorous modelling of the soil media beneath the structure to get better estimates of the structural response.

Fig. 12.

Modal shapes of the tower on silty clay.

Fig. 13.

Modal shapes of the tower on day.

TABLE IV.

FREQUENCIES OF THE SOIL-STRUCTURE SYSTEM

The dynamic response of the tower for both the assumed soil properties is given in Figure 14. The dynamic deflections are observed to be more than twice the static deflection of 0.079 3 m reported by Blumberg and Strader (1969). As seen from Table V, the dynamic magnification factor for the force record given in Figure 11 is found to be 1.89 for a tower with fixed base; but the factor may be found to be higher if a larger length of record with numerous peaks is processed, and if the peak force occurs farther from the initial build-up period. The ratio η, between the dynamic deflections of the actual and the average force records is less than 1 since the force record given in Figure 2 has only a single peak in the initial build-up period of forces.

Fig. 14.

Dynamic displacements in structure and soil.

TABLE V.

RATIO OF DEFLECTIONS OF THE STRUCTURE

Static deflection =0.0741 m

Table VI shows the effect of modification of the force record (Fig. 11) on the tower response. The dynamic magnification factor and the differential deflection coefficient α₂ change from 1.89 to 2.99 and 1.00 to 1.585 respectively. This is a reasonable quantitative indication of the influence of the frequency content of the forcing functions. It can also be observed from Table V that the magnification factor is dependent on the soil properties of the foundation. It is therefore hardly necessary to stress the considerable need for more studies of this nature to evaluate the adequacy of recent code specifications for ice forces on structures.

TABLE VI.

INFLUENCE OF THE FORCE-TIME RECORD ON DEFLECTION OF THE TOWER

The method formulated in the paper, for the generation of artificial ice-force records, gives realistic estimates of dynamic ice forces imposed on off-shore structures. By the proper choice of the ice-floe thickness variation shown in Figure 8 and the other parameters of the ice-force generation model, a combination of force situations contributing to the largest magnification factor can be obtained. The in-built design safety factor for the worst possible case can thus be obtained with better accuracy by using artificially generated force records.

Acknowledgements

The support of the investigation by Imperial Oil Co. Grant No. 04-2026 and N.R.C. Grant No. A8119 is gratefully acknowledged. The authors wish to express their gratitude to Dr R. T. Dempster, Dean of the Faculty of Engineering and Applied Science, and Dr A. A. Bruneau, Vice-President of Professional Schools and Community Services, Memorial University of Newfoundland, St. John’s, Newfoundland, for their keen interest and encouragement. Thanks are due to Professor D. Dunsiger for very valuable discussions on spectral analysis and the use of the Hewlett Packard 5451B Fourier Analyser.

Discussion

J. SCHWARZ: AS we know from Peyton’s test records, the highest ice pressure occurs just before the ice-floe comes to rest. In this case no further oscillations occur. How do you consider this case?

D. V. REDDY: The ice conditions in Cook Inlet, Alaska, over several seasons indicated successive impact by ice-floes with fairly constant velocities. Peyton’s records showed a clear “ratcheting” phenomenon establishing the validity of the ice model of Reference MatlockMatlock and others (1971) consisting of elastic-plastic vertical cantilevers mounted on a rigid base. A flexible structure is deformed by the impacting ice-floe and successive impacts will subject it to oscillation. Other ice environments may produce different kinds of response but these need to be studied.

References

Afanas’yev, V. P. 1972. Davleniye l’da na vertikal’nyye pregrady [Ice pressure on vertical structures]. Transportnoye Stroitel’stvo 1972, Vyp, 3, p. 47–4 .8. [English translation: Canada. National Research Council. Technical Translation 1708, 1973.]Google Scholar

Bathe, K.-J.and Wilson, E. L. 1973. Stability and accuracy analysis of direct integration methods. International Journal of Earthquake Engineering and Structural Dynamics Vol. 1, No. 3, p. 283–91 .CrossRef Google Scholar

Bathe, K.-J.and others. 1974. SAP IV: a structural analysis program for static and dynamic response of linear systems, by K.-J. Bathe, E. L. Wilson,F. E. Peterson. Berkeley, California, Earthquake Engineering Research Center, College of Engineering, University of California. (Report No. EERC 73-11. Revised.)Google Scholar

Bilello, M. A., and Bates, R. E. 1971. Ice thickness observations, North American Arctic and Subarctic, 1966-67, 1967-68. U.S. Cold Regions Research and Engineering Laboratory. Special Report 43, Pt. V.Google Scholar

Bilello, M. A.and Bates, R. E. 1975. Ice thickness observations, North American Arctic and Subarctic, 1970-71, 1971-72. U.S. Cold Regions Research and Engineering Laboratory. Special Report 43, Pt. iVII.Google Scholar

Blenkarn, K. A. 1970. Measurement and analysis of ice forces on Cook Inlet structures. Offshore Technology Conference, April 22-34 [1970]—Houston, Texas. Preprints Vol. 2, p. 365–78 .Google Scholar

Blumberg, R.and Strader, N.R. II. 1969. Dynamic analysis of offshore structures. Offshore Technology Conference, May 18-21, 1969, Houston, Texas. [Preprint] Vol. 2, p. 107–26 .Google Scholar

Chakrabarti, P.and Chopra, A. K. 1972. Earthquake response of gravity dams including reservoir interaction effects. Berkeley, California, Earthquake Engineering Research Center, College of Engineering, University of California. (Report No. EERC 72-6.)Google Scholar

Chakrabarti, P.and Chopra, A. K. 1973. Earthquake analysis of gravity dams including hydrodynamic interaction. Earthquake Engineering and Structural Dynamics Vol. 2, No. 2, p. 143–58 .CrossRef Google Scholar

Clough, R. W.and Penzien, J. 1975. Dynamics of structures. New York, etc., McGraw-Hill.Google Scholar

Durning, P. J.and Engle, D. 1976. Vibration instrumentation system measures an offshore platform’s response to dynamic loads. ASCE/EMD specialty conference. Dynamic response of structures: instrumentation, testing methods and system identification, March 30 and 31, 1976, University of California, Los Angeles p. 143—58 .Google Scholar

Hochtief., 1975. Hochtief, Aktiengesellschaft fur Hoch-und Tiefbauten vorm. Gebr. Helfmann. Hochtief systems. Essen, Hochtief.Google Scholar

Housner, G. W. 1959. Behavior of structures during earthquakes. Proceedings of the American Society of Civil Engineers. Journal of the Engineering Mechanics Division Vol. 85, No. EM4, p. 109–29 .Google Scholar

Jazrawi, W.and Davies, J. F. Unpublished. A monopod drilling system for the Canadian Beaufort Sea. [Paper presented at Ice Tech Symposium, Montreal, 1975.]Google Scholar

Khanna, J. 1969. Elastic soil-structure interaction. In Proceedings of the fourth World Conference on Earthquake Engineering, Santiago de Chile, January 13-18, 1969. Vol. 3. Santiago de Chile, Asociaci?n Chilena de Sismología e Ingenería Antisísmica, p. A-6 143–52 .)Google Scholar

King, R. 1974. Hydroelastic models and offshore design. Civil Engineering Vol. 69, October 1974, p. 43–47 .Google Scholar

Kivisild, H. R. 1969. Ice impact on marine structures. In Ice seminar. A conference sponsored by the Petroleum Society of C.I.M, Calgary, Alberta, May 6-7, 1968. [Montreal], Canadian Institute of Mining and Metallurgy, p. 77 82. (Special Volume 10.))Google Scholar

Liaw, C.-Y.and Chopra, A. K. 1973. Earthquake response of axisymmetric tower structures surrounded by water. Berkeley, California, Earthquake Engineering Research Center, College of Engineering, University of California. (Report No. EERC 73-25.)Google Scholar

Liaw, C.-Y.and Chopra, A. K. 1974. Dynamics of towers surrounded by water. Earthquake Engineering and Structural Dynamics Vol. 3, No. 1, p. 33–49 .CrossRef Google Scholar

Liaw, C.-Y.and Reimer, R. B. 1975. Hydrodynamic interaction effects on the cylindrical legs of deepwater platforms. Proceedings. Offshore Technology Conference, May 5-8, 1975, Houston, Texas Vol. 2, p. 777–86 .Google Scholar

Matlock, H.and others. 1971. Analytical model for ice-structure interaction, by H. Matlock, W. P. Dawkins and J. J. Panak. Proceedings of the American Society of Civil Engineers. Journal of the Engineering Mechanics Division Vol. 97, No. EM4, p. 1083–92 .Google Scholar

Neill, C. R. [1971.] Ice pressure on bridge piers in Alberta, Canada. In [International Association of Hydraulic Research.] I.A.H.R. symposium: ice and its action on hydraulic structures, Reykjavik, Iceland, [8]-10 September 1970. [Delft, International Association of Hydraulic Research], paper 6.1.)Google Scholar

Neill, C. R. [1972.] Force fluctuations during ice-floe impact on piers. In [International Association of Hydraulic Research.] I.A.H.R. symposium. Ice and its action on hydraulic structures, Leningrad, USSR, 26-29 September 1972. [Leningrad, International Association of Hydraulic Research], [Vol. 1], p. 44–50 .)Google Scholar

Neill, C. R. Unpublished. Selected aspects of ice forces on piers and piles. [Paper presented at Canadian Society for Civil Engineering. Atlantic Region Hydrotechnical Conference, Fredericton, New Brunswick, 1974.]Google Scholar

Norwegian contractors. 1975. Norwegian contractors: A/F Heyer-Ellefsen, F. Selmer, Thor. Furuholmen. Oslo, A/F Høyer-Ellefsen, etc.Google Scholar

Offshore. 1975. Gravity structures gain acceptance. Offshore July 1975, p. 83–98 .Google Scholar

Peyton, H. R. 1966. Sea ice strength. College, Alaska, Geophysical Institute, University of Alaska. (UAGR-182.)Google Scholar

Peyton, H. R. 1968[a], Ice and marine structures. Ocean Industry Vol. 3, No. 3, p. 40–44 ; No. 9, p. 59-65 ; No. 12, p. 51-58.Google Scholar

Peyton, H. R. 1968[b]. Sea ice forces. Canada. National Research Council. Associate Committee on Geotechnical Research. Technical Memorandum No. 92, p. 117–24 .Google Scholar

Reddy, D. V.and Cheema, P. S. 1974. Response of an offshore structure to random ice forces. Ocean ’74. IEEE international conference on engineering in the ocean environment. Record Vol. 1, p. 84–88 .Google Scholar

Reddy, D. V.and others. 1975. Stochastic response of a three-dimensional offshore tower to ice forces, [by] D. V. Reddy, P. S. Cheema, A. S. J. Swamidas, A. K. Haldar. In Frankenstein, G. E.. ed. Proceedings, third International Symposium on lee Problems, 18-21 August 1975, Hanover, New Hampshire. [Hanover, N.H.], International Association of Hydraulic Research, p. 499–514 .)Google Scholar

Reddy, D. V.and others. [^e 1976.] Ice force response spectrum modal analysis of offshore towers, [by] D. V, Reddy, P, S. Cheema and A. S. J. Swamidas. In Proceedings of the third International Conference on Port and Ocean Engineering under Arctic Conditions held in Fairbanks, Alaska, University of Alaska, 11-15 August 1975. Vol. 2. Fairbanks, Alaska, University of Alaska. Institute of Marine Science, p. 887–910.)Google Scholar

Rech, N. 1972. Impact of an ice floe against a sloping face. Technical University of Denmark. Institute of Hydrodynamics and Hydraulic Engineering. Progress Report 26, p. 23–28 .Google Scholar

Richart, F. E. jr.,and others. 1970. Vibrations of soils and foundations, [by] F. E. Richart, Jr., J. R. Hall, Jr., R. D. Woods. Englewood Cliffs, N.J., Prentice-Hall Inc. (Prentice-Hall International Series in Theoretical and Applied Mechanics.)Google Scholar

Ruiz., P.and Penzien, J. 1969. Probabilistic study of the behavior of structures during earthquakes. Berkeley, California. Earthquake Engineering Research Center, College of Engineering, University of California.Google Scholar

Schwarz, J. 1970. Treibeisdruck auf Pfähle. Mitteilungen des Franzius-Instituts, TU Hannover Ht. 34, p. 89–193 .Google Scholar

Schwarz, J. [1971.] The pressure of floating ice-fields on piles. In [International Association of Hydraulic Research.] I.A.H.R. symposium: ice and its action on hydraulic structures, Reykjavik, Iceland, [8]-10 September 1970. [Delft, International Association of Hydraulic Research], paper 6.3.)Google Scholar

Selmer 1974. Selmer tripod. Oslo, F. Selmer.Google Scholar

Selmer. 1975. Selmer tripod 300. Oslo, P. Selmer.Google Scholar

Sundararajan., C.and Reddy, D. V. [1974.] Stochastic analysis of ice-structure interaction. In Karlsson, T.. ed. Proceedings, the second International Conference on Port and Ocean Engineering under Arctic Conditions. … Reykjavik, Iceland, August 27-30, 1973. [Reykjavik, University of Iceland], p. 345–53 .)Google Scholar

Swamidas, A. S. J.and Reddy, D. V. 1976. Dynamic ice-structure interaction of offshore monopod towers by hybridisation of the computer programs EATSW and SAP-IV. First SAP Users’ Conference and Compute Workshop. June 7-11, 1976, University of Southern California, Department of Civil Engineering, Los Angeles Session IV. paper 4.Google Scholar

Utt, M. E.and others. 1976. Estimation of the foundation condition of a fixed platform by measurement of dynamic response, by M. E. Utt, P. J. Darning, F. C. Duthweiler and D. D. Engle. Proceedings. Offshore Technology Conference, May 3-6 [1976] — Houston, Texas Vol. 3, p. 65–74 .Google Scholar

Vandiver, J. K. 1975. Detection of structural failure on fixed platforms by measurement of dynamic response, Proceedings. Offshore Technology Conference, May 5–8, 1975, Houston, Texas Vol. 2, p. 343–52 .Google Scholar

Wilson, E. L. 1972. Solid SAP: a static analysis program for three dimensional solid structures. Berkeley, California, Structures and Materials Research, Dept. of Civil Engineering, Structural Engineering Laboratory. University of California. (UC SESM 71-19. Revised.)Google Scholar

Zienkiewicz, O. C. 1971. The finite element method in engineering science. London, McGraw-Hill.Google Scholar

Fig. 1.Union Oil monopod platform at Cook Inlet, Alaska (photograph by a member of Brown & Root Company, reproduced by permission from Marine Technology Society Journal, Vol. 8, No. 7, 1974).