# Research on Runway Blockage Probability through Simulation

Authors: Wan-Li Wei; Ye-Qing Chen; Zhen-Qing Wang; Jian-Hui Wang; Xin-Jian Ren; Shang Ma
Journal: Shock and Vibration
 (2022)

---

## Abstract

The damage and blockade of airport runways are important objectives in modern warfare. The minimum takeoff and landing window is regarded as the damage criterion of an airport runway. On the basis of Mathematica, a simulation is performed by judging whether the inequality group has a solution domain. An improved algorithm for judging whether an airport has a minimum takeoff and landing window is proposed. The algorithm is used to calculate and analyze the blockade probability of an ordinary missile to an airport and the influence of the performance of a guided submunition on the blockade of an airport runway. Results show that increasing the number of submunition, reducing the circular error probable of submunition, and increasing the number of subprojectiles can increase the blocking capability of submunition, while increasing the dispersion radius of subprojectiles will first increase the blocking probability and then decrease it. Through simulation, the inflection point of the blocking efficiency of submunitions is found. This article improves the judgment formula of takeoff and landing window, Compared with the window scanning method, the proposed method solves the problems of error caused by excessively large scanning step and extremely long calculation time caused by considerably small scanning step.


---

## Body

## 1. Introduction

In modern wars, the guarantee of air power plays an important role in national defense. As an important part of an airport, a runway is more obvious than other targets, which makes it the most vulnerable area. Another important issue is the availability of an area for aircraft to take off after the runway is attacked. Many criteria exist for evaluating missile damage to a target, which are generally based on the damage mechanism of a shock wave and debris caused by the missile to the target; the corresponding damage model is established, and the damage level is determined in accordance with damage criteria [1].For the airport runway, the crater caused by the missile is the most important factor. Airport runways use the minimum takeoff and landing window as a criterion to judge whether an airport is blocked. The takeoff and landing window is the smallest intact rectangular area required for a fixed-wing aircraft to take off and land, which means that there should be no craters in this area. After an attack, if no area can be found within a certain period to meet the requirements for takeoff, then the runway is considered to have been blocked successfully. In the past, most research objects were ordinary missiles [2, 3]. With the development of technology, the damage effect of a guided submunition is enhanced [4], and the blockade efficiency is increased. Many factors influence the damage effect of a guided submunition on a target [5, 6], such as the shooting accuracy of the submunition, the dispersion radius of subprojectiles, and the size of a crater caused by subprojectiles. In particular, the dispersion radius of subprojectiles and shooting accuracy of submunition have a great impact on the airport blockade probability. Because airports have many sizes, the runway parameters also should be considered when analyzing the influencing factors of a guided submunition on airport damage [7]. At present, when analyzing the damage of a guided submunition to an airport, the Monte Carlo method [8] is often used to simulate and calculate the impact point [9, 10]. Three methods are adopted to judge whether a minimum takeoff and landing window exists; they are the pixel simulation method [11], random sampling method, and region search method [12]. Prabhakar [13] developed software to evaluate the damage of a guided submunition to an airport and other facilities. They considered many factors that affect the damage of a guided submunition to a target, but their analysis of the standard of airport damage and influencing factors is insufficient. Xiaomei Wang [11] used the area search method to study the impact point of a submunition and the blockade probability of an airport runway, but they did not research the situation in which a plane takes off with an angle to the runway. Junfei Xu [12] used inequality to judge whether the impact point of subprojectiles has an influence on the takeoff and landing window. However, an analysis of the probability of blocking the runway is lacking, especially about the dispersion radius of subprojectiles, the number of subprojectiles, the size of the runway, and the takeoff and landing window. The inequality they used can be further optimized, and the core method is still a region search method. At present, most region search methods have the problems of an excessively large scanning step size resulting in an error and an extremely small scanning step size resulting in long calculation time.Because the past research about the existing form of takeoff and landing window has imperfection and the judgment formula also has some optimization space, this paper optimizes the judgment formula of takeoff and landing window. In order to solve the inaccurate calculation problem of window scanning method, in this article, finding the takeoff and landing window is transformed into judging whether the inequality system has a solution domain. This research solves the problems of region search methods, and the solution domain diagram is processed to quickly judge whether there is a solution domain by judging the image quality. In this way, it is more convenient to simulate the blocking probability. The calculation of blocking probability is no longer affected by scanning step size. The simulation software is Mathematica [14]. The influence of various factors on the blockade probability is systematically analyzed, and a formula is provided to judge rapidly whether an airport can be blocked. This method can effectively simulate and analyze the blockade probability of common missiles, guided submunitions, and other weapons, which lays a foundation for the damage simulation research of other weapons in the future.


## 2. Blockade Criteria for Airport Runways

As the most important infrastructure of an airport, a runway undertakes the tasks of taxiing, taking off and landing, and parking. An airport runway is generally rectangular, and different layout designs are adopted in accordance with different types of airports. Airport runways are composed of main and auxiliary runways, and the length and width of runways should meet the requirements of safe takeoff and taxiing. The main factors affecting the probability of runway blockade are the characteristics of airport runways and damage weapons. The characteristics of airport runways include the strength and thickness of runway surfaces and the size of runways [15]. When the damage weapon is shrapnel, the characteristics of the weapon include the number of submunition, the crater area, the shooting accuracy of the submunition (circular error probable (CEP)), and the dispersion radius of subprojectiles. Generally, the takeoff and landing window is used as the criterion of airport damage, which is related to such factors as the taxiing distance and the track width of the aircraft.Assuming that the runway length isl, the runway width is h, and the crater radius is r, the required length and width of the takeoff and landing window are m and n, respectively. Figure 1 shows the impact range of an effective antirunway projectile. S is the takeoff and landing window required by aircraft. The width of regions ②, ④, ⑤, and ⑦ is r, and regions ①, ③, ⑥, and ⑧ are a quarter circle of radius r. When a missile hits these areas, the takeoff and landing window will be affected by the projectile. At this time, the area shown by S will fail. Therefore, judging whether a runway is blocked can be realized by judging whether area S exists on the runway that makes all missiles hit outside the area as shown in Figure 1.Figure 1
Impact range of an effective antirunway projectile.After a runway is damaged, an aircraft can take off via two ways: one is to take off along the runway length, and the other is to take off at a certain angle with the runway length, as shown in Figure2. Thus, research on damage criteria should consider the existence of a certain angle between the takeoff and landing window and the runway length.Figure 2
Two ways of takeoff and landing.


## 3. Simulation Principle

The established coordinate system is shown in Figure3, and the coordinate axis is set at the end of the runway. The coordinates of projectile Pi are (ei, fi). In accordance with the above analysis, when the projectile falls on the rectangles EFIJ and KLGH and four quarter circles centered on A, B, C, and D, a takeoff and landing window exists.Figure 3
Verification graph of the minimum takeoff and landing window.In accordance with the vector rule, when the projectile falls outside the rectangle EFIJ, the condition is transformed into the system of inequalities, which is shown as follows:(1)JI⟶×JP⟶⋅FE⟶×FP⟶≤0orEJ⟶×EP⟶⋅IF⟶×IP⟶≤0.When the projectile falls outside the rectangle KLGH,(2)KH⟶×KP⟶⋅GL⟶×GP⟶≤0orLK⟶×LP⟶⋅HG⟶×HP⟶≤0.For the convenience of calculation, when a projectile falls on the outside of four quarter circles, it can be regarded as falling on the outside of four circles withA, B, C, and D as the center. The redundant part coincides with the rectangles EFIJ and KLGH, and it does not affect the result. In this case, the system of inequalities is(3)PA⟶≥r,PB⟶≥r,PC⟶≥r,PD⟶≥r.At the same time, the end pointsA, B, C, and D of the takeoff and landing window should be within the runway. Combined with the above analysis, let the coordinates of point A be (x, y) and the angle between AB and X axis be z such that the coordinates of the other points are B (x + m·cosz, y + m·sinz), C (x + m·cosz-n·sinz, y + m·sinz + n·cosz), and D (x-n·sinz, y + n·cosz). Each end point should satisfy the system of inequalities, which is shown as follows:(4)0≤x≤l,0≤y≤h,0≤x+m⋅cosz≤l,0≤y+m⋅sinz≤h,0≤x+m⋅cosz−n⋅sinz≤l,0≤y+m⋅sinz+n⋅cosz≤h,0≤x−n⋅sinz≤l,0≤y+n⋅cosz≤h.We transform the problem into whether the solution domains ofx, y, and z satisfying the abovementioned inequality system exist. The data are integrated into the inequality system to obtain the solution domain satisfying the inequality system.When the inequality system has a solution domain, the runway blockade fails. If no rectangular region satisfies the condition, then the runway blockade is successful. Figure4 shows the flow of analyzing whether the takeoff and landing window exists. Different numbers of inequalities will be generated in accordance with the number of projectiles, which can be processed using Mathematica.Figure 4
The flow of analyzing whether the takeoff and landing window exists.Mathematica is scientific computing software, which combines numerical and symbolic computing engine, graphics system, programming language, text system, and advanced connection with other applications [14]. In accordance with the actual situation, the airport runway and missile data are substituted into the inequality, then the solution domain of the inequality system is obtained. The solution domain of the inequality system is output in the form of image by using RegionPlot3D code to determine rapidly whether the runway is blocked successfully and save analysis time, as shown in Figure 5.Figure 5
Solution domain graph of the inequality system. (a) Existence of a solution domain and (b) no solution domain.
(a)(b)When the blocking probability needs to be simulated for many times, ParallelTable code is used to complete the simulation work. An image recognition method is used for statistics to improve the work efficiency. The axis is removed from the solution domain diagram and processed into the form shown in Figure6. The quality of each image is judged using the ImageMeasurements code. When the system of inequalities has a solution domain, that is, the runway has a takeoff and landing window, the image quality is not 1. When the system of inequalities does not have a solution domain, that is, no takeoff and landing window exists in the runway, the image is blocked and the image quality is 1. The runway blockade probability can be rapidly calculated through the image quality.Figure 6
Solution domain graph of the inequality system after processing. (a) Existence of a solution domain and (b) no solution domain.
(a)(b)


## 4. Simulation of Damage to Airport Runways by Conventional Missiles

Tables1 and 2 show the typical data of airport runways and the takeoff and landing windows for some aircraft.Table 1
Typical data of foreign airport runways.
Runway categoryAircraft typeRunway lengthl (m)Runway widthh (m)Level 1Heavy bomber2500–500060–100Level 2Medium bomber250045–60Fighter-bomber200045Level 3Fighter1800–200040Level 4Trainer<180030Table 2
Takeoff and landing windows for some airplanes.
Aircraft typeIDFMirage 2000-5F-16Length of a takeoff and landing window m (m)800457350Width of a takeoff and landing window n (m)202020An airport adopts a runway with a length of 3000 m and a width of 60 m, and the minimum takeoff and landing window is 400 m in length and 20 m in width. Given that the crater radius is related to the charge quantity of a missile and the strength and thickness of each layer of an airport, the simulated crater radius is randomly distributed in the range of 4–6 m. The missiles are 10, 20, and 30, and the drop points are randomly and evenly distributed within the airport runway. The simulation times are 100.The missile landing point (ei, fi) is(5)ei=l⋅ui,fi=h⋅ui,where u is a uniformly distributed random number between (0, 1).Table3 presents that the blockage probability of ordinary missiles to an airport is low, and a large number of missiles are needed to complete the blockade of the airport. The simulation results show that the uniform distribution of impact points on the airport is an ideal state, but, in fact, it is difficult to achieve due to the problem of missile accuracy. Therefore, to save the economic cost and achieve the tactical purpose rapidly and efficiently, countries began to develop submunition and other weapons to block airports.Table 3
Simulation result.
Number of missiles102030Blockage probability3 (%)39 (%)88 (%)


## 5. Damage of Submunitions to Airport Runways

For submunitions, the main factors affecting the probability of blocking runways are the characteristics of runways and submunitions. The factors that influence the blocking probability of submunitions include the number of submunitions, the crater area, the shooting accuracy of submunitions, and the dispersion radius of subprojectiles [15]. Generally, CEP is used to express the accuracy, which is the sum of the systematic error and the dispersion error of the impact point. CEP is defined as [12](6)P=12πσxσy∫∫x2+y2≤R12exp−12x−μx2σx2−y−μy2σy2dxdy=0.5.In the formula,P is the hit probability of shrapnel, R1 is the radius of the dispersion circle; μx and μy are the estimation of systematic error; and σx and σy are the estimation of random error.The actual coordinates of a submunition are as follows [16]:(7)X=Xm+u⋅CEP/1.1774,Y=Ym+v⋅CEP/1.1774,where (x, y) are the actual impact points of the submunition, u and v are two random numbers that are independently subject to normal distribution in (0, 1), and (Xm, Ym) are the aiming points of the submunition.The coordinates of a subprojectile (ei, fi) are as follows:(8)ei=X+R⋅j⋅cosk⋅2π,fi=Y+R⋅j⋅sink⋅2π,where j and k are uniformly distributed random numbers between (0, 1) and R is the scattering radius.As shown in Figure7, the aiming points of the submunition are evenly distributed along the center line of the runway, d is the distance among the aiming points of the submunition, and the distance between the aiming points of the submunition at both ends of the runway and the boundary of the runway is d/2.Figure 7
Aiming points of a submunition.Therefore, the coordinates of the preaiming points of the submunition should be(9)Xi=d2+i−1⋅d,Yi=h2,d=lS.In the formula,l is the length of the runway, h is the width of the runway, and S is the number of submunitions. Figure 8 shows the simulation of the drop point of subprojectiles.Figure 8
Simulation of the drop point of subprojectiles.Some subprojectiles will be dispersed out of the runway during the process of submunition dispersion, and these subprojectiles will not affect the runway. A schematic of the impact point of subprojectiles is shown in Figure9. The subprojectiles should be filtered to accelerate the calculation and improve the simulation efficiency. The coordinate system in Figure 3 is used when the subprojectile coordinates (ei, fi) satisfy the following formula:(10)−r≤ei≤l+r,−r≤fi≤h+r.Figure 9
Schematic of the impact point of subprojectiles.At this time, the subprojectiles will damage the runway, which are called effective subprojectiles. All subprojectiles are screened, and all valid ones are recorded. The coordinates of effective subprojectiles are integrated into the inequality system to judge whether a solution domain exists, then the blockage probability is obtained.In accordance with the analysis, the simulation process is as follows:(1)
The runway data and the minimum takeoff and landing window size are input(2)
The number of submunition, CEP, the dispersion radius of subprojectiles, and the size of the subprojectile crater are input(3)
The impact points of submunitions and subprojectiles are generated, and the invalid subprojectiles are filtered out(4)
A solution domain graph of the inequality system is generated, the existence of a solution domain is judged using image quality, and the blockage probability is recordedThe specific simulation process is shown in Figure10.Figure 10
Simulation flowchart.The same airport and submunition data as those used by Xiaomei Wang [11] are used to verify the accuracy of the simulation. Suppose that the CEP of a submunition is 100 m, the dispersion radius R is 150 m, every submunition contains 120 subprojectiles, and the radius of the subprojectile crater is 2 m. The submunition dispersion method presents a uniform distribution in the circle. A runway with a length of 3000 m and a width of 60 m is adopted for the target airport. The minimum takeoff and landing window has a size of 400 m × 20 m. The simulation times is 1000.Table4 shows the comparison of simulation results. The blockage probability obtained using the method in this article is generally smaller than that obtained using the method of Xiaomei Wang [11]. This paper considers the situation in which the takeoff and landing window has a certain angle with the runway, and no omission exists due to the problem of scanning step size. Consequently, the result calculated in this paper is smaller than that in Reference [11].Table 4
Comparison of simulation results.
Simulation methodNumber of submunitions5 (%)6 (%)7 (%)8 (%)9 (%)Method of this article1.610.120.131.242.3Method of Xiaomei Wang [11]1.81125.73646Figure11 clearly shows that in the case of five or six submunitions damaging the runway, the difference in blockage probability between the two methods is minimal. The main reason is that the distance among the aiming points is excessively large and the blockage probability is extremely small. When the number of submunitions is increased to 7, 8, and 9, the difference in blockage probability between the two methods becomes evident. The reason is that with the decrease in the distance among the aiming points, the blockage probability decreases due to the existence of an oblique takeoff and landing window. Hence, the gap between the two methods increases.Figure 11
Comparison of simulation results.


### 5.1. Simulation Analysis of Submunition Blockage Probability

#### 5.1.1. Influence of the Number of Submunitions on the Blockage Probability

To study the influence of the number of submunitions on the airport blockage probability, in accordance with the above data, the influence of the number of submunitions on the airport blockage probability is studied by changing the number of submunitions. The simulation results are shown in Figure12.Figure 12
Influence of the number of submunitions on the blockage probability.From the figure, increasing the number of submunitions can increase the blockage probability, which is approximately linear at the beginning. When the number of submunitions increases to a certain extent, the effect of increasing the number of submunitions on improving the blockage probability decreases. The turning point is at the place where the number of submunitions is 11. At this time, the dispersion radius of subprojectiles is 150 m, which is close to 136 m half of the distance among the preaiming points of the submunitions. The coverage area of subprojectiles has overlapped. This condition is speculated as the reason for the decline in blockage efficiency, which will be discussed later.


#### 5.1.2. Influence of the Takeoff and Landing Window Length and Runway Length on the Blockage Probability

Owing to aircraft performance and other reasons, different aircraft need diverse takeoff and landing window lengths. Simulation analysis is carried out to study the influence of the takeoff and landing window length on the blockage probability. The runway length is 3000 m, and the number of submunitions is 6. Other data are the same as above. The simulation results are shown in Figure13.Figure 13
Influence of different takeoff and landing window lengths on the blockage probability.When the length of the takeoff and landing window increases, the blockage probability of the same number of submunitions increases. When the length of the takeoff and landing window is 250–400 m, the blocking probability increases more gently than that when the length is 400–800 m. Given that few subprojectiles landed on the runway and the probability of blockade is small, the difference caused by the change in takeoff and landing windows is minimal.Although the blockage probability of the same number of submunitions under different runway lengths is different, when the number of submunitions is changed such that the preaiming points of the submunitions are the same, the blockade probability is the same. To verify the analysis, the runway length is set to 4000 m, and other data are the same as above. The simulation results are shown in Table5.Table 5
Blockage probability simulation at 4000 m runway length.
Number of submunitions (pieces)89101112Blockage probability (%)10.121.223.435.938.8The simulation results are transformed into variables with the preaiming point spacing of submunitions. The data in Figure11 are processed and plotted in Figure 14.Figure 14
Influence of the preaiming point spacing of submunitions on the blockage probability in different runway lengths.When the preaiming points of submunitions are evenly distributed along the runway axis, the damage probability of different numbers of submunitions to a 4000 m runway falls near the damage probability curve of submunitions to a 3000 m runway. Therefore, the effect of the number of submunitions on the blockage probability of runways is essentially the distance among the preaiming points of submunitions. When we obtain the blockage probability curve of a certain type of submunition to an airport with different numbers, we can apply the curve to judge the blockage probability of this type of submunition to the airport with other lengths by conversion. Then, we can simplify the calculation of the blockage probability of the submunition to the airport runway.


#### 5.1.3. Influence of the Number of Subprojectiles on the Blockage Probability

Six submunitions and different numbers of subprojectiles are used to simulate the impact of the number of subprojectiles on the blockage probability. Other data are the same as above. The simulation results are shown in Table6.Table 6
Influence of the number of subprojectiles on the blockage probability.
Subprojectile number (pieces)406080100120140160180200Blockage probability (%)0.335.67.810.112.715.118.722.2Subprojectile number (pieces)210220240260280300320340360Blockage probability (%)23.725.428.329.830.931.331.531.832.2Figure15 depicts that when the number of subprojectiles is 20–40, the blockage probability increases slowly. Given that the density of subprojectiles in the dispersion range is small, only a few subprojectiles fall on the runway, and increasing the number of subprojectiles cannot produce evident changes in the blockade probability. When the number of subprojectiles is 80–240, the impact of the number of subprojectiles on the blockage probability tends to increase. When the number of subprojectiles is more than 240, the curve tends to be flat. The analysis shows that when a large number of subprojectiles fall on the runway, no takeoff and landing window exists in the coverage area of the submunition, and the window exists in the area not covered by the submunition. Thus, the number of bullets should be reasonably selected to save the economic cost and ensure the blockade efficiency.Figure 15
Influence of the number of subprojectiles on the blockage probability.


#### 5.1.4. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability

Six submunitions are used, the number of subprojectiles is 120, and different dispersion radii of subprojectiles are chosen to simulate. The simulation results are shown in Table7.Table 7
Influence of the dispersion radius of subprojectiles on the blockage probability.
Dispersion radius of subprojectiles (m)255075100125150175200225250275300Blockage probability (%)000.42.67.310.119.52733.53737.838.9Dispersion radius of subprojectiles (m)325350375400425450475500525550575600Blockage probability (%)39.638.237.335.132.729.827.926.423.621.818.816.7Figure15 depicts that when the dispersion radius of subprojectiles increases from 75 m to 250 m, the blocking probability increases rapidly. The probability of the subprojectiles blocking the runway is affected by the coverage area. With the increase in the dispersion radius of subprojectiles, the probability of the subprojectiles landing on the airport runway increases. The coverage area also increases, which makes the blockade probability increase rapidly. When the radius increases from 250 m to 350 m, the change in the blockage probability is relatively gentle. At this time, the dispersion radius of subprojectiles begins to overlap, and the blocking efficiency decreases. When the radius increases from 350 m to 600 m, the density of subprojectiles and the number of subprojectiles falling on the runway decrease, which reduces the blockage probability. As a result, the influence of the dispersion radius of subprojectiles on the blockage probability has an extreme value. When optimizing the blocking performance, we should consider the influence of the dispersion radius of subprojectiles, and the dispersion radius should not be increased blindly. A certain subprojectile density should be ensured.


#### 5.1.5. Influence of CEP and Crater Size of Subprojectiles on the Blockage Probability

We select 10 submunitions to study the impact of the CEP and the size of the crater produced by subprojectiles on the blockage probability. Other data are the same as above. The simulation data are shown in Table8.Table 8
Influence of the CEP and crater size of subprojectiles on the blockage probability.
Radius of crater/mCEP100 (%)90 (%)80 (%)70 (%)60 (%)50 (%)40 (%)r = 149.363.575.48994.299.1100r = 25464.778.789.396.399.2100r = 357.270.681.591.897.699.3100r = 46072.683.692.797.799.4100The dispersion of submunitions at the preview point will produce the largest impact area on the runway. The CEP accuracy affects the probability of the submunitions falling at the preview point. In accordance with the data, with the decrease in CEP, the blockage probability of submunitions increases remarkably. Increasing the crater size of subprojectiles will also increase the blockage probability of submunitions, but the effect is unobvious compared with the effect when reducing CEP. With the increase in accuracy, when the CEP of submunitions is reduced to a certain range, the change in the influence area of submunitions on the runway at the actual dispersion point is reduced. The change in the blockage probability will also be reduced because the dispersion range of submunitions is circular. Although improving the accuracy of submunitions is more effective than increasing the explosive capability of subprojectiles, with the development of today’s airport runway repair technology, the repair time of a large crater is much longer than that of a small crater. Subprojectiles still need a certain power to ensure that the damage can produce sufficient blocking time to meet the tactical requirements.


#### 5.1.6. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability under Different Numbers of Subprojectiles

Six submunitions are used, and the numbers of subprojectiles are 120, 210, and 300. The remaining data are the same as above to study the influence of the dispersion radius of subprojectiles on the blockage probability under different numbers of subprojectiles. The simulation results are shown in Figure16.Figure 16
Influence of the dispersion radius of subprojectiles on the blockage probability under different numbers of subprojectiles.In accordance with the simulation results, when the dispersion radius is 50 m and the number of subprojectiles is 300, the blockage probability is 0.3%. When the number of other subprojectiles is 0, increasing the number of subprojectiles can produce blocking at a lower dispersion radius. In the case of different numbers of subprojectiles, the gentle change range of the blockage probability of submunitions is different. Increasing the number of subprojectiles can increase the gentle change range, but the starting point of the gentle rise range is the same. The analysis shows that increasing the number of subprojectiles can ensure that a certain number of subprojectiles falls on the runway such that the subprojectile density can meet the requirements of maintaining or increasing the blockage probability. When the number of subprojectiles is extremely low, the blockage probability caused by increasing the dispersion radius of subprojectiles will decrease even if it does not reach the flat change range. Therefore, to ensure the stability of the blockage of the airport runway and reduce the impact caused by the fluctuation in the dispersion radius, we can select the appropriate number of submunitions to produce a certain gentle change range and choose a dispersion radius in the gentle change range to ensure the blockage of the airport runway.


#### 5.1.7. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability under Different CEP

To research the influence of the dispersion radius of subprojectiles on the blockage probability under different CEP, four submunitions are used, the number of subprojectiles is 150, and CEP is set to 1, 25, and 50 m. Other data are the same as above. The simulation results are shown in Figure17.Figure 17
Influence of the dispersion radius of subprojectiles on the blockage probability under different CEP.From the simulation results, changing the CEP of submunitions cannot change the trend of the effect exerted by the dispersion radius of subprojectiles on the blockage probability. The same gentle change range occurs under different CEP, and changing CEP only changes the blockage probability. This result shows that the influences of CEP and projectile radius are independent of each other. The analysis of Figures14, 18, and 19 indicates that the starting point of the gentle change interval in Figures 14 and 18 is 250 m; the starting point of the gentle rise interval in Figure 19 is 375 m, which is half of the distance among the preaiming points of shrapnel (d/2). Accordingly, the next research is to determine the reason through simulation.Figure 18
Influence of the dispersion radius of subprojectiles on the blockage probability.Figure 19
Influence of the CEP and crater size of subprojectiles on the blockage probability.


#### 5.1.8. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability under Different Preaiming Point Spacing of Submunitions

Four, six, and eight submunitions are selected for simulation. In accordance with the previous analysis, CEP does not affect the range of gentle change in the blockage probability. Given that the blockage probability of four submunitions is extremely small when CEP = 100 m, four submunitions are studied under the condition of CEP = 50 m, and six and eight submunitions are studied under the condition of CEP = 100 m. The simulation results are shown in Figure20.Figure 20
Influence of the dispersion radius of subprojectiles on the blockage probability under different preaiming point spacing of submunitions.The starting point of the range at which the blockage probability of submunitions changes gently is at half of the distance among the preaiming points of submunitions. The analysis shows that when the dispersion radius of submunitions is half of the distance among the preaiming points of submunitions, the dispersion ranges of subprojectiles at the preaiming points begin to overlap. As the dispersion radius continues to rise, the overlapping ranges continue to increase, the blockage efficiency begins to decline, and the changing trend of the blockage probability varies from rapid rise to slow rise. From the above analysis, we can conclude that the influence of the dispersion radius of subprojectiles on the starting point of the flat change interval is only related to the preaiming points of submunitions and has no connection with CEP and the number of subprojectiles.


### 5.2. Fitting Formula for Blocking Probability

In wars, judging whether airports can be blocked rapidly is sometime necessary. Therefore, a formula should be provided to calculate the probability of airport blockade accurately. When the dispersion radius of subprojectiles is greater than half of the distance among the preaiming points of submunitions (d/2), the blocking efficiency will decrease. Hence, the fitting range is selected in the area when the dispersion radius of subprojectiles is less than d/2. The runway length can be solved by converting it into the distance among the preaiming points of submunitions.In the above analysis, a coupling effect exists between the number of subprojectiles and the dispersion radius of subprojectiles. The number of subprojectiles can affect the flat range of the dispersion radius of subprojectiles in the blockage probability. Therefore, the following research is performed to verify whether the number of subprojectiles has an impact on the area before the starting point of the interval (d/2).The blockage probability when the dispersion radius of subprojectiles is less than 250 m (d/2) in Figure 16 is selected, and the processed data are presented in Table 9.Table 9
Processed data.
Dispersion radius of subprojectiles (m)Blockage probability (%)120 pieces/210 pieces210 pieces/300 pieces120 pieces/300 pieces750.2500.7620.1901000.5000.7540.3771250.6460.7060.4561500.3910.7870.3081750.5150.7580.3902000.5030.8510.4282250.5310.8410.4472500.4890.9000.440When the dispersion radius of subprojectiles is less than 250 m (d/2), a minimal difference in the ratio of the blockage probability of the dispersion radius of subprojectiles exists for each number of subprojectiles. The coupling effect between the number of subprojectiles and the dispersion radius of subprojectiles is small. In accordance with the previous analysis, the influence of various factors on the blocking probability can be expressed in the form of power function. Accordingly, the fitting formula has the following form:(11)P=A⋅Sα1⋅Cα2⋅Rα3⋅sα4⋅rα5⋅mα6+B.In the formula,P is the blockage probability, S is the number of submunitions, C is the CEP, R is the dispersion radius of subprojectiles, s is the number of bullets, r is the radius of crater, m is the length of the takeoff and landing window, and A, B, a1, a2, a3, a4, and a5 are the coefficients.The fitting results are as follows:(12)P=3.9175×10−6⋅S1.6⋅C−0.05⋅R1.0582⋅s0.3786⋅r0.05⋅m1.11−33.598.The correlation coefficient is 0.93. Given that the blockage probability of submunitions to airport runways is affected by the characteristics of runways and submunitions, various influencing factors exist, and the fitting formula has a certain error compared with simulation. Thus, in accordance with the needs of judging time and accuracy, the fitting formula or simulation times should be selected reasonably to achieve the corresponding demands.


## 5.1. Simulation Analysis of Submunition Blockage Probability

### 5.1.1. Influence of the Number of Submunitions on the Blockage Probability

To study the influence of the number of submunitions on the airport blockage probability, in accordance with the above data, the influence of the number of submunitions on the airport blockage probability is studied by changing the number of submunitions. The simulation results are shown in Figure12.Figure 12
Influence of the number of submunitions on the blockage probability.From the figure, increasing the number of submunitions can increase the blockage probability, which is approximately linear at the beginning. When the number of submunitions increases to a certain extent, the effect of increasing the number of submunitions on improving the blockage probability decreases. The turning point is at the place where the number of submunitions is 11. At this time, the dispersion radius of subprojectiles is 150 m, which is close to 136 m half of the distance among the preaiming points of the submunitions. The coverage area of subprojectiles has overlapped. This condition is speculated as the reason for the decline in blockage efficiency, which will be discussed later.


### 5.1.2. Influence of the Takeoff and Landing Window Length and Runway Length on the Blockage Probability

Owing to aircraft performance and other reasons, different aircraft need diverse takeoff and landing window lengths. Simulation analysis is carried out to study the influence of the takeoff and landing window length on the blockage probability. The runway length is 3000 m, and the number of submunitions is 6. Other data are the same as above. The simulation results are shown in Figure13.Figure 13
Influence of different takeoff and landing window lengths on the blockage probability.When the length of the takeoff and landing window increases, the blockage probability of the same number of submunitions increases. When the length of the takeoff and landing window is 250–400 m, the blocking probability increases more gently than that when the length is 400–800 m. Given that few subprojectiles landed on the runway and the probability of blockade is small, the difference caused by the change in takeoff and landing windows is minimal.Although the blockage probability of the same number of submunitions under different runway lengths is different, when the number of submunitions is changed such that the preaiming points of the submunitions are the same, the blockade probability is the same. To verify the analysis, the runway length is set to 4000 m, and other data are the same as above. The simulation results are shown in Table5.Table 5
Blockage probability simulation at 4000 m runway length.
Number of submunitions (pieces)89101112Blockage probability (%)10.121.223.435.938.8The simulation results are transformed into variables with the preaiming point spacing of submunitions. The data in Figure11 are processed and plotted in Figure 14.Figure 14
Influence of the preaiming point spacing of submunitions on the blockage probability in different runway lengths.When the preaiming points of submunitions are evenly distributed along the runway axis, the damage probability of different numbers of submunitions to a 4000 m runway falls near the damage probability curve of submunitions to a 3000 m runway. Therefore, the effect of the number of submunitions on the blockage probability of runways is essentially the distance among the preaiming points of submunitions. When we obtain the blockage probability curve of a certain type of submunition to an airport with different numbers, we can apply the curve to judge the blockage probability of this type of submunition to the airport with other lengths by conversion. Then, we can simplify the calculation of the blockage probability of the submunition to the airport runway.


### 5.1.3. Influence of the Number of Subprojectiles on the Blockage Probability

Six submunitions and different numbers of subprojectiles are used to simulate the impact of the number of subprojectiles on the blockage probability. Other data are the same as above. The simulation results are shown in Table6.Table 6
Influence of the number of subprojectiles on the blockage probability.
Subprojectile number (pieces)406080100120140160180200Blockage probability (%)0.335.67.810.112.715.118.722.2Subprojectile number (pieces)210220240260280300320340360Blockage probability (%)23.725.428.329.830.931.331.531.832.2Figure15 depicts that when the number of subprojectiles is 20–40, the blockage probability increases slowly. Given that the density of subprojectiles in the dispersion range is small, only a few subprojectiles fall on the runway, and increasing the number of subprojectiles cannot produce evident changes in the blockade probability. When the number of subprojectiles is 80–240, the impact of the number of subprojectiles on the blockage probability tends to increase. When the number of subprojectiles is more than 240, the curve tends to be flat. The analysis shows that when a large number of subprojectiles fall on the runway, no takeoff and landing window exists in the coverage area of the submunition, and the window exists in the area not covered by the submunition. Thus, the number of bullets should be reasonably selected to save the economic cost and ensure the blockade efficiency.Figure 15
Influence of the number of subprojectiles on the blockage probability.


### 5.1.4. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability

Six submunitions are used, the number of subprojectiles is 120, and different dispersion radii of subprojectiles are chosen to simulate. The simulation results are shown in Table7.Table 7
Influence of the dispersion radius of subprojectiles on the blockage probability.
Dispersion radius of subprojectiles (m)255075100125150175200225250275300Blockage probability (%)000.42.67.310.119.52733.53737.838.9Dispersion radius of subprojectiles (m)325350375400425450475500525550575600Blockage probability (%)39.638.237.335.132.729.827.926.423.621.818.816.7Figure15 depicts that when the dispersion radius of subprojectiles increases from 75 m to 250 m, the blocking probability increases rapidly. The probability of the subprojectiles blocking the runway is affected by the coverage area. With the increase in the dispersion radius of subprojectiles, the probability of the subprojectiles landing on the airport runway increases. The coverage area also increases, which makes the blockade probability increase rapidly. When the radius increases from 250 m to 350 m, the change in the blockage probability is relatively gentle. At this time, the dispersion radius of subprojectiles begins to overlap, and the blocking efficiency decreases. When the radius increases from 350 m to 600 m, the density of subprojectiles and the number of subprojectiles falling on the runway decrease, which reduces the blockage probability. As a result, the influence of the dispersion radius of subprojectiles on the blockage probability has an extreme value. When optimizing the blocking performance, we should consider the influence of the dispersion radius of subprojectiles, and the dispersion radius should not be increased blindly. A certain subprojectile density should be ensured.


### 5.1.5. Influence of CEP and Crater Size of Subprojectiles on the Blockage Probability

We select 10 submunitions to study the impact of the CEP and the size of the crater produced by subprojectiles on the blockage probability. Other data are the same as above. The simulation data are shown in Table8.Table 8
Influence of the CEP and crater size of subprojectiles on the blockage probability.
Radius of crater/mCEP100 (%)90 (%)80 (%)70 (%)60 (%)50 (%)40 (%)r = 149.363.575.48994.299.1100r = 25464.778.789.396.399.2100r = 357.270.681.591.897.699.3100r = 46072.683.692.797.799.4100The dispersion of submunitions at the preview point will produce the largest impact area on the runway. The CEP accuracy affects the probability of the submunitions falling at the preview point. In accordance with the data, with the decrease in CEP, the blockage probability of submunitions increases remarkably. Increasing the crater size of subprojectiles will also increase the blockage probability of submunitions, but the effect is unobvious compared with the effect when reducing CEP. With the increase in accuracy, when the CEP of submunitions is reduced to a certain range, the change in the influence area of submunitions on the runway at the actual dispersion point is reduced. The change in the blockage probability will also be reduced because the dispersion range of submunitions is circular. Although improving the accuracy of submunitions is more effective than increasing the explosive capability of subprojectiles, with the development of today’s airport runway repair technology, the repair time of a large crater is much longer than that of a small crater. Subprojectiles still need a certain power to ensure that the damage can produce sufficient blocking time to meet the tactical requirements.


### 5.1.6. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability under Different Numbers of Subprojectiles

Six submunitions are used, and the numbers of subprojectiles are 120, 210, and 300. The remaining data are the same as above to study the influence of the dispersion radius of subprojectiles on the blockage probability under different numbers of subprojectiles. The simulation results are shown in Figure16.Figure 16
Influence of the dispersion radius of subprojectiles on the blockage probability under different numbers of subprojectiles.In accordance with the simulation results, when the dispersion radius is 50 m and the number of subprojectiles is 300, the blockage probability is 0.3%. When the number of other subprojectiles is 0, increasing the number of subprojectiles can produce blocking at a lower dispersion radius. In the case of different numbers of subprojectiles, the gentle change range of the blockage probability of submunitions is different. Increasing the number of subprojectiles can increase the gentle change range, but the starting point of the gentle rise range is the same. The analysis shows that increasing the number of subprojectiles can ensure that a certain number of subprojectiles falls on the runway such that the subprojectile density can meet the requirements of maintaining or increasing the blockage probability. When the number of subprojectiles is extremely low, the blockage probability caused by increasing the dispersion radius of subprojectiles will decrease even if it does not reach the flat change range. Therefore, to ensure the stability of the blockage of the airport runway and reduce the impact caused by the fluctuation in the dispersion radius, we can select the appropriate number of submunitions to produce a certain gentle change range and choose a dispersion radius in the gentle change range to ensure the blockage of the airport runway.


### 5.1.7. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability under Different CEP

To research the influence of the dispersion radius of subprojectiles on the blockage probability under different CEP, four submunitions are used, the number of subprojectiles is 150, and CEP is set to 1, 25, and 50 m. Other data are the same as above. The simulation results are shown in Figure17.Figure 17
Influence of the dispersion radius of subprojectiles on the blockage probability under different CEP.From the simulation results, changing the CEP of submunitions cannot change the trend of the effect exerted by the dispersion radius of subprojectiles on the blockage probability. The same gentle change range occurs under different CEP, and changing CEP only changes the blockage probability. This result shows that the influences of CEP and projectile radius are independent of each other. The analysis of Figures14, 18, and 19 indicates that the starting point of the gentle change interval in Figures 14 and 18 is 250 m; the starting point of the gentle rise interval in Figure 19 is 375 m, which is half of the distance among the preaiming points of shrapnel (d/2). Accordingly, the next research is to determine the reason through simulation.Figure 18
Influence of the dispersion radius of subprojectiles on the blockage probability.Figure 19
Influence of the CEP and crater size of subprojectiles on the blockage probability.


### 5.1.8. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability under Different Preaiming Point Spacing of Submunitions

Four, six, and eight submunitions are selected for simulation. In accordance with the previous analysis, CEP does not affect the range of gentle change in the blockage probability. Given that the blockage probability of four submunitions is extremely small when CEP = 100 m, four submunitions are studied under the condition of CEP = 50 m, and six and eight submunitions are studied under the condition of CEP = 100 m. The simulation results are shown in Figure20.Figure 20
Influence of the dispersion radius of subprojectiles on the blockage probability under different preaiming point spacing of submunitions.The starting point of the range at which the blockage probability of submunitions changes gently is at half of the distance among the preaiming points of submunitions. The analysis shows that when the dispersion radius of submunitions is half of the distance among the preaiming points of submunitions, the dispersion ranges of subprojectiles at the preaiming points begin to overlap. As the dispersion radius continues to rise, the overlapping ranges continue to increase, the blockage efficiency begins to decline, and the changing trend of the blockage probability varies from rapid rise to slow rise. From the above analysis, we can conclude that the influence of the dispersion radius of subprojectiles on the starting point of the flat change interval is only related to the preaiming points of submunitions and has no connection with CEP and the number of subprojectiles.


## 5.1.1. Influence of the Number of Submunitions on the Blockage Probability

To study the influence of the number of submunitions on the airport blockage probability, in accordance with the above data, the influence of the number of submunitions on the airport blockage probability is studied by changing the number of submunitions. The simulation results are shown in Figure12.Figure 12
Influence of the number of submunitions on the blockage probability.From the figure, increasing the number of submunitions can increase the blockage probability, which is approximately linear at the beginning. When the number of submunitions increases to a certain extent, the effect of increasing the number of submunitions on improving the blockage probability decreases. The turning point is at the place where the number of submunitions is 11. At this time, the dispersion radius of subprojectiles is 150 m, which is close to 136 m half of the distance among the preaiming points of the submunitions. The coverage area of subprojectiles has overlapped. This condition is speculated as the reason for the decline in blockage efficiency, which will be discussed later.


## 5.1.2. Influence of the Takeoff and Landing Window Length and Runway Length on the Blockage Probability

Owing to aircraft performance and other reasons, different aircraft need diverse takeoff and landing window lengths. Simulation analysis is carried out to study the influence of the takeoff and landing window length on the blockage probability. The runway length is 3000 m, and the number of submunitions is 6. Other data are the same as above. The simulation results are shown in Figure13.Figure 13
Influence of different takeoff and landing window lengths on the blockage probability.When the length of the takeoff and landing window increases, the blockage probability of the same number of submunitions increases. When the length of the takeoff and landing window is 250–400 m, the blocking probability increases more gently than that when the length is 400–800 m. Given that few subprojectiles landed on the runway and the probability of blockade is small, the difference caused by the change in takeoff and landing windows is minimal.Although the blockage probability of the same number of submunitions under different runway lengths is different, when the number of submunitions is changed such that the preaiming points of the submunitions are the same, the blockade probability is the same. To verify the analysis, the runway length is set to 4000 m, and other data are the same as above. The simulation results are shown in Table5.Table 5
Blockage probability simulation at 4000 m runway length.
Number of submunitions (pieces)89101112Blockage probability (%)10.121.223.435.938.8The simulation results are transformed into variables with the preaiming point spacing of submunitions. The data in Figure11 are processed and plotted in Figure 14.Figure 14
Influence of the preaiming point spacing of submunitions on the blockage probability in different runway lengths.When the preaiming points of submunitions are evenly distributed along the runway axis, the damage probability of different numbers of submunitions to a 4000 m runway falls near the damage probability curve of submunitions to a 3000 m runway. Therefore, the effect of the number of submunitions on the blockage probability of runways is essentially the distance among the preaiming points of submunitions. When we obtain the blockage probability curve of a certain type of submunition to an airport with different numbers, we can apply the curve to judge the blockage probability of this type of submunition to the airport with other lengths by conversion. Then, we can simplify the calculation of the blockage probability of the submunition to the airport runway.


## 5.1.3. Influence of the Number of Subprojectiles on the Blockage Probability

Six submunitions and different numbers of subprojectiles are used to simulate the impact of the number of subprojectiles on the blockage probability. Other data are the same as above. The simulation results are shown in Table6.Table 6
Influence of the number of subprojectiles on the blockage probability.
Subprojectile number (pieces)406080100120140160180200Blockage probability (%)0.335.67.810.112.715.118.722.2Subprojectile number (pieces)210220240260280300320340360Blockage probability (%)23.725.428.329.830.931.331.531.832.2Figure15 depicts that when the number of subprojectiles is 20–40, the blockage probability increases slowly. Given that the density of subprojectiles in the dispersion range is small, only a few subprojectiles fall on the runway, and increasing the number of subprojectiles cannot produce evident changes in the blockade probability. When the number of subprojectiles is 80–240, the impact of the number of subprojectiles on the blockage probability tends to increase. When the number of subprojectiles is more than 240, the curve tends to be flat. The analysis shows that when a large number of subprojectiles fall on the runway, no takeoff and landing window exists in the coverage area of the submunition, and the window exists in the area not covered by the submunition. Thus, the number of bullets should be reasonably selected to save the economic cost and ensure the blockade efficiency.Figure 15
Influence of the number of subprojectiles on the blockage probability.


## 5.1.4. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability

Six submunitions are used, the number of subprojectiles is 120, and different dispersion radii of subprojectiles are chosen to simulate. The simulation results are shown in Table7.Table 7
Influence of the dispersion radius of subprojectiles on the blockage probability.
Dispersion radius of subprojectiles (m)255075100125150175200225250275300Blockage probability (%)000.42.67.310.119.52733.53737.838.9Dispersion radius of subprojectiles (m)325350375400425450475500525550575600Blockage probability (%)39.638.237.335.132.729.827.926.423.621.818.816.7Figure15 depicts that when the dispersion radius of subprojectiles increases from 75 m to 250 m, the blocking probability increases rapidly. The probability of the subprojectiles blocking the runway is affected by the coverage area. With the increase in the dispersion radius of subprojectiles, the probability of the subprojectiles landing on the airport runway increases. The coverage area also increases, which makes the blockade probability increase rapidly. When the radius increases from 250 m to 350 m, the change in the blockage probability is relatively gentle. At this time, the dispersion radius of subprojectiles begins to overlap, and the blocking efficiency decreases. When the radius increases from 350 m to 600 m, the density of subprojectiles and the number of subprojectiles falling on the runway decrease, which reduces the blockage probability. As a result, the influence of the dispersion radius of subprojectiles on the blockage probability has an extreme value. When optimizing the blocking performance, we should consider the influence of the dispersion radius of subprojectiles, and the dispersion radius should not be increased blindly. A certain subprojectile density should be ensured.


## 5.1.5. Influence of CEP and Crater Size of Subprojectiles on the Blockage Probability

We select 10 submunitions to study the impact of the CEP and the size of the crater produced by subprojectiles on the blockage probability. Other data are the same as above. The simulation data are shown in Table8.Table 8
Influence of the CEP and crater size of subprojectiles on the blockage probability.
Radius of crater/mCEP100 (%)90 (%)80 (%)70 (%)60 (%)50 (%)40 (%)r = 149.363.575.48994.299.1100r = 25464.778.789.396.399.2100r = 357.270.681.591.897.699.3100r = 46072.683.692.797.799.4100The dispersion of submunitions at the preview point will produce the largest impact area on the runway. The CEP accuracy affects the probability of the submunitions falling at the preview point. In accordance with the data, with the decrease in CEP, the blockage probability of submunitions increases remarkably. Increasing the crater size of subprojectiles will also increase the blockage probability of submunitions, but the effect is unobvious compared with the effect when reducing CEP. With the increase in accuracy, when the CEP of submunitions is reduced to a certain range, the change in the influence area of submunitions on the runway at the actual dispersion point is reduced. The change in the blockage probability will also be reduced because the dispersion range of submunitions is circular. Although improving the accuracy of submunitions is more effective than increasing the explosive capability of subprojectiles, with the development of today’s airport runway repair technology, the repair time of a large crater is much longer than that of a small crater. Subprojectiles still need a certain power to ensure that the damage can produce sufficient blocking time to meet the tactical requirements.


## 5.1.6. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability under Different Numbers of Subprojectiles

Six submunitions are used, and the numbers of subprojectiles are 120, 210, and 300. The remaining data are the same as above to study the influence of the dispersion radius of subprojectiles on the blockage probability under different numbers of subprojectiles. The simulation results are shown in Figure16.Figure 16
Influence of the dispersion radius of subprojectiles on the blockage probability under different numbers of subprojectiles.In accordance with the simulation results, when the dispersion radius is 50 m and the number of subprojectiles is 300, the blockage probability is 0.3%. When the number of other subprojectiles is 0, increasing the number of subprojectiles can produce blocking at a lower dispersion radius. In the case of different numbers of subprojectiles, the gentle change range of the blockage probability of submunitions is different. Increasing the number of subprojectiles can increase the gentle change range, but the starting point of the gentle rise range is the same. The analysis shows that increasing the number of subprojectiles can ensure that a certain number of subprojectiles falls on the runway such that the subprojectile density can meet the requirements of maintaining or increasing the blockage probability. When the number of subprojectiles is extremely low, the blockage probability caused by increasing the dispersion radius of subprojectiles will decrease even if it does not reach the flat change range. Therefore, to ensure the stability of the blockage of the airport runway and reduce the impact caused by the fluctuation in the dispersion radius, we can select the appropriate number of submunitions to produce a certain gentle change range and choose a dispersion radius in the gentle change range to ensure the blockage of the airport runway.


## 5.1.7. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability under Different CEP

To research the influence of the dispersion radius of subprojectiles on the blockage probability under different CEP, four submunitions are used, the number of subprojectiles is 150, and CEP is set to 1, 25, and 50 m. Other data are the same as above. The simulation results are shown in Figure17.Figure 17
Influence of the dispersion radius of subprojectiles on the blockage probability under different CEP.From the simulation results, changing the CEP of submunitions cannot change the trend of the effect exerted by the dispersion radius of subprojectiles on the blockage probability. The same gentle change range occurs under different CEP, and changing CEP only changes the blockage probability. This result shows that the influences of CEP and projectile radius are independent of each other. The analysis of Figures14, 18, and 19 indicates that the starting point of the gentle change interval in Figures 14 and 18 is 250 m; the starting point of the gentle rise interval in Figure 19 is 375 m, which is half of the distance among the preaiming points of shrapnel (d/2). Accordingly, the next research is to determine the reason through simulation.Figure 18
Influence of the dispersion radius of subprojectiles on the blockage probability.Figure 19
Influence of the CEP and crater size of subprojectiles on the blockage probability.


## 5.1.8. Influence of the Dispersion Radius of Subprojectiles on the Blockage Probability under Different Preaiming Point Spacing of Submunitions

Four, six, and eight submunitions are selected for simulation. In accordance with the previous analysis, CEP does not affect the range of gentle change in the blockage probability. Given that the blockage probability of four submunitions is extremely small when CEP = 100 m, four submunitions are studied under the condition of CEP = 50 m, and six and eight submunitions are studied under the condition of CEP = 100 m. The simulation results are shown in Figure20.Figure 20
Influence of the dispersion radius of subprojectiles on the blockage probability under different preaiming point spacing of submunitions.The starting point of the range at which the blockage probability of submunitions changes gently is at half of the distance among the preaiming points of submunitions. The analysis shows that when the dispersion radius of submunitions is half of the distance among the preaiming points of submunitions, the dispersion ranges of subprojectiles at the preaiming points begin to overlap. As the dispersion radius continues to rise, the overlapping ranges continue to increase, the blockage efficiency begins to decline, and the changing trend of the blockage probability varies from rapid rise to slow rise. From the above analysis, we can conclude that the influence of the dispersion radius of subprojectiles on the starting point of the flat change interval is only related to the preaiming points of submunitions and has no connection with CEP and the number of subprojectiles.


## 5.2. Fitting Formula for Blocking Probability

In wars, judging whether airports can be blocked rapidly is sometime necessary. Therefore, a formula should be provided to calculate the probability of airport blockade accurately. When the dispersion radius of subprojectiles is greater than half of the distance among the preaiming points of submunitions (d/2), the blocking efficiency will decrease. Hence, the fitting range is selected in the area when the dispersion radius of subprojectiles is less than d/2. The runway length can be solved by converting it into the distance among the preaiming points of submunitions.In the above analysis, a coupling effect exists between the number of subprojectiles and the dispersion radius of subprojectiles. The number of subprojectiles can affect the flat range of the dispersion radius of subprojectiles in the blockage probability. Therefore, the following research is performed to verify whether the number of subprojectiles has an impact on the area before the starting point of the interval (d/2).The blockage probability when the dispersion radius of subprojectiles is less than 250 m (d/2) in Figure 16 is selected, and the processed data are presented in Table 9.Table 9
Processed data.
Dispersion radius of subprojectiles (m)Blockage probability (%)120 pieces/210 pieces210 pieces/300 pieces120 pieces/300 pieces750.2500.7620.1901000.5000.7540.3771250.6460.7060.4561500.3910.7870.3081750.5150.7580.3902000.5030.8510.4282250.5310.8410.4472500.4890.9000.440When the dispersion radius of subprojectiles is less than 250 m (d/2), a minimal difference in the ratio of the blockage probability of the dispersion radius of subprojectiles exists for each number of subprojectiles. The coupling effect between the number of subprojectiles and the dispersion radius of subprojectiles is small. In accordance with the previous analysis, the influence of various factors on the blocking probability can be expressed in the form of power function. Accordingly, the fitting formula has the following form:(11)P=A⋅Sα1⋅Cα2⋅Rα3⋅sα4⋅rα5⋅mα6+B.In the formula,P is the blockage probability, S is the number of submunitions, C is the CEP, R is the dispersion radius of subprojectiles, s is the number of bullets, r is the radius of crater, m is the length of the takeoff and landing window, and A, B, a1, a2, a3, a4, and a5 are the coefficients.The fitting results are as follows:(12)P=3.9175×10−6⋅S1.6⋅C−0.05⋅R1.0582⋅s0.3786⋅r0.05⋅m1.11−33.598.The correlation coefficient is 0.93. Given that the blockage probability of submunitions to airport runways is affected by the characteristics of runways and submunitions, various influencing factors exist, and the fitting formula has a certain error compared with simulation. Thus, in accordance with the needs of judging time and accuracy, the fitting formula or simulation times should be selected reasonably to achieve the corresponding demands.


## 6. Conclusion

The proposed method has the advantages of fast judgment speed and accurate results. This article considers the existence of an oblique takeoff and landing window. Compared with the window scanning method, the proposed method solves the problems of error caused by excessively large scanning step and extremely long calculation time caused by considerably small scanning step. When studying the attack capability of submunitions to airports, we should consider various factors, such as the dispersion radius of subprojectiles, CEP, the number of submunitions, the crater area, and the economic cost. In accordance with the simulation analysis, the following conclusions are drawn:(1)
The starting point of the smooth range of the radius of subprojectile scattering coverage circle in the blockage probability is 1/2 of the distance among the preaiming points of submunitions. It is not influenced by the CEP of submunitions and the number of subprojectiles. When the distance is less than this value, the blocking efficiency will decrease.(2)
Through simulation, the inflection point of the blocking efficiency of submunitions is found. The research results provide a reference for the design of submunition and the optimization of strike scheme. When the number of subprojectiles is excessively low, the blockage probability caused by increasing the dispersion radius of subprojectiles will decrease even if it does not reach the inflection point. The dispersion radius and the number of submunitions should be considered comprehensively. The subprojectiles should be designed with reasonable dispersion density.(3)
This article improves the judgment formula of takeoff and landing window. The proposed method has a wide range of applicability and feasibility. It can judge the impact of ordinary missiles, submunitions, and other warheads on airport runway damage. A fitting formula for the probability of submunition blocking the airport is given. The proposed method and the formula can provide help for the study of the best depth of explosion, the best number of missiles, the dispersion radius of subprojectiles, CEP, and the best size of bullet craters.[17–24]



---

*Source: 5631800-2022-04-29.xml*