matlab的编辑器函数,matlab 在一个编辑器求多个函数的积分，刚学还不太会

最新推荐文章于 2023-11-22 21:46:06 发布

宝之大者

最新推荐文章于 2023-11-22 21:46:06 发布

阅读量1.5k

点赞数

文章标签： matlab的编辑器函数

本帖最后由 jelly_mu 于 2018-4-4 21:59 编辑

想求齿轮啮合刚度的曲线

%这是我编辑的脚本文件这是程序代码

z1 = 25; %齿数

z2 = 35;

m = 4; %模数

a = deg2rad(20); %齿轮压力角

c = 0.25; %顶隙系数

ha = 1; %齿顶高系数

x = 0; %变位系数

E = 2*10^9; %弹性模量

v = 0.3; %泊松比

L = 40; %齿宽

rb1 = m*z1*cos(a);

rb2 = m*z2*cos(a); %基圆半径

ra1 = m*(z1+2*ha);

ra2 = m*(z2+2*ha); %齿顶圆半经

inva = tan(a)-a; %渐开线函数

S = (pi/2+2*x*tan(a))*m; %分度圆齿厚

Sb1 = cos(a)*(S+m*z1*inva); %基圆齿厚

Sb2 = cos(a)*(S+m*z2*inva);

kh = pi*E*L/(4*(1-v^2)); %赫兹接触刚度

alphak2 = Sb1/(2*rb1);

rho1 = 34.2020;

alphak1 = rho1/rb1-alphak2;

alphak4 = Sb2/(2*rb2);

rho2 = 47.8828;

alphak3 = rho2/rb2-alphak4;

syms alpha1 beta1 gamma1 alpha2 beta2 gamma2;

f1 = 3*(1+cos(alphak1)*((alphak2-alpha1)*sin(alpha1)-cos(alpha1)))^2*(alphak2-alpha1)*cos(alpha1) ...

/(2*E*L*(sin(alpha1)+(alphak2-alpha1)*cos(alpha1)))^3;

f2 = 1.2*(1+v)*(alphak2-beta1)*cos(beta1)*(cos(alphak1))^2*alphak1 ...

/(E*L*(sin(beta1)+(alphak2-beta1)*cos(beta1)));

f3 = (alphak2-gamma1)*cos(gamma1)*(sin(alphak1))^2/(2*E*L*(sin(gamma1)+(alphak2-gamma1)*cos(gamma1)));

f4 = 3*(1+cos(alphak3)*((alphak4-alpha2)*sin(alpha2)-cos(alpha2)))^2*(alphak4-alpha2)*cos(alpha2) ...

/(2*E*L*(sin(alpha2)+(alphak4-alpha2)*cos(alpha2)))^3;

f5 = 1.2*(1+v)*(alphak4-beta2)*cos(beta2)*(cos(alphak3))^2*alphak3 ...

/(E*L*(sin(beta2)+(alphak4-beta2)*cos(beta2)));

f6 = (alphak4-gamma2)*cos(gamma2)*(sin(alphak3))^2/(2*E*L*(sin(gamma2)+(alphak4-gamma2)*cos(gamma2)));

inv1 = int(f1,alpha1,-alphak1,alphak2);

inv2 = int(f2,beta1,-alphak1,alphak2);

inv3 = int(f3,gamma1,-alphak1,alphak2);

inv4 = int(f1,alpha2,-alphak3,alphak4);

inv5 = int(f2,beta2,-alphak3,alphak4);

inv6 = int(f3,gamma2,-alphak3,alphak4);

kb1 = 1/inv1; %弯曲刚度

ks1 = 1/inv2; %剪切刚度

kc1 = 1/inv3; %轴向压缩刚度

kb2 = 1/inv4; %弯曲刚度

ks2 = 1/inv5; %剪切刚度

kc2 = 1/inv6; %轴向压缩刚度

K = 1/(1/kh+1/kb1+1/kb2+1/ks1+1/ks2+1/kc1+1/kc2)

plot(Kt)

输出结果怎么是这，不是我想要的啊，图形也不知道怎么画了。

>> Kt

K =

1/(int(-(740145109493709560989073156557*cos(beta1)*((39*beta1)/25 - 27306594293067297/450359962737049600))/(2535301200456458802993406410752*(80000000000*sin(beta1) - 80000000000*cos(beta1)*(beta1 - 700169084437623/18014398509481984))), beta1, -2928266506797653/9007199254740992, 700169084437623/18014398509481984) + int(-(7351342885398277*cos(gamma1)*(gamma1 - 700169084437623/18014398509481984))/(72057594037927936*(160000000000*sin(gamma1) - 160000000000*cos(gamma1)*(gamma1 - 700169084437623/18014398509481984))), gamma1, -2928266506797653/9007199254740992, 700169084437623/18014398509481984) + int(-(3*cos(alpha1)*(alpha1 - 700169084437623/18014398509481984)*((4267691476218981*cos(alpha1))/4503599627370496 + (4267691476218981*sin(alpha1)*(alpha1 - 700169084437623/18014398509481984))/4503599627370496 - 1)^2)/(160000000000*sin(alpha1) - 160000000000*cos(alpha1)*(alpha1 - 700169084437623/18014398509481984))^3, alpha1, -2928266506797653/9007199254740992, 700169084437623/18014398509481984) - (9705821984532438570026014034160952038482313691*cos(beta1)*((39*beta1)/25 - 27306594293067297/450359962737049600))/(91343852333181432387730302044767688728495783936*(80000000000*sin(beta1) - 80000000000*cos(beta1)*(beta1 - 700169084437623/18014398509481984))) - (96401130640100694239431278718051*cos(gamma1)*(gamma1 - 700169084437623/18014398509481984))/(2596148429267413814265248164610048*(160000000000*sin(gamma1) - 160000000000*cos(gamma1)*(gamma1 - 700169084437623/18014398509481984))) - (39340212588197589*cos(alpha1)*(alpha1 - 700169084437623/18014398509481984)*((4267691476218981*cos(alpha1))/4503599627370496 + (4267691476218981*sin(alpha1)*(alpha1 - 700169084437623/18014398509481984))/4503599627370496 - 1)^2)/(36028797018963968*(160000000000*sin(alpha1) - 160000000000*cos(alpha1)*(alpha1 - 700169084437623/18014398509481984))^3) + 2241217244448703/154742504910672534362390528)

错误使用 plot

数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。

出错 Kt (line 62)

plot(K)

2018-4-4 09:55 上传