Translate this MATLAB code to Excel VBA 

 

if app.DieselButton.Value $==$ true

if string$($app$.$ThrottleTypeSwitch.Value$)==$ "WOT"

if string$($app$.$StrokeSwitch.Value$)=="2"$

if app.NaturallyAspiratedButton.Value $==$ true

$\%--------$DIESEL WOTHROTTLE $2$STROKE NATURALLYASPIRATED

nst$=1$;

Vd$=($en$\_$vol$/$noc$)*(1/1000)$; $\%$displacement volume in m$^3$

Vc $=$ Vd$/($rc$-1)$; $\%$clearence volume in m$^3$

S$=((4*$Vd$)/(($bts$^2)*($pi$)))^(1/3)$; $\%$stroke length in m

B$=$bts$*$S; $\%$bore length in m

V$1=$ Vd $+$ Vc; $\%$v$1$ in m$^3$

V$2=$ Vc; $\%$v$2$ in m$^3$

P$2=$ P$1*$rc$^$k; $\%$p$2$ in kPa

T$2=$ T$1*$rc$^($k$-1)$; $\%$t$2$ in K

m$=($P$1*$V$1)/($R$*$T$1)$; $\%$m in kg

ma$=($af$/($af$+1))*(1-($exr$/100))*$m; $\%$mass of air in kg

mf$=(1/($af$+1))*(1-($exr$/100))*$m; $\%$mass of fuel in kg

mex$=($exr$/100)*$m; $\%$mass of exhaust in kg

T$3=(($mf$*$hv$*($nc$/100))/($m$*$cv$))+$T$2$; $\%$t$3$ in K

P$3=$P$2$; $\%$p$3$ in kPa

V$3=$V$2*($T$3/$T$2)$;

V$4=$V$1$;

T$4=$T$3*(($V$3/$V$4)^($k$-1))$; $\%$t$4$ in K

P$4=$P$3*(($V$3/$V$4)^$k$)$; $\%$p$4$ in kPa

T$\_1=$linspace$($T$1,$T$2,1000)$;

T$\_2=$linspace$($T$3,$T$4,1000)$;

Ap$=($pi$*($B$^2))/4$;

p$=1\mathrm{.}225$;

V$\_12=$linspace$($V$1,$V$2,1000)$;

V$\_34=$linspace$($V$3,$V$4,1000)$;

P$\_12=$P$1.*($V$1./$V$\_12).^$k;

$6$

P$\_34=$P$3.*($V$3./$V$\_34).^$k;

Qin$=$mf$*$hv$*($nc$/100)$;

Qout$=$m$*$cv$*($T$1-$T$4)$;

W$12=($m$*$R$*($T$2-$T$1))/(1-$k$)$;

W$34=($m$*$R$*($T$4-$T$3))/(1-$k$)$;

Wnet$=$W$12+$W$34$;

nt$=(1-($T$1/$T$2))*100$;

imep$=$Wnet$/($V$1-$V$2)$;

Wdoti$=(($Wnet$*($rpm$/60))/$nst$)*$noc;

Up$\_$avg$=2*$S$*($rpm$/60)$;

Wdotb$=($nm$/100)*$Wdoti;

Wi$=($Wdoti$*$nst$)/($rpm$/60)$;

Wb$=$Wi$*($nm$/100)$;

Wf$=$Wi$-$Wb;

Wdotf$=$Wdoti$-$Wdotb;

bmep$=($nm$/100)*$imep;

Torque$=(($bmep$*($en$\_$vol$/1000))/(2*$nst$*$pi$))*1000$;

BSP$=($Wdotb$/$Ap$)*$noc;

OPD$=($Wdotb$/$en$\_$vol$)*$noc;

bsfc$=(($mf$*($rpm$/60)*0\mathrm{.}5)/$Wdotb$)/$noc;

nv$=($ma$/($p$*$Vd$))*100$;

V$5=$V$1$;

V$6=$V$2$;

W$61=$P$1*($V$1-$V$6)$;

W$56=$P$1*($V$6-$V$5)$;

T$=[$T$1$ T$2$ T$3$ T$4$ T$1]$;

P$=[$P$1$ P$2$ P$3$ P$4$ P$1]$;

V$=[$V$1$ V$2$ V$3$ V$4$ V$1$ V$5$ V$6$ V$1]$;

TT$=[$T$2$ T$3]$;

V$23=[$V$2$ V$3]$;

P$23=[$P$2$ P$3]$;

TT$1=[$T$4$ T$1]$;

V$41=[$V$4$ V$1]$;

P$41=[$P$4$ P$1]$;

Bore $=0\mathrm{.}01*$B;

Stroke $=0\mathrm{.}01*$S;

a $=$ Stroke$/2$ ;

r $=$ con$\_$rod$*0\mathrm{.}01$ ;

r$2=$ r$/$a;

figure$(1)$

semilogy$([$V$1$ V$2],[$P$1$ P$1],\text{'}$c$\text{'},$V$\_12,$P$\_12,\text{'}$b$\text{'},$V$23,$P$23,\text{'}$r$\text{'},$V$\_34,$P$\_34,\text{'}$g$\text{'},$V$41,$P$41,\text{'}$y$\text{'})$

$7$

ylim$([1010^4])$

xlabel$(\text{'}$Volume $[$m$^3]\text{'})$

ylabel$(\text{'}$Pressure $[$Pa$]\text{'})$

tetha$2$a$=$linspace$(0,720,5)$;

figure$(2)$

plot$($tetha$2$a$,$P$)$

title$(\text{'}$Crank angle vs pressure'$)$

xlabel$(\text{'}$Crank Angle'$)$

ylabel$(\text{'}$pressure$\text{'})$

tetha$2$b$=$linspace$(0,720,8)$;

figure$(3)$

plot$($tetha$2$b$,$V$)$

title$(\text{'}$Crank angle vs Volume'$)$

xlabel$(\text{'}$Crank Angle'$)$

ylabel$(\text{'}$Volume$\text{'})$

U$=$ones$(1,720)$;

H$=$ones$(1,720)$;

tetha$=(1$:$1$:$720)$;

for tetha$2=1$:$1$:$720$

pistonspeed $=($pi$/2)*$sind$($tetha$2)*(1+($cosd$($tetha$2)/$sqrt$($r$2^2-($sind$($tetha$2)).^2)))*$Up$\_$avg;

pistonposition $=$ a$.*$cosd$($tetha$2)+$ sqrt$($r$^2-($a$^2).*$sind$($tetha$2).^2)$;

U$($tetha$2)=$ pistonspeed;

H$($tetha$2)=$ pistonposition;

end

figure$(4)$

plot$($tetha$,$U$)$

title$(\text{'}$Crank angle vs Piston speed'$)$

xlabel$(\text{'}$Crank Angle'$)$

ylabel$(\text{'}$Piston speed'$)$

figure$(5)$

plot$($tetha$,$H$)$

title$(\text{'}$Crank angle vs Piston Position'$)$

xlabel$(\text{'}$Crank Angle'$)$

ylabel$(\text{'}$Piston Position'$)$

app.OutputValuesTextArea.Value$=$sprintf$([\text{'}$ T$1=\%\mathrm{.}0$f K

T$2=\%\mathrm{.}0$f K

T$3=\%\mathrm{.}0$f K

T$4=\%\mathrm{.}0$f K

$\text{'}...$

$\text{'}$ P$1=\%\mathrm{.}0$f kPa

P$2=\%\mathrm{.}0$f kPa

P$3=\%\mathrm{.}0$f kPa

P$4=\%\mathrm{.}0$f kPa

$\text{'}...$

$\text{'}$ V$1=\%\mathrm{.}7$f m$3$

V$2=\%\mathrm{.}7$f m$3$

V$3=\%\mathrm{.}7$f m$3$

V$4=\%$