④不同位置的火花塞内置电阻对火花电流的影响
function Ig=huohuasaidianliu3(p,Vs,l,Ro,Co,aer,lg,lg1,lg2,R,R1,R2,Cr,Cq) dt=0.002;
t=-0.002:dt:80;Xo=110;
Cw=Co*l;Rw=Ro*l;K=(2*aer/p).*(Vs/lg).^2;
Ipr=1/4.*K.*Vs.*exp((K.*t-Xo)./2).*(1+exp((K.*t-Xo)./2)).^(-1.5); Ir=(Vs./R).*exp(-t./(R.*(Cq+(Cw+Cr)/2)));
Iw=(Vs/(Rw+R)).*exp(-t./((Rw+R)*(Co+Cw+Cr/2))); Ig=Ipr-Ir-Iw; K1=(2*aer/p).*(Vs/lg1).^2;
Ipr1=1/4.*K1.*Vs.*exp((K1.*t-Xo)./2).*(1+exp((K1.*t-Xo)./2)).^(-1.5); Ir1=(Vs./R1).*exp(-t./(R1.*(Cq+(Cw+Cr)/2)));
Iw1=(Vs/(Rw+R1)).*exp(-t./((Rw+R1)*(Co+Cw+Cr/2))); Ig1=Ipr1-Ir1-Iw1; K2=(2*aer/p).*(Vs/lg2).^2;
Ipr2=1/4.*K2.*Vs.*exp((K2.*t-Xo)./2).*(1+exp((K2.*t-Xo)./2)).^(-1.5); Ir2=(Vs./R2).*exp(-t./(R2.*(Cq+(Cw+Cr)/2)));
Iw2=(Vs/(Rw+R2)).*exp(-t./((Rw+R2)*(Co+Cw+Cr/2))); Ig2=Ipr2-Ir2-Iw2;
plot(t,Ig,'r-',t,Ig1,'b--',t,Ig2,'g-','linewidth',2); xlabel('t/us');ylabel('A');
legend('火花塞内置电阻距离放电电极7cm','火花塞内置电阻距离放电电极4cm','火花塞内置电阻距离放电电极3cm'); end
⑤不同阻值的火花塞内置电阻对火花电流的影响
function Ig=huohuasaidianliu2(p,Vs,l,Ro,Co,aer,lg,R,R1,R2,Cr,Cq) dt=0.002;
t=-0.002:dt:120;Xo=120;
Cw=Co*l;K=(2*aer/p).*(Vs/lg).^2;Rw=Ro*l;
Ipr=1/4.*K.*Vs.*exp((K.*t-Xo)./2).*(1+exp((K.*t-Xo)./2)).^(-1.5);
32
Ir=(Vs./R).*exp(-t./(R.*(Cq+(Cw+Cr)/2)));
Iw=(Vs/(Rw+R)).*exp(-t./((Rw+R)*(Co+Cw+Cr/2))); Ig=Ipr-Ir-Iw;
Ir1=(Vs./R1).*exp(-t./(R1.*(Cq+(Cw+Cr)/2)));
Iw1=(Vs/(Rw+R1)).*exp(-t./((Rw+R1)*(Co+Cw+Cr/2))); Ig1=Ipr-Ir1-Iw1;
Ir2=(Vs./R2).*exp(-t./(R2.*(Cq+(Cw+Cr)/2)));
Iw2=(Vs/(Rw+R2)).*exp(-t./((Rw+R2)*(Co+Cw+Cr/2))); Ig2=Ipr-Ir2-Iw2;
plot(t,Ig,'r-',t,Ig1,'b--',t,Ig2,'g-','linewidth',2); xlabel('t/us');ylabel('A');
legend('火花塞电阻为5k欧姆','火花塞电阻为15k欧姆','火花塞电阻为20k欧姆'); end
⑥火花电流
function Ig=huohuasaidianliu01(p,Vs,l,Ro,Co,aer,lg,R,Cr,Cq) Xo=120;
dt=0.002;t=-0.002:dt:120; Cw=Co*l;
Rw=Ro*l;K=(2*aer/p).*(Vs/lg).^2;
Ipr=1/4.*K.*Vs.*exp((K.*t-Xo)./2).*(1+exp((K.*t-Xo)./2)).^(-1.5); Ir=(Vs./R).*exp(-t./(R.*(Cq+(Cw+Cr)/2)));
Iw=(Vs/(Rw+R)).*exp(-t./((Rw+R)*(Co+Cw+Cr/2))); Ig=Ipr-Ir-Iw;
plot(t,Ig,'r-','linewidth',2);
33