致谢
27
参考文献
[1]邹逢兴，张湘平等.电磁兼容技术[M].北京：国防工业出版社.2005 [2]国家标准 GB/T4365-2003.电工术语电磁兼容
[3]Canavero F.Numerical simulation for early EMC design of cars [C].4th conference on electromagnetic compatibility.Brugge.Beligum.2000.9:32-39
[4]Neumayer R.Stelzer A.Continuous simulation of system-level automobile EMC problems[C].IEEE International Symposium on Electromagnetic Compatibility.2003.(8):409-413
[5]Jens Benecke，Stefan Dickmann Analytical HF Model for Multipole DC Motors[C]. Proceedings of the 18th international Zurich Symposium on Electromagnetic Compatibility，EMC.2007.1:201-204
[6]吴元杰，徐博候，徐航.火花塞点火过程非线性模型及参数识别.内燃机学报.1999 [7]汪泉弟，刘春艳，俞集辉.汽车火花点火系统电磁干扰的抑制方法[J].重庆大学学报(自然科学报).2007.30（7）：46-49
[8]LI Xu,YU Ji-hui,WANG Quan-di,LI Yong-ming.Noise Spark Current cause by Capacitive Discharge in an Automobile Secondary Ignition System[C].IEEE international symposium on Microwave,Antenna,Propagation and EMC Technologies.2007. 8:1319-1322
[9]秦曾衍，左公宁，王永荣，吴泓，孙广生，孙鹞鸿.高压强脉冲放电及其应用.北京：北京工业大学出版社.2000
[10](苏)米夏兹TA.大功率毫微秒脉冲的产生。北京：原子能出版社.1982 [11]江贤祚.天线原理[M].北京：北京航空航天大学出版社.1993 [12]赵福堂.汽车电器与电子设备。北京：北京理工大学出版社.1997 [13]关文达.汽车构造.北京：机械工业出版社.2005 [14]刘新亮.汽车电器的电磁兼容.汽车电器.2001 [15]邱关源，罗先觉.电路.北京：高等教育出版社.2006 [16]周守昌.电路原理.北京：高等教育出版社.1999
[17]Fujiwara O.Calculation of Ignition Noise Level Caused by Plug Gap Breakdown. IEEE Trans on Electromagnetic Compatibility.1982,24(2):26-32 [18]姜述刚等.汽车电磁兼容性综述.汽车技术.1996
[19]刘昆敏.汽车点火系电磁干扰的抑制.武汉汽车工业大学.1995
28
[20]张春化，赛小平.汽车电器与电路.北京：人民邮电出版社.2003 [21]关文达.汽车构造.北京：机械工业出版社.2005
[22]徐立.我国汽车电磁兼容技术发展状况.安全与电磁兼容.2003 [23]马西奎.电磁场理论及应用.西安：西安交通大学出版社.2000 [24]王长清.现代计算机电磁学基础.北京.北京大学出版社.2005
29
附录
?不同阻值的火花塞对比
function [G1,G2]=myfile1(R1,R,C,Cr,Cq,Cs,Ro,Co,Cp) dt=1;f=0.1:dt:1000; Cu=Cq+Cr+Cs; Wrq=Co/(Ro*Cu.^2); Wq=1/(Cu*R);Wd=1/(C*R); Wr=Co/(Ro*Cs.^2);W=2*pi*f;
K=(1+W.^2./(Wd.^2)).*(Wr./(2.*W)+(1-sqrt(Wr./(2.*W)))).^2;
M=(1+(W./Wd).*sqrt(Wrq./(2.*W))).^2+(1+sqrt(Wrq./(2.*W))).^2.*(W./Wq).^2; G1=(Cs/(Cr+Cp)).*sqrt(K)./sqrt(M); Wq1=1/(Cu*R1);Wd1=1/(C*R1);
K1=(1+W.^2./(Wd1.^2)).*(Wr./(2.*W)+(1-sqrt(Wr./(2.*W)))).^2;
M1=(1+(W./Wd1).*sqrt(Wrq./(2.*W))).^2+(1+sqrt(Wrq./(2.*W))).^2.*(W./Wq1).^2; G2=(Cs/(Cr+Cp)).*sqrt(K1)./sqrt(M1); plot(f,G1,'b-',f,G2,'g--','linewidth',2); axis([0,130,0,0.05]);
xlabel('f/MHz');ylabel('G(w)');
legend('火花塞电阻为10k欧姆','火花塞电阻为5k欧姆'); end
?电阻型火花塞与电感型火花塞对比
function [G1,G2]=myfile(R,L,C,Cr,Cq,Cs,Ro,Co,Cp) dt=1;f=0.1:dt:1000;
Cu=Cq+Cr+Cs;Wrq=Co/(Ro*Cu.^2); Wq=1/(Cu*R);Wd=1/(C*R); Wr=Co/(Ro*Cs.^2);W=2*pi*f;
K=(1+W.^2./(Wd.^2)).*(Wr./(2.*W)+(1-sqrt(Wr./(2.*W)))).^2;
M=(1+(W./Wd).*sqrt(Wrq./(2.*W))).^2+(1+sqrt(Wrq./(2.*W))).^2.*(W./Wq).^2; G1=(Cs/(Cr+Cp)).*sqrt(K)./sqrt(M);
A=(1-W.^2.*(Cu*L)+Cu*R.*sqrt(Wrq.*W./2)+W.^2.*Cu.*L.*sqrt(Wrq./(2.*W))).^2;
30
B=(Cu*Ro+Cu*R.*sqrt(Wrq.*W./2)-W.^2.*Cu.*L.*sqrt(Wrq./(2.*W))).^2; N=((1-L*C.*W.^2).^2+C.^2*R.^2.*W.^2).*(Wr./(2.*W)+(1-sqrt(Wr./(2.*W)))).^2; G2=(Cs/(Cr+Cp)).*sqrt(N)./sqrt(A+B); plot(f,G1,'r-',f,G2,'b--','linewidth',2); xlabel('f/MHz');ylabel('对比');
legend('电阻型火花塞','电感型火花塞');axis([0,1000,0,0.25]); end
?火花塞和高压点火线的内置电阻对火花电流的影响
function Ig=huohuasaidianliu0(p,Vs,l,l1,l2,Ro,Co,aer,lg,R,R1,R2,Cr,Cq) dt=0.002;
t=-0.002:dt:200;Xo=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; Cw=Co*l1;Cw=Co*l2;
Rw1=Ro*l1;K=(2*aer/p).*(Vs/lg).^2;Ir1=(Vs./R1).*exp(-t./(R1.*(Cq+(Cw+Cr)/2))); Iw1=(Vs/(Rw1+R1)).*exp(-t./((Rw1+R1)*(Co+Cw+Cr/2))); Ig1=Ipr-Ir1-Iw1;
Rw2=Ro*l2;K=(2*aer/p).*(Vs/lg).^2;Ir2=(Vs./R2).*exp(-t./(R2.*(Cq+(Cw+Cr)/2))); Iw2=(Vs/(Rw2+R2)).*exp(-t./((Rw2+R2)*(Co+Cw+Cr/2))); Ig2=Ipr-Ir2-Iw2;
plot(t,Ig,'r',t,Ig1,'c-',t,Ig2,'b--','linewidth',2); xlabel('t/us');ylabel('A'); axis([0,90,-8,30]);
legend('高压点火线有电阻火花塞无电阻','高压点火线无电阻火花塞有电阻','高压点火线有电阻火花塞有电阻'); end
31