ニッケル基合金部品の寿命推定方法
【課題】多結晶のニッケル基合金からなるニッケル基合金部品について高精度に寿命を推定することができるニッケル基合金部品の寿命推定方法を提供する。
【解決手段】多結晶ニッケル基合金から形成されたニッケル基合金部品１の寿命推定方法であって、ニッケル基合金部品１の各部位における使用時の使用温度の推定値と応力の推定値とに基づいて、相対的に損傷の度合いが大きい最損傷部位Ｗを特定する第１工程と、最損傷部位Ｗについて１または複数の手法により破断までの破断時間を求める第２工程と、を有することを特徴とする。
【解決手段】多結晶ニッケル基合金から形成されたニッケル基合金部品１の寿命推定方法であって、ニッケル基合金部品１の各部位における使用時の使用温度の推定値と応力の推定値とに基づいて、相対的に損傷の度合いが大きい最損傷部位Ｗを特定する第１工程と、最損傷部位Ｗについて１または複数の手法により破断までの破断時間を求める第２工程と、を有することを特徴とする。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ニッケル基合金部品の寿命推定方法に関するものである。
【背景技術】
【０００２】
ニッケル基合金は、高温下における強度と耐クリープ性に優れ、ガスタービン、ジェットエンジン、化学プラント等に使用されている。このような性質を有するニッケル基合金であっても、高温・高応力に曝される使用環境下では、経年劣化による疲労破壊やクリープ破壊等が生じる。従って、このような高温・高応力下で使用されるニッケル基合金部品の余寿命を正確且つ定量的に予測することは、ニッケル基合金部品の検査や交換等の時期を計画する上で非常に重要である。
【０００３】
例えば下記特許文献１には、クリープ試験終了後のニッケル基合金試料から、（１）γ´相（析出物）の円相当直径の平均値とラーソン・ミラー・パラメータとの関係、（２）ラーソン・ミラー・パラメータとラフト化比率（γ´相の応力負荷を受ける方向のアスペクト比と無負荷のアスペクト比の比）との関係を予め求めておき、使用時間が既知の寿命推定対象におけるγ´相の円相当直径及び（１）の関係からラーソン・ミラー・パラメータを求め、このラーソン・ミラー・パラメータ、寿命推定対象を測定して求めたラフト化比率及び（２）の関係に基づいてクリープ破断寿命比（使用時間／未使用状態からの破断時間）を推定する技術が開示されている。
【特許文献１】特開平８−２９７０７９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ところで、例えば論文「三浦信祐、他２名、“多結晶Ｎｉ基超合金，ＩＮ−１００，のクリープ変形により形成されるγ´相のラフト構造”、鉄と鋼、社団法人日本鉄鋼協会、２００４年８月１日、Ｖｏｌ．９０、Ｎｏ．８、ｐ．５６７−５７３」に記載されているように、多結晶ニッケル基合金の結晶方位は各結晶粒によって異なり、またγ´相のラフト化は結晶粒の結晶方位によって異なる。
従って、上記従来技術は、単結晶ニッケル基合金には適用できるものの、多結晶ニッケル基合金には適用することができない。すなわち、多結晶ニッケル基合金では各結晶粒の結晶方位が異なるので、上記従来技術のように各結晶粒の結晶方位の違いを考慮しないで一律に測定した応力負荷面と応力負荷無し面の長径及び短径に基づいてラフト化比率を求める手法は、多結晶ニッケル基合金に適用した場合に大きな誤差を含むものとなる。
【０００５】
本発明は上述した事情に鑑みてなされたものであり、多結晶ニッケル基合金からなるニッケル基合金部品について迅速かつ高精度に寿命を推定することができるニッケル基合金部品の寿命推定方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明は、上記課題を解決するために以下の手段を採用する。
すなわち、本発明は、第１の解決手段として、多結晶ニッケル基合金から形成されたニッケル基合金部品の寿命推定方法であって、前記ニッケル基合金部品の各部位における使用時の使用温度の推定値と応力の推定値とに基づいて、相対的に損傷の度合いが大きい最損傷部位を特定する第１工程と、前記最損傷部位について１または複数の手法により破断までの破断時間を求める第２工程とを有する、という手段を採用する。
【０００７】
また、本発明は、第２の解決手段として、上記第１の解決手段における前記第１工程が、前記ニッケル基合金部品の使用時の応力軸方向に垂直な断面において前記使用温度の推定値が最も高い垂直断面最損傷部位を特定する第１Ａ工程と、前記垂直断面最損傷部位の応力軸方向に位置する前記各部位のうち、前記使用温度の推定値と、前記応力の推定値と、前記ニッケル基合金部品の使用開始から経過した使用時間と、前記ニッケル基合金部品の同等部品について予め取得されたクリープマスターカーブとに基づいて、前記最損傷部位を特定する第１Ｂ工程と、からなるという手段を採用する。
【０００８】
また、本発明は、第３の解決手段として、上記第１又は第２の解決手段における前記ニッケル基合金部品が前記ニッケル基合金からなる基材の表面に金属コーティング層が形成されたものである場合、前記第１工程では、前記ニッケル基合金部品と同等な複数の試料について、前記ニッケル基合金部品の使用によって前記基材に形成される拡散層の厚さ（拡散層厚）とラーソン・ミラー・パラメータとの関係をそれぞれ求める第１Ｃ工程と、各試料に関する前記関係に基づいて拡散層厚とラーソン・ミラー・パラメータとの関係を示す第１の近似式を求める第１Ｄ工程と、前記第１の近似式とラーソン・ミラー・パラメータの定義式との第１の連立方程式に、前記各部位における拡散層厚の測定値と、前記使用時間とを代入する第１Ｅ工程と、から前記使用温度の推定値を求める、という手段を採用する。
また、本発明は、第４の解決手段として、上記第３の解決手段における前記第２工程では、前記第１Ｃ工程から第１Ｅ工程によって求められた前記使用温度の推定値を、前記クリープマスターカーブに当てはめることにより前記最損傷部位の破断時間を求める、という手段を採用する。
【０００９】
また、本発明は、第５の解決手段として、上記第１から４のうちいずれかの解決手段における前記第２工程が、異なる試験温度でクリープ試験を行った前記ニッケル基合金部品と同等な複数の試料について、前記試料の応力軸に直交する直交面における複数の結晶粒の結晶方位と結晶粒の結晶格子面に沿って析出した複数の析出物の厚さとをそれぞれ計測すると共にラーソン・ミラー・パラメータを求める第２Ａ工程と、前記結晶粒の結晶方位と析出物の厚さとに基づいて析出物の真の厚さ（析出物厚）を求める第２Ｂ工程と、前記複数の試料のラーソン・ミラー・パラメータと、前記結晶方位と前記応力軸とのなす負荷角度の余弦値との関係の直線性が良好な領域における特定の前記負荷角度の前記析出物厚と、に基づいて第２の近似式を求める第２Ｃ工程と、前記第２の近似式とラーソン・ミラー・パラメータの定義式との第２の連立方程式に、前記最損傷部位から得られる前記特定の負荷角度と同一の負荷角度における析出物厚と、前記ニッケル基合金部品の使用開始から経過した使用時間と、を代入することにより前記最損傷部位における使用温度の推定値を求める第２Ｄ工程と、このようにして求めた前記最損傷部位の使用温度の推定値を、前記クリープマスターカーブに当てはめることにより前記最損傷部位における破断時間を求める第２Ｅ工程と、に基づいて前記最損傷部位の破断時間を求める、という手段を採用する。
【００１０】
また、本発明は、第６の解決手段として、上記第１から５の解決手段における前記第２工程では、前記最損傷部位から作成した試料に対してクリープ試験を行って破断時間を求める、という手段を採用する。
また、本発明は、第７の解決手段として、上記第１から６の解決手段における前記ニッケル基合金部品がガスタービンの動翼または静翼である、という手段を採用する。
【発明の効果】
【００１１】
本発明によれば、ニッケル基合金部品の各部位のうち相対的に損傷の度合いが大きい最損傷部位を特定し、最損傷部位について１または複数の手法により破断時間（寿命）の推定を行うので、最損傷部位以外の部位における破断時間を推定する必要が無く、クリープ破壊の生じる可能性が最も高い部位に集中して１または複数の手法で寿命の推定を行うことができ、迅速且つ高精度に寿命を推定することが可能となる。
【００１２】
また、前記第１工程は、垂直断面最損傷部位の応力軸方向に位置する各部位のうち、使用温度の推定値と、応力の推定値と、ニッケル基合金部品の使用開始から経過した使用時間と、前記ニッケル基合金部品の同等部品について予め取得されたクリープマスターカーブとに基づいて最損傷部位を特定するので、最損傷部位を迅速に特定することができ、より迅速且つ高精度にニッケル基合金部品の寿命の推定を行うことができる。
【発明を実施するための最良の形態】
【００１３】
以下、本発明の一実施形態について、図面を参照して説明する。本実施形態では、ガスタービンにおける圧縮機やタービンの構成部品となる動翼（ガスタービン動翼）を寿命推定対象とする。
【００１４】
図１は、使用開始から所定時間が経過したガスタービン動翼１を示す図であり、図２は、使用開始から所定時間が経過したガスタービン動翼１の要部断面図である。ガスタービン動翼１は、図１に示すように翼型断面を有する棒状部品であり、タービン軸（回転軸）に放射状に取り付けられて使用される。また、ガスタービン動翼１は、図２に示すように、多結晶ニッケル基合金からなる基材１ａ上に拡散層１ｂ、コーティング層１ｃ及び酸化膜層１ｄが形成されている。ガスタービン動翼１はタービン回転軸に放射状に取り付けられて使用に供されるので、応力軸Ｌは、図示するようにガスタービン動翼１の長さ方向となる。
【００１５】
本来（使用開始前）のガスタービン動翼１は、基材１ａの表面にアルミニウム（Ａｌ）からなるコーティング層１ｃが形成されたものであるが、高温環境下での使用によってコーティング層１ｃのアルミニウムが基材１ａ内に拡散することによって拡散層１ｂが形成されると共に、酸化によってコーティング層１ｃの表面に酸化膜層１ｄが形成される。
【００１６】
本実施形態は、ガスタービン動翼１の全体を評価対象とするのではなく、ガスタービン動翼１についてクリープ損傷が最も激しい部位（最損傷部位）を最初に特定し、当該最損傷部位について寿命の推定を行うものである。
より具体的には、ガスタービン動翼１に作用する遠心力Ｆの方向（応力軸Ｌの方向）に垂直な断面Ｓｒにおけるガスタービン動翼１の表面（翼状表面）の使用時における温度分布（使用温度の分布）を推定し、この断面Ｓｒにおいて使用温度が最も高い部位を垂直断面最損傷部位Ｂ１として特定する。すなわち、断面Ｓｒは応力軸Ｌに垂直なので、断面Ｓｒにおけるガスタービン動翼１の表面の各部位には同一の応力が作用すると考えることができ、よって上記表面の各部位のうち使用温度が最も高い部位Ｂ１は、断面Ｓｒにおいてクリープ損傷が最も激しい部位（垂直断面最損傷部位Ｂ１）とみなすことができる。
【００１７】
次に、垂直断面最損傷部位Ｂ１の応力軸方向におけるガスタービン動翼１の表面についての使用温度の分布と使用時応力の分布とをそれぞれ推定する。すなわち、遠心力Ｆの方向に垂直な各断面におけるガスタービン動翼１の表面の使用温度の分布の状態はほぼ同様のものとなるので、垂直断面最損傷部位Ｂ１の遠心力方向に位置する各部位が各断面において最も使用温度が高い部位であるとみなすことができる。また、上述の通り各断面において生じる応力はその断面において同一と考えられるから、垂直断面最損傷部位Ｂ１の遠心力Ｆの方向に位置する各部位が各断面においてクリープ損傷が最も激しい部位であるとみなすことができる。従って、この垂直断面最損傷部位Ｂ１の遠心力方向に位置する各部位の中で最もクリープ損傷が激しい部位が最損傷部位Ｗとなる。
【００１８】
そして、垂直断面最損傷部位Ｂ１の遠心力方向に位置する各部位における使用時応力（推定値）から、予めガスタービン動翼１と同等の試料から取得したクリープ破断データを参照してラーション・ミラー・パラメータＰを求めて、このラーション・ミラー・パラメータＰと、垂直断面最損傷部位Ｂ１の遠心力方向に位置する各部位における使用温度（推定値）とを下記ラーソン・ミラー・パラメータ定義式（１）に代入することにより各破断時間を算出し、この破断時間が最も短い部位をガスタービン動翼１の全体において損傷が最も激しい部位Ｗ（最損傷部位）と特定する。そして、本実施形態では、この最損傷部位Ｗのみについて破断時間（余寿命）の推定を行う。
【００１９】
まず、上述のように、初めにガスタービン動翼１に作用する遠心力方向に垂直な断面Ｓｒの温度分布を求める。本実施形態においては、基材１ａ上に形成される拡散層１ｂの厚さｄ（拡散層厚）から断面Ｓｒの温度分布を求める。
具体的には、以下に説明するように、ラーソン・ミラー・パラメータＰの定義式（１）と、使用開始前のガスタービン動翼１と同等な複数の試料に関する熱時効試験（熱処理）の試験結果から同定されたラーソン・ミラー・パラメータＰＥと拡散層１ｂの厚さｄＥ（拡散層厚）との近似式（２）と、Ｐ＝ＰＥとして上記定義式（１）と近似式（２）とを連立させて得られる関係式（３）に基づいてガスタービン動翼１の使用温度Ｔの推定値を求め、この使用温度Ｔの推定値をガスタービン動翼１の使用温度と初期破断時間とに相当する関係（クリープマスターカーブ）に当てはめることにより破断時間（寿命）を推定するものである。なお、各変数の添字Ｅは、試料に関するものを意味するが、関係式（３）においてはガスタービン動翼１の拡散層厚ｄを代入するので、ｄＥをｄとしている。
Ｐ＝Ｔ・（Ｃ＋ｌｏｇｔ） （１）
ＰＥ＝（ｄＥ−β１）／α１ （２）
Ｔ＝（ｄ−β１）／{α１（Ｃ＋ｌｏｇｔ）} （３）
【００２０】
ラーソン・ミラー・パラメータは、周知のように、金属部品のクリープ破壊における温度と時間との関係を示すパラメータとして、上記定義式（１）で表される。この定義式（１）において、Ｔは使用温度（絶対温度）、ｔは使用時間、Ｃは材料固有の材料定数である。
【００２１】
また、近似式（２）におけるα１，β１は、上記試料に関する実験結果に基づいて同定される定数である。本実施形態では、以下に説明するように、定数α１，β１をより正確に同定するために、実際に使用された複数のガスタービン動翼１から取得した拡散層厚ｄの分布範囲に該当する熱時効試験（熱処理）の試験結果に限定して材料定数Ｃを最適化する。
【００２２】
上記熱時効試験（熱処理）では、上記複数の試料について試験温度ＴＥと試験時間ｔＥを変えた熱時効試験（熱処理）を行う。試験温度ＴＥは、大まかに把握されるガスタービン動翼１の実際の使用温度範囲内及びその近傍における複数温度、例えば４つの試験温度Ｔ１，Ｔ２，Ｔ３，Ｔ４が設定され、一方、試験時間ｔＥについては、ガスタービン動翼１の実際の使用時間範囲内における複数の時間、例えば３つの試験時間ｔ１，ｔ２，ｔ３が設定される。すなわち、熱時効試験（熱処理）では、各試料について４つの試験温度Ｔ１，Ｔ２，Ｔ３，Ｔ４と３つの試験時間ｔ１，ｔ２，ｔ３とが設定される。
【００２３】
そして、上記定義式（１）における使用温度Ｔに試験温度Ｔ１，Ｔ２，Ｔ３，Ｔ４をそれぞれ代入すると共に使用時間ｔに試験時間ｔ１，ｔ２，ｔ３をそれぞれ代入することにより、各試料に関するラーソン・ミラー・パラメータＰＥを求める。ラーソン・ミラー・パラメータＰＥの算出に際しては、別途クリープ試験等で材料定数Ｃを算出してこの値を設定する。
【００２４】
続いて、上記各試料を研磨して内部を露出させ、顕微鏡を用い拡散層厚ｄＥを計測する。この拡散層厚ｄＥは、顕微鏡から出力される試料の画像について画像処理装置を用いて画像処理することによって、あるいは検査員が試料の画像を目視することによって計測される。なお、顕微鏡は拡散層厚ｄＥを計測するのに適切な分解能を有していればよく、例えば光学顕微鏡や走査型電子顕微鏡を用いることができる。
【００２５】
図３は、このようにして求められた各試料のラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥとの関係を示す点（計測点）を片対数グラフ上に示した図である。この図では、試料毎に異なる表記方法で計測点を示している。本実施形態では、このような各試料に関する計測点の全てに例えば最小二乗法を適用することにより、全ての試料に関するラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥとの関係を示す特性直線を求める。このようにして各試料に関する熱時効試験（熱処理）の試験結果から得られた特性直線は、上述した近似式（２）において定数α１、β１が同定されたものである。
【００２６】
ここで、図３に破線で示すように、同一温度における試料の計測点は良好な直線性を示す。つまり各温度については近似性の高い特性直線を求めることができるが、全ての温度の計測点を１本の特性直線を求めるための要素とした場合には、異なる温度域が含まれるので、近似性の高い１本の特性直線を求めることは困難である。
【００２７】
本実施形態では、このような事情を回避して、より近似性の高い特性直線Ｎ、つまり近似式（２）を求めるために、上記と同様に顕微鏡を用いて実際に使用された複数のガスタービン動翼１から拡散層厚ｄの分布範囲Ｈを予め取得し、各試料に関するラーソン・ミラー・パラメータＰＥを求める際に、分布範囲Ｈに該当する試料の熱時効試験（熱処理）の試験結果（計測点）のみに限定して材料定数Ｃを最適化する。すなわち、近似式（２）の近似性が最も高くなるように、図３に示すように各試料に関する計測点を分布範囲Ｈ内に狭く限定し、この上で各試料に関するラーソン・ミラー・パラメータＰＥを求める際の材料定数Ｃを最適化する。なお、関係式（３）における材料定数Ｃも上記最適化したものを用いる。
【００２８】
図４は、このような材料定数Ｃを最適化処理前の各試料のラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥとの関係を示す計測点を片対数グラフ上に示した図である。この図４に示すように、材料定数Ｃについて最適化処理を行わない場合、各試料の計測点に関する特性直線は、分布範囲Ｈ内において比較的乖離しており、各試料の全計測点から求められる１本の特性直線は近似性が低いものとなる。
【００２９】
本実施形態では、近似式（２）を特定することによって得られた近似性の高い関係式（３）に、ガスタービン動翼１に作用する遠心力方向に垂直な断面Ｓｒにおける各周縁の拡散層厚ｄを上記と同様に顕微鏡観察により求め、この拡散層厚ｄを関係式（３）に代入することにより、断面Ｓｒにおけるガスタービン動翼１の表面（翼状表面）の使用時における温度分布（使用温度Ｔの分布）を求める。
例えば、図１において、断面Ｓｒの縁に沿って拡散層厚ｄを計測していき、この計測された拡散層厚ｄを順次関係式（３）に代入して、断面Ｓｒにおけるガスタービン動翼１の表面（翼状表面）の使用時における温度分布（使用温度Ｔの分布）を推定し、この断面Ｓｒにおいて使用温度が最も高い部位を垂直断面最損傷部位Ｂ１として特定する。
【００３０】
次に、垂直断面最損傷部位Ｂ１を特定した後に、ガスタービン動翼１の垂直断面最損傷部位Ｂ１を含む垂直断面最損傷部位Ｂ１の遠心力方向の各部位における使用温度Ｔの推定値と応力σの推定値を求める。図１においては、遠心力方向に沿う全ての部位のうち部位Ｂ１１〜Ｂ１４を例示として示すが、遠心力方向に沿う全ての部位について使用温度Ｔを推定して構わない。ここで、使用温度Ｔの推定値については、断面Ｓｒの場合と同様に拡散層厚ｄと近似式（３）から、応力σについては、その部位に生じる遠心力Ｆ及びその部位から遠心力方向側端面までの棒の自重並びに断面Ｓｒの断面積に基づいて求めている。なお、遠心力方向に垂直な各断面の断面積は一定としている。
【００３１】
図５は、上記ガスタービン動翼１と同等な動翼にクリープ試験を行って得られたクリープマスターカーブを示しており、試料に負荷した応力σとラーソン・ミラー・パラメータＰとの関係に相当する特性直線を示している。なお、このクリープマスターカーブのラーソン・ミラー・パラメータＰは、未使用状態から破断までの使用時間ｔ（初期破断時間）と使用温度Ｔとから求めている。
【００３２】
この図５に示すクリープマスターカーブに、上記のようにして求めた各部位における応力σの推定値を当てはめることにより、各部位におけるラーソン・ミラー・パラメータＰを求め、この各ラーソン・ミラー・パラメータＰと各部位の使用温度Ｔとラーソン・ミラー・パラメータの定義式（１）より未使用状態からの初期破断時間を求める。さらに、初期破断時間から使用時間ｔを減算することにより、各部位の破断までの破断時間（余寿命）を算出する。そして、各部位の中で破断時間（余寿命）が最も短い部位を最損傷部位Ｗとして特定する。
例えば関係式（３）により、ある部位Ｂ１１における応力σの推定値σ１１が求まれば、クリープマスターカーブにより、応力σ１１に対応するラーソン・ミラー・パラメータＰ１１を求めることができる。このラーソン・ミラー・パラメータＰ１１と、部位Ｂ１１における使用温度Ｔの推定値と、ラーソン・ミラー・パラメータの定義式（１）に代入することにより、部位Ｂ１１における初期破断時間を求めることができ、この初期破断時間の推定値から使用時間を減算することにより、破断時間（余寿命）を推定することができる。
【００３３】
上記のように最損傷部位Ｗを特定した後、拡散層１ｂからの寿命推定手法、析出物評価からの寿命推定手法、クリープ試験からの寿命推定手法のうちいずれかを行う。但し、破断時間（余寿命）の推定結果は、拡散層１ｂからの寿命推定手法、析出物評価からの寿命推定手法、クリープ試験に基づく寿命推定手法の順に早く求めることができるが、逆にその結果の信頼性は順に低いものとなっている。そのため、ガスタービン動翼１の作動環境や効率、使用時間等によっていずれかの上記寿命推定手法が適宜選択されて破断時間（余寿命）が求められる。
【００３４】
図６は、上記ガスタービン動翼１と同等な動翼にクリープ試験を行って得られたクリープマスターカーブを示しており、上記ガスタービン動翼１の未使用状態から破断までの使用時間ｔ（初期破断時間）と使用温度Ｔとの関係に相当する特性曲線を示している。
【００３５】
以下、各寿命推定手法について説明する。
まず、初めに拡散層１ｂからの寿命推定手法について説明する。
拡散層１ｂからの寿命推定手法は、上記のように求めた最損傷部位Ｗにおける拡散層厚ｄを測定値から求めた使用温度Ｔの推定値をクリープマスターカーブに当てはめることによってガスタービン動翼１の未使用状態から破断までの使用時間ｔ（初期破断時間）の推定値を求める。そして、ガスタービン動翼１の使用開始から現在までの使用時間ｔを上記初期破断時間の推定値から減算することによりガスタービン動翼１の破断時間（余寿命）を求める。例えば関係式（３）により、最損傷部位Ｗの使用温度ＴＳが求まると、クリープマスターカーブにより使用温度ＴＳにおける初期破断時間ｔＯの推定値が求められる。そして、この初期破断時間ｔＯの推定値からガスタービン動翼１の使用時間ｔＳを減算すれば、最損傷部位Ｗの破断時間（余寿命）ｔＮが推定される。
【００３６】
続いて、析出物評価からの寿命推定手法について説明する。
図７は、使用開始から所定時間が経過したガスタービン動翼１の上記応力軸Ｌと直行する断面Ｓｒにおける結晶粒と結晶方位の関係を示す測定結果である。この図７に示すように、ガスタービン動翼１の応力軸Ｌに直交する断面Ｓｒには、多数の結晶粒Ｇが分布している。また、これら結晶粒Ｇは、図示する立方体が示すように互いに異なる結晶方位を有している。
なお、ガスタービン動翼１は、使用開始から時間が経過するに従って、未使用時には立方体状であった結晶粒内の析出物（γ´相）が板状となって、結晶粒の結晶格子面に沿って粗大化することが知られている。
【００３７】
以下に説明するように、ラーソン・ミラー・パラメータＰの定義式（１）と、使用開始前のガスタービン動翼１と同等な複数の試料に関するクリープ試験の試験結果から同定されたラーソン・ミラー・パラメータＰＥと析出物の真の厚さ（析出物厚ＤＥ）との近似式（４）と、Ｐ＝ＰＥとして上記定義式（１）と近似式（４）とを連立させて得られる関係式（５）に基づいてガスタービン動翼１の使用温度Ｔの推定値を求め、この使用温度Ｔの推定値をガスタービン動翼１の使用温度Ｔと初期破断時間との関係を示す特性図（図６に示すクリープマスターカーブ）に当てはめることにより破断時間（余寿命）を推定するものである。
Ｐ＝Ｔ・（Ｃ＋ｌｏｇｔ） （１）
ＰＥ＝（ＤＥ−β２）／α２ （４）
Ｔ＝（Ｄ−β２）／{α２（Ｃ＋ｌｏｇｔ）} （５）
【００３８】
また、本実施形態では、上記近似式（４）を求めるに際して、その近似性を高めるために、各々の結晶粒の結晶方位を考慮することによって、正確な析出物厚ＤＥを求めると共に、析出物厚ＤＥと結晶方位の関係から定数α２，β２をより正確に同定する。なお、各変数の添字Ｅは、上記と同様に試料に関するものを意味するが、関係式（５）においてはガスタービン動翼１の析出物厚Ｄを代入するので、ＤＥをＤと表記している。
【００３９】
上記クリープ試験では、上記複数の試料について、同一の引張（負荷）応力σＥについて異なる試験温度ＴＥが設定される。試験温度ＴＥは、ガスタービン動翼１の実際の使用温度範囲内における複数温度、例えば３つの試験温度Ｔ４，Ｔ５，Ｔ６が設定され、一方、引張応力については、例えばガスタービン動翼１に作用する実際の負荷応力、例えばガスタービン動翼１に加わる遠心力Ｆを断面Ｓｒの断面積で除した値が設定される。
【００４０】
そして、上記定義式（１）における使用温度Ｔに試験温度Ｔ４，Ｔ５，Ｔ６（Ｔ４＜Ｔ５＜Ｔ６）をそれぞれ代入すると共に、使用時間ｔに各試料の試験時間ｔ４，ｔ５，ｔ６（ｔ４＜ｔ５＜ｔ６）をそれぞれ代入することにより、各試料に関するラーソン・ミラー・パラメータＰＥを求める。
【００４１】
図８は、上記試料における析出物の測定厚ＭＥと析出物厚ＤＥとの関係を示す模式図であり、符号ＧＥは結晶粒を、またＥＥは析出物をそれぞれ示している。上述したように各試料についてラーソン・ミラー・パラメータＰＥを求めて、各試料に作用する応力の応力軸ＬＥに直交する断面Ｓａを露出させて、走査型電子顕微鏡を用いることにより当該断面Ｓａの各析出物ＥＥについて測定厚ＭＥをそれぞれ測定すると共に、走査型電子顕微鏡にＥＢＳＰ（Electron Back Scattering Pattern）検出器を接続することにより上記断面Ｓａの各結晶粒ＧＥについて結晶方位ＤＩＥ（結晶格子面Ｓｂの法線方向）をそれぞれ測定する。
【００４２】
図９は、上記図８における結晶粒ＧＥの結晶方位ＤＩＥ（結晶格子面Ｓｂの法線方向）と応力軸ＬＥとのなす角度θＥ（負荷角度）を示す模式図である。
各結晶粒ＧＥの断面Ｓａは、応力軸ＬＥに直交する面であるが、結晶粒ＧＥの結晶方位ＤＩＥは、応力軸ＬＥに関係なく個々の方向となる。すなわち、負荷角度θＥは、結晶粒ＧＥ毎に異なった値となる。従って、各結晶粒ＧＥについて得られた測定厚ＭＥに同じく各結晶粒ＧＥについて得られた負荷角度θＥの余弦値cosθＥ（負荷角余弦値）を乗算することにより、各結晶粒ＧＥについて結晶方位ＤＩＥの方向における析出物ＥＥの真の厚さ（析出物厚ＤＥ）を求める。
【００４３】
以上の処理によって、試料毎にラーソン・ミラー・パラメータＰＥが求められると共に、各試料の断面Ｓａにおける複数の結晶粒ＧＥについて結晶方位ＤＩＥが求められ、また複数の結晶粒ＧＥの結晶格子面Ｓｂに沿って各々粗大化する複数の析出物ＥＥについて析出物厚ＤＥが求められる。
【００４４】
図１０は、上記負荷角余弦値cosθＥと析出物厚ＤＥの関係を示したグラフである。
図１０（ａ）は、試験温度Ｔ５、試験時間ｔ４でクリープ試験を行った場合の関係であり、負荷角余弦値cosθＥが大きい結晶粒ＧＥほど析出物厚ＤＥが直線状に大きくなっている。一方、図１０（ｂ）は、試験温度Ｔ５、試験時間ｔ６でクリープ試験を行った場合の関係であり、負荷角余弦値cosθＥが略０．８以上の領域では直線性が崩れている。
【００４５】
この直線性の崩れは、負荷角度θＥが０°に近いので、結晶格子面Ｓｂの受ける負荷応力が引張応力の大きさに近く比較的に大きいことと、長時間の試験によって析出物ＥＥの崩壊が始まったためであると推測される。なお、試験温度Ｔ６、試験時間ｔ４の場合についても、負荷角余弦値cosθＥが１近傍においては直線性が崩れる傾向が認められる。
【００４６】
一方、図１１は、負荷角余弦値cosθＥが０．８の場合と１．０の場合とについて、各試料のラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係を示す点（計測点）を片対数グラフ上に示した特性図である。この図において、白丸は負荷角余弦値cosθＥが１．０の場合の計測点を示し、黒丸は負荷角余弦値cosθＥが０．８の場合を示し、また曲線状の破線は白丸の計測点に最小二乗法を適用することにより求めた近似曲線を、直線状の実線は黒丸の計測点に最小二乗法を適用することにより求めた近似直線を示している。
【００４７】
上記近似曲線は、負荷角余弦値cosθＥが１．０の場合には、ラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係が直線近似することができないと共にその近似性が低いことを示し、また上記近似直線は、負荷角余弦値cosθＥが０．８の場合には、ラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係を直線近似することができると共にその近似性が高いことを示している。
【００４８】
このような図１０及び図１１の特性を総合的に解釈すると、ラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係を高精度に直線近似する、つまり上述した近似性の高い近似式（４）を求める（定数α２、β２を高精度に同定する）ためには、負荷角余弦値cosθＥと析出物厚ＤＥとの関係が直線近似することが可能な負荷角余弦値cosθＥの領域、つまり負荷角余弦値cosθＥが１．０近傍の領域を除外した領域において近似式（４）を求める必要がある。
【００４９】
すなわち、本実施形態では、負荷角余弦値cosθＥの値として、図１０（ｂ）において直線性が良好な負荷角余弦値cosθＥが０．８以下の領域の析出物厚ＤＥをのうち最も大きな値（例えば０．８）を限定的に採用することにより、近似式（４）の定数α２、β２を高精度に同定する。
【００５０】
そして、Ｐ＝ＰＥとして、上記のようにして求めた近似式（４）とラーソン・ミラー・パラメータＰの定義式（１）とを連立させて得られる関係式（５）に、ガスタービン動翼１の使用時間ｔ及び上記のように限定的に採用した負荷角度θＥ（例えば、cosθＥ＝０．８）と同一の負荷角度となる最損傷部位Ｗにおける析出物の析出物厚Ｄを代入することにより、最損傷部位Ｗにおける使用温度Ｔの推定値を求める。このガスタービン動翼１の最損傷部位Ｗにおける析出物の析出物厚Ｄは、垂直な断面（例えば断面Ｓｒ）を露出させて、負荷角度θＥと同一の負荷角度となる結晶粒Ｇから直接的に求められる析出物厚Ｄを用いてもよいし、図１０のように各結晶粒Ｇから析出物厚Ｄと負荷角度の余弦値との関係を求めて、この関係から負荷角度θＥと同一の負荷角度となる析出物厚Ｄを求め、これを用いてもよい。
【００５１】
そして、図６に示すクリープマスターカーブに上記使用温度Ｔの推定値を当てはめることにより、ガスタービン動翼１の未使用状態から破断までの使用時間ｔ（初期破断時間）を求め、この初期破断時間ｔからガスタービン動翼１の使用開始から現在までの使用時間ｔを減算することにより最損傷部位Ｗの破断時間（余寿命）を求める。
【００５２】
続いて、クリープ試験の寿命推定手法について説明する。
クリープ試験の寿命評価では、ガスタービン動翼１の最損傷部位Ｗから複数の試料を作成して、拡散層厚ｄより推定された使用温度Ｔ及び計算で求められた応力σ、あるいはこの使用温度Ｔ及び応力σに近い値の試験温度及び試験応力でクリープ試験を行い、残りの破断時間を特定する。
【００５３】
このような本実施形態によれば、損傷が激しくクリープ破壊が生じる可能性が最も高いと思われる最損傷部位Ｗを特定して、最損傷部位Ｗに対して集中的に破断時間（余寿命）を推定するので、他の部位の破断時間（余寿命）を推定しなくてもよく、迅速且つ高精度に寿命の推定を行うことができる。
【００５４】
また、計測された拡散層厚ｄを近似式（２）に代入して求めた使用温度Ｔの推定値に基づいて垂直断面最損傷部位Ｂ１と特定し、垂直断面最損傷部位Ｂ１を含んで遠心力方向に平行に沿った各部位のうち、計算で求められた応力σと使用温度Ｔの推定値とに基づいて、最損傷部位Ｗを特定するので、最損傷部位Ｗの特定を迅速且つ高精度に特定することができる。
【００５５】
また、本実施形態の拡散層１ｂからの寿命推定手法によれば、近似式（２）、つまりラーソン・ミラー・パラメータＰと拡散層厚ｄとの関係が複数の試料に関する熱時効試験（熱処理）の試験結果から高精度に特定されるので、関係式（３）と拡散層厚ｄとに基づいてガスタービン動翼１の各部位の使用温度Ｔをより迅速且つ精度良く推定することが可能である。よって当該使用温度Ｔの推定値に基づいてガスタービン動翼１の破断時間（余寿命）を高精度に推定することができる。
【００５６】
また、本実施形態の析出物評価からの寿命推定手法よれば、近似式（４）、つまりラーソン・ミラー・パラメータＰＥと析出物（γ´相）の析出物厚ＤＥとの関係が、複数の試料に関するクリープ試験の試験結果から高精度に特定される。すなわち、応力軸ＬＥと直交する断面Ｓａにおける結晶粒ＧＥの析出物ＥＥについて測定された測定厚ＭＥを結晶格子面が接する結晶粒ＧＥの結晶方位ＤＩＥから求められた負荷角余弦値cosθＥを用いて析出物厚ＤＥ（析出物ＥＥの真の厚さ）を求め、当該析出物厚ＤＥに基づいて近似式（４）を求めるので、使用温度Ｔの推定値を精度良く求めることが可能であり、よって破断時間（余寿命）の推定値を精度良く求めることができる。
【００５７】
また、本実施形態の析出物評価からの寿命推定手法よれば、負荷角余弦値cosθＥと析出物厚ＤＥとの関係の直線性が良好な負荷角余弦値cosθＥの領域の析出物厚ＤＥのみに限定して近似式（４）を特定するので、使用温度Ｔの推定値をさらに精度良く求めることが可能であり、よって破断時間（余寿命）の推定値をさらに精度良く求めることができる。
【００５８】
また、本実施形態のクリープ試験の寿命推定手法によれば、最損傷部位Ｗと特定して行うので試料を少なくすることができると共に、信頼性が最も高い余寿命を求めることができる。
【００５９】
なお、本発明は上記実施形態に限定されるものではなく、例えば以下の変形例が考えられる。
（１）上記最損傷部位Ｗの寿命評価は、上記３つのうちいずれか一つでなく、複数を組み合わせて用いてもよい。例えば、コーティングからの寿命推定手法により、迅速な寿命推定結果を得て、この結果に基づいてガスタービンを稼動させつつクリープ試験による寿命推定手法を行って信頼性の高い結果を求めてもよい。
（２）上記実施形態では、断面Ｓｒにおける最も高い温度の垂直断面最損傷部位Ｂ１を求め、垂直断面最損傷部位Ｂ１の応力軸方向の各部位における使用温度及び応力を求めたが、これは必須の構成要件ではなく、最損傷部位Ｗを合理的かつ迅速に特定するためにこのような構成にしたものであるから、その近傍の部位についての破断時間ｔから最損傷部位Ｗを特定してもよい。
また、最損傷部位Ｗは単数のものだけではなく、損傷度が同様の場合は複数となる場合もあり得る。
（３）上記実施形態では、拡散層厚ｄから使用温度を推定したが、例えばＦＥＭ解析を行ってガスタービン動翼１の各部位の使用温度を求めてもよい。また、応力分布についてもＦＥＭ解析の結果を用いてよい。この場合は、アルミニウム合金がコーティングされていないニッケル基合金部品についても本発明を適用することができる。
（４）上記実施形態では、ラーソン・ミラー・パラメータＰＥの算出に際して、別途クリープ試験等で材料定数Ｃを算出してこの値を設定したが、一般的に用いられる「２０」に設定してもよい。
【図面の簡単な説明】
【００６０】
【図１】本発明の一実施形態において、使用開始からある程度時間が経過したガスタービン動翼１を示す図である。
【図２】本発明の一実施形態において、使用開始からある程度時間が経過したガスタービン動翼１の内部状態を示す要部断面図である。
【図３】本発明の一実施形態において、材料定数Ｃを調整した後の各試料のラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥの関係を示すグラフである。
【図４】本発明の一実施形態において、材料定数Ｃを調整する前の各試料のラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥとの関係を示すグラフである。
【図５】本発明の一実施形態において、ガスタービン動翼１の応力σとラーソン・ミラー・パラメータＰとの関係を示すグラフである。
【図６】本発明の一実施形態において、ガスタービン動翼１の未使用状態から破断までの使用時間ｔ（初期破断時間）と使用温度Ｔとの関係を示すグラフである。
【図７】本発明の一実施形態において、ガスタービン動翼１の応力軸Ｌとの直行面における結晶粒界と結晶方位の関係を示す概念図である。
【図８】本発明の一実施形態において、析出物の測定厚ＭＥと析出物厚ＤＥの関係を示す図である。
【図９】本発明の一実施形態において、結晶方位ＤＩＥと応力軸ＬＥの為す負荷角度θＥを示す図である。
【図１０】本発明の一実施形態において、負荷角余弦値cosθＥと析出物厚ＤＥの関係の一例を示した図である。
【図１１】本発明の一実施形態において、負荷角余弦値cosθＥが０．８の場合と１．０の場合とについて、各試料のラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係を示す図である。
【符号の説明】
【００６１】
１…ガスタービン動翼（ニッケル基合金部品）
１ａ…基材
１ｂ…拡散層
Ｂ１…垂直断面最損傷部位
Ｇ，ＧＥ…結晶粒
ＥＥ…析出物
Ｓａ…所定面
Ｓｂ…結晶格子面
Ｌ，ＬＥ…応力軸
Ｄ，ＤＥ…析出物厚
ＤＩＥ…結晶方位
ｄ，ｄＥ…拡散層厚
Ｗ…最損傷部位
【技術分野】
【０００１】
本発明は、ニッケル基合金部品の寿命推定方法に関するものである。
【背景技術】
【０００２】
ニッケル基合金は、高温下における強度と耐クリープ性に優れ、ガスタービン、ジェットエンジン、化学プラント等に使用されている。このような性質を有するニッケル基合金であっても、高温・高応力に曝される使用環境下では、経年劣化による疲労破壊やクリープ破壊等が生じる。従って、このような高温・高応力下で使用されるニッケル基合金部品の余寿命を正確且つ定量的に予測することは、ニッケル基合金部品の検査や交換等の時期を計画する上で非常に重要である。
【０００３】
例えば下記特許文献１には、クリープ試験終了後のニッケル基合金試料から、（１）γ´相（析出物）の円相当直径の平均値とラーソン・ミラー・パラメータとの関係、（２）ラーソン・ミラー・パラメータとラフト化比率（γ´相の応力負荷を受ける方向のアスペクト比と無負荷のアスペクト比の比）との関係を予め求めておき、使用時間が既知の寿命推定対象におけるγ´相の円相当直径及び（１）の関係からラーソン・ミラー・パラメータを求め、このラーソン・ミラー・パラメータ、寿命推定対象を測定して求めたラフト化比率及び（２）の関係に基づいてクリープ破断寿命比（使用時間／未使用状態からの破断時間）を推定する技術が開示されている。
【特許文献１】特開平８−２９７０７９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ところで、例えば論文「三浦信祐、他２名、“多結晶Ｎｉ基超合金，ＩＮ−１００，のクリープ変形により形成されるγ´相のラフト構造”、鉄と鋼、社団法人日本鉄鋼協会、２００４年８月１日、Ｖｏｌ．９０、Ｎｏ．８、ｐ．５６７−５７３」に記載されているように、多結晶ニッケル基合金の結晶方位は各結晶粒によって異なり、またγ´相のラフト化は結晶粒の結晶方位によって異なる。
従って、上記従来技術は、単結晶ニッケル基合金には適用できるものの、多結晶ニッケル基合金には適用することができない。すなわち、多結晶ニッケル基合金では各結晶粒の結晶方位が異なるので、上記従来技術のように各結晶粒の結晶方位の違いを考慮しないで一律に測定した応力負荷面と応力負荷無し面の長径及び短径に基づいてラフト化比率を求める手法は、多結晶ニッケル基合金に適用した場合に大きな誤差を含むものとなる。
【０００５】
本発明は上述した事情に鑑みてなされたものであり、多結晶ニッケル基合金からなるニッケル基合金部品について迅速かつ高精度に寿命を推定することができるニッケル基合金部品の寿命推定方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明は、上記課題を解決するために以下の手段を採用する。
すなわち、本発明は、第１の解決手段として、多結晶ニッケル基合金から形成されたニッケル基合金部品の寿命推定方法であって、前記ニッケル基合金部品の各部位における使用時の使用温度の推定値と応力の推定値とに基づいて、相対的に損傷の度合いが大きい最損傷部位を特定する第１工程と、前記最損傷部位について１または複数の手法により破断までの破断時間を求める第２工程とを有する、という手段を採用する。
【０００７】
また、本発明は、第２の解決手段として、上記第１の解決手段における前記第１工程が、前記ニッケル基合金部品の使用時の応力軸方向に垂直な断面において前記使用温度の推定値が最も高い垂直断面最損傷部位を特定する第１Ａ工程と、前記垂直断面最損傷部位の応力軸方向に位置する前記各部位のうち、前記使用温度の推定値と、前記応力の推定値と、前記ニッケル基合金部品の使用開始から経過した使用時間と、前記ニッケル基合金部品の同等部品について予め取得されたクリープマスターカーブとに基づいて、前記最損傷部位を特定する第１Ｂ工程と、からなるという手段を採用する。
【０００８】
また、本発明は、第３の解決手段として、上記第１又は第２の解決手段における前記ニッケル基合金部品が前記ニッケル基合金からなる基材の表面に金属コーティング層が形成されたものである場合、前記第１工程では、前記ニッケル基合金部品と同等な複数の試料について、前記ニッケル基合金部品の使用によって前記基材に形成される拡散層の厚さ（拡散層厚）とラーソン・ミラー・パラメータとの関係をそれぞれ求める第１Ｃ工程と、各試料に関する前記関係に基づいて拡散層厚とラーソン・ミラー・パラメータとの関係を示す第１の近似式を求める第１Ｄ工程と、前記第１の近似式とラーソン・ミラー・パラメータの定義式との第１の連立方程式に、前記各部位における拡散層厚の測定値と、前記使用時間とを代入する第１Ｅ工程と、から前記使用温度の推定値を求める、という手段を採用する。
また、本発明は、第４の解決手段として、上記第３の解決手段における前記第２工程では、前記第１Ｃ工程から第１Ｅ工程によって求められた前記使用温度の推定値を、前記クリープマスターカーブに当てはめることにより前記最損傷部位の破断時間を求める、という手段を採用する。
【０００９】
また、本発明は、第５の解決手段として、上記第１から４のうちいずれかの解決手段における前記第２工程が、異なる試験温度でクリープ試験を行った前記ニッケル基合金部品と同等な複数の試料について、前記試料の応力軸に直交する直交面における複数の結晶粒の結晶方位と結晶粒の結晶格子面に沿って析出した複数の析出物の厚さとをそれぞれ計測すると共にラーソン・ミラー・パラメータを求める第２Ａ工程と、前記結晶粒の結晶方位と析出物の厚さとに基づいて析出物の真の厚さ（析出物厚）を求める第２Ｂ工程と、前記複数の試料のラーソン・ミラー・パラメータと、前記結晶方位と前記応力軸とのなす負荷角度の余弦値との関係の直線性が良好な領域における特定の前記負荷角度の前記析出物厚と、に基づいて第２の近似式を求める第２Ｃ工程と、前記第２の近似式とラーソン・ミラー・パラメータの定義式との第２の連立方程式に、前記最損傷部位から得られる前記特定の負荷角度と同一の負荷角度における析出物厚と、前記ニッケル基合金部品の使用開始から経過した使用時間と、を代入することにより前記最損傷部位における使用温度の推定値を求める第２Ｄ工程と、このようにして求めた前記最損傷部位の使用温度の推定値を、前記クリープマスターカーブに当てはめることにより前記最損傷部位における破断時間を求める第２Ｅ工程と、に基づいて前記最損傷部位の破断時間を求める、という手段を採用する。
【００１０】
また、本発明は、第６の解決手段として、上記第１から５の解決手段における前記第２工程では、前記最損傷部位から作成した試料に対してクリープ試験を行って破断時間を求める、という手段を採用する。
また、本発明は、第７の解決手段として、上記第１から６の解決手段における前記ニッケル基合金部品がガスタービンの動翼または静翼である、という手段を採用する。
【発明の効果】
【００１１】
本発明によれば、ニッケル基合金部品の各部位のうち相対的に損傷の度合いが大きい最損傷部位を特定し、最損傷部位について１または複数の手法により破断時間（寿命）の推定を行うので、最損傷部位以外の部位における破断時間を推定する必要が無く、クリープ破壊の生じる可能性が最も高い部位に集中して１または複数の手法で寿命の推定を行うことができ、迅速且つ高精度に寿命を推定することが可能となる。
【００１２】
また、前記第１工程は、垂直断面最損傷部位の応力軸方向に位置する各部位のうち、使用温度の推定値と、応力の推定値と、ニッケル基合金部品の使用開始から経過した使用時間と、前記ニッケル基合金部品の同等部品について予め取得されたクリープマスターカーブとに基づいて最損傷部位を特定するので、最損傷部位を迅速に特定することができ、より迅速且つ高精度にニッケル基合金部品の寿命の推定を行うことができる。
【発明を実施するための最良の形態】
【００１３】
以下、本発明の一実施形態について、図面を参照して説明する。本実施形態では、ガスタービンにおける圧縮機やタービンの構成部品となる動翼（ガスタービン動翼）を寿命推定対象とする。
【００１４】
図１は、使用開始から所定時間が経過したガスタービン動翼１を示す図であり、図２は、使用開始から所定時間が経過したガスタービン動翼１の要部断面図である。ガスタービン動翼１は、図１に示すように翼型断面を有する棒状部品であり、タービン軸（回転軸）に放射状に取り付けられて使用される。また、ガスタービン動翼１は、図２に示すように、多結晶ニッケル基合金からなる基材１ａ上に拡散層１ｂ、コーティング層１ｃ及び酸化膜層１ｄが形成されている。ガスタービン動翼１はタービン回転軸に放射状に取り付けられて使用に供されるので、応力軸Ｌは、図示するようにガスタービン動翼１の長さ方向となる。
【００１５】
本来（使用開始前）のガスタービン動翼１は、基材１ａの表面にアルミニウム（Ａｌ）からなるコーティング層１ｃが形成されたものであるが、高温環境下での使用によってコーティング層１ｃのアルミニウムが基材１ａ内に拡散することによって拡散層１ｂが形成されると共に、酸化によってコーティング層１ｃの表面に酸化膜層１ｄが形成される。
【００１６】
本実施形態は、ガスタービン動翼１の全体を評価対象とするのではなく、ガスタービン動翼１についてクリープ損傷が最も激しい部位（最損傷部位）を最初に特定し、当該最損傷部位について寿命の推定を行うものである。
より具体的には、ガスタービン動翼１に作用する遠心力Ｆの方向（応力軸Ｌの方向）に垂直な断面Ｓｒにおけるガスタービン動翼１の表面（翼状表面）の使用時における温度分布（使用温度の分布）を推定し、この断面Ｓｒにおいて使用温度が最も高い部位を垂直断面最損傷部位Ｂ１として特定する。すなわち、断面Ｓｒは応力軸Ｌに垂直なので、断面Ｓｒにおけるガスタービン動翼１の表面の各部位には同一の応力が作用すると考えることができ、よって上記表面の各部位のうち使用温度が最も高い部位Ｂ１は、断面Ｓｒにおいてクリープ損傷が最も激しい部位（垂直断面最損傷部位Ｂ１）とみなすことができる。
【００１７】
次に、垂直断面最損傷部位Ｂ１の応力軸方向におけるガスタービン動翼１の表面についての使用温度の分布と使用時応力の分布とをそれぞれ推定する。すなわち、遠心力Ｆの方向に垂直な各断面におけるガスタービン動翼１の表面の使用温度の分布の状態はほぼ同様のものとなるので、垂直断面最損傷部位Ｂ１の遠心力方向に位置する各部位が各断面において最も使用温度が高い部位であるとみなすことができる。また、上述の通り各断面において生じる応力はその断面において同一と考えられるから、垂直断面最損傷部位Ｂ１の遠心力Ｆの方向に位置する各部位が各断面においてクリープ損傷が最も激しい部位であるとみなすことができる。従って、この垂直断面最損傷部位Ｂ１の遠心力方向に位置する各部位の中で最もクリープ損傷が激しい部位が最損傷部位Ｗとなる。
【００１８】
そして、垂直断面最損傷部位Ｂ１の遠心力方向に位置する各部位における使用時応力（推定値）から、予めガスタービン動翼１と同等の試料から取得したクリープ破断データを参照してラーション・ミラー・パラメータＰを求めて、このラーション・ミラー・パラメータＰと、垂直断面最損傷部位Ｂ１の遠心力方向に位置する各部位における使用温度（推定値）とを下記ラーソン・ミラー・パラメータ定義式（１）に代入することにより各破断時間を算出し、この破断時間が最も短い部位をガスタービン動翼１の全体において損傷が最も激しい部位Ｗ（最損傷部位）と特定する。そして、本実施形態では、この最損傷部位Ｗのみについて破断時間（余寿命）の推定を行う。
【００１９】
まず、上述のように、初めにガスタービン動翼１に作用する遠心力方向に垂直な断面Ｓｒの温度分布を求める。本実施形態においては、基材１ａ上に形成される拡散層１ｂの厚さｄ（拡散層厚）から断面Ｓｒの温度分布を求める。
具体的には、以下に説明するように、ラーソン・ミラー・パラメータＰの定義式（１）と、使用開始前のガスタービン動翼１と同等な複数の試料に関する熱時効試験（熱処理）の試験結果から同定されたラーソン・ミラー・パラメータＰＥと拡散層１ｂの厚さｄＥ（拡散層厚）との近似式（２）と、Ｐ＝ＰＥとして上記定義式（１）と近似式（２）とを連立させて得られる関係式（３）に基づいてガスタービン動翼１の使用温度Ｔの推定値を求め、この使用温度Ｔの推定値をガスタービン動翼１の使用温度と初期破断時間とに相当する関係（クリープマスターカーブ）に当てはめることにより破断時間（寿命）を推定するものである。なお、各変数の添字Ｅは、試料に関するものを意味するが、関係式（３）においてはガスタービン動翼１の拡散層厚ｄを代入するので、ｄＥをｄとしている。
Ｐ＝Ｔ・（Ｃ＋ｌｏｇｔ） （１）
ＰＥ＝（ｄＥ−β１）／α１ （２）
Ｔ＝（ｄ−β１）／{α１（Ｃ＋ｌｏｇｔ）} （３）
【００２０】
ラーソン・ミラー・パラメータは、周知のように、金属部品のクリープ破壊における温度と時間との関係を示すパラメータとして、上記定義式（１）で表される。この定義式（１）において、Ｔは使用温度（絶対温度）、ｔは使用時間、Ｃは材料固有の材料定数である。
【００２１】
また、近似式（２）におけるα１，β１は、上記試料に関する実験結果に基づいて同定される定数である。本実施形態では、以下に説明するように、定数α１，β１をより正確に同定するために、実際に使用された複数のガスタービン動翼１から取得した拡散層厚ｄの分布範囲に該当する熱時効試験（熱処理）の試験結果に限定して材料定数Ｃを最適化する。
【００２２】
上記熱時効試験（熱処理）では、上記複数の試料について試験温度ＴＥと試験時間ｔＥを変えた熱時効試験（熱処理）を行う。試験温度ＴＥは、大まかに把握されるガスタービン動翼１の実際の使用温度範囲内及びその近傍における複数温度、例えば４つの試験温度Ｔ１，Ｔ２，Ｔ３，Ｔ４が設定され、一方、試験時間ｔＥについては、ガスタービン動翼１の実際の使用時間範囲内における複数の時間、例えば３つの試験時間ｔ１，ｔ２，ｔ３が設定される。すなわち、熱時効試験（熱処理）では、各試料について４つの試験温度Ｔ１，Ｔ２，Ｔ３，Ｔ４と３つの試験時間ｔ１，ｔ２，ｔ３とが設定される。
【００２３】
そして、上記定義式（１）における使用温度Ｔに試験温度Ｔ１，Ｔ２，Ｔ３，Ｔ４をそれぞれ代入すると共に使用時間ｔに試験時間ｔ１，ｔ２，ｔ３をそれぞれ代入することにより、各試料に関するラーソン・ミラー・パラメータＰＥを求める。ラーソン・ミラー・パラメータＰＥの算出に際しては、別途クリープ試験等で材料定数Ｃを算出してこの値を設定する。
【００２４】
続いて、上記各試料を研磨して内部を露出させ、顕微鏡を用い拡散層厚ｄＥを計測する。この拡散層厚ｄＥは、顕微鏡から出力される試料の画像について画像処理装置を用いて画像処理することによって、あるいは検査員が試料の画像を目視することによって計測される。なお、顕微鏡は拡散層厚ｄＥを計測するのに適切な分解能を有していればよく、例えば光学顕微鏡や走査型電子顕微鏡を用いることができる。
【００２５】
図３は、このようにして求められた各試料のラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥとの関係を示す点（計測点）を片対数グラフ上に示した図である。この図では、試料毎に異なる表記方法で計測点を示している。本実施形態では、このような各試料に関する計測点の全てに例えば最小二乗法を適用することにより、全ての試料に関するラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥとの関係を示す特性直線を求める。このようにして各試料に関する熱時効試験（熱処理）の試験結果から得られた特性直線は、上述した近似式（２）において定数α１、β１が同定されたものである。
【００２６】
ここで、図３に破線で示すように、同一温度における試料の計測点は良好な直線性を示す。つまり各温度については近似性の高い特性直線を求めることができるが、全ての温度の計測点を１本の特性直線を求めるための要素とした場合には、異なる温度域が含まれるので、近似性の高い１本の特性直線を求めることは困難である。
【００２７】
本実施形態では、このような事情を回避して、より近似性の高い特性直線Ｎ、つまり近似式（２）を求めるために、上記と同様に顕微鏡を用いて実際に使用された複数のガスタービン動翼１から拡散層厚ｄの分布範囲Ｈを予め取得し、各試料に関するラーソン・ミラー・パラメータＰＥを求める際に、分布範囲Ｈに該当する試料の熱時効試験（熱処理）の試験結果（計測点）のみに限定して材料定数Ｃを最適化する。すなわち、近似式（２）の近似性が最も高くなるように、図３に示すように各試料に関する計測点を分布範囲Ｈ内に狭く限定し、この上で各試料に関するラーソン・ミラー・パラメータＰＥを求める際の材料定数Ｃを最適化する。なお、関係式（３）における材料定数Ｃも上記最適化したものを用いる。
【００２８】
図４は、このような材料定数Ｃを最適化処理前の各試料のラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥとの関係を示す計測点を片対数グラフ上に示した図である。この図４に示すように、材料定数Ｃについて最適化処理を行わない場合、各試料の計測点に関する特性直線は、分布範囲Ｈ内において比較的乖離しており、各試料の全計測点から求められる１本の特性直線は近似性が低いものとなる。
【００２９】
本実施形態では、近似式（２）を特定することによって得られた近似性の高い関係式（３）に、ガスタービン動翼１に作用する遠心力方向に垂直な断面Ｓｒにおける各周縁の拡散層厚ｄを上記と同様に顕微鏡観察により求め、この拡散層厚ｄを関係式（３）に代入することにより、断面Ｓｒにおけるガスタービン動翼１の表面（翼状表面）の使用時における温度分布（使用温度Ｔの分布）を求める。
例えば、図１において、断面Ｓｒの縁に沿って拡散層厚ｄを計測していき、この計測された拡散層厚ｄを順次関係式（３）に代入して、断面Ｓｒにおけるガスタービン動翼１の表面（翼状表面）の使用時における温度分布（使用温度Ｔの分布）を推定し、この断面Ｓｒにおいて使用温度が最も高い部位を垂直断面最損傷部位Ｂ１として特定する。
【００３０】
次に、垂直断面最損傷部位Ｂ１を特定した後に、ガスタービン動翼１の垂直断面最損傷部位Ｂ１を含む垂直断面最損傷部位Ｂ１の遠心力方向の各部位における使用温度Ｔの推定値と応力σの推定値を求める。図１においては、遠心力方向に沿う全ての部位のうち部位Ｂ１１〜Ｂ１４を例示として示すが、遠心力方向に沿う全ての部位について使用温度Ｔを推定して構わない。ここで、使用温度Ｔの推定値については、断面Ｓｒの場合と同様に拡散層厚ｄと近似式（３）から、応力σについては、その部位に生じる遠心力Ｆ及びその部位から遠心力方向側端面までの棒の自重並びに断面Ｓｒの断面積に基づいて求めている。なお、遠心力方向に垂直な各断面の断面積は一定としている。
【００３１】
図５は、上記ガスタービン動翼１と同等な動翼にクリープ試験を行って得られたクリープマスターカーブを示しており、試料に負荷した応力σとラーソン・ミラー・パラメータＰとの関係に相当する特性直線を示している。なお、このクリープマスターカーブのラーソン・ミラー・パラメータＰは、未使用状態から破断までの使用時間ｔ（初期破断時間）と使用温度Ｔとから求めている。
【００３２】
この図５に示すクリープマスターカーブに、上記のようにして求めた各部位における応力σの推定値を当てはめることにより、各部位におけるラーソン・ミラー・パラメータＰを求め、この各ラーソン・ミラー・パラメータＰと各部位の使用温度Ｔとラーソン・ミラー・パラメータの定義式（１）より未使用状態からの初期破断時間を求める。さらに、初期破断時間から使用時間ｔを減算することにより、各部位の破断までの破断時間（余寿命）を算出する。そして、各部位の中で破断時間（余寿命）が最も短い部位を最損傷部位Ｗとして特定する。
例えば関係式（３）により、ある部位Ｂ１１における応力σの推定値σ１１が求まれば、クリープマスターカーブにより、応力σ１１に対応するラーソン・ミラー・パラメータＰ１１を求めることができる。このラーソン・ミラー・パラメータＰ１１と、部位Ｂ１１における使用温度Ｔの推定値と、ラーソン・ミラー・パラメータの定義式（１）に代入することにより、部位Ｂ１１における初期破断時間を求めることができ、この初期破断時間の推定値から使用時間を減算することにより、破断時間（余寿命）を推定することができる。
【００３３】
上記のように最損傷部位Ｗを特定した後、拡散層１ｂからの寿命推定手法、析出物評価からの寿命推定手法、クリープ試験からの寿命推定手法のうちいずれかを行う。但し、破断時間（余寿命）の推定結果は、拡散層１ｂからの寿命推定手法、析出物評価からの寿命推定手法、クリープ試験に基づく寿命推定手法の順に早く求めることができるが、逆にその結果の信頼性は順に低いものとなっている。そのため、ガスタービン動翼１の作動環境や効率、使用時間等によっていずれかの上記寿命推定手法が適宜選択されて破断時間（余寿命）が求められる。
【００３４】
図６は、上記ガスタービン動翼１と同等な動翼にクリープ試験を行って得られたクリープマスターカーブを示しており、上記ガスタービン動翼１の未使用状態から破断までの使用時間ｔ（初期破断時間）と使用温度Ｔとの関係に相当する特性曲線を示している。
【００３５】
以下、各寿命推定手法について説明する。
まず、初めに拡散層１ｂからの寿命推定手法について説明する。
拡散層１ｂからの寿命推定手法は、上記のように求めた最損傷部位Ｗにおける拡散層厚ｄを測定値から求めた使用温度Ｔの推定値をクリープマスターカーブに当てはめることによってガスタービン動翼１の未使用状態から破断までの使用時間ｔ（初期破断時間）の推定値を求める。そして、ガスタービン動翼１の使用開始から現在までの使用時間ｔを上記初期破断時間の推定値から減算することによりガスタービン動翼１の破断時間（余寿命）を求める。例えば関係式（３）により、最損傷部位Ｗの使用温度ＴＳが求まると、クリープマスターカーブにより使用温度ＴＳにおける初期破断時間ｔＯの推定値が求められる。そして、この初期破断時間ｔＯの推定値からガスタービン動翼１の使用時間ｔＳを減算すれば、最損傷部位Ｗの破断時間（余寿命）ｔＮが推定される。
【００３６】
続いて、析出物評価からの寿命推定手法について説明する。
図７は、使用開始から所定時間が経過したガスタービン動翼１の上記応力軸Ｌと直行する断面Ｓｒにおける結晶粒と結晶方位の関係を示す測定結果である。この図７に示すように、ガスタービン動翼１の応力軸Ｌに直交する断面Ｓｒには、多数の結晶粒Ｇが分布している。また、これら結晶粒Ｇは、図示する立方体が示すように互いに異なる結晶方位を有している。
なお、ガスタービン動翼１は、使用開始から時間が経過するに従って、未使用時には立方体状であった結晶粒内の析出物（γ´相）が板状となって、結晶粒の結晶格子面に沿って粗大化することが知られている。
【００３７】
以下に説明するように、ラーソン・ミラー・パラメータＰの定義式（１）と、使用開始前のガスタービン動翼１と同等な複数の試料に関するクリープ試験の試験結果から同定されたラーソン・ミラー・パラメータＰＥと析出物の真の厚さ（析出物厚ＤＥ）との近似式（４）と、Ｐ＝ＰＥとして上記定義式（１）と近似式（４）とを連立させて得られる関係式（５）に基づいてガスタービン動翼１の使用温度Ｔの推定値を求め、この使用温度Ｔの推定値をガスタービン動翼１の使用温度Ｔと初期破断時間との関係を示す特性図（図６に示すクリープマスターカーブ）に当てはめることにより破断時間（余寿命）を推定するものである。
Ｐ＝Ｔ・（Ｃ＋ｌｏｇｔ） （１）
ＰＥ＝（ＤＥ−β２）／α２ （４）
Ｔ＝（Ｄ−β２）／{α２（Ｃ＋ｌｏｇｔ）} （５）
【００３８】
また、本実施形態では、上記近似式（４）を求めるに際して、その近似性を高めるために、各々の結晶粒の結晶方位を考慮することによって、正確な析出物厚ＤＥを求めると共に、析出物厚ＤＥと結晶方位の関係から定数α２，β２をより正確に同定する。なお、各変数の添字Ｅは、上記と同様に試料に関するものを意味するが、関係式（５）においてはガスタービン動翼１の析出物厚Ｄを代入するので、ＤＥをＤと表記している。
【００３９】
上記クリープ試験では、上記複数の試料について、同一の引張（負荷）応力σＥについて異なる試験温度ＴＥが設定される。試験温度ＴＥは、ガスタービン動翼１の実際の使用温度範囲内における複数温度、例えば３つの試験温度Ｔ４，Ｔ５，Ｔ６が設定され、一方、引張応力については、例えばガスタービン動翼１に作用する実際の負荷応力、例えばガスタービン動翼１に加わる遠心力Ｆを断面Ｓｒの断面積で除した値が設定される。
【００４０】
そして、上記定義式（１）における使用温度Ｔに試験温度Ｔ４，Ｔ５，Ｔ６（Ｔ４＜Ｔ５＜Ｔ６）をそれぞれ代入すると共に、使用時間ｔに各試料の試験時間ｔ４，ｔ５，ｔ６（ｔ４＜ｔ５＜ｔ６）をそれぞれ代入することにより、各試料に関するラーソン・ミラー・パラメータＰＥを求める。
【００４１】
図８は、上記試料における析出物の測定厚ＭＥと析出物厚ＤＥとの関係を示す模式図であり、符号ＧＥは結晶粒を、またＥＥは析出物をそれぞれ示している。上述したように各試料についてラーソン・ミラー・パラメータＰＥを求めて、各試料に作用する応力の応力軸ＬＥに直交する断面Ｓａを露出させて、走査型電子顕微鏡を用いることにより当該断面Ｓａの各析出物ＥＥについて測定厚ＭＥをそれぞれ測定すると共に、走査型電子顕微鏡にＥＢＳＰ（Electron Back Scattering Pattern）検出器を接続することにより上記断面Ｓａの各結晶粒ＧＥについて結晶方位ＤＩＥ（結晶格子面Ｓｂの法線方向）をそれぞれ測定する。
【００４２】
図９は、上記図８における結晶粒ＧＥの結晶方位ＤＩＥ（結晶格子面Ｓｂの法線方向）と応力軸ＬＥとのなす角度θＥ（負荷角度）を示す模式図である。
各結晶粒ＧＥの断面Ｓａは、応力軸ＬＥに直交する面であるが、結晶粒ＧＥの結晶方位ＤＩＥは、応力軸ＬＥに関係なく個々の方向となる。すなわち、負荷角度θＥは、結晶粒ＧＥ毎に異なった値となる。従って、各結晶粒ＧＥについて得られた測定厚ＭＥに同じく各結晶粒ＧＥについて得られた負荷角度θＥの余弦値cosθＥ（負荷角余弦値）を乗算することにより、各結晶粒ＧＥについて結晶方位ＤＩＥの方向における析出物ＥＥの真の厚さ（析出物厚ＤＥ）を求める。
【００４３】
以上の処理によって、試料毎にラーソン・ミラー・パラメータＰＥが求められると共に、各試料の断面Ｓａにおける複数の結晶粒ＧＥについて結晶方位ＤＩＥが求められ、また複数の結晶粒ＧＥの結晶格子面Ｓｂに沿って各々粗大化する複数の析出物ＥＥについて析出物厚ＤＥが求められる。
【００４４】
図１０は、上記負荷角余弦値cosθＥと析出物厚ＤＥの関係を示したグラフである。
図１０（ａ）は、試験温度Ｔ５、試験時間ｔ４でクリープ試験を行った場合の関係であり、負荷角余弦値cosθＥが大きい結晶粒ＧＥほど析出物厚ＤＥが直線状に大きくなっている。一方、図１０（ｂ）は、試験温度Ｔ５、試験時間ｔ６でクリープ試験を行った場合の関係であり、負荷角余弦値cosθＥが略０．８以上の領域では直線性が崩れている。
【００４５】
この直線性の崩れは、負荷角度θＥが０°に近いので、結晶格子面Ｓｂの受ける負荷応力が引張応力の大きさに近く比較的に大きいことと、長時間の試験によって析出物ＥＥの崩壊が始まったためであると推測される。なお、試験温度Ｔ６、試験時間ｔ４の場合についても、負荷角余弦値cosθＥが１近傍においては直線性が崩れる傾向が認められる。
【００４６】
一方、図１１は、負荷角余弦値cosθＥが０．８の場合と１．０の場合とについて、各試料のラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係を示す点（計測点）を片対数グラフ上に示した特性図である。この図において、白丸は負荷角余弦値cosθＥが１．０の場合の計測点を示し、黒丸は負荷角余弦値cosθＥが０．８の場合を示し、また曲線状の破線は白丸の計測点に最小二乗法を適用することにより求めた近似曲線を、直線状の実線は黒丸の計測点に最小二乗法を適用することにより求めた近似直線を示している。
【００４７】
上記近似曲線は、負荷角余弦値cosθＥが１．０の場合には、ラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係が直線近似することができないと共にその近似性が低いことを示し、また上記近似直線は、負荷角余弦値cosθＥが０．８の場合には、ラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係を直線近似することができると共にその近似性が高いことを示している。
【００４８】
このような図１０及び図１１の特性を総合的に解釈すると、ラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係を高精度に直線近似する、つまり上述した近似性の高い近似式（４）を求める（定数α２、β２を高精度に同定する）ためには、負荷角余弦値cosθＥと析出物厚ＤＥとの関係が直線近似することが可能な負荷角余弦値cosθＥの領域、つまり負荷角余弦値cosθＥが１．０近傍の領域を除外した領域において近似式（４）を求める必要がある。
【００４９】
すなわち、本実施形態では、負荷角余弦値cosθＥの値として、図１０（ｂ）において直線性が良好な負荷角余弦値cosθＥが０．８以下の領域の析出物厚ＤＥをのうち最も大きな値（例えば０．８）を限定的に採用することにより、近似式（４）の定数α２、β２を高精度に同定する。
【００５０】
そして、Ｐ＝ＰＥとして、上記のようにして求めた近似式（４）とラーソン・ミラー・パラメータＰの定義式（１）とを連立させて得られる関係式（５）に、ガスタービン動翼１の使用時間ｔ及び上記のように限定的に採用した負荷角度θＥ（例えば、cosθＥ＝０．８）と同一の負荷角度となる最損傷部位Ｗにおける析出物の析出物厚Ｄを代入することにより、最損傷部位Ｗにおける使用温度Ｔの推定値を求める。このガスタービン動翼１の最損傷部位Ｗにおける析出物の析出物厚Ｄは、垂直な断面（例えば断面Ｓｒ）を露出させて、負荷角度θＥと同一の負荷角度となる結晶粒Ｇから直接的に求められる析出物厚Ｄを用いてもよいし、図１０のように各結晶粒Ｇから析出物厚Ｄと負荷角度の余弦値との関係を求めて、この関係から負荷角度θＥと同一の負荷角度となる析出物厚Ｄを求め、これを用いてもよい。
【００５１】
そして、図６に示すクリープマスターカーブに上記使用温度Ｔの推定値を当てはめることにより、ガスタービン動翼１の未使用状態から破断までの使用時間ｔ（初期破断時間）を求め、この初期破断時間ｔからガスタービン動翼１の使用開始から現在までの使用時間ｔを減算することにより最損傷部位Ｗの破断時間（余寿命）を求める。
【００５２】
続いて、クリープ試験の寿命推定手法について説明する。
クリープ試験の寿命評価では、ガスタービン動翼１の最損傷部位Ｗから複数の試料を作成して、拡散層厚ｄより推定された使用温度Ｔ及び計算で求められた応力σ、あるいはこの使用温度Ｔ及び応力σに近い値の試験温度及び試験応力でクリープ試験を行い、残りの破断時間を特定する。
【００５３】
このような本実施形態によれば、損傷が激しくクリープ破壊が生じる可能性が最も高いと思われる最損傷部位Ｗを特定して、最損傷部位Ｗに対して集中的に破断時間（余寿命）を推定するので、他の部位の破断時間（余寿命）を推定しなくてもよく、迅速且つ高精度に寿命の推定を行うことができる。
【００５４】
また、計測された拡散層厚ｄを近似式（２）に代入して求めた使用温度Ｔの推定値に基づいて垂直断面最損傷部位Ｂ１と特定し、垂直断面最損傷部位Ｂ１を含んで遠心力方向に平行に沿った各部位のうち、計算で求められた応力σと使用温度Ｔの推定値とに基づいて、最損傷部位Ｗを特定するので、最損傷部位Ｗの特定を迅速且つ高精度に特定することができる。
【００５５】
また、本実施形態の拡散層１ｂからの寿命推定手法によれば、近似式（２）、つまりラーソン・ミラー・パラメータＰと拡散層厚ｄとの関係が複数の試料に関する熱時効試験（熱処理）の試験結果から高精度に特定されるので、関係式（３）と拡散層厚ｄとに基づいてガスタービン動翼１の各部位の使用温度Ｔをより迅速且つ精度良く推定することが可能である。よって当該使用温度Ｔの推定値に基づいてガスタービン動翼１の破断時間（余寿命）を高精度に推定することができる。
【００５６】
また、本実施形態の析出物評価からの寿命推定手法よれば、近似式（４）、つまりラーソン・ミラー・パラメータＰＥと析出物（γ´相）の析出物厚ＤＥとの関係が、複数の試料に関するクリープ試験の試験結果から高精度に特定される。すなわち、応力軸ＬＥと直交する断面Ｓａにおける結晶粒ＧＥの析出物ＥＥについて測定された測定厚ＭＥを結晶格子面が接する結晶粒ＧＥの結晶方位ＤＩＥから求められた負荷角余弦値cosθＥを用いて析出物厚ＤＥ（析出物ＥＥの真の厚さ）を求め、当該析出物厚ＤＥに基づいて近似式（４）を求めるので、使用温度Ｔの推定値を精度良く求めることが可能であり、よって破断時間（余寿命）の推定値を精度良く求めることができる。
【００５７】
また、本実施形態の析出物評価からの寿命推定手法よれば、負荷角余弦値cosθＥと析出物厚ＤＥとの関係の直線性が良好な負荷角余弦値cosθＥの領域の析出物厚ＤＥのみに限定して近似式（４）を特定するので、使用温度Ｔの推定値をさらに精度良く求めることが可能であり、よって破断時間（余寿命）の推定値をさらに精度良く求めることができる。
【００５８】
また、本実施形態のクリープ試験の寿命推定手法によれば、最損傷部位Ｗと特定して行うので試料を少なくすることができると共に、信頼性が最も高い余寿命を求めることができる。
【００５９】
なお、本発明は上記実施形態に限定されるものではなく、例えば以下の変形例が考えられる。
（１）上記最損傷部位Ｗの寿命評価は、上記３つのうちいずれか一つでなく、複数を組み合わせて用いてもよい。例えば、コーティングからの寿命推定手法により、迅速な寿命推定結果を得て、この結果に基づいてガスタービンを稼動させつつクリープ試験による寿命推定手法を行って信頼性の高い結果を求めてもよい。
（２）上記実施形態では、断面Ｓｒにおける最も高い温度の垂直断面最損傷部位Ｂ１を求め、垂直断面最損傷部位Ｂ１の応力軸方向の各部位における使用温度及び応力を求めたが、これは必須の構成要件ではなく、最損傷部位Ｗを合理的かつ迅速に特定するためにこのような構成にしたものであるから、その近傍の部位についての破断時間ｔから最損傷部位Ｗを特定してもよい。
また、最損傷部位Ｗは単数のものだけではなく、損傷度が同様の場合は複数となる場合もあり得る。
（３）上記実施形態では、拡散層厚ｄから使用温度を推定したが、例えばＦＥＭ解析を行ってガスタービン動翼１の各部位の使用温度を求めてもよい。また、応力分布についてもＦＥＭ解析の結果を用いてよい。この場合は、アルミニウム合金がコーティングされていないニッケル基合金部品についても本発明を適用することができる。
（４）上記実施形態では、ラーソン・ミラー・パラメータＰＥの算出に際して、別途クリープ試験等で材料定数Ｃを算出してこの値を設定したが、一般的に用いられる「２０」に設定してもよい。
【図面の簡単な説明】
【００６０】
【図１】本発明の一実施形態において、使用開始からある程度時間が経過したガスタービン動翼１を示す図である。
【図２】本発明の一実施形態において、使用開始からある程度時間が経過したガスタービン動翼１の内部状態を示す要部断面図である。
【図３】本発明の一実施形態において、材料定数Ｃを調整した後の各試料のラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥの関係を示すグラフである。
【図４】本発明の一実施形態において、材料定数Ｃを調整する前の各試料のラーソン・ミラー・パラメータＰＥと拡散層厚ｄＥとの関係を示すグラフである。
【図５】本発明の一実施形態において、ガスタービン動翼１の応力σとラーソン・ミラー・パラメータＰとの関係を示すグラフである。
【図６】本発明の一実施形態において、ガスタービン動翼１の未使用状態から破断までの使用時間ｔ（初期破断時間）と使用温度Ｔとの関係を示すグラフである。
【図７】本発明の一実施形態において、ガスタービン動翼１の応力軸Ｌとの直行面における結晶粒界と結晶方位の関係を示す概念図である。
【図８】本発明の一実施形態において、析出物の測定厚ＭＥと析出物厚ＤＥの関係を示す図である。
【図９】本発明の一実施形態において、結晶方位ＤＩＥと応力軸ＬＥの為す負荷角度θＥを示す図である。
【図１０】本発明の一実施形態において、負荷角余弦値cosθＥと析出物厚ＤＥの関係の一例を示した図である。
【図１１】本発明の一実施形態において、負荷角余弦値cosθＥが０．８の場合と１．０の場合とについて、各試料のラーソン・ミラー・パラメータＰＥと析出物厚ＤＥとの関係を示す図である。
【符号の説明】
【００６１】
１…ガスタービン動翼（ニッケル基合金部品）
１ａ…基材
１ｂ…拡散層
Ｂ１…垂直断面最損傷部位
Ｇ，ＧＥ…結晶粒
ＥＥ…析出物
Ｓａ…所定面
Ｓｂ…結晶格子面
Ｌ，ＬＥ…応力軸
Ｄ，ＤＥ…析出物厚
ＤＩＥ…結晶方位
ｄ，ｄＥ…拡散層厚
Ｗ…最損傷部位
【特許請求の範囲】
【請求項１】
多結晶ニッケル基合金から形成されたニッケル基合金部品の寿命推定方法であって、
前記ニッケル基合金部品の各部位における使用時の使用温度の推定値と応力の推定値とに基づいて、相対的に損傷の度合いが大きい最損傷部位を特定する第１工程と、
前記最損傷部位について１または複数の手法により破断までの破断時間を求める第２工程と、
を有することを特徴とするニッケル基合金部品の寿命推定方法。
【請求項２】
前記第１工程は、
前記ニッケル基合金部品の使用時の応力軸方向に垂直な断面において前記使用温度の推定値が最も高い垂直断面最損傷部位を特定する第１Ａ工程と、
前記垂直断面最損傷部位の応力軸方向に位置する前記各部位のうち、前記使用温度の推定値と、前記応力の推定値と、前記ニッケル基合金部品の使用開始から経過した使用時間と、前記ニッケル基合金部品の同等部品について予め取得されたクリープマスターカーブとに基づいて、前記最損傷部位を特定する第１Ｂ工程と、
からなることを特徴とする請求項１に記載のニッケル基合金部品の寿命推定方法。
【請求項３】
前記ニッケル基合金部品が前記ニッケル基合金からなる基材の表面に金属コーティング層が形成されたものである場合、
前記第１工程では、
前記ニッケル基合金部品と同等な複数の試料について、前記ニッケル基合金部品の使用によって前記基材に形成される拡散層の厚さ（拡散層厚）とラーソン・ミラー・パラメータとの関係をそれぞれ求める第１Ｃ工程と、
各試料に関する前記関係に基づいて拡散層厚とラーソン・ミラー・パラメータとの関係を示す第１の近似式を求める第１Ｄ工程と、
前記第１の近似式とラーソン・ミラー・パラメータの定義式との第１の連立方程式に、前記各部位における拡散層厚の測定値と、前記使用時間とを代入する第１Ｅ工程と、
から前記使用温度の推定値を求めることを特徴とする請求項１又は２に記載のニッケル基合金部品の寿命推定方法。
【請求項４】
前記第２工程では、前記第１Ｃ工程から第１Ｅ工程によって求められた前記使用温度の推定値を、前記クリープマスターカーブに当てはめることにより前記最損傷部位の破断時間を求めることを特徴とする請求項３に記載のニッケル基合金部品の寿命推定方法。
【請求項５】
前記第２工程は、異なる試験温度でクリープ試験を行った前記ニッケル基合金部品と同等な複数の試料について、前記試料の応力軸に直交する直交面における複数の結晶粒の結晶方位と結晶粒の結晶格子面に沿って析出した複数の析出物の厚さとをそれぞれ計測すると共にラーソン・ミラー・パラメータを求める第２Ａ工程と、
前記結晶粒の結晶方位と析出物の厚さとに基づいて析出物の真の厚さ（析出物厚）を求める第２Ｂ工程と、
前記複数の試料のラーソン・ミラー・パラメータと、前記結晶方位と前記応力軸とのなす負荷角度の余弦値との関係の直線性が良好な領域における特定の前記負荷角度の前記析出物厚と、に基づいて第２の近似式を求める第２Ｃ工程と、
前記第２の近似式とラーソン・ミラー・パラメータの定義式との第２の連立方程式に、前記最損傷部位から得られる前記特定の負荷角度と同一の負荷角度における析出物厚と、前記ニッケル基合金部品の使用開始から経過した使用時間と、を代入することにより前記最損傷部位における使用温度の推定値を求める第２Ｄ工程と、
このようにして求めた前記最損傷部位の使用温度の推定値を、前記クリープマスターカーブに当てはめることにより前記最損傷部位における破断時間を求める第２Ｅ工程と、
に基づいて前記最損傷部位の破断時間を求めることを特徴とする請求項１から４のうちいずれか一項に記載のニッケル基合金部品の寿命推定方法。
【請求項６】
前記第２工程では、前記最損傷部位から作成した試料に対してクリープ試験を行って破断時間を求めることを特徴とする請求項１から５のうちいずれか一項に記載のニッケル基合金部品の寿命推定方法。
【請求項７】
前記ニッケル基合金部品は、ガスタービンの動翼または静翼であることを特徴とする請求項１から６のうちいずれか一項に記載のニッケル基合金部品の寿命推定方法。
【請求項１】
多結晶ニッケル基合金から形成されたニッケル基合金部品の寿命推定方法であって、
前記ニッケル基合金部品の各部位における使用時の使用温度の推定値と応力の推定値とに基づいて、相対的に損傷の度合いが大きい最損傷部位を特定する第１工程と、
前記最損傷部位について１または複数の手法により破断までの破断時間を求める第２工程と、
を有することを特徴とするニッケル基合金部品の寿命推定方法。
【請求項２】
前記第１工程は、
前記ニッケル基合金部品の使用時の応力軸方向に垂直な断面において前記使用温度の推定値が最も高い垂直断面最損傷部位を特定する第１Ａ工程と、
前記垂直断面最損傷部位の応力軸方向に位置する前記各部位のうち、前記使用温度の推定値と、前記応力の推定値と、前記ニッケル基合金部品の使用開始から経過した使用時間と、前記ニッケル基合金部品の同等部品について予め取得されたクリープマスターカーブとに基づいて、前記最損傷部位を特定する第１Ｂ工程と、
からなることを特徴とする請求項１に記載のニッケル基合金部品の寿命推定方法。
【請求項３】
前記ニッケル基合金部品が前記ニッケル基合金からなる基材の表面に金属コーティング層が形成されたものである場合、
前記第１工程では、
前記ニッケル基合金部品と同等な複数の試料について、前記ニッケル基合金部品の使用によって前記基材に形成される拡散層の厚さ（拡散層厚）とラーソン・ミラー・パラメータとの関係をそれぞれ求める第１Ｃ工程と、
各試料に関する前記関係に基づいて拡散層厚とラーソン・ミラー・パラメータとの関係を示す第１の近似式を求める第１Ｄ工程と、
前記第１の近似式とラーソン・ミラー・パラメータの定義式との第１の連立方程式に、前記各部位における拡散層厚の測定値と、前記使用時間とを代入する第１Ｅ工程と、
から前記使用温度の推定値を求めることを特徴とする請求項１又は２に記載のニッケル基合金部品の寿命推定方法。
【請求項４】
前記第２工程では、前記第１Ｃ工程から第１Ｅ工程によって求められた前記使用温度の推定値を、前記クリープマスターカーブに当てはめることにより前記最損傷部位の破断時間を求めることを特徴とする請求項３に記載のニッケル基合金部品の寿命推定方法。
【請求項５】
前記第２工程は、異なる試験温度でクリープ試験を行った前記ニッケル基合金部品と同等な複数の試料について、前記試料の応力軸に直交する直交面における複数の結晶粒の結晶方位と結晶粒の結晶格子面に沿って析出した複数の析出物の厚さとをそれぞれ計測すると共にラーソン・ミラー・パラメータを求める第２Ａ工程と、
前記結晶粒の結晶方位と析出物の厚さとに基づいて析出物の真の厚さ（析出物厚）を求める第２Ｂ工程と、
前記複数の試料のラーソン・ミラー・パラメータと、前記結晶方位と前記応力軸とのなす負荷角度の余弦値との関係の直線性が良好な領域における特定の前記負荷角度の前記析出物厚と、に基づいて第２の近似式を求める第２Ｃ工程と、
前記第２の近似式とラーソン・ミラー・パラメータの定義式との第２の連立方程式に、前記最損傷部位から得られる前記特定の負荷角度と同一の負荷角度における析出物厚と、前記ニッケル基合金部品の使用開始から経過した使用時間と、を代入することにより前記最損傷部位における使用温度の推定値を求める第２Ｄ工程と、
このようにして求めた前記最損傷部位の使用温度の推定値を、前記クリープマスターカーブに当てはめることにより前記最損傷部位における破断時間を求める第２Ｅ工程と、
に基づいて前記最損傷部位の破断時間を求めることを特徴とする請求項１から４のうちいずれか一項に記載のニッケル基合金部品の寿命推定方法。
【請求項６】
前記第２工程では、前記最損傷部位から作成した試料に対してクリープ試験を行って破断時間を求めることを特徴とする請求項１から５のうちいずれか一項に記載のニッケル基合金部品の寿命推定方法。
【請求項７】
前記ニッケル基合金部品は、ガスタービンの動翼または静翼であることを特徴とする請求項１から６のうちいずれか一項に記載のニッケル基合金部品の寿命推定方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２００９−７４８６８（Ｐ２００９−７４８６８Ａ）
【公開日】平成２１年４月９日（２００９．４．９）
【国際特許分類】
【出願番号】特願２００７−２４２９６７（Ｐ２００７−２４２９６７）
【出願日】平成１９年９月１９日（２００７．９．１９）
【出願人】（０００００００９９）株式会社ＩＨＩ (5,014)
【Ｆターム（参考）】
【公開日】平成２１年４月９日（２００９．４．９）
【国際特許分類】
【出願日】平成１９年９月１９日（２００７．９．１９）
【出願人】（０００００００９９）株式会社ＩＨＩ (5,014)
【Ｆターム（参考）】
[ Back to top ]