CN104298848B - A kind of analysis and appraisal procedure for rain bullet self-destroying function reliability - Google Patents
A kind of analysis and appraisal procedure for rain bullet self-destroying function reliability Download PDFInfo
- Publication number
- CN104298848B CN104298848B CN201410335207.4A CN201410335207A CN104298848B CN 104298848 B CN104298848 B CN 104298848B CN 201410335207 A CN201410335207 A CN 201410335207A CN 104298848 B CN104298848 B CN 104298848B
- Authority
- CN
- China
- Prior art keywords
- self
- formula
- destruction unit
- destruction
- coyote hole
- 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.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of analysis for rain bullet self-destroying function reliability and appraisal procedure, this method is simple to operate, can fast and effeciently realize the analysis and assessment of rain bullet self-destroying function reliability；Task pattern of the invention according to the type self-destroyed rain bullet self-destroying function of gunpowder, propose the rain bullet system self-destruction fail-safe analysis assessment models that self-destruction is realized by fixed gunpowder self-destruction unit and motion gunpowder self-destruction unit, the simulation model of motion gunpowder self-destruction unit reliability is constructed with classical interior ballistics principle, and then the reliability of rain bullet self-destroying function is designed, analyzed and assessed.
Description
Technical field
The present invention relates to a kind of analysis of reliability and appraisal procedure, more particularly to be exclusively used in gunpowder type self-destroyed artificial for one kind
The analysis of the reliability of anti-hail bullet self-destroying function and appraisal procedure.
Background technology
Self-destruction is hail-suppression rocket bullet of increasing rain artificially（Abbreviation rain bullet）One of the mode of body safe landing, self-destruction is typically
Completed by self-destruction subsystem responsible；Traditional self-destruction subsystem includes head self-destruction unit, middle part self-destruction unit and afterbody self-destruction
Unit；Self-destruction unit in head is fixed on the front end of body in advance；And middle part self-destruction unit and afterbody self-destruction unit are then in self-destruction
Specified location is moved in journey；But can middle part self-destruction unit and afterbody self-destruction unit reliably be transported to specified location,
Can body is then directly connected to fully achieve self-destruction；Therefore, to middle part self-destruction unit and the motion process of afterbody self-destruction unit
Fail-safe analysis is carried out, is the key of correctly predicted rain bullet self-destroying function reliability；
At present, the country is not yet studied rain bullet self-destroying function reliability design and analysis and appraisal procedure, therefore also
There is no definite method to be determined rain bullet self-destroying function reliability；Therefore, we are according to passing through fixed gunpowder self-destruction unit
The task pattern of the type self-destroyed rain bullet self-destroying function of self-destruction is realized with motion gunpowder self-destruction unit, it is proposed that gunpowder is type self-destroyed
Rain bullet fail-safe analysis appraisal procedure, to design a kind of analysis of self-destruction rain bullet self-destroying function reliability, assessment models, by this
Model carries out design and analysis and the assessment of rain bullet self-destroying function reliability.
The content of the invention
The present invention is in order to make up the deficiencies in the prior art there is provided a kind of analysis for rain bullet self-destroying function reliability with commenting
Estimate method, this method is simple to operate, can fast and effeciently realize the analysis of the type self-destroyed rain bullet self-destroying function reliability of gunpowder with commenting
Estimate.
For achieving the above object, the present invention is adopted the following technical scheme that：
A kind of analysis and appraisal procedure for rain bullet self-destroying function reliability, it is characterised in that：Including walking in detail below
Suddenly：
Step 1: building the logic diagram model of rain bullet self-destroying function subsystem reliability according to operation principle；
Step 2: analyzing influence realizes the key factor of rain bullet self-destroying function reliability, and collect priming system component
Reliability data；
Step 3: self-destruction unit and afterbody self-destruction unit need starting mix gaseous-pressure to promote motion in the middle part of rain bullet
To the operation principle of specified location, based on classical interior ballistics theory, the mathematical modeling of the rain bullet self-destruction course of work is set up；
Step 4: utilizing work of the monte carlo method to middle part self-destruction unit and afterbody self-destruction unit motion to specified location
Make process and carry out analog simulation, solve middle part self-destruction unit and afterbody self-destruction unit motion to the reliable number of degrees of specified location
According to；
Can to increasing rain artificially rail-proof bomb self-destroying function Step 5: combining the rain bullet self-destroying function reliability block diagram model built
Analyzed and assessed by property；
Further, in step one, the reliability block diagram model includes gunpowder extension unitR 1, head self-destruction unitR T, middle part self-destruction unitR C, afterbody self-destruction unitR A；The gunpowder extension unitR 1With head self-destruction unitR TInput end signal
Connection；The head self-destruction unitR TOutput end and middle part self-destruction unitR CInput end signal connection；The middle part self-destruction list
MemberR COutput end and afterbody self-destruction unitR AInput end signal connection；
Wherein, gunpowder extension unitR 1Including electric igniterR 11, electric igniterR 12, electric igniterR 17, electric igniterR 18、
Flame agent Dark Ignition hasR 13, flame agent Dark Ignition toolR 14, flame agentR 15, one-level Delay TubeR 16, electric Delay TubeR 19, starting mixR 110；
The electric igniterR 11Have with flame agent Dark IgnitionR 13Signal is connected；The electric igniterR 12Have with flame agent Dark IgnitionR 14Letter
Number connection；The flame agent Dark Ignition toolR 13, flame agent Dark Ignition toolR 14With flame agentR 15Signal is connected；The flame agentR 15With
One-level Delay TubeR 16Signal is connected；The one-level Delay TubeR 16, electric Delay TubeR 19And starting mixR 110Signal is connected；The electricity
Igniter headR 11Respectively with electric igniterR 17, electric igniterR 18Signal is connected；
According to the reliability block diagram, set the Reliable Mathematics model of rain bullet self-destroying function as：
R=R 1 R T R C R A Formula（1）
The reliability of gunpowder extension unit, its formula of reliability are can obtain by self-destruction subsystem reliability logic diagram
For：
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)
WhereinR 11、R 12、R 17、R 18The reliability of electric igniter is expressed as,R 13、R 14Be expressed as flame agent Dark Ignition tool can
By degree,R 15The reliability of flame agent is expressed as,R 16The reliability of one-level Delay Tube is expressed as,R 19It is expressed as the reliable of electric Delay Tube
Degree,R 110It is expressed as the reliability of starting mix；The reliability of middle part self-destruction unitR CAnd the reliability of afterbody self-destruction unitR AIt is then
The reliability of the two unit priming systems is moved to the synthesis of the reliability of designated displacement, i.e.,
R C= R C1× R CL Formula (3)
R A= R A1× R AL Formula (4)
Wherein,R C1For the reliability of middle part self-destruction unit priming system,R A1For the reliability of afterbody self-destruction unit priming system.
Further, in step one, the operation principle realized according to rain bullet self-destroying function enters reliability model structure
Row simplifies, and the model includes body, the first coyote hole, the second coyote hole, middle part self-destruction body, afterbody self-destruction body, middle card and nozzle hole；
First coyote hole is arranged on body foremost, and is set close to middle part self-destruction body；Second coyote hole is arranged on middle part self-destruction body
With the middle part of afterbody self-destruction body；The middle card is arranged on the middle part of body；The nozzle hole is arranged on the tail end of body；
Further, in step 2, the key factor of rain bullet self-destroying function reliability, which is analyzed, to be realized to influence, it is main
To include the following aspects：
（1）, the reliability block diagram set up by step one determine influence realize rain bullet self-destroying function reliability it is crucial because
Element is reliability of two self-destruction unit motions to specified location；
（2）, the reliability model set up according to step one clearly require the priming system component reliability data of collection；
Further, in step 3, constructed mathematical modeling is based primarily upon hypothesis below：
（1）, whether priming system transmittance process at different levels reliable in self-destruction unit, whether self-destruction body brisance meets requirement,
Body can be scraped completely i.e. when specified location is detonated to be determined by the reliability of priming system element；
（2）, self-destruction unit motion to specified location process meet the related hypothesis of classical interior ballistics model；
（3）, due to connect two coyote holes small aperture it is minimum, its circulated gases is minimum with respect to for coyote hole gas gross
Therefore assume in aperture without gas flowing；
（4）, self-destruction unit is by the process that starting mix is pushed to specified location does not consider the viscosity of gas, gas is in body
The friction of wall, combustion gas uses indirect method amendment to the heat transfer of chamber wall；
（5）, powder burning and two self-destruction unit motions carry out all under average pressure；
（6）, in view of it is middle block middle part self-destruction unit when middle part self-destruction unit speed it is minimum, to ballistic impact also pole
It is small, therefore be approximately considered in the heat produced by this process on the second coyote hole without influence；
（7）, due to starting mix burning and self-destruction unit motion process time it is very short, much smaller than the extension time of Delay Tube,
Therefore its self-destruction time and height are only relevant with the outer trajectory that delay time and the rain bullet of Delay Tube are launched, and its influence is ignored herein.
Constructed mathematical modeling mainly includes the following six stage, specific as follows：
1st, the first stage, from two coyote holes（That is the first coyote hole, the second coyote hole）Gunpowder start to be burned to afterbody self-destruction list
First setting in motion, this process is coyote hole constant volume burning process；The involved formula of whole process is as follows：
Combinatorial formula（5）
This is a composite class formula, wherein, first formula is the second coyote hole combustion gas growth equation formula, second formula
For the second coyote hole burning velocity equation, the 3rd formula is the second coyote hole pressure equation formula；
In above formula,p 2For the pressure of the second coyote hole,V 0' be the second coyote hole original volume；ω 2For the second coyote hole
Explosive payload.With starting pressure P0For initial value, the second coyote hole powder burning situation is solved；
2nd, second stage, is transported since afterbody self-destruction unit and the time of releasing of middle part self-destruction unit are to middle part self-destruction unit
Dynamic, this process is the burning of the second coyote hole and external expansion work, the first coyote hole constant volume burning process；Whole process it is involved
Formula is as follows：
Combinatorial formula (6)
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, and second formula is second
Coyote hole combustion gas growth equation formula, the 3rd formula is the first coyote hole burning velocity equation, and the 4th formula fires for the second coyote hole
Gas growth equation formula, the 5th formula is the first coyote hole pressure equation formula, and the 6th formula is the second coyote hole equation of gas state
Formula, the 7th formula is afterbody self-destruction unit momentum equation, and the 8th formula is afterbody self-destruction unit motion equation, the 9th
Individual formula is the second coyote hole energy conservation equation formula, and the tenth formula is calculates end condition equation, and the 11st formula is
The calculation equation of the enthalpy of second coyote hole unit mass gas outflow；
In above formula,p 1For the pressure of the first coyote hole,V 0For the original volume of the first coyote hole；ω 1For the first coyote hole
Explosive payload,w 2The gaseous mass percentage revealed for the second coyote hole due to the gap between afterbody self-destruction unit and body,T 2For
Gas temperature in two coyote holes,l 2Displacement for afterbody self-destruction unit with respect to middle part self-destruction unit,φ 2For afterbody self-destruction unit
Secondary work(design factor,m 2For the quality of afterbody self-destruction unit,v 2For the speed of afterbody self-destruction unit,S 2For afterbody self-destruction unit
Sectional area,p 0For external atmospheric pressure,R g =R/M, M are the molal weight of starting mix burning gases；
3rd, the phase III, from middle part self-destruction unit setting in motion untill the second coyote hole gunpowder completely burning, this stage two
Individual coyote hole gunpowder burns 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, and second formula is first
Coyote hole burning velocity equation, the 3rd formula is the second coyote hole combustion gas growth equation formula, and the 4th formula is the first coyote hole gas
Body equation of state, the 5th formula is the second coyote hole gas equation, and the 6th formula is the dynamic of middle part self-destruction unit
Equation is measured, the 7th formula is middle part self-destruction unit motion equation, and the 8th formula is the momentum side of afterbody self-destruction unit
Formula, the 9th formula is afterbody self-destruction unit motion equation, and the tenth formula is the energy conservation equation formula of the first coyote hole,
11st formula is the energy conservation equation formula of the second coyote hole, and the 12nd formula flows out for the first coyote hole unit mass gas
Enthalpy calculation equation；
In above formula,w 1For gas matter of first coyote hole due to the gap leakage between middle part self-destruction unit and body
Measure percentage,T 1For the gas temperature in the first coyote hole,l 1Displacement for middle part self-destruction unit with respect to body,φ 1For middle part self-destruction
The secondary work(design factor of unit,m 1For the quality of middle part self-destruction unit,v 1For the speed of middle part self-destruction unit,S 1For middle part certainly
The sectional area of unit is ruined,h 1The enthalpy flowed out for the first coyote hole unit mass gas；
4th, fourth stage, untill being burned to the burning completely of the first coyote hole starting mix completely from the second coyote hole starting mix, this rank
The second coyote hole expansion work of section, the burning of the first coyote hole 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 velocity equation, and second formula is first
Coyote hole gas equation, the 3rd formula is the second coyote hole gas equation, and the 4th formula is the first coyote hole
Energy conservation equation formula, the 5th formula is the energy conservation equation formula of the second coyote hole, and the 6th formula is middle part self-destruction unit
Momentum equation, the 7th formula is middle part self-destruction unit motion equation, and the 8th formula is the dynamic of afterbody self-destruction unit
Equation is measured, the 9th formula is afterbody self-destruction unit motion equation；
5th, the 5th stage, from the burning completely of the first coyote hole starting mix to middle part self-destruction unit motion to body middle part specific bit
Put, now two coyote hole expansion works；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, and second formula is second
Coyote hole gas equation, the 3rd formula is the energy conservation equation formula of the first coyote hole, and the 4th formula is the second coyote hole
Energy conservation equation formula, the 5th formula is the momentum equation of middle part self-destruction unit, and the 6th formula is middle part self-destruction list
First equation of motion, the 7th formula is the momentum equation of afterbody self-destruction unit, and the 8th formula is transported for afterbody self-destruction unit
Dynamic equation；The motion conditions of middle part self-destruction unit can be calculated by the solution of equation group；
6th, the 6th stage from middle part self-destruction unit motion to body middle part specified location to afterbody self-destruction unit motion to bullet
Self-destruction unit is no longer transported because the associated mechanisms of body medium position are blocked in the middle part of body afterbody specified location moment, this stage
Arrive, the first coyote hole gas is no longer revealed, afterbody self-destruction unit is still moved under gaseous-pressure promotion；Whole process it is involved
Formula it is as follows：
Combinatorial formula (10)
In this combinatorial formula, first formula is the second coyote hole gas equation, and second formula is afterbody
The equation of momentum of self-destruction unit, the 3rd formula is afterbody self-destruction unit motion equation, and the 4th formula is the energy of the second coyote hole
Conservation equation formula is measured, the motion conditions of afterbody self-destruction unit are gone out by solving equations；
In above-mentioned six equation groups：
（1）, the first coyote hole flow into the second coyote hole gas-flow equation：
Formula（11）
In formula,φ 3For Flux Loss coefficient,S 1' it is middle part self-destruction unit and the section in body gap
Product,m 1' be the first coyote hole gaseous mass；
（2）, the second coyote hole flow into extraneous gas-flow equation：
Formula（12）
In formula,φ 4For Flux Loss coefficient,S 2' it is afterbody self-destruction unit and the section in body gap
Product,m 2' be the second coyote hole gaseous mass；
Further, in step 4, considering in the middle part of influence and after the principal element of afterbody self-destruction unit motion, to utilize
Monte carlo method is to middle part self-destruction unit motion to middle card position, the reliability of afterbody self-destruction unit motion to nozzle hole position
Analyzed, comprised the following steps that：
A, self-destroying function is realized, it is necessary to which card position in body can be pushed into by middle part self-destruction unit by meeting starting mix, simultaneously
Afterbody self-destruction unit is pushed into nozzle hole position, if specified location is to middle part self-destruction unit bottom under original state in the middle part of body
Distance isY 1, the distance of afterbody specified location to original state lower tail self-destruction unit bottom isY 2；
In actual production process, due to accuracy of manufacture influence,Y 1WithY 2To obey the parameter of certain regularity of distribution, and it is overallY 1Distribution function bef 1y (y), it is overallY 2Distribution function bef 2y (y)；
B, according to above-mentioned mathematical modeling by emulation middle part self-destruction unit can be drawn with respect to body displacementL 1A sample
And afterbody self-destruction unit is with respect to body displacementL 2A sample, andL 1Distribution function bef 1Z (z), it is overallL 2Distribution function
Forf 2Z (z)；Middle part, afterbody self-destruction unit can be reached in the middle part of body, afterbody specified location needs to meetL 1≥Y 1, andL 2≥Y 2This
One condition；
C, determine after the distribution functions of parameters, then can determine that specified location in the middle part of middle part self-destruction unit motion to body
Reliability, calculation formula is：
Formula（13）
The reliability that afterbody self-destruction unit motion to afterbody is specified, calculation formula is：
Formula（14）
Self-destruction unit motion to the topmost factor of specified location is influenceed to be filled for the gunpowder of the first coyote hole and the second coyote hole
Doseω 1、ω 2And the quality of middle part self-destruction unit and afterbody self-destruction unitm 1、m 2；
In a practical situation, explosive payloadω i (i=1,2) and each self-destruction element qualitym j (j=1,2) there is certain value
Scope 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 Form j Nominal value, △m j Form j Tolerance；
D, the mathematical modeling based on interior ballistics theory according to foundation, can obtain self-destruction unit with respect to body displacement and two
The chamber charge and the related function of two self-destruction element qualities：
L 1=D(ω i ,m j ) formula（17）
L 2=D(ω i ,m j ) formula（18）
Formula（17）It is the displacement with middle part self-destruction unit with respect to bodyL 1For the equation group of unknown function, formula（18）It is afterbody
Displacement of the self-destruction unit with respect to bodyL 2For the equation group of unknown function；
E, to obtainL 1、L 2Distribution function, emulated by Monte Carlo, according to dress of the actual conditions to two coyote holes
The requirement that the regularity of distribution and quality inspection that the quality of dose and middle part self-destruction unit and afterbody self-destruction unit is obeyed by it are checked on is carried out
Random sampling；Such askDuring secondary sampling, extract on requestω ik (i=1,2) andm jk (j=1,2)；Willω ik Withm jk Substitution (17),
（18）Formula（Solve equation group）The operation of the displacement of middle part self-destruction unit body relative with afterbody self-destruction unit can be respectively obtained
Sample value：
l 1k =D(ω 1k ,ω 2k ;m 1k , m 2k ) formula（19）
l 2k =D(ω 1k ,ω 2k ;m 1k , m 2k ) formula（20）
After simulation model runs n times, shift value of N number of middle part self-destruction unit with respect to body can obtainl 1k (k=1,2,…,
N) and N number of afterbody self-destruction unit with respect to body shift valuel 2k (k=1,2,…,N)；Respectively with N number ofl 1k Value is used as middle part self-destruction
Unit is with respect to body displacementL 1A sample, with N number ofl 2k Value is as afterbody self-destruction unit with respect to body displacementL 2A sample
This, draws the frequency histogram of two parameters, and rough estimate or assumes totalityL 1Distribution functionf 1Z (z), it is overallL 2's
Distribution functionf 2Z (z)；Assume to test by sample data pair, to judge the correctness assumed；It is final to determine distribution functionf 1Z (z) and distribution functionf 2Z (z), by（13）（14）Formula can determine that rain bullet middle part and afterbody self-destruction unit motion to specified location
Reliability；
Further, in step e, the process of Monte Carlo emulation is comprised the following steps that：
Step（One）, clear and definite sampling parametricω 1、ω 2Withm 1、m 2And simulation times N；
Step（Two）, using Monte Carlo to parameter carry out random sampling；
Step（Three）, by parameterω 1、ω 2Withm 1、m 2And simulation times N substitutes into formula（17）、（18）Solved；
Step（Four）, outputl 1k 、l 2k And it is stored in arrayL 1、L 2；N during this >=k；
Step（Five）, output arrayL 1、L 2, and draw distribution function；
Compared with prior art, using such scheme, the beneficial effects of the invention are as follows：The present invention is according to the type self-destroyed rain of gunpowder
Play the task pattern of self-destroying function, it is proposed that self-destruction is realized by fixed gunpowder self-destruction unit and motion gunpowder self-destruction unit
Rain bullet system self-destruction fail-safe analysis assessment models, with classical interior ballistics principle construct motion gunpowder self-destruction unit can
It is designed, analyzes and assesses by the simulation model of property, and then to the reliability of rain bullet self-destroying function.
Brief description of the drawings
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 invention；
Fig. 4 is Monte Carlo simulation flow schematic diagrames in the present invention.
Embodiment
The present invention is further described with reference to the accompanying drawings and examples.
As shown in figure 1, a kind of analysis and appraisal procedure for rain bullet self-destroying function reliability, it is characterised in that：Including
Step in detail below：
Step 1: building the logic diagram model of rain bullet self-destroying function subsystem reliability according to operation principle；
Step 2: analyzing influence realizes the key factor of rain bullet self-destroying function reliability, and collect priming system component
Reliability data；
Step 3: self-destruction unit and afterbody self-destruction unit need starting mix gaseous-pressure to promote motion in the middle part of rain bullet
To the operation principle of specified location, based on classical interior ballistics theory, the mathematical modeling of the rain bullet self-destruction course of work is set up；
Step 4: utilizing work of the monte carlo method to middle part self-destruction unit and afterbody self-destruction unit motion to specified location
Make process and carry out analog simulation, solve middle part self-destruction unit and afterbody self-destruction unit motion to the reliable number of degrees of specified location
According to；
Can to increasing rain artificially rail-proof bomb self-destroying function Step 5: combining the rain bullet self-destroying function reliability block diagram model built
Analyzed and assessed by property；
Preferably, as shown in Fig. 2 in step one, the reliability block diagram model includes gunpowder extension unitR 1, head
Self-destruction unitR T, middle part self-destruction unitR C, afterbody self-destruction unitR A；The gunpowder extension unitR 1With head self-destruction unitR TIt is defeated
Enter end signal connection；The head self-destruction unitR TOutput end and middle part self-destruction unitR CInput end signal connection；In described
Portion's self-destruction unitR COutput end and afterbody self-destruction unitR AInput end signal connection；
Wherein, gunpowder extension unitR 1Including electric igniterR 11, electric igniterR 12, electric igniterR 17, electric igniterR 18、
Flame agent Dark Ignition hasR 13, flame agent Dark Ignition toolR 14, flame agentR 15, one-level Delay TubeR 16, electric Delay TubeR 19, starting mixR 110；
The electric igniterR 11Have with flame agent Dark IgnitionR 13Signal is connected；The electric igniterR 12Have with flame agent Dark IgnitionR 14Letter
Number connection；The flame agent Dark Ignition toolR 13, flame agent Dark Ignition toolR 14With flame agentR 15Signal is connected；The flame agentR 15With
One-level Delay TubeR 16Signal is connected；The one-level Delay TubeR 16, electric Delay TubeR 19And starting mixR 110Signal is connected；The electricity
Igniter headR 11Respectively with electric igniterR 17, electric igniterR 18Signal is connected；
According to the reliability block diagram, set the Reliable Mathematics model of rain bullet self-destroying function as：
R=R 1 R T R C R A Formula（1）
The reliability of gunpowder extension unit, its formula of reliability are can obtain by self-destruction subsystem reliability logic diagram
For：
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)
WhereinR 11、R 12、R 17、R 18The reliability of electric igniter is expressed as,R 13、R 14Be expressed as flame agent Dark Ignition tool can
By degree,R 15The reliability of flame agent is expressed as,R 16The reliability of one-level Delay Tube is expressed as,R 19It is expressed as the reliable of electric Delay Tube
Degree,R 110It is expressed as the reliability of starting mix；The reliability of middle part self-destruction unitR CAnd the reliability of afterbody self-destruction unitR AIt is then
The reliability of the two unit priming systems is moved to the synthesis of the reliability of designated displacement, i.e.,
R C= R C1× R CL Formula (3)
R A= R A1× R AL Formula (4)
Wherein,R C1For the reliability of middle part self-destruction unit priming system,R A1For the reliability of afterbody self-destruction unit priming system.
Preferably, as shown in figure 3, in step one, the operation principle realized according to rain bullet self-destroying function, by reliability mould
Type structure is simplified, the model include body 5, the first coyote hole 1, the second coyote hole 3, middle part self-destruction body 2, afterbody self-destruction body 4,
Middle card 6 and nozzle hole 7；First coyote hole 1 is arranged on body 5 foremost, and is set close to middle part self-destruction body 2；Second medicine
Room 3 is arranged on the middle part of middle part self-destruction body 2 and afterbody self-destruction body 4；The middle card 6 is arranged on the middle part of body 5；The spray
The mouth of pipe 7 is arranged on the tail end of body 5；
Preferably, in step 2, the key factor of rain bullet self-destroying function reliability, which is analyzed, to be realized to influence, mainly
Including the following aspects：
（1）, the reliability block diagram set up by step one determine influence realize rain bullet self-destroying function reliability it is crucial because
Element is reliability of two self-destruction unit motions to specified location；
（2）, the reliability model set up according to step one clearly require the priming system component reliability data of collection；
Preferably, in step 3, constructed mathematical modeling is based primarily upon hypothesis below：
（1）, whether priming system transmittance process at different levels reliable in self-destruction unit, whether self-destruction body brisance meets requirement,
Body can be scraped completely i.e. when specified location is detonated to be determined by the reliability of priming system element；
（2）, self-destruction unit motion to specified location process meet the related hypothesis of classical interior ballistics model；
（3）, due to connect two coyote holes small aperture it is minimum, its circulated gases is minimum with respect to for coyote hole gas gross
Therefore assume in aperture without gas flowing；
（4）, self-destruction unit is by the process that starting mix is pushed to specified location does not consider the viscosity of gas, gas is in body
The friction of wall, combustion gas uses indirect method amendment to the heat transfer of chamber wall；
（5）, powder burning and two self-destruction unit motions carry out all under average pressure；
（6）, in view of it is middle block middle part self-destruction unit when middle part self-destruction unit speed it is minimum, to ballistic impact also pole
It is small, therefore be approximately considered in the heat produced by this process on the second coyote hole without influence；
（7）, due to starting mix burning and self-destruction unit motion process time it is very short, much smaller than the extension time of Delay Tube,
Therefore its self-destruction time and height are only relevant with the outer trajectory that delay time and the rain bullet of Delay Tube are launched, and its influence is ignored herein.
Constructed mathematical modeling mainly includes the following six stage, specific as follows：
1st, the first stage, from two coyote holes（That is the first coyote hole, the second coyote hole）Gunpowder start to be burned to afterbody self-destruction list
First setting in motion, this process is coyote hole constant volume burning process；The involved formula of whole process is as follows：
Combinatorial formula（5）
This is a composite class formula, wherein, first formula is the second coyote hole combustion gas growth equation formula, second formula
For the second coyote hole burning velocity equation, the 3rd formula is the second coyote hole pressure equation formula；
In above formula,p 2For the pressure of the second coyote hole,V 0' be the second coyote hole original volume；ω 2For the second coyote hole
Explosive payload.With starting pressure P0For initial value, the second coyote hole powder burning situation is solved；
2nd, second stage, is transported since afterbody self-destruction unit and the time of releasing of middle part self-destruction unit are to middle part self-destruction unit
Dynamic, this process is the burning of the second coyote hole and external expansion work, the first coyote hole constant volume burning process；Whole process it is involved
Formula is as follows：
Combinatorial formula (6)
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, and second formula is second
Coyote hole combustion gas growth equation formula, the 3rd formula is the first coyote hole burning velocity equation, and the 4th formula fires for the second coyote hole
Gas growth equation formula, the 5th formula is the first coyote hole pressure equation formula, and the 6th formula is the second coyote hole equation of gas state
Formula, the 7th formula is afterbody self-destruction unit momentum equation, and the 8th formula is afterbody self-destruction unit motion equation, the 9th
Individual formula is the second coyote hole energy conservation equation formula, and the tenth formula is calculates end condition equation, and the 11st formula is
The calculation equation of the enthalpy of second coyote hole unit mass gas outflow；
In above formula,p 1For the pressure of the first coyote hole,V 0For the original volume of the first coyote hole；ω 1For the first coyote hole
Explosive payload,w 2The gaseous mass percentage revealed for the second coyote hole due to the gap between afterbody self-destruction unit and body,T 2For
Gas temperature in two coyote holes,l 2Displacement for afterbody self-destruction unit with respect to middle part self-destruction unit,φ 2For afterbody self-destruction unit
Secondary work(design factor,m 2For the quality of afterbody self-destruction unit,v 2For the speed of afterbody self-destruction unit,S 2For afterbody self-destruction unit
Sectional area,p 0For external atmospheric pressure,R g =R/M, M are the molal weight of starting mix burning gases；
3rd, the phase III, from middle part self-destruction unit setting in motion untill the second coyote hole gunpowder completely burning, this stage two
Individual coyote hole gunpowder burns 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, and second formula is first
Coyote hole burning velocity equation, the 3rd formula is the second coyote hole combustion gas growth equation formula, and the 4th formula is the first coyote hole gas
Body equation of state, the 5th formula is the second coyote hole gas equation, and the 6th formula is the dynamic of middle part self-destruction unit
Equation is measured, the 7th formula is middle part self-destruction unit motion equation, and the 8th formula is the momentum side of afterbody self-destruction unit
Formula, the 9th formula is afterbody self-destruction unit motion equation, and the tenth formula is the energy conservation equation formula of the first coyote hole,
11st formula is the energy conservation equation formula of the second coyote hole, and the 12nd formula flows out for the first coyote hole unit mass gas
Enthalpy calculation equation；
In above formula,w 1For gas matter of first coyote hole due to the gap leakage between middle part self-destruction unit and body
Measure percentage,T 1For the gas temperature in the first coyote hole,l 1Displacement for middle part self-destruction unit with respect to body,φ 1For middle part self-destruction
The secondary work(design factor of unit,m 1For the quality of middle part self-destruction unit,v 1For the speed of middle part self-destruction unit,S 1For middle part certainly
The sectional area of unit is ruined,h 1The enthalpy flowed out for the first coyote hole unit mass gas；
4th, fourth stage, untill being burned to the burning completely of the first coyote hole starting mix completely from the second coyote hole starting mix, this rank
The second coyote hole expansion work of section, the burning of the first coyote hole 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 velocity equation, and second formula is first
Coyote hole gas equation, the 3rd formula is the second coyote hole gas equation, and the 4th formula is the first coyote hole
Energy conservation equation formula, the 5th formula is the energy conservation equation formula of the second coyote hole, and the 6th formula is middle part self-destruction unit
Momentum equation, the 7th formula is middle part self-destruction unit motion equation, and the 8th formula is the dynamic of afterbody self-destruction unit
Equation is measured, the 9th formula is afterbody self-destruction unit motion equation；
5th, the 5th stage, from the burning completely of the first coyote hole starting mix to middle part self-destruction unit motion to body middle part specific bit
Put, now two coyote hole expansion works；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, and second formula is second
Coyote hole gas equation, the 3rd formula is the energy conservation equation formula of the first coyote hole, and the 4th formula is the second coyote hole
Energy conservation equation formula, the 5th formula is the momentum equation of middle part self-destruction unit, and the 6th formula is middle part self-destruction list
First equation of motion, the 7th formula is the momentum equation of afterbody self-destruction unit, and the 8th formula is transported for afterbody self-destruction unit
Dynamic equation；The motion conditions of middle part self-destruction unit can be calculated by the solution of equation group；
6th, the 6th stage from middle part self-destruction unit motion to body middle part specified location to afterbody self-destruction unit motion to bullet
Self-destruction unit is no longer transported because the associated mechanisms of body medium position are blocked in the middle part of body afterbody specified location moment, this stage
Arrive, the first coyote hole gas is no longer revealed, afterbody self-destruction unit is still moved under gaseous-pressure promotion；Whole process it is involved
Formula it is as follows：
Combinatorial formula (10)
In this combinatorial formula, first formula is the second coyote hole gas equation, and second formula is afterbody
The equation of momentum of self-destruction unit, the 3rd formula is afterbody self-destruction unit motion equation, and the 4th formula is the energy of the second coyote hole
Conservation equation formula is measured, the motion conditions of afterbody self-destruction unit are gone out by solving equations；
In above-mentioned six equation groups：
（1）, the first coyote hole flow into the second coyote hole gas-flow equation：
Formula（11）
In formula,φ 3For Flux Loss coefficient,S 1' it is middle part self-destruction unit and the section in body gap
Product,m 1' be the first coyote hole gaseous mass；
（2）, the second coyote hole flow into extraneous gas-flow equation：
Formula（12）
In formula,φ 4For Flux Loss coefficient,S 2' it is afterbody self-destruction unit and the section in body gap
Product,m 2' be the second coyote hole gaseous mass；
Preferably, in step 4, considering in the middle part of influence and after the principal element of afterbody self-destruction unit motion, to utilize illiteracy
Special Carlow method is entered to middle part self-destruction unit motion to middle card position, afterbody self-destruction unit motion to the reliability of nozzle hole position
Row analysis, is comprised the following steps that：
A, self-destroying function is realized, it is necessary to which card position in body can be pushed into by middle part self-destruction unit by meeting starting mix, simultaneously
Afterbody self-destruction unit is pushed into nozzle hole position, if specified location is to middle part self-destruction unit bottom under original state in the middle part of body
Distance isY 1, the distance of afterbody specified location to original state lower tail self-destruction unit bottom isY 2；
In actual production process, due to accuracy of manufacture influence,Y 1WithY 2To obey the parameter of certain regularity of distribution, and it is overallY 1Distribution function bef 1y (y), it is overallY 2Distribution function bef 2y (y)；
B, according to above-mentioned mathematical modeling by emulation middle part self-destruction unit can be drawn with respect to body displacementL 1A sample
And afterbody self-destruction unit is with respect to body displacementL 2A sample, andL 1Distribution function bef 1Z (z), it is overallL 2Distribution function
Forf 2Z (z)；Middle part, afterbody self-destruction unit can be reached in the middle part of body, afterbody specified location needs to meetL 1≥Y 1, andL 2≥Y 2This
One condition；
C, determine after the distribution functions of parameters, then can determine that specified location in the middle part of middle part self-destruction unit motion to body
Reliability, calculation formula is：
Formula（13）
The reliability that afterbody self-destruction unit motion to afterbody is specified, calculation formula is：
Formula（14）
Self-destruction unit motion to the topmost factor of specified location is influenceed to be filled for the gunpowder of the first coyote hole and the second coyote hole
Doseω 1、ω 2And the quality of middle part self-destruction unit and afterbody self-destruction unitm 1、m 2；
In a practical situation, explosive payloadω i (i=1,2) and each self-destruction element qualitym j (j=1,2) there is certain value
Scope 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 Form j Nominal value, △m j Form j Tolerance；
D, the mathematical modeling based on interior ballistics theory according to foundation, can obtain self-destruction unit with respect to body displacement and two
The chamber charge and the related function of two self-destruction element qualities：
L 1=D(ω i ,m j ) formula（17）
L 2=D(ω i ,m j ) formula（18）
Formula（17）It is the displacement with middle part self-destruction unit with respect to bodyL 1For the equation group of unknown function, formula（18）It is afterbody
Displacement of the self-destruction unit with respect to bodyL 2For the equation group of unknown function；
E, to obtainL 1、L 2Distribution function, emulated by Monte Carlo, according to dress of the actual conditions to two coyote holes
The requirement that the regularity of distribution and quality inspection that the quality of dose and middle part self-destruction unit and afterbody self-destruction unit is obeyed by it are checked on is carried out
Random sampling；Such askDuring secondary sampling, extract on requestω ik (i=1,2) andm jk (j=1,2)；Willω ik Withm jk Substitution (17),
（18）Formula（Solve equation group）The operation of the displacement of middle part self-destruction unit body relative with afterbody self-destruction unit can be respectively obtained
Sample value：
l 1k =D(ω 1k ,ω 2k ;m 1k , m 2k ) formula（19）
l 2k =D(ω 1k ,ω 2k ;m 1k , m 2k ) formula（20）
After simulation model runs n times, shift value of N number of middle part self-destruction unit with respect to body can obtainl 1k (k=1,2,…,
N) and N number of afterbody self-destruction unit with respect to body shift valuel 2k (k=1,2,…,N)；Respectively with N number ofl 1k Value is used as middle part self-destruction
Unit is with respect to body displacementL 1A sample, with N number ofl 2k Value is as afterbody self-destruction unit with respect to body displacementL 2A sample
This, draws the frequency histogram of two parameters, and rough estimate or assumes totalityL 1Distribution functionf 1Z (z), it is overallL 2's
Distribution functionf 2Z (z)；Assume to test by sample data pair, to judge the correctness assumed；It is final to determine distribution functionf 1Z (z) and distribution functionf 2Z (z), by（13）（14）Formula can determine that rain bullet middle part and afterbody self-destruction unit motion to specified location
Reliability；
Preferably, as shown in Fig. 2 in step e, the process of Monte Carlo emulation is comprised the following steps that：
Step（One）, clear and definite sampling parametricω 1、ω 2Withm 1、m 2And simulation times N；
Step（Two）, using Monte Carlo to parameter carry out random sampling；
Step（Three）, by parameterω 1、ω 2Withm 1、m 2And simulation times N substitutes into formula（17）、（18）Solved；
Step（Four）, outputl 1k 、l 2k And it is stored in arrayL 1、L 2；N during this >=k；
Step（Five）, output arrayL 1、L 2, and draw distribution function；
Determine in the middle part of rain bullet and after afterbody self-destruction unit motion to the reliability of specified location, with reference to rain bullet self-destroying function
Reliability logic block diagram can solve the reliability of rain bullet self-destroying function, complete the reliability assessment to rain bullet self-destroying function；Knot
Close reliability design of the assessment result to rain bullet self-destroying function and propose related opinions；
Task pattern of the invention according to the type self-destroyed rain bullet self-destroying function of gunpowder, it is proposed that pass through fixed gunpowder self-destruction
Unit and motion gunpowder self-destruction unit realize the rain bullet system self-destruction fail-safe analysis assessment models of self-destruction, with classical interior ballistics
Learn principle and construct the simulation model of motion gunpowder self-destruction unit reliability, and then the reliability of rain bullet self-destroying function is set
Meter, analysis are with assessing；
The present invention is not limited to above-mentioned specific embodiment, one of ordinary skill in the art from above-mentioned design,
Without performing creative labour, made a variety of conversion are within the scope of the present invention.
Claims (1)
1. a kind of analysis and appraisal procedure for rain bullet self-destroying function reliability, it is characterised in that：Including step in detail below：
Step 1: building the logic diagram model of rain bullet self-destroying function subsystem reliability according to operation principle；
Step 2: analyzing influence realizes the key factor of rain bullet self-destroying function reliability, and collect the reliable of priming system component
Degrees of data；
Step 3: self-destruction unit and afterbody self-destruction unit need starting mix gaseous-pressure to promote to descend to move to finger in the middle part of rain bullet
The operation principle put is positioned, based on classical interior ballistics theory, the mathematical modeling of the rain bullet self-destruction course of work is set up；
Step 4: using monte carlo method to the worked of middle part self-destruction unit and afterbody self-destruction unit motion to specified location
Cheng Jinhang analog simulations, solve middle part self-destruction unit and afterbody self-destruction unit motion to the reliability data of specified location；
Step 5: combining the rain bullet self-destroying function reliability block diagram model built to rail-proof bomb self-destroying function reliability of increasing rain artificially
Analyzed and assessed；
In step one, the reliability block diagram model includes gunpowder extension unit R1, head self-destruction unit RT, middle part self-destruction list
First RC, afterbody self-destruction unit RA；The gunpowder extension unit R1With head self-destruction unit RTInput end signal connection；The head
Portion's self-destruction unit RTOutput end and middle part self-destruction unit RCInput end signal connection；The middle part self-destruction unit RCOutput
End and afterbody self-destruction unit RAInput end signal connection；
Wherein, gunpowder extension unit R1Including electric igniter R11, electric igniter R12, electric igniter R17, electric igniter R18, flame agent
Dark Ignition has R13, flame agent Dark Ignition tool R14, flame agent R15, one-level Delay Tube R16, electric Delay Tube R19, starting mix R110；It is described
Electric igniter R11Has R with flame agent Dark Ignition13Signal is connected；The electric igniter R12Has R with flame agent Dark Ignition14Signal connects
Connect；The flame agent Dark Ignition has R13, flame agent Dark Ignition tool R14With flame agent R15Signal is connected；The flame agent R15With one-level
Delay Tube R16Signal is connected；The one-level Delay Tube R16, electric Delay Tube R19With starting mix R110Signal is connected；The electric ignition
Head R11Respectively with electric igniter R17, electric igniter R18Signal is connected；
According to the reliability block diagram, set the Reliable Mathematics model of rain bullet self-destroying function as：
R=R1RTRCRAFormula (1)
The reliability of gunpowder extension unit is can obtain by self-destruction subsystem reliability logic diagram, its formula of reliability is：
R1={ 1- { 1- [1- (1-R11R13)×(1-R12R14)]×R15×R16}×{1-[1-(1-R17)×(1-R18)]×R19}}
×R110
Formula (2)
Wherein R11、R12、R17、R18It is expressed as the reliability of electric igniter, R13、R14It is expressed as the reliable of flame agent Dark Ignition tool
Degree, R15It is expressed as the reliability of flame agent, R16It is expressed as the reliability of one-level Delay Tube, R19The reliability of electric Delay Tube is expressed as,
R110It is expressed as the reliability of starting mix；The reliability R of middle part self-destruction unitCAnd the reliability R of afterbody self-destruction unitABe then this two
The reliability of individual unit priming system is moved to the synthesis of the reliability of designated displacement, i.e.,
RC=RC1×RCLFormula (3)
RA=RA1×RALFormula (4)
Wherein, RC1For the reliability of middle part self-destruction unit priming system, RA1For the reliability of afterbody self-destruction unit priming system；
According to rain bullet self-destroying function realize operation principle, reliability model structure is simplified, the model including body,
First coyote hole, the second coyote hole, middle part self-destruction body, afterbody self-destruction body, middle card and nozzle hole；First coyote hole is arranged on body
Foremost, and close to middle part self-destruction body set；Second coyote hole is arranged on the middle part of middle part self-destruction body and afterbody self-destruction body；Institute
State the middle middle part blocked and be arranged on body；The nozzle hole is arranged on the tail end of body；
In step 2, the key factor of rain bullet self-destroying function reliability, which is analyzed, to be realized to influence, mainly including following several
Individual aspect：
(1), the reliability block diagram set up by step one determines that influence realizes that the key factor of rain bullet self-destroying function reliability is
Two self-destruction unit motions to specified location reliability；
(2), the reliability model set up according to step one clearly requires the priming system component reliability data of collection；In step
In rapid three, constructed mathematical modeling is based primarily upon hypothesis below：
(1), whether priming system transmittance process at different levels are reliable in self-destruction unit, and whether self-destruction body brisance meets requirement, that is, exists
Specified location, which can scrap body completely when detonating, to be determined by the reliability of priming system element；
(2), self-destruction unit motion to the process of specified location meets the related hypothesis of classical interior ballistics model；
(3), because the small aperture for connecting two coyote holes is minimum, its circulated gases is minimum with respect to for coyote hole gas gross therefore false
If without gas flowing in aperture；
(4), self-destruction unit does not consider the sticky of gas, gas to body inwall by the process that starting mix is pushed to specified location
Friction, combustion gas uses indirect method amendment to the heat transfer of chamber wall；
(5), powder burning and two self-destruction unit motions are carried out all under average pressure；
(6) it is also minimum to ballistic impact, in view of middle part self-destruction unit speed is minimum during self-destruction unit in the middle part of middle block, therefore
It is approximately considered in the heat produced by this process on the second coyote hole without influence；
(7), because starting mix burning and self-destruction unit motion process time are very short, much smaller than the extension time of Delay Tube, therefore its
Self-destruction time and height are only relevant with the outer trajectory that delay time and the rain bullet of Delay Tube are launched, and its influence is ignored herein；
Constructed mathematical modeling mainly includes the following six stage, specific as follows：
1), the first stage, since the gunpowder of two coyote holes (i.e. the first coyote hole, the second coyote hole) being burned to afterbody self-destruction unit opens
Begin to move, this process is coyote hole constant volume burning process；The involved formula of whole process is as follows：
This is a composite class formula, wherein, first formula is the second coyote hole combustion gas growth equation formula, and second formula is the
Two coyote hole burning velocity equations, the 3rd formula is the second coyote hole pressure equation formula；
In above formula, p2For the pressure of the second coyote hole, V0' be the second coyote hole original volume；ω2For the dress of the second coyote hole
Dose；With starting pressure P0For initial value, the second coyote hole powder burning situation is solved；
2), second stage, 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 burning of the second coyote hole and external expansion work, the first coyote hole constant volume burning process；The involved formula of whole process
It is as follows：
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, and second formula is the second coyote hole
Combustion gas growth equation formula, the 3rd formula is the first coyote hole burning velocity equation, and the 4th formula is given birth to for the second coyote hole combustion gas
Into equation, the 5th formula is the first coyote hole pressure equation formula, and the 6th formula is the second coyote hole gas equation, the
Seven formula are afterbody self-destruction unit momentum equation, and the 8th formula is afterbody self-destruction unit motion equation, the 9th public affairs
Formula is the second coyote hole energy conservation equation formula, and the tenth formula is calculates end condition equation, and the 11st formula is second
The calculation equation of the enthalpy of coyote hole unit mass gas outflow；
In above formula, p1For the pressure of the first coyote hole, V0For the original volume of the first coyote hole；ω1For the powder charge of the first coyote hole
Amount, w2For gaseous mass percentage of second coyote hole due to the gap leakage between afterbody self-destruction unit and body, T2For the second medicine
Indoor gas temperature, l2Displacement for afterbody self-destruction unit with respect to middle part self-destruction unit, φ2For the secondary of afterbody self-destruction unit
Work(design factor, m2For the quality of afterbody self-destruction unit, v2For the speed of afterbody self-destruction unit, S2For cutting for afterbody self-destruction unit
Area, p0For external atmospheric pressure, Rg=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 untill the second coyote hole gunpowder completely burning, two medicines of this stage
Room gunpowder burns and expansion work；The involved formula of whole process is as follows：
In this combinatorial formula, first formula is the first coyote hole combustion gas growth equation formula, and second formula is the first coyote hole
Burning velocity equation, the 3rd formula is the second coyote hole combustion gas growth equation formula, and the 4th formula is the first coyote hole gas shape
State equation, the 5th formula is the second coyote hole gas equation, and the 6th formula is the momentum side of middle part self-destruction unit
Formula, the 7th formula is middle part self-destruction unit motion equation, and the 8th formula is the momentum equation of afterbody self-destruction unit,
9th formula is afterbody self-destruction unit motion equation, and the tenth formula is the energy conservation equation formula of the first coyote hole, the tenth
One formula is the energy conservation equation formula of the second coyote hole, and the 12nd formula is the enthalpy that the first coyote hole unit mass gas flows out
Calculation equation；
In above formula, w1For gaseous mass hundred of first coyote hole due to the gap leakage between middle part self-destruction unit and body
Fraction, T1For the gas temperature in the first coyote hole, l1Displacement for middle part self-destruction unit with respect to body, φ1For middle part self-destruction unit
Secondary work(design factor, m1For the quality of middle part self-destruction unit, v1For the speed of middle part self-destruction unit, S1For middle part self-destruction list
The sectional area of member, h1The enthalpy flowed out for the first coyote hole unit mass gas；
4), fourth stage, untill being burned to the burning completely of the first coyote hole starting mix completely from the second coyote hole starting mix, this stage the
Two coyote hole expansion works, the burning of the first coyote hole and expansion work；The involved formula of whole process is as follows：
In this combinatorial formula, first formula is the first coyote hole burning velocity equation, and second formula is the first coyote hole
Gas equation, the 3rd formula is the second coyote hole gas equation, and the 4th formula is the energy of the first coyote hole
Conservation equation formula, the 5th formula is the energy conservation equation formula of the second coyote hole, and the 6th formula is the dynamic of middle part self-destruction unit
Equation is measured, the 7th formula is middle part self-destruction unit motion equation, and the 8th formula is the momentum side of afterbody self-destruction unit
Formula, the 9th formula is afterbody self-destruction unit motion equation；
5), the 5th stage, from the burning completely of the first coyote hole starting mix to middle part self-destruction unit motion to body middle part specified location,
Now two coyote hole expansion work；The involved formula of whole process is as follows：
In this combinatorial formula, first formula is the first coyote hole gas equation, and second formula is the second coyote hole
Gas equation, the 3rd formula is the energy conservation equation formula of the first coyote hole, and the 4th formula is the energy of the second coyote hole
Conservation equation formula is measured, the 5th formula is the momentum equation of middle part self-destruction unit, and the 6th formula is transported for middle part self-destruction unit
Dynamic equation, the 7th formula is the momentum equation of afterbody self-destruction unit, and the 8th formula is afterbody self-destruction unit motion side
Formula；The motion conditions of middle part self-destruction unit can be calculated by the solution of equation group；
6), the 6th stage from middle part self-destruction unit motion to body middle part specified location to afterbody self-destruction unit motion to body tail
Self-destruction unit is no longer transported to because the associated mechanisms of body medium position are blocked in the middle part of portion's specified location moment, this stage,
First coyote hole gas is no longer revealed, and afterbody self-destruction unit is still moved under gaseous-pressure promotion；Whole process it is involved
Formula is as follows：
In this combinatorial formula, first formula is the second coyote hole gas equation, and second formula is afterbody self-destruction
The equation of momentum of unit, the 3rd formula is afterbody self-destruction unit motion equation, and the 4th formula is kept for the energy of the second coyote hole
Permanent equation, the motion conditions of afterbody self-destruction unit are gone out by solving equations；
In above-mentioned six equation groups：
(1), the first coyote hole flows into the gas-flow equation of the second coyote hole：
In formulaφ3For Flux Loss coefficient, S1' it is middle part self-destruction unit and the sectional area in body gap, m1′
For the gaseous mass of the first coyote hole；
(2), the second coyote hole flows into extraneous gas-flow equation：
In formulaφ4For Flux Loss coefficient, S2' it is afterbody self-destruction unit and the sectional area in body gap, m2′
For the gaseous mass of the second coyote hole；
In step 4, considering in the middle part of influence and after the principal element of afterbody self-destruction unit motion, to utilize monte carlo method
To middle part self-destruction unit motion to middle card position, afterbody self-destruction unit motion to the reliability of nozzle hole position is analyzed, tool
Body step is as follows：
A, self-destroying function is realized, it is necessary to card position in body can be pushed into by middle part self-destruction unit by meeting starting mix, while by tail
Portion's self-destruction unit is pushed into nozzle hole position, if distance of the specified location to middle part self-destruction unit bottom under original state in the middle part of body
For Y1, the distance of afterbody specified location to original state lower tail self-destruction unit bottom is Y2；
In actual production process, due to accuracy of manufacture influence, Y1With Y2To obey the parameter of certain regularity of distribution, and totality Y1's
Distribution function is f1y(y), overall Y2Distribution function be f2y(y)；
B, according to above-mentioned mathematical modeling by emulation middle part self-destruction unit can be drawn with respect to body displacement L1A sample and afterbody
Self-destruction unit is with respect to body displacement L2A sample, and L1Distribution function be f1Z(z), overall L2Distribution function be f2Z
(z)；Middle part, afterbody self-destruction unit can be reached in the middle part of body, afterbody specified location need to meet L1≥Y1, and L2≥Y2This
Part；
C, determine after the distribution functions of parameters, then can determine that specified location in the middle part of middle part self-destruction unit motion to body can
By degree, calculation formula is：
The reliability that afterbody self-destruction unit motion to afterbody is specified, calculation formula is：
Self-destruction unit motion is influenceed to the powder charge amount that the topmost factor of specified location is the first coyote hole and the second coyote hole
ω1、ω2And the quality m of middle part self-destruction unit and afterbody self-destruction unit1、m2；
In a practical situation, explosive payload ωiAnd each self-destruction element quality m (i=1,2)j(j=1,2) there is certain value model
Enclose and statistical nature：ωip-△ωi≤ωi≤ωip+△ωiFormula (15)
mjp-△mj≤mj≤mjp+△mjFormula (16)
ω in formulaipFor ωiNominal value, △ ωiFor ωiTolerance, mjpFor mjNominal value, △ mjFor mjTolerance；
D, the mathematical modeling based on interior ballistics theory according to foundation, can obtain self-destruction unit with respect to body displacement and two coyote holes
Explosive payload and the related function of two self-destruction element qualities：
L1=D (ωi,mj) formula (17)
L2=D (ωi,mj) formula (18)
Formula (17) is the displacement L with respect to body with middle part self-destruction unit1For the equation group of unknown function, formula (18) is afterbody self-destruction
Displacement L of the unit with respect to body2For the equation group of unknown function；
E, to obtain L1、L2Distribution function, emulated by Monte Carlo, according to explosive payload of the actual conditions to two coyote holes
The requirement that the regularity of distribution and quality inspection obeyed with the quality of middle part self-destruction unit and afterbody self-destruction unit by it are checked on carries out random
Sampling；Such as in kth time sampling, ω is extracted on requestikAnd m (i=1,2)jk(j=1,2)；By ωikAnd mjkSubstitution (17),
(18) formula (solving equation group) can respectively obtain the operation of the displacement of middle part self-destruction unit body relative with afterbody self-destruction unit
Sample value：
l1k=D (ω1k,ω2k；m1k,m2k) formula (19)
l2k=D (ω1k,ω2k；m1k,m2k) formula (20)
After simulation model runs n times, shift value l of N number of middle part self-destruction unit with respect to body can obtain1k(k=1,2 ..., N)
And N number of afterbody self-destruction unit is with respect to the shift value l of body2k(k=1,2 ..., N)；Respectively with N number of l1kValue is used as middle part self-destruction list
The relative body displacement L of member1A sample, with N number of l2kValue is as afterbody self-destruction unit with respect to body displacement L2A sample,
Draw the frequency histogram of two parameters, and rough estimate or assume overall L1Distribution function f1Z(z), overall L2Distribution
Function f2Z(z)；Assume to test by sample data pair, to judge the correctness assumed；It is final to determine distribution function f1Z
And distribution function f (z)2Z(z), can determine that in the middle part of rain bullet by (13) (14) formula and afterbody self-destruction unit motion to specified location can
By degree；
In step e, the process of Monte Carlo emulation is comprised the following steps that：
Step (1), clear and definite sampling parametric ω1、ω2And m1、m2And simulation times N；
Step (2), using Monte Carlo to parameter carry out random sampling；
Step (3), by parameter ω1、ω2And m1、m2And simulation times N substitutes into formula (17), (18) and solved；
Step (4), output l1k、l2kAnd it is stored in array L1、L2；N >=k during this；
Step (5), output array L1、L2, and draw its distribution function.
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 CN104298848A (en) | 2015-01-21 |
| CN104298848B true 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) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
-
2014
- 2014-07-15 CN CN201410335207.4A patent/CN104298848B/en not_active Expired - Fee Related
Non-Patent Citations (4)
| Title |
|---|
| BL-1A增雨防雹火箭系统安全使用探讨;卢培玉 等;《山东气象》;20100331;第30卷(第1期);全文 * |
| 基于蒙特卡洛仿真的液压系统动态可靠性;苏春 等;《东南大学学报》;20060531;第36卷(第3期);第三段 * |
| 增雨防雹火箭弹的工作原理和技术参数;魏旭辉 等;《沙漠与绿洲气象》;20090831;第3卷;全文 * |
| 增雨防雹火箭的可靠性分配及其单元可靠性增长试验;王铮;《固体火箭技术》;20060228;第29卷(第2期);第1、2、3部分 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104298848A (en) | 2015-01-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103630363B (en) | Test method and test device for simulating high altitude ignition ability of turbine engine | |
| CN104298848B (en) | A kind of analysis and appraisal procedure for rain bullet self-destroying function reliability | |
| CN106052491B (en) | A kind of simulated explosion shock wave and high speed fragments group's connected load Experiment on Function method | |
| CN103267644A (en) | Engine performance simulation method | |
| CN105446167B (en) | Hypersonic scramjet engine real-time model, emulation mode | |
| CN107944150A (en) | Interior outflow multicomponent chemical reaction flow field integration numerical simulator | |
| CN106599400A (en) | Fast calculation and dynamic simulation method of aircraft tail flame infrared radiation | |
| CN109376445A (en) | Gas-turbine unit starts modeling method | |
| Li et al. | A numerical study of flow field evolution in a pulsed detonation engine | |
| CN109408915A (en) | Solid-rocket scramjet engine Combustion Flow Field emulation mode | |
| Wessam et al. | Firing precision evaluation for unguided artillery projectile | |
| CN106323661A (en) | Testing method of tunnel smoke vent opening concentrated smoke evacuation performance testing device | |
| CN106844820B (en) | Carbon dioxide inflow parameter determination method based on deep space exploration test | |
| CN104897407A (en) | Engine stability scaling method | |
| CN109738196B (en) | Ramjet performance space-ground conversion method based on ground direct connection test | |
| CN208283323U (en) | Gas in a kind of simulation coal mine down-hole tunnel, coal-dust explosion experimental provision | |
| CN109870080A (en) | The smoke screen of explosive type IR flare smoke grenade generates model | |
| Paxson | A simplified model for detonation based pressure-gain combustors | |
| US6519554B1 (en) | Computer implemented system and method for evaluating gas generator launchers | |
| CN205991881U (en) | Tail room for engine testsand gas extraction system | |
| CN205991883U (en) | Tail room for engine testsand gas extraction system | |
| CN106198030A (en) | Tail room for engine testsand gas extraction system | |
| Wilhelm et al. | Flight simulation of a hybrid projectile to estimate the impact of launch angle on range extension | |
| Thomas | The 9/11 truth movement: The top conspiracy theory, a decade later | |
| Cao | Parametric cycle analysis of continuous rotating detonation ejector-augmented rocket engine |
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 |
Granted publication date: 20170728 Termination date: 20210715 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |