流体工学 終末速度 ホームページプログラミング
非線形方程式を逐次代入法で解く。 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>