JP2017007207A

JP2017007207A - 可塑化シミュレーション装置、その可塑化シミュレーション方法および可塑化シミュレーションプログラム

Info

Publication number: JP2017007207A
Application number: JP2015124948A
Authority: JP
Inventors: 一馬中川; Kazuma Nakagawa; 克之荒木; Katsuyuki Araki; 公一尾崎; Koichi Ozaki
Original assignee: Japan Steel Works Ltd
Current assignee: Japan Steel Works Ltd
Priority date: 2015-06-22
Filing date: 2015-06-22
Publication date: 2017-01-12
Anticipated expiration: 2035-06-22
Also published as: JP6185515B2

Abstract

【課題】極めて高い精度で可塑化能力を算出することができるシミュレーション装置、可塑化シミュレーション方法および可塑化シミュレーションプログラムを提供する。
【解決手段】スクリュ式の可塑化装置におけるシリンダ内の樹脂材料がソリッドベッド、メルトフィルム、およびメルトプールの何れか状態にある３ゾーンモデルに基づいて可塑化能力を算出する可塑化シミュレーション装置であって、樹脂材料の樹脂物性と可塑化装置の運転条件と該可塑化装置の構成データとを含むパラメータを取得する取得部１０１と、初期可塑化能力を取得し、初期可塑化能力、パラメータ、質量保存の式およびエネルギー保存の式を用いて固相率およびメルトフィルム厚みを含む物理量を算出し、パラメータおよび物理量に基づいて３次元流動計算を行い、可塑化能力を算出する解析部１０２及び判定部１０３とを備えた。
【選択図】図３

Description

本発明の実施形態は、射出成形機に代表される単軸または多軸スクリュ式の可塑化装置の可塑化行程をシミュレートする可塑化シミュレーション装置、その可塑化シミュレーション方法および可塑化シミュレーションプログラムに関するものである。

従来、この種の可塑化装置において溶融可塑化される樹脂材料の樹脂温度や固相率といった物理量を算出し、可塑化工程をシミュレートするシミュレーション技術が知られている。一般的にこのようなシミュレーションにおいては、単位時間において可塑化される材料の重量を示す可塑化能力（ｋｇ／ｈ）が物理量の算出のために用いられており、この可塑化能力を固定して、または、所定のスクリュ特性式を用いて算出した可塑化能力を用いて、物理量を算出することが広く実施されている（下記、特許文献１、非特許文献１及び２参照）。

特開２００７−００７９５１号公報

日本製鋼所技報６３号ＥｎｇｉｎｅｅｒｉｎｇＰｒｉｎｃｉｐｌｅｓｏｆＰｌａｓｔｉｃａｔｉｎｇＥｘｔｒｕｓｉｏｎ

しかしながら、シミュレーション上ではなく、実際の可塑化能力は数ある諸条件に応じて変動するものであるため、上述したように可塑化能力を固定してしまうと、精度の高い物理量の算出結果を得ることが困難であるという問題があった。一方、スクリュ特性式を用いて諸条件に応じた可塑化能力を算出するようにしても、当該スクリュ特性式に従った計算精度でしか可塑化能力を算出できず、可塑化能力を固定した場合と同様に精度の高い物理量の算出結果を得ることが困難であるという問題があった。

本発明は上述した問題点を解決するためになされたものであり、極めて高い精度で可塑化能力を算出することができるシミュレーション装置、可塑化シミュレーション方法および可塑化シミュレーションプログラムを提供することを目的とする。

上述した課題を解決するため、本発明の一態様は、スクリュ式の可塑化装置におけるシリンダ内の樹脂材料がソリッドベッド、メルトフィルム、およびメルトプールの何れか状態にある３ゾーンモデルに基づいて可塑化能力を算出する可塑化シミュレーション装置であって、前記樹脂材料の樹脂物性と前記可塑化装置の運転条件と該可塑化装置の構成データとを含むパラメータを取得する情報取得部と、初期可塑化能力を取得する初期可塑化能力取得部と、前記初期可塑化能力、前記パラメータ、質量保存の式およびエネルギー保存の式を用いて固相率およびメルトフィルム厚みを含む物理量を算出する物理量算出部と、前記パラメータおよび前記物理量に基づいて３次元流動計算を行い、可塑化能力を算出する可塑化能力算出部とを備える。

また、本発明の一態様は、スクリュ式の可塑化装置におけるシリンダ内の樹脂材料がソリッドベッド、メルトフィルム、およびメルトプールの何れか状態にある３ゾーンモデルに基づいて可塑化能力を算出する可塑化シミュレーション装置が実行する可塑化シミュレーション方法であって、前記樹脂材料の樹脂物性と前記可塑化装置の運転条件と該可塑化装置の構成データとを含むパラメータを取得するステップと、初期可塑化能力を取得するステップと、前記初期可塑化能力、前記パラメータ、質量保存の式およびエネルギー保存の式を用いて固相率およびメルトフィルム厚みを含む物理量を算出するステップと、前記パラメータおよび前記物理量に基づいて３次元流動計算を行い、可塑化能力を算出するするステップとを有する。

また、本発明の一態様は、スクリュ式の可塑化装置におけるシリンダ内の樹脂材料がソリッドベッド、メルトフィルム、およびメルトプールの何れか状態にある３ゾーンモデルに基づいて可塑化能力を算出することをコンピュータに実行させる可塑化シミュレーションプログラムであって、前記コンピュータを、前記樹脂材料の樹脂物性と前記可塑化装置の運転条件と該可塑化装置の構成データとを含むパラメータを取得する情報取得部と、初期可塑化能力を取得する初期可塑化能力取得部と、前記初期可塑化能力、前記パラメータ、質量保存の式およびエネルギー保存の式を用いて固相率およびメルトフィルム厚みを含む物理量を算出する物理量算出部と、前記パラメータおよび前記物理量に基づいて３次元流動計算を行い、可塑化能力を算出する可塑化能力算出部として機能させる。

本発明によれば、極めて高い精度で可塑化能力を算出することができる。

本実施の形態においてシミュレーションされる単軸スクリュ式射出成形機を示す概略断面図である。本実施の形態に係る可塑化シミュレーション装置のハードウェア構成を示すブロック図である。本実施の形態に係る可塑化シミュレーション装置の機能構成を示す機能ブロック図である。３ゾーンモデルを説明するための図である。３ゾーンモデルにおける流動制限を説明するための図である。本実施の形態に係る可塑化解析処理を示すフローチャートである。初期可塑化能力算出処理を示すフローチャートである。２次元流動解析処理を示すフローチャートである。３次元流動解析処理を示すフローチャートである。可塑化能力算出処理を示すフローチャートである。実機での実験結果、可塑化解析処理により算出された算出結果、及びスクリュ特性式のみを用いた一般的な解析処理により算出された算出結果の各可塑化能力を示す図である。各算出結果の可塑化能力と実機での実験結果の可塑化能力との一致度を示す図である。各算出結果の可塑化能力と実機での実験結果の可塑化能力との誤差の割合を示す図である。可塑化シミュレーションプログラムが情報処理装置に適用される場合を示す図である。

以下、本発明の実施の形態について図面を参照しつつ説明する。

本実施の形態においては、図１に示されるような、樹脂材料が投入されるホッパ１と、樹脂材料を溶融するためのシリンダ２と、樹脂材料を溶融混練すると共にシリンダ２先端部（下流側）へ搬送する単軸のスクリュ３とを備えるスクリュ式射出成形機４の定常状態における可塑化解析を行う可塑化シミュレーション装置を例にとり説明を行う。即ち、本実施の形態に係る可塑化シミュレーション装置は、スクリュ３が回転している回転状態の各種物理量を算出するものである。なお、本実施の形態においては、可塑化装置として射出成形機を例に挙げ説明するが、これに限定されるものではなく、スクリュ式の押出機や多軸スクリュ式の可塑化装置においても本発明を適用することが可能である。以下、本実施の形態について図面を参照しつつ、その詳細を説明する。

（装置構成）
図２は、本実施の形態に係る可塑化シミュレーション装置のハードウェア構成を示すブロック図である。図２に示されるように、可塑化シミュレーション装置１０は、ＣＰＵ(
ＣｅｎｔｒａｌＰｒｏｃｅｓｓｉｎｇＵｎｉｔ)１１、記憶部１２、入力部１３、表
示部１４、ＨＤＤ（Ｈａｒｄｄｉｓｋｄｒｉｖｅ）１５、を有する。

ＣＰＵ１１は記憶部１２上に展開されるＯＳ(ＯｐｅｒａｔｉｎｇＳｙｓｔｅｍ)、ＢＩＯＳ(ＢａｓｉｃＩｎｐｕｔ／ＯｕｔｐｕｔＳｙｓｔｅｍ)、アプリケーション等の各種プログラムを実行し、可塑化シミュレーション装置１０の制御を行う。記憶部１２は、所謂ＲＡＭ(ＲａｎｄｏｍＡｃｃｅｓｓＭｅｍｏｒｙ)などの揮発性のメモリであり、実行されるプログラムの作業領域として利用される。

入力部１３は、可塑化シミュレーション装置１０を使用するユーザからの入力（例えば後述する各種パラメータ）を受け付けるものであり、例えば、ディスプレイ上の特定の位置を指定するためのポインティングデバイスであるマウスや、文字または特定の機能等が割り当てられた複数のキーが配列されたキーボードである。

表示部１４は、ＯＳおよびＯＳ上で動作するアプリケーションのＧＵＩ（ＧｒａｐｈｉｃａｌＵｓｅｒＩｎｔｅｒｆａｃｅ）や解析結果を表示するディスプレイ等の出力装置である。このような出力装置により、例えば、解析結果として算出された物理量がヒストグラムといった形で表示される。ＨＤＤ１５は、後述する可塑化解析処理において用いられる各種パラメータや当該処理により算出される各物理量といったデータが格納される、所謂不揮発性の記憶領域である。

（機能構成）
次に、可塑化シミュレーション装置１０の機能構成を説明する。図３は、本実施の形態に係る可塑化シミュレーション装置の機能構成を示す機能ブロック図である。図３に示されるように、可塑化シミュレーション装置１０は、取得部１０１と、解析部１０２と、判定部１０３と、出力部１０４とを、機能として有する。これら機能は、ＣＰＵ１１や記憶部１２等の前述したハードウェア資源が協働することにより実現される。

取得部１０１は、所定の物理量を算出するために必要な各種パラメータを取得するものである。なお、当該パラメータとして、ユーザにより手入力された情報を取得してもよく、解析対象のスクリュ式射出成形機４に対応して固有に紐付けられた情報をＨＤＤ１５から取得するようにしてもよい。解析部１０２は、取得された各種パラメータに基づいてシミュレーションの解析値として可塑化能力や各物理量の算出を行うものである。判定部１０３は、可塑化解析処理における各種判定を行うものであり、出力部１０４は、解析部１０２の算出結果や判定部１０３の判定結果を受けて、これらの結果を表示部１４やＨＤＤ１５へ出力するものである。

（概要）
理解を容易にするために、本実施の形態に係る可塑化解析処理の概要を簡単に説明する。本実施の形態においては、溶融樹脂量（適宜、溶融量とも称する）、樹脂温度、固相率、メルトフィルム厚み、スクリュ動力、圧力等の各物理量の算出をＴａｄｍｏｒモデルに倣い、図４に示されるような３ゾーンモデルを導入してゾーン毎に分けて行われる。図４に示される参照符号５１は、スクリュ３及びスクリュ３のフライト５からなる領域において樹脂材料が固形で存在するソリッドベッドを示し、参照符号５２は、樹脂材料の溶融が促進されるメルトフィルムを示し、参照符号５３は、溶融状態にある樹脂材料である溶融樹脂が滞留するメルトプールを示している。

この３ゾーンモデルにおいて可塑化される樹脂材料は、図５に示されるように各矢印で示される流速ベクトルを持って流動している。なお、樹脂材料はスクリュ３により押し出されるため、ここではシリンダ２にも流速ベクトルが付されている。本来、メルトフィルム５２にある溶融樹脂は、シリンダ２と接触するためにその溶融が促進されるが、そこからメルトプール５３へ流入する流量、換言すると、メルトフィルム５２とメルトプール５３との境界にあるメッシュの流量は、当該メッシュのメルトフィルム５２の溶融樹脂の溶融量（ソリッドベッド５１の溶融量）を超えることはない。しかしながら、従来の解析手法では、計算上、当該流量が当該溶融量を超えてしまう可能性がある。そのため、本実施の形態においては、溶融樹脂の流速ベクトルのうち、メルトフィルム５２からメルトプール５３へ流入する流速ベクトル６０の流量が、メルトフィルム５２の溶融量を超えないよう制限する、流量制限を可塑化能力の算出に組み込んでいる。このように実際には起こりえない事象を制限することで、高い精度での可塑化能力の算出を実現している。

（処理動作）
以下、図面を用いて本実施の形態に係る可塑化シミュレーション装置１０の動作の詳細を説明する。図６は、本実施の形態に係る可塑化解析処理を示すフローチャートである。なお、この処理において算出される各値は、出力部１０４により適宜記憶部１２またはＨＤＤ１５に格納される。そのため、本フロー内においては各値の格納処理は割愛する。図６に示されるように、本実施の形態に係る可塑化解析処理は、先ず、後述する３次元流動解析処理において用いられる初期値の可塑化能力である初期可塑化能力を算出するための初期可塑化能力算出処理が実行される（Ｓ１）。その後、各物理量を算出するために用いられる後述する分割せん断発熱量及び後述する最大流量を算出するための２次元流動解析を行う２次元流動解析処理が実行され（Ｓ２）、初期可塑化能力と、分割せん断発熱量とに基づいて、可塑化能力および各種物理量を算出する３次元解析を行う３次元流動解析処理が実行され（Ｓ３）、本フローは終了となる。以下、これら各処理の詳細を順次説明する。

先ず、初期可塑化能力算出処理の詳細を説明する。図７は、初期可塑化能力算出処理を示すフローチャートである。図７に示されるように、先ず、取得部１０１は、可塑化能力や各物理量を算出するための樹脂物性、運転条件およびスクリュ式射出成形機４の構成データの各種パラメータを取得する（Ｓ１０１，Ｓ１０２，Ｓ１０３）。ステップＳ１０１において取得される樹脂物性のパラメータとしては、例えば、粘度フィッティングによるモデルパラメータ、固体および溶融体の密度、比熱、熱伝導率、融点、溶融熱量などが挙げられる。ステップＳ１０２において取得される運転条件のパラメータとしては、例えば、スクリュ回転数、背圧、シリンダ設定温度、シリンダ設定温度境界位置、スクリュ位置、メッシュ刻み数、計算サイクル数、計量樹脂量、原料樹脂温度などが挙げられる。ステップＳ１０３において取得される構成データのパラメータとしては、例えば、シリンダ径、スクリュ径、スクリュ溝深さ、スクリュリード、フライト幅、フライトクリアランスなどが挙げられる。この構成データは、スクリュ３の先端部分やホッパ２近辺の供給部分な
ど複数の計算領域（スクリュ３の各ポジション）における各値であることが好ましい。

以上の取得されたパラメータに基づいて、解析部１０２は初期可塑化能力を算出し（Ｓ１０４）、本フローは終了となる。ここで算出する初期可塑化能力は、単位時間当たりの樹脂材料（溶融樹脂）の重量（ｋｇ／ｈ）で示されるものであり、下記（１）式または（２）式で表わされるスクリュ特性式を用いて算出される。なお、（１）式は、シリンダ２内で樹脂材料がスクリュ３により牽引されている部分である牽引流部分における計算に用いられ、（２）式はその他の部分における計算に用いられる。

この（１）式および（２）式においては、Ｑ：体積流量、Ａ及びＢ：係数、Ｖ_ｂｚ：樹脂の流れ方向の移動速度、ｗ：スクリュ溝幅、ｈ：スクリュ溝深さ、φ：スクリュへリックス角、η：樹脂粘度、∂Ｐ／∂ｚ：圧力変化量である。また、（２）式における樹脂粘度ηは、下記（３）式および（４）式の粘度式により求められる。

この（３）式および（４）式においては、η：樹脂粘度、Ａ，ａ，ｂ，ｃ：粘度係数、ｇ：せん断速度、Ｔ_ｒ及びｎ：粘度パラメータ、Ｔ：樹脂温度である。これら（１）〜（４）式に、取得部１０１により取得された各種パラメータのうちの必要なパラメータが導入されることで、初期可塑化能力が算出される。ここでの必要なパラメータとは、スクリュ回転数、シリンダ径、スクリュ溝深さ、溶融体密度等である。

次に、前述した２次元流動解析処理の詳細を説明する。図８は、２次元流動解析処理を示すフローチャートである。図８に示されるように、先ず、解析部１０２は、取得部１０１により取得された各種パラメータを記憶部１２またはＨＤＤ１５から取得し（Ｓ２０１）、これらパラメータを用いた一般的な２次元流動計算を行い、溶融樹脂のメルトフィルム５２内における速度分布を算出する（Ｓ２０２）。算出後、解析部１０２は、当該速度分布に基づいて、下記（５）式で示されるメルトフィルム５２内におけるメッシュ毎のせん断速度を算出し、メルトフィルム５２内におけるせん断速度分布を算出する（Ｓ２０３）。

この（５）式においては、γ^・：せん断速度、ｄｕ：メッシュ間速度差、ｄｙ：１メッシュの距離（厚み）、である。せん断速度分布の算出後、このせん断速度分布に基づいて、メルトフィルム５２内におけるメッシュ毎の粘度、せん断発熱量を算出する（Ｓ２０４）。樹脂粘度は２次元解析の一般的な手法により求めてもよく、上記（３）式および（４）式の粘度式により求めてもよい。せん断発熱量は、下記（６）式により求めることができる。

この（６）式においては、Ｑ：せん断発熱量（［Ｗ］）、η：樹脂粘度、γ^・：せん断速度である。粘度、せん断発熱量の算出後、解析部１０２は、これらに基づいてメルトフィルム５２からソリッドベッド５１およびソリッドベッド５１内部の熱流束を算出する（Ｓ２０５）。メルトフィルム５２からソリッドベッド５１への熱流束は下記（７）式、ソリッドベッド５１内部への熱流束は下記（８）式により求めることができる。

（７）式および（８）式においては、ｑ_ｆｓ：メルトフィルムからソリッドベッドへの熱流束（［Ｗ／ｍ^２］）、ｋ_ｍ：溶融体の熱伝導率、Ｔ_ｆ：メルトフィルム温度、Ｔ_ｍ：樹脂融点、δ：メルトフィルム厚み、Ｑ：せん断発熱量、ｑ_ｓｓ：ソリッドベッド内部への熱流束（［Ｗ／ｍ^２］）、ｋ_ｓ：固体の熱伝導率、ｈ：スクリュ溝深さ、Ａ：パラメータ、Ｂ：パラメータ、Ｃ：パラメータである。

熱流束算出後、解析部１０２は、算出したせん断発熱量と熱流束とから、メルトフィルム５２における溶融樹脂の溶融量を算出し、これを単位時間で割ることにより、メルトフィルム５２からメルトプール５３へ流れる溶融樹脂の単位時間当たりの最大流量を算出し（Ｓ２０６）、本フローは終了となる。具体的には、せん断発熱量と熱流束とを用いて、下記（９）式を解くことにより最大流量ｍを算出する（したがって、溶融量はｍｄｔである）。

この（９）式においては、ｍ：最大流量（［ｍ^３／ｓ］）、ｑ_ｆｓ：メルトフィルムからソリッドベッドへの熱流束、ｑ_ｓｓ：ソリッドベッド内部への熱流束、Ｗ_ｓ：ソリッドベッド幅、ｄｚ：流れ方向メッシュ幅、ρ：樹脂密度、Ｃ_ｐ：樹脂比熱、Ｔ_ｍ：樹脂融点、Ｔ_ｓ：ソリッドベッド温度、Ｈ：樹脂融解熱である。

ここで、本実施の形態においては、メルトフィルム５２におけるせん断発熱量のＣ／（Ａ＋Ｂ＋Ｃ）がメルトフィルム５２内での温度上昇に用いられることとする。即ち、複数のメッシュのうち、メルトフィルム５２に対応するメッシュについては、上記（６）式で求めたせん断発熱量にＣ／（Ａ＋Ｂ＋Ｃ）を乗じ、後述する各物理量の算出においてこのＣ／（Ａ＋Ｂ＋Ｃ）としたせん断発熱量を用いる。せん断発熱量をＣ／（Ａ＋Ｂ＋Ｃ）とする理由は、以下のとおりである。

メルトフィルム５２からソリッドベッド５１への熱流束は、上記（７）式により求められる。これによれば、せん断発熱量はＡ／（Ａ＋Ｂ＋Ｃ）がソリッドベッド５１へ流れることがわかる。また、シリンダ２からメルトフィルム５２への熱流束は、下記（１０）式により求められる。これによれば、せん断発熱量はＢ／（Ａ＋Ｂ＋Ｃ）がシリンダ２へ流れることがわかる。よって、残りのＣ／（Ａ＋Ｂ＋Ｃ）のせん断発熱量がメルトフィルム５２内に残り、溶融樹脂の温度上昇に用いられることがわかる。したがって、メルトフィルム５２に対応するメッシュのせん断発熱量をＡ／（Ａ＋Ｂ＋Ｃ）、Ｂ／（Ａ＋Ｂ＋Ｃ）、Ｃ／（Ａ＋Ｂ＋Ｃ）と分割することで、現実に沿った熱量を与えることができ、これを可塑化能力、各物理量の算出に用いることで、精度の高い解析結果を得ることができる。以後、このメルトフィルム５２におけるＡ／（Ａ＋Ｂ＋Ｃ）、Ｂ／（Ａ＋Ｂ＋Ｃ）、Ｃ／（Ａ＋Ｂ＋Ｃ）と分割したせん断発熱量を、分割せん断発熱量と称する。

この（１０）式においては、ｑ_ｂｆ：シリンダからメルトフィルムへの熱流束、ｋ_ｍ：溶融体の熱伝導率、Ｔ_ｆ：メルトフィルム温度、Ｔ_ｍ：樹脂融点、δ：メルトフィルム厚み、Ｑ：せん断発熱量、Ａ：パラメータ、Ｂ：パラメータ、Ｃ：パラメータである。

次に、前述した３次元流動解析処理の詳細を説明する。図９は、３次元流動解析処理を示すフローチャートである。図９に示されるように、先ず、解析部１０２は、取得部１０１により取得された各種パラメータ、算出した初期可塑化能力、分割せん断発熱量を記憶部１２またはＨＤＤ１５から取得し（Ｓ３０１）、取得した初期可塑化能力を第１可塑化能力と設定する（Ｓ３０２）。

設定後、解析部１０２は、これらのデータに基づいて各物理量を算出する（Ｓ３０３）。本実施の形態においては、算出される物理量として溶融樹脂量、樹脂温度、固相率、メルトフィルム５２厚み、スクリュ動力、圧力が挙げられる。これら物理量は全計算領域（メッシュ毎）に渡って算出される。これらの物理量は、第１可塑化能力と、取得部１０１
により取得された各パラメータのうちの必要なパラメータと、質量保存則およびエネルギー保存則を用いて、更にメルトフィルム５２に対しては分割せん断発熱量を別途用いて算出する。具体的には、図４に示されるソリッドベッド５１における各種物理量の算出には下記（１１）式と（１２）式とが用いられる。

この（１１）式および（１２）式においては、ρ：樹脂密度、Ｘ：ソリッドベッド幅、ｈ：スクリュ溝深さ、δ_ｍ：メルトフィルム厚み、Ｈ：エンタルピー、ｖ_ｓｚ：ソリッドベッドの流れ方向移動速度、λ：融解熱、ｍ_ｆ：ソリッドベッドのメルトフィルム側融解質量流速、ｍ_ｐ：ソリッドベッドのメルトプール側融解質量流速、ｋ_ｓ：固体樹脂の熱伝導率、Ｔ_ｓ：ソリッドベッド温度、ｑ_ｆｓ：メルトフィルムからソリッドベッド界面への熱流束、ｑ_ｐｓ：メルトプールからソリッドベッド界面への熱流束、である。

また、図４に示されるメルトフィルム５２における各物理量の算出には下記（１３）式と（１４）式とが用いられる。

この（１３）式および（１４）式においては、ρ：樹脂密度、Ｘ：ソリッドベッド幅、δ_ｍ：メルトフィルム厚み、Ｈ：エンタルピー、ｖ_ｆｚ：メルトフィルム流れ方向速度、ｖ_ｆｘ：メルトフィルムかき出し方向速度、ｑ_ｆｓ：ソリッドベッド界面からソリッドベッド内部への熱流束、ｋ_ｍ：溶融体熱伝導率、Ｔ_ｆ：メルトフィルム温度、ｑ_Ｂｆ：シリンダからメルトフィルムへの熱流束、Ｑ_Ｇ：メルトフィルムでのせん断発熱量、ｍ_ｆ：ソリッドベッドのメルトフィルム側融解質量流速、である。なお、Ｑ_Ｇには分割せん断発熱量が導入される。

また、図４に示されるメルトプール５３における各物理量の算出には下記（１５）式と（１６）式とが用いられる。

この（１５）式および（１６）式においては、ρ：樹脂密度、ｗ：スクリュ溝幅、Ｘ：ソリッドベッド幅、ｈ：スクリュ溝深さ、Ｈ：エンタルピー、ｖ_ｐｚ：メルトプール流れ方向速度、δ_ｍ：メルトフィルム厚み、ｑ_ｐｓ：メルトプールからソリッドベッド界面への熱流束、ｋ_ｍ：溶融体熱伝導率、Ｔ_ｐ：メルトプール温度、ｖ_ｆｘ：メルトフィルムかき出し方向速度、ｑ_Ｂｐ：シリンダからメルトプールへの熱流束、Ｑ’_Ｇ：メルトプールせん断発熱量、ｍ_ｐ：ソリッドベッドのメルトプール側融解質量流速、である。

前述した（１１）式、（１３）式、および（１５）式はエネルギー保存則を示しており、（１２）式、（１４）式および（１６）式は質量保存則を示している。また、質量保存則およびエネルギー保存則をソリッドベッド５１と、メルトフィルム５２と、メルトプール５３とに分けたため、この質量保存則およびエネルギー保存則を満足させるように融解質量流速が導入されている。これらの数式に、取得された各パラメータのうちの必要なパラメータ、第１可塑化能力、分割せん断発熱量等を導入することで、溶融樹脂量、樹脂温度、固相率、メルトフィルム厚み、スクリュ動力、圧力等を算出することができる。

図９に戻り、ステップＳ３０３の各物理量の算出を終えると、解析部１０２は、固相率と、メルトフィルム厚みとから、メルトプール５３の形状メッシュを作成する（Ｓ３０４）。この形状メッシュの作成方法は、例えば有限体積法等の手法を用いて格子メッシュを作成すればよく、一般的な手法であるためここでの説明は省略する。形状メッシュの作成後、解析部１０２は、可塑化能力算出処理を実行する（Ｓ３０５）。可塑化能力算出処理は、可塑化能力の算出に流動計算を用いるものであり、この処理により算出されるメルトフィルム５２からメルトプール５３への溶融樹脂の流出量を第２可塑化能力とするものである。この可塑化能力算出処理についての詳細は後述する。

第２可塑化能力算出後、判定部１０３は、算出した第２可塑化能力が所定の条件を満たすか否かを判定する（Ｓ３０６）。本実施の形態においては、この所定の条件を第１可塑化能力と、第２可塑化能力との変化量が０．０１％以下であることとする。これは、可塑化能力によってメルトプール５３が生成する位置が異なるために実施している処理であり、その変化量が極めて小さくなった場合に定常状態になったと判断できる。なお、変化量は、この値に限定するものではなく、可塑化能力に変化が見られないと判断可能な値であればよく、例えば０．１や０．００５等でもよいが、精度の観点からいえばより小さい値の方が好ましい。

第２可塑化能力が所定の条件を満たさないと判定された場合（Ｓ３０６，ＮＯ）、解析部１０２は、第２可塑化能力を第１可塑化能力として更新して新たな第１可塑化能力とし（Ｓ３０７）、ステップＳ３０３の各種物理量を算出する処理へ移行する。一方、第２可塑化能力が所定の条件を満たすと判定された場合（Ｓ３０６，ＹＥＳ）、解析部１０２は、第２可塑化能力に基づいて、上記ステップＳ３０３の各種物理量を算出する処理と同様の処理を行うことにより、各物理量を算出する（Ｓ３０８）。算出後、出力部１０４は、第２可塑化能力や、これに基づく各物理量をユーザの入力に応じて適宜記憶部１２、表示部１４、ＨＤＤ１５へ出力し（Ｓ３０９）、本フローは終了となる。例えば、出力部１０４は、各物理量をスクリュ３上の各計算領域（ポジション）毎に予め指定されたアウトプットファイルにアウトプットし、表示部１４上へ表示する。

次に、前述した可塑化能力算出処理の詳細を説明する。本実施の形態に係る可塑化能力算出処理は、一般的な非圧縮流動計算に、本実施の形態に係る流動制限を組み込んだものである。図１０は、その可塑化能力算出処理を示すフローチャートである。図１０に示されるように、先ず、解析部１０２は、流動計算を一般化した、保存則の基礎方程式である下記（１７）式および（１８）式に基づいて、仮の速度（メルトフィルム５２からメルトプール５３へ流れる溶融樹脂の流速）を算出する（Ｓ４０１）。算出後、解析部１０２は、算出した仮の速度の圧力補正を行うと共に、下記（１９）式に基づいて樹脂温度を計算する（Ｓ４０２）。ここで圧力補正された仮の速度を補正速度と称して以後説明を行う。

この（１７）式においては、左から、１項目が非定常項、２項目が対流項、３項目が拡散項、４項目が生成項を示している。また、（１８）式および（１９）式においては、ρ：密度、Ｃ_ｐ：比熱、Ｔ：温度、ｔ：時間、ｑ：熱流束、τ：応力、ｖ：速度、Ｐ：圧力である。

圧力補正後、判定部１０３は、流量制限を行う必要があるか否かを判定する（Ｓ４０３）。本実施の形態においては、この判定を、メルトフィルム５２とメルトプール５３の境界のメッシュにおける補正速度が、ステップＳ２の２次元流動解析処理で算出された最大流量に基づいて算出される最大速度（最大流速）を超えるか否かにより判定するようにしている。補正速度が最大速度を超える場合、メルトフィルム５２の溶融量よりメルトフィルム５２からメルトプール５３へ流出する流出量が大きいという矛盾した状態にあると判断され、流量制限を行う必要があると判定される。より具体的には、判定部１０３は、下記（２０）式を用いて最大流量から最大速度を算出し、最大速度＜補正速度であるか否か
を判定する。

この（２０）式においては、ｖ_ｍａｘ：最大速度、ｍ：最大流量、δ：メルトフィルム厚み、Δｚ：３次元メッシュの長さである。

流量制限を行う必要があると判定された場合（Ｓ４０３，ＹＥＳ）、解析部１０２は、補正速度を最大速度とし（Ｓ４０４）、上記（１８）式、（１９）式を用いた収束処理（マスバランス）を実行することにより真の速度を算出する（Ｓ４０５）。算出後、解析部１０２は、真の速度に基づいてメルトフィルム５２とメルトプール５３の境界のメッシュの体積を算出し、これに溶融体（溶融樹脂）の密度に乗じることにより第２可塑化能力を算出し（Ｓ４０６）、本フローは終了となる。一方、流量制限を行う必要がないと判定された場合（Ｓ４０３，ＮＯ）、補正速度をそのままに、真の速度を算出する処理に移行する。

以上に説明した可塑化解析処理により算出された解析結果について、その優位性を図面を用いて簡単に説明する。図１１は、実機での実験結果、可塑化解析処理により算出された算出結果、及びスクリュ特性式のみを用いた一般的な解析処理により算出された算出結果の各可塑化能力を示す図である。図１２は、各算出結果の可塑化能力と実機での実験結果の可塑化能力との一致度を示す図であり、図１３は、各算出結果の可塑化能力と実機での実験結果の可塑化能力との誤差の割合を示す図である。各図面に示されるＰＰは、ポリプロピレン、ＰＥはポリエチレン、ＰＭＭＡはポリメタクリル酸メチル樹脂、ＰＡ６は、ナイロン６、ＰＡ６６はナイロン６６を示している。実験結果は、実際に樹脂材料を溶融可塑化させた際の計測結果から導かれた定常状態の可塑化能力を示し、解析結果Ａは、本実施の形態に係る可塑化解析処理により導かれた定常状態の可塑化能力を示し、解析結果Ｂは従来の計算手法（スクリュ特性式を用いた手法）により導かれた定常状態の可塑化能力を示している。

図１１〜図１３に示されるように、何れのプラスチックにおいても、解析結果Ｂより解析結果Ａが実験結果に近い値を導き出しており、特にＰＡ６においてそれは顕著である。したがって、本実施の形態に係る可塑化解析処理によれば、従来と比較して非常に高い精度で可塑化能力を算出できることがわかる。

本実施の形態によれば、諸条件を反映した３次元の流動計算から予測精度の高い可塑化能力を算出することができる。特に、これに流量制限を組み込むことにより、メルトフィルム５２とメルトプール５３との境界にあるメッシュの流量が、最大流量以下になるように制限することができるため、当該メッシュの流量が最大流量を超えるといった現実に起こりえない事象を防止することができる。したがって、現実に近いより精度の高い可塑化能力を得られ、延いては、高精度の解析結果を得ることができる。更に、分割せん断発熱量に基づいて可塑化能力および各物理量を算出することにより、メルトフィルム５２内の溶融樹脂の温度上昇に用いられるせん断発熱量を適切な値にすることができ、上記の効果をより高めることが可能となる。これは特にスクリュ動力において顕著である。

なお、本実施の形態においては、スクリュ特性式を用いて初期可塑化能力を算出するようにしているが、可塑化能力のあたりをつけて、入力または格納された所定の値を初期可
塑化能力として取得するようにしてもよい。

本発明は、その要旨または主要な特徴から逸脱することなく、他の様々な形で実施することができる。そのため、前述の実施の形態は、あらゆる点で単なる例示に過ぎず、限定的に解釈してはならない。本発明の範囲は、特許請求の範囲によって示すものであって、明細書本文には、何ら拘束されない。更に、特許請求の範囲の均等範囲に属する全ての変形、様々な改良、代替および改質は、全て本発明の範囲内のものである。

また、実施の形態にて述べた可塑化シミュレーション装置１０における各種ステップを、可塑化シミュレーションプログラムとして、図１４に示されるような、コンピュータにより読み取り可能な可搬型の記録媒体８に記憶させ、当該記録媒体８を情報処理装置９に読み込ませることにより、前述した機能を情報処理装置９に実現させることができる。記録媒体８としては、例えば、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤディスク等）、磁気ディスク（ハードディスクドライブ等）、フラッシュメモリ、ＩＣカード、更にネットワークを介することで伝送可能な媒体等、コンピュータで読み取りや実行が可能な全ての媒体が含まれる。

なお、特許請求の範囲に記載の可塑化シミュレーション装置は、例えば、前述の実施の形態における可塑化シミュレーション装置１０である。情報取得部は、例えば取得部１０１であり、初期可塑化能力取得部、物理量算出部、最大流量算出部および溶融量算出部は、例えば解析部１０２であり、可塑化能力算出部は、例えば解析部１０２および判定部１０３である。また、流速判定部は、例えば判定部１０３である。

３スクリュ、４スクリュ式射出成形機、１０可塑化シミュレーション装置、１０１取得部、１０２解析部、１０３判定部。

Claims

スクリュ式の可塑化装置におけるシリンダ内の樹脂材料がソリッドベッド、メルトフィルム、およびメルトプールの何れか状態にある３ゾーンモデルに基づいて可塑化能力を算出する可塑化シミュレーション装置であって、
前記樹脂材料の樹脂物性と前記可塑化装置の運転条件と該可塑化装置の構成データとを含むパラメータを取得する情報取得部と、
初期可塑化能力を取得する初期可塑化能力取得部と、
前記初期可塑化能力、前記パラメータ、質量保存の式およびエネルギー保存の式を用いて固相率およびメルトフィルム厚みを含む物理量を算出する物理量算出部と、
前記パラメータおよび前記物理量に基づいて３次元流動計算を行い、可塑化能力を算出する可塑化能力算出部と
を備えることを特徴とする可塑化シミュレーション装置。
前記パラメータを用いて２次元流動計算を行い、前記樹脂材料の単位時間当たりの溶融量である最大流量を算出する最大流量算出部を更に備え、
前記可塑化能力算出部は、
前記物理量と前記パラメータとを用いて３次元流動計算を行い、所定のメッシュにおける前記樹脂材料の流速を算出する溶融量算出部と、
前記最大流量に基づいて、前記流速が適正な値であるか否かを判定する流速判定部と、
を有し、
前記流速が適正な値でないと判定された場合、前記パラメータおよび前記最大流量に基づいて可塑化能力を算出する
ことを特徴とする請求項１記載の可塑化シミュレーション装置。
前記流速判定部は、前記流速が前記最大流量に基づく最大流速を超えるか否かを判定し、前記流速が前記最大流速を超えると判定した場合、前記流速が適正な値でないと判定することを特徴とする請求項２記載の可塑化シミュレーション装置。
前記最大流量算出部は、前記パラメータを用いて２次元流動計算を行い、せん断速度分布を算出し、
前記可塑化能力算出部は、前記パラメータ、前記最大流量、および前記せん断速度分布に基づいて可塑化能力を算出することを特徴とする請求項２または請求項３記載の可塑化シミュレーション装置。
前記最大流量算出部は、前記パラメータを用いて２次元流動計算を行い、せん断速度分布を算出し、該せん断速度分布に基づいてメルトフィルムにおけるせん断発熱量を算出し、
前記可塑化能力算出部は、前記せん断発熱量に基づいて可塑化能力および各物理量を算出することを特徴とする請求項２〜請求項４の何れか１項に記載の可塑化シミュレーション装置。
スクリュ式の可塑化装置におけるシリンダ内の樹脂材料がソリッドベッド、メルトフィルム、およびメルトプールの何れか状態にある３ゾーンモデルに基づいて可塑化能力を算出する可塑化シミュレーション装置が実行する可塑化シミュレーション方法であって、
前記樹脂材料の樹脂物性と前記可塑化装置の運転条件と該可塑化装置の構成データとを含むパラメータを取得するステップと、
初期可塑化能力を取得するステップと、
前記初期可塑化能力、前記パラメータ、質量保存の式およびエネルギー保存の式を用いて固相率およびメルトフィルム厚みを含む物理量を算出するステップと、
前記パラメータおよび前記物理量に基づいて３次元流動計算を行い、可塑化能力を算出するするステップと
を有する可塑化シミュレーション方法。
スクリュ式の可塑化装置におけるシリンダ内の樹脂材料がソリッドベッド、メルトフィルム、およびメルトプールの何れか状態にある３ゾーンモデルに基づいて可塑化能力を算出することをコンピュータに実行させる可塑化シミュレーションプログラムであって、
前記コンピュータを、
前記樹脂材料の樹脂物性と前記可塑化装置の運転条件と該可塑化装置の構成データとを含むパラメータを取得する情報取得部と、
初期可塑化能力を取得する初期可塑化能力取得部と、
前記初期可塑化能力、前記パラメータ、質量保存の式およびエネルギー保存の式を用いて固相率およびメルトフィルム厚みを含む物理量を算出する物理量算出部と、
前記パラメータおよび前記物理量に基づいて３次元流動計算を行い、可塑化能力を算出する可塑化能力算出部
として機能させるための可塑化シミュレーションプログラム。