CN111460713B - 基于包壳材料在电磁感应加热条件下的温度分布有限元分析法 - Google Patents
基于包壳材料在电磁感应加热条件下的温度分布有限元分析法 Download PDFInfo
- Publication number
- CN111460713B CN111460713B CN202010244135.8A CN202010244135A CN111460713B CN 111460713 B CN111460713 B CN 111460713B CN 202010244135 A CN202010244135 A CN 202010244135A CN 111460713 B CN111460713 B CN 111460713B
- Authority
- CN
- China
- Prior art keywords
- field
- electromagnetic
- heat
- temperature distribution
- vector
- 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
Links
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提供了基于包壳材料在电磁感应加热条件下的温度分布有限元分析法;对LOCA条件下的包壳材料进行热场分析;通过前处理软件SOLIDWORKS,将被研究物体进行3D建模,并将建好的模型导入到多物理场耦合软件COMSOL中;对感应加热线圈进行电磁特性分析,得出电磁场的分布规律云图,并将其转化为热源,耦合到固体传热模块中去,进而得出被研究物体的表面和中心处的温升;本发明目的是研究在不同温度升高情况下,金属体内温度分布情况,并和实验进行对比,最终实现用仿真分析的方法去替代实验测试。
Description
技术领域
本发明涉及到特殊样品的温度测量领域,特别是一种基于有限元分析的包壳材料在电磁感应加热下的温度分析方法。
背景技术
安全壳作为保护反应堆安全的重要部件,研究其在LOCA事故中的温度分布对维护反应堆的安全性至关重要,为了减少对外部环境的放射性释放,必须在事故期间保证安全壳结构的完整性,即在LOCA条件下安全壳所承受的最高压力和温度应均低于设计值。其中,锆合金因其优异的辐照、力学和抗腐蚀性能,已经广泛被应用于安全壳材料当中。然而对锆合金在LOCA条件下服役的相关实验研究受限于实验条件的限制而未能得到充分研究,因此需要运用模拟仿真的方法对合金内部温度场的分布进行计算。
发明内容
本发明是为了解决上述的问题,设计了一种基于有限元分析的在LOCA条件下获取样品温升的方法。
实现上述目的本发明的技术方案为,一种基于包壳材料在电磁感应加热条件下的温度分布有限元分析法。该方法包含如下几个步骤:
步骤一:根据锆合金包壳的具体尺寸,以及周围加热装置的尺寸及线圈的位置排布,在SolidWorks中建立实体3D模型;
步骤二:将3D模型导入到多物理场耦合软件COMSOL中去,在软件中形成几何实体。定义好各种材料的电磁特性及热特征参数。对输入的几何体进行网格剖分,包括选择单元类型,设置网格尺寸,测试网格精度,做好所有的前处理工作;
步骤三:选择磁场和传热场作为计算物理场。选用麦克斯韦方程组分析和研究电磁场的规律
步骤四:为了简化电磁场计算问题,需要定义两个参数(矢量磁势和标量电势)把电场和磁场的变量分离开计算,分别得到一个独立的电场和磁场的偏微分方程;
步骤五:建立材料属性和场量之间的本构方程;
步骤六:定义线圈的属性,包括线圈匝数,横截面积,输入电流大小和方向等;
步骤七:设置电磁场求解边界条件;
步骤八:设置能量守恒方程求解温度场;
其中ρ0,Cp,k分别代表物质的密度,比热容和热导率,热源项也包含了自身的产热,对外的热传导以及辐射散热三个部分;
步骤九:添加电磁热源,表达式如下:
步骤十:设置边界条件,金属的两个侧边采用热通量边界条件,金属上表面设置表面对环境辐射边界条件,其余边界为绝热边界。
步骤十一:求解器设置,磁场采用线圈几何特性求解,电磁场和传热场耦合采用频域-瞬态求解器进行求解。
步骤十二:绘制加热金属的整体温度分布云图,采集探测点的表面的温升数据并绘制线图。
步骤十三:与实验结果进行对比,验证仿真结果的精确性,从而成为代替实验测试的方法。
与现有技术比,本发明具有如下优势:利用本发明的技术方案,基于包壳材料在电磁感应加热条件下的温度分布有限元分析法,能解决现有的在LOCA条件下,传感器探头无法贴近被测物体表面,红外法无法测量被测对象的真实温度的局限。通过有限元仿真计算分析,有效地缩短了表面温升的测量时间,大大地提高了工作效率。
附图说明
本说明书包括以下附图,所示内容分别是:
图1是本发明所述的有限元分析法的流程示意图。
图2是本发明所述的整体三维模型。
图3是本发明所述的电流激励曲线。
图4是本发明所述的样品温度分布云图。
图5是本发明所述的样品表面温升曲线。
图6是本发明所述的仿真温升曲线与实际测量曲线的对比。
图中标识:10-温度检测点;20-线圈;30-被加热金属。
具体实施方式
下面对照附图,通过对实施例的描述,对本发明的具体实施方式作进一步详细的说明,目的是帮助本领域的技术人员对本发明的构思、技术方案有更完整、准确和深入的理解,并有助于其实施。
如图1所示,该方法包含如下几个步骤:
步骤一:根据锆合金包壳的具体尺寸,以及周围加热装置的尺寸及线圈的位置排布,在SolidWorks中建立实体3D模型,3D模型的结构特征为被加热金属设置于线圈中心处;
步骤二:将3D模型导入到多物理场耦合软件COMSOL中去,在软件中形成几何实体,如图2所示。
定义好材料的电磁特性及热特征参数,如下表所示:
SG1 | 4.1*0.01 | 0.041 | 电阻率参数 |
SG2 | 1/0.017241*SG1[MS/m] | 2.3781E6 S/m | 金属电导率 |
A1 | 22.6[W/(m*K)] | 22.6W/(m·K) | 金属导热系数 |
A3 | 6.49[g/cm^3] | 6490kg/m3 | 金属密度 |
k2 | A2/(A3*k1) | 277.13J/(kg·K) | 金属热容(转换后) |
对输入的几何体进行网格剖分,采用网格最大单元大小为0.01m,最小单元大小为1e-4m,完整网格包含487663个域单元、15427个边界元和552个边单元。
步骤三:选择磁场和传热场作为计算物理场。选用麦克斯韦方程组分析和研究电磁场的规律;
步骤四:为了简化电磁场计算问题,需要定义两个参数(矢量磁势和标量电势)把电场和磁场的变量分离开计算,分别得到一个独立的电场和磁场的偏微分方程;
步骤五:建立材料属性和场量之间的本构方程;
步骤六:定义线圈的属性,线圈匝数为3,横截面积为1e-6m2,输入的电流大小随时间的变化如图3所示,则求得的电流密度为:
步骤七:设置电磁场求解边界条件。除了线圈的输入边界条件,其余的边界条件均为磁绝缘:
步骤八:设置能量守恒方程求解温度场;
其中ρ0,Cp,k分别代表物质的密度,比热容和热导率,热源项也包含了自身的产热,对外的热传导以及辐射散热三个部分。
步骤九:添加电磁热源,表达式如下:
步骤十:设置边界条件,金属的两个侧边采用热通量边界条件,金属上表面设置表面对环境辐射边界条件,其余边界为绝热边界。
步骤十一:求解器设置,磁场采用线圈几何特性求解,电磁场和传热场耦合采用频域-瞬态求解器进行求解。
步骤十二:采用金属体初始温度300摄氏度,以300摄氏度每分钟的速度上升到500摄氏度,再以600度每分钟的速度上升到800摄氏度,再以300度每分钟的速度上升到1000摄氏度,再以600度每分钟的速度上升到1200摄氏度。相应云图及温升曲线如图4,图5所示。
步骤十三:与实验结果进行对比,如图6所示,结果显示仿真与实验的误差在3%以内,计算结果非常可靠,具备了周期短,易操作,准确性高的条件,说明该方法完全可以用来替代实验测试。
以上结合附图对本发明进行了示例性描述。显然,本发明具体实现并不受上述方式的限制。只要是采用了本发明的方法构思和技术方案进行的各种非实质性的改进;或未经改进,将本发明的上述构思和技术方案直接应用于其它场合的,均在本发明的保护范围之内。
Claims (1)
1.一种基于包壳材料在电磁感应加热条件下的温度分布有限元分析法,其特征在于,该方法包括如下步骤:
步骤一:根据锆合金包壳的具体尺寸,以及周围加热装置的尺寸及线圈的位置排布,建立实体3D模型;
步骤二:将3D模型导入到多物理场耦合软件COMSOL中去,在软件中形成几何实体;定义好锆合金材料的电磁特性及热特征参数;对输入的几何体进行网格剖分,包括选择单元类型,设置网格尺寸,测试网格精度,做好所有的前处理工作;
步骤三:选择磁场和传热场作为计算物理场,选用麦克斯韦方程组分析和研究电磁场的规律:
ρ——电荷体密度(C/m3);
步骤四:为了简化电磁场计算问题,需要定义两个参数的分别为:矢量磁势和标量电势,把电场和磁场的变量分离开计算,分别得到一个独立的电场和磁场的偏微分方程;
步骤五:建立材料属性和场量之间的本构方程;
步骤六:定义线圈的属性,包括线圈匝数、横截面积、输入电流大小和方向;
步骤七:设置电磁场求解边界条件;
步骤八:设置能量守恒方程求解温度场;
其中ρ0,Cp,k分别代表物质的密度,比热容和热导率,热源项也包含了自身的产热,对外的热传导以及辐射散热三个部分;
步骤九:添加电磁热源,表达式如下:
步骤十:设置边界条件,金属的两个侧边采用热通量边界条件,金属上表面设置表面对环境辐射边界条件,其余边界为绝热边界;
步骤十一:求解器设置,磁场采用线圈几何特性求解,电磁场和传热场耦合采用频域-瞬态求解器进行求解;
步骤十二:绘制加热金属的整体温度分布云图,采集探测点的表面的温升数据并绘制线图;
步骤十三:与实验结果进行对比,验证仿真结果的精确性,从而成为代替实验测试的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010244135.8A CN111460713B (zh) | 2020-03-31 | 2020-03-31 | 基于包壳材料在电磁感应加热条件下的温度分布有限元分析法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010244135.8A CN111460713B (zh) | 2020-03-31 | 2020-03-31 | 基于包壳材料在电磁感应加热条件下的温度分布有限元分析法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111460713A CN111460713A (zh) | 2020-07-28 |
CN111460713B true CN111460713B (zh) | 2022-03-01 |
Family
ID=71685095
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010244135.8A Active CN111460713B (zh) | 2020-03-31 | 2020-03-31 | 基于包壳材料在电磁感应加热条件下的温度分布有限元分析法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111460713B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111931370B (zh) * | 2020-08-06 | 2022-08-02 | 天津大学 | 一种基于comsol的陶瓷绝缘子闪烧方法 |
CN112229146B (zh) * | 2020-10-20 | 2022-05-03 | 西安电子科技大学 | 微波烘干的干燥控制方法、系统、设备、仿真优化及应用 |
CN112432968B (zh) * | 2020-10-21 | 2022-08-30 | 中国核动力研究设计院 | 辐照后反应堆结构材料热导率测试样的制备方法及试样盒 |
CN112784463B (zh) * | 2021-01-29 | 2023-12-19 | 华中科技大学 | 一种井中重力仪系统的内部温度仿真方法 |
CN113094955B (zh) * | 2021-04-12 | 2023-02-24 | 兰州交通大学 | 一种微波加热仿真分析方法 |
CN115906555B (zh) * | 2022-10-12 | 2023-09-19 | 中建二局安装工程有限公司 | 基于comsol的焊接过程中残余应力及其漏磁场的预测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106650093A (zh) * | 2016-12-21 | 2017-05-10 | 浙江大学 | 电磁器件分析中电磁‑热‑应力三场耦合的去耦合计算方法 |
CN107331425A (zh) * | 2017-07-07 | 2017-11-07 | 西安交通大学 | 一种熔融物与液态金属相互作用研究的试验系统及方法 |
CN107767025A (zh) * | 2017-09-20 | 2018-03-06 | 中广核研究院有限公司 | 一种风险指引的核电厂大破口失水事故分析方法 |
CN109243641A (zh) * | 2018-10-18 | 2019-01-18 | 中国核动力研究设计院 | 用于压水堆失水事故的反应堆压力容器实验模拟体 |
CN109440045A (zh) * | 2018-12-29 | 2019-03-08 | 日照钢铁控股集团有限公司 | 一种电磁感应加热线圈超薄防护工艺 |
CN109767853A (zh) * | 2018-12-19 | 2019-05-17 | 西安交通大学 | 模拟核反应堆失水事故后包壳骤冷破裂实验装置及方法 |
CN109781215A (zh) * | 2019-01-17 | 2019-05-21 | 合肥工业大学 | 一种电磁流量传感器多场耦合感应电动势的数值仿真方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100040187A1 (en) * | 2008-08-12 | 2010-02-18 | Searete Llc, A Limited Liability Corporation Of The State Of Delaware | Heat pipe nuclear fission deflagration wave reactor cooling |
-
2020
- 2020-03-31 CN CN202010244135.8A patent/CN111460713B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106650093A (zh) * | 2016-12-21 | 2017-05-10 | 浙江大学 | 电磁器件分析中电磁‑热‑应力三场耦合的去耦合计算方法 |
CN107331425A (zh) * | 2017-07-07 | 2017-11-07 | 西安交通大学 | 一种熔融物与液态金属相互作用研究的试验系统及方法 |
CN107767025A (zh) * | 2017-09-20 | 2018-03-06 | 中广核研究院有限公司 | 一种风险指引的核电厂大破口失水事故分析方法 |
CN109243641A (zh) * | 2018-10-18 | 2019-01-18 | 中国核动力研究设计院 | 用于压水堆失水事故的反应堆压力容器实验模拟体 |
CN109767853A (zh) * | 2018-12-19 | 2019-05-17 | 西安交通大学 | 模拟核反应堆失水事故后包壳骤冷破裂实验装置及方法 |
CN109440045A (zh) * | 2018-12-29 | 2019-03-08 | 日照钢铁控股集团有限公司 | 一种电磁感应加热线圈超薄防护工艺 |
CN109781215A (zh) * | 2019-01-17 | 2019-05-21 | 合肥工业大学 | 一种电磁流量传感器多场耦合感应电动势的数值仿真方法 |
Non-Patent Citations (3)
Title |
---|
LOCA条件下典型包壳材料感应加热温升行为的模拟研究;赵鑫;《工程与试验》;20200315;106-108 * |
Preliminary safety analysis of the PWR with accident-tolerant fuels during severe accident conditions;Xiaoli Wua;《Annals of Nuclear Energy》;20150306;1-13 * |
不同工况下APWR反应堆的数值传热计算;李媛媛;《化工机械》;20120215;69-71,95 * |
Also Published As
Publication number | Publication date |
---|---|
CN111460713A (zh) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111460713B (zh) | 基于包壳材料在电磁感应加热条件下的温度分布有限元分析法 | |
CN104036125B (zh) | 一种油浸式变压器内部温度场的精确计算方法 | |
CN105740536B (zh) | 对直缝焊管中频热处理过程进行动态仿真的模拟方法 | |
Yang et al. | Simulation of edge cracks using pulsed eddy current stimulated thermography | |
Smolka et al. | Experimental validation of the coupled fluid flow, heat transfer and electromagnetic numerical model of the medium-power dry-type electrical transformer | |
Cheng et al. | A combined experimental-numerical method to evaluate powder thermal properties in laser powder bed fusion | |
Zec et al. | Finite element analysis of nondestructive testing eddy current problems with moving parts | |
Chen et al. | Thermal–hydraulic performance of a 5× 5 rod bundle with spacer grid in a nuclear reactor | |
CN106294914A (zh) | 一种考虑非线性各向异性的换流变压器阀侧三维电场仿真方法 | |
CN105183988B (zh) | 一种接地螺栓通电后温度及应力应变有限元计算分析方法 | |
CN103823114A (zh) | 介质材料不同温度下辐射诱导电导率测量装置及方法 | |
CN111597747A (zh) | 一种多点触发三元锂动力电池模组热失控仿真及预测方法 | |
Lin et al. | Distributed current source method for modeling magnetic and eddy-current fields induced in nonferrous metallic objects | |
Feng et al. | Further validation of liquid metal MHD code for unstructured grid based on OpenFOAM | |
Gajbhiye et al. | Validation and verification of a robust 3-D MHD code | |
Chaharmiri et al. | The effect of stepped field shaper on magnetic pressure and radial displacement in electromagnetic inside bead forming: experimental and simulation analyses using MAXWELL and ABAQUS software | |
Nili et al. | Coupled electro-thermo-mechanical simulation for multiple pellet fabrication using spark plasma sintering | |
Malekzadeh et al. | Magnetic field effect on laminar heat transfer in a pipe for thermal entry region | |
CN105550463B (zh) | 钢板电磁感应加热过程温度场的预测方法 | |
Ocilka et al. | Simulation model of induction heating in Comsol multiphysics | |
Belik et al. | Induction Heating of a Fluidized Pebble Bed: Numerical and Experimental Analysis | |
Zhu et al. | A novel nonlinearity marginalization technique for effective solution of induction heating problems by cell method | |
Vyas et al. | Entropy analysis of radiative MHD forced convection flow with weakly temperature dependent convection coefficient in porous medium channel | |
Liu et al. | Coupled modeling of electromagnetic-thermal problem in induction heating process considering material properties | |
Landi et al. | A virtual design approach to simulate the hob energy performance |
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 |