(57)【要約】【課題】亀裂が進展する過程を評価すること。【解決手段】評価対象であるタイヤにおける結合界面の界面ポテンシャルエネルギ関数のパラメータを、引張剥離実験又はせん断剥離実験の少なくとも一方に基づいて決定する（ステップＳ１０３）。次に、前記タイヤを有限個の要素に分割して作成した解析モデル中の結合界面に、決定した界面ポテンシャルエネルギを含む界面要素を作成し、配置する（ステップＳ１０４）。そして、界面ポテンシャルエネルギ関数に基づいて前記界面要素の剛性を計算し（ステップＳ１０５）、計算した剛性を解析モデルの全体連立方程式に代入して（ステップＳ１０６）、所定の解析条件に基づいて前記全体連立方程式を解くことにより、解析を実行する（ステップＳ１０７、ステップＳ１０８）。そして、全体連立方程式を解いて得られた前記界面の物理量の大きさから、亀裂の進展の可否を判定し、タイヤの破壊を評価する（ステップＳ１０９）。【選択図】        図２
2008145200.tif 000002
【特許請求の範囲】
【請求項１】
構造体の任意部位における破損の発生及び破損進展現象を予測するにあたり、
評価対象である前記構造体における結合界面の界面ポテンシャルエネルギ関数のパラメータを、引張剥離実験又はせん断剥離実験の少なくとも一方に基づいて決定する手順と、
前記構造体を有限個の要素に分割して作成した解析モデル中の結合界面に、前記パラメータを求めた前記界面ポテンシャルエネルギを含む界面要素を作成し、配置するとともに、前記界面ポテンシャルエネルギ関数に基づいて前記界面要素の剛性を計算する手順と、
前記剛性を前記解析モデルの全体連立方程式に代入し、所定の解析条件に基づいて前記全体連立方程式を解く手順と、
全体連立方程式を解いて得られた、前記結合界面の物理量の大きさから、亀裂の進展の可否を判定する手順と、
を含むことを特徴とする構造体の破壊予測方法。
【請求項２】
前記界面ポテンシャルエネルギ関数のパラメータは、破壊の実験によって得られる荷重点での荷重−変位間の関係と、破壊の解析によって得られる荷重−変位間の関係とが一致するように決定されることを特徴とする請求項１に記載の構造体の破壊予測方法。
【請求項３】
前記構造体の解析モデルに発生した亀裂が進展した後、新しく生成された亀裂面同士が貫通する場合には、新しく生成された前記亀裂面間に相互接触を設定することを特徴とする請求項１又は２に記載の構造体の破壊予測方法。
【請求項４】
前記破壊の実験及び前記破壊の解析は、引張モード、面内せん断モード、又は面外せん断モードのうち少なくとも一つによる破壊の形態で行うことを特徴とする請求項２又は３に記載の構造体の破壊予測方法。
【請求項５】
引張モードが支配的である破壊の形態である場合、引張モードのみで前記界面ポテンシャルエネルギ関数のパラメータを決定し、
面内せん断モードが支配的である破壊の形態である場合、又は面外せん断モードが支配的である破壊の形態である場合、せん断モードのみで前記界面ポテンシャルエネルギ関数のパラメータを決定することを特徴とする請求項４に記載の構造体の破壊予測方法。
【請求項６】
前記構造体において亀裂の進展する方向が予め分かっている場合、前記亀裂が進展可能な方向に沿ってのみ、前記界面要素を配置することを特徴とする請求項１〜５のいずれか１項に記載の構造体の破壊予測方法。
【請求項７】
前記構造体に存在するすべての異種材料間、又は、前記構造体において指定した異種材料間、又は、前記解析モデルに存在するすべての要素の境界、又は、前記解析モデルにおいて指定した領域における前記要素の境界のうち少なくとも一つに、前記界面要素を作成し、配置することを特徴とする請求項１〜５のいずれか１項に記載の構造体の破壊予測方法。
【請求項８】
請求項１〜７のいずれか１項に記載の構造体の破壊予測方法をコンピュータに実行させることを特徴とする構造体の破壊予測用コンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、構造体の破壊の予測に関する。
【背景技術】
【０００２】
コンピュータの発達にともない、近似解析手法を用いたコンピュータシミュレーションによって、材料や構造体の力学的特性を予測・解析する手法が実用化されるようになっている。このようなコンピュータシミュレーションを用いて構造体の破壊現象を予測することにより、構造体の健全性を評価することが行われている。特許文献１には、Ｊ積分値による破壊パラメータを求めるとともに、破壊を阻止する抗力を求め、亀裂進展を把握することで、タイヤの経時変化を予測する技術が開示されている。
【０００３】
【特許文献１】特開２００５−４７２９５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかし、特許文献１に開示されているような、Ｊ積分値の計算によって亀裂の進展を判定する手法では、亀裂の進展の有無は判別できるが、亀裂が進展した後における亀裂の進展や亀裂が徐々に進展する現象の評価ができない。そこで、この発明は、上記に鑑みてなされたものであって、亀裂が進展した後や亀裂が徐々に進展する過程を評価できる構造体の破壊予測方法及び構造体の破壊予測用コンピュータプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００５】
上述の目的を達成するために、本発明に係る構造体の破壊予測方法は、構造体の任意部位における破損の発生及び破損進展現象を予測するにあたり、評価対象である前記構造体における結合界面の界面ポテンシャルエネルギ関数のパラメータを、引張剥離実験又はせん断剥離実験の少なくとも一方に基づいて決定する手順と、前記構造体を有限個の要素に分割して作成した解析モデル中の結合界面に、前記パラメータを求めた前記界面ポテンシャルエネルギを含む界面要素を作成し、配置するとともに、前記界面ポテンシャルエネルギ関数に基づいて前記界面要素の剛性を計算する手順と、前記剛性を前記解析モデルの全体連立方程式に代入し、所定の解析条件に基づいて前記全体連立方程式を解く手順と、全体連立方程式を解いて得られた、前記結合界面の物理量の大きさから、亀裂の進展の可否を判定する手順と、を含むことを特徴とする。
【０００６】
この構造体の破壊予測方法は、亀裂が進展する結合界面間における結合強度を表すものとして、非線形のばねに相当する界面ポテンシャルエネルギ関数を導入する。そして、界面ポテンシャルエネルギ関数を、前記結合界面に生成される界面要素に設定して、前記界面要素が要素間に生成された解析モデルを所定の解析条件に基づいて解析する。このように、非線形のばねに相当する界面ポテンシャルエネルギ関数を界面要素に設定することで、亀裂の進展を結合界面同士の結合強度に基づいて判断できるので、亀裂が進展した後における亀裂の進展や亀裂が徐々に進展する過程を精度よく再現し、評価できる。また、複雑な補強材を有する複合材料であっても、非線形のばねに相当する界面ポテンシャルエネルギ関数を界面要素に設定すればよいので、複雑な補強材を有する複合材料で構成される構造体の破壊を精度よく再現し、評価できる。
【０００７】
次の本発明に係る構造体の破壊予測方法のように、前記構造体の破壊予測方法において、前記界面ポテンシャルエネルギ関数のパラメータは、破壊の実験によって得られる荷重点での荷重−変位間の関係と、破壊の解析によって得られる荷重−変位間の関係とが一致するように決定することが好ましい。
【０００８】
次の本発明に係る構造体の破壊予測方法のように、前記構造体の破壊予測方法において、前記構造体の解析モデルに発生した亀裂が進展した後、新しく生成された亀裂面同士が貫通する場合には、新しく生成された前記亀裂面間に相互接触を設定することが好ましい。
【０００９】
次の本発明に係る構造体の破壊予測方法のように、前記構造体の破壊予測方法において、前記破壊の実験及び前記破壊の解析は、引張モード、面内せん断モード、又は面外せん断モードのうち少なくとも一つによる破壊の形態で行うことが好ましい。
【００１０】
次の本発明に係る構造体の破壊予測方法のように、前記構造体の破壊予測方法において、引張モードが支配的である破壊の形態である場合、引張モードのみで前記界面ポテンシャルエネルギ関数のパラメータを決定し、面内せん断モードが支配的である破壊の形態である場合、又は面外せん断モードが支配的である破壊の形態である場合、せん断モードのみで前記界面ポテンシャルエネルギ関数のパラメータを決定することが好ましい。
【００１１】
次の本発明に係る構造体の破壊予測方法のように、前記構造体の破壊予測方法において、前記構造体において亀裂の進展する方向が予め分かっている場合、前記亀裂が進展可能な方向に沿ってのみ、前記界面要素を配置することが好ましい。
【００１２】
次の本発明に係る構造体の破壊予測方法のように、前記構造体の破壊予測方法において、前記構造体に存在するすべての異種材料間、又は、前記構造体において指定した異種材料間、又は、前記解析モデルに存在するすべての要素の境界、又は、前記解析モデルにおいて指定した領域における前記要素の境界のうち少なくとも一つに、前記界面要素を作成し、配置することが好ましい。
【００１３】
次の本発明に係る構造体の破壊予測用コンピュータプログラムは、前記構造体の破壊予測方法をコンピュータに実行させることを特徴とする。
【００１４】
この構造体の破壊予測用コンピュータプログラムによれば、前述の構造体の破壊予測方法が、コンピュータを利用して実現できる。
【発明の効果】
【００１５】
この発明に係る構造体の破壊予測方法及び構造体の破壊予測用コンピュータプログラムでは、亀裂が進展した後における亀裂の進展や亀裂が徐々に進展する過程を評価できる。
【発明を実施するための最良の形態】
【００１６】
以下、この発明につき図面を参照しつつ詳細に説明する。なお、この発明を実施するための最良の形態（以下実施形態という）によりこの発明が限定されるものではない。また、下記実施形態における構成要素には、当業者が容易に想定できるもの、あるいは実質的に同一のもの、いわゆる均等の範囲のものが含まれる。なお、本発明は、異種材料の境界における剥離や亀裂の進展のみならず、同種材料中に破壊の起点（例えば亀裂等）が発生し、そこから破壊が進展するような場合についても適用できる。また、以下においては、構造体の例としてタイヤを取り上げるが、本発明は、タイヤ以外の構造物に適用してもよい。なお、本発明は、タイヤのような、複雑な補強材を有する複合材料で構成される構造物の任意部位における破損の発生及び破損の進展現象を数値解析手法で予測する際に好適である。
【００１７】
本実施形態は、構造体の任意部位における破損の発生及び破損進展現象を予測するものであり、次の点に特徴がある。すなわち、亀裂が進展する結合界面間における結合強度を表すものとして、非線形のばねに相当する界面ポテンシャルエネルギを導入する。そして、界面ポテンシャルエネルギを、前記結合界面に生成される界面要素に設定して、前記界面要素が要素間に生成された解析モデルを所定の解析条件に基づいて解析する。
【００１８】
図１−１、図１−２は、本実施形態に係る構造体の破壊予測装置を示す説明図である。この構造体の破壊予測装置５０は、処理部５２と記憶部５４とで構成される。また、この構造体の破壊予測装置５０には、入出力装置５１が接続されており、ここに備えられた入力部５３で、材料定数その他の解析条件を処理部５２へ入力する。ここで、入力部５３には、キーボード、マウス等の入力デバイスを使用することができる。また、解析結果は、表示装置５５に出力することができる。
【００１９】
記憶部５４には、本実施形態に係る構造体の破壊予測方法をコンピュータに実行させるためのコンピュータプログラムが格納されている。ここで、記憶部５４は、ハードディスク装置や光磁気ディスク装置、又はフラッシュメモリ等の不揮発性のメモリ、ＣＤ−ＲＯＭ（Read Only Memory）等のような読み出しのみが可能な記憶媒体や、ＲＡＭ（Random Access Memory）のような揮発性のメモリ、あるいはこれらの組み合わせにより構成することができる。
【００２０】
図１−２に示すように、処理部５２は、解析モデル作成部５２Ｍと、界面設定部５２Ｂと、演算部５２Ｃとを含んで構成される。また、処理部５２は、メモリ及びＣＰＵ（Central Processing Unit）により構成されている。本実施形態に係る構造体の破壊予測方法を実行するにあたっては、処理部５２が前記コンピュータプログラムを当該処理部５２に組み込まれたメモリに読み込んで演算する。その際に処理部５２は、適宜記憶部５４へ演算途中の数値を格納し、また格納した数値を取り出して演算を進める。なお、本実施形態に係る構造体の破壊予測方法を実現するためのコンピュータプログラムの代わりに、専用のハードウェアを用いて本実施形態に係る構造体の破壊予測方法を実現してもよい。
【００２１】
なお、上記コンピュータプログラムは、コンピュータシステムへ既に記録されているコンピュータプログラムとの組み合わせによって、本実施形態に係る構造体の破壊予測方法を実現できるものであってもよい。また、上記コンピュータプログラムをコンピュータ読み取り可能な記録媒体に記録するとともに、当該コンピュータプログラムをコンピュータシステムに読み込ませ、実行することにより本実施形態に係る構造体の破壊予測方法を実行してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器などのハードウェアを含むものとする。次に、本実施形態に係る構造体の破壊予測方法の手順について説明する。
【００２２】
図２は、本実施形態に係る構造体の破壊予測方法の手順を示すフローチャートである。図３−１、図３−２は、本実施形態に係る構造体の破壊予測方法で用いるタイヤの解析モデルを作成する手順を示す説明図である。本実施形態に係る構造体の破壊予測方法を実行するにあたって、ステップＳ１０１において、構造体の破壊予測装置５０の処理部５２が備える解析モデル作成部５２Ｍは、図３−１に示す、評価対象であるタイヤ１を複数の微小要素（要素）１１1、１１2、・・・１１ｎに分割して、図３−２に示す解析モデル１０を作成する。ここで、解析モデル１０の前記微小要素１１1、１１2、・・・１１ｎを、便宜上通常要素という。
【００２３】
次に、ステップＳ１０２において、破壊の評価対象となる領域を抽出する。この領域は、評価対象であるタイヤ１において、破壊を評価する領域であり、亀裂の進展する領域、すなわち、破壊の発生する領域である。例えば、カーカスとサイドトレッドゴムとが接する領域や、ベルトとアンダトレッドゴムとが接する領域等のように、異種材料で形成される領域がある。また、キャップトレッドに設けられるリブやブロックが剥離するような場合、キャップトレッド中に破壊を評価する領域が存在する場合のように、同種材料で形成される領域がある。抽出された破壊の評価対象となる領域は、少なくとも座標データが記憶部５４に格納される。
【００２４】
破壊を評価する領域を抽出することにより、その領域における破壊のモードや材料定数等の、亀裂の進展を予測する際に必要な情報を得ることができる。破壊のモードは、例えば、引張のモード、面内せん断のモード、面外せん断のモード、あるいはこれらを組み合わせたモードがある。材料定数は、亀裂が進展する領域、すなわち破壊を評価する領域における材料の弾性係数や破壊靭性等がある。破壊を評価する領域を設定したら、ステップＳ１０３において、破壊を評価する領域における界面ポテンシャルエネルギ関数を決定する。次に、界面ポテンシャルエネルギ関数の設定手順を説明する。
【００２５】
図４は、本実施形態に係る構造体の破壊予測方法において、界面ポテンシャルエネルギ関数を設定する手順を示すフローチャートである。図５−１は、引張の破壊モードの実験に用いる実験サンプルを示す模式図である。図５−２は、引張の破壊モードの実験に用いる解析モデルを示す模式図である。図６−１は、せん断の破壊モードの実験に用いる実験サンプルを示す模式図である。図６−２は、せん断の破壊モードの実験に用いる解析モデルを示す模式図である。図７は、実験結果及び解析結果の一例を示す説明図である。
【００２６】
構造材料の結合界面が剥離したり切断されたりする過程では、新しく破損表面が生成されるため、エネルギの消費が必要であり、この消費したエネルギを表すポテンシャルエネルギ関数が存在する。このポテンシャルエネルギ関数を界面ポテンシャルエネルギ関数といい、新しく破損界面を生成するために消費されるエネルギを、界面ポテンシャルエネルギという。界面ポテンシャルエネルギ関数を決めるため、本実施形態では、界面ポテンシャルエネルギ関数を界面要素に導入した解析モデルによる解析結果と、試験片の破壊実験（引張剥離実験又はせん断剥離実験）によって得られた実験結果とが一致するように、界面ポテンシャルエネルギ関数のパラメータを定める。
【００２７】
界面ポテンシャルエネルギ関数φは、構造体に作用する荷重をＦ、この荷重による構造体の変形をｕとすると、∂φ／∂ｕ＝［Ｆ］、∂2φ／∂ｕ2＝［Ｋｓ］の関係がある。ここで、Ｋｓは、剥離やせん断等によって破壊する部分における結合界面の剛性を表し、後述する界面要素の剛性である。なお、［Ｆ］、［Ｋｓ］は、いずれも行列式の形で記述されている。
【００２８】
界面ポテンシャルエネルギ関数を設定する際には、ステップＳ２０１において、上記ステップＳ１０２で抽出された破壊の評価対象領域における界面ポテンシャルエネルギ関数φを定義する。次にステップＳ２０２−１において、引張あるいはせん断の破壊実験に用いる試験片を準備する。引張破壊実験に用いる試験片は、図５−１に示す引張用試験片２であり、せん断破壊試験に用いる試験片は、図５−２に示すせん断用試験片３である。
【００２９】
また、ステップＳ２０２−２において、解析モデル作成部５２Ｍは、パラメータ設定用解析モデルを作成する。パラメータ設定用解析モデルは、例えば、有限要素法に基づいて作成される。図５−１に示す引張用試験片２からパラメータ設定用解析モデルを作成する場合、解析モデル作成部５２Ｍは、引張用試験片２を複数の微小要素（要素）１２1、１２2、・・・１２ｎに分割して、図５−２に示す引張破壊パラメータ設定用解析モデル１２を作成する。また、図６−１に示すせん断用試験片３からパラメータ設定用解析モデルを作成する場合、解析モデル作成部５２Ｍは、引張用試験片２を複数の微小要素（要素）１３1、１３2、・・・１３ｎに分割して、図６−２に示すせん断破壊パラメータ設定用解析モデル１３を作成する。なお、ステップＳ２０２−１とステップＳ２０２−２の順序は問わない。
【００３０】
ここで、引張破壊及びせん断破壊について説明する。引張破壊（あるいは引張剥離破壊）は、図５−１に示すような破壊のモードであり、引張用試験片２の一部を固定部２Ｆで固定するとともに、引張用試験片２の一端２Ｔを固定部２Ｆから引き離すように荷重Ｆを与えることによって、亀裂が矢印Ｙの方向に進展して、引張用試験片２の一部が剥離する破壊のモードである。せん断破壊は、図６−１に示すような破壊のモードであり、せん断用試験片３の一部を固定部３Ｆで固定するとともに、せん断用試験片３を切断するようにせん断用試験片３の一端３Ｔに荷重Ｆを与えることによって、せん断用試験片３の一部がせん断（あるいはせん断剥離）する破壊のモードである。
【００３１】
基本的な破壊形式には、引張（開口）破壊のモード、面内せん断破壊のモード及び面外せん断破壊のモードがある。通常の破壊形式は、これら３種類の基本的な破壊のモードに分けて考えることができる。上述した３種類の独立な破壊のモードの少なくとも２つのモードを組み合わせると、複雑な破壊形式を表すことができる。
【００３２】
例えば、破壊が主に引張（開口）破壊のモードに支配される場合、引張破壊（あるいは引張剥離）の実験のみに基づいて界面ポテンシャルエネルギ関数のパラメータを決定することができる。また、破壊が主にせん断（面内せん断、あるいは面外せん断）モードに支配される場合、せん断破壊（あるいはせん断剥離）の実験のみに基づいて界面ポテンシャルエネルギ関数のパラメータを決定することができる。
【００３３】
本実施形態においては、上述した３種類のモードのうち少なくとも一つに基づいて決定する。これによって、界面ポテンシャルエネルギ関数を容易に決定でき、また、均質材のみならず、カーカスやベルトのような繊維補強材料を含んでいる複合材料であっても、界面ポテンシャルエネルギ関数を容易に決定することができる。
【００３４】
図５−２に示す引張破壊パラメータ設定用解析モデル１２は、引張破壊を再現するため、引張用試験片２の結合界面で亀裂が発生、進展して、引張破壊が発生するものとして、引張用試験片２を解析モデル化する。引張破壊パラメータ設定用解析モデル１２は、引張用試験片２の結合界面に界面要素２０を配置して、この界面要素２０に、ステップＳ２０１で定義された界面ポテンシャルエネルギ関数φを設定する。そして、引張破壊パラメータ設定用解析モデル１２の一端１２Ｔに荷重Ｆが与えられて、引張破壊が再現される。
【００３５】
図６−２に示すせん断破壊パラメータ設定用解析モデル１３は、せん断破壊を再現するため、せん断用試験片３の結合界面で亀裂が発生、進展して、せん断破壊が発生するものとして、せん断用試験片３を解析モデル化する。せん断破壊パラメータ設定用解析モデル１３は、せん断用試験片３の結合界面に界面要素２１を配置して、この界面要素２１に、ステップＳ２０１で定義された界面ポテンシャルエネルギ関数φを設定する。そして、せん断破壊パラメータ設定用解析モデル１３の一端１３Ｔに荷重Ｆが与えられて、せん断破壊が再現される。
【００３６】
実際のタイヤの使用状況を考慮するため、引張用試験片２やせん断用試験片３は、新品又は旧品のタイヤからカットして作製することができる。また、評価対象の構造物であるタイヤと同等の材料から作製してもよい。引張破壊パラメータ設定用解析モデル１２やせん断破壊パラメータ設定用解析モデル１３は、引張用試験片２やせん断用試験片３と同じ寸法に設定することが好ましいが、実験の便利性を考慮した上で、任意の寸法を使用してもよい。
【００３７】
引張あるいはせん断の破壊実験に用いる試験片を準備し、パラメータ設定用解析モデルを作成したら、ステップＳ２０３において、パラメータ設定用解析モデルの界面要素に設定した界面ポテンシャルエネルギ関数φのパラメータを設定する。次に、ステップＳ２０４において、構造体の破壊予測装置５０の処理部５２が備える演算部５２Ｃは、設定したパラメータや境界条件（荷重や固定等）で、パラメータ設定用解析モデルを解析する。
【００３８】
引張破壊パラメータ設定用解析モデル１２やせん断破壊パラメータ設定用解析モデル１３に配置する界面要素２０は、亀裂が発生し、進展を開始した後における荷重の状態によっては、新しく生成された亀裂面間で貫通が発生するおそれがある。これを防ぐため、解析を開始する前に、亀裂が進展可能な部分に配置した界面要素の節点同士に、相互接触を設定しておくことが好ましい。
【００３９】
パラメータ設定用解析モデルの解析が終了したら、ステップＳ２０５において、構造体の破壊予測装置５０の処理部５２が備える界面設定部５２Ｂは、実験結果と解析結果との相違が所定の範囲内にあるか否かを判定する。なお、ステップＳ２０５における判定の前に、引張あるいはせん断の破壊実験は予め済ませておき、結果を記憶部５４に格納しておく。ここで、引張あるいはせん断の破壊実験、及びパラメータ設定用解析モデルの解析について説明する。
【００４０】
図７に示すように、本実施形態では、引張用試験片２やせん断用試験片３（以下必要に応じて試験片という）に引張あるいはせん断の荷重Ｆを与え、荷重Ｆと変位ｕとの間の関係を求める。図７の実線ＥＸが、引張あるいはせん断の破壊実験によって得られた荷重−変位の関係である。一方、引張破壊パラメータ設定用解析モデル１２やせん断破壊パラメータ設定用解析モデル１３（以下必要に応じてパラメータ設定用解析モデルという）も、引張あるいはせん断の破壊実験と同じ条件で解析されて、荷重Ｆと変位ｕとの間の関係が求められる。図７の点線ＡＮ１、一点鎖線ＡＮ２、二点鎖線ＡＮ３が、パラメータ設定用解析モデルを解析することによって得られた荷重−変位の関係である。点線ＡＮ１、一点鎖線ＡＮ２、二点鎖線ＡＮ３は、それぞれ界面ポテンシャルエネルギ関数φのパラメータが異なる。
【００４１】
このように、界面ポテンシャルエネルギ関数φのパラメータを異ならせることによって、パラメータ設定用解析モデルの解析結果が異なる。本実施形態では、実験結果との相違が所定の範囲内にある解析結果が得られた界面ポテンシャルエネルギ関数φのパラメータを、界面ポテンシャルエネルギ関数φのパラメータとして用いる。図７に示す例では、例えば、実験結果と解析結果との相違が、荷重Ｆ及び変位ｕの全範囲において所定の閾値ΔＤの範囲内にある場合に、実験結果と解析結果との相違が所定の範囲内にあると判定し、そのときの界面ポテンシャルエネルギ関数φのパラメータを、構造体の破壊予測に用いる界面ポテンシャルエネルギ関数φのパラメータとする。図７に示す例では、二点鎖線ＡＮ３で示す解析結果が実験結果に対して所定の閾値ΔＤの範囲内にあるので、二点鎖線ＡＮ３で示す解析結果が得られたときにおける界面ポテンシャルエネルギ関数φのパラメータを、構造体の破壊予測に用いる界面ポテンシャルエネルギ関数φのパラメータとする。
【００４２】
ステップＳ２０５においてＮｏと判定された場合、すなわち、界面設定部５２Ｂが、実験結果と解析結果との相違が所定の範囲内にないと判定した場合、ステップＳ２０６において界面ポテンシャルエネルギ関数φのパラメータを変更する。そして、ステップＳ２０４に戻って、演算部５２Ｃが、設定したパラメータや境界条件で、パラメータ設定用解析モデルを解析する。なお、変更するのは、界面ポテンシャルエネルギ関数φのパラメータのみであり、他の解析条件は変更しない。
【００４３】
ステップＳ２０５においてＮｏと判定された場合、すなわち、界面設定部５２Ｂが、実験結果と解析結果との相違が所定の範囲内にあると判定した場合、ステップＳ２０７において、そのときの界面ポテンシャルエネルギ関数φのパラメータが、構造体の破壊予測に用いる界面ポテンシャルエネルギ関数φのパラメータとして決定される。そして、ステップＳ２０８において、ステップＳ２０７で決定されたパラメータを用いて、構造体（本実施形態ではタイヤ）の破壊予測に用いる界面ポテンシャルエネルギ関数φが決定される。
【００４４】
ステップＳ２０８において、タイヤの破壊予測に用いる界面ポテンシャルエネルギ関数φが決定されたら、図２に示すステップＳ１０４に進む。ステップＳ１０４では、ステップＳ１０１で作成した解析モデル１０に、界面要素が配置される。界面要素は、解析モデル１０において亀裂の発生及び進展する部分であり、破壊の起点となる部分である。界面要素は、異種材料の境界や、同種材料中に存在する破壊の起点（例えば亀裂等）から破壊が進展するような領域に配置される。
【００４５】
図８は、界面要素を配置する手順の一例を示すフローチャートである。図９−１〜図９−３は、界面要素を配置する手順の一例を示す説明図である。図９−１には、解析対象の構造体であるタイヤ１の一部の断面構成が示してある。タイヤ１は、第１部材Ｍ（例えば、ゴム）と、第２部材Ｗ（例えば、ワイヤ等の繊維材料）とが、境界Ｓで接合されている。このタイヤ１では、破壊の評価対象領域が、境界Ｓに設定されている。均質材料の場合、破壊の評価対象領域は任意の箇所に設定され、前記破壊の評価対象領域に、界面要素が設定される。
【００４６】
ここで、本実施形態においては、解析モデル１０中において、亀裂が発生、進展し、破壊の基点となる領域に界面要素が設定される。すなわち、解析モデル１０において界面要素の部分に亀裂が発生し、進展するように取り扱われ、タイヤ１における破壊の評価対象領域に荷重が作用していない場合、タイヤ１の界面領域は体積を持たないことになる。このため、本実施形態において、解析モデル１０に配置される界面要素は、荷重を受けないときには体積を持たないように設定される。
【００４７】
界面要素を配置するにあたり、ステップＳ３０１において、界面設定部５２Ｂは、図９−２に示すように、解析モデル１０の境界Ｓの部分に境界要素２０ａ、２ｂ・・・２ｄ等を生成する。界面要素２０ａ、２０ｂ・・・２０ｄ等は、境界Ｓ上の共通節点５ａ、５ｂ・・・５ｅを利用して作成する。なお、共通節点５ａ、５ｂ・・・５ｅは、通常要素１１1、１１2、・・・１１ｎを構成する節点でもある。ここで、境界Ｓ上に共通節点５ａ、５ｂ・・・５ｅがない場合、境界Ｓ上に共通節点が生成されるように、解析モデル１０の通常要素を新たに生成しなおしてもよい。
【００４８】
図９−３に示すように、本実施形態において、界面要素２０ａ、２０ｂ・・・２０ｄは、共通節点５ａ、５ｂ・・・５ｅと、共通節点５ａ、５ｂ・・・５ｅと同じ位置（すなわちδ＝０の位置）に、新たに生成した新節点６ａ、６ｂ・・・６ｅとで構成される。なお、δは０でもよく、この場合には、共通節点５ａ、５ｂ・・・５ｅと新節点６ａ、６ｂ・・・６ｅとの座標情報は同じものとなる。なお、界面要素２０ａ、２０ｂ・・・２０ｄを生成しても、第１部材モデルＭ＿ｍと第２部材モデルＷ＿ｍとの厚さＬは、図９−２に示す解析モデル１０Ａにおける第１部材モデルＭ＿ｍと第２部材モデルＷ＿ｍとの厚さＬに等しい。
【００４９】
境界Ｓ内に界面要素２０ａ、・・２０ｄが生成されたら（ステップＳ３０１、図９−３）、ステップＳ３０２において、界面設定部５２Ｂは、生成した界面要素２０ａ、２０ｂ・・・２０ｄに隣接する通常要素（１１Ａ、１１Ｂ、１１ａ、１１ｂ等）の結合情報（座標や力の関係等）を修正する。これによって、図２のステップＳ１０２で抽出した、解析モデル１０中における破壊の評価対象領域に、界面要素２０ａ、２０ｂ・・・２０ｄ（以下、界面要素２０という）が配置される。なお、界面要素を配置する手法はこれに限られず、例えば、図１−１に示す構造体の破壊予測装置５０が備える表示装置５５に解析モデル１０を表示して、入力部５３を用いたマニュアル操作によって配置してもよい。
【００５０】
図１０−１〜図１３−２は、界面要素の種類を示す説明図である。本実施形態において、界面要素２０は、２つの線、あるいは２つの面から構成され、結合界面的力学特性を表すように、また、荷重を受けないときには体積を持たないように設定する。図１０−１に示す界面要素２０Ａは、線要素の界面要素であり、２個の節点６から構成される。図１１−１に示す界面要素２０Ｂは、２次元の界面要素であり、４個の節点６から構成される。図１２−１に示す界面要素２０Ｃは、８個の節点６から構成される３次元の界面要素であり、図１３−１に示す界面要素２０Ｄは、６個の節点６から構成される３次元の界面要素である。界面要素２０Ａ〜２０Ｄは、荷重Ｆが負荷されていない場合には体積が０である。しかし、図３−２に示す解析モデル１０の解析が開始され、荷重Ｆが負荷されると、図１０−２、図１１−２、図１２−２、図１３−２に示すように、その荷重Ｆに応じて変位が発生し、それに応じて体積が増加する。
【００５１】
図１４−１、図１４−２は、解析モデル中に界面要素を配置した例を示す模式図である。図１４−１、図１４−２は、タイヤを解析モデル化したものであり、図１４−１には、タイヤのサイド部分における亀裂の発生、及び任意方向への亀裂の進展を評価するために、解析モデル１０中に界面要素２０を配置した例を示す。図１４−２には、タイヤのベルトエッジ部分における亀裂の発生、及び任意方向への亀裂の進展を評価するために、解析モデル１０のベルト１５の周辺に界面要素２０を配置した例を示す。
【００５２】
図１４−１、図１４−２に示す解析モデル１０では、解析モデル１０に含まれる要素（通常要素）の境界のうち、指定した領域に存在する要素の境界に界面要素２０を配置した例を示す。このように、解析モデル１０に含まれる要素のうち、指定した領域に界面要素２０を配置すれば、界面要素２０を配置する手間を低減でき、また、計算の負荷を低減させることによって計算効率が向上する。なお、解析モデル１０に含まれるすべての要素の境界に界面要素２０を配置してもよく、このようにすれば、解析精度を向上させることができる。また、界面要素２０は、構造体であるタイヤに存在するすべての異種材料間に配置してもよく、タイヤにおいて指定した異種材料間に配置してもよい。
【００５３】
さらに、構造体であるタイヤにおいて、発生した亀裂が進展する方向が予め分かっている場合、亀裂が進展可能な方向に沿ってのみ、界面要素２０を配置してもよい。このようにすれば、亀裂の進展しない方向に界面要素２０を配置する必要はないので、界面要素２０の配置の負荷を低減できる。また、計算の負荷を低減できるので、計算効率が向上する。
【００５４】
上記手法によって解析モデル１０中における破壊の評価対象領域に界面要素を配置したら（ステップＳ１０４：図２）、界面設定部５２Ｂは、図２に示すステップＳ１０５で、界面要素２０の剛性［Ｋｓ］を計算する。ここで、ステップＳ１０４とステップＳ１０５との順序は問わず、界面要素２０の剛性［Ｋｓ］を計算してから、界面要素２０を解析モデル１０内における破壊の評価対象領域へ配置してもよい。
【００５５】
界面要素２０の剛性［Ｋ］を計算したら、ステップＳ１０６において、解析モデル作成部５２Ｍは、図３−２に示す解析モデル１０を表す全体連立方程式（解析モデル全体連立方程式という）に界面要素２０の剛性［Ｋｓ］を代入する。ここで、解析モデル全体連立方程式は、式（１）で表され、また、界面要素の剛性［Ｋｓ］を代入した解析モデル全体連立方程式（修正全体連立方程式という）は、式（２）で表される。ここで、［Ｋ］は解析モデル１０の通常要素の剛性であり、［Ｆ］は解析モデル１０に与えられる荷重であり、［ｕ］は前記荷重によって解析モデル１０に発生する変位である。［Ｋ］、［Ｋｓ］、［Ｆ］、［ｕ］は、いずれも行列式である。
［Ｋ］＝［Ｆ］［ｕ］・・・（１）
｛［Ｋ］＋［Ｋｓ］｝＝［Ｆ］［ｕ］・・・（２）
【００５６】
ステップＳ１０６で修正全体連立方程式を作成したら、ステップＳ１０７で荷重Ｆの大きさや方向等の解析条件が設定される。そして、ステップＳ１０８で、構造体の破壊予測装置５０の処理部５２が備える演算部５２Ｃは、前記解析条件で修正全体連立方程式を解くことにより、解析モデル１０が破壊する過程を解析する。
【００５７】
図１５−１〜図１５−３は、解析中における亀裂の進展を示す説明図である。解析モデル１０の解析中、発生した亀裂が進展すると、図１５−１に示すように、亀裂面２５Ａ、２５Ｂが新たに生成される。このとき、荷重の作用する方向によっては、図１５−２に示すように、新たに生成された亀裂面２５Ａ、２５Ｂ同士が貫通して、正しい解析が行えなくなるおそれがある。このため、本実施形態では、新たに生成された亀裂面２５Ａ、２５Ｂ同士が貫通することを回避するため、解析モデル作成部５２Ｍは、解析モデル１０に発生した亀裂が進展した後、新しく生成された亀裂面２５Ａ、２５Ｂ同士が貫通する場合には、新しく生成された前記亀裂面２５Ａ、２５Ｂ間に相互接触を設定する。これによって、図１５−３に示すように、新しく生成された亀裂面２５Ａ、２５Ｂ同士の貫通を回避できる。
【００５８】
解析が終了したら、ステップＳ１０９で破壊の評価をする。これは、解析モデル１０に配置した界面要素２０を構成する節点間における結合力、応力、ひずみ、反力等の、界面要素２０を構成する節点間の結合の強さを表す物理量に基づいて評価する。例えば、界面要素２０を構成する節点間における結合力の大きさによって、界面要素２０に亀裂が発生したか否かを判定したり、前記結合力の大きさによって、亀裂の進展を評価したりする。ステップＳ１１０では、ステップＳ１０９の破壊の評価結果に基づき、例えば、タイヤ１の耐久性を予測する。
【００５９】
以上、本実施形態では、構造体の任意部位における破損の発生及び破損進展現象を予測するにあたり、亀裂が進展する結合界面間における結合強度を表すものとして、非線形のばねに関数を、前記結合界面に生成される界面要素に設定して、前記界面要素が要素間に生成された解析モデルを所定の解析条件に基づいて解析する。このように、非線形のばねに相当する界面ポテンシャルエネルギ関数を界面要素に設定することで、亀裂の進展を結合界面同士の結合強度に基づいて判断できるので、亀裂が進展した後における亀裂の進展や亀裂が徐々に進展する過程を評価できる。また、複雑な補強材を有する複合材料であっても、非線形のばねに相当する界面ポテンシャルエネルギ関数を界面要素に設定すればよいので、複雑な補強材を有する複合材料で構成される構造体の破壊を簡易に評価できる。
【００６０】
また、本提案の方法により、界面ポテンシャルエネルギ関数を界面要素に導入した解析モデルによる解析結果と、剥離実験又はせん断実験の実験結果に一致するように、基づいて界面ポテンシャルエネルギ関数のパラメータを決定する。これによって、界面結合強度の大きさによって亀裂の発生や進展を判定できるので、亀裂の駆動力である破壊物理量パラメータを計算し、それらと限界値との比較に基づいて亀裂進展判定したり、解析途中で解析モデルを修正したりする必要はない。その結果、簡易な手法で、亀裂の発生や進展を評価することができる。
【産業上の利用可能性】
【００６１】
以上のように、本発明に係る構造体の破壊予測方法及び構造体の破壊予測用コンピュータプログラムは、構造体の破壊のシミュレーションに有用であり、特に、亀裂が進展した後における亀裂の進展や亀裂が徐々に進展する過程を評価することに適している。
【図面の簡単な説明】
【００６２】
【図１−１】本実施形態に係る構造体の破壊予測装置を示す説明図である。
【図１−２】本実施形態に係る構造体の破壊予測装置を示す説明図である。
【図２】本実施形態に係る構造体の破壊予測方法の手順を示すフローチャートである。
【図３−１】本実施形態に係る構造体の破壊予測方法で用いるタイヤの解析モデルを作成する手順を示す説明図である。
【図３−２】本実施形態に係る構造体の破壊予測方法で用いるタイヤの解析モデルを作成する手順を示す説明図である。
【図４】本実施形態に係る構造体の破壊予測方法において、界面ポテンシャルエネルギ関数を設定する手順を示すフローチャートである。
【図５−１】引張の破壊モードの実験に用いる実験サンプルを示す模式図である。
【図５−２】引張の破壊モードの実験に用いる解析モデルを示す模式図である。
【図６−１】せん断の破壊モードの実験に用いる実験サンプルを示す模式図である。
【図６−２】せん断の破壊モードの実験に用いる解析モデルを示す模式図である。
【図７】実験結果及び解析結果の一例を示す説明図である。
【図８】界面要素を配置する手順の一例を示すフローチャートである。
【図９−１】界面要素を配置する手順の一例を示す説明図である。
【図９−２】界面要素を配置する手順の一例を示す説明図である。
【図９−３】界面要素を配置する手順の一例を示す説明図である。
【図１０−１】界面要素の種類を示す説明図である。
【図１０−２】界面要素の種類を示す説明図である。
【図１１−１】界面要素の種類を示す説明図である。
【図１１−２】界面要素の種類を示す説明図である。
【図１２−１】界面要素の種類を示す説明図である。
【図１２−２】界面要素の種類を示す説明図である。
【図１３−１】界面要素の種類を示す説明図である。
【図１３−２】界面要素の種類を示す説明図である。
【図１４−１】解析モデル中に界面要素を配置した例を示す模式図である。
【図１４−２】解析モデル中に界面要素を配置した例を示す模式図である。
【図１５−１】解析中における亀裂の進展を示す説明図である。
【図１５−２】解析中における亀裂の進展を示す説明図である。
【図１５−３】解析中における亀裂の進展を示す説明図である。
【符号の説明】
【００６３】
１  タイヤ
２  引張用試験片
３  せん断用試験片
１０  解析モデル
１２  引張破壊パラメータ設定用解析モデル
１３  せん断破壊パラメータ設定用解析モデル
２０、２０ａ、２０ｂ、２０ｃ、２０ｄ  界面要素
５０  構造物の破壊予測装置
５１  入出力装置
５２  処理部
５２Ｂ  界面設定部
５２Ｃ  演算部
５２Ｍ  解析モデル作成部
５３  入力部
５４  記憶部
５５  表示装置
【図１−１】
2008145200.tif 000003
【図１−２】
2008145200.tif 000004
【図２】
2008145200.tif 000005
【図３−１】
2008145200.tif 000006
【図３−２】
2008145200.tif 000007
【図４】
2008145200.tif 000008
【図５−１】
2008145200.tif 000009
【図５−２】
2008145200.tif 000010
【図６−１】
2008145200.tif 000011
【図６−２】
2008145200.tif 000012
【図７】
2008145200.tif 000013
【図８】
2008145200.tif 000014
【図９−１】
2008145200.tif 000015
【図９−２】
2008145200.tif 000016
【図９−３】
2008145200.tif 000017
【図１０−１】
2008145200.tif 000018
【図１０−２】
2008145200.tif 000019
【図１１−１】
2008145200.tif 000020
【図１１−２】
2008145200.tif 000021
【図１２−１】
2008145200.tif 000022
【図１２−２】
2008145200.tif 000023
【図１３−１】
2008145200.tif 000024
【図１３−２】
2008145200.tif 000025
【図１４−１】
2008145200.tif 000026
【図１４−２】
2008145200.tif 000027
【図１５−１】
2008145200.tif 000028
【図１５−２】
2008145200.tif 000029
【図１５−３】
2008145200.tif 000030