(57)【特許請求の範囲】
【請求項１】
ゴムと、シリカと、これらを結合する界面結合剤とを含むゴム材料を数値解析が可能な要素でモデル化したゴム材料モデルを設定するステップと、
前記ゴム材料モデルに条件を設定して変形計算を行うステップと、
前記変形計算から必要な物理量を取得するステップとを含むゴム材料のシミュレーション方法であって、
前記ゴム材料モデルは、ゴムマトリックスをモデル化したマトリックスモデルと、
前記マトリックスモデル中に配置されかつ前記シリカをモデル化した複数個のシリカモデルと、
前記各シリカモデルの周りを環状に取り囲みかつ前記マトリックスモデルよりも硬い物性が定義された界面モデルとを含み、
前記シリカモデルは、前記界面モデルを介して紐状に連結された複数個を含み、
前記ゴム材料モデルを設定するステップは、前記界面モデルに物性のパラメータを決定するパラメータ決定工程を含み、
前記パラメータ決定工程は、
解析対象となる前記ゴム材料モデルと同じ配合を有する第１の未加硫ゴム組成物と、該第１の未加硫ゴム組成物から界面結合剤を除いた配合の第２の未加硫ゴム組成物とを少なくとも準備する工程と、
前記第１及び第２の未加硫ゴム組成物をそれぞれ溶剤に浸漬することにより、該第１及び第２の未加硫ゴム組成物からそれぞれマトリックスゴムを除去した第１及び第２の残留物を得る工程と、
少なくとも前記第２の残留物のｔａｎδのピーク温度Ｔ２と、前記第１の残留物のｔａｎδのピーク温度Ｔ１とを測定し、これらに基づいて前記ピーク温度の差Ｔ２−Ｔ１を推定する工程と、
前記第１の未加硫ゴム組成物からシリカを除いた配合の基本の加硫ゴム材料を少なくとも含むとともに該基本の加硫ゴム材料とは架橋密度のみが異なる複数種類の加硫ゴム材料を準備する工程と、
前記各加硫ゴム材料のｔａｎδのピーク温度を測定し、前記加硫ゴム材料のｔａｎδのピーク温度と架橋密度との関係を得る工程と、
前記関係から、ｔａｎδのピーク温度が前記基本の加硫ゴム材料のｔａｎδのピーク温度Ｔ３に前記差Ｔ２−Ｔ１を加えた温度に等しい加硫ゴム材料の架橋密度を特定する工程と、
前記特定された架橋密度の加硫ゴム材料の物性に基づいて前記界面モデルの前記パラメータを定義する工程とを含むことを特徴とするゴム材料のシミュレーション方法。
【請求項２】
前記紐状に連結されたシリカモデルは、任意の軸線方向に関して、ゴム材料モデルの領域の全範囲を連続してのびる請求項１記載のゴム材料のシミュレーション方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ゴムと、シリカと、これらを結合する界面結合剤とが配合されたゴム材料の変形を精度良く解析するのに役立つゴム材料のシミュレーション方法に関する。
【背景技術】
【０００２】
従来、タイヤ、スポーツ用品、その他各種の工業製品に使用されているゴム材料には、その機械的特性を向上させるために、カーボンが多用されていた。しかしながら、近年では、カーボンに代えてシリカが多用されつつある。その理由は、シリカ配合はカーボン配合に比べてエネルギーロスが小さいので、例えばタイヤの転がり抵抗を小さくし、燃費性能の向上に寄与するためである。また、シリカは脱石油資源であるため、環境にも優しいフィラーとも言える。従って、シリカが配合されたゴム材料の変形を精度良くコンピュータシミュレーションで解析することは今後のタイヤ等の開発にきわめて有益となる。
【０００３】
従来、コンピュータを用いたゴム材料のシミュレーション方法としては、下記の非特許文献１及び特許文献１などが知られている。特許文献１のものでは、ゴム材料モデルとして、数値解析が可能な要素でモデル化したマトリックスモデルとフィラーモデルとが設定され、フィラーの影響を考慮した有限要素法による変形計算が行われる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許３６６８２３８号公報
【非特許文献１】Ellen M. Arruda and Marry C. Boyce著「 A THREE-DIMENSIONAL CONSTITUTIVE MODEL FOR THE LARGE STRECH BEHAVIOR OF RUBBER ELASTIC MATERIALSS」 Journal of the Mechanics and Physics of Solids Volume 41, Issue 2, Pages 389-412 (February 1993)
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、シリカ配合ゴム材料のシミュレーションを行う際には、図１５に示されるように、ゴムマトリックスを模したマトリックスモデル２１と、その中に分散配置されかつシリカをモデル化したシリカモデル２２と、該シリカモデル２２の表面を環状に取り囲む界面モデル２３とからなるゴム材料モデル２０が採用されていた。
【０００６】
前記マトリックスモデル２１及びシリカモデル２２には、それぞれの物性が予め定義される。例えば、シリカモデル２２は、硬質の弾性体として取り扱われるので、弾性率が設定される。また、マトリックスモデル２１には、前記物性として、フィラーが配合されていないいわゆる純ゴムの引張試験の結果に基づいた物性（応力と伸びとの関係を示す関数）が定められる。
【０００７】
また、シリカ配合ゴム材料では、シリカとゴムとを結合するために、シランカップリング剤等の界面結合剤が配合される。従って、これまでのシミュレーションにおいては、シリカとゴムとが界面結合剤によって結合されていることを前提として、シリカとゴムとの界面に前記界面モデル２３が定義されていた。
【０００８】
しかしながら、このようなゴム材料モデル２０を使用して、例えば引張変形を与えたときのシミュレーションを行った場合、応力−伸びの結果が、実際の引張試験で得られた応力−伸びの結果と大きく異なる場合がある。種々の実験の結果、発明者らは、このような結果のずれは、シリカモデルの配置に原因があることを知見した。
【０００９】
即ち、発明者らは、シリカの配置についてさらに鋭意研究を重ねた結果、界面結合材は、シリカとマトリックスゴムとを接着させるのみならず、シリカ同士を紐状に連結してシリカ粒子の動きを拘束しているとの仮定を打ち立てた。そして、このようなモデルに従ってシミュレーションを行った結果、実際の試験結果と非常に高い相関性を得ることができることを見出し本発明を完成させるに至った。
【００１０】
以上のように、本発明は、精度良くシリカ配合ゴム材の変形を計算しうるシミュレーション方法を提供することを主たる目的としている。
【課題を解決するための手段】
【００１１】
本発明のうち請求項１記載の発明は、ゴムと、シリカと、これらを結合する界面結合剤とを含むゴム材料を数値解析が可能な要素でモデル化したゴム材料モデルを設定するステップと、前記ゴム材料モデルに条件を設定して変形計算を行うステップと、前記変形計算から必要な物理量を取得するステップとを含むゴム材料のシミュレーション方法であって、前記ゴム材料モデルは、ゴムマトリックスをモデル化したマトリックスモデルと、前記マトリックスモデル中に配置されかつ前記シリカをモデル化した複数個のシリカモデルと、前記各シリカモデルの周りを環状に取り囲みかつ前記マトリックスモデルよりも硬い物性が定義された界面モデルとを含み、前記シリカモデルは、前記界面モデルを介して紐状に連結された複数個を含み、前記ゴム材料モデルを設定するステップは、前記界面モデルに物性のパラメータを決定するパラメータ決定工程を含み、前記パラメータ決定工程は、解析対象となる前記ゴム材料モデルと同じ配合を有する第１の未加硫ゴム組成物と、該第１の未加硫ゴム組成物から界面結合剤を除いた配合の第２の未加硫ゴム組成物とを少なくとも準備する工程と、前記第１及び第２の未加硫ゴム組成物をそれぞれ溶剤に浸漬することにより、該第１及び第２の未加硫ゴム組成物からそれぞれマトリックスゴムを除去した第１及び第２の残留物を得る工程と、少なくとも前記第２の残留物のｔａｎδのピーク温度Ｔ２と、前記第１の残留物のｔａｎδのピーク温度Ｔ１とを測定し、これらに基づいて前記ピーク温度の差Ｔ２−Ｔ１を推定する工程と、前記第１の未加硫ゴム組成物からシリカを除いた配合の基本の加硫ゴム材料を少なくとも含むとともに該基本の加硫ゴム材料とは架橋密度のみが異なる複数種類の加硫ゴム材料を準備する工程と、前記各加硫ゴム材料のｔａｎδのピーク温度を測定し、前記加硫ゴム材料のｔａｎδのピーク温度と架橋密度との関係を得る工程と、前記関係から、ｔａｎδのピーク温度が前記基本の加硫ゴム材料のｔａｎδのピーク温度Ｔ３に前記差Ｔ２−Ｔ１を加えた温度に等しい加硫ゴム材料の架橋密度を特定する工程と、前記特定された架橋密度の加硫ゴム材料の物性に基づいて前記界面モデルの前記パラメータを定義する工程とを含むことを特徴する。
【００１２】
また請求項２記載の発明は、前記紐状に連結されたシリカモデルは、任意の軸線方向に関して、ゴム材料モデルの領域の全範囲を連続してのびる請求項１記載のゴム材料のシミュレーション方法である。
【発明の効果】
【００１４】
本発明によれば、ゴム材料モデルは、ゴムマトリックスをモデル化したマトリックスモデルと、該マトリックスモデル中に配置されかつシリカをモデル化した複数個のシリカモデルと、前記各シリカモデルの周りを環状に取り囲みかつ前記マトリックスモデルよりも硬い物性が定義された界面モデルとを含む。そして、前記シリカモデルは、界面モデルを介して紐状に連結された複数個を含む。このようなゴム材料を用いてシミュレーション（変形計算）を行うことにより、実際のシリカ配合ゴム材料の変形挙動と非常に相関の高い計算結果を得ることができる。
【図面の簡単な説明】
【００１５】
【図１】本実施形態で用いたコンピュータ装置の一例を示す斜視図である。
【図２】本実施形態の処理手順を示すフローチャートである。
【図３】ゴム材料モデル（微視構造）の一実施形態を示す線図である。
【図４】（Ａ）は粘弾性材料、（Ｂ）はその分子鎖１構造を説明する線図、（Ｃ）は１本の分子鎖の拡大図、（Ｄ）はセグメントの拡大図である。
【図５】（Ａ）は粘弾性材料の網目構造体を示す斜視図、（Ｂ）は８鎖モデルの一例を示す斜視図である。
【図６】（Ａ）、（Ｂ）は分子鎖の接合点の破断を説明する線図である。
【図７】マトリックスモデル及び界面モデルの物性の一例を示す応力−ひずみ曲線である。
【図８】ゴム材料モデル（微視構造）の他の実施形態を示す線図である。
【図９】シリカを配合していない加硫ゴムのｔａｎδと温度との関係を示すグラフである。
【図１０】変形シミュレーションの手順を示すフローチャートである。
【図１１】均質化法を説明する微視構造と全体構造との関係を示す。
【図１２】（ａ）〜（ｅ）は、ゴム材料モデルの実施例、比較例を示す線図である。
【図１３】実施例の８鎖モデルの構成図である。
【図１４】シミュレーションの結果と実験値とを比較する公称応力−伸び曲線である。
【図１５】従来のシリカ配合ゴム材料のモデルである。
【発明を実施するための形態】
【００１６】
以下、本発明の実施の一形態が図面に基づき説明する。
図１には、本発明のシミュレーション方法を実施するためのコンピュータ装置１が示されている。このコンピュータ装置１は、本体１ａ、キーボード１ｂ、マウス１ｃ及びディスプレイ装置１ｄを含んで構成されている。本体１ａの内部には、ＣＰＵ、ＲＯＭ、作業用メモリー及び磁気ディスク等のなどの大容量記憶装置が設けられる。また、本体１ａには、ＣＤ−ＲＯＭやフレキシブルディスクのドライブ装置１ａ１、１ａ２が設けられる。そして、前記大容量記憶装置には後述する本発明のシミュレーション方法を実行するための処理手順（プログラム）が記憶されている。
【００１７】
本実施形態のシミュレーション方法では、ゴムと、シリカと、これらを結合する界面結合剤とを含むシリカ配合ゴム材料の変形がシミュレーションされる。
【００１８】
前記界面結合材として、本実施形態では、シランカップリング剤が用いられる。シランカップリング剤としては、例えば、ビス（３−トリエトキシシリルプロピル）ポリスルフィド、ビス（２−トリエトキシシリルエチル）ポリスルフィド、ビス（３−トリメトキシシリルプロピル）ポリスルフィド、ビス（２−トリメトキシシリルエチル）ポリスルフィド、ビス（４−トリエトキシシリルブチル）ポリスルフィド、ビス（４−トリメトキシシリルブチル）ポリスルフィドなどが挙げられ、これらのシランカップリング剤は単独で用いてもよく、２種以上を組み合わせて用いても良い。なかでも、シランカップリング剤の添加効果およびコストの両立から、ビス（３−トリエトキシシリルプロピル）ジスルフィドなどが好適に用いられる。
【００１９】
図２には、本実施形態のシミュレーション方法の処理手順の一例が示される。本実施形態では、先ず、ゴム材料のモデルが設定される（ステップＳ１）。
【００２０】
図３には、繰り返し最小単位の微視構造としてのゴム材料モデル２の一例が視覚化して示されている。該ゴム材料モデル２は、解析しようとするシリカ配合ゴム材料の微小領域が、有限個の小さな要素（メッシュ）２ａ、２ｂ、２ｃ…に置き換えられたものである。各要素２ａ、２ｂ、２ｃ…は、数値解析が可能に定義される。
【００２１】
前記数値解析が可能とは、例えば有限要素法、有限体積法、差分法又は境界要素法といった数値解析法により、各要素ないし系全体についての変形計算が可能なことを意味する。具体的には、各要素２ａ、２ｂ、２ｃ…について、座標系における節点座標値、要素形状、材料特性などが定義される。各要素２ａ、２ｂ、２ｃ…には、例えば２次元平面としての三角形ないし四辺形の要素、３次元要素としては、例えば４ないし６面体の要素が好ましく用いられる。これにより、ゴム材料モデル２は、前記コンピュータ装置１にて取り扱い可能な数値データを構成する。
【００２２】
この実施形態のゴム材料モデル２は、後述する変形シミュレーションにおいて平面ひずみ状態の解析、さらに詳しくはｙ軸方向の引張変形シミュレーションが行われる。つまり、本実施形態では、ｚ方向にはひずみを持たない２次元のシミュレーションが行われる。この実施形態において、微視構造としてのゴム材料モデル２は、例えば３００nm×３００nmの正方形である。
【００２３】
また、本実施形態のゴム材料モデル２は、ゴムマトリックス部分がモデル化されたマトリックスモデル３と、このマトリックスモデル３の中に配置されかつシリカがモデル化されたシリカモデル４と、各シリカモデル４の周りを環状に取り囲みかつ前記マトリックスモデル３よりも硬い物性が定義された界面モデル５とを含む。
【００２４】
前記マトリックスモデル３は、ゴム材料モデル２の主要部を構成し、かつ、例えば三角形ないし四辺形の複数個の要素を用いて表現されている。変形計算を行うために、マトリックスモデル３を構成する各要素には、その物性として応力と伸びとの関係を表す関数が定義される。本実施形態のゴム材料のシミュレーション方法では、ゴム弾性応答を表現するために、前記マトリックスモデル３及び界面モデル５のゴム部分は、いずれも分子鎖網目理論に基づいて計算が行われる。
【００２５】
分子鎖網目理論について念のため述べると、図４（Ａ）、（Ｂ）に示されるように、連続体としてのゴム材料ａは、微視構造として、無秩序に配向された分子鎖ｃが接合点ｂで連結された網目構造を持つとの考えを前提とするものである。接合点ｂは、例えば分子間の化学的結合であってそれには架橋点などが含まれる。
【００２６】
分子鎖網目理論では、接合点ｂが原子の揺らぎ周期に対して長時間的には平均位置が変化しないものとし、接合点ｂの回りの摂動が無視される。さらに、二つの接合点ｂ、ｂを両端に持つ分子鎖ｃの端−端ベクトル（end-to-end vector ）は、それが埋め込まれているゴム材料の連続体と共変形するとの仮定が置かれる。
【００２７】
また、１本の分子鎖ｃは、図４（Ｃ）に示されるように、複数のセグメントｄから構成されるものとする。各セグメントｄは、化学的には同図（Ｄ）に示すように、炭素原子が共有結合によって連結した複数個のモノマーｆが連結したものに相当する。個々の炭素原子は、原子同士の結合軸の周りで互いに自由に回転しうるため、セグメントｄは全体として曲がりくねるなど様々な形態をとることができる。モノマーｆの数が十分多ければ、スケーリング則によってセグメントｄの巨視的な性質は変わらないので、一つのセグメントｄは、分子鎖網目理論において繰り返しの最小構成単位として取り扱われる。さらに、図４（Ｃ）に示される二つの接合点ｂ、ｂによって定義される分子鎖の形態は非ガウス統計分布に従うものとされる。従って、二つの接合点ｂ、ｂを結ぶ方向に伸び（ストレッチ）λを加えた場合に生じる応力σは、下式（１）で表すことができる。
【数１】
0005559594.tif 000002
【００２８】
そして、網目構造の全体的な応答特性は、個々の分子鎖の寄与を考慮して得ることができるが、それら全てを正確に考慮して計算するのは数学的に困難である。このため、分子鎖網目理論では、平均化手法が導入される。本実施形態では、８鎖モデルを採用した分子鎖網目理論に準拠している。即ち、この手法では、図５（Ａ）に示されるように、超弾性体であるゴム材料は、巨視的には、微小な８鎖モデルｇ…が集合した立方体状の網目構造体ｈであると仮定される。また、一つの８鎖モデルｇは、図５（Ｂ）に拡大して示されるように、立方体の中心に定められた一つの接合点ｂ１から、各頂点に設けられた８つの各接合点ｂ２にそれぞれ分子鎖ｃがのびているものと仮定して計算が行われる。
【００２９】
また、本実施形態では、分子鎖網目理論に、材料のひずみに応じた接合点ｂの消滅が考慮される。現実のゴム材料において、荷重の負荷における変形過程では、分子鎖の互いに絡み合った部分（即ち、前記接合点ｂ）が大きなひずみによって消滅することが知られている。つまり、接合点ｂの数が減少する。例えば、図６（Ａ）に示されるように、一つの接合点ｂで接合されている分子鎖ｃ１ないしｃ４に矢印方向の引張応力が作用すると、各分子鎖ｃ１ないしｃ４の伸びにより接合点ｂは大きなひずみを受けて消滅する。この結果、図６（Ｂ）に示されるように、これまで２本であった分子鎖ｃ１及びｃ２は、１本の長い分子鎖ｃ５になる。分子鎖ｃ３及びｃ４についても同様である。このような現象は、ゴム材料の負荷変形が進むにつれて逐次発生する。また、接合点ｂの数の減少は、１本の分子鎖ｃに含まれるセグメントの数を増加させることになる。
【００３０】
そして、上記の現象が、図５（Ａ）に示した網目構造体ｈに適用される。該網目構造体ｈは、幅方向、高さ方向及び奥行き方向にそれぞれ８鎖モデルｇがｋ個結合されているものとする。網目構造体ｈに含まれる接合点ｂの総数を「からみ数」として符号ｍで表すと、それは式（２）で表される。同様に、網目構造体ｈに含まれる分子鎖ｃの数（即ち、マトリックスモデル３の単位体積中に含まれる分子鎖の数）ｎは、式（３）で表される。
ｍ＝（ｋ＋１）3＋ｋ3        …式（２）
ｎ＝８ｋ3                   …式（３）
【００３１】
ここで、ｋは十分に大きい数とすると、上式（２）からｋの３次項以外を省略して次式（４）が得られる。さらに式（３）及び（４）の関係から、からみ数ｍは、ｎを用いて式（５）で表される。
ｍ＝２ｋ3   …式（４）
ｍ＝ｎ／４  …式（５）
【００３２】
さらに、網目構造体ｈは、変形の前後で材料の出入りが無いので、そのセグメントの総数ＮA は常に一定と仮定できる。従って、式（６）及び（７）が成り立つ。
ＮA ＝ｎ・Ｎ              …式（６）
Ｎ＝ＮA ／ｎ＝ＮA ／４ｍ  …式（７）
【００３３】
ここで、上述の接合点ｂの消滅は、マトリックスモデル３におけるからみ数ｍを減少させ、ひいては１本の分子鎖ｃに含まれるセグメントの数を増加させることに相当する。従って、８鎖モデルを用いた分子鎖網目理論に、上述の負荷変形（伸び）に伴った接合点の減少を導入するために、前記式（１）における分子鎖１本当たりの平均セグメント数Ｎを、ひずみに関するパラメータλｃに応じて増大させる。このひずみに関するパラメータλｃとしては、例えば、ひずみ、伸び、ひずみ速度又はひずみの１次の不変量などが挙げられる。また、上記負荷変形時とは、微小時間の間でモデルのひずみが増大する変形であり、逆に除荷変形時とは、ひずみが減少する変形とする。
【００３４】
以上より、本実施形態のマトリックスモデル３では、１本の分子鎖当たりの平均セグメント数Ｎは、下式（８）のように、１本の分子鎖当たりの初期セグメント数Ｎ0に、材料のひずみに応じたパラメータｆ（λc）の項が加算された関数として定義される。
Ｎ（λc）＝Ｎ0＋ｆ（λc）  …（８）
【００３５】
従って、上記式（８）で得られる平均セグメント数Ｎを用いた前記式（１）は、マトリックスモデル３のひずみに依存して応力を変化させることができる。そして、後述するゴム材料モデル２の変形シミュレーションでは、マトリックスモデル３の各要素について、負荷変形時においては前記パラメータλc が常時計算され、それらは式（８）に代入される。これにより、当該要素のセグメント数Ｎが常に更新されてシミュレーションに取り込まれる。
【００３６】
次に、前記シリカモデル４は、シリカを四辺形の複数個の要素を用いてモデル化したもので、全体として円形に形成されている。本実施形態において、各シリカモデル４の粒子径（直径）は全て等しく設定されている。なお、三次元モデルの場合、シリカモデル４は、球形にモデル化されるのが望ましい。また、シリカは、直径約１０〜３００nm程度であり、ゴムに比べて非常に硬い粒子からなる。シリカモデル４には、このような解析対象となるシリカの物性とほぼ等しい物性が設定される。即ち、本実施形態において、シリカモデル４は、粘弾性体ではなく弾性体として取り扱われる。また、シリカモデル４の粒子の個数は、例えば、解析対象のゴム材料のシリカ配合量に基づいて決定することができる。
【００３７】
前記界面モデル５は、シリカとマトリックスゴムとを化学的に結合させるシランカップリング剤の働きをシミュレーションに取り込むためにモデル化したものである。本実施形態の界面モデル５は、シリカモデル４の周りを小さい厚さｔで環状に連続して取り囲むように設定されている。従って、界面モデル５の内周面はシリカモデル４の外周面に接触している。本実施形態では、界面モデル５の内周面とシリカモデル４の外周面とは、互いに剥離しない条件が設定されるが、必要に応じて、予め定めた値以上の応力が生じたときに、シリカモデル４と界面モデル５との境界を分離させるような条件が設定されても良い。なお、界面モデル５の外周面の大部分は、マトリックスゴムモデル３に接触している。
【００３８】
前記界面モデル部６の厚さｔは、特に限定されるものではないが、種々の実験結果などに鑑み、シリカモデル４の直径の１０〜３０％程度、より好ましくは１５〜２５％程度に設定されるのが実際のゴム材料と整合する点で望ましい。
【００３９】
また、界面モデル５にも前記式（１）に準じた応力と伸びとの関係が定義される。しかし、界面モデル５で用いられる式（１）の平均セグメント数Ｎは、固定値Ｎ0に設定され、式（８）で採用されるようなｆ（λｃ）の項は含まれない。
【００４０】
また、実際の界面結合材の物性に鑑み、この界面モデル５はマトリックスゴムよりも硬い物性が定義される。本実施形態では、式（８）の平均セグメント数Ｎの初期値Ｎ0が、マトリックスモデル３に適用される初期値Ｎ0よりも小さく設定されている。これにより、界面モデル５は、マトリックスモデル３よりも伸び難く定義される。つまり、界面モデル５には、マトリックスモデル４よりも硬い物性（正確には粘弾性特性）が定義される。ただし、界面モデル５は、シリカモデル４よりは軟らかいのは言うまでもない。
【００４１】
また、界面モデル５には、マトリックスモデル３よりも小さい限界ストレッチが定義されるのが望ましい。実際のシリカ配合ゴム材料では、シリカ粒子間にも多くの界面結合剤が存在していると考えられるため、その部分の架橋密度がゴムマトリクス部分に比べて高いと推察することは合理的である。この状態をシミュレーションに反映させるために、図７に示されるように、界面モデル５については、応力が立ち上がり始めるひずみ（限界ストレッチ）を、マトリックスモデル３のそれよりも小さく設定するのが望ましい。
【００４２】
さらに、ゴム材料モデル２は、界面モデル５を介して紐状に連結された複数個のシリカモデル４（以下、このようなシリカモデルの連結体を「紐状体６」と呼ぶ場合がある。）を含んでいる。発明者らは、実際のシリカ配合ゴム中でのシリカの配列を完全に解明し得てはいないが、透過型電子顕微鏡による測定等の種々の解析でシリカ配置を特定したところ、概ね、シリカ配合ゴムのシリカ粒子は、その約３０％程度が界面結合材を介して他のシリカと紐状に連結されたネットワーク構造を有することを知見した。これは、界面結合剤同士が縮合反応により結合したためと考えられる。また、縮合反応を起こした界面結合剤は、架橋剤として働き、シリカ周辺の架橋密度を高めることも判明している。従って、上述のようなネットワーク構造は、シリカ粒子間の距離を小さく維持してゴムの変形に対する抵抗となり、シリカ粒子自体の補強効果との相乗作用により、ゴム全体をさらに固くする役割を担っていると考えられる。従って、ゴム中のシリカ粒子の配置をゴム材料モデル２に正確に取り込むことは、さらに精度の高いシミュレーションを実現するのに役立つ。この効果については、後の実施例で明らかにする。
【００４３】
本実施形態において、前記紐状体６は、任意の軸線方向に関して、ゴム材料モデル２の領域の全範囲を連続してのびている。即ち、この実施形態では、各紐状体６は、引張シミュレーションの引張方向であるｙ軸方向に沿って直線状に上下にのびる。また、紐状体６は、ｘ軸方向に間隔を開けて複数本が設けられるのが望ましい。なお、紐状体６が「ゴム材料モデル２の領域の全範囲を連続してのびる」とは、紐状体６の両端それぞれの界面モデル５ａ、５ｂが微視構造１の任意の方向（この例ではｙ軸方向）の端部２Ａ、２Ｂに接しており、その間で界面モデル５が途切れることなく連結されていることを意味している。
【００４４】
紐状体６の形状は、上記実施形態に限定されるものではない。例えば、図８に示されるように、紐状体６は、微視構造のセル中のｙ軸方向の全領域を右下がりにのびる第１部分７と、この第１部分７に交差してｙ軸方向の全領域を左下がりにのびる第２部分８とを含む略Ｘ字状のランダムな配置をとることもできる。また、紐状体６は、さらに複雑な屈曲形状をとることもできる。
【００４５】
また、界面モデル５にも、マトリックスゴム３及びシリカモデル４同様、その物性についてのパラメータ（密度、弾性率等）が入力される。しかしながら、界面モデル５の物性を直接測定して調べることは現実的には困難である。従来では、シミュレーションの解析対象となるゴム材料モデルと同等のシリカ配合ゴムと同じような変形挙動（例えば同じような引張試験の結果が得られるような挙動）となるように、界面モデル５の物性に関するパラメータが設定されている。しかしながら、このようなパラメータの入力では、試験条件が異なると、計算精度が大きく悪化する場合がある。そこで、本実施形態では、界面モデル５の物性については、以下で説明されるような工程に従って、界面モデル５の物性のパラメータを決定している。
【００４６】
先ず、解析対象のゴム材料モデルと同じ配合を有する第１の未加硫ゴム組成物と、該第１の未加硫ゴム組成物から界面結合剤のみを除いた配合を有する第２の未加硫ゴム組成物とが少なくとも準備される。なお、第１の未加硫ゴム組成物には、界面結合剤が８ｗｔ％及び硫黄が１ｐｈｒ含まれているものとする。第１の未加硫ゴム組成物と第２の未加硫ゴム組成物とは、界面結合剤の有無のみで異なっており、残余の配合は同一である。本実施形態では、より好ましい態様として、第１の未加硫ゴム組成物とは異なる配合量の界面結合剤を含む複数種類の他の未加硫ゴム組成物も準備される。これらの各ゴム組成物は、バンバリーミキサーなどで十分かつ均一に混練された状態で準備される。
【００４７】
次に、各未加硫ゴム組成物をそれぞれ溶剤に浸漬することにより、少なくとも第１及び第２の未加硫ゴム組成物からそれぞれマトリックスゴムを除去した第１及び第２の残留物を得る工程が行われる。一例として、塊状の前記各未加硫ゴム組成物を例えば１５０メッシュ程度の金網カゴに入れ、これを室温環境下でトルエンに約４８時間浸漬させる。これにより、各組成物からマトリックスゴムがトルエン中に溶出する。従って、金網カゴの第１の未加硫ゴム組成物からは、シリカ、界面結合剤、硫黄及びその他の添加剤を含む第１の残留物が得られる。他方、金網カゴの第２の未加硫ゴム組成物や他の未加硫ゴム組成物からは、シリカ、硫黄及びその他の添加剤を含む第２の残留物が得られる。
【００４８】
なお、加硫ゴムでは、マトリックスゴムを溶剤に溶出させることができないので、上記工程では、加硫前のゴム組成物が用いられる。また、上記第１の残留物は、マトリックスゴムを実質的に又は完全に含まないため、前記ゴム材料モデル２からマトリックスモデル３を除いて得られるシリカモデル４及び界面モデル５の結合物と等価なものと見なすことができる。従って、後述するように第１の残留物を利用することで、より正確な界面モデル５のパラメータを設定することができる。
【００４９】
次に、本実施形態では、各残留物（本実施形態では、第１、第２及び他の残留物）について、例えば粘弾性試験が行われ、それぞれのｔａｎδのピーク温度が測定される。この測定では、例えば金網カゴに残存しているシリカと界面結合剤とを主体的に含む残留物をプレスしてバルク状のサンプルを得、このサンプルから加硫ゴムの粘弾性試験で用いられるのと同様の試験片を切り出して測定が行われる。なお、ｔａｎδの測定条件は、次の通りである。
初期ひずみ：１０％
変形モード：引張
周波数：１０Hz
片振幅：１％
【００５０】
次に、本実施形態では、上記測定されたｔａｎδのピーク温度の結果等に基づいて前記ピーク温度の差Ｔ２−Ｔ１を推定する工程が行われる。図９（ａ）は、各残留物について、ｔａｎδのピーク温度と、界面結合材の配合量との関係を示すグラフである。この実施形態では、得られたデータについて線形近似を行うことにより、残留物のｔａｎδのピーク温度と界面結合剤との関係を示す近似直線が計算されている。近似直線は、ｙ＝０．８１６７ｘ−３８．３９６で表される。近似直線の精度を高めるためにも、本実施形態のように他の残留物を含むのが望ましい。また、前記近似直線を用いると、第２の残留物のｔａｎδのピーク温度Ｔ２と第１の残留物のｔａｎδのピーク温度Ｔ１との差Ｔ２−Ｔ１は、約６．５３３６℃と推定できる。なお、温度の差Ｔ２−Ｔ１を推定する工程は、このような線形の近似直線を用いる場合のみならず、二次以上の近似曲線を用いても良いのは言うまでもない。
【００５１】
次に、前記第１の未加硫ゴム組成物からシリカを除いた配合の基本の加硫ゴム材料が準備される。また、この基本の加硫ゴム材料の他にも、該基本の加硫ゴム材料とは架橋密度のみが異なる複数種類の加硫ゴム材料が準備される。そして、これらの各加硫ゴム材料について、例えば粘弾性試験を行い、前記各加硫ゴム材料のｔａｎδのピーク温度が測定される。測定条件は、上記の通りである。
【００５２】
図９（ｂ）は、各加硫ゴム材料について、ｔａｎδのピーク温度と、架橋密度（これは、硫黄量で便宜的に表される。）との関係を示すグラフである。この実施形態では、得られたデータについて線形近似を行うことにより、加硫ゴム材料のｔａｎδのピーク温度と架橋密度（硫黄量）との関係を示す近似直線が計算されている。近似直線は、ｙ＝３．０１５９ｘ−３７．９０６で表される。なお、この加硫ゴム材料についても、近似直線の精度を高めるために、サンプル数を多く確保することが望ましい。
【００５３】
次に、前記加硫ゴム材料のｔａｎδのピーク温度と架橋密度（硫黄量）との関係から、ｔａｎδのピーク温度が前記基本の加硫ゴム材料（硫黄量：１ｐｈｒ）のｔａｎδのピーク温度Ｔ３に前記差Ｔ２−Ｔ１を加えた温度（Ｔ３＋Ｔ２−Ｔ１）に等しい加硫ゴム材料の架橋密度（硫黄量）が特定される。この実施形態において、基本の加硫ゴム材料は、硫黄量１ｐｈｒであるため、そのｔａｎδのピーク温度Ｔ３は前記近似式から−２８．３５６５℃になる。従って、温度（Ｔ３＋Ｔ２−Ｔ１）は、−２８．３５６５℃になる。従って、この温度を前記近似式に代入して変数ｘを求めると、その値は約３．１７になる。この値は、ｔａｎδのピーク温度が、温度Ｔ３＋Ｔ２−Ｔ１である加硫ゴム材料の硫黄量（ｐｈｒ）と推定（特定）できる。
【００５４】
従って、上記で特定された架橋密度（硫黄量）の加硫ゴム材料の物性を調べ、それに基づいて前記界面モデル５の物性に関する前記パラメータを定義することができる。
【００５５】
このように、本実施形態では、硫黄を１ｐｈｒ含む残留物の界面結合剤を０wt％（界面結合剤の影響なしでマトリックスゴムと性質同等）から８wt％としたときのｔａｎδのピーク温度の上昇量Ｔ２−Ｔ１と、シリカ及び界面結合材を含まずかつ硫黄の配合量が１phrからいくらか増量させたとき（この例では２．１７ｐｈｒの増量させたとき）の加硫ゴムのｔａｎδのピーク温度の上昇量とが一致するとき、その加硫ゴムの物性を界面モデルの物性として定義している。即ち、シリカの界面結合剤に由来している界面モデルの架橋密度を、硫黄量を用いて近似的に見積もっている。種々の実験の結果、界面結合剤を配合することによる残留物のｔａｎδのピーク温度の上昇量と、硫黄を配合することによる加硫ゴムのｔａｎδのピーク温度の上昇量とは相関があることが判明している。従って、このようにしてパラメータを決定することでより精度の良い計算結果が得られる。
【００５６】
よって、上述の加硫ゴム材料が特定できれば、その物性を種々測定し、それらの値に基づいて前記界面モデル５のパラメータを容易に定義することができる。例えば、上記工程で求まった架橋密度（平均架橋密度）となるように、シリカを配合していない架橋ゴムを別途作成し、それについて再度、物性試験を行い、その物性を界面モデル５の物性値のパラメータとして決定することができる。
【００５７】
次に、ゴム材料モデル２を変形させるための変形条件が設定される（ステップＳ２）。本実施形態では、図３のｙ方向に任意の平均ひずみ速度を加えてゴム材料モデル２に引張変形を与える条件が定義される。ただし、変形条件は種々定めうるのは言うまでもない。
【００５８】
次に本実施形態のシミュレーション方法では、上述のように設定されたゴム材料モデル２を用いて変形シミュレーションが行われる（ステップＳ３）。変形シミュレーションの具体的な処理手順は、図１０に示される。変形シミュレーションでは、先ずゴム材料モデル２の各種のデータがコンピュータ装置１に入力される（ステップＳ３１）。入力されるデータには、各要素に定義された節点の位置や材料特性といった情報が含まれる。
【００５９】
コンピュータ装置１では、入力されたデータに基づいて各要素の剛性マトリックスを作成し（ステップＳ３２）、しかる後、全体構造の剛性マトリックスが組み立てられる（ステップＳ３３）。全体構造の剛性マトリックスには、既知節点の変位、節点力が導入され（ステップＳ３４）、剛性方程式の解析が行われる。そして、未知節点変位が決定され（ステップＳ３５）、各要素のひずみ、応力、主応力といった物理量を計算し、出力する（ステップＳ３６ないし３７）。
【００６０】
ステップＳ３８では、計算を終了させるか否かの判定がなされ、否定的である場合には、ステップＳ３２以降が繰り返される。このようなシミュレーション（変形計算）は、例えば有限要素法を用いたエンジニアリング系の解析アプリケーションソフトウエア（例えば米国リバモア・ソフトウェア・テクノロジー社で開発・改良されたＬＳ−ＤＹＮＡ等）を用いて行うことができる。
【００６１】
また、本シミュレーションは、均質化法（漸近展開均質化法）に基づいて行われる。均質化法は、図１１に示されるように、図３に示した微視構造（均質化法では「ユニットセル」とも呼ばれる）を周期的に持っているゴム材料全体Ｍを表現するｘI と、前記微視構造を表現するｙI との独立した２変数が用いられる。微視的スケールと巨視的スケールという異なる尺度の場におけるそれぞれ独立した変数を漸近展開することにより、図３に示した微視構造のモデル構造を反映させたゴム材料全体の平均的な力学応答を近似的に求めることができる。
【００６２】
前記変形計算が行われると、その結果から必要な物理量を取得することができる（ステップＳ４）。物理量としては、シリカ配合ゴムの変形挙動を調べるために、応力−ひずみ曲線が特に有効である。また、前記ゴム材料モデル２の各要素の時系列的な変形状態や物理量の分布を視覚化して表示することもできる。この際、各要素には、応力に応じた着色を施すことが望ましい。
【００６３】
そして、本実施形態では、上述のような界面モデル５を含むゴム材料モデル２を用いてシミュレーションを行うことにより、これまで以上に精度良い計算結果が得られる（後述の実施例参照）。これは、外力が加えられたときのシリカとマトリックスゴムとの界面挙動は勿論のこと、シリカ粒子間の動きが現実のゴム材料の内部での挙動と近似しているためと推測できる。
【００６４】
以上、本発明の実施形態について詳述したが、本発明は、種々の態様に変形して実施することができる。
【実施例】
【００６５】
１）解析モデル
図１２（ａ）〜（ｅ）に示した微視構造のモデル１〜５をそれぞれ周期的に持っているシリカ配合ゴム材料の全体巨視的モデルが設定された。
【００６６】
２）８鎖モデルの具体的構成式
この実施例では、ゴムの粘弾性挙動をさらに正確に記述するために、図１３に示されるように、粘弾性８鎖モデルとダンパーで構成されるモデルとを用いた。実際のゴムの分子鎖は、周囲の分子鎖との摩擦に起因した粘性を持っている。このような摩擦を表現するために、８鎖モデルＡの各分子鎖に、粘性抵抗をもつバネ・ダンパーの標準モデルが導入された。従って、単分子鎖の二つの接合点を結ぶ方向に伸びλc を加えた場合に生じる応力σcは、前記式（１）から次式で表すことができる。
【００６７】
【数２】
0005559594.tif 000003
【００６８】
また，図１３に示す単分子鎖の各要素のストレッチをλα、λβ、λγとすると、λα＝λｃである。その他の添え字α、β、γについても，図１３に示す要素と対応させている．ただし，λｃ＝λβ・λγである。
【００６９】
また、変形前の体積を基準とした単位体積当たりの仕事に相当するひずみエネルギー密度関数Ｗを用いると、応力σｃは、次式のように表される。
【００７０】
【数３】
0005559594.tif 000004
【００７１】
上記２つの式より、恒等的に次式が成り立つ。
【００７２】
【数４】
0005559594.tif 000005
【００７３】
また、８鎖モデルの場合、主ストレッチをλ1、λ2、λ3とすると、分子鎖のストレッチλｃは、√｛（λ12＋λ22＋λ32）／３｝と表すことができるので、下式が成り立つ。
【００７４】
【数５】
0005559594.tif 000006
【００７５】
上記３つの式より、上記８鎖モデルＡの主ストレッチ方向の応力σiAと、ストレッチλiとは次の関係で与えられる。
【数６】
0005559594.tif 000007
【００７６】
また、上記８鎖モデルＢの主ストレッチ方向の応力σiBと、ストレッチλi’とは次の関係で与えられる。
【数７】
0005559594.tif 000008
【００７７】
また、この実施例では、非圧縮性ゴム材料を扱うものとする。従って、非圧縮性を満たすために静水圧ｐを用いる。このとき、上の２つの式を用いると、非圧縮性ゴム材料の構成式は、次のようになる。
σi＝σiA＋σiB−ｐ
【００７８】
また、上記構成式の速度形式は、次のように表すことができ、この実施例ではこの式が用いられる。
【数８】
0005559594.tif 000009
【００７９】
３）マトリックスモデルの１本の分子鎖の平均セグメント数Ｎ
下式が採用された。
Ｎ（λc ）＝Ｎ0＋ｆ（λc）
＝Ｎ0＋ａ0＋ａ1λc＋ａ2λc2
ここで、ｆ（λc）の具体的な関数の決め方の一例について述べる。まず、マトリックスモデルの材料変形に伴うからみ点の変化は、λc（伸び）のみに依存する。従って、最初に、シリカ未充填の純ゴムでの応力−ひずみ曲線をうるシミュレーションを行ない、その結果が実際の純ゴムの試験結果と整合するようにｆ（λc）の係数が決定される。先ず、１回目の負荷変形サイクルでは、負荷時のからみ点数が変化（減少）し、除荷時は、除荷開始時の値から変化しないものとする。そして、シミュレーションの結果と、実際の試験結果とが整合するよう、前記関数ｆ（λc）の係数が決定される。なお、一旦、除荷された後の再負荷時においては、からみ点数の変化は不可逆的なものとする。即ち、再負荷時において、前回の負荷変形サイクルで到達した最大伸び以下の変形領域では、平均セグメント数Ｎは変化しないものとする。つまり、１回目の変形サイクルの除荷と、２回目の変形サイクルの再負荷でのＮの値は同一になる。そして、さらに変形が進み、１回目の変形サイクルで経験した最大伸びを超えた場合、Ｎの値が再び変化（減少）するものとして定められる。なお、本実施形態では上式の定数部などを次のように定めた。
ａ0＝１．２０
ａ1＝−２．１９
ａ2＝０．９９
［要素Ａ］
ＣRβＡ＝０．２２ＭＰａ
初期セグメントＮA0＝１４
総セグメント数ＮAα＝７．５４×１０26
［要素Ｂ］
ＣRαＢ＝０．２２ＭＰａ
初期セグメントＮB0＝１４
［粘弾性要素］
Ｃ1A＝５．０×１０5
Ｃ2A＝−０．５
ｍA＝３．５
Ｃ1D＝３．０×１０5
Ｃ2D＝−０．５
ｍD＝５．５
【００８０】
４）界面モデル
次のように定数等を設定した。
Ｎ＝Ｎ0
初期セグメント数（Ｎ＝Ｎ0）：８．０（固定値）
総セグメント数ＮA：７．５４×１０26
界面部の厚さ：シリカモデルの径の１０％
第１の接続部の幅Ｔ１：シリカモデルの径の２０％
第２の接続部の幅Ｔ２：シリカモデルの径の３５％
ＣαRs ＝０．３８５ＭＰａ
【００８１】
５）シリカモデルなどのパラメータ
粒子径は全て同一
シリカモデルの縦弾性率：１００ＭＰａ
シリカモデルのポアソン比：０．３
シリカモデルの体積含有率μ：２０％
【００８２】
６）他の条件
材料温度：Ｔ＝２９６Ｋ
ｋB ＝１．３８０６６×１０-29
ペナルティ定数：１００
カップリング剤の含有率：８wt％
【００８３】
７）変形条件
上記巨視的モデルに一様な一軸引張変形を発生させるため、図１１のｘ2方向に一定の変形速度１００mm／minを加え、伸びが２０％になるまで変形を与えた。なお上記巨視的モデルは、厚さ方向（図３のＺ軸方向）に変化しないように８鎖モデルを用いた分子鎖網目理論に基づいて計算を行った。
【００８４】
８）計算結果
図１４には、実施例、比較例の各ゴム材料モデルの応力−伸び曲線について、シミュレーションと実験値とが記載されている。図１４から明らかなように、実施例のモデルは、いずれも従来例に比して、実験値と非常に相関が高いことが確認できた。特に、引張方向に沿って直線状でかつゴム材料モデルの領域の全範囲を連続してのびるモデル１及びモデル３の精度が特に良好であることが確認できた。一方、比較例のモデルは、実施例に比べると、実験値に比べて相関が低いことが確認できる。
【符号の説明】
【００８５】
１  コンピュータ装置
２  ゴム材料モデル
３  マトリックスモデル
４  シリカモデル
５  界面モデル
６  紐状体
【図１】
0005559594.tif 000010
【図２】
0005559594.tif 000011
【図３】
0005559594.tif 000012
【図４】
0005559594.tif 000013
【図５】
0005559594.tif 000014
【図６】
0005559594.tif 000015
【図７】
0005559594.tif 000016
【図８】
0005559594.tif 000017
【図９】
0005559594.tif 000018
【図１０】
0005559594.tif 000019
【図１１】
0005559594.tif 000020
【図１２】
0005559594.tif 000021
【図１３】
0005559594.tif 000022
【図１４】
0005559594.tif 000023
【図１５】
0005559594.tif 000024