Microsoft Word - numero_53_art_30_2840 A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 394 Applicability of Vyalov’s equations to ice wall strength estimation A. Kostina, M. Zhelnin, O. Plekhov, I. Panteleev Institute of Continuous Media Mechanics of the Ural Branch of Russian Academy of Science, Russia kostina@icmm.ru, http://orcid.org/0000-0002-5721-3301 zhelnin.m@icmm.ru, http://orcid.org/0000-0003-4498-450X poa@icmm.ru, http://orcid.org/0000-0002-0378-8249 pia@icmm.ru, http://orcid.org/0000-0002-7430-3667 L. Levin, M. Semin Mining institute of the Ural Branch of Russian Academy of Science, Russia aerolog_lev@mail.ru, http://orcid.org/0000-0003-0767-9207 seminma@outlook.com, http://orcid.org/0000-0001-5200-7931 ABSTRACT. A simple analytical relations are commonly used in engineering practice to calculate ice wall thickness. One of them is Vyalov’s relation that takes into account the features of a real technological process of tubing lining and the inelastic deformation associated with frozen soil creep. An estimation of applicability and margin of safety of this equation is an issue of engineering mechanics. In this paper, we propose a mathematical model for description of ice wall deformation under natural external loading and present the results of the computational experiments in which an optimal thickness for the ice wall is determined. Based on this simulation, we modify the existing analytical relation, which makes it possible to calculate the thickness of an ice wall of unlimited height. KEYWORDS: Ice wall; Frozen soil; Mine shaft; Vyalov’s equation. Citation: Kostina, A., Zhelnin, M., Plekhov, O., Panteleev, I., Levin, L., Semin, M., Applicability of Vyalov’s equations to ice wall strength estimation, Frattura ed Integrità Strutturale, 53 (2020) 394-405. Received: 20.05.2020 Accepted: 01.06.2020 Published: 01.07.2020 Copyright: © 2020 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION he ice wall is an important engineering structure the integrity of which determines the safety of sinking operations in vertical mine shafts under construction. The effectiveness of the ice wall as a temporary protective lining depends on the adequate evaluation of its thickness. Optimal assessment of the wall thickness comes from the need to reduce the cost of mining operations at the design stage. This will help to exclude the failure of an ice wall and the breakthrough of groundwater into the mine. Furthermore, the savings in artificial freezing costs may be significant. Engineering calculations for ice walls are currently performed for two limit states - strain condition and strength condition. In strain calculations, it is necessary to determine the minimum thickness of the ice wall at which its deformation in the design phase does not exceed the value acceptable for shaft construction and promotes no dangerous movements of freezing columns. The objective of the limit strength design is to calculate the optimal minimum thickness of the ice wall which T https://youtu.be/LmyKQbPy_Ik A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 395 prevents its failure under external loads, i.e. no cracks occur in the wall. At this thickness value, the ice wall is in the ultimate stress-strain state, i.e., in the state where the stresses at a given instant of time do not exceed the frozen soil strength at the same instant. The most common methods for calculating of ice wall thickness are those proposed by Domke [1] and Lame-Gadolin [1]. The Lame-Godolin formula was obtained in a linear elastic approximation and did not prove effectiveness at soil depths of more than 50 m. The Domke equation takes into account the plastic deformations of frozen soil, the onset of which is determined according to the criterion for the difference between the largest and smallest principal normal stresses. Because the strength of saturated soils increases with an increase in the mechanical pressure up to a certain value, the application of the Mohr-Coulomb and Druker-Prager criteria makes it possible to calculate plastic strain with higher accuracy and to take into account the fact that the frozen soil strength increases with decreasing temperature [2]. The triggering of high hydrostatic pressure leads to melting and failure of the ice contained in the pores. This has given impetus to different modifications of these criteria [3]. The thermo-hydro-dynamical models of soil freezing which take into account plastic strains are presented in [4-7]. In [4-6], the so-called Barcelona Basic Model [8] is used to describe the plastic deformations of saturated frozen soils. According to this model, the stress-strain state of soils is determined by analyzing two state variables – effective stress and a parameter that describes the moisture migration effect. For the non-frozen soil, the Barcelona Basic Model is simplified and reduced to the modified Cam-Clay model. The thermo-hydro-dynamical models based on this approach can be employed to consider the influence of moisture migration and mechanical pressure on the plastic deformation of frozen soils. In [6], the Barcelona Basic Model is generalized to the case of large deformations. Another way [7] has been proposed to derive the constitutive relations for describing the elastoplastic behavior of frozen soils in the framework of a thermodynamic approach. It has been suggested in [9] that the yield surface defined in terms of a double clay hardening model can be used to assess plastic strains. Elastic and plastic strains characterize the instantaneous response of the material to the applied loads. However, due to the pronounced rheological behavior of saturated soils under long-term loads, their deformation increases and strength reduces. This feature of soils is of prime importance for shaft sinking where the loads from the surrounding rocks and ground water are taken over by the lateral surface of a mine during the entire time required for lining construction. In order to calculate the optimal thickness of an ice wall () and to take into account the rheological behavior of frozen soils, S.S.Vyalov suggested the following relations [10, 11]:     1 111 1 ' 1 ,                     mm m p m p h E a k A t T a (1)                          1 2 1 1 2 , m m p m p E a A t T a (2)   3   s p ph E t (3)                              2 0 1 2 0 45 /2 1 0 45 / 2 1 1 1 2 45 / 2 tg p p tg E a c t tg (4) where  is the optimal ice wall thickness, a is the cylinder inner radius, k’ is the coefficient dependent on soil compaction conditions, m, A(tp, ) are the parameters characterizing the rheological behavior of the soil, tp is the loading period, p is the load applied to the ice wall, h is the height of unfixed part of mine shaft, is the maximum radial movement of the frozen A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 396 soil, s(tp) is the ultimate strength of frozen soil to uniaxial compression, c(tp) is the cohesion, and  is the angle of internal friction. Eqns. (1)-(2) estimate for the cylinder wall thickness according to the criterion for the maximum radial movement of its internal wall, and Eqns. (3)-(4) - the ultimate stress criterion. We notice that Eqns. (1) and (3) are obtained on the assumption that the ice wall is a cylinder of finite height, and Eqns. (2) and (4) – on the assumption that the ice wall is a cylinder of unlimited height. It is also worth noting that formula (3) is obtained based on relationship (1) assuming that the frozen soil exhibits ideal plastic behavior described by the von Mises model. In [12], the applicability of formula (1) for calculating the frozen wall thickness was evaluated. The aim of our work is to analyze the correctness of analytical relationships (3) and (4) for calculating the frozen wall thickness according to the ultimate stress criterion. For this purpose, the ice wall thickness calculated by these formulas was compared with the numerical simulation results obtained by the formulation whose geometry is close to the real shaft sinking conditions. Using the obtained data, we propose a modification of formula (4) because it shows better agreement with the numerical simulation results. In addition, we perform a comparative analysis of the values of ice wall thickness calculated by formulas (1)-(4) and by their modifications, which provides guidelines for using each of these formulas. MATHEMATICAL FORMULATION n order to evaluate the correctness of relations (3) and (4), we have computed the problem of deformation of ice wall exposed to the external radial load p applied to the lateral surface of the cylinder and the vertical load p’ acting on the upper end of the cylinder. The calculation scheme of shaft sinking is given in Fig.1. An ice wall is modeled as a two- piece cylinder. One part is a hollow cylinder of a height equal to that of the ice wall section without lining, and another part is a solid section. The soil inside the ice wall is assumed to be a thawed soil that experiences only elastic deformations. The load applied to the lateral surface of the cylinder is calculated by the formula for mining and hydrostatic pressure, which takes into account the influence of cohesive forces and the results of hydrogeological observations. Figure 1: Computational domain The mathematical model is based on the following hypotheses and assumptions: 1) The computational domain has radial symmetry, and therefore calculations are made in axisymmetric formulation. 2) The rocks under consideration are isotropic at macrolevel. 3) Deformations as well as strain increments at each step are small. 4) The time and temperature effects on the strength and elastic characteristics of rocks are taken into account parametrically. 5) The surrounding thawed and burden formations generate constant rock pressure. 6) A linear relation exists between the stress and elastic strain tensors. 7) The minimum value for the ice wall thickness is determined on the assumption that no plastic deformations are present on its internal wall. 8) The Mohr-Coulomb criterion is used as a yield criterion that is commonly applied to describe soil failure (especially, that of cohesionless soil). By this criterion, the rock failure begins when the maximum tangential stress n reaches the critical value which is dependent on the normal stress n in the same area. The mathematical statement of the problem includes equilibrium Eqn. (5), Hooke’s law (6), geometric relation for small strain tensor (7), yield condition (8), and associated plastic flow rule (9): I A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 397   σ 0 (5)  : σ C ε εp (6)  1 2      ε u u T (7) 0     nnF tg c (8)     ε σ  p F (9) where ={r, ,z, rz} is the stress tensor, ={r, ,z, rz} is the strain tensor, С is the tensor of elastic constants which reduces in the isotropic case to two elastic constants (E′ – Young’s modulus,  – Poisson’s ratio), u={ur, uz} is the displacement vector, p={pr, p,pz, prz} is the plastic strain tensor, F is the Mohr-Coulomb yield criterion, n is the maximum tangential stress in the area with a normal n, n is the normal stress operating in the same area,  is the indefinite multiplier determined using Prager’s compatibility conditions, and the dots over the symbols denote time derivatives. The Prager compatibility conditions are written as 0, 0  F , (or 0F (10) 0, 0   F , (and 0F ) (11) The system of Eqns. (5)-(11) is supplemented with boundary conditions that correspond to the calculation scheme from Fig.1. 1 0 zu (12) 2   n σ P (13) 3 '  n σ P (14) where the vectors P={p,0} and P′={p′,0} correspond to the calculated loads obtained by the formulas:  r hp p p (15) 0 0 2 90 90' 2 2 2                      r mp g h tg c tg (16)   h w gwp g h (17) ' '  mp g h (18) In formulas (15) - (18), the following notation is used: m=2000 kg/m3 – average density of the material, h′ – bed rock density, w=1000 kg/m3 – water density, g=10 m/s2 – gravitational constant, hgw=1.5 m – groundwater depth. The load ph is considered only for saturated rocks. A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 398 MATERIALS AND METHOD he allowable ice wall thickness was determined for three typical materials: sand, chalk and clay. We considered the bed rock depth values of 100m, 200m, 300m, and 500m. Tabs. 1 and 2 give the mechanical parameters of the soil in thawed and frozen states which were provided by the Institute of Nature Management of the NAS of Belarus. The parameters of the frozen soil were obtained at T=-80C and for load duration of 12 h. These data were used for calculating the load acting on the ice wall according to formulas (15) – (18). The obtained data are given in Tabs. 3-5. Soil E′, GPa  с, kPa , 0 Sand 0.32 0.3 9.6 30 Chalk 0.3 0.35 1 31.5 Clay 0.364 0.18 105 25 Table 1: Mechanical parameters for the soil in thawed state. Soil E′, GPa  с, MPa , 0 s, MPa Sand 5.3 0.18 6.33 37 4.4 Chalk 2.5 0.15 6.20 24 7 Clay 1.7 0.17 3.78 10 2.2 Table 2: Mechanical parameters for the soil at T=-80C. Load 100 m 200 m 300 m 500 m p, MPa 1.557 3.222 4.889 8.222 p′, MPa 2 4 6 10 Table 3: Loads acting on the frozen sand ice wall Load 100 m 200 m 300 m 500 m p,MPa 0.626 1.253 1.881 3.135 p′,MPa 2 4 6 10 Table 4: Loads acting on the frozen chalk ice wall. Load 100 m 200 m 300 m 500 m p, MPa 0.678 1.490 2.301 3.925 p′, MPa 2 4 6 10 Table 5: Loads acting on the frozen clay ice wall. The finite-element method was used for the solution of boundary value problem. The examined area was divided into rectangular elements. We controlled the convergence of the numerical solution based on a series of calculations on grids with elements of different sizes for the clay layer at a depth of 200m. The sizes of elements in the ice wall area with no mine working support (lining) varied from 1.28 to 0.02m. The relative error erri is defined by the formula T A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 399 min min 100%       i i m m m err (19) where i[0.02; 1.28] indicates the size of the element used to determine the von Mises stress value m on the internal wall of the ice cylinder, and minm is the von Mises stress value on the internal wall of the ice cylinder obtained in the solution of the problem with minimum element size. Fig. 2 shows the effect of size of the finite element of the grid on relative error. The results demonstrate that the optimal size of the element is 0.04m. It is interesting that if the element size is maximum (1.28m), then the relative error does not exceed 8.3%, which is well within the permissible limits of an engineering error. Figure 2: Relative error vs finite element size. RESULTS AND DISCUSSION Comparison of numerical simulation results and analytical estimates he developed mathematical model is used for assessing the optimal ice wall thickness according to the criterion for ultimate stress-strain state. The calculation results are summarized in Tab. 6. When the optimal ice wall thickness is calculated by Eqn. (2), the ultimate uniaxial compression strength values are taken from Tab. 2. The values for the cohesion and internal friction angle used in Eqn. (4) are also given in Tab. 2. The calculation results are summarized in Tabs. 7 and 8. The obtained results demonstrate that, for all materials under study, the design thicknesses obtained by Eqn. (3) exceed those found using Eqn. (4). This can be attributed to the fact that the internal friction is ignored, which increases the strength. In general, the scatter in the estimates obtained by Eqns. (3) and (4) is great (12 times for sand); especially, when the ice wall experiences high loads. h′, m Sand Chalk Clay 100 1.05 0.35 0.75 200 1.75 0.50 1.2 300 2.35 0.75 1.4 500 2.5 1.05 2.6 Table 6: Ice wall thickness obtained via numerical simulation. h′, m Sand Chalk Clay 100 3.065 0.774 2.667 200 6.342 1.550 5.864 300 9.623 2.327 9.058 500 16.183 3.879 15.451 Table 7: Ice wall thickness calculated by formula (3). T A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 400 h′, m Sand Chalk Clay 100 0.304 0.171 0.404 200 0.595 0.340 0.910 300 0.860 0.508 1.441 500 1.330 0.838 2.582 Table 8: Ice wall thickness calculated by formula (4). A comparison of the ice wall thickness values found during the numerical solution of problem (5)-(18) with the results obtained by Eqns. (3) and (4) is illustrated in Fig.3. For all materials under study, Eqn. (3) shows the most distinction from the numerical simulation results (Fig. 3a). The greatest ice wall thickness reserve is seen for clay and sand, especially, for the layers at a depth of 500 m. This is because the ice wall is subjected to maximum loads at this depth. In this case, the deviations from the numerical results obtained for clay and sand are the greatest, 12.8 m and 13.7m, respectively. (a) (b) Figure 3: Difference between the designed ice wall thickness obtained via numerical simulations and the analytical estimations for the ice wall thickness found by Eqns. (3) (a) and (4) (b). It was found that the results of numerical simulations of the ice wall thickness for chalk and sand are higher than the estimations obtained by Eqn. (4) (Fig.3 b). The most distinction from the numerical simulation results is 1.49m for sand at a depth of 300 m and 0.242 for chalk at a depth of 300 m. A comparison of the same results but obtained for clay shows that the numerical calculations give the higher ice wall thickness values for the material at depths of 100-200 m. At depths of 300-500 m, the simulation results are comparable with those obtained by formula (4); the relative error does not exceed 2.8%. Modification of Eqn. (4) Analysis of the obtained results has revealed that Eqn. (4) has the best correlation with the numerical simulation results. This equation was modified here in order to get a more accurate description of the relation between the designed ice wall thickness and the applied load. The modification made it possible to minimize the difference between the calculation results obtained by this formula and the numerical simulation data. Thus, we have                                   2 0 1 2 0 45 /2 1 0 45 / 2 1 1 1 2 45 / 2 tgp tg E k a d c tg (20) where k and d are the piecewise constants 1 1 2 2 2 3   , ,       k p p p k k p p p A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 401 1 1 2 2 2 3   , ,       d p p p d d p p p The intervals for the coefficients for each material were set reasoning from the best agreement with the numerical simulation data. Approximations of the parameters k and d for each material are given in Tab. 9 and 10. The calculated results for the designed ice wall thickness for chalk, clay and sand obtained by new Eqn. (20) are presented in Fig.4. Soil k1 k2 p1, MPa p2, MPa p3, MPa Sand 2.337 0.319 1.557 4.889 8.222 Chalk 1.186 0.909 0.626 1.881 3.135 Clay 0.624 1.052 0.678 2.301 3.925 Table 9: Approximation of the coefficient k. Soil d1 d2 p1, MPa p2, MPa p3, MPa Sand -6.674 5.648 1.557 4.889 8.222 Chalk -0.847 0.760 0.626 1.881 3.135 Clay 2.516 -0.387 0.678 2.301 3.925 Table 10: Approximation of the coefficient d. (a) (b) (c) Figure 4: Approximation of Eqn. (4) (solid lines) by the piecewise-linear functions (dashed lines) for sand (a), chalk (b) and clay (c). A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 402 In addition to approximation (20), it is suggested that formula (4) can be modified as                                   2 0 1 2 0 45 /2 1 0 45 / 2 1 1 1 2 45 / 2 tgp tg E a z c tg (21) The values of parameter z for sand, chalk and clay are given in Tab. 11. The calculated results are obtained for the designed thickness by Eqn. (21) and the numerical simulation results are shown for each material in Fig.5. The approximations are nonlinear, which is caused by the nonlinear dependence of stress on the radial coordinate and, as a consequence, on the desired ice wall thickness. The nonlinear character is determined by the material parameters of the soil, e.g., for chalk and sand, the curve is convex upwards and for clay – downwards. Chalk and sand exhibit similar cohesion values (6.2 and 6.33 MPa), and clay has significantly lower cohesion value (3.78 MPa). It is interesting that almost the same nonlinear behavior is observed for Eqn. (4). Sand Chalk Clay Z 1.140 0.198 0.153 Table 11: Values of the parameter z, entering Eqn. (21). (a) (b) (c) Figure 5: Results of approximation of Eqn. (4) (curve 2) by function (21) (curve 1) for sand (a), chalk (b), clay (c). Markers – numerical solution. A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 403 A comparison of the results of approximation by Eqn. (4) shows that in case of piecewise-linear dependence (20) they are in rather good qualitative and quantitative agreement. A deficiency of this approach is the necessity of calculating four material parameters, as well as the intervals for determining piecewise-linear correction coefficients. Therefore, it is advisable to use formula (21) in order make rapid estimation for the designed ice wall thickness. In this case, it is necessary to determine only one material parameter, whereas formula (20) can be used when there is a need for more accurate calculation. Comparative analysis of Vyalov’s formulas for ultimate displacements and stresses Fig. 6 presents the calculation results for ice wall thickness obtained for sand, chalk and clay by Eqns. (1)-(2) recommended for ultimate displacement calculations and by Eqns. (3)-(4) - for strength calculations. Apart from these formulas, this figure also shows the calculation results obtained by Eqn. (20) and by the equation which is a modification of relation (1) [12]:                         1 111 1 ( ) 1 2 , mm m p m p ha E g p A t T a (22) where the material function g(p) is dependent on the value of the load acting on the ice wall and determined by applying the numerical simulation results. The function for the materials under study was taken from [12]. (a) (b) (c) Figure 6: Dependence of the designed ice wall thickness on the load for sand (a), chalk (b), clay (c): curve 1 – Eqn. (1), curve 2 – Eqn. (2), curve 3 – Eqn. (3), curve 4 – Eqn. (1), curve 5 – Eqn. (20), curve 6 – Eqn. (22). For sand (Fig.6a), the maximum discrepancy between Eqns. (4) and (20) is 1.49m, which corresponds to the load of 889 МPа. Analogous thicknesses are determined using Eqns. (3) and (22) for the load less than 4.889 МPа. Besides, for the load not exceeding 3.222 МPа, the thickness calculated by Eqn. (2) is close to the value obtained by Eqns. (3) and (22). Eqn. (2) A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 404 cannot be used in this case to calculate the ice wall thickness at a depth of 500 m because the base of power function becomes negative. For chalk Fig.6(b), there are two similar data blocks. The first group includes Eqns. (4) and (20), and another – Eqns. (2), (3) and (22). Eqn. (1) yields the maximum ice wall thickness reserve at all loads examined here. The maximum difference between Eqns. (4) and (20) is equal to 0.242 m and is achieved at the load of 1.881 МPа. The ice wall thicknesses which enter the second group have similar values up to the load of 1.881 МPа. At loads more than 2.25 МPa, Eqn. (2) gives values that exceed those calculated by Eqns. (3) and (22). Analogous results are obtained for clay (Fig. 6c). Two groups of equations provide almost similar values for ice wall thickness: (4), (20) and (2), (22). The maximum difference between Eqns. (4) and (20) corresponds to the load of 0.678 МPа and is equal to 0.35m. Eqns. (2) and (22) yield close thickness values at the load less than 2.301 MPa , and Eqn. (2) cannot be used to calculate the ice wall thickness at the load exceeding this value. The maximum margin of safety is obtained using Eqn. (3). On the basis of the obtained results, some conclusions were drawn regarding the applicability of the analyzed relations to ice wall thickness calculations. To perform strength calculations, it is desirable that Eqn. (20) would be used for all soils under study because Eqn. (3) yields excess thickness, and Eqn. (4) is obtained on the assumption of the ultimate state of the ice wall, and therefore it should be used with some prescribed margin of safety. In ultimate displacement calculations for chalk, it is advisable to use Eqn. (22) because Eqn. (1) yields the excess thickness for all loads under consideration. At loads less than 2.25 MPa (this value corresponds to a depth of 300 m), Eqn. (2) can be applied. In ultimate displacements calculations for clay up to the load of 2.301 МPа (this value corresponds to a depth of 300m), relations (2) and (22) can be used. In ice wall thickness calculations for frozen sand at the load no less than 4.889 МPа (this value corresponds to a depth of 300 m), it is desirable that relation (22) is used. At loads less than 3.222 MPa (this value corresponds to a depth of 200 m), Eqn. (2) can also be applied. At high loads that exceed 4.899 for sand and 2.301 MPa for clay (these values correspond to a depth of 300 m), the ice wall thicknesses calculated by Eqns. (1)-(4), (20) and (22) exhibit great scatter, which generates a need for further investigations aimed at modifying the existing constitutive equations. This necessity stems from the fact that at great loads the hydrostatic pressure can have a significant effect on the stress-strain state of rocks. CONCLUSION e have performed a theoretical study to evaluate the applicability of analytical Vyalov’s formulas to the calculation of the ultimate stress state of an ice-rock cylinder of unlimited and finite heights. The numerical results showed that the design values of the ice wall thickness obtained by the formula for a cylinder of limited height have a margin of safety for all considered rocks (sand, chalk, and clay). This is associated with the fact that the internal friction is ignored, which increases the strength. The application of the Coulomb–Mohr criterion made it possible to model the deformation of the ice-rock cylinder. A comparison of the thicknesses obtained numerically with the results found by Eqn. (3) showed that the greatest thickness reserve was observed for clay and sand. Using the results of numerical simulation of the stress-strain state of the ice wall, we suggested two modifications of Eqn. (4). The first variant involves the use of the piecewise-linear relation and enables describing the simulation results both qualitatively and quantitatively. The second is a simple variant and it can be used for making rapid estimation for the designed ice-wall thickness. A comprehensive comparative analysis of six formulas (four Vyalov’s formulas and two modifications) has revealed that in strength calculations it is advisable to use Eqn. (20) for all soils examined here because Eqn. (3) causes the excess margin of thickness to occur, and Eqn. (4) must be used with some prescribed margin of safety. In ultimate displacement calculations for chalk, Eqn. (22) or (2) should be used at the load not exceeding 2.25 MPa, which corresponds to the soil depth of 300 m. In creep calculations for clay at the load of 2.301 MPa, which corresponds to the soil depth of 300 m, relations (2) and (22) can be used. In order to calculate the ice-wall made of frozen sand at the load less than 4.889 MPa (up to the soil depth of 300 m), it is advisable to use relation (22) and relation (2) at the load less than 3.222 MPa (up to the soil depth of 200 m). At high loads that exceed 4.899 MPa for sand and 2.301 MPa for clay there is a need for further investigations so that the existing constitutive equations describing the stress-strain state of rocks can be modified. W A. Kostina et alii, Frattura ed Integrità Strutturale, 53 (2020) 394-405; DOI: 10.3221/IGF-ESIS.53.30 405 ACKNOWLEDGMENTS his research was supported by 17-11-01204 project (Russian Science Foundation). REFERENCES [1] Levin, L. Yu, Semin, M.A. and Plekhov, O.A. (2018). Comparative analysis of existing methods for calculation frozen wall thickness for mine shafts under construction, Bulletin of PNRPU. Construction and Architecture, 9(4), pp. 93-103. DOI: 10.15593/2223-9826/2018.4.09. [2] Lai, Y., Xu, X., Dong, Y. and Li, S. (2013). Present situation and prospect of mechanical research on frozen soils in China, Cold Reg. Sci. Technol., 87, pp. 6-18. DOI: 10.1016/j.coldregions.2012.12.001 [3] Liu, X., Liu, E., Zhang, D., Zhang, G. and Song, B. (2019). Study on strength criterion for frozen soil. Cold Reg. Sci. Technol.,161, pp. 1-20. DOI: 10.1016/j.coldregions.2019.02.009. [4] Nishimura, S., Gens, A., Olivella, S. and Jardine, R. J. (2009). THM-coupled finite element analysis of frozen soil: formulation and application, Geotechnique, 59(3), p. 159. DOI: 10.1680/geot.2009.59.3.159. [5] Ghoreishian, A. S. A., Grimstad, G., Kadivar, M. and Nordal, S. (2016). Constitutive model for rate-independent behavior of saturated frozen soils. Can. Geotech. J., 53(10), pp. 1646-1657. DOI: 10.1139/cgj-2015-0467. [6] Na, S. H. and Sun, W. C. (2017). Computational thermo-hydro-mechanics for multiphase freezing and thawing porous media in the finite deformation range. Comput. Methods. Appl. Mech. Eng.,318, pp. 667-700. DOI: 10.1016/j.cma.2017.01.028. [7] Liu, E., Lai, Y., Wong, H. and Feng, J. (2018). An elastoplastic model for saturated freezing soils based on thermo- poromechanics. Int. J. Plasticity,107, pp. 246-285. DOI: 10.1016/j.ijplas.2018.04.007. [8] Alonso, E. E., Gens, A. and Josa, A. (1990). A constitutive model for partially saturated soils, Géotechnique, 40(3), pp. 405-430. DOI: 10.1680/geot.1990.40.3.405. [9] Liu, E. L. and Xing H. L. (2009). A double hardening thermo-mechanical constitutive model for overconsolidated clays. Acta Geotech., 4(1), pp. 1-6. DOI: 10.1007/s11440-008-0053-4. [10] Vyalov, S.S., Zaretsky, Y.K., Gorodetsky, S.E. (1979). Stability of mine workings in frozen soils, Engineering Geology, 13 (1-4), pp. 339-351. DOI: 10.1016/0013-7952(79)90041-3. [11] Vyalov, S.S. (1986). Rheological fundamentals of soil mechanics. Elsevier, Amsterdam, the Netherlands [12] Zhelnin, M., Kostina, A., Plekhov, O., Panteleev, I. and Levin, L. (2019). Frattura ed Integrità Strutturale, 49, pp. 156- 166. DOI: 10.3221/IGF-ESIS.49.17. T << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice