特開 2019-3301 情報処理装置、情報処理方法、及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】特開2019-3301(P2019-3301A)

(43)【公開日】2019年1月10日

(54)【発明の名称】情報処理装置、情報処理方法、及びプログラム

(51)【国際特許分類】

   G06F 17/50 20060101AFI20181207BHJP

【ＦＩ】

   G06F17/50 612H

【審査請求】未請求

【請求項の数】11

【出願形態】ＯＬ

【全頁数】28

(21)【出願番号】特願2017-115657(P2017-115657)

(22)【出願日】2017年6月13日

(71)【出願人】

【識別番号】316005926

【氏名又は名称】ソニーセミコンダクタソリューションズ株式会社

【住所又は居所】神奈川県厚木市旭町四丁目１４番１号

(74)【代理人】

【識別番号】100121131

【弁理士】

【氏名又は名称】西川 孝

(74)【代理人】

【識別番号】100082131

【弁理士】

【氏名又は名称】稲本 義雄

(72)【発明者】

【氏名】本郷 一泰

【住所又は居所】神奈川県厚木市旭町四丁目１４番１号 ソニーセミコンダクタソリューションズ株式会社内

(72)【発明者】

【氏名】平野 嵩明

【住所又は居所】神奈川県厚木市旭町四丁目１４番１号 ソニーセミコンダクタソリューションズ株式会社内

【テーマコード（参考）】

5B046

【Ｆターム（参考）】

5B046AA08

5B046JA08

(57)【要約】

【課題】亀裂をより正確に予測することができる。
【解決手段】モデル取得部は、モデル生成部または図示せぬ外部機器などから構造体モデルＭ_Dを取得する。亀裂を有さない要素Ｅ０における振幅負荷エネルギーＡは、当該要素Ｅ０を構成する材料に応じて実験的に得られる相当応力σと相当弾性ひずみεとの関係に基づいて設定される。相当弾性ひずみεは亀裂変数φに依存するため、振幅負荷エネルギーＡは亀裂変数φの関数として表現される。亀裂予測部は、この振幅エネルギーに比例する項を有する微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。本開示は、例えば、亀裂の予測を行う亀裂予測装置に適用することができる。
【選択図】図１

【特許請求の範囲】

【請求項1】

所定の構造体に対応する構造体モデルを取得するモデル取得部と、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
を備える情報処理装置。

【請求項2】

前記振幅負荷エネルギーは、繰り返し負荷時の応力の最大値と最小値の差分と、歪みの最大値と最小値の差分とが積算されて設定される
請求項１に記載の情報処理装置。

【請求項3】

前記振幅負荷エネルギーに係る定数は、対応する材料のS-N線図の高サイクル疲労領域の特性線の傾きに応じた値として設定される
請求項１に記載の情報処理装置。

【請求項4】

前記振幅負荷エネルギーに係る定数は、対応する材料界面のS-N線図の高サイクル疲労剥離領域の特性線の傾きに応じた値として設定される
請求項１に記載の情報処理装置。

【請求項5】

前記亀裂予測部は、非周期的な繰り返し負荷に対して、応力時間変化の変曲点毎に解析区間を区切り、区切られた解析区間毎に前記構造体に発生する亀裂を予測する
請求項１に記載の情報処理装置。

【請求項6】

前記亀裂予測部は、前記構造体モデルの各位置に設定され、塑性変形時に散逸されるエネルギーを前記亀裂変数を利用して表現する塑性散逸エネルギーに比例する項をさらに含む前記微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
請求項１に記載の情報処理装置。

【請求項7】

前記塑性散逸エネルギーは、相当応力を相当塑性ひずみの微小増分で積分した量を利用して設定される
請求項６に記載の情報処理装置。

【請求項8】

前記塑性散逸エネルギーは、相当応力と降伏応力との差と、相当塑性ひずみとの積を利用して設定され、相当応力が降伏応力より小さい場合にゼロとされる
請求項６に記載の情報処理装置。

【請求項9】

前記微分方程式は、空間座標の２階微分に比例する拡散項をさらに含む
請求項１に記載の情報処理装置。

【請求項10】

所定の構造体に対応する構造体モデルを取得し、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
情報処理方法。

【請求項11】

所定の構造体に対応する構造体モデルを取得するモデル取得部と、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
して、コンピュータを機能させるプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、情報処理装置、情報処理方法、及びプログラムに関し、特に、亀裂をより正確に予測することができる情報処理装置、情報処理方法、及びプログラムに関する。

【背景技術】

【0002】

一般的に、半導体デバイス等の各種構造体には、その製造過程などにおいて機械的応力等の様々な応力が加わる。構造体にこのような応力が加わると、その内部に亀裂が発生することがある。このような亀裂の発生を防止するために、構造体に予め発生する可能性のある亀裂を予測する技術が用いられる。

【0003】

特許文献１には、亀裂を含む断面上の異種材料の面積比に応じて各材料の応力とひずみの関係を単一の材料の応力のひずみの関係に置き換えたときの材料定数を用いて求めた破壊力学パラメータＪ積分値が亀裂先端の材料の弾塑性破壊靭性値に等しくなるときの荷重を破壊荷重とする技術が開示されている。

【0004】

特許文献２には、複数種類の材料の界面を跨ぐ亀裂を予測可能な技術が開示されている。特許文献２に係る技術では、構造体の内部で仮想的に亀裂を進行させたときのエネルギー解放率を算出し、エネルギー解放率の大きい方向に亀裂が進行するものと予測する。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特開２００９−１６００２８号公報

【特許文献2】特開２０１１−２０４０８１号公報

【発明の概要】

【発明が解決しようとする課題】

【0006】

近年、半導体デバイス等の構造体の多様化に伴い、構造体を構成する材料として金属材料や樹脂材料が広く用いられるようになってきている。したがって、金属材料や樹脂材料を用いて構成された構造体に発生する可能性のある亀裂を正確に予測可能な技術が要求されている。

【0007】

特許文献１に係る技術では、界面をまたぐ亀裂に対応することが難しく、疲労破壊への対応が難しかった。特許文献２に係る技術では、計算負荷が高く、予め定めた特定面以外の面を解析することが困難であった。

【0008】

本開示は、このような状況に鑑みてなされたものであり、亀裂をより正確に予測することができるものである。

【課題を解決するための手段】

【0009】

本技術の一側面の情報処理装置は、所定の構造体に対応する構造体モデルを取得するモデル取得部と、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部とを備える。

【0010】

前記振幅負荷エネルギーは、繰り返し負荷時の応力の最大値と最小値の差分と、歪みの最大値と最小値の差分とが積算されて設定される。

【0011】

前記振幅負荷エネルギーに係る定数は、対応する材料のS-N線図の高サイクル疲労領域の特性線の傾きに応じた値として設定される。

【0012】

前記振幅負荷エネルギーに係る定数は、対応する材料界面のS-N線図の高サイクル疲労剥離領域の特性線の傾きに応じた値として設定される。

【0013】

前記亀裂予測部は、非周期的な繰り返し負荷に対して、応力時間変化の変曲点毎に解析区間を区切り、区切られた解析区間毎に前記構造体に発生する亀裂を予測することができる。

【0014】

前記亀裂予測部は、前記構造体モデルの各位置に設定され、塑性変形時に散逸されるエネルギーを前記亀裂変数を利用して表現する塑性散逸エネルギーに比例する項をさらに含む前記微分方程式を計算することにより、前記構造体に発生する亀裂を予測することができる。

【0015】

前記塑性散逸エネルギーは、相当応力を相当塑性ひずみの微小増分で積分した量を利用して設定される。

【0016】

前記塑性散逸エネルギーは、相当応力と降伏応力との差と、相当塑性ひずみとの積を利用して設定され、相当応力が降伏応力より小さい場合にゼロとされる。

【0017】

前記微分方程式は、空間座標の２階微分に比例する拡散項をさらに含むことができる。

【0018】

本技術の一側面の情報処理方法は、所定の構造体に対応する構造体モデルを取得し、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する。

【0019】

本技術の一側面のプログラムは、所定の構造体に対応する構造体モデルを取得するモデル取得部と、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部として、コンピュータを機能させる。

【0020】

本技術の一側面においては、所定の構造体に対応する構造体モデルが取得される。そして、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂が予測される。

【発明の効果】

【0021】

本技術によれば、高サイクル疲労破壊を予測することができる。この結果、亀裂をより正確に予測することができる。

【0022】

なお、本明細書に記載された効果は、あくまで例示であり、本技術の効果は、本明細書に記載された効果に限定されるものではなく、付加的な効果があってもよい。

【図面の簡単な説明】

【0023】

【図1】本技術の一実施形態に係る亀裂予測装置（情報処理装置）の構成例を示すブロック図である。

【図2】図１の亀裂予測装置の亀裂予測処理について説明するフローチャートである。

【図3】亀裂予測処理において生成される構造体モデルの一例を示す図である。

【図4】繰り返し負荷に対応した対象材料にかかる応力振幅を定義するための図である。

【図5】対象材料のS-N線図を表す図である。

【図6】材料の低サイクル、高サイクルの疲労特性を記述する際の式を説明する図である。

【図7】対象物に対して想定される負荷が、周期的な振動負荷でない場合の例を説明する図である。

【図8】構造体モデルに加える荷重条件の一例を説明するための図である。

【図9】構造体モデルにおける振幅負荷エネルギーの分布の一例を示す図である。

【図10】構造体モデルにおける亀裂変数の分布の一例を示す図である。

【図11】構造体モデルにおいて予測される亀裂の一例を示す図である。

【図12】弾性率と亀裂変数との関係の一例を示す図である。

【図13】障壁エネルギーと亀裂変数との関係の一例を示す図である。

【図14】本技術の他の実施形態に係る亀裂予測装置（情報処理装置）の構成例を示すブロック図である。

【図15】図１４の亀裂予測装置の亀裂予測処理について説明するフローチャートである。

【図16】亀裂予測処理において設定される塑性散逸エネルギーの表現方法の一例を示す図である。

【図17】コンピュータの構成例を示すブロック図である。

【発明を実施するための形態】

【0024】

以下、本開示を実施するための形態（以下実施の形態とする）について説明する。なお、説明は以下の順序で行う。また、図面には、適宜相互に直交するＸ軸、Ｙ軸、及びＺ軸が示されている。Ｘ軸、Ｙ軸、及びＺ軸は全図において共通である。
０．概要
１．第１の実施形態
２．第２の実施形態
３.コンピュータ

【0025】

＜０．概要＞
［亀裂予測方法の概要］
本技術に係る亀裂予測方法（情報処理方法）の概要について説明する。本技術に係る亀裂予測方法では、フェーズフィールド（Phase Field）法の考え方を応用して構造体Ｄに発生する亀裂を予測する。まず、本技術に関連するフェーズフィールド法の考え方による亀裂予測方法について説明する。

【0026】

（フェーズフィールド法の考え方による亀裂予測方法）
構造体Ｄ内のエネルギーＦは、障壁エネルギーｆ_doub、勾配エネルギーｆ_grad、及び弾性エネルギーｆ_elastを利用して、式（１）で表される。

【数1】

【0027】

フェーズフィールド法の考え方によると、式（１）から、次の微分方程式（２）を導出することができる。

【数2】

【0028】

微分方程式（２）の左辺は、易動度Ｍの逆数と、亀裂の発生の有無を表現する亀裂変数φの時間微分と、の積で構成されている。微分方程式（２）の右辺は、空間座標の２階微分の拡散項∇（ξ∇φ）と、障壁エネルギーｆ_doubの微分項と、弾性エネルギーｆ_elastの微分項と、で構成されている。微分方程式（２）では、弾性エネルギーｆ_elastの微分項によって、弾性エネルギーｆ_elastの解放率が表現されている。

【0029】

微分方程式（２）の計算を行うために、まず構造体Ｄの各位置にそれぞれ亀裂変数φが設定される。より詳細には、亀裂を有さない位置と、亀裂を有する位置と、で異なる亀裂変数φが設定される。例えば、亀裂を有さない位置の亀裂変数φが「０」と設定され、亀裂を有する位置の亀裂変数φが「１」と設定される。

【0030】

そして、微分方程式（２）の計算を進めると、時間の経過に伴い、亀裂変数φを「０」と設定した位置のうち、亀裂変数φが「１」以上となる位置が現れる。フェーズフィールド法の考え方による亀裂予測方法では、所定時間が経過した後に、亀裂変数φが「１」以上である位置に亀裂が発生しているものと予測することができる。

【0031】

フェーズフィールド法の考え方による亀裂予測方法では、微分方程式（２）を計算することにより、亀裂を迅速に予測可能である。また、フェーズフィールド法の考え方による亀裂予測方法では、複数種類の材料の界面を跨ぐ亀裂を予測可能であるため、複数の材料で構成された構造体Ｄに発生する亀裂を予測可能である。更に、フェーズフィールド法の考え方による亀裂予測方法では、亀裂の形状に制約がないため、高い汎用性が得られる。

【0032】

フェーズフィールド法の考え方による亀裂予測方法では、微分方程式（２）に含まれる弾性エネルギーｆ_elastの微分項で表現される弾性エネルギーｆ_elastの解放率によって、脆性破壊による亀裂を予測することができる。しかし、フェーズフィールド法の考え方による亀裂予測方法では、微分方程式（２）に塑性変形に対応する項が含まれていないため、塑性変形を伴う疲労破壊による亀裂を予測することができない。

【0033】

したがって、フェーズフィールド法の考え方を利用した亀裂予測方法では、金属材料や樹脂材料などの疲労破壊しやすい材料を用いて構成された構造体Ｄに発生する亀裂を正確に予測することは困難である。

【0034】

そこで、本技術の発明者は、フェーズフィールド法の考え方を応用し、塑性変形時に主に熱として散逸されるエネルギー（以下、「塑性散逸エネルギーｆ_plast」と称する）を含む項を微分方程式（２）に導入することにより、疲労破壊による亀裂を予測可能となることを見出した。以下、本技術に係るフェーズフィールド法の考え方を応用した亀裂予測方法について説明する。

【0035】

（フェーズフィールド法の考え方を応用した亀裂予測方法）
本技術に係るフェーズフィールド法の考え方を応用した亀裂予測方法では、微分方程式（２）に塑性散逸エネルギーｆ_plastの項を導入した、次の微分方程式（３）を利用する。

【数3】

【0036】

微分方程式（３）において、塑性散逸エネルギーｆ_plastは、弾性エネルギーｆ_elastとは異なり、微分項ではない。これは、弾性エネルギーｆ_elastは時間とともに解放されるのに対し、塑性散逸エネルギーｆ_plastは時間とともに蓄積されるためである。微分方程式（３）では、塑性散逸エネルギーｆ_plastを微分項としないことにより、塑性散逸エネルギーｆ_plastの蓄積を表現することができる。

【0037】

このように、微分方程式（３）には、弾性エネルギーｆ_elastの解放率を表現する弾性エネルギーｆ_elastの微分項と、塑性散逸エネルギーｆ_plastの蓄積を表現する塑性散逸エネルギーｆ_plastの項と、が含まれる。したがって、微分方程式（３）を計算することにより、脆性破壊及び疲労破壊の双方を考慮して、亀裂を予測することが可能となる。

【0038】

なお、上述した塑性散逸エネルギーｆ_plastでは、疲労破壊のうち、低サイクル疲労（塑性疲労）を解析することはできる。しかしながら、高サイクル疲労（弾性疲労）解析について対応することが困難であった。

【0039】

そこで、本技術の発明者は、さらに、フェーズフィールド法の考え方を応用し、微分方程式（２）に、高サイクルの繰り返し疲労（負荷）を表すエネルギー（以下、振幅負荷エネルギーＡと称する）を含む項を導入することにより、高サイクル疲労（弾性疲労）による亀裂を予測可能となることを見出した。

【0040】

（フェーズフィールド法の考え方を応用した亀裂予測方法）
本技術に係るフェーズフィールド法の考え方を応用した亀裂予測方法では、微分方程式（２）に、振幅負荷エネルギーＡの項を導入した、次の微分方程式（４）を利用する。

【数4】

【0041】

微分方程式（４）は、微分方程式（３）において、塑性散逸エネルギーｆ_plastを、振幅負荷エネルギーＡに入れ替えた式である。すなわち、微分方程式（４）は、微分方程式（２）に対して、繰り返し負荷に対応する高サイクル疲労解析を行うための振幅負荷エネルギーＡの項を追加したものである。

【0042】

微分方程式（４）において、振幅負荷エネルギーＡは、弾性エネルギーｆ_elastとは異なり、微分項ではない。これは、弾性エネルギーｆ_elastは時間とともに解放されるのに対し、振幅負荷エネルギーＡは時間とともに蓄積されるためである。微分方程式（４）では、振幅負荷エネルギーＡを微分項としないことにより、振幅負荷エネルギーＡの蓄積を表現することができる。

【0043】

このように、微分方程式（４）には、弾性エネルギーｆ_elastの解放率を表現する弾性エネルギーｆ_elastの微分項と、振幅負荷エネルギーＡの蓄積を表現する振幅負荷エネルギーＡの項と、が含まれる。したがって、微分方程式（４）を計算することにより、脆性破壊及び高サイクル疲労破壊の双方を考慮して、亀裂を予測することが可能となる。

【0044】

このため、本技術では、金属材料や樹脂材料などの疲労破壊のうち、高サイクル疲労（弾性疲労）しやすい材料を用いて構成された構造体Ｄに発生する亀裂を正確に予測可能である。また、本技術に係る亀裂予測方法でも、フェーズフィールド法の考え方による亀裂予測方法と同様に、複数種類の材料の界面を跨ぐ亀裂を予測可能であるため、複数の材料で構成された構造体Ｄに発生する亀裂を迅速に予測可能である。更に、本技術に係る亀裂予測方法でも、フェーズフィールド法の考え方による亀裂予測方法と同様に、亀裂の形状に制約がないため、高い汎用性が得られる。

【0045】

＜１．第１の実施形態＞
［亀裂予測装置の詳細］
図１は、本技術の第１の実施形態に係る亀裂予測装置（情報処理装置）１０の構成を示すブロック図である。なお、第１の実施形態においては、疲労破壊のうち弾性疲労（高サイクル疲労）下について述べていく。

【0046】

亀裂予測装置１０は、モデル生成部１１と、モデル取得部１２と、亀裂変数設定部１３と、振幅負荷エネルギー設定部１４と、微分方程式生成部１５と、亀裂予測部１６と、を備える。

【0047】

モデル生成部１１は、構造体Ｄの構成を再現したモデル（構造体モデル）Ｍ_Dを生成する。モデル取得部１２は、モデル生成部１１により生成された構造体モデルＭ_Dを取得する。

【0048】

亀裂変数設定部１３は、モデル取得部１２により取得された構造体モデルＭ_Dの各要素Ｅに亀裂の有無を表現する亀裂変数φを設定する。振幅負荷エネルギー設定部１４は、モデル取得部１２により取得された構造体モデルＭ_Dの各要素Ｅに、振幅負荷エネルギーＡを設定する。

【0049】

微分方程式生成部１５は、亀裂変数設定部１３により設定された亀裂変数φと、振幅負荷エネルギー設定部１４により設定された振幅負荷エネルギーＡとを利用して微分方程式を作成する。

【0050】

亀裂予測部１６は、微分方程式生成部１５により生成された微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。

【0051】

［亀裂予測動作の例］
図２のフローチャートは、亀裂予測装置１０の亀裂予測処理を説明するフローチャートである。図３乃至図１３は、図２に示す各ステップを説明するための図である。以下、図２に沿って、図３乃至図１３を適宜参照しながら、第１の実施形態に係る亀裂予測処理について説明する。

【0052】

（モデル生成ステップ）
ステップＳ１１において、モデル生成部１１は、構造体Ｄの構成を再現したモデル（構造体モデル）Ｍ_Dを生成する。構造体モデルＭ_Dによって任意の構造体Ｄの構成を再現可能である。構造体モデルＭ_Dによってその構成を再現可能な構造体Ｄとしては、例えば、半導体デバイスなどの各種デバイスが挙げられる。

【0053】

亀裂予測方法には、有限要素法（FEM：Finite Element Method）や差分法（FDM：Finite Difference Method）を利用することができる。また、陰解法や陽解法を利用することができる。有限要素法では、任意の形状に対応可能であり、高い汎用性が得られる。差分法では、計算の並列化が容易であり、計算が早いというメリットが得られる。陰解法では、タイムステップを大きくとれるというメリットが得られる。本実施形態では、有限要素法を用いるため、構造体モデルＭ_Dが複数の要素Ｅから構成される。

【0054】

図３は、ステップＳ１１において生成される構造体モデルＭ_Dを例示する図である。図３Ａは構造体モデルＭ_Dの斜視図であり、図３Ｂは構造体モデルＭ_Dの図３ＡのＡ−Ａ'線に沿った断面図である。図３に示す構造体モデルＭ_Dによってその構成が再現されている構造体Ｄでは、概形が立方体であり、上面の中央にＺ軸方向に延びる初期亀裂が形成されている。

【0055】

この場合、構造体モデルＭ_Dでは、Ｙ軸方向上面のＸ軸方向中央部にＺ軸方向に並ぶ５つの要素Ｅが亀裂を有する要素Ｅ１とされ、その他の要素Ｅが亀裂を有さない要素Ｅ０とされる。図３には、亀裂を有する要素Ｅ１が斜線で示され、亀裂を有さない要素Ｅ０が白抜きで示されている。なお、空孔などの自由空間を有する要素Ｅも、亀裂を有する要素Ｅ１と同様に扱うことが好ましい。

【0056】

以下の説明では、一例として図３に示す構造体モデルＭ_Dを用いて説明するが、他の構造体モデルＭ_Dについても同様に扱うことが可能であることは勿論である。

【0057】

なお、構造体モデルＭ_Dが予め用意されている場合などには、本ステップＳ１１は省略しても構わない。

【0058】

（モデル取得ステップ）
ステップＳ１２において、モデル取得部１２は、ステップＳ１１で生成された構造体モデルＭ_Dを取得する。なお、ステップＳ１１を行わない場合には、本ステップＳ１２においては、図示せぬ外部機器などから構造体モデルＭ_Dを取得することができる。

【0059】

（亀裂変数設定ステップ）
ステップＳ１３において、亀裂変数設定部１３は、ステップＳ１２で取得された構造体モデルＭ_Dの各要素Ｅに亀裂の有無を表現する亀裂変数φを設定する。

【0060】

具体的には、構造体モデルＭ_Dにおいて、亀裂を有さない要素Ｅ０と、亀裂を有する要素Ｅ１と、で異なる亀裂変数φが設定される。つまり、亀裂を有さない要素Ｅ０の亀裂変数φが「ｍ」と設定され、亀裂を有する要素Ｅ１の亀裂変数φが「ｍ」とは異なる「ｎ」と設定される。「ｍ」と「ｎ」とではいずれが大きくても構わない。一例として、亀裂を有さない要素Ｅ０の亀裂変数φが「０」と設定され、亀裂を有する要素Ｅ１の亀裂変数φが「１」と設定される。

【0061】

なお、構造体モデルＭ_Dに予め亀裂変数φが設定されている場合などには、本ステップＳ１３は省略しても構わない。

【0062】

（振幅負荷エネルギー設定ステップ）
ステップＳ１４において、振幅負荷エネルギー設定部１４は、ステップＳ１２で取得された構造体モデルＭ_Dの各要素Ｅに、振幅負荷エネルギーＡを設定する。なお、既に亀裂を有する要素Ｅ１では、弾性応答が生じないため、振幅負荷エネルギーＡが蓄積されない。このため、要素Ｅ１の振幅負荷エネルギーＡは「０」と設定される。

【0063】

亀裂を有さない要素Ｅ０における振幅負荷エネルギーＡは、当該要素Ｅ０を構成する材料に応じて実験的に得られる相当応力σと相当弾性ひずみεとの関係に基づいて設定される。相当弾性ひずみεは亀裂変数φに依存するため、振幅負荷エネルギーＡは亀裂変数φの関数として表現される。

【0064】

なお、構造体モデルＭ_Dに予め振幅負荷エネルギーＡが設定されている場合などには、本ステップＳ１４は省略しても構わない。

【0065】

（微分方程式生成ステップ）
ステップＳ１５において、微分方程式生成部１５は、ステップＳ１３で設定した亀裂変数φと、ステップＳ１４で設定した振幅負荷エネルギーＡと、が利用されて微分方程式を作成する。

【0066】

ステップＳ１５において生成される微分方程式としては、一例として、上記の微分方程式（４）が挙げられる。また、ステップＳ１５においては、微分方程式（４）を改良した、次の微分方程式（５）が生成されてもよい。なお、式（５）は、拡散方程式でもある。

【数5】

【0067】

微分方程式（５）には、障壁エネルギーｆ_doubの微分項のフィッティング用定数ｗ_doub、弾性エネルギーｆ_elastの微分項のフィッティング用定数ｗ_elast、及び振幅負荷エネルギーＡの項のフィッティング用定数ｗ_hcが導入されている。これにより、構造体Ｄの構成などに応じて、障壁エネルギーｆ_doubの微分項と、弾性エネルギーｆ_elastの微分項と、振幅負荷エネルギーＡの項と、のそれぞれの重み付けを最適化することが可能である。これにより、構造体Ｄに発生する亀裂を更に正確に予測可能となる。

【0068】

なお、式（５）における振幅負荷エネルギーＡは、繰り返し負荷に対応した対象材料にかかる応力振幅Δσを、図４のグラフで、応力振幅Δσ＝最大応力σ_max−最小応力σ_minと定義したときに、その負荷に対応するエネルギーである。振幅負荷エネルギーＡは、次の式（６）で表される。

【数6】

【0069】

ここで、Δεは、ひずみ振幅であり、最大ひずみε_max−最小ひずみε_minである。なお、このひずみは弾性ひずみである。

【0070】

また、式（５）の振幅負荷エネルギーＡに係るｗ_hcについて、図５を参照して説明する。

【0071】

図５は、対象材料のS-N線図を表している。S-N線図は、縦軸に応力振幅Δσ、横軸に破断までの繰り返し数Nをとって疲労試験結果をグラフで表したものである。実際には、繰り返し数Ｎによって、メカニズムの異なる低サイクル疲労（塑性疲労）領域と高サイクル疲労（弾性疲労）領域に分かれており、ｗ_hcは、高サイクル疲労領域の傾きに対応したパラメータであり、解析上は、亀裂の進行速度を決めるパラメータである。

【0072】

なお、材料の低サイクル、高サイクルの疲労特性を記述する際には、図６に示される数式である次の式（７）でフィッティングされる場合が多い。

【数7】

【0073】

式（７）の右辺の第１項は、低サイクル疲労領域の応力振幅を表し、第２項が、高サイクル疲労領域の応力振幅を表している。その際、ｗ_hcは、関数ｆを用いて、次の式（８）で表される。

【数8】

【0074】

なお、ここで、関数ｆは、線形であってもよいし、非線形であってもよい。

【0075】

異種材料が、物理結合、化学結合、粘着結合、静電結合などで結合しており、その結合力寿命が各材料のバルク内結合力よりも小さい場合、結合部に繰り返し負荷がかかると界面剥離が進行する。その際の特性は、図５および図６で示されたように、界面での疲労破断についても特性線図が得られる。

【0076】

解析モデルの界面上に対応したｗ_hcxを定義することで、高サイクルにおける剥離の進行も解析可能となる。

【0077】

さらに、図７を参照して、対象物に対して想定される負荷が、周期的な振動負荷でない（すなわち、非周期負荷である）場合についての対応を説明する。

【0078】

非周期負荷の例として、対象にかかる応力の時刻歴が、図７に示されるような場合を考える。図７の例においては、応力振幅は、時間によってバラバラであり、また、振幅のピーク間の時間の間隔もバラバラである。

【0079】

このような場合においては、解析手法上、時間軸に沿って、亀裂進行を解析していく手法となっているため、時間軸を応力の変曲点ごとに区切り、その変曲点毎の時間を区間として、その区間の応力振幅Δσを元に微分方程式を用いて解析を行い、亀裂進行判定する。１つの区間の亀裂進行判定後、次の区間の応力振幅Δσを元に微分方程式を用いて解析を行い、亀裂進行判定するという手順で処理を進めるようにする。

【0080】

例えば、図７に示されるように、σ_min1，σ_max1，σ_min2，σ_max2，σ_min3，σ_max3は、それぞれ、応力の変曲点であり、時刻t₀,t₁,t₂,t₃,t₄,t₅における応力である。σ_min1乃至σ_max1の振幅の区間（時間）Ｔ₁（t₀乃至t₁）においては、Δσ１＝σ_max1−σ_min1を元に微分方程式を用いて解析が行われる。次に、σ_max1乃至σ_min2の振幅の区間Ｔ₂（t₁乃至t₂）においては、Δσ２＝σ_max1−σ_min2を元に微分方程式を用いて解析が行われる。σ_min2乃至σ_max2の振幅の区間Ｔ₃（t₂乃至t₃）においては、Δσ３＝σ_max2−σ_min2を元に微分方程式を用いて解析が行われる。次に、σ_max2乃至σ_min3の振幅の区間Ｔ₄（t₃乃至t₄）においては、Δσ４＝σ_max2−σ_min3を元に微分方程式を用いて解析が行われる。σ_min3乃至σ_max3の振幅の区間Ｔ₅（t₄乃至t₅）においては、Δσ５＝σ_max3−σ_min3を元に微分方程式を用いて解析が行われる。

【0081】

以上のように、変曲点ごとに区切り、その区間のΔσを元に微分方程式を用いて解析するようにしたので、どのような負荷波形にも対応して亀裂を予測することができる。

【0082】

なお、微分方程式が予め生成されている場合などには、本ステップＳ１５は省略しても構わない。

【0083】

（亀裂予測ステップ）
ステップＳ１６においては、亀裂予測部１６は、ステップＳ１５により生成された微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。なお、ステップＳ１５を行わない場合には、本ステップＳ１６においては外部機器などから取得した微分方程式を計算することにより構造体Ｄに発生する亀裂が予測される。

【0084】

微分方程式を計算する際に、構造体Ｄに加わる応力をまず再現するために、構造体モデルＭ_Dに対して荷重条件を付与して応力解析を行う。図８は、構造体モデルＭ_Dに付与する荷重条件の一例を示している。図８に示す例では、構造体モデルＭ_Dにおいて、Ｘ軸方向左側の面を固定（拘束）した状態で、Ｘ軸方向右側の面に引張荷重を加える。そして、この荷重条件下で微分方程式を計算することにより、時間の経過が伴う各要素Ｅ０における亀裂変数φの変化が得られる。

【0085】

図９は、図８に示すように構造体モデルＭ_Dに荷重条件を付与した場合の、ある時刻における振幅負荷エネルギーＡの分布を示している。図９には、振幅負荷エネルギーＡの等しい等エネルギー面が破線で示されている。構造体モデルＭ_D内において、等エネルギー面が円弧状に広がっており、亀裂の先端である要素Ｅ１のＹ軸方向下面に近いほど振幅負荷エネルギーＡが大きい。

【0086】

図１０は、図８に示すように構造体モデルＭ_Dに荷重条件を付与した場合の、ある時刻における亀裂変数φの分布を示している。図１０には、亀裂変数φの等しい等亀裂変数面が破線で示されている。構造体モデルＭ_D内において、等亀裂変数面は、亀裂の先端である要素Ｅ１のＹ軸方向下面からＹ軸方向下方に延びる楕円弧状に広がっている。内側の等亀裂変数面ほど亀裂変数φが大きい。

【0087】

ステップＳ１６では、所定時間が経過した後に、亀裂変数φが「１」以上である要素Ｅ０に亀裂が発生するものと予測される。例えば、要素Ｅ１のＹ軸方向下側の３つの要素Ｅ０の亀裂変数φが「１」以上であった場合には、図１１に示すように、当該３つの要素Ｅ０に亀裂が発生したものとして当該３つの要素Ｅ０を要素Ｅ１に変更する。

【0088】

また、微分方程式の計算を進める過程において、亀裂変数φが「１」以上となった要素Ｅ０の振幅負荷エネルギーＡを順次「０」に変更していくことが好ましい。これにより、構造体モデルＭ_Dの亀裂を順次更新しながら微分方程式の計算を進めることができるため、より正確に亀裂を予測することができる。

【0089】

このように、ステップＳ１６においては、微分方程式を計算することによって、所定時間が経過した後の構造体モデルＭ_Dにおける亀裂を有する要素Ｅ１の分布が得られる。そして、構造体モデルＭ_Dにおける亀裂を有する要素Ｅ１の分布によって、構造体Ｄに発生する亀裂を予測することができる。

【0090】

［微分方程式の変形例］
ステップＳ１５で生成される微分方程式は、フェーズフィールド法の考え方に基づいて生成された微分方程式（４），（５）に限定されず、適宜変更可能である。以下、本技術で利用可能な微分方程式の変形例について説明する。

【0091】

１．構造体Ｄを形成する材料に応じたカスタマイズフェーズフィールド法の考え方に基づいて生成された微分方程式（４），（５）では、空間座標の２階微分に比例する拡散項や弾性エネルギーｆ_elastの微分項が含まれているため、多岐にわたる材料で形成された構造体Ｄに適用可能である。つまり、微分方程式（４），（５）では、高い汎用性が得られる。

【0092】

この一方で、構造体Ｄを形成する材料によっては、微分方程式（４），（５）には不要な項が含まれていることになる。したがって、ステップＳ１５では、微分方程式について、構造体Ｄを形成する材料に応じて、不要な項を排除するなどのカスタマイズをすることにより、構造体Ｄに発生する亀裂を迅速かつ的確に予測可能となる。

【0093】

本技術に係る亀裂予測方法は、金属材料や樹脂材料などの高サイクル疲労破壊しやすい材料を用いて構成された構造体Ｄの亀裂を予測可能であればよい。このため、ステップＳ１５で生成する微分方程式は、少なくとも、亀裂変数φの時間微分に比例する項と、ステップＳ１４で設定した振幅負荷エネルギーＡに比例する項と、を含んでいればよい。

【0094】

また、図１の亀裂予測装置１０は、少なくともモデル取得部１２及び亀裂予測部１６を備えていればよい。つまり、ステップＳ１１，Ｓ１３乃至Ｓ１５を実行しない場合には、亀裂予測装置１０がモデル生成部１１、亀裂変数設定部１３、振幅負荷エネルギー設定部１４、及び微分方程式生成部１５を備えていなくても構わない。また、亀裂予測装置１０は、必要に応じて、上記以外の構成を含んでいてもよい。

【0095】

以下、構造体Ｄを形成する材料の一例を挙げ、当該材料に応じてカスタマイズされた微分方程式を例示する。なお、構造体Ｄを形成する材料は、以下に挙げるものに限定されず、任意の材料であってよい。また、各材料に対応する微分方程式も、以下に例示するものに限定されず、任意にカスタマイズ可能である。

【0096】

（ａ）脆性破壊しにくい材料
構造体Ｄを形成する材料が脆性破壊しにくい場合、例えば、振幅負荷エネルギーＡの項以外の項を排除し、高サイクル疲労破壊のみを考慮した、次の微分方程式（９）を用いることができる。

【数9】

微分方程式（９）には、亀裂変数φの時間微分の項、及び振幅負荷エネルギーＡの項のみが含まれる。なお、振幅負荷エネルギーＡの項には、フィッティング用定数ｗ_hcが含まれていてもよい。このように、振幅負荷エネルギーＡの項以外の項を排除して簡単化された微分方程式（９）を用いることにより、計算負荷を大幅に小さくすることができる。

【0097】

（ｂ）弾性率Ｂが異方性を有する材料
構造体Ｄを形成する材料の弾性率Ｂが異方性を有する場合、例えば、弾性率Ｂの異方性を考慮した、次の微分方程式（１０）を用いることができる。

【数10】

【0098】

微分方程式（１０）中、系のエネルギーＦ_sysは、次の式（１１）で表される。

【数11】

【0099】

式（１１）中、勾配エネルギーｆ_gradは、次の式（１２）で表される。式（１１）中、弾性エネルギーｆ_elastは、次の式（１３）で表される。

【数12】

式（１２）中、κは、材料定数を示す。

【0100】

【数13】

式（１３）中、εは垂直ひずみを示す。

【0101】

式（１３）では、弾性率Ｂ_ijklをテンソル（行列）として扱うことにより、弾性率Ｂの異方性を適切に予測結果に反映させることができる。このため、微分方程式（１０）を計算することによって、弾性率Ｂが異方性を有する材料で形成された構造体Ｄに発生する亀裂を的確に予測することができる。

【0102】

なお、構造体Ｄを形成する材料の弾性率Ｂが等方性を有する場合には、弾性エネルギーｆ_elastは、式（１３）に代えて、次の式（１４）で表すことができる。

【数14】

式（１４）中、νはポアソン比を示し、γはせん断ひずみを示す。

【0103】

また、式（１４）中の弾性率Ｂは、例えば、図１２に示すような亀裂変数φに依存する関数とすることができる。図１２に示す関数では、亀裂変数φが大きくなるにつれて、弾性率Ｂが小さくなる。つまり、図１２に示す関数では、振幅負荷エネルギーＡの蓄積に伴って構造体Ｄを形成する材料の弾性が低下することを表現することができる。

【0104】

（ｃ）靭性値が異方性を有する材料
構造体Ｄを形成する材料の靭性値が異方性を有する場合、例えば、拡散項の係数を亀裂変数φの勾配の関数、すなわち界面の法線方向の関数とする、次の微分方程式（１５）を用いることができる。

【数15】

【0105】

微分方程式（１５）における系のエネルギーＦ_sysは、次の式（１６）で表される。

【数16】

【0106】

式（１６）における勾配エネルギーｆ_gradは、次の式（１７）で表される。式（１６）における弾性エネルギーｆ_elastは、次の式（１９）で表される。

【数17】

式（１７）におけるκは、材料定数を示し、次の式（１８）で表される。

【0107】

【数18】

式（１８）において、ａは異方性関数である。

【0108】

【数19】

式（１９）において、νはポアソン比を示し、γはせん断ひずみを示す。

【0109】

式（１８）を用いることによって、界面の向きに応じて拡散係数を変えることができ、亀裂の進行のしやすさを方向によって変化させることができる。これにより、靭性値の異方性を予測結果に適切に反映させることができる。このため、微分方程式（１５）を計算することによって、靭性値が異方性を有する材料で形成された構造体Ｄに発生する亀裂を的確に予測することができる。

【0110】

（ｄ）脆性破壊と高サイクル疲労破壊とが同時に進行する材料
構造体Ｄが脆性破壊と高サイクル疲労破壊とが同時に進行する材料によって形成されている場合、構造体Ｄには脆性破壊と高サイクル疲労破壊との組み合わせによる亀裂が発生する。脆性破壊と高サイクル疲労破壊との組み合わせによる亀裂を予測するために、例えば、微分方程式（２０）を用いることができる。

【数20】

【0111】

微分方程式（２０）における系のエネルギーＦ_sysは、次の式（２１）で表される。

【数21】

【0112】

式（２１）における勾配エネルギーｆ_gradは、次の式（２２）で表される。式（２１）における弾性エネルギーｆ_elastは、次の式（２３）で表される。

【数22】

式（２２）において、κは材料定数を示す。

【0113】

【数23】

式（２３）において、νはポアソン比を示し、εは垂直ひずみを示し、γはせん断ひずみを示す。

【0114】

微分方程式（２０）では、弾性エネルギーｆ_elastの解放率によって脆性破壊を解析することができるとともに、振幅負荷エネルギーＡの蓄積によって高サイクル疲労破壊を解析することもできる。したがって、微分方程式（２０）を計算することによって、構造体Ｄに発生する脆性破壊と高サイクル疲労破壊との組み合わせによる亀裂を予測可能となる。

【0115】

２．界面の安定化
構造体Ｄにおける亀裂を良好に表現するために、構造体モデルＭ_Dにおいて亀裂を有さない要素Ｅ０と亀裂を有する要素Ｅ１との界面を安定化することが好ましい。つまり、亀裂を有さない要素Ｅ０の亀裂変数φが「０」近傍の値であり、亀裂を有する要素Ｅ１の亀裂変数が「１」近傍の値であり、いずれの要素Ｅもなるべく亀裂変数が「０」と「１」との中間の値である状態とならないことが好ましい。

【0116】

構造体モデルＭ_Dにおいて亀裂を有さない要素Ｅ０と亀裂を有する要素Ｅ１との界面を安定化するために、例えば、次の微分方程式（２４）を用いることができる。

【数24】

【0117】

微分方程式（２４）における系のエネルギーＦ_sysは、次の式（２５）で表される。

【数25】

【0118】

式（２５）中、障壁エネルギーｆ_doubは、次の式（２６）で表される。式（２５）中、弾性エネルギーｆ_elastは、次の式（２７）で表される。

【数26】

式（２６）において、κは材料定数を示す。

【0119】

【数27】

式（２７）におけるｅ_doubはエネルギー障壁を示す。

【0120】

上述した式（２７）は、図１３に示すようなダブルウェル関数である。つまり、障壁エネルギーｆ_doubが亀裂変数φ＝０，１において極小値をとる。このため、亀裂変数φは、「０」又は「１」近傍の値をとりやすく、「０」と「１」との中間の値をとりにくい。これにより、構造体モデルＭ_Dにおいて亀裂を有さない要素Ｅ０と亀裂を有する要素Ｅ１との界面が安定化する。

【0121】

また、式（２７）では、構造体Ｄを形成する材料における亀裂の発生しにくさをエネルギー障壁ｅ_doubによって表現することが可能である。つまり、亀裂が発生しやすい材料の場合にはエネルギー障壁ｅ_doubを低くし、亀裂が発生しにくい材料の場合にはエネルギー障壁ｅ_doubを高くすることができる。

【0122】

以上により、本技術の第１の実施形態によれば、構造体に発生する延性破壊のうち、高サイクル疲労（弾性疲労）による亀裂を迅速に予測することができる。

【0123】

なお、第１の実施形態は、構造体に発生する疲労破壊のうち、高負荷の塑性ひずみに入らない疲労（低サイクル疲労）を考慮しない場合、低負荷の高サイクル疲労（弾性疲労）しか考慮しない場合について説明した。ただし、材料によっては、構造体に発生する疲労破壊として、低負荷の高サイクル疲労（弾性疲労）だけでなく、高負荷の塑性ひずみ（低サイクル疲労）を考慮が必要な場合があり、この場合について、次の第２の実施形態で説明する。

【0124】

＜２．第２の実施形態＞
第２の実施形態に係る技術として、低負荷の高サイクル疲労（弾性疲労）および高負荷の低サイクル疲労（塑性疲労）を含む構造体に発生する疲労破壊を予測可能な方法について説明する。なお、以降は、第１の実施形態との差分に焦点をあてて説明する。

【0125】

本技術の第２の実施形態においては、微分方程式（４）に、塑性散逸エネルギーｆ_plastの項を導入した、次の微分方程式（２８）を利用する。

【数28】

【0126】

微分方程式（２８）は、微分方程式（４）において、塑性散逸エネルギーｆ_plastの項を追加した式である。また、微分方程式（２８）は、微分方程式（３）において、振幅負荷エネルギーＡの項を追加した式であるともいえる。すなわち、微分方程式（２８）は、微分方程式（２）に対して、繰り返し負荷に対応する低サイクル疲労解析を行うための塑性散逸エネルギーｆ_plastと高サイクル疲労解析を行うための振幅負荷エネルギーＡの項を追加したものである。低負荷の高サイクル疲労（弾性疲労）および高負荷の塑性ひずみ（低サイクル疲労）の両方含む場合は、この式（２８）を用いて、予測が行われる。

【0127】

微分方程式（２８）において、塑性散逸エネルギーｆ_plastおよび振幅負荷エネルギーＡは、弾性エネルギーｆ_elastとは異なり、微分項ではない。これは、弾性エネルギーｆ_elastは時間とともに解放されるのに対し、塑性散逸エネルギーｆ_plastおよび振幅負荷エネルギーＡは時間とともに蓄積されるためである。微分方程式（２８）では、塑性散逸エネルギーｆ_plastおよび振幅負荷エネルギーＡを微分項としないことにより、塑性散逸エネルギーｆ_plastおよび振幅負荷エネルギーＡの蓄積を表現することができる。

【0128】

このように、微分方程式（２８）には、弾性エネルギーｆ_elastの解放率を表現する弾性エネルギーｆ_elastの微分項と、塑性散逸エネルギーｆ_plastの蓄積を表現する塑性散逸エネルギーｆ_plastの項と、振幅負荷エネルギーＡの蓄積を表現する振幅負荷エネルギーＡの項と、が含まれる。したがって、微分方程式（２８）を計算することにより、脆性破壊及び疲労破壊（高低サイクル疲労含む）の双方を考慮して、亀裂を予測することが可能となる。

【0129】

このため、本技術では、金属材料や樹脂材料などの低サイクル疲労（塑性疲労）および高サイクル疲労（弾性疲労）を含む疲労破壊しやすい材料を用いて構成された構造体Ｄに発生する亀裂を正確に予測可能である。また、本技術に係る亀裂予測方法でも、フェーズフィールド法の考え方による亀裂予測方法と同様に、複数種類の材料の界面を跨ぐ亀裂を予測可能であるため、複数の材料で構成された構造体Ｄに発生する亀裂を迅速に予測可能である。更に、本技術に係る亀裂予測方法でも、フェーズフィールド法の考え方による亀裂予測方法と同様に、亀裂の形状に制約がないため、高い汎用性が得られる。

【0130】

［亀裂予測装置の詳細］
図１４は、本技術の第２の実施形態に係る亀裂予測装置（情報処理装置）１０１の構成を示すブロック図である。亀裂予測装置１０１は、モデル生成部１１と、モデル取得部１２と、亀裂変数設定部１３と、エネルギー設定部１１１と、微分方程式生成部１１２と、亀裂予測部１６と、を備える。

【0131】

図１４の亀裂予測装置１０１は、振幅負荷エネルギー設定部１４が、エネルギー設定部１１１に入れ替わった点と、微分方程式生成部１５が微分方程式生成部１１２に入れ替わった点とが、図１の亀裂予測装置１０と異なっている。図１４の亀裂予測装置１０１は、モデル生成部１１と、モデル取得部１２と、亀裂変数設定部１３と、亀裂予測部１６とを備える点が図１の亀裂予測装置１０と共通している。

【0132】

すなわち、エネルギー設定部１１１は、構造体モデルＭ_Dの各要素Ｅに、振幅負荷エネルギーＡと塑性散逸エネルギーｆ_plastとをそれぞれ設定する。

【0133】

微分方程式生成部１１２は、亀裂変数設定部１３により設定された亀裂変数φと、エネルギー設定部１１１により設定された振幅負荷エネルギーＡと、エネルギー設定部１１１により設定された塑性散逸エネルギーｆ_plastとを利用して微分方程式を作成する。

【0134】

したがって、図１４の亀裂予測部１６は、微分方程式生成部１１２により生成された微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。

【0135】

［亀裂予測動作の例］
図１５のフローチャートは、図１４の亀裂予測装置１０１の亀裂予測処理を説明するフローチャートである。以下、図１５に沿って、第２の実施形態に係る亀裂予測処理について説明する。なお、図１５のステップＳ１１１乃至Ｓ１１３、およびＳ１１６は、図２のステップＳ１１乃至Ｓ１３、およびＳ１６と基本的に同様の処理を行うため、その詳細な説明は省略する。

【0136】

（モデル生成ステップ）
ステップＳ１１１において、モデル生成部１１は、構造体Ｄの構成を再現したモデル（構造体モデル）Ｍ_Dを生成する。

【0137】

（モデル取得ステップ）
ステップＳ１１２において、モデル取得部１２は、ステップＳ１１１で生成された構造体モデルＭ_Dを取得する。

【0138】

（亀裂変数設定ステップ）
ステップＳ１１３において、亀裂変数設定部１３は、ステップＳ１１２で取得された構造体モデルＭ_Dの各要素Ｅに亀裂の有無を表現する亀裂変数φを設定する。

【0139】

（振幅負荷エネルギー設定ステップ）
ステップＳ１１４において、エネルギー設定部１１１は、ステップＳ１１２で取得された構造体モデルＭ_Dの各要素Ｅに、塑性散逸エネルギーｆ_plastと振幅負荷エネルギーＡとを設定する。なお、ここでは、図２のステップＳ１４との差分である、塑性散逸エネルギーｆ_plastの設定についてのみ説明する。

【0140】

すなわち、既に亀裂を有する要素Ｅ１では、塑性変形が生じないため、塑性散逸エネルギーｆ_plastが蓄積されない。このため、要素Ｅ１の塑性散逸エネルギーｆ_plastは「０」と設定される。

【0141】

亀裂を有さない要素Ｅ０における塑性散逸エネルギーｆ_plastは、当該要素Ｅ０を構成する材料に応じて実験的に得られる相当応力σと相当塑性ひずみε_pとの関係に基づいて設定される。相当塑性ひずみε_pは亀裂変数φに依存するため、塑性散逸エネルギーｆ_plastは、亀裂変数φの関数として表現される。

【0142】

図１６は、ステップＳ１１４において要素Ｅ０に設定される塑性散逸エネルギーｆ_plastの表現方法の一例を示す図である。図１６は、構造体Ｄを形成する材料から得られる相当応力−相当塑性ひずみ線図の一例を示している。図１６では、縦軸が相当応力σを示し、横軸が相当塑性ひずみε_pを示している。また、図１６には降伏応力σ_Yが示されている。

【0143】

図１６に示す相当応力−相当塑性ひずみ線図が得られる材料は、相当応力σが降伏応力σ_Y未満の領域で弾性変形し、相当応力σが降伏応力σ_Y以上の領域で塑性変形する。塑性散逸エネルギーｆ_plastは、相当応力σが降伏応力σ_Y以上のとき、当該材料の塑性変形により主に熱エネルギーとして散逸されるエネルギーを表現するものである。

【0144】

塑性散逸エネルギーｆ_plastは、例えば、図１６Ａ及び図１６Ｂに斜線で示す領域の面積として定義することが可能である。図１６Ａに斜線で示す領域の面積は、相当応力σを相当塑性ひずみε_pの微小増分で積分した量を利用して、例えば、次の式（２９）によって算出可能である。

【数29】

【0145】

また、図１６Ｂに斜線で示す領域の面積は、相当応力σと降伏応力σ_Yとの差と、相当塑性ひずみε_pと、の積を利用して、例えば、次の式（３０）によって算出可能である。

【数30】

【0146】

ただし、式（３０）において、相当応力σが降伏応力σ_Yより小さい場合の塑性散逸エネルギーｆ_plastをゼロとする。

【0147】

塑性散逸エネルギーｆ_plastの表現方法は、構造体Ｄを形成する材料や物理現象などに応じて、正確に亀裂を予測できるように使い分けることが可能である。なお、塑性散逸エネルギーｆ_plastを表現する関数は、式（２９）及び式（３０）に限定されず、相当応力σと相当塑性ひずみε_pとの関係に基づいて適宜作成可能である。

【0148】

なお、構造体モデルＭ_Dに振幅負荷エネルギーＡおよび塑性散逸エネルギーｆ_plastが予め設定されている場合などには、本ステップＳ１１４は省略しても構わない。

【0149】

（微分方程式生成ステップ）
ステップＳ１１５において、微分方程式生成部１１２は、ステップＳ１１３で設定した亀裂変数φと、ステップＳ１１４で設定された振幅負荷エネルギーＡと、ステップＳ１１４で設定された塑性散逸エネルギーｆ_plastが利用されて微分方程式を作成する。

【0150】

ステップＳ１１５において生成される微分方程式としては、一例として、上記の微分方程式（２８）が挙げられる。また、ステップＳ１１５においては、微分方程式（２８）を改良した、次の微分方程式（３１）が生成されてもよい。なお、式（３１）は、拡散方程式でもある。

【数31】

【0151】

微分方程式（３１）には、障壁エネルギーｆ_doubの微分項のフィッティング用定数ｗ_doub、弾性エネルギーｆ_elastの微分項のフィッティング用定数ｗ_elast、塑性散逸エネルギーｆ_plastの項のフィッティング用定数ｗ_plast、及び振幅負荷エネルギーＡの項のフィッティング用定数ｗ_hcが導入されている。これにより、構造体Ｄの構成などに応じて、障壁エネルギーｆ_doubの微分項と、弾性エネルギーｆ_elastの微分項と、塑性散逸エネルギーｆ_plastの項と、振幅負荷エネルギーＡの項と、のそれぞれの重み付けを最適化することが可能である。これにより、構造体Ｄに発生する亀裂を更に正確に予測可能となる。

【0152】

ここで、第１の実施形態において、振幅負荷エネルギーＡについて上述された微分方程式の変形例、例えば、式（９）、式（１０）、式（１５）、式（２０）、および式（２４）などは、振幅負荷エネルギーＡを塑性散逸エネルギーｆ_plastに変えて、同様に適用される。

【0153】

なお、微分方程式が予め生成されている場合などには、本ステップＳ１１５は省略しても構わない。

【0154】

（亀裂予測ステップ）
ステップＳ１１６において、亀裂予測部１６は、ステップＳ１１５で生成された微分方程式を計算することにより、構造体Ｄに発生する亀裂を予測する。なお、ステップＳ１１５を行わない場合には、本ステップＳ１１６においては外部機器などから取得した微分方程式を計算することにより構造体Ｄに発生する亀裂が予測される。

【0155】

以上により、本技術の第２の実施形態によれば、構造体に発生する疲労破壊（低サイクル疲労（塑性疲労）および高サイクル疲労（弾性疲労））による亀裂を迅速に予測することができる。

【0156】

すなわち、本技術の第２の実施形態によれば、脆性亀裂、低サイクル疲労（塑性疲労）による亀裂、高サイクル繰り返し疲労による亀裂など、すべての亀裂現象が一気通貫で解析可能となる。すなわち、本技術によれば、構造体Ｄに発生する亀裂を迅速かつより正確に予測可能となる。

【0157】

[その他の実施形態]
以上、本技術の実施形態について説明したが、本技術は上述の実施形態にのみ限定されるものではなく、本技術の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。

【0158】

例えば、上記実施形態では構造体モデルＭ_Dの要素Ｅを１次要素としたが、必要に応じて、構造体モデルＭ_Dの要素Ｅを２次要素としてもよい。この場合、構造体モデルＭ_Dの各要素Ｅ内における亀裂変数φの分布を考慮することができるため、構造体Ｄの亀裂を更に正確に予測可能となる。

【0159】

以上のように、本技術によれば、亀裂有無を表現する場の亀裂変数φを振幅負荷エネルギーに比例する項をソース項とした拡散方程式で時間発展させた方程式を数値的に解くようにした。

【0160】

これにより、界面をまたぐような亀裂にも対応し、亀裂形状も自由に変化できるフェーズフィールド法の利点を生かしながら、金属／樹脂の弾性域での高サイクル疲労破壊を予測可能になる。

【0161】

さらに、本技術によれば、亀裂有無を表現する場の亀裂変数φを塑性散逸エネルギーに比例する項と、亀裂有無を表現する場の亀裂変数φをソース項とした拡散方程式で時間発展させた方程式を数値的に解くようにした。

【0162】

これにより、界面をまたぐような亀裂にも対応し、亀裂形状も自由に変化できるフェーズフィールド法の利点を生かしながら、金属／樹脂の塑性域での低サイクル繰り返し疲労破壊および弾性域での高サイクル疲労破壊を予測可能になる。これにより、亀裂をより正確に予測することができる。

【0163】

＜３．コンピュータ＞
上述した一連の処理は、ハードウエアにより実行させることもできるし、ソフトウエアにより実行させることもできる。一連の処理をソフトウエアにより実行する場合には、そのソフトウエアを構成するプログラムが、コンピュータにインストールされる。ここでコンピュータには、専用のハードウエアに組み込まれているコンピュータや、各種のプログラムをインストールすることで、各種の機能を実行することが可能な、例えば汎用のパーソナルコンピュータ等が含まれる。

【0164】

図１７は、上述した一連の処理をプログラムにより実行するコンピュータのハードウエアの構成例を示すブロック図である。

【0165】

図１７に示されるコンピュータにおいて、CPU（Central Processing Unit）３０１、ROM（Read Only Memory）３０２、RAM（Random Access Memory）３０３は、バス３０４を介して相互に接続されている。

【0166】

バス３０４にはまた、入出力インタフェース３０５も接続されている。入出力インタフェース３０５には、入力部３０６、出力部３０７、記憶部３０８、通信部３０９、およびドライブ３１０が接続されている。

【0167】

入力部３０６は、例えば、キーボード、マウス、マイクロホン、タッチパネル、入力端子などよりなる。出力部３０７は、例えば、ディスプレイ、スピーカ、出力端子などよりなる。記憶部３０８は、例えば、ハードディスク、RAMディスク、不揮発性のメモリなどよりなる。通信部３０９は、例えば、ネットワークインタフェースよりなる。ドライブ３１０は、磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリなどのリムーバブルメディア３１１を駆動する。

【0168】

以上のように構成されるコンピュータでは、CPU３０１およびバス３０４を介して、RAM３０３にロードして実行することにより、上述した一連の処理が行われる。RAM３０３にはまた、CPU３０１が各種の処理を実行する上において必要なデータなども適宜記憶される。

【0169】

コンピュータ（CPU３０１）が実行するプログラムは、例えば、パッケージメディア等としてのリムーバブルメディア３１１に記録して適用することができる。その場合、プログラムは、リムーバブルメディア３１１をドライブ３１０に装着することにより、入出力インタフェース３０５を介して、記憶部３０８にインストールすることができる。

【0170】

また、このプログラムは、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線または無線の伝送媒体を介して提供することもできる。その場合、プログラムは、通信部３０９で受信し、記憶部３０８にインストールすることができる。

【0171】

その他、このプログラムは、ROM３０２や記憶部３０８に、あらかじめインストールしておくこともできる。

【0172】

また、本技術の実施の形態は、上述した実施の形態に限定されるものではなく、本技術の要旨を逸脱しない範囲において種々の変更が可能である。

【0173】

なお、以上において、１つの装置（または処理部）として説明した構成を分割し、複数の装置（または処理部）として構成するようにしてもよい。逆に、以上において複数の装置（または処理部）として説明した構成をまとめて１つの装置（または処理部）として構成されるようにしてもよい。また、各装置（または各処理部）の構成に上述した以外の構成を付加するようにしてももちろんよい。さらに、システム全体としての構成や動作が実質的に同じであれば、ある装置（または処理部）の構成の一部を他の装置（または他の処理部）の構成に含めるようにしてもよい。つまり、本技術は、上述した実施の形態に限定されるものではなく、本技術の要旨を逸脱しない範囲において種々の変更が可能である。

【0174】

以上、添付図面を参照しながら本開示の好適な実施形態について詳細に説明したが、本開示はかかる例に限定されない。本開示の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本開示の技術的範囲に属するものと了解される。

【0175】

なお、本技術は以下のような構成も取ることができる。
（１） 所定の構造体に対応する構造体モデルを取得するモデル取得部と、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
を備える情報処理装置。
（２） 前記振幅負荷エネルギーは、繰り返し負荷時の応力の最大値と最小値の差分と、歪みの最大値と最小値の差分とが積算されて設定される
前記（１）に記載の情報処理装置。
（３） 前記振幅負荷エネルギーに係る定数は、対応する材料のS-N線図の高サイクル疲労領域の特性線の傾きに応じた値として設定される
前記（１）または（２）に記載の情報処理装置。
（４） 前記振幅負荷エネルギーに係る定数は、対応する材料界面のS-N線図の高サイクル疲労剥離領域の特性線の傾きに応じた値として設定される
前記（１）または（２）に記載の情報処理装置。
（５） 前記亀裂予測部は、非周期的な繰り返し負荷に対して、応力時間変化の変曲点毎に解析区間を区切り、区切られた解析区間毎に前記構造体に発生する亀裂を予測する
前記（１）乃至（４）のいずれかに記載の情報処理装置。
（６） 前記亀裂予測部は、前記構造体モデルの各位置に設定され、塑性変形時に散逸されるエネルギーを前記亀裂変数を利用して表現する塑性散逸エネルギーに比例する項をさらに含む前記微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
前記（１）乃至（５）のいずれかに記載の情報処理装置。
（７） 前記塑性散逸エネルギーは、相当応力を相当塑性ひずみの微小増分で積分した量を利用して設定される
前記（６）に記載の情報処理装置。
（８） 前記塑性散逸エネルギーは、相当応力と降伏応力との差と、相当塑性ひずみとの積を利用して設定され、相当応力が降伏応力より小さい場合にゼロとされる
前記（６）または（７）に記載の情報処理装置。
（９） 前記微分方程式は、空間座標の２階微分に比例する拡散項をさらに含む
前記（６）乃至（８）のいずれかに記載の情報処理装置。
（１０） 所定の構造体に対応する構造体モデルを取得し、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
情報処理方法。
（１１） 所定の構造体に対応する構造体モデルを取得するモデル取得部と、
前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
して、コンピュータを機能させるプログラム。

【符号の説明】

【0176】

１０ 亀裂予測装置（情報処理装置）， １１ モデル生成部， １２ モデル取得部，１３ 亀裂変数設定部， １４ 振幅負荷エネルギー設定部，１５ 微分方程式生成部， １６ 亀裂予測部，１１０ 亀裂予測装置， １１１ エネルギー設定部， １１２ 微分方程式生成部 ，Ｍ_D 構造体モデル， Ｅ，Ｅ０，Ｅ１ 要素

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】