ホーム > 特許ランキング > ポリプラスチックス株式会社
知財求人 - 知財ポータルサイト「IP Force」
特開2024-46200樹脂成形品のボイドを予測する方法、樹脂成形品のボイドを低減する方法、樹脂成形品のボイドを予測するためのプログラム及び樹脂成形品のボイドを低減するためのプログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024046200
(43)【公開日】2024-04-03
(54)【発明の名称】樹脂成形品のボイドを予測する方法、樹脂成形品のボイドを低減する方法、樹脂成形品のボイドを予測するためのプログラム及び樹脂成形品のボイドを低減するためのプログラム
(51)【国際特許分類】
B29C 45/76 20060101AFI20240327BHJP
【ＦＩ】
B29C45/76
【審査請求】未請求
【請求項の数】12
【出願形態】ＯＬ
(21)【出願番号】P 2022151444
(22)【出願日】2022-09-22
(71)【出願人】
【識別番号】390006323
【氏名又は名称】ポリプラスチックス株式会社
(74)【代理人】
【識別番号】110002077
【氏名又は名称】園田・小林弁理士法人
(72)【発明者】
【氏名】▲浜▼野 裕輔
(72)【発明者】
【氏名】青木 現
(72)【発明者】
【氏名】宮崎 晃弘
(72)【発明者】
【氏名】天野 雄太
【テーマコード（参考）】
4F206
【Ｆターム（参考）】
4F206AM23
4F206JA07
4F206JL02
4F206JL09
4F206JQ90
(57)【要約】 （修正有）
【課題】熱可塑性樹脂を金型に射出成形してなる射出成形品の内部に生じるボイドの発生挙動を予測する方法を提供する。
【解決手段】前記射出成形品を複数の要素に分割した解析用モデルを作成することと、熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めることと、前記温度分布を用いて、離型後の冷却過程における温度分布を算出することと、前記温度分布及び前記圧力分布から、予め測定された前記熱可塑性樹脂の、弾性定数の温度依存性データ、熱膨張率のデータ及びＰＶＴデータを使用して、前記解析用モデルの弾性定数分布及び温度荷重分布を算出することと、前記離型後の冷却過程における前記弾性定数分布及び前記温度荷重分布を用いた構造解析によって、解析用モデルの各要素に発生するひずみを算出することと、前記ひずみからボイドの発生場所及び／又はボイド量を予測することとを含む方法。
【選択図】図２
【解決手段】前記射出成形品を複数の要素に分割した解析用モデルを作成することと、熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めることと、前記温度分布を用いて、離型後の冷却過程における温度分布を算出することと、前記温度分布及び前記圧力分布から、予め測定された前記熱可塑性樹脂の、弾性定数の温度依存性データ、熱膨張率のデータ及びＰＶＴデータを使用して、前記解析用モデルの弾性定数分布及び温度荷重分布を算出することと、前記離型後の冷却過程における前記弾性定数分布及び前記温度荷重分布を用いた構造解析によって、解析用モデルの各要素に発生するひずみを算出することと、前記ひずみからボイドの発生場所及び／又はボイド量を予測することとを含む方法。
【選択図】図２
【特許請求の範囲】
【請求項1】
熱可塑性樹脂を金型に射出成形してなる射出成形品の内部に生じるボイドの発生挙動を予測する方法であって、
前記射出成形品を複数の要素に分割した解析用モデルを作成するステップ（Ｓ１）と、
熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）と、
前記温度分布を用いて、離型後の冷却過程における前記解析用モデルの温度分布を算出するステップ（Ｓ３）と、
前記温度分布及び前記圧力分布から、予め測定された前記熱可塑性樹脂の、弾性定数の温度依存性データ、熱膨張率のデータ及びＰＶＴデータを使用して、前記解析用モデルの弾性定数分布及び温度荷重分布を算出するステップ（Ｓ４）と、
前記離型後の冷却過程における前記弾性定数分布及び前記温度荷重分布を用いた構造解析によって、前記解析用モデルの各要素に発生するひずみを算出するステップ（Ｓ５）と、
前記ひずみからボイドの発生場所及び／又はボイド量を予測するステップ（Ｓ６）と、を含む方法。
前記射出成形品を複数の要素に分割した解析用モデルを作成するステップ（Ｓ１）と、
熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）と、
前記温度分布を用いて、離型後の冷却過程における前記解析用モデルの温度分布を算出するステップ（Ｓ３）と、
前記温度分布及び前記圧力分布から、予め測定された前記熱可塑性樹脂の、弾性定数の温度依存性データ、熱膨張率のデータ及びＰＶＴデータを使用して、前記解析用モデルの弾性定数分布及び温度荷重分布を算出するステップ（Ｓ４）と、
前記離型後の冷却過程における前記弾性定数分布及び前記温度荷重分布を用いた構造解析によって、前記解析用モデルの各要素に発生するひずみを算出するステップ（Ｓ５）と、
前記ひずみからボイドの発生場所及び／又はボイド量を予測するステップ（Ｓ６）と、を含む方法。
【請求項2】
前記熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）は、前記熱可塑性樹脂を前記金型のゲートからキャビティへ射出してから離型する工程において前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）である、請求項１に記載の方法。
【請求項3】
前記解析用モデルの弾性定数分布及び温度荷重分布を算出するステップ（Ｓ４）において、前記弾性定数分布及び前記温度荷重分布はゲートシールを始点とし、射出開始から設定された時間が経過するのを終点として算出する、請求項１に記載の方法。
【請求項4】
前記ひずみを算出するステップ（Ｓ５）において、前記射出成形品が保圧冷却された離型時を始点とし、射出開始から設定された時間が経過するのを終点として前記温度分布から算出した温度依存性を考慮した、前記弾性定数分布及び前記温度荷重分布を用いて算出する、請求項１に記載の方法。
【請求項5】
前記ゲートシールが、前記解析用モデルについての流動解析において、前記射出成形品の重量が最大であって、かつ前記ゲートシールまでの時間が最小となるように、成形条件が設定されることによって決定される、請求項３に記載の方法。
【請求項6】
前記ゲートシールが、ゲート中心部の温度が流動停止温度Ｔｓに達する時点である、請求項３に記載の方法。
【請求項7】
前記流動停止温度Ｔｓは、前記ＰＶＴデータに対し、２－ｄｏｍａｉｎ Ｔａｉｔ ＰＶＴモデルのデータフィッティング係数ｂ５及びｂ６を用いて、圧力Ｐについて、Ｔｓ＝ｂ５＋ｂ６×Ｐとする、請求項６に記載の方法。
【請求項8】
前記流動停止温度Ｔｓは、比熱測定において１℃／ｍｉｎから５０℃／ｍｉｎにて冷却した際の変曲点とする、請求項６に記載の方法。
【請求項9】
前記熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）において、交流定常法（ＩＳＯ２２００７－６）と呼ばれる熱伝導率測定手法にて求めた熱伝導率を計算に用いる、請求項１に記載の方法。
【請求項10】
熱可塑性樹脂を金型に射出成形してなる射出成形品の内部に生じるボイドを低減する方法であって、請求項１に記載された方法によって予測されたボイド量に対し、設計、成形条件及び成形材料のうちの１つ以上を変化させたときの予測結果を比較して、ボイド量が所定の量以下に低減されるまでこれを繰り返す、方法。
【請求項11】
コンピュータに、請求項１に記載の方法を実行させるプログラム。
【請求項12】
コンピュータに、請求項１０に記載の方法を実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、樹脂成形品の成形不良発生の予測に関する。
【背景技術】
【0002】
熱可塑性樹脂を用いた複雑な形状の部品製造において射出成形が用いられる。その成形条件又は製品形状によっては、樹脂成形品に「ヒケ（成形品表面に生じる窪み）」又は「ボイド（成形品内部に生じる空洞）」と呼ばれる成形不良が生じることがある。
【0003】
ヒケ及びボイドは熱可塑性樹脂の射出成形において、溶融状態で射出された熱可塑性樹脂が冷却され固化する過程で生じる。特に結晶性樹脂の場合、金型充填直後はランダム状態だった分子鎖が、結晶化により配向（折りたたまれて整列）し、その結果、金型充填直後の体積（金型寸法）よりも体積が減少（収縮）することによってヒケ又はボイドが生じる。
【0004】
これらの成形不良が発生すると、製品の寸法精度低下（例えば、気密用途の部品では窪みによって相手側部材と接するシール面に隙間ができる）又は強度低下（ボイドを起点として破壊が生じやすくなる）が起こり得る。そこでヒケ及びボイドの抑制に関する技術向上が求められている。
【0005】
ヒケ及びボイドを抑制する対策としては、実際に成形した成形品を確認して成形条件を変更する、又は、成形品のゲート又は肉厚の設計変更等をする。しかし、それらに掛る時間及び費用が膨大になることから、近年では、流動解析ソフトウェアを用いた射出成形シミュレーションによりヒケ及び／又はボイドの発生を予測し、製品形状及び成形条件の適正化ができないか検討されてきた。
特許文献１では、流動解析ソフトにより得られたデータを構造解析ソフトによる歪み解析に応用し、ボイドの発生を予測する。
特許文献１では、流動解析ソフトにより得られたデータを構造解析ソフトによる歪み解析に応用し、ボイドの発生を予測する。
【先行技術文献】
【特許文献】
【0006】
【特許文献1】特開第２００９－２３３８８２号公報
【発明の概要】
【発明が解決しようとする課題】
【0007】
既存の手法では、樹脂材料による射出成形プロセスにおいて、冷却後の成形品のボイドの発生が十分な精度で予測できないことがあった。
本開示の目的は、樹脂材料による射出成形プロセスによる樹脂成形品のボイドの発生挙動を精度よく予測する方法を提供することである。この課題を解決することにより、ボイドの発生しない成形品を得るための製品形状設計、金型設計、成形条件設定、成形材料を設計段階で事前に想定することが可能となり、製品化を効率的に行うことができる。
本開示の目的は、樹脂材料による射出成形プロセスによる樹脂成形品のボイドの発生挙動を精度よく予測する方法を提供することである。この課題を解決することにより、ボイドの発生しない成形品を得るための製品形状設計、金型設計、成形条件設定、成形材料を設計段階で事前に想定することが可能となり、製品化を効率的に行うことができる。
【課題を解決するための手段】
【0008】
本開示は以下の態様を含む。
【0009】
［１］熱可塑性樹脂を金型に射出成形してなる射出成形品の内部に生じるボイドの発生挙動を予測する方法であって、
前記射出成形品を複数の要素に分割した解析用モデルを作成するステップ（Ｓ１）と、
熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）と、
前記温度分布及び前記圧力分布を用いて、離型後の冷却過程における前記解析用モデルの温度分布を算出するステップ（Ｓ３）と、
前記温度分布及び前記圧力分布から、予め測定された前記熱可塑性樹脂の、弾性定数の温度依存性データ、熱膨張率のデータ及びＰＶＴデータを使用して、前記解析用モデルの弾性定数分布及び温度荷重分布を算出するステップ（Ｓ４）と、
前記離型後の冷却過程における前記弾性定数分布及び前記温度荷重分布を用いた構造解析によって、前記解析用モデルの各要素に発生するひずみを算出するステップ（Ｓ５）と、
前記ひずみからボイドの発生場所及び／又はボイド量を予測するステップ（Ｓ６）と、を含む方法。
前記射出成形品を複数の要素に分割した解析用モデルを作成するステップ（Ｓ１）と、
熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）と、
前記温度分布及び前記圧力分布を用いて、離型後の冷却過程における前記解析用モデルの温度分布を算出するステップ（Ｓ３）と、
前記温度分布及び前記圧力分布から、予め測定された前記熱可塑性樹脂の、弾性定数の温度依存性データ、熱膨張率のデータ及びＰＶＴデータを使用して、前記解析用モデルの弾性定数分布及び温度荷重分布を算出するステップ（Ｓ４）と、
前記離型後の冷却過程における前記弾性定数分布及び前記温度荷重分布を用いた構造解析によって、前記解析用モデルの各要素に発生するひずみを算出するステップ（Ｓ５）と、
前記ひずみからボイドの発生場所及び／又はボイド量を予測するステップ（Ｓ６）と、を含む方法。
【0010】
［２］前記熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）は、前記熱可塑性樹脂を前記金型のゲートからキャビティへ射出してから離型する工程において前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）である、［１］に記載の方法。
【0011】
［３］前記解析用モデルの弾性定数分布及び温度荷重分布を算出するステップ（Ｓ４）において、前記弾性定数分布及び前記温度荷重分布はゲートシールを始点とし、射出開始から設定された時間が経過するのを終点として算出する、［１］又は［２］に記載の方法。
【0012】
［４］前記ひずみを算出するステップ（Ｓ５）において、前記射出成形品が保圧冷却された離型時を始点とし、射出開始から設定された時間が経過するのを終点として前記温度分布から算出した温度依存性を考慮した、前記弾性定数分布及び前記温度荷重分布を用いて算出する、［１］から［３］のいずれか一項に記載の方法。
【0013】
［５］前記ゲートシールが、前記解析用モデルについての流動解析において、前記射出成形品の重量が最大であって、かつ前記ゲートシールまでの時間が最小となるように、成形条件が設定されることによって決定される、［３］に記載の方法。
【0014】
［６］前記ゲートシールが、ゲート中心部の温度が流動停止温度Ｔｓに達する時点である、［３］に記載の方法。
【0015】
［７］前記流動停止温度Ｔｓは、前記ＰＶＴデータに対して、２－ｄｏｍａｉｎ Ｔａｉｔ ＰＶＴモデルのデータフィッティング係数ｂ５及びｂ６を用いて、圧力Ｐについて、Ｔｓ＝ｂ５＋ｂ６×Ｐとする、［６］に記載の方法。
【0016】
［８］前記流動停止温度Ｔｓは、比熱測定において１℃／ｍｉｎから５０℃／ｍｉｎにて冷却した際の変曲点とする、［６］に記載の方法。
【0017】
［９］前記熱可塑性樹脂を成形する工程における前記解析用モデルの温度分布及び圧力分布を求めるステップ（Ｓ２）において、交流定常法（ＩＳＯ２２００７－６）と呼ばれる熱伝導率測定手法にて求めた熱伝導率を計算に用いる、［１］から［８］のいずれか一項に記載の方法。
【0018】
［１０］熱可塑性樹脂を金型に射出成形してなる射出成形品の内部に生じるボイドを低減する方法であって、［１］から［９］のいずれか一項に記載された方法によって予測されたボイド量に対し、設計、成形条件及び成形材料のうちの１つ以上を変化させたときの予測結果を比較して、ボイド量が所定の量以下に低減されるまでこれを繰り返す、方法。
【0019】
［１１］コンピュータに、［１］から［９］のいずれか一項に記載の方法を実行させるプログラム。
【0020】
［１２］コンピュータに、［１０］に記載の方法を実行させるプログラム。
【図面の簡単な説明】
【0021】
【発明を実施するための形態】
【0022】
図１は、射出成形の各工程を説明するための図である。図１は一つの金型に注目して、成形の時系列を示している。
射出成形は、金型キャビティ内にゲートを通じて樹脂を射出することから開始する。射出直後は樹脂の送出速度が制御されているが、例えば予め設定された量の９９％の樹脂が充填された時に、樹脂圧力の制御（保圧）に切り替えられる（Ｖ－Ｐ切替）。切り替え後は保圧して射出が続けられる。
射出成形は、金型キャビティ内にゲートを通じて樹脂を射出することから開始する。射出直後は樹脂の送出速度が制御されているが、例えば予め設定された量の９９％の樹脂が充填された時に、樹脂圧力の制御（保圧）に切り替えられる（Ｖ－Ｐ切替）。切り替え後は保圧して射出が続けられる。
【0023】
ここで、安定した射出成形品を得るための指標の一つとしてゲートシールがある。ゲートシールとはゲート部の樹脂が固化し、流動停止する現象であり、このようにゲート部の樹脂が固化し、流動停止する時間はゲートシール時間と呼ばれる。ゲートシール時間前に保圧を停止すると、溶融樹脂がゲートを通じて射出成形機側に逆流するため充填不良や重量減少を生じる。一方ゲートシール時間以降に保圧を停止してもゲートは固化しているため樹脂が逆流することはなく、安定した射出成形品を得ることができるため、成形現場では必ず測定する指標である。
【0024】
ゲートシール時間は、成形品の重量が最大であって、かつゲートシールまでの時間が最小になるように成形条件を設定することによって決定してもよい。成形品の重量及びゲートシールまでの時間は実験（計量）によって取得してもよい。成形品の重量及びゲートシールまでの時間は実験に代えて流動解析によりシミュレーションで求めてもよい。
あるいは、ゲート中心部の温度が流動停止温度Ｔｓに達する時点をゲートシール時間としてもよい。ゲート中心部の温度が流動停止温度Ｔｓに達する時点は実験によって取得してもよいし、シミュレーションで求めてもよい。
流動停止温度Ｔｓは、比熱測定において１℃／ｍｉｎから５０℃／ｍｉｎにて冷却した際の変曲点としてもよい。
または、流動停止温度Ｔｓは、前記ＰＶＴデータに対し、２－ｄｏｍａｉｎ Ｔａｉｔ ＰＶＴモデルのデータフィッティング係数ｂ５及びｂ６を用いて、圧力Ｐについて、Ｔｓ＝ｂ５＋ｂ６×Ｐとしてもよい。ここで、２－ｄｏｍａｉｎ Ｔａｉｔ ＰＶＴモデルは、温度Ｔ及び圧力Ｐにおける比容積ｖ（Ｔ，Ｐ）が、定数Ｃと指数関数Ｂ（Ｔ）により、ｖ（Ｔ，Ｐ）＝ｖ０（Ｔ）［１－Ｃ×ｌｎ（１＋Ｐ／Ｂ（Ｔ））］＋ｖｔ（Ｔ，Ｐ）のように表されるとするものである。とくに、ｖ０（Ｔ）はＴの１次式であって、所定の温度Ｔｔに対してｖｔ（Ｔ，Ｐ）は高温側（Ｔ＞Ｔｔ）では零であるが、低温側（Ｔ＜Ｔｔ）ではＴとＰの指数関数になる。そして、このモデルにおいて、ＰＶＴデータからＴｔ＝ｂ５＋ｂ６×Ｐのように、Ｐの１次式でフィッティングされる。
なお、ｂ５はＰＶＴデータ（図４参照）の変曲点になる。ｂ５を境に傾きが変化しており、固相－液相転移が起こっている。
あるいは、ゲート中心部の温度が流動停止温度Ｔｓに達する時点をゲートシール時間としてもよい。ゲート中心部の温度が流動停止温度Ｔｓに達する時点は実験によって取得してもよいし、シミュレーションで求めてもよい。
流動停止温度Ｔｓは、比熱測定において１℃／ｍｉｎから５０℃／ｍｉｎにて冷却した際の変曲点としてもよい。
または、流動停止温度Ｔｓは、前記ＰＶＴデータに対し、２－ｄｏｍａｉｎ Ｔａｉｔ ＰＶＴモデルのデータフィッティング係数ｂ５及びｂ６を用いて、圧力Ｐについて、Ｔｓ＝ｂ５＋ｂ６×Ｐとしてもよい。ここで、２－ｄｏｍａｉｎ Ｔａｉｔ ＰＶＴモデルは、温度Ｔ及び圧力Ｐにおける比容積ｖ（Ｔ，Ｐ）が、定数Ｃと指数関数Ｂ（Ｔ）により、ｖ（Ｔ，Ｐ）＝ｖ０（Ｔ）［１－Ｃ×ｌｎ（１＋Ｐ／Ｂ（Ｔ））］＋ｖｔ（Ｔ，Ｐ）のように表されるとするものである。とくに、ｖ０（Ｔ）はＴの１次式であって、所定の温度Ｔｔに対してｖｔ（Ｔ，Ｐ）は高温側（Ｔ＞Ｔｔ）では零であるが、低温側（Ｔ＜Ｔｔ）ではＴとＰの指数関数になる。そして、このモデルにおいて、ＰＶＴデータからＴｔ＝ｂ５＋ｂ６×Ｐのように、Ｐの１次式でフィッティングされる。
なお、ｂ５はＰＶＴデータ（図４参照）の変曲点になる。ｂ５を境に傾きが変化しており、固相－液相転移が起こっている。
【0025】
保圧に引き続き、予め定められた冷却時間の間、金型内で射出成形品を冷却する（金型内冷却）。冷却時間を経過すると、型が開かれ、射出成形品を金型から突出して離型する。その後、金型は再び締められて金型キャビティを形成して、次の射出成形品を得るための成形に移る。突出された射出成形品を金型の外で冷却する（金型外冷却）ことができる。金型外冷却は、例えば射出成形品が室温になるまで続けることができる。
【0026】
実際の射出成形品においては、特に製品が肉厚の場合、ボイドは製品が金型外へ取り外された後（離型後）発生している場合がある。従来の流動解析ソフトでは金型内での温度と圧力しか計算ができないため、金型外における温度推移を計算して実際の現象を考慮することができない。
そのため、従来の流動解析ソフトを使用する手法では、ボイドの発生を精度よく予測できない場合があると考えられる。
さらに、金型外における温度推移を計算してそこからボイドの発生を予測する手法は知られていない。
そのため、従来の流動解析ソフトを使用する手法では、ボイドの発生を精度よく予測できない場合があると考えられる。
さらに、金型外における温度推移を計算してそこからボイドの発生を予測する手法は知られていない。
【0027】
本開示では、金型から離型した後の冷却を考慮した解析を行う。それにより得られた温度データから、弾性定数分布及び体積収縮率分布を得て、構造解析（ひずみ解析）に応用する連成解析により、実際の製品における結果と遜色ないレベルでボイドの発生場所と量の予測ができることが見出された。なお、本開示では金型の内外で、空冷によって冷却することを詳細に説明するが、冷却手段は空冷に限らない。
【0028】
以下、本開示の実施形態にかかるボイドの発生挙動を予測する方法について詳細に説明する。なお、本開示は以下の実施形態に限定されない。
（一実施形態）
（一実施形態）
【0029】
（ボイドの発生挙動を予測する方法）
本実施形態のボイドの発生挙動を予測する方法の一例について、図２を参照して詳細に説明する。この方法は、流動解析の結果を構造解析に応用する連成解析を行う。
図２のフローチャートに示すように、本実施形態のボイドの発生挙動を予測する方法の一例は、解析用モデルの作成（Ｓ１）、熱可塑性樹脂を射出して成形する工程（射出成形工程）における温度分布及び圧力分布の算出（Ｓ２）、離型後の冷却過程における温度分布の算出（金型外冷却解析）（Ｓ３）、弾性定数分布及び温度荷重分布の算出（Ｓ４）、解析用モデルの各要素に発生するひずみの算出（Ｓ５）、及び、ボイドの発生場所及び／又はボイド量の予測（Ｓ６）の各ステップを含む。
本実施形態のボイドの発生挙動を予測する方法の一例について、図２を参照して詳細に説明する。この方法は、流動解析の結果を構造解析に応用する連成解析を行う。
図２のフローチャートに示すように、本実施形態のボイドの発生挙動を予測する方法の一例は、解析用モデルの作成（Ｓ１）、熱可塑性樹脂を射出して成形する工程（射出成形工程）における温度分布及び圧力分布の算出（Ｓ２）、離型後の冷却過程における温度分布の算出（金型外冷却解析）（Ｓ３）、弾性定数分布及び温度荷重分布の算出（Ｓ４）、解析用モデルの各要素に発生するひずみの算出（Ｓ５）、及び、ボイドの発生場所及び／又はボイド量の予測（Ｓ６）の各ステップを含む。
【0030】
（解析用モデルの作成（Ｓ１））
射出成形品を複数の要素に分割した解析用モデルを作成するステップ（Ｓ１）は射出成形品の形状を微小な要素に分割して、シミュレーションの実行に必要なモデルを作成する（例として図６参照）。例えば、３次元形状測定、又は、ＣＡＤシステム等により、射出成形品の形状（これは射出成形品の設計上の形状又は金型の形状等とすることができる。金型の形状としては、ランナー及びゲートの位置、数、大きさなどの条件も含む）を計算機に取り込む。ついで、要素分割プリプロセッサ等で計算機に取り込んだ形状を複数の３次元要素に分割して、解析用モデルを作成する。
なお流動解析と構造解析の連成解析に際して、流動解析用モデルとそれとは別の構造解析用モデルを用意して解析を実行することもあり得る。しかし、本実施形態では流動解析から構造解析へ同じ解析用モデルを引き継いで解析する例を詳しく説明する。
射出成形品を複数の要素に分割した解析用モデルを作成するステップ（Ｓ１）は射出成形品の形状を微小な要素に分割して、シミュレーションの実行に必要なモデルを作成する（例として図６参照）。例えば、３次元形状測定、又は、ＣＡＤシステム等により、射出成形品の形状（これは射出成形品の設計上の形状又は金型の形状等とすることができる。金型の形状としては、ランナー及びゲートの位置、数、大きさなどの条件も含む）を計算機に取り込む。ついで、要素分割プリプロセッサ等で計算機に取り込んだ形状を複数の３次元要素に分割して、解析用モデルを作成する。
なお流動解析と構造解析の連成解析に際して、流動解析用モデルとそれとは別の構造解析用モデルを用意して解析を実行することもあり得る。しかし、本実施形態では流動解析から構造解析へ同じ解析用モデルを引き継いで解析する例を詳しく説明する。
【0031】
（射出成形工程における温度分布及び圧力分布の算出（Ｓ２））
熱可塑性樹脂を射出して成形する工程（射出成形工程）における解析用モデルの温度分布及び圧力分布を算出するステップ（Ｓ２）では、作成した解析用モデルを使用して、射出成形の流動解析（シミュレーション）を実行する。流動解析によって、熱可塑性樹脂を射出して成形する工程（金型内冷却工程も含みうる）における解析用モデルの各要素の温度及び圧力を算出する。これにより、射出成形工程における解析用モデルの温度分布及び圧力分布を算出する。
流動解析を実行する際には、射出成形に使用する熱可塑性樹脂についての物性値を入力する。流動解析に使用する物性値には、圧力、体積、及び温度の関係を表すデータ（以下、「ＰＶＴデータ」と称する。例として図４参照。）、熱伝導率データ、及び、比熱データ等がある。
なお、熱可塑性樹脂について、比熱は示差走査熱量計（ＤＳＣ）により測定することができる。熱可塑性樹脂の熱伝導率は交流定常法（ＩＳＯ ２２００７－６）にて求めた熱伝導率を用いることができる。
熱伝導率測定手法には、交流定常法、熱線法やホットディスク法などがある。交流定常法は良好な精度を与える。
熱可塑性樹脂を射出して成形する工程（射出成形工程）における解析用モデルの温度分布及び圧力分布を算出するステップ（Ｓ２）では、作成した解析用モデルを使用して、射出成形の流動解析（シミュレーション）を実行する。流動解析によって、熱可塑性樹脂を射出して成形する工程（金型内冷却工程も含みうる）における解析用モデルの各要素の温度及び圧力を算出する。これにより、射出成形工程における解析用モデルの温度分布及び圧力分布を算出する。
流動解析を実行する際には、射出成形に使用する熱可塑性樹脂についての物性値を入力する。流動解析に使用する物性値には、圧力、体積、及び温度の関係を表すデータ（以下、「ＰＶＴデータ」と称する。例として図４参照。）、熱伝導率データ、及び、比熱データ等がある。
なお、熱可塑性樹脂について、比熱は示差走査熱量計（ＤＳＣ）により測定することができる。熱可塑性樹脂の熱伝導率は交流定常法（ＩＳＯ ２２００７－６）にて求めた熱伝導率を用いることができる。
熱伝導率測定手法には、交流定常法、熱線法やホットディスク法などがある。交流定常法は良好な精度を与える。
【0032】
次に、熱可塑性樹脂の流動解析を行うための解析条件を入力する。成形条件としては、樹脂（シリンダー）温度、金型温度、射出速度、保圧力及び保圧時間などを含む。なお、金型温度は冷却水温度もしくはヒーター設定温度と同様である。
また、成形条件には流動解析の始点と終点の指定を含む。流動解析の始点は例えば金型のゲートからキャビティに熱可塑性樹脂を射出し始めた時点とすることができる。
流動解析の終点は、溶融した熱可塑性樹脂を金型キャビティへ射出後、冷却時間が終了し、離型前までとすることができる。特に離型の時点までとすることができる。流動解析における離型の時点は射出開始から設定されたある時間（実工程における離型までの時間との比較で決めてもよい）が経過した時点として指定し得る。
また、成形条件には流動解析の始点と終点の指定を含む。流動解析の始点は例えば金型のゲートからキャビティに熱可塑性樹脂を射出し始めた時点とすることができる。
流動解析の終点は、溶融した熱可塑性樹脂を金型キャビティへ射出後、冷却時間が終了し、離型前までとすることができる。特に離型の時点までとすることができる。流動解析における離型の時点は射出開始から設定されたある時間（実工程における離型までの時間との比較で決めてもよい）が経過した時点として指定し得る。
【0033】
（金型外冷却解析（Ｓ３））
金型外冷却解析のステップ（Ｓ３）では、算出した温度分布を用いて離型後の冷却過程における、解析用モデルの各要素の温度分布を算出する。このステップについては、構造解析ソフトを使用して解析することができる。
金型外冷却解析の始点は離型時となる。金型内での解析の場合、射出成形品に加えて、金型自体をモデル化し、成形品と金型との相互作用を考慮するため、金型の有限要素モデルを用いる必要がある。一方、金型外での解析の場合、金型自体のモデルは必要ないため、射出成形品のモデル自体及び境界条件を大きく変える必要があり、金型内での解析の延長としても、解析条件の設定が極めて困難になる。
金型外冷却解析の終点は、例えば射出成形品が金型外で室温になる時点とすることができる。「射出成形品が金型外で室温になる時点」については、例えば射出開始から実工程で計測して得られた射出成形品が金型外で室温となる時間で指定してもよい。あるいは、例えば金型外冷却解析の中で射出成形品の温度が室温になる時点を求めてもよい。一般的には射出開始から設定されたある時間を経過した時点を金型外冷却解析の終点とすることができる。
金型外冷却解析のステップ（Ｓ３）では、算出した温度分布を用いて離型後の冷却過程における、解析用モデルの各要素の温度分布を算出する。このステップについては、構造解析ソフトを使用して解析することができる。
金型外冷却解析の始点は離型時となる。金型内での解析の場合、射出成形品に加えて、金型自体をモデル化し、成形品と金型との相互作用を考慮するため、金型の有限要素モデルを用いる必要がある。一方、金型外での解析の場合、金型自体のモデルは必要ないため、射出成形品のモデル自体及び境界条件を大きく変える必要があり、金型内での解析の延長としても、解析条件の設定が極めて困難になる。
金型外冷却解析の終点は、例えば射出成形品が金型外で室温になる時点とすることができる。「射出成形品が金型外で室温になる時点」については、例えば射出開始から実工程で計測して得られた射出成形品が金型外で室温となる時間で指定してもよい。あるいは、例えば金型外冷却解析の中で射出成形品の温度が室温になる時点を求めてもよい。一般的には射出開始から設定されたある時間を経過した時点を金型外冷却解析の終点とすることができる。
【0034】
（弾性定数分布及び温度荷重分布の算出（Ｓ４））
弾性定数分布及び温度荷重分布を算出するステップ（Ｓ４）では、算出された温度分布（熱可塑性樹脂を成形する工程における温度分布及び離型後の冷却過程における温度分布）と圧力分布（熱可塑性樹脂を成形する工程における圧力分布）を用いて射出成形品の弾性定数分布と温度荷重分布（温度差の分布）を求める。特に温度分布及び圧力分布から、計算機のための変換プログラムによって、射出成形品の弾性定数分布と温度荷重分布を算出することができる。また、温度荷重分布とともに体積収縮率分布も算出し得る。
弾性定数分布及び温度荷重分布を算出するステップ（Ｓ４）では、算出された温度分布（熱可塑性樹脂を成形する工程における温度分布及び離型後の冷却過程における温度分布）と圧力分布（熱可塑性樹脂を成形する工程における圧力分布）を用いて射出成形品の弾性定数分布と温度荷重分布（温度差の分布）を求める。特に温度分布及び圧力分布から、計算機のための変換プログラムによって、射出成形品の弾性定数分布と温度荷重分布を算出することができる。また、温度荷重分布とともに体積収縮率分布も算出し得る。
【0035】
（弾性定数分布）
射出成形品の弾性定数分布は、予め取得した成形材料（熱可塑性樹脂）の弾性定数（一般には応力とひずみの比例係数）の温度依存性データを、流動解析により得た温度分布に当てはめることにより、求めることができる。
図３は、ある熱可塑性樹脂の弾性定数の温度依存性データを示している。横軸は温度であり、縦軸は弾性定数（ヤング率）である。横軸上の温度Ｔｃは、熱可塑性樹脂の転移温度を示す。転移温度とは、固化領域と溶融領域とに分ける温度のことであり、ノーフロー温度、固体化温度、固化温度、又は、固液転移温度とも呼ばれている。
図３のように、熱可塑性樹脂の弾性定数の温度依存性データが得られていれば、温度分布から弾性定数の分布を得ることができる。
なお、弾性定数は、単軸応力に対する弾性率（ヤング率）、ポアソン比及びせん断弾性率の３種類の値を含むことができる。これら３種類の値の弾性定数の使用により、より精度の高い解析が可能になる。
射出成形品の弾性定数分布は、予め取得した成形材料（熱可塑性樹脂）の弾性定数（一般には応力とひずみの比例係数）の温度依存性データを、流動解析により得た温度分布に当てはめることにより、求めることができる。
図３は、ある熱可塑性樹脂の弾性定数の温度依存性データを示している。横軸は温度であり、縦軸は弾性定数（ヤング率）である。横軸上の温度Ｔｃは、熱可塑性樹脂の転移温度を示す。転移温度とは、固化領域と溶融領域とに分ける温度のことであり、ノーフロー温度、固体化温度、固化温度、又は、固液転移温度とも呼ばれている。
図３のように、熱可塑性樹脂の弾性定数の温度依存性データが得られていれば、温度分布から弾性定数の分布を得ることができる。
なお、弾性定数は、単軸応力に対する弾性率（ヤング率）、ポアソン比及びせん断弾性率の３種類の値を含むことができる。これら３種類の値の弾性定数の使用により、より精度の高い解析が可能になる。
【0036】
（ポアソン比分布）
射出成形品のポアソン比分布は、予め取得した成形材料（熱可塑性樹脂）の弾性定数（ヤング率）の温度依存性データおよびＰＶＴデータから求めた体積弾性率から計算した値を、流動解析により得た温度分布に当てはめることにより、求めることができる。
射出成形品のポアソン比分布は、予め取得した成形材料（熱可塑性樹脂）の弾性定数（ヤング率）の温度依存性データおよびＰＶＴデータから求めた体積弾性率から計算した値を、流動解析により得た温度分布に当てはめることにより、求めることができる。
【0037】
（せん断弾性率分布）
射出成形品のせん断弾性率分布は、弾性定数（ヤング率）の温度依存性データおよびポアソン比の温度依存性データから計算した値を、流動解析により得た温度分布に当てはめることにより、求めることができる。
射出成形品のせん断弾性率分布は、弾性定数（ヤング率）の温度依存性データおよびポアソン比の温度依存性データから計算した値を、流動解析により得た温度分布に当てはめることにより、求めることができる。
【0038】
（体積収縮率分布）
射出成形品の体積収縮率分布は、予め取得（実測）した成形材料（熱可塑性樹脂）の圧力、体積、及び温度の関係を表すデータ（「ＰＶＴデータ」）を、流動解析により得た温度分布と圧力分布データに当てはめることにより、求めることができる。
図４は、ある熱可塑性樹脂のＰＶＴデータを示している。横軸は温度（単位は℃）であり、縦軸は密度の逆数、すなわち、比容積（単位はｃｍ３／ｇ）である。図４に示されているのは圧力（Ｐ）が５０ＭＰａの場合の温度と比容積の関係である。５０ＭＰａ以外のいくつかの圧力についても、温度と比容積の関係が予め実測されているとする。実測値のない圧力における温度と比容積の関係は、実測値のある圧力における温度と比容積の関係から、例えば内挿で取得できる。
図４では、圧力（Ｐ）が５０ＭＰaの場合に、ゲート部の熱可塑性樹脂が固化し、流動停止する時間（「ゲートシール時間」と言う）において、熱可塑性樹脂の温度が２００℃で、比容積が約０．８２ｃｍ３／ｇであることを示す。同様に、熱可塑性樹脂による射出成形品の温度が金型温度となる時点において、熱可塑性樹脂の温度が４０℃で、比容積が約０．７０ｃｍ３／ｇであることを示す。
射出成形品の体積収縮率分布は、予め取得（実測）した成形材料（熱可塑性樹脂）の圧力、体積、及び温度の関係を表すデータ（「ＰＶＴデータ」）を、流動解析により得た温度分布と圧力分布データに当てはめることにより、求めることができる。
図４は、ある熱可塑性樹脂のＰＶＴデータを示している。横軸は温度（単位は℃）であり、縦軸は密度の逆数、すなわち、比容積（単位はｃｍ３／ｇ）である。図４に示されているのは圧力（Ｐ）が５０ＭＰａの場合の温度と比容積の関係である。５０ＭＰａ以外のいくつかの圧力についても、温度と比容積の関係が予め実測されているとする。実測値のない圧力における温度と比容積の関係は、実測値のある圧力における温度と比容積の関係から、例えば内挿で取得できる。
図４では、圧力（Ｐ）が５０ＭＰaの場合に、ゲート部の熱可塑性樹脂が固化し、流動停止する時間（「ゲートシール時間」と言う）において、熱可塑性樹脂の温度が２００℃で、比容積が約０．８２ｃｍ３／ｇであることを示す。同様に、熱可塑性樹脂による射出成形品の温度が金型温度となる時点において、熱可塑性樹脂の温度が４０℃で、比容積が約０．７０ｃｍ３／ｇであることを示す。
【0039】
射出成形品の体積収縮率分布は、つぎのように算出する。
まず、流動解析により得た温度分布データと圧力分布データから、収縮前のある時点での温度（Ｔ１）と圧力（Ｐ１）を求める。これらを予め取得した熱可塑性樹脂のＰＶＴデータに当てはめることにより、収縮前の時点での体積（Ｖ１）に換算する。
図４のゲートシール時間の例では、Ｔ１＝２００℃、Ｐ１＝５０ＭＰａであって、比容積は約０．８２ｃｍ３／ｇである。つまり、この時点での単位重量（１ｇ）あたりの体積は、Ｖ１＝約０．８２ｃｍ３である。
同様に、流動解析により得た温度分布データと圧力分布データから、収縮後のある時点での温度（Ｔ２）と圧力（Ｐ２）を求める。これをＰＶＴデータに当てはめることにより、収縮後の時点での体積（Ｖ２）に換算する。
なお、収縮後のある時点については、射出開始から設定されたある時間（例えば実工程における、射出成形品が金型温度となる時間との比較で決めてもよい）が経過した時点として指定し得る。
図４の金型温度まで冷却された例では、Ｔ２＝４０℃、Ｐ２＝５０ＭＰａであって、比容積は約０．７０ｃｍ３／ｇである。つまり、この時点での単位重量（１ｇ）あたりの体積は、Ｖ２＝約０．７０ｃｍ３である。
そして収縮前後の各要素の体積（Ｖ１、Ｖ２）から、各要素の体積収縮率（（収縮前の体積－収縮後の体積）÷（収縮前の体積）、つまり、（Ｖ１－Ｖ２）／Ｖ１）の分布を求める。つまり、圧力が５０ＭＰａで一定（Ｐ１＝Ｐ２＝５０ＭＰａ）であって、Ｔ１＝２００℃からＴ２＝４０℃まで冷却するときの体積収縮率は、（Ｖ１－Ｖ２）／Ｖ１＝（０．８２－０．７０）／０．８２＝０．１４６、つまり、１４．６％と算出できる。
収縮の前後で圧力が一定でなくても、流動解析によって、収縮前の各要素の温度（Ｔ１）と圧力（Ｐ１）と、収縮後の各要素の温度（Ｔ２）と圧力（Ｐ２）が算出されている。収縮前後の圧力（Ｐ１、Ｐ２）における温度（Ｔ１、Ｔ２）と比容積の関係（ＰＶＴデータ）が予め取得されていれば、各要素の体積収縮率（Ｖ１－Ｖ２）／Ｖ１は算出できる。
このように、ＰＶＴデータと、温度分布及び圧力分布から体積収縮率分布を算出できる。
まず、流動解析により得た温度分布データと圧力分布データから、収縮前のある時点での温度（Ｔ１）と圧力（Ｐ１）を求める。これらを予め取得した熱可塑性樹脂のＰＶＴデータに当てはめることにより、収縮前の時点での体積（Ｖ１）に換算する。
図４のゲートシール時間の例では、Ｔ１＝２００℃、Ｐ１＝５０ＭＰａであって、比容積は約０．８２ｃｍ３／ｇである。つまり、この時点での単位重量（１ｇ）あたりの体積は、Ｖ１＝約０．８２ｃｍ３である。
同様に、流動解析により得た温度分布データと圧力分布データから、収縮後のある時点での温度（Ｔ２）と圧力（Ｐ２）を求める。これをＰＶＴデータに当てはめることにより、収縮後の時点での体積（Ｖ２）に換算する。
なお、収縮後のある時点については、射出開始から設定されたある時間（例えば実工程における、射出成形品が金型温度となる時間との比較で決めてもよい）が経過した時点として指定し得る。
図４の金型温度まで冷却された例では、Ｔ２＝４０℃、Ｐ２＝５０ＭＰａであって、比容積は約０．７０ｃｍ３／ｇである。つまり、この時点での単位重量（１ｇ）あたりの体積は、Ｖ２＝約０．７０ｃｍ３である。
そして収縮前後の各要素の体積（Ｖ１、Ｖ２）から、各要素の体積収縮率（（収縮前の体積－収縮後の体積）÷（収縮前の体積）、つまり、（Ｖ１－Ｖ２）／Ｖ１）の分布を求める。つまり、圧力が５０ＭＰａで一定（Ｐ１＝Ｐ２＝５０ＭＰａ）であって、Ｔ１＝２００℃からＴ２＝４０℃まで冷却するときの体積収縮率は、（Ｖ１－Ｖ２）／Ｖ１＝（０．８２－０．７０）／０．８２＝０．１４６、つまり、１４．６％と算出できる。
収縮の前後で圧力が一定でなくても、流動解析によって、収縮前の各要素の温度（Ｔ１）と圧力（Ｐ１）と、収縮後の各要素の温度（Ｔ２）と圧力（Ｐ２）が算出されている。収縮前後の圧力（Ｐ１、Ｐ２）における温度（Ｔ１、Ｔ２）と比容積の関係（ＰＶＴデータ）が予め取得されていれば、各要素の体積収縮率（Ｖ１－Ｖ２）／Ｖ１は算出できる。
このように、ＰＶＴデータと、温度分布及び圧力分布から体積収縮率分布を算出できる。
【0040】
（温度荷重）
各要素の体積収縮率を予め取得した熱可塑性樹脂の体積膨張率で割ることにより温度差に換算することで、温度荷重（温度差の分布）を求める。体積膨張率は元の体積Ｖが温度上昇ΔＴによってΔＶだけ変化した際の関係式ΔＶ／Ｖ＝βΔＴの係数βである。なお一般的には体積膨張率はＰＶＴデータより求めるが、熱可塑性樹脂が等方性の場合、体積膨張率は線膨張率の３倍の値であるから、ＰＶＴデータより求めた体積膨張率に代えて、予め測定した線膨張率から取得してもよい。体積膨張率と線膨張率を熱膨張率と総称する。熱膨張率（体積膨張率及び線膨張率）は温度依存性を有してもよい。
なお、温度荷重は「収縮前後の温度差」を意味するものであり、収縮率分布を求める過程で利用した「収縮前のある時点での温度」と「収縮後のある時点での温度」の差とは区別される。流動解析で得た温度データから差を求めた後者の温度差は、単に「温度のみ」を計算したものである。それに対して、前者の温度荷重は、一度ＰＶＴデータにより「温度と圧力」から体積を求めた上で体積膨張率を用いて温度に戻した値であることから、「温度のみ」ではなく「温度と圧力」を用いて計算したものである。よって、温度荷重は、実際の成形において樹脂が受ける圧力の影響を考慮した「温度差」となっている点で、収縮前のある時点と収縮後のある時点での温度のみの差とは異なり、温度荷重を使用する方がより精度の高い解析結果が得られる。
各要素の体積収縮率を予め取得した熱可塑性樹脂の体積膨張率で割ることにより温度差に換算することで、温度荷重（温度差の分布）を求める。体積膨張率は元の体積Ｖが温度上昇ΔＴによってΔＶだけ変化した際の関係式ΔＶ／Ｖ＝βΔＴの係数βである。なお一般的には体積膨張率はＰＶＴデータより求めるが、熱可塑性樹脂が等方性の場合、体積膨張率は線膨張率の３倍の値であるから、ＰＶＴデータより求めた体積膨張率に代えて、予め測定した線膨張率から取得してもよい。体積膨張率と線膨張率を熱膨張率と総称する。熱膨張率（体積膨張率及び線膨張率）は温度依存性を有してもよい。
なお、温度荷重は「収縮前後の温度差」を意味するものであり、収縮率分布を求める過程で利用した「収縮前のある時点での温度」と「収縮後のある時点での温度」の差とは区別される。流動解析で得た温度データから差を求めた後者の温度差は、単に「温度のみ」を計算したものである。それに対して、前者の温度荷重は、一度ＰＶＴデータにより「温度と圧力」から体積を求めた上で体積膨張率を用いて温度に戻した値であることから、「温度のみ」ではなく「温度と圧力」を用いて計算したものである。よって、温度荷重は、実際の成形において樹脂が受ける圧力の影響を考慮した「温度差」となっている点で、収縮前のある時点と収縮後のある時点での温度のみの差とは異なり、温度荷重を使用する方がより精度の高い解析結果が得られる。
【0041】
（解析用モデルの各要素に発生するひずみの算出（Ｓ５））
解析用モデルの各要素に発生するひずみを算出するステップ（Ｓ５）では、離型後の冷却過程における弾性定数分布（特にヤング率分布、ポアソン比分布及びせん断弾性率分布）と、体積収縮率分布及び／又は温度荷重分布（体積収縮率分布と温度荷重分布のうちの少なくとも１つ）を用いた構造解析を行い、解析用モデルの各要素に発生するひずみ（解析用モデルにおけるひずみ量の分布）を算出する。
なお、離型後の温度分布を求める際に、離型直後の初期温度を求める必要があり、金型内における温度分布にて初期温度を求めている。
ひずみの算出の時間的始点は、例えば、保圧冷却された離型の時点とすることができる。ひずみの算出の時間的終点は、例えば、射出成形品が室温となった時点とすることができる。あるいは、ひずみの算出の時間的終点は、射出開始から設定されたある時間（例えば実工程との比較により設定してもよい）を経過した時点とすることができる。
解析用モデルの各要素に発生するひずみを算出するステップ（Ｓ５）では、離型後の冷却過程における弾性定数分布（特にヤング率分布、ポアソン比分布及びせん断弾性率分布）と、体積収縮率分布及び／又は温度荷重分布（体積収縮率分布と温度荷重分布のうちの少なくとも１つ）を用いた構造解析を行い、解析用モデルの各要素に発生するひずみ（解析用モデルにおけるひずみ量の分布）を算出する。
なお、離型後の温度分布を求める際に、離型直後の初期温度を求める必要があり、金型内における温度分布にて初期温度を求めている。
ひずみの算出の時間的始点は、例えば、保圧冷却された離型の時点とすることができる。ひずみの算出の時間的終点は、例えば、射出成形品が室温となった時点とすることができる。あるいは、ひずみの算出の時間的終点は、射出開始から設定されたある時間（例えば実工程との比較により設定してもよい）を経過した時点とすることができる。
【0042】
（ボイドの発生場所及び／又はボイド量の予測（Ｓ６））
ボイドの発生場所及び／又はボイド量を予測するステップ（Ｓ６）においては、解析用モデルにおけるひずみ量の分布に基づいて変形量を算出してボイドの発生を予測する。
ひずみ量の分布からまずヒケの発生を予測して、ヒケが発生した形状についての解析をすることによって、ボイドの発生を予測することができる。
またはヒケ発生の予測をせずに、例えば、解析用モデルの各要素におけるひずみが、樹脂材料ごとに予め実測で求めたボイド発生の閾値以上となる場合に、その要素でボイドが発生すると予測することができる。
従来の解析ソフトでは、弾性定数（ヤング率、ポアソン比及びせん断弾性率）及び線膨張率の温度依存性を考慮しないことから、常温の弾性定数（一つの値）しか入力できないため、温度が変化した場合の各部の弾性定数の変化の考慮が不足しており、ボイド発生について正確な予測ができなかった。これに対して、本開示に係る改良手法では、流動解析で得た温度分布と圧力分布とに基づいて、各要素の温度変化による弾性定数分布（特にヤング率分布、ポアソン比分布及びせん断弾性率分布）に加え、体積収縮率分布より計算した温度荷重を用いて成形品各部のひずみ量が算出できるため、予測精度が向上している。
ボイドの発生場所及び／又はボイド量を予測するステップ（Ｓ６）においては、解析用モデルにおけるひずみ量の分布に基づいて変形量を算出してボイドの発生を予測する。
ひずみ量の分布からまずヒケの発生を予測して、ヒケが発生した形状についての解析をすることによって、ボイドの発生を予測することができる。
またはヒケ発生の予測をせずに、例えば、解析用モデルの各要素におけるひずみが、樹脂材料ごとに予め実測で求めたボイド発生の閾値以上となる場合に、その要素でボイドが発生すると予測することができる。
従来の解析ソフトでは、弾性定数（ヤング率、ポアソン比及びせん断弾性率）及び線膨張率の温度依存性を考慮しないことから、常温の弾性定数（一つの値）しか入力できないため、温度が変化した場合の各部の弾性定数の変化の考慮が不足しており、ボイド発生について正確な予測ができなかった。これに対して、本開示に係る改良手法では、流動解析で得た温度分布と圧力分布とに基づいて、各要素の温度変化による弾性定数分布（特にヤング率分布、ポアソン比分布及びせん断弾性率分布）に加え、体積収縮率分布より計算した温度荷重を用いて成形品各部のひずみ量が算出できるため、予測精度が向上している。
【0043】
さらに、本開示の方法では、金型外での冷却時における解析用モデルの各要素の温度分布及び圧力分布の算出すること（Ｓ３）により、金型外冷却を考慮したボイドの発生予測を行うことができるため、予測精度がさらに向上している。
【0044】
次に、本実施形態に係る方法（以下、単に「改良手法」とも言う）による実施例及び従来手法による比較例を示し、本開示の実施形態に係るヒケ及びボイドの予測について具体的に説明する。なお、本開示はこれらの実施例に限定されない。
【0045】
（製品形状及び成形材料）
図５に実施例に係る、ボルト状の厚肉でリブを有する成形品形状の３次元ＣＡＤモデルを示す。実成形品の成形材料（熱可塑性樹脂）はポリプラスチックス製ジュラコン（登録商標）ＰＯＭ Ｍ９０-４４（無充填材)である。なお、成形材料の比熱は示差走査熱量計（ＤＳＣ）により測定し、熱伝導率を交流定常法（ＩＳＯ２２００７－６）と呼ばれる熱伝導率測定手法にて求めた。
図５に示すように、六角形フランジ部の幅は２０ｍｍで厚さは６ｍｍである。また、円柱の高さは１５ｍｍである。六角形フランジ部の一辺にゲートがある。
図５に実施例に係る、ボルト状の厚肉でリブを有する成形品形状の３次元ＣＡＤモデルを示す。実成形品の成形材料（熱可塑性樹脂）はポリプラスチックス製ジュラコン（登録商標）ＰＯＭ Ｍ９０-４４（無充填材)である。なお、成形材料の比熱は示差走査熱量計（ＤＳＣ）により測定し、熱伝導率を交流定常法（ＩＳＯ２２００７－６）と呼ばれる熱伝導率測定手法にて求めた。
図５に示すように、六角形フランジ部の幅は２０ｍｍで厚さは６ｍｍである。また、円柱の高さは１５ｍｍである。六角形フランジ部の一辺にゲートがある。
【0046】
なお、実施例及び比較例における流動解析には、Autodesk社製 Moldflow（登録商標） Insight 2019.0.5 （３次元ソリッドモデル）ビルド20180921.0959_C70L71（以下、Moldflowと称する）を用いた。
実施例における金型外冷却解析及び構造解析には株式会社アライドエンジニアリング製 Adventure Cluster2021（以下、ADVCと称する）を用いた。
実施例における弾性定数分布及び体積収縮率分布の算出（Ｓ５）には、Visual Basic（登録商標）による自作の変換プログラムを使用した。
実施例における金型外冷却解析及び構造解析には株式会社アライドエンジニアリング製 Adventure Cluster2021（以下、ADVCと称する）を用いた。
実施例における弾性定数分布及び体積収縮率分布の算出（Ｓ５）には、Visual Basic（登録商標）による自作の変換プログラムを使用した。
【0047】
（実施例）
実成形品の成形条件と本手法による解析条件を表１に示した。実施例において、熱可塑性樹脂の物性値は、実成形品の熱可塑性樹脂の物性値と同じにしている。
実成形品の成形条件と本手法による解析条件を表１に示した。実施例において、熱可塑性樹脂の物性値は、実成形品の熱可塑性樹脂の物性値と同じにしている。
【表1】
【0048】
また、図６に実施例に使用した解析用モデルの断面図を示す。
【0049】
実施例においては射出、保圧及び金型内冷却の解析（図２のＳ２）を行い、金型外冷却解析（図２のＳ３）による温度分布から、変換プログラムによって、射出成形品（内部及び表面）の弾性定数分布（特にヤング率分布、ポアソン比分布及びせん断弾性率分布）と体積収縮率分布から計算した温度荷重を得て（図２のＳ４）、構造解析を行い、ひずみを算出して、ボイドの発生場所及びボイド量の予測を行った（図２のＳ５とＳ６）。
【0050】
（比較例）
従来手法に係る比較例では、実施例と同じ解析用モデルに対して、実施例と同じ熱可塑性樹脂の物性値（弾性定数等の温度依存性を除く）及び同じ解析条件を使用して、流動解析によって算出した体積収縮率分布の値からボイドの予測をした。前述のように、従来手法では、弾性定数（ヤング率、ポアソン比及びせん断弾性率）、及び線膨張率は温度依存性を考慮していないため、温度が変化した場合の各部の弾性定数の変化を考慮していない。つまり、比較例では金型内冷却の解析精度が劣る結果となっている。
そして、比較例ではボイドの予測のための連成解析は行わずに、流動解析から得た体積収縮率分布からボイドの予測をした。
従来手法に係る比較例では、実施例と同じ解析用モデルに対して、実施例と同じ熱可塑性樹脂の物性値（弾性定数等の温度依存性を除く）及び同じ解析条件を使用して、流動解析によって算出した体積収縮率分布の値からボイドの予測をした。前述のように、従来手法では、弾性定数（ヤング率、ポアソン比及びせん断弾性率）、及び線膨張率は温度依存性を考慮していないため、温度が変化した場合の各部の弾性定数の変化を考慮していない。つまり、比較例では金型内冷却の解析精度が劣る結果となっている。
そして、比較例ではボイドの予測のための連成解析は行わずに、流動解析から得た体積収縮率分布からボイドの予測をした。
【0051】
（参考例）
さらに参考例では、実施例と同じ解析用モデルに対して、実施例と同じ熱可塑性樹脂の物性値（弾性定数等の温度依存性も含む）及び同じ解析条件を使用して、連成解析を行い、構造解析で得られたひずみからボイドの予測をした。参考例における解析条件は表１の実成形品の成形条件と同じにしている。また、参考例において、熱可塑性樹脂の物性値は、実成形品の熱可塑性樹脂の物性値と同じにしている。
なお、参考例では金型外冷却の解析を行っていない。金型内冷却の解析だけ行っている。
さらに参考例では、実施例と同じ解析用モデルに対して、実施例と同じ熱可塑性樹脂の物性値（弾性定数等の温度依存性も含む）及び同じ解析条件を使用して、連成解析を行い、構造解析で得られたひずみからボイドの予測をした。参考例における解析条件は表１の実成形品の成形条件と同じにしている。また、参考例において、熱可塑性樹脂の物性値は、実成形品の熱可塑性樹脂の物性値と同じにしている。
なお、参考例では金型外冷却の解析を行っていない。金型内冷却の解析だけ行っている。
【0052】
（ボイド予測の比較）
図７に、実成形品、従来手法による比較例、金型外解析を行わない参考例及び改良手法による実施例における、ボイド発生の比較を示す。
図７のＡからＤ（以下図７Ａから図７Ｄと称する）には、実成形品のＸ線ＣＴ、流動解析のみの従来手法による比較例、流動解析と構造解析の連成解析を行うが、金型外冷却解析は行わない参考例、及び、改良手法による実施例をそれぞれ示す。
図７に、実成形品、従来手法による比較例、金型外解析を行わない参考例及び改良手法による実施例における、ボイド発生の比較を示す。
図７のＡからＤ（以下図７Ａから図７Ｄと称する）には、実成形品のＸ線ＣＴ、流動解析のみの従来手法による比較例、流動解析と構造解析の連成解析を行うが、金型外冷却解析は行わない参考例、及び、改良手法による実施例をそれぞれ示す。
【0053】
比較例では成形品内部の全体にボイドが発生しており、実成形品でのボイド発生箇所との乖離が大きくなる。
実施例では、実成形品と解析のボイド発生箇所が近づき、実際のボイド発生現象を連成解析で精度よく表現できるようになっている。
実施例では、実成形品と解析のボイド発生箇所が近づき、実際のボイド発生現象を連成解析で精度よく表現できるようになっている。
【0054】
（ボイドを低減する方法）
本開示には、熱可塑性樹脂を金型に射出成形してなる射出成形品の内部に生じるボイドを低減する方法も含まれる。
すなわち、上述したボイドの発生挙動を予測する方法によって、予測されたボイド量に対し、設計、成形条件及び成形材料のうちの１つ以上を変化させたときの予測結果を比較して、ボイド量が所定の量以下に低減されるまでこれを繰り返すようにする。
本開示には、熱可塑性樹脂を金型に射出成形してなる射出成形品の内部に生じるボイドを低減する方法も含まれる。
すなわち、上述したボイドの発生挙動を予測する方法によって、予測されたボイド量に対し、設計、成形条件及び成形材料のうちの１つ以上を変化させたときの予測結果を比較して、ボイド量が所定の量以下に低減されるまでこれを繰り返すようにする。
【0055】
本開示は、上述の実施形態に限定されるものではなく、上述の構成に対して、構成要素の付加、削除又は転換を行った様々な変形例も含むものとする。また、各実施例が様々に組み合わせることが可能である。
特に、本開示は、形状、材料又は条件について、上記の実施形態又は実施例に限定されて解釈されるべきではない。
特に、本開示は、形状、材料又は条件について、上記の実施形態又は実施例に限定されて解釈されるべきではない。
【0056】
さらに本開示に係るボイドの発生挙動を予測する方法、又は、本開示に係るボイドを低減する方法を１又は複数のプロセッサに実行させるためのプログラムも本開示に含まれる。当該プログラムは、コンピュータ読み取り可能で非一時的な（ｎｏｎ－ｔｒａｎｓｉｔｏｒｙ）記憶媒体に記録されて提供されてよい。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】