(57)【要約】【課題】従来の有限要素モデルを用いた計算方法に比べて効率よく短時間に計算できる複合材料の材料定数の算出方法及び複合材料中の材料の体積分率の算出方法を提供する。【解決手段】複合材料の材料定数の算出方法では、母相中に、第１の材料が内側球状粒子として分散し、この内側球状粒子の周りに、第２の材料が被覆して構成された外側殻状層を有する仮想複合材料を定め、この仮想複合材料の有効材料定数を未知数とする非線形方程式を用意する。次に、方程式を解いて、仮想複合材料の有効材料定数を求める。方程式は、外側殻状層を囲む周囲の材料定数を、求めるべき前記仮想複合材料の有効の材料定数として定める。この方程式を用いて複合材料中の材料の体積分率を算出する。【選択図】  図１
2010250796.tif 000002
【特許請求の範囲】
【請求項１】
力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、複合材料の力学的有効材料定数をコンピュータが算出する方法であって、
前記母相中に、前記第１の材料が既知の体積分率で内側球状粒子として分散し、この内側球状粒子の周囲に前記第２の材料を被覆層として所定の厚さで被覆して構成された外側殻状層を有する仮想複合材料を、前記複合材料として定めることにより、この仮想複合材料の有効材料定数を未知数とする非線形方程式を用意するステップと、
前記用意した非線形方程式を解いて、前記仮想複合材料の有効材料定数を、前記複合材料の有効材料定数として求めるステップと、を有し、
前記非線形方程式は、前記仮想複合材料中の前記外側殻状層を囲む周囲の材料定数を、算出すべき前記仮想複合材料の有効材料定数として定めることにより得られる再帰型非線型方程式である、算出方法。
【請求項２】
前記仮想複合材料中の変位は、前記内側球状粒子と前記被覆層との境界で、表面張力及び変位が連続し、前記外側殻状層と前記母相との境界で、表面張力及び変位が連続するように、定めることにより、前記非線形方程式が用意される、請求項１に記載の算出方法。
【請求項３】
前記母相の材料定数をＣAと表し、前記第１の材料の材料定数をＣBと表し、前記第２の材料の材料定数をＣCと表し、前記第１の材料の体積分率をｖBと表し、前記第２の材料の体積分率をｖCと表し、前記仮想複合材料の有効材料定数をＣ*と表したとき、
前記非線形方程式は下記式で表され、
下記式中のＡB（ＣB，ＣC，Ｃ＊）及びＡC（ＣB，ＣC，Ｃ＊）はそれぞれ、前記仮想複合材料中の前記内側球状粒子を囲む周囲の材料定数を、前記第２の材料の材料定数ＣCとし、前記外側殻状層を囲む周囲の材料定数を、算出すべき前記仮想複合材料の有効材料定数として定めることにより得られ、前記材料定数Ｃ＊に関して非線形な式で定まる比例定数である、請求項１または２に記載の算出方法。
Ｃ*  ＝  ＣA＋ｖB・（ＣB−ＣA）・ＡB（ＣB，ＣC，Ｃ＊）
＋ｖC・（ＣC−ＣA）・ＡC（ＣB，ＣC，Ｃ＊）
【請求項４】
応力及びひずみをせん断成分と静水圧成分に分解し、前記静水圧成分に対応する前記仮想複合材料中の変位をｕi（i＝１、２または３）としたとき、
前記変位ｕiは、前記内側球状粒子及び前記外側殻状層の内部において、前記内側球状粒子及び前記外側殻状層の中心を基準とする位置ｘi（i＝１、２または３）に比例し、
前記変位ｕiは、前記内側球状粒子及び前記外側殻状層の外部において、前記中心を基準とする位置ｘi（i＝１、２または３）に比例する項と、前記中心からの距離の３乗に反比例し、前記位置ｘi（i＝１、２または３）に比例する項を有するように定めることにより、前記非線形方程式が用意される、請求項１〜３のいずれか１項に記載の算出方法。
【請求項５】
応力及びひずみを、せん断成分と静水圧成分に分解し、前記せん断成分に対応する前記仮想複合材料中の変位をｕi（i＝１、２または３）としたとき、
前記変位ｕiを、前記内側球状粒子及び前記外側殻状層の内部及び外部において、前記内側球状粒子及び前記外側殻状層の中心からの距離の０乗、距離の２乗、距離の−３乗及び距離の５乗のそれぞれに比例する項の加算式で定めることにより、前記非線形方程式が用意される、請求項１〜３のいずれか１項に記載の算出方法。
【請求項６】
前記非線形方程式は、ニュートン・ラフソン法を用いて、解が収束するように解かれる、請求項１〜５のいずれか１項に記載の算出方法。
【請求項７】
力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、前記第１の材料または前記第２の材料の体積分率をコンピュータが算出する方法であって、
前記複合材料の有効材料定数を実験結果から特定するステップと、
前記母相中に、前記第１の材料が内側球状粒子として分散し、この内側球状粒子の周囲に前記第２の材料を被覆層として所定の厚さで被覆して構成された外側殻状層を有する仮想複合材料を、前記複合材料として定めることにより、前記第１の材料または前記第２の材料の体積分率を未知数とする非線形方程式を用意するステップと、
前記用意した非線形方程式を解いて、前記第１の材料または前記第２の材料の体積分率を求めるステップと、を有し、
前記非線形方程式は、前記仮想複合材料中の前記外側殻状層を囲む周囲の材料定数を、求めるべき前記仮想複合材料の有効材料定数として定めることにより得られる再帰型非線型方程式である、算出方法。
【請求項８】
前記仮想材料中の変位は、前記内側球状粒子と前記被覆層との境界で、表面張力及び変位が連続し、前記外側殻状層と前記母相との境界で、表面張力及び変位が連続するように、定めることにより、前記非線形方程式が用意される、請求項７に記載の算出方法。
【請求項９】
前記母相の材料定数をＣAと表し、前記第１の材料の材料定数をＣBと表し、前記第２の材料の材料定数をＣCと表し、前記第１の材料の体積分率をｖBと表し、前記第２の材料の体積分率をｖCと表し、前記仮想複合材料の有効材料定数をＣ*と表したとき、
前記非線形方程式は下記式で表され、
下記式中のＡB（ＣB，ＣC，Ｃ＊）及びＡC（ＣB，ＣC，Ｃ＊）はそれぞれ、前記仮想複合材料中の前記内側球状粒子を囲む周囲の材料定数を、前記第２の材料の材料定数ＣCとし、前記外側殻状層を囲む周囲の材料定数を、算出すべき前記仮想複合材料の有効材料定数として定めることにより得られ、前記材料定数Ｃ＊に関して非線形な式で定まる比例定数である、請求項７または８に記載の算出方法。
Ｃ*  ＝  ＣA＋ｖB・（ＣB−ＣA）・ＡB（ＣB，ＣC，Ｃ＊）
＋ｖC・（ＣC−ＣA）・ＡC（ＣB，ＣC，Ｃ＊）
【請求項１０】
応力及びひずみを、せん断成分と静水圧成分に分解し、前記静水圧成分に対応する前記仮想材料中の変位をｕi（i＝１、２または３）としたとき、
前記変位ｕiは、前記内側球状粒子及び前記外側殻状層の内部において、前記内側球状粒子及び前記外側殻状層の中心を基準とする位置ｘi（i＝１、２または３）に比例し、
前記変位ｕiは、前記内側球状粒子及び前記外側殻状層の外部において、前記中心を基準とする位置ｘi（i＝１、２または３）に比例する項と、前記中心からの距離の３乗に反比例し、前記位置ｘi（i＝１、２または３）に比例する項を有するように定めることにより、前記非線形方程式が用意される、請求項７または８に記載の算出方法。
【請求項１１】
応力及びひずみを、せん断成分と静水圧成分に分解し、前記せん断成分に対応する前記仮想複合材料の変位をｕi（i＝１、２または３）としたとき、
前記変位ｕiを、前記内側球状粒子及び前記外側殻状層の内部及び外部において、前記内側球状粒子及び前記外側殻状層の中心からの距離の０乗、距離の２乗、距離の−３乗及び距離の５乗のそれぞれに比例する項の加算式で定めることにより、前記非線形方程式が用意される、請求項７または８に記載の算出方法。
【請求項１２】
前記非線形方程式は、２分法を用いて解かれる請求項７〜１２のいずれか１項に記載の体積分率の算出方法。
【請求項１３】
力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、複合材料の力学的有効材料定数を算出する、コンピュータの実行可能なプログラムを記録した記録メディアであって、
前記プログラムは、
前記母相中に、前記第１の材料が既知の体積分率で内側球状粒子として分散し、この内側球状粒子の周囲に前記第２の材料を被覆層として所定の厚さで被覆して構成された外側殻状層を有する仮想複合材料を、前記複合材料として定めることにより、この仮想複合材料の有効材料定数を未知数とする非線形方程式を用意する手順と、
前記用意した非線形方程式を解いて、前記仮想複合材料の有効材料定数を、前記複合材料の有効材料定数として求める手順と、をコンピュータに実行させ、
前記非線形方程式は、前記仮想複合材料中の前記外側殻状層を囲む周囲の材料定数を、求めるべき前記仮想複合材料の有効材料定数として定めることにより得られる再帰型非線型方程式である、記録メディア。
【請求項１４】
前記仮想複合材料中の変位は、前記内側球状粒子と前記被覆層との境界で、変位及び表面張力が連続し、前記外側殻状層と前記母相との境界で、変位及び表面張力が連続するように、定めることにより、前記非線形方程式が用意される、請求項１３に記載の記録メディア。
【請求項１５】
前記母相の材料定数をＣAと表し、前記第１の材料の材料定数をＣBと表し、前記第２の材料の材料定数をＣCと表し、前記第１の材料の体積分率をｖBと表し、前記第２の材料の体積分率をｖCと表し、前記仮想複合材料の有効材料定数をＣ*と表したとき、
前記非線形方程式は下記式で表され、
下記式中のＡB（ＣB，ＣC，Ｃ＊）及びＡC（ＣB，ＣC，Ｃ＊）はそれぞれ、前記仮想複合材料中の前記内側球状粒子を囲む周囲の材料定数を、前記第２の材料の材料定数ＣCとし、前記外側殻状層を囲む周囲の材料定数を、算出すべき前記仮想複合材料の有効材料定数として定めることにより得られ、前記材料定数Ｃ＊に関して非線形な式で定まる比例定数である、請求項１３または１４に記載の記録メディア。
Ｃ*  ＝  ＣA＋ｖB・（ＣB−ＣA）・ＡB（ＣB，ＣC，Ｃ＊）
＋ｖC・（ＣC−ＣA）・ＡC（ＣB，ＣC，Ｃ＊）
【請求項１６】
力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、前記第１の材料または前記第２の材料の体積分率を算出する、コンピュータの実行可能なプログラムを記録した記録メディアであって、
前記プログラムは、
前記複合材料の有効材料定数を実験結果から特定する手順と、
前記母相中に、前記第１の材料が内側球状粒子として分散し、この内側球状粒子の周囲に前記第２の材料を被覆層として所定の厚さで被覆して構成された外側殻状層を有する仮想複合材料を、前記複合材料として定めることにより、前記第１の材料または前記第２の材料の体積分率を未知数とする非線形方程式を用意する手順と、
前記用意した非線形方程式を解いて、前記第１の材料の体積分率を求める手順と、コンピュータに実行させ、
前記非線形方程式は、前記仮想複合材料中の前記外側殻状層を囲む周囲の材料定数を、求めるべき前記仮想複合材料の有効材料定数として定めることにより得られる再帰型非線型方程式である、記録メディア。
【請求項１７】
前記仮想複合材料中の変位は、前記内側球状粒子と前記被覆層との境界で、変位及び表面張力が連続し、前記外側殻状層と前記母相との境界で、変位及び表面張力が連続するように、定めることにより、前記非線形方程式が用意される、請求項１６に記載の記録メディア。
【請求項１８】
前記母相の材料定数をＣAと表し、前記第１の材料の材料定数をＣBと表し、前記第２の材料の材料定数をＣCと表し、前記第１の材料の体積分率をｖBと表し、前記第２の材料の体積分率をｖCと表し、前記仮想複合材料の有効材料定数をＣ*と表したとき、
前記非線形方程式は下記式で表され、
下記式中のＡB（ＣB，ＣC，Ｃ＊）及びＡC（ＣB，ＣC，Ｃ＊）はそれぞれ、前記仮想複合材料中の前記内側球状粒子を囲む周囲の材料定数を、前記第２の材料の材料定数ＣCとし、前記外側殻状層を囲む周囲の材料定数を、算出すべき前記仮想複合材料の有効材料定数として定めることにより得られ、前記材料定数Ｃ＊に関して非線形な式で定まる比例定数である、請求項１６または１７に記載の記録メディア。
Ｃ*  ＝  ＣA＋ｖB・（ＣB−ＣA）・ＡB（ＣB，ＣC，Ｃ＊）
＋ｖC・（ＣC−ＣA）・ＡC（ＣB，ＣC，Ｃ＊）
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、複合材料の力学的有効材料定数を算出する方法と、力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、複合材料における前記第１の材料または前記第２の材料の体積分率を算出する方法と、これらの方法をコンピュータに実行させるプログラムを記録した記録メディアに関する。
【背景技術】
【０００２】
従来より、所定の材料が母相中に含まれた複合材料の力学特性を精度よく推定することが盛んに試みられている。これは、複合材料が所望の特性になるように、どのような力学特性を持つ材料を、どの程度の体積分率で含めればよいか、実験上で探索する替わりに、コンピュータ上で効率よく探索するためである。この探索の結果、所望の材料に関する早期の配合設計が可能となる。
【０００３】
このような状況下、下記特許文献１には、複数の微小要素からなり３次元的に不均質な変形特性を有するミクロ構造が１方向にのみ周期的に配置されるマクロ構造体の構造解析方法が記載されている。当該文献では、マクロ構造体における周期性の単位であるユニットセルを特定し、ユニットセルの均質な材料特性を有するものとして均質化弾性係数を求める。この後、マクロ構造体を均質化弾性特性を有するものとしてモデル化し、周期性配置の１方向に沿った任意の位置における変形を求めるマクロスケール解析を行う。得られた周期性配置の１方向に沿った任意の位置における変形を、その位置におけるユニットセルを構成する各微小要素のそれぞれに与えてそれらの局所的応答を得るローカル解析を行う。
この構造解析方法により、断面内で不均質なマクロ構造の構造計算時間を短縮することができるとされている。
【０００４】
しかし、上記構造解析方法は、微小要素を用いて構成した有限要素モデルを用いる方法である。このため、上記方法は、モデルの作成及び計算に多くの時間を要し、迅速な初期設計及び初期開発において、有効な手段となりにくい、といった問題があった。
【０００５】
一方、複合材料の力学特性を求めるために、従来より、バネ及びダッシュポットを用いた古典的解析モデルも用いられる。このモデルは、計算時間が短く効率的であるが、複合材料のミクロ状態を考慮したモデルでないため、得られる結果の情報が少なく、精度も低い、といった問題もある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００７−１２２２４２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
そこで、本発明は、上記問題点を解決するために、従来の有限要素モデルを用いた計算方法に比べて効率よく短時間に計算できる複合材料の材料定数の算出方法及び複合材料中の材料の体積分率の算出方法を提供するとともに、上記方法をコンピュータで実行するプログラムを記録した記録メディアを提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的は、力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、複合材料の力学的有効材料定数をコンピュータが算出する、以下の方法で達成することができる。
すなわち、方法は、
（Ａ）前記母相中に、前記第１の材料が既知の体積分率で内側球状粒子として分散し、この内側球状粒子の周囲に前記第２の材料を被覆層として所定の厚さで被覆して構成された外側殻状層を有する仮想複合材料を、前記複合材料として定めることにより、この仮想複合材料の有効材料定数を未知数とする非線形方程式を用意するステップと、
（Ｂ）前記用意した非線形方程式を解いて、前記仮想複合材料の有効材料定数を、前記複合材料の有効材料定数として求めるステップと、を有する。
（Ｃ）その際、前記非線形方程式は、前記仮想複合材料中の前記外側殻状層を囲む周囲の材料定数を、算出すべき前記仮想複合材料の有効材料定数として定めることにより得られる再帰型非線型方程式である。
【０００９】
上記方法において、前記仮想複合材料中の変位は、前記内側球状粒子と前記被覆層との境界で、変位及び表面張力が連続し、前記外側殻状層と前記母相との境界で、変位及び表面張力が連続するように、定めることにより、前記非線形方程式が用意されることが好ましい。      
【００１０】
また、前記母相の材料定数をＣAと表し、前記第１の材料の材料定数をＣBと表し、前記第２の材料の材料定数をＣCと表し、前記第１の材料の体積分率をｖBと表し、前記第２の材料の体積分率をｖCと表し、前記仮想複合材料の有効材料定数をＣ*と表したとき、
前記非線形方程式は下記式で表されことが好ましい。  
下記式中のＡB（ＣB，ＣC，Ｃ＊）及びＡC（ＣB，ＣC，Ｃ＊）はそれぞれ、前記仮想複合材料中の前記内側球状粒子を囲む周囲の材料定数を、前記第２の材料の材料定数ＣCとし、前記外側殻状層を囲む周囲の材料定数を、算出すべき前記仮想複合材料の有効材料定数として定めることにより得られ、前記材料定数Ｃ＊に関して非線形な式で定まる比例定数である。
Ｃ*  ＝  ＣA＋ｖB・（ＣB−ＣA）・ＡB（ＣB，ＣC，Ｃ＊）
＋ｖC・（ＣC−ＣA）・ＡC（ＣB，ＣC，Ｃ＊）
【００１１】
さらに、上記目的は、力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、前記第１の材料または前記第２の材料の体積分率をコンピュータが算出する、以下の方法で達成することができる。
（Ｄ）前記複合材料の有効材料定数を実験結果から特定するステップと、
（Ｅ）前記母相中に、前記第１の材料が内側球状粒子として分散し、この内側球状粒子の周囲に前記第２の材料を被覆層として所定の厚さで被覆して構成された外側殻状層を有する仮想複合材料を、前記複合材料として定めることにより、前記第１の材料または前記第２の材料の体積分率を未知数とする非線形方程式を用意するステップと、
（Ｆ）前記用意した非線形方程式を解いて、前記第１の材料または前記第２の材料の体積分率を求めるステップと、を有する。
（Ｇ）その際、前記非線形方程式は、前記仮想複合材料中の前記外側殻状層を囲む周囲の材料定数を、求めるべき前記仮想複合材料の有効材料定数として定めることにより得られる再帰型非線型方程式である。
【００１２】
上記方法において、前記仮想材料中の変位は、前記内側球状粒子と前記被覆層との境界で、表面張力及び変位が連続し、前記外側殻状層と前記母相との境界で、表面張力及び変位が連続するように、定めることにより、前記非線形方程式が用意されることが好ましい。
【００１３】
また、前記母相の材料定数をＣAと表し、前記第１の材料の材料定数をＣBと表し、前記第２の材料の材料定数をＣCと表し、前記第１の材料の体積分率をｖBと表し、前記第２の材料の体積分率をｖCと表し、前記仮想複合材料の有効材料定数をＣ*と表したとき、
前記非線形方程式は下記式で表されることが好ましい。
下記式中のＡB（ＣB，ＣC，Ｃ＊）及びＡC（ＣB，ＣC，Ｃ＊）はそれぞれ、前記仮想複合材料中の前記内側球状粒子を囲む周囲の材料定数を、前記第２の材料の材料定数ＣCとし、前記外側殻状層を囲む周囲の材料定数を、算出すべき前記仮想複合材料の有効材料定数として定めることにより得られ、前記材料定数Ｃ＊に関して非線形な式で定まる比例定数である。
Ｃ*  ＝  ＣA＋ｖB・（ＣB−ＣA）・ＡB（ＣB，ＣC，Ｃ＊）
＋ｖC・（ＣC−ＣA）・ＡC（ＣB，ＣC，Ｃ＊）
【００１４】
さらに、上記目的は、力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、複合材料の力学的の有効材料定数を、上述の算出方法を用いて算出する、コンピュータに実行させるプログラムを記録した記録メディアによって、達成することができる。
同様に、力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、前記第１の材料または前記第２の材料の体積分率を、上述の算出方法により算出する、コンピュータに実行させるプログラムを記録した記録メディアによって、達成することができる。
【００１５】
上述の算出方法に共通する再帰型非線型方程式は、後述する式（６）で表される方程式に基づく。具体的には、再帰型非線形方程式に用いる式（６）中の比例定数Ａは、仮想複合材料のひずみ場をナビエの方程式を利用して求めることにより定まる。そのとき、仮想複合材料中の球状粒子を囲む周囲の材料定数を、仮想複合材料の有効材料定数として定める。すなわち、セルフコンシステント近似を用いる。
【発明の効果】
【００１６】
本発明は、母相中に所定の材料が球状粒子として分散した仮想複合材料を、複合材料として定め、さらに、用いる非線形方程式は、仮想複合材料中の外側殻状層を囲む周囲の材料定数を、求めるべき前記仮想複合材料の有効材料定数として定めることにより得られる、解析的な再帰型非線型方程式である。このため、複合材料の有効材料定数の算出及び複合材料中の材料の体積分率の算出を効率よく短時間に行うことができる。
【図面の簡単な説明】
【００１７】
【図１】本発明の複合材料の材料定数の算出方法及び複合材料中の材料の体積分率の算出方法を実施する計算装置の一例の概略の構成を示す図である。
【図２】（ａ）〜（ｃ）は、複合材料の応力ひずみを説明する図である。
【図３】本発明の複合材料の材料定数の算出方法及び複合材料中の材料の体積分率の算出方法に用いる仮想複合材料を説明する図である。
【図４】本発明の複合材料の材料定数の算出方法の一例のフローを示す図である。
【図５】本発明の複合材料中の材料の体積分率の算出方法の一例のフローを示す図である。
【図６】複合材料の有限要素モデルを説明する図である。
【発明を実施するための形態】
【００１８】
以下、本発明の複合材料の材料定数の算出方法、及び複合材料中の材料の体積分率の算出方法を、下記図面に示す実施形態に基いて詳細に説明する。
図１は、複合材料の材料定数の算出方法、及び複合材料中の材料の体積分率の算出方法を実施する計算装置１０を示す。
【００１９】
計算装置１０は、ＣＰＵ１２、ＲＯＭ１４、ＲＡＭ１６及び入出力ポート１８を備えるコンピュータにより構成されている。計算装置１０は、ＲＯＭ１４に記録されたプログラムを起動することにより、条件設定モジュール２０と、非線形方程式設定モジュール２２と、方程式解法モジュール２４と、収束解判定モジュール２６と、結果処理モジュール２８と、を形成する。すなわち、上記モジュールは、ソフトウェアを起動することによって構成される各部分である。計算装置１０のモジュールは、例えば、数式処理ソフトウェアを用いて形成される。
計算装置１０は、入出力ポート１８を介して、記憶装置３０、マウスやキーボード等の入力装置３２、及び、プリンタやモニタ等の出力装置３４に接続されている。
【００２０】
計算装置１０は、
（１）力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、複合材料の力学的有効材料定数を算出する第１の処理と、
（２）力学的材料定数が既知の母相中に、力学的材料定数が既知の第１の材料及び第２の材料を含む複合材料に関して、前記第１の材料または前記第２の材料の体積分率を算出する第２の処理と、を選択的に実行する。ここで、複合材料の力学的有効材料定数とは、複合材料全体の力学的材料定数である。複合材料全体の力学材料定数とは、複合材料を均質な材料としてみた時の力学材料定数をいう。以降、力学的材料定数を材料定数という。
【００２１】
条件設定モジュール２０は、第１の処理及び第２の処理のどちらの処理を行うかを設定するとともに、第１の処理あるいは第２の処理に必要な情報を揃えて値を設定する部分である。
このモジュール２０は、第１の処理では、母相の材料定数と、各材料の材料定数、及び各材料の体積分率の各数値を設定する。材料定数は、例えば、ラメ定数λ，μの組、あるいは、ラメ定数λ，μの一方と体積弾性率Ｋの組を含む。体積弾性率Ｋは、ラメ定数λ，μとの間で、Ｋ＝λ＋２／３・μにより関係付けられるので、例えば、体積弾性率Ｋとラメ定数μを用いて材料定数を表すことができる。
また、このモジュール２０は、第２の処理では、実験により取得した複合材の有効材料定数、各材料の材料定数、及び前記材料のうちのいくつかの体積分率の数値を設定する。
これらの数値は、入力装置３２を介したオペレータの指示入力により、あるいは記憶装置３０に予め記憶された情報を呼び出すことにより、設定される。
【００２２】
非線形方程式設定モジュール２２は、条件設定モジュール２０で設定された材料定数及び体積分率の値を用いて、第１の処理あるいは第２の処理に応じて、非線形方程式を用意する部分である。このモジュールは、第１の処理及び第２の処理のいずれの場合も、複数の材料を母相中に含む複合材料を想定したとき、この複合材料の有効材料定数を算出するための再帰型非線形方程式、すなわち、後述する式（６）を呼び出し、この方程式の各係数に、材料定数及び体積分率に基づいた値を与える。これにより、解くべき未知数からなる方程式が用意される。
【００２３】
方程式解法モジュール２４は、第１の処理の場合、算出すべき複合材料の有効材料定数に値を付与したとき、この値の付与により、用意された再帰型非線形方程式を用いたニュートン・ラフソン法の漸化式に従って、収束解に近い複合材料の材料定数の値を再帰的に算出する部分である。すなわち、方程式解法モジュール２４は、最初に複合材料の有効材料定数の値として初期値（ｎ＝１）を設定し、この値を用いて、漸化式にしたがって、収束解に近い複合材料の有効材料定数の値（ｎ＝２）を再帰的に求める。この値は、後述する収束解判定モジュール２６に送られ、判定の対象となる。判定の結果、解が収束しない場合、算出された材料定数の値を用いて、上記漸化式に従って、さらに収束解に近い複合材料の有効材料定数の値（ｎ＝３）を再帰的に求める。このようにして、方程式解法モジュール２４は、収束解判定モジュール２６において、収束解が見出されるまで、複合材料の有効材料定数の値を再帰的に求める処理を繰り返す。
【００２４】
また、第２の処理の場合、方程式解法モジュール２４は、下記に示す関数ｆ（ｘ）（ｘは求めようとする所定の材料の体積分率）を定め、この関数ｆ（ｘ）の上限値ｘ＝ｘ1における値ｆ（ｘ1）と下限値ｘ＝ｘ2における値ｆ（ｘ2）との間で積ｆ（ｘ1）・ｆ（ｘ2）を算出し、この結果を収束判解判定モジュール２８に送る。
【００２５】
関数ｆ（ｘ）＝（実験により取得した複合材料の有効材料定数）−（再帰型非線形方程式（６）で求める体積分率ｘにおける複合材料の有効材料定数）
【００２６】
このような積ｆ（ｘ1）・ｆ（ｘ2）の算出は、収束解が見出されるまで、上限値ｘ＝ｘ1と下限値ｘ＝ｘ2における値を変更しながら繰り返し求められる。
【００２７】
収束解判定モジュール２６は、方程式解法モジュール２４の算出結果に基づいて収束解か否かを判定する部分である。
第１の処理の場合、繰返回数（ｎ＋１）で求めた複合材料の有効材料定数の値と、繰返回数ｎでもとめた複合材料の有効材料定数の値との差分の絶対値が予め設定された閾値を下回る場合、繰返回数（ｎ＋１）における複合材料の有効材料定数の値が収束解であるとされる。それ以外は、収束解は見出されないとして、繰返回数（ｎ＋１）で求めた値を用いて方程式解法モジュール２４で計算が行われるように、指示される。上述したように、再帰型非線型方程式は、係数等に具体的な数値が付与された方程式なので、この非線形方程式の微分した関数も既知である。したがって、収束解判定モジュール２６は、公知のニュートン・ラフソン法を用いて収束解を見出すことができる。
第２の処理の場合、ｆ（ｘ1）・ｆ（ｘ2）の値が、負か否かが調べられ、後述するように、２分法を用いて、上限値及び下限値の値が再設定される。
【００２８】
結果処理モジュール２８は、収束解判定モジュール２６で得られた収束解を用いて、複合材料の有効ヤング率やせん断剛性を算出し、あるいは、求めようとする材料の体積分率を取得する部分である。
【００２９】
出力装置３４は、結果処理モジュール２８で得られた各情報をプリント出力し、あるいは画面表示を行う部分である。
記憶装置３０には、予め設定された複合材料や複合材料に含まれる材料の材料定数等が記録保持されているデータベースを備えている。
【００３０】
このような処理は、複合材料の有効材料定数の算出を、以下に示す再帰型非線形方程式を用いて解析的に行うことにより実現される。以下、この再帰型非線形方程式について詳細に説明する。
【００３１】
図２（ａ）〜（ｃ）は、複合材料を説明する図である。
以下で説明する応力、ひずみ及び材料定数は、それぞれ２階テンソル、２階テンソル、及び４階テンソルであるが、わかり易く説明するためにスカラーとして説明する。
複合材料は、図２（ａ）に示すように、母相であるエポキシ樹脂（以降、エポキシという）ＰにフィラーＦ及びウレタン材料Ｕが含まれて構成されている。その際、フィラーＦの周りを一定の厚さで被覆するようにウレタンＵが覆われている。ここで、ｖP，ｖF，ｖUは、エポキシＰ、フィラーＦ及びウレタンＵの体積分率であり、ｖP＋ｖF＋ｖU＝１である。
【００３２】
この複合材料では、図２（ａ）に示すように、フィラーＦ及びウレタンＵに作用する応力σF，σUとエポキシＰに作用する応力σPとが体積分率を重み付け係数として重み付け加算されて、複合材料全体の平均応力σTとなる。これに対応して、図２（ｂ）に示すように、フィラーＦ及びウレタンＵに発生するひずみεF，εUとエポキシＰに発生するひずみεPとが体積分率を重み付け係数として重み付け加算されて、複合材料全体の平均ひずみεTとなる。平均応力σT及び平均ひずみεTは、下記式（１）、（２）で表される。
【００３３】
【数１】
2010250796.tif 000003
【数２】
2010250796.tif 000004
【００３４】
一方、複合材料全体における平均応力σTとひずみεTとの関係を規定する材料定数Ｃ＊は、図２（ｃ）に示すように、下記式（３）で表される。すなわち、材料定数Ｃ＊を持つ均一な材質の仮想複合材料が想定される。
【数３】
2010250796.tif 000005
【００３５】
ここで、ひずみに注目し、ウレタンＵに発生するひずみεUと複合材料全体に発生するひずみεTとの間の関係、及びフィラーＦに発生するひずみεFと複合材料全体に発生するひずみεTとの間の関係を、下記式（４），（５）で表す。
【００３６】
【数４】
2010250796.tif 000006
【数５】
2010250796.tif 000007
【００３７】
式（４）は、ウレタンＵ内の平均ひずみであるεUがεTの関数として表され、εUが比例定数ＡUによりεTと関係付けられることを表している。式（５）は、フィラーＦ内の平均ひずみであるεFがεTの関数として表され、εFが比例定数ＡFによりεTと関係付けられることを表している。
このような比例定数ＡF，ＡUは、エポキシＰ中に、フィラーＦが既知の体積分率で内側球状粒子として分散し、この内側球状粒子の周囲にウレタンＵを被覆層として所定の厚さで被覆して構成された外側殻状層を有する仮想複合材料を複合材料として定め、さらに、外側殻状層を囲む周囲の材料定数を、算出すべき仮想複合材料の有効材料定数Ｃ＊として定めることにより求めることができる。すなわち、比例定数ＡU，ＡFを、ウレタンＵの材料定数と、フィラーＦの材料定数と、均一な材質の仮想複合材料の有効材料定数Ｃ＊とを用いて表すことができる。ウレタンＵの材料定数をＣUとすると、ＡU＝ＡU（ＣU，ＣF，Ｃ＊）と表され、ＡF＝ＡF（ＣU，ＣF，Ｃ＊）と表される。ＡU＝ＡU（ＣU，ＣF，Ｃ＊）及びＡF＝ＡF（ＣU，ＣF，Ｃ＊）の求め方は後述する。ＡU＝ＡU（ＣU，ＣF，Ｃ＊）及びＡF＝ＡF（ＣU，ＣF，Ｃ＊）は、ＣU，ＣF，Ｃ＊の複雑な式で表された非線形な式となっている。ここで、ＣU，ＣF，Ｃ＊は、例えば、ウレタンＵ及びフィラーＦの体積弾性率ＫUとラメ定数μU、 体積弾性率ＫFとラメ定数μF、及び仮想複合材料の体積弾性率Ｋ＊及びラメ定数μ＊を表す。このような比例定数を用いて、式（１）〜（４）を整理すると、下記式（６）が導かれる。ここにＣPは、エポキシＰの材料定数である。
【００３８】
【数６】
2010250796.tif 000008
【００３９】
式（６）は、左辺の材料定数Ｃ＊を求める式であり、右辺にあるＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）は、ＣU，ＣF，Ｃ＊に関して非線形な式となっている。このため、式（６）は、材料定数Ｃ＊に関して、再帰型非線形方程式となっている。
【００４０】
このような比例定数ＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）は、図３に示すように、以下の構成の外側殻状層が配置されたモデルを仮想し、この外側殻状層を囲む周囲の材料定数を、求めるべき仮想複合材料の有効材料定数Ｃ＊として定めることにより得られる。外側殻状層は、材料定数Ｃ＊を持つ均一な材質の仮想複合材料中にあるフィラーＦの内側球状粒子の周りに、一定の厚さで被覆層としてウレタンＵが覆うように構成されている。
すなわち、ＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）は、体積力が作用せず、無限遠で印加されるひずみε∞ij（ｉ，ｊは、１〜３の自然数）が存在する条件下、公知のナビエの方程式を解くことにより得られる。なお、無限遠で印加されるひずみε∞ijは下記式（７）のように、第１項の静水圧成分と第２項のせん断成分とに分解できる。このため、ＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）についても静水圧成分と第２項のせん断成分に分けて算出することができる。
【００４１】
【数７】
2010250796.tif 000009
【００４２】
以下、ＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）の第１項の静水圧成分と第２項のせん断成分の算出について説明する。
【００４３】
（静水圧成分に基づくＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）の算出）
まず、一般に原点で対称な単一球状粒子が母相にある場合を考える。このとき、式（７）中の第１項の静水圧成分に関するナビエの方程式を満足する変位ｕi（ｉは１〜３の自然数）は、下記式（８）のように表される。これは、式中の左辺と右辺のテンソルの階数を一致させる必要性から一意的に定まる。
ここで、球状粒子における原点の変位ｕiが有限の値を持ち、ひずみ成分が無限遠で式（７）中の第１項のひずみになる条件と、球体粒子の表面における変位ｕiの連続性及び表面張力の連続性の条件とを与えることにより、球体粒子内部の変位ｕiと球体粒子外部の被覆層での変位ｕiを求めることができる。具体的には、球体粒子内部（フィラーＦ）の変位ｕiは下記式（９）で、球体粒子外部（ウレタンＵ）の被覆層での変位ｕiは下記式（１０）で表される。ここで、ｕiFは、球状粒子内部であるフィラーF内の変位であり、ｕiUは、球状粒子外部であるウレタンＵ内の変位である。ａ，ｂは、フィラーＦ及びウレタンＵの半径を表す。
【００４４】
【数８】
2010250796.tif 000010
【数９】
2010250796.tif 000011
【数１０】
2010250796.tif 000012
【００４５】
すなわち、母相に１つの球状粒子がある空間内の各方向の変位ｕiは、球状粒子の内部（フィラーＦ）において、球状粒子の中心を基準として定められる位置ｘi（i＝１、２または３）に比例し、球状粒子の外部（ウレタンＵ）の被覆層において、球状粒子の中心を基準として定められる位置ｘi（i＝１、２または３）に比例する項と、球状粒子の中心からの距離の３乗に反比例し、位置ｘiに比例する項とを有するように定められる。
そのとき、上記式（９）、（１０）で表される球状粒子内部及び被覆層における変位から、ひずみを求めることにより、比例定数ＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）の静水圧成分を求めることができる。下記式（１１）は、体積弾性率Ｋに関する比例定数ＡF（ＣU，ＣF，Ｃ＊）の静水圧成分の式を例示している。同様に、下記式（１２）は、体積弾性率Ｋに関する比例定数ＡU（ＣU，ＣF，Ｃ＊）の静水圧成分の式を例示している。ここで、ｔはウレタンとフィラーの半径の比、すなわちｔ＝ｂ／ａを表し、ＫU，ＫF，Ｋ*は、ウレタンＵの体積弾性率，フィラーＦの体積弾性率及び算出すべき体積弾性率を表し、μU，μF，μ*は、ウレタンＵのラメの定数，フィラーＦのラメの定数及び算出すべきラメの定数を表す。
【００４６】
【数１１】
2010250796.tif 000013
【数１２】
2010250796.tif 000014
【００４７】
（せん断成分に基づくＡ（ＣB，Ｃ＊）の算出）
式（７）中の第２項のせん断成分に関するナビエの方程式を満足する変位ｕi（ｉは１〜３の自然数）は、母相に１つの球状粒子がある空間上の変位を用いて定めることができる。具体的には、母相に１つの球状粒子がある空間上の変位は、この球状粒子の内部及び外部を問わず、球状粒子の中心からの距離の０乗、距離の２乗、距離の−３乗及び距離の５乗のそれぞれに比例する項の加算式で定められる。具体的には、変位ｕiを、下記式(１３)のように、式中の無限遠ひずみ成分に比例するように定める。この式は、テンソルの階数をナビエの方程式上で対応させる必要性から一意的に定められたものである。ここで、ｘｊ，ｘkは、３次元座標におけるｘiと異なる座標成分を表す。また、ｆ（ｒ）、ｇ（ｒ）は球状粒子の中心からの距離ｒのみによる関数を表す。
【００４８】
【数１３】
2010250796.tif 000015
【００４９】
このとき、ナビエの方程式を満足するｆ（ｒ）、ｇ（ｒ）をｆ（ｒ）∝ｒｌ、ｇ（ｒ）∝ｒｌ（ｌは整数）としてｒの次数ｌを次元解析により求めることにより、ｒの次数ｌは０，２，−３，５となる。したがって、ナビエの方程式を満足する変位ｕiは、球状粒子の内部及び外部を問わず、球状粒子の中心からの距離の０乗、距離の２乗、距離の−３乗及び距離の５乗のそれぞれに比例する項の加算式により表される。例えば、次数ｌ＝２のときの式は、下記式（１４）に示すように表される。勿論、式（１４）中のλ、μは、球状粒子の内部では、球状粒子の内部の材料定数が、球状粒子の外部では、球状粒子の外部の材料定数が用いられる。
【００５０】
【数１４】
2010250796.tif 000016
【００５１】
すなわち、変位ｕiは、球状粒子の内部及び外部において、球状粒子の中心からの距離の０乗、距離の２乗、距離の−３乗及び距離の５乗のそれぞれに比例する項の加算式として表すことができる。静水圧成分の導出と同様に、母相中に１つの球状粒子がある空間内の変位ｕiについて、フィラーＦで形成される内側球状粒子と、この内側球状粒子の周囲にウレタンＵを被覆層として所定の厚さｔで被覆して構成された外側殻状層とにそれぞれ適用して、内側球状粒子と、内側球状粒子と同心球殻を成す外側殻状層とで作られる各位置での応力場、ひずみ場及び変位ｕiを求めることができる。具体的には、図３に示すフィラーＦで形成される内側球状粒子の境界上での変位及び表面張力の連続性の条件と、外側殻状層（ウレタンＵの被覆層）の境界上での変位及び表面張力の連続性の条件と、を用いて、変位ｕiを一意的に求めることができる。さらに、この変位ｕiからひずみを求めることにより、比例定数ＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）のせん断成分を求めることができる。求めた比例定数ＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）のせん断成分の式は、極めて複雑であるため省略するが、少なくともＣU，ＣF，Ｃ＊に関して非線形となっている。
【００５２】
以上より、比例定数ＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）が定められ、この比例定数ＡU（ＣU，ＣF，Ｃ＊），ＡF（ＣU，ＣF，Ｃ＊）が、材料定数Ｃ＊を求める式（６）に代入されて、再帰型非線形方程式が導かれる。この再帰型非線形方程式は、母相中に、フィラーＦが既知の体積分率で内側球状粒子として分散し、この内側球状粒子の周囲にウレタンＵが被覆層として所定の厚さ（ｂ−ａ）で被覆して構成された外側殻状層を有する仮想複合材料を定め、さらに、仮想複合材料中の外側殻状層を囲む周囲の材料定数を、算出すべき仮想複合材料の有効材料定数として定めることにより、実現したものである。
この仮想複合材料を利用した本発明の算出方法による複合材料の有効材料定数の計算結果は、後述するように、従来の長時間の計算を要する有限要素モデルを用いた場合に近く、精度が高いものとなっている。
【００５３】
このような再帰型非線形方程式を用いて複合材料の有効材料定数Ｃ＊を算出する方法について説明する。
【００５４】
図４は、上記材料定数Ｃ＊を算出するフローを説明する図である。
まず、条件設定モジュール２０にて、複合材料の有効材料定数Ｃ＊を算出する第１の処理が設定される。この設定は、入力装置３２を介してオペレータから指示入力されて行われる。ここで、材料定数Ｃ＊は、算出すべき未知数である体積弾性率Ｋefとラメ定数μefを代表して表している。以降、ＫP，ＫU，ＫFは、それぞれエポキシＰ、ウレタンＵ、フィラーＦの体積弾性率であり、μP，μU，μFは、それぞれエポキシＰ、ウレタンＵ、フィラーＦのラメ定数であり、ｖU，ｖFは、それぞれウレタンＵ、フィラーＦの体積分率である、とする。
【００５５】
まず、エポキシＰ，ウレタンＵ及びフィラーＦの材料定数ＣP，ＣU，ＣFと、ウレタンＵ及びフィラーＦの体積分率ｖu，ｖF，の値が、記憶装置３０中のデータベースから呼び出されて取得される（ステップＳ１０）。材料定数と体積分率の値は、非線形方程式設定モジュール２２に送られて、上述の式（６）の各係数の値が定められ、材料定数Ｃ＊のみの再帰型非線形方程式に修正される。すなわち、方程式が用意される。
【００５６】
次に、方程式解法モジュール２４にて、材料定数Ｃ＊の未知数であるＫef,μefの初期値が設定される（ステップＳ２０）。初期値を設定する理由は、式（６）が非線形であるため、解が収束するまで繰り返すためである。初期値はどのような値であってもよいが、例えば、母相であるエポキシＰの材料定数ＣPとウレタンＵの材料定数ＣUとフィラーＦの材料定数ＣFの、体積分率による加重平均値を選択するとよい。この初期値を、Ｋef,＝Ｋ1，μef＝μ1と定める。さらに、この初期値をそれぞれ、Ｋn，μn（ｎ＝１）と表す。
次に、方程式解法モジュール２４にて、上述の式（６）を修正した下記式（１５）を定義する。具体的には、式（６）の左辺を右辺に移動し、その時の右辺をＦ（Ｃ＊)として定義する。
【００５７】
【数１５】
2010250796.tif 000017
【００５８】
さらに、下記式（１６）にしたがって、Ｆ（Ｃ＊)を用いてＣ(n+1)が定義される。ここで、式（１６）は、Ｆ（Ｃ＊)＝０となる収束解を求めるためにニュートン・ラフソン法を適用した式である。ここで、Ｃ(n)，Ｃ(n+1)は体積弾性率Ｋn，Ｋn+1及びラメの定数μn，μn+1を代表して表す符号である。式（１６）のＣ(n)に、体積弾性率Ｋn，ラメの定数μnが与えられ、式（１６）に従って体積弾性率Ｋn+1及びラメの定数μn+1が算出され取得される（ステップＳ３０）。ここで、式（１６）中のＦ´（Ｃ(n)）はヤコビのテンソルであり、Ｆ´（Ｃ(n)）-1はＦ´（Ｃ(n)）の逆テンソルを表す。すなわち、Ｆ´（Ｃ(n)）-1・Ｆ´（Ｃ(n)）＝Ｉ（Ｉは単位テンソル）である。
【００５９】
【数１６】
2010250796.tif 000018
【００６０】
算出されたＫn+1，μn+1と、この算出に用いられたＫn，μnが収束解判定モジュール２６に送られて、それぞれの差分Ｋn+1−Ｋn，μn+1−μnの絶対値が求められる。この絶対値が予め設定された閾値ε1，ε2と比較される（ステップＳ４０）。比較の結果、差分Ｋn+1−Ｋn，μn+1−μnの絶対値がいずれも閾値ε1，ε2より小さい条件を満たす場合（Ｙｅｓの場合）、求めるべき仮想複合材料の有効材料定数Ｃ＊であるＫef,μefは、それぞれＫn+1，μn+1に決定される（ステップＳ５０）。
一方、比較の結果、上記条件を満たさない場合（Ｎｏの場合）、ｎ＋１をｎとし（ステップＳ６０）、さらに、式（１６）を用いてＫn+1，μn+1を算出するためにステップＳ３０に戻る。こうして、ステップＳ４０における条件を満足するまで、ステップＳ３０，Ｓ４０，Ｓ６０が繰り返される。
【００６１】
ステップＳ５０において求められた仮想複合材料の有効材料定数Ｃ＊は、エポキシＰに、フィラーＦとウレタンＵを含む複合材料の有効材料定数Ｃ＊として定められる。
【００６２】
次に、このように定められた複合材料の有効材料定数Ｃ＊、具体的には体積弾性率Ｋ*及びラメの定数μ＊を用いてヤング率Ｅ*が算出される。ヤング率Ｅ*は、下記式（１７）に従って算出され、プリンタやディスプレイ等の出力装置３４に出力される。
【００６３】
【数１７】
2010250796.tif 000019
【００６４】
算出されたヤング率Ｅ*は、後述するように、有限要素モデルを用いて数値計算をおこなって得られるヤング率Ｅとよく対応した値となっている。このため、本発明の複合材料の材料定数の算出方法は有効であることがわかる。また、本発明の複合材料の材料定数の算出方法は、式（６）に示される式を用いて解析的に算出できるので、有限要素モデルを用いた計算に比べて短時間に算出結果を得ることができ、効率的である。
【００６５】
本発明の複合材料の材料定数の算出方法で用いる式（６）は、解析式であるので、例えば、所定の材料の体積分率を解くべき未知数とすることもできる。このとき、複合材料の有効材料定数は、予め実験等により求めておく必要がある。このような材料の体積分率の算出は、上述した第２の処理に対応する。この第２の処理の具体的な方法を以下説明する。
【００６６】
図５は、第２の処理のフローを示す図である。エポキシＰに、フィラーＦとウレタンＵを含む複合材料において、フィラーＦの体積分率ｖFを算出する場合を想定する。
【００６７】
まず、条件設定モジュール２０において、エポキシＰ、ウレタンＵ、フィラーＦの材料定数、具体的には、ＫP，μP，ＫU，μU，ＫF，μF及びウレタンＵの体積分率ｖUが取得される（ステップＳ１００）。これらの材料定数等は、記憶装置３０中のデータベースから呼び出されて取得される。あるいは、入力装置３２から入力される。
次に、複合材料の力学的実験による計測結果より、複合材料全体の材料定数である有効材料定数、具体的にはＫ＊，μ＊が特定され取得される（ステップＳ１１０）。計測結果は、予め記憶装置３０に記憶されており、条件設定モジュール２０に計測結果が呼び出されて、複合材料の有効材料定数が特定される。さらに、非線形方程式設定モジュール２２では、取得された材料定数や体積分率の各値が式（６）の各係数に代入されて、未知数が体積分率ｖFの方程式が用意される。この場合、関数ｆ（ｘ）における体積分率ｘは、フィラーＦの体積分率であり、ｆ（ｘ）で用いる式（６）は、エポキシＰにフィラーＦが分散し、このフィラーＦの周りに一定の厚さでウレタンＵが被覆された構成の仮想複合材料の有効材料定数を求める式である。
【００６８】
この後、フィラーＦの体積分率ｖFの上限値ｘ1及び下限値ｘ2が設定される（ステップＳ１２０）。上限値及び下限値は、後述する体積分率ｖFの算出を、２分法を用いて行うために設定される。この設定は、入力装置３２を通して、オペレータの指示入力に基づいて設定されてもよいし、予め設定されたデフォルト値を用いてもよい。
【００６９】
次に、方程式解法モジュール２４において、上限値ｘ1及び下限値ｘ2が用いられて、ｆ（ｘ1）・ｆ（ｘ2）の値が負か否かが判定される（ステップＳ１３０）。ｆ（ｘ1）・ｆ（ｘ2）の値が正の場合、上限値ｘ1及び下限値ｘ2が変更される（ステップＳ１４０）。変更の方法は、特に制限されないが、ｘ1及びｘ2は体積分率であるので、上限値ｘ1はより大きく、下限値ｘ2はより小さくなるように設定されるとよい。
【００７０】
ステップＳ１３０における判定で肯定される場合（Ｙｅｓの場合）、公知の２分法に基づいて、以下の処理が行われる。
すなわち、まず、ｘ3＝（ｘ1＋ｘ2）／２が求められ（ステップＳ１５０）、さらに、ｆ（ｘ1）・ｆ（ｘ3）の値が負か否かが判定される（ステップＳ１６０）。判定の結果、ｆ（ｘ1）・ｆ（ｘ3）が負である場合（Ｙｅｓの場合）、下限値ｘ2がｘ3に変更され（ステップＳ１７０）、ｆ（ｘ1）・ｆ（ｘ3）が負でない場合（Ｎｏの場合）、上限値ｘ1がｘ3に変更される（ステップＳ１８０）。
この後、上限値ｘ1と下限値ｘ2との差分の絶対値が予め設定された閾値ε5より小さいか否かが判定される（ステップＳ１９０）。絶対値が閾値ε5より小さい場合（Ｙｅｓの場合）、上限値ｘ1または下限値ｘ2が算出すべき体積分率ｖFとして定める（ステップＳ２００）。ステップＳ１９０において、絶対値が閾値ε5より小さくない場合（Ｎｏの場合）、ステップＳ１５０に戻る。こうして、ステップＳ１９０の判定で、絶対値が閾値ε5より小さくなるまで、ステップＳ１５０、Ｓ１６０、Ｓ１７０、Ｓ１８０を繰り返す。
【００７１】
以上の方法により、フィラーＦの体積分率ｖFは算出される。算出したフィラーＦの体積分率ｖFは、出力装置３４に出力される。
なお、算出すべき体積分率は、フィラーＦを対象とする他、ウレタンＵを対象とすることもできる。また、算出すべき体積分率は、複数の種類の材料を対象としてもよい。
複合材料の材料定数の算出方法及び体積分率の算出方法では、母相に２つの異なる材料を含む複合材料を用いて説明したが、母相に含まれる材料は３つ以上であってもよい。母相は、エポキシに限られず、金属やセラミック等の無機材料であってもよい。また、母相に含まれる材料は、母相を強化する強化材料に限られず、上述したようにウレタン等の軟質材料であってもよい。
また、複合材料の母相に含まれる材料は、必ずしも粒子として分散しているものである必要は無く、母相に相溶するものを仮想複合材料として定めてもよい。
【００７２】
また、複合材料の材料定数の算出方法及び体積分率の算出方法は、コンピュータ上でプログラムを実行することにより実現することもできる。このプログラムはコンピュータで読み取り可能な記録メディアに記録されたものである。この記録メディアは、通信回線を利用してダウンロードされたプログラムを記録したものも含まれる。
【００７３】
次に、上述の複合材料の有効材料定数の算出方法の有効性を説明する。
上記実施形態では、エポキシＰ、フィラーＦ及びウレタンＵからなる複合材料を用いて説明したが、以下の例ではゴムＲ、シリカＳ、及びバウンドラバーＢからなるゴム材料を用いた。
（１）ゴム：１００ｇ
ＳＢＲ（スチレン−ブタジエンゴム）：    ７５ｇ
ＢＲ（ブタジエンゴム）：                ２５ｇ
（２）シリカ：３０ｇ
（３）その他：
カップリング剤                          ２．４ｇ
アロマオイル：                        ３４．１２ｇ
亜鉛：                                  ２．５ｇ
ステアリン酸：                        １．５ｇ
硫黄：                                  １．４７ｇ
加硫促進剤：                              ３．７ｇ
パラフィンワックス：                      １ｇ
その他薬品：                              ４ｇ
作製した上記ゴム材料では、シリカ粒子の周囲を取り巻くようにバウンドラバー層が形成されている。したがって、ゴム材料は、ゴムを母相とし、この母相中にシリカが球状粒子として分散し、このシリカの周りを被覆層としてバウンドラバー層が覆う構成となっている。
【００７４】
上記ゴム材料におけるヤング率及びポアソン比と、体積分率（％）は以下の通りである。なお、ヤング率及びポアソン比から、公知の式により、ラメ定数及び体積弾性率に変換し、変換したラメ定数及び体積弾性率を用いた。
【表１】
2010250796.tif 000020
【００７５】
本発明の複合材料の材料定数の算出方法を用いたゴム材料全体の算出結果のヤング率は１．４６３(ＧＰａ)であった。一方、有限要素モデルを用いた複合材料全体の解析により得られたヤング率は１．４３１（ＧＰａ）であった。
これより、本発明の方法により得られるヤング率は、有限要素モデルを用いて得られるヤング率と極めて近い値を示し、本発明の方法が有効であることがわかる。
なお、有限要素モデルを用いた材料定数の算出は、図６に示すような、ゴムＲを母相とし、このゴムＲ中に球状のシリカＳが体積分率８．８６％含まれ、バウンドラバーＢがシリカＳの周りに体積分率８．８６％含まれて構成される複合材料を想定した有限要素モデルを作成して、複合材料の有効材料定数を算出した。図中下端部を拘束し上端部に引っ張り変位を与えて応力−ひずみの初期傾きを調べてヤング率とした。なお、作成した有限要素モデルは、要素数２５６００個、節点数２５７６１個の２次元平面応力要素モデルである。解析は、非線形有限要素法ソフトウェアＡＢＡＱＵＳ（商品名）を用いた。
【００７６】
以上、本発明の複合材料の材料定数の算出方法、複合材料中の材料の体積分率の算出方法及び記録メディアについて詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。
【符号の説明】
【００７７】
１０  計算装置
１２  ＣＰＵ
１４  ＲＯＭ
１６  ＲＡＭ
１８  入出力ポート
２０  条件設定モジュール
２２  非線形方程式設定モジュール
２４  方程式解法モジュール
２６  収束解判定モジュール
２８  結果処理モジュール
３０  記憶装置
３２  入力装置
３４  出力装置
【図１】
2010250796.tif 000021
【図２】
2010250796.tif 000022
【図３】
2010250796.tif 000023
【図４】
2010250796.tif 000024
【図５】
2010250796.tif 000025
【図６】
2010250796.tif 000026