CN104298848A - Method for analyzing and evaluating reliability of self-destroying function of cannonball for artificial precipitation - Google Patents

Method for analyzing and evaluating reliability of self-destroying function of cannonball for artificial precipitation Download PDF

Info

Publication number
CN104298848A
CN104298848A CN201410335207.4A CN201410335207A CN104298848A CN 104298848 A CN104298848 A CN 104298848A CN 201410335207 A CN201410335207 A CN 201410335207A CN 104298848 A CN104298848 A CN 104298848A
Authority
CN
China
Prior art keywords
self
formula
destruction unit
coyote hole
destruction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410335207.4A
Other languages
Chinese (zh)
Other versions
CN104298848B (en
Inventor
刘卫东
黎济东
占多产
金卫平
高拥军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanchang University
Original Assignee
Nanchang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanchang University filed Critical Nanchang University
Priority to CN201410335207.4A priority Critical patent/CN104298848B/en
Publication of CN104298848A publication Critical patent/CN104298848A/en
Application granted granted Critical
Publication of CN104298848B publication Critical patent/CN104298848B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for analyzing and evaluating the reliability of the self-destroying function of a cannonball for artificial precipitation. By the adoption of the method for analyzing and evaluating the reliability of the self-destroying function of the cannonball for artificial precipitation, operation is easy, and the reliability of the self-destroying function of the cannonball for artificial precipitation can be rapidly and effectively analyzed and evaluated. According to the method for analyzing and evaluating the reliability of the self-destroying function of the cannonball for artificial precipitation, according to the work task mode of the self-destroying function of the gunpowder self-destroying type cannonball for artificial precipitation, an analyzing and evaluating model for the self-destroying reliability of a self-destroying cannonball system for artificial precipitation is achieved through a fixed gunpowder self-destroying unit and a movable gunpowder self-destroying unit, a simulation model for the reliability of the movable gunpowder self-destroying unit is established according to the classic interior ballistics principle, and then the reliability of the self-destroying function of the cannonball for artificial precipitation is designed, analyzed and evaluated.

Description

A kind of analysis for rain bullet self-destroying function reliability and appraisal procedure
Technical field
The present invention relates to a kind of analysis and appraisal procedure of reliability, particularly relate to a kind of analysis and the appraisal procedure that are exclusively used in the reliability of the type self-destroyed rail-proof bomb self-destroying function of increasing rain artificially of gunpowder.
Background technology
Self-destruction is one of the mode of body safe landing of hail-suppression rocket bullet (be called for short rain bullet) of increasing rain artificially, and self-destruction is normally completed by self-destruction subsystem responsible; Traditional self-destruction subsystem comprises head self-destruction unit, middle part self-destruction unit and afterbody self-destruction unit; Head self-destruction unit is fixed on the front end of body in advance; Middle part self-destruction unit and afterbody self-destruction unit then move to assigned address in self-destruction process; But can middle part self-destruction unit and afterbody self-destruction unit reliably be transported to assigned address, be then directly connected to body and realize self-destruction completely; Therefore, carrying out fail-safe analysis to the motion process of middle part self-destruction unit and afterbody self-destruction unit, is the key of correct Prediction rain bullet self-destroying function reliability;
At present, domesticly not yet rain bullet self-destroying function reliability design and analysis and appraisal procedure to be studied, therefore also just do not have definite method to determine rain bullet self-destroying function reliability; For this reason, we are according to the task pattern being realized the type self-destroyed rain bullet self-destroying function of self-destruction by fixing gunpowder self-destruction unit and motion gunpowder self-destruction unit, propose gunpowder type self-destroyed rain bullet fail-safe analysis appraisal procedure, design a kind of analysis of self-destruction rain bullet self-destroying function reliability, assessment models, carried out design and analysis and the assessment of rain bullet self-destroying function reliability by this model.
Summary of the invention
The present invention is in order to make up the deficiencies in the prior art, and provide a kind of analysis for rain bullet self-destroying function reliability and appraisal procedure, the method is simple to operate, fast and effeciently can realize analysis and the assessment of gunpowder type self-destroyed rain bullet self-destroying function reliability.
For achieving the above object, the present invention adopts following technical scheme:
For analysis and the appraisal procedure of rain bullet self-destroying function reliability, it is characterized in that: comprise following concrete steps:
Step one, build the logic diagram model of rain bullet self-destroying function subsystem reliability according to principle of work;
Step 2, analyzing influence realize the key factor of rain bullet self-destroying function reliability, and collect the reliability data of priming system components and parts;
Step 3, to need according to self-destruction unit in the middle part of rain bullet and afterbody self-destruction unit starting mix gaseous-pressure to promote to descend to move to the principle of work of assigned address, based on classical interior ballistics theory, set up the mathematical model of the rain bullet self-destruction course of work;
Step 4, utilize monte carlo method to carry out analog simulation to middle part self-destruction unit and afterbody self-destruction unit motion to the course of work of assigned address, solve middle part self-destruction unit and the afterbody self-destruction unit motion reliability data to assigned address;
The rain bullet self-destroying function reliability block diagram model that step 5, combination build carries out analyzing to rail-proof bomb self-destroying function reliability of increasing rain artificially and assesses;
Further, in step one, described reliability block diagram model comprises gunpowder extension unit r 1, head self-destruction unit r t, middle part self-destruction unit r c, afterbody self-destruction unit r a; Described gunpowder extension unit r 1with head self-destruction unit r tinput end signal connect; Described head self-destruction unit r toutput terminal and middle part self-destruction unit r cinput end signal connect; Described middle part self-destruction unit r coutput terminal and afterbody self-destruction unit r ainput end signal connect;
Wherein, gunpowder extension unit r 1comprise electric igniter r 11, electric igniter r 12, electric igniter r 17, electric igniter r 18, flame agent Dark Ignition tool r 13, flame agent Dark Ignition tool r 14, flame agent r 15, one-level Delay Tube r 16, electric Delay Tube r 19, starting mix r 110; Described electric igniter r 11with flame agent Dark Ignition tool r 13signal connects; Described electric igniter r 12with flame agent Dark Ignition tool r 14signal connects; Described flame agent Dark Ignition tool r 13, flame agent Dark Ignition tool r 14all with flame agent r 15signal connects; Described flame agent r 15with one-level Delay Tube r 16signal connects; Described one-level Delay Tube r 16, electric Delay Tube r 19all and starting mix r 110signal connects; Described electric igniter r 11respectively with electric igniter r 17, electric igniter r 18signal connects;
According to this reliability block diagram, the Reliable Mathematics model of setting rain bullet self-destroying function is:
r= r- 1 r t r c r aformula (1)
Can be obtained the fiduciary level of gunpowder extension unit by self-destruction subsystem reliability logic diagram, its formula of reliability is:
R 1＝{1-{1-[1-(1- R 11 ?R 13)?×?(1- R 12 ?R 14)]× ?R 15× ?R 16}×{1-[1-(1- R 17)?×(1- R 18)]× R 19}}× R 110
Formula (2)
Wherein r 11, r 12, r 17, r 18be expressed as the fiduciary level of electric igniter, r 13, r 14be expressed as the fiduciary level of flame agent Dark Ignition tool, r 15be expressed as the fiduciary level of flame agent, r 16be expressed as the fiduciary level of one-level Delay Tube, r 19be expressed as the fiduciary level of electric Delay Tube, r 110be expressed as the fiduciary level of starting mix; The reliability of middle part self-destruction unit r cand the reliability of afterbody self-destruction unit r a--then that the fiduciary level of these two unit priming systems moves to the comprehensive of the fiduciary level of designated displacement with it, namely
r c= r c1× r cLformula (3)
r a= r a1× r aLformula (4)
Wherein, r c1for the fiduciary level of middle part self-destruction unit priming system, r a1for the fiduciary level of afterbody self-destruction unit priming system.
Further, in step one, according to the principle of work that rain bullet self-destroying function realizes, reliability model structure simplified, described model comprises body, the first coyote hole, the second coyote hole, middle part self-destruction body, afterbody self-destruction body, middle card and nozzle hole; Described first coyote hole is arranged on body foremost, and arranges near middle part self-destruction body; Second coyote hole is arranged on the middle part of middle part self-destruction body and afterbody self-destruction body; Described middle card is arranged on the middle part of body; Described nozzle hole is arranged on the tail end of body;
Further, in step 2, key factor impact being realized to rain bullet self-destroying function reliability is analyzed, and mainly comprises the following aspects:
(1) it is the reliabilities of two self-destruction unit motions to assigned address that the reliability block diagram, set up by step one determines to affect the key factor realizing rain bullet self-destroying function reliability;
(2) reliability model, set up according to step one clearly needs the priming system component reliability data of collecting;
Further, in step 3, constructed mathematical model is mainly based on following hypothesis:
(1), in self-destruction unit whether priming system transmittance process at different levels is reliable, and whether self-destruction body brisance meets the demands, and namely body can be scraped completely when assigned address detonates and need be determined by the reliability of priming system element;
(2), self-destruction unit motion to the process of assigned address meets the relevant hypothesis of classical interior ballistics model;
(3), due to the small aperture being communicated with two coyote holes minimum, its circulated gases minimum relative to coyote hole gas gross thus hypothesis aperture in flow without gas;
(4), by the process that starting mix is pushed to assigned address, self-destruction unit do not considered that the viscosity of gas, gas are to the friction of body inwall, and combustion gas adopts indirect method correction to the heat transfer of chamber wall;
(5), powder burning and two self-destruction unit motions all carry out under mean pressure;
(6) in the middle part of, when considering that middle card blocks middle part self-destruction unit, self-destruction unit speed is minimum, also minimum to ballistic impact, therefore approximate think the heat that produces in this process to the second coyote hole without impact;
(7), due to starting mix burning and self-destruction unit motion process time very short, much smaller than the defer time of Delay Tube, thus its self-destruction time and height only relevant with the outer trajectory that defer time and the rain bullet of Delay Tube are launched, ignoring it at this affects.
Constructed mathematical model mainly comprises following six stages, specific as follows:
1, the first stage, to take fire afterbody self-destruction unit setting in motion from the gunpowder of two coyote holes (i.e. the first coyote hole, the second coyote hole), this process is coyote hole constant volume combustion process; The involved formula of whole process is as follows:
combinatorial formula (5)
This is a composite class formula, and wherein, first formula is the second coyote hole combustion gas growth equation formula, and second formula is the second coyote hole burning rate equation, and the 3rd formula is the second coyote hole pressure equation formula;
In above formula, p 2be the pressure of the second coyote hole, v 0' be the original volume of the second coyote hole; ω 2it is the explosive load of the second coyote hole.With starting pressure P 0for initial value, solve the second coyote hole powder burning situation;
2, subordinate phase, from afterbody self-destruction unit and the time of releasing of middle part self-destruction unit to middle part self-destruction unit setting in motion, this process is the second coyote hole burning and external expansion work, the first coyote hole constant volume combustion process; The involved formula of whole process is as follows:
combinatorial formula (6)
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, second formula is the second coyote hole combustion gas growth equation formula, 3rd formula is the first coyote hole burning rate equation, 4th formula is the second coyote hole combustion gas growth equation formula, 5th formula is the first coyote hole pressure equation formula, 6th formula is the second coyote hole gas equation, 7th formula is afterbody self-destruction unit momentum equation, 8th formula is afterbody self-destruction unit motion equation, 9th formula is the second coyote hole energy conservation equation formula, tenth formula is for calculating end condition equation, 11 formula is the calculation equation of the enthalpy that the second coyote hole unit mass gas flows out,
In above formula, p 1be the pressure of the first coyote hole, v 0it is the original volume of the first coyote hole; ω 1be the explosive load of the first coyote hole, w 2be the gaseous mass percentage that the second coyote hole is revealed due to the gap between afterbody self-destruction unit and body, t 2be the gas temperature in the second coyote hole, l 2for the displacement of afterbody self-destruction unit relative middle part self-destruction unit, φ 2for the secondary merit design factor of afterbody self-destruction unit, m 2for the quality of afterbody self-destruction unit, v 2for the speed of afterbody self-destruction unit, s 2for the sectional area of afterbody self-destruction unit, p 0for external atmospheric pressure, r g =R/M, M are the molal weight of starting mix burning gases;
3, the phase III, from middle part self-destruction unit setting in motion to the second coyote hole gunpowder burns completely, this stage, two coyote hole gunpowder all burnt and expansion work; The involved formula of whole process is as follows:
combinatorial formula (7)
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, second formula is the first coyote hole burning rate equation, 3rd formula is the second coyote hole combustion gas growth equation formula, 4th formula is the first coyote hole gas equation, 5th formula is the second coyote hole gas equation, 6th formula is the momentum equation of middle part self-destruction unit, 7th formula is middle part self-destruction unit motion equation, 8th formula is the momentum equation of afterbody self-destruction unit, 9th formula is afterbody self-destruction unit motion equation, tenth formula is the energy conservation equation formula of the first coyote hole, 11 formula is the energy conservation equation formula of the second coyote hole, 12 formula is the calculation equation of the enthalpy that the first coyote hole unit mass gas flows out,
In above formula, w 1be the gaseous mass percentage that the first coyote hole is revealed due to the gap between middle part self-destruction unit and body, t 1be the gas temperature in the first coyote hole, l 1for the displacement of the relative body of middle part self-destruction unit, φ 1for the secondary merit design factor of middle part self-destruction unit, m 1for the quality of middle part self-destruction unit, v 1for the speed of middle part self-destruction unit, s 1for the sectional area of middle part self-destruction unit, h 1it is the enthalpy that the first coyote hole unit mass gas flows out;
4, fourth stage, is burned to till the first coyote hole starting mix burns completely from the second coyote hole starting mix, this stage second coyote hole expansion work completely, the first coyote hole burning and expansion work; The involved formula of whole process is as follows:
combinatorial formula (8)
In this combinatorial formula, first formula is the first coyote hole burning rate equation, second formula is the first coyote hole gas equation, 3rd formula is the second coyote hole gas equation, 4th formula is the energy conservation equation formula of the first coyote hole, 5th formula is the energy conservation equation formula of the second coyote hole, 6th formula is the momentum equation of middle part self-destruction unit, 7th formula is middle part self-destruction unit motion equation, 8th formula is the momentum equation of afterbody self-destruction unit, 9th formula is afterbody self-destruction unit motion equation,
5, five-stage, burns to assigned address, now two coyote hole expansion work in the middle part of middle part self-destruction unit motion to body from the first coyote hole starting mix completely; The involved formula of whole process is as follows:
combinatorial formula (9)
In this combinatorial formula, first formula is the first coyote hole gas equation, second formula is the second coyote hole gas equation, 3rd formula is the energy conservation equation formula of the first coyote hole, 4th formula is the energy conservation equation formula of the second coyote hole, 5th formula is the momentum equation of middle part self-destruction unit, 6th formula is middle part self-destruction unit motion equation, 7th formula is the momentum equation of afterbody self-destruction unit, and the 8th formula is afterbody self-destruction unit motion equation; The motion conditions that can calculate middle part self-destruction unit is solved by system of equations;
6, the 6th stage from assigned address in the middle part of middle part self-destruction unit motion to body to afterbody self-destruction unit motion to afterbody assigned address moment, in the middle part of this stage, self-destruction unit is no longer transported to because the associated mechanisms of body medium position blocks, first coyote hole gas is no longer revealed, the still motion under gaseous-pressure promotes of afterbody self-destruction unit; The involved formula of whole process is as follows:
combinatorial formula (10)
In this combinatorial formula, first formula is the second coyote hole gas equation, second formula is the equation of momentum of afterbody self-destruction unit, 3rd formula is afterbody self-destruction unit motion equation, 4th formula is the energy conservation equation formula of the second coyote hole, is gone out the motion conditions of afterbody self-destruction unit by solving equations;
In above-mentioned six system of equations:
(1), the first coyote hole flows into the gas-flow equation of the second coyote hole:
formula (11)
In formula , φ 3for Flux Loss coefficient, s 1' be the sectional area in middle part self-destruction unit and body gap, m 1' be the gaseous mass of the first coyote hole;
(2), the second coyote hole flows into extraneous gas-flow equation:
formula (12)
In formula , φ 4for Flux Loss coefficient, s 2' be the sectional area in afterbody self-destruction unit and body gap, m 2' be the gaseous mass of the second coyote hole;
Further, in step 4, after the principal element considering impact middle part and afterbody self-destruction unit motion, utilize monte carlo method to middle part self-destruction unit motion to middle card position, afterbody self-destruction unit motion is analyzed to the reliability of nozzle hole position, and concrete steps are as follows:
A, will realize self-destroying function, middle part self-destruction unit can be pushed into card position in body by demand fulfillment starting mix, afterbody self-destruction unit is pushed into nozzle hole position simultaneously, if in the middle part of body under assigned address to original state the distance of middle part self-destruction unit bottom be y 1, afterbody assigned address to the distance of original state lower tail self-destruction unit bottom is y 2;
In actual production process, due to manufacturing accuracy impact, y 1with y 2for obeying the parameter of certain regularity of distribution, and totally y 1distribution function be f 1 y ( y), totally y 2distribution function be f 2 y ( y);
B, middle part self-destruction unit relative body displacement can be drawn according to above-mentioned mathematical model by emulation l 1a sample and the relative body displacement of afterbody self-destruction unit l 2a sample, and l 1distribution function be f 1 z ( z), totally l 2distribution function be f 2 z ( z); Middle part, afterbody self-destruction unit can arrive in the middle part of body, afterbody assigned address need meet l 1>= y 1, and l 2>= y 2this condition;
C, determine parameters distribution function after, then can determine the fiduciary level of assigned address in the middle part of middle part self-destruction unit motion to body, computing formula is:
formula (13)
The fiduciary level that afterbody self-destruction unit motion is specified to afterbody, computing formula is:
formula (14)
Affect the powder charge amount that self-destruction unit motion to the topmost factor of assigned address is the first coyote hole and the second coyote hole ω 1, ω 2and the quality of middle part self-destruction unit and afterbody self-destruction unit m 1, m 2;
In a practical situation, explosive load ω i ( i=1,2) and each self-destruction element quality m j ( j=1,2) all there is certain span and statistical nature: ω ip -△ ω i ≤ ω i ≤ ω ip + △ ω i formula (15)
m jp -△ m j ≤ m j ≤ m jp + △ m j formula (16)
In formula ω ip for ω i nominal value, △ ω i for ω i tolerance, m jp for m j nominal value, △ m j for m j tolerance;
D, according to set up the mathematical model based on interior ballistics theory, the function that the displacement of self-destruction unit relative body is relevant to two the chamber charges and two self-destruction element qualities can be obtained:
l 1= d( ω i , m j ) formula (17)
l 2= d( ω i , m j ) formula (18)
Formula (17) is with the displacement of the relative body of middle part self-destruction unit l 1for the system of equations of unknown function, formula (18) is the displacement of the relative body of afterbody self-destruction unit l 2for the system of equations of unknown function;
E, for obtaining l 1, l 2distribution function, emulated by Monte Carlo, carry out random sampling according to the requirement that the quality of actual conditions to the explosive load of two coyote holes and middle part self-destruction unit and afterbody self-destruction unit is checked on by its regularity of distribution of obeying and quality inspection; As kduring secondary sampling, extract on request ω ik ( i=1,2) and m jk ( j=1,2); Will ω ik with m jk substitute into (17), (18) formula (i.e. solving equation group) can obtain respectively in the middle part of the operation sample value of self-destruction unit and the displacement of the relative body of afterbody self-destruction unit:
l 1 k = d( ω 1 k , ω 2 k ; m 1 k , m 2 k ) formula (19)
l 2 k = d( ω 1 k , ω 2 k ; m 1 k , m 2 k ) formula (20)
After realistic model runs N time, the shift value of the relative body of N number of middle part self-destruction unit can be obtained l 1 k ( k=1,2 ..., N) and the shift value of the relative body of N number of afterbody self-destruction unit l 2 k (k=1,2 ..., N); Respectively with N number of l 1 k value is as the relative body displacement of middle part self-destruction unit l 1a sample, with N number of l 2 k value is as the relative body displacement of afterbody self-destruction unit l 2a sample, draw the frequency histogram of two parameters, and guestimate or suppose overall l 1distribution function f 1 z ( z), totally l 2distribution function f 2 z ( z); By sample data, hypothesis is tested, to judge the correctness supposed; Finally determine distribution function f 1 z ( z) and distribution function f 2 z ( z), can determine in the middle part of rain bullet by (13) (14) formula and afterbody self-destruction unit motion to the fiduciary level of assigned address;
Further, in step e, the process that Monte Carlo emulates, concrete steps are as follows:
Step (one), clear and definite sampling parametric ω 1, ω 2with m 1, m 2and simulation times N;
Step (two), Monte Carlo is utilized to carry out random sampling to parameter;
Step (three), by parameter ω 1, ω 2with m 1, m 2and simulation times N substitutes into formula (17), (18) solve;
Step (four), output l 1 k , l 2k and stored in array l 1, l 2; N in this process>= k;
Step (five), output array l 1, l 2, and draw distribution function;
Compared with prior art, adopt such scheme, the invention has the beneficial effects as follows: the present invention is according to the task pattern of gunpowder type self-destroyed rain bullet self-destroying function, propose the rain bullet system self-destruction fail-safe analysis assessment models being realized self-destruction by fixing gunpowder self-destruction unit and motion gunpowder self-destruction unit, use classical interior ballistics principle to construct the realistic model of motion gunpowder self-destruction unit reliability, and then design and analysis and assessment are carried out to the reliability of rain bullet self-destroying function.
Accompanying drawing explanation
Fig. 1 is appraisal procedure schematic flow sheet of the present invention;
Fig. 2 is reliability block diagram model structure schematic diagram of the present invention;
Fig. 3 is reliability model structural representation of the present invention;
Fig. 4 is Monte Carlo simulation flow schematic diagram in the present invention.
Embodiment
Below in conjunction with drawings and Examples, the present invention is further described.
As shown in Figure 1, a kind of analysis for rain bullet self-destroying function reliability and appraisal procedure, is characterized in that: comprise following concrete steps:
Step one, build the logic diagram model of rain bullet self-destroying function subsystem reliability according to principle of work;
Step 2, analyzing influence realize the key factor of rain bullet self-destroying function reliability, and collect the reliability data of priming system components and parts;
Step 3, to need according to self-destruction unit in the middle part of rain bullet and afterbody self-destruction unit starting mix gaseous-pressure to promote to descend to move to the principle of work of assigned address, based on classical interior ballistics theory, set up the mathematical model of the rain bullet self-destruction course of work;
Step 4, utilize monte carlo method to carry out analog simulation to middle part self-destruction unit and afterbody self-destruction unit motion to the course of work of assigned address, solve middle part self-destruction unit and the afterbody self-destruction unit motion reliability data to assigned address;
The rain bullet self-destroying function reliability block diagram model that step 5, combination build carries out analyzing to rail-proof bomb self-destroying function reliability of increasing rain artificially and assesses;
Preferably, as shown in Figure 2, in step one, described reliability block diagram model comprises gunpowder extension unit r 1, head self-destruction unit r t, middle part self-destruction unit r c, afterbody self-destruction unit r a; Described gunpowder extension unit r 1with head self-destruction unit r tinput end signal connect; Described head self-destruction unit r toutput terminal and middle part self-destruction unit r cinput end signal connect; Described middle part self-destruction unit r coutput terminal and afterbody self-destruction unit r ainput end signal connect;
Wherein, gunpowder extension unit r 1comprise electric igniter r 11, electric igniter r 12, electric igniter r 17, electric igniter r 18, flame agent Dark Ignition tool r 13, flame agent Dark Ignition tool r 14, flame agent r 15, one-level Delay Tube r 16, electric Delay Tube r 19, starting mix r 110; Described electric igniter r 11with flame agent Dark Ignition tool r 13signal connects; Described electric igniter r 12with flame agent Dark Ignition tool r 14signal connects; Described flame agent Dark Ignition tool r 13, flame agent Dark Ignition tool r 14all with flame agent r 15signal connects; Described flame agent r 15with one-level Delay Tube r 16signal connects; Described one-level Delay Tube r 16, electric Delay Tube r 19all and starting mix r 110signal connects; Described electric igniter r 11respectively with electric igniter r 17, electric igniter r 18signal connects;
According to this reliability block diagram, the Reliable Mathematics model of setting rain bullet self-destroying function is:
r= r- 1 r t r c r aformula (1)
Can be obtained the fiduciary level of gunpowder extension unit by self-destruction subsystem reliability logic diagram, its formula of reliability is:
R 1＝{1-{1-[1-(1- R 11 ?R 13)?×?(1- R 12 ?R 14)]× ?R 15× ?R 16}×{1-[1-(1- R 17)?×(1- R 18)]× R 19}}× R 110
Formula (2)
Wherein r 11, r 12, r 17, r 18be expressed as the fiduciary level of electric igniter, r 13, r 14be expressed as the fiduciary level of flame agent Dark Ignition tool, r 15be expressed as the fiduciary level of flame agent, r 16be expressed as the fiduciary level of one-level Delay Tube, r 19be expressed as the fiduciary level of electric Delay Tube, r 110be expressed as the fiduciary level of starting mix; The reliability of middle part self-destruction unit r cand the reliability of afterbody self-destruction unit r a--then that the fiduciary level of these two unit priming systems moves to the comprehensive of the fiduciary level of designated displacement with it, namely
r c= r c1× r cLformula (3)
r a= r a1× r aLformula (4)
Wherein, r c1for the fiduciary level of middle part self-destruction unit priming system, r a1for the fiduciary level of afterbody self-destruction unit priming system.
Preferably, as shown in Figure 3, in step one, according to the principle of work that rain bullet self-destroying function realizes, reliability model structure simplified, described model comprises body 5, first coyote hole 1, second coyote hole 3, middle part self-destruction body 2, afterbody self-destruction body 4, middle card 6 and nozzle hole 7; Described first coyote hole 1 is arranged on body 5 foremost, and arranges near middle part self-destruction body 2; Second coyote hole 3 is arranged on the middle part of middle part self-destruction body 2 and afterbody self-destruction body 4; Described middle card 6 is arranged on the middle part of body 5; Described nozzle hole 7 is arranged on the tail end of body 5;
Preferably, in step 2, key factor impact being realized to rain bullet self-destroying function reliability is analyzed, and mainly comprises the following aspects:
(1) it is the reliabilities of two self-destruction unit motions to assigned address that the reliability block diagram, set up by step one determines to affect the key factor realizing rain bullet self-destroying function reliability;
(2) reliability model, set up according to step one clearly needs the priming system component reliability data of collecting;
Preferably, in step 3, constructed mathematical model is mainly based on following hypothesis:
(1), in self-destruction unit whether priming system transmittance process at different levels is reliable, and whether self-destruction body brisance meets the demands, and namely body can be scraped completely when assigned address detonates and need be determined by the reliability of priming system element;
(2), self-destruction unit motion to the process of assigned address meets the relevant hypothesis of classical interior ballistics model;
(3), due to the small aperture being communicated with two coyote holes minimum, its circulated gases minimum relative to coyote hole gas gross thus hypothesis aperture in flow without gas;
(4), by the process that starting mix is pushed to assigned address, self-destruction unit do not considered that the viscosity of gas, gas are to the friction of body inwall, and combustion gas adopts indirect method correction to the heat transfer of chamber wall;
(5), powder burning and two self-destruction unit motions all carry out under mean pressure;
(6) in the middle part of, when considering that middle card blocks middle part self-destruction unit, self-destruction unit speed is minimum, also minimum to ballistic impact, therefore approximate think the heat that produces in this process to the second coyote hole without impact;
(7), due to starting mix burning and self-destruction unit motion process time very short, much smaller than the defer time of Delay Tube, thus its self-destruction time and height only relevant with the outer trajectory that defer time and the rain bullet of Delay Tube are launched, ignoring it at this affects.
Constructed mathematical model mainly comprises following six stages, specific as follows:
1, the first stage, to take fire afterbody self-destruction unit setting in motion from the gunpowder of two coyote holes (i.e. the first coyote hole, the second coyote hole), this process is coyote hole constant volume combustion process; The involved formula of whole process is as follows:
combinatorial formula (5)
This is a composite class formula, and wherein, first formula is the second coyote hole combustion gas growth equation formula, and second formula is the second coyote hole burning rate equation, and the 3rd formula is the second coyote hole pressure equation formula;
In above formula, p 2be the pressure of the second coyote hole, v 0' be the original volume of the second coyote hole; ω 2it is the explosive load of the second coyote hole.With starting pressure P 0for initial value, solve the second coyote hole powder burning situation;
2, subordinate phase, from afterbody self-destruction unit and the time of releasing of middle part self-destruction unit to middle part self-destruction unit setting in motion, this process is the second coyote hole burning and external expansion work, the first coyote hole constant volume combustion process; The involved formula of whole process is as follows:
combinatorial formula (6)
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, second formula is the second coyote hole combustion gas growth equation formula, 3rd formula is the first coyote hole burning rate equation, 4th formula is the second coyote hole combustion gas growth equation formula, 5th formula is the first coyote hole pressure equation formula, 6th formula is the second coyote hole gas equation, 7th formula is afterbody self-destruction unit momentum equation, 8th formula is afterbody self-destruction unit motion equation, 9th formula is the second coyote hole energy conservation equation formula, tenth formula is for calculating end condition equation, 11 formula is the calculation equation of the enthalpy that the second coyote hole unit mass gas flows out,
In above formula, p 1be the pressure of the first coyote hole, v 0it is the original volume of the first coyote hole; ω 1be the explosive load of the first coyote hole, w 2be the gaseous mass percentage that the second coyote hole is revealed due to the gap between afterbody self-destruction unit and body, t 2be the gas temperature in the second coyote hole, l 2for the displacement of afterbody self-destruction unit relative middle part self-destruction unit, φ 2for the secondary merit design factor of afterbody self-destruction unit, m 2for the quality of afterbody self-destruction unit, v 2for the speed of afterbody self-destruction unit, s 2for the sectional area of afterbody self-destruction unit, p 0for external atmospheric pressure, r g =R/M, M are the molal weight of starting mix burning gases;
3, the phase III, from middle part self-destruction unit setting in motion to the second coyote hole gunpowder burns completely, this stage, two coyote hole gunpowder all burnt and expansion work; The involved formula of whole process is as follows:
combinatorial formula (7)
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, second formula is the first coyote hole burning rate equation, 3rd formula is the second coyote hole combustion gas growth equation formula, 4th formula is the first coyote hole gas equation, 5th formula is the second coyote hole gas equation, 6th formula is the momentum equation of middle part self-destruction unit, 7th formula is middle part self-destruction unit motion equation, 8th formula is the momentum equation of afterbody self-destruction unit, 9th formula is afterbody self-destruction unit motion equation, tenth formula is the energy conservation equation formula of the first coyote hole, 11 formula is the energy conservation equation formula of the second coyote hole, 12 formula is the calculation equation of the enthalpy that the first coyote hole unit mass gas flows out,
In above formula, w 1be the gaseous mass percentage that the first coyote hole is revealed due to the gap between middle part self-destruction unit and body, t 1be the gas temperature in the first coyote hole, l 1for the displacement of the relative body of middle part self-destruction unit, φ 1for the secondary merit design factor of middle part self-destruction unit, m 1for the quality of middle part self-destruction unit, v 1for the speed of middle part self-destruction unit, s 1for the sectional area of middle part self-destruction unit, h 1it is the enthalpy that the first coyote hole unit mass gas flows out;
4, fourth stage, is burned to till the first coyote hole starting mix burns completely from the second coyote hole starting mix, this stage second coyote hole expansion work completely, the first coyote hole burning and expansion work; The involved formula of whole process is as follows:
combinatorial formula (8)
In this combinatorial formula, first formula is the first coyote hole burning rate equation, second formula is the first coyote hole gas equation, 3rd formula is the second coyote hole gas equation, 4th formula is the energy conservation equation formula of the first coyote hole, 5th formula is the energy conservation equation formula of the second coyote hole, 6th formula is the momentum equation of middle part self-destruction unit, 7th formula is middle part self-destruction unit motion equation, 8th formula is the momentum equation of afterbody self-destruction unit, 9th formula is afterbody self-destruction unit motion equation,
5, five-stage, burns to assigned address, now two coyote hole expansion work in the middle part of middle part self-destruction unit motion to body from the first coyote hole starting mix completely; The involved formula of whole process is as follows:
combinatorial formula (9)
In this combinatorial formula, first formula is the first coyote hole gas equation, second formula is the second coyote hole gas equation, 3rd formula is the energy conservation equation formula of the first coyote hole, 4th formula is the energy conservation equation formula of the second coyote hole, 5th formula is the momentum equation of middle part self-destruction unit, 6th formula is middle part self-destruction unit motion equation, 7th formula is the momentum equation of afterbody self-destruction unit, and the 8th formula is afterbody self-destruction unit motion equation; The motion conditions that can calculate middle part self-destruction unit is solved by system of equations;
6, the 6th stage from assigned address in the middle part of middle part self-destruction unit motion to body to afterbody self-destruction unit motion to afterbody assigned address moment, in the middle part of this stage, self-destruction unit is no longer transported to because the associated mechanisms of body medium position blocks, first coyote hole gas is no longer revealed, the still motion under gaseous-pressure promotes of afterbody self-destruction unit; The involved formula of whole process is as follows:
combinatorial formula (10)
In this combinatorial formula, first formula is the second coyote hole gas equation, second formula is the equation of momentum of afterbody self-destruction unit, 3rd formula is afterbody self-destruction unit motion equation, 4th formula is the energy conservation equation formula of the second coyote hole, is gone out the motion conditions of afterbody self-destruction unit by solving equations;
In above-mentioned six system of equations:
(1), the first coyote hole flows into the gas-flow equation of the second coyote hole:
formula (11)
In formula , φ 3for Flux Loss coefficient, s 1' be the sectional area in middle part self-destruction unit and body gap, m 1' be the gaseous mass of the first coyote hole;
(2), the second coyote hole flows into extraneous gas-flow equation:
formula (12)
In formula , φ 4for Flux Loss coefficient, s 2' be the sectional area in afterbody self-destruction unit and body gap, m 2' be the gaseous mass of the second coyote hole;
Preferably, in step 4, after the principal element considering impact middle part and afterbody self-destruction unit motion, utilize monte carlo method to middle part self-destruction unit motion to middle card position, afterbody self-destruction unit motion is analyzed to the reliability of nozzle hole position, and concrete steps are as follows:
A, will realize self-destroying function, middle part self-destruction unit can be pushed into card position in body by demand fulfillment starting mix, afterbody self-destruction unit is pushed into nozzle hole position simultaneously, if in the middle part of body under assigned address to original state the distance of middle part self-destruction unit bottom be y 1, afterbody assigned address to the distance of original state lower tail self-destruction unit bottom is y 2;
In actual production process, due to manufacturing accuracy impact, y 1with y 2for obeying the parameter of certain regularity of distribution, and totally y 1distribution function be f 1 y ( y), totally y 2distribution function be f 2 y ( y);
B, middle part self-destruction unit relative body displacement can be drawn according to above-mentioned mathematical model by emulation l 1a sample and the relative body displacement of afterbody self-destruction unit l 2a sample, and l 1distribution function be f 1 z ( z), totally l 2distribution function be f 2 z ( z); Middle part, afterbody self-destruction unit can arrive in the middle part of body, afterbody assigned address need meet l 1>= y 1, and l 2>= y 2this condition;
C, determine parameters distribution function after, then can determine the fiduciary level of assigned address in the middle part of middle part self-destruction unit motion to body, computing formula is:
formula (13)
The fiduciary level that afterbody self-destruction unit motion is specified to afterbody, computing formula is:
formula (14)
Affect the powder charge amount that self-destruction unit motion to the topmost factor of assigned address is the first coyote hole and the second coyote hole ω 1, ω 2and the quality of middle part self-destruction unit and afterbody self-destruction unit m 1, m 2;
In a practical situation, explosive load ω i ( i=1,2) and each self-destruction element quality m j ( j=1,2) all there is certain span and statistical nature: ω ip -△ ω i ≤ ω i ≤ ω ip + △ ω i formula (15)
m jp -△ m j ≤ m j ≤ m jp + △ m j formula (16)
In formula ω ip for ω i nominal value, △ ω i for ω i tolerance, m jp for m j nominal value, △ m j for m j tolerance;
D, according to set up the mathematical model based on interior ballistics theory, the function that the displacement of self-destruction unit relative body is relevant to two the chamber charges and two self-destruction element qualities can be obtained:
l 1= d( ω i , m j ) formula (17)
l 2= d( ω i , m j ) formula (18)
Formula (17) is with the displacement of the relative body of middle part self-destruction unit l 1for the system of equations of unknown function, formula (18) is the displacement of the relative body of afterbody self-destruction unit l 2for the system of equations of unknown function;
E, for obtaining l 1, l 2distribution function, emulated by Monte Carlo, carry out random sampling according to the requirement that the quality of actual conditions to the explosive load of two coyote holes and middle part self-destruction unit and afterbody self-destruction unit is checked on by its regularity of distribution of obeying and quality inspection; As kduring secondary sampling, extract on request ω ik ( i=1,2) and m jk ( j=1,2); Will ω ik with m jk substitute into (17), (18) formula (i.e. solving equation group) can obtain respectively in the middle part of the operation sample value of self-destruction unit and the displacement of the relative body of afterbody self-destruction unit:
l 1 k = d( ω 1 k , ω 2 k ; m 1 k , m 2 k ) formula (19)
l 2 k = d( ω 1 k , ω 2 k ; m 1 k , m 2 k ) formula (20)
After realistic model runs N time, the shift value of the relative body of N number of middle part self-destruction unit can be obtained l 1 k ( k=1,2 ..., N) and the shift value of the relative body of N number of afterbody self-destruction unit l 2 k (k=1,2 ..., N); Respectively with N number of l 1 k value is as the relative body displacement of middle part self-destruction unit l 1a sample, with N number of l 2 k value is as the relative body displacement of afterbody self-destruction unit l 2a sample, draw the frequency histogram of two parameters, and guestimate or suppose overall l 1distribution function f 1 z ( z), totally l 2distribution function f 2 z ( z); By sample data, hypothesis is tested, to judge the correctness supposed; Finally determine distribution function f 1 z ( z) and distribution function f 2 z ( z), can determine in the middle part of rain bullet by (13) (14) formula and afterbody self-destruction unit motion to the fiduciary level of assigned address;
Preferably, as shown in Figure 2, in step e, the process that Monte Carlo emulates, concrete steps are as follows:
Step (one), clear and definite sampling parametric ω 1, ω 2with m 1, m 2and simulation times N;
Step (two), Monte Carlo is utilized to carry out random sampling to parameter;
Step (three), by parameter ω 1, ω 2with m 1, m 2and simulation times N substitutes into formula (17), (18) solve;
Step (four), output l 1 k , l 2k and stored in array l 1, l 2; N in this process>= k;
Step (five), output array l 1, l 2, and draw distribution function;
Determine in the middle part of rain bullet and after afterbody self-destruction unit motion to the fiduciary level of assigned address, the reliability logic block diagram in conjunction with rain bullet self-destroying function can solve the fiduciary level of rain bullet self-destroying function, completes the reliability assessment to rain bullet self-destroying function; In conjunction with assessment result, related opinions is proposed to the reliability design of rain bullet self-destroying function;
The present invention is according to the task pattern of gunpowder type self-destroyed rain bullet self-destroying function, propose the rain bullet system self-destruction fail-safe analysis assessment models being realized self-destruction by fixing gunpowder self-destruction unit and motion gunpowder self-destruction unit, use classical interior ballistics principle to construct the realistic model of motion gunpowder self-destruction unit reliability, and then design and analysis and assessment are carried out to the reliability of rain bullet self-destroying function;
The present invention is not limited to above-mentioned concrete embodiment, and those of ordinary skill in the art is from above-mentioned design, and without performing creative labour, done all conversion, all drop within protection scope of the present invention.

Claims (7)

1., for analysis and the appraisal procedure of rain bullet self-destroying function reliability, it is characterized in that: comprise following concrete steps:
Step one, build the logic diagram model of rain bullet self-destroying function subsystem reliability according to principle of work;
Step 2, analyzing influence realize the key factor of rain bullet self-destroying function reliability, and collect the reliability data of priming system components and parts;
Step 3, to need according to self-destruction unit in the middle part of rain bullet and afterbody self-destruction unit starting mix gaseous-pressure to promote to descend to move to the principle of work of assigned address, based on classical interior ballistics theory, set up the mathematical model of the rain bullet self-destruction course of work;
Step 4, utilize monte carlo method to carry out analog simulation to middle part self-destruction unit and afterbody self-destruction unit motion to the course of work of assigned address, solve middle part self-destruction unit and the afterbody self-destruction unit motion reliability data to assigned address;
The rain bullet self-destroying function reliability block diagram model that step 5, combination build carries out analyzing to rail-proof bomb self-destroying function reliability of increasing rain artificially and assesses.
2. a kind of analysis for rain bullet self-destroying function reliability according to claim 1 and appraisal procedure, it is characterized in that: in step one, described reliability block diagram model comprises gunpowder extension unit r 1, head self-destruction unit r t, middle part self-destruction unit r c, afterbody self-destruction unit r a; Described gunpowder extension unit r 1with head self-destruction unit r tinput end signal connect; Described head self-destruction unit r toutput terminal and middle part self-destruction unit r cinput end signal connect; Described middle part self-destruction unit r coutput terminal and afterbody self-destruction unit r ainput end signal connect;
Wherein, gunpowder extension unit r 1comprise electric igniter r 11, electric igniter r 12, electric igniter r 17, electric igniter r 18, flame agent Dark Ignition tool r 13, flame agent Dark Ignition tool r 14, flame agent r 15, one-level Delay Tube r 16, electric Delay Tube r 19, starting mix r 110; Described electric igniter r 11with flame agent Dark Ignition tool r 13signal connects; Described electric igniter r 12with flame agent Dark Ignition tool r 14signal connects; Described flame agent Dark Ignition tool r 13, flame agent Dark Ignition tool r 14all with flame agent r 15signal connects; Described flame agent r 15with one-level Delay Tube r 16signal connects; Described one-level Delay Tube r 16, electric Delay Tube r 19all and starting mix r 110signal connects; Described electric igniter r 11respectively with electric igniter r 17, electric igniter r 18signal connects;
According to this reliability block diagram, the Reliable Mathematics model of setting rain bullet self-destroying function is:
r= r- 1 r t r c r aformula (1)
Can be obtained the fiduciary level of gunpowder extension unit by self-destruction subsystem reliability logic diagram, its formula of reliability is:
R 1＝{1-{1-[1-(1- R 11 ?R 13)?×?(1- R 12 ?R 14)]× ?R 15× ?R 16}×{1-[1-(1- R 17)?×(1- R 18)]× R 19}}× R 110
Formula (2)
Wherein r 11, r 12, r 17, r 18be expressed as the fiduciary level of electric igniter, r 13, r 14be expressed as the fiduciary level of flame agent Dark Ignition tool, r 15be expressed as the fiduciary level of flame agent, r 16be expressed as the fiduciary level of one-level Delay Tube, r 19be expressed as the fiduciary level of electric Delay Tube, r 110be expressed as the fiduciary level of starting mix; The reliability of middle part self-destruction unit r cand the reliability of afterbody self-destruction unit r a--then that the fiduciary level of these two unit priming systems moves to the comprehensive of the fiduciary level of designated displacement with it, namely
r c= r c1× r cLformula (3)
r a= r a1× r aLformula (4)
Wherein, r c1for the fiduciary level of middle part self-destruction unit priming system, r a1for the fiduciary level of afterbody self-destruction unit priming system.
3. a kind of analysis for rain bullet self-destroying function reliability according to claim 1 and appraisal procedure, it is characterized in that: the principle of work realized according to rain bullet self-destroying function, reliability model structure simplified, described model comprises body, the first coyote hole, the second coyote hole, middle part self-destruction body, afterbody self-destruction body, middle card and nozzle hole; Described first coyote hole is arranged on body foremost, and arranges near middle part self-destruction body; Second coyote hole is arranged on the middle part of middle part self-destruction body and afterbody self-destruction body; Described middle card is arranged on the middle part of body; Described nozzle hole is arranged on the tail end of body.
4. a kind of analysis for rain bullet self-destroying function reliability according to claim 1 and appraisal procedure, is characterized in that: in step 2, and key factor impact being realized to rain bullet self-destroying function reliability is analyzed, and mainly comprises the following aspects:
(1) it is the reliabilities of two self-destruction unit motions to assigned address that the reliability block diagram, set up by step one determines to affect the key factor realizing rain bullet self-destroying function reliability;
(2) reliability model, set up according to step one clearly needs the priming system component reliability data of collecting.
5. a kind of analysis for rain bullet self-destroying function reliability according to claim 1 and appraisal procedure, is characterized in that: in step 3, and constructed mathematical model is mainly based on following hypothesis:
(1), in self-destruction unit whether priming system transmittance process at different levels is reliable, and whether self-destruction body brisance meets the demands, and namely body can be scraped completely when assigned address detonates and need be determined by the reliability of priming system element;
(2), self-destruction unit motion to the process of assigned address meets the relevant hypothesis of classical interior ballistics model;
(3), due to the small aperture being communicated with two coyote holes minimum, its circulated gases minimum relative to coyote hole gas gross thus hypothesis aperture in flow without gas;
(4), by the process that starting mix is pushed to assigned address, self-destruction unit do not considered that the viscosity of gas, gas are to the friction of body inwall, and combustion gas adopts indirect method correction to the heat transfer of chamber wall;
(5), powder burning and two self-destruction unit motions all carry out under mean pressure;
(6) in the middle part of, when considering that middle card blocks middle part self-destruction unit, self-destruction unit speed is minimum, also minimum to ballistic impact, therefore approximate think the heat that produces in this process to the second coyote hole without impact;
(7), due to starting mix burning and self-destruction unit motion process time very short, much smaller than the defer time of Delay Tube, thus its self-destruction time and height only relevant with the outer trajectory that defer time and the rain bullet of Delay Tube are launched, ignoring it at this affects;
Constructed mathematical model mainly comprises following six stages, specific as follows:
1), the first stage, to take fire afterbody self-destruction unit setting in motion from the gunpowder of two coyote holes (i.e. the first coyote hole, the second coyote hole), this process is coyote hole constant volume combustion process; The involved formula of whole process is as follows:
combinatorial formula (5)
This is a composite class formula, and wherein, first formula is the second coyote hole combustion gas growth equation formula, and second formula is the second coyote hole burning rate equation, and the 3rd formula is the second coyote hole pressure equation formula;
In above formula, p 2be the pressure of the second coyote hole, v 0' be the original volume of the second coyote hole; ω 2it is the explosive load of the second coyote hole;
With starting pressure P 0for initial value, solve the second coyote hole powder burning situation;
2), subordinate phase, from afterbody self-destruction unit and the time of releasing of middle part self-destruction unit to middle part self-destruction unit setting in motion, this process is the second coyote hole burning and external expansion work, the first coyote hole constant volume combustion process; The involved formula of whole process is as follows:
combinatorial formula (6)
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, second formula is the second coyote hole combustion gas growth equation formula, 3rd formula is the first coyote hole burning rate equation, 4th formula is the second coyote hole combustion gas growth equation formula, 5th formula is the first coyote hole pressure equation formula, 6th formula is the second coyote hole gas equation, 7th formula is afterbody self-destruction unit momentum equation, 8th formula is afterbody self-destruction unit motion equation, 9th formula is the second coyote hole energy conservation equation formula, tenth formula is for calculating end condition equation, 11 formula is the calculation equation of the enthalpy that the second coyote hole unit mass gas flows out,
In above formula, p 1be the pressure of the first coyote hole, v 0it is the original volume of the first coyote hole; ω 1be the explosive load of the first coyote hole, w 2be the gaseous mass percentage that the second coyote hole is revealed due to the gap between afterbody self-destruction unit and body, t 2be the gas temperature in the second coyote hole, l 2for the displacement of afterbody self-destruction unit relative middle part self-destruction unit, φ 2for the secondary merit design factor of afterbody self-destruction unit, m 2for the quality of afterbody self-destruction unit, v 2for the speed of afterbody self-destruction unit, s 2for the sectional area of afterbody self-destruction unit, p 0for external atmospheric pressure, r g =R/M, M are the molal weight of starting mix burning gases;
3), the phase III, from middle part self-destruction unit setting in motion to the second coyote hole gunpowder burns completely, this stage, two coyote hole gunpowder all burnt and expansion work; The involved formula of whole process is as follows:
combinatorial formula (7)
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, second formula is the first coyote hole burning rate equation, 3rd formula is the second coyote hole combustion gas growth equation formula, 4th formula is the first coyote hole gas equation, 5th formula is the second coyote hole gas equation, 6th formula is the momentum equation of middle part self-destruction unit, 7th formula is middle part self-destruction unit motion equation, 8th formula is the momentum equation of afterbody self-destruction unit, 9th formula is afterbody self-destruction unit motion equation, tenth formula is the energy conservation equation formula of the first coyote hole, 11 formula is the energy conservation equation formula of the second coyote hole, 12 formula is the calculation equation of the enthalpy that the first coyote hole unit mass gas flows out,
In above formula, w 1be the gaseous mass percentage that the first coyote hole is revealed due to the gap between middle part self-destruction unit and body, t 1be the gas temperature in the first coyote hole, l 1for the displacement of the relative body of middle part self-destruction unit, φ 1for the secondary merit design factor of middle part self-destruction unit, m 1for the quality of middle part self-destruction unit, v 1for the speed of middle part self-destruction unit, s 1for the sectional area of middle part self-destruction unit, h 1it is the enthalpy that the first coyote hole unit mass gas flows out;
4), fourth stage, be burned to till the first coyote hole starting mix burns completely completely from the second coyote hole starting mix, this stage second coyote hole expansion work, the first coyote hole burning and expansion work; The involved formula of whole process is as follows:
combinatorial formula (8)
In this combinatorial formula, first formula is the first coyote hole burning rate equation, second formula is the first coyote hole gas equation, 3rd formula is the second coyote hole gas equation, 4th formula is the energy conservation equation formula of the first coyote hole, 5th formula is the energy conservation equation formula of the second coyote hole, 6th formula is the momentum equation of middle part self-destruction unit, 7th formula is middle part self-destruction unit motion equation, 8th formula is the momentum equation of afterbody self-destruction unit, 9th formula is afterbody self-destruction unit motion equation,
5), five-stage, burn to assigned address, now two coyote hole expansion work in the middle part of middle part self-destruction unit motion to body completely from the first coyote hole starting mix; The involved formula of whole process is as follows:
combinatorial formula (9)
In this combinatorial formula, first formula is the first coyote hole gas equation, second formula is the second coyote hole gas equation, 3rd formula is the energy conservation equation formula of the first coyote hole, 4th formula is the energy conservation equation formula of the second coyote hole, 5th formula is the momentum equation of middle part self-destruction unit, 6th formula is middle part self-destruction unit motion equation, 7th formula is the momentum equation of afterbody self-destruction unit, and the 8th formula is afterbody self-destruction unit motion equation; The motion conditions that can calculate middle part self-destruction unit is solved by system of equations;
6), the 6th stage from assigned address in the middle part of middle part self-destruction unit motion to body to afterbody self-destruction unit motion to afterbody assigned address moment, in the middle part of this stage, self-destruction unit is no longer transported to because the associated mechanisms of body medium position blocks, first coyote hole gas is no longer revealed, the still motion under gaseous-pressure promotes of afterbody self-destruction unit; The involved formula of whole process is as follows:
combinatorial formula (10)
In this combinatorial formula, first formula is the second coyote hole gas equation, second formula is the equation of momentum of afterbody self-destruction unit, 3rd formula is afterbody self-destruction unit motion equation, 4th formula is the energy conservation equation formula of the second coyote hole, is gone out the motion conditions of afterbody self-destruction unit by solving equations;
In above-mentioned six system of equations:
(1), the first coyote hole flows into the gas-flow equation of the second coyote hole:
formula (11)
In formula , φ 3for Flux Loss coefficient, s 1' be the sectional area in middle part self-destruction unit and body gap, m 1' be the gaseous mass of the first coyote hole;
(2), the second coyote hole flows into extraneous gas-flow equation:
formula (12)
In formula , φ 4for Flux Loss coefficient, s 2' be the sectional area in afterbody self-destruction unit and body gap, m 2' be the gaseous mass of the second coyote hole.
6. a kind of analysis for rain bullet self-destroying function reliability according to claim 1 and appraisal procedure, it is characterized in that: in step 4, after the principal element considering impact middle part and afterbody self-destruction unit motion, utilize monte carlo method to middle part self-destruction unit motion to middle card position, afterbody self-destruction unit motion is analyzed to the reliability of nozzle hole position, and concrete steps are as follows:
A, will realize self-destroying function, middle part self-destruction unit can be pushed into card position in body by demand fulfillment starting mix, afterbody self-destruction unit is pushed into nozzle hole position simultaneously, if in the middle part of body under assigned address to original state the distance of middle part self-destruction unit bottom be y 1, afterbody assigned address to the distance of original state lower tail self-destruction unit bottom is y 2;
In actual production process, due to manufacturing accuracy impact, y 1with y 2for obeying the parameter of certain regularity of distribution, and totally y 1distribution function be f 1 y ( y), totally y 2distribution function be f 2 y ( y);
B, middle part self-destruction unit relative body displacement can be drawn according to above-mentioned mathematical model by emulation l 1a sample and the relative body displacement of afterbody self-destruction unit l 2a sample, and l 1distribution function be f 1 z ( z), totally l 2distribution function be f 2 z ( z); Middle part, afterbody self-destruction unit can arrive in the middle part of body, afterbody assigned address need meet l 1>= y 1, and l 2>= y 2this condition;
C, determine parameters distribution function after, then can determine the fiduciary level of assigned address in the middle part of middle part self-destruction unit motion to body, computing formula is:
formula (13)
The fiduciary level that afterbody self-destruction unit motion is specified to afterbody, computing formula is:
formula (14)
Affect the powder charge amount that self-destruction unit motion to the topmost factor of assigned address is the first coyote hole and the second coyote hole ω 1, ω 2and the quality of middle part self-destruction unit and afterbody self-destruction unit m 1, m 2;
In a practical situation, explosive load ω i ( i=1,2) and each self-destruction element quality m j ( j=1,2) all there is certain span and statistical nature: ω ip -△ ω i ≤ ω i ≤ ω ip + △ ω i formula (15)
m jp -△ m j ≤ m j ≤ m jp + △ m j formula (16)
In formula ω ip for ω i nominal value, △ ω i for ω i tolerance, m jp for m j nominal value, △ m j for m j tolerance;
D, according to set up the mathematical model based on interior ballistics theory, the function that the displacement of self-destruction unit relative body is relevant to two the chamber charges and two self-destruction element qualities can be obtained:
l 1= d( ω i , m j ) formula (17)
l 2= d( ω i , m j ) formula (18)
Formula (17) is with the displacement of the relative body of middle part self-destruction unit l 1for the system of equations of unknown function, formula (18) is the displacement of the relative body of afterbody self-destruction unit l 2for the system of equations of unknown function;
E, for obtaining l 1, l 2distribution function, emulated by Monte Carlo, carry out random sampling according to the requirement that the quality of actual conditions to the explosive load of two coyote holes and middle part self-destruction unit and afterbody self-destruction unit is checked on by its regularity of distribution of obeying and quality inspection; As kduring secondary sampling, extract on request ω ik ( i=1,2) and m jk ( j=1,2); Will ω ik with m jk substitute into (17), (18) formula (i.e. solving equation group) can obtain respectively in the middle part of the operation sample value of self-destruction unit and the displacement of the relative body of afterbody self-destruction unit:
l 1 k = d( ω 1 k , ω 2 k ; m 1 k , m 2 k ) formula (19)
l 2 k = d( ω 1 k , ω 2 k ; m 1 k , m 2 k ) formula (20)
After realistic model runs N time, the shift value of the relative body of N number of middle part self-destruction unit can be obtained l 1 k ( k=1,2 ..., N) and the shift value of the relative body of N number of afterbody self-destruction unit l 2 k (k=1,2 ..., N); Respectively with N number of l 1 k value is as the relative body displacement of middle part self-destruction unit l 1a sample, with N number of l 2 k value is as the relative body displacement of afterbody self-destruction unit l 2a sample, draw the frequency histogram of two parameters, and guestimate or suppose overall l 1distribution function f 1 z ( z), totally l 2distribution function f 2 z ( z); By sample data, hypothesis is tested, to judge the correctness supposed; Finally determine distribution function f 1 z ( z) and distribution function f 2 z ( z), can determine in the middle part of rain bullet by (13) (14) formula and afterbody self-destruction unit motion to the fiduciary level of assigned address.
7. a kind of analysis for rain bullet self-destroying function reliability according to claim 6 and appraisal procedure, is characterized in that: in step e, the process that Monte Carlo emulates, and concrete steps are as follows:
Step (one), clear and definite sampling parametric ω 1, ω 2with m 1, m 2and simulation times N;
Step (two), Monte Carlo is utilized to carry out random sampling to parameter;
Step (three), by parameter ω 1, ω 2with m 1, m 2and simulation times N substitutes into formula (17), (18) solve;
Step (four), output l 1 k , l 2k and stored in array l 1, l 2; N in this process>= k;
Step (five), output array l 1, l 2, and draw its distribution function.
CN201410335207.4A 2014-07-15 2014-07-15 A kind of analysis and appraisal procedure for rain bullet self-destroying function reliability Expired - Fee Related CN104298848B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410335207.4A CN104298848B (en) 2014-07-15 2014-07-15 A kind of analysis and appraisal procedure for rain bullet self-destroying function reliability

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410335207.4A CN104298848B (en) 2014-07-15 2014-07-15 A kind of analysis and appraisal procedure for rain bullet self-destroying function reliability

Publications (2)

Publication Number Publication Date
CN104298848A true CN104298848A (en) 2015-01-21
CN104298848B CN104298848B (en) 2017-07-28

Family

ID=52318571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410335207.4A Expired - Fee Related CN104298848B (en) 2014-07-15 2014-07-15 A kind of analysis and appraisal procedure for rain bullet self-destroying function reliability

Country Status (1)

Country Link
CN (1) CN104298848B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107742007A (en) * 2017-09-18 2018-02-27 中国人民解放军海军工程大学 The computational methods of sheet metal ballisticslimited velocity under a kind of positive penetration of flat nose low speed
CN107742006A (en) * 2017-09-18 2018-02-27 中国人民解放军海军工程大学 The computational methods of sheet metal ballisticslimited velocity under a kind of positive penetration of tack hollow projectile cartridge low speed

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
卢培玉 等: "BL-1A增雨防雹火箭系统安全使用探讨", 《山东气象》 *
王铮: "增雨防雹火箭的可靠性分配及其单元可靠性增长试验", 《固体火箭技术》 *
苏春 等: "基于蒙特卡洛仿真的液压系统动态可靠性", 《东南大学学报》 *
魏旭辉 等: "增雨防雹火箭弹的工作原理和技术参数", 《沙漠与绿洲气象》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107742007A (en) * 2017-09-18 2018-02-27 中国人民解放军海军工程大学 The computational methods of sheet metal ballisticslimited velocity under a kind of positive penetration of flat nose low speed
CN107742006A (en) * 2017-09-18 2018-02-27 中国人民解放军海军工程大学 The computational methods of sheet metal ballisticslimited velocity under a kind of positive penetration of tack hollow projectile cartridge low speed
CN107742007B (en) * 2017-09-18 2021-05-18 中国人民解放军海军工程大学 Method for calculating limit speed of sheet steel trajectory under low-speed penetration of flush bomb
CN107742006B (en) * 2017-09-18 2021-05-18 中国人民解放军海军工程大学 Method for calculating limit speed of sheet steel trajectory under low-speed forward penetration of flat-head hollow bullet

Also Published As

Publication number Publication date
CN104298848B (en) 2017-07-28

Similar Documents

Publication Publication Date Title
CN105446167B (en) Hypersonic scramjet engine real-time model, emulation mode
CN104298848A (en) Method for analyzing and evaluating reliability of self-destroying function of cannonball for artificial precipitation
CN110287599A (en) A kind of propellant powder generates the construction method of plasma rule model
CN106599352B (en) Reliability analysis method for airplane telex control system
Cai et al. Adaptive simulations of cavity-based detonation in supersonic hydrogen–oxygen mixture
Wessam et al. Firing precision evaluation for unguided artillery projectile
Brophy et al. Detonation initiation improvements using swept-ramp obstacles
Paxson A simplified model for detonation based pressure-gain combustors
McNesby et al. Optical measurement of peak air shock pressures following explosions
US6519554B1 (en) Computer implemented system and method for evaluating gas generator launchers
Derk et al. High-pressure burning rate measurements by direct observation
Omang et al. Height of burst explosions: a comparative study of numerical and experimental results
CN109738196B (en) Ramjet performance space-ground conversion method based on ground direct connection test
Cao et al. A novel launching system applying a relay chamber technology and its optimization
M. Rashad et al. Two-phase flow interior ballistics model of naval large caliber guided projectile gun system
CN202900202U (en) Low-minute level twin-core cord type delay detonating mechanism
Wilhelm et al. Flight simulation of a hybrid projectile to estimate the impact of launch angle on range extension
Cao Parametric cycle analysis of continuous rotating detonation ejector-augmented rocket engine
Abou-Elela et al. Ballistic analysis of a projectile provided with base bleed unit
Tao et al. Output reliability assessment of nuclear initiating explosive device based on bootstrap method
Serre et al. The French PROMETHEE program on hydrocarbon fueled dual mode ramjet-Status in 2001
CN203584388U (en) Simulated perforating bullet and composite perforator
CN109342067A (en) A kind of solid rocket ramjet test macro
Viguier et al. Test of the Lapcat II small scale flight experiment model in the ONERA F4 wind tunnel
Acharya et al. Comprehensive Three Dimensional Mortar Interior Ballistics Model for 120mm Mortar System with Experimental Validation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170728

Termination date: 20210715