CN105572720A - Flight space atmosphere neutron fluence rate multi-model composite calculating method and apparatus thereof - Google Patents

Abstract

The invention provides a flight space atmosphere neutron fluence rate multi-model composite calculating method and an apparatus thereof. The method comprises the following steps of acquiring a task profile and an energy threshold value of current electronic equipment, wherein the task profile comprises a flight height, a latitude and a longitude of the current electronic equipment at each task point; according to the task profile and the energy threshold value of the current electronic equipment, calculating a first atmosphere neutron fluence rate of each task point; according to the energy threshold value of the current electronic equipment, using a Boeing model to calculate a second atmosphere neutron fluence rate of each task point and using an NASA model to calculate a third atmosphere neutron fluence rate of each task point; calculating a mean value of the first atmosphere neutron fluence rate, the second atmosphere neutron fluence rate and the third atmosphere neutron fluence rate of each task point respectively and taking the mean value as an atmosphere neutron fluence rate of the task point. By using the method and the apparatus, the atmosphere neutron fluence rate can be accurately calculated and an important basis is provided for airborne electronic equipment protection and evaluation.

Description

Flying area atmospheric neutron fluence rate multi-model composite computing method and device

Technical field

The present invention relates to microelectronics technology, particularly relate to a kind of flying area atmospheric neutron fluence rate multi-model composite computing method and device.

Background technology

At present, mainly contain the method that two kinds obtain atmospheric neutron fluence rate in the world, the Boeing model of Ge Shi Boeing exploitation, another is NASA-Langley model, is called for short NASA model.

Based on the value actually of the 1-10MeV atmospheric neutron fluence rate carried out the sixties in the world, Boeing develops initial simplification Boeing model.This model hypothesis 1-10MeV atmospheric neutron fluence rate can be decomposed into three factors, and a factor changes along with the change of height, and a factor changes along with the change of latitude, and a factor considers the energy of neutron.The people such as A.Taber, by carrying out matching to a large amount of detection data, have shown that 1-10MeVBoeing empirical model is used for predicting the atmospheric neutron flux of differing heights and latitude.Its concrete formula is as follows:

dN/dE＝0.3459E ^-0.9219×exp[-0.01522·(lnE) ²]

N(E)dE＝26E ^-1.16±0.2×exp[-(0.0069x)]·dE

φ _1-10(ω _Lat)＝0.6252exp{-0.461[cos(2×ω _Lat)] ²-0.94cos(2×ω _Lat)+0.252}

In formula: N, φ are neutron flux;

E is neutron energy;

X is atmospheric depth, g/cm2;

ω Lat is latitude.

U.S. NASA is also by having researched and proposed NASA-Langley model to atmospheric neutron distribution situation.NASA model is called as AIR model, is to utilize the exploitation of the flight survey data of the sixties to the seventies to form.In NASA model, atmospheric neutron flux affects by three major parameters, is respectively: atmospheric density (g/cm2), vertical cut rigidity (GV) and sun environmental baseline.

NASA model is prediction 1-10MeV atmospheric neutron flux empirical model more accurately, and this model basic calculating formula is as follows:

${\mathrm{\φ}}_{1-10}(x,R,C)=f(R,C)\·\exp (-\frac{x}{\mathrm{\λ}})\·F(R,C)\·\exp (-\frac{x}{\mathrm{\Λ}});$

Wherein:

F(R,C)＝(Λ/λ)·f(R,C)·exp(x _m/Λ-x _m/λ)；

f(R,C)＝exp(250/λ)φ(250,R,C)

$\begin{array}{c}\mathrm{\φ}(250,R,C)=0.17+[0.787+0.035(C-100)]\exp (-\frac{R^2}{20})\\ +[-0.107-0.0265(C-100)+0.612\exp \left(\frac{C-100}{3.73}\right)]\exp (-\frac{R^2}{139.2})\end{array},;$

λ＝165+2R；

x _m＝50+ln{2000+exp[-2(C-100)]}；

In formula, the atmospheric neutron flux that φ 1-10 (x, R, C) is 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 is thick relevant to height A, and A unit is feet:

$x=1033\exp \{-[0.04534-(1.17E-9)\×{\left|\frac{A-1.05E5}{1000}\right|}^{3.58}]\×\frac{A}{1000}\}$

The factor affecting atmospheric neutron fluence rate mainly contains 4: atmospheric density, cut-off rigidity, energy range, solar activity.Wherein atmospheric density is by highly determining, cut-off rigidity is that longitude and latitude determines.

Boeing model mainly considers height and the impact of latitude on atmospheric neutron fluence rate, have ignored the impact of longitude, therefore applies the atmospheric neutron fluence rate that Boeing model obtains accurate not.

Compared with Boeing model, NASA model is comparatively accurate.But NASA model is only applicable to lower than 20km height.And due to the impact of parameter, require that atmospheric density is lower than 250g/cm2, inapplicable and floor level.

Summary of the invention

(1) technical matters that will solve

The present invention proposes a kind of flying area atmospheric neutron fluence rate multi-model composite computing method and device, for solving in Boeing model, to calculate atmospheric neutron fluence rate accurate not, and NASA model be only applicable to lower than 20km height and in the inapplicable problem of floor level, reach the object of accurate Calculation atmospheric neutron fluence rate more, for the protection of air environment and evaluation provide important evidence.

(2) technical scheme

For solving the problems of the technologies described above, the invention provides a kind of flying area atmospheric neutron fluence rate multi-model composite computing method, comprising the following steps:

Obtain mission profile and the energy threshold of current electronic device, described mission profile comprises flying height, latitude, the longitude of current electronic device at each task point;

The first air 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 described current electronic device, adopt Boeing model to calculate the second air neutron fluence rate of each task point, adopt NASA model to calculate the third-largest gas neutron fluence rate of each task point;

Calculate the average of the first air neutron fluence rate of each task point, the second air neutron fluence rate and the third-largest gas neutron fluence rate respectively, as the atmospheric neutron fluence rate of this task point.

Wherein, the energy threshold of described acquisition current electronic device, comprising:

Obtain the energy threshold of crucial Sensitive Apparatus in current electronic device, as the energy threshold of current electronic device.

Wherein, the described mission profile according to current electronic device and energy threshold calculate the first air neutron fluence rate of each task point, specifically comprise:

Show according to global longitude and latitude is corresponding with cut-off rigidity, search the cut-off rigidity that the longitude and latitude at each task point place in described mission profile is corresponding under 12.2km altitudes;

The cut-off rigidity current according to each task point, utilizes the fluence rate multiple table of cut-off rigidity and 12.2km height, searches the first fluence rate multiple of each task point ground, fluence rate relative New York fluence rate when 12.2km height;

According to described first fluence rate multiple, utilize height latitude multiple table, search the second fluence rate multiple of ground, the fluence rate of each task point when flying height relative New York fluence rate;

The energy calculating each task point according to the second fluence rate multiple is greater than fluence rate during 10MeV;

According to the energy threshold of described current electronic device, utilize energy multiple table to carry out energy threshold correction, obtain the first air neutron fluence rate of each task point.

Wherein, the described energy calculating each task point according to the second fluence rate multiple is greater than fluence rate during 10MeV, specifically comprises:

Acquisition ground, New York energy is greater than standard fluence rate during 10MeV;

The fluence rate when energy of described each task point is greater than 10MeV is calculated according to described standard fluence rate and described second fluence rate multiple.

Wherein, the average of described the first air neutron fluence rate, the second air neutron fluence rate and the third-largest gas neutron fluence rate that calculate each task point respectively, specifically comprises:

Extract the most bad value of the first air neutron fluence rate of each task point, representative value and mean value;

Extract the most bad value of the second air neutron fluence rate of each task point, representative value and mean value;

Extract the most bad value of the third-largest gas neutron fluence rate of each task point, representative value and mean value;

Calculate the average of the average of the most bad value of each task point, the average of representative value and mean value respectively, as the atmospheric neutron fluence rate of this task point.

The invention allows for a kind of flying area atmospheric neutron fluence rate multi-model composite computing device, described device comprises:

Acquisition module, for obtaining mission profile and the energy threshold of current electronic device, described mission profile comprises flying height, latitude, the longitude of current electronic device at each task point;

First computing module, for calculating the first air neutron fluence rate of each task point according to the mission profile of current electronic device and energy threshold;

Second computing module, for the energy threshold according to described current electronic device, adopts Boeing model to calculate the second air neutron fluence rate of each task point, adopts NASA model to calculate the third-largest gas neutron fluence rate of each task point;

Mean value computation module, for calculating the average of the first air neutron fluence rate of each task point, the second air neutron fluence rate and the third-largest gas neutron fluence rate respectively, as the atmospheric neutron fluence rate of this task point.

Wherein, the energy threshold of described current electronic device is the energy threshold of crucial Sensitive Apparatus in current electronic device.

Wherein, described first computing module specifically comprises:

Cut-off rigidity searches unit, for showing according to global longitude and latitude is corresponding with cut-off rigidity, searches the cut-off rigidity that the longitude and latitude at each task point place in described mission profile is corresponding under 12.2km altitudes;

First fluence rate multiple searches unit, for according to each task point current by rigidity, utilize the fluence rate multiple table of cut-off rigidity and 12.2km height, search the first fluence rate multiple of each task point ground, fluence rate relative New York fluence rate when 12.2km height;

Second fluence rate multiple searches unit, for according to described first fluence rate multiple, utilizes height latitude multiple table, searches the second fluence rate multiple of ground, the fluence rate of each task point when flying height relative New York fluence rate;

Condition computing unit, the energy for calculating each task point according to the second fluence rate multiple is greater than fluence rate during 10MeV;

Energy threshold amending unit, for the energy threshold according to described current electronic device, utilizes energy multiple table to carry out energy threshold correction, obtains the first air neutron fluence rate of each task point.

Wherein, described condition computing unit specifically comprises:

Obtain subelement, for obtaining standard fluence rate when ground, New York energy is greater than 10MeV;

Computation subunit, the energy for calculating described each task point according to described standard fluence rate and described second fluence rate multiple is greater than fluence rate during 10MeV.

Wherein, described mean value computation module specifically comprises:

First extraction unit, for extracting the most bad value of the first air neutron fluence rate of each task point, representative value and mean value;

Second extraction unit, for extracting the most bad value of the second air neutron fluence rate of each task point, representative value and mean value;

3rd extraction unit, for extracting the most bad value of the third-largest gas neutron fluence rate of each task point, representative value and mean value;

Average calculation unit, for calculating the average of the average of the most bad value of each task point, the average of representative value and mean value respectively, as the atmospheric neutron fluence rate of this task point.

(3) beneficial effect

The flying area atmospheric neutron fluence rate multi-model composite computing method adopting the present invention to propose and device, specify that cut-off rigidity is the key factor affecting atmospheric neutron fluence rate, instead of single latitude, the latitude impact revised in Boeing model is longitude, latitude, solve NASA model in the inapplicable problem of floor level, each different-energy threshold value is on the impact of atmospheric neutron fluence rate by quantifying, more accurately calculates the atmospheric neutron fluence rate of electron device at task point place.

Accompanying drawing explanation

Can understanding the features and advantages of the present invention clearly by reference to accompanying drawing, accompanying drawing is schematic and should not be construed as and carry out any restriction 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 flow diagram that the embodiment of the present invention one proposes;

Fig. 2 is a kind of flying area atmospheric neutron fluence rate multi-model composite computing apparatus module figure that the embodiment of the present invention two proposes.

Embodiment

For making the object of the embodiment of the present invention, technical scheme and advantage clearly, below in conjunction with the accompanying drawing in the embodiment of the present invention, technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.

Air environment with the complicated microelectronic component of storage organization will inevitably meet with the high energy atmospheric neutron of about every square centimeter of 300 ~ 18000 1MeV ~ 1000MeV per hour in the place environment of flying height (3000 ~ 20000 meters).These high-energy neutrons can penetrate cabin skin, beat on the kernel instruction control module or critical data storage unit of air environment, produce soft error and hard fault, cause navigating (navigation neceiver), radar-probing system (Active Phased Array Radar), data network (the AFDX network switch), communication (optical fiber/bus), high speed computer system, avionic device, engine (FADEC), telex system, automatic Pilot technology, flight alarm, display screen, there is blank screen in other flight system etc. containing electron device, crash, reset, restart, loss of data, the safety hazards such as order loss.In order to set up protection and appraisement system, quantitatively characterizing must be carried out to the harm of atmospheric neutron single particle effect.And atmospheric neutron fluence rate is the important parameter of harm quantitatively characterizing.To this, 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 flow diagram that the embodiment of the present invention one proposes, and as shown in Figure 1, comprises the following steps:

S101 obtains mission profile and the energy threshold of current electronic device, and described mission profile comprises flying height, latitude, the longitude of current electronic device at each task point.

In the present embodiment, obtain the energy threshold of current electronic device, comprising: the energy threshold obtaining crucial Sensitive Apparatus in current electronic device, as the energy threshold of current electronic device.

The mission profile proposed in the embodiment of the present invention is operation flight height, latitude, the parameter such as longitude and flight beginning and ending time of current electronic device.Because the factor affecting atmospheric neutron fluence rate comprises cut-off rigidity, and cut-off rigidity is determined jointly by longitude and latitude, therefore, the present invention is longitude, latitude two parameters by the latitude parameter optimized in Boeing model, can improve the computational accuracy of atmospheric neutron fluence rate.

S102 calculates the first air neutron fluence rate of each task point according to the mission profile of current electronic device and energy threshold;

S103, according to the energy threshold of described current electronic device, adopts Boeing model to calculate the second air neutron fluence rate of each task point, adopts NASA model to calculate the third-largest gas neutron fluence rate of each task point;

S104 calculates the average of the first air neutron fluence rate of each task point, the second air neutron fluence rate and the third-largest gas neutron fluence rate respectively, as the atmospheric neutron fluence rate of this task point.

In embodiments of the present invention, the first air neutron fluence rate of each task point is calculated according to the mission profile of current electronic device and energy threshold, specifically comprise: show according to global longitude and latitude is corresponding with cut-off rigidity, as shown in table 1, search the cut-off rigidity that the longitude and latitude at each task point place in described mission profile is corresponding under 12.2km altitudes; According to each task point current by rigidity, utilize the fluence rate multiple table of cut-off rigidity and 12.2km height, as shown in table 2, search the first fluence rate multiple of each task point ground, fluence rate relative New York fluence rate when 12.2km height; According to described first fluence rate multiple, utilize height latitude multiple table, as shown in table 3, search the second fluence rate multiple of ground, the fluence rate of each task point when flying height relative New York fluence rate; The energy calculating each task point according to the second fluence rate multiple is greater than fluence rate during 10MeV; According to the energy threshold of described current electronic device, utilize energy multiple table, as shown in table 4, carry out energy threshold correction, obtain the first air neutron fluence rate of each task point, comprise the most bad value, representative value (height 12200,45 degree) and mean value.

The global longitude and latitude of table 1 is corresponding with cut-off rigidity to be shown (1)

The global longitude and latitude of table 1 is corresponding with cut-off rigidity to be shown (2)

Show global longitude and latitude is corresponding with cut-off rigidity, north latitude is just, south latitude is negative; Longitude is expressed as east longitude 0-360 °, is 360-x ° of W.As 30S, 120W, be-30N, 240E.

The fluence rate multiple table of table 2 cut-off rigidity 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

Table 3 height latitude multiple table (1)

Ground fluence rate multiple table, each high latitude of table 3 relative New York (2)

Table 4 energy threshold multiple table

Sequence number	Energy (MeV)	Neutron fluence rate (#/cm2h)	With 1-10MeV multiple proportion
				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 calculating each task point according to the second fluence rate multiple proposed in embodiments of the present invention is greater than fluence rate during 10MeV, specifically comprises:

Acquisition ground, New York energy is greater than standard fluence rate during 10MeV; The fluence rate when energy of described each task point is greater than 10MeV is calculated according to described standard fluence rate and described second fluence rate multiple.

In embodiments of the present invention, step S103 specifically comprises: adopt Boeing model to calculate the above atmospheric neutron fluence rate of energy threshold of each task point, draws the most bad value, representative value (height 12200,45 degree) and mean value.And adopt NASA model to calculate the above atmospheric neutron fluence rate of energy threshold of each task point.Find out the most bad value, representative value (height 12200,45 degree) and mean value.

In the embodiment of the present invention, step S104, specifically comprises:

Extract the most bad value of the first air neutron fluence rate of each task point, representative value and mean value;

Extract the most bad value of the second air neutron fluence rate of each task point, representative value and mean value;

Extract the most bad value of the third-largest gas neutron fluence rate of each task point, representative value and mean value;

Calculate the average of the average of the most bad value of each task point, the average of representative value and mean value respectively, as the atmospheric neutron fluence rate of this task point.

Flying area provided by the invention atmospheric neutron fluence rate multi-model composite computing method, by longitude, the latitude determination cut-off rigidity of task point, latitude impact in Boeing model is modified to the cut-off rigidity determined by longitude, latitude, and be applicable to any height, solve NASA model in the inapplicable problem of floor level, each different-energy threshold value is on the impact of atmospheric neutron fluence rate by quantifying, more accurately calculates the atmospheric neutron fluence rate of electron device at task point place.

Fig. 2 is a kind of flying area atmospheric neutron fluence rate multi-model composite computing apparatus module figure that the embodiment of the present invention two proposes, and as shown in Figure 2, described device comprises:

Acquisition module 201, for obtaining mission profile and the energy threshold of current electronic device, described mission profile comprises flying height, latitude, the longitude of current electronic device at each task point.In the present embodiment, the energy threshold of described current electronic device is the energy threshold of crucial Sensitive Apparatus in current electronic device.

First computing module 202, for calculating the first air neutron fluence rate of each task point according to the mission profile of current electronic device and energy threshold;

Second computing module 203, for the energy threshold according to described current electronic device, adopts Boeing model to calculate the second air neutron fluence rate of each task point, adopts NASA model to calculate the third-largest gas neutron fluence rate of each task point;

Mean value computation module 204, for calculating the average of the first air neutron fluence rate of each task point, the second air neutron fluence rate and the third-largest gas neutron fluence rate respectively, as the atmospheric neutron fluence rate of this task point.

The first computing module 202 in the present embodiment specifically comprises:

Cut-off rigidity searches unit, for showing according to global longitude and latitude is corresponding with cut-off rigidity, searches the cut-off rigidity that the longitude and latitude at each task point place in described mission profile is corresponding under 12.2km altitudes;

First fluence rate multiple searches unit, for according to each task point current by rigidity, utilize the fluence rate multiple table of cut-off rigidity and 12.2km height, search the first fluence rate multiple of each task point ground, fluence rate relative New York fluence rate when 12.2km height;

Second fluence rate multiple searches unit, for according to described first fluence rate multiple, utilizes height latitude multiple table, searches the second fluence rate multiple of ground, the fluence rate of each task point when flying height relative New York fluence rate;

Condition computing unit, the energy for calculating each task point according to the second fluence rate multiple is greater than fluence rate during 10MeV;

Energy threshold amending unit, for the energy threshold according to described current electronic device, utilizes energy multiple table to carry out energy threshold correction, obtains the first air neutron fluence rate of each task point.

Condition computing unit in the present embodiment specifically comprises:

Obtain subelement, for obtaining standard fluence rate when ground, New York energy is greater than 10MeV;

Computation subunit, the energy for calculating described each task point according to described standard fluence rate and described second fluence rate multiple is greater than fluence rate during 10MeV.

Mean value computation module 204 in the present embodiment specifically comprises:

First extraction unit, for extracting the most bad value of the first air neutron fluence rate of each task point, representative value and mean value;

Second extraction unit, for extracting the most bad value of the second air neutron fluence rate of each task point, representative value and mean value;

3rd extraction unit, for extracting the most bad value of the third-largest gas neutron fluence rate of each task point, representative value and mean value;

Average calculation unit, for calculating the average of the average of the most bad value of each task point, the average of representative value and mean value respectively, as the atmospheric neutron fluence rate of this task point.

The flying area atmospheric neutron fluence rate multi-model composite computing method that the present invention proposes and device, specify that cut-off rigidity is the key factor affecting atmospheric neutron fluence rate, instead of single latitude, the latitude parameter optimized in Boeing model is longitude, latitude two parameters, and the method can calculate atmospheric neutron fluence rate more accurately.

The invention solves NASA model in the inapplicable problem of floor level, be suitable for and various height.

The present invention adds the impact of energy range.Boeing model and NASA mainly calculate the atmospheric neutron fluence rate that energy is 1-10MeV, and consider that the energy threshold of different components is different, the threshold value of some devices may be 4MeV, 5MeV, and the method is by adding the impact of energy range, computing power can be greater than 1MeV, be greater than 2MeV to the neutron fluence rate being greater than the different-energy scopes such as 10MeV, make result of calculation more accurate.

The present invention is by computation of mean values, and improve computational accuracy, user can choose the atmospheric neutron fluence rate of the most applicable data as task point by rule of thumb.

Through the above description of the embodiments, those skilled in the art can be well understood to the present invention can by hardware implementing, and the mode that also can add necessary general hardware platform by software realizes.Based on such understanding, technical scheme of the present invention can embody with the form of software product, it (can be CD-ROM that this software product can be stored in a non-volatile memory medium, USB flash disk, portable hard drive etc.) in, comprise some instructions and perform method described in each embodiment of the present invention in order to make a computer equipment (can be personal computer, server, or the network equipment etc.).

It will be appreciated by those skilled in the art that accompanying drawing is the schematic diagram of a preferred embodiment, the module in accompanying drawing or flow process might not be that enforcement the present invention is necessary.

It will be appreciated by those skilled in the art that the module in the device in embodiment can carry out being distributed in the device of embodiment according to embodiment description, also can carry out respective change and be arranged in the one or more devices being different from the present embodiment.The module of above-described embodiment can merge into a module, also can split into multiple submodule further.

Be only several specific embodiment of the present invention above, but the present invention is not limited thereto, the changes that any person skilled in the art can think of all should fall into protection scope of the present invention.

Claims (10)

1. a flying area atmospheric neutron fluence rate multi-model composite computing method, is characterized in that, comprising:

Obtain mission profile and the energy threshold of current electronic device, described mission profile comprises flying height, latitude, the longitude of current electronic device at each task point;

The first air 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 described current electronic device, adopt Boeing model to calculate the second air neutron fluence rate of each task point, adopt NASA model to calculate the third-largest gas neutron fluence rate of each task point;

Calculate the average of the first air neutron fluence rate of each task point, the second air neutron fluence rate and the third-largest gas neutron fluence rate respectively, as the atmospheric neutron fluence rate of this task point.

2. method according to claim 1, is characterized in that, the energy threshold of described acquisition current electronic device, comprising:

Obtain the energy threshold of crucial Sensitive Apparatus in current electronic device, as the energy threshold of current electronic device.

3. method according to claim 1, is characterized in that, the described mission profile according to current electronic device and energy threshold calculate the first air neutron fluence rate of each task point, specifically comprise:

Show according to global longitude and latitude is corresponding with cut-off rigidity, search the cut-off rigidity that the longitude and latitude at each task point place in described mission profile is corresponding under 12.2km altitudes;

The cut-off rigidity current according to each task point, utilizes the fluence rate multiple table of cut-off rigidity and 12.2km height, searches the first fluence rate multiple of each task point ground, fluence rate relative New York fluence rate when 12.2km height;

According to described first fluence rate multiple, utilize height latitude multiple table, search the second fluence rate multiple of ground, the fluence rate of each task point when flying height relative New York fluence rate;

The energy calculating each task point according to the second fluence rate multiple is greater than fluence rate during 10MeV;

According to the energy threshold of described current electronic device, utilize energy multiple table to carry out energy threshold correction, obtain the first air neutron fluence rate of each task point.

4. according to the method that claim 3 is stated, it is characterized in that, the described energy calculating each task point according to the second fluence rate multiple is greater than fluence rate during 10MeV, specifically comprises:

Acquisition ground, New York energy is greater than standard fluence rate during 10MeV;

The fluence rate when energy of described each task point is greater than 10MeV is calculated according to described standard fluence rate and described second fluence rate multiple.

5. method according to claim 1, is characterized in that, the average of described the first air neutron fluence rate, the second air neutron fluence rate and the third-largest gas neutron fluence rate that calculate each task point respectively, specifically comprises:

Extract the most bad value of the first air neutron fluence rate of each task point, representative value and mean value;

Extract the most bad value of the second air neutron fluence rate of each task point, representative value and mean value;

Extract the most bad value of the third-largest gas neutron fluence rate of each task point, representative value and mean value;

Calculate the average of the average of the most bad value of each task point, the average of representative value and mean value respectively, as the atmospheric neutron fluence rate of this task point.

6. a flying area atmospheric neutron fluence rate multi-model composite computing device, it is characterized in that, described device comprises:

Acquisition module, for obtaining mission profile and the energy threshold of current electronic device, described mission profile comprises flying height, latitude, the longitude of current electronic device at each task point;

First computing module, for calculating the first air neutron fluence rate of each task point according to the mission profile of current electronic device and energy threshold;

Second computing module, for the energy threshold according to described current electronic device, adopts Boeing model to calculate the second air neutron fluence rate of each task point, adopts NASA model to calculate the third-largest gas neutron fluence rate of each task point;

Mean value computation module, for calculating the average of the first air neutron fluence rate of each task point, the second air neutron fluence rate and the third-largest gas neutron fluence rate respectively, as the atmospheric neutron fluence rate of this task point.

7. device according to claim 6, is characterized in that, the energy threshold of described current electronic device is the energy threshold of crucial Sensitive Apparatus in current electronic device.

8. device according to claim 6, is characterized in that, described first computing module specifically comprises:

Cut-off rigidity searches unit, for showing according to global longitude and latitude is corresponding with cut-off rigidity, searches the cut-off rigidity that the longitude and latitude at each task point place in described mission profile is corresponding under 12.2km altitudes;

First fluence rate multiple searches unit, for according to each task point current by rigidity, utilize the fluence rate multiple table of cut-off rigidity and 12.2km height, search the first fluence rate multiple of each task point ground, fluence rate relative New York fluence rate when 12.2km height;

Second fluence rate multiple searches unit, for according to described first fluence rate multiple, utilizes height latitude multiple table, searches the second fluence rate multiple of ground, the fluence rate of each task point when flying height relative New York fluence rate;

Condition computing unit, the energy for calculating each task point according to the second fluence rate multiple is greater than fluence rate during 10MeV;

Energy threshold amending unit, for the energy threshold according to described current electronic device, utilizes energy multiple table to carry out energy threshold correction, obtains the first air neutron fluence rate of each task point.

9. device according to claim 8, is characterized in that, described condition computing unit specifically comprises:

Obtain subelement, for obtaining standard fluence rate when ground, New York energy is greater than 10MeV;

Computation subunit, the energy for calculating described each task point according to described standard fluence rate and described second fluence rate multiple is greater than fluence rate during 10MeV.

10. device according to claim 6, is characterized in that, described mean value computation module specifically comprises:

First extraction unit, for extracting the most bad value of the first air neutron fluence rate of each task point, representative value and mean value;

Second extraction unit, for extracting the most bad value of the second air neutron fluence rate of each task point, representative value and mean value;

3rd extraction unit, for extracting the most bad value of the third-largest gas neutron fluence rate of each task point, representative value and mean value;

Average calculation unit, for calculating the average of the average of the most bad value of each task point, the average of representative value and mean value respectively, as the atmospheric neutron fluence rate of this task point.