Radiation effect shielding calculation method and device based on three-dimensional spacecraft model

Abstract

A radiation effect shielding calculation method based on a three-dimensional spacecraft model is provided. The method includes obtaining the three-dimensional spacecraft model; dividing a 4π space irradiation environment into meshes to obtain space irradiation vectors corresponding to the meshes respectively; setting one or more analysis points at each target model component; inputting the space irradiation vectors to the three-dimensional spacecraft model to obtain a three-dimensional distribution of actual shielding at each analysis point; obtaining a three-dimensional distribution of an equivalent aluminum shielding thicknesses according to equivalent thickness conversion for different materials; calculating residual irradiation effect values of the three-dimensional distribution of the equivalent aluminum shielding thicknesses of each analysis point to obtain corresponding detailed irradiation data; and performing information post processing and integrating in all directions to obtain a corresponding total radiation effect value.

Claims

1. A radiation effect shielding calculation method based on a three-dimensional spacecraft model, comprising: step 1, obtaining the three-dimensional spacecraft model; step 2, dividing a 4π space irradiation environment where the three-dimensional spacecraft model is located, into meshes to obtain space irradiation vectors corresponding to the meshes respectively and taking the space irradiation vectors as whole-space irradiation vectors; wherein the step 2 comprises: selecting a mesh division from a group consisting of an equal-angle mesh division with equal angles in longitude and latitude directions, an equal-area mesh division with equal angles in a longitude direction and equal heights in a latitude direction and an intelligent mesh division with equal angles in a longitude direction and equal heights in a latitude direction as well as mesh optimization in two poles; setting numbers of longitudinal divisions and latitudinal divisions in the mesh division to divide the 4π space irradiation environment where the three-dimensional spacecraft model is located, into the meshes; setting one or more analysis points at each of target model components in the three-dimensional spacecraft model by one of importing position coordinates of the one or more analysis points in batches and using a mouse to click positions at the target model components; using a ray tracing method, using each of the analysis points of the target model components as a starting point, assigning the meshes irradiation vectors that pass through geometrical centers of the meshes corresponding thereto and cover model entities, with each the analysis point as a center and outward directions, of the three-dimensional spacecraft model, to obtain the space irradiation vectors corresponding to the meshes respectively, and taking the space irradiation vectors as the whole-space irradiation vectors corresponding to each the analysis point; step 3, inputting the whole-space irradiation vectors corresponding to each the analysis point to the three-dimensional spacecraft model to perform shielding calculation by a shielding thickness calculation method, selecting one of a slant method, a norm method and a weight method according to a type of radiation particles to obtain a three-dimensional distribution of actual shielding thicknesses at each the analysis point, and performing an equivalent thickness conversion for different materials to obtain a three-dimensional distribution of equivalent aluminum shielding thicknesses including different directions and solid angles of each the analysis point; wherein the step 3 comprises: inputting the whole-space irradiation vectors corresponding to each the analysis point to the three-dimensional spacecraft model to perform the shielding calculation by the shielding thickness calculation method, obtaining position coordinates (x.sub.i, y.sub.i, z.sub.i) of an intersection of an ith space irradiation vector and the target model component in the three-dimensional spacecraft model, calculating a distribution d={d.sub.1, d.sub.2, d.sub.3, . . . , d.sub.n} of actual intersection thicknesses of actual rays in directions of the space irradiation vectors and the target model components in the three-dimensional spacecraft model, wherein each the actual intersection thickness is calculated by the following formula: $d_{i} = \underset{j = 1}{\overset{k}{.Math.}} \sqrt{(x_{j 1}^{2} - x_{j 2}^{2})}$ wherein d.sub.i represents the actual intersection thickness of the actual ray in the direction of the ith space irradiation vector and the target model components in the three-dimensional spacecraft model; k represents a quantity of the target model components, which intersect with the ith space irradiation vector; x.sub.j1 represents the position coordinates (x.sub.i1, y.sub.i1, z.sub.i1) of the intersection of the actual ray incident on a jth target model component; x.sub.j2 represents the position coordinates (x.sub.j2, y.sub.j2, z.sub.j2) of the intersection of the actual ray emergent on the jth target model component; selecting the one of the slant method, the norm method and the weight method according to the type of radiation particles to work out the distribution d={d.sub.1, d.sub.2, d.sub.3, . . . , d.sub.n} of actual intersection thickness to obtain the three-dimensional distribution of actual shielding thicknesses at each the analysis point; wherein the slant method comprises: obtaining a three-dimensional distribution d.sub.z={d.sub.1, d.sub.2, d.sub.3, . . . , d.sub.n} of actual shielding thicknesses at each the analysis point by the actual rays propagating along a straight line in the target model components; wherein the norm method comprises: obtaining a three-dimensional distribution d.sub.f={d.sub.1×cos α.sub.1, d.sub.2×cos α.sub.2, d.sub.3×cos α.sub.3, . . . , d.sub.n×cos α.sub.i} of actual shielding thicknesses at each the analysis point by the actual rays propagating along a normal direction of an incident point into a surface when the actual rays shoot into the target model components; wherein α.sub.i represents an incident angle, and i=1, . . . , n; wherein the weight method comprises: obtaining a three-dimensional distribution $d_{q} = {d_{1} \times \frac{1}{2} (\cos α_{11} + \cos α_{12}), d_{2} \times \frac{1}{2} (\cos α_{21} + \cos α_{22}), d_{3} \times \frac{1}{2} (\cos α_{31} + \cos α_{32}), .Math., d_{n} \times \frac{1}{2} (\cos α_{i 1} + \cos α_{i 2})}$ angle and α.sub.i2 represents an emergent angle, and i=1, . . . , n; converting the actual shielding thickness at each the analysis point to the equivalent aluminum shielding thicknesses according to material properties of model components to as the equivalent aluminum shielding thicknesses of the space irradiation vectors respectively, and combining the equivalent aluminum shielding thicknesses of the space irradiation vectors corresponding to each the analysis point to obtain the three-dimensional distribution of the equivalent aluminum shielding thicknesses including the different directions and the solid angles of each the analysis point; step 4, calculating residual irradiation effect values of the three-dimensional distribution of the equivalent aluminum shielding thicknesses of each the analysis point according to an acquired radiation effect curve to obtain detailed irradiation data of a three-dimensional distribution of irradiation effect values of each the analysis point; and step 5, performing information post processing on the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of each the analysis point to generate and display a two-dimensional cloud image, a three-dimensional cloud image, a three-dimensional ray diagram and a detailed three-dimensional distribution statement of radiation effect, and integrating the detailed radiation data of the three-dimensional distribution of the irradiation effect values of each the analysis point in all directions to obtain a total radiation effect value of each the analysis point to evaluate an irradiation effect of each the analysis point.

2. The radiation effect shielding calculation method based on the three-dimensional spacecraft model according to claim 1, wherein the step 1 comprises: filtering and simplifying the target model components of the three-dimensional spacecraft model and obtaining the target model components after filtering and simplifying, comprising: obtaining a three-dimensional spacecraft model file in one of a step format file and an iges format file; extracting information of the three-dimensional spacecraft model file to obtain model information, wherein the model information comprises units, densities, volumes, masses and names of the target model components; checking unit consistency of the target model components of the model information, and performing interference check to determine whether the model entities overlap; modifying the model information to make the units of the target model components be consistent with one another, and the model entities do not overlap; presetting density threshold values, volume threshold values, mass threshold values and name threshold values, and determining whether the densities, the units, the volumes, the masses and the names are coved in presetting threshold value ranges respectively; adjusting one of the densities, the units, the volumes, the masses and the names until that the densities, the units, the volumes, the masses and the names in the model information are within the presetting threshold value ranges respectively, when the one of the densities, the units, the volumes, the masses and the names in the model information exceeds the presetting threshold value ranges respectively; completing the filtering and simplifying the target model components of the three-dimensional spacecraft model when the densities, the units, the volumes, the masses and the names in the model information are within the presetting threshold value ranges respectively; setting specific material properties of the target model components of the three-dimensional spacecraft model by one of material leading-in and interface setting to determine the material properties of the model components; and obtaining the three-dimensional spacecraft model.

3. The radiation effect shielding calculation method based on the three-dimensional spacecraft model according to claim 1, wherein the step 4 comprises: step 4.1, reading the acquired radiation effect curve and extracting radiation effect depths; step 4.2, comparing the three-dimensional distribution of the equivalent aluminum shielding thicknesses of each the analysis point with the radiation effect depths, performing interpolation calculation on the irradiation effect values corresponding to the radiation effect depths, taking weight according to a proportion of a mesh area corresponding to the space irradiation vectors in a total mesh area to obtain radiation effect values of the space radiation vector and using the radiation effect values as the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of the analysis point; and step 4.3, repeating the step 4.1 and the step 4.2 to obtain the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of each the analysis point.

4. The radiation effect shielding calculation method based on the three-dimensional spacecraft model according to claim 1, wherein the step 5 comprises: displaying the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of each the analysis points in a form of color temperature image to obtain the two-dimensional cloud image and the three-dimensional cloud image; displaying the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of each the analysis point in a form of ray colored by color temperature to create the three-dimensional ray diagram; highlighting and cruising displaying the three-dimensional ray diagram according to a designed longitude and latitude range to display the equivalent aluminum shielding thickness and radiation effect values in the designed longitude and latitude range correspondingly; generating the detailed three-dimensional distribution statement of radiation effect according to positions of the meshes, the equivalent aluminum shielding thicknesses, and the radiation effect values; and integrating the irradiation effect values of each the analysis point in all directions to obtain the total irradiation effect value of each the analysis point to evaluate the irradiation effect of each the analysis point.

5. A radiation effect shielding calculation device based on a three-dimensional spacecraft model, comprising: a processor; and a storage, connected with the processor, wherein the storage is stored with computer programs; and wherein the processor is configured to, when the computer programs are executed by the processor, implement the radiation effect shielding calculation method based on the three-dimensional spacecraft model according to claim 1.

Description

BRIEF DESCRIPTION OF A DRAWING

(1) FIG. 1 is a flow diagram of a radiation effect shielding calculation method based on a three-dimensional spacecraft model according to an embodiment of the disclosure.

(2) FIG. 2 is a schematic structural diagram of the three-dimensional spacecraft model according to the embodiment of the disclosure.

(3) FIG. 3 is a schematic structural diagram of full space irradiation vectors according to an embodiment of the disclosure.

(4) FIG. 4A is a schematic structural diagram of an equal-angle mesh division according to an embodiment of the disclosure.

(5) FIG. 4B is a schematic structural diagram of an equal-area mesh division according to an embodiment of the disclosure.

(6) FIG. 4C is a schematic structural diagram of an intelligent mesh division according to an embodiment of the disclosure.

(7) FIG. 5 is a schematic structural diagram of analysis points in the three-dimensional spacecraft model according to an embodiment of the disclosure.

(8) FIG. 6 is a schematic structural diagram of a two-dimensional cloud image of a shielding thickness according to an embodiment of the disclosure.

(9) FIG. 7 is a schematic structural diagram of a three-dimensional distribution of actual shielding thicknesses according to an embodiment of the disclosure.

(10) FIG. 8 is a schematic structural diagram of a radiation effect curve according to the embodiment of the disclosure.

(11) FIG. 9 is a schematic structural diagram of a three-dimensional ray diagram according to an embodiment of the disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(12) The disclosure is further described combined with the attached drawing and embodiments.

(13) The disclosure provides a radiation effect shielding calculation method based on a three-dimensional spacecraft model. The method realizes three-dimensional model reading, model checking and modification, model simplification and material setting, radiation effect depth curve import, mesh division, three-dimensional shielding thickness analysis, three-dimensional dose analysis and three-dimensional dose display.

(14) Reference to FIG. 1, the method includes:

(15) step 1, obtaining a three-dimensional spacecraft model;

(16) intelligently checking and modifying the three-dimensional spacecraft model; filtering and simplifying the three-dimensional spacecraft model; and setting material properties;

(17) in an illustrated embodiment of the disclosure, obtaining a three-dimensional spacecraft model file in step format or iges format, which is based on a free CAD software;

(18) extracting information of the three-dimensional spacecraft model file to obtain model information, which includes densities, units, volumes, masses and names of model components;

(19) intelligently checking density accuracy and unit consistency of the model components of the model information; performing interference check to determine whether model entities overlap; and intelligently modifying the model information to make the densities of the model components be accurate, the units of the model components be consistent, and the model entities do not overlap;

(20) according to needs of a user, respectively presetting threshold values for density, volume, mass and name; and determining whether the densities, the units, the volumes, the masses and the names are coved in presetting threshold value ranges respectively;

(21) if any one of the densities, the units, the volumes, the masses and the names in the model information exceeds the presetting threshold value ranges respectively, correspondingly adjusting until that the units, the volumes, the masses and the names in the model information are within the presetting threshold value ranges respectively;

(22) if the densities, the units, the volumes, the masses and the names in the model information are within the presetting threshold value ranges respectively, completing the filtering and simplifying the model components of the three-dimensional spacecraft model;

(23) setting specific material properties of the model components of the three-dimensional spacecraft model by material leading-in or interface setting to determine the material properties of the model components;

(24) thereby obtaining the three-dimensional spacecraft model;

(25) step 2, dividing a 4π space irradiation environment, where the three-dimensional spacecraft model is located, into meshes to obtain space irradiation vectors corresponding to the meshes respectively; and taking the space irradiation vectors as whole-space irradiation vectors;

(26) setting one or more analysis points at each target model component in the three-dimensional spacecraft model;

(27) in an illustrated embodiment of the disclosure, according to needs of calculation accuracy and efficiency, selecting an equal-angle mesh division with equal angles in longitude and latitude directions or an equal-area mesh division with equal angles in a longitude direction and equal heights in a latitude direction or an intelligent mesh division with equal angles in a longitude direction and equal heights in a latitude direction as well as mesh optimization in two poles, combined with setting numbers of longitudinal divisions and latitudinal divisions in the mesh division to divide the 4π space irradiation environment, where the three-dimensional spacecraft model is located, into the meshes;

(28) setting the one or more analysis points at the each target model component in the three-dimensional spacecraft model by importing position coordinates of the one or more analysis points in batches or using a mouse to click positions at the target model components;

(29) using a ray tracing method, using each analysis point of the target model components as a starting point, assigning the meshes irradiation vectors that pass through geometrical centers of the meshes corresponding thereto and cover the model entities, with each analysis point as a center and outward directions to obtain the space irradiation vectors corresponding to the meshes respectively; taking the space irradiation vectors as the whole-space irradiation vectors corresponding to each analysis point; marking a position of each mesh while the mesh division to locate the position of the each mesh for subsequent calculation;

(30) step 3, inputting the whole-space irradiation vectors corresponding to each analysis point to the three-dimensional spacecraft model to perform shielding calculation by a shielding thickness calculation method; selecting a slant method or a norm method or a weight method according to a type of radiation particles to obtain a three-dimensional distribution of actual shielding thicknesses at each analysis point; performing an equivalent thickness conversion for different materials to obtain a three-dimensional distribution of equivalent aluminum shielding thicknesses including different directions and solid angles of each analysis point; in an illustrated embodiment of the disclosure, inputting the whole-space irradiation vectors corresponding to each analysis point to the three-dimensional spacecraft model to performing the shielding calculation by the shielding thickness calculation method, obtaining position coordinates (x.sub.i, y.sub.i, z.sub.i) of an intersection of an ith space irradiation vector and the target model component in the three-dimensional spacecraft model;

(31) calculating a distribution d={d.sub.1, d.sub.2, d.sub.3, . . . , d.sub.n} of actual intersection thicknesses of actual rays in directions of the space irradiation vectors and the target model components in the three-dimensional spacecraft model, which is calculated by the following formula:
d.sub.i=Σ.sub.j=1.sup.k√{square root over ((x.sub.j1.sup.2−x.sub.j2.sup.2))}

(32) in an illustrated embodiment of the disclosure, d.sub.i representing the actual intersection thickness of the actual ray in the direction of the ith space irradiation vector and the target model components in the three-dimensional spacecraft model; k representing a quantity of the target model components, which intersect with the ith space irradiation vector; x.sub.j1 representing the position coordinates (x.sub.i1, y.sub.i1, z.sub.i1) of the intersection of the actual incident ray on the jth target model component; x.sub.j2 represents the position coordinates (x.sub.j2, y.sub.j2, z.sub.j2) of the intersection of the actual emergent ray on a jth target model component;

(33) selecting the slant method or the norm method or the weight method according to the type of radiation particles to work out the distribution d={d.sub.1, d.sub.2, d.sub.3, . . . , d.sub.n} of actual intersection thickness; thereby obtaining the three-dimensional distribution of actual shielding thicknesses at each analysis point;

(34) in an illustrated embodiment of the disclosure, the slant method including: obtaining a three-dimensional distribution d.sub.z={d.sub.1, d.sub.2, d.sub.3, . . . d.sub.n} of actual shielding thicknesses at each analysis point by the actual rays propagating along a straight line in the target model components;

(35) the norm method including: obtaining a three-dimensional distribution d.sub.f={d.sub.1×cos α.sub.1, d.sub.2×cos α.sub.2, d.sub.3×cos α.sub.3, d.sub.n×cos α.sub.i} of actual shielding thicknesses at each analysis point by the actual rays propagating along a normal direction of an incident point into a surface when the actual rays shoot into the target model components; of which α.sub.i represents an incident angle, and i=1, . . . , n;

(36) the weight method including: obtaining a three-dimensional distribution

(37) $d_{q} = {d_{1} \times \frac{1}{2} (\cos α_{11} + \cos α_{12}), d_{2} \times \frac{1}{2} (\cos α_{21} + \cos α_{22}), d_{3} \times \frac{1}{2} (\cos_{31} + \cos α_{32}), .Math., d_{n} \times \frac{1}{2} (\cos α_{i 1} + \cos α_{i 2})}$
of actual shielding thicknesses at each analysis point; of which α.sub.i1 represents an incident angle and α.sub.i2 represents an emergent angle, and i=1, . . . , n;

(38) converting the actual shielding thicknesses at each analysis point to the equivalent aluminum shielding thicknesses according to material properties of model components to as the equivalent aluminum shielding thicknesses of the space irradiation vectors respectively; and combining the equivalent aluminum shielding thicknesses of the space irradiation vectors corresponding to each analysis point to obtain the three-dimensional distribution of the equivalent aluminum shielding thicknesses including the different directions and the solid angles of each analysis point;

(39) step 4, calculating residual irradiation effect values of the three-dimensional distribution of the equivalent aluminum shielding thicknesses of each analysis point according to an acquired radiation effect curve to obtain detailed irradiation data of a three-dimensional distribution of irradiation effect values of each analysis point;

(40) in an illustrated embodiment of the disclosure, reading the acquired radiation effect curve and extracting radiation effect depths;

(41) comparing the three-dimensional distribution of the equivalent aluminum shielding thicknesses of each analysis point with the radiation effect depths; performing interpolation calculation on radiation effect values corresponding to the radiation effect depths; taking weight according to a proportion of a mesh area corresponding to the space irradiation vectors in a total mesh area to obtain radiation effect values of the space radiation vector; and using the radiation effect values as the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of the analysis point;

(42) repeating the above process to obtain the detailed irradiation data of the distribution of the three-dimensional irradiation effect value at each analysis point;

(43) step 5, performing information post processing on the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of each analysis point to generate and display a two-dimensional cloud image, a three-dimensional cloud image, a three-dimensional ray diagram and a detailed three-dimensional distribution statement of radiation effect; and integrating the detailed radiation data of the three-dimensional distribution of the irradiation effect values of each analysis point in all directions to obtain a total radiation effect value of each analysis point to evaluate an irradiation effect of each analysis point;

(44) in an illustrated embodiment of the disclosure, displaying the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of each analysis point in color temperature images to obtain the two-dimensional cloud image and the three-dimensional cloud image;

(45) displaying the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of each analysis point in color-temperature colored rays to create the three-dimensional ray diagram;

(46) highlighting and cruising displaying the three-dimensional ray diagram according to a designed longitude and latitude range to display the equivalent aluminum shielding thickness and radiation effect values in the designed longitude and latitude range correspondingly;

(47) generating the detailed three-dimensional distribution statement of radiation effect according to positions of the meshes, the equivalent aluminum shielding thicknesses, and the radiation effect values;

(48) integrating the irradiation effect values of each analysis point in all directions to obtain the total irradiation effect value of each analysis point to evaluate the irradiation effect of each analysis point.

(49) The disclosure also provides a radiation effect shielding calculation apparatus based on a three-dimensional spacecraft model, including:

(50) a model acquisition module, which is used to obtain the three-dimensional spacecraft model, intelligently check and modify the three-dimensional spacecraft model, filter and simplify the three-dimensional spacecraft model and set material properties;

(51) a space division module, which is used to divide a 4π space irradiation environment, where the three-dimensional spacecraft model is located, into meshes to obtain space irradiation vectors corresponding to the meshes respectively, take the space irradiation vectors as whole-space irradiation vectors, and set one or more analysis points at each target model component in the three-dimensional spacecraft model;

(52) a shielding thickness acquisition module, which is used to input the whole-space irradiation vectors to the three-dimensional spacecraft model to perform shielding calculation by a shielding thickness calculation method, select one of a slant method, a norm method and a weight method according to a type of radiation particles to obtain a three-dimensional distribution of actual shielding thicknesses at each analysis point, and performing an equivalent thickness conversion for different materials to obtain a three-dimensional distribution of equivalent aluminum shielding thicknesses including different directions and solid angles of each analysis point;

(53) an irradiation effect calculation module, which is used to calculate residual irradiation effect values of the three-dimensional distribution of equivalent aluminum shielding thicknesses of each analysis point according to an acquired radiation effect curve to obtain detailed irradiation data of a three-dimensional distribution of irradiation effect values of each analysis point;

(54) an evaluation module, configured to perform information post processing on the detailed irradiation data of the three-dimensional distribution of the irradiation effect values of each analysis point to generate and display a two-dimensional cloud image, a three-dimensional cloud image, a three-dimensional ray diagram and a detailed three-dimensional distribution statement of radiation effect, and integrate the detailed radiation data of the three-dimensional distribution of the irradiation effect values of each analysis point in all directions to obtain a total irradiation effect value of each analysis point, so as to evaluate an irradiation effect of each analysis point.

(55) The disclosure also provides a radiation effect shielding calculation device based on a three-dimensional spacecraft model, which includes a processor and a storage connected with the processor.

(56) The storage is used to store computer programs.

(57) The processor is used to implement the radiation effect shielding calculation method based on the three-dimensional spacecraft model, when the computer programs are executed by the processor.

(58) It should be specific that the above embodiments only aim to describe the technical solutions of the disclosure but not to limit. Although the above embodiments describe the disclosure in detail, it should be understood by those skilled in the technical field that modification or equivalent replacement of the technical solutions of the disclosure cannot depart from the spirit and scope of the technical solutions of the disclosure, which shall be covered by the scope of the disclosure.

Radiation effect shielding calculation method and device based on three-dimensional spacecraft model

Assignee

Inventors

Cpc classification

Classification Explorer

G06F2119/08

PHYSICS

Classification Explorer

G06F30/15

PHYSICS

Classification Explorer

G06F30/20

PHYSICS

Classification Explorer

G06F2111/10

PHYSICS

International classification

Classification Explorer

G06F30/20

PHYSICS

Classification Explorer

G06F30/15

PHYSICS

Abstract

Claims

Description