The purpose of this study was to analyze the anisotropic behaviour of layered soil at small strain level. This thesis comprised experimental, numerical, and theoretical aspects and provided an insight to the concept of the ‘equivalent’ homogeneous cross-anisotropic material indicating that the layered soil can be replaced by an equivalent transversely isotropic material under certain conditions. In the experimental program, this study extended Vaid’s work (1971) to conduct K0-compression tests in a rigid triaxial cell with a flexible lateral boundary for the determination of the K0-value of normally consolidated soil. The error induced by the compliance of the cell-water system that always existed in the previous study was successively eliminated by attaching a compliance correction system (i.e., GDS controller) to the triaxial cell. Three stress path tests (i.e., the K0-compression test, the plane strain compression test, and the hydrostatic pressure compression test) were conducted consecutively in a rigid triaxial cell for a layered soil specimen to determine the elastic cross-anisotropic properties of the equivalent homogeneous material. The applicability of the proposed approach was demonstrated by conducting the tests on two types of soil (i.e., the homogeneous soil and the layered soil). A three-parameter constitutive equation for describing the soil’s cross-anisotropic elastic behaviour was modified in a basic FEM program. Based on the material properties determined by the experimental study as an input, numerical simulations (i.e., the numerical K0-test) using FEM were conducted to compare numerical results with the test results. In the theoretical part, this study made use of the Reuss and the Voigt approximations and proposed a simple, yet physically meaningful, approach to determine the equivalent cross-anisotropic elastic properties of a multilayered medium. To simplify the exposition, a multilayered medium with two constituent materials that were both isotropic was examined. / Thesis / Master of Applied Science (MASc)

Since the last three decades, several studies have been conducted related to improvement in bearing capacity of pavements, embankments, and shallow foundations resting on geosynthetic reinforced soil. Most of the work has been carried out on single layer soil e.g., sand or clay layer only. Very few studies are available on a double layer soil system; but no study is available on the local soil of Carbondale, Illinois. The present study investigates the physical and engineering properties of a local soil and commonly available sand and improvement in the bearing capacity of a local soil for a rectangular footing by replacing top of the local soil with sand layer and placing geogrids at different depths. Seven tests on the model footing were performed to establish the load versus settlement curves of unreinforced and reinforced soil supporting a rectangular foundation. The improvement in bearing capacity is compared with the bearing capacity of the local soil and double layer unreinforced soil system. The test results focus on the improvement in bearing capacity of local soil and double layer unreinforced soil system in non-dimensional form i.e., BCR (Bearing Capacity Ratio). The results obtained from the present study show that bearing capacity increases significantly with the increasing number of geogrid layers. The bearing capacity for double layer soil increases, by placing three inch sand layer at the top of local soil, was not significant. The bearing capacity of the local soil increased at an average of 7% with three inches sand layer. The bearing capacity for the double layer soil increases with an average of 16.67% using one geogrid layer at interface of soils (i.e., local soil and sand) with u/B equal to 0.67. The bearing capacity for the double layer soil increases with an average of 33.33% while using one geogrid in middle of sand layer having u/B equal to 0.33. The improvement in bearing capacity for double layer soil maintaining u/B equal to 0.33 and h/B equal to 0.33; for two, three and four number geogrid layer were 44.44%, 61.11%, 72.22%, respectively. The results obtained from this research work may be useful for the specific condition or similar type of soil available anywhere to improve the bearing capacity of soil for foundation and pavement design.

No description available.

Neste trabalho, é desenvolvido um programa de computador para a análise estática e tridimensional de problemas de interação solo-estrutura. O programa permite considerar várias camadas de solo, cada qual com características físicas diferentes. Sobre este solo, o qual pode conter estacas, podem ser apoiados diversos tipos de estruturas, tais como placas e até um edifício. Todos os materiais considerados são homogêneos, isotrópicos, elásticos e lineares. O solo tridimensional é modelado com o método dos elementos de contorno (MEC), empregando as soluções fundamentais de Kelvin e uma técnica alternativa na consideração do maciço não-homogêneo. Esta técnica, que é uma contribuição original deste trabalho, é baseada no relacionamento das soluções fundamentais de deslocamento dos diferentes domínios, permitindo que sejam analisados como um único sólido sem a necessidade de equações de equilíbrio e compatibilidade. Isso reduz o sistema de equações final e melhora a precisão dos resultados, conforme comprovado nos exemplos apresentados. Para reduzir o custo computacional sem prejudicar a precisão dos resultados, é utilizada uma malha de elementos de contorno infinitos (ECI) nas bordas da malha de ECs para modelar o comportamento das variáveis de campo em longas distâncias. A formulação do ECI mapeado utilizado é outra contribuição original deste trabalho, sendo baseado em um EC triangular. É demonstrado por meio de exemplos que tal formulação é eficiente para a redução de malha, contribuindo de forma significativa na redução do custo computacional. Todas as estruturas que interagem com o solo, incluindo as de fundação, são simuladas empregando o método dos elementos finitos (MEF). Cada estaca é modelada como uma linha de carga empregando um único elemento finito com 14 parâmetros nodais, o qual utiliza funções de forma do quarto grau para aproximar os deslocamentos horizontais, do terceiro grau para as forças horizontais e deslocamentos verticais, do segundo grau para as forças cisalhantes verticais e constantes para as reações da base. Este elemento é empregado em outros trabalhos, no entanto os autores utilizam as soluções fundamentais de Mindlin na consideração da presença da estaca no solo. Desta forma, a formulação desenvolvida neste trabalho com as soluções fundamentais de Kelvin pode ser considerada mais uma contribuição original. No edifício, que pode incluir um radier como estrutura de fundação, são utilizados dois tipos de EFs. Os pilares e vigas são simulados com elementos de barra, os quais possuem dois nós e seis graus de liberdade por nó. As lajes e o radier são modelados empregando elementos planos, triangulares e com três nós. Nestes EFs triangulares são superpostos efeitos de membrana e flexão, totalizando também seis graus de liberdade por nó. O acoplamento MEC/MEF é feito transformando as cargas de superfície do MEC em carregamentos nodais reativos no MEF. Além de exemplos específicos nos Capítulos teóricos, um Capítulo inteiro é dedicado a demonstrar a abrangência e precisão da formulação desenvolvida, comparando-a com resultados de outros autores. / In this work, a computer code is developed for the static analysis of three-dimensional soil-structure interaction problems. The program allows considering a layered soil, which may contain piles. This soil may support several structures, such as shells or even an entire building. All materials are considered homogeneous, isotropic, elastic and linear. The three-dimensional soil is modeled with the boundary element method (BEM), employing Kelvin fundamental solutions and an alternative multi-region technique. This technique, which is an original contribution of this work, is based on relating the displacement fundamental solution of the different domains, allowing evaluating them as an unique solid and not requiring compatibility or equilibrium equations. In such a way, the final system of equations is reduced and more accurate results are obtained, as demonstrated in the presented examples. In order to reduce the computational cost maintaining the accuracy, an infinite boundary element (IBE) mesh is employed at the BE mesh limits to model the far field behavior. The mapped IBE utilized, based on a triangular EC, is another original contribution of this work. In the presented examples it is demonstrated that this IBE formulation is efficient for mesh reduction, implying on a significant computational cost reduction. All structures that interact with the soil, including the foundations, are simulated with de finite element method (FEM). The piles are modeled using a one-dimensional 14 parameter finite element, with forth degree shape functions for horizontal displacement approximation, third degree shape functions for horizontal forces and vertical displacement, second degree shape functions for vertical share force, and constant for the base reaction. This element is employed in other works, however the authors utilize Mindlin fundamental solutions for the pile presence consideration in the soil. In such a way, the formulation developed in this work with Kelvin fundamental solutions may be considered one more original contribution. The building, which may include a radier as a foundation structure, is modeled using two types os FEs. Piles and beams are simulated using bar FEs with two nodes and six degrees of freedom per node. The radier and pavements are modeled employing plane triangular three-node FEs. In these FEs plate and membrane effects are superposed, totalizing six degrees of freedom per node. FEM/BEM coupling is made by transforming the BEM tractions in nodal reactions in the FEM. Even though specific examples are presented in the theoretical Chapters, a role Chapter is dedicated for demonstrating the formulation accuracy and coverage. In most examples, the results are compared with the ones obtained by other authors.

Natural coarse textured soils comprise a significant portion (approximately 20%) of the area to be mined at Suncor, Syncrude (aurora mine), Albian/Shell, and CNRL mines in the Alberta’s oil sands (Macyk, 2006). Although similar in soil textural classifications, the undisturbed areas support a range of ecosite types which exhibit different moisture regimes, suggesting that there are natural mechanisms controlling the plant available water sufficient for forest development. The global objective of this study was to evaluate the potential for textural variability to enhance water storage in coarse textured soil. The observations of the infiltration and drainage behaviour of natural and reclaimed coarse-texture soils in this study have demonstrated that this potential exists and can be applied in reclamation design to achieve the ranges of soil water storage needed to establish different ecosites. Field based infiltration and drainage testing, pit excavation and sampling have been completed on 14 sites (7 natural and 7 reclaimed). Bulk saturated hydraulic conductivity and field capacity were estimated for each of the 14 sites based on the field test results. The observed transient water dynamics give an indication of the effect of layering on these material properties. Laboratory analysis of water content (650 samples), particle size (650 samples), water retention (35 samples), organic carbon (100 samples) as well as calibration of field instrumentation were completed on a large number of samples (approximate values shown in brackets above) across all sites. The laboratory analysis was used to characterize textural variability (mean and standard deviation of the particle diameter) for the layered sites and estimate the soil water retention curve (SWRC) relationships for the range of soil textures encountered at the study sites. Pedotransfer functions (PTFs) were used to investigate if there were significant differences in the residual sum of squares between estimated and measured SWRCs. The measured organic carbon was used to aid in estimating permanent wilting point (WP) used in the calculation of the available water holding capacity (AWHC) of all profiles. An investigation into the calibration of the moisture capacitance probe (MCP) was undertaken as part of a comparison of the measured and simulated volumetric water content (VWC) profiles. Water storage at the cessation of drainage was related to the soil texture and textural variability as measured in the laboratory. Sites with more textural variability generally stored more water for plant use. There appeared to be a limit to what can be considered ‘useful’ textural variability. If adjacent soil layers had too extreme a contrast in texture and therefore hydraulic conductivity, unstable/preferential flow (i.e. bypassing of some of the water and nutrients from plant roots) occurred. The total porosity calculated from field samples was often higher than the maximum measured VWC in each layer which may be indicative of one or more factors that resulted in less than full saturation being attained within the targeted 1 m depth of saturation during the test. Some of these factors include: errors in sampling leading to an overestimate of total porosity; lateral flow along textural interfaces; air entrapment within the rapidly advancing wetting front; unstable/preferential flow as a result of the high contrast in hydraulic conductivity (fine over coarse) between adjacent layers (i.e. Ks Ratio >20) or where tests were conducted on slopes (i.e. funnel flow). This latter case was common at the reclaimed sites. A modelling study of one uniform (SV10) and one layered (NLFH1) natural site was conducted. The models were built by incorporating soil properties of the layers in the various soil profiles as estimated from field and/or laboratory testing. This study offers a comparison between various PTFs and their ability to capture the soil-water storage/dynamics during infiltration and drainage testing. The Arya PTF gave a better estimation of the laboratory measured SWRCs. However, when modeling the measured infiltration and drainage testing for the relatively uniform site SV10, the Arya PTF and Modified Kovacs (MK) PTF performed similarly. The Arya PTF performing slightly better for the infiltration phase and the MK PTF performing slightly better for the drainage phase. Both PTFs gave a reasonable estimation of water storage but the MK PTF gave a better estimation of the water storage with time as compared to the Arya PTF. For the highly layered site NLFH1, neither model performed well. The Arya PTF gave a substantially better estimation of the infiltration phase and gave the better estimation of the magnitude of water storage with time, the MK PTF performed marginally better for the drainage phase and gave a better estimation of the shape of the water storage with time. Generally, the study showed that the replication of the profile water storage requirements for the layered natural ecosites (‘b’ and ‘d’ ecosites) has been achieved and can be achieved by layering (or even mixing) available coarse textured reclamation materials. This study has indicated that replicating the highly uniform ecosites (‘a’ ecosites) is where the bigger challenge lies in reclamation. Reclaiming with a diversity of target ecosites is essential to achieving the pre-disturbance land capability standard that the mine operators are bound by. The temptation may exist to simply condone reclamation that has met or exceeded the pre-exisiting land capability. However, problems with ground water recharge and regional water distribution are likely to arise if large areas of lower functioning ecosites are replaced with higher functioning ecosites.

Neste trabalho, é desenvolvido um programa de computador para a análise estática e tridimensional de problemas de interação solo-estrutura. O programa permite considerar várias camadas de solo, cada qual com características físicas diferentes. Sobre este solo, o qual pode conter estacas, podem ser apoiados diversos tipos de estruturas, tais como placas e até um edifício. Todos os materiais considerados são homogêneos, isotrópicos, elásticos e lineares. O solo tridimensional é modelado com o método dos elementos de contorno (MEC), empregando as soluções fundamentais de Kelvin e uma técnica alternativa na consideração do maciço não-homogêneo. Esta técnica, que é uma contribuição original deste trabalho, é baseada no relacionamento das soluções fundamentais de deslocamento dos diferentes domínios, permitindo que sejam analisados como um único sólido sem a necessidade de equações de equilíbrio e compatibilidade. Isso reduz o sistema de equações final e melhora a precisão dos resultados, conforme comprovado nos exemplos apresentados. Para reduzir o custo computacional sem prejudicar a precisão dos resultados, é utilizada uma malha de elementos de contorno infinitos (ECI) nas bordas da malha de ECs para modelar o comportamento das variáveis de campo em longas distâncias. A formulação do ECI mapeado utilizado é outra contribuição original deste trabalho, sendo baseado em um EC triangular. É demonstrado por meio de exemplos que tal formulação é eficiente para a redução de malha, contribuindo de forma significativa na redução do custo computacional. Todas as estruturas que interagem com o solo, incluindo as de fundação, são simuladas empregando o método dos elementos finitos (MEF). Cada estaca é modelada como uma linha de carga empregando um único elemento finito com 14 parâmetros nodais, o qual utiliza funções de forma do quarto grau para aproximar os deslocamentos horizontais, do terceiro grau para as forças horizontais e deslocamentos verticais, do segundo grau para as forças cisalhantes verticais e constantes para as reações da base. Este elemento é empregado em outros trabalhos, no entanto os autores utilizam as soluções fundamentais de Mindlin na consideração da presença da estaca no solo. Desta forma, a formulação desenvolvida neste trabalho com as soluções fundamentais de Kelvin pode ser considerada mais uma contribuição original. No edifício, que pode incluir um radier como estrutura de fundação, são utilizados dois tipos de EFs. Os pilares e vigas são simulados com elementos de barra, os quais possuem dois nós e seis graus de liberdade por nó. As lajes e o radier são modelados empregando elementos planos, triangulares e com três nós. Nestes EFs triangulares são superpostos efeitos de membrana e flexão, totalizando também seis graus de liberdade por nó. O acoplamento MEC/MEF é feito transformando as cargas de superfície do MEC em carregamentos nodais reativos no MEF. Além de exemplos específicos nos Capítulos teóricos, um Capítulo inteiro é dedicado a demonstrar a abrangência e precisão da formulação desenvolvida, comparando-a com resultados de outros autores. / In this work, a computer code is developed for the static analysis of three-dimensional soil-structure interaction problems. The program allows considering a layered soil, which may contain piles. This soil may support several structures, such as shells or even an entire building. All materials are considered homogeneous, isotropic, elastic and linear. The three-dimensional soil is modeled with the boundary element method (BEM), employing Kelvin fundamental solutions and an alternative multi-region technique. This technique, which is an original contribution of this work, is based on relating the displacement fundamental solution of the different domains, allowing evaluating them as an unique solid and not requiring compatibility or equilibrium equations. In such a way, the final system of equations is reduced and more accurate results are obtained, as demonstrated in the presented examples. In order to reduce the computational cost maintaining the accuracy, an infinite boundary element (IBE) mesh is employed at the BE mesh limits to model the far field behavior. The mapped IBE utilized, based on a triangular EC, is another original contribution of this work. In the presented examples it is demonstrated that this IBE formulation is efficient for mesh reduction, implying on a significant computational cost reduction. All structures that interact with the soil, including the foundations, are simulated with de finite element method (FEM). The piles are modeled using a one-dimensional 14 parameter finite element, with forth degree shape functions for horizontal displacement approximation, third degree shape functions for horizontal forces and vertical displacement, second degree shape functions for vertical share force, and constant for the base reaction. This element is employed in other works, however the authors utilize Mindlin fundamental solutions for the pile presence consideration in the soil. In such a way, the formulation developed in this work with Kelvin fundamental solutions may be considered one more original contribution. The building, which may include a radier as a foundation structure, is modeled using two types os FEs. Piles and beams are simulated using bar FEs with two nodes and six degrees of freedom per node. The radier and pavements are modeled employing plane triangular three-node FEs. In these FEs plate and membrane effects are superposed, totalizing six degrees of freedom per node. FEM/BEM coupling is made by transforming the BEM tractions in nodal reactions in the FEM. Even though specific examples are presented in the theoretical Chapters, a role Chapter is dedicated for demonstrating the formulation accuracy and coverage. In most examples, the results are compared with the ones obtained by other authors.

Worldwide there is considerable concern over the effects of human activities on the quantity and quality of freshwater. Measurement of infiltration behaviour will be important for improving freshwater management. This study identifies that New Zealand has a sporadic history of measuring soil water movement attributes on a limited number of soil types, although the current practical demand should be large for management of irrigation, dairy farm effluent disposal, as well as municipal / domestic waste- and storm-water disposal. Previous research has demonstrated that infiltration behaviour is governed by the interplay between numerous mechanisms including hydrophobicity and preferential flow, the latter being an important mechanism of contaminant leaching for many NZ soils. Future characterisation will need to recognise the dynamic nature of these interactions, and be able to reliably characterise the key infiltration mechanisms. Since macropores are responsible for preferential flow, it is critical that infiltration studies use a representative sample of the macropore network. The aim of this project was to study the mechanisms governing the infiltration behaviour of a layered soil in large (50 x 70 cm) monolith lysimeters, where the connectivity of the macropore network remains undisturbed. Four lysimeters of the Gorge silt loam were collected, a structured soil with four distinct layers. On each lysimeter there were four separate infiltration experiments, with water applied under suctions of 0, 0.5, 1, and 1.5 kPa by a custom-built tension infiltrometer. Each lysimeter was instrumented with 30 tensiometers, located in arrays at the layer boundaries. There was also a field experiment using ponded dye infiltration to visually define preferential flowpaths. Analysis of dye patterns, temporal variability in soil matric potential (Ψm), and solute breakthrough curves all show that preferential flow is an important infiltration mechanism. Preferential flowpaths were activated when Ψm was above -1.5 kPa. During saturated infiltration, at least 97% of drainage was through the ‘mobile’ pore volume of the lysimeter (θm), estimated among the lysimeters at 5.4 – 8.7 % of the lysimeter volume. Early-time infiltration behaviour did not show the classical square-root of time behaviour, indicating sorptivity was not the governing mechanism. This was consistent across the four lysimeters, and during infiltration under different surface imposed suctions. The most likely mechanism restricting sorptivity is weak hydrophobicity, which appears to restrict infiltration for the first 5 – 10 mm of infiltration. Overall, the Gorge soil’s early-time infiltration behaviour is governed by the dynamic interaction between sorptivity, hydrophobicity, the network of air-filled pores, preferential flow and air encapsulation. Long-time infiltration behaviour was intimately linked to the temporal dynamics of Ψm, which was in turn controlled by preferential flow and soil layer interactions. Preferential flowpaths created strong inter-layer connectivity by allowing an irregular wetting front to reach lower layers within 2 – 15 mm of infiltration. Thereafter, layer interactions dominate infiltration for long-time periods, as Ψm in soil layers with different K(Ψm) relationships self-adjusts to try to maintain a constant Darcy velocity. An important finding was that Ψm rarely attained the value set by the tension infiltrometer during unsaturated infiltration. The results show that ‘true’ steady-state infiltration is unlikely to occur in layered soils. A quasi-steady state was identified once the whole column had fully wet and layer interactions had settled to where Ψm changes occurred in unison through each soil layer. Quasi-steady state was difficult to identify from just the cumulative infiltration curve, but more robustly identified as when infiltration matched drainage, and Ψm measurements showed each layer had a stable hydraulic gradient. I conclude that the in-situ hydraulic conductivity, K(Ψm), of individual soil layers can be accurately and meaningfully determined from lysimeter-scale infiltration experiments. My results show that K(Ψm) is different for each soil layer, and that differences are consistent among the four lysimeters. Under saturated flow the subsoil had the lowest conductivity, and was the restricting layer. Most interestingly this pattern reversed during unsaturated flow. As Ψm decreased below -0.5 to -1 kPa, the subsoil was markedly more conductive, and the topsoil layers became the restricting layers. All four soil layers demonstrate a sharp decline in K(Ψm) as Ψm decreases, with a break in slope at ~ -1 kPa indicating the dual-permeability nature of all layers.

The problem of pile - soil interaction is examined in the Thesis at hand by means of both theoretical analyses and experimental investigations. Pile foundations in seismically prone areas are subjected to both direct loading, such as axial and lateral forces imposed at their heads, resulting from a phenomenon known as inertial interaction, and indirect loading along their body, such as imposed displacements due to the passage of various types of seismic waves, resulting from a phenomenon known as kinematic interaction. Along this vein, a family of analytical models of the Tajimi type are presented in the framework of linear elastodynamic theory to explore the effects of axial and lateral pile - soil interaction in homogeneous and inhomogeneous soil under static and dynamic (kinematic and inertial) loading. Apart from simplified two-dimensional models of the Baranov - Novak type, few analytical solutions are available to tackle these problems in three dimensions, the majority of which are restricted to the analysis of an elastic half space under static conditions. The proposed models are based on a continuum solution pioneered by Japanese investigators (notably Matuso & Ohara and Tajimi) in the 1960’s. In the realm of this approach the soil is modelled as a continuum, while the pile is conveniently modelled as a rod or a beam by strength-of-materials theory. Displacements and stresses are expressed through Fourier series in terms of the natural modes of the soil medium. Fundamental to the analysis presented in this study is that the influence of horizontal soil displacement on axial pile response and vertical displacement on lateral response, respectively, are negligible. However, their effect on stresses is not negligible which differentiates the proposed models from the classical Tajimi solutions in which the aforementioned displacements are set equal to zero. The above approximations are attractive, as they lead to a straightforward uncoupling of the equations of motions, even in inhomogeneous media, unlike the classical elastodynamic theory where the uncoupling is generally impossible in presence of inhomogeneity. Although approximate, the proposed models are advantageous over available analytical models and rigorous numerical schemes, as they require relative simple computations and provide excellent predictions of pile response at the frequency ranges of interest in earthquake engineering and geotechnics. In addition, they are advantageous over existing simplified analytical approaches of the Winkler type, as they are more accurate, self - standing, free of empirical constants and provide more realistic simulation of the problem. The main advantage over numerical methods (finite and boundary elements) lies in the derivation of the solution in closed form and the elucidation of complex mechanisms related to the dynamic interaction phenomenon, such as radiation damping and wave propagation in in homogeneous media. The main goal of the theoretical effort lies in the derivation of solutions in closed - form for: (i) the static stiffness and the dynamic impedances (dynamic stiffness and damping coefficients) at the pile head, (ii) translational and rotational kinematic response factors (pile head displacement or rotation over free-field response), (iii) actual, depth- dependent, Winkler moduli (spring and damping coefficients), (iv) corresponding average, depth- independent, Winkler moduli to match the pile head stiffness. In addition, simple approximate formulae for Winkler moduli to be used in engineering practice are proposed, to improve the predictions of Winkler models. Pile-to-pile interaction is investigated on the basis of the superposition method for axially loaded piles. Closed-form expressions for attenuation functions are derived to be used individually or in conjunction with more elaborate methods providing more accurate predictions for static and dynamic interaction factors to assess the vertical stiffness of pile groups. New dimensionless frequency ratios controlling pile response are introduced. Finally, new solutions are added in the context of analytical Winkler models for investigating the behaviour of piles under kinematic loading due to vertically-propagating S waves. Emphasis is given on the influence of boundary conditions of the pile. With reference to kinematic pile bending, insight into the physics of the problem is gained through a rigorous superposition scheme involving an infinitely-long pile excited kinematically, and a pile of finite length excited by a concentrated force and a moment at the tip. Contrary to the classical elastodynamic theory where pile response is governed by six dimensionless ratios, in the realm of Winkler theory three only ratios suffice to fully describe the interaction problem, from which the mechanical slenderness and the effective dimensionless frequency are introduced for the first time. The selection of an appropriate value for the Winkler modulus in the accuracy of the kinematic Winkler model is demonstrated. The theoretical results are compared to new experimental data obtained from a series of tests on piles carried out on scaled models performed on the shaking table at University of Bristol Laboratory (BLADE) within the framework of the Seismic Engineering Research Infrastructures (SERIES) program, sponsored by FP7, and contribute in the investigation of pile - soil interaction. / Στην παρούσα διατριβή εξετάζεται το πρόβλημα της αλληλεπίδρασης πασσάλου - εδάφους μέσω συνδυασμένης θεωρητικής ανάλυσης και πειραματικής διερεύνησης. Οι πάσσαλοι, ως μέσο θεμελίωσης σε σεισμογόνες περιοχές, υπόκεινται σε άμεση φόρτιση στην κεφαλή, μέσω δυνάμεων και ροπών, ως αποτέλεσμα του φαινομένου της αδρανειακής αλληλεπίδρασης, αλλά και σε έμμεση φόρτιση σε όλο τους το μήκος, μέσω επιβαλλόμενων εδαφικών μετακινήσεων, ως αποτέλεσμα του φαινομένου της κινηματικής αλληλεπίδρασης. Στην κατεύθυνση αυτή παρουσιάζεται η ανάπτυξη οικογένειας αναλυτικών προσομοιωμάτων τύπου Tajimi, στο πλαίσιο της γραμμικής ελαστοδυναμικής θεωρίας, για τη διερεύνηση της αξονικής και πλευρικής αλληλεπίδρασης πασσάλου - εδάφους σε ομοιογενείς και ανομοιογενείς εδαφικούς σχηματισμούς, υπό στατική και δυναμική φόρτιση κινηματικής και αδρανειακής μορφής. Εκτός από απλοποιημένα διδιάστατα προσομοιώματα τύπου Baranov - Novak, ελάχιστες αναλυτικές λύσεις είναι διαθέσιμες σε τρεις διαστάσεις, οι περισσότερες των οποίων περιορίζονται στην ανάλυση ελαστικού ημίχωρου υπό στατικές συνθήκες. Τα προτεινόμενα προσομοιώματα βασίζονται σε μια πρωτοποριακή λύση συνεχούς μέσου (κατά Matsuo & Ohara και Tajimi) η οποία αναπτύχθηκε στη δεκαετία του 1960, αλλά δεν επεκτάθηκε ουσιαστικά μέχρι την παρούσα εργασία. Στο πλαίσιο αυτής της προσέγγισης το έδαφος προσομοιώνεται ως συνεχές μέσο και ο πάσσαλος ως ράβδος ή δοκός σύμφωνα με τη τεχνική θεωρία της κάμψης (παραδοχή επιπεδότητας των διατομών), ενώ οι μετακινήσεις και οι τάσεις εκφράζονται μέσω αναπτυγμάτων Fourier σε όρους των φυσικών ιδιομορφών του εδαφικού μέσου. Θεμελιώδης παραδοχή της προτεινόμενης ανάλυσης είναι ότι η επιρροή της οριζόντιας εδαφικής μετακίνησης στην αξονική απόκριση του πασσάλου, αλλά και η επιρροή της κατακόρυφης μετακίνησης στην πλευρική απόκριση θεωρούνται αμελητέες, ωστόσο η επίδρασή τους στις τάσεις είναι μη μηδενική, πράγμα που τις διαφοροποιεί από τις κλασσικές λύσεις τύπου Tajimi στις οποίες οι ανωτέρω μετακινήσεις μηδενίζονται. Οι ανωτέρω προσεγγίσεις κρίνονται ως ιδιαίτερα ελκυστικές καθώς οδηγούν στην άμεση απόζευξη των εξισώσεων της κίνησης, ακόμη και σε ανομοιογενή μέσα, αντίθετα με την κλασσική ελαστοδυναμική θεωρία, η απόζευξη είναι γενικώς αδύνατη παρουσία εδαφικής ανομοιογένειας. Παρά τον προσεγγιστικό τους χαρακτήρα, τα προτεινόμενα αναλυτικά προσομοιώματα πλεονεκτούν ως προς διαθέσιμα αναλυτικά προσομοιώματα και αυστηρά αριθμητικά σχήματα, καθώς απαιτούν σχετικά απλούς υπολογισμούς και παρέχουν εξαιρετικές προβλέψεις της απόκρισης του πασσάλου για το εύρος συχνοτήτων που παρουσιάζει ενδιαφέρον στη σεισμική μηχανική και τα γεωτεχνικά. Επιπρόσθετα, υπερτερούν ως προς υφιστάμενες αναλυτικές προσεγγίσεις τύπου Winkler, καθώς είναι ακριβέστερα, αυτόνομα, απαλλαγμένα από εμπειρικές σταθερές και προσφέρουν ρεαλιστικότερη προσομοίωση του προβλήματος. Το κύριο πλεονέκτημα έναντι των αριθμητικών μεθόδων (πεπερασμένα και συνοριακά στοιχεία) έγκειται στην εξαγωγή της λύσης σε κλειστή μορφή και στη διερεύνηση πολύπλοκων φαινομένων που σχετίζονται με την αλληλεπίδραση πασσάλου - εδάφους, όπως αυτό της απόσβεσης ακτινοβολίας και της διάδοσης κυμάτων στο έδαφος - ειδικά παρουσία ανομοιογένειας. Ο κύριος στόχος της θεωρητικής διερεύνησης υλοποιείται με την εξαγωγή λύσεων σε κλειστή μορφή για: (i) τη στατική και δυναμική στιφρότητα και απόσβεση στην κεφαλή του πασσάλου, (ii) τους συντελεστές κινηματικής απόκρισης σε μετάθεση και στροφή, (iii) τους πραγματικούς, συναρτήσει του βάθους, συντελεστές Winkler (συντελεστής στιφρότητας ελατηρίων και συντελεστής απόσβεσης), (iv) τους αντίστοιχους μέσους, ανεξάρτητους από το βάθος, συντελεστές Winkler. Επιπρόσθετα, παρουσιάζονται ακριβέστερες των διαθέσιμων στη βιβλιογραφία απλές προσεγγιστικές σχέσεις για τον υπολογισμό του συντελεστή Winkler με σκοπό τη βελτίωση της ακρίβειας των προσομοιωμάτων Winkler. Διερευνάται η αλληλεπίδραση πασσάλου προς πάσσαλο στην περίπτωση αξονικά φορτισμένων πασσάλων με βάση την αρχή της επαλληλίας. Εξάγονται λύσεις σε κλειστή μορφή για τις συναρτήσεις εξασθένισης ώστε να χρησιμοποιηθούν αυτόνομα ή σε συνδυασμό με πιο εκλεπτυσμένες λύσεις δίνοντας με στόχο ακριβέστερες προβλέψεις για τους συντελεστές αλληλεπίδρασης, οδηγώντας έτσι σε πιο ρεαλιστικές εκτιμήσεις της κατακόρυφης στιφρότητας ομάδας πασσάλων. Εισάγονται νέοι αδιάστατοι λόγοι συχνοτήτων που καθορίζουν την απόκριση του πασσάλου. Τέλος, παρουσιάζονται νέες λύσεις σε αναλυτικά προσομοιώματα Winkler για τη διερεύνηση της συμπεριφοράς πασσάλων υποκείμενων σε φόρτιση λόγω της κατακόρυφης διάδοσης διατμητικών κυμάτων στο έδαφος, με έμφαση στην επίδραση των οριακών συνθηκών του προβλήματος. Σε αντίθεση με την κλασσική ελαστοδυναμική θεωρία που η απόκριση του πασσάλου καθορίζεται από έξι αδιάστατους λόγους, στο πλαίσιο της θεωρίας Winkler επαρκούν μόνο τρεις για την πλήρη περιγραφή της αλληλεπίδρασης πασσάλου - εδάφους, εκ των οποίων η μηχανική λυγηρότητα και η ενεργός αδιάστατη συχνότητα παρουσιάζονται για πρώτη φορά. Καταδεικνύεται η σημασία επιλογής της κατάλληλης τιμής του συντελεστή Winkler στην ακρίβεια των εν λόγω προσομοιωμάτων. Προτείνεται σύστημα υπέρθεσης που αποτελείται από ένα απειρομήκη πάσσαλο που διεγείρεται κινηματικά και έναν πάσσαλο πεπερασμένου μήκους που υπόκειται σε αδρανειακή φόρτιση για τη διαλεύκανση της λειτουργίας του θεμελιώδους μηχανισμού που καθορίζει την κινηματική κάμψη του πασσάλου. Τα θεωρητικά αποτελέσματα συγκρίνονται με νέα πειραματικά δεδομένα από σειρά δοκιμών σε διάταξη πασσάλων υπό κλίμακα που εκτελέστηκαν στο σεισμικό προσομοιωτή του Πανεπιστήμιου του Bristol στο πλαίσιο του Ευρωπαϊκού Προγράμματος SERIES, το οποίο χρηματοδοτήθηκε από το κοινοτικό πλαίσιο FP7 που συμβάλλουν στην περαιτέρω διερεύνηση του φαινομένου της κινηματικής αλληλεπίδρασης εδάφους - πασσάλου.

Worldwide there is considerable concern over the effects of human activities on the quantity and quality of freshwater. Measurement of infiltration behaviour will be important for improving freshwater management. This study identifies that New Zealand has a sporadic history of measuring soil water movement attributes on a limited number of soil types, although the current practical demand should be large for management of irrigation, dairy farm effluent disposal, as well as municipal / domestic waste- and storm-water disposal. Previous research has demonstrated that infiltration behaviour is governed by the interplay between numerous mechanisms including hydrophobicity and preferential flow, the latter being an important mechanism of contaminant leaching for many NZ soils. Future characterisation will need to recognise the dynamic nature of these interactions, and be able to reliably characterise the key infiltration mechanisms. Since macropores are responsible for preferential flow, it is critical that infiltration studies use a representative sample of the macropore network. The aim of this project was to study the mechanisms governing the infiltration behaviour of a layered soil in large (50 x 70 cm) monolith lysimeters, where the connectivity of the macropore network remains undisturbed. Four lysimeters of the Gorge silt loam were collected, a structured soil with four distinct layers. On each lysimeter there were four separate infiltration experiments, with water applied under suctions of 0, 0.5, 1, and 1.5 kPa by a custom-built tension infiltrometer. Each lysimeter was instrumented with 30 tensiometers, located in arrays at the layer boundaries. There was also a field experiment using ponded dye infiltration to visually define preferential flowpaths. Analysis of dye patterns, temporal variability in soil matric potential (Ψm), and solute breakthrough curves all show that preferential flow is an important infiltration mechanism. Preferential flowpaths were activated when Ψm was above -1.5 kPa. During saturated infiltration, at least 97% of drainage was through the ‘mobile’ pore volume of the lysimeter (θm), estimated among the lysimeters at 5.4 – 8.7 % of the lysimeter volume. Early-time infiltration behaviour did not show the classical square-root of time behaviour, indicating sorptivity was not the governing mechanism. This was consistent across the four lysimeters, and during infiltration under different surface imposed suctions. The most likely mechanism restricting sorptivity is weak hydrophobicity, which appears to restrict infiltration for the first 5 – 10 mm of infiltration. Overall, the Gorge soil’s early-time infiltration behaviour is governed by the dynamic interaction between sorptivity, hydrophobicity, the network of air-filled pores, preferential flow and air encapsulation. Long-time infiltration behaviour was intimately linked to the temporal dynamics of Ψm, which was in turn controlled by preferential flow and soil layer interactions. Preferential flowpaths created strong inter-layer connectivity by allowing an irregular wetting front to reach lower layers within 2 – 15 mm of infiltration. Thereafter, layer interactions dominate infiltration for long-time periods, as Ψm in soil layers with different K(Ψm) relationships self-adjusts to try to maintain a constant Darcy velocity. An important finding was that Ψm rarely attained the value set by the tension infiltrometer during unsaturated infiltration. The results show that ‘true’ steady-state infiltration is unlikely to occur in layered soils. A quasi-steady state was identified once the whole column had fully wet and layer interactions had settled to where Ψm changes occurred in unison through each soil layer. Quasi-steady state was difficult to identify from just the cumulative infiltration curve, but more robustly identified as when infiltration matched drainage, and Ψm measurements showed each layer had a stable hydraulic gradient. I conclude that the in-situ hydraulic conductivity, K(Ψm), of individual soil layers can be accurately and meaningfully determined from lysimeter-scale infiltration experiments. My results show that K(Ψm) is different for each soil layer, and that differences are consistent among the four lysimeters. Under saturated flow the subsoil had the lowest conductivity, and was the restricting layer. Most interestingly this pattern reversed during unsaturated flow. As Ψm decreased below -0.5 to -1 kPa, the subsoil was markedly more conductive, and the topsoil layers became the restricting layers. All four soil layers demonstrate a sharp decline in K(Ψm) as Ψm decreases, with a break in slope at ~ -1 kPa indicating the dual-permeability nature of all layers.