Hostname: page-component-586b7cd67f-t8hqh Total loading time: 0 Render date: 2024-11-25T02:10:06.936Z Has data issue: false hasContentIssue false

Initial experiments on the effects of particles at grain boundaries on the anelasticity and creep behavior of granular ice

Published online by Cambridge University Press:  14 September 2017

Min Song
Affiliation:
Thayer School of Engineering, Dartmouth College, HanoverNH 03755-8000, USA E-mail: Min.Song@Dartmouth.edu
David M. Cole
Affiliation:
US Army Cold Regions Research and Engineering Laboratory, 72 Lyme Road, HanoverNH 03755-1290, USA
Ian Baker
Affiliation:
Thayer School of Engineering, Dartmouth College, HanoverNH 03755-8000, USA E-mail: Min.Song@Dartmouth.edu
Rights & Permissions [Opens in a new window]

Abstract

Experimental observations of the influence of particles at grain boundaries on the anelasticity and creep behavior of granular fresh-water ice are presented. Ice with particle contents of 0–4 wt.% was investigated under both reversed direct-stress and creep loading conditions at –12˚C. The results show that the particles decreased the grain-boundary relaxation by suppressing grain-boundary sliding at higher frequencies (10–1 to 101 Hz). In addition, the modulus increased by up to 30%, and the internal friction decreased by up to 30% at a frequency of 1 Hz. Staged creep tests showed that the particles affected the creep rate substantially. The minimum creep rate of ice containing 1 wt.% particles is 40% higher than that of particle-free ice, indicating that mechanisms besides dislocation glide aid the creep deformation.

Type
Research Article
Copyright
Copyright © The Author(s) [year] 2004

Introduction

The mechanical properties of fresh-water polycrystalline ice have been investigated by a number of workers (e.g. Reference GlenGlen, 1955; Reference DuvalDuval, 1978; Reference JonesJones, 1982; Reference Duval and AshbyDuval and others; 1983; Reference WeertmanWeertman, 1983; Reference JackaJacka, 1984a, Reference Jackab; Reference JackaJacka and Maccagnan, 1984). Recent decades have seen substantial progress in the understanding of the processes involved in glacier flow (Reference NyeNye, 1957; Reference RigsbyRigsby, 1958; Reference RaymondRaymond, 1973; Reference Russell-Head and BuddRussell-Head and Budd, 1979; Reference Paterson and BuddPaterson and Budd, 1982; Reference PatersonPaterson, 1983; Reference Shoji and LangwayShoji and Langway, 1985). However, a full understanding of the constitutive behavior of clean and sediment-laden ice has proven elusive. Understanding how fine particles affect the mechanical behavior of ice is important for the study of glacier flow since low concentrations of debris are common in glacier ice. Although the importance of dislocations in ice deformation is well recognized, a dislocation-based model has been problematic because of the difficulties in quantifying dislocation density. Recently, a method combining cyclic loading and creep experiments has been developed that provides a way to track the dislocation density as a function of the particle content, stress/strain path and thermal history (Reference ColeCole, 1995; Reference Cole and DurellCole and Durell, 2001). Using this approach, Cole and Reference Cole and DurellDurell (2001) found that the flow rate is dislocation drag-controlled for both sea ice and freshwater ice, so that the strain rate is a linear function of the dislocation density. An n = 3 power-law behavior was observed when the applied creep stress controls the effective dislocation density, and n = 1 was found when the dislocation density is effectively independent of stress.

The anelasticity behavior of both fresh-water ice and sea ice has also been widely studied (e.g. Reference VaudreyVaudrey, 1977; Reference DuvalDuval, 1978; Reference Vassoille, Maï and PerezVassoille and others, 1978; Reference ColeCole, 1990, Reference Cole1995, Reference Cole1998; Reference Cole and DurellCole and Durell, 1995). In general, ice exhibits a high level of anelasticity under a wide range of conditions. Previous work (Reference ColeCole, 1995) indicated that both mobile dislocation processes and grain-boundary sliding contribute significantly to ice anelasticity. But how fine particles affect anelasticity, or mobile dislocation processes and grain-boundary sliding, is still unknown. Some studies (Reference Mosher and RajMosher and Raj, 1974; Reference Mori, Kos and MonzenMori and others, 1983; Reference ShigenakaShigenaka and others, 1983) on metals have shown that precipitates at grain boundaries will block grain-boundary sliding and, thus, inhibit grain-boundary relaxation processes. However, there are no studies of the effects of particles at grain boundaries on the anelasticity of polycrystalline ice. One reason is the difficulty in producing ice with the particles spread uniformly along the grain boundaries.

A large number of measurements have been made on the strength of ice with solid inclusions such as sand, silt and clay. However, understanding of the effect of impurities on the rate-controlling deformation process is still incomplete. Reference Hooke, Dahlin and KauperHooke and others (1972) found a high variability in the strain rate for a range of sand concentrations from 0 to 10 vol.%, Reference Baker and Gerberichwhile Baker and Gerberich (1979) found that at temperatures from –20 to –5˚C solid inclusions, in concentrations from 1.3 to 6.6 vol.%, increased the creep strain rate. Later, Reference Shoji and LangwayShoji and Langway (1985) found an increased flow rate for dirty basal ice from Camp Century, Greenland, compared with clean basal ice. Reference Holdsworth and BullHoldsworth and Bull (1970) also observed an enhanced flow rate in the basal amber-colored ice of Meserve Glacier, Antarctica. However, these studies did not clarify the effect of impurities at different microstructural locations, i.e. whether having particles located at grain boundaries or in the grain interiors has the same effect on the creep rate.

As a first step towards understanding the importance of particle location to the flow processes, this study focuses on the effects of particles located at grain boundaries on both the anelasticity and creep behavior of granular ice.

Experimental

Granular fresh-water ice specimens and debris-laden ice specimens (0.1 wt.%, 0.5 wt.%, 1 wt.%, 4 wt.%) with grain-sizes of 3–5mm containing silt-sized particles (average particle size of 0.05 mm) along the grain boundaries were prepared by the method invented by Reference ColeCole (1979). The specimen dimensions were 127mm in length and 50.4 mm in diameter.

Anelasticity experiments employed the reversed direct-stress testing method presented by Reference ColeCole and Durell (1995), and the procedures described by Reference ColeCole (1993). The command signal was a zero-mean-stress sinusoid. A load cell mounted in line with the specimen monitored the axial load, and two displacement transducers mounted on the specimen were used to monitor the displacements. Digitized records provided 200 points per load cycle regardless of frequency. The raw data are smoothed using procedures described by Reference Cole and DurellCole and Durell (1995), and the internal friction and elastic modulus were calculated from the smoothed points.

Staged creep tests, which consisted of alternating periods of uniaxial cyclic and creep loading, were used to examine the creep behavior. As a first step to determine the effects of particles along the grain boundaries, a rather high stress level of 1.45 MPa was used for creep tests at –12˚C, so that a relationship between creep and particle content could be found quickly. Low-stress-level tests are being performed at the time of writing. A set of load cycles was applied prior to creep deformation to assess the initial state of the specimens. Then a constant creep load was applied, followed by load removal, anelastic recovery and sets of loading cycles. This sequence was repeated several times for a given creep stress until the total viscous strain was in the range 0.06–0.08. Because the dislocation density is directly related to the hysteresis loop area, the cyclic loading response is used to track changes in the dislocation density after each creep loading. The cyclic loading frequency used was 0.1 Hz, which was sufficiently low to produce measurable anelastic strain, but high enough to avoid cyclic viscous strain.

The specimens before and after creep were cut into thin sections for microstructural analysis using optical microscopy with polarized light.

Results and Discussion

Anelasticity

Figure 1 shows the relationship between modulus, internal friction and frequency of cyclic loading for particle-free ice and the debris-laden ice with the particles at the grain boundaries. It is clearly seen that the modulus, which increases with increasing frequency, is similar (variation <10%) for all specimens at the low frequencies of 0.001 and 0.01 Hz, but an increase in modulus can be seen for the debris-laden ice when the loading frequency is ∼0.1Hz. At 1 Hz, the modulus of the debris-laden ice is 30% higher than that of particle-free ice. Consistent with the increases in modulus, the internal friction of debris-laden ice decreases compared with particle-free ice.

Fig. 1. Modulus and internal friction as function of cyclic loading frequency for particle-free ice and debris-laden ice (before creep).

It is known that dislocation relaxation can be an order of magnitude greater than grain-boundary relaxation, and that the grain-boundary relaxation peak occurs at ∼5Hz while the dislocation relaxation peak occurs at about ∼5×10–4Hz at –12˚C (Reference ColeCole, 1995). Thus, the modulus increase at 1 Hz in the debris-laden ice is due to the particles decreasing the grain-boundary relaxation. Reference ShigenakaShigenaka and others (1983) also observed this phenomenon in metals (Cu-Fe and Cu-SiO2). When the particles are located at grain boundaries, they retard grain-boundary sliding through the development of internal stresses, which decrease the free sliding distance and the relaxation time, thus increasing the modulus and decreasing the internal friction. The physical basis for dislocation relaxation at low frequency is the oscillatory motion of basal plane dislocations under a sinusoidally varying stress. This will not be affected by the particles at the grain boundaries. Hence, the modulus for particle-free ice and debris-laden ice at low frequencies (0.001 and 0.01 Hz) is unaffected by grain-boundary particles.

Creep

Figure 2 shows the relationship between the minimum creep rate and the particle content. It can be seen that the minimum creep rate increases for low particle contents up to 1 wt.%, after which it decreases. Figure 3 shows the relationship between the creep rate and viscous strain for fresh-water ice and debris-laden ice with 0.5 wt.%, 1 wt.% and 4 wt.% particles. All the curves show a similar trend: the creep rate initially decreases (primary creep regime), reaches a minimum near 1% strain (secondary creep) and subsequently increases as it enters tertiary creep. It is clear that the strain rate for debris-laden ice with 1 wt.% particles at the grain boundaries is higher than the strain rate for both particle-free ice and debris-laden ice with 0.5 wt.% and 4 wt.% particles at the grain boundaries. Reference Hooke, Dahlin and KauperHooke and others (1972) also found an increased creep rate for ice containing a low content of dispersed fine sand. Figure 4 shows the calculated dislocation density for particle-free ice and debris-laden ice during creep, calculated using the model developed by Cole and Reference Cole and DurellDurell (2001).

Fig. 2. Minimum strain rate vs debris content for specimens with particles at the grain boundaries.

Fig. 3. Creep rate vs viscous strain for particle-free ice and debris-laden ice with particles at the grain boundaries (stress level 1.45 MPa).

Fig. 4. Calculated dislocation density vs viscous strain.

(1)

where ω is the angular frequency, τ is the relaxation time, is the loss compliance, αd is the peak width parameter for the dislocation relaxation, δDd is the strength of the relaxation, K is the restoring stress constant, ρ is the dislocation density, Ω is the orientation factor and b is the Burgers vector. We see that except for some fluctuations during primary creep, the dislocation density increases as the creep rate increases. The calculated mobile dislocation densities of the specimens with 0.5 wt.% and 1 wt.% particles are lower than those of both the particle-free specimen and the specimen with 4 wt.% particles. The reason for this is still unclear. Figure 5 shows the fitted linear relationship between the dislocation density and the creep rate. The fluctuation in dislocation density at low strains could be due to dynamic recovery during primary creep. When the creep rate subsequently increases, the dislocation density also increases. Interestingly, for debris-laden ice with 1 wt.% particles, the dislocation density is lower but the creep rate is higher than for particle-free ice and debris-laden ice with 4 wt.% particles. This indicates that another creep mechanism is possibly aiding creep besides dislocation glide for the 1 wt.% debris-laden ice. One possibility may be that the particles at the grain boundaries act as preferential nucleation sites for recrystallization. This will increase the overall nucleation rate, and many small new grains will be nucleated along the grain boundary, increasing the grain-boundary area. These new recrystallized grains have no particles along their boundaries, and the boundaries are thus free to slide. Reference Baker and GerberichBaker and Gerberich (1979) observed this phenomenon for ice with dispersed solid inclusions.

Fig. 5. Relationship between dislocation density and viscous strain. The lines were fitted using linear regression.

Microstructure

Figure 6 shows the microstructures of ice before and after creep deformation for debris-laden ice with 1 wt.% particles. It is clearly seen that dynamic recrystallization occurred during creep, a phenomenon observed in all the specimens. Hence, the average grain-size decreased after creep, compared with the unstrained grain-size. Figure 7 shows how the content of particles at the grain boundary affected recrystallization. Numerous new grains were nucleated at the grain boundaries and grew by consuming the deformed grains. It is seen that the recrystallized grain-size of the debris-laden ice with 1 wt.% particles is smaller than that of the other specimens. Thus, the debris-laden ice with 1 wt.% particles has a larger grain-boundary area, and these new grain boundaries have no particles to inhibit boundary sliding.

Fig. 6. Microstructures (a) before and (b) after staged-creep test for ice with 1 wt.% particles (strain 0.07).

Fig. 7. Microstructures after staged-creep test for debris-laden ice with different particle contents. (a) Particle-free ice (strain 0.06); (b) 0.5 wt.% particles ice (strain 0.075); (c) 1 wt.% particles ice (strain 0.07); and (d) 4 wt.% particles ice (strain 0.075).

Summary and Conclusions

A combination of cyclic and constant compression loading was used to study the anelasticity and creep behavior of debris-laden ice and particle-free ice. The results show that the particles increase the modulus at a frequency of 1 Hz but do not affect it at low frequencies. The particles, which are located at the grain boundaries, inhibit grain-boundary sliding during the cyclic loading, which decreases the grain-boundary relaxation time or even stops grain-boundary relaxation completely at 1 Hz.

Staged creep tests show that the dislocation density increases with increasing creep rate for all the specimens. This phenomenon indicates that dislocation glide is the major mechanism of creep for the prevailing conditions. However, the high creep rate observed for the 1 wt.% debris-laden ice without a corresponding increase in dislocation density, compared with particle-free ice, shows that another creep mechanism might also aid creep deformation. Small new recrystallized grains increase the grain-boundary area, which may increase the boundary-sliding contribution to creep.

Acknowledgements

This research was supported by the US National Science Foundation Office of Polar Programs, Arctic Natural Sciences Program (OPP 011737; program manager J. Dionne). We thank G. Durell for his valuable assistance in developing the creep and cyclic loading equipment.

References

Baker, R.W. and Gerberich, W.W. 1979. The effect of crystal size and dispersed-solid inclusions on the activation energy for creep of ice. J. Glaciol., 24(90), 179194.Google Scholar
Cole, D.M. 1979. Preparation of polycrystalline ice specimens for laboratory experiments. Cold Reg. Sci. Technol., 1(2), 153159.Google Scholar
Cole, D.M. 1990. Reversed direct-stress testing of ice: initial experimental results and analysis. Cold Reg. Sci. Technol., 18(3), 303321.Google Scholar
Cole, D.M. 1993. The effect of creep on the constitutive behavior of saline ice at low temperature. In Ice mechanics. Vol. AMD-163,. New York, American Society of Mechanical Engineers. Applied Mechanics Division, 261271.Google Scholar
Cole, D.M. 1995. A model for the anelastic straining of saline ice subjected to cyclic loading. Philos. Mag. A, 72(1), 231248.Google Scholar
Cole, D.M. 1998.Modeling the cyclic loading response of sea ice. Int. J. Solids Struct, 35(31–32), 40674075.Google Scholar
Cole, D.M. and Durell, G.D. 1995. The cyclic loading of saline ice. Philos. Mag. A, 72(1), 209229.Google Scholar
Cole, D.M. and Durell, G.D. 2001. A dislocation-based analysis of strain history effects in ice. Philos. Mag. A, 81(7), 18491872.Google Scholar
Duval, P. 1978. Anelastic behaviour of polycrystalline ice. J. Glaciol., 21(85), 621628.Google Scholar
Duval, P., Ashby, M.F. and I. Anderman. 1983. Rate-controlling processes in the creep of polycrystalline ice. J. Phys. Chem., 87(21), 40664074.Google Scholar
Glen, J.W. 1955. The creep of polycrystalline ice. Proc. R. Soc. London, Ser. A, 228(1175), 519538.Google Scholar
Holdsworth, G. and Bull, C. 1970. The flow of cold ice: investigations on Meserve Glacier, Antarctic. International Association of Scientific Hydrology Publication, 86 (Symposium at Hanover 1968 – Antarctic Glaciological Exploration (ISAGE)), 204216.Google Scholar
Hooke, R. LeB., Dahlin, B.B. and Kauper, M.T. 1972. Creep of ice containing dispersed fine sand. J. Glaciol., 11(63), 327336.Google Scholar
Jacka, T.H. 1984a. Laboratory studies on relationships between ice crystal size and flow rate. Cold Reg. Sci. Technol.,, 10(1), 3142.Google Scholar
Jacka, T.H. 1984b. The time and strain required for development of minimum strain rates in ice. Cold Reg. Sci. Technol.,, 8(3), 261268.CrossRefGoogle Scholar
Jacka, T.H. and M. Maccagnan. 1984. Ice crystallographic and strain rate changes with strain in compression and extension. Cold Reg. Sci. Technol., 8(3), 269286.Google Scholar
Jones, S.J. 1982. The confined compressive strength of polycrystalline ice. J. Glaciol., 28(98), 171177.Google Scholar
Mori, T., Kos, M. and Monzen, R. 1983. Particle blocking in grain boundary sliding and associated internal friction. Acta Metall., 31(2), 275283.Google Scholar
Mosher, D.R. and Raj, J. (1974) Use of the internal friction technique to measure rates of grain boundary slidin. Acta Metall., 2(2), 21469 Google Scholar
Nye, J.F. 1957. The distribution of stress and velocity in glaciers and ice-sheets. Proc. R. Soc. London, Ser. A, 239(1216), 113133.Google Scholar
Paterson, W. S. B. 1983. Deformation within polar ice sheets: an analysis of the Byrd Station and Camp Century borehole-tilting measurements. Cold Reg. Sci. Technol., 8(2), 165179.Google Scholar
Paterson, W. S. B. and Budd, W.F. 1982. Flow parameters for ice sheet modelling. Cold Reg. Sci. Technol., 6(2), 175177.Google Scholar
Raymond, C.F. 1973. Inversion of flow measurements for stress and rheological parameters in a valley glacier. J. Glaciol., 12(64), 1944.CrossRefGoogle Scholar
Rigsby, G.P. 1958. Effect of hydrostatic pressure on velocity of shear deformation on single ice crystals. J. Glaciol., 3(24), 273278.271–272.Google Scholar
Russell-Head, D.S. and Budd, W.F. 1979. Ice-sheet flow properties derived from bore-hole shear measurements combined with ice-core studies. J. Glaciol., 24(90), 117130.CrossRefGoogle Scholar
Shigenaka, N., R. Monzen and T. Mori. 1983. Internal friction caused by sliding on grain boundaries with second phase particles Cu-Fe and Cu-SiO2 . Acta Metall., 31(12), 20872093.Google Scholar
Shoji, H. and Langway, C.C. Jr., 1985. Comparison of mechanical test on the Dye-3, Greenland ice core and artificial laboratory ice. Ann. Glaciol.,, 6, 305.Google Scholar
Vassoille, R., Maï, C. and Perez, J. 1978. Inelastic behaviour of ice Ih single crystals in the low-frequency range due to dislocations. J. Glaciol., 21(85), 375384.Google Scholar
Vaudrey, K.D. (1977) Ice engineering – study of related properties of floating sea-ice sheets and summary of elastic and viscoelastic analyses. Port Hueneme, CA Naval Construction Battalion Center. Civil Engineering Laboratory. (Technical Report R 860).Google Scholar
Weertman, J. (1983) Creep deformation of ic. Annu. Rev. Earth Planet. Sci., 1(1), 1215 Google Scholar
Figure 0

Fig. 1. Modulus and internal friction as function of cyclic loading frequency for particle-free ice and debris-laden ice (before creep).

Figure 1

Fig. 2. Minimum strain rate vs debris content for specimens with particles at the grain boundaries.

Figure 2

Fig. 3. Creep rate vs viscous strain for particle-free ice and debris-laden ice with particles at the grain boundaries (stress level 1.45 MPa).

Figure 3

Fig. 4. Calculated dislocation density vs viscous strain.

Figure 4

Fig. 5. Relationship between dislocation density and viscous strain. The lines were fitted using linear regression.

Figure 5

Fig. 6. Microstructures (a) before and (b) after staged-creep test for ice with 1 wt.% particles (strain 0.07).

Figure 6

Fig. 7. Microstructures after staged-creep test for debris-laden ice with different particle contents. (a) Particle-free ice (strain 0.06); (b) 0.5 wt.% particles ice (strain 0.075); (c) 1 wt.% particles ice (strain 0.07); and (d) 4 wt.% particles ice (strain 0.075).