(57)【要約】【課題】解析精度を維持しながら解析の高速化を図ることができるき裂進展解析方法を提供する。【解決手段】き裂進展解析方法は、領域分割ステップ、剛性マトリクス作成ステップ、弾塑性解析ステップ、および弾性解析ステップを有する。領域分割ステップは、弾塑性解析が行われる弾塑性解析領域と、弾性解析が行われる弾性解析領域とに解析モデルの領域を分割する。剛性マトリクス作成ステップは、弾塑性解析領域に弾塑性マトリクスを作成するとともに、弾性解析領域に弾性マトリクスを作成する。弾塑性解析ステップは、弾塑性マトリクスの計算を行う。弾性解析ステップは、弾塑性解析領域における境界の変位に弾性解析領域における境界の変位を一致させた後、弾性解析領域における境界の反力に弾塑性解析領域における境界の反力が一致するまで、弾性マトリクスの計算を繰り返す。【選択図】図１
2016180631.tif 000002
【特許請求の範囲】
【請求項１】
弾塑性解析が行われる弾塑性解析領域と、弾性解析が行われる弾性解析領域とに、解析モデルの領域を分割する領域分割ステップと、
前記弾塑性解析領域に弾塑性マトリクスを作成するとともに、前記弾性解析領域に弾性マトリクスを作成する剛性マトリクス作成ステップと、
前記弾塑性マトリクスの計算を行う弾塑性解析ステップと、
前記弾塑性解析領域における境界の変位に前記弾性解析領域における境界の変位を一致させた後、前記弾塑性解析領域における境界の反力に前記弾性解析領域における境界の反力が一致するまで、前記弾性マトリクスの計算を繰り返して行う弾性解析ステップと、
前記弾塑性解析ステップおよび前記弾性解析ステップにおいて得られた計算結果に基づいて前記解析モデルの形状を更新する形状更新ステップと、
を有することを特徴とするき裂進展解析方法。
【請求項２】
前記領域分割ステップ、前記剛性マトリクス作成ステップ、前記弾塑性解析ステップ、前記弾性解析ステップ、および前記形状更新ステップを有するステップ本体を繰り返して行うことを特徴とする請求項１記載のき裂進展解析方法。
【請求項３】
前記ステップ本体のうち２回目以降に行われるものにおける前記領域分割ステップは、その直前に行われる前記ステップ本体の計算結果から算出されるミーゼスの相当応力に基づいて、前記弾塑性解析領域と前記弾性解析領域とに前記解析モデルの領域を分割することを特徴とする請求項２記載のき裂進展解析方法。
【請求項４】
前記ステップ本体よりも前に初期ステップを有し、
前記初期ステップは、
前記解析モデルの全体に弾塑性マトリクスを作成する初期剛性マトリクス作成ステップと、
前記初期剛性マトリクス作成ステップにおいて作成された前記弾塑性マトリクスの計算を行う初期弾塑性解析ステップと、
を有することを特徴とする請求項２または３記載のき裂進展解析方法。
【請求項５】
前記初期ステップは、さらに、
前記解析モデルの初期形状および前記解析モデルにおけるき裂の初期形状を入力する形状入力ステップと、
前記解析モデルの材料物性を入力する材料物性入力ステップと、
前記解析モデルの境界条件を入力する境界条件入力ステップと、
を有することを特徴とする請求項４記載のき裂進展解析方法。
【請求項６】
前記ステップ本体のうち１回目に行われるものにおける前記領域分割ステップは、前記初期ステップの計算結果から算出されるミーゼスの相当応力に基づいて、前記弾塑性解析領域と前記弾性解析領域とに前記解析モデルの領域を分割することを特徴とする請求項４または５記載のき裂進展解析方法。
【請求項７】
弾塑性解析が行われる弾塑性解析領域と、弾性解析が行われる弾性解析領域とに、解析モデルの領域を分割する領域分割ステップと、
前記弾塑性解析領域に弾塑性マトリクスを作成するとともに、前記弾性解析領域に弾性マトリクスを作成する剛性マトリクス作成ステップと、
前記弾塑性マトリクスの計算を行う弾塑性解析ステップと、
前記弾塑性解析領域における境界の変位に前記弾性解析領域における境界の変位を一致させた後、前記弾塑性解析領域における境界の反力に前記弾性解析領域における境界の反力が一致するまで、前記弾性マトリクスの計算を繰り返して行う弾性解析ステップと、
前記解析モデルの形状を更新する形状更新ステップと、
を有することを特徴とするき裂進展解析プログラム。
【請求項８】
き裂の形状を測定する測定装置と、
請求項７記載のき裂進展解析プログラムを有する計算機と、
を有するき裂進展解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、き裂進展解析方法、き裂進展解析プログラムおよびき裂進展解析装置に関する。
【背景技術】
【０００２】
有限要素法により原子力プラント機器などの大型構造物におけるき裂進展解析を行う場合、数メートル規模の解析モデルに対して、数ｍｍのき裂を仮定して解析を行うことがある。このような場合、き裂の周辺のみが弾塑性変形を示し、その他の大部分については弾性変形を示す。現状、一部のみが弾塑性変形を示すにもかかわらず、全体として複雑な弾塑性解析を行っており、解析時間が膨大となっている。解析時間を短くするために解析モデルの要素を粗くすることが考えられるが、解析モデルの要素を粗くすると解析精度が低下する。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０００−３３９２９７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明が解決しようとする課題は、解析精度を維持しながら解析の高速化を図ることができるき裂進展解析方法を提供することである。
【課題を解決するための手段】
【０００５】
実施形態のき裂進展解析方法は、弾塑性解析が行われる弾塑性解析領域と、弾性解析が行われる弾性解析領域とに、解析モデルの領域を分割する領域分割ステップと、前記弾塑性解析領域に弾塑性マトリクスを作成するとともに、前記弾性解析領域に弾性マトリクスを作成する剛性マトリクス作成ステップと、前記弾塑性マトリクスの計算を行う弾塑性解析ステップと、前記弾塑性解析領域における境界の変位に前記弾性解析領域における境界の変位を一致させた後、前記弾塑性解析領域における境界の反力に前記弾性解析領域における境界の反力が一致するまで、前記弾性マトリクスの計算を繰り返して行う弾性解析ステップと、前記弾塑性解析ステップおよび前記弾性解析ステップにおいて得られた計算結果に基づいて前記解析モデルの形状を更新する形状更新ステップと、を有することを特徴とする。
【発明の効果】
【０００６】
実施形態のき裂進展解析方法、き裂進展解析プログラムおよびき裂進展解析装置は、弾性解析が行われる弾性解析領域と、弾塑性解析が行われる弾塑性解析領域とに解析モデルの領域を分割する。これにより、解析精度を維持しながら解析の高速化を図ることができる。
【図面の簡単な説明】
【０００７】
【図１】実施形態のき裂進展解析方法を示したフローチャートである。
【図２】実施形態のき裂進展解析方法の解析対象となる解析モデルの一例である。
【図３】実施形態のき裂進展解析装置を示す構成図である。
【発明を実施するための形態】
【０００８】
以下、本発明を実施するための形態について説明する。
【０００９】
（第１の実施形態）
図１は、実施形態のき裂進展解析方法を示したフローチャートである。また、図２は、実施形態のき裂進展解析方法の解析対象となる解析モデルの一例である。
【００１０】
実施形態のき裂進展解析方法は、き裂の進展を解析するために使用されるものであり、原子力プラント機器などの大型構造物におけるき裂の進展を解析するために好適に使用されるものである。
【００１１】
まず、図２に示される解析モデル１０について説明する。
解析モデル１０は、実施形態のき裂進展解析方法の解析対象となる試験片の解析モデルである。解析モデル１０は、上下対称であり、上下方向の中心から上側となる第１の部分１１と、下側となる第２の部分１２とを有する。なお、図２では、第１の部分１１のみを実線で示し、かつ剛性マトリクスを付与している。
【００１２】
第１の部分１１と第２の部分１２とは、一方の端部側（図中、右側）が接続部１３により接続されており、他方の端部側（図中、左側）に空間部１４が設けられている。第１の部分１１、第２の部分１２には、それぞれ、解析モデル１０に荷重を加えるときに使用される第１の孔部１１１、第２の孔部１２１が設けられている。
【００１３】
荷重は、第１の孔部１１１と第２の孔部１２１とが互いに離間する方向に印加される。すなわち、第１の孔部１１１には、図中上側方向に荷重が印加され、第２の孔部１２１には、図中下側方向に荷重が印加される。これにより、接続部１３における空間部１４側の部分にき裂が発生し、このき裂が空間部１４とは反対側に向かって進展する。
【００１４】
このような解析モデル１０については、き裂の周辺の領域、例えば、図中、黒枠内の領域が弾塑性解析の対象となる弾塑性解析領域２１となる。また、それ以外の領域、例えば、図中、黒枠外の領域が弾性解析の対象となる弾性解析領域２２となる。通常、弾塑性解析領域２１、弾性解析領域２２は、解析の進行と共にその領域の位置や大きさが変化する。
【００１５】
次に、図１に示される各ステップについて説明する。
実施形態のき裂進展解析方法は、初期ステップＳ１０と、ステップ本体Ｓ２０と、出力ステップＳ３０とを有する。初期ステップＳ１０は、き裂進展解析の最初に行われるステップ群である。ステップ本体Ｓ２０は、初期ステップＳ１０に続けて行われるステップ群である。出力ステップＳ３０は、ステップ本体Ｓ２０に続けて行われる。通常、ステップ本体Ｓ２０は２回以上繰り返して行われ、これ以外の初期ステップＳ１０および出力ステップＳ３０は１回のみ行われる。
【００１６】
まず、初期ステップＳ１０について説明する。
初期ステップＳ１０は、形状入力ステップＳ１１、材料物性入力ステップＳ１２、境界条件入力ステップＳ１３、初期剛性マトリクス作成ステップＳ１４、および初期弾塑性解析ステップＳ１５を有する。以下、初期ステップＳ１０の各ステップについて具体的に説明する。
【００１７】
形状入力ステップＳ１１では、解析モデル１０の初期形状、当該解析モデル１０におけるき裂の初期形状および初期位置、当該解析モデル１０における要素の形状および分割数などを入力する。具体的には、解析対象となる構造物の初期形状、当該構造物におけるき裂の初期形状および初期位置、当該構造物における要素の形状および分割数などを入力する。構造物としては、原子力プラント機器などの大型構造物が好適なものとして挙げられる。要素の形状として、四面体一次要素、四面体二次要素、六面体一次要素、六面体二次要素などが挙げられる。
【００１８】
材料物性入力ステップＳ１２では、解析モデル１０の材料物性、すなわち解析対象となる構造物の材料物性を入力する。材料物性として、ヤング率、耐力、加工硬化指数、ポアソン比などが挙げられる。
【００１９】
境界条件入力ステップＳ１３では、解析モデル１０の境界条件、すなわち解析対象となる構造物の境界条件を入力する。境界条件として、力、変位、応力、およびひずみなどが挙げられる。
【００２０】
初期剛性マトリクス作成ステップＳ１４では、解析モデル１０の全体に弾塑性マトリクスを作成する。
【００２１】
初期弾塑性解析ステップＳ１５では、初期剛性マトリクス作成ステップＳ１４において作成された弾塑性マトリクスの計算（弾塑性解析）を行う。この計算（弾塑性解析）では、弾塑性マトリクスにおける力、変位、応力、およびひずみなどを計算する。
【００２２】
次に、ステップ本体Ｓ２０について説明する。
ステップ本体Ｓ２０は、領域分割ステップＳ２１、剛性マトリクス作成ステップＳ２２、弾塑性解析ステップＳ２３、弾性解析ステップＳ２４、形状更新ステップＳ２５、回数計測ステップＳ２６を有する。通常、ステップ本体Ｓ２０は、２回以上繰り返して行われる。以下、ステップ本体Ｓ２０の各ステップについて具体的に説明する。
【００２３】
領域分割ステップＳ２１では、弾塑性解析が行われる弾塑性解析領域２１と、弾性解析が行われる弾性解析領域２２とに、解析モデル１０の領域を分割する。繰り返して行われるステップ本体Ｓ２０のうち１回目に行われるステップ本体Ｓ２０における領域分割ステップＳ２１は、例えば、初期ステップＳ１０における計算結果と、材料物性入力ステップＳ１２で入力された耐力とを利用して行われる。
【００２４】
具体的には、１回目に行われるステップ本体Ｓ２０における領域分割ステップＳ２１については、初期ステップＳ１０における計算結果から、下記式（１）によりミーゼスの相当応力σVMを算出する。そして、このミーゼスの相当応力σVMと、材料物性入力ステップＳ１２で入力された耐力とを利用して、弾塑性解析領域２１と弾性解析領域２２とに解析モデル１０の領域を分割する。例えば、ミーゼスの相当応力σVMが、材料物性入力ステップＳ１２で入力された耐力の８割以上である節点を含む領域を弾塑性解析領域２１とし、その他の領域を弾性解析領域２２とする。
【００２５】
【数１】
2016180631.tif 000003
…（１）
【００２６】
剛性マトリクス作成ステップＳ２２では、弾塑性解析領域２１に弾塑性マトリクスを作成し、弾性解析領域２２に弾性マトリクスを作成する。
【００２７】
弾塑性解析ステップＳ２３では、弾塑性解析領域２１に作成された弾塑性マトリクスの計算（弾塑性解析）を行う。この計算（弾塑性解析）では、弾塑性マトリクスにおける力、変位、応力、およびひずみなどを計算する。
【００２８】
弾性解析ステップＳ２４では、弾性解析領域２２に作成された弾性マトリクスの計算（弾性解析）を行う。弾性マトリクスの計算（弾性解析）は、以下のようにして行われる。
【００２９】
まず、弾性解析領域２２における境界の変位を、弾塑性解析領域２１における境界の変位に一致させる。すなわち、弾塑性解析領域２１と弾性解析領域２２とでは、解析方法が異なることから、力や変位の計算結果が不連続となる。そこで、弾性解析領域２２における境界の変位として、既に得られた弾塑性解析領域２１における境界の変位を使用する。これにより、弾塑性解析領域２１と弾性解析領域２２との境界における計算結果が連続となり、解析精度が維持される。
【００３０】
変位を一致させた後、弾性解析領域２２における境界の反力が弾塑性解析領域２１における境界の反力に等しくなるまで、弾性解析領域２２に作成された弾性マトリクスについて計算（弾性解析）を繰り返す。この計算（弾性解析）では、弾性マトリクスにおける力、変位、応力、およびひずみなどを計算する。
【００３１】
形状更新ステップＳ２５では、弾塑性解析ステップＳ２３および弾性解析ステップＳ２４において得られた計算結果に基づいて、解析モデル１０の形状を更新する。
【００３２】
形状更新ステップＳ２５の後、再び、領域分割ステップＳ２１、剛性マトリクス作成ステップＳ２２、弾塑性解析ステップＳ２３、弾性解析ステップＳ２４、および形状更新ステップＳ２５を有するステップ本体Ｓ２０が繰り返して行われる。
【００３３】
ここで、２回目以降に行われるステップ本体Ｓ２０、特に領域分割ステップＳ２１については、初期ステップＳ１０における計算結果を利用する代わりに、その直前に行われたステップ本体Ｓ２０の計算結果が利用される。例えば、ｎ回目に行われるステップ本体Ｓ２０における領域分割ステップＳ２１では、ｎ−１回目に行われるステップ本体Ｓ２０における計算結果が利用される。
【００３４】
２回目以降に行われるステップ本体Ｓ２０における領域分割ステップＳ２１についても、ミーゼスの相当応力σVMと、材料物性入力ステップＳ１２で入力された耐力とを利用して、弾塑性解析領域２１と弾性解析領域２２とに解析モデル１０の領域を分割（再分割）する。
【００３５】
具体的には、その直前に行われるステップ本体Ｓ２０における計算結果から、式（１）によりミーゼスの相当応力σVMを算出する。そして、このミーゼスの相当応力σVMと、材料物性入力ステップＳ１２で入力された耐力とを利用して、弾塑性解析領域２１と弾性解析領域２２とに解析モデル１０の領域を分割する。例えば、ミーゼスの相当応力σVMが、材料物性入力ステップＳ１２で入力された耐力の８割以上である節点を含む領域を弾塑性解析領域２１とし、その他の領域を弾性解析領域２２とする。
【００３６】
２回目以降に行われるステップ本体Ｓ２０におけるその他のステップ、すなわち、剛性マトリクス作成ステップＳ２２、弾塑性解析ステップＳ２３、弾性解析ステップＳ２４、形状更新ステップＳ２５は、１回目に行われるステップ本体Ｓ２０における各ステップと同様に行われる。
【００３７】
ステップ本体Ｓ２０の繰り返し回数は、解析対象となる構造物の現在までの使用状況や将来の使用状況などを考慮して適宜設定することができる。例えば、解析対象となる構造物が原子力プラント機器である場合、現在までの運転状況や将来の運転状況などに応じて繰り返し回数を設定する。
【００３８】
繰り返し回数は、各ステップ本体Ｓ２０の最後に行われる回数計測ステップＳ２６により計測される。計測された繰り返し回数が既定値未満の場合、再び、ステップ本体Ｓ２０における一連のステップが行われる。一方、計測された繰り返し回数が既定値に達した場合、出力ステップＳ３０に進んで最終的な計算結果を出力する。計算結果は、例えば、コンター図などとして出力される。計算結果は、例えば、原子力プラント機器などの余寿命の評価に使用される。
【００３９】
以上説明したように、実施形態のき裂進展解析方法によれば、弾塑性解析が行われる弾塑性解析領域と、弾性解析が行われる弾性解析領域とに解析モデル１０の領域を分割する。これにより、き裂の周辺などの弾塑性変形領域についてのみ弾塑性解析が行われ、その他の弾性変形領域については弾性解析が行われる。通常、解析モデル１０の大部分は弾性解析領域となり、かつ弾性解析の計算時間は短いことから、全体として解析の高速化を図ることができる。
【００４０】
また、実施形態のき裂進展解析方法によれば、き裂の周辺などの高精度の計算が必要となる弾塑性変形領域については弾塑性解析により高精度の計算が行われる。また、実施形態のき裂進展解析方法によれば、全体として解析の高速化が図られていることから、従来の方法のように解析の高速化を図るために剛性マトリクスにおける要素を粗くする必要がない。これにより、全体として、解析精度も維持することができる。
【００４１】
（第２の実施形態）
図３は、実施形態のき裂進展解析装置を示す構成図である。
実施形態のき裂進展解析装置４０は、き裂の進展を解析する計算機４１と、き裂の形状を測定する探傷装置４２とを有する。なお、図３には、探傷装置４２により、き裂の形状が測定される測定対象５０を合わせて図示している。測定対象５０としては、原子力プラント機器などが挙げられる。
【００４２】
計算機４１は、例えば、入力部４１１、表示部４１２、および記憶部４１３を有する。入力部４１１は、計算機４１に情報を入力するものであり、例えば、キーボードなどである。表示部４１２は、計算機４１に入力される情報、計算機４１において得られた情報を表示するものであり、例えば、モニタなどである。記憶部４１３は、き裂進展解析プログラムを記憶するものであり、例えば、ＨＤＤなどである。なお、き裂進展解析プログラムは、上述したき裂進展解析方法と同様のステップを有するものである。記憶部４１３は、さらに測定対象５０の情報を含むデータベース４１４を有する。測定対象５０が原子力プラント機器である場合、データベース４１４には、原子力プラント機器の形状、材料物性値、運転状況などの情報が含まれる。
【００４３】
探傷装置４２は、測定対象５０におけるき裂５１の形状を測定するための測定部４２１と、この測定部４２１の制御を行うとともに、測定部４２１からの測定情報が入力される装置本体４２２とを有する。
【００４４】
以下、き裂進展解析装置４０を使用したき裂進展解析について説明する。
まず、探傷装置４２により測定対象５０におけるき裂５１の形状を測定する。測定されたき裂５１の形状は、探傷装置４２から計算機４１に入力され、き裂進展解析プログラムにおけるき裂の初期形状、すなわち、初期ステップＳ１０の形状入力ステップＳ１１におけるき裂の初期形状として利用される。
【００４５】
別途、データベース４１１を利用して、き裂進展解析プログラムの実施に必要な情報を入力する。例えば、データベース４１１から測定対象５０の形状、すなわち解析に使用する解析モデルの形状を選択し、初期ステップＳ１０の形状入力ステップＳ１１における解析モデルの初期形状を入力する。同様に、データベース４１１から測定対象５０の材料物性値、すなわち解析に使用する解析モデルの材料物性値を選択し、初期ステップＳ１０の材料物性入力ステップＳ１２における材料物性値を入力する。また、データベース４１１から測定対象５０の境界条件、すなわち解析に使用する解析モデルの境界条件を選択し、初期ステップＳ１０の境界条件入力ステップＳ１３における境界条件を入力する。
【００４６】
入力後、計算機４１において、初期ステップＳ１０の初期剛性マトリクス作成ステップＳ１４により、入力された情報に基づいて解析モデルの全体に弾塑性マトリクスが作成される。その後、初期ステップＳ１０の初期弾塑性解析ステップＳ１５により、弾塑性マトリクスの計算（弾塑性解析）が行われる。
【００４７】
次いで、ステップ本体Ｓ２０の領域分割ステップＳ２１により、初期弾塑性解析ステップＳ１５により得られた計算結果に基づいて、弾塑性解析領域２１と弾性解析領域２２とに解析モデルの領域が分割される。そして、ステップ本体Ｓ２０の剛性マトリクス作成ステップＳ２２により、弾塑性解析領域２１に弾塑性マトリクスが作成されると共に、弾性解析領域２２に弾性マトリクスが作成される。
【００４８】
さらに、ステップ本体Ｓ２０の弾塑性解析ステップＳ２３により、弾塑性解析領域２１に作成された弾塑性マトリクスの計算（弾塑性解析）が行われる。また、ステップ本体Ｓ２０の弾性解析ステップＳ２４により、弾性解析領域２２に作成された弾性マトリクスの計算（弾性解析）が行われる。そして、ステップ本体Ｓ２０の形状更新ステップＳ２５により、弾塑性解析ステップＳ２３および弾性解析ステップＳ２４により得られた計算結果に基づいて解析モデルの形状が更新される。
【００４９】
ステップ本体Ｓ２０、すなわち、領域分割ステップＳ２１、剛性マトリクス作成ステップＳ２２、弾塑性解析ステップＳ２３、および弾性解析ステップＳ２４は、所定の回数が繰り返して行われる。ステップ本体Ｓ２０の繰り返し回数は、各ステップ本体Ｓ２０の最後に行われる回数計測ステップＳ２６により計測される。所定の繰り返し回数に達した後、出力ステップＳ３０により、最終的な計算結果が出力される。計算結果は、例えば、コンター図などとして表示部４１２に出力される。計算結果は、例えば、原子力プラント機器などの余寿命の評価に使用される。
【００５０】
以上説明したように、実施形態のき裂進展解析装置４０によれば、測定対象５０におけるき裂の形状を反映させたき裂進展解析を行うことができる。また、実施形態のき裂進展解析装置４０によれば、計算機４１と探傷装置４２とからなる簡単な構造のために持ち運びが容易であり、測定対象５０が設けられている現地において解析を行うことができる。また、実施形態のき裂進展解析装置４０によれば、き裂進展解析の計算結果を利用することにより、測定対象５０における余寿命の評価を行うことができる。例えば、原子力プラント機器などの余寿命の評価を行うことができる。
【００５１】
以上、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【００５２】
Ｓ１０…初期ステップ、Ｓ１１…形状入力ステップ、Ｓ１２…材料物性入力ステップ、Ｓ１３…境界条件入力ステップ、Ｓ１４…初期剛性マトリクス作成ステップ、Ｓ１５…初期弾塑性解析ステップ、Ｓ２０…ステップ本体、Ｓ２１…領域分割ステップ、Ｓ２２…剛性マトリクス作成ステップ、Ｓ２３…弾塑性解析ステップ、Ｓ２４…弾性解析ステップ、Ｓ２５…形状更新ステップ、Ｓ２６…回数計測ステップ、Ｓ３０…出力ステップ。
１０…解析モデル、１１…第１の部分、１１１…孔部、１２…第２の部分、１２１…孔部、１３…接続部、１４…空間部、２１…弾塑性解析領域、２２…弾性解析領域、４０…き裂進展解析装置、４１…計算機、４１１…入力部、４１２…表示部、４１３…記憶部、４１４…データベース、４２…探傷装置、４２１…測定部、４２２…装置本体、５０…測定対象、５１…き裂。
【図１】
2016180631.tif 000004
【図２】
2016180631.tif 000005
【図３】
2016180631.tif 000006