CN105572720B - Flying area atmospheric neutron fluence rate multi-model composite computing method and device - Google Patents
Flying area atmospheric neutron fluence rate multi-model composite computing method and device Download PDFInfo
- Publication number
- CN105572720B CN105572720B CN201410550850.9A CN201410550850A CN105572720B CN 105572720 B CN105572720 B CN 105572720B CN 201410550850 A CN201410550850 A CN 201410550850A CN 105572720 B CN105572720 B CN 105572720B
- Authority
- CN
- China
- Prior art keywords
- fluence rate
- task point
- neutron fluence
- rate
- value
- 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.)
- Active
Links
Landscapes
- Measurement Of Radiation (AREA)
Abstract
The present invention provides a kind of flying area atmospheric neutron fluence rate multi-model composite computing method and devices, including：The mission profile and energy threshold of current electronic device are obtained, the mission profile includes current electronic device in the flying height of each task point, latitude, longitude；The first atmospheric neutron fluence rate of each task point is calculated according to the mission profile of current electronic device and energy threshold；According to the energy threshold of the current electronic device, the second air neutron fluence rate of each task point is calculated using Boeing model, and the third air neutron fluence rate of each task point is calculated using NASA models；Calculate separately the first atmospheric neutron fluence rate of each task point, the mean value of the second air neutron fluence rate and third air neutron fluence rate, the atmospheric neutron fluence rate as the task point.The present invention can more be precisely calculated atmospheric neutron fluence rate, and important evidence is provided for the protection and evaluation of air environment.
Description
Technical field
The present invention relates to microelectronics technology more particularly to a kind of flying area atmospheric neutron fluence rate multi-model are compound
Computational methods and device.
Background technology
Currently, in the world there are mainly two types of obtain atmospheric neutron fluence rate method, one be Boeing exploitation wave
Sound model, the other is NASA-Langley models, abbreviation NASA models.
Value actually based on the 1-10MeV atmospheric neutron fluence rates that the sixties in the world carries out, Boeing develop just
The simplification Boeing model of beginning.Model hypothesis 1-10MeV atmospheric neutrons fluence rate can be decomposed into three factors, a factor
Change with the variation of height, a factor changes with the variation of latitude, and a factor considers the energy of neutron.
A.Taber et al. has shown that 1-10MeV Boeing empirical models are used for predicting not by being fitted to a large amount of detection datas
With the atmospheric neutron flux of height and latitude.Its specific formula is as follows：
DN/dE=0.3459E-0.9219×exp[-0.01522·(lnE)2]
N (E) dE=26E-1.16±0.2×exp[-(0.0069x)]·dE
φ1-10(ωLat)=0.6252exp { -0.461 [cos (2 × ωLat)]2-0.94cos(2×ωLat)+0.252}
In formula：N, φ is neutron flux；
E is neutron energy；
X is atmospheric depth, g/cm2；
ω Lat are latitude.
U.S. NASA has researched and proposed NASA-Langley models also by atmospheric neutron distribution situation.NASA models
It is referred to as AIR models, is developed using the flight survey data of the sixties to the seventies.In NASA models, atmospheric neutron
Flux is influenced by three major parameters, respectively：Atmospheric density (g/cm2), vertical cut rigidity (GV) and sun environmental condition.
NASA models are a prediction more accurate empirical models of 1-10MeV atmospheric neutrons flux, which counts substantially
It is as follows to calculate formula：
Wherein：
F (R, C)=(Λ/λ) f (R, C) exp (xm/Λ-xm/λ)；
F (R, C)=exp (250/ λ) φ (250, R, C)
λ=165+2R；
xm=50+ln { 2000+exp [- 2 (C-100)] }；
In formula, φ 1-10 (x, R, C) are the atmospheric neutron flux of 1-10MeV；
X is atmospheric density, g/cm2；R is cut-off rigidity；
C is solar activity constant；
λ, xm and Λ are the intermediate parameters calculated.
Wherein, air thickness is related to height A, and A units are feet：
The factor for influencing atmospheric neutron fluence rate mainly has 4：Atmospheric density, cut-off rigidity, energy range, the sun are lived
It is dynamic.By highly determining, cut-off rigidity determines wherein atmospheric density for longitude and latitude.
Boeing model mainly considers influence of the height with latitude to atmospheric neutron fluence rate, has ignored the influence of longitude, because
This is inaccurate using the atmospheric neutron fluence rate that Boeing model obtains.
Compared with Boeing model, NASA models are more accurate.But NASA models are only applicable to be less than 20km height.And by
In the influence of parameter, it is desirable that atmospheric density is less than 250g/cm2, not applicable and ground level.
Invention content
(1) technical problems to be solved
The present invention proposes a kind of flying area atmospheric neutron fluence rate multi-model composite computing method and device, for solving
Certainly calculating atmospheric neutron fluence rate is inaccurate in Boeing model and NASA models are only applicable to less than 20km height and on ground
The not applicable problem of face height, achievees the purpose that more to accurately calculate atmospheric neutron fluence rate, is the protection of air environment
Important evidence is provided with evaluation.
(2) technical solution
In order to solve the above technical problems, the present invention provides a kind of flying area atmospheric neutron fluence rate multi-model composite computing
Method includes the following steps：
The mission profile and energy threshold of current electronic device are obtained, the mission profile includes current electronic device every
The flying height of a task point, latitude, longitude；
The first atmospheric neutron fluence of each task point is calculated according to the mission profile of current electronic device and energy threshold
Rate；
According to the energy threshold of the current electronic device, in the second air that each task point is calculated using Boeing model
Sub- fluence rate calculates the third air neutron fluence rate of each task point using NASA models；
In the first atmospheric neutron fluence rate, the second air neutron fluence rate and third air that calculate separately each task point
The mean value of sub- fluence rate, the atmospheric neutron fluence rate as the task point.
Wherein, the energy threshold for obtaining current electronic device, including：
Obtain the energy threshold of crucial Sensitive Apparatus in current electronic device, the energy threshold as current electronic device.
Wherein, the mission profile and energy threshold according to current electronic device calculates the first air of each task point
Neutron fluence rate specifically includes：
According to global longitude and latitude table corresponding with cut-off rigidity, the longitude and latitude where each task point in the mission profile is searched
Degree corresponding cut-off rigidity under 12.2km altitudes；
According to the current cut-off rigidity of each task point, using the fluence rate multiple table of cut-off rigidity and 12.2km height,
Search each task point first fluence rate multiple of fluence rate with respect to New York ground fluence rate in 12.2km height；
Each task point is searched in flying height using height latitude multiple table according to the first fluence rate multiple
Fluence rate with respect to New York ground fluence rate the second fluence rate multiple；
Fluence rate of the energy of each task point more than 10MeV when is calculated according to the second fluence rate multiple；
According to the energy threshold of the current electronic device, energy threshold amendment is carried out using energy multiple table, is obtained every
First atmospheric neutron fluence rate of a task point.
Wherein, described that fluence rate of the energy of each task point more than 10MeV when, tool are calculated according to the second fluence rate multiple
Body includes：
Obtain standard fluence rate when New York ground energy is more than 10MeV；
The energy that each task point is calculated according to the standard fluence rate and the second fluence rate multiple is more than
Fluence rate when 10MeV.
Wherein, the first atmospheric neutron fluence rate for calculating separately each task point, the second air neutron fluence rate and
The mean value of third air neutron fluence rate, specifically includes：
Extract most bad value, representative value and the average value of the first atmospheric neutron fluence rate of each task point；
Extract most bad value, representative value and the average value of the second air neutron fluence rate of each task point；
Extract most bad value, representative value and the average value of the third air neutron fluence rate of each task point；
Mean value, the mean value of the mean value of representative value and average value for calculating separately the most bad value of each task point, as this
The atmospheric neutron fluence rate of task point.
The invention also provides a kind of flying area atmospheric neutron fluence rate multi-model composite computing device, described device packets
It includes：
Acquisition module, mission profile and energy threshold for obtaining current electronic device, the mission profile include working as
Preceding electronic equipment is in the flying height of each task point, latitude, longitude；
First computing module, for calculating each task point according to the mission profile and energy threshold of current electronic device
First atmospheric neutron fluence rate；
Second computing module is calculated each for the energy threshold according to the current electronic device using Boeing model
Second air neutron fluence rate of task point calculates the third air neutron fluence rate of each task point using NASA models；
Mean value computation module, the first atmospheric neutron fluence rate, the second atmospheric neutron for calculating separately each task point
The mean value of fluence rate and third air neutron fluence rate, the atmospheric neutron fluence rate as the task point.
Wherein, the energy threshold of the current electronic device is the energy cut-off of crucial Sensitive Apparatus in current electronic device
Value.
Wherein, first computing module specifically includes：
Cut-off rigidity searching unit, for according to global longitude and latitude table corresponding with cut-off rigidity, searching the mission profile
In longitude and latitude where each task point corresponding cut-off rigidity under 12.2km altitudes；
First fluence rate multiple searching unit, for according to each task point it is current by rigidity, utilize cut-off rigidity
With the fluence rate multiple table of 12.2km height, the opposite New York ground fluence of each task point fluence rate in 12.2km height is searched
First fluence rate multiple of rate；
Second fluence rate multiple searching unit, for according to the first fluence rate multiple, using height latitude multiple table,
Search second fluence rate multiple of fluence rate of each task point in flying height with respect to New York ground fluence rate；
Condition computing unit, when for calculating the energy of each task point according to the second fluence rate multiple more than 10MeV
Fluence rate；
Energy threshold amending unit, for according to the energy threshold of the current electronic device, using energy multiple table into
Row energy threshold amendment obtains the first atmospheric neutron fluence rate of each task point.
Wherein, the condition computing unit specifically includes：
Subelement is obtained, for obtaining standard fluence rate of the New York ground energy more than 10MeV when；
Computation subunit, for calculating each task according to the standard fluence rate and the second fluence rate multiple
The energy of point is more than fluence rate when 10MeV.
Wherein, the mean value computation module specifically includes：
First extraction unit, the most bad value of the first atmospheric neutron fluence rate for extracting each task point, representative value with
Average value；
Second extraction unit, the most bad value of the second air neutron fluence rate for extracting each task point, representative value with
Average value；
Third extraction unit, the most bad value of the third air neutron fluence rate for extracting each task point, representative value with
Average value；
Average calculation unit, the mean value of the most bad value for calculating separately each task point, the mean value of representative value and flat
The mean value of mean value, the atmospheric neutron fluence rate as the task point.
(3) advantageous effect
Using flying area atmospheric neutron fluence rate multi-model composite computing method and device proposed by the present invention, specify
An important factor for cut-off rigidity is influence atmospheric neutron fluence rate, rather than single latitude, correct the latitude shadow in Boeing model
It is longitude, latitude to ring, and solves the problems, such as that NASA models are not applicable in ground level, by quantifying each different-energy threshold value pair
The influence of atmospheric neutron fluence rate is more precisely calculated atmospheric neutron fluence rate of the electronic device at task point.
Description of the drawings
The features and advantages of the present invention can be more clearly understood by reference to attached drawing, attached drawing is schematically without that should manage
Solution is carries out any restrictions to the present invention, in the accompanying drawings：
Fig. 1 is a kind of flying area atmospheric neutron fluence rate multi-model composite computing method that the embodiment of the present invention one proposes
Flow chart；
Fig. 2 is a kind of flying area atmospheric neutron fluence rate multi-model composite computing device that the embodiment of the present invention two proposes
Module map.
Specific implementation mode
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
The every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.
Air environment with storage organization complexity microelectronic component flying height (3000~20000 meters) from
The about high energy air of 300~18000 1MeV~1000MeV every square centimeter per hour will necessarily be met in right space environment
Neutron.These high-energy neutrons can penetrate cabin skin, beat kernel instruction control unit or critical data in air environment
In storage unit, soft error and hard fault are generated, leads to navigate (navigation neceiver), radar-probing system (active phased array thunder
Up to), data network (the AFDX network switch), communication (optical fiber/bus), high speed computer system, avionic device, start
Machine (FADEC), telex system, automatic Pilot technology, flight alarm, display screen, other flight systems etc. containing electronic device
There is blank screen, crash, resets, restarts, loss of data, ordering the safety hazards such as loss.In order to establish protection and appraisement system,
Atmospheric neutron single particle effect must be endangered and carry out quantitatively characterizing.And atmospheric neutron fluence rate is endanger quantitatively characterizing important
Parameter.In this regard, the present invention proposes a kind of method that can obtain atmospheric neutron fluence rate.
Fig. 1 is a kind of flying area atmospheric neutron fluence rate multi-model composite computing method that the embodiment of the present invention one proposes
Flow chart, as shown in Figure 1, including the following steps：
S101 obtains the mission profile and energy threshold of current electronic device, and the mission profile includes current electronic device
In the flying height of each task point, latitude, longitude.
In the present embodiment, the energy threshold of current electronic device is obtained, including：It obtains crucial sensitive in current electronic device
The energy threshold of device, the energy threshold as current electronic device.
The mission profile proposed in the embodiment of the present invention be the operation flight height of current electronic device, latitude, longitude and
The parameters such as flight beginning and ending time.Factor due to influencing atmospheric neutron fluence rate includes cut-off rigidity, and cut-off rigidity is by longitude
It is codetermined with latitude, therefore, the present invention is two longitude, latitude parameters, energy by the latitude parameter optimized in Boeing model
Enough improve the computational accuracy of atmospheric neutron fluence rate.
S102 calculates the first atmospheric neutron of each task point according to the mission profile and energy threshold of current electronic device
Fluence rate；
S103 calculates the second largest of each task point according to the energy threshold of the current electronic device using Boeing model
Gas neutron fluence rate calculates the third air neutron fluence rate of each task point using NASA models；
S104 calculates separately the first atmospheric neutron fluence rate of each task point, the second air neutron fluence rate and the third-largest
The mean value of gas neutron fluence rate, the atmospheric neutron fluence rate as the task point.
In embodiments of the present invention, each task point is calculated according to the mission profile of current electronic device and energy threshold
First atmospheric neutron fluence rate, specifically includes：According to global longitude and latitude table corresponding with cut-off rigidity, as shown in table 1, described in lookup
Longitude and latitude in mission profile where each task point corresponding cut-off rigidity under 12.2km altitudes；According to each task
Point it is current by rigidity, using the fluence rate multiple table of cut-off rigidity and 12.2km height, as shown in table 2, search each
The business point first fluence rate multiple of fluence rate with respect to New York ground fluence rate in 12.2km height；According to first fluence rate
Multiple as shown in table 3, searches fluence rate of each task point in flying height with respect to New York using height latitude multiple table
Second fluence rate multiple of ground fluence rate；When calculating the energy of each task point more than 10MeV according to the second fluence rate multiple
Fluence rate；According to the energy threshold of the current electronic device, using energy multiple table, as shown in table 4, energy threshold is carried out
Correct, obtain the first atmospheric neutron fluence rate of each task point, including most bad value, representative value (12200,45 degree of height) with put down
Mean value.
The global longitude and latitude of table 1 table (1) corresponding with cut-off rigidity
The global longitude and latitude of table 1 table (2) corresponding with cut-off rigidity
In global longitude and latitude table corresponding with cut-off rigidity, north latitude is just that south latitude is negative；Longitude is expressed as east longitude 0-
360 °, as 360-x ° of W.Such as 30S, 120W, as -30N, 240E.
The fluence rate multiple table of 2 cut-off rigidity of table and 12.2km height
| Cut-off rigidity | Multiple | Cut-off rigidity | Multiple | Cut-off rigidity | Multiple | Cut-off rigidity | Multiple |
| 0.00 | 561.70 | 2.20 | 506.89 | 8.38 | 193.2 | 13.43 | 111.25 |
| 0.01 | 561.70 | 2.30 | 503.85 | 8.99 | 174.36 | 13.50 | 110.34 |
| 0.02 | 561.70 | 2.32 | 497.65 | 9.23 | 167.52 | 13.57 | 109.77 |
| 0.03 | 561.70 | 2.40 | 493.72 | 9.56 | 163.33 | 13.62 | 109.49 |
| 0.05 | 561.70 | 2.56 | 483.35 | 9.57 | 162.04 | 13.76 | 108.20 |
| 0.06 | 561.70 | 2.73 | 463.07 | 9.66 | 161.49 | 13.77 | 108.10 |
| 0.08 | 561.70 | 2.75 | 463.07 | 10.17 | 153.78 | 13.80 | 107.74 |
| 0.09 | 561.70 | 2.85 | 454.95 | 10.55 | 146.52 | 13.88 | 107.01 |
| 0.13 | 561.70 | 3.01 | 447.68 | 10.81 | 142.74 | 13.92 | 106.75 |
| 0.14 | 561.70 | 3.05 | 441.25 | 10.97 | 140.69 | 14.06 | 105.50 |
| 0.17 | 561.69 | 3.20 | 429.35 | 10.98 | 140.12 | 14.10 | 105.15 |
| 0.19 | 561.70 | 3.30 | 419.22 | 11.10 | 138.42 | 14.19 | 104.37 |
| 0.20 | 561.70 | 3.56 | 399.60 | 11.23 | 132.08 | 14.20 | 104.20 |
| 0.24 | 561.70 | 3.68 | 385.15 | 11.25 | 139.55 | 14.35 | 103.09 |
| 0.25 | 561.70 | 3.71 | 387.28 | 11.32 | 134.72 | 14.37 | 102.84 |
| 0.27 | 561.70 | 4.18 | 356.77 | 11.59 | 131.82 | 14.46 | 102.01 |
| 0.28 | 561.70 | 4.28 | 346.05 | 11.61 | 131.57 | 14.64 | 100.62 |
| 0.36 | 561.70 | 4.30 | 360.00 | 11.66 | 130.16 | 14.65 | 100.54 |
| 0.38 | 561.70 | 4.33 | 347.28 | 11.67 | 130.16 | 14.71 | 100.05 |
| 0.42 | 561.70 | 4.55 | 333.40 | 11.72 | 131.18 | 14.94 | 98.24 |
| 0.52 | 561.70 | 4.71 | 329.90 | 12.22 | 123.64 | 14.95 | 98.16 |
| 0.55 | 561.70 | 4.90 | 313.20 | 12.35 | 122.38 | 15.16 | 96.56 |
| 0.68 | 561.69 | 5.02 | 311.59 | 12.36 | 122.50 | 15.37 | 95.00 |
| 0.74 | 561.67 | 5.37 | 290.39 | 12.39 | 121.59 | 16.00 | 90.21 |
| 0.78 | 561.63 | 5.43 | 287.54 | 12.67 | 119.06 | 17.00 | 82.41 |
| 0.79 | 561.67 | 5.58 | 280.60 | 12.72 | 119.06 | ||
| 0.85 | 561.50 | 5.83 | 272.62 | 12.73 | 118.84 | ||
| 0.89 | 561.48 | 6.04 | 260.48 | 12.87 | 116.93 | ||
| 0.93 | 561.39 | 6.10 | 262.51 | 12.88 | 116.82 | ||
| 0.99 | 560.73 | 6.63 | 239.35 | 12.91 | 116.30 | ||
| 1.08 | 560.16 | 6.78 | 233.48 | 12.99 | 115.68 |
| Cut-off rigidity | Multiple | Cut-off rigidity | Multiple | Cut-off rigidity | Multiple | Cut-off rigidity | Multiple |
| 1.10 | 559.93 | 6.86 | 232.13 | 13.02 | 115.37 | ||
| 1.15 | 560.04 | 7.13 | 222.10 | 13.06 | 114.96 | ||
| 1.19 | 558.02 | 7.17 | 230.46 | 13.07 | 114.66 | ||
| 1.44 | 551.49 | 7.41 | 213.65 | 13.20 | 113.35 | ||
| 1.69 | 540.63 | 7.57 | 211.90 | 13.22 | 113.35 | ||
| 1.85 | 530.70 | 7.73 | 204.62 | 13.26 | 112.96 | ||
| 1.86 | 534.99 | 7.89 | 196.71 | 13.31 | 112.66 | ||
| 2.19 | 514.32 | 8.03 | 195.44 | 13.38 | 111.74 |
3 height latitude multiple table (1) of table
3 each high latitude of table is with respect to New York ground fluence rate multiple table (2)
4 energy threshold multiple table of table
| Serial number | Energy (MeV) | Neutron fluence rate (#/cm2h) | With 1-10MeV multiple proportions |
| 1 | >1 | 9200 | 2.88 |
| 2 | >2 | 8300 | 2.59 |
| 3 | >3 | 7700 | 2.41 |
| 4 | >4 | 7300 | 2.28 |
| 5 | >5 | 7000 | 2.19 |
| 6 | >6 | 6800 | 2.13 |
| 7 | >7 | 6500 | 2.03 |
| 8 | >8 | 6400 | 2.00 |
| 9 | >9 | 6200 | 1.94 |
| 10 | >10 | 6000 | 1.88 |
| 11 | 1-10 | 3200 | 1.00 |
The energy that each task point is calculated according to the second fluence rate multiple proposed in embodiments of the present invention is more than 10MeV
When fluence rate, specifically include：
Obtain standard fluence rate when New York ground energy is more than 10MeV；According to the standard fluence rate and described second
Fluence rate multiple calculates fluence rate of the energy of each task point more than 10MeV when.
In embodiments of the present invention, step S103 is specifically included：The energy cut-off of each task point is calculated using Boeing model
It is worth the above atmospheric neutron fluence rate, obtains most bad value, representative value (12200,45 degree of height) and average value.And use NASA models
Calculate the above atmospheric neutron fluence rate of energy threshold of each task point.Find out most bad value, representative value (12200,45 degree of height)
With average value.
In the embodiment of the present invention, step S104 is specifically included：
Extract most bad value, representative value and the average value of the first atmospheric neutron fluence rate of each task point；
Extract most bad value, representative value and the average value of the second air neutron fluence rate of each task point；
Extract most bad value, representative value and the average value of the third air neutron fluence rate of each task point；
Mean value, the mean value of the mean value of representative value and average value for calculating separately the most bad value of each task point, as this
The atmospheric neutron fluence rate of task point.
Flying area atmospheric neutron fluence rate multi-model composite computing method provided by the invention, passes through the warp of task point
Degree, latitude determine cut-off rigidity, and the latitude influence in Boeing model is modified to the cut-off rigidity determined by longitude, latitude, and
And it is suitable for any height, solve the problems, such as that NASA models are not applicable in ground level, by quantifying each different-energy threshold
It is worth the influence to atmospheric neutron fluence rate, is more precisely calculated atmospheric neutron fluence rate of the electronic device at task point.
Fig. 2 is a kind of flying area atmospheric neutron fluence rate multi-model composite computing device that the embodiment of the present invention two proposes
Module map, as shown in Fig. 2, described device includes：
Acquisition module 201, mission profile and energy threshold for obtaining current electronic device, the mission profile include
Current electronic device is in the flying height of each task point, latitude, longitude.In the present embodiment, the energy of the current electronic device
Measure the energy threshold that threshold value is crucial Sensitive Apparatus in current electronic device.
First computing module 202, for calculating each task according to the mission profile and energy threshold of current electronic device
First atmospheric neutron fluence rate of point；
Second computing module 203 is calculated every for the energy threshold according to the current electronic device using Boeing model
Second air neutron fluence rate of a task point calculates the third air neutron fluence rate of each task point using NASA models；
Mean value computation module 204, for calculating separately in the first atmospheric neutron fluence rate of each task point, the second air
The mean value of sub- fluence rate and third air neutron fluence rate, the atmospheric neutron fluence rate as the task point.
The first computing module 202 in the present embodiment specifically includes：
Cut-off rigidity searching unit, for according to global longitude and latitude table corresponding with cut-off rigidity, searching the mission profile
In longitude and latitude where each task point corresponding cut-off rigidity under 12.2km altitudes；
First fluence rate multiple searching unit, for according to each task point it is current by rigidity, utilize cut-off rigidity
With the fluence rate multiple table of 12.2km height, the opposite New York ground fluence of each task point fluence rate in 12.2km height is searched
First fluence rate multiple of rate；
Second fluence rate multiple searching unit, for according to the first fluence rate multiple, using height latitude multiple table,
Search second fluence rate multiple of fluence rate of each task point in flying height with respect to New York ground fluence rate；
Condition computing unit, when for calculating the energy of each task point according to the second fluence rate multiple more than 10MeV
Fluence rate；
Energy threshold amending unit, for according to the energy threshold of the current electronic device, using energy multiple table into
Row energy threshold amendment obtains the first atmospheric neutron fluence rate of each task point.
Condition computing unit in the present embodiment specifically includes：
Subelement is obtained, for obtaining standard fluence rate of the New York ground energy more than 10MeV when；
Computation subunit, for calculating each task according to the standard fluence rate and the second fluence rate multiple
The energy of point is more than fluence rate when 10MeV.
Mean value computation module 204 in the present embodiment specifically includes：
First extraction unit, the most bad value of the first atmospheric neutron fluence rate for extracting each task point, representative value with
Average value；
Second extraction unit, the most bad value of the second air neutron fluence rate for extracting each task point, representative value with
Average value；
Third extraction unit, the most bad value of the third air neutron fluence rate for extracting each task point, representative value with
Average value；
Average calculation unit, the mean value of the most bad value for calculating separately each task point, the mean value of representative value and flat
The mean value of mean value, the atmospheric neutron fluence rate as the task point.
Flying area atmospheric neutron fluence rate multi-model composite computing method and device proposed by the present invention, specifies cut-off
An important factor for rigidity is influence atmospheric neutron fluence rate, rather than single latitude, the latitude parameter optimized in Boeing model are
Two longitude, latitude parameters, this method can more accurately calculate atmospheric neutron fluence rate.
The present invention solves the problems, such as that NASA models are not applicable in ground level, is applicable in and various height.
The influence of energy range is added in the present invention.Boeing model is mainly calculated with NASA in the air that energy is 1-10MeV
Sub- fluence rate, and in view of the energy threshold of different components is different, the threshold values of some devices may be 4MeV, 5MeV, and the party
Method can be more than 1MeV with computing capability, be more than the different-energies models such as 2MeV to more than 10MeV by the influence of addition energy range
The neutron fluence rate enclosed keeps result of calculation more accurate.
The present invention improves computational accuracy, user can choose most suitably used data as task by rule of thumb by calculating mean value
The atmospheric neutron fluence rate of point.
Through the above description of the embodiments, those skilled in the art can be understood that the present invention can lead to
Hardware realization is crossed, the mode of necessary general hardware platform can also be added to realize by software.Based on this understanding, this hair
Bright technical solution can be expressed in the form of software products, which can be stored in a non-volatile memories
In medium (can be CD-ROM, USB flash disk, mobile hard disk etc.), including some instructions are used so that a computer equipment (can be
Personal computer, server or network equipment etc.) execute method described in each embodiment of the present invention.
It will be appreciated by those skilled in the art that attached drawing is the schematic diagram of a preferred embodiment, the module in attached drawing or stream
Journey is not necessarily implemented necessary to the present invention.
It will be appreciated by those skilled in the art that the module in device in embodiment can describe be divided according to embodiment
It is distributed in the device of embodiment, respective change can also be carried out and be located in one or more devices different from the present embodiment.On
The module for stating embodiment can be merged into a module, can also be further split into multiple submodule.
Disclosed above is only several specific embodiments of the present invention, and still, the present invention is not limited to this, any ability
What the technical staff in domain can think variation should all fall into protection scope of the present invention.
Claims (8)
1. a kind of flying area atmospheric neutron fluence rate multi-model composite computing method, which is characterized in that including：
The mission profile and energy threshold of current electronic device are obtained, the mission profile includes current electronic device at each
The flying height of business point, latitude, longitude；
The first atmospheric neutron fluence rate of each task point is calculated according to the mission profile of current electronic device and energy threshold；
According to global longitude and latitude table corresponding with cut-off rigidity, searches the longitude and latitude in the mission profile where each task point and exist
Corresponding cut-off rigidity under 12.2km altitudes；
It is searched using the fluence rate multiple table of cut-off rigidity and 12.2km height according to the current cut-off rigidity of each task point
The each task point first fluence rate multiple of fluence rate with respect to New York ground fluence rate in 12.2km height；
According to the first fluence rate multiple note of each task point in flying height is searched using height latitude multiple table
Second fluence rate multiple of the dose rate with respect to New York ground fluence rate；
Fluence rate of the energy of each task point more than 10MeV when is calculated according to the second fluence rate multiple；
According to the energy threshold of the current electronic device, energy threshold amendment is carried out using energy multiple table, obtains each
First atmospheric neutron fluence rate of business point；
According to the energy threshold of the current electronic device, the second atmospheric neutron that each task point is calculated using Boeing model is noted
Dose rate calculates the third air neutron fluence rate of each task point using NASA models；
Calculate separately the first atmospheric neutron fluence rate, the second air neutron fluence rate and third the atmospheric neutron note of each task point
The mean value of dose rate, the atmospheric neutron fluence rate as the task point.
2. according to the method described in claim 1, it is characterized in that, it is described obtain current electronic device energy threshold, including：
Obtain the energy threshold of crucial Sensitive Apparatus in current electronic device, the energy threshold as current electronic device.
3. the method stated according to claim 1, which is characterized in that described to calculate each task point according to the second fluence rate multiple
Energy is more than fluence rate when 10MeV, specifically includes：
Obtain standard fluence rate when New York ground energy is more than 10MeV；
When calculating the energy of each task point more than 10MeV according to the standard fluence rate and the second fluence rate multiple
Fluence rate.
4. according to the method described in claim 1, it is characterized in that, first atmospheric neutron for calculating separately each task point
The mean value of fluence rate, the second air neutron fluence rate and third air neutron fluence rate, specifically includes：
Extract most bad value, representative value and the average value of the first atmospheric neutron fluence rate of each task point；
Extract most bad value, representative value and the average value of the second air neutron fluence rate of each task point；
Extract most bad value, representative value and the average value of the third air neutron fluence rate of each task point；
Mean value, the mean value of the mean value of representative value and average value for calculating separately the most bad value of each task point, as the task
The atmospheric neutron fluence rate of point.
5. a kind of flying area atmospheric neutron fluence rate multi-model composite computing device, which is characterized in that described device includes：
Acquisition module, mission profile and energy threshold for obtaining current electronic device, the mission profile include current electricity
Sub- equipment is in the flying height of each task point, latitude, longitude；
First computing module, for calculating the first of each task point according to the mission profile and energy threshold of current electronic device
Atmospheric neutron fluence rate；
First computing module specifically includes：
Cut-off rigidity searching unit, it is every in the mission profile for according to global longitude and latitude table corresponding with cut-off rigidity, searching
Longitude and latitude where a task point corresponding cut-off rigidity under 12.2km altitudes；
First fluence rate multiple searching unit, for according to each task point it is current by rigidity, using cut-off rigidity with
The fluence rate multiple table of 12.2km height searches the opposite New York ground fluence rate of each task point fluence rate in 12.2km height
The first fluence rate multiple；
Second fluence rate multiple searching unit, for being searched using height latitude multiple table according to the first fluence rate multiple
Second fluence rate multiple of fluence rate of each task point in flying height with respect to New York ground fluence rate；
Condition computing unit, the energy for calculating each task point according to the second fluence rate multiple are more than fluence when 10MeV
Rate；
Energy threshold amending unit carries out energy for the energy threshold according to the current electronic device using energy multiple table
Threshold value amendment is measured, the first atmospheric neutron fluence rate of each task point is obtained；
Second computing module calculates each task for the energy threshold according to the current electronic device using Boeing model
Second air neutron fluence rate of point, the third air neutron fluence rate of each task point is calculated using NASA models；
Mean value computation module, the first atmospheric neutron fluence rate, the second air neutron fluence for calculating separately each task point
The mean value of rate and third air neutron fluence rate, the atmospheric neutron fluence rate as the task point.
6. device according to claim 5, which is characterized in that the energy threshold of the current electronic device is Current electronic
The energy threshold of crucial Sensitive Apparatus in equipment.
7. device according to claim 5, which is characterized in that the condition computing unit specifically includes：
Subelement is obtained, for obtaining standard fluence rate of the New York ground energy more than 10MeV when；
Computation subunit, for calculating each task point according to the standard fluence rate and the second fluence rate multiple
Energy is more than fluence rate when 10MeV.
8. device according to claim 5, which is characterized in that the mean value computation module specifically includes：
First extraction unit, the most bad value of the first atmospheric neutron fluence rate for extracting each task point, representative value with it is average
Value；
Second extraction unit, the most bad value of the second air neutron fluence rate for extracting each task point, representative value with it is average
Value；
Third extraction unit, the most bad value of the third air neutron fluence rate for extracting each task point, representative value with it is average
Value；
Average calculation unit, mean value, the mean value of representative value and the average value of the most bad value for calculating separately each task point
Mean value, the atmospheric neutron fluence rate as the task point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410550850.9A CN105572720B (en) | 2014-10-16 | 2014-10-16 | Flying area atmospheric neutron fluence rate multi-model composite computing method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410550850.9A CN105572720B (en) | 2014-10-16 | 2014-10-16 | Flying area atmospheric neutron fluence rate multi-model composite computing method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105572720A CN105572720A (en) | 2016-05-11 |
| CN105572720B true CN105572720B (en) | 2018-08-03 |
Family
ID=55883061
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410550850.9A Active CN105572720B (en) | 2014-10-16 | 2014-10-16 | Flying area atmospheric neutron fluence rate multi-model composite computing method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105572720B (en) |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2910638B1 (en) * | 2006-12-21 | 2009-02-13 | Thales Sa | METHOD FOR SIMULATING FAILURE RATE OF AN ELECTRONIC EQUIPMENT DUE TO NEUTRON RADIATION |
| CN102540241A (en) * | 2010-12-17 | 2012-07-04 | 中国辐射防护研究院 | Method for measuring neutron dose equivalent by utilizing semiconductor array |
| CN202584750U (en) * | 2012-02-17 | 2012-12-05 | 中国原子能科学研究院 | Neutron fluence automatic measuring system for physical start of reactor |
-
2014
- 2014-10-16 CN CN201410550850.9A patent/CN105572720B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105572720A (en) | 2016-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103675834B (en) | A kind of indoor satellite signal simulation system | |
| CN111538043B (en) | Method and device for eliminating non-through-view satellite in complex environment and electronic equipment | |
| CN110606212A (en) | Approach landing method, device, equipment and storage medium | |
| Ercha et al. | Ionospheric TEC data assimilation and now-casting system over China | |
| CN105676103B (en) | Obtain the method and device in atmospheric neutron single particle effect Sensitive Apparatus sensitivity section | |
| CN103440408A (en) | Rapid analysis method for landing safety probability of lunar probe | |
| KR101537013B1 (en) | Apparatus and method of track generation for global navigation satellite system simulator | |
| CN105572720B (en) | Flying area atmospheric neutron fluence rate multi-model composite computing method and device | |
| CN105572717B (en) | Avionic device atmospheric neutron fluence rate calculation method and system | |
| CN109521442A (en) | One kind being based on the quick cloth station method of satellite-based augmentation system | |
| CN104925271B (en) | The defining method reentering standard trajectory instruction of lift formula aircraft | |
| Jin et al. | Ionospheric anomalies during the March 2013 geomagnetic storm from BeiDou Navigation Satellite System (BDS) observations | |
| CN105676016B (en) | The method and device in neutron single-particle effect device sensitivity section is obtained using BGR | |
| CN105572718B (en) | Flying area atmospheric neutron fluence rate calculation method and system | |
| CN105676017B (en) | The method and device in single particle effect device sensitivity section is obtained using test data | |
| CN105572719B (en) | A kind of calculation method and system of 12.2km height atmospheric neutron fluence rate | |
| Xue et al. | A novel LAAS pseudo-range error over-bound method based on improved pseudo-range error distribution model | |
| CN104281659A (en) | GIS (geographic information system)-platform based dynamic loading method for real-time weather of graticule data | |
| CN103675845A (en) | Single-carrier multiple-antenna occultation signal generation system | |
| CN116415429B (en) | Arm length accurate estimation method for satellite formation | |
| Damaceno et al. | Regional short-term forecasting model to predict ionospheric scintillation and TEC at low latitudes | |
| Hanif et al. | Implementation of real-time kinematic data to determine the ionospheric total electron content | |
| CN107313766A (en) | Attitude data bearing calibration and device | |
| Shao et al. | A particle swarm optimization based algorithm for the calculation of user differential range error | |
| Wang et al. | Availability Analysis of GNSS RAIM Based on Maximum Undetectable Position Bias |
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 |