CN109783869B - 一种预测反应堆压力容器焊缝热老化晶界p偏析的方法 - Google Patents
一种预测反应堆压力容器焊缝热老化晶界p偏析的方法 Download PDFInfo
- Publication number
- CN109783869B CN109783869B CN201811545856.1A CN201811545856A CN109783869B CN 109783869 B CN109783869 B CN 109783869B CN 201811545856 A CN201811545856 A CN 201811545856A CN 109783869 B CN109783869 B CN 109783869B
- Authority
- CN
- China
- Prior art keywords
- grain boundary
- segregation
- pressure vessel
- reactor pressure
- grain
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
Abstract
本发明涉及一种预测反应堆压力容器焊缝热老化晶界P偏析的方法,将时间参数t代入反应堆压力容器焊缝热老化晶界P偏析预测模型,根据结果确定反应堆压力容器焊缝由热老化导致的晶界P浓度的变化情况。本发明能够利用三维原子探针测试方法获得的数据计算预测模型中的各个参数,最终获得用以预测压力容器材料脆化程度的模型。采用该方法可以进一步验证模型在理论上的准确性,加强预测结果与材料自身微观特征之间的关联性。通过将预测得到的P浓度与实验测得的结果进行对比,可以不断的对模型进行修正,使得预测结果更加准确。
Description
技术领域
本发明反应堆压力容器安全评估技术领域,具体涉及一种预测反应堆压力容器焊缝热老化晶界P偏析的方法。
背景技术
反应堆压力容器(RPV)作为压水堆核电站最重要的关键设备,包裹着反应堆活性区和其他必要设备且在反应堆全寿期内都无法更换,其性能的可靠性直接决定了核电站的服役年限,关系到核电站运行的经济性与安全性。由于RPV长期处在高温、高压和较强辐射场条件下运行,其老化脆化行为非常复杂,主要包括热老化脆化和辐照脆化。一旦在服役过程中RPV发生了严重的脆化,将直接威胁到反应堆的运行安全,则必须提前关闭核电站,这无疑对核电厂的经济效益造成严重的影响,而焊缝作为RPV的薄弱环节,更是直接影响到整个RPV的安全可靠性。而P偏析引起的晶界结合力下降是导致Ni-Cr-Mo-V钢焊缝在热老化过程中发生脆化的主要原因,因此,P偏析的定量分析是保证反应堆运行安全的重要途径之一。
目前,对于P偏析的预测可以采用实验的方法测得不同热老化时间样品晶界处的P浓度,将P浓度与服役时间画图,而后利用含有未知参数的模型对P浓度-服役时间图进行拟合,获得各参数方法,从而对模型进行完善,进而利用该模型对P偏析进行预测。但这种做法并不是直接计算参数,实验较为繁琐。
发明内容
针对现有技术中存在的缺陷,本发明的目的是提供一种预测反应堆压力容器焊缝热老化晶界P偏析的方法,该技术方案能够实现对反应堆压力容器焊缝热老化晶界P偏析的预测,有助于根据其服役时间来预测脆化情况,保证运行安全。
本发明的技术方案如下:
一种预测反应堆压力容器焊缝热老化晶界P偏析的方法,将时间参数t代入反应堆压力容器焊缝热老化晶界P偏析预测模型,根据结果确定反应堆压力容器焊缝由热老化导致的晶界P浓度的变化情况;
所述反应堆压力容器焊缝热老化晶界P偏析预测模型由以下步骤得到:
(1)对晶界处进行三维原子探针实验,得到元素i在垂直于晶界方向上的一维浓度分布图;
(2)对一维浓度分布图进行积分后获得元素i的积累成分曲线;
(3)根据所述积累成分曲线数据来确定晶界厚度δ,利用公式计算初始晶界富集系数s0和富集系数s;
(4)将晶界厚度δ、初始晶界富集系数s0和富集系数s代入下列公式以得到反应堆压力容器焊缝热老化晶界P偏析预测模型:
上述公式中,Cb(t)为服役时间t时的晶界P浓度;
R为晶粒半径;
δ为晶界厚度;
D为P扩散自由能。
进一步地,上述的预测反应堆压力容器焊缝热老化晶界P偏析的方法,根据所述积累成分曲线数据来确定晶界厚度δ,利用公式计算初始晶界富集系数s0和富集系数s为:
其中,初始晶界富集系数s0:
富集系数s:
确定晶界厚度δ:取lα与lGB的交点B和取lβ与lGB的交点A之间的部分确定为晶界,则AB部分对应的厚度定义为晶界厚度δ;
上述公式中,和分别为元素在晶界两侧晶粒α和β内的原子浓度,为元素i在晶界处的浓度,其中为积累成分曲线在晶粒α内的线性拟合线lα的斜率,为积累成分曲线在晶粒β内的线性拟合线lβ的斜率,为积累成分曲线在晶界处部分的线性拟合线lGB的斜率。
进一步地,上述的预测反应堆压力容器焊缝热老化晶界P偏析的方法,所述晶界厚度δ、初始晶界富集系数s0和富集系数s的最终取值为分别多次计算后的平均值。
进一步地,上述的预测反应堆压力容器焊缝热老化晶界P偏析的方法,在VVER-1000型反应堆压力容器焊缝热老化晶界P偏析预测时,将时间参数t代入的反应堆压力容器焊缝热老化晶界P偏析预测模型为:
本发明的有益效果如下:
本发明能够利用三维原子探针测试方法获得的数据计算预测模型中的各个参数,最终获得用以预测压力容器材料脆化程度的模型。采用该方法可以进一步验证模型在理论上的准确性,加强预测结果与材料自身微观特征之间的关联性。通过将预测得到的P浓度与实验测得的结果进行对比,可以不断的对模型进行修正,使得预测结果更加准确。
附图说明
图1为本发明中得到反应堆压力容器焊缝热老化晶界P偏析预测模型的流程图。
图2为本发明积累成分曲线示意图。
图3为本发明具体实施例中,VVER-1000型RPV用Ni-Cr-Mo-V钢焊缝材料晶界P元素积累成分曲线。
图4为本发明的模型计算值与三维原子探针测量值比较示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提供了一种预测反应堆压力容器焊缝热老化晶界P偏析的方法,将时间参数t代入反应堆压力容器焊缝热老化晶界P偏析预测模型,根据结果确定反应堆压力容器焊缝由热老化导致的晶界P浓度的变化情况。
其中,所述反应堆压力容器焊缝热老化晶界P偏析预测模型由以下步骤得到:
S1、对晶界处进行三维原子探针实验,得到元素i在垂直于晶界方向上的一维浓度分布图;
S2、对一维浓度分布图进行积分后获得元素i的积累成分曲线(如图2所示);
S3、根据所述积累成分曲线数据来确定晶界厚度δ,利用公式计算初始晶界富集系数s0和富集系数s;
S4、将晶界厚度δ、初始晶界富集系数s0和富集系数s代入下列公式以得到反应堆压力容器焊缝热老化晶界P偏析预测模型:
上述公式中,Cb(t)为服役时间t时的晶界P浓度;
R为晶粒半径;
δ为晶界厚度;
其中α是一个系数,没有具体的物理意义;
D为P扩散自由能。
上述方法中,根据所述积累成分曲线数据来确定晶界厚度δ,利用公式计算初始晶界富集系数s0和富集系数s为:
其中,初始晶界富集系数s0:
富集系数s:
确定晶界厚度δ:取lα与lGB的交点B和取lβ与lGB的交点A之间的部分确定为晶界,则AB部分对应的厚度定义为晶界厚度δ;
上述公式中,和分别为元素在晶界两侧晶粒α和β内的原子浓度,为元素i在晶界处的浓度,其中为积累成分曲线在晶粒α内的线性拟合线lα的斜率,为积累成分曲线在晶粒β内的线性拟合线lβ的斜率,为积累成分曲线在晶界处部分的线性拟合线lGB的斜率。
为了提高准确率,所述晶界厚度δ、初始晶界富集系数s0和富集系数s的最终取值为分别多次计算后的平均值。
实施例
图3所示为作为实施例的某VVER-1000型RPV用Ni-Cr-Mo-V钢焊缝材料晶界P元素积累成分曲线。对晶粒α和β内的积累成分曲线进行线性拟合分别得到lα和lβ,由lα和lβ的斜率得到和分别为0.0028at.%和0.0081at.%。对积累成分曲线中间部分进行线性拟合后得到拟合线lGB,由lGB的斜率得到为0.118at.%。将lα与lGB的交点B和lβ与lGB的交点A之间的部分确定为晶界,根据原始数据获得AB部分对应的晶界厚度δ为3.4nm。则富集系数和晶界厚度如下:
δ=3.4nm
将计算得到的s0、s、δ带入到预测模型中得到:
需要说明的是模型Cb(t)中所使用的晶界厚度δ和晶界富集系数s、s0是利用该计算方法在获得大量数据以后平均得到的。
将得到的预测模型结果与实验数据得到的结果进行比较,比较结果如图4所示,从图4中可以看到预测模型与实验点2符合较好,验证了参数计算方法与模型本身的可靠性,需要说明的是实验点1参与模型计算,不进行比较。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (4)
1.一种预测反应堆压力容器焊缝热老化晶界P偏析的方法,其特征在于,将时间参数t代入反应堆压力容器焊缝热老化晶界P偏析预测模型,根据结果确定反应堆压力容器焊缝由热老化导致的晶界P浓度的变化情况;
所述反应堆压力容器焊缝热老化晶界P偏析预测模型由以下步骤得到:
(1)对晶界处进行三维原子探针实验,得到元素i在垂直于晶界方向上的一维浓度分布图;
(2)对一维浓度分布图进行积分后获得元素i的积累成分曲线;
(3)根据所述积累成分曲线数据来确定晶界厚度δ,利用公式计算初始晶界富集系数s0和富集系数s;
(4)将晶界厚度δ、初始晶界富集系数s0和富集系数s代入下列公式以得到反应堆压力容器焊缝热老化晶界P偏析预测模型:
上述公式中,Cb(t)为服役时间t时的晶界P浓度;
R为晶粒半径;
δ为晶界厚度;
D为P扩散自由能。
3.如权利要求2所述的预测反应堆压力容器焊缝热老化晶界P偏析的方法,其特征在于,所述晶界厚度δ、初始晶界富集系数s0和富集系数s的最终取值为分别多次计算后的平均值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811545856.1A CN109783869B (zh) | 2018-12-17 | 2018-12-17 | 一种预测反应堆压力容器焊缝热老化晶界p偏析的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811545856.1A CN109783869B (zh) | 2018-12-17 | 2018-12-17 | 一种预测反应堆压力容器焊缝热老化晶界p偏析的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109783869A CN109783869A (zh) | 2019-05-21 |
CN109783869B true CN109783869B (zh) | 2020-08-21 |
Family
ID=66497186
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811545856.1A Active CN109783869B (zh) | 2018-12-17 | 2018-12-17 | 一种预测反应堆压力容器焊缝热老化晶界p偏析的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109783869B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101030467A (zh) * | 2005-03-23 | 2007-09-05 | 信越化学工业株式会社 | 功能梯度稀土永磁体 |
CN101166839A (zh) * | 2005-08-30 | 2008-04-23 | 株式会社开滋 | 青铜类低铅合金 |
CN102099869A (zh) * | 2008-08-12 | 2011-06-15 | 三菱重工业株式会社 | 原子能设备的运转方法 |
CN102930167A (zh) * | 2012-11-05 | 2013-02-13 | 中国原子能科学研究院 | 反应堆压力容器辐照脆化预测评估方法 |
CN103572201A (zh) * | 2013-11-18 | 2014-02-12 | 中国原子能科学研究院 | 低活化铁素体-马氏体钢表面粉末包埋渗铝及后处理工艺 |
CN103710626A (zh) * | 2013-12-26 | 2014-04-09 | 南阳汉冶特钢有限公司 | Cr-Mo系ASTM A387 Gr22CL2压力容器钢及生产方法 |
CN104404296A (zh) * | 2014-11-24 | 2015-03-11 | 安徽科技学院 | 一种核反应堆用铅铋合金的制备方法 |
EP3091060A1 (en) * | 2011-09-22 | 2016-11-09 | Saint-Gobain Cristaux et Détecteurs | Scintillation compound including a rare earth element and a process of forming the same |
JP2017055211A (ja) * | 2015-09-08 | 2017-03-16 | キヤノン株式会社 | 撮像制御装置およびその制御方法、情報処理装置及びその制御方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2876299B1 (fr) * | 2004-10-12 | 2007-01-12 | Solvay | Membranes echangeuses d'ions structurees dans l'epaisseur et procede pour la fabrication de ces membranes |
US7892915B1 (en) * | 2006-03-02 | 2011-02-22 | National Semiconductor Corporation | High performance SiGe:C HBT with phosphorous atomic layer doping |
JP2014077642A (ja) * | 2012-10-09 | 2014-05-01 | Jfe Steel Corp | 鋼材のhic感受性の評価方法およびそれを用いた耐hic性に優れたラインパイプ用高強度厚鋼板の製造方法 |
CN106378459B (zh) * | 2016-09-13 | 2019-01-29 | 安泰核原新材料科技有限公司 | 一种高硼不锈钢中子吸收材料及其制备方法 |
-
2018
- 2018-12-17 CN CN201811545856.1A patent/CN109783869B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101030467A (zh) * | 2005-03-23 | 2007-09-05 | 信越化学工业株式会社 | 功能梯度稀土永磁体 |
CN101166839A (zh) * | 2005-08-30 | 2008-04-23 | 株式会社开滋 | 青铜类低铅合金 |
CN102099869A (zh) * | 2008-08-12 | 2011-06-15 | 三菱重工业株式会社 | 原子能设备的运转方法 |
EP3091060A1 (en) * | 2011-09-22 | 2016-11-09 | Saint-Gobain Cristaux et Détecteurs | Scintillation compound including a rare earth element and a process of forming the same |
CN102930167A (zh) * | 2012-11-05 | 2013-02-13 | 中国原子能科学研究院 | 反应堆压力容器辐照脆化预测评估方法 |
CN103572201A (zh) * | 2013-11-18 | 2014-02-12 | 中国原子能科学研究院 | 低活化铁素体-马氏体钢表面粉末包埋渗铝及后处理工艺 |
CN103710626A (zh) * | 2013-12-26 | 2014-04-09 | 南阳汉冶特钢有限公司 | Cr-Mo系ASTM A387 Gr22CL2压力容器钢及生产方法 |
CN104404296A (zh) * | 2014-11-24 | 2015-03-11 | 安徽科技学院 | 一种核反应堆用铅铋合金的制备方法 |
JP2017055211A (ja) * | 2015-09-08 | 2017-03-16 | キヤノン株式会社 | 撮像制御装置およびその制御方法、情報処理装置及びその制御方法 |
Non-Patent Citations (5)
Title |
---|
"Assessment of segregation kinetics in water-moderated reactors pressure vessel steels under long-term operation";E.A. Kuleshova et al;《Journal of Nuclear Materials》;20160509;第110-122页 * |
"Effect of subgrain structure on the kinetics of phosphorus segregation in grain boundaries";M.V. Sorokin et al;《Materials Letters》;20150603;第151-154页 * |
"Thermodynamic and Kinetic Modeling of Grain Boundary Equilibrium Segregation of P in α-Fe";Y. Yang et al;《Calphad》;20171231;第1-15页 * |
"反应堆压力容器用Ni-Cr-Mo-V钢焊缝热老化脆化行为研究";王成龙等;《原子能科学技术》;20180320(第07期);第1243-1249页 * |
"热老化RPV模型钢中溶质析出激励的原子尺度模拟研究";贾丽霞;《原子能科学技术》;20141031;第48卷;第535-539页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109783869A (zh) | 2019-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Holt | In-reactor deformation of cold-worked Zr–2.5 Nb pressure tubes | |
US20140092934A1 (en) | Method and System for Evaluating Creep Damage of High Temperature Component | |
Chopra et al. | Review of the margins for ASME code fatigue design curve-effects of surface roughness and material variability. | |
CN102967512A (zh) | 基于非对称循环应力控制加载的低周疲劳寿命预测方法 | |
EP4012723A1 (en) | Method for monitoring irradiation embrittlement of nuclear power plant reactor pressure vessel | |
KR101447012B1 (ko) | 금속합금 구조재의 응력부식균열에 대한 비파괴적 평가방법 및 그 금속구조재의 잔여수명 평가방법 | |
Rahimi et al. | A new method for predicting susceptibility of austenitic stainless steels to intergranular stress corrosion cracking | |
US5901071A (en) | Method of evaluating corrosion resistance of metal material, method of designing alloy of high corrosion resistance, method of diagnosing corroded state of metal material, and method of operating plant | |
Ganesh Kumar et al. | High temperature tensile properties of 316LN stainless steel investigated using automated ball indentation technique | |
Sahu et al. | The effect of grain boundary structure on sensitization behavior in a nickel-based superalloy | |
McRae et al. | Diffusivity of hydrogen isotopes in the alpha phase of zirconium alloys interpreted with the Einstein flux equation | |
Bosch et al. | Microstructure, mechanical properties and IASCC susceptibility of stainless steel baffle bolts after 30 years of operation in a PWR | |
CN117723446B (zh) | 基于扩散多元节及skpfm测量氢扩散系数的方法 | |
CN109783869B (zh) | 一种预测反应堆压力容器焊缝热老化晶界p偏析的方法 | |
CN109800458B (zh) | 一种评估反应堆安全运行的方法 | |
JP5483385B2 (ja) | 原子力プラントの管理方法 | |
Muránsky et al. | Assessment of welding-induced plasticity via electron backscatter diffraction | |
Alsit et al. | Experimental and numerical investigation of stress corrosion cracking in AISI 4340 steel under near neutral solution conditions | |
RU2534045C1 (ru) | Способ прогнозирования ресурсоспособности сталей корпусов реакторов ввэр-1000 | |
Leary et al. | Critical Review of Strain Measures for Characterisation of Fatigue Damage in ASME Section III Fatigue Assessments | |
Limbadri et al. | Effect of processing routes on orientation-dependent tensile flow behavior of Zircaloy-4 at elevated temperatures | |
Smith et al. | Assessment of mechanical properties and microstructure characterizing techniques in their ability to quantify amount of cold work in 316L alloy | |
Shabakhty et al. | Investigation on fatigue damage calibration factors in offshore structures | |
Coult et al. | Results from environmentally-assisted short crack fatigue testing on austenitic stainless steels | |
Kim et al. | Ex situ and in situ characterization of stress corrosion cracking of nickel-base alloys at high temperature |
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 |