计算机辅助公差设计.docVIP

计算机辅助公差设计 已知：h=300.01mm,d=300.02mm ，=900.05 试编写matlab程序，并求： F的尺寸及其公差 h,d, 的标准差对的敏感度和贡献率，画出直方图。 解： syms f h d th td tj tf dh dd dj vf gh gd gj theta=sym(theta) f=d/(2*sin(theta/2))+d/2-h %所求尺寸F dh=diff(f,h) %F对H的敏感度 dd=diff(f,d) %F对D的敏感度 dj=diff(f,theta) %F对theta的敏感度 vf=sqrt((dh*th/3)^2+(dd*td/3)^2+(dj*tj*pi/540)^2) tf=3*vf %F的偏差 gh=(dh*th/3)^2/((tf/3)^2)*100 %H对F的贡献度 gd=(dd*td/3)^2/((tf/3)^2)*100 %D对F的贡献度 gj=(dj*tj*pi/540)^2/((tf/3)^2)*100 %theta对F的贡献度 %代入数据进行计算 h=30 d=30 theta=pi/2 th=0.01 td=0.02 tj=0.000278*pi f=eval(f) dh=eval(dh) dd=eval(dd) dj=eval(dj) vf=eval(vf) tf=eval(tf) gh=eval(gh) gd=eval(gd) gj=eval(gj) MATLAB计算结果如下： theta = theta f = d/2 - h + d/(2*sin(theta/2)) dh = -1 dd = 1/(2*sin(theta/2)) + 1/2 dj = -(d*cos(theta/2))/(4*sin(theta/2)^2) vf = (th^2/9 + (td^2*(1/(2*sin(theta/2)) + 1/2)^2)/9 + (pi^2*d^2*tj^2*cos(theta/2)^2)/(4665600*sin(theta/2)^4))^(1/2) tf = 3*(th^2/9 + (td^2*(1/(2*sin(theta/2)) + 1/2)^2)/9 + (pi^2*d^2*tj^2*cos(theta/2)^2)/(4665600*sin(theta/2)^4))^(1/2) gh = (100*th^2)/(9*(th^2/9 + (td^2*(1/(2*sin(theta/2)) + 1/2)^2)/9 + (pi^2*d^2*tj^2*cos(theta/2)^2)/(4665600*sin(theta/2)^4))) gd = (100*td^2*(1/(2*sin(theta/2)) + 1/2)^2)/(9*(th^2/9 + (td^2*(1/(2*sin(theta/2)) + 1/2)^2)/9 + (pi^2*d^2*tj^2*cos(theta/2)^2)/(4665600*sin(theta/2)^4))) gj = (pi^2*d^2*tj^2*cos(theta/2)^2)/(46656*sin(theta/2)^4*(th^2/9 + (td^2*(1/(2*sin(theta/2)) + 1/2)^2)/9 + (pi^2*d^2*tj^2*cos(theta/2)^2)/(4665600*sin(theta/2)^4))) h = 30 d = 30 theta = 1.5708 th = 0.0100 td = 0.0200 tj = 8.7336e-004 f = 6.2132 dh = -1 dd = 1.2071 dj = -10.6066 vf = 0.0087 tf = 0.0261 gh = 14.6441 gd = 85.3521 gj = 0.0038 由计算可得如下结果： F=6.2132±0.0261 F对H的敏感度 dh=-1 F对D的敏感度 dd=1.2071 F对theta的敏感度 dj=-10.6066 F的偏差 tf=0.0261 H对F的贡献度