Microsoft Word - numero_49_art_17_2462 M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 156 Focused on Russian mechanics contributions for Structural Integrity Numerical analysis of application limits of Vyalov’s formula for an ice-soil wall thickness M. Zhelnin, A. Kostina, O. Plekhov, I. Panteleev Institute of Continuous Media Mechanics of the Ural Branch of Russian Academy of Science, Russia zhelnin.m@icmm.ru, http://orcid.org/0000-0003-4498-450X kostina@icmm.ru, http://orcid.org/0000-0002-5721-3301 poa@icmm.ru, http://orcid.org/0000-0002-0378-8249 pia@icmm.ru, http://orcid.org/0000-0002-7430-3667 L. Levin Mining institute of the Ural Branch of Russian Academy of Science, Russia aerolog_lev@mail.ru ABSTRACT. Artificial ground freezing is an efficient technique which allows one to sink the mine shafts under complex hydrogeological conditions. The aim of artificial ground freezing is a creation of a temporary wall of frozen soil around the intended excavation. To estimate an optimal thickness of an ice-soil wall, the Vyalov’s formula is widely used by engineers. The article is devoted to analysis of Vyalov’s formula on the basis of the numerical simulation. The numerical simulation has been conducted by the finite element method. For the simulation of a stress-strain state of an ice-soil wall a new computational scheme has been proposed. The scheme is based on Vyalov’s design layout for a vertical shaft sinking and takes into account a soil layer beyond the excavation bottom. A mechanical behavior of frozen soil is described by Vyalov’s constitutive relations. As a result, it has been shown that values of the wall thickness given by Vyalov’s formula do not agree with the ones obtained by the numerical simulation. In order to conform results given by Vyalov’s formula and the numerical simulation, two modifications of the formula have been proposed. KEYWORDS. Artificial ground freezing; Mine shaft; Ice-soil wall: Inelastic deformation; Numerical simulation. Citation: Zhelnin, M., Kostina, A., Plekhov, O., Panteleev, I., Levin, L., Numerical analysis of application limits of Vyalov’s formula for an ice-soil wall thickness, Frattura ed Integrità Strutturale, 49 (2019) 156-166. Received: 04.04.2019 Accepted: 22.05.2019 Published: 01.07.2019 Copyright: © 2019 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 t the present time mineral deposits characterize of an increase of a depth of an occurrence and a complication of hydrogeological conditions of building of a vertical mine shaft for their elaboration. Artificial ground freezing (AGF) is an efficient technique which allows sinking of the mine shafts through weak, unstable, water saturation A M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 157 soil stratums. The aim of AGF is a creation of a temporary wall of frozen soil around the intended excavation. Ice in soil pores changes its mechanical and hydrodynamical properties, so at certain temperature the ice-soil wall can withstand rock pressure from unfrozen rock mass and eliminate groundwater filtration. The freezing of ground is carried out by a series of wells that are drilled around the excavation. Due to circulation of liquid refrigerant inside the pipes, heat extracts from the near-wellbore domain and the ice-soil retaining structure is established. Strength and reliability of the ice-soil wall determine safety of a shaft sinking and geotechnical works. Mechanical behavior of the wall depends on properties of soil that its formed. Presence of ice and unfrozen water in frozen soil leads to occurrence of rheological properties that characterize development of slow deformations (creep) and losing strength during a prolonged load action [1]. As a result, the ice-soil wall could significantly deform under rock pressure [2]. A lot of experimental studies have been carried out to study time-dependent deformation of frozen soil under uniaxial and multiaxial compression tests [1-10]. Creep deformation of frozen soil depends on long-term strength limit [3]. If the stress does not exceed the long-term strength limit, deformation rate tends to zero and failure time is infinity. In the case evolution of deformation is characterized as attenuate creep. If the stress is larger than the long-term strength, non- attenuate creep develops. In general, the creep process can be divided on three stage: unsteady creep with a decreasing deformation rate, viscoplastic flow with an almost constant rate, and progressive flow with an increasing rate after that failure arises. In [1, 3] it is noted that the long-term strength limit of frozen soil significantly depends on temperature of frozen soil. In [4] from uniaxial creep tests of warm ice-rich sand it has been established that creep characteristics are affected by small temperature variation. In [5, 6] on the basis of triaxial tests it has been shown that a decrease of temperature leads to an improvement of strength properties. In [7] a study of an influence of thermal gradient on the frozen sand has been carried out. It has been determined that a decrease of thermal gradient leads to an increase of the long-term strength limit. In [8] it is noted that rheological properties of frozen soils also depend on ice and unfrozen water saturations, grain-size distribution, mineral content. In [9] it has been concluded that cryostructure has an influence on creep behavior of frozen soil. In [10] based on triaxial creep tests it has been shown that coarse-grained content in frozen silty clay causes a rise of it strength. For describing the phenomena of creep of frozen soils various constitutive relations have been developed. To reflect changes in soils that arises during creep process at the microscale into the macroscale, micromechanics [11], thermodynamics [12, 13] and damage mechanics [10] are used. However, to determine material parameters for this models difficult and costly experimental studies have to be conducted. As a result, for engineering purposes macroscale phenomenological models are widely used [14]. In [3, 14, 15] extensive reviews of various models describing creep under complex stress state and various stages of deformation are presented. It could be divided in two types: models that take into account an influence of confining pressure on creep deformation and models in that the influence is supposed to be negligible. In [16, 17] Nishihara and Burgers models have been improved to a new creep constitutive models that describe viscoplasticity deformation of frozen soils under high confining pressure. One of basic model for analysis of creep deformation of ice-soil retaining structure is Vyalov’s constitutive relations [1,2]. The relations are based on the Norton-Bailey power law in that a modulus depends on time and temperature. An influence of confining pressure on creep of frozen soil is neglected. The advantage of the model is that material parameters for describing deformation in complex stress state can be estimated from simple uniaxial loading tests. Experimental verification of the Vyalov’s relations has been performed in uniaxial and multiaxial creep tests in wide range of temperature and ice saturation of frozen soils [1, 2, 18]. In [19] experimental values of creep deformation of beams of frozen soil have been compared to values obtained by a finite element modeling. In [20] experimental and numerical studies of interaction of footings with frozen sand have been performed. In [21, 22] an improvement of Vyalov’s relations has been performed by consideration of an influence of volumetric ice content on creep deformation. In [1, 23] on the basis of Vyalov’s relations an analytical formula for an estimation of an optimal thickness of an ice-soil wall from a condition of maximal creep deformation has been derived. Nowadays Vyalov’s formula is widely applied by engineers for design parameters of an ice-soil retaining structure and determination of AGF regimes [24]. However, analytical derivation of the formula is accompanied by a series of strong assumptions that could not agree with real conditions of a shaft sinking. In [25] deformation of an ice-soil wall of Callovian sandy loam has been estimated by numerical modelling with using Vyalov’s relations for various depths of a shaft sinking. It has been shown that deformation of an ice-soil wall with thickness estimated by Vyalov’s formula exceeds admissible values for depths of the sinking more than 100 meters. The present work is devoted to analysis of Vyalov’s formula on the basis of numerical modeling. The simulation of a sinking of a vertical shaft under AGF is performed according to Vyalov’s design layout of a shaft sinking [1] modified by consideration of a soil layer beyond the excavation bottom. It is considered three type of soils: chalk, sand and clay that M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 158 are prevalent on potash deposits in the Petrikov Site in Belarus. Results of the numerical simulation is used for modifications of Vyalov’s formula to reconcile values of the wall thickness obtained by the formula to the numerical ones. MATHEMATICAL MODELS FOR DEFORMATION OF AN ICE-SOIL RETAINING STRUCTURE n [1] a design layout for a vertical shaft sinking with using AGF has been proposed Fig. 1. The ice-soil retaining structure is a hollow cylinder of inner radius a and outer radius b that encloses unfrozen soil. As a shaft is sunk an inner surface of the cylinder is exposed. The shaft sinking is performed gradually with an installation of a shaft tubbing lining. Therefore, only a height h of the inner surface is unsupported. This part of the surface is deformed under rock pressure p during the time prt that requests for building of the lining. At that the upper end of the cylinder is supposed to be fixed by the tubbing lining, as it prevents of deformation of the soil to the cylinder center. Figure 1: Vyalov’s design layout for a vertical shaft sinking with applied of AGF. According to the scheme in [16] the following system of equation in axial-symmetric configuration for describing a stress- state state of an ice-soil wall has been proposed:            0 rzr r r r z , (1)      0 rz rz z r , (2) , , , , zr r r r z rz uu u u r r z z              (3) with boundary conditions          0r rz zr a r a z h , (4)    r r b p , (5)     0 0z zz z h u u , (6)   0r z hu , (7) where ( , , )r z – cylindrical coordinates, u – displacement vector, σ - Cauchy stress tensor, ε – strain tensor, a , b – inner radius and outer radius of the ice-soil cylinder. The geometry of the ice-soil cylinder is presented in Fig 1. As it can I M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 159 be seen from (1)-(3) axial stress is supposed to be constant along the height of the cylinder and axial strain is supposed to be constant along the radius of the cylinder. The Vyalov’s constitutive relations for frozen soil is generally accepted. The relationship between stress and strain are written as   ( , ) meff effA T t , (8) where   2 dev( ):dev( ) 3 eff ε ε is the equivalent strain,   3 dev( ):dev( ) 2 eff σ σ is the equivalent stress, dev – deviatoric part of a tensor,      1 2( , ) 3 ( 1) m A T t T t – deformation coefficient depending on temperature T of the frozen soil and time t of loading action, m ,  ,  – material constants. It is assumed that frozen soil is isotropic, homogenous material. The volumetric creep strain of frozen soil is taken to be zero: tr 0ε , (9) where tr – trace of a tensor. The hydrostatic stress trσ has no an influence on the strain. A contribution of the elastic strain is neglected. An additional Hencky’s relations are taken into the account: tr     σ( ), , ,i i i r z (10) 2rz rz  , (11) From engineering practice and experiments on a physical modeling of an operation of an ice-soil wall it is well known that under an action of rock-pressure radial strain on the inner surface of the cylinder could attain significant values that threaten of failure of freezing pipes and deviation on designed characteristics of the mine shaft. Therefore, the problem of estimation of the ice-soil wall thickness E b a  under that the radial displacement ru on the inner surface do not exceed a limiting value  is stated:   r r au . (12) On the basis of analytical solution of the problem (1)-(11) taking into account the condition (12) in [16] the following estimation of an optimal thickness E has been obtained: 1 1 1(1 ) 1 1 ( , ) m m m pr m ph E a A T t a                . (13) The formula (13) is named Vyalov’s formula. To analytically solve the problem (1)-(11) a series of assumptions have been purposed. The key assumption relates to expression for distribution of shear strain rz on the upper end of the ice-soil cylinder. According to the assumption a distribution of the shear strain is approximated by the following expression: rz a rh    . (14) The assumption allows significantly simplifying mathematical treatment and obtaining formula convenient for engineering calculations. Also it should be noted that the problem (1)-(11) has been solved for small neighborhoods of the ends of the ice-soil cylinder. M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 160 In order to analysis Vyalov’s formula a numerical simulation of stress-strain state of the ice-soil cylinder has been carried out. Vyalov’s design layout has been modified by taking into account a soil layer beyond the bottom of the excavation as usual it performs in modelling of vertical shafts (Fig. 2). The layer also consists of the hollow ice-soil cylinder and unfrozen soil inside it. The system of equations considered in the simulation is written as: div 0σ , (15)   T 1 [grad ( grad ) ] 2 ε u u , (16) where grad – the gradient operator, div – the divergence operator. Boundary conditions are identical to (5), (6). The other boundaries of the cylinder are free. Figure 2: A scheme of an ice-soil wall that is used in the numerical simulation. The constitutive relations for describing inelastic strain of the ice-soil cylinder under the rock pressure coincide with Vyalov’s relations. An additional it is taken into account elastic part of strain. The unfrozen soil it is supposed to be isotropic, homogenous material that can undergo only elastic deformation. The simulation is performed in axial-symmetric configuration with using the finite element method. The geometrical domain is divided on mapped elements of the second order. The rock pressure p acting on the ice-soil retaining structure is estimated by the standard engineering formula [26]:  l wp p p , (17)                  2Tan 2 Tan 2 2 l sp Hg C , (18) w w wp gH , (19) where lp – effective lateral pressure, wp – hydrostatic pressure, H – depth of bottom of the soil stratum within that the ice-soil cylinder is placed, l – average density of soils in the rock mass where a vertical shaft is constructed,  – friction angle of unfrozen soil, C – cohesion coefficient of unfrozen soil, w – density of groundwater, wH – groundwater hydraulic head. RESULTS OF THE NUMERICAL SIMULATION he numerical simulation has been conducted for three types of soils: clay, sand and chalk. Elastic and rheological properties of the soils have been determined from uniaxial compression tests of cylindrical samples of core from the Petrikov Site of Starobinsky potash deposit in Belarus. This region is elaborated by Belaruskali. The T M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 161 determined mechanical properties for frozen soils are listed in Tab. 1 and for unfrozen soils in Tab. 2. According to design documentations for building of a vertical mine shaft accepted by Belaruskali, minimal temperature of frozen soil in an ice-soil wall must not exceed 8T C   , so the properties are determined at this temperature. The deformation coefficient is estimated for 24prt hour. Technical parameters for a vertical shaft sinking with using AGF given by the documentation are presented in Tab. 3. The numerical simulation has been carried out for the depth H of soil stratums occurrence within the range from 100 m to 500 m. Soil A , MPa K , GPa G , GPa m Clay 9.37 0.86 0.73 0.47 Sand 6.65 2.76 2.25 0.3 Chalk 15.90 1.19 1.09 0.48 Table 1: Mechanical properties of the frozen soils at temperature 8  T C . A – deformation coefficient in the constitutive relation (8), K – bulk modulus, G – shear modulus, m – degree in in the constitutive relation (8). Soil K , GPa G , GPa C , kPa  ,  Clay 0.19 0.15 105 25 Sand 0.27 0.12 9.6 30 Chalk 0.33 0.11 1 31.5 Table 2: Mechanical properties of the unfrozen soils. K – bulk modulus, G – shear modulus, C – cohesion coefficient,  – friction angle. a , m h , m  , m 5.25 5 0.1 Table 3: Technical parameters for a vertical shaft sinking. In Fig. 3 - 5 distributions of radial displacement ru in an excavation obtained by the numerical simulation are presented. The thickness S of the ice-soil wall is estimated by Vyalov’s formula (13). For the depth up to 300 m Vyalov’s formula gives an overvalued estimation of optimal thickness of the ice-soil wall independent on considered soil. In chalk stratum the obtained displacements do not exceed the admissible value for all considered range of the depths. At the depth 500 m the radial displacement of the ice-soil wall of frozen clay exceeds the admissible value by 1 cm. At the same time the displacement of the wall of frozen sand at this depth exceeds the admissible value more than two times. Qualitatively the distributions for soils under consideration coincide. A maximum value of the displacement attains on the inner surface of the ice-soil cylinder. Thus it can be concluded that Vyalov’s formula gives an overvalued estimation of the optimal thickness of the wall of quasi-brittle carbonate rocks such as chalk. For more ductile soils such as clay and sand Vyalov’s formula gives an overvalued estimation of the optimal thickness of the wall for the depths less than 300 m and an undervalued the one for the large depths. Distributions of shear strain rz in depending on radius on the top end of the ice-soil cylinder are presented in Fig 6. It can be seen that for soils under consideration values of the strain obtained by the numerical simulation significantly differ from the ones given by the Eqn. (14). Normalization of the distribution on the maximum value it allows one to conclude that there is a qualitative discrepancy between the numerical and analytical distributions (Fig 7). M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 162 a b Figure 3: Distribution of radial displacement ru in excavation in a clay stratum at depth H 300 m (a) and 500 m (b). a b Figure 4: Distribution of radial displacement ru in excavation in a sand stratum at depth H 300 m (a) and 500 m (b). a b Figure 5: Distributions of radial displacement ru in excavation in a chalk stratum at depth H 300 m (a) and 500 m (b). a b c Figure 6: Distributions of shear strain rz in depending on radius on the top end of the ice-soil cylinder of frozen clay (a), sand (b), chalk (c). Black curve is given by the expression (14). Other curves are obtained by the numerical simulation. M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 163 a b c Figure 7: Normalized distributions of shear strain rz in depending on radius on the top end of the ice-soil cylinder of frozen clay (a), sand (b), chalk (c). Black line is given by the expression (14), other lines are obtained by the numerical simulation. MODIFICATIONS OF VYALOV’S FORMULA he results of the numerical simulation presented in the previous section show that for the thickness of the ice-soil wall given by Vyalov’s formula the radial displacement either is less than the admissible value or exceeds it. Thus to reconcile Vyalov’s formula to the results of the numerical simulation it has to be modified. From the plots of the shear strain obtained by the numerical simulation (Fig. 6, Fig. 7) it can be seen that the all curves converge to a limit curve under the increasing of the depth. As a consequence, it can be concluded that at the large depth H the shear strain rz on the top end of the ice-soil cylinder do not reflect a relationship between a maximal value of the radial displacement ru of the inner surface of the cylinder and rock pressure p . Thus a modification of Eqn. (14) does not allow one to modify Vyalov’s formula that it coincide with the result of the numerical modeling. Fig. 8 shows estimations of the thickness E of the ice-soil wall obtained by Vyalov’s formula and the numerical simulation on the basis of the condition (12) for clay, sand and chalk stratums. The rock pressure p corresponds to the depth H within the range from 100 m to 500 m. Figure 8: Estimations of the thickness E of the ice-soil wall obtained by Vyalov’s formula (dashed curve) and the numerical modeling (solid line) based on the condition of limiting radial strain for clay, sand and chalk stratums. By analysis of the presented data it has been noted that for a certain value of the rock pressure Vyalov’s formula can be conformed to the results of the numerical simulation by a simple way: 1 2 VyalovE E  , (20) T M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 164 where VyalovE – a value of the thickness given by Vyalov’s formula (13), 1E – a modified value. In Fig. 9 (a) for the soils under consideration values obtained with using the Eqn. (20) is compared to the ones given by the numerical modeling. The presented plots show that independently on soil under consideration the Eqn. (20) can be applied when VyalovE is less than 11 m. For the larger values of VyalovE the Eqn. (20) gives an undervalued estimation of the wall thickness in comparison to the results of the numerical simulation. a b Figure 9: (a) Estimations of the thickness E of the ice-soil wall given by the Eqn. (20) (solid line) and the results of the numerical simulation (markers; the abscissa of a marker is a value given by Vyalov’s formula, the ordinate is a value obtained by the numerical simulation.). (b) Estimations of the thickness E of the ice-soil wall given by the Eqn. (20) (solid line) and the results of the numerical modeling (markers). A value VyalovE depends on a value of the rock pressure acting on the ice-soil wall. Therefore, it can be deduced that the Eqn. (20) is correct within a range of the rock pressure. In order to describe the dependency of the thickness of the ice- soil wall obtained by the numerical simulation within the all range of rock pressure the following modification of Vyalov’s formula has been proposed:                1 1 1 2 (1 )1 1 ( ) 1 2 ( , ) m m m pr m ph E a g p A T t a . (21) Here the function  ( )g g p is a quadratic function 2 2 1 0( )g p p p     , (22) where coefficients i , ( 1, 2, 3)i  , are depended on considered frozen soil. In Tab. 4 values of the coefficients for clay, sand and chalk are listed. Fig. 9 (b) shows estimations of the thickness of the ice-soil wall of frozen clay, sand and chalk given by the Eqn. (21). It can be seen that the proposed equation allows one to describe the results of the numerical simulation both qualitatively and quantitatively. Soil 131 10  , Pa-2 7 2 10  , Pa-1 3 Clay 1.95 -6.02 1.31 Sand 0.51 -2.16 1.20 Chalk 1.00 -4.75 1.42 Table 4: Material coefficients  ,i ( 1, 2, 3)i  , for the Eqn. (22). M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 165 CONCLUSIONS his work is devoted to analysis of Vyalov’s formula for determination of an optimal thickness of an ice-soil retaining structure on the basis of a numerical simulation. For the numerical simulation of a stress-strain state of the ice-soil wall a new computational scheme has been proposed. The scheme is based on Vyalov’s design layout for a vertical shaft sinking and takes into account a soil layer beyond the excavation bottom. The layer consists of a hollow ice-soil cylinder and unfrozen soil inside it. A mechanical behavior of frozen soil is described in according to Vyalov’s constitutive relations. In comparison with an analytical approach used for derivation of the formula in the numerical simulation it has not to make additional assumptions concerning to distributions of stress and strain. As a result of the numerical simulation it has been shown that if a depth of soil stratum occurrence is less than 300 m, Vyalov’s formula gives an overvalued estimation of the wall thickness for the all soils under consideration. For chalk, radial displacement of the wall is less than an admissible value within all range of the studied depths. For sand and clay, the displacement exceeds the admissible value if the depth is large than 300 m. Thus, in quasi-brittle carbonate rock stratum an ice-soil retaining structure with a wall thickness estimated by Vyalov’s formula has excess margin of safety. In ductile soil stratums, that are occurrence at large depth, Vyalov’s formula gives an understated estimation of the wall thickness. It has been studied of an analytical expression for distribution of shear strain on the top end of the cylinder. For soils and depths under consideration the distributions given by the expression do not coincide qualitatively and quantitatively with the ones obtained by the numerical simulation. Since Vyalov’s formula is derived on the basis of the analytical expression, this is a reason due that estimates of the wall thickness obtained by Vyalov’s formula do not agree to the numerical results. In order to conform estimations of the wall thickness given by the Vyalov’s formula to the results of the numerical simulations, two modifications of the formula have been proposed. The first modification could be useful for engineering computations. According to the modification if a value of a thickness of the ice-soil wall given by Vyalov’s formula is less than 11 m than it has to be divided on 2 to coincide with the value given by the numerical simulation. This rule is fulfilled for the all studied soils. The second modification intends including in Vyalov's formula a quadratic function on a rock pressure that acts on an ice-soil wall. Coefficients of the function depends on soil from that consists of the ice-soil wall. The modification allows one to describe the results of the numerical simulation qualitatively and quantitatively within all range of the rock pressure. However, to determine the coefficients of the function the numerical simulations has to be conducted for typical stratums of an elaborated deposit. Thus, it can be concluded that for soil stratums occurrence at small depths a use of Vyalov’s formula in design computations causes overvalued costs on the AGF process. To prevent dangerous breakdown at large depths of a shaft sinking, estimates given by Vyalov’s formula should be corrected by the numerical simulation. ACKNOWLEDGMENTS his research was supported by 17-11-01204 project (Russian Science Foundation). REFERENCES [1] 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. [2] Vyalov, S.S. (1986). Rheological fundamentals of soil mechanics. Elsevier, Amsterdam, the Netherlands [3] Lai, Y., Xu, X., Dong, Y., Li, S. (2013). Present situation and prospect of mechanical research on frozen soils in China, Cold Regions Science and Technology, 87, pp. 6-18. DOI: 10.1016/j.coldregions.2012.12.001. [4] Yang, Y., Lai, Y., Chang, X. (2010). Experimental and theoretical studies on the creep behavior of warm ice-rich frozen sand. Cold Regions Science and Technology, 63(1-2), 61-67. DOI: 10.1016/j.coldregions.2010.04.011. [5] Zhou Z., Ma W., Zhang S., Du H., Mu Y., Li G. (2016). Multiaxial creep of frozen loess, Mechanics of Materials, 95, pp. 172-191. DOI: 10.1016/j.mechmat.2015.11.020. T T M. Zhelnin et alii, Frattura ed Integrità Strutturale, 49 (2019) 156-166; DOI: 10.3221/IGF-ESIS.49.17 166 [6] Esmaeili-Falak M., Katebi H., Javadi A. (2018). Experimental Study of the Mechanical Behavior of Frozen Soils-A Case Study of Tabriz Subway, Periodica Polytechnica Civil Engineering, 62(1), pp. 117-125. DOI: 10.3311/PPci.10960. [7] Zhao, X., Zhou, G. (2013). Experimental study on the creep behavior of frozen clay with thermal gradient. Cold Regions Science and Technology, 86, 127-132. DOI: 10.1016/j.coldregions.2012.10.012. [8] Wang D.Y., Ma W., Wen Z., Chang X.X. (2008). Study on strength of artificially frozen soils in deep alluvium, Tunnelling and underground space technology, 23(4), pp. 381-388. DOI: 10.1016/j.tust.2007.06.010. [9] Bray, M. T. (2012). The influence of cryostructure on the creep behavior of ice-rich permafrost. Cold regions science and technology, 79, 43-52. DOI: 10.1016/j.coldregions.2012.04.003. [10] Hou, F., Lai, Y., Liu, E., Luo, H., & Liu, X. (2018). A creep constitutive model for frozen soils with different contents of coarse grains. Cold Regions Science and Technology, 145, 119-126. DOI: 10.1016/j.coldregions.2017.10.013. [11] Lackner R., Pichler C., Kloiber A. (2008). Artificial ground freezing of fully saturated soil: viscoelastic behavior, Journal of engineering mechanics, 134(1), pp. 1-11. DOI: 10.1061/(ASCE)0733-9399(2008)134:1(1). [12] Andersland O.B., Akili W. (1967). Stress effect on creep rates of a frozen clay soil. Geotechnique 17(1), pp. 27–39. DOI: 10.1680/geot.1967.17.1.27. [13] Fish A.M. (1984). Thermodynamic model of creep at constant stress and constant strain rate, Cold Regions Science and Technology, 9(2), pp. 143-161. DOI: 10.1016/0165-232X(84)90006-5. [14] Qi J., Wang S., Yu F. (2013). A Review on Creep of Frozen Soils, In Constitutive Modeling of Geomaterials. Springer, Berlin, Heidelberg, pp. 129–133. DOI: 10.1007/978-3-642-32814-5_6. [15] Razbegin V.N., Vyalov S.S., Maksimyak R.V., Sadovskii A.V., (1996). Mechanical properties of frozen soils. Soil Mechanics and Foundation Engineering 33(2), pp. 35–45., DOI: 10.1007/BF02354292. [16] Li, D. W., Fan, J. H., Wang, R. H. (2011). Research on visco-elastic-plastic creep model of artificially frozen soil under high confining pressures. Cold Regions Science and Technology, 65(2), 219-225. DOI: 10.1016/j.coldregions.2010.08.006 [17] Wang, S., Qi, J., Yin, Z., Zhang, J., & Ma, W. (2014). A simple rheological element based creep model for frozen soils. Cold Regions Science and Technology, 106, 47-54. DOI: 10.1016/j.coldregions.2014.06.007. [18] Klein J. (1981). Finite element method for time-dependent problems of frozen soils, International Journal for Numerical and Analytical Methods in Geomechanics, 5(3), pp. 263-283. DOI: 10.1002/nag.1610050304. [19] Klein J., Jessberger H.L. (1979). Creep stress analysis of frozen soils under multiaxial states of stress. Engineering Geology, 13(1-4), pp. 353-365. DOI: 10.1016/0013-7952(79)90042-5. [20] Puswewala U.G.A., Rajapakse, R.K.N.D., Domaschuk L., Lach R. P. (1992). Finite element modelling of pressuremeter tests and footings on frozen soils. International journal for numerical and analytical methods in geomechanics, 16(5), pp. 351-375. DOI: 10.1002/nag.1610160505. [21] Arenson, L. U., Springman, S. M. (2005). Mathematical descriptions for the behaviour of ice-rich frozen soils at temperatures close to 0 C. Canadian Geotechnical Journal, 42(2), 431-442. DOI: 10.1139/t04-109. [22] Bray, M. T. (2013). Secondary creep approximations of ice-rich soils and ice using transient relaxation tests. Cold Regions Science and Technology, 88, 17-36. [23] Vyalov S.S., Gorodetsky, S.E., Zaretsky Yu.K., Gmoshinsky V.G, Grigoreva V.G., Pekarskaya N.K., Shusherina E.P., (1962). Strength and Creep of Frozen Soils and Design of Ice—Soil Retaining Structures. U.S.S.R. Acad. Press, Moscow (in Russian). [24] Andersland O.B., Ladanyi B., 2013. An introduction to frozen ground engineering. Chapman and Hall, New York. [25] Kostina A., Zhelnin M., Plekhov O., Panteleev I., Levin L. (2018). Creep behavior of ice-soil retaining structure during shaft sinking, Procedia Structural Integrity, 13, pp. 1273-1278. DOI: 10.1016/j.prostr.2018.12.260. [26] Terzaghi K., Peck R.B., Mesri G. (1996). Soil mechanics in engineering practice. John Wiley & Sons, New York. << /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