特開 2022-118541 繊維強化樹脂成形品の繊維配向状態予測方法、装置、プログラム及び記録媒体

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

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

(11)【公開番号】P2022118541

(43)【公開日】2022-08-15

(54)【発明の名称】繊維強化樹脂成形品の繊維配向状態予測方法、装置、プログラム及び記録媒体

(51)【国際特許分類】

   B29C 70/10 20060101AFI20220805BHJP

   B29C 70/42 20060101ALI20220805BHJP

   B29C 45/76 20060101ALI20220805BHJP

   G06F 30/10 20200101ALI20220805BHJP

   G06F 30/23 20200101ALI20220805BHJP

   G01N 25/02 20060101ALN20220805BHJP

   B29K 105/10 20060101ALN20220805BHJP

【ＦＩ】

B29C70/10

B29C70/42

B29C45/76

G06F17/50 680Z

G06F17/50 612H

G01N25/02 Z

B29K105:10

【審査請求】未請求

【請求項の数】14

【出願形態】ＯＬ

(21)【出願番号】P 2021015125

(22)【出願日】2021-02-02

(71)【出願人】

【識別番号】520167645

【氏名又は名称】東レエンジニアリングＤソリューションズ株式会社

(71)【出願人】

【識別番号】390026538

【氏名又は名称】ダイキョーニシカワ株式会社

(71)【出願人】

【識別番号】000003137

【氏名又は名称】マツダ株式会社

(74)【代理人】

【識別番号】110001427

【氏名又は名称】弁理士法人前田特許事務所

(72)【発明者】

【氏名】岡田 有司

(72)【発明者】

【氏名】末吉 耕平

(72)【発明者】

【氏名】石田 元伸

(72)【発明者】

【氏名】黒川 大輔

(72)【発明者】

【氏名】古川 智司

(72)【発明者】

【氏名】田中 宣隆

(72)【発明者】

【氏名】鈴木 広之

(72)【発明者】

【氏名】田中 慶和

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

2G040

4F205

4F206

5B046

5B146

【Ｆターム（参考）】

2G040AA01

2G040BA02

2G040BA22

4F205AA11

4F205AB25

4F205AD16

4F205AR20

4F205HA12

4F205HA27

4F205HA34

4F205HA36

4F205HB01

4F205HK19

4F205HL15

4F206AA11

4F206AB11

4F206AB25

4F206AD16

4F206AM23

4F206AR20

4F206JA07

4F206JL09

5B046AA05

5B046JA09

5B146DJ03

5B146DJ07

(57)【要約】

【課題】簡単な手法で繊維配向の予測精度を向上できる繊維強化樹脂成形品の繊維配向状態予測方法、装置、プログラム及び記録媒体を提供する。
【解決手段】繊維強化樹脂成形品の繊維配向状態予測方法は、繊維強化樹脂の射出成形品に含まれる繊維の配向状態をコンピュータシミュレーションにより予測する方法であって、解析モデル作成工程と、樹脂の流速を算出する流動解析工程と、樹脂の流速に基づき、繊維配向テンソルを算出する繊維配向テンソル算出工程と、を備える。繊維配向テンソル算出工程で、樹脂におけるキャビティに接触する表面から所定深さまでの領域を最表層部、最表層部よりも深い領域を内側部としたときに、内側部の繊維配向テンソルに基づいて、最表層部の繊維配向テンソルを算出する。
【選択図】図２

【特許請求の範囲】

【請求項1】

繊維強化樹脂の射出成形品に含まれる繊維の配向状態をコンピュータシミュレーションにより予測する方法であって、
鋳型のキャビティの形状データを複数の微小要素に分割して解析モデルを作成する解析モデル作成工程と、
前記キャビティに前記繊維強化樹脂を射出したときの、該繊維強化樹脂に含まれる樹脂の流速を算出する流動解析工程と、
前記樹脂の前記流速に基づき、前記繊維の配向状態を示す繊維配向テンソルを算出する繊維配向テンソル算出工程と、を備え、
前記繊維配向テンソル算出工程で、前記樹脂における前記キャビティに接触する表面から所定深さまでの領域を最表層部、該最表層部よりも深い領域を内側部としたときに、前記内側部の繊維配向テンソルに基づいて、前記最表層部の繊維配向テンソルを算出する
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測方法。

【請求項2】

請求項１において、
前記繊維配向テンソル算出工程で、前記最表層部の繊維配向テンソルは、前記内側部のフローフロント部位における繊維配向テンソルの平均値として算出される
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測方法。

【請求項3】

請求項２において、
前記内側部は、前記最表層部に接する表層部と、該表層部の内側の中心部と、を備え、
前記平均値は、少なくとも前記表層部の繊維配向テンソルと前記中心部の繊維配向テンソルとの平均値である
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測方法。

【請求項4】

請求項２又は請求項３において、
前記平均値は、前記内側部を深さ方向に７以上１８以下の層に分割し、該層の各々について得られた繊維配向テンソルの平均値である
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測方法。

【請求項5】

請求項１～４のいずれか一において、
前記最表層部における前記所定深さは、０．２ｍｍ以下である
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測方法。

【請求項6】

請求項１～５のいずれか一において、
前記解析モデルの前記微小要素は、２．５次元薄肉シェル要素である
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測方法。

【請求項7】

繊維強化樹脂の射出成形品に含まれる繊維の配向状態をコンピュータシミュレーションにより予測する装置であって、
鋳型のキャビティの形状データを複数の微小要素に分割して解析モデルを作成する解析モデル作成手段と、
前記キャビティに前記繊維強化樹脂を射出したときの、該繊維強化樹脂に含まれる樹脂の流速を算出する流動解析手段と、
前記樹脂の前記流速に基づき、前記繊維の配向状態を示す繊維配向テンソルを算出する繊維配向テンソル算出手段と、を備え、
前記繊維配向テンソル算出手段は、前記樹脂における前記キャビティに接触する表面から所定深さまでの領域を最表層部、該最表層部よりも深い領域を内側部としたときに、前記内側部の繊維配向テンソルに基づいて、前記最表層部の繊維配向テンソルを算出する
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測装置。

【請求項8】

請求項７において、
前記最表層部の繊維配向テンソルは、前記内側部のフローフロント部位における繊維配向テンソルの平均値として算出される
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測装置。

【請求項9】

請求項８において、
前記内側部は、前記最表層部に接する表層部と、該表層部の内側の中心部と、を備え、
前記平均値は、少なくとも前記表層部の繊維配向テンソルと前記中心部の繊維配向テンソルとの平均値である
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測装置。

【請求項10】

請求項８又は請求項９において、
前記平均値は、前記内側部を深さ方向に分割してなる７以上１８以下の層の各々について得られた繊維配向テンソルの平均値である
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測装置。

【請求項11】

請求項７～１０のいずれか一において、
前記最表層部における前記所定深さは、０．２ｍｍ以下である
ことを特徴とする繊維強化樹脂の繊維配向状態予測装置。

【請求項12】

請求項７～１１のいずれか一において、
前記解析モデルの前記微小要素は、２．５次元薄肉シェル要素である
ことを特徴とする繊維強化樹脂成形品の繊維配向状態予測装置。

【請求項13】

請求項１～６のいずれか一に記載された繊維強化樹脂成形品の繊維配向状態予測方法をコンピュータに実行させるためのプログラム。

【請求項14】

請求項１３に記載されたプログラムを記録したコンピュータ読み取り可能な記録媒体。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、繊維強化樹脂成形品の繊維配向状態予測方法及びその装置、該方法をコンピュータに実行させるためのプログラム、並びに該プログラムを記録したコンピュータ読み取り可能な記録媒体に関する。

【背景技術】

【0002】

従来、繊維強化樹脂の射出成形品における製品設計等の容易化を目的として、鋳型内の繊維強化樹脂の流動固化挙動をＣＡＥ（Computer－Aided－Engineering）を用いて解析することが行われている（例えば、特許文献１，２参照）。

【0003】

繊維強化樹脂では、繊維の配向状態に異方性が生じると、成形品の収縮量、弾性率等にも異方性が生じ、延いては成形品の寸法精度、強度等に影響する。従って、流動樹脂中の繊維配向を精度よく予測することは、射出成形品の寸法精度、強度等を予測する上で重要である（例えば、特許文献３参照）。

【先行技術文献】

【特許文献】

【0004】

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

【特許文献2】特開２０１５－１８９１１７号公報

【特許文献3】特開２０１４－２２６８７１号公報

【発明の概要】

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

【0005】

ところで、従来のＣＡＥ解析では、流動樹脂中の繊維配向挙動、すなわち繊維の向き及びその分布度合いは、流動解析により得られる樹脂の流速に基づいて算出している。

【0006】

一方、射出成形品の繊維配向の実測結果では、最表層部は流動部に比べ異方性が低下する傾向にある。しかし、従来式により計算した場合、最表層部ではせん断流の影響が大きくなり、樹脂の流動方向に繊維が強く配列するという結果が得られ、繊維配向の予測精度が低下する。

【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】

【図1】本開示に係る繊維配向状態予測装置の構成例を示す図。

【図2】本開示に係る繊維配向状態予測方法を説明するためのフローチャート。

【図3】樹脂の２種類の流れと繊維の配向との関係を説明するための図。

【図4】フローフロント部位を含む先端領域における樹脂の流れと繊維の配向との関係を説明するための図。

【図5】試験片を厚さ方向に９層に分割した様子を示す図。

【図6】図５の層Ｌ２（表層部）及び層Ｌ５（中心部）のＸ線－ＣＴ画像。

【図7】繊維配向テンソルの定義を説明するための図。

【図8】図５の層Ｌ２（表層部）における角度毎の相対繊維量のヒストグラム。

【図9】図５の層Ｌ５（中心部）における角度毎の相対繊維量のヒストグラム。

【図10】図５の各層における繊維配向テンソルの実測値を示すグラフ。

【図11】実現象のメカニズムと、従来のＣＡＥによる繊維配向予測方法のメカニズムとの比較を示す図。

【図12】最表層部の繊維配向テンソルの実測値と流動部の繊維配向テンソルの平均値との関係を示すグラフ。

【図13】試験片の解析モデルを示す図。

【図14】厚さ１．５ｍｍの試験片における繊維配向テンソルの実測値並びに実施例及び比較例の予測値を示すグラフ。

【図15】厚さ２．０ｍｍの試験片における繊維配向テンソルの実測値並びに実施例及び比較例の予測値を示すグラフ。

【図16】厚さ２．５ｍｍの試験片における繊維配向テンソルの実測値並びに実施例及び比較例の予測値を示すグラフ。

【図17】厚さ３．０ｍｍの試験片における繊維配向テンソルの実測値並びに実施例及び比較例の予測値を示すグラフ。

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

【0038】

以下、本開示の実施形態を図面に基づいて詳細に説明する。以下の好ましい実施形態の説明は、本質的に例示に過ぎず、本開示、その適用物或いはその用途を制限することを意図するものでは全くない。

【0039】

＜繊維配向状態予測装置＞
図１に、本開示の繊維強化樹脂成形品の繊維配向状態予測装置１００（以下、「予測装置１００」ともいう。）の構成例を示す。予測装置１００は、繊維強化樹脂の射出成形品に含まれる繊維の配向状態をコンピュータシミュレーションにより予測する装置であり、コンピュータ１１０を基本構成とするＣＡＥ（Ｃｏｍｐｕｔｅｒ Ａｉｄｅｄ Ｅｎｇｉｎｅｅｒｉｎｇ）システムである。

【0040】

予測装置１００は、記憶部１２０と、プロセッサ１３０と、を備える。また、予測装置１００は、例えばディスプレイ等からなる表示部１４０、キーボード等からなる入力部１５０、及び各種記録媒体１７０に保存された情報を取得するための読取部１６０等を備える。記憶部１２０及び／又は記録媒体１７０には、演算処理用のプログラム及び各種解析用データ等の情報が格納される。プロセッサ１３０は、記憶部１２０に格納された上記情報、入力部１５０を介して入力された情報、及び読取部１６０を介して記録媒体１７０から取得した情報等に基づいて、各種演算処理を行う。

【0041】

予測装置１００は、解析モデル作成手段１３１により、繊維強化樹脂の流路となる鋳型のキャビティを定義した３Ｄ ＣＡＤデータ等の形状データを複数の微小要素に分割して解析モデルを作成する。微小要素を作成する解析モデル作成手段１３１としては、例えば東レエンジニアリングＤソリューションズ株式会社製の３Ｄ ＴＩＭＯＮ（登録商標）－Ｐｒｅ／Ｐｏｓｔ、株式会社エヌ・エス・ティ製のＦＥＭＡＰ（登録商標）、エムエスシーソフトウェア株式会社製のＰａｔｒａｎ（登録商標）、Ａｌｔａｉｒ社製のＨｙｅｐｅｒ ｍｅｓｈ（登録商標）等のＣＡＥプリプロセッサを使用できる。

【0042】

樹脂及び繊維等の種類、配合等に関する材料特性データ、射出速度及び樹脂温度等の成形条件が記載された境界条件データ等に基づき、解析条件設定手段１３２で解析条件を設定する。そして、流動解析手段１３３で、キャビティに材料を射出したときの樹脂の挙動を解析する流動解析を実行する。流動解析により、微小要素ごとの樹脂の流速の情報を含む各種データが得られる。なお、本明細書において、「流速」の語は、「流れの速度」と同義であり、流れの速さ及び方向の両者を含む概念として使用している。解析条件設定手段１３２及び流動解析手段１３３としては、例えば東レエンジニアリングＤソリューションズ株式会社製の３Ｄ ＴＩＭＯＮ（登録商標）等の射出成形ＣＡＥソフトウェアを使用できる。

【0043】

繊維配向テンソル算出手段１３４は、流動解析手段１３３により算出された樹脂の流速に基づき、繊維の配向状態を示す繊維配向テンソルを算出する。繊維配向テンソル算出手段１３４としては、例えば東レエンジニアリングＤソリューションズ株式会社製の３Ｄ ＴＩＭＯＮ（登録商標）－ＦＩＢＥＲ等の繊維配向解析ソルバーを活用できる。

【0044】

＜射出成形品＞
解析対象の射出成形品としては、特に限定する意図ではないが、例えば自動車用部品、ロケット、航空機等の部品、スポーツ用品等が挙げられる。好ましくは、車両の内外装部材等の板状の射出成形品が挙げられる。

【0045】

樹脂としては、特に限定されるものではなく、周知の樹脂を対象とすることができるが、具体的には例えばポリプロピレン樹脂、ポリカーボネート樹脂、ポリエステル樹脂、ポリアミド（ＰＡ）樹脂等が挙げられる。これらの樹脂は１種又は２種以上が混合されて用いられ得る。繊維としては、特に限定されるものではなく、周知の繊維を用いることができるが、具体的には例えばガラス繊維、炭素繊維、セルロースナノ繊維等が挙げられる。これらの繊維は１種又は２種以上が混合されて用いられ得る。繊維径、繊維長等は、特に限定されるものではなく、一般的に用いられる条件とすることができる。成形品中における繊維の含有量は、特に限定されるものではなく、一般的な含有量とすることができるが、例えば１質量％以上５０質量％以下、好ましくは５質量％以上３０質量％以下である。また、成形品は、成形性、強度、意匠性、機能性等の向上の観点から、５質量％程度のフィラー、顔料、染料、耐衝撃性改良剤、ＵＶ吸収剤等の添加材等を含有してもよい。これらの添加材は単独で又は複数種添加され得る。

【0046】

なお、上述の解析モデルの微小要素としては、特に限定されるものではなく、ソリッド要素、シェル要素等の周知の要素を採用できるが、板状の射出成形品を解析対象とする場合には、２．５次元薄肉シェル要素を採用することが望ましい。形状データを薄い板状のシェル要素に分割することで、計算量を抑え、簡単且つ短時間で繊維配向の予測を精度よく行うことができる。

【0047】

＜繊維配向状態予測方法＞
図２は、本開示に係る繊維強化樹脂成形品の繊維配向状態予測方法（以下、「予測方法」ともいう。）の実施の手順を示すフローチャートである。予測方法は、解析モデル作成工程Ｓ１と、解析条件設定工程Ｓ２と、流動解析工程Ｓ３と、繊維配向テンソル算出工程Ｓ４と、を備える。

【0048】

まず、解析モデル作成工程Ｓ１において、上述のごとく、解析モデル作成手段１３１により、３次元ＣＡＤ等を用いて作成した鋳型のキャビティの形状データ等を数値解析用の微小要素に分割し、解析モデルを作成する。

【0049】

続いて、解析条件設定工程Ｓ２において、材料特性データ、境界条件データ等の解析条件を設定する。そして、流動解析工程Ｓ３において、流動解析を行い、樹脂の流速を算出する。

【0050】

次に、繊維配向テンソル算出工程Ｓ４において、流動解析の結果得られた樹脂の流速に基づき、繊維配向テンソルを算出する。そして、繊維強化樹脂の充填が完了するまで流動解析Ｓ３及び繊維配向テンソル算出工程Ｓ４の工程が繰り返される（工程Ｓ５）。

【0051】

繊維配向テンソル算出工程Ｓ４では、例えば下記式（１）で示されるＡｄｖａｎｉ－Ｔｕｃｋｅｒの式を用いて繊維の配向状態が予測される（特許文献３参照）。

【0052】

【数1】

【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】

１００ 繊維配向状態予測装置、予測装置
１２０ 記憶部
１３０ プロセッサ
１３１ 解析モデル作成手段
１３２ 解析条件設定手段
１３３ 流動解析手段
１３４ 繊維配向テンソル算出手段
１７０ 記録媒体
２００ 樹脂
２１０ 先端領域
２１１ フローフロント部位
２１２ 流動部（内側部）
２１３ 中心部
２１５ 表層部
２１７ 最表層部
Ｃ キャビティ
Ｍ 鋳型
Ｓ１ 解析モデル作成工程
Ｓ２ 解析条件設定工程
Ｓ３ 流動解析工程
Ｓ４ 繊維配向テンソル算出工程

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】