#### Abstract

It is so important to consider the passive defense problem in any places there have been attacks by varies kinds of military threats and terrorists. It is certain that social security is related to overcoming on these perils and protection from country. Vital facilities are one of examples that should be protected. Vital facilities include roads, bridges, transmission lines, and telecom and media network. With attention to the intense dependent to export and transmit of oil and gas and with consideration of this point that many places are full of gas and oil resource, the protection of these lines is very important. In recent years, occurrence of varies kinds of terrorist accidents in relation to important structures in all the world causes that the explosion loads have special attention. Explosion can generate much damage with vibration in vast soil media. Thus, it is important to predict the dynamic impact load and its treatment response. With attention to regardable development of numerical methods in recent decades, it is possible to investigate the explosion effects on surface and underground structures. In this research, the newest applied method modeling of the explosion phenomenon has been investigated and comprehensive information has been earned. In this investigation, problem of explosion wave’s propagation effects on buried pipes simulated by ABAQUS/CAE 6.10-1 was studied based on the finite element method. Surface explosion effects on gas buried pipe lines and their dynamic response have been investigated depending on properties and their characteristics. The variation of buried pipe depth effects and variation effects in soil properties around pipe in different cases has been considered, and the results are here. The results showed that in buried pipes under surface explosions, displacements, major stresses, and strains decrease in clay, dense, and loose sands with increase of buried depth. These results obtain that because of increase of closuring of pipes in soil when internal friction angle increases for a kind of soil, the stress on pipe rim will decrease also. It was also observed that the pipe performance in clay and loose sands is better than that in compacted sand, respectively.

#### 1. Introduction

Vital arteries include roads, stairs, tunnels, transmission lines (water, oil, and gas), and communication and media networks. If one of the vital arteries is damaged, malfunctioning, urban activities, or relief work will be paralyzed during the crisis, thus increasing the loss of life and property. If the threatening factors of vital arteries are divided into two categories, natural and human factors, earthquakes, and explosions can be mentioned as examples of them, respectively. Over recent years, various terrorist incidents on important structures worldwide have caused special attention to be paid to explosive devices, in large cities that use underground spaces for highways, tunnels, and underground pipes. Therefore, it is essential to predict the dynamic impact loads and investigate the behavioral response of the structures. The reliable and economical design of structures requires a better understanding of the sophisticated and practical parameters of these structures. Due to the necessity and importance of these vital arteries, and if there is proper information, we will be able to predict the possible effects of the explosion on them according to the specifications of the existing pipelines; in addition, if necessary, we should take steps to cure them. Also, with such information, the design of pipelines in the future can be more precise and safer, and the considerations related to the passive defense of vital arteries can be considered more effective. In recent years, there have been significant advances in developing numerical and laboratory methods for studying such systems. In the present study, comprehensive information will be obtained by examining the latest methods used to model the explosion phenomenon and the responses obtained from underground structures.

##### 1.1. Modeling the Behavior of Materials

The simulation of the nonlinear part of the soil behavior has been performed using the modified Drucker–Prager plastic processor (cap) model. We also need to specify the level of yield, flow law, and hardening law. The rupture level of the Drucker–Prager is obtained from the following equation [1]:

and represent the angle of friction of the material and its bond, respectively, and can be defined as a function of temperature, , and others. *T* and *P* are equivalent to the amount of deviation stress and compressive stress, respectively [2] Table 1.

According to Iran’s 2800 regulations, the selected base soil is type IV, which is considered soft soil. Then, by changing the effective parameters in the interaction and response of the structure, about 50 analyses with different parameters have been performed, which have been then displayed under different diagrams.

In addition, the used mechanical properties for simulation of steel pipeline can be found in Figure 1.

#### 2. Numerical Solution Methods

##### 2.1. Motion Equations

The most common solution to the problem of soil-structure interaction is the analysis based on the substructure method. In this method, the linear problem of soil-structure interaction is separated into a series of simpler subproblems, and then the results are combined using the principle of superposition. In the substructure isolation method [3], the whole soil-structure system, shown in Figure 2, is divided into three substructures. Substructure I includes the free area of the structure, substructure II includes the volume of soil removed, and substructure III includes the surface structure and its foundation. In the substructure isolation method, it is assumed that soil and structure interaction occurs only at the common boundary of the substructures (soil contact surface). The equation of motion of the substructures shown in Figure 2 can be written in matrix form (1). For harmonic excitation, mixed force vectors, and frequency displacement, the equation of motion for the soil-structure interaction system will be (8), when indices I, II, and III are related to the three substructures and indices *i*, , and *s*, respectively, are related to the degrees of freedom corresponding to the nodes at the boundary of soil and structure, removed soil volume, and the upper part of the structure (Figure 3):

##### 2.2. The Site Response

For solving the problem of the construction response, it is necessary to make and solve the problem of eigenvalue for the model. In calculations related to volumetric waves from the submatrices calculated from the characteristics of each layer, they are used to form eigenvalue equations. In Kiran and Manoj’s study [4], based on the construction model with horizontal layers and the assumption of linear deformation changes within each layer, the eigenvalue problem can be divided into two separate algebraic eigenvalue problems. These include the Riley wave, and the other one is for the movement of the Love wave. The equation of the eigenvalue for the motion of the Riley wave can be written in the following matrix form using the isolated soil model. There are two degrees of freedom in this model at each boundary between the two layers; thus, each *n*-layer system will have 2*n* degrees of freedom. In the previous equation, is the mode vibration angular frequency, *K* is the eigenvalue, and is the corresponding 2*n* component eigenvector. The dimensions of the matrices [*A*], [*B*], [*G*], and [*M*] are , and these matrices are obtained by the summation of the matrices related to sublayers. If the depth of the *j* layer from the top is and also the volumetric mass, shear modulus, and its constant are , , and , respectively, these matrices will be as follows:

Eigenvalue equations:

Eigenvalue equation using numerical techniques:

Matrices and are the continuous and concentrated mass matrices, respectively. Using the numerical techniques suggested by Was, equation (5) of the eigenvalue can be obtained. By solving the equation, 2*n* Riley mode, as well as 2*n* wavenumber, is obtained. These values will be used to calculate the conditions of energy-absorbing boundaries in wave deformation motions on the system model screen. Based on the layered soil model, the eigenvalue problem for the motion of the Love wave can be written in the following form. In this waveform, at the boundary of each layer, there is only one degree of freedom. The following matrices , , and are derived from the following matrices.

Shi [5] wrote the equation of motion for SV inclined waves using the *n*-layer soil system. The equation of motion of the soil system, which is affected by SV waves, can be written as (18). Matrices , and are obtained from the summation of the submatrices defined in equations (6) and (7). The matrix is obtained from the summation of the following submatrices, , defined as equation (19). If and are the binding and shear modulus of the *J* layer, respectively, the vector is a two-component vector, which defines the load vector at the bottom of the layer (rock bed). From the solution of equation (7), the displacement vector is obtained. At any distance *x*, the free-field motion can be achieved using equation (7) and is the coefficient of modal participation that is calculated by applying the control motion at a point:

The semi-infinite semispace can be modeled by two methods, variable depth and viscous boundary at the base. Substructure methods are valid only for linear analyses. However, soils exhibit nonlinear behavior depends on strain in response to dynamic loading. Nonlinear soil behavior can be considered using the equivalent linear method proposed by Syed and Idris [6] and others. In this way, soil nonlinear characteristics are estimated by equivalent linear characteristics, including shear modulus and damping coefficient.

##### 2.3. The Impedance

In the substructure isolation method, the equations of motion of the SSI system, including the impedance matrix , are shown in equation (8). In this method, the impedance matrix is calculated only for the boundary nodes (*i* nodes in Figure 3). At each analysis frequency, the impedance matrix is calculated by the dynamic compliance matrix. In three-dimensional problems, estimating the dynamic compliance matrix is reduced to the problem of finding the horizontal layer system response to loading at the boundary of the layers. After calculating the elements of mass and stiffness matrices, the equation of motion becomes (9). When C is the dynamic stiffness matrix , *R* is the impedance matrix of the energy-absorbing boundaries. Indices *c* and *p* are related to the degrees of freedom on the boundary line and the model environment, and and are the corresponding displacement sizes. The displacement of the outer nodes of model C is obtained from equation (10), where the index *m* is related to the Fourier harmonic degree. The vector under the coefficients of the modal participation depends on the 3*n* mode of displacement of the *n*-layer soil system:

Furthermore, matrix is a matrix that is a function of the radial distance from the axis of the model, the eigenvalue, the eigenvector, and the *m* degree function from the second type of Stanković functions [7]. Equation (11) is used with the modal participation vector to calculate the displacement at any point at an *r* radius distance from the model axis. At any frequency of analysis, a dynamic compliance matrix is a matrix for a system with *i* interactional nodes within the free field environment. The direct impedance calculation method demands to calculate the compliance matrix for all the nodes involved. The impedance matrix is then calculated by inverting the dynamic compliance matrix:

##### 2.4. Structural Analysis

In this section, the structural and removed soil characteristics, used in the motion equation coefficient matrix equations (12), are calculated including the components *A*, *B*, and *T*. The removed soil and structure are modeled by standard finite element models, and then their dynamic characteristics are calculated. Within the soil-structure system, in the event of stable dynamic excitation, equation (12) is formed and dissolved at separate selected harmonic frequencies. In the case of harmonic excitation of results, the harmonic mixed conversion acceleration functions indicate the response of the whole system to the harmonic input motion at the control point. Transient motions, such as earthquakes, are analyzed using separated Fourier conversion technique. Utilizing these techniques, the base input motion that specified at N separate points is evenly distributed over the *T* period. The final result is obtained after summation in the following way:

For a single harmonic input, *f* separate values at *L* intervals can be calculated by Fourier inverse converting over *D*. Complete solution demands to form a linear equation system, and then it is solved for all FFT frequencies. This procedure requires considerable time and computational work. A high-level frequency and also an efficient method can be utilized so that the mixed response range, *M*, and *Z* are selected and calculated at several critical frequencies besides response values at other FFT frequencies. It can also be identified by interpolation [8].

##### 2.5. Soil and Structure Interaction

Within studying the interaction problem, the dynamic response of the structure is affected by the interactions of the structure, the foundation, and the surrounding soil. There are two general substructural and direct methods available for analyzing interaction issues [9]. In the substructure method, the entire infinite environment is modeled; after calculating the soil impedance matrix, it will be appended to the dynamic stiffness of the structure. This process will make the problem more straightforward to figure out and more comfortable with calculating. Under the direct method, only a part of the soil environment, which is close to the structure, is modeled [9]. In this method, complex geometric shapes, changes in soil properties, and nonlinear behavior of the environment can be considered. Inside the analysis of the time domain, nonlinear behavior of the soil and structure materials can also be considered.

This study operates the conventional model in defining the interaction of contact surfaces, the Columbus friction model. The model has been applied using the contact element, which will be explained below. The Columbus friction model determines the frictional behavior between contact surfaces utilizing a friction coefficient *µ*. The default coefficient of friction is zero, so the tangential displacement will be zero until the surface tension exceeds a critical value for shear stress, and this depends on the vertical contact compression and is defined according to the following equation:where implies the coefficient of friction and is the contact compression between the two surfaces. This equation determines the extreme value of shear stress for the contact surfaces involved in the collision. As long as the shear stress between them is equal to the frictional shear stress, the contact levels will not slide upon each other, as displayed in the following behavioral diagram (Figure 4).

Soil and structure interaction is significant, particularly for large and massive structures that have been built in soft soil. When affected by a dynamic load, a dynamic reaction of the soil and structure complex signifies a function of the dynamic characteristics, induced forces, and stimuli, and the dynamic model of the system which involved the dynamic model of the structure linking with the dynamic model of the environment. Interaction is a boundary value problem.

Therefore, it requires an infinite environment model; however, infinite modeling of the environment is not achievable. Hence, some terms must be considered that by a partial model of the infinite environment, the state of radiation (i.e., the waves are not reflected from the infinite environment) is satisfied. If the force, displacement, and dynamic stiffness upon the infinite boundary are displayed with *P* (*x*, *t*) and *U* (*x*, *t*) *K* (*x*, *t*), respectively, the following equation is settled:

The purpose of the study soil and structural interaction is to obtain the dynamic stiffness and shape of the wave motion in the meeting points of soil and structure. The direct solution method for studying soil and structure interaction is based on the formulation of finite elements for soil and structure. This feature emphasizes the importance of the method in the topic of the interaction of soil and structure. Another significant point is the ability of the finite element method to consider nonlinear effects, both material and geometric. This makes the method one of the most potent numerical methods in structural analysis.

##### 2.6. Loading Surface Blast

When an explosion occurs under ideal conditions of the theory, in the vicinity of the solid surface, waved shape will be hemispherical (Figure 5). (It should be noted that during an explosion within the air, the wave front is spherical.) Therefore, the energy of the blast wave is concentrated into a smaller area. For this purpose, all the relations obtained for blasting in the air can be utilized by doubling the weight of the explosives. However, in cases where the ground surface is firm and impermeable, laboratory tests are recommended.

The parameters of the blast wavefront are significant. The analytical solution of these parameters was first expressed by Huganiot and Rankine to describe shocks in the ideal gas. These equations are expressed for the velocity of the wavefront and maximum dynamic compression as follows [10]: , and are the static overcompression, atmospheric compression, and velocity of sound, respectively, within the air under atmospheric compression. The scaling distance, *Z*, is defined as follows: [*R*3] is the actual distance from the center of the explosive to the target point, and W is the charge mass in kilograms of TNT. Utilizing *Z* allows us to have a concise and practical expression of the blast wave for various situations [10]:

An equivalent TNT mass is required before the parameters can be extracted for an explosion. There are several ways to express equivalent TNT, but the simplest is the ratio of the specific mass energy of actual explosives to the specific mass energy of TNT. The specific gravity energy of TNT is equal to 6700 kJ/kg [11]. Other parameters of the blast wave include the positive phase duration and the positive impulse, which is the area under the compression-time diagram from time to the end of the positive phase [11]. A standard method for extracting blast wave parameters is to practice graphs provided in some sources, such as TM5-1300 [10]:

There are many relations to maximum overcompression due to conventional (chemical) explosions based on researchers’ explosive applied. Kiran’s relation to the traditional chemical explosion is given [4]: maximum overcompression curves (as well as the curves of other parameters of the blast wave) exist as a function of the distance measured in 1300-500 related to the defense organization. And the maximum overcompression can be obtained by utilizing them (Figure 6):

When an explosive device explodes on the surface of the earth or at a very short distance from it, the blast wave will ideally have a hemispherical wavefront, which is different from the spherical wavefront obtained in an aerial explosion [10]. Therefore, the released energy is applied to a smaller surface, and corrective coefficients must be applied to the air blast equations to obtain the overcompression in the surface explosion. Within 1963, Newmark proposed the following relation for calculating the maximum overcompression in chemical surface explosions:where *P*_{s} is the overcompression, *W* is the equivalent weight of TNT (tons), and *R* is the distance from the center of the explosion on the ground to the desired point in terms of meters (Figure 7).

In studies conducted by various researchers such as Bulson, Bashara, and Henrich on the case of blasts, experiments were reported based on the amount of explosive equivalent TNT; therefore, it seems essential to provide a correlation to equate other explosives with TNT. The mass of explosives can be converted to TNT using the combustion temperature [12]:

In these relations, , , , , and are equivalent to TNT weight, TNT combustion temperature, the explosive combustion temperature, the weight of the explosive, and the energy equivalent of a nuclear explosion (e.g., the amount of TNT *T*, which produces energy equivalent to the energy of a nuclear blast), respectively. Also, is a factor related to the amount of actual energy released as a blast wave (a blast wave of nuclear explosives accounts for 70% of total energy, and the rest of it is wasted in the form of light and radiation) Table 2 [13].

##### 2.7. Study of the Propagation of Blast Waves in the Soil Environment

The compressive waves created due to the explosion propagate in the soil environment, and colliding with the soil causes the soil to crumble and flow [14]. After the explosion, the resulting compressive wave propagates rapidly at a velocity close to sound velocity (about 300 m/s). This wave puts a compression of about 1000 pounds per square inch. Other letters are waves of tension or shock. The boundary impact between two different materials on how the wave propagates can be significant. Estimating the relative amplitude and directions in which the waves are reflected or refracted at the boundaries is available, utilizing the relations obtained from the theory of elasticity. Using the following equation also holding wave velocity and mass density, we can obtain the relation between and where *c*, *ρ*, and *A*_{i} are the velocity of the wave, density, and the amplitude of the initial stress wave, respectively. Moreover, is the amplitudes of the reflected stress wave:

The behavior of soil under dynamic loading is one of the most exciting topics for engineers in mines, buildings, and defense structures. Generally, the soil is a three-phase mixture of solid mineral particles, water, and the air. Connected or separated solid particles form the soil skeleton. Water and air are located among the cavities between solid particles of soil. Soil is the three-phase mixture; that is why it is challenging to predict the deformation of soils; not only is the skeletal structure complex, but the properties of each component are fundamentally different. The different deformation characteristics of each phase and soil-structure make the soil deformation mechanism strongly dependent on the ratio of components in the soil and loading conditions. As a result, the mechanism of deformation of saturated and unsaturated soils is different, and such a difference in dynamic loading will be much more noticeable than static loading.

Both the mechanisms of skeletal deformation and the deformation of all soil phases affect simultaneously. However, in different loading stages and depending on the different ratios of components in the soil, one of the mechanisms is dominant so that the other mechanisms can be ignored [15]. For dry soils, which include a large amount of air and have a low percentage of water, the first mechanism prevails when subjected to static loading or slow dynamic loading. The initial compaction of the air is extreme. However, by increasing compression, the bonds between the soil particles are deformed; furthermore, the soil skeleton is damaged. Such soil is considered dense soil. Then, the second mechanism becomes more critical while the first mechanism gradually disappears [16].

The extraordinary dependence of soil properties on loading conditions makes it tough to develop a unique model for the deformation of soils under dynamic loading, particularly for high explosion rate loading, which varies significantly in different mass charge situations under loading conditions. Within an explosion, the soil close to the mass charge is deeply compacted. Hence, a high-stress wave is generated and propagates outward into the soil. As the distance from the charge increases, the stress wave decreases rapidly, and as a result, the soil compaction decreases. Therefore, in areas close to the charge, the second mechanism is predominant. In contrast, with increasing distance from the mass charge, the first mechanism becomes more and more significant and finally dominates the deformation of the soil. The extent of the nearby area depends on the fundamental characteristics of the soil. Soil deformation mechanisms under static and transient loading have been studied over the years, and many models have been proposed to describe soil behavior under such loads.

However, insufficient research studies have been conducted on soil behavior under transient loads, such as explosion loads, due to these complexities. Modeling the soil response to an extremely variable pressure requires a robust soil model because new behavioral models typically cover a small range of deformations. Besides, in such high-stress conditions, in the vicinity of the charge, it is illogical to consider the rigid materials for the soil particles considered in ordinary soil dynamics. Instead, the solid phase change should be regarded in the soil behavior model. Based on Kandaur’s conceptual analysis of soil deformation, a three-phase behavioral model for shock loading was developed by Wang (2003) to provide all of these requirements. In this model, the soil is considered a three-phase system consisting of solid particles, water, and air, in which the climate fills the solid particles form the soil skeleton, and the space between the particles is filled by water and air. The solid phase is considered plastic. Since the blast load duration is concise and there is not enough time for the water and air to escape from the soil particles, the relative displacement of the air-water and soil skeleton is neglected. Therefore, the model is suitable for describing soil behavior under normal blast loading [17].

In Figure 8, elements *A*, *B*, and *C* represent the deformation of the solid particles, water, and air, respectively, and elements *D* and *E* expose the friction and adhesive strength between the solid particles. The bonds between solid particles have been represented by a series of strings. Elements *A*, *B*, and *C* reflect the second deformation mechanism, and elements *D* and *E* are related to the first deformation mechanism. Based on the preceding modeling, a soil model for shock loading is obtained, which includes the following: equation of state (EOS), stress-strain relation, resistance model, and failure model for soil skeleton [18]. Different techniques are available in the technical literature to model a blast wave on buried structures. These impacts are in the form of displacement waves, acceleration, or compression. Naturally, the source of all methods is the same and is accomplished by solving the shock wave propagation equations in the soil environment. Of course, these relations are also interchangeable. The maximum amount of particle displacement (*x*) in terms of meters, at a distance of *R* (meters) from the source of the explosion, entirely or relatively buried, can be estimated from relation (22), where is the mass of the explosive material in kilograms, *c* is the seismic wave velocity in meters per second, is the coupling coefficient, calculated from Figure 9, and *n* is the dimensionless reduction coefficient; from Table 3, it is obtained that [19]

**(a)**

**(b)**

Also, the maximum velocity of particles in meters per second at a distance of *R* (meters) from the explosion source was calculated from equation (23). Moreover, the maximum pressure of the free field was determined based on equation (24). Ρ is the density of the soil (kg/m^{3}), C is the load wave velocity (m/s), and P_0 is in pascals. The loading wave velocity *C* depends on the velocity of the seismic wave and the maximum particle velocity. *C* is high at short distances due to the extraordinary velocity of the particles, but its value decreases as the seismic wave rises with increasing distance. Also, its quantity is never less than the seismic wave velocity [20]:

Also, the time of continuous explosion compression on underground structures is obtained from equation (26) in which is the specific momentum of the free zone, and placing the quantities of displacement and velocity of the particles in the above relation reaches the following relation for the time of total continuity equation (27). Inside the above relation, transpires in seconds, *R* stands in meters, and *C* is in meters per second. However, to make modeling the analysis more comfortable, the reflection pressure drop has been omitted. The relationship of the blast pressure on the pipe in time has become

#### 3. Modeling and Conclusion

##### 3.1. Sensitivity Analysis to Control Model Dimensions

Since the infinite soil environment is modeled in a finite form, the question arises about the dimensions of the soil model. Soil model dimensions should be chosen so that, while small, it has the least impact on the results; to minimize the time of analysis to get accurate results. Because the larger the model, the more accurate the results, but due to the larger dimensions of the model, the time of analysis is extended once the analysis of the model may take days. For this purpose, three models were made according to Table 4.

The above diagrams show the acceleration and wave velocity at the farthest node of the defined models (corner points) for the same loading and boundary conditions. According to this, it can be concluded based on obtained results (i.e., soil with dimensions of 50 100 100 × 25 m). The velocity and acceleration of the blast wave resulting from the explosion reaching the farthest node tend to zero. As the dimensions grow, only the analysis time increases, but there is no change in the results and outputs. Thus, the model that was considered for this project, according to the sensitivity analysis, is model 3, i.e., soil with dimensions of 100 × 50 × 25 meters, Figures 10 and 11.

##### 3.2. Finite Element Model

The geometry of the model consists of two parts, soil and pipe. The dimensions of each are as follows: soil in the form of a cube with dimensions of 100 × 50 × 25 meters, and a steel pipe with a diameter of 1 meter, a length of 100 meters, and a thickness of 3/14 mm. After determining the geometry of the model, the specifications of the materials used should be defined. Soil and steel specifications are shown in Tables 5 and 6, respectively. For steel materials, the yield stress is equal to 0.5% strain, and ultimate tensile strength is equal to a 19% strain. The Drucker–Prager plastic model was also used to define soil behavior. All members are homogeneous. For pipe modeling, shell elements have been used, and for soil modeling, solid elements have been used.

Two methods can be applied in modeling soil-structure interaction. One is soil replacement modeling with springs in three directions (two directions perpendicular to the pipe and one direction tangent to the pipe). Another way is to build a full model, including soil, pipes, and their interaction in three dimensions. In this study, due to the closest approximation to reality, the second method was utilized. In this modeling, the interaction between soil and pipe elements within the software is modeled on surface-to-surface mode, between the outer surface of the pipe and the inner surface of the soil and by the penalty method. In this method, the normal contact of two levels with the hard contact method is considered to prevent the surfaces from sinking into each other. Also, the tangential interaction of the two surfaces with a friction coefficient of 0.5 that was suggested by Vasouras and Kolbadi [21, 22] is regarded in the analyzes. This number is a suitable value for the interaction of soil with metal structures. However, in a part of the study, the effect of the quantity of this coefficient of friction on the behavior of the pipe was investigated, and the results are presented below. Also, at this level, damping contact was considered with a coefficient of 0.05.

Since the peripheral soil is semi-infinite, it must be modeled so that all its properties, including its being semi-infinity, are considered. Based on the previously performed sensitivity analysis (Figure 12), we saw that the velocity and acceleration of the blast wave that reaches the farthest node tend to zero [23, 24]. When the wave reaches the farthest nodes from all sides, it is damped. Therefore, there is no need to use energy absorbing and damping boundaries on the edges of the soil environment. Thus, symmetry conditions can be utilized on the lateral surfaces of the soil model without any challenges [25]. Of course, it should be noted that these conditions, i.e., not using adsorbent and damping boundaries on the lateral surfaces of the soil environment, are appointed just because of the type of loading, i.e., explosion, as well as the dimensions considered for this model, and cannot be generalized to other situations [26]. For example, the conditions for earthquake loading are different. If the adsorbent and damping boundaries are not used on the lateral surfaces of the soil environment, the wave reflected from the structure by these boundaries is reflected in the environment and causes errors in the results [27]. Therefore, absorbing boundaries must be applied to prevent the reflection of the waves. The soil floor also has encastre boundary conditions.

**(a)**

**(b)**

**(c)**

##### 3.3. Loading

The loads applied to the pipe include the load due to the blast wave, the internal pressure of the pipe, and the overburden soil weight of the pipe. The blast load is applied as a compression-time series for an equivalent TNT mass of 10, 20, and 30 kg at a point on the soil based on the equations presented equation (22). The place of the surface load application was a hole dug in the shape of a hemisphere with a radius of 1 meter on the soil surface. The internal pressure of the pipe was also applied as a uniform pressure on the inner surface of the pipe. Nine-node cubic elements were used in meshing the soil volume, which is smaller in sensitive areas in Figure 13.

**(a)**

**(b)**

**(c)**

##### 3.4. Type of Analysis

Knowing the nature of each of these methods and gaining experience will be simple to choose the solution. In the explicit method, the results at each moment (*x*_{n} + 1) are obtained directly from the results before (*x*_{n}). In this way, new situations are calculated by considering the velocity and acceleration of the elements at the moment *n* and placing them in the following equation:

On the other hand, in the implicit method, the size of Δ*t* does not matter much. In this method, the equations of position , velocity , and acceleration are solved simultaneously by iterative methods. In this way, there is no trace of speed and acceleration in the nodes. In high-speed physical phenomena, such as blast forces or collision forces that inflict a significant load on a structure in a short period, solver convergence is implicitly impossible. Besides, observation of the reaction of the structure in small time intervals is considered. The following is the process of spreading the stress wave propagation caused by the explosion in the soil and pipe under the blast load in Figure 14.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

#### 4. Review of Analysis Results

To analyze the sensitivity of the parameters affecting the performance of the pipe, the models and the results of their analysis are examined separately. The parameters that underlie the judgment are stress, strain, and displacement. The stress understudy is phonemic stress obtained from the square root of the sum of the principal stresses. Displacement also means the maximum absolute displacement of one of the pipe nodes.

##### 4.1. The Impact of Pipe Burial Depth

Blast loading was performed on three types of soil (loose sand, hard sand, and clay) and its impact at three different depths. Depth of burial pipe tension, strain, and maximum displacement of the pipe is significantly reduced. The effect of increasing the burial depth of pipes on their behavior against explosion can be mentioned for various reasons, such as increasing the pipe blockage in the soil, increasing the distance of the pipe to the explosion place, and increasing the stiffness of the soil (Table 7). Therefore, it seems that the performance of buried pipes is better in clay soils and loose sands, respectively, especially with increasing pipe burial depth; this is evident. The phenomenon due to the higher seismic wave velocity in hard sandy soils is directly related to the modulus of elasticity of the soil and increases the blast wave maximum pressure (Figures 15–17).

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

##### 4.2. The Impact of Explosive Weight

Based on the relations presented to calculate the mass charge, the weight of the mass charge directly affects the loads caused by the explosion. Therefore, according to the results obtained from the analysis of the models given in Figures 18–20, it can be stated that increasing the weight of charge mass, strain and maximum displacement of the pipe will increase. Here, too, the performance of pipes buried in loose soils and sands, particularly clay soils, seems to be better (Table 8).

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

##### 4.3. The Impact of Blast Distance from the Pipe Axis

Based on the results obtained from the analysis of the models given in Figure 21, it can be said that by increasing the blast distance, strain and maximum displacement of the pipe will be significantly reduced. This decreasing trend tends to zero after specific distances (depending on the conditions), where it occurs after a distance of 5 meters (Table 9).

**(a)**

**(b)**

**(c)**

##### 4.4. The Impact of Internal Pipe Pressure

According to the results of the analysis of the models, which are given in Table 10 and 11, it can be said that the internal pressure has little impact on the pipe response practically.

##### 4.5. The Impact of Friction Coefficient between Soil and Pipe

According to the analysis results, which are given within Tables 12 and 13, it can be recognized that practically, the friction coefficient between the soil and the pipe has a negligible impact on the pipe response, due to the propagation of the wave, which is perpendicular to the axis of the pipe. Hence, the impact it has on the pipe is more inclined to normal type than the contact. So, the normal interaction of soil and pipes is much more critical than their contact interaction.

##### 4.6. Pipe Diameter Impact

As expected, the larger the diameter of the pipe, the better the pipe performance as the cross-sectional moment of inertia increases based on Figure 22 and Table 14.

**(a)**

**(b)**

**(c)**

##### 4.7. Pipe Thickness Impact

As expected, the thicker the pipe, the better the performance of the pipe as the moment of inertia increases as can be found in Figure 23 and Table 15.

**(a)**

**(b)**

**(c)**

##### 4.8. The Impact of Soil Internal Friction Coefficient

Based on the results obtained from the analysis of the models given in Figure 24, it can be stated that by increasing the internal friction coefficient of the soil, the maximum strain and displacement of the pipe decrease (Table 16).

**(a)**

**(b)**

**(c)**

##### 4.9. Soil Specific Gravity Impact

Based on the results of the analyzes given in Figure 25, it can be understood that by increasing the specific gravity of the soil, the maximum stress and strain of the pipe are significantly reduced (Table 17). However, there is no noticeable change in maximum pipe displacement.

**(a)**

**(b)**

**(c)**

#### 5. Conclusions

Summary of the achievements of studying the performance of buried pipelines under the effect of explosive loading is presented in this section. Then, based on the author’s experience, to reduce ambiguities, conduct experiments, furthermore analyze some of the obtained results, and also present a series of suggestions and approaches for using the results of this research, some ideas are mentioned in the suggestions section.

As can be seen, from the problem under study, first, a suitable model was constructed in three-dimensional form using the finite element method. Then, in the analysis, by changing the parameters that seemed to be important in the behavior of the pipe under the effect of the surface explosion, we examined the sensitivity of the pipe response to the change of the desired parameter. After analyzing and extracting the results, these results were presented in charts, which were based on judgments and conclusions about pipe behavior. At the end of the analysis used in this study, results were obtained to understand better the behavior and optimal design of buried pipelines under the impact of surface explosions, which are mentioned as follows:(1)Blast wave arrival velocity depends on the velocity of the soil environment compressive seismic wave. The arrival time of the wave in clay soils is longer than other soils. The blast wave response time to the structure is a significant factor for the response of underground structures such as buried pipes.(2)Due to the low mass of the pipe and the high dependence of its behavior on the strains that occur in the surrounding soil, it can be stated that the entering velocity spectrum of the pipe is much more critical than the acceleration spectrum.(3)As the pipe depth increases, the maximum stress, strain, and displacement of the pipe decrease almost linearly. Buried pipes seem to perform better in clay soils and loose sands, respectively, especially while the burial depth of the pipe increases. The reason is that this phenomenon is the higher seismic wave velocity in stiff sandy soils, which is directly related to the modulus of soil elasticity and increases the maximum blast wave compression.(4)By increasing the blast distance from the tension pipe axis, the maximum strain and displacement of the pipe are significantly reduced. This decreasing trend tends to zero after specific distances (according to environmental conditions, pipes, and explosions); here, it has occurred after 5 meters.(5)The internal pressure of the pipe and the coefficient of interaction friction between the soil and the pipe have a negligible effect on the response of the pipe.(6)The larger the diameter and thickness of the pipe, the better the performance of the pipe with increasing the cross-sectional moment of inertia. The relation between the results and the conversions in these parameters is approximately linear.(7)By increasing the coefficient of internal friction of the soil, the maximum strain and displacement of the pipe decrease. The effectiveness of this parameter is significant compared to other parameters.(8)While the specific gravity of the soil increases, the maximum pipe stress and strain decrease considerably, but there is no noticeable change in the maximum pipe displacement.(9)According to the previous results, in the case of buried shell structures, such as pipes, it can be understood that their performance is more affected by the deformation of the structure (such as pipe bending) than the effects of soil-structure interaction.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon reasonable request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

The authors of this article are grateful for the valuable guidance of Mr. Abbas Aghasi.