Disclosure of Invention
The invention aims to provide a method for constructing an anti-ship missile warhead explosion scaling model, which is characterized in that based on the size of an equivalent naked explosive of a warhead, the ultimate expansion radius of the warhead is obtained according to energy conservation calculation, and Poisson ratio parameters of materials are obtained according to material experiments at different temperatures, so that the material selection is determined; and performing equivalent scaling calculation according to the parameters, and constructing an anti-ship missile warhead explosion scaling model according to the equivalent scaling calculation result. The invention can simplify the missile model in the simulation, simultaneously keep the original characteristics, weaken the influence of some unnecessary factors on the simulation prediction, and ensure that the simulation prediction is simpler and more convenient and the result is more accurate. The invention can save the experiment cost, improve the personnel safety, and provide a convenient, rapid and reliable reference basis for the design and optimization of the warhead or the protection structure of the related weapon and the like.
The purpose of the invention is realized by the following technical scheme.
The invention discloses a method for constructing an anti-ship missile warhead explosion scaling model, which comprises the following steps of:
after the explosive explosion of the warhead, the internal energy and the kinetic energy of detonation products and the kinetic energy of high-speed fragments are finally converted into quasi-static air pressure to act on a structure, and the size of the equivalent bare explosive is calculated by measuring the quasi-static pressure;
step two, neglecting energy loss caused by the streaming effect of the explosive products, and calculating to obtain the ultimate expansion radius of the warhead through the energy conservation principle;
step three, determining the radius-thickness ratio of the shell according to the failure mode of the warhead;
selecting a proper material model according to material fracture experiments under different temperatures, different strain rates and various stresses;
step five, determining materials according to the equivalent bare charge size obtained in the step one, the limit expansion radius obtained in the step two, the diameter-thickness ratio of the shell obtained in the step three and the material determined in the step four, and according to the characteristic parameters of the explosive load in the cabin, such as the overpressure peak value delta p, the impulse i and the quasi-static pressure p of the shock waveqsAnd the explosive characteristics, air characteristics, structure, material and position parameters are combined to complete the establishment of the warship-resisting missile warhead explosion scaling model through impulse similarity law, proportion law and dimensional analysis.
Further comprises the following steps: and (4) predicting the mechanical property of the explosion of the anti-ship missile warhead explosion scaling model according to the anti-ship missile warhead explosion scaling model constructed in the step five, and supporting the design and optimization of the warhead or protective structure of the related weapon.
The method for calculating the equivalent bare charge size by measuring the quasi-static pressure comprises the following steps of;
in the formula: mEBCCalculating equivalent TNT equivalent, kg, for the shock wave pressure load of the warhead; mchargeIs the TNT equivalent of the actual loading of the warhead, kg; mcaseIs the warhead shell mass, kg.
High speed initial speed V of fragment0The calculation formula is as follows:
in the formula: m issThe warhead housing mass, and E the characteristic energy of the warhead charge.
The characteristic energy E of the charge of the warhead is:
d is the explosive loading detonation velocity of the warhead;
further, the required payload equivalent bare charge size W in step one is:
in the formula: efkFor total kinetic energy of the fragments, EtotalThe total internal energy of the warhead is charged.
Step two, the limit expansion radius R of the warhead0Solved by the following equation:
in the formula: rhofAnd hfThe warhead material density and the shell thickness, ρeIs the charge density, r0Is the initial radius of the warhead housing.
Step three, determining the radius-thickness ratio of the shell according to the failure mode of the warhead:
according to the stress state and the microscopic damage fracture mechanism at the later stage of the deformation of the expansion ring in the warhead,
and the generated stretching and shearing failure forms to determine the required diameter-thickness ratio in the third step:
in the formula: r and R are the outer diameter and the inner diameter of the expanded warhead respectively. Typical failures include ductile necking failures, composite failures, and shear failures.
Selecting a proper material model according to material fracture experiments under different temperatures, different strain rates and various stresses;
further, material deformation is generally divided into four typical phases, an overall plastic phase, a stable necking phase, a local necking development phase, and a final fracture phase. And obtaining the elastic limit, the tensile strength and the yield stress parameters of the material according to the deformation of the four stages, thereby calculating the Young modulus or Poisson ratio of the material and further determining the material required by the step four.
Step five, determining materials according to the equivalent bare charge size obtained in the step one, the limit expansion radius obtained in the step two, the diameter-thickness ratio of the shell obtained in the step three and the material determined in the step four, and according to the characteristic parameters of the explosive load in the cabin, such as the overpressure peak value delta p, the impulse i and the quasi-static pressure p of the shock waveqsAnd the explosive characteristics, air characteristics, structure, material and position parameters are constructed by an impulse similarity law, a proportional law and a dimensional analysis method according to the explosive compression ratio model of the warship-resisting missile warhead.
(1) Explosive parameters: mass W of explosive and density rho of explosiveeEnergy e released by unit mass of explosive explosion and adiabatic index gamma of expansion producte;
(2) Air parameters: initial state pressure p0Air density ρaAir adiabatic index gammaa;
(3) Structural parameters are as follows: the length l, the width b and the height h of the cabin;
(4) position parameters: the distance R between the explosion-facing surface and the explosive.
The length L, the mass M and the time T are taken as basic dimensions, and the dimensions of each physical parameter are shown in Table 1.
TABLE 1 characteristic dimension of explosive load in cabin
Further, according to Π's law, the expression of shock wave pressure and related factors of influence is:
Δp=f(W,ρe,e,γe,p0,γa,l,b,h,ρa,R) (7)
removing influence of factors such as dimensionless quantity and the like, and simplifying into:
Δp=f(W,ρe,e,p0,l,b,h,ρa,R) (8)
with W, rhoeAnd e is a basic quantity, (8) expressed as a dimensionless formula
The explosive types, the warhead shell materials and the air parameters of the scaling model and the original model are the same, and then
(γe,γa,ρe,e,p0,ρa)=const (10)
Then equation (9) is simplified to
According to the similarity rate, if the blast overpressure delta p of the warhead scaling model is required to be reducedmExplosion overpressure delta p of original modelpEqual, then:
in the formula, a subscript p represents a prototype, and a subscript m represents a model. From the law of proportionality, the pressure satisfies the following relationship
Δpp=Δpm (13)
The impulse similarity is:
likewise, it can be derived from the law of proportionality
Obtained by adopting similar analysis method and dimensional analysis method
pqsp=pqsm (16)
And completing the construction of the anti-ship missile warhead explosion scaling model.
Advantageous effects
1. Accuracy and authenticity in anti-ship missile simulation analysis are difficult to guarantee, but full-size anti-ship missiles are high in experiment cost, high in difficulty and long in time consumption. The invention discloses a method for constructing an anti-ship missile warhead explosion scaling model, which obtains an equivalent model of an anti-ship missile warhead explosive and an equivalent model of a shell through theoretical derivation calculation, realizes the test of replacing a full-size missile through the test of the scaling model, reduces the test cost and difficulty and shortens the test period.
2. The invention discloses a method for constructing an anti-ship missile warhead explosion scaling model, which is characterized in that on the basis of the size of an equivalent naked explosive of a warhead, the ultimate expansion radius of the warhead is obtained through calculation according to energy conservation, and Poisson ratio parameters of materials are obtained through material experiments at different temperatures, so that the material selection is determined; and then, carrying out equivalent scaling calculation according to the parameters, and constructing an anti-ship missile warhead explosion scaling model according to the equivalent scaling calculation result.
3. The anti-ship missile warhead explosion scaling model construction method disclosed by the invention can simplify the missile model in simulation, simultaneously keep the original characteristics, weaken the influence of some unnecessary factors on simulation prediction, make the simulation prediction more convenient and make the result more accurate. The invention can save the experiment cost, improve the personnel safety, and provide a convenient, rapid and reliable reference basis for the design and optimization of the warhead or the protection structure of the related weapon and the like.
Detailed Description
The specific steps of the scaling model can be intuitively understood by combining the attached figure 1, and the construction of the scaling model of the explosive missile in the warhead cabin of the anti-ship missile can be completed by implementing the steps. The invention is explained below with reference to the example of a semi-piercing single-vessel missile, which is merely illustrative and not restrictive.
In the embodiment, a scaling model of an explosive missile in a warhead cabin of a semi-penetration single-ship missile is constructed, firstly, a series of required parameters are obtained through measurement and experiments, and the calculation and the description are carried out step by step through five steps. Firstly, air is simplified into an ideal medium, no dissipation exists in the propagation of sound waves in the air, the sound velocity and the density of the air are always kept unchanged, and the sound velocity and the density are respectively assumed to be 346m/s and 1.185kg/m3. Models with the scaling ratios of 1/2, 1/3, 1/4, 1/5, 1/6, 1/8 and 1/10 are respectively established by referring to the original models. Take a scale model of 1/3 scale as an example.
As shown in fig. 1, the method for constructing the anti-ship missile warhead explosion scaling model disclosed in this embodiment specifically includes the following steps:
step 1, after the explosive explosion of the warhead, the internal energy and the kinetic energy of detonation products and the kinetic energy of high-speed fragments are finally converted into quasi-static air pressure to act on a structure, and the equivalent bare explosive size is calculated by measuring the quasi-static pressure; the mass of the fragments of the warhead explosion is 0.055kg, the initial speed of the fragments of the warhead explosion is 2000m/s, and the flying angle of the fragments of the warhead is 3.5 degrees, so that the equivalent TNT equivalent of the warhead is calculated to be 200kg, and the equivalent TNT of the shock wave pressure load is calculated to be 125 kg.
Step 2, assuming that the initial expansion speed is 100m/s in the expansion fragmentation process, wherein the whole expansion process comprises four stages, (1) an integral plasticity stage (0-18.03 mu s): as the expansion ring expands, the whole body deforms uniformly, and generates larger plasticity, and the plastic temperature rise can reach 204K. (2) Stable necking stage (18.03-25.2 μ s): the plastic deformation of local area is increased sharply, and a series of necking is formed; plastic deformation of the non-necked region ceases and the temperature no longer increases (3) the local necking development phase (25.2-29.4 μ s): over time, the partly necked region develops rapidly, the temperature rise further accelerating; a part of the necking is then inhibited (stressed necking) by the unloading of the carriers formed by the adjacent rapidly developing necking zone and finally becomes residual necking in the fragment; (4) in the final fragmentation phase (after 29.4 μ s), local fractures are formed, and the fragmentation formation phase. Neglecting energy loss caused by the streaming effect of the explosive product, and calculating to obtain the ultimate expansion radius of the warhead to be 0.8m according to the energy conservation principle;
step 3, determining that the diameter-thickness ratio of the shell is 0.33 by using the failure mode of the warhead mainly as shear failure and the damage mode as composite damage, wherein the ductile reaming damage and the shear damage of the material both exist;
and 4, selecting a material model at room temperature, obtaining the length of the half-penetration single-ship missile of 0.9m, the Young modulus of 200GPa and the Poisson ratio of 0.25 according to the given material parameters, setting and importing the parameters into ABAQUS for simulation analysis and test.
And 5, substituting the obtained parameters into an impulse similarity law formula and a proportional formula, and reducing the model to 0.2m by calculation, wherein the size is reduced by three times compared with the previous size, and the error is about 3.5%. The experiment results of the scale models of all proportions are shown in figure 2, and the experiment results are basically consistent with the simulation results. In consideration of the precision error of experimental test, the calculation result is well matched with the experimental value in the literature, and the requirement of engineering application can be met. Theoretically, the scaling rule provided by the method can be suitable for model tests under any scaling proportion, and the larger the scaling multiple is, the more the cost is saved. However, in actual tests and simulation analysis, as the scaling factor is increased, the external parameters are changed greatly, which will affect the test precision and increase the error. And with the increase of the scaling factor, the size of the model is gradually reduced, the difficulty of processing and manufacturing is increased, and parameters such as the strength of the material are influenced, so that an appropriate scaling ratio is selected in an actual test.
Step 6: and (5) predicting the mechanical property of the explosion of the anti-ship missile warhead explosion scaling model according to the anti-ship missile warhead explosion scaling model established in the step 5, and supporting the design and optimization of the warhead or the protection structure of the related weapon.
The invention can simplify more complex experiments, and the prior scaling configuration inspection and relevant experimental parameter collection are used for other researches and the like.
The above detailed description is intended to illustrate the objects, aspects and advantages of the present invention, and it should be understood that the above detailed description is only exemplary of the present invention and is not intended to limit the scope of the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.