JP2019003301A - 情報処理装置、情報処理方法、及びプログラム - Google Patents

情報処理装置、情報処理方法、及びプログラム Download PDF

Info

Publication number
JP2019003301A
JP2019003301A JP2017115657A JP2017115657A JP2019003301A JP 2019003301 A JP2019003301 A JP 2019003301A JP 2017115657 A JP2017115657 A JP 2017115657A JP 2017115657 A JP2017115657 A JP 2017115657A JP 2019003301 A JP2019003301 A JP 2019003301A
Authority
JP
Japan
Prior art keywords
crack
energy
differential equation
term
variable
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2017115657A
Other languages
English (en)
Inventor
本郷　一泰
Kazuyasu Hongo
一泰 本郷
嵩明 平野
Takaaki Hirano
嵩明 平野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sony Semiconductor Solutions Corp
Original Assignee
Sony Semiconductor Solutions Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sony Semiconductor Solutions Corp filed Critical Sony Semiconductor Solutions Corp
Priority to JP2017115657A priority Critical patent/JP2019003301A/ja
Priority to PCT/JP2018/020652 priority patent/WO2018230339A1/ja
Priority to US16/619,123 priority patent/US11531801B2/en
Publication of JP2019003301A publication Critical patent/JP2019003301A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0033—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining damage, crack or wear
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/32—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/006—Crack, flaws, fracture or rupture
    • G01N2203/0062—Crack or flaws
    • G01N2203/0064—Initiation of crack
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/006—Crack, flaws, fracture or rupture
    • G01N2203/0062—Crack or flaws
    • G01N2203/0066—Propagation of crack
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0073—Fatigue
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/0202—Control of the test
    • G01N2203/0212—Theories, calculations
    • G01N2203/0214—Calculations a priori without experimental data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

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

Description

本開示は、情報処理装置、情報処理方法、及びプログラムに関し、特に、亀裂をより正確に予測することができる情報処理装置、情報処理方法、及びプログラムに関する。
一般的に、半導体デバイス等の各種構造体には、その製造過程などにおいて機械的応力等の様々な応力が加わる。構造体にこのような応力が加わると、その内部に亀裂が発生することがある。このような亀裂の発生を防止するために、構造体に予め発生する可能性のある亀裂を予測する技術が用いられる。
特許文献１には、亀裂を含む断面上の異種材料の面積比に応じて各材料の応力とひずみの関係を単一の材料の応力のひずみの関係に置き換えたときの材料定数を用いて求めた破壊力学パラメータＪ積分値が亀裂先端の材料の弾塑性破壊靭性値に等しくなるときの荷重を破壊荷重とする技術が開示されている。
特許文献２には、複数種類の材料の界面を跨ぐ亀裂を予測可能な技術が開示されている。特許文献２に係る技術では、構造体の内部で仮想的に亀裂を進行させたときのエネルギー解放率を算出し、エネルギー解放率の大きい方向に亀裂が進行するものと予測する。
特開２００９−１６００２８号公報 特開２０１１−２０４０８１号公報
近年、半導体デバイス等の構造体の多様化に伴い、構造体を構成する材料として金属材料や樹脂材料が広く用いられるようになってきている。したがって、金属材料や樹脂材料を用いて構成された構造体に発生する可能性のある亀裂を正確に予測可能な技術が要求されている。
特許文献１に係る技術では、界面をまたぐ亀裂に対応することが難しく、疲労破壊への対応が難しかった。特許文献２に係る技術では、計算負荷が高く、予め定めた特定面以外の面を解析することが困難であった。
本開示は、このような状況に鑑みてなされたものであり、亀裂をより正確に予測することができるものである。
本技術の一側面の情報処理装置は、所定の構造体に対応する構造体モデルを取得するモデル取得部と、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部とを備える。
前記振幅負荷エネルギーは、繰り返し負荷時の応力の最大値と最小値の差分と、歪みの最大値と最小値の差分とが積算されて設定される。
前記振幅負荷エネルギーに係る定数は、対応する材料のS-N線図の高サイクル疲労領域の特性線の傾きに応じた値として設定される。
前記振幅負荷エネルギーに係る定数は、対応する材料界面のS-N線図の高サイクル疲労剥離領域の特性線の傾きに応じた値として設定される。
前記亀裂予測部は、非周期的な繰り返し負荷に対して、応力時間変化の変曲点毎に解析区間を区切り、区切られた解析区間毎に前記構造体に発生する亀裂を予測することができる。
前記亀裂予測部は、前記構造体モデルの各位置に設定され、塑性変形時に散逸されるエネルギーを前記亀裂変数を利用して表現する塑性散逸エネルギーに比例する項をさらに含む前記微分方程式を計算することにより、前記構造体に発生する亀裂を予測することができる。
前記塑性散逸エネルギーは、相当応力を相当塑性ひずみの微小増分で積分した量を利用して設定される。
前記塑性散逸エネルギーは、相当応力と降伏応力との差と、相当塑性ひずみとの積を利用して設定され、相当応力が降伏応力より小さい場合にゼロとされる。
前記微分方程式は、空間座標の２階微分に比例する拡散項をさらに含むことができる。
本技術の一側面の情報処理方法は、所定の構造体に対応する構造体モデルを取得し、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する。
本技術の一側面のプログラムは、所定の構造体に対応する構造体モデルを取得するモデル取得部と、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部として、コンピュータを機能させる。
本技術の一側面においては、所定の構造体に対応する構造体モデルが取得される。そして、前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂が予測される。
本技術によれば、高サイクル疲労破壊を予測することができる。この結果、亀裂をより正確に予測することができる。
なお、本明細書に記載された効果は、あくまで例示であり、本技術の効果は、本明細書に記載された効果に限定されるものではなく、付加的な効果があってもよい。
本技術の一実施形態に係る亀裂予測装置（情報処理装置）の構成例を示すブロック図である。 図１の亀裂予測装置の亀裂予測処理について説明するフローチャートである。 亀裂予測処理において生成される構造体モデルの一例を示す図である。 繰り返し負荷に対応した対象材料にかかる応力振幅を定義するための図である。 対象材料のS-N線図を表す図である。 材料の低サイクル、高サイクルの疲労特性を記述する際の式を説明する図である。 対象物に対して想定される負荷が、周期的な振動負荷でない場合の例を説明する図である。 構造体モデルに加える荷重条件の一例を説明するための図である。 構造体モデルにおける振幅負荷エネルギーの分布の一例を示す図である。 構造体モデルにおける亀裂変数の分布の一例を示す図である。 構造体モデルにおいて予測される亀裂の一例を示す図である。 弾性率と亀裂変数との関係の一例を示す図である。 障壁エネルギーと亀裂変数との関係の一例を示す図である。 本技術の他の実施形態に係る亀裂予測装置（情報処理装置）の構成例を示すブロック図である。 図１４の亀裂予測装置の亀裂予測処理について説明するフローチャートである。 亀裂予測処理において設定される塑性散逸エネルギーの表現方法の一例を示す図である。 コンピュータの構成例を示すブロック図である。
以下、本開示を実施するための形態（以下実施の形態とする）について説明する。なお、説明は以下の順序で行う。また、図面には、適宜相互に直交するＸ軸、Ｙ軸、及びＺ軸が示されている。Ｘ軸、Ｙ軸、及びＺ軸は全図において共通である。
０．概要
１．第１の実施形態
２．第２の実施形態
３.コンピュータ
＜０．概要＞
［亀裂予測方法の概要］
本技術に係る亀裂予測方法（情報処理方法）の概要について説明する。本技術に係る亀裂予測方法では、フェーズフィールド（Phase Field）法の考え方を応用して構造体Ｄに発生する亀裂を予測する。まず、本技術に関連するフェーズフィールド法の考え方による亀裂予測方法について説明する。
（フェーズフィールド法の考え方による亀裂予測方法）
構造体Ｄ内のエネルギーＦは、障壁エネルギーｆdoub、勾配エネルギーｆgrad、及び弾性エネルギーｆelastを利用して、式（１）で表される。
Figure 2019003301
フェーズフィールド法の考え方によると、式（１）から、次の微分方程式（２）を導出することができる。
Figure 2019003301
微分方程式（２）の左辺は、易動度Ｍの逆数と、亀裂の発生の有無を表現する亀裂変数φの時間微分と、の積で構成されている。微分方程式（２）の右辺は、空間座標の２階微分の拡散項∇（ξ∇φ）と、障壁エネルギーｆdoubの微分項と、弾性エネルギーｆelastの微分項と、で構成されている。微分方程式（２）では、弾性エネルギーｆelastの微分項によって、弾性エネルギーｆelastの解放率が表現されている。
微分方程式（２）の計算を行うために、まず構造体Ｄの各位置にそれぞれ亀裂変数φが設定される。より詳細には、亀裂を有さない位置と、亀裂を有する位置と、で異なる亀裂変数φが設定される。例えば、亀裂を有さない位置の亀裂変数φが「０」と設定され、亀裂を有する位置の亀裂変数φが「１」と設定される。
そして、微分方程式（２）の計算を進めると、時間の経過に伴い、亀裂変数φを「０」と設定した位置のうち、亀裂変数φが「１」以上となる位置が現れる。フェーズフィールド法の考え方による亀裂予測方法では、所定時間が経過した後に、亀裂変数φが「１」以上である位置に亀裂が発生しているものと予測することができる。
フェーズフィールド法の考え方による亀裂予測方法では、微分方程式（２）を計算することにより、亀裂を迅速に予測可能である。また、フェーズフィールド法の考え方による亀裂予測方法では、複数種類の材料の界面を跨ぐ亀裂を予測可能であるため、複数の材料で構成された構造体Ｄに発生する亀裂を予測可能である。更に、フェーズフィールド法の考え方による亀裂予測方法では、亀裂の形状に制約がないため、高い汎用性が得られる。
フェーズフィールド法の考え方による亀裂予測方法では、微分方程式（２）に含まれる弾性エネルギーｆelastの微分項で表現される弾性エネルギーｆelastの解放率によって、脆性破壊による亀裂を予測することができる。しかし、フェーズフィールド法の考え方による亀裂予測方法では、微分方程式（２）に塑性変形に対応する項が含まれていないため、塑性変形を伴う疲労破壊による亀裂を予測することができない。
したがって、フェーズフィールド法の考え方を利用した亀裂予測方法では、金属材料や樹脂材料などの疲労破壊しやすい材料を用いて構成された構造体Ｄに発生する亀裂を正確に予測することは困難である。
そこで、本技術の発明者は、フェーズフィールド法の考え方を応用し、塑性変形時に主に熱として散逸されるエネルギー（以下、「塑性散逸エネルギーｆplast」と称する）を含む項を微分方程式（２）に導入することにより、疲労破壊による亀裂を予測可能となることを見出した。以下、本技術に係るフェーズフィールド法の考え方を応用した亀裂予測方法について説明する。
（フェーズフィールド法の考え方を応用した亀裂予測方法）
本技術に係るフェーズフィールド法の考え方を応用した亀裂予測方法では、微分方程式（２）に塑性散逸エネルギーｆplastの項を導入した、次の微分方程式（３）を利用する。
Figure 2019003301
微分方程式（３）において、塑性散逸エネルギーｆplastは、弾性エネルギーｆelastとは異なり、微分項ではない。これは、弾性エネルギーｆelastは時間とともに解放されるのに対し、塑性散逸エネルギーｆplastは時間とともに蓄積されるためである。微分方程式（３）では、塑性散逸エネルギーｆplastを微分項としないことにより、塑性散逸エネルギーｆplastの蓄積を表現することができる。
このように、微分方程式（３）には、弾性エネルギーｆelastの解放率を表現する弾性エネルギーｆelastの微分項と、塑性散逸エネルギーｆplastの蓄積を表現する塑性散逸エネルギーｆplastの項と、が含まれる。したがって、微分方程式（３）を計算することにより、脆性破壊及び疲労破壊の双方を考慮して、亀裂を予測することが可能となる。
なお、上述した塑性散逸エネルギーｆplastでは、疲労破壊のうち、低サイクル疲労（塑性疲労）を解析することはできる。しかしながら、高サイクル疲労（弾性疲労）解析について対応することが困難であった。
そこで、本技術の発明者は、さらに、フェーズフィールド法の考え方を応用し、微分方程式（２）に、高サイクルの繰り返し疲労（負荷）を表すエネルギー（以下、振幅負荷エネルギーＡと称する）を含む項を導入することにより、高サイクル疲労（弾性疲労）による亀裂を予測可能となることを見出した。
（フェーズフィールド法の考え方を応用した亀裂予測方法）
本技術に係るフェーズフィールド法の考え方を応用した亀裂予測方法では、微分方程式（２）に、振幅負荷エネルギーＡの項を導入した、次の微分方程式（４）を利用する。
Figure 2019003301
微分方程式（４）は、微分方程式（３）において、塑性散逸エネルギーｆ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に予め振幅負荷エネルギーＡが設定されている場合などには、本ステップＳ１４は省略しても構わない。
（微分方程式生成ステップ）
ステップＳ１５において、微分方程式生成部１５は、ステップＳ１３で設定した亀裂変数φと、ステップＳ１４で設定した振幅負荷エネルギーＡと、が利用されて微分方程式を作成する。
ステップＳ１５において生成される微分方程式としては、一例として、上記の微分方程式（４）が挙げられる。また、ステップＳ１５においては、微分方程式（４）を改良した、次の微分方程式（５）が生成されてもよい。なお、式（５）は、拡散方程式でもある。
Figure 2019003301
微分方程式（５）には、障壁エネルギーｆdoubの微分項のフィッティング用定数ｗdoub、弾性エネルギーｆelastの微分項のフィッティング用定数ｗelast、及び振幅負荷エネルギーＡの項のフィッティング用定数ｗhcが導入されている。これにより、構造体Ｄの構成などに応じて、障壁エネルギーｆdoubの微分項と、弾性エネルギーｆelastの微分項と、振幅負荷エネルギーＡの項と、のそれぞれの重み付けを最適化することが可能である。これにより、構造体Ｄに発生する亀裂を更に正確に予測可能となる。
なお、式（５）における振幅負荷エネルギーＡは、繰り返し負荷に対応した対象材料にかかる応力振幅Δσを、図４のグラフで、応力振幅Δσ＝最大応力σmax−最小応力σminと定義したときに、その負荷に対応するエネルギーである。振幅負荷エネルギーＡは、次の式（６）で表される。
Figure 2019003301
ここで、Δεは、ひずみ振幅であり、最大ひずみεmax−最小ひずみεminである。なお、このひずみは弾性ひずみである。
また、式（５）の振幅負荷エネルギーＡに係るｗhcについて、図５を参照して説明する。
図５は、対象材料のS-N線図を表している。S-N線図は、縦軸に応力振幅Δσ、横軸に破断までの繰り返し数Nをとって疲労試験結果をグラフで表したものである。実際には、繰り返し数Ｎによって、メカニズムの異なる低サイクル疲労（塑性疲労）領域と高サイクル疲労（弾性疲労）領域に分かれており、ｗhcは、高サイクル疲労領域の傾きに対応したパラメータであり、解析上は、亀裂の進行速度を決めるパラメータである。
なお、材料の低サイクル、高サイクルの疲労特性を記述する際には、図６に示される数式である次の式（７）でフィッティングされる場合が多い。
Figure 2019003301
式（７）の右辺の第１項は、低サイクル疲労領域の応力振幅を表し、第２項が、高サイクル疲労領域の応力振幅を表している。その際、ｗhcは、関数ｆを用いて、次の式（８）で表される。
Figure 2019003301
なお、ここで、関数ｆは、線形であってもよいし、非線形であってもよい。
異種材料が、物理結合、化学結合、粘着結合、静電結合などで結合しており、その結合力寿命が各材料のバルク内結合力よりも小さい場合、結合部に繰り返し負荷がかかると界面剥離が進行する。その際の特性は、図５および図６で示されたように、界面での疲労破断についても特性線図が得られる。
解析モデルの界面上に対応したｗ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の微分項が含まれているため、多岐にわたる材料で形成された構造体Ｄに適用可能である。つまり、微分方程式（４），（５）では、高い汎用性が得られる。
この一方で、構造体Ｄを形成する材料によっては、微分方程式（４），（５）には不要な項が含まれていることになる。したがって、ステップＳ１５では、微分方程式について、構造体Ｄを形成する材料に応じて、不要な項を排除するなどのカスタマイズをすることにより、構造体Ｄに発生する亀裂を迅速かつ的確に予測可能となる。
本技術に係る亀裂予測方法は、金属材料や樹脂材料などの高サイクル疲労破壊しやすい材料を用いて構成された構造体Ｄの亀裂を予測可能であればよい。このため、ステップＳ１５で生成する微分方程式は、少なくとも、亀裂変数φの時間微分に比例する項と、ステップＳ１４で設定した振幅負荷エネルギーＡに比例する項と、を含んでいればよい。
また、図１の亀裂予測装置１０は、少なくともモデル取得部１２及び亀裂予測部１６を備えていればよい。つまり、ステップＳ１１，Ｓ１３乃至Ｓ１５を実行しない場合には、亀裂予測装置１０がモデル生成部１１、亀裂変数設定部１３、振幅負荷エネルギー設定部１４、及び微分方程式生成部１５を備えていなくても構わない。また、亀裂予測装置１０は、必要に応じて、上記以外の構成を含んでいてもよい。
以下、構造体Ｄを形成する材料の一例を挙げ、当該材料に応じてカスタマイズされた微分方程式を例示する。なお、構造体Ｄを形成する材料は、以下に挙げるものに限定されず、任意の材料であってよい。また、各材料に対応する微分方程式も、以下に例示するものに限定されず、任意にカスタマイズ可能である。
（ａ）脆性破壊しにくい材料
構造体Ｄを形成する材料が脆性破壊しにくい場合、例えば、振幅負荷エネルギーＡの項以外の項を排除し、高サイクル疲労破壊のみを考慮した、次の微分方程式（９）を用いることができる。
Figure 2019003301
 微分方程式（９）には、亀裂変数φの時間微分の項、及び振幅負荷エネルギーＡの項のみが含まれる。なお、振幅負荷エネルギーＡの項には、フィッティング用定数ｗhcが含まれていてもよい。このように、振幅負荷エネルギーＡの項以外の項を排除して簡単化された微分方程式（９）を用いることにより、計算負荷を大幅に小さくすることができる。
（ｂ）弾性率Ｂが異方性を有する材料
構造体Ｄを形成する材料の弾性率Ｂが異方性を有する場合、例えば、弾性率Ｂの異方性を考慮した、次の微分方程式（１０）を用いることができる。
Figure 2019003301
微分方程式（１０）中、系のエネルギーＦsysは、次の式（１１）で表される。
Figure 2019003301
式（１１）中、勾配エネルギーｆgradは、次の式（１２）で表される。式（１１）中、弾性エネルギーｆelastは、次の式（１３）で表される。
Figure 2019003301
 式（１２）中、κは、材料定数を示す。
Figure 2019003301
 式（１３）中、εは垂直ひずみを示す。
式（１３）では、弾性率Ｂijklをテンソル（行列）として扱うことにより、弾性率Ｂの異方性を適切に予測結果に反映させることができる。このため、微分方程式（１０）を計算することによって、弾性率Ｂが異方性を有する材料で形成された構造体Ｄに発生する亀裂を的確に予測することができる。
なお、構造体Ｄを形成する材料の弾性率Ｂが等方性を有する場合には、弾性エネルギーｆelastは、式（１３）に代えて、次の式（１４）で表すことができる。
Figure 2019003301
 式（１４）中、νはポアソン比を示し、γはせん断ひずみを示す。
また、式（１４）中の弾性率Ｂは、例えば、図１２に示すような亀裂変数φに依存する関数とすることができる。図１２に示す関数では、亀裂変数φが大きくなるにつれて、弾性率Ｂが小さくなる。つまり、図１２に示す関数では、振幅負荷エネルギーＡの蓄積に伴って構造体Ｄを形成する材料の弾性が低下することを表現することができる。
（ｃ）靭性値が異方性を有する材料
構造体Ｄを形成する材料の靭性値が異方性を有する場合、例えば、拡散項の係数を亀裂変数φの勾配の関数、すなわち界面の法線方向の関数とする、次の微分方程式（１５）を用いることができる。
Figure 2019003301
微分方程式（１５）における系のエネルギーＦsysは、次の式（１６）で表される。
Figure 2019003301
式（１６）における勾配エネルギーｆgradは、次の式（１７）で表される。式（１６）における弾性エネルギーｆelastは、次の式（１９）で表される。
Figure 2019003301
 式（１７）におけるκは、材料定数を示し、次の式（１８）で表される。
Figure 2019003301
 式（１８）において、ａは異方性関数である。
Figure 2019003301
 式（１９）において、νはポアソン比を示し、γはせん断ひずみを示す。
式（１８）を用いることによって、界面の向きに応じて拡散係数を変えることができ、亀裂の進行のしやすさを方向によって変化させることができる。これにより、靭性値の異方性を予測結果に適切に反映させることができる。このため、微分方程式（１５）を計算することによって、靭性値が異方性を有する材料で形成された構造体Ｄに発生する亀裂を的確に予測することができる。
（ｄ）脆性破壊と高サイクル疲労破壊とが同時に進行する材料
構造体Ｄが脆性破壊と高サイクル疲労破壊とが同時に進行する材料によって形成されている場合、構造体Ｄには脆性破壊と高サイクル疲労破壊との組み合わせによる亀裂が発生する。脆性破壊と高サイクル疲労破壊との組み合わせによる亀裂を予測するために、例えば、微分方程式（２０）を用いることができる。
Figure 2019003301
微分方程式（２０）における系のエネルギーＦsysは、次の式（２１）で表される。
Figure 2019003301
式（２１）における勾配エネルギーｆgradは、次の式（２２）で表される。式（２１）における弾性エネルギーｆelastは、次の式（２３）で表される。
Figure 2019003301
 式（２２）において、κは材料定数を示す。
Figure 2019003301
 式（２３）において、νはポアソン比を示し、εは垂直ひずみを示し、γはせん断ひずみを示す。
微分方程式（２０）では、弾性エネルギーｆelastの解放率によって脆性破壊を解析することができるとともに、振幅負荷エネルギーＡの蓄積によって高サイクル疲労破壊を解析することもできる。したがって、微分方程式（２０）を計算することによって、構造体Ｄに発生する脆性破壊と高サイクル疲労破壊との組み合わせによる亀裂を予測可能となる。
２．界面の安定化
構造体Ｄにおける亀裂を良好に表現するために、構造体モデルＭDにおいて亀裂を有さない要素Ｅ０と亀裂を有する要素Ｅ１との界面を安定化することが好ましい。つまり、亀裂を有さない要素Ｅ０の亀裂変数φが「０」近傍の値であり、亀裂を有する要素Ｅ１の亀裂変数が「１」近傍の値であり、いずれの要素Ｅもなるべく亀裂変数が「０」と「１」との中間の値である状態とならないことが好ましい。
構造体モデルＭDにおいて亀裂を有さない要素Ｅ０と亀裂を有する要素Ｅ１との界面を安定化するために、例えば、次の微分方程式（２４）を用いることができる。
Figure 2019003301
微分方程式（２４）における系のエネルギーＦsysは、次の式（２５）で表される。
Figure 2019003301
式（２５）中、障壁エネルギーｆdoubは、次の式（２６）で表される。式（２５）中、弾性エネルギーｆelastは、次の式（２７）で表される。
Figure 2019003301
 式（２６）において、κは材料定数を示す。
Figure 2019003301
 式（２７）におけるｅdoubはエネルギー障壁を示す。
上述した式（２７）は、図１３に示すようなダブルウェル関数である。つまり、障壁エネルギーｆdoubが亀裂変数φ＝０，１において極小値をとる。このため、亀裂変数φは、「０」又は「１」近傍の値をとりやすく、「０」と「１」との中間の値をとりにくい。これにより、構造体モデルＭDにおいて亀裂を有さない要素Ｅ０と亀裂を有する要素Ｅ１との界面が安定化する。
また、式（２７）では、構造体Ｄを形成する材料における亀裂の発生しにくさをエネルギー障壁ｅdoubによって表現することが可能である。つまり、亀裂が発生しやすい材料の場合にはエネルギー障壁ｅdoubを低くし、亀裂が発生しにくい材料の場合にはエネルギー障壁ｅdoubを高くすることができる。
以上により、本技術の第１の実施形態によれば、構造体に発生する延性破壊のうち、高サイクル疲労（弾性疲労）による亀裂を迅速に予測することができる。
なお、第１の実施形態は、構造体に発生する疲労破壊のうち、高負荷の塑性ひずみに入らない疲労（低サイクル疲労）を考慮しない場合、低負荷の高サイクル疲労（弾性疲労）しか考慮しない場合について説明した。ただし、材料によっては、構造体に発生する疲労破壊として、低負荷の高サイクル疲労（弾性疲労）だけでなく、高負荷の塑性ひずみ（低サイクル疲労）を考慮が必要な場合があり、この場合について、次の第２の実施形態で説明する。
＜２．第２の実施形態＞
第２の実施形態に係る技術として、低負荷の高サイクル疲労（弾性疲労）および高負荷の低サイクル疲労（塑性疲労）を含む構造体に発生する疲労破壊を予測可能な方法について説明する。なお、以降は、第１の実施形態との差分に焦点をあてて説明する。
本技術の第２の実施形態においては、微分方程式（４）に、塑性散逸エネルギーｆplastの項を導入した、次の微分方程式（２８）を利用する。
Figure 2019003301
微分方程式（２８）は、微分方程式（４）において、塑性散逸エネルギーｆ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の微小増分で積分した量を利用して、例えば、次の式（２９）によって算出可能である。
Figure 2019003301
また、図１６Ｂに斜線で示す領域の面積は、相当応力σと降伏応力σYとの差と、相当塑性ひずみεpと、の積を利用して、例えば、次の式（３０）によって算出可能である。
Figure 2019003301
ただし、式（３０）において、相当応力σが降伏応力σYより小さい場合の塑性散逸エネルギーｆplastをゼロとする。
塑性散逸エネルギーｆplastの表現方法は、構造体Ｄを形成する材料や物理現象などに応じて、正確に亀裂を予測できるように使い分けることが可能である。なお、塑性散逸エネルギーｆplastを表現する関数は、式（２９）及び式（３０）に限定されず、相当応力σと相当塑性ひずみεpとの関係に基づいて適宜作成可能である。
なお、構造体モデルＭDに振幅負荷エネルギーＡおよび塑性散逸エネルギーｆplastが予め設定されている場合などには、本ステップＳ１１４は省略しても構わない。
（微分方程式生成ステップ）
ステップＳ１１５において、微分方程式生成部１１２は、ステップＳ１１３で設定した亀裂変数φと、ステップＳ１１４で設定された振幅負荷エネルギーＡと、ステップＳ１１４で設定された塑性散逸エネルギーｆplastが利用されて微分方程式を作成する。
ステップＳ１１５において生成される微分方程式としては、一例として、上記の微分方程式（２８）が挙げられる。また、ステップＳ１１５においては、微分方程式（２８）を改良した、次の微分方程式（３１）が生成されてもよい。なお、式（３１）は、拡散方程式でもある。
Figure 2019003301
微分方程式（３１）には、障壁エネルギーｆ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 構造体モデル， Ｅ，Ｅ０，Ｅ１ 要素

Claims (11)

  1. 所定の構造体に対応する構造体モデルを取得するモデル取得部と、
    前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
    を備える情報処理装置。
  2. 前記振幅負荷エネルギーは、繰り返し負荷時の応力の最大値と最小値の差分と、歪みの最大値と最小値の差分とが積算されて設定される
    請求項１に記載の情報処理装置。
  3. 前記振幅負荷エネルギーに係る定数は、対応する材料のS-N線図の高サイクル疲労領域の特性線の傾きに応じた値として設定される
    請求項１に記載の情報処理装置。
  4. 前記振幅負荷エネルギーに係る定数は、対応する材料界面のS-N線図の高サイクル疲労剥離領域の特性線の傾きに応じた値として設定される
    請求項１に記載の情報処理装置。
  5. 前記亀裂予測部は、非周期的な繰り返し負荷に対して、応力時間変化の変曲点毎に解析区間を区切り、区切られた解析区間毎に前記構造体に発生する亀裂を予測する
    請求項１に記載の情報処理装置。
  6. 前記亀裂予測部は、前記構造体モデルの各位置に設定され、塑性変形時に散逸されるエネルギーを前記亀裂変数を利用して表現する塑性散逸エネルギーに比例する項をさらに含む前記微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
    請求項１に記載の情報処理装置。
  7. 前記塑性散逸エネルギーは、相当応力を相当塑性ひずみの微小増分で積分した量を利用して設定される
    請求項６に記載の情報処理装置。
  8. 前記塑性散逸エネルギーは、相当応力と降伏応力との差と、相当塑性ひずみとの積を利用して設定され、相当応力が降伏応力より小さい場合にゼロとされる
    請求項６に記載の情報処理装置。
  9. 前記微分方程式は、空間座標の２階微分に比例する拡散項をさらに含む
    請求項１に記載の情報処理装置。
  10. 所定の構造体に対応する構造体モデルを取得し、
    前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する
    情報処理方法。
  11. 所定の構造体に対応する構造体モデルを取得するモデル取得部と、
    前記構造体モデルの各位置に設定され、亀裂の有無を表現する亀裂変数の時間微分に比例する項と、前記構造体モデルの各位置に設定され、繰り返し負荷時にかかるエネルギーを前記亀裂変数を利用して表現する振幅負荷エネルギーに比例する項とを含む微分方程式を計算することにより、前記構造体に発生する亀裂を予測する亀裂予測部と
    して、コンピュータを機能させるプログラム。
JP2017115657A 2017-06-13 2017-06-13 情報処理装置、情報処理方法、及びプログラム Pending JP2019003301A (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2017115657A JP2019003301A (ja) 2017-06-13 2017-06-13 情報処理装置、情報処理方法、及びプログラム
PCT/JP2018/020652 WO2018230339A1 (ja) 2017-06-13 2018-05-30 情報処理装置、情報処理方法、及びプログラム
US16/619,123 US11531801B2 (en) 2017-06-13 2018-05-30 Information processing device, information processing method, and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017115657A JP2019003301A (ja) 2017-06-13 2017-06-13 情報処理装置、情報処理方法、及びプログラム

Publications (1)

Publication Number Publication Date
JP2019003301A true JP2019003301A (ja) 2019-01-10

Family

ID=64659158

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017115657A Pending JP2019003301A (ja) 2017-06-13 2017-06-13 情報処理装置、情報処理方法、及びプログラム

Country Status (3)

Country Link
US (1) US11531801B2 (ja)
JP (1) JP2019003301A (ja)
WO (1) WO2018230339A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109783947A (zh) * 2019-01-21 2019-05-21 南京大学 一种采水型地裂缝数值模拟方法
JP2020177374A (ja) * 2019-04-17 2020-10-29 横浜ゴム株式会社 複合材料の解析方法、及び複合材料の解析用コンピュータプログラム

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111650062B (zh) * 2020-06-30 2023-05-26 武汉钢铁有限公司 金属材料的等幅总应变-寿命曲线的分段式测试方法
CN113849935B (zh) * 2021-09-29 2024-06-04 哈尔滨理工大学 一种相场法中与区域无关积分求解方法
CN114169209B (zh) * 2021-12-24 2023-07-28 哈尔滨工业大学 一种提取非均匀材料混合型弥散裂纹尖端参量的数值方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4378506B2 (ja) * 2004-07-09 2009-12-09 株式会社産学連携機構九州 疲労き裂成長曲線の推定法、推定プログラムおよび推定装置
US7289913B2 (en) * 2004-10-08 2007-10-30 Jentek Sensors, Inc. Local feature characterization using quasistatic electromagnetic sensors
US8209133B2 (en) * 2007-08-16 2012-06-26 Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College Rapid determination of fatigue failure based on temperature evolution
JP2009160028A (ja) 2007-12-28 2009-07-23 Sammy Corp 入賞口ユニット及び遊技機
JP4938695B2 (ja) * 2008-01-23 2012-05-23 富士通株式会社 き裂進展評価装置及びき裂進展評価方法
US9243985B2 (en) * 2009-10-05 2016-01-26 Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College Fracture fatigue entropy determination
JP2011204081A (ja) 2010-03-26 2011-10-13 Fujitsu Ltd 半導体設計プログラム、半導体設計方法および半導体設計装置
US8380447B2 (en) * 2010-04-14 2013-02-19 Simmonds Precision Products, Inc. Quantification of condition indicators in the presence of synchronous noise
US9581570B2 (en) * 2011-02-10 2017-02-28 University Of South Carolina Determination of the remaining life of a structural system based on acoustic emission signals
EP2924596A4 (en) * 2013-05-10 2016-09-14 Nippon Steel & Sumitomo Metal Corp DEFORMATION ANALYSIS DEVICE, DEFORMATION ANALYSIS PROCESS AND PROGRAM
US20150330950A1 (en) * 2014-05-16 2015-11-19 Eric Robert Bechhoefer Structural fatigue crack monitoring system and method
WO2017077668A1 (ja) * 2015-11-05 2017-05-11 ソニー株式会社 情報処理装置、情報処理方法、及びプログラム
EP3229006A1 (en) * 2016-04-08 2017-10-11 Ansaldo Energia Switzerland AG Method for determination of fatigue lifetime consumption of an engine component

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109783947A (zh) * 2019-01-21 2019-05-21 南京大学 一种采水型地裂缝数值模拟方法
JP2020177374A (ja) * 2019-04-17 2020-10-29 横浜ゴム株式会社 複合材料の解析方法、及び複合材料の解析用コンピュータプログラム
JP7310249B2 (ja) 2019-04-17 2023-07-19 横浜ゴム株式会社 複合材料の解析方法、及び複合材料の解析用コンピュータプログラム

Also Published As

Publication number Publication date
US20200117846A1 (en) 2020-04-16
US11531801B2 (en) 2022-12-20
WO2018230339A1 (ja) 2018-12-20

Similar Documents

Publication Publication Date Title
WO2018230339A1 (ja) 情報処理装置、情報処理方法、及びプログラム
Pant et al. Appropriate viscous damping for nonlinear time‐history analysis of base‐isolated reinforced concrete buildings
Jensen et al. Model-reduction techniques for Bayesian finite element model updating using dynamic response data
US6301970B1 (en) Cumulative damage model for structural analysis of filed polymeric materials
Zordan et al. Improved equivalent viscous damping model for base-isolated structures with lead rubber bearings
Alegre et al. Some aspects about the crack growth FEM simulations under mixed-mode loading
JP6408856B2 (ja) 高分子材料のシミュレーション方法
JP6944736B1 (ja) 疲労寿命予測方法、疲労寿命予測装置、疲労寿命予測プログラムおよび記憶媒体
de Larrard et al. Effect of the Young modulus variability on the mechanical behaviour of a nuclear containment vessel
JP2011243197A (ja) 数値的モデルを用いて高分子材料の材料特性をシミュレートする方法およびシステム
CN112580233B (zh) 一种考虑晶粒尺寸的金属强韧性能预测方法
Barbu et al. High cycle fatigue simulation: A new stepwise load-advancing strategy
Flicek et al. Determination of the shakedown limit for large, discrete frictional systems
WO2017077668A1 (ja) 情報処理装置、情報処理方法、及びプログラム
Bhamare et al. A multi-temporal scale approach to high cycle fatigue simulation
Wang Investigating some technical issues on cohesive zone modeling of fracture
Bhattacharyya et al. A model reduction technique in space and time for fatigue simulation
do Rêgo et al. A methodology for fretting fatigue life estimation using strain-based fracture mechanics
JP2005519300A (ja) 異方性材料からなる構成部材の弾塑性挙動を求める方法、およびこの方法の使用法
EP2787457B1 (en) Contact simulation method for rubber material
Gozin et al. Quarter elliptical crack growth using three dimensional finite element method and crack closure technique
JP2019144765A (ja) 柱状ハニカム構造体の最大引張応力推定方法及び最大引張応力推定装置
Gurve et al. Modal and harmonic analysis of stiffened plate using First-order shear deformation theory
JPWO2011108468A1 (ja) 材料定数推定システム及び材料定数推定方法
Kurt et al. Three-dimensional non-planar crack growth analysis using enriched finite elements