Loading...
 
[Show/Hide Left Column]
 

A Review of Ice Loading and the Evolution of the Finnish-Swedish Ice Class Rules

By Kaj Riska, DSc., Director, Helsinki, ILS Oy, Professor, Norwegian University for Science and Technology

By Jorma Kämäräinen, DSc., Senior Maritime Inspector, Finnish Transport Safety Agency

Shorten version. Full article is available as PDF attachment at the end of this page

DRAFT TO THE SNAME ANNUAL MEETING IN 2011

ABSTRACT

The Finnish-Swedish Ice Class Rules (FSICR) have their origin in rules given in 1890. Since that time, the rules have evolved and, at present, the FSICR can be considered an „industry standard‟ for designing ships for first-year ice environments. The evolution of the rules has mostly been based on feedback from ships, in the form of ice damage, though new research results have also influenced rule development. In this article, the background to the hull rules in the FSICR is described and analyzed. First, it makes it clear that the rules are an integral part of the winter navigation system for Finland and Sweden. The rule formulations are described and great attention is given to analyzing the design point in the rules – the design point includes the limit state used and the frequency of reaching it. In order to analyze the design point, statistics of ice loads are investigated, based mainly on the measurements carried out in the Baltic. The analysis of the design point in the FSICR suggests that it entails reaching the yield limit once in a winter or the plastic limit once in a ship‟s lifetime. The strength of the plating and frames is also quite balanced.

Keywords: Icebreaker, ice load, ice rules, ice load measurements, plastic response

1. INTRODUCTION

Navigation in ice-covered waters used to be a seasonal adventure during the sailing ship era. With sail powered ships, ice could not be broken, only avoided. With the advent of steel hulls, power by steam, and, later, diesel engines and propulsion by propellers, ice could be forced to a certain extent. The first real winter navigation started in the Great Lakes in the mid-19th century. It quickly became evident that more economical and safe navigation would be achieved if special purpose icebreakers were used to assist the less ice-capable merchant ships. The first steam-powered icebreaking ship was the City Ice Boat No. 1, built by the city of Philadelphia in 1837. She was a wooden paddle steamer intended to break ice in the harbor. The first European steam-powered icebreaker, as well as the first ever metal-hull icebreaker, was the Russian Pilot, built in 1864. Its bow had been designed with some ice-breaking capability (the stem angle was 20°). These and other early icebreakers were able to extend the summer season, but it took several decades, to the end of the 19th century, before navigation through ice – winter navigation – started to have economical importance.

When winter navigation commenced on an economically significant scale, it became necessary to regulate the construction and use of ships intended for ice. The first rules for winter navigation were issued in Finland in 1890 – these were called the Imperial Statutes for construction and fitting out ships for winter navigation (see, e.g., Siivonen 1977), as Finland was part of Russia at that time. Since that time, many updates to the Finnish ice class rules have been made: updates that have been necessary in view of the feedback from the earlier rules in the form of data about ice damage, or when new technical solutions have been introduced. The Finnish ice class rules have and still are an example for classification societies and, at present, these rules can be considered an industry standard for ships navigating in sea areas in which there is only first-year ice.

Each set of Finnish ice class rules (since 1932, there have been more than one strengthening level, ice class, in the ice class rules) is based on some idea about ice loading as well as some acceptable strength level. As knowledge about ice loads improves and feedback from ice damage accumulates, changes to the rule foundation warrant rule updates. In this paper, the knowledge background that forms the foundation for the present Finnish-Swedish Ice Class Rules (FSICR) is described. The aim is to make the rule background transparent and indicate how the new knowledge could be incorporated into the rule updates. The FSICR include requirements for hull and machinery scantlings as well as for the minimum propulsion power of ships, but the present article is restricted to describing only the hull rules.

The aim of the article is to describe the design point of the hull structures in the FSICR. The design point includes a description of the loading, especially frequency and magnitude, and allowable response. As the FSICR are an integrated part of Baltic winter navigation, the background to the FSICR is described first. After this, the rules themselves are described in short, followed by a lengthy description of the factors included in determining the rule ice loading. The structural response to be included in the design point is described before the design point that is finally selected (and some alternatives to it) is touched upon. The article tries to show that putting together ice class rules is a holistic process in which many different parts of the rules interact. Thus, it is not possible just to change one part of the rules that is deemed incorrect without a complete revision or at least a review.

2. BACKGROUND TO THE FINNISH-SWEDISH ICE CLASS RULES

Finland is one of two countries in the world whose ports are all ice bound every (average) winter – Estonia being the other. Fluent, efficient, and economical winter navigation is thus a prerequisite of all economic activity in Finland. Even Sweden is less dependent on winter navigation, as ports south of Stockholm only experience minor hindrance from ice.

2.1 Winter Navigation System

The Finnish and Swedish winter navigation system is based on: a) icebreakers that escort the merchant fleet through the worst ice conditions; b) an ice-strengthened merchant fleet; and c) actions of maritime authorities (rules, regulations, and fees). As the Baltic has seasonal ice cover, the ice season lasts from three to five months annually, and all Finnish ports are icebound every average winter, ships visiting Finland and northern Sweden year-round must be ice capable. The merchant ships sailing in the northern Baltic compete with open water ships during a large part of the year, thus the ice strengthening and ice performance must not have a great effect on their competitiveness, vis-à-vis the open water ships. This has led to the present balance between the required ice capability of merchant ships and the number of icebreakers escorting them, and the safety level of winter navigation. The interaction between the different parts of winter navigation makes it a system.

2.2 Elements of the Winter Navigation System

The maritime authorities set the FSICR and the requirements for ice classes (traffic restrictions) as well as the fairway dues. The Finnish-Swedish Ice Class Rules include requirements for the ship hull and machinery strength as well as some requirement for ship performance in ice. The requirement for ice performance interacts with the icebreaker escort service: if the merchant ships had lower ice performance, they would need more icebreaker escorts. More icebreakers would then be needed to serve all the ships – as all ships fulfilling the traffic restrictions are given icebreaker escort, if needed – and this would lead to increased costs. The increased costs would be covered by higher fairway dues. The Finnish fairway dues are scaled with the ice class: higher ice class ships pay lower fairway dues.

In order to reduce the escort times of merchant ships and to make winter navigation as efficient as possible, requirements are set for ice classes and for a minimum amount of cargo. These restrictions are set for each port and are increased as the winter progresses, following the development of the ice cover. The rationale of the minimum amount of cargo is that if each ship only carries a small amount of cargo, more ships are needed for a certain amount of trade and thus more icebreaker escorts are required per transported ton. The traffic restrictions are declared roughly one week before they come into force, and they are published on the home page of the Finnish Transport Agency (www.liikennevirasto.fi). The traffic restrictions are also shown in the daily ice charts issued by the Finnish (and at the same time Swedish) Ice Service; see Fig. 3. The ice charts contain information about the ice cover, traffic restrictions, and icebreakers active in different sea areas.

The fairway dues are scaled according to the ice class, and the amount of the fairway due is set by the net tonnage (NT) of each ship. There is a maximum amount of fairway due (at present 98,400 €). The fairway due is paid for each port visit, though cargo ships do not need to pay the fairway due more than 10 times per year (30 times per year for passenger ships).

2.3 Evolution of the FSICR

The first Finnish rules for winter navigation were given in 1890. These rules only contained requirements for ship equipment (e.g., how many pairs of skates and skis had to be on board) and the general and structural arrangements on the ship. The first Finnish ice class rules for shipping were given in 1920. These rules started the long tradition of the so-called „percentage rules.‟ In these rules, the scantlings were set as some relative increase in the open water scantlings. If the relative increase in the plate thickness and the frame section modulus are the same, the percentage practice leads to a greater increase in plate strength, as plate thickness is proportional to the square root of the strength, whereas the frame section modulus is linearly dependent on strength. This works against a sound structural strength hierarchy. The first Finnish ice class rules were given in 1932, when three ice classes (IA, IB, and IC) were introduced, as well as ice class II corresponding to open water ships and ice class III corresponding to barges. At this time, the ice classes were also coupled with the fairway dues. Only minor corrections to these rules have been made in four decades, and the last change was the introduction of ice class IA Super in 1965.

The year-round trade to northern Finnish ports (Kemi, Oulu, and Raahe, as well as the Swedish port of Luleå) increased greatly in the 1960s. The increased trade caused three major changes to the Finnish winter navigation system. The first one was that ships built according to the percentage rules were found to be too weak – in the sense that the amount of ice damage exceeded the economic pain threshold. This led to a large ice damage survey and complete revision of the ice class rules; see Johansson (1967). For the first time, the damage survey gave an idea of the magnitude of the ice pressure, with the cost of fixing the load height – this was chosen to be the same as the maximum level of ice thickness, i.e., 800 mm. The rules themselves, based on the ice damage survey, were published in 1971.  

The second major change to the Finnish winter navigation system was a result of an agreement between Finland and Sweden to develop the winter navigation system together. Icebreaker services in the Bothnian Bay and Sea would also be organized together. As the ice class rules are an essential part of the winter navigation system, the rules were to be given jointly; thus the 1971 set of rules is called the Finnish-Swedish Ice Class Rules. The third major change to the winter navigation system was the realization that far more powerful icebreakers were needed. The result was the development of Urho class icebreakers, which have almost twice the power of the (at that time) previous icebreakers. The Urho class, delivered between 1974 and 1977, can be considered the last classic icebreaker class with two bow and two stern propellers.

The next revision of the FSICR occurred in 1985 when the hull rules were changed. The main change was to the ice load height, based on extensive ice load measurements on merchant ships and icebreakers; see for example, Vuorio et al. (1979), or Kujala and Vuorio (1986). The measurements showed that the load height compatible with the measured ice load and structural response is much less than the ice thickness. The ice load height was made the ice class factor instead of ice pressure. The idea here was that ice is similar, even if the thickness is somewhat greater, and thus the ice pressure, being an ice material constant, is the same for each ice class, but as the ice thickness related to each ice class is different, the load heights are consequently different for each ice class. While the load height changed, the line load determined from the ice damage data remained constant in the rules.

Recent measurements indicate that the ice edge is failing in a fashion to create a very narrow contact between the structure and the ice edge. This produces a very small contact height; see, for example, the first observation of this type of contact in Riska et al. (1990). The problem, for design purposes, of this line-like contact is that while the local pressures are high (at least 60 MPa has been measured in a very small area; see Frederking and Sudom 2008), the measured ice load obtained from controlled tests divided by the load length along the contact line, q, is relatively small compared with the measured values from ship tests. It has also been observed that when ice is relatively warm (close to melting point) or the relative velocity between the ice edge and structure is low (some cm/s), the ice does not fail in a manner that produces the line but rather so that the contact is on the whole height of the ice sheet. At the same time, the line load q is larger than in the brittle case, in which the line forms; see Sodhi et al. (1998). The design case is not clear, but the suggestions for future research given toward the end of the paper have now already been anticipated.

The 1932 ice class rules already contained a requirement for performance in ice – this was given as the required propulsion power (PD) as a function of the product of ship length and breadth (this product was changed to displacement in the 1971 rules), with different constants for each ice class. In the 1985 rules, this power requirement was amended so that the constants in the previous equation included some dependency on ship beam and bow shape. The ice performance requirement was changed drastically in 1999 (and slightly amended in 2002) when an ice performance requirement (instead of power requirement) was introduced. The requirement states that the ship must be able to maintain a speed of at least 5 knots in a brash ice channel the thickness of which is a class factor (1 m for the ice class IA). If the designer cannot show compliance with this, the rules give an equation for the rule channel resistance (RCH) and propulsion power. The 1999 rules give quite high propulsion power for large ships, see Fig. 6, and thus the rules contain a proviso for this kind of case. In fact, the rules state that the designer can show compliance with the ice performance requirement by, e.g., conducting ice model tests.

The FSICR were further updated with regard to the ice water lines (2006): the ship draught must be between the Upper and Lower Ice Water Line (UIWL and LIWL) when operating in ice. These waterlines may deviate from the ones given in the draught mark. In 2008, the new machinery rules were introduced. This time, the whole structure of the machinery rules was changed. Since the 1971 rules, the machinery rules have been based on the concept of ice torque Q = m·DP2 where m is a class factor and DP the propeller diameter. The new rules only give the ice loads on the propeller and the designer must calculate the propeller scantlings himself/herself.

The most recent change in the FSICR has been to streamline the hull rules. This was done in 2010. When the rules were updated in 1985, several factors were introduced, the background to which is somewhat obscure. The aim in 1985 was not to change the strength level for transverse frames (the new idea of load height had a large impact on longitudinally framed structures). These factors were now put into the context of the newer research results. These 2010 Finnish-Swedish Ice Class Rules form the basis of this paper, and the background to the hull rule formulation is described.

This short elaboration on the evolution of the Finnish-Swedish Ice Class Rules shows that the strength level included in the rules is under constant scrutiny – this scrutiny is provided by the approximately 12,000 ship visits to Finnish ports every winter. Ice damage surveys like the one by Hänninen (2005) suggest that some damage occurs to ship hulls, propellers, and rudders every winter. This indicates that the strength level is not set too high, and as the amount of damage is small it can be considered as not being too low either. This validation of the FSICR every winter since the present strength level was first introduced for transverse frames (other parts have been updated) in 1971 makes the Finnish-Swedish Ice Class Rules very robust and reliable.

3. STRUCTURE OF FINNISH-SWEDISH ICE CLASS RULES

Four ice classes are defined in the Finnish-Swedish Ice Class Rules. These are in order of strength from high to low: IA Super, IA, IB, and IC. The strength level in each ice class corresponds roughly to the loading from a certain level of ice thickness: for IA Super 1.0 m down to IC 0.4 m. The ice thickness for IA Super is higher than the maximum level of ice thickness observed in the Baltic outside the fast ice zone. The design point for each ice class is a collision with a channel edge (as ships are escorted in the worst ice conditions) or with a consolidated layer of older ridges. As the consolidated layer of ridges may be 1.8 times thicker than the level ice at the same location, and the maximum average level ice thickness in the middle of the sea basins is about 60 cm, this results in about 1.0-m-thick ice. It should be noted that this design scenario does not state ship speed – it is considered that no speed restrictions should exist, as this would handicap much of the navigation in ice. It is still somewhat unclear, however, which ship-ice interaction scenario causes the highest loads. This uncertainty in the design scenario is discussed briefly in Section 4.5.

The Finnish-Swedish Ice Class Rules (TRAFI 2010) consist of three main parts: rules for the performance of ships in ice, rules for hull strength, and rules for machinery strength. There are also regulations for ice class draught, rudder and steering arrangements, and some miscellaneous machinery requirements. The topic in this paper is the regulations concerning hull strength, and the quantities included in the hull rules are now described.

3.1 Ice Pressure and Ice Load

The definition of the ice load is the most important part of the hull rules. The ice load in the FSICR is, in principle, defined so that the ice pressure is constant for all classes (nominal ice pressure p0) and the load height h is the class factor (from 0.35 m for ice class IA Super to 0.22 m for ice class IC). The magnitude of the nominal ice pressure is 5.6 MPa.

3.2 Scantlings

Once the ice load is specified and the limit state is defined, the scantlings can be calculated. The limit state used in the Finnish-Swedish Ice Class Rules is the yield limit – consequently, only the elastic response of the structures needs to be derived. The plate thickness equations are based on a similar equation used to design car decks under tire loading. The frame equations are based on simple beam formulation.

3.3 Some Observations

The structure of the Finnish-Swedish Ice Class Rules is quite simple and the flow of the calculations is easy to follow and perform. Some observations can be made, however, concerning the structure of the rules and the main points of the rule rationale. These observations include:

a) The ice load is independent of the hull shape. The simplicity of the formulation is the main reason for this, and not much knowledge exists on the effect of the hull shape. The load is constant along the whole bow and can be qualitatively justified by the fact that two physical effects influence the load: speed of indentation into the ice and frame angle. The indentation speed depends on the projection of speed on the shell normal. This decreases when moving from the stem towards the bow shoulder area. The frame normal angle (frame inclination on the vertical plane including the frame normal) is usually greatest at the bow and decreases towards the shoulder area. The ice load increases with the indentation speed, whereas the load decreases with the frame normal angle – these effects have opposite trends along the bow waterline.

b) According to some other ice class rules, the longitudinal location of the structures influences the required scantlings. This is not the case in the FSICR and similar argumentation as that used above.

c) The ship size description includes both the power and displacement through the factor k. This factor could be called an „aggressiveness factor‟ as it describes the ship inertia and instantaneous speed. The drawback is that there is no theoretical justification for the use of this factor. This question is investigated briefly in Section 4.2.

d) The design point includes the yield as the limit state. If the loading for plating and framing had similar return periods (see Chapter 6), this would induce an unsafe structural strength hierarchy. Plating and frames would be of similar strength and, as frames have less plastic reserve than plating, under ultimate loads the frames would collapse, first leading to greater damage than just plating failing. This is corrected in the present rules, however, with different safety factors for plating and frames.

Now that the structure of the Finnish-Swedish Ice Class Rules is clear, it is time to turn to the background of these hull rules. Here, only the most important elements of the background are described. The aim is to expose the knowledge basis for the different formulations that have been incorporated into the rules.

4. DESCRIPTION OF ICE LOAD

The ice load description is the core of any hull design for ice. The principles of the local ice load description in the Finnish-Swedish Ice Class Rules were described in Chapter 3 and a sketch of these was shown in Fig. 8. Here, the different aspects that have been included in the determination of rule loads are described with the aim of the rule structure becoming more transparent.

4.1 Ice Load Patch Quantities

For structural design purposes, the ice load is commonly assumed to be described by uniform ice pressure, termed pav, on a rectangular load patch of height h and length l. Thus, the total force is F = pav·h·l. The design ship-ice interaction in the Baltic is a collision with a level ice edge, a channel edge when the ship is escorted or with the consolidated layer of an ice ridge. In all these cases, the contact area, i.e., the load patch, is quite narrow in height: in other words, the aspect ratio of the load patch h/l is small. This load description suggests a useful load quantity mentioned already, that of the line load q = p·h. The usefulness of the line load stems from the fact that many structural members are not sensitive to load height but rather to line load. Furthermore, the ice load measurements that are described in more depth in Section 4.5 give the ice load acting on one frame, which is F = p·h·s = q·s (assuming, of course, that the load length is more than one frame spacing s). Thus the line load q can be obtained directly from measurements. The other quantities from the ice load measurements to be used in this analysis are the frame stress on top of the frame flange σFR and the plate stress in the horizontal direction σPL.

The problem with the ice load is that it includes three quantities that are not well known individually. The only quantity that can be known with some confidence is the line load. This line load has been obtained directly from the full-scale measurements and by calculating the ice load from observed ice damage. The insight into the rule on ice load is made even more difficult as different structural members are not sensitive to the total load but to ice pressure (plating), line load (transverse frames), or line load along the whole span in the case of longitudinal frames. An investigation of the ship-ice interaction mechanics suggest that the total ice load F in a collision with ice is mostly dependent on the collision speed and ice thickness (or ice mass in the case of individual ice floes). The local ice pressure p is not dependent – at least not in the same way as the total force – on the ice thickness or collision speed. The indentation rate does influence the contact pressure but only below relatively low indentation rates (some cm/s). Here, the aim is not to discuss the various theories for ice loads or ice pressure but to look at what the rationale is for the rule ice load in the FSICR.

4.2 Ice Load Dependency on Ship Main Particulars

The ice pressure in the FSICR depends on the ship propulsion power and displacement through a factor k, which is given in Eq. (2). The power to be used in the calculation of the factor k is the actual power delivered continuously to the ship propellers (or propulsion). There has been some confusion concerning this actual power and the required rule power, but the situation is clear: the power to be used in calculations is the delivered power PD. The factor k accounts for the possibility of colliding with ice at high speed – thus the propulsion power – and the possibility of penetrating severe ice by using ship inertia. The latter scenario may occur in a channel in which there are thick side ridges, including a consolidated layer. A smaller vessel will not reach the side ridges as she is pushed back to the channel by ice, but a larger ship has a large inertia that will carry her to the consolidated layer of the side ridges.

The background to the quantity k is in the ice damage studies carried out in the late 1960s when the year-round navigation to all Finnish ports increased sharply. The ice damage study was carried out by calculating the strength of ship shell structures in the case of undamaged ships and estimating the ice load causing the observed damages (Johansson 1967). The calculation gave the line load accurately for transverse frames, but as the load height was assumed to be 800 mm, the ice pressure values deduced were low.

The rule ice pressure on the high and low ends of the ship size scale has been criticized: for the smaller vessels, the low end of the pressure – 1.29 MPa in the bow region – is considered low and the fact that there is no upper limit is considered unrealistic this was corrected in the 2010 rules by introducing a ceiling. At the high end, the bow ice pressure would exceed the nominal ice pressure p0 = 5.6 MPa when k = 80.3. This corresponds roughly to a large SUEZMAX tanker of ice class IA Super (Δ = 170 000 dwt, PD = 30 MW).

In order to investigate the size dependency of ice pressure, ice damage data from 2003 (Hänninen 2005) is analyzed. This is the last comprehensive survey of ice damage, a survey that includes ships built according to the 1985 rules. Fig. 10 shows the results of the ice damage survey carried out during winter 2003. The relative abundance of ice damage on different sizes of ships must be balanced with the relative number of this size of ship in traffic. Small ships (less than 2000 dwt) and large ships (more than 20,000 dwt) seem to be more susceptible to ice damage, as their total share of ice damage is 48%. Their share of ship traffic is only 19% however.

> 20,000 dwt 12.3%

10,000-20,000 dwt 2.6%

4000-10,000 dwt 1.7%

2000-4000 dwt 2.2%

< 2000 dwt 6.7%

These figures suggest that there is a slight increase in the tendency for larger vessels to be less susceptible to ice damage, but the accuracy is not high in view of the small amount of ice damage to larger ships. Before a definite conclusion can be drawn from the damage data, the ships that have no ice class or were damaged below the ice belt should be removed from the statistics. Even these data seem to show some bias towards the smaller end of the ship size, as two of the four damaged smaller vessels, measured with dwt, are actually quite large ships, passenger car ferries with a small deadweight. Thus, the small relative increase in the amount of damage in the small ship category is misleading, as half of the ships with a small dwt are actually quite large ships, similar in displacement to bulk carriers of about 6 … 10,000 dwt. The larger ships seem to suffer slightly more ice damage. The relative and weighted amounts of ice damage for different sizes of ships.

Ship dwt

Relative amount of all hull ice damage

Number of damaged ships in each size category

< 2000

30.8%

3.1%

2000-4000

0

0

4000-10 000

38.5%

1.3%

10 000-20 000

15.4%

1.2%

> 20 000

15.4%

5.8%

The dwt is not a good measure of ship size. Thus, a similar analysis of the ice damage frequency should be carried out for the modified data using the factor k. There is no distribution of ships visiting Finnish ports divided according to the factor k however. Thus, only the relative number of damaged ships can be given; see Table 3. As the average value of the factor k can be roughly estimated from the fleet, giving a value of about 7, it can be concluded that there seems to be no trend in the dependency on the factor k.

It is also possible to look for the size dependency in the data from the ice load measurements. There are data from long-term (measurements covering at least one whole winter), full-scale ice load measurements in the Baltic from four ships: IB Sisu, MT Kemira, MS Arcturus, and MT Kashira. A collection of the data from these ships was made by Hänninen (2002). MT Kemira is described in Section 4.5 and, for the purpose here, it suffices to describe the other ships briefly. IB Sisu is a 16.8 MW icebreaker operating in the northern Baltic. The measurements were taken in winters from 1979 to 1985. MS Arcturus is a roro ship operating between Central Europe and the Gulf of Finland. Her ice class is IA Super. Three winters were covered in the local hull ice load measurements: 1985, 1986, and 1988. MT Kashira is a tanker operating in the Gulf of Finland and the Arctic (only Baltic measurements are considered here). Her ice class is UL, corresponding roughly to IA Super. The measurement years were 1985-1990. In this context, the frame ice load measurement results are investigated. The ships are arranged in order of severity of the ice conditions encountered (IB Sisu at the top). Furthermore, the ice load values given include the measured maximum and predicted most probable value in 1000 days. Taking into account the different areas of operation, no trend vs. k in the measured or predicted maxima can be discerned.

The conclusion from the analysis of the size effect is that there is no strong suggestion to make changes to the present rule formulation. The only suggestion that can be made at this time is that a ceiling on the ice pressure is introduced. It is not likely that the frame ice load would increase beyond the measured maximum of both IB Sisu and MT Kemira, i.e., about 2000 kN/m. With a maximum rule load height of 35 cm, this means that the maximum pressure could be set as the nominal ice pressure p0 = 5.6 MPa. This decision is also justified by the fact that ships with a very large k factor value are large tankers (or bulkers) and are thus not likely to navigate in ice for as long or as aggressively as IB Sisu or MT Kemira.

4.3 Load Length Dependency

The ice pressure in the rule formulation depends on the load length. This is taken into account in the 2010 FSICR with a coefficient ca given in Eq. (3), but it is repeated here for clarity

, max 1.0 and min 0.35; l0 = 0.6 m. (7) aallc0

This coefficient is stated in the rules to account for “the probability that the full length of the area under consideration will be under pressure at the same time.” It greatly resembles the dependency of the average ice pressure on the load area, the formulation of which is used in offshore applications in particular – but for ships in first year ice, the load is more line-like and thus the dependency on the load length is used rather than the dependency on area. As the upper limit of the coefficient ca is set to be reached with length l0 = 0.6 m, it is clear that the coefficient is only important for longitudinal framing systems. As explained in Section 3.1, each structural member has an associated load length. The plate response in longitudinal framing systems is sensitive to load length only up to a length of roughly twice the frame spacing (roughly the distance s from the point under study). Thus, the length used with the plating design in longitudinally framed structures is l0 = 1.7s.

The present form of the length-related coefficient is based on a set of measurements from ships in the Arctic and the Baltic (see Frederking & Kubat 2005 for icebreakers Louis S. St. Laurent, Oden, and Polar Sea; Kujala and Vuorio 1986 for icebreaker Sisu; Kujala 1991 for ice damage data). All these data include simultaneous measurements (except, of course, the data based on ice damage data) of ice load on different lengths.

4.4 Equivalent Ice Pressure for Plating

Ice load measurements in the Baltic have indicated that the maximum stress in a frame, the maximum stress at the adjacent plate field, and the ice force acting on the same frame are not compatible with each other if the ice pressure distribution is assumed uniform. It has been noted that especially the plate stress is lower than the calculated value using the measured ice load, assuming the ice pressure to be uniformly distributed (and a reasonable load height). This incompatibility has been associated with the flexibility of plating reducing the ice pressure in the middle of the frame span; see Vuorio et al. (1979). The phenomenon has been explained by the similarity with stiffened plating resting on an elastic foundation (Varsta 1983); if a framed structure is pressed against an elastic or Winkler (i.e., foundation sustaining only pressure, not shear forces) foundation, the contact pressure varies so that the pressure is lowest at the mid-span between the frames. This analogy from elastic foundations may be valid for transversely framed structures for which the line-like load continues over several frame spaces.

4.5 Statistics of Ice Loads

The local ice loading process on the ship shell consists of separate peaks. Each of the peaks represents an impact with an ice edge at a separate cusp. This ice edge shape is formed by the breaking pattern by repeated breaking of these cusps. When the ship hits an ice cusp, the ice edge is first crushed and the contact force then increases with the contact area. The maximum force reached is determined by the bending strength of the ice cover. At any fixed location on the hull, the loading process, consisting of separate peaks, seems quite stochastic; this process has been described numerically by Su et al. (2010). The stochastic ice loading process must be described by statistical distributions – and it is natural to use the load distribution that is the distribution of the load peaks. Often, however, the maximum load values are described by the maximum value observed within some selected time period. Most of the Baltic ship ice load measurements that have been carried out for this period are either 12 h or 24 h.

The analysis of the design ice loading requires knowledge about the statistical properties of ice loads. The theoretical research into this matter is just commencing, so only knowledge that exists comes from measurements. The initial theoretical analysis (see, e.g., Su et al. 2010) suggests that the origin of the statistics of the local ice loading can be divided into two sources: external and internal. The external source is the variability of encountered ice conditions. The internal source is created by the variability of the breaking pattern of ice; this breaking pattern causes variations in the ice load even if the ice strength and thickness parameters would be constant. Here, however, only measured data are considered and the considerations are made somewhat more concrete by considering the measured data from MT Kemira.

MT Kemira is a chemical carrier that operates between the ports Kokkola and Uusikaupunki in Finland and central European ports. She visits Finland often, with a rotation of about one week. Her ice class is IA Super with a length of Lpp = 105.0 m, beam B = 17.5 m, power 3400 kW, and the deadweight is 5800 dwt. This ship was selected for measurements as she navigates regularly in the iciest conditions in the Baltic – thus her operational spectrum can be assumed to represent a typical Baltic merchant ship of her size and shape, naturally. Kemira‟s ice load data are described by Kujala (1989), Gyldén and Riska (1989), and Muhonen (1991, 1992).

MT Kemira was instrumented to measure the ice load on one frame as well as the plate and frame stresses (labeled PL and FN). The frame load measurement was carried out using shear strain gauges attached at the neutral axis of the frame: the difference between the two shear stresses on the same frame is proportional to the load on the frame between these gauges.

The measurement system does not record the time histories of the signals but, instead, calculates two histograms for each measurement channel and measurement period. One histogram consists of the digitized samples from each measurement channel with a sampling rate of 100 Hz. The other histogram includes the ice force peaks, with each peak calculated using the Rayleigh separation for the peaks (with 25% of the peak value as the separator). These histograms were stored after each 12 h period throughout the winter for the years 1985-1991.

The data used here are derived from the maxima of each 12 h peak distribution. Only the maxima from periods when the ship was sailing in ice are included. Before using the data, some processing has to be done. The first assumption was to omit the smallest peaks, as these most probably contain some electronic noise as well as very small ice-induced responses. Specifically, peaks smaller than 20 kN are ignored. The other assumption was to collate the data so that from each 12 h period, the maximum frame load of the measured peaks at the same frame was taken to the data. Thus, the data consist of MAXFFR1,FFR2, MAXFFR6,FFR7, and MAXFFR11,FFR12, where FFR stands for frame load in the gauge denotation used in the measurements.

The advantage of using the Kemira data from these seven years is that they integrate different winters (three severe, two average, and two mild) as well as the ice thickness development during each winter. The disadvantages include the dependency on MT Kemira’s size, hull shape, and navigation spectrum (MT Kemira mostly visited the ports in Uusikaupunki and Kokkola in Finland). The fact that during each 12 h period the ship did not navigate in ice the whole time – in fact, the ice navigation was less than 12 h for most of the 12 h periods (see Muhonen 1992, Appendix 3). This effect has been investigated by Hänninen (2002): the effect was noted to be quantitatively important but did not change the data qualitatively.

In order to have some insight into the data, the statistical distributions are investigated briefly. The basic distribution is the distribution of peaks, fP(q), from a time period of 12 h. The bow frame load data from 3/16/1991 (A.M.) are used for this demonstration, data from Muhonen (1992). On this day, the ship was operating through the Quark having left Kokkola at midnight. The maximum undeformed level ice thickness during this time and in this area was about 50 cm.

This short analysis of the maxima from a single measurement campaign raised some questions concerning the extrapolation of loads and the connection of the estimated pdfs to ice conditions. It would be interesting to study the relationships between the peak distribution, distribution of maxima from certain time periods, and the method of extrapolation to lifetime values. This falls outside the focus of the present article, however, and warrants a study in itself.

To conclude the study of the statistics on ice loads, the ice load measurement results from measurements carried out with four different vessels in the Baltic are described; see Hänninen (2002). The results of a statistical curve fit using the Gumbel I extreme value distribution on the data from these four ships are shown in Fig. 21. The load q is presented versus the return period T of the load. As the statistical distribution used is Gumbel I, the most probable maximum load from period T is given in Eq. It is noteworthy that the values of the constant α, which are proportional to the standard deviation of the distributions, are very close to each other. The values from the northern Baltic (Sisu and Kemira) and values from the Gulf of Finland (Arcturus and Kashira) are also very close to each other, respectively. These observations would warrant closer investigation.

5. STRUCTURAL RESPONSE

The response in terms of the maximum stress of plating and frames is investigated briefly in order to analyze the design point using different limit states in Chapter 6. Several studies have been carried out investigating the plastic and ultimate response of ship hull structures under ice loads – one of the first ones is Johansson (1967). The other studies include studies in the project SAFEICE (Kujala, ed. 2007), and damage studies Varsta et al. (1978), Ranki (1986), Hayward (2001), Daley (2002a,b), Valkonen (2006), and Kaldasaun (2010). The development of formulations for the elastic and plastic response of plating and frames would be a main topic in itself; here, only some formulations are selected for use in investigating the design point.

5.1 Plate Response

In order to calculate the elastic response of plating, the formula developed for the Finnish-Swedish ice class rules (Trafi 2010) can be used.

5.2 Framing Response

The equations corresponding to the plating response, especially for the plastic frame response, are more complicated than the corresponding ones for plating. The different boundary conditions of continuous beams and end brackets cause variation as well as different deformation mechanisms. If only pure bending deformation is taken into account, the formulation is often too simple, as at least the shear deformation should also be taken into account. There is no clear formulation for the inclusion of shear deformation, and as there is quite great uncertainty in describing the ice loads, the frame response in this study is obtained by only describing the simple bending response.

A reference for this formulation is Hayward (2007). The frame response formulation is simple, as it does not consider the shear deformation or the material post-yield behavior. In order to check the validity of the frame plastic formulation, the results in Eqs. (27-30) are compared with results from FE calculations conducted by Det Norske Veritas (Holtsmark & Strömme 2004).

6. DESIGN POINT IN THE ICE RULES

6.1 Description of the Design Point

The aim of the structural response formulations in Chapter 5 is to derive a relationship between the limit response, the scantlings, and the load. Here, only a transversely framed hull structure is investigated. The limit response that is investigated includes both elastic and plastic limit states.

Definition of the limit states for plating and frames

Limit state (label)

Plating

Frames

Elastic (Y)

Stress reaching the yield stress σy somewhere in the plate

Stress reaching the yield stress σy somewhere in the frame

Plastic (P)

Stress distribution reaching full plasticity somewhere in the plate; Permanent deformation still zero

2-hinge formation at the frame supports

Ultimate (U)

Permanent deformation (wP) reaching a specified value

3-hinge formation at the frame supports and the mid-span

The reason the permanent deformation is specified for plating but not for frames is that no easy formulations are available for frame deformation, and the simple formulation gives deformation that is too small. Some formulations for frame deformation exist, like the one given in NORSOK N-004 (Appendix A, section 3, eq. A.3.28), but it was decided to use only loading for simple 2- and 3-hinge formulations for this scoping study. These also ignore shear deformations. The plastic reserve due to membrane stresses is large for plating, and simple formulations for the permanent deformation of patch-loaded plates exist (see, for example, Hayward, 2007)

Three different hull areas are investigated: bow, midship, and stern, using the measured loading as representative for these areas. This is not strictly correct, as the bow, as well as the stern area, contains a range of hull angles. The stern measurement position is where the frame angle is almost vertical – thus these loads can be considered conservative for what influences the hull shape. The bow load measurement position is where the waterline angle is α = 19o and the frame angle is β = 24o. This position is at about the centre of the bow area length and thus the measured loads can be considered characteristic of the whole Kemira bow, on average.

The ice class of MT Kemira is IA Super. The calculated scantlings are compared with the minimum values required by the Finnish-Swedish Ice Class Rules (Trafi 2010) and the ice class PC6 of the International Association of Classification Societies Polar Class Rules (IACS 2007). These two classes are generally considered as approximately equivalent – this equivalency is one of the motives of this study, as the limit state in the design point in the Finnish-Swedish Ice Class Rules is the yield (Y) limit state and in the IACS rules it is a plastic limit state (either fully plastic or ultimate, it is not clear which one is used). The required scantlings of MT Kemira for these two classes are given in Table 10. It should be noted that the required section modulus for the ice class PC6 is a plastic section modulus while for ice class IA Super, it is elastic. The plate thickness given is the net thickness.

6.2 Design Point Calculations

When the load statistics and strength equation have been derived, the required scantlings can be determined versus the return period using the limit state as a parameter. This plot for plate net thickness for the bow is shown in Fig. 24 – using MT Kemira as an example. This plot also shows the required net plate thickness according to ice class IA Super and PC6. Similar plots can be made for the midbody and stern hull regions.

6.3 Design Point

The calculations shown in Section 6.2 can be used to compare the scantlings required by the ice classes IA Super and PC6. This is done by determining the required plate thickness or frame section modulus and comparing these with those thicknesses or frame section moduli corresponding to the different limit states once in the ship‟s lifetime (1250 days in ice per lifetime). It is immediately clear that the built scantlings correspond closely to the plastic limit state, whereas the scantlings given by the IA Super requirements correspond to some permanent deformation of the plating and to about 2-hinge formation for framing. The plate thickness for the midbody and stern hull regions required by the ice class PC6 seem to be somewhat low as dents as deep as the plate thickness are allowed.

7. FUTURE CONSIDERATIONS

All ice class rules, including the Finnish-Swedish Ice Class Rules, need updating when a new insight into the ice loading, structural response, or design point is obtained. A sound development environment would be one in which the rules were critically appraised quite often by including new knowledge of the loading and investigation, if all possible structural responses had been taken into account. The feedback from the built ships is also important and a systematic way to collect the feedback in the form of ice damage should exist. If it relied only on the feedback from the practice (ice damage) and this feedback urged changes in the loading, structural response formulations, or design point, then the changes would be somewhat arbitrary – at worst this could lead to „design by disaster.‟ Here, only two possible items to be taken into account in the future development of the FSICR are mentioned: that of the hull shape in determining the ice loads and the validation of the design point.

7.1 Hull Shape Effect on Ice Loads

A recent research project, SAFEICE, produced two slightly different results from the influence of hull angles. The theoretical calculations given in Valanto (2005) suggested that the influence of hull lines is, at most, quite small, as there are two competing physical processes that influence the hull angle dependency. The waterline angle α influences the indentation speed into ice and the normal frame angle βn influences the amount of the total force that is required to bend the ice. The hull ice load measurements on a model scale (Izumiyama 2005) suggested that the ice load is inversely proportional to βn.

In order to investigate the influence of hull lines on the contact force, some formulations to account for the hull angles are investigated here. Four different ships, from which there are full-scale data (IB Sisu, MV Arcturus, MT Kemira, and MT Kashira), are used as a comparison.

There is much scatter in shape factors. All give similar results for MV Arcturus and MT Kashira, but the two rule formulations give much deviation from the other factors for IB Sisu. As the calculated force is based on the MT Kemira results, it would be surprising if the values obtained for IB Sisu were close to the actual measured value, as IB Sisu is navigating much more aggressively than MT Kemira. Thus, if IB Sisu were navigated as MT Kemira, the expected forces would be much less than the actual observed value in IB Sisu measurements. This gives plausibility to the factors given by Kujala (1994) and Varsta (1984). It is difficult to draw conclusions however.

7.2 Strength Level Validation

The required strength level of the ice class rules – or any structural rules for that matter – is not something exact that can be calculated from the first principles. It depends on how much damage, in a very general sense, is accepted. How close the different limit states are in terms of loading also influences the selection of the required strength level. If there is not much plastic reserve after the yield point, i.e., the load causing collapse is close to the load causing the first yield, then the required strength level must be set quite high. Methods have been developed to evaluate the probability of failure or structural reliability, but all these require knowledge about the ice load statistics and this is only known about in a rudimentary fashion. The only method left to assess the adequate strength level is feedback from past experience.

The feedback from past experience must be collected in an organized fashion. One way to do this is to look at accidents that happen during the winter season compared with other seasons. As the winter season is about four months long, there seems to be a slightly increased risk of accidents during the winter. Unfortunately, the HELCOM database does not allow the investigation if ice cover was the cause of the accident – but judging from the locations of the accidents, ice rarely caused any major accidents. This may be a matter of reporting.

Three major Baltic ice damage data collections have been carried out: Johansson (1967), Kujala (1991), and Hänninen (2005). The first one of these resulted in the ice class rules of 1971. It also defined the strength level that has basically been followed since, even if several modifications to the loading have been made. The next damage data collection occurred at the end of the 1980s and thus not many ships built under the 1985 FSICRs were included. This damage data collection indicated that longitudinally framed ships are quite susceptible to ice damage – this drawback was corrected in the 1985 rules when the load height was reduced while keeping the ice loading constant.

Finally, the winter 2003 ice damage data collection covered a winter that was more severe than almost ten preceding ones. Quite a large amount of damage of different kinds was observed, with more than 100 cases. Many cases related to practices that had become slacker during the preceding mild winters however. During the winter, 27 cases of hull damage were observed; see Fig. 28. This number represents about 0.2% of all the port calls during the winter. Hull ice damage cannot thus be considered very frequent. This indicates that the strength level of the rules is at an acceptable level.

The final matter to be investigated is the advantage that can be gained by having an ice class. The traffic restrictions to ports in the Gulf of Finland often vary from country to country. For the present winter, the ice class requirement for Russia and Estonia has been the minimum ice class of IC, whereas to Finnish ports the requirement has been IA. Is there an increased risk of ice damage using lower ice class tonnage?

This question can be answered by looking at the ice damage data from winter 2003, Hänninen (2005). These data only include ships bound to and from Finnish ports (100 cases of ice damage, 10,000 port calls with about 1000 different ships). The ice damage data suggest that the majority of ships damaged had an ice class of IA or IA Super (59% of all cases). If the ice damage is calculated per port call, then the probability of ice damage for IA and IA Super ships is 0.7%, whereas the same figure for II and IC ice classes is 6.0%. The probability of ice damage is about nine times greater for low ice class ships: ice class surely has an impact on susceptibility to ice damage.

8. CONCLUSION

The Finnish rules for winter navigation have a long history stemming back to 1890 when the first rules were given. These rules, currently called the Finnish-Swedish Ice Class Rules (FSICR), are an integral part of the winter navigation system. Finnish trade is highly dependent on a marine connection, as about 90% of the trade is transported by sea and, as all Finnish ports are ice bound in an average winter, it is clear that winter navigation is crucial to Finland. This situation is less dominant for Sweden, as Sweden has ice-free ports on its west coast.

Winter navigation depends on ice-strengthened merchant ships to navigate efficiently without undue stoppages in traffic, which is often line traffic with fixed schedules. In order to achieve a winter navigation system that runs smoothly, all ships that fulfill the requirements for a minimum amount of cargo and an ice class are given an icebreaker escort to and from the Finnish and Swedish ports. Furthermore, the safety of shipping is ensured by an adequate strength level given in the ice class rules, and the need for icebreakers is optimized by requiring some minimum ice-going performance from merchant ships in the ice class rules.

Winter navigation – or rather year-round navigation – is considered to have started when the steamer SS Express II sailed on her maiden voyage on 12/15/1877 from the Finnish port of Hanko to Stockholm in Sweden. Since those days, winter navigation has expanded, and since the late 1960s all Finnish and Swedish ports have been kept open. In 1966, Finland and Sweden agreed to develop the winter navigation system together – this resulted in, among other things, the first Finnish-Swedish Ice Class Rules in 1971 and an order for an icebreaker series, the Urho class. Three of these icebreakers were built for Sweden and two for Finland.

The 1971 ice class rules were based on an extensive ice damage survey and resulted in the required strength level for the hull and machinery being established. These strength levels (each ice class has its own strength level) have essentially remained unaltered since then. The 1971 ice class rules also defined the ice load explicitly for the first time. Although the rules have seen several changes, the rule structure is similar, even in the newest 2010 rules. Experiences of ships built for the Finnish-Swedish ice classes have made several changes necessary to the rules since 1971. In 1985, the new knowledge about ice loads was incorporated into the rules and, at the same time, the strength level for longitudinally framed ships was made adequate. In 1999 and, slightly modified, in 2002, the new requirements for ice-going performance were given. The machinery rules were amended in 2008 and, finally, in 2010, the hull rules were streamlined to remove any inconsistencies that had remained in the load and strength definitions.

The Finnish-Swedish Ice Class Rules can be considered very robust, as every winter about 10,000 ship visits test the rule strength and ice performance levels. The present winter (2011) has been a little more severe than the long-term average. Some rudder and propeller damage has been observed, but the feeling is that the amount of ice damage is not alarmingly high. The ice performance of ships, on the other hand, is an issue that has to be investigated, as the traffic to the northernmost ports in Finland and Sweden ran very slowly for several weeks in winter 2011. Shipowners claim that they have lost several hundred ship-days in delays. Again, feedback from the existing winter navigation system may give an initiative for changes. This quick accommodation of the ice class rules to the „facts of life‟ is the strength of the Finnish-Swedish Ice Class Rules.

Bibliography



     

    Daley, C. 2002a: Derivation of Plastic Framing Requirements for Polar Ships. Marine Structures, Vol. 15, pp. 544-559.

    Daley, C. 2002b: Application of Plastic Framing Requirements for Polar Ships. Structures, Vol. 15, pp. 533-542.

    Frederking, R. & Sudom, D. 2008: Local Ice Pressure Distributions during 1990 Hobson's Choice Ice Island Multi-Year Ice Indentation Tests. Proc. 19th IAHR Symposium on Ice, Vol. 2, pp 815-827, Vancouver, B.C. Canada, 2008.

    Frederking, R. & Kubat, I. 2005: Definition and Application of Pressure-area Relations and Line Load Distributions. Deliverable D3-2 from the project SAFEICE. In Kujala, ed. (2007), Vol. 1, pp. 134-150.

    Gyldén, R. & Riska, K. 1989: Ice Load Measurements onboard the MS Kemira, Winter 1990. Helsinki University of Technology, Laboratory of Naval Architecture and Marine Engineering, Report M-93, 13 p. + app.

    Hayward, R. 2001: Plastic Response of Ship Shell Plating Subjected to Loads of Finite Height. ME thesis, Faculty of Engineering and Applied Science, Memorial University of Newfoundland, Canada, 113 p. + app.

    Hayward, R. 2007: Principles of Plastic Design. Report D8-4 from the SAFEICE project, 21 p. (see also Kujala, ed. 2007).

    Holtsmark, G. & Strömme, H. 2004: The Capacity of Panel Stiffeners Subjected to Lateral Pressure Loads. Det Norske Veritas, Technical Report No. 2004-0168, 39 p. 43

    Hänninen, S. 2002: The Use of Statistical Methods in Determination of Design Ice Load on Ship Hull Frame in the Baltic Sea in Finnish. M.Sc. thesis, Helsinki University of Technology, Department of Mechanical Engineering, 100 p. + app.

    Hänninen, S. 2005: Incidents and Accidents in Winter Navigation in the Baltic Sea, Winter 2002-2003. Winter Navigation Research Board, Research Report No. 54, Helsinki, 39 p.

    IACS 2007: Requirements Concerning Polar Class. International Association of Classification Societies, 38 p.

    Izumiyama, K. 2005: Description of Local Ice Loading and Hull Area Factors based on Model Test Results. Deliverable D3-1 from the SAFEICE project. In Kujala, ed. (2007), Vol. 1, pp. 115-133.

    Johansson, B. 1967: On the Ice Strengthening of Ship Hulls. International Ship Building Progress, vol. 14, pp. 231-245.

    Kaldasaun, J. 2010: Risk-Based Approach for Structural Design of Ice-Strengthened Vessels Navigating in the Baltic Sea. M.Sc. thesis, Aalto University School of Science and Technology, Department of Applied Mechanics, 83 p. + app.

    Kujala, P. & Vuorio, J. 1986: Results and Statistical analysis of Ice Load Measurements on Board Icebreaker

    Sisu in Winters 1979 to 1985. Winter Navigation Research Board, Research Report No 43, Helsinki. 52 p. + app. 73 p.

    Kujala, P. 1989: Results of Long-Term Ice Load Measurements onboard Chemical Tanker

    Kemira in the Baltic Sea during the Winters 1985 to 1988. Winter Navigation Research Board, Research Report No. 47, 55 p. + app.

    Kujala, P. 1991: Damage Statistics of Ice-Strengthened Ships in the Baltic Sea. Winter Navigation Research Board, Research Report No. 50, Helsinki, 66 p.

    Kujala, P. 1994: On the Statistics of Ice Loads on Ship Hull in the Baltic. Acta Polytechnica Scandinavica, Mechanical Engineering Series No. Me 116, 98 p.

    Kujala, P. ed. 2007: Increasing the Safety of Icebound Shipping, Final Scientific Report, Vol. 2. Helsinki University of Technology, Ship Laboratory, Report M-302, 347 p.

    Lensu, M. 2002: Short Term Prediction of Ice Loads Experienced by Ice Going Ships. Helsinki University of Technology, Ship Laboratory, Report M-269, 59 p.

    Lindholm, J-E., Riska, K., & Joensuu, A. 1990: Contact between Structure and Ice – Results from Ice Crushing Tests with Flexible Indentor. Helsinki University of Technology, Laboratory of Naval Architecture and Marine Engineering, Report M-101, Otaniemi, 30 p. + app.

    Muhonen, A. 1991: Ice Load Measurements onboard the MS Kemira, Winter 1990. Helsinki University of Technology, Laboratory of Naval Architecture and Marine Engineering, Rpt M-109, 24 p. + app.

    Muhonen, A. 1992: Ice Load Measurements onboard the MS Kemira, Winter 1991. Helsinki University of Technology, Laboratory of Naval Architecture and Marine Engineering, Rpt M-121, 26 p. + app.

    NORSOK Standard N-004 2004: Design of Steel Structures. Standards Norway, 287 p.

    Ranki, E. 1986: Determination of Ice Loads from the Permanent Deformations of Shell Structure in Ships in Finnish. Lic. Tech. thesis, Helsinki University of Technology.

    Riska, K., Rantala, H., & Joensuu, A. 1990. Full Scale Observations of Ship-Ice Contact. Helsinki University of Technology, Lab. of Naval Architecture and Marine Eng., Report M-97, Espoo, 1990, 54 p. + 293 app.

    Riska, K. & Windeler, M. 1997: Ice-Induced Stresses in the Shell Plating of Ice-Going Vessels. Helsinki University of Technology, Arctic Offshore Research Centre, Rpt M-219, Otaniemi, Finland, 34 p.

    Riska, K., Uto, S., & Tuhkuri, J. 2002: Pressure Distribution and Response of Multiplate Panels 44

    under Ice Loading. Cold Regions Science and Technology, 34(2002), pp. 209-225.

    Riska, K. 2006a: Factors Influencing the Power Requirement in the Finnish-Swedish Ice Class Rules. Report to the Finnish Maritime Administration, to be published in the Winter Navigation Research Board report series, 63 p.

    Riska, K. 2006b: Ice Classification of Large Vessels. Conference on the Winter Navigation on the Baltic Sea, ICEDAY 2006, 9-10 February, 2006, Kemi.

    Riska, K. 2007: Application of the SAFEICE Project Results in Developing the Finnish Swedish Ice Class Rules. Deliverable D7-3 from the SAFEICE project. In Kujala, ed (2007), 27 p.

    Riska, K. 2011: Design Point in Ice Class Rules. Report to the Finnish Transport Safety Agency, to be published in Winter Navigation Research Board report series, 34 p.

    RMRS 2007: Rules for the Classification and Construction of Sea-Going Ships. Russian Maritime Register of Shipping, Vol. 1, Saint Petersburg, 461 p.

    Siivonen, O. 1977: The Development of Finnish Ice Class Rules. Ice, Ships and Winter Navigation Symposium, Oulu 16-17.12.1977, pp. 154-171.

    Sodhi, D, Takeuchi, T., Nakazawa, N., Akagawa, S., & Saeki, H. 1998: Medium-scale indentation tests on sea ice at various speeds. Cold Regions Science and Technology, Vol. 28, pp. 161-182.

    Su, B., Riska, K., & Moan, T. 2010: Numerical Simulation of Local Ice Loads in Uniform and Randomly Varying Ice Conditions. Cold Regions Science and Technology, Vol. 65, pp. 145-159.

    Trafi 2010: Ice Class Regulations 2010: “Finnish-Swedish Ice Class Rules 2010.” Finnish Transport Safety Agency, 23.11.2010 TRAFI/31298/03.04.01/2010, 48 p.

    Tuhkuri, J. 1993: Laboratory Tests of Ship Structures under Ice Loading, Vol. 1. Helsinki University of Technology, Ship Laboratory, Report M-166, Otaniemi, 171 p.

    Uto, S. 2000: Influence of Plate Rigidity on Ice Loading under Line-Like Contact between Ice and Stiffened Ship Hull Structure. Report M-254, Arctic Offshore Research Centre, Helsinki University of Technology, 47 p.

    Valanto, P. 2005: Spatial Distribution of Numerically Predicted Ice Loads on Ship Hulls in Level Ice. Deliverable D6-3 from the project SAFEICE. In Kujala, ed. (2007), Vol. 2, pp. 25-51.

    Valkonen, J. 2006: Determination of Ship Ice Load from Hull Ice Damages. M. Sc. thesis, Helsinki University of Technology, Ship Laboratory, 98 p. + app.

    Varsta, P., Droumev, I. & Hakala, M. 1978: On Plastic Design of an Ice-Strengthened Frame. Winter Navigation Research Board, Report No. 27, 54 p.

    Varsta, P. 1983: On the Mechanics of Ice Load on Ships in Level Ice in the Baltic Sea. Technical Research Centre of Finland, Publication No. 11, 91 p.

    Varsta, P. 1984. Determination of Ice Loads Semi-empirically. Ship Strength and Winter Navigation, VTT Symposium 52, pp. 177-182.

    Vuorio, J., Riska, K., & Varsta, P. 1979: Long-Term Measurements of Ice Pressure and Ice-Induced Stresses on the Icebreaker SISU in Winter 1978. Winter Navigation Research Board, Research Report no. 28, 50 p.