特開 2023-37317 ＰＶＴ特性計算モデル推定システムおよび方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2023037317

(43)【公開日】2023-03-15

(54)【発明の名称】ＰＶＴ特性計算モデル推定システムおよび方法

(51)【国際特許分類】

   B29C 45/76 20060101AFI20230308BHJP

   G06F 30/27 20200101ALN20230308BHJP

   G06F 113/22 20200101ALN20230308BHJP

   G06F 111/10 20200101ALN20230308BHJP

【ＦＩ】

B29C45/76

G06F30/27

G06F113:22

G06F111:10

【審査請求】未請求

【請求項の数】10

【出願形態】ＯＬ

(21)【出願番号】P 2021143980

(22)【出願日】2021-09-03

(71)【出願人】

【識別番号】000005108

【氏名又は名称】株式会社日立製作所

(74)【代理人】

【識別番号】110000279

【氏名又は名称】弁理士法人ウィルフォート国際特許事務所

(72)【発明者】

【氏名】小林 漢

(72)【発明者】

【氏名】八木 大介

(72)【発明者】

【氏名】中土 裕樹

【テーマコード（参考）】

4F206

5B146

【Ｆターム（参考）】

4F206AM23

4F206JP18

5B146AA06

5B146AA10

5B146BA04

5B146DC03

(57)【要約】

【課題】樹脂材料のＰＶＴ特性計算モデルを容易に推定できるようにする。
【解決手段】樹脂材料の樹脂材料計算モデルを推定するＰＶＴ特性計算モデル推定システム１において、制御装置１１と、制御装置１１により使用される記憶装置１３とを含み、演算装置１１は、前記第１樹脂材料で成形した際のプロセスデータと成形品品質と、少なくとも１つ以上の所定の関係式とに基づいて、前記第１樹脂材料計算モデルを推定し、前記所定の関係式は、前記第１樹脂材料と同一樹脂種である第２樹脂材料のＰＶＴ特性の実測により得られたＰＶＴ特性の計算モデルである第２樹脂材料計算モデルと、前記第２樹脂材料を成形した際のプロセスデータと成形品品質との関係を示す式であり、前記第２樹脂材料計算モデルと、前記第２樹脂材料を用いて成形した際のプロセスデータと成形品品質とが紐づけられて記憶装置１３に蓄積されたデータに基づいて生成される。
【選択図】図７

【特許請求の範囲】

【請求項1】

第１樹脂材料のＰＶＴ特性の計算モデルである第１樹脂材料計算モデルを推定するＰＶＴ特性計算モデル推定システムであって、
プロセッサと、前記プロセッサに接続される記憶装置とを含み、
前記プロセッサは、前記第１樹脂材料で成形した際のプロセスデータと成形品品質と、少なくとも１つ以上の所定の関係式とに基づいて、前記第１樹脂材料計算モデルを推定し、
前記所定の関係式は、
前記第１樹脂材料と同一樹脂種である第２樹脂材料のＰＶＴ特性の実測により得られたＰＶＴ特性の計算モデルである第２樹脂材料計算モデルと、前記第２樹脂材料を成形した際のプロセスデータと成形品品質との関係を示す式であり、前記第２樹脂材料計算モデルと、前記第２樹脂材料を用いて成形した際のプロセスデータと成形品品質とが紐づけられて前記記憶装置に蓄積されたデータに基づいて生成される、
ＰＶＴ特性計算モデル推定システム。

【請求項2】

請求項１に記載のＰＶＴ特性計算モデル推定システムであって、
前記第１樹脂材料と、前記第２樹脂材料とは、同一の型番の樹脂材料である
ＰＶＴ特性計算モデル推定システム。

【請求項3】

請求項１記載のＰＶＴ特性計算モデル推定システムであって、
前記所定の関係式は、
前記第２樹脂材料計算モデルの係数を出力パラメータとし、基準成形条件により前記第２樹脂材料で成形した際のプロセスデータと前記成形品品質とを入力パラメータとする回帰式であって、機械学習を用いて生成される、
ＰＶＴ特性計算モデル推定システム。

【請求項4】

請求項３記載のＰＶＴ特性計算モデル推定システムであって、
前記所定の関係式の生成に使用される前記プロセスデータは、前記第２樹脂材料の所定の単位毎に、各所定の単位間で共通の基準成形条件により成形した際の、金型内センサの測定値の特徴量に基づいて生成される、
ＰＶＴ特性計算モデル推定システム。

【請求項5】

請求項４記載のＰＶＴ特性計算モデル推定システムであって、
前記所定の単位は、材料サプライヤから供給されるリサイクル材の納品単位である、
ＰＶＴ特性計算モデル推定システム。

【請求項6】

請求項４又は請求項５に記載のＰＶＴ特性計算モデル推定システムであって、
前記特徴量には、圧力、温度のうち少なくともいずれか一つが含まれる、
ＰＶＴ特性計算モデル推定システム。

【請求項7】

請求項６記載のＰＶＴ特性計算モデル推定システムであって、
前記特徴量には、前記第２樹脂材料の圧力及び／又は温度の状態情報が含まれており、
前記第２樹脂材料の前記状態情報は、
前記金型内センサの測定値のピーク値、射出開始からピークまでの前記測定値の積分値、射出開始から型開きまでの前記測定値の積分値、前記測定値の最大微分値のうちの少なくともいずれか一つに基づいて計算される、
ＰＶＴ特性計算モデル推定システム。

【請求項8】

請求項３～請求項７のいずれか一項に記載のＰＶＴ特性計算モデル推定システムであって、
前記成形品品質は、成形品重量、成形品寸法のうち少なくともいずれか一つが含まれる、
ＰＶＴ特性計算モデル推定システム。

【請求項9】

請求項１～請求項８のいずれか一項に記載のＰＶＴ特性計算モデル推定システムであって、
前記プロセッサは、
複数のクライアント計算機からネットワークを介して発信される、前記ＰＶＴ特性が実測された前記第２樹脂材料のＰＶＴ特性の計算モデルと、前記第２樹脂材料を用いて成形した際のプロセスデータと成形品品質と、が紐づいたデータを前記記憶装置に保存し、前記記憶装置に保存したデータを基に前記所定の関係式を生成し、前記関係式を前記ネットワークを介して前記クライアント計算機に送信する、
ＰＶＴ特性計算モデル推定システム。

【請求項10】

第１樹脂材料のＰＶＴ特性の計算モデルである第１樹脂材料計算モデルを推定するＰＶＴ特性計算モデル推定システムよるＰＶＴ特性計算モデル推定方法であって、
ＰＶＴ特性計算モデル推定システムは、
前記第１樹脂材料で成形した際のプロセスデータと成形品品質と、少なくとも１つ以上の所定の関係式とに基づいて、前記第１樹脂材料計算モデルを推定し、
前記所定の関係式は、
前記第１樹脂材料と同一樹脂種である第２樹脂材料のＰＶＴ特性の実測により得られたＰＶＴ特性の計算モデルである第２樹脂材料計算モデルと、前記第２樹脂材料を成形した際のプロセスデータと成形品品質との関係を示す式であり、前記第２樹脂材料計算モデルと、前記第２樹脂材料を用いて成形した際のプロセスデータと成形品品質とが紐づけられて記憶装置に蓄積されたデータに基づいて生成される、
ＰＶＴ特性計算モデル推定方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、樹脂材料のＰＶＴ特性の計算モデルである樹脂材料計算モデルを推定する技術に関する。

【背景技術】

【0002】

射出成形プロセスにおいて、樹脂材料の物性ばらつきは成形品品質に影響を与える。特に比容積の圧力・温度依存性（ＰＶＴ特性）は成形品の寸法に大きく影響するため、ＰＶＴ特性が異なる材料を成形する際は成形条件の調整が必要となる。ここで成形条件の調整は容易な作業ではなく、熟練作業者の技術と時間を要する。このため、成形条件の調整を支援するために、射出成形ＣＡＥを用いたシミュレーションが活用されている。しかし、射出成形ＣＡＥを使用するためには、樹脂材料のＰＶＴ特性の測定・計算モデル構築が必要となり、専用の比容積測定装置を用いた複数回の試験を行わなければならない。このため、測定・モデル構築の低コスト化が課題である。この課題に対して、従来、専用装置を使わずに射出成形機および金型の情報からＰＶＴ特性を測定するための方法が検討されている（特許文献１）。

【0003】

特許文献１に記載の方法では、射出成形機から金型に樹脂を充填する経路を製品部とＰＶＴ特性を測定するためのプランジャユニット部で切り替えられるようにし、プランジャユニット部に樹脂を充填した際の樹脂温度、樹脂圧、プランジャ位置を検出し、ある樹脂温度・圧力における比容積を測定する。また、プランジャユニット部に供給した樹脂を供給側に逆送することにより、複数の樹脂温度・圧力条件下における比容積測定処理を繰り返して行い、ＰＶＴ特性を取得することができる。

【先行技術文献】

【特許文献】

【0004】

【特許文献1】特開２００６－１３７０５７号公報

【発明の概要】

【発明が解決しようとする課題】

【0005】

射出成形における成形品品質は、樹脂の材料特性の影響を受ける。例えば、同一種類（例えば、ポリプロピレン、略称：ＰＰ）の材料であっても、グレード単位やサプライヤ単位毎にＰＶＴ特性が異なるため、同じ成形条件で成形しても、金型内の収縮挙動が異なり、成形品品質が変動する。

【0006】

近年、プラスチック廃棄物による海洋汚染問題や、中国や東南アジアによるプラスチック廃棄物の輸入禁止措置を受けて、プラスチックのリサイクル材の活用にこれまで以上に大きな注目が集まっている。欧州を中心に一部の地域でバージン材の使用に対する課税や法規制も検討されており、プラスチックを用いた製品を製造するメーカにとってリサイクル材の活用は喫緊の課題となっている。しかし、リサイクル材は、成形時の熱履歴、使用時環境による劣化、およびリサイクル時の異物混入や熱履歴によって、バージン材と比較してＰＶＴ特性のばらつきが大きいため、バージン材以上に成形品品質のばらつきも大きくなる。

【0007】

図１は、後述する実施例の優位性を明らかにするための比較例を示すグラフであり、従来技術ではない。図１は、材料の種類は同じ（ＰＰ）かつ同じ型番であり、製造ロットが異なる３つのリサイクル材、ロットＡ、ロットＢ、ロットＣに対する、圧力水準が５０ＭＰａの場合の比容積の温度依存性を示している。図１において横軸は、樹脂温度を示し、縦軸は比容積を示している。ここで、ロットとは、サプライヤからの納品時期の違いを意味している。図１の各点プロットは、専用の比容積測定装置で得られた実測値であり、曲線プロットは、実測値に対して、数１に示す２－ｄｏｍａｉｎ Ｔａｉｔ ＰＶＴモデルとよばれる計算モデルをフィッティングした結果である。図１は、圧力一水準(５０ＭＰａ)における温度依存性のみを示しているが、圧力水準が変わった場合でもロット間における比容積の絶対値、温度依存変化率の大小関係は保たれている。図１より、同種の材料であってもロットが異なるとＰＶＴ特性が大きく変動することが確認できる。

【0008】

【数1】

【0009】

図２は、樹脂材料のＰＶＴ特性の測定に用いる、一般的なロードセル方式の測定装置の概要図である。この測定装置では、炉体を含む温度調節機構が備わった試験槽に被測定対象（試料）の樹脂を投入し、樹脂に対してピストンで圧力を加えた時の体積を測定することによって、ある温度・圧力値（温度及び圧力値）における比容積を測定することができる。

【0010】

次に、図３を用いて、一般的なＰＶＴ特性計算モデルの構築方法を説明する。図２に示す装置を用いて、任意の温度・圧力条件（温度及び圧力条件）下での比容積の測定を、所定の回数分異なる温度・圧力条件で繰り返し行う。その後、樹脂材料のＰＶＴ特性を表す数理モデルを選択し、比容積の測定点に対してフィッティングを施すことで、ＰＶＴ特性計算モデルを構築する。

【0011】

図４に構築したＰＶＴ特性計算モデルの一例を示す。図４に示すＰＶＴ特性計算モデルは、樹脂に与える圧力を固定した状態で温度を変化させて測定した比容積に対して、２－ｄｏｍａｉｎ Ｔａｉｔ ＰＶＴモデルでフィッティングした結果である。グラフより、数理モデルが実測比容積をよく表現していることが確認できる。

【0012】

特許文献１に記載の方法では、図２に示す測定装置を使わず、金型にＰＶＴ特性を取得するための専用ユニットを取り付けることが前提となっている。従って、この方法では、測定用の専用ユニットと、製品を得るための金型とを接続および切り替えるための追加開発が必要である。また、ＰＶＴ特性を得るためには、同一樹脂に対して複数の樹脂温度・圧力条件下で測定を行う必要がある。ここで、樹脂は与えられた温度や圧力の影響により、各種物性が変化する性質を持つ。この方法では、ＰＶＴ特性を得るために専用ユニットから樹脂を逆送し、同一樹脂を繰り返し使用しながら複数の温度・圧力条件下で測定を行うため、樹脂特性が劣化する虞がある。従って、この方法では、複数条件での測定を繰り返すための工数が必要であるとともに、被測定樹脂の劣化による測定精度が低下し、得られたＰＶＴ特性を用いた成形では、成形品品質が低下する可能性も考えられる。

【0013】

また、射出成形ＣＡＥと連携するためには、一つの樹脂材料のＰＶＴ特性を測定した後、ＰＶＴ特性を記述する計算モデルを用いて実測値に対してフィッティングを施し、材料の計算モデルを構築する工程も必要となる。

【0014】

本発明は、上記課題に鑑みなされたものであり、その目的は、樹脂材料のＰＶＴ特性計算モデルを容易に推定することができる技術を提供することにある。

【課題を解決するための手段】

【0015】

上記目的を達成するため、一観点に係るＰＶＴ特性計算モデル推定システムは、第１樹脂材料のＰＶＴ特性の計算モデルである第１樹脂材料計算モデルを推定するＰＶＴ特性計算モデル推定システムであって、プロセッサと、前記プロセッサにより使用されるメモリとを含み、前記プロセッサは、前記第１樹脂材料で成形した際のプロセスデータと成形品品質と、少なくとも１つ以上の所定の関係式とに基づいて、前記第１樹脂材料計算モデルを推定し、前記所定の関係式は、前記第１樹脂材料と同一樹脂種である第２樹脂材料のＰＶＴ特性の実測により得られたＰＶＴ特性の計算モデルである第２樹脂材料計算モデルと、前記第２樹脂材料を成形した際のプロセスデータと成形品品質との関係を示す式であり、前記第２樹脂材料計算モデルと、前記第２樹脂材料を用いて成形した際のプロセスデータと成形品品質とが紐づけられて前記メモリに蓄積されたデータに基づいて生成される。

【発明の効果】

【0016】

本発明によれば、樹脂材料のＰＶＴ特性計算モデルを容易に推定することができる。

【図面の簡単な説明】

【0017】

【図1】図１は、比較例としての、納品時期が異なるリサイクル材に対するＰＶＴ特性図である。

【図2】図２は、ＰＶＴ特性の測定装置の構成図である。

【図3】図３は、ＰＶＴ特性計算モデル構築方法を説明する図である。

【図4】図４は、ＰＶＴ特性計算モデルを示す図である。

【図5】図５は、実施例１に係るＰＶＴ特性計算モデル推定システムの機能ブロック図である。

【図6】図６は、実施例１に係る射出成形機の概念図である。

【図7】図７は、実施例１に係るＰＶＴ特性計算モデル推定システムの実現に使用することができる計算機の構成図である。

【図8】図８は、実施例１に係る射出成形におけるプロセスデータの取得に用いる金型の概念図である。

【図9】図９は、実施例１に係る特徴量を登録する処理のブロック図である。

【図10】図１０は、実施例１に係る射出成形時におけるセンサによる取得データの時系列変化を示すグラフである。

【図11】図１１は、実施例１に係る回帰モデルの学習用データセットの構成図である。

【図12】図１２は、実施例１に係る回帰モデル学習処理のフローチャートである。

【図13】図１３は、実施例１に係るＰＶＴ特性計算モデル係数推定処理のフローチャートである。

【図14】図１４は、実施例１に係るＰＶＴ特性計算モデル係数の推定用データセットの構成図である。

【図15】図１５は、実施例１に係る推定されたＰＶＴ特性計算モデルを説明する図である。

【図16】図１６は、実施例２に係るＰＶＴ特性計算モデル推定システムの機能ブロック図である。

【図17】図１７は、実施例１及び実施例２に係るＰＶＴ特性計算モデルデータベースの構成図である。

【発明を実施するための形態】

【0018】

実施形態について、図面を参照して説明する。なお、以下に説明する実施形態は特許請求の範囲に係る発明を限定するものではなく、また実施形態の中で説明されている諸要素及びその組み合わせの全てが発明の解決手段に必須であるとは限らない。

【0019】

本実施形態では、樹脂材料のＰＶＴ特性計算モデルを推定するシステムを提供する。このシステムは、例えば、ＰＶＴ特性が既知である第２樹脂材料のＰＶＴ特性計算モデルと、第２樹脂材料を用いて各所定の単位間で共通の基準射出成形条件で成形した際のセンサデータと成形品品質との関係を示す所定の関係式を、予め得られている、第２樹脂材料計算モデルと、第２樹脂材料を用いて成形した際のプロセスデータと成形品品質と、が紐づけられて蓄積されたデータに基づいて予め生成しておき、第１樹脂材料を用いて基準射出成形条件で成形した際のプロセスデータと、成形品品質と、前記生成した所定の関係式とから、前記第１樹脂材料のＰＶＴ特性計算モデルを推定する。

【0020】

本実施形態によれば、ＰＶＴ特性が既知である材料のＰＶＴ特性計算モデルと、前記ＰＶＴ特性が既知である樹脂材料を用いた成形から予め取得したセンサデータおよび成形品品質とから、少ない工数で容易に所定の樹脂材料のＰＶＴ特性計算モデルを得ることができる。これにより、例えばリサイクル材のような同一型番であっても納品単位（ロット）毎でのＰＶＴ特性ばらつきが大きいような材料に対して、納品単位毎のＰＶＴ特性計算モデルを得るための物性試験が不要となり、ＰＶＴ特性を基にした納品単位毎の射出成形条件の調整が可能になり、納品単位間のばらつき抑制およびリサイクル材の成形品品質の向上を実現できる。

【実施例0021】

図５は、実施例１に係るＰＶＴ特性計算モデル推定システムの機能ブロック図である。ＰＶＴ特性計算モデル推定システム１は、例えば製造実行システム２と、製造工場３と、学習及び計算モデル生成システム４を含む。以下に述べるＰＶＴ特性計算モデル推定システム１の各機能の一部または全部は、ソフトウェアとして構成してもよいし、ソフトウェアとハードウェアとの協働として実現してよいし、固定的な回路を有するハードウェアを用いて実現してもよい。製造実行システム２、および製造工場３の有する機能の少なくとも一部を、オペレータが実行してもよい。

【0022】

製造実行システム２は、製造工場３に対して生産実行を指示するシステムである。製造実行システム２は、予め設定しておいた基準となる射出成形条件（以下、基準成形条件とする。）に基づいた製造条件で生産する指示を製造工場３に送る。製造条件には、例えば、生産（射出成形）に用いる射出成形機を特定する情報、生産に使用する金型を特定する情報、生産に使用する材料を特定する情報、生産する成形品の数量、生産時期、要求品質などが含まれてもよい。

【0023】

製造実行システム２について詳細に説明する。製造実行システム２は、例えば、製造条件決定部２１と、製造実行指示部２２とを備える。

【0024】

製造条件決定部２１は、基準成形条件に基づいて、上述の製造条件を決定する機能である。製造条件決定部２１は、製造条件に関する情報を、製造実行指示部２２に渡す。なお、製造条件に関する情報は、製造工場３を介して学習及び計算モデル生成システム４に送信される。製造条件に関する情報は、金型、射出成形機および材料に関する所定の情報を含むことができる。

【0025】

所定の情報は、例えば、金型の容量、金型のランナー構成を含んでもよく、さらに、使用する樹脂材料のサプライヤ単位の製造番号（型番）およびサプライヤごとの納品番号（例えば、ロット番号）等を含んでいてもよい。

【0026】

製造実行指示部２２は、製造工場３に対して製造実行を指示する機能である。なお、製造実行を生産と呼ぶこともできる。製造実行指示には、例えば、基準成形条件での成形の要求が含まれる。本実施形態では、製造実行指示部２２は、製造条件決定部２１から送信された所定の情報を製造工場３に送信する。所定の情報は、製造工場３を介して学習及び計算モデル生成システム４に入力され、樹脂材料のサプライヤ単位の製造番号およびサプライヤごとの納品番号毎に、ＰＶＴ特性計算モデルを推定するために使用される。

【0027】

製造工場３を説明する。製造工場３は、製造実行システム２からの製造実行指示を受けて、射出成形プロセス３２を実行する。以下、射出成形を「ＩＭ」と略記する場合がある。

【0028】

製造工場３は、例えば、製造実行部３１と、複数台の射出成形機５０（図６で後述）と、複数台の金型５０９（図６で後述）と、成形品品質検査部３３とを有する。以下、成形品品質検査部３３を品質検査部３３と略記する場合がある。

【0029】

製造実行部３１は、製造実行指示部２２から基準成形条件での成形の指示を受け取ると、指示された金型と材料との組み合わせに対して、予め指定された基準成形条件を射出成形機５０に入力することにより射出成形プロセス３２を実行する。すなわち、射出成形プロセス３２は、基準成形条件に従って射出成形するプロセスである。

【0030】

品質検査部３３は、射出成形プロセス３２で得られた成形品の品質を検査する機能である。成形品品質は、例えば、寸法、重量に基づいて評価する。成形品の品質検査は、自動的に行われてもよいし、検査員により手動で行われてもよいし、半自動で行われてもよい。

【0031】

なお、品質検査部３３で検査する成形品品質を、成形品の反り量、バリ、傷、光沢、色彩に基づいて評価してもよい。品質検査部３３は、製造実行部３１が受け取った製造条件と、成形品品質の検査結果とを紐づけて、学習及び計算モデル生成システム４の成形品品質記録部４０４に出力する。

【0032】

なお、本実施例において、金型を使用した際の、使用材料ごとの材料特性に係わる情報は、予め製造工場３が保有する各射出成形機および金型に搭載されたセンサ３４により、射出成形機の所定の位置及び金型内の所定の位置おける物理量を測定して、学習及び計算モデル生成システム４のセンサ情報記録部４０１に出力される。金型に搭載されたセンサの詳細に関しては、図８を用いて後述する。ここで、材料特性に係わる情報とは、例えば、材料の流動性や材料物性、およびこれらに相関する物理量を意味する。

【0033】

ここで、射出成形機の所定の位置とは、例えば、ノズル先端部などである。金型内の所定の位置とは、例えば、金型の樹脂流入口などである。物理量には、例えば、樹脂の圧力、樹脂の温度、樹脂の速度、樹脂の材料物性、および金型の開き量（型開き量）が含まれる。材料物性とは、例えば、樹脂の密度、樹脂の粘度、樹脂の繊維長の分布（強化繊維含有量材料の場合）などである。このうち、材料の流動性に最も相関する物理量は、樹脂の粘度であるが、材料の流動性に相関する物理量としては、その他の圧力、温度、および速度から算出した流動性に相関する量を用いることもできる。

【0034】

学習及び計算モデル生成システム４は、例えば、センサ情報記録部４０１、特徴量抽出部４０２、特徴量データベース４０３、成形品品質記録部４０４、ＰＶＴ特性計算モデル係数記憶部４０５、連結処理部４０６、学習用データベース４０７、回帰モデル学習部４０８、回帰モデル保存部４０９、学習済みモデル記憶部４１０、回帰モデル読み出し部４１１、推定用データベース４１２、ＰＶＴ特性計算モデル推定部４１３、ＰＶＴ特性計算モデル保存部４１４、及びＰＶＴ特性計算モデルデータベース４１５を含む。

【0035】

一般に、計算モデルとは、ある現象を数学的に抽象化した数理モデルを指し、ＰＶＴ特性計算モデルとは、樹脂材料のＰＶＴ特性を表す数理モデルを指す。また、本実施例の説明においてＰＶＴ特性計算モデルは、所定の材料単位毎に数理モデルのパラメータ（係数）が決定された計算モデルのことを指す。

【0036】

センサ情報記録部４０１は、製造工場３の射出成形プロセス３２中においてセンサ３４により取得される、金型内の所定の位置における物理量を記録する。特徴量抽出部４０２は、センサ情報記録部４０１で一時記録された物理量に対して特徴量抽出を行い、所定の材料単位と、射出成形機と金型との組み合わせと、抽出した特徴量とを紐づけて特徴量データベース４０３に記録する。ここで、特徴量は物理量そのままでもよく、物理量及び特徴量は、プロセスデータに対応する。また、所定の材料単位とは、例えば、材料サプライヤごとの樹脂材料（リサイクル材も含む）の型番ごとの単位または、同一の型番のロットごとの単位等の材料を区別する単位を意味する。また、特徴量データベース４０３に記録された、材料単位と、射出成形機と金型との組み合わせと、抽出した特徴量とを紐づけたデータ群のことを、本実施例では、特徴量データセットという。特徴量抽出部４０２の処理に関しては、図６および図７を用いて後述する。

【0037】

ここで、射出成形機と金型との組み合わせを固定すると共に、射出成形プロセス３２における成形条件を固定した状態において、所定の材料単位だけが変化した場合、センサ３４から得られる物理量から抽出される特徴量は、材料単位間の材料情報（例えば、流動性や物性値）の変化の影響を強く受ける。

【0038】

このため、材料単位毎に固有の材料情報を特徴量間の差として記録することが可能となる。ここで、センサ３４から得られる物理量から抽出される特徴量は、射出成形機及び金型の機差による影響を受けるため、特徴量データセットは、射出成形機と金型の組み合わせごとに、材料単位毎に固有の材料情報を特徴量間の差として保存する。本実施例では、特徴量データセットは、全て、射出成形機と金型の組み合わせが固定されているものとして説明する。

【0039】

成形品品質記録部４０４は、品質検査部３３により検査された成形品品質と、製造条件（所定の材料単位と、射出成形機と金型との組み合わせ）と、を紐づけて記録する。

【0040】

ＰＶＴ特性計算モデル係数記憶部４０５には、所定の材料単位毎に予め測定された樹脂材料のＰＶＴ特性により構築されたＰＶＴ特性計算モデルの各係数が、所定の材料単位の情報と紐づいて記憶されている。本実施例において、ＰＶＴ特性計算モデルには、例えば２－ｄｏｍａｉｎ－ｔａｉｔモデル、Ｓｐｅｎｃｅｒ－Ｇｌｉｍｏｒｅモデル、Ｍｏｄｉｆｉｅｄ－Ｃｅｌｌモデル、Ｓｉｍｈａ－Ｓｏｍｃｙｎｓｋｙモデルなどの計算モデルを用いてよい。

【0041】

連結処理部４０６は、特徴量データベース４０３から取得した材料情報と紐づいた特徴量データセットと、成形品品質記録部４０４から取得した成形品品質及び製造条件と、ＰＶＴ特性計算モデル係数記憶部４０５から取得したＰＶＴ特性計算モデル（第２樹脂材料計算モデル）の各係数とを、材料情報を結合キーとして、特徴量データセットと、成形品品質と、ＰＶＴ特性計算モデルの各係数と、製造に使用した材料とが対応したデータセット（以下、このデータセットを学習用データセットという）を生成し、この学習用データセットを学習用データベース４０７に記録する。

【0042】

また、連結処理部４０６は、ＰＶＴ特性計算モデル係数記憶部４０５に係数が記憶されていない樹脂材料（例えば、ＰＶＴ特性が未知である樹脂材料：第１樹脂材料）の材料情報に対応する、特徴量データベース４０３から取得した材料情報と紐づいた特徴量データセットと、成形品品質記録部４０４から取得した成形品品質及び製造条件とを、材料情報を結合キーとして、特徴量データセットと、成形品品質と、製造に使用した材料とが対応したデータセット（以下、このデータセットを推定用データセットという）を生成し、この推定用データセットを推定用データベース４１２に記録する。この推定用データセットは、ＰＶＴ特性計算モデル係数を含んでいない。

【0043】

回帰モデル学習部４０８は、学習用データベース４０７から学習用データセットを取得し、特徴量データセットと成形品品質とを説明変数とし、ＰＶＴ特性計算モデル係数の内の１つを目的変数とした回帰モデルを用いて、説明変数から目的変数を予測する回帰モデルを学習し、学習済み回帰モデルを生成する。

【0044】

一般に回帰モデルとは、説明変数（Ｘ）から目的変数（ｙ）を予測するモデル（ｙ＝ｆ（Ｘ））のことを指し、回帰モデル内のパラメータは学習用データによって決定される。本実施例の説明において「回帰モデル」というときは、回帰モデル全般のことを指し、「学習済み回帰モデル」というときは、学習用データによって、モデルのパラメータが決定された回帰モデルのことを指す。

【0045】

本実施例において、回帰モデルには、例えば線形回帰、リッジ回帰、サポートベクタマシン、ニューラルネットワーク、ランダムフォレスト回帰等の回帰モデル、またはこれらを組み合わせた回帰モデルを用いてよい。また、使用する回帰モデルが、例えばニューラルネットワークのように複数の目的変数を持つことができる場合は、目的変数としてＰＶＴ特性計算モデルの係数の中から複数選定してもよい。

【0046】

回帰モデル保存部４０９は、回帰モデル学習部４０８において生成された学習済み回帰モデルを、学習済み回帰モデル記憶部４１０に記録する。回帰モデル読み出し部４１１は、学習済み回帰モデル記憶部４１０から学習済み回帰モデルを取得し、ＰＶＴ特性計算モデル推定部４１３へ入力する。

【0047】

推定用データベース４１２は、推定用データセットを記録する。

【0048】

ＰＶＴ特性計算モデル推定部４１３は、ＰＶＴ特性が未知である樹脂材料のＰＶＴ特性計算モデルを推定する。ＰＶＴ特性計算モデル推定部４１３は、推定用データベース４１２から読み出した推定用データセットと、回帰モデル読み出し部４１１から読み出した学習済み回帰モデル（例えば、同一の型番のＰＶＴ特性計算モデルが既知のロットの情報に基づいて学習された回帰モデル）とに基づいて、ＰＶＴ特性が未知である所定の材料単位の樹脂材料のＰＶＴ特性計算モデル（第１樹脂材料計算モデル）を推定し、推定したＰＶＴ特性計算モデル（推定ＰＶＴ特性計算モデル）をＰＶＴ特性計算モデル保存部４１４へ入力する。

【0049】

ＰＶＴ特性計算モデル保存部４１４は、入力された推定ＰＶＴ特性計算モデルをＰＶＴ特性計算モデルデータベース４１５に格納する。

【0050】

ＰＶＴ特性計算モデルデータベース４１５は、ＰＶＴ特性計算モデル係数記憶部４０５に格納されている既知のＰＶＴ特性計算モデルと、ＰＶＴ特性計算モデル保存部４１４により格納されたＰＶＴ特性計算モデル（推定ＰＶＴ特性計算モデル）とを、対応する材料単位と、ＰＶＴ特性計算モデルが既知であるか又は推定されたものであるかを示す情報とに紐づいたデータとして記録する。図１７は、本実施例のＰＶＴ特性計算モデルデータベースの構成を示す。図１７に示すように、ＰＶＴ特性計算モデルデータベース４１５は、所定の材料単位を示す樹脂材料の型番及びロットと、既知又は推定の情報と、ＰＶＴ特性計算モデルの係数とを格納する。

【0051】

図７は、ＰＶＴ特性計算モデル推定システム１の実現に使用することができる計算機１０の構成例を示す。ここでは、一つの計算機１０によりＰＶＴ特性計算モデル推定システム１を実現する場合を説明するが、これに限らず、複数の計算機を連携させることにより一つまたは複数のＰＶＴ特性計算モデル推定システム１を構築することもできる。また、上述の通り、製造実行システム２、および製造工場３は、専用のソフトウェアやハードウェアを用いず、各機能の一部または全部をオペレータが実施することで、ＰＶＴ特性計算モデル推定システム１を実現することもできる。

【0052】

計算機１０は、例えば、演算装置１１、メモリ１２、記憶装置１３、入力装置１４、出力装置１５、通信装置１６、媒体インターフェース部１７を備えており、それら各装置１１～１７は通信経路ＣＮ１により接続されている。通信経路ＣＮ１は、例えば、内部バス、ＬＡＮ（Ｌｏｃａｌ Ａｒｅａ Ｎｅｔｗｏｒｋ）などである。

【0053】

演算装置１１は、例えばマイクロプロセッサなどから構成されており、プロセッサの一例である。演算装置１１は、マイクロプロセッサに限らず、例えば、ＤＳＰ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ）、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）などを含んでもよい。演算装置１１は、記憶装置１３に記録されたコンピュータプログラムをメモリ１２に読み出して実行することにより、ＰＶＴ特性計算モデル推定システム１としての各機能２１、２２、３１、４０１～４１５、６０を実現する。

【0054】

記憶装置１３は、コンピュータプログラムとデータとを記憶する装置であり、例えば、フラッシュメモリまたはハードディスクなどの書き換え可能な記憶媒体を有する。記憶装置１３には、オペレータにＧＵＩ（Ｇｒａｐｈｉｃａｌ Ｕｓｅｒ Ｉｎｔｅｒｆａｃｅ）を提供するＧＵＩ部６０を実現するためのコンピュータプログラムと、上述した各機能２１、２２、３１、４０１～４１５を実現するためのコンピュータプログラムとが格納される。

【0055】

入力装置１４は、オペレータが計算機１０に情報を入力する装置である。入力装置１４としては、例えば、キーボード、タッチパネル、マウスなどのポインティングデバイス、音声指示装置（いずれも不図示）などがある。出力装置１５は、計算機１０が情報を出力する装置である。出力装置１５としては、例えば、ディスプレイ、プリンタ、音声合成装置（いずれも不図示）などがある。

【0056】

通信装置１６は、外部の情報処理装置と計算機１０とを通信経路ＣＮ２を介して通信させる装置である。外部の情報処理装置としては、図示せぬ計算機のほかに、外部記憶装置１９がある。計算機１０は外部記憶装置１９に記録されたデータ（計算機固有情報、生産実績など）およびコンピュータプログラムを読み込むことができる。計算機１０は、記憶装置１３に記憶されたコンピュータプログラムおよびデータの全部または一部を、外部記憶装置１９に送信して記憶させることもできる。

【0057】

媒体インターフェース部１７は、外部記録媒体１８に読み書きする装置である。外部記録媒体１８としては、例えば、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ Ｓｅｒｉａｌ Ｂｕｓ）メモリ、メモリカード、ハードディスクなどがある。外部記録媒体１８から記憶装置１３に記憶されたコンピュータプログラムおよびデータの全部または一部を外部記録媒体１８に転送して記憶させることもできる。

【0058】

図６は、射出成形機５０の概念図である。図６を用いて、射出成形プロセスの各過程を説明する。本実施例において、成形現象とは、射出成形プロセスにおいて生じる一連の現象を示す。本実施例では、射出成形プロセスを、計量および可塑化過程と、射出および保圧過程と、冷却過程と、取出過程とに大別する。

【0059】

計量および可塑化過程では、可塑化用モータ５０１を駆動力としてスクリュー５０２を後退させ、ホッパー５０３から樹脂ペレット５０４をシリンダ５０５内へ供給する。そして、ヒータ５０６による加熱とスクリュー５０２の回転とにより、樹脂を可塑化させて均一な溶融状態とする。スクリュー５０２の背圧および回転数の設定により、溶融樹脂の密度と強化繊維の破断度合いとが変化し、これらの変化は成形品品質に影響を与える。

【0060】

射出および保圧過程では、射出用モータ５０７を駆動力としてスクリュー５０２を前進させ、ノズル５０８を介して溶融樹脂を金型５０９内へ射出する。金型５０９内に射出された溶融樹脂には、金型５０９の壁面からの冷却と、流動に起因するせん断発熱とが並行して作用する。すなわち溶融樹脂は、冷却作用と加熱作用を受けながら金型５０９内を流動する。金型５０９を閉じておく力である型締め力が小さい場合は、溶融樹脂の固化後に微小な金型開きが生じてしまい、その微小な隙間により成形品品質が影響を受ける。

【0061】

冷却過程では、一定温度に保持された金型５０９により、溶融樹脂が固化温度以下に冷却される。この冷却過程において発生する残留応力は、成形品の品質に影響を与える。残留応力は、金型内での流動により生じる材料物性の異方性、保圧による密度分布、成形収縮率の不均等に伴い発生する。

【0062】

取出過程では、金型５０９を開閉するモータ５１１を駆動力として型締機構５１２を駆動させることにより、金型５０９を開く。そして、突き出し用モータ５１３を駆動力としてエジェクタ機構５１４を駆動させることにより、固化した成形品を金型５０９から取り出す場合において、十分な突き出し力が成形品に均等に作用しなかったときには、成形品に残留応力が残ってしまい、成形品の品質に影響する。

【0063】

射出成形機５０において、ロードセル５１０による圧力値が、入力された成形条件内の圧力値へ近づくように圧力制御される。シリンダ５０５の温度は、複数のヒータ５０６により制御される。スクリュー５０２の形状とシリンダ５０５の形状とノズル５０８の形状とによって、射出成形機毎に異なる圧力損失が生じる。これにより、金型５０９の樹脂流入口における圧力は、射出成形機５０に入力された成形条件に示される圧力よりも低い値となる。さらに、ヒータ５０６の配置とノズル部における樹脂のせん断発熱とに起因して、金型５０９の樹脂流入口における樹脂温度は、射出成形機５０に入力された成形条件に示される樹脂温度と異なる場合がある。

【0064】

射出機構の構成（スクリュー５０２の形状、シリンダ５０５の形状、ノズル５０８の形状、ヒータ５０６の配置等）は、射出成形機ごとに異なり、これが機差となり成形品品質に影響を与える場合がある。

【0065】

成形品品質は、形状特性（重量、寸法、厚さ、ヒケ、バリ、反りなど）と、外観不良などの表面特性（ウェルド、シルバー、焼け、白化、傷、気泡、剥離、フローマーク、ジェッティング、色・光沢など）と、機械的・光学特性（引張強度、耐衝撃特性、透過率など）で評価される。

【0066】

本実施例においては、評価する成形品の品質として、例えば、重量、寸法とを用いている。成形品形状にショートショットやオーバーパック等の不良が無い場合、成形品の重量は金型容積に対する樹脂材料の比容積によって変化すると考えられるため、重量はＰＶＴ特性と相関の大きい物理量である。また、同じく上述の不良が無い場合、成形品の寸法は樹脂の収縮率によって変化すると考えらえるため、寸法もＰＶＴ特性と相関の大きい物理量である。

【0067】

形状特性は、射出および保圧過程と冷却過程とにおける、圧力および温度の履歴と型締力とに強い相関がある。表面特性は、発生する現象に対してそれぞれ発生要因が異なるが、例えばフローマークおよびジェッティングは、射出過程における樹脂の温度と速度に強い相関がある。機械的および光学的特性は、例えば引張強度の場合、破壊試験での評価が必要となるため、重量などの相関する他の品質指標で評価されることが多い。

【0068】

成形条件には、射出成形プロセスの各過程に対応したパラメータが設定される。計量および可塑化過程については、計量位置、サックバック、背圧、背圧速度、および回転数などが設定される。射出および保圧過程については、圧力と温度と時間と速度とがそれぞれ設定される。射出および保圧過程については、射出と圧力とを切り替えるスクリュー位置（ＶＰ切替位置）と、金型５０９の型締め力も設定される。冷却過程については、保圧後の冷却時間が設定される。温度に関するパラメータとして、複数のヒータ５０６の温度、および金型５０９を冷却するための冷媒の温度および流量などが設定される。

【0069】

図８は、本実施例において、射出成形におけるプロセスデータの取得に用いる金型の概要を示す。図８において、図８（Ａ）は、金型７０の製品部の上面図を示し、図８（Ｂ）は、金型７０の製品部の側面図を示し、図８（Ｃ）は、金型のランナー部７５の上面図を示す。この金型７０は、ランナー部７５から５点のピンゲート方式で製品部へ樹脂が流入する構造である。実際の成形実験では、製品部のセンサ配置部７１に、金型内の樹脂の温度を測定するための樹脂用温度センサ７２と、樹脂の圧力を測定する圧力センサ７３と、金型の温度を測定する金型用温度センサ７４とを配置し、ランナー部７５のセンサ配置部７６に、樹脂用温度センサ７２と、圧力センサ７３とを配置し、これらの時間変化を取得した。これら樹脂用温度センサ７２と、圧力センサ７３と、金型用温度センサ７４とが、センサ３４に対応する。成形に使用する樹脂材料は、ポリプロピレン（ＰＰ）を用いた。射出成形機は、最大型締力１５０ｔおよびスクリュー径４４ｍｍの電動射出成形機を用いた。

【0070】

図９は、センサ３４から得られた物理量に対して特徴量抽出部４０２が特徴量抽出を行い、その結果を特徴量データベース４０３へ記録することで、材料情報を取得する方法の例を示すブロック図である。図９に示す材料情報の取得情報は、所定の位置に所定の物理量を計測するセンサが設けられた「センサ付き金型」、「センサ内臓金型」のいずれかを用いることにより実現される。

【0071】

まず任意の材料６０１に対して、固定された成形条件である基準成形条件６０２を、実際の射出成形機６０３へ入力することにより、金型内の所定部位における物理量を取得する。ここで、射出成形機６０３は、図６で述べた射出成形機５０に対応する。また、基準成形条件は、図５で述べた射出成形プロセス３２を実行する際に射出成形機に入力される条件に対応する。

【0072】

金型内の所定部位における物理量は、使用する材料自体の材料情報、金型および射出成形機固有の機差、成形条件の影響をうける。このため、金型と射出成形機の組み合わせごとに基準成形条件６０２を設定することにより、機差および成形条件の影響を抑制して、材料固有の材料情報を物理量の特徴量として特徴量データベース６１０に記録できる。なお、基準成形条件は、金型と射出成形機の組み合わせごとに変更してもよい。

【0073】

実際の射出成形機６０３において成形現象を取得するためには、金型内センサ６０６を用いる。金型６０４内の任意の位置に金型内センサ６０６を配置することにより、金型６０４内の成形現象を直接測定して、材料情報と相関した物理量の実測値６０８を取得することができる。成形品６０５の品質は、製品品質検査６０７により取得され、成形品品質記録部６１１に記録される。

【0074】

得られた物理量の実測値から、特徴量が抽出される（６０９）。得られた物理量の実測値は、いずれも射出成形プロセス中の時間変化として取得されるため、直接評価することは難しい。そこで、本実施例では、物理量の時間変化から、材料情報と相関する特徴量を取得することにより、材料情報の定量的評価した。材料間で共通した基準成形条件で成形することによって、材料間の特徴量を比較することで、材料間の材料情報が比較可能となる。

【0075】

図１０および図１１を用いて、図４で述べたＰＶＴ特性を得るための実験例の測定結果を説明する。図１０は、納品時期の異なる３つのＰＰの材料ロットＡ、Ｂ、Ｃそれぞれに対して、基準成形条件で射出成形した際のランナーのセンサ配置部７６における、圧力センサ７３および樹脂用温度センサ７２の時系列データを示す。図１０に示すように、センサ取得波形は、同一成形条件であっても、ロット毎に時系列データの変動が異なっており、圧力センサ７３の時系列データが、材料固有の材料情報、例えば図１に示す所定の材料単位間のＰＶＴ特性のばらつきの影響を受けていることが確認できる。

【0076】

図１１には、図１０に示す材料ロット毎の金型内の異なる位置に取り付けた複数の圧力センサおよび樹脂温度センサ波形から抽出された特徴量の結果例が含まれている。図１１における説明変数の一部分が特徴量の結果例である。特徴量としては、圧力センサおよび樹脂温度センサのセンサ値のピーク値、ピーク値までの積分値、ピーク値からの積分値等がある。なお、図示していないが、圧力センサおよび樹脂温度センサのセンサ値の最大微分値、ピーク値および位置に設置したセンサ値を特徴量としてもよい。

【0077】

図１１によれば、特徴量がロット間でばらついていることから、抽出した特徴量が材料情報の影響を受けていることが確認できる。また、同じ材料ロットでもセンサ取り付け位置によって特徴量が変化していることが確認できる。このように、ロット毎に特徴量抽出処理を行い、ロット毎に特徴量が紐づいたデータセットを特徴量データセットといい、特徴量データベース４０３に記録される。本実施例では、特徴量データセットは、図１１に示すように、圧力センサおよび樹脂用温度センサ、金型用温度センサの複数センサのそれぞれの特徴量を紐づけたデータセットとしている。

【0078】

以下、物理量を測定する金型の部位、樹脂情報に相関する物理量のパラメータおよび特徴量について説明する。

【0079】

まず、物理量を測定する金型内の部位（以下、測定部位）について説明する。いずれの金型構造においても、測定部位は、少なくとも金型内の樹脂流入口からキャビティ内に至るまでのスプルー部あるいはランナー部を含むことが好ましい。

【0080】

スプルー部およびランナー部に加えて、例えば、キャビティ内のゲート直下部、樹脂合流部（ウェルド部）、流動末端部などのように、特徴的な流動が観測されうる部位を測定部位として使用してもよい。この場合、取り付け位置が異なる複数のセンサにより得られる物理量から、より高精度に材料固有のＰＶＴ特性に相関する物理量を求めることができる。

【0081】

例えば、基準成形条件での成形において、成形中の樹脂の圧力・温度は金型内の位置によって異なる。そのため、取り付け位置が異なる複数のセンサによって、１つの成形条件で異なる圧力・温度条件下における樹脂の物理量を測定することが可能である。これにより、より高精度に材料固有のＰＶＴ特性に相関する物理量を求めることができる。

【0082】

なお、金型構造と測定する物理量とによって、適切な測定部位は異なる。金型開き量以外の物理量では、いずれの金型構造であっても、可能であるならスプルー部を測定部位とするのが好ましい。なお、本明細書において「好ましい」という表現は、何らかの有利な効果を期待できるという意味で使用しているにすぎず、その構成が必須であることを意味していない。

【0083】

サイドゲート、ジャンプゲート、サブマリンゲート、およびバナナゲートでは、スプルー部直下のランナー部や、ゲート直前のランナー部などにセンサ配置する。ピンゲートの場合、３プレート構造となるため、センサ配置には工夫が必要だが、スプルー部直下のランナー部などにセンサを配置する。ピンゲートの場合、キャビティには繋がらないダミーのランナーを測定用に設けて測定部位としてもよい。測定専用の部位を設けることにより、金型設計の自由度が向上する。ピンゲートの場合、キャビティには繋がらないダミーのランナーを測定用に設けて測定部位としてよい。測定専用の部位を設けることにより、金型設計の自由度が向上する。フィルムゲートやファンゲートの場合、ゲート部に流入する前のランナー部にセンサを設置する。

【0084】

上述の物理量として測定するパラメータについて説明する。本実施例では、所定の材料単位のＰＶＴ特性計算モデルの推定を行うために、少なくとも圧力および温度を測定する。圧力と温度の測定には、例えば、金型内圧力センサ、金型表面温度センサ、樹脂温度センサなどを用いることができる。樹脂温度センサには、熱電対などの接触式温度センサ、または赤外線放射温度計などの非接触式温度センサのいずれかまたは両方を用いることができる、圧力と温度とのいずれの物理量も、射出成形プロセス中の時間変化を記録する。

【0085】

上述の物理量の特徴量について説明する。本実施例においては、例えば、圧力の最大値（時間変化のピーク値）と積分値と、温度の最大値（時間変化のピーク値）とを特徴量として用いることができる。また、圧力の時間変化に対して、時間微分値の最大値を取得することも有効である。圧力の時間微分値の最大値は、材料の瞬間粘度と相関がある。また、圧力の積分値は、射出過程と保圧過程とで分けて算出してもよい。保圧過程における圧力の積分値はゲートシール後の無圧力下における比容積の温度依存変化と相関が大きい物理量であると考えられる。

【0086】

図１１を用いて、学習用データベース４０７に記録される学習用データセットを説明する。学習用データセットは、図１１に示すように、特徴量データセットと、成形品品質と、ＰＶＴ特性計算モデルの各係数（目的変数に対応）と、製造条件決定部２１から取得された製造条件と、所定の材料単位とが、紐づけられたデータである。学習用データセットは、学習用データベース４０７では、射出成形機と金型の組み合わせに対応付けられて、射出成形機と金型の組み合わせ毎に記録される。

【0087】

本実施例では、学習用データセットとしては、前述したＰＶＴ特性が既知である材料ロットＡ，Ｃを使用し、成形品品質として重量、寸法で評価し、ＰＶＴ特性計算モデルとして２－ｄｏｍａｉｎ－Ｔａｉｔモデルを用いた。

【0088】

図１２を用いて、学習及び計算モデル生成システム４における回帰モデル学習処理について説明する（回帰モデル学習フローチャート）。回帰モデル学習処理では、まず、回帰モデル学習部４０８は、学習用データベース４０７から学習用データセットを読み出し、学習用データセットの中の特徴量データセットと、成形品品質とについて標準化処理を行う（Ｓ１０１）。

【0089】

図１１は、回帰モデル学習処理で用いた射出成形機と金型の組み合わせ毎に記録される学習用データセットの一例である。学習用データセットは、材料ロットＡ、Ｃの各ロット名、各センサにおける圧力・樹脂温度のピーク値、ピーク値までの積分値、ピーク値からの積分値と成形品重量、寸法、および各ロットのＰＶＴ特性計算モデルの係数が紐づいている。

【0090】

標準化処理とは、対象とする学習用データセットの列ごとに、データ群の平均値および標準偏差を算出し、各列のデータに対して、その列の平均値を引き、その列の標準偏差で割る処理である。標準化処理により各列のデータ群は平均が０、標準偏差が１になるため、列間の単位の違いによる影響を除くことができる。一般に、標準化処理を行うことにより、回帰モデルの学習精度を向上させる効果がある。

【0091】

次に、学習及び計算モデル生成システム４の回帰モデル学習部４０８は、学習用データセットのＰＶＴ特性計算モデル係数を読み込む（Ｓ１０２）。次いで、回帰モデル学習部４０８は、目的変数（出力パラメータ）をＰＶＴ特性計算モデル係数の１つとし、説明変数（入力パラメータ）を、標準化処理後の特徴量データセット及び成形品品質とした、回帰モデルを作成して、この回帰モデルを機械学習により学習することにより学習済み回帰モデルを生成し、このような学習済み回帰モデルを、使用するＰＶＴ特性計算モデルの係数分繰り返して、各係数のそれぞれに対応する学習済み回帰モデルを生成する（Ｓ１０３）。これにより、使用するＰＶＴ特性計算モデルの係数分の学習済み回帰モデルが生成されることとなる。

【0092】

本実施例では、ＰＶＴ特性計算モデルとして、２－ｄｏｍａｉｎ－ｔａｉｔモデルを用いた例を示しているが、ＰＶＴ特性計算モデルには、例えば、Ｓｐｅｎｃｅｒ－Ｇｌｉｍｏｒｅモデル、Ｍｏｄｉｆｉｅｄ－Ｃｅｌｌモデル、Ｓｉｍｈａ－Ｓｏｍｃｙｎｓｋｙモデルなどの計算モデルを用いてよい。

【0093】

また、本実施例では、学習済み回帰モデルを例えば、リッジ回帰を用いて生成したが、線形回帰、サポートベクタマシン、ニューラルネットワーク、ランダムフォレスト回帰等の他の回帰モデル、またはこれらを組み合わせた回帰モデルを用いて学習済み回帰モデルを生成するようにしてもよい。また、使用する回帰モデルが、例えばニューラルネットワークのように複数の目的変数を持つことができる場合は、回帰モデルにおける目的変数として、ＰＶＴ特性計算モデルの複数の係数を同時に選定してもよい。

【0094】

次いで、学習及び計算モデル生成システム４は、標準化処理に係る情報および学習済み回帰モデルを学習済みモデル記憶部４１０に保存する（Ｓ１０４）。

【0095】

図１３、図１４を用いて、学習及び計算モデル生成システム４におけるＰＶＴ特性計算モデル係数推定処理を説明する（ＰＶＴ特性計算モデル係数推定フローチャート）。学習及び計算モデル生成システム４のＰＶＴ特性計算モデル推定部４１３は、推定用データベース４１２からＰＶＴ特性計算モデルを推定する対象の材料（基本的には、ＰＶＴ特性計算モデルの係数が未知である材料）の推定用データセット（例えば、特徴量データセット及び成形品品質）を読み込み、推定用データセットについて標準化処理を行う（Ｓ１０５）。標準化処理は、ステップＳ１０１で算出した各列の平均値、標準偏差を用いて同様な処理を行ってもよい。

【0096】

次に、学習及び計算モデル生成システム４のＰＶＴ特性計算モデル推定部４１３は、回帰モデル読み出し部４１１により学習済み回帰モデル記憶部４１０から対応する学習済み回帰モデルを読み出し（Ｓ１０６）、読み出した学習済み回帰モデルの説明変数（入力パラメータ）に対して標準化処理後の推定用データセットの値を入力することで、推定用データセットに対してＰＶＴ特性計算モデルの係数を生成する（Ｓ１０７）。

【0097】

図１４に推定用データセットの一例を示す。本実施例では、回帰モデルにリッジ回帰を用いたため、ＰＶＴ特性計算モデルの係数のそれぞれに対して、ステップＳ１０６及びＳ１０７の処理を繰り返し、推定する対象の材料のＰＶＴ特性計算モデルのすべての係数を生成した。このように生成した係数をＰＶＴ特性計算モデルの係数とすることにより、推定対象の材料のＰＶＴ特性計算モデル（推定ＰＶＴ特性計算モデル）が特定できる。

【0098】

次いで、学習及び計算モデル生成システム４のＰＶＴ特性計算モデル保存部４１４が、標準化処理に係る情報、推定用データセット、及び推定ＰＶＴ特性計算モデルを、ＰＶＴ特性計算モデルデータベース４１５に保存する（Ｓ１０８）。

【0099】

このようにして、本実施例では、推定対象の所定単位の樹脂材料について、この樹脂材料のＰＶＴ特性を細かく測定することなく、容易且つ迅速にＰＶＴ特性計算モデルを推定することができる。例えば、同一の型番におけるロット単位のＰＶＴ特性計算モデルを容易に推定することができる。

【0100】

図１５は、本実施例において、生成した推定ＰＶＴ特性計算モデルの妥当性評価結果を示す。本評価は、本実施例で生成した学習済み回帰モデルの有効性を確かめるために行ったものであり、ＰＶＴ特性の計算モデル推定システム１を実行するうえで必ずしも実施する必要はない。以下、本評価の方法について説明する。

【0101】

本評価では、推定したＰＶＴ特性計算モデルの妥当性を確認するために、ＰＶＴ特性計算モデルが既知である材料ロットＡ、Ｂ、Ｃにおける材料ロットＢをＰＶＴ特性モデルが未知である材料ロットと仮定した。具体的には、材料ロットＡ、Ｃの学習用データセットで学習済み回帰モデルを作成し、この学習済み回帰モデルと、材料ロットＢの推定用データセットを用いて、材料ロットＢのＰＶＴ特性計算モデルを推定し、推定したＰＶＴ特性計算モデルと、既知である材料ロットＢのＰＶＴ特性とを比較評価を行った。

【0102】

図１５は、推定された材料ロットＢのＰＶＴ特性計算モデルと、公知の方法（温度を３０℃から２５０℃まで５℃刻み、圧力水準５０、１００、１５０ＭＰａ）で測定された材料ロットＢのＰＶＴ特性との比較を示すグラフである。図１５のグラフは、横軸が温度であり、縦軸が比容積であり、図中の線プロットが本実施例の方法により推定した材料ロットＢのＰＶＴ特性計算モデルを示し、点プロットが公知の方法で測定された材料ロットＢのＰＶＴ特性（比容積の実測値）である。図１５の凡例に示す線・点の形状がそれぞれの圧力水準に対応する。

【0103】

図１５に示すグラフでは、点プロットに対して線プロットの一致度が高いほど、推定したＰＶＴ特性計算モデルの妥当性が高いことを意味する。実測値に対する計算モデルの妥当性を定量評価する代表的な指標としては、決定係数が使用され、決定係数が１に近づくほど、計算モデルの妥当性が高いと言える。図１５に示すグラフにおいては、測定時の圧力水準毎に決定係数を計算すると、いずれの圧力水準に対しても、測定したＰＶＴ特性と推定ＰＶＴ特性計算モデルとの決定係数は０．９以上であり、推定したＰＶＴ特性計算モデルの妥当性が高いことが確認できる。本評価によって、射出成形プロセス３２における基準成形条件による成形から得られる特徴量データおよび成形品品質を用いた学習済み回帰モデルより、ＰＶＴ特性が未知の材料のＰＶＴ特性計算モデルを推定する方法が、良好な方法であることが確認できる。

【実施例0104】

次に、実施例２に係るＰＶＴ特性計算モデル推定システムについて説明する。なお、実施例１に係るＰＶＴ特性計算モデル推定システムと同様な構成については、同一の符号を付す。実施例２に係るＰＶＴ特性計算モデル推定システムは、実施例１における計算機１０の機能を、インターネットを介して接続された１以上のクライアント８（クライアント計算機）とクラウドサーバ８０とで実現するようにしたシステムである。クライアント８は、データ送信部８０１と、取得データで記憶部８０２とを含み、クラウドサーバ８０上で構築された学習及び計算モデル生成システム４に接続可能である。

【0105】

クラウドサーバ８０には、学習及び計算モデル生成システム４が構築されている。クラウドサーバ８０は、ＰＶＴ特性計算モデルデータベース４１５から所定の材料単位のＰＶＴ特性計算モデルをクライアント８に提供する機能と、クライアント８から提供された成形品品質とセンサデータとを用いて、学習済みモデル記憶部４１０に記憶された学習済み回帰モデルから対象の材料単位についてのＰＶＴ特性計算モデルを推定し、クライアント８に推定ＰＶＴ特性計算モデルを提供する機能と、クライアント８から提供された成形品品質とセンサ情報とＰＶＴ特性計算モデルとを用いて、学習用データベース４０７の学習用データセットを更新し、学習用データセットにより新たに回帰モデルを学習することで学習済み回帰モデルを更新する機能と、を有する。

【0106】

以下、図１６と図１７を用いて、実施例２に係るＰＶＴ特性計算モデル推定システムのそれぞれの機能について説明する。

【0107】

図１６を用いて、ＰＶＴ特性計算モデルをクライアント８に提供する機能を説明する。
クライアント８は、クラウドサーバ８０からＰＶＴ特性計算モデルデータベース４１５の所定の材料単位毎のＰＶＴ特性計算モデルを取得可能であり、取得したＰＶＴ特性計算モデルを取得データ記憶部８０２に保存することができる。取得データ記憶部８０２に記憶されたＰＶＴ特性計算モデルは、例えば射出成形ＣＡＥ等に用いてもよい。

【0108】

図１６と図１３を用いて、クライアント８から提供された成形品品質とセンサデータとを用いて、学習済み回帰モデルからＰＶＴ特性計算モデルを推定し、クライアント８に推定ＰＶＴ特性計算モデルを提供する機能を説明する。クライアント８は、クラウドサーバ８０に対して、データ送信部８０１により、推定対象の材料の所定単位の情報と、対象の材料を用いて基準成形条件で成形した際の、成形品品質と、センサデータとを送信する。

【0109】

クラウドサーバ８０の学習及び計算モデル生成システム４は、データ送信部８０１から送信された、推定対象の材料の所定単位の情報と成形品品質とを成形品品質記録部４０４に格納し、センサデータをセンサ情報記録部４０１に渡し、ＰＶＴ特性計算モデル係数推定処理（Ｓ１０５～Ｓ１０８）を実行する。この時、射出成形機による機差を考慮するため、学習用データベース４０７における、クライアント８から送信された、基準成形条件の金型と射出成形機との組み合わせと一致する学習用データセットが使用される。

【0110】

ＰＶＴ特性計算モデル係数推定処理（Ｓ１０５～Ｓ１０８）によってＰＶＴ特性計算モデル保存部４１４に記録された推定ＰＶＴ特性計算モデルは、ＰＶＴ特性計算モデルデータベース４１５に格納される。クライアント８は、クラウドサーバ８０から、ＰＶＴ特性計算モデルデータベース４１５に格納されている、情報を送信した材料についての推定されたＰＶＴ特性計算モデルを取得して、取得データ記憶部８０２に保存することができる。

【0111】

図１６と図１２を用いて、クライアント８から提供された、成形品品質と、センサデータと、ＰＶＴ特性計算モデルとを用いて、クラウドサーバ８０上の学習用データベース４０７を更新し、新たに回帰モデルを学習することで学習済み回帰モデルを更新する機能を説明する。クライアント８は、クラウドサーバ８０の学習及び計算モデル生成システム４に、任意の所定単位の材料情報と、対象の材料を用いて基準成形条件で成形した際の、成形品品質と、センサ情報と、ＰＶＴ特性計算モデルとを送信する。

【0112】

学習及び計算モデル生成システム４は、データ送信部８０１から送信された、所定単位の材料情報及び成形品品質を成形品品質記録部４０４に格納し、センサデータをセンサ情報記録部４０１に格納し、ＰＶＴ特性計算モデルをＰＶＴ特性計算モデル係数記憶部４０５に格納し、学習済み回帰モデル構築処理（Ｓ１０１～Ｓ１０４）を実行する。

【0113】

この時、クラウドサーバ８０は、学習用データベース４０７に、クライアント８から送信された材料情報によって新たに構築された学習用データセットをマージして、学習用データベース４０７を更新してもよい。なお、クラウドサーバ８０は、クライアント８から送信された、基準成形条件の金型と射出成形機との組み合わせ毎に学習用データベース４０７を構築する。

【0114】

学習済み回帰モデル構築処理（Ｓ１０１～Ｓ１０４）において、更新された学習用データベース４０７を用いて生成された学習済み回帰モデルを、回帰モデル保存部４０９が学習済みモデル記憶部４１０に記録してもよい。また、回帰モデル読み出し部４１１は、更新前後の学習済み回帰モデルを選択する機能を有してもよい。

【0115】

上述のクライアント８から送信された材料情報を基に学習済み回帰モデルを更新する機能によって、１以上のクライアント８からのデータを用いて効率よく回帰モデルの学習に必要となる学習用データセットを拡充可能であるため、学習済み回帰モデルの構築・更新のための工数を低減することが可能となる。また、一般に、機械学習による回帰モデルの学習において、学習用データベースのサンプル数が多い程、回帰モデルの精度が良くなることが知られており、上述のように学習用データセットを拡充可能であるために、学習済み回帰モデルの精度を向上することができる。

【0116】

なお、本発明は、上述の実施例に限定されるものではなく、本発明の趣旨を逸脱しない範囲で、適宜変形して実施することが可能である。

【0117】

例えば、上記実施例では、同一型番の材料の各ロットを単位として、樹脂材料のＰＶＴ特性の実測により得られたＰＶＴ特性の計算モデルと、樹脂材料を成形した際のプロセスデータと、成形品品質との関係を示す式とを紐づけて記憶装置に格納し、同一型番の材料の所定のロットについてのＰＶＴ特性計算モデルを推定するようにしていたが、例えば、同一樹脂種の樹脂材料を単位として、樹脂材料のＰＶＴ特性の実測により得られたＰＶＴ特性の計算モデルと、樹脂材料を成形した際のプロセスデータと、成形品品質との関係を示す式とを紐づけて記憶装置に格納し、同一樹脂種の所定の樹脂材料についてのＰＶＴ特性計算モデルを推定するようにしてもよい。

【0118】

また、上記実施例において、プロセッサ等の演算装置が行っていた処理の一部又は全部を、専用のハードウェア回路で行うようにしてもよい。また、上記実施例におけるプログラムは、プログラムソースからインストールされてよい。プログラムソースは、プログラム配布サーバ又は記憶メディア（例えば可搬型の記憶メディア）であってもよい。