The internal structure of ice sheets reflects the history of flow and deformation experienced by the ice mass. Flow and deformation are controlled by processes occurring within the ice mass and at its boundaries, including surface accumulation or ablation, ice rheology, basal topography, basal sliding, and basal melting or freezing. The internal structure and basal environment of ice sheets is studied with ice-penetrating radar. Recently, radar observations in Greenland and Antarctica have imaged large englacial structures rising from near the bed that deform the overlying stratigraphy into anticlines, synclines, and overturned folds. The mechanisms that may produce these structures include basal freeze-on, travelling slippery patches at the ice base, and rheological contrasts within the ice column. In this thesis, I explore the setting and mechanisms that produce large basal stratigraphic structures inside ice sheets. First, I use radar data to map subglacial hydrologic networks that deliver meltwater uphill towards freeze-on structures in East Antarctica. Next, I use a thermomechanical flowline model to demonstrate that trains of alternating slippery and sticky patches can form underneath ice sheets and travel downstream over time. The disturbances to the ice flow field produced by these travelling patches produce stratigraphic folds resembling the observations. I then examine the overturned folds produced by a single travelling sticky patch using a kinematic flowline model. This model is used to interpret stratigraphic measurements in terms of the dynamic properties of basal slip. Finally, I use a simple local one-dimensional model to estimate the thickness of basal freeze-on that can be produced based on the supply of available meltwater, the thermal boundary conditions, ice sheet geometry, and the ice flow regime.

As new materials are developed and further concerns on green alternatives and serviceability arise, understanding material behavior during the entire span of their lifetime becomes crucial to engineering applications. Moreover, many problems display a significant dependence to time and loading effects which, by varying across multiple time scales, require material models that incorporate these effects into any valid characterization and prediction. This dissertation aims at proposing a new approach to analyze and predict viscoelastic materials that deteriorate during multiple loading conditions. The model is constructed from mechanical and mathematical basis while satisfying physical laws. In this work, the proposed constitutive law is used for the analysis of the mechanical properties of ice. The mechanical behavior, biaxial envelop and multiple loading types demonstrate the validity of the model when compared to experimental results and other ice models available in the literature. A rigorous calibration scheme for the viscoelastic and damage parameters is also presented. Moreover, as material deterioration or damage is modeled in standard Finite Elements software, it is commonly known that computational results can be dependent on the spatial discretization or mesh. That is, damage zone and energy dissipation are dependent on the selection of the mesh yielding a disappearing damage zone and energy dissipation upon refinement. This non-physical behavior is corrected by the novel regularization approach proposed in this document, which introduces a length scale of the material and produces results that are no longer sensitive to the mesh selection. The nonlocal damage model is finally used in the analysis of asphalt concrete viscoelastic behavior and cracking prediction. As presented in the ice case, a rigorous calibration approach is presented first followed by the validation to experimental data available in the literature under different loading conditions. The coupled viscoelastic and damage model is compared to other model and their Finite Elements implementations are highlighted in terms of computational efficiency. A nonlinear coupled system for solving this problem is programmed as a User Element in a commercial Finite Element analysis software.

The lateral "corners" where Kamb and Whillans Ice Streams (KIS and WIS) discharge into the Ross Ice Shelf share common geometries and ice mechanical settings. At both corners of the now-stagnant KIS outlet, shear margins of apparently different ages confine regions with a relatively flat, smooth surface expression. These features are called the "Duckfoot" on the northern, right-lateral side and the "Goosefoot" on the other. It has been suggested, on evidence found in ice internal layers, that the flat ice terrains on KIS were afloat in the recent past, at a time when the ice stream grounding line was upstream of its present location. The overdeepening in the bed just upstream of the KIS grounding line supports this view of the past geometry. The right-lateral margin at the outlet of the currently active WIS, the location of Subglacial Lake Englehardt (SLE), appears to have many similarities with the right lateral margin of KIS, though with a less developed looking inboard margin. This paper presents a mechanical analysis using surface and bed topography and velocity datasets comparing the Duckfoot flat ice terrain with the terrain around Subglacial Lake Englehardt. At both locations mechanical thinning along shear margins and lows in the bed topography redirects basal water routing towards the features. Here, I consider the history of these features and their role in ice stream variability by comparison of the relict and modern features and via numerical modeling of ice shelf grounding and ungrounding in response to variations in ice flow. We propose two scenarios for the development of flat ice terrains/subglacial lakes at the outlets of ice streams. In the first, development of a lake in the hydraulic potential low along a shear margin forces a margin jump as shearing develops along the inboard shore of the margin lake. This thesis presents evidence for an inboard (relative to the main outboard shear margin) zone of shear along the inboard shoreline of SLE, suggesting that subglacial lakes along shear margins are capable of facilitating shear margin jumps. In the second, grounding line advance around a relative low in the bed, creating adjacent margins along the lakeshores, forms a remnant lake. Discerning which of these scenarios is appropriate at the KIS outlet has implications for understanding the history of the ice stream grounding line. An ice flow model is used to place these local conditions in a regional context by studying the effect of internal perturbations, such as ice rise stagnation or inward margin jumps, on grounding line position. Bathymetry is important in determining ice stream flow in the ways that might not be otherwise realized in 1-D flow model studies. In the numerical modeling experiments, grounding line advance across the KIS outlet is mediated by the overdeepening in the bed and proceeds not in the direction of ice flow but transverse to flow. This finding adds complexity to both a flowline view of grounding line migration and the theory that grounding lines are unstable in the presence of inward sloping bed topography.

No description available.

Thesis (M.Eng.)--Memorial University of Newfoundland, 2003. / Bibliography: leaves 119-127.

The evaluation of design iceberg impact loads for offshore structures and the influence of hydrodynamic effects on impact loads are examined. Important hydrodynamic effects include iceberg added mass, wave-induced oscillatory iceberg motions, and the influence of a large structure on the surrounding flow field and subsequent velocities of approaching icebergs. The significance of these phenomena has been investigated using a two-body numerical diffraction model and through a series of experiments modelling the drift of various sized icebergs driven by waves and currents approaching a large offshore structure. Relevant findings from the hydrodynamic studies have been incorporated into two probabilistic models which can be used to determine design iceberg collision events with a structure based on either iceberg kinetic energy upon impact or global sliding force acting on the structure. Load exceedence probabilities from the kinetic energy and sliding force models are evaluated using the second-order reliability method. Output from the probabilistic models can be used to determine design collision parameters and to assess whether more sophisticated modelling of various impact processes is required. The influence of the structure on velocities of approaching icebergs is shown to be significant when the structure horizontal dimension is greater than twice the iceberg dimension. As expected, wave-induced oscillatory motions dominate the collision velocity for smaller icebergs but have a negligible effect on velocity for larger icebergs. / Applied Science, Faculty of / Civil Engineering, Department of / Graduate

Thèse (D.Ing.) -- Université du Québec à Chicoutimi, 2005. / Bibliogr.: f. 133-137. Document électronique également accessible en format PDF. CaQCU

Les phénomènes liés au comportement à la rupture de la glace sous impact sont fréquents dans le génie civil, pour les structures offshore, et les processus de dégivrage. Pour réduire les dommages causés par l'impact de la glace et optimiser la conception des structures ou des machines, l'étude sur le comportement à la rupture dynamique de la glace sous impact est nécessaire. Ces travaux de thèse portent donc sur la propagation dynamique des fissures dans la glace sous impact. Une série d'expériences d'impact est réalisée avec un dispositif de barres de Hopkinson. La température est contrôlée par une chambre de refroidissement. Le processus dynamique de la rupture de la glace est enregistré avec une caméra à grande vitesse et ensuite analysé par des méthodes d'analyse d'images. La méthode des éléments finis étendus complète cette analyse pour évaluer la ténacité dynamique. Au premier abord, le comportement dynamique de la glace sous impact est étudié avec des échantillons cylindriques afin d'établir la relation contrainte-déformation dynamique qui sera utilisée dans les simulations numériques plus tard. Nous avons observé de multi-fissuration dans les expériences sur les échantillons cylindriques mais son étude est trop difficile à mener. Pour mieux comprendre la propagation des fissures dans la glace, des échantillons rectangulaires avec une pré-fissure sont employés. En ajustant la vitesse d'impact on aboutit à la rupture des spécimens avec une fissure principale à partir de la pré-fissure. L'histoire de la propagation de fissure et de sa vitesse sont évaluées par analyse d'images basée sur les niveaux de gris et par corrélation d'images. La vitesse de propagation de la fissure principale est identifiée dans la plage de 450 à 610 m/s ce qui confirme les résultats précédents. Elle varie légèrement au cours de la propagation, dans un premier temps elle augmente et se maintient constante ensuite et diminue à la fin. Les paramètres obtenus expérimentalement, tels que la vitesse d'impact et la vitesse de propagation de fissure, sont utilisés pour la simulation avec la méthode des éléments finis étendus. La ténacité d'initiation dynamique et la ténacité dynamique en propagation de fissure sont déterminées lorsque la simulation correspond aux expériences. Les résultats indiquent que la ténacité dynamique en propagation de fissure est linéaire vis à vis de la vitesse de propagation et semble indépendante de la température dans l'intervalle -15 à -1 degrés. / The phenomena relating to the fracture behaviour of ice under impact loading are common in civil engineering, for offshore structures, and de-ice processes. To reduce the damage caused by ice impact and to optimize the design of structures or machines, the investigation on the dynamic fracture behaviour of ice under impact loading is needed. This work focuses on the dynamic crack propagation in ice under impact loading. A series of impact experiments is conducted with the Split Hopkinson Pressure Bar. The temperature is controlled by a cooling chamber. The dynamic process of the ice fracture is recorded with a high speed camera and then analysed by image methods. The extended finite element method is complementary to evaluate dynamic fracture toughness at the onset and during the propagation. The dynamic behaviour of ice under impact loading is firstly investigated with cylindrical specimen in order to obtain the dynamic stress-strain relation which will be used in later simulation. We observed multiple cracks in the experiments on the cylindrical specimens but their study is too complicated. To better understand the crack propagation in ice, a rectangular specimen with a pre-crack is employed. By controlling the impact velocity, the specimen fractures with a main crack starting from the pre-crack. The crack propagation history and velocity are evaluated by image analysis based on grey-scale and digital image correlation. The main crack propagation velocity is identified in the range of 450 to 610 m/s which confirms the previous results. It slightly varies during the propagation, first increases and keeps constant and then decreases. The experimentally obtained parameters, such as impact velocity and crack propagation velocity, are used for simulations with the extended finite element method. The dynamic crack initiation toughness and dynamic crack growth toughness are determined when the simulation fits the experiments. The results indicate that the dynamic crack growth toughness is linearly associated with crack propagation velocity and seems temperature independent in the range -15 to -1 degrees.