(57)【特許請求の範囲】
【請求項１】
メッシュモデルを用いて成形型のキャビティ内の充填材を含有した溶融樹脂を流動解析して成形品表面から微小量突出する２次ウエルドラインの発生を予測する樹脂成型品の流動解析方法であって、
成形工程初期に発生する１次ウエルドラインの発生位置を演算する第１のステップと、
前記演算された１次ウエルドライン発生位置に基づいて、前記樹脂成型品の板厚方向の平均繊維配向と、前記溶融樹脂の体積収縮率とを用いて２次ウエルドラインの突出高さを演算する第２のステップとを備え、
２次ウエルドラインの発生位置におけるキャビティの板厚方向の幅をＴ１ｍｏｌｄ、一般部位置におけるキャビティの板厚方向の幅をＴ２ｍｏｌｄ、２次ウエルドラインの発生位置における体積収縮率をＶ１、一般部位置における体積収縮率をＶ２、２次ウエルドラインの発生位置における充填材の板厚方向に対する平均繊維配向度をｆ１、一般部位置における充填材の板厚方向に対する平均繊維配向度をｆ２、２次ウエルドラインの突出高さをｈとしたとき、ｈを、
ｈ＝（１−Ｖ１）・Ｔ１ｍｏｌｄ・ｆ１−（１−Ｖ２）・Ｔ２ｍｏｌｄ・ｆ２
の関係式に基づき演算することを特徴とする樹脂成型品の流動解析方法。
【請求項２】
前記第１のステップと第２のステップとの間において、前記溶融樹脂の流動方向の平均繊維配向度と前記溶融樹脂の平均移動速度とに基づいて２次ウエルドラインの発生位置を演算する第３のステップを備えたことを特徴とする請求項１に記載の樹脂成型品の流動解析方法。
【請求項３】
前記第３のステップと第２のステップとの間において、前記２次ウエルドライン発生位置の少なくとも一部の周辺の解析メッシュ単位を小さい解析メッシュ単位に変更する第４のステップを備えたことを特徴とする請求項２に記載の樹脂成型品の流動解析方法。
【請求項４】
メッシュモデルにより成形型内の充填材を含有した溶融樹脂を流動解析して成形品表面から突出する２次ウエルドラインの発生を予測する樹脂成型品の流動解析装置であって、
成形工程初期に発生する１次ウエルドラインの発生位置を演算する１次ウエルドライン作成部と、
前記演算された１次ウエルドライン発生位置に基づき前記溶融樹脂の体積収縮率と、前記樹脂成型品の板厚方向の平均繊維配向とを用いて２次ウエルドラインの高さを演算する２次ウエルドライン作成部とを備え、
前記２次ウエルドライン作成部は、２次ウエルドラインの発生位置におけるキャビティの板厚方向の幅をＴ１ｍｏｌｄ、一般部位置におけるキャビティの板厚方向の幅をＴ２ｍｏｌｄ、２次ウエルドラインの発生位置における体積収縮率をＶ１、一般部位置における体積収縮率をＶ２、２次ウエルドラインの発生位置における充填材の板厚方向に対する平均繊維配向度をｆ１、一般部位置における充填材の板厚方向に対する平均繊維配向度をｆ２、２次ウエルドラインの突出高さをｈとしたとき、ｈを、
ｈ＝（１−Ｖ１）・Ｔ１ｍｏｌｄ・ｆ１−（１−Ｖ２）・Ｔ２ｍｏｌｄ・ｆ２
の関係式に基づき演算することを特徴とする樹脂成型品の流動解析装置。
【請求項５】
メッシュモデルにより成形型内の充填材を含有した溶融樹脂を流動解析して成形品表面から突出する２次ウエルドラインの発生を予測する樹脂成型品の流動解析装置のコンピュータを、
成形工程初期に発生する１次ウエルドラインの発生位置を演算する１次ウエルドライン作成手段と、
前記演算された１次ウエルドライン発生位置に基づき前記溶融樹脂の体積収縮率と、前記樹脂成型品の板厚方向の平均繊維配向とを用いて２次ウエルドラインの高さを演算する２次ウエルドライン作成手段であって、
２次ウエルドラインの発生位置におけるキャビティの板厚方向の幅をＴ１ｍｏｌｄ、一般部位置におけるキャビティの板厚方向の幅をＴ２ｍｏｌｄ、２次ウエルドラインの発生位置における体積収縮率をＶ１、一般部位置における体積収縮率をＶ２、２次ウエルドラインの発生位置における充填材の板厚方向に対する平均繊維配向度をｆ１、一般部位置における充填材の板厚方向に対する平均繊維配向度をｆ２、２次ウエルドラインの突出高さをｈとしたとき、ｈを、
ｈ＝（１−Ｖ１）・Ｔ１ｍｏｌｄ・ｆ１−（１−Ｖ２）・Ｔ２ｍｏｌｄ・ｆ２
の関係式に基づき演算する２次ウエルドライン作成手段と、
して機能させるための流動解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、樹脂成型品の流動解析技術に関し、特に樹脂成型品の板厚方向の平均繊維配向度と樹脂の体積収縮率とを用いた樹脂成型品の流動解析方法、流動解析装置及び流動解析プログラムに関する。
【背景技術】
【０００２】
従来、樹脂の射出成形は、成形型のキャビティ内へ溶融樹脂を射出して１工程で複雑な形状の樹脂成型品を成形できるため、幅広い分野において樹脂成型品の製造に用いられている。自動車分野では、バンパ等の強度が要求される部品にも補強用の充填材（ガラスファイバ、タルク等）を含有した樹脂成型品が多く採用されている。この射出成形により成形された樹脂成型品には、ウエルドラインと呼ばれるノッチ状の線状痕や微小量突出した線状痕が成形品表面に発生することがあり、特に充填材として扁平状のタルク等を混合した場合にはウエルドラインが発生し易く外観上の不具合となっていた。そこで、射出成形プロセスの流動シミュレーションを行い、ウエルドラインの発生位置とその量とをコンピュータ上で予測し、その予測値に基づき射出ゲートの位置、射出タイミング、射出圧力等を調節してウエルドライン発生位置を変更する等の対策が行われている。
【０００３】
特許文献１に記載された樹脂成型品の流動解析方法は、少なくとも１つの射出ゲートから金型内部へ樹脂を射出する解析モデルを用いて、各ゲートから射出され対向方向から流動してきた流動樹脂が相互に接触（会合）するとき、接触部に生じる１次ウエルドライン上に複数の仮想粒子を生成し、これら仮想粒子の移動経路を求めてその移動距離を算出し流動樹脂の一方が他方の流動樹脂内へ侵入した潜り込み距離を求め、潜り込み距離から樹脂成型品の表面に生じる表面凹凸部（２次ウエルドライン）の発生位置と表面凹凸量（高さ・深さ）を予測している。この流動解析方法では、潜り込み距離と表面凹凸量との相関式を樹脂の種類毎に予め試作試験等により求めておき、この相関式と仮想粒子の予測潜り込み距離とにより表面凹凸量を予測する。
【０００４】
特許文献２に記載された樹脂成型品の品質予測方法は、メッシュモデルを用いて溶融樹脂の流動解析を行い、流動解析結果に基づき、２次ウエルドラインの発生位置を抽出し、ウエルドライン周辺での樹脂の流れ方向に沿った移動量の経時的変化を推定し、２次ウエルドライン周辺での樹脂の流動層厚み方向の経時的変化を推定し、流れ方向経時変化と厚み方向経時変化とに基づき樹脂会合部であるウエルド界面形状を推定して成型品表面を評価している。この品質予測方法では、ウエルド界面が途中で折り返さない場合に比べてウエルド界面が途中で折り返す場合に小さくなる値、また、ウエルド界面の厚さ方向中央部での曲率が大きい程小さくなる値となる評価関数を予め試作試験等により求めておき、この評価関数とウエルド界面形状とにより成型品表面の品質を予測する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０００−３４３５７５号公報
【特許文献２】特開２００８−２０７４４０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
成形時に成形型のキャビティ内において２つ以上の溶融樹脂の流動層が会合する樹脂会合部の表面に、ノッチ状の１次ウエルドラインと、この１次ウエルドラインから所定距離離隔した位置に微小量突出する２次ウエルドラインとが夫々発生する。そこで、本発明者は、１次ウエルドラインと２次ウエルドラインとの発生メカニズムを検討した結果、以下の点に着眼した。
図１５に示すように、会合前段階では、溶融樹脂内にタルク等の扁平状の充填材を含有する場合、充填材は剪断応力の発生方向に応じて配向する特性を有している。それ故、剪断応力の大きな溶融樹脂と金型との境界部分の充填材は流動方向（横方向）に配向し、剪断応力の小さな中央部分の充填材はキャビティの板厚方向（縦方向）に配向する。
【０００７】
溶融樹脂の流動層が会合する会合段階では、夫々の流動層の進行方向先端部が衝突してキャビティの板厚方向に樹脂会合部を形成する。この樹脂会合部周辺の充填材は縦配向し、流動層と金型との境界部分の充填材は樹脂会合部（流動層の進行方向）側程キャビティの板厚方向内側へ傾斜状になるよう傾斜配向している。それ故、樹脂会合部の表面にはＶノッチ状に凹入した１次ウエルド部が形成され、この１次ウエルド部が線状に連続して１次ウエルドラインが形成されている。
【０００８】
会合段階後の内部流動段階では、未硬化の板厚方向内側の内部流動層が高圧力側から低圧力側へ移動する。樹脂会合部の移動に伴い、縦配向した充填材もその数を略維持しながら低圧力側へ移動する。夫々の内部流動層は溶融樹脂の硬化が進むにつれて流動量が減少し、流動層の溶融樹脂はキャビティの板厚方向外側から固化層を形成している。
【０００９】
夫々の内部流動層の圧力が均衡する圧力均衡段階では、内部流動層及びウエルド界面の移動が終了する。圧力均衡段階における樹脂会合部周辺には高い密度で縦配向した充填材が存在し、また、金型との境界部分の充填材は樹脂会合部側程キャビティの板厚方向外側へ傾斜状になるよう傾斜配向している。それ故、樹脂会合部周辺の樹脂の容量は少なく、樹脂会合部以外の部分（一般部分）の樹脂の容量は大きいため、換言すれば、圧力均衡段階における樹脂会合部は縦配向した充填材により一般部分に比べて樹脂が収縮し難いため、型開き後、溶融樹脂の硬化時における樹脂会合部周辺の収縮量は小さく、樹脂会合部以外の部分の収縮量は大きくなる。即ち、樹脂表面には一般部分に比べて微小量突出した２次ウエルド部が形成され、この２次ウエルド部が線状に連続して２次ウエルドラインが形成されている。以上により、１次ウエルドラインよりも２次ウエルドラインが製品品質に与える影響が大きく、ウエルドラインは樹脂会合部における充填材の板厚方向に対する配向傾向（配向度）によってその突出量が異なることに着目した。
【００１０】
特許文献１の流動解析方法では、潜り込み距離から樹脂成型品の表面に生じる表面凹凸部の発生位置と表面凹凸量を相関式に基づき予測している。しかし、潜り込み距離と表面凹凸量との相関式は、予め試作試験等により求められた相関式であり、樹脂の種類、樹脂溶融温度、充填材の含有率等の射出条件毎に数多くの相関式を準備する必要がある。しかも、樹脂に含有されるタルクの配向は、潜り込み距離から直接的に演算するものではなく、相関式の準備に多大な労力を費やすにも拘わらず、実際の樹脂成形品に即した高い精度の表面凹凸量を得ることができない虞がある。
【００１１】
特許文献２の品質予測方法では、樹脂会合部近傍とそれ以外の一般部分で板厚方向のタルク配向が相違することに着目し、成形時における樹脂会合部形状を推定して評価関数に基づき成型品表面を評価している。しかし、この評価関数は、特許文献１の流動解析方法と同様に、予め試作試験等により求められた評価関数であり、樹脂の種類、樹脂溶融温度、充填材の含有率等の射出条件毎に数多くの評価関数を準備する必要がある。しかも、ウエルド界面形状を推定するに留まり、樹脂に含有されるタルクの配向度に基づく表面凹凸量を演算するものではないため、高精度の表面凹凸量や表面凹凸量に基づく品質予測を得ることは難しい。
【００１２】
本発明の目的は、樹脂成型品の２次ウエルドラインの突出高さを予測でき、高い精度の品質評価を可能にできる樹脂成型品の流動解析方法、流動解析装置及び流動解析プログラム等を提供することである。
【課題を解決するための手段】
【００１３】
請求項１の樹脂成型品の流動解析方法は、メッシュモデルを用いて成形型のキャビティ内の充填材を含有した溶融樹脂を流動解析して成形品表面から微小量突出する２次ウエルドラインの発生を予測する樹脂成型品の流動解析方法であって、成形工程初期に発生する１次ウエルドラインの発生位置を演算する第１のステップと、前記演算された１次ウエルドライン発生位置に基づいて、前記樹脂成型品の板厚方向の平均繊維配向度と、前記溶融樹脂の体積収縮率とを用いて２次ウエルドラインの突出高さを演算する第２のステップとを備え、２次ウエルドラインの発生位置におけるキャビティの板厚方向の幅をＴ１ｍｏｌｄ、一般部位置におけるキャビティの板厚方向の幅をＴ２ｍｏｌｄ、２次ウエルドラインの発生位置における体積収縮率をＶ１、一般部位置における体積収縮率をＶ２、２次ウエルドラインの発生位置における充填材の板厚方向に対する平均繊維配向度をｆ１、一般部位置における充填材の板厚方向に対する平均繊維配向度をｆ２、２次ウエルドラインの突出高さをｈとしたとき、ｈを、  
ｈ＝（１−Ｖ１）・Ｔ１ｍｏｌｄ・ｆ１−（１−Ｖ２）・Ｔ２ｍｏｌｄ・ｆ２
の関係式に基づき演算することを特徴としている。
【００１４】
この流動解析方法においては、演算された１次ウエルドライン発生位置に基づき２次ウエルドラインを予測するため、実際に発生する２次ウエルドラインに対して高精度の２次ウエルドラインの発生予測ができる。しかも、ガラスファイバ等強化繊維の水平方向の並び状態を示す指標である平均繊維配向をタルク等の扁平状の充填材の配向に代用できるため、充填材の配向を演算により容易に得ることができる。
【００１５】
請求項２の発明は、請求項１の発明において、前記第１のステップと第２のステップとの間において、前記溶融樹脂の流動方向の平均繊維配向度と前記溶融樹脂の平均移動速度とに基づいて２次ウエルドラインの発生位置を演算する第３のステップを備えたことを特徴としている。
【００１６】
請求項３の発明は、請求項２の発明において、前記第３のステップと第２のステップとの間において、前記２次ウエルドライン発生位置の少なくとも一部の周辺の解析メッシュ単位を小さい解析メッシュ単位に変更する第４のステップを備えたことを特徴としている。
【００１７】
【００１８】
請求項４の樹脂成型品の流動解析装置は、メッシュモデルにより成形型内の充填材を含有した溶融樹脂を流動解析して成形品表面から突出する２次ウエルドラインの発生を予測する樹脂成型品の流動解析装置であって、成形工程初期に発生する１次ウエルドラインの発生位置を演算する１次ウエルドライン作成部と、前記演算された１次ウエルドライン発生位置に基づき前記溶融樹脂の体積収縮率と、前記樹脂成型品の板厚方向の平均繊維配向度とを用いて２次ウエルドラインの高さを演算する２次ウエルドライン作成部とを備え、前記２次ウエルドライン作成部は、２次ウエルドラインの発生位置におけるキャビティの板厚方向の幅をＴ１ｍｏｌｄ、一般部位置におけるキャビティの板厚方向の幅をＴ２ｍｏｌｄ、２次ウエルドラインの発生位置における体積収縮率をＶ１、一般部位置における体積収縮率をＶ２、２次ウエルドラインの発生位置における充填材の板厚方向に対する平均繊維配向度をｆ１、一般部位置における充填材の板厚方向に対する平均繊維配向度をｆ２、２次ウエルドラインの突出高さをｈとしたとき、ｈを、  
ｈ＝（１−Ｖ１）・Ｔ１ｍｏｌｄ・ｆ１−（１−Ｖ２）・Ｔ２ｍｏｌｄ・ｆ２
の関係式に基づき演算することを特徴としている。
【００１９】
請求項５の流動解析プログラムは、メッシュモデルにより成形型内の充填材を含有した溶融樹脂を流動解析して成形品表面から突出する２次ウエルドラインの発生を予測する樹脂成型品の流動解析装置のコンピュータを、成形工程初期に発生する１次ウエルドラインの発生位置を演算する１次ウエルドライン作成手段と、前記演算された１次ウエルドライン発生位置に基づき前記溶融樹脂の体積収縮率と、前記樹脂成型品の板厚方向の平均繊維配向度とを用いて２次ウエルドラインの高さを演算する２次ウエルドライン作成手段であって、２次ウエルドラインの発生位置におけるキャビティの板厚方向の幅をＴ１ｍｏｌｄ、一般部位置におけるキャビティの板厚方向の幅をＴ２ｍｏｌｄ、２次ウエルドラインの発生位置における体積収縮率をＶ１、一般部位置における体積収縮率をＶ２、２次ウエルドラインの発生位置における充填材の板厚方向に対する平均繊維配向度をｆ１、一般部位置における充填材の板厚方向に対する平均繊維配向度をｆ２、２次ウエルドラインの突出高さをｈとしたとき、ｈを、
ｈ＝（１−Ｖ１）・Ｔ１ｍｏｌｄ・ｆ１−（１−Ｖ２）・Ｔ２ｍｏｌｄ・ｆ２
の関係式に基づき演算する２次ウエルドライン作成手段として機能させることを特徴としている。
【発明の効果】
【００２０】
請求項１の発明によれば、扁平状の充填材の配向と相関の高い平均繊維配向により充填材の板厚方向の配向度を演算するため、２次ウエルドラインの発生位置における充填材の板厚方向の配向度を精度よく演算できる。充填材の板厚方向の配向度と溶融樹脂の体積収縮率とを用いて、２次ウエルドラインの発生位置と一般部位置との溶融樹脂の収縮量を演算でき、樹脂成型品の一般部に対する２次ウエルドラインの突出高さを精度よく予測できる。しかも、一般部から突出した２次ウエルドラインの突出高さを充填材の板厚方向の配向度に基づき直接的に予測するため、品質評価の精度を高くできる。
また、樹脂成型品の２次ウエルドラインの突出高さを充填材の板厚方向の平均繊維配向度と溶融樹脂の体積収縮率とを用いて精度よく演算することができる。
【００２１】
請求項２の発明によれば、溶融樹脂の流動方向の平均繊維配向度と溶融樹脂の平均移動速度とに基づいて樹脂会合面の発生位置を演算するため、１次ウエルドライン発生位置から内部流動層のウエルド境界面の移動を演算でき、実際に発生する２次ウエルドラインに対して解析結果としての２次ウエルドラインの発生位置を高い精度で予測できる。しかも、２次ウエルドラインの発生位置の演算精度が増加するため、突出高さ演算の精度を増すことができる。
【００２２】
請求項３の発明によれば、２次ウエルドライン発生位置の少なくとも一部の周辺の解析メッシュ単位を小さい解析メッシュ単位に変更するため、流動解析速度の低下を抑制しつつ、一般部から突出した２次ウエルドラインの突出高さを精度よく演算できる。
【００２３】
請求項４の発明によれば、樹脂成型品の流動解析装置において請求項１と略同様の効果を奏する。
請求項５の発明によれば、樹脂成型品の流動解析装置のコンピュータに流動解析プログラムを実行させることにより、請求項１と略同様の効果を奏する。
【図面の簡単な説明】
【００２４】
【図１】本発明の実施例１の樹脂成形品に係るバンパの斜視図である。
【図２】本発明の実施例１に係る流動解析装置のブロック図である。
【図３】メッシュモデルを示す図である。
【図４】１次ウエルドライン作成工程の説明図である。
【図５】２次ウエルドライン作成工程前半の説明図である。
【図６】２次ウエルドライン作成工程後半の説明図である。
【図７】リメッシュ工程の説明図である。
【図８】平均繊維配向画面を示す図である。
【図９】体積収縮率画面を示す図である。
【図１０】２次ウエルドライン突出高さとタルクとの関係を示す図である。
【図１１】流動解析処理の処理手順を示すフローチャートである。
【図１２】２次ウエルド外観評価点と２次ウエルドライン突出高さとの相関関係を示すグラフである。
【図１３】２次ウエルド指数差と２次ウエルドライン突出高さとの相関関係を示すグラフである。
【図１４】溶融樹脂に係る温度と比容積との相関関係を示す図である。
【図１５】１次ウエルドラインと２次ウエルドラインとの説明図である。
【発明を実施するための形態】
【００２５】
以下、本発明を実施するための形態について実施例に基づいて説明する。
【実施例１】
【００２６】
以下、本発明の実施例について図１〜図１４に基づいて説明する。
本実施例では、図１に示すように、一方向に長く延びた平板状の部材、例えば自動車前部のバンパＷを扁平状のタルク（滑石）を含有した合成樹脂材料を用いて射出成形法により製造する場合の流動解析方法について説明する。成形型ＭのキャビティＣ（図１０参照）には、複数のゲートＧ１〜Ｇ７等が設けられ、溶融樹脂が各ゲートからキャビティＣに射出される。尚、ゲート数は限定されるものではなく、樹脂成型品の形状や寸法等に応じて適宜設定することができる。
【００２７】
図２に示すように、流動解析装置１は、例えばコンピュータシステムによって構成され、入力装置２と、出力装置３と、コンピュータ本体４等を備えている。流動解析装置１は、ＣＤ−ＲＯＭ、ＩＣカード等の記録媒体に情報を記録若しくは再生する記録再生装置や、他のコンピュータシステム等とＬＡＮ等の通信装置を介してプログラムやデータを送受信することも可能である。
【００２８】
入力装置２は、キーボードやマウス等を備え、情報を入力可能に構成されている。出力装置３は、ディスプレイやＣＲＴ等の表示装置、プリンタ等の印刷装置等を備え、各種情報を出力可能に構成されている。コンピュータ本体４は、ＣＰＵ５と、ＲＯＭ６と、ＲＡＭ７と、入力インターフェース８と、出力インターフェース９等を有し、これらは通信バスにより送受信可能に接続されている。
【００２９】
ＲＯＭ６は、半導体メモリ及び大容量のハードディスク装置等によって構成され、ＣＰＵ５を駆動制御するメイン制御プログラム、出力装置３に各種情報を表示させる表示制御プログラムを格納している。更に、ＲＯＭ６は、流動解析用プログラム、１次ウエルドライン作成用プログラム、２次ウエルドライン作成用プログラムを格納している。
ＲＡＭ７には、流動解析結果、流動解析に必要な各種パラメータ、射出圧力等の成形条件、流動解析実行中の変数や各種グラフ等を一時的に記憶している。
【００３０】
ＣＰＵ５は、ＲＯＭ６に格納された各種プログラムを実行し、入力装置２、出力装置３、入出力インターフェース８，９、記録再生装置及び通信装置等を制御している。ＣＰＵ５には、流動解析部１０と、１次ウエルドライン作成部１１（１次ウエルドライン作成手段）と、２次ウエルドライン作成部１２（２次ウエルドライン作成手段）等が設けられている。
【００３１】
流動解析部１０は、図３に示すように、形状データを任意の大きさの微小要素（メッシュ）に分割して、溶融樹脂の流動解析を行うためのメッシュモデル２０を設定し、解析条件に基づき時間経過毎の速度情報や充填時間情報等の解析結果を演算可能に構成されている。メッシュモデル２０は、所定の大きさの第１解析メッシュ単位である三角形状のメッシュで構成される。
【００３２】
形状データは、ＣＡＤシステム（図示略）により作成され、入力装置２から入力される。解析条件は、ゲートから射出する溶融樹脂の射出圧力、射出速度、成形型の締付け力、成形型やゲートの温度、ゲート位置等の成形条件等であり、これら解析条件は入力装置２から入力され、ＲＡＭ７に記憶される。解析結果は、ゲートから射出された溶融樹脂が各メッシュの計算点を通過するときの流速情報や溶融樹脂が各メッシュに到達するまでの到達時間情報等である。速度情報は、各メッシュの重心位置における平均移動速度と溶融樹脂の流動方向を備えた速度ベクトルｖで表される。平均移動速度は、形状データのメッシュを板厚方向に複数層の解析要素に分割し、各解析要素における移動速度の平均値である。流動解析部１０で演算された解析結果は、ＲＡＭ７に記憶され、出力装置３に画面表示される。
【００３３】
流動解析部１０は、リメッシュ部１３と、平均繊維配向情報作成部１４と、体積収縮率情報作成部１５とを備えている。リメッシュ部１３は、メッシュモデル２０の特定領域における解析メッシュ単位を第１解析メッシュ単位よりも小さな第２解析メッシュ単位の三角形状のメッシュにリメッシュし、特定領域（リメッシュ領域２１）の流動解析を行う（図７参照）。これにより、多大な解析時間をかけることなく、特定領域（リメッシュ領域２１）の詳細な流動解析を行うことができる。
【００３４】
平均繊維配向情報作成部１４は、後述する平均繊維配向に基づきメッシュモデル２０の指定部位におけるタルクの板厚方向の配向度（配向度Ｘ３）を表示するための画面表示用情報を作成可能に構成されている。体積収縮率情報作成部１５は、溶融樹脂の比容積（ｃｍ３／ｇ）に基づきメッシュモデル２０の指定部位における体積収縮率を表示するための画面表示用情報を作成可能に構成されている。図１４に示すように、樹脂材料は、溶融樹脂に作用する圧力に応じて温度と比容積との相関関係が一義的に設定されるため、溶融樹脂の温度とその圧力を設定することで樹脂密度を演算することができる。それ故、樹脂の種類毎に樹脂温度と比容積との相関関係グラフをＲＡＭ７に格納している。平均繊維配向情報作成部１４と体積収縮率情報作成部１５とにより作成された画面表示用情報は、ＲＡＭ７に記憶され、出力装置３によりタルクの板厚方向の配向度と体積収縮率が画面表示される。
【００３５】
１次ウエルドライン作成部１１は、成形工程初期に発生する１次ウエルドラインの発生位置を演算可能に構成されている。図４に示すように、１次ウエルドライン作成部１１は、流動解析部１０により演算された速度情報に基づき、対向方向から流動してきた異なる２つ以上の溶融樹脂流動層の速度ベクトルｖが最初に会合（衝突）した樹脂会合位置（１次ウエルド発生起点）、及び前記樹脂会合部に連続して各溶融樹脂流動層の速度ベクトルｖが会合した後続位置を、最初に会合した時点からの経過時間に応じて演算し、後続の会合位置の座標を算出している。更に、１次ウエルドライン作成部１１は、１次ウエルド発生起点ＬＳ１と演算された後続の会合位置を結び１次ウエルドラインＬ１を設定し、１次ウエルドライン発生位置の座標をＲＡＭ７に記憶させる。１次ウエルドライン作成部１１により演算された１次ウエルドライン発生位置は、出力装置３により形状データ上に画面表示される。
【００３６】
２次ウエルドライン作成部１２は、樹脂成型品の流動方向の平均繊維配向と溶融樹脂の平均移動速度とに基づいて２次ウエルドラインＬ２の発生位置を演算可能に構成されている。２次ウエルドライン作成部１２は、平均繊維配向演算部１６と、２次ウエルド指数差演算部１７とを備えている。尚、平均繊維配向は、形状データのメッシュを板厚方向に複数層の解析要素に分割し、各解析要素における平均繊維配向率の平均値である。
【００３７】
異なる２つ以上の溶融樹脂流動層が会合した場合、樹脂会合部は初期に会合した１次ウエルドライン発生位置から内部流動層の圧力差に応じて流動層の圧力が均衡する位置（２次ウエルドライン発生位置）に移動する。各内部流動層の圧力が均衡した場合、溶融樹脂流動層の速度ベクトルｖは合流して進路を大きく変更する。それ故、図５に示すように、２次ウエルドライン作成部１２では、１次ウエルド発生起点ＬＳ１の速度ベクトルｖの会合部を追跡し、夫々の溶融樹脂流動層の速度ベクトルｖが合流して進行方向が安定した位置を２次ウエルド発生起点ＬＳ２として抽出し、２次ウエルド発生起点ＬＳ２の座標を演算している。
【００３８】
平均繊維配向演算部１６は、溶融樹脂に含有されたタルクの板厚方向に対する配向度（板厚方向に対する配向傾向）を演算可能に構成されている。本実施例では、タルクの配向傾向とガラス繊維等の配向傾向とが近似できることに着目し、ガラス繊維等の流動方向の平均繊維配向度をタルクの流動方向の配向度、ガラス繊維等の板厚方向の平均繊維配向度をタルクの板厚方向の配向度に代用している。
【００３９】
次式（１）で表すことができる公知の三次元メッシュの繊維配向式、例えばFolgar-Tucker繊維配向式等により平均繊維配向を演算し、この平均繊維配向の演算結果に基づき、タルクの各配向度を予測している。
【数１】
0005573276.tif 000002
［ｆ］は固有値と固有ベクトルで表される繊維配向式、Ｘ１は最大主配向方向の平均繊維配向度（０<Ｘ１<１）、Ｘ２は最大主配向方向に対して垂直方向の平均繊維配向度（０<Ｘ２<１）、Ｘ３は板厚方向の平均繊維配向度（０<Ｘ３<１）、ｘ，ｙ，ｚはタルク重心の絶対座標である。
【００４０】
２次ウエルドライン作成部１２は、流動方向の配向度Ｘ１に基づき２次ウエルドライン発生位置を演算可能に構成されている。つまり、図６に示すように、メッシュモデル２０上に演算された配向度Ｘ１を表示し、２次ウエルド発生起点ＬＳ２から配向度Ｘ１を追跡して２次ウエルドライン発生位置を抽出し、その座標をＲＡＭ７に記憶させる。２次ウエルドライン作成部１２により演算された２次ウエルドライン発生位置は、出力装置３により形状データ上に画面表示される。
【００４１】
２次ウエルド指数差演算部１７は、溶融樹脂の体積収縮率と配向度Ｘ３とを用いて一般部（２次ウエルドライン以外の領域）に対する２次ウエルドラインの突出高さ（２次ウエルド指数差）を演算可能に構成されている。成形工程初期に樹脂会合部において板厚方向に配向したタルクは、後続の内部流動層の移動に伴い１次ウエルドライン発生位置から２次ウエルドライン発生位置へ移動する。このとき、溶融樹脂の収縮後の樹脂成型品の板厚Ｔは、次式（２）により表すことができる。
【数２】
0005573276.tif 000003
尚、Ｔｍｏｌｄは成形型ＭのキャビティＣの板厚方向の幅、Ｖｒは溶融樹脂の板厚方向の樹脂収縮率、ａはタルク配向度に対する収縮係数、Ｆは板厚方向のタルク配向度である。
【００４２】
図１０に示すように、２次ウエルドラインの発生位置における板厚をＴ１、一般部位置における板厚をＴ２、２次ウエルドラインの発生位置における板厚方向の樹脂収縮率をＶｒ１、一般部位置における板厚方向の樹脂収縮率をＶｒ２、２次ウエルドラインの発生位置における板厚方向のタルク配向度をＦ１、一般部位置における板厚方向のタルク配向度をＦ２、２次ウエルドラインの発生位置におけるキャビティの板厚方向の幅をＴ１ｍｏｌｄ、一般部位置におけるキャビティの板厚方向の幅をＴ２ｍｏｌｄ、２次ウエルド指数差をｈとしたとき、次式（３）により表すことができる。
【数３】
0005573276.tif 000004
【００４３】
ここで、板厚方向の樹脂収縮率Ｖｒは樹脂が三次元的に収縮することから算出が難しく、また、タルク配向度Ｆは現時点適した計算式が存在しないため、次式（４）に示すように、樹脂収縮率Ｖｒに代えて体積収縮率Ｖ、タルク配向度Ｆに代えて平均繊維配向の板厚方向の平均繊維配向度ｆを夫々代用している。
【数４】
0005573276.tif 000005
【００４４】
式（４）の係数を省略した近似式による演算結果と試作試験等による実測値とを相関分析した結果、高い相関関係が判明したため、２次ウエルド指数差演算部１７では、次式（５）を用いて２次ウエルド指数差を演算している。この式（５）や三次元メッシュの繊維配向式（１）等はＲＡＭ７に予め記憶されている。
【数５】
0005573276.tif 000006
【００４５】
図１１のフローチャートに基づき、本流動解析方法の処理手順を説明する。尚、Ｓｉ（ｉ＝１，２，…）は各処理ステップを示す。
条件取得工程では、ＣＰＵ５が記憶装置５から溶融樹脂の射出圧力や射出速度等の解析条件、繊維配向式（１）、ウエルド指数差式（５）等本流動解析に必要な情報を取得する（Ｓ１）。
【００４６】
メッシュモデル設定行程では、バンパＷの形状データに基づき作成されたメッシュモデル２０を設定する（Ｓ２）。メッシュモデル２０は、第１メッシュ単位の三角形状のメッシュで分割されている（図３参照）。ＣＰＵ５では、Ｓ１で取得された解析条件を用いて、経過時間に応じた溶融樹脂流動の挙動を解析し、各流動層の流速情報や到達時間情報等の解析結果の履歴を取得する。
【００４７】
１次ウエルドライン作成工程では、Ｓ２で取得された解析結果の履歴に基づき、バンパＷに発生が予測される１次ウエルドラインＬ１を作成する（Ｓ３）。各流動層の速度ベクトルｖの挙動に基づき１次ウエルド発生起点ＬＳ１と１次ウエルドラインＬ１を設定し、１次ウエルドライン発生位置の座標を取得する（図４参照）。
【００４８】
平均繊維配向演算工程では、溶融樹脂に含有されたタルクについて平均繊維配向を式（１）を用いて演算し、配向度Ｘ１、配向度Ｘ２、配向度Ｘ３を算出する（Ｓ４）。
２次ウエルドライン作成工程では、１次ウエルド発生起点ＬＳ１における速度ベクトルｖに基づき樹脂会合部の移動を追跡して２次ウエルド発生起点ＬＳ２を抽出する（Ｓ５）（図５参照）。メッシュモデル２０上に配向度Ｘ１を表示し、２次ウエルド発生起点ＬＳ２から配向度Ｘ１を追跡して２次ウエルドラインＬ２を抽出し、２次ウエルドライン発生位置の座標を取得する（図６参照）。
【００４９】
次に、２次ウエルドライン発生位置の少なくとも一部の周辺領域についてリメッシュ領域２１を作成済みか否か判定し、リメッシュ領域２１を作成していない場合、リメッシュ工程（Ｓ６）を行う。Ｓ６では、図７に示すように、Ｓ５で取得された２次ウエルドライン発生位置の座標に基づき、メッシュモデル２０において２次ウエルドラインＬ２の周辺領域を第２メッシュ単位の三角形状のメッシュで分割したリメッシュ領域２１を設定し、リメッシュ領域２１を用いて平均繊維配向について再度解析を行う。リメッシュ後の解析結果は速度ベクトルの表示が小さく確認が困難であるため、本実施例では、２次ウエルドライン発生位置の座標はリメッシュ工程前の解析結果を用いている。尚、２次ウエルドライン発生位置の座標について、再度解析を行うことも可能である。
【００５０】
２次ウエルド指数差演算工程では、一般部位置に対する２次ウエルドラインの突出高さに相当する２次ウエルド指数差を演算する（Ｓ７）。
まず、２次ウエルドラインＬ２と２次ウエルドラインＬ２から所定距離離隔した一般部の位置について、夫々の板厚方向の平均繊維配向度と体積収縮率とを設定する。
図８に示すように、タルクの配向度Ｘ３に基づく平均繊維配向画面を表示し、２次ウエルドラインＬ２上の点Ｐ１における配向度Ｘ３と点Ｐ１から流動方向に対して垂直方向に所定距離離隔した一般部位置の点Ｐ２，Ｐ３における配向度Ｘ３を取得する。同様に、図９に示すように、体積収縮率に基づく体積収縮率画面を表示し、２次ウエルドラインＬ２上の点Ｐ１における体積収縮率と点Ｐ１から流動方向に対して垂直方向に所定距離離隔した一般部位置の点Ｐ２，Ｐ３における体積収縮率を取得する。尚、平均繊維配向画面や体積収縮率画面を表示する際、２次ウエルドラインＬ２を同じ画面上に表示することも可能である。
【００５１】
次に、取得された各位置の配向度Ｘ３及び体積収縮率を式（５）に代入して一般部位置に対する２次ウエルドラインＬ２の突出高さに相当する２次ウエルド指数差を演算する。この演算結果に基づき、２次ウエルド指数差が所定値よりも大きな場合、樹脂成型品の外観への影響が大きいと判定している。
【００５２】
図１２，図１３に基づき、実際の２次ウエルドラインの突出高さと２次ウエルド指数差との相関関係を説明する。試作試験によりタルクを含有する樹脂成形品のサンプルを複数製作し、これらサンプルについて、２次ウエルドラインの突出高さの実測値、２次ウエルド外観評価点数及び２次ウエルド指数差を求めた。
【００５３】
図１２に示すように、各サンプルについて、横軸を２次ウエルド外観評価点、縦軸を２次ウエルドラインの突出高さの実測値としたグラフを作成した。２次ウエルド外観評価点と２次ウエルドラインの突出高さとは一次関数の関係式で表すことができる。また、２次ウエルド外観評価点がα以上のとき、つまり、２次ウエルドラインの突出高さがβ以下のとき、樹脂成型品の外観は良好であった。
【００５４】
次に、図１３に示すように、各サンプルについて、横軸を２次ウエルド指数差、縦軸を２次ウエルドラインの突出高さの実測値としたグラフを作成した。２次ウエルド指数差と２次ウエルドラインの突出高さとは一次関数の関係式で表すことができる。以上により、２次ウエルド指数差と２次ウエルド外観評価点とは相関することが分かる。それ故、２次ウエルドラインの突出高さがβ以下のとき、つまり、２次ウエルド指数差がγ以下のとき、樹脂成型品の外観は良好と判定することができる。
【００５５】
次に、本流動解析方法、流動解析装置１及び流動解析プログラムの作用、効果について説明する。
この本流動解析方法は、メッシュモデル２０を用いて成形型のキャビティ内の溶融樹脂を流動解析して成形品表面から微小量突出する２次ウエルドラインＬ２の発生を予測する樹脂成型品の流動解析方法であって、成形工程初期に発生する１次ウエルドラインＬ１の発生位置を演算する１次ウエルドライン作成工程と、前記演算された１次ウエルドライン発生位置に基づいて、前記樹脂成型品Ｗの板厚方向の平均繊維配向度Ｘ３と、前記溶融樹脂の体積収縮率Ｖとを用いて２次ウエルドラインＬ２の突出高さｈを演算する２次ウエルド指数差演算工程と、を備えている。
【００５６】
この本流動解析方法によれば、扁平状のタルク配向と相関の高い平均繊維配向によりタルクの板厚方向の配向度Ｘ３を演算するため、２次ウエルドラインＬ２の発生位置におけるタルクの板厚方向の配向度を精度よく演算できる。配向度Ｘ３と溶融樹脂の体積収縮率Ｖとを用いて、２次ウエルドラインＬ２の発生位置と一般部位置との溶融樹脂の収縮量を演算でき、樹脂成型品の一般部位置に対する２次ウエルドラインの突出高さｈを精度よく予測できる。しかも、一般部から突出した２次ウエルドラインＬ２の突出高さをタルクの板厚方向の配向度に基づき直接的に予測するため、品質評価の精度を高くできる。
【００５７】
前記１次ウエルドライン作成工程と２次ウエルド指数差演算工程との間において、平均繊維配向から求めた流動方向の配向度Ｘ１と前記溶融樹脂の平均移動速度とに基づいて２次ウエルドラインの発生位置を演算する２次ウエルドライン作成工程を備えるため、１次ウエルドライン発生位置から内部流動層の樹脂会合部の移動を演算でき、実際に発生する２次ウエルドラインＬ２に対して２次ウエルドラインＬ２の発生位置を高い精度で予測できる。しかも、２次ウエルドラインＬ２の発生位置の演算精度が増加するため、突出高さｈの演算精度を増すことができる。
【００５８】
２次ウエルドライン作成工程と２次ウエルド指数差演算工程との間において、前記２次ウエルドライン発生位置の少なくとも一部の周辺の解析メッシュ単位を小さい解析メッシュ単位に変更するリメッシュ領域設定工程を備えたため、流動解析速度の低下を抑制しつつ、一般部から突出した２次ウエルドラインＬ２の突出高さｈを精度よく演算できる。
【００５９】
前記溶融樹脂はタルクを含有し、２次ウエルドラインＬ２の発生位置の点Ｐ１におけるキャビティの板厚方向の幅をＴ１ｍｏｌｄ、一般部位置の点Ｐ２（又はＰ３）におけるキャビティの板厚方向の幅をＴ２ｍｏｌｄ、２次ウエルドラインの発生位置の点Ｐ１における体積収縮率をＶ１、一般部位置の点Ｐ２（又はＰ３）における体積収縮率をＶ２、２次ウエルドラインの発生位置の点Ｐ１におけるタルクの板厚方向に対する平均繊維配向度をｆ１、一般部位置の点Ｐ２（又はＰ３）におけるタルクの板厚方向に対する平均繊維配向度をｆ２、２次ウエルドラインの突出高さをｈとしたとき、ｈを、  
ｈ＝（１−Ｖ１）・Ｔ１ｍｏｌｄ・ｆ１−（１−Ｖ２）・Ｔ２ｍｏｌｄ・ｆ２
の関係式に基づき演算するため、樹脂成型品Ｗの２次ウエルドラインＬ２の突出高さｈをタルクの板厚方向の平均繊維配向度と溶融樹脂の体積収縮率とを用いて精度よく演算することができる。
【００６０】
この本流動解析装置１は、メッシュモデル２０により成形型内の溶融樹脂を流動解析して成形品表面から突出する２次ウエルドラインＬ２の発生を予測する樹脂成型品の流動解析装置であって、成形工程初期に発生する１次ウエルドラインＬ１の発生位置を演算する１次ウエルドライン作成部１１と、前記演算された１次ウエルドライン発生位置に基づき前記溶融樹脂の体積収縮率Ｖと、前記樹脂成型品の板厚方向の平均繊維配向度Ｘ３とを用いて２次ウエルドラインの高さｈを演算する２次ウエルドライン作成部１２と、を備えているため、前述と略同様の効果を奏することができる。
【００６１】
記憶装置５に記憶された制御プログラムは、メッシュモデル２０により成形型内の溶融樹脂を流動解析して成形品表面から突出する２次ウエルドラインＬ２の発生を予測する樹脂成型品の流動解析装置１のコンピュータ本体４を、成形工程初期に発生する１次ウエルドラインＬ１の発生位置を演算する１次ウエルドライン作成部１１と、前記演算された１次ウエルドライン発生位置に基づき前記溶融樹脂の体積収縮率Ｖと、前記樹脂成型品の板厚方向の平均繊維配向度Ｘ３とを用いて２次ウエルドラインＬ２の高さｈを演算する２次ウエルドライン作成部１２として機能させるため、前述と略同様の効果を奏することができる。
【００６２】
次に、前記実施例を部分的に変更した変形例について説明する。
１〕前記実施例においては、樹脂に含有した充填材にタルクを用いた例を説明したが、少なくとも、溶融樹脂の流動に伴い配向する傾向を有する楕円状や扁平状の充填材を含有する樹脂であれば何れの樹脂成形品の流動解析にも適用可能である。また、樹脂成形品は自動車のバンパに限ることなく、板状の成形品であれば何れの樹脂成形品の流動解析にも適用可能である。
【００６３】
２〕前記実施例においては、２次ウエルド指数差演算工程の前にリメッシュ工程を行い再度解析する例を説明したが、初期のメッシュモデルのメッシュ単位が解析精度に適合している場合は、リメッシュ工程を省略することも可能である。このとき、最初に演算した板厚方向の平均繊維配向度等の演算結果を２次ウエルド指数差演算工程で使用する。
【００６４】
３〕前記実施例においては、一般部に対する２次ウエルドラインの突出高さである２次ウエルド指数差を予測する例を説明したが、単に２次ウエルドラインの突出高さである２次ウエルド指数を予測することも可能であり、２次ウエルド指数により品質評価を行うことも可能である。
【００６５】
４〕前記実施例においては、流動解析装置が１次ウエルド発生起点における速度ベクトルに基づき樹脂会合部の移動を追跡して２次ウエルド発生起点ＬＳ２を抽出する例を説明したが、画面上の速度ベクトルに基づき人手により２次ウエルド発生起点ＬＳ２を抽出することも可能である。
【００６６】
５〕その他、当業者であれば、本発明の趣旨を逸脱することなく、前記実施例に種々の変更を付加した形態で実施可能であり、本発明はそのような変更形態も包含するものである。
【産業上の利用可能性】
【００６７】
本発明は、樹脂成型品の流動解析において、樹脂成型品の板厚方向の平均繊維配向度と樹脂の体積収縮率とを用いて２次ウエルドラインの突出高さを演算することにより、扁平状の充填材を含有した樹脂でも、高い精度の品質評価を可能にできる。
【符号の説明】
【００６８】
１          流動解析装置
４          コンピュータ
１０        流動解析部
１１        １次ウエルドライン作成部
１２        ２次ウエルドライン作成部
１３        リメッシュ部
１６        平均繊維配向演算部
１７        ２次ウエルド指数演算部
２０        メッシュモデル
Ｌ１        １次ウエルドライン
Ｌ２        ２次ウエルドライン
Ｘ１        （最大主配向方向）平均繊維配向度
Ｘ３        （板厚方向）平均繊維配向度
【図１】
0005573276.tif 000007
【図２】
0005573276.tif 000008
【図３】
0005573276.tif 000009
【図４】
0005573276.tif 000010
【図５】
0005573276.tif 000011
【図６】
0005573276.tif 000012
【図７】
0005573276.tif 000013
【図８】
0005573276.tif 000014
【図９】
0005573276.tif 000015
【図１０】
0005573276.tif 000016
【図１１】
0005573276.tif 000017
【図１２】
0005573276.tif 000018
【図１３】
0005573276.tif 000019
【図１４】
0005573276.tif 000020
【図１５】
0005573276.tif 000021