Numerical Simulations of the Drifting Ice Sheets Collision with the Bridge Pier

2013 ◽  
Vol 368-370 ◽  
pp. 1383-1386
Author(s):  
Lian Zhen Zhang ◽  
Wei Xiong

The drifting ice sheets impact with the bridge pier and other hydraulic structures in the rivers, which may damage even cause collapse of the structures. In this paper, the FEM software package LS-DYNA was used to performed the numerical simulations of the collision process of the ice sheets and the bridge piers to make clear the interaction between them and to understand the failure mechanism of the ice sheet. The elastic strain-stress model with von mises failure criterion was used to describe the ice material. The brittle damage material model was used to describe the concrete pier. Three types thickness of ice sheets were performed at various velocity of the ice sheet respectively. The impact process of every case were displayed and the time history curve of the collision force were given out. The simulations results show that the peak value of the collision force time history curve increases with the velocity of the sheet firstly and then decreases with the velocity of the ice sheet. There is one critical velocity which relate to the compressive strength of the ice sheet. The simulation result were also compared with the different bridge design code, which show that the code result is more conservative in bridge design.

2016 ◽  
Vol 12 (12) ◽  
pp. 2195-2213 ◽  
Author(s):  
Heiko Goelzer ◽  
Philippe Huybrechts ◽  
Marie-France Loutre ◽  
Thierry Fichefet

Abstract. As the most recent warm period in Earth's history with a sea-level stand higher than present, the Last Interglacial (LIG,  ∼  130 to 115 kyr BP) is often considered a prime example to study the impact of a warmer climate on the two polar ice sheets remaining today. Here we simulate the Last Interglacial climate, ice sheet, and sea-level evolution with the Earth system model of intermediate complexity LOVECLIM v.1.3, which includes dynamic and fully coupled components representing the atmosphere, the ocean and sea ice, the terrestrial biosphere, and the Greenland and Antarctic ice sheets. In this setup, sea-level evolution and climate–ice sheet interactions are modelled in a consistent framework.Surface mass balance change governed by changes in surface meltwater runoff is the dominant forcing for the Greenland ice sheet, which shows a peak sea-level contribution of 1.4 m at 123 kyr BP in the reference experiment. Our results indicate that ice sheet–climate feedbacks play an important role to amplify climate and sea-level changes in the Northern Hemisphere. The sensitivity of the Greenland ice sheet to surface temperature changes considerably increases when interactive albedo changes are considered. Southern Hemisphere polar and sub-polar ocean warming is limited throughout the Last Interglacial, and surface and sub-shelf melting exerts only a minor control on the Antarctic sea-level contribution with a peak of 4.4 m at 125 kyr BP. Retreat of the Antarctic ice sheet at the onset of the LIG is mainly forced by rising sea level and to a lesser extent by reduced ice shelf viscosity as the surface temperature increases. Global sea level shows a peak of 5.3 m at 124.5 kyr BP, which includes a minor contribution of 0.35 m from oceanic thermal expansion. Neither the individual contributions nor the total modelled sea-level stand show fast multi-millennial timescale variations as indicated by some reconstructions.


2021 ◽  
Author(s):  
Sam Sherriff-Tadano ◽  
Ayako Abe-Ouchi ◽  
Akira Oka ◽  
Takahito Mitsui ◽  
Fuyuki Saito

Abstract. Glacial periods undergo frequent climate shifts between warm interstadials and cold stadials on a millennial time-scale. Recent studies have shown that the duration of these climate modes varies with the background climate; a colder background climate and lower CO2 generally results in a shorter interstadial and a longer stadial through its impact on the Atlantic Meridional Overturning Circulation (AMOC). However, the duration of stadials was shorter during the Marine Isotope Stage 3 (MIS3) compared with MIS5, despite the colder climate in MIS3, suggesting potential control from other climate factors on the duration of stadials. In this study, we investigated the role of glacial ice sheets. For this purpose, freshwater hosing experiments were conducted with an atmosphere–ocean general circulation model under MIS5a, MIS3 and MIS3 with MIS5a ice sheet conditions. The impact of ice sheet differences on the duration of the stadials was evaluated by comparing recovery times of the AMOC after freshwater forcing was reduced. Hosing experiments showed a slightly shorter recovery time of the AMOC in MIS3 compared with MIS5a, which was consistent with ice core data. We found that larger glacial ice sheets in MIS3 shortened the recovery time. Sensitivity experiments showed that stronger surface winds over the North Atlantic shortened the recovery time by increasing the surface salinity and decreasing the sea ice amount in the deepwater formation region, which set favourable conditions for oceanic convection. In contrast, we also found that surface cooling by larger ice sheets tended to increase the recovery time of the AMOC by increasing the sea ice thickness over the deepwater formation region. Thus, this study suggests that the larger ice sheet in MIS3 compared with MIS5a could have contributed to the shortening of stadials in MIS3, despite the climate being colder than that of MIS5a, when the effect of surface wind played a larger role.


2012 ◽  
Vol 6 (6) ◽  
pp. 4897-4938 ◽  
Author(s):  
S. Charbit ◽  
C. Dumas ◽  
M. Kageyama ◽  
D. M. Roche ◽  
C. Ritz

Abstract. Since the original formulation of the positive-degree-day (PDD) method, different PDD calibrations have been proposed in the literature in response to the increasing number of observations. Although these formulations provide a satisfactory description of the present-day Greenland geometry, they have not all been tested for paleo ice sheets. Using the climate-ice sheet model CLIMBER-GRISLI coupled with different PDD models, we evaluate how the parameterization of the ablation may affect the evolution of Northern Hemisphere ice sheets in the transient simulations of the last glacial cycle. Results from fully coupled simulations are compared to time-slice experiments carried out at different key periods of the last glacial period. We find large differences in the simulated ice sheets according to the chosen PDD model. These differences occur as soon as the onset of glaciation, therefore affecting the subsequent evolution of the ice system. To further investigate how the PDD method controls this evolution, special attention is given to the role of each PDD parameter. We show that glacial inception is critically dependent on the representation of the impact of the temperature variability from the daily to the inter-annual time scale, whose effect is modulated by the refreezing scheme. Finally, an additional set of sensitivity experiments has been carried out to assess the relative importance of melt processes with respect to initial ice sheet configuration in the construction and the evolution of past Northern Hemisphere ice sheets. Our analysis reveals that the impacts of the initial ice sheet condition may range from quite negligible to explaining about half of the LGM ice volume depending on the representation of stochastic temperature variations which remain the main driver of the evolution of the ice system.


1983 ◽  
Vol 4 ◽  
pp. 180-187 ◽  
Author(s):  
B. Michel ◽  
D. Blanchet

The problem of a floating ice sheet hitting a structure with a vertical face appears to be a simple one but, in fact, has only been solved for a limited number of cases. Research work on this question usually reports on an indentation coefficient which relates the average pressure on the indenter to the uniaxial crushing strength of the ice. Very few tests have been made in the brittle range of ice failure. In this particular area of study, this paper reports on 27 tests that were conducted in a cold-room water basin where controlled S2floating ice sheets were produced with a surface area of 4 × 4 m, three sides being fully restrained and the other, freely float! no, being submitted to the impact of the moving indenter. All tests were carried out at computed indentation rates varying from 0.017 to 0.34 s-1. In this range this ice would normally be considered to act as a brittle material. The thickness of the ice sheets varied from 1.2 to 9.0 cm and the indenter width from 5 cm to 1 m. Overall, the aspect ratio relating these two parameters could be varied from 0.5 to 83.Results have shown that for aspect ratios <5, there was an important oscillatory effect which caused the formation of pi asti fi ed triangles in front of the indenter, increasing its resistance as it would under ductile conditions. Because of this plastification, an extrusion effect appeared in front of the indenter as the broken ice crystals were blown up and down in front of the fast-moving indenter. The theory of plasticity which gives an indentation coefficient of 2.97 seems to apply in this case. Another mode of failure which occurred with aspect ratios 5 was cleavage in the plane of the ice sheet which also gives a higher indentation coefficient for S2ice, but of the same order of magnitude as previously.For intermediate values of the aspect ratio, between 5 and 20, the theory of elasticity used by Michel (1978) seems to apply well. Shear cracks are first formed on both sides of the square indenter and control the maximum pressure when they propagate inside forming big triangles in front of it.Finally, for aspect ratios ~>20, buckling of the ice occurs, either after or at the same time as the formation of wedges, together with a reduction in the indentation coefficient to a value close to that given by the theory of buckling of a truncated 45° wedge with a hinged edge.


2020 ◽  
Author(s):  
Thomas Kleiner ◽  
Jeremie Schmiedel ◽  
Angelika Humbert

<p>Ice sheets constitute the largest and most uncertain potential source of future sea-level rise. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) brings together a consortium of international ice sheet and climate models to explore the contribution from the Greenland and Antarctic ice sheets to future sea-level rise.</p> <p>We use the Parallel Ice Sheet Model (PISM, pism-docs.org) to carry out spinup and projection simulations for the Antarctic Ice Sheet. Our treatment of the ice-ocean boundary condition previously based on 3D ocean temperatures (initMIP-Antarctica) has been adopted to use the ISMIP6 parameterisation and 3D ocean forcing fields (temperature and salinity) according to the ISMIP6 protocol.</p> <p>In this study, we analyse the impact of the choices made during the model initialisation procedure on the initial state. We present the AWI PISM results of the ISMIP6 projection simulations and investigate the ice sheet response for individual basins. In the analysis, we distinguish between the local and non-local ice shelf basal melt parameterisation.</p>


2020 ◽  
Author(s):  
Felicity McCormack ◽  
Roland Warner ◽  
Adam Treverrow ◽  
Helene Seroussi

<p>Viscous deformation is the main process controlling ice flow in ice shelves and in slow-moving regions of polar ice sheets where ice is frozen to the bed. However, the role of deformation in flow in ice streams and fast-flowing regions is typically poorly represented in ice sheet models due to a major limitation in the current standard flow relation used in most large-scale ice sheet models – the Glen flow relation – which does not capture the steady-state flow of anisotropic ice that prevails in polar ice sheets. Here, we highlight recent advances in modeling deformation in the Ice Sheet System Model using the ESTAR (empirical, scalar, tertiary, anisotropic regime) flow relation – a new description of deformation that takes into account the impact of different types of stresses on the deformation rate. We contrast the influence of the ESTAR and Glen flow relations on the role of deformation in the dynamics of Thwaites Glacier, West Antarctica, using diagnostic simulations. We find key differences in: (1) the slow-flowing interior of the catchment where the unenhanced Glen flow relation simulates unphysical basal sliding; (2) over the floating Thwaites Glacier Tongue where the ESTAR flow relation outperforms the Glen flow relation in accounting for tertiary creep and the spatial differences in deformation rates inherent to ice anisotropy; and (3) in the grounded region within 80km of the grounding line where the ESTAR flow relation locally predicts up to three times more vertical shear deformation than the unenhanced Glen flow relation, from a combination of enhanced vertical shear flow and differences in the distribution of basal shear stresses. More broadly on grounded ice, the membrane stresses are found to play a key role in the patterns in basal shear stresses and the balance between basal shear stresses and gravitational forces simulated by each of the ESTAR and Glen flow relations. Our results have implications for the suitability of ice flow relations used to constrain uncertainty in reconstructions and projections of global sea levels, warranting further investigation into using the ESTAR flow relation in transient simulations of glacier and ice sheet dynamics. We conclude by discussing how geophysical data might be used to provide insight into the relationship between ice flow processes as captured by the ESTAR flow relation and ice fabric anisotropy.</p>


2008 ◽  
Vol 22 (09n11) ◽  
pp. 1361-1368 ◽  
Author(s):  
GOUPING REN ◽  
ZHU LI

The impact test on short concrete filled steel tubular column was conducted through DHR-9401 dropped hammer tester. Based on analysis on recorded time-history curve of impact force, the relations of impact force with respect to confining effect coefficient and impact speed are obtained. So are done that of the impact duration. By use of ANSYS/LS-DYNA, values of impact peak force in relation with those of impact speed were computed in the case of unitary material model and composite material model respectively. The simulation results show that peak force-speed curve of unitary model has better description of test data than that of composite one. Critical energy is found to increase linearly with the steel ratio when steel tube and concrete remain unchanged.


2013 ◽  
Vol 871 ◽  
pp. 240-246
Author(s):  
Tae Ho Yang ◽  
Young Shin Lee ◽  
Yoon Jae Kim ◽  
Tae Hyeong Kim ◽  
Chang Won Shul ◽  
...  

Polyurethane (PU) S80A was used as the material of the elastomer of the shock programmer in this paper. To validate Ogden hyper-elastic material model in simulation, the small impact test was performed. As the comparison for the time history of the acceleration between the impact test and simulation was performed. Using the cylindrical shock programmer, the constant used in Ogden hyper-elastic material model was calculated. The wave shape of the acceleration was obtained with the noised sign. To clearly obtain the wave shape of the acceleration the cylindrical shock programmer, the truncated conical shock programmer was used. Using the Ogden hyper-elastic material model, design variables of the shock programmer with the truncated conical shape was studied. Using the shock programmer with truncated conical shape the range on the level and time duration of the acceleration in simulation was from 494.9 m/s2 to 10941 m/s2 and from 1.3 msec to 23.5 msec, respectively.


2021 ◽  
Author(s):  
Nathaelle Bouttes ◽  
Didier Roche ◽  
Fanny Lhardy ◽  
Aurelien Quiquet ◽  
Didier Paillard ◽  
...  

<p>The last deglaciation is a time of large climate transition from a cold Last Glacial Maximum at 21,000 years BP with extensive ice sheets, to the warmer Holocene 9,000 years BP onwards with reduced ice sheets. Despite more and more proxy data documenting this transition, the evolution of climate is not fully understood and difficult to simulate. The PMIP4 protocol (Ivanovic et al., 2016) has indicated which boundary conditions to use in model simulations during this transition. The common boundary conditions should enable consistent multi model and model-data comparisons. While the greenhouse gas concentration evolution and orbital forcing are well known and easy to prescribe, the evolution of ice sheets is less well constrained and several choices can be made by modelling groups. First, two ice sheet reconstructions are available: ICE-6G (Peltier et al., 2015) and GLAC-1D (Tarasov et al., 2014). On top of topographic changes, it is left to modelling groups to decide whether to account for the associated bathymetry and land-sea mask changes, which is technically more demanding. These choices could potentially lead to differences in the climate evolution, making model comparisons more complicated.</p><p>We use the iLOVECLIM model of intermediate complexity (Goosse et al., 2010) to evaluate the impact of different ice sheet reconstructions and the effect of bathymetry changes on the global climate evolution during the Last deglaciation. We test the two ice sheet reconstructions (ICE-6G and GLAC-1D), and have implemented changes of bathymetry and land-sea mask. In addition, we also evaluate the impact of accounting for the Antarctic ice sheet evolution compared to the Northern ice sheets only.</p><p>We show that despite showing the same long-term changes, the two reconstructions lead to different evolutions. The bathymetry plays a role, although only few changes take place before ~14ka. Finally, the impact of the Antarctic ice sheet is important during the deglaciation and should not be neglected.</p><p>References</p><p>Goosse, H., et al., Description of the Earth system model of intermediate complexity LOVECLIM version 1.2, Geosci. Model Dev., 3, 603–633, https://doi.org/10.5194/gmd-3-603-2010, 2010</p><p>Ivanovic, R. F., et al., Transient climate simulations of the deglaciation 21–9 thousand years before present (version 1) – PMIP4 Core experiment design and boundary conditions, Geosci. Model Dev., 9, 2563–2587, https://doi.org/10.5194/gmd-9-2563-2016, 2016</p><p>Peltier, W. R., Argus, D. F., and Drummond, R., Space geodesy constrains ice age terminal deglaciation: The global ICE-6G_C (VM5a) model, J. Geophys. Res.-Sol. Ea., 120, 450–487, doi:10.1002/2014JB011176, 2015</p><p>Tarasov,L.,  et al., The global GLAC-1c deglaciation chronology, melwater pulse 1-a, and a question of missing ice, IGS Symposium on Contribution of Glaciers and Ice Sheets to Sea-Level Change, 2014</p>


1983 ◽  
Vol 4 ◽  
pp. 180-187 ◽  
Author(s):  
B. Michel ◽  
D. Blanchet

The problem of a floating ice sheet hitting a structure with a vertical face appears to be a simple one but, in fact, has only been solved for a limited number of cases. Research work on this question usually reports on an indentation coefficient which relates the average pressure on the indenter to the uniaxial crushing strength of the ice. Very few tests have been made in the brittle range of ice failure. In this particular area of study, this paper reports on 27 tests that were conducted in a cold-room water basin where controlled S2 floating ice sheets were produced with a surface area of 4 × 4 m, three sides being fully restrained and the other, freely float! no, being submitted to the impact of the moving indenter. All tests were carried out at computed indentation rates varying from 0.017 to 0.34 s-1. In this range this ice would normally be considered to act as a brittle material. The thickness of the ice sheets varied from 1.2 to 9.0 cm and the indenter width from 5 cm to 1 m. Overall, the aspect ratio relating these two parameters could be varied from 0.5 to 83.Results have shown that for aspect ratios <5, there was an important oscillatory effect which caused the formation of pi asti fi ed triangles in front of the indenter, increasing its resistance as it would under ductile conditions. Because of this plastification, an extrusion effect appeared in front of the indenter as the broken ice crystals were blown up and down in front of the fast-moving indenter. The theory of plasticity which gives an indentation coefficient of 2.97 seems to apply in this case. Another mode of failure which occurred with aspect ratios 5 was cleavage in the plane of the ice sheet which also gives a higher indentation coefficient for S2 ice, but of the same order of magnitude as previously.For intermediate values of the aspect ratio, between 5 and 20, the theory of elasticity used by Michel (1978) seems to apply well. Shear cracks are first formed on both sides of the square indenter and control the maximum pressure when they propagate inside forming big triangles in front of it.Finally, for aspect ratios ~>20, buckling of the ice occurs, either after or at the same time as the formation of wedges, together with a reduction in the indentation coefficient to a value close to that given by the theory of buckling of a truncated 45° wedge with a hinged edge.


Sign in / Sign up

Export Citation Format

Share Document