JP2020016542A - Probabilistic destructive evaluation device, method, and program - Google Patents

Probabilistic destructive evaluation device, method, and program Download PDF

Info

Publication number
JP2020016542A
JP2020016542A JP2018139480A JP2018139480A JP2020016542A JP 2020016542 A JP2020016542 A JP 2020016542A JP 2018139480 A JP2018139480 A JP 2018139480A JP 2018139480 A JP2018139480 A JP 2018139480A JP 2020016542 A JP2020016542 A JP 2020016542A
Authority
JP
Japan
Prior art keywords
size
crack
probability distribution
time
measured
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
JP2018139480A
Other languages
Japanese (ja)
Inventor
陽平 大野
Yohei Ono
陽平 大野
板谷　雅雄
Masao Itaya
雅雄 板谷
琢矢 小川
Takuya Ogawa
琢矢 小川
貴広 林
Takahiro Hayashi
貴広 林
利之 斎藤
Toshiyuki Saito
利之 斎藤
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.)
Toshiba Corp
Toshiba Energy Systems and Solutions Corp
Original Assignee
Toshiba Corp
Toshiba Energy Systems and Solutions Corp
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 Toshiba Corp, Toshiba Energy Systems and Solutions Corp filed Critical Toshiba Corp
Priority to JP2018139480A priority Critical patent/JP2020016542A/en
Publication of JP2020016542A publication Critical patent/JP2020016542A/en
Pending legal-status Critical Current

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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

To provide a probabilistic destructive evaluation technology that can perform a soundness evaluation having a reduced excessive maintainability even when the development of an analysis code is insufficient.SOLUTION: In this technology, a theory size a(t) of a crack proceeding is calculated on the basis of a probability distribution 21of an initial size of the crack existing in a reactor core internal structure before an in-service period t, an actual measurement size 23of the crack measured by non-destructive inspection at j-th time tduring the in-service period is kept, a shift amount μbetween the theory size a(t) and the actual measurement size 23at j-th time tis calculated, a probability distribution 21of the initial size is modified on the basis of the shift amount μ, the theory size a(t) of the crack proceeding with lapse time is re-calculated on the basis of a probability distribution 31of the modified initial size.SELECTED DRAWING: Figure 2

Description

本発明の実施形態は、中性子照射を受ける構造物の確率論的破壊評価技術に関する。   Embodiments of the present invention relate to a probabilistic fracture assessment technique for a structure that is subjected to neutron irradiation.

金属材料は、中性子照射を受けると破壊靱性値が低下することが知られている。中性子照射を受ける構造物として、例えば原子力プラントの原子炉圧力容器内に設置された炉内構造物が挙げられる。   It is known that the fracture toughness value of a metallic material decreases when it is subjected to neutron irradiation. As a structure to be subjected to neutron irradiation, there is, for example, an in-reactor structure installed in a reactor pressure vessel of a nuclear power plant.

原子力プラントは、高経年化に伴い、保全の重要性が高まっている。このため、供用期間中検査（ＩＳＩ：In-service Inspection）と呼ばれる定期検査を定期的に行い、設備の維持に関する技術基準への適合状態が満足されていることの確認が行なわれる。そして、このＩＳＩにおいて、き裂が発見された場合は、このき裂の進展を予測し設備の構造上の健全性評価が行なわれる。   As nuclear power plants age, the importance of maintenance is increasing. For this reason, a periodic inspection called an in-service inspection (ISI: In-service Inspection) is periodically performed, and it is confirmed that the state of conformity with the technical standards for maintenance of the equipment is satisfied. If a crack is found in the ISI, the propagation of the crack is predicted and the structural integrity of the facility is evaluated.

ところで、原子力プラントの炉内構造物に対する従来の健全性評価は、決定論的な評価が行われてきた。しかし、この決定論的な評価による場合、き裂の大きさや材料の破壊強度といった複数の影響パラメータに関し、それぞれ最も厳しい側のデータを採用して評価を行う。このため決定論的な評価は、全体の安全率が高くなる一方で、過度に保守的な評価となってしまう。   By the way, the conventional soundness evaluation of the internal structure of a nuclear power plant has been deterministic. However, in the case of this deterministic evaluation, evaluation is performed by using data on the most severe side for each of a plurality of influence parameters such as a crack size and a material breaking strength. For this reason, deterministic evaluation is an excessively conservative evaluation while increasing the overall safety factor.

そこで、確率論的破壊力学（ＰＦＭ：Probabilistic Fracture Mechanics）の理論を導入して原子力プラントの健全性評価を行うことが検討されている。ＰＦＭでは、健全性評価に関わるほとんどの影響パラメータを確率変数とし、その固有の確率分布に基づいて評価を行う。このためＰＦＭでは、影響パラメータのばらつきや不確かさを考慮して、より現実に近い評価モデルを確立し健全性評価を行なえる。   Therefore, it has been studied to evaluate the integrity of a nuclear plant by introducing the theory of probabilistic fracture mechanics (PFM). In PFM, most of the influence parameters related to soundness evaluation are used as random variables, and evaluation is performed based on the unique probability distribution. For this reason, the PFM can establish a more realistic evaluation model and perform soundness evaluation in consideration of the variation and uncertainty of the influence parameter.

そしてＰＦＭでは、決定論のように最も厳しい側のデータを用いて解析する必要が無くなり、極端に保守的な評価と成ることが避けられる。さらにＰＦＭは、決定論的な評価に比べ、非破壊検査や保全計画等で得られる多数の影響パラメータの効果を、健全性評価に論理的に反映させることができる。   In PFM, there is no need to perform analysis using the data on the most severe side as in determinism, and extremely conservative evaluation can be avoided. Furthermore, PFM can logically reflect the effects of a large number of influence parameters obtained in non-destructive inspection, maintenance planning, and the like in the soundness evaluation as compared with deterministic evaluation.

特開２００７−１９８８３８号公報JP 2007-198838 A

しかし、ＰＦＭの解析コードでは、確率変数として定義したパラメータの確率密度関数を数式としてインプットする必要がある。しかし、この確率密度関数を決定するためには膨大なサンプルデータが必要となるため、多くのパラメータにおいて確率密度関数を精度良く決定することが困難である。   However, in the PFM analysis code, it is necessary to input a probability density function of a parameter defined as a random variable as a mathematical expression. However, since a large amount of sample data is required to determine the probability density function, it is difficult to accurately determine the probability density function for many parameters.

ところで、原子力プラントの炉内構造物は、その内側に配置された炉心燃料から中性子照射を受ける。この炉内構造物の内表面に照射された中性子束は、構造物の材料内部で減衰する。このため、炉内構造物の中性子照射量は、内表面が大きく外表面が小さくなるように分布し、この中性子照射量の分布に依存して、炉内構造物の破壊靭性値の低下が進行していく。   By the way, a reactor internal structure of a nuclear power plant is subjected to neutron irradiation from a core fuel disposed inside. The neutron flux irradiated on the inner surface of the reactor internal structure is attenuated inside the material of the structure. For this reason, the neutron irradiation dose in the reactor internal structure is distributed so that the inner surface is large and the outer surface is small, and the fracture toughness value of the furnace internal structure decreases depending on the distribution of the neutron irradiation amount. I will do it.

照射誘起応力腐食割れ（ＩＡＳＣＣ：Irradiation Assisted Stress Corrosion Cracking）の進展に着目した炉内構造物の健全性評価に関しては、この中性子照射量の分布を考慮して、加速するき裂進展速度を決定する必要がある。しかし、炉内構造物の中性子束分布を考慮したＰＦＭの解析コードの開発は未整備である。   Regarding the soundness evaluation of reactor internals focusing on the development of Irradiation Assisted Stress Corrosion Cracking (IASCC), determine the accelerated crack growth rate in consideration of the distribution of neutron irradiation. There is a need. However, the development of PFM analysis code considering the neutron flux distribution in the reactor internals has not been developed yet.

本発明の実施形態はこのような事情を考慮してなされたもので、解析コードの開発が未整備であっても、過剰な保守性を解消した健全性評価を実施できる確率論的破壊評価技術を提供することを目的とする。   The embodiment of the present invention has been made in view of such circumstances, and even if the development of the analysis code is not yet developed, a probabilistic destruction evaluation technology capable of performing a soundness evaluation that eliminates excessive maintainability. The purpose is to provide.

確率論的破壊評価装置において、供用期間前の炉内構造物に存在するき裂の初期サイズの確率分布に基づいて、供用時間の経過とともに進展するき裂の理論サイズを計算する第１計算部と、供用期間中の第ｊ時点（ｊ：自然数）において非破壊検査により測定されたき裂の実測サイズを保持するデータ保持部と、前記ｊ時点における前記理論サイズと前記実測サイズとの間のシフト量を演算する演算部と、前記シフト量に基づいて前記初期サイズの確率分布を修正する修正部と、前記修正された初期サイズの確率分布に基づいて、経過時間とともに進展するき裂の理論サイズを再計算する第２計算部と、を備えることを特徴とする。   A first calculation unit that calculates a theoretical size of a crack that propagates with a lapse of service time based on a probability distribution of an initial size of a crack existing in a reactor internal structure before a service period in the stochastic fracture evaluation device. A data holding unit for holding the measured size of the crack measured by the nondestructive inspection at the j-th time point (j: natural number) during the service period; and a shift between the theoretical size and the measured size at the j-time point. A calculation unit for calculating the amount, a correction unit for correcting the probability distribution of the initial size based on the shift amount, and a theoretical size of a crack that progresses with elapsed time based on the corrected probability distribution of the initial size. And a second calculating unit for recalculating.

本発明の実施形態により、解析コードの開発が未整備であっても、過剰な保守性を解消した健全性評価を実施できる確率論的破壊評価技術が提供される。   According to the embodiment of the present invention, a probabilistic destruction evaluation technique capable of performing a soundness evaluation that eliminates excessive maintainability even if development of an analysis code is undeveloped is provided.

本発明の実施形態に係る確率論的破壊評価装置を示すブロック図。FIG. 1 is a block diagram showing a probabilistic fracture evaluation device according to an embodiment of the present invention. プラントの供用時間に対し、進展するき裂のサイズをシミュレーションしたグラフ。The graph which simulated the size of the growing crack with respect to the operation time of a plant. 計算部及び演算部で用いられる演算式の説明図。FIG. 4 is an explanatory diagram of an arithmetic expression used in a calculation unit and an arithmetic unit. 解析対象として二次元き裂を有する平板構造物の概念図。The conceptual diagram of the flat-plate structure which has a two-dimensional crack as an analysis object. 実施形態に係る確率論的破壊評価方法及び確率論的破壊評価プログラムを説明するフローチャート。4 is a flowchart illustrating a probabilistic destruction evaluation method and a probabilistic destruction evaluation program according to the embodiment.

以下、本発明の実施形態を添付図面に基づいて説明する。図１は本発明の実施形態に係る確率論的破壊評価装置１０を示すブロック図である。図２はプラントの供用時間（横軸）に対し進展するき裂のサイズ（縦軸）をシミュレーションしたグラフである。   Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings. FIG. 1 is a block diagram showing a probabilistic fracture evaluation device 10 according to an embodiment of the present invention. FIG. 2 is a graph simulating the crack size (vertical axis) that grows with respect to the operating time of the plant (horizontal axis).

図１に示すように（適宜、図２参照）、確率論的破壊評価装置１０は、供用期間前ｔ0の炉内構造物に存在するき裂の初期サイズの確率分布２１0に基づいて供用時間ｔの経過とともに進展するき裂の理論サイズａ0(ｔ)を計算する第１計算部１１と、供用期間中の第ｊ時点ｔj（ｊ：自然数）において非破壊検査により測定されたき裂の実測サイズ２３jを保持するデータ保持部２０と、ｊ時点ｔjにおける理論サイズａ0(ｔ)と実測サイズ２３jとの間のシフト量μjを演算する演算部１６と、このシフト量μjに基づいて初期サイズの確率分布２１0を修正する修正部１５と、この修正された初期サイズの確率分布３１0に基づいて経過時間とともに進展するき裂の理論サイズａj(ｔ)を再計算する第２計算部１２と、を備えている。なお、供用開始後であっても初回の検査前、例えば供用開始後1年の時点をｔ0として、このｔ0時点におけるき裂サイズを初期サイズと読替えてもよい。 As shown in FIG. 1 (as appropriate, see FIG. 2), the probabilistic failure assessment device 10, in service on the basis of the probability distribution 21 0 of the initial size of a crack is present in the furnace structure of prior service period t 0 A first calculation unit 11 for calculating a theoretical size a 0 (t) of a crack that grows with the passage of time t; and a crack measured by a nondestructive inspection at a j-th time point t j (j: natural number) during the service period. a data holding unit 20 for holding the actual size 23 j, an arithmetic unit 16 for calculating a shift amount mu j between the theoretical size at time j t j a 0 (t) and the actual size 23 j, the shift amount a correction unit 15 for correcting the probability distribution 21 0 of the initial size based on mu j, the modified initial size of the probability distribution 31 0 theory crack to progress over time based on the size a j (t) And a second calculator 12 for recalculating. As initial test before even after the start of service, for example, the time point of start of service one year after t 0, may be read as the initial size of the crack size can in this t 0 point.

さらに確率論的破壊評価装置１０は、再計算された理論サイズａj(ｔ)の確率分布の有効上限値を時間軸方向に外挿した外挿線３２jが予め定められた許容限界値２４に達する供用限界時点Ｌjを認定する認定部１８を備えている。 Further, the probabilistic fracture evaluation apparatus 10 sets an extrapolation line 32 j obtained by extrapolating the recalculated effective upper limit value of the probability distribution of the theoretical size a j (t) in the time axis direction to a predetermined allowable limit value 24. A certification unit 18 for certifying the service limit time point Lj at which the service limit time is reached.

さらに確率論的破壊評価装置１０において、データ保持部２０は、供用期間中の第ｊ＋１時点ｔj+1において非破壊検査により測定されたき裂の実測サイズ２３j+1を保持し、演算部１６は第ｊ＋１時点ｔj+1における実測サイズ２３j+1と前回の第ｊ時点ｔjで再計算された理論サイズａj(ｔ)との間のシフト量μj+1を再演算し、修正部１５はこの再演算されたシフト量μj+1に基づいて前回の第ｊ時点ｔjで修正した初期サイズの確率分布３１0をさらに修正し、第２計算部１２はｊ＋１時点ｔj+1で修正された初期サイズの確率分布３３0に基づいて経過時間とともに進展するき裂の理論サイズａj+1(ｔ)を再計算する。 Further, in the probabilistic fracture evaluation apparatus 10, the data holding unit 20 holds the actual size 23 j + 1 of the crack measured by the nondestructive inspection at the j + 1 time point t j + 1 during the service period, and is recalculated shift amount mu j + 1 between the re-calculated theoretical size a j (t) in the (j + 1) time t j actual size 23 in + 1 j + 1 and the previous j-th time point t j, correction unit 15 will further modify the probability distribution 31 0 of the initial size that is fixed in the j-th time point t j of previous based on the shift amount mu j + 1 which is the re-calculation, the second calculating unit 12 j + 1 time t j recalculate the theory of crack progressing size a j + 1 (t) over time based on a probability distribution 33 0 of the modified initial size at +1.

き裂の実測サイズ２３jは、供用期間中の第ｊ時点ｔjにおいて非破壊検査により測定されるものである。この非破壊検査を行う装置としては、渦電流探傷装置、超音波探傷装置等の非破壊検査に一般的に用いられるものである。 The measured size 23 j of the crack is measured by a non-destructive inspection at the j-th time point t j during the service period. As an apparatus for performing this nondestructive inspection, an apparatus generally used for nondestructive inspection such as an eddy current flaw detector and an ultrasonic flaw detector.

なお、ある時点（図では時点ｔ1）の非破壊検査において、検査装置の検出下限２５を下回る等して、き裂が測定されない場合がある。このような場合、修正部１５はこの時点（図では時点ｔ1）において初期サイズの確率分布２１0の修正を実行せずに、次の時点（図では時点ｔj）で測定されたき裂の実測サイズ２３jに基づいて最初の確率分布２１0の修正を実行する。 In a non-destructive inspection at a certain time point (time point t 1 in the figure), the crack may not be measured due to a lower limit of the detection lower limit 25 of the inspection device. In this case, the correction unit 15 does not execute the correction of the probability distribution 21 0 of the initial size at this point (time point in Figure t 1), the following point Taki measured at (time t j in the figure) Cracks executing the correction of the first probability distribution 21 0 based on the measured size 23 j.

データ保持部２０には、き裂の初期サイズの確率分布２１0、中性子束の確率分布２２、非破壊検査により測定されたき裂の実測サイズ２３j（ｊ：自然数）及び許容限界値２４が、情報データとして保持されている。また、図示を省略するが、データ保持部２０は、応力σmの確率分布も保持する場合もある。 The data holding unit 20 stores a probability distribution 21 0 of the initial size of the crack, a probability distribution 22 of the neutron flux, an actually measured size 23 j (j: natural number) of the crack measured by the nondestructive inspection, and an allowable limit value 24. It is held as information data. Although not shown, the data holding unit 20 may also hold the probability distribution of the stress σ m .

確率論的破壊力学（ＰＦＭ：Probabilistic Fracture Mechanics）の計算は、通常モンテカルロ法を用いて行う。まず、さまざまな形状・寸法の初期サイズのき裂を乱数により発生させる。そして、この発生させたき裂について、評価対象の劣化モードに対応したき裂進展則に従って破壊評価を行う。   Calculation of stochastic fracture mechanics (PFM: Probabilistic Fracture Mechanics) is usually performed using the Monte Carlo method. First, cracks of various sizes and initial sizes are generated by random numbers. Then, a fracture evaluation is performed on the generated crack in accordance with a crack growth rule corresponding to the degradation mode to be evaluated.

図３は、計算部１１，１２で用いられる演算式の説明図である。
き裂の初期サイズの確率分布２１0（図２）は、図３の式（１）で与えられる。ここで、aはき裂サイズ、λsはパラメーター（確率変数）である。このように、小さいサイズのき裂は大きいサイズのき裂よりも存在確率が高い。ここで、０．０５ｍｍ以下のき裂は進展しない可能性が示唆されているため、初期サイズの確率分布２１0に下限を設けている。パラメータであるλsの確率密度関数は図３の式（２）のガンマ分布で与えられる。式（２）中のＫおよびθは任意の定数であり、ガンマ関数Γ(ｋ)は図３の式（３）で与えられる。 FIG. 3 is an explanatory diagram of arithmetic expressions used in the calculation units 11 and 12.
The probability distribution 21 0 (FIG. 2) of the initial size of the crack is given by equation (1) in FIG. Here, a is a crack size, and λ s is a parameter (random variable). Thus, small-sized cracks have a higher probability of existence than large-sized cracks. Here, the following crack 0.05mm for may not progress has been suggested is provided with a lower limit to the probability distribution 21 0 of the initial size. The probability density function of the parameter λ s is given by the gamma distribution of equation (2) in FIG. K and θ in Expression (2) are arbitrary constants, and the gamma function Γ (k) is given by Expression (3) in FIG.

膜応力σｍおよび曲げ応力σｂ等の外部応力は、確率変数として表される。膜応力σmの確率密度関数は図３の式（４）で与えられる。ここで、σは標準偏差、μは中央値である。 External stresses such as the film stress σ m and the bending stress σ b are represented as random variables. The probability density function of the film stress σ m is given by equation (4) in FIG. Here, σ is the standard deviation, and μ is the median.

炉内構造物の板厚方向における中性子束分布は、図３の式（５）のような指数関数で与えられる。ここで、ａは板厚方向の深さ、ｈおよびｆ０は定数である。中性子束の確率分布２２は正規分布で表わされ、式（５）の中性子束ｆtを中央値とし、中央値と任意の定数の積を標準偏差とする。 The neutron flux distribution in the thickness direction of the reactor internals is given by an exponential function as shown in Equation (5) in FIG. Here, a is the thickness direction of the depth, h, and f 0 is a constant. Probability distribution 22 of the neutron flux is represented by a normal distribution, the neutron flux f t of formula (5) and the median, the product of the median and arbitrary constant and the standard deviation.

炉内構造物の一例である上部格子板や炉心シュラウドでは、炉心で発生する中性子の照射により材料の劣化が生じて破壊靱性値が低下する。特に上部格子板では中性子照射量が大きく、高経年化対策における健全性評価が重要視されている。   In an upper lattice plate and a core shroud, which are examples of the reactor internals, the material is deteriorated by irradiation of neutrons generated in the core, and the fracture toughness value is reduced. In particular, the neutron irradiation amount is large in the upper lattice plate, and the importance of soundness evaluation in measures against aging is emphasized.

上部格子板の外表面は、材料中の中性子束の減衰により、中性子照射量が内表面に比べて小さいため、破壊靱性値の低下は内表面より小さい。その結果、上部格子板には、内表面から外表面にかけて破壊靱性値が減少する分布が存在することになる。   On the outer surface of the upper lattice plate, the amount of neutron irradiation is smaller than that of the inner surface due to the attenuation of the neutron flux in the material. As a result, the upper lattice plate has a distribution in which the fracture toughness value decreases from the inner surface to the outer surface.

この上部格子板の内表面にき裂が生じたと仮定し、外部からの荷重によって上部格子板が破壊するか否かを評価する場合、き裂先端の破壊靱性値の分布を考慮して評価する必要がある。また、照射誘起応力腐食割れ（ＩＡＳＣＣ）によるき裂の進展を考慮する場合においても、き裂の進展速度は中性子照射量に依存するため、中性子束の減衰による中性子照射量の分布を考慮して評価する必要がある。   Assuming that a crack has occurred on the inner surface of this upper lattice plate, when evaluating whether or not the upper lattice plate is broken by an external load, the evaluation is made in consideration of the distribution of the fracture toughness value of the crack tip. There is a need. Also, when considering crack propagation due to irradiation-induced stress corrosion cracking (IASCC), the crack growth rate depends on the neutron dose, so the distribution of the neutron dose due to the attenuation of the neutron flux is taken into account. Need to be evaluated.

図４は、解析対象の例としての二次元き裂５１を有する平板構造物５０の概念図である。この平板構造物５０に、引張荷重Ｆと、曲げモーメントＭが付加されている場合を示している。引張荷重Ｆは、二次元き裂５２に垂直の方向に付加されている。二次元き裂５２に対する荷重の作用の形式は、モードＩの引張形式、モードＩＩのせん断形式、モードＩＩＩの面外せん断形式に分類される。なお、き裂に対する荷重の作用の仕方として分類される、モードＩの引張形式、モードＩＩのせん断形式、モードＩＩＩの面外せん断形式のうち、図４は、モードＩの場合を例示しているが、それ以外のモードについても検討される。   FIG. 4 is a conceptual diagram of a flat plate structure 50 having a two-dimensional crack 51 as an example of an analysis target. A case where a tensile load F and a bending moment M are added to the flat plate structure 50 is shown. The tensile load F is applied in a direction perpendicular to the two-dimensional crack 52. The type of action of the load on the two-dimensional crack 52 is classified into a mode I tension type, a mode II shear type, and a mode III out-of-plane shear type. FIG. 4 exemplifies the case of mode I among the mode I tensile mode, mode II shear mode, and mode III out-of-plane shear mode classified as a method of applying a load to the crack. However, other modes are considered.

図３の式（６）（７）は、応力腐食割れの進展予測を表す式である。応力腐食割れのき裂進展速度は、図３の式（６）で表される。ここで、ａはき裂サイズ、ｔは時間、Ｋ(ｔ)は任意時間の応力拡大係数、Ｄとｎは材料定数である。そして応力拡大係数は、図３の式（７）で表される。ここで、σは応力、Ｍは補正係数である。このような応力腐食割れの進展予測式は、き裂に対する荷重の作用の仕方として分類されるモードに対応したものが適用される。なお記載を省略するが、応力腐食割れの進展予測を表す式（６）（７）には、中性子束の確率分布２２の寄与因子も含まれている。   Equations (6) and (7) in FIG. 3 are equations for predicting the development of stress corrosion cracking. The crack growth rate of stress corrosion cracking is represented by equation (6) in FIG. Here, a is a crack size, t is time, K (t) is a stress intensity factor at an arbitrary time, and D and n are material constants. The stress intensity factor is expressed by equation (7) in FIG. Here, σ is a stress, and M is a correction coefficient. For such a stress corrosion crack growth prediction formula, one corresponding to a mode classified as a method of applying a load to a crack is applied. Although not shown, the equations (6) and (7) representing the prediction of the development of stress corrosion cracking also include the contributing factors of the probability distribution 22 of the neutron flux.

第１計算部１１で計算されるき裂の理論サイズａ0(ｔ)は、き裂の初期サイズの確率分布２１0に基づいて供用時間ｔを変数として計算され、任意の時間ｔにおけるき裂サイズの確率分布（図示される符号２１1，２１j，２１j+1はその一部）を与える。また第２計算部１２で計算されるき裂の理論サイズａj(ｔ)は、修正された初期サイズの確率分布３１0に基づいて供用時間ｔを変数として計算され、任意の時間ｔにおけるき裂サイズの確率分布（図示される符号３１j，３１j+1はその一部）を与える。同様に第２計算部１２で計算されるき裂の理論サイズａj+1(ｔ)は、修正された初期サイズの確率分布３３0に基づいて供用時間ｔを変数として計算され、任意の時間ｔにおけるき裂サイズの確率分布（図示される符号３３j+1はその一部）を与える。 Theoretical size a 0 of a crack is calculated by the first calculation unit 11 (t) is calculated as a variable serviced time t based on the probability distribution 21 0 of the initial size of the crack, cracks at any time t The probability distribution of the size (the reference numerals 21 1 , 21 j , and 21 j + 1 are some of them) is given. The second calculation unit 12 theoretical size of Ki is calculated crack in a j (t) is calculated as a variable serviced time t on the basis of the probability distribution 31 0 of the modified initial size, can at any time t The probability distribution of the crack size (the reference numerals 31 j and 31 j + 1 are a part thereof) is given. Similarly theory crack which is calculated by the second calculating section 12 size a j + 1 (t) is calculated serviced time t as a variable based on a probability distribution 33 0 of the modified initial size, any time The probability distribution of the crack size at t (33 j + 1 shown is a part thereof) is given.

演算部１６は、図２に示される時点ｔ1のように、非破壊検査の検出下限２５等の理由により、実測サイズ２３jが測定されずデータ保持部２０に存在しない場合は、シフト量の演算を行わない。そして、データ保持部２０に実測サイズ２３jが登録された後は、これを取得し、初演算であるときは、第１計算部１１から理論サイズａ0(ｔ)を取得する。そして、時間ｔjにおける理論サイズａ0(ｔ)と実測サイズ２３jとのシフト量μ1を演算する。 When the actual measurement size 23 j is not measured and does not exist in the data holding unit 20 due to the detection lower limit 25 of the nondestructive inspection or the like, as in the time t 1 shown in FIG. Do not perform calculations. After the actual size 23 j is registered in the data holding unit 20, the actual size 23 j is acquired. If the actual operation is the first operation, the theoretical size a 0 (t) is acquired from the first calculation unit 11. Then, the shift amount μ 1 between the theoretical size a 0 (t) and the measured size 23 j at time t j is calculated.

そして演算部１６における二回目以降のシフト量μj+1の演算は、データ保持部２０から実測サイズ２３j+1を取得し、第２計算部１２から理論サイズａj(ｔ)を取得する。そして、時間ｔj+1における理論サイズａj(ｔ)と実測サイズ２３j+1とのシフト量μj+1を演算する。 In the second and subsequent calculations of the shift amount μ j + 1 by the calculation unit 16, the measured size 23 j + 1 is obtained from the data holding unit 20 and the theoretical size a j (t) is obtained from the second calculation unit 12. . Then, the shift amount μ j + 1 between the theoretical size a j (t) at the time t j + 1 and the measured size 23 j + 1 is calculated.

修正部１５は、演算部１６から取得したシフト量μjに基づいて、初期サイズの確率分布２１0を修正して確率分布３１0に更新する。なお、次回演算されたシフト量μj+1による更新対象は、前回の時点ｔjに更新された初期サイズの確率分布３１0であり、これを確率分布３３0に更新する。 Correction unit 15, based on the shift amount mu j acquired from the operation unit 16, and updates the probability distribution 31 0 modify the probability distribution 21 0 of the initial size. Incidentally, updated by the next calculated shift amount mu j + 1 is the probability distribution 31 0 of the initial size that is updated in the previous time point t j, update this to the probability distribution 33 0.

この初期サイズの確率分布２１0の更新は、図３の式（８）で表されるベイズの定理を用いる。ここで、aは初期き裂サイズ、Ｘは検査結果、Ｐ(ａ)は初期き裂サイズａが得られる確率、Ｐ(Ｘ｜ａ)は尤度関数（初期き裂サイズがａの場合に検査結果Ｘが得られる確率）、Ｐ(ａ｜Ｘ)は検査結果Ｘが得られた場合のき裂サイズのａが得られる確率である。 Updating of the probability distribution 21 0 of the initial size, using Bayes' theorem represented by Equation 3 (8). Here, a is the initial crack size, X is the inspection result, P (a) is the probability of obtaining the initial crack size a, and P (X | a) is the likelihood function (when the initial crack size is a, The probability of obtaining the inspection result X) and P (a | X) are the probability of obtaining the crack size a when the inspection result X is obtained.

修正部１５で修正・更新された初期サイズの確率分布３１0，３３0は、第２計算部１２に送られる。そしてこの第２計算部１２において、供用時間ｔを変数としてき裂の理論サイズａj(ｔ)，ａj+1(ｔ)が再計算される。 Probability distribution 31 0, 33 0 of revision and update initial size correction unit 15 is sent to the second computing section 12. Then, in the second calculation unit 12, the theoretical sizes a j (t) and a j + 1 (t) of the crack are recalculated using the service time t as a variable.

有効上限値外挿部１７（図１）は、図２に示す理論サイズａ0(ｔ)，ａj(ｔ)，ａj+1(ｔ)の各々において、それぞれの確率分布の有効上限値を決定し、さらに連続して外挿させた外挿線３２0，３２j，３２j+1を生成する。この有効上限値は、確率分布の有意水準から決定されるが、その決定方法に限定はない。 The effective upper limit extrapolation unit 17 (FIG. 1) calculates the effective upper limit of the probability distribution for each of the theoretical sizes a 0 (t), a j (t), and a j + 1 (t) shown in FIG. Is determined, and extrapolated lines 32 0 , 32 j , and 32 j + 1 which are successively extrapolated are generated. The effective upper limit is determined from the significance level of the probability distribution, but there is no limitation on the method of determining the upper limit.

認定部１８（図１）は、理論サイズａj(ｔ)の有効上限値の外挿線３２jを時間軸方向に外挿した値が予め定められた許容限界値２４に達する供用限界時点Ｌjを認定する。許容限界値２４は、炉内構造物に存在するき裂が進展してこの炉内構造物が破損に至る臨界レベルである。定期検査が繰り返される度に、供用限界時点Ｌ0，Ｌj，Ｌj+1が延長され、健全性評価における過剰な保守性が緩和されていくのが判る。 The certifying unit 18 (FIG. 1) determines the service limit time point L at which the value obtained by extrapolating the extrapolation line 32 j of the effective upper limit value of the theoretical size a j (t) in the time axis direction reaches the predetermined allowable limit value 24. j is certified. The allowable limit value 24 is a critical level at which a crack existing in the furnace internal structure grows and the furnace internal structure is damaged. It can be seen that the service limit time points L 0 , L j , L j + 1 are extended each time the periodic inspection is repeated, and excessive maintainability in the soundness evaluation is eased.

図５は実施形態に係る確率論的破壊評価方法及び確率論的破壊評価プログラムを説明するフローチャートである（適宜、図２参照）。
ます供用期間前ｔ0の炉内構造物に存在するき裂の初期サイズの確率分布２１0のデータを取得する（Ｓ１１）。次にこの確率分布２１0のデータに基づいて供用時間ｔの経過とともに進展するき裂の理論サイズａ0(ｔ)を計算する（Ｓ１２）。 FIG. 5 is a flowchart illustrating a probabilistic fracture evaluation method and a probabilistic fracture evaluation program according to the embodiment (see FIG. 2 as appropriate).
Service period to get the data of the probability distribution 21 0 of the initial size of the crack that is present in the furnace structure before t 0 (S11). Then calculated theoretical size a 0 of a crack to progress (t) with the lapse of-service time t on the basis of the data of the probability distribution 21 0 (S12).

ｊ＝１として（Ｓ１３）、供用期間中の第ｊ時点ｔjにおいて測定された非破壊検査によるき裂の実測サイズ２３jを保持する（Ｓ１４ Ｙｅｓ）。なお、非破壊検査機器の検出下限２５にかかる等して、き裂の実測サイズ２３jを保持できなかった場合は、次回の第ｊ＋１時点ｔj+1の非破壊検査で再試行する（Ｓ１４ Ｎｏ）。 As j = 1 (S13), the actual size 23 j of the crack measured by the non-destructive inspection measured at the j-th time point t j during the service period is held (S14 Yes). If the actually measured crack size 23 j cannot be retained due to the detection lower limit 25 of the non-destructive inspection device or the like, a retry is performed in the next non-destructive inspection at the j + 1 time point t j + 1 (S 14). No).

ｊ時点ｔjにおける理論サイズａ0(ｔ)と実測サイズ２３jとの間のシフト量μjを演算し（Ｓ１５）、このシフト量μjに基づいて初期サイズの確率分布２１0を修正する（Ｓ１６）。そして、修正された初期サイズの確率分布３１0に基づいて経過時間とともに進展するき裂の理論サイズａj(ｔ)を再計算する（Ｓ１７）。 calculates the shift amount mu j between the theoretical size a 0 (t) and the actual size 23 j at time j t j (S15), corrects the probability distribution 21 0 of the initial size based on the shift amount mu j (S16). Then, recalculating the modified initial size of the probability distribution 31 0 theory of a crack develops over time based on the size a j (t) (S17) .

再計算された理論サイズａj(ｔ)の確率分布の有効上限値を時間軸方向に外挿し、外挿線３２jを生成する（Ｓ１８）。この外挿線３２jが予め定められた許容限界値２４に達する供用限界時点Ｌjを認定する（Ｓ１９）。認定された供用限界時点Ｌjに基づいて、プラント供用をこのまま継続することが適合と判定されれば、（Ｓ１４）に戻って次の時点で非破壊検査が実施され（Ｓ２０ Ｎｏ）、不適合と判定されれば確率論的破壊評価の処理は終了する（Ｓ２０ Ｙｅｓ、ＥＮＤ）。 The effective upper limit of the recalculated probability distribution of the theoretical size a j (t) is extrapolated in the time axis direction to generate an extrapolated line 32 j (S18). The extrapolation line 32 j is certified serviced limit when L j to reach the permissible limit value 24 predetermined (S19). Based on the approved serviced limit point L j, if it is determined that the match is to continue plant serviced this state, (S14) to return the non-destructive inspection is performed at the next time point (S20 No), incompatible with If it is determined, the probabilistic destruction evaluation process ends (S20 Yes, END).

以上述べた少なくともひとつの実施形態の確率論的破壊評価装置によれば、供用期間中に実施した非破壊検査の測定結果を利用して、解析コードの開発が未整備であっても、過剰な保守性を解消した健全性評価を実施することが可能となる。   According to the probabilistic destruction evaluation device of at least one embodiment described above, even if the development of the analysis code is not yet developed, the excessive use of the measurement result of the nondestructive inspection performed during the service period It is possible to carry out soundness evaluation that eliminates maintainability.

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。   Although several embodiments of the present invention have been described, these embodiments are provided by way of example and are not intended to limit the scope of the invention. These embodiments can be implemented in other various forms, and various omissions, replacements, changes, and combinations can be made without departing from the gist of the invention. These embodiments and their modifications are included in the scope and gist of the invention, and are also included in the invention described in the claims and equivalents thereof.

確率論的破壊評価装置の構成要素は、コンピュータのプロセッサで実現することも可能であり、確率論的破壊評価プログラムにより動作させることが可能である。   The components of the stochastic destruction evaluation device can be realized by a processor of a computer, and can be operated by a probabilistic destruction evaluation program.

以上説明した確率論的破壊評価装置は、専用のチップ、ＦＰＧＡ（Field Programmable Gate Array）、ＧＰＵ（Graphics Processing Unit）、又はＣＰＵ（Central Processing Unit）などのプロセッサを高集積化させた制御装置と、ＲＯＭ（Read Only Memory）やＲＡＭ（Random Access Memory）などの記憶装置と、ＨＤＤ（Hard Disk Drive）やＳＳＤ（Solid State Drive）などの外部記憶装置と、ディスプレイなどの表示装置と、マウスやキーボードなどの入力装置と、通信Ｉ／Ｆとを、備えており、通常のコンピュータを利用したハードウェア構成で実現できる。   The probabilistic destruction evaluation device described above is a control device in which a processor such as a dedicated chip, an FPGA (Field Programmable Gate Array), a GPU (Graphics Processing Unit), or a CPU (Central Processing Unit) is highly integrated; Storage devices such as ROM (Read Only Memory) and RAM (Random Access Memory); external storage devices such as HDD (Hard Disk Drive) and SSD (Solid State Drive); display devices such as displays; mice and keyboards And a communication I / F, and can be realized by a hardware configuration using a normal computer.

また確率論的破壊評価装置で実行されるプログラムは、ＲＯＭ等に予め組み込んで提供される。もしくは、このプログラムは、インストール可能な形式又は実行可能な形式のファイルでＣＤ−ＲＯＭ、ＣＤ−Ｒ、メモリカード、ＤＶＤ、フレキシブルディスク（ＦＤ）等のコンピュータで読み取り可能な記憶媒体に記憶されて提供するようにしてもよい。   Further, a program executed by the probabilistic destruction evaluation device is provided by being incorporated in a ROM or the like in advance. Alternatively, the program is provided as a file in an installable format or an executable format stored in a computer-readable storage medium such as a CD-ROM, a CD-R, a memory card, a DVD, and a flexible disk (FD). You may make it.

また、本実施形態に係る確率論的破壊評価装置で実行されるプログラムは、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせて提供するようにしてもよい。また、装置は、構成要素の各機能を独立して発揮する別々のモジュールを、ネットワーク又は専用線で相互に接続し、組み合わせて構成することもできる。   In addition, the program executed by the probabilistic destruction evaluation apparatus according to the present embodiment may be stored on a computer connected to a network such as the Internet, and provided by being downloaded via the network. In addition, the device may be configured by combining different modules that independently perform the functions of the constituent elements with each other via a network or a dedicated line and combining them.

１０…確率論的破壊評価装置、１１…第１計算部、１２…第２計算部、１５…シフト修正部（修正部）、１６…シフト量演算部（演算部）、１７…有効上限値外挿部、１８…供用限界時点認定部（認定部）、２０…データ保持部、２１…き裂の初期サイズの確率分布、２２…中性子束の確率分布、２３…き裂の実測サイズ、２４…許容限界値、２５…検出下限、３１，３３…理論サイズの確率分布、３２…外挿線。   DESCRIPTION OF SYMBOLS 10 ... Probabilistic destruction evaluation apparatus, 11 ... 1st calculation part, 12 ... 2nd calculation part, 15 ... Shift correction part (correction part), 16 ... Shift amount calculation part (calculation part), 17 ... Outside effective upper limit Insertion part, 18 ... service limit time certification part (certification part), 20 ... data holding part, 21 ... probability distribution of initial size of crack, 22 ... probability distribution of neutron flux, 23 ... measured size of crack, 24 ... Tolerable limit value, 25: lower detection limit, 31, 33: probability distribution of theoretical size, 32: extrapolation line.

Claims (7)

供用期間前の炉内構造物に存在するき裂の初期サイズの確率分布に基づいて、供用時間の経過とともに進展するき裂の理論サイズを計算する第１計算部と、
供用期間中の第ｊ時点（ｊ：自然数）において非破壊検査により測定されたき裂の実測サイズを保持するデータ保持部と、
前記ｊ時点における前記理論サイズと前記実測サイズとの間のシフト量を演算する演算部と、
前記シフト量に基づいて前記初期サイズの確率分布を修正する修正部と、
前記修正された初期サイズの確率分布に基づいて、経過時間とともに進展するき裂の理論サイズを再計算する第２計算部と、を備えることを特徴とする確率論的破壊評価装置。 A first calculation unit that calculates a theoretical size of a crack that propagates with a lapse of service time, based on a probability distribution of an initial size of a crack existing in the in-furnace structure before the service period;
A data holding unit that holds an actually measured size of the crack measured by the nondestructive inspection at a j-th time point (j: a natural number) during the service period;
A calculating unit that calculates a shift amount between the theoretical size and the measured size at the time point j;
A correction unit that corrects the probability distribution of the initial size based on the shift amount,
A second calculating unit that recalculates the theoretical size of the crack that grows with the lapse of time based on the corrected probability distribution of the initial size. 請求項１に記載の確率論的破壊評価装置において、
前記再計算された前記理論サイズの確率分布の有効上限値を時間軸方向に外挿した値が予め定められた許容限界値に達する供用限界時点を認定する認定部を備える確率論的破壊評価装置。 The stochastic fracture evaluation device according to claim 1,
Probabilistic destruction evaluation device comprising a certification unit for certifying a service limit point at which a value obtained by extrapolating in the time axis direction the effective upper limit value of the probability distribution of the recalculated theoretical size reaches a predetermined allowable limit value . 請求項１又は請求項２に記載の確率論的破壊評価装置において、
前記データ保持部は、供用期間中の第ｊ＋１時点において非破壊検査により測定されたき裂の実測サイズを保持し、
前記演算部は、前記第ｊ＋１時点における前記実測サイズと前回の第ｊ時点ｔjで再計算された理論サイズとの間のシフト量を再演算し、
前記修正部は、前記再演算されたシフト量に基づいて、前回の第ｊ時点で修正した初期サイズの確率分布をさらに修正し、
前記第２計算部は、前記ｊ＋１時点で修正された初期サイズの確率分布に基づいて、経過時間とともに進展するき裂の理論サイズを再計算する確率論的破壊評価装置。 The stochastic fracture evaluation device according to claim 1 or 2,
The data holding unit holds the measured size of the crack measured by the non-destructive inspection at the j + 1 time point during the service period,
The calculation unit may re-calculates the shift amount between the actual size and re-calculated theoretical size at the time j t j of previous in the (j + 1) th time point,
The correcting unit further corrects the probability distribution of the initial size corrected at the previous j-th point based on the recalculated shift amount,
The stochastic fracture evaluation device, wherein the second calculation unit recalculates the theoretical size of a crack that grows with time based on the probability distribution of the initial size corrected at the time point j + 1. 請求項３に記載の確率論的破壊評価装置において、
前記第ｊ時点の非破壊検査において前記き裂が測定されなかった場合は、
前記修正部は、前記ｊ時点において前記修正を実行せずに、前記ｊ＋１時点で測定された前記き裂の実測サイズに基づいて前記修正を実行する確率論的破壊評価装置。 The stochastic destruction evaluation device according to claim 3,
When the crack is not measured in the non-destructive inspection at the j-th time,
The stochastic fracture evaluation device, wherein the correction unit performs the correction based on an actually measured size of the crack measured at the time j + 1 without performing the correction at the time j. 請求項１から請求項４のいずれか１項に記載の確率論的破壊評価装置において、
前記第１計算部及び／又は前記第２計算部は、前記炉内構造物における中性子束の確率分布にも基づいて、前記理論サイズを計算する確率論的破壊評価装置。 The stochastic fracture evaluation device according to any one of claims 1 to 4,
The first calculation unit and / or the second calculation unit calculates a theoretical size based on a probability distribution of a neutron flux in the reactor internals. 供用期間前の炉内構造物に存在するき裂の初期サイズの確率分布に基づいて、供用時間の経過とともに進展するき裂の理論サイズを計算するステップと、
供用期間中の第ｊ時点（ｊ：自然数）において非破壊検査により測定されたき裂の実測サイズを保持するステップと、
前記ｊ時点における前記理論サイズと前記実測サイズとの間のシフト量を演算するステップと、
前記シフト量に基づいて前記初期サイズの確率分布を修正するステップと、
前記修正された初期サイズの確率分布に基づいて、経過時間とともに進展するき裂の理論サイズを再計算するステップと、を含むことを特徴とする確率論的破壊評価方法。 Based on the probability distribution of the initial size of the cracks present in the furnace internals before the service period, calculating the theoretical size of the cracks that propagate with the passage of service time,
Holding the measured size of the crack measured by the nondestructive inspection at the j-th time point (j: natural number) during the service period;
Calculating a shift amount between the theoretical size and the measured size at the time point j;
Modifying the probability distribution of the initial size based on the shift amount,
Recalculating the theoretical size of a crack that propagates with the passage of time based on the modified probability distribution of the initial size. コンピュータに、
供用期間前の炉内構造物に存在するき裂の初期サイズの確率分布に基づいて、供用時間の経過とともに進展するき裂の理論サイズを計算するステップ、
供用期間中の第ｊ時点（ｊ：自然数）において非破壊検査により測定されたき裂の実測サイズを保持するステップ、
前記ｊ時点における前記理論サイズと前記実測サイズとの間のシフト量を演算するステップ、
前記シフト量に基づいて前記初期サイズの確率分布を修正するステップ、
前記修正された初期サイズの確率分布に基づいて、経過時間とともに進展するき裂の理論サイズを再計算するステップ、を実行させることを特徴とする確率論的破壊評価プログラム。 On the computer,
Calculating the theoretical size of a crack that grows with the passage of service time, based on the probability distribution of the initial size of the cracks present in the furnace internals before the service period,
Holding the measured size of the crack measured by the nondestructive inspection at the j-th time point (j: natural number) during the service period;
Calculating a shift amount between the theoretical size and the measured size at the time point j;
Modifying the probability distribution of the initial size based on the shift amount,
Recalculating the theoretical size of a crack that grows with time based on the corrected initial size probability distribution.
JP2018139480A 2018-07-25 2018-07-25 Probabilistic destructive evaluation device, method, and program Pending JP2020016542A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018139480A JP2020016542A (en) 2018-07-25 2018-07-25 Probabilistic destructive evaluation device, method, and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018139480A JP2020016542A (en) 2018-07-25 2018-07-25 Probabilistic destructive evaluation device, method, and program

Publications (1)

Publication Number Publication Date
JP2020016542A true JP2020016542A (en) 2020-01-30

Family

ID=69579500

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018139480A Pending JP2020016542A (en) 2018-07-25 2018-07-25 Probabilistic destructive evaluation device, method, and program

Country Status (1)

Country Link
JP (1) JP2020016542A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022038840A1 (en) * 2020-08-17 2022-02-24 株式会社日立製作所 Crack progress evaluation device and crack progress evaluation program
JP7512073B2 (en) 2020-04-28 2024-07-08 株式会社東芝 System and method for evaluating and analyzing soundness

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7512073B2 (en) 2020-04-28 2024-07-08 株式会社東芝 System and method for evaluating and analyzing soundness
WO2022038840A1 (en) * 2020-08-17 2022-02-24 株式会社日立製作所 Crack progress evaluation device and crack progress evaluation program

Similar Documents

Publication Publication Date Title
US9792555B2 (en) Probabilistic modeling and sizing of embedded flaws in ultrasonic nondestructive inspections for fatigue damage prognostics and structural integrity assessment
Soliman et al. Life-cycle management of fatigue-sensitive structures integrating inspection information
Kim et al. Cost-based optimum scheduling of inspection and monitoring for fatigue-sensitive structures under uncertainty
Liu et al. Utility and information analysis for optimum inspection of fatigue-sensitive structures
US10275546B2 (en) Method of predicting crack growth and information processing device therefor
Shekari et al. A predictive approach to fitness-for-service assessment of pitting corrosion
Guan et al. Probabilistic modeling and sizing of embedded flaws in ultrasonic non-destructive inspections for fatigue damage prognostics and structural integrity assessment
El Hajj et al. A condition‐based deterioration model for the stochastic dependency of corrosion rate and crack propagation in corroded concrete structures
Keprate et al. Fatigue and fracture degradation inspection of offshore structures and mechanical items: the state of the art
Yuan et al. Probabilistic assessment for concrete spalling in tunnel structures
Jena et al. Differential evolution: an inverse approach for crack detection
Wang et al. Probabilistic analysis of corrosion of reinforcement in RC bridges considering fuzziness and randomness
JP2020016542A (en) Probabilistic destructive evaluation device, method, and program
Šomodíková et al. fib models for modeling of chloride ion ingress and concrete carbonation: Levels of assessment of input parameters
Uchida et al. Improvement of plant reliability based on combining of prediction and inspection of crack growth due to intergranular stress corrosion cracking
JP6746512B2 (en) Destruction evaluation analysis device, destructive evaluation analysis system and destructive evaluation analysis method
Khaleel et al. A model for predicting vessel failure probabilities including the effects of service inspection and flaw sizing errors
JP2017167095A (en) Crack evaluation device and crack evaluation method
JP6636891B2 (en) Destruction evaluation analyzer, destruction evaluation system, and destruction evaluation method
Clayton et al. Thick concrete specimen construction, testing, and preliminary analysis
Innal et al. Treatment of uncertainty in probabilistic risk assessment using Monte Carlo analysis
Li et al. A novel computational approach for assessing system reliability and damage detection delay: Application to fatigue deterioration in offshore structures
Nikbin A unified European approach to high temperature defect assessment code and its incorporation in a knowledge base system
Qin Probabilistic modeling and Bayesian inference of metal-loss corrosion with application in reliability analysis for energy pipelines
JP2005242712A (en) Method and program for working out maintenance plan of plant