Introduction
Many workers have studied the uniaxial mechanical behaviour of ice but few have bothered with more complicated stress systems, even though most engineering, and glacier-flow, situations involve complex stress fields. Reference Parameswaran and JonesRigsby (1958) studied the effect of hydrostatic pressure on the shear creep behavior of four single crystals and concluded that if the temperature was measured from the pressure-melting point, hydrostatic pressure had no effect on the creep strain-rate. In other words, the creep strain-rate was increased by the hydrostatic pressure, at constant temperature. Reference Haefeli, Haefeli, Jaccard and de QuervainHaefeli and others (1968) deformed three polycrystalline ice samples in compressive creep with a 30 MN m–2 superimposed hydrostatic pressure and they also showed that the creep-rate was increased by the hydrostatic pressure. At higher strain-rates, comparable to those used in the present study, Reference Goughnour and AnderslandGoughnour and Andersland (1968) noticed a slight increase in the strength of pure ice samples with a confining pressure of 0.7 MN m–2. Reference Simonson, Simonson, Jones and GreenSimonson and others ( [1975]) conducted Triaxial tests on both ice and a sand-ice mixture up to 200 MN m–2 at –10°C and found a “monotonic decrease in strength to the highest pressure tested”. There have been several studies on various sand-ice and permafrost systems (Reference Chamberlain, Chamberlain, Groves and PerhamChamberlain and others, 1972; Reference Alkire and AnderslandAlkire and Andersland, 1973; Reference Smith and CheathamSmith and Cheatham, 1975; Reference ParameswaranParameswaran, 1980; Reference Parameswaran and JonesParameswaran and Jones, 1981). The object of this work was to study the effect of confining pressure on the compressive behaviour of ice over a wide range of strain-rates, from the normal creep range to the brittle fracture range. This paper reports all the results obtained at –11 ± 1°C, in the strain-rate range 10–7 to 10–1 s–1 Preliminary data were given in Reference JonesJones (1978).
Sample preparation
Details of the method of sample preparation were given in Reference JonesJones (1978). Random polycrystalline samples with grain-size less than 1.0 mm diameter were produced; these had a density of 0.913 ±0.001 Mg m–3 at –10°C which indicates the presence of some air bubbles. These bubbles were too small to be seen individually with the naked eye but gave the ice a slightly hazy appearance. Figure 1 shows a typical thin section through crossed polaroids. The average grain diameter of this sample is 0.8 mm as determined by counting the number of intercepts of grain boundaries with a single straight line. The sample dimensions for the mechanical tests were 20 mm diameter by 60 mm long. Stainless-steel caps were frozen onto both ends of the ice, using a special jig to maintain alignment of the sample and end caps. These end caps prevented any sideways motion of the ice sample while it was in the testing machine.
Method
Details of the high-pressure cell and general method of testing are given in Reference JonesJones (1978). In brief, the ice was deformed at a constant strain-rate (strictly constant rate of deformation) in compression on an Instron model 1116 tester with a superimposed hydrostatic pressure of up to 85 MN m–2. For comparison, unconfined tests were carried out at atmospheric pressure (0.1 MN m–2). For strain-rates less than 5 ×10–3 s–1 the load–time curve was recorded on the Instron chart recorder. Above this strain-rate the load was recorded on a galvanometric type recorder (Kyowa Electronic Instrument Co. Rapicorder RMV-300E) which had a much faster response time than the standard Instron chart recorder. Reference JonesJones (1978) explains the precautions that were necessary, and the method by which the load-time curves were converted to true-stress-true-strain curves.
Results and discussion
Figure 2 shows some of the stress-strain data obtained at a strain-rate of 5.5 ×10–3 s–1. The stress plotted in Figure 2 is the differential stress i.e. σ 1 – σ 3, the amount by which the axial stress departs from the confining pressure. At this strain-rate the unconfined test gave an elastic deformation followed by brittle failure, as shown in the Figure. The application of a 10 MN m–2 confining pressure suppressed this brittle failure and the associated cracking, and allowed the stress to rise much higher than in the unconfined test, leading to plastic yielding and flow. The yield stress, or strength, thus defined, was greater for confining pressures of 34 MN m–2, but for higher confining pressures the strength value fell, as shown by the curve for an 85 MN m–2 confining pressure in Figure 2. If these yield stresses, or brittle failure stresses, as the case may be, are then plotted against confining pressure, curves such as those shown in Figure 3 are obtained for different strain-rates. As can be seen by comparing these curves, the effect of the confining pressure at low strain-rates is much less than that at high strain-rates. This is because, at low strain-rates, the ice deforms plastically in the unconfined tests, as well as in the confined ones, with little or no cracking. At high strain-rates the confining pressure is needed to suppress the cracking which would otherwise lead to failure.
Figure 3 also shows that at the high strain-rates, the rise in strength with confining pressure is rapid, with the bulk of the increase occurring in the first 10 MN m–2 confining pressure, whereas the decrease in strength after the maximum is reached is relatively gradual. At the low strain-rates (≈10–6 s–1) the maximum is reached at a confining pressure of about 10 – 20 MN m–2 whereas at the high strain-rates (≈10–3 s–1) the peak occurs at about 40 MN m–2. However, there is too much scatter in the data to determine the exact dependence of peak strength on confining pressure. For a temperature of –11°C, the melting point would be reached at a confining pressure of about 100 MN m–2, above which the strength would have to be zero.
This decrease in strength with increasing confining pressure, after the maximum has been reached, is probably due to an increase in the amount of liquid at the grain-boundary junctions. Reference Barnes, Barnes, Tabor and WalkerBarnes and others (1971) have shown that above –10°C grain-boundary sliding and the presence of liquid at the triple junctions of the grain boundaries serve primarily as accommodation processes for the high stresses generated by basal glide within the grains. In the present experiment the confining pressure causes a lower melting point so that more triple-junction melting occurs, and hence a lower stress is required to maintain the deformation rate. If this is so, we would not expect to see the same softening if the tests were conducted at temperatures of –20°C or lower at the same confining pressures, because triple-junction melting is not thought to be significant below –10°C.
The peak strength obtained at each strain-rate, from curves such as those shown in Figure 3, was plotted as a function of strain-rate, as shown in Figure 4. The unconfined, atmospheric-pressure results are also shown for comparison. The confined strengths are always greater than the unconfined, but the two sets of data converge clearly at low strain rates. At high strain-rates, above 10–2 s–1, the unconfined strengths decrease rapidly, whereas the confined strengths continue to increase, although the rate of increase with confining pressure is less. However, this latter result may be an artefact of the experiment because at high speeds (greater than 1.4 ×10–2 s–1) the load capacity of the Instron testing machine was limited to 5 000 kg, which in turn limited the confining hydrostatic pressure to 50 MN m–2, and it is possible that a larger confining pressure would have given a greater strength than that shown in Figure 4. However, there is certainly no decrease in strength in the confined tests above 10–2 s–1.
Two straight lines and one curve, have been drawn on Figure 4. The straight lines are the best fit to the unconfined and confined data for strain-rates between 10–6 and 10–3 s–1. The data for lower strain-rates are less accurate because only a few tests were carried out, and at higher strain-rates the data clearly deviate from a straight line, as discussed above. A straight-line relationship on the double logarithmic scale of Figure 4 implies ∝σ n where n is the slope of this line. For the confined tests n = 3.95 and for the unconfined tests n = 5.04. This value of 5.04 for the unconfined tests is in agreement with other workers at high strain-rates (Reference RamseierRamseier, 1976), who has analyzed his results in a similar way. The curved line drawn on Figure 4 was generated using the equation of Reference Barnes, Barnes, Tabor and WalkerBarnes and others (1971):
and using mean values of A′(2.58 ×1010 s–1), α(0.266 MN–1 m2), n(3.02), and Q(78.1 kJ mol–1) derived from their creep results between –8°C and –45°C. The curve is, therefore, not a best fit to my data but, rather, it shows that the present unconfined constant strain-rate data agree very well with the Barnes and others constant-stress results. This supports the idea of a correspondence principle outlined by Reference Hooke and HookeHooke and others (1980). It also shows that it is largely a matter of personal preference whether one uses a power law or a hyperbolic sine law to express empirically the dependence of stress on strain-rate, at least over this range of strain-rates.
At low strain-rates below about 10–6 s–1, where there is no cracking evident, most workers have found n to be nearer to 3. My value for the confined tests 3.95 falls half-way between this and my value of 5.04 for the unconfined tests. Since the two straight lines converge at low strain-rates -and even cross if, in creep, hydrostatic pressure increases the flow rate slightly (Reference Parameswaran and JonesRigsby, 1958; Reference Haefeli, Haefeli, Jaccard and de QuervainHaefeli and others, 1968) – it is not possible for the n value of the unconfined tests to be less than the n-value for the confined tests. Therefore, at such low strain-rates, the w-value for unconfined tests will not fall below 3.95, unless the confined tests also give an n less than 3.95 in this strain-rate range.
Although a straight line has been drawn through the unconfined data of Figure 4, it can be seen that between the strain rates of 10–5 and 10–4 s–1 there appears to be a “plateau” in which the stress increases very little. Many tests were carried out specifically in this strain-rate range, as shown by the numbers on Figure 4, and this “plateau” seems genuine. The Instron machine speed was checked and was correct. No other workers have commented on this “plateau” and the reason for it is not clear. It does not mark the transition from no-cracking to cracking in the ice, since this occurs at a lower strain-rate (between 10–7 and 10–6 s–1). There does not appear to be a similar plateau in the confined tests.
Above 5 ×10–3 s–1, the unconfined tests show a dramatic reduction in strength with further increase in strain-rate, whereas the confined strength continues to increase. This reduction in unconfined strength is due to a combination of two effects. First, end effects at high strain rates are important (Reference Hawkes and MellorHawkes and Mellor, 1970) and slight irregularities in the flatness of the ends or in the alignment of the end caps can cause stress concentrations and lead to premature cracking and failure. However, because ice is a relatively soft material, this may not be so important as the second factor, which is that the Instron testing machine itself is not infinitely rigid. This leads to the release of strain energy stored in the machine when the test specimen reaches the limit of its resistance. The loading platens accelerate violently and the sample is destroyed. When there is a confining pressure present, it acts as a stiffening element on the testing machine and this may be the reason why the confining pressure tests in Figure 4 do not show a reduction in strength at the high strain rates (Reference Sinha and FrederkingSinha and Frederking, 1979). Also, other workers, Reference MichelMichel (1978), have observed the reduction in unconfined compressive strength beginning at a strain-rate of about 5 ×10–4 s rather than the value of 10–2 s–1 found here. This is, presumably, because their testing machine was less rigid, and their end effects more important, than in this study. It is reasonable to expect that another study, using a more rigid testing machine, would not show the reduction until even higher stresses were reached.
As mentioned earlier, the samples contained a small amount of air, about 0.6% by volume. At –11°C and confining pressures above 10.5 MN m–2 this will form a clathrate structure within the ice. Such clathrates are thought to occur in the deep Antarctic ice (Reference MillerMiller, 1969). If the mechanical properties of this air clathrate had been very different from those of ice, we might have seen a significant change in yield stress at a confining pressure of about 10 MN m–2. However, no such change was noticed in Figure 3 and Figure 4 from which we conclude that a small amount of air clathrate does not significantly affect the flow of ice. We intend to carry out a further study on the effect of air bubbles on the strength of ice.