Hostname: page-component-cd9895bd7-fscjk Total loading time: 0 Render date: 2024-12-23T11:27:00.624Z Has data issue: false hasContentIssue false

Fracture Criterion for Ice Using a Dislocation Model

Published online by Cambridge University Press:  20 January 2017

V.R. Parameswaran*
Affiliation:
Division of Building Research, National Research Council of Canada, Ottawa, Ontario K1A 0R6, Canada
Rights & Permissions [Opens in a new window]

Abstract

The propensity for crack formation on basal and prismatic planes and twin nucleation on two pyramidal planes in hexagonal ice is examined on the basis of fracture theories and anisotropic elastic theories of dislocations. Calculations show that the ease of crack formation is almost the same in the basal and prismatic planes, whereas it is almost impossible to produce mechanical twinning in ice by stress alone.

Résumé

Résumé

Sur la base des théories de la fracture et des théories élastiques anisotropiques des dislocations on examine les possibilités pour la formation de fissures sur les plans de base et sur les faces pyramidales ainsi que la nucléation des macles sur deux faces pyramidales dans de la glace hexagonale. Les calculs montrent que la formation de fissure est presque aussi facile sur le plan de base et sur les faces prismatiques, alors qu’il est presque impossible de produire des macles dans la glace par le seul exercice d’un effort mécanique.

Zusammenfassung

Zusammenfassung

Auf der Grundlage von Bruchtheorien und Theorien zur anisotropen elastischen Versetzung wird die Neigung zur Rissbildung an Grund- und prismatischen Flächen und zur Zwillingskernbildung an zwei pyramidischen Flächen in hexagonalem Eis untersucht. Berechnungen zeigen, dass sich Risse ebenso leicht in Grundflächen wie in prismatischen Flächen bilden können, dass es aber fast unmöglich ist, mechanische Zwillingsbildung im Eis allein durch Druckspannung zu erzeugen.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1982

Introduction

To formulate the failure criterion of a solid it is important to know the elastic strain energy E of the solid containing cracks and the crack extension force G, which is equal to the negative of the derivative of the elastic energy with respect to crack length. In considering stress fields due to cracks in fracture mechanics, it is customary to distinguish the three cases shown in Figure 1 (Reference Bilby and EshelbyBilby and Eshelby, 1968):

Fig. 1. The three modes of crack propagation in a solid.

The elastic stress field and the energy of such cracks in terms of externally applied stress fields have been derived by several workers: Reference StrohStroh (1958) has formally treated the problem for a general anisotropic medium using the method of Fourier transforms and dual integral equations, and Reference Barnett and AsaroBarnett and Asaro (1972) have used the method of double-ended pile-up of continuously distributed, infinitesimal dislocations, to solve the problem of a slit-like crack in an arbitrarily anisotropic elastic medium under uniform stress. Reference YooYoo (1979, 1980) has used these solutions to calculate the elastic energy of slit cracks in hexagonal crystals and to develop a dislocation model for twinning and fracture in metals of hexagonal close-packed structure.

In this paper the above solutions are used for ice having a hexagonal structure to calculate the ratios of the critical stresses for crack initiation on the basal and prismatic planes, as well as the ratios of the surface energies associated with cracks. The ratios of the critical stresses for basal cleavage and {1012} twin and for prismatic cleavage and {1011} twin are also calculated. From these ratios a prediction is made of the tendency for twinning and cleavage in ice.

Background of the calculations

Reference YooYoo (1979, 1980) derived the equations for critical stresses for crack and twin initiation from the Reference Bilby and EshelbyBilby and Eshelby (1968) theory. By considering a crack as a double-ended pile-up of dislocations, the critical stresses for formation of the three modes of cracks shown in Figure 1 are as follows (Reference YooYoo, 1979): for a Mode I crack under a uniaxial stress σ 22 :

(1)

for a Mode II crack under a shear stress σ 21 :

(2)

and for a Mode III crack under a shear stress σ 23:

(3)

The subscripts 1, 2, and 3 correspond to the cartesian coordinates x, y, and z. The elastic compliance factor Bij is the symmetrical second-rank tensor defined by Reference StrohStroh (1958), depending on the elastic constants Cij and crack orientation; i, j take any of the values 1, 2, and 3; 2a is the length of the crack. The surface energy γ associated with the crack plane x 2 =0, is given by (Reference YooYoo, 1979)

(4)

where is the transformed elastic stiffness constant with respect to the crack plane (hkl) having interplanar spacing d.

For an orthotropic medium (a solid with three mutually perpendicular planes of symmetry) Bij s are the inverse of the energy factors Kij for a dislocation, and expressions for these in terms of standard elastic constants have been obtained by Reference Savin, Savin, Chernov and StrokovaSavin and others (1976). For a dislocation lying on the basal plane with Burgers vector , they show that

(5)
(6)

and

(7)

where

(8)

From Equation (4), the ratio of the surface energies associated with cracks in the basal and prismatic planes is given by (Reference YooYoo, 1979)

(9)

where for {1010} planes and 2 for {1020} planes, c and a are the lattice parameters for the hexagonal lattice, in the present case of ice. Equations (1) to (9) can be used to give the ratio of the critical stress for crack extension in the basal and prismatic planes as (Reference YooYoo, 1979)

(10)
(11)

The subscripts ∥ and ⊥ indicate whether crack extension occurs parallel or perpendicular to the c-axis.

To calculate the critical stress for twin nucleation, Reference Yoo, Haasen, Gerold and KostorzYoo (1980) treats a micro-twin as a slit crack, for simplicity, neglecting the finite thickness of the twin and the welded boundary conditions for twin interfaces. Although far from rigorous, this treatment enables one to estimate the propensity for twin nucleation versus crack propagation. Two cases are considered: (a) basal cleavage versus {1012} twin, and (b) prismatic cleavage versus {1011>} twin.

From Equation (1) the critical stress for crack extension on the basal plane can be written as

(12)

where f c is the total inelastic resistance to crack extension; this includes both the lattice resistance to dislocation motion and the surface energy associated with the crack forces. Considering a micro-twin as a Mode II slit crack, the condition for twin extension from the Bilby–Eshelby theory, is given by Reference Yoo, Haasen, Gerold and KostorzYoo (1980) as

(13)

where denotes the stress components with respect to the twin axis x′i inclined at an angle θ with the original coordinate axis xi , and f t is the total inelastic resistance to twin extension.

For a uniaxial stress σ22,

(14)

For a {1012} ‹1012› twin, which is the most common type in hexagonal close-packed crystals (see Fig. 2a),

(15)

Using Equations (14) and (15) in (13), the critical stress for twin propagation is given by

(16)

where β = c/a. For the isotropic case, B′ 12 = 0 and B′ 11 = B22; hence, from Equations (12) and (16), the ratio of the critical stresses for basal cleavage and {1012} ‹1011› twin is

(17)

Next, considering the ratio of the critical stresses for prismatic cleavage and a {

} twin, for a {
} ‹
› twin (see Fig. 2b):
(18)

Using the conditions of Equations (18) and (14) in (13), the critical stress for twin propagation can be derived as (Reference Yoo, Haasen, Gerold and KostorzYoo, 1980)

(19)

Analogous to Equation (12), the critical normal stress for extension of a Mode I crack on the {

} prism plane is
(20)

From Equations (19) and (20), and assuming B′ 12 = 0 for the isotropic case, the ratio of the critical stresses for prismatic cleavage and a {

} 〈101 2〉 twin is;
(21)

Fig. 2. Geometry of a {1012} ‹ 1011› twin (a), and a{1011} ‹101 2 › twin (b).

Equations (9) to (11), (17), and (21) may now be used to calculate the surface energies of cracks in the basal and prism planes, and the ratios of the critical stresses for crack formation and twinning at various temperatures in ice. For a hexagonal lattice, the non-zero elastic constants are

The temperature dependence of these constants for ice has been given by Reference DantlDantl (1968), and the temperature dependence of the lattice parameters c and a for ice by Reference LonsdaleLonsdale (1958).

Results and discussions

Values of the elastic constants Cij and the lattice parameters c and a calculated for various temperatures from Reference DantlDantl’s (1968) and Reference LonsdaleLonsdale’s (1958) equations are shown in Table I. Figure 3a and b shows the variation of Cij, c, a, and c/a with temperature. Using these values in Equations (9) to (11), the ratios of surface energies (γ b/γ p) and the critical stresses for crack formation (σ b/σ p) in the basal and prismatic planes were calculated. Assuming f c = f t (f c and f t are the lattice resistance to crack and twin propagation, and these include surface energies) in Equations (17) and (21), the ratios of the critical stresses for crack and twin propagation, namely σ c/ σ t, were also calculated. All these values are shown in Table II. The value of γbp{1010} is about 1.03 showing that the probability or ease of crack formation is almost the same in the basal and {1010} prismatic planes. The ratio of the surface energies in the basal and the {1120} planes, however, is about 0.89, indicating that crack formation in the {1120} plane is more difficult than that in the basal or {1010} planes. The value of the ratio σ b/ σ p of critical stresses for crack extension in the basal and prismatic planes is almost equal to 1 for all temperatures, for cracks extending parallel or perpendicular to the c-axis.

Table I. Values of Elastic Constants Cij (units: 103 MPa) and lattice parameters c and a (units: 10–10 m)

Table II. Values of Energy Factors Kij = 1/Bij (units: 103 MPa) and ratios of surface energy and critical stresses for crack and twin propagation

Fig. 3. (a) Variation of elastic constants Cij of ice with temperature, (b). Variation of lattice parameters c and a with temperature.

These results suggest that, under a uniaxial stress, cracks could propagate in ice with equal ease, either along the basal or along the prism planes, although the probability of their occurrence on the basal and the {1010} planes is slightly more than that along the {1120} plane. (In this paper, crack propagation along planes at other angles to the c-axis, such as the pyramidal planes, has not been considered.) This coincides with observations made by Reference GoldGold (1972), who carried out detailed and extensive studies on various aspects of cracking activity in ice. He found that in polycrystalline ice subjected to uniaxial compressive stress at temperatures between –5 and –30°C, about 25% of the observed cracks are along the grain boundaries, 42% in the basal plane, and 20% perpendicular to the basal plane; the rest were irregular. He also observed that the most probable direction of crack propagation is parallel to the compressive stress direction. Reference Parameswaran and JonesParameswaran and Jones (1975) tested single-crystal and polycrystalline ice at – 196°C (77 K) under uniaxial compression. No particular orientation dependence was observed for crack propagation. The modes of failure were in general: (a) axial cleavage with one or more major cracks splitting the sample along the loading direction, and (b) cataclasis or general internal crumbling by formation of multiple cracks. The cracks did not follow any particular pattern or crystallographic plane.

Considering the ratio σ c t of the critical stress for basal cleavage and {1012} twin, or prismatic cleavage and {1011} twin, it may be seen from Table II that σ c/σ t has values of 0.5 and 0.41 respectively for the two cases. This shows that mechanical twinning is almost impossible in ice, because the critical stress required for twin initiation is twice that for cleavage crack propagation. This explains why no one so far has observed mechanical twinning in ice crystals, although several workers have studied the deformation behaviour of ice under various conditions. This does not preclude the existence of growth twins in ice. The peculiar shapes of snow crystals (grown under both natural and artificial conditions) observed by investigators such as Reference Bentley and HumphreysBentley and Humphreys (1931), Reference NakayaNakaya (1954), Reference KikuchiKikuchi (1970), Reference Kikuchi and HoganKikuchi and Hogan (1976), and Reference YamashitaYamashita (1971), have recently been explained as twins by Reference Kobayashi, Kobayashi, Furukawa, Kikuchi and UyedaKobayashi and others (1976) on the basis of the generalized concept of coincidence-site lattice (CSL) theory, Reference Kobayashi and FurukawaKobayashi and Furukawa (1975) have also discussed the formation of the twelve-branched snow crystals observed by Reference DoiDoi (1832, 1840), Reference Bentley and HumphreysBentley and Humphreys (1931), and Reference NakayaNakaya (1954), in terms of rotation twinning.

From Equation (1) a theoretical value of the critical stress for crack extension can also be estimated. The critical stress intensity factor or strain-energy release rate for ice was found to depend on temperature and strain-rate, and under conditions of crack arrest associated with thermal shock Reference GoldGold (1963) calculated the strain-energy release rate to be about 300 ergs cm-2. This means that

or

where γ is again the surface energy and P is the external work done. Reference GoldGold (1972) uses a value for γ w= 100 ergs/cm-2 (10–7 MN m–1). One can also calculate the width 2a of a critical crack by considering it as a pile-up of n freely slipping edge dislocations (Reference BulloughBullough, 1964)

(22)

where μ is the shear modulus (C 66 for ice), and b is the magnitude of the Burgers vector (4.523 × 10–10 m). The width of the critical crack depends on the choice of n, the number of dislocations in a pile-up. Reference BulloughBullough (1964) suggested a value of n = 300–700 and Reference Friedel and VassamilletFriedel (1964) suggested that n could be 102 to 103. Using the values of μ (=C 66) and v from Table I, and assuming γ= 100 ergs cm–2 (1 – 10–7 MN m–1), 2a=0.23 mm at –196°C and 0.2 mm at –10°C for n= 500, and 2a=0.92 mm at –196°C and 0.8 mm at –10°C for n=103. These values are in good agreement with the experimentally observed value of 1 mm for the first large crack in ice subjected to a compressive stress (Reference GoldGold, 1972).

Using these values of a in Equation (1), the critical stress values calculated for crack formation at –10°C were 2.49 and 1.25 MPa, respectively. These agree very well with the observed failure stresses for ice measured under tension, namely, 1 to 3 MPa (Reference MichelMichel, 1978). The critical stress values calculated from Equation (1) for ice at –196°C, using n = 500 and 1000, were 2.52 MPa and 1.26 MPa, respectively. These values are about one order of magnitude smaller than the fracture stresses observed by Reference Parameswaran and JonesParameswaran and Jones (1975) Footnote * . The values of failure stress under compression measured at –196°C were in the range of 24 to 39 MPa for natural as well as artificially grown single and polycrystalline ice. They observed, however, that the first crack appeared in the very early stages of deformation in the stress range 0.1–0.5 MPa.

Conclusion

A dislocation model based on theories of fracture and anisotropic elasticity has been used to calculate the ratio of surface energies and the critical stresses for formation of cracks in the basal and prismatic planes of hexagonal ice. These calculations show that crack formation occurs with almost equal ease in both basal and prismatic planes. The calculations also show that it is almost impossible to produce mechanical twinning in ice by stress alone, because the critical stress for twin initiation is twice that required for cleavage crack formation, and this explains why no one so far has observed mechanical twins in ice.

Acknowledgement

Sincere thanks are extended to Dr L. W. Gold for helpful discussions. This paper is a contribution from the Division of Building Research, National Research Council of Canada, and is published with the approval of the Director of the Division.

Footnotes

page note page 168 note * The numbers along the ordinates in figures 4, 5, and 9 of this reference are incorrect. The correct values of fracture stress in figures 4 and 5 vary between 10 and 100 MN m–2 (or MPa). not 100 to 1 000 MN m–1. In figure 9 the fracture stresses are in units of MN m–1 (or MPa), not 104 MN m–2.)

References

Barnett, D. M. and Asaro, R. J. 1972. The fracture mechanics of slit-like cracks in anisotropic elastic media. Journal of Mechanics and Physics of solids, Vol. 20, No. 6, p. 353–66.Google Scholar
Bentley, W.A. and others Humphreys, W. J. 1931. Snow crystals. New York, McGraw-Hill Book Co., Inc. Google Scholar
Bilby, B. A. and others Eshelby, J. D. 1968. Dislocations and the theory of fracture, (In Liebowitz, H., ed. Fracture: anadvanced treatise. New York. Academic Press, Vol. 1, p. 99182.)Google Scholar
Bullough, R. 1964. The cracked dislocation under tension. Philosophical Magazine, Eighth Ser., Vol. 9 No. 102 p. 915–25.Google Scholar
Dantl, G. 1968. Die elastischen Moduln von Eis-Einkrislallen Physik der kondensierten Materie. Bd. 7, Ht. 5 p. 390–97..Google Scholar
Doi, S. 1832. Sekka zusetsu [Illustrations of snow crystals], 2 vols. [Reprinted, with Doi (1840). in the series Nihon Kagakukoten Zensho (Tokyo, Asahi Shinbun-sha). [Vol.] 6, No. 19, p. 629–84..]Google Scholar
Doi, T. 1840.>. sekka zusetsu [Supplementary volume of illustrations of snow crystals], 2 vols.Google Scholar
Friedel, L. W. 1964.>. Dislocations. English edition, based on a translation by Vassamillet, L.F. Oxford,. Pergamon Press. .+Dislocations.+English+edition,+based+on+a+translation+by+Vassamillet,+L.F.+Oxford,.+Pergamon+Press.>Google Scholar
Gold, L. W. 1963. Crack formation in ice plates by thermal shock. Canadian Journal of Physics Vol. 41, No. 10, p. 1712–28.CrossRefGoogle Scholar
Gold, L. W. 1972.>. The failure process in columnar-grained ice. Canada. National Research Council. Division of Building Research. Technical Paper No. 369 Google Scholar
Kikuchi, k. 1970. Peculiar shapes of solid precipitation observed at Syowa Station, Antarctica. Journal of the Meteorological Society of Japan, 48 No. 3, p. 243–50.Google Scholar
Kikuchi, K. and Hogan, A. W. 1976. Snow crystal observations in summer season at Amundsen–Scott South Pole Station, Antarctica. Journal of the Faculty of Science, Hokkaido University. Ser. 7, vol. 5 No. 1 p. 120.Google Scholar
Kobayashi, T. and Furukawa, Y. 1975. On twelve-branched snow crystals. Journal of Crystal Growth.vol. 28, No. 1 p. 2128.Google Scholar
Kobayashi, T and others 1976.>. On twinned structures in snow crystals. [by] Kobayashi, T and Furukawa, Y and Kikuchi, and Uyeda, H. Journal of Crystal Growth. vol. 32, No. 2, p. 233–49.CrossRefGoogle Scholar
Lonsdale, K. 1958. The structure of ice. Proceedings of the Royal Society of London, ser. A, Vol. 247, no. 1251, p. 424–34.Google Scholar
Michel, B. 1978.>. The strength or polycrystalline ice. Canadian Journal of Civil Engineering, Vol. 5, no. 3 p. 285300.Google Scholar
Nakaya, U. 1954.>. Snow crystals. Cambridge, Mass., Harvard University Press. Google Scholar
Parameswaran, V., and Jones, S. J. 1975. Brittle fracture of ice at 77 k. Journal of Glaciology, vol. 14 No. 71 p. 305–15.Google Scholar
Savin, M. M. and others 1976. Energy factor of dislocations in hexagonal crystals, by Savin, M. M. Chernov, V. M. and Strokova, A. M. Physica Status Solidi(a). vol. 35, No. 2 p. 747–54.Google Scholar
Stroh, A. 1958 Dislocations and cracks in anisotropic elasticity. Philosophical Magazine Eighth Ser. vol. 3 No. 30 p. 625–46.Google Scholar
Yamashita, A. N. 1971. Skeleton ice crystals of non-hexagonal shape grown in free fall. Journal of the Meteorological Society of Japan, vol. 49 No. 309 p. 215–31.Google Scholar
Yoo, M. H. 1979. The elastic energy of slit cracks in hexagonal crystals. Scripta Metallurgica, vol. 13 No. 2 p. 131–36.CrossRefGoogle Scholar
Yoo, M. H. 1980. A dislocation model for twinning and fracture and its application to H.C.P metals. (In Haasen, P., and others, ed. Strength of metals and alloys. Proceedings of the 5th international conference, Aachen, Federal Republic of Germany, August 27–31,1979. Edited by Haasen, P. Gerold, V. and Kostorz, G. Oxford etc., Pergamon Press, Vol. 2, p. 825–30.)Google Scholar
Figure 0

Fig. 1. The three modes of crack propagation in a solid.

Figure 1

Fig. 2. Geometry of a {1012} ‹1011› twin (a), and a{1011} ‹1012› twin (b).

Figure 2

Table I. Values of Elastic Constants Cij (units: 103 MPa) and lattice parameters c and a (units: 10–10 m)

Figure 3

Table II. Values of Energy Factors Kij = 1/Bij (units: 103 MPa) and ratios of surface energy and critical stresses for crack and twin propagation

Figure 4

Fig. 3. (a) Variation of elastic constants Cij of ice with temperature, (b). Variation of lattice parameters c and a with temperature.