陕西科技大学毕业设计（论文）说明书 第32页
0.0384 ans = 0.9993
效果图如图3－4所示：
得到
Kv?f?v??-1.4090?10-10v6?3.3099?10?8v5?3.0013?10?6v4?1.3206?10v3?0.0030?v2?0.0384?v?0.9993?4
对于精度为7级的齿轮而言，经表面硬化处理由《机械设计》表10－3知齿间载荷分布系数取为KH??KF??1.2
由《机械设计》表10－4查得7级精度，小齿轮相对支承非对称布置时，
*KH??1.12?0.18?1?0.6?d2??d2?0.23?10?3b?1.408?0.23?10?3b，齿高h?2ham?2m得
b/h?b故由《机械设计》图10－13得齿向载荷分布系数 2mKF??1.12?0.18?1?0.6?d2??d2?0.23?10?3b?1.408?0.23?10?3b
Ft?
2T d
待入各参数值，得到以下四个不等式：
2T??xxn?1.5?Kv?53??K??x1??zhou1??x12?1?x5x3?60000?g16?X??474.5??540
x1?x5x3?x122T??xxn?1.5?Kv?53??K??x1??zhou2??x12?1?x5x3x12?60000?g17?X??474.5??522.5x1?x5x3x12?x12n???xx?64x??i?2T12??K??x2??zhou2??1.5?Kv??1?x6x4?x12??60000?????g18?X??474.5??540ix2?x6x4?x12n???xx?64x??i?2T12??K??x2??zhou3??1.5?Kv??1?i?x12??60000?xx64??x12??g19?X??474.5??522.5iix2?x6x4?x12x12
（5） 齿轮的弯曲应力应不大于其许用值
由《机械设计》公式10－5知，齿根危险截面的弯曲强度条件式为
?F?2KTY1FaYSa???F?
?dm3z12齿形系数YFa及应力校正系数YSa可由《机械设计》表10－5知与齿轮齿数有关可以通过Matlab曲线拟合的方法进行公式化；
YSa的拟合源程序：
z=[17,18,19,20,21,22,23,24,25,26,27,28,29,30,35,40,45,50,60,70,80,90,100,150,200];
ysa=[1.52,1.53,1.54,1.55,1.56,1.57,1.575,1.58,1.59,1.595,1.60,1.61,1.62,1.625,1.65,1.67,1.68,1.70,1.73,1.75,1.77,1.78,1.79,1.83,1.865,];
y2=polyfit(z,ysa,5); y2 x=17:1:200;
subplot(2,2,2);plot(z,ysa);
subplot(2,2,4);plot(z,ysa,'ro',x,y2(1).*x.^5+y2(2).*x.^4+y2(3).*x.^3+y2(4).*x.^2+y2(5).*x+y2(6)); x=[17,18,19,20,21,22,23,24,25,26,27,28,29,30,35,40,45,50,60,70,80,90,100,150,200] y2(1).*x.^5+y2(2).*x.^4+y2(3).*x.^3+y2(4).*x.^2+y2(5).*x+y2(6) y2(1)
陕西科技大学毕业设计（论文）说明书 第34页
y2(2) y2(3) y2(4) y2(5) y2(6)
得到系数为： ans = 2.4308e-011 ans = -1.3553e-008 ans = 2.8989e-006 ans = -3.0694e-004 ans = 0.0178 ans = 1.2953 得到拟合公式为：
YSa?2.4308?10?11x5?1.3553?10?8x4?2.8989?10?6x3?3.0694?10x?0.0178x?1.2953?42
YFa拟合的源程序：
data=[17,2.97; 18,2.91; 19,2.85; 20,2.80; 21,2.76; 22,2.72; 23,2.69; 24,2.65; 25,2.62; 26,2.60; 27,2.57; 28,2.55; 29,2.53; 30,2.52; 35,2.45; 40,2.40; 45,2.35;
50,2.32; 60,2.28; 70,2.24; 80,2.22; 90,2.20; 100,2.18; 150,2.14; 200,2.12 ];
init_lambda=[0,0];
lambda=fminsearch('fun_e3',init_lambda,[],data); x=data(:,1); y=data(:,2);
A=[exp(lambda(1)*x) exp(lambda(2)*x)]; a=A\\y;
estimated_y=a(1)*exp(lambda(1)*x)+a(2)*exp(lambda(2)*x) subplot(2,2,1);plot(x,y);
subplot(2,2,3);plot(x,y,'ro',x,estimated_y,'b-') lambda(1) lambda(2) a(1) a(2)
得到系数为： ans = -4.8639e-004 ans = -0.0794 ans = 2.3135 ans = 2.5255 得到拟合公式为：
YFa?2.3135?e
?4.8639?10?4x?2.5255?e?0.0794x
陕西科技大学毕业设计（论文）说明书 第36页
效果图如图3－5所示：
待入各参数值可得到以下四个不等式：
??xxn?3?Kv?53??K??x1?Tzhou1?YFa?x3??YSa?x3??60000?g20?X???303.57x1x52x3n??xxx5312?x123?Kv??60000??g21?X??n??xx?64x123?Kv??60000??g22?X??n??xx64?x123?Kv??60000??g23?X??????K??x1??Tzhou2YFa?x3x12??YSa?x3x12?????238.862x1x5x3x12????K??x2??Tzhou2YFa?x4??YSa?x4?????303.57x2x62x4???i??i???K??x2??Tzhou3YFa?x4??YSa?x4???x12??x12????238.86
?i?x2x62?x4??x12?（6） 轴的弯扭强度校核计算
由《机械设计》公式15－50知，轴的弯扭合成强度条件为
?M???T??ca????4???W2W????22M2???T?W2????1?