Dimensioning of Room and Pillars in Horizontal Layers

An optimal dimensioning of the room and pillars of the square or rectangular section is nding adequate transversal section support for the pillar; in that sense, all forces acting on the pillar must be evaluated. e pillar has its ends subjected to normal (σ) and shear stresses (τ). In addition, the pillar can also be subjected to dynamic stresses from the seismic waves generated by rock blasting. us, when calculating the dimensions of a mine pillar, the stresses imposed by seismic waves must be considered. e combination of lithostatic and dynamic stresses generated by blasting can aect the pillar, generating regions, or areas that will exhibit the formation of cracks, microcracks, chipping, and other detrimental processes for its stability. In that regard, it is appropriate to assess the dynamic eects of the energy delivered in each wave pulse so that these dynamic quantities can be included in the dimensioning of the pillars. is paper presents the physical and mathematical formulation describing the dimensioning issues of room and pillars and its solution using numerical modeling. e numerical modeling is performed using nonlinear mathematical programming and simulation with FLAC 3D (Fast Lagrangian Analysis of Continua in 3D). In this study, it was proven that the stress generated by the seismic wave represents a dynamic stress of nearly 5.5MPa. 'Dimensioned pillars considering wave stresses allows it to recover an average 1.6% less when it is compared with pillars only under overlying rock mass'.


Introduction
Mining is an activity whose purpose is to extract minerals that society needs and uses economically. Mining requires that the nancial resources be e ectively allocated [1]. Mining operations include the activities of drilling, rock blasting, loading, and transportation [2][3][4]. Rock excavation is the rst activity in mineral production. It can be achieved using mechanical methods or with the use of explosives. Traditionally, at industrial-scale production, explosives are used for rock excavation. e explosives must be distributed in holes systematically in the rock mass to be blasted. is task is made to fragment the rock by releasing the explosive's energy. It follows the principles of physic dynamics, cinematic, and the expansion of hypercompressed gases. Not all the energy released by the explosives is used in the fragmentation process, resulting in collateral e ects such as y rock, airblast, and seismic waves.
As a product of the detonations, the seismic waves will result in the loading and unloading of the rock mass in a repetitive manner, similar to the conditions of a body subjected to laboratory tests of simple cyclical compression [5,6]. e repetitive application of stresses can be detrimental for the geotechnical parameters of the materials (Young's modulus (E), shear modulus (G), constant K, and yield stress (σ e )). Additionally, the magnitude of the tensile strength (σ t ), compression (σ c ) of the rock mass, and the factor of safety (FS) of mine structures are constantly subjected to variations as a consequence of the abovementioned dynamic stresses [7].
Regarding a mine that uses the room and pillar mining method, the mine pillar is a structure subjected to lithostatic stresses (primary stresses), corresponding to the weight of a prismatic volume of an overlying rock that acts on its top [8].
e conventional sizing of a mine pillar does not include the effect of the wave stress produced by rock blasting on these structures. For this reason, the methodologies that allow integrating these dynamic loads in the sizing of mine pillars must be adopted. e purpose of finding the pillar's optimal dimensions is to obtain the safest geometry for the room and pillars system with the maximum recovery of the minerals in the deposit; hence, a detailed evaluation of the dynamic stresses acting on the pillar must be performed; accordingly, when dynamic stress is simplified under the premise "its value is negligible," such assumption characterizes an error in the dimensioning of a pillar, since pillar stays overloaded (underdimensioned because of simplified estimates to find its dimensions). ereby, there is a risk that pillar is working in transient regime close to its yield stregth limit, and 'any increasing' stress on the pillar would drive its strength parameters to break. erefore, it is necessary to carry out rock mass monitoring with geophones/seismograph for measuring dynamic problems whose aim is controlling the energy released by rock blasting.

Methodology
is paper uses a novel methodology to include the stresses generated by blasting around the pillars for its dimensioning. e methodology can be summarized as follows: (1 )Analysis and evaluation of recorded seismograph information, lithological information, and geotechnical data of the area under study are achieved (2 )Characterization of the dynamic properties of the rock mass using the Fourier Transform to understand the mechanism of energy dissipation in the rock mass domain includes the application of conservation theories of energy, mass, as well as determination of the particle velocity responsible for displacements and deformation, and contrasting this information via mathematical models and simulation with a commercial software (I-Blast V.8) (3) Formulation and development of equations integrate the sizing mine pillar model (4) Design one algorithm that provides the solution for dimensioning room and pillars through nonlinear mathematical programming and simulation with FLAC 3D (5) Finally, there is analysis of results

Seismic Waves Generated by Rock Blasting.
e study of oscillations and vibrations is a fundamental part of dynamic studies in physics [9]. In general, the vibrational wave is attenuated and can be characterized as a damped vibration where a portion of the energy delivered by the wave is 'absorbed' by the rock material [6]. e deformation produced in a material is proportional to the applied tension and can be represented by the postulates of Hooke's law (σ � Eε, with σ being the stress, and E and ε are Young's modulus and unit strain, respectively). e detonation of an explosive charge in the rock mass in a blast hole imposes a dynamic loading in the vicinity of the structures, and, depending on the magnitude of these dynamic loads, the rock mass will begin to fail [10].

Conservation of Movement, Conservation of Mass, and
Energy. e work developed by force 'F' is equivalent to the variation of energy (dW � − dU); such a statement is supported by the reduced equation of Bernoulli: where V, P, ρ, and K are the velocity, pressure, density, and constant, respectively. e gases responsible for the fracture are dissipated by exerting an internal pressure in the middle; from the Bernoulli equation, it is verified that the pressure at any point in the rock mass can be found once the pressure (inside the hole) and the velocity at an initial point are known.
Tesarik and Hustrulid [11], used hydrodynamics concepts to assess the peak particle velocity (theoretical velocity) at a distance from the origin due to an explosive charge. is methodology is based on the continuity equation (∇•(ρv) � − (dρ/dt), where ρ is the density, t is the time, and v is the velocity vector. e explosive energy fragments the rock in the process, resulting in a material with the behavior of a flow (by an area element [DϵR 3 ]). e flow of the material can be described as work with a domain (D), and the flow problem can be analyzed using the Stokes theorem (the first following equation) and Gren's Potential theorem (the second following equation). By this way, the equilibrium relationships for the system in motion are expressed by the movement conservation equations (the third following equation) [3,[12][13][14].
where the variable τ represents the stress acting on the elementary unit, ρ is the volumetric mass (density), u, v, an d w are the velocity vectors, g is the gravitational component, ∇ is the gradient operator, and x, y, and z are the planar coordinates. e conservation of mass is applied to fluids (flows) in motion, and its generalized equation is expressed through the integral surface of a vector function ( is an elementary area). is analysis will consider some basic criteria accepted by Tesarik and Hustrulid (2009): irrotational flow (∇ × V → � 0) and rock incompressibility (∇• V → � 0 or zρ/zt � 0) have a behavior similar to an inviscid liquid.
Particles that carry the seismic wave are subjected to stages of excitement, and the wave passage is responsible for the potential energy of deformation and the kinetic energy [9]. e quantity of energy released in each wave pulse is related proportionally to the amplitude. Its magnitude depends on the parameters selected in the drilling operation, rock mass properties, and blasting sequencing [15]. PAN-EIRO states that the liberated energy is suddenly in a short interval of time in the form of high pressures and temperatures, which are sufficiently big enough to produce the fragmentation of the material, and the residual stresses dissipated in the vicinity produce the vibrations in the ground [16]. e internal energy (U) is generated from the sudden disintegration of the explosive, characterized as an "adiabatic process," which is transformed into energy work and is manifested by rock fragmentation and subproducts. e blasted rocks under this kinetic process follow the principles of the Newtonian flow. Considering the criterion of continuity and conservation of mass and energy, the equations derived from Navier-Stokes allow finding an implicit way through which the energy dissipates by the wave in the rock volume. is law describes the movement of the flow in the middle, and, for that, the following premises for the flow in question must be assumed: density (ρ) and viscosity (μ) are constant.

Vibrations and Stress on the Ground.
After detonating the explosive material in the blast hole, a set of stress waves are generated, traveling through the medium as vibrations. Many researches proposed equations (Duvall and Fongelson, 1962, Ambraseys and Hendrom, 1968, and Indian Standard Predictor, 1973) to estimate those particle velocities. e velocity can be measured using seismographs [17], and the signal information obtained from the seismograph records provides relevant data about the behavior of the ground.
e 'maximum standby charge method or scaled distance' is amply accepted to predict a vibration level at a certain point. One of the methodologies to determine the PPV, as well as the characteristics of the ground, is expressed by Devine's equation [3].
With PPV being the peak particle velocity, D is the distance between the location of the explosive charge and the point of interest, Q is the amount of standby charge (explosive charge), and α and K are the attenuation and the ground coefficients, respectively. In addition, for determining the PPV in a short field (close to the explosive charge), the Holmberg and Persson equations can be used, and the PPV may also be determined by the hydrodynamic theory as mentioned in Section 2.2; this criterion has a physical-mathematical basis in fluid dynamics [18].
Zhang admits to an empirical relation to find the deformations made by a wave from a distance 'r' in the detonation point where K is the rock mass constant, B takes the value of 3 (however, this value has a mathematical basis and can be found from the rock mass characterization and the function decoded via FFT), the constant r h is the radius of the hole, and a constant β is the attenuation coefficient [19].

Explosive Energy.
e shock of explosive energy relates to the detonation velocity of the explosive through the following expression: With P d , being the detonation pressure (GPa), ρe is explosive density; VOD is explosive detonation velocity (m/ s); and c indicates the specific heat ratio of the gases produced by the detonation [20]. On the other hand, the energy delivered by the explosive to fragment the rock is a function of the pressure mentioned above; the explosive density is corrected in the Chapman-Jouguet plane (ρ CJ � (4/3)ρe) and Gurney's factor (X � 0.68 − 0.21ρe) for explosive mixtures with ANFO (Ammonium Nitrate-Fuel Oil) and X � 0.2 + 0.37ρe for mixtures with emulsion [19].

Frequency Analysis of Wave Signals.
e frequency analysis of the seismic waves can be done using the Fourier Transform.
e Fourier Transform allows converting the time-domain waveform in the frequency-domain [21]. is is done using where x(m) is the magnitude data (frequency information); x(n) is the waveform signal in a discrete form; w, m, n, t, and N are the angular frequency; and m is the component associated with each sample spectrum, n is the component associated with the parcels that characterize the series of coefficients, t is the time, and N is the amount of data.

Assessment of the Stresses Generated by Energy Dissipation through Waves on Support Structures: Mine Pillar.
A general equation for the propagation of a plane wave (on the x axis) in a rock medium, continuous and isotropic: is equation is a simplification of a 3D wave, whose expression is defined by Deformations (ε) are proportional to the wave velocity in the rock mass, which is ε � ε(ξ) with 'ξ' being the wave velocity (v) at point of the analysis; additionally, the seismic energy can be expressed as a function of the velocity. erefore, it is possible to write an equation that represents the behavior of the disturbance based on the following mathematical model: where U o is the fundamental amplitude, k is the constant attenuation, c is the peak particle velocity, and x is the distance [11].
In the same order of ideas, the paper in [9] establishes a relation of energy associated with the deformation of the material and the acting stress is differential expression is associated with variable amplitude; the author considers a solution based on a complex conjugate model to represent the decay of the wave with time and formulates the following function: z � Ae (pt+∝)j and its implicit solution: (− p 2 + jpc + ω o )Ae (pt+∝)j � 0; the constants A and ∝ allow adjusting the solution to the values of displacement and velocity, while p is a value associated with variable angular frequency, c is a constant related to the viscosity of the material, and ω o is the natural frequency.
According to the expressions proposed by the authors mentioned above, it is evident that the variable amplitude can be written as follows: ese different relations and the formulation of Devine can be equated to define the magnitude of the stress of the acting wave that produces the deformation.
On the other hand, Goodman sustains that knowing Young's modulus values and Poisson relation (μ), the deformations are defined for a state of three-dimensional stresses in an elastic solid: isotropic and linear [5]. From the concepts of stress-deformation and the theories by the above-mentioned researchers, it is possible to infer the following equation, which describes the state of stress produced by the passage of the wave in the rock mass. Although the design of the following equation is to consider Hooke's law, it has complex solution elements (attenuation coefficient); its solucion depends on the elastic parameters of the rock mass, expression representing the displacemente component (z/z(x ′ )), time, and the position of the pulse (v).
Note that, in a more restricted concept in (9), the variable x seems to be the position of the pulse responsible for the deformation (different from the general concept of 'wave path'), and, in this understanding, x is replaced by the affine x', not to create ambiguity in the understanding; thus, x' becomes an 'independent' and differentiable variable. At the same time, it is a function linked to the dimensions of the spans and the pillars. Analogously, the variable t can also be conveniently rearranged as an expression of velocity (speed of the wave in the process of attenuation and the dimensions of the pillar and spans).

Dynamic Stress on the Mine Pillar.
e function u(x, t � (x � /v)) defines the behavior of the 'displacements' of the particles when subjected to dynamic stresses, and the variables x, x´, t, and v express the pulse position, time, and speed.
ere is a function zu(x, t)/zx that expresses the deformations (ε(ξ) � (zs − zx)) associated with the variable tension (∆σ � zσ); such relations are linked implicitly to the velocity of the wave dissipated in the rock mass. e generalized expression of Hooke's law allows us to relate the 'stresses and strains' since the size must meet the stability criteria: On the other hand, (10) can be written as a function of the characteristics of the traveling wave propagating in the rock mass. e ∆σ component represents the dynamic stress, and its equation is a complex derivate expression, whose solution is composed of a real and an imaginary part In an attempt to facilitate the presentation, the math presentation of (11) is summarized throughout symbolic operators (a, A, B, C, D, E, and F). 4 Advances in Materials Science and Engineering e geometric variables (W p , L p , W o , and L o ) are conjugated of the tributary area dimensions; E m , is the rock mass elasticity modulus; ρ rm is rock mass density; ρ e is density of explosive; E U is the useful energy provided after the detonation of the explosive, f h is a factor associated with the amount of charge (ANFO) per hole, r is the radius of the hole, and μ is called Poisson's ratio. Figure 1 is a schematization of the process of mineral extraction by the room and pillar method, in which the origin of the coordinates (O) is located in the center of the room section (W o H p ), x, u(x � 0, t)). e dynamic tension (σ o ) active in the pillar ( Figure 1) is defined by the quantity of charge delivered in each pulse of energy when the wave runs through the distance 'x-vt' in Figure 2 on the horizontal axis OX.
erefore, the value of σ o, must be added in the lithostatic stress acting on the pillar (σ p ).

Case Study
e sizing problem of the pillars presented in this document refers to a real situation that happened in a mining company situated in Brazil. e company performs mechanized mining of mineral resources. In this study, an emphasis has been put on the room and pillar sizing, considering the stresses delivered by the dissipated waves of explosive detonations; after these produced energies are quantified, they will soon be incorporated in the formulation that integrates the algorithm to find optimal dimensions of the section of pillar support. It is opportune to quote the resistance criterion of Hoek and Brown's 2002 version, which considers the empirical factor 'D'; such a factor 'D' correlates the damage produced by the excavation of rocks with explosives and the resistance of the rock mass in an empirical way. However, with the support of physical tools-mathematics by what is cited in the methodology-the tools of the wave are studied and included in the sizing process of the room and mine pillars. us, the sizing problem is patterned and the results obtained represent information with bigger accuracy in contrast with simplified models or empirical foundations, in which the wave voltages are disregarded.
Regarding the characteristics of the deposit, the ore layer is at a depth between 300 and 350 meters and has a specific weight of 0.027 MN/m 3 to a barren rock (overburden), whose mineral occurrence is given in layers of 2 m of thickness in a horizontal position. e roof that comprises the spans shows a sedimentary structure (slab) of 1 m of thickness, and the mineralized rock is a quartzite; in Table 1, the lithological materials as well as the geotechnical data for these materials are broken down. e mining company in question uses the ANFO (total energy of the explosive: 912 kcal/kg) explosive, which has a specific mass of 900 kg/m 3 . e approximate loading ratio of 3.3 kg/m 3 (4.5 kg/hole) for a stope whose average section is 4.3 m × 2 m, with holes of 45 mm in diameter and 45 holes, is distributed in the blasting stope.

Analysis of Ground
Characteristics. For this case study, the seismograph data was recorded in a stope, with a standard section of 4 m × 2 m; the sensors were installed at distances between 30 and 75 m from the stope. Table 2 includes the information corresponding to (1) 'peak velocity (V)' corresponding to the channel 'L.' e L channel represents seismographic data with higher intensity information and, therefore, the wave energy is increased and (2) scaled distance is illustrated in the same table (SD (D/ �� Q √ ), with D being the distance between the stope and the sensor � 70 m, Q � 70 kg). e parameters for the rock mass are 'K � 153.04 and α � − 1.463'. e information for K and α Advances in Materials Science and Engineering is used to analyze the wave velocity in a short eld in the software I-Blast V. 8.

Analysis of Seismographic Data through Fourier
Transform.
e information in Table 3 includes the result of the analysis of the seismograph records after being analyzed using (7). e table includes (a) phase angle (φ) − (π/3), (b) amplitude (Uo) ≈ 0.45 m, and the fundamental frequency (f ) 48.7 Hz with the dominant angular frequency (ω) 306 rad/s.

Analysis of the Velocities Responsible for the Stresses in the Pillar.
e software I-Blast V.8 simulates the wave propagation, derived from information contained in Section 3.1 and Table 2. In this way, the stress found in the centroid of the pillar is an implicit function of liberated energies, which are expressed in scalar magnitudes of velocity. It must be highlighted that the simulation provides the peak velocity of a particle; therefore, these wave magnitudes represent the maximum value of the wave pulse velocity.

Lithostatic Stresses and E ective Stresses on the Pillar.
ese are stresses related to the lithological load exerted by the rock mass on the pillar that works as a support at the applied depth. Table 4 presents the average lithostatic stress, as well as the pillar e ective stress (transversal section of pillar support ≈ 10.3 m 2 and tributary area ≈ 72 m 2 ). Such information corresponds to the historical data provided by the geomechanical department.

Approach to the Problem.
e problem of pillar dimensioning, including the dynamic stresses, was solved using numerical modeling through nonlinear mathematical programming (PNL). It is possible to establish those recovery functions of deposits whose variables are the geometric unknowns, which must be respected in order to guarantee the mining arrangement. e combination of the recovery functions and the geometric restrictions must be written in a function of structural interdependence that characterizes the algorithm as a whole, which will allow obtaining the sought answers. In addition, the information found in the mathematical model, in which the dynamic forces are incorporated, will be evaluated by simulation.

Math Formulation.
Betts shows a solution set, in which the structure establishes the objective function and the restrictions to be respected [22]; it is as follows: evaluate the n-vector: x T [x 1 , x 2 , . . ., x n ], to maximize (or minimize) the function F F(x i ).
Subject to 'm' restrictions, C L ≤ C(x) ≤ C U . e limits of coherence: x L ≤ x ≤ x U ∀x ? R n /R ⟶ C r , where x T is the n-dimensional vector of the project variables; F(x) is the objective function to be optimized; C(x) corresponds to the restrictions that condition the solution; the values of x represent the permissible magnitudes associated with the coherence of the model and decide the space of solutions in the domain R n , whose values for the variables are of the real type within the closed-convergent bound that encompasses the solution, and the expressions are continuous and di erentiable functions (C r ).

Objective Function.
e math expression of the objective to be reached must involve the corresponding variables to the geometric data of the mining arrangement (see Figure 3); in this way, the recovery of the function (Rec) is expressed by Rec Rec (W p , W o , L p , and L o ), with W p , W o , L o , and L p being the pillar width, span width, pillar length, and span length, respectively.
Considering a regular arrangement and a uniform distribution of the pillars in the area corresponding to the mineral block, as well as assuming that the entire thickness (H p ) of the ore layer is mined, the recovery function  (objective function) is defined as follows:

Tributary Area eory (TAT).
is method is based on the proposition of Newton's third law [6,[24][25][26]. For this, the strength of pillar equals the weight of the volume of overlying prismatic-shaped rock, overlying the vertical reaction exerted over an influential area (At) or "tributary" (such dimensions include the pillar itself and its surroundings going up to half of the adjacent spans) (Figure 3). e average stress that acts over the pillar, as given by , is compared with the pillar strength (σ rp ) increased by a safety factor; the average stress must be less than or equal to the strength of the pillar [27].
In the formulation above, it is assumed that the pillar section (A p ) is subjected to an uniform stress; this assumption takes the following premises: (1) there is flat and extensive topography in relation to depth and (2) there is no redistribution of stresses in the mine confines, and it is assumed that all pillars are under the same axial load [12].

Pillar Strength.
e magnitude of pillar strength is bound to the shape effect [28]; in this understanding, the expressions that give the pillar strength must integrate the geometric variables of the pillar and soon be equated based on the TAT, characterizing a limit restriction.
On the other hand, the requested resistance in the pillar must consider the safety margin (FS pilar ≈ FS) so that the structure does not have to face the ultimate stress regime (yield stress). Given the above, the equation that satisfies this requirement is σ rp � (σ v (At/Ap) + ∆σ)FS pillar .

Advances in Materials Science and Engineering
A lot of research proposed mathematical relations deduced from retroanalysis and by laboratory tests. e two expressions that are broadly accepted and used for pillar design projects are the equation of Obert and Duvall [29] as in the rst following equation and the equation of Salomon and Munro (1967) as the second following equation (14), well known for its application: With σ H/D 1 being the strength of a rock volume of unit dimensions, W p /H p is the relation of pillar slenderness (W p , H p being pillar width and height, respectively). It is worth highlighting that (11) has been countersigned through algebraic proof by Pariseau [24]; such a demonstration must be sustained according to Mohr-Coulomb's rupture criterion.
Regarding pillars with a rectangular support cross section, Mark and Chasse (1997) proposed an expression that involves the pillar edge dimensions, shown in the following equation. is expression is a derivative obtained from the Bieniawski equation:

Strength of Foundations.
e acting load (lithostatic weight) at the top of the pillar is 'passed on' (equilibrium principle) to the oor; consequently, there is a demand for strength from the loaded pillar. It is veri ed that the oor is charged by the 'retransmitted charge' and the charge is due to the mass of the pillar itself. In order to calculate the resistance capacity of the foundation to hold loads due to 'overlying weight', it is advised to use the classic expression applicable to shallow foundations, as proposed by Terzaghi [6]: Q b is the load capacity of the foundation in units of force/area φ is the angle of friction of the foundation rock c is the cohesion of the foundation rock c is the speci c weight of the foundation rock W p is the width of the pillar Sc = 1.0 − 0.4(W P /L P ) and Sq = 1.0 + sinφ(Wp/Lp), with L p being the length of the pillar erefore, the equation that characterizes this restriction is written asσ p ≤ Q b .

Stability of the Layers at Compose the Roof of the Room Spans.
In the estimated width of the spans, two considerations are made [2,26]: (1) the access widths must guarantee the operability of equipment and production systems as well as the supply of services (e.g., air for ventilation) and (2) the structure and stability of the roof must be the same as the other parts.
In another view, the roof can be formed by a slab or a set of strati ed and sedimentary slabs, arranged in layers, and can be determined as a structure exion through concepts of an elastic beam [30]. Goodman uses the mechanic concept of the materials, speci cally the elastic beam theory, to analyze behavior [5]; Pariseau formulates that the slab behavior submitted to pure exion considering the plane strain theory [24]. For that, the restriction of stability cited by the aforementioned authors is written in the following manner [30], in which the analysis of the 'n' slabs willing to ex under an equivalent charge 'Q' is as follows: where Q is a load value (stress) distributed by unit of length; E i and E j are the modulus of elasticity of the rst slab and the upper slabs, respectively; and h i and h j are their respective thickness.
e Q value is a function of the 'n' "stable" slabs: , it follows that the second layer (superior) is loading at rst (bottom) If Q(1) its smaller than Q (2), it follows that the second layer is supporting part of the weight of the rst is analysis must be realized for the entire structure of the 'n' layers that make up the ceiling (1, 2, 3, . . ., n), and, in the absence of tension between the layers Q(n) and Q(n + 1), there is decoupling between the superior and bottom slabs. e equation for the restriction of width (or length) of the span is de ned as follows:  Advances in Materials Science and Engineering Here, h is the size of the layer; σ t is the tensile strength; Q is the stress exerted by the set of slabs that will inflect and FS slab is the security factor for the slabs. By what was already mentioned, the minimum span coherence constraint is written like this: Min(W o , L o ) ≥ L min ; L min is the minimum width, usually defined as a function of the largest equipment that will travel inside the mine, with a gap on each side of the road.

Numerical Modeling and Simulation.
e mining area under study is in a domain presenting 3 types of lithology. e information and input data (geometric arrangement) are provided by the optimization through the algorithm (as described in Section 3.1); therefore, the formulations corresponding to the hydrodynamic theory are placed in the model script.
e information for the modeling is the following: (1) For square section: Items 'a' refer to pillars affected by wave energy, whereas pillars in item 'b' do not consider dynamic energy (for both sections).

Formulation of the Optimal
Pillar strength constraint on the 'x' axis is Pillar constraint on the 'y' axis is Foundation stability constraint is Stability constraint of the slab on the 'x' axis is Stability constraint of the slab on the 'y' axis is Coherence constraint is

Simulation with FLAC 3D (Fast Lagrangian Analysis of Continua in 3D
). It should be highlighted that the Mohr-Coulomb failure criterion was adopted for assessing this model (see Figure 4). e values related to the geometric data are found in the algorithm, as previously described in Section 3.5.7; therefore, the geotechnical information and the wave equations are inserted in the script (Sections 3 and 3.2). e FISH language of FLAC 3D allows creating routines to generate the model interactively; this simplicity facilitates scenario analysis, as well as providing data and obtaining answers quickly. e routine that characterizes the numerical model is established as follows:   programming, as given in Section 3.6. e objective function structure with the restrictive equations represents the sentences that compose the algorithm (model), whose equations depict physical behavior for the mine pillar under static load and dynamic stress. e analytical model described was analyzed in the Mathcad V. 14.0 software and the subsequent evaluation with the FLAC 3D software, for which the variables mentioned above were considered (Sections 3.5.2 and 3.5.7).
On the other hand, the responsible velocity for the dynamic stresses corresponding to the absolute values was analyzed in the centroid of the pillar support section, and    these found velocities were analyzed using the numeric methods given by the Newton theory of uxions and contrasted with the I-Blast V. 8 software, as detailed in the Methodology Section. In Table 5 are the optimal values for the dimensions of the geometric arrangement of the pillar and a discriminated analysis for a pillar with square and rectangular cross section support. As shown in Figure 5, the dynamic stresses due to the seismic wave represent an additional maximum load of 5.99 MPa; therefore, the pillars designed considering these dynamic stresses are more robust; this represents on average 3.7 tons more ore mass contained in each pillar compared to the pillars conventionally designed.

Conclusion
As presented in Table 5, the results corresponding to the recovery factors for the mining arrangement, when considering the wave stress, are relatively inferior in relation to the conventional dimensioning, in which the dynamic stress delivered by wave pulses is neglected. is di erence is due to the fact that the conventional sizing follows simpli ed criteria; hence, its attractive recovery (recovery factor) is detrimental to the pillar quality to respond to stress request. e dimensioned pillars, in a conventional manner, will o er an average recovery of 1.6% more compared to pillar dimensioning, when wave stresses are considered. e percentage of extra recovery represents a value of 3.7 tons on average per pillar; this mineral mass is an apparent gain in the dimensioned pillars.
As shown in Section 4, the dynamic stresses are around 4-5.99 MPa, a value under the mathematical model; although these stresses are much lower than the strength of the pillar, when added to the lithostatic stresses, they could create plasticization processes in several portions in the weakened structure of the pillar. us, compressive stress values were also obtained (Figures 6(a) and 6(b)) under the values obtained from the algorithm using Mathcad (see Table 5) in both sections (square x rectangular).
All velocities were calculated through I-Blast software and hydrodynamic theory application using Mathcad algorithm. Both methodologies demonstrated coherence. Pwave velocity: 424.5 mm/s and 444.64 mm/s (I-Blast and hydrodynamic theory, respectively). Such velocities refer to information on the pillar centroid. Furthermore, velocities answers were contrasted and veri ed with heat map (Figure 7, energy released). Both cases also provided similar information 1 Value of 2 for the ANFO explosive agent.
Pillars with rectangular support cross sections represent the best support option as compared to square section pillars, since the load is better distributed in the former. e wave stress (P-wave) and the dimensions of the widths of the room spans are inversely proportional; these results are in agreement with the velocitythat is contrasted by the simulation and the numerical analysis of the potential velocity. In both cases (support), there is a di erence of 0.73% and 0.66% when the wave e ect is or is not considered severally (see Table 5).
Data Availability e information and data that made this study possible will be made available by the authors upon reasonable request.

Conflicts of Interest
e authors declare that they have no con icts of interest.