CN103149581A - Computing method of displacement absorbed doses of space three-dimensional structure device - Google Patents
Computing method of displacement absorbed doses of space three-dimensional structure device Download PDFInfo
- Publication number
- CN103149581A CN103149581A CN2013100249236A CN201310024923A CN103149581A CN 103149581 A CN103149581 A CN 103149581A CN 2013100249236 A CN2013100249236 A CN 2013100249236A CN 201310024923 A CN201310024923 A CN 201310024923A CN 103149581 A CN103149581 A CN 103149581A
- Authority
- CN
- China
- Prior art keywords
- space
- dimensional structure
- displacement
- absorbed dose
- computing method
- 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.)
- Pending
Links
Images
Landscapes
- Measurement Of Radiation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a computing method of displacement absorbed doses of a space three-dimensional structure device and relates to the computing method of displacement absorbed doses. The method solves the technical problems, existing in the present displacement absorbed doses, of being high in cost, bad in accuracy. The computing method comprises a step of determining a radiation-receiving energy spectrum of the space three-dimensional structure device; a step of calculating a distribution curve of absorbed doses following with depth of aluminum; a step of ensuring quality of all members of the space three-dimensional structure device; a step of ensuring calculating points of the members; a step of ensuring a plume in a subdivision mode by regarding the calculating points as poles of ball coordination; a step of equating materials on the plume to the aluminum and calculating a displacement absorbed dose valve of the plume; a step of adding up all the displacement absorbed doses values of the plume in order to gain the displacement absorbed doses value of the calculating points; a step of repeating from the fifth step to the seventh step in order to gain the displacement absorbed doses values of all the calculating points of the space three-dimensional structure device. The calculation of the displacement absorbed doses is conducted under an irradiation condition of space rail charge particles.
Description
Technical field
The present invention relates to the computing method of displacement absorbed dose.
Background technology
Accurately calculate spacecraft and device at rail displacement absorbed dose, it is one of important means of estimating spacecraft and device lifetime/reliability, especially to the device of displacement damage sensitivity (as, solar cell, CCD device, electric light and photoelectric device etc.) all the more so.Comparatively deep in the calculating of rail ionization absorbed dose to spacecraft and device in the world at present.For example, be based on the earliest one dimension Monte Carlo method and develop and be applied to SHIELDOSE on engineering
[1]Program.This program can the Computing One-Dimensional housing, the target absorbed dose of unlimited thickness of slab and the protection of limited thickness of slab, and wherein the material of protective layer mainly is thought of as Al, and the material of target can be Si, SiO
2, H
2O, Al etc.This program convenience of calculation, but can not carry out the calculating of 3 D complex structure absorbed dose, can not the displacement calculating absorbed dose.Subsequently, developed NOVICE, MCNP and the Geant4 supervisor of developing based on three-dimensional Monte Carlo method.The computing method of these programs can consider that the ionization absorbed dose of labyrinth calculates, but because the convergence of Monte Carlo computing method is poor especially, if so structure more complicated, quite consuming time, sometimes may all can not get result preferably by some months, this is very unrealistic on scientific research, especially engineering.These years recently, the U.S. and European countries have developed a kind of new ionization absorbed dose computing method---fan-shaped Meshing Method, at first this method is carried out the 3 D complex structure grid and is divided, then divide for every kind of grid, database in conjunction with Monte Carlo method is calculated calculates.The result of calculation of the method is comparatively accurate, and has greatly improved counting yield, is convenient to practical implementation.
Compare with the fast development of ionization absorbed dose computing method, the computing method development of 3 D complex displacement structure absorbed dose is relatively slow.In the calculating of rail displacement absorbed dose, be mainly the displacement damage of estimating spacecraft and device.What be used for the earliest the evaluation displacement damage is equivalent fluence method, and equivalent fluence method is the space service life assessment technique that is used for estimating monocrystalline silicon and unijunction solar cell that proposes in the phase at the beginning of the eighties in last century.This technology can be reinforced the damage of device resist displacement reference frame is provided.Its main shortcoming is: need a large amount of ground simulation tests, and test is to carry out in the mode of ex situ, this just need to provide a large amount of tested device spares, and the test assessment cost is high, and owing to being non-in site measurement, its accuracy is also relatively poor.Development subsequently be that the essence of the NIEL method in NRL laboratory is displacement absorbed dose method, but it does not consider that the displacement absorbed dose is with the distribution of the degree of depth.At present also be not specifically designed in the world the 3 D complex structure devices in the computing method of rail displacement absorbed dose.
Summary of the invention
The method cost that the present invention will solve existing displacement calculating absorbed dose is high, the technical matters of poor accuracy, and the computing method of space three-dimensional structure device displacement absorbed dose are provided.
The computing method of space three-dimensional structure device displacement absorbed dose of the present invention are carried out according to the following steps:
One, press military service track and the military service period of space three-dimensional structure instrument, the power spectrum of the radiation that mensuration space three-dimensional structure device is accepted, the power spectrum of radiation is radiation belt of the earth proton power spectrum, radiation belt of the earth electronic energy spectrum and solar cosmic ray proton power spectrum;
The power spectrum of the radiation that two, obtains according to step 1 utilizes Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve D~t;
Three, determine material, density and the composition of all members in the space three-dimensional structure instrument;
Four, determine the calculation level of member, this calculation level is the arbitrfary point in the space three-dimensional structure instrument;
Five, the limit take calculation level as spherical coordinates, determine spherical coordinate system, spherical coordinate system (r, θ, φ), get θ ∈ 0~π, φ ∈ 0~2 π carries out subdivision for the structure of space three-dimensional structure instrument to θ and φ, determines position and the quantity of polar curve r according to the cell surface of institute's subdivision;
Six, be aluminium with the equivalence of the material on each polar curve, and calculate the thickness t of equivalent aluminium, and according to the D in step 2~t curve, calculate the displacement absorbing agent value on each polar curve;
Seven, the displacement absorbing agent value on all polar curves is added and, obtain the displacement absorbed dose of this calculation level;
Eight, repeating step five to seven, obtain the displacement absorbing agent value of all calculation levels in the space three-dimensional structure instrument.
The present invention finds out the calculation level in the space three-dimensional structure instrument one by one, take each calculation level as the utmost point heart, calculate the equivalent thickness of aluminium of all materials that pass of polar curve of this utmost point heart, obtain the displacement absorbed dose of this utmost point heart place calculation level, method is simple, and accurately, the inventive method can be for all spacecraft military service tracks and complicated topology layout, and the displacement absorbed dose that carries out under the radiation environment of space orbit charged particle calculates.
Description of drawings
Fig. 1 tests the structural representation of 3 d structure model figure in;
Fig. 2 is geostationary orbit in test one, utilizes the AE8 radiation belt electronic model of NASA to calculate acquisition radiation belt electronics differential energy spectrum;
Fig. 3 is geostationary orbit in test one, utilizes the AP8 radiation belt proton model of NASA to calculate acquisition radiation belt proton differential energy spectrum;
Fig. 4 is geostationary orbit in test one, utilizes the JPL91 model to calculate the track solar cosmic ray proton power spectrum that obtains;
Fig. 5 be the test one in the displacement absorbed dose with aluminium depth profile curve D~t.Wherein a is that the displacement absorbed dose of terrestrial radiation having electronic is with the relation curve of aluminium protective thickness; B is that the displacement absorbed dose of radiation belt of the earth proton is with the relation curve of aluminium protective thickness; C is that the displacement absorbed dose of solar cosmic ray proton is with the relation curve of aluminium protective thickness; D is that total displacement absorbed dose is with the relation curve with the aluminium protective thickness;
Embodiment
Embodiment one: present embodiment the computing method of space three-dimensional structure device displacement absorbed dose carry out according to the following steps:
One, press military service track and the military service period of space three-dimensional structure instrument, the power spectrum of the radiation that mensuration space three-dimensional structure device is accepted, the power spectrum of radiation is radiation belt of the earth proton power spectrum, radiation belt of the earth electronic energy spectrum and solar cosmic ray proton power spectrum;
The power spectrum of the radiation that two, obtains according to step 1 utilizes Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve D~t;
Three, determine material, density and the composition of all members in the space three-dimensional structure instrument;
Four, determine the calculation level of member, this calculation level is the arbitrfary point in the space three-dimensional structure instrument;
Five, the limit take calculation level as spherical coordinates, determine spherical coordinate system, spherical coordinate system (r, θ, φ), get θ ∈ 0~π, φ ∈ 0~2 π carries out subdivision for the structure of space three-dimensional structure instrument to θ and φ, determines position and the quantity of polar curve r according to the cell surface of institute's subdivision;
Six, be aluminium with the equivalence of the material on each polar curve, and calculate the thickness t of equivalent aluminium, and according to the D in step 2~t curve, calculate the displacement absorbing agent value on each polar curve;
Seven, the displacement absorbing agent value on all polar curves is added and, obtain the displacement absorbed dose of this calculation level;
Eight, repeating step five to seven, obtain the displacement absorbing agent value of all calculation levels in the space three-dimensional structure instrument.
Present embodiment is found out the calculation level in the space three-dimensional structure instrument one by one, take each calculation level as the utmost point heart, calculate the equivalent thickness of aluminium of all materials that pass of polar curve of this utmost point heart, obtain the displacement absorbed dose of this utmost point heart place calculation level, method is simple, the inventive method can be for all spacecraft military service tracks and complicated topology layout, and the displacement absorbed dose that carries out under the radiation environment of space orbit charged particle calculates.
Embodiment two: what present embodiment was different from embodiment one is: in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision be: interval 0~π of θ is divided into 200~1000 parts, 0~2 π of φ is divided into 200~2000 parts.Other are identical with embodiment one.
Embodiment three: what present embodiment was different from embodiment one is: in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision be: interval 0~π of θ is divided into 250 parts, 0~2 π of φ is divided into 250 parts.Other are identical with embodiment one.
Embodiment four: what present embodiment was different from one of embodiment one to three is: in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision be: interval 0~π of θ is divided into 400 parts, 0~2 π of φ is divided into 300 parts.Other are identical with one of embodiment one to three.
With following verification experimental verification beneficial effect of the present invention:
Test one: the computing method of the space three-dimensional structure device displacement absorbed dose of test one are carried out according to the following steps:
One, press military service track and the military service period of space three-dimensional structure instrument, the power spectrum of the radiation that mensuration space three-dimensional structure device is accepted, the power spectrum of radiation is radiation belt of the earth proton power spectrum, radiation belt of the earth electronic energy spectrum and solar cosmic ray proton power spectrum;
The power spectrum of the radiation that two, obtains according to step 1 utilizes Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve D~t;
Three, determine material, density and the composition of all members in the space three-dimensional structure instrument;
Four, determine the calculation level of member, this calculation level is the arbitrfary point in the space three-dimensional structure instrument;
Five, the limit take calculation level as spherical coordinates, determine spherical coordinate system, spherical coordinate system (r, θ, φ), get θ ∈ 0~π, φ ∈ 0~2 π carries out subdivision for the structure of space three-dimensional structure instrument to θ and φ, and interval 0~π of θ is divided into 300 parts, 0~2 π of φ is divided into 300 parts, determines position and the quantity of polar curve r according to the cell surface of institute's subdivision;
Six, be aluminium with the equivalence of the material on each polar curve, and calculate the thickness t of equivalent aluminium, and according to the D in step 2~t curve, calculate the displacement absorbing agent value on each polar curve;
Seven, the displacement absorbing agent value on all polar curves is added and, obtain the displacement absorbed dose of this calculation level;
Eight, repeating step five to seven, obtain the displacement absorbing agent value of all calculation levels in the space three-dimensional structure instrument.
In this test the schematic diagram of space three-dimensional structure instrument as shown in Figure 1,
Determine orbit parameter according to the space three-dimensional structure device in rail military service task in this test procedure one, comprise orbit altitude, orbit inclination, excentricity etc., can spectrum model be calculated the power spectrum of the radiation that the track differential energy spectrum accepts by track.Fig. 2 is geostationary orbit, utilizes the AE8 radiation belt electronic model of NASA to calculate acquisition radiation belt electronics differential energy spectrum; Fig. 3 is geostationary orbit, utilizes the AP8 radiation belt proton model of NASA to calculate acquisition radiation belt proton differential energy spectrum; Fig. 4 is geostationary orbit, utilizes the JPL91 model to calculate the track solar cosmic ray proton power spectrum that obtains;
Take the power spectrum of step 1 gained as data source, utilize Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve in step 3, concrete outcome as shown in Figure 5.Wherein a is that the displacement absorbed dose of terrestrial radiation having electronic is with the relation curve of aluminium protective thickness; B is that the displacement absorbed dose of radiation belt of the earth proton is with the relation curve of aluminium protective thickness; C is that the displacement absorbed dose of solar cosmic ray proton is with the relation curve of aluminium protective thickness; D is that total displacement absorbed dose is with the relation curve with the aluminium protective thickness; Curve d is 10 at the thickness of protective layer Al
-8~10
-3In time, overlap with curve b, and curve d is 2 * 10 at the thickness of protective layer Al
-3~4 * 10
0In time, overlap with curve a.
For space three-dimensional structure instrument shown in Figure 1, define its inner required calculation level, concrete as " A " point in Fig. 1.As seen from the figure, the safeguard structure placement differences of " A " some surrounding is very large, calculates if use Equivalent A l spherical shell, and its error is very large.And its actual labyrinth situation is considered in the calculating of this test one, and result of calculation is comparatively accurate.For " A " calculation level, carry out three-dimensional ball subdivision.Carry out three-dimensional ball subdivision with around calculation level by position angle and polar angle, calculate the corresponding equivalent material of each little grid direction and density.
With the result of calculation combination in step 2 and six, obtain the displacement absorbed dose of the calculation level that defines in the 3 D complex structure.During concrete calculating, need each little grid integration in the hope of last result of calculation.In this example, result of calculation is 6rad (Si).Thus result also as can be known, its equivalent protective thickness is larger, greater than 1g/cm
2
Test in one be with electron device under geostationary orbit, utilize labyrinth displacement absorbed dose computing method, calculate its inner displacement absorbed dose.For other labyrinth, can realize identical target according to the process of the inventive method.The inventive method can be for all spacecraft military service tracks and complicated topology layout, and the displacement absorbed dose that carries out under the radiation environment of space orbit charged particle calculates.
Claims (4)
1. the computing method of space three-dimensional structure device displacement absorbed dose is characterized in that the method carries out according to the following steps:
One, press military service track and the military service period of space three-dimensional structure instrument, the power spectrum of the radiation that mensuration space three-dimensional structure device is accepted, the power spectrum of radiation is radiation belt of the earth proton power spectrum, radiation belt of the earth electronic energy spectrum and solar cosmic ray proton power spectrum;
The power spectrum of the radiation that two, obtains according to step 1 utilizes Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve D~t;
Three, determine material, density and the composition of all members in the space three-dimensional structure instrument;
Four, determine the calculation level of member, this calculation level is the arbitrfary point in the space three-dimensional structure instrument;
Five, the limit take calculation level as spherical coordinates, determine spherical coordinate system, spherical coordinate system (r, θ, φ), get θ ∈ 0~π, φ ∈ 0~2 π carries out subdivision for the structure of space three-dimensional structure instrument to θ and φ, determines position and the quantity of polar curve r according to the cell surface of institute's subdivision;
Six, be aluminium with the equivalence of the material on each polar curve, and calculate the thickness t of equivalent aluminium, and according to the D in step 2~t curve, calculate the displacement absorbing agent value on each polar curve;
Seven, the displacement absorbing agent value on all polar curves is added and, obtain the displacement absorbed dose of this calculation level;
Eight, repeating step five to seven, obtain the displacement absorbing agent value of all calculation levels in the space three-dimensional structure instrument.
2. the computing method of space three-dimensional structure device displacement absorbed dose according to claim 1, it is characterized in that in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision being: interval 0~π of θ is divided into 200~1000 parts, 0~2 π of φ is divided into 200~2000 parts.
3. the computing method of space three-dimensional structure device displacement absorbed dose according to claim 1, it is characterized in that in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision being: interval 0~π of θ is divided into 250 parts, 0~2 π of φ is divided into 250 parts.
4. the computing method of space three-dimensional structure device displacement absorbed dose according to claim 1, it is characterized in that in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision being: interval 0~π of θ is divided into 400 parts, 0~2 π of φ is divided into 300 parts.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013100249236A CN103149581A (en) | 2013-01-23 | 2013-01-23 | Computing method of displacement absorbed doses of space three-dimensional structure device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013100249236A CN103149581A (en) | 2013-01-23 | 2013-01-23 | Computing method of displacement absorbed doses of space three-dimensional structure device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103149581A true CN103149581A (en) | 2013-06-12 |
Family
ID=48547763
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013100249236A Pending CN103149581A (en) | 2013-01-23 | 2013-01-23 | Computing method of displacement absorbed doses of space three-dimensional structure device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103149581A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106404810A (en) * | 2016-10-19 | 2017-02-15 | 哈尔滨工业大学 | Displacement and radiation damage equivalence evaluation method for aromatic polymer insulating material |
CN111581719A (en) * | 2020-04-28 | 2020-08-25 | 北京天工科仪空间技术有限公司 | Radiation effect calculation method, device and equipment based on spacecraft three-dimensional shielding |
CN113543615A (en) * | 2021-06-29 | 2021-10-22 | 中国科学院长春光学精密机械与物理研究所 | Irradiation protection method for space electronic equipment |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2420861A1 (en) * | 2010-08-20 | 2012-02-22 | Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO | A radiation dose meter for measuring radiation dose in an external magnetic field |
CN102564931A (en) * | 2010-12-31 | 2012-07-11 | 中国航天科技集团公司第五研究院第五一○研究所 | Method for simulation test about comprehensive space radiation effect of surface function material for spacecraft |
-
2013
- 2013-01-23 CN CN2013100249236A patent/CN103149581A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2420861A1 (en) * | 2010-08-20 | 2012-02-22 | Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO | A radiation dose meter for measuring radiation dose in an external magnetic field |
CN102564931A (en) * | 2010-12-31 | 2012-07-11 | 中国航天科技集团公司第五研究院第五一○研究所 | Method for simulation test about comprehensive space radiation effect of surface function material for spacecraft |
Non-Patent Citations (1)
Title |
---|
李兴冀: "星用双极型器件带电粒子辐照效应及损伤机理", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106404810A (en) * | 2016-10-19 | 2017-02-15 | 哈尔滨工业大学 | Displacement and radiation damage equivalence evaluation method for aromatic polymer insulating material |
CN106404810B (en) * | 2016-10-19 | 2019-04-09 | 哈尔滨工业大学 | The shifted radiation equivalent damage evaluation method of aromatic polymer insulating materials |
CN111581719A (en) * | 2020-04-28 | 2020-08-25 | 北京天工科仪空间技术有限公司 | Radiation effect calculation method, device and equipment based on spacecraft three-dimensional shielding |
CN111581719B (en) * | 2020-04-28 | 2024-03-26 | 北京天工科仪空间技术有限公司 | Radiation effect calculation method, device and equipment based on spacecraft three-dimensional shielding |
CN113543615A (en) * | 2021-06-29 | 2021-10-22 | 中国科学院长春光学精密机械与物理研究所 | Irradiation protection method for space electronic equipment |
CN113543615B (en) * | 2021-06-29 | 2022-11-01 | 中国科学院长春光学精密机械与物理研究所 | Irradiation protection method for space electronic equipment |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shikaze et al. | Measurements of 0.2–20 GeV/n cosmic-ray proton and helium spectra from 1997 through 2002 with the BESS spectrometer | |
Mitrofanov et al. | Dynamic albedo of neutrons (DAN) experiment onboard NASA’s Mars Science Laboratory | |
Cârloganu et al. | Towards a muon radiography of the Puy de Dôme | |
Ackermann et al. | Searches for cosmic-ray electron anisotropies with the Fermi Large Area Telescope | |
CN103116176A (en) | Ground equivalent fluence calculating method for electronic component charged particle irradiation effect | |
Schraube et al. | Experimental verification and calculation of aviation route doses | |
Velinov et al. | Model for induced ionization by galactic cosmic rays in the Earth atmosphere and ionosphere | |
Yáñez et al. | A method to measure the integral vertical intensity and angular distribution of atmospheric muons with a stationary plastic scintillator bar detector | |
CN103149581A (en) | Computing method of displacement absorbed doses of space three-dimensional structure device | |
Mitrica et al. | A mobile detector for measurements of the atmospheric muon flux in underground sites | |
Samnøy et al. | Microdosimetry with a 3D silicon on insulator (SOI) detector in a low energy proton beamline | |
CN112051601A (en) | Virtual source principle-based source boundary parameter Monte Carlo inversion technology | |
Gustafsson et al. | PHITS simulations of the Matroshka experiment | |
Tran et al. | A novel silicon microdosimeter using 3D sensitive volumes: Modeling the response in neutron fields typical of aviation | |
CN110057843A (en) | Radioactive source shipping container gamma shielding method for testing performance and system based on G (E) function | |
Kurochkin et al. | Study of the radiation environment caused by galactic cosmic rays at flight altitudes, at the summit of the Zugspitze and at PTB Braunschweig | |
Oláh et al. | Improvement of cosmic-ray muography for Earth sciences and civil engineering | |
CN107545098B (en) | Compound device space continuous spectrum proton displacement damage ground equivalent evaluation method | |
CN105701327A (en) | Method for assessing radiation dose rate of earth-orbiting satellite | |
Panuntun et al. | Validation of HPGe efficiency modeling by MCNP5 for volume source | |
Prokopovich et al. | Comparison of SOI microdosimeter and tissue equivalent proportional counter measurements at the CERF facility | |
Lv et al. | Research on optimization design of radiation dose shield hardening for aerospace components | |
Chen et al. | Monte Carlo simulation for background noise study of the Wudalianchi volcanic density imaging | |
Mitrica | 20 years of cosmic muons research performed in IFIN-HH | |
Mitrica et al. | Studies of the high energy muon flux in Unirea salt mine from Slanic Prahova |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20130612 |