Flat Simulator Multi (Ver.2.0.0)
改良成果資料

2021/6
株式会社HASL
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

はじめにFlat Simulator Multi (Ver.2.0.0) では，多層流動に対する解析機能の計算安定性と用途拡大を目的として，2.5D FEM 定式化の改良/拡張に取り組みました．従来の多層流動に対する主な改良点は，以下に示す通りです．
1) 多層流体を対象とした発達流れの定式化と完全発達流れの定式化を総合する
ことで，計算安定性及び収束性の向上させるとともに計算負荷を低減させました．
2) 発達流れの定式化に，従来，解析対象外とされていた多層流体の界面包み込み現象を疑似的に表現可能な計算モデルを実装しました．
以下に各改良成果とその運用方法について解説します．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

11. 発達多層流れ1),2)と完全発達多層流れ3)の統合化解法
発達流れとは，流動方向に対して変化する流動状態，完全発達流れとは，変化しない流動状態をそれぞれ意味します．多層流動において界面形状が空間的に変化する状態は発達流れとして取り扱うのが妥当です．一方，従来の2.5D FEMは，完全発達流れを前提とした定式化を採用しています．Flat Simulator Multi (Ver.1.0.0)では，多層界面の空間的な変化を定量化するために2.5D FEMの定式化を改良し，発達流れの解析へ適用範囲を拡張しています．発達流れと完全発達流れの定式化の差は，表1に示す様に，多層流体が形成する界面の法線方向に対する応力の釣り合い式の表現の差に表されます．これらの定式化の詳細については，文献1)を参照して下さい．

表1 多層流体の界面上境界条件

1) 谷藤眞一郎，依藤大輔，鬼防崇，辰巳昌典：成形加工，‘定常非ニュートン粘性多層流動に対する一般化2.5次元有限要素法による定式化
第1報：定式化と理論検証’ 33(2),60(2021)
2) 株式会社HASL，谷藤眞一郎，多層流体解析プログラム，多層流体解析システム，及び多層流体解析方法，特許第6547157号，2019-07-24
3) Puissant, S., Demay, Y., Vergnes, B. and Agassant, J.F. : Polym. Eng. Sci.,34(3),201(1994)

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

2Flat Simulator Multi (Ver.1.0.0)には，発達流れの2.5D FEM 定式化に基づいて，多層流体の全流動領域を解析可能とする解析機能を実装しました．しかしながら，実用性を向上させるために，定式化を更に改良し，計算安定性や収束性を改善させる課題が残されていました．
一般的に，多層流動は，上流側で発達流れ状態であっても下流側で完全発達状態に漸近する傾向を示します．Flat
Simulator Multi (Ver.2.0.0)では，多層流動のこのような特徴に着目し，全流動領域に対する解析モデルを発達流れの定式化を適用する上流側モデルと完全発達流れの定式化を適用する下流側モデルに分離して解析する機能を実装しています．すなわち，完全発達流れの2.5D FEM 定式化に基づく解析機能，及び上流側モデルと下流側モデルを連結するインターフェイス機能が新規に実装されました．モデル分割に伴い，計算規模が縮小します．また，発達流れの定式化を適用する解析領域を限定することで計算収束性の向上が期待できます．
発達流れと完全発達流れの定式化では，多層流体が形成する界面の計算方法が異なります．発達流れの定式化では，前述したように，界面上法線方向の応力の釣り合い式：

 pl  2l

1 Dhl
1 Dhl 1  pl 1  2l 1
hl Dt
hl 1 Dt

(1)を利用して界面を計算します．
一方，完全発達流れでは，多層流体の圧力は共通と見なされ，層 l の流束 i 成分 qli は，粘度や流路形状を利用して計算される流動コンダクタンス Si を利用して次式で表されます．

qli   Sl

p
xi

(2)

圧力 p は，層毎に肉厚方向に平均化した連続方程式に(2)式を代入して得られる方程式を利用して求めます．
すなわち，
3
qli
  p 


 Sl
0
xi 
i 1 xi
i 1 xi 
3

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

3(3)多層流動において，最も流量の多い層 m を基準層とします．(3)式は，基準層 m を利用して以下に示す様に書き直されます．

  p  3   Slp  3p  qlm
S

S

S
0
 l
 

  m

m
xi  i 1 xi  Sm
xi  i 1
xi xi
i 1 xi 
3(4)

ここで，χqlm は，層 l と m の流量比であり，

 qlm 

Sl
Sm

(5)

と定義します．
(4)式は，基準層 m の平均流速ベクトルと流量比の等値面の法線ベクトルが，互いに直交することを意味します．
すなわち，基準層の平均流速ベクトルから求められる流線に沿って，流量比は一定となります．
流動コンダクタンスは，層厚に依存する関数であり，(4)式から求められる流量比 ˆ qlmを満足するように以下に示す
ニュートンラフソン法を利用して層厚 hl を計算します．
k 1
l

h h
k
l



ˆ qlm   qlm (hlk 1 ) qlm / hh  hlk 1

ここで， 上付き添え字 k は反復計算回数を表します．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

4

for l  1 ~ n, (l  m)

(6)

2. 疑似界面包み込みモデル4),5)
界面包み込み現象は，図1に示す様に低粘性流体が高粘性流体を包み込む現象として知られています．この現象は，矩形断面形状の流路を上流側に配置するフィードブロックタイプ多層コートハンガーダイ内の多層流動状態を予測する上で重視されています．

図1 界面包み込み現象と座標系の定義(η1/η2<1)4) 谷藤眞一郎，依藤大輔，鬼防崇，辰巳昌典：成形加工，‘定常非ニュートン粘性多層流動に対する一般化2.5次元有限要素法による定式化
第2報：疑似界面包み込みモデルの実装と実験検証’，プラスチック成形加工学会誌投稿済み（査読終了アクセプト，2021/12月号掲載予定）
5) 株式会社HASL，谷藤眞一郎，多層流体解析プログラム，多層流体解析システム，及び多層流体解析方法，特許出願済(2021/2/10)（審査中）

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved5粘性流動解析の範疇では，この界面包み込み現象を正確に表現する難しく，その定量化に際して粘弾性解析が必要とされています．また，この現象は，流動方向に対して界面形状が変化する発達流れの状態で発現します．一般的に，粘弾性解析は粘性解析と比較して計算負荷が増加し，解析の難度も高くなります．また，現状では，3D粘弾性流体解析法が進歩していますが，その適用範囲は基礎的な解析モデルに限定され，実用の多層押出ダイを対象とした解析事例は，数少ない状況です．運用が難しく，計算時間が長くなる粘弾性流体解析は，原理的に可能であっても工学支援ツールとして手軽に活用することは難しいと考えます．Flat Simulator Multi (Ver.2.0.0)では，界面包み込み現象
を解析対象とし，且つ計算負荷の増加を避け，運用が容易な解析技術を構築することを目的として，疑似的に当現象を表現可能な計算モデルを開発しました．以下に疑似界面包み込みモデルの理論背景について解説します．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved6粘弾性効果を近似的に考慮するために，解析的な取り扱いが容易であるCEF(Criminale Ericksen Filbey)モデル6)を採用します．CEFモデルでは，速度勾配テンソルLや変形速度テンソルDを利用した単純な表式で粘弾性余剰応
力テンソルτが次式に示す様に表されます．
τ  2 D  1 D  4 2 D  D(7)

ここで，変形速度テンソルD の upper convective time derivative は，次式で与えられます．



D

D
 u  D  L  D  D  LT
t(8)

図1内に示すように座標系を定義し，矩形断面に垂直なx方向(Machine Direction：MD)の流速成分uが，MDに一様になる一定粘度ηの定常流動を想定します．y方向(肉厚方向)の流速成分vやz方向（TD：Transverse Direction）の流速成分wの絶対値は，uの絶対値と比較して十分小さいとして無視します．

6) Criminale, Jr. W.O., Ericksen, J.L. and Filbey, Jr. G. L. : Arch. Rat. Mech. Anal., 1,410(1985)

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved7このような条件下では，(7)式より，以下に示すように余剰応力成分を求めることができます．

 xx   1  2  xy2  xz2
 xy  xy xz  xz
2
 yy   2xy    
2 xy xz
 yz
    2
2 xz
 zz

ここで，


(9)

u

 xy  y  u xz z

(10)としました．
2.5D FEMで採用する潤滑近似では，TDに垂直な側壁(図1内のz=±W/2の壁面)の影響が無視され，

xz  0 と見なされます．この場合，(9)式より，各種法線応力差係数の定義式
 N1   xx   yy   1xy2
2
 N 2   yy   zz   2xy(11)が導かれます．
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved8第2法線応力差は，第１法線応力差と比較して小さく5)，第2応力差係数を第1法線応力差係数を用いて実験的に，

2 0.15 ~ 0.10
1

(12)

と表現されます．
(9)式を利用して，運動方程式は，次式で表されます．

p  xx  xy  xz





 xy   xz 
x
x
y
z
y
z

(13)

p  xy  yy  yz






 2xy2   2xy xz 
y
x
y
z
y
z

(14)

p  xz  yz  zz





  2xyxz   2xz2
z
x
y
z
y
z

(15)









(13)式のMDの運動方程式は，粘性流体の運動方程式と一致します．一方，(14)式と(15)式の運動方程式に着目すると，第2法線応力差係数を含む速度勾配項が，矩形断面内の圧力勾配に寄与することが分かります．完全発達状態にある粘性流動では，MDのみに圧力勾配が発生し，断面内の圧力は一定です．第2法線応力差を考慮することで，MDの流速成分の速度勾配が，MDと垂直な断面内に圧力勾配を発生させる効果を示すのは，粘性流体と異なる特性と言えます．これらの圧力勾配は，近似的に無視した断面内の二次流れを発現させる要因になります．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

92.5D FEM では，層厚が薄いことを理由として各層内の肉厚方向に対する圧力勾配を無視する近似が採用されています．一方，TDに対する圧力勾配を表現する(15)式は，肉厚方向に対して平均化することで，2.5D FEMによる発達流れの定式化の拡張に利用可能です．第2法線応力差係数と流路肉厚が，TDに対して一定であることを前提とし，(15)式から，以下に示すように肉厚方向に平均化された圧力勾配の評価式が得らます．
H /2
p
1
p
 2   H /2 2 dy xz dy z
H  H/ 2 z
H z   H/ 2


(16)

定常単層粘性流動に対しては，(16)式を解析的に積分することが可能です．図1に示すような矩形断面内定常単層流動の(13)式の理論解は，次式で表されます．


 (2n  1) z 
cosh
2
2
n





1 p 1  2  H   4 H(1)H cos  (2n  1) y    y    3 
u  y, z  3 x  2 H
 2    n 1 (2n  1) cosh  (2n  1) W 






2H




(17)式よりひずみ速度は，次式で与えられます．


 (2n  1) z sinh nu
4 H 1 p(1)H
 cos  (2n  1) y  xz  2z  x  n 1 (2n  1) 2H
 (2n  1) W 


cosh




2
H





Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

10

(18)

(17)

(16)式に含まれるひずみ速度の2乗の肉厚平均値は，(18)式を利用し，



H /2

2
xz

1
2
2




dy


g f ( z)xzw
H  H/2(19)と表されます．ここで，

w 

W p
2 x

(20)

は，H→∞ とした際に側壁上(z=±W/2)で計算されるひずみ速度の絶対値です．また，形状因子 gf は次式で表されます．

 (2n  1) z  sinh 
2  
32  H 
1
H

 

g f ( z)  4   
4
(2
n

1)

W
  W  n 1 (2n  1)  cosh 
2H




2

(21)

(16)式をTD(z方向）に積分すると，
p ( z )  p (0)   2 xz2   2w2 g f ( z )(22)が得られ，第2法線応力差係数が負であることに留意すると，圧力の肉厚平均値は，TDに対して低下することが分
かります．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

11第2法線応力差係数を0として扱う従来の2.5D FEM では，(22)式に示す様に断面内の圧力は，中心部圧力の肉厚平均値に一致するとされています．第2法線応力差係数を考慮することで表現されるTDに対する圧力勾配を定式化に反映させるために，(22)式が多層流体の各層のTDに対する圧力勾配を表現すると仮定し，(1)式に示した界面法線方向の応力の釣り合い式を以下に示すように補正します．

 pl  2l

1 Dhl
1 Dhl 1
  2 l2   pl 1  2l 1
  2 l21
hl Dt
hl 1 Dt(23)

ここで，各層のひずみ速度の肉厚平均値は，次式で与えられます．

2



2
l

 W pl  g f (z) 2l x 

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

12

(24)

形状因子gf は，図2に示す様にアスペクト比：H/Wに依存して変化します．アスペクト比が1以下の条件下では，
H/W=1，すなわち正方断面流路の場合が最大値となり，第2法線応力差の影響が顕著になります．このように，従来無視されていた第2法線応力差係数と形状因子で表現される側壁効果を考慮することで，多層流体が形成する界面の包み込み現象を疑似的に表現することが可能になります．

図2 形状因子gf の TD座標(z/W) 及びアスペクト比(H/W)依存性

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

133. Flat Simulator Multi (Ver.2.0.0)の運用方法
図3 に示す様にFlat Simulator Multi (Ver.2.0.0)は，発達流れの定式化に基づく解析プログラム（DevelopingSolver.exe)と完全発達流れに基づく解析プログラム(FullyDevelopedSolver.exe)から構成されます．前述した疑似界面包み込みモデル(PEM : Pseudo Encapsulation Model)は，DevelopingSolverに実装されており，解析モデルに予め適用領域を規定して運用します．各解析モデルは，単独での運用あるいは連携した運用が可能です．図3に示す様に単独で運用する場合，多層合流部を含む上流側の領域の解析には，DevelopingSolverを適用します．一方，下流側の多層完全発達流れ領域の解析にはFullyDevelopedSolverを適用します．一方，多層押出ダイの全域を対象とする場合，図4に示す様に解析モ
デルを上流側と下流側に分離し，各領域にそれぞれDevelopingSolverとFullyDevelopedSolverを適用することが合理的です．この場合，図5に示す様に上流側の発達多層流動解析を先行して実施し，その結果を継続して実施する下流側の完全発達多層流動解析の流入口境界条件として考慮します．

発達多層流動解析プログラム：
DevelopingSolver.exe & PEM model

完全発達多層流動解析プログラム：
FullyDevelopedSolver.exe

発達多層流動解析
結果

完全発達多層流動解析結果

図3 Flat Simulator Multi (Ver.2.0.0) 構成解析プログラムの単独運用

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

14Inlet601) 単層流路領域:(DevelopingSolver適用)

合流断面

A’

A

2) 多層（フィードブロック）領域：(DevelopingSolver & PEM適用)
B’

モデル結合断面 B

303) コートハンガー領域(FullyDevelopedSolver適用)
C’

Outlet600

図4 Flat Simulator Multi (Ver.2.0.0) 構成解析プログラムの解析対象領域（モデルをB-B’断面で分離し，上流側にDevelopingSolver，下流側にFullyDevelopedSolverを適用）
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

15

1st step : 上流側発達多層流動領域の解析

2nd step : 下流側完全発達多層流動領域の解析

発達多層流動解析プログラム：
DevelopingSolver.exe & PEM model

完全発達多層流動解析プログラム：
FullyDevelopedSolver.exe

発達多層流動解析
結果：xxx.fluxoutlet

上流側発達多層流動解析
結果を反映した下流側完全発達多層流動解析結果

図5 Flat Simulator Multi (Ver.2.0.0) 構成解析プログラムの連携運用

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

163.1 発達多層流動解析プログラム(Developing Solver) の運用方法
発達多層流動解析プログラムの運用方法は，Ver.1.0.0システムと同様ですが，幾つかの改良を施しています．以下に関心が高いフィードブロックタイプコートハンガーダイを解析対象とする運用例について解説します．

1) フィードブロック流路モデル作成
発達多層流動解析プログラムは，上流側の単層流動領域及び合流直後のフィードブロック流路内の発達多層流動領域を解析対象とします．一方，下流側コートハンガーダイ内の多層流動は，後述する完全多層流動解析プログラムが解析対象とします．ここでは，発達多層流動解析プログラムが解析対象とする三層フィードブロック流路モデルの作成方法について解説します．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

17最初にテンプレートを利用して中間層(層番号2)のフィードブロック単層流路モデルを作成します．本生成例では，マニフ
ォールド断面形状を円型状としてモデル化していますが，他形状も自由に作成可能です．

図6 中間層(層2)のフィードブロック単層流路モデル

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

18最初にテンプレートを利用して中間層(層番号2)のフィードブロック単層流路モデルを作成します．
モデル保存ボタンを押してメッシュ情報を保存します．モデルは，流入口を起点（y=0)として，MD(-y)方向に作成されます．メッシュセーブフォームのMD最小座標は，流出位置，最大座標は，流入位置を意味します．下流側コートハンガーダ
イモデルが，y=0を起点として-y方向に作成されていることと合流後の多層流路の長さを考慮し，並進移動量を設定してメッシュ情報を保存します．本作成例では，y=-125~0で作成された元モデルをdy=250の並進移動した位置に保存します．
すなわち，保存後のモデルは，y=125~250の範囲に定義されます．

図7 メッシュ情報の並進移動保存（中間層）
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved19

表層(層番号1,2)の単層流路形状が中間層と等しい場合，紙面と垂直方向に並進移動した位置にメッシュ情報を保存します．本生成例では，下表層(層番号1)のメッシュ情報は，(dx,dy,dz) : (0,250,-80)，上表層(層番号2)のメッシュ情報は，(dx,dy,dz) : (0,250,80)の位置に並進移動して保存します．従来のVer.1.0.0 システムでは，各メッシュ情報を保存後，メッシ
ュ生成を再度行って，情報を保存していましたが，Ver.2.0.0 システムでは，流路形状が等しい場合，メッシュの再生成は不要であり，メッシュ保存ボタンを連続的に押すことで異なる並進位置に複数のメッシュ情報の保存が可能です．各情報を保存後，システムを一旦終了します．

図8 メッシュ情報の並進移動保存（下表層）Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

20システムを再起動し，中間層(層番号2)のメッシュ情報をインポートした後，修正/2.5Dメッシュ追加をプルダウン選択します．
表示されるフォームにおいてAppend import チェックボックスをチェック状態とし，下表層(層番号1)と上表層(層番号3)のメッシ
ュ情報を順にインポートします．

図9 メッシュ情報のAppend import

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved21フォーム内のMesh modify ボタンを押し，Rotation ラジオボタンをONとし，Box pick チェックボックスをチェック状態として，下図に示す様に下表層のランナーを除く計算要素をボックスピック選択します．次にNode focus ボタンを押し，選択領域の上左単の節点を一回目のマウスクリック選択，上右単の節点を二回目のマウスクリックで選択します．この操作により，選択要素の上端片が回転軸として選択されます．軸方向は，一回目クリック位置（左）から二回目クリック位置（右）方向に設定されています．
右ねじの法則に従って回転角度を指定します．回転軸に対して時計周りを正，半時計周りを負とします．本回転例では，下表
層は，反時計回り，すなわち-45°，上表層は，時計回り45°に回転させます．回転角度を設定以後，Modifyボタンを押すこと
で，選択領域の計算要素が指定回転軸周りに回転します．

図10 選択要素の回転移動
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

22

各表層の計算要素を下図に示す様に軸周りに回転した後，Mesh export ボタンを押してメッシュを保存します．作業を継続する際にシステムを再起動する必要があるため，メッシュ情報は，この段階において適当な名称を付け必ず保存して下さい．

図11 表層要素の回転移動例

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

23下流側に装着されるコートハンガーダイに対するメッシュ情報を予め作成しておき，当メッシュ情報をインポートします．
以下に示す操作では，上流側フィード単層流路と下流側コートハンガーダイ多層流動領域を連結するフィード多層流動領域を作成します．

図12 下流側コートハンガーダイメッシュ情報のAppend import

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

24フィード多層流動領域のコーナに下図に示す様にポイントを作成します．ポイントアイコンをマウスクリック選択し，NodeFocus Click チェックボックスをチェック状態として下図に示す様にモデルのコーナをマウスクリック選択し，8 ポイントを作成
します．作成後のポイント情報は，Geometry exportボタンを押して適当な名称を付け保存します．

ポイントアイコン

3
5
1

4
7

6

8

2

Node Focus
Click チェックボ
ックス

Geometry export
ボタン

図13 コーナポイントの作成
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved25次にラインアイコンをマウスクリック選択し，2ポイントを結ぶラインを作成します．Point 1,2の順にマウスクリック選択するとpoint 1，2を結ぶラインが作成されます．この機能を利用し，下図に示す様に3つの4角形領域を作成します．作成後のライン情報は，Geometry export ボタンを押して適当な名称を付けファイルに保存します．

図14 領域の輪郭形状を表現するラインの作成

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

26

システムを再起動し，上流側フィード単層流路のメッシュ情報を読み込み，修正/2.5Dメッシュ追加をプルダウン選択して表示されるフォームにおいて，Geometry import ボタンを押し，前述の操作で作成したライン情報を読み込みます．

図15 上流側フィード単層流路メッシュ情報とライン情報の読み込み
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved27

メッシュ生成アイコンをマウスクリック選択し，要素分割領域の横（幅）方向と縦（高さ）方向の分割数と肉厚を指定します．プ
ロパティ番号は，1以外を設定して下さい（1はマニフォールド部で既に利用しているため）．Line Focus Click チェックボックスをチェック状態とし，下図に示す様に四角形領域の構成辺を下右上左の順序でマウスクリック選択します．マウスクリック選択時に交換を促すメッセージが表示されたらOKボタンを押します．最後のラインを選択した後，正常にループが形成されたというメッセージが出力されるのでOKボタンを押して下さい．この場合は，Generateボタンを押すことで，正常にメッシュを生成可能です．ループがクローズしていないというメッセージの場合は，メッシュの生成は不可です．この場合は，ポイントやライン
情報が連結していないという意味ですので，これらの情報の確認と再生成が必要です．縦（高さ）方向の要素分割数は任意ですが，横（幅）方向の要素分割数は，上流側フィード単層流路と下流側コートハンガーダイのメッシュ情報と等しく設定する必要があります．

メッシュ生成アイコン

図16 四角形要素分割領域の選択
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

28

以上の操作を繰り返し，下図に示す様に連結流路領域の全てをメッシュ分割します．本生成例では，合流領域の肉厚は20mm，要素分割数は20×40，上流側単層流路の肉厚は1mm，分割数は20×10としています．分割後は，隣接共有辺上に重複して生成されている節点をマージするために，Node mergeボタンを押します．節点マージ後，Mesh export ボタンを押して適当な名称を付けメッシュ情報をファイルに保存します．

図17 要素分割と節点のマージ
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved29

次にMesh modify ボタンを押し，Layer configuration ラジオボタンをONとして層構成情報を設定します．本作成例では，Maximum layer number ： 最大層数は3と入力します．YZボタンを押してモデルを側面から観察し，Box pick チェックボックスをチェック状態とし，層構成を設定する要素を選択します．選択後，Layer configuration ：層構成を設定し，Modify ボタンを押すことで層構成情報が設定されます．本作成例では，下図に示す様に４つの領域に層構成を設定しています．層構成情報の設定後，Mesh export ボタンを押して適当な名称を付けメッシュ情報を保存します．多層解析用メッシュファイルの拡張子は，
multimshに自動的に変更されます．

領域①
層構成：1

領域②
層構成：2

領域③
層構成：3

領域④
層構成：1,2,3

図18 層構成情報の設定
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved30次に修正/2.5D メッシュ加工をプルダウン選択し，表示されるフォームでフレーム固定チェックボックスをチェック状態とし，ボックスピックラジオボタンをONとして要素全体をボックスピック選択します．温度境界条件変更ラジオボタンを押し，熱伝達係数と環境温度を設定して更新ボタンを押します．

図19 温度境界条件の設定
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved31必要に応じて要素肉厚情報を変更します．肉厚を変更する要素を選択します．肉厚ラジオボタンをONとし，区分線形補間ラジオ
ボタンをONとして表示されるフォームで，MD補間ラジオボタンをONとするとMD（高さ方向）の肉厚分布を区分線形補間で設定可能です．高さ比率と肉厚の関係を設定し，追加ボタンを押します．ここで高さ比率0は，選択要素の下端（最下流側），1は上端（最上流側）を意味します．更新ボタンを押すことで，下図に示す選択要素の肉厚のMD（高さ方向）の分布が更新されます．

図20 MD方向肉厚分布の変更
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

32

Ver.2.0.0に新しく実装された疑似界面包み込みモデル（PEM : Pseudo Encapsulation Model)を適用する領域には，プ
ロパティ番号を999に変更します．すなわち，プロパティ番号999の流路では，第2法線応力差がもたらす界面包み込み現象が疑似的に表現されます．PEMを適用する領域は，矩形断面を特徴とする流路に限定されます．合流直後の領域
は，流動状態の変化が激しく，流路アスペクト比(流路肉厚H/流路幅W)が小さため，PEMの適用領域から除外します．
ケースバイケースですが，合流位置より数メッシュ（本設定例では5メッシュ）より下流側の矩形断面流路を表現する計算要素にプロパティ番号999を設定します．PEMを適用する要素選択後，Property 番号ラジオボタンをONとし，プロパティ番号のテキストボックスに999を入力し，更新ボタンを押します．

図21 疑似界面包み込みモデル適用領域（プロパティ番号999)の設定
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved33ツール/Boundary condition set をプルダウン選択し，Qアイコンをマウスクリック選択します．3か所の流入口をマウスクリッ
ク選択して，決定ボタンを押し，流量にダミーパラメータ1，温度拘束条件同時設定パネルの設定チェックボックスをチェック状態として，設定ボタンを押します．

図22 流入口境界条件の設定

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved34

P アイコンをマウスクリック選択します．流出口をボックスピック選択して，決定ボタンを押し，圧力を0，温度拘束条件同時
設定パネルの設定チェックボックスを非チェック状態として，設定ボタンを押します．以上の操作で解析モデルの作成作業は終了です．ファイル/メッシュファイルエクスポートをプルダウン選択し，適当な名称を付けてメッシュ情報をファイルに保存します．発達多層流動解析では，当ファイルをメッシュファイルとして利用します．

図23 流出口境界条件の設定

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved35多層フィードブロックタイプコートハンガーダイの解析では，合流部断面形状（フィードブロックプロファイル）の変化が界面形成状況に及ぼす影響の分析が重視されます．2.5D 解析では，前述の手順で作成したFEMモデルを再生成し直す必要は無く，共用可能なメリットがあります．前述した操作で作成した解析モデルのメッシュファイルをインポートし，修正/2.5D メッ
シュ加工をプルダウン選択します．以下に示す操作では，合流位置から上流側3メッシュ分の流路肉厚をTD（幅方向）に対して指定した分布に変更する方法について解説します．各層のみを単独表示すると合流位置が不明になるため，下図に示す様にフレームチェックボックスをチェック状態とし，合流位置の目印となる要素をマウスクリック選択します．

図24 フィードブロックプロファイルの設定方法（合流位置確認用の要素選択）

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved36中間層(層番号2)の肉厚分布を設定するには，層表示番号を2として層表示変更ボタンを押します．当操作により表層(層番号1,3)の情報はマスクされます．先に選択した目印を基準に上流側3メッシュ分の要素をボックスピックで選択します．肉厚ラジオ
ボタンをONとし，区分線形補間ラジオボタンもONとします．表示される区分線形補間設定フォームにおいて，TD補間ラジオボ
タンをONとし，幅比率と肉厚の関係を追加ボタンを押してリストボックスに登録します．この登録情報は，保存/読込ボタンを押すことでファイルへの保存，ファイルからの読み込みが可能です．幅比率０はセンター，１はエッジを表し，更新ボタンを押すことでセンターに対して対称に肉厚分布が設定されます．本解析モデルでは，高粘性層である中間層の層厚が界面包み込み現象によって薄肉化するため，その不均一性を相殺するようにエッジ部のクリアランスを0.5mm拡大しています．

図25 フィードブロックプロファイルの設定方法（中間層フードブロックプロファイルの変更）
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

37

同様の操作に従って，表層(層番号1,3)の肉厚は，界面包み込み現象を相殺させるために中間層とは逆にエッジ部のクリアラ
ンスを下図に示す様に縮小させています．肉厚変更後は，元のファイル名とは別名でメッシュ情報を保存します．2.5D FEMではメッシュの再生成を行うこと無く，TDに対する肉厚分布情報を変更することでフィードブロックプロファイルの変化が多層界面形状に及ぼす影響を分析可能です．

図26 フィードブロックプロファイルの設定方法（表層フードブロックプロファイルの変更）
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved382) 多層押出解析条件の設定
発達多層流動解析プログラムを利用する場合，Solver selection パネルでDeveloping ラジオボタンをONとします．完全発達多層流動解析と比較して収束性が悪く，数百～数千回の非ニュートン反復計
算を要します．また，妥当な解析結果を得るために圧力差寄与係数や界面計算緩和係数を調整する必要があります．これらのパラメータは，1以下の正値に設定します．値か小さいほど，計算安定性が向上しますが，収束解を得るまでの反復計算回数が増加します．

多層解析用メッシュ
ダミーファイル名（存在しないファイル名を指定）
多層条件設定ボタン
値が大きいほど，圧力差が形状変化に及ぼす影響
大（最大値1:0.01程度が適正値，問題によっては更に小さく設定する必要あり）
多層界面計算の未緩和係数（最大値1:0.01程度が適正値，問題によっては更に小さく設定する必要あり）

非ニュートン反復計算回数は数百～
数千回の設定で十分

Developing ラジオボタンをon

図27 発達多層流動解析の基本的な条件設定(プロジェクト名:Upperfeedblock)
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved39多層押出解析で採用する押出量条件や材料物性ファイルは，多層条件設定ボタンを押して表示される下記フォームに設定します．物性データファイルと押出量，及び第2法線応力差係数を層数分設定し，Ver.1.0.0と同様，Evolution parameter を反復計算回数の増加に伴って1に漸近するように設定します．粘度比や流量比が大幅に異なるような非線形性の強い問題の場合は，当パラメータの増加を緩やかにするよう調整を要します．Evolution parameter の増加に伴って，多層流体間の相互作用が厳密に考慮されます．Evolution parameter を1とするのが理想的ですが，非線形性の強い問題の場合は，1に到達する前に解析不能になる場
合があります．そのような場合には，結果が得られる条件内で解析結果を評価する妥協が必要です．

物性データファイルを選択
設定（層数分）第2法線応力差係数を設定（層数分：負の値に設定）

押出量物性を設定（層数分）

Evolution パラメータをを設定

図28 発達多層流動解析の条件設定フォーム
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

40Ver.2.0.0 で新たに入力データとして追加された第2法線応力差に依存して界面形状が変化します．第2法線応力差の精密な測定は難しく，(12)式に示したように，大雑把に第1法線応力差の-10~-15%程度に設定されることが多いです．第1法
線応力差や第2法線応力差は，温度とひずみ速度に依存して変化します．一般的に第2法線応力差係数の絶対値は，ひずみ速度や温度の上昇に伴って低下する傾向を示すと言われています．粘度と同様，これらの非線形性を考慮することに技
術的な障害はありませんが，第2法線応力差の収集が難しいことと PEM が様々な近似に立脚していることを理由として，
Ver.2.0.0では，第2法線応力差係数を層毎に一定のパラメータ（負の値）に簡略化し，実測される層厚のフィッティングパラメ
ータとみなしています．200℃の温度条件下のLDPEに対して報告されている第1法線応力差は，ひずみ速度0.1~10 s-1の範囲において，103～104 Pa の範囲で変化します．仮に界面上のひずみ速度を10 s-1，第1法線応力差を104 Pa，第2法線応力差が第1法線応力差の-15%として，第2法線応力差係数を試算すると-15 Pa・s2となります．解析例では，中間層の界面上のひずみ速度が,およそ5 s-1であったため，第2法線応力差係数を -60 Pa・s2と試算しています．このように第2法線応力
差係数の絶対値は，1~102のオーダーに設定することを推奨します．一般的に，多層流動では，安定性を確保するために中間層の粘度を高めた層配置が採用されることが多く，低粘性層の表層が，中間層を包み込む現象が問題視されています．
粘度の高低と第2法線応力差係数の絶対値の高低は，対応関係にあると考えられますが，極端に高粘性層の第2法線応力差係数の絶対値が，低粘性層のよりも大きい場合には，高粘性層が低粘性層を逆に包み込み，実測とは逆の傾向が示されます．第2法線応力係数の絶対値が大きいほど，界面形状の変化が大きくなりますが，第2法線応力係数の比率も界面形成状態に強く影響します．第2法線応力差を適切に設定しないと，界面が押しつぶされ，解析不能に陥ります．妥当な解析結果を得るには，第2法線応力差を適切に調整するためのケーススタディを要しますが，計算時間は短く，1ケース当たり，数分～数十分の計算時間内に解析結果が得られます．
物性データファイルと第2法線応力差係数，及び押出量を層数分設定し，Evolution parameter を設定した後，確定ボタンを押します．閉じるボタンを押して当ホームを閉じ，解析実行タブメニューで保存ボタンを押して解析条件をファイル保存します．実行ボタンを押すことで発達多層流動解析プログラムが起動し，多層熱流動解析が行われます．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved413) 多層熱流動解析結果のポスト処理
コンター図の作画方法は従来通り，表示層番号を設定し，作画内容をコンボボックスで選択します．Ver.2.0.0では，各層の肉厚を3D表示する機能が新規実装されました．コンターコントロールパネル内の肉厚を反映チェックボックスをチェック状態として描画ボタンを押すと下図に示す様に各層の肉厚を反映した3D化表示が可能です．

肉厚を反映
チェックボックス

表層(層番号1)

中間層(層番号2)
図29 各層の肉厚を反映した3D化表示

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved423.2 完全発達多層流動解析プログラム(Fully Developed Solver) の運用方法
完全発達多層流動解析プログラムは，下流側コートハンガーダイ内の多層流動状況を効率的に定量化可能です．当解析プログラムは，単独運用の他，前述したように上流側発達多層流動解析結果を解析モデルの流入境界条件に反映した連携運用が可能です．以下に，当解析モデルの運用法について解説します．
完全発達多層流動解析プログラムを利用する場合，Solver Selection パネルにおいてFully developed ラジオ
ボタンをONとします．この際，Developing → Fullydeveloped パネルが表示されます．同パネル内のInterface チェックボックスをチェック状態とすると，連携ファイル(拡張子:fluxoutlet)の設定テキストボックスとSelectボタンが表示されます．この場合，連携ファイルを利用した連携運用となります．単独運用の場合，Interfaceチェックボッ
クスを非チェック状態とします．

図30 完全発達流動解析プログラムの選択と上流側発達多層流動解析結果連携ファイルの選択

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

433.2.1 完全発達多層流動解析プログラム(Fully Developed Solver) の単独運用
単独運用の場合，以下に示す手順に従って操作を進めます．

1) コートハンガーダイモデル作成
多層コートハンガーダイを作成する場合，Multi layer runner set チェックボックスをチェック状態とし，上流側フィートブロック矩形断面流路の半幅と肉厚をmm単位で指定します．
Multi layer runner set
チェックボックス

矩形断面流路半幅
矩形断面流路肉厚

図31 多層コートハンガーダイ基本形状設定フォーム
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

44従来の方法に従って，2.5D 解析用メッシュを作成し，メッシュ情報とデータを保存します．次に，修正/2.5Dメッシュ追加をプルダウン選択し，表示されるフォームでLayer configuration ラジオボタンをonとします．モデル全体をボックスピックで選択し，最大層数と層構成番号を設定し，Modifyボタンを押しで多層解析用メッシュに変更します．変更後，close ボタンを押し，当フォームを閉じます．

図32 多層コートハンガーダイ解析用メッシュへの変更

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved45

多層解析用メッシュには，部分的に温度境界条件が設定されていない領域があるため，修正/2.5Dメ
ッシュ加工をプルダウン選択して表示されるフォームで，温度境界条件を設定します．フレーム固定チェックボックスをチェックし，ボックスピックで全要素を選択します．温度境界条件変更ラジオボタンをONとし，熱伝達係数と環境温度を設定し，更新ボタンを押します．メッセージボックスのOKボタンを押した後，Close ボタンを押して当フォームを閉じます．

図33 多層コートハンガーダイ解析用メッシュの温度境界条件設定
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved46

ツール/Boundary Condition Set をプルダウン選択し，Qアイコンをマウスクリックします．ボックスピッ
クで流入境界に位置する節点を一括選択し，決定ボタンを押します．流量規定タブメニューでは，仮の流量1cm3/secを設定し，温度拘束チェックボックスをチェック状態として，境界温度を設定し，設定ボタ
ンを押します．以上の操作で解析上必要になる流入口境界条件が設定されます．流出口には，デフォルト状態で圧力0既定/断熱境界条件が設定されています．

図34 多層コートハンガーダイ解析用メッシュの流入口境界条件設定
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

47仮流量1を設定した節点には，解析条件として設定される押出流量が均一に配分されます．意図的に流量配分を指定することも可能です．流量配分を指定する場合は，流量を設定する節点を選択した
後，流量規定タブメニューに新規実装された分布設定ボタンを押します．

分布設定ボタン

図35 流量配分用分布設定ボタン
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

48表示された流量配分設定フォームにおいて，流量配分比の幅比率依存性を区分線形分布で指定します．ここで，幅比率は，0がセンター，1がエッジを表します．流量配分比は，センターに対して左右対称に設定されます．幅比率と流量配分比のテキストボックスに，設定値を入力し，追加ボタンを押すことで，設定値がリストボックスに登録されます．２点以上登録すると分布状態がグラフプロットされます．ポイントドラッグ移動チェックボックスをチェック状態とす
るとグラフプロットされているポイントをマウスドラッグで移動可能になります．登録情報は，読込/保存ボタンを押すことで読込/保存が可能です．設定層に流量配分を設定する層番号を入力し，設定ボタンを押すことで流量配分比が設定されます．流量配分比は，全層に対して設定する必要があります．下記の例では，層数：3とし，表層：1,3には，エッジ部がセンター部より流量配分比が高い設定，中間層：2には，均一な流量配分比を設定しています．
表層:1,3の不均一流量配分比設定

中間層:2の均一流量配分比設定

設定ボタン

流量配分比を設定する
層番号を指定図36 各層毎の流量配分比の設定
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved49層 l の節点 i に配分される流量 Qi lは．次式で表されます．

w
Qil  n i Qexl
 wi

(25)

k 1

ここで，n は選択節点数，Wi は，節点 i の流量配分比，Qexl は，層 l の全押出流量です．
前述した通り，完全発達流れの定式化では，押出流量の最も多い基準層 m に対する相対層 l の流量比は，
基準層の流線に沿って一定となります．流量配分比を設定しない場合（均一流量配分比の場合）には，流量比
は，解析モデル全体にわたって，解析条件として設定する各層の全流量の比率に等しくなります．従って，流量配分比を設定しない場合には，流量比分布の計算は省略され，解析結果として該当する情報も出力されません．一方，流量配分比を不均一に設定する場合には，(25)式で計算される各層の流量比を流入口境界条件として考慮し，流量比分布の計算が行われ，解析結果として流量比分布解析結果が出力されます．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

50以上で多層解析用メッシュの生成作業は終了です．最後に，ファイル/メッシュファイルエクスポートを
プルダウン選択して表示さっれるファイルセーブダイアログボックスにおいて適当なファイル名を設定
してメッシュ情報を保存します．多層解析用メッシュファイルの拡張子は，multimsh です．

図37 多層コートハンガーダイ解析用メッシュの保存
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved512) 多層押出解析条件の設定
完全発達多層流動解析プログラムを利用する場合，Solver selection パネルでFully developed ラジ
オボタンをONとします．発達多層流動解析と比較して収束性が良好で，数十回の非ニュートン反復計
算で収束解が得られます．また，発達多層流動解析で調整を要する圧力差寄与係数や界面計算緩和係数はダミーパラメータとなり，運用も容易です．

多層解析用メッシュ
ダミーファイル名（存在しないファイル名を指定）

非ニュートン反復計算回数は数十回の設定で十分

多層条件設定ボタン

ダミーパラメータ（1に設定）
Fully developed ラジオボタンをon

図38 完全発達多層流動解析の基本的な条件設定
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved52

多層押出解析で採用する押出量条件や材料物性ファイルは，多層条件設定ボタンを押して表示される下記フォー
ムに設定します．発達多層流動解析で調整を要するEvolution parameter もダミーパラメータとなります．Iteration 1,
Evolution parameter 1と設定し，追加ボタンを押してリストボックスにダミーパラメータを登録して下さい．物性データフ
ァイルと押出量を層数分設定した後，確定ボタンを押します．閉じるボタンを押して当ホームを閉じ，解析実行タブメ
ニューで保存ボタンを押して解析条件をファイル保存します．実行ボタンを押すことで完全発達多層流動解析プログ
ラムが起動し，多層熱流動解析が行われます．

物性データファイルを選択
設定（層数分）押出量物性を設定（層数分）

Evolution パラメータをダミーパラメータとして設定

図39 完全発達多層流動解析の条件設定フォーム
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved533) 多層熱流動解析結果のポスト処理
コンター図の作画方法は従来通り，表示層番号を設定し，作画内容をコンボボックスで選択します．流入口に流量比分布を設定した完全発達多層熱流動解析では，相対層（基準層を除く層）の流量比分布の図化出力項目が追加されました．基準層や発達多層流動解析では，流量比分布は出力されません．

作画項目の選択
表示層番号の設定

図40 コンター図の作画例（流量比分布コンター図）
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved54グラフ図の作画方法も従来通りです．主要な評価項目である流出口における計算量をグラフプロットするには，
2.5D多層流動解析結果(拡張子25dmultirst)ファイルを読み込み，グラフコントロールフォームで，Outlet ラジオボタン
をONとして，流出口の節点を一括選択します．出力層番号を設定し，作画項目をコンボボックスで選択します．X方向のラジオボタンをONとした後，描画ボタンを押すことで選択された作画項目がグラフプロットされます．

Outlet ラジオボタン

出力層番号の設定
作画項目選択コンボボックス

描画ボタン

Sort X方向ラジオボタン

図41グラフ図の作画例（流出口における表層1の層厚分布グラフ図）
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved55計算量の肉厚方向の分布をグラフプロットする場合は，マウスクリックラジオボタンをONとし，作画対象とする節点をマウスクリック選択します．作画項目選択コンボボックスで層～の作画項目を選択し，描画ボタンを押すことで，選択された作画項目の肉厚方向の分布がグラス図出力されます．

図42 グラフ図の作画例（センター，エッジ部における層流速分布グラフ図）
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

56

3.2.2 発達多層流動解析プログラム(Developing Solver）と完全発達多層流動解析プログラム(Fully Developed Solver) の連携運用
本節では，フィードブロックタイプコートハンガーダイをサンプルとして，多層流動解析プログラムの連携運用例を示します．解析対象は，下図に示す2.5D FEM モデルです．

100 mmCA

A’ ：多層流体合流位置

B

B’ ：モデル結合位置1000 mm

C’ ：流出位置
流出口リップクリアランス2 mm

2.5D FEM モデル
図43 フィードブロックタイプ多層コートハンガーダイ解析モデル
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

57

3D 化モデル

解析モデルは，図43内B-B’で連結されています．上流側モデルの解析に対しては，発達多層流動解析プログラム(Developing Solver)を適用しました．多層押出条件は，以下に示す通りです．第2法線応力差係数を設定し，上流側
フィードブロック流路内の界面包み込み現象を疑似的に表現しています．

図44 上流側モデル多層流動解析条件
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

58

また，上流側モデルでは，合流部断面の流路クリアランスを中間層4mm，表層1mmとしたflat 条件(メッシュファイル名:Upperfeedblock.multimsh)と下図に示す肉厚分布を設定したvariated条件(メッシュファイル名:Upperfeedblockvariated.multimsh)の比較解析を行いました．

表層フィードブロックプロファイル

中間層フィードブロックプロファイル

図45 Variated 条件の合流部肉厚分布（page 37,38参照）

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved59最初に，以下に示す計算コントロールデータを利用し，上流側モデルに対するflat 条件とvariated 条件の2ケースの解析を行いました．これらの解析用データ一式は，テストサンプルとしてご提供するフォルダー：testsample 内に保存されていますので，動作確認や別モデル作成時のご参考資料として下さい．

Flat条件：Upperfeedblock.fcal

Variated条件：Upperfeedblockvariated.fcal

図46 上流側モデル解析用計算コントロールファイルの内容
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

60次に，下流側モデルの完全発達多層流動の解析を行いました．Flat条件とvariated条件は，上流側モデルの解析結果として出力されるUpperfeedblock.fluxoutlet とUpperfeedblockvariated.fluxoutlet をマウスクリック選択するのみで自動的に連携解析が行われます．すなわち，これらのインターフェイスファイルの差を除き，解析条件やコートハンガーダイのメッシュ情報(Tdietest.multimsh)などの情報は両ケース共通です．完全発達多層流動解析プログラムは，上流側モデルで計算されている流出口の流束や温度分布を下流側モデルの流入口の熱流動条件として自動的に設定します．

Interface チェックボックス
をチェックInterface チェックボックス
をチェックInterface ファイル名：
Upperfeedblockvariated.fluxoutlet

Interface ファイル名：
Upperfeedblock.fluxoutlet

Variated条件：T-dievariated.fcal

Flat条件：T-dieflat.fcal

図47 下流側モデル解析用計算コントロールファイルの内容
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

61

図48に示す様に発達多層流れ定式化にPEMを導入することにより，流動(MD)方向に対して界面包み込み現象が促進される傾向や合流断面形状の変化が層厚分布に及ぼす影響が表現可能です．定性的に予測される通り，合流断面形状を変化させたVariated条件の解析結果では，界面包み込み現象によって表層エッジ部の表層厚が厚くなる傾向がFlat条件と比較して抑制されています．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

62

mm

表層

mm

Variated

Flat

中間層

Variated

Flat

図48 上流側フィードブロックモデル内の層厚分布予測結果の比較
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

63図49に完全発達多層流動解析プログラムを利用して計算された下流側コートハンガーダイ内の流量比分布の予測結果を示します．各モデルの流入口には，flat条件とvariated条件を採用した解析結果を利用して熱流動境界条件が設定されています．前述した通り，完全発達多層流動の定式化においては，基準層の流線に沿って流量比が一定になります（page 4,(4)式参照）．また，解析領域内の流量比分布の予測結果を利用して層厚分布の計算が可能です（page 4,(6)式参照）．完全発達多層流動解析プログラムの運用は容易です．また，計算時間も短く，これらの解析結果は，1ケース当たり30sec以内（非ニュートン反復計
算回数：20回）に得られます．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

64

Flat

Variated
図49 下流側コートハンガーダイモデル内の流量比分布予測結果の比較

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved65工学的に関心があるのは，流出口の多層フィルムの層厚分布です．シミュレーション結果では，ダイ内で流動状態にある層厚分布が予測されますので，一定速度で巻き取られる多層フィルムの層厚分布とは異なります．フィルム巻き取り時（フィルムキャスト）におけるネックインやエッジビーズによって多層フィルムの層厚分布は，変化しますが，
平面伸長変形が支配的なフィルム中央部では，質量保存則を利用してフィルム肉厚を流束（単位幅当たりの流量）に換算した情報が，シミュレーション結果との対比に適する情報と考えられます．流束は，実測位置でのフィルム層厚に巻き取り速度を乗じることで計算されます．
シミュレーションで予測される流束は，下図に示す様にインターフェイスファイルとして利用する拡張子fluxoutletの出力ファイル内に格納されています．
層数

出口で計算される各層流量(単位cc/s)

入口に設定される各層の流量と出口で計算される各層の流量を比較する
ことで質量保存則の表現精度を確認できます（非圧縮条件では，設定流量＝計算流量）

・・・

各層の設定流量(単位cc/s)

流出口に配置される節点数

節点流束＝節点流量/メッシュ幅（単位cm2/s)

層番号

節点のTD座標（単位cm)

節点流量（単位cm3/s)

節点流量を計算しているメッシュ
幅(単位cm)

図50 Fluxoutファイルの内容(T-dieflat.fluxoutlet)
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved66

拡張子fluxoutletの出力ファイルに出力されている節点のTD座標はcm単位です．モデル情報が直接出力されていますので，グラフ図の作画の際に適当にソートして下さい．また，流束の単位はcm2/sですので，必要に応じてSI単位に換算して下さい(座標は10倍でmm，流束は100倍でmm2/sに換算されます）．図50にFluxoutletファイルの情報をエクセルを利用してソート，SI単位換算した情報のグラフプロット図を示します．表層では，流束についても界面包み込み現象によってエッジ部が支配的になる傾向が示されています．
Variated条件では，この傾向が抑制され，最大流束と最小流束の差が低減し，多層フィルム層厚の均一性が向上する傾向が示されています．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved67Flux rate [mm2/s]

Flux rate [mm2/s]

TD position [mm]TD position [mm]

表層流束

中間層流束

図51 下流側コートハンガーダイモデル流出口における流束分布予測結果の比較

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved68本解析例では，上流側と下流側にモデルを分割し，それぞれのモデルに対して解析結果が出力されていますが，
Ver.2.0.0では，各解析結果を統合して表示する機能が実装されました．特に連携解析で得られた下流側の解析結果を先に読み込み，次いで上流側解析結果を読み込むと，上流側圧力解析結果に下流側モデルの流入口で計算されれている圧力が加算されます．非圧縮性解析では，一定圧力の加減操作によって，流動状態や界面形成状態は変化しないため，このような圧力更新が可能です．本解析結果の処理では，先に下流側コートハンガーダイモデルの解析
結果を読み込みます．解析結果ファイル名の右側に追加されたAppendチェックボックスをONとして，次に上流側フィードブロックモデルの解析結果を読み込むと下図に示す様に両モデルが統合表示されます．

Appendチェックボックス

図52 解析結果の統合表示（Append import）

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved69この統合表示機能を利用した温度と圧力の全体解析結果をそれぞれ図53と図54に示します．合流部断面形状の変化は，温度や圧力分布の予測結果にも影響を及ぼすことが示されています．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

70

MPa

Variated

Flat

図53 圧力分布解析結果の比較

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

71

℃

Variated

Flat

図54 温度分布解析結果の比較

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved72