CN108345026A - A method of calculating power spectrum after charged particle protective layer - Google Patents

A method of calculating power spectrum after charged particle protective layer Download PDF

Info

Publication number
CN108345026A
CN108345026A CN201810136597.0A CN201810136597A CN108345026A CN 108345026 A CN108345026 A CN 108345026A CN 201810136597 A CN201810136597 A CN 201810136597A CN 108345026 A CN108345026 A CN 108345026A
Authority
CN
China
Prior art keywords
protective materials
energy
proton
spectrum
ion
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
CN201810136597.0A
Other languages
Chinese (zh)
Other versions
CN108345026B (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201810136597.0A priority Critical patent/CN108345026B/en
Publication of CN108345026A publication Critical patent/CN108345026A/en
Application granted granted Critical
Publication of CN108345026B publication Critical patent/CN108345026B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/36Measuring spectral distribution of X-rays or of nuclear radiation spectrometry

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Measurement Of Radiation (AREA)

Abstract

A method of power spectrum after charged particle protective layer being calculated, the present invention relates to the methods of power spectrum after calculating charged particle protective layer.The purpose of the present invention is to solve the prior arts there are the calculating speed of power spectrum after orbit of charged particle protective layer is slow, calculates that the time is long and space particle radiation can make the short service life of electronic device, the problem of bringing about great losses.Process is:Determine the differential energy spectrum of track incoming particle;Determine protective materials and protective materials thickness;Determine the incident angle of incoming particle;If incoming particle is electronics, E ' is calculated1With the relationship of range and protective materials thickness;If incoming particle is proton or ion, E ' is calculated2、E′3、E′4With E '5Respectively with the relationship of range and protective materials thickness;Calculate the inelastic cross section of incoming particle and protective materials;Calculate separately the differential energy spectrum after electronics, proton and ion protective layer.The present invention is used for space environment effect, nuclear science and applied technical field.

Description

A method of calculating power spectrum after charged particle protective layer
Technical field
The present invention relates to the methods of power spectrum after calculating charged particle protective layer, belong to space environment effect, nuclear science and answer Use technical field.
Background technology
Spacecraft is run in space orbit, it will is influenced by a variety of environmental factors:As particle radiation, microgravity and Elemental oxygen etc..With the rapid development of electronic technology, the application of electronic component is more and more extensive.Satellite, spaceship and On space shuttle, extensively using electronic device, to realize various functions.For electronic device, particle radiation is most fatal ring Border factor.The early stage spacecraft of transmitting once repeatedly failed because of damage.With the development of science and technology, although the spacecrafts such as satellite The example for generating critical failure is fewer and fewer, but still happens occasionally.Space particle radiation can make the short service life of satellite, from And it brings about great losses.Even if using the shadow for the methods of encapsulating and adding protection and cannot guarantee that electronic device avoids radiation completely It rings.
The content being concerned about the most in aerospace engineering relates generally to two aspects:First, quick and precisely assessing sensing unit (electronics Component, thermal control coating, conducting wire, optical component etc.) radiation absorbed dose that is subject to of place;Second is that for the protection of Sensitive Apparatus Design finds a kind of fast appraisement method.
But the prior art is there are the calculating speed of power spectrum after orbit of charged particle protective layer is slow, it is long to calculate the time, and Space particle radiation can make the short service life of electronic device, the problem of bringing about great losses.
Invention content
The purpose of the present invention is to solve the prior art, there are the calculating speeds of power spectrum after orbit of charged particle protective layer Slowly, it calculates that the time is long and space particle radiation can make the short service life of electronic device, the problem of bringing about great losses, And propose it is a kind of calculating charged particle protective layer after power spectrum method.
It is a kind of calculating charged particle protective layer after the method detailed process of power spectrum be:
Step 1: determining the differential energy spectrum f (E ' of track incoming particle1)、f(E′2)、f(E′3)、f(E′4) and f (E '5);
Step 2: determining protective materials and protective materials thickness t;
Step 3: determining the incident angle α of incoming particle;
Step 4: if incoming particle is electronics, according to the protective materials thickness t that step 2 determines, rail in step 1 is calculated Differential energy spectrum f (the E ' of road incident electron1) in ENERGY E '1With the relationship of range and protective materials thickness;
Step 5:
If 1) incoming particle is proton, what the protective materials thickness t determined according to step 2 and step 3 determined enters radion The incident angle α of son, calculates separately the differential energy spectrum f (E ' of step 1 middle orbit incident proton2) in ENERGY E '2With f (E '3) in ENERGY E '3With the relationship of range and protective materials thickness;
If 2) incoming particle is ion, what the protective materials thickness t determined according to step 2 and step 3 determined enters radion The incident angle α of son, calculates separately the differential energy spectrum f (E ' of step 1 middle orbit incident ion4) in ENERGY E '4With f (E '5) in ENERGY E '5With the relationship of range and protective materials thickness;
Step 6: calculating the inelastic cross section of incoming particle and protective materials;
Step 7: if incoming particle is electronics, according to the differential energy spectrum f (E ' of step 11), the protective materials of step 2 it is thick The ENERGY E that degree t, step 4 obtain '1With the inelastic cross section σ of step 6, the differential after radiation belt of the earth electronic protective layers is calculated Power spectrum;
Step 8:
1) if incoming particle is proton, according to the differential energy spectrum f (E ' of step 12), the protective materials thickness t of step 2, The ENERGY E that the incident angle α of step 3, step 5 obtain '2With the inelastic cross section σ of step 6, radiation belt of the earth proton is calculated Differential energy spectrum after protective layer;
According to the differential energy spectrum f (E ' of step 13), the protective materials thickness t of step 2, step 3 incident angle α, step Rapid five obtain ENERGY E '3With the inelastic cross section σ of step 6, the differential energy spectrum after solar cosmic ray proton protective layer is calculated;
2) if incoming particle is ion, according to the differential energy spectrum f (E ' of step 14), the protective materials thickness t of step 2, The ENERGY E that the incident angle α of step 3, step 5 obtain '4With the inelastic cross section σ of step 6, solar cosmic ray ion is calculated Differential energy spectrum after protective layer;
According to the differential energy spectrum f (E ' of step 15), the protective materials thickness t of step 2, step 3 incident angle α, step Rapid five obtain ENERGY E '5With the inelastic cross section σ of step 6, the differential energy spectrum after galactic comic ray ion protective layer is calculated.
Beneficial effects of the present invention are:
Spacecraft of the present invention is mainly comprehensive by spaces such as the proton of different-energy, electronics and heavy ions with radiation-sensitive member Close the influence of radiation environment factor;The present invention is based on the original power spectrums of the original charged particle of common track in engineering, are distinguishing just On the basis of ion and anion, for different overcoat thickness, achieve the purpose that quickly to calculate power spectrum after protective layer, step letter Single, easily operated, calculating speed is fast, is greatly lowered and calculates the time, extends the working life contracting of electronic device, reduces damage It loses, is of great significance to the in-orbit reliability of spacecraft, there is apparent advantage and be widely applied foreground.
Fig. 4 gives after 1mm, 2mm, 3mm Al protective materials, power spectrum after the protective layer of terrestrial radiation having electronic.By scheming As it can be seen that with the increase of overcoat thickness, power spectrum gradually decreases after protective layer.It is 10 in terrestrial radiation having electronic-1When MeV, 3mmAl Afterwards after power spectrum, 2mmAl after power spectrum and 1mmAl power spectrum 102-104MeV-1.cm-2.s-1Between, original power spectrum is 106-107MeV-1.cm-2.s-1Between;Fig. 5 gives after 1mm, 2mm, 3mm Al protective materials, energy after the protective layer of radiation belt of the earth proton Spectrum.As seen from the figure, with the increase of overcoat thickness, power spectrum gradually decreases after low energy proton protective layer, and energetic portions variation is not Greatly.
Description of the drawings
Fig. 1 is that incident electron vertical incidence of the present invention calculates path schematic diagram, L1For Actual path, L2To calculate path;
Fig. 2 is that incident ion oblique incidence of the present invention calculates path schematic diagram, L3For path, t is thickness, and α is incidence angle Degree;
Fig. 3 is the original power spectrum schematic diagram of track;
Fig. 4 is power spectrum schematic diagram of calculation result of the proton after protective layer;
Fig. 5 is power spectrum schematic diagram of calculation result of the electronics after protective layer.
Specific implementation mode
Specific implementation mode one:The method detailed process of power spectrum after a kind of calculating charged particle protective layer of present embodiment For:
Spacecraft of the present invention is mainly comprehensive by spaces such as the proton of different-energy, electronics and heavy ions with radiation-sensitive member Close the influence of radiation environment factor.The present invention is based on the original power spectrums of the original charged particle of common track in engineering, are distinguishing just On the basis of ion and anion, for different overcoat thickness, achieve the purpose that quickly to calculate power spectrum after protective layer.
It is according to the present invention a kind of based on the original power spectrum of the common original charged particle of track, application in engineering Including material, device, electronic system and structure.The technology is characterized in that, based on common track grandfather tape electrochondria in engineering The original power spectrum of son, calculates the ions energy spectrum after different type and different-thickness protective layer.
Different types of incoming particle, especially positively charged charged particle (ion) and electronegative charged particle (electricity Son), in the transport process of material and device, generate degree of injury and path difference.The study found that electronics is based primarily upon vertically Incidence, as shown in Figure 1;It is directly incident that ion is based primarily upon practical incident path mode, as shown in Figure 2.
Step 1: determining the differential energy spectrum f (E ' of track incoming particle1)、f(E′2)、f(E′3)、f(E′4) and f (E '5);
Step 2: determining that protective materials and protective materials thickness t (artificially select protective materials, each protective materials has it The thickness of itself);
Step 3: determining the incident angle α (artificial settings) of incoming particle;
Step 4: if incoming particle is electronics, range is unrelated with the thickness of protective materials, is determined according to step 2 anti- Protective material thickness t calculates the differential energy spectrum f (E ' of step 1 middle orbit incident electron1) in ENERGY E '1With range and protective materials The relationship of thickness;
Step 5:
1) if incoming particle is proton, range is related with the thickness of protective materials, the protection material determined according to step 2 The incident angle α for expecting the incoming particle that thickness t and step 3 are determined, calculates separately the differential energy of step 1 middle orbit incident proton Compose f (E '2) in ENERGY E '2With f (E '3) in ENERGY E '3With the relationship of range and protective materials thickness;
2) if incoming particle is ion, range is related with the thickness of protective materials, the protection material determined according to step 2 The incident angle α for expecting the incoming particle that thickness t and step 3 are determined, calculates separately the differential energy of step 1 middle orbit incident ion Compose f (E '4) in ENERGY E '4With f (E '5) in ENERGY E '5With the relationship of range and protective materials thickness;
Step 6: calculating the inelastic cross section of incoming particle and protective materials;
Step 7: if incoming particle is electronics, range is unrelated with the thickness of protective materials.According to the differential energy of step 1 Compose f (E '1), the protective materials thickness t of step 2, step 4 obtain ENERGY E '1With the inelastic cross section σ of step 6, calculate ground Differential energy spectrum after spherical radiation having electronic protective layer;
Step 8:
1) if incoming particle is proton, range is related with the thickness of protective materials, according to the differential energy spectrum f of step 1 (E′2), the protective materials thickness t of step 2, the incident angle α of step 3, step 5 obtain ENERGY E '2It is non-with step 6 Elastic cross-section σ calculates the differential energy spectrum after radiation belt of the earth proton protective layer;
According to the differential energy spectrum f (E ' of step 13), the protective materials thickness t of step 2, step 3 incident angle α, step Rapid five obtain ENERGY E '3With the inelastic cross section σ of step 6, the differential energy spectrum after solar cosmic ray proton protective layer is calculated;
2) if incoming particle is ion, range is related with the thickness of protective materials, according to the differential energy spectrum f of step 1 (E′4), the protective materials thickness t of step 2, the incident angle α of step 3, step 5 obtain ENERGY E '4It is non-with step 6 Elastic cross-section σ calculates the differential energy spectrum after solar cosmic ray ion protective layer;
According to the differential energy spectrum f (E ' of step 15), the protective materials thickness t of step 2, step 3 incident angle α, step Rapid five obtain ENERGY E '5With the inelastic cross section σ of step 6, the differential energy spectrum after galactic comic ray ion protective layer is calculated.
Differential energy spectrum f (E ' are found out in step 11)、f(E′2)、f(E′3)、f(E′4) and f (E '5), according to f (E '1)、f (E′2)、f(E′3)、f(E′4) and f (E '5) energy and quantity of errant incoming particle can be obtained, what step 4 and step 5 were found out E′1、E′2、E′3、E′4、E′5It is energy, the f (E ' of step 1 is corresponded to the energy1)、f(E′2)、f(E′3)、f(E′4) and f (E′5), so that it may the quantity of step 4 and step 5 to energy is obtained respectively, Step 7: the f (E ' in step 81)、f(E′2)、f (E′3)、f(E′4) and f (E '5) it is that step 1 is found out, Step 7: the S (E ' in step 81)、S(E′2)、S(E′3)、S(E′4) With S (E '5) in E '1、E′2、E′3、E′4With E '5It is that step 4 and step 5 are found out.
The present disclosure applies equally to space other radiation conditions, such as Jupiter, Saturn, Mars, the moon etc., Jupiter, Saturn Equally it is terrestrial radiation having electronic, radiation belt of the earth proton, solar cosmic ray proton, solar cosmic ray ion and galactic comic ray Ion;The moon only has solar cosmic ray ion and galactic comic ray ion;The computational methods of power spectrum are the same as phase of the present invention after protective layer Together.
Specific implementation mode two:The present embodiment is different from the first embodiment in that:Rail is determined in the step 1 The differential energy spectrum of road incoming particle;Detailed process is:
Satellite orbital altitude and inclination of satellite orbit are determined according to satellite orbit;
Determined respectively according to satellite orbital altitude and inclination of satellite orbit terrestrial radiation having electronic, radiation belt of the earth proton, Solar cosmic ray proton, solar cosmic ray ion and galactic comic ray ion;
AE8 or AE9 models are selected according to terrestrial radiation having electronic, obtain the differential energy spectrum f (E ' of track incident electron1);
AP8 or AP9 models are selected according to radiation belt of the earth proton, obtain the differential energy spectrum f (E ' of track incident proton2);
ESP models are selected according to solar cosmic ray proton, obtain the differential energy spectrum f (E ' of track incident proton3);
CREME96 models are selected according to solar cosmic ray ion, obtain the differential energy spectrum f (E ' of track incident ion4);
CREME96 models are selected according to galactic comic ray ion, obtain the differential energy spectrum f (E ' of track incident ion5);
Particle includes ion, electronics and proton.
Other steps and parameter are same as the specific embodiment one.
Specific implementation mode three:The present embodiment is different from the first and the second embodiment in that:If in the step 4 Incoming particle is electronics, then range is unrelated with the thickness of protective materials, according to the protective materials thickness t that step 2 determines, is calculated Differential energy spectrum f (the E of step 1 middle orbit incident electron1') in ENERGY E1' the relationship with range and protective materials thickness, it is specific public Formula is:
Wherein, Ri(E1) it is residual range of the radiation belt of the earth incident electron in protective materials, pass through EGS, GEANT etc. Software is calculated;RiFor residual range;E1For energy (radiation belt of the earth electron energy).
Other steps and parameter are the same as one or two specific embodiments.
Specific implementation mode four:Unlike one of present embodiment and specific implementation mode one to three:The step 5 1) if incoming particle is proton in, and range is related with the thickness of protective materials, the protective materials thickness t determined according to step 2 The incident angle α of the incoming particle determined with step 3, calculates separately the differential energy spectrum f of step 1 middle orbit incident proton (E′2) in ENERGY E '2With f (E '3) in ENERGY E '3With the relationship of range and protective materials thickness, specific formula is:
Wherein, Ri(E2) it is residual range of the radiation belt of the earth incident proton in protective materials, pass through SRIM, GEANT etc. Software is calculated;RiFor residual range;E2For energy (radiation belt of the earth proton energy);
Ri(E3) it is residual range of the solar cosmic ray incident proton in protective materials, pass through the softwares such as SRIM, GEANT It is calculated;RiFor residual range;E3For energy (solar cosmic ray proton energy);
2) if incoming particle is ion, range is related with the thickness of protective materials, the protection material determined according to step 2 The incident angle α for expecting the incoming particle that thickness t and step 3 are determined, calculates separately the differential energy of step 1 middle orbit incident ion Compose f (E '4) in ENERGY E '4With f (E '5) in ENERGY E '5With the relationship of range and protective materials thickness, specific formula is:
Wherein, Ri(E4) it is residual range of the solar cosmic ray incident ion in protective materials, pass through SRIM, GEANT etc. Software is calculated;RiFor residual range;E4For energy (solar cosmic ray ion energy);
Ri(E5) it is residual range of the galactic comic ray incident ion in protective materials, pass through the softwares such as SRIM, GEANT It is calculated;RiFor residual range;E5For energy (galactic comic ray ion energy).
Other steps and parameter are identical as one of specific implementation mode one to three.
Specific implementation mode five:Unlike one of present embodiment and specific implementation mode one to four:The step 6 The middle inelastic cross section for calculating incoming particle and protective materials, specific formula are:
Wherein, N is avogadros constant, and A is protective materials atomic number, and M is protective materials atomic weight.
Other steps and parameter are identical as one of specific implementation mode one to four.
Specific implementation mode six:Unlike one of present embodiment and specific implementation mode one to five:The step 7 If incoming particle is electronics in, and range is unrelated with the thickness of protective materials.According to the differential energy spectrum f (E ' of step 11), step The ENERGY E that two protective materials thickness t, step 4 obtain '1With the inelastic cross section σ of step 6, terrestrial radiation having electronic is calculated Differential energy spectrum after protective layer, specific formula are:
Wherein, S (E '1) and S (E1) be energy be E '1And E1Prevention sheet of the terrestrial radiation having electronic in protective materials Neck, is calculated by softwares such as EGS, GEANT;S(E′1) obtained according to step 4 ENERGY E '1It obtains.
Other steps and parameter are identical as one of specific implementation mode one to five.
Specific implementation mode seven:Unlike one of present embodiment and specific implementation mode one to six:The step 8 1) if incoming particle is proton in, and range is related with the thickness of protective materials, according to the differential energy spectrum f (E ' of step 12), step The ENERGY E that rapid two protective materials thickness t, the incident angle α of step 3, step 5 obtain2' the inelastic cross section with step 6 σ, calculates the differential energy spectrum after radiation belt of the earth proton protective layer, and specific formula is:
Wherein, S (E '2) and S (E2) be energy be E '2And E2Prevention sheet of the radiation belt of the earth proton in protective materials Neck, is calculated by softwares such as SRIM, GEANT;S(E′2) obtained according to step 5 ENERGY E '2It obtains;
According to the differential energy spectrum f (E ' of step 13), the protective materials thickness t of step 2, step 3 incident angle α, step Rapid five obtain ENERGY E '3With the inelastic cross section σ of step 6, the differential energy spectrum after solar cosmic ray proton protective layer is calculated, Specifically formula is:
Wherein, S (E '3) and S (E3) be energy be E '3And E3Prevention sheet of the solar cosmic ray proton in protective materials Neck, is calculated by softwares such as SRIM, GEANT;S(E′3) obtained according to step 5 ENERGY E '3It obtains;
2) if incoming particle is ion, range is related with the thickness of protective materials, according to the differential energy spectrum f of step 1 (E′4), the protective materials thickness t of step 2, the incident angle α of step 3, step 5 obtain ENERGY E '4It is non-with step 6 Elastic cross-section σ, calculates the differential energy spectrum after solar cosmic ray ion protective layer, and specific formula is:
Wherein, S (E '4) and S (E4) be energy be E '4And E4Prevention sheet of the solar cosmic ray ion in protective materials Neck, is calculated by softwares such as SRIM, GEANT;S(E′4) obtained according to step 5 ENERGY E '4It obtains;
According to the differential energy spectrum f (E ' of step 15), the protective materials thickness t of step 2, step 3 incident angle α, step Rapid five obtain ENERGY E '5With the inelastic cross section σ of step 6, the differential energy spectrum after galactic comic ray ion protective layer is calculated, Specifically formula is:
Wherein, S (E '5) and S (E5) be energy be E '5And E5Prevention sheet of the galactic comic ray ion in protective materials Neck, is calculated by softwares such as SRIM, GEANT;S(E′5) obtained according to step 5 ENERGY E '5It obtains.
Other steps and parameter are identical as one of specific implementation mode one to six.
Beneficial effects of the present invention are verified with following embodiment:
Embodiment one:
The quick calculation method process of power spectrum is after a kind of charged particle protective layer of the present embodiment:
For LEO tracks (800km, 98 °), Al protective materials calculates separately the electronics and proton spectrum after protective layer.LEO The original differential energy spectrum of track terrestrial radiation having electronic and proton is as shown in Figure 3.Fig. 4 gives through 1mm, 2mm, 3mm Al protection After material, power spectrum after the protective layer of terrestrial radiation having electronic.As seen from the figure, with the increase of overcoat thickness, power spectrum after protective layer It gradually decreases.Fig. 5 gives after 1mm, 2mm, 3mm Al protective materials, power spectrum after the protective layer of radiation belt of the earth proton.By Figure as it can be seen that with overcoat thickness increase, power spectrum gradually decreases after low energy proton protective layer, and energetic portions variation is little.
The present invention can also have other various embodiments, without deviating from the spirit and substance of the present invention, this field Technical staff makes various corresponding change and deformations in accordance with the present invention, but these corresponding change and deformations should all belong to The protection domain of appended claims of the invention.

Claims (7)

1. a kind of method of power spectrum after calculating charged particle protective layer, it is characterised in that:The method detailed process is:
Step 1: determining the differential energy spectrum f (E of track incoming particle1′)、f(E2′)、f(E3′)、f(E4') and f (E5′);
Step 2: determining protective materials and protective materials thickness t;
Step 3: determining the incident angle α of incoming particle;
Step 4: if incoming particle is electronics, according to the protective materials thickness t that step 2 determines, calculates step 1 middle orbit and enter Differential energy spectrum f (the E of radio1') in ENERGY E1' the relationship with range and protective materials thickness;
Step 5:
1) if incoming particle is proton, the incoming particle that the protective materials thickness t determined according to step 2 and step 3 determine Incident angle α calculates separately the differential energy spectrum f (E of step 1 middle orbit incident proton2') in ENERGY E2' and f (E3') in energy E3' the relationship with range and protective materials thickness;
2) if incoming particle is ion, the incoming particle that the protective materials thickness t determined according to step 2 and step 3 determine Incident angle α calculates separately the differential energy spectrum f (E of step 1 middle orbit incident ion4') in ENERGY E4' and f (E5') in energy E5' the relationship with range and protective materials thickness;
Step 6: calculating the inelastic cross section of incoming particle and protective materials;
Step 7: if incoming particle is electronics, according to the differential energy spectrum f (E of step 11'), the protective materials thickness t of step 2, The ENERGY E that step 4 obtains1' inelastic cross section the σ with step 6, calculates the differential energy after radiation belt of the earth electronic protective layers Spectrum;
Step 8:
1) if incoming particle is proton, according to the differential energy spectrum f (E of step 12'), the protective materials thickness t of step 2, step 3 Incident angle α, the obtained ENERGY E of step 52' inelastic cross section the σ with step 6, calculates radiation belt of the earth proton protective layer Differential energy spectrum afterwards;
According to the differential energy spectrum f (E of step 13'), the protective materials thickness t of step 2, the incident angle α of step 3, step 5 Obtained ENERGY E3' inelastic cross section the σ with step 6, calculates the differential energy spectrum after solar cosmic ray proton protective layer;
2) if incoming particle is ion, according to the differential energy spectrum f (E of step 14'), the protective materials thickness t of step 2, step 3 Incident angle α, the obtained ENERGY E of step 54' inelastic cross section the σ with step 6, calculates solar cosmic ray ion protective layer Differential energy spectrum afterwards;
According to the differential energy spectrum f (E of step 15'), the protective materials thickness t of step 2, the incident angle α of step 3, step 5 Obtained ENERGY E5' inelastic cross section the σ with step 6, calculates the differential energy spectrum after galactic comic ray ion protective layer.
2. according to claim 1 it is a kind of calculate charged particle protective layer after power spectrum method, it is characterised in that:The step The differential energy spectrum of track incoming particle is determined in one;Detailed process is:
Satellite orbital altitude and inclination of satellite orbit are determined according to satellite orbit;
Terrestrial radiation having electronic, radiation belt of the earth proton, the sun are determined respectively according to satellite orbital altitude and inclination of satellite orbit Ultra rays proton, solar cosmic ray ion and galactic comic ray ion;
AE8 or AE9 models are selected according to terrestrial radiation having electronic, obtain the differential energy spectrum f (E of track incident electron1′);
AP8 or AP9 models are selected according to radiation belt of the earth proton, obtain the differential energy spectrum f (E of track incident proton2′);
ESP models are selected according to solar cosmic ray proton, obtain the differential energy spectrum f (E of track incident proton3′);
CREME96 models are selected according to solar cosmic ray ion, obtain the differential energy spectrum f (E of track incident ion4′);
CREME96 models are selected according to galactic comic ray ion, obtain the differential energy spectrum f (E of track incident ion5′);
Particle includes ion, electronics and proton.
3. according to claim 2 it is a kind of calculate charged particle protective layer after power spectrum method, it is characterised in that:The step If incoming particle is electronics in four, according to the protective materials thickness t that step 2 determines, step 1 middle orbit incident electron is calculated Differential energy spectrum f (E1') in ENERGY E1' the relationship with range and protective materials thickness, specific formula are:
Wherein, Ri(E1) it is residual range of the radiation belt of the earth incident electron in protective materials, pass through EGS, GEANT software meter It obtains;RiFor residual range;E1For energy.
4. according to claim 3 it is a kind of calculate charged particle protective layer after power spectrum method, it is characterised in that:The step If 1) incoming particle is proton in five, the incoming particle that the protective materials thickness t determined according to step 2 and step 3 determine Incident angle α calculates separately the differential energy spectrum f (E of step 1 middle orbit incident proton2') in ENERGY E2' and f (E3') in energy E3' the relationship with range and protective materials thickness, specific formula are:
Wherein, Ri(E2) it is residual range of the radiation belt of the earth incident proton in protective materials, pass through SRIM, GEANT software meter It obtains;RiFor residual range;E2For energy;
Ri(E3) it is residual range of the solar cosmic ray incident proton in protective materials, it is calculated by SRIM, GEANT software It arrives;RiFor residual range;E3For energy;
2) if incoming particle is ion, the incoming particle that the protective materials thickness t determined according to step 2 and step 3 determine Incident angle α calculates separately the differential energy spectrum f (E of step 1 middle orbit incident ion4') in ENERGY E4' and f (E5') in energy E5' the relationship with range and protective materials thickness, specific formula are:
Wherein, Ri(E4) it is residual range of the solar cosmic ray incident ion in protective materials, pass through SRIM, GEANT software meter It obtains;RiFor residual range;E4For energy;
Ri(E5) it is residual range of the galactic comic ray incident ion in protective materials, it is calculated by SRIM, GEANT software It arrives;RiFor residual range;E5For energy.
5. according to claim 4 it is a kind of calculate charged particle protective layer after power spectrum method, it is characterised in that:The step The inelastic cross section of incoming particle and protective materials is calculated in six, specific formula is:
Wherein, N is avogadros constant, and A is protective materials atomic number, and M is protective materials atomic weight.
6. according to claim 5 it is a kind of calculate charged particle protective layer after power spectrum method, it is characterised in that:The step If incoming particle is electronics in seven, according to the differential energy spectrum f (E of step 11'), the protective materials thickness t of step 2, step 4 Obtained ENERGY E1' inelastic cross section the σ with step 6, calculates the differential energy spectrum after radiation belt of the earth electronic protective layers, specifically Formula is:
Wherein, S (E1') and S (E1) be energy be E1' and E1Stopping power of the terrestrial radiation having electronic in protective materials, lead to EGS, GEANT software is crossed to be calculated.
7. according to claim 6 it is a kind of calculate charged particle protective layer after power spectrum method, it is characterised in that:The step If 1) incoming particle is proton in eight, according to the differential energy spectrum f (E of step 12'), the protective materials thickness t of step 2, step The ENERGY E that three incident angle α, step 5 obtain2' inelastic cross section the σ with step 6 calculates the protection of radiation belt of the earth proton Differential energy spectrum after layer, specific formula are:
Wherein, S (E2') and S (E2) be energy be E2' and E2Stopping power of the radiation belt of the earth proton in protective materials, lead to SRIM, GEANT software is crossed to be calculated;
According to the differential energy spectrum f (E of step 13'), the protective materials thickness t of step 2, the incident angle α of step 3, step 5 Obtained ENERGY E3' inelastic cross section the σ with step 6, calculates the differential energy spectrum after solar cosmic ray proton protective layer, specifically Formula is:
Wherein, S (E3') and S (E3) be energy be E3' and E3Stopping power of the solar cosmic ray proton in protective materials, lead to SRIM, GEANT software is crossed to be calculated;
2) if incoming particle is ion, according to the differential energy spectrum f (E of step 14'), the protective materials thickness t of step 2, step 3 Incident angle α, the obtained ENERGY E of step 54' inelastic cross section the σ with step 6, calculates solar cosmic ray ion protective layer Differential energy spectrum afterwards, specific formula are:
Wherein, S (E4') and S (E4) be energy be E4' and E4Stopping power of the solar cosmic ray ion in protective materials, lead to SRIM, GEANT software is crossed to be calculated;
According to the differential energy spectrum f (E of step 15'), the protective materials thickness t of step 2, the incident angle α of step 3, step 5 Obtained ENERGY E5' inelastic cross section the σ with step 6, calculates the differential energy spectrum after galactic comic ray ion protective layer, specifically Formula is:
Wherein, S (E5') and S (E5) be energy be E5' and E5Stopping power of the galactic comic ray ion in protective materials, lead to SRIM, GEANT software is crossed to be calculated.
CN201810136597.0A 2018-02-09 2018-02-09 Method for calculating energy spectrum of charged particle after protective layer Active CN108345026B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810136597.0A CN108345026B (en) 2018-02-09 2018-02-09 Method for calculating energy spectrum of charged particle after protective layer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810136597.0A CN108345026B (en) 2018-02-09 2018-02-09 Method for calculating energy spectrum of charged particle after protective layer

Publications (2)

Publication Number Publication Date
CN108345026A true CN108345026A (en) 2018-07-31
CN108345026B CN108345026B (en) 2021-06-15

Family

ID=62960144

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810136597.0A Active CN108345026B (en) 2018-02-09 2018-02-09 Method for calculating energy spectrum of charged particle after protective layer

Country Status (1)

Country Link
CN (1) CN108345026B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102426382A (en) * 2011-09-01 2012-04-25 舟山雷大电子科技有限公司 Measuring equipment of electron accelerator output beam parameter
CN103116176A (en) * 2013-01-23 2013-05-22 哈尔滨工业大学 Ground equivalent fluence calculating method for electronic component charged particle irradiation effect
CN103983874A (en) * 2014-05-16 2014-08-13 中国科学院新疆理化技术研究所 Electron accelerator and method for achieving electron beam low-fluence environment
WO2014184589A1 (en) * 2013-05-16 2014-11-20 Ibex Innovations Ltd X-ray detector apparatus
WO2015050141A1 (en) * 2013-10-03 2015-04-09 日立アロカメディカル株式会社 Semiconductor radiation detector, nuclear medicine diagnostic device using same, and method for producing semiconductor radiation detector
CN104570044A (en) * 2014-12-10 2015-04-29 北京永新医疗设备有限公司 Real-time detecting system and method for radioactive materials carried by personnel
CN106442295A (en) * 2016-10-19 2017-02-22 哈尔滨工业大学 Ground simulation test method of a spatial comprehensive environment effect of polymeric material used for spacecraft
CN106483061A (en) * 2016-10-19 2017-03-08 哈尔滨工业大学 A kind of method that cascade and point defect are produced based on charged particle irradiation

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102426382A (en) * 2011-09-01 2012-04-25 舟山雷大电子科技有限公司 Measuring equipment of electron accelerator output beam parameter
CN103116176A (en) * 2013-01-23 2013-05-22 哈尔滨工业大学 Ground equivalent fluence calculating method for electronic component charged particle irradiation effect
WO2014184589A1 (en) * 2013-05-16 2014-11-20 Ibex Innovations Ltd X-ray detector apparatus
WO2015050141A1 (en) * 2013-10-03 2015-04-09 日立アロカメディカル株式会社 Semiconductor radiation detector, nuclear medicine diagnostic device using same, and method for producing semiconductor radiation detector
CN103983874A (en) * 2014-05-16 2014-08-13 中国科学院新疆理化技术研究所 Electron accelerator and method for achieving electron beam low-fluence environment
CN104570044A (en) * 2014-12-10 2015-04-29 北京永新医疗设备有限公司 Real-time detecting system and method for radioactive materials carried by personnel
CN106442295A (en) * 2016-10-19 2017-02-22 哈尔滨工业大学 Ground simulation test method of a spatial comprehensive environment effect of polymeric material used for spacecraft
CN106483061A (en) * 2016-10-19 2017-03-08 哈尔滨工业大学 A kind of method that cascade and point defect are produced based on charged particle irradiation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
方绍辉: "叠层式固体径迹探测器测量快中子能谱影响因素的研究", 《中国核物理学会 会议论文集》 *
李兴冀: "电子辐照LDPE_MWCNTs复合材料的熔融与结晶行为", 《物理学报》 *

Also Published As

Publication number Publication date
CN108345026B (en) 2021-06-15

Similar Documents

Publication Publication Date Title
Messenger et al. Modeling solar cell degradation in space: A comparison of the NRL displacement damage dose and the JPL equivalent fluence approaches
Picozza et al. Scientific goals and in-orbit performance of the high-energy particle detector on board the CSES
Ya'acob et al. Mitigation of space radiation effects on satellites at Low Earth Orbit (LEO)
CN106568702A (en) Ground simulation test method for space comprehensive environmental effect of exposed material for medium and high orbit spacecrafts
Conrad et al. Collateral damage to satellites from an EMP attack
Summers et al. Low energy proton-induced displacement damage in shielded GaAs solar cells in space
Plis et al. Effect of simulated GEO environment on the properties of solar panel coverglasses
CN115186468A (en) Real-time quantitative representation method, device and system for space comprehensive environment
Tverskaya Diagnostics of the magnetosphere based on the outer belt relativistic electrons and penetration of solar protons: A review
CN108345026A (en) A method of calculating power spectrum after charged particle protective layer
Okamoto et al. Gamma-ray irradiation effects on CdTe solar cell dosimeter
Durcanin et al. The definition of the low Earth orbital environment and its effect on thermal control materials
CN113740241A (en) Space environment synergistic effect cutting method of nonmetallic material for spacecraft
Kim et al. Initial operation and preliminary results of the instrument for the study of stable/storm-time space (ISSS) on board the next generation small satellite-1 (NEXTSat-1)
Altıngün et al. A simulation study for the expected performance of Sharjah-Sat-1 payload improved X-Ray Detector (iXRD) in the orbital background radiation
Cooper et al. Local time asymmetry of drift shells for energetic electrons in the middle magnetosphere of Saturn
Imaizumi et al. Displacement damage dose analysis of the output characteristics of In0. 5Ga0. 5P and Cu (In, Ga)(S, Se) 2 solar cells irradiated with alpha ray simulated helium ions
Yamaoka et al. Solar Neutron and Gamma-ray Monitor on the ChubuSat-2 Satellite
Palmer et al. Gamma-ray continuum and line observations of SN 1987A
Reitz Cosmic rays, nuclear reactors, and manned space systems
Claeys et al. Radiation environments and component selection strategy
Okamoto et al. Effects of He ion and electron beam irradiation on a CdTe solar cell dosimeter in a substrate configuration
Brucker et al. Successful large-scale use of CMOS devices on spacecraft traveling through intense radiation belts
Panasyuk Empirical models of terrestrial trapped radiation
Fischi et al. Digital Tracker Array Radiation Test Report

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