KnE Materials Science | IV Sino-Russian ASRTU Symposium on Advanced Materials and Processing Technology (ASRTU) | pages: 78-83

1. Introduction

The Northern Sea Route (NSR) undergoes significant changes in the ice conditions every year. Drifting sea ice is characterized by the following parameters: age; thickness; dimensions; concentration; extent of failure, and other. Ice crystals consist of fresh ice, brine, insoluble salts, and air inclusions. Deformation and fracture of ice are due to its structure, temperature, salinity, and type of deformation. Bending failure of ice is typical for ships and hydraulic structures. The work presents the study of changes in the physical and mechanical properties of sea ice in bending.

2. Ice Physics

The annual change in the temperature of the ice thickness T(z) is shown in Fig. 1a. In summer, the ice temperature is practically constant being close to the average value T¯:T(z)T¯ = const. Z coordinate reports with the bottom of the ice surface (Fig. 1a) with a thickness H.

Figure 1

Characteristics of sea ice: (a) temperature variation of sea ice T over the horizon z throughout the year (I...XII – months of the year); (b) changing the Young's modulus E Y (1) and deformation modulus at a bend E B (2) of sea ice by the brine content.


Ice is characterized by salinity: S(z)S¯ = const.

The average salinity S¯ for the seas (NSR) is defined by the following Ryvlin's expression [1]:



Multi-year ice has a salinity S¯=36 ‰. Salinity is also determined by the relative content of the brine in the volume of ice υ p described by formula of Frankenstein and Garner [2]:



Ice is separated into two layers: the upper water-snow has a granular structure, while the lower core layer consists of congelation prisms or fibrillation.

3. The Mechanical Properties of Ice Under Load

Deformation and fracture of ice is described in frame of kinetic and mechanical approaches. The kinetic one describes the development of microscopic damage (defects) basing on studies of Griffith, Irwin, and others. For example, this approach is applied for study of bending and compression of ice by Vershinin [2]. A more common approach uses mechanical stress state diagrams σε (Fig. 2a).

Figure 2

Ice tension diagrams: (a) schematic diagrams σ - ϵ under bending of ice; (b) marine congelation ice in compression (small samples; σ B ,// - loading direction relative to the horizon, the temperature indicated).


Ice deforms elastically at low T¯ and elastoplastically at increasing T¯ (Fig. 2b). In order to take it into account, a simple Hooke's law is replaced by more complex rheological models [3]. At the beginning of the ice deformation, its behavior is described by the elastic modulus (Young's) E Y . With increasing load, secant modulus E e is used (Fig. 2a). The extreme value E i - is deformation module E B = E i m i n when ice is breaking. The relative value of EB/EYσBσY characterizes the degree of ice elasticity. Here, σ Y is the ultimate tensile stress of ice, when it is fully elastically deformed (Fig. 2a). Respectively, 1-EB/EY=σY-σBσY is the measure of ice plasticity.

Figure 3

The cellular model of sea ice and the tensile strength of: (a) ice model of the cellular structure (cell brine); (b) limit σ B bending strength consoles.

Figure 4

Dependence of the relation of deformation module of sea ice at a bend to the limit strength depending on the temperature. Salinity of the ice is 3 ‰.


Estimates made by authors show an elastic behavior of the ice at a bend near the neutral sheet and increase of the ice plasticity near the surface (Fig. 4).

4. The Elastic Characteristics of Ice in Bending

The initial deformation of ice is characterized by the elastic modulus E Y . Reliable way to determine E Y is a seismic dynamic method. In this case, the speed of bending waves is influenced by all features of ice construction [2]. These field experiments are approximated by the following dependence:



Here, the Young's modulus for the freshwater ice EY(0,T¯)9,21(1-0,0015T¯) (GPa).

Simplified example in the study of the bending deformation is the equating of modules (intersecting modules at fracture point) in compression and tension E p = E 1 = E B. The closest case to the interaction of vessels with ice and sloped waterworks are experiments on the destruction of sawn consoles in ice. Under the destruction on a "top-down" the module E B is defined by the following expression [4]:



Dependences of the elastic characteristics of the E Y and E B are shown in Fig. 1b.

5. Results and Discussion

When bending, it is believed that the break-up occurs when the upper limit of tensile strength layer σ 2 is achieved. A direct test of the large ice samples tensile is technologically impossible. Therefore, this test is replaced by the destruction of the console afloat with the definition of conditional flexural strength σ B = σ 2 = σ p under the assumption of the equality of the deformation modules E B = E 1 = E 2. Strength limit σ B is determined using the results of the theory of bending beams [5]:


(5) where ρ - bending radius of neutral console layer.

Presentation of σBS,¯T¯ data is carried out using a cellular model of ice (Fig. 3a) [1]. Cells with brine of relative volume υ p pass through the volume of freshwater ice. Therefore, we obtain the following expression [1]:


(6) where σB(0,T¯) – strength of fresh ice, υ 0 - limiting amount of brine, containing which the ice has no strength (σ B→0); κ - exponent, depending on the shape of the cell section.

Eq. (6) is universal for the entire range of υ p . For a small amount of brine, strength of the cellular structure is more accurately defined by the formula [6]:


(7) where the exponent n is determined from experimental data.

Tensile strength of fresh ice, determined using results of Butyagin [7], Lavrov [8], and others, is (MPa):



The authors come to the following final formulas for the flexural strength:

for γ p <0, 06⋯0, 07


(9) for γ p >0, 07



Results of σB(S¯,T¯) calculations using (8-10) are shown in Fig. 3b.

Defined values of the physical and mechanical characteristics E Y , E B, and σ B of ice can be directly used in the calculation of the ice propulsion of ships and ice interaction with hydraulic structures. Despite the stability of the concept of modulus flexural E B, it will allow to determine reliably the value of the rigidity of the cylindrical field of ice as a plate:


(11) where μ=0,333+6,105·10-2exp0,182T¯- Poisson coefficient of sea ice [2].

The ratio (E B/σ B), in accordance with (5), is related to the radius of curvature of the ice neutral layer in the place of destruction by the following formula:



The nature of the changes of this ratio is shown in Fig. 4.

In the simulation of traffic movements and interactions of hydraulic engineering constructions with ice in the ice model basin, the fulfillment of the condition (E B/σ B) = i d e m is one of the criteria of dynamic similarity of the experiment [9]. The calculated value (E B/σ B) indicates its growth with decreasing temperature T¯ (Fig. 4). This is contrary to the judgment of the opposite tendencies of its measurement [10]. Character of changes of (E B/σ B) indicated in Fig. 4 is defined by the calculation in Table 1.

Table 1

Character of changes (E B/σ B) in bending.

The average temperature of the ice T¯ ,C
-20 -15 -11 -7 -3
υ p at s¯=3 / 8,97 11,43 15,01 22,68 50,78
EB(0,003;T¯),GPa 2,07 2 1,9 1,7 1,15
σB(0,003;T¯),MPa 0,88 0,80(5) 0,76 0,57(5) 0,30(5)
(E B/σ B)*10−3 2,35 2,48 2,5 2,96 3,77
Figure 5

Dependence of the secant modulus (b) and the relative stress (c) sea-ice of congelation at σ¯=-11 C at a bend. (a) Change in temperature with the thickness. s¯ - salinity of the ice.


Numerically, it is explained by the fact that the gradient (negative) dEBdυp is much less than the corresponding gradient dσBdυp . Studies of stressed state of ice in bending show that fracture occurs in the middle of ice stretching zone, rather than on the surface (Fig. 5). This makes the relation (E B/σ B) an artificial characteristic of the physical and mechanical properties of ice. It would be more correct to consider it as the relative value of the curvature radius of the neutral layer in the place of destruction in bending in accordance with (11).

6. Conclusion

Results of research lead to the following conclusions:

1. Strength limit σ B and deformation module EB of ice under bending are conditional mechanical characteristics of sea ice. These characteristics do not reflect the actual stress state of ice destruction at the time.

2. At the same time, module EB allows estimating reliably the cylindrical stiffness of ice field.

3. The ratio of the module to the strength limit EB/σ B determines the relative radius of curvature of neutral layer in the place of ice destruction. It is shown that this ratio increases with the increase of ice temperature.

4. Dependencies for determining of EB and σ B of sea ice that generalize the experimental data are obtained.



Doronin YP., Kheysin DE., Gidrometeoizdat, Year: 1975,


Vershinin SA., Impact of ice on marine structure shelf. Results of science and technology, Series Watercraft , Year: 1988, Volume: 13, MoscowVINITI


Koshkin SV., Generalized model of rheological behavior of ice under load, in. Materials scientific-technical conference2015 Page: 238-240.


Taranukha N. A., Koshkin S. V., Review and research of physical and mechanical properties for sea ice, 24th International Ocean and Polar Engineering Conference, ISOPE 2014 BusanJune 2014kor Page: 1093-1095.


Pisarenko GS., Yakovlev AP, Matveev VV: Handbook of material resistance , Year: 1988, KievNaukova Dumka


Loset S., Shhiner KN., Gudmenstad O., Hoyland K., Effects of ice on marine and coastal structure , Year: 2010, St.PetersburgLan


Butyagin IP., The strength of ice and ice cover , Year: 1966, Novosibirsk, NaukaSiberian Division


Lavrov VV., Deformation and strength of ice , Year: 1969, LeningradGidrometeoizdat


Ionov BP., Gramuzov EM: Ice propulsion of ships. St, Petersburg: Shipbuilding , Year: 2001,


Sasonov KE., Theoretical foundation of navigation in ice. St.Peterburg: Krylov State Reserch Centre , Year: 2010, Krylov State Reserch CentreTheoretical foundation of navigation in ice. St.Peterburg



  • Downloads 32
  • Views 108



ISSN: 2519-1438