(57)【要約】【課題】亀裂をより正確に予測することができる。【解決手段】モデル取得部は、モデル生成部または図示せぬ外部機器などから構造体モデルＭDを取得する。亀裂を有さない要素Ｅ０における振幅負荷エネルギーＡは、当該要素Ｅ０を構成する材料に応じて実験的に得られる相当応力σと相当弾性ひずみεとの関係に基づいて設定される。相当弾性ひずみεは亀裂変数φに依存するため、振幅負荷エネルギーＡは亀裂変数φの関数として表現される。亀裂予測部は、この振幅エネルギーに比例する項を有する微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。本開示は、例えば、亀裂の予測を行う亀裂予測装置に適用することができる。【選択図】図１
2019003301.tif 000002
【特許請求の範囲】
【請求項１】
所定の構造体に対応する構造体モデルを取得するモデル取得部と、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
を備える情報処理装置。
【請求項２】
前記振幅負荷エネルギーは、繰り返し負荷時の応力の最大値と最小値の差分と、歪みの最大値と最小値の差分とが積算されて設定される
請求項１に記載の情報処理装置。
【請求項３】
前記振幅負荷エネルギーに係る定数は、対応する材料のS-N線図の高サイクル疲労領域の特性線の傾きに応じた値として設定される
請求項１に記載の情報処理装置。
【請求項４】
前記振幅負荷エネルギーに係る定数は、対応する材料界面のS-N線図の高サイクル疲労剥離領域の特性線の傾きに応じた値として設定される
請求項１に記載の情報処理装置。
【請求項５】
前記亀裂予測部は、非周期的な繰り返し負荷に対して、応力時間変化の変曲点毎に解析区間を区切り、区切られた解析区間毎に前記構造体に発生する亀裂を予測する
請求項１に記載の情報処理装置。
【請求項６】
前記亀裂予測部は、前記構造体モデルの各位置に設定され、塑性変形時に散逸されるエネルギーを前記亀裂変数を利用して表現する塑性散逸エネルギーに比例する項をさらに含む前記微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
請求項１に記載の情報処理装置。
【請求項７】
前記塑性散逸エネルギーは、相当応力を相当塑性ひずみの微小増分で積分した量を利用して設定される
請求項６に記載の情報処理装置。
【請求項８】
前記塑性散逸エネルギーは、相当応力と降伏応力との差と、相当塑性ひずみとの積を利用して設定され、相当応力が降伏応力より小さい場合にゼロとされる
請求項６に記載の情報処理装置。
【請求項９】
前記微分方程式は、空間座標の２階微分に比例する拡散項をさらに含む
請求項１に記載の情報処理装置。
【請求項１０】
所定の構造体に対応する構造体モデルを取得し、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
情報処理方法。
【請求項１１】
所定の構造体に対応する構造体モデルを取得するモデル取得部と、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
して、コンピュータを機能させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、情報処理装置、情報処理方法、及びプログラムに関し、特に、亀裂をより正確に予測することができる情報処理装置、情報処理方法、及びプログラムに関する。
【背景技術】
【０００２】
一般的に、半導体デバイス等の各種構造体には、その製造過程などにおいて機械的応力等の様々な応力が加わる。構造体にこのような応力が加わると、その内部に亀裂が発生することがある。このような亀裂の発生を防止するために、構造体に予め発生する可能性のある亀裂を予測する技術が用いられる。
【０００３】
特許文献１には、亀裂を含む断面上の異種材料の面積比に応じて各材料の応力とひずみの関係を単一の材料の応力のひずみの関係に置き換えたときの材料定数を用いて求めた破壊力学パラメータＪ積分値が亀裂先端の材料の弾塑性破壊靭性値に等しくなるときの荷重を破壊荷重とする技術が開示されている。
【０００４】
特許文献２には、複数種類の材料の界面を跨ぐ亀裂を予測可能な技術が開示されている。特許文献２に係る技術では、構造体の内部で仮想的に亀裂を進行させたときのエネルギー解放率を算出し、エネルギー解放率の大きい方向に亀裂が進行するものと予測する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００９−１６００２８号公報
【特許文献２】特開２０１１−２０４０８１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
近年、半導体デバイス等の構造体の多様化に伴い、構造体を構成する材料として金属材料や樹脂材料が広く用いられるようになってきている。したがって、金属材料や樹脂材料を用いて構成された構造体に発生する可能性のある亀裂を正確に予測可能な技術が要求されている。
【０００７】
特許文献１に係る技術では、界面をまたぐ亀裂に対応することが難しく、疲労破壊への対応が難しかった。特許文献２に係る技術では、計算負荷が高く、予め定めた特定面以外の面を解析することが困難であった。
【０００８】
本開示は、このような状況に鑑みてなされたものであり、亀裂をより正確に予測することができるものである。
【課題を解決するための手段】
【０００９】
本技術の一側面の情報処理装置は、所定の構造体に対応する構造体モデルを取得するモデル取得部と、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部とを備える。
【００１０】
前記振幅負荷エネルギーは、繰り返し負荷時の応力の最大値と最小値の差分と、歪みの最大値と最小値の差分とが積算されて設定される。
【００１１】
前記振幅負荷エネルギーに係る定数は、対応する材料のS-N線図の高サイクル疲労領域の特性線の傾きに応じた値として設定される。
【００１２】
前記振幅負荷エネルギーに係る定数は、対応する材料界面のS-N線図の高サイクル疲労剥離領域の特性線の傾きに応じた値として設定される。
【００１３】
前記亀裂予測部は、非周期的な繰り返し負荷に対して、応力時間変化の変曲点毎に解析区間を区切り、区切られた解析区間毎に前記構造体に発生する亀裂を予測することができる。
【００１４】
前記亀裂予測部は、前記構造体モデルの各位置に設定され、塑性変形時に散逸されるエネルギーを前記亀裂変数を利用して表現する塑性散逸エネルギーに比例する項をさらに含む前記微分方程式を計算することにより、前記構造体に発生する亀裂を予測することができる。
【００１５】
前記塑性散逸エネルギーは、相当応力を相当塑性ひずみの微小増分で積分した量を利用して設定される。
【００１６】
前記塑性散逸エネルギーは、相当応力と降伏応力との差と、相当塑性ひずみとの積を利用して設定され、相当応力が降伏応力より小さい場合にゼロとされる。
【００１７】
前記微分方程式は、空間座標の２階微分に比例する拡散項をさらに含むことができる。
【００１８】
本技術の一側面の情報処理方法は、所定の構造体に対応する構造体モデルを取得し、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する。
【００１９】
本技術の一側面のプログラムは、所定の構造体に対応する構造体モデルを取得するモデル取得部と、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部として、コンピュータを機能させる。
【００２０】
本技術の一側面においては、所定の構造体に対応する構造体モデルが取得される。そして、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂が予測される。
【発明の効果】
【００２１】
本技術によれば、高サイクル疲労破壊を予測することができる。この結果、亀裂をより正確に予測することができる。
【００２２】
なお、本明細書に記載された効果は、あくまで例示であり、本技術の効果は、本明細書に記載された効果に限定されるものではなく、付加的な効果があってもよい。
【図面の簡単な説明】
【００２３】
【図１】本技術の一実施形態に係る亀裂予測装置（情報処理装置）の構成例を示すブロック図である。
【図２】図１の亀裂予測装置の亀裂予測処理について説明するフローチャートである。
【図３】亀裂予測処理において生成される構造体モデルの一例を示す図である。
【図４】繰り返し負荷に対応した対象材料にかかる応力振幅を定義するための図である。
【図５】対象材料のS-N線図を表す図である。
【図６】材料の低サイクル、高サイクルの疲労特性を記述する際の式を説明する図である。
【図７】対象物に対して想定される負荷が、周期的な振動負荷でない場合の例を説明する図である。
【図８】構造体モデルに加える荷重条件の一例を説明するための図である。
【図９】構造体モデルにおける振幅負荷エネルギーの分布の一例を示す図である。
【図１０】構造体モデルにおける亀裂変数の分布の一例を示す図である。
【図１１】構造体モデルにおいて予測される亀裂の一例を示す図である。
【図１２】弾性率と亀裂変数との関係の一例を示す図である。
【図１３】障壁エネルギーと亀裂変数との関係の一例を示す図である。
【図１４】本技術の他の実施形態に係る亀裂予測装置（情報処理装置）の構成例を示すブロック図である。
【図１５】図１４の亀裂予測装置の亀裂予測処理について説明するフローチャートである。
【図１６】亀裂予測処理において設定される塑性散逸エネルギーの表現方法の一例を示す図である。
【図１７】コンピュータの構成例を示すブロック図である。
【発明を実施するための形態】
【００２４】
以下、本開示を実施するための形態（以下実施の形態とする）について説明する。なお、説明は以下の順序で行う。また、図面には、適宜相互に直交するＸ軸、Ｙ軸、及びＺ軸が示されている。Ｘ軸、Ｙ軸、及びＺ軸は全図において共通である。
０．概要
１．第１の実施形態
２．第２の実施形態
３.コンピュータ
【００２５】
＜０．概要＞
［亀裂予測方法の概要］
本技術に係る亀裂予測方法（情報処理方法）の概要について説明する。本技術に係る亀裂予測方法では、フェーズフィールド（Phase Field）法の考え方を応用して構造体Ｄに発生する亀裂を予測する。まず、本技術に関連するフェーズフィールド法の考え方による亀裂予測方法について説明する。
【００２６】
（フェーズフィールド法の考え方による亀裂予測方法）
構造体Ｄ内のエネルギーＦは、障壁エネルギーｆdoub、勾配エネルギーｆgrad、及び弾性エネルギーｆelastを利用して、式（１）で表される。
【数１】
2019003301.tif 000003
【００２７】
フェーズフィールド法の考え方によると、式（１）から、次の微分方程式（２）を導出することができる。
【数２】
2019003301.tif 000004
【００２８】
微分方程式（２）の左辺は、易動度Ｍの逆数と、亀裂の発生の有無を表現する亀裂変数φの時間微分と、の積で構成されている。微分方程式（２）の右辺は、空間座標の２階微分の拡散項∇（ξ∇φ）と、障壁エネルギーｆdoubの微分項と、弾性エネルギーｆelastの微分項と、で構成されている。微分方程式（２）では、弾性エネルギーｆelastの微分項によって、弾性エネルギーｆelastの解放率が表現されている。
【００２９】
微分方程式（２）の計算を行うために、まず構造体Ｄの各位置にそれぞれ亀裂変数φが設定される。より詳細には、亀裂を有さない位置と、亀裂を有する位置と、で異なる亀裂変数φが設定される。例えば、亀裂を有さない位置の亀裂変数φが「０」と設定され、亀裂を有する位置の亀裂変数φが「１」と設定される。
【００３０】
そして、微分方程式（２）の計算を進めると、時間の経過に伴い、亀裂変数φを「０」と設定した位置のうち、亀裂変数φが「１」以上となる位置が現れる。フェーズフィールド法の考え方による亀裂予測方法では、所定時間が経過した後に、亀裂変数φが「１」以上である位置に亀裂が発生しているものと予測することができる。
【００３１】
フェーズフィールド法の考え方による亀裂予測方法では、微分方程式（２）を計算することにより、亀裂を迅速に予測可能である。また、フェーズフィールド法の考え方による亀裂予測方法では、複数種類の材料の界面を跨ぐ亀裂を予測可能であるため、複数の材料で構成された構造体Ｄに発生する亀裂を予測可能である。更に、フェーズフィールド法の考え方による亀裂予測方法では、亀裂の形状に制約がないため、高い汎用性が得られる。
【００３２】
フェーズフィールド法の考え方による亀裂予測方法では、微分方程式（２）に含まれる弾性エネルギーｆelastの微分項で表現される弾性エネルギーｆelastの解放率によって、脆性破壊による亀裂を予測することができる。しかし、フェーズフィールド法の考え方による亀裂予測方法では、微分方程式（２）に塑性変形に対応する項が含まれていないため、塑性変形を伴う疲労破壊による亀裂を予測することができない。
【００３３】
したがって、フェーズフィールド法の考え方を利用した亀裂予測方法では、金属材料や樹脂材料などの疲労破壊しやすい材料を用いて構成された構造体Ｄに発生する亀裂を正確に予測することは困難である。
【００３４】
そこで、本技術の発明者は、フェーズフィールド法の考え方を応用し、塑性変形時に主に熱として散逸されるエネルギー（以下、「塑性散逸エネルギーｆplast」と称する）を含む項を微分方程式（２）に導入することにより、疲労破壊による亀裂を予測可能となることを見出した。以下、本技術に係るフェーズフィールド法の考え方を応用した亀裂予測方法について説明する。
【００３５】
（フェーズフィールド法の考え方を応用した亀裂予測方法）
本技術に係るフェーズフィールド法の考え方を応用した亀裂予測方法では、微分方程式（２）に塑性散逸エネルギーｆplastの項を導入した、次の微分方程式（３）を利用する。
【数３】
2019003301.tif 000005
【００３６】
微分方程式（３）において、塑性散逸エネルギーｆplastは、弾性エネルギーｆelastとは異なり、微分項ではない。これは、弾性エネルギーｆelastは時間とともに解放されるのに対し、塑性散逸エネルギーｆplastは時間とともに蓄積されるためである。微分方程式（３）では、塑性散逸エネルギーｆplastを微分項としないことにより、塑性散逸エネルギーｆplastの蓄積を表現することができる。
【００３７】
このように、微分方程式（３）には、弾性エネルギーｆelastの解放率を表現する弾性エネルギーｆelastの微分項と、塑性散逸エネルギーｆplastの蓄積を表現する塑性散逸エネルギーｆplastの項と、が含まれる。したがって、微分方程式（３）を計算することにより、脆性破壊及び疲労破壊の双方を考慮して、亀裂を予測することが可能となる。
【００３８】
なお、上述した塑性散逸エネルギーｆplastでは、疲労破壊のうち、低サイクル疲労（塑性疲労）を解析することはできる。しかしながら、高サイクル疲労（弾性疲労）解析について対応することが困難であった。
【００３９】
そこで、本技術の発明者は、さらに、フェーズフィールド法の考え方を応用し、微分方程式（２）に、高サイクルの繰り返し疲労（負荷）を表すエネルギー（以下、振幅負荷エネルギーＡと称する）を含む項を導入することにより、高サイクル疲労（弾性疲労）による亀裂を予測可能となることを見出した。
【００４０】
（フェーズフィールド法の考え方を応用した亀裂予測方法）
本技術に係るフェーズフィールド法の考え方を応用した亀裂予測方法では、微分方程式（２）に、振幅負荷エネルギーＡの項を導入した、次の微分方程式（４）を利用する。
【数４】
2019003301.tif 000006
【００４１】
微分方程式（４）は、微分方程式（３）において、塑性散逸エネルギーｆplastを、振幅負荷エネルギーＡに入れ替えた式である。すなわち、微分方程式（４）は、微分方程式（２）に対して、繰り返し負荷に対応する高サイクル疲労解析を行うための振幅負荷エネルギーＡの項を追加したものである。
【００４２】
微分方程式（４）において、振幅負荷エネルギーＡは、弾性エネルギーｆelastとは異なり、微分項ではない。これは、弾性エネルギーｆelastは時間とともに解放されるのに対し、振幅負荷エネルギーＡは時間とともに蓄積されるためである。微分方程式（４）では、振幅負荷エネルギーＡを微分項としないことにより、振幅負荷エネルギーＡの蓄積を表現することができる。
【００４３】
このように、微分方程式（４）には、弾性エネルギーｆelastの解放率を表現する弾性エネルギーｆelastの微分項と、振幅負荷エネルギーＡの蓄積を表現する振幅負荷エネルギーＡの項と、が含まれる。したがって、微分方程式（４）を計算することにより、脆性破壊及び高サイクル疲労破壊の双方を考慮して、亀裂を予測することが可能となる。
【００４４】
このため、本技術では、金属材料や樹脂材料などの疲労破壊のうち、高サイクル疲労（弾性疲労）しやすい材料を用いて構成された構造体Ｄに発生する亀裂を正確に予測可能である。また、本技術に係る亀裂予測方法でも、フェーズフィールド法の考え方による亀裂予測方法と同様に、複数種類の材料の界面を跨ぐ亀裂を予測可能であるため、複数の材料で構成された構造体Ｄに発生する亀裂を迅速に予測可能である。更に、本技術に係る亀裂予測方法でも、フェーズフィールド法の考え方による亀裂予測方法と同様に、亀裂の形状に制約がないため、高い汎用性が得られる。
【００４５】
＜１．第１の実施形態＞
［亀裂予測装置の詳細］
図１は、本技術の第１の実施形態に係る亀裂予測装置（情報処理装置）１０の構成を示すブロック図である。なお、第１の実施形態においては、疲労破壊のうち弾性疲労（高サイクル疲労）下について述べていく。
【００４６】
亀裂予測装置１０は、モデル生成部１１と、モデル取得部１２と、亀裂変数設定部１３と、振幅負荷エネルギー設定部１４と、微分方程式生成部１５と、亀裂予測部１６と、を備える。
【００４７】
モデル生成部１１は、構造体Ｄの構成を再現したモデル（構造体モデル）ＭDを生成する。モデル取得部１２は、モデル生成部１１により生成された構造体モデルＭDを取得する。
【００４８】
亀裂変数設定部１３は、モデル取得部１２により取得された構造体モデルＭDの各要素Ｅに亀裂の有無を表現する亀裂変数φを設定する。振幅負荷エネルギー設定部１４は、モデル取得部１２により取得された構造体モデルＭDの各要素Ｅに、振幅負荷エネルギーＡを設定する。
【００４９】
微分方程式生成部１５は、亀裂変数設定部１３により設定された亀裂変数φと、振幅負荷エネルギー設定部１４により設定された振幅負荷エネルギーＡとを利用して微分方程式を作成する。
【００５０】
亀裂予測部１６は、微分方程式生成部１５により生成された微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。
【００５１】
［亀裂予測動作の例］
図２のフローチャートは、亀裂予測装置１０の亀裂予測処理を説明するフローチャートである。図３乃至図１３は、図２に示す各ステップを説明するための図である。以下、図２に沿って、図３乃至図１３を適宜参照しながら、第１の実施形態に係る亀裂予測処理について説明する。
【００５２】
（モデル生成ステップ）
ステップＳ１１において、モデル生成部１１は、構造体Ｄの構成を再現したモデル（構造体モデル）ＭDを生成する。構造体モデルＭDによって任意の構造体Ｄの構成を再現可能である。構造体モデルＭDによってその構成を再現可能な構造体Ｄとしては、例えば、半導体デバイスなどの各種デバイスが挙げられる。
【００５３】
亀裂予測方法には、有限要素法（FEM：Finite Element Method）や差分法（FDM：Finite Difference Method）を利用することができる。また、陰解法や陽解法を利用することができる。有限要素法では、任意の形状に対応可能であり、高い汎用性が得られる。差分法では、計算の並列化が容易であり、計算が早いというメリットが得られる。陰解法では、タイムステップを大きくとれるというメリットが得られる。本実施形態では、有限要素法を用いるため、構造体モデルＭDが複数の要素Ｅから構成される。
【００５４】
図３は、ステップＳ１１において生成される構造体モデルＭDを例示する図である。図３Ａは構造体モデルＭDの斜視図であり、図３Ｂは構造体モデルＭDの図３ＡのＡ−Ａ'線に沿った断面図である。図３に示す構造体モデルＭDによってその構成が再現されている構造体Ｄでは、概形が立方体であり、上面の中央にＺ軸方向に延びる初期亀裂が形成されている。
【００５５】
この場合、構造体モデルＭDでは、Ｙ軸方向上面のＸ軸方向中央部にＺ軸方向に並ぶ５つの要素Ｅが亀裂を有する要素Ｅ１とされ、その他の要素Ｅが亀裂を有さない要素Ｅ０とされる。図３には、亀裂を有する要素Ｅ１が斜線で示され、亀裂を有さない要素Ｅ０が白抜きで示されている。なお、空孔などの自由空間を有する要素Ｅも、亀裂を有する要素Ｅ１と同様に扱うことが好ましい。
【００５６】
以下の説明では、一例として図３に示す構造体モデルＭDを用いて説明するが、他の構造体モデルＭDについても同様に扱うことが可能であることは勿論である。
【００５７】
なお、構造体モデルＭDが予め用意されている場合などには、本ステップＳ１１は省略しても構わない。
【００５８】
（モデル取得ステップ）
ステップＳ１２において、モデル取得部１２は、ステップＳ１１で生成された構造体モデルＭDを取得する。なお、ステップＳ１１を行わない場合には、本ステップＳ１２においては、図示せぬ外部機器などから構造体モデルＭDを取得することができる。
【００５９】
（亀裂変数設定ステップ）
ステップＳ１３において、亀裂変数設定部１３は、ステップＳ１２で取得された構造体モデルＭDの各要素Ｅに亀裂の有無を表現する亀裂変数φを設定する。
【００６０】
具体的には、構造体モデルＭDにおいて、亀裂を有さない要素Ｅ０と、亀裂を有する要素Ｅ１と、で異なる亀裂変数φが設定される。つまり、亀裂を有さない要素Ｅ０の亀裂変数φが「ｍ」と設定され、亀裂を有する要素Ｅ１の亀裂変数φが「ｍ」とは異なる「ｎ」と設定される。「ｍ」と「ｎ」とではいずれが大きくても構わない。一例として、亀裂を有さない要素Ｅ０の亀裂変数φが「０」と設定され、亀裂を有する要素Ｅ１の亀裂変数φが「１」と設定される。
【００６１】
なお、構造体モデルＭDに予め亀裂変数φが設定されている場合などには、本ステップＳ１３は省略しても構わない。
【００６２】
（振幅負荷エネルギー設定ステップ）
ステップＳ１４において、振幅負荷エネルギー設定部１４は、ステップＳ１２で取得された構造体モデルＭDの各要素Ｅに、振幅負荷エネルギーＡを設定する。なお、既に亀裂を有する要素Ｅ１では、弾性応答が生じないため、振幅負荷エネルギーＡが蓄積されない。このため、要素Ｅ１の振幅負荷エネルギーＡは「０」と設定される。
【００６３】
亀裂を有さない要素Ｅ０における振幅負荷エネルギーＡは、当該要素Ｅ０を構成する材料に応じて実験的に得られる相当応力σと相当弾性ひずみεとの関係に基づいて設定される。相当弾性ひずみεは亀裂変数φに依存するため、振幅負荷エネルギーＡは亀裂変数φの関数として表現される。
【００６４】
なお、構造体モデルＭDに予め振幅負荷エネルギーＡが設定されている場合などには、本ステップＳ１４は省略しても構わない。
【００６５】
（微分方程式生成ステップ）
ステップＳ１５において、微分方程式生成部１５は、ステップＳ１３で設定した亀裂変数φと、ステップＳ１４で設定した振幅負荷エネルギーＡと、が利用されて微分方程式を作成する。
【００６６】
ステップＳ１５において生成される微分方程式としては、一例として、上記の微分方程式（４）が挙げられる。また、ステップＳ１５においては、微分方程式（４）を改良した、次の微分方程式（５）が生成されてもよい。なお、式（５）は、拡散方程式でもある。
【数５】
2019003301.tif 000007
【００６７】
微分方程式（５）には、障壁エネルギーｆdoubの微分項のフィッティング用定数ｗdoub、弾性エネルギーｆelastの微分項のフィッティング用定数ｗelast、及び振幅負荷エネルギーＡの項のフィッティング用定数ｗhcが導入されている。これにより、構造体Ｄの構成などに応じて、障壁エネルギーｆdoubの微分項と、弾性エネルギーｆelastの微分項と、振幅負荷エネルギーＡの項と、のそれぞれの重み付けを最適化することが可能である。これにより、構造体Ｄに発生する亀裂を更に正確に予測可能となる。
【００６８】
なお、式（５）における振幅負荷エネルギーＡは、繰り返し負荷に対応した対象材料にかかる応力振幅Δσを、図４のグラフで、応力振幅Δσ＝最大応力σmax−最小応力σminと定義したときに、その負荷に対応するエネルギーである。振幅負荷エネルギーＡは、次の式（６）で表される。
【数６】
2019003301.tif 000008
【００６９】
ここで、Δεは、ひずみ振幅であり、最大ひずみεmax−最小ひずみεminである。なお、このひずみは弾性ひずみである。
【００７０】
また、式（５）の振幅負荷エネルギーＡに係るｗhcについて、図５を参照して説明する。
【００７１】
図５は、対象材料のS-N線図を表している。S-N線図は、縦軸に応力振幅Δσ、横軸に破断までの繰り返し数Nをとって疲労試験結果をグラフで表したものである。実際には、繰り返し数Ｎによって、メカニズムの異なる低サイクル疲労（塑性疲労）領域と高サイクル疲労（弾性疲労）領域に分かれており、ｗhcは、高サイクル疲労領域の傾きに対応したパラメータであり、解析上は、亀裂の進行速度を決めるパラメータである。
【００７２】
なお、材料の低サイクル、高サイクルの疲労特性を記述する際には、図６に示される数式である次の式（７）でフィッティングされる場合が多い。
【数７】
2019003301.tif 000009
【００７３】
式（７）の右辺の第１項は、低サイクル疲労領域の応力振幅を表し、第２項が、高サイクル疲労領域の応力振幅を表している。その際、ｗhcは、関数ｆを用いて、次の式（８）で表される。
【数８】
2019003301.tif 000010
【００７４】
なお、ここで、関数ｆは、線形であってもよいし、非線形であってもよい。
【００７５】
異種材料が、物理結合、化学結合、粘着結合、静電結合などで結合しており、その結合力寿命が各材料のバルク内結合力よりも小さい場合、結合部に繰り返し負荷がかかると界面剥離が進行する。その際の特性は、図５および図６で示されたように、界面での疲労破断についても特性線図が得られる。
【００７６】
解析モデルの界面上に対応したｗhcxを定義することで、高サイクルにおける剥離の進行も解析可能となる。
【００７７】
さらに、図７を参照して、対象物に対して想定される負荷が、周期的な振動負荷でない（すなわち、非周期負荷である）場合についての対応を説明する。
【００７８】
非周期負荷の例として、対象にかかる応力の時刻歴が、図７に示されるような場合を考える。図７の例においては、応力振幅は、時間によってバラバラであり、また、振幅のピーク間の時間の間隔もバラバラである。
【００７９】
このような場合においては、解析手法上、時間軸に沿って、亀裂進行を解析していく手法となっているため、時間軸を応力の変曲点ごとに区切り、その変曲点毎の時間を区間として、その区間の応力振幅Δσを元に微分方程式を用いて解析を行い、亀裂進行判定する。１つの区間の亀裂進行判定後、次の区間の応力振幅Δσを元に微分方程式を用いて解析を行い、亀裂進行判定するという手順で処理を進めるようにする。
【００８０】
例えば、図７に示されるように、σmin1，σmax1，σmin2，σmax2，σmin3，σmax3は、それぞれ、応力の変曲点であり、時刻t0,t1,t2,t3,t4,t5における応力である。σmin1乃至σmax1の振幅の区間（時間）Ｔ1（t0乃至t1）においては、Δσ１＝σmax1−σmin1を元に微分方程式を用いて解析が行われる。次に、σmax1乃至σmin2の振幅の区間Ｔ2（t1乃至t2）においては、Δσ２＝σmax1−σmin2を元に微分方程式を用いて解析が行われる。σmin2乃至σmax2の振幅の区間Ｔ3（t2乃至t3）においては、Δσ３＝σmax2−σmin2を元に微分方程式を用いて解析が行われる。次に、σmax2乃至σmin3の振幅の区間Ｔ4（t3乃至t4）においては、Δσ４＝σmax2−σmin3を元に微分方程式を用いて解析が行われる。σmin3乃至σmax3の振幅の区間Ｔ5（t4乃至t5）においては、Δσ５＝σmax3−σmin3を元に微分方程式を用いて解析が行われる。
【００８１】
以上のように、変曲点ごとに区切り、その区間のΔσを元に微分方程式を用いて解析するようにしたので、どのような負荷波形にも対応して亀裂を予測することができる。
【００８２】
なお、微分方程式が予め生成されている場合などには、本ステップＳ１５は省略しても構わない。
【００８３】
（亀裂予測ステップ）
ステップＳ１６においては、亀裂予測部１６は、ステップＳ１５により生成された微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。なお、ステップＳ１５を行わない場合には、本ステップＳ１６においては外部機器などから取得した微分方程式を計算することにより構造体Ｄに発生する亀裂が予測される。
【００８４】
微分方程式を計算する際に、構造体Ｄに加わる応力をまず再現するために、構造体モデルＭDに対して荷重条件を付与して応力解析を行う。図８は、構造体モデルＭDに付与する荷重条件の一例を示している。図８に示す例では、構造体モデルＭDにおいて、Ｘ軸方向左側の面を固定（拘束）した状態で、Ｘ軸方向右側の面に引張荷重を加える。そして、この荷重条件下で微分方程式を計算することにより、時間の経過が伴う各要素Ｅ０における亀裂変数φの変化が得られる。
【００８５】
図９は、図８に示すように構造体モデルＭDに荷重条件を付与した場合の、ある時刻における振幅負荷エネルギーＡの分布を示している。図９には、振幅負荷エネルギーＡの等しい等エネルギー面が破線で示されている。構造体モデルＭD内において、等エネルギー面が円弧状に広がっており、亀裂の先端である要素Ｅ１のＹ軸方向下面に近いほど振幅負荷エネルギーＡが大きい。
【００８６】
図１０は、図８に示すように構造体モデルＭDに荷重条件を付与した場合の、ある時刻における亀裂変数φの分布を示している。図１０には、亀裂変数φの等しい等亀裂変数面が破線で示されている。構造体モデルＭD内において、等亀裂変数面は、亀裂の先端である要素Ｅ１のＹ軸方向下面からＹ軸方向下方に延びる楕円弧状に広がっている。内側の等亀裂変数面ほど亀裂変数φが大きい。
【００８７】
ステップＳ１６では、所定時間が経過した後に、亀裂変数φが「１」以上である要素Ｅ０に亀裂が発生するものと予測される。例えば、要素Ｅ１のＹ軸方向下側の３つの要素Ｅ０の亀裂変数φが「１」以上であった場合には、図１１に示すように、当該３つの要素Ｅ０に亀裂が発生したものとして当該３つの要素Ｅ０を要素Ｅ１に変更する。
【００８８】
また、微分方程式の計算を進める過程において、亀裂変数φが「１」以上となった要素Ｅ０の振幅負荷エネルギーＡを順次「０」に変更していくことが好ましい。これにより、構造体モデルＭDの亀裂を順次更新しながら微分方程式の計算を進めることができるため、より正確に亀裂を予測することができる。
【００８９】
このように、ステップＳ１６においては、微分方程式を計算することによって、所定時間が経過した後の構造体モデルＭDにおける亀裂を有する要素Ｅ１の分布が得られる。そして、構造体モデルＭDにおける亀裂を有する要素Ｅ１の分布によって、構造体Ｄに発生する亀裂を予測することができる。
【００９０】
［微分方程式の変形例］
ステップＳ１５で生成される微分方程式は、フェーズフィールド法の考え方に基づいて生成された微分方程式（４），（５）に限定されず、適宜変更可能である。以下、本技術で利用可能な微分方程式の変形例について説明する。
【００９１】
１．構造体Ｄを形成する材料に応じたカスタマイズフェーズフィールド法の考え方に基づいて生成された微分方程式（４），（５）では、空間座標の２階微分に比例する拡散項や弾性エネルギーｆelastの微分項が含まれているため、多岐にわたる材料で形成された構造体Ｄに適用可能である。つまり、微分方程式（４），（５）では、高い汎用性が得られる。
【００９２】
この一方で、構造体Ｄを形成する材料によっては、微分方程式（４），（５）には不要な項が含まれていることになる。したがって、ステップＳ１５では、微分方程式について、構造体Ｄを形成する材料に応じて、不要な項を排除するなどのカスタマイズをすることにより、構造体Ｄに発生する亀裂を迅速かつ的確に予測可能となる。
【００９３】
本技術に係る亀裂予測方法は、金属材料や樹脂材料などの高サイクル疲労破壊しやすい材料を用いて構成された構造体Ｄの亀裂を予測可能であればよい。このため、ステップＳ１５で生成する微分方程式は、少なくとも、亀裂変数φの時間微分に比例する項と、ステップＳ１４で設定した振幅負荷エネルギーＡに比例する項と、を含んでいればよい。
【００９４】
また、図１の亀裂予測装置１０は、少なくともモデル取得部１２及び亀裂予測部１６を備えていればよい。つまり、ステップＳ１１，Ｓ１３乃至Ｓ１５を実行しない場合には、亀裂予測装置１０がモデル生成部１１、亀裂変数設定部１３、振幅負荷エネルギー設定部１４、及び微分方程式生成部１５を備えていなくても構わない。また、亀裂予測装置１０は、必要に応じて、上記以外の構成を含んでいてもよい。
【００９５】
以下、構造体Ｄを形成する材料の一例を挙げ、当該材料に応じてカスタマイズされた微分方程式を例示する。なお、構造体Ｄを形成する材料は、以下に挙げるものに限定されず、任意の材料であってよい。また、各材料に対応する微分方程式も、以下に例示するものに限定されず、任意にカスタマイズ可能である。
【００９６】
（ａ）脆性破壊しにくい材料
構造体Ｄを形成する材料が脆性破壊しにくい場合、例えば、振幅負荷エネルギーＡの項以外の項を排除し、高サイクル疲労破壊のみを考慮した、次の微分方程式（９）を用いることができる。
【数９】
2019003301.tif 000011
微分方程式（９）には、亀裂変数φの時間微分の項、及び振幅負荷エネルギーＡの項のみが含まれる。なお、振幅負荷エネルギーＡの項には、フィッティング用定数ｗhcが含まれていてもよい。このように、振幅負荷エネルギーＡの項以外の項を排除して簡単化された微分方程式（９）を用いることにより、計算負荷を大幅に小さくすることができる。
【００９７】
（ｂ）弾性率Ｂが異方性を有する材料
構造体Ｄを形成する材料の弾性率Ｂが異方性を有する場合、例えば、弾性率Ｂの異方性を考慮した、次の微分方程式（１０）を用いることができる。
【数１０】
2019003301.tif 000012
【００９８】
微分方程式（１０）中、系のエネルギーＦsysは、次の式（１１）で表される。
【数１１】
2019003301.tif 000013
【００９９】
式（１１）中、勾配エネルギーｆgradは、次の式（１２）で表される。式（１１）中、弾性エネルギーｆelastは、次の式（１３）で表される。
【数１２】
2019003301.tif 000014
式（１２）中、κは、材料定数を示す。
【０１００】
【数１３】
2019003301.tif 000015
式（１３）中、εは垂直ひずみを示す。
【０１０１】
式（１３）では、弾性率Ｂijklをテンソル（行列）として扱うことにより、弾性率Ｂの異方性を適切に予測結果に反映させることができる。このため、微分方程式（１０）を計算することによって、弾性率Ｂが異方性を有する材料で形成された構造体Ｄに発生する亀裂を的確に予測することができる。
【０１０２】
なお、構造体Ｄを形成する材料の弾性率Ｂが等方性を有する場合には、弾性エネルギーｆelastは、式（１３）に代えて、次の式（１４）で表すことができる。
【数１４】
2019003301.tif 000016
式（１４）中、νはポアソン比を示し、γはせん断ひずみを示す。
【０１０３】
また、式（１４）中の弾性率Ｂは、例えば、図１２に示すような亀裂変数φに依存する関数とすることができる。図１２に示す関数では、亀裂変数φが大きくなるにつれて、弾性率Ｂが小さくなる。つまり、図１２に示す関数では、振幅負荷エネルギーＡの蓄積に伴って構造体Ｄを形成する材料の弾性が低下することを表現することができる。
【０１０４】
（ｃ）靭性値が異方性を有する材料
構造体Ｄを形成する材料の靭性値が異方性を有する場合、例えば、拡散項の係数を亀裂変数φの勾配の関数、すなわち界面の法線方向の関数とする、次の微分方程式（１５）を用いることができる。
【数１５】
2019003301.tif 000017
【０１０５】
微分方程式（１５）における系のエネルギーＦsysは、次の式（１６）で表される。
【数１６】
2019003301.tif 000018
【０１０６】
式（１６）における勾配エネルギーｆgradは、次の式（１７）で表される。式（１６）における弾性エネルギーｆelastは、次の式（１９）で表される。
【数１７】
2019003301.tif 000019
式（１７）におけるκは、材料定数を示し、次の式（１８）で表される。
【０１０７】
【数１８】
2019003301.tif 000020
式（１８）において、ａは異方性関数である。
【０１０８】
【数１９】
2019003301.tif 000021
式（１９）において、νはポアソン比を示し、γはせん断ひずみを示す。
【０１０９】
式（１８）を用いることによって、界面の向きに応じて拡散係数を変えることができ、亀裂の進行のしやすさを方向によって変化させることができる。これにより、靭性値の異方性を予測結果に適切に反映させることができる。このため、微分方程式（１５）を計算することによって、靭性値が異方性を有する材料で形成された構造体Ｄに発生する亀裂を的確に予測することができる。
【０１１０】
（ｄ）脆性破壊と高サイクル疲労破壊とが同時に進行する材料
構造体Ｄが脆性破壊と高サイクル疲労破壊とが同時に進行する材料によって形成されている場合、構造体Ｄには脆性破壊と高サイクル疲労破壊との組み合わせによる亀裂が発生する。脆性破壊と高サイクル疲労破壊との組み合わせによる亀裂を予測するために、例えば、微分方程式（２０）を用いることができる。
【数２０】
2019003301.tif 000022
【０１１１】
微分方程式（２０）における系のエネルギーＦsysは、次の式（２１）で表される。
【数２１】
2019003301.tif 000023
【０１１２】
式（２１）における勾配エネルギーｆgradは、次の式（２２）で表される。式（２１）における弾性エネルギーｆelastは、次の式（２３）で表される。
【数２２】
2019003301.tif 000024
式（２２）において、κは材料定数を示す。
【０１１３】
【数２３】
2019003301.tif 000025
式（２３）において、νはポアソン比を示し、εは垂直ひずみを示し、γはせん断ひずみを示す。
【０１１４】
微分方程式（２０）では、弾性エネルギーｆelastの解放率によって脆性破壊を解析することができるとともに、振幅負荷エネルギーＡの蓄積によって高サイクル疲労破壊を解析することもできる。したがって、微分方程式（２０）を計算することによって、構造体Ｄに発生する脆性破壊と高サイクル疲労破壊との組み合わせによる亀裂を予測可能となる。
【０１１５】
２．界面の安定化
構造体Ｄにおける亀裂を良好に表現するために、構造体モデルＭDにおいて亀裂を有さない要素Ｅ０と亀裂を有する要素Ｅ１との界面を安定化することが好ましい。つまり、亀裂を有さない要素Ｅ０の亀裂変数φが「０」近傍の値であり、亀裂を有する要素Ｅ１の亀裂変数が「１」近傍の値であり、いずれの要素Ｅもなるべく亀裂変数が「０」と「１」との中間の値である状態とならないことが好ましい。
【０１１６】
構造体モデルＭDにおいて亀裂を有さない要素Ｅ０と亀裂を有する要素Ｅ１との界面を安定化するために、例えば、次の微分方程式（２４）を用いることができる。
【数２４】
2019003301.tif 000026
【０１１７】
微分方程式（２４）における系のエネルギーＦsysは、次の式（２５）で表される。
【数２５】
2019003301.tif 000027
【０１１８】
式（２５）中、障壁エネルギーｆdoubは、次の式（２６）で表される。式（２５）中、弾性エネルギーｆelastは、次の式（２７）で表される。
【数２６】
2019003301.tif 000028
式（２６）において、κは材料定数を示す。
【０１１９】
【数２７】
2019003301.tif 000029
式（２７）におけるｅdoubはエネルギー障壁を示す。
【０１２０】
上述した式（２７）は、図１３に示すようなダブルウェル関数である。つまり、障壁エネルギーｆdoubが亀裂変数φ＝０，１において極小値をとる。このため、亀裂変数φは、「０」又は「１」近傍の値をとりやすく、「０」と「１」との中間の値をとりにくい。これにより、構造体モデルＭDにおいて亀裂を有さない要素Ｅ０と亀裂を有する要素Ｅ１との界面が安定化する。
【０１２１】
また、式（２７）では、構造体Ｄを形成する材料における亀裂の発生しにくさをエネルギー障壁ｅdoubによって表現することが可能である。つまり、亀裂が発生しやすい材料の場合にはエネルギー障壁ｅdoubを低くし、亀裂が発生しにくい材料の場合にはエネルギー障壁ｅdoubを高くすることができる。
【０１２２】
以上により、本技術の第１の実施形態によれば、構造体に発生する延性破壊のうち、高サイクル疲労（弾性疲労）による亀裂を迅速に予測することができる。
【０１２３】
なお、第１の実施形態は、構造体に発生する疲労破壊のうち、高負荷の塑性ひずみに入らない疲労（低サイクル疲労）を考慮しない場合、低負荷の高サイクル疲労（弾性疲労）しか考慮しない場合について説明した。ただし、材料によっては、構造体に発生する疲労破壊として、低負荷の高サイクル疲労（弾性疲労）だけでなく、高負荷の塑性ひずみ（低サイクル疲労）を考慮が必要な場合があり、この場合について、次の第２の実施形態で説明する。
【０１２４】
＜２．第２の実施形態＞
第２の実施形態に係る技術として、低負荷の高サイクル疲労（弾性疲労）および高負荷の低サイクル疲労（塑性疲労）を含む構造体に発生する疲労破壊を予測可能な方法について説明する。なお、以降は、第１の実施形態との差分に焦点をあてて説明する。
【０１２５】
本技術の第２の実施形態においては、微分方程式（４）に、塑性散逸エネルギーｆplastの項を導入した、次の微分方程式（２８）を利用する。
【数２８】
2019003301.tif 000030
【０１２６】
微分方程式（２８）は、微分方程式（４）において、塑性散逸エネルギーｆplastの項を追加した式である。また、微分方程式（２８）は、微分方程式（３）において、振幅負荷エネルギーＡの項を追加した式であるともいえる。すなわち、微分方程式（２８）は、微分方程式（２）に対して、繰り返し負荷に対応する低サイクル疲労解析を行うための塑性散逸エネルギーｆplastと高サイクル疲労解析を行うための振幅負荷エネルギーＡの項を追加したものである。低負荷の高サイクル疲労（弾性疲労）および高負荷の塑性ひずみ（低サイクル疲労）の両方含む場合は、この式（２８）を用いて、予測が行われる。
【０１２７】
微分方程式（２８）において、塑性散逸エネルギーｆplastおよび振幅負荷エネルギーＡは、弾性エネルギーｆelastとは異なり、微分項ではない。これは、弾性エネルギーｆelastは時間とともに解放されるのに対し、塑性散逸エネルギーｆplastおよび振幅負荷エネルギーＡは時間とともに蓄積されるためである。微分方程式（２８）では、塑性散逸エネルギーｆplastおよび振幅負荷エネルギーＡを微分項としないことにより、塑性散逸エネルギーｆplastおよび振幅負荷エネルギーＡの蓄積を表現することができる。
【０１２８】
このように、微分方程式（２８）には、弾性エネルギーｆelastの解放率を表現する弾性エネルギーｆelastの微分項と、塑性散逸エネルギーｆplastの蓄積を表現する塑性散逸エネルギーｆplastの項と、振幅負荷エネルギーＡの蓄積を表現する振幅負荷エネルギーＡの項と、が含まれる。したがって、微分方程式（２８）を計算することにより、脆性破壊及び疲労破壊（高低サイクル疲労含む）の双方を考慮して、亀裂を予測することが可能となる。
【０１２９】
このため、本技術では、金属材料や樹脂材料などの低サイクル疲労（塑性疲労）および高サイクル疲労（弾性疲労）を含む疲労破壊しやすい材料を用いて構成された構造体Ｄに発生する亀裂を正確に予測可能である。また、本技術に係る亀裂予測方法でも、フェーズフィールド法の考え方による亀裂予測方法と同様に、複数種類の材料の界面を跨ぐ亀裂を予測可能であるため、複数の材料で構成された構造体Ｄに発生する亀裂を迅速に予測可能である。更に、本技術に係る亀裂予測方法でも、フェーズフィールド法の考え方による亀裂予測方法と同様に、亀裂の形状に制約がないため、高い汎用性が得られる。
【０１３０】
［亀裂予測装置の詳細］
図１４は、本技術の第２の実施形態に係る亀裂予測装置（情報処理装置）１０１の構成を示すブロック図である。亀裂予測装置１０１は、モデル生成部１１と、モデル取得部１２と、亀裂変数設定部１３と、エネルギー設定部１１１と、微分方程式生成部１１２と、亀裂予測部１６と、を備える。
【０１３１】
図１４の亀裂予測装置１０１は、振幅負荷エネルギー設定部１４が、エネルギー設定部１１１に入れ替わった点と、微分方程式生成部１５が微分方程式生成部１１２に入れ替わった点とが、図１の亀裂予測装置１０と異なっている。図１４の亀裂予測装置１０１は、モデル生成部１１と、モデル取得部１２と、亀裂変数設定部１３と、亀裂予測部１６とを備える点が図１の亀裂予測装置１０と共通している。
【０１３２】
すなわち、エネルギー設定部１１１は、構造体モデルＭDの各要素Ｅに、振幅負荷エネルギーＡと塑性散逸エネルギーｆplastとをそれぞれ設定する。
【０１３３】
微分方程式生成部１１２は、亀裂変数設定部１３により設定された亀裂変数φと、エネルギー設定部１１１により設定された振幅負荷エネルギーＡと、エネルギー設定部１１１により設定された塑性散逸エネルギーｆplastとを利用して微分方程式を作成する。
【０１３４】
したがって、図１４の亀裂予測部１６は、微分方程式生成部１１２により生成された微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。
【０１３５】
［亀裂予測動作の例］
図１５のフローチャートは、図１４の亀裂予測装置１０１の亀裂予測処理を説明するフローチャートである。以下、図１５に沿って、第２の実施形態に係る亀裂予測処理について説明する。なお、図１５のステップＳ１１１乃至Ｓ１１３、およびＳ１１６は、図２のステップＳ１１乃至Ｓ１３、およびＳ１６と基本的に同様の処理を行うため、その詳細な説明は省略する。
【０１３６】
（モデル生成ステップ）
ステップＳ１１１において、モデル生成部１１は、構造体Ｄの構成を再現したモデル（構造体モデル）ＭDを生成する。
【０１３７】
（モデル取得ステップ）
ステップＳ１１２において、モデル取得部１２は、ステップＳ１１１で生成された構造体モデルＭDを取得する。
【０１３８】
（亀裂変数設定ステップ）
ステップＳ１１３において、亀裂変数設定部１３は、ステップＳ１１２で取得された構造体モデルＭDの各要素Ｅに亀裂の有無を表現する亀裂変数φを設定する。
【０１３９】
（振幅負荷エネルギー設定ステップ）
ステップＳ１１４において、エネルギー設定部１１１は、ステップＳ１１２で取得された構造体モデルＭDの各要素Ｅに、塑性散逸エネルギーｆplastと振幅負荷エネルギーＡとを設定する。なお、ここでは、図２のステップＳ１４との差分である、塑性散逸エネルギーｆplastの設定についてのみ説明する。
【０１４０】
すなわち、既に亀裂を有する要素Ｅ１では、塑性変形が生じないため、塑性散逸エネルギーｆplastが蓄積されない。このため、要素Ｅ１の塑性散逸エネルギーｆplastは「０」と設定される。
【０１４１】
亀裂を有さない要素Ｅ０における塑性散逸エネルギーｆplastは、当該要素Ｅ０を構成する材料に応じて実験的に得られる相当応力σと相当塑性ひずみεpとの関係に基づいて設定される。相当塑性ひずみεpは亀裂変数φに依存するため、塑性散逸エネルギーｆplastは、亀裂変数φの関数として表現される。
【０１４２】
図１６は、ステップＳ１１４において要素Ｅ０に設定される塑性散逸エネルギーｆplastの表現方法の一例を示す図である。図１６は、構造体Ｄを形成する材料から得られる相当応力−相当塑性ひずみ線図の一例を示している。図１６では、縦軸が相当応力σを示し、横軸が相当塑性ひずみεpを示している。また、図１６には降伏応力σYが示されている。
【０１４３】
図１６に示す相当応力−相当塑性ひずみ線図が得られる材料は、相当応力σが降伏応力σY未満の領域で弾性変形し、相当応力σが降伏応力σY以上の領域で塑性変形する。塑性散逸エネルギーｆplastは、相当応力σが降伏応力σY以上のとき、当該材料の塑性変形により主に熱エネルギーとして散逸されるエネルギーを表現するものである。
【０１４４】
塑性散逸エネルギーｆplastは、例えば、図１６Ａ及び図１６Ｂに斜線で示す領域の面積として定義することが可能である。図１６Ａに斜線で示す領域の面積は、相当応力σを相当塑性ひずみεpの微小増分で積分した量を利用して、例えば、次の式（２９）によって算出可能である。
【数２９】
2019003301.tif 000031
【０１４５】
また、図１６Ｂに斜線で示す領域の面積は、相当応力σと降伏応力σYとの差と、相当塑性ひずみεpと、の積を利用して、例えば、次の式（３０）によって算出可能である。
【数３０】
2019003301.tif 000032
【０１４６】
ただし、式（３０）において、相当応力σが降伏応力σYより小さい場合の塑性散逸エネルギーｆplastをゼロとする。
【０１４７】
塑性散逸エネルギーｆplastの表現方法は、構造体Ｄを形成する材料や物理現象などに応じて、正確に亀裂を予測できるように使い分けることが可能である。なお、塑性散逸エネルギーｆplastを表現する関数は、式（２９）及び式（３０）に限定されず、相当応力σと相当塑性ひずみεpとの関係に基づいて適宜作成可能である。
【０１４８】
なお、構造体モデルＭDに振幅負荷エネルギーＡおよび塑性散逸エネルギーｆplastが予め設定されている場合などには、本ステップＳ１１４は省略しても構わない。
【０１４９】
（微分方程式生成ステップ）
ステップＳ１１５において、微分方程式生成部１１２は、ステップＳ１１３で設定した亀裂変数φと、ステップＳ１１４で設定された振幅負荷エネルギーＡと、ステップＳ１１４で設定された塑性散逸エネルギーｆplastが利用されて微分方程式を作成する。
【０１５０】
ステップＳ１１５において生成される微分方程式としては、一例として、上記の微分方程式（２８）が挙げられる。また、ステップＳ１１５においては、微分方程式（２８）を改良した、次の微分方程式（３１）が生成されてもよい。なお、式（３１）は、拡散方程式でもある。
【数３１】
2019003301.tif 000033
【０１５１】
微分方程式（３１）には、障壁エネルギーｆdoubの微分項のフィッティング用定数ｗdoub、弾性エネルギーｆelastの微分項のフィッティング用定数ｗelast、塑性散逸エネルギーｆplastの項のフィッティング用定数ｗplast、及び振幅負荷エネルギーＡの項のフィッティング用定数ｗhcが導入されている。これにより、構造体Ｄの構成などに応じて、障壁エネルギーｆdoubの微分項と、弾性エネルギーｆelastの微分項と、塑性散逸エネルギーｆplastの項と、振幅負荷エネルギーＡの項と、のそれぞれの重み付けを最適化することが可能である。これにより、構造体Ｄに発生する亀裂を更に正確に予測可能となる。
【０１５２】
ここで、第１の実施形態において、振幅負荷エネルギーＡについて上述された微分方程式の変形例、例えば、式（９）、式（１０）、式（１５）、式（２０）、および式（２４）などは、振幅負荷エネルギーＡを塑性散逸エネルギーｆplastに変えて、同様に適用される。
【０１５３】
なお、微分方程式が予め生成されている場合などには、本ステップＳ１１５は省略しても構わない。
【０１５４】
（亀裂予測ステップ）
ステップＳ１１６において、亀裂予測部１６は、ステップＳ１１５で生成された微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。なお、ステップＳ１１５を行わない場合には、本ステップＳ１１６においては外部機器などから取得した微分方程式を計算することにより構造体Ｄに発生する亀裂が予測される。
【０１５５】
以上により、本技術の第２の実施形態によれば、構造体に発生する疲労破壊（低サイクル疲労（塑性疲労）および高サイクル疲労（弾性疲労））による亀裂を迅速に予測することができる。
【０１５６】
すなわち、本技術の第２の実施形態によれば、脆性亀裂、低サイクル疲労（塑性疲労）による亀裂、高サイクル繰り返し疲労による亀裂など、すべての亀裂現象が一気通貫で解析可能となる。すなわち、本技術によれば、構造体Ｄに発生する亀裂を迅速かつより正確に予測可能となる。
【０１５７】
[その他の実施形態]
以上、本技術の実施形態について説明したが、本技術は上述の実施形態にのみ限定されるものではなく、本技術の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。
【０１５８】
例えば、上記実施形態では構造体モデルＭDの要素Ｅを１次要素としたが、必要に応じて、構造体モデルＭDの要素Ｅを２次要素としてもよい。この場合、構造体モデルＭDの各要素Ｅ内における亀裂変数φの分布を考慮することができるため、構造体Ｄの亀裂を更に正確に予測可能となる。
【０１５９】
以上のように、本技術によれば、亀裂有無を表現する場の亀裂変数φを振幅負荷エネルギーに比例する項をソース項とした拡散方程式で時間発展させた方程式を数値的に解くようにした。
【０１６０】
これにより、界面をまたぐような亀裂にも対応し、亀裂形状も自由に変化できるフェーズフィールド法の利点を生かしながら、金属／樹脂の弾性域での高サイクル疲労破壊を予測可能になる。
【０１６１】
さらに、本技術によれば、亀裂有無を表現する場の亀裂変数φを塑性散逸エネルギーに比例する項と、亀裂有無を表現する場の亀裂変数φをソース項とした拡散方程式で時間発展させた方程式を数値的に解くようにした。
【０１６２】
これにより、界面をまたぐような亀裂にも対応し、亀裂形状も自由に変化できるフェーズフィールド法の利点を生かしながら、金属／樹脂の塑性域での低サイクル繰り返し疲労破壊および弾性域での高サイクル疲労破壊を予測可能になる。これにより、亀裂をより正確に予測することができる。
【０１６３】
＜３．コンピュータ＞
上述した一連の処理は、ハードウエアにより実行させることもできるし、ソフトウエアにより実行させることもできる。一連の処理をソフトウエアにより実行する場合には、そのソフトウエアを構成するプログラムが、コンピュータにインストールされる。ここでコンピュータには、専用のハードウエアに組み込まれているコンピュータや、各種のプログラムをインストールすることで、各種の機能を実行することが可能な、例えば汎用のパーソナルコンピュータ等が含まれる。
【０１６４】
図１７は、上述した一連の処理をプログラムにより実行するコンピュータのハードウエアの構成例を示すブロック図である。
【０１６５】
図１７に示されるコンピュータにおいて、CPU（Central Processing Unit）３０１、ROM（Read Only Memory）３０２、RAM（Random Access Memory）３０３は、バス３０４を介して相互に接続されている。
【０１６６】
バス３０４にはまた、入出力インタフェース３０５も接続されている。入出力インタフェース３０５には、入力部３０６、出力部３０７、記憶部３０８、通信部３０９、およびドライブ３１０が接続されている。
【０１６７】
入力部３０６は、例えば、キーボード、マウス、マイクロホン、タッチパネル、入力端子などよりなる。出力部３０７は、例えば、ディスプレイ、スピーカ、出力端子などよりなる。記憶部３０８は、例えば、ハードディスク、RAMディスク、不揮発性のメモリなどよりなる。通信部３０９は、例えば、ネットワークインタフェースよりなる。ドライブ３１０は、磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリなどのリムーバブルメディア３１１を駆動する。
【０１６８】
以上のように構成されるコンピュータでは、CPU３０１およびバス３０４を介して、RAM３０３にロードして実行することにより、上述した一連の処理が行われる。RAM３０３にはまた、CPU３０１が各種の処理を実行する上において必要なデータなども適宜記憶される。
【０１６９】
コンピュータ（CPU３０１）が実行するプログラムは、例えば、パッケージメディア等としてのリムーバブルメディア３１１に記録して適用することができる。その場合、プログラムは、リムーバブルメディア３１１をドライブ３１０に装着することにより、入出力インタフェース３０５を介して、記憶部３０８にインストールすることができる。
【０１７０】
また、このプログラムは、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線または無線の伝送媒体を介して提供することもできる。その場合、プログラムは、通信部３０９で受信し、記憶部３０８にインストールすることができる。
【０１７１】
その他、このプログラムは、ROM３０２や記憶部３０８に、あらかじめインストールしておくこともできる。
【０１７２】
また、本技術の実施の形態は、上述した実施の形態に限定されるものではなく、本技術の要旨を逸脱しない範囲において種々の変更が可能である。
【０１７３】
なお、以上において、１つの装置（または処理部）として説明した構成を分割し、複数の装置（または処理部）として構成するようにしてもよい。逆に、以上において複数の装置（または処理部）として説明した構成をまとめて１つの装置（または処理部）として構成されるようにしてもよい。また、各装置（または各処理部）の構成に上述した以外の構成を付加するようにしてももちろんよい。さらに、システム全体としての構成や動作が実質的に同じであれば、ある装置（または処理部）の構成の一部を他の装置（または他の処理部）の構成に含めるようにしてもよい。つまり、本技術は、上述した実施の形態に限定されるものではなく、本技術の要旨を逸脱しない範囲において種々の変更が可能である。
【０１７４】
以上、添付図面を参照しながら本開示の好適な実施形態について詳細に説明したが、本開示はかかる例に限定されない。本開示の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本開示の技術的範囲に属するものと了解される。
【０１７５】
なお、本技術は以下のような構成も取ることができる。
（１）  所定の構造体に対応する構造体モデルを取得するモデル取得部と、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
を備える情報処理装置。
（２）  前記振幅負荷エネルギーは、繰り返し負荷時の応力の最大値と最小値の差分と、歪みの最大値と最小値の差分とが積算されて設定される
前記（１）に記載の情報処理装置。
（３）  前記振幅負荷エネルギーに係る定数は、対応する材料のS-N線図の高サイクル疲労領域の特性線の傾きに応じた値として設定される
前記（１）または（２）に記載の情報処理装置。
（４）  前記振幅負荷エネルギーに係る定数は、対応する材料界面のS-N線図の高サイクル疲労剥離領域の特性線の傾きに応じた値として設定される
前記（１）または（２）に記載の情報処理装置。
（５）  前記亀裂予測部は、非周期的な繰り返し負荷に対して、応力時間変化の変曲点毎に解析区間を区切り、区切られた解析区間毎に前記構造体に発生する亀裂を予測する
前記（１）乃至（４）のいずれかに記載の情報処理装置。
（６）  前記亀裂予測部は、前記構造体モデルの各位置に設定され、塑性変形時に散逸されるエネルギーを前記亀裂変数を利用して表現する塑性散逸エネルギーに比例する項をさらに含む前記微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
前記（１）乃至（５）のいずれかに記載の情報処理装置。
（７）  前記塑性散逸エネルギーは、相当応力を相当塑性ひずみの微小増分で積分した量を利用して設定される
前記（６）に記載の情報処理装置。
（８）  前記塑性散逸エネルギーは、相当応力と降伏応力との差と、相当塑性ひずみとの積を利用して設定され、相当応力が降伏応力より小さい場合にゼロとされる
前記（６）または（７）に記載の情報処理装置。
（９）  前記微分方程式は、空間座標の２階微分に比例する拡散項をさらに含む
前記（６）乃至（８）のいずれかに記載の情報処理装置。
（１０）  所定の構造体に対応する構造体モデルを取得し、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
情報処理方法。
（１１）  所定の構造体に対応する構造体モデルを取得するモデル取得部と、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
して、コンピュータを機能させるプログラム。
【符号の説明】
【０１７６】
１０  亀裂予測装置（情報処理装置），  １１  モデル生成部，  １２  モデル取得部，１３  亀裂変数設定部，  １４  振幅負荷エネルギー設定部，１５  微分方程式生成部，  １６  亀裂予測部，１１０  亀裂予測装置，  １１１  エネルギー設定部，  １１２  微分方程式生成部  ，ＭD  構造体モデル，  Ｅ，Ｅ０，Ｅ１  要素
【図１】
2019003301.tif 000034
【図２】
2019003301.tif 000035
【図３】
2019003301.tif 000036
【図４】
2019003301.tif 000037
【図５】
2019003301.tif 000038
【図６】
2019003301.tif 000039
【図７】
2019003301.tif 000040
【図８】
2019003301.tif 000041
【図９】
2019003301.tif 000042
【図１０】
2019003301.tif 000043
【図１１】
2019003301.tif 000044
【図１２】
2019003301.tif 000045
【図１３】
2019003301.tif 000046
【図１４】
2019003301.tif 000047
【図１５】
2019003301.tif 000048
【図１６】
2019003301.tif 000049
【図１７】
2019003301.tif 000050