10 OPEN "PREME. DAT" FOR OUTPUT AS位
20 'こここここここここ=(Program Name :PACRYME)==二二=二一
30 '一一一一一一一一一( r-lll1容積理論による拡散係数の推算)一一一
50 '本木本本木本木本本木本木本木本木本本本木本木本木本木本木本本木本字本木本木本木本木本本木本木本木本木本本本木本本本木本本本本木本本木'
60 '木本木本木本 nrun =山ßUJく準の数 料本'
70 '木本料水木 m = パラメータの数 料本'
80 '本木林本木 n = データの数 料本'
90 '本木本木林 Bl，B2， B3， B4 =溶媒の比容積定数(m3/kg単位でう 料本' 100 '本木林林 C11， C12， C13 =高分子の比容積定数(m3/kg単位で、) 料本' 110 '本木林料 KAI = Floryのχノマラメータ[-J 料本' 120 '本木林料 V1HATS = O[KJで、の溶媒の比容積[m3/kgJ 料本' 130 '料本木林 V2HATS = O[KJで、の高分子の比容積[m3/kgJ 料本' 140 '本木林料 GAJ.伽1 =溶媒のオーバーラップファクター[-] 林本' 150 '本木林料GAMM2 = 高分子のオーバーラップファクター[-J 料本' 160 '本木林料 DATD(i) =拡散係数データ[m2/sJ 林本'
170 '料本木林 DAIT(i) =温度データ[KJ 料本'
180 '料本木林 DATW1(i) =溶媒の重量分率データ[-] 料本'
190 '料本木林 DO =定数Do 料本'
200 '本木林料 ED = エネルギーノミラメータ E 木林'
210 '料本木林GUZAIこど 本木本'
220 '料本木林 RG = ガス定数(=8.3144 [J/molKJ) 本料'
230 '料本木林 添字l 二溶媒 2 = 高分子 本木本'
240 '本木林料 CALT = 相関結果をグラフに表示する時の温度[KJ 本木本'
250 '料料料 M1 =溶媒の分子量 本木本'
260 '料本木林 ECOH 二溶媒の凝集エネルギー口/molJ 料本' 270 '本木本木林 VB =溶媒のがI�点分子容[m3/molJ 本木本' 280 '木本本木本木本木本木本木本木本木本本木本木本木本木本木本本木本木本本本本木本本木本木本本木本木本本木本木本本木本木本木本本本木本本本' 290 '一一一一一一一一一一一一一一(自己列の宣言)一一一一
300 DIM DATD(200)，CALD(200)， GCALD(100， 100)， DATDL(200)，DATDH(200) 310 DIM DATW1(200)， CALW1(200)， DATW1L(200)，DATW1H(200)
320 DIM DA汀(200)，CALT(10)，DAη1'(200) 330 DIM G1 %(100)， G2%(100)， RK%(20)
360 RG=8. 3144
370 '一一一一一一一一一(データの読み込み)一一一一一 380 READ NR閃
390郎AD N
400 READ M1， ECOH， VB 410 READ B1， B2， B3，B4 420 READ C11，C12，C13 430 READ V2HATS，GA刷2 440 FOR 1二1 TO NRUN 450 READ CALT(1) 460 NEXT 1
470 FOR 1=1 TO N
480 READ DATD(1)，DATW1(1)，DAη(1) 490 DAm(1)二DATT(1)
500阻XT 1
510 '--_一一一一一ー一一一(ノミラメータの計算)一一一 520 DO=l. 835E-10刈1-7.753E-9
530 ED=O. 1252本ECOH-2537 540 KA1 =2. 901-7. 160E-5本ECOH 550 GUZAI=l. 671E3本VB十O.2133
560 V1HATS=(0.5955E3本VBtO.011912)/M1 570 GA刷1=0.003391E似VB十O. 2581
580 '一一一一一一一一一一一(設定パラメータの印刷 )---590 PR1NT"Program Name : PACRYME"
600 LPR1NT"Program Name : PACRYME"
610 PRINT 620 LPR1NT
630 PR1NT"System:Acryl adhesivetMethyl acetate"
640 LPR1NT"System:Acry ladhesivetMethyl acetate"
650 PR1NT 660 LPRINT
670 PR1NT"相関式の型D二00*(1-ゆ1)�2*(1-2本χ本ゆ1) 本exp(-(wl本Vlsト打w2*V2s)/VFH本γ)) "
680 LPR1NT"相関式の型D二Do*(l-Ø 1)�2*(1-2本χ本ゆ1) 本exp(-(wl*V1s十ど本w2本V2s)/VF取り) "
690 PRINT 700 LPR1NT
710 PRINT"データの数n= ";N :PRINT 720 LPRINT"データの数n= ";N :LPRINT
730 PRINT"溶媒の分子量M1=";M1"溶媒の凝集エネルギ- ECOH二";ECOH :PRINT 740 LPRINT"溶媒の分子量M1=";M1"溶媒の凝集エネルギーECOH=";ECOH:LPRINT 750 PRINT"溶媒の沸点分子容VB=";VB :PRINT
760 LPRINT"溶媒の沸点分子容VB=";VB :LPRINT
770 PRINT"溶媒の比容積定数B1=";B1" B2= ";B2" B3= ";B3" B4二";B4 780 LPRINT"溶媒の上防予積定数B1=";B1" B2= ";B2" B3= ";B3" B4= ";B4 790 PRINT"高づナ子の比零積定数C11=";C11" C12ご";C12" C13=";C13 800 LPRINT"高分子の比容積定数C11=";C11" C12二 ";C12"C13=";C13 810 PRINT"溶媒のO[KJでの比容積V1HATS 二
V1HATS" 高分子のO[KJでの比容積V2臥TS = ";V2HATS 820 LPRINT"溶媒のO[KJでの比容積 V1HATS 二
V1HATS" 高分子のO[KJでの比容積V2HATS二";V2臥TS 830 PRINT"溶媒のoverlap factor GA削1 = ";GAMM1
， 高分子のoverlap factor 印刷2二 ";GA紬12 840 LPRINT"溶媒のoverlap factor GA1凹1 = ";GAMM1
' 高分子のoverlap factor GAMM2 = ";GA1凹2 850 PRINT"χパラメータKAI=" ; KAI : PRINT 860 LPRINT"χノミラメータKAI二"; KAI : LPRINT
870 PRINT"ノミラメータDO=";DO"ノミラメータE=";ED :PRINT 880 LPRINT"ノミラメータDOニ";DO"ノミラメータE二";ED :LPRINT 890 PRINT"ノミラメータ f GUZAI =" ; GUZAI : PRINT
900 LPRINT"ノミラメータ f GUZAIニ";GUZAI :LPRINT 910 INPUT" Yキーを押して下さし\0 " ; Y$
920 ' 一一一一一一一一一一一一一一( 拡散係数の計算 〉一一一 930 PRINT " No. T[KJ W1 [ -J D. exp[m2/ s J
D. cal[m2/sJ Error(%)"
940 LPRINT " No. T[KJ W1[-J D. exp[m2/sJ
D. cal[m2/sJ Error( %)"
950 ' 一一一一一一( 惟算されたパラメータでの拡散係数の計算 )一一 960 ER1=0
970 FOR 1=1 TO N
980 V1=1/(B1-B以(DATT(I)-273. 15)-B3/(B4-(DAl寸(1)-273. 15))) 990 V2二1/(C11+C12本(DAη(1)-273.15)十C13本(DATT(I)-273. 15)�2)
1020 F刊かDATW1(1)本(V1-V1HATS)/G川仰い(1-DATW1(1))本(V2-V2臥TS)/G AMM2 1030 Q1=D01*(1-FA11)円以(1-2*KA1*FA11)
1040 CALD(1)ニQ1本EXP(-(DATW1(1)本V1臥TS十GUZA1*(1-DATW1(1))本V2HATS)/F刊の 1050 ERO=(DATD(1)-CALD(1))/DATD(1)本100
1060 ER1=ER1十店S(ERO)
1070 PR1NT US1NG " 桝井 柑井.柑 拡柑井 井井.非井押井一
柑.柑非材一… 柑.柑ザー�"; 1， DAη(1)， DATW1(1)，DATD(1)，C札D(1)，ERO 1080 LPR1NT US1NG " 桝井 柑井.柑 拡桝井 柑.柑非井井……
柑.柑柑井…… 柑.仰が…�";1， DAη(1)， DATW1(1)，DATD(1)，C札D(1)，ERO 1090 WR1TE非 2，DATT(1)，DATW1(1)，DATD(1)， CALD(1)，ERO
1100 NEXT 1 1110 TER=ER1/N 1120 PR1NT "
1130 LPR1NT "
1140 PR1NT US1NG"average error(%)二柑拡非開……";TER 1150 LPR1NT US1NG"average error(%)=柑北井柑……" ;TER 1160 PR1NT
1170 LPR1NT 1180 CLOSE
1220 '一一一一一一一一一一( グラフ表示用の拡散係数の計算 )一一一
1230 '一一一一一一一一一一(データの最大、 最小値のチェッ ク)一一一
1240 FOR 1=1 TO N 1250 DATW1H(1)=DATW1(1) 1260 DATW1L(1)=DATW1(1) 1270 DATDH(1)二DATD(1) 1280 DATDL(1)=DATD(1) 1290肥XT 1
1300 FOR 1二1 TO N-1
1310 1F DATW1L(1)くこDATW1L(1十1) THEN DATW1L(1十1)ニDATW1L(1) 1320悶XT 1
1330 FOR 1=1 TO N-1
1340 1F DATW1H(I)>二DATW1H(1十1) THEN DATW1H(I十1)=DATW1H(1) 1350悶XT 1
1360 FOR 1=1 TO N-1
1370 1F DATDL(1)く二DATDL(It1) THEN DATDL(1t1)=DATDL(1) 1380阻XT 1
1390 FOR 1=1 TO N-1
1400 IF DATDH(I)>二DATDH(It1) THEN DATDH(It1)=DATDH(I) 1410 NEXT 1
1420 W1胤XD=DATW1H(N) : W1MIND=DATW1L(N) 1430 D臥XD=DATDH(N) : DMIND=DATDL(N)
1440 PRINT" W1min二 " ;W1MIND"[-] W1max = " ;W1刷XD" [-] "
1450 LPR1NT" W1min = ";W1MIND"[-] W1max二" ;W1MAXD" [一 ]"
1460 PRINT" Dmin二" ;DM1ND"[-J Dmax二" ;D胤XD"[-] "
1470 LPRINT" Dmin = " ;DMIND" [ -J Dmax = " ;D胤XD" [一 ] "
1480 PR1NT 1490 LPR1NT
1500 1NPUT" グラフ表示をします。 yキーを押して下さし\0 " ; Y$
1510 1F Y$ご" Y" OR Y$=" y" THENωTO 1520 ELSE 1500
1520 1NPUT" 酢酸メチルの重量分率の表示の最小値を数値で、入れて下さい";W1M1N 1530 1NPUT" 酢酸メチルの主主分率の表示の技大他を数値で、入れて下さい" ;W1胤X 1540 '本木本木本本木本木本木本木本木本木本本木本木本木本本木本木本木本木本本木本本本本木本木本本木本本木本本木本木本本木本本木本本*本木本' 1550 INPUT" 拡散係数のぷ示の段小値を数値で-入れて下さい" ;DM1N
1560 1NPUT" 拡散係数の表示の最大値を数値で、入れて下さい";D臥X 1570 '一一一一一一一一一一一(グラフ表示用の拡散係数の計算〉一一 1580問=100
1590 GH=CW1胤X-W1M1N)/NH 1600 FOR ]=1 TO NRUN 1610 FOR 1=0 TO問 1620 CALW1(1)ごGH本I十W1MIN 1630 DATT(1)ごCALT(])
1640 V1=1/CB1-B2本(DAIT(I)-273.15)-B3/(B4-(DAIT(1)-273. 15))) 1650 V2=1/(C11tC12本(DAIT(I)-273.15)十C13本(DAIT(I)-273. 15)�2) 1660 D01=DO本EXP(-ED/RG/DAl寸(1))
1670 FAl1=CALW1(I)本V1/(CALW1(I)本V1t(1-CALW1(I))本V2)
1680 FVHG=CALW1(I)本(V1-V1HATS)/G品仰い(1-CALW1(I))本(V2-V2HATS)/GA削2 1690 Q1=D01*(1-FAl1)�2*(1-2*KAI*FAl1)
1700 GCALD(I， ])=Q1*EXP(-(CALW1(I)本V1HATS十GUZAI*(1-CALW1(I))本V2HATS)/FVHG) 1710肥XT 1
1720悶XT ] 1730 SCREEN 3 1740 C凶l
1770 SX1=180 :SY1二300 :SX2=580 : SY2=50 1780 DX=50
1790 XSCALE二(SX2-SX1)/(W1臥X-W1M1N) :YSCALE=(SY1-SY2)/(F1X(L∞(D臥X) /L∞(10非))-F1X(L∞(DM1N)/LOG(10非))十1)
1800 '一一一一一一(グラフの軸表示 ，スケール表示)--1810 L1NE(SX1， SY1)-(SX2，SY2)，4，B
1820 GET@(10， 10)一(10，20)，G1%
1830 GET@(10， 10)-(20， 10)， G2%
1840 NXS=F1X((SX2-SX1)/DX) : NYSニF1X(LOG(D胤X)/LOG(10非))-F1X(LOG(DM1N) /L∞(10井))は
1850 DY=(SY1-SY2)/NYS 1860 FOR 1=0 TO NXS
1870 PUT@(SX1tDX本1，SYl-5)， G1%，PSET， 0，4 1880 NEXT 1
1890 FOR L=l TO NYS 1900 FOR-1二1 TO 10
1910 DYY=(LCお(10*1)/LC氾(10井)-1)本DY 1920 SYYY=SY1-(DYY十DY*(L-1)) 1930 PUT@(SX1，SYYY)，G2%，PSET， 0，4 1940 NEXT 1 :NEXT L
1950 '一一一一一一一一一一(実験値のプロット)一一
1960 1F LCお(DM1N)/LOG(10非)くF1X(LOG(DM1N)/LCお(10井)) THEN NYSS=
F1X(LC氾(DM1N)/LOG(10非))-1
1970 IF LOG(DMIN)/LCお(10非)=FIX(LCお(DMIN)/LOG(10井)) THEN NYSS=
F1X(LOG(DM1N)/LCお(10井)) 1980 111=1
1990 FOR 1=1 TO N
2000 X=(DATW1(1)-W1M1N)ばSCALEtSX1 : Y=SY1-(LCお(DATD(1))/LOG(10約一NYSS) 本YSCALE
2010 CIRCLE(X， Y)，4， 7 : PAINT(X， Y)， 111，7 2020 1F 1=NωTO 2040
2030 IF DA廿T(1) く DA廿T(1十1) THEN 111二111t1 2040 NEXT 1
2050 '一一一一一一一一一一
(決定したパラメータによる計算線を描く)--2060 FOR ]=1 TO NRUN 2070 FOR 1=1 TO問
2080 XX(I， ])=(CALW1(1)-W1MIN)刈SCALEtSX1 :YY(I， ])=SY1-(L∞(GCALD( 1， ]))
/L∞(10約一NYSS)本YSCALE 2090 NEXT 1
2100 NEXT J
2110 FOR Jご1 TO NRUN 2120 FOR 1=1 TO NH-1
2130 L1NE (XX(I， J)， YY(I， J))-(XX(I+1， J)， YY(I十1，J))， 4 2140 NEXT 1
2150 NEXT J
2160 FOR 1=0 TO NXS
2170 XW1(I)=(W1MAX-W1M1N)/NXS本I+W1M1N 2180 NEXT 1
2190 '一一一一一一一一一(機軸の目盛の表示)一一一
2200 FOR 1=0 TO NXS 2210 NS=XW1(1) 2220ωSUB 2750 2230 FOR.K=l TO LNS 2240 KX=LNS本8-K刈
2250 PUT@(SX1+DX本1-KX， SY1 +20)，臥NJI(KCODE(K)) 2260 NEXT K
2270 NEXT 1
2280 '一一一一 一一一一一一(縦軸の目盛の表示)
---2290 FOR 1=0 TO NYS 2300 NS=10�(NYSS+1) 2310ωSUB 2750 2320 FOR K二1 TO LNS 2330 KY二LNS本8-K本8
2340 PUT@(SX1-KY-20，SY1-DY本1-10)，臥NJ1(KCODE(K)) 2350 NEXT K
2360 NEXT 1
2370 '一一一一一一一一一(横軸の説明)一一一一一一一一一一一一
2380 A$="W1[ -]"
2390 LNS=LEN(A$) 2400 FOR ト1 TO LNS 2410 CC=ASC(MID$(A$， 1， 1))
2420 PUT@((SX1十SX2)/2-LNS本8+1本8+10，SY1十40)，臥NJI(CC)
2450 LNS=LEN(A$) 2460 FOR 1=1 TO LNS 2470 CC=ASC(M 1 D$(A$， 1 ，1))
2480 PUT@((SX1tSX2)/2-LNS本8十I本8t70，SYl-280)，臥NJ1(CC) 2490悶XT 1
2500 '一一一一一一一一一一一(横軸に漢字を表示する)一一一一 2510 AA$="質量分率"
2520 KLNS=KLEN(AA$， 2) 2530 FOR 1 =1 TO KLNS
2540 AAA$=J 1 S$(KMID$(AA$， 1 t 1， 1)) 2550 AB=V札("朗勺AAA$)
2560 PUT@((SX1tSX2)/2-KLNS本16十i本16-50，SY1ト40)，KANJ 1(AB)，PSET，2， 0 2570 NEXT 1
2580 '一一一一一一一一一(縦軸に漢字を表示する)一一一一一
2590 BB$二"相互拡散係数"
2600 YAXIS$=BB$
2610阻.3Y二KLEN(YAX 1 S$，2) 2620 FOR 1 =1 TO阻.3Y 2630ωSUB 2820
2640 PUT@(SX1-100， (SY1tSY2)/2t阻.3Y*16-16刈)，RK%， PSET， 2， 0 2650 NEXT 1
2660 '一一一一一一一一一一一(縦軸に1 byte文字を書く)一一一
2670 CC$二"D[m2/s]"
2680 YAX1 S$二CC$
2690阻.3YY= LEN (Y AX 1 S$) 2700 FOR 1 =1 TO KLSYY 2710ωSUB 2980
2720 PUT@(SX1-100， (SY1tSY2)/2-阻.3YY本8-8料ト50)，RK%，PSET， 7， 0 2730 NEXT 1
2740 END
2750 '一一一一一一一一一一(数字を漢字コー ドに変換するサフゃルーチン)一一一
2760 NS$=STR$(NS) : LNS二LEN(NS$)
2770 1 F NS)=O THEN LNS=LNS-1 : NS$ニRIGHT$(NS$，LNS) 2780 FOR K=l TO LNS
2790 KCODE(K)=ASC(MID$(NS$，K， 1)) 2800 NEXT K
2810郎TぼN
2820 '一一一一一一一一一一(文字を90度回転するサフ守ルーチン)一一
2830 '一一一一一一一ー一一一(2byte系全角文字の場合)一一一一一一一一 2840 K$=J IS$(KMID$(Y AXIS$， 1 t 1， 1))
2850 PUT@(O， 0)， KANJI(VAL("組"tK$))，PSET 2860 RK9Vo(0)=16 :郎9Vo(1)=16
2870 FOR XX=15 TO 0 STEP -1 2880 Sl=O :S2=0
2890 FOR YY=O TO 8
2900 W1=POINT(XX， YY) :昭二POINT(XX，YY十8) 2910 IF W1二7 THEN Sl=Slt2�(7-YY)
2920 IF W2=7 THEN S2=S2t2�(7-YY) 2930 NEXT YY
2940 RK9Vo(17-XX)二VAL("組"十郎X$(256本S2tS1)) 2950 NEXT XX
2960 PUT@(O， 0)， KANJI(VAL("組"十K$))，PSET， 0， 0 2970郎T四N
2980 '一一一一一一一一一一(文字を90度回転するサブルーチン)一一
2990 '一一一一一一一一一一一(lby.te系文字の場合)一一一 3000 KYY=ASC(KMID$(YAXIS$， 1， 1))
3010 PUT@(O， 0)，臥NJI(KYY)，PSET 3020 RK9Vo(0)=16 :郎9Vo(1)=8 3030 FOR XX=7 TO 0 STEP -1 3040 Sl=O :S2=0
3050 FOR yy=o TO 7
3060 W1二POINT(XX，YY) : W2=POINT(XX， YYt8) 3070 IF W1=7 THEN Sl=Sl十2�(7-YY)
3080 IF W2=7 THEN S2=S2t2�(7-YY) 3090限XT YY
3100 RK9Vo(9-XX)=VAL("組"十HEX$(256本S2tS1)) 3110 NEXT XX
3120 PUT@(O， 0)，臥NJI(VAL("組"十K$))，PSET，0， 0 3130肥T限N
4600 DATA 1
4602 ' 一一一一一一一一一一一一一一(データの数)--4605 DATA 17
4610 '一一一一一一(溶媒の分子量， 沸点分子容， 凝集エネルギー)一一 -4620 DATA 74.08，26390，8. 31E-5
4670 '一一一一一一一一一一一一(比容積定数 溶媒)一一 一一一 4680 DATA 1012. 8， 1. 08， 15000， 280
4690 '一一一一一一一一一一一(比容積定数
高分子)一一一一一一一一一一一一一-4700 DATA 1000， -7. 257E-1， 6. 09E-4
4710 '一一一一(高分子の臨界比空孔容積， オーバーラップファクター)一一一一一一'
4720 DATA O. 955E-3， 1.300
4750 '一一一一一一一一一一一(グラフの表示温度)一一一 4760 DATA 313. 15
4770 '一一一一一一一一一一一一(拡散係数，質量分率，温度)一一一一一 4780 DATA 1. 79E-11， O. 0127， 313.15
4790 DATA 2. 09E-11， O. 0173， 313.15 4800 DATA 2. 13E-11， O. 0222， 313. 15 4810 DATA 2. 38E-11， O. 0249，313.15 4820 DATA 2. 51E-11， 0.0304，313.15 4830 DATA 2. 82E-11， O. 0337， 313.15 4840 DATA 2. 82E-11， O. 0340， 313.15 4850 DATA 3. 06E-11， O. 0400， 313.15 4860 DATA 3. 16E-11， O. 0443， 313.15 4870 DATA 3. 66E-11， O. 0558， 313.15 4880 DATA 4. 77E-11， O. 0784， 313.15 4890 DATA 5. 55E-11， O. 0888， 313.15 4900 DATA 5. 89E-11， 0.1010，313.15 4910 DATA 5. 87E-11， 0.1060，313.15 4920 DATA 6. 41E-11， O. 1180，313. 15 4930 DATA 6. 66E-11， O. 1200，313. 15 4940 DATA 7. 28E-11， 0.1310，313.15
謝 辞
本研究を遂行するにあたり、終始適切なご指導ご鞭捲を賜りました九州大学 大学院工学研究科化学システム工学専攻 荒井康彦教授に深く感謝の意を表 します。本論文の審査にあたり、数々の有益なご助言を賜りました九州大学大 学院工学研究科化学システム工学専攻 船津和守教授、九州大学大学院工学研 究科材料物性工学専攻 梶山千里教授、九州大学大学院工学研究科応用物理学 専攻 根本紀夫教授に深く御礼申しあげます。
本研究を行うにあたり、的確なご助言とご指導を賜りました山口大学工学部 佐野雄二名誉教授、山口大学工学部応用化学工学科 山本修一助教授に深く感 謝し1たします。また終始お励まし頂いた山口大学工学部応用化学工学科 熊田 誠教授に深く感謝し1たします。
本研究に対して適切なご助言、多大なご協力を頂きました宇部工業高等専門 学校物質工学科 福地賢治教授、九州大学大学院工学研究科化学システム工学 専攻 岩井芳夫助教授に深く感謝いたします。
本研究を遂行するにあたり、 アクリル粘着剤試料の提供と貯蔵弾性率の測定 をして頂いた積水化学工業株式会社 堀芳樹氏に深く感謝し、たします。さらに、
積水化学工業株式会社時代から今日まで終始お励まし頂いた古田武人氏に感 謝いたします。
本論文の資料作成にご協力頂きました九州大学大学院工学研究科化学シス テム工学専攻 米津節子助手に深く感謝し1たします。 また、会議室利用などで 大変お世話になりました九州大学大学院工学研究科化学システム工学専攻化 学プロセス教室の事務室の皆様に感謝し1たします。
本研究の研究クゃループとして、測定ならびに計算にご協力を頂いた山口大学 工学部応用化学工学科の卒業生諸君に感謝いたします。
最後に、健康管理に気を配ってくれた妻 留美子と両親ならびに心の支えと なってくれた長女 美保に感謝し1たします。