(57)【要約】【課題】発泡流路構造の適正化および発泡材料注入量の適正化を図るための発泡流動挙動を解析できる計算方法、計算プログラムを開発する。【解決手段】密度を、流体が注入口を通過してからの経過時間項および肉厚項を含む関数として入力して、ある時間において計算された密度を連続の式および運動方程式に代入することにより流動速度および圧力を求める。ここで、粘度は時間項を含む関数とする。このとき、計算の収束を判定し、収束しない場合は、境界条件および成形条件を修正する。また、計算が収束した場合には、密度および圧力分布の判定を行い、密度および圧力分布が設計許容範囲からはずれる場合には、モデル形状、物性値または成形条件を修正する計算手法を用いた解析プログラムによって、発泡流路構造の適正化および発泡材料注入量の適正化を図る。【選択図】図３
2006018616.tif 000002
【特許請求の範囲】
【請求項１】
解析対象製品の発泡充填部分のデータを記憶装置から読み込み、当該データに基づいて３次元ソリッド要素に分解処理をし、注入口から注入された発泡材料の密度変化を伴う流動過程を計算する際に、
密度が時間変化を伴う関数として表される連続の式、運動方程式を、前記３次元ソリッド要素に基づいて演算処理し、
注入口を通過する発泡材料に時間の履歴を残し、前記の２式に含まれる前記発泡材料の密度が注入口を通過してからの経過時間項を含む関数（密度の関数式）として入力され、各タイムステップ毎に前記密度の関数式に基づいて算出された密度を、前記連続の式および運動方程式に代入して、少なくとも密度の演算出力を行なうことを特徴とする有限要素解析法。
【請求項２】
解析対象製品の発泡充填部分のデータを記憶装置から読み込み、当該データに基づいて３次元ソリッド要素に分解処理をし、注入口から注入された発泡材料の粘度変化を伴う流動過程を計算する際に、
粘度が時間変化を伴う関数として表される連続の式、運動方程式を、前記３次元ソリッド要素に基づいて演算処理し、
注入口を通過する発泡材料に時間の履歴を残し、前記の２式に含まれる前記発泡材料の粘度が注入口を通過してからの経過時間項を含む関数（粘度の関数式）として入力され、各タイムステップ毎に前記粘度の関数式に基づいて算出された粘度を、前記連続の式および運動方程式に代入して、少なくとも粘度の演算出力を行なうことを特徴とする有限要素解析法。
【請求項３】
請求項１記載の解析方法であって、発泡材料を注入している時間においては、注入口における密度が一定値であり、前記発泡材料が注入口を通過してからの経過時間を用いて計算された密度分布の演算出力を行うことを特徴とした有限要素解析法。
【請求項４】
請求項１〜２いずれか記載の解析方法であって、注入口を通過してからの経過時間分布の演算出力を行うことを特長した有限要素解析法。
【請求項５】
請求項２記載の解析方法であって、発泡材料を注入している時間においては、注入口における粘度が一定値であり、前記発泡材料が注入口を通過してからの経過時間を用いて計算された粘度分布の演算出力を行うことを特徴とした有限要素解析法。
【請求項６】
請求項１〜２記載の密度または粘度の関数式が、時間と肉厚を含む関数として表され、任意の座標の肉厚が数値として入力され、粘度または密度の計算を行うことを特徴とする有限要素解析法。
【請求項７】
請求項１〜２記載の密度または粘度の関数式が、時間または肉厚を含む関数として表され、前記肉厚が他のCADソフトなどで算出された値として入力され、粘度または密度の計算を行うことを特徴とする有限要素解析法。
【請求項８】
請求項１〜２記載の解析に入力される粘度が、せん断速度または密度を含む関数として表され、粘度または密度の計算を行うことを特徴とする有限要素解析法。
【請求項９】
請求項１〜８いずれか記載の解析法を用いたことを特徴とする有限差分３次元発泡流動解析プログラムまたは有限要素３次元発泡流動解析プログラムを格納したことを特徴とする記録媒体。
【請求項１０】
請求項９記載の３次元発泡流動解析プログラムによって得られる密度を含む結果を用いて、発泡材料の原液を注入する位置、最終充填場所、注入量、または発泡樹脂流路の部品配置、または成形条件、または熱伝導率などの発泡材料の物性値を決めたことを特徴とする発泡製品の製造方法。
【発明の詳細な説明】
【背景技術】
【０００１】
本発明は、プラスチック発泡成形加工技術に係り、冷蔵庫や建設材料、自動車等に使用される断熱材、衝撃吸収材等の発泡部品を発泡成形する際の３次元発泡流動解析方法に関する。
時間を変数とする発泡材料の密度の関数を導入し、密度が時間と共に減少する発泡挙動を解析可能な解析プログラムが、特許文献１（特開2001-318909号公報）、特許文献２（特開2003-91561号公報）に記載されている。
【０００２】
特に、特許文献１では、発泡材料全体を均一の密度とみなし、さらに、その密度として、発泡原料を攪拌した発泡材料を出すノズルを最初に出た発泡材料のノズルを出てからの経過時間で算出した密度を用いている。
【０００３】
また、特許文献２では、特許文献１の技術に加え、肉厚の変動によって発泡材料の密度が変化することを考慮した関数を用いて発泡材料の発泡流動解析を行うことが記載されている。
【０００４】
【特許文献１】特開2001-318909号公報
【特許文献２】特開2003-91561号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
実際の製造プロセスでは、ノズル近傍では実質的に一定な密度のままであり、時間とともに密度が変動することはない。  
また、最初にノズルから出た発泡材料の影響により、その後にノズルから出る発泡材料の流動性が低下する。1秒程度までの発泡流動解析であれば、近似値として有用であると考えるが、それ以上の発泡流動解析では現実との差が大きくなりすぎる。
【０００６】
かかる理由から特許文献１及び２による解析は精度が低かった。
【０００７】
本発明の目的は、発泡流動解析の精度を向上することにある。
【０００８】
また、特許文献２には、発泡樹脂の密度を区分けして表示しようとするものが記載されているが、肉厚によって区分けしているので、肉厚一定では特許文献１と同じ出力となってしまう問題がある。
【０００９】
本発明の他の目的は、ノズルから出る発泡材料の密度を一定値とすることによって、肉厚一定であっても、発泡密度の大きさによる区分けが可能な発泡流動解析を実現することにある。
【課題を解決するための手段】
【００１０】
上記課題を解決するため、本発明は、密度が注入口を通過してからの時間経過を含む関数として入力されることを特徴とし、このタイムステップ毎の密度変化を連続の式および運動方程式に代入して計算することにより、圧力、流動速度、密度分布などの結果を得ることができる３次元発泡流動解析法または発泡流動解析プログラムを提供する。なお、実際のウレタン材料などの発泡現象において、粘度も時間項を含む硬化反応を伴って変化する挙動を示すので、詳細な発泡流動解析を行うために、密度だけでなく粘度も、注入口を通過してからの時間経過を含む関数として入力し、タイムステップ毎の密度および粘度変化を連続の式、運動方程式に代入して計算することにより、圧力、流動速度、密度分布などの結果を得る３次元発泡流動解析法または発泡流動解析プログラムを提供することもできる。
【発明の効果】
【００１１】
本発明によれば、発泡流動解析システムの解析精度を向上することができる。
【発明を実施するための最良の形態】
【００１２】
以下、添付の図面を参照しながら、本発明に係る実施の形態について説明する。
【実施例１】
【００１３】
まず、解析対象となる発泡成形工程を図１を用いて説明する。
【００１４】
発泡成形工程は、発泡材料充填装置のノズルから発泡型へ充填することにより、なされる。
【００１５】
発泡材料充填装置は、発泡原料を格納した複数のタンクと、そのタンクから出る原料の量を制御するバルブと、そのバルブを介して繋がって、複数のタンクから流れ込んだ発泡原料を混合した発泡材料をノズルから出すミキシングヘッドとを有している。
【００１６】
本実施例では、シクロペンタン(C5H10)発泡剤を充填したポリオール１を一つのタンクに格納し、イソシアネート２を他のタンクに格納し、それらの２液をミキシングヘッド３によって攪拌し、発泡型４の中に２液を攪拌することにより形成した発泡材料５を充填することにより、発泡させるものとした。
【００１７】
次に、この発泡製品の発泡成形工程における発泡材料の発泡流動解析システムについて説明する。  
発泡流動解析システムは、図５示すハードウェア構成で後述する図３のフローを備えたソフトウェアが実行されることにより、機能する。
【００１８】
具体的には、計算装置６、記録装置10(ハードディスク、MOなど)を備えた計算装置７、この２つの計算装置を繋ぐＬＡＮ８、計算装置７が備える表示装置９を備えている。計算装置６で作成したCADデータをLAN８を介して計算装置７に転送するようにしておく。計算装置７に転送されたCADデータは、計算装置７の記録装置10(ハードディスク、MOなど)に記録される。計算機７は、図3で示すフローチャートに従い計算を実行し、結果を記録装置10に記録した後、表示装置９に結果を表示する。
【００１９】
図示してはいないが、計算装置６及び７には、当然キーボードやマウス等の入力デバイスを備えている。
【００２０】
次に、図３のフローチャートに沿って発泡流動解析プログラムの処理を説明する。
【００２１】
まず、モデル形状作成ステップ５０１では、オペレータによって入力装置を介して特定された解析対象モデル、つまり、発泡材料が充填される空間の形状データを記憶装置１０から読み出す。具体的には、発泡材料の流動部分の初期設計モデルのＣＡＤデータを記憶装置から読み込む。
【００２２】
次に、３次元ソリッド要素作成のステップ５０２では、モデル形状作成ステップ５０１で読み込んだデータを用いて、発泡材料が充填される断熱部分の形状を複数の特定空間（３次元ソリッドの有限要素）に分解し、各有限要素の形状データを作成する。
【００２３】
次に、流体の物性値入力ステップ５０３では、解析を行う発泡材料の物性値である密度式（５）、粘度式(６)を入力するように、オペレータに催促する表示を行い、入力装置からこれらのデータを受け付ける。
【００２４】
次に、境界条件、成形条件入力ステップ５０４において、オペレータに対して３次元ソリッド要素内に流体が注入する際の初期速度、初期密度、初期粘度、注入場所、肉厚の入力をするように、オペレータに催促する表示を行い、入力装置からこれらのデータを受け付ける。受け付けた注入場所のデータとモデル形状作成ステップ５０１で読み込んだ形状データから、注入口の断面積を算出する。
【００２５】
次に、オペレータからの解析開始の指示と初期時間増分を受け付ける。
【００２６】
ステップ５０５として、この指示に基づいて、記録装置に格納された連続の式(１)および運動方程式(２)〜(４)を呼び出し、これまで入力を受け付けた、初期時間増分、初期速度、初期密度、初期粘度、粘度式(６)、密度式(５)を代入し、発泡材料が流動していく、速度、圧力、密度および粘度を計算する。この計算結果と、初期時間増分と、注入口に最も近い有限要素の位置と対応つけて記憶装置に保存する。
【００２７】
次にステップ５０６で、ステップ５０５の計算を繰り返し行うとともに、代入していった時間増分（初期時間増分を含む）を加算していく。つまり、注入口を発泡材料が通過してからその３次元ソリッド要素に到達する経過時間を算出する。なお、初期時間増分は最初の計算にのみ使用し、２回目の計算からは収束安定性を確保できる時間増分を適宜算出し、初期時間増分の代わりに代入するようにする。
【００２８】
また、２回目の計算からは、注入口に最も近い有限要素（第１の有限要素）だけでなく、発泡材料が流動していく方向（速度で示される方向）に隣接する有限要素群（第２の有限要素）に対してもおこなう。第１の有限要素に対して計算する場合には、計算で新たに求めた時間増分を用いる。第２の有限要素に対して計算をする場合には、第１の有限要素の出力である、速度、密度、粘度を、初期速度、初期密度、初期粘度を入力値として用いる。但し、密度と粘度を計算するときには、注入口を発泡材料が通過してからその有限要素に到達する経過時間（初期時間増分＋計算で求めた時間増分）を用いて計算する。
【００２９】
２回目の出力としては、有限要素の位置（第１の有限要素と第２の有限要素）、注入口を発泡材料が通過してからその有限要素に到達する経過時間（第１の有限要素にとっては計算で求めた時間増分であり、第２の有限要素にとっては初期時間増分と計算で求めた時間増分の和）、速度、密度、粘度を対応付けて記憶しておく。１回目の出力と２回目の出力として保存しているデータテーブルを図４として示す。ここで、要素Aが第一の要素を示しており、初期時間増分Δt１において計算された速度VA1、密度ρ１、粘度ρ１を算出する。第２の有限要素Ｂ〜Ｆは、これらの値を初期値として、第2回目の計算を行う。ここで、第２の有限要素Ｂ〜Ｆの速度は各要素によって異なるが、入り口を通過してからの時間はΔt1+Δt2で等しいので、式(５)（６）により、肉厚値が等しい場合には第２の有限要素Ｂ〜Ｆの密度および粘度は等しくなる。
【００３０】
この繰り返しは、各有限要素の体積をあらかじめ記憶しておき、計算に用いた有限要素の体積の和を求め、形状データから求めたモデルの発泡体積と一致したときに終了する。
【００３１】
ステップ５０７で、収束の判定を行う。収束の判定手法は、圧力とあらかじめ定めておいた圧力範囲とを対比し、範囲内にある場合を収束として判定する。収束しない場合には、ステップ５０１〜５０４のいづれかに戻る。この際、オペレータに入力を促し、どのステップに戻るかを決める。
【００３２】
ステップ５０８では、注入口を発泡材料が通過してからその有限要素に到達する経過時間と、密度、粘度を用いて表示を行う。表示の例を図6〜9に示す。
【００３３】
発泡解析の表示例を図６に示す。これは、図２に示した形状を用いて、前記した密度式（５）および粘度式（６）を、連続の式(１)、運動方程式(２)〜(４)と連動させることにより発泡流動解析を行った肉厚中心部分のYZ平面における結果である。発泡領域12の寸法は、高さ500(mm)、幅250(mm)、肉厚40(mm)とし、注入場所１１から初期密度ρ０＝1000(kg/m3)の発泡材料５をZ軸のプラス方向に流量0.2（kg/s）で注入し続け、式（５）におけるＡ＝１とし、式（６）における初期粘度η0を１(Pa・s)、材料固有の定数F＝1.２、te＝60とし、重力はＺ軸のマイナス方向に加えた。ここでは、時間変化に伴う材料の発泡流動過程を示しており、10(s)、20(s)、30(s)における発泡材料５の充填量を示している。このように、発泡領域12に注入された発泡材料５が、Ｚ方向のプラス方向に発泡流動している過程が解析できている。
【００３４】
また、この解析による出力結果の一例として、30(s)後における注入口を通過してからの経過時間の分布を図７に示す。図中の線が、経過時間が等しい位置を示しており、入り口を通過している流体の経過時間が0(ｓ)、Z軸プラス方向の最初に注入された流体の経過時間が30(ｓ)の値を示している。更に、30(s)後における発泡密度の分布を図８に示す。発泡密度は式（５）に示すように、注入口を通過してからの経過時間の関数となっている。ここで、密度が等しい位置を図中の線で示しており、入り口における密度が初期密度1000(kg/m3)で最も高く、最初に注入された発泡材料の先端(Z軸プラス方向)の密度が最も低くなる分布を示している。また、30(s)後における粘度の分布を図９に示す。粘度は式（６）に示すように、注入口を通過してからの経過時間の関数となっている。ここで、粘度が等しい位置を図中の線で示しており、入り口における粘度が初期粘度1(Pa・s)で最も低く、最初に注入された発泡材料の先端(Z軸プラス方向)の粘度が最も高くなる分布を示している。
【００３５】
また、以上示した解析により、発泡材料５のやわらかさを評価するため、弾性率を密度の式として入力し、求められた弾性率から発泡材料５のやわらかさを求めることもできる。
【００３６】
また、上述した、ステップ５０３で入力する粘度に関して、式（６）で示した時間の関数の他に、式（５）で示す密度式において、発泡後の最終密度をρminとし、この最終密度を密度の最小値とすると、粘度は式（７）で表すことが出来る。ここで、η0：初期温度、γ：せん断速度、ｃ、ｎ：材料固有の係数を示す。
η=（η0 (ρ（ｔ）/ρmin)C  ）／γｎ        (７)
このように、粘度は発泡後の最終密度を基準とした発泡の進行率と、せん断速度の関係として表すことが出来る。但し、本解析に用いる粘度は式（７）で示される関数だけではなく、時間項を含む関数など任意の値を用いることが出来る。
【００３７】
ここで、上述した式１から式６について説明する。
【００３８】
一般的に、ウレタン発泡成形における初期密度から発泡後の最終密度との倍率を示す発泡倍率に影響を与える主因子は治具温度と肉厚であることが知られており、この治具温度と肉厚を考慮した発泡挙動を解析するためには、解析モデルを３次元ソリッド要素に分割し、連続の式(１)、運動方程式(２)〜(４)によって、流動速度、圧力、温度を求める手法を用いる。ここで、ρ；密度、ｕ；ｘ方向速度、υ；ｙ方向速度、ω；ｚ方向速度、Ｐ；圧力、ｔ；時間、η；粘度、を示している。
【００３９】
これらの式をソフトウェア上で実現した場合には、「ρの時間変化を含む関数」、「粘度η」は入力値として使用される。
【００４０】
【数１】
2006018616.tif 000003
【００４１】
【数２】
2006018616.tif 000004
【００４２】
【数３】
2006018616.tif 000005
【００４３】
【数４】
2006018616.tif 000006
【００４４】
ここで、密度変化のない液体を扱う場合には、∂ρ/∂t=0となる。しかし、冷蔵庫ウレタン材料等の発泡材料を扱う場合には、∂ρ/∂t≠0となり、前記した発泡倍率に影響を与える主因子である治具温度と肉厚を考慮した発泡材料の密度変化を別途求めることが必要となる。
【００４５】
このとき、一つの治具の場所による温度バラツキおよび発泡成形毎の治具温度のバラツキを無視すると、可視化実験などの結果を用いて、ある治具温度における密度の変化を数式化することにより、密度が時間項および肉厚項を含む関数として表記できる。図２に発泡型4に囲まれた発泡領域12の断面図を示す。ここで、注入口11から流量Ｑで初期密度ρ０、初期粘度η0の発泡材料5を一定時間流入した場合の密度変化関数を式（５）として、粘度変化関数を式（６）として示す。ここで、ｔ：発泡材料が入り口を通過してからの経過時間、Ｈ；製品肉厚、ρ0；初期密度、η0；初期粘度、Ａ、F、te；材料固有の定数とする。
ρ(ｔ)＝ρ0／(Ａ×ｔ2／Ｈ＋１)                  (５)
η = η0 ( ( 1+ｔ／te )／( 1-ｔ／te ) )F          (６)            
この式（５）、(６)によって求められる各タイムステップにおける密度および粘度を、式(１)〜(４)に代入することにより、発泡時の流動速度、圧力、温度などの結果を求めることができる。この密度式（５）を用いると、材料が流動しながら発泡する過程を解析できる。
【００４６】
ここで、肉厚Ｈは、肉厚を構成するある壁面からもう一面の壁面への垂線として計算されるか、もしくは座標ごとの肉厚値として任意に入力することが出来るものとする。また、肉厚の入力方法として、他のCADソフトなどで算出された値を、各要素ごとの肉厚値に変換して入力することもできる。
【００４７】
なお、式（５）（６）において、密度および粘度は時間の関数として表したが、本発明は、これだけに限定するものではなく、使用する発泡材料の特性に応じて、時間項および肉厚項を含む任意の関数(圧力、粘度などの変数も考慮に入れた関数としても表すことができるものとする。) に対応した式で表すことができるものとする。
【図面の簡単な説明】
【００４８】
【図１】ポリオールとイソシアネートの２液混合による発泡過程
【図２】解析検討に用いた発泡型に囲まれた発泡領域の断面図
【図３】発泡材料の注入口を通過してからの経過時間項および肉厚項を含む関数を入力する場合の発泡流動解析フローチャート
【図４】時間増分ごとの各要素における計算項目
【図５】発泡解析を行うハードウェア構成図
【図６】発泡解析結果（時間変化に伴う材料の発泡流動過程）
【図７】発泡解析結果（発泡材料の注入口を通過してからの経過時間分布）
【図８】発泡解析結果（密度分布）
【図９】発泡解析結果（粘度分布）
【符号の説明】
【００４９】
１…ポリオール  ２…イソシアネート  ３…ミキシングヘッド  ４…発泡型
５…発泡材料    ６…計算装置  ７…計算装置  ８…LAN
９…表示装置    １０…記録装置    １１…注入口    １２…発泡領域  
【図１】
2006018616.tif 000007
【図２】
2006018616.tif 000008
【図３】
2006018616.tif 000009
【図４】
2006018616.tif 000010
【図５】
2006018616.tif 000011
【図６】
2006018616.tif 000012
【図７】
2006018616.tif 000013
【図８】
2006018616.tif 000014
【図９】
2006018616.tif 000015