CN108563839B - 一种核设施退役模型程式化仿真方法 - Google Patents

一种核设施退役模型程式化仿真方法 Download PDF

Info

Publication number
CN108563839B
CN108563839B CN201810243219.2A CN201810243219A CN108563839B CN 108563839 B CN108563839 B CN 108563839B CN 201810243219 A CN201810243219 A CN 201810243219A CN 108563839 B CN108563839 B CN 108563839B
Authority
CN
China
Prior art keywords
model
programming
physical object
stylized
constructing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810243219.2A
Other languages
English (en)
Other versions
CN108563839A (zh
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 Engineering University
Original Assignee
Harbin Engineering University
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 Engineering University filed Critical Harbin Engineering University
Priority to CN201810243219.2A priority Critical patent/CN108563839B/zh
Publication of CN108563839A publication Critical patent/CN108563839A/zh
Application granted granted Critical
Publication of CN108563839B publication Critical patent/CN108563839B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明的目的在于提供一种核设施退役模型程式化仿真方法，用3dsMax软件根据确定的核设施参数构建模型文件，并将模型文件保存为3DS格式；导入3DS格式的模型文件，获得模型参数；构建模型包围盒；根据实体对象包围盒构建程式化模型；根据实体对象顶点在程式化模型内的个数确定实体对象外部程式化模型；根据实体对象的模型参数构建实体对象内部程式化模型；使用组合几何方法构建实体对象的程式化模型并写为标准输入文件；将标准输入文件输入到辐射剂量计算程序内获得辐射场剂量分布。本发明实现了高效的将复杂核设施3dsMax模型简化为可计算辐射剂量分布的程式化模型，实现了对内部为空的壳状核设施模型的辐射场剂量分布计算。

Description

一种核设施退役模型程式化仿真方法
技术领域
本发明涉及的是一种仿真方法，具体地说是核设施仿真方法。
背景技术
核设施维修与退役是核设施生命周期的重要环节。在制定核设施维修或退役策略时，为降低辐射对工作人员的伤害，需要准确地了解核设施辐射剂量分布。在核退役辐射仿真中，为获得准确的辐射场剂量分布结果，就必须对复杂的核设施模型进行简化计算。
目前，在核退役虚拟辐射场计算的所有方法中，蒙特卡罗方法与点核积分方法应用广泛，且都可以通过标准输入文件计算辐射场剂量分布。国外已有的软件一般采用直接构建输入文件的方法计算核设施辐射场剂量分布，在处理内部为空的核设施模型时将模型转换为非空的模型，剂量分布计算误差大。
综上所述，开发出一种准确、高效的可将内部为空的核设施模型简化为可计算辐射剂量分布的程式化模型的仿真方法对核退役仿真具有重大的实际意义。
发明内容
本发明的目的在于提供将核设施模型简化为程式化模型，并生成标准输入文件的一种核设施退役模型程式化仿真方法。
本发明的目的是这样实现的：
本发明一种核设施退役模型程式化仿真方法，其特征是：
(1)用3dsMax软件根据确定的核设施参数构建模型文件，并将模型文件保存为3DS格式；
(2)导入3DS格式的模型文件，获得模型参数；
(3)构建模型包围盒；
(4)根据实体对象包围盒构建程式化模型；
(5)根据实体对象顶点在程式化模型内的个数确定实体对象外部程式化模型；
(6)根据实体对象的模型参数构建实体对象内部程式化模型；
(7)使用组合几何方法构建实体对象的程式化模型并写为标准输入文件；
(8)将标准输入文件输入到辐射剂量计算程序内获得辐射场剂量分布。
本发明还可以包括：
1、所述的构建模型包围盒包括：
遍历模型所有顶点，根据模型顶点在x、y、z方向上的最大最小值构建模型场景立方体包围盒，遍历每个实体对象的所有顶点，构建实体对象立方体包围盒，尺寸为△x×△y×△z，包围盒中心点为O(xo,yo,zo)。
2、所述的根据实体对象包围盒构建程式化模型包括，
根据实体对象包围盒的中心点与尺寸，构建基于包围盒的椭球体、柱体和立方体的程式化模型：
(1)椭球体方程为：
Figure BDA0001605845620000021
当△x＝△y＝△z时，方程为球体方程；
(2)柱体的母线平行于坐标轴：
母线平行z轴的柱体方程：
Figure BDA0001605845620000022
当△x＝△y时，方程为圆柱体方程；
(3)立方体方程为：
Figure BDA0001605845620000023
3、所述的根据实体对象顶点在程式化模型内的个数确定实体对象外部程式化模型包括：
计算椭球体和柱体程式化模型包含顶点的个数，选择满足顶点数目要求且顶点最多的模型作为实体对象的外部程式化模型，如果顶点个数相同，则选择椭球体作为实体对象的外部程式化模型，如果椭球体和柱体程式化模型都不满足数目要求，则选择立方体作为实体对象的外部程式化模型。
4、所述的根据实体对象的模型参数构建实体对象内部程式化模型包括：
导入实体对象的模型参数，如果该实体对象为“SL”，内部为空的壳状结构，则根据该实体对象的厚度与外部程式化模型构建实体对象的内部程式化模型。
5、所述的使用组合几何方法构建实体对象的程式化模型并写为标准输入文件包括：
使用组合几何方法，利用布尔运算将所有实体对象外部程式化模型、内部程式化模型组合成为实体对象的程式化模型，并将实体对象的程式化模型与对应的模型参数写为标准输入文件；
组合几何方法的过程如下：
(1)为外部程式化模型、内部程式化模型编号；
(2)在内部程式化模型编号前附上“+”，外部程式化模型编号前附上“-”，其中“-”表示在基本体内部，“+”表示在基本体外部；
(3)将同一实体对象的外部程式化模型、内部程式化模型用布尔逻辑符“AND”连接构成实体对象的程式化模型，将不同实体对象的程式化模型用布尔逻辑符“OR”连接构成场景的程式化模型，其中“AND”表示交集，“OR”表示并集。
本发明的优势在于：
1、本发明实现了高效的将复杂核设施3dsMax模型简化为可计算辐射剂量分布的程式化模型；
2、本发明实现了对内部为空的壳状核设施模型的辐射场剂量分布计算。
附图说明
图1为本发明的流程图；
图2a为根据实体对象包围盒构建程式化模型(柱体)，图2b为根据实体对象包围盒构建程式化模型(椭球体)，图2c为根据实体对象包围盒构建程式化模型(立方体)；
图3为使用组合几何方法构建实体对象的程式化模型。
具体实施方式
下面结合附图举例对本发明做更详细地描述：
结合图1-3，本发明采用3dsMax软件建模与C++语言编程，主要功能为：构建确定参数的核设施模型，使用程式化建模和组合几何方法，将核设施3dsMax模型简化为程式化模型，该程式化模型生成的标准输入文件可用点核积分方法和MCNP计算辐射场剂量分布。整个发明包括3dsMax建模、核设施模型程式化、生成标准输入文件三个模块。
本发明采用以下技术方案：
1、用3dsMax软件根据确定的核设施参数构建模型，并将文件保存为3DS格式。
1)在3dsMax软件中直接为模型对象命名。
源项类模型命名以SO为区分符开头，命名格式为：
SO,名称
举例：SO,Box1
内部非空的实体类屏蔽模型命名以SH为区分符开头，命名格式为：
SH,名称
举例：SH,Box2
内部为空的壳类屏蔽模型命名以SL为区分符开头，命名格式为：
SL,壳最厚处的厚度,名称
举例：SL,4,Sphere1
2)在3dsMax材质编辑器内对材质球命名，将材质球导入到使用该材质的实体上。
源项类命名格式为：
能量,活度,±份额,其他信息
举例：0.052,3.27E+06,-1,U
屏蔽类材质命名格式为：
密度,原子序数,±份额,其他信息
举例：2.702,13,-1.0,Al
其中“+”代表原子个数的份额，“-”代表质量份额。
2、导入3DS核设施模型文件，获得模型参数。
导入3DS文件，读取3DS文件顶点列表，具体的算法如下：
1)定义3个动态指针，分别指向3个动态数组，用来存储模型的所有顶点，定义一个变量count，用来保存模型的顶点总数；
2)将所有的顶点读入到指定的内存；
3)删除指针变量，释放内存。
3、构建模型包围盒。
遍历模型所有顶点，根据模型顶点在x、y、z方向上的最大最小值构建模型场景立方体包围盒。遍历每个实体对象的所有顶点，构建实体对象立方体包围盒，尺寸为△x×△y×△z，包围盒中心点为O(xo,yo,zo)。
4、根据实体对象包围盒构建程式化模型。
如图2所示，根据实体对象包围盒的中心点与尺寸，构建基于包围盒的椭球体、柱体和立方体的程式化模型。
1)椭球体方程为：
Figure BDA0001605845620000051
当△x＝△y＝△z时，方程为球体方程。
2)柱体的母线平行于坐标轴。下面给出了母线平行z轴的柱体方程：
Figure BDA0001605845620000052
当△x＝△y时，方程为圆柱体方程。母线平行x轴、y轴的柱体方程与此类似。
3)立方体方程为：
Figure BDA0001605845620000053
5、根据实体对象顶点在程式化模型内的个数确定实体对象外部程式化模型。
实验证明，如果实体对象85％以上的顶点都在程式化模型内，则该程式化模型能较好的拟合实体对象。计算椭球体和柱体4种程式化模型包含顶点的个数，选择满足顶点数目要求且顶点最多的模型作为实体对象的外部程式化模型。如果顶点个数相同，则选择椭球体作为实体对象的外部程式化模型。如果椭球体和柱体程式化模型都不满足数目要求，则选择立方体作为实体对象的外部程式化模型。
6、根据实体对象的模型参数构建实体对象内部程式化模型。
导入实体对象的模型参数，如果该实体对象为“SL”，是内部为空的壳状结构，则根据该实体对象的厚度与外部程式化模型构建实体对象的内部程式化模型。
7、使用组合几何方法构建实体对象的程式化模型并写为标准输入文件。
使用组合几何方法，利用布尔运算将所有实体对象外部程式化模型、内部程式化模型组合成为实体对象的程式化模型，并将实体对象的程式化模型与对应的模型参数写为标准输入文件。如图3所示，组合几何方法的主要过程如下：
1)为外部程式化模型、内部程式化模型编号。
2)在内部程式化模型编号前附上“+”，外部程式化模型编号前附上“-”。其中“-”表示在基本体内部，“+”表示在基本体外部。
3)将同一实体对象的外部程式化模型、内部程式化模型用布尔逻辑符“AND”连接构成实体对象的程式化模型，将不同实体对象的程式化模型用布尔逻辑符“OR”连接构成场景的程式化模型。其中“AND”表示交集(通常省略AND)，“OR”表示并集。
8、将标准输入文件输入到辐射剂量计算程序内获得辐射场剂量分布。

Claims (5)

1.一种核设施退役模型程式化仿真方法，其特征是：
(1)用3dsMax软件根据确定的核设施参数构建模型文件，并将模型文件保存为3DS格式；
(2)导入3DS格式的模型文件，获得模型参数；
(3)构建模型包围盒；
(4)根据实体对象包围盒构建程式化模型；
(5)根据实体对象顶点在程式化模型内的个数确定实体对象外部程式化模型；
(6)根据实体对象的模型参数构建实体对象内部程式化模型；
(7)使用组合几何方法构建实体对象的程式化模型并写为标准输入文件；
(8)将标准输入文件输入到辐射剂量计算程序内获得辐射场剂量分布；
所述的根据实体对象包围盒构建程式化模型包括，
根据实体对象包围盒的中心点与尺寸，构建基于包围盒的椭球体、柱体和立方体的程式化模型：
a、椭球体方程为：
Figure FDA0003416296070000011
当Δx＝Δy＝Δz时，方程为球体方程；
b、柱体的母线平行于坐标轴：
母线平行z轴的柱体方程：
Figure FDA0003416296070000012
当Δx＝Δy时，方程为圆柱体方程；
c、立方体方程为：
Figure FDA0003416296070000021
2.根据权利要求1所述的一种核设施退役模型程式化仿真方法，其特征在于：所述的构建模型包围盒包括：
遍历模型所有顶点，根据模型顶点在x、y、z方向上的最大最小值构建模型场景立方体包围盒，遍历每个实体对象的所有顶点，构建实体对象立方体包围盒，尺寸为Δx×Δy×Δz，包围盒中心点为O(xo,yo,zo)。
3.根据权利要求1所述的一种核设施退役模型程式化仿真方法，其特征在于：所述的根据实体对象顶点在程式化模型内的个数确定实体对象外部程式化模型包括：
计算椭球体和柱体程式化模型包含顶点的个数，选择满足顶点数目要求且顶点最多的模型作为实体对象的外部程式化模型，如果顶点个数相同，则选择椭球体作为实体对象的外部程式化模型，如果椭球体和柱体程式化模型都不满足数目要求，则选择立方体作为实体对象的外部程式化模型。
4.根据权利要求1所述的一种核设施退役模型程式化仿真方法，其特征在于：所述的根据实体对象的模型参数构建实体对象内部程式化模型包括：
导入实体对象的模型参数，如果该实体对象为“SL”，内部为空的壳状结构，SL为内部为空的壳类屏蔽模型命名的区分符开头，则根据该实体对象的厚度与外部程式化模型构建实体对象的内部程式化模型。
5.根据权利要求1所述的一种核设施退役模型程式化仿真方法，其特征在于：所述的使用组合几何方法构建实体对象的程式化模型并写为标准输入文件包括：
使用组合几何方法，利用布尔运算将所有实体对象外部程式化模型、内部程式化模型组合成为实体对象的程式化模型，并将实体对象的程式化模型与对应的模型参数写为标准输入文件；
组合几何方法的过程如下：
(1)为外部程式化模型、内部程式化模型编号；
(2)在内部程式化模型编号前附上“+”，外部程式化模型编号前附上“-”，其中“-”表示在基本体内部，“+”表示在基本体外部；
(3)将同一实体对象的外部程式化模型、内部程式化模型用布尔逻辑符“AND”连接构成实体对象的程式化模型，将不同实体对象的程式化模型用布尔逻辑符“OR”连接构成场景的程式化模型，其中“AND”表示交集，“OR”表示并集。
CN201810243219.2A 2018-03-23 2018-03-23 一种核设施退役模型程式化仿真方法 Active CN108563839B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810243219.2A CN108563839B (zh) 2018-03-23 2018-03-23 一种核设施退役模型程式化仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810243219.2A CN108563839B (zh) 2018-03-23 2018-03-23 一种核设施退役模型程式化仿真方法

Publications (2)

Publication Number Publication Date
CN108563839A CN108563839A (zh) 2018-09-21
CN108563839B true CN108563839B (zh) 2022-04-05

Family

ID=63532880

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810243219.2A Active CN108563839B (zh) 2018-03-23 2018-03-23 一种核设施退役模型程式化仿真方法

Country Status (1)

Country Link
CN (1) CN108563839B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110489789B (zh) * 2019-07-10 2023-04-18 哈尔滨工程大学 一种核设施退役环境中放射性气体扩散评估方法
CN111161416B (zh) * 2019-12-11 2023-08-29 北京互时科技股份有限公司 根据模型形状信息精确调整模型显示优先级的方法和系统
CN112288874B (zh) * 2020-10-20 2023-03-24 哈尔滨工程大学 基于cad模型与布尔运算的伽马辐射建模计算仿真方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106297935A (zh) * 2016-11-01 2017-01-04 中国核动力研究设计院 一种废物坑内低放废物的回取方法
CN106934082A (zh) * 2015-12-30 2017-07-07 核工业北京地质研究院 基于德尔菲法的高放废物地质处置库场址适宜性评价方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2011326330B2 (en) * 2010-11-08 2014-12-18 Ezymine Pty Limited A collision avoidance system and method for human commanded systems
JP2014137061A (ja) * 2013-01-18 2014-07-28 Kazuo Ogami 火力発電
CN104377693A (zh) * 2014-11-11 2015-02-25 国家电网公司 一种发电生产模拟模型
CN105956781A (zh) * 2016-05-12 2016-09-21 国网甘肃省电力公司电力科学研究院 电网设备全寿命周期费用计算方法
CN107026450B (zh) * 2016-09-19 2019-10-15 国网重庆市电力公司电力科学研究院 基于遗传算法的b2g和v2g配电网功率控制系统及方法
CN107194103B (zh) * 2017-06-07 2020-11-10 哈尔滨工程大学 一种基于3dsMax的核设施模型辐射场剂量仿真方法
CN107330187B (zh) * 2017-06-30 2020-11-03 哈尔滨工程大学 一种核设施退役辐射场剂量分布仿真方法
CN107330186B (zh) * 2017-06-30 2020-06-16 哈尔滨工程大学 一种3dsMax核设施模型辐射场剂量分布仿真方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934082A (zh) * 2015-12-30 2017-07-07 核工业北京地质研究院 基于德尔菲法的高放废物地质处置库场址适宜性评价方法
CN106297935A (zh) * 2016-11-01 2017-01-04 中国核动力研究设计院 一种废物坑内低放废物的回取方法

Also Published As

Publication number Publication date
CN108563839A (zh) 2018-09-21

Similar Documents

Publication Publication Date Title
CN108563839B (zh) 一种核设施退役模型程式化仿真方法
CN107194103B (zh) 一种基于3dsMax的核设施模型辐射场剂量仿真方法
CN107330186B (zh) 一种3dsMax核设施模型辐射场剂量分布仿真方法
CN109190144A (zh) 一种任意形状放射源辐射屏蔽计算仿真方法
WO2016119537A1 (zh) 基于深度剥离的核辐射屏蔽计算网格生成方法和系统
CN109344459B (zh) 一种蒙特卡罗程序的可视建模与转换方法及系统
CN108427861A (zh) 一种基于mpt工具包构建材料周期性多晶结构的方法
CN106126928B (zh) 获得固态及液态氟盐热中子散射数据库的方法及数据库
CN108804865A (zh) 一种用于模拟等离子体中氘氚聚变反应的方法、系统、计算机设备及存储介质
CN113672849A (zh) 一种基于轴向通量展开的一步法输运计算方法及系统
CN113609099A (zh) 一种基于蒙特卡罗方法制作聚变堆多群屏蔽数据库的方法
Martz Mcnp6 unstructured mesh initial validation and performance results
Aarnio et al. A long writeup of the FLUKA82 program
CN108470079A (zh) 航天发射场涉核操作辐射安全评估仿真方法
Lyytinen Parametric studies of breeding blanket thickness of a HELIAS stellarator using Monte Carlo neutron transport code Serpent2
Wang et al. Development of a hybrid parallelism Monte Carlo transport middleware on mesh geometry
US20200027003A1 (en) Package method of radioactive dismantled parts
Hart Automated Doppler broadening of cross sections for neutron transport applications
Chapoutier et al. AREVA developments for an efficient and reliable use of Monte Carlo codes for radiation transport applications
Ponomarev et al. NASA-developed ProE-based tool for the ray-tracing of spacecraft geometry to determine radiation doses and particle fluxes in habitable areas of spacecraft and in the human body
Shriwise Geometry Query Optimizations in CAD-based Tessellations for Monte Carlo Radiation Transport
CN113901029B (zh) 一种支持堆芯多专业计算软件的统一建模方法及系统
US20220284142A1 (en) Systems and methods for computer-aided design (cad) exchange
Baxter Nonlinear Solution of the Time Eigenvalue of a Fast Burst Reactor Using the Finite Volume Method
Phillips et al. Solving the Discretised Boltzmann Transport Equations using Neural Networks: Applications in Neutron Transport

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