WO2020071237A1 - 水素脆化進行度の評価方法及び水素脆化進行度評価装置 - Google Patents

水素脆化進行度の評価方法及び水素脆化進行度評価装置

Info

Publication number
WO2020071237A1
WO2020071237A1 PCT/JP2019/037857 JP2019037857W WO2020071237A1 WO 2020071237 A1 WO2020071237 A1 WO 2020071237A1 JP 2019037857 W JP2019037857 W JP 2019037857W WO 2020071237 A1 WO2020071237 A1 WO 2020071237A1
Authority
WO
WIPO (PCT)
Prior art keywords
hydrogen embrittlement
progress
transition probability
time
rupture
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.)
Ceased
Application number
PCT/JP2019/037857
Other languages
English (en)
French (fr)
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.)
NTT Inc
Original Assignee
Nippon Telegraph and Telephone 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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2020550356A priority Critical patent/JP7108214B2/ja
Priority to US17/275,144 priority patent/US20220042901A1/en
Publication of WO2020071237A1 publication Critical patent/WO2020071237A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N17/00—Investigating resistance of materials to the weather, to corrosion, or to light
    • G01N17/006—Investigating resistance of materials to the weather, to corrosion, or to light of metals
    • E—FIXED CONSTRUCTIONS
    • E04—BUILDING
    • E04C—STRUCTURAL ELEMENTS; BUILDING MATERIALS
    • E04C5/00—Reinforcing elements, e.g. for concrete; Auxiliary elements therefor
    • E04C5/01—Reinforcing elements of metal, e.g. with non-structural coatings
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/20—Metals
    • G01N33/204—Structure thereof, e.g. crystal structure

Definitions

  • the present invention relates to a method and an apparatus for evaluating the degree of progress of hydrogen embrittlement.
  • the conventional method of evaluating the hydrogen embrittlement resistance properties of steel materials based on the rupture time is mainly for the qualitative evaluation between steel materials, and it is possible to know the degree of hydrogen embrittlement progress in the course of fracture. Did not.
  • the present invention has been made in view of the above, and an object of the present invention is to know the degree of hydrogen embrittlement progress in the course of a steel material being broken by hydrogen embrittlement.
  • the method for evaluating the degree of hydrogen embrittlement progress includes a first step of obtaining a cumulative distribution function of rupture time from hydrogen embrittlement rupture test data measured under test conditions assuming reinforcing bars in prestressed concrete; A second step of calculating the transition probability from the cumulative distribution function of the rupture time, and treating the transition probability as a hydrogen embrittlement progress degree And a third step of determining
  • the hydrogen embrittlement progress evaluation apparatus is a cumulative distribution function calculation unit for obtaining a cumulative distribution function of rupture time from hydrogen embrittlement rupture test data measured under test conditions assuming rebar in prestressed concrete, hydrogen
  • the embrittlement rupture is regarded as a stochastic process in which the deterioration state sequentially changes and leads to rupture, a transition probability calculation unit that calculates the transition probability from the cumulative distribution function of the rupture time, and hydrogen embrittlement by integrating the transition probability with time.
  • FIG. 3 is a diagram in which hydrogen embrittlement rupture time data of FIG. 2 is plotted on Weibull probability paper. It is a graph which shows the time change of transition probability. 4 is a graph showing a change over time in the degree of hydrogen embrittlement. 5 is a graph showing a temporal change of a transition probability obtained by another method from the graph of FIG. 4. It is a graph which shows the time change of the hydrogen embrittlement progress calculated using the transition probability calculated
  • FIG. 9 is a graph showing transition probabilities obtained from the degree of hydrogen embrittlement progression in FIG. 8.
  • FIG. 10 is a diagram illustrating a cumulative rupture probability predicted using the transition probabilities in FIG. 9 and a result when an experiment is actually performed under the same conditions. It is a flowchart which shows the flow of a process of the hydrogen embrittlement progress evaluation apparatus in this Embodiment. It is a figure showing an example of hardware constitutions of a hydrogen embrittlement progress evaluation device.
  • FIG. 1 is a functional block diagram illustrating a configuration of a hydrogen embrittlement progress evaluation apparatus 1 according to the present embodiment.
  • the hydrogen embrittlement progress evaluation apparatus 1 shown in FIG. 1 includes a cumulative distribution function calculator 11, a transition probability calculator 12, and a hydrogen embrittlement progress calculator 13.
  • Each unit included in the hydrogen embrittlement progress evaluation apparatus 1 may be configured by a computer including an arithmetic processing unit, a storage device, and the like, and the processing of each unit may be executed by a program.
  • This program is stored in a storage device included in the hydrogen embrittlement progress evaluation apparatus 1, and can be recorded on a recording medium such as a magnetic disk, an optical disk, or a semiconductor memory, or can be provided through a network.
  • the storage device 2 and the display device 3 are connected to the hydrogen embrittlement progress evaluation device 1.
  • the hydrogen embrittlement progress evaluation device 1 may include a storage device 2 and a display device 3.
  • the cumulative distribution function calculation unit 11 calculates the cumulative distribution function of the rupture time from a large number of hydrogen embrittlement rupture test data measured under the test conditions assuming the reinforcing steel in the prestressed concrete.
  • the hydrogen embrittlement rupture test data is stored in the storage device 2.
  • the cumulative distribution function calculation unit 11 reads out the hydrogen embrittlement fracture test data from the storage device 2 and calculates the cumulative distribution function of the fracture time.
  • the hydrogen embrittlement rupture time data shown in FIG. 2 is obtained by repeating a test of measuring the rupture time by charging hydrogen to the steel material while applying a constant tensile stress to the high-strength steel material, and obtaining a large number of rupture times. is there.
  • the rupture times are arranged in ascending order and plotted using the average rank method.
  • the test conditions are as follows.
  • the test piece used was a round bar having a length of 50 cm and a diameter of 7 mm.
  • the tensile stress was 0.9 times the tensile strength (0.9 ⁇ b).
  • a cathode charging method in which a steel material was immersed in an aqueous electrolyte solution and a negative potential was applied was used.
  • electrolyte aqueous solution a 1 mol / L sodium bicarbonate aqueous solution + 0.13 mol / L thiocyanate ammonium aqueous solution was used.
  • the applied potential is -1V @ vs. SSE was used.
  • the cumulative distribution function F (t) is estimated using, for example, general Weibull probability paper.
  • the cumulative distribution function F (t) was calculated for the hydrogen embrittlement rupture time data of FIG. 2 using Weibull probability paper as shown in FIG. 3, the following equation (1) was obtained.
  • the transition probability calculating unit 12 regards the hydrogen embrittlement rupture as a stochastic process in which the deterioration state sequentially transitions to the fracture state, and calculates the transition probability thereof. .
  • the inventors have found that by making the following two assumptions, it is possible to calculate the transition probability of the deterioration state in hydrogen embrittlement, and to estimate the degree of hydrogen embrittlement from the transition probability. .
  • Transition probability is the probability that the deterioration state changes only once per unit time, and represents the average value of the rate of increase in deterioration. That is, it is considered that the value H (t) obtained by integrating the transition probability (rate of increase in deterioration) with time represents the degree of progress of hydrogen embrittlement.
  • transition of the degradation state due to hydrogen embrittlement can be considered to be the same mechanism as the transition of the equipment to the spare equipment due to the failure of the equipment in the standby system. That is, the transition probability ⁇ (t) of hydrogen embrittlement can be obtained by using the expressions (1) and (2).
  • equation (2) cannot be solved analytically, but can be accurately approximated by setting the transition probability ⁇ (t) to the following equation (3).
  • Equation (4) shows a graph of the change over time of the transition probability expressed by the equation (4).
  • n may be any sufficiently large value.
  • the hydrogen embrittlement progress calculating unit 13 calculates the hydrogen embrittlement progress by integrating the transition probability with time. There is no particular definition of the method, and a general method may be used.
  • the hydrogen embrittlement progress calculating unit 13 causes the display device 3 to display the obtained hydrogen embrittlement progress.
  • FIG. 5 shows the time change of the hydrogen embrittlement progress H (t) calculated by the hydrogen embrittlement progress evaluation apparatus 1 of the present embodiment.
  • the hydrogen embrittlement progress evaluation apparatus 1 associates the hydrogen embrittlement rupture test data used in the calculation of the hydrogen embrittlement progress with the life data of the steel material of the actual equipment, and calculates the actual hydrogen embrittlement progress from the obtained hydrogen embrittlement progress.
  • transition probability ⁇ (t) can be obtained more accurately by using the following method.
  • the cumulative distribution function F (t) of the rupture time in equation (1) By substituting the cumulative distribution function F (t) of the rupture time in equation (1) into equation (2) and rearranging the equation for t, the following equation (6) is obtained.
  • a value from 1 to 100 is sequentially substituted for the hydrogen embrittlement progression H (t) in the equation (6), and the value of t corresponding to the value of H (t) is obtained, whereby the change over time of H (t) is obtained.
  • Ask for. ⁇ (t) can be approximately expressed by the slope of H (t).
  • the slope of H (t) is obtained by the difference method using the following equation (7), and is set as ⁇ (t).
  • t1 and t2 are arbitrary times, and t1 ⁇ t2.
  • ⁇ (t) is directly calculated without using an approximation formula, so that ⁇ (t) can be obtained more accurately.
  • FIG. 6 shows a comparison between ⁇ (t) obtained by the above-described method of calculating the slope of H (t) and a method using the approximation formula (4).
  • the result obtained by the method of obtaining the slope of H (t) is shown by a plot, and the result obtained by the method using the approximate expression of Expression (4) is shown by a solid line. From FIG. 6, it can be seen that, particularly in the latter half of the test, the method of obtaining the slope of H (t) has a remarkable difference in accuracy as compared with the method using the approximate expression of Expression (4).
  • FIG. 7 shows the time change of the hydrogen embrittlement progression H (t) calculated using ⁇ (t) obtained by using the method for obtaining the slope of H (t) in the hydrogen embrittlement progression evaluation apparatus 1. .
  • the hydrogen embrittlement when environmental conditions change, for example, the hydrogen embrittlement when the environmental conditions change such that the hydrogen generation condition due to the corrosion reaction or the progress rate of the hydrogen embrittlement changes according to the season. It can predict the degree of progress.
  • the degree of hydrogen embrittlement progress H (t) is obtained for different environmental conditions by the method described above.
  • the degree of hydrogen embrittlement in the case where electrochemical hydrogen charging is performed at 0.025 mA / mm 2 simulating a state in which hydrogen embrittlement progresses rapidly is indicated by a dashed line.
  • the broken line indicates the degree of hydrogen embrittlement when electrochemical hydrogen charging is performed at 0.015 mA / mm 2 simulating a slow state.
  • H (t) indicated by a solid line changes in the same manner as the hydrogen embrittlement progression of 0.025 mA / mm 2 until 240 minutes, and the hydrogen embrittlement of 0.015 mA / mm 2 after 240 minutes. It changes as well as the degree of progress.
  • ⁇ Whether the prediction of hydrogen embrittlement progression when environmental conditions change is correct can be confirmed by a probabilistic simulation (Monte Carlo simulation).
  • the cumulative rupture probability is obtained by repeating the simulation from the start of the test to the rupture under the condition that the rupture occurs.
  • the cumulative fracture probability obtained by the simulation is compared with the result obtained when an experiment is actually performed under the same conditions. If the two values agree with each other, it can be said that the prediction of the degree of hydrogen embrittlement progression is correct.
  • FIG. 9 shows ⁇ (t) obtained from the slope of H (t) in FIG.
  • FIG. 10 shows the cumulative fracture probability obtained by repeating the above simulation 1000 times using ⁇ (t) in FIG. 9 and the result of an experiment actually performed under the same conditions.
  • the simulation result is indicated by a dotted line
  • the experimental result is indicated by a plot.
  • the cumulative rupture probability obtained by the experiment and the cumulative rupture probability obtained by the simulation using ⁇ (t) are almost the same, and it can be seen that the experimental result can be predicted by the simulation. Therefore, it can be said that the degree of hydrogen embrittlement progress when environmental conditions change can be predicted correctly.
  • FIG. 11 is a flowchart showing the flow of processing of the hydrogen embrittlement progress evaluation apparatus 1 in the present embodiment. It is assumed that the hydrogen embrittlement rupture test data measured under the test conditions assuming the reinforcing steel in the prestressed concrete has been acquired and stored in the storage device 2.
  • the cumulative distribution function calculator 11 calculates a cumulative distribution function of the rupture time from the hydrogen embrittlement test data (step S11).
  • the transition probability calculation unit 12 regards the hydrogen embrittlement rupture as a probability process in which the deterioration state sequentially transitions to the fracture, and calculates the transition probability from the cumulative distribution function of the rupture time (step S12).
  • the hydrogen embrittlement progress calculating unit 13 calculates the hydrogen embrittlement progress by integrating the transition probability with time (step S13).
  • the renewal time estimation unit uses the hydrogen embrittlement progress degree based on the life data in the actual equipment of the reinforcing bar used for measuring the hydrogen embrittlement fracture test data. Estimate the renewal time of the actual equipment from the
  • a cumulative distribution function of rupture time is calculated from a large number of hydrogen embrittlement rupture test data measured under test conditions assuming a reinforcing bar in prestressed concrete, and hydrogen embrittlement is calculated.
  • the hydrogen embrittlement progress evaluation apparatus 1 described above includes, for example, a central processing unit (CPU) 901, a memory 902, a storage 903, a communication device 904, and an input device 905 as shown in FIG. , An output device 906 can be used.
  • the hydrogen embrittlement progress evaluation apparatus 1 is realized by the CPU 901 executing a predetermined program loaded on the memory 902.

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Architecture (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Ecology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Environmental Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
  • Investigating And Analyzing Materials By Characteristic Methods (AREA)

Abstract

プレストレストコンクリート内の鉄筋を想定した試験条件により測定された、多数の水素脆化破断試験データから破断時間の累積分布関数を計算し、水素脆化破断を劣化状態が逐次推移して破断状態に至る確率過程とみなし、破断時間の累積分布関数から推移確率を計算し、推移確率を時間で積分して水素脆化進行度を求める。

Description

水素脆化進行度の評価方法及び水素脆化進行度評価装置
　本発明は、水素脆化進行度の評価方法および装置に関する。
　近年、建築資材としてプレストレストコンクリート内の鉄筋に使用される高強度鋼材の需要が高まっている。高強度鋼材は水素脆化感受性が高く、鋼材中に経時的に水素が含まれることにより延性が失われ、強度が著しく低下する。この現象は水素脆化と呼ばれている（非特許文献１参照）。鋼材の水素脆化に関し、鋼材の性能として耐水素脆化特性の評価が行われている。例えば、鋼材に一定の引張応力を付与（定荷重試験）しながら鋼材中に電気化学的に水素チャージし、鋼材の破断時間を測定する評価方法が広く用いられている（非特許文献２参照）。
白神哲夫、「鉄鋼材料における水素脆化」、材料と環境、公益社団法人腐食防食学会、２０１１年、第６０巻、第５号、pp. 236-420 ２０％チオシアン酸アンモニウム溶液中でのＰＣ鋼材の水素脆化試験方法、ＪＳＣＥ　Ｓ　１２０１、社団法人腐食防食協会、２０１２年 塩見弘著、「信頼性工学入門」、改訂二版、丸善株式会社、１９７１年、p. 103
　実際に設備等に適用される鋼材について、鋼材が水素脆化によって破断に至る途中の水素脆化進行度を知ることができれば、鋼材の使用期間中にどの程度水素脆化が進行しているのかを予想でき、鋼材が破断する前の適切なタイミングで設備の更改を行う判断が可能となる。
　しかしながら、従来用いられてきた破断時間による鋼材の耐水素脆化特性の評価方法は、鋼材間の定性的評価が主な目的であり、破断に至る途中の水素脆化進行度を知ることはできなかった。
　本発明は、上記に鑑みてなされたものであり、鋼材が水素脆化によって破断に至る途中の水素脆化進行度を知ることを目的とする。
　本発明に係る水素脆化進行度の評価方法は、プレストレストコンクリート内の鉄筋を想定した試験条件により測定された水素脆化破断試験データから破断時間の累積分布関数を求める第１工程と、水素脆化破断を劣化状態が逐次推移して破断に至る確率過程とみなし、前記破断時間の累積分布関数から推移確率を計算する第２工程と、前記推移確率を時間で積分して水素脆化進行度を求める第３工程と、を有することを特徴とする。
　本発明に係る水素脆化進行度評価装置は、プレストレストコンクリート内の鉄筋を想定した試験条件により測定された水素脆化破断試験データから破断時間の累積分布関数を求める累積分布関数計算部と、水素脆化破断を劣化状態が逐次推移して破断に至る確率過程とみなし、前記破断時間の累積分布関数から推移確率を計算する推移確率計算部と、前記推移確率を時間で積分して水素脆化進行度を求める水素脆化進行度計算部と、を有することを特徴とする。
　本発明によれば、鋼材が水素脆化によって破断に至る途中の水素脆化進行度を知ることができる。
本実施形態の水素脆化進行度評価装置の構成を示す機能ブロック図である。 水素脆化破断時間データの一例を示す図である。 図２の水素脆化破断時間データをワイブル確率紙へプロットした図である。 推移確率の時間変化を示すグラフである。 水素脆化進行度の時間変化を示すグラフである。 図４のグラフと別の方法で求めた推移確率の時間変化を示すグラフである。 別の方法で求めた推移確率を用いて算出した水素脆化進行度の時間変化を示すグラフである。 異なる環境条件の水素脆化進行度を組み合わせて求めた水素脆化進行度を示すグラフである。 図８の水素脆化進行度から求めた推移確率を示すグラフである。 図９の推移確率を用いて予測した累積破断確率と、実際に同条件で実験を行った場合の結果を示す図である。 本実施の形態における水素脆化進行度評価装置の処理の流れを示すフローチャートである。 水素脆化進行度評価装置のハードウェア構成の一例を示す図である。
　以下、本発明の実施の形態について図面を用いて説明する。
　（水素脆化進行度評価装置の構成）
　図１は、本実施の形態における水素脆化進行度評価装置１の構成を示す機能ブロック図である。同図に示す水素脆化進行度評価装置１は、累積分布関数計算部１１、推移確率計算部１２、及び水素脆化進行度計算部１３を備える。水素脆化進行度評価装置１が備える各部は、演算処理装置、記憶装置等を備えたコンピュータにより構成して、各部の処理がプログラムによって実行されるものとしてもよい。このプログラムは水素脆化進行度評価装置１が備える記憶装置に記憶されており、磁気ディスク、光ディスク、半導体メモリ等の記録媒体に記録することも、ネットワークを通して提供することも可能である。
　水素脆化進行度評価装置１には、記憶装置２および表示装置３が接続される。水素脆化進行度評価装置１が記憶装置２および表示装置３を備えてもよい。
　累積分布関数計算部１１は、プレストレストコンクリート内の鉄筋を想定した試験条件により測定された、多数の水素脆化破断試験データから破断時間の累積分布関数を計算する。水素脆化破断試験データは、記憶装置２に記憶させておく。累積分布関数計算部１１は、記憶装置２から水素脆化破断試験データを読み出して破断時間の累積分布関数を計算する。
　水素脆化破断時間データは、例えば図２に示すようなものを用いる。図２に示した水素脆化破断時間データは、高強度鋼材に一定の引張応力を付与しながら鋼材中に水素チャージして破断時間を測定する試験を繰り返し、多数の破断時間を取得したものである。図２では、破断時間を小さい方から順に並べ平均ランク法を用いてプロットしている。試験条件は以下の通りである。試験片は長さ５０ｃｍで直径７ｍｍの丸棒平滑材を用いた。引張応力は引張強さの０．９倍の応力（０．９σｂ）とした。水素チャージ方法は鋼材を電解質水溶液に浸漬させて負電位を印加する陰極チャージ法を用いた。電解質水溶液には１ｍｏｌ／Ｌの炭酸水素ナトリウム水溶液＋０．１３ｍｏｌ／Ｌのチオシアン線アンモニウム水溶液を用いた。印加電位には－１Ｖ　ｖｓ．　ＳＳＥを用いた。試験数はＮ＝５５とした。これらの試験条件はプレストレストコンクリート内の鉄筋を想定しており、鋼材中に含まれる水素量を実際の環境よりも極端に増やして水素脆化を加速させている。
　水素脆化破断時間データから破断時間の累積分布関数を求めるには、特に方法の規定はなく、一般的な方法を用いれば良い。水素脆化の破断時間はワイブル分布に従うと考えられることから、例えば一般的なワイブル確率紙を用いて累積分布関数Ｆ（ｔ）を推定する。図２の水素脆化破断時間データに対し、図３に示すようにワイブル確率紙を用いて累積分布関数Ｆ（ｔ）を計算したところ、以下の（１）式が得られた。
Figure JPOXMLDOC01-appb-M000001
　推移確率計算部１２は、累積分布関数計算部１１の計算した累積分布関数を用いて、水素脆化破断を劣化状態が逐次推移して破断状態に至る確率過程とみなしてその推移確率を計算する。
　発明者らは、以下の２つの仮定を置くことで、水素脆化における劣化状態の推移確率を計算することが可能となり、推移確率から水素脆化の進行度を推定可能であることを見出した。
　仮定１：水素脆化では劣化が時間とともに蓄積し、一定値を超えたところで破断が発生する。
　仮定２：ある劣化状態から次の劣化状態への変化は推移確率λ(t)に従って確率的に発生する。
　推移確率は単位時間あたりに劣化状態が１回だけ推移する確率であり、これは劣化の増加速度の平均値を表している。すなわち、推移確率（劣化の増加速度）を時間で積分した値Ｈ（ｔ）が水素脆化の進行度を表していると考えられる。
　劣化が確率的に進行し、破断に至るという考え方は金属材料の疲労破壊においても用いられているが、疲労破壊の場合は推移確率λを時間変化しない一定値としているため、水素脆化のように推移確率（劣化の増加速度）が時間変化すると考えられる現象に対してはそのまま適用することができない。一方、信頼性工学の分野では、装置が故障した際に順次予備の装置に切り替わっていく待機系と呼ばれるシステムにおいて、装置１個の故障率λ(ｔ)とシステム全体が故障する累積分布関数Ｆ(ｔ)の関係は、予備の装置がｎ個の場合は下記の（２）式で表されることが知られている（非特許文献３参照）。
Figure JPOXMLDOC01-appb-M000002
　水素脆化における劣化状態の推移は、待機系において装置が故障し予備の装置へ推移することと同様の機構であると考えることができる。すなわち、（１）式と（２）式を用いることで、水素脆化の推移確率λ(ｔ)を求めることができる。また、累積分布関数Ｆ（ｔ）がワイブル分布に従うとき、（２）式は解析的に解くことができないが、推移確率λ（ｔ）を以下の（３）式とすることで精度良く近似可能となる。
Figure JPOXMLDOC01-appb-M000003
　上記の方法に従い、劣化状態が１００回推移したところで破断が発生する（ｎ＝１００）とし、（２）式に（３）式を代入し、（１）式の累積分布関数Ｆ（ｔ）に適合するように（３）式のＡ，Ｂ，Ｃの値を求めると、推移確率λ（ｔ）は以下の（４）式のように計算される。図４に（４）式で示される推移確率の時間変化のグラフを示す。
Figure JPOXMLDOC01-appb-M000004
　なお、破断が発生するｎの値は任意に取っているが、ｎの値が十分に大きければｎを変えても推移確率λ（ｔ）の曲線はほぼ相似形となることから、破断に至るまでの水素脆化進行度を知る上ではｎは十分に大きな任意の値で良いこととなる。
　水素脆化進行度計算部１３は、推移確率を時間で積分して水素脆化進行度を求める。特に方法の規定はなく、一般的な方法を用いれば良い。水素脆化進行度計算部１３は、求めた水素脆化進行度を表示装置３に表示させる。
　（４）式から水素脆化進行度Ｈ（ｔ）を計算すると、水素脆化進行度Ｈ（ｔ）は以下の（５）式のように計算される。
Figure JPOXMLDOC01-appb-M000005
　図５に、本実施形態の水素脆化進行度評価装置１により計算した水素脆化進行度Ｈ（ｔ）の時間変化を示す。
　劣化状態が１００回推移したところ（Ｈ（ｔ）＝１００）で破断が発生するとしているため、Ｈ（ｔ）＝１００となるｔ（ｔ＝２．６時間）は平均破断時間を表している。また、図５から破断の９０％の水素脆化進行度（Ｈ（ｔ）＝９０）となる時間ｔを読み取るとｔ＝１．４時間である。
　水素脆化進行度評価装置１は、水素脆化進行度の計算に用いた水素脆化破断試験データと実設備の鋼材の寿命データとを対応付けておき、求めた水素脆化進行度から実設備の更改時期を推定する更改時期推定部を備えてもよい。例えば、今回用いた水素脆化破断時間データの取得に使用した鋼材の実環境での平均寿命が２０年の場合、Ｈ（ｔ）＝１００となるｔ＝２．６時間が実環境での２０年に相当する。水素脆化進行度が破断に至る９０％で実設備を更改する場合、Ｈ（ｔ）＝９０となるｔ＝１．４時間に相当する実環境での約１１年で実設備を更改すればよいことになる。
　（推移確率の別の求め方）
　推移確率λ（ｔ）は以下の方法を用いればさらに正確に求めることができる。（１）式の破断時間の累積分布関数Ｆ（ｔ）を（２）式に代入し、ｔについて式を整理すると以下の（６）式が得られる。
Figure JPOXMLDOC01-appb-M000006
　この（６）式の水素脆化進行度Ｈ（ｔ）に１から１００までの値を順次代入し、Ｈ（ｔ）の値と対応するｔの値を求めることでＨ（ｔ）の経時変化を求める。λ（ｔ）は近似的にＨ（ｔ）の傾きで表すことができる。以下の（７）式を用いて差分法によりＨ（ｔ）の傾きを求めてλ（ｔ）とする。
Figure JPOXMLDOC01-appb-M000007
　ここで、ｔ１、ｔ２は任意の時間であり、ｔ１＜ｔ２である。この方法は近似式を用いずに直接λ（ｔ）を算出しているため、より正確にλ（ｔ）を求めることができる。
　上記のＨ（ｔ）の傾きを求める方法と（４）式の近似式を用いた方法のそれぞれにより求めたλ（ｔ）を比較すると図６のようになる。図６では、Ｈ（ｔ）の傾きを求める方法で求めた結果をプロットで示し、（４）式の近似式を用いた方法で求めた結果を実線で示している。図６から、特に試験後半においては、Ｈ（ｔ）の傾きを求める方法は（４）式の近似式を用いた方法に比べて精度の違いが顕著であることがわかる。
　図７に、水素脆化進行度評価装置１においてＨ（ｔ）の傾きを求める方法を用いて求めたλ（ｔ）を用いて計算した水素脆化進行度Ｈ（ｔ）の時間変化を示す。図７から、Ｈ（ｔ）＝１００となる平均破断時間ｔは２．５時間であることがわかる。また、Ｈ（ｔ）＝９０となる時間ｔは１．２時間であることがわかる。
　（環境条件が変化する場合）
　本実施形態では、環境条件が異なる水素脆化進行度を組み合わせることで、例えば腐食反応による水素発生条件または水素脆化の進行速度が季節により変化するような、環境条件が変化する場合の水素脆化進行度を予測できる。
　まず、これまでに説明した方法により、異なる環境条件について水素脆化進行度Ｈ（ｔ）を求めておく。図８では、水素脆化の進行が速い状態を模擬する０．０２５ｍＡ／ｍｍ２で電気化学的な水素チャージを行った場合の水素脆化進行度を一点鎖線で示し、水素脆化の進行が遅い状態を模擬する０．０１５ｍＡ／ｍｍ２で電気化学的な水素チャージを行った場合の水素脆化進行度を破線で示している。
　試験開始時は水素チャージが０．０２５ｍＡ／ｍｍ２で行われ、試験開始から２４０分経過後に水素チャージが０．０１５ｍＡ／ｍｍ２で行われるという環境条件を考えた場合、Ｈ（ｔ）は図８の実線に示すような挙動となる。図８の例では、実線で示すＨ（ｔ）は２４０分までは０．０２５ｍＡ／ｍｍ２の水素脆化進行度と同様に変化し、２４０分以降は０．０１５ｍＡ／ｍｍ２の水素脆化進行度と同様に変化する。水素脆化進行度計算部１３が、異なる環境条件に基づいて求めた水素脆化進行度を組み合わせることで、環境条件が変化する場合の水素脆化進行度が求められる。
　環境条件が変化する場合の水素脆化進行度の予測が正しいかは、確率的シミュレーション（モンテカルロシミュレーション）により確認することができる。時間ｔからΔｔ経過する間に劣化状態が次の状態へ遷移する確率Ｐは、λ（ｔ）を用いてＰ＝λ（ｔ）Δｔとして表される。例えば、Ｐ＝０．０１となるようΔｔを適宜設定し、時間ｔからΔｔ経過する間に劣化状態がＰ＝０．０１で次の劣化状態に推移し、劣化状態が１００に達すると破断が発生するという条件で、試験開始から破断までシミュレーションを繰り返し行うことで累積破断確率が求められる。シミュレーションにより求めた累積破断確率と、実際に同条件で実験を行った場合の結果を比較し、両者が一致すれば水素脆化進行度の予測は正しいと言える。
　図８のＨ（ｔ）の傾きから求めたλ（ｔ）を図９に示す。図９のλ（ｔ）を用いて上記のシミュレーションを１０００回繰り返し行って求めた累積破断確率と、実際に同条件で実験を行った場合の結果を図１０に示す。図１０では、シミュレーション結果を点線で示し、実験結果をプロットで示している。実験により求めた累積破断確率とλ（ｔ）を用いたシミュレーションにより求めた累積破断確率は概ね一致し、シミュレーションによって実験結果を予測できていることがわかる。したがって、環境条件が変化する場合の水素脆化進行度は正しく予測できていると言える。
　このように、環境条件が変化しない場合の水素脆化進行度をいくつか把握しておけば、環境条件が変化する場合の水素脆化進行度を予測できるようになる。実際の環境では、環境条件が変化する場合が多いことから、この方法を用いれば簡易な実験によって実際の環境での水素脆化進行度の予測が可能となる。
　（水素脆化進行度評価装置の動作）
　次に、本実施の形態における水素脆化進行度評価装置１の処理の流れについて説明する。
　図１１は、本実施の形態における水素脆化進行度評価装置１の処理の流れを示すフローチャートである。プレストレストコンクリート内の鉄筋を想定した試験条件により測定された水素脆化破断試験データは取得済みであり、記憶装置２に記憶されているとする。
　累積分布関数計算部１１は、水素脆化試験データから破断時間の累積分布関数を求める（ステップＳ１１）。
　推移確率計算部１２は、水素脆化破断を劣化状態が逐次推移して破断に至る確率過程とみなし、破断時間の累積分布関数から推移確率を計算する（ステップＳ１２）。
　水素脆化進行度計算部１３は、推移確率を時間で積分して水素脆化進行度を求める（ステップＳ１３）。
　水素脆化進行度評価装置１が更改時期推定部を備える場合、更改時期推定部は、水素脆化破断試験データの測定に使用した鉄筋の実設備における寿命データに基づいて、水素脆化進行度から実設備の更改時期を推定する。
　以上説明したように、本実施の形態によれば、プレストレストコンクリート内の鉄筋を想定した試験条件により測定された、多数の水素脆化破断試験データから破断時間の累積分布関数を計算し、水素脆化破断を劣化状態が逐次推移して破断状態に至る確率過程とみなし、破断時間の累積分布関数から推移確率を計算し、推移確率を時間で積分して水素脆化進行度を求めることで、鋼材が水素脆化によって破断に至る途中の水素脆化進行度を知ることができる。これによって、事前に鋼材の水素脆化進行度を評価しておけば、鋼材の使用期間中にどの程度水素脆化が進行するのかを予想することができ、鋼材を使用した設備の更改の判断材料とすることができる。
　上記説明した水素脆化進行度評価装置１には、例えば、図１２に示すような、中央演算処理装置（ＣＰＵ）９０１と、メモリ９０２と、ストレージ９０３と、通信装置９０４と、入力装置９０５と、出力装置９０６とを備える汎用的なコンピュータシステムを用いることができる。このコンピュータシステムにおいて、ＣＰＵ９０１がメモリ９０２上にロードされた所定のプログラムを実行することにより、水素脆化進行度評価装置１が実現される。
　１…水素脆化進行度評価装置
　１１…累積分布関数計算部
　１２…推移確率計算部
　１３…水素脆化進行度計算部
　２…記憶装置
　３…表示装置

Claims (10)

  1. 　プレストレストコンクリート内の鉄筋を想定した試験条件により測定された水素脆化破断試験データから破断時間の累積分布関数を求める第１工程と、
    　水素脆化破断を劣化状態が逐次推移して破断に至る確率過程とみなし、前記破断時間の累積分布関数から推移確率を計算する第２工程と、
    　前記推移確率を時間で積分して水素脆化進行度を求める第３工程と、
    　を有することを特徴とする水素脆化進行度の評価方法。
  2. 　前記第２工程では、前記推移確率をλ（ｔ）＝Ａｔ－Ｂ＋Ｃと近似して求めることを特徴とする請求項１に記載の水素脆化進行度の評価方法。
  3. 　前記第２工程では、前記水素脆化進行度のそれぞれに対応する時間の値を求め、前記水素脆化進行度の変化の割合を前記推移確率とすることを特徴とする請求項１に記載の水素脆化進行度の評価方法。
  4. 　前記第３工程では、さらに、異なる環境条件に基づいて求めた水素脆化進行度を組み合わせて環境条件が変化する場合の水素脆化進行度を求めることを特徴とする１ないし３のいずれかに記載の水素脆化進行度の評価方法。
  5. 　前記水素脆化破断試験データの測定に使用した鉄筋の実設備における寿命データに基づいて、前記水素脆化進行度から実設備の更改時期を推定する第４工程を有することを特徴とする請求項１ないし４のいずれかに記載の水素脆化進行度の評価方法。
  6. 　プレストレストコンクリート内の鉄筋を想定した試験条件により測定された水素脆化破断試験データから破断時間の累積分布関数を求める累積分布関数計算部と、
    　水素脆化破断を劣化状態が逐次推移して破断に至る確率過程とみなし、前記破断時間の累積分布関数から推移確率を計算する推移確率計算部と、
    　前記推移確率を時間で積分して水素脆化進行度を求める水素脆化進行度計算部と、
    　を有することを特徴とする水素脆化進行度評価装置。
  7. 　前記推移確率計算部は、前記推移確率をλ（ｔ）＝Ａｔ－Ｂ＋Ｃと近似して求めることを特徴とする請求項６に記載の水素脆化進行度評価装置。
  8. 　前記推移確率計算部は、前記水素脆化進行度のそれぞれに対応する時間の値を求め、前記水素脆化進行度の変化の割合を前記推移確率とすることを特徴とする請求項６に記載の水素脆化進行度評価装置。
  9. 　前記水素脆化進行度計算部は、さらに、異なる環境条件に基づいて求めた水素脆化進行度を組み合わせて環境条件が変化する場合の水素脆化進行度を求めることを特徴とする６ないし８のいずれかに記載の水素脆化進行度評価装置。
  10. 　前記水素脆化破断試験データの測定に使用した鉄筋の実設備における寿命データに基づいて、前記水素脆化進行度から実設備の更改時期を推定する更改時期推定部を有することを特徴とする請求項６ないし９のいずれかに記載の水素脆化進行度評価装置。
PCT/JP2019/037857 2018-10-01 2019-09-26 水素脆化進行度の評価方法及び水素脆化進行度評価装置 Ceased WO2020071237A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2020550356A JP7108214B2 (ja) 2018-10-01 2019-09-26 水素脆化進行度の評価方法及び水素脆化進行度評価装置
US17/275,144 US20220042901A1 (en) 2018-10-01 2019-09-26 Hydrogen Embrittlement Progress Evaluation Method and Hydrogen Embrittlement Progress Evaluation Device

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2018-186637 2018-10-01
JP2018186637 2018-10-01

Publications (1)

Publication Number Publication Date
WO2020071237A1 true WO2020071237A1 (ja) 2020-04-09

Family

ID=70054622

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/037857 Ceased WO2020071237A1 (ja) 2018-10-01 2019-09-26 水素脆化進行度の評価方法及び水素脆化進行度評価装置

Country Status (3)

Country Link
US (1) US20220042901A1 (ja)
JP (1) JP7108214B2 (ja)
WO (1) WO2020071237A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022097205A1 (ja) * 2020-11-04 2022-05-12 日本電信電話株式会社 破断確率推定方法

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020261422A1 (ja) * 2019-06-26 2020-12-30 日本電信電話株式会社 鉄筋の水素脆化破断リスクの評価方法
CN116296798B (zh) * 2023-02-15 2024-10-08 中国铁道科学研究院集团有限公司铁道建筑研究所 一种基于dic技术测定受拉混凝土本构关系的方法及系统
KR20240149991A (ko) 2023-04-07 2024-10-15 주식회사 피레타 듀얼 감도 방식의 ae 신호감지를 이용한 수소 이송관 안정성 측정장치
CN121024910B (zh) * 2025-10-31 2026-01-30 厚普智慧物联科技有限公司 一种适用于氢能压缩机的状态监测及故障预测系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004340817A (ja) * 2003-05-16 2004-12-02 Mitsubishi Heavy Ind Ltd 金属の水素脆化予測法及び水素脆化予測法を用いた金属部材設計方法
US20120024077A1 (en) * 2010-07-29 2012-02-02 Kyushu University, National University Corporation Inclusion rating method
JP2017219532A (ja) * 2016-06-03 2017-12-14 日本電信電話株式会社 水素脆化特性評価試験における試験片面積の算出方法および試験セルサイズの算出方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002365202A (ja) * 2001-06-07 2002-12-18 Kawasaki Heavy Ind Ltd ボルトの水素脆化割れ予測方法
JP2004309197A (ja) 2003-04-03 2004-11-04 Sumitomo Metal Ind Ltd 耐遅れ破壊性の評価方法
JP5986126B2 (ja) * 2014-03-05 2016-09-06 日本電信電話株式会社 水素脆化評価方法
JP6225087B2 (ja) 2014-08-26 2017-11-01 日本電信電話株式会社 材料評価のための前処理方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004340817A (ja) * 2003-05-16 2004-12-02 Mitsubishi Heavy Ind Ltd 金属の水素脆化予測法及び水素脆化予測法を用いた金属部材設計方法
US20120024077A1 (en) * 2010-07-29 2012-02-02 Kyushu University, National University Corporation Inclusion rating method
JP2017219532A (ja) * 2016-06-03 2017-12-14 日本電信電話株式会社 水素脆化特性評価試験における試験片面積の算出方法および試験セルサイズの算出方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
IZUMI, MASANORI ET AL.: "Transactions of the Architectural Institute of Japan", CONCRETE FRACTURE THEORY BASED ON STOCHASTIC PROCESS THEORY, vol. 287, 30 January 1980 (1980-01-30), pages 1 - 13, Retrieved from the Internet <URL:https://jstage.jst.go.jp/article/aijsaxx/287/0/287_KJ00003748923/_article/-char/ja> [retrieved on 20191028] *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022097205A1 (ja) * 2020-11-04 2022-05-12 日本電信電話株式会社 破断確率推定方法
JPWO2022097314A1 (ja) * 2020-11-04 2022-05-12
WO2022097314A1 (ja) * 2020-11-04 2022-05-12 日本電信電話株式会社 破断確率推定方法
JP7464888B2 (ja) 2020-11-04 2024-04-10 日本電信電話株式会社 破断確率推定方法

Also Published As

Publication number Publication date
JP7108214B2 (ja) 2022-07-28
JPWO2020071237A1 (ja) 2021-09-02
US20220042901A1 (en) 2022-02-10

Similar Documents

Publication Publication Date Title
JP7108214B2 (ja) 水素脆化進行度の評価方法及び水素脆化進行度評価装置
Lin et al. Integrating random shocks into multi-state physics models of degradation processes for component reliability assessment
CN110197288B (zh) 故障影响下设备的剩余使用寿命预测方法
US20090271235A1 (en) Apparatus and method for generating survival curve used to calculate failure probability
JP6352851B2 (ja) 水素脆化感受性評価方法
CN102980806A (zh) 一种预测多步加载条件下金属材料低周疲劳寿命的方法
KR20150057553A (ko) 피로수명 예측방법
CN109117566A (zh) 一种基于产品寿命预测模型的混合维修规划方法
CN109522517B (zh) 一种电表失效率计算方法和系统
JP5064259B2 (ja) 劣化評価装置及び劣化評価方法
JP6470716B2 (ja) 水素脆化特性評価試験における試験片面積の算出方法および試験セルサイズの算出方法
Chen et al. Stochastic modelling and lifecycle performance assessment of bond strength of corroded reinforcement in concrete
JP7029061B2 (ja) 亀裂等発生検知方法、亀裂等発生検知装置及び亀裂等発生検知プログラム
JP2003004599A (ja) 構造物の疲労損傷度の評価方法及び表示方法
JP6940774B2 (ja) 水素侵入挙動推定方法、水素侵入挙動推定装置及び水素侵入挙動推定プログラム
Bayane et al. Quantification of the conditional value of SHM data for the fatigue safety evaluation of a road viaduct
JP7464888B2 (ja) 破断確率推定方法
US20190250890A1 (en) Testing of physical random number generators
CN117725664A (zh) 开孔板连接件的极限承载力计算方法、装置、设备及介质
CN119670882B (zh) 一种桥梁时变可靠度评估方法、装置、设备及介质
CN109992875B (zh) 一种切换设备剩余寿命的确定方法及系统
US11946856B2 (en) Method for estimating steel rupture starting point, device for estimating steel rupture starting point, and program for estimating steel rupture starting point
KR20180078112A (ko) 저장 장치의 수명 예측 장치 및 방법
Al-Osta et al. Probability of failure of corroding reinforced concrete columns under eccentric loading
Nepal et al. Gamma process modelling for lifecycle performance assessment of corrosion affected concrete structures

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19868799

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2020550356

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19868799

Country of ref document: EP

Kind code of ref document: A1