Microsoft Word - 205 completed IIUM Engineering Journal, Vol. x, No. x, 200x Ahmadi 17 REMARKS ON DEFORMATION IN SNOW WITH RESIDUAL STRESS AND GRAVITATIONAL FORCE BY DETONATION IMPACT SIAVOSH AHMADI Department of Mathematics, Tehran north branch, Islamic Azad University Tehran, Iran. s_ahmadi@iau-tnb.ac.ir ABSTRACT: Estimation of strain (deformation) and strain energy at the time of an explosion has been undertaken in the present investigation. A basic concept of explosive pressure which is required for the estimation of strain and strain energy has been drawn from the underwater explosion with suitable modification to the present studies. One of the most significant aspects of the present investigation is to augment the model by accounting it for gravitational force and residual stresses which are of realistic in nature. Studies consist of computational aspect by observing the variation of strain and strain energy with respect to material parameters such as Elastic modulus (E), poisons ratio (ν), density (ρ), and flow parameters such as peak pressure (Pm), sound velocity (Vs), and length of snow slab (L). One of the interesting observation found in the computational aspect is, the deformation is found to decrease as the velocity of sound (Vs) in snow increases which appears to be strange but found factual. The model has been compared with the other existing papers in the literature and found that, the present model yields to better results. ABSTRAK: Anggaran regangan (canggaan) dan tenaga regangan pada masa letupan diambil kira dalam kajian ini. Konsep asas tekanan letupan yang diperlukan untuk anggaran regangan dan tenaga regangan diperolehi daripada letupan bawah air dengan pengubahsuaian yang sesuai dalam kajian terkini. Antara aspek penting dalam kajian ini adalah penambahan model dengan mengambil kira daya graviti dan tegasan sisa yang sebenar. Kajian terdiri daripada aspek pengiraan, dengan memerhatikan variasi regangan dan tenaga regangan terhadap parameter bahan seperti modulus anjal (E), nisbah Poisson’s (ν), ketumpatan (ρ), dan parameter aliran seperti tekanan puncak (Pm), halaju bunyi (Vs), dan panjang papak salji (L). Antara pemerhatian menarik yang ditemui dalam aspek pengiraan ialah, canggaan didapati berkurangan ketika halaju bunyi (Vs) salji bertambah. Walaupun ini nampak asing, tetapi ianya terbukti benar. Model dibandingkan dengan kertas kerja lain dan didapati bahawa, model terkini memberikan keputusan yang lebih tepat. KEYWORDS: avalanche dynamics; fracture mechanics; detonation pressure 1. INTRODUCTION Avalanches in snow bound hilly regions are basically due to instability caused in the snow pack lying on the slopes. The causes for such instability are many. It may be due to pressure metamorphism [1] or due to temperature metamorphism [2] or it may be due to wind induced drift phenomenon [3-4]. The studies on avalanches will not only help in forecasting of it (much before the eventuality to occur) but also in taking precaution for the mitigation aspect (preserving either the structure or the lives once the forces IIUM Engineering Journal, Vol. x, No. x, 200x Ahmadi 18 (avalanche) can be estimated beforehand). Also, studies help in minimizing the property damage as result of an avalanche. In view of its importance, studies pertaining to avalanches are undertaken in the present investigations. Fracture mechanics deals with the propagation of crack as a result of an impact on the material and the relations there on stress and strain with material properties [Elastic modulus (E), Poison’s ratio (ν), density (ρ), cross sectional area (A), crack length (a), and so forth] [5-6]. Since instability in snowpack ultimately results into an avalanche, it is understood that, instability may induce cracks in the snowpack before sliding from the top hill. This crack, it is believed that, has links with fracture mechanics for studying stress and strain aspect of snow as a material. Basic concepts of fracture mechanics have been used in analyzing the avalanche related aspects. Avalanches are both natural and artificial one. In a natural avalanche, the release of snow mass from the top of the hill is by the natural phenomenon which is described in the earlier paragraphs. Whereas in artificial ones, the release of an avalanche, is basically due to an explosion in the snowpack from the remote location. This method of artificial triggering is by and large gaining popularity in snow bound hilly regions in the Western Countries since much of the habitation is located in the western part of the globe. Moreover this method is safe and economical one when compared to the destruction (avalanche) it may cause when released naturally. The studies in the present investigations have been basically focused onto the estimation of important flow parameters [velocity and acceleration] during the explosion in the snowpack. In addition to the focused parts, the studies have been also concentrated onto the estimation of other parameters such as strain and strain energy parameters. As is well-known that strain and strain energy is one of the important parameters in the estimation of avalanche formations as a result of an explosion. It is generally noticed during the artificial triggering that, not all explosive detonation results into an avalanche. Though there are various other parameters also required to be estimated ( metamorphic state of snow pack, aging of snow pack, slope aspect and so forth) before going for an avalanche forecasting however, it is generally believed that strain (εy) and strain energy (Es) are the most effective parameters to estimate and quantify for it during the explosion aspect. In view of its importance of these parameters, the present paper has been stressed on strain (εy) and strain energy (Es) parameters estimations in addition to velocity and acceleration aspects. 2. ANALYSIS The model has been aimed at first obtaining important flow parameters [velocity, maximum velocity, and acceleration] before going for the determination of strain (εy) and strain energy (Es). 3. ESTIMATION OF FLOW PARAMETERS The relation between the pressure and velocity at the time of an explosion on the snowpack (lying on slope of hill) is given by [7] 2.P(t) – ρ.Vs.V = m dV/dt (1) where P(t) is the pressure, Vs is velocity of sound in medium, V is the velocity of projectile, ρ is the density of the water, m is the mass per unit area (water).This relation has been modified to our present requirement of snow as a material. The proposed formula is written as : 2.P(t) – ρ.Vs.V + σg – σr = m dV/dt (2) IIUM Engineering Journal, Vol. x, No. x, 200x Ahmadi 19 where σg , σr ρ, and V are gravitational and residual stress, density of snow, and velocity of the snow pack respectively. P(t) is taken to be of the form from Ref. [8] as P(t) = Pm .exp(-t/T) (3) where Pm is the peak pressure, t is the time, and T is the time decay factor. This peak pressure is due to an explosion. Schewizer [9] has correlated the term (σg – σr) to snow density (ρ), gravity (g), and to slab depth by (Fig. 1) σg / σr =1.5, σ = σg – σr =(1/3) ρ.g.H.sinθ (4) where H is snow slab depth and θ is the slope angle. This relation [eqn.(4)] has been successfully employed in our earlier studies [10]. Fig. 1: Flow geometry. In order to solve eqn.(2), we need to have one condition (initial or boundary) and this has been taken to be of the from V(0) = 0. Using this condition eqn.(2) simplifies to V(t) = [2.Pm.T/( ρ.(L - Vs.T))].{exp(-Vs/L) – exp(-t/T)} + [(g.H.sinθ)/3].{1- exp(-Vs.t/L)} (5) having known velocity, maximum velocity of the snowpack (Vmax) soon after the explosion can be estimated by equating its gradient to zero. That is dV/dt = 0. Taking differential of (5) with respect to ‘t’ and equating it to zero we get the time to achieve this maximum velocity. tVmax = L.T.ln{[σ.( T.Vs - L) + 2.Pm.Vs.T ] / (2.Pm.L)}/(Vs.T - L) (6) Substituting this ‘tVmax’ for ‘t’ in equation (5), we get maximum velocity (Vmax) as: Vmax = [2.Pm.T/(ρ.[L-Vs.T])].{{[σ.(Vs.T-L) + 2Pm.Vs.T]/(2.Pm.L)} V s .T/(L-V s .T) – { [σ.(Vs.T-L) + 2Pm.Vs.T]/(2.Pm.L)} L/(L-V s .T) }} + σ/(ρ.Vs) [1-{[σ.(Vs.T-L) + 2Pm.Vs.T]/(2.Pm.L)} V s .T/(L-V s .T) }] (7) Knowing velocity from equation (5), acceleration can be obtained by taking differentiation with respect to “t”. That is a(t) = dV(t)/dt = [2.Pm.T/( ρ.(L - Vs.T))].{(-Vs/L)exp(-Vst/L)+(1/T) exp(-t/T)} + [(g.H.sinθ)/3].{(Vs/L exp(-Vs.t/L)} (8) Estimation of Strain and Strain Energy From the classical theory of elasticity [5,6], deformation energy (Wdef), stress (σyy) in y direction, and strain (εy) are related by IIUM Engineering Journal, Vol. x, No. x, 200x Ahmadi 20 Wdef = ∫ σyy dεy (9) where σyy and εy are normal stress and normal strain respectively. σyy and εy are related by [5-6] εx = [σxx - ν (σyy + σzz)] / E (10) εy = [σyy - ν(σxx + σzz)] / E (11) where E, ν, σxx , σyy, σzz , εx are Elastic modulus, Poison’s ratio, normal stress in x, y and z directions, and normal strain in the x-direction respectively. Assuming σzz and εx both equal to zero, eqn.(10) and (11) simplifies to εy =(1/E) σyy (1-ν 2 ) (12) σyy = E.εy/(1-ν 2 ) (13) we also know that Kinetic energy (KE) [Wk] is related by KE = Wk = ½ m.V 2 max (14) It is assumed that, explosion energy is same as the deformation (strain) energy and is equated to kinetic energy. we get ½ ρ.V2max = {E/[2.(1-ν 2 )].} ε2y (15) substituting for Vmax from equation (7) in eqn.(15), we get deformation (normal strain(εy)) as: εy = [ρ.(1-ν 2 )/E] 0.5 . [2.Pm.T/(ρ.[L-Vs.T])].{{[σ.(Vs.T-L) + 2Pm.Vs.T]/(2.Pm.L)} V s .T/(L-V s .T) – {[σ.(c.T-L) + 2.Pm.Vs.T]/(2.Pm.L)} L/(L-V s .T) }} + σ/(ρ.Vs) [1-{[σ.(Vs.T-L) + Pm.Vs.T]/(2.Pm.L)} V s .T/(L-V s .T) }] (16) substituting εy in equation (13) and using eqn.(9) we get strain energy (Es) which is assumed to be same as deformation energy (Wdef) as: Es = {E/[2.(1-ν 2 )].} ε2y (17) 4. RESULTS AND DISCUSSION Flow parameters velocity (V), acceleration (a), maximum velocity (Vmax), strain (εy), and strain energy (Es) have been computed and their variation have been observed with material parameters such as density, velocity of sound in the medium, elastic modulus (E), and Poison’s ratio. The computed values have been shown in Fig. 2-6. The present computed values have been compared with the model having no gravitational and residual stress [10]. In order to compute flow parameters we need data on peak pressure Pm, ρ, E, T, t, ν, εy, and the same has been taken from Ref [11-12]. These data have been shown in Table 1. The results indicate, the present computed values are lower (except for density, and velocity of sound Vs) in comparison to without gravitation and residual stresses [the solid line in Fig. 2-6 indicate for without gravitation and residual stresses]. The reason for such lower values could be attributed due to inclusions of the more realistic parameters in the present investigations which resulted into more accurate values. One of the IIUM Engineering Journal, Vol. x, No. x, 200x Ahmadi 21 interesting observations of velocity computation is, the present variations of V with Vs is found to be increasing and realistic whereas without gravitational and residual terms is found decreasing. The reason for this discrepancy could be basically due to subtraction of velocity of sound from peak pressure terms [eqn.(2)]. The time delay factor T did not find much effect on the maximum velocity computations [Fig. 3].The results on acceleration [Fig. 4] is found to be very distinct when compared with and without gravitational and residual stresses. The justification for this discrepancy could be given in similar conjunction with that of velocity computations. The variation of acceleration with Pm [Fig. 4] also found not to yield significant variations in comparison to earlier model [10]. The results on strain (εy) with Vs, L (slab length) and T indicates that, the addition of gravitation and residual increases the strain component whereas with Pm and ρ it decreases [Fig. 5]. The variation of strain energy (Es) with E, ν, and with εy shows that, strain energy increases with increase in E, ν and εy [Fig. 6]. Also it is observed that, the present computed values of Es with εy are found to be greater than that of without gravitational and residual stress. All the computed results [Fig. 2-6] have found to agree with the physics of flows and deformation property. For instance, increase in strain with peak pressure and decrease with elastic modulus adheres to results from classical physics. Table 1: Data of flow parameters. Parameters Typical values Range Peak pressure : Pm 300 pa 200-300 Distance of blast center: r 90m 90-9000 Time decay : T 3 s 2-15 Arrival time : t 0.1 s 0-2 Density : ρ 200 kg/m3 100-400 Slab length : L 4 m 2-8 Sound speed in snow : Vs 300 m/s 100-300 Elastic modulus : E 10 Mpa 0.5-10 Poison ratio : ν 0.2 0.1-0.3 Strain : εy 5.7×10 -5 2×10 -5 – 10 - 4 IIUM Engineering Journal, Vol. x, No. x, 200x Ahmadi 22 Fig. 2: Variation of velocity (V) with flow parameters (Vs, L). Fig. 3: Variation of maximum velocity (Vmax) with flow parameters (Vs,L). Fig. 4: Variation of acceleration (a) with flow parameters (Vs, L). Fig. 5: variation of strain (εy) with flow parameters (Vs, L). 100 120 140 160 180 200 220 240 260 280 300 0.005 0.01 0.015 0.02 0.025 0.03 100 120 140 160 180 200 220 240 260 280 300 9.2 9.4 9.6 9.8 10 10.2 V Vs with gravitational and residual forces(rigth column) without gravitational and residual forces(left column) 2 3 4 5 6 7 8 9.4 9.6 9.8 x 10 -3 2 3 4 5 6 7 8 9.8 10 10.2 V L 100 120 140 160 180 200 220 240 260 280 300 0 0.01 0.02 0.03 100 120 140 160 180 200 220 240 260 280 300 0.01 0.02 0.03 0.04 Vmax Vs 100 120 140 160 180 200 220 240 260 280 300 0 0.05 0.1 a V s 0 0.005 2 3 4 5 6 7 8 0 0.01 0.02 a L 0.005 0.01 100 120 140 160 180 200 220 240 260 280 300 0 1 2 x 10 -4 100 120 140 160 180 200 220 240 260 280 300 0 1 2 x 10 -4 εγ Vs 2 3 4 5 6 7 8 4.2 4.22 4.24 4.26 4.28 4.3 4.32 4.34 4.36 x 10 -5 2 3 4 5 6 7 8 5.64 5.66 5.68 5.7 5.72 5.74 5.76 5.78 5.8 x 10 -5 εγ L IIUM Engineering Journal, Vol. x, No. x, 200x Ahmadi 23 Fig. 6: Variation of strain energy with strain, elastic modulus. 5. CONCLUSION The present studies deal with the estimation of strain and strain energy at the time of an explosion as these parameters plays significant role in forecasting of an avalanche. The model has been developed with the basic concepts from underwater explosions with suitable modification to the present model. The flow variables (velocity, maximum velocity, and acceleration, strain and strain energy) have been computed for various material parametric effects (density, peak pressure, slab length, elastic modulus, poisons ratio, velocity of sound in medium) and found that the present computed results yield better results when compared with the earlier models. Also, the model observed to yield the results which could be quantifiable with the physics of flow and deformation theory from classical theories. The results need to be quantifiable with the avalanche formation which is our continued effort. REFERENCES [1] Palmer,A.G and J.R.Rice.The growth of slip surfaces in the progressive failure of over- consolidated clay.Proc.R.Sac.London.Scr.A, 332 (1591),518-527,1973. [2] P.S.Negi, R.K Aggrawal, D N Sethi and Uttam Chand, “Measurment of thermal conductivity of seasonal snow by a pc based transient thermal probe method”, Proceedings International Symposium on Snow Monitoring and Avalanches(ISSMA-2004), Snow and Avalanche Study Establishment(SASE), Manali, India, 67-70(2004). [3] RONALD I, PERLA M, AND MARTINELLI.Jr,”AVALANCHE HAND BOOK” U.S.Dep.Agric.,Agric.Handb.489,254p,1978. [4] PAUL M.B.FOHN,”SNOW TRANSPORT OVER MOUNTAIN CRESTS”, journal of Glaciology,Vol 26,469-492,(1980) [5] Dowling. N.E. : Mechanical behavior of materials (1993).M. Young, The Technical Writer's Handbook. Mill Valley, CA: University Science, 1989. [6] Kannien, M.F. and Popelar, C.H.: Advanced Fracture Mechanics Oxford University Press (1985). [7] Hollyer RS. Direct shock-wave damage to merchant ships from non-contact derwater explosion. Trans SNAME 1959;67:773–84. [8] C.F. Hunga,*, P.Y. Hsua, J.J. Hwang-Fuub: Elastic shock response of an air-backed plate to underwater explosion(in detonation file). Trans SNAME 2005:151-158. [9] Schweizer J. Review of Dry Snow Slab Avalanche Release. J. Cold Region Science and Technology, 30, 43-57 (1999). 3 3.5 4 4.5 5 5.5 6 x 10 -5 0.008 0.01 0.012 0.014 0.016 0.018 3 3.5 4 4.5 5 5.5 6 x 10 -5 2 4 6 8 10 12 x 10 -3 left column rigth column Es ε γ 0.1 0.15 0.2 0.25 0.3 0.35 0.0166 0.0168 0.017 0.0172 0.0174 0.0176 0.0178 0.018 0.0182 0.0184 0.1 0.15 0.2 0.25 0.3 0.35 3.75 3.8 3.85 3.9 3.95 4 4.05 4.1 4.15 4.2 x 10 -3 Es ν strain = 0.00003 strain = 0.00005 IIUM Engineering Journal, Vol. x, No. x, 200x Ahmadi 24 [10] Siavosh Ahmadi and R.N. Pralhad : Modeling propagation of explosive shock waves in snow on hilly top and prediction towards an avalanche release, Proceedings34 National Society For Fluid Mechanics and Fluid Power, BITS Mesra Ranchi, Dec. 10-12, 2007. [11] Donald Albert: Reduction of blast noise by a snow cover, Primary subject classification: 21.3.4; Secondary subject classification: 72.8, 2002. [12] Virendra Kumar and H. Muthurajan : “Artificial triggering of avalanches using explosives”, Proceedings International Symposium on Snow Monitoring and Avalanches(ISSMA-2004): Snow and Avalanche Study Establishment(SASE), Manali, India., 137-144, 2004.