<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJG</journal-id><journal-title-group><journal-title>Open Journal of Geology</journal-title></journal-title-group><issn pub-type="epub">2161-7570</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojg.2016.63012</article-id><article-id pub-id-type="publisher-id">OJG-64251</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Numerical Modeling of Explosive Loading on Strong Rock Mass with Discontinuity System in Depth 25 m
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>abi</surname><given-names>Nakhaei</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mohsen</surname><given-names>Jami</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Homan</surname><given-names>Najafizade</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Farhad</surname><given-names>Jahani Adimi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Seyed</surname><given-names>Hossein Mirzeinaly Yazdi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mohsen</surname><given-names>Pourkermani</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Nazanin</surname><given-names>Sarhadi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Geology, Zahedan Branch, Islamic Azad University, Zahedan, Iran</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>nabi.nakhaei@gmail.com(AN)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>07</day><month>03</month><year>2016</year></pub-date><volume>06</volume><issue>03</issue><fpage>117</fpage><lpage>130</lpage><history><date date-type="received"><day>11</day>	<month>December</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>4</month>	<year>March</year>	</date><date date-type="accepted"><day>7</day>	<month>March</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Effect of Explosion on underground structures where drilled on the stony materials can be affected by explosive strength and rock mass properties of tunnel. In this paper, it has been studied on effect of explosion of 10 tons TNT on the round tunnel (diameter: 10 m) in 25 m depth in the mass rocks (RMR &gt; 90) to estimate stability of tunnel to construct underground structure. Regarding to the studied stones, geo-mechanical indices of rock mass have been found to use as strength criteria in UDEC. In the next step, tunnel in real dimensions and founded indices is modeled with UDEC and analyzed stably. Explosion wave indices are found by UFC Instruction and other experimental relations to administer force from explosion wave on the model. Finally, on the base of FISH, indices of explosion wave are administrated in UDEC and tunnel is analyzed dynamically. According to modeling in depth of 25 m, amount of displacement of drilled tunnel in mass rocks (RMR &lt; 90) is very lower than standard criterion and the tunnel is very stable. Development
   
  of plastic zone, vertical and horizontal displacement of rock mass around the tunnel, speed variations and stress from explosion wave are modeled and recorded as the graphs and figures.
 
</p></abstract><kwd-group><kwd>Underground Structure</kwd><kwd> Plastic Zone</kwd><kwd> Round Tunnel</kwd><kwd> UDEC Software</kwd><kwd> UFC Instruction</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Today, Passive Defense is an important case in army. Passive defense is a civil defense and refers to a series of steps to use weapons do not need to run it possible financial damage to equipment and facilities critical military and civilian casualties or the extent of the damage and prevent losses may be reduced to a minimum. One of the most important principles of Passive Defense is to construct safe structures. Some places such as unity tunnels, underground hospitals and mines can be used as the safe structures for Passive Defense. It is necessary to study on these places and estimate their stability. On the one hand, explosion waves have caused dynamical stress on rock masses around the underground places and led to failure and convergence and new forms around the rock mass that it can be covered stability of underground spaces. Then, it is necessary to study on explosion and its effects on underground structures and stability of these places from explosive loads.</p><p>The most important aims of this paper include: study on the important factors in seismic loading effect on stability of underground structures, effect of depth on the rocks with discontinuity system and analyze tunnel stability and amount of wave damping. Generally, stability of these structures against explosive loads from bombs and rockets should be studied before using them as the passive defense structures. On the one side, destroying from rocks and bomb is one of the most complex mechanical processes that it has been studied carefully during 2 centuries ago. At present, three methods can be used to evaluate these cases: finding experimental relations on the basis of repeated tests, developing simple dynamic splash methods, penetration and explosion by using laws of motion and materials and advanced numerical methods to model explosion load effects on underground structures. Methods based on experimental relations are very expensive and especial for each property. Then, they cannot be used for all conditions. Methods based on engineering models are simple with some information. They are ideal to evaluate and study on results from explosion loads. But they are not precise to analyze because of using making simple hypothesizes.</p><p>On the other side, advanced numerical methods have considered all cases and provided better results. There are some problems to use numerical methods such as need to precise input information, time to compute, high computational cost, precisely understanding solutions of problem. Because of increasing computational power of computers and saving basic information, researchers and scientists would like to use numerical methods recently.</p><p>An explosion is a rapid increase in volume and release of energy in an extreme manner, usually with the generation of high temperatures. TNT is an explosive material and it is stable and pure chemically and it can be compared with various explosive materials [<xref ref-type="bibr" rid="scirp.64251-ref1">1</xref>] . As mentioned above, an explosion is a rapid increase in volume and release of energy in an extreme manner. Based upon it, explosive materials have various explosive specifications regarding type and methods of combination. It consists of:</p><p>Released energy, fast explosion that fast explosion of TNT is about 2000 f/s.</p><p>There are many important variables that effect on exploring rocks such as properties of stone and explosive materials and explosive patterns. Each explosion has two parameters: Pressure and Volume. Two parameters of distributing different components of energy in explosion process are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>On the basis place of exploring to the earth, Explosions can be categorized in three groups [<xref ref-type="bibr" rid="scirp.64251-ref2">2</xref>] :</p><p>Explosions that placed up to the earth.</p><p>Explosions that placed on the earth (model of this paper).</p><p>Explosions that placed under the earth.</p></sec><sec id="s2"><title>2. Background</title><p>Sheng-tian and Zhi-yuan (1993) have described dynamic interaction of underground structures and local rocks by using Bending Equations of thick and thin plates [<xref ref-type="bibr" rid="scirp.64251-ref3">3</xref>] . Also, Chen and Zhao (1998) have studied on wave propagation process from explosion in discontinuous and jointed rock masses by using UDEC and AUTODYN [<xref ref-type="bibr" rid="scirp.64251-ref4">4</xref>] . Morris et al. 2004 have presented results from loading of tunnels placed on Nevada Site (they are placed in granite and tuff areas with maximum stresses) by using method of three dimensional separated elements [<xref ref-type="bibr" rid="scirp.64251-ref5">5</xref>] . Lu (2005) has studied on effect of different parameters including geological conditions and density of explosive material on explosion in Swedish tunnels and modeled these results by using Artificial Neural Network [<xref ref-type="bibr" rid="scirp.64251-ref6">6</xref>] . Heuze and Morrise (2006) have performed laboratory and desert tests to model effect of explosion in jointed rock mass [<xref ref-type="bibr" rid="scirp.64251-ref7">7</xref>] . Jiao et al. (2007) has studied on effect of wave propagation from explosion by using Discontinuous Deformation Analysis (DDA) in jointed rock mass [<xref ref-type="bibr" rid="scirp.64251-ref8">8</xref>] . Mortazavi (2003) has evaluated role of strength properties of discontinuities on stability of underground spaces under dynamic load by using Discontinuous Deformation Analysis [<xref ref-type="bibr" rid="scirp.64251-ref9">9</xref>] . Also, Mirzeynali and Bejarivalam (2011) have studied on numerical modeling of explosive loading on weak continuous rock mass by using Finite Difference Method with FLAC (A 2 Dimentional</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Computing explosion wave parameters [<xref ref-type="bibr" rid="scirp.64251-ref11">11</xref>] </title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1210380x7.png"/></fig><p>Software) [<xref ref-type="bibr" rid="scirp.64251-ref10">10</xref>] .</p></sec><sec id="s3"><title>3. Method of Computing Explosion Wave Parameters</title><p>Today, both experimental and numerical methods can be used to compute wave parameters. In this paper, we use experimental methods and UFC Instruction [<xref ref-type="bibr" rid="scirp.64251-ref11">11</xref>] .</p><p>According to UFC Instruction, <xref ref-type="fig" rid="fig1">Figure 1</xref> should be used to compute wave parameters.</p><p>where</p><p>W: weight of explosive material (on the basis of pond (Ib)).</p><p>P: Pressure from explosion (on the basis of PSI).</p><p>R: Distance from explosion place (on the basis of ft.).</p><p>T: Time (on the basis of ms).</p></sec><sec id="s4"><title>4. Modeling</title><p>This rock mass is very strong (RMR &gt; 90). Complete Mechanical Specifications of joints and the rocks have been studied in <xref ref-type="table" rid="table1">Table 1</xref> [<xref ref-type="bibr" rid="scirp.64251-ref2">2</xref>] .</p><sec id="s4_1"><title>4.1. Computing Wave Parameters from Explosion</title><p>Both experimental and numerical methods have been used to compute wave specifications from exploring 10 tons TNT.</p><p>Based upon it, necessary specifications to find wave coefficients are as follow as:</p><p>Distance from Explosion Place (R): place of measuring wave parameters has been placed at a distance of 1.5 m from explosive material, then</p><disp-formula id="scirp.64251-formula1"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1210380x8.png"  xlink:type="simple"/></disp-formula><p>Weight of Explosive Material (W): 10 tons TNT.</p><disp-formula id="scirp.64251-formula2"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1210380x9.png"  xlink:type="simple"/></disp-formula><p>Then, Parameter Z will be computed to be used UFC Instruction that is equal to:</p><disp-formula id="scirp.64251-formula3"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1210380x10.png"  xlink:type="simple"/></disp-formula><p>At the next stage, pressure from explosion is about psi13250 that is equal to 90.1 as following:</p><disp-formula id="scirp.64251-formula4"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1210380x11.png"  xlink:type="simple"/></disp-formula><p>Finally, refer to UFC Instruction (<xref ref-type="fig" rid="fig1">Figure 1</xref>) to compute Parameters of Explosive Wave.</p></sec><sec id="s4_2"><title>4.2. Tunnel Modeling by UDEC</title><p>Parameters from explosion wave have been modeled in a tunnel, 25 m in depth. Then, Tunnel Stability will be studied in this depth. After primary modeling in UDEC, the drawn model will be solved because unbalancing powers should be zero (<xref ref-type="fig" rid="fig2">Figure 2</xref>).</p><p>After wave damping from explosion, developing plastic area, displacement cantores and speed variations are according to Figures 3-5.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Mechanical specifications of rock mass and joints to model [<xref ref-type="bibr" rid="scirp.64251-ref2">2</xref>] </title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Mechanical Properties of Rocks</th><th align="center" valign="middle" >Mechanical Properties of Joints</th></tr></thead><tr><td align="center" valign="middle" >Dens: 0.002700</td><td align="center" valign="middle" >Jkn: 10e6</td></tr><tr><td align="center" valign="middle" >G: 3e3</td><td align="center" valign="middle" >Jks: 10e6</td></tr><tr><td align="center" valign="middle" >K: 5e3</td><td align="center" valign="middle" >Jten: 100</td></tr><tr><td align="center" valign="middle" >Ucs: 100</td><td align="center" valign="middle" >Coh: 100</td></tr><tr><td align="center" valign="middle" >Prat: 0.25</td><td align="center" valign="middle" >Fric: 40</td></tr><tr><td align="center" valign="middle" >Coh: 20</td><td align="center" valign="middle" >Spacing: 30 cm = 0.3 m</td></tr><tr><td align="center" valign="middle" >Fric: 25</td><td align="center" valign="middle" >-</td></tr><tr><td align="center" valign="middle" >RMR &gt; 90</td><td align="center" valign="middle" >-</td></tr></tbody></table></table-wrap><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Variations of unbalancing powers in the model</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1210380x12.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Developing a plastic area in a tunnel (depth: 25)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1210380x13.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Displacement cantors variations in a tunnel (depth: 25)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1210380x14.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Cantors variations in a tunnel (depth: 25)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1210380x15.png"/></fig><p>Stress changes of the explosion in the ceiling, walls and floor of the tunnel is in accordance with the images below.</p><disp-formula id="scirp.64251-formula5"><graphic  xlink:href="http://html.scirp.org/file/1-1210380x16.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.64251-formula6"><graphic  xlink:href="http://html.scirp.org/file/1-1210380x17.png"  xlink:type="simple"/></disp-formula><p>Following graph shows displacement variations (vertical and horizontal) in ceiling, walls and floor of tunnel, respectively.</p><disp-formula id="scirp.64251-formula7"><graphic  xlink:href="http://html.scirp.org/file/1-1210380x18.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.64251-formula8"><graphic  xlink:href="http://html.scirp.org/file/1-1210380x19.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.64251-formula9"><graphic  xlink:href="http://html.scirp.org/file/1-1210380x20.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.64251-formula10"><graphic  xlink:href="http://html.scirp.org/file/1-1210380x21.png"  xlink:type="simple"/></disp-formula><p>Displacement from an explosion in a tunnel (depth: 25)</p><disp-formula id="scirp.64251-formula11"><graphic  xlink:href="http://html.scirp.org/file/1-1210380x22.png"  xlink:type="simple"/></disp-formula><p>Horizontal Displacement</p><disp-formula id="scirp.64251-formula12"><graphic  xlink:href="http://html.scirp.org/file/1-1210380x23.png"  xlink:type="simple"/></disp-formula><p>Vertical Displacement</p></sec></sec><sec id="s5"><title>5. Conclusions</title><p>Results from this modeling are as following:</p><p>1) Developing is limited in depth of 15 m.</p><p>2) Minimum overload altitude is 20 m for underground structure to resist against explosion.</p><p>3) Optimal depth and stability of a tunnel is more than 25 m regarding this rock mass.</p><p>4) Studying on Tunnel Behaviors such as displacement and strain on the ceiling indicate that the tunnel will be stable if it is placed in every depth and under dynamical load equal to 10 Tons TNT but optimal depth to drill a tunnel with minimum strain is 20 - 30 m.</p></sec><sec id="s6"><title>Cite this paper</title><p>NabiNakhaei,MohsenJami,HomanNajafizade,Farhad JahaniAdimi,Seyed Hossein MirzeinalyYazdi,MohsenPourkermani,NazaninSarhadi, (2016) Numerical Modeling of Explosive Loading on Strong Rock Mass with Discontinuity System in Depth 25 m. Open Journal of Geology,06,117-130. doi: 10.4236/ojg.2016.63012</p></sec><sec id="s7"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.64251-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Smith, P.D. (1994) Blasting and Ballistic Loading of Structures. Butterworth-Heineman Ltd.</mixed-citation></ref><ref id="scirp.64251-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Arab, M.R. and Mirzeynali, S.H. (2001) Numerical Modeling of Explosive Loading on Strong Rock Mass with Discontinuity System. Malek Ashtar Industry University.</mixed-citation></ref><ref id="scirp.64251-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Yang, S.T. and Cao, Z.Y. (1983) An Analytical Solution for Underground Structure Rock-Dynamic Interaction. Applied Mathematics &amp; Mechanics, 4.</mixed-citation></ref><ref id="scirp.64251-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Chen, S.G. and Zhao, J. (1998) A Study of UDEC Modeling for Blast Wave Propagation in Jointed Rock Masses. International Journal of Rock Mechanics and Mining Sciences, 35, 93-99.  
http://dx.doi.org/10.1016/S0148-9062(97)00322-7</mixed-citation></ref><ref id="scirp.64251-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Morris, J.P., Rubin, M.B., Blair, S.C., Glenn, L.A. and Heuze, F.E. (2004) Simulations of Underground Structures Subjected to Dynamic Loading Using the Distinct Element Method. Engineering Computations, 21, 384-408.  
http://dx.doi.org/10.1108/02644400410519848</mixed-citation></ref><ref id="scirp.64251-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Lu, Y. (2005) Underground Blast Induced Ground Shock and Its Modeling Using Artificial Neural Network. Journal of Computers and Geotechnics, 32, 164-178. http://dx.doi.org/10.1016/j.compgeo.2005.01.007</mixed-citation></ref><ref id="scirp.64251-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Heuze, F.E. and Morris, J.P. (2006) Insights into Ground Shock in Jointed Rocks and the Response of Structures There-In. International Journal of Rock Mechanics and Mining Sciences, 44, 647-676.  
http://dx.doi.org/10.1016/j.ijrmms.2006.09.011</mixed-citation></ref><ref id="scirp.64251-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Jiao, Y.Y., Zhang, X.L., Zhao, J. and Liu, Q.S. (2007) Viscous Boundary of DDA for Modeling Stress Wave Propagation in Jointed Rock. International Journal of Rock Mechanics and Mining Sciences, 44, 1070-1076.  
http://dx.doi.org/10.1016/j.ijrmms.2007.03.001</mixed-citation></ref><ref id="scirp.64251-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Mortazavi, A. (2008) Numerical Analysis of Strength Properties of Discontinues on Stability of Underground Spaces under Dynamic Load. The Sixth Tunnel Conference, Civil Engineering College of Tehran Science and Indust University.</mixed-citation></ref><ref id="scirp.64251-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Bajarivalam, H. and Mirzeynali, S.H. (2010) Modeling of Explosive Loading on Weak Rock Mass with Continuity System. Malek Ashtar Industry University.</mixed-citation></ref><ref id="scirp.64251-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">TM 5-1300 (UFC 3-340-02) US (1990) Army Corps of Engineers, Structures to Resist the Effects of Accidental Explosions. US Army Corps of Engineers, Washington DC.</mixed-citation></ref></ref-list></back></article>