Over the past ten years considerable attention has been paid to the development of a continuum model of large-scale sea-ice deformation (see recent reviews by Reference Hibler and ColbeckHibler, 1980, and Reference Coon and PritchardCoon, [c1980]). The most commonly used constitutive law is derived from a yield curve using a “flow rule” (e.g. Reference CoonCoon, 1974; Reference RothrockRothrock, 1975; Reference HiblerHibler, 1977, Reference Hibler1979; Reference Ralston and PritchardRalston [c1980]). The purpose of this note is to point out an alternative derivation for the constitutive law which, while still using a yield curve, avoids using the flow rule and thereby avoids the mathematically imposed relationship between dilatation and pressure-dependent shear strength which the flow rule implies.
Consider the tear-shaped yield of Reference RothrockRothrock (1975) shown in Figure 1. In region A, the yield curves widen indicating an increase in the shear strength with increasing isotropic compressive stress
. This is a realistic feature of the model but, as is well known, the use of the flow rule(with f(σ1, σ2)= 0 on the yield curve) in this region also implies a definite rate of steady divergence. It is certainly true that some materials do dilate initially when they are sheared, but there is no physical law which would require continuing dilatation for all materials whose yield curves diverge.
In region B, the yield curve narrows and cuts across the isotropic stress axis. With the flow rule in mind, this geometry of yield curve is chosen to describe the tendency of the ice to converge under conditions of large compressive isotropic stress. Unfortunately, this geometry carries with it the prediction that the shear strength will weaken with increasing isotropic stress.
It should be clear that to retain full flexibility in choosing a constitutive law, the mathematically convenient flow rule must be avoided. It cannot, after all, be derived rigorously for pressure-dependent materials (Reference HunterHunter, 1976), and its implications are unphysical. In its place we propose a more general procedure starting with the Reiner–Rivlin equation. To illustrate this procedure, consider a two-dimensional continuum with the following properties.
-
(1) No equilibium pressure—the ice has no tendency to expand of its own accord (although it might do so because of random wind or current forcing).
-
(2) No ability to support tension—the ice is presumed to be already cracked in several directions.
-
(3) Resistance to deviatoric shearing is proportional to the isotropic compressive stress holding the ice floes together, and independent of the strain-rate and ice thickness—this is analogous to the Coulomb law for granular materials (Reference CoonCoon, 1974).
-
(4) The ice is nearly non-divergent until a critical value of isotropic stress is reached. Choosing this critical stress proportional to the ice thickness is consistent with the ice-ridging energy analysis of Reference RothrockRothrock (1975). We further assume, for the present purposes, that this process is independent of on-going shearing.
-
(5) The ice is horizontally isotropic.
The Reiner–Rivlin equation is the most general local constitutive relation that satisfies condition (5). In two dimensions it reduces to
where the S i are any of a number of scalar state variables and θ 1, θ 2 are the two invariants of the strain-rate tensor
and
The functions
and are to be chosen to satisfy conditions (1)–(4). This procedure is simplified if we split the stress tensor into its isotropic and deviatoric parts in the usual way.Then
where
and
. The pressure, defined asis
while the deviatoric stress
Condition (4) can be satisfied by choosing the viscosity function γ to be any function which is zero during ice divergence, but which allows the stress to rise quickly to a limiting value (dependent on ice thickness) for ice convergence. An example of such a function is
where C(h) is a thickness-dependent yield stress for the onset of rapid convergence. With the constant n taken to be large, Equation (11) is nearly unaffected by the choice of the normalizing factor
, and the dilatation rate θ 1 becomes vanishingly small unless the isotropic stress is close to C. Note that (1) is satisfied by Equation (11) asThe desired behavior for the deviatoric stress (condition (3)) can be obtained by choosing
The crucial step here is using the pressure p as a state variable in the Reiner–Rivlin constitutive relation. This choice is physically and mathematically correct as the pressure is a scalar invariant of the system. The fact that it makes the constitutive relation transcendental (i.e. the viscosity function now depends on a stress) causes no difficulty, either for statically determinate or indeterminate problems. This is so because the constitutive relation must in any case be solved simultaneously with the other field equations.
For incompressible materials this method is the only way of representing Coulomb behavior within the Reiner–Rivlin system as the rate-of-divergence
contains no information about the level of isotropic stress. For compressible materials other procedures could be used, but this is probably the most convenient way as it explicitly displays the influence of isotropic stress on shear strength. Note also that there is no conflict between the use of pressure as a state variable in μ, and the lack of volumetric elasticity implied by condition (1).An example of Equation (13) which satisfies condition (3) is
If again n is taken to be large, the constant D may be interpreted as a coefficient of sliding or internal friction. The magnitude of D determines the angle between the yield curve and the isotropic stress axis. We require D ≤ 1 in order to satisfy condition (2).
The choice of how large to make n in Equations (11) and (14) depends on how sharp a transition exists between rigid and plastic behavior. For large n, a description in terms of a yield curve (Fig. 2) is appropriate. Note that in this model, the effect of ice thickness on shear strength is indirect. Under compression, the thicker ice can withstand a larger pressure and therefore, through its Coulomb behavior, have an increased shear strength.
The use of this constitutive relation can be illustrated by examining the (statically determinate) problem of a belt of sea ice blown against a straight infinite coastline by a wind stress with on-shore and long-shore components τ ⊥ and τ ∥ respectively. (This problem is equivalent to the problem of a viscous or granular material resting on an inclined plane under the influence of gravity. In both problems, the normal and tangential stress increase proportionately away from the free boundary. The tilt angle of the plane is equivalent to the angle of the wind-stress vector). Neglecting the Coriolis force and the frictional forces beneath the ice, the stress supported by the ice is
where x is measured shoreward of the ice edge (see Fig. 3). The value of σ yy is determined by the following argument. Along a very long coast even a slight long-shore divergence would produce large long-shore velocities which would be resisted by friction. Thus locally, σ yy will take on a value equal to σ xx . A similar result would follow if the shoreline were slightly curved forming a large closed basin. This stress field could be represented on Figure 2 as a straight line emanating from the origin at an angle ϕ from the isotropic stress axis. In the case of small ϕ (nearly on-shore winds) the stress trajectory remains in the nearly rigid field until the normal stress condition for the initial thickness (h 0) is exceeded at position x 1.
From Equations (9), (11), and (15)
Beyond that point the ice thickness increases according to
As the wind stress angle ϕ increases, the stress trajectory eventually coincides with the shearing yield curve ϕ = D. Under these conditions shearing motion begins simultaneously all across the ice belt, in both the thickened and unthickened regions. (This is analogous to the angle-of-repose in the inclined plane problem). This result is very sensitive to choice of yield-curve shape. If the product Dp in Equation (14) were replaced by a more general function g(p), the curvature of this function (and therefore the curvature of the yield curve) would determine whether shearing would be initiated at the shore or at the ice edge.
The use of the flow rule in this problem would lead to a somewhat different result. When the on-shore wind-angle ϕ increased sufficiently to produce shearing motion, the flow rule predicts that the ice would begin to expand outwards away from the coast. This expansion would continue indefinitely, or at least until the ice became so thin that it could not withstand the existing isotropic stresses.
Summary
The object of this note is to suggest a procedure whereby the transition between nearly rigid and nearly plastic sea-ice behavior can be represented mathematically. The method is to use the Reiner–Rivlin equation with pressure included as a state variable and the two viscosity functions chosen to satisfy various empirical constraints. The “flow rule” is not appropriate for this kind of problem as sea ice is a compressible medium. The relationship between shear strength and compressibility predicted by the flow rule is not physically or empirically based, although it may in some circumstances appear reasonable.
The particular viscosity functions in Equations (11) and (14) are derived from conditions (1)–(4), and the solutions to the ice-belt problem are presented for illustrative purposes only. They are not meant to be representative of actual sea-ice dynamics. The author hopes that with a more careful and extensive list of empirical and physical constraints, a useful constitutive relation could be derived.
Acknowledgements
My interest in this problem stems from several stimulating conversations with W. D. Hibler and R. E. Moritz.