(57)【要約】【課題】熱硬化性樹脂中へ粒子を配合する工程での粒子の攪拌挙動を迅速かつ高精度に予測する。【解決手段】熱硬化樹脂への粒子分散の設計支援装置がモデル作成部１２と流動解析部１３から構成され、モデル作成部１２は各版容器と液状樹脂部の形状を定義し、流動解析部１３内のクラスター内粒子数解析部１４では液状熱硬化性樹脂中に投入されたクラスター内粒子数の時間変化率を計算し、微小時間経過後の前記クラスター内粒子数の増分あるいは減分を近似的に求め、新しい時刻におけるクラスター内粒子数を逐次求めていく。流動解析部１３内の熱硬化性樹脂粘度解析部１５では粘度算出式を用いてクラスター内粒子数解析部１４中の粘度を計算する。また、流動解析部１３内の温度、せん断速度解析部１６では樹脂の温度、せん断速度分布を算出してクラスター内粒子数解析部１４に代入し、攪拌層内のクラスター内粒子数の分布と時間変化を逐次計算する。【選択図】図１
2010184365.tif 000002
【特許請求の範囲】
【請求項１】
モデル作成部と流動解析部から構成される熱硬化性樹脂への粒子分散の設計支援装置であって、前記モデル作成部は攪拌容器と容器内の液状樹脂部の形状を定義し、前記流動解析部は、クラスター内粒子数解析部と熱硬化性樹脂粘度解析部と温度、せん断速度解析部とを備え、前記クラスター内粒子数解析部は、前記熱硬化性樹脂中に投入された粒子の凝集体であるクラスターを定義し、前記クラスター内粒子数解析部は、クラスター内粒子数の時間変化率を粒子の凝集項とせん断破壊項の差によって計算し、差分式により微小時間経過後の前記クラスター内粒子数の増分あるいは減分を近似的に求め、新しい時刻における前記クラスター内粒子数を逐次求めていくモデルを用いてクラスター内粒子数を解析し、前記熱硬化性樹脂粘度解析部は、前記凝集項ならびに前記せん断破壊項には粘度項が含まれ、前記粘度は反応の進行により値が変化する状態を表せる粘度算出式を用いて粘度を算出し、温度、せん断速度解析部は、前記流体と熱の移動を記述する質量・運動量・エネルギーの保存方程式と組み合わせて、有限差分法、あるいは有限要素法を用いて数値解析し、前記攪拌槽モデル部内の樹脂の温度、粘度、せん断速度分布を算出して前記凝集項とせん断破壊項に代入して前記攪拌層モデル部内のクラスター内粒子数の分布と時間変化を逐次計算することを特徴とする熱硬化性樹脂への粒子分散の設計支援装置。
【請求項２】
請求項１に記載の設計支援装置であって、前記粘度算出式は温度と時間の関数となる式を用いることを特徴とする熱硬化性樹脂への粒子分散の設計支援装置。
【請求項３】
請求項１に記載の設計支援装置であって、前記粘度算出式は温度と反応率の関数となる式を用いるとともに反応率を算出するための反応速度式を併せて用いることを特徴とする熱硬化性樹脂への粒子分散の設計支援装置。
【請求項４】
コンピュータにより熱硬化性樹脂への粒子分散の設計を支援する方法であって、モデル作成ステップと流動解析ステップからなり、前記モデル作成ステップにおいて攪拌容器と容器内の液状樹脂部の形状を定義するとともに、前記流動解析ステップでは樹脂中に投入された粒子の凝集体であるクラスターを定義し、前記クラスター内粒子数の時間変化率を粒子の凝集項とせん断破壊項の差によって計算し、差分式により微小時間経過後の前記クラスター内粒子数の増分あるいは減分を近似的に求め、新しい時刻における前記クラスター内粒子数を逐次求めていくモデルを用い、前記凝集項ならびに前記せん断破壊項には粘度項が含まれ、前記粘度は反応の進行により値が変化する状態を表せる粘度算出式を用い、流体と熱の移動を記述する質量・運動量・エネルギーの保存方程式と組み合わせて、有限差分法、あるいは有限要素法を用いて数値解析し、前記攪拌槽モデル部内の樹脂の温度、粘度、せん断速度分布を算出して前記凝集項とせん断破壊項に代入して前記攪拌層モデル部内のクラスター内粒子数の分布と時間変化を逐次計算することを特徴とする熱硬化性樹脂への粒子分散の設計を支援する方法。
【請求項５】
コンピュータにより熱硬化性樹脂への粒子分散の設計を支援するプログラムであって、モデル作成ステップと流動解析ステップからなり、前記モデル作成ステップにおいて攪拌容器と容器内の液状樹脂部の形状を定義するとともに、前記流動解析ステップでは樹脂中に投入された粒子の凝集体であるクラスターを定義し、前記クラスター内粒子数の時間変化率を粒子の凝集項とせん断破壊項の差によって計算し、差分式により微小時間経過後の前記クラスター内粒子数の増分あるいは減分を近似的に求め、新しい時刻における前記クラスター内粒子数を逐次求めていくモデルを用い、前記凝集項ならびに前記せん断破壊項には粘度項が含まれ、前記粘度は反応の進行により値が変化する状態を表せる粘度算出式を用い、流体と熱の移動を記述する質量・運動量・エネルギーの保存方程式と組み合わせて、有限差分法、あるいは有限要素法を用いて数値解析し、前記攪拌槽モデル部内の樹脂の温度、粘度、せん断速度分布を算出して前記凝集項とせん断破壊項に代入して前記攪拌層モデル部内のクラスター内粒子数の分布と時間変化を逐次計算することを特徴とする熱硬化性樹脂への粒子分散の設計を支援するプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、熱硬化性樹脂への粒子分散の設計支援技術に関する。
【背景技術】
【０００２】
熱硬化性樹脂はその硬化物が接着性、機械的強度、電気絶縁性、耐薬品性などに優れた性質を有しているため、電気機器の構造絶縁材料や電子、機構部品の封止材料として広く用いられている。この中でモータや発電機、変圧器、開閉器のような高電圧機器においては電気絶縁性、機械的特性、信頼性確保のためにコイルなどの電流が流れる箇所とその周りの筐体との間に樹脂を注入して固める方法が用いられる。
【０００３】
近年、これらの機器は高出力化のために発熱量が増え、小型化のために使用樹脂量が減少し、樹脂による強度確保が困難な状況になりつつある。このため高強度樹脂の開発と実用化が急務になっている。樹脂の強度を飛躍的に向上させる手段として有力なのは、樹脂にナノサイズの粒子を配合させ、分子レベルの補強効果を付与するナノコンポジット材料である。この配合工程では、攪拌槽の中の液状の熱硬化性樹脂にナノサイズの粒子を投入し、羽根の回転や遠心力を与えて粒子を樹脂中に分散させる手法が用いられる。
【０００４】
しかし、粒子が非常に小さく、粒子数が非常に多いため、粒子は極めて凝集しやすく、均一に分散させるための条件設定が非常に難しい状況になっている。攪拌装置仕様や条件を合理的に決定するためには攪拌挙動のシミュレーション技術の構築が重要となってきた。特に熱硬化性樹脂は反応の進行により攪拌中に粘度が複雑に変化するため、適切な粒子の挙動モデルと連立させて、数値解析により温度、粘度、流速などの物理量の変化を計算し、各種条件と粒子の分散状態の関連を事前に予測することが必要となる。
【０００５】
本発明は熱硬化性樹脂への粒子分散の設計支援行う装置、方法、およびコンピュータプログラムを対象とする。
【０００６】
特許文献１は反応速度モデルを基にして熱硬化性樹脂の流動挙動と硬化後の残留ひずみまでを一貫して解析する装置、方法である。特許文献２もやはり反応速度モデルを基にして熱硬化性樹脂の流動挙動を解析する装置、方法である。特許文献３は所定の空間内で運動する粉体の挙動を解析する装置、方法に係り、現像、転写などのプロセスを含む電子写真プロセスにおける現像剤の挙動を解析する装置、方法とコンピュータプログラムである。特許文献４は容器内における粒子の微小時間経過毎の位置を求めるシミュレーション装置、方法に関するものであり、特に電子写真プロセスでの粉体の混合攪拌時の静電気容量をシミュレートする装置と方法である。特許文献５も特許文献４と同じである。特許文献６は電子写真技術を用いた画像形成装置に関連し、容器内における粒子位置を微小時間毎に求めるシミュレーション装置と記録媒体である。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００６−２０５７４０号公報
【特許文献２】特開平１１−２３２２５０号公報
【特許文献３】特開２００６−３３０１５６号公報
【特許文献４】特開平１０−２６０１５９号公報
【特許文献５】特開平１０−２６０１６０号公報
【特許文献６】特開２００３−２２３０４９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
特許文献１と２は熱硬化性樹脂の解析モデルとしてはかなり厳密であるが、樹脂は均質体としての取り扱いとなる。すなわち、樹脂中に充填材としての粒子が含まれていてもその粒子は時間や場所によらず常に均一の分布を保っているという前提がある。したがって、最初は凝集している粒子が樹脂中へ分散していく時間変化は取り扱えないという問題がある。
【０００９】
特許文献３は粉体が流体中の拡散していく挙動をシミュレーションするが、熱硬化性樹脂特有の粘度変化は取り扱っていない。特許文献４と特許文献５は静電気を持つ粉体の空気中での攪拌時の挙動を解析するもので熱硬化性樹脂中への粒子の分散とは物理モデルが異なっている。
【００１０】
特許文献６は攪拌時の粒子の接触など粒子間の相互作用を考慮して解析する手法であるが、熱硬化性樹脂中への粒子の分散とは物理モデルが異なっている。以上述べたように、熱硬化性樹脂中への粒子の分散を予測する手法がこれまでになく、粒子配合工程では試行錯誤を繰り返すのが現状となっている。
【００１１】
本発明は上記事情に鑑みてなされたものであり、本発明の目的は熱硬化性樹脂への粒子の分散挙動を事前に予測し、この結果を基に攪拌装置の設計や条件設定を合理的に行うことにある。
【課題を解決するための手段】
【００１２】
上記問題を解決するために、本発明の設計支援装置は、攪拌槽のモデル作成部と流動解析部から構成される。前記モデル作成部において前記攪拌槽の形状や回転状態を定義する。前記流動解析部では液状熱硬化性樹脂中に投入された粒子の凝集体であるクラスターを定義し、前記クラスター内粒子数の時間変化率を粒子の凝集項とせん断破壊項の差によって計算し、差分式により微小時間経過後の前記クラスター内粒子数の増分あるいは減分を近似的に求め、新しい時刻における前記クラスター内粒子数を逐次求めていくモデルを用いる。前記凝集項ならびに前記せん断破壊項には粘度項が含まれ、前記粘度は反応の進行により値が変化する状態を表せる粘度算出式を用いる。
【００１３】
さらに、流体と熱の移動を記述する質量・運動量・エネルギーの保存方程式と組み合わせて、有限差分法、あるいは有限要素法を用いて数値解析し、前記攪拌槽モデル部内の樹脂の温度、粘度、せん断速度分布を算出して前記凝集項とせん断破壊項に代入して前記攪拌層モデル部内のクラスター内粒子数の分布と時間変化を逐次計算する。
【発明の効果】
【００１４】
本発明によれば、クラスター内粒子数の変化は樹脂の熱・流動解析で得られる任意の場所と時間における温度、粘度、せん断速度などの値から計算できる。すなわち、粒子形状や数を幾何学的にモデル化せずに目的の計算が可能になるので、迅速なシミュレーションを行うことができる。また、熱硬化性樹脂特有の反応の進行による粘度変化を考慮するので高精度のシミュレーションが可能となる。
【図面の簡単な説明】
【００１５】
【図１】図１は本発明の実施例１の設計支援装置の概略構成図である。
【図２】図２は図１に示す設計支援装置のハードウエア構成例を示す図である。
【図３】図３は本発明の実施例１の３次元流動解析処理を説明するためのフロー図である。
【図４】図４は本発明の実施例１の粘度式の等温粘度変化を示す図である。
【図５】図５は本発明の実施例１の非等温粘度変化を取り扱う説明図である。
【図６】図６は本発明の実施例１の攪拌槽の温度と攪拌限界粘度到達時間の関係の計算例である。
【図７】図７は攪拌温度７５℃における時間と粘度比の関係の計算例である。
【図８】図８はせん断速度をパラメータとしたときの時間とクラスター内粒子数の関係の計算例である。
【図９】図９は粒子径をパラメータとしたときの時間とクラスター内粒子数の関係の計算例である。
【図１０】図１０は熱硬化性樹脂の反応率α、粘度ηの変化を説明する図である。
【図１１】図１１は本発明の実施例２の設計支援装置の概略構成図である。
【図１２】図１２は本発明の実施例２の３次元流動解析処理を説明するためのフロー図である。
【発明を実施するための形態】
【００１６】
以下、本発明の実施形態について図面を参照して説明する。
【実施例１】
【００１７】
図１は本発明の実施例１の設計支援装置の概略構成図である。図示するように、本実施例の設計支援装置はＧＵＩ（Graphical User Interface）部１１とモデル作成部１２、流動解析部１３を有する。流動解析部は１３はクラスター内粒子数解析部１４と熱硬化性樹脂粘度解析部１５ならびに温度、せん断速度解析部１６に分かれ、お互いに解析データの授受を行いながら、それぞれ、モデル作成部１２で作成したモデルに対応した解析を行う。
【００１８】
ＧＵＩ部１１は、画面表示およびキーボード、マウス等の入力装置を介して、ユーザより各種指示や情報などの入力を受け付けたり、３次元流動解析の結果を表示したりする。また、モデル作成部１２は、ＧＵＩ部１１を介して受け付けたユーザの指示に従い、設計支援の対象とする攪拌容器と容器内の液状樹脂部の形状データ（モデルデータ）を作成する。
【００１９】
流動解析部１３は、解析対象領域で流動中の熱硬化性樹脂への粒子の分散状態の変化を解析する（３次元流動解析）。流動解析部１３は、クラスター内粒子数解析部１４、熱硬化性樹脂粘度解析部１５、及び温度、せん断速度解析部１６から構成される。
【００２０】
クラスター内粒子数解析部１４では樹脂中に投入された粒子の凝集体であるクラスターを定義し、前記クラスター内粒子数の時間変化率を粒子の凝集項とせん断破壊項の差によって計算し、差分式により微小時間経過後の前記クラスター内粒子数の増分あるいは減分を近似的に求め、新しい時刻における前記クラスター内粒子数を逐次求める計算を行う。
【００２１】
熱硬化性樹脂粘度解析部１５は前記粒子の凝集項中の分散媒（粒子を投入する前の溶液）粘度と前記せん断破壊項中のスラリー（粒子を投入後の溶液）粘度を計算する。分散媒粘度の計算では熱硬化性樹脂の反応の進行による粘度変化が計算できる粘度算出式を用いる。スラリー粘度の計算では上記分散媒粘度と粒子含有率、前記クラスター内の粒子数、クラスターの空隙率などの関数となる粘度式を用いる。なお、前記分散媒粘度算出式は温度と時間の関数形、温度と反応率の関数形の２種類があり、後者は反応率算出用の反応速度式も含まれる。粘度計算結果はクラスター内粒子解析部１４と温度、せん断速度解析部１６に送られ、それぞれ値がセットされる。
【００２２】
温度、せん断速度解析部１６は、流体と熱の移動を記述する質量・運動量・エネルギーの保存方程式からなり、熱硬化性樹脂粘度解析部から送られた粘度データは運動量保存方程式中にセットされるとともに、有限差分法、あるいは有限要素法を用いて数値解析し、温度とせん断速度の計算結果は前記クラスター内粒子解析部１４に送られる。これにより、新しい時刻でのクラスター内粒子数が求められる。以上の手順は微小時間経過ごとにモデル作成部１２内の液状樹脂部全域で行われ、解析を続行する。
【００２３】
上記構成の設計支援装置は、例えば図２に示すような、ＣＰＵ２１と、メモリ２２と、ＨＤＤ等の外部記憶装置２３と、ＣＤ−ＲＯＭやＤＶＤ−ＲＯＭ等の可搬性を有する記憶媒体２４から情報を読み出す読取装置２５と、キーボードやマウスなどの入力装置２６と、ＣＲＴやＬＣＤなどの表示装置２７と、インターネットなどのネットワーク通信を行なうための通信装置２８とを備えた一般的なコンピュータシステムにより構成される。
【００２４】
本発明の設計支援装置は、図２のコンピュータシステムにおいて、ＣＰＵ２１がメモリ２２上にロードされた所定のプログラム（モデル作成部１２を実現する３Ｄ−ＣＡＤ、ＣＡＭあるいはＣＡＥプログラム、熱硬化性樹脂流動解析部１３を実現する３次元流動解析プログラム）を実行することで実現できる。これらのプログラムは、読取装置２５を介して記憶媒体２４から、あるいは、通信装置２８を介してインターネットなどからの通信媒体から、メモリ２２に直接ロードしてもよいし、あるいは、一旦、外部記憶装置２３にダウンロードしてから、メモリ２２にロードしてもよい。
【００２５】
図３は本発明の実施例１の設計支援装置による３次元流動解析処理を説明するためのフロー図である。流動解析部１３ではＧＵＩ部１１を介してユーザより３次元攪拌解析に必要な、熱硬化性樹脂の物性値を受け取る（Ｓ２０１）。本実施例では後述する粘度計算式中の係数、比熱、密度、熱伝導率などである。
【００２６】
次に、熱硬化性樹脂流動解析部１３ではＧＵＩ部１１を介してユーザより３次元流動解析のための各種条件（境界条件、解析条件、および初期条件）を受け取る（Ｓ２０２）。各種条件には、攪拌槽の温度、形状データ、回転速度、使用樹脂量、攪拌される粒子の大きさ、粒子の含有率などが含まれる。
【００２７】
モデル作成部１２は、ＧＵＩ部１１を介してユーザより３次元流動解析指示を受け着ける。それから、指定されたモデルデータを、モデル作成部１２から取り込む（Ｓ２０３）。次に、モデル作成部１２は、取り込んだデータが特定する樹脂領域を３次元解析領域に設定する。そして、この３次元解析領域を複数の３次元ソリッド要素に分割する際の条件（分割数や要素サイズなど）を、ＧＵＩ部１１を介してユーザより受け取る（Ｓ２０４）。そして、受け付けた分割条件に従い、３次元解析領域を複数の３次元ソリッド領域に分割する（Ｓ２０５）。
【００２８】
次に、流動解析部１３の中のクラスター内粒子数解析部１４において初期条件から計算できる固定パラメータの値を計算する（Ｓ２０６）。具体的な内容については後述する。次にクラスターない粒子解析部１４で用いる比較用クラスター内粒子数ｎ０とクラスター切断回数Ｎｂをセットする（Ｓ２０７）。この具体的な内容は後述する。
【００２９】
以上の準備が終わるとクラスター内粒子数解析部１４の中のクラスター内粒子数変化率ｄｎ／ｄｔの計算を行う（Ｓ２０８）。この式は下記の形を用いる。
ｄｎ／ｄｔ＝ＢｇＴＮ０／η＋Ｓｇφｎγ−Ｓｈ（ｄ０３ｎ／Ｎｂ）（ｎ／（１−ε）−１）ηｓγ２                                  （式１）
【００３０】
ここで、ｎ：クラスター内粒子数、ｔ：時間、Ｂｇ：ブラウン凝集速度定数とボルツマン定数を含む係数、Ｔ：温度、Ｎ０：解析領域内の粒子数、η：分散媒（粒子を含まない熱硬化性樹脂単体）の粘度、Ｓｇ：せん断凝集速度定数を含む係数、φ：固体体積分率、γ：せん断速度、Ｓｈ：粒子間結合エネルギーを含む係数、ｄ０：粒子直径、Ｎｂ：クラスター切断回数、ε：クラスター内ボイド率、ηｓ：スラリー（粒子を含む熱硬化性樹脂）の粘度  である。
【００３１】
式１の右辺第１項はブラウン運動による粒子の凝集項、第２項はせん断速度が働くとクラスターが変形して隣のクラスターと結合しやすくなることを表すせん断凝集項、第３項はせん断速度が働くとクラスターにせん断エネルギーが働きクラスターが破壊して分散が進行することを表すせん断破壊項である。数１の意味は右辺第１項と第２項の和が右辺第３項よりも大きくなればｄｎ／ｄｔは正となり、クラスター内の粒子数が増加する方向、すなわち、凝集が進むことを示す。反対に負になればクラスター内の粒子数が減少する方向、すなわち、分散が進むことを示す。数１においてクラスター内ボイド率εは次式で計算される。
ε＝εｍａｘ（１−ｎ（−０．４））            （式２）
【００３２】
ここで、εｍａｘ：クラスター内最大空隙率である。クラスターは粒子の凝集体であり、粒子間には空隙部が存在し、ｎ＝１（完全分散）でε＝０、ｎが増加するとεも増加する傾向を持つ。数２はこれを表せる形になっている。クラスター内粒子数ｎは後述する方法により設定する。また、固体体積分率φは次式で計算される。
φ＝φ０／（１−ε（１−φ０）                （式３）
【００３３】
ここで、φ０：凝集なしの状態での固体体積分率である。粒子が凝集したクラスターが存在すると空隙部の存在によりφはφ０よりも大きくなる。数３を用いるとその状態が計算できる。数１の解析領域内の粒子数Ｎ０は樹脂に投入する粒子の体積と粒子径ｄ０から計算する。また、粒子体積と樹脂体積の関係から粒子の占める割合であるφ０が求められる。通常、攪拌開始直前には粒子がクラスターに分かれており、その大きさからクラスター内粒子数ｎを計算する。また、クラスター切断回数Ｎｂは１回の切断により２つに分裂することを繰り返したとしてＮ０とｎの関係から数学的に求めることができる。数１の溶媒粘度ηは後述する熱硬化性樹脂用の粘度計算式から求める。一方、スラリー粘度ηｓは次式で求められる。
ηｓ＝ηｅｘｐ（２．５φ／（１−φ／φｍ））  （式４）
【００３４】
ここで、φｍ：最大個体体積分率である。数４は固体体積分率φが大きくなり、φｍに近づくにつれて粒子を含む熱硬化性樹脂の粘度が急激に上昇する状態を表すことができる。また、数１のせん断速度γは後述する粘度計算式と質量、運動量、エネルギー保存式を用いて求められる。また、せん断発熱などにより温度が変化する場合の温度も求められる。
【００３５】
以上をまとめると、Ｓ２０２の受付条件は数１の中のＢｇ、Ｔ、Ｓｇ、Ｓｈ、ｄ０、数２のεｍａｘ：、式４のφｍとなる。また、Ｓ２０６で計算される値は数１のＮ０、φ、ｎ、Ｎｂ、εとなる。また、Ｓ２０２で初期のｎがｎ０として、Ｎｂはそのままセットされる。
【００３６】
次に、熱硬化性樹脂粘度解析部１５の具体例について説明する。熱硬化性樹脂粘度解析部１５では熱硬化性樹脂用粘度計算式を用いて時間ｔを初期時間にセットする。それから３次元ソリッド要素毎に時間ｔにおける熱硬化性樹脂分散媒の粘度を算出する（Ｓ２０９）。
なお、等温粘度式は数１〜数４で示される。
η＝η０（（１＋ｔ／ｔ０）／（１−ｔ／ｔ０））Ｃ  （式５）
η０＝ａ  ｅｘｐ（ｂ／Ｔ）                        （式６）
ｔ０＝ｄ  ｅｘｐ（ｄ／Ｔ）                        （式７）
Ｃ＝ｆ／Ｔ−ｇ                                    （式８）
【００３７】
ここで、ηは分散媒粘度、ｔは時間、Ｔは温度、η０は初期分散媒粘度、ｔ０はゲル化時間、Ｃは粘度上昇係数、ａ，ｂ，ｄ，ｅ，ｆ，ｇは材料の固有係数である。図４はこの粘度式の等温特性を示したものである。各温度Ｔにおいて分散媒粘度は時間ｔで初期粘度となり、時間の経過とともに反応の進行により分散媒粘度が増大し、ゲル化時間において分散媒粘度は無限大となる。また、温度が高くなるにしたがって、初期分散媒粘度は低くなり、ゲル化時間は短くなる。
【００３８】
一方、非等温状態での粘度変化は以下の手法で予測できる。
式１で、
μ＝（η／η０）（１／Ｃ）                      （式９）
τ＝ｔ／ｔ０                                （式１０）
と置くと次式が得られる。
μ＝（１＋τ）／（１−τ）                  （式１１）
【００３９】
ここで、μは無次元分散媒粘度、τは無次元時間である。μ−τ特性曲線を図５の（ａ）に示す。いま、状態１で時間ｔ１、温度Ｔ１とし、ここからそれぞれ微小量Δｔ、ΔＴ変化し、状態２でｔ２，Ｔ２になったとする。これにより、μ−τ特性曲線上ではτ１からΔτ変化し、μ１がμ２になる。すなわち、時間と温度が同時に変化する現象を一本の曲線上の変化として取り扱うことができる。Δτを微小量とすれば近似的に次式が成立する。
μ２＝μ１＋（ｄμ／ｄτ）μ１Δτ
＝μ１＋２Δτ／（（１−τ１）２）      （式１２）
【００４０】
また、Δτは近似的に次式で求めることができる。
Δτ＝（∂τ／∂ｔ）τ１Δｔ＋（∂τ／∂Ｔ）τ１ΔＴ
＝Δｔ／（ｄ  ｅｘｐ（ｅ／Ｔ１））＋ΔＴ  ｅ  τ１／（Ｔ１２）  （式１３）
【００４１】
式１３を式１２に代入すれば既知の値を用いてμ２が求められる。式１１により、
τ２＝（μ２−１）／（μ２＋１）         （式１４）
となる。
【００４２】
式９から、状態２の粘度は次式で求められる。
η２＝η０（Ｔ２）μ２Ｃ（Ｔ２）              （式１５）
【００４３】
以上の手順をτ＝０から１とみなせる無次元時刻まで繰り返せば、図５−（ｂ）に示す初期状態からゲル化に至るまでの実際の分散媒粘度変化が計算できる。このようにして得られた粘度は逐次クラスター内粒子数解析部１４に送られ、数１の分散媒粘度ηがセットされる。一方、数４からスラリー粘度ηｓが求められ数１のηｓがセットされる。
【００４４】
次に温度、せん断速度解析部１６の具体例について説明する。ここでは質量、運動量、エネルギ保存方程式を用いて、要素毎に、時間ｔにおける温度、速度、せん断速度などを算出する（Ｓ２１０）。
【００４５】
なお、質量保存方程式は、式１６で表される。
（∂ρ／∂ｔ）  ＋ρ（▽・ｖ）＝０                        （式１６）
また、運動量保存方程式は、数１７で示される。
ρ（∂ｖ／∂ｔ  ＋  ｖ・▽ｖ）＝−▽ｐ  ＋▽・τ＋ρｇ    （式１７）
また、エネルギ保存方程式は、数１８で示される。
ρＣｐ（∂Ｔ／∂ｔ  ＋  ｖ・▽Ｔ）＝λ▽２Ｔ  ＋τ：▽ｖ  （式１８）
【００４６】
ここで、ρは密度、ｔは時間、▽はナブラ演算子、ｖは速度ベクトル、ｐは圧力、τは偏差応力テンソル、ｇは重力ベクトル、Ｃｐは定圧比熱、Ｔは温度、λは熱伝導率である。なお、τは粘度ηｓと速度勾配であるせん断速度を用いて計算される。式１６〜式１８は３次元の偏微分方程式であり、厳密解は求められないので、有限差分法あるいは有限要素法などの数値解析手法により、温度、速度、せん断速度などの近似解が求められる。
【００４７】
ここで計算されたせん断速度γの値はクラスター内粒子数解析部１４の中の数１にセットされる。また、温度は次の時刻での数１の温度Ｔ、分散媒粘度η、スラリー粘度ηｓの計算用にセットされる。以上で初期状態における数１の右辺にすべて値が代入される。
【００４８】
次に、クラスター内粒子数解析部１４において数１の計算結果からｄｎ／ｄｔ＜０か？という判断が行われる（Ｓ２１１）。これがＮｏになれば凝集進行あるいは分散が進まない状態となるので、分散解析は終了となる。Ｙｅｓになれば分散が進行しているので次のη／η０＜所定値か？という判定に進む（Ｓ２１２）。
【００４９】
ここで、η０は式６に記載した分散媒初期温度である。温度一定の場合は反応の進行によりηは単調増加、温度上昇を伴う場合は一旦減少した後増加する。粒子の攪拌工程で時間を要すると反応が進行しすぎて粘度上昇し、実用に適さなくなるため、ここでは初期粘度に対する粘度の比であるη／η０を計算し、所定値（例えば１．５程度の値）と比較し、Ｎｏならその時間で解析終了となる。Ｙｅｓの場合はまだ攪拌可能なので時間をタイムステップΔｔ進める（Ｓ２１３）。ここで新しい時刻でのクラスター内粒子数ｎを計算する（Ｓ２１４）。
【００５０】
ここでは、式１を次のように差分形式に書き換える。
Δｎ／Δｔ＝ＢｇＴＮ０／η＋Ｓｇφｎγ−Ｓｈ（ｄ０３ｎ／Ｎｂ）（ｎ／（１−ε）−１）ηｓγ２                      （式１９）
ここで、Δｎ：新しい時刻でのｎの近似変化分である。これから次式が得られる。
Δｎ＝（ＢｇＴＮ０／η＋Ｓｇφｎγ−Ｓｈ（ｄ０３ｎ／Ｎｂ）（ｎ／（１−ε）−１）ηｓγ２  ）Δｔ                  （式２０）
【００５１】
右辺はすべて値が入力されており、式２０によりΔｔ経過後のΔｎが近似的に計算できる。分散進行の場合はΔｎは負になっている。最初に計算済みの初期クラスター内粒子数ｎにこのΔｎを加えると新しい時刻におけるｎが計算できる。
【００５２】
次にクラスター内粒子数解析部１４においてｎ≒ｎ０／２か？という判定が行われる（２１５）。この判定はクラスター内粒子数が比較用クラスター内粒子数ｎ０の半分に達したら、クラスターが新たに二つに分裂し、クラスター切断回数Ｎｂを１つ増やすために用いられる。数２０をそのまま計算すると小数も出てくるのでｎがほぼｎ０／２とみなせる状態での判定方式としている。
【００５３】
最初はＳ２０７において初期ｎが比較用クラスター内粒子数ｎ０になっており、まず、この値と比較される。ここでＮｏの場合はまたＳ２０７に戻り、ｎ０は初期ｎのまま、Ｎｂも初期にセットした値のままとなる。Ｙｅｓの場合はクラスターが分裂した状態なので、ｎ０／２を新しいｎ０にセットしなおす。また、切断回数が１回増えているので、Ｎｂ＋１を新しいクラスター切断回数としてＮｂにセットしなおす。この手順を繰り返し、Ｓ２１１またはＳ２１２の判定でＮｏとなる状態で計算を終了する。
【００５４】
次に、具体的な解析例を説明する。図６はあるエポキシ樹脂の粘度パラメータを用いたときの攪拌槽の温度とη／η０が１．５に達する時間ｔの関係の計算結果である。攪拌槽の温度が高くなるほど反応の進行が早いため短い時間で粘度が上昇する。ここではＳ２１２の所定値は１．５とした。この場合、各温度においてグラフの下側が攪拌可能時間となる。
【００５５】
図７は、攪拌槽の温度を７５℃とした場合の時間ｔとη／η０の関係の計算結果である。この温度は後述する攪拌解析で設定した温度条件である。時間の経過とともに反応が進行するためη／η０は初期の１より増加し、約１６０００ｓで１．５に達する。したがって、この温度条件では１６０００ｓまでは攪拌可能時間となる。
【００５６】
図８は、攪拌槽を回転させ、攪拌槽内でせん断速度γがもっとも大きかった場所におけるクラスター内粒子数の変化を計算した結果である。ここでは４種類の回転数での結果をまとめて示した。それぞれの条件でのせん断速度γは２００，５００，１０００，２０００ｓ−１となっている。ここでは図６の特性のエポキシ樹脂に粒子径ｄ０が１２ｎｍの粒子をφ０＝０．０２の配合量、攪拌槽の温度７５℃で攪拌する条件を用いた。γが２００ｓ−１ではｎは時間の経過とともに緩やかに減少し、γの増加に伴いｎは小さくなる。
【００５７】
また、γが大きくなると、式１の特性により攪拌開始直後に急激にｎが小さくなり時間の経過とともにｎは飽和する。図８での攪拌時間は３６００ｓであり、攪拌限界時間である１６０００ｓより小さいのでこの時間内では問題なく攪拌できる。なお、攪拌槽内では場所によりせん断速度γは異なるが、攪拌はせん断速度のもっとも大きい場所を起点として対流輸送により速やかに行われるので、各時刻において場所によるｎの違いは非常に小さい。
【００５８】
図９も攪拌槽を回転させ、攪拌槽内の一箇所におけるクラスター内粒子数の変化を計算した結果である。ここでは粒子径ｄ０が５，１２，２５，５０ｎｍの４種類を用いたときの結果をまとめて示した。ここでは、図６の特性のエポキシ樹脂に粒子をφ０＝０．０２の配合量、攪拌槽の温度７５℃で攪拌し、γが１０００ｓ−１となる回転数での条件を用いた。Ｄ０が小さいと凝集しやすいのでｎの減少は少ないがｄ０の増加とともにｎは急激に減少する。Ｄ０＝５０ｎｍではほぼｎ＝１近くになっており、完全分散に近い状態が計算されている。
【００５９】
このように、本発明では分散に影響するパラメータを詳細に机上検討でき、粒子分散のプロセス設計を合理的に行うことができる。
【実施例２】
【００６０】
次に、本発明の実施例２の設計支援装置及び３次元流動解析処理について説明する。実施例１では粘度は温度と時間の関数として表し、ゲル化時刻までの粘度変化を計算している。このような熱硬化性樹脂の物性変化は反応の進行によるものであり、反応率の変化を表す式を基礎にして粘度計算や流動解析を行なうのが本発明の実施例２である。
【００６１】
図１０に、熱硬化性樹脂の等温状態での物性値の変化を示す。熱硬化性樹脂は時間の経過とともに反応が進行し、反応率αと分子量が増えていく。これにより、等温状態では液体の粘度ηが上昇する。流動が可能なのはまだ粘度が低い領域にあるときである。粘度は指数関数的に増加を続け、反応率がゲル化反応率αｇｅｌに達すると粘度は無限大になり、ゲル化が起きる。また、αｌｉｍは攪拌限界反応率であり、αｇｅｌより小さい値に設定しておく。実施例２では各時刻ごとに計算されるαの値をαｌｉｍと比較することにより、実施例１でのＳ２１２に相当する判定ができる。
【００６２】
図１１に、本発明の実施例２の設計支援装置の概略図を示す。実施例２においても設計支援装置の概略は図１と同じである。異なる点は、熱硬化性樹脂粘度解析部１５に反応率の計算を行う反応率計算部１７が加わり、温度、せん断速度解析部１６に反応による発熱速度の計算を行う発熱速度計算部１８が加わることである。実施例２において設計支援装置の実現手段も図２と同じである。
【００６３】
図１２は、本発明の実施例２の設計支援装置の３次元流動解析処理を説明するフロー図である。流動解析部１３ではＧＵＩ部１１を介してユーザより３次元攪拌解析に必要な、熱硬化性樹脂の物性値を受け取る（Ｓ３０１）。本実施形態では後述する熱反応式や粘度計算式中の係数、比熱、密度、熱伝導率などである。
【００６４】
次に、熱硬化性樹脂流動解析部１３では、ＧＵＩ部１１を介してユーザより３次元流動解析のための各種条件（境界条件、解析条件、および初期条件）を受け取る（Ｓ３０２）。各種条件には、攪拌槽の温度、形状データ、回転速度、使用樹脂量、攪拌される粒子の大きさ、粒子の含有率などが含まれる。
【００６５】
モデル作成部１２は、ＧＵＩ部１１を介して、ユーザより３次元流動解析指示を受け着ける。それから指定されたモデルデータを、モデル作成部１２から取り込む（Ｓ３０３）。次に、モデル作成部１２は、取り込んだデータが特定する樹脂領域を３次元解析領域に設定する。そして、この３次元解析領域を複数の３次元ソリッド要素に分割する際の条件（分割数や要素サイズなど）を、ＧＵＩ部１１を介してユーザより受け取る（Ｓ３０４）。そして、受け付けた分割条件に従い、３次元解析領域を複数の３次元ソリッド領域に分割する（Ｓ３０５）。Ｓ３０６，Ｓ３０７，Ｓ３０８の内容はそれぞれ実施例１のＳ２０６，Ｓ２０７，Ｓ２０８と同じである。
【００６６】
次に、流動解析部１３の熱硬化性樹脂粘度解析部１５の具体例について説明する。まず、時間ｔを初期時間にセットする。それから、熱硬化性樹脂用の熱反応式と温度条件を用いて３次元ソリッド要素毎に時間ｔにおける反応率と発熱速度を算出する（Ｓ３０９）。
熱反応式は式２１〜式２５で表される。
【００６７】
∂α／∂ｔ＝（Ｋ１＋Ｋ２αＭ）（１−α）Ｎ      （式２１）
Ｋ１＝Ｋａ  ｅｘｐ（−Ｅａ／Ｔ）                （式２２）
Ｋ２＝Ｋｂ  ｅｘｐ（−Ｅｂ／Ｔ）                （式２３）
α＝Ｑ／Ｑ０                                    （式２４）
∂Ｑ／∂ｔ＝Ｑ０（Ｋ１＋Ｋ２αＭ）（１−α）Ｎ  （式２５）
【００６８】
ここで、αは反応率、ｔは時間、Ｔは温度、∂α／∂ｔは反応速度、Ｋ１、Ｋ２は温度の関数で表される係数、Ｎ、Ｍ、Ｋａ  Ｋｂ  Ｅａ  Ｅｂは材料の固有係数、Ｑは時刻ｔまでの発熱量、Ｑ０は反応終了までの総発熱量、∂Ｑ／∂ｔは発熱速度を示している。そのうち、Ｎ、Ｍ、Ｋａ  、Ｋｂ  、Ｅａ  、Ｅｂ、Ｑ０はステップＳ３０１で受け付けた熱硬化性樹脂の物性値である。また、温度ＴはステップＳ３０２で受け付けた成形条件である。数２１〜数２５に初期状態から微小時間Δｔ経過毎の温度Ｔを逐次代入していけば反応速度、発熱速度の時間変化が計算でき、反応速度を微小時間経過Δｔ毎に時間方向に近似積分すれば反応率の時間変化が計算できる。
【００６９】
次に、粘度式と反応率、温度条件を用いて要素毎に時間ｔにおける分散媒粘度を算出する（Ｓ３１０）。
【００７０】
また、粘度式は式２６〜式２８で表される。
η＝η０（（１＋α／αｇｅｌ）／（１−α／αｇｅｌ））Ｃ  （式２６）
η０＝ａ  ｅｘｐ（ｂ／Ｔ）                                （式２７）
Ｃ＝ｆ／Ｔ−ｇ                                            （式２８）
【００７１】
ここで、ηは分散媒粘度、Ｔは温度、η０は初期分散媒粘度、αは反応率、αｇｅｌはゲル化時の反応率、Ｃは粘度上昇係数であり、ａ，ｂ，ｆ，ｇならびにαｇｅｌは材料の固有係数である。ａ，ｂ，ｆ，ｇならびにαｇｅｌはステップＳ３０１で受け付けた熱硬化性樹脂の物性値である。ａ，ｂ，ｆ，ｇ、αｇｅｌの値と温度、ならびにＳ３０９で計算されたαの値を式２６〜式２８に代入すれば分散媒粘度が計算できる。
【００７２】
この手法を用い、微小時間Δｔ変化毎に、温度条件とそのときの反応率を逐次代入していけば、等温状態では図４と、非等温状態では図５−（ｂ）と同様の熱硬化性樹脂特有の分散媒の粘度変化が計算できる。分散媒粘度ηを式４に代入してスラリー粘度ηｓが計算できる。Ｓ３１０で計算されたηとηｓはＳ３０８に送られ、数１中に値がセットされる。
【００７３】
次に、質量、運動量、エネルギー保存方程式を用いて、要素毎に、時間ｔにおける温度、速度、せん断速度などを算出する（Ｓ３１１）。なお、質量、運動量の保存方程式は、実施例１の式１６、式１７と同じになる。
また、エネルギー保存方程式は、式２９で示される。
ρＣｐ（∂Ｔ／∂ｔ  ＋  ｖ・▽Ｔ）＝λ▽２Ｔ  ＋τ：▽ｖ＋ρ（ｄＱ／ｄｔ）
（式２９）
【００７４】
ここで、ρは密度、Ｃｐは定圧比熱、Ｔは温度、ｔは時間、ｖは速度ベクトル、▽はナブラ演算子、λは熱伝導率、τは偏差応力テンソル、Ｑは発熱量である。式２９では実施例１の式１４に発熱速度ｄＱ／ｄｔが加わっている。この発熱速度はＳ３０９で計算されており、これを用いて式２９により温度計算を行う。すなわち、実施例２では熱硬化性樹脂の反応発熱を含んだ解析が可能となり、より正確に解析ができる。ここで計算された温度とせん断速度はＳ３０８に送られ式１にセットされる。
【００７５】
Ｓ３１２の判定内容は実施例１のＳ３１１と同じである。Ｓ３１３では、Ｓ３０９で計算された反応率αが攪拌限界反応率αｌｉｍと比較され、αがαｌｉｍに達した場合は解析が終了となる。αがαｌｉｍによりも小さい場合には時間ｔをタイムステップΔｔ進める（Ｓ３１４）。Ｓ３１５，Ｓ３１６，Ｓ３１７は、それぞれ実施例１のＳ２１４，Ｓ２１５，Ｓ２１６と同じになる。
【００７６】
以上述べた手法により、反応発熱の影響を考慮した上で、粒子の分散に影響するパラメータを詳細に机上検討でき、粒子分散のプロセス設計を合理的に行うことができる。
【産業上の利用可能性】
【００７７】
本発明によれば、クラスター内粒子数の変化は樹脂の熱・流動解析で得られる任意の場所と時間における温度、粘度、せん断速度などの値から計算できる。すなわち、粒子形状や数を幾何学的にモデル化せずに目的の計算が可能になるので、迅速なシミュレーションを行うことができる。また、熱硬化性樹脂特有の反応の進行による粘度変化を考慮するので高精度のシミュレーションが可能となり、実験前段階での粒子攪拌プロセスの最適化を図れるため、産業上の利用可能性は極めて高い。
【符号の説明】
【００７８】
１１    ＧＵＩ部
１２    モデル作成部
１３    流動解析部
１４    クラスター内粒子数解析部
１５    熱硬化性樹脂粘度解析部
１６    温度、せん断速度解析部
１７    反応率計算部
１８    発熱速度計算部
２１    ＣＰＵ
２２    メモリ
２３    外部記憶装置
２４    記憶媒体
２５    読取装置
２６    入力装置
２７    表示装置
２８    通信装置
【図１】
2010184365.tif 000003
【図２】
2010184365.tif 000004
【図３】
2010184365.tif 000005
【図４】
2010184365.tif 000006
【図５】
2010184365.tif 000007
【図６】
2010184365.tif 000008
【図７】
2010184365.tif 000009
【図８】
2010184365.tif 000010
【図９】
2010184365.tif 000011
【図１０】
2010184365.tif 000012
【図１１】
2010184365.tif 000013
【図１２】
2010184365.tif 000014