(57)【要約】【課題】  格子要素法を用いた定常的なスナップショット解析であっても、溶融樹脂材料の自由表面形状を再現できる流動挙動予測装置、その方法およびプログラムを提供する。【解決手段】  スクリュ式の混練装置により混練される流体状の材料の流動挙動を予測する流動挙動予測装置１０であって、材料の物性と、可塑化装置の操作条件および構成データとを使用して、格子要素法に基づいて複数の要素により構築されるスクリュ形状データを作成し、充満状態にある該材料の定常状態における流動解析を行い、少なくとも複数の要素を構築する各節点における圧力値を算出する解析部１０２と、算出された各節点の圧力値に基づいて、各節点に材料が存在するか否かを判定する判定部１０３と、材料が存在すると判定された第１節点と、存在しないと判定された第２節点との間に自由表面位置を設定し、該自由表面位置を連結させて前記材料の自由表面形状を作成する評価部１０４とを備えた。【選択図】図３
2016078397.tif 000002
【特許請求の範囲】
【請求項１】
スクリュ式の混練装置により混練される流体状の材料の流動挙動を予測する流動挙動予測装置であって、
前記材料の物性と、前記混練装置の操作条件および構成データとを使用して、格子要素法に基づいて複数の要素により構築されるスクリュ形状データを作成し、充満状態にある該材料の定常状態における流動解析を行い、少なくとも前記複数の要素を構築する各節点における圧力値を算出する圧力値算出部と、
算出された各節点の圧力値に基づいて、各節点に前記材料が存在するか否かを判定する節点判定部と、
前記材料が存在すると判定された第１節点と、存在しないと判定された第２節点との間に自由表面位置を設定し、該自由表面位置を連結させて前記材料の自由表面形状を作成する形状作成部と
を備えることを特徴とする流動挙動予測装置。
【請求項２】
前記判定部は、節点の圧力値が正の値である場合に該節点が第１節点であると判定し、節点の圧力値が負の値である場合に該節点が第２節点であると判定することを特徴とする請求項１記載の流動挙動予測装置。
【請求項３】
前記判定部は、節点の圧力値が大気圧より高い場合に該節点が第１節点であると判定し、節点の圧力値が大気圧より低い場合に該節点が第２節点であると判定することを特徴とする請求項１記載のことを特徴とする請求項１記載の流動挙動予測装置。
【請求項４】
作成された自由表面形状が連続的であるか否かを判定し、該自由表面形状が連続的でないと判定された場合、連続的でないと判定された箇所に該当する要素に樹脂が存在しているとして該自由表面形状を修正する修正部を更に備えることを特徴とする請求項１〜請求項３の何れか１項に記載の流動挙動予測装置。
【請求項５】
スクリュ式の混練装置により混練される流体状の材料の流動挙動を予測する流動挙動予測装置が行う流動挙動予測方法であって、
前記材料の物性と、前記混練装置の操作条件および構成データとを使用して、格子要素法に基づいて複数の要素により構築されるスクリュ形状データを作成し、充満状態にある該材料の定常状態における流動解析を行い、少なくとも前記複数の要素を構築する各節点における圧力値を算出するステップと、
算出された各節点の圧力値に基づいて、各節点に前記材料が存在するか否かを判定するステップと、
前記材料が存在すると判定された第１節点と、存在しないと判定された第２節点との間に自由表面位置を設定し、該自由表面位置を連結させて前記材料の自由表面形状を作成するステップと
を備えることを特徴とする流動挙動予測方法。
【請求項６】
スクリュ式の混練装置により混練される流体状の材料の流動挙動を予測する処理をコンピュータに実行させるための流動挙動予測プログラムであって、
前記コンピュータを、
前記材料の物性と、前記混練装置の操作条件および構成データとを使用して、格子要素法に基づいて複数の要素により構築されるスクリュ形状データを作成し、充満状態にある該材料の定常状態における流動解析を行い、少なくとも前記複数の要素を構築する各節点における圧力値を算出する圧力値算出部と、
算出された各節点の圧力値に基づいて、各節点に前記材料が存在するか否かを判定する節点判定部と、
前記材料が存在すると判定された第１節点と、存在しないと判定された第２節点との間に自由表面位置を設定し、該自由表面位置を連結させて前記材料の自由表面形状を作成する形状作成部
として機能させるための流動挙動予測プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、可塑化装置、特に樹脂等の材料を収容するシリンダ内に１つまたは複数のスクリュが回動自在に設けられたスクリュ式の押出成形機といった混練装置における、当該材料の溶融状態における流動解析（樹脂流動解析）を行うと共にその流動挙動を予測する流動挙動予測装置、その流動挙動予測方法および流動挙動予測プログラムに関するものである。
【背景技術】
【０００２】
従来、上記のような押出成形機により溶融可塑化される樹脂材料の流動挙動を予測する手法として、有限要素法や有限体積法、粒子法など様々な手法が用いられている。有限要素法や有限体積法は、流動場を規格化された格子状の要素に分割し、それぞれの要素内に存在する節点を計算点とし、当該計算点の物理量を連立方程式にて解く手法である。これらの演算手法は総じて格子要素法と称されている。
【０００３】
この格子要素法による解析は、樹脂材料の完全溶融および完全充満状態（シリンダ内に樹脂が充満した状態）を前提としたスクリュの混練性能予測や樹脂材料の発熱量、滞留時間などの物理量予測に適しており、ニュートン流体および純粘性非ニュートン流体を対象とした演算・評価技術は成熟の域に達している（下記、特許文献１および非特許文献１参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平１０−２４４５７９号公報
【非特許文献】
【０００５】
【非特許文献１】Ｋ．Ｆｕｎａｔｓｕ，  Ｓ．Ｋｉｈａｒａ，  Ｍ．Ｍｉｙａｚａｋｉ，  Ｓ．Ｋａｔｓｕｋｉ  ａｎｄ  Ｔ．Ｋａｊｉｗａｒａ，  Ｐｏｌｙｍ．Ｅｎｇ．Ｓｃｉ．，４２，７０７（２００２）
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、上述したように格子要素法は、シリンダ内に樹脂が充満した充満状態を前提としたものであるため、シリンダ内に樹脂が充満しておらず、溶融した樹脂材料が流体となって自由表面を形成する非充満状態の定常的なスナップショット解析による予測は困難であるという問題があった。
【０００７】
本発明は上述した問題点を解決するためになされたものであり、格子要素法を用いた定常的なスナップショット解析であっても、流体である樹脂材料の自由表面形状を再現できる流動挙動予測装置、その流動挙動予測方法および流動挙動予測プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
上述した課題を解決するため、本発明の一態様は、スクリュ式の混練装置により混練される流体状の材料の流動挙動を予測する流動挙動予測装置であって、前記材料の物性と、前記混練装置の操作条件および構成データとを使用して、格子要素法に基づいて複数の要素により構築されるスクリュ形状データを作成し、充満状態にある該材料の定常状態における流動解析を行い、少なくとも前記複数の要素を構築する各節点における圧力値を算出する圧力値算出部と、算出された各節点の圧力値に基づいて、各節点に前記材料が存在するか否かを判定する節点判定部と、前記材料が存在すると判定された第１節点と、存在しないと判定された第２節点との間に自由表面位置を設定し、該自由表面位置を連結させて前記材料の自由表面形状を作成する形状作成部とを備える。
【０００９】
また、本発明の一態様は、スクリュ式の混練装置により混練される流体状の材料の流動挙動を予測する流動挙動予測装置が行う流動挙動予測方法であって、前記材料の物性と、前記混練装置の操作条件および構成データとを使用して、格子要素法に基づいて複数の要素により構築されるスクリュ形状データを作成し、充満状態にある該材料の定常状態における流動解析を行い、少なくとも前記複数の要素を構築する各節点における圧力値を算出するステップと、算出された各節点の圧力値に基づいて、各節点に前記材料が存在するか否かを判定するステップと、前記材料が存在すると判定された第１節点と、存在しないと判定された第２節点との間に自由表面位置を設定し、該自由表面位置を連結させて前記材料の自由表面形状を作成するステップとを備える。
【００１０】
また、本発明の一態様は、スクリュ式の混練装置により混練される流体状の材料の流動挙動を予測する処理をコンピュータに実行させるための流動挙動予測プログラムであって、前記コンピュータを、前記材料の物性と、前記混練装置の操作条件および構成データとを使用して、格子要素法に基づいて複数の要素により構築されるスクリュ形状データを作成し、充満状態にある該材料の定常状態における流動解析を行い、少なくとも前記複数の要素を構築する各節点における圧力値を算出する圧力値算出部と、算出された各節点の圧力値に基づいて、各節点に前記材料が存在するか否かを判定する節点判定部と、前記材料が存在すると判定された第１節点と、存在しないと判定された第２節点との間に自由表面位置を設定し、該自由表面位置を連結させて前記材料の自由表面形状を作成する形状作成部として機能させる。
【発明の効果】
【００１１】
本発明によれば、格子要素法を用いた定常的なスナップショット解析であっても、溶融樹脂材料の自由表面形状を再現することができる。
【図面の簡単な説明】
【００１２】
【図１】本実施の形態に係る流動挙動予測装置の予測対象となるスクリュ式押出機におけるスクリュ混練場を模式的に示す平面図である。
【図２】本実施の形態に係る流動挙動予測装置のハードウェア構成を示すブロック図である。
【図３】本実施の形態に係る流動挙動予測装置の機能構成を示す機能ブロック図である。
【図４】８節点６面体要素にてマーチングキューブ法を適用した自由表面形状の描画パターンを示す図である。
【図５】本実施の形態に係る流動挙動予測処理を示すフローチャートである。
【図６】流動解析処理を示すフローチャートである。
【図７】解析評価処理を示すフローチャートである。
【図８】二軸スクリュ押出成形機の流動解析によるスクリュ表面の樹脂圧力分布結果を示す図である。
【図９】二軸スクリュ押出成形機の流動解析によるシリンダ表面の樹脂圧力分布結果を示す図である。
【図１０】図８のＡ−Ａ断面図である。
【図１１】二軸スクリュ押出成形機の流動解析による溶融樹脂材料の充満状態を描画した結果を示す図である。
【図１２】図８のＡ−Ａ断面における溶融樹脂材料の充満状態を描画した結果を示す図である。
【図１３】流動挙動予測プログラムが情報処理装置に適用される場合を示す図である。
【発明を実施するための形態】
【００１３】
以下、本発明の実施の形態について図面を参照しつつ説明する。
【００１４】
本実施の形態においては、図１に示されるような、樹脂材料を収容するシリンダ１と、樹脂材料を溶融混練すると共にシリンダ１の先端部（下流側）へ搬送する２つのスクリュ２とを備える二軸スクリュ押出成形機３における、そのスクリュ混練場（スクリュ２及びシリンダ１内に形成される流路空間）にて溶融可塑化される樹脂材料の樹脂流動挙動を予測する流動挙動予測装置を例にとり説明を行う。本実施の形態においては、二軸スクリュ押出成形機を例に挙げ説明するが、これに限定されるものではなく、スクリュ式の可塑化装置、混練装置、単軸または３軸以上のスクリュ式の押出成形機においても本発明を適用することができる。以下、本実施の形態について図面を参照しつつ、その詳細を説明する。
【００１５】
（装置構成）
図２は、本実施の形態に係る流動挙動予測装置のハードウェア構成を示すブロック図である。図２に示されるように、流動挙動予測装置１０は、ＣＰＵ(Central Processing Unit)１１、記憶部１２、入力部１３、表示部１４、ＨＤＤ（Hard Disk Drive）１５、を有する。
【００１６】
ＣＰＵ１１は記憶部１２上に展開されるＯＳ(Operating System)、ＢＩＯＳ(Basic Input/Output System)、アプリケーション等の各種プログラムを実行し、流動挙動予測装置１０の制御を行う。記憶部１２は、所謂ＲＡＭ(Random Access Memory)などの揮発性のメモリであり、実行されるプログラムの作業領域として利用される。
【００１７】
入力部１３は、流動挙動予測装置１０を使用するユーザからの入力を受け付けるものであり、例えば、ディスプレイ上の特定の位置を指定するためのポインティングデバイスであるマウスや、文字または特定の機能等が割り当てられた複数のキーが配列されたキーボードである。
【００１８】
表示部１４は、ＯＳおよびＯＳ上で動作するアプリケーションのＧＵＩ（Graphical User Interface）、後述する図８〜図１２に示される解析結果および描画データを表示するディスプレイ等の出力装置である。
【００１９】
ＨＤＤ１５は、後述する流動挙動予測処理において用いられる各種パラメータや当該処理により算出される各種物理量といったデータが格納される、所謂不揮発性の記憶領域である。
【００２０】
（機能構成）
次に、流動挙動予測装置１０の機能構成を説明する。図３は、本実施の形態に係る流動挙動予測装置の機能構成を示す機能ブロック図である。図３に示されるように、流動挙動予測装置１０は、取得部１０１と、解析部１０２と、判定部１０３と、評価部１０４と、出力部１０５とを、機能として有する。これら機能は、ＣＰＵ１１や記憶部１２等の前述したハードウェア資源が協働することにより実現される。
【００２１】
取得部１０１は、所定の物理量を算出するために必要な各種パラメータを取得するものである。なお、当該パラメータは、ユーザにより手入力された情報を取得してもよく、予測対象の二軸スクリュ押出成形機３に対応してＨＤＤ１５から取得するようにしてもよい。
【００２２】
解析部１０２は、取得部１０１により取得された各種パラメータに基づいて３次元のスクリュデータの作成や樹脂材料の流動挙動を解析するものであり、本実施の形態においては、前述した格子要素法を用いて樹脂材料の流動挙動を解析する樹脂流動解析を行う。判定部１０３は、流動挙動予測処理に係る各種判定処理を行うものである。判定処理としては、流動挙動予測処理に係る収束判定や、後述する格子要素の各節点に樹脂材料が存在するか否かの判定等が挙げられる。
【００２３】
評価部１０４は、樹脂流動解析結果を評価するものであり、具体的には、当該評価として流体となった樹脂材料の自由表面形状を作成する。出力部１０５は、作成された自由表面形状を含めて流動挙動予測処理の結果を、例えば圧力分布図等として表示部１４やＨＤＤ１５へ出力するものである。
【００２４】
（流動挙動予測処理の概要）
以下に、前述した流動挙動予測処理の理解を容易にするため、その概要を説明する。本実施の形態における流動挙動予測処理は、差分法を適用しない定常解析、即ち、時間のファクタを除いた定常状態を想定した樹脂流動解析を、格子要素法を用いて行うものである。具体的には、一般的な流動解析と同様、演算対象の演算領域であるスクリュ混練場（流路空間）の何れかの部位に圧力の境界条件を設定して３次元のスクリュ形状を具体的に再現し、その流路空間を三次元形状の格子要素に分割する。この分割した複数の格子要素それぞれを構築する節点を計算点として所定の方程式を解くことにより、各計算点における物理量を算出し、この物理量に基づいて樹脂材料の流動挙動を予測する。
【００２５】
なお、格子要素に分割する場合、その格子要素の形状（メッシュ）として、三角形の面から構成される四面体要素を採用してもよいが、四角形の面から構成される六面体要素を採用することが好ましい。これは、四面体要素はコンピュータ上でソフトウェアが自動的に演算領域を要素分割する際に採用されることが多いが、演算精度は六面体要素を採用した方が高いためである。したがって、本実施の形態においては、六面体要素を採用して格子要素への分割を行う。なお、他の多面体要素を採用してもよいことは言うまでもない。
【００２６】
このような格子要素法を用いた流動解析において、設定した圧力境界条件が１０ＭＰａなどの非常に高い値である場合は、演算領域全体の流体圧力が正の値を有する解析結果が得られるケースが多い。一方で、圧力境界条件の値が１ＭＰａ以下などの比較的大気圧に近い場合は、演算領域内の圧力結果が正の値と負の値の分布を示す。これは、与えた圧力境界条件が実現象に沿った現実的な値であっても、その結果は非現実的な負圧領域が存在することを意味する。即ち、この解析結果は演算領域内の数値バランスを保つために出力された数値であって、実際の現象と対比すると、負圧の領域は流体が存在しない空間であると判断することができる。つまり、流体が完全に充満しているとの前提で実施される格子要素法を適用した解析でも、負圧の流体が存在すると出力された領域は流体が存在しない空間であるとみなすことができる。本発明者らは、この考えに基づき、圧力値が正の値と負の値との境界面を流体の自由表面とすることにより、格子要素法を用いた定常的なスナップショット解析であっても、溶融樹脂材料の自由表面形状を再現できることを見出した。この考え方は解析の前提条件（完全充満）と結果の処理（非充満状態の判定）とが理論的に矛盾するため、同様の考え方に基づいた考察例はこれまでにはない。
【００２７】
このような自由表面形状を描画する手法としては、図４に示されるようなマーチングキューブ法（例えば、Ｗ．  Ｌｏｒｅｎｓｅｎ  ａｎｄ  Ｈ．  Ｃｌｉｎｅ，  Ｃｏｍｐｕｔｅｒ  Ｇｒａｐｈｉｃｓ  Ｖｏｌ．２１，  Ｎｏ．４，  ｐｐ．１６３‐１６９  （１９８７）．参照）に基づくことが好ましい。具体的には、前述したように六面体要素へ分割されると、格子要素の辺の頂点は８カ所となる。その頂点（節点）を計算点とし、この８つの計算点それぞれが正の圧力値を有する３２通りのパターンを考えることができる。それぞれのパターンで格子要素内にて自由表面を表現するには、マーチングキューブ法によれば１５通りのパターンのうちのいずれかで描画することが可能となる。この描画法で描かれた流体の自由表面および充満状態は自由表面近傍の流体圧力がわずかに正の値を有する比較的低い圧力であるが、混練用途のスクリュで構成される流体が完全充満し混練が促進される混練領域では数ＭＰａの高い圧力が作用する妥当な結果が表示できる。
【００２８】
また、負圧の演算領域は実際には０Ｐａ以下を取り得る領域場ではないため、その圧力は０と判断することができる。そのため、結果出力の際に圧力分布を表示させる場合、流体が存在しない空間領域は０Ｐａと表示することが可能である。ただし、この負圧領域の圧力値の修正は、実際の解析演算中に実施すると物理量の収支バランスが崩壊するため定常状態の収束解が得られなくなる。そのため、演算終了後に出力された結果のデータを修正することが好ましい。
【００２９】
（処理動作）
次に、流動挙動予測装置１０の前述した各機能により実行される流動挙動予測処理について、図面を用いてその詳細を説明する。図５は、本実施の形態に係る流動挙動予測処理を示すフローチャートである。図５に示されるように、先ず、取得部１０１は、後述する流動解析処理による物理量の算出に用いられる各種パラメータを取得する（Ｓ１，Ｓ２，Ｓ３）。取得方法としては、例えば流動挙動予測装置１０のユーザによる入力部１３を介した入力により、取得部１０１が各種パラメータを取得してもよく、予めＨＤＤ１５内に格納された各種パラメータを読み出して取得するようにしてもよい。
【００３０】
ステップＳ１で取得されるパラメータは、格子要素法において用いる樹脂物性に関するものであり、例えば、粘度フィッティングによる粘度モデル式のパラメータ等が挙げられる。ステップＳ２で取得されるパラメータは、二軸スクリュ押出成形機３の操作条件に関するパラメータであり、例えば、樹脂材料の押出量、スクリュ２の回転数、二軸スクリュ押出成形機３における樹脂材料の流出圧力または流入圧力、流体の流出量または流速、外部からの伝熱量（熱解析を実施する場合）などが挙げられる。ステップＳ３で取得されるパラメータは、二軸スクリュ押出成形機３の構成データに関するものであり、例えば、シリンダ径、スクリュ形状、スクリュ溝深さなどが挙げられる。この構成データは、スクリュ２の先端部分や図示しないホッパ近辺の供給部分など複数の演算領域における各値であることが好ましい。
【００３１】
各種パラメータの取得後、解析部１０２及び判定部１０３により、スクリュ形状を作成してその演算領域における樹脂流動解析を行う流動解析処理が実行され（Ｓ４）、流動解析処理後、評価部１０４により、前述した解析結果（計算点の物理量）を評価して自由表面形状を作成する解析評価処理が実行される（Ｓ５）。解析評価処理後、出力部１０５は評価結果を表示部１４に出力し（Ｓ６）、本フローは終了となる。なお、評価結果としては各計算点の充満率、圧力値、速度等の物理量が数値や分布図形式等で出力される。次に、流動解析処理、解析評価処理について具体的に説明する。
【００３２】
先ず、ステップＳ４の流動解析処理について説明する。図６は、流動解析処理を示すフローチャートである。流動解析処理は、一般的に用いられている格子要素法に基づいて、スクリュ２及び樹脂材料が収容されたシリンダ１内に形成されるスクリュ混練場（流路空間）を対象に樹脂材料の流動解析を行う処理であり、先ず解析部１０２により、３次元のスクリュ形状の再現が行われる（Ｓ４０１）。簡単に説明すると、一般的な樹脂流動解析と同様に、演算領域、即ちシリンダ１内の最上流部（流体流入部）に流体の流入量あるいは流速を与え、当該領域内のいずれかの部位（一般的には流出口）に圧力の境界条件を設定する。なお、この圧力境界条件は取得部１０１により予め取得されるようにしてもよい。この圧力境界条件に加え、スクリュの回転数や、熱解析を実施する場合には外部からの伝熱量の設定が施される。これらの設定値をもとに、スクリュ形状を具体的に再現する。スクリュ形状の詳細な再現方法は、一般的な手法であるため、その詳細は省略する。
【００３３】
スクリュ形状の再現後、解析部１０２は、その流路空間を三次元形状の格子要素、即ち六面体要素へと分割する（Ｓ４０２）。分割後、解析部１０２は、計算点に与える物理量を算出するために初期圧力値としてＰ０＝０と設定し（Ｓ４０３）、剪断速度や粘度等の物理量を算出する。例えば、シリンダ１とスクリュ２との間の流路空間でスクリュ２が回転することによって定められる剪断速度は、以下の（１）式により求めることが好ましい。
【００３４】
【数１】
2016078397.tif 000003
【００３５】
この（１）式においては、Ｄ：シリンダ径、Ｈ：溝深さ、Ｎ：スクリュ回転数、である。
【００３６】
また、粘度は樹脂材料を対象とした非ニュートン流体のモデル式を適用して算出することが好ましい。非ニュートンのモデル式は複数種類が提案されているが、例えば、そのうちの代表モデルであるＰｏｗｅｒ−ｌａｗモデルを例にすると、以下の（２）式により求めることができる。
【００３７】
【数２】
2016078397.tif 000004
【００３８】
この（１）式においては、η：粘度、ｋ，ｎ：物性パラメータ、である。
【００３９】
計算点に付与する物理量の算出後、解析部１０２は、算出された剪断速度や粘度等の物理量を用いてＮａｖｉｅｒ−Ｓｔｏｋｅｓ方程式により各計算点の圧力値、速度、滞留時間、充満率等の物理量を算出する（Ｓ４０５）。なお、計算点の物理量を求められる方程式であれば、何れの方程式を用いてもよい。
【００４０】
算出後、判定部１０３は、収束判定のために、各計算点における圧力値の差分であるΔＰをΔＰ＝Ｐ−Ｐ０により求める（Ｓ４０６）。ＰはＮａｖｉｅｒ−Ｓｔｏｋｅｓ方程式により算出された圧力値である。各計算点の圧力値差分ΔＰの算出後、判定部１０３は、全ての計算点の圧力値差分ΔＰが所定の収束判定値未満であるか否かを判定する（Ｓ４０７）。収束判定値としては、本実施の形態のような高粘度の樹脂流動解析では、比較的圧力値の精度が求められるため、小数点以下の値を設定することが好ましく、演算時間の観点から０．１〜０．００１Ｐａの範囲で設定することがより好ましい。
【００４１】
全ての計算点の圧力値差分ΔＰが判定収束値未満でないと判定された場合（Ｓ４０７，ＮＯ）、判定部１０３は、Ｐ０＝ΔＰとし（Ｓ４０８）、再度ステップＳ４０４の物理量の算出処理を行う。一方、全ての計算点の圧力値差分ΔＰが判定収束値未満であると判定された場合（Ｓ４０７，ＹＥＳ）、本フローは終了となる。
【００４２】
次に、ステップＳ５の解析評価処理について説明する。図７は、解析評価処理を示すフローチャートである。解析評価処理は、前述した自由表面形状を形成するための処理であり、先ず評価部１０４により、各計算点のうちの何れか１つが選択される（Ｓ５０１）。選択後、判定部１０３は、選択した計算点の収束判定後における圧力値が０未満、即ち負の値か否かを判定する（Ｓ５０２）。圧力値が負の値であると判定された場合（Ｓ５０２，ＹＥＳ）、評価部１０４は当該圧力値を０とし（Ｓ５０３）、判定部１０３による全ての計算点を選択したか否かの判定に移行する（Ｓ５０４）。
【００４３】
一方、圧力値が負の値でないと判定された場合（Ｓ５０２，ＮＯ）、そのままステップＳ５０４の判定処理に移行する。全ての計算点を選択したと判定された場合（Ｓ５０４，ＹＥＳ）、評価部１０４は、圧力値が正の値である計算点と負の値である計算点との中間に対し、自由表面が存在する位置を示す自由表面位置を設け、これら計算点を連結し、自由表面形状を作成する（Ｓ５０５）。これにより、圧力値が正の値にある等値面を求めることができるため、ステップＳ６の評価結果出力の処理において、出力部１０５にてその等値面を三次元的に描画することにより流体である溶融樹脂材料の自由表面の表現を行うことができる。また、圧力分布も現実的な非充満状態の描画が可能となる。一方、全ての計算点を選択していないと判定された場合（Ｓ５０４，ＮＯ）、評価部１０４は、ステップＳ５０１の計算点の選択処理に移行し、選択していない計算点を選択する。
【００４４】
自由表面形状作成後、判定部１０３は、自由表面形状が連続的か否か、換言すると、自由表面形状が格子要素間で連続性を有しているか否かを判定する（Ｓ５０６）。実際の解析では、本来の自由表面とは別に樹脂材料が存在している内部に圧力が低い領域が生じ、気泡のような空間（自由表面）が形成される場合がある。そのような場合、本来の自由表面とその内部の空間のどちらが実際の自由表面であるかを認識し、正確な描画を行う必要がある。本実施の形態においては、自由表面形状が連続的か否かを判定し、連続的でないと判定された場合、即ち気泡のような空間が形成されていると判定された場合、自由表面を修正することにより、極めて正確な自由表面形状を得ることを可能としている。なお、「自由表面形状が連続的か否か」の判定は、その空間領域がシリンダ要素に接しているか否かを基準とすることが好ましい。これは、本来の自由表面の外側の気体存在領域がシリンダ表面と接触しており、樹脂内部で生じた単なる気泡である場合はシリンダ表面には接触していないとの考えによるものである。
【００４５】
自由表面形状が連続的であると判定された場合（Ｓ５０６，ＹＥＳ）、本フローは終了となる。一方、自由表面形状が連続的でないと判定された場合（Ｓ５０６，ＮＯ）、評価部１０４は、連続的でない箇所の格子要素が全て前述したような気泡といった空間領域であると判断し、当該格子要素が樹脂材料で満たされているとして自由表面形状を修正する自由表面形状修正処理を実行し（Ｓ５０７）、本フローは終了となる。
【実施例】
【００４６】
本実施の形態に係る流動挙動予測装置１０により適切に溶融樹脂材料の自由表面形状が描画できるか、下記の条件を与えて流動挙動予測処理を実施した。解析対象としてφ６９ｍｍの図１に示されるような同方向回転型二軸押出機にて、押出量を２００ｋｇ／ｈ、スクリュ回転数を２００ｒｐｍとし、ポリプロピレン溶融体（ＭＦＲ＝４．０）を流動させた条件とした。スクリュ形状はフライトリード６５ｍｍのフルフライトを２１１．２５ｍｍの長さとし、出口断面の平均樹脂圧力を０．５ＭＰａと設定した。図８〜図１２に当該条件下における流動挙動予測処理の結果を示す。
【００４７】
図８は、二軸スクリュ押出成形機３の樹脂流動解析によるスクリュ表面の樹脂圧力分布結果を示す図であり、図９は、二軸スクリュ押出成形機３の樹脂流動解析によるシリンダ表面の樹脂圧力分布結果を示す図である。これらの分布図は、本実施例に係る流動挙動予測処理結果であるため、当然負の値を有する圧力値は存在していない。スクリュ２が回転すると樹脂材料を下流側へ搬送させる力が生じるため、図８及び図９からわかるとおり下流側（図右側）へかけて圧力値が次第に上昇している妥当な結果が示されている。
【００４８】
図１０は、図８のＡ−Ａ断面図である。スクリュ回転により溶融樹脂材料を搬送させようとする表面（Ｐｕｓｈｉｎｇ側）近傍の圧力値が上昇し、その背面（Ｌｅａｄｉｎｇ側）は圧力値が低下して０Ｐａとなっていることが表現できている。
【００４９】
図１１は、二軸スクリュ押出成形機の流動解析による溶融樹脂材料の充満状態を描画した結果を示す図である。これは、正の圧力値を有する計算点からマーチングキューブ法を用いて溶融樹脂材料の存在領域を図示したものである。図８および図５と比較して分かるとおり、正の圧力値を示す領域に樹脂が存在する描画となっており、その結果は実際の現象を鑑みても十分に妥当と判断できる。
【００５０】
図１２は、図８のＡ−Ａ断面における溶融樹脂材料の充満状態を描画した結果を示す図である。フライトのＰｕｓｈｉｎｇ側にのみ樹脂が存在しており、その充満状態および自由表面形状は格子要素のサイズを考慮すると妥当である。なお、さらに細かく要素分割を行うとより実際の挙動に近い自由表面状態が描画できることが示唆される。
【００５１】
以上の結果から、本実施の形態に係る流動挙動予測処理のような、格子要素により演算された圧力結果をもとに樹脂材料の充満状態を描画する手段では、現実的に十分な表現が可能であることが示された。また、他の解析と比較しても自由表面形状を明確に再現できており、その解析及び評価時間も差分法を適用しない定常解析を格子要素法で行っているため、最長でも１０分程度の演算時間で迅速に実施することができた。
【００５２】
例えば、自由表面が形成される非充満状態の予測とその表現を行うための一手法としてＶＯＦ（Volume Of Fluid）法を適用した解析が知られている。ＶＯＦ法は、初期値として演算対象系全体（演算領域全体）の流体（溶融樹脂材料）の充満率を設定するか、または、流体の供給量に従って演算対象系の最上流の境界面から流体を徐々に供給する演算条件のもとで、時間ステップを進行させた繰り返し計算による非定常状態での非充満状態を予測することが可能である。しかしながら、前述したＶＯＦ法を適用した解析は、高粘性流体である溶融樹脂材料では解析自体の収束安定性が悪く、またＶＯＦ法の難点である気液界面、即ち自由表面形状の判別が不明確であり、さらには演算時間が必要な非定常解析を行わなければならないため、本実施の形態のように定常的なスナップショット解析による予測をしたい場合、非定常解析を進めることで流動状態が安定した定常状態を導き出す必要があり、当該予測をリーズナブルに行うことが困難である。
【００５３】
ＶＯＦ法と同様に、非充満状態の予測とその表現を行う一手法として粒子法を適用した解析が知られている。粒子法は、液体の流動性を予測する手段としてＭＰＳ（Moving Particle Simulation）法やＳＰＨ（Smoothed Particle Hydrodynamic）法のいずれかを用いることが一般的である。これらの演算手法は、いずれも計算点である粒子そのものの物理量を差分法で解くことで、溶融状態にある樹脂材料の粘性流動挙動を表現することができるものである。また、これらの手法は、計算点が時間ステップごとに移動する離散要素解法であることから、結果として自由表面形状の表現も可能であり、速度場や圧力場を解くと同時にシリンダ内の樹脂材料の非充満挙動が予測でき、溶融樹脂材料を対象とした解析自体の収束安定性がＶＯＦ法と比較して良いという利点がある。しかしながら、この粒子法を適用した解析も非定常挙動の予測手段であるため、ＶＯＦ法と同様、定常状態の評価を行うには演算対象系全体の物理量が安定するまでの長い時間の計算が必要である。
【００５４】
即ち、押出成形機内のスクリュ混練場における樹脂の非充満状態を定常的なスナップショット解析にて予測したい場合、非定常解析を進めることで流動状態が安定した定常状態を導き出す必要があるため、これらの演算手法を適用すると少なくとも数時間、延いては数日を要し、一般的な設計検証に必要とされる数十秒から数分の演算時間で非充満状態の予測を行うことが困難である。一方、本実施の形態によれば、前述したように明確に自由表面形状を再現でき、且つ、その演算時間も極めて短時間とすることが可能となる。
【００５５】
なお、本実施の形態においては、流動挙動解析処理にて圧力値が負の値にある計算点を０としたが、これに限定されるものではない。自由表面の形状は大気圧下での表面張力によって形成されるとの考えに基づいた場合、自由表面を形成する境界圧力を０Ｐａとせず１０１３２５Ｐａの大気圧としてもよい。この場合、計算点の圧力値が大気圧より高いか低いかを判定基準とした描画とすることが好ましい。
【００５６】
また、本実施の形態においては、圧力値に基づく収束判定を行っているがこれに限定されるものではなく、時間以外のファクタ、例えば速度や充満率といった物理量に基づいて収束判定を行ってもよい。
【００５７】
本発明は、その要旨または主要な特徴から逸脱することなく、他の様々な形で実施することができる。そのため、前述の実施の形態は、あらゆる点で単なる例示に過ぎず、限定的に解釈してはならない。本発明の範囲は、特許請求の範囲によって示すものであって、明細書本文には、何ら拘束されない。更に、特許請求の範囲の均等範囲に属する全ての変形、様々な改良、代替および改質は、全て本発明の範囲内のものである。
【００５８】
また、実施の形態にて述べた流動挙動予測装置１０における各種ステップを、可塑化シミュレーションプログラムとして、図１３に示されるような、コンピュータにより読み取り可能な可搬型の記録媒体８に記憶させ、当該記録媒体８を情報処理装置９に読み込ませることにより、前述した機能を情報処理装置９に実現させることができる。記録媒体８としては、例えば、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤディスク等）、磁気ディスク（ハードディスクドライブ等）、フラッシュメモリ、ＩＣカード、更にネットワークを介することで伝送可能な媒体等、コンピュータで読み取りや実行が可能な全ての媒体が含まれる。
【００５９】
なお、特許請求の範囲に記載の流動挙動予測装置は、例えば、前述の実施の形態における流動挙動予測装置１０である。圧力値算出部は、例えば解析部１０２であり、節点判定部は、例えば判定部１０３である。形状作成部は、例えば評価部１０４であり、修正部は、例えば判定部１０３および評価部１０４である。
【符号の説明】
【００６０】
１０  流動挙動予測装置、１０２  解析部、１０３  判定部、１０４  評価部。
【図１】
2016078397.tif 000005
【図２】
2016078397.tif 000006
【図３】
2016078397.tif 000007
【図５】
2016078397.tif 000008
【図６】
2016078397.tif 000009
【図７】
2016078397.tif 000010
【図１３】
2016078397.tif 000011
【図４】
2016078397.tif 000012
【図８】
2016078397.tif 000013
【図９】
2016078397.tif 000014
【図１０】
2016078397.tif 000015
【図１１】
2016078397.tif 000016
【図１２】
2016078397.tif 000017