CN112797856A - Method for rapidly evaluating position load of minimum risk bomb of transport aircraft - Google Patents

Method for rapidly evaluating position load of minimum risk bomb of transport aircraft Download PDF

Info

Publication number
CN112797856A
CN112797856A CN202110131766.3A CN202110131766A CN112797856A CN 112797856 A CN112797856 A CN 112797856A CN 202110131766 A CN202110131766 A CN 202110131766A CN 112797856 A CN112797856 A CN 112797856A
Authority
CN
China
Prior art keywords
load
equivalent
mach
overpressure
explosive
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
CN202110131766.3A
Other languages
Chinese (zh)
Other versions
CN112797856B (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.)
Civil Aviation University of China
Original Assignee
Civil Aviation University of China
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 Civil Aviation University of China filed Critical Civil Aviation University of China
Priority to CN202110131766.3A priority Critical patent/CN112797856B/en
Publication of CN112797856A publication Critical patent/CN112797856A/en
Application granted granted Critical
Publication of CN112797856B publication Critical patent/CN112797856B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42—AMMUNITION; BLASTING
    • F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
    • F42B35/00—Testing or checking of ammunition

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses a method for quickly evaluating the position load of a minimum risk bomb of a transport plane, which comprises the following steps of firstly, placing a detonation center at a selected position; secondly, calculating pressure pulses acting on internal points of the fuselage; finally, the dynamic load is converted into an equivalent static load for predicting structural failure caused by explosion. According to the invention, the equivalent static pressure is obtained through calculation, the scale of structural damage caused by explosive explosion can be estimated through the value, and the estimated value is utilized to assist in evaluating whether the residual structure can continuously bear the flight load.

Description

Method for rapidly evaluating position load of minimum risk bomb of transport aircraft
Technical Field
The invention relates to the field of rapid evaluation of minimum risk bombs, in particular to a rapid evaluation method for position loads of minimum risk bombs of a transportation type airplane.
Background
The federal aviation administration of the united states issued a FAR 25-127 amendment on 28/10/2008, requiring manufacturers to design a "minimum risk bomb location" on the aircraft for placement of discovered suspect devices. Through the comprehensive design of the structure and the system, the key structure and the system related to the safety of the airplane can be protected to the maximum extent once the key structure and the system are exploded, so that the safety level of the airplane is effectively improved.
The validity of the aircraft minimum risk bomb position (LRBL) does not need to be verified experimentally, i.e. the validity of the position does not need to be verified by exploding the aircraft with "real bombs". However, depending on the design, other experiments related to LRBL are required, so how to establish a method for quickly evaluating the location load of the minimum risk bomb is an effective means for supporting the LRBL design and compliance verification.
Disclosure of Invention
Aiming at the existing problems, the invention provides a method for quickly evaluating the position load of the minimum risk bomb of the transport aircraft.
The technical solution for realizing the purpose of the invention is as follows:
the method for rapidly evaluating the position load of the minimum risk bomb of the transport-type airplane is characterized by comprising the following steps of:
step 1: placing the detonation center at a selected position, recording coordinates of the position as (a, b, c), defining parameters of the airplane body, and initializing the parameters: radius R, length L, thickness t and explosive equivalent w;
step 2: based on the coordinates and initialization parameters of the calculated point in the fuselage, the relevant parameters of the shock wave acting on the point in the fuselage and the load holding time t are calculated0And the shock wave related parameters include: mach number M of incident wavexReflected overpressure prMach number M of Mach wavesMach overpressure pmsAnd maximum Mach angle of reflection betamax；
And step 3: by calculating the maximum overpressure pmaxThe dynamic coefficient y is used for converting the dynamic load generated by explosion into the equivalent static load under the equivalent quantity and predicting the structural damage caused by explosion;
and 4, step 4: according to the maximum overpressure p obtainedmaxAnd the dynamic coefficient y, calculating to obtain the equivalent static pressure based on the detonation distance and the explosive equivalent;
and 5: the stress condition of the aircraft structure after the bomb explosion can be solved according to the obtained equivalent static pressure, and then whether the residual strength of the aircraft structure can continuously bear the flight load or not is evaluated in an auxiliary mode.
Further, the incident wave Mach number M obtained in the step 2xReflection overpressure prMach number M of Mach wavesMach overpressure pmsAngle of maximum Mach reflection betamaxAnd time of load t0Respectively as follows:
Figure BDA0002925605260000021
wherein, according to the cube root law of proportionality: if two same explosive charges are provided, the geometrical shapes of the charges are similar to each other, the charges have different sizes, and when the charges are exploded in the same atmosphere, the charges are at the same proportional distance
Figure BDA0002925605260000031
Similar shock waves are generated, D is a space distance, and w is the equivalent weight of the explosive;
Figure BDA0002925605260000032
Figure BDA0002925605260000033
Ms＝Mx/sinβ (4)，
wherein beta is the contained angle between the ray that sends out from the initiation point and the fuselage profile intersection point tangent line, b is the horizontal radial distance of explosive and equivalent static load point to:
β＝90-arctan(R/b) (5)，
Figure BDA0002925605260000034
Figure BDA0002925605260000035
Figure BDA0002925605260000036
further, the maximum overpressure p of step 3maxComprises the following steps:
Figure BDA0002925605260000037
further, modeling the aircraft structure as a one-dimensional system with the natural frequency as a characteristic, wherein the static load in step 3 is the static load which generates the same deformation as the dynamic explosion, and estimating the natural frequency T of the structure, wherein the estimation formula is as follows:
Figure BDA0002925605260000038
wherein f is the cycle frequency;
and:
Figure BDA0002925605260000039
wherein t is the skin thickness, C is the frequency constant, k is the material factor, and l is the distance between the reinforcement frames.
Further, step 4 calculates the equivalent static pressure P for each of the explosive distance and the explosive equivalentEThe formula of (1) is:
PE＝pmaxy (12)。
compared with the prior art, the method has the following beneficial effects:
firstly, the damage range of the explosion shock wave obtained by calculation on the aircraft skin can meet the requirement of any precision;
secondly, the invention can quickly provide the equivalent static pressure calculated by each explosion distance and the explosive equivalent, estimate the structural damage scale caused by explosive explosion according to the obtained equivalent static pressure, and assist in evaluating whether the residual structure can continuously bear the flight load by using the estimated value.
Drawings
FIG. 1 is a process flow diagram of an evaluation method of the present invention;
FIG. 2 is a schematic diagram of a Mach number versus proportional distance fit curve in the present invention;
FIG. 3 is t0T and pE/pmaxThe graph of (a), namely a curve diagram of the kinetic coefficient y;
Detailed Description
In order to make those skilled in the art better understand the technical solution of the present invention, the following further describes the technical solution of the present invention with reference to the drawings and the embodiments.
Referring to the attached fig. 1, the method for rapidly evaluating the location load of the minimum risk bomb of the transport-type airplane, provided by the invention, comprises the following operation steps:
step 1: placing the detonation center at a selected position, recording coordinates of the position as (a, b, c), defining parameters of the airplane body, and initializing the parameters: radius R, length L, thickness t and explosive equivalent w;
step 2: based on the coordinates and initialization parameters of the calculated point in the fuselage, the relevant parameters of the shock wave acting on the point in the fuselage and the load holding time t are calculated0And the shock wave related parameters include: mach number M of incident wavexReflected overpressure prMach number M of Mach wavesMach overpressure pmsAnd maximum Mach angle of reflection betamax；
And step 3: by calculating the maximum overpressure pmaxThe dynamic coefficient y is used for converting the dynamic load generated by explosion into the equivalent static load under the equivalent quantity and predicting the structural damage caused by explosion;
wherein, the dynamic load refers to that the explosive load is in a change state (dynamic pressure), and the static load refers to a static load (static pressure); the maximum overpressure is the maximum value of dynamic pressure, and the power coefficient y is a parameter for converting the maximum value of dynamic pressure into static pressure;
and 4, step 4: according to the maximum overpressure p obtainedmaxAnd the dynamic coefficient y, calculating to obtain the equivalent static pressure based on the detonation distance and the explosive equivalent;
and 5: the stress condition of the aircraft structure after the bomb explosion can be solved according to the obtained equivalent static pressure, and then whether the rest structure of the aircraft structure can continuously bear the flight load or not is evaluated in an auxiliary mode.
The equivalent static pressure is acted on the airplane structure, so that the stress of the airplane structure after the bomb explodes can be solved, and the residual strength can be evaluated in an auxiliary manner;
preferably, in step 2, the proportional distance z of any position of the fuselage is solved, and an explicit Shocks in air annex table is fitted to obtain the proportional distance z and the mach number MxThe Mach number M of the incident wave is obtained by the piecewise function betweenxComprises the following steps:
Figure BDA0002925605260000061
wherein, according to Hopkinson's law of proportionality, also called cube-root's law of proportionality: if two same explosive charges are provided, the geometrical shapes of the charges are similar to each other, the charges have different sizes, and when the charges are exploded in the same atmosphere, the charges are at the same proportional distance
Figure BDA0002925605260000062
Generating similar shock waves, wherein D is the space distance, and w is the equivalent weight of the explosive;
according to the obtained incident wave Mach number MxObtaining the overpressure at this point, the reflected overpressure prMach number M of Mach wavesMach overpressure pmsMaximum Mach reflectionAngle betamaxAnd time of load t0Respectively as follows:
Figure BDA0002925605260000063
Figure BDA0002925605260000064
Ms＝Mx/sinβ (4)，
wherein, beta is 90-arctan (R/b) which represents the included angle between a ray drawn from the initiation point and the tangent of the intersection point of the fuselage outline, and b is the horizontal radial distance between the explosive and the equivalent static load point;
Figure BDA0002925605260000065
Figure BDA0002925605260000066
Figure BDA0002925605260000067
preferably, in step 3, said maximum overpressure pmaxComprises the following steps:
Figure BDA0002925605260000068
Figure BDA0002925605260000071
at this time, a static load which is the same as the static load generated by the dynamic explosion in the process of generating deformation is determined, the natural frequency T of the structure must be estimated, and since the deformation of the structure under the dynamic load is related to the natural frequency of the structure (the ratio of the explosion action time to the natural frequency is linearly related to the ratio of the equivalent static load to the dynamic load), the natural frequency of the structure is considered when the dynamic load is equivalent to the static load, and the calculation formula is as follows:
Figure BDA0002925605260000072
wherein f is the cycle frequency;
and:
Figure BDA0002925605260000073
wherein t is the skin thickness, C is the frequency constant, k is the material factor, and l is the distance between the reinforcement frames.
According to the abscissa as t0T, s, ordinate pE/pmaxI.e. the curve of the kinetic coefficient y, the fitted piecewise function is:
Figure BDA0002925605260000074
in FIG. 3, t is0T and pE/pmaxThe abscissa of the graph of (a) is denoted as s, the ordinate is denoted as y, y is linearly related to x, the resulting curve is a curve of the kinetic coefficient y, so the value of s can be determined by a series of equations preceding equation (13), pmaxIt can also be found that P can be obtained according to the formula (13)E。
Further, an equivalent static pressure P can be calculated for each of the explosive distance and the explosive equivalentEComprises the following steps:
PE＝pmaxy (12)。
examples
Assuming that the explosive equivalent is 80g, the skin thickness is 2mm, the distance l between the reinforcing frames is 300mm, and the explosion distance is the incident wave Mach number MxReflection overpressure prMach number M of Mach wavesMach overpressure pms，
Maximum Mach angle of reflection betamaxAnd time of load t0Respectively as follows:
Mx＝20.6266-44.3921z+28.3631z2＝3.4681；
inputting R and b values by a program according to the parameters of the airplane and the coordinate points for calculating the static load,
and calculating: β -90-arctan (R/b) -1.2558;
Figure BDA0002925605260000081
Figure BDA0002925605260000082
Ms＝Mx/sinβ＝Inf；
Figure BDA0002925605260000083
Figure BDA0002925605260000084
Figure BDA0002925605260000085
in this case, beta is not more than betamaxMaximum overpressure pmaxComprises the following steps:
pmax＝pr＝76.7230Psi；
at this time, the natural frequency T of the structure is:
Figure BDA0002925605260000091
s＝t0/T＝0.0641；
the kinetic coefficient y is obtained by equation (13):
y＝0.02+s＝0.0841；
converting the dynamic load into equivalent static load for predicting structural damage caused by explosion and equivalent static pressure PEComprises the following steps:
PE＝pmaxy＝6.4563Psi。
according to the obtained equivalent static pressure PEThe size of the structural damage caused by the explosive blast is estimated and this estimate is then used to assist in assessing whether the remaining structure can continue to withstand the flight loads.
Those not described in detail in this specification are within the skill of the art. Although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that various changes in the embodiments and modifications of the invention can be made, and equivalents of some features of the invention can be substituted, and any changes, equivalents, improvements and the like, which fall within the spirit and principle of the invention, are intended to be included within the scope of the invention.

Claims (5)

1. The method for rapidly evaluating the position load of the minimum risk bomb of the transport-type airplane is characterized by comprising the following steps of:
step 1: placing the detonation center at a selected position, recording coordinates of the position as (a, b, c), defining parameters of the airplane body, and initializing the parameters: radius R, length L, thickness t and explosive equivalent w;
step 2: based on the coordinates and initialization parameters of the calculated point in the fuselage, the relevant parameters of the shock wave acting on the point in the fuselage and the load holding time t are calculated0And the shock wave related parameters include: mach number M of incident wavexReflected overpressure prMach number M of Mach wavesMach overpressure pmsAnd maximum Mach angle of reflection betamax；
And step 3: by calculating the maximum overpressure pmaxThe dynamic coefficient y is used for converting the dynamic load generated by explosion into the equivalent static load under the equivalent quantity and predicting the structural damage caused by explosion;
step (ii) of4: according to the maximum overpressure p obtainedmaxAnd the dynamic coefficient y, calculating to obtain the equivalent static pressure based on the detonation distance and the explosive equivalent;
and 5: the stress condition of the aircraft structure after the bomb explosion can be solved according to the obtained equivalent static pressure, and then whether the rest structure of the aircraft structure can continuously bear the flight load or not is evaluated in an auxiliary mode.
2. The method for rapidly estimating the location load of the minimum risk bomb of a transport-type aircraft according to claim 1, wherein the incident wave Mach number M obtained in step 2xReflection overpressure prMach number M of Mach wavesMach overpressure pmsAngle of maximum Mach reflection betamaxAnd time of load t0Respectively as follows:
Figure FDA0002925605250000011
wherein, according to the cube root law of proportionality: if two same explosive charges are provided, the geometrical shapes of the charges are similar to each other, the charges have different sizes, and when the charges are exploded in the same atmosphere, the charges are at the same proportional distance
Figure FDA0002925605250000012
Similar shock waves are generated, D is a space distance, and w is the equivalent weight of the explosive;
Figure FDA0002925605250000021
Figure FDA0002925605250000022
Ms＝Mx/sinβ (4)，
wherein beta is the contained angle between the ray that sends out from the initiation point and the fuselage profile intersection point tangent line, b is the horizontal radial distance of explosive and equivalent static load point to:
β＝90-arctan(R/b) (5)，
Figure FDA0002925605250000023
Figure FDA0002925605250000024
Figure FDA0002925605250000025
3. method for the rapid assessment of the location load of a minimum risk bomb of a transport aircraft according to claim 1, characterised in that the maximum overpressure p of step 3 ismaxComprises the following steps:
Figure FDA0002925605250000026
4. the method for rapidly estimating the location load of the minimum risk bomb of transportation aircraft as claimed in claim 1, wherein the aircraft structure is modeled as a one-dimensional system characterized by natural frequency, the static load in step 3 is the static load which generates the same deformation as the dynamic explosion, the natural frequency T of the structure is estimated, and the estimation formula is as follows:
Figure FDA0002925605250000031
wherein f is the cycle frequency;
and:
Figure FDA0002925605250000032
wherein t is the skin thickness, C is the frequency constant, k is the material factor, and l is the distance between the reinforcement frames.
5. The method for rapidly estimating the position load of the minimum risk bomb of transport-type airplane according to claim 1, wherein the step 4 is to calculate the equivalent static pressure P for each detonation distance and explosive equivalentEThe formula of (1) is:
PE＝pmaxy (12)。
CN202110131766.3A 2021-01-30 2021-01-30 Method for rapidly evaluating position load of minimum risk bomb of transport aircraft Active CN112797856B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110131766.3A CN112797856B (en) 2021-01-30 2021-01-30 Method for rapidly evaluating position load of minimum risk bomb of transport aircraft

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110131766.3A CN112797856B (en) 2021-01-30 2021-01-30 Method for rapidly evaluating position load of minimum risk bomb of transport aircraft

Publications (2)

Publication Number Publication Date
CN112797856A true CN112797856A (en) 2021-05-14
CN112797856B CN112797856B (en) 2023-01-03

Family

ID=75813097

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110131766.3A Active CN112797856B (en) 2021-01-30 2021-01-30 Method for rapidly evaluating position load of minimum risk bomb of transport aircraft

Country Status (1)

Country Link
CN (1) CN112797856B (en)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB719771A (en) * 1950-01-18 1954-12-08 British Electricity Authority Improvements in or relating to tele-communication equipment
US20030112107A1 (en) * 2001-12-04 2003-06-19 Nmr Holdings No. 2 Pty Ltd Asymmetric tesseral shim coils for magnetic resonance
US20040061895A1 (en) * 2002-09-27 2004-04-01 Brother Kogyo Kabushiki Kaisha Image forming device and image forming method
CN101151504A (en) * 2005-04-08 2008-03-26 独立行政法人产业技术综合研究所 Blasting treating method
US20110168004A1 (en) * 2009-10-20 2011-07-14 Henegar Douglas W System and method for mitigating and directing an explosion aboard an aircraft
CN204713434U (en) * 2015-04-02 2015-10-21 中国人民解放军空军勤务学院 Aerial bomb ammunition carrier
CN106197184A (en) * 2016-08-01 2016-12-07 西北工业大学 A kind of civil aircraft main cabin portable directional explosion-proof device
RU2649999C1 (en) * 2017-04-17 2018-04-06 Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом") Method of estimation of fougasseness characteristics in air explosion of a moving test object (variants)
CN109581381A (en) * 2018-11-28 2019-04-05 中国民航大学 Enhanced turbulent flow detection method based on the vertical load factor
CN109780956A (en) * 2019-03-06 2019-05-21 西安近代化学研究所 A kind of cumulative pressure release civil aircraft main cabin directional explosion-proof device
US20210016900A1 (en) * 2018-03-29 2021-01-21 Japan Aerospace Exploration Agency Method of designing a shape of an airframe of a supersonic aircraft, production method of a supersonic aircraft, and supersonic aircraft

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB719771A (en) * 1950-01-18 1954-12-08 British Electricity Authority Improvements in or relating to tele-communication equipment
US20030112107A1 (en) * 2001-12-04 2003-06-19 Nmr Holdings No. 2 Pty Ltd Asymmetric tesseral shim coils for magnetic resonance
US20040061895A1 (en) * 2002-09-27 2004-04-01 Brother Kogyo Kabushiki Kaisha Image forming device and image forming method
CN101151504A (en) * 2005-04-08 2008-03-26 独立行政法人产业技术综合研究所 Blasting treating method
US20110168004A1 (en) * 2009-10-20 2011-07-14 Henegar Douglas W System and method for mitigating and directing an explosion aboard an aircraft
CN204713434U (en) * 2015-04-02 2015-10-21 中国人民解放军空军勤务学院 Aerial bomb ammunition carrier
CN106197184A (en) * 2016-08-01 2016-12-07 西北工业大学 A kind of civil aircraft main cabin portable directional explosion-proof device
RU2649999C1 (en) * 2017-04-17 2018-04-06 Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом") Method of estimation of fougasseness characteristics in air explosion of a moving test object (variants)
US20210016900A1 (en) * 2018-03-29 2021-01-21 Japan Aerospace Exploration Agency Method of designing a shape of an airframe of a supersonic aircraft, production method of a supersonic aircraft, and supersonic aircraft
CN109581381A (en) * 2018-11-28 2019-04-05 中国民航大学 Enhanced turbulent flow detection method based on the vertical load factor
CN109780956A (en) * 2019-03-06 2019-05-21 西安近代化学研究所 A kind of cumulative pressure release civil aircraft main cabin directional explosion-proof device

Also Published As

Publication number Publication date
CN112797856B (en) 2023-01-03

Similar Documents

Publication Publication Date Title
CN104568613B (en) A kind of vessel cabin detonation inside Structure Damage evaluation method based on equivalent scale model
Buzyurkin et al. Determination and verification of Johnson–Cook model parameters at high-speed deformation of titanium alloys
Guzas et al. Air blast load generation for simulating structural response
Schonberg et al. Penetration and ricochet phenomena in oblique hypervelocity impact
Jenq et al. Simulation of a rigid plate hit by a cylindrical hemi-spherical tip-ended soft impactor
Fu et al. Scaling the response of armor steel subjected to blast loading
Hedayati et al. Foam-core effect on the integrity of tailplane leading edge during bird-strike event
Wang et al. Experimental and numerical study of tire debris impact on fuel tank cover based on coupled Eulerian-Lagrangian method
CN112797856B (en) Method for rapidly evaluating position load of minimum risk bomb of transport aircraft
Zhang et al. Study on ballistic penetration resistance of titanium alloy TC4, Part II: Numerical analysis
Yu et al. Analytical determination of the critical impact location for wing leading edge under birdstrike
Izadifard et al. Blastwave parameters assessment at different altitude using numerical simulation
Zhu et al. Axial distribution of fragments from the dynamic explosion fragmentation of metal shells
Weber et al. Numerical simulation of a 100-ton ANFO detonation
Konokman et al. Computational and experimental study of high-speed impact of metallic Taylor cylinders
Zhang et al. Application of numerical methods for crashworthiness investigation of a large aircraft wing impact with a tree
Barnat Experimental and numerical study of influence of incidence angle of shock wave created by explosive charge on the steel plate
Warsiyanto Dynamic response of the wing during UAV collision
Jha et al. Bird strike damage and analysis of UAV’s airframe
Hou et al. Dynamic response and damage of the fan blades during UAV ingestion into an aero-engine
Masi Blast actions from high explosives. Studies on their simulation and effects
Guida et al. Evaluation and validation of multi-physics FE methods to simulate bird strike on a wing leading edge
Wu et al. Safety assessment of aircraft panel under the impact load by tire fragment based on thermal–mechanical effect
Grujicic et al. Smart-fuze design and development based on computational analysis of warhead/urban-target interactions
Wang et al. An integral numerical analysis of impact of a commercial aircraft on nuclear containment

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant