エアフォイルの余寿命を判定する方法400であって、処理システムで前記エアフォイルに関連する少なくとも1つの振動応答パラメータを受信するステップ402と、亀裂を生じた少なくとも1つのエアフォイルの有限要素モデルを処理するステップであって、前記亀裂を生じたエアフォイルの有限要素モデルにアクセスするステップ404と、破壊力学に基づく有限要素解析を用いて、前記亀裂を生じたエアフォイルの有限要素モデルのモーダル応力強度係数SIFを計算するステップ406と、少なくとも部分的に、前記モーダルSIF及び少なくとも1つの振動応答パラメータに基づいて振動SIFを計算するステップ408とを含む、亀裂を生じた少なくとも1つのエアフォイルの有限要素モデルを処理するステップと、少なくとも部分的に、前記振動SIFに基づいて前記エアフォイルの余寿命インジケータを計算するステップ410とを含む方法。前記少なくとも1つの振動応答パラメータが、エアフォイルの加振周波数、エアフォイルの振動振幅、及びエアフォイルの構造減衰のうちの少なくとも1つを含む、請求項1に記載の方法。前記余寿命インジケータを計算するステップが、少なくとも1つの材料特性パラメータを考慮にいれるステップを更に含む、請求項1に記載の方法。前記余寿命インジケータが、亀裂の成長速度、サイクルカウント判定、及び亀裂のサイズ対サイクルカウント判定のうちの少なくとも1つを含む、請求項1に記載の方法。エアフォイルの余寿命を判定するシステム108であって、前記エアフォイルに関連する少なくとも1つの振動応答パラメータを受信する402受信機モジュールと、亀裂を生じた少なくとも1つのエアフォイルの有限要素モデルを処理する有限要素解析エンジン206とを備え、前記亀裂を生じたエアフォイルの有限要素モデルを処理するステップが、前記亀裂を生じたエアフォイルの有限要素モデルにアクセスするステップ404と、破壊力学に基づく有限要素解析を用いて、前記亀裂を生じたエアフォイルの有限要素モデルのモーダル応力強度係数SIFを計算するステップ406と、少なくとも部分的に、前記モーダルSIF及び少なくとも1つの振動応答パラメータに基づいて振動SIFを計算するステップ408とを含み、且つ少なくとも部分的に、前記振動SIFに基づいて前記エアフォイルの余寿命インジケータを計算する410判定モジュール208を備えるシステム。前記少なくとも1つの振動応答パラメータが、エアフォイルの加振周波数、エアフォイルの振動振幅、及びエアフォイルの構造減衰のうちの少なくとも1つを含む、請求項5に記載のシステム。前記余寿命インジケータが、亀裂の成長速度、サイクルカウント判定、及び亀裂のサイズ対サイクルカウント判定のうちの少なくとも1つを含む、請求項5に記載のシステム。エアフォイルの余寿命を判定するために、コンピュータ実行可能命令でエンコードされた非一時的コンピュータ読取り可能媒体を備えるコンピュータプログラム製品であって、前記コンピュータ実行可能命令は、実行されると、1つ又は複数のプロセッサが、前記エアフォイルに関連する少なくとも1つの振動応答パラメータを受信し402、亀裂を生じた少なくとも1つのエアフォイルの有限要素モデルを処理し、前記亀裂を生じたエアフォイルの有限要素モデルを処理するステップが、前記亀裂を生じたエアフォイルの有限要素モデルにアクセスするステップ404と、破壊力学に基づく有限要素解析を用いて前記亀裂を生じたエアフォイルの有限要素モデルのモーダル応力強度係数SIFを計算するステップ406と、少なくとも部分的に、前記モーダルSIF及び少なくとも1つの振動応答パラメータに基づいて振動SIFを計算するステップ408とを含む、前記亀裂を生じたエアフォイルの有限要素モデルを処理するステップと、少なくとも部分的に、前記振動SIFに基づいて前記エアフォイルの余寿命インジケータを計算するステップ410とを含む、コンピュータプログラム製品。前記少なくとも1つの振動応答パラメータが、エアフォイルの加振周波数、エアフォイルの振動振幅、及びエアフォイルの構造減衰のうちの少なくとも1つを含む、請求項8に記載のコンピュータプログラム製品。前記余寿命インジケータが、亀裂の成長速度、サイクルカウント判定、及び亀裂のサイズ対サイクルカウント判定のうちの少なくとも1つを含む、請求項8に記載のコンピュータプログラム製品。従って、既知の解決策に関連する上記及びその他の欠点を克服する方法及びシステムが必要とされる。本発明の実施形態は、共振条件又はそれに近い条件で、振動レベルの上昇による高いサイクル疲労破損を受けるリスクにあるエアフォイルの余寿命を判定するシステム及び方法を提供する。一実施形態によれば、エアフォイルの余寿命を判定する方法は、エアフォイルに関連する少なくとも1つの振動応答パラメータを処理システムで受信するステップを含む。この方法は、亀裂を生じた少なくとも1つのエアフォイルの有限要素モデルを処理する。亀裂を生じたエアフォイルの有限要素モデルの処理プロセスは、亀裂を生じたエアフォイルの有限要素モデルにアクセスし、破壊力学に基づく有限要素解析を用いて亀裂を生じたエアフォイルの有限要素モデルのモーダル応力強度係数SIFを計算し、少なくとも部分的にモーダルSIFに基づき、且つ少なくとも1つの振動応答パラメータに基づいて振動SIFを計算する。一実施形態によれば、エアフォイルの余寿命を判定するシステムは、エアフォイルに関連する少なくとも1つの振動応答パラメータを受信するための受信機モジュールを含む。システムは更に、亀裂を生じた少なくとも1つのエアフォイルの有限要素モデルを処理するための有限要素解析エンジンを含む。有限要素解析エンジンは、亀裂を生じたエアフォイルの有限要素モデルにアクセスすることによって、亀裂を生じたエアフォイルの有限要素モデルを処理し、破壊力学に基づく有限要素解析を用いて亀裂を生じたエアフォイルの有限要素モデルのモーダル応力強度係数SIFを計算し、少なくとも部分的にモーダルSIFに基づき、且つ少なくとも1つの振動応答パラメータに基づいて振動SIFを計算する。システムは、少なくとも部分的に振動SIFに基づいてエアフォイルの余寿命インジケータを計算するための判定モジュールを含む。別の実施形態によれば、エアフォイルの余寿命を判定するためのコンピュータ実行可能命令でエンコードされた非一時的コンピュータ読取可能媒体を備えるコンピュータプログラム製品が提供される。コンピュータ実行可能命令は実行されると、1つ以上のプロセッサがエアフォイルに関連する少なくとも1つの振動応答パラメータを受信するようにさせる。コンピュータ実行可能命令は更に、1つ以上のプロセッサが亀裂を生じた少なくとも1つのエアフォイルの有限要素モデルを処理するようにさせる。処理プロセスには亀裂を生じたエアフォイルの有限要素モデルにアクセスし、破壊力学に基づく有限要素解析を用いて亀裂を生じたエアフォイルの有限要素モデルのモーダル応力強度係数SIFを計算し、少なくとも部分的にモーダルSIFに基づき、且つ少なくとも1つの振動応答パラメータに基づいて振動SIFを計算するステップが含まれる。次いでコンピュータ実行可能命令は、1つ以上のプロセッサが少なくとも部分的に振動SIFに基づいてエアフォイルの余寿命インジケータを計算するようにさせる。上記及びその他の利点並びに特徴は、添付図面を参照した好ましい実施形態の以下の詳細な説明によってより容易に理解されよう。本発明の実施形態は、共振条件又はそれに近い条件で高いサイクル疲労を受けるリスクにあるエアフォイルの余寿命を判定する方法、システム及びコンピュータプログラム製品を提供する。一実施形態では、方法はエアフォイルの余寿命を計算するために物理学に基づく三次元有限要素解析を利用する。図1は、様々な実施形態が動作し得る例示的環境100を示す。環境100はターボ機械102、ブレード監視システム104、フィールドデータストア106、及び余寿命判定システム108を含む。ターボ機械102は、ガスタービンエンジン、スチームタービンなどを含み得るが、それらに限定されない。ブレード監視システム104は、圧縮機のエアフォイルなどのターボ機械のエアフォイルの動作状態を監視する。ブレード監視システム104は、例えばブレードの振動振幅、ブレードの振動周波数及びエアフォイルの構造減衰などの振動応答パラメータを監視し得る。ブレード監視システム104は更に、エアフォイルの亀裂データをも監視し得る。亀裂データは、亀裂の位置、亀裂のサイズ、亀裂伝搬方向、又はそれらの組み合わせを含み得る。動作状態及び亀裂データに加えて、ブレード監視システム104は、エアフォイルの共振周波数などのエアフォイルの構造特性をも監視し得る。ブレード監視システム104は、動作状態、亀裂データ及びエアフォイルの構造特性を監視するセンサを含み得る。あるいは、亀裂データ及びエアフォイルの構造特性を試験及び検査システムから取り込んでも良く、又は技術者からの入力として受け取っても良い。その代わりに、又はそれに加えて、亀裂データ及びエアフォイルの構造特性をフィールドデータストア106から受け取る。フィールドデータストア106は、幾つかのターボ機械102のフィールドデータを記憶し得る。フィールドデータストア106は、ターボ機械102の検査から得られた亀裂データを記憶し得る。本発明の一実施形態では、超音波検査技術、又は光学検査技術を用いて亀裂を検知し、亀裂データを取得する。例として特定の検査技術が存在し、本発明の実施形態の範囲を限定するものではないことが理解されよう。フィールドデータストア106は更に、エアフォイルの幾何形状、エアフォイルのブレード寸法、エアフォイル内のブレード数、定格動作速度、定格動作負荷などを含むがこれらに限定されないターボ機械102の設計パラメータをも記憶し得る。一実装形態では、フィールドデータストア106は、ターボ機械102を配備し得る施設にローカルなものである。別の実装形態では、フィールドデータストア106は、ターボ機械102のメーカーの施設などの中央位置にあり、適宜のネットワークを経て余寿命判定システム108に接続される。更に、余寿命判定システム108を、同じ位置又は別の位置に配備された複数のターボ機械102に接続する。あるいは、各ターボ機械102を専用の余寿命判定システム108に接続しても良い。図2は、一実施形態による例示的な余寿命判定システム108のブロック図である。余寿命判定システム108は、立体幾何モジュール202、有限要素モデラ204、有限要素解析エンジン206、及び判定モジュール208を含む。余寿命判定システム108は、1つ以上の余寿命インジケータを提示するために出力デバイス210に結合される。余寿命判定システム108は更に、受信機モジュール図2には図示せずに結合される。受信機モジュールは、亀裂データ、ターボ機械の設計パラメータ、及び動作状態を受信し得る。受信機モジュールは動作状態を受信するためにブレード監視システム104に結合され、且つ亀裂データ及びターボ機械の設計パラメータを受信するためにフィールドデータストア106に結合される。受信機モジュールは、UART、USART、シリアルポート、コントローラエリアネットワークバス、IEEE802.15イーサネット商標インターフェース、ワイヤレステレメトリネットワークなどを含むがこれらに限定されない適宜のインターフェースを用いてブレード監視システム104及びフィールドデータストア106と通信し得る。立体幾何モジュール202は、ターボ機械102のエアフォイルの立体幾何モデルを生成する。立体幾何モジュール202は、エアフォイルの形状、エアフォイルの寸法及びターボ機械の圧縮機又はタービン内のエアフォイルの数などのターボ機械102の設計パラメータを受け入れる。次いで立体幾何モジュール202は、設計パラメータに基づいてエアフォイルの立体幾何モデルを生成する。設計パラメータに加えて、立体幾何モジュール202は、亀裂を生じたエアフォイルの立体幾何モデルを生成するために亀裂データを受け入れる。立体幾何モジュール202は、亀裂のサイズや亀裂伝搬方向などの異なる構造的損傷レベルの亀裂を生じた幾つかのエアフォイルについて立体幾何モデルを生成する。亀裂を生じたエアフォイルの立体幾何モデルは、亀裂を生じた現場の個別のエアフォイルに対応する。次いで、立体幾何モジュール202は立体幾何モデルを有限要素モデラ204に送る。有限要素モデラ204は、立体幾何モデルをメッシュ分割して有限要素モデルを生成する。有限要素モデラ204は、例えば四面体要素などの多面体要素を用いて有限要素モデルを生成する。しかし、それらに限定されないが六面体、八面体、十二面体、及び二十面体要素などのその他の適宜の多面体要素を用いて有限要素モデルを生成しても良い。メッシュ密度は、余寿命判定システム108の必要な、且つ利用可能な計算能力に応じて選択される。一実装形態では、有限要素モデラ204は、ターボ機械の圧縮機又はタービンの単一のエアフォイルの有限要素モデルを生成する。ロータスロットと接触する接合面での有限要素モデルのノードに非線形接触動力学を適用する。このような非線形接触動力学は、高精細解像度で物理を捕捉し、平均SIFを正確に計算し得る。あるいは、エアフォイルの接合面がいかなるマイクロスリップ運動もなくロータスロットと完全に接触するという前提の下で、ロータスロットと接触する接合面のノードに固定境界条件を適用しても良い。このような前提は、高精細度を犠牲にすれば計算が簡単で単純な有限要素モデルを可能にする。次いで有限要素モデラ204は、有限要素モデルを有限要素解析エンジン206に転送し得る。次いで有限要素解析エンジン206は、破壊力学に基づくアルゴリズムを用いて応力強度係数SIFを計算する。有限要素解析エンジン206は、静的応答パラメータについて平均SIFを計算し、且つ振動応答モードについてモーダルSIF言い換えるとピーク交番SIFを計算する。有限要素解析エンジン206は、J積分法、M積分法、又は変位外挿法などの破壊力学を用いて平均SIF及びモーダルSIFを計算する。本明細書には破壊力学に基づくアルゴリズムの特定例が記載されるが、有限要素解析エンジン206は、既知のいずれの破壊力学に基づくアルゴリズムを用いて平均SIF及びモーダルSIFを計算しても良い。次いで有限要素解析エンジン206は、共振周波数及びそれに近い周波数での高サイクル条件での振動に起因する応力を解明するために、振動応答パラメータを考慮に入れる。有限要素解析エンジン206は、モーダルSIF及び振動応答パラメータに基づいて振動SIFを計算する。一実装形態では、有限要素解析エンジン206は、一定間隔で計測されたエアフォイルの先端振動振幅及びエアフォイルの構造減衰を用いて、モーダルSIFを共振周波数における振動SIFに変換する。次いで有限要素解析エンジン206は、周波数離調パラメータを考慮に入れつつ、適宜の一自由度系の強制応答方程式を用いて近共振周波数における振動SIFを計算する。周波数離調パラメータを計算するため、有限要素解析エンジン206はモーダル解析を行って、損傷していないエアフォイルの共振速度、及び亀裂を生じたエアフォイルの共振速度を特定する。亀裂を生じたエアフォイルの共振周波数を計算する方程式の例は以下の通りである。共振速度対亀裂サイズの曲線によって周波数離調パラメータが得られる。次いで有限要素解析エンジン206は、亀裂を生じたエアフォイルの振動振幅を計算し得る。亀裂を生じたエアフォイルの振動振幅を計算する方程式の例は以下の通りである。亀裂を生じたエアフォイルの振動振幅x’0及び計測されたエアフォイルの振動振幅を用いて以下のようなスケーリング係数を計算し得る。次いで有限要素解析エンジン206は、モーダルSIFをスケーリング係数と乗算することによって振動SIFを計算し得る。エアフォイルの振動振幅、及びエアフォイルの応力応答は、例えばエアフォイルの共振周波数、エアフォイルの構造減衰、加振力の振幅、及び振動モードの節直径などの要因によって決まる。エアフォイルの振動振幅に影響を及ぼす別の要因は、ノズルの非対称性、伴流の形状変化、漏流、及び継手などのターボ機械構造内の乱流などの要因に起因するターボ機械を通る流体の不均一な流れである。多数の要因によってエアフォイルの振動振幅の正確な評価が非常に困難になり、複雑な計算上の流体力学や強制応答解析を要する場合が多い。このような技術には膨大な計算量と時間を要する。測定された振動応答パラメータ及びモーダルSIFに基づいて振動SIFを計算することで、これらの欠点が克服される。次いで有限要素解析エンジン206は、平均SIF及び振動SIFに基づいて有効SIFを計算する。平均SIFはターボ機械のエアフォイルのへの静的負荷を解明し、振動SIFはターボ機械のエアフォイルへの動的負荷を解明する。有限要素解析エンジン206は、従来の物理学に基づく破壊力学技術を用いて有効SIFを計算する。判定モジュール208は有効SIFの値を受け取り、有効SIFに基づいて亀裂成長速度を計算する。判定モジュール208は、振動応力率を考慮に入れて△Keffを計算するためにウォーカーの方程式を用いる。ウォーカーの方程式は、負荷率Rを用いて平均応力効果を組み入れるモデルである。次いで判定モジュール208は、パリスの法則を用いて亀裂成長速度を計算する。パリスの方程式の一般形式は以下の通りである。判定モジュール208は更に、亀裂成長速度のシグモイド関数表現を用いても良い。シグモイド関数表現は亀裂成長速度に対するSIFをプロットし、エアフォイルが故障する前に、亀裂が伝搬し得ないSIFの範囲、亀裂がゆっくりと伝搬するSIFの範囲、及び亀裂が急激に伝搬するSIFの範囲を特定する。次いで判定モジュール208は、亀裂のサイズ対サイクルカウントの曲線を得るために、様々なレベルの構造的損傷について亀裂の成長速度を数値積分する。亀裂のサイズ対サイクルカウントの曲線は、ターボ機械のエアフォイルが破局的な損傷を受ける前に、所定の亀裂サイズについて監視された動的負荷条件の下でのターボ機械のエアフォイルの余寿命を示す。更に、エアフォイルが破局的な損傷を受ける時期を判定するために、破壊力学方法及び幾何的方法を導入しても良い。1つの破壊力学的方法は、亀裂を生じたエアフォイルのKmaxがエアフォイルの材料の臨界SIF以上である場合に、亀裂を生じたエアフォイルが急激な亀裂の成長を受け、妨げとなる破局的な故障を示すという事実に基づくものである。1つの幾何的方法は、亀裂がエアフォイルの亀裂のないリガメントを通った場合に、妨げとなる破局的な損傷を予測する。破壊力学的方法と幾何的方法のいずれか一方又は両方を用いて、妨げとなる破局的な故障を特定しても良い。実際に故障するかなり前に、ターボ機械のエアフォイルの妨げとなる故障を特定するための、慎重な余寿命インジケータを提供するように判定モジュール208を調整する。次いで出力デバイス210が余寿命インジケータを提示する。出力デバイス210は、適当なディスプレイで技術者に余寿命インジケータを提示する。あるいは、出力デバイス210は、それに限定されないがIEEE802.15イーサネット商標回線、電話回線、無線ネットワークなどの適当なネットワークを用いて余寿命インジケータを中央監視ステーションに送信しても良い。一実施形態では、本明細書に記載の個々のモジュールは、非一時的コンピュータ読取り可能媒体に記憶されたコンピュータ実行可能命令として実装される。非一時的コンピュータ読み取り可能媒体は、ランダムアクセスメモリRAM、読み出し専用メモリROM、フラッシュメモリ、又はいずれかの適当な記憶装置である。コンピュータ読取り可能媒体は、余寿命判定システム108の全体的な動作を制御するためのオペレーティングシステム及びその他のソフトウエアを記憶させる。図3は、一実施形態による亀裂を生じたエアフォイルの有限要素モデルを生成するために例示的なプロセス300を示すフローチャートである。ステップ302で、立体幾何モジュール202は、少なくとも1つのエアフォイル設計パラメータを受信する。エアフォイル設計パラメータには、エアフォイルの形状、エアフォイルの寸法、及び例えばターボ機械の圧縮機又はタービン内のエアフォイルの数などのパラメータは含まれる。ステップ304で、立体幾何モジュール202は、亀裂開始位置、亀裂伝搬方向、亀裂のサイズ及び亀裂の形状を含むがこれらに限定されない少なくとも1つの亀裂パラメータを受信する。一実装形態では、立体幾何モジュール202はフィールドデータストア106から亀裂パラメータを受信する。一実施形態では、亀裂パラメータには、亀裂を生じたエアフォイルの検査による実際のフィールドデータが含まれる。ステップ306で、立体幾何モジュール202は、少なくとも1つのエアフォイル設計パラメータに基づいて損傷していないエアフォイルの立体幾何モデルを生成する。立体幾何モジュール202は、自動化された立体幾何モデル化アルゴリズムを用いてエアフォイルの立体幾何モデルを生成する。あるいは、フィールドデータストア106が、様々な異なるターボ機械のエアフォイルの立体幾何モデルを記憶していても良い。立体幾何モジュール202は、余寿命が判定されるエアフォイルの立体幾何モデルを直接受信する。次いで立体幾何モジュール202は、立体幾何モデルを有限要素モデラ204に転送する。ステップ308で、有限要素モデラ204は、損傷していないエアフォイルの立体幾何モデルに基づいて、少なくとも部分的に、1つのエアフォイルの有限要素モデルを生成する。有限要素モデラ204は、例えば四面体要素などの適当な多面体要素を用いて立体幾何モデルをメッシュ分割する。余寿命判定システム108が必要とする正確さ及び精度によって、メッシャが使用するメッシュサイズを調整する。ステップ310で、有限要素モデラ204は、亀裂パラメータ及び損傷していないエアフォイルの有限要素モデルに基づいて、亀裂を生じた少なくとも1つのエアフォイルの有限要素モデルを生成する。一実装形態では、有限要素モデラ204は、各々が現場で検査された亀裂を生じた別個のエアフォイルに対応する、幾つかの亀裂を生じたエアフォイルの有限要素モデルを生成する。有限要素モデラ204は、精度を高めるため、高解像度、つまり小さなメッシュサイズで、亀裂領域を再びメッシュ分割する。有限要素解析エンジン206は、亀裂を生じたエアフォイルの有限要素モデルを処理して、平均SIF、モーダルSIF、振動SIF、及び有効SIFを計算する。有効SIFは亀裂の成長速度を計算し、その後、エアフォイルの余寿命を計算するために用いられる。有限要素解析エンジン206は更に、損傷していないエアフォイルの有限要素モデルを処理して、損傷していないエアフォイルの共振周波数を特定する。図4は、一実施形態によるエアフォイルの余寿命を判定する例示的プロセス400を示すフローチャートである。ステップ402で、有限要素解析エンジン206は、エアフォイルに関連する振動応答パラメータを受信する。振動応答パラメータには、それらに限定されないが、エアフォイルの加振周波数、エアフォイルの振動振幅、及びエアフォイルの構造減衰が含まれる。次いで有限要素解析エンジン206は、現場で検査された亀裂を生じた個別のエアフォイルに対応する、亀裂を生じたエアフォイルの各々の有限要素モデルを処理する。ステップ404~408は、有限要素解析エンジン206が亀裂を生じた各々のエアフォイルの有限要素モデルについて処理する処理ステップである。ステップ404で、有限要素解析エンジン206は、亀裂を生じたエアフォイルの有限要素モデルにアクセスする。有限要素解析エンジン206は、有限要素モデラ204から亀裂を生じたエアフォイルの有限要素モデルにアクセスする。あるいは、有限要素解析エンジン206は、有限要素モデラ204によって生成された亀裂を生じたエアフォイルの有限要素モデルが記憶されているデータストアから亀裂を生じたエアフォイルの有限要素モデルにアクセスしても良い。ステップ406で、有限要素解析エンジン206は、破壊力学技術に基づく物理学を用いた有限要素解析により、亀裂を生じたエアフォイルの有限要素モデルのモーダル応力強度係数SIFを計算する。ステップ408で、有限要素解析エンジン206は、モーダルSIF及び振動応答パラメータに基づいて振動SIFを計算する。言い換えると、有限要素解析エンジン206は、共振又は近共振周波数でエアフォイルが受けた振動応答を因数分解して、SIFを計算する。測定された振動応答パラメータを利用することで、普通ならば大幅に大きい計算力を要し、強制応答解析、又は計算上の流体力学などの従来の技術を用いて大幅に長い時間を要するはずの振動SIFの計算が簡潔になり、計算時間が短縮される。この時点で、有限要素解析エンジン206は、亀裂を生じた様々なエアフォイル有限要素モデルに対応する振動SIFの表を有している。ステップ410で、有限要素解析エンジン206は、少なくとも部分的に振動SIFに基づいてエアフォイルの余寿命インジケータを計算する。特定の一実装形態では、有限要素解析エンジン206は、従来の破壊力学技術を用いて振動SIF及び平均SIFに基づいて有効SIFを計算する。次いで有限要素解析エンジン206は、有効SIFに基づいて亀裂の成長速度を計算する。有限要素解析エンジン206は、それに限定されないがターボ機械の点火開始回数及びターボ機械の点火時間などのミッションサイクルデータを考慮に入れつつ亀裂の成長速度を数値積分する。次いで有限要素解析エンジン206は、亀裂のサイズ対ミッションサイクルカウントの曲線に基づいて余寿命インジケータを計算する。余寿命インジケータの計算については図2に関連して詳細に記載されている。本発明の実施形態を特定のハードウエア及びソフトウエアのコンポーネントに関連して記載したが、異なる組み合わせのハードウエア及び/又はソフトウエアのコンポーネントを使用することもでき、且つハードウエアに実装されたものとして記載された特定の動作をソフトウエアに実装しても良く、又はその逆でも良いことが当業者には理解されよう。また、当業者には別の実施形態も明らかであろう。本発明を説明目的のために幾つかの実施形態に関連して説明した。この説明から、本発明は記載の実施形態に限定されるものではなく、添付の請求項の趣旨と範囲によってのみ限定される修正及び変更を伴って実施しても良いことを当業者には理解されよう。