A method for constructing a three-dimensional static-dynamic blast power field of a preformed fragment warhead

Abstract

The application provides a prefabricated fragment warhead three-dimensional static and dynamic blast power field construction method, comprising the following steps: step one, prefabricated fragment warhead structure modeling; step two, fragment sampling; step three, fragment static blast initial velocity calculation and correction; step four, fragment static blast scattering direction angle calculation; step five, fragment dynamic blast initial velocity and dynamic blast scattering direction angle calculation; step six, prefabricated fragment warhead three-dimensional static and dynamic blast power field construction: steps two to five are performed for each fragment until the dynamic blast initial velocity and dynamic blast scattering direction angle of all fragments are calculated, and then the prefabricated fragment warhead three-dimensional static and dynamic blast power field under the specified terminal velocity is obtained. The application comprehensively considers the influence of the difference between fragments, the leakage of detonation products, the end effect, randomness and the like. It is verified that the application can accurately construct the three-dimensional fragment field of a complex prefabricated fragment warhead containing multiple types and large magnitude fragments, thereby providing support for warhead power and damage efficiency evaluation.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of explosives and damage, and relates to a prefabricated fragment warhead, in particular to a method for constructing a three-dimensional static and dynamic blast power field of a prefabricated fragment warhead. BACKGROUND

[0002] The blast-kill warhead is widely used in the battlefield, and the accurate calculation of the blast field is of great significance for the damage efficiency evaluation and actual operation of the ammunition. The fragment is one of the core damage elements of the blast-kill warhead, and the prefabricated fragment is a typical representative.

[0003] The engineering calculation method of the initial speed of the prefabricated fragment currently has the following problems: ①Researchers often take cylindrical or flat charging as the research object, pay attention to the relationship between the shape and mass of the prefabricated fragment and the size of the fragment ejection initial speed, but ignore the influence of the fragment arrangement and the shape of the warhead shell on the ejection of the prefabricated fragment; ②The randomness of the fragment driving and releasing process is not considered, so that the initial speed and ejection direction of the fragments at the same position along the axis of the bomb are the same, which is inconsistent with the actual situation; ③The engineering algorithm can only give one-dimensional distribution of the initial speed of the fragment along the axis of the bomb, but the calculation of the damage efficiency of the ammunition needs to consider the three-dimensional interaction between the fragment field and the target, and the conversion method from one dimension to three dimensions is not perfect; ④For the warhead with complex structure, multiple types of fragments and large number of fragments (more than 10w), it is difficult to meet the data requirements of traditional engineering modeling by relying on experiments and simulation.

[0004] Therefore, it is urgent to provide a three-dimensional prefabricated fragment field calculation method considering the structure of the bomb and the randomness of the initial driving of the fragment, to support the calculation of the damage efficiency of the blast-kill ammunition. SUMMARY

[0005] In view of the deficiencies in the prior art, the purpose of the present application is to provide a method for constructing a three-dimensional static and dynamic blast power field of a prefabricated fragment warhead, to solve the technical problem that the calculation accuracy of the damage efficiency of the blast-kill ammunition in the prior art needs to be further improved.

[0006] In order to solve the above technical problems, the present application adopts the following technical solutions:

[0007] A method for constructing a three-dimensional static and dynamic blast power field of a prefabricated fragment warhead, which is performed according to the following steps:

[0008] Step 1, modeling of the prefabricated fragment warhead structure:

[0009] Modeling the prefabricated fragment warhead structure according to the prefabricated fragment warhead structure to obtain the shell contour curve.

[0010] Step 2, fragment sampling:

[0011] For different types of fragments, the single fragment is taken as the object, and the specific position of the fragment in the prefabricated fragment warhead is determined by sampling.

[0012] Step three, calculation and correction of the static explosion initial velocity of the fragment:

[0013] Step 301, calculation of the average initial velocity v0 of the fragment.

[0014] Step 302, based on v0, the end effect is corrected to obtain v 0x1 .

[0015] Step 303, based on v 0x1 , the leakage problem of the detonation product is corrected to obtain v 0x2 .

[0016] Step 304, based on v 0x2 , the speed deviation caused by the mass difference of the fragment is corrected to obtain v 0x .

[0017] Wherein:

[0018] v0 is the average initial velocity of the fragment, in m / s;

[0019] v 0x1 is the speed of the fragment at different positions from the explosion end surface after correction of the end effect, in m / s;

[0020] v 0x2 is the static explosion initial velocity of the fragment after correction of the leakage problem, in m / s;

[0021] v 0x is the static explosion initial velocity of the fragment after correction, in m / s.

[0022] Step four, calculation of the scattering direction angle of the static explosion fragment:

[0023] Step 401, according to the shell shape curve obtained in step one and the specific position of the fragment in the prefabricated fragment warhead obtained in step two, the outer normal direction of the shell at the corresponding position is calculated, and the included angle between the outer normal of the shell surface of the prefabricated fragment warhead and the positive direction of the axis of the shell is represented as .

[0024] Step 402, based on obtained in step 401, the direction deflection angle θ s of the fragment is calculated according to the Shapiro formula.

[0025] Step 403, calculation of the scattering direction angle of the fragment under the treatment of the Shapiro formula

[0026] In the formula,

[0027] is the direction angle of the fragment dispersion under the Shapiro formula processing, and the unit is rad;

[0028] is the direction angle of the fragment dispersion under the Shapiro formula processing, and the unit is rad;

[0029] is the direction angle of the fragment dispersion under the Shapiro formula processing, and the unit is rad; s is the direction angle of the fragment dispersion under the Shapiro formula processing, and the unit is rad.

[0030] Step five, the calculation of the initial velocity and the direction angle of the dynamic explosion of the fragment:

[0031] On the basis of step four, based on the vector superposition principle, the superposition result of the terminal speed vector of the ammunition and the speed vector of the fragment is calculated, and the initial velocity and the direction angle of the dynamic explosion of the fragment are obtained.

[0032] Step six, the construction of the three-dimensional static and dynamic explosion power field of the prefabricated fragment warhead:

[0033] Steps two to five are performed for each fragment until the initial velocity and the direction angle of the dynamic explosion of all fragments are calculated, and then the three-dimensional static and dynamic explosion power field of the prefabricated fragment warhead under the specified terminal speed is obtained.

[0034] Compared with the prior art, the present application has the following technical effects:

[0035] (I) The present application comprehensively considers the influence of the fragment arrangement form, the mass difference of the fragment, the structure of the warhead shell and other factors on the fragment field, so that the accuracy of the calculation result is higher.

[0036] (II) The present application considers the randomness of the driving and releasing process of the fragment, and can better reflect the real scattering situation of the fragment field.

[0037] (III) The present application calculates based on a single fragment, and finally constructs a three-dimensional static and dynamic explosion power field of the prefabricated fragment warhead, which can meet the calculation input demand of various damage efficiency algorithms.

[0038] (IV) The present application comprehensively considers the influence of the difference between the types of fragments, the leakage of detonation products, the end effect and randomness. It has been verified that the present application can accurately construct a three-dimensional fragment field of a complex prefabricated fragment warhead containing multiple types and large magnitudes of fragments, thereby providing support for the power and damage efficiency evaluation of the warhead. BRIEF DESCRIPTION OF DRAWINGS

[0039] Figure 1 is a prefabricated fragment warhead structure schematic diagram of the embodiment of the present application.

[0040] Figure 2 is a warhead coordinate system construction schematic diagram of the embodiment of the present application.

[0041] Figure 3 is a discrete collection result diagram of the warhead shell shape curve of the embodiment of the present application.

[0042] Figure 4 is a modeling schematic diagram of the fragment distribution area of the embodiment of the present application.

[0043] Figure 5 is a three-dimensional space distribution (left) and ground distribution (right) diagram of the fragment field (static explosion-5m explosion height-90° drop angle) of the embodiment of the present application.

[0044] Figure 6 is a three-dimensional space distribution (left) and ground distribution (right) diagram of the fragment field (static explosion-5m explosion height-75° drop angle) of the embodiment of the present application.

[0045] Figure 7 is a three-dimensional space distribution (left) and ground distribution (right) diagram of the fragment field (340m / s drop speed-5m explosion height-90° drop angle) of the embodiment of the present application.

[0046] Figure 8 is a three-dimensional space distribution (left) and ground distribution (right) diagram of the fragment field (340m / s drop speed-5m explosion height-75° drop angle) of the embodiment of the present application.

[0047] Figure 9 is a diagram of the relationship between the initial speed of each fragment and the scattering direction angle of the fragment field under the static explosion condition of the embodiment of the present application.

[0048] Figure 10 is a diagram of the relationship between the initial speed of each fragment and the scattering direction angle of the fragment field under the 340m / s drop speed condition of the embodiment of the present application.

[0049] The meanings of the various reference numerals in the drawings are as follows: 1 - charge, 2 - fragment layer, 3 - outer shell.

[0050] The specific content of the present application is further explained and described in detail below in combination with embodiments. DETAILED DESCRIPTION

[0051] It should be noted that the components and devices in the present application, if not specifically stated, all adopt components and devices known in the prior art.

[0052] The specific embodiments of the present application are given below, and it should be noted that the present application is not limited to the following specific embodiments, and any equivalent transformation made on the basis of the technical solutions of the present application falls within the protection scope of the present application.

[0053] Embodiment:

[0054] The present embodiment gives a method for constructing a three-dimensional static and dynamic explosion power field of a prefabricated fragment warhead, which is performed according to the following steps:

[0055] Step one, prefabricated fragment warhead structure modeling:

[0056] According to the prefabricated fragment warhead structure, the shell shape curve of the prefabricated fragment warhead structure is obtained.

[0057] In this embodiment, the prefabricated fragment warhead structure is as shown in Figure 1 The total weight of the warhead is 90 kg, the total length is 623 mm, and the maximum outer diameter is 213 mm; the charge is TNT, the charge length is 500 mm, and the charge amount is 20 kg. 6000 pieces of 7mm prefabricated fragments and 21000 pieces of 5mm prefabricated fragments are loaded.

[0058] Step one includes the following steps:

[0059] Step 101, establish the coordinate system of the prefabricated fragment warhead:

[0060] As shown in Figure 2 , a cylindrical coordinate system is established with the center of the bullet head as the origin, and the prefabricated fragment warhead central axis, the prefabricated fragment warhead radial direction and the prefabricated fragment warhead circumferential azimuth angle as the reference. In the cylindrical coordinate system, the prefabricated fragment warhead central axis is the x-axis, and the prefabricated fragment warhead radial direction is the y-axis.

[0061] In this embodiment, the red dot in Figure 2 is the detonation point.

[0062] Step 102, obtain the prefabricated fragment warhead structure parameters:

[0063] In the cylindrical coordinate system established in step 101, the data is read by taking points to obtain the length of the prefabricated fragment warhead and the charge length, the average charge diameter, the fragment layer thickness, the shell thickness, the x-coordinate of the detonation point and the shell shape curve.

[0064] In this embodiment, the length of the prefabricated fragment warhead is 623mm, the charge length is 500mm, the average charge diameter is 146mm, the fragment layer thickness is 22.7mm, the shell thickness is 10mm, and the x-coordinate of the detonation point is 61.37mm.

[0065] Further in step 102, points are taken on the outer surface of the prefabricated fragment warhead shell to realize the discretization and digitization of the shell shape curve; points are taken densely at places where the shape changes obviously, and a small amount of points are taken at places where the shape changes slowly, as shown in Figure 3 .

[0066] Step 103, modeling of the fragment distribution area:

[0067] According to the arrangement of the fragments, the fragments are partitioned along the bullet axis to obtain the number of fragments in each region.

[0068] In this specific embodiment, such as Figure 4 As shown, the 5mm fragments are divided into one region with an x-coordinate interval of (221, 514)mm, and the number of fragments in the interval is 21,000. The 7mm fragments are divided into two regions with x-coordinate intervals of (63, 221)mm and (514, 560)mm, respectively.

[0069] The ratio of interval lengths is approximated as the ratio of the number of fragments in each region, resulting in 4647 and 1353 fragments in the two regions, respectively. Alternatively, the number of fragments can be converted in other ways or directly assigned according to the detailed design parameters, depending on the actual situation.

[0070] Step 2, fragment sampling:

[0071] For different types of fragments, the specific location of each fragment within the pre-fragmented warhead is determined by sampling, using a single fragment as the target.

[0072] Step two includes the following steps:

[0073] Step 201: For fragments of the same type, random sampling is performed between 0° and 360° in the circumferential direction to determine the circumferential scattering direction angle of the fragments. The circumferential scattering direction angle is the azimuth angle of the fragments in the cylindrical coordinate system.

[0074] Step 202: Based on the total number of fragments of the same type and the number of fragments in each region, randomly sample to determine the region to which the fragment belongs.

[0075] Specifically, in this embodiment, for 7mm fragments, with a total of 6000 fragments, sampling is performed between 1 and 6000. If the result is between 1 and 4647, the fragment belongs to the first region; otherwise, it belongs to the second region. For 5mm fragments, since there is only one region, no sampling is required, or only formal sampling is needed.

[0076] Step 203: Based on the coordinate range of the region, randomly sample and determine the coordinates of the fragment along the projectile axis.

[0077] In this specific embodiment, if a 7mm fragment belongs to the first region, and the x-coordinate interval of this region is (221, 514)mm, then random sampling is performed between 221 and 514 to obtain the coordinates of the fragment.

[0078] Step 204: Based on the number of layers of fragments in the area, randomly sample to determine which layer the fragments are arranged in.

[0079] In this embodiment, a certain 7mm fragment belongs to the first region, the fragments in this region are arranged in three layers, and then random sampling can be performed between 0 and 3; if the result is [0, 1), the fragment comes from the first layer, if the result is [1, 2), the fragment comes from the first layer, and if the result is [2, 3], the fragment comes from the third layer.

[0080] Step three, calculation and correction of the initial velocity of the fragment static explosion:

[0081] Step 301, calculate the average initial velocity v0 of the fragment; the calculation method in step 301 is as follows:

[0082]

[0083] In the formula:

[0084] v0 is the average initial velocity of the fragment, and the unit is m / s;

[0085] is the Gurney constant;

[0086] D e is the detonation velocity, and the unit is m / s, for TNT, take 6900 m / s;

[0087] β is the load coefficient, that is, the ratio of the mass of the charge per unit length to the mass of the shell per unit length (here, the mass of the shell is approximated by the total mass of the prefabricated fragment warhead minus the mass of the charge) ;

[0088] M e is the mass of the charge, and the unit is kg;

[0089] L e is the length of the charge, and the unit is m;

[0090] M W is the total mass of the prefabricated fragment warhead, and the unit is kg;

[0091] L W is the length of the prefabricated fragment warhead, and the unit is m.

[0092] Step 302, on the basis of v0, the end effect is corrected to obtain v 0x1 ; the correction method of step 302 is as follows:

[0093]

[0094] In the formula:

[0095] v 0x1 is the velocity of the fragment at different positions from the detonation end face after correction of the end effect, and the unit is m / s;

[0096] v0 is the average initial velocity of the fragment, and the unit is m / s;

[0097] x is the distance from the fragment to the end surface of initiation, in m; obtained according to the coordinates of the fragment and the initiation point;

[0098] d is the average charge diameter, in m;

[0099] L e is the charge length, in m;

[0100] A, B, C, and D are all constant parameters.

[0101] In this embodiment, for the case of one-end initiation, the constant parameters A, B, C, and D are respectively taken as 0.361, -2.3617, 0.192, and -3.03.

[0102] Step 303, on the basis of v 0x1 , the leakage problem of the detonation product is corrected to obtain v 0x2 ; the correction method of step 303 is as follows:

[0103] v 0x2 = v 0x1 ·modify1

[0104] In the formula:

[0105] v 0x2 is the corrected initial velocity of the fragment, in m / s;

[0106] v 0x1 is the velocity of the fragment at different positions from the end surface of initiation after correction of the end effect, in m / s;

[0107] modify1 is a leakage correction coefficient.

[0108] In this embodiment, for the leakage correction coefficient modify1, the inner layer fragment takes 0.8, the next inner layer takes 0.85, and the other layers take 0.9.

[0109] Step 304, on the basis of v 0x2 , the speed deviation caused by the mass difference of the fragments is corrected to obtain v 0x ; the correction method of step 304 is as follows:

[0110] v 0x = v 0x2 ·modify2

[0111] In the formula:

[0112] v 0x is the corrected initial velocity of the fragment, in m / s;

[0113] v0x2 The corrected fragment static explosion initial velocity is m / s.

[0114] modify2 is a mass correction coefficient.

[0115] In the embodiment, the distribution of loading energy is not uniform during the process of loading the fragments by the detonation products, and the constraint of the materials such as the shell and the adhesive makes the velocity difference of the fragments not too large, so the simple energy distribution method cannot be used to calculate the difference between the initial velocities of the fragments.

[0116] In the embodiment, the research results show that the mass correction coefficient modify2 can be 0.95 for large mass fragments, 1.05 for small mass fragments, and 1 for medium mass fragments, and ideal correction results can be obtained. In the embodiment, the correction coefficient of the 7mm fragment is 0.95, and the correction coefficient of the 5mm fragment is 1.05.

[0117] In step 305, the normal sampling is used to process the randomness of the initial velocity of the fragment.

[0118] In the embodiment, the calculated v 0x is used as the mean value, and 60.79m / s is used as the standard deviation (satisfying the probability of 90% that the initial velocity of the fragment falls within the range of v 0x ±100m / s), a normal distribution is constructed, the static explosion initial velocity of the fragment is sampled, and the randomness of the initial velocity of the fragment is superimposed. In actual implementation, the researchers can also adjust the parameters of the normal distribution according to the specific circumstances.

[0119] Step four, calculation of the direction angle of the static explosion scattering of the fragment:

[0120] In step 401, the outer normal direction of the shell at the corresponding position is calculated according to the shell shape curve obtained in step one and the specific position of the fragment in the preformed fragment warhead obtained in step two, and the included angle between the outer normal of the shell surface of the preformed fragment warhead and the positive direction of the axis of the shell is denoted as

[0121] In step 402, the direction deflection angle θ of the fragment is calculated according to the Shapiro formula on the basis of s The calculation method in step 402 is as follows:

[0122]

[0123] In the formula, θ

[0124] θ s is the direction deflection angle of the fragment, and the unit is rad;

[0125] v 0x ​V0 is the static explosion initial velocity of the fragment, unit: m / s;

[0126] φ1 is the angle between the outer normal of the preformed fragment warhead shell surface and the positive direction of the shell axis, unit: rad;

[0127] φ2 is the angle between the propagation direction of the detonation wave and the positive direction of the shell axis when the detonation wave reaches a certain position, unit: rad; In this case, the angle between the line connecting the initiation point and the fragment position and the shell axis is taken;

[0128] x is the distance of the fragment from the initiation end surface, unit: m;

[0129] d is the average charge diameter, unit: m / s;

[0130] D e is the detonation velocity, unit: m / s.

[0131] Step 403, calculate the fragment scattering direction angle under Shapiro formula processing The calculation method is as follows:

[0132]

[0133] In this embodiment, the direction of the shell head is taken as the positive direction, and “-” is taken when the fragment position is in front of the initiation point, and “+” is taken when the fragment position is behind the initiation point. In this embodiment, the fragments are all located behind the initiation point, so all take “+”.

[0134] Step 404, use normal sampling to process the randomness of the fragment scattering direction.

[0135] In this embodiment, a normal distribution with mean μ and standard deviation σ is constructed, and the static explosion scattering direction angle of the fragment is obtained by random sampling The calculation method of μ and σ is as follows:

[0136]

[0137] σ=kθ s

[0138] In the formula:

[0139] k is the scattering control coefficient, and in this embodiment, k is taken as 0.304, which satisfies the probability that the fragment scattering direction angle falls between and is about 90%.

[0140] Step five, calculate the dynamic explosion initial velocity and dynamic explosion scattering direction angle of the fragment:

[0141] On the basis of step four, based on the vector superposition principle, the superposition result of the shell terminal velocity vector and the fragment velocity vector is calculated to obtain the dynamic explosion initial velocity and the dynamic explosion scattering direction angle of the fragment.

[0142] If only the fragment static blast field is studied, this step can be skipped.

[0143] Step six, three-dimensional static-dynamic blast field construction of prefabricated fragment warhead:

[0144] Steps two to five are performed for each fragment until the dynamic blast initial velocity and dynamic blast scattering direction angle of all fragments are calculated, and the three-dimensional static-dynamic blast field of the prefabricated fragment warhead under the specified terminal velocity is obtained.

[0145] In this embodiment, the three-dimensional space distribution and ground distribution of the fragment field under the conditions of static blast-5m burst height-90° fall angle, static blast-5m burst height-75° fall angle, 340m / s fall velocity-5m burst height-90° fall angle, and 340m / s fall velocity-5m burst height-75° fall angle are calculated, and are shown in Figure 5 、 Figure 6 、 Figure 7 、 Figure 8 respectively. The corresponding relationship between the initial velocity and scattering direction angle of the fragments under static blast and 340m / s fall velocity is also given, as shown in Figure 9 、 Figure 10 . The dark data points in the figure are 5mm fragments, and the light data points are 7mm fragments.

[0146] This embodiment aims to explain the invention in detail and does not involve actual prefabricated fragment warheads. The invention has been verified by experiments and simulation data, and the calculation deviation can be controlled within 20%, which confirms that the method can be effectively applied to the accurate construction of the fragment field of prefabricated fragment warheads and the damage effectiveness evaluation of kill-blast ammunition.

Claims

1. A method for constructing a three-dimensional static-dynamic blast power field of a preformed fragment warhead, characterized in that, The method is performed according to the following steps: Step one, prefabricated fragment warhead structure modeling: According to the prefabricated fragment warhead structure, the prefabricated fragment warhead structure is modeled, and the shell shape curve is obtained; Step two, fragment sampling: For different types of fragments, taking a single fragment as the object, the specific position of the fragment in the prefabricated fragment warhead is determined by sampling; Step two includes the following steps: Step 201, for fragments of the same type, randomly sample between 0-360° in the circumferential direction to determine the circumferential scattering direction angle of the fragment; Step 202, according to the total number of fragments of the same type and the number of fragments in each region, randomly sample to determine the region to which the fragment belongs; Step 203, according to the coordinate range of the region, randomly sample to determine the coordinate of the fragment along the axis of the bomb; Step 204, according to the number of layers of the fragment in the region, randomly sample to determine the layer in which the fragment is arranged; Step three, fragment static explosion initial velocity calculation and correction: Step 301, calculate the average initial velocity of the fragments ; Step 302, on the basis of the end effect is corrected to obtain ; Step 303, on the basis of , the problem of detonation product leakage is corrected to obtain ; The correction method of step 303 is as follows: In the formula: Fragment static explosion initial velocity after correction for leakage problem, in m / s; is a leakage correction coefficient; Step 304, on the basis of the speed deviation caused by the mass difference of the fragments is corrected to obtain ; The correction method of step 304 is as follows: In the formula: To amend the completed fragment static initial velocity, unit: m / s; is the mass correction factor; Step four, fragment static explosion scattering direction angle calculation: Step 401, according to the shell shape curve obtained in step one and the specific position of the fragment in the preformed fragment warhead obtained in step two, the outer normal direction of the shell at the corresponding position is calculated, so that the angle between the outer normal of the preformed fragment warhead shell surface and the positive direction of the projectile axis is obtained is performed Step 402, on the basis of the result of step 401, the deflection angle of the fragment direction is calculated according to Shapiro formula ; Step 403, the deflection angle of the fragment direction is calculated according to the formula of the deflection angle of the fragment direction in the case of the impact of the missile with the target ; Step 403, calculate the fragment dispersion direction angle under the Shapiro formula processing , ; In the formula: is the deflection angle of the fragment, in radians; Step five, fragment dynamic explosion initial velocity and dynamic explosion scattering direction angle calculation: On the basis of step four, based on the principle of vector superposition, the superposition result of the terminal velocity vector of the ammunition and the velocity vector of the fragment is calculated to obtain the dynamic explosion initial velocity and the dynamic explosion scattering direction angle of the fragment; Step six, prefabricated fragment warhead three-dimensional static and dynamic explosion power field construction: Steps two to five are performed for each fragment until the dynamic explosion initial velocity and the dynamic explosion scattering direction angle of all fragments are calculated, and then the prefabricated fragment warhead three-dimensional static and dynamic explosion power field under the specified terminal velocity is obtained.

2. The method of constructing the three-dimensional static-dynamic blast power field of the preformed fragment warhead according to claim 1, characterized in that, Step one includes the following steps: Step 101, establish the prefabricated fragment warhead coordinate system: A cylindrical coordinate system is established with the center of the bomb head as the origin, the central axis of the prefabricated fragment warhead, the radial direction of the prefabricated fragment warhead, and the azimuth angle of the prefabricated fragment warhead as the reference. In the cylindrical coordinate system, the central axis of the prefabricated fragment warhead is the x-axis, and the radial direction of the prefabricated fragment warhead is the y-axis; Step 102, obtain the prefabricated fragment warhead structure parameters: The data is read in the column coordinate system established in step 101, and the length of the preformed fragment warhead and the charge length, average charge diameter, fragment layer thickness, shell thickness, and initiation point are obtained x coordinates and shell shape curve; Step 103, fragment distribution region modeling: Divide the fragments along the axis to obtain the number of fragments in each region.

3. The method of claim 1, wherein the method is characterized by: The calculation method in step 301 is as follows: In the formula: Vavg = 800 m / s for the average initial velocity of the fragments; Gurney's constant; D is the detonation velocity, in m / s; is the load factor.

4. The method of claim 3, wherein the method is characterized by: The correction method of step 302 is as follows: In the formula: The end effect is corrected for the fragment velocity at different positions from the initiation end face, in units of m / s; x D is the distance from the end surface of the initiation to the fragment, in meters; d D is the average charge diameter in m; L is the charge length in m; 、 、 、 are constant parameters.

5. The method of constructing the three-dimensional static-dynamic blast power field of the preformed fragment warhead according to claim 1, characterized in that, The calculation method in step 402 is as follows: In the formula: is the deflection angle of the fragment, in radians; To amend the completed fragment static initial velocity, unit: m / s; is the angle between the direction of the detonation wave propagation and the positive direction of the axis of the projectile at the moment when the detonation wave reaches the point of interest, in radians. D is the distance from the end surface of the initiation to the fragment, in meters; D is the average charge diameter in m / s; D is the detonation velocity, in m / s.

Images