流体工学 終末速度 ホームページプログラミング 


非線形方程式を逐次代入法で解く。 termvel.html

スクリプト部分の記述

(VBScript の式の部分)
 <script language="VBScript"><!--
Dim dpm,dp,rof,rop,mu,ZEP,ZX0,NL,X,ZX1,ZX2,ZXD,ZX3,Rep,vt,CR,G
Sub A_onClick
dpm=Valdpm.value
dp=dpm/1000
Valdp.value=dp
rof=Valrof.value
rop=Valrop.value
mu=Valmu.value
NL =1 
ZEP = .000001
ZX0=100
ZXD=1
Do While ABS(ZXD) > ZEP
X = ZX0
velocity
ZX1 = G
X = ZX1
velocity
ZX2 = G
ZXD = ZX2 - 2 * ZX1 + ZX0
ZX3 = ZX2
ZX3 = ZX0 - (ZX1 - ZX0) * (ZX1 - ZX0) / ZXD
ZX0 = ZX3
NL=NL+1
Loop
Valvt.value=ZX3
ValRep.value=Rep
ValNL.value=NL
End Sub
Sub velocity
vt=X
Rep=vt*dp*rof/mu
If Rep>0.035 Then
CR=(0.55+4.8/(Rep^0.5))^2
Else
CR=24/Rep
End If
G=(4*9.807*dp*(rop-rof)/(3*rof*CR))^0.5
End Sub
--></script>