[Truncated abstract] The demand for offshore wind turbines is increasing in densely populated areas, such as Europe. These constructions are typically founded on a gravity foundation or a large 'mono pile'. Gravity foundations can only be used at locations where strong soils exist and water depths are limited. Costs associated with a 'mono pile' type foundation contribute to a very large percentage of the total investment costs. This research, therefore, focuses upon a different foundation for offshore wind turbines, namely suction caissons beneath a tripod. This foundation can be used in all kinds of soil types and is cheaper than the 'mono pile' foundation, both in the amount of steel used and installation costs. Cheaper foundations can contribute to a more competitive price for offshore wind energy in comparison with other energy resources. To date, there have been relatively few studies to investigate the behaviour of this type of foundation during the installation process and during operational and ultimate loading for seabed conditions comprising dense sand. Two types of investigations were performed during this research to determine the behaviour of suction caissons beneath a tripod. Firstly, an existing computer program was extended to predict the typical loading conditions for a tripod foundation. Secondly, centrifuge tests on small scale suction caissons were performed to investigate the behaviour during the installation and loading phases. The computer program developed helped to quantify the likely ranges of environmental loading on an offshore wind turbine. For a typical 3 MW wind turbine of 90 m height, the vertical load is low at around 7 MN. During storm conditions the horizontal hydrodynamic load can be in the order of 4 MN. During normal working conditions the horizontal aerodynamic loads can reach 0.4 MN, but can increase to 1.2 MN when the pitch system malfunctions and gusts reach 30 m/s. This aerodynamic load will result in a very large contribution to the overturning moment, due to the high action point of this load. When the wind turbine is placed on top of a tripod, these large moments are counteracted by a push-pull system. ... The development of differential pressure was found to depend on the soil permeability, the extraction speed and a consolidation effect. During cyclic loading no obvious signs of a decrease in resistance were observed. During very fast cyclic loading differential pressures developed, which could increase the drained frictional resistance by approximately 40%. All centrifuge tests results were used to develop methods to predict or back calculate the installation process of suction caissons in sand and layered soil, and the behaviour during tensile and cyclic loading. These methods all use the cone resistance as the main input parameter and predict the force (or required suction) as a function of time, for a given rate of pumping or uplift displacement, in addition to the variation of suction with penetration (or force with uplift displacement). These new methods provide a useful tool in designing a reliable foundation for offshore wind turbines consisting of a tripod arrangement of suction caissons embedded in dense sand.

The research carried out in this thesis has concentrated on the application of numerical solutions to geotechnical penetration problems in offshore engineering. Several important issues closely relevant to deep-water oil and gas developments were investigated, covering installation of suction caisson foundations, interpretation of fullflow penetrometers and shallow penetration of a cylindrical object (submarine pipeline or T-bar), all in clayey sediments such as are often encountered in deep-water sites. These problems are commonly characterised by large vertical movements of structural elements relative to the seabed. A large deformation finite element method was adopted and further developed to simulate these challenging problems, referred to as Remeshing and Interpolation Technique with Small Strain. In this approach, a sequence of small strain Lagrangian increments, remeshing and interpolation of stresses and material properties are repeated until the required displacement has been reached. This technique is able to model relative motion between the penetrating objects and the soil, which is critical for evaluating soil heave inside the caissons, the effect of penetration-induced remoulding on the resistance of full-flow penetrometers, and influence of soil surface heave on the embedment of pipelines. '...' Simple expressions were presented allowing the resistance factors for the T-bar and ball penetrometers to be expressed as a function of the rate and strain-softening parameters. By considering average strength conditions during penetration and extraction of these full-flow penetrometers, an approximate expression was derived that allowed estimation of the hypothetical resistance factor with no strain-softening, and hence an initial estimate of the stain-rate dependency of the soil. Further simulations of cyclic penetration tests showed that a cyclic range of three diameters of the penetrometers was sufficient to avoid overlap of the failure mechanism at the extremes and mid-point of the cyclic range. The ball had higher resistance factors compared with the T-bar, but with similar cyclic resistance degradation curves, which could be fitted accurately by simple expressions consistent with the strain-softening soil model adopted. Based on the curve fitting, more accurate equations were proposed to deduce the resistance factor with no strain-softening, compared with that suggested previously based on the resistances measured in the first cycle of penetration and extraction. The strain-rate dependency was similar in intact or post-cyclic soil for a given rate parameter. The resistance factor for the post-cyclic condition was higher than that for the initial conditions, to some degree depending upon soil sensitivity and brittleness parameter. For the shallow penetration of a cylindrical object, the penetration resistance profile observed from centrifuge model tests was very well captured by the numerical simulation. The mechanism of shear band shedding was reproduced by the numerical technique, although the frequency of the shear band generation and the exact shape of the heave profile were not correctly captured, which were limited by the simple strainsoftening soil model adopted.

Caisson or pier foundations are encountered as part of the foundation system of tall structures such as bridges, transmission towers, heliostats, etc, and correspond to rigid blocks of length-to-diameter (D/B) ratio on the order of D/B = 2-6. As a result of their geometry and stiffness characteristics, the mechanisms of load transfer from the superstructure to the surrounding soil and their kinematic response to seismic wave propagation are governed by a complex stress distribution at the pier-soil interface, which cannot be adequately represented by means of simplified Winkler models for shallow foundations or flexible piles. Continuum model solutions, such as 3D finite elements (FE) cannot be employed frequently in practice for the design of non-critical facilities due to the cost and effort associated with these analyses. The objective of this work is to develop a Winkler-type model for the analysis of transversely-loaded caissons, which approximately accounts for all the main soil resistance mechanisms mobilized, while retaining the advantages of simplified methodologies for design at intermediate levels of target accuracy. Investigation of the governing load-transfer mechanisms and development of complex spring functions is formulated on the basis of 3D FE simulations. Initially, the soil-structure stiffness matrix is computed by subjecting the pier to transverse static and dynamic loading at the top, and numerically estimating the response. Complex frequency-dependent functions are next developed for the spring constants by equating the stiffness matrix terms to the analytical expressions developed for the four-spring model. Sensitivity analyses are conducted for optimization of the truncated numerical domain size, finite element size and far-field dynamic boundary conditions to avoid spurious wave reflections. Simulations are next conducted to evaluate the transient response of the foundation subjected to vertically propagating shear waves, and results are compared to the response predicted by means of the 4-spring model. Finally, the applicability of the method is assessed for soil profiles with depth-varying properties. While the methodology developed is applicable for linear elastic media with no material damping, the expressions of complex spring functions may be extended to include hysteretic damping, nonlinear soil behavior and soil-foundation interface separation, as shown in the conclusions.