Hostname: page-component-78c5997874-4rdpn Total loading time: 0 Render date: 2024-11-17T15:03:21.476Z Has data issue: false hasContentIssue false

On the strain-rate sensitivity of columnar ice

Published online by Cambridge University Press:  20 January 2017

M. E. Manley
Affiliation:
Thayer School of Engineering, Dartmouth College, Hanover, New Hampshire 03755, U.S.A.
E. M. Schulson
Affiliation:
Thayer School of Engineering, Dartmouth College, Hanover, New Hampshire 03755, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

A power law relation between stress and strain rate of the form was used to describe the response to strain rate of S1 ice loaded across the columns at –10° C. The rate exponent, n, decreased with increasing strain from about 4.6 at an observed peak on the load displacement curve to approximately 2.6 at a shortening of 2%. Analysis of these results and of the results of other authors on different forms of ice deformed at the same temperature suggests that the power law exponent, n is proportional to Fc/Fg The parameter Fc/Fg is the far-field basal dislocation climb force divided by the glide force.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1997

During an investigation of the ductile compressive failure of S1 ice (designation of Reference Michel, and Ramseier,Michel and Ramseir, 1971), we performed a series of experiments to determine the strain-rate sensitivity, n, of the flow stress

where B is a constant and έ is the strain rate. Instead of obtaining one value of about n = 3 as one might expect from Glen’s law for power law creep of polycrystals (Reference Duval,, Ashby, and Andermn,Duval and others, 1983; Reference Weertman,Weertman, 1983), we obtained several values whose range seemed too large to ascribe to experimental scatter. This result was reminiscent of earlier observations by Reference Frederking,Frederking (1977) and by Reference Schulson, and Buck,Schulson and Buck (1995). It also brought to mind Reference Weertman,Weertman’s (1983) compilation. Ia an attempt to rationalize this behavior, we found a relationship between the strain-rate sensitivity and the ratio of the climb force to the glide force on basal dislocations. Hence this note.

The Sl ice was grown unidirectional as a sheet from fresh water. The procedure and the material are fully described elsewhere (Reference Manley,Manley. 1996). The ice consisted of columnar-shaped grains about 8 mm in diameter whose crystallographic c axes were parallel to the long axes of the columns to within ±5° . Plate-like specimens (159 mm x 159 mm x 25 mm) were milled from the sheet, with the columns perpendicular to the large faces. Five specimens were then compressed uniaxially across the columns (i.e. along a direction more or less within the basal plane) using a servo-hydraulic loading frame, at a temperature of –10 ± 0.2° C. Each specimen was initially strained at a constant but different rate (

and 10-4s-1) to beyond the peak in the stress–strain curve; i.e. to a strain of about 0.005 as measured using a displacement gauge mounted across the ice. The peak occurred at a strain of about 0.003. Subsequently, and without unloading, the strain rate was either instantaneously increased or decreased by a factor of ten and the ice was shortened an additional 0.2–0.4%, after which the strain rate was instantaneously restored to its initial value and the ice was shortened by this amount again. This cycle was repeated twice before the specimen was unloaded. The strain-rate sensitivity was computed from the peak in the flow stress in the usual manner (i.e. from the slope of a plot of log flow stress vs log strain rate) and from the flow stresses just before (σ1), and just after (σ2) the change in strain rate through the relationship
.

The ice kinked in every experiment, owing to the load being applied in the basal plane. An analysis is given elsewhere (Reference Manley, and Schulson,Manley and Schulson, 1997). This feature, however, did not prevent the rate sensitivity from being determined. The stress–strain curves are given by Reference Manley,Manley (1996). Upon analysis, they showed that the sensitivity decreased from n = 4.6 at an axial strain of 0.003 to n = 2.6 at a strain of 0.02 (Fig. 1). The vertical bars denote the spread in values obtained from the different specimens.

Fig. 1. The exponent, n, from the relation

vs strain, for S1 fresh-water ice deformed at –10° C. The star (*) was obtained from a six tests (Reference Manley,Manley, 1996) from 10−6 to 10−4s−1 where the log of the peak compressive stress was plotted vs the log of the strain rate. All of the other data were obtained from strain-rate change tests.

To rationalize this dependence, we assumed that the ice deformed by both the glide of dislocations in the basal plane and the climb of dislocations out of the basal plane, following Reference Duval,, Ashby, and Andermn,Duval and others (1983) and Reference Weertman,Weertman (1983). We then calculated the ratio of the climb force, F e, to the glide force. F g, on the dislocations within the kink band, by knowing the lattice rotation from the optical analysis given elsewhere (Reference Manley, and Schulson,Manley and Schulson, 1997) and by making the transformation given in the Appendix. We found that n increased more or less linearly with the ratio F e/F g (Fig. 2a). We made a similar calculation for S2 fresh-water ice (where the c axes are randomly oriented within a plane more or less perpendicular to the long axis of the columnar grains (Reference Michel, and Ramseier,Michel and Ramseier. 1971: Reference Michel, and Ramseier,Michel. 1978) deformed at –10° C by compressing biaxially across the long axis of the columns (Reference Schulson, and Buck,Schulson and Buck, 1995), as shown in the Appendix (Fig. 5). The results fit the same curve (Fig. 2a). Finally, we made the analysis for S2 saline ice of 4–5 ppt salinity, also loaded biaxailly across the columns at –10° C (Reference Kuchn,Kuchn, 1993; Reference Schulson, and Nickolayev,Schulson and Nickolayev, 1995). Again, n increased more or less linearly with increasing F e/F g (Fig. 2b).

Fig. 2. Plots of n values at –10° C vs F e/F g for (a) fresh-water columnar ice and (b) saline columnar ice. Although five points are shown the present work in (Fig. 1, only three points are shown here, because we determined the kink reorientation for the S1 ice in only three cases (Reference Manley, and Schulson,Manely and Schulson, 1997).

Thus, it appears that to a first approximation n=n o+kF c/F g where n o is about 2 and 2.8 for fresh-water and salt-water ice, respectively, and к is around 0.2 for both materials at –10° C.

We suggest that this relationship reflects a transition from essentially two-dimensional (glide) to three-dimensional (glide plus climb) dislocation motion as n increases from 2 to 3. The higher n values may reflect a non-linear dependence of dislocation velocity on stress. It is interesting to note, however, that our criterion predicts that n = 2.2 for uniaxially loaded single crystals whose basal planes are inclined 45° to the direction of loading. This is in agreement with the value n = 2.1 obtained by Reference Kuchn, and Glen,Homer and Glen (1978) who crept such crystals at –10°C.

Hopefully, this note will stimulate further discussion.

Acknowledgements

We would like to thank E. Gratz and J. Mellon for their help with laboratory equipment. Also, we acknowledge the helpful comments by the anonymous reviewers. The work was supported by the U.S. Army Research Office, grant No. DAAL03-92-G-0187.

Appendix Calculations of F c/F g

Consider an edge dislocation within a volume element (Fig. 3). The climbing force is related to the normal stress by F c=bσn, where b is the Burger’s vector. The glide force is related to the shear stress by F g=bσs. So

(A1)

Calculation for S1 ice loaded uniaxially

Assuming a uniform applied stress (Fig. 4), the stresses on a volume element containing a basal edge dislocation on the

plane are

and

where θ is the tilting angle of the basal planes within the kink band (Reference Manley, and Schulson,Manley and Schulson, 1997). Using Equation (Al), the ratio is

The ratio decreases as θ increases from near zero to π/4. The ratio approaches infinity θ = 0. At this point basal glide cannot occur.

Fig. 3. Sketch of glide and climb forces on an edge dislocation.

Fig. 4. Sketch showing reorientation of basal slip plane within kink.

Calculation for S2 ice loaded biaxially

For this case the basal planes (shown as lines crossing the columns; Fig. 5) of individual columns are oriented randomly with respect to the principal stresses (the columns are pointing out of the page in the figure). For an individual column having its basal plane oriented at some angle α with respect to the maximum principal stress axis (σι) the ratio is

For a given angle (>0) the function increases as R 21 increases to near unity. The ratio approaches infinity at =0, π/2 and also when R 21=1. A global ratio for a given specimen was found by assuming the average angle is α = 19.5° . A single crystal placed at this angle with respect to the loading direction would have the same average Schmidt factor as the average factor (0.32) for S2 ice.

Fig. 5. Sketch showing across-column biaxial loading of S2 ice. Traces of the basal planes are shown.

References

Duval,, P., Ashby,, M. F. and Andermn,, I. 1983. Rate-controlling processes in the creep of polycrystalline ice. J. Phys. Chem., 87(21), 40664074.Google Scholar
Frederking,, R. 1977. Plane-strain compressive strength of columnar-grained and granular-snow ice. J. Glaciol., 18(80), 505516.Google Scholar
Homer,, D. R. and Glen,, J. W. 1978. The creep activation energies of ice. J. Glaciol., 21(85), 429444.Google Scholar
Kuchn,, G. A. and Glen,, J. W. 1978. The creep activation energies of ice. J. Glaciol., 21(85), 429444.Google Scholar
Kuchn,, G. A. 1993. Laboratory-grown saline ice — a first-year sea ice analog; the uniaxial compressive properties and the unaixial tensile properties after compressive pre-strain. (M.S. thesis, Dartmouth College, Hanover, NH.)Google Scholar
Manley,, M. E. 1996. Deformation in S1 and S2 ice loaded across the columns. (M.Sc. thesis, Dartmouth College, Hanover, NH.)Google Scholar
Manley,, M. E. and Schulson,, E. M. 1997. Kinks and cracks in S1 ice under across-column compression. Philos. Mag. Lett., 75(2), 8390.Google Scholar
Michel,, B. and Ramseier,, R. O. 1971. Classification of river and lake ice. Can. Geotech. J., 8(1), 3645.Google Scholar
Schulson,, E. M. and Buck,, S. E. 1995. The ductile-to-brittle transition and ductile failure envelopes of orthotropic ice under biaxial compression. Acta Metall. Materialia., 43(10), 36613668.Google Scholar
Schulson,, E. M. and Nickolayev,, O. Y. 1995. Failure of columnar saline ice under biaxial compression: failure envelopes and the brittle-to-ductile transition. J. Geophys. Res., 100(B11), 22,38322,400,400. (Correction: J. Geophys. Res., 101(B3), 1996, 5659.)Google Scholar
Weertman,, J. 1983. Creep deformation of ice. Annu. REv. Earth Planet. Sci., 11, 215240.Google Scholar
Figure 0

Fig. 1. The exponent, n, from the relation vs strain, for S1 fresh-water ice deformed at –10° C. The star (*) was obtained from a six tests (Manley, 1996) from 10−6 to 10−4s−1 where the log of the peak compressive stress was plotted vs the log of the strain rate. All of the other data were obtained from strain-rate change tests.

Figure 1

Fig. 2. Plots of n values at –10° C vs Fe/Fg for (a) fresh-water columnar ice and (b) saline columnar ice. Although five points are shown the present work in (Fig. 1, only three points are shown here, because we determined the kink reorientation for the S1 ice in only three cases (Manely and Schulson, 1997).

Figure 2

Fig. 3. Sketch of glide and climb forces on an edge dislocation.

Figure 3

Fig. 4. Sketch showing reorientation of basal slip plane within kink.

Figure 4

Fig. 5. Sketch showing across-column biaxial loading of S2 ice. Traces of the basal planes are shown.