JP3994709B2 - 発泡材料の3次元発泡解析方法及びそれを記録した記録媒体並びにそれを用いた発泡材料の断熱性評価方法 - Google Patents
発泡材料の3次元発泡解析方法及びそれを記録した記録媒体並びにそれを用いた発泡材料の断熱性評価方法 Download PDFInfo
- Publication number
- JP3994709B2 JP3994709B2 JP2001282549A JP2001282549A JP3994709B2 JP 3994709 B2 JP3994709 B2 JP 3994709B2 JP 2001282549 A JP2001282549 A JP 2001282549A JP 2001282549 A JP2001282549 A JP 2001282549A JP 3994709 B2 JP3994709 B2 JP 3994709B2
- Authority
- JP
- Japan
- Prior art keywords
- density
- foam material
- foaming
- equation
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 238000005187 foaming Methods 0.000 title claims description 76
- 239000000463 material Substances 0.000 title claims description 39
- 238000004458 analytical method Methods 0.000 title claims description 38
- 238000009413 insulation Methods 0.000 title claims description 21
- 238000011156 evaluation Methods 0.000 title description 2
- 239000006261 foam material Substances 0.000 claims description 53
- 238000004364 calculation method Methods 0.000 claims description 30
- 239000006260 foam Substances 0.000 claims description 28
- 238000000034 method Methods 0.000 claims description 28
- 238000002347 injection Methods 0.000 claims description 21
- 239000007924 injection Substances 0.000 claims description 21
- 238000004134 energy conservation Methods 0.000 claims description 11
- 239000007787 solid Substances 0.000 claims description 11
- RGSFGYAAUTVSQA-UHFFFAOYSA-N Cyclopentane Chemical compound C1CCCC1 RGSFGYAAUTVSQA-UHFFFAOYSA-N 0.000 claims description 8
- 239000012948 isocyanate Substances 0.000 claims description 8
- 150000002513 isocyanates Chemical class 0.000 claims description 8
- 229920005862 polyol Polymers 0.000 claims description 8
- 150000003077 polyols Chemical class 0.000 claims description 8
- 239000004088 foaming agent Substances 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- DMEGYFMYUHOHGS-UHFFFAOYSA-N heptamethylene Natural products C1CCCCCC1 DMEGYFMYUHOHGS-UHFFFAOYSA-N 0.000 claims description 4
- 238000013461 design Methods 0.000 description 14
- JOYRKODLDBILNP-UHFFFAOYSA-N Ethyl urethane Chemical compound CCOC(N)=O JOYRKODLDBILNP-UHFFFAOYSA-N 0.000 description 11
- 238000005206 flow analysis Methods 0.000 description 11
- 239000007788 liquid Substances 0.000 description 11
- 238000010097 foam moulding Methods 0.000 description 10
- 239000012530 fluid Substances 0.000 description 8
- 238000000465 moulding Methods 0.000 description 8
- 238000009826 distribution Methods 0.000 description 7
- 238000002156 mixing Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 238000001125 extrusion Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 230000020169 heat generation Effects 0.000 description 3
- 239000011810 insulating material Substances 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 239000011347 resin Substances 0.000 description 3
- 229920005989 resin Polymers 0.000 description 3
- 238000001816 cooling Methods 0.000 description 2
- 238000007872 degassing Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012950 reanalysis Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000002076 thermal analysis method Methods 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 229920005830 Polyurethane Foam Polymers 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 238000009835 boiling Methods 0.000 description 1
- 230000005587 bubbling Effects 0.000 description 1
- 239000004035 construction material Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000003050 experimental design method Methods 0.000 description 1
- 239000012774 insulation material Substances 0.000 description 1
- 239000002984 plastic foam Substances 0.000 description 1
- 239000011496 polyurethane foam Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 229920001187 thermosetting polymer Polymers 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/22—Moulding
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
【発明の属する技術分野】
本発明は、プラスチック発泡成形加工技術に係り、冷蔵庫や建設材料、自動車等に使用される断熱材、衝撃吸収材等を発泡成形する際の3次元流動解析方法及びその方法を実行させるプログラムに関する。
【0002】
【従来の技術】
有限差分3次元流動解析または有限要素3次元流動解析方法に関して、密度変化のある流体の解析を行う場合には、圧力による圧縮性を考慮する手法が用いられている。この圧縮性を考慮した流体解析の例として、特開平7-334484号公報、特開平6-187321号公報があげられるが、密度が時間項を含む関数として増加する発泡挙動を伴う解析方法または解析プログラムの例は報告されていない。
【0003】
【発明が解決しようとする課題】
近年の3D-CAD・CAM・CAE環境にあっては、発泡流動に伴う圧力を考慮した部品配置の適正化など3D-CAEを用いたコンカレントな設計支援が、開発期間短縮、コスト低減のために必要不可欠である。しかし、冷蔵庫の断熱材として用いられているウレタン材料などの発泡挙動を解析する場合には、従来の圧力による圧縮性を考慮した圧縮性気体の解析方法では、そのウレタン材料等の発泡過程の物理的性質が、圧縮性気体の性質とは大きく異なり、正確に対応できない問題があった。また、冷蔵庫の断熱部分は、冷蔵庫の内装意匠や、冷却部や電源回路等の存在によって複雑な形状をしており、設計に際しては、断熱部分を充填するために必要なウレタン材料の最適量や、ウレタン材料の適切な注入個所を決めなければならない。しかし、従来の解析方法では発泡現象に対しては十分な解析ができず、製品の断熱部分の隅々にまで十分に発泡が行き渡らないといった製品不良が生じる場合もあった。冷蔵庫の場合では、発泡後の密度が大きく、発泡倍率が低いと冷却効率低下の原因となり、エネルギー消費効率の低下にもつながっていた。
【0004】
そこで、本発明は、ウレタン材料等の発泡現象における多種の発泡要因が発泡倍率に与える影響を明確にし、発泡倍率に最も影響を与える要因である治具温度および製品肉厚に対応した密度変化を考慮して解析する必要性を見出した。加えて、ウレタン発泡は熱反応を伴うので、発熱式をエネルギ保存式に代入し、熱反応率による密度変化を考慮する必要がある。
【0005】
よって、連続の式、運動方程式およびエネルギ保存式によって計算される3次元ソリッド要素を用いた流動解析法であって、密度が上記の時間および肉厚を含む関数、または熱反応の関数として入力されことを特徴とする3次元発泡流動解析法、その方法を実行させるプログラム及びそのプログラムをが格納された記録媒体を提供することを目的とする。そして、この発泡解析方法を用いて、製品の適切な設計を支援することを目的とする。
【0006】
【課題を解決するための手段】
上記課題を解決するため、本発明は、密度が時間および肉厚を含む関数、または熱反応の関数として入力されることを特徴とし、このタイムステップ毎の密度変化を連続の式、運動方程式およびエネルギ保存式に代入して計算することにより、圧力、流動速度、温度、密度分布などの結果を得ることができる3次元発泡流動解析法または発泡流動解析プログラムを提供する。なお、実際のウレタン材料などの発泡現象において、粘度は時間項を含む硬化反応を伴い、密度は時間変化に加え、粘度および圧力などによって変化する挙動を示すので、詳細な発泡流動解析を行うために、密度および粘度を、時間項を含む関数として入力し、タイムステップ毎の密度および粘度変化を連続の式、運動方程式およびエネルギ保存式に代入して計算することにより、圧力、流動速度、温度、密度分布などの結果を得る3次元発泡流動解析法または発泡流動解析プログラムを提供することもできる。また、冷蔵庫などの断熱材としての性能を評価するために、密度の関数として熱伝導率を求めることもできる。
【0007】
【発明の実施の形態】
以下、添付の図面を参照しながら、本発明に係る実施の一形態について説明する。
【0008】
まず、本実施の形態に係る発泡成形工程を図1に示す。これは、シクロペンタン(C5H10)発泡剤を充填したポリオール1とイソシアネート2の2液をミキシングヘッド3によって攪拌し、発泡型4の中に2液を攪拌した発泡材料5を充填することにより、発泡させる工程の例を示したものである。ここで、発泡に用いる2液をポリオール1とイソシアネート2としたが、本発明はこれだけに限定されるものではなく、発泡型4も任意の形状とできるものとする。
【0009】
この発泡成形において、治具温度などの成形条件、発泡材料が満たされる製品肉厚などの因子により、発泡挙動が大きく異なる。ここでは、各因子が、発泡材料の初期の体積から発泡成形後の体積の倍率で表す発泡倍率に与える影響を明確にするため、実験計画法を用いた検討を行った。発泡に影響を与える因子は、
(A)ガス抜き場所12の寸法
(B)発泡材料の充填場所
(C)治具温度
(D)液体温度
(E)ポリオール1とイソシアネート2の比率
(F)ミキシングヘッドの注入圧力
(G)発泡剤充填量
(H)製品肉厚
とし、因子(A)は2水準、(B)〜(H)は3水準の条件を変更することにより検討した。なお、因子(A)(B)(H)は冷蔵庫の形状設計および治具形状設計によって変更できる範囲で条件を変更し、因子(C)〜(G)は冷蔵庫の量産成形においてバラツキが予想される範囲内で条件を変更した。発泡に用いた治具形状13は図2に示す形状とする。結果を図3に示す。ここでは、発泡倍率を感度(dB)で表しており、感度が大きいほど、発泡倍率が大きいことを示している。このように、(C)治具温度および(H)製品肉厚の因子が発泡倍率に与える影響が支配的であることが分かる。
【0010】
このような治具温度と肉厚を考慮した発泡挙動を解析するためには、解析モデルを3次元ソリッド要素に分割し、連続の式(1)、運動方程式(2)〜(4)、エネルギ保存式(5)によって、流動速度、圧力、温度を求める手法を用いる。ここで、ρ;密度、u;x方向速度、υ;y方向速度、ω;z方向速度、T;温度、P;圧力、t;時間、η;粘度、Cp;定圧比熱、β;体積膨張係数、λ;熱伝導率を示している。
【0011】
これらの式をソフトウェア上で実現した場合には、「ρの時間変化を含む関数」、「粘度η」、「定圧比熱Cp」、「体積膨張係数β」は入力値として使用される。
【0012】
【数1】
【数2】
【数3】
【数4】
【数5】
ここで、密度変化のない液体を扱う場合には、∂ρ/∂t=0となる。しかし、冷蔵庫ウレタン材料等の発泡材料を扱う場合には、∂ρ/∂t≠0となり、前記した発泡倍率に影響を与える主因子である治具温度と肉厚を考慮した発泡材料の密度変化を別途求めることが必要となる。
【0013】
このとき、一つの治具の場所による温度バラツキおよび発泡成形毎の治具温度のバラツキを無視すると、可視化実験などの結果を用いて、ある治具温度における密度の変化を数式化することにより、密度が時間項および肉厚項を含む関数として表記できる。式(6)、(7)は、初期密度を1とし、5秒後から密度が時間の2乗に反比例し、肉厚に比例して発泡する材料の密度を示している。なお、H;製品肉厚を示している。
【0014】
【数6】
ρ=1 、(0≦t<5) (6)
【数7】
ρ(t)=H/((t−5)2+1)、(t≧5) (7)
この式(6)、(7)によって求められる各タイムステップにおける密度を、式(1)〜(5)に代入することにより、発泡時の流動速度、圧力、温度などの結果を、肉厚Hを考慮した上で求めることができる。この密度式(6)、(7)を用いると、材料が流動しながら発泡する過程を解析できる。
【0015】
ここで、図4に示す複数の6面体の要素により発泡充填部分が形成されている場合に、要素14における式(7)で示した肉厚Hは、要素14を含んで発泡充填部分の肉厚を構成するX,Y,Z方向の壁面(境界)間の距離Lx,Ly,Lzの最小値として表すことができる。また、図5に示すように、発泡型4や発泡充填部分が曲面形状を有する場合の肉厚Hは、解析対象個所の曲面形状の要素14を含む発泡充填部分の肉厚を構成するX,Y,Z方向の壁面(境界)間の距離Lx,Ly,Lzの最小値とすることで実現できる。図5の場合は、結果としてはLyが最小値(肉厚H)となるが、この算出方法としては、発泡充填部分のY方向でみた場合の二壁面の要素15および要素16を選定し、要素15の発泡充填部分を成し且つ他要素と接続していない境界の中央から、要素16の発泡充填部分を成し且つ他要素と接続していない境界への垂線17の長さ(Ly)として各要素の肉厚Hを求めることもできる。
【0016】
なお、式(7)において、密度は時間の2乗に反比例し、肉厚に比例する関数として表したが、本発明は、これだけに限定するものではなく、使用する発泡材料の特性に応じて、時間項および肉厚項を含む任意の関数(圧力、温度、粘度などの変数も考慮に入れた関数としても表すことができるものとする。)に対応した式で表すことができるものとする。
【0017】
例えば、式(7)においては、一つの治具の場所による温度バラツキおよび発泡成形毎の治具温度のバラツキを無視して密度の変化式を求めたが、本発明はこれだけに限定されるものではなく、密度を時間項、肉厚項および治具温度の関数とすることにより、熱計算で求めた治具の場所毎の温度に対応した発泡過程を計算することもできるものとする。
【0018】
一般的に、発泡倍率が高く、密度が小さい場合には熱伝導率が低くなり、発泡倍率が低く、密度が大きい場合には熱伝導率が高くなるので、上記で示した熱計算を行う場合には、発泡材料の熱伝導率を密度の関数として表すこともできるものとする。
【0019】
また、樹脂材料を発泡させながら押出し成形する場合には、押出し成形用ノズルと押出されてからの大気圧(圧力=0)間の圧力差が、材料を発泡させる主要因となり、樹脂粘度および表面張力が発泡形状を決める場合には、密度を圧力、樹脂粘度および表面張力を含む関数としても表すことができる。
【0020】
密度を、時間項と肉厚項を含む関数として入力する場合の発泡解析及び製品設計支援のフローチャートを図6に示す。ここでは、冷蔵庫の断熱部材として、ウレタン発泡材料を使用する場合を例に説明する。
【0021】
まず、モデル形状作成のステップ601では、冷蔵庫の断熱部分の初期設計モデルをCADデータ等として記憶装置から読み込む。次に、3次元ソリッド要素作成のステップ602では、読み込まれた初期設計モデルのCADデータに基づき、発泡材料が充填される断熱部分の形状を複数の特定形状(3次元ソリッドの有限要素)に分解する。次に、流体の物性値入力ステップ603では、解析を行う発泡材料の物性値である比熱、粘度式、密度式および熱伝導率式などを入力する。ここで、発泡材料毎に特有の密度の時間項および肉厚項を含む任意の関数を入力できるものとする。
【0022】
次に、境界条件、成形条件入力ステップ604において、3次元ソリッド要素内に流体が注入する際の初期速度、初期注入量、初期温度、注入場所および金型温度などの条件を入力する。ステップ605にて、ステップ603で入力した密度の時間変化式に従って、各計算過程に対応した密度の値を求め、この密度の値からステップ606にて熱伝導率を求めた後、ステップ605で求めた密度を用いてステップ607にて連続の式(1)および運動方程式(2)〜(4)により、流体速度および圧力などを計算する。更に、温度計算を行う場合には、各計算過程に於いて、連続の式(1)、運動方程式(2)〜(4)にエネルギ保存式(5)を連動させて計算させることもできる(ステップ608)。これらの式(1)〜(5)を用いた計算を収束判定(ステップ609)し、収束しない場合には、ステップ604の境界条件、成形条件およびステップ502の要素分割などを修正し、計算が収束する条件を選定する。計算が収束する場合には、冷蔵庫の断熱部分についての計算結果を記憶装置に記録し、この結果を計算機に表示させ、熱伝導率が部分的に高い部分から冷蔵庫内部の熱が外部に漏洩することなどを防止するために、ステップ610において、密度、熱伝導率などの適正判断を行う。このとき、密度、熱伝導率が部分的に高くなる場合には、冷蔵庫の断熱部分のモデル形状、または使用する発泡材料、または発泡材料の注入口、注入速度などの境界条件、成形条件を変更することにより、再計算を行う。なお、このステップ610における判定は、解析結果を人が判断し、条件などを変更して再解析を行なうこともできるし、密度の最大許容値を解析時に入力しておき、解析により求められた密度がこの許容値以上であれば、プログラムが自動的に注入口などを変更し、再解析することもできる。
【0023】
以上で示したステップにより、試作成形の前段階で、解析を用いた冷蔵庫の断熱部分形状および材料の注入口などの設計支援を行うことができる。
【0024】
また、以上では、密度を時間項と肉厚項の関数として表した例を示しており、肉厚をB,C,D(B<C<D)と変更した場合の比容積(密度の逆数)と時間の実験により求めた関係の一例を図7に示す。このように、肉厚が大きい場合には、発泡後の比容積が大きくなり、発泡が終了するまでの時間は短くなる、という発泡挙動の詳細が分かるが、肉厚条件毎に比容積の時間変化が異なり、密度の時間変化の定式化が困難である。よって、図3に示した発泡後の容積に影響を与える主因子である肉厚(B,C,D) および治具温度(F,G,H)を変更した場合の密度と時間の関係を定式化するため、無次元化した比容積と無次元化した時間条件の関係を図8に示す。
【0025】
ここで、ρ:発泡中の密度、ρi:発泡前の液体密度、ρmin:発泡後の最終密度、t:時間、t0:発泡しなくなる時間とした場合の無次元化した比容積を式(8)に、無次元化した時間を式(9)に示す。
【0026】
【数8】
無次元化比容積=((1/ρ)-(1/ρi))/((1/ρmin) -(1/ρi)) (8)
【数9】
無次元化時間=t/t0 (9)
このように、無次元化した比容積と無次元化した時間条件の関係は、一本の線で表されることから、発泡材料の条件によらないマスターカーブが求められる。また、図7で示したように、発泡後の密度および発泡しなくなる時間は、肉厚などの関数として表すことができ、式(8)で示した無次元化比容積は、式(9)で示した無次元化時間の関数として求められるので、密度は時間項を含む関数として表すことができる。
【0027】
このように無次元化した比容積と無次元化した時間条件から求められた密度の時間項を含む関数を、図6に示す603に示す密度式として入力することもできる。なお、発泡後の密度は、肉厚、治具温度、発熱量などの任意の関数として、発泡しなくなる時間は、肉厚、治具温度、粘度などの任意の関数として表すことができるものとする。
【0028】
また、図1で示したポリオール1とイソシアネート2から成るウレタン発泡材料は、粘度が時間の経過により高くなる熱硬化性の特性を示すので、密度の時間変化とともに、粘度の時間項を含む関数式を入力して計算することもできる。
【0029】
ここで、治具温度および肉厚が一定の条件で発泡成形した後の初期充填場所からの発泡材料の流動長さにおける密度を図9に示す。このように、発泡材料の流動は初期充填場所からの流動長さが長くなるに従い、密度が大きくなる性質がある。式(10)に初期充填場所からの流動長さと密度の関係を示す。ここでは、図9に示すように直線近似しており、Lは初期充填場からの流動長さを、Aは図9に示す近似直線の傾き(定数)を、Bは初期充填場所における密度(定数)を、ρ1は発泡成形終了後の密度を示している。
【0030】
【数10】
ρ1=AL+B (10)
式(7)で示した密度の関数は、可視化実験などにより、肉厚が等しければ、場所によらず密度は一定であると仮定している。しかし、実際には図9に示すように、場所により密度差があり、また発泡中の場所毎の密度差の測定は、大変困難である。よって、式(10)で示す発泡後の流動長と密度の関係から式(7)で示す発泡流動中の密度ρ(t)を補正する手法を用いる。
【0031】
冷蔵庫発泡成形において、注入場所を各種変更した検討を行うので、流動長で密度を補正するのは困難である。よって、各要素への発泡材料が流入する時間を用いて式(7)に示す密度を補正する手法を用いる。
【0032】
ここで、解析に用いる形状の肉厚;H1,幅;L1,注入量;Gとすると、発泡後の密度ρ1と時間tの関係は式(11)で示される。
【0033】
【数11】
ρ1= (G*A/(H1*L1))*(1/(H/((t-5)2+1)) +B (11)
よって、ある要素Bに時間t1において流入した発泡材料の要素Bにおける発泡後の密度は、式(12)で表される。なお、t1:ある要素への発泡材料の流入時間とする。
【0034】
【数12】
ρ1(t1)= (G*A/(H1*L1))*(1/(H/((t1−5)2+1))+B (t≧5)(12)
また、式(6)(7)で示した密度式は、場所による密度差を考慮していないので、全体の平均値を示す値である。ここで、各タイムステップにおける密度と流動長の関係は、発泡後の関係と等しくなると考えることにより、式(13)に、式(6)(7)で示した密度の時間変化式を式(12)で示したρ1(t1)により補正した密度の式を示す。なお、ρ3:発泡後の平均密度とする。
【0035】
【数13】
ρ(t)=H/((t-5)2+1)+ρ1(t1)-ρ3 (t≧5) (13)
この式(13)で示す密度式を図6で示したフローチャートのステップ603において入力することもできるものとする。
【0036】
以上では、発泡後の流動長さと密度の関係から、ソリッド3次元要素への流入時間を用いて密度変化式を補正した式を示すが、本発明はこれだけに限定されるものではなく、初期充填部分からの流動長Lなどを用いて、密度を補正することもできるものとする。
【0037】
冷蔵庫に用いるポリウレタン発泡は、シクロペンタン(C5H10)などの発泡剤を充填したポリオール1とイソシアネート2の混合による熱反応により、液体の発泡剤が気化することによって生じる。よって、発熱式をエネルギ方程式に代入することにより、二液混合時の熱反応を表現できる。発熱式を代入したエネルギ方程式を式(14)、熱反応率の時間変化式を式(15)、反応率の時間変化式を式(16)、温度の係数となる関数式を式(17)、熱反応率の式を式(18)に示す。なお、α:反応率、Q:単位質量当りの発熱量、d Q/dt:単位質量当りの発熱速度、Q0:単位質量当りの総発熱量、K(T):温度の係数となる関数、N、Kα、Ea:材料固有の係数を示す。
【0038】
【数14】
【数15】
dQ/dt= Q0K(T)(1-α(t))N (15)
【数16】
dα(t)/dt =K(T) (1-α(t))N (16)
【数17】
K(T) = Kαexp(-Ea/T) (17)
【数18】
α(t) = Q/Q0 (18)
上記のように、発泡中の熱反応を考慮することにより、発泡密度は熱反応率の関数として示すことができる。ここで、シクロペンタン(C5H10)発泡剤の沸点を49℃とすると、発泡材料の温度が49℃以下では、発泡はおこらないと考えられる。式(19)(20)は発泡材料の温度が49℃以下では、密度は発泡前の液体密度と等しく、49℃以上では、密度は熱反応率の関数として表している。なお、ρi:発泡前の液体密度、αf:49℃に達した時間における反応率、ρmin:発泡後の最小密度(飽和最小密度)を示す。
【0039】
【数19】
ρ=ρi (19)
【数20】
ρ=1/((α-αf)(1/ρmin)+1/ρi) (T≧49℃) (20)
また、一般的に、発泡倍率が高く、密度が小さい場合には熱伝導率が低くなり、発泡倍率が低く、密度が大きい場合には熱伝導率が高くなるので、発泡材料の熱伝導率を密度の関数として表すことができる。式(21)は熱伝導率の密度関数を示す。なお、J、M;材料固有の係数を示す。
【0040】
【数21】
λ=J(ρ/ρi)M (21)
ここで、粘度を反応率と温度の関数として、式(22)(23)(24)に示す。なお、η0:初期温度、αg:ゲル化した時の反応率、a、b、f、g:材料固有の係数を示す。
【0041】
【数22】
η=η0 ((1+α/αg)/(1-α/αg))C (22)
【数23】
η0=aexp(b/T) (23)
【数24】
C=f/T‐g (24)
なお、以上で示した発熱式、熱反応式、密度式、粘度式はこれだけに限定されるものではなく、任意の式を用いることができる。例えば、粘度式として、時間と温度の式を用いることもできるものとする。
【0042】
式(14)〜(24)で示した熱反応を考慮した場合のフローチャートを図10に示す。ここでは、冷蔵庫の断熱部材として、ウレタン発泡材料を使用する場合を例に説明する。
まず、モデル形状作成のステップ1001では、冷蔵庫の断熱部分の初期設計モデルをCADデータ等として記憶装置から読み込む。次に、3次元ソリッド要素作成のステップ1002では、読み込まれた初期設計モデルのCADデータに基づき、発泡材料が充填される断熱部分の形状を複数の特定形状(3次元ソリッドの有限要素)に分解する。次に、流体の物性値入力ステップ1003では、解析を行う発泡材料の物性値である比熱、熱反応式、発熱式、密度式、熱伝導率式、粘度式などを入力する。ここで、熱反応式、発熱式、密度式、熱伝導率式、粘度式は式(12)〜(22)で示した式を入力できるものとする。
【0043】
次に、境界条件、成形条件入力ステップ1004において、3次元ソリッド要素内に流体を注入する際の初期速度、初期注入量、初期温度、注入場所および金型温度などの条件を入力する。ステップ1005にて、ステップ1003で入力した密度の各タイムステップ毎の値を求め、この密度値と1つ前のタイムステップにおいて、連続の式、運動方程式および発熱式を入力したエネルギ保存式によって計算した圧力、温度および反応率を用いて、時間変化する密度を計算し、この密度値を用いてステップ1006にて熱伝導率を求める。
【0044】
ステップ1003で入力した密度式に従って、各計算過程に対応した密度の値を求めた後、この密度を用いてステップ1007にて連続の式(1)および運動方程式(2)〜(4)により、流体速度および圧力などを計算する。更に、ステップ1008にて、各計算過程に於いて、連続の式(1)、運動方程式(2)〜(4)に発熱式を代入したエネルギ保存式(10)を連動させて計算する。このとき、熱伝導率はステップ1006で求めた値を用いる。これらの式を用いた計算を収束判定(ステップ1009)し、収束しない場合には、ステップ1004の境界条件、成形条件およびステップ1002の要素分割などを修正し、計算が収束する条件を選定する。計算が収束する場合には、冷蔵庫の断熱部分についての計算結果を記憶装置に記録し、この結果を計算機に表示させ、熱伝導率が部分的に高い部分から冷蔵庫内部の熱が外部に漏洩することなどを防止するために、ステップ1010において、密度、熱伝導率などの適正判断を行う。このとき、密度、熱伝導率が部分的に高くなる場合には、冷蔵庫の断熱部分のモデル形状、または使用する発泡材料、または発泡材料の注入口、注入速度などの境界条件、成形条件を変更することにより、再計算を行う。なお、このステップ1010における判定は、解析結果を人が判断し、条件などを変更して再解析を行なうこともできるし、密度の最大許容値を解析時に入力しておき、解析により求められた密度がこの許容値以上であれば、プログラムが自動的に注入口などを変更し、再解析することもできる。
【0045】
以上で示したステップにより、試作成形の前段階で、解析を用いた冷蔵庫の断熱部分形状および材料の注入口などの設計支援を行うことができる。
【0046】
なお、式(1)〜(5)および(14)で示した連続の式、運動方程式およびエネルギ保存式は直行座標表示したが、円筒座標および球面座標など任意の座標表示できるものとする。また、式(1)〜(5)および(14)で示した連続の式、運動方程式およびエネルギ保存式において、重力項など任意に省略することができるものとする。
【0047】
ここで、解析を行うハードウェアの構成図を図11に示す。ここでは、計算装置6で作成したCADデータをLAN8を介して計算装置7に転送して計算する例を示している。計算装置7に転送されたCADデータは、計算装置7の記録装置10(ハードディスク、MOなど)に記録され、図2〜4で示したフローチャートに従い、計算を実行し、結果を記録装置10に記録した後、表示装置9に結果を表示できる構成となっている。
【0048】
発泡解析の一例を図12に示す。これは、図2に示した発泡治具形状13を用いて、前記した密度式(6)(7)を、連続の式(1)、運動方程式(2)〜(4)と連動させることにより発泡流動解析を行った結果である。注入場所11から発泡材料5を注入し、動粘度を15(cm2/s)、初期充填量300(cm3)、重力はZ軸のマイナス方向に加えた。ここでは、時間変化に伴う材料の発泡流動過程を示しており、0.2(s)、20(s)、30(s)、40(s)における発泡材料5の充填量を示している。このように、重力により、Z方向のマイナス方向の治具13に注入された発泡材料5が、Z方向のプラス方向に発泡流動している過程が解析できている。
【0049】
また、この解析による出力結果の一例として、密度分布を図13に示す。ここでは、前記した密度式(6)(13)を用いて、流動長による密度変化を求めており、図中の線が、密度の0.14〜0.16(g/cm3) (0.05(g/cm3)間隔)の値を示しており、密度が小さい場所の熱伝導率が低く、密度が大きい場所の熱伝導率が高いことが、式(21)に示した密度と熱伝導率の関係式から求めることができる。また、ここで示した密度の解析結果により、発泡材5のやわらかさを評価するため、弾性率を密度の式として入力し、求められた弾性率から発泡材5のやわらかさを求めることもできる。つまり、前記した図13において、密度が小さい場所は、弾性率が低く、やわらかく、密度が大きい場所は、弾性率が高く、かたくなることが、密度と弾性率の関数から求めることができる。
【0050】
また、以上で示した解析により熱伝導率が求められた発泡材5の断熱性をシミュレーションにより評価するための一例を図14に示す。ここでは、6面体を構成する断熱壁18〜23が、熱伝導率分布が計算された発泡材5で仕切られた境界条件を用いている。これは、XY平面を成す断熱壁18、19、YZ平面を成す断熱壁20、21、XZ平面を成す断熱壁22および発泡材5により囲まれた空間24の温度と、XY平面を成す断熱壁18、19、YZ平面を成す断熱壁20、21、XZ平面を成す断熱壁23および発泡材5により囲まれた空間25の温度に差を付け、空間24および25の温度変化を計算することにより、発泡材5の断熱性を計算できる事例を示している。このA-A断面のシミュレーション解析結果の一例を図15に示す。ここでは、空間24を冷蔵庫内の温度が低い領域、空間25を外気温とした場合の一定時間後の温度分布を示している。ここで、外気温をT2(一定値)と考えると、一定時間後には空間24の温度がT1となり、この温度T1は断熱材5の熱伝導率によって求められる値である。このように、発泡解析で求められた熱伝導率の値を用いて、熱解析を行うことにより、発泡材の断熱性能が計算できる。
【0051】
以上では、単純化したパネルの例を示したが、本解析を冷蔵庫の断熱層に用いることにより、複数個設置する材料注入場所の適正化、材料注入量の適正化、使用する発泡材料の適正化、断熱層形状の適正化などを試作成形をしないで行えるので、有効な設計支援ツールとして活用できる。また、断熱性を評価する場合にも、発泡材で囲まれた空間と外部との温度の計算により、断熱性を評価することもできる。
【0052】
以上、密度の時間変化を伴う計算方法を示したが、本発明は、密度の時間変化を入力して、計算する3D発泡流動解析プログラムおよびプログラムの記録媒体にも適用できるものとし、本プログラムを用いて、発泡材料充填量の適正化、発泡流動に伴う圧力を考慮した部品配置の適正化などを行った製品の製造方法にも適用することができる。なお、以上では、連続の式、運動方程式、エネルギ保存式を用いた解析フローチャートを示したが、本発明はこれだけに限定されるものではなく、連続の式、運動方程式を用いて、速度および圧力を求める解析を行うこともできる。
【0053】
【発明の効果】
本発明の実施形態により、発泡充填部の形状を考慮した正確な発泡解析を実現することができ、発泡材の弾性評価や断熱性評価による製品設計の支援を行うことができる。
【図面の簡単な説明】
【図1】 ポリオールとイソシアネートの2液混合による発泡過程
【図2】 発泡倍率に各因子が与える影響を明確にするための発泡実験治具
【図3】 各因子が発泡倍率に与える影響
【図4】 各要素における肉厚の算出方法(1)
【図5】 各要素における肉厚の算出方法(2)
【図6】 密度の時間項および肉厚項を含む関数を入力する場合の流動解析フローチャート
【図7】 肉厚条件を変更した場合の比容積(密度の逆数)と時間の関係
【図8】 肉厚および治具温度条件を変更した場合の無次元化した比容積と無次元化した時間の関係
【図9】 初期充填場所からの流動長さと密度の関係
【図10】 密度式が発熱項を含む関数として入力する場合の流動解析フローチャート
【図11】 発泡解析を行うハードウェア構成図
【図12】 発泡解析結果(時間変化に伴う材料の発泡流動過程)
【図13】 発泡解析結果(発泡材料の密度分布)
【図14】 発泡材の断熱性評価解析方法
【図15】 発泡解析で求めた密度分布を用いた熱解析結果
【符号の説明】
1…ポリオール 2…イソシアネート 3…ミキシングヘッド 4…発泡型
5…発泡材料 6…計算装置 7…計算装置 8…LAN
9…表示装置 10…記録装置 11…注入場所 12…ガス抜き場所
13…発泡治具形状 14…要素A 15…要素B 16…要素C 17…要素B外側境界の中央から要素Cへの垂線 18…断熱壁(XY平面、Zマイナス方向)
19…断熱壁(XY平面、Zプラス方向) 20…断熱壁(YZ平面、Xプラス方向)
21…断熱壁(YZ平面、Xマイナス方向) 22…断熱壁(XZ平面、Yマイナス方向)
23…断熱壁(XZ平面、Yプラス方向) 24…発泡材で仕切られた空間(1)
25…発泡材で仕切られた空間(2)
Claims (5)
- 所定の空間形状データを記憶装置から計算装置に取り込み、この空間を複数の特定空間である3次元ソリッド有限要素空間に分割し、時間項及び有限要素空間の肉厚項とを含む発泡材料の密度式及びその密度の項を含む熱伝導率式を計算装置に入力し、更に、有限要素空間に注入される発泡材料の少なくとも初期速度、初期注入量、注入場所の情報を入力し、前記密度式から発泡材料の注入開始からの時間経過に対応した発泡材料の密度を算出し、この密度を用いて前記熱伝導率式から発泡材料の熱伝導率を算出し、得られた密度及び熱伝導率を連続の式及び運動方程式に代入して発泡材料の少なくとも流動速度と圧力を計算することを特徴とする発泡材料の3次元発泡解析方法。
- 所定の空間形状データを記憶装置から計算装置に取り込み、この空間を複数の特定空間である3次元ソリッド有限要素空間に分割し、有限要素空間に注入される発泡材料の少なくとも比熱及び熱反応率を含む熱反応式、発熱式、反応率項を含む密度式、その密度の項を含む熱伝導率式及び粘度式を計算装置に入力し、更に、有限要素空間に注入される発泡材料の少なくとも初期速度、初期注入量、注入場所の情報を入力し、前記熱反応式及び発熱式から求められる発泡材の反応率を用いて発泡材料の密度を算出し、この密度を用いて前記熱伝導率式から発泡材料の熱伝導率を算出し、得られた密度及び熱伝導率を連続の式、運動方程式及びエネルギ保存式に代入して発泡材料の少なくとも流動速度、圧力及び温度を計算することを特徴とする発泡材料の3次元発泡解析方法。
- 前記発泡材料はシクロペンタン発泡剤を含むポリオール及びイソシアネート混合材料であることを特徴とする請求項1または2に記載の発泡材料の3次元発泡解析方法。
- 請求項1または2に記載の発泡材料の3次元発泡解析方法により算出された発泡材料の熱伝導率を用いて、該発泡材料を挟んで隣接する2空間の温度変化を計算することを特徴とする発泡材料の断熱性評価方法。
- 請求項1または2に記載の発泡材料の3次元発泡解析方法における各手順を記録したことを特徴とする記録媒体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001282549A JP3994709B2 (ja) | 2001-09-18 | 2001-09-18 | 発泡材料の3次元発泡解析方法及びそれを記録した記録媒体並びにそれを用いた発泡材料の断熱性評価方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001282549A JP3994709B2 (ja) | 2001-09-18 | 2001-09-18 | 発泡材料の3次元発泡解析方法及びそれを記録した記録媒体並びにそれを用いた発泡材料の断熱性評価方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003091561A JP2003091561A (ja) | 2003-03-28 |
JP3994709B2 true JP3994709B2 (ja) | 2007-10-24 |
Family
ID=19106184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001282549A Expired - Lifetime JP3994709B2 (ja) | 2001-09-18 | 2001-09-18 | 発泡材料の3次元発泡解析方法及びそれを記録した記録媒体並びにそれを用いた発泡材料の断熱性評価方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3994709B2 (ja) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4419719B2 (ja) * | 2004-07-02 | 2010-02-24 | 株式会社日立製作所 | 発泡樹脂材料の3次元流動解析方法 |
JP2008033380A (ja) * | 2006-07-26 | 2008-02-14 | Hitachi Appliances Inc | 製品の断熱性能の解析方法、および製品の断熱性能の解析プログラム |
JP4807246B2 (ja) * | 2006-12-12 | 2011-11-02 | トヨタ自動車株式会社 | ガス溶解度予測方法並びに発泡性樹脂の流動解析方法及びプログラム |
JP5492046B2 (ja) * | 2010-10-14 | 2014-05-14 | 株式会社日立製作所 | 多孔質体内の樹脂のボイド成長解析方法 |
JP5739838B2 (ja) * | 2012-05-31 | 2015-06-24 | 株式会社日立製作所 | 樹脂加圧注型装置および加圧注入方法 |
CN110000293B (zh) * | 2019-03-07 | 2020-08-07 | 南京航空航天大学 | 一种考虑非均匀温度场的复杂曲面件热成形模具设计方法 |
-
2001
- 2001-09-18 JP JP2001282549A patent/JP3994709B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2003091561A (ja) | 2003-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20050114104A1 (en) | Method and apparatus for modeling injection of a fluid in a mold cavity | |
KR100704809B1 (ko) | 3차원 발포 해석 방법, 그것을 이용한 제품 설계 지원 방법및 그것들을 기록한 기록 매체 | |
Alfaro et al. | Meshless methods with application to metal forming | |
US8095348B2 (en) | Support system, support method and support program of resin molded article | |
Ilinca et al. | Three‐dimensional finite element solution of gas‐assisted injection moulding | |
JP3994709B2 (ja) | 発泡材料の3次元発泡解析方法及びそれを記録した記録媒体並びにそれを用いた発泡材料の断熱性評価方法 | |
Tie et al. | Three-dimensional finite element method for the filling simulation of injection molding | |
Tao et al. | A non-iterative direct-forcing immersed boundary method for thermal discrete unified gas kinetic scheme with Dirichlet boundary conditions | |
Zhou et al. | Numerical filling simulation of injection molding based on 3D finite element model | |
JP3804400B2 (ja) | 3次元発泡解析方法、それを用いた製品設計支援方法及びそれらを記録した記録媒体 | |
JP2008033380A (ja) | 製品の断熱性能の解析方法、および製品の断熱性能の解析プログラム | |
Marin et al. | Influences of the mesh in the CAE simulation for plastic injection molding | |
JP2011022660A (ja) | 流体の数値計算装置 | |
Shen et al. | The thermochemical non-equilibrium scale effects of the high enthalpy nozzle | |
US7010474B1 (en) | Performance predictor for fabricated parts | |
Rayner | Multigrid flow solutions in complex two‐dimensional geometries | |
Mukutmoni et al. | Role of Accurate numerical simulation of brake cooldown in brake design process | |
Li et al. | Three‐dimensional simulations of non‐isothermal flow for gas penetration in complex cavity during gas assisted injection moulding process | |
Isaev et al. | Testing of numerical methods, convective schemes, algorithms for approximation of flows, and grid structures by the example of a supersonic flow in a step-shaped channel with the use of the CFX and fluent packages | |
US11651124B2 (en) | Anti-warping design method for resin molded article, recording medium, and anti-warping design device for resin molded article | |
Li et al. | Adaptive coupled arbitrary Lagrangian–Eulerian finite element and meshfree method for injection molding process | |
Polynkin et al. | 3d simulation of gas assisted injection molding analysis of primary and secondary gas penetration and comparison with experimental results | |
JP3023969B2 (ja) | 冷熱サイクル構造体温度の解析方法および金型装置系の設計装置 | |
Wang et al. | FEM of Gas‐Assisted Injection Molding Based on 3D Model | |
Kim et al. | Numerical simulation of PU foaming flow in a refrigerator cabinet |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041012 |
|
RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20060419 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070508 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070619 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070710 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070723 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100810 Year of fee payment: 3 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 3994709 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100810 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100810 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110810 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120810 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130810 Year of fee payment: 6 |
|
EXPY | Cancellation because of completion of term |