CN112446147A - 一种提高压紧力裕量的燃料组件压紧系统分析方法 - Google Patents

一种提高压紧力裕量的燃料组件压紧系统分析方法 Download PDF

Info

Publication number
CN112446147A
CN112446147A CN202011314003.4A CN202011314003A CN112446147A CN 112446147 A CN112446147 A CN 112446147A CN 202011314003 A CN202011314003 A CN 202011314003A CN 112446147 A CN112446147 A CN 112446147A
Authority
CN
China
Prior art keywords
fuel assembly
pressing force
force
hold
parameters
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
Application number
CN202011314003.4A
Other languages
English (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.)
Nuclear Power Institute of China
Original Assignee
Nuclear Power Institute of China
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 Nuclear Power Institute of China filed Critical Nuclear Power Institute of China
Priority to CN202011314003.4A priority Critical patent/CN112446147A/zh
Publication of CN112446147A publication Critical patent/CN112446147A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/70Machine learning, data mining or chemometrics
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Software Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明涉及核反应堆燃料组件技术领域,具体公开了一种提高压紧力裕量的燃料组件压紧系统分析方法。该方法包括:根据燃料组件在堆芯中的受力分析,建立燃料组件压紧力裕量模型;根据压紧力影响参数类型,确定压紧弹簧变形量,获得各种工况下的燃料组件压紧力;获取燃料组件压紧力裕量;在规定的次数内,重复计算获得燃料组件压紧力和压紧力裕量,并在给定置信度和容许区间情况下,对各种工况下的压紧力和压紧力裕量的计算结果进行数据统计。本方法将压紧力影响参数按照一定的概率分布进行取值,相对于确定论分析方法,不用各变量的极限值,可解决影响参数取极限值带来的保守性问题。

Description

一种提高压紧力裕量的燃料组件压紧系统分析方法
技术领域
本发明属于核反应堆燃料组件技术领域,具体涉及一种提高压紧力裕量的燃料组件压紧系统分析方法。
背景技术
燃料组件压紧系统是核反应堆燃料组件的重要部件,其主要功能是补偿燃料组件的高度差以及燃料组件辐照生长和松弛引起的高度变化,为燃料组件提供合适的轴向压紧力,以防止燃料组件跳起,从而保证燃料组件的结构完整性。压紧力过小会使燃料组件脱离堆芯下板的支撑,压紧力过大又易使燃料组件在运行期间发生变形,可能导致控制棒组件落棒不畅。在国家核安全局核电厂安全审查中,燃料组件压紧系统性能是关注的重点之一。
目前,一般采用确定论方法开展燃料组件压紧系统分析计算。采用确定论方法计算时,对各种影响压紧力的参数均取极限值,并且对不确定性进行算术累加。事实上,各种影响压紧力的参数均同时处于极限值的概率是非常小的。因此,采用确定论方法验证计算得到的最小压紧力过于保守,往往使得压紧力裕量不足,计算结果难以满足燃料组件设计准则要求。为更准确地评估燃料组件压紧力,保证合适压紧力的同时又减小压紧力对燃料组件弯曲变形的影响,需要采用一种更合理的分析方法开展燃料组件压紧系统分析。
发明内容
本发明的目的在于提供一种提高压紧力裕量的燃料组件压紧系统分析方法,克服现有技术中确定论分析方法的缺点,解决燃料组件压紧系统的压紧力计算过于保守导致压紧力裕量不足的问题。
本发明的技术方案如下:一种提高压紧力裕量的燃料组件压紧系统分析方法,该方法包括:
S1、根据燃料组件在堆芯中的受力分析,建立燃料组件压紧力裕量模型;
S2、根据压紧力影响参数类型,确定压紧弹簧变形量,获得各种工况下的燃料组件压紧力;
S3、对燃料组件的重力、水力学作用力以及浮力随机取值,并根据燃料组件压紧力裕量模型,获取燃料组件压紧力裕量;
S4、在规定的次数内,重复计算获得燃料组件压紧力和压紧力裕量,并在给定置信度和容许区间情况下,对各种工况下的压紧力和压紧力裕量的计算结果进行数据统计。
所述S1中燃料组件压紧力裕量模型具体为:
燃料组件压紧力裕量满足以下要求:
M=Fr+W-Fh-A≥0
其中,M为给定置信度下的燃料组件压紧力裕量统计值,Fr为给定置信度下的压紧力统计值,W为燃料组件重力,Fh为燃料组件水力学作用力,A为燃料组件浮力。
所述S2中获得各种工况下燃料组件压紧力具体包括:
S2.1、确定燃料组件压紧力影响参数类型;
S2.1.1、梳理确定燃料组件压紧力影响参数;
S2.1.2、对燃料组件压紧力影响参数属于确定值参数或随机性参数进行分类判断;
S2.2、确定燃料组件压紧力影响中的随机性参数的不确定性;
S2.3、根据压紧力影响参数,计算获得压紧弹簧变形量及压紧力。
所述S2.3中获得压紧弹簧变形量及压紧力的具体步骤为:
S2.3.1、获取与燃料组件弹簧变形量相关参数的随机取值,进行弹簧变形量计算;
S2.3.2、对燃料组件辐照生长及弹簧辐照松弛进行随机取值,并获取各种工况下的压紧系统压紧力。
所述S2.3.1中与燃料组件弹簧变形量相关的参数具体包括:堆腔高度、管座高度、弹簧片高度及厚度;获取每个工况下的所述的堆腔高度、管座高度、弹簧片高度及厚度参数,对上述参数进行随机取值,根据各参数对随机取值,计算获得弹簧变形量。
所述S2.1.2中燃料组件压紧力影响参数属于确定值参数或随机性参数具体为;
将影响燃料组件压紧力中的对应工况的温度参数、材料的热膨胀系数、板弹簧特性曲线划分为确定值参数;将堆腔高度参数、燃料组件结构参数、燃料组件重量及水力载荷数据则归为随机性参数。
所述S3中获取燃料组件压紧力裕量具体包括:
对燃料组件重力W、燃料组件水力学作用力Fh、为燃料组件浮力A进行随机取值,并根据S2中获得压紧力,通过燃料组件压紧力裕量计算公式,获得燃料组件压紧力裕量。
所述S4具体包括:
在设定的计算次数内,获得给定置信度和容许区间下的压紧力统计值,在给定置信度95%和容许区间95%情况下,对各种工况下的压紧力和压紧力裕量的计算结果进行数据统计。
所述S4中对压紧力和压紧力裕量的计算结果进行数据统计具体公式为:
计算结果的算术平均值:
Figure BDA0002790736210000041
计算结果的标准差:
Figure BDA0002790736210000042
计算输出的最小值:
Figure BDA0002790736210000043
计算输出的最大值:
Figure BDA0002790736210000044
式中,n为设定的计算次数,data为用于统计的计算结果数据,k为给定的置信度和容许区间下双侧统计容许区间系数。
所述S2.2中确定燃料组件压紧力影响中的随机性参数的不确定性是根据燃料组件压紧力影响中的随机性参数的名义值、公差以及随机数分布类型,采用现有随机数生成算法,获得各随机性参数的随机数。
本发明的显著效果在于:本发明所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,将压紧力影响参数按照一定的概率分布进行取值,通过给定的置信度水平进行结果统计,计算结果更接近实际值,因而可获得较高的压紧力裕量。模拟计算次数越多,计算得到的结果越精确。相对于确定论分析方法,本发明不用各变量的极限值,可解决影响参数取极限值带来的保守性问题。
附图说明
图1为本发明所述的一种提高压紧力裕量的燃料组件压紧系统分析方法流程图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
实施例1
如图1所示,一种提高压紧力裕量的燃料组件压紧系统分析方法,该方法具体包括:
S1、根据燃料组件在堆芯中的受力分析,建立燃料组件压紧力裕量模型;
根据燃料组件在堆芯中的受力情况,建立燃料组件压紧力裕量满足以下要求:
M=Fr+W-Fh-A≥0
其中,M为给定置信度下的燃料组件压紧力裕量统计值,Fr为给定置信度下的压紧力统计值,W为燃料组件重力,Fh为燃料组件水力学作用力,A为燃料组件浮力;
S2、根据压紧力影响参数类型,确定压紧弹簧变形量,获得各种工况下的压紧系统的压紧力;
S2.1、确定燃料组件压紧力影响参数类型;
S2.1.1、梳理确定燃料组件压紧力影响参数;
S2.1.2、对燃料组件压紧力影响参数属于确定值参数或随机性参数进行分类判断;
将影响燃料组件压紧力中的对应工况的温度参数、材料的热膨胀系数、板弹簧特性曲线划分为确定值参数;将堆腔高度参数、燃料组件结构参数、燃料组件重量及水力载荷数据则归为随机性参数;
S2.2、确定燃料组件压紧力影响中的随机性参数的不确定性;
根据燃料组件压紧力影响中的随机性参数的名义值、公差以及随机数分布类型,采用现有随机数生成算法,获得各随机性参数的随机数;
各随机性参数的随机数分布类型根据实际情况选取,例如可选取正态分布或均匀分布,其中,当一个输入参数不能确定为正态分布,从保守上考虑则可将该参数确定为均匀分布;
S2.3、根据压紧力影响参数,计算获得压紧弹簧变形量及压紧力;
S2.3.1、获取与燃料组件弹簧变形量相关参数的随机取值,进行弹簧变形量计算;
获取每个工况下的堆腔高度、管座高度、弹簧片高度及厚度参数,对上述参数进行随机取值,根据各参数对随机取值,计算获得弹簧变形量;
S2.3.2、对燃料组件辐照生长及弹簧辐照松弛进行随机取值,并获取各种工况下的压紧系统压紧力;
S3、对燃料组件的重力、水力学作用力以及浮力随机取值,并根据燃料组件压紧力裕量模型,获取燃料组件压紧力裕量;
对燃料组件重力W、燃料组件水力学作用力Fh、为燃料组件浮力A进行随机取值,并根据S2中获得压紧力,通过燃料组件压紧力裕量计算公式,获得燃料组件压紧力裕量;
S4、在规定的次数内,重复计算获得燃料组件压紧力和压紧力裕量,并在给定置信度和容许区间情况下,对各种工况下的压紧力和压紧力裕量的计算结果进行数据统计;
在设定的计算次数内,获得给定置信度和容许区间下的压紧力统计值,在给定置信度(一般为95%)和容许区间(容许区间选定的总体比例一般为95%)情况下,对各种工况下的压紧力和压紧力裕量的计算结果进行数据统计;
各种工况下的压紧力和压紧力裕量的数据统计计算方法如下:
计算结果的算术平均值:
Figure BDA0002790736210000061
计算结果的标准差:
Figure BDA0002790736210000062
计算输出的最小值:
Figure BDA0002790736210000063
计算输出的最大值:
Figure BDA0002790736210000064
式中,n为设定的计算次数,data为用于统计的计算结果数据,k为给定的置信度和容许区间下双侧统计容许区间系数,其取值可参考GB3359-1982。
实施例2
本实施例与实施例1不同之处在于,压紧力裕量M值也为给定置信度下的统计值。通过对燃料组件重力W、燃料组件水力学作用力Fh以及燃料组件浮力A进行随机取值,然后参与给定置信度下的压紧力计算,当达到给定的计算次数时,就可以获得给定置信度和容许区间下的压紧力裕量统计值;该压紧力裕量统计值应大于实施例1中的压紧力裕量M值。

Claims (10)

1.一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,该方法包括:
S1、根据燃料组件在堆芯中的受力分析,建立燃料组件压紧力裕量模型;
S2、根据压紧力影响参数类型,确定压紧弹簧变形量,获得各种工况下的燃料组件压紧力;
S3、对燃料组件的重力、水力学作用力以及浮力随机取值,并根据燃料组件压紧力裕量模型,获取燃料组件压紧力裕量;
S4、在规定的次数内,重复计算获得燃料组件压紧力和压紧力裕量,并在给定置信度和容许区间情况下,对各种工况下的压紧力和压紧力裕量的计算结果进行数据统计。
2.根据权利要求1所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,所述S1中燃料组件压紧力裕量模型具体为:
燃料组件压紧力裕量满足以下要求:
M=Fr+W-Fh-A≥0
其中,M为给定置信度下的燃料组件压紧力裕量统计值,Fr为给定置信度下的压紧力统计值,W为燃料组件重力,Fh为燃料组件水力学作用力,A为燃料组件浮力。
3.根据权利要求1所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,所述S2中获得各种工况下燃料组件压紧力具体包括:
S2.1、确定燃料组件压紧力影响参数类型;
S2.1.1、梳理确定燃料组件压紧力影响参数;
S2.1.2、对燃料组件压紧力影响参数属于确定值参数或随机性参数进行分类判断;
S2.2、确定燃料组件压紧力影响中的随机性参数的不确定性;
S2.3、根据压紧力影响参数,计算获得压紧弹簧变形量及压紧力。
4.根据权利要求3所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,所述S2.3中获得压紧弹簧变形量及压紧力的具体步骤为:
S2.3.1、获取与燃料组件弹簧变形量相关参数的随机取值,进行弹簧变形量计算;
S2.3.2、对燃料组件辐照生长及弹簧辐照松弛进行随机取值,并获取各种工况下的压紧系统压紧力。
5.根据权利要求4所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,所述S2.3.1中与燃料组件弹簧变形量相关的参数具体包括:堆腔高度、管座高度、弹簧片高度及厚度;获取每个工况下的所述的堆腔高度、管座高度、弹簧片高度及厚度参数,对上述参数进行随机取值,根据各参数对随机取值,计算获得弹簧变形量。
6.根据权利要求3所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,所述S2.1.2中燃料组件压紧力影响参数属于确定值参数或随机性参数具体为;
将影响燃料组件压紧力中的对应工况的温度参数、材料的热膨胀系数、板弹簧特性曲线划分为确定值参数;将堆腔高度参数、燃料组件结构参数、燃料组件重量及水力载荷数据则归为随机性参数。
7.根据权利要求1所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,所述S3中获取燃料组件压紧力裕量具体包括:
对燃料组件重力W、燃料组件水力学作用力Fh、为燃料组件浮力A进行随机取值,并根据S2中获得压紧力,通过燃料组件压紧力裕量计算公式,获得燃料组件压紧力裕量。
8.根据权利要求1所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,所述S4具体包括:
在设定的计算次数内,获得给定置信度和容许区间下的压紧力统计值,在给定置信度95%和容许区间95%情况下,对各种工况下的压紧力和压紧力裕量的计算结果进行数据统计。
9.根据权利要求1或8所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,所述S4中对压紧力和压紧力裕量的计算结果进行数据统计具体公式为:
计算结果的算术平均值:
Figure FDA0002790736200000031
计算结果的标准差:
Figure FDA0002790736200000032
计算输出的最小值:
Figure FDA0002790736200000033
计算输出的最大值:
Figure FDA0002790736200000034
式中,n为设定的计算次数,data为用于统计的计算结果数据,k为给定的置信度和容许区间下双侧统计容许区间系数。
10.根据权利要求3所述的一种提高压紧力裕量的燃料组件压紧系统分析方法,其特征在于,所述S2.2中确定燃料组件压紧力影响中的随机性参数的不确定性是根据燃料组件压紧力影响中的随机性参数的名义值、公差以及随机数分布类型,采用现有随机数生成算法,获得各随机性参数的随机数。
CN202011314003.4A 2020-11-20 2020-11-20 一种提高压紧力裕量的燃料组件压紧系统分析方法 Pending CN112446147A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011314003.4A CN112446147A (zh) 2020-11-20 2020-11-20 一种提高压紧力裕量的燃料组件压紧系统分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011314003.4A CN112446147A (zh) 2020-11-20 2020-11-20 一种提高压紧力裕量的燃料组件压紧系统分析方法

Publications (1)

Publication Number Publication Date
CN112446147A true CN112446147A (zh) 2021-03-05

Family

ID=74737205

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011314003.4A Pending CN112446147A (zh) 2020-11-20 2020-11-20 一种提高压紧力裕量的燃料组件压紧系统分析方法

Country Status (1)

Country Link
CN (1) CN112446147A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106384610A (zh) * 2016-11-18 2017-02-08 中国核动力研究设计院 一种燃料组件压紧系统载荷评价方法
CN110402467A (zh) * 2017-02-27 2019-11-01 泰拉能源公司 用于为核反应堆建模的系统和方法
US20200176139A1 (en) * 2018-11-29 2020-06-04 Elsoltec Fluid sloshing behavior analysis apparatus and method thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106384610A (zh) * 2016-11-18 2017-02-08 中国核动力研究设计院 一种燃料组件压紧系统载荷评价方法
CN110402467A (zh) * 2017-02-27 2019-11-01 泰拉能源公司 用于为核反应堆建模的系统和方法
US20200176139A1 (en) * 2018-11-29 2020-06-04 Elsoltec Fluid sloshing behavior analysis apparatus and method thereof

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
SY JEON等: "A Hold-down Margin Assessment using Statistical Method for the PWR Fuel Assembly", 《TRANSACTIONS OF THE KOREAN NUCLEAR SOCIETY AUTUMN MEETING》 *
张林等: "燃料组件压紧部件分析研究", 《核动力工程》 *
曾孝敏等: "基于蒙特卡洛方法的燃料组件压紧力研究", 《中国核电》 *
朱发文等: "基于蒙特卡洛方法的燃料组件压紧力计算模型研究", 《核动力工程》 *

Similar Documents

Publication Publication Date Title
CN113076684B (zh) 一种核反应堆堆芯调棒过程瞬态参数智能计算方法
EP2077561A1 (en) Thermal limit analysis with hot-channel model for boiling water reactors
CN105260941A (zh) 一种供应侧参与新能源调峰的技术经济评价方法
CN107229771B (zh) 进行核燃料板弹簧压紧力模拟测定的方法
US6377908B1 (en) Method for optimizing transfer function outputs
Franzén Evaluation of Fuel Assembly Bow Penalty Peaking Factors for Ringhals 3: Based on a Cycle Specific Core Water Gap Distribution
CN113919204B (zh) 一种多状态制造系统可用性综合重要度分析方法
CN112446147A (zh) 一种提高压紧力裕量的燃料组件压紧系统分析方法
CN111274624B (zh) 一种基于rbf代理模型的多工况异形节点拓扑优化设计方法
CN112651153A (zh) 一种确定晶体塑性有限元模型材料参数的方法
CN107221367A (zh) 一种换棒组件再入堆的安全评价方法
CN112613156A (zh) 一种精细化燃料棒性能分析方法
CN107591847B (zh) 一种利用变参数的方式调节水电机组agc的方法
KR102572085B1 (ko) 확률론적 분석방법을 이용한 천연가스 공급 안정성 평가방법
CN112149242A (zh) 一种考虑应力松弛和辐照影响的堆内构件压紧弹簧疲劳可靠性评估方法
CN112100867B (zh) 电网负荷预测方法
CN107403027A (zh) 一种建立计及检修影响的电气设备故障率修正模型的方法
CN102184316A (zh) 一种自动选定变力弹簧型号及确定该型号信息的方法
JP2021039024A (ja) 炉心計算方法、炉心計算プログラムおよび炉心計算装置
Forget et al. Integral Full Core Multi-Physics PWR Benchmark with Measured Data
CN111308977A (zh) 暂降及短时中断严重度评估与治理的工业过程建模方法
CN107918822B (zh) 一种评估在运核电厂汽轮机出力提升能力的方法
CN112418533B (zh) 一种清洁能源电量分解预测方法
CN116933930B (zh) 一种区间降雨影响下的电站坝前水位预测方法
CN116307406B (zh) 一种基于物联网的能效管理方法及系统

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210305